1. Introduction
The ultimate load capacity of a ship is a crucial parameter for assessing its safety and reliability and has been extensively investigated by scholars. To date, numerous well-established methodologies have been developed to compute the static ultimate load capacity of ship structures. Smith [
1] proposed an incremental-iteration method for calculating the ultimate load capacity of a hull girder. In this method, the cross-section is divided into different types of calculation units, each with its own specific stress–strain curve. By gradually increasing the curvature of the cross-section, calculating the stress of each unit, and integrating the stress of all units on the cross-section, one can obtain the bending moment-curvature curve where the peak value of the curve is the static ultimate load capacity of the hull girder. The curvature of the cross section refers to the angle at which the cross section of the hull girder rotates around the neutral axis. Because of its simplicity and high calculation accuracy, the incremental-iteration method has become a standard method for calculating the static ultimate load capacity of hull girders in classification societies. Tanaka et al. [
2] extended the incremental-iteration method and established a method for calculating the ultimate load capacity of ship under combined bending and torsion, considering the influence of shear stress. Ueda and Sherif [
3] proposed a numerical method called the Idealized Structural Unit Method (ISUM) for calculating the ultimate load capacity of ship structures. The structural units in the ISUM are larger than the elements in the finite element method (FEM), which greatly reduces the computation time [
4]. Lindemann and Kaeding [
5] used the ISUM to study the ultimate load capacity of the stiffened plates under lateral pressure and in-plane compression. Underwood et al. [
6] proposed a new ISUM that can assess the static ultimate load capacity of damaged structures through collapse analysis of the stiffened plates.
In recent years, with the rapid improvement in computer performance, the nonlinear FEM has been widely used to compute the ultimate load capacity of ship structures. The nonlinear FEM is characterized by changes in structural stiffness in the FE analysis. There are three kinds of nonlinearity: geometric nonlinearity, material nonlinearity and boundary nonlinearity. The term “geometric nonlinearity” mainly refers to large deformations. The term “material nonlinearity” mainly refers to the consideration of the yield strength of the material. The term “boundary nonlinearity” mainly refers to changes in boundary conditions encountered during the analysis. It is common to encounter both geometric nonlinearity and material nonlinearity when analyzing ultimate load capacity. Thus, FE analysis is nonlinear. The ultimate load capacity of ship structures can be predicted by empirical formulas based on the numerous results obtained by the FEM. Zhang and Khan [
7] proposed a semianalytical empirical formula and then performed numerical simulation of stiffened plates using FEM to determine the coefficients of the empirical formula. Kim et al. [
8] proposed an empirical formula to predict the ultimate load capacity of the stiffened plates subjected to longitudinal compression through numerical calculations of the stiffened plates with T-bar and flat-bar stiffeners. Xu et al. [
9] studied the effects of lateral pressure and stiffener types on the collapse behavior of the stiffened plates through nonlinear FEM and derived an empirical formula for predicting the ultimate load capacity of the stiffened plates under the combined actions of axial compression and different levels of lateral pressure based on the results of the numerical calculation.
In all of the above methods, the ultimate load capacity of ship structures is evaluated under static or quasistatic conditions. However, when a ship encounters extreme conditions such as large freak waves or underwater bubbles, the amplitude of the load on the ship is very large and the load duration can be reduced to the order of milliseconds [
10,
11]. Obviously, the load on ship structures cannot be considered as a static or quasistatic load in such situations. In particular, the incident report from
MOL Comfort indicates that dynamic loads may be one of the causes of ship structure collapse [
12]. Thus, it is meaningful to study the dynamic ultimate load capacity of ship structures. Yamada [
13] studied the effect of strain rate on the ultimate load capacity of container ships by using a half-sinusoidal dynamic load. The numerical calculation revealed that the strain rate has a significant effect on the ultimate load capacity and that the ultimate load capacity could be increased by 10–20% when the strain rate was considered. Jagite et al. [
14] found that in calculating the dynamic ultimate load capacity of ship structures subjected to wave loads, the effect of strain rate can be ignored for container ships. Yang et al. [
15] studied the dynamic ultimate load capacity of rectangular plates under axial compression. These authors used the FE results for numerous ship plates to derive an empirical formula for predicting the dynamic ultimate compressive capacity of ship plates. Yang and Wang [
16] studied the dynamic buckling of the stiffened plates under in-plane impact loads by theoretical derivation and considered the influence of rotational constraint stiffness on the dynamic response of the stiffened plates. Paik [
17] studied the dynamic ultimate load capacity of plates under axial compression loads with different loading speeds. The experimental results show that as loading speed increases, the ultimate compressive capacity of the plate also increases gradually. Liu et al. [
18] studied the dynamic failure of rectangular plates under lateral impact via impact experiments and FE analysis. Jagite et al. [
19] studied the dynamic ultimate load capacity of the stiffened plates under axial compression and lateral load in a parameterized way. Load values obtained by parameterization are more consistent with those measured under real-world conditions. The limitations of the existing strain rate model are also discussed.
The current method for determining the dynamic ultimate load capacity of ship structures relies primarily on the concept of dynamic buckling. Ari-Gur and Simonetta [
20] defined the dynamic failure criterion as the case in which a small increase in the dynamic load amplitude causes a sudden increase in the dynamic response of the structure. Xiong et al. [
21] studied the dynamic ultimate compressive load capacity of the stiffened plates based on this criterion and derived a reasonably accurate empirical formula. Yang et al. [
22] proposed a one-time thickness-deformation method for determining the dynamic ultimate load capacity by studying a container ship’s bow. Their formula defines the dynamic ultimate load capacity as the amplitude of the dynamic load at which the maximum deflection of the structure is equal to one times plate thickness. Yang and Wang [
16] proposed a new formula for determining the dynamic critical buckling load. According to this formula, when the slope of the load-displacement curve at any given point reaches 35, the corresponding load at that point is the dynamic critical buckling load. There are other methods for determining the dynamic ultimate load capacity, but they are similar to those already described and will not be discussed further here. However, these methods all essentially use the deformation of the structure to determine the dynamic ultimate load capacity. There are drawbacks to applying these methods in some cases, as will be described below. Therefore, we hope to develop a novel method for determining the dynamic ultimate load capacity of ship structures.
Considering the context described above, this study uses a nonlinear FEM to investigate the dynamic ultimate load capacity of the stiffened plates as one of the primary components of ship structures. This work proposes a novel “two-step” approach to determine the dynamic load capacity of the stiffened plates based on the static ultimate load capacity after impact. This method is theoretically applicable to a variety of ship structures. Subsequently, the failure mode of stiffened plates under impact load is discussed. Considering that real ship structures may experience continuous impacts with different load amplitudes, this work includes an analysis of the dynamic ultimate load capacity of the stiffened plates subjected to continuous impact load, as well as of the influence of the impact load cycles and impact load sequence on the dynamic ultimate load capacity of the stiffened plate. Finally, the applicability of the two-step approach to hull girders is studied.
3. Numerical Results and Discussion
The purpose of this study is to propose a new method for determining the dynamic ultimate load capacity of ship structures. In order to make this method more universal, five different sizes of the stiffened panels were designed according to
Table 1, as shown in
Table 3. Zhang and Khan [
34] found that the plate slenderness ratio
ranges from 1.0 to 4.5 and that the stiffener slenderness ratio
ranges from 0.15 to 0.95 on ship structures. As shown in
Table 3, the range of
and
in the designed model are both within the range reported by Zhang and Khan.
3.1. Static Ultimate Load Capacity Results
Before calculating the dynamic ultimate load capacity, the static ultimate load capacity and the first-order vibration period of each model were calculated first. The impact load is determined by the impact duration and the impact load amplitude. Therefore, the static ultimate load capacity and the first-order vibration period can be used as the reference values for the amplitude of the impact load and the impact duration, respectively. Then, the nondimensional impact load amplitude and the nondimensional impact duration can be defined as and , respectively.
The accuracy of FEM can be verified by comparing the FE results with those of empirical formula, as discussed by Xiong et al. [
21]. The results for static ultimate load capacity were verified by comparison with the empirical formula derived by Xu et al. [
9]. This formula is relatively accurate for predicting the static ultimate load capacity of T-bar stiffened plates. The empirical formula is as follows:
The results for static ultimate load capacity of FEM, along with the empirical formula and errors, are shown in
Table 4. All error values are within 10%. The maximum error is 6.41%. Xiong et al. [
21] compared the static ultimate load capacity obtained by FEM and the empirical formula. In that study, the maximum error was 7.14%, and all errors were within 10%. Thus, the results of the FEM calculation in this study are sufficiently precise.
3.2. The Determination Criterion of the Dynamic Ultimate Load Capacity
The dynamic ultimate load capacity of each stiffened plate model was computed based on the first vibration period, considering 5 different impact durations: 0.2T0, 0.5T0, T0, 2T0 and 5T0. According to the calculation procedure of the two-step approach, the maximum value of the dynamic response of the stiffened plate was recorded in step 1, and then the static ultimate load capacity after impact was recorded in step 2. The peak value of the dynamic response is the maximum axial displacement of the loaded edge, i.e., . The maximum axial displacement of the stiffened plate under static ultimate load capacity is . Thus, the nondimensional maximum axial displacement is defined as . The static ultimate load capacity of the stiffened plate after impact is denoted as . The nondimensional static ultimate load capacity after impact is defined as .
Model 3 was taken as an example to illustrate the calculation process in the two-step approach. The process of applying the two-step approach is similar to that used for other models. When the impact duration
, the results obtained by using the two-step approach are shown in
Figure 10. The axial displacement of the stiffened plate over time subjected to impact loads with different amplitudes in step 1 is shown in
Figure 10a. The axial displacement of the loaded edge can also be called “end shortening”. As shown in the figure, the stiffened plate vibrates after impact and then tends to stabilize. As load amplitude increases, the maximum axial displacement also increases gradually. When
, the stiffened plate undergoes almost no plastic deformation after impact. However, when
, the stiffened plate undergoes plastic deformation, and as the impact load amplitude increases, the plastic deformation of the stiffened plate increases gradually.
When the vibration of the stiffened plate becomes stable, the calculation in step 1 ends. The calculation in step 2 then begins, and forced displacement is applied at both ends of the stiffened plate to calculate the static ultimate load capacity.
Figure 10b shows the compression load of the stiffened plate over time after impacts of different amplitudes. Note that step 1 was carried out in Dynamic-Implicit, with a total calculation time of 0.3 s, which is the real time. By contrast, step 2 was carried out in Static-General; in this step, the time is not real and defaults to 1 s. The maximum compressive load is the static ultimate load capacity after impact. It shown in
Figure 10b that when
, the static ultimate load capacity after impact
does not decrease and is equal to the static ultimate load capacity
. However, when
, the structure undergoes plastic deformation, which indicates that low levels of plastic deformation have no influence on the static ultimate load capacity of the stiffened plates. As the impact load amplitude increases, the static ultimate load capacity after impact decreases obviously at
.
The curve of the nondimensional maximum axial displacement versus nondimensional load amplitude and the curve of nondimensional static ultimate load capacity after impact versus nondimensional load amplitude from Model 3 can be obtained by gradually, iteratively increasing the amplitude of the impact load for each impact duration. The results are shown in
Figure 11.
By observing the trend in the maximum axial displacement versus the impact load amplitude in
Figure 11a, we find that the maximum axial displacement gradually increases as the impact load amplitude increases and that the slope of the curve increases continuously. However, it shown in
Figure 11b that the static ultimate load capacity after impact eventually stabilizes as the impact load amplitude increases, i.e., the static ultimate load capacity will not decrease after impact. Then, as the amplitude of the impact load continues to increase, the static ultimate load capacity after impact begins to decrease. The slope of the curve in
Figure 11b is close to zero at the beginning, then begins to decrease; eventually, the slope becomes almost constant.
Therefore, we can find a point on the curve after which the slope of the curve remains constant, i.e., where the second derivative is 0, and define of the point as the dynamic ultimate load capacity. The dynamic ultimate load capacity is denoted as .
According to the above-described method of determining dynamic ultimate load capacity, the dynamic ultimate load capacity of Model 3 under different impact durations can be obtained, as shown in
Table 5. The percentage reduction of static ultimate load capacity after impact is given in the table. After the stiffened plate is subjected to the impact load whose amplitude is the dynamic ultimate load capacity, its static ultimate load capacity is reduced by 1.4–2.9%.
3.3. Dynamic Ultimate Load Capacity of Stiffened Plates
According to the procedure of the two-step approach, the calculation results from the other models are shown in
Figure 12. The dynamic ultimate load capacity of each model under different impact durations is shown in
Figure 13. As the impact duration increases, the dynamic ultimate load capacity of the stiffened plate gradually decreases. By observing the calculation results of Model 1 in
Figure 13, we find that the dynamic ultimate load capacity is obviously smaller when
. Moreover, we observe the curve of the nondimensional maximum axial displacement versus nondimensional load amplitude of the stiffened plate in
Figure 12a, and we find that the nondimensional maximum axial displacement corresponding to
is the largest when
. The cause of this phenomenon may be that higher-order deformation modes are activated under short impact duration.
Figure 14 shows the displacement distribution at dynamic ultimate state of Model 1 with
and the vibration mode with 103.36 Hz. As shown in
Figure 14, the displacement distribution and vibration mode of the stiffened plate are similar. And the vibration period is about
when frequency is equal to 103.36 Hz, which is close to the impact duration
. Therefore, when
, the impact load activates the high-order vibration mode of the stiffened plate, causing the resonance of the stiffened plate, thus the axial displacement of the stiffened plate increases significantly.
Table 6 shows the percentage reduction of static ultimate load capacity after impacts of different duration. The static ultimate load capacity of after impact all models is reduced by no more than 3%. It indicates that the static ultimate load capacity of the stiffened plate will decrease slightly when the stiffened plate is subjected to the impact load whose amplitude is the dynamic ultimate load capacity. We propose another method to determine the dynamic ultimate load capacity, that is, the static ultimate load capacity is reduced by 3% after the structure is subjected to the impact load, then the amplitude of this impact load is the dynamic ultimate load capacity. This method is suitable for the case where ship structures are subjected to impact loads of complex shapes.
3.4. Dynamic Failure Mode of Stiffened Plates
The failure mode of the stiffened plate can be divided into 6 types [
34], as follows: mode I, overall failure mode; mode II, failure of the local plate between stiffeners; mode III, beam-column-type failure of the stiffeners with the attached plate; mode IV, local buckling of the stiffener web; mode V, tripping of the stiffener; mode VI, overall yielding of the stiffened plate. When stiffened plates are subjected to dynamic impact loads, the failure mode will change.
Figure 15 shows the stress diagram of Model 1 in the ultimate state. As shown in
Figure 15a–e, the failure mode of the stiffened plate at the dynamic ultimate state changes as the impact duration increases. The failure mode of the stiffened plate is mainly failure of the local plate (mode II) when
. The deformation of the stiffener is relatively small. As the impact duration increases, both the tripping of the stiffener and the failure mode of the local plate arise when
, i.e., the failure modes are mainly mode II and mode V. The beam-column-type failure of the stiffened plate (mode III) is the main failure mode when
. In addition, the failure mode at this time is the same as the static failure mode of the stiffened plate in
Figure 15f.
Figure 16 shows the stress diagram of Model 3 in the ultimate state. It can also be seen from
Figure 16a–e that as the dynamic impact duration increases, the failure mode of the stiffened plate changes. The failure mode of the stiffened plate is mainly failure of local plate (mode II) when
. The deformation of the stiffener is very small. when
, the stiffener deformation increases, several half waves appear and the main failure modes are mode II and mode V. When
, the tripping of the stiffener (mode V) is the main failure mode. The failure mode in
is the same as the static failure mode of the stiffened plate in
Figure 16f.
In summary, the impact duration can affect the dynamic failure mode of the stiffened plate. The dynamic failure mode of the stiffened plate for is the same as the static failure mode.
3.5. Effect of the Number of Impact Load Cycles on Dynamic Ultimate Load Capacity
Figure 17 shows the schematic diagram of the impact load cycles: one cycle, two cycles and three cycles. Using the two-step approach, the dynamic ultimate load capacity of Model 2 under different numbers of impact load cycles is calculated when
.
Figure 18 shows the calculation results obtained by the two-step approach. As shown in
Figure 18a, the greater the number of impact cycles, the larger the maximum axial displacement of the stiffened plate under an impact load of the same amplitude. This effect occurs because the greater the number of impacts, the larger the load impulse applied to the stiffened plate and the larger the maximum axial displacement. Based on the results shown in
Figure 18b, the dynamic ultimate load capacity under one impact, two impacts and three impacts can be determined as
,
and
, respectively. Therefore, the dynamic ultimate load capacity decreases as the number of impact cycles increases.
3.6. Effect of the Impact Load Sequence on Dynamic Ultimate Load Capacity
Real ship structures may be subjected to continuous impact loads with different amplitudes. Therefore, it is meaningful to study the dynamic ultimate load capacity of ship structures under continuous impact. We assume that the stiffened plate is impacted three times continuously, each time with different amplitudes. As shown in
Figure 19, there are a total of six different sequences of the impact load s. The duration of each impact is denoted as
. The three amplitudes of the impact load are
,
and
, respectively. For ease of recording, the sequence of the impact load s is written as 1, 2, 3. For example, the sequence of the impact load amplitudes in
Figure 19a is
,
,
, which can be denoted as “A-123”.
With Model 2 as the case study and static ultimate load capacity as the reference value, the dynamic ultimate load capacity is calculated by using the two-step approach when
. The calculation results are shown in
Figure 20. As higher-order modes may be activated, as observed in the course of continuous impacts, the curve of the nondimensional static ultimate load capacity after impact versus nondimensional impact load impact becomes non-smooth. The dynamic ultimate load capacity is derived from the percentage reduction in the static ultimate load capacity after the impact. According to the above discussion, the static ultimate load capacity after impact cannot be reduced by more than 3%, so we have chosen 3% reduction the criterion for determining the dynamic load capacity. Thus, results for different impact load sequences can be obtained, as shown in
Table 7. The 3% reduction is based on the nondimensional results from the five models of stiffened plates under five different impact times. Therefore, it is reasonable to take a 3% reduction as the criterion for determining the dynamic ultimate load capacity. When the last impact load is the one with the highest amplitude, i.e., when the load sequence is “A-123”, the dynamic ultimate load capacity of the stiffened plate is the smallest. Because the stiffened plate has been subjected to two impact loads with different amplitudes and then subjected to an impact load with a higher amplitude, the deformation and plastic region of the stiffened plate will be significantly larger than in other cases, reducing the dynamic ultimate load capacity. For comparison, when the load sequence is “A-321”, the highest-amplitude load occurs first. After that impact, the smaller loads are not able to expand the plastic region further. Additionally, plastic regions have already been generated when the smaller loads occur, resulting in smaller deformations and a greater ultimate load capacity.
3.7. Application of the Two-Step Approach to a Hull Girder
In order to illustrate the applicability of the two-step approach to other ship structures, the two-step approach was used to calculate the dynamic ultimate load capacity of a hull girder.
Figure 21 shows the schematic diagram of a single-frame model of a container-ship hull girder. The two reference points (RP-1 and RP-2) and all nodes at both ends of the model are coupled by coupling constraints. Both the boundary conditions and the impact bending moment are applied to the reference points. One end of the model is fixed, and the other end is subjected to impact bending moments. The shape of the impact bending moment is also half-sinusoidal, as shown in
Figure 8. According to the natural vibration period of the single-frame model,
is selected as the calculation case and the dynamic ultimate load capacity of the hull girder in the sagging state is calculated.
The procedure for calculating the dynamic ultimate load capacity of the hull girder by the two-step approach is similar to that used for the stiffened plates. First, the static ultimate load capacity and the static ultimate rotation angle of the hull girder are calculated. Thus, the nondimensional maximum dynamic response of the hull girder (the rotation angle of the cross section of the hull girder during impact.) and nondimensional amplitude of the impact bending moment can be defined as and , respectively. The static ultimate load capacity of the hull girder after impact is denoted as .
Figure 22 shows the results of calculating the dynamic ultimate load capacity of the hull girder by the two-step approach. From the trend shown by the curve in
Figure 22b, it can be determined that the dynamic ultimate load capacity when
is
and that the static ultimate load capacity after the impact is reduced by 5%. In practice, ship designers can determine the percentage reduction in static ultimate load capacity after impact (e.g., 8% or 10%) according to their experience and needs. Therefore, the two-step approach can also be used to determine the dynamic ultimate load capacity of the hull girder.
4. Conclusions
This study proposes a novel “two-step” approach for determining the dynamic ultimate load capacity. Taking a stiffened plate as a case study, the nonlinear FEM-based two-step approach is discussed in detail. Moreover, the applicability of the two-step approach to a hull girder is verified. The main conclusions of this study are summarized as follows:
(1) A criterion for the evaluation of the dynamic ultimate load capacity of ship structures is proposed in this paper based on the curve of nondimensional static ultimate load capacity after impact versus nondimensional impact load amplitude. These data can yield a reasonable determination of dynamic ultimate load capacity.
(2) The high-order deformation modes will be activated when . Then, the deformation of the stiffened plate will increase and the dynamic ultimate load capacity of the stiffened plate will decrease.
(3) The static ultimate load capacity after the stiffened plate is subjected to an impact load whose amplitude is the dynamic ultimate load capacity is reduced by less than 3%.
(4) The failure mode of the stiffened plate will change as the impact duration changes. When , the dynamic failure mode of the stiffened plate is the same as the static failure mode.
(5) The dynamic ultimate load capacity of the stiffened plate decreases as the number of impact cycles increases, while the impact duration remains constant.
(6) When the impact load with the highest amplitude is the last in the load sequence, the dynamic ultimate load capacity of the stiffened plate is the smallest; when the impact load with the highest amplitude is first in the load sequence, the dynamic ultimate load capacity of the stiffened plate is the greatest.