*HEADING
 Buckling of car pillar with spot welds - standard 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
**************************************************
*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
*ELSET, ELSET=WPIL2,GENERATE
 109, 47109, 1000
 110, 47110, 1000
 111, 47111, 1000
 112, 47112, 1000
 113, 47113, 1000
 114, 47114, 1000
 115, 47115, 1000
 116, 47116, 1000
 117, 47117, 1000
 118, 47118, 1000
 119, 47119, 1000
 120, 47120, 1000
*ELSET, ELSET=BPIL2,GENERATE
 200, 47200, 1000
 201, 47201, 1000
 202, 47202, 1000
 203, 47203, 1000
 204, 47204, 1000
 209, 47209, 1000
 210, 47210, 1000
 211, 47211, 1000
 212, 47212, 1000
 213, 47213, 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=BPIL2
 BPIL2, SPOS
*SURFACE,TYPE=ELEMENT,NAME=WPIL2
 WPIL2, SPOS
*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
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,
*FASTENER, INTERACTION NAME=SPOTWELDS, ELSET=FASTCONN,PROPERTY=WELDPROPS
 BPIL2, WPIL2
*FASTENER PROPERTY, NAME=WELDPROPS
 0.001,
*CONTACT PAIR,INTERACTION=INT
 BPIL2, WPIL2
*SURFACE INTERACTION,NAME=INT
*CONSTRAINT CONTROLS, NO CHECKS
**************************************************
*STEP,INC=200,NLGEOM=YES
*STATIC
 3.e-4,1.0
*BOUNDARY
 60000, 1,1,   0.
 60000, 2,2, -.25
 60000, 3,3,  -.02
 60000, 4,4,  -.07
 60000, 5,5,   0.
 60000, 6,6, -.785
********************
*OUTPUT, FIELD, FREQ=5
*NODE OUTPUT
 U, 
*ELEMENT OUTPUT
 S,PE
*ELEMENT OUTPUT, ELSET=FASTCONN, directions=yes
 CTF,CEF,CUE
*OUTPUT,HISTORY,FREQ=5
*ELEMENT OUTPUT, ELSET=FASTCONN
CTF,
*ENERGY OUTPUT, VAR=PRESELECT
*END STEP