*HEADING
ADAPTIVE MESHING EXAMPLE
3D SPHERICAL INDENTATION
Units - N, mm, sec
*RESTART,TIMEMARKS=YES,WRITE,NUM=1
*NODE,NSET=ALLN
*NSET,NSET=N1
*NSET,NSET=N10
*NSET,NSET=N401
*NSET,NSET=N410
*NFILL,NSET=TOP2D
N1,N401,40,10
*NFILL,NSET=BOT2D
N10,N410,40,10
*NFILL,NSET=HEAD
TOP2D,BOT2D,9,1
*NCOPY, SHIFT, OLD SET=HEAD, NEW SET=TAIL,
CHANGE NUMBER=16400
*NFILL, NSET=NALL
HEAD, TAIL,   40,  410
*ELEMENT,TYPE=C3D8R
*ELGEN,ELSET=BLANK
*NCOPY, SHIFT, OLD SET=TOP2D, NEW SET=TOPLAST,
CHANGE NUMBER=16400
*NFILL, NSET=TOP
TOP2D, TOPLAST,   40,  410
*NCOPY, SHIFT, OLD SET=BOT2D, NEW SET=BOTLAST,
CHANGE NUMBER=16400
*NFILL, NSET=BOT
BOT2D, BOTLAST,   40,  410
*INCLUDE,INPUT=ale_indent_sphelset.inp
*NODE,NSET=NOUT
*ELEMENT,TYPE=MASS,ELSET=PMASS
*MASS,ELSET=PMASS
0.2,
*SOLID SECTION,ELSET=BLANK,MATERIAL=FOAM,CONTROLS=SECT
*SECTION CONTROLS,HOURGLASS=ENHANCED,
KINEMATICS=ORTHOGONAL,NAME=SECT
*MATERIAL,NAME=FOAM
*ELASTIC
 7.5E3, 0.0
*CRUSHABLE FOAM,HARDENING=ISOTROPIC
*CRUSHABLE FOAM HARDENING
*DENSITY
 60.E-9,
*BOUNDARY
BOT,1,3,
100000,1,
100000,3,
100000,4,6,
*AMPLITUDE,NAME=RAMPP,TIME=TOTAL TIME,
DEFINITION=SMOOTH STEP
0.0,0.0,0.06,-250.,
*ELSET,ELSET=SMALL,GEN
*NSET,NSET=SMALL
*SURFACE,TYPE=ELEMENT,NAME=TARGET
UPPER,S5
*SURFACE,TYPE=REVOLUTION,NAME=PUNCH
START,100.,80.
LINE, 100.,0.
CIRCL,0.,-100.,0.,0.
*RIGID BODY,REF NODE=100000,
ANALYTICAL SURFACE =PUNCH
*STEP
*DYNAMIC,EXPLICIT
,0.06
*BOUNDARY,AMPLITUDE=RAMPP
100000,2,2,1.
*SURFACE INTERACTION,NAME=IMP_TARG
*CONTACT PAIR,INTERACTION=IMP_TARG
PUNCH,TARGET
*FILE OUTPUT,NUMBER INTERVAL=1, TIMEMARKS=YES
*EL FILE, ELSET=SMALL
 LE,
*NODE FILE, NSET=NOUT
 U,RF
*ENERGY FILE
*****
*NSET, NSET=QA_TEST_REFN
NOUT,
*ELSET, ELSET=QA_TEST
SMALL,
*NSET, NSET=QA_TEST
SMALL,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST
PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U,
*NODE OUTPUT, NSET=QA_TEST_REFN
RF,
*OUTPUT, HIST, FREQ=9999
*ENERGY OUTPUT, VAR=PRESELECT
******
*ADAPTIVE MESH, ELSET=BLANK, MESH SWEEPS=3
*** output for figures
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT
PEEQ,
*NODE OUTPUT
U,
*END STEP