<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing with OASIS Tables v3.0 20080202//EN" "journalpub-oasis3.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:oasis="http://docs.oasis-open.org/ns/oasis-exchange/table" xml:lang="en" dtd-version="3.0" article-type="research-article">
  <front>
    <journal-meta><journal-id journal-id-type="publisher">WES</journal-id><journal-title-group>
    <journal-title>Wind Energy Science</journal-title>
    <abbrev-journal-title abbrev-type="publisher">WES</abbrev-journal-title><abbrev-journal-title abbrev-type="nlm-ta">Wind Energ. Sci.</abbrev-journal-title>
  </journal-title-group><issn pub-type="epub">2366-7451</issn><publisher>
    <publisher-name>Copernicus Publications</publisher-name>
    <publisher-loc>Göttingen, Germany</publisher-loc>
  </publisher></journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.5194/wes-8-575-2023</article-id><title-group><article-title>Prognostics-based adaptive control strategy  for lifetime control of wind turbines</article-title><alt-title>Prognostics-based adaptive control strategy  for lifetime control of wind turbines</alt-title>
      </title-group><?xmltex \runningtitle{Prognostics-based adaptive control strategy  for lifetime control of wind turbines}?><?xmltex \runningauthor{E. Kipchirchir et al.}?>
      <contrib-group>
        <contrib contrib-type="author" corresp="yes" rid="aff1">
          <name><surname>Kipchirchir</surname><given-names>Edwin</given-names></name>
          <email>edwin.kipchirchir@uni-due.de</email>
        <ext-link>https://orcid.org/0000-0002-5567-7225</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff2">
          <name><surname>Do</surname><given-names>M. Hung</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff3">
          <name><surname>Njiri</surname><given-names>Jackson G.</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Söffker</surname><given-names>Dirk</given-names></name>
          
        </contrib>
        <aff id="aff1"><label>1</label><institution>Chair of Dynamics and Control, University of Duisburg-Essen, <?xmltex \hack{\break}?>Lotharstr. 1-21, 47057, Duisburg, Germany</institution>
        </aff>
        <aff id="aff2"><label>2</label><institution>School of Heat Engineering and Refrigeration, Hanoi University of Science and Technology, <?xmltex \hack{\break}?> No. 1, Dai Co Viet, Hai Ba Trung, Hanoi, Vietnam</institution>
        </aff>
        <aff id="aff3"><label>3</label><institution>Department of Mechatronic Engineering, Jomo Kenyatta University of Agriculture and Technology, 62000-00200, Nairobi, Kenya</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">Edwin Kipchirchir (edwin.kipchirchir@uni-due.de)</corresp></author-notes><pub-date><day>24</day><month>April</month><year>2023</year></pub-date>
      
      <volume>8</volume>
      <issue>4</issue>
      <fpage>575</fpage><lpage>588</lpage>
      <history>
        <date date-type="received"><day>23</day><month>November</month><year>2021</year></date>
           <date date-type="rev-request"><day>30</day><month>November</month><year>2021</year></date>
           <date date-type="rev-recd"><day>10</day><month>February</month><year>2023</year></date>
           <date date-type="accepted"><day>24</day><month>March</month><year>2023</year></date>
      </history>
      <permissions>
        <copyright-statement>Copyright: © 2023 Edwin Kipchirchir et al.</copyright-statement>
        <copyright-year>2023</copyright-year>
      <license license-type="open-access"><license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><self-uri xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023.html">This article is available from https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023.html</self-uri><self-uri xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023.pdf">The full text article is available as a PDF file from https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023.pdf</self-uri>
      <abstract><title>Abstract</title>

      <p id="d1e126">Variability in wind profiles in both space and time is responsible for fatigue loading in wind turbine components. Advanced control methods for mitigating structural loading in these components have been proposed in previous works. These also incorporate other objectives like speed and power regulation for above-rated wind speed operation. In recent years, lifetime control and extension strategies have been proposed to guarantee power supply and operational reliability of wind turbines. These control strategies typically rely on a fatigue load evaluation criteria to determine the consumed lifetime of these components, subsequently varying the control set point to guarantee a desired lifetime of the components. Most of these methods focus on controlling the lifetime of specific structural components of a wind turbine, typically the rotor blade or tower. Additionally, controllers are often designed to be valid about specific operating points and hence exhibit deteriorating performance in varying operating conditions. Therefore, they are not able to guarantee a desired lifetime in varying wind conditions. In this paper an adaptive lifetime control strategy is proposed for controlled aging of rotor blades to guarantee a desired lifetime while considering damage accumulation level in the tower. The method relies on an online structural health monitoring system to vary the lifetime controller gains based on a state-of-health (SoH) measure by considering the desired lifetime at every time step. For demonstration, a 1.5 MW National Renewable Energy Laboratory (NREL) reference wind turbine is used. The proposed adaptive lifetime controller regulates structural loading in the rotor blades to guarantee a predefined damage level at the desired lifetime without sacrificing the speed regulation performance of the wind turbine. Additionally, a significant reduction in the tower fatigue damage is observed.</p>
  </abstract>
    </article-meta>
  </front>
<body>
      

<sec id="Ch1.S1" sec-type="intro">
  <label>1</label><title>Introduction</title>
      <p id="d1e138">Growing demand for wind energy has led to the development of large wind turbines. However, these turbines are less tolerant to system performance degradation and faults <xref ref-type="bibr" rid="bib1.bibx13" id="paren.1"/>. To ensure utility-scale wind turbines operate with respect to their design lifetime, advanced control strategies have been developed in recent years to reduce structural loading of blades and the tower. Most of these incorporate additional objectives such as power optimization and rotor speed regulation. The objective of the lifetime control of wind turbines using prognostics-based load mitigation strategies has become more important in recent years. Most of the proposed methods focus on controlling the lifetime of one structural component of a wind turbine, typically the rotor blade or the tower, without considering the fatigue damage level in other components. These lifetime controllers are also designed to be valid regarding specific operating points.</p>
      <?pagebreak page576?><p id="d1e144">A control strategy for extending the maintenance interval of wind turbine blades under assumed crack initiation, detected using a data filtering algorithm, is proposed <xref ref-type="bibr" rid="bib1.bibx6" id="paren.2"/>. In <xref ref-type="bibr" rid="bib1.bibx7" id="text.3"/> and <xref ref-type="bibr" rid="bib1.bibx23" id="text.4"/>, a set of switching controllers with varying degrees of load mitigation are engaged sequentially based on the accumulated damage obtained from an online damage evaluation model to extend the lifetime of rotor blades. An adaptive lifetime controller is proposed in <xref ref-type="bibr" rid="bib1.bibx9" id="text.5"/> to guarantee the desired lifetime of the tower. Depending on the damage accumulation and the predicted lifetime provided by an online damage evaluation model, the weights of the lifetime controller are varied. However, in <xref ref-type="bibr" rid="bib1.bibx6" id="text.6"/>, <xref ref-type="bibr" rid="bib1.bibx7" id="text.7"/>, <xref ref-type="bibr" rid="bib1.bibx23" id="text.8"/>, and <xref ref-type="bibr" rid="bib1.bibx9" id="text.9"/>, fatigue damage is considered in only one turbine component. The lifetime controllers used are not adaptive to varying wind conditions. In recent times, resilient control has been proposed in <xref ref-type="bibr" rid="bib1.bibx1" id="text.10"/>, <xref ref-type="bibr" rid="bib1.bibx4" id="text.11"/>, <xref ref-type="bibr" rid="bib1.bibx12" id="text.12"/>, and <xref ref-type="bibr" rid="bib1.bibx15" id="text.13"/> to minimize the effect of unanticipated faults or unexpected dynamics to maintain the operation of a wind turbine within a limited degradation tolerance bound. However, resilient control does not address the problem of controlling life consumption in wind turbine components to avoid early fatigue failures. Although new concepts like operational modal analysis (OMA), which relies on measurement data to analyze vibrating structures, are becoming the industry standard for condition monitoring and diagnosis especially for offshore wind turbines <xref ref-type="bibr" rid="bib1.bibx18 bib1.bibx5 bib1.bibx11 bib1.bibx24" id="paren.14"/>, these concepts have yet to be integrated for prognosis and lifetime control of wind turbines.</p>
      <p id="d1e188">In this work an adaptive lifetime control strategy is proposed for controlling the aging of rotor blades to guarantee a desired lifetime while considering damage accumulation level in the tower. A robust disturbance accommodating control (RDAC) proposed in <xref ref-type="bibr" rid="bib1.bibx10" id="text.15"/> is used for rotor speed regulation and load mitigation in the tower, while a prognostics-based adaptive independent pitch control (aIPC), which adapts to wind speed variation, is used for lifetime control of rotor blades. By monitoring the accumulated damage using an online structural health evaluation model, the load mitigation level in the blades is controlled by varying the control gains in the respective IPC controllers based on a threshold evaluation of the estimated lifetime. As an improvement to the approaches in the aforementioned contributions, the proposed adaptive lifetime control strategy regulates fatigue loading in the rotor blades to reach a predefined damage limit at the desired lifetime with subsequent reduction in tower damage accumulation. This is realized without trade-off in speed and power regulation performance. Although the proposed lifetime control strategy is applied to wind turbine operations in above-rated operation, this concept can be extended for lifetime control in below-rated operation with suitable objectives including reducing structural loads while ensuring maximum power extraction as proposed in <xref ref-type="bibr" rid="bib1.bibx9" id="text.16"/>.</p>
      <p id="d1e197">The paper is organized as follows. In Sect. 2, a theoretical background on wind turbine health monitoring is given. In Sect. 3, the design of the primary RDAC for rotor speed regulation and tower load mitigation, as well as the prognostics-based aIPC lifetime controller for controlled aging of rotor blades, is outlined. The proposed prognostics-based adaptive lifetime control strategy, which incorporates the primary and lifetime controllers, as well as an online damage evaluation model, is described in Sect. 4. In Sect. 5, simulation results based on the performance evaluation of the proposed prognostics-based adaptive lifetime control strategy on a reference wind turbine are discussed. Lastly, a summary and conclusions are given in Sect. 6.</p>
</sec>
<sec id="Ch1.S2">
  <label>2</label><title>Wind turbine health monitoring</title>
      <p id="d1e208">Wind speed variability subjects wind turbine components like blades and the tower to cyclic loading. This causes damage to be accumulated in these components over time causing gradual degradation until failure occurs. Therefore, structural health monitoring of wind turbines is important in preventing the occurrence of fatigue failure before reaching the related design lifetime. Information on the damage evolution in a component can be utilized as a health indicator for failure detection, as well as for developing control measures to guarantee the desired lifetime. This section outlines the methods used for estimating the damage accumulation in wind turbine components.</p>
<sec id="Ch1.S2.SS1">
  <label>2.1</label><title>Evaluation of damage accumulation</title>
      <?pagebreak page577?><p id="d1e218">A wind turbine endures varying and complex load conditions over its lifetime. Fatigue analysis is therefore important in determining the consumed lifetime of its components. Component degradation starts at microscale as micro-cracks resulting from irreversible changes in the microstructure and propagates gradually until it fails. Assumptions of underlying damage evolution laws are often used to estimate the actual damage level, as well as to predict the remaining useful life (RUL) of a component. Component-specific high-cycle fatigue experiments are used to generate  <inline-formula><mml:math id="M1" display="inline"><mml:mi>S</mml:mi></mml:math></inline-formula>–<inline-formula><mml:math id="M2" display="inline"><mml:mi>N</mml:mi></mml:math></inline-formula> curves (Wöhler curve), which describe the relationship between applied stress amplitude <inline-formula><mml:math id="M3" display="inline"><mml:mi>S</mml:mi></mml:math></inline-formula> and the number of load cycles <inline-formula><mml:math id="M4" display="inline"><mml:mi>N</mml:mi></mml:math></inline-formula> that would cause failure. This forms the basis for the mathematical relation for fatigue analysis in wind turbine components expressed as
            <disp-formula id="Ch1.E1" content-type="numbered"><label>1</label><mml:math id="M5" display="block"><mml:mrow><mml:msup><mml:mi>s</mml:mi><mml:mi>m</mml:mi></mml:msup><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mi>K</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
          where <inline-formula><mml:math id="M6" display="inline"><mml:mi>s</mml:mi></mml:math></inline-formula> denotes the stress range amplitude and <inline-formula><mml:math id="M7" display="inline"><mml:mi>m</mml:mi></mml:math></inline-formula> the Wöhler exponent (typically 3 for steel materials like the tower and 10 for composites like the blade; <xref ref-type="bibr" rid="bib1.bibx25" id="altparen.17"/>). The material parameter of fatigue damage at failure <inline-formula><mml:math id="M8" display="inline"><mml:mi>K</mml:mi></mml:math></inline-formula>, e.g., ultimate tensile strength, is related to the number of load cycles <inline-formula><mml:math id="M9" display="inline"><mml:mi>N</mml:mi></mml:math></inline-formula>.</p>
      <p id="d1e300">Wind turbine components are designed for a service lifetime of at least 20 years according to the International Electrotechnical Commission (IEC) standard, with these structural components facing roughly between <inline-formula><mml:math id="M10" display="inline"><mml:mrow><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mn mathvariant="normal">8</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M11" display="inline"><mml:mrow><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mn mathvariant="normal">9</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> fatigue load cycles <xref ref-type="bibr" rid="bib1.bibx30" id="paren.18"/>. The component lifetime is typically arrived at using the projected number of fatigue cycles and average wind conditions it will encounter in its lifetime. Additionally, the IEC standard specifies that a wind turbine component should be designed to maintain its structural integrity in case it experiences 50-year extreme wind events during its lifetime.</p>
      <p id="d1e328">Fatigue damage in components can be assessed using linear damage accumulation theory based on Miner's rule or nonlinear damage accumulation theories <xref ref-type="bibr" rid="bib1.bibx29" id="paren.19"/>. Due to its simplicity, Miner's rule <xref ref-type="bibr" rid="bib1.bibx20" id="paren.20"/> is widely used. Wind speed variability induces a varying-amplitude load spectrum on wind turbine components. To use Miner's rule, the complex spectrum of varying load is often transformed using the rain-flow counting (RFC) algorithm first proposed by <xref ref-type="bibr" rid="bib1.bibx19" id="text.21"/> into simple uniform loading, from which stress range histograms can be extracted and used to assess the accumulated damage. A schematic of this procedure is shown in Fig. <xref ref-type="fig" rid="Ch1.F1"/>.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F1"><?xmltex \currentcnt{1}?><?xmltex \def\figurename{Figure}?><label>Figure 1</label><caption><p id="d1e345">Conventional fatigue load estimation <xref ref-type="bibr" rid="bib1.bibx25" id="paren.22"/>.</p></caption>
          <?xmltex \igopts{width=236.157874pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f01.png"/>

        </fig>

      <p id="d1e357"><?xmltex \hack{\newpage}?>By combining RFC and Miner's rule, damage accumulation <inline-formula><mml:math id="M12" display="inline"><mml:mrow><mml:msub><mml:mi>D</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is calculated as
            <disp-formula id="Ch1.E2" content-type="numbered"><label>2</label><mml:math id="M13" display="block"><mml:mrow><mml:msub><mml:mi>D</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:msub><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>n</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>n</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msubsup><mml:mi>s</mml:mi><mml:mi>i</mml:mi><mml:mi>m</mml:mi></mml:msubsup></mml:mrow><mml:mi>K</mml:mi></mml:mfrac></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
          where <inline-formula><mml:math id="M14" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula> denotes the total number of related stress range histograms, <inline-formula><mml:math id="M15" display="inline"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> the incremental damage at the <inline-formula><mml:math id="M16" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula>th stress range histogram, <inline-formula><mml:math id="M17" display="inline"><mml:mrow><mml:msub><mml:mi>n</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> the number of applied load cycles in each histogram bin, <inline-formula><mml:math id="M18" display="inline"><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> the number of cycles to failure at the <inline-formula><mml:math id="M19" display="inline"><mml:mi>i</mml:mi></mml:math></inline-formula>th stress range histogram, and <inline-formula><mml:math id="M20" display="inline"><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> the applied load amplitude in each histogram bin. With continual load application, damage in a component progresses from an undamaged state <inline-formula><mml:math id="M21" display="inline"><mml:mrow><mml:msub><mml:mi>D</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> to the point it is considered to have reached its end of life when the accumulated damage <inline-formula><mml:math id="M22" display="inline"><mml:mrow><mml:msub><mml:mi>D</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>. In this case, the component is considered to have exhausted its structural reserves. Although other cycle counting algorithms including level crossing counting, peak counting, and simple range counting exist, RFC algorithms are the most widely applied for fatigue analysis <xref ref-type="bibr" rid="bib1.bibx21" id="paren.23"/>.</p>
</sec>
<sec id="Ch1.S2.SS2">
  <label>2.2</label><title>Online rain-flow counting</title>
      <p id="d1e571">Most standard RFC algorithms generate equivalent load cycles from complex load spectra by pairing local minimum and maximum points using the three-point counting rule. Therefore, the entire load history is needed beforehand for the equivalent cycles to be generated. This process is computationally inefficient because the algorithm has to process all the stored loading data. Therefore, standard RFC  cannot be used for real-time monitoring or control of life consumption of a component  <xref ref-type="bibr" rid="bib1.bibx21" id="paren.24"/>.</p>
      <p id="d1e577">In <xref ref-type="bibr" rid="bib1.bibx21" id="text.25"/>, a real-time implementation of the RFC algorithm is proposed. By employing a three-point counting rule recursively, the extremal points of time-series loading data are processed and stored in two flexible stacks as they occur to pick out the full and half cycles. For each identified cycle and using Miner's rule, the life consumption of a component is then calculated and incremented online. This allows for the online determination of the consumed life of a component, as well as implementation of lifetime control. In this paper, the online damage evaluation algorithm <xref ref-type="bibr" rid="bib1.bibx21" id="paren.26"/> is adopted for evaluating the accumulated damage in rotor blades and the tower. This information is then used to adapt the lifetime controller to guarantee a predefined service life of the wind turbine components.</p>
</sec>
</sec>
<sec id="Ch1.S3">
  <label>3</label><title>Control strategy for load mitigation and speed regulation</title>
      <p id="d1e595">In this paper, the robust disturbance accommodating controller (RDAC) <xref ref-type="bibr" rid="bib1.bibx10" id="paren.27"/>, proposed for rotor speed regulation and mitigation of tower fore–aft bending moments, is extended to include an adaptive independent<?pagebreak page578?> pitch controller (aIPC), which is used as a dynamic lifetime controller for reducing blade flapwise bending moments in a wind turbine operating in an above-rated wind speed region. In this section the description of the reference wind turbine (RWT) is outlined. Additionally, the description of the adaptive robust observer-based controller, which is adapted for lifetime control, is summarized.</p>
<sec id="Ch1.S3.SS1">
  <label>3.1</label><title>Wind turbine model description</title>
      <p id="d1e608">A 1.5 MW WindPACT reference wind turbine developed by the National Renewable Energy Laboratory (NREL) <xref ref-type="bibr" rid="bib1.bibx26" id="paren.28"/>, which is available in fatigue, aerodynamics, structures, and turbulence (FAST) design code <xref ref-type="bibr" rid="bib1.bibx17" id="paren.29"/>, is chosen as the test bed for the design and simulation of the proposed adaptive lifetime control strategy. This onshore wind turbine model was developed based on a real-life commercial wind turbine used in the WindPACT program. The specifications of this turbine are summarized in Table <xref ref-type="table" rid="Ch1.T1"/>. It is a three-bladed, upwind, horizontal-axis wind turbine, having 24 degrees of freedom (DoFs) describing its flexibility. However, a few DoFs are enabled to obtain reduced-order linear-time-invariant (LTI) models used for controller design.</p>

<?xmltex \floatpos{t}?><table-wrap id="Ch1.T1"><?xmltex \currentcnt{1}?><label>Table 1</label><caption><p id="d1e622">The 1.5 MW WindPACT reference wind turbine specifications.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="3">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="right"/>
     <oasis:colspec colnum="3" colname="col3" align="left"/>
     <oasis:thead>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">Parameter</oasis:entry>
         <oasis:entry colname="col2">Value</oasis:entry>
         <oasis:entry colname="col3">Unit</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1">Rated rotor speed</oasis:entry>
         <oasis:entry colname="col2">20.463</oasis:entry>
         <oasis:entry colname="col3">rpm</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Hub height</oasis:entry>
         <oasis:entry colname="col2">84.288</oasis:entry>
         <oasis:entry colname="col3">m</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Cut-in, rated, cut-out wind speed</oasis:entry>
         <oasis:entry colname="col2">4, 12, 25</oasis:entry>
         <oasis:entry colname="col3">m s<inline-formula><mml:math id="M23" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Gearbox ratio</oasis:entry>
         <oasis:entry colname="col2">87.965</oasis:entry>
         <oasis:entry colname="col3">–</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Blade radius</oasis:entry>
         <oasis:entry colname="col2">35</oasis:entry>
         <oasis:entry colname="col3">m</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Rated power</oasis:entry>
         <oasis:entry colname="col2">1.5</oasis:entry>
         <oasis:entry colname="col3">MW</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Blade pitch range</oasis:entry>
         <oasis:entry colname="col2">0–90</oasis:entry>
         <oasis:entry colname="col3"><inline-formula><mml:math id="M24" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Pitch rate</oasis:entry>
         <oasis:entry colname="col2">10</oasis:entry>
         <oasis:entry colname="col3"><inline-formula><mml:math id="M25" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula> s<inline-formula><mml:math id="M26" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Optimal tip-speed ratio (<inline-formula><mml:math id="M27" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">opt</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">7.0</oasis:entry>
         <oasis:entry colname="col3">–</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Maximum power coefficient (<inline-formula><mml:math id="M28" display="inline"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mi mathvariant="normal">max</mml:mi></mml:msub></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">0.5</oasis:entry>
         <oasis:entry colname="col3">–</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Optimum pitch angle (<inline-formula><mml:math id="M29" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">β</mml:mi><mml:mi mathvariant="normal">opt</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">2.6</oasis:entry>
         <oasis:entry colname="col3"><inline-formula><mml:math id="M30" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula></oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table><?xmltex \gdef\@currentlabel{1}?></table-wrap>

      <p id="d1e868">The nonlinear generalized equation of motion for the wind turbine is expressed as
            <disp-formula id="Ch1.E3" content-type="numbered"><label>3</label><mml:math id="M31" display="block"><mml:mrow><mml:mi>M</mml:mi><mml:mo>(</mml:mo><mml:mi>q</mml:mi><mml:mo>,</mml:mo><mml:mi>u</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo><mml:mover accent="true"><mml:mi>q</mml:mi><mml:mo mathvariant="normal">¨</mml:mo></mml:mover><mml:mo>+</mml:mo><mml:mi>f</mml:mi><mml:mo>(</mml:mo><mml:mi>q</mml:mi><mml:mo>,</mml:mo><mml:mover accent="true"><mml:mi>q</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mo>,</mml:mo><mml:mi>u</mml:mi><mml:mo>,</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
          where <inline-formula><mml:math id="M32" display="inline"><mml:mi>M</mml:mi></mml:math></inline-formula> denotes the mass matrix containing inertia and mass components and <inline-formula><mml:math id="M33" display="inline"><mml:mi>f</mml:mi></mml:math></inline-formula> is a nonlinear function of the enabled DoFs <inline-formula><mml:math id="M34" display="inline"><mml:mi>q</mml:mi></mml:math></inline-formula> and their first derivative <inline-formula><mml:math id="M35" display="inline"><mml:mover accent="true"><mml:mi>q</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover></mml:math></inline-formula>, as well as the control input <inline-formula><mml:math id="M36" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula>, the disturbance input <inline-formula><mml:math id="M37" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, and time <inline-formula><mml:math id="M38" display="inline"><mml:mi>t</mml:mi></mml:math></inline-formula>. The nonlinear model Eq. (<xref ref-type="disp-formula" rid="Ch1.E3"/>) available in FAST is linearized with respect to an operating point in the above-rated region. By enabling the DoFs, which capture the most important wind turbine dynamics of interest, and specifying the operating point defined by a constant hub-height wind speed, pitch angle, and rotor speed, linearization is carried out numerically in FAST yielding periodic (azimuth-dependent) matrices of LTI models.</p>
</sec>
<sec id="Ch1.S3.SS2">
  <label>3.2</label><title>Controller for load mitigation and speed regulation</title>
      <p id="d1e1005">An adaptive robust observer-based controller, which in combination with an online damage evaluation model used for lifetime control of wind turbine components, is briefly outlined.</p>
<sec id="Ch1.S3.SS2.SSS1">
  <label>3.2.1</label><title>Robust disturbance accommodating controller</title>
      <p id="d1e1015">The RDAC, proposed in previous work <xref ref-type="bibr" rid="bib1.bibx10" id="paren.30"/>, is briefly outlined for principal understanding. To obtain an LTI model for controller design while avoiding unnecessary complexity in the linear model, six DoFs including tower fore–aft bending in first mode; drivetrain torsional displacement; blades 1, 2, and 3 flapwise displacements in first mode; variable speed generator; and drivetrain rotational flexibility are chosen. These DoFs capture the most important dynamics corresponding to the desired closed-loop performance with respect to load mitigation in wind turbine blades and the  tower, as well as generator speed regulation assuming also a flexible drivetrain. The linear model is obtained by linearizing the nonlinear model Eq. (<xref ref-type="disp-formula" rid="Ch1.E1"/>) with respect to an operating point in the above-rated wind speed region defined by a steady hub-height wind speed of <inline-formula><mml:math id="M39" display="inline"><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mi mathvariant="normal">op</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">18</mml:mn></mml:mrow></mml:math></inline-formula> m s<inline-formula><mml:math id="M40" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>, a pitch angle of <inline-formula><mml:math id="M41" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">β</mml:mi><mml:mi mathvariant="normal">op</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">20</mml:mn></mml:mrow></mml:math></inline-formula><inline-formula><mml:math id="M42" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>, and a rotor speed of <inline-formula><mml:math id="M43" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mi mathvariant="normal">op</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">20.463</mml:mn></mml:mrow></mml:math></inline-formula> rpm. The states <inline-formula><mml:math id="M44" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> used for controller design are
              <disp-formula id="Ch1.E4" content-type="numbered"><label>4</label><mml:math id="M45" display="block"><mml:mrow><mml:mi>x</mml:mi><mml:mo>=</mml:mo><mml:mfenced close="]" open="["><mml:mtable class="array" columnalign="center"><mml:mtr><mml:mtd><mml:mtext>tower-top fore–aft displacement</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>drivetrain torsional displacement</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>blade 1 flapwise displacement</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>blade 2 flapwise displacement</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>blade 3 flapwise displacement</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>generator speed</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>tower fore–aft velocity</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>drivetrain torsional velocity</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>blade 1 flapwise velocity</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>blade 2 flapwise velocity</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>blade 3 flapwise velocity</mml:mtext></mml:mtd></mml:mtr></mml:mtable></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p>
      <p id="d1e1153">The obtained reduced-order LTI model is expressed in state-space form as
              <disp-formula id="Ch1.E5" content-type="numbered"><label>5</label><mml:math id="M46" display="block"><mml:mtable class="split" rowspacing="0.2ex" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mover accent="true"><mml:mi>x</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mi>A</mml:mi><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:mi>B</mml:mi><mml:mi>u</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>B</mml:mi><mml:mi>d</mml:mi></mml:msub><mml:mi>d</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mi>y</mml:mi><mml:mo>=</mml:mo><mml:mi>C</mml:mi><mml:mi>x</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
            where <inline-formula><mml:math id="M47" display="inline"><mml:mrow><mml:mi>A</mml:mi><mml:mo>,</mml:mo><mml:mi>B</mml:mi><mml:mo>,</mml:mo><mml:msub><mml:mi>B</mml:mi><mml:mi>d</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math id="M48" display="inline"><mml:mi>C</mml:mi></mml:math></inline-formula> denote the state-space system; <inline-formula><mml:math id="M49" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> the control input, which is the collective pitch angle; <inline-formula><mml:math id="M50" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> the states; <inline-formula><mml:math id="M51" display="inline"><mml:mi>d</mml:mi></mml:math></inline-formula> the wind disturbance; and <inline-formula><mml:math id="M52" display="inline"><mml:mi>y</mml:mi></mml:math></inline-formula> the measured outputs, which include rotor speed and tower-top fore–aft bending moment.</p>
      <?pagebreak page579?><p id="d1e1261">The model Eq. (<xref ref-type="disp-formula" rid="Ch1.E5"/>) is augmented with a pitch actuator model, which accounts for the slow pitch actuator dynamics. To counteract wind disturbance effects, the model is extended with an assumed step disturbance waveform <xref ref-type="bibr" rid="bib1.bibx27 bib1.bibx28" id="paren.31"/>, which approximates sudden uniform rotor effective wind velocity fluctuations. To meet the rotor speed regulation objective with zero steady-state tracking error, the model is further extended with a partial integral action.</p>
      <p id="d1e1269">To ensure closed-loop system stability, robustness, and optimality, a mixed-sensitivity <inline-formula><mml:math id="M53" display="inline"><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mi mathvariant="normal">∞</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> norm of the closed loop transfer function is used as a cost function to optimize the disturbance accommodating controller (DAC) parameters including observer gain <inline-formula><mml:math id="M54" display="inline"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, state controller <inline-formula><mml:math id="M55" display="inline"><mml:mrow><mml:msub><mml:mi>K</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, disturbance rejection controller <inline-formula><mml:math id="M56" display="inline"><mml:mrow><mml:msub><mml:mi>K</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, and the integral gain <inline-formula><mml:math id="M57" display="inline"><mml:mrow><mml:msub><mml:mi>K</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> in a single step. The mixed sensitivity H<inline-formula><mml:math id="M58" display="inline"><mml:msub><mml:mi/><mml:mi mathvariant="normal">∞</mml:mi></mml:msub></mml:math></inline-formula> optimization problem is formulated as
              <disp-formula id="Ch1.E6" content-type="numbered"><label>6</label><mml:math id="M59" display="block"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mo>*</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:munder><mml:mi mathvariant="normal">argmin</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mo>∈</mml:mo><mml:mi mathvariant="script">R</mml:mi></mml:mrow></mml:munder><mml:msub><mml:mfenced open="∥" close="∥"><mml:mtable class="matrix" columnalign="center" framespacing="0em"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mi>S</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mi>R</mml:mi><mml:mi>S</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msub><mml:mi>T</mml:mi></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced><mml:mi mathvariant="normal">∞</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
            where <inline-formula><mml:math id="M60" display="inline"><mml:mrow><mml:msup><mml:mi>R</mml:mi><mml:mo>*</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> denotes the optimized controller and <inline-formula><mml:math id="M61" display="inline"><mml:mi mathvariant="script">R</mml:mi></mml:math></inline-formula> a set of controllers <inline-formula><mml:math id="M62" display="inline"><mml:mi>R</mml:mi></mml:math></inline-formula> that stabilize the plant. The weighting functions <inline-formula><mml:math id="M63" display="inline"><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M64" display="inline"><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math id="M65" display="inline"><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> are introduced to ensure desired robust performance, while <inline-formula><mml:math id="M66" display="inline"><mml:mi>S</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M67" display="inline"><mml:mrow><mml:mi>R</mml:mi><mml:mi>S</mml:mi></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math id="M68" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> denote the related sensitivity, control effort, and complementary sensitivity functions, respectively. The problem to find an optimal RDAC, RDAC<inline-formula><mml:math id="M69" display="inline"><mml:msup><mml:mi/><mml:mo>*</mml:mo></mml:msup></mml:math></inline-formula>, is formulated as
              <disp-formula id="Ch1.E7" content-type="numbered"><label>7</label><mml:math id="M70" display="block"><mml:mrow><mml:msup><mml:mi mathvariant="normal">RDAC</mml:mi><mml:mo>*</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:munder><mml:mi mathvariant="normal">argmin</mml:mi><mml:mrow><mml:mi mathvariant="normal">RDAC</mml:mi><mml:mo>∈</mml:mo><mml:mi mathvariant="script">R</mml:mi><mml:mi mathvariant="script">D</mml:mi><mml:mi mathvariant="script">A</mml:mi><mml:mi mathvariant="script">C</mml:mi></mml:mrow></mml:munder><mml:mo>‖</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mi mathvariant="normal">zd</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>P</mml:mi><mml:mo>,</mml:mo><mml:mi>R</mml:mi><mml:mi>D</mml:mi><mml:mi>A</mml:mi><mml:mi>C</mml:mi><mml:mo>)</mml:mo><mml:msub><mml:mo>‖</mml:mo><mml:mi mathvariant="normal">∞</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
            where <inline-formula><mml:math id="M71" display="inline"><mml:mrow><mml:mi mathvariant="script">R</mml:mi><mml:mi mathvariant="script">D</mml:mi><mml:mi mathvariant="script">A</mml:mi><mml:mi mathvariant="script">C</mml:mi></mml:mrow></mml:math></inline-formula> denotes a set of controllers RDAC that stabilize the generalized plant <inline-formula><mml:math id="M72" display="inline"><mml:mi>P</mml:mi></mml:math></inline-formula>, and <inline-formula><mml:math id="M73" display="inline"><mml:mrow><mml:msub><mml:mi>G</mml:mi><mml:mi mathvariant="normal">zd</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the transfer function from the exogenous inputs <inline-formula><mml:math id="M74" display="inline"><mml:mi>d</mml:mi></mml:math></inline-formula>  to the controlled outputs <inline-formula><mml:math id="M75" display="inline"><mml:mi>z</mml:mi></mml:math></inline-formula>.</p>
      <p id="d1e1599">Nonsmooth H<inline-formula><mml:math id="M76" display="inline"><mml:msub><mml:mi/><mml:mi mathvariant="normal">∞</mml:mi></mml:msub></mml:math></inline-formula> synthesis proposed in <xref ref-type="bibr" rid="bib1.bibx2" id="text.32"/>, used for problems with structural and stability constraints, is applied to find an optimal controller RDAC<inline-formula><mml:math id="M77" display="inline"><mml:msup><mml:mi/><mml:mo>*</mml:mo></mml:msup></mml:math></inline-formula> with robust gains <inline-formula><mml:math id="M78" display="inline"><mml:mi>L</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M79" display="inline"><mml:mi>K</mml:mi></mml:math></inline-formula> for tower load mitigation and rotor speed regulation. It is implemented in MATLAB using the <italic>hinfStruct</italic> command <xref ref-type="bibr" rid="bib1.bibx3" id="paren.33"/>. In Fig. <xref ref-type="fig" rid="Ch1.F2"/> the application of the RDAC to the 1.5 MW NREL RWT is shown. An actuator transfer function is included in the generalized plant <inline-formula><mml:math id="M80" display="inline"><mml:mi>P</mml:mi></mml:math></inline-formula>, to account for the blade pitch actuator dynamics. Hub-height wind disturbance <inline-formula><mml:math id="M81" display="inline"><mml:mi>d</mml:mi></mml:math></inline-formula> excites the wind turbine dynamics in above-rated operation. Measurement outputs including rotor speed <inline-formula><mml:math id="M82" display="inline"><mml:mi mathvariant="italic">ω</mml:mi></mml:math></inline-formula> and tower fore–aft bending moment <inline-formula><mml:math id="M83" display="inline"><mml:mi mathvariant="italic">ζ</mml:mi></mml:math></inline-formula> are fed to the RDAC, which generates a collective pitch angle <inline-formula><mml:math id="M84" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula> as a control signal for regulating rotor speed at the rated value and for reducing tower fore–aft bending moment oscillations. The RDAC is robust against modeling errors and wind disturbances.
The desired trade-off between robust stability and performance is achieved by choosing suitable weighting functions <inline-formula><mml:math id="M85" display="inline"><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">11</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M86" display="inline"><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">12</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math id="M87" display="inline"><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>. To effect rotor speed response and ensure robustness against wind disturbances, <inline-formula><mml:math id="M88" display="inline"><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">11</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is designed as an inverted low-pass filter. To reduce the first mode of tower fore–aft oscillation, <inline-formula><mml:math id="M89" display="inline"><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">12</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is designed as an inverted notch filter centered at 2.56 rad s<inline-formula><mml:math id="M90" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>. To reduce controller activity at high frequencies thereby increasing robustness, <inline-formula><mml:math id="M91" display="inline"><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is chosen as the high-pass filter.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F2"><?xmltex \currentcnt{2}?><?xmltex \def\figurename{Figure}?><label>Figure 2</label><caption><p id="d1e1763">RDAC for wind turbines.</p></caption>
            <?xmltex \igopts{width=236.157874pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f02.png"/>

          </fig>

      <p id="d1e1772">Both objectives of rotor speed regulation and tower load reduction for wind turbines operating in an above-rated wind speed region are met while ensuring robustness against modeling errors and wind disturbances. However, RDAC<inline-formula><mml:math id="M92" display="inline"><mml:msup><mml:mi/><mml:mo>*</mml:mo></mml:msup></mml:math></inline-formula> is only valid within its design operating point and suffers performance deterioration outside this envelop. Additionally, its control input signal is a collective pitch angle, hence cannot be applied for reducing blade oscillations due to vertical wind shear, which can only be achieved through IPC.</p>
</sec>
<sec id="Ch1.S3.SS2.SSS2">
  <label>3.2.2</label><title>Adaptive independent pitch controller</title>
      <p id="d1e1792">This controller is desired to counteract periodic aerodynamic loading of the rotor blades due to vertical wind shear. It is designed to reduce 1P (once per revolution; 0.333 Hz) blade flapwise oscillations and is adaptive to change in the operating point due to horizontal wind speed fluctuations. Five IPC controllers, each designed to be operational over a particular wind speed bin in the above-rated wind speed region, together with a switching mechanism based on the incoming wind speed, are used to realize aIPC. The linear models, used for designing respective IPC controllers, are extracted from the nonlinear wind turbine model Eq. (<xref ref-type="disp-formula" rid="Ch1.E1"/>) at different operating points as shown in Table <xref ref-type="table" rid="Ch1.T2"/>. It is important to note that the steady wind speeds (together with associated pitch angle and rated rotor speed) are only used to define operating points for extracting linear state-space models used for designing each of the IPC controllers. However, a stochastic wind profile is used for excitation of the closed-loop dynamic response of the wind turbine. Predefined wind speed bins are only used for thresholding based on the incoming hub-height stochastic wind speed to establish the appropriate IPC controller to be utilized in continuous operation.</p>

<?xmltex \floatpos{t}?><table-wrap id="Ch1.T2" specific-use="star"><?xmltex \currentcnt{2}?><label>Table 2</label><caption><p id="d1e1802">Design operating points for the IPC controllers.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="5">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="center"/>
     <oasis:colspec colnum="3" colname="col3" align="center"/>
     <oasis:colspec colnum="4" colname="col4" align="center"/>
     <oasis:colspec colnum="5" colname="col5" align="center"/>
     <oasis:thead>
       <oasis:row>
         <oasis:entry colname="col1">IPC</oasis:entry>
         <oasis:entry colname="col2">Wind speed bin</oasis:entry>
         <oasis:entry colname="col3">Steady wind speed</oasis:entry>
         <oasis:entry colname="col4">Blade pitch angle</oasis:entry>
         <oasis:entry colname="col5">Rotor speed</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">controller</oasis:entry>
         <oasis:entry colname="col2">[m s<inline-formula><mml:math id="M93" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>]</oasis:entry>
         <oasis:entry colname="col3">[m s<inline-formula><mml:math id="M94" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>]</oasis:entry>
         <oasis:entry colname="col4">[<inline-formula><mml:math id="M95" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>]</oasis:entry>
         <oasis:entry colname="col5">[rpm]</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1">1</oasis:entry>
         <oasis:entry colname="col2">12–15</oasis:entry>
         <oasis:entry colname="col3">14</oasis:entry>
         <oasis:entry colname="col4">13.10</oasis:entry>
         <oasis:entry colname="col5">20</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">2</oasis:entry>
         <oasis:entry colname="col2">15–17</oasis:entry>
         <oasis:entry colname="col3">16</oasis:entry>
         <oasis:entry colname="col4">16.75</oasis:entry>
         <oasis:entry colname="col5">20</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">3</oasis:entry>
         <oasis:entry colname="col2">17–19</oasis:entry>
         <oasis:entry colname="col3">18</oasis:entry>
         <oasis:entry colname="col4">19.83</oasis:entry>
         <oasis:entry colname="col5">20</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">4</oasis:entry>
         <oasis:entry colname="col2">19–21</oasis:entry>
         <oasis:entry colname="col3">20</oasis:entry>
         <oasis:entry colname="col4">22.47</oasis:entry>
         <oasis:entry colname="col5">20</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">5</oasis:entry>
         <oasis:entry colname="col2">21–25</oasis:entry>
         <oasis:entry colname="col3">22</oasis:entry>
         <oasis:entry colname="col4">24.84</oasis:entry>
         <oasis:entry colname="col5">20</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table><?xmltex \gdef\@currentlabel{2}?></table-wrap>

      <p id="d1e1986">To capture the most important dynamics  with respect to blade load mitigation and speed regulation, four DoFs including blade flapwise displacement for each blade in the first mode and variable speed generator are chosen. Correspondingly, seven states <inline-formula><mml:math id="M96" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> are included in the linear model used for designing the aIPC controller. To capture periodicity due to vertical wind shear, 24 equispaced azimuth positions are selected for linearization. To integrate this<?pagebreak page580?> periodicity in controller design, multi-blade coordinate (MBC) transformation proposed in <xref ref-type="bibr" rid="bib1.bibx8" id="text.34"/> is used to transform blade dynamics from the rotating to the non-rotating frame. The MBC-transformed reduced-order models are then averaged to obtain a weakly periodic LTI model described in state-space form as
              <disp-formula id="Ch1.E8" content-type="numbered"><label>8</label><mml:math id="M97" display="block"><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mover accent="true"><mml:mi>x</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mi>A</mml:mi><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:mi>B</mml:mi><mml:mi>u</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>B</mml:mi><mml:mi>d</mml:mi></mml:msub><mml:mi>d</mml:mi></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mi>y</mml:mi><mml:mo>=</mml:mo><mml:mi>C</mml:mi><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:mi>v</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
            where <inline-formula><mml:math id="M98" display="inline"><mml:mrow><mml:mi>A</mml:mi><mml:mo>,</mml:mo><mml:mi>B</mml:mi><mml:mo>,</mml:mo><mml:msub><mml:mi>B</mml:mi><mml:mi>d</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math id="M99" display="inline"><mml:mi>C</mml:mi></mml:math></inline-formula> denote the state-space system, <inline-formula><mml:math id="M100" display="inline"><mml:mrow><mml:mi>u</mml:mi><mml:mo>=</mml:mo><mml:mo>[</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">3</mml:mn></mml:msub><mml:msup><mml:mo>]</mml:mo><mml:mi>T</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> denotes the perturbed independent pitch angles, and <inline-formula><mml:math id="M101" display="inline"><mml:mi>d</mml:mi></mml:math></inline-formula> is the wind disturbance. The measurements <inline-formula><mml:math id="M102" display="inline"><mml:mi>y</mml:mi></mml:math></inline-formula>, which include the blade root flapwise bending moment for each blade, are assumed to be distorted with noise <inline-formula><mml:math id="M103" display="inline"><mml:mi>v</mml:mi></mml:math></inline-formula>.</p>
      <p id="d1e2140">Using a linear quadratic Gaussian (LQG) control method, Eq. (<xref ref-type="disp-formula" rid="Ch1.E8"/>) is used to design an observer-based controller. The full-state feedback controller <inline-formula><mml:math id="M104" display="inline"><mml:mi>K</mml:mi></mml:math></inline-formula> is designed using a linear quadratic regulator (LQR) technique by minimizing the cost function:
              <disp-formula id="Ch1.E9" content-type="numbered"><label>9</label><mml:math id="M105" display="block"><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mi mathvariant="normal">QR</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mi mathvariant="normal">∞</mml:mi></mml:munderover><mml:mo>(</mml:mo><mml:msup><mml:mi>x</mml:mi><mml:mi>T</mml:mi></mml:msup><mml:mi>Q</mml:mi><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:msup><mml:mi>u</mml:mi><mml:mi>T</mml:mi></mml:msup><mml:mi>R</mml:mi><mml:mi>u</mml:mi><mml:mo>)</mml:mo><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
            while solving the algebraic Riccati equation (ARE) <inline-formula><mml:math id="M106" display="inline"><mml:mrow><mml:msup><mml:mi>A</mml:mi><mml:mi>T</mml:mi></mml:msup><mml:mi>P</mml:mi><mml:mo>+</mml:mo><mml:mi>P</mml:mi><mml:mi>A</mml:mi><mml:mo>-</mml:mo><mml:mi>P</mml:mi><mml:mi>B</mml:mi><mml:msup><mml:mi>R</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup><mml:msup><mml:mi>B</mml:mi><mml:mi>T</mml:mi></mml:msup><mml:mi>P</mml:mi><mml:mo>+</mml:mo><mml:mi>Q</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula>, assuming (<inline-formula><mml:math id="M107" display="inline"><mml:mrow><mml:mi>A</mml:mi><mml:mo>,</mml:mo><mml:mi>B</mml:mi></mml:mrow></mml:math></inline-formula>) is fully controllable. Here <inline-formula><mml:math id="M108" display="inline"><mml:mi>Q</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M109" display="inline"><mml:mi>R</mml:mi></mml:math></inline-formula> denote the state and control input weighting matrices, respectively, whose elements are tuned to achieve the desired dynamic response with respect to blade load mitigation and rotor speed regulation, while <inline-formula><mml:math id="M110" display="inline"><mml:mi>P</mml:mi></mml:math></inline-formula> is the solution to the ARE. To implement optimal full-state feedback control <inline-formula><mml:math id="M111" display="inline"><mml:mrow><mml:mi>u</mml:mi><mml:mo>=</mml:mo><mml:mi>K</mml:mi><mml:mover accent="true"><mml:mi>x</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> using estimated states <inline-formula><mml:math id="M112" display="inline"><mml:mover accent="true"><mml:mi>x</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover></mml:math></inline-formula>, a Kalman state estimator is used to design the observer gain <inline-formula><mml:math id="M113" display="inline"><mml:mi>L</mml:mi></mml:math></inline-formula> by minimizing the state estimation covariance error <inline-formula><mml:math id="M114" display="inline"><mml:mrow><mml:mi>E</mml:mi><mml:mo>(</mml:mo><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>-</mml:mo><mml:mover accent="true"><mml:mi>x</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mo>)</mml:mo><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>-</mml:mo><mml:mover accent="true"><mml:mi>x</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:msup><mml:mo>)</mml:mo><mml:mi>T</mml:mi></mml:msup><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:math></inline-formula> while solving the ARE <inline-formula><mml:math id="M115" display="inline"><mml:mrow><mml:mi>A</mml:mi><mml:msub><mml:mi>P</mml:mi><mml:mi>f</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mi>f</mml:mi></mml:msub><mml:msup><mml:mi>A</mml:mi><mml:mi>T</mml:mi></mml:msup><mml:mo>-</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mi>f</mml:mi></mml:msub><mml:msup><mml:mi>C</mml:mi><mml:mi>T</mml:mi></mml:msup><mml:msubsup><mml:mi>R</mml:mi><mml:mi>f</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msubsup><mml:mi>C</mml:mi><mml:msub><mml:mi>P</mml:mi><mml:mi>f</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>Q</mml:mi><mml:mi>f</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula>, assuming (<inline-formula><mml:math id="M116" display="inline"><mml:mi>A</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M117" display="inline"><mml:mi>C</mml:mi></mml:math></inline-formula>) is fully observable. Here, <inline-formula><mml:math id="M118" display="inline"><mml:mrow><mml:msub><mml:mi>Q</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M119" display="inline"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are process disturbance and measurement noise covariance matrices, respectively, while <inline-formula><mml:math id="M120" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the solution to the ARE.</p>
      <p id="d1e2479">Figure <xref ref-type="fig" rid="Ch1.F3"/> illustrates the implementation of one of the five IPC controllers. The wind profile <inline-formula><mml:math id="M121" display="inline"><mml:mi>d</mml:mi></mml:math></inline-formula> excites the dynamics of the wind turbine in the above-rated wind speed region. The perturbed blade root flapwise bending moment measurements <inline-formula><mml:math id="M122" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mi>y</mml:mi></mml:mrow></mml:math></inline-formula> are transformed from the rotating to the fixed coordinate frame of controller design using an inverse MBC transformation matrix <inline-formula><mml:math id="M123" display="inline"><mml:mrow><mml:mi>T</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ψ</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, which relies on real-time rotor azimuth angle measurements <inline-formula><mml:math id="M124" display="inline"><mml:mi mathvariant="italic">ψ</mml:mi></mml:math></inline-formula>. The perturbed independent pitch angles <inline-formula><mml:math id="M125" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">β</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are obtained by transforming the control input <inline-formula><mml:math id="M126" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> back to the rotating frame using the MBC transformation matrix <inline-formula><mml:math id="M127" display="inline"><mml:mrow><mml:mi>T</mml:mi><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ψ</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>. By summing <inline-formula><mml:math id="M128" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">β</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and the collective pitch angle <inline-formula><mml:math id="M129" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">β</mml:mi><mml:mi>c</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> from the RDAC, the IPC signal <inline-formula><mml:math id="M130" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">β</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is obtained.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F3"><?xmltex \currentcnt{3}?><?xmltex \def\figurename{Figure}?><label>Figure 3</label><caption><p id="d1e2600">Independent pitch controller.</p></caption>
            <?xmltex \igopts{width=236.157874pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f03.png"/>

          </fig>

      <p id="d1e2609">The five IPC controllers are designed following this procedure, each at a predefined operating point, to cover the entire range of operation in the above-rated regime. A switching mechanism is implemented to activate each controller at a predefined operating range based on the incoming hub-height wind speed for system excitation using a stochastic wind profile.</p>
</sec>
</sec>
</sec>
<sec id="Ch1.S4">
  <label>4</label><?xmltex \opttitle{Control of wind turbine lifetime: an illustrative example using the 1.5\,MW NREL reference wind turbine}?><title>Control of wind turbine lifetime: an illustrative example using the 1.5 MW NREL reference wind turbine</title>
      <p id="d1e2623">To control the lifetime consumption in wind turbine blades, the adaptive robust observer-based controller (RDAC<inline-formula><mml:math id="M131" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula>aIPC), implemented using two control loops, is combined with an online damage evaluation model as shown in Fig. <xref ref-type="fig" rid="Ch1.F4"/>. A wind profile excites the wind turbine dynamics in the above-rated regime. The RDAC <xref ref-type="bibr" rid="bib1.bibx10" id="paren.35"/>, which is robust against modeling errors generates the<?pagebreak page581?> primary collective pitch controller (CPC) signal for rotor speed regulation and tower load mitigation, while aIPC is used as the lifetime controller to dynamically control the damage accumulation of the rotor blades. The IPC angles are perturbed with respect to the CPC signal from RDAC, forming the control input <inline-formula><mml:math id="M132" display="inline"><mml:mi>u</mml:mi></mml:math></inline-formula> to the wind turbine.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F4" specific-use="star"><?xmltex \currentcnt{4}?><?xmltex \def\figurename{Figure}?><label>Figure 4</label><caption><p id="d1e2647">Prognostics-based adaptive lifetime control.</p></caption>
        <?xmltex \igopts{width=426.791339pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f04.png"/>

      </fig>

      <p id="d1e2656">The blade-root flapwise bending moment measurements <inline-formula><mml:math id="M133" display="inline"><mml:mi>y</mml:mi></mml:math></inline-formula> are logged into memory during simulation. The online damage evaluation model based on the real-time implementation of the RFC algorithm <xref ref-type="bibr" rid="bib1.bibx21" id="paren.36"/> calculates the accumulated damage at every time step <inline-formula><mml:math id="M134" display="inline"><mml:mrow><mml:msub><mml:mi>D</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. The estimated lifetime of the blade <inline-formula><mml:math id="M135" display="inline"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">e</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>  used as a state-of-health (SoH) indicator is calculated as
          <disp-formula id="Ch1.E10" content-type="numbered"><label>10</label><mml:math id="M136" display="block"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">e</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>D</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:msub><mml:mi>D</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M137" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> denotes the current time step, while <inline-formula><mml:math id="M138" display="inline"><mml:mrow><mml:msub><mml:mi>D</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> denotes the accumulated damage at the design lifetime. At every time step <inline-formula><mml:math id="M139" display="inline"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, the estimated RUL can be calculated as
          <disp-formula id="Ch1.E11" content-type="numbered"><label>11</label><mml:math id="M140" display="block"><mml:mrow><mml:mi>R</mml:mi><mml:mi>U</mml:mi><mml:mi>L</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">e</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mfenced close=")" open="("><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>D</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>D</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p>
      <p id="d1e2814">Based on the threshold evaluation of <inline-formula><mml:math id="M141" display="inline"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">e</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, the load mitigation level in the respective IPC controllers is controlled by selecting the appropriate gains <inline-formula><mml:math id="M142" display="inline"><mml:mi>L</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M143" display="inline"><mml:mi>K</mml:mi></mml:math></inline-formula> every 10 ms, which is the time interval chosen for lifetime threshold evaluation. For illustrative purposes a lifetime of 600 s is chosen.  Three threshold levels are set such that if <inline-formula><mml:math id="M144" display="inline"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">e</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is below the lower limit of the desired lifetime (<inline-formula><mml:math id="M145" display="inline"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">e</mml:mi></mml:msub><mml:mspace linebreak="nobreak" width="0.125em"/><mml:mo>&lt;</mml:mo><mml:mspace width="0.125em" linebreak="nobreak"/><mml:mn mathvariant="normal">580</mml:mn></mml:mrow></mml:math></inline-formula>), maximum gains of respective IPC controllers are selected to increase the blade load mitigation level. If <inline-formula><mml:math id="M146" display="inline"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">e</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> falls within a range of the desired lifetime (<inline-formula><mml:math id="M147" display="inline"><mml:mrow><mml:mn mathvariant="normal">580</mml:mn><mml:mo>≤</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">e</mml:mi></mml:msub><mml:mo>≥</mml:mo><mml:mn mathvariant="normal">620</mml:mn></mml:mrow></mml:math></inline-formula>), optimum gains, which strike a balance between load mitigation and speed regulation, are selected. On the other hand, if the value of <inline-formula><mml:math id="M148" display="inline"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">e</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is higher than the desired lifetime (<inline-formula><mml:math id="M149" display="inline"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">e</mml:mi></mml:msub><mml:mspace linebreak="nobreak" width="0.125em"/><mml:mo>&gt;</mml:mo><mml:mspace width="0.125em" linebreak="nobreak"/><mml:mn mathvariant="normal">620</mml:mn></mml:mrow></mml:math></inline-formula>) and hence blade load mitigation level can be compromised, minimum gains are chosen.</p>
      <p id="d1e2929">It is important to note that two levels of switching are implemented as illustrated in Fig. <xref ref-type="fig" rid="Ch1.F5"/>. The first level, used for switching between different IPC controllers, is defined based on the incoming hub-height stochastic wind speed. Highly uncertain wind turbine anemometer measurements should suffice, as strict accuracy is not required for switching. Predefined wind speed bins are used for thresholding and activating a suitable IPC controller from the designed bank of controllers. This ensures that an appropriate IPC controller is used for the prevailing operating conditions. The second level of switching relies on a lifetime estimate of the blades obtained from the online damage evaluation model to adapt both the full-state feedback gain <inline-formula><mml:math id="M150" display="inline"><mml:mi>K</mml:mi></mml:math></inline-formula> and observer gain <inline-formula><mml:math id="M151" display="inline"><mml:mi>L</mml:mi></mml:math></inline-formula> of the IPC controller activated in the first level of switching to achieve the targeted lifetime. The combined switching constitutes aIPC lifetime control.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F5"><?xmltex \currentcnt{5}?><?xmltex \def\figurename{Figure}?><label>Figure 5</label><caption><p id="d1e2950">The aIPC switching implementation.</p></caption>
        <?xmltex \igopts{width=170.716535pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f05.png"/>

      </fig>

<?xmltex \hack{\newpage}?>
</sec>
<sec id="Ch1.S5">
  <label>5</label><title>Results and discussion</title>
      <p id="d1e2970">This section presents and discusses the simulation results obtained from evaluating the adaptive lifetime control strategy using the 1.5 MW NREL RWT in FAST design code. Following the IEC 61400-1 recommendation for fatigue load evaluation, a 600 s stochastic wind profile generated using TurbSim software <xref ref-type="bibr" rid="bib1.bibx16" id="paren.37"/> is used for closed-loop system excitation. The full-field IEC Kaimal type A wind profile shown in Fig. <xref ref-type="fig" rid="Ch1.F6"/>a has a mean hub-height wind speed of 18 m s<inline-formula><mml:math id="M152" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>, a turbulence intensity of 16 % at 15 m s<inline-formula><mml:math id="M153" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>, and vertical wind shear with a power-law exponent of 0.2. Although such a high wind speed has a low occurrence probability, it drives the dynamics of the wind turbine from near-rated to cut-off wind speeds. Therefore, it is useful for demonstrating the performance of the proposed control strategy over a wide range of wind turbine operations.  While blade edgewise (E–W) and tower side–side (S–S) bending moments contribute to the total fatigue damage of the respective components, in this contribution, blade flapwise (F–W) and tower fore–aft (F–A) bending moments are chosen since they are the main structural loads that drive fatigue damage of respective components in above-rated turbine operation. This sufficiently demonstrates the application of lifetime estimation of wind turbine components as a state-of-health indicator to establish a trade-off between load mitigation and speed regulation to  guarantee a given damage at a desired lifetime.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F6" specific-use="star"><?xmltex \currentcnt{6}?><?xmltex \def\figurename{Figure}?><label>Figure 6</label><caption><p id="d1e3004">Adaptive lifetime control performance.</p></caption>
        <?xmltex \igopts{width=497.923228pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f06.png"/>

      </fig>

      <p id="d1e3013">The performance of the lifetime control scheme in different blade load mitigation scenarios is shown in Fig. <xref ref-type="fig" rid="Ch1.F6"/>b. As shown, the adaptive lifetime control strategy controls the damage accumulation in the blades to reach the predefined damage limit at the  desired lifetime of 600 s. While the control strategy with maximum load mitigation achieves the same desired result, the lifetime control scheme spreads the incremental damage accumulation over the entire operation window by dynamically switching between the different load mitigation levels.</p>
      <p id="d1e3019">The baseline controller (without lifetime control) used for comparisons is RDAC combined with aIPC without lifetime control, whereby switching in aIPC is only based on incoming wind speed. To evaluate the proposed controller in varying wind fields, six profiles with mean wind speeds of 18 and 14 m s<inline-formula><mml:math id="M154" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>, each having three seeds, are used. The three wind profiles having a mean wind speed of 18 m s<inline-formula><mml:math id="M155" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> are shown in Fig. <xref ref-type="fig" rid="Ch1.F7"/>a. A comparison in blade F–W bending moment load mitigation performance is shown in Fig. <xref ref-type="fig" rid="Ch1.F7"/>b. On average (for the three wind fields), an adaptive lifetime controller achieves a 9.39 % reduction in standard deviation compared with the baseline controller. Additionally, there is significant reduction in the accumulated damage as shown in Fig. <xref ref-type="fig" rid="Ch1.F7"/>c.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F7" specific-use="star"><?xmltex \currentcnt{7}?><?xmltex \def\figurename{Figure}?><label>Figure 7</label><caption><p id="d1e3054">Blade fatigue load mitigation for 18 m s<inline-formula><mml:math id="M156" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> wind.</p></caption>
        <?xmltex \igopts{width=497.923228pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f07.png"/>

      </fig>

      <p id="d1e3075">The performance of the adaptive lifetime control strategy in mitigating tower loads is also evaluated. As illustrated in Fig. <xref ref-type="fig" rid="Ch1.F8"/>a, a significant reduction in tower F–A oscillation is observed, with the average standard deviation reducing by<?pagebreak page582?> 6.58 %. A reduction in tower damage accumulation can be seen in Fig. <xref ref-type="fig" rid="Ch1.F8"/>b. This shows that the lifetime control of blades, which reduces 1P fatigue loads, leads to reduced damage accumulation in the tower due to 3P fatigue loads. To evaluate the performance of the proposed lifetime control scheme in load mitigation in different tower and blade load channels, damage equivalent load (DEL) analysis is carried out using MLife software <xref ref-type="bibr" rid="bib1.bibx14" id="paren.38"/>. Based on fatigue analysis carried out with results obtained using the 10 min stochastic wind profile shown in Fig. <xref ref-type="fig" rid="Ch1.F7"/>a, the lifetime controller reduces DELs in the blade F–W and tower F–A as illustrated in Fig. <xref ref-type="fig" rid="Ch1.F8"/>c. No noticeable change in blade E–W DEL is achieved. However, a slight increase is seen in tower S–S DEL, which is attributed to a slight increase in pitch activity for improved load reduction.</p>

      <?xmltex \floatpos{p}?><fig id="Ch1.F8" specific-use="star"><?xmltex \currentcnt{8}?><?xmltex \def\figurename{Figure}?><label>Figure 8</label><caption><p id="d1e3091">Tower fatigue load mitigation and DEL analysis for 18 m s<inline-formula><mml:math id="M157" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> wind.</p></caption>
        <?xmltex \igopts{width=497.923228pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f08.png"/>

      </fig>

      <p id="d1e3112">Despite the adaptive lifetime controller achieving improved performance in reducing damage accumulation in both rotor blade and tower, it is important to ascertain that this does not compromise the speed and power regulation performance. To illustrate this, the rotor speed and generator power are evaluated as shown in Fig. <xref ref-type="fig" rid="Ch1.F9"/>. With lifetime control, no notable change is realized with respect to rotor speed regulation. Although a slight increase of 8.7 % in power standard deviation is noted, the generated power fluctuates within acceptable limits. The mean power is identical at 1560.13   and 1559.97 kW for the baseline and lifetime controllers, respectively. Improvement in both load mitigation and rotor speed regulation is achieved with insignificant additional pitch activity as illustrated in Fig. <xref ref-type="fig" rid="Ch1.F9"/>c. The average total pitch travel marginally increases by 0.13 %.</p>

      <?xmltex \floatpos{p}?><fig id="Ch1.F9" specific-use="star"><?xmltex \currentcnt{9}?><?xmltex \def\figurename{Figure}?><label>Figure 9</label><caption><p id="d1e3122">Speed and power regulation performance and pitch actuator usage for 18 m s<inline-formula><mml:math id="M158" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> wind.</p></caption>
        <?xmltex \igopts{width=497.923228pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f09.png"/>

      </fig>

      <p id="d1e3143">Given that an 18 m s<inline-formula><mml:math id="M159" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> wind field realization has a low occurrence probability since wind turbines spend most of the time operating in near-rated wind conditions, the proposed adaptive lifetime control strategy is evaluated using a near-rated wind profile. For this, three IEC type C stochastic wind field realizations shown in Fig. <xref ref-type="fig" rid="Ch1.F10"/>a, each with a mean speed of 14 m s<inline-formula><mml:math id="M160" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> and a turbulence intensity of 12 % at 15 m s<inline-formula><mml:math id="M161" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>, are used. Fatigue load mitigation performance of the proposed lifetime controller in the blades is evaluated against the baseline controller as illustrated in Fig. <xref ref-type="fig" rid="Ch1.F10"/>b. The lifetime controller achieves a 10.1 % reduction in standard deviation in  blade F–W bending moment. Additionally, significant reduction in the accumulated damage is achieved.</p>

      <?xmltex \floatpos{p}?><fig id="Ch1.F10" specific-use="star"><?xmltex \currentcnt{10}?><?xmltex \def\figurename{Figure}?><label>Figure 10</label><caption><p id="d1e3188">Blade fatigue load mitigation for 14 m s<inline-formula><mml:math id="M162" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> wind.</p></caption>
        <?xmltex \igopts{width=497.923228pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f10.png"/>

      </fig>

      <p id="d1e3209">A reduction in tower fore–aft loading and damage accumulation is realized as shown in Fig. <xref ref-type="fig" rid="Ch1.F11"/>a. The standard deviation in tower loading is reduced by 11.2 %.  Fatigue analysis is carried out using simulation results based on the wind profiles shown in Fig. <xref ref-type="fig" rid="Ch1.F10"/>a. As illustrated in Fig. <xref ref-type="fig" rid="Ch1.F11"/>c, the adaptive lifetime controller achieves DEL reduction in all load channels except the slight increase in tower S–S DEL.</p>

      <?xmltex \floatpos{p}?><fig id="Ch1.F11" specific-use="star"><?xmltex \currentcnt{11}?><?xmltex \def\figurename{Figure}?><label>Figure 11</label><caption><p id="d1e3220">Tower fatigue load mitigation and DEL analysis for 14 m s<inline-formula><mml:math id="M163" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> wind.</p></caption>
        <?xmltex \igopts{width=497.923228pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f11.png"/>

      </fig>

      <p id="d1e3242">Speed and power regulation performance of the proposed controller is also evaluated at near-rated wind conditions as shown in Fig. <xref ref-type="fig" rid="Ch1.F12"/>. With lifetime control, improvement is realized in both speed and power regulation, with the standard deviation in rotor speed and generated power reducing<?pagebreak page583?> by 1.3 % and 1.2 %, respectively. The mean power is identical at 1553.5  and 1553.73 kW for the baseline and lifetime controllers, respectively. Improvement in both load mitigation and speed regulation comes with insignificant additional pitch activity as illustrated in Fig. <xref ref-type="fig" rid="Ch1.F12"/>c. The average total pitch travel increases slightly by 0.5 %. Therefore, the proposed adaptive lifetime control strategy performs well in near-rated wind conditions.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F12" specific-use="star"><?xmltex \currentcnt{12}?><?xmltex \def\figurename{Figure}?><label>Figure 12</label><caption><p id="d1e3251">Speed and power regulation performance and pitch actuator usage for 14 m s<inline-formula><mml:math id="M164" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> wind.</p></caption>
        <?xmltex \igopts{width=497.923228pt}?><graphic xlink:href="https://wes.copernicus.org/articles/8/575/2023/wes-8-575-2023-f12.png"/>

      </fig>

</sec>
<sec id="Ch1.S6" sec-type="conclusions">
  <label>6</label><title>Summary and conclusion</title>
      <p id="d1e3280">In this paper, a prognostics-based adaptive control strategy for lifetime control of wind turbines is presented. A robust disturbance accommodating controller (RDAC) designed using mixed sensitivity <inline-formula><mml:math id="M165" display="inline"><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mi mathvariant="normal">∞</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> control is used as the primary controller for mitigating tower loads and regulating rotor speed using a CPC signal. On the other hand, an aIPC controller designed using the LQG control method is used as a lifetime controller. The gains of each of its five IPC controllers are adapted based on the state of health of the rotor blades obtained using an online damage evaluation model to strike a compromise between lifetime control through load mitigation and speed regulation.</p>
      <p id="d1e3294">Through simulation using a 1.5 MW<?pagebreak page586?> wind turbine model, it is demonstrated that the adaptive lifetime control strategy controls the damage accumulation in the blades to guarantee a given damage limit at the desired lifetime. A reduction in accumulated damage in the tower is also realized. This can potentially be used for optimizing maintenance scheduling in wind farms by synchronizing aging of wind turbine components, hence reducing operation and maintenance costs and increasing operational reliability. Fatigue analysis indicates a reduction in DELs in most load channels. This improvement is realized without compromise in the speed and power regulation performance. The lifetime controller achieves these results without significant increase in pitch actuator duty cycle. In the future, adaptive lifetime control based on nonlinear damage accumulation models will be considered. Additionally, the use of new concepts for state-of-health indicators such as a change in modal parameters for structural health monitoring will be explored.</p>
</sec>

      
      </body>
    <back><notes notes-type="codeavailability"><title>Code availability</title>

      <p id="d1e3301">The code is not publicly available and cannot be shared.</p>
  </notes><notes notes-type="dataavailability"><title>Data availability</title>

      <p id="d1e3307">Simulation data were obtained using Fatigue Aerodynamics Structures and Turbulence  version 7 (FAST v7), a free wind turbine simulation tool developed by National Renewable Energy Laboratory (NREL). For additional information see the FAST user guide  (<uri>https://www.nrel.gov/wind/nwtc/fastv7.html</uri>, <xref ref-type="bibr" rid="bib1.bibx22" id="altparen.39"/>).</p>
  </notes><notes notes-type="authorcontribution"><title>Author contributions</title>

      <p id="d1e3319">DS and JGN proposed the original idea of combining an online damage evaluation model with a switching IPC load  mitigation control scheme for the lifetime extension of wind turbine blades. MHD and DS extended the work of JGN by developing an adaptive lifetime controller (RDAC) to control life consumption in the tower. Based on these ideas and with supervision from DS, EK developed an adaptive lifetime controller (aIPC), which is adaptive to wind speed changes and in which the gains of each IPC controller are adapted to dynamically control the damage accumulation of rotor blades while monitoring the life consumption of the tower. EK evaluated the prognostics scheme (RDAC<inline-formula><mml:math id="M166" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula>aIPC) on a 1.5 MW RWT by running simulations and analyzing obtained results. With valuable input from DS, EK wrote the manuscript. All authors contributed to this work from concept to manuscript stage.</p>
  </notes><notes notes-type="competinginterests"><title>Competing interests</title>

      <p id="d1e3332">The contact author has declared that none of the authors has any competing interests.</p>
  </notes><notes notes-type="disclaimer"><title>Disclaimer</title>

      <p id="d1e3338">Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p>
  </notes><ack><title>Acknowledgements</title><p id="d1e3344">This work is partly supported through a scholarship awarded to Edwin Kipchirchir by the German Academic Exchange Service (DAAD) in cooperation with the Ministry of Education of Kenya for his PhD study at the Chair of Dynamics and Control, UDE, Germany. We also thank Paul Fleming and the three anonymous reviewers who helped to improve the manuscript.</p></ack><notes notes-type="reviewstatement"><title>Review statement</title>

      <p id="d1e3349">This paper was edited by Paul Fleming and reviewed by three anonymous referees.</p>
  </notes><ref-list>
    <title>References</title>

      <ref id="bib1.bibx1"><?xmltex \def\ref@label{{Acho et~al.(2016)}}?><label>Acho et al.(2016)</label><?label Acho2016?><mixed-citation>Acho, L., Rodellar, J., Tutivén, C., and Vidal, Y.: Passive fault tolerant
control strategy in controlled wind turbines, in: 2016 3rd Conference on
Control and Fault-Tolerant Systems (SysTol),   IEEE, Barcelona,
Spain, 7–9 September 2016, 636–641, <ext-link xlink:href="https://doi.org/10.1109/SYSTOL.2016.7739820" ext-link-type="DOI">10.1109/SYSTOL.2016.7739820</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bibx2"><?xmltex \def\ref@label{{Apkarian and Noll(2006)}}?><label>Apkarian and Noll(2006)</label><?label Apkarian2006?><mixed-citation>Apkarian, P. and Noll, D.: Nonsmooth H<inline-formula><mml:math id="M167" display="inline"><mml:msub><mml:mi/><mml:mi mathvariant="normal">∞</mml:mi></mml:msub></mml:math></inline-formula> synthesis, IEEE T. Automat. Contr., 51, 71–86, <ext-link xlink:href="https://doi.org/10.1109/TAC.2005.860290" ext-link-type="DOI">10.1109/TAC.2005.860290</ext-link>, 2006.</mixed-citation></ref>
      <ref id="bib1.bibx3"><?xmltex \def\ref@label{{Apkarian and Noll(2017)}}?><label>Apkarian and Noll(2017)</label><?label Apkarian2017?><mixed-citation>Apkarian, P. and Noll, D.: The H<inline-formula><mml:math id="M168" display="inline"><mml:msub><mml:mi/><mml:mi mathvariant="normal">∞</mml:mi></mml:msub></mml:math></inline-formula> control problem is solved,
Aerospace Lab, 13,  1–11, <ext-link xlink:href="https://doi.org/10.12762/2017.AL13-01" ext-link-type="DOI">10.12762/2017.AL13-01</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bibx4"><?xmltex \def\ref@label{{Azizi et~al.(2019)Azizi, Nourisola, and Shoja-Majidabad}}?><label>Azizi et al.(2019)Azizi, Nourisola, and Shoja-Majidabad</label><?label Azizi2019?><mixed-citation>Azizi, A., Nourisola, H., and Shoja-Majidabad, S.: Fault tolerant control of
wind turbines with an adaptive output feedback sliding mode controller,
Renew. Energ., 135, 55–65, <ext-link xlink:href="https://doi.org/10.1016/j.renene.2018.11.106" ext-link-type="DOI">10.1016/j.renene.2018.11.106</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bibx5"><?xmltex \def\ref@label{{Bajri{\'{c}} et~al.(2018)Bajri{\'{c}}, Jan, and
R{\"{u}}dinger}}?><label>Bajrić et al.(2018)Bajrić, Jan, and
Rüdinger</label><?label Bajric2017?><mixed-citation>Bajrić, A., Jan, H., and Rüdinger, F.: Evaluation of damping estimates
by automated operational modal analysis for offshore wind turbine tower
vibrations, Renew. Energ., 116, 153–163,
<ext-link xlink:href="https://doi.org/10.1016/j.renene.2017.03.043" ext-link-type="DOI">10.1016/j.renene.2017.03.043</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bibx6"><?xmltex \def\ref@label{{Beganovic et~al.(2015)Beganovic, Njiri, Rothe, and
S{\"{o}}ffker}}?><label>Beganovic et al.(2015)Beganovic, Njiri, Rothe, and
Söffker</label><?label Beganovic2015?><mixed-citation>Beganovic, N., Njiri, J. G., Rothe, S., and Söffker, D.: Application of
Diagnosis and Prognosis to Wind Turbine System Based on Fatigue Load, in:
2015 IEEE Conference on Prognostics and Health Management (PHM),
IEEE, Austin, Texas, 22–25 June 2015, 1–6, <ext-link xlink:href="https://doi.org/10.1109/ICPHM.2015.7245031" ext-link-type="DOI">10.1109/ICPHM.2015.7245031</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bibx7"><?xmltex \def\ref@label{{Beganovic et~al.(2018)Beganovic, Njiri, and
S{\"{o}}ffker}}?><label>Beganovic et al.(2018)Beganovic, Njiri, and
Söffker</label><?label Beganovic2018?><mixed-citation>Beganovic, N., Njiri, J. G., and Söffker, D.: Reduction of structural
loads in wind turbines based on adapted control strategy concerning online
fatigue damage evaluation models,  Energies,  11,
3429, <ext-link xlink:href="https://doi.org/10.3390/en11123429" ext-link-type="DOI">10.3390/en11123429</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bibx8"><?xmltex \def\ref@label{{Bir(2010)}}?><label>Bir(2010)</label><?label Bir2010?><mixed-citation>Bir, G. S.: User's guide to MBC3: Multi-blade coordinate transformation code
for 3-bladed wind turbine, Tech. rep., National Renewable Energy Laboratory
(NREL),  <uri>https://www.nrel.gov/docs/fy10osti/44327.pdf</uri> (last access: 19 April 2023), 2010.</mixed-citation></ref>
      <ref id="bib1.bibx9"><?xmltex \def\ref@label{{Do and S{\"{o}}ffker(2020)}}?><label>Do and Söffker(2020)</label><?label Do2020a?><mixed-citation>Do, M. H. and Söffker, D.: Wind turbine lifetime control using structural
health monitoring and prognosis, IFAC-PapersOnLine, 53, 12669–12674,
<ext-link xlink:href="https://doi.org/10.1016/j.ifacol.2020.12.1847" ext-link-type="DOI">10.1016/j.ifacol.2020.12.1847</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bibx10"><?xmltex \def\ref@label{{Do and S{\"{o}}ffker(2022)}}?><label>Do and Söffker(2022)</label><?label Do2021?><mixed-citation>Do, M. H. and Söffker, D.: Wind turbine robust disturbance accommodating
control using non-smooth H<inline-formula><mml:math id="M169" display="inline"><mml:msub><mml:mi/><mml:mi mathvariant="normal">∞</mml:mi></mml:msub></mml:math></inline-formula> optimization, Wind Energy, 25,
107–124, <ext-link xlink:href="https://doi.org/10.1002/we.2663" ext-link-type="DOI">10.1002/we.2663</ext-link>,
2022.</mixed-citation></ref>
      <ref id="bib1.bibx11"><?xmltex \def\ref@label{{Dong et~al.(2018)Dong, Lian, Wang, Yu, and Zhao}}?><label>Dong et al.(2018)Dong, Lian, Wang, Yu, and Zhao</label><?label Dong2017?><mixed-citation>Dong, X., Lian, J., Wang, H., Yu, T., and Zhao, Y.: Structural vibration
monitoring and operational modal analysis of offshore wind turbine
structure, Ocean Eng., 150, 280–297,
<ext-link xlink:href="https://doi.org/10.1016/j.oceaneng.2017.12.052" ext-link-type="DOI">10.1016/j.oceaneng.2017.12.052</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bibx12"><?xmltex \def\ref@label{{{El Maati} and {El Bahir}(2020)}}?><label>El Maati and El Bahir(2020)</label><?label Elmaati2020?><mixed-citation>El Maati, Y. A. and El Bahir, L.: Optimal fault tolerant control of
large-scale wind turbines in the case of the pitch actuator partial faults,
Complexity, 2020, 1–17, <ext-link xlink:href="https://doi.org/10.1155/2020/6210407" ext-link-type="DOI">10.1155/2020/6210407</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bibx13"><?xmltex \def\ref@label{{Gao and Liu(2021)}}?><label>Gao and Liu(2021)</label><?label Gao2021?><mixed-citation>Gao, Z. and Liu, X.: An overview on fault diagnosis, prognosis and resilient
control for wind turbine systems, Processes, 9, 300,
<ext-link xlink:href="https://doi.org/10.3390/pr9020300" ext-link-type="DOI">10.3390/pr9020300</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bibx14"><?xmltex \def\ref@label{{Hayman(2012)}}?><label>Hayman(2012)</label><?label Hayman2012?><mixed-citation>Hayman, G.: MLife theory manual for version 1.00, Tech. rep., National
Renewable Energy Laboratory (NREL), Golden, Colorado,
<uri>https://www.nrel.gov/wind/nwtc/assets/pdfs/mlife-theory.pdf</uri> (last access: 19 April 2023),
2012.</mixed-citation></ref>
      <ref id="bib1.bibx15"><?xmltex \def\ref@label{{Jain and Yam{\'{e}}(2020)}}?><label>Jain and Yamé(2020)</label><?label Jain2020?><mixed-citation>Jain, T. and Yamé, J.: Health-aware fault-tolerant receding horizon
control of wind turbines, Control Eng. Pract., 95, 104236,
<ext-link xlink:href="https://doi.org/10.1016/j.conengprac.2019.104236" ext-link-type="DOI">10.1016/j.conengprac.2019.104236</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bibx16"><?xmltex \def\ref@label{{Jonkman and Kilcher(2012)}}?><label>Jonkman and Kilcher(2012)</label><?label Jonkman2012?><mixed-citation>Jonkman, B. J. and Kilcher, L.: TurbSim user's guide: version 1.06.00, Tech.
rep., National Renewable Energy Laboratory (NREL), Golden, Colorado, USA,
<uri>https://www.nrel.gov/wind/nwtc/assets/pdfs/turbsim.pdf</uri> (last access: 19 April 2023), 2012.</mixed-citation></ref>
      <ref id="bib1.bibx17"><?xmltex \def\ref@label{{Jonkman and {Buhl}(2005)}}?><label>Jonkman and Buhl(2005)</label><?label Jonkman2005?><mixed-citation>Jonkman, J. M. and Buhl Jr., M. L.: FAST user's guide, Tech. Rep.
NREL/EL-500-29798, National Renewable Energy Laboratory (NREL), Golden,
Colorado, USA, <uri>https://www.nrel.gov/docs/fy06osti/38230.pdf</uri> (last access: 19 April 2023),
2005.</mixed-citation></ref>
      <ref id="bib1.bibx18"><?xmltex \def\ref@label{{Kim et~al.(2019)Kim, Kim, and Choe}}?><label>Kim et al.(2019)Kim, Kim, and Choe</label><?label Kim2019?><mixed-citation>Kim, H.-C., Kim, M.-H., and Choe, D.-E.: Structural health monitoring of
towers and blades for floating offshore wind turbines using operational modal
analysis and modal properties with numerical-sensor signals, Ocean
Eng., 188, 106226, <ext-link xlink:href="https://doi.org/10.1016/j.oceaneng.2019.106226" ext-link-type="DOI">10.1016/j.oceaneng.2019.106226</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bibx19"><?xmltex \def\ref@label{{Matsuishi and Endo(1968)}}?><label>Matsuishi and Endo(1968)</label><?label Matsuishi1968?><mixed-citation>
Matsuishi, M. and Endo, T.: Fatigue of metals subjected to varying stress,
Japan Society of Mechanical Engineers, 68, 37–40, 1968.</mixed-citation></ref>
      <ref id="bib1.bibx20"><?xmltex \def\ref@label{{Miner(1945)}}?><label>Miner(1945)</label><?label Miner1945?><mixed-citation>Miner, M. A.: Cumulative damage in fatigue, J. Appl. Mech.,
12,  A159–A164, <ext-link xlink:href="https://doi.org/10.1115/1.4009458" ext-link-type="DOI">10.1115/1.4009458</ext-link>, 1945.</mixed-citation></ref>
      <ref id="bib1.bibx21"><?xmltex \def\ref@label{{Musallam and Johnson(2012)}}?><label>Musallam and Johnson(2012)</label><?label Musallam2012?><mixed-citation>Musallam, M. and Johnson, C. M.: An efficient implementation of the rainflow
counting algorithm for life consumption estimation, IEEE T.
Reliab., 61, 978–986, <ext-link xlink:href="https://doi.org/10.1109/TR.2012.2221040" ext-link-type="DOI">10.1109/TR.2012.2221040</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bibx22"><?xmltex \def\ref@label{National Renewable Energy Laboratory(2023)}?><label>National Renewable Energy Laboratory(2023)</label><?label NREL2023?><mixed-citation>National Renewable Energy Laboratory (NREL):   FAST v7, <uri>https://www.nrel.gov/wind/nwtc/fastv7.html</uri>, last access: 19 April 2023.</mixed-citation></ref>
      <ref id="bib1.bibx23"><?xmltex \def\ref@label{{Njiri et~al.(2019)Njiri, Beganovic, Do, and S{\"{o}}ffker}}?><label>Njiri et al.(2019)Njiri, Beganovic, Do, and Söffker</label><?label Njiri2019?><mixed-citation>Njiri, J. G., Beganovic, N., Do, M. H., and Söffker, D.: Consideration of
lifetime and fatigue load in wind turbine control, Renew. Energ., 131,
818–828, <ext-link xlink:href="https://doi.org/10.1016/j.renene.2018.07.109" ext-link-type="DOI">10.1016/j.renene.2018.07.109</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bibx24"><?xmltex \def\ref@label{{Pegalajar-Jurado and Bredmose(2019)}}?><label>Pegalajar-Jurado and Bredmose(2019)</label><?label Pegalajar2019?><mixed-citation>Pegalajar-Jurado, A. and Bredmose, H.: Reproduction of slow-drift motions of a
floating wind turbine using second-order hydrodynamics and operational modal
analysis, Mar. Struct., 66, 178–196,
<ext-link xlink:href="https://doi.org/10.1016/j.marstruc.2019.02.008" ext-link-type="DOI">10.1016/j.marstruc.2019.02.008</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bibx25"><?xmltex \def\ref@label{{Ragan and Manuel(2007)}}?><label>Ragan and Manuel(2007)</label><?label Ragan2007?><mixed-citation>Ragan, P. and Manuel, L.: Comparing estimates of wind turbine fatigue loads
using time-domain and spectral methods, Wind Eng., 31, 83–99,
<ext-link xlink:href="https://doi.org/10.1260/030952407781494494" ext-link-type="DOI">10.1260/030952407781494494</ext-link>, 2007.</mixed-citation></ref>
      <ref id="bib1.bibx26"><?xmltex \def\ref@label{{Rinker and Dykes(2018)}}?><label>Rinker and Dykes(2018)</label><?label Rinker2018?><mixed-citation>Rinker, J. and Dykes, K.: WindPACT reference wind turbines, Tech. rep.,
National Renewable Energy Laboratory (NREL), Golden, Colorado, USA,
<uri>https://www.nrel.gov/docs/fy18osti/67667.pdf</uri> (last access: 19 April 2023), 2018.</mixed-citation></ref>
      <ref id="bib1.bibx27"><?xmltex \def\ref@label{{Wright(2003)}}?><label>Wright(2003)</label><?label Wright2004?><mixed-citation>Wright, A. D.: Modern control design for flexible wind turbines, University
of Colorado at Boulder, Golden, Colorado, USA,
<uri>https://www.proquest.com/openview/f8f93112e807853ee3cbc16d4c8a64e7/1?pq-origsite=gscholar&amp;cbl=18750&amp;diss=y</uri> (last access: 19 April 2023),
2003.</mixed-citation></ref>
      <ref id="bib1.bibx28"><?xmltex \def\ref@label{{Wright and Fingersh(2008)}}?><label>Wright and Fingersh(2008)</label><?label Wright2008?><mixed-citation>Wright, A. D. and Fingersh, L. J.: Advanced control design for wind turbines
part I: Control design, implementation, and initial tests, Tech. Rep.
NREL/TP-500-42437, National Renewable Energy Laboratory (NREL), Golden,
Colorado, USA,  <uri>https://www.osti.gov/biblio/927269</uri> (last access: 19 April 2023), 2008.</mixed-citation></ref>
      <ref id="bib1.bibx29"><?xmltex \def\ref@label{{Yuan et~al.(2015)Yuan, Li, Huang, Zhu, and Gao}}?><label>Yuan et al.(2015)Yuan, Li, Huang, Zhu, and Gao</label><?label Yuan2014?><mixed-citation>Yuan, R., Li, H., Huang, H. Z., Zhu, S. P., and Gao, H.: A nonlinear fatigue
damage accumulation model considering strength degradation and its
applications to fatigue reliability analysis, Int. J. Damage Mech.,
24,  646–662,
<ext-link xlink:href="https://doi.org/10.1177/1056789514544228" ext-link-type="DOI">10.1177/1056789514544228</ext-link>, 2015.</mixed-citation></ref>
      <?pagebreak page588?><ref id="bib1.bibx30"><?xmltex \def\ref@label{{Ziegler et~al.(2018)Ziegler, Gonzalez, Rubert, Smolka, and
Melero}}?><label>Ziegler et al.(2018)Ziegler, Gonzalez, Rubert, Smolka, and
Melero</label><?label Ziegler2018?><mixed-citation>Ziegler, L., Gonzalez, E., Rubert, T., Smolka, U., and Melero, J. J.: Lifetime
extension of onshore wind turbines: A review covering Germany, Spain,
Denmark, and the UK, Renewable and Sustainable Energy Reviews, 82,
1261–1271, <ext-link xlink:href="https://doi.org/10.1016/j.rser.2017.09.100" ext-link-type="DOI">10.1016/j.rser.2017.09.100</ext-link>, 2018.</mixed-citation></ref>

  </ref-list></back>
    <!--<article-title-html>Prognostics-based adaptive control strategy  for lifetime control of wind turbines</article-title-html>
<abstract-html/>
<ref-html id="bib1.bib1"><label>Acho et al.(2016)</label><mixed-citation>
      
Acho, L., Rodellar, J., Tutivén, C., and Vidal, Y.: Passive fault tolerant
control strategy in controlled wind turbines, in: 2016 3rd Conference on
Control and Fault-Tolerant Systems (SysTol),   IEEE, Barcelona,
Spain, 7–9 September 2016, 636–641, <a href="https://doi.org/10.1109/SYSTOL.2016.7739820" target="_blank">https://doi.org/10.1109/SYSTOL.2016.7739820</a>, 2016.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>Apkarian and Noll(2006)</label><mixed-citation>
      
Apkarian, P. and Noll, D.: Nonsmooth H<sub>∞</sub> synthesis, IEEE T. Automat. Contr., 51, 71–86, <a href="https://doi.org/10.1109/TAC.2005.860290" target="_blank">https://doi.org/10.1109/TAC.2005.860290</a>, 2006.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>Apkarian and Noll(2017)</label><mixed-citation>
      
Apkarian, P. and Noll, D.: The H<sub>∞</sub> control problem is solved,
Aerospace Lab, 13,  1–11, <a href="https://doi.org/10.12762/2017.AL13-01" target="_blank">https://doi.org/10.12762/2017.AL13-01</a>, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>Azizi et al.(2019)Azizi, Nourisola, and Shoja-Majidabad</label><mixed-citation>
      
Azizi, A., Nourisola, H., and Shoja-Majidabad, S.: Fault tolerant control of
wind turbines with an adaptive output feedback sliding mode controller,
Renew. Energ., 135, 55–65, <a href="https://doi.org/10.1016/j.renene.2018.11.106" target="_blank">https://doi.org/10.1016/j.renene.2018.11.106</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>Bajrić et al.(2018)Bajrić, Jan, and
Rüdinger</label><mixed-citation>
      
Bajrić, A., Jan, H., and Rüdinger, F.: Evaluation of damping estimates
by automated operational modal analysis for offshore wind turbine tower
vibrations, Renew. Energ., 116, 153–163,
<a href="https://doi.org/10.1016/j.renene.2017.03.043" target="_blank">https://doi.org/10.1016/j.renene.2017.03.043</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>Beganovic et al.(2015)Beganovic, Njiri, Rothe, and
Söffker</label><mixed-citation>
      
Beganovic, N., Njiri, J. G., Rothe, S., and Söffker, D.: Application of
Diagnosis and Prognosis to Wind Turbine System Based on Fatigue Load, in:
2015 IEEE Conference on Prognostics and Health Management (PHM),
IEEE, Austin, Texas, 22–25 June 2015, 1–6, <a href="https://doi.org/10.1109/ICPHM.2015.7245031" target="_blank">https://doi.org/10.1109/ICPHM.2015.7245031</a>, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>Beganovic et al.(2018)Beganovic, Njiri, and
Söffker</label><mixed-citation>
      
Beganovic, N., Njiri, J. G., and Söffker, D.: Reduction of structural
loads in wind turbines based on adapted control strategy concerning online
fatigue damage evaluation models,  Energies,  11,
3429, <a href="https://doi.org/10.3390/en11123429" target="_blank">https://doi.org/10.3390/en11123429</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>Bir(2010)</label><mixed-citation>
      
Bir, G. S.: User's guide to MBC3: Multi-blade coordinate transformation code
for 3-bladed wind turbine, Tech. rep., National Renewable Energy Laboratory
(NREL),  <a href="https://www.nrel.gov/docs/fy10osti/44327.pdf" target="_blank"/> (last access: 19 April 2023), 2010.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>Do and Söffker(2020)</label><mixed-citation>
      
Do, M. H. and Söffker, D.: Wind turbine lifetime control using structural
health monitoring and prognosis, IFAC-PapersOnLine, 53, 12669–12674,
<a href="https://doi.org/10.1016/j.ifacol.2020.12.1847" target="_blank">https://doi.org/10.1016/j.ifacol.2020.12.1847</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>Do and Söffker(2022)</label><mixed-citation>
      
Do, M. H. and Söffker, D.: Wind turbine robust disturbance accommodating
control using non-smooth H<sub>∞</sub> optimization, Wind Energy, 25,
107–124, <a href="https://doi.org/10.1002/we.2663" target="_blank">https://doi.org/10.1002/we.2663</a>,
2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>Dong et al.(2018)Dong, Lian, Wang, Yu, and Zhao</label><mixed-citation>
      
Dong, X., Lian, J., Wang, H., Yu, T., and Zhao, Y.: Structural vibration
monitoring and operational modal analysis of offshore wind turbine
structure, Ocean Eng., 150, 280–297,
<a href="https://doi.org/10.1016/j.oceaneng.2017.12.052" target="_blank">https://doi.org/10.1016/j.oceaneng.2017.12.052</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>El Maati and El Bahir(2020)</label><mixed-citation>
      
El Maati, Y. A. and El Bahir, L.: Optimal fault tolerant control of
large-scale wind turbines in the case of the pitch actuator partial faults,
Complexity, 2020, 1–17, <a href="https://doi.org/10.1155/2020/6210407" target="_blank">https://doi.org/10.1155/2020/6210407</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>Gao and Liu(2021)</label><mixed-citation>
      
Gao, Z. and Liu, X.: An overview on fault diagnosis, prognosis and resilient
control for wind turbine systems, Processes, 9, 300,
<a href="https://doi.org/10.3390/pr9020300" target="_blank">https://doi.org/10.3390/pr9020300</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>Hayman(2012)</label><mixed-citation>
      
Hayman, G.: MLife theory manual for version 1.00, Tech. rep., National
Renewable Energy Laboratory (NREL), Golden, Colorado,
<a href="https://www.nrel.gov/wind/nwtc/assets/pdfs/mlife-theory.pdf" target="_blank"/> (last access: 19 April 2023),
2012.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>Jain and Yamé(2020)</label><mixed-citation>
      
Jain, T. and Yamé, J.: Health-aware fault-tolerant receding horizon
control of wind turbines, Control Eng. Pract., 95, 104236,
<a href="https://doi.org/10.1016/j.conengprac.2019.104236" target="_blank">https://doi.org/10.1016/j.conengprac.2019.104236</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>Jonkman and Kilcher(2012)</label><mixed-citation>
      
Jonkman, B. J. and Kilcher, L.: TurbSim user's guide: version 1.06.00, Tech.
rep., National Renewable Energy Laboratory (NREL), Golden, Colorado, USA,
<a href="https://www.nrel.gov/wind/nwtc/assets/pdfs/turbsim.pdf" target="_blank"/> (last access: 19 April 2023), 2012.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>Jonkman and Buhl(2005)</label><mixed-citation>
      
Jonkman, J. M. and Buhl Jr., M. L.: FAST user's guide, Tech. Rep.
NREL/EL-500-29798, National Renewable Energy Laboratory (NREL), Golden,
Colorado, USA, <a href="https://www.nrel.gov/docs/fy06osti/38230.pdf" target="_blank"/> (last access: 19 April 2023),
2005.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>Kim et al.(2019)Kim, Kim, and Choe</label><mixed-citation>
      
Kim, H.-C., Kim, M.-H., and Choe, D.-E.: Structural health monitoring of
towers and blades for floating offshore wind turbines using operational modal
analysis and modal properties with numerical-sensor signals, Ocean
Eng., 188, 106226, <a href="https://doi.org/10.1016/j.oceaneng.2019.106226" target="_blank">https://doi.org/10.1016/j.oceaneng.2019.106226</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>Matsuishi and Endo(1968)</label><mixed-citation>
      
Matsuishi, M. and Endo, T.: Fatigue of metals subjected to varying stress,
Japan Society of Mechanical Engineers, 68, 37–40, 1968.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>Miner(1945)</label><mixed-citation>
      
Miner, M. A.: Cumulative damage in fatigue, J. Appl. Mech.,
12,  A159–A164, <a href="https://doi.org/10.1115/1.4009458" target="_blank">https://doi.org/10.1115/1.4009458</a>, 1945.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib21"><label>Musallam and Johnson(2012)</label><mixed-citation>
      
Musallam, M. and Johnson, C. M.: An efficient implementation of the rainflow
counting algorithm for life consumption estimation, IEEE T.
Reliab., 61, 978–986, <a href="https://doi.org/10.1109/TR.2012.2221040" target="_blank">https://doi.org/10.1109/TR.2012.2221040</a>, 2012.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib22"><label>National Renewable Energy Laboratory(2023)</label><mixed-citation>
      
National Renewable Energy Laboratory (NREL):   FAST v7, <a href="https://www.nrel.gov/wind/nwtc/fastv7.html" target="_blank"/>, last access: 19 April 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib23"><label>Njiri et al.(2019)Njiri, Beganovic, Do, and Söffker</label><mixed-citation>
      
Njiri, J. G., Beganovic, N., Do, M. H., and Söffker, D.: Consideration of
lifetime and fatigue load in wind turbine control, Renew. Energ., 131,
818–828, <a href="https://doi.org/10.1016/j.renene.2018.07.109" target="_blank">https://doi.org/10.1016/j.renene.2018.07.109</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib24"><label>Pegalajar-Jurado and Bredmose(2019)</label><mixed-citation>
      
Pegalajar-Jurado, A. and Bredmose, H.: Reproduction of slow-drift motions of a
floating wind turbine using second-order hydrodynamics and operational modal
analysis, Mar. Struct., 66, 178–196,
<a href="https://doi.org/10.1016/j.marstruc.2019.02.008" target="_blank">https://doi.org/10.1016/j.marstruc.2019.02.008</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib25"><label>Ragan and Manuel(2007)</label><mixed-citation>
      
Ragan, P. and Manuel, L.: Comparing estimates of wind turbine fatigue loads
using time-domain and spectral methods, Wind Eng., 31, 83–99,
<a href="https://doi.org/10.1260/030952407781494494" target="_blank">https://doi.org/10.1260/030952407781494494</a>, 2007.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib26"><label>Rinker and Dykes(2018)</label><mixed-citation>
      
Rinker, J. and Dykes, K.: WindPACT reference wind turbines, Tech. rep.,
National Renewable Energy Laboratory (NREL), Golden, Colorado, USA,
<a href="https://www.nrel.gov/docs/fy18osti/67667.pdf" target="_blank"/> (last access: 19 April 2023), 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib27"><label>Wright(2003)</label><mixed-citation>
      
Wright, A. D.: Modern control design for flexible wind turbines, University
of Colorado at Boulder, Golden, Colorado, USA,
<a href="https://www.proquest.com/openview/f8f93112e807853ee3cbc16d4c8a64e7/1?pq-origsite=gscholar&amp;cbl=18750&amp;diss=y" target="_blank"/> (last access: 19 April 2023),
2003.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib28"><label>Wright and Fingersh(2008)</label><mixed-citation>
      
Wright, A. D. and Fingersh, L. J.: Advanced control design for wind turbines
part I: Control design, implementation, and initial tests, Tech. Rep.
NREL/TP-500-42437, National Renewable Energy Laboratory (NREL), Golden,
Colorado, USA,  <a href="https://www.osti.gov/biblio/927269" target="_blank"/> (last access: 19 April 2023), 2008.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib29"><label>Yuan et al.(2015)Yuan, Li, Huang, Zhu, and Gao</label><mixed-citation>
      
Yuan, R., Li, H., Huang, H. Z., Zhu, S. P., and Gao, H.: A nonlinear fatigue
damage accumulation model considering strength degradation and its
applications to fatigue reliability analysis, Int. J. Damage Mech.,
24,  646–662,
<a href="https://doi.org/10.1177/1056789514544228" target="_blank">https://doi.org/10.1177/1056789514544228</a>, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib30"><label>Ziegler et al.(2018)Ziegler, Gonzalez, Rubert, Smolka, and
Melero</label><mixed-citation>
      
Ziegler, L., Gonzalez, E., Rubert, T., Smolka, U., and Melero, J. J.: Lifetime
extension of onshore wind turbines: A review covering Germany, Spain,
Denmark, and the UK, Renewable and Sustainable Energy Reviews, 82,
1261–1271, <a href="https://doi.org/10.1016/j.rser.2017.09.100" target="_blank">https://doi.org/10.1016/j.rser.2017.09.100</a>, 2018.

    </mixed-citation></ref-html>--></article>
