Links
Comment on page

Stability

The aim of stability analysis is to determine if a design is statically and dynamically stable.

The pitching moment coefficient and the concept of the neutral point

The moment about the y-axis of the body-fixed coordinate system affects the stability of the design. To have a design that is "statically stable" the moment about the y-axis should be restoring. What is meant by a restoring moment?
Consider the image below. As the angle of attack is increased, the moment about the y-axis
MyM_{y}
becomes more and more negative. This means that the design has a tendency to return to its original angle of attack. This tendency is called "static stability".
Static stability tendency
Since the moment about the y-axis
MyM_{y}
depends on the center of gravity location, moving the center of gravity location forward or aft will affect the stability of the design. The image blow shows the pitching moment coefficient
CmC_{m}
as a function of the angle of attack. For a stable design the derivative of
CmC_{m}
with respect to
α\alpha
must be negative
Cmα<0\frac{\partial C_{m}}{\partial\alpha}<0
.
Pitching moment coefficient versus angle of attack
Movement of the center of gravity front to aft affects the value of
Cmα\frac{\partial C_{m}}{\partial\alpha}
. The derivative is referred to as the pitch static stability derivative and must be negative if the design is to be stable. It is interesting to determine as the center of gravity is moved at what point the value of
Cmα\frac{\partial C_{m}}{\partial\alpha}
becomes 0. This point (or position of the center of gravity) is referred to as the neutral point.
If we have the values of
CmC_{m}
and
CZC_{Z}
at two different angles of attack we can determine the location of the neutral point using the following equation:
Cmcg=CmxcgcrefCZ.C_{m}^{cg}=C_{m}-\frac{x_{cg}}{c_{ref}}C_{Z}\text{.}
The above equation calculates the pitching moment about an arbitrary point
xcgx_{cg}
along the x-axis of the body-fixed coordinate system. Since we are interested in finding at what point
Cmα\frac{\partial C_{m}}{\partial\alpha}
becomes zero we take the derivative of the equation with respect to the angle of attack and equate it to zero.
Cmcgα=CmαxcgcrefCZα=0\frac{\partial C_{m}^{cg}}{\partial \alpha}=\frac{\partial C_{m}}{\partial \alpha}-\frac{x_{cg}}{c_{ref}}\frac{\partial C_{Z}}{\partial \alpha}=0
Rearranging the equation for
xcgx_{cg}
gives:
xcg=xnp=crefCmαCZαΔCmΔCZCmα=5Cmα=0CZα=5CZα=0.x_{cg}=x_{np}=-c_{ref}\frac{\frac{\partial C_{m}}{\partial \alpha}}{\frac{\partial C_{Z}}{\partial \alpha}}\approx\frac{\Delta C_{m}}{\Delta C_{Z}}\approx\frac{C_{m}\rvert_{\alpha=5}-C_{m}\rvert_{\alpha=0}}{C_{Z}\rvert_{\alpha=5}-C_{Z}\rvert_{\alpha=0}}\text{.}
Solving the above equation will give the location of the center of gravity for which
Cmα=0\frac{\partial C_{m}}{\partial\alpha}=0
. This location is referred to as the neutral point.
If the the center of gravity is behind the neutral point, the design will be statically unstable. If the center of gravity is in front of the neutral point, the design will be statically stable.
The static margin
SMSM
is a measure of how statically stable the design is. The static margin can be calculated from the from the following equation:
SM=xcgxnpcrefSM=\frac{x_{cg}-x_{np}}{c_{ref}}

Example

The following example shows how to determine the neutral point location of a fixed-wing UAV with reference area
𝑆ref𝑆_{ref}
of
0.7729 m20.7729\text{ }m^2
and reference chord
crefc_{ref}
of
0.2544 m20.2544\text{ }m^2
. A body-fixed coordinate system is placed at the nose. The reference speed is
Vref=25 m/s|V_{ref} |=25\text{ }m/s
and the reference density, pressure and viscosity are evaluated at an altitude of
h=0 mh=0\text{ }m
(ISA).
Both fixed-wake and free-wake simulations of the fixed-wing UAV were performed at
α=0\alpha=0
and
α=5\alpha=5
deg.
Fixed-wake simulation of the UAV at zero angle of attack
Free-wake simulation of the UAV at zero angle of attack
Fixed-wake simulation of the UAV at 5 deg angle of attack
Free-wake simulation of the UAV at 5 deg angle of attack
The table below summarises the values of the
CZC_{Z}
and
CmC_{m}
coefficients obtained from the simulations.
Simulation type
α\alpha^\circ
CZC_{Z}
CmC_{m}
Fixed-wake
0
−0.5357
−1.1156
Fixed-wake
5
−1.0417
−2.2843
Free-wake
0
−0.2774
−0.6383
Free-wake
5
−1.0280
−2.2641
The neutral point is calculated from the previously derived expression:
xnp=crefΔCmΔCZ=crefCmα=5Cmα=0CZα=5CZα=0x_{np}=-c_{ref}\frac{\Delta C_{m}}{\Delta C_{Z}}=-c_{ref}\frac{C_{m}\rvert_{\alpha=5}-C_{m}\rvert_{\alpha=0}}{C_{Z}\rvert_{\alpha=5}-C_{Z}\rvert_{\alpha=0}}
For the fixed-wake simulations the neutral point is located at
x=0.5866 mx=-0.5866\text{ }m
whereas for the free-wake simulations the neutral point is located at
x=0.5501 mx=-0.5501\text{ }m
. The absolute difference amounts to approximately 6.5 percent.
The free-wake simulations give more accurate prediction for the neutral point location as they account for the wake deformation. This is confirmed by forced-pitch oscillation simulations. The neutral point locations obtained from forced-pitch oscillation simulations and from free-wake simulations agree well.
Once the neutral point location is determined, the location of the center of gravity can be calculated using the following equation:
xcg=xnp+SMcref,x_{cg}=x_{np}+SMc_{ref},
where
SMSM
is the static margin. The greater the static margin the more statically stable the UAV is in pitch. For a 10% static margin (
SM=0.1SM=0.1
)
xcgx_{cg}
is
0.5247 m-0.5247\text{ }m
.