*heading
MISES PLASTICITY WITH COMBINED ISOTROPIC-KINEMATIC HARDENING
WITH MULTIPLE BACKSTRESSES,
UNIAXIAL TENSION (LOAD CONTROL)TEST; DIRECTION 1, C3D8
*part, name=part-1
*node,nset=alln
 1,0.,0.,0.
 2,1.,0.,0.
 3,1.,1.,0.
 4,0.,1.,0.
 5,0.,0.,1.
 6,1.,0.,1.
 7,1.,1.,1.
 8,0.,1.,1.
*element,type=c3d8r,elset=alle
 1,1,2,3,4,5,6,7,8
*solid section,elset=alle,material=alle, controls=ctrl
*end part
*assembly, name=assembly-1
*instance, name=part-1-1, part=part-1
*end instance
*end assembly
*section controls, name=ctrl, hourglass=enhanced
*material,name=alle
*elastic
 200.e3,.3
*density
 7.82000e-09,
*plastic, hardening=combined, data type=parameters, numberbackstress=3
 200., 2.222e4, 34.65
*cyclic hardening, parameters
 200., 2000., 0.25
*boundary
 assembly-1.part-1-1.1,1,3
 assembly-1.part-1-1.4,1,1
 assembly-1.part-1-1.4,3,3
 assembly-1.part-1-1.5,1,1
 assembly-1.part-1-1.5,2,2
 assembly-1.part-1-1.8,1,1
**step 1 tension***
*amplitude,name=smoothing_curve,definition=smooth step
 0., 0., 10e-3, 1.
*step,nlgeom=no
*dynamic,explicit
 ,10e-3
*dload,amplitude=smoothing_curve
 assembly-1.part-1-1.1, p4, -701. 
*output, field
*element output, elset=assembly-1.part-1-1.alle
 S, LE, PE, PEEQ, ALPHA
*node output, nset=assembly-1.part-1-1.alln
 U, RF
*output, history, variable=preselect
*element output, elset=assembly-1.part-1-1.alle
 S, PE, LE, PEEQ, ALPHA, ALPHAP
*file output, number interval = 10
*el file, elset=assembly-1.part-1-1.alle
 PEEQ, ALPHA, ALPHAP
*end step