*HEADING 
 SUPERPLASTIC FORMING OF BOX WITH 
 MEMBRANES - AUTOMATIC LOADING
*RESTART, WRITE, FREQUENCY=30
**
** PLATE DEFINITION
**
*NODE
1,     4.0, 20.0625,  -4.0
23,   26.0, 20.0625,  -4.0 
3201,  4.0, 20.0625, -36.0
3223, 26.0, 20.0625, -36.0 
*NGEN, NSET=EDGE3
1,3201,100
*NGEN, NSET=EDGE5
23,3223,100
*NFILL, NSET=PLATE
EDGE3,EDGE5,22,1
*NSET, NSET=EDGE1, GENERATE
1,23,1
*NSET, NSET=EDGE7, GENERATE
3201,3223,1
*NSET, NSET=CENTER
1,
*ELEMENT, TYPE=M3D4, ELSET=PLATE1
1,1,2,102,101
*ELGEN, ELSET=PLATE1
1,21,1,1,31,100,100
*NSET,ELSET=PLATE1, NSET=NCONT
*ELEMENT, TYPE=M3D4, ELSET=PLATE2
22,22,23,123,122
3101,3101,3102,3202,3201
*ELGEN, ELSET=PLATE2
22,32,100,100
3101,21,1,1
*ELSET,ELSET=PLATE
PLATE1,PLATE2
*MEMBRANE SECTION, ELSET=PLATE, MATERIAL=SUPRAL
.125, 
**
** MATERIAL IS CLOSE TO SUPRAL100 AT 470C
**
*MATERIAL, NAME=SUPRAL
*ELASTIC
10.3E6, .34
*CREEP, LAW=TIME
1.48E-9, 2., 0.
**
** CONTACT surface
**
*NODE, NSET=DIE
 10000, 0.0, 0.0, 0.0
*RIGID BODY, ELSET=ERIGID, REFNODE=10000
*NODE, INPUT=superplasticbox_node.inp
*ELEMENT,TYPE=R3D3,ELSET=ERIGID,
INPUT=superplasticbox_element.inp
*SURFACE,NAME=DIE
 ERIGID,SPOS
*SURFACE,NAME=SLAVES
PLATE,SNEG
*CONTACT PAIR,INTERACTION=DIE_NODE,TYPE=SURFACE,MINIMUM DISTANCE=NO
 SLAVES,DIE
*SURFACE INTERACTION,NAME=DIE_NODE
**
** BOUNDARY AND INITIAL CONDITIONS
**
*BOUNDARY
EDGE1,3
EDGE7,2,3
EDGE3,1
EDGE5,1,2
10000,1,6
*INITIAL CONDITIONS, TYPE=STRESS, UNBALANCED STRESS=STEP
PLATE, 10.0, 10.0
*AMPLITUDE,DEFINITION=SOLUTION DEPENDENT,NAME=AUTO
1.,0.1,1000.
*NSET,NSET=NSELECT
101,121,1001,1101,1507,2509
1,113,1211,1308,1408,1508,1509,1705
**
** STEP 1
**
*STEP, NLGEOM, UNSYMM=YES
*STATIC
2.E-3,1.0,
*DLOAD
PLATE,P,-0.4 
*CONTROLS, PARAMETER=FIELD
1.0e-5,
*CONTACT PRINT, FREQUENCY=100
*CONTACT FILE, FREQUENCY=100, NSET=NSELECT
*EL PRINT, ELSET=PLATE , FREQUENCY=100
S, E
CE, 
SINV, 
*PRINT, CONTACT=YES
*NODE PRINT, NSET=EDGE3, FREQUENCY=100
U, 
*NODE FILE, NSET=EDGE3, FREQUENCY=100
U, 
*OUTPUT,FIELD, Frequency=100,variable=preselect
*NODE OUTPUT,NSET=EDGE3
U, 
*END STEP
**
** STEP 2
**
*STEP, INC=500, NLGEOM, unsymm=yes
*VISCO, CETOL=0.01
0.2, 2000.0
*DLOAD,AMPLITUDE=AUTO
PLATE,P,-0.4
*CREEP STRAIN RATE CONTROL, ELSET=PLATE, 
AMPLITUDE=AUTO 
0.02 , 
*CONTROLS, PARAMETER=FIELD
1.0e-5,
** TO WRITE THE AUTOMATIC SOLUTION CONTROL 
** VARIABLES AMPCU AND RATIO TO THE RESULTS
** FILE EVERY INCREMENT SUCH FILE HAS TO BE
** ACTIVATED
*NODE FILE, NSET=CENTER, FREQUENCY=1
U, 
*OUTPUT,FIELD, OP=ADD,Frequency=1
*NODE OUTPUT,NSET=CENTER
U, 
*END STEP