Theodorsen Function

where the generalized forces are given in Eqs. (4.22). The function C(k) is a complex-valued function of the reduced frequency k, given by

where W,"(k) are Hankel functions of the second kind, which can be expressed in terms of

Bessel functions of the first and second kind, respectively, as Hf\k) = J„(k) - iY„(k).

The function C(k) is called Theodorsen's function and is plotted in Fig. 4.9. Note that C(k) is real and equal to unity for the steady case (i.e., for k — 0). As k increases, one finds that the imaginary part increases in magnitude while the real part decreases. As k tends to infinity, C(k) approaches 1 /2. However, for practical situations k does not exceed values of the order of unity. Hence, the plot in Fig. 4.9 only extends to k — 1. When any harmonic function is multiplied by C(k), its magnitude is reduced and a phase lag is introduced.

A few things are noteworthy concerning Eqs. (4.78). First, in Theodorsen's theory the lift-curve slope is equal to lit. Thus, the first of the two terms in the lift is the circulatory

Figure 4.9 Plot of the real and imaginary parts of C(k) for k varying from zero, where C(k) = 1, to unity.

lift without the effect of shed vortices multiplied by C(k). The multiplication by C(k) is a consequence of the theory having taken into account the effect of shed vorticity. The second term in the lift as well as the pitching moment are noncirculatory, depending on the acceleration and angular acceleration of the airfoil. The circulatory lift is the more significant of the two terms in the lift.

For steady flow the circulatory lift is linear in the angle of attack, but for unsteady flow there is no single angle of attack since the flow direction varies along the chordline as the result of the induced flow. However, just so we can discuss the concept for unsteady flow, it is possible to introduce a so-called effective angle of attack. For simple harmonic motion it can be inferred from Theodorsen's theory that an effective angle of attack is a = C(k)

As we will show in Section 4.5.2 by comparison with the finite-state aerodynamic model introduced therein, a is the angle of attack measured at the three-quarter chord based on an averaged value of the induced flow. Recall that in steady-flow aerodynamics, the angle of attack is the pitch angle 0. Here, however, a depends on 0 as well as on /;, 0, and k. Because of these additional terms and because of the behavior of C(k), we expect changes in magnitude and phase between 0 and a carrying over into changes in the magnitude and phase of the lift and pitching moment relative to that of 0. Indeed, the function C(k) is sometimes called the lift-deficiency function because it reduces the magnitude of unsteady lift relative to steady lift. It also introduces an important phase shift between the peak values of pitching oscillations and corresponding oscillations in lift and pitching moment. An approximation of Theodorsen theory in which C(k) is set equal to unity is called a "quasi-steady" thin-airfoil theory. Such an approximation has value only for cases in which k is restricted to be very small.

Theodorsen's theory may be used in classical flutter analysis. There the reduced frequency of flutter is not known a priori. One can find k at the flutter condition using the method described in Section 4.3. Theodorsen's theory may also be used in the k and p-k methods, as described in Sections 4.4.1 and 4.4.2, respectively.

4.5.2 Finite-State Unsteady Thin-Airfoil Theory of Peters et al.

Although the Theodorsen theory is an excellent choice for classical flutter analysis, there are situations in which an alternative approach is needed. First, frequently one needs to calculate the modal damping in subcritical flight conditions. Second, there is a growing interest in the active control of flutter, and design of controllers requires that the system be represented in state-space form. To meet these requirements, one needs to represent the actual aerodynamic loads (which are in the frequency domain in Theodorsen's theory) in terms of time-domain differential equations. Finite-state theories approximate the actual infinite-state aerodynamic model to within engineering accuracy. One such approach is the finite-state, induced-flow theory for inviscid, incompressible flow of Peters et al.

Consider a typical section of a rigid, symmetric wing, shown in Fig. 4.2, and the additional vectorial directions defined in Fig. 4.10. To begin the presentation of this theory, we first relate the three sets of unit vectors:

1. a set fixed in the inertial frame, ii and ij, such that the air is flowing at velocity

2. a set fixed in the wing, bi and b2, with bi directed along the zero-lift line toward the leading edge and bj perpendicular to bi,

3. a set ai and kn associated with the local relative wind vector at the three-quarter chord, such that aj is along the relative wind vector and a2 is perpendicular to it, in the assumed direction of the lift.

The relationships among these unit vectors can be simply stated as and bi b2

Induced-flow theories approximate the effects of shed vortices based on changes they cause in the flow field near the airfoil. Thus, the velocity field near the airfoil consists of the

^ lift direction

^ lift direction zero-lift line direction of the velocity of T relative to wind

Figure 4.10 Schematic showing geometry of the zero-lift line, relative wind, and lift directions.

zero-lift line direction of the velocity of T relative to wind

Figure 4.10 Schematic showing geometry of the zero-lift line, relative wind, and lift directions.

freestream velocity plus an additional component to account for the induced flow. Although the induced flow varies throughout the flow field, we will approximate its value near the airfoil as an average value along the chordline. Thus, the local inertial wind velocity is written approximately as — Ui\ + /.obi, where /,« is the average induced flow (perpendicular to the airfoil zero-lift line). According to classical thin-airfoil theory, one should calculate the angle of attack using the instantaneous relative wind velocity vector as calculated at T. To represent the relative wind velocity vector at T, one can write the relative wind vector (i.e., the velocity of the wing with respect to the air) as Wai and set it equal to the inertial velocity of T minus the inertial air velocity, that is,

Wai = \'T ~ (-Uii + A0b2) = Yp + Uii - A0i>2, where vr is the inertial velocity of the three-quarter chord, given by vT — vp + Ob'i x rPT, and rpx is the position vector from P to T. From Fig. 4.2 one finds that r pt —

Thus,

+3 -1

Responses

• WEGAHTA
How to plot theodorsen function?
6 years ago
• nebay
What is the lift deficiency function?
1 year ago