*HEADING
 NONSTRUCTURAL MASS VERIFICATION - SINGLE ELEMENT DECKS WITH "MASS"
 UNITS.
Verifies nonstructural "total mass" distribution in proportion
to the element mass.
The elements in ELSET "noscale" have density that is 8 times that
of those in elset "scale", but both have equal volume.  A third
of the "noscale" elset mass is given as nonstructural mass 
for "element mass" proportional distribution.  This will raise the 
effective density in both elsets by factor 35/27 while still
leaving them in 8:1 ratio.  The EDT of "noscale" elements 
should be sqrt(8) times that of "scale" elements as the density
ratio has not changed.
*PREPRINT, MODEL=YES
*RESTART,WRITE
*********************
*****B21*************
*NODE,NSET=ALL-NODES
1,   0.0, 1.
2,   2.0, 1.
101, 0.0, 0.
102, 2.0, 0.
*NSET,NSET=DISPLACEMENT
2, 102
*NSET,NSET=ROTATE-Z
2, 102
*NSET,NSET=FIXED-NODES
1,101
*ELEMENT,TYPE=B21,ELSET=REFERENCE1
1, 1,2
*ELEMENT,TYPE=B21,ELSET=TEST1
101, 101,102
*ELSET,ELSET=NOSCALE-ELEMENTS
REFERENCE1,
*ELSET,ELSET=SCALE-ELEMENTS
TEST1 
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=REFERENCE1,MATERIAL=NOSCALE-MATERIAL
 1,1,0.5,-0.5
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=TEST1,MATERIAL=SCALE-MATERIAL
 1,1,0.5,-0.5
*********************
*****B22*************
*NODE,NSET=ALL-NODES
401,   0.0, 1.
402,   2.0, 1.
403,   1.0, 1.
501, 0.0, 0.
502, 2.0, 0.
503, 1.0, 0.
*NSET,NSET=DISPLACEMENT
402, 502
*NSET,NSET=DISPLACEMENT_mid
403, 503
*NSET,NSET=ROTATE-Z
402, 502
*NSET,NSET=ROTATE-Z_mid
403, 503
*NSET,NSET=FIXED-NODES
401,501
*ELEMENT,TYPE=B22,ELSET=REFERENCE1a
401, 401,403,402
*ELEMENT,TYPE=B22,ELSET=TEST1a
501, 501,503,502
*ELSET,ELSET=NOSCALE-ELEMENTS
REFERENCE1a,
*ELSET,ELSET=SCALE-ELEMENTS
TEST1a,
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=REFERENCE1a,MATERIAL=NOSCALE-MATERIAL
 1,1,0.5,-0.5
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=TEST1a,MATERIAL=SCALE-MATERIAL
 1,1,0.5,-0.5
*********************
****B31*************
*NODE,NSET=ALL-NODES
1001,   0.0, 5.
1002,   2.0, 5., 2.
1101, 0.0, 4.
1102, 2.0, 4., 2.
*NSET,NSET=DISPLACEMENT
1002, 1102
*NSET,NSET=ROTATE-Y
1002, 1102
*NSET,NSET=ROTATE-Z
1002, 1102
*NSET,NSET=FIXED-NODES
1001,1101
*ELEMENT,TYPE=B31,ELSET=REFERENCE2
1001, 1001,1002
*ELEMENT,TYPE=B31,ELSET=TEST2
1101, 1101,1102
*ELSET,ELSET=NOSCALE-ELEMENTS
REFERENCE2,
*ELSET,ELSET=SCALE-ELEMENTS
TEST2,
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=REFERENCE2,MATERIAL=NOSCALE-MATERIAL
 1,1,0.5,-0.5
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=TEST2,MATERIAL=SCALE-MATERIAL
 1,1,0.5,-0.5
*************************************
****B32*************
*NODE,NSET=ALL-NODES
1401,   0.0, 5.
1402,   2.0, 5., 2.
1403,   1.0, 5., 1.
1501, 0.0, 4.
1502, 2.0, 4., 2.
1503, 1.0, 4., 1.
*NSET,NSET=DISPLACEMENT
1402, 1502
*NSET,NSET=DISPLACEMENT_mid
1403, 1503
*NSET,NSET=ROTATE-Y
1402, 1502
*NSET,NSET=ROTATE-Z
1402, 1502
*NSET,NSET=ROTATE-Y_mid
1403, 1503
*NSET,NSET=ROTATE-Z_mid
1403, 1503
*NSET,NSET=FIXED-NODES
1401,1501
*ELEMENT,TYPE=B32,ELSET=REFERENCE2a
1401, 1401,1403,1402
*ELEMENT,TYPE=B32,ELSET=TEST2a
1501, 1501,1503,1502
*ELSET,ELSET=NOSCALE-ELEMENTS
REFERENCE2a,
*ELSET,ELSET=SCALE-ELEMENTS
TEST2a,
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=REFERENCE2a,MATERIAL=NOSCALE-MATERIAL
 1,1,0.5,-0.5
*BEAM SECTION,SECTION=TRAPEZOID,ELSET=TEST2a,MATERIAL=SCALE-MATERIAL
 1,1,0.5,-0.5
*************************************
*****S3R *******
*NODE,NSET=ALL-NODES
2001,   5.0, 5.
2002,   7.0, 5.
2003,   7.0, 8.
2004,   5.0, 8.
2101, 5.0, 0.
2102, 7.0, 0.
2103, 7.0, 3.
2104, 5.0, 3.
*NSET,NSET=DISPLACEMENT
2002, 2102
*NSET,NSET=ROTATE-Y
2002, 2102
*NSET,NSET=FIXED-NODES
2001,2101,2004,2104
*ELEMENT,TYPE=S3R,ELSET=REFERENCE3
2001, 2001,2002,2003
2002, 2001,2003,2004
*ELEMENT,TYPE=S3R,ELSET=TEST3
2101, 2101,2102,2103
2102, 2101,2103,2104
*ELSET,ELSET=NOSCALE-ELEMENTS
REFERENCE3,
*ELSET,ELSET=SCALE-ELEMENTS
TEST3,
***NODAL THICKNESS
**2001, .2
**2002, .2
**2003, .4
**2004, .4
**2101, .2
**2102, .2
**2103, .4
**2104, .4
*SHELL SECTION,
ELSET=REFERENCE3,MATERIAL=NOSCALE-MATERIAL
.5, 
*SHELL SECTION,
ELSET=TEST3,MATERIAL=SCALE-MATERIAL
.5, 
*****************************
*****S4R WITH REBAR *********
*NODE,NSET=ALL-NODES
3001,   5.0, 15.
3002,   7.0, 15.
3003,   7.0, 18.
3004,   5.0, 18.
3101, 5.0, 10.
3102, 7.0, 10.
3103, 7.0, 13.
3104, 5.0, 13.
*NSET,NSET=DISPLACEMENT
3002, 3102
*NSET,NSET=ROTATE-Y
3002, 3102
*NSET,NSET=FIXED-NODES
3001,3101,3004,3104
*ELEMENT,TYPE=S4R,ELSET=REFERENCE4
3001, 3001,3002,3003,3004
*ELEMENT,TYPE=S4R,ELSET=TEST4
3101, 3101,3102,3103,3104
*ELSET,ELSET=NOSCALE-ELEMENTS
REFERENCE4,
*ELSET,ELSET=SCALE-ELEMENTS
TEST4,
*SHELL SECTION,ELSET=REFERENCE4,MATERIAL=NOSCALE-MATERIAL
.5, 
*REBAR LAYER
REB4,.5, .5,0.2,NOSCALE-MATERIAL,0.0,1
*SHELL SECTION,ELSET=TEST4,MATERIAL=SCALE-MATERIAL
.5, 
*REBAR LAYER
REB4,.5, .5,0.2,SCALE-MATERIAL,0.0,1
*************************
*****SAX1 WITH REBAR ****
*NODE,NSET=ALL-NODES
4001,   10.0, 5.
4002,   12.0, 5.
4101, 10.0, 0.
4102, 12.0, 0.
*NSET,NSET=DISPLACEMENT
4002, 4102
*NSET,NSET=ROTATE-Z
4002, 4102
*NSET,NSET=FIXED-NODES
4001,4101
*ELEMENT,TYPE=SAX1,ELSET=REFERENCE5
4001, 4001,4002
*ELEMENT,TYPE=SAX1,ELSET=TEST5
4101, 4101,4102
*ELSET,ELSET=NOSCALE-ELEMENTS
REFERENCE5,
*ELSET,ELSET=SCALE-ELEMENTS
TEST5,
*SHELL SECTION,ELSET=REFERENCE5,MATERIAL=NOSCALE-MATERIAL
.5,
*REBAR LAYER, GEOMETRY=ANGULAR
REB5,.25,3.581,0.2,NOSCALE-MATERIAL,0.0,1 
*SHELL SECTION,ELSET=TEST5,MATERIAL=SCALE-MATERIAL
.5,
*REBAR LAYER, GEOMETRY=ANGULAR
REB5,.25,3.581,0.2,SCALE-MATERIAL,0.0,1  
*********************************
*****SAX1 **
*NODE,NSET=ALL-NODES
5001,   10.0, 15.
5002,   12.0, 15.
5101, 10.0, 10.
5102, 12.0, 10.
*NSET,NSET=DISPLACEMENT
5002, 5102
*NSET,NSET=ROTATE-Z
5002, 5102
*NSET,NSET=FIXED-NODES
5001,5101
*ELEMENT,TYPE=SAX1,ELSET=REFERENCE6
5001, 5001,5002
*ELEMENT,TYPE=SAX1,ELSET=TEST6
5101, 5101,5102
*ELSET,ELSET=NOSCALE-ELEMENTS
REFERENCE6,
*ELSET,ELSET=SCALE-ELEMENTS
TEST6,
***NODAL THICKNESS
**5001, .2
**5002, .4
**5101, .2
**5102, .4
*SHELL SECTION,
ELSET=REFERENCE6,MATERIAL=NOSCALE-MATERIAL
.5, 
*SHELL SECTION,
ELSET=TEST6,MATERIAL=SCALE-MATERIAL
.5, 
******************************
** REMAINDER OF DECK SAME FOR ALL TESTS
******************************
*MATERIAL,NAME=NOSCALE-MATERIAL
*ELASTIC
1.E6, .25
*DENSITY
1.E-4, 
*MATERIAL,NAME=SCALE-MATERIAL
*ELASTIC
1.E6, .25
*DENSITY
0.125E-4, 
*BOUNDARY
FIXED-NODES,1,6
**************************
*AMPLITUDE,NAME=SMOOTH,DEFINITION=SMOOTH STEP
0.,0.,  .00015, 1.
**************************
*NSET,NSET=QA_TEST
FIXED-NODES,
*ELSET,ELSET=QA_TEST
SCALE-ELEMENTS,
NOSCALE-ELEMENTS,
*NONSTRUCTURAL MASS, UNITS=TOTAL MASS, ELSET=QA_TEST
0.6923333e-2
*STEP
*DYNAMIC,EXPLICIT
,.00015
*BOUNDARY,AMP=SMOOTH
DISPLACEMENT, 2,2, .05
DISPLACEMENT_mid, 2,2, .025
ROTATE-Y, 5,5, .05
ROTATE-Y_mid, 5,5, .025
ROTATE-Z, 6,6, .05
ROTATE-Z_mid, 6,6, .025
**********************
*FILE OUTPUT,NUMBER INTERVAL=4,TIMEMARKS=YES
*NODE FILE, NSET=FIXED-NODES
RF,
** RF not expected to be noisy.
*EL FILE
EMSF, EDT
** EMSF, EDT not expected to be noisy.
***********************
*OUTPUT, FIELD,NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST
EMSF, EDT
*NODE OUTPUT,NSET=QA_TEST
RF,
*****
***OUTPUT,HISTORY,TIME INTERVAL=0.
***NODE OUTPUT,NSET=FIXED-NODES
**RF,
***ELEMENT OUTPUT, ELSET=SCALE-ELEMENTS
**EDT,
***ELEMENT OUTPUT, ELSET=NOSCALE-ELEMENTS
**EDT,
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=3.75E-05
*PRINT,DMASS=YES
*END STEP