1. Introduction
It is common for stiffened and unstiffened panels to be used for primary structural supporting members in the field of ocean and shore technology (OST), including shipping and offshore industries. In general, mild (MS24) and high tensile (AH32 or AH 36) steel grades are being used as construction materials for voyages in the Southern Sea Route (SSR). In the case of the Northern Sea Route (NSR), which was opened due to the global warming effect, the different grades of the steels, i.e., B, D, E, or F, are recommended to be used [1,2,3].
A wide range of studies in assessing and predicting the structural condition of intact and damaged structures have been conducted by many researchers for the robust design of ships and offshore structures. In particular, the finite element method (FEM), one of the famous numerical methods, is considered a powerful technique for solving issues in various fields, such as engineering and medical. Moreover, the computational fluid dynamics (CFD)-based and fluid-structure interaction (FSI)-based numerical simulations are also getting more and more popular in structural design, following the development of computer technology [4].
The experimental and analytical methods, meanwhile, are also considered useful ways to resolve engineering issues, provided that there is adequate financial support and time frame given, subject to physical limitation, such as equipment size and testing area issues. In the case of analytical methods, while it provides an absolute solution, it has the limitation of not covering complex geometries with actual environmental conditions [5].
In this regard, the design formulation or the empirical formulation approach is studied [6] and adopted by major classification societies, such as Lloyd’s Register (LR), the American Bureau of Shipping (ABS), and Det Norske Veritas Germanischer Lloyd (DNV GL), and followed by major shipyards in general. A number of empirical formulations and simplified techniques were proposed by various researchers in terms of ultimate limit state of intact and damaged hull girders [7,8,9,10,11,12,13,14], unstiffened panel (plate) [15,16,17], stiffened panel [6,18,19,20,21,22], and many others.
Recently, historical and technical reviews on existing empirical formulations in predicting the ultimate limit state (ULS) have been conducted for unstiffened and stiffened panels [23,24]. In the case of stiffened panel, which is made of steel and aluminium, most of the existing empirical formulations [6,21,25,26,27,28,29,30] were developed as a function of two parameters. such as plate slenderness ratioβand column slenderness ratio (λ), by adopting the plate stiffener combination (PSC) model.
More recently, Kim et al. [21] addressed that ultimate compressive strength of stiffened panel tends to fluctuate in the lower range of the column of slenderness ratio (λ ) as shown in Figure 1. This means that stiffened panel has a higher nonlinearity with the combination of plate and stiffener geometries. Therefore, simplified or single line shape empirical formulation, such as a function ofβandλ, may not be accurate in representing the ultimate strength behaviour of stiffened panel. Moreover, additional parameters should be considered in predicting the accurate ultimate compressive strength behaviour of stiffened panel.
In this regard, Kim et al. [6] introduced a refined empirical formulation in predicting the ultimate compressive strength of T-bar type stiffened panel, as shown in Figure 2, by using a conventional data processing technique with four (4) parameters, such as plate slenderness ratio (β), column slenderness ratio (λ), web slenderness ratio (hw/tw), and moment of inertia of stiffener to moment of inertia of plate ratio in z-direction (vertical) (Isz/Ipz ) From the wide range of numerical simulations (in total 10,500 scenarios of T-bar stiffened panel), an advanced empirical formulation by data processing technique has been proposed [6] with reliable accuracy of ULS, compared by ANSYS FE numerical simulation results (R2 = 0.98). Details may be referred to Kim et al. [6].
Most of the existing empirical formulations can be used in predicting the ULS of T-bar stiffened panel under longitudinal compression. However, there are limited studies conducted on flat- and angle-bar stiffened panel. In this regard, 540 numerical simulations were conducted by the ANSYS non-linear finite element method (NLFEM) to obtain the ultimate strength of flat-bar type stiffened panel under longitudinal compression. The ULS results of flat-bar stiffened panel were utilised as the input for the data processing. In addition, the applicability of the empirical formulation proposed by Kim et al. [6] was also tested as to whether it can be fitted for flat-bar or otherwise.
Finally, the accuracy of the refined empirical formulation for flat-bar stiffened panel obtained in this study was verified by statistical analysis. The applicability of the outcome obtained from this study was verified by ANSYS FE numerical simulations, as well as existing empirical formulations for flat-bar stiffened panel, by adopting single line shape formulations [21,25,26,28,29].
A useful outcome was achieved in predicting the ULS of flat-bar stiffened panel, which is one of the primary structural components of ships and offshore structures. 2. Brief Review of the Existing Formulations
As mentioned earlier, recently Zhang [23] and Kim et al. [24] provided detailed technical reviews on existing design and empirical formulations to predict the ultimate strength of the stiffened panel. The details of the same can be found in the articles mentioned above.
In this study, the representative existing empirical formulations are addressed in this section. The formulations introduced in this paper are used for comparison with FE numerical simulation by statistical analysis in the following section. In general, the empirical formulations in predicting the ultimate limit state (ULS) of the stiffened panel formulates as a function of plate slenderness ratio (β) and column slenderness ratio (λ), as shown in Equation (1).
σxuσYeq.=f(β, λ)
whereσxuis the ultimate compressive strength in x-axis (under longitudinal compression),σYeq.is the equivalent yield strength of plate and stiffener,βis the plate slenderness ratio (=(bp/tp)⋅(σYp/E)),λis the column slenderness ratio (=[L/(π⋅r)]⋅(σYeq./E)),σYp= is the yield strength of plate,Eis the Young’s modulus,Lis the length of the stiffened panel,ris the radius of gyration (=(I/A)), andIandAare the moment of inertia and sectional area for the plate-stiffener combination (PSC) model’s section, respectively.
By assuming the basic format of empirical formulations, illustrated in Equation (1), a number of empirical formulations (or closed form shape formulations) were proposed by experimental and numerical methods. The common design formulations used in a shipyard are presented in Equations (2a–c). Details on existing design formulations can be found in the research by Paik [16].
Euler formulation
σEσYeq.={1.0 for λ ≤1.01/λ2 for λ >1.0
Johnson–Ostenfeld (J–O) formulation
σcrσYeq.={ 1/λ2 for σE/σYeq. ≤0.51−λ2/4 for σE/σYeq. >0.5
Perry–Robertson (P–R) formulation under axial compression
σxuσYeq.=12(1+1+ηλ2)− 14(1+1+ηλ2)2−1λ2
Lin [25] also suggested the generalised shape of the empirical formulation to predict the ultimate strength of the stiffened panel subject to longitudinal compression, as shown in Equation (3a). Thereafter, Paik and Thayamballi [26] proposed the revised empirical formulation, based on the test database collected. In addition, they set the upper limit of the empirical formulation, of which the ultimate compressive strength of stiffened panel may not exceed the elastic buckling stress (σYeq./λ2 ), as shown in Equation (3b). The coefficients consisting of Equation (3a,b) are summarised in Table 1.
General shape of the empirical formulation
σxuσYeq.=1c1+c2 λ2+c3 β2+c4 λ2 β2+c5 λ4
where c1 – c5 are the coefficients, which can be found in Table 1.
Upper limit of empirical formulation
σxuσYeq.=1c1+c2 λ2+c3 β2+c4 λ2 β2+c5 λ4≤1λ2
More recent empirical formulations are also summarised in Equations (4)–(6). Zhang and Khan [28] proposed Equation (4) with a limitation of the range of column slenderness ratio (λ≤2 ). Kim et al. [21] also suggested that empirical formulation should be based on numerical simulations, which allows to cover the wide range of the column slenderness ratio. Recently, Xu et al. [29] proposed empirical formulation to predict the ULS of all types of stiffened panel, i.e., T-bar, angle-bar, and flat-bar subject to longitudinal compression, as well as lateral pressure. Among others, the empirical formulation for flat bar applied in this study is only presented in Equation (6).
Z-K formulation [28]
σxuσYeq.=1β0.281.0+λ3.2 for λ≤2 range only
Kim’s formulation [21]
σxuσYeq.=10.8884+eλ2+10.4121+eβ
Xu’s formulation [29] for flat bar
σxuσYeq.=1(1.127−4.915λ+0.49β+0.773λβ+10.075λ2−0.109β2−0.14λ2 β2−7.089λ3+0.04β3+0.01λ3 β3+1.564λ4)≤1λ2
Details of existing empirical formulations and their technical reviews can be found in the research by Zhang [23] and Kim et al. [24]. It is highlighted that the existing empirical formulations in Equations (2)–(6) are presented as a single line shape equation, and this is one of the important reasons why advanced empirical formulation is required for predicting more accurate ULS results.
Recently, Mei and Wang [31] also proposed single line shape empirical formulation in predicting ULS of stiffened panel, which is a similar shape that proposed by Lin [25] and Paik-Thayamballi [26]. However, they limited the maximum order of plate and column slenderness ratios as second and third, respectively. The obtained FE results were limited to propose empirical formulation, so the outcome was not compared in this study.
Lastly, the common structural rule (CSR) method proposed by International Association of Classification Societies (IACS) [32] may also be considered as industrial standard for practical purposes. In April 2006, CSR for oil tankers [33] and bulk carriers [34] was proposed. With regards to design rule and methodology, detailed comparison with historical summaries was documented by Jiameng et al. [35] and Kim et al. [21]. Therefore, it will not be discussed further here.
3. Ultimate Strength Calculations by the Non-Linear Finite Element Method (NLFEM)
Limit state design (LSD), also known as load and resistance factor design (LRFD), is now a well-known design method in the field of structural engineering. The LSD [16] includes ultimate limit state (ULS), fatigue limit state (FLS), accidental limit state (ALS), and serviceability limit state (SLS). Among others, a number of studies have been conducted on ULS-based design, technique, and its applications in terms of ULS application to stiffened panels [36,37,38,39], scaling effect [40,41], stiffened panels with opening [42], dynamic ULS [43], low temperature effect [14], and prediction of ULS by an artificial neural network (ANN) [44]. The FLS of an offshore riser by an ANN and simplified method [45,46], the ALS and ULS of a non-ice class aged ship [4,47], a steel plated structure [3], a ship’s hull [2], and Floating Liquefied Natural Gas (FLNG) [48,49] are also investigated.
3.1. Selection of Scenarios of Flat-Bar Stiffened Panel
Recently, Kim et al. [24] conducted a wide range of technical reviews on existing empirical formulation in predicting the ultimate strength of stiffened panel subject to longitudinal compression. In addition, they tested the accuracy of the existing empirical formulations by conducting 10,500 cases of T-bar stiffened panel numerical simulations using ANSYS NLFEM. The total number of stiffened panel scenarios were selected, as shown in Equation (7).
1a × 1b × 7tp × 10hw × 6tw × 5bf × 5tf = 10,500 scenarios
whereais the plate length,bis the plate breadth,tpis the plate thickness,hwis the web height,twis the web thickness, is the flange breadth, andtfis the flange thickness.
In this study, the flat-bar stiffened panel, which consisted of plate and web, was targeted, which meant that flat-bar stiffened panel did not have flange. In this regard, the 10,500 scenarios, as mentioned above in Equation (7), could be reduced to 420 cases by neglecting the flange part (bfandtf), so the scenarios of flange breadth (bf) and flange thickness (tf ) could be not considered in this study. More details of the 420 scenarios, including material and geometric properties, can found in Equation (8) and Table 2.
1a × 1b × 7tp × 10hw × 6tw = 420 scenarios
3.2. Structural Modelling
It is recognised that ultimate strength behaviour of stiffened and unstiffened panels varies depending on structural modelling technique, which means that the application of the appropriate FE structural modelling technique is essentially required in order to obtain realistic outcomes. In particular, the effect of assumed boundary condition, model size in longitudinal and transverse directions, material modelling technique, mesh size, initial imperfections, such as initial deflection and welding-induced residual stress, and many other elements should be carefully taken into consideration [50].
International Ship and Offshore Structures Congress (ISSC) [7] conducted a wide range of parametric studies on ultimate strength of stiffened panel by considering the changes of geometries. They studied the effect of model size on ULS by selecting one bay/one span and two bay/two span (1/2-1-1/2 model) stiffened panel models, as shown in Figure 3a,b. The details of boundary conditions for both models shown in Figure 3a,b are summarised in Table 3a,b, respectively. As expected, the one bay/one span model resulted in the overestimation of the ULS value, rather than the two bay/two span model. This was due to the effect of sideways deformation of the stiffeners located at the transverse frames, which were not allowed for the one bay/one span model. This was also caused by the effect of boundary condition assumed in both models. The imperfection sensitivity and geometric effects in stiffened panel was also studied by Ahmer Wadee and Farsi [51].
Based on the findings by ISSC [7], the two bay/two span model, as shown in Figure 3b and Table 3b was adopted in this study with the average level of initial deflection for plate and initial distortions for stiffeners. Details on initial deflection can be found in the Appendix. The welding-induced residual stress effect was not considered in this study. It was reported that 10–13% decrement of ULS of stiffened panel is expected to be achieved due to the effect of welding-induced residual stress [52]. The number of mesh in plate and web part is 10 and 6, respectively, based on mesh convergence studies [7,21].
With regards to structural modelling technique, some recent studies can be referred to in terms of scaling effect and [40,41] geometric evaluation [53]. In the case of element type, the Shell181 element with four nodes and six degrees of freedom per node were adopted. In addition, the classical metal plasticity model in the ANSYS material library was used in this study. The applied material model was based on von Mises yield criterion to define isotropic yielding and isotropic hardening behaviour. In the present study, the bilinear type material curve with no tangential slope (no hardening effect) was adopted. For example, the yield stress will not change as plastic strain increased.
ISSC [54] also reported the effect of hardening on the ultimate strength of structures. The hardening behaviour appears only after material yields, when the stiffened panel undergoes large deformations and plastic strains appear. However, a bilinear curve with no hardening may help to secure the additional safety margin. Based on this, ISSC [7] also adopted a bilinear shape material curve, which consists of the combination of material elastic modulus and yield strength.
For the FE simulation, the arc-length method was employed in conjunction with the modified Newton Raphson method in standard and modified forms, which enabled us to obtain an accurate load-shortening curve. In addition, all simulations were conducted based on the above method, with a large displacement option in the ANSYS NLFEM. The above FE simulation setting may be capable and was recommended for implementing buckling analysis of the complex structural shape. 3.3. Structural Analysis and Results
Prior to conducting numerical simulation, the results of FE technique validation with ISSC [7] is shown in Figure 4a. A total of 420 numerical simulations by ANSYS were conducted to obtain ULS of flat-bar stiffened panel under longitudinal compression. As shown in Figure 4b, the ULS trends were plotted based onλandβ. As expected, ULS tended to decrease when theλincreased. When the plate slenderness ratio increased or the plate was getting thinner, the ULS tended to maintain a general trend. This is t discussed further in the following section, in comparison with empirical formulation.
4. Development of Empirical Formulation and Verification of its Applicability 4.1. Empirical Formulation for Flat-Bar Stiffened Panel
Recently, an advanced empirical formulation was conducted by Kim et al. [6] to predict the ultimate strength of the T-bar type stiffened panel under longitudinal compression. A wide range of the numerical simulations of 10,500 cases in total were conducted by considering the changes of the geometric properties in terms of plate and stiffener.
This confirmed that the fluctuation behaviour of the ULS was found in the small range of the column slenderness ratio, as illustrated in Figure 1. In order to predict the ULS of stiffened panel more accurately, Kim et al. [6] additionally considered two important parameters, i.e., web slenderness ratio,hw/tw, and moment of inertia of stiffener to moment of inertia of plate ratio in z-direction,Ipz/Isz=(hw tw3 + tf bf3)/(tp b3), in addition to the two original basic parameters, i.e.,βandλ, as shown in Equation (1).
The numerical simulation results obtained were analysed by a data processing technique, and the polynomial function shape empirical formulation [6] is presented in Equation (9). By adopting the proposed empirical formulation, 15 coefficients for flat-bar were newly obtained based on FE numerical simulation in this study. Table 4 shows the 15 coefficients for T-bar and flat-bar, consisting of the empirical formulations. The data processing procedure proposed by Kim et al. [6] is shown in Figure 5.
σxuσYeq=[c0+(c1+c2λ+c3β+c4hwtw+c5IpzIsz)λ+(c6+c7β+c8hwtw+c9IpzIsz)1β +(c10+c11hwtw+c12IpzIsz)hwtw +(c13+c14IpzIsz)IpzIsz]≤1.0
The ULS results obtained by ANSYS FE simulation in Figure 4b were directly compared with results obtained by the empirical formulas in Figure 6a–h. In general, the empirical formulation shows relatively good agreement with ANSYS results, based on R2 results. This can also be confirmed by mean and coefficient of variation (COV) values in Table 5 and Figure 7a. When theβ increased, the ULS tended to have more of a general shape heading to the southeast direction, as shown in Figure 6e–h. In this study, the general shape represented a tendency to regularly decrease diagonally. In particular, this tendency could be observed when the plate was thin. It seems that the buckling of the plate element was affecting the overall collapse behaviour of the stiffened panel. As shown in Figure 6b–d, when the plate slenderness ratio was less than 1.8, which is generally considered a thick plate, ULS tended to fluctuate greatly, based on the variation of the stiffener size. A good understanding of accuracy was required for eachβwhen using the empirical formulation. In particular, this formulation may need a thin plate, having a large value ofβ , as shown in Figure 6e–h.
In this study, we verified that the applicability of the empirical formulation in Equation (9) can be extended to flat-bar stiffened panel with a new set of coefficients as summarised in Table 4. This shows relatively good agreement with ANSYS results, with the maximum and minimum range of the R2 values (0.8881 ≤ R2 ≤ 0.9435). However, it can be further improved by studying the phenomenon of the flat-bar stiffened panel under longitudinal compression in future. Particularly, the collapse behaviour of the plate under longitudinal compression should be studied.
4.2. Statistical Analysis for Verification of Developed Empirical Formulation
The ULS results obtained by various methods, such as numerical simulations [55,56] and empirical formulations [21,25,26,28,29], together with the proposed refined empirical formulation in this study, as well as ALPS/ULSAP [56], which is a computer program for ultimate limit state assessment for stiffened panel (ULSAP), were plotted in Appendix B and presented in Figure A1a–g. Details on ALPS/ULSAP can be found in Appendix C. In the case of Zhang and Khan [28], they limited it within the range ofλ≤2 . The detailed comparisons were conducted by statistical analysis, as summarised in Table 5. The statistical analysis results are also represented in Figure 7a–h. As expected, design formulations, such as J-O, P-R, and Euler tended to overestimate the ULS by about 55–65%, compared to the ANSYS FE numerical simulation, referring to the mean value in the “ALL” column shown in Table 5 (1.5463 ≤ mean ≤ 1.6539 and 0.1922 ≤ COV ≤ 0.1932).
In the case of the empirical formulations, improved results of mean and COV were observed, compared to design formulations (1.1225 ≤ mean ≤ 1.3922 and 0.1395 ≤ COV ≤ 0.1662). Most of the existing empirical formulations slightly overestimated the ULS values. On the other hand, ALPS/ULSAP, which is considered am analytical solution, underestimated ULS values by about 17–18%, compared to ANSYS numerical simulation results (mean = 0.8260 and COV = 0.4046). In particular, severe underestimation was observed when the column slenderness ratio (λ) between 0.1 and 0.3. If this range ofλwas excluded, this study showed that the mean and COV values were significantly improved to 0.9912 and 0.1389, respectively. In this study, we selected a reliable but limited range of flat-bar type stiffened panel. The empirical formulation proposed by this study provided well-fitted ULS results with ANSYS FE simulations (mean = 1.0024 and COV = 0.0583).
In summary, design formulations, which are generally adopted in shipbuilding, overestimate ULS values rather than empirical formulations and analytical solutions (ALPS/ULSAP). Most of the existing empirical formulations show good agreement with the refined FEM results by ANSYS. However, single line shaped empirical formulations still have limitations for predicting the ULS values accurately. The ALPS/ULSAP, which is an analytical method-based solution, can be considered a reliable way to take into account the nonlinearity of the ULS values. However, ALPS/ULSAP is recommended only when theλis greater than 0.3. An additional advantage of ALPS/ULSAP is that it enables robust design through pessimistic analysis results.
Lastly, the proposed empirical formulation in this study considered two more parameters mentioned in Section 4.1, so that it enabled us to predict ULS values and its nonlinearities accurately.
4.3. Additional Considerations
With regards to the effect of material grade on ultimate strength of stiffened panel, Li et al. [57] provided the investigation results using the ABAQUS numerical simulation code by using the arc-length method and allowing large displacement of the structures. They only modelled the single plate-stiffener combination (PSC) model, but the results may be useful for the comparison of different yield strength effects. They adopted the hungry horse mode, based on information by ISSC [58], and considered four steel grades, such as 235, 315, 355, and 390 MPa, of yield strength. They concluded that increasing the yield strength can effectively improve the ULS of stiffened panels, except for those slender members with large column slenderness ratio. They also mentioned that ULS is dispersed between 0.2 and 0.8 when column slenderness ratio ranges between 0.4 and 0.8.
Prior to this study, Kim et al. [21] stated this fluctuation behaviour in small column slenderness ratios. Li et al. [57] also mentioned that the maximum reduction due to welding-induced residual stress effect might be 10%. It was also stated by Khan and Zhang [52] that 10–13% of ultimate strength decreased due to the effect of compressive residual stress. The detailed investigation on the effect of geometrical imperfections can be found in the research by ISSC [7,54]. With regards to initial imperfection, ISSC [59] conclude that less than 5% of the ultimate strength may be influenced by ultimate strength.
5. Conclusions
In this study, the refined empirical formulation was proposed to predict ultimate strength performance or ultimate limit state (ULS) of flat-bar type steel stiffened panel under longitudinal compression. In total, 420 cases of numerical simulations by the ANSYS non-linear finite element method (NLFEM) were conducted and used as input data to develop an empirical formulation. The formulation obtained showed good agreement with ANSYS results in general (0.8881 ≤ R2 ≤ 0.9485), as shown in Table 6. In conclusion, it was verified that the obtained empirical formulation obtained was well fitted with ANSYS numerical simulation results (R2 = 0.9435).
The detailed results are summarised as follows. 5.1. Findings
-
When the plate slenderness ratio (β ) increased, the ULS tended to be a generalised shape, which represented a tendency to regularly decrease diagonally, as shown in Figure 6e–h. This may have been caused by the buckling of the plate element, which affected the overall collapse behaviour of the stiffened panel. In addition, this trend was observed when the plate was considered thin (1.8≤β)
-
As represented in Figure 6b–d, when the plate got thick, the ultimate limit state (ULS or ultimate strength) tended to fluctuate greatly depending on the variation of the stiffener size.
-
Two parameters, i.e., column slenderness ratio (λ) and plate slenderness ratio (β ), were considered as the main parameters of the existing empirical formulations. As indicated by Kim et al. [21,24], single line-shaped existing empirical formulations may not be able to implement the fluctuation behaviour of ULS. This means that the additional parameters should be considered in predicting the accurate ULS in the region of fluctuation.
-
The applicability of the refined empirical formulation proposed by Kim et al. [6] was tested by statistical analysis. It was confirmed that proposed empirical formulation can be applied to flat-type stiffened panel with 14 modified coefficients.
-
The empirical formulation proposed by this study was based on the ANSYS numerical simulation results with assumed scenarios in Table 2 and boundary conditions in Table 3b. It is well recognised that data processing depends on the input data. This means that other types of input data, e.g., ULS values obtained by an experimental or analytical method, would provide slightly different final outcomes. Nevertheless, the numerical simulation results assumed by a simply supported boundary condition with average level initial deflection may help designers in the robust design of ships and offshore structures by maintaining the additional structural safety margin.
- This study only considers the prediction of ULS of steel stiffened panel, i.e., a ship’s deck or upper side shell stiffened panel under longitudinal compression. Other types of applied loadings, such as biaxial compression and lateral pressure, should also be taken into consideration in future.
- With regard to initial imperfection, initial deflection of plate and initial distortion of stiffener are only considered in this study, while welding-induced residual stress was not considered.
- In addition, other types of stiffener (angle-bar type) should also be studied to develop the empirical formulation.
Coefficient | Lin [26] | P–T* [26] |
---|---|---|
C1 | 0.960 | 0.995 |
C2 | 0.765 | 0.936 |
C3 | 0.176 | 0.170 |
C4 | 0.131 | 0.188 |
C5 | 1.046 | −0.067 |
* Paik and Thayamballi (P–T).
(a) | |||||
Material Properties | |||||
Yield Strength | Elastic Modulus | Poisson’s Ratio | |||
Plate and Stiffener | 315 MPa | 205.8 GPa | 0.3 | ||
Geometric Properties (Unit = mm) | |||||
Plate | Length | 4150 | |||
breadth | 830 | ||||
Thickness | 9.5, 11, 14, 16, 21.5, 32.5, and 44.5 | ||||
Stiffener | Web height | 200, 284, 300, 360, 425, 460, 500, 700, 800, and 1,000 | |||
Web thickness | 10, 11.5, 12.5, 13.5, 20, and 28 | ||||
Summary | |||||
Plate slenderness ratio | 3.4181, 2.9620, 2.3194, 2.0295, 1.5103, 0.9991, and 0.7297 | ||||
Column slenderness ratio (PSC * model) | 420 cases | ||||
(b) | |||||
Geometrical Variable | No. of Scenarios | Selected Scenarios | |||
Plate length (a) | 1 | 4150 mm | |||
Plate width (b) | 1 | 830 mm | |||
Plate thickness (tp) | 7 | 44.5, 32.5, 21.5, 16, 14, 11, 9.5 mm | |||
Web height (hw) | 10 | 200, 284, 300, 360, 425, 460, 500, 700, 800, 1000 mm | |||
Web thickness (tw) | 6 | 10, 11.5, 12.5, 13.5, 20, 28 mm | |||
Terms | Coefficients | |
---|---|---|
T-Bar [6] | Flat-Bar (Present) | |
C0 | −0.1449 | −1.5721 |
C1 | 2.9787 | 5.6591 |
C2 | −2.6098 | −3.7336 |
C3 | −0.2418 | −0.6934 |
C4 | 1.2374 × 10−3 | −1.8581 × 10−2 |
C5 | 1.3470 × 10−2 | 1.7858 × 10−2 |
C6 | 0.8841 | 1.3546 |
C7 | −0.3361 | −0.3482 |
C8 | 1.5975 × 10−3 | −1.9443 × 10−3 |
C9 | 2.7745 × 10−3 | 0.8850 × 10−3 |
C10 | −7.5919 × 10−3 | 1.8299 × 10−2 |
C11 | 3.2442 × 10−5 | −1.2316 × 10−4 |
C12 | 4.9670 × 10−5 | 1.4994 × 10−4 |
C13 | 1.3267 × 10−2 | −1.8752 × 10−4 |
C14 | −5.4149 × 10−5 | −1.6306 × 10−5 |
Note: The permissions of reproduction for this T-bar information [6] in Equation (9) and Table 4 are given by Elsevier.
Existing Methods | Plate Slenderness Ratio (β) | ||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
0.7297 | 0.9991 | 1.5103 | 2.0295 | 2.3194 | 2.9520 | 3.4181 | ALL | ||||||||||
Mean | COV | Mean | COV | Mean | COV | Mean | COV | Mean | COV | Mean | COV | Mean | COV | Mean | COV | ||
Design Formulations | J-O | 1.2958 | 0.1634 | 1.3206 | 0.1453 | 1.3994 | 0.1084 | 1.5453 | 0.1151 | 1.6289 | 0.1267 | 1.7750 | 0.1462 | 1.8591 | 0.1577 | 1.5463 | 0.1932 |
P-R | 1.3047 | 0.1777 | 1.3337 | 0.1618 | 1.4122 | 0.1179 | 1.5567 | 0.1136 | 1.6393 | 0.1224 | 1.7825 | 0.1398 | 1.8640 | 0.1511 | 1.5562 | 0.1924 | |
Euler | 1.4282 | 0.2393 | 1.4344 | 0.2071 | 1.4955 | 0.1327 | 1.6413 | 0.1037 | 1.7277 | 0.1067 | 1.8818 | 0.1212 | 1.9684 | 0.1316 | 1.6539 | 0.1922 | |
Empirical Formulations | Lin | 1.1532 | 0.1314 | 1.1473 | 0.1189 | 1.1349 | 0.1148 | 1.1436 | 0.1383 | 1.1392 | 0.1491 | 1.0939 | 0.1637 | 1.0456 | 0.1723 | 1.1225 | 0.1452 |
P-T | 1.2112 | 0.1528 | 1.1933 | 0.1432 | 1.1587 | 0.1161 | 1.1508 | 0.1270 | 1.1415 | 0.1404 | 1.0909 | 0.1621 | 1.0409 | 0.1743 | 1.1410 | 0.1530 | |
Z-K | 1.2987 | 0.1466 | 1.3149 | 0.1354 | 1.2455 | 0.1047 | 1.2691 | 0.1198 | 1.2899 | 0.1324 | 1.3159 | 0.1517 | 1.3240 | 0.1626 | 1.2938 | 0.1395 | |
Kim | 1.1448 | 0.1465 | 1.0990 | 0.1353 | 1.0730 | 0.1080 | 1.1186 | 0.1336 | 1.1487 | 0.1493 | 1.1936 | 0.1717 | 1.2150 | 0.1840 | 1.1418 | 0.1562 | |
Xu | 1.3531 | 0.1767 | 1.3212 | 0.1333 | 1.2891 | 0.1125 | 1.3440 | 0.1253 | 1.3871 | 0.1386 | 1.4839 | 0.1652 | 1.5668 | 0.1824 | 1.3922 | 0.1662 | |
Present | 0.9857 | 0.0904 | 1.0363 | 0.0712 | 0.9909 | 0.0565 | 0.9886 | 0.0306 | 0.9958 | 0.0299 | 1.0064 | 0.0364 | 1.0128 | 0.0474 | 1.0024 | 0.0583 | |
Analytical Solution | ALPS/ULSAP | 0.8704 | 0.3479 | 0.9012 | 0.3456 | 0.8965 | 0.3843 | 0.8527 | 0.3963 | 0.7940 | 0.4180 | 0.7342 | 0.4396 | 0.7331 | 0.4618 | 0.8260 | 0.4046 |
* Johnson and Ostenfeld (J-O); Perry and Robertson (P-R) Lin = Lin (Lin) [25]; Paik and Thayamballi (P-T) [26], Zhang and Khan (Z-K) [28], Kim et al. (Kim) [21], and Xu et al. (Xu) [29].
Plate Slenderness Ratio (β) | R2 |
---|---|
0.7297 | 0.9347 |
0.9991 | 0.9180 |
1.5103 | 0.8958 |
2.0295 | 0.9392 |
2.3194 | 0.9485 |
2.9620 | 0.9296 |
3.4181 | 0.8881 |
Total | 0.9435 |
Author Contributions
Conceptualization, D.K.K.; methodology, D.K.K., S.Y.Y.; software, H.L.L.; validation D.K.K. and H.L.L.; formal analysis, D.K.K., S.Y.Y., H.L.L., and N.-K.C.; investigation, N.-K.C.; writing-original draft preparation, D.K.K. and S.Y.Y.; writing-review and editing, N.-K.C.; All authors have read and agreed to the published version of the manuscript.
Funding
This study was supported by the Research Program funded by the SeoulTech (Seoul National University of Science and Technology.
Acknowledgments
Authors appreciate kind supports from POSTECH (Korea), UTP (Malaysia), and Newcastle University (UK).
Conflicts of Interest
The authors declare no conflict of interest.
Appendix A. Initial Imperfections
Initial imperfections including initial distortion (initial deflection of plate and initial distortion of stiffener) and welding induced residual stress should be carefully taken into account for the ultimate strength analysis of structures. In this study, we only considered initial distortions.
The three types of initial distortions considered in this study are summarised as follows. In general, it can be categorised as plate initial deflection (wopl), the column type initial distortion of the stiffener (woc), and the sideways initial distortion of the stiffener (wos ). In case of plate initial deflection, there are various modes such as hungry horse mode, mountain mode, spoon mode, sinusoidal mode, and buckling mode [60,61]. Among them, buckling mode is adopted in this study.
1. Buckling mode initial deflection of plate:
wopl=Aosin(mπxa)sin(πyb)for plate;
2. Column type distortion of stiffener:
woc=Bosin(πxa)sin(πyB)for stiffener (Column-type);
3. Sideways initial distortion of stiffener:
wos=Cozhwsin(πxa)for stiffener (Sideways).
where, Ao, Bo, and Co = coefficients of the initial distortion (0.1β2t, 0.0015a, and 0.0015 respectively), m = buckling mode of the plate which provides minimum integer satisfyinga/b≤m(m+1) , a = plate length, b = plate breadth, B = breadth of stiffened panel (1 bay-1 span), hw = web height.
Appendix B. Comparison of ULS Results with Existing Methods
Jmse 08 00605 g0a1a 550Jmse 08 00605 g0a1b 550
Figure A1.Comparison of ULS results.
Figure A1.Comparison of ULS results.
Jmse 08 00605 g0a1b
Appendix C. ALPS/ULSAP Program
The ALPS/ULSAP program, developed for assessing ultimate limit state of stiffened panel is calculating its ultimate strength based on assumed six (6) types of collapse modes as shown in Figure A2. Once the ultimate strength values for six different modes are calculated, the lowest ultimate strength value will then be selected as ULS of stiffened panel. In this regards, ALPS/ULSAP gives relatively lower value of ultimate strength than other methods but it provides additional safe margins. This ALPS/ULSAP results are also presented in Figure A1a-g.
In this study, stiffener-induced collapse by beam-column type collapse (collapse mode 3) and stiffener-induced collapse by web buckling (mode 4) were observed from the most of the scenarios while overall collapse mode (mode 1) is also detected in typical scenario.
Details may be found in the ISSC [7] report, as well as by following references [60,62,63]. In addition, ALPS/ULSAP program is also available using MAESTRO modeller together with ALPS/HULL program [63].
Jmse 08 00605 g0a2 550
Figure A2.ALPS/ULSAP method and assumed collapse mode for ULS calculation, ALPS = nonlinear Analysis of Large Plated Structures; ULSAP = Ultimate Limit State Assessment of stiffened Panel) [64].
Figure A2.ALPS/ULSAP method and assumed collapse mode for ULS calculation, ALPS = nonlinear Analysis of Large Plated Structures; ULSAP = Ultimate Limit State Assessment of stiffened Panel) [64].
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2020. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
This research aims to study the ultimate limit state (ULS) behaviour of stiffened panel under longitudinal compression by a non-linear finite element method (NLFEM). There are different types of stiffeners mainly being used in shipbuilding, i.e., T-bar, flat-bar, and angle-bar. However, this research focuses on the ultimate compressive strength behaviour of flat-bar stiffened panel. A total of 420 reliable scenarios of flat-bar stiffened panel were selected for numerical simulation by the ANSYS NLFEM. The ultimate strength behaviours obtained were used as data for the development of closed form shape empirical formulation. Recently, our group proposed an advanced empirical formulation for T-bar stiffened panel, and the applicability of the proposed formulation to flat-bar stiffened panel is confirmed by this study. The accuracy of the empirical formulation obtained for flat-bar stiffened panel was validated by finite element (FE) simulation results of statistical analysis (R2 = 0.9435). The outcome obtained will be useful for ship structural designers in predicting the ultimate strength performance of flat-bar type stiffened panel under longitudinal compression.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer