*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
 100., 1100., 32.2, 2.0
 0., 0., 0., 0.
 0., 0., 0., 1100.
***********
***********
** Wave parameters (Airy):
**
**     Wave amplitude - 50 (height = 100)
**     Wave period    - 15
**     Phase angle    - 0
**     Wave direction - (1,0,0)
**
*WAVE,TYPE=AIRY,WAVE PERIOD
 50,15.,0., 1.,0.,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
*******************************************
**DESCRIBE SURFACE SECTION FOR AQUA WAVE VISUALIZATION
*******************************************
*NODE
 9001, -200, -150, 1100
 9030, 200, -150, 1100
 9031, -200, -130, 1100
 9060, 200, -130, 1100
 9061, -200, -110, 1100
 9090, 200, -110, 1100
 9091, -200, -90, 1100
 9120, 200, -90, 1100
 9121, -200, -70, 1100
 9150, 200, -70, 1100
 9151, -200, -50, 1100
 9180, 200, -50, 1100
 9181, -200, -30, 1100
 9210, 200, -30, 1100
 9211, -200, -10, 1100
 9240, 200, -10, 1100
 9241, -200, 10, 1100
 9270, 200, 10, 1100
 9271, -200, 30, 1100
 9300, 200, 30, 1100
 9301, -200, 50, 1100
 9330, 200, 50, 1100
 9331, -200, 70, 1100
 9360, 200, 70, 1100
 9361, -200, 90, 1100
 9390, 200, 90, 1100
 9391, -200, 110, 1100
 9420, 200, 110, 1100
*NGEN, NSET=NAQUA
 9001, 9030, 1
 9031, 9060, 1
 9061, 9090, 1
 9091, 9120, 1
 9121, 9150, 1
 9151, 9180, 1
 9181, 9210, 1
 9211, 9240, 1
 9241, 9270, 1
 9271, 9300, 1
 9301, 9330, 1
 9331, 9360, 1
 9361, 9390, 1
 9391, 9420, 1
*ELEMENT, TYPE=SFM3D4R, ELSET=AQUAVIZ
 9001, 9001, 9002, 9032, 9031
*ELGEN, ELSET = AQUAVIZ
 9001, 29, 1, 1, 13, 30, 30
*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
 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
*OUTPUT, FIELD, NUMBER INTERVAL = 50
*NODE OUTPUT
U
*END STEP