Airfoil Analysis
Airfoil Analysis
Airfoil Analysis
NASA
Brooks (ref. 7), Montoya (ref. 8) and Yu (ref. 9), among others, using the
classical Euler-Bernoulli theory allowing for various complexities. Using the
Timoshenko beam theory, which allows for shear deflection and rotary inertia
effects, and incorporating further complexities, the coupled equations of
motion were derived by Carnegie (ref. 10), Rao and Rao (ref. 11) and
Subrahmanyam, Kulkarni and Rao (ref. 12) among others. Recent research in helicopter and turbo-prop blades indicates that the linear equations of motion are
inadequate for establishing flutter boundaries. Geometric nonlinear theory is
to be used with various degrees of complexity. Some such theories were presented by Hodges and Dowel 1 (ref. 13), Rosen and Friedmann (ref. 14), and Kaza
and Kvaternik (ref. 15) among others.
Several methods of solution of vibration problems are available. Broadly,
these are classified as belonging to either the continuum model approach or to
the discrete model approach. A comprehensive review of the work done using the
various methods in blade vibration area can be found in the review articles by
Rao (ref. 16) and Leissa (ref. 17) among others.
An examination of the published literature reveals that among the various
methods used for solution of the differential equations of motion, the firstorder finite-difference method has attracted the greatest attention while
relatively few works exist which use higher-order finite-differences. Even in
those few works using second-order central differences either the approximations made to eliminate fictitious stations were not good enough to give
accuracy superior to first-order theory results (ref. 18) or a method of using
forward or backward differences at beam boundaries was used to avoid the fictitious stations (ref. 19). An improvement to the second order finite-difference
method was made by the present authors and this refined method was applied to
uncoupled vibration analysis of tapered cantilever beams in (ref. 20), and to
pretwisted tapered blades executing coupled bending-bending vibrations in
(ref. 21). However, the applicability of the method for vibration analysis of
pretwisted asymmetric blades has not been investigated.
Another area, currently receiving considerable attention, is the utilization
of mixed or hybrid variational techniques for vibration analysis. Among these
mixed variational principles, the Reissner method (ref. 22) has received the
greatest attention. For solution of the dynamic vibratory problems, the usual
way of formulating the dynamic principle is to consider the difference of the
kinetic energy and the Reissner functional, consisting of contributions from
potential and complementary energies together with the work done by the body
forces and surface tractions (ref. 23), and to minimize the same according to
the Ritz process. One thus has the flexibility to vary displacements and
stresses independently and can obtain simultaneously good distributions of all
the independent parameters. Another more flexible way of formulating the
dynamic variational functional is demonstrated by Dean and Plass (ref. 24)
where the kinetic energy effects were incorporated by using linear momentum
density in conjunction with the Reissner functional. This way to include
angular momentum and rotary inertia effects was demonstrated by the present
first author's earlier work (refs. 25 and 26), and it was observed that both
mixed variational principles (refs. 22 to 24) lead to identical results if
identical shape functions and beam kinematics were used in the analysis. However, attempts to incorporate coupling between the flexural and torsional
degrees of freedom were not made so far. An important advantage offered by
the Dean and Plass method is that one can handle the gyroscopic forces in the
vibratory problem in a straightforward manner. Since this variational principle contains products of linear or angular momenta and the corresponding
velocities in the kinetic energy equivalent, and since the gyroscopic forces
also are functions of velocities, the variational functional possesses consistently first order time derivative functions. Unlike the Reissner or potential
energy methods, which require either a transformation procedure or other more
complicated methods for solution (refs. 27 and 28), the Dean and Plass method
leads directly to a standard eigenvalue problem.
The objectives of the present work are: (1) to apply the first, and second
order finite-difference methods, the Dean and Plass method and the Reissner
method for the general case of vibration with coupling between bending and
torsional degrees of freedom; (2) to assess the computational advantages and
disadvantages of each method; and (3) to study the effects of pretwist and
asymmetry in coupling the component modes. Numerical results from all four
methods will be presented and comparisons will be made to those from other
existing theoretical and experimental results (refs. 29 to 33).
SYMBOLS
B,C,D
matrices
A.j,B..,...I'..
B-j
torsional stiffness
C,
warping rigidity
c.f.
Young's modulus
f-,g.,h-
shape functions
modulus of rigidity
h
I , ,
I , ,
I , ,
Ixx
I
Icf
IP
IR
ID,IR
i,j
k
L
M ,M
length of blade
bending moments
n
Pn
r* ,ry
rx,rv etc
r
xo'rvo
coordinates
C(
rcf
SI
_
TI-
t
T
TQ9
time
kinetic energy
twisting moment
U.j
UVA,U..,U,
y i.
v
displacement field
displacement in x, y, and z directions
volume
P.
^xl'^vl
surface fractions
directions, respectively
PQ0
x,y
x,,y,
^'Y'^l'^1
x,x, ,y,y, ;xx,yy coordinates axes through centroid and center of flexure
respectively
zz
longitudinal axis
z
e. .
'J
e
p
ii
warping function
Y
n
( )'
(*)
pretwist angle
axial fractional length z/L
prime denotes differentiation with respect to z
dot over a parameter represents differentiation with respect
to time, t
successive differentiation with respect to z and t
(*)'
; u
+rU > Avi
i 'J
i etc>
J
-
x,y,e
3<|)
r 3Ui
relative amplitudes
DEVELOPMENT OF FREQUENCY EQUATIONS
(2)
The coefficient aj_, b^ ... hj_, 32, t>2 ^2 a3 ^3 are all functions
of n and are presented in appendix B.
Q
Q
Q
Q
Q
Q
E
<V
=p
Q
ti
E
Q
<V
(4)
In the preceding equation, A,, B, , D and E. are all square matrices of order
(nxn), {x-j}, {y-j} and {e-j} are column matrices containing the linear and
torsional displacements of the n-stations and F, G, H, I, J, K and L are
diagonal. The matrix Q is a null matrix. Further, each submatrix A, B,
or D has a band width of five for the first-order theory and seven for secondorder theory. The submatrix has a band width of three for first-order theory
and five for second-order theory. The matrices A, E, , D and are
nonsymmetric in first and second order finite-difference applications. For
brevity elements of these matrices are not presented, but one can easily develop
these matrices with the help of references 20 and 21.
Frequency Equation from the Dean and PI ass Principle
The Dean and Plass dynamic variational functional is developed in
appendix A and has the functional relationship of the form
3 , p }M }M }j ,n,t)
(5)
The parameters in the functional given by equation (5) are independent of each
other and depend only on n and t.
The following shape functions in series form (refs. 23 and 25) are assumed
for solution of the coupled bending-bending-torsion vibration problem:
Z A.
2
3 i
n - R2 n + R3 n ) n
G(t) =
PX1 =
B. f, (n) H(t)
*
i=l
A. f.(n)
6(t)
y= E
c1 f1 (n) G(t)
Pyl - E
D. f . (n) H(t)
k
k
1
Mxx = E E.
R
(1-n)
G(t)
=
E
4
i=l 1 H
i=l
My = E
e-
E.g.(n)
6(t)
'n
Fi 91 (n) 6(t)
k
.
k
1
E G.
n
(1-R
n)
6(t)
=
E
b5
i=l 1
i=l
k
TQ = E
k
Pe = E
G-h^n)
G(t)
n n
H. g. (n) 6(t)
^ h1 (n) H(t)
(6)
The parameters RI, R;? . . . R$t G(t), H(t) etc. are defined in appendix B.
The shape functions given in equation (2) satisfy the boundary conditions
(ref. 12) applicable to a cantilever beam fixed at n = 0 and free at n = 1.
Substituting equations (6) and their appropriate derivatives into equation
(A13), eliminating the time dependence according to the averaging procedure
"'Pn
lDdt
(7)
and applying the Ritz process to minimize Ip with respect to the arbitrary
parameters, that is,
8l
a!D aID
o
=
=
~3AT "aBT ' * ' TTT= 0>
k = If2f
----K>
^8^
A + Pn JB = 0
In equation (9), A and B are symmetric square matrices.
of these matrices are not presented here.
(9)
ds
R -Iff
(10)
-rv
2jr/P
v v -v
- (I
2MI IMI 1.I
4 *
x y xlyl
,2
2E
+ | |MI
i
y xlxl dz dt
(11)
Using similar shape functions as were used for the Dean and Plass method and
proceeding on the same lines as described earlier, one can obtain the frequency
equation which can be written as
C~ + P^D
n~ = 0
(12)
Matrices C and JD are symmetric, the elements of which are not presented here
to save space.
METHOD OF SOLUTION
The eigenvalue problems given by equations (4), (9) and (12) were solved
by using a standard eigenvalue extraction routine. This library routine determines all the eigenvalues and the associated eigenvectors. For the eigenvalue
problem formulation using mixed variational approaches, certain integrations
were to be performed numerically. This was accomplished by using a 15 point
Gaussian quadrature formula. The computer programs developed in FORTRAN language were run on an IBM 370 computer at the NASA Lewis Research Center. The
lowest ten coupled frequencies and associated mode shapes were determined for
studying the effects of asymmetry and pretwist in coupling the modes using
typical airfoil blade data (ref. 2). Results are presented and discussed in
what follows.
RESULTS AND DISCUSSION
pretwist angles. The modes were then classified as belonging to a basic mode
category by examining the associated mode shapes in relation to either the
uncoupled mode shapes or to the coupled bending-bending mode shapes. These
results are tabulated in tables VI to VIII. Thus, in table VI, for blades with
rxo.= 0.193 mm, ryo = 0 are presented, in table VII those for rxo = 0,
r = 1.1938 mm are shown while in table VIII, the results of general asymmetric airfoil blade with rxo = 0.193 mm, ryo = 1.1938 mm are listed for
various angles for pretwist. It may be notea that the frequencies are not
necessarily in an increasing order as the mode number is increased since the
modes are classified as belonging to a basic mode category. Only one set of
coupled bending-bending-torsion mode shapes are presented in figure 4 for the
case Y = 89 for brevity. The general coupling trends as observed from these
results are presented below.
(1) For pretwisted blades of doubly symmetric cross section, strong
coupling exists between the second uncoupled flapwise mode (F2) and fundamental
uncoupled chordwise mode (Si). As pretwist increases, the lower frequency
decreases and the higher frequency increases to form the corresponding coupled
frequencies. Similar coupling between the uncoupled modes S2 and F5 exists
but the magnitude of frequency variations are more severe for the set (F2, SI)
than for the set (S2, F5). This frequency variation trend is consistent with
earlier observations reported in references 23 and 35. The coupled bendingbending mode shapes obtained in the present investigation are compared to those
in reference 33. Consistent trends have been observed in all the results
compared.
(2) For untwisted blades with asymmetry about only one plane, coupling
between bending in one principal direction and torsion occurs. Results for
such a case with ry0 = 0 are shown in column 3 of table VI, and for the case
r = 0 in table Vn. From the results presented in table VI, one can
observe that the frequencies corresponding to (F2, Tl) and (F3, T2) couple in
such a way that the lower uncoupled frequency in each pair is reduced and the
higher one is increased to form the corresponding coupled bending-torsion frequencies. Although not shown here, similar coupling trends exist for the pairs
(F4, T3) and (F5, T4). Similar coupling trends are observed for the coupled
modes having coupling between chordwise bending and torsional motions. This
coupling trend is shown by the frequency pair (SI, Tl) in table VII. Although
not shown, the frequency pair (T3, S2) also shows similar coupling trends.
(3) For pretwisted blades with asymmetry about only one plane, the coupling trends are combinations of those discussed in (1) and (2) above. A careful examination of the results presented in tables VI and VII for various
values of pretwist angles in relation to the corresponding ones in table V
indicates that a specific trend in the three close frequency values corresponding to F2, SI and Tl modes exists. The frequencies corresponding to the pretwisted symmetric blade listed in table V will be further coupled because of
asymmetry in such a way that the lower two frequencies will be reduced and the
highest one will be increased to result in the coupled bending-bending-torsion
frequencies. Other secondary coupling trends are observed for the frequency
pair (T2, F4) when the frequencies are very close to one another, in table
VII, and for all pretwist angles in table VI. The reason for this is that the
flap-torsion coupling is more important for blades having a symmetry about the
chord (table VI) than for those with asymmetry about the chord (table VII).
11
14
15
APPENDIX A
DEVELOPMENT OF THE DEAN AND PLASS VARIATIONAL FUNCTIONAL
The Dean and Plass dynamic variational functional (ref. 24) can be
written as
(A2)
Uy = y + x.e = y1 + x^
(A3)
(A4)
For the computation of the kinetic energy, one can disregard the effects of
longitudinal inertia effects in the total kinetic energy without great loss of
accuracy in the final results. It will be assumed that the longitudinal
inertia effects are negligible and thus the linear momentum densities about
x and y directions take the following forms:
Px = P Ux = P (xj - y^e)
(A5)
Py = P Uy = P (yx + 8le)
(A6)
Equations (A5) and (A6) contain terms associated with linear velocities x^
and y\ together with the angular velocity e. One can explicitly define
the linear and angular momenta by defining the linear momenta about x^X}, y^y
centroidal directions and angular momentum due to elastic twist independently
as follows:
Pxl = PX!
16
(A7)
(A8)
yl -
(A9)
Thus
.t. dv =/
U
/[pxl
- p^e) ( Xl -
Yle)+
+ pxje) (y x + x^) dA dz
(P
ML
= {L
[PxlA*l
+ P
yl A*l
(A10)
+ P A
e J dz
and
j0n?-/
v
'2
P A
xl
+ P
2
ylA
22
P A
+ e
dz
(All)
The Reissner functional, IR$ can be reduced to the following form for the
coupled bending-bending-torsion case neglecting shear deformation, body forces,
surface fractions and warping rigidity effects:
dz
M x - T
V
^
2E v(Ixlxl
, ,1ylyl
, , - Ixlyl
, .)
y
(A12)
Making use of equations (A10) to (A12) in equation (Al), the Dean and Plass
dynamic variational functional becomes
x + r e) + P -.AVJ
(y
y ' yl
P 0 Ae -
M 2 I y, y, - 2M M.I . . + M 2 I.
" +Mx " - T e
' +7
e r + x vlvl
v xlvl
v xlxl
-~ xx y
xiyi
y
xixi
My
T
xlxlylyl
dz
(A13)
Treating Pxi, Pyi, P0, Mx, My, T0, x, y and e as independent variables
(independent of one another but dependent on z and t) and performing the variational process (ref. 12)
=0
17
(A14)
\
xl =
(*+ ry)
=p
*l
yl = p
e = >< I X lxl
=
+ I
ylyl)e/A='V/A
i i - My ]
_ (Mx Iylyl
(A15)
(
x Vxlxl - Vxlyl*
e = T0J2C
3Z
It
Tt
2
(A16)
3Z
-TT- (M ) = 0 at z = 0, L or x prescribed
T~
(^Xww) = 0 at z = 0, L or y prescribed
oZ
(A17)
T9 = 0 at z = 0, L or e prescribed
Equations (A15) relate the linear and angular momenta to respective linear or
angular velocities, curvatures to bending moments and rate of twist to twisting
moment and torsional rigidity. These relations can alternatively be looked
upon as stress-strain relations. Equations (A16) represent the conditions of
equilibrium or motion. The first two of equations (A16) relate the time rate
changes of linear momenta per unit length to corresponding rates of loading
(forces per unit length). The last of equations (A16) equates the time rate of
change of angular momentum per unit length to the twisting moment (torque) per
unit length. Equations (A17) are the boundary conditions associated with the
present formulation. Assuming harmonic solutions for eliminating the time
18
dz
El xlxl
xlyl
,2
d2x
El
dz
dz
Iz'
xlyl dz
,x + PAPnrye
(A18)
- pIcfPne
- - pArxPn*
Equations (A18) agree with those derived earlier by Carnegie (ref. 10) for the
specialized case considered here.
19
APPENDIX B
SHAPE FUNCTIONS, GEOMETRIC PROPERTIES, AND COEFFICIENTS
OF DIFFERENTIAL EQUATIONS
Shape functions:
fi(n) = n1 (RI n - R2 n2 + RS n3)
9i = 4 (1-n)1
hi = n1 (1 - RS n)
G(t) = sin Pt
H(t) = cos Pt
Geometric properties:
rx = rxo cosYn - ryosin Yn
ry = ryo cos Yn - i"xosln Yn
P
!xlxl = !yy sin Yn + IXx
Ixlyl.,LaCiisin
sin
C2
(Iyy - IXx)
sn
0*2 = Iyy
62
f2
91
PAL
/E
91
02 = pAL4 (ryo
o cos w + r
sin
a3 = 1; d3 = P [I
A (r
yo )]/C
= pA (ryo cos
sin
Yn)/C; h = L/n
21
REFERENCES
1.
Duncan, W. J., Ellis, D. L., and Scruton, C.: The Flexural Center and
the Center of Twist of an Elastic Cylinder. Philos. Mag., Series 7,
vol. 16, no. 104, 1933, pp. 201-235.
2.
3.
4.
Wagner, H.: Torsion and Buckling of Open Section. NACA TM-807, 1936.
5.
Garland, C. F.: The Normal Modes of Vibration of Beams having Noncollinear Elastic and Mass Axes. J. Appl. Mech., vol. 7, no. 3, Sept. 1940,
pp. A97-105.
6.
7.
8.
Montoya, J.: Coupled Bending and Torsional Vibrations of a Twisted Rotating Blade. Brown Boveri Rev., vol. 53, no. 3, 1966, pp. 216-230.
9.
Yu, Y. Y.: Variational Equation of Motion for Coupled Flexure and Torsion
of Bars of Thin-walled Open Section Including Thermal Effect. J. Appl.
Mech., vol. 38, June 1971, pp. 502-506.
10. Carnegie, W.: A Note One the Application of the Variational Method to
Derive the Equations of Dynamic Motion of a Pretwisted Cantilever Blade
Mounted on the Periphery of a Rotating Disc Allowing for Shear Deflection,
Rotary Inertia and Torsion Bending. Bull. Mech. Eng. Educ., vol. 5, 1966,
pp. 221-223.
11. Rao, D. K.; and Rao, J. S.: Equations of Motion of Rotating Pretwisted
Blades in Bending-Bending-Torsion with Effects of Warping, Shear, Rotary
Inertia etc. Proc. Silver Jubilee Cong. Aero. Soc. India. Paper No. 4.3,
1973.
12. Subrahmanyam, K. B.; Kulkarni, S. V.; and Rao, J. S.: Application of the
Reissner Method to Derive the Coupled Bending-Torsion Equations of Dynamic
Motion of Rotation Pretwisted Cantilever Blading with Allowance for Shear
Deflection, Rotary Inertia, Warping and Thermal Effects. J. Sound Vib.,
vol. 84, no. 2, 1982, pp. 223-240.
13. Hodges, D. H.; and Dowell, E. H.: Nonlinear Equations of Motion for the
Elastic Bending and Torsion of Twisted Nonuniform Rotor Blades. NASA
TN-D-7818, 1974.
22
14. Rosen, A.; and Friedmann, P. P.: Nonlinear Equations of Equilibrium for
Elastic Helicopter on Wind Turbine Blades Undergoing Moderate Deformation.
UCLA-EN6-7718, University of California, Los Angeles, School of Engineering and Applied Science Report, 1978.
15. Kaza, K. R. V.; and Kvaternik, R. G.: Nonlinear Aeroelastic Equations
for Combined Flapwise Bending, Chordwise Bending, Torsion and Extension
of Twisted Non-Uniform Rotor Blades in Forward Flight. NASA TM-74059,
1977.
16. Rao, J. S.: Turbomachine Blade Vibration. Shock Vib. Dig., vol. 15,
no. 5, 1983, pp. 3-9.
17. Leissa, Arthur: Vibrational Aspects of Rotating Turbomachinery Blades.
Appl. Mech. Rev., vol. 34, no. 5, May 1981, pp. 629-635.
18. Greenwood, Donald T.: The Use of Higher Order Difference Methods in Beam
Vibration Analysis. NASA TN-D-964, 1961.
19. Salvadori, M. G.: Numerical Computation of Buckling Loads by Finite
Differences. Proc. Am. Soc. Civ. Eng. vol. 75, no. 10, Dec. 1949,
pp. 1441-1475.
20. Subrahmanyam, K. B.; and Kaza, K. R. V.: An Improved Finite-Difference
Analysis of Uncoupled Vibrations of Tapered Cantilever Beams. NASA
TM-83495, Sept. 1983.
21. Subrahmanyam, K. B.; and Kaza, K. R. V.: An Improved Finite-Difference
Vibration Analysis of Pretwisted, Tapered Beams. Proceedings, 12th Southeastern Conference on Theoretical and Applied Mechanics, vol. 1, Eng.
Ext. Serv., Sen. Eng., Auburn Univ., 1984, pp. 118-126.
22. Reissner, E.: On a Variational Theorem in Elasticity, J. Math. Phys.
vol. 29, 1950, pp. 90-95.
23. Subrahmanyam, K. B.; and Rao, J. S.: Coupled Bending-Bending Vibrations
of Pretwisted Tapered Cantilever Beams Treated by the Reissner Method.
J. Sound Vib., vol. 82, no. 4, 1982, pp. 577-592.
24. Dean, Thomas Scott; and Plass, H. J. Jr.: A Dynamic Variational Principle
for elastic Bodies, and its Application to Approximations in Vibration
Problems. Developments in Mechanics. Vol. 3, Part 2, T. C. Huang, and
M. W. Johnson, Jr. eds., John Wiley and Sons, Inc., 1967, pp. 107-118.
25. Kumar, D.; Kulkarni, S. V.; and Subrahmanyam, K. B.: Uncoupled Vibrations
of Tapered Cantilever Beams Treated by the Dean and Plass Dynamic Variational Principle. J. Sound Vib., vol. 79, no. 4, 1981, pp. 609-615.
26. Subrahmanyam, K. B.; Kulkarni, S. V.; and Rao, P. M.: Dean and Plass
Method Calculations of Flexural Frequencies of Temeshenko Beams. J. Sound
Vib., vol. 81, no. 1, Mar. 1982, pp. 141-146. . "
23
27. Leissa, A. W.; and Co, Chi-min: Coriolis Effects on the Vibrations of
Rotating Beams and Plates. Proceedings, 12th Southeastern Conference on
Theoretical and Applied Mechanics, vol. 2, Eng.Ext.'Sen/.', Sch. Eng.,
Auburn Univ., 1984, pp. 508-513.
28. Dugundji, J.; and Wendell, J. H.: Some Analysis Methods for Rotating
Systems with Periodic Coefficients. AIAA J., vol. 21, no. 6, 1983,
pp. 890-897.
29. Carnegie, W.: Vibrations of Pretwisted Cantilever Blading. Proc. Inst.
Mech. Eng., Vol. 173, No. 12, 1959, pp. 343-374.
30. Carnegie, W.; Dawson, B.; and Thomas, J.: Vibration Characteristics of
Cantilever Blading. Proc. Inst. Mech. Eng., Vol. 180, pt. 31, 1965-66,
pp. 71-89.
31. Carnegie, W.; and Dawson, B.: Vibration Characteristics of Straight
. Blades of Asymmetrical Aerofoil Cross-Section. Aeronaut. Q., Vol. 20,
1968, pp. 178-190.
32. Rao, J. S.; and Carnegie, W.: Solution of the Equations of Motion of
Coupled Bending-Bending-Torsion Vibrations of Turbine Blades by the Method
of Ritz-Galerkin. Int. 0. Mech. Sci., Vol. 12, 1970, pp. 875-882.
33. Dawson, B.; and Carnegie, W.: Modal Curves of Pretwisted Beams of Rectangular Cross-Section. J. Mech. Eng. Sci., Vol. 11, No. 1, 1969, pp. 1-13.
34. Richardson, L. F.: The Approximate Arithmetical Solution by Finite
Differences of Physical Problems involving Differential Equations, with
an Application to the Stresses in a Masonry Dam. Philos. Trans. R. Soc.
London, Ser. A, vol. 210, no. IX, 1911, pp. 307-357.
35. Rosard, D. D.: Natural Frequencies of Twisted Cantilever Beams. J. Appl.
Mech., vol. 20, no. 2, 1953, pp. 241-244.
36. The International Mathematical and Statistical Library (ISML), Houston,
Texas, 1975.
24
=
=
=
=
E
P
rxo
ryo
=
=
=
=
Pretwist angle,
Y>
deg
0
15
30
45
60
75
89
Torsional rigidity, C
Nm2/rad
lb in.2/rad
9.1265
9.1759
9.3269
9.6282
10.0440
10.6470
11.5600
3180.1684
3197.4
3250.0
3355.0
3500.0
3710.0
4028.0
0*
-30'
FIRST- AND
Mode
number
n o 10
20
30
40
50
60
10
20
30
Extrapolated
1
2
3
4
5
6
7
95.929
581.449
836.642
1074.560
1557.277
2875.489
3064.882
96.534
599.808
841.330
1076.166
1659.921
3087.379
3201.482
96.648
603.358
842.204
1076.468
1680.532
3091.582
3269.552
96.687
604.613
842.510
1076.574
1687.877
3093.051
3294.057
96.706
605.196
842.652
1076.623
1691.301
3093.731
3305.523
96.716
6D5.513
842.729
1076.650
1693.168
3094.100
3311.786
96.706
603.475
842.658
1076.654
1672.018
3093.880
3207.516
96.734
605.911
842.874
1076.704
1694.623
3094.862
3313.419
96.737
606.141
842.895
1076.709
1696.647
3094.923
3322.705
96.738
606.204
842.901
1076.710
1697.190
3094.937
3325.192
1
2
3
4
5
103.644
384.472
1032.234
1187.743
1187.743
2622.674
3112.428
98.249
482.100
988.711
1109.888
1548.318
3089.437
3151.108
98.552
501.378
980.034
1104.447
1607.632
3129.428
3213.081
97.324
508.273
976.911
1102.855
1628.241
3133.785
3245.109
97.221
511.493
975.446
1102.162
1637.768
3135.325
3260.536
97.165
513.251
974.645
1101.796
1642.942
3136.082
3269.030
101.665
464.308
1011.948
1134.652
1490.641
3074.468
3136.994
98.171
507.565
981.844
1105.379
1632.749
3139.478
3272.776
97.554
513.490
976.653
1102.703
1647.241
3138.599
3286.343
97.378
515.149
975.145
1101.983
1651.192
3138.141
3290.033
1
1
2
3
4
5
6
7
94.042
822.909
1079.098
96.664
96.738
96.739
573.170 599.444 606.259
842.904
842.363 842.903
1076.708 1076.708 1076.711
3183.905 1596.621 1629.061
5216.526 3107.792 3095.274
5293.266 3210.000
5
96.739
606.219
842.904
1076.711
1697.967
3051.492
3095.062
96.739
606.235
842.904
1076.711
1696.414
3094.935
3328.162
= -30*
1
1
2
3
4
5
6
7
94.514
623.304
1072.483
97.041
96.997
97.041
515.031 517.281
495.286
962.517 972.845
960.023
1097.665 1096.317 1100.997
3180.700 1565.746 1600.834
3693.199 3159.649 3066.255
5559.758 3187.625
5
97.041
517.262
972.777
1100.970
1655.388
3031.560
3154.065
6
97.041
517.267
972.809
1100.990
1653.893
3137.568
3291.277
Mode
number,
classification
1
2
3
4
5
6
7
8
9
10
(Fl)
(F2)
(SI)
(Tl)
(F3)
(T2)
(F4)
(T3)
(S2)
(F5)
Uncoupled
bending and
torsion modes
r
xo = ryo = Y =
96.738
606.158
864.603
1051.080
1696.742
3153.229
3323.060
5255.263
5417.605
5488.150
Y = 45
Y = 60
Y = 89
98.417
454.448
1136.247
a
!079.585
1600.250
a
3238.745
3234.039
a
5397.786
5149.207
6086.129
99.319
403.476
1250.120
a
l!02.668
1562.417
a
3307.992
3176.773
a
5513.196
5047.868
6417.764
100.255
330.908
1282.718
a
l!82.920
1674.196
a
3548.748
3071.047
a
5914.446
4815.582
7315.336
= 30
97.558
513.908
1014.142
a
!062.557
1648.125
a
3187.661
3281.225
a
5312.648
5242.015
5815.830
Mode
number
1
2
3
4
5
6
7
Uncoupled
frequency
rxo = ryo = Y = 0
96.738
606.158
864.603
1051.080
1696.742
3153.229
3323.060
(Fl)
(F2)
(SI)
(Tl)
(F3)
(T2)
(F4)
Coupled
bendingtorsion
frequency,
Coupled bending-bending-torsion
frequency
Y = 0
Y = 30
Y = 45
Y = 60
Y = 89
96.737
606.142
a
864.603
1051.095
1696.647
3153.270
3322.728
97.558
513.891
1013.977
1062.743
1648.040
3187.303
3281.327
98.417
454.432
1136.514
1079.348
1600.157
3245.592
3226.973
99.319
403.461
1250.296
1102.544
1562.299
3309.104
3175.471
100.255
330.898
1283.036
1182.623
1674.142
3549.324
3070.251
1
2
3
4
5
6
7
Uncoupled
frequency...^.,.;.
r
xo = ryo = Y = 0
96.738 (Fl)
606.158 (F2)
864.603 (SI)
1051.080 (Tl)
1696.742 (F3)
3153.229 (T2)
3323.060 (F4)~
Coupled
bend ingtorsion
frequency,
Y =0
a
96.738
606.158
842.893
1076.696
a
!696.742
3094.865
a
3323.060
a
Coupled bending-bending-torsion
frequency
= 30
Y = 45
Y = 60
Y = 89
97.556
513.342
974.417
1105.563
1647.418
3131.506
3290.388
98.412
454.127
1170.423
1047.110
1599.494
3270.061
3166.320
99.311
403.354
1264.208
1088.446
1561.908
3304.629
3154.937
100.246
330.907
1289.706
1174.069
1671.711
3542.868
3069.589
1
2
3
4
5
6
7
Uncoupled
frequency
r
xo = ryo = Y = 0
Y = 0
Y = 30
Y = 45
Y = 60
Y = 89
96.738 (Fl)
606.158 (F2)
864.603 (SI)
1051.080 (Tl)
1696.742 (F3)
3153.229 (T2
3323.060 (F4)
96.737
606.141
842.895
1076.709
1696.647
3094.922
3322.705
97.557
513.162
972.178
1108.465
1647.424
3124.353
3294.018
98.413
453.978
1174.357
1043.971
1599.445
3277.243
3154.608
99.313
403.259
1267.441
1085.937
1562.017
3308.751
3145.672
100.246
330.896
1292.688
1170.910
1672.784
3539.662
3067.755
10
20
30
40
50
60
831
550
736
484
966
838
856
5390
17014
1
2
3
4
5
6
59
215
490
924
1494
2193
84
348
847
1639
2881
4410
Figure L - Cross-section of asymmetrical airfoil blade showing position of the centroid and center-of-flexure.
1300 r
1200
1100
1000
ANTICLOCKWISE PRETWIST
(POSITIVE?)
- CLOCKWISE PRETWIST
(NEGATIVE y)
TRANSFORMATION METHOD
(NEGATIVE y>
EXPERIMENTAL (NEGATIVE?)
ALSO FOR
ANTICLOCKWISE
PRETWIST-v
>o
90
(a) First, second, third, and fourth coupled mode frequencies.
3600
Of
=>
3500
3400
1500
20
40
60
1.0
RELATIVE
AMPLITUDE
.8
.6
.4
.2
0
-.2
-.4
-.6
-.8
-1.0
.4
.6
.8
1.0
.2
.4
.6
.8
1.0
Figure 3. - Coupled bending-bending, and uncoupled torsion mode shapes of a pretwisted blade of
doubly symmetric cross section. rxo ryo 0, y 89.
-1.
1.0
.8
.6
.4
Q-
.2
0
-.2
-.4
-.6
-.8
-1.0
-1.0
.2
.4
.6
.8
AXIAL FRACTIONAL LENGTH
1.0
RELATIVE
AMPLITUDE
-1.0
UJ
a:
-1.0
.2
.4
.6
.8
1.0
0
.2
AXIAL FRACTIONAL LENGTH
.4
.6
.8
1.0
RELATIVE
AMPLITUDE
.8
.6
.4
o.
.2
0
-.2
-.4
-.6
-.8
-1.0
-.6
-.8
-1.0
.2
.4
.6
.8
AXIAL FRACTIONAL LENGTH
1.0
1. Report No.
NASA TM-83735
5. Report Date
505-33-42
B. Performing Organization Report No.
7. Author(s)
E-2175
Technical Memorandum
K. B. Subrahmanyam, NBKR Institute of Science and Technology, Mechanical Engineering Dept., Vidyanagar 524 413, India and NRC-NASA Research Associate; K. R.
V. Kaza, Lewis Research Center. Prepared for the XVIth International Congress of
Theoretical and Applied Mechanics cosponsored by the International Union of Theoretical and Applied Mechanics (IUTAM) and the Technical University of Denmark,
Lyngby, Denmark, August 19-25, 1984.
16. Abstract
Blade vibration
Reissner method
Dean and Plass method
Finite- difference methods
19. Security Classlf. (of this report)
Unclassified
Unclassified - unlimited
STAR Category 39
Unclassified
*For sale by the National Technical Information Service. Springfield, Virginia 22161
22. Price*
Washington, D.C.
20546
Official Butineu
Pomgt md FMS Ptd
National Aeronautics and
Space AdmtniMrttion
NASA-451
fVIASA
POSTMASTF.R:
U Undeliverihlc (Srnion I S*
PollI Mtnutl) lt Nut Rrlurn