Multi-Objective Optimization of Spatially Truss Structures Based on Node Movement

: This paper discusses the solutions for topology optimization of spatially discrete structures. The optimization objects are the structural weight and the maximum displacement. The optimization variables include structural node coordinates, and the improved MOEA (Multi-objective Evolutionary Algorithm) method is used to optimize the structure. The innovation of this study is that it breaks through the shortage of constant node position in the optimization thought of traditionally discrete structure in the “Ground Structure Approach” and uses the coordinate of the node as the optimization variable for the optimization calculation. The result is not a single one but a set of optimal solutions through the evolution (i.e., Pareto optimal solutions); on this basis, the most suitable solution can be found according to the boundary conditions or other related requirements. Using the C# language to compile the calculation program, ANSYS ﬁnite element software is used to analyze the structure, and the Pareto front surface was automatically drawn to determine the optimal layout form of the discrete structure. The analysis results show that the improved MOEA method can provide an e ﬀ ective method to solve such optimization problems.


Introduction
Optimizing a discrete structure is a multi-objective optimization problem with some constraining conditions that include not only strength issues but also weight and rigidity factors. As early as 1904, Michell [1] researched the structural design of stress constraint under a single load using an analytic method and proposed a theory of truss optimization design that was called Michell theory. In a subsequent study, Prager [2] found that the Michell structure was a more flexible structure than others. In the single loading case, the bar direction is distributed along the main stress. In 1964, Dorn et al. [3] proposed a "Ground Structure Approach", which introduced the numerical method into the optimization region for the first time. The "Ground Structure" was connected by many components, which included all loading points and supporting points. The principle of this method was that: on the basis of the ground structure model, some optimization algorithm (mathematical programming approach or criterion method) can be used to omit needless bars from the ground structure according to regulations and constraints. For example, the sectional area of a bar attaches to the zero or floor level, which is unnecessary, and ultimately, the remaining bars form the optimization topology of the structure. Then, there are the other numerical algorithms, such as the weight method, constraint method [4], simplex method [5] and gradient method [6].
In recent years, the artificial intelligence algorithm has been widely used for solving the truss optimization problem. Sun et al. [7] listed all connection types of truss structure first and then set a bound of threshold value for a section surface on the truss member. Subsequently, the section surface of the truss member was used as a design variable, which was optimized by the GA (Genetic Algorithm), and in this method, some trusses with section area that was less than the floor threshold level were omitted to achieve topological optimization. Ioannis et al. [8] used the PSO (Particle Swarm Optimization) algorithm for topological optimization of a truss structure, and this theory translated the structure optimization problem constraint with stress and displacement into an unconstrained optimization issue by the penalty function. Under this condition, the optimization problem of truss structure constraint with stress and displacement was finally solved. In recent years, Free Gradient Methods for short (FGMs ) have been applied to many industrial applications in aerospace and civil engineering [9].
However, regardless of using the traditional numerical optimization method or a new artificial intelligence optimization method, there is no fundamental change in the optimization variables. Thus, most of the topological optimization of the truss structure is based on the ground structure, and the result of the optimal solution is actually a subset of given layouts. Thus, the range of the ground structure adoption is defined from the beginning of the optimization, which leads to an approximate optimal solution and loss of the optimal solution because the true value is not included. Therefore, the real topological optimization is supposed to search optimization in any given region, which means that researchers should consider the effect of the ultimate topological form on node movement of the structure.
For optimization goals, one of the biggest features of the evolutionary algorithm search solution is that it is a group search algorithm that can deal with the parallel solving of the problem of multi-objective optimization with two or more goals. In 1967, Rosenberg proposed a solution to a multi-objective optimization problem using a genetic search algorithm in his doctoral thesis, while VEGA (Vector Evaluated Genetic Algorithm) did not appear until 1985, which was the first multi-objective evolutionary algorithm (MOEA) [10]. Then, many researchers improved this method later. Fonseca and Fleming [11] proposed multi-objective optimization GAs, Srivivas and Deb [12] put forward a non-dominated sorting GA, and these methods were called MOEA based on Pareto. San et al. [13] employed NURBS (Non-Uniform Rational B-Splines) functions to represent the geometrical shape and then optimized a class of structures under mass constraints and self-weight. Mroginski et al. [14] proposed a selective GA for the multi-objective design of trusses, and their approach involved sensitivity analysis of each variable to reduce the search space and to improve the convergence rate.
A multi-objective optimization problem is different from single objective optimization; it is not a single answer but a group of solutions, called the Pareto frontier. When there are multiple targets, there are usually solutions that cannot be easily compared. This solution is referred to as the Pareto optimal solution. In this paper, to solve the problem of discrete optimization, node movement in the structure was considered. Using a C# language self-programming optimization program, ANSYS finite element software analysis and calculations were used to draw the current Pareto optimal frontier in real time for structural optimization. Subsequently, the weight and maximum displacement of the discrete structure were used as the optimization goal, and the space truss and coordinates of the tower nodes were used as optimization variables to determine the optimal topological mode for this kind of discrete structure.

Basic Concept of MOEA
A multi-objective optimization problem is different from single objective optimization, which is not a single answer but a group of solutions, called the Pareto frontier [15]. The conception of Pareto was suggested by French economist V. Pareto. First, the Pareto frontier is defined as one or some individuals in an optimization region that had no absolute good individual that was better than any Appl. Sci. 2020, 10,1964 3 of 17 other one. "Absolutely good" means that at least one goal is good, and the other goals are not bad. Then, the multi-objective problem requests a group of equilibrium solutions, not a single optimal solution.
For a given issue, the Pareto optimal solution and frontier is definite, and an MOEA is used to solve the multi-objective optimization problem with the view of obtaining the Pareto optimal set and then confirming the final solution according to other information about the problem (such as boundary conditions, etc.). In many test functions, the Pareto frontier of the problem is known, and it can be used to test the algorithm. Nevertheless, the Pareto frontier in practical problems is unknown. (Figure 1 shows hollow circles that indicate a set of solutions, and solid square spots indicate the Pareto frontier as the optimal solution set with the boundary conditions of 0.6 < f 1 < 1.0 and 0.6 < f 2 < 1.0; the blue circle region shows the selected optimal solutions.) Appl. Sci. 2020, 10, x FOR PEER REVIEW  3 of 17   Then, the multi-objective problem requests a group of equilibrium solutions, not a single optimal  solution. For a given issue, the Pareto optimal solution and frontier is definite, and an MOEA is used to solve the multi-objective optimization problem with the view of obtaining the Pareto optimal set and then confirming the final solution according to other information about the problem (such as boundary conditions, etc.). In many test functions, the Pareto frontier of the problem is known, and it can be used to test the algorithm. Nevertheless, the Pareto frontier in practical problems is unknown. (Figure 1 shows hollow circles that indicate a set of solutions, and solid square spots indicate the Pareto frontier as the optimal solution set with the boundary conditions of 0.6 < f1 < 1.0 and 0.6 < f2 < 1.0; the blue circle region shows the selected optimal solutions.)

Improvement to MOEA
There are two main points in the improved MOEA approach of this paper: 1. It abandons the traditional optimization algorithm, which only considers a single target of structure weight in the optimization process, and adopts two optimization targets of weight and maximum displacement for control. 2. The optimization variable is not the traditional bar section, but instead uses the truss node coordinates as the variable. 3. The fitness function is not calculated by the objective function directly, but by indirect connection of the Pareto strength and the objective function value. The Pareto strength increases during the evolutionary process (shown in Figure 2); this is the individual solution with low density, and then the total disaggregation advances to the Pareto optimal frontier through an evolutionary process to finally obtain the representative optimal solution set. Because of the direct relevance of the fitness function and Pareto strength, the improved algorithm ensures the convergent solution set as a group of optimal solutions [16,17]; moreover, the density degree of the unit in the solution space is considered so that this method can maintain population diversity and uniform distribution of disaggregation and prevent the disaggregation from converging in a certain part of the Pareto optimal frontier [18,19].

Improvement to MOEA
There are two main points in the improved MOEA approach of this paper: 1. It abandons the traditional optimization algorithm, which only considers a single target of structure weight in the optimization process, and adopts two optimization targets of weight and maximum displacement for control. 2. The optimization variable is not the traditional bar section, but instead uses the truss node coordinates as the variable. 3. The fitness function is not calculated by the objective function directly, but by indirect connection of the Pareto strength and the objective function value. The Pareto strength increases during the evolutionary process (shown in Figure 2); this is the individual solution with low density, and then the total disaggregation advances to the Pareto optimal frontier through an evolutionary process to finally obtain the representative optimal solution set. Because of the direct relevance of the fitness function and Pareto strength, the improved algorithm ensures the convergent solution set as a group of optimal solutions [16,17]; moreover, the density degree of the unit in the solution space is considered so that this method can maintain population diversity and uniform distribution of disaggregation and prevent the disaggregation from converging in a certain part of the Pareto optimal frontier [18,19].

Mathematical Model
Multi-objective optimization is different from single-objective optimization; multi-objective optimization usually has multiple optimal solutions that meet the conditions. In order to explain how to deal with multiple optimal solutions, the mathematical model of multi-objective optimization problem is given.
Assuming there are m optimization objectives, they are minimized as shown in Equation (3): (3)

Mathematical Model
Multi-objective optimization is different from single-objective optimization; multi-objective optimization usually has multiple optimal solutions that meet the conditions. In order to explain how to deal with multiple optimal solutions, the mathematical model of multi-objective optimization problem is given.

Pareto Dominate Relations
When there are multiple targets, there are usually solutions that cannot be easily compared. This solution is usually the Pareto optimal solution. The definition of decision space and target space are given here.
Definition (1): definition of decision space: Definition (2): definition of target space: The objective function f(x) maps Ω to the set Π.
Definition (3): definition of dominate relations: Let x a and x b be two different solutions. If x a dominates x b , the two conditions of Equations (6) and (7) must be satisfied.
For all subgoals, x a is no worse than x b : At least there are subgoals that make xa better than xb: where m is the number of subgoals; then x a x b " "represents dominate relations. The distribution of each individual in the population is shown by a Pareto strength value, indicating the number of individual i in the evolutionary population P. The calculation formula is shown as Equation (8): Based on the Pareto strength value, the fitness value of an individual represents the sum of the strength values of the dominant individual in the evolutionary population. The calculation formula is shown as Equation (9): If there is the same original fitness value between individuals, the density information of individuals is combined into the final fitness. The density information here refers to the application of the adjacent individual strategy, that is, the distance between individual i and other individuals in the evolutionary population P, then calculated and arranged in increasing order. The calculation formula is shown as Equation (10): where δ k i where is the distance between individual i and the kth adjacent individual. When the distance between individual i and the kth adjacent individual is larger, the individual distribution is better; the smaller the fitness value, the more likely the individual is to be selected for the next evolution.
Considering these factors, the C# language was adopted to compile the MOEA optimization program based on the Pareto strength value order, that each individual in the population is distributed with a strength value, representing the number of individuals dominated by individual in the evolutionary population. The fitness value of an individual represents the sum of the strength values of the dominant individual in the evolutionary population. The smaller the fitness value, the more likely the individual is to be selected for the next evolution. The ANSYS scripting language APDL was used to write the analytical procedure for the truss finite element in this study. During the counting process, two objective functions (weight and maximum displacement of the structure) were calculated by the finite element method, and this method adopted an alternative pattern to update in real time by way of transferring the value to the optimization module through the connector.

Numerical Method Verification
Deb, K. [20] indicated that there are two difficult key problems in multi-objective optimization: 1. how to achieve optimal population converging to the Pareto frontier and 2. how to ensure specimen diversity in the population. Thus, there are various optimal solutions that can be found, such as convex or concave (shown in Figure 3) or discrete or heterogeneous distributions. For this Pareto optimal solution with such characteristics, many researchers have designed a corresponding test function, which is confined to two objectives for each one.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 6 of 17 maximum displacement of the structure) were calculated by the finite element method, and this method adopted an alternative pattern to update in real time by way of transferring the value to the optimization module through the connector.

Numerical Method Verification
Deb, K. [20] indicated that there are two difficult key problems in multi-objective optimization: 1. how to achieve optimal population converging to the Pareto frontier and 2. how to ensure specimen diversity in the population. Thus, there are various optimal solutions that can be found, such as convex or concave (shown in Figure 3) or discrete or heterogeneous distributions. For this Pareto optimal solution with such characteristics, many researchers have designed a corresponding test function, which is confined to two objectives for each one.  To test the veracity of the algorithm in this study, there is an article that lists the two normative testing functions that are able to reflect the essential feature of multi-objective problems [21]. MOP1 and MOP2 test functions can examine the algorithm and procedure, which represent two familiar conditions in the structure optimization as convex (MOP1) and concave (MOP2) of the Pareto frontier, respectively.

Function Test
MOP1 is shown as Equation (11), which has one optimization variable x, and two optimization objectives, ƒ1 and ƒ2, and this optimal problem has a convex Pareto frontier: The operating parameters of MOP1 are shown in Table 1.  To test the veracity of the algorithm in this study, there is an article that lists the two normative testing functions that are able to reflect the essential feature of multi-objective problems [21]. MOP1 and MOP2 test functions can examine the algorithm and procedure, which represent two familiar conditions in the structure optimization as convex (MOP1) and concave (MOP2) of the Pareto frontier, respectively.

Test Function MOP1
MOP1 is shown as Equation (11), which has one optimization variable x, and two optimization objectives, ƒ 1 and ƒ 2 , and this optimal problem has a convex Pareto frontier: MOP1 : The operating parameters of MOP1 are shown in Table 1.  (12)) has a concave Pareto frontier: MOP2 : The operating parameters of MOP2 are shown in Table 2.

Test Results
The numerical experiment for the improved MOEA was completed in the organizational procedure above, and the evolutionary process and final optimal consequence are shown in Figure 3. The results showed that two types of test functions (convex and concave) obtained the Pareto solution set with symmetrical distributions and good results.

Space Truss Optimization
The initial model of the truss was 3m long and was divided into 6 sections, and the physical dimensions of each component of the truss are shown in Figure 4. Moreover, the truss section is in the form of an equilateral triangle, as shown in Figure 5. According to the requirements of the truss structure, the angle between the two bars should be controlled from 30 • to 50 • . This structure adopts the same section form bar, which does not cross, and they are divided into two forms: a left oblique and a right oblique form. The boundary condition is fixed at one end and free at the other end, and the free end is subjected to a downward concentrated force and the force size is 1000 N. The truss material is steel structure, and beam189 of ANSYS software is used for analysis and calculation.
In the process of truss optimization calculation, since the Beam element is adopted in the truss, each member is divided into three elements. In the whole optimization process, it takes 10 seconds to calculate each optimization step and about 5 minutes to calculate 50 steps.  (12)) has a concave Pareto frontier: MOP2: The operating parameters of MOP2 are shown in Table 2.

Test Results
The numerical experiment for the improved MOEA was completed in the organizational procedure above, and the evolutionary process and final optimal consequence are shown in Figure  3. The results showed that two types of test functions (convex and concave) obtained the Pareto solution set with symmetrical distributions and good results.

Space Truss Optimization
The initial model of the truss was 3m long and was divided into 6 sections, and the physical dimensions of each component of the truss are shown in Figure 4. Moreover, the truss section is in the form of an equilateral triangle, as shown in Figure 5. According to the requirements of the truss structure, the angle between the two bars should be controlled from 30° to 50°. This structure adopts the same section form bar, which does not cross, and they are divided into two forms: a left oblique and a right oblique form. The boundary condition is fixed at one end and free at the other end, and the free end is subjected to a downward concentrated force and the force size is 1000 N. The truss material is steel structure, and beam189 of ANSYS software is used for analysis and calculation.
In the process of truss optimization calculation, since the Beam element is adopted in the truss, each member is divided into three elements. In the whole optimization process, it takes 10 seconds to calculate each optimization step and about 5 minutes to calculate 50 steps.  The truss nodes are considered to move along the horizontal chord (z-direction: color number direction as shown in Figure 6). The specific process is as follows: 1. We made the chord node move to the optimal location first, and then connected the diagonal bar. 2. In the optimization procedure, the truss forms were classified into two types (Figure 7): the first type: with a pitch number equal to the upper and lower chords, and the second type: unequal type: the top and bottom strings are the same length, but the top string is more knotted than the bottom string, which is staggered. 3. An "on-off" pattern was used, i.e., nodes were dispersed into four independent dots with a continuous moving range around the half-section, i.e., eight equal diversion points of the truss node (Figure 6.) so that it was moving the down cord node. This not only avoids the connection of cross-node number, but also improves the calculation efficiency and makes the optimized result modular. The truss nodes are considered to move along the horizontal chord (z-direction: color number direction as shown in Figure 6). The specific process is as follows: 1. We made the chord node move to the optimal location first, and then connected the diagonal bar. 2. In the optimization procedure, the truss forms were classified into two types (Figure 7): the first type: with a pitch number equal to the upper and lower chords, and the second type: unequal type: the top and bottom strings are the same length, but the top string is more knotted than the bottom string, which is staggered. 3. An "on-off" pattern was used, i.e., nodes were dispersed into four independent dots with a continuous moving range around the half-section, i.e., eight equal diversion points of the truss node (Figure 6.) so that it was moving the down cord node. This not only avoids the connection of cross-node number, but also improves the calculation efficiency and makes the optimized result modular. The truss nodes are considered to move along the horizontal chord (z-direction: color number direction as shown in Figure 6). The specific process is as follows: 1. We made the chord node move to the optimal location first, and then connected the diagonal bar. 2. In the optimization procedure, the truss forms were classified into two types (Figure 7): the first type: with a pitch number equal to the upper and lower chords, and the second type: unequal type: the top and bottom strings are the same length, but the top string is more knotted than the bottom string, which is staggered. 3. An "on-off" pattern was used, i.e., nodes were dispersed into four independent dots with a continuous moving range around the half-section, i.e., eight equal diversion points of the truss node ( Figure 6.) so that it was moving the down cord node. This not only avoids the connection of cross-node number, but also improves the calculation efficiency and makes the optimized result modular.   The analytical result involves optimization of the discrete variables in the space truss structures, where the variable is the z direction coordinate of the truss node and its locator variable is eight positions. The individual number of MOEA population is 100; then, the maximum sample size N is set to 30, and the final result is obtained by iterative computations, from the 5th step, 16th step, 25th step, 31st step, 42nd step, and to the 50th step (Figure 8a-f).

1.74
Step-50 This method has the following characteristics. 1. During the first few steps of evolution, there were relatively few Pareto front points, and there were often obvious inflection points. 2. As the number of iterations increased, the points on the front edge of Pareto became more evenly distributed, and the inflection point positions became more blurred. 3. As the optimization process continued, the Pareto frontier did not stop the forward leading until a satisfying result was obtained (as shown in Figure 9). All data points resulting from this calculation and the total Pareto frontier curve are shown in Figure 10. This method has the following characteristics. 1. During the first few steps of evolution, there were relatively few Pareto front points, and there were often obvious inflection points. 2. As the number of iterations increased, the points on the front edge of Pareto became more evenly distributed, and the inflection point positions became more blurred. 3. As the optimization process continued, the Pareto frontier did not stop the forward leading until a satisfying result was obtained (as shown in Figure 9). All data points resulting from this calculation and the total Pareto frontier curve are shown in Figure 10.

Displacement (mm)
This method has the following characteristics. 1. During the first few steps of evolution, there were relatively few Pareto front points, and there were often obvious inflection points. 2. As the number of iterations increased, the points on the front edge of Pareto became more evenly distributed, and the inflection point positions became more blurred. 3. As the optimization process continued, the Pareto frontier did not stop the forward leading until a satisfying result was obtained (as shown in Figure 9). All data points resulting from this calculation and the total Pareto frontier curve are shown in Figure 10. The truss shape with typical evolution characteristics is shown in Figure 7, and the result indicates that truss shape appeared disordered and had a poor regularity during the early stage of evaluation. When there were more than 30 steps in the evaluation, four small pyramid structures on the two sides with a vertical web member disappeared from the initial model, while the vertical web member on the surface was maintained. Moreover, the diagonal web member presents a topological configuration from end to end. With an increase in evolution steps, a rectangular pyramid also shows the variation trend of dimensional unification ( Figure 11 shows the truss topological graph that corresponds to the optimization procedure in Figure 8). were relatively few Pareto front points, and there were often obvious inflection points. 2. As the number of iterations increased, the points on the front edge of Pareto became more evenly distributed, and the inflection point positions became more blurred. 3. As the optimization process continued, the Pareto frontier did not stop the forward leading until a satisfying result was obtained (as shown in Figure 9). All data points resulting from this calculation and the total Pareto frontier curve are shown in Figure 10. The truss shape with typical evolution characteristics is shown in Figure 7, and the result indicates that truss shape appeared disordered and had a poor regularity during the early stage of evaluation. When there were more than 30 steps in the evaluation, four small pyramid structures on the two sides with a vertical web member disappeared from the initial model, while the vertical web member on the surface was maintained. Moreover, the diagonal web member presents a topological configuration from end to end. With an increase in evolution steps, a rectangular pyramid also shows the variation trend of dimensional unification ( Figure 11 shows the truss topological graph that corresponds to the optimization procedure in Figure 8). The truss shape with typical evolution characteristics is shown in Figure 7, and the result indicates that truss shape appeared disordered and had a poor regularity during the early stage of evaluation. When there were more than 30 steps in the evaluation, four small pyramid structures on the two sides with a vertical web member disappeared from the initial model, while the vertical web member on the surface was maintained. Moreover, the diagonal web member presents a topological configuration from end to end. With an increase in evolution steps, a rectangular pyramid also shows the variation trend of dimensional unification ( Figure 11 shows the truss topological graph that corresponds to the optimization procedure in Figure 8 In practical structural design, the chosen modal is not the top result in the optimal computation but a relatively good result according to the particular needs, such as site condition, user demand, economy, beauty, and so on. Based on the rule of truss evolution, the topological form is good, which has a side that is a symmetrical rectangular pyramid, surface as a vertical web member, and diagonal web member from end to end. According to the loading symmetry, practical engineering should choose the truss of the upper chord with more half-sections, as in Figure 12. As a result, this truss has an optimizing index that is closer to the Pareto frontier, which is a reasonable choice. In practical structural design, the chosen modal is not the top result in the optimal computation but a relatively good result according to the particular needs, such as site condition, user demand, economy, beauty, and so on. Based on the rule of truss evolution, the topological form is good, which has a side that is a symmetrical rectangular pyramid, surface as a vertical web member, and diagonal web member from end to end. According to the loading symmetry, practical engineering should choose the truss of the upper chord with more half-sections, as in Figure 12. As a result, this truss has an optimizing index that is closer to the Pareto frontier, which is a reasonable choice.  Figure 12. Suggested selection point.

Space Tower Optimization
The form for a tower structure is shown in Figure 13, whose internode length of horizontal and vertical lengths are both 3 m. There is a 10 kN horizontal load on left side node (as shown in Figure  14). Under this load, the structure was optimized and analyzed to minimize the total structure weight and the maximum structural displacement. There were three boundary conditions. (a) The four angular points of each part of the tower were defined as the control nodes, which allowed the displacement value to move in the horizontal direction from 0 to 1.5m for intra-movement; moreover, the maximum displacement distance was half the side length. After control nodes were located, other nodes were obtained by interpolating, then connected to the vertical members. (b) The coordinate of the base position did not change.

Space Tower Optimization
The form for a tower structure is shown in Figure 13, whose internode length of horizontal and vertical lengths are both 3 m. There is a 10 kN horizontal load on left side node (as shown in Figure 14). Under this load, the structure was optimized and analyzed to minimize the total structure weight and the maximum structural displacement. There were three boundary conditions. (a) The four angular points of each part of the tower were defined as the control nodes, which allowed the displacement value to move in the horizontal direction from 0 to 1.5m for intra-movement; moreover, the maximum displacement distance was half the side length. After control nodes were located, other nodes were obtained by interpolating, then connected to the vertical members. (b) The coordinate of the base position did not change.  Figure 12. Suggested selection point.

Space Tower Optimization
The form for a tower structure is shown in Figure 13, whose internode length of horizontal and vertical lengths are both 3 m. There is a 10 kN horizontal load on left side node (as shown in Figure  14). Under this load, the structure was optimized and analyzed to minimize the total structure weight and the maximum structural displacement. There were three boundary conditions. (a) The four angular points of each part of the tower were defined as the control nodes, which allowed the displacement value to move in the horizontal direction from 0 to 1.5m for intra-movement; moreover, the maximum displacement distance was half the side length. After control nodes were located, other nodes were obtained by interpolating, then connected to the vertical members. (b) The coordinate of the base position did not change.  As shown by the continuous evolution diagram of the tower frame, all the nodes were constantly moving during the early stage of the optimization calculation, but the movement was more chaotic (as shown in Figure 15). When optimal computation accessed 30 steps, the nodes close to the base are moved away from the section center, meaning the tower was sturdier. However, the closer the node was to the top, the closer to the section center, meaning the tower was thinner, which is also conformed to the load-bearing characteristics of the cantilever. When the optimal computation achieved 50 steps, the tower weight was 10,000 kg at the maximum displacement of 35 mm, compared with the first generation for the same displacement, which lost a weight of 4%. In addition, the overall configuration showed that the points between the nodes were smoother. As shown by the continuous evolution diagram of the tower frame, all the nodes were constantly moving during the early stage of the optimization calculation, but the movement was more chaotic (as shown in Figure 15). When optimal computation accessed 30 steps, the nodes close to the base are moved away from the section center, meaning the tower was sturdier. However, the closer the node was to the top, the closer to the section center, meaning the tower was thinner, which is also conformed to the load-bearing characteristics of the cantilever. When the optimal computation achieved 50 steps, the tower weight was 10,000 kg at the maximum displacement of 35 mm, compared with the first generation for the same displacement, which lost a weight of 4%. In addition, the overall configuration showed that the points between the nodes were smoother. Appl. Sci. 2020, 10, x FOR PEER REVIEW 14 of 17 Step-50 evolution diagram. Figure 16 shows the Pareto frontier curve drawn using the optimization step corresponding to Figure 15. With continuous deepening of evolution, the data points became smoother and more uniform, and the optimization effect became more obvious; thus, the method is feasible, and the results are credible.  Figure 16 shows the Pareto frontier curve drawn using the optimization step corresponding to Figure 15. With continuous deepening of evolution, the data points became smoother and more uniform, and the optimization effect became more obvious; thus, the method is feasible, and the results are credible.

Conclusions
(1) In this study, an improved MOEA optimization method was adopted to realize optimization of spatial discrete structures using moving bar node coordinates.
(2) In contrast to the traditional optimization, which uses weight as the only object, the improved MOEA multi-objective optimization method is a combination of the structural weight and the maximum displacement of the spatially discrete structure. This method uses the strength rank idea and obtains a group of optimal solutions, not an individual value. Then, the optimal case can be confirmed according to other constraint conditions, such as the boundary value and the user's preferences.
(3) This algorithm was tested using standard mathematical functions, and the results indicate that the rate of convergence is fast and that no inferior solution evenly distributes in the solution

Conclusions
(1) In this study, an improved MOEA optimization method was adopted to realize optimization of spatial discrete structures using moving bar node coordinates. (2) In contrast to the traditional optimization, which uses weight as the only object, the improved MOEA multi-objective optimization method is a combination of the structural weight and the maximum displacement of the spatially discrete structure. This method uses the strength rank idea and obtains a group of optimal solutions, not an individual value. Then, the optimal case can be confirmed according to other constraint conditions, such as the boundary value and the user's preferences.
(3) This algorithm was tested using standard mathematical functions, and the results indicate that the rate of convergence is fast and that no inferior solution evenly distributes in the solution space, which maintains good diversity without converging to a local part. Thus, we used the combination of this method with ANSYS to analyze a three-dimensional truss and tower structure, and the results confirmed that the modified MOEA algorithm appears to be feasible and superior with respect to the optimization of the discrete structure design.

Conflicts of Interest:
The authors declare no conflict of interest.