*HEADING
: DYNAMIC PRESSURE AND CLOSED-END BUOYANCY [PIPE31]
AIRY WAVES AND *DYNAMIC PROCEDURE
 -- 50 ELEMENTS WITH PB AND PARTIAL SUBMERGENCE
 -- 1 ELEMENT WITH TSB
EPR SOLUTION: 8307 LBS., TOTAL MAX. VERTICAL LOAD
ABAQUS SOLUTION: 8310 LBS., TOTAL MAX. VERTICAL LOAD
***********
** /Post journal file:
**
** 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
***********
***********
** Wave parameters (Airy):
**
**     Wave amplitude - 50 (height = 100)
**     Wave period    - 15
**     Phase angle    - 0
**     Wave direction - (1,0,0)
**
*WAVE,TYPE=AIRY,WAVE PERIOD
***********
***********
*NODE
*NODE,NSET=N_TSB
*NGEN,NSET=N_PB
*ELEMENT,TYPE=PIPE31
*ELGEN,ELSET=E_PB
*BOUNDARY
 N_PB,1,6
*******************************************
**DESCRIBE SURFACE SECTION FOR AQUA WAVE VISUALIZATION
*******************************************
*NODE
*NGEN, NSET=NAQUA
*ELEMENT, TYPE=SFM3D4R, ELSET=AQUAVIZ
*ELGEN, ELSET = AQUAVIZ
*SURFACE SECTION, ELSET = AQUAVIZ, AQUAVISUALIZATION = YES
*******************************************
*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
*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
*OUTPUT, FIELD, NUMBER INTERVAL = 50
*NODE OUTPUT
U
*END STEP