*HEADING
  Paremetric study of prestrain: uniaxial strain.
*PARAMETER
 major=2.0
 beta=-1.0
 minor=beta*major
 a0=-1.0
 a=1.0
 t0=0.0
 t0_major=-0.05
 td_major=-1.0/major
 td_minor=-1.0/minor
*PREPRINT,MODEL=YES
*NODE,NSET=NALL
*ELEMENT,TYPE=S4R,ELSET=EALL
*ORIENTATION,NAME=ORI
*SHELL SECTION,ELSET=EALL,MATERIAL=MAT,ORIENTATION=ORI
1.0
*MATERIAL,NAME=MAT
*ELASTIC
2.648E+10,0.167
*DENSITY
2400
*PLASTIC
 88352278.13	,0
 128052248.4	,0.01
 150009607.2	,0.02
 166128903.5	,0.03
 179154007.3	,0.04
 190217647.5	,0.05
 199909364	,0.06
 208579231.4	,0.07
*DAMAGE INITIATION, CRITERION=MK, FNN = 0.0, FNT= 0.0, FEQ= 0.0, NUMBER IMPERFECTION=100
 0.9999, 0.0
*DAMAGE EVOLUTION,TYPE=ENERGY, SOFTENING=EXPONENTIAL
 0.01
*NSET,NSET=BOTTOM
*NSET,NSET=LEFT
*NSET,NSET=TOP
*NSET,NSET=RIGHT
*BOUNDARY,TYPE=VELOCITY
BOTTOM,2,2
LEFT,1,1
NALL,4,6
*AMPLITUDE,NAME=AMP_MAJOR,DEFINITION=DECAY
<a0>,<a>,<t0_major>,<td_major>
*AMPLITUDE,NAME=AMP_MINOR,DEFINITION=DECAY
<a0>,<a>,<t0>,<td_minor>
*AMPLITUDE,NAME=AMP_INITIAL,DEFINITION=DECAY
<a0>,<a>,<t0>,<td_major>
*STEP
*DYNAMIC,EXPLICIT
,0.05
*BOUNDARY,AMPLITUDE=AMP_INITIAL
TOP,2,2,1.0
*BOUNDARY,TYPE=VELOCITY
RIGHT,1,1,0.0
*OUTPUT,FIELD,NUMBER INTERVAL=20, TIME MARKS=YES
*NODE OUTPUT,NSET=NALL
 U,
*OUTPUT,HISTORY,TIME INTERVAL=0.002
*ELEMENT OUTPUT,ELSET=EALL
SDEG,STATUS,MKCRT
LEP,PEEQ
*END STEP
*STEP
*DYNAMIC,EXPLICIT
,1.0
*BOUNDARY,AMPLITUDE=AMP_MAJOR
TOP,2,2,1.0
*BOUNDARY,AMPLITUDE=AMP_MINOR
RIGHT,1,1,1.0
*END STEP