1/*!
2 * \file   norton.mtest
3 * \brief
4 * \author Helfer Thomas
5 * \brief  09 avril 2013
6 */
7
8@ModellingHypothesis 'Axisymmetrical';
9
10@Behaviour<aster> '@MFrontAsterBehavioursBuildPath@' 'asternorton';
11@MaterialProperty<constant> 'YoungModulus'     150.e9;
12@MaterialProperty<constant> 'PoissonRatio'       0.3;
13@Real 'A' 8.e-67;
14@Real 'E' 8.2;
15@MaterialProperty<constant> 'NortonCoefficient' 'A';
16@MaterialProperty<constant> 'NortonExponent' 'E';
17
18@Real 'srr' 20e6;
19@ImposedStress 'SRR' 'srr';
20// Initial value of the elastic strain
21@Real 'EELRR0' 0.00013333333333333333;
22@Real 'EELZZ0' -0.00004;
23@InternalStateVariable 'ElasticStrain' {'EELRR0','EELZZ0','EELZZ0',0.};
24// Initial value of the total strain
25@Strain {'EELRR0','EELZZ0','EELZZ0',0.};
26// Initial value of the total stresses
27@Stress {'srr',0.,0.,0.};
28
29@ExternalStateVariable 'Temperature' 293.15;
30
31@Times {0.,3600 in 20};
32
33// tests on strains
34// note: EquivalentViscoplasticStrain is known at 1.e-12 (defaut value
35// for @StrainEpsilon), thus we may expect the strain to be known at
36// '3.6*1.e-9'. If pratice, things are a bit better but not much
37// better.
38@Test<function> 'ERR' 'EELRR0+A*SRR**E*t'     1.e-9;
39@Test<function> 'EZZ' 'EELZZ0-0.5*A*SRR**E*t' 1.e-10;
40@Test<function> 'ETT' 'EELZZ0-0.5*A*SRR**E*t' 1.e-10;
41@Test<function> 'ERZ' '0.'                    1.e-10;
42// tests on internal state variables
43@Test<function> 'ElasticStrainRR' 'EELRR0'  1.e-12;
44@Test<function> 'ElasticStrainZZ' 'EELZZ0'  1.e-12;
45@Test<function> 'ElasticStrainTT' 'EELZZ0'  1.e-12;
46@Test<function> 'ElasticStrainRZ' '0.'      1.e-12;
47@Test<function> 'EquivalentViscoplasticStrain' 'A*SRR**E*t' 1.e-12;
48// this test is a bit paranoiac since SRR is imposed
49@Test<function> 'SRR' 'SRR'     1.e-3;
50// check that the mechanical equilibrium is satisfied
51@Test<function> 'SZZ' '0.'      1.e-3;
52@Test<function> 'STT' '0.'      1.e-3;
53@Test<function> 'SRZ' '0.'      1.e-3;
54