1 /*!
2 * \file   norton.mtest
3 * \brief
4 * \author Helfer Thomas
5 * \brief  09 avril 2013
6 */
7
8@ModellingHypothesis 'Axisymmetrical';
9@AccelerationAlgorithm 'AlternateSecant';
10
11@Real 'Y'  150.e9;
12@Real 'nu' 0.3;
13@Behaviour<umat> '@MFrontCastemBehavioursBuildPath@' 'umatnorton';
14@MaterialProperty<constant> 'YoungModulus'      'Y';
15@MaterialProperty<constant> 'PoissonRatio'      'nu';
16@Real 'A' 8.e-67;
17@Real 'E' 8.2;
18@MaterialProperty<constant> 'NortonCoefficient' 'A';
19@MaterialProperty<constant> 'NortonExponent' 'E';
20
21@Real 'err' 1.5e-4;
22@Real 'srr' 'YoungModulus*err';
23@ImposedStrain 'ERR' 'err';
24// Initial value of the elastic strain
25@Real 'EELRR0' 'err';
26@Real 'EELZZ0' '-PoissonRatio*err';
27@InternalStateVariable 'ElasticStrain' {'EELRR0','EELZZ0','EELZZ0',0.};
28// Initial value of the total strain
29@Strain {'EELRR0','EELZZ0','EELZZ0',0.};
30// Initial value of the total stresses
31@Stress {'srr',0.,0.,0.};
32
33@ExternalStateVariable 'Temperature' 293.15;
34
35@Times {0.,3600 in 100};
36
37// // tests on internal state variables
38@Test<function> 'ElasticStrainRR' '(srr/Y)/((1+(E-1)*A*Y*t*srr**(E-1))**(1./(E-1)))'      1.5e-6;
39@Test<function> 'ElasticStrainTT' '-nu*(srr/Y)/((1+(E-1)*A*Y*t*srr**(E-1))**(1./(E-1)))'  1.5e-6;
40@Test<function> 'ElasticStrainZZ' '-nu*(srr/Y)/((1+(E-1)*A*Y*t*srr**(E-1))**(1./(E-1)))'  1.5e-6;
41@Test<function> 'ElasticStrainRZ' '0.'      1.e-12;
42@Test<function> 'EquivalentViscoplasticStrain' '(srr/Y)*(1-1/((1+(E-1)*A*Y*t*srr**(E-1))**(1./(E-1))))' 1.5e-6;
43// // this test is a bit paranoiac since SRR is imposed
44@Test<function> 'SRR' 'srr/((1+(E-1)*A*Y*t*srr**(E-1))**(1./(E-1)))' 1.e6;
45// check that the mechanical equilibrium is satisfied
46@Test<function> 'SZZ' '0.'      1.e-3;
47@Test<function> 'STT' '0.'      1.e-3;
48@Test<function> 'SRZ' '0.'      1.e-3;
49