Gust Alleviation of a Large Aircraft with a Passive Twist Wingtip

: This paper presents an investigation into the gust response and wing structure load alleviation of a 200-seater aircraft by employing a passive twist wingtip (PTWT). The research was divided into three stages. The first stage was the design and analysis of the baseline aircraft, including aerodynamic analysis, structural design using the finite element (FE) method and flutter analysis to meet the design requirements. Dynamic response analysis of the aircraft to discrete (one-cosin) gust was also performed in a range of gust radiances specified in the airworthiness standards. In the second stage, a PTWT of a length of 1.13 m was designed with the key parameters determined based on design constraints and, in particular, the aeroelastic stability and gust response. Subsequent gust response analysis was performed to evaluate the effectiveness of the PTWT for gust alleviation. The results show that the PTWT produced a significant reduction of gust-induced wingtip deflection by 21% and the bending moment at the wing root by 14% in the most critical flight case. In the third stage, effort was made to study the interaction and influence of the PTWT on the symmetric and unsymmetrical manoeuvring of the aircraft when ailerons were in operation. The results show the that PTWT influence with a reduction of the aircraft


Introduction
Aircraft will be inevitably exposed to atmospheric disturbances or gust of different degrees of severity during flight.The gust produces extra aerodynamic forces on the structure, especially the wing, and leads to the dynamic responses of the aircraft with undesired transverse movement.For some aircraft, the gust load is greater than the manoeuvring load and produces the most critical load case in the structure design.For civil airliners, especially with a high aspect ratio wing, the gust response causes concerns for the airframe safety, riding comfort and handling quality.To ensure aircraft safety, airworthiness authorities have specified a typical gust model and parameters as part of the certification requirements for structure design.Modern aircraft structures have been significantly improved by employing composite materials and advanced alloys.Nevertheless, gust load will cause an increase of the airframe weight, reduction of the structural fatigue life and higher maintenance cost.It is therefore desirable to alleviate the gust load by means of an active or passive control system in the first place.
In fact, various active and passive control technologies have been developed for flutter suppression and gust alleviation, and some been applied to aircraft [1][2][3][4][5].Comparing with the active control system, a passive control device is usually simpler and more reliable, although, maybe, less effective.This current investigation is focused on the effectiveness of gust alleviation by a specially-designed passive twist wingtip (PTWT).This particular PTWT has been studied in previous research for large unmanned air vehicles, including a flying wing aircraft [6,7] and jointed-wing sensor craft [8][9][10][11][12] of very flexible and high aspect ratio wings.The results have shown that the gust load has been significantly reduced by employing such a passive gust load alleviation device with optimized design parameters.This current investigation is aimed at evaluating the feasibility and effectiveness of the PTWT for gust response alleviation of a large civil aircraft of a high aspect ratio wing and conventional configuration.The example aircraft taken in the case study is a 200-seater airliner designed in an MSc group design project at Cranfield University in 2013 [20].The technical details are presented in Section 3. The study was conducted in three stages.
In the first stage, to ensure that the analysis and results are representable in general cases, the design and modelling of the example aircraft was carried out to meet the airworthiness requirements.Starting with the aerodynamic loading analysis, the composite wing structure was modelled and analysed by using the finite element (FE) method.Subsequent modal and flutter analysis was carried out to satisfy the aeroelastic design constraint.For this baseline aircraft model, gust response analysis to a discrete (one-cosin) gust was performed in a range of gust velocities and gradients specified in the airworthiness European Aviation Safety Agency (EASA) CS-25 [13].
In the second stage, an FE model of the passive twist wingtip (PTWT) of 1.13 m in length was created to replace the baseline wingtip.The key design parameters of the PTWT were determined for maximum gust alleviation subject to the required aeroelastic stability and normal operation at cruising.Gust response analysis of the wing and the whole aircraft with free heave motion was subsequently performed to evaluate the effectiveness of the PTWT in the specified gust range.The analysis results have shown that a significant reduction of gust response by employing the PTWT has been achieved in terms of the wing deflection and bending moment.
In the third stage, the investigation was extended to the influence of the PTWT on the control effectiveness of the aircraft in symmetric and unsymmetrical manoeuvrability.Although additional aerodynamic force will be produced mainly in the wing sections where ailerons were operated, the PTWT will react to the operation to some extent and have a potentially negative effect.In the case that the negative effect is unacceptable, a PTWT locking system would be necessary.The investigation has shown that the PTWT is a reliable and very effective approach to alleviate gust load and, therefore, has the potential for application.

Design Requirements and Procedure
Based on the conceptual design of an aircraft, the most critical load conditions are identified for the preliminary design of structures.For a large civil aircraft, the design has to comply with the CS-25 airworthiness requirements.In this investigation, the buckling analysis method [14] under practical design constraints [15] has been used in the stiffened skin panel design.For a wing made of carbon/epoxy composite materials, the usual laminate strength criteria in terms of the failure index (FI) are used for strength evaluation.One of the three widely used interactive failure criteria, the Tsai-Wu criterion, was also employed, since the other two would give similar results with negligible differences.In addition, the maximum strain of the structural laminate components is limited to ε = 0.0035 (3500 με) under the damage tolerance constraint.This figure was used as a general guideline for composite structure design in a world-leading civil aircraft industry.
Based on the preliminary design, the wing structure was modelled and analysed using the FE method, implemented in the Patran/Nastran package.After the static analysis to ensure the structural stiffness and strength requirements are satisfied, the modal and flutter analysis and dynamic response of the wing to 1-cosin discrete gust was conducted in a range of gust lengths specified in CS-25.The original 200-seater airliner having a conventional wing without the PTWT is taken as the baseline model.The PTWT FE model is then created to replace the baseline wingtip.Firstly, the PTWT effectiveness on gust alleviation was evaluated only for the wing clamped at the root.Subsequently, the gust alleviation was analysed for the whole aircraft with the fuselage and tail modelled as rigid bodies.Finally, concerning the potential influence of the PTWT on the flight performance, analysis was conducted on the symmetrical and unsymmetrical manoeuvring of the aircraft.

Aeroelastic Analysis
The aeroelastic analysis for large aircraft is normally performed in a generalized coordinate system in the frequency domain.The aeroelastic governing equations in matrix form can be expressed as below [16].
where Mhh, Dhh and Khh are the generalized mass, damping and stiffness matrix, respectively; M is the Mach number; k is the reduced frequency; {qh} is the generalized coordinate vector; and are the real and imaginary parts of the generalized unsteady aerodynamic matrix = [Ф] T {Lhh}[Ф], where {Lhh} is the unsteady aerodynamic force vector and [Ф] the modal matrix of the structure in spatial coordinates; { } = [Ф] T {LG} is the generalized aerodynamic force due to gust, where {LG} is the gust-induced aerodynamic force vector based on the gust model with details presented in the section below.
The unsteady aerodynamic force is calculated by using the doublet-lattice method [17,18] in Nastran.In the aerodynamic modelling, the airfoil profile effect is ignored, since it has been considered in the steady aerodynamic calculation.The wing surface is divided into planar trapezoidal segments, each having a constant dihedral with the sides parallel to the flow direction.These segments are further divided into panels with aspect ratios between 1/3 and 3.0.Higher density mesh is taken near the wing edges and those regions where discontinuous downwash and pressures of large gradients occur.The structural FE stress and modal analysis is performed by using Nastran.When both aerodynamic grids and structural meshes are created, a surface spline method was used to interpolate the aerodynamic grid points to structural mesh points for aeroelastic analysis.
For flutter analysis, the {QG} term in Equation ( 1) is zero, and the equation is then solved to determine the flutter speed by using the "pk" method.For gust response analysis with non-zero {QG}, the {qh} is solved and transformed back to the displacement in space coordinates.Both calculations use the aeroelastic module built in to Nastran.

Gust Model and Load
The air turbulence in real life encountered in flight is usually random and continuous in a short period.For large aircraft design purposes, two types of air turbulence models are accepted and specified in the airworthiness regulation CS-25.One is a continuous model of random characteristics represented by the power spectrum density and another simplified discrete gust model, such as the 1-cosine model.The external gust aerodynamic force and response depends on the input gust model.In this study, the main concern is the relative gust response of an aircraft with the PTWT compared with the baseline rather than the absolute value of the gust response.The 1-cosin discrete gust model is thus adapted in this current investigation.The gust velocity profile is expressed below [13]: where Uds is a specified design gust velocity, s is the distance of the aircraft penetrated in the gust and H is the gust gradient distance, which is the distance parallel to the airplane's flight path for the gust to reach its peak velocity.According to CS-25, the gust distance H is in the range from 9 to 107 m, including a typical value 12.5 ̅ , where ̅ is the wing mean chord.For an aircraft flying at speed V and encountering a discrete gust modelled by Equation ( 2), the gust penetrating distance is s = Vt, and Equation ( 2) can be rewritten as 0.5Uds (1 − cos(πVt/H)) or 0.5Uds(1 − cos(ωt)) in the time-domain.From this expression, an equivalent gust frequency for the dynamic response analysis can be calculated as ω = πV/H in radians or ω = V/2H in Hz.The design gust velocity (Uds), as defined in CS-25, can be calculated as follows: where the reference gust velocity Uref or gust intensity depends on the equivalent airspeed (EAS) and altitude and is specified in CS-25.At speeds between stall speed and cruise speed VC, Uref = 17.03 m/s at sea level; linearly decreasing to 13.41 m/s from sea level to 15,000 ft; further linearly decreasing to 6.36 m/s from 15,000 to 60,000 ft.At dive speed VD, Uref takes half of the specified value for VC.
The Fg is defined as the flight profile alleviation factor and expressed below.It increases linearly from sea level to 1.0 at the maximum operating altitude Zmo as follows [13]: According to the gust velocity, the aerodynamic force due to the gust can be calculated from: where ρ is the air density, s is the reference area of the wing, a1 is the lift curve slope of the wing and V0 is the equivalent flight speed.For a symmetric vertical gust, the gust alleviating factor is defined as: where = ̅ is the mass parameter, m is the mass of the aircraft and ̅ is the wing mean chord.
Substituting the following expression of aerodynamic force into Equation ( 5), Equation ( 5) can be rewritten in proportion to the steady aerodynamic force as: The gust-induced aerodynamic force LG was transformed to the generalized form {QG}, as defined earlier, and used to calculate the structural response in terms of deflection based on Equation (1).

Aircraft Data and Aerodynamic Load
The 200-seater airliner as illustrated in Figure 1a is taken as the baseline airframe.This aircraft has an H-tail, and twin open-rotor engines are mounted above the rear fuselage.This configuration makes the high aspect ratio wing clean for high aerodynamic efficiency, but also sensitive to gust load.The cruise speed (EAS) is Mach 0.78 or 231 m/s at an altitude of 35,000 ft.Some additional technical data of the aircraft design are listed in Table 1.The aerodynamic force distributions over the wing in trimmed level flight and under gust load without and with the PTWT in action are illustrated in Figure 1b.The wing airfoil chosen for the aircraft is a supercritical profile to keep natural laminar flow in cruise.To identify the most critical loading case for wing structure design, the initial analysis of aerodynamic loads in the combination of various manoeuvre cases was carried out based on lifting line theory and elliptical approximation of the spanwise lift distribution [19].For validation purposes, numerical simulation of the aerodynamic loads was subsequently performed by using the Fluent Package based on computational fluid dynamics (CFD) method.In the steady aerodynamic analysis, the aeroelastic effect was ignored.Figure 2 shows the aerodynamic pressure on a 2D wing section and the whole 3D wing at an angle of attack (AoA) = 5.54°.The difference of the resulting wingtip deflections due to the aerodynamic forces by the lifting line and CFD methods is within 5%, with the CFD results slightly greater [20].In this result, only the aerodynamic washout due to the spanwise jig angle was counted in the 3D wing shape.The aeroelastic washout due to structural deformation was not considered at this stage, since the influence would not be significant due to the small sweep angle of the wing.
The most critical load on the wing for structure design was calculated using the lifting line method by combining the aerodynamic forces for each section in the worst, but possible, conditions.It is usually represented in the form of shear force and the associated bending moment plot along the wingspan, as shown in Figure 3.The torque distribution along the wing span is similar to the bending moment and results in a maximum value of 1,000 kNm at the wing root.To simplify the aerodynamic analysis, the winglet was ignored in the current study.This led to the wingtip aerodynamic performance and, consequently, the effectiveness of PTWT gust alleviation being underestimated.

Structural Layout and Initial Design
The conventional structure layout of the fuselage and wing of the aircraft design has been adapted as illustrated in Figure 4a.The chordwise dimension of the flaps and aileron is 30% chord.The spanwise position and dimension of the aileron and elevator are also shown in Figure 4a.The primary structure of the composite wing is made of the upper and lower stiffened skin covers and 21 ribs supported by a front spar and rear spar located at 15% and 65% chord, respectively.Top-hat stringers and integral blade stringers with a spacing of 150 mm were selected to reinforce the upper and lower skin covers, respectively, as illustrated in Figure 4b.An intermediate modulus carbon/epoxy composite prepreg has been chosen with material modulus E1 and strength Xt and Xc in the fibre direction, E2, Yt and Yc in the off-fibre direction and shear modulus G12, strength S and Poisson's ratio v12, presented in Table 2.  Based on the loading and structural layout, an initial structural design with quasi-isotropic layup (50/40/10) (i.e., 50%, 40% and 10% of 0°, ±45° and 90° layers) for the skin and different layup for spars and stiffeners [17] was carried out.Buckling analysis of the stiffened composite skin panels was carried out by employing the finite strip method in Engineering Science Data Unit (ESDU) 03001 [14], together with the practical design rules given by Niu [15].As a result, the upper skin panel thickness varies from 22 mm inboard to 3 mm at the wingtip, as shown in Figure 5a.The buckling load factors are between 1.05 and 1.5.The lower skin thickness varies from 17 mm inboard to 3 mm at the wingtip with a failure index below 0.8 under an ultimate load factor of 3.75 along the wing span.Taking the outer wing box structure as an example, the wing component mass breakdown of a total mass of 1436 kg is shown in Figure 5b.

Structural FE Model and Stress Analysis
Based on the initial design, the FE model of the composite wing structure was created, using Patran/Nastran.In the FE model, the skins, ribs and spar webs are modelled using 18,500 shell elements for the outer wing; the stringers and spar capes are modelled by using beam elements.The FE analysis results show a maximum displacement of 1.78 m at the wingtip under the static ultimate load, as shown in Figure 3.The associated failure index (FI) of the upper skin panel is shown in Figure 6a, and the strain distribution of the upper and lower skin covers is shown in Figure 6b.
The results show that the maximum FI is less than 1.0 and the strain less than 3500 µε (micro-strain) in all regions of the wing box.It is observed that the most critical FI and strain close to the limit occur only at the wing root, especially around the rear spar, where high stress concentration exists due to the clamped boundary condition at discrete nodes.The results indicate that the design of the wing structure has met the strength and buckling requirements.

Aeroelastic and Gust Response Analysis
Based on the FE model, modal analysis was carried out as part of the aeroelastic analysis procedure.The first five frequencies in both maximum take-off mass (MTOM) and operational empty fuel mass (OEM) are listed in Table 3. Due to the inertia effect of the fuel in the wing box, the frequencies in MTOM are smaller than that in the OEM case.In the subsequent aeroelastic analysis, the five dominant bending and torsion modes were considered in this preliminary design stage.The results are presented in the V-g plot, shown in Figure 7a, indicating that the flutter speed is Vf = 580 m/s for MTOM and Vf = 820 m/s for the OEM case, which keeps a large safety margin from the limit of 283 m/s according to the CS-25 regulation.In the gust response analysis, four different values of gust gradient distance were selected from the lowest 9 m to the highest 107 m, including the typical value 12.5c = 48.5 m, according to the CS-25 specified range.The standard structural damping coefficient 0.03 was taken into account [22] in the model as expressed in Equation (1).Corresponding to the four gust gradient distances, the flight profile alleviation factor Fg = 0.974 is obtained from Equation (4), and the design gust velocity and gust equivalent frequencies in the most critical case at sea level and cruise speed were calculated and listed in Table 4.It is noted from Tables 3 and 4 that the gust load falls into the frequency range of the first two bending modes of the wing structure in the MTOM case and the first mode in the OEM case.The gust in the specified range is likely to cause a large gust response at the structure's natural frequency.For example, the gust gradient H = 25 m and frequency 3.5 Hz are identical to the first bending mode of the wing in the MTOM case, which could be critical.Taking the MTOM case at sea level, the gust response in terms of the wingtip deflection has been calculated for all four gust gradients and frequencies listed in Table 4. From the results, as shown in Figure 7b, it is noted, however, that the maximum wingtip deflection of 1.59 m actually corresponds to the gust gradient H = 107 m associated with the greatest gust velocity.The gust response is in addition to the static deflection in the 1g level flight.Comparing with the 1.78 m wingtip deflection and 3500 µε (micro-strain) under the ultimate load, this result indicates that the gust load reaches nearly 90% of the most critical ultimate load.In addition to the lifting force in trimmed level flight, the maximum aerodynamic force on the wing would exceed the ultimate load by 15.7%.This raises a serious concern for the structural integrity.The problem could be resolved either by the cost of the weight penalty or a load alleviation technology.
The gust response results for the four gust gradients H in different altitudes and cruise speeds VC and the design speed for maximum gust intensity VB were also calculated.A summary of the maximum wingtip displacements is listed in Table 5.From Tables 4 and 5 and Figure 7b, it is noted that the most critical condition for the aircraft in MTOM is the gust gradient H = 107 m at VC and sea level corresponding to the maximum gust velocity Uds and response.

The PTWT Key Parameters and Influence on Gust Response
The PTWT is a separate wing segment mounted at the wingtip through a shaft and torque spring connected to the wing front spar, as illustrated in Figure 8a.For the shaft located in front of the aerodynamic centre, the gust-induced aerodynamic force will cause the PTWT twist nose to be down, as illustrated in Figure 8a.This twist produces a negative AoA and a negative aerodynamic force on the PTWT that leads to gust alleviation of the wing.There are five key design parameters that affect the PTWT design and performance, the dimension, the shaft position in chordwise direction, the torque spring stiffness and the mass and location of the centre of gravity (CG).In this current study, a small PTWT of 1.13 m in the same length as the wingtip chord was taken.The rest of the three parameters will be determined subject to the design constraints and, in particular, the PTWT twist angle and aeroelastic stability.The basic requirement for the torque spring is set to be stiff enough to maintain the designed aerodynamic angle of 3.7° under the 1g level flight lift in cruise.Under extra lift due to gust, the maximum twist angle of the PTWT is limited to −10° to avoid stall.The shaft working as the PTWT front spar and also the spring supporter is of adequate bending and torsion stiffness.
In the previous study on a similar passive gust alleviation device for a flying wing aircraft [7], the torque spring stiffness between 10 and 40 kNm and the shaft location in the range of 5%-35% chord was evaluated.It was noted that the required spring stiffness was greater due to the larger wing span of 63 m and the longer device of 1.85 m compared with the aircraft in the current study.In the previous study, the shaft location was set at 5% chord from the leading edge for the maximum gust alleviation effect.
In the current PTWT design, the shaft connected to the front spar was located at 10% chord from the leading edge.The spring stiffness was considered in the range of 0.1-10 kNm/rad.The PTWT is made of sandwich composite skin reinforced by the spar and has a total mass of 30 kg.Since the PTWT mass is mainly in the front of the mid-chord dominated by the shaft and spring, its CG is located at 5% chord behind the shaft (15% from leading edge).Based on the initial design, an FE model of the PTWT was created at the wingtip to replace the original fixed wingtip.In the gust response analysis, the control surfaces were locked in neutral position.The gust response of the aircraft in terms of wingtip deflection was calculated in the most critical condition within a specified range of PTWT spring stiffness to determine the minimum stiffness.For the case of the wing clamped at the root, the gust response results are presented in Figure 8b.Within the range of spring stiffness from K = 0.1 to 10 kNm/rad, the PTWT torsional frequency increases from 0.034 to 0.173 Hz.This result indicates that coupling can be avoided between the PTWT torsion mode and the first bending mode of the wing.Figure 8b also shows that the hardest spring of 10 kNm/rad would result in a gust response of 1.57 m of wingtip deflection, which is close to the result for the baseline aircraft, as shown in Figure 7b.The minimum gust response of 1.18 m deflection is associated with the softest spring with the PTWT twist angle restricted to −10°.
The aerodynamic model of the PTWT was also created as a separate surface at the wingtip in the gust response analysis.Figure 8b shows that the gust response is most sensitive to the spring stiffness between 1 and 2 kNm/rad.For a stiffness of 2 kNm/rad or greater, the PTWT twist angle in response to the specified gust is below −5°.In this case, the gust alleviation can reach 7.5% with the maximum wingtip deflection reduced to 1443 from 1560 mm for the baseline aircraft, as shown in Figure 7b.If the minimum stiffness of 0.1 kNm/rad were taken, the PTWT twist angle in response to the gust would reach −10°.In this case, 25% gust alleviation can be obtained with the wingtip deflection reduced from 1560 to 1170 mm.To keep the PTWT in neutral position at cruise speed, however, the minimum stiffness of 1 kNm/rad is required and, therefore, used in the PTWT design.The spring stiffness of 1 kNm/rad can be obtained by using a steel helical spring of five coils with coil diameter of 20 mm and a spring diameter of 100 mm.
To ensure that the torsional stiffness and mass satisfy the aeroelastic constraint, flutter analysis of the wing with PTWT using different parameters was performed.Figure 9a,b shows the results of flutter speed variation within a range of stiffness and mass, respectively.A critical flutter speed of 283 m/s, which is 15% greater than the dive speed, is marked to indicate the aeroelastic design requirement according to CS-25.The results show that the flutter speed is 300 m/s for the wing with the PTWT of 30 kg mounted on the shaft with a spring stiffness of 1 kNm/rad.This PTWT design satisfies the aeroelastic stability requirement.Figure 9a,b also indicates that the flutter speed will be increased by either increasing the torsional stiffness or reducing the mass of the PTWT based on the current design.Based on the PTWT design with a spring stiffness of 1 kNm/rad, the gust response has been calculated in the four gust gradient cases and different flight conditions.The results of the gust response reduction of the wing achieved by employing the PTWT are listed in Table 6.From the results, it is noted that the PTWT is generally more effective at larger gust gradients or lower gust frequencies.

Gust Response of the Whole Aircraft
In the above study, the wing was clamped at the root, so that the analysis was focused on the PTWT effectiveness during the wing loading without the aircraft body motion.In such a case, the gust response is relatively large and overestimated compared to real life.To calculate the gust response of the whole aircraft, the fuselage and tail are modelled as a rigid body with free motion.For symmetric motion, the degrees of freedom of the aircraft model are set free in transverse and pitching directions.The gust response was measured in terms of displacement at the aircraft CG and the wingtip deflection relative to the CG.Firstly, the gust response of the baseline aircraft was calculated in the most critical gust case, as described in the previous section, with the result shown in Figure 10.Comparing with the previous results of the wing alone absorbing all of the gust energy, as shown in Figure 7b, the gust load on the wing was largely reduced with the wingtip deflection reduced by 48% from 1.59 to 0.82 m, as shown in Figure 10b.This is because the gust energy was absorbed by the wing elastic deformation and also the potential energy of the aircraft with an upward heave motion over 8 m, as shown in Figure 10a.Comparing with the wingtip deflection and strain under the static ultimate load, this result indicates that the gust load in real lift actually reaches 46% of the ultimate load and puts the structural integrity at risk.
When the PTWT was employed, the gust-induced wingtip deflection relative to the CG was calculated and is shown in Figure 11a.In comparison with the baseline results in Figure 10, it can be seen that a significant reduction of wingtip deflection by 21.2% from 0.83 to 0.65 m has been achieved.Figure 11b shows the elastic twist angle of the wing and the PTWT twist angle reaching almost −9° relative to the wing in response to the gust load.
Based on the analysis, the bending moment at the wing root was also reduced significantly by 10.6%.A summary of the aircraft gust alleviation results at sea level for both the MTOM and OEM cases by employing the PTWT in the same most critical gust H = 107 m, described in the previous section, is presented in Table 7.The results indicate that the gust-induced extra load on the wing can be reduced from 46% of the ultimate load (1.75 g) to 36.5% (1.39 g) by employing the PTWT.Comparing with the baseline aircraft where the maximum load in the most critical gust exceeded the design ultimate load by 15.7%, the PTWT has reduced the maximum aerodynamic load below the ultimate load.

The PTWT Influence on the Flight Manoeuvre Performance
Due to the nature of the PTWT in gust alleviation, it could have a negative influence on the flight performance when manoeuvring is required by operating the control surfaces.An investigation was carried out on the evaluation of the aircraft symmetric and unsymmetrical manoeuvrability with the PTWT unlocked.
In the symmetric manoeuvre simulation, a step input from the aircraft elevator was considered.In response to the elevator rotation, the aircraft would produce an incidence change that can be expressed by an exponential decay time function = 1 −  , where a decay constant λ = 4 was taken [20].For this aircraft, the angle of incidence converged to 5° after 1 s time associated with a steady pitch rate, as shown in Figure 12a.The corresponding variation of the PTWT twist angle over 4 s period is shown in Figure 12b.In response to the elevator input and the PTWT reaction, the normal velocity of the aircraft in MTOM increased from the level flight speed and converged to a stable value, as shown in Figure 13a.Comparing with the baseline aircraft, the normal velocity was reduced by 1.7% due to the PTWT-induced drag.The corresponding reduction of acceleration was 2.7% only.Figure 13b shows that the aircraft heave motion as a result of the incidence change was reduced by only 3% due to the PTWT.The normal velocity converged to a constant value, and the heave motion kept increasing, because no reverse elevator input was set to recover the aircraft status.The results indicate that the PTWT influence on the symmetric manoeuvrability is negligible.This is mainly because the PTWT twist angle, as shown in Figure 12b, has little influence on the increase of the aircraft incidence and lift due to the elevator action, as shown in Figure 12a.A summary of the PTWT influence on the symmetric manoeuvring measured as the percentage reduction of the aircraft normal acceleration, velocity and movement is listed in Table 8.In the unsymmetrical rolling manoeuvre, an input aileron rotation at a rate of 6°/s was considered in an exponential time function, as shown in Figure 14a.In response to the aileron rotation, the PTWT had a twist angle up to −6° in the opposite direction, as shown in Figure 14a.Due to the PTWT effect, the rolling moment of the aircraft is reduced by 20.5% compared with the baseline aircraft, as shown in Figure 14b.Consequently, the roll rate and the resulting roll angle were reduced by 11.6% and 10.6%, respectively, as shown in Figure 15.It is noted that the significant reduction of the roll moment is due to the negative effect of PTWT on both wings.If the PTWT nose-up rotation were restricted, the PTWT influence on the wing with the aileron rotating upward would be eliminated and the influence on the aircraft rolling moment would be reduced significantly.In addition, the associated roll acceleration was reduced by 13%.A summary of the PTWT influence on the unsymmetrical manoeuvring measured as the percentage reduction of the aircraft roll acceleration, rate and angle after four seconds is listed in Table 9.

Table 9.
The PTWT influence on the aircraft unsymmetrical manoeuvring.

Conclusions
The feasibility and effectiveness of employing a passive twist wingtip (PTWT) for gust alleviation has been evaluated through a case study of a 200-seater passenger aircraft with a high aspect ratio wing in this investigation.The baseline airframe, especially the composite wing structure, has been modelled and analysed at the detailed design level subject to the design criteria in line with CS-25.In the most critical gust load case, the maximum aerodynamic load on the wing can exceed the design ultimate load by 15.7% and cause concern for the structural integrity.In the PTWT design, the key parameters were determined based on the gust response and flutter analysis of the wing to achieve maximum effectiveness and aeroelastic stability.
The results show that significant gust response alleviation can be achieved by employing the PTWT in the whole range of gust gradients specified in the airworthiness standards.In the most critical case, the gust-induced wingtip deflection and bending moment at the wing root can be reduced by 21% and 11%, respectively.Consequently, the maximum load, including the worst gust load, on the wing can be reduced below the design ultimate load by employing the PTWT.
The PTWT influence on the symmetric manoeuvring of the aircraft in terms of normal velocity and heave motion is negligible.However, the PTWT influence on unsymmetrical manoeuvring of the aircraft raises a concern with the roll moment and roll rate reduced by 20.5% and 12%, respectively.Since the PTWT reduces the rolling control effectiveness, it makes the system more conservative rather than detrimental.To offset the negative effect, the input aileron rotation rate or angle need to be increased to enhance the performance.Alternative measures, such as a PTWT locking system, which is not an ideal option, could also be used during manoeuvring operation.The investigation has shown that the PTWT is a reliable and effective means for gust alleviation and has potential application for such a large civil aircraft.

Notation a1
Lift parameters MTOM, MLM and MZFM represent the maximum take-off mass, maximum landing mass and maximum zero fuel mass of the aircraft, respectively.

Figure 1 .
Figure 1.(a) The aircraft configuration and (b) the aerodynamic force distributions in three cases.PTWT, passive twist wingtip.

Figure 2 .
Figure 2. (a) Aerodynamic pressure over a 2D wing section and (b) the whole 3D wing.

Figure 3 .
Figure 3. (a) The most critical shear force and (b) bending moment distribution along the wing span.

Figure 4 .
Figure 4. (a) Structural layout of the aircraft and (b) the cross-section view of the wing box.

Figure 5 .
Figure 5. (a) Skin cover thickness and (b) mass breakdown of outer wing box from the initial design.

Figure 6 .
Figure 6.(a) The failure index (FI) results of the upper skin; and (b) strain plot of the lower and upper skins of the wing.

Figure 7 .
Figure 7. (a) The flutter V-g plot for MTOM and (b) gust response at four gust gradients in the MTOM case.

Figure 8 .
Figure 8.(a) The PTWT and (b) the wing gust response for a range of torsion stiffness.

Figure 9 .
Figure 9. (a) The wing flutter speed vs. the PTWT shaft torsion stiffness and (b) the mass.

Figure 10 .
Figure 10.(a) Gust response at the centre of gravity (CG) and (b) the wingtip deflection of the baseline aircraft.

Figure 11 .
Figure 11.(a) The wingtip deflection with and without PTWT and (b) the wing and PTWT twist angle.

Figure 12 .
Figure 12.(a) The aircraft incidence and (b) the PTWT twist angle in response to an elevator step input.

Figure 13 .
Figure 13.(a) The normal velocity and (b) heave motion of the aircraft with and without the PTWT.

Figure 14 .Figure 15 .
Figure 14.(a) The input aileron rotation angle and (b) the corresponding rolling moment of the aircraft.

Table 1 .
Technical data for the aircraft design.

Table 3 .
The first a few frequencies in the full fuel and empty fuel cases.

Table 4 .
The design gust velocity and frequency against the gust gradient distances.

Table 5 .
The gust response in different flight conditions.

Table 6 .
The gust alleviation achieved by the PTWT in different flight conditions.

Table 7 .
The gust response reduction with the PTWT against the baseline aircraft.

Table 8 .
The PTWT influence on the aircraft symmetric manoeuvring.