*HEADING
*PARAMETER
***major can be interpreted as the strain rate
 major= 2.0
 beta= -0.6
 minor=beta*major
 t_step1 = 0.250724
 a0=-1.0
 t0=0.0
 a_step1=1.0
 td_major_step1=-1.0/major
 td_minor_step1=-1.0/minor
 a_major_step2= exp(t_step1*major)
 a_minor_step2= exp(t_step1*minor)
 major_step2 = 2.0
 beta_step2 = 0.3
 minor_step2 = beta_step2 * major_step2
 td_major_step2 = -1.0/major_step2
 td_minor_step2 = -1.0/minor_step2
*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_step1>,<t0>,<td_major_step1>
*AMPLITUDE,NAME=AMP_MINOR,DEFINITION=DECAY
<a0>,<a_step1>,<t0>,<td_minor_step1>
*AMPLITUDE,NAME=AMP_MAJOR2,DEFINITION=DECAY
<a0>,<a_major_step2>,<t0>,<td_major_step2>
*AMPLITUDE,NAME=AMP_MINOR2,DEFINITION=DECAY
<a0>,<a_minor_step2>,<t0>,<td_minor_step2>
*STEP
*DYNAMIC,EXPLICIT
,<t_step1>
*BOUNDARY,AMPLITUDE=AMP_MAJOR
TOP,2,2,1.0
*BOUNDARY,AMPLITUDE=AMP_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_MAJOR2
TOP,2,2,1.0
*BOUNDARY,AMPLITUDE=AMP_MINOR2
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