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