*HEADING
THERMAL FORMING OF AN ALUMINUM SHEET                   
S.I Units, Kg, M, s, N
**
*RESTART, WRITE, FREQUENCY=20, OVERLAY
**
** Rigid body node definitions
*NODE,NSET=RIGID
       1,   -1.800000,    0.000000,    0.000000
      17,    1.800000,    0.000000,    0.000000
     901,    0.000000,    0.000000,   13.000000
     681,   -1.300000,   10.300000,    0.000000
     697,    1.300000,   10.300000,    0.000000
     902,    0.000000,   10.300000,    6.000000
*NGEN,NSET=END1,LINE=C
1,17,1,901
*NGEN,NSET=END2,LINE=C
681,697,1,902
*NFILL,NSET=RIGID
END1,END2,40,17
**
** Plate node definitions
*NODE,NSET=DEFORM
   10001,   -1.000000,   10.150000,   0.1725
   10017,    1.000000,   10.150000,   0.1725
   10681,   -1.500000,    0.150000,   0.1725
   10697,    1.500000,    0.150000,   0.1725
*NGEN,NSET=SHORT
10001,10017
*NGEN,NSET=LONG
10681,10697
*NFILL,NSET=DEFORM
SHORT,LONG,40,17
**
** Rigid body reference node
*NODE,NSET=FIX
99999,       0.,          0.,         0.
**
** Plate element definitions
*ELEMENT, TYPE=S4R, ELSET=SHEET
   10001,   10001,   10002,   10019,   10018
*ELGEN,ELSET=SHEET
   10001,  16,   1,   1,  40,  17,  16
*ELSET,ELSET=SHORT,GENERATE
 10001,10016,1
**
** Rigid surface element definitions
*ELEMENT, TYPE=R3D4, ELSET=CONTACT
       1,   1,   2,  19,  18
*ELGEN,ELSET=CONTACT
       1,  16,   1,   1,  40,  17,  16
**
** Plate 10mm in thickness
*SHELL SECTION, ELSET=SHEET, MATERIAL=ALUMINUM
         0.01,       5
**
*MATERIAL, NAME=ALUMINUM
** 
*ELASTIC
7.919714E10, 0.3, 292.
*DENSITY
2.9e3, 
*EXPANSION
8.367e-6, 
**
*PLASTIC
1.e8,0,293.
1.e3,0,673.
**
*RIGID BODY, ELSET=CONTACT, REF NODE=99999
**
*SURFACE, NAME=M1
CONTACT, SPOS
*SURFACE, NAME=S1
SHEET, SPOS
**
*CONTACT PAIR, INTERACTION=I1
      S1,      M1
*SURFACE INTERACTION, NAME=I1
*FRICTION, SLIP TOLERANCE=0.02
0.1, 
**
** Nodes down the middle of the model on which 
** boundary conditions are placed. 
*NSET,NSET=MIDDLE
   10009,   10026,   10043,   10060,   
   10077,   10094,   10111,   10128,
   10145,   10162,   10179,   10196,   
   10213,   10230,   10247,   10264,
   10281,   10298,   10315,   10332,   
   10349,   10366,   10383,   10400,
   10417,   10434,   10451,   10468,   
   10485,   10502,   10519,   10536,
   10553,   10570,   10587,   10604,   
   10621,   10638,   10655,   10672,
   10689,
**
**
*BOUNDARY
10009,2,2,0.0
MIDDLE,1,1,0.0
99999, 1,6,0.0
**
** Temperature amplitude card for use 
** during forming step
**
*AMPLITUDE, NAME=TEMP_PROF, VALUE=ABSOLUTE
    0,293.,  3600,633. ,  7200,323.
**
** Initial ambient temperature for the plate
**
*INITIAL CONDITIONS, TYPE=TEMPERATURE
DEFORM,293.,0.
**
** ---------------------------------------------
**
** Step 1 - Apply gravity loading together with 
** 1 Bar pressure loading on plate to force it 
** into the die.
*STEP, INC=10000, NLGEOM
*STATIC
       0.1,          1.,       1.E-10
**
** pressure + gravity
**
*DLOAD, OP=NEW
SHEET,GRAV,9.81, 0., 0., -1.
SHEET, P,      1.e5
**
*CONTACT CONTROLS,MASTER=M1, SLAVE=S1,STABILIZE,TANGENT FRACTION=0.
,,0.5       
**
*NODE PRINT, FREQ=0
*EL PRINT, FREQ=0
*CONTACT PRINT, FREQ=999
*NODE FILE,NSET=SHORT,FREQ=999
U, 
*EL FILE, ELSET=SHORT, FREQ=999 
S,E,PE
*OUTPUT,FIELD,FREQ=20
*NODE OUTPUT
U,
*ELEMENT OUTPUT
S,E,PEEQ
*END STEP
**
** ---------------------------------------------
**
** Step 2 - Now perform the forming step. 
** The temperature follows the profile on the 
** *AMPLITUDE card
*STEP, INC=10000, NLGEOM
*STATIC
100.,7200.,1e-1,5.e2
** 
*TEMPERATURE, OP=NEW, AMPLITUDE=TEMP_PROF
DEFORM,1.,1.,1.,1.,1.
**
*END STEP
**
** --------------------------------------------
**
** Step 3 - Now release 1 Bar pressure 
** (but maintain gravity load)
*STEP, INC=10000, NLGEOM
*STATIC
0.1, 1. 
**
*DLOAD, OP=NEW
SHEET,GRAV,9.81, 0., 0., -1.
**
*END STEP