The impact and sensitivity analysis of beam and stringer for stiffness center of composite wing structure

<!-- metadata commented in wiki content


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<big>'''The impact and sensitivity analysis of beam and stringer for stiffness center of composite wing structure'''</big></div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Shengjun Qiao<sup>1</sup><span id="fnc-1"></span><span style="text-align: center; font-size: 75%;">[[#fn-1|<sup>1</sup>]]</span>, Ruixue Ma<sup>1</sup>, Jin Jiao<sup>1</sup>, Xiangteng Ma<sup>2</sup>, Xing Liu<sup>1</sup></div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">1. School of Aircraft Engineering, Xi′ an Aeronautical University, Xi′an 710077, China.</span></div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">2. Department of Engineering Mechanics, Northwestern Polytechnical University, Xi′an 710129, China.</span></div>
-->

'''Abstract: '''Wing stiffness center should be determined firstly for structure detail design. The present study focused on the impact analysis of beam and stringer for wing cross-section stiffness center based on thin wall structure mechanics theory. In order to discuss the impact of beam and stringer on stiffness center, the sensitivity formulas for stiffness center of wing cross-section were derived and expressed in terms of the dimension and layout of beams and stringers. The results indicated that the structural layouts of beam and stringer were important influencing factors in stiffness detail design of full composite wing structure. The research results can provide an important reference for the stiffness design and aeroelastic design of the full composite wing.

'''Key words: '''Composite wing, Stiffness center, Layout, Dimension, Sensitivity

=1 Introduction=

In structure detail design of aircraft, the stiffness center of the wing structure is an important factor which must be considered. Not only the aerodynamic force was impacted by stiffness center of composite wing, but also the distribution of static and the deformation of aeroelastic [1]. When the cross sections from wing root to wing tip are selected, the wing stiffness center axis can be determined by calculating stiffness center of each cross-section.

At present, many researchers focus on studying the stiffness center of the wing section. They are based on structural mechanics theory and use model simplification methods to determine the stiffness center position [2-5], which is efficient and utility in engineering design. Based on the N-S equation, Ye et al. [6] not only analyzed the airfoil flow field characteristics, but also studied the influence of different positions of the elastic axis on the aerodynamic characteristics. Cesnik et al. [7] carried out the optimization design of the internal layout of the blade structure by controlling the position of the elastic axis of the helicopter blade, and in order to maximize the static twist actuation of the blade. The effect of the elastic axis was investigated, with the elastic axis closing to the airfoil midpoint, the pitch amplitude increases slightly while the plunge amplitude fluctuates with an extremum [8].

Sensitivity is also the derivative information, which reflects the impact of design variable on objective function or constraint function [9]. As an important research field in engineering, analysis was used previously in shape optimization of continuum structures and structure response [10-11]. Later according to different design requirements, sensitivity formulas were derived individually based on shape design variations or structure parameters [12-15], structure static responses, dynamic responses [16-21], and linear or nonlinear functions [22-25]

Before starting the calculation of the wing stiffness center, the paper needs to simplify the cross-sectional structure of the wing using an equivalent method.Firstly, based on the simplified cross-sectional structure of the wing, the calculation formula of the cross-sectional stiffness center is derived by using the related theory of thin-walled structural mechanics.Secondly, the impact analysis of dimension and layout for beam and stringer on the stiffness center of the wing was accomplished. Finally, sensitivity analysis was investigated on the influencing factors of beam and stringer dimension and layout for the wing stiffness center. The research results can provide certain guidance in the process of stiffness design and aeroelastic design of composite wing.

=2 Cross-section stiffness center=

Wing stiffness centers were calculated for different cross-sections selected from root to tip. However, in the process of stiffness centers calculation, in order to facilitate calculation and programming, it is necessary to simplify the area and materials of different structures in the cross section. The beam edges and stringers are respectively equivalent to a point on the section, and different materials are equivalent to the same material. The simplified model for stiffness center calculation is shown in Fig. 1.
[[File:Fig. 1 .jpg|centre|500x500px|
'''Fig. 1 Simplified wing cross-section model'''
|thumb]]

2.1 Skin area

In order to obtain the stiffness data and the stiffness center position required for the aeroelastic analysis, the cross-sectional area of the skin is taken as the effective area, and the definite integral is used in the calculation process. The formula for calculating the cross-sectional area A of the skin in the interval [a, b] can be written as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">{A}_{skin}=\int_{a}^{b}tds</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}


Where, ''A<sub>skin</sub> ''represents the cross-sectional area of the skin and ''t'' represents the thickness of the skin.

2.2 Stringer area

The stringer mainly bears the axial force caused by the bending moment of the wing. When establishing the stiffness center calculation model, the section is generally defined on the midplane of the skin, so that the cross-sectional area of the stringer should be converted. The simplified stringer is a rod with a concentrated area, and its entire area is concentrated on a point in the skin plane. The cross-sectional area of the stringer after conversion can be written as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">{A}_{stringer}={A}_{0}\frac{{{h}_{1}}^{2}}{{{h}_{0}}^{2}}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}


Where, ''A<sub>stringer</sub>'' and ''A<sub>0</sub>''are reduced and initial area, respectively. ''h<sub>1</sub> ''is the distance between skin and neutral surface of airfoil cross section, ''h<sub>0 </sub>''is the distance between stringer centroid and neutral surface of airfoil cross section.

2.3 stiffness center calculation

Cross-sectional area equivalent method based on sections 2.1 and 2.2 , in order to calculate the stiffness center for wing cross sections, first of all, we need to calculate the total area A of the components of the wing section under normal stress, and the centroid coordinates of the section denoted by <math display="inline">{\epsilon }_{0}</math> and <math display="inline">{\eta }_{0}</math>.

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>A={\oint {\varphi }_it_ids+\sum_i{\varphi }_iA_i}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\xi }_0=\frac{\oint {\varphi }_it_ix_ids+\sum_i{\varphi }_iA_ix_i}{A}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\eta }_0=\frac{\oint {\varphi }_it_iy_ids+\sum_i{\varphi }_iA_iy_i}{A}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}


Where  <math>{\varphi }_i</math> is reduction factor,  <math>t_i</math> is the thickness of skin which number is  <math>i</math> ,  <math>A</math> is the sum area of all component cross section,  <math>{\xi }_0</math> and  <math>{\eta }_0</math> are  <math>x</math> coordinate value and  <math>y</math> coordinate value of centroid, respectively.

According to coordinate transfer theory, in the area inertia coordinate y and z with individually representing wing chord and thickness orientation, the moment of inertia for every reduced cross sections on wing can be written as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| [[Image:Draft_qiao_302098793-image13.png|216px]]
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| [[Image:Draft_qiao_302098793-image14.png|216px]]
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
|}


The Eq (7) is the moment of the cross sections in the inertia axis, and then stiffness center of every cross section on wing can be calculated by Eq (8) according to the thin wall structural mechanics theory.

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| [[Image:Draft_qiao_302098793-image15.png|210px]]
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
|}


where  [[Image:Draft_qiao_302098793-image16.png|18px]] is twice the area of enclosed cross sections, and  [[Image:Draft_qiao_302098793-image17.png|12px]] is the distance from the centroid to the tangent of each component.

=3 The influence analysis of stiffness center on beam=

As the main frame, the wing beam plays an important role in wing structure. The wing stiffness center is impacted by beam dimensions and layouts. The influencing degree can be found eventually by comparing the different stiffness centers, which were gained by changing beam dimensions and layouts. Wing finite element model was shown in Fig. 2.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:Fig. 2 Wing finite element model.jpg|centre|500x500px|'''Fig. 2 Wing finite element model'''|thumb]]</div>

3.1 Beam dimension

Stiffness center from wing root to tip is different. In order to research the impact of the beams, three wing stiffness centers are analyzed using three different wing models, which are original dimensions, enlarged beam dimensions and reduced beam dimensions at the same time, as shown in Table 1. The wing is divided into three sections according to different beam dimensions and numbers as shown in Fig. 2.

<div id="_GoBack" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 1 Beam dimensions'''</div>

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;text-align: center;"|Model
|  style="border-top: 1pt solid black;text-align: center;"|Section A
|  style="border-top: 1pt solid black;text-align: center;"|Section B
|  style="border-top: 1pt solid black;text-align: center;"|Section C
|-
|  style="border-top: 1pt solid black;text-align: center;"|Original
| style="border-top: 1pt solid black;text-align: center;" |3
|  style="border-top: 1pt solid black;text-align: center;"|2
|  style="border-top: 1pt solid black;text-align: center;"|2
|-
|  style="text-align: center;"|Enlarged
|  style="text-align: center;"|3.5
|  style="text-align: center;"|2.5
|  style="text-align: center;"|2.5
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|Reduced
|  style="border-bottom: 1pt solid black;text-align: center;"|2.5
|  style="border-bottom: 1pt solid black;text-align: center;"|1.5
|  style="border-bottom: 1pt solid black;text-align: center;"|1.5
|}<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:Fig. 3 Wing cross sections.jpg|centre|500x500px|'''Fig. 3 Wing cross sections'''|thumb]][[File:Draft_qiao_302098793-image20.png|centre|450x450px|'''Fig. 4 Stiffness center of enlarged beam and original'''|thumb]]</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:Draft_qiao_302098793-image21.png|centre|thumb|450x450px|'''Fig. 5 Stiffness center of reduced beam and original model''']]</div>

Wing cross section is shown in Fig. 3. Stiffness center of enlarged beam dimension is compared with the original by linear fitting as shown in Fig. 4. Similarly, the reduced beam dimension is also compared with the original by linear fitting as shown in Fig. 5. Spacing of front beam and rear beam in the position of wing root is 560.606 mm, which is 197.08 mm in the position of wing tip. If  [[Image:Draft_qiao_302098793-image22.png|12px]] denotes the distance between stiffness center and front beam, and  [[Image:Draft_qiao_302098793-image23.png|12px]] denotes the distance between front beam and rear beam, then the value of C (C= [[Image:Draft_qiao_302098793-image22.png|12px]] / [[Image:Draft_qiao_302098793-image23.png|12px]] ) in the positions of wing root and wing tip under different models can be obtained, as shown in Table 2. It is indicated that the enlarged and reduced beam dimensions have small impact comparing with the original model.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 2 The values of C in the positions of wing root and wing tip for different models'''</div>

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Original model (%)
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Enlarged beam dimension (%)
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Reduced beam dimension (%)
|-
|  style="border-top: 1pt solid black;text-align: center;"|Wing root
|  style="border-top: 1pt solid black;text-align: center;"|44.7721
|  style="border-top: 1pt solid black;text-align: center;"|45.8022
|  style="border-top: 1pt solid black;text-align: center;"|46.3281
|-
|  style="border-bottom: 2pt solid black;text-align: center;"|Wing tip
|  style="border-bottom: 2pt solid black;text-align: center;"|43.4955
|  style="border-bottom: 2pt solid black;text-align: center;"|42.4467
|  style="border-bottom: 2pt solid black;text-align: center;"|43.2312
|}

3.2 Beam layout

This section focused on the impact analyze of wing cross sections stiffness center on beam by comparing stiffness centers of different models. The percentage curves of the stiffness center with respect to the change of front beam and rear beam are shown in Figs. 6 and 7 respectively.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:Draft_qiao_302098793-image24.png|centre|450x450px|'''Fig. 6 The percent curve of the stiffness center according to the changing of front beam'''|thumb]][[File:Draft_qiao_302098793-image25.png|centre|450x450px|'''Fig. 7 The percent curve of the stiffness center according to the changing of back beam'''|thumb]]</div>

The position of the stiffness center is impacted by the layout of the front beam and rear beam. The positions of them were defined by percent of the distance between them and the front border with respect to the total wing width. When changing the percent of the front beam positions from 15.343% to 21.98% of total wing width, the percent of the stiffness center positions varies from 40.176% to 42.128% correspondingly. Similarly, when changing the percent of the rear beam position from 65.288% to 71.926% of total wing width, the percent of the stiffness center positions varies from 40.061% to 41.976% correspondingly.

=4 The influence analysis of stiffness center on stringer=

4.1 Stringer dimension

The change in the cross-sectional dimension of the stringer will affect the cross-sectional structure of the wing and at the same time affect the stiffness centerof the cross-sectional wing. In order to research the impact of the stringers, three wing stiffness centers were analyzed using three different wing models, which are original dimensions, enlarged stringer dimensions and reduce stringer dimensions at the same time, as shown in Table 3.  Two Stringer cross sections are shown in Fig. 8.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 3 three stringer dimensions'''</div>

{| style="width: 100%;border-collapse: collapse;" 
|-
|  rowspan='2' style="border-top: 2pt solid black;text-align: center;"|models
|  colspan='4'  style="border-top: 2pt solid black;text-align: center;"|Top wing in section A
|  colspan='4'  style="border-top: 2pt solid black;text-align: center;"|Down wing in section A
|  colspan='4'  style="border-top: 2pt solid black;text-align: center;"|wing in section B and C
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|W
|  style="border-bottom: 1pt solid black;text-align: center;"|H
|  style="border-bottom: 1pt solid black;text-align: center;"|t
|  style="border-bottom: 1pt solid black;text-align: center;"|W1
|  style="border-bottom: 1pt solid black;text-align: center;"|W
|  style="border-bottom: 1pt solid black;text-align: center;"|H
|  style="border-bottom: 1pt solid black;text-align: center;"|t
|  style="border-bottom: 1pt solid black;text-align: center;"|W1
|  style="border-bottom: 1pt solid black;text-align: center;"|W
|  style="border-bottom: 1pt solid black;text-align: center;"|H
|  style="border-bottom: 1pt solid black;text-align: center;"|t1
|  style="border-bottom: 1pt solid black;text-align: center;"|t2
|-
|  style="border-top: 1pt solid black;text-align: center;"|original
|  style="border-top: 1pt solid black;text-align: center;"|40
|  style="border-top: 1pt solid black;text-align: center;"|30
|  style="border-top: 1pt solid black;text-align: center;"|3
|  style="border-top: 1pt solid black;text-align: center;"|25
|  style="border-top: 1pt solid black;text-align: center;"|30
|  style="border-top: 1pt solid black;text-align: center;"|30
|  style="border-top: 1pt solid black;text-align: center;"|3
|  style="border-top: 1pt solid black;text-align: center;"|25
|  style="border-top: 1pt solid black;text-align: center;"|15
|  style="border-top: 1pt solid black;text-align: center;"|25
|  style="border-top: 1pt solid black;text-align: center;"|2
|  style="border-top: 1pt solid black;text-align: center;"|2
|-
|  style="text-align: center;"|enlarge
|  style="text-align: center;"|45
|  style="text-align: center;"|35
|  style="text-align: center;"|3.5
|  style="text-align: center;"|30
|  style="text-align: center;"|25
|  style="text-align: center;"|25
|  style="text-align: center;"|3.5
|  style="text-align: center;"|30
|  style="text-align: center;"|18
|  style="text-align: center;"|30
|  style="text-align: center;"|2.5
|  style="text-align: center;"|2.5
|-
|  style="border-bottom: 2pt solid black;text-align: center;"|reduce
|  style="border-bottom: 2pt solid black;text-align: center;"|35
|  style="border-bottom: 2pt solid black;text-align: center;"|25
|  style="border-bottom: 2pt solid black;text-align: center;"|2.5
|  style="border-bottom: 2pt solid black;text-align: center;"|20
|  style="border-bottom: 2pt solid black;text-align: center;"|20
|  style="border-bottom: 2pt solid black;text-align: center;"|20
|  style="border-bottom: 2pt solid black;text-align: center;"|2.5
|  style="border-bottom: 2pt solid black;text-align: center;"|20
|  style="border-bottom: 2pt solid black;text-align: center;"|12
|  style="border-bottom: 2pt solid black;text-align: center;"|20
|  style="border-bottom: 2pt solid black;text-align: center;"|1.5
|  style="border-bottom: 2pt solid black;text-align: center;"|1.5
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
{|
|-
|[[File:Draft_qiao_302098793-image26.png|centre|300x300px|(a)|thumb]]
| [[Image:Draft_qiao_302098793-image27.png|center|196x196px|(b)|thumb]]
|}
'''Fig. 8 Stringer cross-section    (a) Section A  (b) Section B and C'''
</div>

These stiffness centers of the two wing models with enlarged stringer dimensions and original are compared, as shown in Fig. 9. Similarly, these stiffness centers of the two wing models with reduced stringer dimensions and original model are also compared in Fig. 10.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:Draft_qiao_302098793-image28.png|centre|450x450px|'''Fig. 9 Stiffness center of enlarged stringer and original model'''|thumb]][[File:Draft_qiao_302098793-image29.png|centre|450x450px|'''Fig. 10 Stiffness center of reduced stringer and original model'''|thumb]]</div>

The values of C in the positions of wing root and wing tip for different models can be obtained, as shown in Table 4. It is indicated that enlarged and reduced stringer dimensions have little impact comparing with the original model.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 4 The value of C in the positions of wing root and wing tip under different models'''</div>

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Enlarged stringer dimension (%)
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Reduced stringer dimension (%)
|-
|  style="border-top: 1pt solid black;text-align: center;"|Wing root
|  style="border-top: 1pt solid black;text-align: center;"|45.779
|  style="border-top: 1pt solid black;text-align: center;"|46.553
|-
|  style="border-bottom: 2pt solid black;text-align: center;"|Wing tip
|  style="border-bottom: 2pt solid black;text-align: center;"|44.1105
|  style="border-bottom: 2pt solid black;text-align: center;"|43.3256
|}


4.2 Stringer layout

Wing stiffness centers were different for the different stringer layouts. In order to study the influence of the structural layout of the stringer on the stiffness center, it is necessary to keep the layout and dimensions of other cross-sectional structures unchanged to achieve the purpose of the research. Wing cross section is shown in Fig. 3, different stiffness centers are gained when changing the middle group stringer position in the four models, which are shown in Fig. 11.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:Draft_qiao_302098793-image30.png|centre|450x450px|'''Fig. 11 The percent curve of the stiffness center according to middle group stringer on different stringer spacing'''|thumb]]</div>

When changing the percentage of middle stringer positions of total wing width are from 37.232% to 49.828%, the percentage of the variation range of stiffness center positions in total wing width are individually shown in Table 5. It is indicated that stiffness center position has a great change for different models with different spacings of stringers, but has small impact on the spacing of stringers.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 5 the percents of the variation range of stiffness center positions for four models'''</div>

{| style="width: 100%;border-collapse: collapse;" 
|-
|  rowspan='2' style="border-top: 2pt solid black;text-align: center;"|Stringer spacing
|  colspan='2'  style="border-top: 2pt solid black;text-align: center;"|Variation range of stiffness center (%)
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|beginning
|  style="border-bottom: 1pt solid black;text-align: center;"|end
|-
|  style="border-top: 1pt solid black;text-align: center;"|104 mm
|  style="border-top: 1pt solid black;text-align: center;"|38.199
|  style="border-top: 1pt solid black;text-align: center;"|43.7
|-
|  style="text-align: center;"|116 mm
|  style="text-align: center;"|38.248
|  style="text-align: center;"|43.667
|-
|  style="text-align: center;"|128 mm
|  style="text-align: center;"|38.283
|  style="text-align: center;"|43.626
|-
|  style="border-bottom: 2pt solid black;text-align: center;"|138 mm
|  style="border-bottom: 2pt solid black;text-align: center;"|38.305
|  style="border-bottom: 2pt solid black;text-align: center;"|43.567
|}

=5 Sensitivity analysis for beam and stringer=

Sensitivity is the derivative information, which reflects the impact of design variable on objective function or constraint function. The sensitivity formula for wing cross sections stiffness is derived according to the thickness  [[Image:Draft_qiao_302098793-image31.png|12px]] of beam and stringer, which is written as
[[File:(9).png|centre|thumb|629x629px|(9)]]
with

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:Draft_qiao_302098793-image33.png|centre|thumb|248x248px]]</div>

As the distance  [[Image:Draft_qiao_302098793-image34.png|12px]] between wing beam or middle group stringer and inertia axis of wing cross sections, the sensitivity formula for wing cross sections stiffness is derived as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
|[[File:Draft_qiao_302098793-image35.png|centre|thumb|507x507px]]
|}
(10)
|}


with

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:10-2.png|centre|thumb]]</div>

Mean value sensitivity reflects the influencing degree of the average value of variable for structure. According to the formulation (4), when enlarging or reducing the thicknesses of beam and stringer at the same time, the results of mean sensitivity for the stiffness center of wing cross sections for the thickness of beam and stringer are obtained in Table 6.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 6 Mean value sensitivity for thickness of beams and stringers'''</div>

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Variables
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Beams
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Stringers
|-
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|Mean value sensitivity
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|1.691
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.0352
|}

When researching the impact of the front beam and rear beam layouts for stiffness center, it is necessary to keep the positions of other structure while only change the position of front beam or rear beam. According to the Eq (5), the results of mean sensitivity of beam positions of stiffness center are obtained in Table 7.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 7 Mean value sensitivity for layout of beams'''</div>

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Variables
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Front beam
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Rear beam
|-
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|Mean value sensitivity
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.2906
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.288
|}

The distribution of the stringers for wing cross section was shown in Fig. 2. In order to research the impact of different stringer spacing for stiffness center, 4 models with different stringer spacing were obtained. Change the position of stringers at the same time when middle group stringer is the design variable of sensitivity, the results of mean value sensitivity for stringer positions about stiffness center for different models are obtained in Table 8.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 8 Mean sensitivity for layout of stringers on different models'''</div>

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;text-align: center;"|Variables
|  colspan='4'  style="border-top: 2pt solid black;text-align: center;"|Middle group stringer
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|Models
|  style="border-bottom: 1pt solid black;text-align: center;"|1
|  style="border-bottom: 1pt solid black;text-align: center;"|2
|  style="border-bottom: 1pt solid black;text-align: center;"|3
|  style="border-bottom: 1pt solid black;text-align: center;"|4
|-
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|Mean value sensitivity
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.4462
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.4413
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.4328
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.4174
|}

We concluded that it has small impact by changing the beam and stringer thickness for position of stiffness center on full composite wing structure. However, it is large to the changing of the positions of front beam and rear beam. Comparing 4 different models with only different stringer spacing, it is found that the positions of wing cross sections stiffness center are all impacted largely for the 4 models.

=6 Conclusion=

This paper studies the influence of the size and layout of the stringers and beams on the stiffness center of the wing cross section. The sensitivity formulas for stiffness center of wing cross sections were derived and expressed in terms of variations with dimension and position of beam and stringer. The important conclusions obtained through the research are as follows:

* According to the impact analysis of beam and stringer for wing cross sections stiffness center, it is found that the stiffness centers of full composite wing are remarkable sensitive to structure layout of beam and stringer, which has small impact for the change of stringers spacing.

* Based on the formula of stiffness center, sensitivity formulas for stiffness center were derived, which was used accurately in analyzing the impact of beam and stringer. This effect of the structure layout of beam and stringer are considered as an important factor in stiffness design of full composite wing structure.

* The influence analysis and sensitivity of stiffness are very useful in the design of full composite wing structure.

==Acknowledgments==

This work is funded by the Scientific Research Foundation for Doctors, the Natural Science Basic Research Plan in Shaanxi Province of China (No: 2019JQ-912) and the Natural Science Foundation of Xi′ an Aeronautical University (No: 2018KY1226).

==Reference==

[1] D Guan, K Zhong. Aeroelastic tailoring applied to composite wing [J]. Acta Aeronautica Et Astronautica Sinica, 1989, 10(5): 221-226.

[2] S H Huo, F S Wang, Z Yuan, Z F Yue. Composite wing elastic axis for aeroelasticity optimization design [J]. Aircraft Engineering and Aerospace Technology, 2013, 85(1): 10-15.

[3] R J Guyan. Reduction of stiffness and mass matrices [J]. AIAA Journal, 1965, 13(2): 380-385.

[4] J Shawn. Advanced flexible wing technology assessment for transport applications [D]. Massachusetts Institute of Technology, U, S, 1995

[5] Y C LIU, G N HUANG, L J WANG. Computation of stiffness characteristic of a high aspect ratio wing [J]. Airplane engineering, 2006(4): 21-24.

[6] Z Y Ye, Y J Xie, J Wu. The effects of wind-tunnel model vibration on flow field and aerodynamics of an airfoil [J]. Engineering Mechanics, 2009, 26(4): 240-245.

[7] C E S Cesnik, J Mok, A S Parikh. Optimization design framework for integrally twisted helicopter blades [J]. AIAA 2004-1761, 2004(4): 1-15.

[8] D C Li, J W Xiang. Nonlinear aeroelastic analysis of airfoil using quasi-analytical approach [J]. Acta Aeronautica Et Astronautica Sinica, 2007, 28(5): 1080-1084.

[9] R T Haftka, H M Addlman. Recent developments in structural sensitivity analysis [J]. Structural Optimization, 1989, 3(1): 137-152.

[10] S Y Wang, Y Sun, R H Gallagher. Sensitivity analysis in shape optimization of continuum structures [J]. Computers and Structures, 1985, 20:855-867.

[11] P C Pandey, P Bakshi. Analytical response sensitivity using hybrid finite elements [J]. Computers and Structures, 1999, 70: 525-534.

[12] A Tafreshi. Shape design sensitivity analysis with respect to the positioning of features in composite structures using the boundary element method [J]. Engineering Analysis with Boundary Elements, 2006, 30: 1-13.

[13] J H Choi. Shape design sensitivity analysis for stability of elastic continuum structures [J]. International Journal of Solids and Structures, 2007, 44: 1593-1607.

[14] J H Choi. Shape design sensitivity analysis and optimization of general plane arch structures [J]. Finite Elements in Analysis and Design, 2002, 39: 119-136.

[15] C C Antonio, L N Hoffbauer. Uncertainty analysis based on sensitivity applied to angle-ply composite structures [J]. Reliability Engineering and System Safety, 2007, 92: 1353-1362.

<span id='OLE_LINK1'></span>[[#cite-OLE_LINK1|[16]]] T D Hien, M Kleiber. Stochastic structural design sensitivity of static response [J]. Computers and Structures, 1991, 38(5/6): 659-667.

[17] O Sergeyey, Z Mroz. Sensitivity analysis and optimal design of 3D frame structures for stress and frequency constraints [J]. Computers and Structures, 2000. 76: 167-185.

[18] Z D Wang, S H Chen. Sensitivity analysis of eigenmodes and dynamic responses for intelligent structures [J]. Finite Elements in Analysis and Design, 1999, 33: 71-81.

[19] Y Zhao, F Y Wang. Sensitivity analysis for structures subjected to stationary random excitations [J]. Procedia Engineering, 2011, 14: 2475-2480.

[20] M Ohsaki. Design sensitivity analysis and optimization for nonlinear buckling of finite-dimensional elastic conservative structures [J]. Comput. Methods Appl. Mech. Engrg., 2005, 194: 3331-3358.

[21] E P Petrov. Analysis of sensitivity and robustness of forced response for nonlinear dynamic structures [J]. Mechanical Systems and Signal Processing, 2009, 23: 68-86.

[22] S Etienne, D Pelletier. A general approach to sensitivity analysis of fluid-structure interactions [J]. Journal of Fluids and Structures, 2005, 21: 169-186.

[23] M G Kim, H L Jang, S Cho. Adjoint design sensitivity analysis of reduced atomic systems using generalized langevin equation for lattice structures [J]. Journal of computational physics, 2013, 240: 1-19.

[24] J M Pajot, K Maute. Analytical sensitivity analysis of geometrically nonlinear structures based on the co-rotational finite element method [J]. Finite Elements in Analysis and Design, 2006, 42: 900-913.

[25] H A Mang, G Hofinger, X Jia. On the interdependency of primary and initial secondary equilibrium paths in sensitivity analysis of elastic structures [J]. Comput. Methods Appl. Mech. Engrg., 2011, 200: 1558-1567.

----

Corresponding author.

''E-mail address:'' [mailto:qiaoshengjun123@163.com qiaoshengjun123@163.com]