*HEADING
/Explicit ADAPTIVE MESHING verification (Units: N,m)
Material Testing for Drucker-Prager Plasticity
Two step analysis to test effect of mesh smoothing
*RESTART,WRITE,NUMBER=10
*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
*NSET,NSET=QA_TEST
5,105,205,
*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
*ELSET,ELSEt=QA_TEST
16,32,13,29,
****************************************
*SOLID SECTION,MATERIAL=MAT1,ELSET=E3D,CONTROL=SECT
0.12, 
*SECTION CONTROLS, NAME=SECT, HOURGLASS=STIFFNESS
*MATERIAL,NAME=MAT1
*ELASTIC
3.0E6,0.3
*DRUCKER PRAGER
10.,1.,10.
*DRUCKER PRAGER HARDENING
.15E5, 0.0
0.6E5,  1.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, 2, 0.0
LEFT, 1, 1, 0.0
RIGHT,1, 1, 0.0
BACK, 3, 3, 0.0
FRONT, 3, 3, 0.0
BOT, 1, 1, 0.0
*STEP
*DYNAMIC,EXPLICIT
,1.1
*BOUNDARY, TYPE=DISPLACEMENT, AMPLITUDE=AMP
1,2,2,0.
2,2,2,-.20
3,2,2,-.40
4,2,2,-.60
5,2,2,-.80
101,2,2,0.
102,2,2,-.20
103,2,2,-.40
104,2,2,-.60
105,2,2,-.80
201,2,2,0.
202,2,2,-.20
203,2,2,-.40
204,2,2,-.60
205,2,2,-.80
*FILE OUTPUT, NUM=1, TIMEMARKS=NO
*NODE FILE
 U,
*EL FILE 
PEEQ,
*ENERGY FILE
ALLIE,ALLKE,ALLWK,ETOTAL,ALLPD
*OUTPUT,FIELD,NUMBER INTERVAL=3,VARIABLE=PRESELECT
*OUTPUT,FIELD,NUMBER INTERVAL=3
*NODE OUTPUT,NSET=QA_TEST
 U
*ELEMENT OUTPUT,ELSET=QA_TEST
 PEEQ
*END STEP
*STEP
*DYNAMIC,EXPLICIT
,1.
*ADAPTIVE MESH,ELSET=E3D,FREQUENCY=1,INITIAL MESH SWEEPS=1
*END STEP