1/*! 2 * \file castemincompressibleplanestressgornetdesmoratbehaviourbiaxialtesting.mtest 3 * \brief 4 * \author Helfer Thomas 5 * \brief 09 avril 2013 6 */ 7 8@AccelerationAlgorithm 'AlternateSecant'; 9@ModellingHypothesis 'PlaneStress'; 10@Behaviour<umat> '@MFrontCastemBehavioursBuildPath@' 'umatincompressibleplanestressgornetdesmoratbehaviour'; 11 12@Real 'C_1' '0.183'; 13@Real 'C_2' '0.0034'; 14@Real 'E' '3*(2*(C_1+C_2))'; 15 16@MaterialProperty<constant> 'YoungModulus' '100*E'; 17@MaterialProperty<constant> 'PoissonRatio' 0.49999; 18@MaterialProperty<constant> 'H1' '0.142236'; 19@MaterialProperty<constant> 'H2' '1.5854659e-2'; 20@MaterialProperty<constant> 'H3' '3.4946541e-4'; 21 22@ExternalStateVariable 'Temperature' {0:293.15,3600.:800}; 23 24@ImposedDeformationGradient 'FXX' 1.; 25@ImposedDeformationGradient 'FYY' 1.; 26@ImposedDeformationGradient 'FXY' {0:0,1:3.}; 27@ImposedDeformationGradient 'FYX' 0; 28 29@Times {0.,1 in 10}; 30 31@Test<function> {'SXX':'2*H1*exp(H3*FXY**4)*FXY**2'} 1.e-12; 32@Test<function> {'SYY':'-6*H2/sqrt(FXY**2+3)*FXY**2'} 1.e-12; 33@Test<function> {'SXY':'2*sqrt(2)*(H1*exp(H3*FXY**4)+3*H2/sqrt(FXY**2+3))*FXY'} 1.e-12; 34 35