*HEADING
WHIPPING ANALYSIS: EXPLICIT
*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
10001,50.0,15.0,-15.0
10002,50.0,3.536,-3.536
*nset,nset=source
10001
*nset,nset=standoff
10002
*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
**
*Incident Wave Interaction Property, name=UNDEX, type=SPHERE
1500., 1025.
*UNDEX Charge Property
52.1E6, 0.09E-3, 0.00, 0.00, 1045.0E6, 1.30
1630.0, 100.0, 9.81E0, 1.01325E5, 0.0, 2.25
0.6,,0.0,0.0
65.0, 0.0, 0.0, 1.0
**
*IMPEDANCEPROPERTY,NAME=HARDREFL
0.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,
*amplitude,defi=tab,name=pointy,time=total
-1,0,   0,0,   0.1,1, 0.2,0,
**
*ELSET,ELSET=ELOUT
51,
*NSET,NSET=NOUT
51,
**
** STEP 1
** ------
*STEP
*DYNAMIC,EXPLICIT,direct,initial=yes
1.e-9,5.e-8,
*DIAGNOSTICS, DEFORMATION SPEED CHECK=DETAIL
**
** 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
*Incident Wave Interaction, UNDEX, property=UNDEX
LOADSURF,source,standoff,1.
*INCIDENTWAVEREFLECTION
65.00, 0.0,0.0,1.0, -1.0
**
*OUTPUT, FIELD, NUM=0
*OUTPUT,HISTORY, FREQUENCY=10
*NODE OUTPUT,NSET=NOUT
V2,V3
*ELEMENT OUTPUT,ELSET=ELOUT
1,2,3,4,5
E11
*END STEP
**
** STEP 2
** ------
*STEP
*DYNAMIC,EXPLICIT
,4.95,
**
** 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
*Incident Wave Interaction, UNDEX, property=UNDEX
LOADSURF,source,standoff,1.
*INCIDENTWAVEREFLECTION
65.00, 0.0,0.0,1.0, -1.0
**
*OUTPUT, FIELD, NUM=0
*OUTPUT,HISTORY, FREQUENCY=3
*NODE OUTPUT,NSET=NOUT
V2,V3
*ELEMENT OUTPUT,ELSET=ELOUT
1,2,3,4,5
E11,
*END STEP