Tools › Aerodynamics › Rotors and Propellers › Fan Blade Stress


Fan Blade Stress and Frequencies

Centrifugal force, flapwise and lagwise bending moments, combined stress and the first bending frequencies of a ducted-fan blade from its planform, section, material and operating point. The aerodynamic loads along the span come from the same blade-element momentum model as the Ducted Fan Performance tool; the blade sections are modeled as solid or thin-shell ellipses of the entered chord and thickness (their area and second moments by the site's section-geometry library) and pitched with the local blade angle, so the thrust and in-plane moments resolve into flap and lag bending about the section axes. The tool reports the root loads, the stress along the span with the critical station, margins against yield and ultimate with a factor of safety, the blade mass, and the first flap and lag frequencies at rest and at speed (Rayleigh with Southwell's centrifugal stiffening) on a Campbell diagram against the per-rev harmonics.

  • Spanwise thrust and in-plane loads from blade-element momentum theory on the entered blade

  • Centrifugal force, flap and lag bending moments and the combined stress at every station; the critical station and its margins

  • Section as a solid or thin-shell ellipse of the local chord and thickness, pitched with the blade angle

  • Blade mass, first flap and lag bending frequencies at rest and at speed, Campbell diagram with 1P to 6P lines

Fan Blade Stress and Frequencies


Unit System: /

Rotor and Operating Point







Blade Stations (hub to tip)

1
2
3
4
5

Section and Material










Airfoil




















Resolution


×

Fan Blade Stress and Frequencies

Blade Loads, Stress and Frequencies

Loads

A fan blade is a rotating cantilever. Three loads act on it: the centrifugal force of its own mass, which pulls along the span and grows from zero at the tip to its full value at the root; the thrust distribution, which bends the blade out of the rotor plane (flap); and the in-plane force that reacts the torque, which bends it in the plane of rotation (lag). The thrust and in-plane loads along the span come from the same blade-element momentum model as the Ducted Fan Performance tool, so the stress analysis uses the loads the aerodynamic design implies. The centrifugal force at a station is the integral of ρm A Ω² r over the blade outboard of it; the bending moments are the integrals of the outboard loads times their arm.

Sections

The section at each station is modeled as an ellipse whose axes are the local chord and thickness, either solid or a thin shell of constant skin thickness. Its area and second moments are computed exactly by the site's section-geometry library (Green's theorem on the polygon), the same code that serves the Section Properties and Beam Stress tools. An ellipse has about 20 percent more area than a real airfoil of the same chord and thickness and a comparable flapwise second moment, so the stresses are a first estimate; for a detailed check, draw the real section in the Beam Stress tool with the loads reported here.

The section sits at the local pitch θ, so the flap moment (vector in the rotor plane) and the lag moment (vector along the axis) are resolved into moments about the chord and about the thickness axis before the bending stress is formed. The stress σ = CF/A + Mu v / Iu − Mv u / Iv is evaluated around the ellipse and its largest tensile value taken; the station with the largest value is the critical one, usually near the root for a tapered blade but sometimes outboard where the thickness falls faster than the moment.

Margins

The operating-point loads are limit loads. The yield margin is written at limit load times the yield factor FSyld (default 1.15, the value several certification agencies require for no detrimental permanent set; 1.0 for a plain limit check), the ultimate margins at limit load times the factor of safety FS (default 1.5), following aircraft practice. The overspeed case multiplies the centrifugal stress by the square of the overspeed factor while holding the aerodynamic bending: a conservative bound on the stress of a runaway rotor, since the aerodynamic loads would in fact also rise.

Frequencies

The first flap and lag bending frequencies are found by Rayleigh's method with the cantilever shape ψ = 1 − cos(π(r − Rh)/(2L)), which satisfies the clamped root. At rest, ω0² = ∫EI ψ″² dr / ∫m ψ² dr with the flapwise or lagwise second moment. Rotation stiffens the blade through the centrifugal tension: Southwell's relation gives ωflap² = ω0² + K Ω² with K = ∫CF ψ′² dr / (Ω² ∫m ψ² dr), and for the lag mode the in-plane centrifugal softening reduces the coefficient by one. Rayleigh's frequencies are upper bounds, typically within a few percent for the first mode.

The Campbell diagram plots both frequencies against rotor speed together with the per-rev lines 1P to 6P. A crossing near the operating speed is a resonance: 1P is excited by any inflow asymmetry (a crosswind, an inlet distortion), and the duct sees the blade-passing order. The tool flags a first mode within 10 percent of an integer per-rev at the operating speed.

Not included

Torsion and the pitch-flap coupling of a swept or skewed blade, the root attachment and its stress concentration, fatigue (the alternating 1P bending of a fan in a crosswind), composite ply-level allowables, and the aerodynamic damping of the modes.




References

Johnson, W. (1980). Helicopter Theory. Princeton University Press, Princeton, NJ.

Bielawa, R. L. (2006). Rotary Wing Structural Dynamics and Aeroelasticity, 2nd ed. American Institute of Aeronautics and Astronautics, Reston, VA.

Southwell, R. V., and Gough, B. S. (1921). The Free Transverse Vibration of Airscrew Blades. Reports and Memoranda No. 766, Aeronautical Research Committee, London.

Glauert, H. (1935). "Airplane Propellers." In Durand, W. F. (ed.), Aerodynamic Theory, Vol. IV, Division L, Julius Springer, Berlin, pp. 169–360.

Bruhn, E. F. (1973). Analysis and Design of Flight Vehicle Structures. Jacobs Publishing, Carmel, IN.


DISCLAIMER: The methods and calculation on aNOWysis.com are provided as means to educate the reader. There is no warranty, expessed or implied, as to the relevance or accuracy of the methods or calculations. All calculations should be checked by a qualified analyst for appropriateness and accuracy. We will not be held liable for any damages, injuries or deaths caused by using methods or calculations from this site.