*HEADING
BEAM TEST PROBLEM
SECTION CONTROLS USED (HOURGLASS=ENHANCED)
*NODE
101,0.,0.
109,.4,0.
501,0.,.1
509,.4,.1
*NGEN,NSET=BOTTOM
101,109
*NGEN,NSET=TOP
501,509
*NFILL
BOTTOM,TOP,4,100
*NSET,NSET=PINNED
301,
*NSET,NSET=SYMMETRY,GENERATE
109,509,100
*ELEMENT,TYPE=CPE4R,ELSET=SOLID
101,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
10105,.2,-0.2
10106,.2,-0.2
10110,.4,-0.2
10501,0.,-0.1
10505,.2,-0.1
10506,.2,-0.1
10510,.4,-0.1
*NGEN,NSET=BOTL
10101,10105
*NGEN,NSET=BOTR
10106,10110
*NGEN,NSET=TOPL
10501,10505
*NGEN,NSET=TOPR
10506,10510
*NFILL
BOTL,TOPL,4,100
*NFILL
BOTR,TOPR,4,100
*NSET,NSET=PINNED
10301,
*NSET,NSET=SYMMETRY,GENERATE
10110,10510,100
*ELEMENT,TYPE=CPE4R,ELSET=SOLID2
10101,10101,10102,10202,10201
*ELGEN,ELSET=SOLID2
10101,4,1,1,4,100,100
*ELEMENT,TYPE=CPE4R,ELSET=SOLID2
10105,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
*NSET,NSET=PR,GEN
10106,10506,100
*MPC
PIN,PL,PR
*NODE
20101,0.,-0.4
20105,.2,-0.4
20106,.2,-0.4
20110,.4,-0.4
20501,0.,-0.3
20505,.2,-0.3
20506,.2,-0.3
20510,.4,-0.3
*NGEN,NSET=BOTLL
20101,20105
*NGEN,NSET=BOTRR
20106,20110
*NGEN,NSET=TOPLL
20501,20505
*NGEN,NSET=TOPRR
20506,20510
*NFILL
BOTLL,TOPLL,4,100
*NFILL
BOTRR,TOPRR,4,100
*NSET,NSET=PINNED
20301,
*NSET,NSET=SYMMETRY,GENERATE
20110,20510,100
*ELEMENT,TYPE=CPE4R,ELSET=SOLID3
20101,20101,20102,20202,20201
*ELGEN,ELSET=SOLID3
20101,4,1,1,4,100,100
*ELEMENT,TYPE=CPE4R,ELSET=SOLID2
20105,20106,20107,20207,20206
*ELGEN,ELSET=SOLID2
20105,4,1,1,4,100,100
*ELSET,ELSET=TOP
20401,20402,20403,20404,20405,20406,20407,20408
*ELSET,ELSET=PLL,GEN
20104,20504,100
*ELSET,ELSET=PRR,GEN
20105,20505,100
*SOLID SECTION,ELSET=SOLID,MATERIAL=MAT,CONTROL=A
.025, 
*SOLID SECTION,ELSET=SOLID2,MATERIAL=MAT,CONTROL=A
.025, 
*SOLID SECTION,ELSET=SOLID3,MATERIAL=MAT,CONTROL=A
.025, 
*SECTION CONTROLS, NAME=A, HOURGLASS=ENHANCED
*MATERIAL,NAME=MAT
*ELASTIC,TYPE=ISO
1.E9,0.,
*DENSITY
1000., 
*BOUNDARY
SYMMETRY,1,1
PINNED,1,2,0.
*RESTART,WRITE,FREQ=20
*NSET, NSET=NHIST
309, 10310, 20310
*NSET, NSET=NSUP
301,109,209,309,409,509
*SURFACE,TYPE=ELEMENT,NAME=M3R
PRR,S4
*SURFACE,TYPE=ELEMENT,NAME=M3L
PLL,S2
*NSET,NSET=QA_TEST_NHIST
NHIST,
*CONTACT PAIR,TIED,ADJUST=0.0,INTERACTION=INTER
M3L,M3R
*SURFACE INTERACTION,NAME=INTER
**
*STEP,NLGEOM=YES
*DYNAMIC,DIRECT,APPLICATION=TRANSIENT FIDELITY
2.E-4,20.E-3
*DLOAD
TOP,P3,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