/*! * \file castemnortonrk2-rk54.mtest * \brief * \author Helfer Thomas * \brief 09 avril 2013 */ @ModellingHypothesis 'Axisymmetrical'; @PredictionPolicy 'ElasticPrediction'; @AccelerationAlgorithm 'AlternateSecant'; @Behaviour '@MFrontCastemBehavioursBuildPath@' 'umatnortonrk4'; @MaterialProperty 'YoungModulus' 150.e9; @MaterialProperty 'PoissonRatio' 0.3; @Real 'E' 8.2; @MaterialProperty 'E[0]' 'E'; @MaterialProperty 'E[1]' 'E'; @MaterialProperty 'E[2]' 'E'; @MaterialProperty 'E[3]' 'E'; @MaterialProperty 'E[4]' 'E'; @MaterialProperty 'E[5]' 'E'; @MaterialProperty 'E[6]' 'E'; @MaterialProperty 'E[7]' 'E'; @MaterialProperty 'E[8]' 'E'; @MaterialProperty 'E[9]' 'E'; @MaterialProperty 'E[10]' 'E'; @MaterialProperty 'E[11]' 'E'; @MaterialProperty 'E[12]' 'E'; @MaterialProperty 'E[13]' 'E'; @MaterialProperty 'E[14]' 'E'; @MaterialProperty 'E[15]' 'E'; @MaterialProperty 'E[16]' 'E'; @MaterialProperty 'E[17]' 'E'; @MaterialProperty 'E[18]' 'E'; @MaterialProperty 'E[19]' 'E'; @MaterialProperty 'E[20]' 'E'; @MaterialProperty 'E[21]' 'E'; @MaterialProperty 'E[22]' 'E'; @MaterialProperty 'E[23]' 'E'; @MaterialProperty 'E[24]' 'E'; @MaterialProperty 'E[25]' 'E'; @MaterialProperty 'E[26]' 'E'; @MaterialProperty 'E[27]' 'E'; @MaterialProperty 'E[28]' 'E'; @MaterialProperty 'E[29]' 'E'; @MaterialProperty 'E[30]' 'E'; @MaterialProperty 'E[31]' 'E'; @MaterialProperty 'E[32]' 'E'; @MaterialProperty 'E[33]' 'E'; @MaterialProperty 'E[34]' 'E'; @MaterialProperty 'E[35]' 'E'; @MaterialProperty 'E[36]' 'E'; @MaterialProperty 'E[37]' 'E'; @MaterialProperty 'E[38]' 'E'; @MaterialProperty 'E[39]' 'E'; @MaterialProperty 'E[40]' 'E'; @MaterialProperty 'E[41]' 'E'; @MaterialProperty 'E[42]' 'E'; @MaterialProperty 'E[43]' 'E'; @MaterialProperty 'E[44]' 'E'; @MaterialProperty 'E[45]' 'E'; @MaterialProperty 'E[46]' 'E'; @MaterialProperty 'E[47]' 'E'; @MaterialProperty 'E[48]' 'E'; @MaterialProperty 'E[49]' 'E'; @MaterialProperty 'E[50]' 'E'; @MaterialProperty 'E[51]' 'E'; @MaterialProperty 'E[52]' 'E'; @MaterialProperty 'E[53]' 'E'; @MaterialProperty 'E[54]' 'E'; @MaterialProperty 'E[55]' 'E'; @MaterialProperty 'E[56]' 'E'; @MaterialProperty 'E[57]' 'E'; @MaterialProperty 'E[58]' 'E'; @MaterialProperty 'E[59]' 'E'; @MaterialProperty 'E[60]' 'E'; @MaterialProperty 'E[61]' 'E'; @MaterialProperty 'E[62]' 'E'; @MaterialProperty 'E[63]' 'E'; @MaterialProperty 'E[64]' 'E'; @MaterialProperty 'E[65]' 'E'; @MaterialProperty 'E[66]' 'E'; @MaterialProperty 'E[67]' 'E'; @MaterialProperty 'E[68]' 'E'; @MaterialProperty 'E[69]' 'E'; @MaterialProperty 'E[70]' 'E'; @MaterialProperty 'E[71]' 'E'; @MaterialProperty 'E[72]' 'E'; @MaterialProperty 'E[73]' 'E'; @MaterialProperty 'E[74]' 'E'; @MaterialProperty 'E[75]' 'E'; @MaterialProperty 'E[76]' 'E'; @MaterialProperty 'E[77]' 'E'; @MaterialProperty 'E[78]' 'E'; @MaterialProperty 'E[79]' 'E'; @MaterialProperty 'E[80]' 'E'; @MaterialProperty 'E[81]' 'E'; @MaterialProperty 'E[82]' 'E'; @MaterialProperty 'E[83]' 'E'; @MaterialProperty 'E[84]' 'E'; @MaterialProperty 'E[85]' 'E'; @MaterialProperty 'E[86]' 'E'; @MaterialProperty 'E[87]' 'E'; @MaterialProperty 'E[88]' 'E'; @MaterialProperty 'E[89]' 'E'; @MaterialProperty 'E[90]' 'E'; @MaterialProperty 'E[91]' 'E'; @MaterialProperty 'E[92]' 'E'; @MaterialProperty 'E[93]' 'E'; @MaterialProperty 'E[94]' 'E'; @MaterialProperty 'E[95]' 'E'; @MaterialProperty 'E[96]' 'E'; @MaterialProperty 'E[97]' 'E'; @MaterialProperty 'E[98]' 'E'; @MaterialProperty 'E[99]' 'E'; @Real 'srr' 20.e6; @ImposedStress 'SRR' 'srr'; // Initial value of the elastic strain @Real 'EELRR0' 0.00013333333333333333; @Real 'EELZZ0' -0.00004; @InternalStateVariable 'ElasticStrain' {'EELRR0','EELZZ0','EELZZ0',0.}; // Initial value of the total strain @Strain {'EELRR0','EELZZ0','EELZZ0',0.}; // Initial value of the stresses @Stress {'srr',0.,0.,0.}; @ExternalStateVariable 'Temperature' 293.15; @Real 'A' 8.e-67; @ExternalStateVariable 'A[0]' 'A/100'; @ExternalStateVariable 'A[1]' 'A/100'; @ExternalStateVariable 'A[2]' 'A/100'; @ExternalStateVariable 'A[3]' 'A/100'; @ExternalStateVariable 'A[4]' 'A/100'; @ExternalStateVariable 'A[5]' 'A/100'; @ExternalStateVariable 'A[6]' 'A/100'; @ExternalStateVariable 'A[7]' 'A/100'; @ExternalStateVariable 'A[8]' 'A/100'; @ExternalStateVariable 'A[9]' 'A/100'; @ExternalStateVariable 'A[10]' 'A/100'; @ExternalStateVariable 'A[11]' 'A/100'; @ExternalStateVariable 'A[12]' 'A/100'; @ExternalStateVariable 'A[13]' 'A/100'; @ExternalStateVariable 'A[14]' 'A/100'; @ExternalStateVariable 'A[15]' 'A/100'; @ExternalStateVariable 'A[16]' 'A/100'; @ExternalStateVariable 'A[17]' 'A/100'; @ExternalStateVariable 'A[18]' 'A/100'; @ExternalStateVariable 'A[19]' 'A/100'; @ExternalStateVariable 'A[20]' 'A/100'; @ExternalStateVariable 'A[21]' 'A/100'; @ExternalStateVariable 'A[22]' 'A/100'; @ExternalStateVariable 'A[23]' 'A/100'; @ExternalStateVariable 'A[24]' 'A/100'; @ExternalStateVariable 'A[25]' 'A/100'; @ExternalStateVariable 'A[26]' 'A/100'; @ExternalStateVariable 'A[27]' 'A/100'; @ExternalStateVariable 'A[28]' 'A/100'; @ExternalStateVariable 'A[29]' 'A/100'; @ExternalStateVariable 'A[30]' 'A/100'; @ExternalStateVariable 'A[31]' 'A/100'; @ExternalStateVariable 'A[32]' 'A/100'; @ExternalStateVariable 'A[33]' 'A/100'; @ExternalStateVariable 'A[34]' 'A/100'; @ExternalStateVariable 'A[35]' 'A/100'; @ExternalStateVariable 'A[36]' 'A/100'; @ExternalStateVariable 'A[37]' 'A/100'; @ExternalStateVariable 'A[38]' 'A/100'; @ExternalStateVariable 'A[39]' 'A/100'; @ExternalStateVariable 'A[40]' 'A/100'; @ExternalStateVariable 'A[41]' 'A/100'; @ExternalStateVariable 'A[42]' 'A/100'; @ExternalStateVariable 'A[43]' 'A/100'; @ExternalStateVariable 'A[44]' 'A/100'; @ExternalStateVariable 'A[45]' 'A/100'; @ExternalStateVariable 'A[46]' 'A/100'; @ExternalStateVariable 'A[47]' 'A/100'; @ExternalStateVariable 'A[48]' 'A/100'; @ExternalStateVariable 'A[49]' 'A/100'; @ExternalStateVariable 'A[50]' 'A/100'; @ExternalStateVariable 'A[51]' 'A/100'; @ExternalStateVariable 'A[52]' 'A/100'; @ExternalStateVariable 'A[53]' 'A/100'; @ExternalStateVariable 'A[54]' 'A/100'; @ExternalStateVariable 'A[55]' 'A/100'; @ExternalStateVariable 'A[56]' 'A/100'; @ExternalStateVariable 'A[57]' 'A/100'; @ExternalStateVariable 'A[58]' 'A/100'; @ExternalStateVariable 'A[59]' 'A/100'; @ExternalStateVariable 'A[60]' 'A/100'; @ExternalStateVariable 'A[61]' 'A/100'; @ExternalStateVariable 'A[62]' 'A/100'; @ExternalStateVariable 'A[63]' 'A/100'; @ExternalStateVariable 'A[64]' 'A/100'; @ExternalStateVariable 'A[65]' 'A/100'; @ExternalStateVariable 'A[66]' 'A/100'; @ExternalStateVariable 'A[67]' 'A/100'; @ExternalStateVariable 'A[68]' 'A/100'; @ExternalStateVariable 'A[69]' 'A/100'; @ExternalStateVariable 'A[70]' 'A/100'; @ExternalStateVariable 'A[71]' 'A/100'; @ExternalStateVariable 'A[72]' 'A/100'; @ExternalStateVariable 'A[73]' 'A/100'; @ExternalStateVariable 'A[74]' 'A/100'; @ExternalStateVariable 'A[75]' 'A/100'; @ExternalStateVariable 'A[76]' 'A/100'; @ExternalStateVariable 'A[77]' 'A/100'; @ExternalStateVariable 'A[78]' 'A/100'; @ExternalStateVariable 'A[79]' 'A/100'; @ExternalStateVariable 'A[80]' 'A/100'; @ExternalStateVariable 'A[81]' 'A/100'; @ExternalStateVariable 'A[82]' 'A/100'; @ExternalStateVariable 'A[83]' 'A/100'; @ExternalStateVariable 'A[84]' 'A/100'; @ExternalStateVariable 'A[85]' 'A/100'; @ExternalStateVariable 'A[86]' 'A/100'; @ExternalStateVariable 'A[87]' 'A/100'; @ExternalStateVariable 'A[88]' 'A/100'; @ExternalStateVariable 'A[89]' 'A/100'; @ExternalStateVariable 'A[90]' 'A/100'; @ExternalStateVariable 'A[91]' 'A/100'; @ExternalStateVariable 'A[92]' 'A/100'; @ExternalStateVariable 'A[93]' 'A/100'; @ExternalStateVariable 'A[94]' 'A/100'; @ExternalStateVariable 'A[95]' 'A/100'; @ExternalStateVariable 'A[96]' 'A/100'; @ExternalStateVariable 'A[97]' 'A/100'; @ExternalStateVariable 'A[98]' 'A/100'; @ExternalStateVariable 'A[99]' 'A/100'; @Real 'tmax' 2500; @Times {0.,100. in 100, 'tmax' in 1000}; // tests on strains // note: EquivalentViscoplasticStrain is known at 1.e-12 (defaut value // for @StrainEpsilon), thus we may expect the strain to be known at // '3.6*1.e-9'. If pratice, things are a bit better but not much // better. @Test 'ERR' 'EELRR0+A*SRR**E*t' 1.e-9; @Test 'EZZ' 'EELZZ0-0.5*A*SRR**E*t' 1.e-10; @Test 'ETT' 'EELZZ0-0.5*A*SRR**E*t' 1.e-10; @Test 'ERZ' '0.' 1.e-10; // tests on internal state variables @Test 'ElasticStrainRR' 'EELRR0' 1.e-12; @Test 'ElasticStrainZZ' 'EELZZ0' 1.e-12; @Test 'ElasticStrainTT' 'EELZZ0' 1.e-12; @Test 'ElasticStrainRZ' '0.' 1.e-12;