*HEADING
  Paremetric study of the effect of prestrain on fld curve:biaxial strain
*PARAMETER
 major_initial=2.0
 minor_initial=2.0
 major=2.0
 beta=-1.0
 minor=beta*major
 a0=-1.0
 a_initial=1.0
 t0=0.0
 a_major=exp(0.05*major_initial)
 a_minor=exp(0.05*minor_initial)
 td_major_initial=-1.0/major_initial
 td_minor_initial=-1.0/minor_initial
 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_major>,<t0>,<td_major>
*AMPLITUDE,NAME=AMP_MINOR,DEFINITION=DECAY
<a0>,<a_minor>,<t0>,<td_minor>
*AMPLITUDE,NAME=INITIAL_MAJOR,DEFINITION=DECAY
<a0>,<a_initial>,<t0>,<td_major_initial>
*AMPLITUDE,NAME=INITIAL_MINOR,DEFINITION=DECAY
<a0>,<a_initial>,<t0>,<td_minor_initial>
*STEP
*DYNAMIC,EXPLICIT
,0.05
*BOUNDARY,AMPLITUDE=INITIAL_MAJOR
TOP,2,2,1.0
*BOUNDARY,AMPLITUDE=INITIAL_MINOR
RIGHT,1,1,1.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