*HEADING
Multiple step analysis of cantilever beam
*RESTART, WRITE
*NODE, NSET=BUILTIN
1, 0.0,0.0,0.0
41,0.0,0.5,0.0
*NGEN, NSET=BUILTIN
1,41,10
*NCOPY,CHANGE NUMBER=9, OLD SET=BUILTIN, NEW SET=TIP, SHIFT
1.0,0.0,0.0
0.0,0.0,0.0,0.0,1.0,0.0,0.0
*NFILL, NSET=ALL
BUILTIN,TIP,9,1
*NSET, NSET=BUILTIN
1,11,21,31,41
*NSET, NSET=TIP
10,20,30,40,50
*ELEMENT, TYPE=CPS4R, ELSET=EALL
1,1,2,12,11
*ELGEN, ELSET=EALL
1,4,10,1,9,1,4
*SOLID SECTION,ELSET=EALL,MATERIAL=MAT1
*MATERIAL,NAME=MAT1
*ELASTIC
200.E9,0.3
*PLASTIC
380.E6,0.00
580.E6,0.35
780.E6,0.50
*DENSITY
7800.,
*AMPLITUDE, NAME=RAMP
0.0,0.0,1.0,1.0
** step 1
*STEP, NLGEOM
Apply downward load to tip of beam
*DYNAMIC
0.1,1.0
*BOUNDARY
BUILTIN,1,6
*CLOAD, AMPLITUDE=RAMP
TIP,2,-5.E6
*OUTPUT, FIELD
*ELEMENT OUTPUT
S,EE,LE,PE,PEEQ
*NODE OUTPUT
U,
*NODE FILE, NSET=TIP
U,
*EL FILE, ELSET=EALL, FREQ=999
S,EE,LE,PE
*MONITOR, DOF=2, NODE=50
*END STEP
** step 2
*STEP
First springback analysis
*STATIC
0.1,1.0
*BOUNDARY
BUILTIN,1,6
*CLOAD, OP=NEW
*END STEP
** step 3
*STEP, NLGEOM
Apply upward load to tip of beam
*DYNAMIC
0.001,0.01
*BOUNDARY
BUILTIN,1,6
*CLOAD, AMPLITUDE=RAMP
TIP,2,5.E7
*END STEP
** step 4
*STEP, NLGEOM
Second springback analysis
*STATIC
0.1,1.0
*BOUNDARY
BUILTIN,1,6
*CLOAD, OP=NEW
*END STEP
** step 5
*STEP, NLGEOM
Pull tip of beam
*STATIC
0.1,1.0
*BOUNDARY
BUILTIN,1,6
*CLOAD, AMPLITUDE=RAMP,OP=NEW
TIP,1,0.5E7
*END STEP
** step 6
*STEP, NLGEOM
Third springback analysis
*STATIC
0.1,1.0
*BOUNDARY
BUILTIN,1,6
*CLOAD, OP=NEW
*END STEP