*HEADING
WHIPPING ANALYSIS
*RESTART,WRITE,FREQUENCY=20000000
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
****************************************
*Part,name=radius03
*NODE, NSET=NALL
 50, 49.0
 51, 50.0
*ELEMENT,TYPE=B31,ELSET=BMS101
 50, 50, 51
*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
0.0,0.0,0.0,5.0,5.0,0.0,0.0,-5.0,
-5.0,0.0
*TRANSVERSE SHEAR STIFFNESS
6.727E10,6.727E10
*DAMPING,ALPHA=0.005,BETA=0.005
**
** THIS IS FOR A BEAM WITH CIRCULAR CROSS SECTION, OF 0.3m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
1025.0,0.0,0.0,0.3,1.0,0.0
*ELEMENT,TYPE=MASS,ELSET=M1
2050, 50
2051, 51
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius1
*NODE, NSET=NALL
 50, 49.0
 51, 50.0
*ELEMENT,TYPE=B31,ELSET=BMS101
 50, 50, 51
*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
0.0,0.0,0.0,5.0,5.0,0.0,0.0,-5.0,
-5.0,0.0
*TRANSVERSE SHEAR STIFFNESS
6.727E10,6.727E10
*DAMPING,ALPHA=0.005,BETA=0.005
**
** THIS IS FOR A BEAM WITH CIRCULAR CROSS SECTION, OF 1m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
1025.0, 0.0,0.0,1.0,1.0,0.0
*ELEMENT,TYPE=MASS,ELSET=M1
2050, 50
2051, 51
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius3
*NODE, NSET=NALL
 50, 49.0
 51, 50.0
*ELEMENT,TYPE=B31,ELSET=BMS101
 50, 50, 51
*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
0.0,0.0,0.0,5.0,5.0,0.0,0.0,-5.0,
-5.0,0.0
*TRANSVERSE SHEAR STIFFNESS
6.727E10,6.727E10
*DAMPING,ALPHA=0.005,BETA=0.005
**
** THIS IS FOR A BEAM WITH CIRCULAR CROSS SECTION, OF 3m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
1025.0, 0.0,0.0,3.0,1.0,0.0
*ELEMENT,TYPE=MASS,ELSET=M1
2050, 50
2051, 51
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius5
*NODE, NSET=NALL
 50, 49.0
 51, 50.0
*ELEMENT,TYPE=B31,ELSET=BMS101
 50, 50, 51
*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
0.0,0.0,0.0,5.0,5.0,0.0,0.0,-5.0,
-5.0,0.0
*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
1025.0, 0.0,0.0,5.0,1.0,0.0
*ELEMENT,TYPE=MASS,ELSET=M1
2050, 50
2051, 51
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius7
*NODE, NSET=NALL
 50, 49.0
 51, 50.0
*ELEMENT,TYPE=B31,ELSET=BMS101
 50, 50, 51
*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
0.0,0.0,0.0,5.0,5.0,0.0,0.0,-5.0,
-5.0,0.0
*TRANSVERSE SHEAR STIFFNESS
6.727E10,6.727E10
*DAMPING,ALPHA=0.005,BETA=0.005
**
** THIS IS FOR A BEAM WITH CIRCULAR CROSS SECTION, OF 7m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
1025.0, 0.0,0.0,7.0,1.0,0.0
*ELEMENT,TYPE=MASS,ELSET=M1
2050, 50
2051, 51
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius9
*NODE, NSET=NALL
 50, 49.0
 51, 50.0
*ELEMENT,TYPE=B31,ELSET=BMS101
 50, 50, 51
*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
0.0,0.0,0.0,5.0,5.0,0.0,0.0,-5.0,
-5.0,0.0
*TRANSVERSE SHEAR STIFFNESS
6.727E10,6.727E10
*DAMPING,ALPHA=0.005,BETA=0.005
**
** THIS IS FOR A BEAM WITH CIRCULAR CROSS SECTION, OF 9m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
1025.0, 0.0,0.0,9.0,1.0,0.0
*ELEMENT,TYPE=MASS,ELSET=M1
2050, 50
2051, 51
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius17
*NODE, NSET=NALL
 50, 49.0
 51, 50.0
*ELEMENT,TYPE=B31,ELSET=BMS101
 50, 50, 51
*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
0.0,0.0,0.0,5.0,5.0,0.0,0.0,-5.0,
-5.0,0.0
*TRANSVERSE SHEAR STIFFNESS
6.727E10,6.727E10
*DAMPING,ALPHA=0.005,BETA=0.005
**
** THIS IS FOR A BEAM WITH CIRCULAR CROSS SECTION, OF 17m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
1025.0, 0.0,0.0,17.0,1.0,0.0
*ELEMENT,TYPE=MASS,ELSET=M1
2050, 50
2051, 51
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius35
*NODE, NSET=NALL
 50, 49.0
 51, 50.0
*ELEMENT,TYPE=B31,ELSET=BMS101
 50, 50, 51
*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
0.0,0.0,0.0,5.0,5.0,0.0,0.0,-5.0,
-5.0,0.0
*TRANSVERSE SHEAR STIFFNESS
6.727E10,6.727E10
*DAMPING,ALPHA=0.005,BETA=0.005
**
** THIS IS FOR A BEAM WITH CIRCULAR CROSS SECTION, OF 35m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
1025.0, 0.0,0.0,35.0,1.0,0.0
*ELEMENT,TYPE=MASS,ELSET=M1
2050, 50
2051, 51
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Assembly,name=beams
*instance,name=radius03,part=radius03
0,-1,-1
*end instance
*instance,name=radius1,part=radius1
0,0,0
*end instance
*instance,name=radius3,part=radius3
0,1,1
*end instance
*instance,name=radius5,part=radius5
0,2,2
*end instance
*instance,name=radius7,part=radius7
0,3,3
*end instance
*instance,name=radius9,part=radius9
0,4,4
*end instance
*instance,name=radius17,part=radius17
0,5,5
*end instance
*instance,name=radius35,part=radius35
0,6,6
*end instance
****************************************
*ELSET,ELSET=ELOUT
radius03.50,
radius1.50,
radius3.50,
radius5.50,
radius7.50,
radius9.50,
radius17.50,
radius35.50,
*NSET,NSET=QA_TEST
radius03.50,
radius1.50,
radius3.50,
radius5.50,
radius7.50,
radius9.50,
radius17.50,
radius35.50,
*SURFACE, NAME=LOADSURF
radius03.BMS101,S1
radius1.BMS101,S1
radius3.BMS101,S1
radius5.BMS101,S1
radius7.BMS101,S1
radius9.BMS101,S1
radius17.BMS101,S1
radius35.BMS101,S1
*end assembly
**************************
**
** 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
50.0,3.536,-3.536
50.0,15.0,-15.0
*INCIDENTWAVE FLUID PROPERTY
2.30625E9,1025.0,0.0
**
** 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
1.8,
** ------
*STEP,INC=99999,NLGEOM
SHOCK PHASE: Small increments to capture shock pulse and its reflection
*DYNAMIC,ALPHA=-0.05,INITIAL=YES,direct
8.0E-5,0.02,
**
** 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
41.0000, 0.0,-1.0,1.0, 1.0
**
*OUTPUT, HISTORY, FREQUENCY=10
*NODE OUTPUT,NSET=QA_TEST
V2,V3
*OUTPUT,FIELD, FREQUENCY=0
*END STEP
**
*STEP,INC=9999999,NLGEOM
BUBBLE PHASE: Automatic increments 
*DYNAMIC,ALPHA=-0.02,INITIAL=YES,HAFTOL=2.e9
1.0E-4,1.98,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
41.0000, 0.0,-1.0,1.0, 1.0
**
*OUTPUT, HISTORY, FREQUENCY=0
*NODE OUTPUT,NSET=QA_TEST
V2,V3
*OUTPUT,FIELD, FREQUENCY=0
*END STEP