Ground Structures-Based Topology Optimization of a Morphing Wing Using a Metaheuristic Algorithm

: This paper presents multi-objective topology and sizing optimization of a morphing wing structure. The purpose of this paper is to design a new aircraft wing structure with a tapered shape for ribs, spars, and skins including a torsion beam for external actuating torques, which is anticipated to modify the aeroelastic characteristic of the aircraft wing using multi-objective optimization. Two multi-objective topology optimization problems are proposed employing ground element structures with high- and low-grid resolutions. The design problem is to minimize mass, maximize difference of lift effectiveness, and maximize the buckling factor of an aircraft wing subject to aeroelastic and structural constraints including lift effectiveness, critical speed, and buckling factors. The design variables include aircraft wing structure dimensions and thickness distribution. The proposed optimization problems are solved by an efﬁcient multi-objective metaheuristic algorithm while the results are compared and discussed. The Pareto optimal fronts obtained for all tests were compared based on a hypervolume metric. The objective function values for Case I and Case II at 10 selected optimal solutions exhibit a range of structural mass as 115.3216–411.6250 kg, 125.0137–440.5869 kg, lift effectiveness as 1.0514–1.1451, 1.0834–1.1639 and bucking factor as 38.895–1133.1864 Hz, 158.1264–1844.4355 Hz, respectively. The best results reveal unconventional aircraft wing structures that can be manufactured using additive manufacturing. This research is expected to serve as a foundation for future research into multi-objective topology optimization of morphing wing structures based on the ground element framework. Two multi-objective topology optimization investigated by employing ground with high- and low-grid resolutions. The design problem has the objectives of weight reduction, difference of lift effectiveness maximization, and buckling factor maximization of an aircraft wing subject to aeroelastic and structural constraints including lift effectiveness, critical speed, and buck- ling factors. The design variables include aircraft wing structure dimensions and thickness distribution and internal skeletal rod sizing and topology. The MM-IPDE optimizer is used to solve the problems where the optimum results show that the topology optimiza- tion problem having higher ground structure resolution gives the higher range of the ob- jective function values. Furthermore, these ground structure designs with high resolution accomplish superior buckling factor and lift effectiveness and slightly better maximum transverse displacement at the same structural mass while the lower ground structure resolution gives the better critical speed at the same structural mass. Moreover, the higher ground structure resolution problem offers more diverse optimal solutions and topologies of the wing structure. In case of testing, the new structure is acceptable for use as a morph- ing aircraft structure. The result shows that the new structure can change the wing box shape and lift effectiveness upon the applied torques. The new structure concept is said to be acceptable to be used as a morphing aircraft wing. resolution of the ground structure for the morphing wing on the resulting optimal solutions is investigated. Two multi-objective topology optimization problems are proposed and investigated by employing ground element structures with high- and low-grid resolutions. The design problem has the objectives of weight reduction, difference of lift effectiveness maximization, and buckling factor maximization of an aircraft wing subject to aeroelastic and structural constraints including lift effectiveness, critical speed, and buckling factors. The design variables include aircraft wing structure dimensions and thickness distribution and internal skeletal rod sizing and topology. The MM-IPDE optimizer is used to solve the problems where the optimum results show that the topology optimization problem having higher ground structure resolution gives the higher range of the objective function values. Furthermore, these ground structure designs with high resolution accomplish superior buckling factor and lift effectiveness and slightly better maximum transverse displacement at the same structural mass while the lower ground structure resolution gives the better critical speed at the same structural mass. Moreover, the higher ground structure resolution problem offers more diverse optimal solutions and topologies of the wing structure. In case of testing, the new structure is acceptable for use as a morphing aircraft structure. The result shows that the new structure can change the wing box shape and lift effectiveness upon the applied torques. The new structure concept is said to be acceptable to be used as a morphing aircraft wing.


Introduction
Weight reduction in conjunction with increased structural and aero-elastic performances such as flutter speed, buckling is vital for competitiveness in the aircraft industry. Therefore, the design and development process of aircraft is essential to boost the potential of aircraft wing structures. However, with conventional wing design standards that pose limitations of fixed structure and control surfaces, it is very difficult to achieve maximum performance.
Hence, to reach the highest performance of the aircraft over the flight envelope, the notion of morphing aircraft is becoming interesting for researchers. This is because adjusting aircraft structures can improve flying performance, control authority, fuel consumption, and multi-mission capabilities to some extent [1,2]. However, due to the limitation of cost structures can be used in practice. The unconventional structural design looks possible for manufacturing [22] while it has repeatability to produce [30].
As mentioned above, the key of the morphing wings is the internal structure which can be classified as a compliant mechanism [17,20,21]. The internal structure of morphing wings needs to be optimized to meet the best possible wing performance while the optimization problem can be both single and multi-objective (MO) to minimize the structural mass and/or maximizing aircraft tall flight performance subjected to several constraints such as flutter speed, flexibility, buckling factor, etc., [31]. The topology and size of the structure are the most popular design variables considered so far in the literature. Over the past decade, the TO of morphing wings has been successfully presented based on both concepts of ground structures [17] and ground elements [20] topology.
Optimization techniques based on metaheuristic or evolutionary algorithms, on the other hand, have also been successfully applied for both single-objective [32,33] and MO problems [34,35]. Although numerous investigations on TO of morphing wings have been studied, most of the work focuses on investigating the performance of an aeroelasticity analysis tool, an optimization solver tool, and an optimum structure. It was found in the literature that there is rarely a work-study on the effect of the resolution of the ground structure or ground element on the design results and computational time.
Therefore, this work concentrates on this research gap, aiming to investigate the effect of ground structure topology resolution on wing morphology while also presenting a novel design for a morphed wing structure. The innovative morphed wings investigated employ a torsion beam for actuating and skeletal elements as internal structure instead of conventional wing ribs and inside spars as shown in Figure 1. The optimization problem is posed to minimize wings mass, maximize difference of lift effectiveness, and maximize buckling factor subjected to aeroelastic and structural constraints including lift effectiveness, critical speed, and buckling factors. The design variables are thickness for ribs, spars, and skins, distribution thickness function decision, and topology and sizes of skeletal elements. Two MO problems are proposed based on skeletal ground structures with high-and low-grid resolutions. The proposed optimization problems are solved by multi-objective metaheuristic with iterative parameter distribution estimation (MM-IPDE) while the results obtained from the skeletal ground element structures with high-and low-grid resolutions are compared and discussed. The remaining of this paper is divided into five sections. Section 2 shows the details of the aircraft wing model and its aeroelastic analysis while the details of the experimental setup are presented in Section 3. The results and discussion, and conclusions are detailed in Sections 4 and 5, respectively. Therefore, the design results are proved that these peculiar structures can be used in practice. The unconventional structural design looks possible for manufacturing [22] while it has repeatability to produce [30]. As mentioned above, the key of the morphing wings is the internal structure which can be classified as a compliant mechanism [17,20,21]. The internal structure of morphing wings needs to be optimized to meet the best possible wing performance while the optimization problem can be both single and multi-objective (MO) to minimize the structural mass and/or maximizing aircraft tall flight performance subjected to several constraints such as flutter speed, flexibility, buckling factor, etc., [31]. The topology and size of the structure are the most popular design variables considered so far in the literature. Over the past decade, the TO of morphing wings has been successfully presented based on both concepts of ground structures [17] and ground elements [20] topology.
Optimization techniques based on metaheuristic or evolutionary algorithms, on the other hand, have also been successfully applied for both single-objective [32,33] and MO problems [34,35]. Although numerous investigations on TO of morphing wings have been studied, most of the work focuses on investigating the performance of an aeroelasticity analysis tool, an optimization solver tool, and an optimum structure. It was found in the literature that there is rarely a work-study on the effect of the resolution of the ground structure or ground element on the design results and computational time.
Therefore, this work concentrates on this research gap, aiming to investigate the effect of ground structure topology resolution on wing morphology while also presenting a novel design for a morphed wing structure. The innovative morphed wings investigated employ a torsion beam for actuating and skeletal elements as internal structure instead of conventional wing ribs and inside spars as shown in Figure 1. The optimization problem is posed to minimize wings mass, maximize difference of lift effectiveness, and maximize buckling factor subjected to aeroelastic and structural constraints including lift effectiveness, critical speed, and buckling factors. The design variables are thickness for ribs, spars, and skins, distribution thickness function decision, and topology and sizes of skeletal elements. Two MO problems are proposed based on skeletal ground structures with highand low-grid resolutions. The proposed optimization problems are solved by multi-objective metaheuristic with iterative parameter distribution estimation (MM-IPDE) while the results obtained from the skeletal ground element structures with high-and low-grid resolutions are compared and discussed. The remaining of this paper is divided into five sections. Section 2 shows the details of the aircraft wing model and its aeroelastic analysis while the details of the experimental setup are presented in Section 3. The results and discussion, and conclusions are detailed in Sections 4 and 5, respectively.

Aircraft Wing Model
In this work, a new concept of morphing wing structure is proposed via employing skeletal elements as part of the internal structure instead of conventional wing rib and

Aircraft Wing Model
In this work, a new concept of morphing wing structure is proposed via employing skeletal elements as part of the internal structure instead of conventional wing rib and spar while the actuator is modelled as a torsion beam. The Goland wing is a rectangular-shaped frame that is comprised of leading, trailing, and center edge spar including11 ribs as shown  Figure 2a. It is selected for the present case study and is modelled as an Aluminum wing box (7075 Aluminum alloy) while the material properties of the wing are shown in Table 1. The chord length and semi-span of the wing are 1.216 m and 6.096 m., respectively while the wing thickness is 0.0508 m. The comparison of the internal wing structure of the original Goland wing and the proposed internal concept structure are shown the Figure 2. spar while the actuator is modelled as a torsion beam. The Goland wing is a rectangularshaped frame that is comprised of leading, trailing, and center edge spar including11 ribs as shown in Figure 2a. It is selected for the present case study and is modelled as an Aluminum wing box (7075 Aluminum alloy) while the material properties of the wing are shown in Table 1. The chord length and semi-span of the wing are 1.216 m and 6.096 m., respectively while the wing thickness is 0.0508 m. The comparison of the internal wing structure of the original Goland wing and the proposed internal concept structure are shown the Figure 2.

Aerodynamic Model
In this paper, a quasi-steady vortex ring method [17] is used for aerodynamic analysis. The aerodynamic forces acting on the structure can be represented as follows: where k = Lω V is the reduced frequency, L is a semi chord length, V is wind velocity, ω is a circulation frequency, q is dynamic pressure, [S] is the diagonal matrix of panel areas, [AIC] is the aerodynamic influence coefficient matrix, [G] is the transformation matrix, and {u} is a vector of structural displacements of the wing finite element model. It should be noted that more accurate aeroelastic analysis can be obtained by using the unsteady aerodynamic analysis such as the doublet lattice method, however, it has been shown in our previous study that the quasi-steady vortex ring method gives acceptable results.

Aeroelasticity Analysis
The wing aeroelasticity model considered herein is established by the mutual interaction of three forces i.e., inertial, elastic, and aerodynamic forces. The governing equation can be written as: where [M], [D], [K] depicts the mass, damping, and stiffness matrices of the wing structure.
[Ad] and [Ak] are aerodynamic damping and stiffness matrices respectively. The latter two matrices are the cause of fluid/structure interaction, which will modify the system damping and stiffness while they are dependent on velocity.

Aerodynamic Model
In this paper, a quasi-steady vortex ring method [17] is used for aerodynamic analysis. The aerodynamic forces acting on the structure can be represented as follows: where k = Lω V is the reduced frequency, L is a semi chord length, V is wind velocity, ω is a circulation frequency, q is dynamic pressure, [S] is the diagonal matrix of panel areas, [AIC] is the aerodynamic influence coefficient matrix, [G] is the transformation matrix, and {u} is a vector of structural displacements of the wing finite element model. It should be noted that more accurate aeroelastic analysis can be obtained by using the unsteady aerodynamic analysis such as the doublet lattice method, however, it has been shown in our previous study that the quasi-steady vortex ring method gives acceptable results.

Aeroelasticity Analysis
The wing aeroelasticity model considered herein is established by the mutual interaction of three forces i.e., inertial, elastic, and aerodynamic forces. The governing equation can be written as: depicts the mass, damping, and stiffness matrices of the wing structure.
[A d ] and [A k ] are aerodynamic damping and stiffness matrices respectively. The latter two matrices are the cause of fluid/structure interaction, which will modify the system damping and stiffness while they are dependent on velocity.

Flutter Analysis
The flutter speed is a speed at which the dynamic system becomes unstable. Such a speed can be determined by sweeping the value of wind velocity V from lower to higher. With a given wind speed, a state-space matrix can be formed while the real parts of its eigenvalues are negative if the wing aeroelastic system is stable. The flutter speed is found when one of the real parts of the eigenvalues is zero.

Divergence Analysis
The divergence is usually represented as a speed when the aerodynamic load overcomes the structural restoration, considered at the steady-state condition. The governing equation can be obtained from Equation (2) by removing the terms related to time written as The system leads to an eigenvalue problem. Since the aerodynamic stiffness matrix can be singular, it is more useful to modify the problem as where λ = 1 q . Once the eigenvalue problem (4) is solved, the largest value of λ will give the lowest divergence speed.

Lift Effectiveness
The lift effectiveness is defined as the ratio of lift force on the flexible structure to its rigid counterpart, therefore, can be computed as: where q = 1 2 ρ air V 2 is the dynamic pressure, ρ air is the air density, S is the matrix of panel areas, and α is the vector of the panels' angles of attack.
[AIC] F is the flexible surface aerodynamic influence coefficient matrix. It is alternatively called the ratio of cruise shape lift to jig shape lift.

Buckling Analyses
Buckling due to aerodynamic loads is crucial for aircraft wing design as most parts of the wing are plate-like. In this work, linear buckling analysis is employed. The buckling phenomenon takes place when the work done by in-plane stress due to bending displacement on the wing overcomes its elastic potential energy. Such a situation leads to an eigenvalue problem. [ where λ is a buckling factor, [K] depicts the structural stiffness matrix, and [K G ] illustrates a global geometrical matrix.

The Proposed Multiobjective Ground Structure Topology Optimization Problem
Two cases of MO ground structure TO design problem of the newly proposed concept of Goland wing are presented based on low (Case I) and high (Case II) grid resolutions. Moreover, 18 and 54 pieces of skeletal rods distributed in the wing internal structure are considered to be a ground structure for Case I and Case II respectively as demonstrated in Figures 3 and 4.
The MO optimization problem is posed to minimize structural mass, maximize buckling factor, and maximize the difference of lift effectiveness due to counterclockwise and clockwise twists of the beam subject to aeroelasticity constraints while the design variables are sizing of ribs, spars, skins and torsion beam, distribution thickness function decision, topology, and sizing of the skeletal rods as presented in Figure 3. The first objective function is set to have a low inertia structure while the second objective is set to enhance the structural static stability. The third objective is the difference between the wing lift effective-Metals 2021, 11, 1311 6 of 17 ness due to counterclockwise and clockwise twist on the beam. This indicator determines the ability to control lift force on the wing through control actuation from twisting the beam. Thus, this wing works with the concept of morphing wing. The investigated MO optimization problem can be expressed as follows: subject to  The MO optimization problem is posed to minimize structural mass, maximize buckling factor, and maximize the difference of lift effectiveness due to counterclockwise and clockwise twists of the beam subject to aeroelasticity constraints while the design variables are sizing of ribs, spars, skins and torsion beam, distribution thickness function decision, topology, and sizing of the skeletal rods as presented in Figure 3. The first objective function is set to have a low inertia structure while the second objective is set to enhance  The MO optimization problem is posed to minimize structural mass, maximize buckling factor, and maximize the difference of lift effectiveness due to counterclockwise and clockwise twists of the beam subject to aeroelasticity constraints while the design variables are sizing of ribs, spars, skins and torsion beam, distribution thickness function decision, topology, and sizing of the skeletal rods as presented in Figure 3. The first objective The above expression includes numerous variables where x is a vector of design variables having lower and upper bounds as x l and x u followed by η L is the wing lift effectiveness. η L,al,min is the minimum allowable value of lift effectiveness. η L,al,max is the maximum allowable value of lift effectiveness. V f is a critical wind speed herein is a flutter speed. V f ,al is an allowable value of critical wind speed. λ al is the permissible magnitude of buckling factor. λ is the minimum magnitude of buckling factor. x 1 is the diameter of the torsion beam. x 2 is the distribution thickness function of the root front spar. x 3 is the thickness of the root front spar. x 4 is the thickness of the tip front spar. x 5 is the distribution thickness function of the root rear spar. x 6 is the thickness of the root rear spar. x 7 is the thickness of the tip rear spar. x 8 is the distribution function of ribs. x 9 is the thickness of leading tip ribs. x 10 is the thickness of tailing tip ribs. x 11 is the distribution thickness function of lower skin. x 12 is the thickness of root lower skin. Note that when the decision variable for beam diameter is equal to zero it means that it will be deleted leading to the variation of wing topology. For the elements of {x 2, x 5, x 8, x 11, x 14 } as decision variables for thickness distribution values, one of the six distribution functions cases (as shown in Figure 5) is selected.

Numerical Simulation Model
The new design concept of the aircraft wing structure is demonstrated using the combination of MATLAB R2016b and ANSYS 2020 R2 Student Version software. The optimi-

Numerical Simulation Model
The new design concept of the aircraft wing structure is demonstrated using the combination of MATLAB R2016b and ANSYS 2020 R2 Student Version software. The optimization process consists of running the MM-IPDE algorithm while function evaluations are achieved using finite element analysis (FEA) for static, buckling, and modal analyses. Aerodynamic loads are computed using the vortex ring method (equivalent to the vortex lattice method) as MATLAB codes. During the flight, the wing is subject to static aerodynamic loads at cruise speed, which means the stress and buckling constraints due to such applied loads are taken into consideration. The static aeroelastic phenomena, divergence speed, and lift effectiveness are computed in MATLAB by extracting a structural global stiffness matrix from ANSYS. Meanwhile, flutter analysis is carried out using a quasi-steady approach meaning that the effect of unsteady aerodynamic is excluded. It should be noted that this work is an initial study. If the more realistic wing is to be synthesized, the more accurate unsteady aerodynamic such as the doublet lattice method should be employed. For FEA with ANSYS, the shell element SHELL181 is used for wing skins, tip chord, and spars whereas all of the skeletal parts are modelled with the three-node line element BEAM188. The element size for automatic meshing is set to be 0.5 mm. Material properties of the aluminum alloy materials are those embedded in the ANSYS Library. The boundary condition set for the wing structure simulation is shown in Figure 6. A ± 1000 N-m torque is applied at the right end surface of the torsion beam while the left end of the wing surface is kept fixed. Aerodynamic force calculated based on the air density (ρ air = 1.2 kg/m 3 ) and free stream velocity (=40 m/s) are applied over the wing surface.

Result and Discussion
The new aircraft wing structural optimization solves the problem using a computer with the following specification: AMD Ryzen 7 3700X with Radeon Graphics 3.00 GHz, 32.00 GB, 64-bit Window 10 operating system. Having performed the optimization runs of the proposed MOTO challenge based on the MM-IPDE, the hypervolume indicator is used to measure the quality of the Pareto front. Figure 8 shows the search history based on the hypervolume values for both cases of low ( Figure 8a) and high (Figure 8b) grid resolutions. The figure indicates that the results converged after 100 and 80 iterations for the design Case I and Case II respectively. This implies that the problem with a higher number of design variables or high resolution of ground structure has a faster convergence rate.

Result and Discussion
The new aircraft wing structural optimization solves the problem using a computer with the following specification: AMD Ryzen 7 3700X with Radeon Graphics 3.00 GHz, 32.00 GB, 64-bit Window 10 operating system. Having performed the optimization runs of the proposed MOTO challenge based on the MM-IPDE, the hypervolume indicator is used to measure the quality of the Pareto front. Figure 8 shows the search history based on the hypervolume values for both cases of low ( Figure 8a) and high (Figure 8b) grid resolutions. The figure indicates that the results converged after 100 and 80 iterations for the design Case I and Case II respectively. This implies that the problem with a higher number of design variables or high resolution of ground structure has a faster convergence rate.         For the maximum transverse displacement, constraint minimum and maximum values are 0.0032751 and 0.014869 mm. for Case. I. while they are 0.0043311 and 0.013977 mm. for Case. II. Figure 10 demonstrates the plots of structural mass versus buckling, lift effectiveness, maximum transverse displacement of the 10 selected solutions for both cases. According to the figure, at the same structural mass, Case II realizes superior buckling and lift effectiveness with slightly better maximum transverse displacement than Case I, whereas Case I exhibited greater critical speed.   Figure 10 demonstrates the plots of structural mass versus buckling, lift effectiveness, maximum transverse displacement of the 10 selected solutions for both cases. According to the figure, at the same structural mass, Case II realizes superior buckling and lift effectiveness with slightly better maximum transverse displacement than Case I, whereas Case I exhibited greater critical speed. With a color bar, Figures 11 and 12 illustrate the topologies and sizes of solutions no. 1 and no. 10 of the 10 selected optimal locations. For both cases, solution no. 1 indicates the points where the least structural mass was obtained, whereas solution no. 10 corresponds to the points where the largest buckling factor was achieved. With a color bar, Figures 11 and 12 illustrate the topologies and sizes of solutions no. 1 and no. 10 of the 10 selected optimal locations. For both cases, solution no. 1 indicates the points where the least structural mass was obtained, whereas solution no. 10 corresponds to the points where the largest buckling factor was achieved.
Moreover, for both solutions, it is also evident from Figure 11 that the skeletal rods for no. 10 are larger than that of no. 1, with considerable variation in thicknesses and skeletal rod diameters at the left end (supported) and minor variation at the right end. Figure 12 depicts the topology altering in Case II against the ground structure of Solution 1, which has the lowest structural mass and buckling factor. The topology of the structure for a solution no. 10, which has the maximum structural mass and buckling factor, is comparable to its ground structure. The variation in thicknesses and skeletal rod diameters are high at the left end (supported) and low at the right end for solution no. 1. For solution no. 10, the thicknesses and skeletal rod diameters have the same size with a larger size than that of no.1.
Overall, it was discovered that the TO issue with better grid resolution acquired a wider range of objective function values when a single MO optimization was performed. Moreover, the higher grid resolution offers the more diverse optimal solutions and topologies of the structure which will be selected through decision-making for the next design steps.
FOR PEER REVIEW 14 of 19 Moreover, for both solutions, it is also evident from Figure 11 that the skeletal rods for no. 10 are larger than that of no. 1, with considerable variation in thicknesses and skeletal rod diameters at the left end (supported) and minor variation at the right end. Figure 12 depicts the topology altering in Case II against the ground structure of Solution 1, which has the lowest structural mass and buckling factor. The topology of the structure for a solution no. 10, which has the maximum structural mass and buckling factor, is comparable to its ground structure. The variation in thicknesses and skeletal rod diameters are high at the left end (supported) and low at the right end for solution no. 1. For solution no. 10, the thicknesses and skeletal rod diameters have the same size with a larger size than that of no.1.
Overall, it was discovered that the TO issue with better grid resolution acquired a wider range of objective function values when a single MO optimization was performed. Moreover, the higher grid resolution offers the more diverse optimal solutions and topologies of the structure which will be selected through decision-making for the next design steps.  One optimal design solution is chosen to analyze the wing capability of being morphed. The comparison between jig and cruise shapes of the wing with various values of actuating torque are given in Table 4. It can be seen that when the twisting torque is in a counterclockwise direction, the lift effectiveness becomes higher as higher torque is applied. The wing's front view shows that the cruise shape is bent upward. On the other hand, twisting the beam in a clockwise direction leads to lower wing lift effectiveness while the wing is bent downward from its front view. The maximum range of the lift effectiveness value is [0.98322, 1.1957]. This implies that the present morphing wing concept has a potential to vary the lift force in the interval of [0.98322 LR, 1.1957 LR] where LR is the total lift force on the wing jig shape at cruise speed. That means aircraft lateral/directional motion can be controlled through structural flexibility without using wing ailerons. However, more investigation in various aspects is required before it can lead to a prototype. One optimal design solution is chosen to analyze the wing capability of being morphed. The comparison between jig and cruise shapes of the wing with various values of actuating torque are given in Table 4. It can be seen that when the twisting torque is in a counterclockwise direction, the lift effectiveness becomes higher as higher torque is applied. The wing's front view shows that the cruise shape is bent upward. On the other hand, twisting the beam in a clockwise direction leads to lower wing lift effectiveness while the wing is bent downward from its front view. The maximum range of the lift effectiveness value is [0.98322, 1.1957]. This implies that the present morphing wing concept has a potential to vary the lift force in the interval of [0.98322 L R , 1.1957 L R ] where L R is the total lift force on the wing jig shape at cruise speed. That means aircraft lateral/directional motion can be controlled through structural flexibility without using wing ailerons. However, more investigation in various aspects is required before it can lead to a prototype.

Conclusions
In this work, a new concept of morphing wing structure design is successfully presented based on ground structure topology optimization. Moreover, the effect of the resolution of the ground structure for the morphing wing on the resulting optimal solutions is investigated. Two multi-objective topology optimization problems are proposed and investigated by employing ground element structures with high-and low-grid resolu-

Conclusions
In this work, a new concept of morphing wing structure design is successfully presented based on ground structure topology optimization. Moreover, the effect of the resolution of the ground structure for the morphing wing on the resulting optimal solutions is investigated. Two multi-objective topology optimization problems are proposed and investigated by employing ground element structures with high-and low-grid resolutions. The design problem has the objectives of weight reduction, difference of lift effectiveness maximization, and buckling factor maximization of an aircraft wing subject to 0.98322 −1000

Conclusions
In this work, a new concept of morphing wing structure design is successfully presented based on ground structure topology optimization. Moreover, the effect of the resolution of the ground structure for the morphing wing on the resulting optimal solutions is investigated. Two multi-objective topology optimization problems are proposed and investigated by employing ground element structures with high-and low-grid resolutions. The design problem has the objectives of weight reduction, difference of lift effectiveness maximization, and buckling factor maximization of an aircraft wing subject to aeroelastic and structural constraints including lift effectiveness, critical speed, and buckling factors. The design variables include aircraft wing structure dimensions and thickness distribution and internal skeletal rod sizing and topology. The MM-IPDE optimizer is used to solve the problems where the optimum results show that the topology optimization problem having higher ground structure resolution gives the higher range of the objective function values. Furthermore, these ground structure designs with high resolution accomplish superior buckling factor and lift effectiveness and slightly better maximum transverse displacement at the same structural mass while the lower ground structure resolution gives the better critical speed at the same structural mass. Moreover, the higher ground structure resolution problem offers more diverse optimal solutions and topologies of the wing structure. In case of testing, the new structure is acceptable for use as a morphing aircraft structure. The result shows that the new structure can change the wing box shape and lift effectiveness upon the applied torques. The new structure concept is said to be acceptable to be used as a morphing aircraft wing.

Data Availability Statement:
The authors confirm that the data supporting the findings of this study are available within the article.