Optimization of Rigidly Supported Guyed Masts

A technique for simultaneous topology, shape, and sizing optimization of tall guyed masts is presented in this paper. The typical scheme of telecommunication masts that are triangular in plan and supported by a certain number of guys’ levels is examined.The objective function is the total mass of the mast, including the mass of the guys. The mast structure is optimized for self-weight and wind loading that is evaluated according to Eurocodes. The nonlinear behaviour of the guyed mast is simplified idealizing the nonlinear guys as approximate boundary conditions for the mast. A comparison of the simplified solution with the results of nonlinear analysis with Ansys shows small discrepancies that are on the safe side. The constraints involve all strength, local and global stability, and slenderness requirements. The optimization problem is solved using evolutionary algorithm with original genome repair procedure. As an example, a typical 96m guyed broadcasting antenna for mobile-phone networks was designed employing the proposed optimization technique and taking up to 10 guys’ levels. It is shown that the optimal mast is supported by 3–5 guys’ levels. The optimal ranges of all remaining design variables were also obtained.


Introduction
In recent years the development of telecommunication networks has stimulated the amount of design and construction of tall steel-guyed masts.The mast structures are usually produced and constructed in certain quantities as typical structures depending on the terrain conditions and required area of antennas.Therefore, research, innovation, and cost optimization of guyed masts are a relevant engineering problem.
The guyed masts are supported by prestressed guy cables attached to the shaft of the mast at different levels.These cables exhibit nonlinear behaviour and are very sensitive to their geometry and prestress levels; therefore, despite the elastic behaviour of the mast shaft, these structures should be analysed by a nonlinear analysis program.The shaft is usually composed of a number of sections of identical structure and dimensions; one section comprises shaft legs, stiffeners between legs of the shaft, and oblique bracing elements.One example of triangular mast and the structure of the shaft are shown in Figures 1 and 2.
The problem of analysis of a guyed mast is even more complicated since the dominant loading on the mast is turbulent wind, and the dynamic effects cannot be neglected in a reliable analysis.
Those nonlinear prestressed structures were analysed in a number of research papers, for example, Gantes et al. [1] and Smith [2].Different analytical and computational methods for evaluation of the strain-stress behaviour were suggested by Wahba et al. [3], He et al. [4], and Gioffrè et al. [5].Several researchers pursued improvements of mast structure seeking the least possible weight (synonymous with cost) of the structure (Jasim and Galeb [6]).In most cases refinement of the structure was understood exclusively as selection of the geometrical scheme of shaft elements, that is, leg, bracing, and stiffener members, and the dimensioning of cross sections.However, the maximum effectiveness of those structures can be achieved by simultaneously tuning all the topological, geometric, and physical parameters.
The topology of the guyed mast is identified by the number of the guys' clusters, the number of guys' anchors to the ground, the geometric scheme of a typical shaft section, 2 and the shaft-ground connection scheme.Some of these parameters were chosen before the optimization.For example, it is obvious (and it is explicitly shown in [7]) that the mast of the triangular horizontal plan is lighter than the square plan mast.Thus, only the triangular mast is considered.As to the shaft-ground connection scheme, two schemes prevail in the engineering practice: pinned at the foundation and fixed.

Advances in Civil Engineering
In the pinned masts the stresses in the leg are distributed more evenly.In fixed masts the horizontal displacements of the shaft cause significant bending moments at the support, and the maximum stresses develop in the lowest sections of the mast.Such uneven distribution of stresses is not rational and results in a substantially heavier structure compared with the pinned masts.However, pinned masts have inherent deficiencies.Firstly, the shaft must be strutted at the assembly stage of mast construction.Secondly, the hinged connection at the support does not provide the torsional stiffness of the mast; therefore, either the antitwist tackle ("mounting star") or sophisticated support construction that does not transfer the bending moment but assures the required torsional stiffness must be provided [8].Due to these drawbacks of the pinned scheme, in this paper the mast scheme that is fixed at the foundation is considered.
The shape of the mast is defined by the vertical locations of the guy attachments to the shaft, the horizontal locations of the guy anchors on the ground, the width of the shaft, and the height of a typical section.The sizing parameters include cross sections of all structural elements of the shaft and cross sections of guys.All leg, bracing, and stiffener elements are round solid steel bars and are of constant diameter along the mast height.In addition, the pretension forces have to be chosen in guys at all levels.
The aim of this paper is twofold: to suggest the technique enabling the use of a classical EA with special genome repair procedure for optimization of guyed masts and to obtain the optimal scheme of a typical 96 m tall mast.The suggested technique should yield the optimal, that is, lightest, scheme for a mast with "one button click" using the following initial data for the structure: In mathematical terms, the mast optimization problem is a mixed integer global optimization problem, since the landscape of objective function is not known and cannot be obtained in a closed form.The results of numerical experiments presented in this paper definitely show that the problem is really multimodal.Since the number of design parameters is large, we use the stochastic optimization algorithms, which do not require the sensitivity information.The stochastic algorithms do not guarantee the global solution of a problem, but in the engineering practice it is more important to find a rational solution that is better than the solution currently known.
In [9,10] we optimized the triangular mast with two clusters of guys using genetic algorithms and simulated annealing and revealed the dominant parameters for optimal design.In [11] the typical 96 m broadcasting antenna for mobile-phone networks with an arbitrary number of guys' clusters employing an evolutionary algorithm (EA) from Matlab [12] was optimized.Comparison of optimized mast structure with the typical industrial design revealed possibilities for significant reduction of steel expenditure.However, the simplified algorithm for evaluation of computational length between adjacent clusters of guys led to a so-called "checkerboard effect": finer structural elements conditioned better objective function values, and some threshold values were needed to prevent excessive refinement of the structure.In this paper the exact computational lengths were used for checking the global stability constraints, and a special genome repair procedure was employed in EA.It is shown that the mass of the mast is significantly reduced by installing at least three guys' clusters: the reduction of mass is up to 40% for a three-guy-level mast compared to a two-level one, but starting from the five-level mast the mass slightly increases.Also, the optimal ranges of the determinant design parameters are discussed.

Problem Formulation
2.1.Idealizations.The structural behaviour of guyed masts is extremely complicated.The guys exhibit geometrically nonlinear behaviour and more at low pretensioning levels.Increasing the pretension forces decreases the nonlinearity and enhances lateral stiffness but simultaneously increases compressive loads and herewith the buckling probability for the mast.The mast itself can be also geometrically nonlinear due its slender structure.In addition to the substantial wind loads, the loads due to the self-weight of the mast with all auxiliary equipment and possible icing should be accounted for.
The use of optimization algorithms with tens of design parameters supposes analysing the mast structure and obtaining the objective function for any number of times.Therefore, a very fast and reliable analysis tool is a pure necessity, and we restrict the analysis only to the linear stage, approximately substituting the guys by the springs and compression forces on the shaft of the mast.The horizontal stiffness of springs is proportional to the displacements at the guy attachment node, and the vertical compression forces depend first of all on the prestress in the guy.The geometrically nonlinear verification solution of the same mast structure via finite element package Ansys proves these approximate assumptions are always on the safe side.In the Ansys solution the guys are represented by cable elements; each guy is divided into 10 elements.In Figure 3 both solutions for the horizontal displacements in the opposite direction to the wind are compared for the obtained best mast scheme: the 5 guys' cluster mast.All geometrical data are presented in Table 2.The discrepancies in horizontal displacements are higher at the middle levels of the mast, but the differences in displacements of nodes at the mast top do not exceed 2% (Figure 3).Since higher displacements obtained by linearized simplified solution also cause higher stresses in the rigidly supported mast, the constraints (see ( 6)-( 8), Section 2.2) slightly narrow the search space of optimization problem to the safer side.More important here is the coincidence of displacements at the mast top: analysis of the intermediate results of optimization reveals that the constraint most difficult to satisfy is the given maximum lateral displacement of the mast top.
The next simplifying assumption concerns the turbulent wind loading on the structure.Instead of dynamic loading, we evaluate the corresponding static wind loads by the patchload method (Eurocode 1, Part 1-4 [13], and Eurocode 3, Part 3-1 [14]) which adequately mimics the dynamic behaviour of the structure and solve the static problem.According to the Eurocodes, the designer has to consider several other loadings in addition to the wind loading of different directions, including icing and several of their combinations.In this paper we choose to optimize the mast for the most critical and determinant case of wind loading, when the direction of the wind is at a right angle to one side of the mast.The loads due to the self-weight of the structure and equipment are accounted for, but not the ice loads.Simultaneous icing and extreme winds are more topical for Nordic and mountainous countries only.However, our meshing program has the straightforward capability of including icing loads of the required type and class according to the European standard on atmospheric icing [15].
Subsequently, after rounding the obtained values of continuous design parameters to suitable values (e.g., geometrical dimensions, to in-stock values), the designer may verify the scheme employing the more accurate nonlinear analysis and taking into account other loading cases.Therefore, the optimization technique presented can be useful during the preliminary mast design stage, supplying the designer with hints for the topology and shape of the mast, as well as geometrical data on all structural elements.

Optimization Problem and Technique.
The complete set of programs consists of four independent parts: the optimization algorithm, meshing program that renders the set of design parameters to the discretized model of the mast The genetic algorithm with preprocessor for genome repair is used as the primary optimization engine.After the best solution from the set of EA solutions is selected, the pattern search algorithm [16] is employed to explore the neighbourhood of the solution in depth and (probably) improve it.
The meshing program is used as a preprocessor for the finite element program.It inputs the set of design parameters produced by the optimization algorithm together with data which do not vary in the optimization process (e.g., mast height, wind velocity, and material characteristics) and prepares the complete dataset for the finite element program.All loadings including self-weight and wind loads on the structure are recalculated using new values of all geometrical parameters.
A fast finite element program is required; hence an original problem-oriented Fortran program was developed for the static analysis.All the structural elements of the shaft are represented by beam elements.Instead of guys at the guy attachment nodes, the horizontal stiffness boundary conditions and corresponding vertical compression forces are applied.
After the finite element analysis, the postprocessor program evaluates the value of the objective function, verifies all constraints, and transmits the objective function value to the optimization engine, which prepares a new estimate of design parameters set on its basis.
In mathematical terms, the optimization problem is formulated as follows: subject to (i) structure equilibrium constraints, (ii) strength constraints in all structural elements of the shaft and guys, (iii) local stability constraints in all leg and bracing elements, (iv) global stability constraints between the clusters of guys and between the lower cluster and mast foundation, (v) slenderness constraints in all leg and bracing elements, (vi) lateral stability constraint of the mast top.
() in ( 1) is a nonlinear objective function of continuous and integer variables :R  → R,  is the number of design parameters , and  ∈ R  is a feasible region of design parameters.The global minimum should be found.No assumptions on unimodality are included into formulation of the problem; that is, many local minima may exist.
The total mass of the mast, including the mass of the guys, is considered in the objective function.Since the material of the guys is more expensive, the mass of the guys is premultiplied by a factor given by the program user.The cost of guy anchors and the cost of connections can also be included into the objective function in a similar manner.
The complete set of design parameters, along with their characteristics and bounds, is listed in Section 3, Table 1.
All strength, stability, and slenderness constraints are assessed according to Eurocode 3, Part 3-1 [14], and Eurocode 3, Part 1-1 [17].The overall structure equilibrium constraints are checked solving the static problem with the finite element program.The typical 2-node beam elements with 12 degrees of freedom are used to represent the structural elements of the shaft.The main statics equation is where []  is the stiffness matrix, {}  are nodal displacements, and {}  are the active forces of the ensemble of elements.The element stiffness matrix can be found in many textbooks, for example, [18].The stress resultants in all finite elements according to [14] are where   is stress due to the static wind load on the whole mast height plus the dead weight of the mast structure and   is cumulative stress due to the pseudo-dynamic loads on the patches of the mast according to [13]; the number of load cases for   guys' clusters is 2(  + 1) + 1.Thus, (2) has to be solved for 2(  + 1) + 2 times.The influence of guys on the behaviour of the shaft is modelled by linear stiffness boundary conditions.The total stiffness estimation of the guys in a guys' cluster when the wind blows toward one guy is given by Gantes et al. [1].However, this case is not critical for the mast fixed at the foundation.The shaft of the mast receives larger displacements in the opposite wind direction.Now two stressed guys absorb all wind action, while the third guy is slack and provides only compression force mainly due its pretension force.In this case we obtain different stiffness' expression of one windward guy: where   is the distance of the guy foundation from the axis of the shaft,   is guy cross section area,  is Young's modulus of guy material,  is the length of the guy,  is the dead weight of the guy per unit length, and   is the pretension force in the guy.
The compressive effects of one set of guys are idealized by the vertical compressive forces applied to the guy attachment nodes: where   is the number of guys in one guys' cluster and H is the height of the mast.Our assumptions (4) and ( 5) are checked via nonlinear modelling of the mast with guys and are of sufficient precision (Figure 3).The spring stiffness according to (4) is slightly less than the stiffness of the real guy but the error is to the safe side.The first of the inequality constraints, the strength constraints, are checked in all structural elements of the mast, accounting for both the axial force N and bending moment  for the allowable stress : The stability requirements for the leg and bracing members are posed according to Eurocode 3, Part 3-1 [14]: where  is the cross section area of the member and   is the steel yield point (dependent on the diameter of the member).According to [14], the partial factor of resistance of the members to member buckling  1 is taken to be 1.0, and reduction factor coefficient  is evaluated via complex engineering formulas depending on the nondimensional slenderness of structural elements.The global stability requirements between adjacent clusters of guys are expressed in the form of a mast segment buckling condition due to the equivalent compression force   and the equivalent bending moment   in the cross section of the mast: where   is the reduction factor for the relevant buckling mode,   is the cross section axial resistance, and  , is the cross section moment resistance (based on either the plastic, elastic, or effective section modulus, depending on classification).Details on the evaluation of these coefficients may be found in Eurocode 3, Part 1-1 [14].
The slenderness constraints for the mast segments between clusters of guys and slenderness constraints for bracing members, correspondingly, are Obtaining the buckling length of a mast segment needed for assessment of the stability constraint takes significant numerical effort.In [11] we used an approximate value of buckling length that was equal to the segment length; this is adequately accurate for a mast supported by 1-4 guys' clusters.For a mast supported by a larger number of guys' clusters it led to unrealistic designs, when the mass of the mast steadily decreases with the increasing number of guys' clusters.In this research the ith segment buckling length is obtained with an empirical computational length coefficient: where   is the coefficient of the segment support stiffness with regard to the shaft stiffness , dependent on the total stiffness provided by the guys ]  : where   is the length of the ith segment between adjacent guys' clusters.The coefficient  1 for the segment between foundation and the first guys' level is obtained with a more awkward procedure: where with   being the number of segments.Finally, the lateral displacement  of the mast top is constrained to one percent of the height of the mast.
The classical real coded genetic algorithm for solving mixed integer optimization problems MI-LXPM [19] was employed as a primary optimization engine.The meshing, FEM analysis, and postprocessor programs are connected to the optimization algorithm as "black boxes."The main genetic parameters, population size, number of generations, probabilities of mutation and crossover, and so on, are tuned to the problem solved.The complete list of genetic parameters is presented in Section 3.
The list of design parameters along with their lower and upper bounds is presented in Table 1.Some of the parameters are interdependent (e.g., the guys cannot intersect or be attached to a nonexisting shaft section); therefore, in the initial stage of optimization-random generation of individuals-only the viable or repaired individuals are included into the initial population.For example, at the random generation of one particular individual of the initial population the following data was obtained: the number of guys' clusters   = 3, number of sections   = 50 (let the upper bound for   be 100), and levels of guys' clusters  1 = 10,  2 = 60, and  3 = 80.In order to repair this infeasible individual, the levels are scaled down proportionally:  1 = 5,  2 = 30, and  3 = 40.
Different conditions hold for the invalid individuals of the next generations obtained due to the mutation or crossover.Since the main idea of optimization technique to keep all constituent programs independent, the obtained infeasible guys' levels are repaired but the individuals with intersecting guys are retained.However, with an additional sufficient penalty given by the preprocessor of the finite element program, the finite element analysis for these individuals is not executed.The penalty is proportional to the number of such violations and acts like a "dead" penalty, fairly reducing the chances of infeasible individuals surviving the selection.
After the best EA solution is selected, the neighbourhood of the solution is explored with the pattern search (PS) algorithm.For optimization problems with a larger number of design parameters, it slightly improves the solution by 3-6%.

Numerical Examples and Discussion
As an example, we optimize the 96 m tall guyed telecommunication mast that is produced by the steel industries.The industrial mast of triangular horizontal plan has four guys' clusters anchored at two points 32 and 59.5 m from the shaft.The width of shaft is 1.20 m.All the leg, bracing, and stiffener elements are round, solid steel bars.At the radii of the structural elements 0.02 m (leg), 0.0125 m (bracing), and 0.007 m (guy), the total mass of the shaft is 10410 kg, while the mass of all guys is 1230 kg.
The optimal mast scheme is sought retaining the scheme of an industrial mast but varying the topology, shape, and sizing parameters of the mast.The range of guys' clusters from 1 to 10 is explored.Constant data specific to this mast which do not vary in the optimization process are as follows: The design parameters along with their characteristics (dimensions in m, forces in N) are listed in Table 1.
We treat all the cross section characteristics as continuous variables.After the solution is obtained, the designer should round the radii of structural elements to the safe side and to the available assortment or in-stock values.Consequently, the modified mast structure will be slightly heavier.
The following genetic parameters of EA were chosen after a few numerical experiments: (vi) Mutation type: power mutation [19] (vii) Penalty type: static, proportional to the extent of constraint violation Since the EA is a stochastic algorithm whose solution generally cannot be repeated, we take 30 independent numerical experiments.In this case the median solution stabilizes and additional numerical experiments do not change the best and median results significantly.After the EA solution, the best obtained result is refined with pattern search.
Solution time depends on the number of finite element mesh nodes that in turn is determined by the number of typical sections   :  =   ⋅ 6 + 3. Thirty runs of the program package usually take from 300 to 1600 min.
The final results of mast optimization successively with 1, 2, . .., 10 guys' clusters, that is, the total mass of structure ( + in Table 2), width of shaft, and radii of all structural elements after solution with the EA and the subsequent PS, are provided in Table 2.In order to reveal the influence of PS on the quality of the result, the value of total mass obtained only using EA (  in Table 2) is also given.The mass is given in kg, while the geometrical parameters are in m.Also, some relevant data on the optimization process-the total number of objective function evaluations (  ) in 30 independent numerical experiments and the total time of computations  (, in h)-are listed.The pretension forces in the guys of different clusters are not provided in the table due to lack of space; generally, the pretension forces are highest at lowest guys and gradually diminish in higher clusters.The pretension forces (also in the form of prestress) for the masts of 2-5 guys' clusters are provided in Table 3.Some additional information on the optimal pretensioning of guys may be found also in [10].The best scheme obtained is depicted in Figure 4.
From the engineering point of view, the most relevant ones are the schemes with 2-5 guys' clusters.Figure 5 definitely shows that the total mass of the mast reduces significantly starting from 3 guys' clusters: the mass of a 3-level mast is 37% less than the mass of a 2-level mast.Introducing more than 4 guys' clusters does not produce a large effect-the mass of a 5-level mast is only 7% less.Thus, the rational scheme for the analysed masts should possess 3 to 5 guys' clusters.
A closer view of the best three results obtained for the masts with 2 to 5 guys' clusters (Table 3) reveals that the optimal parameters, particularly geometrical, are distributed in rather narrow bounds.In many cases the differences between corresponding parameters do not exceed 10%.One of these determinant geometrical parameters is the width of the shaft.For masts clamped at the foundation this parameter conditions two important factors: the bending stiffness and the magnitude of axial forces in the legs of the shaft.The axial forces in the legs diminish with the increasing distance between the legs, at the expense of increasing bending stiffness.Large bending stiffness is not acceptable, since even small shaft displacements cause significant bending moments at the foundation.Due to this, greater guy stiffness is required in the lower guys.
The results also indicate that a small number of guys' clusters determine very high prestress levels in the lower guys, for example, for 2 guys' clusters the values of prestress exceed 300 MPa.The same holds for the 5-level mast: the global stability condition does not allow for diminishing bending stiffness, and this, in turn, requires larger prestress forces in the lower guys for alignment of stresses.A large prestress level has an adverse influence on metal fatigue due to the vibrations of cables.
As to the tilt angle of bracing elements, in all of the best solutions for a small number of guys' clusters it equals 35.5 degrees-this corresponds well to the optimal values of tilt angle reported in [8].This can be explained by the fact that the determinant factor in such structures is the bending of the shaft.In larger numbers of guys' clusters the shear force which determines the stresses in bracing elements is relatively higher.
Another interesting point is the fact that one of the most stable parameters is the tilt angle of the uppermost guys: independently of the number of guys' clusters, its value is close to 50 degrees.

Conclusions
A technique for simultaneous topology/shape/sizing optimization of guyed telecommunication mast is suggested.It comprises four independent parts: the optimization algorithm, finite element linear static analysis program, meshing preprocessor, and postprocessor for evaluation of constraints according to the Eurocodes and the objective function-the total mass of the shaft and guys.With "one button click" the programs yield an optimal scheme for the mast, including all dimensions of structural elements and pretension forces in the guys.The whole optimization loop takes up to 1600 min on a typical PC and provides the designer with an initial suggestion for all elements of the mast structure, thus accelerating the design process.
Using the proposed technique, the typical 96 m steelguyed mast carrying a standard antenna cluster was optimized.The optimization of the mast with different sets of design parameters definitely shows that the most relevant schemes of the mast have three to five guys' clusters, with the optimal scheme being the mast with five guys' clusters.Introducing more guys' clusters leads to a slight increase in the mass of the mast.
Analysis of the three best solutions for each 3-to 5-guycluster mast scheme reveals some stable ranges of design parameters.The optimal number of typical sections along the height of the mast is around 57 for a 2-level mast, 73 for a 3level mast, 81 for a 4-level mast, and 85-90 for a 5-level mast.The optimal width of the mast side is around 1.18, 0.87, 0.72, and 0.67 m, correspondingly.The attachment levels of each guys' cluster and guys' anchoring distances from the shaft also are very close for all of the best solutions.
As stated in prior research [11], the comparison of the optimized 96 m tall mast with the typical industrial design shows possibilities for significant reduction in expenditures on steel.

Figure 1 :Figure 2 :
Figure 1: Typical mast with two clusters of guys.Horizontal plan of the mast.
(i) Height of the mast (ii) The geometrical scheme of the typical mast section (iii) The geometrical dimensions of auxiliary equipment and antennas (iv) Material data of structural elements and guys (Young's moduli, specific weights; material is treated as isotropic) (v) Maximum allowable stresses in elements and guys (vi) Maximum allowable deflection at the top of the mast (vii) Loading data (viii) Lower and upper limits for radii of leg, bracing, stiffener members, and guys

Figure 3 :
Figure 3: Comparison of the idealized solution with the Ansys solution.
(i) Height of the mast: 96 m (ii) Area of antennas distributed in the top 15% of mast height: 10.34 m 2 (iii) Diameter of auxiliary equipment along the whole mast height: 0.05 m (iv) Wind velocity: 25 m/sec (v) Allowable steel stresses: 2.35 × 10 8 N/m 2 (vi) Specific weight of steel: 7.70 × 10 4 N/m 3 (vii) Young's modulus of steel elements: 2.10 × 10 11 N/m 2 (viii) Young's modulus of guys: 1.50 × 10 11 N/m 2 (ix) Allowable guy stresses: 1.00 × 10 9 N/m 2 (x) Weight factor for the mass of guys: 3 (i) Population size: 200 individuals (ii) Number of generations: 200 (iii) Stopping criterion: the first met: the maximum number of generations: 200; stall generations (i.e., if the weighted average change in the objective function value over stall generations is less than function tolerance: 1 − 6): 100; objective function tolerance (i.e., cumulative change in the objective function value over stall generations is less than given function tolerance)(iv) Crossover probability: 0.8 (v) Crossover type: Laplace crossover[19]

Figure 4 :
Figure 4: The best scheme obtained for the mast.

4 Figure 5 :
Figure 5: Dependence of the mast mass on the number of guys' clusters.

Table 1 :
Set of design parameters.

Table 3 :
The main optimal parameters of the best 3 solutions for 2 to 5 guys' clusters masts (m: total mass of the mast, kg; : side of the mast, m;   : number of sections; : tilt angle of bracing elements, deg;   : attachment level of the ith guys' cluster;   : anchoring distance of the ith guys' cluster, m;   : tilt angle of the guy at the ith level, deg;   : prestress force at the ith guy level, kN;  2 : prestress at the ith guy level, MPa;  ,, : radii of the leg, bracing, and guy elements, mm).