*HEADING
: DYNAMIC PRESSURE AND CLOSED-END BUOYANCY [PIPE21]
AIRY WAVES AND *DYNAMIC, EXPLICIT 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.,1000.
  51,0.,1175.
*NODE,NSET=N_TSB
 500,0.,1000.
 501,0.,1001.
*NGEN,NSET=N_PB
 1,51,1
*NSET, NSET=N_PB_END
 1
*ELEMENT,TYPE=PIPE21
 1,1,2
*ELGEN,ELSET=E_PB
 1,50,1,1
*BOUNDARY
 N_PB,1,6
*RESTART,WRITE
*BEAM SECTION,SECTION=PIPE,MATERIAL=MAT1,ELSET=E_PB
 1.0,0.25
 0.,-1.,0.
*ELEMENT,TYPE=B21,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, 
*AMPLITUDE, NAME=GRADUAL, DEFINITION=SMOOTH STEP
 0.0, 0.0, 1.0, 1.0, 15.0, 1.0 
** -------------------------------------------
**     Dynamic pressure test.
** -------------------------------------------
*STEP
 Dynamic pressure test.
*DYNAMIC, EXPLICIT
 ,15
*DLOAD, AMPLITUDE=GRADUAL
 E_PB,PB, 1.0  , 1.1283792
**el,PB, amag0, Diameter (2/sqrt(pi) = area of 1)
*CLOAD, AMPLITUDE=GRADUAL
 500,  TSB,  1.0, 1.0, 0.,-1.
**node,TSB,amag0,area,  normal
**
*OUTPUT, FIELD
*NODE OUTPUT
 RF
**
*OUTPUT, HISTORY, VARIABLE=PRE
**
*OUTPUT, HISTORY
*NODE OUTPUT,NSET=N_PB_END
 RF, 
*NODE OUTPUT,NSET=N_TSB
 RF, 
*END STEP