*HEADING
WHIPPING ANALYSIS
*RESTART,WRITE,FREQUENCY=20000
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
****************************************
*Part,name=radius03
*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 0.3m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
*ELEMENT,TYPE=MASS,ELSET=M1
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius1
*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 1m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
*ELEMENT,TYPE=MASS,ELSET=M1
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius3
*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 3m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
*ELEMENT,TYPE=MASS,ELSET=M1
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius5
*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
*End Part
****************************************
*Part,name=radius7
*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 7m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
*ELEMENT,TYPE=MASS,ELSET=M1
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius9
*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 9m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
*ELEMENT,TYPE=MASS,ELSET=M1
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius17
*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 17m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
*ELEMENT,TYPE=MASS,ELSET=M1
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Part,name=radius35
*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 35m RADIUS,
** FULLY IMMERSED IN A FLUID WITH DENSITY = 1025.
**
*BEAM FLUID INERTIA, FULL
*ELEMENT,TYPE=MASS,ELSET=M1
*MASS,ELSET=M1
 1.0E4
*End Part
****************************************
*Assembly,name=beams
*instance,name=radius03,part=radius03
*end instance
*instance,name=radius1,part=radius1
*end instance
*instance,name=radius3,part=radius3
*end instance
*instance,name=radius5,part=radius5
*end instance
*instance,name=radius7,part=radius7
*end instance
*instance,name=radius9,part=radius9
*end instance
*instance,name=radius17,part=radius17
*end instance
*instance,name=radius35,part=radius35
*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
*INCIDENTWAVE FLUID PROPERTY
**
*IMPEDANCEPROPERTY,NAME=HARDREFL
**
** 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
*DYNAMIC,ALPHA=0.0,HAF=1.0e8,INITIAL=YES
2.0E-3,0.2,1.0E-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, HARDREFL
**
*OUTPUT, HISTORY, FREQUENCY=10
*NODE OUTPUT,NSET=QA_TEST
V2,V3
*OUTPUT,FIELD, FREQUENCY=0
*END STEP
**
*STEP,INC=99999,NLGEOM
*DYNAMIC,ALPHA=0.0,HAF=1.0e7,INITIAL=YES
2.0E-3,1.9,1.0E-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, HARDREFL
**
*OUTPUT, HISTORY, FREQUENCY=0
*OUTPUT,FIELD, FREQUENCY=0
*END STEP