*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 - sc6r, sc8r
*PREPRINT,ECHO=YES,MODEL=YES,HISTORY=YES
*restart, write
**--------------------------------------------- Part1
*PART,NAME=PART-1
  *NODE
   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
   5,  0.0, 0.0, 1.0
   6,  1.0, 0.0, 1.0
   7,  1.0, 1.0, 1.0
   8,  0.0, 1.0, 1.0
**
   11,  0.0, 0.0, 2.0
   12,  1.0, 0.0, 2.0
   13,  1.0, 1.0, 2.0
   14,  0.0, 1.0, 2.0
   15,  0.0, 0.0, 3.0
   16,  1.0, 0.0, 3.0
   17,  1.0, 1.0, 3.0
   18,  0.0, 1.0, 3.0  
  *ELEMENT,TYPE=SC8R
   1, 1, 2, 3, 4, 5, 6, 7, 8
  *ELEMENT, TYPE=SC6R
   11, 11, 12, 13, 15, 16, 17
   12, 11, 13, 14, 15, 17, 18
  *ELSET,ELSET=EA
   1, 11, 12
  *SHELL SECTION,ELSET=EA,MATERIAL=GLASS_EPOXY,ORIENT=RECT,controls=EC-1
   1.0,
  *NSET,NSET=fix1
   1, 4, 5, 8, 11, 14, 15, 18
  *NSET,NSET=fix2
   1, 2, 5, 6, 11, 12, 15, 16
  *NSET,NSET=fix3
   1, 2, 3, 4, 11, 12, 13, 14
  *nset, nset=move
   2, 3, 6, 7, 12, 13, 16, 17
  *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,
**-----------------------------------------------
*Section Controls, name=EC-1, hourglass=ENHANCED
*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=0.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
*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
**---------------------------------------------
*nset,nset=qa_test
qa_test_fix1, app_bc
*STEP, nlgeom=yes
*Dynamic, Explicit, scale factor=0.5
 ,0.2
*BOUNDARY, amplitude=linear_amp
 ASSEMBLY-1.PART-1-1.move, 1,1, 0.02
*Boundary
 ASSEMBLY-1.PART-1-1.fix1, 1,3
 ASSEMBLY-1.PART-1-1.fix2, 2,3
 ASSEMBLY-1.PART-1-1.fix3, 2,3
 ASSEMBLY-1.PART-1-1.move, 2,3
*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,
 STATUS,
*ELEMENT OUTPUT,ELSET=qa_test_ea
 ELDMD,
 DAMAGESHR,
*node output, nset=qa_test_fix1
 RF1,
*node output, nset=app_bc
 U1,
*OUTPUT,FIELD
*node output, nset=qa_test
 U,
*Output, history
*Energy Output, variable=preselect
*END STEP