0
Research Papers

A Continuum Based Three-Dimensional Modeling of Wind Turbine Blades

[+] Author and Article Information
Ahmed H. Bayoumy

Graduate Student
Mechanical Design & Production,
Engineering Department,
Faculty of Engineering,
Cairo University,
Giza 12613, Egypt
e-mail: a.hamdy85@ymail.com

Ayman A. Nada

Assistant Professor
Mechanical Engineering Department,
Benha Institute of Technology,
Benha University,
Benha 13512, Egypt
e-mail: arobust@tedata.net.eg

Said M. Megahed

Professor
Mechanical Design & Production,
Engineering Department,
Faculty of Engineering,
Cairo University,
Giza 12613, Egypt
e-mail: smegahed@cu.edu.eg

1Corresponding author.

2Present address: College of Engineering, Jazan University, Jazan P.O.706, KSA.

Contributed by the Design Engineering Division of ASME for publication in the Journal of Computational and Nonlinear Dynamics. Manuscript received December 5, 2011; final manuscript received August 18, 2012; published online October 30, 2012. Assoc. Editor: Khaled E. Zaazaa.

J. Comput. Nonlinear Dynam 8(3), 031004 (Oct 30, 2012) (14 pages) Paper No: CND-11-1237; doi: 10.1115/1.4007798 History: Received December 05, 2011; Revised August 18, 2012

Accurate modeling of large wind turbine blades is an extremely challenging problem. This is due to their tremendous geometric complexity and the turbulent and unpredictable conditions in which they operate. In this paper, a continuum based three dimensional finite element model of an elastic wind turbine blade is derived using the absolute nodal coordinates formulation (ANCF). This formulation is very suitable for modeling of large-deformation, large-rotation structures like wind turbine blades. An efficient model of six thin plate elements is proposed for such blades with non-uniform, and twisted nature. Furthermore, a mapping procedure to construct the ANCF model of NACA (National Advisory Committee for Aeronautics) wind turbine blades airfoils is established to mesh the geometry of a real turbine blade. The complex shape of such blades is approximated using an absolute nodal coordinate thin plate element, to take the blades tapering and twist into account. Three numerical examples are presented to show the transient response of the wind turbine blades due to gravitational/aerodynamics forces. The simulation results are compared with those obtained using ANSYS code with a good agreement.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.

References

Kalleϕe, B. S., 2007, “Aeroservoelasticity of Wind Turbines,” Ph.D. thesis, Technical University of Denmark.
Holm-Jϕrgensen, K., and Nielsen, S., 2009, “System Reduction in Multibody Dynamics of Wind Turbines,” Multibody Syst. Dyn., 21(2), pp.147–165. [CrossRef]
ShabanaA. A., 2005, Dynamics of Multibody Systems, 3rd ed., Cambridge University, Cambridge, England.
Nada, A. A., Hussein, B. A., Megahed, S. M., and Shabana, A. A., 2010, “Use of the Floating Frame of Reference Formulation in Large Deformation Analysis: Experimental and Numerical Validation,” Proc. Inst. Mech. Eng., Part K: J. Multibody Dyn., 224, pp. 45–58. [CrossRef]
Dufva, K., and Shabana, A. A., 2005, “Analysis of Thin Plate Structure Using the Absolute Nodal Coordinate Formulation,” Proc. Inst. Mech. Eng., Part K: J. Multibody Dyn., 219, pp. 345–355. [CrossRef]
Dufva, K. E., Kerkkänen, K. S., Maqueda, L., and Shabana, A. A., 2007, “Nonlinear Dynamics of Three-dimensional Belt Drives Using the Finite-Element Method,” Nonlinear Dyn., 48(4), pp. 449–466. [CrossRef]
Yoo, W. S., Lee, J. H., Park, S. J., Sohn, J. H., Pogorelov, D., and Dmitrochenko, O., 2004, “Large Deflection Analysis of a Thin Plate: Computer Simulations and Experiments,” Multibody Syst. Dyn., 11, pp. 185–208. [CrossRef]
Vallejo, D. G., Sugiyama, H., and Shabana, A. A., 2005, “Finite Element Analysis of the Geometric Stiffening Effect Using the Absolute Nodal Coordinate Formulation,” Proceedings of IDETC/CIE, ASME 2005 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference, Sept. 24–28, Long Beach, California.
Nada, A. A., and El-Assal, A. M., 2012, “Absolute Nodal Coordinate Formulation of Large-Deformation Piezoelectric Laminated Plates,” Nonlinear Dyn., 67(4), pp. 2441–2454. [CrossRef]
Gerstmayr, J., and Shabana, A. A., 2006, “Analysis of Thin Beams and Cables Using the Absolute Nodal Coordinate Formulation,” Nonlinear Dyn., 45, pp. 109–130. [CrossRef]
Shabana, A. A., 2010, Computational Continuum Mechanics, Cambridge University, Cambridge, England.
Yakoub, Y. R., and Shabana, A. A., 2001, “Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Implementation and Applications,” ASME J. Mech. Des., 123, pp. 614–621. [CrossRef]
Shabana, A. A., and Yakoub, R. Y., 2001, “Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Theory,” ASME J. Mech. Des., 123(4), pp. 606–613. [CrossRef]
Vallejo, D. G., Mayo, J., Escalona, J. L., and Dominguez, J., 2004, “Efficient Evaluation of the Elastic Forces and the Jacobian in the Absolute Nodal Coordinate Formulation,” Nonlinear Dyn., 35(4), pp. 313–329. [CrossRef]
Hussein, B. A., Sugiyama, H., and Shabana, A. A., 2007, “Coupled Deformation Modes in the Finite Element Absolute Nodal Coordinate Formulation,” 12th IFToMM World Congress, Besançon, France, June 18–21.
Dmitrochenko, O. N., Hussein, B. A., and Shabana, A. A., 2009, “Coupled Deformation Modes in the Large Deformation Finite Element Analysis: Generalization,” J. Comput. Nonlinear Dyn., 4(2), p. 021002. [CrossRef]
Hussein, B. A., Negrut, D., and Shabana, A. A., 2008, “Implicit and Explicit Integration in the Solution of the Absolute Nodal Coordinate Differential/Algebraic Equations,” Nonlinear Dyn., 54(4), pp. 283–296. [CrossRef]
Hussein, B. A., and Shabana, A. A., 2011, “Sparse Matrix Implicit Numerical Integration of the Stiff Differential/Algebraic Equations: Implementation,” Nonlinear Dyn., 65, pp. 369–382. [CrossRef]
Sugiyama, H., Escalona, J., and Shabana, A. A., 2003, “Formulation of Three-Dimensional Joint Constraints Using the Absolute Nodal Coordinates,” Nonlinear Dyn., 33, pp. 167–195. [CrossRef]
Vallejo, D. G., Valverde, J., and Dominguez, J., 2005, “An Internal Damping Model for the Absolute Nodal Coordinate Formulation,” Nonlinear Dyn., 42, pp. 347–369. [CrossRef]
Abdel-Nasser, M., and Shabana, A., 2011, “A Nonlinear Visco-elastic Constitutive Model for Large Rotation Finite Element Formulations,” Multibody Syst. Dyn., 26(1), pp. 57–79. [CrossRef]
Dreese, J., 2000, Aero Basics and Design FOIL User Guide, Capitola, California.
Carmichael, R. L., 2001, “Algorithm for Calculating Coordinates of Cambered NACA Airfoils at Specified Chord Locations,” AIAA, pp. 5235-1–5235-8.
Farin, G., 2001, Curves and Surfaces for CAGD A Practical Guide, 5th ed., Morgan Kaufmann Publishers, San Francisco, CA.
El Chazly, N. M., 1993, “Static and Dynamic Analysis of Wind Turbine Blades Using the Finite Element Method,” Renewable Energy, 3(6–7), pp. 705–724. [CrossRef]
Bathe, K. J., 1996, Finite Element Procedures, Prentice-Hall, Engelwood Cliffs, NJ.
Sereshk, M. V., and Salimi, M., 2011, “Comparison of Finite Element Method Based on Nodal Displacement and Absolute Nodal Coordinate Formulation (ANCF) in Thin Shell Analysis,” Int. J. Numer. Methods Biomed. Eng., 27(8), pp. 1185–1198. [CrossRef]
Sugiyama, H., Koyama, H., and Yamashita, H., 2010, “Gradient Deficient Curved Beam Element Using the Absolute Nodal Coordinate Formulation,” J. Comput. Nonlinear Dyn., 5, p. 021001. [CrossRef]
Jureczko, M., Pawlak, M., and Mezyk, A., 2005, “Optimisation of Wind Turbine Blades,” J. Mater. Process. Technol., 167(2–3), pp. 463–471. [CrossRef]
Houghton, E. L., and Carpenter, P. W., 2003, Aerodynamics for Engineering Students, 5th ed., Butterworth-Heinemann, London.
Shabana, A. A., 2010, Computational Dynamics, 3rd ed., Wiley, New York.
ANSYS User's Manual, 2001, 12th ed., SAS IP, Canonsburg.
Cheung,Y. K., Zhang,Y. X., and Wanji, C., 2000, “The Application of a Refined Non-conforming Quadrilateral Plate Bending Element in Thin Plate Vibration and Stability Analysis,” Finite Elem. Anal. Design, 34(2), pp. 175–191. [CrossRef]
Hua, X., Rui, G. U., Jin, J., Liu,Y., Yi, M. A., Cong, Q., and Zheng, Y., 2010, “Numerical Simulation and Aerodynamic Performance Comparison Between Seagull Aerofoil and NACA 4412 Aerofoil Under Low-Reynolds,” Adv. Nat. Sci., 3(2), pp. 244–250.
Rajakumar, S., and Ravindran, D., 2010, “Computational Fluid Dynamics Of Wind Turbine Blade at Various Angles of Attack and Low Reynolds Number,” Int. J. Eng. Sci. Technol., 2(11), pp. 6474–6484.

Figures

Grahic Jump Location
Fig. 1

Airfoil shape parameters

Grahic Jump Location
Fig. 3

Four nodes plate element

Grahic Jump Location
Fig. 5

Constructing airfoil using ANCF (2 elements)

Grahic Jump Location
Fig. 6

NACA and ANCF airfoils (2 elements)

Grahic Jump Location
Fig. 7

Constructing airfoil using ANCF (6 elements)

Grahic Jump Location
Fig. 8

ANCF models of NACA airfoils, red lines represent NACA code profile

Grahic Jump Location
Fig. 9

Nonuniform wind turbine blade

Grahic Jump Location
Fig. 10

WTB NACA 4412, α = 5 (deg)

Grahic Jump Location
Fig. 11

WTB NACA 8612, α = 8 (deg)

Grahic Jump Location
Fig. 12

Twisted wind turbine blade

Grahic Jump Location
Fig. 13

WTB: twisted NACA 4412, β = 30 (deg)

Grahic Jump Location
Fig. 14

WTB: front view: twisted NACA 4412, β = 30 (deg)

Grahic Jump Location
Fig. 15

WTB: twisted NACA 8612, β = 30 (deg)

Grahic Jump Location
Fig. 16

WTB: front view: twisted NACA 8612, β = 30 (deg)

Grahic Jump Location
Fig. 17

Aerodynamic forces

Grahic Jump Location
Fig. 18

FEM blade model with ANSYS (nontapered blade)

Grahic Jump Location
Fig. 19

Transverse deflection of tip edge point

Grahic Jump Location
Fig. 20

FEM blade model with ANSYS (tapered blade)

Grahic Jump Location
Fig. 21

Transverse deflection of tip edge point

Grahic Jump Location
Fig. 22

ANSYS meshing domains of blade-wall and fluid

Grahic Jump Location
Fig. 23

Transient response of aerodynamic forces with ANSYS

Grahic Jump Location
Fig. 24

Comparison of aerodynamics responses

Grahic Jump Location
Fig. 25

WTB model with multisections along the span length

Grahic Jump Location
Fig. 26

Comparison of different models along the span length

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In