*heading
ply fabric damage with JC 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 = 20.678
 hardC = 133.983
 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.
*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= TABULAR
  2.0678E+01,  0.0000E+00
  6.0155E+01,  2.0408E-03
  6.5940E+01,  4.0816E-03
  6.9709E+01,  6.1224E-03
  7.2572E+01,  8.1633E-03
  7.4908E+01,  1.0204E-02
  7.6894E+01,  1.2245E-02
  7.8629E+01,  1.4286E-02
  8.0176E+01,  1.6327E-02
  8.1575E+01,  1.8367E-02
  8.2854E+01,  2.0408E-02
  8.4034E+01,  2.2449E-02
  8.5131E+01,  2.4490E-02
  8.6157E+01,  2.6531E-02
  8.7121E+01,  2.8571E-02
  8.8031E+01,  3.0612E-02
  8.8894E+01,  3.2653E-02
  8.9715E+01,  3.4694E-02
  9.0498E+01,  3.6735E-02
  9.1247E+01,  3.8776E-02
  9.1964E+01,  4.0816E-02
  9.2654E+01,  4.2857E-02
  9.3317E+01,  4.4898E-02
  9.3957E+01,  4.6939E-02
  9.4575E+01,  4.8980E-02
  9.5172E+01,  5.1020E-02
  9.5751E+01,  5.3061E-02
  9.6312E+01,  5.5102E-02
  9.6857E+01,  5.7143E-02
  9.7386E+01,  5.9184E-02
  9.7901E+01,  6.1224E-02
  9.8402E+01,  6.3265E-02
  9.8890E+01,  6.5306E-02
  9.9366E+01,  6.7347E-02
  9.9831E+01,  6.9388E-02
  1.0028E+02,  7.1429E-02
  1.0073E+02,  7.3469E-02
  1.0116E+02,  7.5510E-02
  1.0159E+02,  7.7551E-02
  1.0200E+02,  7.9592E-02
  1.0241E+02,  8.1633E-02
  1.0281E+02,  8.3673E-02
  1.0320E+02,  8.5714E-02
  1.0358E+02,  8.7755E-02
  1.0396E+02,  8.9796E-02
  1.0433E+02,  9.1837E-02
  1.0469E+02,  9.3878E-02
  1.0505E+02,  9.5918E-02
  1.0540E+02,  9.7959E-02
  1.0575E+02,  1.0000E-01
  1.0736E+02,  1.1000E-01
  1.1057E+02,  1.3224E-01
  1.1337E+02,  1.5449E-01
  1.1586E+02,  1.7673E-01
  1.1812E+02,  1.9898E-01
  1.2017E+02,  2.2122E-01
  1.2207E+02,  2.4347E-01
  1.2384E+02,  2.6571E-01
  1.2549E+02,  2.8796E-01
  1.2703E+02,  3.1020E-01
  1.2850E+02,  3.3245E-01
  1.2988E+02,  3.5469E-01
  1.3120E+02,  3.7694E-01
  1.3246E+02,  3.9918E-01
  1.3366E+02,  4.2143E-01
  1.3481E+02,  4.4367E-01
  1.3592E+02,  4.6592E-01
  1.3699E+02,  4.8816E-01
  1.3801E+02,  5.1041E-01
  1.3901E+02,  5.3265E-01
  1.3996E+02,  5.5490E-01
  1.4089E+02,  5.7714E-01
  1.4179E+02,  5.9939E-01
  1.4267E+02,  6.2163E-01
  1.4352E+02,  6.4388E-01
  1.4434E+02,  6.6612E-01
  1.4515E+02,  6.8837E-01
  1.4593E+02,  7.1061E-01
  1.4669E+02,  7.3286E-01
  1.4744E+02,  7.5510E-01
  1.4817E+02,  7.7735E-01
  1.4888E+02,  7.9959E-01
  1.4957E+02,  8.2184E-01
  1.5025E+02,  8.4408E-01
  1.5092E+02,  8.6633E-01
  1.5157E+02,  8.8857E-01
  1.5221E+02,  9.1082E-01
  1.5284E+02,  9.3306E-01
  1.5346E+02,  9.5531E-01
  1.5406E+02,  9.7755E-01
  1.5466E+02,  9.9980E-01
  1.5524E+02,  1.0220E+00
  1.5581E+02,  1.0443E+00
  1.5637E+02,  1.0665E+00
  1.5693E+02,  1.0888E+00
  1.5747E+02,  1.1110E+00
  1.5801E+02,  1.1333E+00
  1.5854E+02,  1.1555E+00
  1.5906E+02,  1.1778E+00
  1.5957E+02,  1.2000E+00
*rate dependent, type= JOHNSON COOK
0.5, 0.05
*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
*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