*HEADING
 Buckling of car  with spot welds
************************************************
**
*NODE,NSET=W0
*NCOPY,CHANGE=50000,OLD=W0,NEW=W50,SHIFT
*NGEN,LINE=P,NSET=WPILA
*NGEN,LINE=P,NSET=WPILB
*NGEN,LINE=P,NSET=WPILC
*NGEN,LINE=P,NSET=WPILD
*NGEN,LINE=P,NSET=WPILE
*NGEN,LINE=P,NSET=WPILF
*NGEN,LINE=P,NSET=WPILG
*NGEN,LINE=P,NSET=WPILH
*NFILL,NSET=WPIL
 WPILA,WPILB,4,1
 WPILB,WPILC,4,1
 WPILC,WPILD,4,1
 WPILD,WPILE,4,1
 WPILE,WPILF,4,1
 WPILF,WPILG,4,1
 WPILG,WPILH,4,1
************************************************
*ELEMENT,TYPE=S4R,ELSET=WPIL
*ELGEN,ELSET=WPIL
*ELGEN,ELSET=WPIL1
************************************************
** box
*NODE,NSET=B0
*NCOPY,CHANGE=50000,OLD=B0,NEW=B50,SHIFT
*NGEN,LINE=P,NSET=BPILA
*NGEN,LINE=P,NSET=BPILB
*NGEN,LINE=P,NSET=BPILC
*NGEN,LINE=P,NSET=BPILD
*NGEN,LINE=P,NSET=BPILE
*NGEN,LINE=P,NSET=BPILF
*NFILL,NSET=BPIL
 BPILA,BPILB,1,1
 BPILB,BPILC,4,1
 BPILC,BPILD,4,1
 BPILD,BPILE,4,1
 BPILE,BPILF,1,1
************************************************
*ELEMENT,TYPE=S4R,ELSET=BPIL
*ELGEN,ELSET=BPIL
*ELGEN,ELSET=BPIL1
************************************************
** roof
*NODE,NSET=ROOF
*NCOPY,CHANGE=0,OLD=ROOF,NEW=ROOF,SHIFT
*ELEMENT,ELSET=ROOF,TYPE=R3D4
*ELGEN,ELSET=ROOF
*ELEMENT,TYPE=MASS,ELSET=REF
*MASS,ELSET=REF
 1.E-3,
************************************************
**
*SHELL SECTION,ELSET=WPIL,MATERIAL=STEEL,
SECTION INTEGRATION=GAUSS
 0.002,
*SHELL SECTION,ELSET=BPIL,MATERIAL=STEEL,
SECTION INTEGRATION=GAUSS
 0.003,
************************************************
**
*MATERIAL,NAME=STEEL
*DENSITY
 7850.,
*ELASTIC
 206.8E9,0.3
*PLASTIC
170.0E6,      0.0000000E+00
180.0E6,      1.7205942E-03
190.0E6,      3.8296832E-03
200.0E6,      6.3897874E-03
210.0E6,      9.4694765E-03
220.0E6,      1.3143660E-02
230.0E6,      1.7493792E-02
240.0E6,      2.2608092E-02
************************************************
*NSET,NSET=BASE,GEN
*BOUNDARY
 BASE,ENCASTRE
************************************************
*SURFACE,TYPE=ELEMENT,NAME=BPIL
 BPIL1,SPOS
*SURFACE,TYPE=ELEMENT,NAME=WPIL
 WPIL1,SPOS
*RIGID BODY,ELSET=ROOF,REF=60000
********************
*NSET,NSET=WELDA,GEN
*NSET,NSET=WELDB,GEN
*NSET,NSET=WELDS
 WELDA,WELDB
*SURFACE INTERACTION,NAME=WELDS
*FRICTION,TAUMAX=1.4E8
 .3,
*BOND
 WELDA, 3000., 1800., 0.5E-3, 2.E-3, 0., 0.
 WELDB, 4000., 2300., 0.5E-3, 0., 0.3E-3, 0.3E-3
********************
*STEP
*DYNAMIC,EXPLICIT
 ,10.E-3
*RESTART,WRITE,NUM=2, TIMEMARKS=NO
********************
***CONTACT PAIR
** BPIL,
********************
***NSET,NSET=WELDA,GEN
** 5203,45203,10000
***NSET,NSET=WELDB,GEN
** 5211,45211,10000
***NSET,NSET=WELDS
** WELDA,WELDB
***CONTACT PAIR,INTER=WELDS,WEIGHT=0.
** BPIL, WPIL
*CONTACT
*CONTACT INCLUSIONS, ALL ELEMENT BASED
*CONTACT PAIR,INTER=WELDS,WEIGHT=0.
 BPIL, WPIL
*BOUNDARY,TYPE=VELOCITY
 60000, 1,1,   0.
 60000, 2,2, -25.
 60000, 3,3,  -2.
 60000, 4,4,  -7.
 60000, 5,5,   0.
 60000, 6,6, -78.5
********************
*NSET, NSET=QA_TEST
WELDS,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=2
*NODE OUTPUT
U,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT,HISTORY,TIME INTERVAL=2.E-5
*NODE OUTPUT,NSET=WELDS
U,
*CONTACT OUTPUT, NSET=WELDS
BONDSTAT,BONDLOAD
*NODE OUTPUT,NSET=ROOF
RF,U
*ENERGY OUTPUT
ALLIE,ALLKE,ALLWK,ALLVD,ETOTAL,ALLFD,ALLSE,
ALLAE,ALLCD,ALLPD
*FILE OUTPUT,TIMEMARKS=YES,NUM=1
*NODE FILE, NSET=QA_TEST
 U,
*ENERGY FILE
*END STEP