Next Article in Journal
Can a 16th Century Shipwreck Be Considered a Mercury Source in the 21st Century?—A Case Study in the Azores Archipelago (Portugal)
Next Article in Special Issue
Epoxy Resin Biocomposites Reinforced with Flax and Hemp Fibers for Marine Applications
Previous Article in Journal
Experimental and Numerical Investigation of Cavity Structure Forced Water Exit from Calm Water at Constant Lifting Velocity
Previous Article in Special Issue
Multi-Scale Analysis for Assessing the Impact of Material Composition and Weave on the Ultimate Strength of GFRP Stiffened Panels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Finite Element Simulations of Novel Submersible Design Based on the ASME Design-by-Analysis Approach

1
Maritime Department, University of Dubrovnik, 20000 Dubrovnik, Croatia
2
Marine and Energy Solutions DIV d.o.o, 10000 Zagreb, Croatia
3
Faculty of Mechanical Engineering and Naval Architecture, University of Zagreb, 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2023, 11(2), 275; https://doi.org/10.3390/jmse11020275
Submission received: 19 November 2022 / Revised: 4 January 2023 / Accepted: 13 January 2023 / Published: 25 January 2023

Abstract

:
Tourist submersibles are a new class of underwater vehicles where new materials and innovative structural concepts are used. Because of the lack of experience with these innovative vessels, prescribed rules of classification societies are not directly applicable in their structural design. Thus, a novel design-by-analysis approach, as given in the Boiler and Pressure Vessel Code of ASME, is employed for dimensioning structural components. Design-by-analysis stands for a consistent set of loads, boundary conditions, modelling principles, failure modes, and safety factors if the finite element analysis is performed instead of using rule formulas. The approach has been used until now in the structural design and optimization of unconventional pressure vessels, while in this paper it is for the first time employed in structural design of the underwater vehicle. This study can serve as guidance for future applications of design-by-analysis in the development of the innovative structures exposed predominantly to the pressure loading.

1. Introduction

Submarines have traditionally been constructed of metallic materials and used by navies for military purposes. This has had a profound influence on their design—the crew is confined within the submarine hull having no direct visual contact with the surroundings. Such a design decision is reasonable since military submarines operate at ocean depths much greater than 200 m by which more than 99% of sunlight is absorbed [1].
When a submarine is built for tourist purposes, one of the main goals is to provide passengers with a view of the underwater world. Thus, the hull must be fitted with openings allowing as broad and as clear view as possible since unobstructed views from a manned submersible is a major selling point for any non-military application. This means that two major challenges must be overcome. The first is related to assessing and taking care of the hull stiffness lost by the openings and the second one is related to selecting a material that combines good visibility and sufficient strength, among other properties.
The first successful submarine concept with hull openings was that of Auguste Piccard in the 1960s. He constructed and used a conical frustum window made of acrylic plastic [2]. Due to the window’s shape, external pressure is distributed more favourably along the window’s connection with the hull structure and, as a result, this type of window is still widely used and investigated, especially for the greatest of ocean depths. Two recent examples are a study on failure modes of acrylic conical frustum windows of the Jialong, a Chinese submersible intended for hadal zones [3], and a study of Pranesh et al. [4] on design of conical frustum fillets to prevent stress concentration.
The problem with the window in the context of a tourist submarine is, however, that its thickness increases rapidly with an increase in the diameter, even at a moderate depth [2], which limits the size of the opening. To solve this problem, a spherical window has been proposed [2]. It allows for greater diameters while preserving good visibility. However, this produces only a hull with “peepholes” rather than offering a panoramic view from the hull interior. To achieve that, the hull itself must be made of a transparent material.
The appropriate material for this purpose is polymethyl methacrylate (PMMA) or acrylic plastic, which has been proved in a monumental study of Jerry Stachiw and his co-workers [2]. The first submersibles having the entire hull made of acrylic plastics and intended for tourist purposes seems to be those of Hyco Technologies of Vancouver, Canada. They developed two concepts in the mid-to-late 1980s [5]. The first, named Gemini, represented a submersible that comprised three acrylic spheres embedded into the framework made of stainless steel [5]. Although the spheres provide an unobstructed view of almost 360°, the concept has some serious drawbacks, especially if the vessel is intended for a much larger number of people. The potential capacity is low since the size of an acrylic sphere is limited (due to manufacturing reasons [2]). Furthermore, it has poor hydrodynamical properties and difficult communication between the spheres (in the case of the Gemini, communication was possible by means of “transfer tunnels” fitted between the spheres [5]).
The second concept of Hyco Technologies, called Aries [5], emerged in the 1980s and was patented in 1990 [6]. It is based on the conventional, cylindrical submarine hull shape but since building the entire hull of acrylic is technologically challenging [2], a solution to the problem was found by fitting the cylindrical acrylic sections between supporting steel rings connected to an external steel structure (the exo-structure). Thus, the steel structure can take most of the load while the acrylic cylinders can provide the unobstructed view of the surroundings. The first submersible of this concept, called Seabus, was built in 1992 and was classed by the American Bureau of Shipping (ABS). To the best of our knowledge, there has been only one more submersible of this concept built so far, called Deepview.
Submersibles of the “Seabus concept”, as it will be referred to here, bring several advantages compared with other similar submarines. They offer an almost unobstructed view from the hull interior. Compared with the spherical type of submersibles, the communication between all parts of the hull is seamless. The modular concept of their hull offers great flexibility in terms of their size and, consequently, number of passengers. The innovative structural arrangement of the Seabus concept, however, brings many design challenges in assessing the structural behaviour of the submarine. In this regard, classification societies rules can be considered. About 2/3 of manned submersibles are classed by three major classification societies, ABS (21%), Det Norske Veritas (DNV) (13%), and Lloyd’d Register (LR) (31%) [7]. ABS is the only one of the three that has rules dedicated to civilian submersibles called Rules for Building and Classing Underwater Vehicles, Systems, and Hyperbaric Facilities [8]. The rules are concerned with the assessment of the metallic parts of a submersible’s structure while the acrylic parts are assessed according to ASME’s code on pressure vessels for human occupancy. Other classification societies offer naval rules that focus on military submarines and thus tend to have more conservative criteria resulting in heavier structures [9].
Although the ABS rules are better suited than the naval rules for the submersible concept dealt with here, a structural design approach based entirely on classification societies’ rules would still include all the well-known cons of such an approach, namely, its inability to achieve an optimum design and its questionable applicability to novel vessel concepts [10].
The limitations of such an approach can be overcome if the structural behaviour is determined by means of a direct strength analysis method (e.g., the finite element method, FEM). This allows for lower structural safety factors since the structural behaviour is assessed more accurately. In addition, the analysis can readily be incorporated into a broader design framework, often referred to as the rational-based design [11], which, together with an optimization tool, can produce an optimum structural arrangement.
To leverage the benefits of the direct strength analysis, an approach by the American Society of Mechanical Engineers (ASME), called the design-by-analysis (DBA), is investigated in this paper. The methodology is a part of ASME’s Boiler and Pressure Vessel Code (BPVC) [12].
The code that first appeared in the 1960s was intended for nuclear facilities. Subsequently, new parts have been added facilitating its application in a broad range of structures that need to withstand a pressure gradient directed inwards (external pressure) or outwards (internal pressure). The pressure hull of the submersible will, in this regard, be considered as a pressure vessel under the external pressure and application of the ASME’s code on the metallic part of the structure will be investigated in the following sections.
There are two ways of demonstrating that a pressure vessel fulfils the BPVC criteria [12]. The first one is the design-by-formula (DBF) approach, which shares all the pros and cons of the approach based entirely on classification societies’ rules [13]. The other approach called design-by-analysis (DBA), undertaken here, is based on direct strength analysis by means of FEM and can potentially result in lighter and safer structures, which has been demonstrated in many studies that deal with shore-based pressure vessels.
In their analysis of a vertical pressure vessel with skirt support, Diamantoudis and Kermanidis demonstrated [14] that the maximum allowable pressure on the structure is as much as 40% higher for mild steel and 70% for a high-strength steel when the DBA approach is used instead of the DBF approach. An even greater difference is reported in [15] where the allowable pressure is 2.32 times greater when the DBA approach is used for the mild steel and 2.5 times greater when the DBA approach is used for the high-strength steel compared with the results obtained by means of the DBF. The main differences between the two studies are that the elastic-plastic method of [12] is used in [14], while the limit state method is used in [15] and the results are compared at different positions, which might explain the discrepancy.
Similar conclusions can be derived from studies concerned with marine structures. In [16], the structural behaviour of high-strength steel hemispheres is determined by means of physical experiments and numerical models. The obtained critical collapse pressures are then compared with the pressures determined according to the rules of classification societies, including the Underwater Vehicles Rules of ABS [8]. The pressures determined based on FEA differ up to 8% from the experimental results, while the critical pressures obtained by means of ABS’s rules are between 28% and 36% greater than the experimental results, which implies the margin available for potential savings when the direct strength assessment methods are used.
An analogous study [17] on ring-stiffened steel cylinders, conducted by the same research group, showed that the collapse pressure determined based on FEA differ from the experiments by up to 14%, while the ABS rules give a conservative collapse pressure assessment of as much as 29%.
An interesting study of Pranesh et al. [18] is concerned with structural design and analysis of a spherical steel submersible for three persons. The authors calculate the collapse external pressure for the submersible of an internal diameter of 2100 mm and shell thickness of 25 mm. The pressure is 8.54 MPa and 11.47 MPa obtained according to the ABS’s rules and nonlinear FEM, respectively, for the same safety factor of 1.5. The potential material and weight gain can be seen more directly in the part of the study where the pressure is held constant while the required shell thickness is calculated. The required shell thickness for the collapse pressure of 7.5 MPa and the stated factor of safety is 22.15 mm according to the ABS’s rules and 17 mm according to the nonlinear FEA.
Two recent studies have used ASME’s BPVC in submarine/submersible design. Ma et al. [19] used the DBF approach of ASME’s BPVC in the design of an innovative subsea shuttle tanker concept. The tanker, similar to the Seabus concept dealt with here, has metallic fore and aft peaks while the central portion of the hull is flooded when the tanker is submerged. The central portion comprises cargo, compensation, and other tanks that are designed according to the DBF approach. The peaks are designed according to the DNV’s naval code.
In the study of Sohn et al. [20], ASME’s BPVC for composites is used in assessing the collapse pressure of a submarine hull made of composite materials. The results were then compared with FEA results and the design-based collapse pressure was significantly lower for all test configurations than the pressure determined based on the direct strength calculations.
Within the context of the rationally based design, the DBA of a pressure vessel can be paired with a structural optimization tool leading to significant savings in terms of structural weight and the material used. In [21], a steel internal pressure vessel is optimized by connecting ANSYS, structural analysis software, with an optimizer in Matlab. The structural analysis is based on linear elastic FEM where the optimization problem is defined as follows. Thicknesses of the pressure vessel shell are defined as design variables that are varied to minimize structural weight while the optimization constraints are based on the ASME BPVC requirements for DBA following the same procedure that will be demonstrated in this paper. It is demonstrated that the final structure has the structural weight reduced by almost 40%.
The paper is structured as follows. After the Introduction, Section 2 elaborates on the proposed procedure and shows its departure from the rulebook approach of ABS and presents relevant characteristics of a tourist submersible concept with emphasis on the pretension force. Section 3 demonstrates the workings of the procedure on structural analysis of the submersible. For the sake of completeness, Section 4 provides a brief discussion on the main aspects of non-metallic parts of the pressure hull, which are not dealt with in detail in this research. The paper ends with relevant conclusions.
DBA has been used so far in the structural design of unconventional offshore pressure vessels, while in the present paper it is used for the first time for dimensioning of structural components of the underwater vehicle. The main issue in the structural design of innovative structures is to find the right balance between applied design loads, finite element mesh size, and safety factors for different failure modes, which is consistently proposed within the DBA approach. This study may therefore be used as a guidance and incentive for future application of DBA in the development of innovative structures predominantly exposed to pressure loading.

2. Design Requirements of the Acrylic Hull Submersible

This chapter consists of two subchapters. In the first subchapter, the proposed computational procedure based on the design-by-analysis (DBA) approach of ASME is described. In the second subchapter, the design of a passenger submersible used to demonstrate the workings of the procedure is presented.

2.1. A Computational Procedure for Structural Analysis of a Passenger Submersible Based on the Design-by-Analysis (DBA) Approach by ASME

Since both submersibles with cylindrical acrylic hulls mentioned in the Introduction have been classed by ABS, the procedure prescribed by the stated classification society will be referred to as the standard procedure for the submersible design. The procedure (Figure 1) starts with ABS’s rules [22] that require, among other things, performing “pressure vessel stress analysis” as given in “Structural evaluation of metallic pressure boundary components (Chapter 6) of [8], “analysis of lifting load and stresses induced in the hull” and “window calculations” in compliance with “Structural evaluation of windows and viewports” (Chapter 7) of [8]. In the context of the vessel described in this paper, it means that two kinds of analyses must be performed. First, all metallic parts of the pressure hull and supporting sub- and super-structure (exo-structure) must be analysed for two load cases: (1) the structure is submerged and behaves like a pressure vessel and (2) the structure is lifted to be either removed from the water or launched. The second type of analysis is related to evaluation of cylindrical acrylic hull segments and is briefly discussed in the paper.
Based on the foregoing, it can be stated that the standard design procedure starts with ABS Rules for Conditions of Classification [22] followed by the Rules for Underwater Vehicles [8] as presented in Figure 1. ABS rules [8], however, leave open an option for “other recognised standards”, which is the primary aim of this paper. Thus, structural evaluation of the metallic portion of the pressure hull and the exo-structure will be conducted by means of the ASME DBA procedure [12], which is denoted by a dashed arrow in Figure 1.
The DBA of ASME is concerned with providing protection against four failure modes: plastic collapse, local failure, buckling, and cyclic loading [12].

2.1.1. Protection against Failure from Plastic Collapse

When the structure is subjected to loads, strains and stresses arise within the structure. As the load increases, the stresses reach the yield limit of the structure, and the first plasticity regions emerge. Followed by further increase in loading, plasticity regions can become so large that the structure loses its load-carrying ability, i.e., plastic collapse occurs. A designer has three methods [12] to evaluate protection against this failure mode: the elastic stress analysis method, the limit-load method, and the elastic-plastic stress analysis method.
The present paper demonstrates the first of the three methods. It is the most widely used since it is based on a simple elastic structural analysis that can easily be performed by means of the finite element method (FEM) in few steps as presented in Figure 1.
The first step consists of identifying all relevant loads on the structure and deriving corresponding load cases (i.e., grouping loads that act simultaneously on the structure) together with adequate boundary conditions. In the problem considered in this paper, the loads are caused by hydrostatic pressure from the surrounding water and self-weight of the structure.
Stresses obtained after structural analysis will depend on the type of finite element (FE) used. The ASME’s code [12] allows for 3D continuum mechanics (“solid”) finite elements and (2D) shell FEs where the choice of the FE will affect the stress categorisation required by the code.
Thus, the code [12] distinguishes between the following stress categories (relevant for the protection against plastic collapse approach based on the elastic stress analysis method): general membrane stress (Pm), local membrane stress (PL), and bending stress (Pb). If shell FEs are used, the stress categorisation is straightforward since a shell FE can output membrane and bending stress separately [12,23] (unless specified otherwise) when the shell is built out of a single homogenous isotropic material and a linear static behaviour is considered. If a 3D continuum mechanics (“solid”) FE is used, a process called linearisation [12] must be performed to obtain the stresses from a general stress tensor. Nowadays, there are many stress categorisation algorithms that can readily be integrated into commercial FEA programs, although the process is rather cumbersome and is still an active field of research [24].
Once the stresses are categorised into membrane and bending stress, it remains to distinguish between general and local membrane stress. The difference is given in the code where it is stated that if a membrane stress surpasses 1.1 of the limit stress (to be clarified in a subsequent paragraph) and when it is confined in a region not greater than R t (where R is the mean radius of the structure and t is its thickness), the membrane stress can be categorised as the local membrane stress. Otherwise, it is considered as the general membrane stress.
The next step consists of calculating equivalent stresses for each category. In this approach, equivalent stresses based on the von Mises yield criterion are required [12].
The last step includes determining allowable stresses and comparing them with the equivalent stresses for each of the categories. There are two allowable stresses prescribed by the code [12]: allowable stresses (denoted by S) used to evaluate general membrane stress and allowable stress (denoted by SPL) used to evaluate local membrane and local membrane and bending stress taken together. Thus, the stress limits can be written as follows:
P m < S
P L < S P L
P L + P b < S P L
S and SPL are functions of the material’s yield stress and are related in the following manner: S P L = 1.5 S .

2.1.2. Protection against Failure from Buckling

In addition to the evaluation of protection against plastic collapse, a design factor for protection against buckling collapse should be satisfied for compressive stress field. The design factor depends on the type of buckling analysis performed (i.e., bifurcation buckling analysis or elastic–plastic collapse analysis). Buckling design factors in the present study are based on “Type 1” factor level, since an elastic buckling analysis is performed without geometrical nonlinearities. Thus, based on [12] a minimum safety factor against buckling reads ϕB = 2/βcr. The capacity reduction factor βcr for unstiffened and ring stiffened cylinders and cones under external pressure reads 0.8, while for spherical shells and spherical, torispherical, and elliptical heads under external pressure βcr are equal to 0.124 [12]. As the boundary conditions are critical for the buckling analysis, a special attention must be paid to avoid that boundary conditions prevent some realistic buckling modes.

2.1.3. Protection against Failure from Cyclic Loading

Protection against failure from cyclic loading consists of assessing the structure for fatigue and ratchetting. These load effects result from three additional stress categories besides the bending stress (Pb) elaborated in Section 2.1.1. These are the primary local membrane (PL), the secondary membrane and bending stress (Q), and peak stress (F). These stress components occur due to structural discontinuities (e.g., change in plate thickness), thermal gradients due to temperature changes, stress concentrations (e.g., junctions, welds, sharp edges), etc. Thus, the total stress to be considered in the fatigue analysis is [12]:
P L + P b + Q + F
In general, the assessment should be carried out if the structure is subjected to cyclic loading [12]. However, if a set of prescribed criteria is fulfilled, no assessment is required. Thus, before starting the fatigue assessment, the fatigue analysis screening is performed (Figure 1).
There are three approaches to the screening. The first one is based on the “previous experience with comparable equipment” [12] where the designer can prove that the fatigue assessment is not necessary if there is a sufficiently long history of use of similar equipment in similar loading conditions. The other two approaches include a series of steps that must be followed if relevant data for the first approach are not available, which is the case here.
There are two methods available in this regard: method A and B. The former one is used in this research because it is applicable to materials with tensile strength below 552 MPa. The goal of the method is to determine the total number of cycles of all alternating loads that act on the structure and compare it with a relevant limit (1000 for the integral and 400 for the nonintegral construction). Since the total number of cycles in this research is 25,000 (see Section 3.2.4), the fatigue assessment must be performed.
The fatigue assessment [12] consists of steps given in Figure 1. The appropriate stresses can be obtained by means of either an elastic or an elastic–plastic analysis. The former will be used here similarly to the plastic collapse evaluation of the previous subchapter.
The first step in the assessment procedure is to determine load history, i.e., to identify all alternating loads in the operating condition. The second step is concerned with determining relevant stress–strain cycles from the load histogram of the previous step. Several cycle counting methods are available (e.g., the rainflow method of ASTM). The total number of the cycles for a structural component of interest is denoted by M.
Once the cycles are determined, a structural analysis is conducted for each cycle (k) to obtain corresponding stresses. This step is similar to the analysis conducted in the plastic collapse assessment and the same tools can be used (e.g., FEM). It is, however, important to note that, as mentioned in the first step, the loads correspond to the operating conditions rather than to the design conditions. The step is conducted for each cycle determined in the previous step.
The fourth step is concerned with determining the alternating stress amplitudes. Two expressions are provided in [12] depending on whether the stresses arise only from mechanical loads (which is the case in this research), or they are caused by temperature changes as well.
Based on the stress amplitudes determined in the previous step, the maximum permissible number of cycle repetitions (Nk) can be obtained from diagrams provided in [12]. This allows us to calculate the fatigue damage in the current (k-th) cycle:
D f , k = n k N k
where nk is the number of repetitions of the k-th cycle.
The structure is considered safe to fatigue if the accumulated fatigue damage for all cycles is less than 1 [12]:
D f = k = 1 M D f , k 1
The number of cycles in a ship design is typically determined by the histogram of wave-induced stresses. Since there is only one bin of the histogram in this case (as will be shown in Section 3.2.4), we can write M = 1. By introducing this and (5) to (6), we obtain the fatigue assessment criterion for the structure at hand:
D f = n 1 N 1 1
where n1 and N1 denote the number of repetitions of the first (and only) cycle/loading condition and the permissible number of cycle repetitions, respectively.

2.2. Submersible Design Details

The previously described computational procedure is demonstrated on a problem of structural analysis of a passenger submersible. The main particulars of the submersible are given in Table 1, while its general arrangement is presented in Figure 2.
While tourist submersibles are not intended for large speeds, their hydrodynamics resistance is still important to ensure enough battery power for sufficient cruising autonomy and to resist sea currents. Thus, hydrodynamic analysis of the bare hull of the tourist submarine moving in forward and transverse directions is performed [25]. External structures required to ensure the structural integrity of acrylic hulls considerably increases hydrodynamic resistance, as described by Gatin et al. [26]. Additionally, all hydrodynamical aspects, such as submersible resistance, seakeeping, and propeller performance are planned to be tested in the towing tank and cavitation tunnel.
The structure of the tourist submersible consists of a pressure hull and an exo-structure. The pressure hull is composed of the steel heads and steel transverse rings that support the acrylic cylindrical sections. The aft head is completely made of steel, and it is partially cylindrical, conical, and spherical at the tip. The fore end is made partially of steel (cylindrical and conical parts), while the foremost part is a hemispherical acrylic window for the piloting of the submarine.
As the mechanical properties of the acrylic are poor and rather sensitive to the axial stresses and temperature variation, the basic design requirement is that acrylic cylinders should not contribute to the structural strength during the diving of the submarine. Instead, the exo-structure is designed to bear longitudinal loads produced by hydrostatic pressure acting on the submarine heads, thus preventing axial load transfer to the acrylic cylinders. The exo-structure is bolted to the aft and fore heads and transverse rings, thus enabling maintenance and replacement of the acrylic cylinders which carry only local loads and transmit them to the metallic structure through the supporting rings. Tie rods connecting aft and fore steel heads are used for fastening cylindrical acrylic sections. Although the role of the tie rods needs to be investigated further, they are used here to keep the acrylic cylinders in place during manufacturing and to account for shortening of the acrylic cylinders with a decrease in temperature.
The important aspect of the submersible structural design is the pretension force in tensioning rods. Although determination of the tension force in the tie rods will be dealt with in a separate paper, it will be demonstrated here how a preliminary assessment of the force in the tie rods can be obtained by means of the theory of elasticity and requirements of [27], which is concerned with the design of acrylic hull sections.
The tie rods have two purposes here. The first one is to keep the acrylic cylinders in place during manufacturing and the second one is to account for shortening of the acrylic cylinders with a decrease in ambient temperature.
Once the hull sections are manufactured and mounted to the exo-structure, their position is fixed with respect to the supporting rings. Between the cylinders and the rings lies a joint that consists of a seal (that prevents water from entering the hull) and two gaskets (made of chloroprene rubber), which provide a seating for the cylinders. Design of the joint can be seen in Figure 3 and is based on [2]. As the temperature decreases, the length of the acrylic cylinder will decrease as well. This could lead to misalignment of the parts in the joint. To account for this, the tie rods can be pre-tensioned such that they keep the cylinders always adjacent to the side gasket. In order to achieve this, it is important to: (1) calculate longitudinal displacements at the ends of the cylinders due to their thermal dilatation, (2) calculate the pressure on the side gasket that causes the displacement obtained in 1, and (3) determine the compression force that produces the pressure from the previous step (magnitude of the tension force in the tie rods is equal to the compression force that keeps the gasket in place).
The procedure presented here is based on the linear theory of elasticity with the following assumptions: the deformation is considered small, the strains infinitesimal, and the material behaviour is linear elastic. These assumptions do not hold since the strain is almost 5% (as will be seen) and the behaviour of chloroprene rubber (CR) is, generally, non-linear elastic [28]. Thus, a more accurate approach would be based on the theory of hyperelasticity considering viscoelastic properties of the rubber. The procedure presented here is proposed for a first assessment because of the following: (1) behaviour of CR is almost linear up to approximately 0.2 engineering strain in pure shear, (2) although the problem at hand is axisymmetric, it can be shown that a plane strain analysis of the problem is sufficiently accurate (the difference in stress resultants is less than 5%) [29]. Since a state of strain and stress in a pure shear/planar tension test is similar to the state of stress and strain under plane strain assumptions, the material behaviour as observed in the test can be used in this calculation. Thus, the modulus of elasticity of CR will be held constant and calculated according to data given in [30] for the pure shear test up to 0.08 engineering strain.
Shortening of the cylinder due to temperature decrease can be expressed as:
Δ L = α · L · Δ ϑ = 0.89   mm
where α = 7.4 × 10−5 mm/ (mm °C) is a coefficient of linear thermal expansion at 20 °C calculated according to data obtained from [27], L is length of the cylinder (Table 2), and Δθ = −5 °C is the temperature decrease. Both temperature values are selected only to illustrate the procedure and will depend on ambient temperature in areas where the submersible is assembled and operated, respectively.
Strain in the side gasket as a consequence of ΔL:
ε = Δ L 2 L = 0.044
Stress in the gasket due to strain ε:
σ = E ε = 0.29   MPa
where E = 6.56 MPa is the modulus of elasticity of chloroprene rubber (CR).
Total compression force needed to produce stress σ in the side gasket (It is assumed that: (1) only a contact pressure acts along the side gasket edge adjacent to the acrylic cylinder and the edge adjacent to the steel ring, (2) the side gasket can expand freely in the other two directions. Thus, a stress component perpendicular to the edge must be equal to the contact pressure that produces the compression force of equation (11)):
F = σ · A = 314.79   kN
where A is the area of the side gasket in contact with the acrylic cylinder. The area can be calculated using d and h given in Table 2. Tension force in a tie rod:
F T R = F n = 39.35   kN
where n is the number of tie rods. Based on the result obtained in (12), the tension force per tie rod used in direct strength calculations in this paper is 40 kN.
In the present work, structural components such as transverse rings, fore and aft end section, and battery cylinders are analysed as separate models for yielding and buckling. Next, the steel part of the submersible structure is studied as a single unit. In the case of the yielding and buckling check, all models are analysed for the test depth of 65 m while the fatigue analysis of the complete model is performed for design depth of 50 m. Furthermore, the vessel is analysed for the three lifting conditions as well.

3. FEM Structural Analysis

The elastic stress analysis is performed by employing the FEM approach with 2D shell elements. The analysed submersible was designed based on ASTM 516 Gr65 steel which, according to the ASME in Section 2, Part D, has the maximum allowable stress of 161 MPa. Furthermore, a high tensile steel (ASTM 841 B) is used for lifting lugs which are modelled by 3D solid elements to account for high stresses which occur due to emergency lifting conditions. Therefore, for the components which are modelled by shell elements and mild steel, two allowable primary membrane stresses, depending on their location, have been conservatively established at S = 150 MPa and SPL = 1.5S = 225 MPa. As described in Section 2.1, membrane stresses are relevant in most parts of the structure while local primary membrane stresses are evaluated in the vicinity of the structural discontinuities. Steel parts are modelled using shell elements with mesh size t × t (t is the thickness of the structural element), while longitudinal rods are modelled by circular bar elements.
As the employed analysis is predominantly linear and the mesh convergence study is not a mandatory requirement in DBA, this was not performed within the scope of the present analysis. Furthermore, it is known that the stress concentration is decreasing for coarser mesh size and that the accuracy of the geometry modelling is reducing. The applied plate elements with mesh size t × t are well known in the maritime industry, and it is considered as the standard for fatigue analysis of thin-walled structures [31]. Such very fine mesh encompasses all geometry stress concentrations and finer mesh than t × t is not recommended in the offshore engineering practice. An alternative to such very fine mesh of plate elements is to use solid elements with at least three elements across the plate thickness. As the plate thickness is generally relatively small (about 20 mm), such a solid element mesh is considered as rather impractical for modelling the whole submarine structure. Nevertheless, in the present study solid elements are used in model areas where credible stress distribution cannot be obtained by shell elements, such as for the lifting lugs for the considered lifting conditions.
FEM analysis is performed for individual components and hull assembly for different load cases. This is necessary to determine possible stress increase due to interactions between individual components. In most cases, loads correspond to the test depth pressure of 65 m except for fatigue assessment which is analysed for the design depth of 50 m. Boundary conditions are employed to have minimal impact on the results while restraining rigid body motion. Preliminary results of the structural analysis and some complementary considerations in addition to those given in the following subsections may be found in the conference paper [28].

3.1. FEM Analysis of Submersible Structural Components

Individual components of transverse rings, fore and aft end cupolas, as well as battery cylinders are analysed as separate models. As previously emphasized, all employed boundary conditions are applied in such a way that they have minimal impact on the results while restraining the component rigid body motion.
Thus, the FE model of the ring is restrained in x-direction at supports of the longitudinal rods since that location is bolted to the exo-structure. Furthermore, four additional springs are used in the mid-span of the flanges to restrain the rigid body motion of the ring for the buckling analysis. Springs in the vertical, centreline plane have negligible stiffness in the y axis, while the same negligible stiffness is considered in the horizontal symmetrical plane in the z direction. For the component analysis of the fore and aft cupola, fixed boundary conditions are applied on the upper and lower fins since these components are bolted to the exo-structure at those locations. Finally, for the analysis of the battery cylinder, the translations in all three axes as well as rotation about the x-axis are restrained at the fore cylinder end, while transverse and vertical motions together with x-axis rotation are restrained at the aft cylinder end [32].
Yielding and buckling analyses are performed for the test depth of 65 m. Thus, the cylindrical battery pod is loaded with the external sea pressure of 6.5 bar (0.65 MPa) on the outside housing. In the case of transverse rings and aft and fore cupolas, an additional pressure of 8.8 MPa, which is the equivalent load, transferred from the acrylic surface of the half cylindrical section is considered. It is conservatively assumed that the whole pressure from the acrylic cylinders is transferred to the inner flange of the analysed component as presented in Figure 4.
For the analysis of steel heads, the preload is modelled by longitudinal forces of 40 kN per bar as calculated in the previous section and the corresponding pressure of 0.277 MPa in the opposite direction on the part of the head. Thus, the equilibrium of forces acting on the head is ensured. The results of the analysed components which satisfy previously described criteria are summarized in Table 3 while Figure 5 shows buckling results for the cylinder batter pod.

3.2. FEM Analysis of Submersible Assembly

For the structural assembly of the submersible, the yielding failure is determined for two load cases, i.e., submerged and lifting conditions while the buckling load factor is calculated for the submerged condition only. Radial and angular deformation analysis is performed for the test submerged condition according to the safety standard for pressure vessels for human occupancy—PVHO-1 [27]. Finally, fatigue analysis is evaluated for the design case of the submerged condition, i.e., 50 m.
The connection between the steel hull and the exo-structure is achieved by bolts which are modelled by beam elements. Nine bolts are used for the connection between one fin and the supporting structure while four bolts connect one lug of the web ring with the exo-structure.

3.2.1. Yielding Analysis

Pressure loads and preloads for the structural assembly are applied in an analogous way for the transverse rings and steel heads. As previously described, transverse rings are loaded with the external sea pressure of 6.5 bar (0.65 MPa) on the outside, exposed flange, while the inner flange is loaded with 8.8 MPa. In the case of aft and fore heads, the hydrostatic pressure load of 0.65 MPa is combined with the pressure due to preload (0.277 MPa) that acts on the web of the first ring and corresponding forces of 40 kN which acts on the longitudinal rods. The external pressure of 0.65 MPa also acts on the cylindrical battery housing as presented in Figure 6.
Boundary conditions for the global analysis of the whole submarine are applied on the tip of the fore and aft heads. All three motions and rotation of the x-axis are restrained at the fore end, while the transverse and vertical motions as well as the x-axis rotation are restrained at the aft head. Loads and boundary conditions are presented in Figure 6.
The results of the yielding analysis are similar to the component analysis results. The majority of structural elements satisfy the design requirement of maximum membrane von Mises stress of 150 MPa, while only local exceedance in the flanges of transverse rings connected to the heads are noticed. The maximum local primary von Mises stress of 178 MPa is well below permissible local value of 225 MPa. The resulting von Mises membrane stresses from the analysis of the structural assembly are presented in Figure 7.

3.2.2. Radial and Angular Deformation Analysis

Radial and angular deformation analysis is performed for large openings according to safety standards for pressure vessels for human occupancy [27]. Based on the PVHO-1 requirements for large openings, a maximum radial deformation rmax of the window seat at maximum internal or external pressure must be less than 0.002·Di, where Di in our case reads:
  • Internal diameter of the acrylic cylinder Di = 2360 mm; rmax = 4.72 mm.
  • Internal diameter of the acrylic fore sphere Di = 2001 mm; rmax = 4.00 mm.
In the case of angular deformation, the maximum allowable deformation should be less than 0.5°. Radial deformation can be calculated based on deformations in y and z directions, which are presented in Figure 8. The results are extracted for nodes with maximum translation on the inner flange diameter of each ring. The first node represents the location with the maximum translation on the connection between the web and flange, while Node 2 represents the corresponding node on the edge of the flange (100 mm from Node 1). Based on these deformation analysis results, it can be concluded that all maximum radial deformations on each window are less than the PHVO-1 requirements. The maximum radial deformation is obtained for the central ring as expected while the maximum angular deformations are calculated for the connection with the aft and fore cupola. All deformation requirements are satisfied.

3.2.3. Buckling Analysis

The lowest global buckling failure mode of the assembled structure is shown in Figure 9, where an adequate buckling load factor of 4.67 is achieved.

3.2.4. Fatigue Analysis

Fatigue analysis is performed for the design depth of 50 m where transverse rings are loaded with the external sea pressure of 5.0 bar (0.5 MPa) on the outside, exposed flange, while the inner flange is loaded with a pressure of 6.8 MPa, which is transferred from acrylic cylinders. Additionally, in the case of the aft and fore heads, the hydrostatic pressure load of 0.5 MPa is combined with the pressure due to the preload (0.277 MPa) that acts on the web of the first/last ring as presented in Figure 10. The same boundary conditions are assumed for the yielding analysis at the test depth.
As the yielding analysis revealed potential stress concentration locations, the fatigue analysis is performed for those details as presented in Figure 11. Two details are analysed:
  • Connection of the aft head and aft cupola ring.
  • Connection of the forehead and fore cupola ring.
The fatigue analysis is performed based on effective equivalent stress ranges considering fatigue strength reduction and fatigue penalty factor. The equivalent stress ranges are calculated based on plate bottom stress components since higher values of von Mises stress were obtained for the bottom plate compared with the top plate. The extracted plate bottom components are extrapolated by the hot spot stress approach and from extrapolated components the effective equivalent stress ranges are calculated as 57.7 MPa.
Based on the calculated effective equivalent stress ranges and effective alternating equivalent stress amplitude, the permissible cycle life based on the ASME fatigue curve of the adequate material can be determined. From Figure 12 it can be concluded that the permissible number of cycle repetitions is slightly below 2.2 × 105. Since the number of cycle repetitions is 25,000, the fatigue assessment criterion (accumulated fatigue damage) according to (7) reads:
D f = 2.5 · 10 4 2.2 · 10 5 1
Therefore, the design of both details is safe.

3.2.5. Lifting Analysis

Following the preliminary calculations, the design with eight lugs is proposed to meet the requirements of lifting conditions. Generally, lifting conditions are analysed to ensure that:
  • Stresses in the structures are acceptable;
  • Relative displacements at the ends of the cylindrical acrylic sections are satisfactorily low and will not cause damage.
First, a normal (static) lifting condition is analysed where a nonlinear analysis is employed since the nonlinear gap elements simulate contact between lifting cables and lifting lugs. During the iterative procedure of nonlinear analysis, the gap elements employ a large (almost infinite) stiffness in the compression condition and zero stiffness in tension. This ensures that the pulling force of the lifting cables is distributed on the nodes in the lifting ears.
Loads are simulated by a submersible weight that reads W = 156.4 t and a gravity constant g = 9.81 m/s2. The model is restrained in the transverse direction on the tip of the aft head and in longitudinal and transversal direction on the tip of the forehead as presented in Figure 13. Moreover, the vertical motion is restrained in the centre of the lifting ears to simulate the tension of lifting cables.
The results are presented in the following figures where Figure 14 shows the deformed model with von Mises stresses. It can be observed that maximum stresses of 46.8 MPa in shell elements and 69.9 MPa in solid elements are well below permissible values. Furthermore, displacements during lifting of the vessel can be seen in Figure 15. The maximum vertical displacement can be determined as 1.1 mm while relative displacements at the ends of the first cylindrical sections (between the first and second ring) are calculated to be 0.8 mm, which is found to be satisfactory by the criteria for acrylic windows. Namely, such a reasonable level of the relative displacements should assure that acrylic cylinders will not be damaged during frequent lifting operations with the submersible.
Besides lifting in the normal condition during launching, it is necessary to evaluate two additional lifting conditions according to ABS, which can be found in Section 6 of Underwater Rules [8]. The first one is the lifting condition from the seabed in the case when one of the battery pods is damaged.
The mass of the flooded battery housing reads 3346 kg which corresponds to the force of 32.8 kN. Since it is assumed that the forward battery pod on the port side is damaged, the upper and lower lifting lugs in the fore part of the submersible are analysed. In the first case, the lifting force due to the damaged battery pod acts on the upper, port side lug, while in the second case, the same force acts on the lower, starboard lifting lug. In both cases, the external pressure of 0.65 MPa on the test depth is considered as well. Boundary conditions for this lifting condition are simulated using spring to ground elements with very low stiffness of 1 N/mm to prevent rigid body motion. The loads when lifting force acts on the lower starboard lug are presented in Figure 16.
Results from the nonlinear analysis are presented in Figure 17 where mid-plate von Mises stresses are shown when the lifting force acts on the lower starboard lug. The maximum stress of 174 MPa is in same range and below permissible values as the one obtained in the global strength analysis for the test depth presented in Figure 7. Moreover, von Mises stresses around the loaded lifting lug are well below permissible values. Comparable results are obtained in the case when the lifting force acts on the upper portside lug as well.
The third analysed lifting condition is the launching/recovering condition at rough seas where the primary lifting lugs and lifting attachments must be analysed for forces of at least 2 g vertical (1 g static plus 1 g dynamic), 1 g transverse, and 1 g longitudinal, acting simultaneously under the most severe loading condition [8]. Thus, the resultant acceleration for the emergency lifting condition reads 25.03 m/s2.
The results from the nonlinear analysis are presented in Figure 18 and Figure 19 where a deformed model with von Mises stresses is shown. The maximum mid-plate von Mises stress of 243 MPa is recorded as may be seen in Figure 18a. Based on these results, designers reinforced the critical areas with high-tensile steel with permissible values of 332 MPa (80% of yield strength). Figure 18b shows von Mises stresses for lifting lugs, which are modelled by solid elements and high-tensile steel. Based on the ABS Rules [8] for the allowable stress of contact problem, 100% of the yield strength of the material (415 MPa) can be used as the limit value. Thus, according to the results in Figure 18b, all solid elements satisfy the given criteria. Furthermore, the total translation can be seen in Figure 19, where it can be concluded that the model rotates around the x axis.

4. Discussion about Non-Metallic Parts of the Pressure Hull

There are two major non-metallic structural groups: the acrylic cylinders and parts of the exo-structure such as soft ballast tanks. The latter parts are not subjected to differential pressures and thus are not required to be included in the pressure hull design. Acrylic cylinders, however, are designed as per the PVHO-1 design standard which does not allow for direct structural assessment of windows [27]. As the stiffness of the acrylic material is almost negligible compared with the steel, the influence of the acrylic cylinders on the structural strength of the pressure hull is rather small. For that reason, acrylic cylinders can be omitted from the analysis of the steel pressure hull without affecting the accuracy of the analysis. Since this paper deals with the analysis of metallic parts of the submersible’s structure, only basic aspects of design and analysis of polymeric parts will be given here. Due to the importance of acrylic cylinders, these aspects must be considered even when the designer is concerned only with the metallic part of the structure and it can be assessed by means of simple formulae and diagrams, at least as a preliminary assessment.
Acrylic cylindrical hull sections are fitted between supporting rings (Figure 2). Due to the structural arrangement where the exo-structure with the fore and aft ends takes most of the loads, the acrylic cylinders are loaded only locally, i.e., external hydrostatic pressure that acts on them is transmitted to the steel rings which are rigidly connected with the exo-structure. Thus, the pressure causes the cylinders to bend between the rings. Since the cylinders are not free to expand in the longitudinal direction, this will typically cause axial stress [2].
The analysis of the problem can be simplified by observing that cylinders, the loads, and the supports are axisymmetric. This means that the problem can be analysed as an axisymmetric problem instead of as a 3D problem. Not only does this facilitate a potential thorough analysis, but it also enables the use of simple analytical formulae that can be found in many engineering handbooks [33].
Besides bending, the cylinders will exhibit significant thermal deformations too. This is a consequence of the thermal properties of acrylic plastic. For example, the coefficient of linear thermal expansion will increase as much as 16% from 4 °C to 27 °C based on data from [27]. This will lead to a variation in the length of the cylinders of a couple of millimetres. If unaccounted for, significant axial forces will arise, which can lead to fatigue failure [2].
Fatigue can also arise due to axial force variation [2]. In the submersible design proposed and analysed in this paper, the exo-structure ideally takes all global loads, which also includes axial force that arises due to hydrostatic pressure that acts on the fore and aft ends of the submarine. Since the ends are fixed to the exo-structure and very stiff, maximum horizontal deformation at the connection with the cylinders is 0.53 mm, which can be obtained from the structural analysis demonstrated in Section 3.2.1 and adjusted for the operational depth of 50 m. This is very close to the thermal deformation obtained in [9] and should be accounted for by the deformation of the side gasket. Structural analysis of a joint between the steel ring and the cylinder is a complex problem that will be dealt with in a separate paper.
The third potential failure mode of the acrylic cylinders is buckling [2]. The cylinders under external pressure are much less susceptible to buckling failure than cylinders under internal pressure [2]. To assess the critical external pressure that causes buckling, Stachiw proposes a set of diagrams giving the critical pressure as a function of thickness to diameter and length to diameter ratios [2]. Based on the diagrams, the critical pressure for a cylinder with dimensions specified in Table 2 is well above the hydrostatic pressure at the design depth.
The second non-metallic critical part, i.e., the joint between the cylinders and the steel rings, as shown in Figure 3, is composed of a seal and three gaskets. Following a joint design proposed in [2], the seal is made of nitrile butadiene rubber (NBR) and the gaskets are made of chloroprene rubber (CR). Both materials are elastomers, polymeric materials that can undergo large elastic deformations. Their behaviour is hyperelastic and viscoelastic [28] and can be analysed only by means of numerical methods that can account for large deformations, finite strains, and geometrical and material nonlinearity.
Although the material behaviour is, in general, non-linear elastic, both materials show an almost linear elastic behaviour up to an engineering strain of about 0.2 for CR and about 0.3 for NBR for a strain state relevant in this analysis [29] with a modulus of elasticity of about 6.5 MPa. This approximation can be used for simple, preliminary assessments of their structural behaviour as demonstrated in a simple calculation of tensile force in the tie rods. Limitations and underlying assumptions of such an analysis are given in Section 2.2 and should always be observed. Additionally, the seal joint between the cylinders and the steel rings are planned to be experimentally tested to obtain material multi-axial mechanical properties.

5. Conclusions

This paper describes and demonstrates the procedure based on the ASME design-by-analysis approach for the structural design of metallic pressure hull elements of the tourist submersible. Design-by-analysis is a recently introduced approach in structural design practice for pressure vessels, while the application for underwater vehicles described in the present paper is a novelty and main contribution of the study. Based on the proposed approach, a novel touristic submersible, with a capacity of 48 passengers at nominal diving depth of 50 m, was analysed. Furthermore, insight into the analysis of non-metallic parts of the pressure hull as well as preliminary pretension force calculation were given. As the stiffness of the non-metallic parts of the submarine is almost negligible compared with the steel, the influence of the acrylic cylinders on the structural strength of the pressure hull is not considered.
According to the design-by-analysis approach, the pressure vessel component should be evaluated for different failure modes, i.e., plastic collapse, local failure, buckling, and cyclic loading. These failure modes were investigated for the individual components and for the structural assembly as well. Radial and angular deformation analyses were performed for connections of acrylic cylinders according to safety standards for pressure vessels for human occupancy. Finally, important lifting conditions were investigated to optimize the location and size of the lifting lugs. As the lifting analysis revealed the location of high-stress concentration, the lifting ears were reinforced using high-tensile steel. All calculations for operational and design depth were performed using pretension force in tensioned rods that enable safe fastening of the acrylic cylinders with thermal variations and diving depth.
The results of the finite element method revealed that protection against plastic collapse of the pressure hull design is satisfactory, as the primary membrane von Mises stresses were always below 150 MPa, while local membrane stresses satisfy the criteria of 225 MPa. Buckling factors for all relevant metallic structures, i.e., transverse web ring, battery cylinder pod, and fore and aft cupola were above minimum design values. Next, the fatigue life for operational depth and 25,000 cycles was found to be adequate according to the relevant S-N curve where the connection between the aft/fore head and aft/fore cupola ring were analysed as the most dangerous locations for potential fatigue failure. Finally, it was shown that the lifting operations can be performed safely since stresses and deformation for various lifting cases are at satisfactory low levels after reinforcing the lifting ears.
As demonstrated in the study, the design-by-analysis procedure can be recommended for the design of metallic parts of innovative underwater vehicles and other structures exposed predominantly to the pressure loading.

Author Contributions

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

Funding

European Union from the European Regional Development Fund within the Operational Program “Competitiveness and Cohesion 2014-202”, project KK.01.2.1.02.0339—Development of the multipurpose luxury touristic and research submarine.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The design of the submarine that is the subject of this study was conducted by Marine and Energy Solutions DIV d.o.o, Zagreb, Croatia. The project is co-financed by the European Union from the European Regional Development Fund within the Operational Program “Competitiveness and Cohesion 2014–2020”, project KK.01.2.1.02.0339—Development of the multipurpose luxury touristic and research submarine. The content of the publication is the sole responsibility of the project partner Faculty of Mechanical Engineering and Naval Architecture, University of Zagreb, Croatia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. How Far Does Light Travel in the Ocean? Available online: https://oceanservice.noaa.gov/facts/light_travel.html (accessed on 26 June 2022).
  2. Stachiw, J.D. Handbook of Acrylics for Submersibles, Hyperbaric Chambers and Aquaria; Best Publishing Company: North Plam Beach, FL, USA, 2003. [Google Scholar]
  3. Du, Q.; Hu, Y.; Cui, W. Safety assessment of the acrylic conical frustum viewport structure for a deep-sea manned submersible. Ships Offshore Struct. 2017, 12(S1), 221–229. [Google Scholar] [CrossRef]
  4. Pranesh, S.B.; Kumar, D.; Subramanian, V.A.; Sathianarayanan, D.; Ramadass, G.A. Numerical and experimental study on the safety of viewport window in a deep sea manned submersible. Ships Offshore Struct. 2020, 15, 769–779. [Google Scholar] [CrossRef]
  5. Pope, G.T. Touring the Deep. Pop. Mech. 1988, 10, 67–70, 102. [Google Scholar]
  6. Forman, W.R. Submersible Observation Vessel. U.S. Patent US4928614, 29 May 1990. [Google Scholar]
  7. Kohnen, W. Manned Underwater Vehicles 2017–2018 Global Industry Overview; Marine Technology Society (MTS): Washington, DC, USA, 2018. [Google Scholar]
  8. Rules for Building and Classing Underwater Vehicles, Systems and Hyperbaric Facilities; ABS: Spring, TX, USA, 2021.
  9. Xing, Y.; Santoso, T.A.D.; Ma, Y. Technical-Economic Feasibility Analysis of Subsea Shuttle Tanker. J. Mar. Sci. Eng. 2022, 10, 20. [Google Scholar] [CrossRef]
  10. Hughes, O.W.; Paik, J.K. Ship Structural Analysis and Design; SNAME: New York, NY, USA, 2010. [Google Scholar]
  11. Payer, H.G.; Schellin, T.E. A class society’s view on rationally based ship structural design. Ships Offshore Struc. 2013, 8, 319–336. [Google Scholar] [CrossRef]
  12. ASME Boiler & Pressure Vessel Code, Section VIII, Division 2; ASME: New York, NY, USA, 2019.
  13. The Design-by-Analysis Manual; European Commission: Petten, The Netherlands, 1999.
  14. Diamantoudis, A.T.; Kermanidis, T. Design by analysis versus design by formula of high strength pressure vessels: A comparative study. Int. J. Press. Vessel. Pip. 2005, 82, 43–50. [Google Scholar] [CrossRef]
  15. Karthikeyan, K.M.B.; Balasubramanian, T.; Bruce, R.A.; Premkumar, P. Pressure Vessel Design by Design by Analysis Route. In Proceedings of the IOP Conference Series: Materials Science and Engineering, Chennai, India, 20 February 2020. [Google Scholar]
  16. Cho, S.-R.; Muttaqie, T.; Lee, S.H.; Paek, J.; Sohn, J.M. Ultimate Strength Assessment of Steel-Welded Hemispheres under External Hydrostatic Pressure. J. Mar. Sci. Appl. 2020, 19, 615–633. [Google Scholar] [CrossRef]
  17. Cho, S.-R.; Muttaqie, T.; Do, Q.T.; Kim, S.; Kim, S.M.; Han, D.-H. Experimental investigations on the failure modes of ring-stiffened cylinders under external hydrostatic pressure. Int. J. Nav. Archit. 2018, 10, 711–729. [Google Scholar] [CrossRef]
  18. Pranesh, S.B.; Sathianarayanan, D.; Ramadass, G.A. Design standards for steel spherical pressure hull for a manned submersible. J. Ocean Eng. Mar. Energy 2022, 8, 137–151. [Google Scholar] [CrossRef]
  19. Ma, Y.; Xing, Y.; Ong, M.C.; Hemmingsen, T.H. Baseline design of a subsea shuttle tanker system for liquid carbon dioxide transportation. Ocean Eng. 2021, 240, 109891. [Google Scholar] [CrossRef]
  20. Sohn, J.M.; Hirdaris, S.; Romanoff, J.; Kim, S.J. Development of Numerical Modelling Techniques for Composite Cylindrical Structures under External Pressure. J. Mar. Sci. Eng. 2022, 10, 466. [Google Scholar] [CrossRef]
  21. Zhou, M.; Patel, A.; Wang, B.; Jin, W. Design Optimization of Pressure Vessel in Compliance With Elastic Stress Analysis Criteria for Plastic Collapse Using an Integrated Approach. J. Press. Vessel. Technol. 2021, 143, 011303. [Google Scholar] [CrossRef]
  22. Rules for Conditions of Classification; ABS: Spring, TX, USA, 2022.
  23. Prebeg, P.; Palaversa, M.; Andric, J.; Tomicic, M. Adaptation of FEM-based open-source software for ship structural analysis. Ships Offshore Struct. 2022, 1–12. [Google Scholar] [CrossRef]
  24. Duan, C.; Ding, L.; Lu, M. Discussion on the Implementation of the Primary Structure Method in Design by Analysis. Appl. Mech. Mater. 2017, 853, 341–345. [Google Scholar] [CrossRef]
  25. Saghi, H.; Parunov, A.; Mikulic, A. Resistance Coefficient Estimation for a Submarine’s Bare Hull Moving in Forward and Transverse Directions. Appl. Sci. 2022, 12, 10953. [Google Scholar] [CrossRef]
  26. Gatin, I.; Čokić, J.; Romić, D.; Parunov, J. CFD Study on the Influence of Exostructure Elements on the Resistance of a Submarine. J. Mar. Sci. Eng. 2022, 10, 1542. [Google Scholar] [CrossRef]
  27. ASME PVHO-1 Safety Standard for Pressure Vessels for Human Occupancy; ASME: New York, NY, USA, 2019.
  28. Bergstrom, J. Mechanics of Solid Polymers: Theory and Computational Modeling; William Andrew-Elsevier: San Diego, CA, USA, 2015. [Google Scholar]
  29. Palaversa, M.; Parunov, J. FEA of Hyperelastic Structures: A Case from the Submarine Design. In Proceedings of the Sorta 2022, Malinska, Croatia, 7 September 2022. [Google Scholar]
  30. Kim, B.; Lee, S.B.; Lee, J.; Cho, S.; Park, H.; Yeom, S.; Park, S.H. A Comparison Among Neo-Hookean Model, Mooney-Rivlin Model and Ogden Model for Chloroprene Rubber. Int. J. Precis. Eng. Manuf. 2012, 13, 759–764. [Google Scholar] [CrossRef]
  31. Det Norke Veritaas. Fatigue Assessment of Ship Structures, Classification Notes No. 30.7; Det Norke Veritaas: Høvik, Norway, June 2010. [Google Scholar]
  32. Ćorak, M.; Šperanda, Z.; Čokić, J.; Parunov, J. Structural analysis of tourist submarine with acrylic hull. In Sustainable Development and Innovations in Marine Technologies, 1st ed.; Ergin, S., Soares, C.G., Eds.; CRC Press: London, UK, 2022. [Google Scholar] [CrossRef]
  33. Young, W.C.; Budynas, R.G. Roark’s Formulas for Stress and Strain, 7th ed.; McGraw-Hill: New York, NY, USA, 2002. [Google Scholar]
Figure 1. The standard and the design-by-analysis procedure.
Figure 1. The standard and the design-by-analysis procedure.
Jmse 11 00275 g001
Figure 2. General view of the transparent acrylic submersible hull.
Figure 2. General view of the transparent acrylic submersible hull.
Jmse 11 00275 g002
Figure 3. Scheme of a joint between the acrylic cylinder and the supporting steel ring (dimensions in mm).
Figure 3. Scheme of a joint between the acrylic cylinder and the supporting steel ring (dimensions in mm).
Jmse 11 00275 g003
Figure 4. Loads on the transverse ring (pressure in MPa).
Figure 4. Loads on the transverse ring (pressure in MPa).
Jmse 11 00275 g004
Figure 5. Buckling mode of the battery cylinder pod.
Figure 5. Buckling mode of the battery cylinder pod.
Jmse 11 00275 g005
Figure 6. Loads and boundary conditions on the submarine structural assembly (force in N, pressure in MPa).
Figure 6. Loads and boundary conditions on the submarine structural assembly (force in N, pressure in MPa).
Jmse 11 00275 g006
Figure 7. Yielding results of the structural assembly (stresses in MPa).
Figure 7. Yielding results of the structural assembly (stresses in MPa).
Jmse 11 00275 g007
Figure 8. Deformation (in mm): (a) y direction (T2), (b) z direction (T3).
Figure 8. Deformation (in mm): (a) y direction (T2), (b) z direction (T3).
Jmse 11 00275 g008
Figure 9. Buckling mode for structural assembly.
Figure 9. Buckling mode for structural assembly.
Jmse 11 00275 g009
Figure 10. Loads and boundary conditions for fatigue analysis (pressure in MPa).
Figure 10. Loads and boundary conditions for fatigue analysis (pressure in MPa).
Jmse 11 00275 g010
Figure 11. Plate bottom von Mises stresses for design depth and detail of the connection of the aft head and the first ring (stresses in MPa).
Figure 11. Plate bottom von Mises stresses for design depth and detail of the connection of the aft head and the first ring (stresses in MPa).
Jmse 11 00275 g011
Figure 12. Fatigue life of the analysed detail [12].
Figure 12. Fatigue life of the analysed detail [12].
Jmse 11 00275 g012
Figure 13. Boundary conditions and gravity loads for emergency lifting condition.
Figure 13. Boundary conditions and gravity loads for emergency lifting condition.
Jmse 11 00275 g013
Figure 14. The von Mises stresses in MPa of the vessel for the lifting condition; (a) plate elements, (b) solid elements.
Figure 14. The von Mises stresses in MPa of the vessel for the lifting condition; (a) plate elements, (b) solid elements.
Jmse 11 00275 g014
Figure 15. Vertical displacements (in mm) of the assembled structure for the lifting condition.
Figure 15. Vertical displacements (in mm) of the assembled structure for the lifting condition.
Jmse 11 00275 g015
Figure 16. Loads and vertical displacements of the assembled structure for the emergency lifting in damaged condition (force in N, pressure in MPa).
Figure 16. Loads and vertical displacements of the assembled structure for the emergency lifting in damaged condition (force in N, pressure in MPa).
Jmse 11 00275 g016
Figure 17. Mid-plate von Mises stresses (in MPa) for emergency lifting in the damaged condition.
Figure 17. Mid-plate von Mises stresses (in MPa) for emergency lifting in the damaged condition.
Jmse 11 00275 g017
Figure 18. High-stress area in lugs for launching/recovering in rough seas: (a) plate elements with mid-plate von Mises stress above 192 MPa; and (b) solid elements with von Mises stresses above 230 MPa.
Figure 18. High-stress area in lugs for launching/recovering in rough seas: (a) plate elements with mid-plate von Mises stress above 192 MPa; and (b) solid elements with von Mises stresses above 230 MPa.
Jmse 11 00275 g018
Figure 19. Mid-plate von Mises stresses for launching/recovering in rough seas.
Figure 19. Mid-plate von Mises stresses for launching/recovering in rough seas.
Jmse 11 00275 g019
Table 1. The particulars of the submersible.
Table 1. The particulars of the submersible.
FeatureValue
Test depth65 m
Length overall25.09 m
Beam4.750 m
Draft3.420 m
Dry weight155 tons
Hull acrylic outer diameter2.64 m
Forward speed (max)2.5 knots
Passengers48
Table 2. Dimensions of a cylindrical hull section and the side gasket.
Table 2. Dimensions of a cylindrical hull section and the side gasket.
DimensionValue
Length of the cylinder (L)2.4 m
Inner diameter of the cylinder (d)2.34 m
Thickness of the cylinder (h)140 mm
Thickness of the side gasket (t)10 mm
Table 3. Structural component results.
Table 3. Structural component results.
ComponentVon Mises Membrane Stresses (MPa)Buckling Factor
Battery cylinder pod138.5
Fore cupola16431.5
Aft cupola16331.5
Transverse ring1064.5
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

Ćorak, M.; Šperanda, Z.; Čokić, J.; Palaversa, M.; Parunov, J. Finite Element Simulations of Novel Submersible Design Based on the ASME Design-by-Analysis Approach. J. Mar. Sci. Eng. 2023, 11, 275. https://doi.org/10.3390/jmse11020275

AMA Style

Ćorak M, Šperanda Z, Čokić J, Palaversa M, Parunov J. Finite Element Simulations of Novel Submersible Design Based on the ASME Design-by-Analysis Approach. Journal of Marine Science and Engineering. 2023; 11(2):275. https://doi.org/10.3390/jmse11020275

Chicago/Turabian Style

Ćorak, Maro, Zdenko Šperanda, Juvel Čokić, Marin Palaversa, and Joško Parunov. 2023. "Finite Element Simulations of Novel Submersible Design Based on the ASME Design-by-Analysis Approach" Journal of Marine Science and Engineering 11, no. 2: 275. https://doi.org/10.3390/jmse11020275

APA Style

Ćorak, M., Šperanda, Z., Čokić, J., Palaversa, M., & Parunov, J. (2023). Finite Element Simulations of Novel Submersible Design Based on the ASME Design-by-Analysis Approach. Journal of Marine Science and Engineering, 11(2), 275. https://doi.org/10.3390/jmse11020275

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