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