*HEADING
      BEAM TEST PROBLEM -- 3D
 SECTION CONTROLS USED (HOURGLASS=STIFFNESS)
*NODE
101,0.,0.,0.
109,.4,0.,0.
501,0.,.1,0.
509,.4,.1,0.
1101,0.,0.,-.025
1109,.4,0.,-.025
1501,0.,.1,-.025
1509,.4,.1,-.025
*NGEN,NSET=BOTTOM
101,109
*NGEN,NSET=TOP
501,509
*NGEN,NSET=BOTTOM
1101,1109
*NGEN,NSET=TOP
1501,1509
*NFILL,NSET=ALL
BOTTOM,TOP,4,100
*NSET,NSET=PINNED
301,1301
*NSET,NSET=SYMMETRY,GENERATE
109,509,100
*NSET,NSET=SYMMETRY,GENERATE
1109,1509,100
*ELEMENT,TYPE=C3D8R,ELSET=SOLID
101, 1101,1102,1202,1201,101,102,202,201
*ELGEN,ELSET=SOLID
101, 8,1,1, 4,100,100
*ELSET,ELSET=TOP
401,402,403,404,405,406,407,408
**
*NODE
10101,0.,-0.2,0.
10105,.2,-0.2,0.
10106,.2,-0.2,0.
10110,.4,-0.2,0.
10501,0.,-0.1,0.
10505,.2,-0.1,0.
10506,.2,-0.1,0.
10510,.4,-0.1,0.
11101,0.,-0.2,-.025
11105,.2,-0.2,-.025
11106,.2,-0.2,-.025
11110,.4,-0.2,-.025
11501,0.,-0.1,-.025
11505,.2,-0.1,-.025
11506,.2,-0.1,-.025
11510,.4,-0.1,-.025
*NGEN,NSET=BOTL
10101,10105
*NGEN,NSET=BOTR
10106,10110
*NGEN,NSET=TOPL
10501,10505
*NGEN,NSET=TOPR
10506,10510
*NGEN,NSET=BOTL
11101,11105
*NGEN,NSET=BOTR
11106,11110
*NGEN,NSET=TOPL
11501,11505
*NGEN,NSET=TOPR
11506,11510
*NFILL,NSET=ALL
BOTL,TOPL,4,100
BOTR,TOPR,4,100
*NSET,NSET=PINNED
10301,11301
*NSET,NSET=SYMMETRY,GENERATE
10110,10510,100
*NSET,NSET=SYMMETRY,GENERATE
11110,11510,100
*ELEMENT,TYPE=C3D8R,ELSET=SOLID2
10101, 11101,11102,11202,11201,10101,10102,10202,10201
*ELGEN,ELSET=SOLID2
10101, 4,1,1, 4,100,100
*ELEMENT,TYPE=C3D8R,ELSET=SOLID2
10105, 11106,11107,11207,11206,10106,10107,10207,10206
*ELGEN,ELSET=SOLID2
10105, 4,1,1, 4,100,100
*ELSET,ELSET=TOP
10401,10402,10403,10404,10405,10406,10407,10408
*NSET,NSET=PL,GEN
10105,10505,100
11105,11505,100
*NSET,NSET=PR,GEN
10106,10506,100
11106,11506,100
*MPC
PIN,PL,PR
**
*NODE
20101,0.,-0.4,0.
20105,.2,-0.4,0.
20106,.2,-0.4,0.
20110,.4,-0.4,0.
20501,0.,-0.3,0.
20505,.2,-0.3,0.
20506,.2,-0.3,0.
20510,.4,-0.3,0.
21101,0.,-0.4,-.025
21105,.2,-0.4,-.025
21106,.2,-0.4,-.025
21110,.4,-0.4,-.025
21501,0.,-0.3,-.025
21505,.2,-0.3,-.025
21506,.2,-0.3,-.025
21510,.4,-0.3,-.025
*NGEN,NSET=BOTLL
20101,20105
*NGEN,NSET=BOTRR
20106,20110
*NGEN,NSET=TOPLL
20501,20505
*NGEN,NSET=TOPRR
20506,20510
*NGEN,NSET=BOTLL
21101,21105
*NGEN,NSET=BOTRR
21106,21110
*NGEN,NSET=TOPLL
21501,21505
*NGEN,NSET=TOPRR
21506,21510
*NFILL,NSET=ALL
BOTLL,TOPLL,4,100
BOTRR,TOPRR,4,100
*NSET,NSET=PINNED
20301,21301
*NSET,NSET=SYMMETRY,GENERATE
20110,20510,100
*NSET,NSET=SYMMETRY,GENERATE
21110,21510,100
*ELEMENT,TYPE=C3D8R,ELSET=SOLID3
20101, 21101,21102,21202,21201,20101,20102,20202,20201
*ELGEN,ELSET=SOLID3
20101, 4,1,1, 4,100,100
*ELEMENT,TYPE=C3D8R,ELSET=SOLID3
20105, 21106,21107,21207,21206,20106,20107,20207,20206
*ELGEN,ELSET=SOLID3
20105, 4,1,1, 4,100,100
*ELSET,ELSET=TOP
20401,20402,20403,20404,20405,20406,20407,20408
*ELSET,ELSET=PLL,GEN
20104,20404,100
*ELSET,ELSET=PRR,GEN
20105,20405,100
*SOLID SECTION,ELSET=SOLID,MATERIAL=MAT,CONTROL=A
*HOURGLASS STIFFNESS
5.E8
*SOLID SECTION,ELSET=SOLID2,MATERIAL=MAT,CONTROL=A
*HOURGLASS STIFFNESS
5.E8
*SOLID SECTION,ELSET=SOLID3,MATERIAL=MAT,CONTROL=A
*HOURGLASS STIFFNESS
5.E8
*SECTION CONTROLS,NAME=A,HOURGLASS=STIFFNESS
0.5
*MATERIAL,NAME=MAT
*ELASTIC,TYPE=ISO
1.E9,0.,
*DENSITY
1000., 
*BOUNDARY
SYMMETRY,1
PINNED,1
PINNED,2
ALL,3
*NSET, NSET=NHIST
309, 10310, 20310
*NSET, NSET=NSUP
301,1301,
109,209,309,409,509,
1109,1209,1309,1409,1509
*RESTART,WRITE,FREQ=20
*SURFACE,TYPE=ELEMENT,NAME=LE
PLL,S4
*SURFACE,TYPE=ELEMENT,NAME=LR
PRR,S6
*TIE, NAME=TIE,POSITION TOLERANCE=.00625E-3
LE,LR
*NSET,NSET=QA_TEST_NHIST
NHIST,
**
*STEP,NLGEOM=YES
*DYNAMIC,DIRECT
2.E-4,20.E-3
*DLOAD
TOP,P5,720000.
***OUTPUT,HISTORY,FREQ=1
***NODE OUTPUT,NSET=NHIST
**U2,
***NODE OUTPUT,NSET=NSUP
**RF,
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLSE,ALLAE,ALLVD,ALLCD,ALLPD,ALLFD,ALLWK,ETOTAL
*OUTPUT, FIELD, FREQ=20
*NODE OUTPUT, NSET=QA_TEST_NHIST
U, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT, FREQ=20
*END STEP