*HEADING
 ABAQUS/Explicit ADAPTIVE MESHING verification (Units: N,m)
 Material Testing for crushable foam plasticity with
 isotropic hardening
 Two step analysis to test effect of mesh smoothing
*NODE
   1, 0., 0.
   5, 4., 0.
  21, 0., 4.
  25, 4., 4.
 101, 0., 0., 2.
 105, 4., 0., 2.
 121, 0., 4., 2.
 125, 4., 4., 2.
 201, 0., 0., 4.
 205, 4., 0., 4.
 221, 0., 4., 4.
 225, 4., 4., 4.
*NGEN,NSET=BOT
   1,   5, 1
 101, 105, 1
 201, 205, 1
*NGEN,NSET=TOP
  21,  25, 1
 121, 125, 1
 221, 225, 1
*NGEN,NSET=BOT1
 1, 5, 1
*NGEN,NSET=BOT2
 101, 105, 1
*NGEN,NSET=BOT3
 201, 205, 1
*NGEN,NSET=TOP1
 21, 25, 1
*NGEN,NSET=TOP2
 121, 125, 1
*NGEN,NSET=TOP3
 221, 225, 1
*NGEN,NSET=RIGHT
   5,  25, 5
 105, 125, 5
 205, 225, 5
*NGEN,NSET=LEFT
   1,  21, 5
 101, 121, 5
 201, 221, 5
*NFILL,NSET=NALL
 BOT1, TOP1, 4, 5
 BOT2, TOP2, 4, 5
 BOT3, TOP3, 4, 5
*NODE
   7, 1.25, 1.25
   9, 3.25, 1.25
  17, 1.25, 3.25
  19, 3.25, 3.25
 107, 1.25, 1.25, 2.
 109, 3.25, 1.25, 2.
 117, 1.25, 3.25, 2.
 119, 3.25, 3.25, 2.
 207, 1.25, 1.25, 4.
 209, 3.25, 1.25, 4.
 217, 1.25, 3.25, 4.
 219, 3.25, 3.25, 4.
*NSET,NSET=FRONT,GEN
 1, 25, 1
*NSET,NSET=BACK,GEN
 201, 225, 1
*ELEMENT,TYPE=C3D8R
 1, 1, 2, 7, 6, 101, 102, 107, 106
*ELGEN,ELSET=E3D
 1, 4, 1, 1, 4, 5, 4, 2, 100, 16
*NSET,NSET=QA_TEST_NALL
 NALL,
*ELSET,ELSEt=QA_TEST_EALL
 E3D,
**
*SOLID SECTION,MATERIAL=MAT1,ELSET=E3D,CONTROL=SECT
 0.12,
*SECTION CONTROLS,NAME=SECT,HOURGLASS=STIFFNESS
*MATERIAL,NAME=MAT1
*ELASTIC
 3.E6, 0.
*CRUSHABLE FOAM,HARDENING=ISOTROPIC
 1.1, 0.2983
** plastic Poisson's ratio is chosen in such a way 
** that the plastic flow is associative.
*CRUSHABLE FOAM HARDENING
 220000.0,  0.0
 246512.1,  0.1
 272938.0,  0.2
 299019.2,  0.3
 324549.1,  0.4
 349345.3,  0.5
 373259.8,  0.6
 396170.9,  0.7
 418011.1,  0.8
 438721.0,  0.9
 458268.2,  1.0
 493843.0,  1.2
 524836.3,  1.4
 551528.4,  1.6
 574307.3,  1.8
 593594.8,  2.0
 629359.6,  2.5
 651994.7,  3.0
 683338.7,  5.0
 688332.2, 10.0
*DENSITY
 500.,
**
*AMPLITUDE,TIME=TOTAL TIME,NAME=AMP,DEFINITION=SMOOTH STEP
 0.0, 0.0, 1.0, 1.0, 4.0, 1.0
*BOUNDARY
 TOP,   2
 LEFT,  1
 RIGHT, 1
 BACK,  3
 FRONT, 3
 BOT,   1
**
*STEP
*DYNAMIC,EXPLICIT
 ,1.1
*BOUNDARY,TYPE=DISPLACEMENT,AMPLITUDE=AMP
   1, 2, 2,  0.
   2, 2, 2, -0.20
   3, 2, 2, -0.40
   4, 2, 2, -0.60
   5, 2, 2, -0.80
 101, 2, 2,  0.
 102, 2, 2, -0.20
 103, 2, 2, -0.40
 104, 2, 2, -0.60
 105, 2, 2, -0.80
 201, 2, 2,  0.
 202, 2, 2, -0.20
 203, 2, 2, -0.40
 204, 2, 2, -0.60
 205, 2, 2, -0.80
*FILE OUTPUT,NUM=1,TIMEMARKS=YES
*NODE FILE
 U,
*EL FILE 
 PEEQ,
*OUTPUT,FIELD,NUMBER INTERVAL=3,VARIABLE=PRESELECT,TIMEMARKS=YES
*OUTPUT,FIELD,NUMBER INTERVAL=1
*NODE OUTPUT,NSET=QA_TEST_NALL
 U,
*ELEMENT OUTPUT,ELSET=QA_TEST_EALL
 PEEQ,
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIME INTERVAL=1.1
**
*END STEP
*STEP
*DYNAMIC,EXPLICIT
 ,1.
*ADAPTIVE MESH,ELSET=E3D,FREQUENCY=1,INITIAL MESH SWEEPS=1
*END STEP