*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
1,0.,0.
2,1.0,0.0
3,1.0,1.0
4,0,1.0
*ELEMENT,TYPE=S4R,ELSET=EALL
1,1,2,3,4
*ORIENTATION,NAME=ORI
1.0,0.0,0.0,0.0,1.0,0.0
3,0.0
*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
 216453882.2	,0.08
 223689315.1	,0.09
 230397958.6	,0.1
 236663789.5	,0.11
 242551356.8	,0.12
 248111454.9	,0.13
 253384846.5	,0.14
 258404792.4	,0.15
 263198822.1	,0.16
 267790005.9	,0.17
 272197888.4	,0.18
 276439187	,0.19
 280528323.8	,0.2
 284477836.2	,0.21
 288298698.7	,0.22
 292000578.3	,0.23
 295592038.5	,0.24
 299080704.9	,0.25
 302473400.5	,0.26
 305776257.1	,0.27
 308994807.5	,0.28
 312134063.3	,0.29
 315198579.3	,0.3
 318192508.9	,0.31
 321119651.4	,0.32
 323983491.4	,0.33
 326787233.8	,0.34
 329533833.9	,0.35
 332226022.8	,0.36
 334866330.1	,0.37
 337457104.3	,0.38
 340000529.4	,0.39
 342498640.5	,0.4
 344953337.6	,0.41
 347366397.4	,0.42
 349739483.8	,0.43
 352074157.9	,0.44
 354371886.2	,0.45
 356634048.4	,0.46
 358861944.3	,0.47
 361056799.8	,0.48
 363219772.8	,0.49
 365351958.3	,0.5
 367454392.5	,0.51
 369528057.6	,0.52
 371573885.4	,0.53
 373592760.6	,0.54
 375585524.3	,0.55
 377552976.5	,0.56
 379495879.2	,0.57
 381414958.7	,0.58
 383310907.6	,0.59
 385184387.3	,0.6
 387036029.6	,0.61
 388866438.6	,0.62
 390676192.4	,0.63
 392465844.5	,0.64
 394235925.2	,0.65
 395986943	,0.66
 397719385.8	,0.67
 399433721.9	,0.68
 401130401.3	,0.69
 402809856.3	,0.7
 404472502.8	,0.71
 406118741	,0.72
 407748956.3	,0.73
 409363519.7	,0.74
 410962788.8	,0.75
 412547108.7	,0.76
 414116812	,0.77
 415672219.9	,0.78
 417213642.5	,0.79
 418741379.3	,0.8
 420255719.8	,0.81
 421756943.9	,0.82
 423245322.5	,0.83
 424721117.5	,0.84
 426184582.5	,0.85
 427635963.2	,0.86
 429075497.5	,0.87
 430503416.1	,0.88
 431919942.4	,0.89
 433325293.3	,0.9
 434719679.2	,0.91
 436103304.3	,0.92
 437476366.6	,0.93
 438839058.6	,0.94
 440191567.4	,0.95
 441534074.6	,0.96
 442866756.7	,0.97
 444189785.4	,0.98
 445503327.8	,0.99
 446807546.3	,1
*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
1,2
*NSET,NSET=LEFT
1,4
*NSET,NSET=TOP
3,4
*NSET,NSET=RIGHT
2,3
*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