*HEADING
 Import with Hashin's composite damage model; from /Xpl to /Std (Explicit analysis)
 uniaxial tension and compression (off-axis angle 0 degree)
 Shell elements - s4, s4r, s3r
*PREPRINT,ECHO=YES,MODEL=YES,HISTORY=YES
*restart, write
**--------------------------------------------- Part1
*PART,NAME=PART-1
  *NODE
   1,  0.0, 0.0
   2,  0.1, 0.0
   3,  0.1, 0.1
   4,  0.0, 0.1
   11, 0.0, 0.2
   12, 0.1, 0.2
   13, 0.1, 0.3
   14, 0.0, 0.3
   21, 0.0, 0.4
   22, 0.1, 0.4
   23, 0.1, 0.5
   24, 0.0, 0.5   
  *ELEMENT,TYPE=S4R, ELSET=EA1
   1, 1, 2, 3, 4
  *ELEMENT, TYPE=S4, ELSET=EA2
   10, 11, 12, 13, 14
  *ELEMENT, TYPE=S3R, ELSET=EA3
   20, 21, 22, 23
   21, 21, 23, 24
  *ELSET, ELSET=EA
   EA1, EA2, EA3
  *SHELL SECTION,ELSET=EA1,MATERIAL=GLASS_EPOXY,ORIENT=RECT
   1.,
  *SHELL SECTION,ELSET=EA2,MATERIAL=GLASS_EPOXY,ORIENT=RECT
   1.,
  *SHELL SECTION,ELSET=EA3,MATERIAL=GLASS_EPOXY,ORIENT=RECT
   1.,
  *NSET,NSET=fix1
   1, 4, 11, 14, 21, 24
  *NSET,NSET=fix2
   1, 2, 11, 12, 21, 22
  *NSET,NSET=fix3
   1, 2, 3, 4, 11, 12, 13, 14,
   21, 22, 23, 24
  *nset, nset=move
   2, 3, 12, 13, 22, 23
  *ORIENTATION,NAME=RECT
   1.0, 0.0, 0.0, 0.0, 1.0, 0.0
   3,0.0
*END PART
**--------------------------------------------- Assembly-1
*ASSEMBLY,NAME=ASSEMBLY-1
*INSTANCE,NAME=PART-1-1,PART=PART-1
*END INSTANCE
*END ASSEMBLY
**--------------------------------------------- End Assembly
*elset, elset=qa_test_ea
 PART-1-1.EA,
*nset, nset=qa_test_fix1
 ASSEMBLY-1.PART-1-1.fix1,
*nset, nset=app_bc
 ASSEMBLY-1.PART-1-1.move,
*MATERIAL,NAME=GLASS_EPOXY
*Density
 1900.0,
*ELASTIC,TYPE=LAMINA
 53.8e9,17.9e9,0.25,8.96e9,8.96e9,6.88e9
**<E11>,<E22>,<nu12>,<G12>,<G13>,<G23>
*DAMAGE INITIATION,CRITERION=HASHIN,ALPHA=1.0
 1034e6,1034e6,27.6e6,138e6,41.4e6,69e6
** <X_T>,<X_C>,<Y_T>,<Y_C>,<S_L>,<S_T>
*DAMAGE EVOLUTION, TYPE=ENERGY, SOFTENING=LINEAR
 20e6,20e6,4.0e4,1.0e6
**DAMAGE STABILIZATION
** 0.0001,0.0001,0.0001,0.0001
*amplitude, name=linear_amp, time=total time
 0.0, 0.0, 0.05, 1.0, 0.15, -2.0, 0.2, 0.0,
 0.25, 2.0, 0.35, -3.0, 0.4, 0.0
**--------------------------------------------
*STEP, nlgeom=yes
*Dynamic, Explicit, scale factor=0.5
 ,0.2
*BOUNDARY, amplitude=linear_amp
 ASSEMBLY-1.PART-1-1.move, 1,1, 0.005
*Boundary
 ASSEMBLY-1.PART-1-1.fix1, 1,1
 ASSEMBLY-1.PART-1-1.fix2, 2,2
 ASSEMBLY-1.PART-1-1.fix3, 3,6
*OUTPUT, FIELD
*ELEMENT OUTPUT, elset=PART-1-1.EA
 S,E,
 DAMAGEFT,DAMAGEFC,
 DAMAGEMT,DAMAGEMC,
 DAMAGESHR,
*NODE OUTPUT
 U,
*OUTPUT,HISTORY, time interval=1.0e-3
*ELEMENT OUTPUT,ELSET=PART-1-1.EA
 ELDMD,
 HSNFTCRT,HSNMTCRT,HSNFCCRT,HSNMCCRT,
 DAMAGEFT,DAMAGEMT,DAMAGESHR,DAMAGEFC,DAMAGEMC,
*ELEMENT OUTPUT,ELSET=qa_test_ea
 ELDMD,
 DAMAGESHR,
*node output, nset=qa_test_fix1
 RF1,
*node output, nset=app_bc
 U1,
*Output, history
*Energy Output, variable=preselect
*END STEP