Next Article in Journal
Data Fusion for Smart Civil Infrastructure Management: A Conceptual Digital Twin Framework
Previous Article in Journal
Study of Mechanical Response of Tunnels Crossing Active Faults in Different Burial Depths
Previous Article in Special Issue
Non-Linear Analysis and Retrofitting by Steelwork of a Precast RC Warehouse
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Non-Linear Modeling of New Wind Girder Used for Bolted Tanks

1
Department of Applied Mechanics, Faculty of Mechanical Engineering, VSB—Technical University of Ostrava, 17. listopadu 2172/15, 708 00 Ostrava, Czech Republic
2
WITKOWITZ ENVI a.s., Ruska 1142/30, 703 00 Ostrava, Czech Republic
*
Author to whom correspondence should be addressed.
Buildings 2023, 13(11), 2724; https://doi.org/10.3390/buildings13112724
Submission received: 8 September 2023 / Revised: 22 October 2023 / Accepted: 27 October 2023 / Published: 29 October 2023
(This article belongs to the Collection Non-linear Modelling and Analysis of Buildings)

Abstract

:
Large-capacity bolted cylindrical tanks for liquid storage are used in many applications. The tanks are made of thin steel sheets that are connected by bolts. A common problem associated with tanks is deforming under extreme loads. Adding wind girders to the tank increases the tank’s buckling capacity, which is defined as the limit load at which the structure loses stability. The girders are usually placed in the horizontal joints of the tank wall. The girders are bent from standard or non-standard steel bars with a uniform cross-section. This type of design is difficult to produce, especially with large profiles or large curvatures, to avoid distortion of the cross-section during bending. Furthermore, the girders are customized to the given openings and curvature for various tank diameters. The resulting solution is then uneconomical and more complicated to store. This paper deals with the design and non-linear modeling of a new shape of wind girder for bolted tanks that eliminates the above-mentioned disadvantages. To analyze the new shape of the girder, a non-linear numerical model of an open-topped tank with various dimensions is designed to study its buckling capacity.

1. Introduction

In many technical applications, above-ground steel tanks are used for the storage of various liquids or bulk materials. Such tanks have a relatively simple geometry. They consist of steel sheets that are bolted or welded together. Tanks are produced with various types of roofs or open-topped. Commonly, there is a risk of buckling under load due to tanks being made of thin steel sheets. For this reason, wind girders are mounted to the tanks. In the case of open-topped tanks, the risk of buckling is more significant due to negative pressures (suction), as mentioned by the authors Uematsu et al. [1]. The girders are used to increase the tank’s stiffness and buckling capacity. In some cases, the girder is mounted only on the top of the tank. However, using more grinders is practical for increasing buckling capacity. The number of girders and their placement are determined by the design according to the specific type and location of the tank. It is necessary to provide the tank wall with a special coating, e.g., enameling, anti-corrosion coatings, or use stainless steel as a material. Potential corrosion also significantly reduces buckling capacity, which has been investigated by the authors Zhao and Lin [2] and Shokrzadeh and Sohrabi [3]. Another risk of buckling occurs during the construction of the tank, as discussed by Jaca and Godoy [4], and therefore, recommended principles during construction should be followed.
Bolted tanks consist of sheets with drilled holes, which are also used for mounting the above-mentioned wind girders onto the tank wall. Each country defines the standards for the design of cylindrical above-ground tanks. The European single market is subjected to Eurocodes. For tank design, it is possible to use Eurocode 3 Part 4-2: Tanks [5] and for loads Eurocode 1 Part 4: Loads on tanks and reservoirs [6] together with others. AWWA D103-09 [7] is the most commonly used standard for the design of bolted tanks in the American continent. When designing tanks, it is necessary to take into account in particular (I) the load from the storage medium, (II) the wind load, (III) the dead load, (IV) the snow load for roofed tanks, (V) the operating load, and (VI) possible seismic effects. In this paper, the authors focused on the use of Eurocodes and the wind load (II) because wind girders are mainly used to prevent buckling under wind action.

1.1. Procedures for the Evaluation of Tank Buckling Capacity

The load-bearing capacity of the tanks is mainly affected by the hydrostatic pressure load and seismic effects at low wind speeds—i.e., a full tank. Wind load greatly affects buckling capacity, especially for higher wind speeds and empty tanks. The bottom section of the tank usually has thicker sheets due to the higher hydrostatic pressure load, so the wind girders reinforce mainly the top section.
The wind load can be implemented in the numerical model using an equivalent axisymmetric pressure-idealized (see Figure 1a) or the real pressure distribution based on the relevant standard, testing, or numerical simulation (see Figure 1b). The Eurocode EN 1991-1-4 [8] specifies a coefficient k W to convert the real pressure distribution to the axisymmetric pressure. The substitution of this distribution is used by Uematsu et al. [1]. The authors Zhao and Lin [2] and Burgos et al. [9] studied both types of wind pressure distribution, and according to their study, axisymmetric pressure is conservative for the buckling of tanks. However, this simplified pressure distribution is only used for analytical designs or initial analyses.
Tanks can be classified as thin-walled structures; therefore, membrane theory can be used to evaluate tank limit states. Membrane theory is an analysis that predicts the behavior of thin-walled structures under continuous loading, assuming only membrane forces in equilibrium with the external load. The assumptions of membrane theory are published in the literature [10] and cannot account for many real cases. In current practice, numerical models using the shell theory with different accuracy are used. Types of the shell theories specified in EN 1993-1-6 [11] are listed in Table 1. Four limit states, which should be considered when evaluating tanks according to EN 1993-1-6, are
  • Plastic limit state;
  • Cyclic plasticity limit state;
  • Buckling limit state;
  • Fatigue limit state.
Because this paper is focused on wind girders, the authors mainly dealt with the buckling of tanks. According to EN 1993-1-6, buckling is defined as “the limit state at which the structure suddenly loses stability under membrane compression or shear”. For these tank analyses, numerical models using the finite element method (FEM) are most commonly used. The least accurate and unsuitable in relation to the stability of thin-walled structures is linear elastic shell analysis (LA). The implementation of LA using a numerical model was studied by Kminek et al. [12]. The most accurate results are given by geometrically and materially non-linear analysis with imperfections (GMNIA).
Elkholy et al. [13] investigated the use of six different shell elements in Ansys Mechanical (Ansys, Canonsburg, PA, USA) for tank analysis considering structure–fluid–soil interaction. They used the experimental model and theoretical solution obtained by Mazuch et al. [14] to validate the numerical results. Two cases, an empty and a filled tank, were tested. The performed analyses showed that the number of elements significantly affects the accuracy of obtained results.
A parametric study for four different large-capacity open-topped tanks was presented by Pan and Liang [15]. They performed linear and non-linear numerical analyses with imperfections. The results show that the main effect on the buckling under wind action is due to the pressure of the windward side, and the use of girders can significantly increase buckling capacity. Chen and Rotter [16] performed an analysis of cylindrical shells with uniform thickness. They considered a stocky cylinder, an intermediate cylinder, and a slender cylinder. In their paper, they found that different dimensionless length parameters play the main role in buckling capacity, and plasticity has a minor role for cylinders with R / t 500. Cao and Zhao [17] carried out a parametric study on the buckling capacity of filled silos with various slendernesses. They applied an axisymmetric pressure, which allowed them to simplify the model of the tank. They performed six different linear and non-linear analyses to study buckling capacity. The analyses showed that material non-linearity strongly influences the behavior of all selected tank slendernesses. The distribution of the wall thickness over the height of the tank also has a significant additional influence on buckling capacity. However, the study by Zhao et al. [18] described that material non-linearities have only a relevant effect when modeling the buckling of a filled tank. In the case of an empty tank, this effect is negligible. However, the implemented imperfections have a significant effect, which reduces buckling capacity in both cases. For this reason, linear bifurcation analysis (eigenvalue buckling) is suitable as an initial indicator of buckling capacity. Eigenvalue buckling shape should be used as imperfections for the following non-linear analysis, which is also recommended by Maraveas et al. [19].
Dynamic analysis of open-topped tanks has been discussed by Chiang and Guzey [20], who performed a modal analysis and obtained mode shapes and natural frequencies. Modal analysis of tanks is important to determine a critical wind speed. The critical wind speed occurs when the vortex shedding frequency is close to the natural frequency of the tank. To estimate the vortex shedding frequency, the Strouhal number is used. Strouhal number is a dimensionless number, which gives the ratio of vortex shedding frequency multiplied by characteristic dimension and flow velocity. A more detailed description of vortex shedding and the Strouhal number for cylinders is given in the publication by the authors Alam and Zhou [21]. An experimental determination of the Strouhal number in a wind tunnel was performed by the authors Ellingsen et al. [22].
Shokrzadeh and Sohrabi [3] investigated the effect of corrosion on buckling. In their paper, they considered both cases of wind pressure distribution and observed that in the case of non-uniform pressure distribution, the effect of corrosion is more significant. The distribution of the wind pressure acting on the tank is also influenced by whether it is an isolated tank or several tanks located close to each other. The effect of several tanks placed behind each other with different distances between them has been studied by Godoy [23] or also Burgos et al. [24].

1.2. Description of Wind Girders

In engineering applications, stiffeners are divided into vertical and horizontal, as shown in Figure 2. Vertical stiffeners (see Figure 2a) are mainly used to increase stiffness in the lower part of the tank or areas with high axial loads. However, these stiffeners are not used very often because of their efficiency. Horizontal stiffeners (Figure 2b) are used to reinforce the tank due to wind action and to prevent buckling. Therefore, they are called wind girders. Wind girders usually take the form of round-edged steel bars of various dimensions, as illustrated in detail A of Figure 2c. These girders are mounted on the top of the tank and eventually between selected rows of the tank wall. The location and dimensions of the girders are specified in the design process. The use of appropriate girders, including their correct placement, increases the buckling capacity of the tanks, which was also confirmed by the authors Pan and Liang [25]. Standard EN ISO 28765 [26] defines the minimum moment of inertia to be used for the upper girder and also for the intermediate girders, depending on their spacing, the loads, and the dimensions of the tank. In the case of conventional girders, the moment of inertia (as well as the stiffness) of the girder is uniform around the circumference of the tank. The dimensions of the girders can be thus easily determined. The most common placement of girders is the upper part because the thinnest sheets are usually placed from the upper edge.
Conventional girders with uniform cross-sections are not very efficient in terms of production technology, especially for large cross-sections of girders or large curvature of tanks. Uematsu et al. [1] have analyzed the stiffness of girders and they show that there is a limiting “stiffness” value of girders. When increasing the stiffness of the girder above the limit value, there is no longer a noticeable increase in buckling capacity. The aspect ratio of the tank and whether the bottom of the tank affects the overall stiffness also matter when designing girders, which was discussed by Bu and Qian [27]. The use of intermediate girders changed the single-mode shapes during buckling and also increased buckling capacity. The study of the geometrical parameters of the girders was carried out by Bu and Qian [28] and Lewandowski et al. [29]. The authors’ findings can be summarized that the design of the dimensions of the girders specified in the standards is conservative. And detailed non-linear analysis is recommended for better use of the material, taking into account all possible factors.

1.3. Possibilities of New Girders Design

As mentioned above, wind girders have a significant effect on the behavior of tanks, especially for buckling. The stiffness requirements for wind girders in open-top tanks were published by the authors Zebyek and Topkaya [30]. They propose a minimum requirement for the second moment of area with respect to negative internal pressure (suction). The possibility of replacing the girder with various modifications and the usage of boundary conditions was studied by Burgos et al. [9]. By designing the optimal type and placing them correctly, buckling capacity can be greatly increased. The main idea of the new shape of the girder is the use of a simple profile (segment) with a rectangular cross-section, which is intricately curved around the circumference of the tank. The new girder design is composed of many segments mounted together around the circumference of the tank. The bend of the segment easily adapts to different tank diameters, and these new profiles are easy to store for use on any tank. For this reason, it is easier to manufacture and customize. The total cost of tank construction (after accounting for production costs) is thus reduced. The new girder shape is used experimentally by several companies, but it is not possible to get details of the geometry, overlapping of the sheets, and design data. The new shape also needs to be customized to the sheet metal punching and designed to the producer’s requirements.
The behavior of the new girders, especially ones with variable moment of inertia around the circumference of the tank, is required to be analyzed using a numerical model. The girders are mainly designed to resist wind load. This means that it is necessary to assess the condition in the case of an empty tank and the application of external pressure. The axisymmetric pressure around the circumference of the tank is conservative for buckling; therefore, this pressure distribution is considered in this paper.
This paper aims at the non-linear modeling of a new type of horizontal wind girder concerning the above-mentioned disadvantages of profiles with uniform cross-sections. For this purpose, a numerical model of the tank with the new type of girder is designed. The numerical model consists of buckling analyses (using GMNIA) and modal analyses for various tank dimensions. By using a more efficient and reliable shape of the girder, the cost of the tanks is reduced through easier production technology. Preliminary results were published by the first author of this paper in his master’s thesis [31].

2. Materials and Methods

This section is divided into 3 parts. In the first part, a description of the design of the new wind girder (NWG) is presented. In the second part, the numerical model of the tank is designed. This part focuses on the material model, boundary conditions, and solution assumptions. To study the behavior of the new design, buckling capacity of the tank with NWG is compared with the tank with round-edged equal leg steel angle bar EN 10056-1 [32] - L60 × 60 × 6 - EN 10025-2 [33]—corresponds to a = 60 mm, t g = 6 mm (L60) (see Figure 2c). Replacement of L60 is selected based on experience of the company WITKOWITZ ENVI a.s. The same material-steel S235JR EN 10025-2—is considered for NWG and L60. In the third part, a parametric study of low-height to high tanks is presented.

2.1. Design of New Wind Girder

The innovation in NWG compared to the conventional girder is the use of two differently curved thin sheets (V-segment and P-segment) adapted to the tank wall (tank radius is R). A shape of NWG formed by segments is shown in Figure 3a. The illustrated holes are intended for mounting bars or segments to the tank wall. The displayed configuration is repeated around the circumference of the tank, and thus the whole girder is formed. The shape of a conventional girder is shown in Figure 3b. The bar of the given length is bent to the required curvature, and then it is mounted to the tank using the depicted openings.
The detail of NWG is shown in Figure 4. The V-shape segment of height h (Item 1) is further supplemented by the plate (P-segment) of the same thickness t and width w (Item 2). The segments are further bolted (Item 3) together and attached to the tank wall (Item 4) to form a trapezoidal girder shape. All these segments have a simple rectangular cross-section, but by overlapping each other, the thickness is strengthened. This effect leads to a local increase in the moment of inertia. According to the designed geometry, the thickness is increased up to four times in some regions.
During manufacturing, distortion of the V-segment and the P-segment does not occur as in conventional profiles. In the case of insufficient rigid girders in the design of the structure, the stiffness of the girder can be easily increased by modifying the selected parameters (especially h and t), and thus a higher buckling capacity can be achieved. Therefore, it is possible to manufacture many segments of NWG according to tank size for efficient mass production.
Standards for designing bolted tanks define the minimum moment of inertia or the section modulus of the wind girders [7,26]. The main difference in the new shape is the fact that NWG does not have a uniform moment of inertia in the circumferential direction like conventional profiles. For this reason, it is advisable to investigate the behavior of the girder by using a numerical model for a reliable design. The proposed basic dimensions of the segment forming NWG are listed in Table 2. Dimensions are designed using a numerical model and adapted to the manufacturer’s structural capabilities.

2.2. Design of Non-Linear Numerical Model

Large-capacity tanks can be considered as thin-walled structures, satisfying the condition R / t w 10, where R is the radius of the tank, and t w is the thickness of the tank wall. For buckling capacity study, geometrically (large deflections and strains) and materially non-linear analysis with imperfections (GMNIA) under external pressure is performed. To implement the imperfections, linear elastic bifurcation analysis (LBA) was used. Displacements in each node were exported, and further, this mesh was imported as an initial shape in GMNIA. The imperfection amplitude was applied to u i = 0.5 t w . The influence of the imperfection amplitude has been studied by many authors [2,15,24,34]. In this paper, a sensitivity study of imperfection amplitude was performed for three variants with NWG (TK1, TK4, and TK7), which varied a lot from the model with L60. Results of this study are presented in Section 3.5. Further, a modal analysis (MA) of all investigated tank types was performed. The mode shapes and natural frequencies are important factors as they can affect the possible vibrations of the tanks during airflow and vortex shedding. The natural frequencies of the tank should be sufficiently different from the vortex shedding frequency. The authors of this paper have many years of experience in the application of numerical methods to solve the behavior of complex structures [35,36,37,38,39]. For this reason, the authors used FEM and designed a parametric numerical model of the tank in ANSYS Workbench 2022 R2 (Ansys, Canonsburg, PA, USA).
A 3D numerical model combines two types of finite elements (shell and beam elements are used). A 4-node shell element SHELL181 with linear displacement approximation using 6 degrees of freedom per node (3 displacements and 3 rotations) is used for the discretization of the tank wall. The girders are modeled using the 2-node beam element BEAM188 with linear displacement approximation using 6 degrees of freedom per node (3 displacements and 3 rotations). Because applied element types (shells and beams) use the same number of degrees of freedom, the connection of these elements is achieved through the node merge. To simplify a girder model, different cross-sections are defined for each part with different shapes. For the V-segment are defined rectangular cross-sections from 50 × 4 mm to 50 × 12 mm according to different thicknesses of the shape in Figure 4. Thus, a different number of the overlapped sheets is taken into account, and it is not required to use the mutual contacts of the individual segments. The P-segment is modeled as an additional thickness of the tank wall so that a rigid connection between the segment and the tank is defined. Due to low stiffness at the top edge of the tank, for the upper girder, the number of P-segments is doubled. For numerical model with L60 is defined as a beam with a uniform cross-section 60 × 60 × 6 according to Figure 2c.
Because in this paper a parametric study was carried out for different tank sizes (see Section 2.4), the finite element mesh parameters are given for multiple models. A sensitivity analysis of the mesh size influence on buckling capacity was performed. Therefore, the shell element edge length in the circumferential direction is set to 32.5 mm for all variants. Element size in axial direction varies for each variant and is between 25 mm for small tanks and 150 mm for large tanks. In the case of L60, the beam element length is the same as the shell element length in the circumferential direction. In the case of NWG, the beam element length varies. The connection part to the tank wall of the V-segment is divided into two elements, and the element length of the remaining parts depends on the size of the overlap of sheets and dimensions of the segment (see Figure 4). However, the preferred length of these parts is set to 32.5 mm. The total number of nodes of the finite element mesh of the model with NWG ranged from 28,815 to 598,640 nodes. For cases with L60, the number of nodes is a bit lower due to simpler geometry, but similar to variants with NWG. The finite element meshes for a model with an NWG and a diameter of 4.3 m are shown in Figure 5. Detail A shows detailed connections of nodes of tank wall and girder, which was performed by node merge in the appropriate places. One can also observe the single cross-sections of NWG, which are discretized by beam elements.
The tank wall and girders materials are considered to be ideally isotropic and homogeneous. The yield stress is typically greater for the tank wall material than for the girder material because the tank wall is usually more loaded. For girders, the yield stress and the tangent modulus play a minor role. A bilinear elastoplastic material model with isotropic hardening was used for both materials. Relation between true stress σ True and deformation ϵ is in Figure 6. The elastic part is defined by Young’s modulus E and Poisson’s ratio μ . The plastic part is from yield stress σ y to ultimate stress σ u l t with tangent modulus E T . The tangent modulus E T is calculated from the elongation and ultimate stress. Young’s modulus and Poisson’s ratio are taken as standard values for conventional structural steels specified in Eurocode EN 1993-1-1 [40]. For the tank wall, KOSMALT E 300 T (K300T) is selected with the mechanical properties given in the U. S. Steel Košice catalog [41]. The girder material is chosen S235JR EN 10025-2 (S235) with the mechanical properties given in Eurocode EN 1993-1-1 [40]. The material constants used for the tank wall and girders are given in Table 3.
To determine buckling capacity, axisymmetric pressure distribution around the tank circumference (see Figure 1a) is gradually applied. The axisymmetrical load pressure is prescribed in the range from p = 0.2–50 kPa depending on the dimensions of the tanks. The full Newton–Raphson method is applied to solve the non-linear model. The number of minimum executed substeps of the Newton–Raphson method is set to 100 and the maximum substeps to 100,000 for all variants. Maximum displacement corresponding to the current pressure level is determined in each substep. From this relation (displacement vs. pressure) are determined the slopes k and their relative deviations
Δ k i = k i k i + 1 k i + 1 .
At the moment when deviation Δ k is not positive or the method fails to converge, the buckling capacity of the tank is assumed to be reached. The critical pressure corresponds to buckling capacity, which is obtained as the previous substep. Typically, the tanks are anchored to a concrete foundation; therefore, fixed support is used at the bottom of the tank to eliminate displacement and rotation of the nodes at the bottom edge of the tank. Only the empty tank case is considered in the whole study, as in this case, there is a higher risk of buckling due to external pressure. The self-load of the tank and the girders is neglected as it has only a minimal effect on the given configuration. Applied boundary conditions are shown in Figure 7. The maximum value of pressure for this variant was prescribed to 50 kPa (see Figure 7a). The bottom edge of the tank wall was fixed as is shown in Figure 7b.

2.3. Analytical Models According to European Standards

This section presents a view on the issue of tank buckling concerning two methods of critical pressure assessment according to European standards. The following simplified analytical models were used in the next section to compare with numerical models. The first way is to evaluate the buckling according to the EN 1993-1-6 [11]. In this evaluation, the wind and subsequent external pressure have a major role in circumferential stress. According to this standard, the critical pressure using membrane theory
p σ θ , R d · t w R
can be determined from the limiting value of the circumferential stress σ θ , R d (mainly a function of tank geometrical dimensions and position of girders, production quality, or size of imperfections and material mechanical properties). For comparison purposes, quality class A according to EN 1993-1-6 was selected in the results presented in Section 3.4.
Another possibility is the use of EN ISO 28765, where the limit external pressure is explicitly defined using Donnell’s theory. This limit value is a function of Young’s modulus, Poisson’s ratio, wall thickness, distance between girders, and tank radius. However, the size of imperfections does not appear as a separate parameter that could be easily changed in this formula.

2.4. Parametric Study

A parametric study is performed on selected tanks to analyze the behavior of the girders under axisymmetrical pressure acting on an empty tank. Three different diameters D and three different tank heights H for each diameter are taken, as listed in Table 4. Tank wall thickness t w = 3 mm is considered. The range of tanks D = 4.3–42 m and H = 2.9–17.2 m is chosen based on the selection from the production list of the company WITKOWITZ ENVI a.s. The considered tank dimensions with mesh statistics and number of segments (V-segment + P-segments) per girder are listed in Table 4.
To analyze NWG behavior, the numerical model including L60 girders is designed. The model has the same mesh for the tank wall but the mesh of the girder is different due to geometry. For the reference numerical model, girders with a standardized profile are used and placed in the same positions as the new girder shape. The girders are placed above the top ring of the tank and between the middle rings. This distribution is only for comparison purposes and does not respect the design distances between the girders or dimensions according to the relevant standards. The critical pressure and associated shape from GMNIA are evaluated in both models. From the modal analysis, the first natural frequency with the corresponding mode shape is evaluated.

3. Results

This section summarizes the results of the behavior of NWG by the parametric study using GMNIA and MA. The calculated results are divided into three categories according to the height of the tanks with NWG compared with L60. The evaluated parameters are: critical pressure for the model with L60 p L 60 , first natural frequency for the model with L60 f L 60 , critical pressure for the model with NWG p N W G , first natural frequency for the model with NWG f N W G , maximum equivalent stress (HMH) in tank wall σ L 60 , w and girder σ L 60 , G for model with L60 and maximum equivalent stress (HMH) in tank wall σ N W G , w and girder σ N W G , G for model with NWG. In addition, the relative deviation between the standardized profile and the new shape for both variables is listed.

3.1. Application of New Wind Girder to Low-Height Tanks

The first evaluated set is for the low-height tanks with height H = 2.9 m (TK1, TK2, and TK3). The results obtained for TK1 presented as shape under critical pressure from GMNIA and also the first mode shape from MA for both types of girders are shown in Figure 8.
The results show that shapes obtained from GMNIA are slightly different and total deformation is significantly higher for the NWG model. The value of critical pressure obtained from the GMNIA is also higher for the model with L60 and is given in Table 5. Mode shapes from MA are quite different for both models too. This is also associated with the values of the first natural frequencies, where for the model with L60, the first natural frequency is f L 60 = 39.111 Hz, and the model with NWG is much lower f N W G = 28.035 Hz.
In the case of the larger diameter tanks (TK2 and TK3), the natural frequencies of both models decrease and also get closer to each other. The results for the TK3 variant are shown in Figure 9. As one can observe from the figure, the mode shapes are already comparable, although the natural frequency is slightly higher in the case of the model with NWG. Shapes at buckling from GMNIA are similar to TK1, but maximum deformation in this case is close to each other. Maximum deformation is also more concentrated in the upper part of the model with NWG, similar to TK1. The comparison of critical pressure and first natural frequency plotted as a function of tank diameter for the low-height tanks is shown in Figure 10. From the plotted comparison, it is obvious that both the critical pressure and the first natural frequency decrease with increasing diameter, and the biggest differences are in variant with diameter D = 4.3 m (TK1).

3.2. Application of New Wind Girders to Middle-Height Tanks

For the middle-height tanks H = 8.6 m (TK4, TK5, and TK6), the first mode shapes and the shapes from GMNIA are comparable, and only the values of critical pressure and natural frequencies are varied and are published in Table 6. The biggest relative deviation of critical pressure is in TK4, but the first natural frequency differs the most in TK6. For the variant TK5, the first mode shape from MA and the shape under critical pressure are shown in Figure 11.
Figure 12 presents critical pressure and the first natural frequency of the models with both types of girders. As with the low-height tanks, there is an obvious trend that as the diameter of the tank increases, both evaluated values decrease.

3.3. Application of New Wind Girders to High Tanks

The results for the high tanks H = 17.2 m (TK7, TK8, and TK9) are similar to those for the middle-height tanks. The shapes from the individual analyses are comparable but critical pressure values, maximum displacements, and natural frequencies vary. As can be observed from Figure 13, the first natural frequency is higher for a model with NWG for all three variants in this evaluated set. However, the critical pressure varies for each variant and the biggest difference is again for a variant with a small diameter (TK7). The mode shapes and shapes under critical pressure for variant TK7 are shown in Figure 14 and the corresponding values are listed in Table 7. For this variant, maximum displacement is higher for a model with NWG, although the shapes are similar.

3.4. Comparison of Results Obtained from LBA and GMNIA with Analytical Models

In this section, the comparison of results from LBA and GMNIA with analytical models is presented. Because LBA is not suitable for thin shell structures, this paper focused on GMNIA, which is more precise. But to compare the results obtained from GMNIA, the above-mentioned methods were also used. Critical pressure from LBA was obtained as eigenvalue, which was multiplied by the defined initial load. Values are listed in Table 8, where p L B A , L 60 is critical pressure for the model with L60 and p L B A , N W G is critical pressure for the model with NWG. Critical pressure from analytical models p E u r o c o d e and p I S O were determined based on standards described in Section 2.3. However, the different types of girders have not been taken into account in these analytical models. In the last two columns in Table 8 is the ratio between GMNIA and LBA for both numerical models. As one can observe, the critical pressures obtained from GMNIA were lower for all variants. When comparing the analytical solution and LBA, one can see that the critical pressure is significantly lower for both methods. An interesting comparison is between the analytical solution and GMNIA. The critical pressure values according to EN ISO 28765 ( p I S O ) were very close to the numerical results for both types of girders for most variants. In some cases, the critical pressure was even higher, so the standard gives less conservative results. The second analytical approach (EN 1993-1-6) is rather more conservative, which may be due to a more detailed implementation of imperfections.

3.5. Geometrical Imperfection Study

In this section, results of the influence of imperfection sizes are presented for three variants (TK1, TK4, and TK7), which differed greatly from the L60 model. These variants include a sample from the low-height, middle-height, and high tanks. The magnitudes of the imperfection amplitudes were selected based on the thickness of the tank wall. Three sizes of amplitudes were considered 0.25 · t w , 0.5 · t w , and 1 · t w and calculated critical pressures are listed in Table 9. The dependence of the critical pressure on the imperfection amplitude is plotted in Figure 15 for both models.
As can be observed from Table 9, increasing the imperfection amplitude decreases the value of critical pressure. The maximum deviation from the imperfection amplitude used in this paper (0.5 · t w ) was about 12%. Figure 16 shows load–deflection curves for variant TK7 for the L60 model and the NWG model. It can be seen from the results that the NWG model is more sensitive to the size of imperfections.

3.6. Evaluation of NWG Results

For some variants (TK3, TK5, TK6, and TK9), buckling capacity is higher for models with NWG. The smallest difference is for variant TK5, where the relative deviation is only 1.5%. On the other hand, the biggest difference is in variant TK6, where the relative deviation is 2.5%. It means that for these variants, the difference in buckling capacity is negligible. An important observation is the higher buckling capacity of a model with L60 for the rest of the variants (TK1, TK2, TK4, TK7, and TK8). The biggest deviations are in variants with small diameters (TK1, TK4, and TK7), specifically 25.1% for variant TK1. The critical pressure values are closest to each other for TK8 with a relative deviation of 7.6%. For these variants, it would be possible to increase the thickness of NWG segments from which the girder is made to increase buckling capacity. The obtained results show the variable behavior of the girders for different tank variants.
Maximum equivalent stresses in tank walls were comparable for both models for most of the variants, except small-diameter tanks (TK1 and TK4). In girders, the stress results were quite different. These differences were probably most affected by the influence of the variable stiffness of NWG and stress concentration. It means that NWG is formed by segments, and in certain places, the thickness is only 4 mm. However, in all variants, the maximum equivalent stress (HMH) did not exceed the yield stress when the load reached the evaluated critical value. So, evaluated failures were due to large deflection or zero tangent stiffness.
Load–deflection curves for two middle-height variants (TK4 and TK6) are plotted in Figure 17. One can observe that for variant, TK6 load–deflection curves were very comparable. It also corresponds to the value of critical pressures, which were almost equal. On the other hand, for variant TK4, it is evident that the curves for the L60 model are above the NWG model. Displacements were smaller at higher levels of load, and critical pressure was larger. So, in this case, the L60 model seems to be more rigid.
Also important are the results of MA, in which the values of the first natural frequency differ significantly. Especially for TK1, where the value of the first natural frequency is considerably lower for the model with NWG. The relative deviation is 28.3%, the absolute deviation is 11.1 Hz and the first mode shape also differs markedly, which is a significant difference for this type of construction. This may be due to insufficient stiffness of NWG in the upper part. If the first natural frequency had to be increased to achieve the same mode shape as the model with L60, the tank wall had to be additionally reinforced or different dimensions of NWG had to be used. On the other hand, the biggest difference is for TK9, where the first natural frequency is almost 50% higher for the NWG model. However, first mode shapes do not differ much in this case. For the remaining variants, mode shapes are already comparable to each other, and only first natural frequencies differ.
Figure 10, Figure 12 and Figure 13 show the trend of the studied variables depending on the tank diameter for both types of girders. From the figures, it can be seen that the natural frequencies of the model with NWG come out worse for low tanks. It is also noticeable that as the tank diameter increases, NWG performs better compared with L60. This reflects the better behavior for the given design of NWG for large tanks (TK3, TK6, and TK9). This makes the use of NWG easier to manufacture, especially for large-diameter tanks, where they represent a significant part of the total mass. On the other hand, L60 shows better behavior for tanks with a diameter of 4.3 m.

4. Discussion

Wind girders are mainly used to increase buckling capacity. The new design of a girder in the form of segments (Figure 3a) can be used to replace the conventional girders (Figure 3b). As the shape of the girder is designed, the basic geometrical parameters given in Table 2 can be modified to achieve higher stiffness of the girder and consequently increase the overall buckling capacity for individual tank applications. The main difference with commonly used girders is their changing stiffness in the circumferential direction, so the numerical model is needed for further evaluation. As the study is used to compare NWG with L60, the values of the loading pressure do not correspond to the specific location of the tank and the real wind distribution according to EN 1993-1-4 [8]. For a concrete application, it would be advisable to implement this and assess the influence of these specific parameters for a given location.
The presented analyses are focused on the non-linear buckling analysis of empty tanks. To assess the stress–strain behavior, it would be possible to design a more precise model and evaluate the stresses and strains in the girders and the tank wall. Only tanks without roofs are included in the analyses, as there is a higher risk of buckling due to negative pressures acting on the inner tank wall. In the case of tanks with a fixed roof, there is partial stiffening of the roof structure at the top edge, and as the tank is closed, no suction from wind occurs in the inner part of the tank. The effect of internal pressures is also considered in publications [1,25].
A study was carried out in the field of testing large-capacity tanks. Based on the obtained information, it is not possible to test the analyzed tanks even in the largest wind tunnels, e.g., in the “Jules Verne climatic wind tunnel” [42] in France. Testing smaller tanks is possible but very expensive. In the case of testing scaled models, it is difficult to perform experiments that are valid for application to analyzed tanks. This is mainly due to the availability of very thin plates and small girders and bolt connections. For these reasons, the experiment of each model is not implemented in this paper. Simplified analytical models according to EN 1993-1-6 and EN ISO 28765 were used to compare with the numerical results (LBA and GMNIA).
In this paper, the seismic load is not considered, which also has a major influence on the tank configuration if the tank is located in the seismic area. Numerical evaluation of tanks excited by seismic loads was carried out by authors Spritzer and Guzey [43]. To model the seismic load, it would be useful to design a more precise 3D model with the replacement of the fluid with finite elements and simulate the problem as a fluid–structure interaction (FSI) as presented in the publication of the authors Compagnoni and Curadelli [44]. In this case, it would be appropriate to investigate both roofed and open-topped tanks and also determine the height of the sloshing wave during the earthquake and the forces acting on the foundation. For roofed tanks with an insufficient height of the water level from the top edge (so-called freeboard), the roof could be damaged due to the effect of the sloshing wave, which is described in more detail in References [45,46]. In terms of the assessment of NWG, the influence of the occurring hydrodynamic effects and the vibration of the tank on the behavior of the wind girders could be investigated.
These analyses are focused on bolted tanks (not welded), and for this reason, it is recommended to have a more precise model with substitution of the real bolts, e.g., by beam elements, and also evaluate the stress in the bolts and the tank wall due to the net section. In this way, it would be possible to design the ideal bolt spacing and size for each tank design under known loads. To perform a more detailed analysis, it is also advisable to consider a full shell model by replacing the single sheets with shells and implementing the contacts between the single segments making up the girder and also between the girder and tank wall. In this paper, simplification is applied and the nodes at the contacts between the girder and the tank are merged. However, for tank analyses, a combined model with shell and beam elements is often used, for example, by Zeybek et al. [47], who dealt with the evaluation of stress resultants in the girders.
As the tank problem is a symmetrical thin-walled structure, it is not recommended to apply a linear buckling analysis. The resulting eigenvalues do not exactly predict the real global buckling and collapse of the tank. Therefore, the linear buckling analysis is preferably used as an input preliminary shape to the following more precise GMNIA, as is also recommended by Maraveas et al. [19]. In this analysis, the imperfections that are present in real structures are already implemented. Similarly, imperfections are included in the analytical equations for a simplified assessment of buckling, given in EN 1993-1-6 [11]. A sensitivity analysis of the imperfections can be performed, and the size of the imperfections modified. The imperfection amplitude can affect buckling capacity, which is also seen in the publication of the authors Chiang and Guzey [20].
The parametric study is performed for tank variants with uniform wall thickness t w = 3 mm. This is not very common in industrial applications, especially for larger tanks. The resulting values of the critical pressure are only used for a clearer comparison of the different types of girders and to realize a parametric study, so the thickness variation in each row is not considered. Three different tank diameters and three heights for each diameter are selected. The values in Figure 10, Figure 12, and Figure 13 show the non-linear trends. To smooth the curves, more values in the diameter interval D = 4.3–42 m should be used. Furthermore, an axisymmetric pressure is applied as a load over the entire circumference and height of the tank, which is conservative to buckling. However, when designing tanks according to real local conditions, the non-uniform distribution of wind pressure according to EN 1994-1-4 [8] is usually applied, where the effect of the negative pressure occurring in open-topped tanks should be included as well.
In Table 10, the main advantages of both wind girders are summarized. Based on the obtained results, NWG has more benefits compared to conventional profiles. Except mass per length is larger for NWG compared to L60, it may not mean that the total price is also higher.

5. Conclusions

This paper has presented the design and numerical analysis of a new wind girder (NWG) and the subsequent comparison of the designed shape with the model using the smallest used round-edged equal leg steel angle bar EN 10056-1 - L60 × 60 × 6 - EN 10025-2-S235JR (L60). The geometric design is based on a research study and on tank dimensions taken from technical applications. After designing the basic dimensions of the girder listed in Table 2, a parametric study of nine tank variants is carried out to compare buckling capacity and modal properties of tanks with conventional profile L60.
The shape of NWG (Figure 3a) is composed of segments with a rectangular cross-section that are bolted to the tank wall and to each other. This shape with a rectangular cross-section is easy to produce and also can be easily bent to the required curvature compared to the commonly used profiles. The designed girder can be quickly modified while retaining the geometry concept by only changing the thickness of the used plates. This increases their moment of inertia and thus buckling capacity of the whole tank. As listed in Table 11, the total mass of NWG is slightly higher than L60. However, due to the simpler cross-section of NWG, the production costs are lower. The main disadvantage of the new shape compared to the conventional girders is the variable moment of inertia around the circumference of the tank. This causes a different behavior, and it is not easy to design the dimensions of the girders using analytic approaches based on the comparison of the moment of inertia or the section modulus. The mentioned standards specify a minimum moment of inertia (or section modulus) for the upper and intermediate girders; therefore, the design of NWG should be verified by the non-linear numerical model. Based on the obtained results, NWG can replace L60; however, a more detailed numerical analysis is required for tank design to verify the application.
A non-linear analysis with imperfections (GMNIA) was used to compare the buckling capacity of girders. The imperfections were taken from the linear bifurcation analysis (LBA). The girders were analyzed on several tank variants of different diameters and heights. The results show that the designed shape of the wind girder can be used instead of the conventional girders. As the results obtained for the analyzed variants are relatively varied, it is advisable to assess the specific tank configuration with a given load using numerical analysis. For some variants, the new design gives a better result, while in other cases, the currently used girders appear to be better. According to the achieved results for the given configuration and chosen variants, NWG seems to be better for tanks with larger diameters. This may be due to the fact that for tanks with larger diameters more segments are used, whereas for small-diameter tanks, the number of segments is insufficient. Different punching and segment geometries could be used to achieve a larger number of V-segments in these cases. It is evident from the modal analysis that, mainly in the cases of low tanks, the first natural frequency for the tanks with the new girder design is lower than for the conventional profiles. In one variant (TK1), the first mode shape was also different and indicates that, in some cases, the tank is not adequately stiffened and a more rigid profile should be used. The lower first natural frequency may also have an influence on the formation of the resonance under wind or seismic action. Conclusions can be summarized in the following terms.
  • The cross-section of NWG segments is very simple and can be easily produced by common manufacturers.
  • Although the total mass of NWG is slightly higher than for L60, due to the obtained results, NWG could be compared with the larger cross-section profiles usually used for tanks with larger diameters and possibly replace them.
  • Although the total mass of NWG is slightly higher than for L60, the total price should be less because of the simpler cross-section and therefore lower material prices.
  • The geometry is designed so that NWG can be used for different tank curvatures and stiffness requirements. By modifying the basic dimensions of segments, the stiffness is increased and overall buckling capacity is also increased.
  • Because NWG is composed of simple segments with various curvatures, it is convenient for storage and assembly.
  • An important difference of NWG compared to conventional girders is the variable moment of inertia around the circumference of the tank, which requires a more detailed assessment using GMNIA.

Author Contributions

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

Funding

This work was supported by Specific Research “Application of Modern Computational and Experimental Approaches in Applied Mechanics” (SP2023/027).

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
3DThree-dimensional
FEMFinite element method
FSIFluid-structure interaction
GMNAGeometrically and materially non-linear analysis
GMNIAGeometrically and materially non-linear analysis with imperfections
GNAGeometrically non-linear elastic analysis
GNIAGeometrically non-linear elastic analysis with imperfections
K300TKOSMALT E 300 T
L60Round-edged equal leg steel angle bar
EN 10056-1 - L60 × 60 × 6 - EN 10025-2 - S235JR
LBALinear elastic bifurcation analysis
MAModal analysis
MNAMaterially non-linear analysis
NWGNew wind girder
P-segmentPlate segment of new wind girder
S235S235JR EN 10025-2
TKTANK
V-segmentV-shape segment of new wind grider

References

  1. Uematsu, Y.; Yamaguchi, T.; Yasunaga, J. Effects of wind girders on the buckling of open-topped storage tanks under quasi-static wind loading. Thin-Walled Struct. 2018, 124, 1–12. [Google Scholar] [CrossRef]
  2. Zhao, Y.; Lin, Y. Buckling of cylindrical open-topped steel tanks under wind load. Thin-Walled Struct. 2014, 79, 83–94. [Google Scholar] [CrossRef]
  3. Shokrzadeh, A.R.; Sohrabi, M.R. Buckling of ground based steel tanks subjected to wind and vacuum pressures considering uniform internal and external corrosion. Thin-Walled Struct. 2016, 108, 333–350. [Google Scholar] [CrossRef]
  4. Jaca, R.C.; Godoy, L.A. Wind buckling of metal tanks during their construction. Thin-Walled Struct. 2010, 48, 453–459. [Google Scholar] [CrossRef]
  5. EN 1993-4-2; Eurocode 3: Design of Steel Structures—Part 4-2: Tanks. The European Committee for Standardization: Brussels, Belgium, 2008.
  6. EN 1991-4; Eurocode 1: Actions on structures—Part 4: Silos and Tanks. The European Committee for Standardization: Brussels, Belgium, 2018.
  7. AWWA D103-09; Factory-Coated Bolted Carbon Steel Tanks for Water Storage. American Water Works Association: Denver, CO, USA, 2009.
  8. EN 1991-1-4; Eurocode 1: Actions on Structures—Part 1-4: General Actions—Wind Loads. The European Committee for Standardization: Brussels, Belgium, 2020.
  9. Burgos, C.A.; Batista-Abreu, J.C.; Calabró, H.D.; Jaca, R.C.; Godoy, L.A. Buckling estimates for oil storage tanks: Effect of simplified modeling of the roof and wind girder. Thin-Walled Struct. 2015, 91, 29–37. [Google Scholar] [CrossRef]
  10. Jawad, M. Theory and Design of Plate and Shell Structures; Springer Science & Business Media: New York, NY, USA, 2012. [Google Scholar]
  11. EN 1993-1-6; Eurocode 3: Design of Steel Structures—Part 1-6: Strength and Stability of Shell Structures. The European Committee for Standardization: Brussels, Belgium, 2008.
  12. Kminek, T.; Marsalek, P.; Karasek, T. Analysis of steel tanks for water storage using shell elements. AIP Conf. Proc. 2019, 2116, 320007. [Google Scholar]
  13. Elkholy, S.A.; Elsayed, A.A.; El-Ariss, B.; Sadek, S.A. Optimal finite element modelling for modal analysis of liquid storage circular tanks. Int. J. Struct. Eng. 2014, 5, 207–241. [Google Scholar] [CrossRef]
  14. Mazuch, T.; Horacek, J.; Trnka, J.; Veselỳ, J. Natural modes and frequencies of a thin clamped–free steel cylindrical storage tank partially filled with water: Fem and measurement. J. Sound Vib. 1996, 193, 669–690. [Google Scholar] [CrossRef]
  15. Pan, J.; Liang, S. A study on the buckling behavior of in-service large open-topped oil-storage tanks. Structures 2021, 29, 211–224. [Google Scholar] [CrossRef]
  16. Chen, L.; Rotter, J.M. Buckling of anchored cylindrical shells of uniform thickness under wind load. Eng. Struct. 2012, 41, 199–208. [Google Scholar] [CrossRef]
  17. Cao, Q.S.; Zhao, Y. Buckling design of large steel silos with various slendernesses. J. Zhejiang Univ. Sci. A 2017, 18, 282–305. [Google Scholar] [CrossRef]
  18. Zhao, Y.; Cao, Q.S.; Su, L. Buckling design of large circular steel silos subject to wind pressure. Thin-Walled Struct. 2013, 73, 337–349. [Google Scholar] [CrossRef]
  19. Maraveas, C.; Balokas, G.A.; Tsavdaridis, K.D. Numerical evaluation on shell buckling of empty thin-walled steel tanks under wind load according to current American and European design codes. Thin-Walled Struct. 2015, 95, 152–160. [Google Scholar] [CrossRef]
  20. Chiang, Y.C.; Guzey, S. Dynamic analysis of aboveground open-top steel tanks subjected to wind loading. Eng. Struct. 2019, 198, 109496. [Google Scholar] [CrossRef]
  21. Alam, M.M.; Zhou, Y. Strouhal numbers, forces and flow structures around two tandem cylinders of different diameters. J. Fluids Struct. 2008, 24, 505–526. [Google Scholar] [CrossRef]
  22. Ellingsen, Ø.M.; Amandolese, X.; Flamand, O.; Hémon, P. Twin Strouhal numbers in pressure loading of circular cylinder at high Reynolds numbers. J. Fluids Struct. 2022, 115, 103782. [Google Scholar] [CrossRef]
  23. Godoy, L.A. Buckling of vertical oil storage steel tanks: Review of static buckling studies. Thin-Walled Struct. 2016, 103, 1–21. [Google Scholar] [CrossRef]
  24. Burgos, C.A.; Jaca, R.C.; Lassig, J.L.; Godoy, L.A. Wind buckling of tanks with conical roof considering shielding by another tank. Thin-Walled Struct. 2014, 84, 226–240. [Google Scholar] [CrossRef]
  25. Pan, J.; Liang, S. Buckling analysis of open-topped steel tanks under external pressure. SN Appl. Sci. 2020, 2, 535. [Google Scholar] [CrossRef]
  26. ISO 28765:2016; Vitreous and Porcelain Enamels—Design of Bolted Steel Tanks for the Storage or Treatment of Water or Municipal or Industrial Effluents and Sludges. International Organization for Standardization: Geneva, Switzerland, 2016.
  27. Bu, F.; Qian, C. On the rational design of the top wind girder of large storage tanks. Thin-Walled Struct. 2016, 99, 91–96. [Google Scholar] [CrossRef]
  28. Bu, F.; Qian, C. A rational design approach of intermediate wind girders on large storage tanks. Thin-Walled Struct. 2015, 92, 76–81. [Google Scholar] [CrossRef]
  29. Lewandowski, M.J.; Gajewski, M.; Gizejowski, M. Numerical analysis of influence of intermediate stiffeners setting on the stability behaviour of thin-walled steel tank shell. Thin-Walled Struct. 2015, 90, 119–127. [Google Scholar] [CrossRef]
  30. Zeybek, Ö.; Topkaya, C. Stiffness requirements for wind girders in open-top cylindrical steel tanks. Thin-Walled Struct. 2022, 176, 109353. [Google Scholar] [CrossRef]
  31. Drahorad, L. Analysis of Steel Tank Supports for Liquid Storage. Master’s Thesis, VSB-Technical University of Ostrava, Ostrava, Czech Republic, 2021. [Google Scholar]
  32. EN 10056-1; Structural Steel Equal and Unequal Leg Angles—Part 1: Dimensions. The European Comittee for Standardization: Brussels, Belgium, 2017.
  33. EN 10025-2; Hot Rolled Products of Structural Steels—Part 2: Technology Delivery Conditions for Non-Alloy Structural Steels. The European Comittee for Standardization: Brussels, Belgium, 2020.
  34. Portela, G.; Godoy, L. Wind pressures and buckling of cylindrical steel tanks with a conical roof. J. Constr. Steel Res. 2005, 61, 786–807. [Google Scholar] [CrossRef]
  35. Horyl, P.; Snuparek, R.; Marsalek, P. Behaviour of frictional joints in steel arch yielding supports. Arch. Min. Sci. 2014, 59, 781–792. [Google Scholar] [CrossRef]
  36. Lesnak, M.; Marsalek, P.; Horyl, P.; Pistora, J. Load-Bearing Capacity Modelling and Testing of Single-Stranded Wire Rope. Acta Montan. Slovaca 2020, 25, 192–200. [Google Scholar]
  37. Rybansky, D.; Sotola, M.; Marsalek, P.; Poruba, Z.; Fusek, M. Study of Optimal Cam Design of Dual-Axle Spring-Loaded Camming Device. Materials 2021, 14, 1940. [Google Scholar] [CrossRef] [PubMed]
  38. Hroncek, J.; Marsalek, P.; Rybansky, D.; Sotola, M.; Drahorad, L.; Lesnak, M.; Fusek, M. Simplified Numerical Model for Determining Load-Bearing Capacity of Steel-Wire Ropes. Materials 2023, 16, 3756. [Google Scholar] [CrossRef]
  39. Rybansky, D.; Marsalek, P.; Sotola, M.; Hroncek, J.; Drahorad, L.; Kusnir, O.; Prokop, J. Design and Behavior of Lightweight Flexible Structure with Spatial Pattern Reducing Contact Surface Fraction. Polymers 2023, 15, 3896. [Google Scholar] [CrossRef]
  40. EN 1993-1-1; Eurocode 3: Design of Steel Structures—Part 1-1: General Rules and Rules for Buildings. The European Committee for Standardization: Brussels, Belgium, 2011.
  41. U. S. Steel Kosice, s.r.o. Steel for Enameling. 2023. Available online: https://www.usske.sk/sk/produkty/ocel-valcovana-za-tepla/ocel-pre-smaltovanie (accessed on 27 July 2023).
  42. The Scientific and Technical Center for Building (CSTB). Jules Verne Climatic Wind Tunnel. 2023. Available online: https://www.cstb.fr/en/test-facilities/jules-verne-climatic-wind-tunnel (accessed on 13 October 2023).
  43. Spritzer, J.; Guzey, S. Nonlinear numerical evaluation of large open-top aboveground steel welded liquid storage tanks excited by seismic loads. Thin-Walled Struct. 2017, 119, 662–676. [Google Scholar] [CrossRef]
  44. Compagnoni, M.E.; Curadelli, O. Experimental and numerical study of the response of cylindrical steel tanks under seismic excitation. Int. J. Civ. Eng. 2018, 16, 793–805. [Google Scholar] [CrossRef]
  45. Phan, H.N.; Paolacci, F.; Mongabure, P. Nonlinear finite element analysis of unanchored steel liquid storage tanks subjected to seismic loadings. In Proceedings of the Pressure Vessels and Piping Conference, Waikoloa, HI, USA, 16–20 July 2017; Volume 58035. [Google Scholar]
  46. Panchal, V.; Jangid, R. Seismic response of liquid storage steel tanks with variable frequency pendulum isolator. KSCE J. Civ. Eng. 2011, 15, 1041–1055. [Google Scholar] [CrossRef]
  47. Zeybek, Ö.; Topkaya, C.; Rotter, J.M. Stress resultants for wind girders in open-top cylindrical steel tanks. Eng. Struct. 2019, 196, 109347. [Google Scholar] [CrossRef]
Figure 1. Scheme of pressure distribution of wind around the tank circumference: (a) idealized; (b) real.
Figure 1. Scheme of pressure distribution of wind around the tank circumference: (a) idealized; (b) real.
Buildings 13 02724 g001
Figure 2. Types of conventional stiffeners: (a) vertical stiffeners; (b) horizontal girders; (c) detail A of horizontal girder.
Figure 2. Types of conventional stiffeners: (a) vertical stiffeners; (b) horizontal girders; (c) detail A of horizontal girder.
Buildings 13 02724 g002
Figure 3. Comparison between the girder shapes: (a) NWG; (b) conventional girder L60.
Figure 3. Comparison between the girder shapes: (a) NWG; (b) conventional girder L60.
Buildings 13 02724 g003
Figure 4. Shape of designed segments forming NWG: Item 1—V-segment; Item 2—P-segment; Item 3—bolted connection; Item 4—tank wall.
Figure 4. Shape of designed segments forming NWG: Item 1—V-segment; Item 2—P-segment; Item 3—bolted connection; Item 4—tank wall.
Buildings 13 02724 g004
Figure 5. Finite element meshes for variants with NWG and diameter D = 4.3 m: (a) low-height tank; (b) middle-height tank; (c) high tank.
Figure 5. Finite element meshes for variants with NWG and diameter D = 4.3 m: (a) low-height tank; (b) middle-height tank; (c) high tank.
Buildings 13 02724 g005
Figure 6. Idealized stress–strain curves for used material model of K300T (blue) and S235 (red).
Figure 6. Idealized stress–strain curves for used material model of K300T (blue) and S235 (red).
Buildings 13 02724 g006
Figure 7. Applied boundary conditions illustrated on the variant with NWG and diameter D = 4.3 m: (a) uniform pressure; and (b) fixed support.
Figure 7. Applied boundary conditions illustrated on the variant with NWG and diameter D = 4.3 m: (a) uniform pressure; and (b) fixed support.
Buildings 13 02724 g007
Figure 8. Tank TK1–shape under critical pressure from GMNIA: (a) L60; (b) NWG; first mode shape from MAl; (c) L60; (d) NWG.
Figure 8. Tank TK1–shape under critical pressure from GMNIA: (a) L60; (b) NWG; first mode shape from MAl; (c) L60; (d) NWG.
Buildings 13 02724 g008
Figure 9. Tank TK3–shape under critical pressure from GMNIA: (a) L60; (b) NWG; first mode shape from MA; (c) L60; (d) NWG.
Figure 9. Tank TK3–shape under critical pressure from GMNIA: (a) L60; (b) NWG; first mode shape from MA; (c) L60; (d) NWG.
Buildings 13 02724 g009
Figure 10. Comparison of low-height tanks: on the left is critical pressure, and on the right is first natural frequency.
Figure 10. Comparison of low-height tanks: on the left is critical pressure, and on the right is first natural frequency.
Buildings 13 02724 g010
Figure 11. Tank TK5—shape under critical pressure from GMNIA: (a) L60; (b) NWG; first mode shape from MA; (c) L60; (d) NWG.
Figure 11. Tank TK5—shape under critical pressure from GMNIA: (a) L60; (b) NWG; first mode shape from MA; (c) L60; (d) NWG.
Buildings 13 02724 g011
Figure 12. Comparison of middle-height tanks: on the left is critical pressure, and on the right is first natural frequency.
Figure 12. Comparison of middle-height tanks: on the left is critical pressure, and on the right is first natural frequency.
Buildings 13 02724 g012
Figure 13. Comparison of high tanks: on the left is critical pressure, and on the right is first natural frequency.
Figure 13. Comparison of high tanks: on the left is critical pressure, and on the right is first natural frequency.
Buildings 13 02724 g013
Figure 14. Tank TK7—shape under critical pressure from GMNIA: (a) L60; (b) NWG; first mode shape from MA; (c) L60; (d) NWG.
Figure 14. Tank TK7—shape under critical pressure from GMNIA: (a) L60; (b) NWG; first mode shape from MA; (c) L60; (d) NWG.
Buildings 13 02724 g014
Figure 15. Dependence of the critical pressure on imperfection amplitude.
Figure 15. Dependence of the critical pressure on imperfection amplitude.
Buildings 13 02724 g015
Figure 16. Load–deflection curves for variant TK7 and different sizes of imperfection.
Figure 16. Load–deflection curves for variant TK7 and different sizes of imperfection.
Buildings 13 02724 g016
Figure 17. Load–deflection curves for variant (a) TK4 and (b) TK6.
Figure 17. Load–deflection curves for variant (a) TK4 and (b) TK6.
Buildings 13 02724 g017
Table 1. Shell theories applicable for numerical simulation of the tanks.
Table 1. Shell theories applicable for numerical simulation of the tanks.
Type of Tank AnalysisShell Theory
Linear elastic shell analysis (LA)Linear bending and stretching
Linear elastic bifurcation analysis (LBA)Linear bending and stretching
Geometrically non-linear elastic analysis (GNA)Non-linear
Materially non-linear analysis (MNA)Linear
Geometrically and materially non-linear analysis (GMNA)Non-linear
Geometrically non-linear elastic analysis with imperfections (GNIA)Non-linear
Geometrically and materially non-linear analysis with imperfections (GMNIA)Non-linear
Table 2. Designed dimensions of V-segment and P-segment forming NWG.
Table 2. Designed dimensions of V-segment and P-segment forming NWG.
Thickness tHeight hWidth w
4 mm150 mm50 mm
Table 3. Material constants used for the tank wall and girders.
Table 3. Material constants used for the tank wall and girders.
Tank WallGirders
Steel gradeK300TS235
Young’s modulus E210,000 MPa210,000 MPa
Poisson’s ratio μ 0.300.30
Yield stress σ y 300 MPa235 MPa
Elongation δ 28%26%
Ultimate stress σ u l t 360 MPa235 MPa
Tangent modulus E T 531 MPa844 MPa
Table 4. Geometry of tank variants selected for the parametric study.
Table 4. Geometry of tank variants selected for the parametric study.
Variant LabelD (m)H (m) t w (m)Volume (m3)Number of Nodes (-)Number of Elements (-)Number of Segments (-)
TK14.32.90.0034228,90026,80070
TK221.42.90.0031043143,900133,400350
TK342.02.90.0034018208,900189,100686
TK44.38.60.00312578,70076,60070
TK521.48.60.0033093392,600382,500350
TK642.08.60.00311,915549,900530,100686
TK74.317.20.003250153,300151,30070
TK821.417.20.0036187417,400407,300350
TK942.017.20.00323,830598,700578,800686
Table 5. Comparison of values obtained for low-height tanks.
Table 5. Comparison of values obtained for low-height tanks.
VariantD p L 60 f L 60 p NWG f NWG p NWG p L 60 p L 60 f NWG f L 60 f L 60 σ L 60 , w σ NWG , w σ L 60 , G σ NWG , G
Label (m) (kPa) (Hz) (kPa) (Hz) (%) (%) (MPa) (MPa) (MPa) (MPa)
TK14.322.50039.11116.85028.035−25.111−28.321671391562
TK221.41.75016.6401.50016.038−14.286−3.6141631617
TK342.00.64011.1870.65012.0161.5637.4041819514
Table 6. Comparison of values obtained for middle-height tanks.
Table 6. Comparison of values obtained for middle-height tanks.
VariantD p L 60 f L 60 p NWG f NWG p NWG p L 60 p L 60 f NWG f L 60 f L 60 σ L 60 , w σ NWG , w σ L 60 , G σ NWG , G
Label (m) (kPa) (Hz) (kPa) (Hz) (%) (%) (MPa) (MPa) (MPa) (MPa)
TK44.38.70010.8747.35010.407−15.517−4.295961281294
TK521.40.5855.3620.5946.4301.49619.9092934517
TK642.00.2003.7620.2054.9052.50030.382161738
Table 7. Comparison of values obtained for high tanks.
Table 7. Comparison of values obtained for high tanks.
VariantD p L 60 f L 60 p NWG f NWG p NWG p L 60 p L 60 f NWG f L 60 f L 60 σ L 60 , w σ NWG , w σ L 60 , G σ NWG , G
Label (m) (kPa) (Hz) (kPa) (Hz) (%) (%) (MPa) (MPa) (MPa) (MPa)
TK74.34.3754.8413.8835.116−11.2495.67911412219119
TK821.40.3602.4260.3333.317−7.63936.7092836519
TK942.00.1121.7210.1142.5691.78649.288141538
Table 8. Values of critical pressure obtained from LBA, GMNIA, and analytical models.
Table 8. Values of critical pressure obtained from LBA, GMNIA, and analytical models.
VariantD p LBA , L 60 p GMNIA , L 60 p G M N I A , L 60 p L B A , L 60 p LBA , NWG p GMNIA , NWG p G M N I A , N W G p L B A , N W G p Eurocode p ISO
Label (m) (kPa) (kPa) (kPa) (kPa) (kPa) (kPa)
TK14.325.83522.5000.87125.57216.8500.65914.40419.717
TK221.42.6131.7500.6702.6371.5000.5691.2971.776
TK342.01.0650.6400.6011.0820.6500.6010.4710.646
TK44.39.1108.7000.9558.9737.3500.8194.7906.557
TK521.40.8390.5850.6970.8450.5940.7030.4310.591
TK642.00.3120.2000.6400.3170.2050.6470.1560.215
TK74.34.4704.3750.9794.3743.8830.8882.3953.279
TK821.40.4080.3600.8830.4090.3330.8120.2150.295
TK942.00.1500.1120.7480.1510.1140.7540.0780.107
Table 9. Critical pressures for variants TK1, TK4, and TK7 for different imperfection amplitudes.
Table 9. Critical pressures for variants TK1, TK4, and TK7 for different imperfection amplitudes.
TK1TK4TK7
u i P L 60  (kPa) P NWG  (kPa) P L 60  (kPa) P NWG  (kPa) P L 60  (kPa) P NWG  (kPa)
0.25 · t w = 0.75 mm23.50017.5008.8007.9704.4004.076
0.50 · t w = 1.50 mm22.50016.8508.7007.3504.3753.884
1.00 · t w = 3.00 mm21.08816.5008.6006.4484.3593.609
Table 10. Comparison of conventional girders and NWG.
Table 10. Comparison of conventional girders and NWG.
Conventional GirdersNew Design of Girders
In some cases difficult to produceEasier production
Difficult adaptation to the tankEasy adaptation to various tank sizes
Need to select from available standardized profilesEasy to change stiffness
Uniform cross-sectionVariable cross-section
Lower massLarger mass (18% more compared to L60)
Table 11. Comparison of mass per length of girders.
Table 11. Comparison of mass per length of girders.
Part of GirderL60NWG
Mass per length of girders (kg/m)6.397.55
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

Drahorad, L.; Marsalek, P.; Hroncek, J.; Rybansky, D.; Sotola, M.; Poruba, Z.; Larys, M. Design and Non-Linear Modeling of New Wind Girder Used for Bolted Tanks. Buildings 2023, 13, 2724. https://doi.org/10.3390/buildings13112724

AMA Style

Drahorad L, Marsalek P, Hroncek J, Rybansky D, Sotola M, Poruba Z, Larys M. Design and Non-Linear Modeling of New Wind Girder Used for Bolted Tanks. Buildings. 2023; 13(11):2724. https://doi.org/10.3390/buildings13112724

Chicago/Turabian Style

Drahorad, Lukas, Pavel Marsalek, Juraj Hroncek, David Rybansky, Martin Sotola, Zdenek Poruba, and Michal Larys. 2023. "Design and Non-Linear Modeling of New Wind Girder Used for Bolted Tanks" Buildings 13, no. 11: 2724. https://doi.org/10.3390/buildings13112724

APA Style

Drahorad, L., Marsalek, P., Hroncek, J., Rybansky, D., Sotola, M., Poruba, Z., & Larys, M. (2023). Design and Non-Linear Modeling of New Wind Girder Used for Bolted Tanks. Buildings, 13(11), 2724. https://doi.org/10.3390/buildings13112724

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