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