*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
*NGEN,NSET=END1,LINE=C
*NGEN,NSET=END2,LINE=C
*NFILL,NSET=RIGID
END1,END2,40,17
**
** Plate node definitions
*NODE,NSET=DEFORM
*NGEN,NSET=SHORT
*NGEN,NSET=LONG
*NFILL,NSET=DEFORM
SHORT,LONG,40,17
**
** Rigid body reference node
*NODE,NSET=FIX
**
** Plate element definitions
*ELEMENT, TYPE=S4R, ELSET=SHEET
*ELGEN,ELSET=SHEET
*ELSET,ELSET=SHORT,GENERATE
**
** Rigid surface element definitions
*ELEMENT, TYPE=R3D4, ELSET=CONTACT
*ELGEN,ELSET=CONTACT
**
** 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
**
**
*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.,293.,293.,293.,293.
**
** ---------------------------------------------
**
** 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,APPROACH
**
*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