*HEADING
  P AND VP TYPE *DLOAD AND *DSLOAD VERIFICATION FOR COH3D8, COH3D6, COH2D4 & COHAX4
*PREPRINT, MODEL=YES
*NODE,NSET=COH3D8
11, 0., 0., 0.
12, 0., 0., 7.
13, 7., 0., 7.
14, 7., 0., 0.
15, 0., 7., 0.
16, 0., 7., 7.
17, 7., 7., 7.
18, 7., 7., 0.
*SYSTEM
0., 10., 0.
*NODE,NSET=COH3D6
21, 0., 0., 0.
22, 0., 0., 7.
23, 7., 0., 7.
24, 7., 0., 0.
25, 0., 7., 0.
26, 0., 7., 7.
27, 7., 7., 7.
28, 7., 7., 0.
*SYSTEM
0., 20., 0.
*NODE,NSET=COH2D4
31, 0., 0., 0.
32, 7., 0., 0.
33, 7., 7., 0.
34, 0., 7., 0.
*SYSTEM
0., 30., 0.
*NODE,NSET=COHAX4
41, 10., 0., 0.
42, 20., 0., 0.
43, 20., 10., 0.
44, 10., 10., 0.
*ELEMENT,TYPE=COH3D8, ELSET=COH3D8
11, 11, 12, 13, 14, 15, 16, 17, 18
*ELEMENT,TYPE=COH3D6, ELSET=COH3D6
21, 21, 22, 23, 25, 26, 27
22, 24, 21, 23, 28, 25, 27
*ELEMENT,TYPE=COH2D4, ELSET=COH2D4
31, 31, 32, 33, 34
*ELEMENT,TYPE=COHAX4, ELSET=COHAX4
41, 41, 42, 43, 44
*NSET, NSET=ALL
COH3D8, COH3D6, COH2D4, COHAX4
*ELSET, ELSET=EALL
COH3D8, COH3D6, COH2D4, COHAX4
*COHESIVE SECTION, THICKNESS=GEOMETRY, MATERIAL=COH, ELSET=EALL, RESPONSE=CONTINUUM
*MATERIAL,NAME=COH
*ELASTIC
3.E6, 0.3
*DENSITY
 10.0, 
*BOUNDARY
 COH3D8,1,3
 COH3D6,1,3
 COH2D4,1,2
 COHAX4,1,2
*NSET, NSET=QA_TEST_ALL
ALL
*ELSET, ELSET=QA_TEST_EALL
EALL
*ELSET, ELSET=COH2D
COH2D4, COHAX4
*SURFACE, NAME=SIDE1_3D
COH3D8, S1
COH3D6, S1
*SURFACE, NAME=SIDE1_2D
COH2D4, S1
*SURFACE, NAME=SIDE1_AX
COHAX4, S1
*SURFACE, NAME=SIDE2_3D
COH3D8, S2
COH3D6, S2
*SURFACE, NAME=SIDE2_2D
COH2D4, S3
*SURFACE, NAME=SIDE2_AX
COHAX4, S3
***************************
** 	STEP-1
***************************
*STEP, NAME=P_DLOAD
*DYNAMIC, EXPLICIT
,1.
*DLOAD
 COH3D8, P1, 3000.
 COH3D8, P2, 3000.
 COH3D6, P1, 3000.
 COH3D6, P2, 3000.
 COH2D, P1, 3000.
 COH2D, P3, 3000.
*OUTPUT, FIELD, NUMBER INTERVAL=1, TIME MARKS=YES
*NODE OUTPUT, NSET=QA_TEST_ALL
 U, RF
*OUTPUT, HISTORY, FREQUENCY=0
*END STEP
***************************
** 	STEP-2
***************************
*STEP, NAME=P_DSLOAD
*DYNAMIC, EXPLICIT
,1.
*DSLOAD
 SIDE1_3D, P, 3000.
 SIDE1_2D, P, 3000.
 SIDE1_AX, P, 3000.
 SIDE2_3D, P, 3000.
 SIDE2_2D, P, 3000.
 SIDE2_AX, P, 3000.
*END STEP
***************************
** 	STEP-3
***************************
*STEP, NAME=VP_DLOAD
*DYNAMIC, EXPLICIT
,1.
*BOUNDARY, OP=MOD, TYPE=VELOCITY
 COH3D8, 2, 2, 100.
 COH3D6, 2, 2, 100.
 COH2D4, 2, 2, 100.
 COHAX4, 2, 2, 100.
*DLOAD
 COH3D8, VP1, 1000.
 COH3D8, VP2, 1000.
 COH3D8, VP3, 1000.
 COH3D8, VP4, 1000.
 COH3D6, VP1, 1000.
 COH3D6, VP2, 1000.
 21, VP3, 1000.
 21, VP4, 1000.
 COH2D, VP1, 1000.
 COH2D, VP2, 1000.
 COH2D, VP3, 1000.
 COH2D, VP4, 1000.
*END STEP
***************************
** 	STEP-4
***************************
*STEP, NAME=VP_DSLOAD
*DYNAMIC, EXPLICIT
,1.
*DSLOAD
 SIDE1_3D, VP, 1000.
 SIDE1_2D, VP, 1000.
 SIDE1_AX, VP, 1000.
 SIDE2_3D, VP, 1000.
 SIDE2_2D, VP, 1000.
 SIDE2_AX, VP, 1000.
*END STEP