*HEADING
BUCKLING OF X SHAPED CROSS SECTION
WITH PENALTY SELF-CONTACT
*PREPRINT,ECHO=NO,MODEL=NO,HISTORY=YES
*NODE
 601,0.,-.2,0.0
 611,0., .2,0.0
 106,0., 0.,-.2
1106,0., 0., .2
**
10601,0.05,-.2,0.0
10611,0.05, .2,0.0
10106,0.05, 0.,-.2
11106,0.05, 0., .2
*NGEN,NSET=A
601,611,1
106,1106,100
*NGEN,NSET=A1
10601,10611,1
10106,11106,100
*NODE
200601,1.,-.2,0.0
200611,1., .2,0.0
200106,1., 0.,-.2
201106,1., 0., .2
**
190601,0.95,-.2,0.0
190611,0.95, .2,0.0
190106,0.95, 0.,-.2
191106,0.95, 0., .2
*NGEN,NSET=B
200601,200611,1
200106,201106,100
*NGEN,NSET=B1
190601,190611,1
190106,191106,100
*NFILL,NSET=ALL1
A1,B1,18,10000
*ELEMENT,TYPE=S4R
  1, 601,602,10602,10601
101, 606,607,10607,10606
201, 106,206,10206,10106
301, 606,706,10706,10606
*ELGEN,ELSET=FIN1
  1, 5,1,1,   20,10000,5
*ELGEN,ELSET=FIN2
101, 5,1,1,   20,10000,5
*ELGEN,ELSET=FIN3
201, 5,100,1, 20,10000,5
*ELGEN,ELSET=FIN4
301, 5,100,1, 20,10000,5
*ELSET,ELSET=SHELLS
FIN1,FIN2,FIN3,FIN4
*SHELL SECTION, SECTION INT=GAUSS,ELSET=SHELLS,MATERIAL=STEEL
.003,
*NODE
500001,0.,-.45,-.45
500002,0., .0,-.45
500003,0., .45,-.45
500004,0.,-.45,0.
500005,0., .0,0.
500006,0., .45,0.
500007,0.,-.45,.45
500008,0., .0,.45
500009,0., .45, .45
599995,0.,0.,0.
*ELEMENT,TYPE=R3D4,ELSET=WALL1
500001, 500001,500002,500005,500004
500002, 500002,500003,500006,500005
500003, 500004,500005,500008,500007
500004, 500005,500006,500009,500008
*NODE
600001,1.,-.45,-.45
600002,1., .0,-.45
600003,1., .45,-.45
600004,1.,-.45,0.
600005,1., .0,0.
600006,1., .45,0.
600007,1.,-.45,.45
600008,1., .0,.45
600009,1., .45, .45
699995,1.,0.,0.
*ELEMENT,TYPE=R3D4,ELSET=WALL2
600001, 600001,600002,600005,600004
600002, 600002,600003,600006,600005
600003, 600004,600005,600008,600007
600004, 600005,600006,600009,600008
*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=FIN3
FIN3,
*SURFACE,TYPE=NODE,NAME=A
 A,
*SURFACE,TYPE=ELEMENT, NAME=WALL1
WALL1,SPOS
*SURFACE,TYPE=ELEMENT, NAME=FIN4
FIN4,
*SURFACE,TYPE=NODE,NAME=B
 B,
*SURFACE,TYPE=ELEMENT, NAME=WALL2
WALL2,SNEG
*SURFACE,TYPE=NODE,NAME=ALL1
 ALL1,
*SURFACE,TYPE=ELEMENT, NAME=FIN1
FIN1,
*SURFACE,TYPE=ELEMENT, NAME=FIN2
FIN2,
*RIGID BODY, ELSET=WALL1, REF NODE=599995
*RIGID BODY, ELSET=WALL2, REF NODE=699995
*NSET,NSET=NRIGID
 599995,699995
*TIE, NAME=TIE1,POSITION TOLERANCE=1.E-4
A, WALL1
*TIE, NAME=TIE2,POSITION TOLERANCE=1.E-4
B, WALL2
*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 PAIR,MECHANICAL CONSTRAINT=PENALTY
FIN1,
FIN2,
FIN3,
FIN4,
*CONTACT PAIR,MECHANICAL CONSTRAINT=PENALTY
FIN1,FIN4
FIN4,FIN2
FIN2,FIN3
FIN3,FIN1
*CONTACT PAIR,MECHANICAL CONSTRAINT=PENALTY
WALL1,ALL1
**
** 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
*END STEP