Wing Structural Model for Overall Aircraft Design of Distributed Electric Propulsion General Aviation and Regional Aircraft
Abstract
:1. Introduction
2. Materials and Methods
2.1. FAST-OAD-GA (Future Aircraft Sizing Tool—Overall Aircraft Design—General Aviation)
2.2. Wing Mass Model
- The spar web withstands the shear force;
- The spar caps supports the bending mometn;
- The skin reacts to the torsion moment;
- The ribs bear the buckling loads.
2.3. Methodology Description
- Load factors: CS 23.321 and CS 23.333 [15] state that the aircraft must comply with the load requirements for any combination of load factor and airspeed within the limits of the flight envelope. This envelope results from the application of the maneuver and gust criteria defined for every weight and flight altitude computed according to CS 23.335, CS 23.337, and CS 23.341. To simplify the process, two constraining weight cases for the wing structure at the typical cruise altitude were selected: the aircraft at MTOW and with minimum fuel in the wings [15]. The choice of considering only the cruise altitude and cruise equivalent airspeed was made since these conditions lead to the highest gust load factor. In addition to that, according to CS 23.305, the wing should also demonstrate compliance to the strength requirements at ultimate loads. These are calculated by applying a safety factor to the maximum/minimum load factors obtained for the aforementioned conditions.
- Lift distribution: The computation of the lift distribution was separated into two parts: the clean lift, which is the lift without accounting for engine wake and the slipstream lift, meaning the additional lift generated by the engine wake. The clean lift distribution is obtained from OpenVSP and was assumed to be of constant shape along the wingspan. It is then scaled to match the lift coefficient required to sustain the extreme load factors. In the computation of the clean lift, the tail down force necessary to maintain the aircraft level was also considered, leading to a increase assumed to balance the weight in cruise (L = W) [12]. The slipstream lift distribution is computed using OpenVSP through a simple actuator-disk engine model to capture the effects of DEP on the aerodynamic loads. For the computation of the slipstream effect it was considered that two types of propellers could be installed on the wing: the High-Lift Propeller (HLP) [16] to provide extra lift at low speed and the Wingtip Propeller (WTP) to provide thrust during cruise and reduce induced drag. In cruise, the HLPs are folded, but at the climb cruise transition, they are still activated and working close to their maximum operating speed. At this time of flight, the potential risk of whirl flutter require setting the torque at the minimum level with a consequent reduction of power [17]. These few seconds of blown wing under cruise conditions are taken into account in the wing sizing process under a steady approximation to ensure conservative results. This lift distribution is then multiplied by the chord at the corresponding station along the span to obtain the aerodynamic loads at each point of the wing. The results of the computation of the lift load for the X-57 are shown in Figure 5.
- Structural load diagrams: five types of loads acting on the wing are considered: the aerodynamic lift load previously obtained, the distribution of the wing and fuel weights (proportional to the local chord), the punctual weight of the engines and the weight of the landing gear. However, when the fuel is placed in the same position as an engine or landing gear housing, its weight is decreased at this section to account for the reduced space available. The shear force diagram is then obtained by the numerical integration of the distributed loads. Similarly, the bending moment diagram can be constructed with the numerical integral of the shear diagram. Pursuing the flexibility of this methodology for hybrid or DEP configurations, the fuel weight is only considered in the applicable cases.
- Mass breakdown of the wing: the mass of each structural component of the wing is proportional to the quantity of material necessary to support the loads, defined as the density of the material multiplied by the volume of the part. The following computations will be done for the two previously mentioned aircraft masses: MTOW and minimum fuel in the wing. Consequently, two values of the mass will be obtained for each element, with the maximum being retained. Hereunder, a non-dimensional reference coordinate system will be used, where the Y axis () is defined positive to the right side of the aircraft, with b the span.
- Web mass in shear : the spar web is a vertical panel designed to react, under assumptions of the presented methodology, to all the shear force of the wing . CS 23.305 stress criteria is applied, which states that the stress on the web must be lower than the maximum shear stress of the material . So, the mass of the web for the whole wing is calculated integrating its section area along its elastic axis which sweep angle is :
- Spar caps mass in bending : the lower and upper spar caps of the main spar react to the bending moment . Similarly to the web sizing, the CS 23.305 stress criteria imposes that the stress on the spar caps cannot surpass the maximum normal stress of the material, , neither in traction () nor in compression (). In general, the resultant shear force of the wing is in the upwards direction. Hence, the upper spar cap will work under compression, while the lower spar cap will be subjected to traction. The consequence is that the upper spar cap will face higher risk of failure under buckling than the lower spar cap [10]. With those considerations, the mass of both spar caps can be calculated with Equation (2), where is the distance between the center of gravity of the spar caps retrieved from [10] and , are the cross-section areas of the upper and lower spar caps, respectively.
- Skin mass in torsion : The skin is a thin panel of a thickness assumed constant that provides the necessary strength to bear the torsion moment. This assumption makes sense from a manufacturing point of view, where varying the thickness for wings of this size represents additional manufacturing complexity that is not desirable for general aviation manufacturers. Besides having a varying skin thickness did not significantly improve the results in [10], hence our choice to keep the wing skin thickness constant. The regulation imposes two sizing criteria: CS 23.455 that limits the deformation of the wing in rolling conditions and one that corresponds to the maximum admissible stress at ultimate loads stated in the CS 23.305. Hence, the skin mass will be determined by the most demanding requirement. In previous versions of the model [18], both criteria were computed for various aircraft. The results revealed that the deformation dimensioning also meets the stress requirements. Consequently, only the CS 23.455 rule was implemented in the current edition of the code. Then, as presented in Equation (4), the skin mass is obtained by integrating the product of the skin thickness, e and the perimeter of the wing-box, l. The objective of this criterion is to guarantee the efficiency of the ailerons control by limiting the roll speed at the design airspeeds . When ailerons are deflected, this norm restricts the wing deformation that changes the angle of attack . This can be approximated by the variation of the twist angle due to the torsion moment produced by the ailerons. The skin weight will then be computed based on the necessary thickness to react that torsion moment. For further detail on the calculation of e, the user is encouraged to see [10] (pp. 84–104).
The mass of the ribs and miscellaneous items are obtained empirically according to the expressions proposed by É. Roux [10] and NASA FLOPS [7], respectively. As for the secondary structure, Torenbeek [2] states that it contributes to around a of the total mass of the wing. Therefore, its weight will be empirically obtained with this proportion. Finally, the total mass of the wing will be the sum of all its components: - Iteration: the method finishes by comparing the computed mass of the wing with the one initially estimated. If these values are different, the tool iterates on the solution with the computed value until the error is lower than . This ensures the obtained wing weight is close enough to the structural weight used to compute relief effects.
2.4. Methodology Validation
2.5. Parametric Analysis
2.5.1. Effect of the Engines Weight Distribution
- The high-lift engines are evenly spaced along the wingspan and it is assumed that the distance between the last HLP and the WTP was greater than the distance between two HLP. The position used for this study is detailed in Table 1.
- The maximum power and mass of the engines are proportionally scaled with the number of HLP. The propeller disc size is also adjusted to have a constant power coefficient. The baseline is defined by the 12 high-lift motors mounted in the X-57 Maxwell, with a propeller diameter of 0.58 m, a shaft power of 10.5 kW and a mass of 6.8 kg per engine.
- The slipstream effect of HLP is not considered for this study. This analysis is focused on the engines weight and not on the aerodynamics though the two parameters will be coupled in an overall aircraft analysis.
- The WTP remains the same in every simulation. This enables to maintain the main source of thrust of the aircraft to counteract the drag.
2.5.2. Effect of Slipstream Consideration
2.5.3. Effect of Engine Placement
3. Results
3.1. Methodology Validation
3.2. Parametric Studies
3.2.1. Effect of the Engines Weight Distribution
3.2.2. Effect of Slipstream Consideration
3.3. Effect of Engine Placement
3.4. Computational Time
4. Discussion
5. Conclusions
Supplementary Materials
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Acknowledgments
Conflicts of Interest
Abbreviations
AR | Aspect Ratio |
DEP | Distributed Electric Propulsion |
EASA | European Union Aviation Safety Agency |
FAST | Future Aircraft Sizing Tool |
GA | General Aviation |
HLP | High-Lift Propeller |
MTOW | Maximum Take-Off Weight |
OAD | Overall Aircraft Design |
OEW | Operating Empty Weight |
OEI | One Engine Inoperative |
WTP | WingTip Propeller |
Appendix A. Statistical Formula Used for Methodology Validation
Appendix B. Reference Technical Data for the X-57 Maxwell and the Beechcraft 76
Parameter | Unit | X-57 Maxwell | Beechcraft 76 | Beechcraft 1900 |
---|---|---|---|---|
MTOW | kg | 1360 | 1769 | 7688 |
MZFW | kg | 1360 | 1381.7 | 6804 |
Wing mass | kg | 166.7 | - | - |
Design fuel | kg | 0 | 387.3 | 884 |
Parameter | Unit | X-57 Maxwell | Beechcraft 76 | Beechcraft 1900 |
---|---|---|---|---|
Fuselage width | m | 1.22 | 1.20 | 1.38 |
Wing span | m | 9.66 | 11.58 | 17.67 |
Aspect ratio | - | 15 | 7.98 | 10.84 |
Wing surface | m | 6.2 | 16.8 | 28.8 |
Root chord | m | 0.74 | 1.48 | 2.22 |
Root thickness | m | 0.088 | 0.21 | 0.4 |
Tip chord | m | 0.52 | 1.42 | 0.93 |
Leading edge sweep | deg | 1.9 | 0.0 | 0.0 |
Aileron chord | %c | 25 | 26 | 0.28 |
Aileron span | %b/2 | 22 | 32 | 0.38 |
Aileron max deflection | deg | 22 | 20 | 20 |
Airfoil thickness ratio | - | 0.12 | 0.15 | 0.18 |
Ribs step | m | 0.6 | 0.6 | 0.6 |
Parameter | Unit | X-57 Maxwell | Beechcraft 76 | Beechcraft 1900 |
---|---|---|---|---|
Cruise altitude | ft | 8000 | 8000 | 20,000 |
Cl | rad | 6.1 | 5.55 | 7.0 |
Cl | - | 4.2 | 1.97 | 2.10 |
Parameter | Unit | X-57 Maxwell | Beechcraft 76 | Beechcraft 1900 |
---|---|---|---|---|
V | m·s | 78.19 | 78.63 | 118.27 |
V | m·s | 97.7 | 110.0 | 163.02 |
V | m·s | 58.17 | 64.75 | 90.38 |
Parameter | Unit | X-57 Maxwell | Beechcraft 76 | Beechcraft 1900 |
---|---|---|---|---|
Cruise engines | nº | 2 | 2 | 2 |
HLP engines | nº | 12 | 0 | 0 |
Cruise engine power | kW | 60.0 | 130.0 | 953 |
HLP power | kW | 10.5 | - | - |
Cruise engines weight | kg | 53.1 | 175.0 | 325.4 |
HLP weight | kg | 6.8 | - | - |
Cruise propeller dia | m | 1.5 | 1.93 | 2.78 |
HLP dia | m | 0.58 | - | - |
Parameter | Unit | X-57 Maxwell | Beechcraft 76 | Beechcraft 1900 |
---|---|---|---|---|
q | lb/ft | 61.48 | 72.44 | - |
V | KEAS | 150 | 153 | - |
N | - | 6.12 | 5.7 | 4.5 |
References
- Moore, K.R.; Ning, A. Distributed electric propulsion effects on existing aircraft through multidisciplinary optimization. In Proceedings of the 2018 AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Kissimmee, FL, USA, 8–12 January 2018; p. 1652. [Google Scholar]
- Ko, A.; Schetz, J.A.; Mason, W.H. Assessment of the potential advantages of distributed-propulsion for aircraft. In Proceedings of the XVI International Symposium on Air Breathing Engines (ISABE), Cleveland, OH, USA, 31 August–5 September 2003. [Google Scholar]
- Gudmundsson, S. General Aviation Aircraft Design: Applied Methods and Procedures; Butterworth-Heinemann: Oxford, UK; Waltham, MA, USA, 2013. [Google Scholar]
- Raymer, D.P. Aircraft Design: A Conceptual Approach; AIAA Education Series; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2012. [Google Scholar]
- Nicolai, L.M.; Carichner, G.E. Fundamentals of Aircraft and Airship Design, Volume 1—Aircraft Design; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2010. [Google Scholar]
- Dorbath, F. A Flexible Wing Modeling and Physical Mass Estimation System for Early Aircraft Design Stages. Ph.D. Thesis, Technische Universität Hamburg, Hamburg, Germany, 2014. [Google Scholar]
- Wells, D.P.; Horvath, B.L.; McCullers, L.A. The Flight Optimization System Weights Estimation Method. 2017. Available online: https://ntrs.nasa.gov/api/citations/20170005851/downloads/20170005851.pdf (accessed on 17 November 2021).
- Bindolino, G.; Ghiringhelli, G.; Ricci, S.; Terraneo, M. Multilevel structural optimization for preliminary wing-box weight estimation. J. Aircr. 2010, 47, 475–489. [Google Scholar] [CrossRef]
- Torenbeek, E. Development and Application of a Comprehensive, Design-Sensitive Weight Prediction Method for Wing Structures of Transport Category Aircraft; Delft University of Technology, Faculty of Aerospace Engineering, Report LR-693; Delft University of Technology: Delft, The Netherlands, 1992. [Google Scholar]
- Roux, E. Modèle de Masse Voilure: Avions de Transport Civil. Ph.D. Thesis, SupAéro-ONERA, Toulouse, France, 2006. [Google Scholar]
- David, C.; Delbecq, S.; Defoort, S.; Schmollgruber, P.; Benard, E.; Pommier-Budinger, V. From FAST to FAST-OAD: An open source framework for rapid Overall Aircraft Design. In IOP Conference Series: Materials Science and Engineering; IOP Publishing: Bristol, UK, 2021; Volume 1024, p. 012062. [Google Scholar]
- Roskam, J. Airplane Design: Part 5—Component Weight Estimation; DARcorporation: Lawrence, KS, USA, 1985. [Google Scholar]
- Litherland, B.; Rieth, K. VSP Aircraft Analysis User Manual; NASA: Washington, DC, USA, 2014. [Google Scholar]
- An, S. Aeroelastic Design of a Lightweight Distributed Electric Propulsion Aircraft with Flutter and Strength Requirements. Ph.D. Thesis, Georgia Institute of Technology, Atlanta, GA, USA, 2015. [Google Scholar]
- EASA. Certification Specifications and Acceptable Means of Compliance for Normal, Utility, Aerobatic, and Commuter Category Aeroplanes CS-23 Amendment 4. July 2015. Available online: https://www.easa.europa.eu/sites/default/files/dfu/CS-23%20Amendment%204.pdf (accessed on 17 November 2021).
- Patterson, M.D. Conceptual Design of High-Lift Propeller Systems for Small Electric Aircraft. Ph.D. Thesis, Georgia Institute of Technology, Atlanta, GA, USA, 2016. [Google Scholar]
- Borer, N.K.; Patterson, M.D. X-57 High-Lift Propeller Control Schedule Development. In Proceedings of the AIAA AVIATION 2020 FORUM, Virtual Event, 15–19 June 2020; p. 3091. [Google Scholar]
- Ruscio, J.P.; Jezegou, J.; Benard, E.; Gomez Pacheco, A.; Laonet, P.; Alonso Castilla, R. Hybrid electric distributed propulsion overall aircraft design process and models for general aviation (FAST GA). In IOP Conference Series: Materials Science and Engineering; IOP Publishing: Bristol, UK, 2021; Volume 1024, p. 012072. [Google Scholar]
- Li, W. Overview of the X-57 Structural Requirements, Modifications, and Airworthiness. 2019. Available online: https://ntrs.nasa.gov/api/citations/20190026541/downloads/20190026541.pdf (accessed on 17 November 2021).
- Moore, M.D. Distributed Electric Propulsion (DEP) Aircraft; NASA Langley Research Center: Hampton, VA, USA, 2012. [Google Scholar]
- Deere, K.A.; Viken, J.K.; Viken, S.; Carter, M.B.; Wiese, M.; Farr, N. Computational analysis of a wing designed for the X-57 distributed electric propulsion aircraft. In Proceedings of the 35th AIAA Applied Aerodynamics Conference, Denver, CO, USA, 5–9 June 2017; p. 3923. [Google Scholar]
- Habermann, A.L. Effects of Distributed Propulsion on Wing Mass in Aircraft Conceptual Design. In Proceedings of the AIAA AVIATION 2020 FORUM, Virtual Event, 15–19 June 2020; p. 2625. [Google Scholar]
- Sinnige, T.; van Arnhem, N.; Stokkermans, T.C.; Eitelberg, G.; Veldhuis, L.L. Wingtip-mounted propellers: Aerodynamic analysis of interaction effects and comparison with conventional layout. J. Aircr. 2019, 56, 295–312. [Google Scholar] [CrossRef] [Green Version]
- Niu, C.; Niu, M.C. Airframe Structural Design: Practical Design Information and Data on Aircraft Structures; Conmilit Press: Hong-Kong, China, 1988. [Google Scholar]
- Megson, T.H.G. Aircraft Structures for Engineering Students; Butterworth-Heinemann: Oxford, UK; Waltham, MA, USA, 2016. [Google Scholar]
Position (s) of the HLP (s) | Position of the WTP | ||||||
---|---|---|---|---|---|---|---|
2 (1 + 1) | 0.33 | - | - | - | - | - | 1.0 |
3 (2 + 1) | 0.25 | 0.5 | - | - | - | - | 1.0 |
4 (3 + 1) | 0.2 | 0.4 | 0.6 | - | - | - | 1.0 |
5 (4 + 1) | 0.17 | 0.35 | 0.5 | 0.65 | - | - | 1.0 |
6 (5 + 1) | 0.143 | 0.286 | 0.429 | 0.572 | 0.715 | - | 1.0 |
7 (6 + 1) | 0.188 | 0.3075 | 0.4269 | 0.5462 | 0.667 | 0.785 | 1.0 |
Estimation Method | Beechcraft 76 | X-57 Maxwell | Beechcraft 1900 | |||
---|---|---|---|---|---|---|
Wing Mass (kg) | Error (%) | Wing Mass (kg) | Error (%) | Wing Mass (kg) | Error (%) | |
Developed model | 182.2 | 0.4 | 174.3 | 4.5 | 722.8 | 0.7 |
Statistical model | 183.0 [4] | Reference | 131.8 [5] | 20.1 | 717.5 [12] | Reference |
Real value | Unknown | - | 166.7 [20] | Reference | Unknown | - |
Wing Structural Component | Beechcraft 76 | X-57 Maxwell | Beechcraft 1900 |
---|---|---|---|
Web | 2.235 kg (1.2%) | 2.409 kg (1.4%) | 15.538 kg (2.1%) |
Lower spar cap | 10.889 kg (6.0%) | 20.549 kg (11.8%) | 64.643 kg (8.8%) |
Upper spar cap | 14.548 kg (8.0%) | 27.455 kg (15.7%) | 86.365 kg (11.8%) |
Skin | 49.73 kg (27.3%) | 68.501 kg (39.3%) | 267.065 kg (36.0%) |
Ribs | 22.128 kg (12.1%) | 3.206 kg (1.8%) | 47.643 kg (6.6%) |
Misc | 37.113 kg (20.4%) | 8.64 kg (5.0%) | 70.863 kg (9.7%) |
Primary structure | 136.643 kg (75.0%) | 130.76 kg (75.0%) | 542.117 kg (75.0%) |
Secondary structure | 45.548 kg (25.0%) | 43.587 kg (25.0 %) | 180.705 kg (25.0%) |
Total | 182.191 kg (100.0%) | 174.347 kg (100.0%) | 722.821 kg (100.0%) |
2 (1 + 1) | 2.34 | 47.10 | 68.50 | 3.21 | 8.64 | 43.26 | 173.05 |
3 (2 + 1) | 2.32 | 46.40 | 43.02 | 172.10 | |||
4 (3 + 1) | 2.31 | 45.95 | 42.87 | 171.48 | |||
5 (4 + 1) | 2.31 | 45.66 | 42.77 | 171.08 | |||
6 (5 + 1) | 2.30 | 45.37 | 42.67 | 170.69 | |||
7 (6 + 1) | 2.28 | 44.58 | 42.40 | 169.6 |
w/slipstream | 2.41 | 48.04 | 68.50 | 3.21 | 8.64 | 43.59 | 174.35 |
w/o slipstream | 2.28 | 44.58 | 42.40 | 169.6 |
Semi-Analytical Model | Statistical Model | |
---|---|---|
Beechcraft 76 | 76 ms | Virtually instantaneous |
X-57 Maxwell | 155 ms | |
Beechcraft 1900 | 87 ms |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Alonso Castilla, R.; Lutz, F.; Jézégou, J.; Bénard, E. Wing Structural Model for Overall Aircraft Design of Distributed Electric Propulsion General Aviation and Regional Aircraft. Aerospace 2022, 9, 5. https://doi.org/10.3390/aerospace9010005
Alonso Castilla R, Lutz F, Jézégou J, Bénard E. Wing Structural Model for Overall Aircraft Design of Distributed Electric Propulsion General Aviation and Regional Aircraft. Aerospace. 2022; 9(1):5. https://doi.org/10.3390/aerospace9010005
Chicago/Turabian StyleAlonso Castilla, Raquel, Florent Lutz, Joël Jézégou, and Emmanuel Bénard. 2022. "Wing Structural Model for Overall Aircraft Design of Distributed Electric Propulsion General Aviation and Regional Aircraft" Aerospace 9, no. 1: 5. https://doi.org/10.3390/aerospace9010005
APA StyleAlonso Castilla, R., Lutz, F., Jézégou, J., & Bénard, E. (2022). Wing Structural Model for Overall Aircraft Design of Distributed Electric Propulsion General Aviation and Regional Aircraft. Aerospace, 9(1), 5. https://doi.org/10.3390/aerospace9010005