*heading
ply fabric damage with power law rate dep.
*preprint, model=yes
*parameter
 E1 = 28.45e3
 E2 = E1
 nu12 = 0.086
 G12 = 1.77e3
 E1c = 27e3
 E2c = E1c
 nu12c = nu12
 sigu1t = 541.5
 sigu2t = sigu1t
 sigu1c = 93
 sigu2c = sigu1c
 sigu12 = 42.191
 Gf1t = sigu1t*sigu1t/E1
 Gf1c = 1e9
 Gf2t = sigu2t*sigu2t/E2
 Gf2c = 1e9
 alpha12 = 0.4
#
# yield0 = 37.6184
 yield0 = 20.678
# hardC = 143.822
 hardC = 133.983
# hardP = 0.299564
 hardP = 0.197272
#
 lDelFlag = 1.0
 sdegMax = 0.9999
 eqpsMax = 0.5
*node, nset=nall
 1, 0.0, 0.0, 0.0
 2, 1.0, 0.0, 0.0
 3, 1.0, 1.0, 0.0
 4, 0.0, 1.0, 0.0
*nset, nset=bot
 1, 2
*nset, nset=top
 3, 4
*nset, nset=left
 1, 4
*nset, nset=right
 2, 3
*element, type=s4r, elset=eall
 1, 1, 2, 3, 4
*orientation, name=45Deg
 1.0, 1.0, 0.0, -1.0, 1.0, 0.0
 3, 0.0
*SECTION CONTROLS, name=sc, MAX DEGRADATION=<sdegMax>, ELEMENT DELETION=YES
*shell section, elset=eall, material=ply_fabric, poisson=material, CONTROLS=sc, orientation=45Deg
 1.0,
*material, name=ply_fabric, SRATE FACTOR=0.1
*density
 1.e-5
*ELASTIC, TYPE=BILAMINA
 <E1>, <E2>, <nu12>,<G12>,<E1c>,<E2c>,<nu12c>
*PLY FABRIC HARDENING, TYPE= JOHNSON COOK
 <yield0>, <hardC>, <hardP>
*rate dependent, type= POWER LAW
0.5,2.
*DAMAGE INITIATION, CRITERION=PLY FABRIC
 <sigu1t>, <sigu1c>, <sigu2t>, <sigu2c>, <sigu12>
*DAMAGE EVOLUTION, TYPE=ENERGY, SOFTENING=EXPONENTIAL 
 <Gf1t>, <Gf1c>, <Gf2t>, <Gf2c>, <alpha12>, 0.955
*PLY FABRIC FAILURE, NUM FIBERS=1
<eqpsMax>, 5.e-2, -5.e-2
*boundary
 nall, 3, 6
 bot, 2, 2
 left, 1, 1
*amplitude, name=amp, definition=smooth step
 0.0, 0.0, 1.0, 0.05
*step
*dynamic, explicit
 ,1.0
*boundary, amplitude=amp
 top, 2, 2, -1.0
 right, 1, 1, 1.0
*output, field, variable=preselect
*element output, elset=eall
 s, ne, sdeg, sdv, DMICRT, status, peeq
 DAMAGEF1T,DAMAGEF1C,DAMAGEF2T,DAMAGEF2C,DAMAGESHR
***edmicrtmax
***dmicrtmax
*output, history, variable=preselect
*element output, elset=eall
 s, ne, le, status, DMICRT, peeq
 DAMAGEF1T,DAMAGEF1C,DAMAGEF2T,DAMAGEF2C,DAMAGESHR
*node output, nset=nall
 rf1, rf2,u
*end step