*HEADING
BUCKLING OF X SHAPED CROSS SECTION
WITH GENERAL CONTACT
*PREPRINT,MODEL=YES,HISTORY=YES,CONTACT=YES
*NODE
**
*NGEN,NSET=A
*NGEN,NSET=A1
*NODE
**
*NGEN,NSET=B
*NGEN,NSET=B1
*NFILL,NSET=ALL1
A1,B1,18,10000
*ELEMENT,TYPE=S4R
*ELGEN,ELSET=FIN1
*ELGEN,ELSET=FIN2
*ELGEN,ELSET=FIN3
*ELGEN,ELSET=FIN4
*ELSET,ELSET=SHELLS
FIN1,FIN2,FIN3,FIN4
*SHELL SECTION, SECTION INT=GAUSS,ELSET=SHELLS,MATERIAL=STEEL
.003,
*NODE
*ELEMENT,TYPE=R3D4,ELSET=WALL1
*NODE
*ELEMENT,TYPE=R3D4,ELSET=WALL2
*MATERIAL,NAME=STEEL
*DENSITY
7850.,
*ELASTIC
200.E9,.3
*PLASTIC
250.E6,0.
700.E6,1.
*BOUNDARY
599995,2,3
599995,4,4
599995,6,6
699995,1,6
*INITIAL CONDITIONS,TYPE=VELOCITY
 599995, 1,1, 50.
 599995, 5,5, 78.54
*RESTART,TIMEMARKS=NO,WRITE,NUM=10
**
*SURFACE,TYPE=ELEMENT, NAME=WALL1
WALL1,SPOS
*SURFACE,TYPE=ELEMENT, NAME=WALL2
WALL2,SNEG
*SURFACE,TYPE=NODE,NAME=A
 A,
*SURFACE,TYPE=NODE,NAME=B
 B,
**
*RIGID BODY, ELSET=WALL1, REF NODE=599995
*RIGID BODY, ELSET=WALL2, REF NODE=699995
*NSET,NSET=NRIGID
*TIE, NAME=TIE1,POSITION TOLERANCE=1.E-4
A, WALL1
*TIE, NAME=TIE2,POSITION TOLERANCE=1.E-4
B, WALL2
*IMPERFECTION, STEP=1, FILE=sscxsec_bkl
*STEP
*DYNAMIC,EXPLICIT
,0.7E-2
*BOUNDARY,TYPE=VEL
** compress 35 percent axially and rotate 31 degrees during the step
599995, 1,1, 50.
599995, 5,5, 78.54
**
*CONTACT
*CONTACT INCLUSIONS,ALL ELEMENT BASED
**
** Output requests for postprocessing
**
*OUTPUT,FIELD,VAR=PRESELECT
*OUTPUT,HISTORY,TIME INTERVAL=2.E-5
*NODE OUTPUT, NSET=NRIGID
U, RF
*ENERGY OUTPUT
ALLIE,ALLKE,ALLWK
**
** Output requests for qa testing
**
*ELSET,ELSET=QA_TESTEL
FIN4,
*NSET,NSET=QA_TESTN,ELSET=FIN4
*NSET,NSET=QA_TESTRB
NRIGID,
*OUTPUT,FIELD,NUMBER INTERVAL=1
*ELEMENT OUTPUT,ELSET=QA_TESTEL
PEEQ,
*NODE OUTPUT,NSET=QA_TESTN
U,
*NODE OUTPUT, NSET=QA_TESTRB
U, RF
*OUTPUT,HISTORY,TIME INTERVAL=1.75E-3
*ENERGY OUTPUT,VAR=PRESELECT
*FILE OUTPUT,NUMBER=4,TIMEMARKS=YES
*ENERGY FILE
*NODE FILE
U,RF
*END STEP