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@Behaviour<abaqus> '@MFrontAbaqusBehavioursBuildPath@' 'ABAQUSORTHOTROPICELASTIC_3D'; 9 10@Real 'E1' '7.8e+10'; 11@Real 'E2' '2.64233e+11'; 12@Real 'E3' '2.64233e+11'; 13@Real 'nu12' '0.0039991'; 14@Real 'nu13' '0.0039991'; 15@Real 'nu23' '0.00736051'; 16@Real 'G12' '4.8e+10'; 17@Real 'G13' '4.8e+10'; 18@Real 'G23' '1.16418e+11'; 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> 'ShearModulus13' 'G13'; 31@MaterialProperty<constant> 'ShearModulus23' 'G23'; 32@MaterialProperty<constant> 'ThermalExpansion1' 'a1'; 33@MaterialProperty<constant> 'ThermalExpansion2' 'a2'; 34@MaterialProperty<constant> 'ThermalExpansion3' 'a3'; 35 36@ExternalStateVariable 'Temperature' {0:293.15,3600.:800}; 37@Real 'ThermalExpansionReferenceTemperature' 293.15; 38@Real 'Tref' 'ThermalExpansionReferenceTemperature'; 39 40@ImposedStrain 'EXX' 0.; 41@ImposedStrain 'EYY' 0.; 42@ImposedStrain 'EZZ' 0.; 43@ImposedStrain 'EXY' 0.; 44@ImposedStrain 'EXZ' 0.; 45@ImposedStrain 'EYZ' 0.; 46 47@Real 'pi' '4*atan(1)'; 48@Real 'phi' 'pi/7'; 49@Real 'c' 'cos(phi)'; 50@Real 's' 'sin(phi)'; 51@RotationMatrix {{'c' ,0,'s'}, 52 {0 ,1,0}, 53 {'-s',0,'c'}}; 54 55@Times {0.,3600 in 20}; 56 57@Real 'S11' '1/E1'; 58@Real 'S22' '1/E2'; 59@Real 'S33' '1/E3'; 60@Real 'S12' '-nu12/E1'; 61@Real 'S13' '-nu13/E1'; 62@Real 'S23' '-nu23/E2'; 63@Real 'dS' 'S11*S22*S33+2*S23*S13*S12-S11*S23*S23-S22*S13*S13-S33*S12*S12'; 64@Real 'C00' '(S22*S33-S23*S23)/dS'; 65@Real 'C01' '(S13*S23-S12*S33)/dS'; 66@Real 'C02' '(S12*S23-S13*S22)/dS'; 67@Real 'C10' '(S13*S23-S12*S33)/dS'; 68@Real 'C11' '(S11*S33-S13*S13)/dS'; 69@Real 'C12' '(S12*S13-S11*S23)/dS'; 70@Real 'C20' '(S12*S23-S13*S22)/dS'; 71@Real 'C21' '(S12*S13-S11*S23)/dS'; 72@Real 'C22' '(S11*S22-S12*S12)/dS'; 73 74@Real 'sig11' '-(C00*a1+C01*a2+C02*a3)'; 75@Real 'sig33' '-(C20*a1+C21*a2+C22*a3)'; 76 77@Test<function> 'SXX' '(c**2*sig11+s**2*sig33)*(Temperature-Tref)' 1.e-3; 78@Test<function> 'SYY' '-(C10*a1+C11*a2+C12*a3)*(Temperature-Tref)' 1.e-3; 79@Test<function> 'SZZ' '(s**2*sig11+c**2*sig33)*(Temperature-Tref)' 1.e-3; 80@Test<function> 'SXZ' 'sqrt(2)*s*c*(-sig11+sig33)*(Temperature-Tref)' 1.e-3; 81@Test<function> 'SXY' '0.' 1.e-3; 82@Test<function> 'SYZ' '0.' 1.e-3; 83