*HEADING
Partial penetration of the element into the sea-bed: 
Elements tested: b21, b21h, b22, b22h, b23, b23h,
		 b31, b31h, b32, b32h, b23, b33h
***********
***********
** 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
***********
***********
*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.
21, 79.8, 0., -1.
22, 79.8, 0.,  1.
23, 79.8, 0., -1.
24, 79.8, 0.,  1.
25, 79.8, 0., -1.
26, 79.8, 0.,  0.
27, 79.8, 0.,  1.
31, 79.8, -1.
32, 79.8,  1.
33, 79.8, -1.
34, 79.8,  1.
35, 79.8, -1.
36, 79.8,  0.
37, 79.8,  1.
*node
991,1.
992,1.
993,1.
994,1.
995,1.
996,1.
981,1.
982,1.
983,1.
984,1.
985,1.
986,1.
*ELEMENT,TYPE=b31
1, 1,2
*element,type=b32
2, 5,6,7
*element,type=b33
3, 3,4
*ELEMENT,TYPE=b21
11, 11,12
*element,type=b22
12, 15,16,17
*element,type=b23
13, 13,14
*ELEMENT,TYPE=b31h
21, 21,22
*element,type=b32h
22, 25,26,27
*element,type=b33h
23, 23,24
*ELEMENT,TYPE=b21h
31, 31,32
*element,type=b22h
32, 35,36,37
*element,type=b23h
33, 33,34
*elset,elset=beams3d
1,2,3,21,22,23
*elset,elset=beams2d
11,12,13,31,32,33
*elset,elset=beams
beams2d,beams3d
*beam section,section=circ,material=mat1,elset=beams2d
1.0, 
0.,0.,-1.

*beam section,section=circ,material=mat1,elset=beams3d
1.0, 
0.,-1.,0.

*material,name=mat1
*elastic
1.e6, 
*mpc
beam,1,991
beam,2,991
beam,3,993
beam,4,993
beam,5,992
beam,6,992
beam,7,992
beam,11,994
beam,12,994
beam,13,996
beam,14,996
beam,15,995
beam,16,995
beam,17,995
beam,21,981
beam,22,981
beam,23,983
beam,24,983
beam,25,982
beam,26,982
beam,27,982
beam,31,984
beam,32,984
beam,33,986
beam,34,986
beam,35,985
beam,36,985
beam,37,985
*BOUNDARY
991, 1,6
992, 1,6
993, 1,6
994, 1,6
995, 1,6
996, 1,6
981, 1,6
982, 1,6
983, 1,6
984, 1,6
985, 1,6
986, 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
beams,fdd,,1.0,1.0,,,zero
*NODE PRINT
 RF, 
*OUTPUT,FIELD
*NODE OUTPUT
 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
beams,fi,,1.0,1.0,1.0
*END STEP