4.1. Bending Mode
As can be seen from the pictures in the validation, when the SHG suffers from sagging damage subjected to the near-field explosion bubble, the two ends move like rigid bodies within a small elastic-plastic limit, and the overall bend is mainly composed of large plastic deformation near the midship, as shown in
Figure 6. Moreover, the crease pattern at the bend position has certain regularity. In the process of sagging damage of the SHG, the upper part of the side plates becomes unstable under compression and gradually expands to form several plastic hinge lines. The deformation range of the lower part of the side plates is generally narrower than that of the upper part, specifically shown as W-shaped creases. The deformation of the bottom plate is relatively simple. The plastic hinge lines extend transverse and straight through the whole width. These plastic hinge lines form the overall bend in the SHG. After the plastic hinge lines are formed, the rotation around the yield area becomes easier, and the overall sagging amplitude is often further increased when subjected to bubble collapse force. For the W-shaped creases’ local deformation mode on the side plates in this paper, an SHG sagging bend was composed of ten plastic hinge lines, eight on the side plates and two on the bottom plate, as shown in
Figure 7.
As the core deformation damage zone, it is important to predict the length of the potential bend zone when facing a near-field underwater explosion load. Therefore, the length of the bend zone, which is not affected by the transverse bulkhead, is analyzed below. Each part of the hull plate rotates around the plastic hinge lines during the bending process. That is, the bending moment plays a major role. In contrast, the energy dissipated by the expansion and contraction of the hull plates is small; therefore, it is ignored.
is used to indicate the length of the bend zone on the undeformed SHG, as shown in
Figure 7.
represents the average bend moment of the SHG during the bending process. The angle that, at one end, rotates around the middle cross-section is recorded as
. It is assumed that the length of each plastic hinge line on the side plates
l has a function related to the length of bend zone
and the height of the side plate
H; then, the length of different plastic hinge lines can be denoted as
. The rotation angle of both sides of each plastic hinge line is positively related to
, assumed to be
. The length of the plastic hinge lines on the bottom plate is the same as the ship width
W, and the rotation angle equals
. As the work of external bend moment equals the sum of dissipated energy of the side and bottom plates’ plastic hinge lines, we have:
where
is the ultimate bend moment per unit length. According to the principle that the
value should make the value of
the minimum, let
; therefore, we have:
According to the geometric relations,
is a constant for the determined plastic hinge lines mode. Therefore,
has a function that only depends on
H. To find the relation between
and
H, keeping the width of SHG unchanged, the sagging damage responses with height–width ratios of 0.4, 0.5, 0.6 and 0.7, under the same load conditions as
Section 3.2, are calculated. The creases are all in a similar W-shaped pattern. The lengths of the bend zones are listed in
Table 3 and fitted in
Figure 8. The three functions in
Figure 8 can fit the relationship between the bend zone’s length and the side plates’ height. In sum, the correlation coefficient of the formula
is closest to 1, and the form is convenient for application; therefore, this linear relation is adopted.
The length data of the bend zone can be obtained from the calculation results of the finite element model using the following method. The top edge rotation angle was plot around the z-axis (UR3), the bottom edge rotation angle was plot around the y-axis (UR2), and their differential curves are shown in
Figure 9. Points 1-5 are the endpoints of plastic hinge lines. The angle at the crease endpoint changes significantly, that is, the peak points of the differential curve.
4.2. Effect of Transverse Bulkheads on Bend Damage
When the transverse bulkheads are located in the potential bend zone, this will obviously have a greater influence on the bend and creases. When the transverse bulkheads are located outside the potential bend zone, this may also have some influence. Therefore, to investigate the transverse bulkhead’s effect on the bend damage characteristics of SHGs, the internal structure was redesigned based on the SHG used in the validation.
Figure 10 presents the SHG configuration diagram. The origin of the coordinate is located in the center of the bottom plate, and the arrangement of the transverse bulkheads is symmetrical about the yoz coordinate plane. By changing the position, eight kinds of SHGs are designed, as shown in
Table 4. The numbering rules are as follows: SHG is the abbreviation of simplified hull girder, and the number in brackets is the dimensionless position of transverse bulkheads
. The calculation formula of
is:
For example, SHG(0) indicates that the transverse bulkhead is located at the mid-section; SHG(1) indicates that the transverse bulkhead coincides with the end plate, i.e., there is no transverse bulkhead. Through the structural response results of SHGs with the staggered and orderly distribution of transverse bulkheads, the bearing effect of transverse bulkheads can be analyzed to some extent.
To study the bend damage patterns, all SHGs were subjected to the same load conditions as the verification experiment. After checking the results, the unbent SHGs were allowed to bend to a certain extent by keeping the charge unchanged and gradually reducing the detonation distance. The final calculation cases are shown in
Table 5.
There are several obvious different bend characteristics in the calculation results. These in clude whether the position is central or off-center, whether there are one or two bends, and whether the size is affected by the transverse bulkheads.
Figure 11 shows the crease evolution process of the bend position with the following characteristics: (a) size is not affected by the transverse bulkheads and there is a single center bend, (b) size is affected by the transverse bulkheads and there is a single center bend, (c) size is affected by the transverse bulkheads and there is a single off-center bend, and (d) size is affected by the transverse bulkheads and there are double off-center bends.
When the transverse bulkheads are located outside the bend zone, it is difficult to effectively limit the crease formation process. The crease expansion is only affected by the hull plates, as shown in
Figure 11a. When the transverse bulkheads are located in the bend zone, the creases are limited between the two transverse bulkheads, the space available for side plate deformation becomes small, and the overall bend resistance is effectively improved, as shown in
Figure 11b. When the transverse bulkhead is located directly above the explosive, the side plates at the corresponding position are firmly fixed, and it is difficult to produce a transverse deformation. Therefore, it is difficult for the two ends of the SHG to rotate with the transverse bulkhead as the center; therefore, the bend position is offset to both sides of the transverse bulkhead. Depending on the load conditions, there can be one or two bends, as shown in
Figure 11c,d. As for
Figure 11d, because the explosive is too close to the water surface, the air is sucked in during the bubble contraction, which destroys the internal pressure state, leading to an early recovery time for sagging amplitude. In general, the relative longitudinal position of the explosive and the transverse bulkhead will affect the bend position and size, and the load strength will affect the bend quantity. In the results of this paper, all the side plate creases are W-shaped. However, more crease patterns may appear under different conditions, which requires further study.
Table 6 lists the damage results for SHGs in different cases. The ratio of transverse bulkhead position to half of the bend zone length was used to characterize the parameter of the relative position of the transverse bulkhead and bend zone, i.e.,
. The calculation formula is
In cases 1 and 2, the SHGs move in whipping mode, so there is no bend zone. The large error in case 3 is due to the influence of the transverse bulkhead, which increases the length of the bend zone by about 12%, indicating that a transverse bulkhead will slightly increase the length of the bend zone. The error in case 9 is large because the transverse bulkheads are located within the predicted bend zone, which restricts the outward development of plastic hinge lines. The same SHG was used in cases 10 and 11. However, the error in case 10 is larger because the load input energy is smaller than in case 11. Part of the load energy is absorbed by the warping deformation of the left hull plate, resulting in a smaller bend zone length. It can be seen from case 11 that when the load is strong enough, the length of the off-center bend zone with only one constrained side can also be well predicted. It can be seen from cases 4-8 that the transverse bulkheads with have little influence on the length of the bend zone. Thus, when the bend zone is less affected by the transverse bulkheads, the prediction error of the formula is less than 5%, which can reflect the potential bend zone length and provide a reference for the anti-impact work of ships.
4.3. Effect of Transverse Bulkheads on Overall Bend Strength
To analyze the effect of the transverse bulkhead position on the overall bend strength of SHGs, the overall response results of eight different SHGs under the same load in cases 1–8 are analyzed below.
Figure 12 shows the contours of different SHGs at various times. The arrows indicate the direction of motion, and the symbols on the curves indicate the position of transverse bulkheads. At an early stage of the explosion, the middle of the SHGs rapidly displaces under the combined action of a spherical shock wave and bubble expansion boundary. Moreover, the displacement of both ends is relatively slow due to the inertial effect, so the overall structure presents hogging bend deformation, as shown in
Figure 12a. In
Figure 12b, the main action time of the shock wave has passed. Under the influence of material resilience, the overall deformation gradually enters the unloading stage, and the bend amplitudes decrease to varying degrees. At this moment, the deflection amplitude of all structures is within 0.0025 m, and the ratio to the total length is 1.67‰, so the SHGs can be consideed close to the horizontal state. In
Figure 12c, all the structures show a sagging bend, and the deflection amplitude significantly varies among the different structures. After the bubble reaches its maximum volume, it starts to rapidly shrink the boundary. All the structures are driven downward by this, which also causes the sagging bend amplitude of SHG(0.2), SHG(0.33), SHG(0.4), SHG(0.6), SHG(0.8), and SHG(1) to expand further and form sagging damage. However, the deflection amplitude of SHG(0) and SHG(0.1) is not obviously affected, as shown in
Figure 12d,e. After the bubble shrinks to the minimum volume, it expands again. The SHGs begin to move upward under the comprehensive action of the expansion bubble boundary, pulsating pressure, jet flow, and restoring force. The sagging amplitudes are restored to different degrees. The deformation of SHG(0) and SHG(0.1) turns into an upward bend, as shown in
Figure 12f. The general response mode of SHG(0) and SHG(0.1) is a whipping motion, and the rest of the SHGs appear to undergo different degrees of sagging damage.
In summary, the response behavior and damage degree of SHGs with different transverse bulkhead configurations are significantly different under the same load. This shows that transverse bulkheads can contribute important longitudinal strength when facing a near-field underwater explosion load, and the effect of transverse bulkheads differs at different positions. When predicting the overall damage behavior of the structure, attention should be paid to the longitudinal position parameters of different transverse bulkheads.
The dimensionless deflection
is introduced as a parameter to characterize the relative bend degree of SHGs.
Figure 13 shows the dimensionless deflection time–history curves of each structure. For these SHGs, if the transverse bulkhead is at the boundary of the predicted bend zone, the corresponding structure number is SHG(0.162). As seen from
Figure 13, there are great differences between the deflection of SHGs with transverse bulkhead position parameters greater than 0.162, indicating that the transverse bulkhead outside the bend zone greatly influences the overall damage. In addition, SHGs with different transverse bulkhead configurations show little difference in hogging behavior but a great difference in sagging behavior. This is because, when bent downward, the upper part of the side plate is compressed, and the lower part is stretched. The bend requires yielding and folding in the upper part to make space. Similarly, in the hogging bend, the lower part of the side plates will appear more yielding and folding. For the SHG structure without the deck used in this paper, the lower part of the side plates is orthogonally connected with the bottom plate to form a strong constraint, which is more difficult to yield than the upper part.
Further, the maximum sagging amplitudes of each structure are given as statistics in
Figure 14. For the SHGs with transverse bulkheads, the response changes from whipping motion to sagging damage as the
increases. With the change in position, the restraining ability of the transverse bulkhead to sagging amplitude is shown to have a three-stage pattern, from strong to weak. The transverse bulkhead in the bend zone can significantly reduce the overall bend amplitude. The transverse bulkhead located outside the bend zone also has a non-negligible effect. For example, SHG(0.8) has a 94.5% change in deflection relative to SHG(1) without transverse bulkheads. The sagging deflection of SHG(1) without a transverse bulkhead is between SHG(0.2) and SHG(0.33) and is not the maximum. This indicates that the transverse bulkhead does not always play a positive role in the bearing near-field explosion bubble load. This is because the transverse bulkhead at a far position will reduce the load-sharing capacity of the structure outside the bend zone, causing the load to become more concentrated in the bend zone. The deflection of SHG(0) is slightly greater than that of SHG(0.1). The reason for this is that the bend position of the transverse bulkhead at the mid-section will be offset, and these two have different deformation modes.
The above results show that the distant transverse bulkhead has a significant influence on the deflection, and it is difficult to properly characterize the overall bend resistance with the parameters of a single section. Therefore, based on the moment of inertia, the equivalent moment of inertia considering the influence of transverse bulkheads at different positions was constructed for evaluation:
where
is the equivalent moment of inertia of SHGs,
is the moment of inertia of the SHGs without the transverse bulkhead section,
is the moment of inertia of the transverse bulkhead section, and
is the weight coefficient of the transverse bulkhead at different positions. For this paper,
= 80.88 cm
,
= 1666.66 cm
. Assuming that the maximum sagging deflection of SHGs is inversely proportional to the equivalent moment of inertia when the load and ship shape is constant, the weight coefficient function of transverse bulkheads at different positions can be calculated according to Formula (
11) and
Figure 14, as shown in
Figure 15, and the expression is:
where
is the boundary position of the bend zone. It can be seen that, when facing a near-field underwater explosion load, the transverse bulkhead in the bend zone makes a significant positive contribution to the overall bend strength. Therefore, it can be considered to properly strengthen the transverse bulkhead in the potential bend zone in the actual ship, such as using high-strength steel. The transverse bulkhead outside the bend zone may play a positive or negative role. Although the contribution coefficient is relatively small, the influence on the equivalent moment of inertia cannot be ignored due to the large difference between
and
. It should be noted that the transverse bulkhead weight coefficient function of Equation (
12) is mainly for the SHG configuration used in this paper. Different SHGs with more transverse bulkheads will be more complicated and require further study.