Dynamic Analysis of an Offshore Wind Turbine Including Soil Effects Procedia Engineering 116 ( 2015 ) 32 – 39 1877-7058 © 2015 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/). Peer- Review under responsibility of organizing committee , IIT Madras , and International Steering Committee of APAC 2015 doi: 10.1016/j.proeng.2015.08.261 ScienceDirect Available online at www.sciencedirect.com 8th International Conference on Asian and Pacific Coasts (APAC 2015) Department of Ocean Engineering, IIT Madras, India. Dynamic Analysis of an Offshore Wind Turbine Including Soil Effects K.A. Abhinav, Nilanjan Saha* aDepartment of Ocean Engineering, IIT Madras, Chennai 600036, India Abstract Offshore wind turbines (OWTs) offer an attractive, sustainable solution to the impending global energy crisis. A major challenge in fixed-bottom OWT design is accounting for soil-structure interaction (SSI) under the influence of random dynamic loading from wind, waves and currents. Usually, SSI is either ignored in OWT studies or is incorporated by means of simplified foundation concepts like the apparent fixity model. OWTs in shallow water depths (less than 30 m) are mostly supported on monopiles - large diameter steel pipe piles driven into the subsoil. Monopiles transfer the dynamic lateral loads into the soil by bending action. The present work deals with the dynamic analysis of the NREL 5MW OWT on a monopile foundation, in Indian waters. It involves parametric studies on various clayey soil profiles - soft, medium stiff and stiff clay. An operational wind speed of 12 m/s and a sea state of 4 m significant wave height and 10 s spectral peak period are considered. The OWT design should ensure that the natural frequency is away from the forcing frequencies of wind, wave and rotor. A water depth of 20 m is considered. Hub-height aerodynamic loads are obtained using the NREL-FAST code, which is based on the blade- element momentum (BEM) theory. The hydrodynamic time domain analyses are performed in the FEM based coupled hydrodynamic - geotechnical software, DNV-GL - USFOS. USFOS makes use of the JONSWAP spectrum to generate irregular waves. Soil is represented by means of p-y, Q-z and t-z curves. Results indicate the significance of including SSI in OWT studies. Variation in response due to change in pile penetration depth and pile diameter are also highlighted. Stiffness of clay is the design driver for OWTs. © 2014 The Authors. Published by Elsevier B.V. Peer-review under responsibility of organizing committee, IIT Madras, and International Steering Committee of APAC 2015 * Corresponding author. Tel.:+91-44-22574827; fax:+91-44-22574802. E-mail address:nilanjan@iitm.ac.in © 2015 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/). Peer- Review under responsibility of organizing committee , IIT Madras , and International Steering Committee of APAC 2015 http://crossmark.crossref.org/dialog/?doi=10.1016/j.proeng.2015.08.261&domain=pdf 33 K.A. Abhinav and Nilanjan Saha / Procedia Engineering 116 ( 2015 ) 32 – 39 Keywords:offshore wind turbine, soil-structure interaction, monopile 1. Introduction Offshore wind turbines are gaining popularity as the need to offset the nonrenewable energy crisis is being recognized. OWTs in shallow waters are usually supported on monopile substructures, which are essentially large diameter piles in the range of 3-6 m. Monopiles transfer lateral loads from wind and waves by bending action, into the soil. Unlike oil and gas structures, OWT systems are subjected to heavy lateral loads, often to the tune of 150% of the vertical loading (Lesny, 2010). This calls for a comprehensive dynamic analysis, especially in soft soils, where the natural frequency of the system could undergo significant variations (Bazeos et al., 2002). Laterally loaded piles are designed using the p y method specified by the API (2000). However, these soil curves are derived from experiments involving small diameter piles and their validity for monopile design is still being debated (Achmus and Abdel-Rahman, 2012). Lesny and Wiemann (2005), Haiderali et al. (2013), Achmus et al. (2009) etc. have made use of quasi-static mudline loads to investigate the response of monopiles. Abhinav and Saha (2015) analysed for a jacket subjected to random wave loads considering soil structure interaction while, Gao et al. (2010) showed that jacket structure can be analysed as a set of monopiles of varying diameters. ABS (2011) followed a quasi-aerostatic-hydrodynamic approach to study OWT response under different water depths. Bisoi and Haldar (2014) analyzed a 2MW OWT under dynamic loads, including soil-structure interaction (SSI) Here, the aerodynamic load at the hub was modelled as a function of the rotor frequency. The present work deals with the analysis of a monopile support structure for the NREL 5MW baseline OWT, under aerodynamic and hydrodynamic loads, for varying soil conditions. The FEM-based code USFOS (SINTEF Group, 2001) is used for hydrodynamic analysis. USFOS implements a nonlinear FE formulation based on Green strain E (valid for all magnitudes of rotation and displacement), where a finite element represents an actual element of the structure. This pre-empts the choice of mesh-size and speeds up computation. The HHT-α method is used for dynamic analysis. Aerodynamic loads are derived using an aeroelastic code, FAST (Jonkman and Buhl, 2005). Three different clayey soil profiles - soft, medium stiff and stiff clay - are considered. The effects of soil stiffness variation, pile penetration depth and pile diameter, on the lateral response and natural frequency of the OWT structure are investigated. Soil stiffness is derived using the soil curves specified by API (2000). Results show significant variations in the response and point towards the need for detailed modelling of soil-structure interaction in the analysis of OWT structures. 2. Numerical modeling 2.1. Structural model The NREL 5MW baseline OWT, conceptualized by Jonkman et al. (2009), is considered in the present study. The numerical model of a monopile supporting the NREL 5MW OWT, in a water depth of 20 m is developed in USFOS. The turbine rotor-nacelle-assembly (RNA) is supported by a steel tower, connected to the monopile through a cylindrical transition piece. The diameter and thickness of the tower varies from 6m and 0.027 m respectively, at the base, to 3.87 m and 0.019 m respectively, at the top. The components of the OWT structure and the corresponding design levels are represented in Figure 1. A higher density value is considered for steel (8500 kg/m3) to account for the absence of bolts, flanges and welds in the model. The monopile has a diameter of 6 m and is founded on uniform stiff clay. Two other soil profiles (medium-stiff clay and soft clay) are also considered for the analysis, to investigate the response of the structure under dynamic loads. The soil profiles are defined in Table 1. Here, γ is the submerged unit weight, Su is the undrained shear strength and 50ε is the strain at 50% of maximum stress in undrained compression test. The depth of embedment is taken as 42m or 7 times the pile diameter. The wall thickness ( t ) of the monopile is defined by API (2000), on the basis of its diameter ( D ), as shown in equation 1. 34 K.A. Abhinav and Nilanjan Saha / Procedia Engineering 116 ( 2015 ) 32 – 39 6.35 100 D t (1) Fig. 1. OWT conceptual and FE model Table 1. Soil profiles Soil type γ (kN/m3) Su (kPa) 50ε Soft clay 8 40 0.020 Medium stiff clay 8 80 0.010 Stiff clay 8 120 0.005 The tower and piles are modelled using 2-noded beam elements. All analyses are performed for an operational hub-height wind speed of 12 m/s and a sea-state defined by a significant wave height ( SH ) of 4 m and wave spectral peak period ( PT ) of 10 s. 2.2. Modelling of soil-structure-interaction Soil-structure interaction is modelled in USFOS using lateral and axial nonlinear spring-to-ground elements along the length of the pile. The properties of the spring elements are derived from the p y , t z , and Q z , curves recommended by API (2000). Pile-soil interaction in the lateral direction is accounted for, by the lateral soil resistance - deflection or p y curves, while that in the vertical direction is modelled using the curves for skin friction ( t z ) and end bearing ( Q z ). More information on the development of cyclic p y curves for soft and stiff clays is described in Reese and van Impe (2011). 35 K.A. Abhinav and Nilanjan Saha / Procedia Engineering 116 ( 2015 ) 32 – 39 2.3. Modelling parameters The OWT is located in a water depth of 20 m in clayey soil. The hub height wind speed is 12 m/s. Simulations are performed for duration of 600 s. A time step value of 0.01 s has been used for dynamic analysis. Soil springs are attached to the pile at a distance of 2 m. A frequency value of 300 is used for the discretizing the JONSWAP spectrum by constant area method. 3. Description of loads OWT are subjected to the action of aerodynamic and hydrodynamic loading. USFOS has coupled hydrodynamic and geotechnical modules. However, the absence of an aerodynamic component requires the use of a separate code for the same. According to Seidel et al. (2009) and Gao et al. (2010), the response of an OWT can be accurately realized using decoupled analyses, if its natural period varies significantly from the predominant wave period. The fundamental vibration mode of the OWT involves the vibration of the tower, indicating the lack of a significant contribution of wave loading, to this mode. This is consistent with the findings in Gao et al. (2010) and points to a lack of strong coupling between wind and wave loads. Hence, a decoupled analysis could be used to model the response of the OWT. The present work ignores the influence of wind loading on the tower of the OWT. The modeling concepts for dynamic loading are discussed below. 3.1. Aerodynamic loads Aerodynamic loads on the NREL 5MW OWT are computed using the FAST code. NREL's TurbSim code (Jonkman, 2009) is used to generate stochastic wind fields corresponding to a hub-height wind speed of 12 m/s. The Kaimal spectrum (Kaimal et al., 1972) is used in conjunction with the Normal Turbulence Model (IEC, 2009) to realize 3-component wind speed vectors encompassing the entire rotor plane. FAST makes use of the blade element momentum theory (Moriarty and Hansen, 2005) to generate time-series of aerodynamic load components at the hub of the OWT. 3.2. Hydrodynamic loads Wave loading on slender structures like piles are calculated using the Morison equation (Chakrabarti, 2005). Here, the total load is obtained as the summation of drag and inertia components, arising from the velocity and acceleration of water particles, respectively. Morison equation gives the hydrodynamic force (f) per unit length of the vertical pile as: 2 1 4 2 . M D D f C C u uu π ρ ρ (2) where D is the diameter of the pile, u is the horizontal water particle velocity, . u is the water particle acceleration, ρ is the density of water and ,M DC C are hydrodynamic coefficients of inertia and drag, respectively. The dot represents acceleration. Long-crested irregular waves for dynamic analyses are modelled using the constant area discretization of the Joint North Sea Wave Project (JONSWAP) spectrum. The JONSWAP spectrum is valid for fetch-limited seas and homogenous wind conditions. 4. Results and Discussion The generation of time series for stochastic wind and wave loading involves the use of pseudo-random number generators for realizing the random phases (Zwick and Muskulus, 2014). The use of different random seeds for the 36 K.A. Abhinav and Nilanjan Saha / Procedia Engineering 116 ( 2015 ) 32 – 39 same wind (or wave) data results in different time-series and hence, statistical uncertainty. This has to be taken care of by multiple simulations, involving different seed numbers. The results present herein are mean values from an ensemble size of 50, involving different seed numbers for generating wind and wave time series. 4.1. Effect of soil variability on natural frequency of OWT As stated in Bhattacharya et al. (2013), a 3-bladed OWT is excited at the following: (a) the rotor frequency (1P), (b) the blade pass frequency (3P), (c) the wave loads and (d) the gusty winds. Economic considerations dictate that the frequency of the OWT lies in between the 1P and the 3P regions. This is termed as the soft-stiff design, which corresponds to a frequency band of 0.20 Hz to 0.35 Hz for the NREL 5MW-OWT. It is observed, from Table 2, that softening soils tends to increase the period of vibration and hence, result in a reduction in natural frequency. In the case of soft-clay, the natural frequency was found to dip into the unsafe 1P region, where the OWT structure becomes vulnerable to resonance. Pile penetration depth of 42 m and diameter of 6 m is considered for analysis. Table 2. Variation of natural frequency Soil type Natural frequency (Hz) Soft clay 0.16 Medium stiff clay 0.21 Stiff clay 0.24 4.2. Effect of soil variability on lateral response of OWT For the considered pile penetration depth of 42 m, medium stiff and stiff clayey soil profiles show a variation of 60% in the lateral displacement at the mudline and 15% at the tower-top, respectively, as shown in Figure 2. The displacement of the OWT in soft clay is exceeds the serviceability limit of 0.2 m at mudline. 37 K.A. Abhinav and Nilanjan Saha / Procedia Engineering 116 ( 2015 ) 32 – 39 Fig. 2. Effect of soil variability Fig. 3. Effect of pile diameter variation 4.3. Effect of pile diameter on lateral response of OWT Three different pile diameters were considered for a penetration depth of 42 m, in stiff clay - 7 m, 6 m and 5 m. The thickness of the piles is defined according to Equation 1. It is observed from Figure 3, that an increase in pile diameter brings about a reduction in the lateral response of the OWT. This is due to the fact that the stiffness obtained from a p y curve is proportional to the diameter of the pile. 4.4. Effect of pile depth on lateral response of OWT Stiff clay exhibits the same displacement profile for pile penetration depths from 42 m to 30 m. Medium stiff clay shows an increase in displacement with reduction in pile penetration depth. These results are shown in Figures 4a and 4b. A marked difference in the displacement of the embedded portion of the pile is also observed. For smaller penetration depths, the response of the OWT in soft clay escalates and this results in failure of the system. 38 K.A. Abhinav and Nilanjan Saha / Procedia Engineering 116 ( 2015 ) 32 – 39 Fig. 4a. Medium stiff clay Fig. 4b. Stiff clay 5. Conclusions Decoupled aero-hydrodynamic simulations on a monopile in clay of varying stiffness underscore the importance of considering SSI in OWT studies. Softer soils increased the time-period of vibration and hence the natural frequency of the OWT shifted to the resonant region. Soil stiffness is the determining factor in the performance of OWTs in clayey soils. Soft clays were found to produce excessive motions that transcend the serviceability limit state, leading to failure. Stiff clays, on the other hand, produced relatively constant response with varying pile depth and diameter. Similar analysis is required to be extended for sandy soils as well as mixed soils, which would form a part of the future work. Acknowledgements The authors would like to thank Dr. Tore Holmas for his help with USFOS. The authors gratefully acknowledge the financial support given by the Earth Systems Science Organization, Ministry of Earth Sciences, Government of India through National Institute of Ocean Technology to conduct part of the research. References ABS, 2011. Design Standards for Offshore Wind Farms, American Bureau of Shipping, Houston, USA Abhinav KA and Saha, N., 2015. Coupled hydrodynamic and geotechnical analysis of jacket offshore wind turbine, Soil Dyanmics and Earthquake Engineering, 73, pp.66-79. Achmus, M., Kuo, Y-S., Abdel-Rahman, K., 2009. Behavior of monopile foundations under cyclic lateral load, Computers and Geotechnics, 36 725–735 Achmus, M., Abdel-Rahman, K.,.2012. Design of Piles for Offshore Wind Energy Foundations with Respect to Horizontal Loading, Proceedings of the Twenty-second International Offshore and Polar Engineering Conference, ISOPE, Rhodes, Greece, 143-150 39 K.A. Abhinav and Nilanjan Saha / Procedia Engineering 116 ( 2015 ) 32 – 39 API, 2000. API-RP-2A - Recommended Practice for Planning, Designing and Constructing Fixed Offshore Platforms - Working Stress Design, American Petroleum Institute, Washington, D.C., USA Bazeos, N., Hatzigeorgiou, G.D., Hondros, I.D., Karamaneas, H., Karabalis, D.L., Beskos, D.E., 2002. Static, seismic and stability analysis of a prototype wind turbine steel tower, Engineering Structures, 24, 1015–1025 Bhattacharya, S., Nikitas, N., Garnsey, J., Alexander, N.A., Cox, J., Lombardi, D., Muir Wood, D., Nash, D.F.T., 2013. Observed dynamic soil–structure interaction in scale testing of offshore wind turbine foundations, Soil Dynamics and Earthquake Engineering, 54, 47-60 Bisoi, S., Haldar, S., 2014. Dynamic analysis of offshore wind turbine in clay considering soil–monopile–tower interaction, Soil Dynamics and Earthquake Engineering, 63, 19 - 35 Chakrabarti, S.K., 2005. Handbook of Offshore Engineering - Volume -I, Elsevier, Amsterdam, The Netherlands Gao, Z., Saha, N., Moan, T. , Amdahl, J., 2010. Dynamic analysis of offshore fixed wind turbines under wind and wave loads using alternative computer codes, Proceedings of the 3rd EAWE Conference, TORQUE 2010: The Science of Making Torque from Wind, Crete, Greece Haiderali, A., Cilingir, U., Madabhushi, G., 2013. Lateral and Axial Capacity of Monopiles for Offshore Wind Turbines, Indian Geotechnical Journal, 43(3):181-194 IEC., 2009. IEC 61400-3 - Wind turbines – Part 3: Design requirements for offshore wind turbines, INTERNATIONAL ELECTROTECHNICAL COMMISSION, Geneva, Switzerland Jonkman, B.J., 2009. TurbSim User's Guide v 1.50, Technical Report NREL/TP-500-46198, National Renewable Energy Laboratory, Golden, CO, USA Jonkman, J.M., Buhl, M.L. Jr., 2005. FAST User's Guide, Technical Report NREL/EL-500-38230, National Renewable Energy Laboratory, Golden, CO, USA Jonkman, J.M., Butterfield, S., Musial, W., Scott, G., 2009. Definition of a 5-MW Reference Wind Turbine for Offshore System Development. Technical Report NREL/TP-500-38060, National Renewable Energy Laboratory, Golden, CO, USA Kaimal, J.C., Wyngaard, J.C., Izumi, Y., Cote, O.R., 1972. Spectral characteristics of surface-layer turbulence, Quarterly Journal of the Royal Meteorological Society, 98, 563–589 Lesny, K., Wiemann, J., 2005. Design aspects of monopiles in German offshore wind farms, In: Frontiers of Offshore Geotechnics, Eds. Gourvenec, S., Cassidy., M, Taylor & Francis Group, London, pp 383-389 Lesny, L., 2010. Foundations for Offshore Wind Turbines, VGE Verlag GmbH, Essen, Germany Moriarty, P. J., Hansen, A.C., 2005. AeroDyn Theory Manual, Technical Report NREL/EL-500-36881, National Renewable Energy Laboratory, Golden, CO, USA Reese, L.C., van Impe, W.F., 2011. Single Piles and Pile Groups Under Lateral Loading, CRC Press/Balkema, Leiden, The Netherlands Seidel, M., Ostermann, F., Curvers, A.P.W.M., Kuhn, M., Kaufer, D., Boker, C., 2009. Validation of offshore load simulations using measurement data from the DOWNVInD project. Conference proceedings of European Offshore Wind, Stockholm, Sweden SINTEF Group, 2001. USFOS-Getting Started, Structural Engineering, Marintek Zwick, D., Muskulus, M., 2014. The simulation error caused by input loading variability in offshore wind turbine structural analysis, Wind Energy (available online) DOI: 10.1002/we.1767