Shock and Vibration
Volume 2015 (2015), Article ID 741986, 12 pages
http://dx.doi.org/10.1155/2015/741986

Research Article

In-Structure Shock of Surface Structures: A Theoretical Approach

Hongyuan Zhou¹ and Xiaojuan Wang²

¹School of Civil and Environmental Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798
²Department of Civil and Environmental Engineering, National University of Singapore Singapore 117576

Received 14 May 2015; Revised 6 July 2015; Accepted 7 July 2015

Academic Editor: Sakdirat Kaewunruen

Copyright © 2015 Hongyuan Zhou and Xiaojuan Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Metal foam cladding emerges as a promising measure for structure protection against close range blast. With proper cladding design, a considerable amount of energy can be absorbed while the load transmitted to the protected structure can be lowered to a controllable level, leading to less damage of the protected structure. However, it is not clear whether the in-structure shock can be mitigated by attaching foam cladding. The present study theoretically investigates the response of a structural member subjected to a close range blast, for both scenarios with and without attaching foam cladding. Subsequently, shock response spectra are employed to assess the in-structure shock of both situations. It is found that applying foam cladding does not effectively mitigate in-structure shock of surface structures.

1. Introduction

The protection of buildings and structures against accidental and terrorist attacks is of interest to scientific and engineering community. A traditional method to mitigate a blast load is to strengthen the structure itself. With thicker structural members and stronger reinforced materials, both the structural stiffness and strength are increased, leading to small deflection subjected to the same blast load. However, the drawbacks of strengthening do exist. With higher rigidity, the incident blast load is violently amplified, resulting in a higher load exerted on structures. Further, the construction and retrofit of such strengthened structures are expensive and time-consuming. Recently, a new protection philosophy emerged: an additional cladding with metal foam core is attached to the target structure, which lowers the load transferred to the protected structure and meanwhile absorbs a considerable amount of energy, when subjected to the same blast load [1]. Sandwich panel behaviors are investigated and summarized [2], while energy absorption of various cellular structures and materials is examined [3]. Field test is conducted for aluminum foam core cladding protected concrete panels, and it is found that the protection is effective [4]. To ensure controllable constant load transmitted to the protected structure, cellular solid with relative density less than 0.2 should be applied [5]. When subjected to the same blast, the load transmitted to the protected structure is lowered while the cladding is crushed, which absorbs a large amount of energy. After that, the damage cladding is easily and quickly replaced with a new one; thus both resiliency and cost efficiency are achieved [6].

Generally, such a protective cladding is made of a face plate and a layer of metal foam or honeycomb as core. Subjected to rapid movement of the face plate induced by blasts, the foam adjacent to the face plate is crushed first; then the densification continues as a plastic shock wave [7]. In the upstream of the shock wave, metal foam remains almost undeformed and the load exerted on the protected structures is the foam plateau stress, which is adjustable by varying the foam relative density and considerably lower than the blast load.

While no doubt the structure protection in terms of strength is of importance, the mitigation effect in terms of vibration can never be ignored, especially when vibration-sensitive device is mounted or attached to the protected structure. From the viewpoint of structural engineering, the strength aspect of protection is a matter of structural limit analysis while the vibration counterpart is an issue of serviceability. In-structure shock is a kind of strong vibration induced by an external blast or ballistic load not damaging or lightly damaging the main structure but causing devices within the structure to fail [8]. In-structure shock of underground structures with rigid body motion taken into account is investigated [9]. Double-layer floor is proposed to successfully mitigate vertical shock of underground structures [10]. Small scale test is conducted to validate the theoretical prediction of in-structure shock of underground structures [11]. Hybrid mounting system is tested [12]. In early years, in-structure shock of surface structures is of minor significance since whether the structure can survive the explosion is the major concern. However, with the development of blast mitigation technology such as consolidation and attaching metal foam cladding, in some occasions, major damages of the protected structures can be avoided. Then the in-structure shock emerges as a problem of substantial significance. In the present study, in-structure shock assessment is performed for a surface structure with and without metal foam cladding subjected to a close range explosion, shown in Figure 1 [6]. Energy method can be used to determine the crater extent and depth. Then the crushing time history of the face plate (and the densified foam) is calculated, taken as the loading time for the protected structure. Subsequently, response of idealized structural member with cladding attached is obtained and compared to that without cladding protection. Finally, the in-structure shock of the structural member with and without cladding is examined.

Figure 1: In-structure shock of surface structure with metal foam cladding subjected to close range blast.

2. Response of Structural Member with Metal Foam Cladding Subjected to a Blast

2.1. Time History of Foam Crushing

To obtain the duration of the load transmitted to the protected structural member from foam densification, crushing of the metal foam cladding subjected to a close range airburst is analyzed first. The cladding attached to the protected structure is a layer of metal foam covered by a thin solid face plate. When subjected to a blast load induced by a close range detonation, the nonuniform pressure will create a bulge, which is deepest in the point originally nearest to the explosion before deformation and its depth decreases towards the bulge periphery to zero. As shown in Figure 1, only the foam under the round bulge is crushed while the other foam remains undeformed [6].

The face plate of a typical blast mitigation cladding is a thin metallic plate. Bending does not contribute to the energy dissipation when the transverse deflection is greater than half of the plate thickness [13]. Subsequently, for sufficiently large deformation of a plate, the influence of the bending is ignored and the response is solely governed by membrane force. In the present study, the face plate is reasonably assumed as a rigid-plastic membrane under a close range explosion, as the face plate deflection is significantly larger than half of its thickness. Take the face plate and the densified form to analyze. Compared to the face plate, the stretching energy dissipation contribution of the metal foam core, both densified and undensified, is relatively small and can be neglected [14]. By equating the work rate by external load, inertia force, and foam resistance to the energy dissipation rate by the membrane stretching, the governing equation for the system consisting of the face plate and densified foam in the blast duration is

where is the nonuniform blast load time history considering fluid-structure interaction; , are the density and thickness of the face plate; , , and are the density, plateau stress, and densification strain of the foam, respectively; is the deflection of the face plate; and and are internal forces (in fact the internal force per unit length, according to the convention in plate theory) of the face plate in radial and circumferential direction. Recall [15]:

is the radius from the bulge center and is the yield stress of the face plate. Equation (1) is recast as

where is the bulge extent [6], where is the peak reflected pressure on face plate defined in (6) later. The first term in the square brackets is the blast load applied on the face sheet, in fact the external driving force. The second term is the inertial force of the system consisting of the face plate and the fully densified foam part, in which is the mass per unit area of the face plate and is the mass per unit area of the densified foam. The third term in the square brackets is the resistance from the dynamic foam crushing, including the foam plateau stress (the quasi-static part) and the dynamic enhancement . Decouple the space and time of the face plate deflection time history:

where is the bulge shape, is the characteristic dimension of the bugle, and . Then the velocity and acceleration of the face sheet are

In the present study, the bell-shaped blast load is approximated as

where is the peak reflected pressure on the face plate and is the blast duration. is the characteristic size of the loading region, the same as the standoff distance between the explosion and face plate. At the load decreases to 1% of the peak value and is neglected. This blast load distribution assumption agrees well with the measured load in test [6].

With (5), (3) is recast as

Substitute (4) and (6) into (7) and integrate it:

In short term,

where

It is worth noting that the coefficients in (9) are simplified by neglecting the trivially small items using the relation . Equation (9) is a nonlinear ordinary differential equation and its solution cannot be analytically obtained due to the presence of the square item. Four-order Runge-Kutta method is adopted to numerically solve the governing equation for .

If the blast load is relatively weak compared to the resistance from the face plate and metal foam, the crushing may cease in the blast duration. However, for a typical close range explosion, the blast intensity is substantially higher than the resistance, where a relatively larger amount of energy is imparted into the cladding resulting in a crushing period longer than the blast duration. Generally, the blast load duration induced by a close range airburst varies in the order of 0.1 milliseconds or less and the structural response may be several times longer. After the blast load, there is no external loading on the cladding. Offset the time by , let , and the governing equation after blast is

The initial condition of the second phase is the terminal response of the blast duration to ensure the continuity of the displacement and velocity:

Again, combine the governing equation and the initial condition for the second phase, and the Runge-Kutta method is employed to solve (11).

When the face plate velocity decreases to zero at , the foam crushing ceases. is determined by intersecting the generalized velocity, , with time axis. Finally the entire foam crushing duration due to the close range explosion is

where is the instant of crushing stops with respect to time starting from .

It is worth noting from coefficient of (9) and (11) that the crushing time depends not only on the blast load and cladding parameters but also on the final bulge depth, which in fact is also a function of the blast and cladding parameters. However, the numerical predictions of (9) and (11) indicate that the crushing time of the cladding is independent of the bulge depth, as shown in Figure 2: when the blast load and cladding parameters are the same, the crushing duration is the same in spite of the different bulge depths arbitrarily assigned. The only purpose of the crushing time history analysis is to determine the total crushing time. Therefore, the bulge depth is not necessarily in its true value and set as unity (although not practical in engineering practice) in the present study.

Figure 2: The total crushing time does not depend on the bulge depth.

According to the shock theory for foam [2, 3, 7], the crushing process in fact is a plastic wave propagation process, shown in Figure 3. Downstream of the shock front, the foam is densified with high stress and strain while, upstream of the shock front, the uncrushed foam is in elastic compression with foam plateau stress [2]. Therefore no matter how strong the load is applied on the face plate, the load transmitted to the protected structure is the foam plateau stress, provided that the foam is not totally crushed. Total crush of the foam layer may result in a load even higher than the blast load on face plate and should be avoided.

Figure 3: Foam crushing in progressive collapse mode and load transmitted to the protected structure.

Consider that a structural member of a building is protected by attaching a metal foam cladding, subjected to a close range spherical detonation. The member is loaded by the foam in a round region due to the axis symmetry of the close spherical explosion, indicated in Figure 1. The boundary of the structural member is connected to adjacent slabs, neither simply supported nor fixed. In the present study, it is considered as simply supported at four edges, leading to a conservative prediction, safe in design.

2.2. Equivalent SDOF System for a Simply Supported Rectangular Structural Member

Subjected to a dynamic load, analytical solution for continuous system is generally difficult and cumbersome to obtain, especially when the load is nonuniformly distributed on a part of the structure. To analyze a continuous system under dynamic load in a relatively simple way with adequate accuracy, a single-degree-of-freedom (SDOF) equivalent method for continuous systems is frequently adopted [16]. In the idealization, the time scale is not altered; therefore the response time history of the idealized system is exactly the same as that of the real structure. Before idealization, the deformation profile of the structural member must be predefined. Various approximate deflection modes can be chosen for the simply supported rectangular plate; in the present study, the fundamental deflection mode is adopted. If higher analysis accuracy is required, higher order deformation modes can be also incorporated, however, with higher analysis cost.

According to Biggs [16], the equivalent mass of the rectangular plate is calculated:

where , , , and are the two edge dimensions, density, and thickness of the structural member, respectively. is the approximate deformation profile, taken as , where the origin of the coordinates is taken as one of the plate corners. The transformation factor for mass is

The subscripts and denote “equivalent” and “real,” respectively. In calculating the equivalent load, double cosine function, instead of double sine function, is employed with coordinate shift while keeping the idealization nature exactly unchanged:

where is the bulge extent. The corresponding real load is

However, (16) cannot be analytically integrated and thus evaluated numerically. The transformation factor for this load is

The transformation factor for stiffness is exactly the same as that for load [16]; therefore the equivalent stiffness is

is the stiffness of the rectangular plate with four edges simply supported which is calculated as follows. Generally, stiffness of a structural member is defined as the force necessary to cause a unit deflection at the point of interest or the ratio of force applied to the deflection caused. For a continuous structural member such as a plate in the present study, to determine the stiffness, the area of loading must be considered. Specifically, a concentrated load on the plate center causes larger deflection at this point, resulting in smaller plate stiffness; when the concentrated force is uniformly distrusted over the whole plate surface with total force exactly the same, the deflection at the plate center is relatively smaller; thus the plate appears stiffer. Therefore, within the elastic limit, although the stiffness is not a function of loading intensity, it is a function of the loading pattern or area in terms of the above-mentioned stiffness definition. Then the plate stiffness with central round loading area in the present study (the most critical case is loaded in the center) should be between the two cases mentioned above. The deflection for a simply supported rectangular plate subjected to uniform load on a circle region around the center is

where

is the bending rigidity of the plate. It is worth noting that the series in (20) converge quickly and only first five terms can achieve adequate accuracy [15]. The real stiffness of the structural member is

Until now, the equivalent mass, load, and stiffness are calculated and the equivalent natural frequency of the idealized SDOF system is

2.3. Response of the Structural Member with Metal Foam Cladding

With metal foam cladding attached, the load exerted on the structural member is the foam plateau stress and the loading duration is the period of cladding crushing. The governing equation for the structural member within blast duration is

The damping ratio is taken as 5% of the critical, typical for surface structures [17]. With initial conditions of zero displacement and velocity, the displacement and acceleration time history of the structural member center are

where

In the postloading phase, the structural member undergoes free vibration until the motion is damped out, which also contributes to the in-structure shock. Offset time by the cladding crushing duration: . The response time history in the free vibration phase is

where

where , are the initial conditions for the free vibration phase ensuing the continuity of the displacement and velocity.

3. Response of Structural Member without Cladding Subjected to a Blast

Without cladding, for the same detonation, the scaled distance is larger due to the absence of the cladding thickness, resulting in a smaller loading intensity and a larger loading area, compared to the load on the face plate.

The real load on the structural member (within the bulge) is

With the same idealization method for the load on structural member as briefed in the preceding section, the equivalent load is

The transformation factor for load is

It should be noted that, in the present study, only the load within the bulge is considered, for both the cases with and without foam cladding attached. In so doing, the response can be compared without considering the load outside the bulge, based on the principle of superposition of elastic systems (the structural member is assumed elastic).

Assume that the structural member is statically loaded with the peak load of time zero, and the static deflection at the center is

where

Then the stiffness is determined:

Subsequently, the natural frequency of the equivalent system is

Although a typical damping does not have significant influence on the structural response within the short blast duration, it has significant influence on the free vibration phase [17]. A typical damping of 5% critical in structural response [17] is used in both phases for consistence. The governing equation for the structural member within the blast duration is

With zero initial conditions, the structural response at the center is

where

It is worth noting that the loading intensity is smaller than that on the face plate due to the larger standoff distance when the cladding is not attached. The response in the postblast phase is a damped vibration:

where , are the initial conditions for the postblast phase ensuring the continuity of the displacement and velocity between the two phases.

4. In-Structure Shock Assessment with a Case Study

A typical event is considered as a case study. Assume a spherical explosive with TNT equivalent charge of 0.157 kg detonated 0.3 m away from the face plate, producing a blast with duration of 0.1 ms and a peak pressure of 8 MPa (directly read from the chart in [18]). The face plate is made of aluminum with thickness 1 mm, while aluminum foam with relative density of 5% is used as cladding core, having a plateau stress of 0.67 MPa. The thickness of the cladding is assumed as 5 cm; thus if the cladding is absent, the standoff distance from the explosion to the structural member is 0.35 m. With the same explosive charge and larger standoff distance, the peak load on the structural member reduced to 5.1 MPa in the present study and the loading area is larger with radius 1.05 m. However, the blast duration is assumed the same as in the case with cladding, 0.1 ms.

By attaching the metal foam cladding, the load experienced by the structure is shifted: without the cladding, the load is a narrow and high triangular pulse while, with cladding, it becomes a triangle with a lower intensity and a longer duration being the foam crushing time. Figure 4 illustrates the blast load applied on the face plate and the load transmitted to the protected structures (the load is defined as resultant force). Although both loads are triangles with different peak values and durations, the mechanisms are different: without cladding, the loading intensity decreases with loading area unchanged while, with cladding, the loading area decreases with the loading intensity unchanged being the foam plateau stress.

Figure 4: Load on structure comparison with and without cladding.

The structural member is a simply supported concrete plate of 2 m by 3 m, with a thickness of 0.2 m. Young’s modulus and Poisson’s ratio are 30 GPa and 0.2, respectively. From (22) and (35), the equivalent circular natural frequencies of the two idealized SDOF systems in the case study are similar: one is 317 rad/s and the other is 318 rad/s. The equivalent mass is exactly the same for both idealizations. For the structural member without cladding, due to the larger standoff distance, the loading area is larger but the bell-shaped loading concentrates in the central area while, for the same structural member with cladding, the loading area is smaller, but the foam crushing induced loading uniformly distributes in this area. Regarding the definition of the effective stiffness, with equivalent stiffness almost the same, the equivalent natural frequencies of the two SDOF systems are similar.

In the case study, when the peak load of the close range airburst is mediate, for instance, 8 MPa, the mitigation effect of the cladding is negligible, as shown in Figure 5. When the peak load increases to a relatively high value by increasing the explosive charge, 20 MPa, the mitigation effect is still not remarkable, indicated in Figure 6. Further, parametric study shows that the foam crushing duration increases with the peak blast load, shown in Figure 7. In addition, the ratio of foam crushing duration to blast duration decreases with increasing blast duration, shown in Figure 8.

(a)

(b)

(c)

Figure 5: Structural member center response comparison with and without cladding. The blast duration is 0.1 ms and peak pressure is 8 MPa on cladding and 5.12 MPa on structure without cladding: (a) displacements; (b) accelerations with initial stage truncated; (c) accelerations in the initial stage.

(a)

(b)

(c)

Figure 6: Structural member center response comparison with and without cladding. The blast duration is 0.1 ms and peak pressure is 20 MPa on cladding and 12.6 MPa on structure without cladding: (a) displacements; (b) accelerations with initial stage truncated; (c) accelerations in the initial stage.

Figure 7: Relationship between the peak pressure and the crushing duration.

Figure 8: Ratio of the crushing time to blast duration.

The responses of a structural member with and without metal foam cladding are analyzed, which are in fact the excitation for the equipment or devices attached to this member. A device, with its support, can be simplified as an SDOF system consisting of a spring, a damper, and a mass. The excitation for this SDOF system is similar to the one frequently used in earthquake engineering, where the time history of base acceleration is the excitation to the attached systems. The response of the SDOF system depends on the characteristics of the excitation acceleration time history, natural frequency, and damping ratio (conventionally taken as zero for short duration load [17]) of the system. Then for a specific device with support, with zero initial condition, the responses such as displacement, velocity, and acceleration time histories can be determined. Subsequently, the maximum displacement, velocity, and acceleration are obtained. If the peak acceleration of certain device is greater than its corresponding tolerance, damage occurs [18].

The mass and stiffness characteristics of a specific device, with its mounting, can be delineated by one parameter: natural frequency or natural period. Therefore an SDOF system with varying natural frequency can represent a series of devices with supports. Shock response spectra are plots of the maximum responses of a series of SDOF systems subjected to the same excitation. When the excitation is prescribed, the maximum displacement, velocity, and acceleration can be read from the tripartite plot in different directions. Figure 9 indicates that when subjected to a close range detonation with peak loading pressure 8 MPa, the in-structure shock experienced by the devices mounted to the structural member is almost the same, regardless whether the metal foam cladding (of thickness 5 cm) is attached or not. When the peak loading is increased to 20 MPa while other parameters are unchanged, the shock level of the device is slightly (not remarkable) smaller for cladding-applied structural member compared to that without cladding, as shown in Figure 10. To sum up, attaching foam cladding to structural member does not mitigate the in-structure shock of surface structures.

Figure 9: Shock response spectra with cladding subjected to a blast load with peak pressure 8 MPa on cladding and 5.122 MPa on structure, with both durations being 0.1 ms.

Figure 10: Shock response spectra with cladding subjected to a blast load with peak pressure 20 MPa on cladding and 12.6 MPa on structure, with both durations being 0.1 ms.

5. Conclusions

A procedure for evaluating the in-structure shock of surface structures is proposed. By attaching a metal foam cladding to a structural member, the loading exerted on the member is altered from a pulse with a higher peak and a shorter duration to a pulse with a lower peak and a longer duration, although the pulses are both triangular. Due to the similar impulses imparted to the structural member and the relatively short loading durations, although elongated in the case with cladding, compared to the natural period of the member, the responses including the displacement and acceleration time histories are similar with and without the cladding. Attaching foam cladding in structural protection against blast does not mitigate in-structure shock.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

P. Zhang, Y. Cheng, and J. Liu, “Numerical analysis of dynamic response of corrugated core sandwich panels subjected to near-field air blast loading,” Shock and Vibration, vol. 2014, Article ID 180674, 16 pages, 2014.

M. F. Ashby, A. Evans, N. A. Fleck, L. G. Gibson, J. W. Hutchinson, and H. N. G. Wadley, Metal Foam: A Design Guide, Butterworth-Heinemann, 2000.

G. X. Lu and T. X. Yu, Energy Absorption of Structures and Materials, CRC Press, Boca Raton, Fla, USA, 2003.

C. Q. Wu, L. Huang, and D. J. Oehlers, “Blast testing of aluminum foam-protected reinforced concrete slabs,” Journal of Performance of Constructed Facilities, vol. 25, no. 5, pp. 464–474, 2011.

L. J. Gibson and M. F. Ashby, Cellular Solids: Structures and Properties, Cambridge University Press, Cambridge, UK, 2nd edition, 1997.

H. Y. Zhou, G. W. Ma, J. D. Li, and Z. Y. Zhao, “Design of metal foam cladding subjected to close-range blast,” Journal of Performance of Constructed Facilities, Article ID 04014110, 2014.

S. R. Reid and C. Peng, “Dynamic uniaxial crushing of wood,” International Journal of Impact Engineering, vol. 19, no. 5-6, pp. 531–570, 1997.

G. W. Ma, H. Y. Zhou, Y. Lu, and K. Chong, “In-structure shock of underground structures: a theoretical approach,” Engineering Structures, vol. 32, no. 12, pp. 3836–3844, 2010.

G. W. Ma, H. Y. Zhou, and K. Chong, “In-structure shock assessment of underground structures with consideration of rigid body motion,” Journal of Engineering Mechanics, vol. 137, no. 12, pp. 797–806, 2011.

H. Zhou and G. Ma, “Double-layer floor to mitigate in-structure shock of underground structures: a conceptual design,” Engineering Structures, vol. 35, pp. 314–321, 2012.

H. Y. Zhou, M. Beppu, G. W. Ma, and Z. Y. Zhao, “In-structure shock of underground structures: a revisit with experimental investigation,” Engineering Structures, vol. 56, pp. 1620–1630, 2013.

S. J. Moon, J. S. Kwak, J. H. Chung et al., “A study on the hybrid mount against vibration and shock for naval ships,” Shock and Vibration, vol. 17, no. 3, pp. 269–283, 2010.

N. Jones, Structural Impact, Cambridge University Press, Cambridge, UK, 1989.

F. Zhu, Z. H. Wang, G. X. Lu, and L. M. Zhao, “Analytical investigation and optimal design of sandwich panels subjected to shock loading,” Materials & Design, vol. 30, no. 1, pp. 91–100, 2009.

S. Timoshenko and S. Woinowsky-Krieger, Theory of Plates and Shells, McGraw-Hill, New York, NY, USA, 2nd edition, 1959.

M. A. Biggs, Introduction to Structural Dynamics, McGraw-Hill, New York, NY, USA, 1964.

R. W. Clough and J. Penzien, Dynamics of Structures, McGraw-Hill, New York, NY, USA, 1993.

U.S. Army Corps of Engineers, Naval Facilities Engineering Command, and and Air Force Civil Engineer Support Agency, “Structures to resist the effects of accidental explosions,” UFC 3-340-02, 2008.