*HEADING
 input deck to test Import with Hashin's composite damage model
 from /Std to /Xpl (Standard analysis)
 uniaxial tension and compression (off-axis angle 45 degree)
 element tested - m3d4r, m3d4, m3d3
*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=M3D4R, ELSET=EA1
   1, 1, 2, 3, 4
  *ELEMENT, TYPE=M3D4, ELSET=EA2
   10, 11, 12, 13, 14
  *ELEMENT, TYPE=M3D3, ELSET=EA3
   20, 21, 22, 23
   21, 21, 23, 24
  *ELSET, ELSET=EA
   EA1, EA2, EA3
  *MEMBRANE SECTION,ELSET=EA1,MATERIAL=GLASS_EPOXY,ORIENT=RECT
   1.,
  *MEMBRANE SECTION,ELSET=EA2,MATERIAL=GLASS_EPOXY,ORIENT=RECT
   1.,
  *MEMBRANE 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,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,
*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, 1.0, 0.35, -2.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.001
*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, freq=1
*ELEMENT OUTPUT,ELSET=qa_test_ea
 ELDMD,
 HSNFTCRT,HSNMTCRT,HSNFCCRT,HSNMCCRT,
 DAMAGEFT,DAMAGEMT,DAMAGESHR,DAMAGEFC,DAMAGEMC,
*node output, nset=qa_test_fix1
 RF1,
*node output, nset=app_bc
 U1,
*Output, history
*Energy Output, variable=preselect
*END STEP