 *HEADING
 Buckling of car pillar with spot welds - explicit analysis
**************************************************
*preprint, model=yes
** w-pillar
*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 pillar
*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
*NSET, NSET=QA_TEST
**************************************************
*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 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
**************************************************
*NSET,NSET=BASE,GEN
*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
*ELEMENT, TYPE=CONN3D2, ELSET=FASTCONN
*ORIENTATION, NAME=GLOBAL
*CONNECTOR SECTION, ELSET=FASTCONN, BEHAVIOR=BEHV,ELIMINATION=NO
CARTESIAN,CARDAN
GLOBAL,
*CONNECTOR BEHAVIOR, NAME=BEHV
*CONNECTOR DERIVED COMPONENT, NAME=normal
*CONNECTOR DERIVED COMPONENT, NAME=normal
*CONNECTOR DERIVED COMPONENT, NAME=shear
*CONNECTOR DERIVED COMPONENT, NAME=shear
*CONNECTOR ELASTICITY, COMPONENT=1
*CONNECTOR ELASTICITY, COMPONENT=2
*CONNECTOR ELASTICITY, COMPONENT=3
*CONNECTOR ELASTICITY, COMPONENT=4
*CONNECTOR ELASTICITY, COMPONENT=5
*CONNECTOR ELASTICITY, COMPONENT=6
*CONNECTOR DAMAGE INITIATION,CRITERION=FORCE
*CONNECTOR POTENTIAL
normal,.02,,MACAULEY
shear,.05,,ABS
*CONNECTOR DAMAGE EVOLUTION,TYPE=MOTION,AFFECTED COMPONENTS,SOFTENING=TABULAR
*CONNECTOR POTENTIAL
*FASTENER, INTERACTION NAME=SPOTWELDS, ELSET=FASTCONN,PROPERTY=WELDPROPS,COUPLING=STRUCTURAL
 BPIL, WPIL
*FASTENER PROPERTY, NAME=WELDPROPS
**************************************************
*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