*HEADING
. . . . . PROJECT
**---------------------------------------------------------
** D E F I N I T I O N ... O F ... N O D E S
**
*NODE 
**
** the following 4 lines define 4 nodes
**
1, 0., 0. 
20, 10., 0.
380, 10.,10.
361, 0.,10.
**
** generate the bottom set of nodes
**
*NGEN, NSET=BOTTOM
**
**(starting node, end node, increment)
**
1,20,1 
**
** generate the top set of nodes
**
*NGEN, NSET=TOP
361,380,1
**
** generate the right set of nodes
**
*NGEN, NSET=RIGHT
20,380,20
**
** generate the left set of nodes
**
*NGEN, NSET=LEFT
1,361,20
**
** generate a set of nodes that includes all nodes by filling
** rows between TOP and BOTTOM 
**
*NFILL,NSET=ALL
**
** starting row, end row, # of rows to generate, increment between 
** nodes of consecutive rows
**
BOTTOM,TOP,18,20
**---------------------------------------------------------
** D E F I N I T I O N.... O F.... E L E M E N T S 
**
*ELEMENT, TYPE=CAX4
1,1,2,22,21 
172,181,182,202,201
**
**(element number, node numbers on the element counterclockwise)
**
** generate element of material A1
**
*ELGEN,ELSET=A1
1,19,1,1,9,20,19
**
** generate element of material A2
**
*ELGEN, ELSET=A2
172,19,1,1,9,20,19
*ELSET, ELSET=EALL
A1,A2
**
** generate all elements
**
**ELGEN,ELSET=ELSET1,ELSET2
**1,19,1,1,19,20,19 
**
** master element,
**# of elements in each row,
**increment in node numbers as you move along the row,
**increment in element numbers as you move along the row,
**# of element rows to be generated,
**increment in node numbers as you move vertically, 
**increment in element numbers as you move vertically.
**---------------------------------------------------------
**
** specify the material used for each element of material A1
**
*MATERIAL,NAME=POWDER1
*ELASTIC, DEPENDENCIES=1
0.1, 0.0001, , 0.5
40000., 0.06, , 0.6
80000., 0.12, , 0.7
120000., 0.18, , 0.8
160000., 0.24, , 0.9
200000., 0.30, , 1.0
*USER DEFINED FIELD
*DEPVAR
1
*PLASTIC
200.
*POROUS METAL PLASTICITY, RELATIVE DENSITY=0.6
1.0,1.0,1.0
*INITIAL CONDITIONS, TYPE=SOLUTION
A1, 0.6
**RELATIVE DENSITY FOR POWDER1
*INITIAL CONDITIONS, TYPE=RELATIVE DENSITY 
ALL, 0.6
*SOLID SECTION, ELSET=A1, MATERIAL=POWDER1
*NODE,nset=mon
9999,0.0,10.0
9998,10.0,0.0
9997,0.0,0.0
**** specify the material used for each element of material A2
**
*MATERIAL,NAME=POWDER2
*ELASTIC, DEPENDENCIES=1
0.2, 0.0001, , 0.5
100000., 0.06, , 0.6
200000., 0.12, , 0.7
300000., 0.18, , 0.8
400000., 0.24, , 0.9
500000., 0.30, , 1.0
*USER DEFINED FIELD
*DEPVAR
1
*PLASTIC
200.
*POROUS METAL PLASTICITY, RELATIVE DENSITY=0.6
1.0,1.0,1.0
*INITIAL CONDITIONS, TYPE=SOLUTION
A2, 0.6
**RELATIVE DENSITY FOR POWDER1
*INITIAL CONDITIONS, TYPE=RELATIVE DENSITY 
ALL, 0.6
*SOLID SECTION, ELSET=A2, MATERIAL=POWDER2
*NODE,nset=mon
9999,0.0,10.0
9998,10.0,0.0
9997,0.0,0.0
**
*SURFACE, TYPE= SEGMENTS, NAME=BOTPUNCH
START,-1.0,0.0
LINE,11.0,0.0
*RIGID BODY, REF NODE=9997,
ANALYTICAL SURFACE=BOTPUNCH
**
*SURFACE, TYPE= SEGMENTS, NAME=SIDE
START,10.0,-1.0
LINE,10.0,11.0
*RIGID BODY, REF NODE=9998,
ANALYTICAL SURFACE=SIDE
**
*SURFACE, TYPE= SEGMENTS, NAME=TOPPUNCH
START,11.0,10.0
LINE,-1.0,10.0
*RIGID BODY, REF NODE=9999,
ANALYTICAL SURFACE=TOPPUNCH
**---------------------------------------------------------
*SURFACE,NAME=POWDER, TYPE=ELEMENTS
EALL
**
** DEFINE SURFACE INTERACTIONS
*CONTACT PAIR, INTERACTION=SMOOTH
POWDER, BOTPUNCH
*CONTACT PAIR, INTERACTION=SMOOTH
POWDER, SIDE
*CONTACT PAIR, INTERACTION=SMOOTH
POWDER, TOPPUNCH
*SURFACE INTERACTION, NAME=SMOOTH
*FRICTION
0.1
**
*** -- STEP 1 -- LOADING
*STEP,NLGEOM, INC=500
*STATIC
0.001,1.75
**
** specify the boundary conditions.
*BOUNDARY
9998, 1, 6, 0.
9997, 1,, 0.
9997, 6,, 0.
9999, 1,, 0.
9999, 6,, 0.
LEFT, 1,,0.
*BOUNDARY, TYPE=VELOCITY
9999, 2,2, -1.0
9997, 2,2, 1.0
**
**
** element set, code of boundary condition,, value
** 
*OUTPUT, FIELD, FREQ=1
*ELEMENT OUTPUT
RD,S,VVF 
*NODE OUTPUT
U, RF
*output, history, freq=1 
*NODE OUTPUT, NSET=MON
U, RF
*END STEP
***
***  -- STEP 2 -- UNLOADING
***
*STEP,NLGEOM
*STATIC
0.001,0.3
**
** specify the boundary conditions.
*BOUNDARY
9998, 1, 6, 0.
9999, 1,, 0.
9999, 6,, 0.
9997, 1,, 0.
9997, 6,, 0.
LEFT, 1,,0.
*BOUNDARY, TYPE=VELOCITY
9999, 2,, +1.5
9997, 2,, 0
**
**
** element set, code of boundary condition,, value
** 
*OUTPUT, FIELD, FREQ=1
*ELEMENT OUTPUT
RD,S,VVF 
*NODE OUTPUT
U, RF
*output, history, freq=1 
*NODE OUTPUT, NSET=MON
U, RF
*END STEP
***
***  -- STEP 3 --- EJECTION
***
*STEP,NLGEOM, INC=500
*STATIC
0.001,1.
**
** specify the boundary conditions.
*BOUNDARY
9998, 1, 6, 0.
9997, 1,, 0.
9997, 6,, 0.
9999, 1,, 0.
9999, 6,, 0.
LEFT, 1,,0.
*BOUNDARY, TYPE=VELOCITY
9999, 2,, +12.
9997, 2,, +12.
**
**
** element set, code of boundary condition,, value
** 
*OUTPUT, FIELD, FREQ=1
*ELEMENT OUTPUT
RD,S,VVF 
*NODE OUTPUT
U, RF
*output, history, freq=1 
*NODE OUTPUT, NSET=MON
U, RF
*END STEP