*HEADING
WHIPPING ANALYSIS: STANDARD
**RESTART,WRITE,overlay
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
*NODE, NSET=NALL
  1,  0.0
  2,  1.0
  3,  2.0
  4,  3.0
  5,  4.0
  6,  5.0
  7,  6.0
  8,  7.0
  9,  8.0
 10,  9.0
 11, 10.0
 12, 11.0
 13, 12.0
 14, 13.0
 15, 14.0
 16, 15.0
 17, 16.0
 18, 17.0
 19, 18.0
 20, 19.0
 21, 20.0
 22, 21.0
 23, 22.0
 24, 23.0
 25, 24.0
 26, 25.0
 27, 26.0
 28, 27.0
 29, 28.0
 30, 29.0
 31, 30.0
 32, 31.0
 33, 32.0
 34, 33.0
 35, 34.0
 36, 35.0
 37, 36.0
 38, 37.0
 39, 38.0
 40, 39.0
 41, 40.0
 42, 41.0
 43, 42.0
 44, 43.0
 45, 44.0
 46, 45.0
 47, 46.0
 48, 47.0
 49, 48.0
 50, 49.0
 51, 50.0
 52, 51.0
 53, 52.0
 54, 53.0
 55, 54.0
 56, 55.0
 57, 56.0
 58, 57.0
 59, 58.0
 60, 59.0
 61, 60.0
 62, 61.0
 63, 62.0
 64, 63.0
 65, 64.0
 66, 65.0
 67, 66.0
 68, 67.0
 69, 68.0
 70, 69.0
 71, 70.0
 72, 71.0
 73, 72.0
 74, 73.0
 75, 74.0
 76, 75.0
 77, 76.0
 78, 77.0
 79, 78.0
 80, 79.0
 81, 80.0
 82, 81.0
 83, 82.0
 84, 83.0
 85, 84.0
 86, 85.0
 87, 86.0
 88, 87.0
 89, 88.0
 90, 89.0
 91, 90.0
 92, 91.0
 93, 92.0
 94, 93.0
 95, 94.0
 96, 95.0
 97, 96.0
 98, 97.0
 99, 98.0
100, 99.0
101,100.0
*ELEMENT,TYPE=B31,ELSET=BMS101
  1,  1,  2
  2,  2,  3
  3,  3,  4
  4,  4,  5
  5,  5,  6
  6,  6,  7
  7,  7,  8
  8,  8,  9
  9,  9, 10
 10, 10, 11
 11, 11, 12
 12, 12, 13
 13, 13, 14
 14, 14, 15
 15, 15, 16
 16, 16, 17
 17, 17, 18
 18, 18, 19
 19, 19, 20
 20, 20, 21
 21, 21, 22
 22, 22, 23
 23, 23, 24
 24, 24, 25
 25, 25, 26
 26, 26, 27
 27, 27, 28
 28, 28, 29
 29, 29, 30
 30, 30, 31
 31, 31, 32
 32, 32, 33
 33, 33, 34
 34, 34, 35
 35, 35, 36
 36, 36, 37
 37, 37, 38
 38, 38, 39
 39, 39, 40
 40, 40, 41
 41, 41, 42
 42, 42, 43
 43, 43, 44
 44, 44, 45
 45, 45, 46
 46, 46, 47
 47, 47, 48
 48, 48, 49
 49, 49, 50
 50, 50, 51
 51, 51, 52
 52, 52, 53
 53, 53, 54
 54, 54, 55
 55, 55, 56
 56, 56, 57
 57, 57, 58
 58, 58, 59
 59, 59, 60
 60, 60, 61
 61, 61, 62
 62, 62, 63
 63, 63, 64
 64, 64, 65
 65, 65, 66
 66, 66, 67
 67, 67, 68
 68, 68, 69
 69, 69, 70
 70, 70, 71
 71, 71, 72
 72, 72, 73
 73, 73, 74
 74, 74, 75
 75, 75, 76
 76, 76, 77
 77, 77, 78
 78, 78, 79
 79, 79, 80
 80, 80, 81
 81, 81, 82
 82, 82, 83
 83, 83, 84
 84, 84, 85
 85, 85, 86
 86, 86, 87
 87, 87, 88
 88, 88, 89
 89, 89, 90
 90, 90, 91
 91, 91, 92
 92, 92, 93
 93, 93, 94
 94, 94, 95
 95, 95, 96
 96, 96, 97
 97, 97, 98
 98, 98, 99
 99, 99,100
100,100,101
*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
2001,  1
2002,  2
2003,  3
2004,  4
2005,  5
2006,  6
2007,  7
2008,  8
2009,  9
2010, 10
2011, 11
2012, 12
2013, 13
2014, 14
2015, 15
2016, 16
2017, 17
2018, 18
2019, 19
2020, 20
2021, 21
2022, 22
2023, 23
2024, 24
2025, 25
2026, 26
2027, 27
2028, 28
2029, 29
2030, 30
2031, 31
2032, 32
2033, 33
2034, 34
2035, 35
2036, 36
2037, 37
2038, 38
2039, 39
2040, 40
2041, 41
2042, 42
2043, 43
2044, 44
2045, 45
2046, 46
2047, 47
2048, 48
2049, 49
2050, 50
2051, 51
2052, 52
2053, 53
2054, 54
2055, 55
2056, 56
2057, 57
2058, 58
2059, 59
2060, 60
2061, 61
2062, 62
2063, 63
2064, 64
2065, 65
2066, 66
2067, 67
2068, 68
2069, 69
2070, 70
2071, 71
2072, 72
2073, 73
2074, 74
2075, 75
2076, 76
2077, 77
2078, 78
2079, 79
2080, 80
2081, 81
2082, 82
2083, 83
2084, 84
2085, 85
2086, 86
2087, 87
2088, 88
2089, 89
2090, 90
2091, 91
2092, 92
2093, 93
2094, 94
2095, 95
2096, 96
2097, 97
2098, 98
2099, 99
2100,100
2101,101
*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
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
0.6,,1.e-12,1.e-12
*ELSET,ELSET=ELOUT
51,
*NSET,NSET=QA_TEST
51,
**
** 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
65.00, 0.0,0.0,1.0, -1.0
**
*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.e10
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
65.00, 0.0,0.0,1.0, -1.0
**
*OUTPUT,HISTORY, FREQUENCY=0
*OUTPUT,FIELD, FREQUENCY=0
*END STEP