Wind power industry is expected to be one of the fastest growing renewable energy sources inthe world. The growth specially focuses on growing industries and markets, because ofeconomical condition for wind power development besides political decisions.According to growth of wind turbine industries, wind turbine blades are growing fast in both sizeand number. The problem that now arises is how to deal with the blades at the end of their lifecycle. This Master Thesis describes existing methods of disposing wind turbine blades.Moreover, the thesis considers alternative method of disposing blades, based on environmentaland safety consideration.

A simplified model is proposed to investigate the stress fields and the strain energy release rate (SERR) associated with cracks in bonded joints in wind turbine blades. The proposed two-dimensional model consists of nonparallel upper and lower shells with adhesive between the shells at the tapered end. Nonlinear finite element analysis (FEA) is performed in a systematic parametric study of material and geo- metric properties. Two failure modes and their locations are predicted at different combinations of parameters: yielding at the outside end of the adhesive and interface cracking at the inside end of the bondline. Effect of the shell curvature on the stress fields is also considered. Based on the classic beam theory and the beam-on-elastic-foundation (BOEF) assumption, stress and displacement fields of the adhesively-bonded joint were determined by a new theoretical model to support the results from the numerical computation. The failure analysis is continued by studying the effects of manufacturing defects in the adhesive bond. Single and multiple voids are embedded to simulate air bubble trapped in the interface. The numerical and analytical studies are conducted to investigate SERR associated with the voids and results are provided to illustrate the effects of void position and void size.

No description available.

A primary motivation for this research comes from the need to improve the ability of polymer-matrix composites to withstand lightning strikes. In particular, we are concerned with lightning strike damage in composite wind turbine blades. The direct effects of lightning strike on polymer-matrix composites often include rapid temperature rise, melting or burning at the lightning attachment points, and mechanical damage due to lightning-induced magnetic force and acoustic shock wave. The lightning strike damage accumulation problem is essentially multiphysic. The lightning plasma channel discharges an electric current up to 200 kA, inducing a severe heat flux at the surface of the composite structure, as well as generating Joule heating through the composite structure. The resulting electro-thermo-mechanical response of the composite structure may include matrix degradation and decomposition, delamination, and fiber breakage and sublimation, thus leading to catastrophic failure. The existing studies related to the lightning strike damage in composites ignored the lightning channel radius expansion during the initial lightning discharge and lacked adequate treatment of material phase transitions. These assumptions significantly simplify the mathematical treatment of the problem and affect the predictive capabilities of the models. Another common feature of these limited studies is that they all focused on carbon-fiber-reinforced polymer-matrix (CFRP) composites, which are electrically conductive. In the present thesis, the thermal responses and thermal ablations in a non-conductive glass-fiber-reinforced polymer-matrix (GFRP) composite wind turbine blade and in a conductive CFRP composite wind turbine blade are studied, respectively. In the case of non-conductive GFRP composite wind turbine blade, prior to the thermal response and thermal ablation analysis, a finite element analysis is performed to calculate the electric field due to lightning stepped leader to estimate the dielectric breakdown of the non-conductive composite wind turbine blade. The estimation of dielectric breakdown is used to determine whether Joule heating needs to be included in the problem formulation. To predict the thermal response and thermal ablation in the composite structure due to lightning strike, a physics-based model describing surface interaction between the lightning channel and the composite structure has been developed. The model consists of: (i) spatial and temporal evolution of the lightning channel as a function of the electric current waveform; (ii) temporary and spatially non-uniform heat flux and current density (in the case of electrically conductive CFRP composite or if dielectric breakdown occurs in the case of non-conductive GFRP composite) generated at the composite structure; and (iii) nonlinear transient heat transfer problem formulation for layered anisotropic composites that includes the moving boundary of the expanding lightning channel and the phase transition moving boundary associated with instantaneous material removal due to sublimation. The model has been employed to investigate the thermal responses and thermal ablations in a GFRP composite laminated panel used in a Sandia 100-meter all-glass baseline wind turbine blade (SNL 100-00) and a typical CFRP composite laminated panel subjected to lightning strike. The temperature-dependent directional material properties for both the GFRP and CFRP composites have been determined in this thesis using a micromechanics approach based on the experimental data for fibers and resin. An integrated Matlab-ABAQUS numerical procedure features the aforementioned aspects (i), (ii), and (iii) of the developed model. The obtained results include the evolution of temperature fields in the composite laminated panel and the progressive shape change of the composite laminated panel due to thermal ablation. The predictions of thermal ablation in the CFRP composite laminated panel are validated by reported experimental results.

Complex systems can be divided into simpler substructures. Determining the properties of each subcomponent by experimental procedures is practical and can serve to verify or calibrate finite element models. In this work, an existing model of a wind turbine blade was improved by use of experimental data. Such a blade is a subpart of a complete wind turbine. For calibration purpose, several material tests were made in order to determine the stiffness and mass properties. Later on, vibration tests of the blades were conducted and compared with simulation results of the improved model. Geometry variability within sets of blades was also studied. The blade twist angles and the center of gravity positions were found to vary moderately, which accounts for differences in blades’ dynamic behavior. Correlations between experimental data and analytical model results were very high for the first eight modeshapes. That is, according to the Model Assurance Criterion the calibrated model achieves a high-quality representation of reality. However, torsional modes in the computer model occur at a higher frequency than the experimental ones. Substructuring of the turbine allows the blades to be modeled and validated independently of the other substructures and can later be incorporated into a complete model of the turbine.

Experimental research has been completed using a custom designed and built 4m diameter wind turbine in a university operated wind facility. The primary goals of turbine testing were to determine the power production of the turbine and to apply the particle image velocimetry (PIV) technique to produce flow visualization images and velocity vector maps near the tip of a blade. These tests were completed over a wide range of wind speeds and turbine blade rotational speeds. This testing was also designed to be a preliminary study of the potential for future research using the turbine apparatus and to outline it's limitations. The goals and results of other large scale turbine tests are also briefly discussed with a comparison outlining the unique aspects of the experiment outlined in this thesis. Power production tests were completed covering a range of mean wind speeds, 6.4 m/s to 11.1 m/s nominal, and rotational rates, 40 rpm to 220 rpm. This testing allowed the total power produced by the blades to be determined as a function of input wind speed, as traditionally found in power curves for commercial turbines. The coefficient of power, Cp, was determined as a function of the tip speed ratio which gave insight into the peak power production of the experimental turbine. It was found, as expected, that the largest power production occurred at the highest input wind speed, 11.1 m/s, and reached a mean value of 3080 W at a rotational rate of 220 rpm. Peak Cp was also found, as a function of the tip speed ratio, to approach 0.4 at the maximum measurable tip speed ratio of 8. Blade element momentum (BEM) theory was also implemented as an aerodynamic power and force prediction tool for the given turbine apparatus. Comparisons between the predictions and experimental results were made with a focus on the Cp power curve to verify the accuracy of the initial model. Although the initial predictions, based on lift and drag curves found in Abbot and Von Doenhoff (1959), were similar to experimental results at high tip speed ratios an extrapolation of the data given by Hoffman et al. (1996) was found to more closely match the experimental results over the full range of tip speed ratios. Finally PIV was used to produce flow visualization images and corresponding velocity maps of the chord-wise air flow over an area at a radius ratio of 0.9, near the tip of a blade. This technique provided insight into the flow over a blade at three different tip speed ratios, 4, 6 and 8, over a range of wind speeds and rotational rates. A discussion of the unique aspects and challenges encountered using the PIV technique is presented including: measuring an unbounded external flow on a rotating object and the turbulence in the free stream affecting the uniform seeding and stability of the flow.

Experimental research has been completed using a custom designed and built 4m diameter wind turbine in a university operated wind facility. The primary goals of turbine testing were to determine the power production of the turbine and to apply the particle image velocimetry (PIV) technique to produce flow visualization images and velocity vector maps near the tip of a blade. These tests were completed over a wide range of wind speeds and turbine blade rotational speeds. This testing was also designed to be a preliminary study of the potential for future research using the turbine apparatus and to outline it's limitations. The goals and results of other large scale turbine tests are also briefly discussed with a comparison outlining the unique aspects of the experiment outlined in this thesis. Power production tests were completed covering a range of mean wind speeds, 6.4 m/s to 11.1 m/s nominal, and rotational rates, 40 rpm to 220 rpm. This testing allowed the total power produced by the blades to be determined as a function of input wind speed, as traditionally found in power curves for commercial turbines. The coefficient of power, Cp, was determined as a function of the tip speed ratio which gave insight into the peak power production of the experimental turbine. It was found, as expected, that the largest power production occurred at the highest input wind speed, 11.1 m/s, and reached a mean value of 3080 W at a rotational rate of 220 rpm. Peak Cp was also found, as a function of the tip speed ratio, to approach 0.4 at the maximum measurable tip speed ratio of 8. Blade element momentum (BEM) theory was also implemented as an aerodynamic power and force prediction tool for the given turbine apparatus. Comparisons between the predictions and experimental results were made with a focus on the Cp power curve to verify the accuracy of the initial model. Although the initial predictions, based on lift and drag curves found in Abbot and Von Doenhoff (1959), were similar to experimental results at high tip speed ratios an extrapolation of the data given by Hoffman et al. (1996) was found to more closely match the experimental results over the full range of tip speed ratios. Finally PIV was used to produce flow visualization images and corresponding velocity maps of the chord-wise air flow over an area at a radius ratio of 0.9, near the tip of a blade. This technique provided insight into the flow over a blade at three different tip speed ratios, 4, 6 and 8, over a range of wind speeds and rotational rates. A discussion of the unique aspects and challenges encountered using the PIV technique is presented including: measuring an unbounded external flow on a rotating object and the turbulence in the free stream affecting the uniform seeding and stability of the flow.

Renewable energy is a serious alternative to deliver the energy needs of an increasing world population and improve economic activity. Wind energy provides better environmental and economic benefits in comparison with the other renewable energy sources. Wind energy is capable of providing 72 TW (TW = 10^12 W) of electric power, which is approximately four and half times the world energy consumption of 15.8 TW as reported in 2006. Since power output extracted from wind turbines is proportional to the square of the blade length and the cube of the wind speed, wind turbine size has grown rapidly in the last two decades to match the increase in power output. As the blade length increases, so does its weight opening up design possibilities to introduce hybrid glass and carbon fiber composite materials as lightweight structural load bearing alternatives. Herein, we investigate the feasibility of introducing modular composite tubulars as well as hybrid sandwich composite skins in the next generation blades. After selecting a target energy output, 8 MW with 80 m blade, airfoil geometry and the layup for the skin as well as internal reinforcements are proposed. They are incorporated into the computational blade via linear shell elements for the skin, and linear beam elements for the composite tubulars to assess the relationship between weight reduction and structural performance. Computational simulations are undertaken to understand the static and dynamic regimes; specifically, displacements, stresses, and vibration modes. The results showed that the composite layers did not exhibit any damage. However, in the balsa core of the sandwich skin, the von Mises stress exceeded its allowable at wind speeds ranging from 11.0 m/sec to 12.6 m/sec. In the blades with composite tubular reinforcement, two different types of damage are observed: a. Stress concentrations at the tubular-skin attachments, and b. Highest von Mises stress caused by the flapping bending moment. The vibration studies revealed a strong coupling mode, bending and twist, at the higher natural frequencies of the blade with tubular truss configuration. The weight saving measures in developing lighter blades in this study did not detract from the blades structural response for the selected load cases.

No description available.

No description available.