Centrifuge experiment and numerical study on the dynamic response of air-backed plate to underwater explosion

This work compares experimental and numerical results concerning the elastic dynamic response of air-backed plate to underwater explosion. Experiments were performed in a centrifuge, both the shock loadings and structure responses were tested, and the bubble oscillation considering centrifugal force was predicted with Geer and Hunter model. The experimental and numerical results illustrates that the centrifugal force have no effect on shockwave for the short duration of action, and with the increase of centrifugal force, the bubble pulse motion cycle, maximum bubble radius, and peak pressure of bubble pulse decreased; both the peak of high-frequency and low-frequency response of the plate decline slightly when alpha damping rising ,the high-frequency response mode almost have no change with different alpha damping ,while, the low-frequency response tends to be an obvious oscillation waveform when the alpha damping is too low; although the difference of the peak acceleration owing to shock wave and bubble reaches to tens of times, the strain responses caused by the two factors are very close, which enucleates that the train response is related to shock energy and has a higher sensitivity to time accumulation; both the peak acceleration and peak strain due to the shock wave remain steady and the two kinds of response by bubble pulse decrease with the centrifugal force increasing; when the target position is more close to the area that the shock wave vertically impact ,the linear correlation between peak velocity and shock factor is much better.


Introduction
Understanding and prediction of the dynamic response of submerged structure subjected to underwater explosion (UNDEX) is an important step for naval ship designers.The UNEDX experiment is the major research method to study the structure responses and failure modes.A prototype experiment is very close to the real experimental condition; however, it also has the defects of high expense, complicated process and poor repeatability.And the available information obtained from a prototype experiment is frequently limited for the instrument used.So, comparing with the prototype experiment, the model experiment is more available to the UNDEX researchers.
By high-speed rotation, the centrifugal force is exerted by centrifuge on the model, and the self-weight stress, deformation and failure mechanism of the model can match the prototype [1], so the centrifuges can not only solve many conventional rock and soil mechanics ,such as the explosive-induced soil ejected kinematic and crater morphology [2], and the variably saturated flow in soil and fractured rock [3], but also directly simulate the damage effect of structures subjected to kinds of dynamic loadings, such as earthquake, explosion and impact [4][5][6], and other engineering problems such as dam construction by directional blasting, mining, slope treatment and foundation treatment [7][8][9].Among them, the centrifugal model study on the dynamic responses and the protection of the structures under the impact of explosion loadings has become a hot topic of the research field of dynamics in recent years.Liu [10] carried out centrifuge experiment on the thin-walled aluminum tubes to investigate the response of transit tunnels in saturated soils under internal blast loading; De [11] simulated the effects of surface explosion on the underground tunnel with the centrifuge and gave potential mitigation measures through numerical analysis; De [12] also used centrifuge to build the physical model to study the effect of an underwater explosion on a buried tunnel and found that the strains in the tunnel increased with water level above the ground surfaces and then remained constant at a certain point ,which was related to the mass of explosives.Lu [13] conducted centrifugal model experiment on the brittle failure of a concrete dam subjected to underwater explosion and put forward the concept of the failure load intensity scale, and that the breaking strength of the brittle material scale is equal to the failure load intensity scale for the first time.
The interaction between the UNEDX shock loading and a structure is very complicated, and usually, the action time of shock wave is very short within a few microseconds while a bubble pulse can usually last for hundreds of microseconds [14].Since the responses of structure subjected to UNDEX shock wave and bubble loading are totally different, only one part of structure response is studied in some researches.However, many experiments have proved the importance to consider both the two factors during an UNDEX analysis [15,16].With the fast development of study method, numerical simulation has become another important and effective approach to study UNDEX problems expect experiment and theoretical analysis.Based on the acoustic-structure coupling technique, the ABAQUS UNDEX Analysis (AUA) method [17] can effectively solve complex UNDEX problems.The advantages of the AUA are that the shock wave and the bubble pulse pressure can be both simplified as the acoustic wave, which means that the two different shock loadings can be unified, and they can be also imposed on the structures by steps according to the experimental conditions.
A plate panel is a basic structural component of ships and submerged structures.In this paper, UNDEX experimental studies carried out on an air-backed rectangular steel plate with different charge weights, standoff distances and centrifugal forces were presented.Small charges and micro detonators were used.The response of the plate panel remained within the elastic range, and no permanent deformation occurred.The numerical analysis of the plate was performed with ABAQUS finite element(FE) code and the results are compared with the experimental data from reference [18] and [19].

Shock wave
During an underwater explosion, there are two main shock loadings, shock wave and bubble pulse.The shock energy delivering to the structure by the underwater explosion is closely related to the charge weight and standoff distance.On the shock wave, Cole had made lots of experimental researches.By analyzing the pressure-time cures of the specific test point, he put forward the empirical formula of UNEDX shock wave based on the similarity theory.The shock pressure at a given point has a sharp peak in time, followed by a decaying exponential function, which can be given by [20]: And the function of peak pressure (MPa), and the time decay constant (μs) are: where, ( ) is the pressure (MPa) at any instantaneous (μs), is the charge weight (kg), is the standoff distance (m), is the shortest radial distance, = ( − )/ is the time delay (μs), is the sound velocity in water (m/s), and , , and are the shock wave parameters of the explosive.

Bubble pulse
On the periodic motion of bubble pulse, the approximate relation of the emergence time of the first of bubble pulse (ms) and the maximum bubble radius (cm) are given by Cole as follows [20]: where is the water depth (m) of the explosive, is of equivalent water depth (m) of atmospheric pressure, which is about 10.33 m, and and are the motion parameters of the bubble.
Moreover, based on the assumption of linear compressible fluid, Geer and Hunter [21] derived the equations of motion of UNDEX bubble, through doubly asymptotic approximation technique (DAA), which simultaneously considers both dilation and translation during bubble oscillation.And the bubble motion can be expressed as: where and respectively represent density and corresponding sound velocity of inner gas product, is the impedance of inner gas and outer fluid, is the hydrostatic pressure, which includes the atmospheric pressure, is the gas product pressure, is the empirical flow drag parameter, and is the correction coefficient to match the experiment result.By giving the initial condition, the five characteristic parameters of bubble oscillation, namely, bubble radius , floating upward displacement , expansion velocity potential , translational potential of outer fluid and translational potential of gas product, can be calculated with Runge-Kutta method.
And the free-field pressure during the bubble oscillation phase can be expressed as [22]: where ( ) and are the second derivative of bubble volume and radius versus time, respectively.

Experiment setup
In this work, the experiments were performed on LXJ-4-450 centrifuge in IWHR laboratory, which has a maximal acceleration of 300 G, effective load of 1.5 t and maximal turning radius of 5.03 m.As Fig. 1 shown, the model box (1.2 m×0.7 m×0.9 m) was put in the basket (1.5 m×1.0 m×1.5 m) and excited centrifugal force along the cantilever direction.The model box was made with alloy plates fastened by bolts, and there was a glass observation window at one of side of the box.The inner walls of the box were coated with absorption materials to reduce the wave reflection.As indicated in Fig. 3, the steel plate of size 0.7 m×0.6 m×0.05 m with an exposed area of 0.6 m×0.65 m was fastened by cement mortar with thickness of 5 cm at the right, left and bottom side.A3 steel was selected as the target plate and the water-storage depth in the model box was 0.6 m.The explosives or micro detonators and the pressure sensor (ICP138) were located along the perpendicular line though the symmetry axis of the plate and at a same water depth.There was a high-speed camera (FASTCAM-ultima APX) with a frequency of 2000 fps and resolution of 1024×1024 put near the observation window to shoot the bubble pulse.The accelerometers (ICPM350B01) and strain gauges (TSJ120-4AA) were installed at the back of the plate (Fig. 4).All the dynamic signals were collected through data-acquisition system with 32 channels, which was put on the top of the shaft of centrifuge.The data-acquisition system was set delay trigger up to 20 ms to ensure obtain the complete measured data containing the initiating time of the charge.

Explosives
In order to ensure precision in small quantity, RDX spherical explosives and micro detonators were used in the experiments (Fig. 2).The RDX spherical explosives are composed by two separated hemispherical explosives, and a 10 cm long detonating fuse with a diameter of 2.6 mm and equivalent TNT weight of 48 mg/m, was detonated by a micro detonator (RDX equivalent 50 mg) to guarantee the reliable central initiation of the charge.Thin plastic sleeve made by 3D printing, hot pyrocondensation (HP) pipe and hydrophobic membrane were adapted to tightly connect and protect the whole detonation device.The density of explosive is 1.65 g/cm 3 , the detonation heat is 5200 J/g, and the detonation velocity is 8160 m/s.Three micro detonators were bounded together to simulate the effect of 150 mg detonator in this work.The explosive is converted to TNT equivalent by a coefficient of 1.58.

Determining shock parameters of shock wave
Seventeen tests, using different amounts of RDX charges and detonators, were performed in this work.The experimental conditions are shown in Table 1, and , , , and separately represent charge weight, centrifugal force, water depth, test distance between the charge and the plate, and standoff distance between the charge and the pressure sensor.Since some errors may occur during experiments due to underwater currents and irregular tension from the signal cables , all the effective standoff distances in each case were modified by the multiplication of the wave-propagating velocity (1498 m/s, 25℃) [21] .The Fig. 5 shows some of the time-history of pressures measured, from which it can be found that the peak pressure increases with the decrease of , and although there is some difference between the peak of S1 and S4 owing to the charge weight or test error, it can also prove that the centrifugal force have no effect on the shock wave, which can be also indicated in Table1.The reason is that the shock wave phase is too short, so the centrifugal (or gravity) effect can be ignored, and in addition, the peak pressure of S6 and S7 illustrate that the small change in water depth almost have no influence on the shock wave test of this work, so, the shock wave parameters can still be dealt with Eqs.(1) to Eqs. (3).
The relationship between peak pressure and weight-distance parameter ⁄ / can be determined through a bi-logarithm linear regression [23], and the regressive curve is = 1.6695 + 4.26415, where, = ln(P) and = ln( ⁄ / ).Fig. 6 shows the experimental data and regression curve.The peak pressure can be calculated by taking exponential of the curve, and the time decay constant can be obtained through a bi-logarithm linear regression on each pressure time-history curve.According to Eqs. ( 2

Centrifugal effect on bubble pulse
Before every test, a scale was put at the same position of the charge and shoot by the high-speed camera.Many data points of every bubble image were depicted according to the scale to get the mean radius of the bubble at a given time, and since the detonation process is very short when compared with the bubble oscillation, the moment of the emergence of the highlight area for the first time was approximately regarded as the beginning time of the bubble oscillation.
Owing to the centrifugal force generated by the high-speed rotation of centrifuge, the gravitational field of fluid and bubble buoyancy change a lot.According to empirical Eq. ( 4) and ( 5), the motion parameters and of bubble are only related to the water depth when the charge type and weight is determined.In this work, different centrifugal force were converted to water depth based on the equivalent conversion of hydrostatic pressure, and with the test data of bubble, the emergence time of the first of bubble pulse and the maximum bubble radius are fitted as follows: From Eqs. ( 12) and ( 13), and the bubble motion image (Fig. 9), it can be directly seen that and decrease with the increase of centrifugal force.And in order to explain this phenomenon in centrifugal model experiment, the radius-time and pressure-time curve of bubble are calculated through Eq. ( 6) to Eq. ( 9).Fig. 7 and Fig. 8 shows the comparison of Geer and Hunter model and experimental data.It can be seen that the change trend of the measured and the theoretically calculated bubble motion are compared well, and peak bubble loadings decline when centrifugal force rising, ignoring the initial charge weight difference.In addition, the theoretical loading curves are more symmetric as a result of assumption of symmetric bubble in Geer and Hunter model.Zamyshlyayev [24] also established a shock loading model dividing the whole explosive process into five phases.On the basis of Zam's theory, the hydrostatic pressure higher, the peak bubble loading higher, which seems to be opposite to our result.And it is because of the assumption of still bubble in the first motion cycle in his model.However, when the centrifugal force increases more, the bubble oscillation become more obvious with the expansion and translational potential changing more sharply, so, the whole bubble displacement in vertical direction increase a lot, and the first bubble motion cycle shorten for the increasing drag force, which can be calculated by the Geer and Hunter model.In each case of this work, the pressure sensor was set up at a specific position, but the bubble center is a moving point source and the distance from the bubble center to the sensor lengthen with the rising centrifugal force, so, although the increase of initial hydrostatic pressure raises the bubble potential energy when bubble shrinking to the minimum radius, the measured loading peak decrease.As a result, in this work, the pressure field caused by bubble oscillation must take consideration of the vertical moving of bubble for the big change of centrifugal force.

Numerical simulation
Based on acoustic-structure coupling technique, the dynamic response of air-backed steel plate subjected to UNDEX explosion considering the centrifugal effect was simulated with explicit analysis on program ABAQUS 16.4.

Analytical model
On account of the fact that the inherent nature of the materials is mainly related to the electromagnetic forces, so the basic material properties in the centrifugal model can be regarded as to be the same with the conventional model [25].The numerical model was established under the same experimental condition (Fig. 10), which was composed by three parts, namely, water, steel plate and cement mortar.Since the response of cement mortar is not the concern of this study, it was treated as discrete rigid body considering its specific shape and meshed with 22650 R3D4 4-node 3-D bilinear rigid quadrilateral elements [26].In the simulation, the steel plate adapts the linear elastic model and the material constant used were summarized in Table 2.The plate part was modeled by 21000 C3D8R 8-node linear brick, reduced integration, hourglass control, solid elements, which is not easy to lead to volume locking and save computing cost, and the steel plate was meshed with five elements in the thickness direction to avoid hourglass [27].The external fluid water was approximate to a regular cuboid, so it was modeled using acoustic medium element AC3D8R, and the material property constant are also described in Table 2. AC3D8R element is a 8-node linear acoustic brick, reduced integration, hourglass control, solid element, which does not have displacement degrees of freedom [28] (Fig. 10).The explicit time integration method is a conditional stable integration and is closely related to the element characteristic length.Since the wavelength of shockwave increase with the frequency decreasing, the maximum permitted element length of fluid element , can be determined by calculating the upper range of shockwave frequency .The general principal is < /( .), where is the wavelength of the shock wave within the element, and the recommended value is ≥ 6. [29].And the water domain was established with 444600 elements in the model.

Boundary condition
The interface between steel plate, cement mortar and water respectively adapt acousticstructure coupling interaction, the interaction between the plate and cement mortar was set as surface to surface contact, and the interfaces between steel plate, cement mortar and the inner walls of model box were separately set as clamped boundary.Although the fluid boundary may lead to the scattering and reflection, which will reinforce or weaken the shock loading, the influence of inner walls and free water surface were ignored in the simulation, considering the small weight of charges and test results.So, the interfaces between fluid and inner walls were set as non-reflective boundary and the reflective effect of free water surface was also not considered.

Input shock loading
Usually, if the bubble oscillates at a very close distance near the surface of the structure or free surface, either the bubble motion cycle or the maximum bubble radius will be much affected , the bubble will be rejected at the expansion stage, and be drawn at the collapse stage, and a water jet will form ,which can be explained by the Bjerkness effect [31].While, in most cases of this work, the test distance and water depth was much bigger than the maximum bubble radius (Fig. 7), namely, / > 2 and / > 2, so, the experiment condition can be approximately treated as mid-field or far-field underwater explosion [32].Mu [33] applied MSC.Dytran to study the characteristics of UNDEX bubble under different boundary conditions, and found that when / > 2, the free surface will not influence the bubble oscillation, so the Geer and Hunter model can be applied to calculate the bubble motion parameters, which also has been proved in section 3.4 and reference [18].From the comparison of the bubble motion tested and calculated, the boundary condition do not affect the bubble oscillation for most experiment cases a lot, so the general input loading-time curves of numerical model were approximately calculated from Eqs. ( 9) and (10).Fig. 11 compares the experimental and theoretical data of shock loading of S9, and considering the obvious bubble flotation, the theoretical vertical displacement-time curve is applied on the source point as a displacement constraint.When the bubble floats to the free water surface, namely the bubble displacement are equal to water depth, the loading input process in FEM analysis will automatically terminate at end time .Fig. 12 shows experimental and theoretical displacement of the bubble of S9, and only the data segment in the black dashed are valid.

Alpha damping effect on the response
Different from the static analysis, the dynamic analysis must consider the effect of structure deformation and motion, such as damping, which refers to the energy dissipation in the actual motion of the structure.Damping exists in most engineering problems, and even though it is only the approximated simulation of the energy absorption features of the structure, rather than the simulation of the physical mechanism bring about this effect, it is very important for obtaining accurate calculation results [34].And the introduction of a reasonable damping value to the numerical simulation can effectively restrict numerical oscillation and increase the physical damping of the whole calculation system to make the numerical results better reflect the real physical phenomenon.
Since the interaction between the plate and UNDEX shock loading is the main concern of this work, the influence of mass ratio damping Alpha is mainly investigated.The Alpha factor define a damping contribution proportionate to element mass matrix, which is related to the absolute nodal velocity of the numerical model, and the effect of the definition of Alpha damping can be compared to the movement of the structure going through a viscous liquid [35], therefore, any nodal movement in the model can lead to the damping force.In the ABAQUS/Explicit analysis, the oscillation of low-frequency response can be effectively adjusted by setting a reasonable Alpha damping [26].Fig. 13 and Fig. 14 show the effect of Alpha damping on the acceleration and strain response of the steel plate, from which it can be seen that the effect is significant.Both the peak of high-frequency responses caused by shock wave impact and low-frequency responses induced by bubble pulse has a little decline with Alpha damping increasing.However, comparing with high-frequency responses, low-frequency responses are more sensitive to the change of Alpha damping.When Alpha damping is very small, the low-frequency responses of acceleration and strain tend to be obvious oscillations, and the acceleration even cannot be distinguished (Fig. 13), while the high-frequency responses almost keep a stable mode of vibration.And by a great deal of trail calculation, the numerical results agree well with the test results when Alpha damping is close to 1000.

Acceleration
Comparison of numerical and tested peak acceleration at various measurement points on the target plate is shown in Table 3. Combined with Fig. 13, it can be seen that the difference between the peak of high-frequency and low-frequency response is very obvious, which is corresponding to the peak loadings and reaches to tens of times.The peak acceleration induced by shock wave is relatively steady, and the peak acceleration caused by bubble pulse decrease gradually for rising centrifugal forces.Owing to the test error, some test and numerical data does not fit very well, and when the standoff distance is very small (15 cm at S10 and 5 cm at S11), the errors between experimental and numerical results all become very big, and it is because that when the charge is very close to the target plate, the formation of bubble jet is greatly affected, which impact the plate like a water hammer [36], and the shock wave and bubble loading will significantly interrelate, so the approximate treatment for free field of a bubble pulse in is not very appropriate again.

Velocity
Although the particle velocity of the steel plate was not measured in this work, it can be obtained through integration the measured accelerations.And since the acceleration peaks caused by shock wave loading are much bigger than that of the bubble loading, only the velocity responses during shock wave phase are discussed.The period of the integration was from the arrival of the shock front to local cavitation occurred [37].In the light of the complex constraint boundary of steel plate and initiation condition of charge, some simplification were adopted in the velocity analysis, including that the transmission of shock wave into the plate are ignored, and that the included angle between the straight line through the source point and a given particle and the normal through the source point is equal to the incident and reflection angle of shock wave, and lastly, the upper surface of the plate part exposed to the air is handled as a clamped boundary.Rajendran [38] have derived the particle velocity of plate under normal incident shock wave, and an angle parameter were introduced in this paper to calculate the velocity considering the impact at any angle of incidence.On the basis of boundary continuity, the velocity of the particle on the incident side is calculated as: And according to Talyor's theory [39], the total pressure in normal direction, which acted on the air-backed plate can be expressed as: where ( ) is the normal velocity of the particle on the plate-water boundary, ( ) and ( ) are the incident and reflected shock wave, respectively, is the incident angle, which is equal to reflection angle, and ( ) is total wall pressure of the plate.The particle velocity can be obtained by substituting Eq. (1) into Eq.( 13) and Eq. ( 14): where is plate mass per unit area, and is the characteristic mass ratio.When ( ) drop to zero, the plate and fluid will be separate and cavitation occur [37], the acceleration process of plate owing to shock wave stop and the normal velocity of plate particle reach to the maximum, given by: The time , when the total wall pressure drops to zero, can be calculated as: And is covert to the theoretical result of Rajendran [38] at = 0.In UNDEX researches, shock factor is always regarded as an effective characteristic parameter indicating shock energy of UNDEX shock wave, and a shock factor is adopted in this paper to evaluate the velocity response, which is defined as [40]: Fig. 15 shows the maximum velocity obtained by integrating measured acceleration, approximated by Eq. ( 18) with peak pressure, and the velocity calculated by the finite element method.All scattered points of peak velocity are linear fitted with the shock factor from the original point, and it can be found that the linear correlation between the peak velocity and the shock factor is much better at the position, where it is much closer to the area the shock wave vertically impact.In particular, when the standoff distance reduces to a certain extent (in Fig. 15(c)), the peak velocity calculated by Eq. ( 18) begins to decrease with shock factor increasing, which is not consistent with the experimental and numerical results, and it is because of the approximation of the incident angle.When the target position is far from the point shock wave vertically incidents upon, the approximated oblique incidence angle is too big and not reasonable again.

Strain
In this work, all the strain gauges were firmly adhered to the steel plate, and the experimental and numerical peak strain (S1, S2, S4 and S5) at various measurement points were indicated in Fig. 16 and Fig. 17.The peak strain caused by shock wave in experiment are relatively decentralized for test errors, and most of the experimental peak strains induced by bubble pulse well match the numerical peaks.From the numerical results, it is can be seen that the peak strains induced by shock wave remain stable when centrifugal force changing, the peak strains due to bubble pulse also fall with the centrifugal force increasing, just like the peak accelerations, however, the peak strains induced by the two factors are very close, which is totally different from the acceleration responses.a) Experimental [19] b) Numerical Fig. 16.Peak strain caused by UNDEX shock wave According to the modal analysis theory on distributed strain measurement [41], the strain response is more similar to the displacement response than the acceleration response, and it is more sensitive to the low-frequency response from the angle of time-frequency domain, which are the responses caused by bubble pulse in this work.Although the pressure of bubble pulse is much lower than that of the shock wave, the duration time of bubble pulse is much higher than that of shock wave, so, even though the high-frequency responses of the strain and acceleration are similar, the change amplitudes of the low-frequency responses of the strain are much higher in contrast to the acceleration responses for the accumulation of time domain.In other words, the strain responses are more sensitive to time accumulation than the acceleration responses.Moreover, from the angle of shock energy, although there is a big difference between the peak pressure of the shock wave and the bubble pulse, the shock energy carried by shock wave and bubble are close [20], so, the difference between the acceleration and strain response can illustrates that the strain response is tightly related to shock energy.And actually, Rajendran has [42] deduced the maximum equivalent strain at one position of the plate in the light of air-backed steel plate impacted only by shock wave, which can be expressed as a function related to peak pressure and decay time, and it can also indirectly elucidate this phenomenon.

Conclusions
In this study, the UNEDX experiments were carried out in a centrifuge to expound the dynamic response of an air-backed steel plate impacted by shock loadings.The bubble oscillation considering centrifugal effect were predicted with Geer and Hunter model, and the experimental response were compared with Finite element analysis under the same plate and shock loadings.According to the experimental and numerical results, we conclude the following: 1) The centrifugal force has no effect on UNDEX shock wave for the short function time, however, with the increase of centrifugal force, the initial hydrostatic pressure increase, the bubble motion cycle and maximum bubble radius decrease, the floating upward of the bubble become more obvious, and the peak pressure of bubble pulse decrease.
2) Both the high-frequency and low-frequency response have a slight decline as the alpha damping rising; the vibration model of high-frequency response remains stable when centrifugal force changes and the low-frequency response tend to be an obvious oscillation waveform as the alpha damping is too low.
3) The difference between the peak accelerations respectively induced by shock wave and bubble pulse are very large ,tending to be tens of times, while, the peak strains caused by the two factors are very close, which illustrates that the strain response is related to shock energy and more sensitive to time accumulation; during the shock wave impact ,both the acceleration and strain response are not influenced by centrifugal force and the two kinds of response peaks decline with centrifugal force increasing.
4) The linear correlation between the peak velocity and the shock factor is much better at the position, where it is closer to the area shock wave vertically impact.

3 Fig. 15 .
Relationships between the peak velocity and shock factor

.
Peak strain caused by bubble pulseJianyu Wu is the principal person to complete this paper, Long and Zhong gave important advice to this paper, and Song and Hu participated in the experiment of this paper.

Table 3 .
Peak acceleration ( ) at various measurement points on the target plate