A.8 Deep drawing of a can bottom

draw_freq_std.inp

*HEADING
Deep drawing of a can bottom -- frequency analysis
ABAQUS/Standard
SI units (kg, m, s, N)
**
** BLANK
**
*NODE,NSET=BLANK
501,.000,.000
601,.045,.000
*NGEN,NSET=BLANK
501,601,1
*NSET,NSET=BSYM
501,
*NSET,NSET=BEND
601,
*ELEMENT,TYPE=SAX1,ELSET=BLANK
501,501,503
*ELGEN,ELSET=BLANK
501,50,2,2
*SHELL SECTION,ELSET=BLANK,MATERIAL=STEEL,
SECTION INTEGRATION=GAUSS
.0005,5
*MATERIAL,NAME=STEEL
*ELASTIC
210.0E9, 0.3
*DENSITY
7800.,
*BOUNDARY
BSYM,1
BSYM,6
BEND,2
*STEP
*FREQUENCY, eigensolver=lanczos
1, , 
*RESTART,WRITE
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*NODE FILE,FREQ=999
U,
*OUTPUT, FIELD
*NODE OUTPUT, NSET=BLANK
U,
*END STEP

draw_stage1_attempt1.inp

*HEADING
Deep drawing of a can bottom -- stage 1, attempt 1
SI units (kg, m, s, N)
**
*NODE,NSET=REFP1
1001,2.6e-2,6.25e-3
*NODE,NSET=REFD1
1101,3.66e-2,-4.25e-3
**
** BLANK
**
*NODE,NSET=BLANK
501,.000,.000
601,.045,.000
*NGEN,NSET=BLANK
501,601,1
*ELEMENT,TYPE=SAX1,ELSET=BLANK
501,501,503
*ELGEN,ELSET=BLANK
501,50,2,2
*SHELL SECTION,ELSET=BLANK,MATERIAL=STEEL,
SECTION INTEGRATION=GAUSS
.0005,5
*NSET,NSET=BSYM
501,
*NSET,NSET=NHIST
521,601
**
** PUNCH1 -- PUNCH FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=PUNCH1, TYPE=SEGMENTS
START,0.032,0.03025
LINE,0.032,0.00625
CIRCL,0.026,0.00025,0.026,0.00625
LINE,0.,0.00025
*RIGID BODY, REF NODE=1001, ANALYTICAL SURFACE=PUNCH1
**
** DIE1 -- DIE FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=DIE1, TYPE=SEGMENTS
START,0.0326,-0.02825
LINE,0.0326,-0.00425
CIRCL,0.0366,-0.00025,0.0366,-0.00425
LINE,0.050,-0.00025
*RIGID BODY, REF NODE=1101, ANALYTICAL SURFACE=DIE1
**
** SURFACE  DEFINITIONS
**
*SURFACE, NAME=BLANK_TOP
BLANK,SPOS
*SURFACE, NAME=BLANK_BOT
BLANK,SNEG
**
** BLANK MATERIAL
**
*MATERIAL,NAME=STEEL
*ELASTIC
210.0E9,0.3
*DENSITY
7800.,
*PLASTIC
  0.912E8, 0.0
  0.131E9, 0.159E-02
  0.171E9, 0.649E-02
  0.211E9, 0.177E-01
  0.251E9, 0.395E-01
  0.291E9, 0.776E-01
  0.331E9, 0.139
  0.391E9, 0.295
**
** AMPLITUDE DEFINITIONS
**
*AMPLITUDE,NAME=FORM1
0.,0.,.033,1.
**
** INITIAL FIXED BOUNDARY CONDITIONS
**
*BOUNDARY
REFP1,1
REFP1,6
REFD1,1,6
BSYM,1
BSYM,6
**
** STEP 1 -- FIRST FORMING STAGE
**
*STEP
*DYNAMIC,EXPLICIT
,.033
**
** Contact Pair Definition
**
*CONTACT PAIR,INTERACTION=INTER
PUNCH1,BLANK_TOP
DIE1,BLANK_BOT
*SURFACE INTERACTION,NAME=INTER
*FRICTION
.1,
*BOUNDARY,TYPE=DISPLACEMENT,AMPLITUDE=FORM1
REFP1,2,2,-.015
*RESTART,WRITE,NUMBER INTERVAL=10
*OUTPUT, FIELD, NUMBER=4
*ELEMENT OUTPUT, ELSET=BLANK
STH, SF, S
*NODE OUTPUT
U,RF 
*OUTPUT, HISTORY, TIME INTERVAL=0.033
*NODE OUTPUT, NSET=NHIST 
U,V
*ENERGY OUTPUT
ALLKE, ALLWK, ALLVD, ALLAE
**
** Output for ABAQUS QA purposes
**
*NSET,NSET=QA_TEST
NHIST,
*OUTPUT, FIELD, NUMBER=1
*NODE OUTPUT, NSET=QA_TEST
U,
*FILE OUTPUT,NUMBER INTERVAL=1
*ENERGY FILE
*NODE FILE
U,RF
*EL FILE
S,
*END STEP

draw_spring1_std.inp

*HEADING
Deep drawing of a can bottom -- first springback
ABAQUS/Standard springback following
draw_stage1_attempt3_sqrt10.inp
SI units (kg, m, s, N)
*IMPORT,STEP=1,INTERVAL=10,STATE=YES,UPDATE=NO
BLANK,
*IMPORT ELSET
BLANK,
*IMPORT NSET
BSYM,NHIST,BLANK
*NSET, NSET=QA_TEST
BLANK,
*ELSET, ELSET=QA_TEST
BLANK,
*BOUNDARY
BSYM,1
BSYM,6
*RESTART,WRITE,FREQ=999
*STEP,NLGEOM
*STATIC
.1,1.
*BOUNDARY,FIXED
BSYM,2,2
*OUTPUT, FIELD, FREQUENCY=10 
*ELEMENT OUTPUT, ELSET=BLANK
STH, SF, S
*NODE OUTPUT
U, 
**
** Output for ABAQUS QA purposes
**
*NODE FILE,FREQ=999
U,
*EL FILE,FREQ=999
S,
*OUTPUT, FIELD,  FREQ=999
*NODE OUTPUT, NSET=QA_TEST
U,
*ELEMENT OUTPUT, ELSET=QA_TEST
S, 
*OUTPUT, HIST,  FREQ=999
*ENERGY OUTPUT, VAR=ALL
**
*END STEP

draw_stage1_attempt2.inp

*HEADING
Deep drawing of a can bottom -- stage 1, attempt 2
SI units (kg, m, s, N)
**
*NODE,NSET=REFP1
1001,2.6e-2,6.25e-3
*NODE,NSET=REFD1
1101,3.66e-2,-4.25e-3
**
** BLANK
**
*NODE,NSET=BLANK
501,.000,.000
601,.045,.000
*NGEN,NSET=BLANK
501,601,1
*ELEMENT,TYPE=SAX1,ELSET=BLANK
501,501,503
*ELGEN,ELSET=BLANK
501,50,2,2
*SHELL SECTION,ELSET=BLANK,MATERIAL=STEEL,SECTION INTEGRATION=GAUSS
.0005,5
*NSET,NSET=BSYM
501,
*NSET,NSET=NHIST
521,601
**
** PUNCH1 -- PUNCH FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=PUNCH1, TYPE=SEGMENTS
START,0.032,0.03025
LINE,0.032,0.00625
CIRCL,0.026,0.00025,0.026,0.00625
LINE,0.,0.00025
*RIGID BODY, REF NODE=1001, ANALYTICAL SURFACE=PUNCH1
**
** DIE1 -- DIE FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=DIE1, TYPE=SEGMENTS
START,0.0326,-0.02825
LINE,0.0326,-0.00425
CIRCL,0.0366,-0.00025,0.0366,-0.00425
LINE,0.050,-0.00025
*RIGID BODY, REF NODE=1101, ANALYTICAL SURFACE=DIE1
**
** SURFACE  DEFINITIONS
**
*SURFACE, NAME=BLANK_TOP
BLANK,SPOS
*SURFACE, NAME=BLANK_BOT
BLANK,SNEG
**
** BLANK MATERIAL
**
*MATERIAL,NAME=STEEL
*ELASTIC
210.0E9,0.3
*DENSITY
7800.,
*PLASTIC
  0.912E8, 0.0
  0.131E9, 0.159E-02
  0.171E9, 0.649E-02
  0.211E9, 0.177E-01
  0.251E9, 0.395E-01
  0.291E9, 0.776E-01
  0.331E9, 0.139
  0.391E9, 0.295
**
** AMPLITUDE DEFINITIONS
**
*AMPLITUDE,NAME=FORM1, DEFINITION=SMOOTH STEP
0.,0.,.033,1.
**
** INITIAL FIXED BOUNDARY CONDITIONS
**
*BOUNDARY
REFP1,1
REFP1,6
REFD1,1,6
BSYM,1
BSYM,6
**
** STEP 1 -- FIRST FORMING STAGE
**
*STEP
*DYNAMIC,EXPLICIT
,.033
**
** Contact Pair Definition
**
*CONTACT PAIR,INTERACTION=INTER
PUNCH1,BLANK_TOP
DIE1,BLANK_BOT
*SURFACE INTERACTION,NAME=INTER
*FRICTION
.1,
*BOUNDARY,TYPE=DISPLACEMENT,AMPLITUDE=FORM1
REFP1,2,2,-.015
*RESTART,WRITE,NUMBER INTERVAL=10
*OUTPUT, HISTORY, TIME INTERVAL=0.033 
*NODE OUTPUT, NSET=NHIST 
U,V
*ENERGY OUTPUT
ALLKE, ALLIE, ALLWK, ALLVD, ALLAE, ALLPD, ETOTAL
*OUTPUT, FIELD, NUMBER=4
*ELEMENT OUTPUT, ELSET=BLANK
STH, SF, S
*NODE OUTPUT
U,RF 
**
** Output for ABAQUS QA purposes
**
*NSET,NSET=QA_TEST
NHIST,
*OUTPUT, FIELD, NUMBER=1
*NODE OUTPUT, NSET=QA_TEST
U,
*FILE OUTPUT,NUMBER INTERVAL=1
*ENERGY FILE
*NODE FILE
U,RF
*EL FILE
S,
*END STEP

draw_stage1_attempt3.inp

*HEADING
Deep drawing of a can bottom -- stage 1, attempt 3
SI units (kg, m, s, N)
**
*NODE,NSET=REFP1
1001,2.6e-2,6.25e-3
*NODE,NSET=REFD1
1101,3.66e-2,-4.25e-3
**
** BLANK
**
*NODE,NSET=BLANK
501,.000,.000
601,.045,.000
*NGEN,NSET=BLANK
501,601,1
*ELEMENT,TYPE=SAX1,ELSET=BLANK
501,501,503
*ELGEN,ELSET=BLANK
501,50,2,2
*SHELL SECTION,ELSET=BLANK,MATERIAL=STEEL,SECTION INTEGRATION=GAUSS
.0005,5
*NSET,NSET=BSYM
501,
*NSET,NSET=NHIST
521,601
**
** PUNCH1 -- PUNCH FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=PUNCH1, TYPE=SEGMENTS
START,0.032,0.03025
LINE,0.032,0.00625
CIRCL,0.026,0.00025,0.026,0.00625
LINE,0.,0.00025
*RIGID BODY, REF NODE=1001, ANALYTICAL SURFACE=PUNCH1
**
** DIE1 -- DIE FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=DIE1, TYPE=SEGMENTS
START,0.0326,-0.02825
LINE,0.0326,-0.00425
CIRCL,0.0366,-0.00025,0.0366,-0.00425
LINE,0.050,-0.00025
*RIGID BODY, REF NODE=1101, ANALYTICAL SURFACE=DIE1
**
** SURFACE  DEFINITIONS
**
*SURFACE, NAME=BLANK_TOP
BLANK,SPOS
*SURFACE, NAME=BLANK_BOT
BLANK,SNEG
**
** BLANK MATERIAL
**
*MATERIAL,NAME=STEEL
*ELASTIC
210.0E9,0.3
*DENSITY
7800.,
*DAMPING,BETA=1.E-7
*PLASTIC
  0.912E8, 0.0
  0.131E9, 0.159E-02
  0.171E9, 0.649E-02
  0.211E9, 0.177E-01
  0.251E9, 0.395E-01
  0.291E9, 0.776E-01
  0.331E9, 0.139
  0.391E9, 0.295
**
** AMPLITUDE DEFINITIONS
**
*AMPLITUDE,NAME=FORM1, DEFINITION=SMOOTH STEP
0.,0.,.033,1.
**
** INITIAL FIXED BOUNDARY CONDITIONS
**
*BOUNDARY
REFP1,1
REFP1,6
REFD1,1,6
BSYM,1
BSYM,6
**
** STEP 1 -- FIRST FORMING STAGE
**
*STEP
*DYNAMIC,EXPLICIT
,.033
**
** Contact Pair Definition
**
*CONTACT PAIR,INTERACTION=INTER
PUNCH1,BLANK_TOP
DIE1,BLANK_BOT
*SURFACE INTERACTION,NAME=INTER
*FRICTION
.1,
*BOUNDARY,TYPE=DISPLACEMENT,AMPLITUDE=FORM1
REFP1,2,2,-.015
*RESTART,WRITE,NUMBER INTERVAL=10
*OUTPUT, FIELD, NUMBER=4
*ELEMENT OUTPUT, ELSET=BLANK
STH, SF, S
*NODE OUTPUT
U,RF 
*OUTPUT, HISTORY, TIME INTERVAL=0.03
*ENERGY OUTPUT
ALLKE, ALLWK, ALLVD, ALLAE
**
** Output for ABAQUS QA purposes
**
*NSET,NSET=QA_TEST
NHIST,
*OUTPUT, FIELD, NUMBER=1
*NODE OUTPUT, NSET=QA_TEST
U,
*FILE OUTPUT,NUMBER INTERVAL=1
*ENERGY FILE
*NODE FILE
U,RF
*EL FILE
S,
*END STEP

draw_stage1_attempt3_sqrt10.inp

*HEADING
Deep drawing of a can bottom -- stage 1, attempt 3
Mass scaling factor = 10 (speedup of sqrt10)
SI units (kg, m, s, N)
**
*NODE,NSET=REFP1
1001,2.6e-2,6.25e-3
*NODE,NSET=REFD1
1101,3.66e-2,-4.25e-3
**
** BLANK
**
*NODE,NSET=BLANK
501,.000,.000
601,.045,.000
*NGEN,NSET=BLANK
501,601,1
*ELEMENT,TYPE=SAX1,ELSET=BLANK
501,501,503
*ELGEN,ELSET=BLANK
501,50,2,2
*SECTION CONTROLS, NAME=ENHANCED,HOURGLASS=ENHANCED
*SHELL SECTION,ELSET=BLANK,MATERIAL=STEEL,SECTION INTEGRATION=GAUSS,
CONTROLS=ENHANCED
.0005,5
*NSET,NSET=BSYM
501,
*NSET,NSET=NHIST
521,601
**
** PUNCH1 -- PUNCH FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=PUNCH1, TYPE=SEGMENTS
START,0.032,0.03025
LINE,0.032,0.00625
CIRCL,0.026,0.00025,0.026,0.00625
LINE,0.,0.00025
*RIGID BODY, REF NODE=1001, ANALYTICAL SURFACE=PUNCH1
**
** DIE1 -- DIE FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=DIE1, TYPE=SEGMENTS
START,0.0326,-0.02825
LINE,0.0326,-0.00425
CIRCL,0.0366,-0.00025,0.0366,-0.00425
LINE,0.050,-0.00025
*RIGID BODY, REF NODE=1101, ANALYTICAL SURFACE=DIE1
**
** SURFACE  DEFINITIONS
**
*SURFACE, NAME=BLANK_TOP
BLANK,SPOS
*SURFACE, NAME=BLANK_BOT
BLANK,SNEG
**
** BLANK MATERIAL
**
*MATERIAL,NAME=STEEL
*ELASTIC
210.0E9,0.3
*DENSITY
7800.,
*DAMPING,BETA=1.E-7
*PLASTIC
  0.912E8, 0.0
  0.131E9, 0.159E-02
  0.171E9, 0.649E-02
  0.211E9, 0.177E-01
  0.251E9, 0.395E-01
  0.291E9, 0.776E-01
  0.331E9, 0.139
  0.391E9, 0.295
**
** AMPLITUDE DEFINITIONS
**
*AMPLITUDE,NAME=FORM1, DEFINITION=SMOOTH STEP
0.,0.,.033,1.
**
** INITIAL FIXED BOUNDARY CONDITIONS
**
*BOUNDARY
REFP1,1
REFP1,6
REFD1,1,6
BSYM,1
BSYM,6
**
** STEP 1 -- FIRST FORMING STAGE
**
**
*ELSET,ELSET=QA_TEST_ALLELEMS
BLANK,
**
*NSET, NSET=QA_TEST_ALLNODES
REFP1, REFD1, BLANK, BLANK
BSYM, NHIST
**
*STEP
*DYNAMIC,EXPLICIT
,.033
*FIXED MASS SCALING, FACTOR=10
**
** Contact Pair Definition
**
*CONTACT PAIR,INTERACTION=INTER
PUNCH1,BLANK_TOP
DIE1,BLANK_BOT
*SURFACE INTERACTION,NAME=INTER
*FRICTION
.1,
*BOUNDARY,TYPE=DISPLACEMENT,AMPLITUDE=FORM1
REFP1,2,2,-.015
*RESTART,WRITE,NUMBER INTERVAL=10
*OUTPUT, HISTORY, TIME INTERVAL=1.E-5 
*NODE OUTPUT, NSET=NHIST 
U,V
*ENERGY OUTPUT
ALLKE, ALLIE, ALLWK, ALLVD, ALLAE, ALLPD, ETOTAL
*OUTPUT, FIELD, NUMBER=4
*ELEMENT OUTPUT, ELSET=BLANK
STH, SF, S
*NODE OUTPUT
U,RF 
**
** Output for ABAQUS QA purposes
**
*FILE OUTPUT,NUMBER INTERVAL=1
*ENERGY FILE
*NODE FILE
U,RF
*EL FILE
S,
*OUTPUT, FIELD, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
U,
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.033
***OUTPUT,HISTORY,VARIABLE=PRESELECT
*END STEP

draw_stage1_attempt3_10.inp

*HEADING
Deep drawing of a can bottom -- stage 1, attempt 3
Mass scaling factor = 100 (speedup of 10)
SI units (kg, m, s, N)
**
*NODE,NSET=REFP1
1001,2.6e-2,6.25e-3
*NODE,NSET=REFD1
1101,3.66e-2,-4.25e-3
**
** BLANK
**
*NODE,NSET=BLANK
501,.000,.000
601,.045,.000
*NGEN,NSET=BLANK
501,601,1
*ELEMENT,TYPE=SAX1,ELSET=BLANK
501,501,503
*ELGEN,ELSET=BLANK
501,50,2,2
*SHELL SECTION,ELSET=BLANK,MATERIAL=STEEL,SECTION INTEGRATION=GAUSS
.0005,5
*NSET,NSET=BSYM
501,
*NSET,NSET=NHIST
521,601
**
** PUNCH1 -- PUNCH FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=PUNCH1, TYPE=SEGMENTS
START,0.032,0.03025
LINE,0.032,0.00625
CIRCL,0.026,0.00025,0.026,0.00625
LINE,0.,0.00025
*RIGID BODY, REF NODE=1001, ANALYTICAL SURFACE=PUNCH1
**
** DIE1 -- DIE FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=DIE1, TYPE=SEGMENTS
START,0.0326,-0.02825
LINE,0.0326,-0.00425
CIRCL,0.0366,-0.00025,0.0366,-0.00425
LINE,0.050,-0.00025
*RIGID BODY, REF NODE=1101, ANALYTICAL SURFACE=DIE1
**
** SURFACE  DEFINITIONS
**
*SURFACE, NAME=BLANK_TOP
BLANK,SPOS
*SURFACE, NAME=BLANK_BOT
BLANK,SNEG
**
** BLANK MATERIAL
**
*MATERIAL,NAME=STEEL
*ELASTIC
210.0E9,0.3
*DENSITY
7800.,
*DAMPING,BETA=1.E-7
*PLASTIC
  0.912E8, 0.0
  0.131E9, 0.159E-02
  0.171E9, 0.649E-02
  0.211E9, 0.177E-01
  0.251E9, 0.395E-01
  0.291E9, 0.776E-01
  0.331E9, 0.139
  0.391E9, 0.295
**
** AMPLITUDE DEFINITIONS
**
*AMPLITUDE,NAME=FORM1, DEFINITION=SMOOTH STEP
0.,0.,.033,1.
**
** INITIAL FIXED BOUNDARY CONDITIONS
**
*BOUNDARY
REFP1,1
REFP1,6
REFD1,1,6
BSYM,1
BSYM,6
**
** STEP 1 -- FIRST FORMING STAGE
**
**
*ELSET,ELSET=QA_TEST_ALLELEMS
BLANK,
**
*NSET, NSET=QA_TEST_ALLNODES
REFP1, REFD1, BLANK, BLANK
BSYM, NHIST
**
*STEP
*DYNAMIC,EXPLICIT
,.033
*FIXED MASS SCALING, FACTOR=100
**
** Contact Pair Definition
**
*CONTACT PAIR,INTERACTION=INTER
PUNCH1,BLANK_TOP
DIE1,BLANK_BOT
*SURFACE INTERACTION,NAME=INTER
*FRICTION
.1,
*BOUNDARY,TYPE=DISPLACEMENT,AMPLITUDE=FORM1
REFP1,2,2,-.015
*RESTART,WRITE,NUMBER INTERVAL=10
*OUTPUT, HISTORY, TIME INTERVAL=1.E-5 
*NODE OUTPUT, NSET=NHIST 
U,V
*ENERGY OUTPUT
ALLKE, ALLIE, ALLWK, ALLVD, ALLAE, ALLPD, ETOTAL
*OUTPUT, FIELD, NUMBER=4
*ELEMENT OUTPUT, ELSET=BLANK
STH, SF, S
*NODE OUTPUT
U,RF 
**
** Output for ABAQUS QA purposes
**
*FILE OUTPUT,NUMBER INTERVAL=1
*ENERGY FILE
*NODE FILE
U,RF
*EL FILE
S,
*OUTPUT, FIELD, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
U,
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.033
***OUTPUT,HISTORY,VARIABLE=PRESELECT
*END STEP

draw_stage1_attempt3_30.inp

*HEADING
Deep drawing of a can bottom -- stage 1, attempt 3
Mass scaling factor = 900 (speedup of 30)
SI units (kg, m, s, N)
**
*NODE,NSET=REFP1
1001,2.6e-2,6.25e-3
*NODE,NSET=REFD1
1101,3.66e-2,-4.25e-3
**
** BLANK
**
*NODE,NSET=BLANK
501,.000,.000
601,.045,.000
*NGEN,NSET=BLANK
501,601,1
*ELEMENT,TYPE=SAX1,ELSET=BLANK
501,501,503
*ELGEN,ELSET=BLANK
501,50,2,2
*SHELL SECTION,ELSET=BLANK,MATERIAL=STEEL,SECTION INTEGRATION=GAUSS
.0005,5
*NSET,NSET=BSYM
501,
*NSET,NSET=NHIST
521,601
**
** PUNCH1 -- PUNCH FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=PUNCH1, TYPE=SEGMENTS
START,0.032,0.03025
LINE,0.032,0.00625
CIRCL,0.026,0.00025,0.026,0.00625
LINE,0.,0.00025
*RIGID BODY, REF NODE=1001, ANALYTICAL SURFACE=PUNCH1
**
** DIE1 -- DIE FOR FORMING STAGE 1 (Analytical Rigid Surface)
**
*SURFACE, NAME=DIE1, TYPE=SEGMENTS
START,0.0326,-0.02825
LINE,0.0326,-0.00425
CIRCL,0.0366,-0.00025,0.0366,-0.00425
LINE,0.050,-0.00025
*RIGID BODY, REF NODE=1101, ANALYTICAL SURFACE=DIE1
**
** SURFACE  DEFINITIONS
**
*SURFACE, NAME=BLANK_TOP
BLANK,SPOS
*SURFACE, NAME=BLANK_BOT
BLANK,SNEG
**
** BLANK MATERIAL
**
*MATERIAL,NAME=STEEL
*ELASTIC
210.0E9,0.3
*DENSITY
7800.,
*DAMPING,BETA=1.E-7
*PLASTIC
  0.912E8, 0.0
  0.131E9, 0.159E-02
  0.171E9, 0.649E-02
  0.211E9, 0.177E-01
  0.251E9, 0.395E-01
  0.291E9, 0.776E-01
  0.331E9, 0.139
  0.391E9, 0.295
**
** AMPLITUDE DEFINITIONS
**
*AMPLITUDE,NAME=FORM1, DEFINITION=SMOOTH STEP
0.,0.,.033,1.
**
** INITIAL FIXED BOUNDARY CONDITIONS
**
*BOUNDARY
REFP1,1
REFP1,6
REFD1,1,6
BSYM,1
BSYM,6
**
** STEP 1 -- FIRST FORMING STAGE
**
**
*ELSET,ELSET=QA_TEST_ALLELEMS
BLANK,
**
*NSET, NSET=QA_TEST_ALLNODES
REFP1, REFD1, BLANK, BLANK
BSYM, NHIST
**
*STEP
*DYNAMIC,EXPLICIT
,.033
*FIXED MASS SCALING, FACTOR=900
**
** Contact Pair Definition
**
*CONTACT PAIR,INTERACTION=INTER
PUNCH1,BLANK_TOP
DIE1,BLANK_BOT
*SURFACE INTERACTION,NAME=INTER
*FRICTION
.1,
*BOUNDARY,TYPE=DISPLACEMENT,AMPLITUDE=FORM1
REFP1,2,2,-.015
*RESTART,WRITE,NUMBER INTERVAL=10
*OUTPUT, HISTORY, TIME INTERVAL=1.E-5 
*NODE OUTPUT, NSET=NHIST 
U,V
*ENERGY OUTPUT
ALLKE, ALLIE, ALLWK, ALLVD, ALLAE, ALLPD, ETOTAL
*OUTPUT, FIELD, NUMBER=4
*ELEMENT OUTPUT, ELSET=BLANK
STH, SF, S
*NODE OUTPUT
U,RF 
**
** Output for ABAQUS QA purposes
**
*FILE OUTPUT,NUMBER INTERVAL=1
*ENERGY FILE
*NODE FILE
U,RF
*EL FILE
S,
*OUTPUT, FIELD, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
U,
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.033
***OUTPUT,HISTORY,VARIABLE=PRESELECT
*END STEP

draw_stage2.inp

*HEADING
Deep drawing of a can bottom -- stage 2
ABAQUS/Explicit forming stage 2 following draw_spring1_std.inp
SI units (kg, m, s, N)
*IMPORT,STATE=YES,UPDATE=NO
BLANK,
*IMPORT NSET
BSYM,NHIST,BLANK
*NSET, NSET=QA_TEST
BLANK,
*ELSET, ELSET=QA_TEST
BLANK,
**
** PUNCH2 -- PUNCH FOR FORMING STAGE 2
**
*NODE,NSET=PUNCH2
201,.023875,.030
202,.023875,.002
232,.025875,.000
262,.031875,.006
263,.031875,.030
298,.025875,.002
299,.025875,.006
1201,.023875,.030
*NGEN,NSET=PUNCH2,LINE=C
202,232,1,298
232,262,1,299
*ELEMENT,TYPE=RAX2,ELSET=PUNCH2
201,201,202
*ELGEN,ELSET=PUNCH2
201,62,1,1
*RIGID BODY,ELSET=PUNCH2,REF NODE=1201
*NSET,NSET=REFP2
1201,
**
** DIE2 -- DIE FOR FORMING STAGE 2
**
*NODE,NSET=DIE2
301,.026,-.010
302,.026, .000
332,.0325, .0065
333,.050, .0065
399,.026, .0065
1301,.026,-.010
1302,.026,-.030
*NGEN,NSET=DIE2,LINE=C
302,332,1,399
*ELEMENT,TYPE=RAX2,ELSET=DIE2
301,301,302
*ELGEN,ELSET=DIE2
301,32,1,1
*RIGID BODY,ELSET=DIE2,REF NODE=1301
*NSET,NSET=REFD2
1301,
**
** DOMER -- DOMER DIE FOR FORMING STAGE 2
**
*NODE,NSET=DOMER
401, .000,   .008
431, .022,   .000
432, .022,  -.010
499, .000,  -.02625
1401,.000,   .008
*NGEN,NSET=DOMER,LINE=C
401,431,1,499
*ELEMENT,TYPE=RAX2,ELSET=DOMER
401,401,402
*ELGEN,ELSET=DOMER
401,31,1,1
*RIGID BODY,ELSET=DOMER,REF NODE=1401
*NSET,NSET=REFDOMER
1401,
**
** Surface Definitions
**
*SURFACE, NAME=BLANK_TOP
BLANK,SPOS
*SURFACE, NAME=BLANK_BOT
BLANK,SNEG
*SURFACE, NAME=PUNCH2
PUNCH2,SNEG
*SURFACE, NAME=DIE2
DIE2,SPOS
*SURFACE, NAME=DOMER
DOMER,SPOS
**
** SPRING AND DASHPOT FOR FORMING STAGE 2
**
*ELEMENT,TYPE=SPRINGA,ELSET=SPRING
1301,1301,1302
*NSET,NSET=SPEND
1302,
*SPRING,ELSET=SPRING

6.3E6,
*ELEMENT,TYPE=DASHPOTA,ELSET=DASHPOT
2301,1301,1302
*DASHPOT,ELSET=DASHPOT

1.E4,
*NSET,NSET=REF
REFP2,REFD2,REFDOMER,SPEND
**
** MASS FOR SPRING IN FORMING STAGE 2
**
*ELEMENT,TYPE=MASS,ELSET=MASS
1302,1301
*MASS,ELSET=MASS
10.,
**
** AMPLITUDE DEFINITIONS
**
*AMPLITUDE,NAME=RIGID,DEFINITION=SMOOTH STEP
0.,0.,1.E-6,1.
*AMPLITUDE,NAME=FORM2,DEFINITION=SMOOTH STEP
0.,0.,.033,1.
**
** STEP 4 -- POSITION RIGID BODIES FOR SECOND PUNCHING STAGE
**
*STEP
*DYNAMIC,EXPLICIT
,1.E-6
*BULK VISCOSITY
0.02, 0.4
*RESTART,WRITE,NUM=1
*BOUNDARY,TYPE=VELOCITY
REFDOMER,1,1,0.
REFDOMER,6,6,0.
REFD2,1,1,0.
REFD2,6,6,0.
SPEND,1,1,0.
SPEND,3,3,0.
REFP2,1,1,0.
REFP2,6,6,0.
*BOUNDARY,AMPLITUDE=RIGID
REFDOMER,2,2,-.02355
REFD2,2,2,-.0153
SPEND,2,2,-.0153
REFP2,2,2,-.01465
*BOUNDARY
BLANK, ENCASTRE
*DIAGNOSTICS,CONTACT INITIAL OVERCLOSURE=DETAIL,DEFORMATION SPEED=DETAIL
*OUTPUT, HISTORY, TIME INTERVAL=1.E-5 
*NODE OUTPUT, NSET=NHIST 
U,V
*NODE OUTPUT, NSET=REF
U, RF
*ENERGY OUTPUT
ALLKE, ALLIE, ALLWK, ALLVD, ALLAE, ALLPD, ETOTAL
*OUTPUT, FIELD, NUMBER=4
*ELEMENT OUTPUT, ELSET=BLANK
STH, SF, S
*NODE OUTPUT
U,
*END STEP
**
** STEP 5 -- ANNEAL
**
*STEP
*ANNEAL
*END STEP
**
** STEP 6 -- SECOND PUNCHING STAGE
**
*STEP
*DYNAMIC,EXPLICIT
,.033
*BULK VISCOSITY
0.02, 0.4
*FIXED MASS SCALING, FACTOR=10, ELSET=BLANK
*DIAGNOSTICS,CONTACT INITIAL OVERCLOSURE=DETAIL,DEFORMATION SPEED=DETAIL
*CONTACT PAIR,INTERACTION=INTER
PUNCH2,BLANK_TOP
DIE2,BLANK_BOT
DOMER,BLANK_BOT
*SURFACE INTERACTION,NAME=INTER
*FRICTION
.1,
*BOUNDARY,TYPE=DISPLACEMENT,AMPLITUDE=FORM2,OP=NEW
REFP2,2,2,-0.010
*BOUNDARY,TYPE=VELOCITY,OP=NEW
REFP2,1,1,0.
REFP2,6,6,0.
REFD2,1,1,0.
REFD2,6,6,0.
SPEND,1,3,0.
REFDOMER,1,6,0.
BSYM,1,1,0.
BSYM,6,6,0.
*MONITOR,NODE=1301,DOF=2
*RESTART,WRITE,NUMBER INTERVAL=10
**
** Output for ABAQUS QA purposes
**
*FILE OUTPUT,NUMBER INTERVAL=1
*ENERGY FILE
*NODE FILE
U,RF
*EL FILE
S,
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIME MARKS=YES
*NODE OUTPUT, NSET=QA_TEST
U,
*ELEMENT OUTPUT, ELSET=QA_TEST
PEEQ,
*OUTPUT,HISTORY,TIME INTER=0.011
*ENERGY OUTPUT, VAR=ALL 
**
*END STEP

draw_spring2_std.inp

*HEADING
Deep drawing of a can bottom -- second springback
ABAQUS/Standard springback following draw_stage2.inp
SI units (kg, m, s, N)
*IMPORT, STATE=YES, UPDATE=NO
BLANK,
*IMPORT ELSET
BLANK,
*IMPORT NSET
BSYM, BLANK
*NSET, NSET=QA_TEST
BLANK,
*ELSET, ELSET=QA_TEST
BLANK,
*BOUNDARY
BSYM,1
BSYM,6
*RESTART,WRITE,FREQ=999
*STEP,NLGEOM
*STATIC
.1,1.
*BOUNDARY,FIXED
BSYM,2,2
*OUTPUT, FIELD, FREQUENCY=10 
*ELEMENT OUTPUT, ELSET=BLANK
STH, SF, S
*NODE OUTPUT
U, 
**
** Output for ABAQUS QA purposes
**
*NODE FILE,FREQ=999
U,
*EL FILE,FREQ=999
S,
*OUTPUT, FIELD,  FREQ=999
*NODE OUTPUT, NSET=QA_TEST
U,
*ELEMENT OUTPUT, ELSET=QA_TEST
S, 
*END STEP