*HEADING
 Buckling of car  with spot welds
************************************************
** 
*NODE,NSET=W0
 101, 0.,   0., 0.065
 105, 0.02, 0., 0.065
 109, 0.02, 0., 0.045
 113, 0.,   0., 0.045
 117, 0.,   0., 0.02
 121, 0.02, 0., 0.02
 125, 0.02, 0., 0.
 129, 0.,   0., 0.
*NCOPY,CHANGE=50000,OLD=W0,NEW=W50,SHIFT
 0.12,0.45,0.
 0.,0.45,0., 0.,0.45,1., -20.
*NGEN,LINE=P,NSET=WPILA
 101,50101,1000,, 0.03,0.2,0.065
*NGEN,LINE=P,NSET=WPILB
 105,50105,1000,, 0.05,0.2,0.065
*NGEN,LINE=P,NSET=WPILC
 109,50109,1000,, 0.05,0.2,0.045
*NGEN,LINE=P,NSET=WPILD
 113,50113,1000,, 0.03,0.2,0.045
*NGEN,LINE=P,NSET=WPILE
 117,50117,1000,, 0.03,0.2,0.02
*NGEN,LINE=P,NSET=WPILF
 121,50121,1000,, 0.05,0.2,0.02
*NGEN,LINE=P,NSET=WPILG
 125,50125,1000,, 0.05,0.2,0.0
*NGEN,LINE=P,NSET=WPILH
 129,50129,1000,, 0.03,0.2,0.0
*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
 101, 101,102,1102,1101
*ELGEN,ELSET=WPIL
 101, 28,1,1, 50,1000,1000
*ELGEN,ELSET=WPIL1
 101, 28,1,1, 48,1000,1000
************************************************
** box 
*NODE,NSET=B0
 200, 0.003,0.,0.03725
 201, 0.003,0.,0.042
 205, 0.025,0.,0.042
 209, 0.025,0.,0.023
 213, 0.003,0.,0.023
 214, 0.003,0.,0.02775
*NCOPY,CHANGE=50000,OLD=B0,NEW=B50,SHIFT
 0.12,0.45,0.
 0.,0.45,0., 0.,0.45,1., -20.
*NGEN,LINE=P,NSET=BPILA
 200,50200,1000,, 0.033,0.2,0.03725
*NGEN,LINE=P,NSET=BPILB
 201,50201,1000,, 0.033,0.2,0.042
*NGEN,LINE=P,NSET=BPILC
 205,50205,1000,, 0.055,0.2,0.042
*NGEN,LINE=P,NSET=BPILD
 209,50209,1000,, 0.055,0.2,0.023
*NGEN,LINE=P,NSET=BPILE
 213,50213,1000,, 0.033,0.2,0.023
*NGEN,LINE=P,NSET=BPILF
 214,50214,1000,, 0.033,0.2,0.02775
*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
 200, 200,201,1201,1200
*ELGEN,ELSET=BPIL
 200, 14,1,1, 50,1000,1000
*ELGEN,ELSET=BPIL1
 200, 14,1,1, 48,1000,1000
************************************************
** roof
*NODE,NSET=ROOF
 60000, 0.0125, 0., 0.0325
*NCOPY,CHANGE=0,OLD=ROOF,NEW=ROOF,SHIFT
 0.12,0.45,0.
 0.,0.45,0., 0.,0.45,1., -20.
*ELEMENT,ELSET=ROOF,TYPE=R3D4
 60101, 49101,49102,50102,50101
 60201, 49201,49202,50202,50201
*ELGEN,ELSET=ROOF
 60101, 28,1,1
 60201, 12,1,1
*ELEMENT,TYPE=MASS,ELSET=REF
 60000, 60000
*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
250.0E6,      2.8581845E-02
260.0E6,      3.5517555E-02
270.0E6,      4.3525275E-02
280.0E6,      5.2722659E-02
290.0E6,      6.3235357E-02
300.0E6,      7.5197279E-02
310.0E6,      8.8750519E-02
320.0E6,      0.1040458
330.0E6,      0.1212430
340.0E6,      0.1405106
350.0E6,      0.1620263
360.0E6,      0.1859779
370.0E6,      0.2125620
380.0E6,      0.2419857
390.0E6,      0.2744660
400.0E6,      0.3102303
410.0E6,      0.3495160
420.0E6,      0.3925720
430.0E6,      0.4396578
440.0E6,      0.4910434
450.0E6,      0.5470111
460.0E6,      0.6078544
470.0E6,      0.6738777
480.0E6,      0.7453985
490.0E6,      0.8227461
500.0E6,      0.9062610
510.0E6,      0.9962980
************************************************
*NSET,NSET=BASE,GEN
 101,129,1
 200,214,1
*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
 5203,45203,10000
*NSET,NSET=WELDB,GEN
 5211,45211,10000
*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,IMPROVED DT METHOD=NO
 ,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