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