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 'AxisymmetricalGeneralisedPlaneStrain'; 9@Behaviour<castem,LogarithmicStrain1D> '@MFrontCastemBehavioursBuildPath@' 'umatorthotropicelastic'; 10 11@Real 'E1' '7.8e+10'; 12@Real 'E2' '2.64233e+11'; 13@Real 'E3' '2.64233e+11'; 14@Real 'nu12' '0.0039991'; 15@Real 'nu13' '0.0039991'; 16@Real 'nu23' '0.00736051'; 17@Real 'G12' '4.8e+10'; 18@Real 'G13' '4.8e+10'; 19@Real 'G23' '1.16418e+11'; 20@Real 'a1' '3.41e-5'; 21@Real 'a2' '2.23e-5'; 22@Real 'a3' '4.63e-5'; 23 24@MaterialProperty<constant> 'YoungModulus1' 'E1'; 25@MaterialProperty<constant> 'YoungModulus2' 'E2'; 26@MaterialProperty<constant> 'YoungModulus3' 'E3'; 27@MaterialProperty<constant> 'PoissonRatio12' 'nu12'; 28@MaterialProperty<constant> 'PoissonRatio13' 'nu13'; 29@MaterialProperty<constant> 'PoissonRatio23' 'nu23'; 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 'ETT' 0.; 40@ImposedStrain 'EZZ' 0.; 41 42@Times {0.,3600 in 20}; 43 44@Real 'S11' '1/E1'; 45@Real 'S22' '1/E2'; 46@Real 'S33' '1/E3'; 47@Real 'S12' '-nu12/E1'; 48@Real 'S13' '-nu13/E1'; 49@Real 'S23' '-nu23/E2'; 50@Real 'dS' 'S11*S22*S33+2*S23*S13*S12-S11*S23*S23-S22*S13*S13-S33*S12*S12'; 51@Real 'C00' '(S22*S33-S23*S23)/dS'; 52@Real 'C01' '(S13*S23-S12*S33)/dS'; 53@Real 'C02' '(S12*S23-S13*S22)/dS'; 54@Real 'C10' '(S13*S23-S12*S33)/dS'; 55@Real 'C11' '(S11*S33-S13*S13)/dS'; 56@Real 'C12' '(S12*S13-S11*S23)/dS'; 57@Real 'C20' '(S12*S23-S13*S22)/dS'; 58@Real 'C21' '(S12*S13-S11*S23)/dS'; 59@Real 'C22' '(S11*S22-S12*S12)/dS'; 60@Test<function> 'SRR' '-(C00*ln(1+a1*(Temperature-Tref))+C01*ln(1+a2*(Temperature-Tref))+C02*ln(1+a3*(Temperature-Tref)))' 1.e-3; 61@Test<function> 'SZZ' '-(C10*ln(1+a1*(Temperature-Tref))+C11*ln(1+a2*(Temperature-Tref))+C12*ln(1+a3*(Temperature-Tref)))' 1.e-3; 62@Test<function> 'STT' '-(C20*ln(1+a1*(Temperature-Tref))+C21*ln(1+a2*(Temperature-Tref))+C22*ln(1+a3*(Temperature-Tref)))' 1.e-3; 63