*HEADING
WHIPPING ANALYSIS: STANDARD
**RESTART,WRITE,overlay
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
*NODE, NSET=NALL
*ELEMENT,TYPE=B31,ELSET=BMS101
*BEAM GENERAL SECTION,SECTION=GENERAL,DENSITY=7850.0,ELSET=BMS101
1.571,19.635,0.0,19.635,39.2709
0.0,1.0,0.0
2.1E11,0.808E11
*SECTION POINTS
*TRANSVERSE SHEAR STIFFNESS
6.727E10,6.727E10
*DAMPING,ALPHA=0.005,BETA=0.005
**
** THIS IS FOR A BEAM WITH CIRCULAR CROSS SECTION, OF 5m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
*ELEMENT,TYPE=MASS,ELSET=M1
*MASS,ELSET=M1
 1.0E4
**************************
*SURFACE, NAME=LOADSURF
BMS101,S1
**
** THIS INCIDENT WAVE PROPERTY INDICATES A
** SOURCE, INITIALLY AT (50,15,-15), WITH A STANDOFF
** (REFERENCE) POINT LOCATED AT (50,3.536,-3.536)
**
*INCIDENT WAVE PROPERTY,NAME=UNDEX,TYPE=S
*INCIDENTWAVE FLUID PROPERTY
**
**
** TO DEFINE HICKS-TYPE BUBBLE MODEL BEHAVIOR, SET A = B = 0,
** ETA ( = WAVE EFFECTS) = 0.0, AND C_d = 2.25
**
*AMPLITUDE,DEFINITION=BUBBLE,NAME=HICKSTYP,TIME=TOTAL
52.1E6, 0.09E-3, 0.00, 0.00, 1045.0E6, 1.30, 1630.0, 100.0, 65.0
1025.0, 1.5E3, 0.0, 0.0, 1.0
9.81E0, 1.01325E5, 0.0, 2.25
0.6,,1.e-12,1.e-12
*ELSET,ELSET=ELOUT
*NSET,NSET=QA_TEST
**
** STEP 1
** ------
*STEP,INC=99999,NLGEOM
SHOCK PHASE: Small increments to capture shock pulse
*DYNAMIC,ALPHA=-0.05,INITIAL=YES,direct
8.0E-5,0.007,
**
** this load definition is for:
** x_standoff = (50,3.536,-3.536)
** x_charge   = (50,15   ,-15)
** ref mag    = 1.0
**
*INCIDENTWAVE,PRESSUREAMPLITUDE=HICKSTYP,PROPERTY=UNDEX
LOADSURF,1.0
*INCIDENTWAVEREFLECTION
**
*OUTPUT, FIELD, FREQUENCY=0
*OUTPUT,HISTORY, FREQUENCY=20
*NODE OUTPUT,NSET=QA_TEST
U,
*END STEP
**
*STEP,INC=9999999,NLGEOM
BUBBLE PHASE: Automatic increments
*DYNAMIC,ALPHA=-0.01,INITIAL=YES,HAFTOL=2.e9
1.0E-4,4.993,2.e-12
**
** this load definition is for:
** x_standoff = (50,3.536,-3.536)
** x_charge   = (50,15   ,-15)
** ref mag    = 1.0
**
*INCIDENTWAVE,PRESSUREAMPLITUDE=HICKSTYP,PROPERTY=UNDEX
LOADSURF,1.0
*INCIDENTWAVEREFLECTION
**
*OUTPUT,HISTORY, FREQUENCY=0
*OUTPUT,FIELD, FREQUENCY=0
*END STEP