3.1. FEM Simulations
FEM simulations have both qualitative and quantitative components. The qualitative component consists in the identification of the most stressed areas, which determines the scheme of the sensors’ arrangement. The quantitative component consists in estimating the value field concerning the primary strain analytical calculations.
In the FEM analysis, the following simplifying hypotheses were considered: the thickness of the material was uniform for the entire thin wall of the sphere and the thickness of the wall was the smallest value measured using an ultrasound technique.
The software used for modeling and simulation was Dassault Systèmes SolidWorks 2014 (Watham, MA, USA) [
21].
It is a well-known fact that SolidWorks represents a complex software tool for design and simulation purposes; to emphasize this aspect, one can quote J. Ed Akin’s work,
Finite Element Analysis Concepts via SolidWorks, 2009: SolidWorks “studies calculate displacements, reaction forces, strains, stresses, failure criterion, factor of safety and error estimates. Available loading conditions include point, line, surface and thermal loads, elastic orthotropic materials are also available. The SW Simulation software also offers several types of nonlinear studies” [
20].
Application of SolidWorks Simulation 2014 allows for several approaches for 3D modeling and simulation:
(1) A first criterion for the selection of meshing is the use of the “part” type elements (individual parts) or “assembly” (entity made of several individual parts) (see
Figure 1 and
Figure 2).
(2) Second, a simplification can be made such that the simulation time is shortened and less computational resources are used by using the symmetry option (circular, in this case) when 3D modeling and meshing the element that is going to be studied.
For our mesh tank simulations, by using a single element (using option “part” without the use of symmetry), one can obtain the model in
Figure 1a. Some of the existing disadvantages are the impossibility of choosing different materials for the elements, namely the spherical tank and its support legs, respectively. Another disadvantage is the existence of some “residues” (see
Figure 1b), as seen on the transversal section of the element from
Figure 1a; one can observe the traces of the intersection of the legs of the spherical tank with the tank itself (in the center). We call these traces 3D modeling residues; these residues will cause problems related to the application of the FEM method (with implications for the difficulty of meshing and, implicitly, for the degree of accuracy).
As a first measure for raising the accuracy of the treatment, one can employ 3D modeling with the use of SolidWorks symmetry. This method allows for the choice of several variants, among which is the algorithmic treatment with the help of circular symmetry, which we consider to be of maximum efficiency in our case in terms of the number of meshing elements used and the computational simulation time. The circular symmetry allowed us to use just a slice (
Figure 2a) of the tank element for meshing and simulation (as represented in
Figure 1c) from a total number of ten slices in which the tank was divided. The input data of the simulation are presented in
Figure 2b,c. The results obtained in this case were expected to be closer to the real values because this method allowed fora larger number of finite elements because the modeled slice portion was smaller compared to the one for the whole tank from
Figure 1a.
According to SolidWorks representation rules (
Figure 2a,b), different red-sized arrows stand for forces or pressure fields and green custom arrows and indicators stand for boundary conditions in terms of degrees of freedom (displacements) as supporting scheme.
In
Figure 2b, the loads coming from internal pressure, weight of the contained LPG, and weight of the tank itself are represented. In
Figure 2c, one can notice the local coordinate system, denoted Coordinate System1 (
Figure 2c), in relation to which the law of linear variation of the load could be modeled according to the actual weight of the LPG stored in the liquid phase by considering an 80% full container (according to the working manual instructions, a tank will be filled only 80% from the total volume, where the rest should remain empty for vaporization phenomena). The FEM parameters are represented in
Figure 2c.
We performed simulations with one slice as represented in
Figure 2a–c. In order to validate and confirm the results obtained with one slice, we used also the option of “assembly” for meshing and simulation. In this case, we used the whole set of ten slices without using the method of simplifying the calculation scheme by using symmetry (see
Figure 2). As noted in
Figure 2c and
Figure 3c, the mesh grid execution time for one slice with options “part” and “symmetry” was approximately 10 s compared to 1.5 min for an assembly composed often slice modules. Moreover, one can see (
Figure 3b) a better 3D modeling quality (as compared to
Figure 1b) due to the total absence of residual artifacts.
The use of the SolidWorks symmetry algorithm for the FEM 3D modeling and simulations of the large structures considered makes this work original; to the best of our knowledge, for complex and large structures, this slicing modeling technique is unheard of in the scientific literature.
The results of the FEM simulations (option “assembly”) are: the values of the equivalent strain field (
Figure 4a), the displacements components field (
Figure 4b), the von Mises equivalent stresses (
Figure 4c), and the maximum von Mises equivalent stresses field (
Figure 4d). In the mechanics of materials, the von Mises yield criterion (also known as the maximum distortion energy criterion) can be formulated in terms of the von Mises stress or equivalent (von Mises) tensile stress, which represents a computable scalar value of stress; a material is said to start yielding when the von Mises stress reaches a value known as the yield strength. The von Mises stress is used to predict the yielding of materials under complex loading from the results of uniaxial tensile stress. Based on a known plane stress state defined by principal normal stresses
the equivalent von Mises stress is given as Equation (1):
3.2. Experimental Method
In what follows, we present the results for the experimental method.
Using the FEM simulations, we determined the zones (the area of intersection of tank with supporting legs) where the resistive strain sensors were to be placed. Using the sensors mounted in the determined places, we obtained the experimental data. In this section, we compare the theoretical von Mises equivalent stress obtained using the FEM method with the experimental ones.
The resistive strain sensor electro-resistive measurements method (measuring electrical resistance variations undergone in a strain gauge sensor grid linked to specific resistive strain sensor bridge equipment) represents one of the most frequently used experimental techniques and targets the real behavior of the material of the mechanical structure that is being investigated. The option to apply this method is based on the fact that research can be carried out on the structure under actual operating conditions.
It is well known that for materials situated in the elastic behavior limit (Hooke’s theory), there is a linear relationship between specific deformations (strains) and stresses. Above this limit, plastic deformation occurs and the relation between the specific deformations and stresses is no longer linear; moreover, the equations that express the stress/strain connection become very complex.
Any specific structural element deforms under a load, i.e., the existence of a normal/tangential field stress (σ, τ). Experimental direct calculation of such tensile parameters is impossible; therefore, in order to address the issue, one can experimentally determine the corresponding strain field and then, based on theoretical relationships (Hooke’s theory) between specific deformations and stresses, establish the stress field values.
The main disadvantage of the experimental method resides in the fact that it does not directly identify the most strained areas of the equipment, a situation that may be overcome with the help of the previous FEM study.
Starting from the geometrical and constructive particularities of the tank, and in order to be able to determine the von Mises stress field state, it has been decided that for each measuring point, a three-directional strain sensor (rosette type) should be used (three-directional strain sensors give the strains according to three directions and this can be used to characterize the principal stresses and their directions; other types of sensors are not adequate for this). The location plan of the resistive strain sensor is shown in
Figure 5.
In
Figure 5, the original constructive design of the tank is presented together with the sensor placements. As observed in the figure, the constructive design is also formed by slices that are not correlated with the ones from the theoretical FEM simulation (the only common part is that both slices include the area of intersection between the sustaining legs and the tank). In
Figure 5, it can be observed that only five sensors were used. This was due to regulatory measurements rules for such tanks and due to the consideration of costs. The policy was to first place a smaller number of sensors, and only if abnormal irregularities were observed, more sensors would be deployed on the tank and the measurements would be repeated. In our case, five sensors were considered enough to be used for the first series of measurements. As determined from the FEM simulations sensors were placed in the area of the intersection between the legs and the tank (as seen in
Figure 5).
In order to perform load-related strain field assessment, a rosette-type sensor with three measurement directions (
Figure 6) was used, with the corresponding temperature characteristics shown in
Figure 7. The sensor manufacturer recommends the use of temperature correction curve 2 from
Figure 7 for a HBM 6/350 CRY81-3L-3M-type sensor (Hottinger Baldwin Messtechnik Gmbh, Darmstadt, Germany).
3.3. Comparison Between the FEM Simulations and Experimental Results
Table 3 indicates the sensor-recorded values (equivalent von Mises stress) with respect to the most strained/stressed point of the supporting legs. The first five columns correspond to the five points from
Figure 5, while the sixth column corresponds to the FEM-calculated values. Characteristics of the nineteen tanks are given in
Table 1. The calculated FEM stress from column 5 was the one obtained with the “assembly” option. Given the symmetry, the stress was the same for all ten FEM slices; therefore, we put just one value in
Table 3. Because we did both with the “assembly” and “part” options, we took the highest stress value from the two FEM simulations to emulate a worst-case scenario; for all cases, although the values were similar, the highest values were given by FEM simulations with the “assembly” option.
From
Table 3, it can be observed that FEM values and experimental values for points P2, P4, and P5 were similar, while P1 and P3 were not. For validating this statistically, we used Tukey’s test. Through this test, which verifies the equality of the means for “
k” selections of possibly different volumes, we aimed to identify the similarities and differences between the values measured experimentally for each supporting leg and the value calculated using the FEM.
The means
were calculated using the Equation (2), which determines
and
:
where “
k” is the number of different volumes involved and
j the number of selected measurement points.
The formula for Tukey’s test is:
where
is the larger of the two means being compared,
is the smaller of the two means being compared, and
SE is the standard error of the sum of the means, which is calculated using Equation (4):
where
S12 and
S22 are the dispersions for the two corresponding selections
and
.
Appendix A contains the raw data and auxiliary information.
The computed “q” value was compared to a certain value obtained from the standardized range distribution “qa” for n–k degrees of freedom, where n = 19, the total number of investigated equipment, and k = 2. One can consider qa to be the critical value. It was considered that the values were different if the computed q value was greater than the critical value qa.
The results are presented in the
Table 4.
The values obtained for
qa with 17 degrees of freedom are presented in
Table 5.
Using Tukey’s test, it was concluded that P2, P4, and P5 value were consistent with the FEM analysis. A significant difference was observed for the values corresponding to points P1 and P3 in relation to the FEM-calculated value. In all cases, the most stressed point was P1, which corresponded to the position of the leg attached to the tank access stairs. Accordingly, the less stressed point was the diametrically opposed one, which was P3 in all cases. Excluding the most stressed point (P1) and the least stressed point (P3), it can be seen that the values for P3, P4, and P5 measured using the tensor-resistive sensors in three directions were within the margin of ±3% compared to the values calculated using the FEM method. This was quite remarkable and doubly validates both the FEM simulations and the experimental method (one validates another).
The difference between the values recorded by sensors at P1 (adjacent to the tank stair access structure) and the opposite one (P3) was quite significant, with this phenomenon being systematically reported in all cases of the investigated spheres. As noted, P1 was the most stressed point while P3 was the least stressed. The access tank stairs was adjacent to P1, placed between P5 and P1. In this context, P3 can be considered the opposite point to P1 (and not also P4), where the position of the access stairs explains why P3 was the less stressed (and not also P4). One result of the study is the fact that the presence of stairs causes a peak stress value point at P1. Such an unconventional structural behavior was caused by specific interface areas (between the stairs and tank sphere) that are prone to stress concentration due to a sudden local increase in terms of general mechanical stiffness caused by the stair’s presence, with a corresponding undesired influence from the point of view of the mechanics of materials.
Hence, from the point of view of the general stress state field, the sensor data in
Table 2 indicate an unacceptable trend concerning the structural fatigue strength (in point P1) with poor life expectancy characteristics.
One can ask why we did not model the stair structures using the FEM. Our starting point for FEM modeling was the European regulatory norms, as defined in EN 13445 (similar in USA and worldwide), which is particular to our tank spheres. This regulation does not include the stairs in the mechanical design of the sphere; therefore, we did not include the stairs in the FEM modeling. However, the experimental results show clearly that the stairs had an influence on the stress behavior of the tank; however this, as noted, is not regulated. We will enlarge this discussion in the next section.