*HEADING
VERIFICATION OF ELEMENTS LOADED WITH VISCOUS PRESSURE (*DSLOAD).
ANALYTICAL RESULTS:
	CPE4R	- 400000
	CPS4R	- 400000	
	C3D8R	- 400000
	CPE3	- 400000	
	CPS3	- 400000	
	C3D4	- 400000	
	C3D6	- 400000	
	S4R	- 400000	
	S3R	- 400000	
	M3D4R	- 400000
        SFM3D4R - 400000	
	M3D3	- 400000
        SFM3D3  - 400000		
	CAX3	- 28900000		
	CAX4R	- 28900000
	SAX1	- 31420000
*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
10601,26.0,0.0
10602,28.0,0.0
10603,28.0,2.0
10604,26.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=sfm3d3
20601,26.0,3.0
20602,28.0,3.0
20603,28.0,5.0
20604,26.0,5.0
*NODE, NSET=REFVEL
99999, 0.0, 0.0, 0.0
*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
**  SURFACE ELEMENTS
*ELEMENT, TYPE=SFM3D4R, ELSET=SURF
 10601, 10601, 10602, 10603, 10604
**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
*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=MASS,ELSET=MASS
99999, 99999
*MATERIAL,NAME=A1ELAS
*ELASTIC
193.1E9,0.3
*DENSITY
7850, 
*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, POISSON=0.3
0.2, 
*MEMBRANE SECTION,ELSET=TRI2DM,MATERIAL=A1ELAS, POISSON=0.3
0.2, 
*MASS, ELSET=MASS
1.0E-6
*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,10301,20301
*NSET,NSET=DIR2
10101,10102,10103,10104,10201,10202
10203,10204,10301,10302,10303,10304,10011,10012,10013,10014
10015,10016,10017,10018,20201,20202,20203,20204,20301,20302
20303,20304,20201,20202,20203,20204,
20011,20012,20013,20014,20015,20016,20017,20018
30011,30012,30013,30014,30015,30016,30017,30018
20101,20102,20103,20104
*NSET,NSET=DIR3
10401,10402,10403,10404,10501,10502,10503,10504
10601,10602,10603,10604,20601,20602,20603,20604
20501,20502,20503,20504,20401,20402,20403,20404
*NSET,NSET=DIR1
10009,10010
*nset,nset=nall
dir1,dir2,dir3
*BOUNDARY
dir2,1,1
dir2,3,3
dir3,1,2
dir1,2,3
REFVEL,1,3
*INITIAL CONDITIONS,TYPE=VELOCITY
DIR2,2,100.
DIR1,1,100.
DIR3,3,100.
*AMPLITUDE, NAME=RAMP
0.,0.,0.0002,1.0
*RESTART,TIMEMARKS=YES,WRITE,NUMBER=1
*SURFACE,TYPE=ELEMENT, NAME=SURF_2D
 10101, S3
 10201, S3
 20102, S3
 20202, S3
*SURFACE,TYPE=ELEMENT, NAME=SURF_3D
 10601, SNEG
 10011, S2
 20013, S3
 20014, S3
 30012, S4
 10401, SNEG
 20401, SNEG
 10501, SNEG
 20402, SNEG
 20501, SNEG
 20601, SNEG
 20502, SNEG
*SURFACE,TYPE=ELEMENT, NAME=SURF_AX
 10301, S3
 20302, S3
 10005, SNEG
**
*NSET, NSET=QA_TEST_ALLNODES
sax1, c3d8r, cpe4r, cps4r
cax4r, s4r, m3d4r, sfm3d4r, cpe3
cps3, c3d4, c3d6, cax3
s3r, m3d3, DIR2, DIR3
DIR1, nall
**
*STEP
*DYNAMIC,EXPLICIT
,0.0002
*DSLOAD, REF NODE=REFVEL
SURF_AX,VP,1000.
SURF_2D,VP,1000.
SURF_3D,VP,1000.
*FILE OUTPUT,TIMEMARKS=YES, NUM=1
*NODE FILE
U,RF
***OUTPUT,HISTORY,TIME INTERVAL=0.0
***ELEMENT OUTPUT,ELSET=EALL
**S,
***NODE OUTPUT,nset=nall
**u,rf
***ENERGY OUTPUT
**allke,allie,allae,allse,allvd,allwk,etotal
*OUTPUT, FIELD, TIMEMARKS=YES, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
U, RF
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.0002
*END STEP