Abstract

Wind turbine blades are adjusted in real-time according to the wind conditions and blade
deformations to improve power generation efficiency. It is necessary to predict and reduce the
aeroelastic deformations of wind turbine blades. In this paper, the equivalent model of the blade is
established by the finite element method (FEM), and the aerodynamic load of the blade is evaluated
based on the blade element momentum (BEM) theory. The aeroelastic coupling model is established
in which the bending-torsion coupling effect of the blade is taken into account. The steady and
dynamic aeroelastic deformations are calculated. The influences of the blade section's shear centre
position and the blade's sweepback design on the deformations are analyzed. The novel approaches
reducing the twist angle of the blade by changing the shear centre position and sweepback of the
blade are presented and proven to be feasible.

Dynamic Analysis of Horizontal Axis

Wind Turbines under Thunderstorm
Downbursts

Abstract
Wind turbines are among the fastest-growing technologies for producing sustainable
energy. Although numerous studies have been carried out to investigate the dynamic
behavior of wind turbines exposed to synoptic winds, the dynamic behavior of such
structures under downbursts is not comprehensively studied and not yet well-
understood. As such, the main objective of the current study is to develop a numerical
model that can be used to investigate the dynamic response of wind turbines subjected
to downbursts taking into account the fluid-structure interaction (FSI) effect. To
achieve this task, two modules are developed in this study and are integrated with the
open-source code, fatigue, aerodynamics, structures, and turbulence (FAST). These
modules introduce a three-dimensional (3D) time history of the downburst wind field
into FAST inflow-wind module and simulate the interaction between the downburst
wind field and the wind turbine structure. The downburst wind field consists of a
moving mean component with a superimposed turbulence. The mean component is
generated using a previously conducted computational fluid dynamics (CFD)
simulation. The turbulence is simulated using the consistent discrete random flow
generation (CDRFG) method as a stochastic process based on the turbulence power
spectral density and the coherence functions pertaining to downbursts. The developed
numerical model is validated using the results of a previously performed experiment
on a wind turbine model sited in microburst-liked winds. A parametric study is then
conducted on a wind turbine model to investigate the dynamic responses of the tower
and blades with and without FSI effect under various downburst configurations. A
comparison between the wind turbine quasistatic and dynamic responses is conducted,
and the dynamic amplification factor (DAF) is then calculated. The aerodynamic
damping of the blades is also estimated and compared with a previously developed
closed form solution.
 Modern horizontal axis wind turbine blades are long, slender, and flexible

structures that can undergo considerable deformation, leading to blade failures

For this reason, it is important to estimate blade .)e.g., blade-tower collision(

behaviors accurately when designing large-scale wind turbines. In this study, a

numerical analysis considering blade torsional degree of freedom, geometric

nonlinearity, and gravity was utilized to examine the effects of these factors on the

aeroelastic blade behavior of a large-scale horizontal axis wind turbine. The results

predicted that flapwise deflection is mainly affected by the torsional degree of

freedom, which causes the blade bending deflections to couple to torsional

deformation, thereby varying the aerodynamic loads through changes in the

effective angle of attack. Edgewise deflection and torsional deformation are

mostly influenced by the periodic gravitational force on the wind turbine blade

Les pales d'éoliennes modernes à axe horizontal sont longues, minces et flexibles

structures pouvant subir des déformations considérables, conduisant à des ruptures de pales

(par exemple, collision lame-tour). Pour cette raison, il est important d'estimer la lame

comportements précis lors de la conception d’éoliennes à grande échelle. Dans cette étude,
un

analyse numérique considérant le degré de liberté de torsion de la pale, géométrique

la non-linéarité et la gravité ont été utilisées pour examiner les effets de ces facteurs sur le

comportement aéroélastique des pales d'une éolienne à axe horizontal à grande échelle. Les
résultats

prédit que la déflexion des volets est principalement affectée par le degré de torsion de

liberté, ce qui amène les déflexions de flexion de la lame à se coupler à la torsion

déformation, faisant ainsi varier les charges aérodynamiques à travers des changements
dans la

angle d'attaque efficace. La déflexion sur chant et la déformation en torsion sont

principalement influencé par la force gravitationnelle périodique sur la pale de l'éolienne

Cet article présente une analyse de la dynamique de la pale flexible d'une HAWT (éolienne à
axe horizontal) à grande échelle. Pour modéliser plus précisément la flexibilité de la pale, le
model «non linéaires d'Euler-Bernoulli est utilisée pour la modélisation structurelle.) La
vibration de flexion de la pale et la déflexion de torsion sont toutes deux prises en compte
lors du calcul des charges aérodynamiques ; ainsi, la théorie BEM (blade element
momentum) utilisée dans cette étude est modifiée. De plus, le modèle de décrochage
dynamique BeL (BeddoeseLeishman) est intégré au modèle modifié par BEM pour étudier les
caractéristiques de décrochage dynamique du profil aérodynamique.

La combinaison de ces deux types de formulations(structurelles et aerodynamique) donne

les équations aéroélastiques intégro-différentielles non linéaires. les équations gouvernantes
seront résolues en utilisant une méthode numérique( methode of line) sous mathematica
pour prédire la réponse aéroélastique du problème.

Les réponses dynamiques dans le domaine temporel de la pale d’éolienne offshore de 5 MW

du NREL (National Renewable Energy Laboratory) sont obtenues. Les résultats de la
simulation indiquent que la vibration et la déformation de la pale ont des effets significatifs
sur les charges aérodynamiques, et que le décrochage dynamique peut provoquer une
fluctuation plus violente des charges aérodynamiques de la pale par rapport au modèle
aérodynamique stable, ce qui peut affecter considérablement l'analyse du spectre de charge
de fatigue de la pale et le conception de la durée de vie en fatigue.

Introduction
À mesure que la capacité de l'éolienne augmente, La taille et la flexibilité des composantes
tel que les pales continuent d'augmenter, et les pales subissent souvent d'importantes
déformations et vibrations, des problèmes géométriques non linéaires évidents se
produisent.

Les pales subissent un couplage aéroélastique dans lequel les forces aérodynamiques et les
réponses dynamiques des composants de l'unité sont couplées les unes aux autres.

Le problème du couplage devient de plus en plus important, la vibration des pales de

l'éolienne et la vibration de l'accouplement dynamique de la tour et de l'arbre principal
affecteront la réponse aéroélastique de l'unité, affectant ainsi la durée de vie en fatigue et la
stabilité aéroélastique de l'unité.

La vibration de la pale la fait bouger Lors d'un déplacement avec un angle d'attaque élevé ou
même un décrochage profond, la bombe aérienne de l'équipage L'analyse nécessite une
analyse précise des caractéristiques instables telles que le décrochage dynamique du profil
de la pale. simulation précise [ 1].

Par conséquent, il est nécessaire de développer un couplage aéroélastique efficace et fiable

pour l'unité. Un modèle d'analyse complet est utilisé pour analyser avec précision la charge
aéroélastique de la pale. Pour la conception optimisée des éoliennes à grande échelle, les
critères de durée de vie en fatigue et de stabilité de fonctionnement de l'unité sont d'une
grande importance [2].
L'effet de couplage aéroélastique des éoliennes est causé par la pression de l'air subie par
l'unité.

Le mode aéroélastique est provoqué par le couplage mutuel de la charge dynamique, de la

déformation structurelle et des vibrations.

Le modèle comprend le calcul de la charge aérodynamique des pales et l'analyse de la

réponse dynamique structurelle. exister

Pièce de modèle aérodynamique, théorie du moment des éléments de pale

momentum, BEM) les résultats des calculs sont plus raisonnables et la solution est rapide.

Dans les éoliennes grand public telles que HAWC2, FAST et GH-Bladed,

Il est largement utilisé dans les logiciels d’analyse élastique [1]. Dans des conditions
instables, la pale L'angle d'attaque du profil aérodynamique varie dans le temps et les
caractéristiques aérodynamiques du profil aérodynamique changent avec l'angle d'attaque.

La courbe aura des caractéristiques d'hystérésis de décrochage dynamique [3], et le

coefficient aérodynamique peut Grands écarts par rapport aux valeurs statiques, comme par
exemple l'éolienne d'essai conjointe conçue par NREL À des vitesses de vent élevées, il existe
un écart de 15 à 20 % entre la valeur de conception et la valeur réelle.

La différence [4] est de réduire la production d'énergie d'environ 20 %. Le retard de la charge

aérodynamique L'hystérésis provoquera à son tour une plus grande vibration des pales, ce
qui entraînera leur durée de vie en fatigue.

raccourcir. Les recherches de SHIPLEY et al. [5] montre que le décrochage dynamique peut
causer relativement Grandes fluctuations de charge aérodynamique. En raison de la
complexité des flux, il est actuellement Théoriquement, il n’existe aucun modèle permettant
de déterminer avec précision la réponse dynamique instable d’un profil aérodynamique.

En pratique, des modèles de calcul semi-empiriques basés sur des résumés de tests sont
souvent utilisés [6-7]. Y compris le modèle ONERA Edlin[8], le modèle à ailes fines
THEODERSON modèle [9] et modèle BEDDOES-LEISHMAN [10-12], etc.

Le modèle B-L a été valorisé et reconnu par les chercheurs [13], Liu Xiong et al. [3, 6, 12]

Ce modèle a été utilisé pour étudier les caractéristiques de décrochage dynamique du profil
aérodynamique et révéler montre ses caractéristiques instables.

Pour le modèle dynamique du système d'éolienne, la construction actuellement

couramment utilisée Les méthodes du module incluent : Méthode de forme de mode
supposé (Forme de mode supposée méthode, AMSM), méthode des éléments finis
(Méthode des éléments finis,FEM) et méthode du système multi-corps (méthode du système
multi-corps,MBS), etc. [14]. Les méthodes MBS modernes prennent en compte la flexibilité
des composants du système Il convient à la construction mécanique de systèmes complexes
comprenant des corps rigides et des corps flexibles.

modèle[15]. MURTAGH et coll. [16-17] énergie éolienne établie grâce à la méthode MBS
L'équation dynamique de l'ensemble du système machine, combinée au modèle
aérodynamique BEM, L'influence de la rotation de l'éolienne sur la vibration de la tour et la
dynamique de l'ensemble de la machine sont analysées numériquement.

force de réponse. Pour les composants avec mouvement spatial et déformation élastique,

ZHAO et coll. [18-19] ont introduit des « super unités » et les ont discrétisées en plusieurs

Un corps rigide relié par des paires cinématiques, des ressorts et des amortisseurs, adapté au
vent Modélisation mécanique de composants tels que les pales de moteur, les tours et les
arbres principaux, tandis que la nacelle,

Les composants tels que les moyeux de roue sont directement traités comme des corps
rigides et grâce à la théorie de la dynamique multi-corps

La dynamique non linéaire du système peut être établie avec moins de degrés de liberté.

Procédure. Li Deyuan et coll. [20] ont utilisé la théorie informatique de la dynamique des
systèmes multi-corps et Méthode de modélisation Roberson-Wittenburg, combinée à
l'aérodynamique BEM modèle, l'équation de couplage aéroélastique multi-corps rigide-
flexible de l'ensemble de la machine est établie, et l'ensemble de la machine est Les modes
de vibration et la réponse dans le domaine temporel du couplage aéroélastique de la
machine ont été analysés numériquement.

simulation.

Compte tenu de l'effet de couplage dynamique de la charge aérodynamique, ce papier)

Introduction de rezai : Les développements et progrès récents dans la technologie de

l’énergie éolienne ont conduit à une croissance constante de la capacité et de la taille des
éoliennes modernes. Cette tendance devrait se poursuivre à l’avenir, notamment dans les
applications offshore [1]. À cet égard, des recherches universitaires et des efforts industriels
sont en cours pour développer une éolienne de 10 MW avec des pales de rotor de 140 m de
diamètre [2,3]. Néanmoins, la question de la limitation du poids est devenue un défi majeur
pour l'agrandissement des éoliennes et a obligé les ingénieurs à concevoir des structures
plus légères qui répondent aux exigences de résistance obligatoires. Ces contraintes ont
conduit à la conception de structures d'éoliennes plus flexibles et, notamment, de pales
d'éoliennes, conçues pour un chargement dynamique. Cette plus grande flexibilité
structurelle pourrait entraîner d'importantes déformations sous charge aérodynamique,
même dans des conditions nominales de fonctionnement [4,5]. Dans ces cas, bien que les
déformations restent faibles, des non-linéarités spéciales (appelées non-linéarités
géométriques) sont introduites dans le système par les grandes rotations résultant des
déformations structurelles élastiques [6]. Des travaux de recherche récents sur le
comportement aéroélastique des éoliennes [4,5,7e10] ont démontré que les non-linéarités
géométriques affectent tous les aspects de ces structures, y compris leurs caractéristiques de
rigidité et de masse ainsi que les chargements aérodynamiques.

Dans la plupart des travaux de recherche antérieurs, les aubes de turbine ont été analysées
structurellement au moyen de modèles de poutres qui fournissent des résultats
suffisamment précis pour un coût de calcul raisonnable. Les résultats analytiques peuvent
être améliorés en utilisant des modèles de poutre d'ordre supérieur tels que VABS [11] ou
des modèles de poutre développés sur la base des théories avancées des composites à
parois minces [12]. Cependant, comme les non-linéarités géométriques affectent
principalement la dynamique globale d'un système, la majorité des travaux de recherche
antérieurs ont utilisé des modèles de poutres d'Euler-Bernoulli non linéaires pour éviter la
complexité supplémentaire généralement introduite par la prise en compte des effets des
non-linéarités géométriques. En utilisant un modèle d'ordre réduit (ROM) dans les directions
volet et bord, Larsen et Nielsen [7] ont étudié les effets des non-linéarités structurelles pour
une pale d'éolienne en rotation soumise à des forces aérodynamiques constantes. Ils ont
ensuite utilisé leur modèle développé pour explorer le rôle du mouvement harmonique du
pied de pale dans les caractéristiques de stabilité des pales en rotation [8]. Wu et coll. [13]
ont utilisé une formulation de poutre géométriquement exacte pour étudier les grandes
déformations d'une pale en rotation au moyen d'un modèle par éléments finis (FEM). Rezaei
et coll. [4] ont utilisé une formule similaire pour modéliser une pale d'éolienne de 5 MW-
NREL et ont étudié les effets des non-linéarités géométriques pour différentes déformations
et conditions de charge. Le modèle développé a été étendu aux analyses aérodynamiques
quasi-stationnaires et instationnaires en tenant compte des non-linéarités associées aux
charges aérodynamiques [5]. Les investigations ont montré l'importance de considérer la
torsion induite par les grandes déformations des pales dans la direction latérale et les effets
de ces déformations sur l'instabilité du flottement des pales en rotation.

Bien que les effets cruciaux des non-linéarités géométriques sur les caractéristiques
dynamiques et aéroélastiques des éoliennes modernes aient été bien démontrés dans des
travaux récents, la plupart des résultats rapportés sont basés sur des modèles monopales qui
ignorent complètement les flexibilités de la tour de l'éolienne et des autres pales. . Par
conséquent, il est essentiel d’étudier les effets de la flexibilité de la tour et les interactions
avec la dynamique des pales du rotor sur le comportement aéroélastique de systèmes
d’éoliennes entiers. Lee et coll. [10] ont présenté un modèle multi-flexible d'une éolienne à
deux pales et ont utilisé le FEM pour obtenir les valeurs propres du système à charge libre.
Dans un autre travail [14], la méthode de synthèse des modes composants et un modèle
linéaire d'éolienne ont été utilisés pour obtenir les réponses dynamiques transitoires.
Thomsen et coll. [15] ont présenté une nouvelle technique expérimentale pour mesurer
l'amortissement aérodynamique des vibrations des pales sur le bord et ont démontré que
l'amortissement total du système est significativement affecté par le mouvement
tourbillonnant du rotor. Ils ont conclu qu’il est essentiel d’analyser l’ensemble d’un système
éolien plutôt qu’une seule pale, en particulier lorsqu’on tente de prédire la stabilité de
l’éolienne. En utilisant un modèle d'éolienne linéaire à trois pales, Hansen et al. [16] ont
dérivé les formes invariantes dans le temps des équations gouvernantes en appliquant la
méthode de transformation de coordonnées multi-lames (également appelée transformation
de Coleman). Ils ont analysé le comportement dynamique d'une pale dans les modes de
tourbillonnement avant et arrière. Bir [17] a utilisé la même technique pour étudier les
tendances de fréquence et d’amortissement d’une pale d’éolienne de 5 MW-NREL à
différentes vitesses de rotation. Dans un autre travail [18] sur un rotor d'éolienne à pales
anisotropes (un rotor à pales différentes pour lequel la transformation de Coleman ne
pouvait pas être appliquée), la méthode plus générale de transformation de Lyapunov-
Floquet a été utilisée pour obtenir l'invariance temporelle

1. STRUCTURAL MODEL
1.3 Blades: Blade in-plane (lead-lag) vibrations defined about zbi axis of the rotating
coordinate system placed at the root of the blade are denoted by θi , out-of-
 plane (flap) vibrations defined about ybi axis are denoted by ϕi and torsional
vibrations defined about xbi axis are denoted by i ε (refer Fig. 2(a)), where i =
1,2,3 denote the blade number. The kinetic energy of the blade is obtained by
integrating the kinetic energy of a small element of length dr located at a radial
distance of r from the blade root, over the full length of the blade (lb). Position
vector of the small element dr is given below.
 (5) where ( ) ( ) ( ) Z i Y i TX i T θ ,T ϕ , ε are the transformation matrices given
in the Appendix; (r,cgy,cgz) are the coordinates of the C.G. of the blade section in
local coordinate system xbiybizbi. Energy expressions for the i th blade are given
below
2. 2. AERODYNAMIC MODEL Beam element momentum (BEM) theory [7] is used to
calculate aerodynamic loads considering wind shear and tower shadow effects. Wind
shear is modeled using wind profile power law using mean wind velocity at nacelle
height and a parameter of 0.2 for the amount of shear. Tower shadow is modeled
using potential theory of flow around a cylinder [7]. The tower diameter variation
along its height is also taken into account in the calculation of the wind velocity near
the tower. Prandtl’s tip loss factor and Glauert corrections are considered in the BEM
theory used in this study. Aerodynamic loads are calculated at the aerodynamic
center of the aerofoil sections of the blade. Velocity triangle without and with
considering blade vibrations are shown in Fig. 3. Figure 3. Velocity triangle at the
blade section (a) without and (b) with considering blade vibrations Expressions for
the inflow angle, angle of attack without considering blade vibrations are given in Eq.
(13) & (14). ( ) ( )        + − Θ = − r a' V a tan w if 1 1 1 0 ω (13) Θaoa = Θif − Θp+t
0 0 (14) where, 0 0 if aoa Θ ,Θ are the inflow angle and angle of attack, ignoring the
blade vibrations; Θ p+t is the sum of blade pitch and section twist angles; Vw and ω
are the wind velocity and rotational frequency of the blade; r is the radial distance of
the blade section from the hub center; a and a’ are the axial and tangential induction
factors. Blade vibration velocities change the relative velocity of wind entering the
blade section as shown in Fig. 3(b). As the structural model of the blade is built with
respect to pitch axis, due to the offset between aerodynamic center and pitch axis,
both the inflow angle and the angle of attack (AOA) depend on torsional vibrations.
Expressions for the inflow angle, angle of attack considering blade vibrations are
given in Eq. (15) & (16). ( ) ( )         + + − + Θ = − 2 1 1 1 1 r a' V V a V tan w if ω
(15) where, ( ) ( ) i AC p t i i i AC p t i i V r a sin V r a cos θ ε ε ϕ ε ε     = − Θ + =
− + Θ + + + 2 1 ( ) aoa if p t i Θ = Θ − Θ + ε + (16) where, if aoa Θ ,Θ are the inflow
angle and angle of attack considering blade vibrations; i i i θ ,ϕ ,ε  are the
velocities of the i th blade in-plane, out-of-plane and torsional vibrations; i ε is the i
th blade torsional vibration; aAC is the distance between the aerodynamic center
(A.C.) and pitch axis (P.A.) (shown in Fig. 3) of the blade section. Change in the AOA
due to blade vibrations can be approximated according to Eq. (17). ( ) ( ) i rel w rel
aoa aoa V V V a V V r a' ε ω − − − + ≈ ∆Θ = Θ − Θ 2 2 0 2 1 0 0 1 1 (17) Lift, drag and
pitching moment coefficients at Θaoa can be expanded using Taylor series expansion
about the 0 Θaoa as given in Eq. (18). ( ) ( ) ( ) ( ) ( ) ( ) (Θ ) ≈ (Θ )+ (Θ )∆Θ Θ ≈ Θ + Θ
∆Θ Θ ≈ Θ + Θ ∆Θ 0 0 0 0 0 0 M aoa M aoa M aoa D aoa D aoa D aoa L aoa L aoa L aoa
C C C ' C C C ' C C C ' (18) where, CL, CD, CM are the lift, drag and moment
coefficients of the aerofoil; CL’, CD’, CM’ are the slopes of the lift, drag and moment
coefficient curves. Aerodynamic loads considering these coefficients are expressed in
Eq. (19)-(21) where, if written in matrix form the coefficients of the variables i i i i
 θ ,ϕ ,ε ,ε  can be separated into matrices CAero and KAero which are known as
aerodynamic damping and stiffness matrices. ( ) L ( aoa ){ ( ) w( ) rel

3. RESULTS AND DISCUSSION Data of the model wind turbine [8] are inserted into the system
of Eq. (12) and eigenvalues are obtained at zero speed and by ignoring gravity terms. Natural
frequencies of the present 13 DOF model are compared with the values predicted using FAST
(aeroelastic simulation tool) [8] in Table 1. As FAST doesn’t model blade torsional vibrations
its natural frequency in the current rigid beam model is compared with the frequency value
predicted in [5] where blade is modeled using FEM. Natural frequencies calculated using
these two models (FEM vs. 13 DOF model) at zero speed match well and thus this 13 DOF
model can be considered as a good representation of dynamic behavior of the wind turbine
structure in the low frequency vibration modes. NREL 5 MW wind turbine operates over a
wind velocity range of 3-25 m/s and its rotational speed increases from 7.35 rpm at a wind
velocity of 3 m/s to 12.1 rpm at a wind velocity of 11.4 m/s and operates at this rotational
speed till 25 m/s wind velocity by changing pitch angle of the blades [8]. Natural frequencies
of the 13 DOF model considering only the rotational speed variation with wind velocities
between 3-25 m/s are shown in Fig. 4. Natural frequencies of the blade DOF change till 11.4
m/s wind velocity due to increase in rotational speed and thereafter remains constant as
there is no change in the rotational speed. Aerodynamic loads calculated for mean wind
velocities in the range between 3-25 m/s at the nacelle are used to determine aeroelastic
natural frequencies of this model which are plotted in Fig. 5. Differences in the natural
frequencies in Fig. 4 & 5 are only due to the aerodynamic coupling considered in the latter
case. Aeroelastic natural frequencies calculated using the FEM model of this wind turbine are
reported in Ref. [10] which are shown in Fig. 6. All vibration modes except tower bending
modes change with an increase in wind velocity. Lag modes in Fig. 5 & 6 follow the same
trend with increasing wind velocities. Flap modes decrease with an increase in wind velocity
in Fig. 6, whereas in Fig. 5, only the collective and regressive flap modes are decreasing with
Table 1. Comparison of natural frequencies Mode Natural frequency (Hz) FEM 13 DOF Model
Tower fore-aft (α) 0.3240 * 0.3237 Tower side-to-side (β) 0.3120 * 0.3291 Yaw motion (γ) -
6.8346 Drivetrain torsion (ψ) 0.6205 * 0.6099 Blade regressive lead-lag (θs) 1.0793 * 1.1170
Blade progressive lead-lag (θc) 1.0898 * 1.1399 Blade collective lead-lag (θ0) - 3.7978 Blade
regressive flap (φs) 0.6664 * 0.6399 Blade progressive flap (φc) 0.6675 * 0.6542 Blade
collective flap (φ0) 0.6993 * 0.6879 Blade regressive torsion (εs) - 5.5705 Blade progressive
torsion (εc) - 5.5738 Blade collective torsion (ε0) 5.6 ~ 5.5708 * From Ref.[8], ~ From Ref.[5]
Figure 4. Natural frequencies variation considering rotational speed changes with wind
velocities using 13 DOF model increase in wind velocity. Aeroelastic flap modes will be
accurately predicted by flexible beam models which captures aerodynamic coupling
accurately. Differences in the aeroelastic natural frequencies reported in Fig. 5 and Fig. 6
arise due to simplification of flexible structure with rigid beam model and can be attributed
to two reasons: centrifugal stiffening and aerodynamic loads. The rigid beam model doesn’t
account for centrifugal stiffening of the rotating beam in bending. In the current model, wind
turbine blades are assumed as rigid beams with torsional spring attached at the blade root. If
the blade is modeled as flexible, lower part of the blade vibrates with high amplitudes as it
experiences higher loads and also it is less stiff when compared to the blade geometry near
the root. So, rigid beam approximation for the blade doesn’t take this into account.