*HEADING
Partial penetration of the element into the sea-bed: 
Elements tested: rb2d2, rb3d3,
		 t2d2, t2d2h, t2d3, t2d3h,
		 t3d2, t3d2h, t3d3, t3d3h
***********
***********
** Aqua parameters:
**
**   Seabed         - 0
**   Still water    - 2
**   Steady current - (1,0,0) No depth variation
**   Density        - 1.99
**   Gravity const  - 32.3
**
*AQUA
 0., 2., 32.2, 1.99
 1., 0., 0., 0.
 1., 0., 0., 2.
***********
***********
** Wave parameters (airy):
**
**     Wave amplitude - 0.1
**     Wave period    - 10
**     Phase angle    - 0
**     Wave direction - (1,0)
**
*wave,type=airy,wave period
.1,10.0,0.,1.,0.0
**.1,80.25,0.,1.,0.0
*NODE,nset=ntr
1,  0., 0., -1.
2,  0., 0.,  1.
3,  0., 0., -1.
4,  0., 0.,  1.
5,  0., .0, -1.
6,  0., .0,  0.
7,  0., .0,  1.
11, 0., -1.
12, 0.,  1.
13,  .0, -1.
14,  .0,  1.
15,  .0, -1.
16,  .0,  0.
17,  .0,  1.
101,  0., 0., -1.
102,  0., 0.,  1.
103,  0., 0., -1.
104,  0., 0.,  1.
105,  0., .0, -1.
106,  0., .0,  0.
107,  0., .0,  1.
111, 0., -1.
112, 0.,  1.
113,  .0, -1.
114,  .0,  1.
115,  .0, -1.
116,  .0,  0.
117,  .0,  1.
*node
991,1.
992,1.
993,1.
994,1.
995,1.
996,1.
997,1
998,1
999,1
1000,1
*ELEMENT,TYPE=rb3D2,ELSET=etr1
1, 1,2
*element,type=rb2d2,elset=etr2
11, 11,12
*element,type=t3d2
3, 3,4
*element,type=t3d3
4, 5,6,7
*element,type=t2d2
13, 13,14
*element,type=t2d3
14, 15,16,17
*element,type=t3d2h
103, 103,104
*element,type=t3d3h
104, 105,106,107
*element,type=t2d2h
113, 113,114
*element,type=t2d3h
114, 115,116,117
*elset,elset=truss
3,4,13,14,103,104,113,114
*elset,elset=etr
etr1,etr2,truss
*rigid body,elset=etr1,ref=991
*rigid body,elset=etr2,ref=992
*solid section,material=mat1,elset=truss
1.0, 
*material,name=mat1
*elastic
1.e6, 
*mpc
beam,3,993
beam,4,993
beam,5,994
beam,6,994
beam,7,994
beam,13,995
beam,14,995
beam,15,996
beam,16,996
beam,17,996
beam,103,997
beam,104,997
beam,105,998
beam,106,998
beam,107,998
beam,113,999
beam,114,999
beam,115,1000
beam,116,1000
beam,117,1000
*BOUNDARY
991, 1,6
992, 1,6
993, 1,6
994, 1,6
995, 1,6
996, 1,6
997,1,6
998,1,6
999,1,6
1000,1,6
**
** Dummy mass element to define mass in the problem
**
*ELEMENT,TYPE=MASS,ELSET=MASS
 99, 9999
*MASS,ELSET=MASS
 1., 
*amplitude,name=zero
0.,0., 1.e-6,0.0
***********************
*** S T E P   1 :  Only drag loads
***		   Analytical result: RF = 0.995 i + 0 j + 0 k
*STEP,inc=100
*static
1.e-10,1.e-10
*dload
etr,fdd,,1.0,1.0,,,zero
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*NODE PRINT
 RF, 
*END STEP
***********************
*** S T E P   2 :  Ramp down the drag load
*step
*static
1.e-10,1.e-10
*dload,op=new
*endstep
***********************
*** S T E P   3 :  Inertial load
**		   Analytical result: RF = 0.389 i + 0 j + 0 k (peak)
*STEP,inc=100
*dynamic
.5,10.5,,.5
*dload
etr,fi,,1.0,1.0,1.0
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*END STEP