Helicopter Flap Dynamics
In this section we provide a derivation of the flapping dynamics for three simple scenarios.
We’ll see why cyclic controls impart max pitch
roughly \(90^o\) ahead of where the blade reaches peak flapping
(e.g. why left lateral cyclic provides max pitch around \(\Psi=0^o\) to
flap the rotor up at \(\Psi=90^o\)). We’ll see how this angle shrinks with the addition of a flap spring or the outward movement of the flap hinge.
Click here for definitions of symbols and abbreviations used here and throughout our site.
First, let’s recall a few equations from Newtonian physics.
The rotational analog of Newton's \(F=m\ddot{x}\) for a blade flapping about a hinge,
with flap angle \(\beta\), is \(M=I \ddot{\beta} \) or $$\begin{equation} \ddot{\beta}=\frac{M}{I} \label{eq:betadd} \end{equation}$$ where \(M\)
is the moment about the flap hinge and \(I\) is the flapping inertia of the blade about the flap hinge.
The centrifugal force
acting on a section of the blade with length \(dr\) and mass \(mdr\) is \(mdr\Omega^2 r\).
This force acts radially outward, away from the shaft. The component normal to the blade, which produces a moment about the hinge,
is \(mdr\Omega^2 r \sin\beta\).
Scaling this by the moment arm (\(r-a\)) about the flap hinge and summing over the blade outboard of the hinge (at \(r=a\)) gives
$$\begin{equation}M_{CF}= -\int_a^R mdr\Omega^2 r(r - a) \beta . \label{eq:cf} \end{equation}$$
We've used the small angle assumption \(\sin\beta \approx \beta\) since flapping is normally limited to just a few degrees. This simplifies the arithmetic significantly.
The flap inertia of the blade about the flap hinge is $$\begin{equation} I = \int_a^R m(r-a)^2dr. \label{eq:inert} \end{equation}$$
If a torsion spring is present at the flap hinge, it provides an additional moment $$\begin{equation}M_k=-k\beta \label{eq:tors} \end{equation}$$ where \(k\) is the stiffness of the spring.
Now we'll use this information to derive the flap dynamics in three simple scenarios.
Scenario 1: Free flapping about a hinge at the center of rotation
This is the simplest case, not requiring excessive arithmetic. A blade is free to flap with no restraint, \(k=0\). The flap hinge is located at the center of rotation, \(a=0\). This is roughly the situation for a teetering rotor. We'll focus on the centrifugal force here. Using equations \ref{eq:betadd} - \ref{eq:inert} above, the dynamics are governed by $$\begin{equation} \ddot{\beta}=\frac{M_{CF}}{I}\approx \frac{-\int_0^R\beta \Omega^2r^2mdr}{\int_0^Rr^2mdr}= -\beta\Omega^2. \label{eq:s1dq} \end{equation}$$
This is analogous to the differential equation for a spring and mass system with the known solution $$\begin{equation} \beta (t)= \beta_m\cos(\Omega t+ \Psi_0) \label{eq:s1} \end{equation}$$ where \(\beta_m\) (the max flapping angle) and \(\Psi_0\) depend on the initial conditions.
In other words, the blade undergoes periodic motion with frequency exactly equal to the frequency of rotation.
This “coincidence” means that the blade returns to the same flap angle after exactly one revolution.
If the flap angle at \(\Psi_1\) is \(\beta_1\) then the flap angle will return to \(\beta_1\) on the next rotor revolution when the blade is at \(\Psi_1\).
In the context of rotors, frequencies are often given in units of “per rev” rather than per second.
What we've just found is that the flap frequency is 1 per rev (denoted 1P).
Scenario 2: The addition of a flap spring
A flap spring provides a moment proportional to the flap angle, \(M=-k\beta\).
Since this increases the moment pushing the blade back to \(\beta =0\), without affecting inertia, it increases \(\ddot{\beta}\) and the flap frequency.
The math is straightforward: add \(-k\beta\) to the numerator of \ref{eq:s1dq} to get $$\begin{equation} \ddot{\beta} \approx -\beta \Omega^2 - \frac{k\beta}{I} = -\beta (\Omega^2 + \frac{k}{I}). \label{eq:s2dq} \end{equation}$$
Hence the frequency increases from \(\Omega\) to \(\omega_k= \sqrt{\Omega^2+\frac{k}{I}}\) and the solution is $$\begin{equation} \beta (t)= \beta_m\cos(\omega_k t+ \Psi_0). \label{eq:s2} \end{equation}$$
Although this frequency is larger than 1P, the rotor undergoes motion that's dominantly 1P. This is because each blade encounters large, 1P aerodynamic forces.
Scenario 3: Offsetting the flap hinge
In this scenario we place the flap hinge a distance \(a>0\) away from the shaft. We consider the flapping \(\beta\) outboard of \(a\) and hence only consider the moment \(M\) about the hinge (due to CF) and the inertia \(I\) about the hinge and outboard of the hinge. We’ll assume the mass per unit length \(m\) is constant along the blade (so that \(m\) can be pulled out of the integrals and cancelled) – even so the arithmetic is more difficult than the prior scenarios. The differential equation is $$\begin{equation} \ddot{\beta} \approx \frac{-\int_a^R\beta\Omega^2r(r-a)mdr}{\int_a^Rm(r-a)^2dr}= -\beta \Omega^2\frac{ \int_a^Rr^2dr – a\int_a^Rrdr }{\int_a^Rr^2dr -2a\int_a^Rrdr - a^2\int_a^Rdr} \end{equation} $$
Now adding \(0=-a\int_a^Rrdr -a^2\int_a^Rdr +a\int_a^Rrdr +a^2\int_a^Rdr\) to the numerator and writing it as the denominator plus the remaining terms this becomes \( \ddot{\beta} = -\beta \Omega^2(1 + s)\). With several more steps of arithmetic you can show that \(s = \frac{3a}{2(R-a)}\). Putting this together we get the differential equation $$\begin{equation} \ddot{\beta} = -\beta \Omega^2(1 + \frac{ 3a }{ 2(R-a) }). \end{equation} $$
Of course this is again a periodic motion given by $$ \begin{equation} \beta (t) = \beta_M \cos ( \Omega \sqrt{1+s} t + \Psi_0). \end{equation}$$
Clearly the frequency increases with \(a\). This is because the "centrifugal moment" due to a section of blade is proportional to \(r(r-a)\) while the inertia due to the same section is proportional to \((r-a)^2\). Hence \(I\) decreases more than \(M\) as the hinge is moved outboard, increasing \(\ddot{\beta}\) and the flap frequency.
Advantage of a flap hinge offset
We’d like to point out a feature of a hinge offset. The centrifugal force acting on a blade created the moment \(M\) described above, but it also applies a force to the hinge radially outward, parallel to the flapped blade. Without a hinge offset, and with roughly opposite flapping on opposite sides of the azimuth (\(\beta (\Psi) = -\beta( \Psi + \pi )\), the force from one blade cancels the force from the blade on the opposite side of the azimuth and therefore can be neglected. However, with a hinge offset, components of these forces parallel to the shaft create a moment about the hub. The moment is exactly as would be desirable to enhance the moment created by tilting the rotor thrust. For example, in picture below, the upward displaced aft blade pulls the aft hinge up (and aft). The the front blade pulls the hinge down (and forward). The forward and backward components cancel, but the up and down components (multiplied by the moment arm \(a\)) add together to create a nose down pitch moment. Furthermore, this moment is not reliant on thrust. This enhances the safety of the helicopter, giving the pilot roll and pitch control in low thrust conditions. Crashes in helicopters with teetering rotors have been attributed to control loss in low thrust conditions.
|