Next Article in Journal
Total Polyphenols and Antioxidant Properties of Selected Fresh and Dried Herbs and Spices
Previous Article in Journal
Active Jet Noise Control of Turbofan Engine Based on Explicit Model Predictive Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Size and Shape Optimization of a Guyed Mast Structure under Wind, Ice and Seismic Loading

by
Raffaele Cucuzza
1,
Marco Martino Rosso
1,
Angelo Aloisio
2,*,
Jonathan Melchiorre
1,
Mario Lo Giudice
1 and
Giuseppe Carlo Marano
1
1
Department of Structural, Geotechnical and Building Engineering, Politecnico di Torino, Corso Duca Degli Abruzzi, 24, 10128 Turin, Italy
2
Civil Environmental and Architectural Engineering Department, Università Degli Studi dell’Aquila, via Giovanni Gronchi n.18, 67100 L’Aquila, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(10), 4875; https://doi.org/10.3390/app12104875
Submission received: 17 April 2022 / Revised: 8 May 2022 / Accepted: 9 May 2022 / Published: 11 May 2022
(This article belongs to the Topic Artificial Intelligence (AI) Applied in Civil Engineering)

Abstract

:
This paper discusses the size and shape optimization of a guyed radio mast for radiocommunications. The considered structure represents a widely industrial solution due to the recent spread of 5G and 6G mobile networks. The guyed radio mast was modeled with the finite element software SAP2000 and optimized through a genetic optimization algorithm (GA). The optimization exploits the open application programming interfaces (OAPI) SAP2000-Matlab. Static and dynamic analyses were carried out to provide realistic design scenarios of the mast structure. The authors considered the action of wind, ice, and seismic loads as variable loads. A parametric study on the most critical design variables includes several optimization scenarios to minimize the structure’s total self-weight by varying the most relevant parameters selected by a preliminary sensitivity analysis. In conclusion, final design considerations are discussed by highlighting the best optimization scenario in terms of the objective function and the number of parameters involved in the analysis.

1. Introduction

Guyed masts are extensively used in the telecommunications industry, and the size, shape, and topology optimization can significantly benefit their transportation and installation. The main loads acting on guyed mast structures arise from wind [1,2], earthquakes [3,4,5,6], sudden rupture of guys [7], galloping of guys [8], and sudden ice shedding from ice-laden guy wires [9].
Their optimization must fulfil several requirements under ultimate and service limit states [10]. Specifically, service limit states are crucial for guyed mast structures due to high-amplitude oscillations caused by their high deformability. In some cases, these vibrations have led to a signal loss caused by excessive displacement and rotation of the antennas and, in other cases, have resulted in permanent deformation or failure. Therefore, size optimization of the guyed mast structure represents a challenging task since the increment of the performance ratio of the materials should be counterbalanced by an adequate lateral stiffness to reduce high-vibration drawbacks [11].
Saxena [12] reported several happenings where heavy icing combined with moderate wind resulted in severe misalignment of towers and complete failure. Novak et al. [13] showed that ice accumulation on some parts of the guy wires and moderate winds could lead to the guy galloping, resulting in unacceptable stress levels throughout the structure. The main topics investigated in the field of guyed structures can be summarized as follows:
  • Structural design. Several researchers investigated the dynamic response of guyed mast structures through experimental tests and numerical modeling to derive design approaches and recommendations [14,15,16]. In particular, there are studies dealing with the dynamic identification and accurate estimate of the wind loads [17,18,19,20,21].
  • Nonlinear dynamics. The proneness to global and local instabilities challenged several scholars to estimate and predict the nonlinear behaviour of guyed masts [22,23,24,25,26].
  • Structural optimization. The need for guyed structures that are easy to install and transport challenged several scholars to optimize their shape in order to reduce the structural mass without reducing the lateral stiffness and prevent instability phenomena [27].
  • Structural control. There are some attempts of control methods to reduce vibrations in mast-like structures [28,29,30]. Among others, Blachowski [31] proposed the use of a hydraulic actuator to control cable forces in guyed masts using Kalman filtering.
This paper tackles the size and shape optimization of guyed mast structures. A video of the considered structure is available in Supplementary Material. Since the first attempts by Bell and Brown [32], many engineers attempted to optimize guyed masts under wind loads using deterministic global optimization algorithms. However, as pointed out by [27], this approach leads to local optimum points, since each design variable was considered separately. Thornton et al. [33] and Uys et al. [34] proposed general procedures for optimizing steel towers under dynamic loads. To the author’s knowledge, Venanzi and Materazzi [35] were the first to implement a multi-objective optimization method for guyed mast structures under wind loads using the stochastic simulated annealing algorithm for size optimization. The objective function implemented by [35] included the sum of the squares of the nodal displacements and the in-plan width of the structure. Zhang and Li [36] attempted to achieve both shape and size optimization in a two-step procedure using the ant colony algorithm (ACA). Cucuzza et al. [37] proposed an alternative approach in which the multi-objective optimization problem has been reduced to a single-objective problem through suitable parameters. Luh and Lin [38] were challenged in achieving the topology, size, and shape optimization of guyed masts using a modification of the binary particle swarm optimization (PSO) and the attractive and repulsive particle swarm optimization.
This paper discusses the size optimization of guyed masts using a genetic algorithm (GA) by considering different design scenarios (e.g., Cucuzza et al. [37] and Manuello et al. [39]). Kaveh and Talatahari [40] noticed that the particle swarm optimization (PSO) is more effective than ACA and the harmony search scheme for optimizing truss structures. However, Deng et al. [41] and Guo and Li [42] were successful in optimizing tapered masts and transmission towers using modifications of genetic algorithms (GA). Moreover, Belevivcius et al. [27] attempted the topology-sizing optimization problem of the guyed mast as a single-level single-objective global optimization problem using GAs.
Therefore, given the numerous successful solutions of guyed masts using GAs, the authors chose to investigate the size optimization of a guyed mast structure using GAs. Following [35], this paper focuses on the size optimization by considering eight possible design scenarios. The purpose of the present paper is two-fold. Firstly, this work aims at achieving a size optimization on a real application case adopting structural verification according to Eurocode 3. During the load evaluation phase, detailed analyses have been conducted, including wind, ice, and seismic actions and the verifications against instabilities. Secondly, the computational intelligence procedure adopted by the authors allowed the investigation of several scenarios simultaneously. As a result, the parameters that mainly affected the design process have been detected to provide preliminary indications to engineers in the practical design of similar structural typologies. Furthermore, the considered case study may represent a benchmark case for validating the reliability and accuracy of alternative numerical approaches. Therefore, the paper is organized as follows. After the case study description and the FE model, the authors introduce the first numerical results and the outcomes of the size and shape optimization.

2. Case Study

The considered structure is a guyed radio mast. It is a thin, slender, vertical structure sustained by tension cables fixed to the ground and typically arranged at 120° between each other.
The main body is a single central column made of tube profiles or truss systems when a high elevation must be reached, see Figure 1. More than one set of cables is placed at different elevations to prevent instability phenomena. Guyed towers are usually built for meteorological purposes or to support radio antennas, such as the one considered in this research. In particular, this structure can be used for a limited time during an event or maintenance of primary transmission towers. Therefore, it is also called a temporary base transceiver station (BTS), typically adopted to supply the immediate service. Sporting events, concerts, motor racing, military camps, and emergency events are typical examples of temporary BTS applications. The BTS is usually mounted on a moveable platform called the shelter.
The considered structure is located in Bassano Del Grappa, in the north of Italy, at a 129 m elevation from the sea level. The surrounding area is low-urbanized, with no relevant obstacles to the wind loads. The total height of the mast is 30.00 m. It is sustained by a central pole where 21 cables are fixed, see Figure 2. Other structural elements with rectangular cross-sections are used to create truss systems connecting cables and the central pole.
The central pole consists of five circular hollow steel profiles with flanged joints and 6 m in length. All connections are bolted, as well as those connecting the cables to the pole. The shelter is a steel box devoted to partially sustaining the structure and hosting electronic equipment. It is usually mounted on a moveable platform.

3. Load Analysis

This section details the loads acting on the structures, from the dead to the variable loads. According to the Italian Standard Regulation NTC2018, the load combinations of the actions have been evaluated at the ultimate limit state (ULS) and, for seismic conditions, at the life safety (LS) limit state. In Appendix A, Table A4 illustrates the most critical combinations for both static and dynamic configurations. Partial safety factors γ and combination coefficients ψ were adopted in order to consider maximization (positive sign) or minimization (negative sign) of effects both for vertical and horizontal actions.

3.1. Dead Loads

The structure is made of steel S355, whose mechanical stress-strain behaviour is depicted in Figure 3, and the characteristics are listed here: f u s = 510 MPa, f y s = 355 MPa, E s = 210,000 MPa, which are the ultimate and yielding stresses and Young’s modulus, respectively.
The cables are made of galvanized steel consisting of 6 strands (216 wires) with an independent metal core (49 wires). The main characteristics are illustrated in Table 1.
The structure investigated consists of a few types of elements, as indicated in Table 2. Dead loads are calculated from the weight per unit length of each member.
The non-structural dead loads originate from the wiring weight and the steel ladder for inspection and maintenance. This load results in 0.3 kN/m. Antennas and parabolas represent the weight of the equipment. Two groups of three antennas are located at 26.00 and 29.25 m in height, with a 120° in mutual spacing. The first one is the model AOC4518R7v06 produced by Huawei®. The second one is the model 6888670N manufactured by Amphenol®. Finally, there are three parabolas located at 23.15 m height, spaced 120° apart from each other, 30 cm in diameter. Table 3 and Table 4 detail the weight of the equipment and the non-structural dead loads.

3.2. Variable Loads

In this section, the detailed load modeling phase, for each variable load considered, is described. With specific reference to the wind action evaluation, the drag and lift forces are calculated according to the CNR-DT 207 R1/2018 [43]. The relationship between inertia and viscous forces, i.e., how wind load impacts to the surface, is taken into account with the Reynold’s number R e with the following expression:
R e ( z ) = l · v m ( z ) ν
where z is the elevation, l is the characteristic length, v m is the averaged wind speed, while ν is the kinematic viscosity of air ( ν = 15 × 10 6 m 2 /s).

3.2.1. Maintenance and Repairing Loads

Following the Italian national recommendations [44], it is supposed that a typical situation of inspection or maintenance is performed by an operator working on the steel ladder. A concentrated load of 120 kg is applied at the top of the tower. Despite that, it is reasonable to believe that the operator could work by using a basket elevator, without loading the structure.

3.2.2. Wind Loads

The wind action was evaluated according to the Italian recommendations in [43]. Firstly, the peak kinetic pressure ( q p ) was evaluated as follows:
q p = 1 2 · ρ · v r 2 · c e ( z )
where p is the kinetic pressure, while:
  • ρ is the air density;
  • v r 2 is the reference wind velocity;
  • c e is the exposure coefficient, varying with the elevation z of the structure.
For this purpose, the equivalent longitudinal or drag forces, f D , and transverse or lift force, f L , are evaluated as follows:
f d r a g = q p ( z ) · l · c d r a g ; f l i f t = q p ( z ) · b · c l i f t
where
  • q p ( z ) is the peak kinetic pressure evaluated at height z;
  • l is the characteristic element size;
  • b is the reference transverse dimension of the section;
  • c d r a g and c l i f t are the longitudinal and transverse dynamic coefficients.
Drag D and Lift L forces are reported in Table A2 and Table A3.

3.2.3. Ice Load

Ice and snow attached to the structural surface can significantly increase the variable loads in flexible and light structures. In particular, the radio mast is very sensitive to changes in the wind-exposed surface. In addition, the ice covering can increase the volume and the surface of the structural elements more than twice due to the low thickness of the central pole. The recommendations in [43] provide several scenarios for ice coverings. In the absence of more detailed evaluations, it is customary to consider an ice sleeve formation that is 12.5 mm thick. After the estimate of the wind loads, the influence of the ice sleeve formation on the structure is considered by assuming an additional exposed surface equal to 15% of the original one.

3.2.4. Seismic Action

Seismic action is evaluated according to the Italian seismic hazard map [44]. A linear dynamic analysis with seismic elastic response spectrum corresponding to the service limit state was carried out. Specifically, seismic actions are considered as acting independently in the X and Y plane directions.
The elastic response spectrum considered in the analysis was calculated by considering the topographic category of the site and geometry of the building (Figure 4). The first 33 vibration modes of the structure are included in the analysis, to reach 85 % of the total participating mass according to the national regulations in [44]. The mass participating ratios are listed in Appendix A.

4. Finite Element Modeling

The structural model was developed using two different element types: beams and cables. Beam elements model the main pole and all structural elements except for the cables. They possess the geometric and material properties of structural elements. The beam elements are used to model the main pole and secondary elements. Moreover, except for the main pole, rotation releases are applied at the ends in order to consider no flexural rigidity, as occurring for trussed structures.
Cable elements are used to simulate the steel ropes. Cable elements undergo large displacements that give rise to geometric nonlinearities. Therefore, the equilibrium of the cables is considered in the deformed configuration using SAP2000. As a result, the structural behaviour of guyed towers can be highly nonlinear, especially for low pre-tension cables, which are prone to large displacements. On the contrary, the nonlinear behaviour becomes less pronounced by increasing the pre-tension, resulting in high compression levels and minor instability effects. This paper considers the envelope of the maximum and minimum responses associated with each load condition.
Figure 5 plots the three modes with a higher mass participation ratio. These are the 10th, 11th, and 12th modes obtained from the dynamic analysis of the mast structure with the dead loads. On the contrary, the first modes arising from the dynamic analysis have lower mass participation factors and are characterized by local deformation of the structural elements. The 10th, 11th, and 12th modes are the first modes exhibiting the global deformation of the mast structure.
X and Y identify the in-plane orthogonal directions. The 10th and 11th modes have an approximate 26% mass participation factor in the Y and X directions, respectively. The natural period is very low and at approximately 0.4 s. The 13th mode is mainly torsional with nearly a 7 and 4% mass participation factor in the X and Y directions.
Figure 6 shows the positive (in dark green and purple) and negative (in red and light green) maximum and minimum envelopes of the axial, shear forces, and bending moments acting on the structural elements. Figure 7 plots the performance ratios of all structural elements except for the cables. The performance ratio is the ratio between the maximum stress in the structural element and the yielding stress. The performance ratios are defined by the colour map next to Figure 7. The plots highlight the presence of a structural element in the first half of the central pole with a high-performance ratio, depicted in yellow. The first section of the central pole has a low performance ratio, lower than 0.25. After the section with a performance ratio in the range 0.4–0.65, the following sections fall in the range 0.25–0.4 and are coloured in green. The top sections of the central pole are not significantly stressed, with a performance ratio of 0–0.25. The bracings have low stress, plotted in cyan, with performance ratios of 0–0.25.
Figure 8 shows the maximum displacements in the X ( u 1 ), Y ( u 2 ) directions and their combination ( u t ) at the service limit state. The maximum displacement is located at the top of the tower, in particular at joint 6 (z = 30.00 m), with a maximum displacement equal to u t = 18.7 mm.

5. Structural Optimization

In optimization problems, the main goal is to find the best conditions in terms of the optimal set of design parameters collected in the design vector x , which minimizes an objective function (OF) f ( x ) [45,46,47]. These problems can be categorized into single-objective or multi-objective based on the number of OFs involved, and a further classification is based on the presence (or not) of constraints [48,49,50]. In the structural optimization field, it is common to deal with constrained optimization, whose general statement is [51]:
min x Ω { f ( x ) } s . t . g q ( x ) 0 q = 1 , , n q h r ( x ) = 0 r = 1 , , n r
where x = { x 1 , , x j , , x n } T is the design vector to be optimized, whose terms are limited into a hyper-rectangular multidimensional box-type search space domain of interest denoted as Ω , given by the Cartesian product of the range of interest of each j-th of each design variable bounded in [ x j l , x j u ] , Ω = [ x 1 l , x 1 u ] × × [ x j l , x j u ] × × [ x n l , x n u ] . The term g q ( x ) in (4) denotes inequality constraints whereas h r ( x ) are equality ones, which further reduce the feasible search space inside Ω . In structural optimization, it is typical to deal with inequality constraints, and a common goal is to minimize the global cost of the structure. Since this involves many terms, the main attempt is minimizing the self-weight of the structure, indirectly connected to material cost, i.e., material usage and natural resources consumption [51]. Several strategies have been developed over the years to handle constraints [52,53,54]. In the present work, the penalty function-based approach was implemented due to its simplicity, allowing converting the problem with OF f ( x ) into a new unconstrained version ϕ ( x ) :
min x Ω { ϕ ( x ) ) } = min x Ω { f ( x ) + H ( x ) }
where H ( x ) is the penalty function. Adopting a static-penalty strategy, H ( x ) , assume this form [55,56]
H s ( x ) = w 1 H N V C ( x ) + w 2 H S V C ( x )
where H N V C is the number of violated constraints and H S V C is the sum of all violations:
H S V C ( x ) = p = 1 n p max { 0 , g p ( x ) }
w 1 and w 2 are the violation control parameters, whose numerical values are assumed equal to w 1 = w 2 = 100 following [55].
In the current study, the authors carried out a parametric study on the design variables of the guyed mast. This fact has led to eight different scenarios, summarized in Table 5. In addition, the starting initial values of the design parameter are listed in Table 6, while the general optimization workflow is illustrated in Figure 9. To compare the results, the focus is related only to the performance ratios PR of the central pole of the guyed radio mast, being the pole the most stressed element. It consists of five segments 6.00 m long with the same cross-section. Thus, starting from the ground level:
  • Pole 1 (0.00 to 6.00 m);
  • Pole 2 (6.00 to 12.00 m);
  • Pole 3 (12.00 to 18.00 m);
  • Pole 4 (18.00 to 24.00 m);
  • Pole 5 (24.00 to 30.00 m).
Starting with a constant diameter of the cross-section for the pole, at the end of the optimization, it is advisable to find a tapered solution following a linear relationship with the height, as represented in Figure 10f. Accordingly, it is possible to shape the pole cross-section with two design variables described by the bottom Φ i and top Φ f diameters. In the following, the different scenarios obtained from the parametric study based on the design variables involved in the optimization problem are described:
  • Scenario A: this scenario involves the diameter Φ , as a sole variable, in the attempt to reduce the material consumption with a constant pole cross-section diameter with the height, as illustrated in Figure 10a.
  • Scenario B: this scenario attempts to refine the previous case by adopting a tapered solution for the pole, by using the bottom Φ i and the top Φ f diameters, as represented in Figure 10b.
  • Scenario C: further improvements are considered concerning scenario B by adding the cable pre-stressing force F as a variable of the optimization, as represented in Figure 10c.
  • Scenario D: further improvements are considered to scenario B by using a unique value for the pole thickness t of the tapered elements of the pole, as represented in Figure 10c.
  • Scenario E: further improvements are considered with respect to scenario B by optimizing both cable pre-stressing force F with a unique value of thickness t for the tapered elements of the pole, as represented in Figure 10e.
  • Scenario F: from the structural analysis of scenario E, it is possible to point out how the linear law for the tapering forces to use a larger section where it is not necessary. Elements 2 and 3 are the most stressed ones. Therefore it is possible to further refine scenario E by considering a thickness value for the intermediate pole elements t i n t e r and a different thickness for the other extremal pole elements t e n d s .
  • Scenario G: in this scenario, the five different thickness values only have been governed for every pole element { t 1 , t 2 , t 3 , t 4 , t 5 } for a constant diameter solution with height, as depicted in Figure 10f.
  • Scenario H: in this last scenario, a complete approach involves both the tapered solution by governing the initial bottom Φ i and the final top Φ f diameters, the five values of thickness for every pole element { t 1 , t 2 , t 3 , t 4 , t 5 } , and even the cable pre-stressing force.

Constraints Involved in the Structural Optimization Problem

The optimization problem statement is reported in (4) and the constraints were treated with the penalty-based approach illustrated in (5), by converting the constrained problem into an equivalent unconstrained one. The resolution of the optimization task considers the structural design assessment required by national and international codes to ensure the safety of constructions. In particular, the structural verifications derive from Eurocode 3 (EN 1993-1-1: 2005) and are referred to the ultimate limit state (ULS). The design verifications include tensile, compression, and buckling verification, and a combined assessment, such as the interaction capacity according to Annex B of the Eurocode 3:
D C = N E d χ y A f y k γ M 1 + k y y M y , E d χ L T W p l , y f y k γ M 1 + k y z M z , E d W p l , z f y k γ M 1 1
D C = N E d χ z A f y k γ M 1 + k z y M y , E d χ L T W p l , y f y k γ M 1 + k z z M z , E d W p l , z f y k γ M 1 1
where D stands for the demand and C stands for the capacity of the structure. Specifically, N E d is the acting axial force, whereas M y , E d and M z , E d represent the acting bending moments in the two principal directions of a planar local reference system centered on the cross section center of gravity. A is the cross section area of the pole, W p l , y and W p l , z are the plastic section modulus in the two principal directions, f y k is the yielding strength of the steel, whereas γ M 1 is the partial safety factor for instability conditions, equal to 1.05 from the Italian National Annex. χ L T is the reduction factor for lateral–torsional buckling, whereas k y y , k y z , k z y , and k z z are interaction factors whose values are derived according to two alternative approaches based on Annex A (method 1)and Annex B (method 2). The global structural deformation referred to the service limit state (SLS) has also been considered by verifying the top displacement of the mast. Specific recommendations for guyed mast structures are missing in national and international codes. Therefore, the authors adopted the suggestions defined in the Italian Technical Code NTC2018 (D.M.17/01/2018) reported in Chapter 4.2.4.2.2 Table 4.2.XIII related to limitations of lateral displacements of steel multi-storey frame structures. These limitations express a threshold condition in terms of the total height of the structure H:
δ SLS , top δ SLS , top , lim = H 500 = 30000 mm 500 = 60 mm
Since this condition is specific for steel multi-storey frame structures, the authors will assume this value as a reasonable choice to ensure service life assessment and preservation of working conditions of the telecommunication guyed mast tower. In the next section, a discussion on the results is carried out.

6. Results and Discussion

The paper compares the outcomes of the size and shape optimization in eight different scenarios, distinguished by different design variables. Scenario A is associated with the worst improvement of the structural performance since a single diameter is used for the central pole. Additionally, industrial steel profiles do not cover all possible ranges of the diameter. Improvements in the structural performance and weight reduction are achieved in the following scenarios when the search space becomes larger by increasing the number of design variables.
Scenario B introduces the tapering of the central pole with a linear variation from the bottom to the top. In this case, the optimal solution is affected by intermediate sections, which are more stressed. Consequently, the end cross-sections are over-estimated. In response to that, Scenario F introduces the linear tapering of the tube thickness t e n d s , t i n t e r to enhance the performance of the optimal solution. Parallelly, in Scenario G, five different thicknesses are adopted ( t 1 , t 2 , t 3 , t 4 , t 5 ), and the results are analogue to case F. Therefore, the thickness of the steel members is a suitable optimization parameter. At the same time, the diameter alone is not capable of returning attractive solutions because a linear interpolation trend is used. In addition, lower and upper limits were imposed for d and t. In particular, for this kind of structure, a minimum diameter d m i n 100 mm and a minimum thickness t m i n 3 mm was imposed.
The cross-section area depends on the square of the thickness. Therefore, small changes in t significantly affect the resulting area. Conversely, if the diameter is the sole search space, despite being tapered linearly with height, even significant modifications may not produce notable improvements. Still, the increment of design variables involved in the structural optimization typically increases the computational efforts. However, the scenario with the highest number of variables was characterized by an average time iteration close to 18s, using a computer with average performance. The computational effort cost of the optimization procedure strongly depends on the machine performance, no convergence issues occur. Table 7 lists the average values of performance ratio obtained from the eight optimization scenarios. All scenarios were collected in terms of number of parameters involved during the analysis. Table 7 proves that the increment in the number of design variables is associated with higher performance ratios. The target of the optimization achieves the best weight reduction, fully exploiting the structural material, without exceeding the ultimate and service limit states. Table 7 lists three sets of performance ratios: the initial one before optimization, the optimized, and the one obtained using commercial steel profiles, called the design performance ratio. The averaged performance ratio is equal to 28% before optimization. It significantly increases from scenario A, nearly 45%, to scenario G with 68%.
Essentially related to PR, mass reduction gives an idea about how much lighter (or heavier) the structure becomes due to the optimization process. It directly provides an estimate of cost savings.
Therefore, the results in Table 8 are consistent with the ones in terms of performance ratios, shown in Table 7.
Figure 11 shows the optimization results for the Scenario G, in term of the performance ratio obtained by averaging the performance ratios for each structural element. The results for all scenarios are reported in Appendix A. Scenario G, depicted below, exhibits higher values of the performance ratios. This fact becomes become more evident for poles 2, 3, and 4. In these cases, the performance ratios, associated with the design solutions, achieved values equal or greater than the optimized one due to the approximation of the design section adopted. In the post-processing phase, in fact, the optimized section chosen by the list of the FE software was manually edited since the structural constraint violation or the maximum performance ratio was not reached during to the optimization process. Moreover, in Table 9, the optimized design section for different independent iterations and the proposed industrial solutions according to product list, provided by the software, are listed. As expected, the mass reduction achieved during the optimization process results higher than the design solution due to the approximation issue. For the proposed scenario, the iteration ( N t r i a l ) that guarantees the best objective function is the second. In Appendix A, the graphical (through histogram charts) and numerical representation (through tables) of the optimization result for each scenario are provided. In order to provide an overview of the objective function trend, the performance ratios and mass reduction for each scenario were collected into Figure 12 and Figure 13. The mentioned values were obtained for each scenario, making an average of the results, before and after optimization, independently, for each steel profile composing the central pole.
Figure 12 highlights an almost monotone increment of the performance ratios to the number of design variables. Interestingly, for a number of variables n ≥ 5, no significant improvements are achieved. Figure 13 emphasizes an important reduction of structural mass as the design variables increase. Once again, n = 5 represents trade-off. If the number of variables exceed 5, no significant improvements are observed.
Figure 12 and Figure 13 show a comparison between each scenario in terms of the average performance ratio and mass reduction, respectively. Figure 12 highlights the difference with the initial state, which has an average performance ratio P R 0 = 25.6%. An evident improvement is achieved for scenarios that include the thickness t as the design variable.
In particular, from Scenarios D, E, F, G, H, the average performance ratios exceed 50%, resulting in a more than 40% difference compared to the initial state. Figure 12 shows that the commercial profiles are sufficient to accommodate the optimized solution. An exception is noticeable in Scenario A because the optimization is performed using just one diameter Φ , which is optimal for a few parts of the structure, while others are “over-fitted”, resulting in a decrease of the performance ratios −28.4% and an increase of structural mass (+173 kg), as shown in Figure 13.
Similarly, a monotonic increment of the structural mass at the end of the optimization process is evident from Figure 13. In this case, the tonnage decreases with the increasing of the parameter’s number. There is an overall mass reduction of about −67.5% (−972 kg) from scenario D to H. In scenarios A, B, and C, the thickness t of structural members is not considered. Therefore, the mass loss is not satisfactory, at about −28.4% (−409 kg). The choice of the best scenario should depend on one of the five situations described above (from D to H) related to the better PRs gain and mass loss.

7. Conclusions

In this paper, a guyed radio mast’s size and shape optimization process was carried out to identify the equilibrium solution that guarantees the lighter optimized model, verifying strength, instability, and deformation requirements. The paper considers a detailed evaluation of the variable loads according to the Eurocodes recommendations. Furthermore, the OAPI was used to perform a structural analysis with the finite element software SAP2000 by considering the non-linearity of the cables. The optimization was carried out using a genetic optimization algorithm. Eight scenarios (labeled from A to H) were investigated, considering different arrangements of the geometric characteristics of the central pole and cables. The input parameters were increased from Scenario A to H to achieve the best fitness value of the self-weight. From Scenario A to H, the mass reduction index generally increased with the computational effort except in scenarios B and E, in which the input parameter did not represent the best vector design for the structural optimization. At this stage, the best design solution was evaluated from the database of cross-sections inside the finite element software. Though Scenario A provides the worst structural solution in terms of objective function, it represents the most convenient optimization strategy due to its low computational effort; on the contrary, Scenario H exhibits the best fitness value with the lowest self-weight, but it represents the most time-consuming solution. The best solution is achieved when the thickness values of each member, which, composed of the central pole, are included in the optimization process. An improvement of the structural behaviour against instability is observed with increasing thickness. This verification is critical for this structure, mainly subjected to normal stresses resulting from self-weight and pre-stressing cable force. The entire optimization process seems to not be sensible to the pole diameter, chosen as the input parameter of the design vector. Although the final results of the FEM analyses are based on the Italian standards, other codes (e.g., Eurocodes, American code, etc.) can be selected from the SAP2000 settings. However, since no detailed analysis was carried out and many standards are based on the semi-probabilistic approach, the final results should be similar, even with different code formulations. Nevertheless, the partial safety factors involved in load combinations remain quite the same from the numerical point of view, regardless of the followed code.
In future developments, the authors will attempt to replace circular hollow sections with built-up steel solutions to achieve the best structural performance and assemblage procedures. Especially for higher structures, guyed radio masts generally consist of a truss skeleton. Another possible development could be a structural optimization for a cable-stayed radio antenna adopting other optimization strategies, such as particle swarm optimization, PSO, and the evolution differential algorithm (EDA), which could be less time-consuming. Finally, it could perform a typological optimization by managing the position of the cable connection, trying to find the best attachment points.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/app12104875/s1.

Author Contributions

Conceptualization, R.C., M.M.R. and J.M.; methodology, R.C., M.M.R., A.A., J.M., M.L.G. and G.C.M.; software, M.L.G., M.M.R. and R.C.; validation, R.C., M.L.G., J.M. and A.A.; formal analysis, R.C., J.M. and A.A.; investigation, A.A. and R.C.; resources, R.C., M.M.R. and G.C.M.; data curation, R.C. and A.A.; writing—original draft preparation, A.A., R.C. and M.M.R.; writing—review and editing, J.M., R.C., M.L.G. and A.A.; visualization, M.M.R., R.C. and A.A.; supervision, G.C.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by project MSCA-RISE-2020 Marie Skłodowska-Curie Research and Innovation Staff Exchange (RISE)—ADDOPTML (ntua.gr).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable comments and suggestions in revising the paper. The authors would like to thank G.C. Marano and the project ADDOPTML for funding/supporting this research.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Drag and lift forces according to [36] at ULS in [Kg/m].
Table A1. Drag and lift forces according to [36] at ULS in [Kg/m].
Wind Action (Drag D, Lift L) at ULS
z (m)Drag_1Lift_1Drag_2Lift_2Drag_2Lift_3
13.926.56.015.26.851.95
23.926.56.015.26.851.95
33.926.56.015.26.851.95
43.926.56.015.26.851.95
53.926.56.015.26.851.95
64.26.976.445.587.342.09
74.447.376.815.97.762.21
84.667.737.156.198.142.32
94.858.057.446.448.482.42
105.038.357.716.688.792.5
115.198.617.966.899.072.58
125.348.868.197.099.332.66
135.489.098.47.279.572.73
145.619.318.67.459.82.79
155.739.518.797.6110.012.85
165.859.78.977.7610.222.91
175.969.889.137.9110.412.97
186.0610.069.298.0410.593.02
196.1610.229.458.1810.763.07
206.2510.389.598.310.923.11
216.3410.539.738.4211.083.16
226.4310.679.868.5411.233.2
236.5110.819.998.6511.383.24
246.5910.9410.118.7511.523.28
256.6711.0710.238.8611.663.32
266.7511.1910.358.9611.793.36
276.8211.3110.469.0511.913.39
286.8911.4310.569.1412.033.43
296.9611.5410.679.2312.153.46
307.0211.6510.779.3212.273.5
Table A2. Drag and lift forces according to [36] at SLS in [Kg/m].
Table A2. Drag and lift forces according to [36] at SLS in [Kg/m].
Wind Action (Drag D, Lift L) at SLS
z (m)Drag_1Lift_1Drag_2Lift_2Drag_2Lift_3
12.293.813.523.054.011.14
22.293.813.523.054.011.14
32.293.813.523.054.011.14
42.293.813.523.054.011.14
52.293.813.523.054.011.14
62.313.843.543.074.041.15
72.323.863.573.094.061.16
82.343.883.583.14.081.16
92.343.893.63.114.11.17
102.353.93.613.124.111.17
112.363.923.623.134.121.17
122.373.933.633.144.131.18
132.373.943.643.154.141.18
142.383.943.643.154.151.18
152.383.953.653.164.161.19
162.393.963.663.174.171.19
172.393.963.663.174.171.19
182.393.973.673.184.181.19
192.43.983.673.184.191.19
202.43.983.683.194.191.19
212.43.993.683.194.21.2
222.43.993.693.194.21.2
232.4143.693.24.211.2
242.4143.73.24.211.2
252.4143.73.24.211.2
262.414.013.73.214.221.2
272.424.013.713.214.221.2
282.424.013.713.214.231.2
292.424.023.713.214.231.21
302.424.023.723.224.231.21
Table A3. Modal participating mass ratios.
Table A3. Modal participating mass ratios.
Modal Participating Mass Ratios
n. ModesPeriod (s)Frequence (Hz)Part. Mass X (%)Part. Mass Y [%]
13.990.2510.00.28
23.990.2510.830.0
33.990.2510.00.55
43.4730.2880.00.25
53.4730.2880.750.0
63.4720.2880.00.5
72.9290.3410.02.52
82.9250.3420.00.05
92.9160.3432.550.0
100.4372.2909.5826.24
110.4342.30426.419.16
120.2064.8537.164.44
130.2034.9345.644.78
140.1556.43727.380.78
150.1446.9353.6726.10
160.1168.5840.140.00
170.1168.5950.000.32
180.1069.4362.359.82
190.05717.4100.030.08
200.05717.4420.010.01
210.05418.5730.480.45
220.05020.1630.030.00
230.04721.4960.000.00
240.04621.5160.000.00
250.03627.9060.270.11
260.03528.3130.040.13
270.03231.2240.010.48
280.03231.7220.200.00
290.03132.5900.010.06
300.02440.8310.000.00
310.02440.8360.000.00
320.02342.60012.320.00
330.02244.5180.0110.59
Table A4. Load combination.
Table A4. Load combination.
Load Combination
ULS Max 1 1.3 · G 1 + 1.5 · G 2 + 1.5 · W i n d 1 + 1.5 · 0.5 · I c e 1 + 1.5 · 0 · Q M
ULS Max1 2 1.3 · G 1 + 1.5 · G 2 + 1.5 · Q M + 1.5 · 0.6 · W i n d 1 + 1.5 · 0.2 · I c e 1
ULS Min 1 1 · G 1 + 0.8 · G 2 + 1.5 · W i n d 1 + 1.5 · 0.5 · I c e 1 + 1.5 · 0 · Q M
Quake 1 E + G 1 + G 2 + 0 · W i n d 1 + 0 · I c e 1 + 0 · Q M
ULS Max 2 1.3 · G 1 + 1.5 · G 2 + 1.5 · W i n d 2 + 1.5 · 0.5 · I c e 2 + 1.5 · 0 · Q M
ULS Max2 1 1.3 · G 1 + 1.5 · G 2 + 1.5 · Q M + 1.5 · 0.6 · W i n d 2 + 1.5 · 0.2 · I c e 2
ULS Min 1 1 · G 1 + 0.8 · G 2 + 1.5 · W i n d 2 + 1.5 · 0.5 · I c e 2 + 1.5 · 0 · Q M
Quake 2 E + G 1 + G 2 + 0.8 · W i n d 2 + 1.5 · I c e 2 + 1.5 · Q M
Figure A1. Scenarios A, B. In blue, orange, and green, the average PRs, respectively, at the initial condition, after optimization, and the design solution.
Figure A1. Scenarios A, B. In blue, orange, and green, the average PRs, respectively, at the initial condition, after optimization, and the design solution.
Applsci 12 04875 g0a1aApplsci 12 04875 g0a1b
Figure A2. Scenarios C, D. In blue, orange, and green, the average PRs, respectively, at the initial condition, after optimization, and the design solution.
Figure A2. Scenarios C, D. In blue, orange, and green, the average PRs, respectively, at the initial condition, after optimization, and the design solution.
Applsci 12 04875 g0a2
Figure A3. Scenarios E, F, H. In blue, orange and green, the average PRs, respectively, at the initial condition, after optimization, and the design solution.
Figure A3. Scenarios E, F, H. In blue, orange and green, the average PRs, respectively, at the initial condition, after optimization, and the design solution.
Applsci 12 04875 g0a3
Table A5. Scenario A results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one, according to the product list.
Table A5. Scenario A results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one, according to the product list.
SCENARIO A—Optimized Solution
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)12112.56000201
Pole 2 (6–12 m)12112.56000201
Pole 3 (12–18 m)12112.56000201
Pole 4 (18–24 m)12112.56000201
Pole 5 (24–30 m)12112.56000201
Total Mass [kg] Σ 1003
Mass variation [kg]−437Mass variation [%]−30.36
SCENARIO A—Design proposed according to the product list
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)139.712.56000235
Pole 2 (6–12 m)139.712.56000235
Pole 3 (12–18 m)139.712.56000235
Pole 4 (18–24 m)139.712.56000235
Pole 5 (24–30 m)139.712.56000235
Total Mass [kg] Σ 1176
Mass variation [kg]−264Mass variation [%]−18.36
Ntrial = 5
Φ opt [mm]OF [kN]
12140.758
12140.758
12140.758
12140.758
12240.849
Table A6. Scenario B results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
Table A6. Scenario B results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
SCENARIO B—Optimized Solution
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)14912.56000252
Pole 2 (6–12 m)13812.56000231
Pole 3 (12–18 m)12612.56000210
Pole 4 (18–24 m)11512.56000189
Pole 5 (24–30 m)10312.56000168
Total Mass [kg] Σ 1051
Mass variation [kg]−389Mass variation [%]−27.02
SCENARIO B—Design proposed according to product list
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)168.312.56000288
Pole 2 (6–12 m)139.712.56000235
Pole 3 (12–18 m)139.712.56000235
Pole 4 (18–24 m)114.312.56000188
Pole 5 (24–30 m)101.612.56000165
Total Mass [kg] Σ 1111
Mass variation [kg]−329Mass variation [%]−22.84
Ntrial = 5; best solutions
Φ i [mm] Φ f [mm]OF [kN]
1489441.248
14610341.466
1489441.248
14610341.466
1499241.230
Table A7. Scenario C results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
Table A7. Scenario C results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
SCENARIO C—Optimized solution
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)14712.56000249
Pole 2 (6–12 m)13612.56000228
Pole 3 (12–18 m)12512.56000208
Pole 4 (18–24 m)11412.56000188
Pole 5 (24–30 m)10312.56000167
Total Mass [kg] Σ 1040
Mass variation [kg]−400Mass variation [%]−27.79
SCENARIO C—Design proposed according to the product list
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)168.312.56000288
Pole 2 (6–12 m)139.712.56000235
Pole 3 (12–18 m)139.712.56000235
Pole 4 (18–24 m)114.3106000154
Pole 5 (24–30 m)101.6106000135
Total Mass [kg] Σ 1048
Mass variation [kg]−392Mass variation [%]−27.22
Ntrial = 5
Φ i [mm] Φ f [mm]F [kN]OF [kN]
152921.841.393
151921.441.339
14992141.230
156922.441.610
147920.841.121
Table A8. Scenario D results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
Table A8. Scenario D results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
SCENARIO D—Optimized solution
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)16166000138
Pole 2 (6–12 m)14766000125
Pole 3 (12–18 m)13366000113
Pole 4 (18–24 m)12066000101
Pole 5 (24–30 m)1066600089
Total Mass [kg] Σ 565
Mass variation [kg]−875Mass variation [%]−60.75
SCENARIO D—Design proposed according to the product list
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)168.366000144
Pole 2 (6–12 m)168.366000144
Pole 3 (12–18 m)139.766000119
Pole 4 (18–24 m)114.36600096
Pole 5 (24–30 m)101.66600085
Total Mass [kg] Σ 588
Mass variation [kg]−853Mass variation [%]−59.20
Ntrial = 5
Φ i [mm] Φ f [mm]t [mm]OF [kN]
16192636.465
146117737.389
16292636.491
16292636.491
16392636.517
Table A9. Scenario E results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
Table A9. Scenario E results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
SCENARIO E—Optimized Solution
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)16566000141
Pole 2 (6–12 m)15066000128
Pole 3 (12–18 m)13566000115
Pole 4 (18–24 m)12166000102
Pole 5 (24–30 m)1066600089
Total Mass [kg] Σ 574
Mass variation [kg]−866Mass variation [%]−60.13
SCENARIO E—Design proposed according to the product list
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)168.366000144
Pole 2 (6–12 m)168.366000144
Pole 3 (12–18 m)139.766000119
Pole 4 (18–24 m)114.36600096
Pole 5 (24–30 m)101.66600085
Total Mass [kg] Σ 588
Mass variation [kg]−853Mass variation [%]−59.20
Ntrial = 5; best solutions
Φ i Φ f tFOF
[mm][mm][mm][kN][kN]
150977.81.337.766
1531126.41.636.964
1659162.336.552
1609171.337.287
1391048.81.338.337
Table A10. Scenario F results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
Table A10. Scenario F results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
SCENARIO F—Optimized Solution
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)1574600091
Pole 2 (6–12 m)14466000122
Pole 3 (12–18 m)1314600075
Pole 4 (18–24 m)1184600067
Pole 5 (24–30 m)1054600060
Total Mass [kg] Σ 415
Mass variation [kg]−1025Mass variation [%]−71.16
SCENARIO F—Design proposed according to the product list
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)168.34600097
Pole 2 (6–12 m)168.356000121
Pole 3 (12–18 m)139.74600080
Pole 4 (18–24 m)114.34600065
Pole 5 (24–30 m)101.64600058
Total Mass [kg] Σ 421
Mass variation [kg]−1019Mass variation [%]−70.75
Ntrial = 3; best solutions
Φ i Φ f t e n d s t i n t e r FOF
[mm][mm][mm][mm][kN][kN]
15592473.235.141
15792460.934.993
15192471.335.058
Table A11. Scenario H results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
Table A11. Scenario H results: optimized solutions for the different independent executions (Ntrial) and the proposed industrial one according to the product list.
SCENARIO H—Optimized Solution
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)1583600069
Pole 2 (6–12 m)14666000124
Pole 3 (12–18 m)1334600076
Pole 4 (18–24 m)1214600069
Pole 5 (24–30 m)1083600047
Total Mass [kg] Σ 385
Mass variation [kg]−1055Mass variation [%]−73.27
SCENARIO H—Design proposed according to the product list
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)168.356000121
Pole 2 (6–12 m)168.34600097
Pole 3 (12–18 m)139.74600080
Pole 4 (18–24 m)139.73600061
Pole 5 (24–30 m)114.33600049
Total Mass [kg] Σ 408
Mass variation [kg]−1032Mass variation [%]−71.65
Ntrial = 3
Φ i Φ i t 1 t 2 t 3 t 4 t 5 FOF
[mm][mm][mm][mm][mm][mm][mm][kN][kN]
164109454330.934.789
16711136433234.839
1589636443234.695

References

  1. Law, S.; Bu, J.; Zhu, X.; Chan, S.L. Wind characteristics of Typhoon Dujuan as measured at a 50 m guyed mast. Wind Struct. 2006, 9, 387–396. [Google Scholar] [CrossRef]
  2. Sparling, B.; Wegner, L. Estimating peak wind load effects in guyed masts. Wind Struct. 2007, 10, 347–366. [Google Scholar] [CrossRef]
  3. Hensley, G.M.; Plaut, R.H. Three-dimensional analysis of the seismic response of guyed masts. Eng. Struct. 2007, 29, 2254–2261. [Google Scholar] [CrossRef]
  4. Amiri, G.G. Seismic sensitivity indicators for tall guyed telecommunication towers. Comput. Struct. 2002, 80, 349–364. [Google Scholar] [CrossRef]
  5. Liu, C.; Fang, D.; Zhao, L. Reflection on earthquake damage of buildings in 2015 Nepal earthquake and seismic measures for post-earthquake reconstruction. Structures 2021, 30, 647–658. [Google Scholar] [CrossRef]
  6. Sun, X.; Tao, X.; Duan, S.; Liu, C. Kappa (k) derived from accelerograms recorded in the 2008 Wenchuan mainshock, Sichuan, China. J. Asian Earth Sci. 2013, 73, 306–316. [Google Scholar] [CrossRef]
  7. Buchholdt, H.; Moossavinejad, S.; Iannuzzi, A. Non-Linear Dynamic Analysis of Guyed Masts Subjected to Wind and Guy R Uptures. Proc. Inst. Civ. Eng. 1986, 81, 353–395. [Google Scholar] [CrossRef]
  8. Ballaben, J.S.; Guzmán, A.M.; Rosales, M.B. Nonlinear dynamics of guyed masts under wind load: Sensitivity to structural parameters and load models. J. Wind Eng. Ind. Aerodyn. 2017, 169, 128–138. [Google Scholar] [CrossRef]
  9. Gerstoft, P.; Davenport, A. A simplified method for dynamic analysis of a guyed mast. J. Wind Eng. Ind. Aerodyn. 1986, 23, 487–499. [Google Scholar] [CrossRef]
  10. de Souza, R.R.; Miguel, L.F.F.; Lopez, R.H.; Miguel, L.F.F.; Torii, A.J. A procedure for the size, shape and topology optimization of transmission line tower structures. Eng. Struct. 2016, 111, 162–184. [Google Scholar] [CrossRef]
  11. Saudi, G. Structural assessment of a guyed mast through measurement of natural frequencies. Eng. Struct. 2014, 59, 104–112. [Google Scholar] [CrossRef]
  12. Saxena, R.; Popplewell, N.; Trainor, P.; Shah, A. Vibrations of complex guyed towers. In Proceedings of the 12th Biennial Conference on Mechanical Vibration and Noise Control, Montreal, QC, Canada, 17–21 September 1989. [Google Scholar]
  13. Novak, M.; Tanaka, H.; Davenport, A.G. Vibration of towers due to galloping of iced cables. J. Eng. Mech. Div. 1978, 104, 457–473. [Google Scholar] [CrossRef]
  14. Wahba, Y.M.; Madugula, M.K.; Monforton, G.R. Shake Table for Dynamic Testing of Guyed Towers. In Building to Last; ASCE: Reston, VA, USA, 1997; pp. 353–357. [Google Scholar]
  15. Madugula, M.K.; Wahba, Y.M.; Monforton, G.R. Dynamic response of guyed masts. Eng. Struct. 1998, 20, 1097–1101. [Google Scholar] [CrossRef]
  16. Luzardo, A.C.; Parnás, V.E.; Rodríguez, P.M. Guy tension influence on the structural behavior of a guyed mast. J. Int. Assoc. Shell Spat. Struct. 2012, 53, 111–116. [Google Scholar]
  17. Davenport, A.; Sparling, B. Dynamic gust response factors for guyed towers. J. Wind Eng. Ind. Aerodyn. 1992, 43, 2237–2248. [Google Scholar] [CrossRef]
  18. Harikrishna, P.; Annadurai, A.; Gomathinayagam, S.; Lakshmanan, N. Full scale measurements of the structural response of a 50 m guyed mast under wind loading. Eng. Struct. 2003, 25, 859–867. [Google Scholar] [CrossRef]
  19. Gioffrè, M.; Gusella, V.; Materazzi, A.; Venanzi, I. Removable guyed mast for mobile phone networks: Wind load modeling and structural response. J. Wind Eng. Ind. Aerodyn. 2004, 92, 463–475. [Google Scholar] [CrossRef]
  20. Clobes, M.; Peil, U. Unsteady buffeting wind loads in the time domain and their effect on the life-cycle prediction of guyed masts. Struct. Infrastruct. Eng. 2011, 7, 187–196. [Google Scholar] [CrossRef]
  21. Pezo, M.L.; Bakić, V.V. Numerical determination of drag coefficient for guyed mast exposed to wind action. Eng. Struct. 2014, 62, 98–104. [Google Scholar] [CrossRef]
  22. Sparling, B.F. The Dynamic Behavior of Guys and Guyed Masts in Turbulent Winds. Ph.D. Thesis, Western University, London, ON, Canada, 1995. [Google Scholar]
  23. Wahba, Y.; Madugula, M.; Monforton, G. Evaluation of non-linear analysis of guyed antenna towers. Comput. Struct. 1998, 68, 207–212. [Google Scholar] [CrossRef]
  24. Madugula, M.K. Dynamic Response of Lattice Towers and Guyed Masts; ASCE Publications: Reston, VA, USA, 2001. [Google Scholar]
  25. Orlando, D.; Gonçalves, P.B.; Rega, G.; Lenci, S. Nonlinear dynamics and instability as important design concerns for a guyed mast. In IUTAM Symposium on Nonlinear Dynamics for Advanced Technologies and Engineering Design; Springer: Berlin/Heidelberg, Germany, 2013; pp. 223–234. [Google Scholar]
  26. Ballaben, J.S.; Rosales, M.B. Nonlinear dynamic analysis of a 3D guyed mast. Nonlinear Dyn. 2018, 93, 1395–1405. [Google Scholar] [CrossRef]
  27. Belevičius, R.; Jatulis, D.; Šešok, D. Optimization of tall guyed masts using genetic algorithms. Eng. Struct. 2013, 56, 239–245. [Google Scholar] [CrossRef]
  28. Gawronski, W.; Bienkiewicz, B.; Hill, R. Wind-induced dynamics of a deep space network antenna. J. Sound Vib. 1994, 178, 67–77. [Google Scholar] [CrossRef]
  29. Fujino, Y.; Warnitchai, P.; Pacheco, B. Active Stiffness Control of Cable Vibration. J. Appl. Mech. 1993, 60, 948–953. [Google Scholar] [CrossRef]
  30. Lacarbonara, W.; Ballerini, S. Vibration mitigation of guyed masts via tuned pendulum dampers. Struct. Eng. Mech. Int. J. 2009, 32, 517–529. [Google Scholar] [CrossRef]
  31. Błachowski, B. Model based predictive control of guyed mast vibration. J. Theor. Appl. Mech. 2007, 45, 405–423. [Google Scholar]
  32. Bell, L.C.; Brown, D.M. Guyed tower optimization. Comput. Struct. 1976, 6, 447–450. [Google Scholar] [CrossRef]
  33. Thornton, C.H.; Joseph, L.; Scarangello, T. Optimization of tall structures for wind loading. J. Wind Eng. Ind. Aerodyn. 1990, 36, 235–244. [Google Scholar] [CrossRef]
  34. Uys, P.; Farkas, J.; Jarmai, K.; Van Tonder, F. Optimisation of a steel tower for a wind turbine structure. Eng. Struct. 2007, 29, 1337–1342. [Google Scholar] [CrossRef]
  35. Venanzi, I.; Materazzi, A. Multi-objective optimization of wind-excited structures. Eng. Struct. 2007, 29, 983–990. [Google Scholar] [CrossRef]
  36. Zhang, Z.Q.; Li, H.N. Two-level optimization method of transmission tower structure based on ant colony algorithm. In Advanced Materials Research; Trans Tech Publications Ltd.: Freienbach, Switzerland, 2011; Volume 243, pp. 5849–5853. [Google Scholar]
  37. Cucuzza, R.; Costi, C.; Rosso, M.M.; Domaneschi, M.; Marano, G.C.; Masera, D. Optimal strengthening by steel truss arches in prestressed girder bridges. In Proceedings of the Institution of Civil Engineers—Bridge Engineering; Thomas Telford Ltd.: London, UK, 2021; pp. 1–21. [Google Scholar] [CrossRef]
  38. Luh, G.C.; Lin, C.Y. Optimal design of truss-structures using particle swarm optimization. Comput. Struct. 2011, 89, 2221–2232. [Google Scholar] [CrossRef]
  39. Manuello Bertetto, A.; Marano, G. Numerical and dimensionless analytical solutions for circular arch optimization. Eng. Struct. 2022, 253, 113360. [Google Scholar] [CrossRef]
  40. Kaveh, A.; Talatahari, S. Particle swarm optimizer, ant colony strategy and harmony search scheme hybridized for optimization of truss structures. Comput. Struct. 2009, 87, 267–283. [Google Scholar] [CrossRef]
  41. Deng, Z.Q.; Zhang, Y.; Huang, H.L.; Li, B. Parametric optimization for a tapered deployable mast in an integrated design environment. In Advanced Materials Research; Trans Tech Publications Ltd.: Freienbach, Switzerland, 2012; Volume 346, pp. 426–432. [Google Scholar]
  42. Guo, H.; Li, Z. Structural topology optimization of high-voltage transmission tower with discrete variables. Struct. Multidiscip. Optim. 2011, 43, 851–861. [Google Scholar] [CrossRef]
  43. delle Ricerche, C.N. Istruzioni per la valutazione delle azioni e degli effetti del vento sulle costruzioni. CNR-DT 2009, 207, 2008. [Google Scholar]
  44. Mordà, N.; Mancini, A. Norme Tecniche per le Costruzioni (NTC 2018) D. Min. Infrastrutture e Trasporti 17 Gennaio 2018; Ministero delle Infrastrutture e dei Trasporti: Roma, Italy, 2018. [Google Scholar]
  45. Melchiorre, J.; Bertetto, A.M.; Marano, G.C. Application of a Machine Learning Algorithm for the Structural Optimization of Circular Arches with Different Cross-Sections. J. Appl. Math. Phys. 2021, 9, 1159–1170. [Google Scholar] [CrossRef]
  46. Rosso, M.M.; Cucuzza, R.; Aloisio, A.; Marano, G.C. Enhanced Multi-Strategy Particle Swarm Optimization for Constrained Problems with an Evolutionary-Strategies-Based Unfeasible Local Search Operator. Appl. Sci. 2022, 12, 2285. [Google Scholar] [CrossRef]
  47. Rosso, M.M.; Cucuzza, R.; Di Trapani, F.; Marano, G.C. Nonpenalty machine learning constraint handling using PSO-svm for structural optimization. Adv. Civ. Eng. 2021, 2021, 6617750. [Google Scholar] [CrossRef]
  48. Rao, S.S. Engineering Optimization: Theory and Practice; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar]
  49. Aloisio, A.; Pasca, D.P.; Battista, L.; Rosso, M.M.; Cucuzza, R.; Marano, G.; Alaggio, R. Indirect assessment of concrete resistance from FE model updating and Young’s modulus estimation of a multi-span PSC viaduct: Experimental tests and validation. Elsevier Struct. 2022, 37, 686–697. [Google Scholar] [CrossRef]
  50. Sardone, L.; Rosso, M.M.; Cucuzza, R.; Greco, R.; Marano, G.C. Computational Design of Comparative models and geometrically constrained optimization of a multi-domain variable section beam based on Timoshenko model. In Proceedings of the EUROGEN2021, 14TH ECCOMAS Thematic Conference on Evolutionary and Deterministic Methods for Design, Optimization and Control, Athens, Greece, 28–30 June 2021. [Google Scholar] [CrossRef]
  51. Christensen, P.W.; Klarbring, A. An Introduction to Structural Optimization; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2008; Volume 153. [Google Scholar]
  52. Coello Coello, C.A. Theoretical and numerical constraint-handling techniques used with evolutionary algorithms: A survey of the state of the art. Comput. Methods Appl. Mech. Eng. 2002, 191, 1245–1287. [Google Scholar] [CrossRef]
  53. Koziel, S.; Michalewicz, Z. Evolutionary Algorithms, Homomorphous Mappings, and Constrained Parameter Optimization. Evol. Comput. 1999, 7, 19–44. [Google Scholar] [CrossRef] [PubMed]
  54. Michalewicz, Z.; Fogel, D. How to Solve It: Modern Heuristics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  55. Parsopoulos, K.; Vrahatis, M. Unified Particle Swarm Optimization for Solving Constrained Engineering Optimization Problems. In International Conference on Natural Computation; Springer: Berlin/Heidelberg, Germany, 2005; Volume 3612, pp. 582–591. [Google Scholar] [CrossRef] [Green Version]
  56. Coello, C. Self-adaptive penalties for GA-based optimization. In Proceedings of the 1999 Congress on Evolutionary Computation-CEC99 (Cat. No. 99TH8406), Washington, DC, USA, 6–9 July 1999; IEEE: Hoboken, NJ, USA, 1999; Volume 1, pp. 573–580. [Google Scholar] [CrossRef]
Figure 1. (a) Render model realized using Tekla Structures. (b) Technical drawing of the structure investigated with dimensions in mm.
Figure 1. (a) Render model realized using Tekla Structures. (b) Technical drawing of the structure investigated with dimensions in mm.
Applsci 12 04875 g001
Figure 2. Pictures and details of the considered structure.
Figure 2. Pictures and details of the considered structure.
Applsci 12 04875 g002
Figure 3. Mechanical stress–strain behaviour of steel S355 implemented in SAP2000.
Figure 3. Mechanical stress–strain behaviour of steel S355 implemented in SAP2000.
Applsci 12 04875 g003
Figure 4. Elastic response spectrum corresponding to the service limit state (SLV), where S a is the spectral acceleration.
Figure 4. Elastic response spectrum corresponding to the service limit state (SLV), where S a is the spectral acceleration.
Applsci 12 04875 g004
Figure 5. (a) Mode 10th—Ts = 0.437 s—mass participation ratio X = 9.6%, Y = 26.2%; (b) Mode 11th—Ts = 0.434 s—mass participation ratio X = 26.4% Y = 9.2%; (c) Mode 12—Ts = 0.206 s—mass participation ratio X = 7.2% Y = 4.4%.
Figure 5. (a) Mode 10th—Ts = 0.437 s—mass participation ratio X = 9.6%, Y = 26.2%; (b) Mode 11th—Ts = 0.434 s—mass participation ratio X = 26.4% Y = 9.2%; (c) Mode 12—Ts = 0.206 s—mass participation ratio X = 7.2% Y = 4.4%.
Applsci 12 04875 g005
Figure 6. (a) Axial force, (b) shear force ( V 2 ), (c) bending moment ( M 2 ), (d) shear force ( V 3 ), (e) bending moment ( M 3 ).
Figure 6. (a) Axial force, (b) shear force ( V 2 ), (c) bending moment ( M 2 ), (d) shear force ( V 3 ), (e) bending moment ( M 3 ).
Applsci 12 04875 g006
Figure 7. Performance ratios of the pole before optimization. Cables are depicted with magenta colour because their performance ratios are not included in the current representation.
Figure 7. Performance ratios of the pole before optimization. Cables are depicted with magenta colour because their performance ratios are not included in the current representation.
Applsci 12 04875 g007
Figure 8. Displacements vs. elevation at the service limit state in two in-plane orthogonal directions ( u 1 , u 2 ) and their combination (tot).
Figure 8. Displacements vs. elevation at the service limit state in two in-plane orthogonal directions ( u 1 , u 2 ) and their combination (tot).
Applsci 12 04875 g008
Figure 9. Workflow of the optimization problem.
Figure 9. Workflow of the optimization problem.
Applsci 12 04875 g009
Figure 10. Parametric study on the design variables involved and representation of the different scenarios described in Table 5.
Figure 10. Parametric study on the design variables involved and representation of the different scenarios described in Table 5.
Applsci 12 04875 g010
Figure 11. Scenario G.—PRs trend. In blue—the performance ratios of each pole before optimization are illustrated, otherwise orange for the optimized solution. In green—PRs at a design configuration according to the product list.
Figure 11. Scenario G.—PRs trend. In blue—the performance ratios of each pole before optimization are illustrated, otherwise orange for the optimized solution. In green—PRs at a design configuration according to the product list.
Applsci 12 04875 g011
Figure 12. In blue, orange, and green, the average PRs, respectively, at the initial condition, after optimization, and design solution.
Figure 12. In blue, orange, and green, the average PRs, respectively, at the initial condition, after optimization, and design solution.
Applsci 12 04875 g012
Figure 13. Increasing the number of design variables, the final mass becomes gradually smaller, until 385 kg (scenario H).
Figure 13. Increasing the number of design variables, the final mass becomes gradually smaller, until 385 kg (scenario H).
Applsci 12 04875 g013
Table 1. Technical specifications of the steel ropes.
Table 1. Technical specifications of the steel ropes.
Steel Ropes (Cables)
Model6 × 36WS + IWRC/265 wires
Construction pattern6 × (14 + (7 + 7) + 7 + 1) + (7 × 7)
Winding directionright cross
Materialgalvanized steel
Resistance1170 N/mm 2 –180 kg/mm 2
Cable diameterWeightAreaWire diameterLoad to failure
[mm][kg/m][mm 2 ][mm][kN]
161.36173.250.91161
181.67212.741.03204
202.02257.321.14252
222.41307.011.26305
Table 2. Computation of the dead loads.
Table 2. Computation of the dead loads.
Computation of Dead Loads
Profile [mm]w [kg/m]Length [m]Wtot [kg]
CircularD168.3 × 12.548651440
D168.3 × 12.5485.652543
Rectangular60 × 40 × 34.353.169124
60 × 40 × 34.351.8971
100 × 40 × 36.130.45617
RopeD161.366712.45351
D161.366715.44363
D161.366724.439300
D161.36675.76324
D161.36678.46335
2651 Kg
Table 3. Weight of equipment, H, W, and D stand for height, width, and depth.
Table 3. Weight of equipment, H, W, and D stand for height, width, and depth.
TypologyModelNoElevation [m]H × W × D [mm]Self-Weight [kg]Clamps [kg]Total [kg]
AntennaAOC4518R7v06329.251509 × 469 × 20639.32 × 5.8153
Antenna6888670N3261997 × 305 × 163322 × 3.9119
Parabolan.d323.15Diameter = 300152.251.6
Table 4. Non-structural dead loads.
Table 4. Non-structural dead loads.
Itemqk [kN/m]Qk [kN]
Steel ladder, other0.3-
Antenna-1.53
Antenna-1.19
Parabolas-0.52
Table 5. Parametric study on the design variables involved and summary of the different scenarios.
Table 5. Parametric study on the design variables involved and summary of the different scenarios.
ScenarioNo. Parameters
A ( Φ )1
B ( Φ i , Φ f )2
C ( Φ i , Φ f , F )3
D ( Φ i , Φ f , t )3
E ( Φ i , Φ f , t , F )4
F ( Φ i , Φ f , t e n d s , t i n t e r , F )5
G ( t 1 , t 2 , t 3 , t 4 , t 5 )5
H ( Φ i , Φ f , t 1 , t 2 , t 3 , t 4 , t 5 )8
Table 6. Total mass of the main pole.
Table 6. Total mass of the main pole.
ParameterMeasureValue
Φ 0 [mm]168.3
t 0 [mm]12.5
L[mm]6000
Mass[kg]288
no elements[-]5
Total Mass[kg]1440
Table 7. Averaged performance ratios obtained in each optimization scenario.
Table 7. Averaged performance ratios obtained in each optimization scenario.
No ParametersPR InitialPR OptimizedPR Design
[%][%][%]
128.045.740.5
239.543.1
350.550.6
454.458
565.860.2
86866
Table 8. Mass values before/after optimization and after proper approximation (design) using commercial steel profiles.
Table 8. Mass values before/after optimization and after proper approximation (design) using commercial steel profiles.
No ParametersInitial Mass [kg]Optimized Mass [kg]Design Mass [kg]
1144010031176
210511111
3803818
4574588
5403453
8385408
Table 9. Scenario G results: optimized solutions for the different independent executions (Ntrial) and proposed industrial one, according to the product list.
Table 9. Scenario G results: optimized solutions for the different independent executions (Ntrial) and proposed industrial one, according to the product list.
SCENARIO G—Optimized Solution
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)168.33600073
Pole 2 (6–12 m)168.34600097
Pole 3 (12–18 m)168.33600073
Pole 4 (18–24 m)168.33600073
Pole 5 (24–30 m)168.33600073
Total Mass [kg] Σ 391
Mass variation [kg]−1050Mass variation [%]−72.88
SCENARIO G—Design proposed according to the product list
Elementd [mm]t [mm]L [mm]Mass [Kg]
Pole 1 (0–6 m)168.34600097
Pole 2 (6–12 m)168.34600097
Pole 3 (12–18 m)168.33600073
Pole 4 (18–24 m)168.33600073
Pole 5 (24–30 m)168.33600073
Total Mass [kg] Σ 414
Mass variation [kg]−1026Mass variation [%]−71.22
Ntrial = 3
t 1 t 2 t 3 t 4 t 5 OF
[mm][mm][mm][mm][mm][kN]
3443334.985
3433334.751
3434334.985
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cucuzza, R.; Rosso, M.M.; Aloisio, A.; Melchiorre, J.; Giudice, M.L.; Marano, G.C. Size and Shape Optimization of a Guyed Mast Structure under Wind, Ice and Seismic Loading. Appl. Sci. 2022, 12, 4875. https://doi.org/10.3390/app12104875

AMA Style

Cucuzza R, Rosso MM, Aloisio A, Melchiorre J, Giudice ML, Marano GC. Size and Shape Optimization of a Guyed Mast Structure under Wind, Ice and Seismic Loading. Applied Sciences. 2022; 12(10):4875. https://doi.org/10.3390/app12104875

Chicago/Turabian Style

Cucuzza, Raffaele, Marco Martino Rosso, Angelo Aloisio, Jonathan Melchiorre, Mario Lo Giudice, and Giuseppe Carlo Marano. 2022. "Size and Shape Optimization of a Guyed Mast Structure under Wind, Ice and Seismic Loading" Applied Sciences 12, no. 10: 4875. https://doi.org/10.3390/app12104875

APA Style

Cucuzza, R., Rosso, M. M., Aloisio, A., Melchiorre, J., Giudice, M. L., & Marano, G. C. (2022). Size and Shape Optimization of a Guyed Mast Structure under Wind, Ice and Seismic Loading. Applied Sciences, 12(10), 4875. https://doi.org/10.3390/app12104875

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