*HEADING
 Import with Hashin's composite damage model; from /Std to /Xpl (Standard analysis)
 uniaxial compression and tension (off-axis angle 45 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,45.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
*DAMAGE STABILIZATION
 0.001,0.001,0.001,0.001
*amplitude, name=linear_amp, time=total time
 0.0, 0.0, 0.05, 1.0, 0.15, -1.0, 0.2, 0.0,
 0.25, 1.0, 0.35, -1.0, 0.4, 0.0
**---------------------------------------------
*STEP, nlgeom=yes
*Static
 0.005, 0.2, ,0.005
*BOUNDARY, amplitude=linear_amp
 ASSEMBLY-1.PART-1-1.move, 1,1, -0.01
*Boundary
 ASSEMBLY-1.PART-1-1.fix1, 1,1
 ASSEMBLY-1.PART-1-1.fix2, 2,2
 ASSEMBLY-1.PART-1-1.fix3, 3,3
*OUTPUT, FIELD
*ELEMENT OUTPUT, elset=PART-1-1.EA
 S,E,
 DAMAGEFT,DAMAGEFC,
 DAMAGEMT,DAMAGEMC,
 DAMAGESHR,
*NODE OUTPUT
 U,
*OUTPUT,HISTORY, freq=1
*ELEMENT OUTPUT,ELSET=qa_test_ea
 ELDMD,
 HSNFTCRT,HSNMTCRT,HSNFCCRT,HSNMCCRT,
 DAMAGEFT,DAMAGEMT,DAMAGESHR,DAMAGEFC,DAMAGEMC,
 STATUS,
*node output, nset=qa_test_fix1
 RF1,
*node output, nset=app_bc
 U1,
*Output, history
*Energy Output, variable=preselect
*END STEP