*HEADING
VERIFICATION OF ELEMENTS LOADED WITH NON-UNIFORM BODY FORCES. 
ANALYTICAL RESULTS:
	CPE4R 	- 800000
	CPS4R 	- 800000
	C3D8R 	- 800000
	CPE3 	- 800000	
	CPS3 	- 800000	
	C3D4 	- 800000	
	C3D6 	- 800000	
	S4R	- 80000
	S3R	- 80000	
	M3D4	- 80000	
	M3D3	- 80000	
	CAX4R	- 57805305
	CAX3	- 57805305
	SAX1	- 628318.53
	T2D2	- 40000
	T3D2	- 40000
*NODE,NSET=T2D2
10005,4.0,0.0
10006,4.0,2.0
*NODE,NSET=T3D2
10007,5.0,0.0
10008,5.0,2.0
*NODE,NSET=SAX1
10009,25.0,0.0
10010,25.0,2.0
*NODE,NSET=C3D8R
10011,7.0,0.0,2.
10012,9.0,0.0,2.
10013,9.0,0.0,0.0
10014,7.0,0.0,0.0
10015,7.0,2.0,2.
10016,9.0,2.0,2.
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
10302,24.,0.0
10303,24.,2.0
10304,22.,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=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.
20012,9.0,3.0,2.
20013,9.0,3.0,0.0
20014,7.0,3.0,0.0
20015,7.0,5.0,2.
20016,9.0,5.0,2.
20017,9.0,5.0,0.0
20018,7.0,5.0,0.0
*NODE,NSET=C3D6
30011,7.0,6.0,2.
30012,9.0,6.0,2.
30013,9.0,6.0,0.0
30014,7.0,6.0,0.0
30015,7.0,8.0,2.
30016,9.0,8.0,2.
30017,9.0,8.0,0.0
30018,7.0,8.0,0.0
*NODE,NSET=CAX3
20301,22.,3.0
20302,24.,3.0
20303,24.,5.0
20304,22.,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
*NODE,NSET=B21
40041, 31.0, 0.0
40042, 32.0, 1.0
*NODE,NSET=B31
40043, 31.0, 0.0, 0.0
40044, 32.0, 1.0, 1.0
*NODE,NSET=B22
40045, 33.0, 0.0
40046, 33.5, 0.5
40047, 34.0, 1.0
*NODE,NSET=B32
40048, 35.0, 0.0, 0.0
40049, 35.5, 0.5, 0.5
40050, 36.0, 1.0, 1.0
**
*ELEMENT,TYPE=T2D2,ELSET=T2D2
10003,10005,10006
*ELEMENT,TYPE=T3D2,ELSET=T3D2
10004,10007,10008
*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=C3D8R,ELSET=C3D8R
10011,10011,10012,10013,10014,10015,10016,10017,10018
*ELEMENT,TYPE=S4R,ELSET=S4R
10401,10401,10402,10403,10404
*ELEMENT,TYPE=SAX1,ELSET=SAX1
10005,10009,10010
*ELEMENT,TYPE=M3D4R,ELSET=M3D4R
10501,10501,10502,10503,10504
*ELEMENT,TYPE=CPE3,ELSET=CPE3
20101,20101,20102,20104
20102,20104,20102,20103
*ELEMENT,TYPE=CAX3,ELSET=CAX3
20301,20301,20302,20304
20302,20304,20302,20303
*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=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,
*ELEMENT,TYPE=B21,ELSET=B21
40041, 40041, 40042
*ELEMENT,TYPE=B31,ELSET=B31
40043, 40043, 40044
*ELEMENT,TYPE=B22,ELSET=B22
40045, 40045, 40046, 40047
*ELEMENT,TYPE=B32,ELSET=B32
40048, 40048, 40049, 40050
**
*SOLID SECTION,ELSET=T2D2,MATERIAL=A1ELAS
0.2,
*SOLID SECTION,ELSET=T3D2,MATERIAL=A1ELAS
0.2,
*SOLID SECTION,ELSET=C3D8R,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=CPE4R,MATERIAL=A1ELAS
2.0,
*SOLID SECTION,ELSET=CPS4R,MATERIAL=A1ELAS
2.0,
*SHELL SECTION,ELSET=S4R,MATERIAL=A1ELAS
0.2,3
*SHELL SECTION,ELSET=SAX1,MATERIAL=A1ELAS
0.02,3
*SHELL SECTION,ELSET=S3R,MATERIAL=A1ELAS
0.2,3
*SOLID SECTION,ELSET=CAX4R,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=CAX3,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=C3D4,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=C3D6,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=CPE3,MATERIAL=A1ELAS
2.0,
*SOLID SECTION,ELSET=CPS3,MATERIAL=A1ELAS
2.0,
*MEMBRANE SECTION,ELSET=M3D4R,MATERIAL=A1ELAS
0.2,
*MEMBRANE SECTION,ELSET=M3D3,MATERIAL=A1ELAS
0.2,
*SOLID SECTION,ELSET=CPS6M,MATERIAL=A1ELAS
2.0
*SOLID SECTION,ELSET=CPE6M,MATERIAL=A1ELAS
2.0
*SOLID SECTION,ELSET=CAX6M,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=C3D10M,MATERIAL=A1ELAS
*BEAM SECTION,ELSET=B21,MATERIAL=A1ELAS,SECTION=CIRC
0.05,
*BEAM SECTION,ELSET=B22,MATERIAL=A1ELAS,SECTION=CIRC
0.05,
*BEAM SECTION,ELSET=B31,MATERIAL=A1ELAS,SECTION=CIRC
0.05,
*BEAM SECTION,ELSET=B32,MATERIAL=A1ELAS,SECTION=CIRC
0.05,
**
*MATERIAL,NAME=A1ELAS
*ELASTIC
193.1E9,0.3
*DENSITY
785,
**
*ELSET,ELSET=EALL
T2D2,T3D2,C3D8R,CPE4R,CPS4R,
S4R,SAX1,S3R,CAX4R,CAX3,C3D4,
C3D6,CPE3,CPS3,M3D4R,M3D3,
CPS6M,CPE6M,CAX6M,C3D10M,
B21,B22,B31,B32
*ELSET,ELSET=E2D
T2D2,
CPE4R,CPS4R,
CPE3,CPS3,
CPS6M,CPE6M
*ELSET,ELSET=BEAM2D
B21,B22
*ELSET,ELSET=E3D
T3D2,
C3D8R,C3D4,C3D6,C3D10M
S4R,S3R,
M3D4R,M3D3
*ELSET, ELSET=BEAM3D
B31,B32
*ELSET,ELSET=AXI
SAX1,CAX4R,CAX3,CAX6M,
**
*NSET,NSET=FIX2
T2D2,SAX1,
CPE4R,CPS4R,CAX4R,
CPE3,CPS3,CAX3,B21,B22
*NSET,NSET=FIX3
T3D2,
C3D8R,C3D4,C3D6,
S3R,S4R,M3D4R,M3D3,
B31,B32
*NSET,NSET=NALL
FIX2,FIX3,
CPS6M,CPE6M,CAX6M,C3D10M
**
*BOUNDARY
FIX2,1,2
FIX3,1,3
**
*NSET,NSET=REACTION
FIX2,FIX3
*NSET,NSET=ACCEL
CPS6M,CPE6M,CAX6M,C3D10M
**
*NSET,NSET=QA_TEST
ACCEL,
*RESTART,TIMEMARKS=YES,WRITE,NUMBER=1
*AMPLITUDE,NAME=RAMP1
0.,0.,0.0005,1.0,0.001,1.0
*AMPLITUDE,NAME=RAMP2
0.,0.,0.0005,1.0,0.001,1.0
*AMPLITUDE,NAME=RAMP3
0.,0.,0.0005,1.0,0.001,1.0
*STEP
*DYNAMIC,EXPLICIT
,0.001
*DLOAD,AMP=RAMP1
E2D,BXNU,100000.
E3D,BXNU,100000.
AXI,BRNU,100000.
BEAM2D,PXNU,25000.
BEAM2D,P2NU,10000.
BEAM3D,PXNU,25000.
BEAM3D,P2NU,10000.
*FILE OUTPUT,NUM=1, TIMEMARKS=YES
*NODE FILE,NSET=REACTION
RF,
*NODE FILE,NSET=ACCEL
U,
********************************************
*OUTPUT,FIELD,NUMBER INTERVAL=1,TIMEMARKS=YES
*NODE OUTPUT,NSET=QA_TEST
U,
*OUTPUT,HISTORY,TIME INTERVAL=0.001
*ENERGY OUTPUT, VAR=PRESELECT
**********************
***OUTPUT,HISTORY,TIME INTERVAL=0.0
***NODE OUTPUT,NSET=NALL
**RF,
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLAE,ALLSE,ALLVD,ALLWK,ETOTAL
*END STEP
*STEP
*DYNAMIC,EXPLICIT
,0.001
*DLOAD,AMP=RAMP2,OP=NEW
E2D,BYNU,100000.
E3D,BYNU,100000.
AXI,BZNU,100000.
BEAM2D,PYNU,25000.
BEAM3D,PYNU,25000.
BEAM3D,P1NU,10000.
*END STEP
*STEP
*DYNAMIC,EXPLICIT
,0.001
*DLOAD,AMP=RAMP3,OP=NEW
E3D,BZNU,100000.
BEAM3D,PZNU,25000.
*END STEP