Hostname: page-component-586b7cd67f-dlnhk Total loading time: 0 Render date: 2024-11-27T20:40:32.947Z Has data issue: false hasContentIssue false

FLIGHT LIMITATIONS IMPOSED ON SINGLE ROTOR AND COAXIAL HELICOPTERS BY THE LIFT EQUATION

Published online by Cambridge University Press:  17 July 2023

B. MALDON
Affiliation:
School of Information and Physical Sciences, University of Newcastle, Callaghan, NSW 2308, Australia; e-mail: [email protected]
MICHAEL H. MEYLAN*
Affiliation:
School of Information and Physical Sciences, University of Newcastle, Callaghan, NSW 2308, Australia; e-mail: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

To compute the maximum speed threshold for helicopters, we model the lift produced by the rotor blades. Using this model, we derive limits for each method traditionally used to alleviate dissymmetry of lift. Additionally, we find the minimum rotor angular velocity required to produce a prescribed lift at a given forward velocity. We derive conditions on the coefficient of lift for helicopter airfoils that maintain altitude. Further considerations are also made with regard to the properties of the air and its effect on helicopter dynamics.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2023. Published by Cambridge University Press on behalf of Australian Mathematical Publishing Association Inc.

1 Introduction

Helicopters remain a vital component of aviation technology, first realised in Igor Sikorsky’s VS300 model, which was the first successful helicopter to feature a main rotor and tail rotor [Reference Fay1]. La Mantia and Dabnichki calculated optimal airfoil thicknesses based on the standard NACA airfoil shapes [Reference La Mantia and Dabnichki4]. Kulchenko et al. applied mathematical models for single-rotor helicopters for autopilot technology in 2014 [Reference Kulchenko, Kostukov, Verevkina and Chufistov3]. Finally, Romani and Casalino elucidated the influence of rotor-induced vortices on the overall production of lift on the helicopter blades [Reference Romani and Casalino10]. The theory of helicopter dynamics encompasses the study of the complex aerodynamic and mechanical principles that govern the flight of helicopters. Helicopters are unique aircraft that rely on the interaction between components, including the main rotor, tail rotor, fuselage and various control systems, to achieve controlled vertical takeoff, landing and manoeuvrability in the air. The aerodynamics of helicopter rotors involve lift, drag, induced flow and rotor wake, which are influenced by blade geometry, airfoil design and rotor speed. Understanding the dynamic behaviour of helicopters requires knowledge of key concepts such as cyclic and collective pitch control, autorotation, ground effect and gyroscopic precession. Numerous scientific studies and engineering advancements have contributed to the development and refinement of helicopter dynamics, many of which can be found in [Reference Johnson2, Reference Leishman5].

Perdomo and Wei studied blade flapping in their paper [Reference Perdomo and Wei9], providing a time-dependent truncated Fourier series to calculate the flapping angle for an airfoil across the rotor disk. Majhi and Ganguli studied dynamic stall and the efficacy of assuming small flapping angles in mathematically modelling helicopters [Reference Majhi and Ganguli7].

Applications for airfoils extend beyond aviation technology for helicopters. Liu et al. recently used a dynamic stall model to study wind turbines for renewable energy applications [Reference Liu, Liang, Li, Godbole and Lu6]. Rozhdestvinsky and Ryzhov completed a survey for flapping-wing propulsors as they apply to aerohydrodynamics in 2003 [Reference Rozhdestvinsky and Ryzhov11].

We define a helicopter’s fundamental components to capture the limiting factors arising from flying helicopters at high forward velocities. Primarily, this paper will study the lift equation as a model for the thrust collectively produced by the helicopter blades.

1.1 Helicopter rotor blades

During operation, helicopter rotor blades spin about the rotor shaft at a typically constant speed $\Omega _{R}$ (usually measured in RPM). We will assume all rotor blades rotate anticlockwise, with the exception of the second rotor disk of a coaxial helicopter or when otherwise noted. Each blade, therefore, has an angular velocity $\Omega $ of

$$ \begin{align*} \Omega = \frac{\pi\Omega_{R}}{30}. \end{align*} $$

Given a helicopter rotor disk of n blades, the azimuth angle $\theta $ for blade k of n at time t is given by

$$ \begin{align*} \theta(t) = \Omega t + \frac{2\pi(k - 1)}{n}, \end{align*} $$

so that at time $t = 0$ , the first helicopter blade in the rotor disk is in the advancing position. Using the azimuth angle, the tip $B_{\text {tip}}$ of blade k of n at time t is therefore given by

$$ \begin{align*} B_{\text{tip}} = R [ \cos(\theta), \sin(\theta), 0 ], \end{align*} $$

where R is the length of the helicopter blade (and, therefore, the radius of the rotor disk).

1.2 Dissymmetry of lift

Under static hover conditions, the lift is generated uniformly across the rotor disk as the blades pass through the air. When helicopter blades are subject to a relative wind (usually induced by the forward velocity of the helicopter), more lift is generated on the advancing side of the rotor disk and less lift is generated on the retreating side. This dissymmetry of lift is detailed in Figure 1 and shows a fundamental problem encountered by single-rotor helicopters. Coaxial helicopters amend this issue with their unique counter-rotating disks.

Figure 1 Diagram showing the direction of rotation with the relative wind.

1.3 Helicopter rotor blade airfoils

To combat the dissymmetry of lift, helicopter rotor blades may rotate about three axes, as shown by Figure 2. The lead-lagging angle $\xi $ , flapping angle $\beta $ and feathering angle $\zeta $ are all used to decrease the lift on the advancing blades and increase lift on the retreating blades.

Figure 2 Diagram of the three angles of rotation for a helicopter rotor blade.

The primary limitations of these measures are found in retreating blade stall and supersonic blade limitations. Figure 3 shows the effect of the retreating blade stall coefficient of lift for a NACA0012 airfoil. Once the angle of attack $\alpha $ approaches its critical value $\alpha _{c}$ , the lift increase plateaus and leads to decreased lift if the angle of attack exceeds the critical value.

Figure 3 Plot of the coefficient of lift against the angle of attack for the NACA0012 airfoil for Reynolds numbers from $5\times 10^5$ to $10^6$ [8]. The x-axis is in degrees.

In the problem of supersonic blade limitations, helicopter rotor blade tips cannot exceed the speed of sound without encountering excessive drag and air resistance liable to break the blades. Mathematically, the angular velocity $\Omega $ of the blades must satisfy

$$ \begin{align*} \Omega R + u < \nu, \end{align*} $$

where $\nu $ is the speed of sound. Rearranging for $\Omega $ , we have

(1.1) $$ \begin{align} \Omega = \eta \frac{\nu - u}{R}, \end{align} $$

where $\eta < 1$ is a tolerance factor. Unless otherwise stated, we will assume $\eta = 1$ for theoretical calculations to obtain the maximum possible lift.

2 Mathematical model

In this section, we model the lift produced by helicopter blade k of n with the standard lift equation

(2.1) $$ \begin{align} \ell_{\text{std}} = \tfrac{1}{2}\rho c v^{2} C_{L}, \end{align} $$

where $\ell _{\text {std}}$ is the lift per unit length, $\rho $ is the air density, c is the chord length, v is the velocity of the blade and $C_{L}$ is the coefficient of lift. In this study, we will consider the lift generated by one blade with azimuth angle $\theta = \Omega t$ . We are assuming uniform chord along blade, straight blades and uniform angle of attack.

2.1 Assumptions and parameter values

We will only consider the lift generated on the disk with polar coordinates. For time-dependent effects during operation, we will use the azimuth angle expression

$$ \begin{align*}\theta(t) = \Omega t + \frac{2\pi(k - 1)}{n}.\end{align*} $$

The coefficient of lift $C_{L}$ is regarded to be a function of angle of attack $\alpha $ only. We will assume $C_{L}$ is either constant or linear with respect to $\alpha $ , based on Figure 3.

Throughout this paper, the rotor blades have a length $3$ m and a chord length of $180$ mm (or $0.18$ m). The default mass of the helicopter is $1000$ kg.

For reference, we will use the NACA0012 airfoil [8] for the coefficient of lift computations. Other airfoils will also be compatible with these expressions, given the lift coefficients as functions of the angle of attack and Reynolds number.

2.2 Standard lift equation

Recall the standard lift equation is given by

$$ \begin{align*} \ell_{\text{std}} = \tfrac{1}{2}\rho c v^{2} C_{L}. \end{align*} $$

To study the effect of a relative wind, we model the linear velocity with the equation

$$ \begin{align*} \nu = \Omega r + u\cos(\theta), \end{align*} $$

where u is the relative wind speed and $\theta $ is the blade azimuth angle. The lift becomes a function of the position on the blade $r \in [0, R]$ and azimuth angle $\theta \in [0, 2\pi )$ .

Figure 4 shows the lift distribution for a helicopter rotor blade in static hover using equation (2.1). In the absence of relative wind, the lift is equal across the rotor disk. The lift is maximized at the tip of the rotor disk due to the maximal linear velocity at the blade tips.

Figure 4 Lift distribution of the rotor disk in standard hover.

Figure 5 shows the lift distribution for a helicopter rotor blade under the relative winds $u = 50\,\text {m}\text {s}^{-1}$ and $u = 120\,\text {m}\text {s}^{-1}$ , using equation (2.1). Without adjusting the angle of attack, the relative wind decreases the lift produced on the retreating side of the rotor disk.

Figure 5 Lift distribution of the rotor disk given a relative wind $u = 50\,\text {m}\text {s}^{-1}$ (left) and $u = 120\,\text {m}\text {s}^{-1}$ (right).

3 Results

3.1 Revised lift equation

In this section, we compute the lift generated by a helicopter rotor blade with the modified equation

$$ \begin{align*} L_{\text{mod}}(u, \theta) = \int_{0}^{R} \frac{\rho c C_{L}}{2}(\Omega r + u\cos(\theta))|\Omega r + u\cos(\theta)| \,dr. \end{align*} $$

Similar expressions can be found in [Reference Johnson2, Reference Leishman5]. Expressing the velocity term as $\nu |\nu |$ rather than $\nu ^{2}$ allows the term to become negative, as experienced by retreating blades at high velocity. Upon explicit calculation,

$$ \begin{align*} L_{\text{mod}}(u, \theta) &= \frac{\rho c C_{L}}{6\Omega}[(\Omega R + u\cos(\theta))^{2} |\Omega R + u\cos(\theta)| - u^{3}\cos^{2}(\theta)|\!\cos(\theta)|]. \end{align*} $$

3.1.1 Retreating blade special case

For retreating blades ( $\theta = 180^{\circ }$ ), the lift generated simplifies to

$$ \begin{align*} L_{\text{mod}, \, R}(u) = \frac{\rho c C_{L}}{6\Omega}[(\Omega R - u)^{2}|\Omega R - u| - u^{3}]. \end{align*} $$

Given $\Omega =(\nu - u)/R$ to ensure maximum RPM, the lift becomes

$$ \begin{align*} L_{\text{mod}, \, R}(u) = \frac{R\rho c C_{L}}{6(\nu - u)}[(\nu - 2u)^{2}|\nu - 2u| -u^{3}]. \end{align*} $$

To find u such that $L_{R}(u) = 0$ , we solve

$$ \begin{align*} [(\nu - 2u)^{2}|\nu - 2u| - u^{3}] = 0. \end{align*} $$

Given $u < \nu /2$ ,

$$ \begin{align*} u = \frac{\nu}{3}. \end{align*} $$

Given $\nu \kern1.3pt{=}\kern1.3pt 343\,\text {m}\text {s}^{-1}$ , we see that traditional helicopters cannot exceed ${\nu /3 \kern1.3pt{\approx}\kern1.3pt 114.3\,\text {m}\text {s}^{-1}}$ . This limit is slightly higher than the known world speed record of $111.35\,\text {m}\text {s}^{-1}$ by the Westland Lynx set in 1986 [Reference Fay1].

3.1.2 Negative lift region

If $u> \nu /3$ , a region of the rotor disk about the tip of the retreating blade will produce negative lift when integrated. To find this region, we solve $L(u, \theta _{c}) = 0$ for $\theta _{c}$ given $u> \nu /3$ . That is, we solve

$$ \begin{align*} \frac{\rho c C_{L}}{6\Omega}[(\Omega R + u\cos(\theta_{c}))^{2}|\Omega R + u\cos(\theta_{c})| - u^{3}\cos^{2}(\theta_{c})|\!\cos(\theta_{c})|] = 0. \end{align*} $$

Letting $\epsilon = \cos (\theta _{c})$ and using $\Omega = (\nu - u)/R$ ,

$$ \begin{align*} (\nu - u(1 - \epsilon))^{2}|\nu - u(1 - \epsilon)| - u^{3}\epsilon^{2}|\epsilon| = 0. \end{align*} $$

Since $\epsilon < 0$ , by expansion and factoring,

$$ \begin{align*} \bigg(\frac{\nu - u}{u}\bigg)(2\epsilon + \nu)\bigg(\epsilon - 2\frac{\nu - u}{u}\bigg) = 0. \end{align*} $$

Therefore, we have $\epsilon = (\nu - u)/2u$ . Upon substitution,

$$ \begin{align*} \theta_{c} = \arccos\bigg(\frac{\nu - u}{2u}\bigg). \end{align*} $$

Therefore, the helicopter blades will produce negative lift for azimuth angles, satisfying

$$ \begin{align*} \theta \in \bigg[\pi - \arccos\bigg(\frac{\nu - u}{2u}\bigg), \, \pi + \arccos\bigg(\frac{\nu - u}{2u}\bigg)\bigg]. \end{align*} $$

3.2 Modified lift equation

To recapture the loss of lift on the retreating rotor blades, we modify the lift equation to be

(3.1) $$ \begin{align} \ell_{\text{mod}} = \tfrac{1}{2}\rho c v|v| C_{L}, \end{align} $$

replacing $v^{2}$ with $v|v|$ . This allows the velocity of the blade to take negative values, which will occur on the retreating side for sufficiently high relative wind speeds. The lift distribution in static hover matches Figure 4 as $v = \Omega r$ . Given a relative wind u, we set $v = \Omega r + u\cos (\theta )$ as before. Figure 6 shows that the relative wind greatly affects the lift produced near the rotor shaft on the retreating side.

Figure 6 Modified lift distribution of the rotor disk given a relative wind $u = 50\,\text {m}\text {s}^{-1}$ (left) and $u = 120\,\text {m}\text {s}^{-1}$ (right).

Figure 7 Plot of the lift generated over time for several forward velocities.

3.2.1 Modified integrated lift equation

Using the same techniques as the standard lift equation, we may calculate the total lift produced by a helicopter rotor blade by integrating the lift over the blade. That is, we compute

$$ \begin{align*} L(\alpha, u, \theta) = \int_{0}^{R} \frac{1}{2}\rho c (\Omega r + u\cos(\theta))|\Omega r + u\cos(\theta)|C_{L}(\alpha) \,dr. \end{align*} $$

Upon computation,

$$ \begin{align*} L(\alpha, u, \theta) = \frac{\rho c C_{L}(\alpha)}{6\Omega}[(\Omega R + u\cos(\theta))^{2}|\Omega R + u\cos(\theta)| - u^{3}\cos^{2}(\theta)|\!\cos(\theta)|]. \end{align*} $$

With $\Omega = (\nu - u)/R$ , the lift becomes

$$ \begin{align*} L(\alpha, u, \theta) = \frac{R\rho c C_{L}(\alpha)}{6(\nu - u)}[(\nu - u + u\cos(\theta))^{2}|\nu - u + u\cos(\theta)| - u^{3}\cos^{2}(\theta)|\!\cos(\theta)|]. \end{align*} $$

Figure 7 plots the lift generated by the blade over time using the usual azimuth angle $\theta (t) = \Omega t$ . The lift oscillates as the blade transitions between advancing and retreating. As the forward velocity increases, the oscillation amplitude increases and the centre decreases. As the velocity may take negative values with the expression $\nu |\nu |$ , we see the special case $u = 120\,\text {m}\text {s}^{-1}$ produces negative lift at the retreating side.

To ensure $L = mg/n$ for vertical hover, the coefficient of lift must satisfy

$$ \begin{align*} C_{L}(\alpha) = \frac{6\Omega m g}{\rho c n}[(\Omega R + u\cos(\theta))^{2}|\Omega R + u\cos(\theta)| - u^{3}\cos^{2}(\theta)|\!\cos(\theta)|]^{-1}. \end{align*} $$

With $\Omega = (\nu - u)/R$ , the coefficient of lift becomes

$$ \begin{align*} C_{L}(\alpha) &= \frac{6\eta(\nu - u) m g}{\rho R c n}\\ &\quad\times [(\eta(\nu - u) + u\cos(\theta))^{2}|\eta(\nu - u) + u\cos(\theta)| - u^{3}\cos^{2}(\theta)|\!\cos(\theta)|]^{-1}. \end{align*} $$

Figure 8 is a plot of the coefficient of lift over time given the standard azimuth angle $\theta (t) = \Omega t$ . For static hover and the relatively low forward velocity $50\,\text {m}\text {s}^{-1}$ , we see the coefficient of lift take reasonable values between $0$ and $1$ . For the forward velocity $120\,\text {m}\text {s}^{-1}$ , we see that the negative lift produced on the retreating side leads to discontinuities the coefficient of lift cannot reconcile. The angle of attack is, therefore, unable to balance the lift on the retreating side for forward velocities $u> \nu /3$ .

Figure 8 Plot of the lift coefficient $C_{L}$ required for vertical hover against time.

3.3 Average modified lift

To better understand the lift performance of the helicopter blade over the entire rotor disk, we also compute the average lift produced over time $\lambda $ with the expression

$$ \begin{align*} L_{\text{mod}, \lambda} = \frac{1}{\lambda}\int_{0}^{\lambda} \int_{0}^{R} \frac{1}{2} \rho c (\Omega r + u\cos(\Omega t))|\Omega r + u\cos(\Omega t) |C_{L}(\alpha) \,dr \,dt. \end{align*} $$

We compute this integral numerically with a standard midpoint Riemann sum with 300 evenly spaced intervals in Figure 9. As the forward velocity increases, the negative lift produced by the retreating side lowers the average lift across the disk. A higher forward velocity results in a more significant difference.

Figure 9 Plot of the average lift $L_{\text {mod}, T}$ over time interval $\lambda $ .

3.4 Average lift across rotor disk

Since the helicopter rotor speed $\Omega $ is dependent on u, the time interval required to consider one revolution is also dependent on u. Therefore, for the standard lift equation, we compute the average value over one particular revolution with the alternative integral

$$ \begin{align*} L_{\text{std}, T}(u) = \frac{1}{2\pi}\int_{0}^{2\pi} \int_{0}^{R} \frac{1}{2}\rho c \bigg(\frac{\nu - u}{R}r + u\cos(\theta)\bigg)^{2}C_{L} \,dr\, d\theta. \end{align*} $$

Upon simplification,

$$ \begin{align*} L_{\text{std}, T}(u) = \frac{\rho c R C_{L}(\alpha)}{12}[5u^2 -4\nu u + 2\nu^2]. \end{align*} $$

For the modified lift, we compute

$$ \begin{align*} L_{\text{mod}, T}(u) = \frac{1}{2\pi}\int_{0}^{2\pi} \int_{0}^{R} \frac{1}{2}\rho c \bigg(\frac{\nu - u}{R}r + u\cos(\theta)\bigg)\bigg| \frac{\nu - u}{R}r + u\cos(\theta) \bigg|C_{L} \,dr\, d\theta. \end{align*} $$

Finally, to reconcile the negative lift from the modified velocity model $\nu |\nu |$ , we also compute the integral

$$ \begin{align*} L_{\max, T} = \frac{1}{2\pi}\int_{0}^{2\pi} \int_{0}^{R} \max\bigg(\frac{1}{2}\rho c \bigg(\frac{\nu - u}{R}r + u\cos(\theta)\bigg) \bigg| \frac{\nu - u}{R}r + u\cos(\theta) \bigg|C_{L}, 0\bigg) \,dr\, d\theta. \end{align*} $$

To better understand the lift differences, we scale the result by the lift obtained in static hover. From Figure 10, we find the models produce similar average lifts for low forward velocity as the effect of dissymmetry of lift is at its lowest influence.

Figure 10 Plot of the normalized average lift over forward velocity for each model.

As the forward velocity increases, the lift produced by the blade decreases as the rotor speed $\Omega $ must decrease. Once the negative lift is produced by the retreating blade at $u = \nu /3$ , we find the lift from the standard and modified model diverge. The standard model increases lift as the forward velocity remains nonnegative. For the modified model, the overall lift decreases as the negative lift on the retreating blades counteracts the positive lift produced by the advancing sides. If the angle of attack can be set to $0^{\circ }$ , the negative effect of the retreating blade diminishes and the lost lift begins to recover.

Since the retreating side of the rotor disk cannot produce any positive lift for velocities above $114.3\,\text {ms}^{-1}$ , we find that only coaxial helicopters are able to rectify the issue. We also note the lift produced in these velocities varies between $50\%\ \mathrm{and}\ 75\%$ of the lift generated in static hover. This is primarily a consequence of the decreased rotor speed $\Omega $ as u increases.

3.5 Minimum rotor speed requirements

To obtain the maximum lift attainable from helicopter blades, we prescribe $\Omega $ with equation (1.1) to the lift equation (3.1). To find the minimum requirements for $\Omega $ to hold a given lift L, we estimate (to enable us to give an approximate analytic formula) the average lift by

$$ \begin{align*} L_{\text{avg}}(u, \Omega) = \frac{1}{4}\bigg[L(u, 0, \Omega) + 2L\bigg(u, \frac{\pi}{2}, \Omega \bigg) + L(u, \pi, \Omega)\bigg], \end{align*} $$

which simplifies to

$$ \begin{align*} L_{\text{avg}}(u, \Omega) = \frac{\rho c}{24\Omega}[(\Omega R + u)^{3} + (\Omega R - u)^{2}|\Omega R - u| + 2\Omega^{3}R^{3} - 2u^{3}]. \end{align*} $$

Solving for $\Omega $ ,

(3.2) $$ \begin{align} \Omega(u, L) = \left\{ \begin{aligned} &2^{-2/3}\dfrac{2^{7/3}L - 2^{1/3}(Rc\rho u^{2}) + ( (Rc\rho u^{2})^{3/2} + K)^{2/3}}{R\sqrt{Rc\rho}((Rc\rho u^{2})^{3/2} + K )^{1/3}}, & u < \sqrt{\frac{3L}{c\rho R}} , \\[.1cm] & -\dfrac{3u}{2R} + \dfrac{\sqrt{3c\rho R(16L + 3c\rho Ru^{2})}}{2c\rho R^{2}}, & u \geq \sqrt{\frac{3L}{c\rho R}}, \end{aligned} \right. \end{align} $$

where $K = \sqrt {-128L^{3} + 96(Rc\rho u^{2})L^{2} - 24(Rc\rho u^{2})^{2}L + 3(Rc\rho u^{2})^{3}}$ .

From Figure 11, we see that the minimum threshold for $\Omega $ decreases as forward velocity u increases. While this will allow the rotor disk to spin at a lower RPM, the maximum threshold imposed by equation (1.1) decreases at a faster rate. Consequently, the interval between the maximum and minimum values for $\Omega $ decreases in length until they coincide at some $u = u_{M}(L)$ dependent on the lift L required by the minimum threshold.

Figure 11 Plot of minimum $\Omega $ for various weights with the theoretical maximum against forward velocity.

By solving $\Omega (u_{M}, L) = (\nu - u_{M})/R$ , we are able to find the maximum allowable speed limit as imposed by the minimum lift L that the rotor disk(s) must generate. Assuming $u_{M}> 3L/c\rho R$ , we find the maximum speed limit is given by

(3.3) $$ \begin{align} u_{M}(L) = \frac{\nu}{4} + \frac{\sqrt{3Rc\rho(3Rc\rho\nu^{2} - 32L)}}{4Rc\rho}. \end{align} $$

In Figure 12, we plot the maximum speed limit as a function of lift required for a variety of helicopter blade lengths. If the forward velocity of a helicopter should exceed this limit (if $u> u_{M}$ ), the rotor angular velocity $\Omega $ must decrease, and the lift produced will not be sufficient to maintain altitude. To compensate for this, other methods for producing lift will be necessary, such as wings or a modified fuselage shape.

Figure 12 Plot of maximum speed $u_{M}$ against lift L for several blade lengths.

We note in equation (3.3) that $u_{M}$ attains its minimum at $L = 3Rc\rho \nu ^{2}/32$ , where $u_{M} = \nu /4$ .

3.6 Lift on the advancing side of the rotor disk

To study the effect of the relative wind on the advancing side of the rotor disk, we calculate the integral

$$ \begin{align*} L_{\text{adv}}(u) = \int_{-\pi/2}^{\pi/2} \int_{0}^{R} \frac{1}{2\pi}\rho c (\Omega r + u\cos(\theta))|\Omega r + u\cos(\theta)| \,dr \,d\theta. \end{align*} $$

Since $\Omega r + u\cos (\theta )> 0$ for advancing rotor blades and the integral is even, we may simplify the expression to

$$ \begin{align*} L_{\text{adv}}(u) = \frac{2}{\pi} \int_{0}^{\pi/2} \int_{0}^{R} \frac{1}{2}\rho c (\Omega r + u\cos(\theta))^{2} \,dr \,d\theta. \end{align*} $$

Upon simplification,

$$ \begin{align*} L_{\text{adv}}(u) = \frac{R\rho c}{12\pi}[3\pi u^{2} + 12\Omega R u + 2\pi(\Omega R)^{2}]. \end{align*} $$

In Figure 13, we plot $L_{\text {adv}}$ against u for the maximum theoretical $\Omega $ prescribed by equation (1.1) and for the minimum values for $\Omega $ required to maintain altitude for $1000$ kg, $2000$ kg and $3000$ kg helicopters. As the forward velocity increases, it reaches the maximum allowable limit $u_{M}$ , at which point the lift generated coincides with the lift obtained by the maximal rotor angular velocity $\Omega $ .

Figure 13 Plot of the average lift over the advancing half of the helicopter rotor disk against forward velocity.

To elucidate the effect of the angle of attack on the advancing side of the blade, we multiply the lift by the coefficient of lift using the maximal values for $\Omega $ . In Figure 14, we see lower angles of attack predictably produce less lift up to the stall threshold $\alpha _{c} \approx 10^{\circ }$ . When the angle of attack exceeds this limit, the lift begins to decrease, as expected with retreating blade stall.

Figure 14 Plot of the average lift over the advancing half of the helicopter rotor disk against forward velocity for several angles of attack.

To additionally incorporate the effect of drag, we alter the integral on the advancing blades with the proportion $(C_L(\alpha ) - C_D(\alpha ))$ . Given that drag increases sharply with the angle of attack, we expect an optimal angle of attack that allows the greatest lift coefficient under the drag coefficient, as shown in Figure 15.

Figure 15 Plot of the lift coefficient and the difference between lift and drag against the angle of attack.

Figure 16 shows the lift on the advancing side with the coefficient of drag included in the computation. Though the effect of the drag is slight, the increased angle of attack magnifies the gap between the lift at $\alpha = 9^{\circ }$ and $\alpha = 12^{\circ }$ .

Figure 16 Plot of the average lift over the advancing half of the helicopter rotor disk against forward velocity for several angles of attack.

3.7 Lift on the retreating side of the rotor disk

To obtain the lift generated on the retreating side of the rotor disk, we calculate

$$ \begin{align*} L_{\text{ret}}(u) = \frac{1}{\pi}\int_{{\pi}/{2}}^{{3\pi}/{2}} \int_{0}^{R} \frac{1}{2}\rho c (\Omega r + u\cos(\theta))|\Omega r + u\cos(\theta)| \,dr \,d\theta. \end{align*} $$

Given that $\Omega r + u\cos (\theta )$ is not necessarily nonnegative nor nonpositive, there is no simplification. To calculate this integral, we use MAPLE’s built-in numerical integration with a partition of $100$ points. In Figure 17, the retreating lift can be seen to rapidly decrease. Without measures to nullify negative lift, the retreating side continues to build negative lift as forward velocity increases. For coaxial helicopters, Figure 13 shows us that the lift increase from the advancing side is more than capable of balancing the lift loss on the retreating side.

Figure 17 Plot of the average lift over the retreating half of the helicopter rotor disk against forward velocity.

3.8 Angle of attack considerations for coaxial helicopters

In our previous mathematical derivations for the lift across the rotor disk (or particular regions of the rotor disk), we have assumed that the coefficient of lift constantly takes its maximum value $C_{L} = 1$ . In this section, we relax this assumption to derive conditions on the coefficient of lift that will enable the lift to be generated equally across the rotor disk.

For low forward velocities, we are able to apply the simple formula

$$ \begin{align*} C_{L}(\theta) = \frac{L}{L_{\text{std}}(u, \theta)}, \end{align*} $$

where $L_{\text {std}}(u, \theta )$ is the lift generated by the blade at forward velocity u and azimuth angle $\theta $ , given by

$$ \begin{align*} L_{\text{std}}(u, \theta) = \frac{\rho c}{6\Omega}[(\Omega R + u\cos(\theta))^{2}|\Omega R + u\cos(\theta)| - u^{3}\cos^{2}(\theta)|\!\cos(\theta)|]. \end{align*} $$

This expression remains valid until the lift generated on the retreating blade is no longer capable of producing sufficient lift. That is, the first threshold is the forward velocity $u_{1}$ such that

$$ \begin{align*} L_{\text{std}}(u_{1}, \pi) = L. \end{align*} $$

Assuming the rotor angular velocity is at its maximum, as shown by equation (1.1), we may simplify the above equation to show $u_{1}$ satisfies

$$ \begin{align*} u_{1}^{3} - \frac{4}{3}\nu u_{1}^{2} - \frac{2}{3}\bigg(\nu^{2} + \frac{L}{\rho c R}\bigg)u_{1} - \frac{\nu}{9}\bigg(\nu^{2} - \frac{6L}{\rho c R}\bigg) = 0. \end{align*} $$

For velocities $u> u_{1}$ , we set the angle of attack on the retreating blades to maximum and increase the angle of attack on the counterpart advancing blade to compensate for the insufficient lift generation.

3.8.1 Lift compensation for retreating blades

As the forward velocity u increases above the first threshold $u_{1}$ defined above, we must produce excess lift on the advancing blades to compensate for the lack of lift produced on the retreating blades. Suppose each blade must generate L lift. If $u < u_{1}$ , $L_{\text {std}}(u, \pi ) < L$ . Therefore, the excess lift E that the advancing blades must produce is given by

$$ \begin{align*} E = L - L_{\text{std}}(u, \pi). \end{align*} $$

In Figure 18, we see the lift compensation E that the advancing blade must produce for its counterpart retreating blade. The curves in Figure 18 are proportional to the maximum lift produced in hover that the highest weight helicopters may attain. Negative values indicate the retreating blade is producing sufficient lift, and compensation is not necessary. As the forward velocity increases, the critical values $E(u) = 0$ coincident with $u_{1}(L)$ are found. The lift compensation required also increases up to the critical value $\nu /3$ .

Figure 18 Plot of the average lift over the retreating half of the helicopter rotor disk against forward velocity.

Therefore, in the case that $u> u_{1}(L)$ , the coefficient of lift will take the form

$$ \begin{align*} C_{L}(\theta) = \begin{cases} \dfrac{2L - \min\bigg(\dfrac{L}{L_{\text{std}}(u, \theta + \pi)}, 1\bigg)L_{\text{std}}(u, \theta + \pi)}{L_{\text{std}}(u, \theta)}, & \theta \in \bigg[0, \dfrac{\pi}{2}\bigg]\cup\bigg[\dfrac{3\pi}{2}, 2\pi\bigg], \\ \min\bigg(\dfrac{L}{L_{\text{std}}(u, \theta)}, 1\bigg), & \theta \in \bigg[\dfrac{\pi}{2}, \dfrac{3\pi}{2}\bigg]. \end{cases} \end{align*} $$

3.8.2 Speed limit imposed on transitioning blades

Consider the lift produced by the blades as they transition between the advancing and retreating halves of the rotor disk (at azimuth angles $90^{\circ }$ and $270^{\circ }$ ). The lift produced on these blades is equal and unaffected by the forward velocity. That is,

$$ \begin{align*} L_{\text{std}}\bigg(u, \frac{\pi}{2}\bigg) = L_{\text{std}}\bigg(u, \frac{3\pi}{2}\bigg) = \frac{\rho c R^{3} \Omega^{2}}{6}. \end{align*} $$

Since the rotor angular velocity $\Omega $ decreases as forward velocity u increases, the lift produced on the transitioning blades decreases. For a given lift L, with the maximal rotor angular velocity $\Omega $ given in equation (1.1),

$$ \begin{align*} \frac{\rho c R}{6}(\nu - u)^{2} = L. \end{align*} $$

Solving for u,

(3.4) $$ \begin{align} u = \nu - \sqrt{\frac{6L}{\rho c R}}. \end{align} $$

We plot the maximum speed limit as a function of lift in Figure 19. As we increase the amount of lift the blade must produce, the maximum forward velocity at which the blades are able to maintain the lift decreases. Therefore, coaxial helicopters are subject to an additional speed limit in equation (3.4) as imposed by the reduced lift on the transitioning blades. Though using four blades per disk reduces the lift demand for each blade, the problem of reduced lift on transitioning blades remains, as seen in Figure 20.

Figure 19 Plot of the maximum speed limit imposed by equation (3.4) against lift.

Figure 20 Plot of the maximum lift produced by the blades at $\theta = 90^{\circ }$ and $\theta = 270^{\circ }$ .

4 Conclusion

From the various models for the lift equation, we conclude that replacing the $\nu ^{2}$ term with $\nu |\nu |$ is recommended to capture the effect of the stall region about the rotor shaft on the retreating side (as seen in Figures 5 and 6). The lift dynamics at higher forward velocities are also better modelled by the modified velocity term. In the ideal lift situation, we find that the retreating blade produces no lift for the forward velocity $u = \nu /3$ , regardless of R. This implies conventional helicopters with one main rotor cannot accelerate to $\nu /3 \approx 114.3\,\text {m}\text {s}^{-1}$ , and struggle to handle dissymmetry of lift at speeds close to $\nu /3$ . We use equation (3.2) to find the minimum rotor angular velocity required to generate L lift. To avoid supersonic blade issues, the rotor speed $\Omega $ must decrease when the forward velocity u increases. Because of this, the overall lift is decreased, as shown in Figure 10. For the proposed coaxial design featuring two counter-rotating blades, no lift will be produced on the retreating blades. This impacts the lift generated by the helicopter, as shown in Figure 10.

Acknowledgements

The authors would like to acknowledge Hyper Q Aerospace for their support. They would also like to thank Graeme Hocking for his input and for his moderation of the related project at MISG 2020.

References

Fay, J., The helicopter, 4th edn (Hippocrene Books, New York, 1987).Google Scholar
Kulchenko, A., Kostukov, V., Verevkina, L. and Chufistov, V., “The features of flight dynamic single rotor helicopter mathematical model: for application in autopilots that based on position-trajectory algorithms”,2014 $1{9}^{\mathrm{th}}$ Int. Conf. on Methods and Models in Automation and Robotics (MMAR) (IEEE, Miedzyzdroje, Poland, 2014) 423428; doi:10.1109/MMAR.2014.6957391.Google Scholar
La Mantia, M. and Dabnichki, P., “Effect of the wing shape on the thrust of flapping wing”, Appl. Math. Model. 35 (2011) 49794990; doi:10.1016/j.apm.2011.04.003.CrossRefGoogle Scholar
Leishman, G. J., Principles of helicopter aerodynamics with CD extra (Cambridge University Press, Cambridge, 2006).Google Scholar
Liu, X., Liang, S., Li, G., Godbole, A. and Lu, C., “An improved dynamic stall model and its effect on wind turbine fatigue load prediction”, Renew. Energy 156 (2020) 117130; doi:10.1016/j.renene.2020.04.040.CrossRefGoogle Scholar
Majhi, J. R. and Ganguli, R., “Helicopter blade flapping with and without small angle assumption in the presence of dynamic stall”, Appl. Math. Model. 34 (2010) 37263740; doi:10.1016/j.apm.2010.02.010.CrossRefGoogle Scholar
Perdomo, O. and Wei, F., “On the flapping motion of a helicopter blade”, Appl. Math. Model. 46 (2017) 299311; doi:10.1016/j.apm.2017.01.055.CrossRefGoogle Scholar
Romani, G. and Casalino, D., “Rotorcraft blade-vortex interaction noise prediction using the lattice-Boltzmann method”, Aerosp. Sci. Technol. 88 (2019) 147157; doi:10.1016/j.ast.2019.03.029.CrossRefGoogle Scholar
Rozhdestvinsky, K. V. and Ryzhov, V. A., “Aerohydrodynamics of flapping-wing propulsors”, Prog. Aerosp. Sci. 39 (2003) 585633; doi:10.1016/S0376-0421(03)00077-0.CrossRefGoogle Scholar
Figure 0

Figure 1 Diagram showing the direction of rotation with the relative wind.

Figure 1

Figure 2 Diagram of the three angles of rotation for a helicopter rotor blade.

Figure 2

Figure 3 Plot of the coefficient of lift against the angle of attack for the NACA0012 airfoil for Reynolds numbers from $5\times 10^5$ to $10^6$ [8]. The x-axis is in degrees.

Figure 3

Figure 4 Lift distribution of the rotor disk in standard hover.

Figure 4

Figure 5 Lift distribution of the rotor disk given a relative wind $u = 50\,\text {m}\text {s}^{-1}$ (left) and $u = 120\,\text {m}\text {s}^{-1}$ (right).

Figure 5

Figure 6 Modified lift distribution of the rotor disk given a relative wind $u = 50\,\text {m}\text {s}^{-1}$ (left) and $u = 120\,\text {m}\text {s}^{-1}$ (right).

Figure 6

Figure 7 Plot of the lift generated over time for several forward velocities.

Figure 7

Figure 8 Plot of the lift coefficient $C_{L}$ required for vertical hover against time.

Figure 8

Figure 9 Plot of the average lift $L_{\text {mod}, T}$ over time interval $\lambda $.

Figure 9

Figure 10 Plot of the normalized average lift over forward velocity for each model.

Figure 10

Figure 11 Plot of minimum $\Omega $ for various weights with the theoretical maximum against forward velocity.

Figure 11

Figure 12 Plot of maximum speed $u_{M}$ against lift L for several blade lengths.

Figure 12

Figure 13 Plot of the average lift over the advancing half of the helicopter rotor disk against forward velocity.

Figure 13

Figure 14 Plot of the average lift over the advancing half of the helicopter rotor disk against forward velocity for several angles of attack.

Figure 14

Figure 15 Plot of the lift coefficient and the difference between lift and drag against the angle of attack.

Figure 15

Figure 16 Plot of the average lift over the advancing half of the helicopter rotor disk against forward velocity for several angles of attack.

Figure 16

Figure 17 Plot of the average lift over the retreating half of the helicopter rotor disk against forward velocity.

Figure 17

Figure 18 Plot of the average lift over the retreating half of the helicopter rotor disk against forward velocity.

Figure 18

Figure 19 Plot of the maximum speed limit imposed by equation (3.4) against lift.

Figure 19

Figure 20 Plot of the maximum lift produced by the blades at $\theta = 90^{\circ }$ and $\theta = 270^{\circ }$.