Assessment of Ship Hull Ultimate Strength under Fire Conditions: The Fire Smith Method Approach

<jats:p>This article presents the fire Smith method, which integrates the plate model and the Smith method, to analyze the impact of fire temperature and extent on the buckling strength of frigates. This investigation focused on the frigate’s buckling strength and how it is affected by a fire, achieved by modifying both deck temperature and area. Leveraging the principle of oxygen consumption, a coefficient “D” was introduced to account for fire temperature and region variations. This enabled the characterization of buckling strength under varying temperature and regional conditions through 64 simulations conducted on a frigate model. The outcomes revealed that the disparity between simulation results and the fire Smith method remained below 10%, establishing a solid basis for engineering assessment. As the high-temperature area decreased and necessitated less oxygen, the ultimate strength initiated a decline. However, upon reaching a certain threshold temperature, the ultimate strength stabilized. Conversely, expansive high-temperature zones caused a decline in ultimate strength, coupled with an increased oxygen requirement. Under consistent oxygen consumption conditions, the rate of ultimate strength reduction intensified. Consequently, these divergent characteristics of ultimate strength in various high-temperature areas underscored that minimizing the expansive fire high-temperature zones can significantly enhance the ship’s fire resistance, safeguarding its structural integrity.</jats:p>

section_content: Ships play an important role in the transportation of goods and upholding national security. In recent years, ship accidents have caused serious economic and property losses. Between 2009 and 2018, 48 bulk carriers with a dead weight exceeding 10,000 tons suffered various accidents that caused more than 100 fatalities [1]. Extensive research efforts have been directed towards enhancing hull safety and reliability, encompassing studies on intact hull ultimate strength [2, 3], damages [4] [5] [6] [7] [8] [9], combined loads [10, 11], and ice-related scenarios [12, 13]. Although the safety of ships has been improved in recent decades, more detailed and advanced research should be considered to ensure that the ultimate strength of ships can withstand loads under various conditions, including fire. 

In recent years, a series of significant ship fires have underscored the criticality of fire safety. Notably, on 12 July 2020, the "USS Bonhomme Richard" of the 7th Fleet of the US Army caught fire while undergoing maintenance work at the dock of the San Diego Naval Base [14]. Similarly, on 17 July 2020, a fire broke out during the repair of the amphibious assault ship "Kearsarge" located at Norfolk [15]. On the evening of 29 March 2023, a passenger ship caught fire near the island of Brasiland in the Philippines, resulting in a devastating loss of at least 29 lives [16]. The frequent occurrence of fire accidents has raised concerns among scholars, who have conducted research on the ultimate strength of ship hulls under fire conditions. Liu [17] has calculated the fire temperature field based on the large eddy model and used a self-developed temperature mapping program to map the temperature load to the structure for thermal coupling calculation. Nassiraei [18] studied the influence of joint geometry and ring plate size on initial stiffness, ultimate strength, and failure mechanism at different high temperatures. At high temperatures, the initial stiffness and ultimate strength of the X-joint will increase significantly when the ring plate is used. Through nonlinear regression analysis, a parameter equation was proposed to determine the high-temperature compressive support effectiveness of X-joints reinforced with ring plates. Nassiraei [19] studied the influence of ring plate size and joint geometry on the static performance of X-joint reinforced by ring plates under high-temperature tensile loads. The failure mechanism of a joint can be significantly improved by using larger ring plates. It was noted that the tensile strength of ringplate-reinforced X-joints diminishes as temperature increases, a trend attributed to the decrease in the elastic modulus of the steel material. Li [20] studied the effect of engine room fires on the ultimate strength of the ship using thermal radiation flux values as boundary conditions. Paik [21, 22] conducted stiffened plate fire tests to measure the ultimate strength under ISO-834 standard [23] conditions. This provides a reference for the failure modes of reinforced panels under fire conditions. Liu [24] studied the temperature distribution and fire response of multi-compartment fire, proposed the predictive model for hull ultimate strength under fire conditions, and obtained a positive correlation between compartment residual strength and the heat release rate. Guo [25] delved into the failure characteristics, stress distribution, and deformation modes of a two-layer cabin under different heat release rates and fire positions. The study predicted residual cabin strength under heat release rates ranging from 0 to 8 MW at different fire positions. 

In order to evaluate the ultimate strength of ship hulls, C. S. Smith proposed the Smith method [2]. Over time, numerous scholars have refined and adapted this method to suit different operational scenarios. Piscopo [26] summarizes the modified incremental iteration method tailored for bulk carriers under alternating loading conditions. Taking ISSC bulk carrier as a reference ship and comparing it with the finite element results from the literature, the effectiveness of this method was verified. Zhao's [27] investigation into large openings in cabins utilized experimental and simulation approaches. By validating the experimental and simulation outcomes through the Smith method, this study highlighted its practical feasibility. Li [28] compared the ultimate strength of a stranded hull by proposing the damage index method and simulation calculation, providing a reference and overview for future research. Xu [29] studied the buckling ultimate strength and proposed an analytical solution to the evaluation method of the hull ultimate strength under extreme wave loads. Zhu [30] studied the residual strength in damaged ships with the vertical and horizontal hull girder loads, as well as under combined stresses. The improved Smith method was used to calculate the residual load of the damaged ship, and the particle swarm optimization method was used to obtain the neutral axis of the damaged section. Putranto's [31] innovative proposition of an equivalent single-layer approach for the hull girder ultimate strength not only enhanced computational speed but also streamlined calculations. Cerik [32] introduced a fast procedure to determine the translation and rotation of the neutral axis at each curvature increment, which was used for a rapid assessment of the hull strength under severe conditions. Guedes Soares [33] studied the elastic-plastic behavior of plates under high temperature, which made a very important reference for the study of mechanical properties of plates at such conditions. Guedes Soares [34] delved into the elasticplastic behavior of plates at high temperatures, in order to comprehend mechanical properties in such conditions. Gordo [35] put forward a method to estimate the ultimate bending moment, considering the strength degradation caused by residual stress caused by welding high temperature, which provides an important reference for the element behavior under high-temperature conditions. Chen [36] proposed a simplified method to study the ultimate strength of welded stiffened plates, which provided a reference for the study of ultimate strength characteristics at high temperature. Teixeira and Guedes Soares [37] introduced the reliability formula of a thermal-insulation-bearing plate under local thermal load. While it is evident that numerous scholars have innovatively proposed improvements to the Smith method for assessing ultimate strength under diverse conditions, it is noteworthy that the application of the fire Smith method remains relatively constrained despite its potential. A first-order second-moment method was proposed to quantify the uncertainty of thermal load and describe the importance of control variables in the limit state function. 

Considering the effect of fire on the ultimate strength of the hull at the beginning of the ship design can reduce the damage of fire to the overall ultimate strength of the hull by using fire-resistant coatings or combustible materials. It is helpful for firefighters to put out the fire in the position that has great influence on the ultimate strength, such as the hull deck. For rapid decision making, it is helpful to rationally allocate firefighters to extinguish fires and realize more effective resource allocation. 

The sinking of the Sheffield [38] after being hit by a flying fish missile is an example that highlights the importance of correctly evaluating the ultimate strength of the hull under fire conditions. This disaster could have been avoided if the impact of fire had been properly considered during the design process. The analysis of the use of aluminum alloy as hull material [39] suggests that the impact of fire may be underestimated during the design process. Using a fire simulation method to evaluate the impact of fire on the ultimate strength of the hull at the beginning of the design process can help to avoid such disasters. 

Analysis of the above research finds that the ultimate strength of ships under fire conditions mainly focuses on two methods: trials and simulations. For the trial method, it can reflect the real fire failure mode, but currently there is limited fire test data, making it difficult to apply to all ship types. For simulation methods, the temperature mapping software needs to be developed, the professional technical level of personnel evaluating the ultimate strength of ships under fire conditions is high, and the work cycle is long. Therefore, a method that supports rapid decision making by command and rescue personnel was needed. However, simulation and trial cannot achieve the goal of simple and rapid assessment of ultimate strength under fire conditions. In summary, to address the issue of simple and rapid assessment of the ultimate strength of ship hulls under fire conditions, a rapid assessment method was proposed.
section_name: Introduction

section_content: The key of the fire Smith method lies in the characteristics of element in the fire stressstrain relationship. With regard to the plate behavior under high temperature, relevant experiments and simulation calculations have been carried out in the past, providing an important reference for the fire Smith method. 

In previous studies, it can be found that the failure mode of steel structure gradually changes from elasticity to plasticity in the process of compression and temperature rise, as shown in Figure 1. Yang's [40, 41] exploration of ultimate strength and failure modes in I-beams at diverse temperatures provided essential insights into the progressive shift from elastic to plastic behavior within steel structures under high temperatures. 

Han [42] made a comprehensive study on the local buckling behavior of an aluminum alloy plate in fire and obtained the typical failure mode under compression, as shown in Figure 2. It can be observed that buckling occurred in the middle of the plate. 

Moreover, Han et al. [43] extended their study to scrutinize the local buckling behavior of high-strength aluminum alloy columns in fire, obtaining the failure mode and ultimate strength curve and using it to evaluate the accuracy of the current design method of aluminum alloy short columns in fire. As shown in Figure 3, it can be found that the failure modes of plates at 20 • C to 400 • C were similar, all of which were buckling failure. Han [42] made a comprehensive study on the local buckling behavior of an aluminum alloy plate in fire and obtained the typical failure mode under compression, as shown in Figure 2. It can be observed that buckling occurred in the middle of the plate. Moreover, Han et al. [43] extended their study to scrutinize the local buckling behavior of high-strength aluminum alloy columns in fire, obtaining the failure mode and ultimate strength curve and using it to evaluate the accuracy of the current design method of aluminum alloy short columns in fire. As shown in Figure 3, it can be found that the failure modes of plates at 20 °C to 400 °C were similar, all of which were buckling failure. Presently, there are few buckling tests on the real ship fire. It is conceivable that there may be an interaction between the plate in the real ship fire, which may be different from Han [42] made a comprehensive study on the local buckling behavior of an aluminum alloy plate in fire and obtained the typical failure mode under compression, as shown in Figure 2. It can be observed that buckling occurred in the middle of the plate. Moreover, Han et al. [43] extended their study to scrutinize the local buckling behavior of high-strength aluminum alloy columns in fire, obtaining the failure mode and ultimate strength curve and using it to evaluate the accuracy of the current design method of aluminum alloy short columns in fire. As shown in Figure 3, it can be found that the failure modes of plates at 20 °C to 400 °C were similar, all of which were buckling failure. Presently, there are few buckling tests on the real ship fire. It is conceivable that there may be an interaction between the plate in the real ship fire, which may be different from Han [42] made a comprehensive study on the local buckling behavior of an aluminum alloy plate in fire and obtained the typical failure mode under compression, as shown in Figure 2. It can be observed that buckling occurred in the middle of the plate. Moreover, Han et al. [43] extended their study to scrutinize the local buckling behavior of high-strength aluminum alloy columns in fire, obtaining the failure mode and ultimate strength curve and using it to evaluate the accuracy of the current design method of aluminum alloy short columns in fire. As shown in Figure 3, it can be found that the failure modes of plates at 20 °C to 400 °C were similar, all of which were buckling failure. Presently, there are few buckling tests on the real ship fire. It is conceivable that there may be an interaction between the plate in the real ship fire, which may be different from Presently, there are few buckling tests on the real ship fire. It is conceivable that there may be an interaction between the plate in the real ship fire, which may be different from the failure mode of a single beam or plate, but it can provide reference to some extent. As is well known, under normal temperature conditions, the hull deck is subjected to buckling failure due to compression [44]. If the failure mode of the hull deck under compression due to high temperature is plastic deformation, before the failure mode of the whole ship fire changes from elasticity to plasticity, if the elastic buckling strength is estimated by plastic yield strength, the ultimate strength of the whole ship in fire will be overestimated. 

Conversely, the buckling force of a plate under compression is less than that of plastic deformation. Therefore, using the elastic buckling strength for estimation could lead to an underestimation of the ultimate strength. However, there is a safety factor built into this approach. Although it may lead to a certain cost in the design process, these expenses are still considerably lower than those incurred due to fire-induced damage to the entire ship. 

In the actual fire event, the development speed is fast, and if the influence of fire is underestimated, it may lead to irreparable cost. Therefore, from an engineering perspective, it is highly important to consider the utilization of elastic buckling as a means to underestimate the ultimate strength on fire. 

Based on the theory of small deflection plates, an analytical solution for the elastic instability load of a single direction uniformly compressed simply supported plate under fire conditions was proposed. The governing equation for in-plane stressed thin rectangular isotropic plate is outlined [45]. The coordinate system of an equal thickness thin plate is shown in Figure 4, and the XY and XZ plane of the plate is the midplane of the plate. For a differential microelement (dx, dy) and (dx, dz) from the plate, the normal and shear stresses acting on each surface are shown in Figure 4. 

by plastic yield strength, the ultimate strength of the whole ship in fire will be overestimated. 

Conversely, the buckling force of a plate under compression is less than that of plastic deformation. Therefore, using the elastic buckling strength for estimation could lead to an underestimation of the ultimate strength. However, there is a safety factor built into this approach. Although it may lead to a certain cost in the design process, these expenses are still considerably lower than those incurred due to fire-induced damage to the entire ship. 

In the actual fire event, the development speed is fast, and if the influence of fire is underestimated, it may lead to irreparable cost. Therefore, from an engineering perspective, it is highly important to consider the utilization of elastic buckling as a means to underestimate the ultimate strength on fire. 

Based on the theory of small deflection plates, an analytical solution for the elastic instability load of a single direction uniformly compressed simply supported plate under fire conditions was proposed. The governing equation for in-plane stressed thin rectangular isotropic plate is outlined [45]. The coordinate system of an equal thickness thin plate is shown in Figure 4, and the XY and XZ plane of the plate is the midplane of the plate. For a differential microelement (dx, dy) and (dx,dz) from the plate, the normal and shear stresses acting on each surface are shown in Figure 4.
section_name: Materials and Methods

section_content: Three fundamental assumptions underpin the small deflection model: (1) The normal strain perpendicular to the midplane is extremely small and can be disregarded. (2) Stress components σz, τzx, and τzy are smaller than σx, σy, and τxy; thus, normal strain εz and shear strain γzx and γzy can be omitted. The bending of thin plates with small deflections can be simplified as a plane stress problem. (3) When thin plate is bent, there is no strain parallel to the midplane at any point in the midplane. 

When solving the small deflection bending problem of thin plates based on displacement, the deflection of thin plates ω is established. An elastic surface an elastic surface differential equation is formulated a small deflection theoretical plate based on geometric, physical, and force equilibrium equations as an unknown function,
section_name: Uniform Compression Thin Plate Model with Four Simply Supported Sides

section_content: Three fundamental assumptions underpin the small deflection model: (1) The normal strain perpendicular to the midplane is extremely small and can be disregarded. (2) Stress components σz, τzx, and τzy are smaller than σx, σy, and τxy; thus, normal strain εz and shear strain γzx and γzy can be omitted. The bending of thin plates with small deflections can be simplified as a plane stress problem. (3) When thin plate is bent, there is no strain parallel to the midplane at any point in the midplane. 

When solving the small deflection bending problem of thin plates based on displacement, the deflection of thin plates ω is established. An elastic surface an elastic surface differential equation is formulated a small deflection theoretical plate based on geometric, physical, and force equilibrium equations as an unknown function,

In the Equation (1), D is the stiffness coefficient of the plate, D = Et 3

12(1-µ 2 ) ; E is the Elastic modulus; T is the thickness; µ is Poisson's ratio; N x1, N y is the buckling load on the midplane of the plate along the x and y axes; and N xy is the shear force per unit length on the midplane of the plate. 

For simple supported plate N y = N xy = 0 under unidirectional uniform compression, then

In Equation ( 3 ), m and n are the half wave numbers in the x and y directions, respectively, both being positive integers and are positive integers, and A mn is undetermined constant. Substituting Equation (3) into Equation (2) yields Equation ( 4 ):

If Equation (4) = 0, each coefficient of the unique condition that Equation ( 4 ) is 0, that is, the instability condition of the plate is

Let k be the buckling coefficient:

In a simple supported plate under uniform unidirectional compression, the Buckling load is the minimum load for the plate to maintain a slightly curved state, assumption n = 1. 

For a simply supported plate subjected to uniform unidirectional compression, the buckling load represents the minimum load required to maintain a slightly curved state. Assuming n = 1, and recognizing that dk/dm = 0, the expression simplifies:

This establishes the smallest buckling coefficient when β = m, and subsequently, k min = 4; therefore,

It can be found that under the condition that the geometric dimensions a and b and the buckling half wave number m and n are determined, the buckling load Nx of the plate varies with the elastic modulus of the plate, then Nx is a function of the elastic modulus E.
section_name: Uniform Compression Thin Plate Model with Four Simply Supported Sides

section_content: We assume that the non-uniform compression plate is also a small deformation plate. Due to the longitudinal force characteristics of the ship, the force in the area further away from the neutral axis increases, while the force in the area closer to the neutral axis decreases. Therefore, the force characteristics of the reinforcing ribs are non-uniform compression. This is shown in Figure 4, where a is the length c and the width. Under the longitudinal load of the ship, the cross-sectional stress of the stiffener along the z-axis direction is linearly distributed. The maximum stress value at the compression edge of the stiffener is σ 1, and the stress at the bottom of the stiffener under compression is σ 2. We define stress gradient coefficient as α 0 = (σ 1σ 2 )/σ 1. Then, the stress at z from the deck σ can be represented as σ = σ 1 (1αz/c). 

We applied simply supported boundary conditions to the loading edge of the longitudinal deck girders and the edge in contact with the deck, while leaving the edge in contact with the flange as a free edge. For the buckling coefficient of a thin plate under non-uniform compression, with three simply supported edges and one free edge, according to the Technical Code of Cold-Formed Thin-Wall Steel Structures [46], the calculation formula is as in Equation ( 10 ):

The European EC3 standard [47] also provides a calculation method for the buckling coefficient of non-uniform compression plates as in Equation ( 11 ):

According to the stress characteristics of the total longitudinal strength of the hull, the area farther away from the neutral axis is more stressed, and the length of the stiffened plate is 0. 02 m. Assuming that the neutral axis is in the center of the ship, the stress at the bottom of the stiffened plate is increased α 2 = α 1 (1 -0. 02/(H/2)) = 0. 985α 1. Then, the stress gradient coefficient α 0 = (σ 1σ 2 )/σ 1 = (1 -0. 985)/1 = 0. 015. According to Equation (10), the buckling coefficient k 2 = 0. 564, and therefore
section_name: Non-Uniform Compression Thin Plate Model with Three Simply Supported Sides and One Free Side

section_content: In terms of the stress characteristics of the flange, while it experiences equal compression as the deck in Section 2. 1, both sides of the flange remain free as they are not supported. Moreover, due to its smaller width-thickness ratio compared to the deck, the thin plate theory is not applicable, and a beam buckling model is more suitable. Therefore, considering its size, constraint form, and stress distribution, the flange is treated as a beam to account for the buckling force. According to Euler's formula,

When the flange width is e, the plate thickness is t, and the flange width is greater than the plate thickness, Equation ( 13 ) can be written as Equation ( 14 ):
section_name: Uniform Compression Beam Model with Two Simply Supported Sides and Two Free Sides

section_content: The hull of Figure 4 comprises H36 steel, characterized by a Young's modulus of 206 GPa and a Poisson's ratio of 0. 3 under normal temperature conditions. To model the evolution of material properties with temperature, one can refer to the recommended data of the American Institute of Steel Construction Specification Committee (AISC) [48]. Figure 5 illustrates the elastic modulus along with its fitting curve, in line with AISC data.
section_name: Ultimate Strength of Uniform and Non-Uniform Compression Plate Models under High-Temperature Conditions

section_content: The hull of Figure 4 comprises H36 steel, characterized by a Young's modulus of 206 GPa and a Poisson's ratio of 0. 3 under normal temperature conditions. To model the evolution of material properties with temperature, one can refer to the recommended data of the American Institute of Steel Construction Specification Committee (AISC) [48]. Figure 5 illustrates the elastic modulus along with its fitting curve, in line with AISC data. According to the curve of the elastic modulus measured by AISC as a function of temperature, the fitting equations for elastic modulus in different temperature ranges can be established: According to the curve of the elastic modulus measured by AISC as a function of temperature, the fitting equations for elastic modulus in different temperature ranges can be established:
section_name: Ultimate Strength of Uniform and Non-Uniform Compression Plate Models under High-Temperature Conditions

section_content: -0. 257T + 252. 87 -0. 385T + 328. 12 -0. 054T + 65. 28 0 < T ≤ 200 200 < T ≤ 600 600 < T ≤ 800 800 < T ≤ 1200 (15) According to Equation (15) and stiffness coefficient of the plate D = Et 3 12(1-µ 2 ) N x1,cr,min as shown in Equation ( 16 ):

Uniform buckling stress σ x1,cr,min is shown in Equations ( 17 ) and ( 18 ):

Similarly to Equation ( 16 ), Nx 2,cr is as shown in Equation ( 19 ):

Non-uniform buckling stress σ x2,cr is shown in Equations ( 20 ) and ( 21 ):

The formula for buckling stress for flange is shown in Equations ( 22 ) and ( 23 ):

C. S. Smith proposed the Smith method to discretize the hull into multiple elements and assumed that these elements did not affect each other [2]. Although the basic structure of the hull is stiffened plates, stiffened plates can be discretized into multiple plate elements, as shown in Figure 6. The key of Smith's method is to obtain the correct stress-strain curve of the element. Taking Reference [43] as an example, it can be found that the force and displacement curve of the plate under pressure at high temperature in Figure 7 shows an approximate linear relationship between force and displacement before the buckling force reaches its maximum value. When the maximum is reached, it approximately presents a "flat area" phenomenon, and with the increase in temperature, the gentler the force is after reaching the maximum. Based on the characteristics of force and displacement of the elements under various temperature conditions in the experiment, to obtain the maximum buckling stress, we utilize Equations ( 24 )-( 26 ) under the assumption of linear stress-strain behavior before reaching the maximum buckling stress. Additionally, after reaching the maximum stress, The key of Smith's method is to obtain the correct stress-strain curve of the element. Taking Reference [43] as an example, it can be found that the force and displacement curve of the plate under pressure at high temperature in Figure 7 shows an approximate linear relationship between force and displacement before the buckling force reaches its maximum value. When the maximum is reached, it approximately presents a "flat area" phenomenon, and with the increase in temperature, the gentler the force is after reaching the maximum. The key of Smith's method is to obtain the correct stress-strain curve of the element. Taking Reference [43] as an example, it can be found that the force and displacement curve of the plate under pressure at high temperature in Figure 7 shows an approximate linear relationship between force and displacement before the buckling force reaches its maximum value. When the maximum is reached, it approximately presents a "flat area" phenomenon, and with the increase in temperature, the gentler the force is after reaching the maximum. Based on the characteristics of force and displacement of the elements under various temperature conditions in the experiment, to obtain the maximum buckling stress, we utilize Equations ( 24 )-( 26 ) under the assumption of linear stress-strain behavior before reaching the maximum buckling stress. Additionally, after reaching the maximum stress, Based on the characteristics of force and displacement of the elements under various temperature conditions in the experiment, to obtain the maximum buckling stress, we utilize Equations ( 24 )-( 26 ) under the assumption of linear stress-strain behavior before reaching the maximum buckling stress. Additionally, after reaching the maximum stress, keep the stress constant while performing calculations at a very small strain ∆ε. 

The operational sequence of the fire Smith method is depicted in Figure 8. First, given the value of element number i, a temperature matrix is constructed according to the temperature distribution. According to Equations ( 24 )-( 26 ), the stress calculation formula of which temperature range the current element is applicable to is determined, and the stress value corresponding to the current temperature is calculated. Assuming the buckling strain value is 0. 005, the stress-strain curve of the current temperature is obtained. The process then proceeds to initialize the curvature and incrementally raise the curvature value. As the curvature value increases, the concurrent strain increment is calculated, in turn determining the current stress value as per the stress-strain relationship, and then calculating the magnitude of the force. The moment is assessed at the neutral axis, followed by a reevaluation of the neutral axis position. Upon detecting the amalgamated force, neutral axis position, curvature magnitude, and unit number that aligns with the predefined criteria, the bending moment value for the unit is output.
section_name: 206

section_content: According to the principle of oxygen consumption, alkane combustibles consume 1 kg of oxygen to release heat, which is approximately 13. 6 MJ/kg [49]. Therefore, it is assumed that under the condition of a closed and insulated compartment, the sum of heat received by the top of an oxygen compartment that consumes the same mass is equal. The
section_name: Simulation Setup

section_content: According to the principle of oxygen consumption, alkane combustibles consume 1 kg of oxygen to release heat, which is approximately 13. 6 MJ/kg [49]. Therefore, it is assumed that under the condition of a closed and insulated compartment, the sum of heat received by the top of an oxygen compartment that consumes the same mass is equal. The purpose of working condition design is to study the ultimate strength characteristics and effects of ship hulls in different regions with the same heat distribution. 

Using the ABAQUS thermal coupling method, the influence of distinct top decks being exposed to oil pool fires is calculated. For smaller-scale oil pools, the fire source area is smaller, the combustion time is longer, and the ignition area is limited, resulting in extended combustion durations, concentrated temperatures, and higher temperature levels. In contrast, larger oil pools exhibit expanded ignition areas, shorter combustion times, dispersed temperatures, and lower temperature levels. To analyze the ultimate strength of the ship under different heating conditions and verify the feasibility of the Smith method under fire conditions, 64 sets of ultimate strength calculations under fire conditions were carried out. We assumed an oxygen consumption coefficient denoted as D; the heating zone width as W as shown in Figure 9 ; and the temperature T, which can be represented as D = W × T. This approach facilitates the exploration of ultimate strength attenuation patterns in ship hulls under differing heating zone widths W and temperatures T. The working condition settings are detailed in Table 1. Figure 9a shows a small-scale oil pool fire. Assuming a heated area length above the fire source of W = 1 m, the ultimate strength characteristics of the hull at different temperatures were studied. Figure 9b shows a medium-sized oil pool fire scenario with the heated area above the fire source at W = 2 m. Under conditions of equivalent oxygen consumption, the temperature of heating the deck with the same heat was half that of W = 1 m. By analogy, Figure 9c,d shows large and larger oil pool fires, assuming that the deck temperatures of W = 3 m and W = 4 m were 1/3 and 1/4 of W = 1 m, respectively, under the same oxygen consumption. The calculation of hull ultimate strength at high temperature was divided into two steps: thermal-mechanical coupling calculation and explicit dynamics calculation. Coupled temp-displacement load step was used. The boundary condition definition is shown in Figure 10. The temperature load was applied above the deck using simply supported boundary conditions at both ends. The element type was S4RT, and the amplitude curves used to define the temperature followed the equally spaced method. 

m. Under conditions of equivalent oxygen consumption, the temperature of heati deck with the same heat was half that of W = 1 m. By analogy, Figure 9c, d The calculation of hull ultimate strength at high temperature was divided into two steps: thermal-mechanical coupling calculation and explicit dynamics calculation. Coupled temp-displacement load step was used. The boundary condition definition is shown in Figure 10. The temperature load was applied above the deck using simply supported boundary conditions at both ends. The element type was S4RT, and the amplitude curves used to define the temperature followed the equally spaced method. After completing the thermal-mechanical coupling calculation, in order to consider the structural deformation caused by temperature stress, the deformed structure was used as the model output for explicit dynamic analysis, and a predefined field method was used to apply the temperature field to the structure. In explicit dynamic analysis, the elements type is S4RT, and the boundary conditions are shown in Figure 11. One end was constrained by a simple support, while the other end was constrained by a coupling method that constrains movement in the Y and Z directions, as well as rotation in the X and Z directions, and applies rotation around the Y axis to the entire plane with a rotation angle of 0. 2°. The ultimate high-temperature strength was obtained by outputting the After completing the thermal-mechanical coupling calculation, in order to consider the structural deformation caused by temperature stress, the deformed structure was used as the model output for explicit dynamic analysis, and a predefined field method was used to apply the temperature field to the structure. In explicit dynamic analysis, the elements type is S4RT, and the boundary conditions are shown in Figure 11. One end was constrained by a simple support, while the other end was constrained by a coupling method that constrains movement in the Y and Z directions, as well as rotation in the X and Z directions, and applies rotation around the Y axis to the entire plane with a rotation angle of 0. 2 •. The ultimate high-temperature strength was obtained by outputting the bending moment value.
section_name: Simulation Setup

section_content: The mesh employed in a finite element model significantly influences the analysis outcomes. With the increase in the number of meshes, the results generally converge. Generally speaking, with the increase in the number of meshes, the calculation results gradually converge but at the same time, and the calculation amount increases, so it is important to correctly grasp the number of meshes and the accuracy of the calculation results. In order to verify the mesh sensitivity of the ultimate strength value of frigates, for the mesh sizes of 0. 08 m, 0. 04 m, 0. 02 m, and 0. 01 m, the total number of meshes were 6058, 15,575, 55,300, and 207,700 respectively. As shown in Figure 12, the ultimate strength of frigates with different mesh sizes was calculated, and it was found that the ultimate strength decreased with the increase in the number of meshes, and the ultimate strengths of different meshes were 3. 465 × 10 7 N•m, 2. 918 × 10 7 N•m, 2. 879 × 10 7 N•m, and 2. 873 × 10 7 N•m, respectively. In order to ensure the accuracy of the calculation results and consider the calculation speed, the frigate's mesh configuration is illustrated in Figure 13, featuring both an overall
section_name: Mesh Sensitivity Analysis

section_content: The mesh employed in a finite element model significantly influences the analysis outcomes. With the increase in the number of meshes, the results generally converge. Generally speaking, with the increase in the number of meshes, the calculation results gradually converge but at the same time, and the calculation amount increases, so it is important to correctly grasp the number of meshes and the accuracy of the calculation results. In order to verify the mesh sensitivity of the ultimate strength value of frigates, for the mesh sizes of 0. 08 m, 0. 04 m, 0. 02 m, and 0. 01 m, the total number of meshes were 6058, 15,575, 55,300, and 207,700 respectively. As shown in Figure 12, the ultimate strength of frigates with different mesh sizes was calculated, and it was found that the ultimate strength decreased with the increase in the number of meshes, and the ultimate strengths of different meshes were 3. 465 × 10 7 N•m, 2. 918 × 10 7 N•m, 2. 879 × 10 7 N•m, and 2. 873 × 10 7 N•m, respectively.
section_name: Mesh Sensitivity Analysis

section_content: The mesh employed in a finite element model significantly influences the analysis outcomes. With the increase in the number of meshes, the results generally converge. Generally speaking, with the increase in the number of meshes, the calculation results gradually converge but at the same time, and the calculation amount increases, so it is important to correctly grasp the number of meshes and the accuracy of the calculation results. In order to verify the mesh sensitivity of the ultimate strength value of frigates, for the mesh sizes of 0. 08 m, 0. 04 m, 0. 02 m, and 0. 01 m, the total number of meshes were 6058, 15,575, 55,300, and 207,700 respectively. As shown in Figure 12, the ultimate strength of frigates with different mesh sizes was calculated, and it was found that the ultimate strength decreased with the increase in the number of meshes, and the ultimate strengths of different meshes were 3. 465 × 10 7 N•m, 2. 918 × 10 7 N•m, 2. 879 × 10 7 N•m, and 2. 873 × 10 7 N•m, respectively. In order to ensure the accuracy of the calculation results and consider the calculation speed, the frigate's mesh configuration is illustrated in Figure 13, featuring both an overall In order to ensure the accuracy of the calculation results and consider the calculation speed, the frigate's mesh configuration is illustrated in Figure 13, featuring both an overall view and an enlarged section. The overall mesh size of the frigate was 0. 02 m, and the mesh type was square with a total of 55,300 elements. 

view and an enlarged section. The overall mesh size of the frigate was 0. 02 m, and the mesh type was square with a total of 55,300 elements.
section_name: Mesh Sensitivity Analysis

section_content: Figure 14 shows typical failure modes under different high temperature conditions. The deformation of the top deck due to thermal expansion at high temperatures resulted in a "wave"-shaped failure pattern. Figure 14a is the failure mode of ship hull under the fire condition of a small-scale oil pool, which corresponds to Figure 9a. It can be found that the top high-temperature area presented a "wave"-shaped deformation under the conditions of fire load and compression load. Figure 14b is the ultimate state of ship fire under the condition of a medium-sized oil pool, which corresponds to Figure 9b. Compared with Figure 14a,b, it can be found that with the increase in the high-temperature area, the "wave-like" deformation area was larger, being caused by the larger structural deformation area caused by high fire temperature, which leads to a more decline in the ultimate strength of a ship. Similarly, Figure 14c,d shows the hull failure mode under the conditions of a large-scale oil pool fire and larger-scale oil pool fire. It can be found that due to the high temperature, the "wave-like" deformation gradually spread to the whole deck, which led to a further decrease in ultimate strength.
section_name: Results

section_content: Figure 14 shows typical failure modes under different high temperature conditions. The deformation of the top deck due to thermal expansion at high temperatures resulted in a "wave"-shaped failure pattern. Figure 14a is the failure mode of ship hull under the fire condition of a small-scale oil pool, which corresponds to Figure 9a. It can be found that the top high-temperature area presented a "wave"-shaped deformation under the conditions of fire load and compression load. Figure 14b is the ultimate state of ship fire under the condition of a medium-sized oil pool, which corresponds to Figure 9b. Compared with Figure 14a,b, it can be found that with the increase in the high-temperature area, the "wave-like" deformation area was larger, being caused by the larger structural deformation area caused by high fire temperature, which leads to a more decline in the ultimate strength of a ship. Similarly, Figure 14c,d shows the hull failure mode under the conditions of a large-scale oil pool fire and larger-scale oil pool fire. It can be found that due to the high temperature, the "wave-like" deformation gradually spread to the whole deck, which led to a further decrease in ultimate strength. view and an enlarged section. The overall mesh size of the frigate was 0. 02 m, and the mesh type was square with a total of 55,300 elements.
section_name: Results

section_content: Figure 14 shows typical failure modes under different high temperature conditions. The deformation of the top deck due to thermal expansion at high temperatures resulted in a "wave"-shaped failure pattern. Figure 14a is the failure mode of ship hull under the fire condition of a small-scale oil pool, which corresponds to Figure 9a. It can be found that the top high-temperature area presented a "wave"-shaped deformation under the conditions of fire load and compression load. Figure 14b is the ultimate state of ship fire under the condition of a medium-sized oil pool, which corresponds to Figure 9b. Compared with Figure 14a,b, it can be found that with the increase in the high-temperature area, the "wave-like" deformation area was larger, being caused by the larger structural deformation area caused by high fire temperature, which leads to a more decline in the ultimate strength of a ship. Similarly, Figure 14c,d shows the hull failure mode under the conditions of a large-scale oil pool fire and larger-scale oil pool fire. It can be found that due to the high temperature, the "wave-like" deformation gradually spread to the whole deck, which led to a further decrease in ultimate strength. Figure 15 shows the ultimate strength of high-temperature regions with lengths of 1 m, 2 m, 3 m, and 4 m under different temperature conditions. It can be seen that with the increase in temperature, the ultimate strength showed a downward trend. When the length of the high-temperature area was 1 m, the ultimate strength of the small-scale oil pool fire in Figure 9a decreased rapidly with the temperature, and then tended to be flat, as shown in Figure 15a. Figure 15b shows the change of the ultimate strength of the large-scale oil pool fire in Figure 9b with temperature. Notably, the onset of decline in the ultimate strength was observed at around 200 • C, marked by an accelerated rate of decrease. Similar to Figure 15c, d, it showed large-scale and larger-scale oil pool fires, that is, when the length of the heating zone was 4 m, the ultimate strength needed more oxygen to start to decline, and it showed a rapid downward trend with the temperature change. Figure 15 shows the ultimate strength of high-temperature regions with lengths of 1 m, 2 m, 3 m, and 4 m under different temperature conditions. It can be seen that with the increase in temperature, the ultimate strength showed a downward trend. When the length of the high-temperature area was 1 m, the ultimate strength of the small-scale oil pool fire in Figure 9a decreased rapidly with the temperature, and then tended be flat, as shown Figure 15a. Figure 15b shows the change of the ultimate strength of the largescale oil pool fire in Figure 9b with temperature. Notably, the onset of decline in the ultimate strength was observed at around 200 °C, marked by an accelerated rate of decrease. Similar to Figure 15c, d, it showed large-scale and larger-scale oil pool fires, that is, when the length of the heating zone was 4 m, the ultimate strength needed more oxygen to start to decline, and it showed a rapid downward trend with the temperature change. To further analyze the ultimate strength characteristics of ships under fire conditions, the maximum ultimate strength at different temperatures was plotted in Figure 16. The fire Smith method was used to calculate and simulate the same temperature conditions as shown in the Figure 16. By comparing the maximum value of the ultimate strength of each working condition calculated by the fire Smith method and the simulation in Figure 16 and Table 2, it was found that the error between the simulation calculation and the fire Smith method was about 9. 5%. The sources of the maximum error value point were the W1 simulation calculation value and the fire smith method when T × W = 16. The reason was that with the increase in temperature, the influence of plasticity on ultimate strength increased gradually, while the smith result of fire with elastic buckling failure was less than that of the simulation calculation. It can provide a reference for the assessment of the ultimate buckling strength of a fire hull. To further analyze the ultimate strength characteristics of ships under fire conditions, the maximum ultimate strength at different temperatures was plotted in Figure 16. The fire Smith method was used to calculate and simulate the same temperature conditions as shown in the Figure 16. By comparing the maximum value of the ultimate strength of each working condition calculated by the fire Smith method and the simulation in Figure 16 and Table 2, it was found that the error between the simulation calculation and the fire Smith method was about 9. 5%. The sources of the maximum error value point were the W 1 simulation calculation value and the fire smith method when T × W = 16. The reason was that with the increase in temperature, the influence of plasticity on ultimate strength increased gradually, while the smith result of fire with elastic buckling failure was less than that of the simulation calculation. It can provide a reference for the assessment of the ultimate buckling strength of a fire hull. As shown in Figure 17, a comparison of maximum ultimate strength values for different high-temperature regions revealed several key characteristics: (1) The ultimate strength of the ship remained relatively stable until the temperature reached 200 • C. For the heating zone W 1 = 1 m, less oxygen was required to reach 200 • C, so less oxygen was required to reduce the ultimate strength. For the case where the heating area W 4 = 4 m, due to the larger heating area, more oxygen was required to reach 200 • C in the entire area, resulting in a decrease in ultimate strength requiring more oxygen. A comparison between two sets of working conditions (T × W = 8) revealed that the ultimate strength for W 1 was 2. 88 × 10 7 N•m, whereas W 4 exhibited an ultimate strength of 2. 57 ×10 7 N•m, resulting in an approximately 11% decrease for W 1. (2) T × W = 12, the heating area W 1 = 1 m, and the ultimate strength remained almost unchanged after reaching 1200 • C. However, when the heating area W 4 = 4 m, the temperature reached 300 • C, and the ultimate strength was still in a downward trend. (3) For T × W = 16, comparing the four different heating regions, it can be found that W 1 had a significant difference in ultimate strength compared to the other three groups due to reaching the temperature threshold. The downward trend of W 2, W 3, and W 4 was similar, with the greater ultimate strength attenuation associated with larger distribution areas. For the same T × W condition, a higher distribution area translated to a more pronounced decrease in ultimate strength. As shown in Figure 17, a comparison of maximum ultimate strength values for different high-temperature regions revealed several key characteristics: (1) The ultimate strength of the ship remained relatively stable until the temperature reached 200 °C. For the heating zone W1 = 1 m, less oxygen was required to reach 200 °C, so less oxygen was required to reduce the ultimate strength. For the case where the heating area W4 = 4 m, due to the larger heating area, more oxygen was required to reach 200 °C in the entire area, resulting in a decrease in ultimate strength requiring more oxygen. A comparison between two sets of working conditions (T × W = 8) revealed that the ultimate strength for W1 was 2. 88 × 10 7 N•m, whereas W 4 exhibited an ultimate strength of 2. 57 ×10 7 N•m, resulting in an approximately 11% decrease for W1. (2) T × W = 12, the heating area W1 = 1 m, and the ultimate strength remained almost unchanged after reaching 1200 °C. However, when the heating area W4 = 4 m, the temperature reached 300 °C, and the ultimate strength was still in a downward trend. (3) For T × W = 16, comparing the four different heating regions, it can be found that W1 had a significant difference in ultimate strength compared to the other three groups due to reaching the temperature threshold. The downward trend of W2, W3, and W4 was similar, the greater ultimate strength attenuation associated with larger distribution areas. For the same T × W condition, a higher distribution area translated to a more pronounced decrease in ultimate strength.
section_name: Results

section_content: On the basis of plate theory and the Smith method, the fire Smith method was proposed. Taking the frigate model as the research object, the simulation calculation was carried out, and it was found that the trend of simulation calculation was consistent with the fire Smith method, so the fire Smith method can be used as a rapid assessment method for hull ultimate strength under fire conditions in engineering. 

Based on the principle of oxygen consumption, the same heat is distributed in different high-temperature regions. Through 64 sets of simulation calculations, the ultimate strength characteristics of different high-temperature regions were found. The study showed that the ultimate strength began to decline with lesser oxygen consumption if the high-temperature region was smaller. Conversely, once the temperature reached the threshold, the ultimate strength remained relatively steady. For larger high-temperature regions, ultimate strength decreased with higher oxygen consumption and with a more pronounced attenuation. Comparing high-temperature regions W2, W3, and W4, it was evident that larger distribution areas led to greater ultimate strength reduction. Therefore,
section_name: Discussion

section_content: On the basis of plate theory and the Smith method, the fire Smith method was proposed. Taking the frigate model as the research object, the simulation calculation was carried out, and it was found that the trend of simulation calculation was consistent with the fire Smith method, so the fire Smith method can be used as a rapid assessment method for hull ultimate strength under fire conditions in engineering. 

Based on the principle of oxygen consumption, the same heat is distributed in different high-temperature regions. Through 64 sets of simulation calculations, the ultimate strength characteristics of different high-temperature regions were found. The study showed that the ultimate strength began to decline with lesser oxygen consumption if the high-temperature region was smaller. Conversely, once the temperature reached the threshold, the ultimate strength remained relatively steady. For larger high-temperature regions, ultimate strength decreased with higher oxygen consumption and with a more pronounced attenuation. Comparing high-temperature regions W 2, W 3, and W 4, it was evident that larger distribution areas led to greater ultimate strength reduction. Therefore, based on the ultimate strength characteristics of different high-temperature regions, controlling the fire scale to a smaller range when oxygen consumption is low can effectively reduce the damage of the fire to the ultimate strength of the ship. 

In previous studies, many scholars have used experimental and simulation methods to study the ultimate strength of ship hulls under fire conditions. They have found that the ultimate strength gradually decreases with the increase in fire temperature and fire power. However, the evaluation speed of experimental and simulation methods is relatively slow, and they cannot provide a rapid assessment of the ultimate strength of ship hulls during fire disasters. Although this article assumes that the buckling strain and minimum buckling stress of the plate under fire conditions lead to errors in the Smith method and simulation method, the ultimate strength of the ship decreased with the increase in fire temperature, which is consistent with previous studies. It is recommended to further explore failure modes and conduct ultimate strength tests under ship fire conditions to both supplement and validate these findings.
section_name: Discussion

section_content: In this study, the fire Smith method was put forward based on the Smith method and plate theory, and the frigate in different high-temperature areas was evaluated based on the oxygen consumption theory. Based on the simulation results and fire Smith method, the following conclusions can be drawn:

(1) In order to quickly evaluate the ultimate strength of ship buckling under fire conditions, a method was proposed based on the Smith method and the plate model. Because there may be a phenomenon that the hull gradually changes from elastic buckling to plasticity, in order to avoid overestimating the ultimate strength of the hull under fire conditions, the elastic buckling equation is used to estimate the ultimate strength value. Although there is a 10% error, it can to some extent avoid overestimating the ultimate strength of the hull before the elastic transition to plasticity, making it safer to evaluate the ultimate strength of the hull in engineering. 

(2) Based on the principle of oxygen consumption, 64 simulations were carried out to study the ultimate strength of different high-temperature regions, and it was found that the high-temperature regions showed a "wavy" failure mode. (3) According to the simulation results, the ultimate strength characteristics of different high-temperature regions were obtained. The smaller the high-temperature region, the less oxygen required for the ultimate strength to begin to decrease. At the same time, the ultimate strength remains basically unchanged after the temperature reaches the threshold. The larger the high-temperature area, the more oxygen is required for the ultimate strength to decrease, but as the oxygen consumption coefficient increases, the ultimate strength increases. (4) Comparing the fire Smith method with the simulation results, we found that the trend was consistent, so the proposed method can be used as an alternative approach for evaluating the ultimate strength of a fire hull.
section_name: Conclusions

section_name: Nomenclature

section_name: Institutional Review Board Statement: Not applicable.

section_content: Funding: This research received no external funding.

section_content: Informed Consent Statement: Not applicable.
section_name: Institutional Review Board Statement: Not applicable.

section_content: Data Availability Statement: Data sharing not applicable. No new data were created or analyzed in this study. Data sharing is not applicable to this article.

section_name: Conflicts of Interest:

section_content: Author Contributions: Conceptualization, J. W. and Z. D. ; methodology, J. W. ; software, J. W. ; validation, J. W. and Z. D. ; formal analysis, J. W. and M. Y. ; investigation, J. W. ; resources, Z. D. ; data curation, X. S. ; writing-original draft preparation, J. W. ; writing-review and editing, J. W. ; visualization, J. W. ; supervision, Z. D., M. Y. and X. S. All authors have read and agreed to the published version of the manuscript.

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: