Next Article in Journal
Evaluation of TMD Performance in Footbridges Using Human Walking Probabilistic Models
Previous Article in Journal
Exploring the Relationship between Preprocessing and Hyperparameter Tuning for Vibration-Based Machine Fault Diagnosis Using CNNs
Previous Article in Special Issue
Combining Computational Fluid Dynamics and Gradient Boosting Regressor for Predicting Force Distribution on Horizontal Axis Wind Turbine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Finite Element Analysis of Wind Turbine Blade Vibrations

1
ACEX Building, Department of Aeronautics, South Kensington Campus, Imperial College London, London SW7 2AZ, UK
2
Department of Theoretical and Applied Mechanics, Taras Shevchenko National University of Kyiv, 01601 Kyiv, Ukraine
3
Department of Machinery, Electronics, and Chemistry, Oslo Metropolitan University, Pilestredet 35, 0166 Oslo, Norway
4
ICP Institute of Computational Physics, ZHAW Zürich University of Applied Sciences, Wildbachstrasse 21, CH-8401 Winterthur, Switzerland
*
Author to whom correspondence should be addressed.
Submission received: 31 January 2021 / Revised: 23 March 2021 / Accepted: 30 March 2021 / Published: 5 April 2021
(This article belongs to the Special Issue Dynamics of Composite Wind Turbine Rotor Blades)

Abstract

:
The article is devoted to the practical problem of computer simulation of the dynamic behaviour of horizontal axis wind turbine composite rotor blades. This type of wind turbine is the dominant design in modern wind farms, and as such its dynamics and strength characteristics should be carefully studied. For this purpose, in this paper the mechanical model of a rotor blade with a composite skin possessing a stiffener was developed and implemented as a finite element model in ABAQUS. On the basis of this computer model, modal analysis of turbine blade vibrations was performed and benchmark cases for the dynamic response were investigated. The response of the system subjected to a uniform underneath pressure was studied, and the root reaction force and blade tip displacement time histories were obtained from the numerical calculations conducted.

1. Introduction

Today, wind energy along with its green counterpart, solar energy, is considered as a predominant source and a pillar of ‘green’ electricity. Currently, the wind power sector is the fastest growing renewable energy sector and contributes significantly to reducing greenhouse gases and to achieving the emission reduction targets set by the Kyoto Protocol and Paris Agreement. Furthermore, the rapidly developing improvements in wind turbines’ design and efficiency increase their economic competitiveness compared to conventional power generating technologies [1].
Over the past decades of the evolution of wind energy technology, the main tendency has been increasing size of wind turbine rotors. This is driven by both economic and technological demands, as larger turbines have enhanced performance and increased production capacity, giving rise to new high-capacity wind power plants [2,3]. In turn, scaling up wind turbines leads to the emergence of new complex problems in their design and operation, for single turbines as well as in wind farm spacing. The solution of these problems requires, along with immoderate experimental investigations, intensive computer simulation of dynamic processes in all their parts, especially in the blades, which in the case of large turbines are now made of composite materials and have a rather complex structural design. Modern turbines have significantly longer blades than their predecessors (tens of meters) and are much more flexible. These blades are designed as long slender composite shell structures with stiffeners, significantly inhomogeneous and becoming more slender towards the free tip. This necessitates a wide range of studies in blade structural dynamics of the performance of such blades to evaluate deformations and stresses in their structure. The possible intensive vibrations in service conditions may not only cause fatigue damage or affect structural strength and stability, but also have significant influence on power production [4]. Another major and important area where simulation of dynamic processes in wind turbines is needed is the vibration control of wind turbines during operation and at extreme conditions [1,5].
For the computer simulation of blades’ dynamic behavior, a number of modern FEM packages like ANSYS or ABAQUS are widely used. They modelling of blades having high accuracy and simulate their dynamic behaviour, taking into account nonlinear geometric and material effects when necessary [6,7].
Simpler beam-type models can be obtained directly using different methods of the blade behaviour approximation. In the past, models based on advanced variants of theory for thin-walled beams have widely used [8,9,10], or an identification approach has been used based on the data obtained from physical or numerical experiments (the latter was considered in [11]).
The aim of the present study is to investigate dynamic characteristics of a horizontal axis wind turbine (HAWT) rotor blade with a stiffener. For this, a corresponding mechanical model was developed and implemented by means of a finite element package ABAQUS [12]. Then a modal analysis of turbine blade vibrations was carried out in order to obtain an eigenfrequency spectrum and shapes of the turbine blade vibration modes, which can be used in transient modal based dynamic analysis and, in future studies, in building different simplified beam-type by identification methods as in [11]. Besides, benchmark cases for the dynamic response of the blade subjected to uniform underneath pressure pulses of different magnitudes were studied in detail to find the loading value that provides a noticeable effect of geometric nonlinearity. The evaluated root reaction forces allow estimation of strength characteristics of the blade in the root section.

2. Methodology

In order to investigate the response of real rotor blades, an ABAQUS shell-type model of the blade was developed to be consistent with industrial demands. Inertial effects due to the rotation were not considered in the present paper but can be taken into account in further investigations. The blade consisted of the skin, a stiffener and a specific attachment zone at the root of the blade. The model incorporated all necessary FEM modeling features for effective simulation of the transient blade response to different types of loading. Blade geometry was a modified version [13] of a 3 MW turbine proposed by [14] and used a scale-up version of the 20.15 m long blade in [15]. The geometry of the model was based on that developed in [11] and is shown in Figure 1.
The rotor blade, which was 44.175 m long, had a circular cross-section at its root attaching the rotor hub to the 1.753-m span location. Then the circular blade cross-sections was gradually transformed to airfoil sections that had NREL airfoil types S818 at the 9.648-m span location, S825 at 32.667-m, and S826 at 41.873-m and 44.175-m (at the blade tip). The shapes of the NREL airfoils used are shown in Figure 1 in [13], and in [16], where data files presenting details of geometry of these airfoils are also given. Table 1 presents the details of the cross-sections framing the blade, and the constitutive parameters of the blade’s skin material which was considered to be a homogeneous unidirectional composite with fibers oriented along the blade as shown in Table 2 [11].
For representation of the considered shell structure, ABAQUS S4R 4-node shell elements with reduced integration and hourglass control were utilized. This element type belongs to the class of general-purpose elements suitable for both thin and thick shells and allowed us to take into account finite membrane strains and arbitrarily large rotations [12]. A brief review of different modern approaches to developing shell elements and comparatively evaluating some of them, including S4 and S4R, can be found in [17].
The model was discretized into finite elements by mapped meshing with element size in z direction about 20 cm, and with 16 subdivisions in transverse direction at the top and bottom sides of the blade (overall more than 7000 elements were used).
At the root of the blade the encastre boundary condition was applied which constrains all active structural degrees of freedom (i.e., displacements and rotations): U 1 = U 2 = U 3 = 0 , U R 1 = U R 2 = U R 3 = 0 .
Some computations at early stage of investigation were carried out using S4 shell elements with full integration which led to more expensive models. However, comparison of results with the ones obtained with S4R elements showed small difference between eigenfrequencies (fractions of a percent), so the S4R element was chosen.
The cross-section of the blade was strengthened by a wooden stiffener (made of longleaf pine [18]) with the values of density, stiffness moduli and Poisson ratios listed in Table 3. The stiffener was chosen to be running along the whole blade as shown in Figure 1 by the red line. The rear view of the blade root and the stiffener is depicted in Figure 2.

3. Results

3.1. Modal Analysis of the Wind Turbine Blade Vibrations

In this subsection modal analysis of a horizontal axis wind turbine (HAWT) composite blade is conducted using the finite element technique. Calculations were performed for three values of the wooden (longleaf pine) stiffener thickness: 100 mm, 150 mm, and 200 mm. Correspondent values of the eigenfrequencies for the eight lowest vibration modes are listed in Table 4, and their shapes are shown in Figure 3.
Analyses of Table 4 and Figure 3 show that, as can be seen from the vibration mode shapes represented in Figure 3, the modes 1, 2, 4 had a flap character, the 3rd mode was mainly lead-lag, and in the modes 5, 6, 8 flap and lead-lag movements were combined. The vibration mode 7 was a twisting mode.
A blade with a wooden stiffener of 150 mm thickness was chosen for the following simulations and analyses.

3.2. Simulation of Blade Dynamics under Step-Wise Pulse Loading

As mentioned above, a blade consisting of a composite skin along with a wooden stiffener of 150 mm thickness was subjected to step-wise uniform underneath pressure of variable magnitudes. The aim of the simulations performed was to find the loading value that provided a noticeable effect of geometric nonlinearity. Two simulations were run for each particular case of the loading amplitude. The first one was performed using the ABAQUS model, taking into account geometric nonlinearity, while the linear elastic ABAQUS model was employed for the second simulation. The computations were done using an implicit direct-integration dynamic procedure provided in ABAQUS/Standard [12] with constant time steps in linear cases, and with variable steps in nonlinear ones. The different time steps were used for calculations. Stable results were obtained with application of the practical convergence technique for time step choice.
Special attention was paid to the investigation of time histories of the blade tip displacement in the y (flap) direction and the z (longitudinal) component of the root reaction force which were obtained from the results of calculations.
Typical results for the blade tip point deflection time history are shown in Figure 4. In the figures in the sequel, deflection levels are in meters while time is in seconds. The results of calculations showed that for up to a magnitude of pressure (in the case considered here 7 × 102 Pa) the influence of nonlinearity was small and may be neglected. However, for larger values of pressure, nonlinear effects such as increasing of amplitude and time lag of the vibrations were observed.
A special procedure was developed to find the corrected distribution of the reaction force results over the root’s circular cross-section because the nodal points of the mesh were not evenly distributed (see red dots in Figure 2). This caused the distribution presented in Figure 5 where the nodal values of reaction forces R z are in Newtons with respect to azimuth angle θ as shown (calculations were performed for pressure loading magnitude 102 Pa, and this graph corresponds to the time t = 0.9 s when the maximum values of the reaction forces were reached).
For this aim, a standard Fourier series approximation was used:
f ( θ ) a 0 2 + n = 1 N a n cos n θ + b n sin n θ , a 0 = 1 π π π f ( θ ) d θ , a n = 1 π π π f ( θ ) cos n θ d θ , b n = 1 π π π f ( θ ) sin n θ d θ .
The typical corrected results for the density of membrane reaction force R ¯ z (N/m) are shown in Figure 6 where the effect of the stiffener is clearly visible as the noticeable rise of the membrane reaction occurred at places of stiffener attachment (90° and −90°). This effect should be taken into account when the strength, fatigue resistance and life prediction of the blade are evaluated.
Graphic representation of the time histories of the reaction force R ¯ z density at 90° (on y-axis, attachment point of the stiffener and skin) obtained by the Fourier approximation of the 9th order (solid line corresponds to the geometrically linear model while dashed line stands for the nonlinear one) for different pressure loading magnitudes is presented in Figure 7.
It can be seen that both the linear and nonlinear models predicted similar responses for the loads which did not exceed a value of 7 × 102 Pa, which implies insignificant geometric nonlinear effects up to this point in the loading space.

3.3. Simulation of Blade Dynamics under Half-Sine Pulse Loading

Here, the results of response simulation for the wind turbine rotor blade subjected to the uniform underneath pressure pulse of the half-sine shape with a duration of 2 s (see Figure 8) are presented. This pulse imitates the load at a wind gust with a single peak. The investigated blade’s composition was a UD composite skin along with a wooden stiffener 150 mm thick.
Displacements of the blade tip are shown in Figure 9 for different load levels of the half-sine pulse pressure load. In this figure, the solid lines correspond to linear simulation while dashed lines are used to depict the case when nonlinear geometrical effects are accounted for.
Nonlinearity effects started to become noticeable from the magnitude of pressure p around 7 × 102 Pa upwards. For the loading magnitude of 103 Pa, the discrepancy between linear and nonlinear cases were more pronounced. It should be mentioned that the nonlinearity level could be estimated by the displacement of nodal points time shift (as well as the displacement extremes). A magnitude of 103 Pa led to the essential movement of the extreme points of graph along time axis to the right. Thus, both displacement extremes and the nodal displacement shift in time if nonlinearity is taken into account. This is due to the dependence of frequency on amplitude and on considering the deformed configuration derived from the equations of motion.
The membrane reaction force densities at the blade root are depicted in Figure 10a–e, which correspond to load magnitudes of 102, 5 × 102, 7 × 102, 7.5 × 102 and 103 Pa, respectively. These results are similar to ones obtained in the previous section.
The results were qualitatively similar to the ones shown in Figure 9 for displacements up to those corresponding to the magnitude of pressure p close to the value 7 × 102 Pa. For larger loading amplitudes, fast oscillations due to nonlinear effects occurred (see graphs in Figure 10d,e. The magnitude of the oscillations increased rapidly with the load. Presence of the effect can be explained by the longitudinal vibration of the blade which was taken into account for the nonlinear blade. Thus this additional “shivering effect” can be attributed to geometrical nonlinearity.
In Figure 11 and Figure 12 the results for simulations up to 10 s for the half-sine pulse with a duration of 2 s with pressure magnitude 103 Pa are presented. The comparison of Figure 9 and Figure 11 shows that the general tendencies were kept: displacement nodal points shifted to the right (as well as the displacement extremes), i.e., the blade demonstrated nonlinear responses of the soft type. The same was observed for the membrane reaction force graphs in Figure 12 (compare with Figure 10e).
To verify the computation procedure, calculations with different time steps were performed. The results agreed with high accuracy. Therefore, the “shivering effect” was physical and is not introduced by computational instability.
It is worth mentioning here that accounting for the geometrical nonlinearity manifests itself in three effects: (i) increasing of the flap deflection up to about 15 per cent at a load magnitude of 103 Pa; (ii) gradually increasing time shift in the nodal displacements and their extremes and (iii) “shivering” in the root membrane forces due to excitation of longitudinal vibrations up to about 20 per cent at the load magnitude 103 Pa.

4. Discussion and Conclusions

In the present study, a shell type mechanical model of a horizontal axis wind turbine (HAWT) rotor blade with a stiffener was developed and implemented in the finite element package ABAQUS. This model allows for detailed investigation of the dynamic response and processes of deformation and evolution of stresses in the rotor blade considered.
Firstly, an eigenvalue modal analysis (linear perturbation analysis) of turbine blade vibrations was carried out. It demonstrated the wide variety of the forms of vibrational motion for the blade, from bending modes in flap and lead-lag directions and twisting modes in their combinations.
Benchmark cases for the dynamic response were investigated. Two types of pulse loading were considered: step-wise, replicating a pulse load as impact or blast, and half-sine shape, which imitates the load at a wind gust. In both cases a uniform pressure pulse was applied to the underneath part of the blade. The computations were performed both in the linear and geometrically nonlinear formulations. Blade tip displacement and root reaction force time histories were obtained from the results of finite element simulation conducted using ABAQUS/Standard. The results of calculations showed that for up to some magnitude of pressure (in the case considered here 7 × 102 Pa) the influence of nonlinearity was small and may be neglected. However, with larger values of pressure, nonlinear effects were observed in the behaviour of reaction forces at the blade’s root, such as gradually increasing amplitude and time lag of the vibrations, and the emergence of high frequency longitudinal vibrations (similar to appearance of overtones in musical instruments).
In general, it was found the geometrical nonlinearity had three main manifestations under the conditions considered: (i) increasing of the flap deflection up to about 15 per cent at the load magnitude 103 Pa; (ii) gradually increasing time shift in the nodal displacements and their extremes and (iii) “shivering” in the root membrane forces due to excitation of longitudinal vibrations up to about 20 per cent at a load magnitude 103 Pa.
The developed finite element model of the blade can be modified to investigate behaviour of the blades made of layered composites including different approaches to damage consideration. Using the results of simulation, both modal analysis and transient dynamics and different types of approximated beam-type models of wind turbine blades on the base of identification techniques can be built [12].

Author Contributions

Conceptualization, N.N., Y.Z. and A.S.F.; methodology, I.G., F.E.M. and Y.Z.; software simulations, N.N., I.G. and F.E.M.; formal analysis of results, N.N., Y.Z., A.S.F.; writing—original draft preparation, N.N., I.G. and F.E.M.; writing—review and editing, A.S.F. and Y.Z.; supervision, A.S.F. and Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to being part of an ongoing study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Staino, A.; Basu, B. Dynamics and control of vibrations in wind turbines with variable rotor speed. Eng. Struct. 2013, 56, 58–67. [Google Scholar] [CrossRef]
  2. Serrano-González, J.; Lacal-Arántegui, R. Technological Evolution of Onshore Wind Turbines—A market-Based Analysis. Wind Energy 2016, 19, 2171–2187. [Google Scholar] [CrossRef] [Green Version]
  3. Wiser, R.H.; Bolinger, M. 2015 Wind Technologies Market Report. U.S. Department of Energy. 2016. Available online: https://emp.lbl.gov/publications/2015-wind-technologies-market-report (accessed on 21 January 2021).
  4. Brondsted, P.; Nijssen, R.P.L. (Eds.) Advances in Wind Turbine Blade Design and Materials; Woodhead Publishing Limited: Cambridge, UK, 2013; p. 477. [Google Scholar]
  5. Staino, A.; Basu, B. Emerging Trends in Vibration Control of Wind Turbines: A Focus on a Dual Control Strategy. Phil. Trans. R. Soc. A 2015, 373, 20140069. [Google Scholar] [CrossRef] [PubMed]
  6. Laird, D.; Montoya, F.; Malcolm, D. Finite element modeling of wind turbine blades. ASME/AIAA Wind Energy Symp. Reno 2005, 9–17. [Google Scholar] [CrossRef]
  7. Tarfaoui, M.; Nachtane, M.; Boudounit, H. Finite element analysis of composite offshore wind turbine blades under operating conditions. J. Therm. Sci. Eng. Appl. 2020, 12, 011001. [Google Scholar] [CrossRef]
  8. Qin, Z.; Librescu, L. Dynamic aeroelastic response of aircraft wings modeled as anisotropic thin-walled beams. J. Aircraft. 2003, 40, 532–543. [Google Scholar] [CrossRef]
  9. Cárdenas, D.; Escárpita, A.A.; Elizalde, H.; Aguirre, J.J.; Ahuett, H.; Marzocca, P.; Probst, O. Numerical validation of a finite element thin-walled beam model of a composite wind turbine blade. Wind Energy 2012, 15, 203–223. [Google Scholar] [CrossRef]
  10. Le, T.B.; Christenson, A.; Calderer, T.; Stolarski, H.; Sotiropoulos, F. A thin-walled composite beam model for light-weighted structures interacting with fluids. J. Fluids Struct. 2020, 95, 102968. [Google Scholar] [CrossRef]
  11. Navadeh, N.; Goroshko, I.O.; Zhuk, Y.A.; Fallah, A.S. An FEM-based AI approach to model parameter identification for low vibration modes of wind turbine composite rotor blades. Europ. J. Comp. Mech. 2017, 26, 541–556. [Google Scholar] [CrossRef]
  12. ABAQUS 6.14 Documentation. 2014. Available online: http://130.149.89.49:2080/v6.14/ (accessed on 21 January 2021).
  13. Chen, K.-N.; Chen, P.-Y. Structural optimization of 3 MW wind turbine blades using a two-step procedure. Int. J. Simul. Multidisci. Des. Optim. 2010, 4, 159–165. [Google Scholar] [CrossRef] [Green Version]
  14. Malcolm, D.J.; Hansen, A.C. WindPACT Turbine Rotor Design Study; Subcontract Report NREL/SR-500-32495; National Renewable Energy Laboratory: Golden, CO, USA, 2006. Available online: https://www.nrel.gov/docs/fy06osti/32495.pdf (accessed on 21 January 2021).
  15. Bir, G.; Migliore, P. Preliminary Structural Design of Composite Blades for Two-and Three-Blade Rotors; Technical Report NREL/TP-500-31486; National Renewable Energy Laboratory: Golden, CO, USA, 2004.
  16. Airfoil Database List (S). Available online: http://airfoiltools.com/search/list?page=s&no=4 (accessed on 21 January 2021).
  17. Laulusa, A.; Bauchau, O.A.; Choi, J.-Y.; Tan, V.B.C.; Li, L. Evaluation of some shear deformable shell elements. Int. J. Solids Struct. 2006, 43, 5033–5054. [Google Scholar] [CrossRef] [Green Version]
  18. Wood Handbook—Wood as an Engineering Material; General Technical Report FPL-GTR-190; U.S. Department of Agriculture, Forest Service, Forest Products Laboratory: Madison, WI, USA, 2010. Available online: https://www.fpl.fs.fed.us/documnts/fplgtr/fpl_gtr190.pdf (accessed on 21 January 2021).
Figure 1. The rotor blade model.
Figure 1. The rotor blade model.
Vibration 04 00020 g001
Figure 2. Root part of the blade.
Figure 2. Root part of the blade.
Vibration 04 00020 g002
Figure 3. Shapes of the turbine blade vibration modes.
Figure 3. Shapes of the turbine blade vibration modes.
Vibration 04 00020 g003aVibration 04 00020 g003b
Figure 4. Displacements of the blade tip point for different load levels (shown in Pa at the lines) for the step-wise pressure loading: solid lines—linear simulation, dashed lines—nonlinear geometrical effects included.
Figure 4. Displacements of the blade tip point for different load levels (shown in Pa at the lines) for the step-wise pressure loading: solid lines—linear simulation, dashed lines—nonlinear geometrical effects included.
Vibration 04 00020 g004
Figure 5. Nodal values of reaction forces (N).
Figure 5. Nodal values of reaction forces (N).
Vibration 04 00020 g005
Figure 6. Density of membrane reaction force in z-direction obtained by the Fourier approximation of the 9th order.
Figure 6. Density of membrane reaction force in z-direction obtained by the Fourier approximation of the 9th order.
Vibration 04 00020 g006
Figure 7. Membrane reaction force density time histories for step-wise pulse loading.
Figure 7. Membrane reaction force density time histories for step-wise pulse loading.
Vibration 04 00020 g007
Figure 8. Half-sine pulse shape.
Figure 8. Half-sine pulse shape.
Vibration 04 00020 g008
Figure 9. Displacement time histories of the blade tip for different levels of the half-sine pulse pressure load (indices near graphs): solid lines—linear simulations, dashed lines—simulation with nonlinear geometrical effects included.
Figure 9. Displacement time histories of the blade tip for different levels of the half-sine pulse pressure load (indices near graphs): solid lines—linear simulations, dashed lines—simulation with nonlinear geometrical effects included.
Vibration 04 00020 g009
Figure 10. Membrane reaction force density time histories for half-sine pulse loading: (a)—102, (b)—5 × 102, (c)—7 × 102, (d)—7.5 × 102 and (e)—103 Pa.
Figure 10. Membrane reaction force density time histories for half-sine pulse loading: (a)—102, (b)—5 × 102, (c)—7 × 102, (d)—7.5 × 102 and (e)—103 Pa.
Vibration 04 00020 g010aVibration 04 00020 g010b
Figure 11. Displacement at the tip of the blade for the 103 Pa magnitude of half-sine pulse pressure load: solid lines—linear simulation, dashed lines—nonlinear geometrical effects included.
Figure 11. Displacement at the tip of the blade for the 103 Pa magnitude of half-sine pulse pressure load: solid lines—linear simulation, dashed lines—nonlinear geometrical effects included.
Vibration 04 00020 g011
Figure 12. Membrane reaction force density time histories for 103 Pa magnitude of half-sine pulse pressure load: solid lines—linear simulation, dashed lines—nonlinear geometrical effects included.
Figure 12. Membrane reaction force density time histories for 103 Pa magnitude of half-sine pulse pressure load: solid lines—linear simulation, dashed lines—nonlinear geometrical effects included.
Vibration 04 00020 g012
Table 1. Section thicknesses along the beam axis (measured from the circular root).
Table 1. Section thicknesses along the beam axis (measured from the circular root).
Number of SectionDistance from the Root in z-Direction (mm)Thickness t (mm)
1175496.02
2964879.04
314,25153.21
429,21213.44
544,17510.34
Table 2. Stiffness moduli and Poisson ratios for the UD composite (1—fibre direction, 2—transverse direction, 3—through-the-thickness direction).
Table 2. Stiffness moduli and Poisson ratios for the UD composite (1—fibre direction, 2—transverse direction, 3—through-the-thickness direction).
ρ (kg/m3)E1E2E3 ν 12 ν 13 ν 23 G12G13G23
(GPa)(GPa)(GPa)(GPa)(GPa)(GPa)
186037.018.018.00.250.250.38.08.06.9
Table 3. Stiffness moduli and Poisson ratios for the longleaf pine.
Table 3. Stiffness moduli and Poisson ratios for the longleaf pine.
ρ (kg/m3)E1E2E3 ν 12 ν 13 ν 23 G12G13G23
(GPa)(GPa)(GPa)(GPa)(GPa)(GPa)
59015.071.53710.82890.3320.3650.3841.070.90420.1808
Table 4. Eigenfrequencies of the rotor blade (UD composite skin and wooden stiffener) vibration modes for different stiffener thicknesses.
Table 4. Eigenfrequencies of the rotor blade (UD composite skin and wooden stiffener) vibration modes for different stiffener thicknesses.
Mode Number NEigenfrequency, Hz
Stiffener Thickness
100 mm150 mm200 mm
10.510220.505290.50076
21.730441.707931.68643
32.169162.102152.04292
44.016883.976063.93519
56.523196.383456.2385
67.090876.986046.90093
77.887667.977028.07623
89.561389.568429.5557
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Navadeh, N.; Goroshko, I.; Zhuk, Y.; Etminan Moghadam, F.; Soleiman Fallah, A. Finite Element Analysis of Wind Turbine Blade Vibrations. Vibration 2021, 4, 310-322. https://0-doi-org.brum.beds.ac.uk/10.3390/vibration4020020

AMA Style

Navadeh N, Goroshko I, Zhuk Y, Etminan Moghadam F, Soleiman Fallah A. Finite Element Analysis of Wind Turbine Blade Vibrations. Vibration. 2021; 4(2):310-322. https://0-doi-org.brum.beds.ac.uk/10.3390/vibration4020020

Chicago/Turabian Style

Navadeh, Navid, Ivan Goroshko, Yaroslav Zhuk, Farnoosh Etminan Moghadam, and Arash Soleiman Fallah. 2021. "Finite Element Analysis of Wind Turbine Blade Vibrations" Vibration 4, no. 2: 310-322. https://0-doi-org.brum.beds.ac.uk/10.3390/vibration4020020

Article Metrics

Back to TopTop