Pressure Hull Design Methods for Unmanned Underwater Vehicles

: This paper describes design methods for the plastic hull of an Unmanned Underwater Vehicle (UUV), with a particular focus on its cylindrical body and nearly spherical domes at the ends. With the proposed approach, the methodologies reported in the literature were compared, and suitable modiﬁcations and improvements were investigated and implemented to extend the classical theories and data to this case study. The investigated underwater vehicle, named FeelHippo, was designed and assembled by the Department of Industrial Engineering of the University of Florence. Its main hull is composed of an extruded PMMA (PolyMethyl MethAcrylate) cylinder and two thermoformed PMMA domes. Breakage of the hull results in destructive phenomena, namely, yielding and buckling. An experimental campaign and FEM (Finite Element Method) analysis were carried out to complete the theoretical study, and the collapse pressures were compared with the derived design values. In conclusion, the proposed innovative method is a lean and effective technique for designing underwater hull domes and predicting the collapse pressures.


Introduction
This paper compares collapse pressure prediction techniques in the literature and presents an innovative hull design strategy starting from an Autonomous Underwater Vehicle (AUV) developed by the Mechatronics and Dynamic Modelling Laboratory (MDM Lab) of the University of Florence. The study is based on theoretical checks, Finite Element Method (FEM) analysis, and experimental tests, and it aims to use the available data to extend the scope of the relationships established in the literature.
In the 1950s [1], the first ROV (Remotely Operated Vehicle) was developed for archaeological research, namely, the POODLE, and for military uses, namely, the Cable-Controlled Underwater Research Vehicle (CURV). Since then, there has been growing interest in the underwater environment and the development of Unmanned Underwater Vehicles (UUVs). The design of underwater vehicles is highly dependent on the sector that will implement it (e.g., [2][3][4]). The shape of the hull, a primary focus of the design, is most often cylindrical [5] or spherical. In underwater robotics [6], the more commonly used geometry is spherical rather than a flat surface because the distribution of stress on the surface is homogeneous, and the hydrodynamic drag is reduced. The hull material is another fundamental aspect to consider since the production process can introduce non-negligible peculiarities [7,8]. In this case, the use of PolyMethyl MethAcrylate (PMMA) as a material allows for an appropriate balance between Yielding is axisymmetric inelastic, while local instability involves an asymmetrical inelastic collapse. The current trend of designing hulls of underwater vehicles is to consider a specific sizing so that only yielding can occur at the design pressure. Because appropriate safety coefficients are used, the critical factors linked to both local and general instability should appear after the operating pressure exceeds the design pressure. In particular, the failure modes due to global instability are very sensitive to manufacturing imperfections [20].
When the external pressure generates compressive stresses that reach high levels, instability buckling phenomena predominate. For basic geometries, the linear buckling load can be determined by analytical and numerical studies [21]. However, the linear buckling load is only an indication of the real buckling resistance of a shell; other important factors, such as geometric imperfections and the effects of the boundary conditions, must be considered [19] to obtain the actual buckling strength. For this reason, this study focused on breakage due to the instability of cylinders. Sizing is carried out to ensure that the structure will not be damaged, regardless of the phenomenon that occurs.
With the application of a uniform pressure range, the collapse of a circular cylinder occurs with an external pressure that is a small fraction compared with the internally applied pressure. This mode of failure is known as shell instability or lobar buckling, and it causes a collapse around its circumference in the form of a number of circumferential waves or lobes [22].
General instability manifests as the physical collapse of the entire shell because of the low strength of the ring-stiffeners. Finally, axisymmetric deformation takes place if the circular cylinder implodes while maintaining its circular shape [22]. Resistance to external pressure is further diminished by initial out-of-circularity [23].
The second geometric shape studied is the dome, which can collapse in different ways but presents a sudden loss of load capacity triggered by buckling. Studies on nonlinearity and marked imperfection sensitivity have been carried out since the study in [24]. Under a uniformly compressed state, both axisymmetric and asymmetrical behaviors can occur, as shown in [25]; in particular, the study focused on buckling and compression. These types of collapse are based on geometric and structural characteristics: for this reason, various breaking prediction theories have been studied to prevent these phenomena, as is the case in the present study. For this geometry, this study builds on the literature and represents a step forward in collapse prediction theories.
A hemispherical shell is able to withstand higher internal pressure than any other geometrical hull with the same wall thickness and radius. Thus, it is a major component of pressure vessel construction and has a smaller surface area per unit volume than any other shape. For complicated structures, such as shells, the modern design technique tends to improve the model investigation since, in most cases, the true behavior of the shell, with the load and the constraints, is not known or very difficult to know. The best approach is to make some assumptions and then verify them by using dedicated tests [26].
Although the literature reports numerous complex theoretical treatments, there is a lack of implementation in real and diverse applications [21], and there are few validations of the proposed theories. The cases presented in the state-of-the-art literature do not fit the results obtained herein (Section 5) because of the non-negligible discrepancy in the constraints and the production process (e.g., [7]). A similar problem, which was found to be inherent in the difficult task of associating a coefficient with a type of constraint, was discussed in [27].
Both the cylinder and domes investigated in this study are made of PMMA, so the objective of this work is also to adapt the classical methods to a different material and, above all, the production method. The most common materials for these types of mobile robots are metal alloys (high-strength steels, aluminum, or titanium) or composites (glass-fiber-reinforced plastic, carbon epoxy, or metal matrix composite) [4,28]. Furthermore, compared with the most common AUVs (e.g., REMUS-100 AUV [29], Fòlaga AUV [30], SPARUS II AUV [31], and Matsya 5.0 [3]), the shape of the designed vehicle investigated in this study is more compact and less torpedo-like.
This study focused on the autonomous vehicle called FeelHippo AUV [32] (shown in Figure 2, Section 2). First, the theoretical critical buckling pressure of the cylinder was calculated using classical formulas [33]; second, the critical buckling pressure was analytically derived and used to determine the optimized size in order to calculate the lowest thickness of the vehicle's central cylinder to address the buckling effect. The available domes were experimentally tested until breakage, and then the dome design pressures were obtained. As described below, additional correlations that were introduced resulted in good agreement between the theoretical results and both the FEM and experimental results.
An overall operational water depth of 30 m was determined. Finally, from the achieved results, a new restyled version of the vehicle was realized.
For the cylindrical geometry, both a theoretical comparison and a simulated validation were performed. The theoretical formulations in the literature were used for verification. The widely applied cylindrical geometry was used to approach the problem of collapse and allowed the research group to become familiar with the phenomenon. For the domes, first, an experimental campaign was carried out. Then, the theories available in the literature were studied: because the correspondence with the proposed case was not satisfactory, some improvements and additions were applied. The results obtained with the new proposed approach are encouraging and allow for the extension of the cases investigated in the literature. Furthermore, the dome geometry is widely used in the field of submarine robotics.
This paper presents the redesign of FeelHippo AUV, particularly its central hull-a cylinder and two flanged domes. The remainder of the paper is organized as follows: Section 2 illustrates the main mechanical and electronic characteristics of the vehicle. Section 3 focuses on the design of the cylindrical component by starting from theory and the state-of-the-art literature and ending with the investigated case and related validation. In Section 4, experimental tests on the dome are described, and Section 5 expands on the design of the domes by applying theoretical treatments and improvements to match the experimental evidence. Finally, Section 6 presents the conclusions and proposes some future developments.
The overall aim of this research is to create a simplified calculation system for the design and sizing of submarine hull parts under pressure and, in particular, to codify a fast and light computational procedure to check the resistance of cylinders and domes. A dedicated experimental campaign enabled the proposal of an extension of classical test cases in the literature, especially for the second geometry. The experimental campaign, in which the collapse of available domes was tested, facilitated the creation of a small dataset from which the theoretical equations were extended and validated. The main contribution of this research is the proposed design method to identify the mechanical resistance of the domes. For the sake of clarity, the classical theoretical theories and thus the FEM were based on the assumption of a perfect sphere, with only the stiffness of the material taken into account. In this work, the thickness variation and the flange constraints were introduced to consider the real shape of the dome and extrapolate suitable corrective coefficients. The proposed approach starts from the geometrical data and enables the identification of the failure mode and the reachable operational water depth.

Case Study: FeelHippo AUV
FeelHippo AUV, shown in Figure 2, is an AUV specially designed to be a development platform to test new payload and new control algorithms for both the MDM Lab research group and its students. This is reflected by the vehicle's involvement in several monitoring tasks and student competitions: SAUC-E (Student Autonomous Underwater Challenge-Europe) in both 2013 and 2016; Eurathlon in 2015; and ERL Emergency Robots in 2017 [34], 2018, and 2019. For this reason, the design specifications for the vehicle are very simple: light, cheap, small, and with an operating depth of a few meters. In the study reported in this paper, the real depth of immersion and navigation in safety were calculated and found to correspond to about 30 m.
The layout of the thrusters imparts FeelHippo AUV with an underactuated motion, with four marine propellers [35] controlling 4 degrees of freedom (DOFs) of the vehicle: two of these are mounted laterally in a V-shape to control the depth and the lateral movement (which is necessary for precise hovering); the other two are mounted in the backside on the horizontal plane, one per side, to control the motion along the longitudinal advancing direction and the yaw angle [36]. The main structure of the vehicle is made of Anticorodal type 6082-T6 aluminum alloy and PMMA, commercially known as Plexiglas R (mechanical characteristics in Tables 1 and 2). In Figure 3, a CAD of the vehicle is shown to highlight the overall dimensions. More details are shown in Figure 4. The vehicle's central body is a plexiglass tube containing the instrumental hardware, and Anticorodal type 6082-T6 constitutes the structural connections, and the two tubes under the main body house the batteries. The body consists of an extruded cylinder and two flanged thermoformed domes. The thermoforming process entails heating a slab of material to 150-160 • C and then deforming it under vacuum by blowing without a counter-mold to maintain good optical quality. As discussed in Section 5, this implies that the thickness is variable and that the shape is not perfectly spherical.
The choice of materials was dictated by compromising between lightness, strength, and the ability to visually check the electronic components. Other external elements were printed in ABS (Acrylonitrile Butadiene Styrene) with a 3D printer from the MDM Lab (Stratasys Dimension Elite).
The differences between the selected materials had a profound influence on the design phase. In fact, the aluminum parts were sized to ensure structural resistance to impacts, which represent another type of possible failure. This choice also ensured that the aluminum parts could endure buckling instability: for this reason, given the difference in the thicknesses of the various elements, this work performed buckling optimization only for the PMMA components (also because this material has a much lower yield strength and Young's modulus). The PMMA domes and the main tube are connected by aluminum rods and flanges, which are equivalent to ring-stiffeners ( Figure 4).
The FeelHippo design is thought to have limited roll and pitch angles during vehicle navigation. This feature is commonly adopted in AUV designs and is obtained by placing its center of gravity below its center of buoyancy. In this case, the AUV has a large and light main body, which contains the mainboards and other hardware, and two smaller and heavy legs, which contain the batteries and balancing weights.
In the electronic architecture, the main computer communicates with other devices by means of serial communication and an internal Ethernet network. Devices connected with serial communication are

•
Radio modem to send short messages from high distances on the surface; • GPS receiver on the surface to quickly localize the vehicle; • FOG (Fiber Optic Gyroscope), a single-axis gyroscope with high precision to improve the pose estimation (particularly the yaw measurement); • Custom board to monitor internal temperature; it also includes a water leakage sensor, battery SOC, and a control custom beacon flash and status LEDs; • USB camera to provide frontal images during navigation; • IMU (Inertial Measurement Unit) to provide data from the 3D internal accelerometer, gyroscope, and magnetometer; • Acoustic modem, which enables a high data transfer rate (13.9 kbit/s) and a high functioning range (3500 m); • Servo driver, an electronic board to control the four outrunner brushless motors.
The other devices, which are connected to each other in a cabled Ethernet network with two switches, are • DVL (Doppler Velocity Log), which provides a 3-axial linear speed and altitude measure referred to the seabed using the Doppler effect and the depth of the vehicle by means of an integrated pressure sensor; • WiFi AccessPoint, which allows for fast, high-band, short-range communication on the surface; • 2D forward-looking sonar that allows the vehicle to see obstacles in situations of poor visibility; • IP cameras to acquire images of the seabed and facilitate 3D reconstruction.
The energy supply that powers the vehicle comprises three LiPo batteries placed inside the aluminum pipes. In addition, other voltage levels are provided by dedicated DC/DC converters. The main voltage value lines are shown in Figure 5, along with the brands and models of the devices.
All these devices are largely in the main plexiglass body, stacked on two plastic planes. Only the IP cameras, 2D forward-looking sonar, and DVL are hung from the external structure.

Design of the Main Cylinder of the AUV Hull
In this section, the equations for local instability are given. The critical buckling pressure equations were derived from the buckling analysis proposed both by R. Von Mises and by D.F. Windenburg and C. Trilling. Particular emphasis was placed on the case of submarines that have a circular cylinder, are closed at the ends, and are deeply submerged in water and thus exposed to pressure from all sides [37]. Since thin cylinders subjected to uniform external pressure collapse in an asymmetrical mode (lobar buckling), the local instability was studied at a fraction of the pressure needed to cause axisymmetric yielding.
Specifically, an analytical reference formulation for the study of short thin tubes that are supported at the ends and subjected to uniform radial and axial pressure was proposed. The formulation was subsequently modified by Von Mises [22]: (1) ν = Poisson's ratio; n = number of circumferential waves lobes (an integer that minimizes the expression).
In a precautionary manner, the difference between the real scenario and the theoretical hypotheses was considered by introducing the corrective factor k [20].
where k varies according to t: k = 0.4 for t < 5 mm k = 0.5 for 5 < t < 7 mm k = 0.6 for t > 7 mm P design = operative pressure at which the vehicle can dive and navigate safely.
The following expression, proposed by D.F. Windenburg and C. Trilling, does not depend on the number of lobes [23]: Its result differs by a maximum of 3.5% from the result of Von Mises's formulation, so it is sufficiently accurate for practical applications.
Ross presented a comparison between the theoretical models, the experimental results, and simulated results [21]. For instance, Figure 6 reports the results for an Anticorodal 6082-T6 tube with the outer diameter 2R + t = 50 mm and t = 1.7 Figure 6. Predicted buckling pressure versus tube length [23].
The literature also includes graphs to calculate the pressure that is actually bearable by the cylinder according to its geometric characteristics and its material. The calculations are based on both a theoretical basis and experimental results; therefore, the results are even more reliable. Figure 7 presents a design chart for thick-walled circular cylinders that collapse under external hydrostatic pressure [22].  In addition, a safety coefficient (SF) < 3 e.g., arbitrarily 1.5 is added to obtain the real design operative pressure (a value of 3 is used for the case in which people are onboard). It can be seen that if 1/λ c < 0.9, the structure will probably fail by elastic instability, followed by axisymmetric deformation. Additionally, there is a link between plastic axisymmetric deformation and inelastic shell instability.
The ability to withstand pressure is reduced if the cylinder is manufactured with an initial out-of-circularity. Comparing experimental and theoretical data reveals that as the length of the cylinder decreases, the theoretical data deviate from the experimental evidence ( Figure 6) because of edge effects, too. Therefore, corrective factors are needed to ensure a certain safety margin.
For the geometry of FeelHippo AUV, in accordance with the aforementioned nomenclature, the cylinder data are summarized in Table 3. First, the theoretical critical buckling pressure (the pressure at which the hull would collapse under ideal conditions) was calculated. Table 4 summarizes the theoretical pressure for the case under examination: the first two columns identify the cylindrical thickness and the number of lobes; the third column reports the values from Equation (1), the fourth column reports the values from Equation (2), the fifth column reports the values from Equation (4), and the last column reports the results of the FEM analysis described below. The values are quite similar among the various formulations. The necessary geometry was obtained by carrying out sizing optimization to perform an FEM comparison and an evaluation of the minimum volume of material-that is, the minimum thickness of the cylinder to withstand the calculated pressure, less than 0.7 MPa. As a result, the axial load was imposed on the frontal section while a uniform external pressure was applied to the structure. The first thickness was set to 10 mm to avoid influencing the test.
The FEM results of the sizing optimization analysis by HyperMesh are shown in Figure 8 (static pressure) and Figure 9 (qualitative, first buckling mode). After several iterations, the solver identified a minimum necessary thickness of 5.08 mm, which can be reasonably approximated to 5 mm.    There is a complete match between the results of the FEM analysis and the theoretical results derived from Von Mises's theory. The number of circumferential lobes is 4, which is in agreement with both Equations (1) and (2).
The last step includes the derivation of the real critical buckling pressure by using both the classical and experimental methods. In sum, the first method uses the theoretical value of the pressure derived by R. Von Mises' or D.F. Windenburg and C. Trilling's equations, with corrections applied by using the corrective factor k. The second (experimental) method corrects the theoretical value through experimental charts (λ c and PKD) and the safety coefficient SF. From the chart in Figure 7, it is possible to calculate the coefficient λ c , particularly λ c = 2.87 (Equation (5)), to find the value of PKD 1.3 and, consequently, P exp = 0.49 MPa (Equation (6)), which is the experimental bearable pressure value. Using a safety factor (Equation (7)) SF = 1.5, the design pressure P design = 0.32 MPa (about 30 m) is obtained. This is the same value that results from using the safety coefficient k = 0.5 and P crit of Equation (1) in Equation (3).
In conclusion, the operative depth obtained for the case under examination is about 30 m. A theoretical investigation was also carried out for other cases that have been tested by other researchers [38]; this renders the validation more credible if the same results are obtained. Table 5 summarizes the buckling pressure data found: specifically, it presents the geometric data of the two aluminum tubes tested by BlueRobotics, their experimentally detected critical buckling pressures, and finally, the theoretical results derived from the FEM analysis and Equation (1). The number of lobes is the same in each of the three cases. It is reasonable to assume that the tested cylinders are almost free from imperfections in addition to having an excellent geometric ratio.

Experimental Tests on Domes
Some experimental tests were carried out in the MDM Lab pressure chamber, as shown in Figure 11. The technical characteristics of the pressure chamber are  Each component, that is, each dome with different features, was tested together with other more resistant components to avoid influencing the test with other collapses. Each test was performed three times, and the same result was obtained in each trial: • two domes were held together by 8 screws, divided by a protective flat seal ( Figure 12); • the hydrostatic pressure was raised step by step and maintained; • the collapse pressure was detected using the manometer (just near the buckling collapse in Figure 13); • after depressurization, each element was checked ( Figure 14 shows an example of yielding collapse occurring close to the flange).   The obtained results are summarized in Table 6 (t s denotes the nominal thickness of the slab, R denotes the mean radius, h denotes the height of the dome, and R f is the external radius of the flange).

Design of Domes
It has been observed that the dome geometry can collapse in different manners depending on certain parameters. First, two equations from the literature are reported [24,39,40].
The first equation calculates P crit , the theoretical critical buckling pressure for a thin sphere: The second equation is related to the critical static compression pressure (based on large-deformation theory of shells) that causes membrane yield in the sphere (P c ): Since these equations represent two different failure modes, each one is valid for a specific field of application, and it is worth noting that their results could not be the same. Compared with the previous notation, it is important to highlight that As one might expect, a hemispherical dome has a better buckling resistance than other kinds of domes. Indeed, hemi-ellipsoidal prolate and hemi-ellipsoidal oblate domes are ten times less robust than hemispherical domes. Experimental tests [21] have shown that the structural failure of domes occurs for values of critical pressure lower than those theoretically predictable by the formulas above. The literature reports many specific correlations but only for very limited fields of application [6]. The aim of this work is to extend the field of application of the available values.
The comparison between the theoretical data and the experimental data from the previous section is not satisfactory (Table 7). Hence, there is a need to introduce corrections to predict the behavior of the case under study. This process starts from the theoretical critical buckling pressure of a hemispherical dome (P H ) and the ideal yield pressure for hemispherical and homogeneous domes (P y ). In this study, the variation in thickness, the aspect ratio, and the type of constraint were taken into account, given their non-negligible influence. The collapse pressures obtained from experimentation are generally less than those calculated using theoretical formulas because of the geometric imperfections. Furthermore, in the case studied, the geometry strongly depends on the manufacturing method.
The chosen domes are made of plastic by thermoforming. The structural characteristics of the objects realized in this type of production process vary according to the geometric characteristics, and they differ considerably from those produced by other production techniques.
Unfortunately, therefore, the experimental theories derived from case studies in the literature [21,41], although suitable for highly regular domes, cannot be directly applied to thermoformed PMMA. Thus, the theories and equations were modified, verified, and compared with the results of the FEM analysis.
In the flange part, the thickness, t s , is the nominal thickness of the plexiglass plate, but the thickness decreases as it approaches the center of the dome, as shown in the chart in Figure 15. Hence, for safety, the thickness used is t H , which is the minimum of the dome (to clarify, in Table 8, the geometrical data are as follows: R is the mean radius, h is the height of the dome, and d is the external radius of the dome; all three variables do not consider the flange). The blue dotted line is the interpolation performed in this study (Equation (10)) [42].
where S 1 , S 2 , h, and d are geometrical quantities and shown in Figure 15; S 1 = t s and S 2 = t H ; and s 1 ,   If the changes in the parameters are considered to be independent of each other in the first approximation, then the variation in the geometry, thickness, and constraints can be analyzed.
It has been observed that an imperfectly spherical shape influences the mode of dome breakage: as the thickness increases, the influence of the membrane buckling behavior decreases and favors yielding, while the point of breakage moves from the top of the dome to the flanged base, close to the constraint. Furthermore, in the thickest domes, the radius of the inner curvature between the dome and the flange greatly influences its breaking behavior under compression, but this geometric parameter depends on the thermoforming process and thickness: the greater the thickness, the greater the radius and, consequently, the greater the mechanical stresses of the section. This is because the dome tends to overturn. On the other hand, as the radius of the external curvature between the flange and the dome decreases, local tensions increase. This justifies a sublinear trend. Other differences, which are discussed in the following part, are the AR (aspect ratio), distance from the unit, the presence of holes in the flange, and an imperfect hemisphere.
The purpose of this study is to derive a simplified and efficient method for the case under study and integrate the problems related to the aspect ratio and constraints into the calculation to identify the limit between compression and buckling.
The formulation of the theoretical critical buckling pressure (Equation (8)) can be corrected using the AR (R/h), which is dependent on the parameters of the dome and non-negligible in this particular case.
In Figure 16, a graph from the literature is given in which the colored dotted lines represent the interpolation performed in this study. Below are the relationships used to interpolate the graph. First is the initial linear part (based on a linear function) (Equation (11)); second is the nonlinear part (based on the power curve) (Equation (12)), and third is the scaling based on R/t H (based on symmetrical sigmoidal) (Equation (13)).
As regards the border between compression and buckling behavior, the parameter λ d is considered. Starting from the state-of-the-art literature [44], λ d (Equation (14)) expresses a ratio between pressures and a relationship between geometric and material features.
where P y = ideal yield pressure of hemispherical and homogeneous shell.
From the experimental results, it was observed that if λ d < 1 (Equation (15) is valid), then the rupture occurs because of the static compression load. On the other hand, if λ d > 1, then the buckling phenomenon occurs (Equation (8) is valid).
It is highlighted that the mechanics of the fracture depend on the type of constraint, the material, its stiffness, and the specific production method.
For example, with consideration of only the constraint, in the particular case under examination, there are holes for connecting screws, and a peak of stress in their surroundings causes a crack and a brittle break that promote the initiation of the collapse.
A corrective factor is introduced ( f y ) to take into account these effects. In this way, the theoretical corrected pressure (P yc ) matches the experimental one. The correction factor is based on a power curve Equation (16) derived from FEM analyses, and it is dependent on the dome thickness.
with y 1 = 0.383 and y 2 = 1.034. The influence of the correction of the coefficient λ dc , which is the same as λ d but with corrected pressures, follows a square root law, while the coefficient f y follows a linear one. Therefore, the transition point of collapse is also affected. As a result of applying this coefficient, the theoretical data coincide with the experimental evidence, as summarized in Table 9. Obviously, similar to the case with the cylinder, a safety factor for the dome (e.g., SFd = 1.25) has to be considered to obtain the pressure design. This starts from the minimum pressure between P designDY for yielding and P designDB for buckling, the pressure at which the vehicle could dive safely.
As discussed, FEM analysis was performed, although with different thicknesses from those experimentally tested (Table 10), to both validate the proposed relationships of buckling pressure and delineate the correction curve for theoretical compression breakage (in the column of the FEM analysis results, only the real collapse pressure is reported; it is a sublinear trend, as is the experimental one). Given the difficulty of identifying the yielding breaking state, a system of constraints suitable for the case under examination was to set up by using experimental data. In particular, the static analysis of the 10 mm thick dome is reported in Figure 17. The results of the intermediate domes were used to define the corrective relationship. The bolts in stainless steel AISI 316, A470 EN ISO 3506, were inserted as connectors and simulated with their pretension clamping force set to 3760 N, which is equivalent to nominal tightening. The lower face of the flange was fixed, while the contact between the dome and the flange was imposed without penetration and with a friction coefficient of 0.5. A tetrahedral solid curvature-based mesh was chosen, with 202,760 elements and 317,983 nodes. The minimum theoretical value between buckling and yielding collapse pressures are in bold and stress the type of break, and the result is the same as the experimental result. Specifically, the collapse occurs because of buckling up to a thickness of 5 mm; beyond this thickness, the collapse is due to yielding. The most significant comparisons of breakage examples are qualitatively reported, and the experimental and FEM results match. In Figure 18, the intermediate (6 mm thick) domes broke in the lowest part because of the yielding phenomenon, and in Figure 19, the thinnest (3 mm thick) domes broke because of membrane buckling in the highest region.  Notably, both the real and the simulated domes are not perfectly hemispherical. With this corrected shape, simulations can be used to predict the real breaking position, as well as the collapse pressure. This characteristic depends directly on the thermoforming process. Moreover, the AR factor and radius of the curvature are very important for numerically reproducing its behavior.
The FEM results well reproduce the experimental data; consequently, they are considered reliable.
In the end, given the correlations reported, a 10 mm thick dome for our vehicle FeelHippo AUV was designed, as shown in Figure 20, in which the left shows the real dome, and the right shows the design chart. The vehicle is thus able to reach a depth of 30 m safely (the less resistant structural component is the cylinder).

Conclusions and Future Developments
This paper presents theoretical and experimental investigations to produce a consistent relationship between both the classic naval design and innovative methods. A good correspondence between the experimental evidence and the optimized simulations is reported.
Regarding the design of the cylindrical tube, the geometry has been widely studied, and the related theory has been largely used; therefore, it is possible to effectively design the component in PMMA, too, even without using FEM analysis in the initial project step.
On the other hand, for the dome geometry, the standard theoretical approach is too generalized, and the specific results are influenced by boundary conditions, such as geometry and the linking constraints. Therefore, the predictive capacity of the algorithms was improved by using the available data in the literature to propose and introduce some correction coefficients. With the coefficients applied, the obtained results are very satisfactory.
Furthermore, because of the buckling instability phenomenon, the formula for sizing the dome was combined with the aspect ratio factor. This was implemented because the high region of the dome is affected by this kind of break, and the breaking mode largely depends on the geometry. The other formula predicts the compression break, which occurs in the area of the dome closest to the constraint. Given the strong dependence on the type of constraint and the production methodology of the dome, which even determines the radii of curvature at the base, an overall coefficient was identified for the specific case of a thermoformed plastic flanged dome fixed to the rest of the vehicle by screws.
As a result of the changes made, structural resistance to buckling and compression seems to be provided in an effective and correct manner for the specific case under study. The new dome was tested, and the results are encouraging.
Possible future developments are related to the evaluation of specific coefficients for each application, individual variation in the thickness, the constraint, and the manufacturing process. This requires a much broader experimental campaign that allows each of the many coefficients to be studied separately. Because of the costs of shell buckling experiments, the effects of the large number of variables to be considered could be analyzed first by a combination of experimental and numerical approaches.
For the thermoforming process, it could also be determined whether it is possible to define an equivalent thickness, starting from the original thickness and taking into account the production method, for other types of plastic material. Finally, the ultimate goal will be to develop a general method and extend it to different applications.