*HEADING
 Test element deletion (s4R elements); default setting
*PREPRINT,ECHO=YES,HISTORY=YES
**--------------------------------------------- Part1
*PART,NAME=PART-1
  *NODE,NSET=NALL
   1,  0.0
   2,  1.0
   3,  2.0
   4,  3.0
**
   5,  0.0, 1.0
   6,  1.0, 1.0
   7,  2.0, 1.0
   8,  3.0, 1.0
**
   9,  0.0, 2.0
  10,  1.0, 2.0
  11,  2.0, 2.0
  12,  3.0, 2.0
**
  13,  0.0, 3.0
  14,  1.0, 3.0
  15,  2.0, 3.0
  16,  3.0, 3.0
  *ELEMENT,TYPE=S4R, ELSET=EALL
   1,  1, 2, 6, 5
   2,  2, 3, 7, 6
   3,  3, 4, 8, 7
   4,  5, 6,10, 9
   5,  6, 7,11,10
   6,  7, 8,12,11
   7,  9,10,14,13
   8, 10,11,15,14
   9, 11,12,16,15
  *ELSET,ELSET=EELAST, GEN
   1,4
   6,9
  *ELSET,ELSET=EDAMAGE
   5,
  *ELSET,ELSET=EA
   1,
  *SHELL SECTION,ELSET=EELAST,MATERIAL=GLASS_EPOXY_DMG,ORIENT=RECT
   1.,
  *SHELL SECTION,ELSET=EDAMAGE,MATERIAL=GLASS_EPOXY_DMG,ORIENT=RECT
   1.,
  *NSET,NSET=fix1
   1,5,9,13
  *NSET,NSET=fix2
   1,2,3,4
  *NSET,NSET=fix3,GEN
   1,16
  *nset, nset=move1
   4,8,12,16
  *nset, nset=move2
   13,14,15,16
  *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_element
 PART-1-1.EDAMAGE,
*MATERIAL,NAME=GLASS_EPOXY
*DENSITY
 1900.0,
*ELASTIC,TYPE=LAMINA
 53.8e9,17.9e9,0.25,8.96e9,8.96e9,6.88e9
*MATERIAL,NAME=GLASS_EPOXY_DMG
*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,   0.0
 1024e6,1024e6,17.6e6,128e6,31.4e6,59e6, 100.0
** <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
** 1.e-4,1.e-4,1.e-4,1.e-4
*amplitude, name=amp1, definition=smooth step, time=total time
 0.0, 0.0, 0.1, 0.1, 0.2, -0.2
**
**    Step 1:  Tension
**
*STEP
*DYNAMIC, EXPLICIT, scale factor=0.5
 , 0.1
*BOUNDARY, amplitude=amp1
 ASSEMBLY-1.PART-1-1.move1, 1,1, 1.0
 ASSEMBLY-1.PART-1-1.move2, 2,2, 1.0
*BOUNDARY
 ASSEMBLY-1.PART-1-1.fix1,  1,1
 ASSEMBLY-1.PART-1-1.nall,  3,6
 ASSEMBLY-1.PART-1-1.fix2,  2,2
*OUTPUT,FIELD
*ELEMENT OUTPUT
 DMICRT,
 E,
 DAMAGEFT,
 DAMAGEFC,
 DAMAGEMT,
 DAMAGEMC,
 DAMAGESHR,
 STATUS,
*NODE OUTPUT
 U,
***OUTPUT, FIELD
***ELEMENT OUTPUT, ELSET=PART-1-1.EDAMAGE
** STATUS,
*OUTPUT,HISTORY
*node output, nset=ASSEMBLY-1.PART-1-1.move1
 U1,
*node output, nset=ASSEMBLY-1.PART-1-1.move2
 U2,
*ELEMENT OUTPUT,ELSET=PART-1-1.EALL
 ELDMD,
 HSNFTCRT,
 HSNFCCRT,
 HSNMTCRT,
 HSNMCCRT,
 S,
 E,
 DAMAGEFT,
 DAMAGEFC,
 DAMAGEMT,
 DAMAGEMC,
 DAMAGESHR,
 STATUS,
*ENERGY OUTPUT, variable=preselect
***output, history, time interval=0.025
***element output, elset=qa_test_element
** DAMAGEFT, DAMAGEMT, DAMAGEFC, DAMAGEMC
*END STEP
**
**    Step 2:  Compression
**
*STEP
*DYNAMIC, EXPLICIT, scale factor=0.5
 , 0.1
***BOUNDARY, amplitude=amp1
** ASSEMBLY-1.PART-1-1.move1, 1,1, 1.0
** ASSEMBLY-1.PART-1-1.move2, 2,2, 1.0
***BOUNDARY
** ASSEMBLY-1.PART-1-1.fix1,  1,1
** ASSEMBLY-1.PART-1-1.fix1,  3,6
** ASSEMBLY-1.PART-1-1.fix2,  2,2
*END STEP