Analysis of the pressure wave parameters caused by TNT underwater explosion forced on the hull of minehunter

In September 2013 Inspectorate of Polish Navy Armament signed a contract to build new ships for combating and setting mines codenamed “Kormoran II”. She has to be adapted to operate in the Polish economic zone, as well as in tactical operations in the Baltic and the North Sea. Basic assumptions tactical – technical shown below [1]:

• Displacement: max: 850 tons;

• Length: 58.00 m;

• Width: 10.30 m;

• Forecastle deck height: 6.40 m;

• The height of the main deck, aft: 4.70 m;

• Draft: 2.60 m;

• Speed: not less than 15 kt;

• Range: not less than 2,500 nautical miles.

Due to the destination of the ship should be examined cases in which the ship may be exposed to the sea mine explosion. It is therefore necessary to distinguish underwater explosion (UNDEX) derived from the enemy charge and the explosion on board in case of a failure while setting up mine warfare. The effects of explosions are difficult to predict. In case of contact explosions which are explosions on board of the ship, it is possible to determine the pressure distribution derived from such an explosion by using the CONWEP procedure which is not described further.

Fig. 1The Kormoran II minehunter designed for the minefield setting

Nowaday, many CAE programs uses CONWEP procedure to simple calculations for free explosions in the air. The value of pressure is determined by the program based on the weight of the load and then using TNT equivalents calculated for the respective loads. A similar situation can be applied to underwater explosions. In order to properly calculation of pressure wave incident on the structure shall be familiar with the phenomenon, analytical description of the pressure wave used by the program and simplifying assumptions.

Underwater explosion is a rapid equilibrium violation caused by detonation of explosive in water. During explosion the charge is converted from solid in to has product of volume equal to volume of explosive material, temperature range T≈ 3300 K and pressure p≈14000 MPa. With so much pressure gas bubble is formed which causes a spherical shock wave that propagates in the medium [3, 4]. The first of the researchers of this phenomenon was the R. H. Cole whose works are the basis for further research. The parameters of the shock wave and the gas bubble parameters are being investigated further. Currently the main of researchers are T. L. Geers, K. S. Hunter and C. K. Park. Most previous publications describes the behavior of gas bubbles, often omitting same shockwave. In this work, omit the description of the bubble in order to focus on the pressure wave from the detonation. This wave in the initial stage moves at a speed of the order of v≈ 5000-8000 m/s. Then, the water molecules act on the water layer adjacent losing speed and travel further the speed of sound in water, which is approximately c0≈ 1500 m/s. Pressure impulse was described by researchers in various ways. It is hard find a universal solution, however, due to the development of measuring instruments in this work decided that for the calculations related to underwater explosion the latest models, ie. model developed by Geers’a and Hunter will be best to use. The following figure shows the profiles described by researchers over the years [5].

Fig. 2The pressure wave profiles described by researchers over the years [5]

Cole formulas shown below. They were determined by the measurement of detonation of 70-136 kg of TNT. Profile of the pressure wave is shown by the following formula:

where: m – charge mass, kg; r – distance of epicenter, m; pmax – peak pressure, MPa; θ – time constant, ms; t – time, ms.

Cole models were presented in 1948. Since that time was corrected by J. S Navagin, W. Stiepanow, and T. L. Geers and K. S. Hunter, the latter formulas are shown below:

where: ac – charge radius, m; Pc= 1.42; vc= 0.992; A= 0.13; B= 0.18.

In the case of explosions of any material other than TNT it is possible to use a TNT equivalent. The equivalent is very difficult to determine because of the multitude of comparative parameters. In practice the heat of explosion, detonation speed and pressure coming from the explosion at a test distance are compared. As the table shows the equivalents determined on the basis of various parameters differ significantly. A good solution is to use average values.

The task of explosions, should be regarded as an appropriate task described by the equation of motion. Recognition of this phenomenon using the finite element method requires the use of an appropriate equation in the form:

where: K – stiffness matrix; M – matrix of inertia (density matrix); C=αM+βK – damping matrix, where α and β are constant coefficients; U,U˙,U¨ – displacement, velocity and acceleration vector; F – load vector, ε˙ – strain rate vector; t – time; R – distance from epicenter; mTNT – TNT charge mass.

In the task load vector F have been set as a time-varying pressure wave in accordance with Eqs. (1), (2) or (3), (4). Other expressions of Eq. (5) in this work are omitted, because they describe the structure reaction while this work refers only to load.

Polish Navy is currently using 4 types of sea mines. Among other, uses OS mine which contains 110 kg TNT (Source: Laboratory of Underwater Weapons Polish Naval Academy). The paper examined a case in which such a charge detonates around the ship of Kormoran II project. Using CAD programs reflected outer shell of the ship, and then implements it into CAE program in which the relevant calculations were carried out.

Charge was placed at a distance of 9 m below the keel of the ship. According to Eq. (3) the pressure on the front of the shock wave is 25.67 MPa. The pressure calculated by the program was 29 MPa. The difference results from the strong relationship between the size and shape of mesh elements. In the paper [5] the effect of mesh size and shape on the propagation of the shock wave were analyzed. It was found that the 4 – nodal tetragonal of size of 0.5 m elements are best to use.

Fig. 3Pressure distribution of non-contact TNT explosion under the keel of the ship

a)

b)

Fig. 4Pressure distribution of non-contact TNT explosion under wet surface of the ship

The pressure distribution of non-contact mine explosion in the acoustic medium in different time points were shown. The pressure wave reaches the ship’s hull and it is reflected resulting in a negative pressure. The decrease of pressure can lead to cavitation under the keel of the ship. In this case, if the ship is treated as a beam supported continuously, due to cavitation a discontinuity of the support will appear and the ship would be load with the forces of gravity. In this work, this phenomenon is omitted.

Study describes the phenomenon of underwater explosion. The model used in the study is the latest model of the shock wave in this type calculations. Problems that results from the application of the acoustic medium and its response for the given load. Also highlights the FEM equation that describes the dynamics task. The complexity of this task depends on the author and expected results. In case of underwater explosions weight loss, elements self – contact or hardening of the material due to the strain rate can be included. The influence of various parameters on the results will be examined in further publications.