Next Article in Journal
The Influence of Pre-Chamber Parameters on the Performance of a Two-Stroke Marine Dual-Fuel Low-Speed Engine
Previous Article in Journal
Uncertainty Analysis and Maneuver Simulation of Standard Ship Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design of Inner Ribs with Unequal Stiffness for Deep-Sea Highly Pressure-Resistant Cylindrical Shells and Utilizing NSGA-2 for Lightweight Optimization

1
Hubei Key Laboratory of Modern Manufacturing Quality Engineering, School of Mechanical Engineering, Hubei University of Technology, Wuhan 430068, China
2
Deep Sea Engineering Division, Institute of Deep Sea Science and Engineering, Chinese Academy of Sciences, Sanya 572000, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(7), 1231; https://doi.org/10.3390/jmse12071231
Submission received: 28 June 2024 / Revised: 18 July 2024 / Accepted: 19 July 2024 / Published: 21 July 2024
(This article belongs to the Section Ocean Engineering)

Abstract

:
For conducting scientific research at depths in the ocean, deep-sea probes are essential pieces of equipment. The cylindrical shell is the most sensible and rational packaging structure for these detectors. New technical challenges for enhancing the pressure resistance and lightweight design of the pressure-resistant cylindrical shell arise from the need to ensure that the detector packaging structure can withstand the immense water pressure at tens of thousands of meters in the underwater environment, while simultaneously reducing the detector packaging structure’s self-weight. This article examines the detection system’s deep-sea pressure-resistant cylindrical shell. To address the issue of the pressure-resistant shell’s insufficient ability to counteract the overall instability caused by the inability to form unstable half-waves in the radial direction when the ring rib pressure-resistant shell experiences it, a design method for the ribs inside the unequal-stiffness pressure-resistant cylindrical shell is suggested. The shell’s instability pressure increases by 9.65 MPa following the stiffness adjustment. Simultaneously, in order to attain even more lightweight optimization, the optimal inner rib section was obtained by applying the orthogonal topology optimization method, which also reduced the weight by 106.8 g and effectively improved the compression stability of the high-pressure cylindrical shell structure. Based on this, key optimization variables were found by performing sensitivity analysis on the cylindrical shell structure’s parameters. Then, with lightweighting as the primary objective, the high-pressure-resistant cylindrical shell’s optimal structural parameters were found using a multi-objective optimization process using the second-generation fast non-dominated genetic algorithm (NSGA-2). This resulted in a weight reduction of 1.2492 kg, or 17.26% of the original pressure-resistant shell. This has led to the development of a lightweight, highly pressure-resistant method for packaging marine exploration equipment structures.

1. Introduction

The ocean is full of resources for human development. However, traditional ocean observation technologies like ship-based observations, ocean buoys, and ocean satellites are no longer able to meet the current needs of ocean exploration due to the urgent need for further in-depth development of underwater resources. For the development of marine resources, high-pressure marine exploration equipment will, therefore, be essential. Submersible detectors come in a variety of shapes and sizes, but they are all made up of three basic parts: a connecting structure, an internal framework platform structure, and a pressure-resistant structure [1,2,3,4,5,6,7]. Considering that it makes up between a quarter and a half of the submersible’s total weight, the pressure-resistant shell structure is one of its most important parts. The primary duty of this structure is to tolerate high water pressure from the outside. To achieve the best pressure resistance performance and floating weight ratio, it is crucial to research and optimize the materials and structural forms of pressure-resistant shells.
Reducing the weight or minimizing the total weight of pressure-resistant shells while keeping the necessary strength and stability for operation is the major goal of optimization.
MingYang et al. [8] used the theory of thin shells to build a mechanical model of spherical pressure-resistant shells, which they then used to combine the Penalty Function Method (PFM) with the Multi-Population Genetic Algorithm (MPGA) to optimize the thickness, strengthening ribs, width of the strengthening ribs, and intersection angle of spherical pressure-resistant shells. A. S. Bhanu Prasan et al. [9] also used the PSO algorithm to optimize the quality of an entire pressure-resistant shell under specified buckling pressure and frequency constraints. Muhammad Imran et al. [10] optimized the pressure-resistant shell of a spherical submersible using the genetic algorithm (GA) in ANSYS, and then performed nonlinear buckling analysis on the pressure-resistant shell using a modified RIKS program in ABAQUS. In order to reach the maximum operating depth, Elsayed Fathallah et al. [11] optimized the entire pressure-resistant shell by reducing the elliptical submersible’s buoyancy coefficient. The design variables that were utilized were the longitudinal beam size, elliptical radius, and thickness of the pressure-resistant shell, all of which were subject to the limitations of failure criteria and shell buckling strength. Muhammad Imran et al. [12] introduced composite materials into cylindrical pressure-resistant shells, aiming to minimize buoyancy coefficients and maximize buckling load coefficients as optimization research objectives, and executed them through a coupled multi-objective genetic algorithm (MOGA). In order to optimize design variables and obtain the optimal cylindrical pressure-resistant shell, Cheng Wang et al. [13] proposed a cylindrical pressure-resistant shell with non-uniform arch ribs. Based on the parameterization model of the response surface method, they established an approximate model for the output response of ship mass. Honglei Liu et al. [14] proposed an adaptive morphogenesis algorithm based on the growth mechanism of leaf veins. They then applied this method to the design of stiffened plate and shell structures, resulting in a distinct stiffener distribution pattern. Paweł Foryś [15] addressed the impact of geometric defects on the form of equilibrium paths, employed an enhanced particle swarm optimization method (MPSO) for numerical optimization, and utilized the finite element method (FEM) to model and solve cylindrical shell structures. Ghasemi et al. [16] investigated the impact of hydrostatic pressure on cylindrical shells and their preparation by combining a genetic algorithm with geometric, structural, stress, and buckling constraints. Additionally, they used penalty functions to get rid of weak models. In the end, they suggested the best model, taking into account the least amount of weight under pressure. An optimization technique based on a free form was developed by Masatoshi et al. [17] and is especially useful for shape optimization in the out-of-plane direction. The buckling coefficient is the optimization’s objective function. The optimization results show that the buckling coefficient can be greatly increased using this method. Wu et al. [18] used the Newton–Lapson iteration method to perform a nonlinear static finite element analysis of axially compressed reinforced cylindrical shell structures with the goal of determining the ultimate load of structural instability. Mahdi et al. [19] optimized the shell and rib parameters under weight and frequency constraints using a genetic algorithm. They also improved the ribs’ material composition and shape, and the outcomes demonstrated that the improved shell had superior vibration characteristics. In order to optimize underwater pressure-resistant shells, Akl et al. [20] used a max–min multi-objective optimization strategy. The optimization objectives included shell vibration, noise, and quality. Sofiyev AH et al. [21] introduced the stability of non-uniform nanocomposite cylindrical shells (INCCSs) under hydrostatic pressure in a hot environment. The effective material properties of non-uniform nanocomposite cylindrical shells were modeled based on extended mixing rules. Based on the effective performance of the materials, the basic relationship and stability equation of the thermal environment were derived. The analytical expressions for the hydrostatic buckling pressure of INCCSs under the framework of FSDT and classical shell theory (CST) were obtained through solutions based on the Galerkin program. Wang M et al. [22] introduced the design and optimization of a pressure cylindrical shell composed of carbon/epoxy resin and glass/epoxy resin composite materials. A method for analyzing composite-pressure cylindrical shells was proposed, and the strength and buckling of composite-pressure cylindrical shells with different rib numbers were analyzed. The effects of thickness and layer angle on critical buckling and strength failure stress were studied. Yang Z et al. [23] used experimental methods to study the buckling pressure and failure strength of non-reinforced fiber-reinforced composite cylindrical shells under external hydrostatic pressure. They collected and discussed strain and pressure data. The results indicate that the improvement in buckling pressure by the ring reinforcement structure is more significant than the increase in failure pressure. The experimental burst pressure of the annular rigid cylindrical shell increased by 23.2%. Zhang X et al. [24] proposed a pressure shell with a corrugated structure to significantly improve its compressive performance. Seven pressure shells were prepared, including six inner corrugated pressure shells and one cylindrical pressure shell, and hydrostatic pressure tests were conducted. The results indicate that the corrugated structure can significantly increase the critical buckling load of the shell. Zhang Y et al. [25] proposed a novel bidirectional corrugated sandwich structure to improve the load-bearing capacity of cylindrical shells. The static and buckling analyses of sandwich shells and unreinforced cylindrical shells with the same volume-to-weight ratio were studied through numerical simulation. The results indicate that the proposed sandwich shell can effectively reduce the ratio of circumferential stress to axial stress from 2 to 1.25, increasing the critical buckling load by approximately 1.63 times. Numerical simulations show that optimizing and adjusting structural parameters can significantly improve the advantages of sandwich shells.
The ultra-high deep-water pressure environment that occurs tens of thousands of meters deep is typically ignored in current research. The standard uniform ribbed method for cylindrical shells may not be able to form unstable half-waves in the radial direction when encountering overall instability in the ultra-deep-water static pressure environment. This could lead to an insufficient ability of the pressure-resistant shell to resist instability. We suggest a technique to modify the internal rib stiffness in order to enhance the pressure-resistant shell’s overall anti-instability capabilities in order to address this general instability scenario. However, in terms of lightweight qualities, conventional, ordinary pressure-resistant materials like aluminum alloy 6061T6, resin-based carbon fiber T700, silicon nitride ceramics, etc. cannot be substituted for the titanium alloy TB9 material used to meet the pressure resistance conditions while minimizing volume at depths of tens of thousands of meters. An orthogonal topology optimization method is proposed to further reduce costs by optimizing the inner rib section after adjustment, given the high density and high cost of titanium alloy TB9 material. By combining the second-generation genetic algorithm, NSGA-2, for additional lightweight design in order to obtain the ideal pressure-resistant shell parameters, a multi-objective optimization design of the pressure-resistant shell is also proposed.
In the context of deep-sea high-pressure cylindrical shell structures, the standard uniform ribbed method for cylindrical shells may encounter issues with unstable half-waves in the axial direction when overall instability arises. Additionally, the method is unable to form unstable half-waves in the radial direction, which leaves the pressure-resistant shell with insufficient ability to withstand instability. The idea behind this article is to improve the overall instability resistance of the pressure-resistant shell by adjusting the uneven stiffness of the ribs inside using a method based on the Liz method. This paper suggests using the second-generation fast non-dominated genetic algorithm (NSGA-2) for the multi-objective optimization design of the pressure shell structure in order to obtain the optimal design parameters for the pressure shell. This will allow for the optimization of the inner rib section while maintaining good anti-instability ability and lead to a lightweight design of the entire device.

2. Pressure-Resistant Cylindrical Shells with Non-Uniform Internal Ribs Reinforced Structurally

2.1. Design of Pressure-Resistant Shell Structure Scheme with Inner Rib Reinforcement

Material and Basic Shape Selection for Reinforced Pressure-Resistant Shell with Internal Ribs

The marine environment is harsher and more complex than the terrestrial environment, and it also has many more unknowns for humans. Consequently, given that the pressure-resistant shell in this article operates in an environment that is 10,000 m below the surface, it is part of a strong, high-pressure working environment. The material used to make the pressure-resistant shell should possess qualities like corrosion resistance, light density, and high strength. The pressure-resistant housing’s primary material, titanium alloy TB9, was selected after pertinent processing technology and operating conditions were taken into account. The parameters of commonly used pressure-resistant materials are shown in Table 1 below.
Current pressure-resistant shells primarily take the shape of regular shapes like ellipsoids, spheres, conics, and cylindrical shapes. Apart from the conventional forms, there exist novel structural variations of pressure-resistant shells, including the ellipsoidal cross-section, lotus root, and corrugated tube.
A spherical pressure-resistant shell’s minimal volume-to-weight ratio and membrane stress are its key benefits, but there are also significant drawbacks. Due to the substantial difference between the spherical pressure-resistant shell’s shape and its streamline, there is a high fluid resistance it encounters when operating underwater, which leads to higher energy consumption. Additionally, the space utilization inside the shell is low and the effective load is small.
With its high space utilization and easy processing, the cylindrical pressure-resistant shell is currently the most popular pressure-resistant shell structure in use worldwide. Based on the dimensions of the sensors needed for this article (as indicated in Figure 1), a cylindrical structure is chosen because the total length of the sensor is approximately 330 mm and it has a long, square shape.

2.2. Stiffness Distribution Design for the Inner Ribs Using Liz Method

2.2.1. A Summary of Inner Rib Rigidity Adjustment Techniques

Currently, pressure-resistant shell reinforced rib sections are primarily designed in rectangular, L-shaped, T-shaped, and I-shaped configurations. The mass of the reinforced rib in the rectangular section is lowest when the section’s moment of inertia is equal. Initially, the shape of the rib section is determined by considering the processing complexity of various types of reinforced ribs. Initially, a rectangle was chosen.
In the study of pressure-resistant shell stiffening technology, the stiffness of each rib is usually considered in an identical manner; this is achieved by selecting the same material and cross-sectional shape for each inner rib. However, there is some logical error in the inner rib positioning. The main reason for this is that the pressure-resistant shell with ring ribs can only form an unstable half-wave in the axial direction when it experiences overall instability; it cannot form an unstable half-wave in the radial direction. This means that the pressure-resistant shell cannot withstand instability very well.
This paper employs the internal rib stiffness allocation technique to address this problem. This technique distributes the stiffness of each rib in accordance with relevant allocation principles, thereby avoiding the instability previously mentioned.
This article’s method for adjusting the stiffness of internal ribs involves modifying the size of the rib section, redistributing the moment of inertia of each rib section, and readjusting the stiffness of every internal rib, all the while maintaining the same number of ribs, rib moments, and total moments of inertia of all the rib sections.

2.2.2. Based on the Internal Rib Stiffness Allocation Using the Ritz Method, a Critical Instability Pressure Model for Pressure-Resistant Shells

The critical pressure formula for pressure-resistant shells with inner rib stiffness allocation is derived by applying the Liz method in conjunction with the conventional critical pressure formula for inner ribs with equal stiffness as a reference.
1.
The stress of equal-stiffness stiffened pressure-resistant shells can be calculated using this formula.
The following is the main formula used to determine the relevant stress when applying ring ribs to a pressure-resistant shell.
At the midpoint of two ribs and its rib edge, the shell’s normal stress [26] is
σ 1 = ( σ 1 ) x = 0 = P R t ( 0.5 ± ε 3 ) σ 1 = ( σ 1 ) x = 0.5 l = P R t ( 0.56 ε 2 )
where
(a)
σ 1 —at the midpoint of two ribs, the shell plate’s axial normal stress;
(b)
σ 1 —the normal stress axially at the rib’s edge;
(c)
P —the pressure that the pressure-resistant shell can withstand from the outside.
σ 2 = P R t ( 1 ε 4 ± μ ε 3 ) σ 2 = P R t ( 1 0.5 μ ) ( 1 ε 1 ) + μ σ 1 σ 2 0 = P R t ( 1 ε 4 )
where
(a)
σ 2 —the normal stress around the circumference of two ribs at their midpoint;
(b)
σ 2 —the normal stress circumferential at the rib’s edge;
(c)
σ 2 0 —the circumferential normal stress at the midpoint of two ribs on the shell plate.
σ f = E w 1 R = σ 2 μ σ 1 = P R t ( 1 0.5 μ ) ( 1 ε 1 )
where σ f —the axial normal stress of ribs.
In accordance with the “Classification and Construction Specification for Submarine Systems and Submarines” [27] requirements, the circumferential normal stress σ 2 0 on the middle surface of the shell plate at the midpoint of the two ribs, the axial normal stress σ f representing the ribs, and the axial normal stress σ 1 at the edge of the ribs must all be verified. The following is the verification formula [26]:
σ 2 0 = P R t ( 1 ε 4 ) 0.85 σ s σ 1 = P R t ( 0.56 ε 2 ) 1.15 σ s σ f = P R t ( 1 ε 1 ) ( 1 0.5 μ ) 0.6 σ s
2.
The Liz method determines the unstable pressure of a pressure-resistant shell and provides instructions for modifying the internal rib stiffness.
The Liz method’s pertinent steps are generally as follows [28]:
(1)
For the internal rib stiffness allocation to the pressure-resistant shell, derive the expression of the instability displacement correlation function. This function must be able to satisfy the boundary conditions and determine the displacement of the ribbed shell deviating from its initial position;
(2)
Find the expression for the strain and its displacement in different directions at any point in the shell;
(3)
Determine how the external force work and the shell’s strain energy are related to each other. Then, utilize this expression to determine the expressions for the shell’s total potential energy after deviation, its strain energy, and the external force work;
(4)
Determine the shell’s stability equation using the principle of minimum potential energy, and then solve for the minimum critical pressure.
3.
We calculate the theoretical critical pressure for pressure-resistant shells that have an internal distribution of rib stiffness.
As illustrated in Figure 2, establish a cylindrical coordinate system for the pressure-resistant shell. The graph’s angle indicates the angle formed by the research point’s diameter plane and the reference diameter plane, while the x-axis direction indicates the pressure-resistant shell’s axial direction.
The deformation at the rigid support and both ends of the shell can be considered continuous when the shell experiences unstable deformation. This continuous deformation can be broken down into axial, tangential, and radial directions [26].
u = A sin n φ cos m π x L v = B sin n φ cos m π x L w = C sin n φ cos m π x L
where
(a)
L —the length of the shell;
(b)
m —the half-wave number of axial instability when the shell is unstable;
(c)
n —the whole-wave number of circumferential instability when the shell is unstable;
(d)
u —the axial displacement of the shell;
(e)
v —the circumferential displacement of the shell;
(f)
w —the radial displacement of the shell; A, B, and C are constants.
With internal rib stiffness allocation, the strain energy of the shell plate for the pressure-resistant shell can be derived as follows, per “Submarine Strength” [27]:
V 1 = π L 4 R D R 2 2 ( 1 μ ) m 2 α 2 + 4 ( 1 μ ) m 2 α 2 n B C + ( m 2 α 2 + n 2 1 ) 2 + 2 ( 1 μ ) m 2 α 2 C 2 V 2 = π L 4 R E t 1 u 2 m 2 α 2 + 1 2 1 μ n 2 A 2 + 1 + μ m n α A B + n 2 + 1 2 ( 1 μ ) m 2 α 2 B 2 + 2 n B C + 2 μ m α A C + C 2 V = V 1 + V 2
where
(a)
V 1 —shell bending strain energy;
(b)
V 2 —strain energy in the middle surface of the shell,
α = π R / L .
Work of the external force U [27]:
U 1 = T 1 2 0 L 0 2 π ( w x ) 2 R d φ d x = π L 4 R T 1 α 2 m 2 C 2 U 2 = T 2 2 0 2 π 0 L x 2 w R d φ d x = π L 4 R T 2 ( n 2 1 ) C 2 U = U 1 + U 2
In order to determine the internal rib stiffness allocation of pressure-resistant shells, assume that the pressure-resistant shell has a total of i ribs [26]. Then, the strain energy is V 3 i , and the stiffness corresponding to the ribs at x i is I i .
V 3 i = E I i 2 R 3 π ( n 2 1 ) 2 C 2 sin 2 m π x i L
It is possible to calculate the strain energy V 3 of the entire pressure-resistant shell rib by adding up the strain energies of each rib [26]:
V 3 = i = 1 L / l 1 V 3 i = π L 4 R 2 E R 2 L ( n 2 1 ) C 2 i = 1 L / l 1 ( I i sin 2 m π x i L )
where L —the total length of the pressure-resistant shell;
l —rib spacing.
The formula for calculating the total potential energy of the pressure-resistant shell with internal rib stiffness adjustment can be obtained by combining Formulas (7)–(9).
= V 1 + V 2 + V 3 U 1 U 2 = π L 4 R D R 2 2 ( 1 μ ) m 2 a 2 B 2 + 4 ( 1 μ ) m 2 a 2 n B C + ( m 2 a 2 + n 2 1 ) 2 + 2 ( 1 μ ) m 2 a 2 C 2 + π L 4 R E t 1 μ 2 m 2 a 2 + 1 2 ( 1 μ ) n 2 A 2 + ( 1 μ ) m n a A B + n 2 1 2 ( 1 μ ) m 2 a 2 B 2 + 2 n B C + 2 μ m a A C + C 2 + π L 4 R 2 E R 2 L ( n 2 1 ) C 2 i = 1 L / l 1 ( I i sin 2 m π x i L ) π L 4 R T 1 α 2 m 2 C 2 π L 4 R T 2 ( n 2 1 ) C 2
By removing the high-order trace from the pressure-resistant shell with the internal rib stiffness allocation calculation formula, the stability equation under a uniform external load of the pressure-resistant shell with internal rib stiffness allocation can be obtained by replacing it with the basic relationship formula of the Liz method [26]:
T 1 m 2 a 2 + T 2 ( n 2 1 ) 2 = D R 2 ( n 2 1 + m 2 a 2 ) + E t m 4 a 4 n 2 + m 2 a 2 + 2 E R 2 L ( n 2 1 ) 2 C 2 i = 1 L / l 1 I i sin 2 m π x i L
Under uniform pressure P, the transverse and longitudinal sections of the pressure-resistant shell experience initial membrane forces denoted by T1 and T2 in the equation above.
It can be roughly described as follows, generally:
T 1 = 0.5 P R T 2 = P R
Formulas (11) and (12) can be combined to obtain the following theoretical critical pressure calculation formula for the overall instability of pressure-resistant shells with internal rib stiffness allocation:
P E = 1 n 2 1 + 0.5 m 2 a 2 D R 3 ( n 2 1 + m 2 a 2 ) 2 + E t R m 4 a 4 n 2 + m 2 a 2 + 2 E R 3 L ( n 2 1 ) 2 i = 1 L / l 1 I i sin 2 m π x i L
where m —the number of half-waves formed upwards when the pressure-resistant shell buckles;
n —the number of half-waves formed by the circumferential instability of the pressure-resistant shell.

2.2.3. Examining the Critical Instability Pressure of Cylindrical Shells with Adjustable Internal Rib Stiffness

The pressure-resistant shell has six evenly spaced ribs and measures L = 350 mm in total length, D = 90 mm in inner diameter, and t = 5 mm in thickness. The rib width is 5 mm, rib spacing is 50 mm, and rib thickness is 3 mm. The chosen material is titanium alloy TB9, with associated characteristics including elastic modulus E = 100 GPa, yield limit δs = 975 MPa, Poisson’s ratio μ = 0.31, and an applied load of 10 MPa.
The main calculation process is as follows:
1. Calculate the total moment of inertia of the ribs.
According to the Chinese CCS submersible design specifications, the formula for calculating the moment of inertia I of the ribs with plates is as follows [26]:
I = I 0 + l t 3 12 + ( y 0 + t 2 ) 2 × l t F l t + F
where
(a)
I 0 —the moment of inertia of each rib section;
(b)
F —the rib cross-sectional area;
(c)
y 0 —the distance from the axis of the rib section to the middle surface of the shell;
(d)
t —the thickness of pressure-resistant shell;
(e)
l —the rib spacing.
The rib section with plates’ moment of inertia can be computed as I = 7166.91 mm4 using Formula (14). Itotal = 43001.46 mm4 is the rib section’s total moment of inertia with six plates.
2. Determine the radial displacement ratio at which each unstable rib occurs.
When an instability occurs in a cylindrical shell, the shell plate’s radial displacement exhibits a trend of progressive decrease on both sides. The formula for radial displacement is as follows:
w = C sin ( n φ ) sin ( m π x L )
3. Reassign the stiffness of each rib.
Redistribute each rib section’s moment of inertia after utilizing Formula (15) to determine the radial displacement ratio between each rib. The following formula applies in this case:
I i = sin m π x i L Σ i = 1 sin m π x i L I t o t a l
In the formula: I i —the stiffness of each rib.
Each rib’s assigned stiffness can be computed after completing the aforementioned three steps. Each rib’s parameters are displayed in Table 2, where changing the rib thickness while keeping the rib width the same allows for the adjustment of each rib section’s moment of inertia.
Through the substitution of the rib section’s moment of inertia from Table 1 into Formula (13), the pressure-resistant shell’s overall instability critical pressure PE1 with inner rib stiffness allocation is 79.68 MPa, while the pressure-resistant shell’s overall instability critical pressure PE2 with equal stiffness is 70.24 MPa.

2.2.4. Analyzing Cylindrical Shells through Simulation with Adjustable Internal Rib Stiffness

The equal-stiffness internal rib shell pressure, PE2, and the overall instability critical pressure, PE1, were obtained in the previous section using the internal rib stiffness adjustment theory’s critical instability pressure calculation formula. ANSYS Workbench 2023 R1 software was used to perform the pertinent simulation calculations discussed in this section. The grid was set as a hexahedral grid with a grid size of 1 mm. The convergence process of the grid was automatically determined by the program, with fixed constraints at the two circular rings on the cylindrical surface. The pressure was set to 10 MPa on both sides of the cylinder. The accuracy of the theoretical formula was confirmed by comparing the obtained results with the theoretical values. At the moment, the most popular approach for researching structural stability is linear buckling analysis. Its primary feature, which is a finite element analysis, is to treat the entire structure as an elastic structure and ignore the original defects of the structure itself [29].
Figure 3 displays the characteristic values of the buckling deformation analysis obtained after the internal rib stiffness of the pressure-resistant shell was adjusted through ANSYS Workbench simulation calculations.
We obtained the Mises stress cloud maps that are pertinent, as illustrated in Figure 4 and Figure 5.
From the above results, it is evident that the maximum stress on the ribs is 96.753 MPa and the maximum stress on the equal-stiffness inner rib pressure-resistant shell is 131.08 MPa. The critical instability pressure in the first-order instability state is 72 MPa, which is represented by the buckling characteristic value of 7.2. The inner rib stiffness-adjusted pressure-resistant shell can withstand stresses up to 131.11 MPa, while the ribs themselves can withstand stresses up to 96.975 MPa. The first-order instability state has a buckling characteristic value of 8.165, which is equivalent to a critical instability pressure of 81.65 MPa. This confirms the accuracy of the theoretical formula for critical instability pressure derived in this article. It can be seen that the critical instability pressure in both simulation results is essentially consistent with the critical instability pressure calculated using the derivative theory formula in the previous section.
In terms of stability, the instability pressure of the inner rib stiffness adjustment method is 9.65 MPa higher than that of the inner rib with equal stiffness, despite the two stiffening methods having similar maximum stresses. Increasing the internal rib stiffness can boost the pressure-resistant shell’s compressive strength and stability, as can be observed.

2.3. Optimization of the Inner Rib Reinforcement Cross-Sectional Orthogonal Topology Using the Variable Density Approach

2.3.1. Design Strategy for Orthogonal Topology Optimization

With the maximum stress of the optimized strengthening rib as the primary optimization objective and the number of iterations needed for topology optimization as the second evaluation indicator, the goal is to obtain the cross-sectional shape of the strengthening rib with the strongest compressive strength. The variable density method is used to optimize the orthogonal topology of the rib’s cross-sectional shape. As a design variable, the range of relative element material density ρ in the optimization region is defined as [0, 1]. The pressure-resistant shell’s minimum weight serves as the objective function, and its maximum stress and maximum strain are set as constraint conditions. Formula (17) [30,31] illustrates how the relative element material density ρ represents the objective function in the mathematical model of the variable density method topology optimization.
f i n d ( ρ ) = ( ρ 1 , ρ 2 , ρ n ) min m ( ρ ) = V ρ d V s . t E ( ρ i ) = ρ i P E 0 V ρ d V M * ( i = 1 , 2 , n ) 0 ρ min < ρ i 1
where
(a)
ρ —the relative density of the unit materials, one of the design variables;
(b)
m ρ —the optimized design area’s quality, represented by the objective function;
(c)
V —the volume occupied by the design area;
(d)
E ρ i —the relative density of unit materials, represented by the relative elastic modulus;
(e)
P —the penalty factor;
(f)
E 0 —the material’s true elastic modulus;
(g)
M * —the optimized quality’s maximum value;
(h)
n —maximize the quantity of limited components within the design domain;
(i)
ρ m i n —the lowest relative density of the material in question.
Grid size t , penalty factor p , and volume factor f are the primary variables influencing the outcomes of topology optimization. We should prevent checkerboard phenomena, porosity, and non-convergence during the topology optimization process and accurately depict the structure’s mechanical characteristics. The selections in this article are 0.3 ≤ f ≤ 0.7, 3 ≤ p ≤ 5, and 0.5 mm ≤ t ≤ 1 mm.
The volume factor f , penalty factor p , and grid size t are the primary factors influencing the optimized object’s performance in terms of pressure resistance, as was observed in the previous step. The three parameters mentioned above are denoted as factors A, B, and C, respectively. It is presumed that there is no interplay among the three factors [32]. As indicated in Table 3, the level of factors for orthogonal topology optimization experiments is established based on the range of values for each factor suggested in the preceding text, as well as comparable research and design experience.
Based on this, the orthogonal experimental table needed for the experiment is shown in Table 4 (note that this article only contains three horizontal factors; therefore, a standard orthogonal table would require the addition of a set of blank factors).

2.3.2. Test Results and Analysis

Test Result

Table 5 displays the pertinent results of the topology optimization based on the variable density method that was carried out using ANSYS Workbench on the nine experiments that were previously mentioned. In topology optimization settings, the mesh is a hexahedral mesh, and the mesh size was automatically selected based on the parameters in Table 5. The mesh convergence process was automatically determined by the program, and the pressure was set to a uniform distribution of 10 MPa on the outer ring of the circular ring. At the same time, cylindrical support was applied to the outer ring, and free constraints were given in the radial direction.

Range Analysis of Experimental Results

The range analysis table for this experiment is displayed in Table 6. In the table, K i i = 1,2 , 3 denotes the average value of the test results corresponding to level i , and k i represents the sum of the test results corresponding to level i in any column. Making a relationship graph between various levels and target parameters to aid in comprehension of how different level factors affect target parameters is also practical, and is shown Figure 6.
The analysis presented above leads to the conclusion that the volume coefficient has the greatest influence on the maximum stress evaluation index, whereas the grid size has the least influence. The number of iteration steps is primarily determined by the penalty factor, with the volume coefficient having the least significant effect. The combinations of A3B1C3 and A1B2C3, respectively, represent the ideal solutions that correspond to the two evaluation indicators of maximum stress and iteration steps. With a maximum stress of 221.72 MPa and 16 iteration steps, A3B1C3 is present in the aforementioned experiment (as indicated by the pertinent simulation figure in Figure 7). Nevertheless, additional testing is required because A1B2C3 was not randomly identified in the orthogonal experiment. The pertinent findings are displayed in Figure 8. The maximum stress following A1B3C3 combination optimization is 285.59 MPa, and 12 iteration steps were used.
When the maximum stress of the A3B1C3 scheme is compared to the A1B2C3 scheme, it is 63.87 MPa less. The optimization process requires four more iterations for the A3B1C3 scheme than for the A1B2C3 scheme. But since achieving high pressure resistance is the primary goal of this article, it was ultimately determined that A3B1C3, with a volume factor of 0.6, a penalty factor of 3, and a grid size of 1 mm, is the best combination for cross-sectional topology optimization.

Analysis of Optimization Results

Subsequent manufacturing processes encounter certain challenges because of the optimized cross-section’s uneven and multiplanar edges. Post-processing procedures are therefore required to guarantee that the optimized model can be produced and processed with ease. Figure 9 displays the model’s processed rib section following imitation optimization. The rib section is referred to in this article as having a bow-like shape because of its shape.
Table 7 displays the information pertaining to the modeling and simulation analysis of the optimized bow section.
The stress dropped from 221.72 MP to 199.51 MP, a decrease of 22.21 MPa, even though the processed section’s mass increased by 5.02 g in comparison to the untreated reinforcing ribs. A weight reduction of 12.43 g (24.89%) and 106.8 g (the total weight of the six ribs) over the original rectangular section allowed for some partial fulfillment of the lightweight requirement.
The internal rib stiffness adjustment method described in the previous section was applied to the obtained section of the arched rib. Sequentially, we numbered the ribs from left to right. B1 is the formula for rib 1’s width, and so forth. Ribs 2, 3, 4, 5, and 6 have widths of b2, b3, b4, b5, and b6, in that order. The different parameters of the rib cross-section will be defined as follows, as shown in Figure 10, where h = b/10 and c = b/5, to guarantee that the cross-sectional shape of the ribs does not change.
A formula was used to determine each rib section’s moment of inertia, and the results are displayed in Table 8.
The three-dimensional schematic diagram is shown in Figure 11.
Finally, as indicated in Table 9, the initial parameters of the complete cylindrical pressure shell structure were obtained.

3. NSGA-2-Based Multi-Objective Parameterization Optimization of Cylindrical Pressure-Resistant Shells

3.1. Sensitivity Analysis of Parameters for Cylindrical Shells with Inner Rib Reinforcement under High Pressure

The primary design variables in this section are the thickness of the pressure-resistant shell (t), the rib thickness (b1, b2, b3, b4, b5, b6), and the rib spacing, L, which are based on the internal rib stiffness adjustment method suggested in the previous text. The rib thickness can be calculated using the internal rib stiffness adjustment method described in the preceding text, which is b1 = b6, b2 = b5, and b3 = b4. This article takes the thickness h of the end caps as an optimization parameter and includes the end caps at both ends of the pressure-resistant shell in the optimized model, in addition to the parameters mentioned above. We simultaneously set the Mises stress S, shell mass M, and buckling characteristic value Q as the output responses.
These are the initial values of the design variables: t = 5 mm, h = 15 mm, b1 = b6 = 9.1 mm, b2 = b5 = 9.91 mm, b3 = b4 = 10.78 mm, L = 52 mm. Table 10 illustrates the range of design variables pertaining to pressure-resistant shells.
For sensitivity analysis, we employed a global sensitivity analysis method based on variance. This paper used an optimal Latin sampling experiment to extract samples consistently and comprehensively. ANSYS Workbench was utilized to conduct simulation analysis after 60 sample points were uniformly extracted for every design variable. The experimental setup for sensitivity analysis was as follows: the mesh was set as a hexahedral mesh with a mesh size of 1 mm. The convergence process of the mesh was automatically determined by the program, and the fixed constraints were at the circular rings on both sides of the cylindrical surface. The pressure was set to 10 MPa on both sides of the cylindrical surface. Figure 12, Figure 13, Figure 14 and Figure 15 display the obtained sensitivity relationships for each design variable. The sensitivity of each design variable to the output is shown in Table 11.
The images and tables above make it easy to determine that, at 0.6983, 0.7105, and 0.7168, respectively, the shell thickness t has the highest sensitivity to stress, buckling, and mass. With a stress ratio of 0.0837, the end cap’s thickness h’s sensitivity to quality is a specific relationship. The most notable rib width is b3, which also explains a portion of the sensitivity in buckling and stress (0.1511 in buckling and 0.0623 in stress). With 0.0937 in stress and 0.0813 in buckling, b2 also has some effect in these areas.

3.2. Multi-Objective Parameterized Optimization Process

Generally speaking, multi-objective optimization is an optimization technique that has several optimization goals that, to varying degrees, directly impact the objective function. Determining the design variables, constraints, and objective function is the first step towards creating a mathematical model of the objective function [33].
Using the sensitivity analysis of each parameter from the previous text as a guide, this article will optimize each parameter by choosing multiple appropriate parameters to be used as design variables for multi-objective optimization.
Meanwhile, the pressure-resistant shell must also adhere to the pertinent clauses of the “Classification and Construction Specification for Submarine Systems and Submarines” and confirm the circumferential normal stress σ 2 0 on the middle surface of the shell plate at the midpoint of the two ribs, the axial normal stress σ f representing the ribs, the critical pressure of instability P c r , and the axial normal stress σ 1 at the edge of the ribs. Thus, the following are the constraint conditions:
y = f ( x 1 , x 2 , x s ) , X = ( x 1 , x 2 , x s ) C S
σ 2 0 = P R t ( 1 ε 4 ) 0.85 σ s σ 1 = P R t ( 0.56 ε 2 ) 1.15 σ s σ f = P R t ( 1 ε 1 ) ( 1 0.5 μ ) 0.6 σ s P c r 1.5 P
Here, we integrate the constraints, goal functions, and design variables that were suggested in the earlier text. For pressure-resistant shells, the following multi-objective optimization model can be created:
g k X h ρ X 0   or 0 , k , ρ 1 , 2 , 3 X = x 1 x 2 x 3 x i x n T n N x i m i n N x i m a x F X = f 1 X f 2 X f 3 X f i X f m X m i n F X = m i n f i X i = 1 , 2 , 3 , m
where
(a)
X —a set of design variables;
(b)
h ρ ( X ) —the function of boundary conditions;
(c)
g ρ ( X ) —the function of constraint conditions;
(d)
F ( X ) —the objective function for optimization.

3.3. Analysis of Optimization Results

Particle swarm optimization, genetic algorithms, ant colony algorithms, and other optimization algorithms are currently popular multi-objective optimization techniques. Apart from the commonly used algorithms previously mentioned, another enhanced version of genetic algorithms is the second-generation fast non-dominated genetic algorithm (NSGA-2). NSGA-2 offers several advantages over traditional optimization algorithms, including quick convergence, quick computation, and good global search performance. For every optimized data set, the second-generation fast non-dominated genetic algorithm will be chosen as the optimization algorithm in this article. The NSGA-2 flowchart is displayed in Figure 16.
A total of thirty sets of sample points were chosen for the experiment through uniform sampling within each design variable’s value range using the best Latin hypercube experiment technique. The 30 sets of sample points were then simulated and examined using ANSYS Workbench. Table 12 displays the pertinent sample points and output response values.
Currently, multi-objective optimization results are typically assessed using four indicators: relative absolute maximum error, relative average error, root mean square error, and coefficient of certainty R 2 . The degree to which variable X in the equation explains variable Y is indicated by the coefficient of certainty, or R 2 . The explanatory power increases with the distance between R 2 and 1. The prediction accuracy is represented by the relative absolute maximum error, relative mean error, and root mean square error. The prediction is more accurate the closer the error value is to 0. The coefficient of certainty R 2 and the relationship are as follows:
R 2 = i = 1 n i ( y i ^ y i ¯ ) 2 i = 1 n i ( y i y i ¯ ) 2
where N —the number of test sample points;
y i ^ —the predicted values of the model;
y i ¯ —the mean of the experimental response value;
y i —the response value of the model.
The comparisons between the simulation values and the expected response values for each output are displayed in Figure 17, Figure 18 and Figure 19.
A variety of indicators for evaluating the approximate model’s accuracy can be obtained using Workbench (see Table 13).
Table 13 shows that the model’s accuracy is good even though the stress values on the two evaluation indicators—relative maximum absolute error and relative average absolute error—are slightly higher but still fall within the permitted range of model error. Stress constraint conditions must be set when optimizing pressure-resistant shells multi-objectively with ANSYS Workbench. In compliance with “Submarine Strength” regulations, submarines delving deeper than two thousand meters are required to have a minimum safety factor greater than 1.5. With a safety factor of η = 1.5, the maximum stress limit established in this article is σ = σ S / η 650 MPa. The target is the maximum buckling characteristic value and the minimum mass. After obtaining the candidate design points (A~E), the pertinent parameters were determined, as displayed in Table 14. Simultaneously, the response values’ trend of changes was also acquired, as Figure 20 illustrates.
Table 15 displays the pertinent findings from the comparison of the five candidate points mentioned above with the original prototype. The quality difference between each candidate point and the prototype is indicated by the symbol M .
It is evident from Table 16. that every design point achieved a weight reduction of more than 1.2 kg in comparison to the original prototype, with candidate point C reaching the greatest weight reduction of 1.2492 kg. Regarding maximum stress, each candidate point’s maximum stress has grown considerably when compared to the original prototype, but it is still quite near to the predetermined 650 MPa, suggesting that all of the material has been used. Every potential point is above 4.3 in the characteristic buckling value, indicating that the stability satisfies the necessary conditions. This article will use candidate point C as the final design scheme for the pressure-resistant shell parameters, even though candidate points A through E meet the various performance requirements of the pressure-resistant shell. Candidate point C achieves the most weight reduction, reaching 1.2492 kg, accounting for 17.26% of the original pressure-resistant shell.
The final values of the various parameters of the pressure-resistant shell obtained are shown in Table 16. It is necessary to round off the various parameters of candidate point C due to the accuracy of the machining. Figure 21 displays the schematic diagram of the pressure-resistant shell structure.

4. Conclusions

In order to increase the critical instability pressure of pressure-resistant shells in settings at tens of thousands of meters of water depth, this article suggests combining the Liz method with an internal rib stiffness adjustment technique. By ensuring that the critical instability pressure met the working conditions, the shape of the inner ribs was optimized using the orthogonal topology optimization method in order to further achieve lightweight design. Based on this, the pressure-resistant shell’s multi-objective optimization design was ultimately optimized using a second-generation fast non-dominated genetic algorithm. The pressure-resistant shell was then further treated to make it lighter. The article’s primary findings are listed below:
(1)
In order to optimize deep-sea high-pressure structural layers, this article suggests a lightweight design approach for pressure-resistant devices in high-pressure deepwater environments that satisfies the non-unstable requirement.
(2)
The determined pressure-resistant shell was subjected to internal rib stiffness adjustment using the proposed internal rib stiffness adjustment method in conjunction with the Liz method. Through the use of ANSYS 2023 R1simulation software, this process was compared and verified with traditional uniform rib addition. There was a 9.65 MPa increase in the critical instability pressure following the modification of the internal rib stiffness and an improvement in the pressure-resistant shell’s overall anti-instability performance.
(3)
Titanium alloy TB9 material is unique in a high-pressure setting. An orthogonal topology optimization method was proposed for the inner rib section in order to reduce material costs. This method yielded the optimal inner rib arch section shape and resulted in a 106.8 g weight reduction. Based on this, a sensitivity analysis was performed on the pressure-resistant shell, and a second-generation fast non-dominated genetic algorithm was applied to carry out additional multi-objective optimization. After determining the pressure-resistant shell’s ideal structural parameters, 1.2492 kg of weight—or 17.26% of the original pressure-resistant shell—was reduced.
This article generally suggests a lightweight design approach for pressure-resistant devices in high-pressure water depth environments that satisfy the stability requirement, as well as an internal rib adjustment method to enhance their anti-instability ability. This offers a fresh method for optimization for lightweighting and strengthening the stability of deep-sea high-pressure-resistant structural layers, creating a process method that is both scientifically and engineeringly significant.

Author Contributions

Conceptualization, Y.H. and H.M.; methodology, X.W.; software, X.W.; validation, Z.L., Y.Y. and H.M.; formal analysis, X.W.; investigation, Z.L.; resources, H.M.; data curation, Y.Y.; writing—original draft preparation, X.W.; writing—review and editing, Y.H.; visualization, Z.L.; supervision, Y.Y.; project administration, H.M.; funding acquisition, Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key R&D Program Projects, grant number 2022YFC3104000.

Institutional Review Board Statement

The study did not require ethical approval.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available on request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ramadass, G.A.; Vedachalam, N.; Ramesh, S.; Sathianarayanan, D.; Subramanian, A.N.; Ramesh, R.; Chowdhury, T.; Pranesh, S.B.; Atmanand, M.A. Challenges in developing deep-water human occupied vehicles. Curr. Sci. 2020, 118, 1687–1693. [Google Scholar] [CrossRef]
  2. Ramadass, G.A.; Ramesh, S.; Vedachalam, N.; Subramanian, A.N.; Sathianarayanan, D.; Ramesh, R.; Harikrishnan, G.; Chowdhury, T.; Jyothi, V.B.N.; Pranesh, S.B.; et al. Unmanned underwater vehicles. Curr. Sci. 2020, 118, 1681–1686. [Google Scholar] [CrossRef]
  3. Joo, M.G.; Qu, Z. An autonomous underwater vehicle as an underwater glider and its depth control. Int. J. Control Autom. Syst. 2015, 13, 1212–1220. [Google Scholar] [CrossRef]
  4. Jaffre, F.; Littlefield, R.; Grund, M.; Purcell, M. Development of a New Version of the REMUS 6000 Autonomous Underwater Vehicle. In Proceedings of the OCEANS 2019—Marseille, Marseille, France, 17–20 June 2019. [Google Scholar]
  5. Purcell, M.; Gallo, D.; Packard, G.; Dennett, M.; Rothenbeck, M.; Sherrell, A.; Pascaud, S. Use of REMUS 6000 AUVs in the search for the Air France Flight 447. In Proceedings of the OCEANS’11 MTS/IEEE KONA, Waikoloa, HI, USA, 19–22 September 2011. [Google Scholar]
  6. Roper, D.; Harris, C.A.; Salavasidis, G.; Pebody, M.; Templeton, R.; Prampart, T.; Kingsland, M.; Morrison, R.; Furlong, M.; Phillips, A.B.; et al. Autosub Long Range 6000: A Multiple-Month Endurance AUV for Deep-Ocean Monitoring and Survey. IEEE J. Ocean. Eng. 2021, 46, 1179–1191. [Google Scholar] [CrossRef]
  7. McPhail, S. Autosub6000: A Deep Diving Long Range AUV. J. Bionic Eng. 2009, 6, 55–62. [Google Scholar] [CrossRef]
  8. Yang, S.; Wang, Y.; Liang, Y.; Wang, S.; Zhang, L.; Yang, M. Optimization design of neutrally buoyant hull for underwater gliders. Ocean. Eng. 2020, 209, 107512. [Google Scholar] [CrossRef]
  9. Prasanna, A.S.; Raju, K.K.S.; Ramji, K.; Satish, P. Free Vibration, Buckling and Design Optimisation of Composite Pressure Hulls. Mater. Today Proc. 2017, 4, 7381–7387. [Google Scholar] [CrossRef]
  10. Imran, M.; Shi, D.; Tong, L.; Waqas, H.M.; Muhammad, R.; Uddin, M.; Khan, A. Design Optimization and Non-Linear Buckling Analysis of Spherical Composite Submersible Pressure Hull. Materials 2020, 13, 2439. [Google Scholar] [CrossRef] [PubMed]
  11. Fathallah, E.; Qi, H.; Tong, L.; Helal, M. Design Optimization of Composite Elliptical Deep-Submersible Pressure Hull for Minimizing the Buoyancy Factor. Adv. Mech. Eng. 2014, 2014, 987903. [Google Scholar] [CrossRef]
  12. Shi, D.; Tong, L.; Elahi, A.; Waqas, H.M.; Uddin, M.; Imran, M. Multi-objective design optimization of composite submerged cylindrical pressure hull for minimum buoyancy factor and maximum buckling load capacity. Def. Technol. 2021, 17, 1190–1206. [Google Scholar]
  13. Wang, C.; Yang, M.; Wang, Y.; Ren, M.; Wang, Z.; Yang, S. Design and optimization of cylindrical hull with non-uniform arch ribs for underwater gliders based on approximate model and experiments. Ocean. Eng. 2022, 259, 111831. [Google Scholar] [CrossRef]
  14. Liu, H.; Li, B.; Yang, Z.; Hong, J. Topology optimization of stiffened plate/shell structures based on adaptive morphogenesis algorithm. J. Manuf. Syst. 2017, 43, 375–384. [Google Scholar] [CrossRef]
  15. Foryś, P. Optimization of cylindrical shells stiffened by rings under external pressure including their post-buckling behaviour. Thin-Walled Struct. 2015, 95, 231–243. [Google Scholar] [CrossRef]
  16. Ghasemi, A.R.; Hajmohammad, M.H. Minimun-Weight Design of Stiffened Shell under Hydrostatic Pressure by Genetic Algorithm; Techno-Press: Daejeon, Republic of Korea, 2015. [Google Scholar]
  17. Shimoda, M.; Okada, T.; Nagano, T.; Shi, J.-X. Free-form optimization method for buckling of hull structures under out-of-plane and in-plane shape variations. Struct. Multidiscip. Optim. 2016, 54, 275–288. [Google Scholar] [CrossRef]
  18. Wu, H.; Yan, Y.; Yan, W.; Liao, B. Adaptive Approximation-based Optimization of Composite Advanced Grid-stiffened Cylinder. Chin. J. Aeronaut. 2010, 23, 423–429. [Google Scholar]
  19. Mehrabani, M.M.; Jafari, A.A.; Azadi, M. Multidisciplinary optimization of a stiffened hull by genetic algorithm. J. Mech. Sci. Technol. 2012, 26, 517–530. [Google Scholar] [CrossRef]
  20. Akl, W.; Ruzzene, M.; Baz, A. Optimal design of underwater stiffened hulls. Struct. Multidiscip. Optim. 2002, 23, 297–310. [Google Scholar] [CrossRef]
  21. Sofiyev, A.H.; Fantuzzi, N. Stability analysis of shear deformable inhomogeneous nanocomposite cylindrical shells under hydrostatic pressure in thermal environment. Materials 2023, 16, 4887. [Google Scholar] [CrossRef]
  22. Wang, M.; Liu, F.; Zhang, J. Collaborative optimization of a composite pressure cylindrical shell based on dynamic penalty functions. Ocean. Eng. 2024, 309, 118542. [Google Scholar] [CrossRef]
  23. Yang, Z.; Zhang, X.; Pan, G.; Xu, Y. Buckling and strength failure for unstiffened and ring-stiffened composite shells. Ocean. Eng. 2023, 278, 114513. [Google Scholar] [CrossRef]
  24. Zhang, X.; Sun, J.; Wu, R.; Peng, Y. Buckling performance of inner corrugated pressure shells under external hydrostatic pressure. Ocean. Eng. 2023, 288, 115963. [Google Scholar] [CrossRef]
  25. Zhang, Y.; Chen, Y.; Yun, L.; Liang, X. Mechanical performance of bio-inspired bidirectional corrugated sandwich pressure shell under external hydrostatic pressure China. Ocean. Eng. 2024, 38, 297–312. [Google Scholar] [CrossRef]
  26. Wang, B.; Yang, M.; Zeng, D.; Hao, P.; Li, G.; Liu, Y.; Tian, K. Post-buckling behavior of stiffened cylindrical shell and experimental validation under non-uniform external pressure and axial compression. Thin-Walled Struct. 2021, 161, 107481. [Google Scholar] [CrossRef]
  27. Wei, R.; Shen, K.; Pan, G.; Li, Z. Experimental study on buckling behavior of composite cylindrical shells with and without circular holes under hydrostatic pressure. Compos. Struct. 2023, 322, 117434. [Google Scholar] [CrossRef]
  28. Cheng, Y. Structural Optimization Design of the Pressure Hull Based on Nonuniform Ring-Frames Theory. Master’s Thesis, Harbin Engineering University, Harbin, China, 2012. [Google Scholar]
  29. Feng, Y. HXJ180 Marine Workover Machine Derrick Lightweighting Study. Master’s Thesis, Yangtze University, Jingzhou, China, 2023. [Google Scholar]
  30. Du, Y.-X.; Wang, W.; Tian, Q.-H. Optimal Designs of Parameters of Structural Topology Optimization by Orthogonal Experiment. Mach. Des. Res. 2011, 27, 14–17. [Google Scholar]
  31. Yang, R.; Tang, P.; Ni, X.; Xiao, L.; Zhao, H. Structural Analysis and Design of an Electric Bus Frame. J. Yancheng Inst. Technol. (Nat. Sci. Ed.) 2022, 35, 30–37. [Google Scholar]
  32. Ren, S.; Gao, A.; Zhang, Y.; Han, W. Finite element analysis and topology optimization of folding mechanism of six-rotor plant protection UAV. J. Chin. Agric. Mech. 2021, 42, 53–58+194. [Google Scholar]
  33. Yang, S. Structure Design and Optimization of Pressure Cabinof Underwater Vehicle. Master’s Thesis, Jiangxi University of Science and Technology, Ganzhou, China, 2020. [Google Scholar]
Figure 1. Schematic diagram of internal sensor size parameters in the detection system.
Figure 1. Schematic diagram of internal sensor size parameters in the detection system.
Jmse 12 01231 g001
Figure 2. Schematic diagram of cylindrical shell’s cylindrical coordinate system.
Figure 2. Schematic diagram of cylindrical shell’s cylindrical coordinate system.
Jmse 12 01231 g002
Figure 3. Eigenvalue buckling analysis deformation diagram. (a) Instability mode of shell with equal stiffness and inner ribs; (b) instability mode of inner rib stiffness adjustment shell characteristics.
Figure 3. Eigenvalue buckling analysis deformation diagram. (a) Instability mode of shell with equal stiffness and inner ribs; (b) instability mode of inner rib stiffness adjustment shell characteristics.
Jmse 12 01231 g003
Figure 4. Mises shell stress cloud maps: (a) Mises stress cloud map of the inner rib cage with the same degree of rigidity; (b) the inner rib stiffness adjustment shell’s Mises stress cloud map.
Figure 4. Mises shell stress cloud maps: (a) Mises stress cloud map of the inner rib cage with the same degree of rigidity; (b) the inner rib stiffness adjustment shell’s Mises stress cloud map.
Jmse 12 01231 g004
Figure 5. Mises stress cloud maps of ribs: (a) Mises stress cloud map of ribs with equal stiffness; (b) Mises stress cloud map for rib stiffness adjustment.
Figure 5. Mises stress cloud maps of ribs: (a) Mises stress cloud map of ribs with equal stiffness; (b) Mises stress cloud map for rib stiffness adjustment.
Jmse 12 01231 g005
Figure 6. Diagrams of relationships between assessment indicators and various levels. (a) The relationship between different levels and maximum stress; (b) the relationship between different levels and iteration steps.
Figure 6. Diagrams of relationships between assessment indicators and various levels. (a) The relationship between different levels and maximum stress; (b) the relationship between different levels and iteration steps.
Jmse 12 01231 g006
Figure 7. Maximum stress cloud and section diagram after optimizing A3B1C3 combination. (a) Maximum stress cloud map; (b) cross-section view.
Figure 7. Maximum stress cloud and section diagram after optimizing A3B1C3 combination. (a) Maximum stress cloud map; (b) cross-section view.
Jmse 12 01231 g007
Figure 8. Maximum stress cloud and section diagram after optimizing A1B2C3 combination. (a) Maximum stress cloud map; (b) cross-section view.
Figure 8. Maximum stress cloud and section diagram after optimizing A1B2C3 combination. (a) Maximum stress cloud map; (b) cross-section view.
Jmse 12 01231 g008
Figure 9. Outline drawings of reinforced rib section after processing: (a) 3D schematic diagram of strengthening ribs; (b) 2D schematic diagram of strengthening ribs.
Figure 9. Outline drawings of reinforced rib section after processing: (a) 3D schematic diagram of strengthening ribs; (b) 2D schematic diagram of strengthening ribs.
Jmse 12 01231 g009
Figure 10. Definition diagram of arched rib parameters.
Figure 10. Definition diagram of arched rib parameters.
Jmse 12 01231 g010
Figure 11. Cross-section diagram of pressure-resistant shell.
Figure 11. Cross-section diagram of pressure-resistant shell.
Jmse 12 01231 g011
Figure 12. Sensitivity relationship curve between design variables and stress.
Figure 12. Sensitivity relationship curve between design variables and stress.
Jmse 12 01231 g012
Figure 13. Sensitivity relationship curve between design variables and buckling characteristic values.
Figure 13. Sensitivity relationship curve between design variables and buckling characteristic values.
Jmse 12 01231 g013
Figure 14. Sensitivity relationship curve between design variables and quality.
Figure 14. Sensitivity relationship curve between design variables and quality.
Jmse 12 01231 g014
Figure 15. Bar chart showing the sensitivity relationship between design variables and various output responses.
Figure 15. Bar chart showing the sensitivity relationship between design variables and various output responses.
Jmse 12 01231 g015
Figure 16. NSGA-2 process diagram.
Figure 16. NSGA-2 process diagram.
Jmse 12 01231 g016
Figure 17. Comparison between predicted and simulated stress values for pressure-resistant shells.
Figure 17. Comparison between predicted and simulated stress values for pressure-resistant shells.
Jmse 12 01231 g017
Figure 18. Comparison between predicted and simulated buckling characteristic values of pressure-resistant shells.
Figure 18. Comparison between predicted and simulated buckling characteristic values of pressure-resistant shells.
Jmse 12 01231 g018
Figure 19. Comparison between predicted and simulated values of pressure shell quality.
Figure 19. Comparison between predicted and simulated values of pressure shell quality.
Jmse 12 01231 g019
Figure 20. The trends of changes between different response values: (a) stress–buckling characteristic value; (b) quality–buckling characteristic value; (c) stress–quality.
Figure 20. The trends of changes between different response values: (a) stress–buckling characteristic value; (b) quality–buckling characteristic value; (c) stress–quality.
Jmse 12 01231 g020
Figure 21. Schematic diagram of pressure-resistant housing structure.
Figure 21. Schematic diagram of pressure-resistant housing structure.
Jmse 12 01231 g021
Table 1. Mechanical parameters of common materials.
Table 1. Mechanical parameters of common materials.
PerformanceStructural Steel Q345Titanium Alloy TB9Aluminum Alloy 6061T6Resin-Based Carbon Fiber T700Silicon Nitride CeramicsGlass
Elastic modulus (GPa)21010069120300100
Density (g/cm3)7.84.52.71.63.22.32
Poisson’s ratio0.280.310.330.30.270.2
Yield limit (MPa)34597528060038201450
Specific strength (MPa/(g/cm3))441721043751194650
Specific rigidity (GPa/(g/cm3))272226759445
Table 2. Distribution of moment of inertia of rib section.
Table 2. Distribution of moment of inertia of rib section.
Rib Number1 and 62 and 53 and 4
Inertia of rib section with plate/mm44902.176923.249675.33
Rib height/mm42.6432.9663.316
Table 3. Orthogonal experimental factor level table for three parameters.
Table 3. Orthogonal experimental factor level table for three parameters.
HorizontalVolume Coefficient (A) Penalty Factor (B) Grid Size (C)
10.42.00.5
20.53.00.75
30.64.01.0
Table 4. Orthogonal topology optimization experimental scheme combination table.
Table 4. Orthogonal topology optimization experimental scheme combination table.
Test NumberVolume Coefficient (A) Penalty Factor (B) Grid Size (C)Blank (D)
10.42.00.5-
20.43.00.75-
30.44.01-
40.52.00.75-
50.53.01-
60.54.00.5-
70.62.01-
80.63.00.5-
90.64.00.75-
Table 5. Comparing cross-sectional topological forms with various combinations of parameters.
Table 5. Comparing cross-sectional topological forms with various combinations of parameters.
Test Number123456789
Combination numberA1B1C1A1B2C3A1B3C3A2B1C2A2B2C3A2B3C1A3B1C3A3B2C1A3B3C2
Cross-sectional figureJmse 12 01231 i001Jmse 12 01231 i002Jmse 12 01231 i003Jmse 12 01231 i004Jmse 12 01231 i005Jmse 12 01231 i006Jmse 12 01231 i007Jmse 12 01231 i008Jmse 12 01231 i009
Maximum stress (MPa)267.74285.59259.47241.85259.24269.01221.72239.14249.34
Iteration steps181312181315161516
Table 6. (a) Range analysis table for maximum stress. (b) Range analysis table for iteration steps.
Table 6. (a) Range analysis table for maximum stress. (b) Range analysis table for iteration steps.
(a)
Experimental IndicatorsTest NumberA Volume CoefficientB Penalty Factor C Mesh SizeD Blank
Maximum stressK1812.8731.31775.89776.32
K2770.1783.97776.78776.32
K3710.2777.82740.43740.46
k1270.93243.77258.63258.77
k2256.7261.32258.93258.77
k3236.73259.27246.81246.82
Range R34.2015.5012.1211.95
Primary and secondary orderA > B > C > D
Excellent levelA3B1C3D3
Optimal combinationA3B1C3
(b)
Experimental IndicatorsTest NumberA Volume CoefficientB Penalty FactorC Mesh SizeD Blank
Iteration stepsK143524847
K246414744
K347434145
k114.3317.331615.67
k215.3313.6715.6714.67
k315.6714.3313.6715
Range R1.333.672.331
Primary and secondary orderB > C > A > D
Excellent levelA3B1C3D3
Optimal combinationA1B2C3
Table 7. Comparison table of data from different cross-sections.
Table 7. Comparison table of data from different cross-sections.
Maximum Stress/MPaMaximum Displacement/mmMass/g
Reinforced ribs with rectangular cross-section160.200.67037.88
Strengthening ribs of the original bow section221.720.10520.43
Reinforced rib of processed arched section199.510.08425.45
Table 8. Related parameters of ribs.
Table 8. Related parameters of ribs.
NumberInertia Moment of Rib Section with Plate/mm4Cross-Section Area/mm2Rib Width/mmRib Height/mm
1 and 65633.38819.679.12.73
2 and 57955.92523.339.912.97
3 and 411118.5327.610.783.23
Table 9. Initial shell-related data.
Table 9. Initial shell-related data.
NameNumberUnit
Shell length350mm
Shell thickness5mm
Inner diameter of shell90mm
Front bulkhead width45mm
Rear bulkhead width45mm
Rib spacing52mm
Number of ribs6Number
Table 10. Design variable range for pressure resistant shell.
Table 10. Design variable range for pressure resistant shell.
Design Variablet/mmh/mmb1/mmb2/mmb3/mmL/mm
Initial value15159.19.9110.7852
Value range7.5~17.510~255~155.5~15.56~1640~58
Table 11. Sensitivity of each design variable to the corresponding output.
Table 11. Sensitivity of each design variable to the corresponding output.
Response Sensitivitytb1b2b3Lh
Stress0.69830.03970.09370.02750.03710.0837
Buckling0.71050.02120.08130.15110.04460.0827
Quality0.71680.02090.02190.0229300.2196
Table 12. Sample points and output response values.
Table 12. Sample points and output response values.
GroupShell Thickness t/mmRib Width b1/mmRib Width b2/mmRib Width b3/mmRib Spacing L/mmEnd Cover Thickness h/mmStress/MPaBuckling Characteristic Value QQuality M/kg
17.756.7512.2510.0845.5522.38935.431.5015.604
28.0812.2513.7511.9252.1517.13878.711.5375.634
38.429.258.588.0848.5521.88850.062.8275.954
48.7514.2510.7513.9249.1524.88871.172.1625.861
59.0812.0814.258.2549.4518.88820.572.1215.859
69.429.425.7514.7548.8522.63775.312.5525.749
79.757.586.9210.2546.4510.13795.272.2525.795
810.0811.428.2510.4240.1520.38858.232.3155.696
910.427.7511.926.7548.2522.13760.143.0965.721
1010.7513.7511.0811.7543.4512.63704.393.836.137
1111.085.088.7514.2542.8524.13688.074.3525.983
1211.428.756.2510.5841.3513.88686.144.7136.213
1311.7510.0814.757.0850.3523.63690.134.2266.574
1412.0813.927.927.7550.6523.38842.724.1476.565
1512.4211.5812.427.2544.3522.88603.134.2746.787
1612.756.9211.2515.9256.6518.63650.924.3387.097
1713.085.2513.5815.2551.5516.38817.614.847.210
1813.426.4212.928.4246.1519.38706.074.5857.340
1913.7511.089.9212.7544.9514.13604.795.7356.803
2014.089.087.4213.5853.0512.38599.775.3927.528
2114.428.256.087.9255.4517.88556.996.1527.001
2214.756.0814.0810.9241.6511.13571.285.9957.569
2315.088.9210.0810.7552.7512.88555.687.1017.579
2415.427.429.4213.4242.5516.13543.096.9957.027
2515.7512.429.756.5846.7517.38556.507.4197.897
2616.0814.757.2513.0844.0519.88563.367.6527.791
2716.4210.4211.4215.5841.9521.38528.728.628.263
2816.7514.427.5815.4247.6524.63523.158.8758.066
2917.089.927.089.7549.7523.13574.6910.1597.962
3017.425.5814.428.9245.2523.88500.1210.3178.286
Table 13. Approximation model accuracy evaluation indicators.
Table 13. Approximation model accuracy evaluation indicators.
Evaluating IndicatorStressBuckling EigenvalueQuality
Determination coefficient R20.996570.999970.99999
Root mean square error0.097750.014271.232 × 10−10
Relative maximum absolute error (%)4.75821.63430.092175
Relative average absolute error (%)4.62450.439330.0096425
Table 14. Candidate points for optimizing pressure-resistant shells.
Table 14. Candidate points for optimizing pressure-resistant shells.
Design VariableShell Thickness t/mmRib Width b1/mmRib Width b2/mmRib Width b3/mmRib Spacing L/mmEnd Cover Thickness h/mm
Initial value159.19.9110.785215
A11.0005.4975.52115.65447.24813.465
B11.0936.9335.67913.27247.32913.145
C11.0199.8165.52512.86145.91413.062
D11.0668.5605.88413.30045.52013.137
E10.93511.2785.61815.46744.53713.113
Table 15. Comparison results between each candidate point and the original prototype.
Table 15. Comparison results between each candidate point and the original prototype.
ResponsesStressBuckling EigenvalueShell Mass (with End Cap)/kg∆M/kg
Original prototype549.335.85627.237
A644.264.36776.01081.2262
B646.734.34645.99881.2382
C649.124.30515.98781.2492
D647.174.34776.00291.2341
E645.884.39686.01161.2254
Table 16. Final parameters of pressure-resistant shell after rounding processing.
Table 16. Final parameters of pressure-resistant shell after rounding processing.
Design VariableShell Thickness t/mmRib Width b1/mmRib Width b2/mmRib Width b3/mmRib Spacing L/mmEnd Cover Thickness h/mm
Numerical value11106134513
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Huang, Y.; Wang, X.; Liu, Z.; You, Y.; Ma, H. Design of Inner Ribs with Unequal Stiffness for Deep-Sea Highly Pressure-Resistant Cylindrical Shells and Utilizing NSGA-2 for Lightweight Optimization. J. Mar. Sci. Eng. 2024, 12, 1231. https://doi.org/10.3390/jmse12071231

AMA Style

Huang Y, Wang X, Liu Z, You Y, Ma H. Design of Inner Ribs with Unequal Stiffness for Deep-Sea Highly Pressure-Resistant Cylindrical Shells and Utilizing NSGA-2 for Lightweight Optimization. Journal of Marine Science and Engineering. 2024; 12(7):1231. https://doi.org/10.3390/jmse12071231

Chicago/Turabian Style

Huang, Yizhe, Xiao Wang, Zhiqiang Liu, Ying You, and Haoxiang Ma. 2024. "Design of Inner Ribs with Unequal Stiffness for Deep-Sea Highly Pressure-Resistant Cylindrical Shells and Utilizing NSGA-2 for Lightweight Optimization" Journal of Marine Science and Engineering 12, no. 7: 1231. https://doi.org/10.3390/jmse12071231

APA Style

Huang, Y., Wang, X., Liu, Z., You, Y., & Ma, H. (2024). Design of Inner Ribs with Unequal Stiffness for Deep-Sea Highly Pressure-Resistant Cylindrical Shells and Utilizing NSGA-2 for Lightweight Optimization. Journal of Marine Science and Engineering, 12(7), 1231. https://doi.org/10.3390/jmse12071231

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop