** Mass scaling verification: Test of special element type with
**                            *FIXED MASS SCALING, FACTOR
*HEADING
MASS SCALING VERIFICATION - SINGLE ELEMENT DECKS WITH FACTOR OPTION
*RESTART,WRITE

****************************************
** Test 1: 2D rotary I and point mass **
****************************************
*NODE,NSET=NALL
1,     0.0,  0.0
2,     2.0,  0.0
101,   0.0,  2.0 
102,   2.0,  2.0
*NSET,NSET=DISP
2, 102,
*NSET,NSET=FIXED
1,101,
*ELEMENT,TYPE=B21,ELSET=BREF2D
1, 1,2
*ELEMENT,TYPE=MASS,ELSET=MREF2D
2, 2
*ELEMENT,TYPE=ROTARYI,ELSET=RIREF2D
3, 2
*ELEMENT,TYPE=B21,ELSET=BTEST2D
101, 101,102
*ELEMENT,TYPE=MASS,ELSET=MTEST2D
102, 102
*ELEMENT,TYPE=ROTARYI,ELSET=RITEST2D
103, 102
*ELEMENT,TYPE=ROTARYI,ELSET=RITEST2D2
104, 102
*MASS,ELSET=MREF2D
5.E-4, 
*MASS,ELSET=MTEST2D
5.E-7, 
*ROTARY INERTIA,ELSET=RIREF2D
0.0,0.0,2.5E-4
*ROTARY INERTIA,ELSET=RITEST2D
0.0,0.0,1.25E-7
*ROTARY INERTIA,ELSET=RITEST2D2
0.0,0.0,1.25E-7
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=BREF2D,MATERIAL=REFER
 1,1,0.5,-0.5
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=BTEST2D,MATERIAL=TEST
 1,1,0.5,-0.5
*ELSET,ELSET=SCALE
BTEST2D,MTEST2D,RITEST2D,RITEST2D2
*************************
** Test 2: 3D rotary I **
*************************
*NODE,NSET=NALL
1001,  4.0, 0.
1002,  6.0, 0.
1101,  4.0, 2.
1102,  6.0, 2.
*NSET,NSET=DISP
1002, 1102,
*NSET,NSET=FIXED
1001,1101,
*ELEMENT,TYPE=B31,ELSET=BREF3D
1001, 1001,1002
*ELEMENT,TYPE=MASS,ELSET=MREF3D
1002, 1002
*ELEMENT,TYPE=ROTARYI,ELSET=RIREF3D
1003, 1002
*ELEMENT,TYPE=B31,ELSET=BTEST3D
1101, 1101,1102
*ELEMENT,TYPE=MASS,ELSET=MTEST3D
1102, 1102
*ELEMENT,TYPE=ROTARYI,ELSET=RITEST3D
1103, 1102
*ORIENTATION,NAME=ORIENT3D,SYSTEM=RECTANGULAR
1.,1.,0.,  -1.,1.,0.
3, 0.
*MASS,ELSET=MREF3D
5.E-4, 
*MASS,ELSET=MTEST3D
5.E-7, 
*ROTARY INERTIA,ELSET=RIREF3D,ORIENTATION=ORIENT3D
2.5E-4, 5.0E-4, 7.5E-4
*ROTARY INERTIA,ELSET=RITEST3D,ORIENTATION=ORIENT3D
2.5E-7, 5.0E-7, 7.5E-7
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=BREF3D,MATERIAL=REFER
 1,1,0.5,-0.5
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=BTEST3D,MATERIAL=TEST
 1,1,0.5,-0.5
*ELSET,ELSET=SCALE
BTEST3D,MTEST3D,RITEST3D
******************
** TEST 3: R2D2 **
******************
*NODE,NSET=NALL
2001,  0.0, 4.
2002,  2.0, 4.
2101,  0.0, 6.
2102,  2.0, 6.
2010,  1.0, 4.5  
2110,  1.0, 6.5  
*NSET,NSET=DISP
2010, 2110,
*ELEMENT,TYPE=R2D2,ELSET=RREF3
2001, 2001,2002
*ELEMENT,TYPE=R2D2,ELSET=RTEST3
2101, 2101,2102
*RIGID BODY,ELSET=RREF3,REFNODE=2010,DENSITY=2.E-4
.5, 
*RIGID BODY,ELSET=RTEST3,REFNODE=2110,DENSITY=2.E-7
.5, 
*ELSET,ELSET=SCALE
RTEST3,
******************
** TEST 4: RAX2 **
******************
*NODE,NSET=NALL
3001,  4.0, 4.
3002,  6.0, 4.
3101,  4.0, 6.
3102,  6.0, 6.
3010,  5.0, 4.5  
3110,  5.0, 6.5  
*NSET,NSET=DISP
3010, 3110
*ELEMENT,TYPE=RAX2,ELSET=RREF4
3001, 3001,3002
*ELEMENT,TYPE=RAX2,ELSET=RTEST4
3101, 3101,3102
*RIGID BODY,ELSET=RREF4,REFNODE=3010,DENSITY=2.E-4
.5, 
*RIGID BODY,ELSET=RTEST4,REFNODE=3110,DENSITY=2.E-7
.5, 
*ELSET,ELSET=SCALE
RTEST4,
******************
** TEST 5: R3D3 **
******************
*NODE,NSET=NALL
4001,   8.0, 4.
4002,  10.0, 4.
4003,   8.0, 5.
4101,   8.0, 7.
4102,  10.0, 7.
4103,   8.0, 8.
4010,   9.0, 4.5
4110,   9.0, 7.5  
*NSET,NSET=DISP
4010,4110
*ELEMENT,TYPE=R3D3,ELSET=RREF5
4001, 4001,4002,4003
*ELEMENT,TYPE=R3D3,ELSET=RTEST5
4101, 4101,4102,4103
*RIGID BODY,ELSET=RREF5,REFNODE=4010,DENSITY=2.E-4
.25, 
*RIGID BODY,ELSET=RTEST5,REFNODE=4110,DENSITY=2.E-7
.25, 
*ELSET,ELSET=SCALE
RTEST5,
******************
** TEST 6: R3D4 **
******************
*NODE,NSET=NALL
5001,  12.0, 4.
5002,  14.0, 4.
5003,  14.0, 6.
5004,  12.0, 6.
5101,  12.0, 7.
5102,  14.0, 7.
5103,  14.0, 9.
5104,  12.0, 9.
5010,  13.0, 4.5
5110,  13.0, 7.5  
*NSET,NSET=DISP
5010,5110
*ELEMENT,TYPE=R3D4,ELSET=RREF6
5001, 5001,5002,5003,5004
*ELEMENT,TYPE=R3D4,ELSET=RTEST6
5101, 5101,5102,5103,5104
*RIGID BODY,ELSET=RREF6,REFNODE=5010,DENSITY=2.E-4
.25, 
*RIGID BODY,ELSET=RTEST6,REFNODE=5110,DENSITY=2.E-7
.25, 
*ELSET, ELSET=DUMMY
9,99,999,9999,99999
*ELSET,ELSET=SCALE
DUMMY, RTEST6,
***************************************
** REMAINDER OF DECK SAME FOR ALL TESTS
***************************************
*MATERIAL,NAME=REFER
*ELASTIC
1.E6, .25
*DENSITY
1.E-4, 
*MATERIAL,NAME=TEST
*ELASTIC
1.E6, .25
*DENSITY
1.E-7, 
*BOUNDARY
FIXED,1,6
*AMPLITUDE,NAME=RAMP
0.,0.,  .01, 1.
**************************
*NSET,NSET=QA_TEST_FIXED
FIXED,
**
*STEP
*DYNAMIC,EXPLICIT
,.00015
*FIXED MASS SCALING, FACTOR=1000.,ELSET=SCALE
*BOUNDARY,AMP=RAMP,TYPE=VELOCITY
DISP, 1,1, 10.
*RESTART,WRITE,NUMBER INTERVAL= 20
*FILE OUTPUT,NUMBER INTERVAL=4,TIMEMARKS=YES
*NODE FILE, NSET=FIXED
RF,
*ENERGY FILE
***OUTPUT,HISTORY,TIME INTERVAL=0.
***NODE OUTPUT,NSET=NALL
**RF,U
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_FIXED
RF,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=3.75E-05
*END STEP