Three-Dimensional CST Parameterization Method Applied in Aircraft Aeroelastic Analysis

Class/shape transformation (CST) method has advantages of adjustable design variables and powerful parametric geometric shape design ability and has been widely used in aerodynamic design and optimization processes.Three-dimensional CST is an extension for complex aircraft and can generate diverse three-dimensional aircraft and the corresponding mesh automatically and quickly. This paper proposes a parametric structural modeling method based on gridding feature extraction from the aerodynamic mesh generated by the three-dimensional CST method.This novel method can create parametric structural model for fuselage and wing and keep the coordination between the aerodynamicmesh and the structuralmesh. Based on the generated aerodynamicmodel and structural model, an automatic process for aeroelastic modeling and solving is presented with the panel method for aerodynamic solver andNASTRAN for structural solver. A reusable launch vehicle (RLV) is used to illustrate the process for aeroelastic modeling and solving. The result shows that this method can generate aeroelastic model for diverse complex three-dimensional aircraft automatically and reduce the difficulty of aeroelastic analysis dramatically. It provides an effective approach to make use of the aeroelastic analysis at the conceptual design phase for modern aircraft.


Introduction
In the wake of requirements for high lift-drag ratio aerodynamic shape and light-weight structure, the aeroelastic phenomena caused by interaction between fluid and structure have a growing influence on the integrated performance of modern aircraft [1,2]. Aeroelastic analysis becomes an important process for modern aircraft design [3,4]. Especially in the conceptual design phase, main performance of an aircraft is determined during this phase and, therefore, how to carry out the aeroelastic analysis quickly and stably to improve the integrated performance and the design rationality of the aircraft scheme will help a lot in the following aircraft design.
Aeroelastic analysis is given attention by lots of researchers [5]. Many tools such as ZAERO [6], ENSAERO [7], and NeoCASS [8] are developed to perform aeroelastic analysis based on the frequency domain analysis method and time domain analysis method and have been widely applied on high aspect ratio wing, unmanned aerial vehicle, and hypersonic aircraft. Although there exist some mature aeroelastic analysis and solving method, the aerodynamic modeling and structural modeling are still a complex and time-consuming process. First, the parametric geometric shape and modeling tools should work together closely and automatically; second, the aerodynamic solver and structural solver should coordinate with each other to ensure that the aeroelastic analysis procedure is executed consistently and accurately. In the conceptual design phase, the design scheme of aircraft usually needs constant modification to improve performance. The size parameters and structural topology of the scheme need to change frequently. It is difficult for the traditional CAD-based geometry modeling method to satisfy the need of rapid geometry iteration and large range modification. Aeroelastic analysis applied in the conceptual design phase faces the following problems: (1) The aerodynamic model and structural model established in many aeroelastic literatures were complicated; the sizing of the aerodynamic shape and structural layout is difficult and time-consuming, which makes it hard to meet the demand for rapid modification at the conceptual design phase.
2 International Journal of Aerospace Engineering (2) Because of the differences between aerodynamic modeling and structural modeling in mechanism, the aerodynamic mesh and structural mesh are not compatible. The data conversion is used additionally to interchange aerodynamic force and structural deformation between the aerodynamic model and the structural model. Extra computation overhead is needed, and some errors are introduced inevitably.
(3) Quantity and position of the structural inner elements, such as beam, bulkhead, spar, and rib, are hard to be modified using the traditional CAD-based geometry modeling method. The structural layout is vested and cannot obtain the optimal scheme by topology optimization.
These problems hinder the application of aeroelastic analysis in the conceptual design phase. More effective aerodynamic modeling and structural modeling method should be studied to simplify and improve the aeroelastic modeling and analysis process. The CST method is an analytical method developed by Kulfan [9,10]. It combines a class function representing a specific class of shapes and a shape function defining the deviation from the class function. These provide an efficient shape parameterization ability on complex geometry using fewer design variables. There are some comparisons with other shape parameterization methods [11,12], showing that the CST method has advantages in smoothness, mathematical efficiency, fitting performances, and Intuitiveness. CST has been widely used to parameterize and optimize two-dimensional airfoil [13,14] and simple 3D aircraft [15,16]. Liu [17,18] presents a multiblock CST method to model the hypersonic aircraft with complex shape, which can join adjacent surfaces smoothly and retain the good properties of the CST method. Leal [19] proposes an aerostructural optimization method for determining in a preliminary manner morphing wing configurations that provide benefits during various disparate flight conditions with CST parameterization.
The authors in [20,21] proposed a universal threedimensional CST method for geometry modeling of complex aircraft. It generates complex three-dimensional geometric shape to support various aircraft aerodynamic shape modeling, which gives a simple and effective way to aerodynamic optimization. This novel three-dimensional CST method is extended in aeroelastic analysis of complex aircraft in this article. A novel parametric structural modeling method is proposed based on gridding feature extraction from the aerodynamic mesh generated by the three-dimensional CST method. An automatic and effective way for aeroelastic modeling and analysis is also established. The structure of this article is as follows. First, the basic principle of the three-dimensional CST method is introduced briefly; then the aircraft characteristic components library is established, including fuselage, wing, and empennage; on this basis, the structural modeling method is presented in detail, and the aeroelastic modeling and analysis process is constructed; finally, a static aeroelastic analysis example is used to verify the proposed aeroelastic modeling and analysis process.

Three-Dimensional CST Method
CST method has an efficient and brief shape description for the two-dimensional airfoils and simple three-dimensional geometries. This section gives an improvement to expand the original two-dimensional CST method to three dimensions. A universal three-dimensional CST method is proposed to provide adjustable parameterization ability to complex threedimensional geometry.

Basic Cross Section
Definition. Most of the aircraft geometry can be described by innumerable cross sections along the axial direction. Combining an analytical function (the class function) with a parametric curve (the shape function), the basic cross section can be defined by the B-spline CST method [22] as follows: where = / . is the lateral coordinate. is the normalized lateral coordinate.
is the total length of the basic cross section in the lateral direction. ( ) is the normal height in the lateral position. Δ ( ) is the eccentric distance in the lateral position. ( ) and ( ) are the class function and the shape function defined in CST method. 1 and 2 are control parameters of the class function. For the symmetrical section, 1 is equal to 2. The class function ( ) is defined as follows: The basic cross section is assembled by the upper ( ) and the inverted lower ( ). If we set Δ ( ) to 0, ignore the shape function, and change the control parameters of the class function simultaneously, different sections can be generated. Various cross section shapes can be generated by changing control parameters of the class function.

B-Spline Shape Function.
The B-spline function can adjust the range of influence by selecting different orders of sub-Bernstein polynomials. It has better local control ability and computation efficiency than Bernstein polynomials. So it has been chosen as the shape function of the CST method in the basic cross section definition. The B-spline function is combined by the massive low-order Bezier curves to approximate the high-order curve. It is defined as follows: where , = 0, 1, . . . , is the weight factor. + 1 order Bspine basic function is defined as the piecewise polynomials at node vector =0 , which can be International Journal of Aerospace Engineering 3 obtained by recurrence formula De Boor-Cox.
The basic cross section defined with B-spline function is as follows: The three-dimensional geometry can be considered as a series of original cross sections arranged parallel along the axial direction. By defining proper cross sections using the method mentioned above, we can get an analytical surface to express geometric shape. In order to introduce an analytical axial rule, the coefficient is replaced by the same B-spline CST definition in formula (5). So the two-dimensional CST function can be expanded to a three-dimensional version, which is shown as follows: where = / . is the axial coordinate. is the normalized axial coordinate. is the total length of the generated surface in the axial direction. Bringing formula (6) to formula (5), we have the definition of the three-dimensional analytical surface.
( , ) upp,low Formula (8) is the analytical expression of the extensional three-dimensional surface, where is the third-dimensional coordinate value along the direction based on the twodimensional normalized coordinate and . ( ) and ( ) are the CST's class functions with control parameter 1, 2, 1 and 2. ( ) and ( ) are the basic functions of the B-spline, which constitute the CST's shape function. , is the discrete control weight factor in the up and low surfaces. Lower case and are the orders of B-spline function, which are expressed as the numbers of lateral and axial control points. The total points of the geometric surface are ( + 1) × ( + 1). More control points mean more design parameters and better parametric geometric shape design ability. When = 0 and = 0, the discrete control weight factor is equal to 1, the shape function ( , ) is also equal to 1, and then entire geometric surface is controlled by the control factors 1, 2, 1, and 2 of the class function. Δ , ( , ) is the eccentric distance in normal position, and the default value is 0.
The geometric surface defined in formula (8) can be described as the third-dimensional coordinate value calculated by the two-dimensional mesh points dispersed as the lateral and axial control points. and are the normalized coordinates defined in [0, 1] × [0, 1]. We also need to define the profile in the lateral direction with respect to the axial coordinate . The profile is expressed the same as the Bspline CST method in formula (5).
is the discrete control weight factor in the lateral coordinate. From the above definition, we obtain the threedimensional analytical surface, which can be transformed into global Cartesian coordinate as formula (10).
where = [0, 1], = [0, 1]. The outward boundary box is used to control the size of the geometric surface, and then we get the full design parameters of the entire geometric surface, which includes the following: Axial and lateral length: , Sectional section control factor: 1, 2 Vertical section control factor: 1, 2 Lateral section control factor: 1, 2 Surface weight factor: , Lateral weight factor: , , 1 , 2 are the size parameters of the geometric surface. 1, 2, 1, 2, 1, 2 are the control parameters of the class function. , and are the discrete control weight factors of the shape function. When , and are equal to 1, the number of total parameters is constant and also minimum. When the , and are expressed as matrix, some control points with the same size of the matrix elements are located in the geometric surface to improve the complexity of the geometric shape. So the size of design parameters could be adjusted dynamically. With these weight factors, the three-dimensional CST method has more flexible parametric geometric shape design ability.

Basic Geometric
Shape. The basic geometric shape is defined by the upper geometric surface and the inverted lower geometric surface. Figure 1 gives some geometric shapes generated by the above three-dimensional CST method with different control parameters. upp, low, and are the matrix generated by , of the upper surface, , of the lower surface, and of the lateral coordinate. The shapes in the same lines are generated with the same value of 1 and 2, but with different 1, 2, 1, and 2. The shapes in the same columns are generated with the same values 1, 2, 1, and 2, but with different 1, 2, and upp. The modified shapes in the third line are generated with different sizes of control weight factors of the upper surface compared with the shapes in the second line.
From the various geometric shapes and cross sections shown in Figure 1, the following features can be found for the three-dimensional CST method: (1) The geometric shape can be expressed as an extruded body made up of a series of basic cross sections with different control factors along the axial direction. The transformation rule of the body is decided by the control factors 1 head , 2 head at head and 1 tail , 2 tail at tail.
(2) The shape profile at the axial and lateral direction is decided by the control factors 1, 2, 1, and 2. It allows continuous control of every single basic cross section. When and are equal to 0 simultaneously, the corresponding cross section is opened. When and are greater than 0, the corresponding cross section is closed to a point. This feature can be used to generate the closed head and the opening body.
(3) upp, low, and are matrix with arbitrary size. They produce a global correction or a local adjustment to the geometric shape. The control ability is decided by the size of the matrix.
The shape profile is parameterized by the section control factors, and the geometric shape can be adjusted globally or locally by the weight factor matrix. These features give the three-dimensional CST method a comprehensive parametric geometric shape design ability with adjustable control parameters. Designers could have a more flexible design pattern to parameterize various geometric model.

Mesh Discretization.
A generic mesh discretization method is utilized to generate the corresponding mesh, which will be used as the aerodynamic mesh and structural mesh. The procedure for the three-dimensional CST method is likely to generate structured mesh surface along the and direction. With these characteristics, the surface is dispersed into some control points uniformly. Then a two-dimensional matrix made up of control points is generated, which can be expressed as follows: Node upp ( , ) = ( ( ) , ( , ) , ( )) , International Journal of Aerospace Engineering 5 where = 1, 2, . . . , − 1, = 1, 2, . . . , − 1. and are the numbers of total points in the and direction. The quadrangle is used as the mesh element consisted with these discrete control points Element ( ) = [Node ( , ) , Node ( , + 1) , The Node( , ) and Element( ) are used to generate the mesh of the analytical geometric shape. Then we get the corresponding quadrilateral mesh.

Three-Dimensional Parametric Geometry Modeling Method
With the above three-dimensional CST method, we get the continuous and smooth geometric shape. But single geometric shape cannot provide enough degree of freedom to generate an entire complex aircraft. Usually a complex aircraft can be split into some standard components. So with the component combination method, a universal three-dimensional parametric geometry modeling method is proposed based on the aircraft characteristic components library. The typical aircraft characteristic components include fuselage, wing, and empennage, which are introduced as follows.

Aircraft Characteristic Components Library
3.1.1. Fuselage. The fuselage can be divided into head, headbody, and head-body-tail, three types according to the structural cabins or aerodynamic shape features. By setting the control factors 1, 2, 1, and 2 to 1 or 0, it is easy to generate the semiclosed head and tail, and the opened body. The control factors 1 and 2 must be kept the same value at the joint position between two parts to ensure surface continuation. The detailed control parameters and the three fuselage types are shown in Figure 2(a The weight factors of the head include The formulas (13) and (14) are the generic forms of a single part of the component. The others are similar too. By varying these control parameters, various geometric shapes can be generated. The parametric geometric shape design ability is dependent on the size of the matrix and . For complex geometric shape, designers can set and to some arbitrary matrix to increase the degree of freedom. But one should note that the number of control parameters also increases as the matrix becomes larger.

Wing.
Wing modeling is similar to the fuselage. The control factors 1, 2, 1, and 2 are set to 0 to describe the opening tip and root of the wing. The sweepback angle and dihedral angle are also parameterized by mesh offset operation as follows: is the sweepback angle and is the dihedral angle. The single wing and double wing are modeled in the components library. The detailed control parameters and the two wing types are shown in Figure 2

3.1.4.
Others. Some typical components are introduced above, which can be used to create many kinds of common aircraft. There may be other requested components such as engine and nozzle. These can be some simplification of the basic three types components library. For example, the nozzle may be some simplification of the single part fuselage with opening head and tail. For other special components, designers can model it using the three-dimensional CST method and package to the aircraft characteristic components library for reuse. ] .

Three-Dimensional Geometry
The rotation matrix and transfer matrix are as follows: With these matrices, the modified control points new are evaluated in formula (18). Then the modified component mesh is generated by these points using the mesh discretization method discussed in Section 2.5.

Universal Structural Modeling Method
Although structural finite element analysis has been widely used in industrial department and research institution, the structure parameterized modeling is still a difficult phase. The most common structural parameterization method is based on parametric CAD model. It is complicated to model complex aircraft and also difficult to support the changing of the structural layout. These limitations cannot satisfy the needs of rapid modification and iterations in the conceptual design phase of the modern aircraft design. This section will provide a detailed illustration of a novel structural modeling method based on the three-dimensional parametric geometry modeling method mentioned above. Structural model is constructed based on the aircraft characteristic components library. Three typical structural models are elucidated below.

Fuselage.
Fuselage is the main body of an aircraft. It is usually designed as a thin-walled structure and constituted by longitudinal stiffeners (like beam and stringer), transverse stiffeners (like bulkhead), and surface skin. A fuselage contains fuel tank, payload, electronic instrument, and other equipment. Also it is as a sole central part connected with wings, empennages, engines, and other components to make up the entire aircraft. The main structural forms of fuselage include girder structure, longeron structure, and monocoque shell structure. Based on the above assumption, all of these structural forms can be simplified as three basic elements: beam/stringer, bulkhead, and skin. The detailed modeling method of these three elements is shown as follows.
4.1.1. Beam/Stringer. Beam and stringer are the longitudinal stiffeners of fuselage structure, used for undertaking axial load coming from fuselage bending, and also used for transferring the load of outer surface skin to bulkhead. Here, the beam and stringer are simplified to one-dimensional beam element. The cross section shape of the beam and the stringer, like I-section, H-section circle-section and so on, can be ignored at the modeling phase and will be considered by the structural solver.
Assuming the surface mesh of fuselage is continuous and compatible among all parts of the component, the discrete nodes of the upper surface or the lower surface can be expressed as where = 1, 2, . . . , , = 1, 2, . . . , . and are the numbers of control points along the and direction. Vector is used to locate the position of the beam and the stringer along transverse direction.
where 0 ≤ 1 < 2 < ⋅ ⋅ ⋅ < ≤ 1. is normalized to [0, 1]. represents the normalized lateral position of the th beam/stringer. For the discrete nodes, it can be handled as the nearest point Node Fuselage ( , ), which has the minimum distance to the position in the th line of the discrete surface consisted of Node Fuselage ( , ). So the node set of the beam/stringer can be defined as Node ( ) = Node ( , ) , = 1, 2, . . . , .

(23)
Then the mesh of beam/stringer can be extracted from the surface mesh of the fuselage. Figure 3 shows the beam/stringer structural layout scheme of a hypersonic aircraft.

Bulkhead.
Bulkhead can be divided into three types: normal bulkhead, reinforced bulkhead, and connected bulkhead. It supports the skin of the fuselage to maintain geometric shape and also undertakes some concentrated load. The definition of bulkhead is similar to the beam/stringer. One-dimensional beam element is selected to simplify the bulkhead model. The cross section shape of the beam also is ignored at the modeling phase and will be considered by the structural solver. Vector is used to locate the position of the bulkhead.   Then the mesh of the bulkhead can be extracted from the surface mesh of the fuselage. Figure 4 shows the bulkhead structural layout scheme of a hypersonic aircraft.

Fuselage Skin.
Fuselage skin is used to maintain the geometric shape of the fuselage. It should be continuous and smooth to support the aerodynamic solver and structural solver. The discrete surface mesh of fuselage is used for both the structural mesh and the aerodynamic mesh. It has the following advantages: (a) the structural model and the aerodynamic model use the same surface mesh coming from mesh discretization in Section 2.5; the data conversion between these models are no longer needed; (b) the mesh of the beam/stringer and the bulkhead are both extracted from the discrete surface mesh; they are coordinated with surface skin mesh naturally. Based on the above definition, the whole structural model of the fuselage can be generated as Figure 5 stringer, and rib are imputed to the wing skin. Based on these assumptions, the wing structure can be simplified to some two-dimensional shell elements. Three typical basic elements are used to construct the wing structure.

Spanwise Elements.
Spanwise elements include spar and stringer. Assuming these elements are arranged along the chordwise direction rigorously. For the discrete nodes Node( , ) of the wing, vector is used to locate the position of spanwise elements.
represents the normalized spanwise position of the th spanwise elements. For the discrete nodes, it can be handled as the nearest point Node Wing ( , ), which has the minimum distance to the position in the th line of the discrete surface consisting of Node Wing ( , ). So the node set of spanwise elements can be defined as where = 1, 2, . . . , . is the number of the spanwise points. is used to define the node distribution in the thickness direction.
where 0 ≤ 1 < 2 < ⋅ ⋅ ⋅ < ≤ 1. is normalized to [0, 1] too. The layouts of the upper surface points and the lower surface points are the same. So by arranging some middle points between the upper surface points and the lower surface points by some node distribution, the total nodes of the wing spanwise section can be generated as linespace ( ) means to arrange points in the regularity of distribution . The element set can be generated by connecting these nodes. Figure 6(a) shows the spanwise structural layout scheme of wing.

Chordwise Elements.
Rib is the chordwise element. Similar to the spanwise elements, vector is used to define the position of the rib, where 0 ≤ 1 < 2 < ⋅ ⋅ ⋅ < ≤ 1. is normalized to [0, 1]. The modeling process is the same as the spanwise elements; the node set can be generated as where = 1, 2, . . . , . is the number of the chordwise points.
is also used to define the node distribution in the thickness direction. So the total nodes of the wing chordwise section can be generated as (32) Figure 6(b) shows the chordwise structural layout scheme of wing.

Wing Skin.
Similar to the fuselage skin, the discrete mesh of the wing is used for both structural mesh and aerodynamic mesh to maintain consistency. Based on the above definition, the entire structural model of the wing can be generated as Figure 7

Others.
Most of the common components have the similar structural layouts to the fuselage and wing. They can be structured by the method mentioned above.

Aeroelastic Modeling and Analysis Process
The process of aeroelastic modeling and analysis includes three steps: three-dimensional CST modeling, aerodynamic/ structural modeling, and aeroelastic analysis. The process of aeroelastic modeling and analysis is shown in Figure 8, and the detailed illustration is as follows.

Three-Dimensional CST Modeling.
Section 3 lists all design parameters of the three common components. There are four types of design parameters available to use to control the geometric shape: layout parameter, size parameter, shape parameter, and local control parameter. The layout parameter and the size parameter are global design parameters. The former changes the position and the posture of the components, and the latter changes the shape size of the components. They can be used to control the global sizing of aircraft. The shape parameter and local control parameter are local design parameters. The former changes the sketch of the main cross section, and the latter adjusts the details of the component. They can be used to further control the detailed shape of the aircraft. Local control parameter can be any size of matrix. The larger the matrix, the better the parametric geometric shape design ability. According to the requirement of the geometry modeling complexity, the proper design parameters can be selected as the design variables; others may stay constant or vary with respect to the selected parameters. With these design variables, geometry model and corresponding mesh of the component can be generated automatically.

Aerodynamic/Structural Modeling.
A complex aircraft is constituted by several components. Applying the first step, the mesh of these components can be generated, respectively, to make up the entire aircraft. For each of these components, the structural mesh also can be obtained in accordance with the structural modeling method introduced in Section 4. But there are two problems that needed to be solved to generate the integrated aerodynamic/structural model.
The first problem is the redundant mesh caused by the mutual nesting between the connected components. This mesh could influence the precision of aeroelastic analysis. In the worst condition, it may cause some error to the aerodynamic solver and the structural solver. PINPOLYHEDRON, an open-source tool, is used to remove this nested mesh. PIN-POLYHEDRON is a C++ code. It provides a novel algorithm to test whether points are inside/outside/on a polyhedron defined by triangular faces and vertices. It can be used for various complicated models such as nonconvex volumes, multimaterial bodies, and so on, and there is no assumption about orientation of the face normal. Above all, the algorithm is very efficient especially for large-scale problems. In this research, the fuselage is as the main polyhedron. If a point of other components is detected inside the main fuselage, this point and the relevant mesh are both removed from the component. Looping all points of the component and removing the nested mesh, then the remaining mesh is taken as the valid mesh. This valid mesh is used as the aerodynamic mesh to evaluate the aerodynamic characteristics and is used as the structural mesh to analyze the structural performance.
The second problem is the connection of the relevant components. For aerodynamic analysis, the panel method is used as the aerodynamic solver, so the connection between components is not necessary. For structure analysis, the connection must be modeled to guarantee the force transmission between the connected components, assuming a structural connection exists only between the bulkhead of the fuselage and the beam/spar of other components. The following method is used to generate these structural connections.
(1) For each bulkhead of fuselage, evaluate the midpoint by counting the average coordinate of the points set on the bulkhead.
(2) For each beam/spar, find the nearest elements to the fuselage, mark as , evaluate the midpoint of , and mark as .
(3) For each beam/spar, find the nearest bulkhead by min | − |, and mark as * .
(4) Sort the points set * on the bulkhead * by counting the distance from , the sorted points set mark as * .
(5) The first points of the sorted points * are used as the connection points set 1 . The connection points set of the beam mark 2 and is the number of the connection points 2 .
(6) Generate some quadrilaterals by connecting the relevant points between the connection points set 1 and connection points set 2 .
(7) Looping all of the beams of the components, generate nodes set and elements set of the structural connection to obtain the connection mesh. Figure 9 shows the structural connection between the fuselage and wing generated by the above method. The red line is the bulkheads of the fuselage. The blue line is the beam/stringer of the fuselage. The quadrilaterals with green edges are the structural connection mesh. This method provides a simple way to generate the structural connection mesh between the main fuselage and the connected components.
This structural connection mesh guarantees the grid continuity for force transmission. However, the compelling connection may reduce the quality of these mesh and influence the connection stiffness between the fuselage and the wing. These influences can be corrected by modifying the material attribute or adding a spring element. In this paper, these impacts are ignored temporarily.

Aeroelastic Analysis.
The aerodynamic solver and the structural solver interplay in the aeroelastic analysis process. AeroCalc, a C++ in-house code with the panel method, is used as the aerodynamic solver. The modified Newtonian impact theory is used to evaluate the windward surface, and the Prandtl-Meyer theory is used to evaluate the leeward surface. NASTRAN is used as the structural solver. The structural mesh and control data are written in BDF format as the script file. Every part of the components has its own property section to define thickness, cross section parameters, and material attributes. The outside surface mesh of aerodynamic model and structural model is consistent, so the force obtained from aerodynamic solver can be applied to the structural model and the displacement obtained from structural solver can be applied to the aerodynamic model without any conversion. These simplify the aeroelastic analysis process tremendously.
The three-dimensional CST modeling method is used as a parametric geometry modeling and mesh generation processor to automatically generate aerodynamic and structural model for aeroelastic analysis. The designer can construct an aircraft geometry shape quickly and automatically obtain the corresponding aeroelastic model. An aeroelastic analysis can be carried out quickly and effectively to give a preliminary evaluation of the aeroelastic effect. This could be very helpful for aircraft preliminary design.
The aerodynamic model and structural model use the same surface mesh. The aerodynamic force can be applied to the structural model and the structural displacement can be added to aerodynamic mesh to generate new aerodynamic model directly. This simplification will limit the mesh adaptation for aerodynamic model and structural model, which may increase the computational cost and reduce analysis precision. But it could save the tedious aeroelastic meshing and modeling process, which is time-consuming, especially for complex aircraft. Also these automatic aeroelastic modeling and analysis process can be used for aeroelastic optimization.   head-body fuselage, a double wing, and a tail, are utilized to build up the RLV geometry. The main parameters of the RLV are listed in Table 1, and the trimmed aerodynamic model is shown in Figure 10(a). Eighteen bulkheads and eight beams are arranged in the head-body fuselage. Five spars and seven ribs are arranged in the double wing. Three spars and seven ribs are arranged in the vertical tail. Ten structural connections are used to connect the fuselage and the wing. Three structural connections are used to connect the fuselage and the tail. The integrated structural model is shown in Figure 10(b). The aluminum alloys are used as structural material for the structural frame and skin, with a density of 2,700 Kg/m 3 and an elastic modulus of 72 GPa. The total number of the aerodynamic elements is 15506, which is controlled by the number of points for mesh discretization in the and direction of each surface. The total number of the structural elements is 18691, including the skin mesh, which is the same as the aerodynamic mesh except the nested mesh, the structural component mesh, and the structural connection mesh.

Static Aeroelastic Analysis Example
The maximum dynamic pressure trajectory point during reentry is chosen as the aeroelastic evaluation point with = 3.6 and = 6.5 degree. The pressure coefficient distribution evaluated by the panel method at the maximum dynamic pressure point is shown in Figure 11(a). The corresponding aerodynamic force is applied to the structure surface directly without interpolation. The structural mesh and the displacement distribution are shown in Figure 11(b). The relative error of the maximum deformation in the direction is used as the convergence criteria during the aeroelastic iterations and it is set to 1 − 5.
The aeroelastic analysis of the RLV converges shortly after five iterations. The iteration history of the lift coefficient , the drag coefficient , and the maximum deformation in the direction are shown in Figure 12. In the final converged condition, the is 1.3 percent lower and the is 0.37 percent lower, which lead to a 0.94 percent decrease in the lift-to-drag ratio than the initial undeformed condition. The maximum structure deformation in the direction occurs at the wingtip position. The maximum deformation at the first structure analysis and the last structure analysis are 0.1287 m and 0.1271 m. The wingtip deformation is shown in Figure 13. The blue line shows the undeformed wing tip, the orange line shows the deformed wing tip after the first structure analysis, and the red line shows the deformed wing tip at the convergence. The change of the wingtip deformation is   not obvious during the aeroelastic iterations, which reaches a stable level after the second step. The three-dimensional CST parameterization method and the aeroelastic analysis process are written in MATLAB code. The geometry modeling costs less than one second, and the generation of aerodynamic mesh and structural mesh cost 5.5 s. The analysis time of aerodynamic analysis and structure analysis per iteration and the total analysis time are shown in Table 2. The whole aeroelastic modeling and analysis process is very fast and efficient. To generate an entire   complex RLV and its analysis models only needs several seconds. This will simplify and shorten the aeroelastic analysis and make it easy to use in the aircraft conceptual design phase.

Conclusion and Future Work
In this paper, a novel aerodynamic and structural modeling method based on the three-dimensional CST is developed to provide a fast and simple way to carry out an entire aeroelastic analysis process. The aerodynamic model is generated by trimming the embedded aircraft characteristic components. The structural model is generated by arranging proper inner parts from the basic mesh topology. The consistency between the aerodynamic model and the structural model is satisfied naturally. Finally, an efficient aeroelastic analysis process is created and tested by a RLV. In conclusion, (1) the three-dimensional parametric geometry modeling method gives a universal way to generate geometry model of common three-dimensional complex aircraft. It inherits the advantages of the original CST method and has fast and stable parametric geometric shape design ability. An aircraft can be modeled with a few control parameters in a few minutes and discretized to surface mesh. More common aircraft characteristic components library will be introduced in future work to give a broader support to various kinds of aircraft; (2) a universal aeroelastic modeling and analysis process is introduced. By using the three-dimensional parametric geometry modeling method, the aerodynamic model and the structural model can be generated together to keep the mesh consistency. Fluid structure interaction (FSI) can be ignored to simplify the aeroelastic analysis process. This simplification may limit the mesh adaptation for aerodynamic model and structural model, but it gives a fast and simple way to carry out an aeroelastic analysis process for complex aircraft, which is very useful in the conceptual design phase; (3) the structural modeling method in this article gives an automatic and fast structural layout parameterization and generation approach. Topology of the structure can be parameterized from a wide range stably. This could be used as parameterized structural model for topology optimization; (4) the generated mesh by the three-dimensional parametric geometry modeling method will be of poor quality at the position where the slope of the surface along axial direction and lateral direction is too large. These are caused by the uniform mesh discretization along the axial direction and lateral direction. These narrow grids may reduce the precision of the aeroelastic analysis. Nonuniform mesh discretization and mesh repairing method will be studied in future work to improve mesh quality.