*HEADING
 ROTATING CYLINDER PROBLEM
**
**  Membrane mesh.
**
*NODE, NSET=QA_TEST_ALLNODES
1001,3.,0.,0.
1010,2.,1.,0.
1101,3.,0.,.25
1110,2.,1.,.25
*NGEN,LINE=C,NSET=M3D
1001,1010,1,,2.,0.,0.
1101,1110,1,,2.,0.,.25
*ELEMENT,TYPE=M3D4R,ELSET=M3D
101,1001,1101,1102,1002
*ELGEN,ELSET=M3D
101, 9,1,1
*NODAL THICKNESS
M3D,0.1
*MEMBRANE SECTION,ELSET=M3D,MATERIAL=MAT,ORIENTATION=M3D,POISSON=MATERIAL,
NODAL THICKNESS
10., 
*EQUATION
2,
1001,1,1.,  1010,2,-1.
2,
1001,2,1.,  1010,1, 1.
2,
1101,1,1.,  1110,2,-1.
2,
1101,2,1.,  1110,1, 1.
*INITIAL CONDITIONS,TYPE=ROTATING VEL
M3D,4000., 0.,0.,0.
2.,0.,0., 2.,0.,1.
*ORIENTATION,NAME=M3D,SYSTEM=CYLINDRICAL
2.,0.,0., 2.,0.,1.
*BOUNDARY
M3D,3,3
**
**  3D brick mesh.
**
*NODE, NSET=QA_TEST_ALLNODES
2001,0.95,1.5,0.
2101,1.05,1.5,0.
2010,0.,2.45,0.
2110,0.,2.55,0.
2099,0.,1.5,0.
2201,0.95,1.5,.25
2301,1.05,1.5,.25
2210,0.,2.45,.25
2310,0.,2.55,.25
2299,0.,1.5,.25
*NGEN,LINE=C,NSET=C3D
2001,2010,1,2099
2201,2210,1,2299
2101,2110,1,2099
2301,2310,1,2299
*ELEMENT,TYPE=C3D8R,ELSET=C3D
201, 2001,2101,2102,2002,2201,2301,2302,2202
*ELGEN,ELSET=C3D
201, 9,1,1
*SOLID SECTION,ELSET=C3D,MATERIAL=MAT
*EQUATION
2,
2001,1,1., 2010,2,-1.
2,
2001,2,1., 2010,1, 1.
2,
2101,1,1., 2110,2,-1.
2,
2101,2,1., 2110,1, 1.
2,
2201,1,1., 2210,2,-1.
2,
2201,2,1., 2210,1, 1.
2,
2301,1,1., 2310,2,-1.
2,
2301,2,1., 2310,1, 1.
*MPC
SLIDER, 2001,2101,2099
SLIDER, 2010,2110,2099
SLIDER, 2201,2301,2299
SLIDER, 2210,2310,2299
*BOUNDARY
C3D,3,3
2099,1,3
2299,1,3
*INITIAL CONDITIONS,TYPE=ROTATING VEL, DEFINITION=NODES
C3D,4000., 0.,0.,0.
2099, 2299
***0.,1.5,0., 0.,1.5,1.
**
**  Shell mesh.
**
*NODE, NSET=QA_TEST_ALLNODES
3001,3.,1.5,0.
3010,2.,2.5,0.
3101,3.,1.5,.25
3110,2.,2.5,.25
*NGEN,LINE=C,NSET=S3D
3001,3010,1,,2.,1.5,0.
3101,3110,1,,2.,1.5,.25
*ELEMENT,TYPE=S4R,ELSET=S3D
301,3001,3101,3102,3002
*ELGEN,ELSET=S3D
301, 9,1,1
*NODAL THICKNESS
S3D, 0.1
*SHELL SECTION, SECTION INT=GAUSS,ELSET=S3D,MATERIAL=MAT,
NODAL THICKNESS
1., 
*EQUATION
2,
3001,1,1.,  3010,2,-1.
2,
3001,2,1.,  3010,1, 1.
2,
3001,6,1.,  3010,6,-1.
2,
3101,1,1.,  3110,2,-1.
2,
3101,2,1.,  3110,1, 1.
*INITIAL CONDITIONS,TYPE=ROTATING VEL
S3D,4000., 0.,0.,0.
2.,1.5,0., 2.,1.5,1.
*BOUNDARY
S3D,3,3
**
*MATERIAL,NAME=MAT
*ELASTIC
1.03E7,.33333
*PLASTIC
4.15E4,0.
521.15E4,10.
*DENSITY
2.508E-4,
**
**
*ELSET,ELSET=QA_TEST_ALLELEMS
M3D, C3D, S3D
**
*NSET, NSET=QA_TEST_ALLNODES
QA_TEST_ALLNODES, M3D, C3D, S3D
**
*STEP
*DYNAMIC,EXPLICIT
,.15E-3
*DLOAD
 C3D,P6,9760.
 M3D,P,-9760.
 S3D,P,-9760.
**
*FILE OUTPUT,NUMBER INTERVAL=2,TIME MARKS=YES
*EL FILE
 PEEQ,
 MISES,
 PE,
 SDV
*NODE FILE
 U,
**
*OUTPUT, FIELD, NUMBER INTERVAL=2, TIMEMARKS=YES
*ELEMENT OUTPUT, ELSET=QA_TEST_ALLELEMS
 PEEQ,
 MISES,
 PE,
 SDV
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
 U,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=7.5E-05
***OUTPUT,HISTORY,VARIABLE=PRESELECT
*END STEP