*HEADING
: DYNAMIC PRESSURE AND CLOSED-END BUOYANCY [PIPE31]
GRIDDED WAVES (LINEAR INTERPOLATION) AND *DYNAMIC PROCEDURE
 -- 50 ELEMENTS WITH PB AND PARTIAL SUBMERGENCE
 -- 1 ELEMENT WITH TSB
EPR SOLUTION: 8307 LBS., TOTAL MAX. VERTICAL LOAD
ABAQUS SOLUTION: 8297 LBS., TOTAL MAX. VERTICAL LOAD
***********
** /Post journal file:
**
** restart,fil=.res
** results,fil=.fil
** read curve,name=rf_pb ,var=rf3,node=n_pb,sum
** read curve,name=rf_tsb,var=rf3,node=500
** disp
** rf_pb
** rf_tsb
**
***********
***********
** Aqua parameters:
**
**   Water depth    - 1000
**   Seabed         - 100
**   Still water    - 1100
**   Steady current - (0,0,0) No steady current
**   Density        - 2.0
**   Gravity const  - 32.2
**
*AQUA
 100., 1100., 32.2, 2.0
 0., 0., 0., 0.
 0., 0., 0., 1100.
***********
***********
** Wave parameters (Gridded):
**
**     Wave amplitude - 50 (height = 100)
**     Wave period    - 15
**     Phase angle    - 0
**     Wave direction - (1,0,0)
**
*WAVE,TYPE=GRIDDED,DATA FILE=gridwave_3d.binary
0.,0.,0.,1.,0.
***********
***********
*NODE
   1,0.,0.,1000.
  51,0.,0.,1175.
*NODE,NSET=N_TSB
 500,0.,0.,1000.
 501,0.,0.,1001.
*NGEN,NSET=N_PB
 1,51,1
*ELEMENT,TYPE=PIPE31
 1,1,2
*ELGEN,ELSET=E_PB
 1,50,1,1
*BOUNDARY
 N_PB,1,6
*RESTART,WRITE,FREQ=10
*BEAM SECTION,SECTION=PIPE,MATERIAL=MAT1,ELSET=E_PB
 1.0,0.25
 0.,-1.,0.
*ELEMENT,TYPE=B31,ELSET=CLOAD
 500,500,501
*BEAM SECTION,SECTION=CIRC,MATERIAL=MAT1,ELSET=CLOAD
 1.0,
 0.,-1.,0.
*BOUNDARY
 500,1,6
*MATERIAL,NAME=MAT1
*ELASTIC
 1.E6,
*DENSITY
 1.0E6, 
** -------------------------------------------
**     Dynamic pressure test.
** -------------------------------------------
*STEP,NLGEOM,INC=1000
 Dynamic pressure test.
*DYNAMIC,DIRECT
 .2,15
*DLOAD
 E_PB,PB, 1.0  , 1.1283792
**el,PB, amag0, Diameter (2/sqrt(pi) = area of 1)
*CLOAD
 500,  TSB,  1.0, 1.0, 0.,0.,-1.
**node,TSB,amag0,area,  normal
*NODE PRINT,TOTAL=YES,NSET=N_PB
 RF, 
*NODE PRINT,NSET=N_TSB
 RF, 
*NODE FILE,FREQ=10
 RF, 
*ENERGY FILE,FREQ=10
*EL PRINT,FREQ=0
*END STEP