Propeller Effects and Elasticity in Aerodynamic Analysis of Small Propeller-Driven Aircraft

1. Introduction

Currently, there is an increasing interest in daily passenger and cargo transport by small aircraft over distances that have been traditionally served by ground transportation vehicles. Although the idea of using small aircraft for transportation is not novel, over recent years there have been significant enhancements in technologies and societal variations that may make these operations become a practical part of our daily life. The convergence of advanced technologies, such as electric propulsion and autonomy, along with new business models, such as mobile application-based ride sharing and network-enabled on-demand services, bring new general aviation markets [1,2,3]. Advanced air mobility (AAM) is the general term defined by NASA for these new small aircraft. AAM services fall into two categories of urban air mobility (UAM) and regional air mobility (RAM) [4,5]. UAM is the transport system that transfers people or goods by air within cities. RAM will be responsible for the transportation of goods and people to rural and remote communities [6]. AAM has the potential to enhance our daily transportation by reducing the travel time, avoiding traffics on the ground and within the cities, smooth the travel between different parts of a city, and advance regional mobility, particularly in areas underserved by the modern air transportation [7].

The propeller effects and the power contribution to the aerodynamics of a small aircraft are essential. From the beginning of the airplane introduction, this has always been important to realize the propeller behavior due to the generation of substantial lateral forces in the presence of side wind [8,9]. Recently, much research has enhanced semi-empirical procedures to study the aerodynamics of sophisticated propeller aircraft and the propeller effects on the stability and handling quality [10,11,12,13]. These research presented an augmented and combined analytical procedures and design data compendia to account for the restrictions of the earlier approaches [8]. Among those, the semi-empirical multidisciplinary analysis program, MAPLA has been developed for the optimization of small conventional and eVTOL aircraft. MAPLA consists of five primary disciplines of weight and balance, aerodynamics, propulsion, stability, and control as well as performance. In a series of previous studies, MAPLA was demonstrated to be able to model the characteristics of small conventional and eVTOL aircraft with an acceptable precision [10,14,15,16].

In the design process of new airplanes, the performance of wing is dominated by the aero-structural-coupled system [17]. Due to the deflection of wings in flight which results from aerodynamic loads that are transferred to the wing structure, the aerodynamic analysis of the wing in flight is different than the rigid wing analysis [18,19]. In order to simulate this behavior, the general approach is to couple the aerodynamic solver and the structural computations. This coupling would need an iteration loop to find the equilibrium solution [20,21,22].

Among aero-structural coupling analysis tools, NASTRAN is a finite element analysis program that was originally developed for NASA [23]. In aero-structural coupling problems, structural and aerodynamic grids are connected by interpolation [24]. This approach enables the independent selection of grid points for both the structural model and the aerodynamic elements of lifting surfaces or bodies. The structural model for a wing may encompass a one-, two-, or three-dimensional array of grid points, while the aerodynamic theory may employ lifting surface theory or strip theory. A versatile interpolation method is available to connect various combinations seamlessly. Any aerodynamic panel or body can be subdivided into subregions for interpolation, utilizing separate functions for each. This interpolation method, known as splining, involves mathematical analysis of beams and plates, including linear splines that generalize infinite beams to accommodate torsional and bending degrees of freedom, surface splines that address solutions for infinite uniform plates, and explicit user-defined interpolation methods [25,26].

Static (or quasi-steady) aeroelasticity investigation is an interdisciplinary field that combines knowledge of aerodynamics, elasticity, and inertial forces [27,28,29]. By understanding and addressing static aeroelastic effects, engineers can ensure the stability and structural integrity of aircraft under various operating conditions, contributing to safer and more reliable flight operations [30]. Structural analysts are primarily concerned with the aerodynamic load redistribution and the resulting internal structural load and stress redistributions [31,32]. They also consider the possibility of static aeroelastic instability, such as divergence. Aerodynamicists and control systems analysts are interested in aerodynamic load redistribution and its effects on aerodynamic stability and control derivatives. The static aeroelastic capability in Nastran addresses these needs by computing aircraft trim conditions, followed by the recovery of structural responses, aeroelastic stability derivatives, and static aeroelastic divergence dynamic pressures [25,33].

The primary objective of this paper is to investigate the significance of propeller effects and elasticity in the aerodynamic analysis of small propeller-driven airplanes. Subsequently, the methodology section elaborates on the proposed procedure in detail, offering a comprehensive insight into the analytical approach. The results pertaining to a twin-engine propeller-driven aircraft are presented meticulously, serving as a case study to illustrate the application of the proposed methodology.

2. Materials and Methods

The methodology section delineates two key analyses: firstly, an examination of propeller effects in rigid aircraft, followed by the integration of elastic wing properties to explore alterations in aerodynamic characteristics. The analysis of propeller effects in rigid aircraft serves as the foundational step, providing insights into the aerodynamic behavior under normal operating conditions. Subsequently, the incorporation of elastic wing properties enriches our understanding by exploring how these variables influence aerodynamic performance. Validation of the studied aircraft is accomplished through a comparative analysis with wind tunnel test results obtained from a similar twin-engine propeller-driven aircraft with a rigid structure. Furthermore, the aerodynamic analysis of the elastic aircraft is conducted utilizing NASTRAN software, complemented by semi-empirical findings derived from the rigid aircraft analysis.

2.1. Rigid Aircraft Aerodynamic Analysis with Propeller Effects

The semi-empirical analysis program MAPLA, initially designed for optimizing small, general aviation aircraft, underwent enhancements for increased efficacy. Its original implementation encompasses five crucial disciplines: aerodynamics, propulsion, performance, weight and balance, as well as stability and control, amalgamating cutting-edge analytical procedures and design data collections into a fully automated method. For this investigation, the propulsion module specifically computes the propeller and power effects on longitudinal and lateral-directional aerodynamic coefficients and stability derivatives of rigid aircraft. Within each subprogram, power-on static stability and control derivatives were initially estimated across various aircraft components, including the wing, fuselage, nacelle, horizontal tail, vertical tail, and high-lift surfaces. The total power-on static stability and control derivatives were subsequently derived by amalgamating these individual contributions. Following this, the dynamic characteristics of the aircraft were estimated using static derivatives. The propulsion module's development was rooted in NASA's work, with enhancements tailored to generalize its application, particularly for small aircraft design and development purposes. In a series of previous investigations, this tool demonstrated its capability to model small aircraft characteristics with commendable precision [8,12,14,15,16,34].

The propeller effects on the lift forces could be divided into two groups such as those created from propeller forces and those created from the propeller slipstream. Hence the aircraft lift could be presented by [34,35]

where propeller forces components are from propeller thrust vector, ${\left(∆{\mathrm{C}}_{\mathrm{L}}\right)}_{\mathrm{T}}$ and propeller normal force, ${\left(∆{\mathrm{C}}_{\mathrm{L}}\right)}_{{\mathrm{N}}_{\mathrm{p}}}$. And propeller slipstream effects contain the lift due to the power-induced change in dynamic pressure, ${\left(∆{\mathrm{C}}_{{\mathrm{L}}_{\mathrm{h}}}\right)}_{∆{\overline{\mathrm{q}}}_{\mathrm{h}}}$and lift change because of the propeller downwash for wing, ${\left(∆{\mathrm{C}}_{\mathrm{L}}\right)}_{∆{\overline{\mathrm{q}}}_{\mathrm{w}}}+{\left(∆{\mathrm{C}}_{\mathrm{L}}\right)}_{{\in }_{\mathrm{p}}}$ and horizontal tail, ${\left(∆{\mathrm{C}}_{{\mathrm{L}}_{\mathrm{h}}}\right)}_{{(∆{\in }_{\mathrm{h}})}_{\mathrm{p}\mathrm{o}\mathrm{w}\mathrm{e}\mathrm{r}}}$.

The same could be done to present the pitching moment coefficient of the aircraft in the presence of propeller effects [34,35]

Where similar to the lift equation, propeller forces components are from propeller thrust vector, ${\left(∆{\mathrm{C}}_{\mathrm{m}}\right)}_{\mathrm{T}}$ and propeller normal force, ${\left(∆{\mathrm{C}}_{\mathrm{m}}\right)}_{{\mathrm{N}}_{\mathrm{p}}}$. And propeller slipstream effects contain the contribution of the power-induced change in dynamic pressure, ${\left(∆{\mathrm{C}}_{{\mathrm{m}}_0}\right)}_{∆{\overline{\mathrm{q}}}_{\mathrm{w}}}$and propeller slipstream-induced dynamic-pressure and angle of attack changes on the wing, ${{\left(∆{\mathrm{C}}_{\mathrm{m}}\right)}_{{\mathrm{w}}_{\mathrm{L}}}=\left(∆{\mathrm{C}}_{\mathrm{m}}\right)}_{∆{\overline{\mathrm{q}}}_{\mathrm{w}}}+{\left(∆{\mathrm{C}}_{\mathrm{m}}\right)}_{{\in }_{\mathrm{p}}}$ , propeller slipstream on nacelle free moments, ${\left(∆{\mathrm{C}}_{\mathrm{m}}\right)}_{{\mathrm{n}}_{\mathrm{p}}}$, and propeller slipstream on dynamic-pressure and downwash on the horizontal tail, ${\left(∆{\mathrm{C}}_{\mathrm{m}}\right)}_{\mathrm{h}}$.

The following components are necessary for the aircraft drag change due to the propeller effects. Firstly, the propeller thrust components parallel to the X-stability axis. Secondly, the change in slipstream stream dynamic pressure. Thirdly, the change in induced drag due to the lift component of the direct propeller forces. And finally, the change in cooling drag. These components are summarized as [34,35]

where $-\mathrm{n}({\mathrm{T}}_{\mathrm{c}}^{\mathrm{\text{'}}}/\mathrm{prop})\cos {\mathsf{\alpha }}_{\mathrm{T}}$ is the component of total thrust parallel to the velocity vector, $∆{\mathrm{C}}_{{\mathrm{D}}_0}$ is the change in profile drag coefficient due to power, $∆{\mathrm{C}}_{{\mathrm{D}}_{\mathrm{i}}}$ is the induced drag coefficient due to power and ${\left(∆{\mathrm{C}}_{\mathrm{D}}\right)}_{\mathrm{c}\mathrm{o}\mathrm{o}\mathrm{l}\mathrm{i}\mathrm{n}\mathrm{g}\mathrm{d}\mathrm{r}\mathrm{a}\mathrm{g}}$ is the change in the cooling system drag coefficient due to power.

The power effects on the side force derivative include the propeller normal force contribution, ${\left(∆{\mathrm{C}}_{{\mathrm{Y}}_{\mathsf{\beta }}}\right)}_{{\mathrm{N}}_{\mathrm{p}}}$, propeller-induced increase in dynamic pressure, ${\left(∆{\mathrm{C}}_{{\mathrm{Y}}_{\mathsf{\beta }}}\right)}_{\mathrm{n}(∆\overline{\mathrm{q}})}$and power-induced sidewash, ${\left(∆{\mathrm{C}}_{{\mathrm{Y}}_{\mathsf{\beta }}}\right)}_{\mathrm{n}\left({\mathsf{\sigma }}_{\mathrm{p}}\right)}$[34,35]

The propeller and power effects on the weathercock stability also include the propeller normal force contribution, ${\left(∆{\mathrm{C}}_{{\mathrm{n}}_{\mathsf{\beta }}}\right)}_{{\mathrm{N}}_{\mathrm{p}}}$, propeller-induced increase in dynamic pressure, ${\left(∆{\mathrm{C}}_{{\mathrm{n}}_{\mathsf{\beta }}}\right)}_{\mathrm{n}(∆\overline{\mathrm{q}})}$and power-induced sidewash, ${\left(∆{\mathrm{C}}_{{\mathrm{n}}_{\mathsf{\beta }}}\right)}_{\mathrm{n}\left({\mathsf{\sigma }}_{\mathrm{p}}\right)}$ [34,36].

Finally, the propeller and power effects on the dihedral derivative include the propeller side force contribution, ${\left(∆{\mathrm{C}}_{{\mathrm{l}}_{\mathsf{\beta }}}\right)}_{{\mathrm{N}}_{\mathrm{p}}}$ and the propeller-induced increase in dynamic pressure and downwash, ${\left(∆{\mathrm{C}}_{{\mathrm{l}}_{\mathsf{\beta }}}\right)}_{\mathrm{w}(∆\overline{\mathrm{q}}+{\mathsf{ϵ}}_{\mathrm{p}})}$ [36].

Change in dihedral derivative due to propeller effects and power on Propeller forces

2.1.1. Rigid Aircraft Analysis Validation

To validate the calculated results, the aerodynamics module outcomes were compared with wind tunnel test data of a twin-engine propeller-driven small aircraft [35,36]. Figure 1a and Figure 1b depict the geometry of the original aircraft model and the modelled aircraft using MAPLA, respectively. Additionally, Table 1 outlines the general characteristics of the small twin-engine propeller airplane and properties of the investigated flight condition.

Figure 1. The twin-engine propeller-driven small airplane under investigation is referenced from NASA’s report [35,36], denoted as a. for the reported aircraft and b. for the modeled aircraft utilizing MAPLA.

Figure 1. The twin-engine propeller-driven small airplane under investigation is referenced from NASA’s report [35,36], denoted as a. for the reported aircraft and b. for the modeled aircraft utilizing MAPLA.

In the subsequent sections, Figure 2a to Figure 2c delineate the longitudinal aerodynamic characteristics of the investigated aircraft across various flight conditions. Additionally, Figure 3a and Figure 3b compare the lateral-directional aerodynamic characteristics of the twin-engine propeller-driven small aircraft using MAPLA with the corresponding data obtained from available wind tunnel tests [35,36]. All results are presented from the power-off conditions where the trust coefficient of the propellers, CT is equal to zero up to CT=0.2 and CT=0.44. The trust coefficient parameter that is counting for propeller effects could be defined using:

$$\mathrm{C}\mathrm{T}=\frac{\mathrm{T}\mathrm{h}\mathrm{r}\mathrm{u}\mathrm{s}\mathrm{t}}{{\overline{\mathrm{q}}}_{\mathrm{\infty }}{\mathrm{S}}_{\mathrm{w}}}$$

(7)

Where ${\overline{\mathrm{q}}}_{\mathrm{\infty }}$ is the dynamic pressure ratio in newton per square meter and ${\mathrm{S}}_{\mathrm{w}}$ is the surface area of the wing in suare meter.

The deviation results of MAPLA from wind tunnel test results are presented in Table 2. With respect to the deviation results presented in Table 1, MAPLA’s results present an accuracy of 5% for lift coefficient values from power-off conditions to the case of CT=0.2 and even higher CT values of 0.44. The same could be seen for the power off conditions for drag coefficient results. However, as power increases the results’ deviation also increases up to 15% for CT=0.2 and 35% for CT=0.44. For the case of pitching moment coefficient results, MAPLA presented higher deviation from the wind tunnel test results starting from 10% for the power-off condition and increasing up to 20% and 50% for CT=0.2 and CT=0.44, respectively. On the other side, for lateral-directional characteristics, MAPLA reported higher deviation from wind tunnel test results where for the weather-cock stability analysis, the reported value for power-off conditions showed a deviation of 35% and this value increased for the higher power conditions up to 40% for CT=0.2 and up to 45% for CT=0.44.This is while, the dihedral effect characteristic showed a better accuracy by 25% deviation from wind-tunnel test results in CT=0, 15% in CT=0.2 and 10% in CT=0.44. Overall, MAPLA showed a better accuracy for longitudinal aerodynamic characteristics compared to the lateral-directional characteristics. The results indicated that MAPLA’s aerodynamics module is able to determine the aerodynamic characteristics of the small aircraft to acceptable accuracy for low fidelity analysis where conceptual and preliminary studies are of interest. However, compared to the other methods such as DATCOM [34], MAPLA’s results are closer to the wind tunnel test data in almost all studied conditions for linear angles of attacks.

2.2. Elastic Wing Integration

In this study, the aeroelastic analysis are based on the NASTRAN approach for static aeroelastic problems and deals with the interaction of aerodynamic and structural forces on a flexible vehicle that results in a redistribution of the aerodynamic loading as a function of airspeed [25]. The static aeroelastic study here addresses the static aeroelastic problem in aircraft trim condition.

Three matrix equations summarize the relationships required to define a set of aerodynamic influence coefficients [25,37]. These are the basic relationships between the lifting pressure and the dimensionless vertical or normal velocity induced by the inclination of the surface to the airstream; i.e., the downwash (or normalwash),

$$\left\{{\mathrm{w}}_{\mathrm{j}}\right\}=\left[{\mathrm{A}}_{\mathrm{j}\mathrm{j}}\right]\{{\mathrm{f}}_{\mathrm{j}}/\mathrm{q}\}$$

(8)

the substantial differentiation matrix of the deflections to obtain downwash,

$$\left\{{\mathrm{w}}_{\mathrm{j}}\right\}=\left[{\mathrm{D}}_{\mathrm{j}\mathrm{k}}^1+\mathrm{i}\mathrm{k}{\mathrm{D}}_{\mathrm{j}\mathrm{k}}^2\right]\left\{{\mathrm{u}}_{\mathrm{k}}\right\}+\left\{{\mathrm{w}}_{\mathrm{j}}^{\mathrm{g}}\right\}$$

(9)

and the integration of the pressure to obtain forces and moments,

$$\left\{{\mathrm{P}}_{\mathrm{k}}\right\}=\left[{\mathrm{S}}_{\mathrm{k}\mathrm{j}}\right]\left\{{\mathrm{f}}_{\mathrm{j}}\right\}$$

(10)

where ${\mathrm{w}}_{\mathrm{j}}$ denotes the downwash, ${\mathrm{w}}_{\mathrm{j}}^{\mathrm{g}}$ is the static aerodynamic downwash and includes the static incidence distribution that may arise from an initial angle of attack, camber, or twist. ${\mathrm{f}}_{\mathrm{j}}$ is the pressure on lifting element j and q is the dynamic pressure. $\mathrm{k}$ is the reduced frequency and equal to $\mathsf{\omega }{\mathrm{b}}_{\mathrm{s}}/\mathrm{V}$ where $\mathsf{\omega }$ is the angular frequency, ${\mathrm{b}}_{\mathrm{s}}$ is the reference semichord and $\mathrm{V}$ is the free-streem velocity. ${\mathrm{A}}_{\mathrm{j}\mathrm{j}}$ is the aerodynamic influence coefficient matrix, a function of Mach number and reduced frequency. ${\mathrm{u}}_{\mathrm{k}}$ and ${\mathrm{P}}_{\mathrm{k}}$ are displacements and forces at aerodynamic grid points, respectively. ${\mathrm{D}}_{\mathrm{j}\mathrm{k}}^1$ and ${\mathrm{D}}_{\mathrm{j}\mathrm{k}}^2$ real and imaginary parts of substantial differentiation matrix, respectively. And ${\mathrm{S}}_{\mathrm{k}\mathrm{j}}$ is the integration matrix.

The three matrices presented as equations 8 through 10 can also be combined to give an aerodynamic influence coefficient matrix [25]:

$$\left[{\mathrm{Q}}_{\mathrm{k}\mathrm{k}}\right]=\left[{\mathrm{S}}_{\mathrm{k}\mathrm{j}}\right]{\left[{\mathrm{A}}_{\mathrm{j}\mathrm{j}}\right]}^{-1}[{\mathrm{D}}_{\mathrm{j}\mathrm{k}}^1+\mathrm{i}\mathrm{k}{\mathrm{D}}_{\mathrm{j}\mathrm{k}}^2]$$

(11)

Structural and aerodynamic grids are connected by interpolation. This allows the independent selection of grid points of the structure and aerodynamic elements of the lifting surfaces/bodies in a manner best suited to the particular theory. The interpolation method is called splining. The structural degrees of freedom are chosen as the independent degrees of freedom and the aerodynamic degrees of freedom are dependent. The splining methods lead to an interpolation matrix, $\left[{\mathrm{G}}_{\mathrm{k}\mathrm{g}}\right]$ hat relates the components of structural grid point deflections $\left\{{\mathrm{u}}_{\mathrm{g}}\right\}$ to the deflections of the aerodynamic grid points $\left\{{\mathrm{u}}_{\mathrm{k}}\right\}$,

$$\left\{{\mathrm{u}}_{\mathrm{k}}\right\}=\left[{\mathrm{G}}_{\mathrm{k}\mathrm{g}}\right]\left\{{\mathrm{u}}_{\mathrm{g}}\right\}$$

(12)

In the case of static aeroelasticity, the downwash relation presented as Equation 8 becomes:

$$\left\{{\mathrm{w}}_{\mathrm{j}}\right\}=\left[{\mathrm{D}}_{\mathrm{j}\mathrm{k}}\right]\left\{{\mathrm{u}}_{\mathrm{k}}\right\}+\left[{\mathrm{D}}_{\mathrm{j}\mathrm{x}}\right]\left\{{\mathrm{u}}_{\mathrm{x}}\right\}+\left\{{\mathrm{w}}_{\mathrm{j}}^{\mathrm{g}}\right\}$$

(13)

where $\left\{{\mathrm{w}}_{\mathrm{j}}\right\}$ is the vector of aerodynamic degrees of freedom, also called angle of attack, $\left\{{\mathrm{u}}_{\mathrm{k}}\right\}$ is the vector of aerodynamic displacements or deformations, $\left\{{\mathrm{u}}_{\mathrm{x}}\right\}$ is the vector of “extra aerodynamic points” used to describe control surface deflections and overall rigid body motions. $\left\{{\mathrm{w}}_{\mathrm{j}}^{\mathrm{g}}\right\}$ is an initial static aerodynamic downwash that includes the static incidence distribution that may arise from an initial angle of attack, camber, or washout (twist). $\left[{\mathrm{D}}_{\mathrm{j}\mathrm{k}}\right]$ denotes the substantial derivative matrix for the aerodynamic displacements. This is basically the ${\mathrm{D}}_{\mathrm{j}\mathrm{k}}^1$ term in Equation 11 where the ${\mathrm{D}}_{\mathrm{j}\mathrm{k}}^2$ term is not used for the extra aerodynamic points [25].

The theoretical aerodynamic pressures could be derived from:

$$\left\{{\mathrm{f}}_{\mathrm{j}}\right\}=\overline{\mathrm{q}}{\left[{\mathrm{A}}_{\mathrm{j}\mathrm{j}}\right]}^{-1}\left\{{\mathrm{w}}_{\mathrm{j}}\right\}$$

(14)

where $\left\{{\mathrm{f}}_{\mathrm{j}}\right\}$ is the vector of pressure on lifting elements. And finally the aerodynamic forces could be presented by

$$\left\{{\mathrm{P}}_{\mathrm{k}}\right\}=\overline{\mathrm{q}}\left[{\mathrm{W}}_{\mathrm{k}\mathrm{k}}\right]\left[{\mathrm{S}}_{\mathrm{k}\mathrm{j}}\right]{\left[{\mathrm{A}}_{\mathrm{j}\mathrm{j}}\right]}^{-1}\left\{{\mathrm{w}}_{\mathrm{j}}\right\}+\mathrm{q}\left[{\mathrm{S}}_{\mathrm{k}\mathrm{j}}\right]\left\{\frac{{\mathrm{f}}_{\mathrm{j}}^{\mathrm{e}}}{\overline{\mathrm{q}}}\right\}$$

(15)

where $\left[{\mathrm{W}}_{\mathrm{k}\mathrm{k}}\right]$ is a matrix of empirical correction factors to adjust each theoretical aerodynamic box lift and moment to agree with experimental data for incidence changes [25,38]. $\left\{\frac{{\mathrm{f}}_{\mathrm{j}}^{\mathrm{e}}}{\overline{\mathrm{q}}}\right\}$ is the vector of experimental pressure coefficients at some reference incidence (e.g., zero angle of attack) for each aerodynamic element.

3. Results

In this section, the results of the rigid airplane will be initially presented, followed by the elastic wing analysis. These analyses will then be integrated to provide comprehensive aerodynamic characteristics of the full aircraft.

3.1. Rigid Aircraft Aerodynamic Analysis

The aerodynamic characteristics of the twin-engine small aircraft, utilizing the geometry outlined in Figure 4 (generated by MAPLA’s geometry module) and detailed in Table 3, have been examined across various flight conditions. This comprehensive analysis aims to better understand the behavior of the investigated airplane, particularly in relation to propeller effects.

The longitudinal aerodynamic characteristics of the twin-engine propeller-driven aircraft, utilizing the geometry outlined in Table 3, are analyzed in Figure 5 across various flight conditions, including power-off, cruise, and different thrust coefficients such as medium (0.1) and higher (0.3) values typical of take-off conditions. The results highlight significant changes attributed to propeller effects in all longitudinal aerodynamic characteristics compared to the power-off condition. Notably, while lift and drag characteristics exhibit relatively minor alterations with increased values, the pitching moment coefficient demonstrates substantial variations due to heightened propeller effects, as indicated by the thrust coefficient.

Figure 6 illustrates the lateral-directional static results of the twin-engine propeller-driven small aircraft in different flight conditions. It is evident that, similar to the longitudinal characteristics, propeller effects intensify the values across all angles of attack.

Figure 7 shows the longitudinal dynamic characteristics results of the twin-engine propeller-driven small aircraft. It is notable that, similar to the trends observed in the static characteristics, the influence of propeller effects amplifies the values across all angles of attack for the longitudinal dynamic characteristics. This indicates a consistent pattern where the propulsion system significantly impacts the aircraft's dynamic behavior, underscoring the importance of considering propeller effects in aerodynamic analyses for accurate performance predictions.

Finally, Figure 8 shows the lateral-directional dynamic characteristics of the of the twin-engine propeller-driven small aircraft. Remarkably, akin to the observations made in the static characteristics, the influence of propeller effects is discernible, intensifying the values across all angles of attack. However, it is noteworthy that the results for Cl_p and Cn_r were more significantly impacted by the propeller effects compared to those of Cl_r and Cn_p. This consistency highlights the significant impact of propeller dynamics on the aircraft's lateral-directional behavior, underscoring the need to account for such effects for comprehensive aerodynamic analysis and design optimization.

3.2. Aircraft Aerodynamic Analysis with Elastic Wing Integration

As previously mentioned, the aeroelastic analysis presented here adopts the NASTRAN approach, primarily focusing on static aeroelastic problems. This approach delves into the intricate interaction between aerodynamic and structural forces acting upon a flexible vehicle, leading to the redistribution of aerodynamic loading relative to airspeed. Specifically, the static aeroelastic study addresses the static aeroelastic problem under aircraft trim conditions [25]. The wing modeling process was conducted using PATRAN, adhering to the methodology proposed by [39]. Figure 9 illustrates the resulting geometry of the wing considered for analysis in the subsequent steps.

In the subsequent phase, the Aero-Structure Coupling module of NASTRAN is employed to generate the coupled aerodynamic and structural models for aeroelastic analysis. The resultant Aero-Structural Coupling model is depicted in Figure 10. Following the creation of the Aero-Structural Coupling model, the Aeroelastic Analysis module is utilized, with the solution type set to Static Aeroelasticity. Specifically, the Flexible Trim method from NASTRAN is employed for this analysis.

Figure 11a presents the pressure distribution across the wing surface in Newton per square meter, emphasizing aerodynamic behavior while assuming rigid components within the wing model considering cruise flight condition at a speed of 90 m/s. Notably, a linear distribution of pressure along the wingspan is evident, accompanied by chordwise nonlinear behavior. Higher pressure values are observed on the front side, contrasting with lower values on the aft side of the wing model. This depiction provides valuable insights into the interaction between aerodynamic forces and the wing structure under conditions of rigidity. In contrast, Figure 11b explores the same analysis while accounting for the presence of elastic components within the wing model. As illustrated, the pressure distribution is influenced by the behavior of elasticity. This representation offers a nuanced understanding of how aerodynamic forces are impacted by the incorporation of elasticity into the wing structure. By comparing these two scenarios, we gain valuable insights into the effects of structural flexibility on aerodynamic behavior and its consequential implications for overall aircraft performance. The resulting aerodynamic characteristics for both rigid and elastic configurations of the total aircraft are consolidated in Table 1, covering the linear range of angle of attack.

Figure 12a illustrates the moment distribution across the wing surface in Newton per meter, assuming rigid components within the wing model considering cruise flight condition at a speed of 90 m/s. Notably, it reveals a linear distribution of moment along the wingspan, accompanied by chordwise nonlinear behavior. Higher moment values are concentrated on the front side, contrasting with lower values on the aft side of the wing model. This depiction provides valuable insights into the interaction between moment changes and the wing structure under conditions of rigidity. In contrast, Figure 12b delves into the same analysis, this time considering the moment distribution of elastic components within the wing model. As depicted, the moment distribution is notably influenced by the behavior of elasticity. By comparing these two scenarios, we gain valuable insights into the effects of structural flexibility on moment behavior and its consequential implications for overall aircraft performance.

Finally, Figure 13a illustrates the force distribution across the wing surface in Newton, emphasizing aerodynamic behavior under the assumption of rigid components within the wing model at a speed of 90 m/s. As with the previous characteristics, a linear distribution of force along the wingspan is evident, accompanied by chordwise nonlinear behavior. Higher force values are concentrated on the front side, while lower values are observed on the aft side of the wing model. This representation offers valuable insights into the interaction between aerodynamic forces and the wing structure under conditions of rigidity. In contrast, Figure 13b delves into the same analysis, this time considering the force distribution of elastic components within the wing model. As depicted, the force distribution is notably influenced by the behavior of elasticity. This representation offers a nuanced understanding of how aerodynamic forces are impacted by the incorporation of elasticity into the wing structure. By comparing these two scenarios, we gain valuable insights into the effects of structural flexibility on force behavior throughout the wing model surface and its consequential implications for overall aircraft performance. The resulting aerodynamic characteristics for both rigid and elastic configurations of the total aircraft are consolidated in Table 1, encompassing the linear range of angle of attack.

In the subsequent step, leveraging MAPLA’s modular capability, the aerodynamic analysis for the wing component was updated using the attained aerodynamic results from NASTRAN. Subsequently, the remaining analysis was conducted by running MAPLA based on the new results for the wing, while considering propeller effects and power contributions. The summary of the results for both rigid and elastic aircraft configurations in different power settings is presented in Table 4, focusing on the linear angles of attack. Furthermore, the presentation of the average deviation from the rigid aircraft configuration offers a more comprehensive understanding, enabling a clearer distinction of the contribution and significance of aeroelasticity in the analysis. Regarding these deviations, it is evident that the aircraft would be notably influenced by elasticity, particularly for longitudinal characteristics, with an increase of nearly 20%. Conversely, for lateral-directional characteristics, this effect appears to be relatively lower, with an estimated impact of approximately 4%. The results achieved in this study are also consistent with the other studies involved with the aeroelasticity aircraft analysis [40].

5. Conclusions

In conclusion, this paper explored the critical influence of propeller effects, power contribution, and elasticity on small propeller-driven aircraft aerodynamics. Through meticulous analysis, a twin-engine propeller-driven aircraft was examined, shedding light on their combined impact. Validation through comparative analysis and NASTRAN software utilization reinforced our findings. Leveraging MAPLA’s capability, the aerodynamic analysis was updated for the wing component, considering propeller effects and power contributions. Results underscore the significance of these factors at various angles of attack. Additionally, deviations from the rigid aircraft configuration highlighted a substantial impact of aeroelasticity, particularly in longitudinal characteristics, with approximately 20% increase, and a lower impact of less than 5% in lateral-directional characteristics. This study advances small propeller-driven aircraft design considerations, offering insights for future research in UAM. By deepening the understanding of aerodynamic complexities, future works can contribute to more efficient and optimized aircraft designs in the evolving aviation landscape.