*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
*NODE,nset=c3d8r
*NODE,nset=cpe4r
*NODE,nset=cps4r
*NODE,nset=cax4r
*NODE,nset=s4r
*NODE,nset=m3d4r
*NODE,nset=sfm3d4r
*node,nset=cpe3
*node,nset=cps3
*node,nset=c3d4
*node,nset=c3d6
*node,nset=cax3
*node,nset=s3r
*node,nset=m3d3
*node,nset=sfm3d3
*NODE, NSET=REFVEL
*ELEMENT, TYPE=CPE4R, ELSET=PE2DELAS
*ELEMENT, TYPE=CPS4R, ELSET=PS2DELAS
*ELEMENT, TYPE=CAX4R, ELSET=CAXELAS
*ELEMENT, TYPE=C3D8R, ELSET=S3DELAS
*ELEMENT, TYPE=S4R, ELSET=SHEELAS
*ELEMENT, TYPE=SAX1, ELSET=SAXIELAS
**  Membrane
*ELEMENT, TYPE=M3D4R, ELSET=MEMBELAS
**  SURFACE ELEMENTS
*ELEMENT, TYPE=SFM3D4R, ELSET=SURF
**Triangular elements
** Plane strain
*ELEMENT,TYPE=CPE3,ELSET=TRI2DPE
*ELEMENT, TYPE=CAX3,ELSET=TRI2DA
** Plane stress
*ELEMENT,TYPE=CPS3,ELSET=TRI2DPS
** 3D Triangular elemets, C3D4
*ELEMENT, TYPE=C3D4, ELSET=TRI3D1
** 3D Triangular elemets, C3D6
*ELEMENT, TYPE=C3D6, ELSET=TRI3D2
** 2D Axisymmetric Triangular elements
** 2D Triangular membrane elements
*ELEMENT,TYPE=M3D3,ELSET=TRI2DM
*ELEMENT,TYPE=SFM3D3,ELSET=TRISURF
** 2D Triangular shell elements
*ELEMENT,TYPE=S3R,ELSET=TRI2DS
*ELEMENT,TYPE=MASS,ELSET=MASS
*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
*NSET,NSET=DIR2
*NSET,NSET=DIR3
*NSET,NSET=DIR1
*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