*HEADING
 Buckling of car pillar with spot welds - explicit analysis
**************************************************
*preprint, model=yes
** w-pillar
*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 pillar
*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
*NSET, NSET=QA_TEST
50113,50114,50115,50116,50117,
**************************************************
*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 INTEGRA=GAUSS
 0.002,
*SHELL SECTION,ELSET=BPIL,MATERIAL=STEEL,
 SECTION INTEGRA=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, 
*SURFACE,TYPE=ELEMENT,NAME=WPIL
 WPIL1, 
*RIGID BODY,ELSET=ROOF,REF=60000
*NSET,NSET=WELDREFNODES,GEN
 5203,45203,10000
 5211,45211,10000
*ELEMENT, TYPE=CONN3D2, ELSET=FASTCONN
1000000,5203,105203
1000001,15203,1015203
1000002,25203,1025203
1000003,35203,1035203
1000004,45203,1045203
2000000,5211,105211
2000001,15211,1015211
2000002,25211,1025211
2000003,35211,1035211
2000004,45211,1045211
*ORIENTATION, NAME=GLOBAL
1.,0.,0.,0.,1.,0.
*CONNECTOR SECTION, ELSET=FASTCONN, BEHAVIOR=BEHV,ELIMINATION=NO
CARTESIAN,CARDAN
GLOBAL,
*CONNECTOR BEHAVIOR, NAME=BEHV
*CONNECTOR DERIVED COMPONENT, NAME=normal
3,
1.0
*CONNECTOR DERIVED COMPONENT, NAME=normal
4,5
.3333,.3333
*CONNECTOR DERIVED COMPONENT, NAME=shear
6,
.5,
*CONNECTOR DERIVED COMPONENT, NAME=shear
1,2
1.0,1.0
*CONNECTOR ELASTICITY, COMPONENT=1
2.E+11,
*CONNECTOR ELASTICITY, COMPONENT=2
2.E+11,
*CONNECTOR ELASTICITY, COMPONENT=3
2.E+11,
*CONNECTOR ELASTICITY, COMPONENT=4
2.E+11,
*CONNECTOR ELASTICITY, COMPONENT=5
2.E+11,
*CONNECTOR ELASTICITY, COMPONENT=6
2.E+11,
*CONNECTOR DAMAGE INITIATION,CRITERION=FORCE
,2.E+5
*CONNECTOR POTENTIAL
normal,.02,,MACAULEY
shear,.05,,ABS
*CONNECTOR DAMAGE EVOLUTION,TYPE=MOTION,AFFECTED COMPONENTS,SOFTENING=TABULAR
1,2,3,4,5,6
0.0,0.0
1.0,1.E-7
*CONNECTOR POTENTIAL
1,
2,
3,
4,
5,
6,
*FASTENER, INTERACTION NAME=SPOTWELDS, ELSET=FASTCONN,PROPERTY=WELDPROPS
 BPIL, WPIL
*FASTENER PROPERTY, NAME=WELDPROPS
 0.001,
**************************************************
*STEP
*DYNAMIC,EXPLICIT
 ,10.E-3
********************
*CONTACT PAIR
 BPIL,
********************
*CONTACT PAIR,INTER=WELDS
 BPIL, WPIL
*SURFACE INTERACTION,NAME=WELDS
*FRICTION
 .3,
********************
*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
********************
*OUTPUT,FIELD,TIME MARKS=YES,NUMBER INTER=50
*NODE OUTPUT
 U, 
*OUTPUT,FIELD, NUMBER=2
*NODE OUTPUT, NSET=QA_TEST
 U,
*OUTPUT,HISTORY
*ELEMENT OUTPUT, ELSET=FASTCONN
CTF,CDIF,CDMG
*ENERGY OUTPUT, VAR=PRESELECT
*END STEP