*HEADING
VERIFICATION OF SURFACES LOADED WITH UNIFORM PRESSURE LOAD (DSLOAD).
ANALYTICAL RESULTS:
	CPE4R	- 400000
	CPS4R	- 400000	
	C3D8R	- 400000
	CPE3	- 400000	
	CPS3	- 400000	
	C3D4	- 400000	
	C3D6	- 400000	
	S4R	- 400000	
	S3R	- 400000	
	M3D4R	- 400000	
	M3D3	- 400000	
	SFM3D4R	- 400000	
	SFM3D3	- 400000		
	CAX3	- 28900000		
	CAX4R	- 28900000
	SAX1	- 31420000
	B21	- 200000
	B31	- 200000
	R2D2	- 200000
	RAX2	- 32670000
	R3D3	- 400000
	R3D4	- 400000
*NODE,nset=b21
10005,4.0,0.0
10006,4.0,2.0
*NODE,nset=b22
10007,5.0,0.0
100071,5.0,1.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=sfm3d4r
10701,30.0,0.0
10702,32.0,0.0
10703,32.0,2.0
10704,30.0,2.0
*NODE,nset=r2d2
10601,3.0,0.0
10602,3.0,2.0
*NODE,nset=rax2
10603,26.0,0.0
10604,26.0,2.0
*NODE,nset=r3d3
10605,27.0,3.0
10606,29.0,3.0
10607,29.0,5.0
10608,27.0,5.0
*NODE,nset=r3d4
10609,27.0,0.0
10610,29.0,0.0
10611,29.0,2.0
10612,27.0,2.0
*node,nset=rbrns
10613,25.0,1.0
10614,26.0,1.0
10615,28.0,1.0
10616,28.0,4.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=sfm3d3
20601,30.0,3.0
20602,32.0,3.0
20603,32.0,5.0
20604,30.0,5.0
**  BEAM ELEMENT
*ELEMENT, TYPE=B21, ELSET=B2D1ELAS
 10003, 10005, 10006
*ELEMENT, TYPE=B22, ELSET=B2D2ELAS
 10004, 10007, 100071, 10008
*ELEMENT, TYPE=CPE4R, ELSET=PE2DELAS
 10101, 10101, 10102, 10103, 10104
*ELEMENT, TYPE=CPS4R, ELSET=PS2DELAS
 10201, 10201, 10202, 10203, 10204
*ELEMENT, TYPE=CAX4R, ELSET=CAXELAS
10301, 10301, 10302, 10303, 10304
*ELEMENT, TYPE=C3D8R, ELSET=S3DELAS
 10011, 10011, 10012, 10013, 10014, 10015, 10016, 10017, 10018
*ELEMENT, TYPE=S4R, ELSET=SHEELAS
 10401, 10401, 10402, 10403, 10404
*ELEMENT, TYPE=SAX1, ELSET=SAXIELAS
 10005, 10009, 10010
**  Membrane
*ELEMENT, TYPE=M3D4R, ELSET=MEMBELAS
 10501, 10501, 10502, 10503, 10504
*ELEMENT, TYPE=SFM3D4R, ELSET=SURF
 10701, 10701, 10702, 10703, 10704
**Triangular elements
** Plane strain
*ELEMENT,TYPE=CPE3,ELSET=TRI2DPE
20101,20101,20102,20104
20102,20104,20102,20103
*ELEMENT, TYPE=CAX3,ELSET=TRI2DA
20301,20301,20302,20304
20302,20304,20302,20303
** Plane stress
*ELEMENT,TYPE=CPS3,ELSET=TRI2DPS
20201,20201,20202,20204
20202,20204,20202,20203
** 3D Triangular elemets, C3D4
*ELEMENT, TYPE=C3D4, ELSET=TRI3D1
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
** 3D Triangular elemets, C3D6
*ELEMENT, TYPE=C3D6, ELSET=TRI3D2
30011,30014,30013,30018,30011,30012,30015
30012,30013,30017,30018,30012,30016,30015
** 2D Axisymmetric Triangular elements
** 2D Triangular membrane elements
*ELEMENT,TYPE=M3D3,ELSET=TRI2DM
20501,20501,20502,20504
20502,20504,20502,20503
** 2D Triangular surface elements
*ELEMENT,TYPE=SFM3D3,ELSET=TRISURF
20601,20601,20602,20604
20602,20604,20602,20603
** 2D Triangular shell elements
*ELEMENT,TYPE=S3R,ELSET=TRI2DS
20401,20401,20402,20404
20402,20404,20402,20403
*ELEMENT,TYPE=R2D2,ELSET=RBR2D2
10601,10601,10602
*ELEMENT,TYPE=RAX2,ELSET=RBRAX2
10602,10603,10604
*ELEMENT,TYPE=R3D3,ELSET=RBR3D3
10603,10605,10606,10608
10604,10608,10606,10607
*ELEMENT,TYPE=R3D4,ELSET=RBR3D4
10605,10609,10610,10611,10612
*MATERIAL,NAME=A1ELAS
*ELASTIC
193.1E9,0.3
*DENSITY
785, 
*BEAM SECTION,SECTION=RECT,ELSET=B2D1ELAS,MATERIAL=A1ELAS
0.2,0.2
*BEAM SECTION,SECTION=RECT,ELSET=B2D2ELAS,MATERIAL=A1ELAS
0.2,0.2
*SOLID SECTION,ELSET=S3DELAS,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=PE2DELAS,MATERIAL=A1ELAS
2.0, 
*SOLID SECTION,ELSET=PS2DELAS,MATERIAL=A1ELAS
2.0, 
*SHELL SECTION,ELSET=SHEELAS,MATERIAL=A1ELAS
0.2, 3
*SHELL SECTION,ELSET=SAXIELAS,MATERIAL=A1ELAS
0.02, 3
*SHELL SECTION,ELSET=TRI2DS,MATERIAL=A1ELAS
0.2, 3
*SOLID SECTION,ELSET=CAXELAS,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=TRI2DA,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=TRI3D1,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=TRI3D2,MATERIAL=A1ELAS
*SOLID SECTION,ELSET=TRI2DPE,MATERIAL=A1ELAS
2.0, 
*SOLID SECTION,ELSET=TRI2DPS,MATERIAL=A1ELAS
2.0, 
*MEMBRANE SECTION,ELSET=MEMBELAS,MATERIAL=A1ELAS
0.2, 
*MEMBRANE SECTION,ELSET=TRI2DM,MATERIAL=A1ELAS
0.2, 
*SURFACE SECTION, ELSET=SURF, DENSITY = 1570
*SURFACE SECTION, ELSET=TRISURF, DENSITY = 1570
*ELSET,ELSET=EALL
10011,10401,10501,20101,20102,20201,20202,20011,20401,20402
10101,10201,20012,20013,20014,20015,30011,30012,20501,10005
10401,10003,10004,10601,10602,10603,10604,10605,10301,20301
*NSET,NSET=BOT
10005,10006,10009,10011,10012,10013,10014,10101,10102
10201,10202,10301,10302,10401,10402,10501,10502
10701,10702,20101,20102,20201,20202
20011,20012,20013,20014,20301,20302,20401,20402,20501,20502
30011,30012,30013,30014,10613,10614,10615,10616,20601,20602
*NSET,NSET=TOP
10007,10008,10010,10015,10016,10017,10018,10103,10104
10203,10204,10303,10304,10403,10404,10503,10504
10703,10704,20103,20104,20203,20204
20015,20016,20017,20018,20303,20304,20403,20404,20503,20504
30015,30016,30017,30018,20603,20604
*BOUNDARY
BOT,1,3
TOP,1,3
*nset,nset=nall
top,bot
*RESTART,WRITE,NUMBER=1
*AMPLITUDE, NAME=RAMP
0.,0.,0.0005,1.0,0.001,1.0
*SURFACE,TYPE=ELEMENT, NAME=SURFACE_2D
10101,s3
10201,s3
20102,s3
20202,s3
*SURFACE,TYPE=ELEMENT, NAME=SURFACE_3D
10011,s2
10401, SNEG
10501, SNEG
10701, SNEG
20013,s3
20014,s3
30012,s4
20502, SNEG
20602, SNEG
20501, SNEG
20601, SNEG
20401, SNEG
20402, SNEG
10603, SNEG
10604, SNEG
10605, SNEG
*SURFACE,TYPE=ELEMENT, NAME=SURFACE_AX
20302,s3
10301,s3
10005,SNEG
10602,SNEG
*SURFACE,TYPE=ELEMENT, NAME=S20101P3
 20101, S3
*SURFACE,TYPE=ELEMENT, NAME=S10101P4
 10101, S4
*SURFACE,TYPE=ELEMENT, NAME=S20011P4
 20011, S4
*SURFACE,TYPE=ELEMENT, NAME=S20201P3
 20201, S3
*SURFACE,TYPE=ELEMENT, NAME=S10201P4
 10201, S4
*SURFACE,TYPE=ELEMENT, NAME=S10301P4
 10301, S4
*SURFACE,TYPE=ELEMENT, NAME=S20301P3
 20301, S3
*SURFACE,TYPE=ELEMENT, NAME=S10011P6
 10011, S6
*SURFACE,TYPE=ELEMENT, NAME=S30011P5
 30011, S5
*SURFACE,TYPE=ELEMENT, NAME=S20014P4
 20014, S4
*SURFACE,TYPE=ELEMENT, NAME=S10005P
 10005, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S10601P
 10601, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S10602P
 10602, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S10603P
 10603, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S10604P
 10604, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S10605P
 10605, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S20502P
 20502, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S20602P
 20602, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S20601P
 20601, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S10401P
 10401, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S10501P
 10501, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S10701P
 10701, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S20501P
 20501, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S20401P
 20401, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S20402P
 20402, SNEG
*SURFACE,TYPE=ELEMENT, NAME=S10003P
 10003, S1
*SURFACE,TYPE=ELEMENT, NAME=S10004P
 10004, S1
*RIGID BODY,ELSET=RBR2D2,REF NODE=10613
*RIGID BODY,ELSET=RBRAX2,REF NODE=10614
*RIGID BODY,ELSET=RBR3D3,REF NODE=10615
*RIGID BODY,ELSET=RBR3D4,REF NODE=10616
*BOUNDARY
10613, 1,6 
10614, 1,6 
10615, 1,6 
10616, 1,6 
**
*NSET, NSET=QA_TEST_ALLNODES
b21, b22, sax1, c3d8r
cpe4r, cps4r, cax4r, s4r
m3d4r, sfm3d4r,  r2d2, rax2, r3d3
sfm3d3, r3d4, rbrns, cpe3, cps3
c3d4, c3d6, cax3, s3r
m3d3, BOT, TOP, nall
**
*STEP
*DYNAMIC,EXPLICIT
,0.001
********LOAD ON HYBRID SURFACES********
*DSLOAD,AMP=RAMP
SURFACE_2D,P,100000
SURFACE_3D,P,100000
SURFACE_AX,P,100000
S10601P ,P,-100000.
**************************
*FILE OUTPUT, NUM=1
*NODE FILE
RF, 
***OUTPUT,HISTORY,TIME INTERVAL=0.0
***NODE OUTPUT,nset=nall
**rf,
*OUTPUT, FIELD, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
RF, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.001
*END STEP
*STEP
*DYNAMIC,EXPLICIT
,0.001
*DSLOAD,AMP=RAMP,OP=NEW
 S10101P4,p, 100000.
 S10201P4,P, 100000.
 S10011P6,p, 100000.
 S10401P, p,-100000.
 S10501P, p,-100000.
 S10701P, p,-100000.
 S20101P3,p, 100000.
 S20201P3,P, 100000.
 S20011P4,P, 100000.
 S20014P4,P, 100000.
 S30011P5,P, 100000.
 S20502P, P,-100000.
 S20602P, P,-100000.
 S20501P, P,-100000.
 S20601P ,P,-100000.
 S20401P, P,-100000.
 S20402P, P,-100000.
 S20301P3,P, 100000.
 S10301p4,P, 100000.
 S10005P, P,-100000.
 S10601P, P,100000.
 S10602P, P,-100000.
 S10603P, P,-100000.
 S10604P, P,-100000.
 S10605P, P,-100000.
 S10003P, P,-100000.
 S10004P, P,-100000.
*FILE OUTPUT, NUM=1
*NODE FILE
RF, 
*OUTPUT, FIELD, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
RF, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.001
*END STEP