Obtaining stability derivatives from forced oscillations
This article explains how to obtain the stability derivatives of a design from forced oscillations
Consider a simple wing-tail design. The design has a wing span
, wing area
, and wing chord of
. The image below shows the mesh of the design.
Mesh of a wing-tail design
A body-fixed coordinate system is placed at the leading edge of the wing. The reference speed is
and the reference density, pressure and viscosity are evaluated at an altitude of
Before performing forced oscillation simulations the location of the neutral point of the wing-tail configuration must be determined. Two simulations at different angles of attack are required to determine the neutral point location. Please refer to the Stability article for more information on how to determine the neutral point.
The table below shows the
values used for determining the location of the neutral point,
The location of the neutral point,
, can be obtained from the following equation:
After substituting the
values from the table into the above equation,
. The minus sign indicates that the neutral point is aft of the origin of the body-fixed coordinate system. The following equation can be used to determine the centre of gravity for a specific static margin,
For a statically stable design the centre of gravity is ahead of the neutral point.
Once the centre of gravity location is obtained, forced oscillation simulations can be performed. Remember to change the x_cg value in the .conf file and to re-run the preprocessor before proceeding. The preprocessor will move the origin of the body-fixed coordinate system to the new location specified by x_cg, y_cg, and z_cg. Please refer to the Configuration file section of the documentation for more information.
A motion file describing the translation and rotation of the body-fixed coordinate system is required to perform a forced oscillation simulation. Please refer to the Custom motion section of the documentation for additional information. To obtain the longitudinal stability derivatives of the design, a forced pitch oscillation motion is required. The pitch angle during the oscillation is defined by the following equation:
is the pitch angle,
is the amplitude of the pitch oscillation,
is the angular frequency of the oscillation (
is the time. In this example
. The pitch rate
can be obtained by differentiating the
equation with respect to time:
Appropriate time step size and number of time steps must be selected to complete at least 2 cycles. In this example
. The image below shows
. Note that
is divided by 100 so the difference between
is easy to see. In the actual motion file
is not divided by 100.
Generated pitch angle and pitch rate
The video below shows the development of the wing-tail configuration wake during the forced oscillation simulation.
Wing-tail design forced pitch oscillation
After running the simulation the aerodynamic coefficients can be plotted against the angle of attack
. The image below shows the
Aerodynamic coefficients versus angle of attack
The locations of the minimum and maximum
are shown in the images. The
coefficient exhibits a quadratic behavior, which is expected.
stability derivatives can be obtained in multiple ways. Here the
stability derivatives are obtained with the least squares method. To obtain the
stability derivatives the following system of equations is defined:
Note that in the above system of equations
is a vector of ones which length is equal to the length of the
vectors. The system is overdetermined - it has more equations than unknowns. However the least squares method can be used to solve the overdetermined system. Solving the system gives the values for
. The same process is repeated for the
A polynomial fitting method is used for
in order to capture its quadratic behavior. A third order polynomial in
is selected. The polynomial is given by the following equation:
Unlike the least squares method the polynomial fitting method cannot be used to determine
can be determined from the following equation:
is the reduced frequency, given by
. The table below shows the obtained stability derivatives:
The wing-tail design is both statically and dynamical stable as
The time history of the aerodynamic coefficients can be used to verify the obtained stability derivatives. The stability derivatives are used to describe the time history of
The above equations can be compared against
obtained from the simulation. The variation of angle of attach
and the pitch rate
with respect to time is known and is given by:
The image below shows the comparison of
Comparison of the aerodynamic coefficients from the simulation and from the obtained aerodynamic derivatives
As seen from the image the aerodynamic coefficients agree well. The minor differences from
are due to the simulation transients (the time it takes for the flow to develop).