** 
** This is a test of the *BEAM ADDED INERTIA for 3d quadratic  beams
*PREPRINT,MODEL=YES,HISTORY=YES
*RESTART,TIMEMARKS=NO,WRITE,NUM=8
*PART, NAME=BEAM1
*NODE,NSET=BEAMS
1, 0.0, 0.0, 0.
2, 5.0, 0.0, 0.
3, 2.5, 0.0, 0.
*ELEMENT,TYPE=B32,ELSET=GEN1
1,1,3,2
*BEAM GENERAL SECTION, SECTION=GENERAL, ELSET=GEN1, DENSITY=1000.0, POISSON=0.3
 0.25,1.5, 0.0, 1.5, 0.010416667
 0.0, 0.0, -1.0
3.E9, 1.1538E9, .001
*BEAM ADDED INERTIA
1700.,.25,.25,0.,10200.,10200.,0.
*END PART
*PART, NAME=BEAM11
*NODE,NSET=BEAMS
101, 0.0, 0.0, 10.
102, 5.0, 0.0, 10.
103, 2.5, 0.0, 10.
10100,0.0,.25,9.75
10200,5.0,.25,9.75
10300,2.5,.25,9.75
*ELEMENT,TYPE=B32,ELSET=GEN2
101,101,103,102
*BEAM GENERAL SECTION, SECTION=GENERAL, ELSET=GEN2, DENSITY=1000.0, POISSON=0.3
 0.25,1.5, 0.0, 1.5, 0.010416667
 0.0, 0.0, -1.0
3.E9, 1.1538E9, .001
*ELEMENT, TYPE=MASS, ELSET=MASS2
1002,10100
1003,10200
*ELEMENT, TYPE=MASS, ELSET=MASS21
1004,10300
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY2
2002,10100
2003,10200
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY21
2004,10300
*MASS, ELSET=MASS2
1416.667
*MASS, ELSET=MASS21
5666.6667
*ROTARY INERTIA, ELSET=ROTARY2
0.,8500,8500.
*ROTARY INERTIA, ELSET=ROTARY21
0.,34000.,34000.
*END PART
*PART, NAME=BEAM2
*NODE,NSET=BEAMS
51, 0.0, 0.5, 0.
52, 5.0, 0.5, 0.
53, 2.5, 0.5, 0.
*ELEMENT,TYPE=B32,ELSET=ARBITRAR1
51,51,53,52
*BEAM SECTION,SECTION=ARBITRARY,ELSET=ARBITRAR1,MATERIAL=TEST1
4,.5,.5,-.5,.5,.1
-.5,-.5,.1
.5,-.5,.1
.5,.5,.1
*BEAM ADDED INERTIA
380.,.25,.25,0.,60.,60.,0.
*END PART
*PART, NAME=BEAM22
*NODE,NSET=BEAMS
151, 0.0, 0.5,10.
152, 5.0, 0.5,10.
153, 2.5, 0.5,10.
15100,0.0,.75,9.75
15200,5.0,.75,9.75
15300,2.5,.75,9.75
*ELEMENT,TYPE=B32,ELSET=ARBITRAR2
151,151,153,152
*BEAM SECTION,SECTION=ARBITRARY,ELSET=ARBITRAR2,MATERIAL=TEST1
4,.5,.5,-.5,.5,.1
-.5,-.5,.1
.5,-.5,.1
.5,.5,.1
*ELEMENT, TYPE=MASS, ELSET=MASS1
1000,15100
1001,15200
*ELEMENT, TYPE=MASS, ELSET=MASS11
1002,15300
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY1
2000,15100
2001,15200
*ELEMENT, TYPE=ROTARYI, ELSET=ROTARY11
2002,15300
*MASS, ELSET=MASS1
316.667
*MASS, ELSET=MASS11
1266,667
*ROTARY INERTIA, ELSET=ROTARY1
0.,50.,50.
*ROTARY INERTIA, ELSET=ROTARY11
0.,200.,200.
*END PART
*ASSEMBLY, NAME=ASSEMBLY
*INSTANCE, NAME=BEAM1_A, PART=BEAM1
 0.0, 0.0, 0.0
*END INSTANCE
*INSTANCE, NAME=BEAM11_A, PART=BEAM11
 0.0, 0.0, 0.0
*END INSTANCE
*INSTANCE, NAME=BEAM2_A, PART=BEAM2
 0.0, 0.0, 0.0
*END INSTANCE
*INSTANCE, NAME=BEAM22_A, PART=BEAM22
 0.0, 0.0, 0.0
*END INSTANCE
*MPC
BEAM,BEAM11_A.101,BEAM11_A.10100
BEAM,BEAM11_A.102,BEAM11_A.10200
BEAM,BEAM11_A.103,BEAM11_A.10300
BEAM,BEAM22_A.151,BEAM22_A.15100
BEAM,BEAM22_A.152,BEAM22_A.15200
BEAM,BEAM22_A.153,BEAM22_A.15300
*END ASSEMBLY
*MATERIAL,NAME=TEST1
*ELASTIC
1.E9,0.
*DENSITY
50., 
*BOUNDARY
ASSEMBLY.BEAM1_A.1,1,6
ASSEMBLY.BEAM11_A.101,1,6
ASSEMBLY.BEAM2_A.51,1,6
ASSEMBLY.BEAM22_A.151,1,6
*NSET, NSET=NHIST
ASSEMBLY.BEAM1_A.2,1,6
ASSEMBLY.BEAM11_A.102,1,6
ASSEMBLY.BEAM2_A.52,1,6
ASSEMBLY.BEAM22_A.152,1,6
*NSET,NSET=QA_TEST_ALLNODES
ASSEMBLY.BEAM1_A.BEAMS,ASSEMBLY.BEAM11_A.BEAMS,
ASSEMBLY.BEAM2_A.BEAMS,ASSEMBLY.BEAM22_A.BEAMS
*AMPLITUDE, NAME=AMP, DEFINITION=SMOOTH STEP
0.0,0.0,.125,1.0
**
*STEP
*DYNAMIC,EXPLICIT
 ,.125
*CLOAD, AMPLITUDE=AMP
NHIST,2,-5.E6
*FILE OUTPUT,NUMBER INTERVAL=4
*NODE FILE, NSET=NHIST
 U,
*OUTPUT, FIELD, NUMBER INTERVAL=4,TIME MARKS=YES
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
 U,
*END STEP