*HEADING
VERIFICATION OF ELEMENTS LOADED WITH GRAVITY LOAD.
ANALYTICAL RESULTS:
	T2D2 	- 31400
	T3D2 	- 31400
	B21  	- 6280
	B31 	- 6280
 	SAX1 	- 19729.202
	MASS	- 7850
	C3D8R 	- 62800
	C3D4 	- 62800
	C3D6	- 62800
	CPE4R	- 62800
	CPE3	- 62800
	CPS4R	- 62800
	CPS3	- 62800
	S4R	- 62800
	S3R	- 62800
	M3D4R	- 62800
	M3D3	- 62800
	CAX4R	- 54255.305
The nonstructural mass is added to the elements to make
sure that the GRAV load is applied on the element mass 
that includes the nonstructural mass contribution. The 
under lying primitive is reduced by reducing the material
density while maintaining the same element total mass.
The RF results should not change. - prasad
*PREPRINT, MODEL=YES
*NODE,NSET=T2D2
10001,2.0,0.0
10002,2.0,2.0
*NODE,NSET=T3D2
10003,3.0,0.0
10004,3.0,2.0
*NODE,NSET=B21
10005,4.0,0.0
10006,4.0,2.0
*NODE,NSET=B31
10007,5.0,0.0
10008,5.0,2.0
*NODE,NSET=SAX1
10009,1.0,0.0
10010,1.0,2.0
*NODE,NSET=C3D8R
10011,7.0,0.0,2.0
10012,9.0,0.0,2.0
10013,9.0,0.0,0.0
10014,7.0,0.0,0.0
10015,7.0,2.0,2.0
10016,9.0,2.0,2.0
10017,9.0,2.0,0.0
10018,7.0,2.0,0.0
*NODE,NSET=CPE4R
10101,10.0,0.0
10102,12.0,0.0
10103,12.0,2.0
10104,10.0,2.0
*NODE,NSET=CPS4R
10201,13.0,0.0
10202,15.0,0.0
10203,15.0,2.0
10204,13.0,2.0
*NODE,NSET=CAX4R
10301,22.0,0.0
10302,24.0,0.0
10303,24.0,2.0
10304,22.0,2.0
*NODE,NSET=S4R
10401,16.0,0.0
10402,18.0,0.0
10403,18.0,2.0
10404,16.0,2.0
*NODE,NSET=M3D4R
10501,19.0,0.0
10502,21.0,0.0
10503,21.0,2.0
10504,19.0,2.0
*NODE,NSET=MASS
10601,6.0,1.0
*NODE,NSET=CPE3
20101,10.0,3.0
20102,12.0,3.0
20103,12.0,5.0
20104,10.0,5.0
*NODE,NSET=CPS3
20201,13.0,3.0
20202,15.0,3.0
20203,15.0,5.0
20204,13.0,5.0
*NODE,NSET=C3D4
20011,7.0,3.0,2.0
20012,9.0,3.0,2.0
20013,9.0,3.0,0.0
20014,7.0,3.0,0.0
20015,7.0,5.0,2.0
20016,9.0,5.0,2.0
20017,9.0,5.0,0.0
20018,7.0,5.0,0.0
*NODE,NSET=C3D6
30011,7.0,6.0,2.0
30012,9.0,6.0,2.0
30013,9.0,6.0,0.0
30014,7.0,6.0,0.0
30015,7.0,8.0,2.0
30016,9.0,8.0,2.0
30017,9.0,8.0,0.0
30018,7.0,8.0,0.0
*NODE,NSET=CAX3
20301,22.0,3.0
20302,24.0,3.0
20303,24.0,5.0
20304,22.0,5.0
*NODE,NSET=S3R
20401,16.0,3.0
20402,18.0,3.0
20403,18.0,5.0
20404,16.0,5.0
*NODE,NSET=M3D3
20501,19.0,3.0
20502,21.0,3.0
20503,21.0,5.0
20504,19.0,5.0
*NODE,NSET=CPS6M
40001,25.0,0.0
40002,26.0,0.0
40003,27.0,0.0
40004,25.0,1.0
40005,26.0,1.0
40006,27.0,1.0
40007,25.0,2.0
40008,26.0,2.0
40009,27.0,2.0
*NODE,NSET=CPE6M
40011,25.0,3.0
40012,26.0,3.0
40013,27.0,3.0
40014,25.0,4.0
40015,26.0,4.0
40016,27.0,4.0
40017,25.0,5.0
40018,26.0,5.0
40019,27.0,5.0
*NODE,NSET=CAX6M
40021,25.0,6.0
40022,26.0,6.0
40023,27.0,6.0
40024,25.0,7.0
40025,26.0,7.0
40026,27.0,7.0
40027,25.0,8.0
40028,26.0,8.0
40029,27.0,8.0
*NODE,NSET=C3D10M
40031,28.0,0.0,0.0
40032,30.0,0.0,0.0
40033,28.0,2.0,0.0
40034,28.0,0.0,2.0
40035,29.0,0.0,0.0
40036,29.0,1.0,0.0
40037,28.0,1.0,0.0
40038,28.0,0.0,1.0
40039,29.0,0.0,1.0
40040,28.0,1.0,1.0
**
*ELEMENT,TYPE=T2D2,ELSET=T2D2
10001,10001,10002
*ELEMENT,TYPE=T3D2,ELSET=T3D2
10002,10003,10004
*ELEMENT,TYPE=B21,ELSET=B21
10003,10005,10006
*ELEMENT,TYPE=B31,ELSET=B31
10004,10007,10008
*ELEMENT,TYPE=SAX1,ELSET=SAX1
10005,10009,10010
*ELEMENT,TYPE=C3D8R,ELSET=C3D8R
10011,10011,10012,10013,10014,10015,10016,10017,10018
*ELEMENT,TYPE=CPE4R,ELSET=CPE4R
10101,10101,10102,10103,10104
*ELEMENT,TYPE=CPS4R,ELSET=CPS4R
10201,10201,10202,10203,10204
*ELEMENT,TYPE=CAX4R,ELSET=CAX4R
10301,10301,10302,10303,10304
*ELEMENT,TYPE=S4R,ELSET=S4R
10401,10401,10402,10403,10404
*ELEMENT,TYPE=M3D4R,ELSET=M3D4R
10501,10501,10502,10503,10504
*ELEMENT,TYPE=MASS,ELSET=MASS
10601,10601
*ELEMENT,TYPE=CPE3,ELSET=CPE3
20101,20101,20102,20104
20102,20104,20102,20103
*ELEMENT,TYPE=CPS3,ELSET=CPS3
20201,20201,20202,20204
20202,20204,20202,20203
*ELEMENT,TYPE=C3D4,ELSET=C3D4
20011,20014,20013,20018,20011
20012,20011,20013,20016,20012
20013,20013,20017,20018,20016
20014,20011,20016,20018,20015
20015,20011,20016,20013,20018
*ELEMENT,TYPE=C3D6,ELSET=C3D6
30011,30014,30013,30018,30011,30012,30015
30012,30013,30017,30018,30012,30016,30015
*ELEMENT,TYPE=CAX3,ELSET=CAX3
20301,20301,20302,20304
20302,20304,20302,20303
*ELEMENT,TYPE=M3D3,ELSET=M3D3
20501,20501,20502,20504
20502,20504,20502,20503
*ELEMENT,TYPE=S3R,ELSET=S3R
20401,20401,20402,20404
20402,20404,20402,20403
*ELEMENT,TYPE=CPS6M,ELSET=CPS6M
40001,40001,40003,40009,40002,40006,40005
40002,40001,40009,40007,40005,40008,40004
*ELEMENT,TYPE=CPS6M,ELSET=CPE6M
40011,40011,40013,40019,40012,40016,40015
40012,40011,40019,40017,40015,40018,40014
*ELEMENT,TYPE=CAX6M,ELSET=CAX6M
40021,40021,40023,40029,40022,40026,40025
40022,40021,40029,40027,40025,40028,40024
*ELEMENT,TYPE=C3D10M,ELSET=C3D10M
40031,40031,40032,40033,40034,40035,40036,40037,
40038,40039,40040,
**
*SOLID SECTION,ELSET=T2D2,MATERIAL=A1ELAS_STRUCT
0.2,
*SOLID SECTION,ELSET=T3D2,MATERIAL=A1ELAS_STRUCT
0.2,
*BEAM SECTION,SECTION=RECT,ELSET=B21,MATERIAL=A1ELAS_STRUCT
0.2,0.2
*BEAM SECTION,SECTION=RECT,ELSET=B31,MATERIAL=A1ELAS_STRUCT
0.2,0.2
*SHELL SECTION,ELSET=S4R,MATERIAL=A1ELAS_STRUCT,POISSON=0.3
0.2,3
*SHELL SECTION,ELSET=S3R,MATERIAL=A1ELAS_STRUCT,POISSON=0.3
0.2,3
*SHELL SECTION,ELSET=SAX1,MATERIAL=A1ELAS_CONT,POISSON=0.3
0.02,3
*SOLID SECTION,ELSET=C3D8R,MATERIAL=A1ELAS_CONT
*SOLID SECTION,ELSET=C3D4,MATERIAL=A1ELAS_CONT
*SOLID SECTION,ELSET=C3D6,MATERIAL=A1ELAS_CONT
*SOLID SECTION,ELSET=CPE4R,MATERIAL=A1ELAS_CONT
2.0,
*SOLID SECTION,ELSET=CPE3,MATERIAL=A1ELAS_CONT
2.0,
*SOLID SECTION,ELSET=CPS3,MATERIAL=A1ELAS_CONT
2.0,
*SOLID SECTION,ELSET=CPS4R,MATERIAL=A1ELAS_CONT
2.0,
*SOLID SECTION,ELSET=CAX4R,MATERIAL=A1ELAS_CONT
*SOLID SECTION,ELSET=CAX3,MATERIAL=A1ELAS_CONT
*MEMBRANE SECTION,ELSET=M3D4R,MATERIAL=A1ELAS_STRUCT,POISSON=0.3
0.2,
*MEMBRANE SECTION,ELSET=M3D3,MATERIAL=A1ELAS_STRUCT,POISSON=0.3
0.2,
*SOLID SECTION,ELSET=CPS6M,MATERIAL=A1ELAS_CONT
2.0
*SOLID SECTION,ELSET=CPE6M,MATERIAL=A1ELAS_CONT
2.0
*SOLID SECTION,ELSET=CAX6M,MATERIAL=A1ELAS_CONT
*SOLID SECTION,ELSET=C3D10M,MATERIAL=A1ELAS_CONT
**
*MATERIAL,NAME=A1ELAS
*ELASTIC
193.1E9,0.3
*DENSITY
7850,
*MATERIAL,NAME=A1ELAS_CONT
*ELASTIC
193.1E9,0.3
*DENSITY
3925.,
*MATERIAL,NAME=A1ELAS_STRUCT
*ELASTIC
193.1E9,0.3
*DENSITY
1962.5,
**
*MASS,ELSET=MASS
785.0,
**
*ELSET,ELSET=EALL
T2D2,T3D2,
B21,B31,
SAX1,
C3D8R,
CPE4R,CPS4R,CAX4R,
S4R,M3D4R,
MASS,
CPE3,CPS3,C3D4,C3D6,CAX3,M3D3,S3R,
CPS6M,CPE6M,CAX6M,C3D10M
*ELSET,ELSET=EALL_NSM
T2D2,T3D2,
B21,B31,
SAX1,
C3D8R,
CPE4R,CPS4R,CAX4R,
S4R,M3D4R,
***MASS,
CPE3,CPS3,C3D4,C3D6,CAX3,M3D3,S3R,
CPS6M,CPE6M,CAX6M,C3D10M
*ELSET,ELSET=EALL_CONT
SAX1,
C3D8R,
CPE4R,CPS4R,CAX4R,
CPE3,CPS3,C3D4,C3D6,CAX3,
CPS6M,CPE6M,CAX6M,C3D10M
*ELSET,ELSET=EALL_AREA
S4R,M3D4R,M3D3,S3R,
*ELSET,ELSET=EALL_LEN1
T2D2,T3D2,
*ELSET,ELSET=EALL_LEN2
B21,B31,
**
*ELSET,ELSET=E2D
T2D2,
B21,
CPE4R,CPS4R,CPE3,CPS3,C3D4,C3D6,
CPS6M,CPE6M
*ELSET,ELSET=E3D
T3D2,
B31,
C3D8R,C3D4,C3D6,C3D10M,
S4R,S3R,
M3D4R,M3D3,
MASS,
**
*ELSET,ELSET=AXI
SAX1,
CAX4R,
CAX3,
CAX6M
**
*NSET,NSET=NALL
T2D2,T3D2,B21,B31,SAX1,
C3D8R,CPE4R,CPS4R,CAX4R,
S4R,M3D4R,MASS,
CPE3,CPS3,C3D4,C3D6,CAX3,
S3R,M3D3,
CPS6M,CPE6M,CAX6M,C3D10M
**
*BOUNDARY
NALL,1,3
**
*RESTART,TIMEMARKS=YES,WRITE,NUMBER=1
*AMPLITUDE,NAME=RAMP
0.,0.,0.00025,1.,0.0005,1.
**
**
*NSET, NSET=QA_TEST_ALLNODES
T2D2, T3D2, B21, B31
SAX1, C3D8R, CPE4R, CPS4R
CAX4R, S4R, M3D4R, MASS
CPE3, CPS3, C3D4, C3D6
CAX3, S3R, M3D3, CPS6M
CPE6M, CAX6M, C3D10M, NALL
**
***
*** Add NonStructural mass to bring the density 
*** effectively upto 7850.0 for all elements.
***
*Nonstructural mass,Elset=EALL_NSM,Units=masspervolume
1962.5,
*Nonstructural mass,elset=EALL_NSM,units=totalmass,
distribution=volumeproportional
0.369567E7
*Nonstructural mass,Elset=EALL_AREA,Units=massperarea
392.5
*Nonstructural mass,Elset=EALL_LEN1,Units=massperlength
392.5
*Nonstructural mass,Elset=EALL_LEN2,Units=massperlength
78.5
*STEP
*DYNAMIC,EXPLICIT
,.0005
*DLOAD,AMP=RAMP
E2D,GRAV,10.,0.,-1.,0.
E3D,GRAV,10.,0.,-1.,0.
AXI,GRAV,10.,0.,-1.,0.
*FILE OUTPUT,NUM=1, TIMEMARKS=YES
*NODE FILE
RF,
***OUTPUT,HISTORY,TIME INTERVAL=0.00005
***NODE OUTPUT,NSET=NALL
**RF,
*OUTPUT, FIELD, NUMBER INTERVAL=1, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
RF,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.0005
*END STEP
**
*STEP
*DYNAMIC,EXPLICIT
,.0005
*DLOAD,AMP=RAMP,OP=NEW
E2D,GRAV,10.,-1.,0.,0.
E3D,GRAV,10.,-1.,0.,0.
AXI,BR,78500.
*FILE OUTPUT,NUM=1, TIMEMARKS=YES
*NODE FILE
RF,
*OUTPUT, FIELD, NUMBER INTERVAL=1, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
RF,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.0005
*END STEP
**
*STEP
*DYNAMIC,EXPLICIT
,.0005
*DLOAD,AMP=RAMP,OP=NEW
E3D,GRAV,10.,0.,0.,-1.
*FILE OUTPUT,NUM=1, TIMEMARKS=YES
*NODE FILE
RF,
*OUTPUT, FIELD, NUMBER INTERVAL=1, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
RF,
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.0005
*END STEP