1@Author Helfer Thomas;
2@Date  09 avril 2013;
3@Description{
4  "This test checks that thermal expansion are correctly taken into "
5  "account for an orthotropic material."
6};
7
8@ModellingHypothesis 'Axisymmetrical';
9
10@Behaviour<aster> '@MFrontAsterBehavioursBuildPath@' 'asterorthotropicelastic';
11
12@Real 'E1'   '7.8e+10';
13@Real 'E2'   '2.64233e+11';
14@Real 'E3'   '2.64233e+11';
15@Real 'nu12' '0.0039991';
16@Real 'nu13' '0.0039991';
17@Real 'nu23' '0.00736051';
18@Real 'G12'  '4.8e+10';
19@Real 'a1'   '3.41e-5';
20@Real 'a2'   '2.23e-5';
21@Real 'a3'   '4.63e-5';
22
23@MaterialProperty<constant> 'YoungModulus1'     'E1';
24@MaterialProperty<constant> 'YoungModulus2'     'E2';
25@MaterialProperty<constant> 'YoungModulus3'     'E3';
26@MaterialProperty<constant> 'PoissonRatio12'    'nu12';
27@MaterialProperty<constant> 'PoissonRatio13'    'nu13';
28@MaterialProperty<constant> 'PoissonRatio23'    'nu23';
29@MaterialProperty<constant> 'ShearModulus12'    'G12';
30@MaterialProperty<constant> 'ThermalExpansion1' 'a1';
31@MaterialProperty<constant> 'ThermalExpansion2' 'a2';
32@MaterialProperty<constant> 'ThermalExpansion3' 'a3';
33
34@ExternalStateVariable 'Temperature' {0:293.15,3600.:800};
35@Real 'ThermalExpansionReferenceTemperature' 293.15;
36@Real 'Tref' 'ThermalExpansionReferenceTemperature';
37
38@ImposedStrain 'ERR' 0.;
39@ImposedStrain 'EZZ' 0.;
40@ImposedStrain 'ETT' 0.;
41@ImposedStrain 'ERZ' 0.;
42
43@Real 'pi'   '4*atan(1)';
44@Real 'phi'  'pi/7';
45@Real 'c'    'cos(phi)';
46@Real 's'    'sin(phi)';
47@RotationMatrix {{'c' ,'s',0},
48                 {'-s','c',0},
49		 {0,0,1}};
50
51@Times {0.,3600 in 20};
52
53@Real 'S11' '1/E1';
54@Real 'S22' '1/E2';
55@Real 'S33' '1/E3';
56@Real 'S12' '-nu12/E1';
57@Real 'S13' '-nu13/E1';
58@Real 'S23' '-nu23/E2';
59@Real 'dS'  'S11*S22*S33+2*S23*S13*S12-S11*S23*S23-S22*S13*S13-S33*S12*S12';
60@Real 'C00' '(S22*S33-S23*S23)/dS';
61@Real 'C01' '(S13*S23-S12*S33)/dS';
62@Real 'C02' '(S12*S23-S13*S22)/dS';
63@Real 'C10' '(S13*S23-S12*S33)/dS';
64@Real 'C11' '(S11*S33-S13*S13)/dS';
65@Real 'C12' '(S12*S13-S11*S23)/dS';
66@Real 'C20' '(S12*S23-S13*S22)/dS';
67@Real 'C21' '(S12*S13-S11*S23)/dS';
68@Real 'C22' '(S11*S22-S12*S12)/dS';
69
70@Real 'sig11' '-(C00*a1+C01*a2+C02*a3)';
71@Real 'sig22' '-(C10*a1+C11*a2+C12*a3)';
72
73@Test<function> 'SRR' '(c**2*sig11+s**2*sig22)*(Temperature-Tref)'    1.e-3;
74@Test<function> 'SZZ' '(s**2*sig11+c**2*sig22)*(Temperature-Tref)'    1.e-3;
75@Test<function> 'SRZ' 'sqrt(2)*s*c*(-sig11+sig22)*(Temperature-Tref)' 1.e-3;
76@Test<function> 'STT' '-(C20*a1+C21*a2+C22*a3)*(Temperature-Tref)'    1.e-3;
77