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
Given a helicopter rotor disk of n blades, the azimuth angle $\theta $ for blade k of n at time t is given by
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
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.
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.
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.
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
where $\nu $ is the speed of sound. Rearranging for $\Omega $ , we have
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
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
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
To study the effect of a relative wind, we model the linear velocity with the equation
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 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.
3 Results
3.1 Revised lift equation
In this section, we compute the lift generated by a helicopter rotor blade with the modified equation
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,
3.1.1 Retreating blade special case
For retreating blades ( $\theta = 180^{\circ }$ ), the lift generated simplifies to
Given $\Omega =(\nu - u)/R$ to ensure maximum RPM, the lift becomes
To find u such that $L_{R}(u) = 0$ , we solve
Given $u < \nu /2$ ,
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
Letting $\epsilon = \cos (\theta _{c})$ and using $\Omega = (\nu - u)/R$ ,
Since $\epsilon < 0$ , by expansion and factoring,
Therefore, we have $\epsilon = (\nu - u)/2u$ . Upon substitution,
Therefore, the helicopter blades will produce negative lift for azimuth angles, satisfying
3.2 Modified lift equation
To recapture the loss of lift on the retreating rotor blades, we modify the lift equation to be
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.
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
Upon computation,
With $\Omega = (\nu - u)/R$ , the lift becomes
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
With $\Omega = (\nu - u)/R$ , the coefficient of lift becomes
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$ .
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
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.
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
Upon simplification,
For the modified lift, we compute
Finally, to reconcile the negative lift from the modified velocity model $\nu |\nu |$ , we also compute the integral
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.
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
which simplifies to
Solving for $\Omega $ ,
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.
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
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.
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
Since $\Omega r + u\cos (\theta )> 0$ for advancing rotor blades and the integral is even, we may simplify the expression to
Upon simplification,
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 $ .
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.
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 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 }$ .
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
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.
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
where $L_{\text {std}}(u, \theta )$ is the lift generated by the blade at forward velocity u and azimuth angle $\theta $ , given by
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
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
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
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$ .
Therefore, in the case that $u> u_{1}(L)$ , the coefficient of lift will take the form
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,
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),
Solving for u,
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.
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.