# Analytical Models for the Response of the Double-Bottom Structure to Underwater Explosion Based on the Wave Motion Theory

Authors: Yingyu Chen; Xiongliang Yao; Wei Xiao
Journal: Shock and Vibration
 (2016)

---

## Abstract

The aim of this paper is to apply the elastic wave motion theory and the classical one-dimensional cavitation theory to analyze the response of a typical double-bottom structure subjected to underwater blast. The section-varying bar theory and the general acoustic impedance are introduced to get the simplified analytical models. The double-bottom structure is idealized by the basic unit of three substructures which include the simple panel, the panel with stiffener (T-shaped), and the panel associated with girder (I-shaped). According to the simplified models, the analytical models for the corresponding substructures are set up. By taking the cavitation effect into account, the process of fluid-structure interaction can be thoroughly understood, as well as the stress wave propagation. Good agreement between the analytical solution and the finite element prediction is achieved. On the other hand, the Taylor predictions for the panel associated with girder (I-shaped) including the effects of cavitation are invalid, indicating a potential field for the analytical method. The validated analytical models are used to determine the sensitivity of structure response to dimensionless geometric parametersα-k, β-k, and γ-k. Based on the dynamic response of the substructures, we establish the approximate analytical models which are able to predict the response of double-bottom structure to underwater explosion.


---

## Body

## 1. Introduction

A considerable amount of literature exists on the feature of underwater explosion load and the dynamic response of ship structure subjected to underwater explosion [1–7]. The dynamic response of double-bottom hull to underwater explosion is of great interest to the defense of military naval vessels. In order to obtain the dynamic response of double-bottom hull, the effects of fluid-structure interaction (FSI) during the loading phase need to be thoroughly understood. The effects of FSI include two aspects: the shock wave propagation in the complex structure and the cavitation effects in the liquid surrounding the concerned structure.First theoretical studies on FSI date back to World War II. Taylor [8] investigated the response of the free-standing air-backed rigid plate subjected to water blast and found the impulse transmitted to the plate could be reduced by reducing the mass of the plate, but the reduction leads to water cavitation at the fluid-structure interface. The development of Taylor’s solution has been done to investigate the response of protective structures [4, 9–12]. When the input wave acts on the fluid-structure interface, the stress wave propagates along the longitudinal direction in structure [13]. By taking longitudinal stress waves into consideration, Jin et al. [14, 15] proposed an analytical model to solve the response of a coated plate to water blast. However, these methods cannot address the cavitation effects.The loads due to cavitation need to be taken into consideration. By treating water as a bilinear elastic medium, Kennard [16] theoretically studied the one-dimensional cavitation phenomena consequent to the shock wave and described the propagating process of the “breaking front” and the “closing front.” Schiffer et al. [4] analytically investigated the one-dimensional response of a spring-supported rigid plate to underwater blast. They focused on the cavitation phenomena in the fluid, including propagation and arrest of breaking fronts and closing fronts. Later on, the structure response and the cavitation process subsequent to underwater blast were experimentally investigated [17–19]. It was found that the breaking front approaching the structure may invert its direction before reaching the fluid-structure interface and the behavior of water cavitation strongly influenced the pressure histories of the structure. Based on linear wave motion theory, Jin et al. [20] investigated a monolithic elastic coating with varying stiffness and thickness with the consideration of FSI effects and cavitation phenomena. The analytical solutions for propagation of breaking fronts and closing front as well as their interactions with the structure under a water blast were found to be in excellent agreement with finite element (FE) predictions.A significant part of the recent literature on blast loading has concentrated on the underwater blast loading of sandwich structure [17, 21–25]. Deshpande and Fleck [21] and Hutchinson and Xue [22] proposed approximate analytical models for the one-dimensional response of sandwich plate subjected to underwater shock load. McMeeking et al. [25] developed an analytical model for wet surface response, including fluid-structure interaction, which can be used for a wide range of core topologies. This model addressed cavitation and incorporated the momentum of reconstituted water attached to the wet face. Yin et al. [11] approached a similar problem by performing FE simulations and developing analytical models. The results showed that the actual incident wave will be enhanced due to the collapse of cavitation bubbles, especially for the soft core.Although considerable efforts have been devoted to understanding the effects of FSI on the 1D response of monolithic plates and sandwich panels, it still remains unclear how the wave propagation and FSI affect the response of a complex structure such as double-bottom structures. Analytical models for study of the double-bottom structures response are rare. When analyzing a tanker grounding accident, the total response of the assembly is obtained through the summation of the responses of all structural members [26, 27]: plates, cruciform, and stiffeners. Based on similar models, a potential alternative simplified calculation model for the total response of the double-bottom can be obtained in two steps. In the first step, we analytically solve the structural response for the individual structural members: the simple panel, the panel with stiffener (T-shaped), and the panel associated with girder (I-shaped). In the second step, the total response of the double-bottom structure can be obtained through the summation of the responses of the three structural members. Based on this, an effective and quick prediction of the response of double-bottom structure can be achieved.In this study, we examined the dynamic response of the free-standing double-bottom structure plates in contact with a liquid on one side to underwater explosion. We proposed analytical models taking into consideration the stress wave propagation in the structure, the cavitation effects in the fluid, and the corresponding effects on the plates. The general acoustic impedance [28, 29] is adopted to model the shape effect in the panel with stiffener (T-shaped) and the panel associated with girder (I-shaped). The theory of cavitation in Kennard’s research [4, 16] is used to model the cavitation effect in the fluid. Based on the wave motion theory [15], the longitudinal stress waves are taken into account to obtain a more precise solution. One important difference between our approach and Jin et al.’s work [14, 20] is that Jin’s theory is usually employed in the q-layer system. We however use it to analyze the double-bottom structure dynamics by introducing the general acoustic impedance and established the relationship between the substructures dynamics and the double-bottom structure dynamics.The outline of the paper is as follows: in Section2, we present the analytical models for substructures; in Section 3, FE models for corresponding substructures are described and the comparisons between the analytical and FE predictions are carried out; in Section 4, by using the validated analytical models, we explore the sensitivity of the substructures response to different parameters; in Section 5, based on the analysis of substructures, we construct the approximate analytical solutions for a double-bottom structure.


## 2. Analytical Models

### 2.1. Model Characteristics

A typical double-bottom structure of a naval ship subjected to underwater explosion is shown in Figure1(a). This structure is slightly idealized when compared to the real structure. The profile of longitudinal stiffener and girder is rectangular and the double-bottom structure’s height is a constant. The stress waves propagating process in the double-bottom structure hull is rather complex. To achieve a simple and effective model, we divide the whole model into three substructures (see Figure 1(b)): the panel (base component of the ship’s hull structure), a panel with longitudinal stiffener (T-shaped), and girders (I-shaped).Figure 1
Sketches of problem geometry for double-bottom structure: (a) typical double-bottom structure to underwater blast; (b) three substructures for the double-bottom structure.
(a)
(b)The profiles for the panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped) are shown in Figure2 and Table 1. The considered variables are dimensions of profile element’s width and thickness of the bottom flange lf,tf, the web dm,tm, and the top flange lb,tb. The profile was characterized with dimensionless geometric parameter lf/tm,dm/tf,lb/tm,dm/lf. The distance between two girders is 1 m.Table 1
Parameter for the substructures (unit: m).
t
f
t
m
t
b
l
f
d
m
l
b
Simple panel
0.01
—
—
—
—
—
Panel with stiffener (T-shaped)
0.01
0.01
—
0.07
0.07
—
Panel with girder (I-shaped)
0.01
0.01
0.01
0.07
0.8
0.07Figure 2
The geometries of the three substructures used in this analysis.
(a)
Panel
(b)
T-shaped plate
(c)
I-shaped plateDue to the varying cross section of the structure, the stress wave will be redistributed since a portion will be reflected and the other will be transmitted at the interface (see Figure3). ρi, ci, and Ai are the density, sound speed, and the area of the section of the medium i, i=1,2, respectively. The medium on the left side of the interface is taken as medium 1 and the other side is medium 2. Therefore, the subscripts of the variables refer to medium 1 and medium 2, for example, ρ1c1A1 and ρ2c2A2. ΔσI,ΔσR, and ΔσT are the disturbance quantity of stress of incident wave, reflected wave, and transmitted wave, respectively.Figure 3
The wave reflection and transmission at the interface.According to the force equilibrium and compatibility at the interface, we can obtain(1)A1ΔσI+ΔσR=A2ΔσT,ΔσIρ1c1-ΔσRρ1c1=ΔσTρ2c2. By denoting the ratio of general acoustic impedance n=(ρ1c1A1)/(ρ2c2A2), (1) can be rewritten as(2)ΔσR=r·ΔσI,r=1-n1+n,ΔσT=s·ΔσI,s=21+n·A1A2,where r is the reflection coefficient and s is the transmission coefficient.In this paper, the general acoustic impedance [28] is introduced to consider the varying section in the structure. Based on this, we can assume that substructures in double-bottom structure subjected to underwater explosion are treated as one-dimensional multilayer system structures to water blast. The analytical method described in [14] is used to analyze the stress wave propagation in this studies. The following assumptions are made: the substructures (panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped)) are treated as multilayer medium with different general acoustic impedance; the pressure wave travels as a longitudinal wave; and the rigid body motion is neglected.Underwater explosions give rise to spherical shockwaves, traveling in water at approximately sonic speed and impinging on structures. At sufficient distance from the point of denotation, spherical shockwaves can be taken as one-dimensional planar wave. Assume that the double-bottom structure is loaded by a planar, exponentially decaying pressure wave with the peak pressurep0 and decay constant θ, traveling towards the structure. Based on the assumptions we listed above, to thoroughly understand the fluid-structure interaction, we proceed to derive the double-bottom structure response for three cases sketched in Figure 4 and described as follows:(1)
A steel panel plate with one medium (mediumsS1) is in contact with water (see Figure 4(a)).(2)
A panel with stiffener (T-shaped) is treated as a structure with two mediums (mediumsS1 and S2) and is in contact with water (see Figure 4(b)).(3)
An panel with girder (I-shaped) is treated as a structure with three mediums (mediumsS1, S2, and S3) and is in contact with water (see Figure 4(c)).Figure 4
Sketches of geometry, reference system, and loading case for (a) the simple panel; (b) the panel with stiffener (T-shaped); (c) the panel with girder (I-shaped).
(a)
(b)
(c)


### 2.2. Fluid-Structure Interaction

According to the reference system shown in Figure4, waves traveling in the positive x direction at an arbitrary time t can be expressed as(3)pinx,t=p0e-t-x/cw/θ,where p0 is the peak pressure and θ is the decay constant. Applying the continuity of the pressure and particle velocity on the wet surface, we can obtain(4)pint+preft=ptrat,pintρwcw-preftρwcw=ptratρscs, where ρw is the water density, cw is the sonic speed of water, and ρs and cs are the density and sonic speed of the structure. pin(t), pref(t), and ptra(t) represent the pressure induced by incident wave, reflected wave, and transmitted wave on the wet surface, respectively.In this paper, we assume that the water cavitation pressure is zero, in line with the assumptions of previous studies in this field [18, 25, 30]. In order to clarify the interaction process, we divide the response process into two stages: (i) the response prior to the onset of first cavitation and (ii) the response subsequent to the first cavitation.


#### 2.2.1. Response of Three Substructures prior to the Onset of First Cavitation

We would like to reemphasize one important difference between our studying and that of Jin et al.’s work [15]. Jin’s theory is usually employed in multilayer system. We developed a method based on Jin’s work to analyze the double-bottom structure dynamics by introducing the general acoustic impedance and established the analytical models between the substructures and the double-bottom structure.(a) Response of the Simple Panel. Consider the problem sketched in Figure 4(a) of a plate in contact with water on the left side. Let us define t2S1=2hS1/cS1 as the time the wave front propagates through a plate with thickness hS1 and back to the wet surface. Assume that the stress waves have been propagated through the steel plates m1 times. Thus, tm1 can be expressed as(5)tm1=m1·t2S1,where m1 can be taken from 0 to infinity. Let t_total be the total calculation time; then, m1 means the value of t_total/t2S1 rounded to the nearest integer. If m1=0, it means that the wave is not transmitted into the steel plate, as shown in Figure 5(a), but it is reflected. If m1≥1, it means that the input wave is transmitted into the steel plate and propagates m1 times in steel plate, as shown in Figure 5(b).Figure 5
The wave propagation process in the panel: (a) the wave not transmitted into the steel plate (b) the wave transmitted into steel plate once or more.
(a)
m
1
=
0
(b)
m
1
>
0Therefore, the reflection pressure on the wet surface can be expressed as(6)pref,m11tm1,t=rws1·pint,form1=0,sws1ss1wrs1wrs1am1-1rs1a·pint-tm1,form1>1,pref1t=∑i=0m1pref,i1ti,t. Here, the subscript ws1 represents the quantity from water to steel plate and s1w indicates the quantity from steel plate to water. Namely, w refers to water, s1 refers to the steel plate for only one medium, and a is for air. rws1 and sws1 are the reflected coefficient and transmitted coefficient from water to steel plate, respectively. The ratio of two different general acoustic impedance mediums is denoted by nws1=(ρwcwAw)/(ρscsAs). Then, we have rws1=(1-nws1)/(1+nws1) and sws1=2(Aw/As1)/(1+nws1).(b) Response of the Panel with Stiffener (T-Shaped). In Figure 4(b), we define t2S1 as the time the input wave propagates once in medium S1 and defined t2S2 as the time the input wave propagates one time in medium S2. hS1 and hS2 are the thickness of medium S1 and medium S2. Then, we have the t2S1=2hS1/cS1, and t2S2=2hS2/cS2. Thus, the specific time tm1,m2 can be expressed as(7)tm1,m2=m1·t2S1+m2·t2S2.The value ofm1 and m2 can be between 0 and infinity. Let t_total be the total calculation time; then, m1 means the value of t_total/t2S1 rounded to the nearest integer and m2 means the value of t_total/t2S2 rounded to the nearest integer. Different cases are shown in Figure 6. If m1=m2=0, it indicates that the input wave is not transmitted into medium S1. If m1≥1 and m2=0, the wave transmits from water into medium S1 and propagates m1 times in medium S1 only. If m1≥1 and m2≥1, it means the wave propagates m1 times in medium S1 and m2 times in medium S2.Figure 6
The wave propagation process in the panel with stiffener (T-shaped): (a)m1=0, m2=0; (b) m1>0, m2=0; (c) m1>0, m2>0.
(a)
(b)
(c)Consider the wave propagation process sketched in Figure6; the reflected pressure wave on the wet surface of the panel with stiffener (T-shaped) can be summarized as (8)pref,m1,m22tm1,m2,t=rws1·pint,form1=m2=0,sws1ss1wrs1wrs1s2m1-1rs1s2·pint-tm1,0,form1≥1,m2=0,∑i=1kfitm1,m2,t,form1≥1,m2≥1,pref2t=∑i=0m1∑j=0m2pref,i,j2ti,j,t, where(9)fktm1,m2,t=Cm2-1k-1Cm1ksws1ss1wss1s2rs2ass2s1krs1wk-1·rs1wrs1s2m1-krs2ars2s1m2-k·pint-tm1,m2,k=minm1,m2,where the superscript k on fk(tm1,m2,t) indicates k times. If k=1, it means that the wave is transmitted once from medium S1 to medium S2 (see Figure 6(c)). The combinatorial number Cm1k=m1!/(m1-k)!k! represents that from m1 times of propagation occurances in mediums S1 the wave transmits k times into medium S2. Cm2-1k-1=(m2-1)!/[(m2-1)-(k-1)]!(k-1)! means that the wave propagating m2 times in medium S2 is divided into k parts.(c) Response of the Panel with Girder (I-Shaped). Consider the problem sketched in Figure 4(c); we assumed that the input wave propagates in media S1, S2, and S3 for m1, m2, and m3 times, respectively. The specific time tm1,m2,m3 can be written as (10)tm1,m2,m3=m1·t2S1+m2·t2S2+m3·t2S3,where m1, m2, and m3 can be taken from 0 to infinity. Let t_total be the total calculation time; then, m1, m2, and m3 mean the value of t_total/t2S1, t_total/t2S2, and t_total/t2S3 rounded to the nearest integer, respectively.We discuss this problem considering four scenarios (see Figure7).Figure 7
The different cases of the input wave propagates in the panel with girder (I-shaped).
(a)
m
1
=
0, m2=0, and m3=0
(b)
m
1
>
0, m2=0, and m3=0
(c)
m
1
>
0, m2>0, and m3=0
(d)
m
1
>
0, m2>0, and m3>0Ifm1=m2=m3=0, it means the input wave does not transmit into medium S1.Ifm1≥1 and m2=m3=0, as shown in Figure 7(b) the wave transmits from water into medium S1 and propagates only in mediums S1 for m1 times. The wave does not transmit into medium S2 or S3.Ifm1≥1, m2≥1, and m3=0, it indicates that the wave is not transmitted into mediums S3. The wave transmission times from medium S1 into medium S2 can be from 1 to k. The number k is the minimum of m1 and m2.Ifm1≥1, m2≥1, and m3≥1, it indicates the input wave propagates through the whole structure. The transmission times from medium S2 into medium S3 can be from 1 to q. The number q is the minimum of m1, m2, and m3.To sum up, the reflected pressure on the wet surface of the panel with girder (I-shaped) at the specific timetm1,m2,m3 can be written as(11)pref,i,j,k3tm1,m2,m3,t=rws1·pint,form1=m2=m3=0,sws1ss1wrs1wrs1s2m1-1rs1s2·pint-tm1,0,0,form1≥1,m2=m3=0,∑i=1kfktm1,m2,0,t,k=minm1,m2,form1≥1,m2≥1,m3=0,∑i=1q∑i=1kfkqtm1,m2,m3,t,q=minm1,m2,m3,form1≥1,m2≥1,m3≥1,pref3t=∑i=0m1∑j=0m2∑k=0m3pref,i,j,k3ti,j,k,t,where(12)fktm1,m2,0,t=Cm2-1k-1Cm1ksws1ss1wss1s2rs2s3ss2s1krs1wk-1rs1wrs1s2m1-krs2s3rs2s1m2-k·pint-tm1,m2,0,fkqtm1,m2,m3,t=Cm3-1q-1Cm2qCm2-1k-1Cm1ksws1ss1wss1s2ss2s1krs1wm1-1rs1s2m1-krs2s1m2-k·ss3s2ss2s3qrs2s3m2-qrs3s2m3-qrs3am3·pint-tm1,m2,m3.The pressure of the wave that transmitsk times from S1 into medium S2 and transmits q times from S2 into medium S3 is denoted as fkq(tm1,m2,m3,t). The combinatorial number Cm2q indicates that from m2 times of propagation occurrences the wave is transmitted q times from S2 into medium S3. Cm3-1q-1 implies that the wave propagating in medium S3 is divided into q parts.As listed above, the plate velocity and pressure for three substructures prior to the onset of first cavitation can be summarized as(13)pwetkt=pint+prefkt,k=1,2,3,(14)vwetkt=pintρwcw-prefktρwcw,k=1,2,3,where the superscript k represents the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively. vwetk(t) is the velocity of the fluid-structure interface.


#### 2.2.2. Response of Three Substructures Subsequent to the First Cavitation

In this section, the analysis is performed with the assumption that the water response can be described by a linear-elastic, reversible pressure versus volumetric strain constitutive relation forp>0. In contrast, we assume that the fluid is unable to sustain any pressure for p<0; then, we obtain(15)p=-ρwcw2εV,forp>0,0,forp≤0.The pressure in the water between the reflected wave front and the wet face is(16)pwkx,t=pint-xcw+prefkt+xcw,k=1,2,3.For-tcw≤x≤0, the pressure given by (16) is valid until the water first cavities. Prior to cavitation, the velocity in the fluid between the reflected wave front and wet surface is presented as(17)vwkt=pint-x/cwρwcw-prefkt+x/cwρwcw,k=1,2,3,where the superscript k represents the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively.The cavitation timetcav(x) at a location x in the water can be obtained by setting the pressure to zero in (16) and solving the implicit Eq. pwk(x,tx)=0. By taking the minimum value of tcav(x), we can obtain the location and time of first cavitation, which are denoted by xc and tc. According to the location and time of first cavitation, we can achieve the pressure and velocity response subsequent to the first cavitation for three cases.Case 1.
Iftc=tcav(0)or  xc=0, it means the first cavitation takes place on the wet surface. A new breaking front forms. In this case, we assume that the pressure on the wet surface is zero and the velocity is given as an average. Therefore, for t≥tc the pressure and velocity are presented as(18)pwetkt=0,vwetkt=∫0tckpwetktdtρskHeqk,where Heqk refer to the corresponding solid plate thickness.Case 2.
Iftc<tcav(0) or xc<0 and xc<xrest<0,tc<trest, it means after the first cavitation t>tc two breaking waves propagate towards opposite directions at a speed in excess of sound speed cw. The breaking front traveling towards the wet surface can be arrested at position xrest and time trest before reaching the fluid-structure interface (see Figure 8).
The conditions for propagation of a breaking front in the water are given by Kennard [16] as(19)pwx=xbf=0;∂vw∂xx=xbf>0;∂p∂xx=xbf<ρwcw∂vw∂xx=xbf,∂p∂xx=xbf≥0propagationinthe+xdirection,∂p∂xx=xbf≤0propagationinthe-xdirection,where xbf is the current position of the breaking front. Once the breaking front traveling in the +x direction has been arrested, if at least one of the conditions fails, a closing front can take place and propagate in the -x direction. Consequently, propagation of this front begins at position xrest and time trest. The pressure and particle velocities at the closing front denoted by pcf and vcf are given by(20)pcf=pcf,-+pcf,+,vcf=vcf,-+vcf,+=-pcf,-+pcf,+ρwcw, where pcf,- and pcf,+ are wave trains approaching the cavitation region and traveling away from the cavitation region, respectively. vcf,- and vcf,+ are particle velocities associated with the wave trains pcf,- and pcf,+, respectively. According to Kennard, we can have the conditions for propagation of a closing front (21)2pcf,->-ρwcwvcavx;ηx,t=∫tbfxt∂vcavx∂xdt,where vcav(x) represents the particle velocity of the cavitated water, evaluated at a position just ahead of the closing front, namely, vcav(xbf)=vw(xbf,tbf). The quantity ηx,t represents a measure of the strain in the cavitation water. The momentum and mass conservation at the closing front in the original coordinates (see Figure 8) are(22)pcf=1-ηρwdxcfdtvcf-vcav,ρwcw2vcf-vcav=ρwcw2η+pcf1-ηdxcfdt.Equations (20) and (22) can be rewritten as(23)dxcfdt=-cwλλ+cw-λη,pcf=1-ηρwcwλ22λ+cw-2λη,where λ is the auxiliary quantity and is expressed as λ=2pcf,-/ρwcw+vcav.
According to Figure8, the solution in Case 2 can be calculated by three steps.
In the first step,t=trest, the closing front pressure pcf is given by the expression prior to the onset of first cavitation. pcf,- and pcf,+ can be written as(24)pcf,-=preft+xcw,pcf,+=pint-xcw.
In the second step,trest<t<trest-xrest/cw, the pcf,- is given by (24). By solving (23) and checking if the conditions in (19) are satisfied, we can obtain wave trains traveling away from the cavitation region pcf,+ and update the current position of closing front xcf. The value of pcf,+ is used to generate a wave train which will act on the plate at tcf-xcf/cw.
In the third step,t≥trest-xrest/cw, the wave trains pcf,+ reaches the wet surface invalidating the expression prior to the onset of first cavitation for the plates velocity. The analysis of the closing front propagation is similar to the process during trest<t<trest-xrest/cw in the second step. pint is given by(25)pint=pcf,+t-trest+xrestcw.
Based on (25), we can compute the plate pressure and velocity expression subsequent to the first cavitation by solving (13) and (14).Figure 8
Schematic illustration of the phenomena of initial cavitation, emergence and propagation of breaking fronts, and development of a closing front.


## 2.1. Model Characteristics

A typical double-bottom structure of a naval ship subjected to underwater explosion is shown in Figure1(a). This structure is slightly idealized when compared to the real structure. The profile of longitudinal stiffener and girder is rectangular and the double-bottom structure’s height is a constant. The stress waves propagating process in the double-bottom structure hull is rather complex. To achieve a simple and effective model, we divide the whole model into three substructures (see Figure 1(b)): the panel (base component of the ship’s hull structure), a panel with longitudinal stiffener (T-shaped), and girders (I-shaped).Figure 1
Sketches of problem geometry for double-bottom structure: (a) typical double-bottom structure to underwater blast; (b) three substructures for the double-bottom structure.
(a)
(b)The profiles for the panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped) are shown in Figure2 and Table 1. The considered variables are dimensions of profile element’s width and thickness of the bottom flange lf,tf, the web dm,tm, and the top flange lb,tb. The profile was characterized with dimensionless geometric parameter lf/tm,dm/tf,lb/tm,dm/lf. The distance between two girders is 1 m.Table 1
Parameter for the substructures (unit: m).
t
f
t
m
t
b
l
f
d
m
l
b
Simple panel
0.01
—
—
—
—
—
Panel with stiffener (T-shaped)
0.01
0.01
—
0.07
0.07
—
Panel with girder (I-shaped)
0.01
0.01
0.01
0.07
0.8
0.07Figure 2
The geometries of the three substructures used in this analysis.
(a)
Panel
(b)
T-shaped plate
(c)
I-shaped plateDue to the varying cross section of the structure, the stress wave will be redistributed since a portion will be reflected and the other will be transmitted at the interface (see Figure3). ρi, ci, and Ai are the density, sound speed, and the area of the section of the medium i, i=1,2, respectively. The medium on the left side of the interface is taken as medium 1 and the other side is medium 2. Therefore, the subscripts of the variables refer to medium 1 and medium 2, for example, ρ1c1A1 and ρ2c2A2. ΔσI,ΔσR, and ΔσT are the disturbance quantity of stress of incident wave, reflected wave, and transmitted wave, respectively.Figure 3
The wave reflection and transmission at the interface.According to the force equilibrium and compatibility at the interface, we can obtain(1)A1ΔσI+ΔσR=A2ΔσT,ΔσIρ1c1-ΔσRρ1c1=ΔσTρ2c2. By denoting the ratio of general acoustic impedance n=(ρ1c1A1)/(ρ2c2A2), (1) can be rewritten as(2)ΔσR=r·ΔσI,r=1-n1+n,ΔσT=s·ΔσI,s=21+n·A1A2,where r is the reflection coefficient and s is the transmission coefficient.In this paper, the general acoustic impedance [28] is introduced to consider the varying section in the structure. Based on this, we can assume that substructures in double-bottom structure subjected to underwater explosion are treated as one-dimensional multilayer system structures to water blast. The analytical method described in [14] is used to analyze the stress wave propagation in this studies. The following assumptions are made: the substructures (panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped)) are treated as multilayer medium with different general acoustic impedance; the pressure wave travels as a longitudinal wave; and the rigid body motion is neglected.Underwater explosions give rise to spherical shockwaves, traveling in water at approximately sonic speed and impinging on structures. At sufficient distance from the point of denotation, spherical shockwaves can be taken as one-dimensional planar wave. Assume that the double-bottom structure is loaded by a planar, exponentially decaying pressure wave with the peak pressurep0 and decay constant θ, traveling towards the structure. Based on the assumptions we listed above, to thoroughly understand the fluid-structure interaction, we proceed to derive the double-bottom structure response for three cases sketched in Figure 4 and described as follows:(1)
A steel panel plate with one medium (mediumsS1) is in contact with water (see Figure 4(a)).(2)
A panel with stiffener (T-shaped) is treated as a structure with two mediums (mediumsS1 and S2) and is in contact with water (see Figure 4(b)).(3)
An panel with girder (I-shaped) is treated as a structure with three mediums (mediumsS1, S2, and S3) and is in contact with water (see Figure 4(c)).Figure 4
Sketches of geometry, reference system, and loading case for (a) the simple panel; (b) the panel with stiffener (T-shaped); (c) the panel with girder (I-shaped).
(a)
(b)
(c)


## 2.2. Fluid-Structure Interaction

According to the reference system shown in Figure4, waves traveling in the positive x direction at an arbitrary time t can be expressed as(3)pinx,t=p0e-t-x/cw/θ,where p0 is the peak pressure and θ is the decay constant. Applying the continuity of the pressure and particle velocity on the wet surface, we can obtain(4)pint+preft=ptrat,pintρwcw-preftρwcw=ptratρscs, where ρw is the water density, cw is the sonic speed of water, and ρs and cs are the density and sonic speed of the structure. pin(t), pref(t), and ptra(t) represent the pressure induced by incident wave, reflected wave, and transmitted wave on the wet surface, respectively.In this paper, we assume that the water cavitation pressure is zero, in line with the assumptions of previous studies in this field [18, 25, 30]. In order to clarify the interaction process, we divide the response process into two stages: (i) the response prior to the onset of first cavitation and (ii) the response subsequent to the first cavitation.


### 2.2.1. Response of Three Substructures prior to the Onset of First Cavitation

We would like to reemphasize one important difference between our studying and that of Jin et al.’s work [15]. Jin’s theory is usually employed in multilayer system. We developed a method based on Jin’s work to analyze the double-bottom structure dynamics by introducing the general acoustic impedance and established the analytical models between the substructures and the double-bottom structure.(a) Response of the Simple Panel. Consider the problem sketched in Figure 4(a) of a plate in contact with water on the left side. Let us define t2S1=2hS1/cS1 as the time the wave front propagates through a plate with thickness hS1 and back to the wet surface. Assume that the stress waves have been propagated through the steel plates m1 times. Thus, tm1 can be expressed as(5)tm1=m1·t2S1,where m1 can be taken from 0 to infinity. Let t_total be the total calculation time; then, m1 means the value of t_total/t2S1 rounded to the nearest integer. If m1=0, it means that the wave is not transmitted into the steel plate, as shown in Figure 5(a), but it is reflected. If m1≥1, it means that the input wave is transmitted into the steel plate and propagates m1 times in steel plate, as shown in Figure 5(b).Figure 5
The wave propagation process in the panel: (a) the wave not transmitted into the steel plate (b) the wave transmitted into steel plate once or more.
(a)
m
1
=
0
(b)
m
1
>
0Therefore, the reflection pressure on the wet surface can be expressed as(6)pref,m11tm1,t=rws1·pint,form1=0,sws1ss1wrs1wrs1am1-1rs1a·pint-tm1,form1>1,pref1t=∑i=0m1pref,i1ti,t. Here, the subscript ws1 represents the quantity from water to steel plate and s1w indicates the quantity from steel plate to water. Namely, w refers to water, s1 refers to the steel plate for only one medium, and a is for air. rws1 and sws1 are the reflected coefficient and transmitted coefficient from water to steel plate, respectively. The ratio of two different general acoustic impedance mediums is denoted by nws1=(ρwcwAw)/(ρscsAs). Then, we have rws1=(1-nws1)/(1+nws1) and sws1=2(Aw/As1)/(1+nws1).(b) Response of the Panel with Stiffener (T-Shaped). In Figure 4(b), we define t2S1 as the time the input wave propagates once in medium S1 and defined t2S2 as the time the input wave propagates one time in medium S2. hS1 and hS2 are the thickness of medium S1 and medium S2. Then, we have the t2S1=2hS1/cS1, and t2S2=2hS2/cS2. Thus, the specific time tm1,m2 can be expressed as(7)tm1,m2=m1·t2S1+m2·t2S2.The value ofm1 and m2 can be between 0 and infinity. Let t_total be the total calculation time; then, m1 means the value of t_total/t2S1 rounded to the nearest integer and m2 means the value of t_total/t2S2 rounded to the nearest integer. Different cases are shown in Figure 6. If m1=m2=0, it indicates that the input wave is not transmitted into medium S1. If m1≥1 and m2=0, the wave transmits from water into medium S1 and propagates m1 times in medium S1 only. If m1≥1 and m2≥1, it means the wave propagates m1 times in medium S1 and m2 times in medium S2.Figure 6
The wave propagation process in the panel with stiffener (T-shaped): (a)m1=0, m2=0; (b) m1>0, m2=0; (c) m1>0, m2>0.
(a)
(b)
(c)Consider the wave propagation process sketched in Figure6; the reflected pressure wave on the wet surface of the panel with stiffener (T-shaped) can be summarized as (8)pref,m1,m22tm1,m2,t=rws1·pint,form1=m2=0,sws1ss1wrs1wrs1s2m1-1rs1s2·pint-tm1,0,form1≥1,m2=0,∑i=1kfitm1,m2,t,form1≥1,m2≥1,pref2t=∑i=0m1∑j=0m2pref,i,j2ti,j,t, where(9)fktm1,m2,t=Cm2-1k-1Cm1ksws1ss1wss1s2rs2ass2s1krs1wk-1·rs1wrs1s2m1-krs2ars2s1m2-k·pint-tm1,m2,k=minm1,m2,where the superscript k on fk(tm1,m2,t) indicates k times. If k=1, it means that the wave is transmitted once from medium S1 to medium S2 (see Figure 6(c)). The combinatorial number Cm1k=m1!/(m1-k)!k! represents that from m1 times of propagation occurances in mediums S1 the wave transmits k times into medium S2. Cm2-1k-1=(m2-1)!/[(m2-1)-(k-1)]!(k-1)! means that the wave propagating m2 times in medium S2 is divided into k parts.(c) Response of the Panel with Girder (I-Shaped). Consider the problem sketched in Figure 4(c); we assumed that the input wave propagates in media S1, S2, and S3 for m1, m2, and m3 times, respectively. The specific time tm1,m2,m3 can be written as (10)tm1,m2,m3=m1·t2S1+m2·t2S2+m3·t2S3,where m1, m2, and m3 can be taken from 0 to infinity. Let t_total be the total calculation time; then, m1, m2, and m3 mean the value of t_total/t2S1, t_total/t2S2, and t_total/t2S3 rounded to the nearest integer, respectively.We discuss this problem considering four scenarios (see Figure7).Figure 7
The different cases of the input wave propagates in the panel with girder (I-shaped).
(a)
m
1
=
0, m2=0, and m3=0
(b)
m
1
>
0, m2=0, and m3=0
(c)
m
1
>
0, m2>0, and m3=0
(d)
m
1
>
0, m2>0, and m3>0Ifm1=m2=m3=0, it means the input wave does not transmit into medium S1.Ifm1≥1 and m2=m3=0, as shown in Figure 7(b) the wave transmits from water into medium S1 and propagates only in mediums S1 for m1 times. The wave does not transmit into medium S2 or S3.Ifm1≥1, m2≥1, and m3=0, it indicates that the wave is not transmitted into mediums S3. The wave transmission times from medium S1 into medium S2 can be from 1 to k. The number k is the minimum of m1 and m2.Ifm1≥1, m2≥1, and m3≥1, it indicates the input wave propagates through the whole structure. The transmission times from medium S2 into medium S3 can be from 1 to q. The number q is the minimum of m1, m2, and m3.To sum up, the reflected pressure on the wet surface of the panel with girder (I-shaped) at the specific timetm1,m2,m3 can be written as(11)pref,i,j,k3tm1,m2,m3,t=rws1·pint,form1=m2=m3=0,sws1ss1wrs1wrs1s2m1-1rs1s2·pint-tm1,0,0,form1≥1,m2=m3=0,∑i=1kfktm1,m2,0,t,k=minm1,m2,form1≥1,m2≥1,m3=0,∑i=1q∑i=1kfkqtm1,m2,m3,t,q=minm1,m2,m3,form1≥1,m2≥1,m3≥1,pref3t=∑i=0m1∑j=0m2∑k=0m3pref,i,j,k3ti,j,k,t,where(12)fktm1,m2,0,t=Cm2-1k-1Cm1ksws1ss1wss1s2rs2s3ss2s1krs1wk-1rs1wrs1s2m1-krs2s3rs2s1m2-k·pint-tm1,m2,0,fkqtm1,m2,m3,t=Cm3-1q-1Cm2qCm2-1k-1Cm1ksws1ss1wss1s2ss2s1krs1wm1-1rs1s2m1-krs2s1m2-k·ss3s2ss2s3qrs2s3m2-qrs3s2m3-qrs3am3·pint-tm1,m2,m3.The pressure of the wave that transmitsk times from S1 into medium S2 and transmits q times from S2 into medium S3 is denoted as fkq(tm1,m2,m3,t). The combinatorial number Cm2q indicates that from m2 times of propagation occurrences the wave is transmitted q times from S2 into medium S3. Cm3-1q-1 implies that the wave propagating in medium S3 is divided into q parts.As listed above, the plate velocity and pressure for three substructures prior to the onset of first cavitation can be summarized as(13)pwetkt=pint+prefkt,k=1,2,3,(14)vwetkt=pintρwcw-prefktρwcw,k=1,2,3,where the superscript k represents the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively. vwetk(t) is the velocity of the fluid-structure interface.


### 2.2.2. Response of Three Substructures Subsequent to the First Cavitation

In this section, the analysis is performed with the assumption that the water response can be described by a linear-elastic, reversible pressure versus volumetric strain constitutive relation forp>0. In contrast, we assume that the fluid is unable to sustain any pressure for p<0; then, we obtain(15)p=-ρwcw2εV,forp>0,0,forp≤0.The pressure in the water between the reflected wave front and the wet face is(16)pwkx,t=pint-xcw+prefkt+xcw,k=1,2,3.For-tcw≤x≤0, the pressure given by (16) is valid until the water first cavities. Prior to cavitation, the velocity in the fluid between the reflected wave front and wet surface is presented as(17)vwkt=pint-x/cwρwcw-prefkt+x/cwρwcw,k=1,2,3,where the superscript k represents the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively.The cavitation timetcav(x) at a location x in the water can be obtained by setting the pressure to zero in (16) and solving the implicit Eq. pwk(x,tx)=0. By taking the minimum value of tcav(x), we can obtain the location and time of first cavitation, which are denoted by xc and tc. According to the location and time of first cavitation, we can achieve the pressure and velocity response subsequent to the first cavitation for three cases.Case 1.
Iftc=tcav(0)or  xc=0, it means the first cavitation takes place on the wet surface. A new breaking front forms. In this case, we assume that the pressure on the wet surface is zero and the velocity is given as an average. Therefore, for t≥tc the pressure and velocity are presented as(18)pwetkt=0,vwetkt=∫0tckpwetktdtρskHeqk,where Heqk refer to the corresponding solid plate thickness.Case 2.
Iftc<tcav(0) or xc<0 and xc<xrest<0,tc<trest, it means after the first cavitation t>tc two breaking waves propagate towards opposite directions at a speed in excess of sound speed cw. The breaking front traveling towards the wet surface can be arrested at position xrest and time trest before reaching the fluid-structure interface (see Figure 8).
The conditions for propagation of a breaking front in the water are given by Kennard [16] as(19)pwx=xbf=0;∂vw∂xx=xbf>0;∂p∂xx=xbf<ρwcw∂vw∂xx=xbf,∂p∂xx=xbf≥0propagationinthe+xdirection,∂p∂xx=xbf≤0propagationinthe-xdirection,where xbf is the current position of the breaking front. Once the breaking front traveling in the +x direction has been arrested, if at least one of the conditions fails, a closing front can take place and propagate in the -x direction. Consequently, propagation of this front begins at position xrest and time trest. The pressure and particle velocities at the closing front denoted by pcf and vcf are given by(20)pcf=pcf,-+pcf,+,vcf=vcf,-+vcf,+=-pcf,-+pcf,+ρwcw, where pcf,- and pcf,+ are wave trains approaching the cavitation region and traveling away from the cavitation region, respectively. vcf,- and vcf,+ are particle velocities associated with the wave trains pcf,- and pcf,+, respectively. According to Kennard, we can have the conditions for propagation of a closing front (21)2pcf,->-ρwcwvcavx;ηx,t=∫tbfxt∂vcavx∂xdt,where vcav(x) represents the particle velocity of the cavitated water, evaluated at a position just ahead of the closing front, namely, vcav(xbf)=vw(xbf,tbf). The quantity ηx,t represents a measure of the strain in the cavitation water. The momentum and mass conservation at the closing front in the original coordinates (see Figure 8) are(22)pcf=1-ηρwdxcfdtvcf-vcav,ρwcw2vcf-vcav=ρwcw2η+pcf1-ηdxcfdt.Equations (20) and (22) can be rewritten as(23)dxcfdt=-cwλλ+cw-λη,pcf=1-ηρwcwλ22λ+cw-2λη,where λ is the auxiliary quantity and is expressed as λ=2pcf,-/ρwcw+vcav.
According to Figure8, the solution in Case 2 can be calculated by three steps.
In the first step,t=trest, the closing front pressure pcf is given by the expression prior to the onset of first cavitation. pcf,- and pcf,+ can be written as(24)pcf,-=preft+xcw,pcf,+=pint-xcw.
In the second step,trest<t<trest-xrest/cw, the pcf,- is given by (24). By solving (23) and checking if the conditions in (19) are satisfied, we can obtain wave trains traveling away from the cavitation region pcf,+ and update the current position of closing front xcf. The value of pcf,+ is used to generate a wave train which will act on the plate at tcf-xcf/cw.
In the third step,t≥trest-xrest/cw, the wave trains pcf,+ reaches the wet surface invalidating the expression prior to the onset of first cavitation for the plates velocity. The analysis of the closing front propagation is similar to the process during trest<t<trest-xrest/cw in the second step. pint is given by(25)pint=pcf,+t-trest+xrestcw.
Based on (25), we can compute the plate pressure and velocity expression subsequent to the first cavitation by solving (13) and (14).Figure 8
Schematic illustration of the phenomena of initial cavitation, emergence and propagation of breaking fronts, and development of a closing front.


## 2.2.1. Response of Three Substructures prior to the Onset of First Cavitation

We would like to reemphasize one important difference between our studying and that of Jin et al.’s work [15]. Jin’s theory is usually employed in multilayer system. We developed a method based on Jin’s work to analyze the double-bottom structure dynamics by introducing the general acoustic impedance and established the analytical models between the substructures and the double-bottom structure.(a) Response of the Simple Panel. Consider the problem sketched in Figure 4(a) of a plate in contact with water on the left side. Let us define t2S1=2hS1/cS1 as the time the wave front propagates through a plate with thickness hS1 and back to the wet surface. Assume that the stress waves have been propagated through the steel plates m1 times. Thus, tm1 can be expressed as(5)tm1=m1·t2S1,where m1 can be taken from 0 to infinity. Let t_total be the total calculation time; then, m1 means the value of t_total/t2S1 rounded to the nearest integer. If m1=0, it means that the wave is not transmitted into the steel plate, as shown in Figure 5(a), but it is reflected. If m1≥1, it means that the input wave is transmitted into the steel plate and propagates m1 times in steel plate, as shown in Figure 5(b).Figure 5
The wave propagation process in the panel: (a) the wave not transmitted into the steel plate (b) the wave transmitted into steel plate once or more.
(a)
m
1
=
0
(b)
m
1
>
0Therefore, the reflection pressure on the wet surface can be expressed as(6)pref,m11tm1,t=rws1·pint,form1=0,sws1ss1wrs1wrs1am1-1rs1a·pint-tm1,form1>1,pref1t=∑i=0m1pref,i1ti,t. Here, the subscript ws1 represents the quantity from water to steel plate and s1w indicates the quantity from steel plate to water. Namely, w refers to water, s1 refers to the steel plate for only one medium, and a is for air. rws1 and sws1 are the reflected coefficient and transmitted coefficient from water to steel plate, respectively. The ratio of two different general acoustic impedance mediums is denoted by nws1=(ρwcwAw)/(ρscsAs). Then, we have rws1=(1-nws1)/(1+nws1) and sws1=2(Aw/As1)/(1+nws1).(b) Response of the Panel with Stiffener (T-Shaped). In Figure 4(b), we define t2S1 as the time the input wave propagates once in medium S1 and defined t2S2 as the time the input wave propagates one time in medium S2. hS1 and hS2 are the thickness of medium S1 and medium S2. Then, we have the t2S1=2hS1/cS1, and t2S2=2hS2/cS2. Thus, the specific time tm1,m2 can be expressed as(7)tm1,m2=m1·t2S1+m2·t2S2.The value ofm1 and m2 can be between 0 and infinity. Let t_total be the total calculation time; then, m1 means the value of t_total/t2S1 rounded to the nearest integer and m2 means the value of t_total/t2S2 rounded to the nearest integer. Different cases are shown in Figure 6. If m1=m2=0, it indicates that the input wave is not transmitted into medium S1. If m1≥1 and m2=0, the wave transmits from water into medium S1 and propagates m1 times in medium S1 only. If m1≥1 and m2≥1, it means the wave propagates m1 times in medium S1 and m2 times in medium S2.Figure 6
The wave propagation process in the panel with stiffener (T-shaped): (a)m1=0, m2=0; (b) m1>0, m2=0; (c) m1>0, m2>0.
(a)
(b)
(c)Consider the wave propagation process sketched in Figure6; the reflected pressure wave on the wet surface of the panel with stiffener (T-shaped) can be summarized as (8)pref,m1,m22tm1,m2,t=rws1·pint,form1=m2=0,sws1ss1wrs1wrs1s2m1-1rs1s2·pint-tm1,0,form1≥1,m2=0,∑i=1kfitm1,m2,t,form1≥1,m2≥1,pref2t=∑i=0m1∑j=0m2pref,i,j2ti,j,t, where(9)fktm1,m2,t=Cm2-1k-1Cm1ksws1ss1wss1s2rs2ass2s1krs1wk-1·rs1wrs1s2m1-krs2ars2s1m2-k·pint-tm1,m2,k=minm1,m2,where the superscript k on fk(tm1,m2,t) indicates k times. If k=1, it means that the wave is transmitted once from medium S1 to medium S2 (see Figure 6(c)). The combinatorial number Cm1k=m1!/(m1-k)!k! represents that from m1 times of propagation occurances in mediums S1 the wave transmits k times into medium S2. Cm2-1k-1=(m2-1)!/[(m2-1)-(k-1)]!(k-1)! means that the wave propagating m2 times in medium S2 is divided into k parts.(c) Response of the Panel with Girder (I-Shaped). Consider the problem sketched in Figure 4(c); we assumed that the input wave propagates in media S1, S2, and S3 for m1, m2, and m3 times, respectively. The specific time tm1,m2,m3 can be written as (10)tm1,m2,m3=m1·t2S1+m2·t2S2+m3·t2S3,where m1, m2, and m3 can be taken from 0 to infinity. Let t_total be the total calculation time; then, m1, m2, and m3 mean the value of t_total/t2S1, t_total/t2S2, and t_total/t2S3 rounded to the nearest integer, respectively.We discuss this problem considering four scenarios (see Figure7).Figure 7
The different cases of the input wave propagates in the panel with girder (I-shaped).
(a)
m
1
=
0, m2=0, and m3=0
(b)
m
1
>
0, m2=0, and m3=0
(c)
m
1
>
0, m2>0, and m3=0
(d)
m
1
>
0, m2>0, and m3>0Ifm1=m2=m3=0, it means the input wave does not transmit into medium S1.Ifm1≥1 and m2=m3=0, as shown in Figure 7(b) the wave transmits from water into medium S1 and propagates only in mediums S1 for m1 times. The wave does not transmit into medium S2 or S3.Ifm1≥1, m2≥1, and m3=0, it indicates that the wave is not transmitted into mediums S3. The wave transmission times from medium S1 into medium S2 can be from 1 to k. The number k is the minimum of m1 and m2.Ifm1≥1, m2≥1, and m3≥1, it indicates the input wave propagates through the whole structure. The transmission times from medium S2 into medium S3 can be from 1 to q. The number q is the minimum of m1, m2, and m3.To sum up, the reflected pressure on the wet surface of the panel with girder (I-shaped) at the specific timetm1,m2,m3 can be written as(11)pref,i,j,k3tm1,m2,m3,t=rws1·pint,form1=m2=m3=0,sws1ss1wrs1wrs1s2m1-1rs1s2·pint-tm1,0,0,form1≥1,m2=m3=0,∑i=1kfktm1,m2,0,t,k=minm1,m2,form1≥1,m2≥1,m3=0,∑i=1q∑i=1kfkqtm1,m2,m3,t,q=minm1,m2,m3,form1≥1,m2≥1,m3≥1,pref3t=∑i=0m1∑j=0m2∑k=0m3pref,i,j,k3ti,j,k,t,where(12)fktm1,m2,0,t=Cm2-1k-1Cm1ksws1ss1wss1s2rs2s3ss2s1krs1wk-1rs1wrs1s2m1-krs2s3rs2s1m2-k·pint-tm1,m2,0,fkqtm1,m2,m3,t=Cm3-1q-1Cm2qCm2-1k-1Cm1ksws1ss1wss1s2ss2s1krs1wm1-1rs1s2m1-krs2s1m2-k·ss3s2ss2s3qrs2s3m2-qrs3s2m3-qrs3am3·pint-tm1,m2,m3.The pressure of the wave that transmitsk times from S1 into medium S2 and transmits q times from S2 into medium S3 is denoted as fkq(tm1,m2,m3,t). The combinatorial number Cm2q indicates that from m2 times of propagation occurrences the wave is transmitted q times from S2 into medium S3. Cm3-1q-1 implies that the wave propagating in medium S3 is divided into q parts.As listed above, the plate velocity and pressure for three substructures prior to the onset of first cavitation can be summarized as(13)pwetkt=pint+prefkt,k=1,2,3,(14)vwetkt=pintρwcw-prefktρwcw,k=1,2,3,where the superscript k represents the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively. vwetk(t) is the velocity of the fluid-structure interface.


## 2.2.2. Response of Three Substructures Subsequent to the First Cavitation

In this section, the analysis is performed with the assumption that the water response can be described by a linear-elastic, reversible pressure versus volumetric strain constitutive relation forp>0. In contrast, we assume that the fluid is unable to sustain any pressure for p<0; then, we obtain(15)p=-ρwcw2εV,forp>0,0,forp≤0.The pressure in the water between the reflected wave front and the wet face is(16)pwkx,t=pint-xcw+prefkt+xcw,k=1,2,3.For-tcw≤x≤0, the pressure given by (16) is valid until the water first cavities. Prior to cavitation, the velocity in the fluid between the reflected wave front and wet surface is presented as(17)vwkt=pint-x/cwρwcw-prefkt+x/cwρwcw,k=1,2,3,where the superscript k represents the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively.The cavitation timetcav(x) at a location x in the water can be obtained by setting the pressure to zero in (16) and solving the implicit Eq. pwk(x,tx)=0. By taking the minimum value of tcav(x), we can obtain the location and time of first cavitation, which are denoted by xc and tc. According to the location and time of first cavitation, we can achieve the pressure and velocity response subsequent to the first cavitation for three cases.Case 1.
Iftc=tcav(0)or  xc=0, it means the first cavitation takes place on the wet surface. A new breaking front forms. In this case, we assume that the pressure on the wet surface is zero and the velocity is given as an average. Therefore, for t≥tc the pressure and velocity are presented as(18)pwetkt=0,vwetkt=∫0tckpwetktdtρskHeqk,where Heqk refer to the corresponding solid plate thickness.Case 2.
Iftc<tcav(0) or xc<0 and xc<xrest<0,tc<trest, it means after the first cavitation t>tc two breaking waves propagate towards opposite directions at a speed in excess of sound speed cw. The breaking front traveling towards the wet surface can be arrested at position xrest and time trest before reaching the fluid-structure interface (see Figure 8).
The conditions for propagation of a breaking front in the water are given by Kennard [16] as(19)pwx=xbf=0;∂vw∂xx=xbf>0;∂p∂xx=xbf<ρwcw∂vw∂xx=xbf,∂p∂xx=xbf≥0propagationinthe+xdirection,∂p∂xx=xbf≤0propagationinthe-xdirection,where xbf is the current position of the breaking front. Once the breaking front traveling in the +x direction has been arrested, if at least one of the conditions fails, a closing front can take place and propagate in the -x direction. Consequently, propagation of this front begins at position xrest and time trest. The pressure and particle velocities at the closing front denoted by pcf and vcf are given by(20)pcf=pcf,-+pcf,+,vcf=vcf,-+vcf,+=-pcf,-+pcf,+ρwcw, where pcf,- and pcf,+ are wave trains approaching the cavitation region and traveling away from the cavitation region, respectively. vcf,- and vcf,+ are particle velocities associated with the wave trains pcf,- and pcf,+, respectively. According to Kennard, we can have the conditions for propagation of a closing front (21)2pcf,->-ρwcwvcavx;ηx,t=∫tbfxt∂vcavx∂xdt,where vcav(x) represents the particle velocity of the cavitated water, evaluated at a position just ahead of the closing front, namely, vcav(xbf)=vw(xbf,tbf). The quantity ηx,t represents a measure of the strain in the cavitation water. The momentum and mass conservation at the closing front in the original coordinates (see Figure 8) are(22)pcf=1-ηρwdxcfdtvcf-vcav,ρwcw2vcf-vcav=ρwcw2η+pcf1-ηdxcfdt.Equations (20) and (22) can be rewritten as(23)dxcfdt=-cwλλ+cw-λη,pcf=1-ηρwcwλ22λ+cw-2λη,where λ is the auxiliary quantity and is expressed as λ=2pcf,-/ρwcw+vcav.
According to Figure8, the solution in Case 2 can be calculated by three steps.
In the first step,t=trest, the closing front pressure pcf is given by the expression prior to the onset of first cavitation. pcf,- and pcf,+ can be written as(24)pcf,-=preft+xcw,pcf,+=pint-xcw.
In the second step,trest<t<trest-xrest/cw, the pcf,- is given by (24). By solving (23) and checking if the conditions in (19) are satisfied, we can obtain wave trains traveling away from the cavitation region pcf,+ and update the current position of closing front xcf. The value of pcf,+ is used to generate a wave train which will act on the plate at tcf-xcf/cw.
In the third step,t≥trest-xrest/cw, the wave trains pcf,+ reaches the wet surface invalidating the expression prior to the onset of first cavitation for the plates velocity. The analysis of the closing front propagation is similar to the process during trest<t<trest-xrest/cw in the second step. pint is given by(25)pint=pcf,+t-trest+xrestcw.
Based on (25), we can compute the plate pressure and velocity expression subsequent to the first cavitation by solving (13) and (14).Figure 8
Schematic illustration of the phenomena of initial cavitation, emergence and propagation of breaking fronts, and development of a closing front.


## 3. Validation

### 3.1. Finite Element Model

As shown in Figure9, the three-dimensional numerical model consists of a water column on the left side of the wet surface. Contact is enforced at the fluid-structure interface, and the free-standing plate is modeled by applying symmetry boundary conditions. The commercial code ABAQUS [31] is used to validate the analytical models described in Section 2.2.Figure 9
A sketch of the boundary conditions in the analysis of free-standing panel with girder (I-shaped) impinged by a water blast.Three-dimensional FE simulations were performed to provide more insight into fluid and structural response. For the corresponding substructures, we consider a homogeneous material withρ=7800Kg/m3, Young’s modulus E=2.11e11Pa, and Poisson’s ratio μ=0.3. The dimensions of the substructures are shown in Table 1. The eight-node three-dimensional brick elements with reduced integration (C3D8R) are used to model the structure. The element meshes of substructures are shown in Figure 10. To capture the stress wave in the structure, the element sizes les1_l, les2_l, and les3_l along the longitudinal direction of the substructures are set to 0.1 mm. The element sizes les1_t, les2_t, and les3_t of 10 mm are used to discretize the substructures along the width, while a single element is used to discretize the substructures along the depth. The elements of the panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped) are 700, 1400, and 9400, respectively.Figure 10
The mesh of fluid and structure: (a) the water column; (b) the panel with stiffener (T-shaped); (c) the panel with girder (I-shaped).
(a)
(b)
(c)
(d)In ABAQUS, by using a Mie-Gruneisen equation of state with a linear Hugoniot relation, we can model the water with the linear-elastic relation described in (15). The density and sound speed of water are set to ρw=1000Kg/m3 and cw=1500m/s, respectively. The bulk modulus of water is Kw=ρwcw2=2.25GPa. The shear modulus and hydrostatic tensile strength of water are set to be zero. Acoustic eight-node three-dimensional brick elements with reduced integration (AC3D8R) are employed to model the water column. To capture the cavitation effect in the water column, the length of the water column is set to 1.2 m. The element size lew_l of the water along the longitudinal direction of the water column is set to 0.3 mm, while the element size lew_t of 10 mm is used to discretize the water along the width, while one single element is used to discretize the water along the depth. The water column comprises 28000 elements.The water cavitation effects are realized by total wave formulation. We set the cavitation limit to be zero and the hydrostatic pressure field in the water column to be sufficiently small. A uniform pressure is imposed on the left surface of the water column. The impulses used all have peak pressurep0=10MPa and the decay constant θ=0.1ms.


### 3.2. Comparison of Analytical and FE Prediction

In this section, we compare the analytical and FE predictions for fluid-structure interface pressure and velocity histories, as well as the motion of breaking and closing fronts, to assess the accuracy of the analytical models.For simplicity, we define nondimensional parameters:t-=t/θ, x-=x/cwθ, p-=p/p0, v-=p/p0/ρwcw, and I-=I/p0θ, where t-, x-, p-, v-, and I- are the nondimensional time, displacement, pressure, velocity, and impulse, respectively.


#### 3.2.1. Closing Front and Breaking Front

Figure11 presents the cavitation time t-c as a function of the location x- and the velocity v-cav in the cavitated water at t-cx for the three substructures. The minimum value on these curves (circled in Figure 11(a)) is the location and time of the onset of first cavitation. Two branches of each curve originated from this point, corresponding to propagation of two breaking fronts traveling away from (left of x-c) and towards (right of x-c) the wet surface. Note that the cavitation time is linear with x- on the left of x-c, except very close to the first cavitation location.Figure 11
Analytical predictions and FE prediction of the spatial variation of (a) the time to cavitation and (b) the velocityv-cav in the cavitated water at t-cx. These results are for cases in Table 1, which are the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped).
(a)
(b)For the simple panel, cavitation initiates at a very small distance from the wet surface and gives rise to a breaking front traveling away from the interface, while the breaking front approaching the interface cannot be arrested before reaching the wet surface. For the panel with stiffener (T-shaped) and the panel with girder (I-shaped), the first cavitation is located within the fluid and the breaking front traveling towards the structure would reverse its motion and become the closing front. Figure11(a) shows the first cavitation time for the water column. For example, the jump in the curve of the panel with girder (I-shaped) indicates that when the first cavitation occurs at x¯=-0.298, the breaking front propagated away from the wet surface. As time elapses, another cavitation takes place at the interface at time t-c0. And this breaking front also propagated away from the wet surface. These two breaking fronts generate from different place, while they travel at the same velocity. The jump means that the cavitation will occur twice at the same displacement and different time. We can see a similar phenomenon in the case of the panel with stiffener (T-shaped).As shown in Figure11(a), the trends of two breaking fronts in the analytical results agree well with FE results. The analytical results in Figure 11(a) only show the first cavitation time for the water column. Thus, when the cavitation occurs at the same displacement and different time, the curve will “jump” without showing the closing front. And this leads to the discrepancy between FE and analytical predictions. From Figure 11(a), we can see that t-c0<trest/θ-xrest/cwθ indicates that the cavitation occurs on the wet surface before the wave train emanating from the closing front reaches the plate.The expressions (13) for the plate velocity and pressure are still valid when t-<t-c0. In Figure 11(b), the velocity v-cav in the cavitated water at t-cx for the three substructures is shown and small difference is found between the analytical solutions and the FE predictions. The particle velocity v-cav0 is the corresponding velocity on the wet surface when cavitation occurs. The value of v-cav0 for the panel, the panel with the stiffener (T-shaped), and the panel with girder (I-shaped) are 1.568, 1.407, and 0.6657, respectively. The variance from the analytical results for the different shape plates is large. This is due to the fact that the panel with girder (I-shaped) has a larger unit mass.


#### 3.2.2. Pressure and Velocity Histories on Wet Surface

Figure12(a) shows nondimensional time histories of nondimensional wet surface pressure p-wet1, while the corresponding time histories of nondimensional velocity v-wet1 are presented in Figure 12(b). As discussed above, in this particular case the first cavitation takes place at a location very close to the wet surface. We can then calculate the analytical solution of the pressure and velocity by using (18) after the wet surface is cavitated. In particular, we investigated the dynamic response of a steel plate by a three-method analytical solution, FEM results, and Taylor’s solution and then did a comparison of the predictions. Both Figures 12(a) and 12(b) show that the analytical solutions are adequately predicted by FE results and Taylor’s solutions [8]. On the whole process, the analytical solution agrees well with the FE prediction for the pressure and velocity. It is to be noted that the small jumps of analytical solutions indicate the reflections and transmissions of the pressure wave. When the stress wave propagation and the cavitation effect are considered, the analytical results are more accurate than Taylor’s results and FE predictions.Figure 12
Nondimensional time histories of (a) nondimensional wet surface pressure and (b) nondimensional velocity for the simple panel.
(a)
(b)To reveal the fluid-structure interaction of the double-bottom structure, the nondimensional time histories of wet surface pressure and velocity for the panel with stiffener (T-shaped) and the panel with girder (I-shaped) are shown in Figures13 and 14, respectively. The FE predictions are presented to validate the analytical solutions. As shown, the analytical solutions agree considerably well with the FE predictions. The oscillations in FE predictions can be explained by the fact that on one hand the FE algorithm for the thin steel plate accounts for the FE predictions that have large numerical diffusions; on the other hand, the FE numerical models are three-dimensional, including the effect of the transverse wave. We note that if cavitation occurs on the wet surface, the velocity of interface would be given as an average. After the wet surface gets cavitated, the process becomes very complex to solve analytically. Thus, the velocity is calculated as an average. Hence, we can conclude that the analytical solutions are validated in these cases.Figure 13
Nondimensional time histories of (a) nondimensional wet surface pressure and (b) nondimensional velocity for the panel with stiffener (T-shaped).
(a)
(b)Figure 14
Nondimensional time histories of (a) nondimensional wet surface pressure and (b) nondimensional velocity for the panel with girder (I-shaped).
(a)
(b)Note that the obvious drops in the analytical solutions indicate the superposition of the reflected wave from the rear surface of the structure and the incident wave. For the panel with stiffener (T-shaped), this reflected wave is the rarefaction wave which causes the great drop in the pressure and the decrement of decay constant (see Figure13(a)). The specific time ticks circled in Figure 13(a), t=0.05, t=0.105, and t=0.1583, are the time ticks of stress wave traveling once, twice, and thrice through the panel with stiffener (T-shaped), namely, the time it takes for the reflected wave from the rear surface of the structure to reach the wet surface. Between two arbitrary time ticks, small oscillations exist. This is as a result of the superposition of the incident wave and the reflected wave from the interface of medium S1 and medium S2. As such, we can characterize the response of the panel with stiffener (T-shaped) by the drops in specific time tick and the small oscillations between the time ticks. The response clearly shows the wave propagation process in T-shaped plate. In Figure 13(b), the characteristics induced by superposition of multiwaves can be expressed as the jumps at specific time ticks and small oscillations between time ticks. Circles in Figure 13(b) are the maximum value of the velocity and the velocity when the wet surface gets cavitated, that is, 1.422 and 1.404. The panel with stiffener (T-shaped) gets its peak velocity before cavitation occurs on the wet surface. As shown, the analytical solution is a little larger than the FE prediction and the differences are reasonably acceptable. As shown in Figure 13, Taylor’s solutions for rigid plate equivalent to panel with stiffener (T-shaped) have similar trend as the analytical and FE prediction. However, Taylor’s solution cannot show the stress wave propagation process in the specific structure and it cannot recognize the minor change in the structure.In Figure14(a), we can see the obvious small jumps at the specific time (t=0.53,t=0.1068) which is the result of the enlargement effect when the stress wave is transmitted from medium S2 to medium S3. The details of the enhancement of pressure will be expressed in Section 4.2. The sharp drop of pressure occurs when the rarefaction wave from the rear surface of the structure reaches the wet surface. Similar to the wave propagation process in the panel with stiffener (T-shaped), the corresponding response of the panel with girder (I-shaped) is characterized by big drops at specific time tick and small oscillations between time ticks. The differences are that the wave travel in the panel with girder (I-shaped) takes longer time to finish one trip through the structure and the small jumps take place due to the existence of the medium S3. In Figure 14(b), we also can describe the trends in nondimensional time history of velocity as jumps at specific time ticks and small oscillations between time ticks, as well as the small drops owing to the existence of the medium S3. The panel with girder (I-shaped) gets its peak velocity before cavitation occurs on the wet surface (see circles in Figure 14(b)). By comparing the velocities in Figures 12(b), 13(b), and 14(b), we can conclude that the velocity of the wet surface is considerably small if the structure is heavy or strong. Correspondingly, the rarefaction wave radiated by the motion of wet surface is small. Thus, the pressure on the wet surface decreases slowly (see Figure 14(a)). As shown in Figure 14, Taylor’s solutions for rigid plate equal to the panel with girder (I-shaped) have a similar trend as the analytical and FE predictions. However, there is a big difference between Taylor’s solution and FE predictions. It can be seen that Taylor’s solution cannot capture the reflected wave from the varying cross section of the structure.


#### 3.2.3. Impulse Histories on Wet Surface

The intensity of loads impacted on the structure consequent to underwater explosion, as quantified by the specific impulse (per unit area) applied on the wet surface, is deduced by solving the following equation:(26)Ikt=∫0tpwetkτdτ,k=1,2,3,where k=1,2,3 refers to the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively.Define the maximum incident impulseI0 (integral of the incident wave) per unit quasistatic state as I0=2p0θ. The severity of structural loading consequent to underwater explosion can be quantified by the maximum nondimensional impulse, I-max=Imax/I0, applied on the load surface of the structure.The analytical solutions for the nondimensional impulseI-k of three substructures agree well with the FE predictions and Taylor’s solutions (see Figure 15). As the thickness of the structure increases, the plate takes longer time to reach its maximum nondimensional impulse. This is due to the fact that the larger thickness of the structure will prolong the propagating time and thus release higher energy to the plate. The comparisons of the analytical results and FE predictions of maximum nondimensional impulse I-maxk for the three substructures are presented in Table 2. We can see that the numerical predictions of I-maxk are little smaller than the analytical results. This is probably caused by the numerical dissipation. As shown in the second column of Table 2, the normalized impulses I-maxk for the three types of plates equal to 6.81%, 12.17%, and 39.87%, respectively, and are all lower than unity. Moreover, comparison between the three plates shows that the specific impulse applied on the wet surface of the panel with girder (I-shaped) is larger than that of the panel with stiffener (T-shaped) and the simple panel, indicating that more momentum is transmitted to the panel with girder (I-shaped). This is due to the fact that the pressure on the wet surface of the panel with girder (I-shaped) decreases slowly and therefore the impulse impacting the plate will increase.Table 2
Analytical and FE predictions of normalized impulse.
Case
Analytical results
FE predictions
Error
Simple panel
0.0681
0.0673
1.25%
The panel with stiffener (T-shaped)
0.1217
0.1184
2.71%
The panel with girder (I-shaped)
0.3987
0.3977
0.25%Figure 15
The analytical solutions and FE predictions of the nondimensional impulse for (a) simple panel; (b) the panel with stiffener (T-shaped); (c) the panel with girder (I-shaped).


## 3.1. Finite Element Model

As shown in Figure9, the three-dimensional numerical model consists of a water column on the left side of the wet surface. Contact is enforced at the fluid-structure interface, and the free-standing plate is modeled by applying symmetry boundary conditions. The commercial code ABAQUS [31] is used to validate the analytical models described in Section 2.2.Figure 9
A sketch of the boundary conditions in the analysis of free-standing panel with girder (I-shaped) impinged by a water blast.Three-dimensional FE simulations were performed to provide more insight into fluid and structural response. For the corresponding substructures, we consider a homogeneous material withρ=7800Kg/m3, Young’s modulus E=2.11e11Pa, and Poisson’s ratio μ=0.3. The dimensions of the substructures are shown in Table 1. The eight-node three-dimensional brick elements with reduced integration (C3D8R) are used to model the structure. The element meshes of substructures are shown in Figure 10. To capture the stress wave in the structure, the element sizes les1_l, les2_l, and les3_l along the longitudinal direction of the substructures are set to 0.1 mm. The element sizes les1_t, les2_t, and les3_t of 10 mm are used to discretize the substructures along the width, while a single element is used to discretize the substructures along the depth. The elements of the panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped) are 700, 1400, and 9400, respectively.Figure 10
The mesh of fluid and structure: (a) the water column; (b) the panel with stiffener (T-shaped); (c) the panel with girder (I-shaped).
(a)
(b)
(c)
(d)In ABAQUS, by using a Mie-Gruneisen equation of state with a linear Hugoniot relation, we can model the water with the linear-elastic relation described in (15). The density and sound speed of water are set to ρw=1000Kg/m3 and cw=1500m/s, respectively. The bulk modulus of water is Kw=ρwcw2=2.25GPa. The shear modulus and hydrostatic tensile strength of water are set to be zero. Acoustic eight-node three-dimensional brick elements with reduced integration (AC3D8R) are employed to model the water column. To capture the cavitation effect in the water column, the length of the water column is set to 1.2 m. The element size lew_l of the water along the longitudinal direction of the water column is set to 0.3 mm, while the element size lew_t of 10 mm is used to discretize the water along the width, while one single element is used to discretize the water along the depth. The water column comprises 28000 elements.The water cavitation effects are realized by total wave formulation. We set the cavitation limit to be zero and the hydrostatic pressure field in the water column to be sufficiently small. A uniform pressure is imposed on the left surface of the water column. The impulses used all have peak pressurep0=10MPa and the decay constant θ=0.1ms.


## 3.2. Comparison of Analytical and FE Prediction

In this section, we compare the analytical and FE predictions for fluid-structure interface pressure and velocity histories, as well as the motion of breaking and closing fronts, to assess the accuracy of the analytical models.For simplicity, we define nondimensional parameters:t-=t/θ, x-=x/cwθ, p-=p/p0, v-=p/p0/ρwcw, and I-=I/p0θ, where t-, x-, p-, v-, and I- are the nondimensional time, displacement, pressure, velocity, and impulse, respectively.


### 3.2.1. Closing Front and Breaking Front

Figure11 presents the cavitation time t-c as a function of the location x- and the velocity v-cav in the cavitated water at t-cx for the three substructures. The minimum value on these curves (circled in Figure 11(a)) is the location and time of the onset of first cavitation. Two branches of each curve originated from this point, corresponding to propagation of two breaking fronts traveling away from (left of x-c) and towards (right of x-c) the wet surface. Note that the cavitation time is linear with x- on the left of x-c, except very close to the first cavitation location.Figure 11
Analytical predictions and FE prediction of the spatial variation of (a) the time to cavitation and (b) the velocityv-cav in the cavitated water at t-cx. These results are for cases in Table 1, which are the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped).
(a)
(b)For the simple panel, cavitation initiates at a very small distance from the wet surface and gives rise to a breaking front traveling away from the interface, while the breaking front approaching the interface cannot be arrested before reaching the wet surface. For the panel with stiffener (T-shaped) and the panel with girder (I-shaped), the first cavitation is located within the fluid and the breaking front traveling towards the structure would reverse its motion and become the closing front. Figure11(a) shows the first cavitation time for the water column. For example, the jump in the curve of the panel with girder (I-shaped) indicates that when the first cavitation occurs at x¯=-0.298, the breaking front propagated away from the wet surface. As time elapses, another cavitation takes place at the interface at time t-c0. And this breaking front also propagated away from the wet surface. These two breaking fronts generate from different place, while they travel at the same velocity. The jump means that the cavitation will occur twice at the same displacement and different time. We can see a similar phenomenon in the case of the panel with stiffener (T-shaped).As shown in Figure11(a), the trends of two breaking fronts in the analytical results agree well with FE results. The analytical results in Figure 11(a) only show the first cavitation time for the water column. Thus, when the cavitation occurs at the same displacement and different time, the curve will “jump” without showing the closing front. And this leads to the discrepancy between FE and analytical predictions. From Figure 11(a), we can see that t-c0<trest/θ-xrest/cwθ indicates that the cavitation occurs on the wet surface before the wave train emanating from the closing front reaches the plate.The expressions (13) for the plate velocity and pressure are still valid when t-<t-c0. In Figure 11(b), the velocity v-cav in the cavitated water at t-cx for the three substructures is shown and small difference is found between the analytical solutions and the FE predictions. The particle velocity v-cav0 is the corresponding velocity on the wet surface when cavitation occurs. The value of v-cav0 for the panel, the panel with the stiffener (T-shaped), and the panel with girder (I-shaped) are 1.568, 1.407, and 0.6657, respectively. The variance from the analytical results for the different shape plates is large. This is due to the fact that the panel with girder (I-shaped) has a larger unit mass.


### 3.2.2. Pressure and Velocity Histories on Wet Surface

Figure12(a) shows nondimensional time histories of nondimensional wet surface pressure p-wet1, while the corresponding time histories of nondimensional velocity v-wet1 are presented in Figure 12(b). As discussed above, in this particular case the first cavitation takes place at a location very close to the wet surface. We can then calculate the analytical solution of the pressure and velocity by using (18) after the wet surface is cavitated. In particular, we investigated the dynamic response of a steel plate by a three-method analytical solution, FEM results, and Taylor’s solution and then did a comparison of the predictions. Both Figures 12(a) and 12(b) show that the analytical solutions are adequately predicted by FE results and Taylor’s solutions [8]. On the whole process, the analytical solution agrees well with the FE prediction for the pressure and velocity. It is to be noted that the small jumps of analytical solutions indicate the reflections and transmissions of the pressure wave. When the stress wave propagation and the cavitation effect are considered, the analytical results are more accurate than Taylor’s results and FE predictions.Figure 12
Nondimensional time histories of (a) nondimensional wet surface pressure and (b) nondimensional velocity for the simple panel.
(a)
(b)To reveal the fluid-structure interaction of the double-bottom structure, the nondimensional time histories of wet surface pressure and velocity for the panel with stiffener (T-shaped) and the panel with girder (I-shaped) are shown in Figures13 and 14, respectively. The FE predictions are presented to validate the analytical solutions. As shown, the analytical solutions agree considerably well with the FE predictions. The oscillations in FE predictions can be explained by the fact that on one hand the FE algorithm for the thin steel plate accounts for the FE predictions that have large numerical diffusions; on the other hand, the FE numerical models are three-dimensional, including the effect of the transverse wave. We note that if cavitation occurs on the wet surface, the velocity of interface would be given as an average. After the wet surface gets cavitated, the process becomes very complex to solve analytically. Thus, the velocity is calculated as an average. Hence, we can conclude that the analytical solutions are validated in these cases.Figure 13
Nondimensional time histories of (a) nondimensional wet surface pressure and (b) nondimensional velocity for the panel with stiffener (T-shaped).
(a)
(b)Figure 14
Nondimensional time histories of (a) nondimensional wet surface pressure and (b) nondimensional velocity for the panel with girder (I-shaped).
(a)
(b)Note that the obvious drops in the analytical solutions indicate the superposition of the reflected wave from the rear surface of the structure and the incident wave. For the panel with stiffener (T-shaped), this reflected wave is the rarefaction wave which causes the great drop in the pressure and the decrement of decay constant (see Figure13(a)). The specific time ticks circled in Figure 13(a), t=0.05, t=0.105, and t=0.1583, are the time ticks of stress wave traveling once, twice, and thrice through the panel with stiffener (T-shaped), namely, the time it takes for the reflected wave from the rear surface of the structure to reach the wet surface. Between two arbitrary time ticks, small oscillations exist. This is as a result of the superposition of the incident wave and the reflected wave from the interface of medium S1 and medium S2. As such, we can characterize the response of the panel with stiffener (T-shaped) by the drops in specific time tick and the small oscillations between the time ticks. The response clearly shows the wave propagation process in T-shaped plate. In Figure 13(b), the characteristics induced by superposition of multiwaves can be expressed as the jumps at specific time ticks and small oscillations between time ticks. Circles in Figure 13(b) are the maximum value of the velocity and the velocity when the wet surface gets cavitated, that is, 1.422 and 1.404. The panel with stiffener (T-shaped) gets its peak velocity before cavitation occurs on the wet surface. As shown, the analytical solution is a little larger than the FE prediction and the differences are reasonably acceptable. As shown in Figure 13, Taylor’s solutions for rigid plate equivalent to panel with stiffener (T-shaped) have similar trend as the analytical and FE prediction. However, Taylor’s solution cannot show the stress wave propagation process in the specific structure and it cannot recognize the minor change in the structure.In Figure14(a), we can see the obvious small jumps at the specific time (t=0.53,t=0.1068) which is the result of the enlargement effect when the stress wave is transmitted from medium S2 to medium S3. The details of the enhancement of pressure will be expressed in Section 4.2. The sharp drop of pressure occurs when the rarefaction wave from the rear surface of the structure reaches the wet surface. Similar to the wave propagation process in the panel with stiffener (T-shaped), the corresponding response of the panel with girder (I-shaped) is characterized by big drops at specific time tick and small oscillations between time ticks. The differences are that the wave travel in the panel with girder (I-shaped) takes longer time to finish one trip through the structure and the small jumps take place due to the existence of the medium S3. In Figure 14(b), we also can describe the trends in nondimensional time history of velocity as jumps at specific time ticks and small oscillations between time ticks, as well as the small drops owing to the existence of the medium S3. The panel with girder (I-shaped) gets its peak velocity before cavitation occurs on the wet surface (see circles in Figure 14(b)). By comparing the velocities in Figures 12(b), 13(b), and 14(b), we can conclude that the velocity of the wet surface is considerably small if the structure is heavy or strong. Correspondingly, the rarefaction wave radiated by the motion of wet surface is small. Thus, the pressure on the wet surface decreases slowly (see Figure 14(a)). As shown in Figure 14, Taylor’s solutions for rigid plate equal to the panel with girder (I-shaped) have a similar trend as the analytical and FE predictions. However, there is a big difference between Taylor’s solution and FE predictions. It can be seen that Taylor’s solution cannot capture the reflected wave from the varying cross section of the structure.


### 3.2.3. Impulse Histories on Wet Surface

The intensity of loads impacted on the structure consequent to underwater explosion, as quantified by the specific impulse (per unit area) applied on the wet surface, is deduced by solving the following equation:(26)Ikt=∫0tpwetkτdτ,k=1,2,3,where k=1,2,3 refers to the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively.Define the maximum incident impulseI0 (integral of the incident wave) per unit quasistatic state as I0=2p0θ. The severity of structural loading consequent to underwater explosion can be quantified by the maximum nondimensional impulse, I-max=Imax/I0, applied on the load surface of the structure.The analytical solutions for the nondimensional impulseI-k of three substructures agree well with the FE predictions and Taylor’s solutions (see Figure 15). As the thickness of the structure increases, the plate takes longer time to reach its maximum nondimensional impulse. This is due to the fact that the larger thickness of the structure will prolong the propagating time and thus release higher energy to the plate. The comparisons of the analytical results and FE predictions of maximum nondimensional impulse I-maxk for the three substructures are presented in Table 2. We can see that the numerical predictions of I-maxk are little smaller than the analytical results. This is probably caused by the numerical dissipation. As shown in the second column of Table 2, the normalized impulses I-maxk for the three types of plates equal to 6.81%, 12.17%, and 39.87%, respectively, and are all lower than unity. Moreover, comparison between the three plates shows that the specific impulse applied on the wet surface of the panel with girder (I-shaped) is larger than that of the panel with stiffener (T-shaped) and the simple panel, indicating that more momentum is transmitted to the panel with girder (I-shaped). This is due to the fact that the pressure on the wet surface of the panel with girder (I-shaped) decreases slowly and therefore the impulse impacting the plate will increase.Table 2
Analytical and FE predictions of normalized impulse.
Case
Analytical results
FE predictions
Error
Simple panel
0.0681
0.0673
1.25%
The panel with stiffener (T-shaped)
0.1217
0.1184
2.71%
The panel with girder (I-shaped)
0.3987
0.3977
0.25%Figure 15
The analytical solutions and FE predictions of the nondimensional impulse for (a) simple panel; (b) the panel with stiffener (T-shaped); (c) the panel with girder (I-shaped).


## 3.2.1. Closing Front and Breaking Front

Figure11 presents the cavitation time t-c as a function of the location x- and the velocity v-cav in the cavitated water at t-cx for the three substructures. The minimum value on these curves (circled in Figure 11(a)) is the location and time of the onset of first cavitation. Two branches of each curve originated from this point, corresponding to propagation of two breaking fronts traveling away from (left of x-c) and towards (right of x-c) the wet surface. Note that the cavitation time is linear with x- on the left of x-c, except very close to the first cavitation location.Figure 11
Analytical predictions and FE prediction of the spatial variation of (a) the time to cavitation and (b) the velocityv-cav in the cavitated water at t-cx. These results are for cases in Table 1, which are the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped).
(a)
(b)For the simple panel, cavitation initiates at a very small distance from the wet surface and gives rise to a breaking front traveling away from the interface, while the breaking front approaching the interface cannot be arrested before reaching the wet surface. For the panel with stiffener (T-shaped) and the panel with girder (I-shaped), the first cavitation is located within the fluid and the breaking front traveling towards the structure would reverse its motion and become the closing front. Figure11(a) shows the first cavitation time for the water column. For example, the jump in the curve of the panel with girder (I-shaped) indicates that when the first cavitation occurs at x¯=-0.298, the breaking front propagated away from the wet surface. As time elapses, another cavitation takes place at the interface at time t-c0. And this breaking front also propagated away from the wet surface. These two breaking fronts generate from different place, while they travel at the same velocity. The jump means that the cavitation will occur twice at the same displacement and different time. We can see a similar phenomenon in the case of the panel with stiffener (T-shaped).As shown in Figure11(a), the trends of two breaking fronts in the analytical results agree well with FE results. The analytical results in Figure 11(a) only show the first cavitation time for the water column. Thus, when the cavitation occurs at the same displacement and different time, the curve will “jump” without showing the closing front. And this leads to the discrepancy between FE and analytical predictions. From Figure 11(a), we can see that t-c0<trest/θ-xrest/cwθ indicates that the cavitation occurs on the wet surface before the wave train emanating from the closing front reaches the plate.The expressions (13) for the plate velocity and pressure are still valid when t-<t-c0. In Figure 11(b), the velocity v-cav in the cavitated water at t-cx for the three substructures is shown and small difference is found between the analytical solutions and the FE predictions. The particle velocity v-cav0 is the corresponding velocity on the wet surface when cavitation occurs. The value of v-cav0 for the panel, the panel with the stiffener (T-shaped), and the panel with girder (I-shaped) are 1.568, 1.407, and 0.6657, respectively. The variance from the analytical results for the different shape plates is large. This is due to the fact that the panel with girder (I-shaped) has a larger unit mass.


## 3.2.2. Pressure and Velocity Histories on Wet Surface

Figure12(a) shows nondimensional time histories of nondimensional wet surface pressure p-wet1, while the corresponding time histories of nondimensional velocity v-wet1 are presented in Figure 12(b). As discussed above, in this particular case the first cavitation takes place at a location very close to the wet surface. We can then calculate the analytical solution of the pressure and velocity by using (18) after the wet surface is cavitated. In particular, we investigated the dynamic response of a steel plate by a three-method analytical solution, FEM results, and Taylor’s solution and then did a comparison of the predictions. Both Figures 12(a) and 12(b) show that the analytical solutions are adequately predicted by FE results and Taylor’s solutions [8]. On the whole process, the analytical solution agrees well with the FE prediction for the pressure and velocity. It is to be noted that the small jumps of analytical solutions indicate the reflections and transmissions of the pressure wave. When the stress wave propagation and the cavitation effect are considered, the analytical results are more accurate than Taylor’s results and FE predictions.Figure 12
Nondimensional time histories of (a) nondimensional wet surface pressure and (b) nondimensional velocity for the simple panel.
(a)
(b)To reveal the fluid-structure interaction of the double-bottom structure, the nondimensional time histories of wet surface pressure and velocity for the panel with stiffener (T-shaped) and the panel with girder (I-shaped) are shown in Figures13 and 14, respectively. The FE predictions are presented to validate the analytical solutions. As shown, the analytical solutions agree considerably well with the FE predictions. The oscillations in FE predictions can be explained by the fact that on one hand the FE algorithm for the thin steel plate accounts for the FE predictions that have large numerical diffusions; on the other hand, the FE numerical models are three-dimensional, including the effect of the transverse wave. We note that if cavitation occurs on the wet surface, the velocity of interface would be given as an average. After the wet surface gets cavitated, the process becomes very complex to solve analytically. Thus, the velocity is calculated as an average. Hence, we can conclude that the analytical solutions are validated in these cases.Figure 13
Nondimensional time histories of (a) nondimensional wet surface pressure and (b) nondimensional velocity for the panel with stiffener (T-shaped).
(a)
(b)Figure 14
Nondimensional time histories of (a) nondimensional wet surface pressure and (b) nondimensional velocity for the panel with girder (I-shaped).
(a)
(b)Note that the obvious drops in the analytical solutions indicate the superposition of the reflected wave from the rear surface of the structure and the incident wave. For the panel with stiffener (T-shaped), this reflected wave is the rarefaction wave which causes the great drop in the pressure and the decrement of decay constant (see Figure13(a)). The specific time ticks circled in Figure 13(a), t=0.05, t=0.105, and t=0.1583, are the time ticks of stress wave traveling once, twice, and thrice through the panel with stiffener (T-shaped), namely, the time it takes for the reflected wave from the rear surface of the structure to reach the wet surface. Between two arbitrary time ticks, small oscillations exist. This is as a result of the superposition of the incident wave and the reflected wave from the interface of medium S1 and medium S2. As such, we can characterize the response of the panel with stiffener (T-shaped) by the drops in specific time tick and the small oscillations between the time ticks. The response clearly shows the wave propagation process in T-shaped plate. In Figure 13(b), the characteristics induced by superposition of multiwaves can be expressed as the jumps at specific time ticks and small oscillations between time ticks. Circles in Figure 13(b) are the maximum value of the velocity and the velocity when the wet surface gets cavitated, that is, 1.422 and 1.404. The panel with stiffener (T-shaped) gets its peak velocity before cavitation occurs on the wet surface. As shown, the analytical solution is a little larger than the FE prediction and the differences are reasonably acceptable. As shown in Figure 13, Taylor’s solutions for rigid plate equivalent to panel with stiffener (T-shaped) have similar trend as the analytical and FE prediction. However, Taylor’s solution cannot show the stress wave propagation process in the specific structure and it cannot recognize the minor change in the structure.In Figure14(a), we can see the obvious small jumps at the specific time (t=0.53,t=0.1068) which is the result of the enlargement effect when the stress wave is transmitted from medium S2 to medium S3. The details of the enhancement of pressure will be expressed in Section 4.2. The sharp drop of pressure occurs when the rarefaction wave from the rear surface of the structure reaches the wet surface. Similar to the wave propagation process in the panel with stiffener (T-shaped), the corresponding response of the panel with girder (I-shaped) is characterized by big drops at specific time tick and small oscillations between time ticks. The differences are that the wave travel in the panel with girder (I-shaped) takes longer time to finish one trip through the structure and the small jumps take place due to the existence of the medium S3. In Figure 14(b), we also can describe the trends in nondimensional time history of velocity as jumps at specific time ticks and small oscillations between time ticks, as well as the small drops owing to the existence of the medium S3. The panel with girder (I-shaped) gets its peak velocity before cavitation occurs on the wet surface (see circles in Figure 14(b)). By comparing the velocities in Figures 12(b), 13(b), and 14(b), we can conclude that the velocity of the wet surface is considerably small if the structure is heavy or strong. Correspondingly, the rarefaction wave radiated by the motion of wet surface is small. Thus, the pressure on the wet surface decreases slowly (see Figure 14(a)). As shown in Figure 14, Taylor’s solutions for rigid plate equal to the panel with girder (I-shaped) have a similar trend as the analytical and FE predictions. However, there is a big difference between Taylor’s solution and FE predictions. It can be seen that Taylor’s solution cannot capture the reflected wave from the varying cross section of the structure.


## 3.2.3. Impulse Histories on Wet Surface

The intensity of loads impacted on the structure consequent to underwater explosion, as quantified by the specific impulse (per unit area) applied on the wet surface, is deduced by solving the following equation:(26)Ikt=∫0tpwetkτdτ,k=1,2,3,where k=1,2,3 refers to the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively.Define the maximum incident impulseI0 (integral of the incident wave) per unit quasistatic state as I0=2p0θ. The severity of structural loading consequent to underwater explosion can be quantified by the maximum nondimensional impulse, I-max=Imax/I0, applied on the load surface of the structure.The analytical solutions for the nondimensional impulseI-k of three substructures agree well with the FE predictions and Taylor’s solutions (see Figure 15). As the thickness of the structure increases, the plate takes longer time to reach its maximum nondimensional impulse. This is due to the fact that the larger thickness of the structure will prolong the propagating time and thus release higher energy to the plate. The comparisons of the analytical results and FE predictions of maximum nondimensional impulse I-maxk for the three substructures are presented in Table 2. We can see that the numerical predictions of I-maxk are little smaller than the analytical results. This is probably caused by the numerical dissipation. As shown in the second column of Table 2, the normalized impulses I-maxk for the three types of plates equal to 6.81%, 12.17%, and 39.87%, respectively, and are all lower than unity. Moreover, comparison between the three plates shows that the specific impulse applied on the wet surface of the panel with girder (I-shaped) is larger than that of the panel with stiffener (T-shaped) and the simple panel, indicating that more momentum is transmitted to the panel with girder (I-shaped). This is due to the fact that the pressure on the wet surface of the panel with girder (I-shaped) decreases slowly and therefore the impulse impacting the plate will increase.Table 2
Analytical and FE predictions of normalized impulse.
Case
Analytical results
FE predictions
Error
Simple panel
0.0681
0.0673
1.25%
The panel with stiffener (T-shaped)
0.1217
0.1184
2.71%
The panel with girder (I-shaped)
0.3987
0.3977
0.25%Figure 15
The analytical solutions and FE predictions of the nondimensional impulse for (a) simple panel; (b) the panel with stiffener (T-shaped); (c) the panel with girder (I-shaped).


## 4. Sensitivity of the Response to Nondimensional Geometric Parameters

Having established the accuracy of the analytical predictions, the analytical models are now used to determine the sensitivity of the plate response to the nondimensional parameters ofα-k=lf/tm, β-k=dm/tf, and γ-k=lb/tm for the panel with stiffener (T-shaped) and the panel with girder (I-shaped), where the superscript k (k=1,2,3) represents the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively.


### 4.1. Response of the Panel with Stiffener (T-Shaped)

α
-
2 is the basic parameter of the panel with stiffener (T-shaped). In this case, to investigate the effect of the single parameter α-2 on the design of the panel with stiffener (T-shaped), we set β-2=7. The time histories of nondimensional wet pressure p-wet2 of different α-2 are presented in Figure 16. For a special case of when α-2 is 1, the problem is the same as the response for the simple panel with a sum thickness of mediums S1 and S2. In this case, only one characteristic is shown, namely, the drops at specific time ticks. When α-2 is greater than 1, the characteristics described in Section 3.2.2 for the panel with stiffener (T-shaped) are shown. When α-2 is less than 10, we find that increasing α-2 will lead to larger drops of the pressure, as a result of the smaller decay constant and the shortened cavitation time on the wet surface. When α-2 is greater than 10, there is little difference between time histories of pressure with different α-2. When α-2→∞, the nondimensional histories of pressure approach those of the simple panel. This can be explained as follows: as α-2 increases the corresponding specific mass (per unit area) decreases dramatically and approaches the specific mass of simple panel, which greatly weakens the effect of the stiffener on the associated simple panel.Figure 16
The nondimensional time histories of the nondimensional pressurep-wet2 for the panel with stiffener (T-shaped) with the parameters (a) α-2 = 1~10, β-2=7; (b) α-2 = 10~100, β-2=7.
(a)
(b)In Figure17, the analytical results of nondimensional maximum impulse I-max2 are plotted against the normalized parameters α-2 for the panel with stiffener (T-shaped). Taylor’s predictions for rigid plates of mass equal to this case are also plotted in Figure 17. Both results predict a similar trend. The results show that α-2 increases rapidly with reduction of I-max2. We note that, for α-2≫10, rate of change of I-max2 with α-2 slows down, and the value of I-max2 drops to 0.15, which is approximately the same value as I-max1 for simple panel. According to the results plotted in Figures 16 and 17, we can conclude that the effect of the stiffener on the associated simple panel is significant only in a very small region (α-2<20). To enhance the significant effect of stiffeners on the associated panel, we can add the number of stiffeners with the same panel.Figure 17
Sensitivity of the nondimensional impulseI-max2 to the nondimensional parameter α-2 of the panel with stiffener (T-shaped).To investigate the effect of the single parameterβ-2 on the dynamic response of the panel with stiffener (T-shaped) we set α-2=15. Figure 18(a) presents the dependence of p-wet2 upon β-2. According to the times the input is wave transmitted through the panel with stiffener (T-shaped) (specific time ticks), the time histories of nondimensional pressure can be divided into few phases. As β-2 increase, the time the wave travels through T-shaped once increases linearly. That is, the time period between two time ticks increases as β-2 increase. The trends in each phase are similar and the plate with small β-2 is easier to get cavitated. Figure 18(b) shows the sensitivity of the nondimensional impulse I-max2 to the normalized parameter β-2. The analytical solution and Taylor’s solution have the same trend. The I-max2 increases linearly with β-2. It means that the larger time period between two time ticks can enlarge the impulse on the wet surface. The trend is the consequence of the fact that the ability of the panel with stiffener (T-shaped) to transport energy increases with increasing dm. Therefore, to enhance the effect of stiffeners on the associated plate, another available way is to increase the value of β-2.Figure 18
(a) The normalized time histories of the nondimensional pressurep-wet2 for the panel with stiffener (T-shaped) with the parameters β-2 and α-2=15; (b) sensitivity of the nondimensional impulse I-max2 to the normalized parameter β-2 for the panel with stiffener (T-shaped), α-2=15.
(a)
(b)


### 4.2. Response of the Panel with Girder (I-Shaped)

Figure19 shows that the nondimensional pressure histories for parameter α-3 for some values between 1 and 70. In these cases, γ-3 equals 2 and β-3 equals 80. It is found that the nondimensional time histories of pressure on the wet surface can be divided into three phases. In the first phase (before the rarefaction wave travels back to the wet surface), the pressure linearly decreases with increasing α-3. In the second phase (after the rarefaction wave reaches the wet surface and before cavitation occurs), the pressure drops sharply; we name it the “cut-off” effect. When α-3 is less than 10, we can observe that the “cut-off” effect is more and more significant as α-3 increases. When α-3 is greater than 10, the “cut-off” effect directly leads to cavitation on the wet surface. The phenomenon of second phase can be explained as follows: due to increasing α-3, the corresponding solid mass (per unit area) decreases, and thus the velocity of wet surface gets larger, resulting in cavitation. Thus, the pressure on the wet surface decreases quickly. In the third phase (after the cavitation occurs), we can see that the cavitation time is shortened when α-3 increases.Figure 19
The normalized time histories of the nondimensional pressurep-wet3 for the panel with girder (I-shaped) with the parameter zone (a) α-3 = 1~10, β-3=80, and γ-3=2; (b) α-3 = 10~70,  β-3=80, and γ-3=2.
(a)
(b)By integrating the pressure histories, we can obtain the analytical maximum nondimensional impulse (see Figure20). It is noted that when α-3 is less than 3, the nondimensional impulse is greater than unity. This is because the total impulse on the wet surface is 2p0θ for the stationary rigid plate case. Therefore, for the heavy plates, the velocity is small and results in I/2p0θ higher than 0.5. Similarly, the significant effect of girders on the associated panel is limited to the parameter α-3. When α-3>10, the maximum impulse exponentially decays.Figure 20
Sensitivity of the nondimensional impulseI-max3 to the normalized parameter α-3 for the panel with girder (I-shaped).As shown in Figure21(a), the γ-3 is set to be 3, and α-3 is 7. The time histories of nondimensional pressure p-wet3 still can be taken as three phases. The phases are divided based on the time the input wave travels once through the panel with girder (I-shaped). As β-3 increases, the time for the wave to propagate once through the panel with girder (I-shaped) increases linearly. The trends in the whole process with different β-3 are the same and share similar characteristics: the great drop occurs at the specific time tick, small oscillations occur between time ticks, and there is a little jump due to medium S3. Looking at the results plotted in Figures 19 and 21(a), we found out that α-3 plays a key role in the first phase and β-3 determines the time ticks of the phase. In the third phase, we see that the plate with small β-3 is easier to get cavitated on wet surface. The corresponding maximum impulse I-max3 is shown in Figure 21(b). The impulse is approximately linear with increasing β-3, indicating that large β-3 increases the effect of girders on the associated panel.Figure 21
(a) The normalized time histories of the nondimensional pressurep-wet3 for the panel with girder (I-shaped) with the parameter β-3. (b) Sensitivity of the nondimensional impulse I-max3 to the normalized parameter β-3 for the panel with girder (I-shaped).
(a)
(b)To investigate the effect ofγ-3, we set α-3 equal to 7 and β-3 equal to 80. The pressure histories for various γ-3 are shown in Figure 22(a). As we can see, γ-3 has a key effect on the jump of the pressure time histories. This stress enhancement occurs when it comes to large impedance ratio on the interface. When the wave transmits from medium S2 to medium S3, due to the varying section, a strengthened reflected wave is generated on the interface of medium S2 and medium S3. And this strengthened reflected wave reaches the wet surface before the rarefaction wave from the rear surface. As such, an interesting phenomenon occurs; there are small jumps before the great drop in pressure, except for the case of γ-3=1. When γ-3=1, the problem becomes the same as that used to solve the response for the panel with stiffener (T-shaped). In this case, the small jump before the big drop in time histories of pressure would never happen. From Figure 22(a) it can be seen that the jump is more prominent with the increasing γ-3. In Figure 22(b), the nondimensional impulse I-max3 is only weakly sensitive to variations of γ-3.Figure 22
(a) The normalized time histories of the nondimensional pressurep-wet3 for the panel with girder (I-shaped) with the parameter γ-3 = 1~10. (b) Sensitivity of the nondimensional impulse I-max3 to the normalized parameter γ-3 for the panel with girder (I-shaped).
(a)
(b)Taylor’s predictions for monolithic plates of equivalent mass are also shown in Figures20, 21(b), and 22(b). When α-3<3 and α-3>10, the differences between the analytical solution and Taylor’s prediction are prominent. Due to the fact that the strengthened reflected wave in the panel with girder (I-shaped) will cause small jumps in time history of pressure, the nondimensional impulse I-max3 of analytical solutions is slightly higher than these of Taylor’s predictions. However, they share a similar trend.


## 4.1. Response of the Panel with Stiffener (T-Shaped)

α
-
2 is the basic parameter of the panel with stiffener (T-shaped). In this case, to investigate the effect of the single parameter α-2 on the design of the panel with stiffener (T-shaped), we set β-2=7. The time histories of nondimensional wet pressure p-wet2 of different α-2 are presented in Figure 16. For a special case of when α-2 is 1, the problem is the same as the response for the simple panel with a sum thickness of mediums S1 and S2. In this case, only one characteristic is shown, namely, the drops at specific time ticks. When α-2 is greater than 1, the characteristics described in Section 3.2.2 for the panel with stiffener (T-shaped) are shown. When α-2 is less than 10, we find that increasing α-2 will lead to larger drops of the pressure, as a result of the smaller decay constant and the shortened cavitation time on the wet surface. When α-2 is greater than 10, there is little difference between time histories of pressure with different α-2. When α-2→∞, the nondimensional histories of pressure approach those of the simple panel. This can be explained as follows: as α-2 increases the corresponding specific mass (per unit area) decreases dramatically and approaches the specific mass of simple panel, which greatly weakens the effect of the stiffener on the associated simple panel.Figure 16
The nondimensional time histories of the nondimensional pressurep-wet2 for the panel with stiffener (T-shaped) with the parameters (a) α-2 = 1~10, β-2=7; (b) α-2 = 10~100, β-2=7.
(a)
(b)In Figure17, the analytical results of nondimensional maximum impulse I-max2 are plotted against the normalized parameters α-2 for the panel with stiffener (T-shaped). Taylor’s predictions for rigid plates of mass equal to this case are also plotted in Figure 17. Both results predict a similar trend. The results show that α-2 increases rapidly with reduction of I-max2. We note that, for α-2≫10, rate of change of I-max2 with α-2 slows down, and the value of I-max2 drops to 0.15, which is approximately the same value as I-max1 for simple panel. According to the results plotted in Figures 16 and 17, we can conclude that the effect of the stiffener on the associated simple panel is significant only in a very small region (α-2<20). To enhance the significant effect of stiffeners on the associated panel, we can add the number of stiffeners with the same panel.Figure 17
Sensitivity of the nondimensional impulseI-max2 to the nondimensional parameter α-2 of the panel with stiffener (T-shaped).To investigate the effect of the single parameterβ-2 on the dynamic response of the panel with stiffener (T-shaped) we set α-2=15. Figure 18(a) presents the dependence of p-wet2 upon β-2. According to the times the input is wave transmitted through the panel with stiffener (T-shaped) (specific time ticks), the time histories of nondimensional pressure can be divided into few phases. As β-2 increase, the time the wave travels through T-shaped once increases linearly. That is, the time period between two time ticks increases as β-2 increase. The trends in each phase are similar and the plate with small β-2 is easier to get cavitated. Figure 18(b) shows the sensitivity of the nondimensional impulse I-max2 to the normalized parameter β-2. The analytical solution and Taylor’s solution have the same trend. The I-max2 increases linearly with β-2. It means that the larger time period between two time ticks can enlarge the impulse on the wet surface. The trend is the consequence of the fact that the ability of the panel with stiffener (T-shaped) to transport energy increases with increasing dm. Therefore, to enhance the effect of stiffeners on the associated plate, another available way is to increase the value of β-2.Figure 18
(a) The normalized time histories of the nondimensional pressurep-wet2 for the panel with stiffener (T-shaped) with the parameters β-2 and α-2=15; (b) sensitivity of the nondimensional impulse I-max2 to the normalized parameter β-2 for the panel with stiffener (T-shaped), α-2=15.
(a)
(b)


## 4.2. Response of the Panel with Girder (I-Shaped)

Figure19 shows that the nondimensional pressure histories for parameter α-3 for some values between 1 and 70. In these cases, γ-3 equals 2 and β-3 equals 80. It is found that the nondimensional time histories of pressure on the wet surface can be divided into three phases. In the first phase (before the rarefaction wave travels back to the wet surface), the pressure linearly decreases with increasing α-3. In the second phase (after the rarefaction wave reaches the wet surface and before cavitation occurs), the pressure drops sharply; we name it the “cut-off” effect. When α-3 is less than 10, we can observe that the “cut-off” effect is more and more significant as α-3 increases. When α-3 is greater than 10, the “cut-off” effect directly leads to cavitation on the wet surface. The phenomenon of second phase can be explained as follows: due to increasing α-3, the corresponding solid mass (per unit area) decreases, and thus the velocity of wet surface gets larger, resulting in cavitation. Thus, the pressure on the wet surface decreases quickly. In the third phase (after the cavitation occurs), we can see that the cavitation time is shortened when α-3 increases.Figure 19
The normalized time histories of the nondimensional pressurep-wet3 for the panel with girder (I-shaped) with the parameter zone (a) α-3 = 1~10, β-3=80, and γ-3=2; (b) α-3 = 10~70,  β-3=80, and γ-3=2.
(a)
(b)By integrating the pressure histories, we can obtain the analytical maximum nondimensional impulse (see Figure20). It is noted that when α-3 is less than 3, the nondimensional impulse is greater than unity. This is because the total impulse on the wet surface is 2p0θ for the stationary rigid plate case. Therefore, for the heavy plates, the velocity is small and results in I/2p0θ higher than 0.5. Similarly, the significant effect of girders on the associated panel is limited to the parameter α-3. When α-3>10, the maximum impulse exponentially decays.Figure 20
Sensitivity of the nondimensional impulseI-max3 to the normalized parameter α-3 for the panel with girder (I-shaped).As shown in Figure21(a), the γ-3 is set to be 3, and α-3 is 7. The time histories of nondimensional pressure p-wet3 still can be taken as three phases. The phases are divided based on the time the input wave travels once through the panel with girder (I-shaped). As β-3 increases, the time for the wave to propagate once through the panel with girder (I-shaped) increases linearly. The trends in the whole process with different β-3 are the same and share similar characteristics: the great drop occurs at the specific time tick, small oscillations occur between time ticks, and there is a little jump due to medium S3. Looking at the results plotted in Figures 19 and 21(a), we found out that α-3 plays a key role in the first phase and β-3 determines the time ticks of the phase. In the third phase, we see that the plate with small β-3 is easier to get cavitated on wet surface. The corresponding maximum impulse I-max3 is shown in Figure 21(b). The impulse is approximately linear with increasing β-3, indicating that large β-3 increases the effect of girders on the associated panel.Figure 21
(a) The normalized time histories of the nondimensional pressurep-wet3 for the panel with girder (I-shaped) with the parameter β-3. (b) Sensitivity of the nondimensional impulse I-max3 to the normalized parameter β-3 for the panel with girder (I-shaped).
(a)
(b)To investigate the effect ofγ-3, we set α-3 equal to 7 and β-3 equal to 80. The pressure histories for various γ-3 are shown in Figure 22(a). As we can see, γ-3 has a key effect on the jump of the pressure time histories. This stress enhancement occurs when it comes to large impedance ratio on the interface. When the wave transmits from medium S2 to medium S3, due to the varying section, a strengthened reflected wave is generated on the interface of medium S2 and medium S3. And this strengthened reflected wave reaches the wet surface before the rarefaction wave from the rear surface. As such, an interesting phenomenon occurs; there are small jumps before the great drop in pressure, except for the case of γ-3=1. When γ-3=1, the problem becomes the same as that used to solve the response for the panel with stiffener (T-shaped). In this case, the small jump before the big drop in time histories of pressure would never happen. From Figure 22(a) it can be seen that the jump is more prominent with the increasing γ-3. In Figure 22(b), the nondimensional impulse I-max3 is only weakly sensitive to variations of γ-3.Figure 22
(a) The normalized time histories of the nondimensional pressurep-wet3 for the panel with girder (I-shaped) with the parameter γ-3 = 1~10. (b) Sensitivity of the nondimensional impulse I-max3 to the normalized parameter γ-3 for the panel with girder (I-shaped).
(a)
(b)Taylor’s predictions for monolithic plates of equivalent mass are also shown in Figures20, 21(b), and 22(b). When α-3<3 and α-3>10, the differences between the analytical solution and Taylor’s prediction are prominent. Due to the fact that the strengthened reflected wave in the panel with girder (I-shaped) will cause small jumps in time history of pressure, the nondimensional impulse I-max3 of analytical solutions is slightly higher than these of Taylor’s predictions. However, they share a similar trend.


## 5. The Interaction between a Double-Bottom Structure and Its Substructures

In Section3, we verified the analytical model for the substructures of the double-bottom structure: the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped). In Section 4, we demonstrated the effects of the nondimensional geometric parameters on the substructures. In this section, we conduct a case study of the substructures, by analytical solution and FE simulation, to obtain the approximate analytical solution to predict the impulse for the complex double-bottom structure system. That is, based on the results from FE method, we try to figure out the interaction between the substructures and a double-bottom structure by using the nondimensional parameters α-k, β-k, and γ-k.It is noted that the parameterα-k for the panel with stiffener (T-shaped) and the panel with girder (I-shaped) is key to relate the response of the double-bottom structure and the substructures, which can be expressed as(27)I-D=I-1·SΓ-∑k=23Ψα-k+∑k=23I-k·Ψα-kSΓ,where the superscript k represent the substructures: the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped), respectively. Γ represents the wet surface of the double-bottom structure and SΓ is the corresponding area. Ψα-k is the function which depends on the parameter α-k. For simplification, we set Ψα-k=α-k·Nk, k=2,3. N is the sum of the corresponding substructures in the double-bottom structure cell. A cell of a double-bottom structure is defined as the structure assembly between two girders (see Figure 23).Figure 23
The sketch of the double-bottom structure cell subjected to underwater explosion.To obtain the solution of (27), an analytical solution and FE simulation are carried out. The FE modeling process is the same as in Section 3.1. Considering the double-bottom structure (see Figure 23), the stiffeners and girders have the same value as shown in Table 1, and the distance between the girders is 1 m. Even though the double-bottom structure is very complex, to validate the effective model, the FE prediction of impulse of double-bottom structure is solved as an average; namely, (28)I-DFEt=∫Γ∫0tp-wetFEτdτds2∫Γds,where Γ represent the wet surface.Figure24(a) shows the comparison between the FE prediction (obtained from (28)) and the analytical solution (obtained from (27)) for the nondimensional impulse maximum I-D,max. In this figure, the FE prediction is a constant. To get the same value as the FE prediction, we have to change the value of  α-k in (27). By comparing the curves in Figure 24(a), we can easily solve (27) by setting α-k=15. The definition of the impulse of a double-bottom structure is the summation of the product of impulses of the substructures and the area percentage of the substructures in a double-bottom structure. As α-k increases, the area percentage of the T-shaped panel and the I-shaped panel in a double-bottom structure increases. Based on the above factors, in Figure 24(a), as α-k increases, the maximum impulse increases.Figure 24
Comparison of analytical solution and FE predictions of (a) the nondimensional maximum impulse and (b) the nondimensional impulse time histories.
(a)
(b)By settingα-k=15, the nondimensional time histories of impulse are showed in Figure 24(b). The analytical predictions agree well with FE predictions, except for the small difference from time 0.4~0.8 of 1.13%. For the difference, the explanations are as follows: the response of the double-bottom is based on the results of substructures; and the cavitation time on the wet surface of substructures is different. In the analytical solution, when cavitation occurs on the wet surface, we assume that the pressure on the wet surface for the corresponding substructures is zero and never changes. After cavitation occurs, the analytical solutions obtained by this assumption cannot present the exact solution in water and thus result in the difference. Actually, the cavitation effect in water is very complex. Compared to the analytical solutions, the FE predictions on water cavitation are more accurate. Though the FE method has the characteristic of being general, the analytical method reveals the underlying physics more clearly and is more convenient to perform a parameter effect analysis. An error of 1.13% is acceptable, and hence we can still make the conclusion that the approximate analytical solution for double-bottom structure is validated in this case. In other words, when we set α-k=15, we can obtain a reasonable impulse for an arbitrary double-bottom structure by (27).In summary, the steps of the analytical method for an arbitrary double-bottom structure in this paper are as follows: the first step is to divide a double-bottom structure into several substructures according to the nondimensional parametersα-k where α-k=15; the second step is to solve the dynamic response of substructures according to the methods in Section 2; the final step is to obtain the dynamic response of double-bottom structure by solving (27) where α-k=15.


## 6. Conclusion

In this paper, the response of double-bottom structure to underwater explosion has been investigated by analytical approaches. We have constructed and validated theoretical models for the substructures (the simple panel, the panel with stiffener (T-shaped), and the panel with girder (I-shaped)) loaded by planar, exponentially decaying underwater shock wave. Based on the dynamic response of the substructures, we established the approximate analytical models which are able to predict the response of the double-bottom structure to an underwater explosion. The shock propagation mechanism in the substructures, FSI effects, and cavitation phenomenon were considered in detail. According to the analysis, the main conclusions are as follows:(1)
The analytical models for substructures loaded by exponentially decaying water blast were validated by comparing their predictions to those obtained from fully coupled 3D FE simulations and good agreements were found. However, Taylor’s predictions for rigid plates of equivalent areal mass including the effects of cavitation fail to capture the reflected wave from the varying cross section of the structure. It shows a potential application field for the analytical method proposed in this paper.(2)
Writing the expression of dynamic response on the wet surface in nondimensional form shows that the response of substructures to underwater blast depends on three parameters, namely,α-k, β-k, and γ-k. α-k is the key parameter defining the effect of the stiffeners or girders on the associated panel. It shows that the effect of the stiffeners or girders on the associated panel is limited to a small region α-k<20. The way to enhance this effect is to increase the number of the corresponding substructures. β-k is the parameter to prolong the wave transmission time through the structure. Increasing the value of β-k is another available way to enhance the effect of the stiffeners or girders on the associated panel.(3)
α-k is key in relating the response of the substructures and the double-bottom structure. It was shown that by taking α-k as 15 very small difference is found between the analytical solution and the FE prediction. The relationship can be expressed as I-D=I-1·SΓ-∑k=23Ψα-k+∑k=23I-k·Ψα-k/SΓ. The parameter β-2 of the panel with stiffener (T-shaped) and β-3 of the panel with girder (I-shaped) play an important role on FSI dynamic response for double-bottom structure.The major limitation of the analytical method in the present paper is the assumption that the structure is completely elastic, neglecting plasticity, buckling, and fracture. When a double-bottom structure subjected to underwater blast, the structure would suffer severe deformation. However, the effects of plasticity and buckling cannot be solved by this method. The spherical fashion of the reflected wave at the abrupt change in cross section and the propagation of flexural wave due to bending cannot be captured. Therefore, more detailed calculations of the plastic wave propagation in 2D and 3D need to be done to reveal the fluid-structure effect on a double-bottom subjected to underwater explosion. This will from part of our future research.Despite these limitations, the current double-bottom model and analyses presented in this paper are useful in understanding the fluid-structure interaction, and the elastic stress wave propagation process in complex structures before the impulse cause the bending or the deformation of the whole structure.



---

*Source: 7368624-2016-06-29.xml*