<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">IJMNTA</journal-id><journal-title-group><journal-title>International Journal of Modern Nonlinear Theory and Application</journal-title></journal-title-group><issn pub-type="epub">2167-9479</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ijmnta.2015.42012</article-id><article-id pub-id-type="publisher-id">IJMNTA-57667</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  A DFIM Sensor Faults Multi-Model Diagnosis Approach Based on an Adaptive PI Multiobserver—Experimental Validation
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>bid</surname><given-names>Aicha</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Benhamed</surname><given-names>Mouna</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Sbita</surname><given-names>Lassaad</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>National Engineering School of Gabes, Tunisia, Photovoltaic, Wind and Geothermal Systems Research, Gabes, Tunisia</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>aicha.abid@gmail.com(BA)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>15</day><month>05</month><year>2015</year></pub-date><volume>04</volume><issue>02</issue><fpage>161</fpage><lpage>178</lpage><history><date date-type="received"><day>19</day>	<month>April</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>23</month>	<year>June</year>	</date><date date-type="accepted"><day>30</day>	<month>June</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  This paper studies the problem of diagnosis strategy for a doubly fed induction motor (DFIM) sensor faults. This strategy is based on unknown input proportional integral (PI) multiobserver. Thecontribution of this paper is on one hand the creation of a new DFIM model based on multi-model approach and, on the other hand, the synthesis of an adaptive PI multi-observer. The DFIM Volt per Hertz drive system behaves as a nonlinear complex system. It consists of a DFIM powered through a controlled PWM Voltage Source Inverter (VSI). The need of a sensorless drive requires soft sensors such as estimators or observers. In particular, an adaptive Proportional-Integral multi-observer is synthesized in order to estimate the DFIM’s outputs which are affected by different faults and to generate the different residual signals symptoms of sensor fault occurrence. The convergence of the estimation error is guaranteed by using the Lyapunov’s based theory. The proposed diagnosis approach is experimentally validated on a 1 kW Induction motor. Obtained simulation results confirm that the adaptive PI multiobserver consent to accomplish the detection, isolation and fault identification tasks with high dynamic performances.
 
</p></abstract><kwd-group><kwd>Diagnosis</kwd><kwd> Doubly Fed Induction Motor</kwd><kwd> Multi-Model Approach</kwd><kwd> Adaptive PI Multi-Observer</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Doubly-fed induction machine (DFIM) has become the most competitive choice in several applications related with renewable energy especially wind energy as a generator or as a motor for industrial applications such as rolling, rail traction or even marine propulsion or pumping. Nevertheless, it arrives that this machine presents an electric or mechanical defect.</p><p>In recent years, the development of fault diagnosis techniques has become an important issue seen the continuing evolution of modern systems complexity and the increasing demand for improving the reliability and security of controlled systems [<xref ref-type="bibr" rid="scirp.57667-ref1">1</xref>] -[<xref ref-type="bibr" rid="scirp.57667-ref5">5</xref>] .</p><p>Observer-based diagnosis method is one of the wide variety of different diagnosis approaches that have already been proposed in the literature, these approaches are based on the use of an adequate model [<xref ref-type="bibr" rid="scirp.57667-ref6">6</xref>] -[<xref ref-type="bibr" rid="scirp.57667-ref10">10</xref>] . It is often difficult to synthesize a sufficient model that is able to take into account the system’s nonlinearity and complexity. Thus many different approaches have been developed to deal with this problem. As a solution, multi-model approach is one of the most widely used modeling techniques. This approach consists in representing the whole behavior of a nonlinear system by a set of simple local models. Generated sub models are then combined using validity function to contribute to the construction of the whole model [<xref ref-type="bibr" rid="scirp.57667-ref11">11</xref>] -[<xref ref-type="bibr" rid="scirp.57667-ref13">13</xref>] . Indeed, the multi- models facilitate the extension of some analysis tools that are developed in linear context to nonlinear context.</p><p>Then, faults diagnosis consists in the use of multiobserver based on decoupled multi-model structure [<xref ref-type="bibr" rid="scirp.57667-ref14">14</xref>] -[<xref ref-type="bibr" rid="scirp.57667-ref19">19</xref>] .</p><p>Within this diagnosis context, the main objective of this paper is the detection and isolation of different types of sensor faults that may affect the doubly fed induction motor.</p><p>As it is the workhorse of industry, the use of the doubly fed induction machine (DFIM) in industrial applications has grown impressively in recent years especially for variable speed applications. Thus, an increasingly growing interest is given to the implementation of a supervision process to ensure a safe application of this machine. As the DFIM is a nonlinear complex system that is subject of load disturbances, it is often difficult to synthesize a single model, therefore, a single observer that is able to detect and isolate the system’s faults. Thus the multi-model approach may be a solution to facilitate diagnosis task.</p><p>This paper treats the different steps for the study of the diagnosis of DFIM sensor faults based on multimodel approach. Starting with the decoupled multi-model modeling of the DFIM Volt per Hertz drive system which consists of a DFIM powered through a controlled PWM Voltage Source Inverter (VSI) in presence of load disturbances, next the design of an adaptive PI multi-observer that will be exploited finally for the faults detection. The adaptive PI multi-observer is synthesized after modification of the PI multi-observer in the case of variable faults. Finally, the proposed diagnosis approach is validated experimentally on an Induction motor.</p><p>Section 2 is dedicated to DFIM modeling in the dq synchronous reference frame. Section 3 deals with the DFIM modeling via decoupled multiple model approach respecting the different steps that lead to the realization of this task. In section 4 an adaptive PI multiobserver is synthesized to be used in section 5 for the detection and isolation of sensor faults that can affect the motor. Section 6 is an experimental validation of the proposed diagnosis method on a squirrel cage induction motor.</p><p>Simulation and experimental results for the DFIM multi-model modeling and the implementation of the diagnosis approach are performed by using the environment MATLAB/SIMULINK.</p></sec><sec id="s2"><title>2. Classic DFIM Modeling</title><p>The DFIM Volt per Hertz drive system consists of a DFIM powered through the grid in stator side and a PWM inverter in the rotor side [<xref ref-type="bibr" rid="scirp.57667-ref1">1</xref>] -[<xref ref-type="bibr" rid="scirp.57667-ref4">4</xref>] .</p><p>The mathematical model of the DFIM is presented here via the dq equations in the reference frame.</p><p>The equations for the stator and rotor voltage can be written as (1)-(4).</p><disp-formula id="scirp.57667-formula1484"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x5.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.57667-formula1485"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x6.png"  xlink:type="simple"/></disp-formula><p>The flux equations can be expressed as</p><disp-formula id="scirp.57667-formula1486"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x7.png"  xlink:type="simple"/></disp-formula><p>The mechanic equation is described as</p><disp-formula id="scirp.57667-formula1487"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x8.png"  xlink:type="simple"/></disp-formula><p>A Simulink model is built using the Equations (1)-(4) and the DFIM parameters in <xref ref-type="table" rid="table1">Table 1</xref>.</p><p>The DFIM model is highly nonlinear since it contains product terms such as speed w with flux φ<sub>d</sub> or φ<sub>q</sub>, consequently with current i<sub>d</sub> or i<sub>q</sub>. Thus to cope with this problem the multi-model approach is next suggested.</p></sec><sec id="s3"><title>3. The DFIM Drive System Modeling via Decoupled Multi-Model Approach</title><p>The multi-model approach consists in representing the system’s behavior with a set of local linear models. Every local model or sub model contributes to this global representation via a validity function which takes its values in {0, 1}.</p><p>The DFIM Modeling with multi-model approach is executed through a sequence of four steps which are clusters estimation, structure identification, parametric identification and local models combination.</p><p>The performance of the cluster estimation depends on the quality of data base which must be rich in information. The process inputs and outputs are acquired after application of a voltage scalar control strategy as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>An excitation produced by applying a variable amplitude high frequency signal at the desired speed loop. The acquisition phase consists in the collection of the DFIM’s output signals; the speed, the rotor currents i<sub>rd</sub> and i<sub>rq</sub>. The rotor frequency fr is considered as the system’s input signal. The multi-model modeling is applied in such a way to create the system’s model in presence of load disturbances, thus a variable Load torque is produced by a pseudo random binary signal.</p><p>Next, the input-output collected data on DFIM are clustered into several groups through a Chui’s clustering</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> The DFIM parameters</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Variable</th><th align="center" valign="middle" >Value</th></tr></thead><tr><td align="center" valign="middle" >R<sub>s</sub></td><td align="center" valign="middle" >0.05 Ω</td></tr><tr><td align="center" valign="middle" >R<sub>r</sub></td><td align="center" valign="middle" >0.38 Ω</td></tr><tr><td align="center" valign="middle" >M<sub>sr</sub></td><td align="center" valign="middle" >47.3 H</td></tr><tr><td align="center" valign="middle" >L<sub>s</sub></td><td align="center" valign="middle" >50 H</td></tr><tr><td align="center" valign="middle" >L<sub>r</sub></td><td align="center" valign="middle" >50 H</td></tr><tr><td align="center" valign="middle" >J<sub>t</sub></td><td align="center" valign="middle" >0.05 Kg・m<sup>2</sup></td></tr><tr><td align="center" valign="middle" >f<sub>t</sub></td><td align="center" valign="middle" >0.003 N・m/rad/s</td></tr><tr><td align="center" valign="middle" >N<sub>p</sub></td><td align="center" valign="middle" >2</td></tr></tbody></table></table-wrap><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> DFIM Volt per Hertz drives system</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x9.png"/></fig><p>algorithm. Then, the structure identification is performed on each cluster using instrumental determinants ration (RDI) method or the general procedure for order estimation while the parameters identification of each sub model are identified using recursive least square (RLS) method. Finally, obtained sub models are combined using the validity concept.</p><p>The different steps of multi-model modeling and implementation are performed thanks to MATLAB/ SIMULINK environment</p><p>The modeling strategy leads to a decoupled multi-model with six sub models which can be presented as follows.</p><disp-formula id="scirp.57667-formula1488"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x10.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-2340179x11.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-2340179x12.png" xlink:type="simple"/></inline-formula>and u are respectively the system’s i<sup>th</sup> state, output and input. Where n<sub>i</sub> = 3, p = 3.</p><disp-formula id="scirp.57667-formula1489"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x13.png"  xlink:type="simple"/></disp-formula><p>V<sub>i</sub> is the i<sup>th</sup> validity value which computes the i<sup>th</sup> sub model’s contribution to the creation of the system’s global output. These functions have the following properties (7).</p><disp-formula id="scirp.57667-formula1490"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x14.png"  xlink:type="simple"/></disp-formula><p>The different obtained DFIM system’s matrixes are expressed as (8)-(11).</p><disp-formula id="scirp.57667-formula1491"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x15.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57667-formula1492"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x16.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57667-formula1493"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x17.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57667-formula1494"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x18.png"  xlink:type="simple"/></disp-formula><p>By exciting the system with a variable input, the modeling results of the speed, the current i<sub>rd</sub> and the current i<sub>rq</sub> are shown in Figures 2-4. We can notice that the multi-model outputs follow with acceptable error the real outputs.</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Speed Multi-model modeling</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x19.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> i<sub>rd</sub> current Multi-model modeling</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x20.png"/></fig><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> i<sub>rq</sub> Multi-model modeling</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x21.png"/></fig><p>The normalized roots mean square modeling error NRMSE for each modeling error; the speed, i<sub>rd</sub> and i<sub>rq</sub> modeling error are calculated and given in <xref ref-type="table" rid="table2">Table 2</xref>.</p><p>The obtained multi-model is compared to a model done by RLS method, <xref ref-type="fig" rid="fig5">Figure 5</xref> approves the efficiency of the multi-model.</p><p>Next this multi-model will be exploited in the multiobserver design.</p></sec><sec id="s4"><title>4. Adaptive Proportional-Integral Multiobserver Designs</title><p>In this section we propose to study in the first part the PI multiobserver then the modified or adaptive PI multiobserver in second part of this section.</p><p>Firstly, the DFIM decoupled multi-model structure is modified in order to take into account the unknown input vector, and then exploited in the rest of this paper in order to conceive the based observer diagnosis’s strategy.</p><disp-formula id="scirp.57667-formula1495"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x22.png"  xlink:type="simple"/></disp-formula><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> NRMSE for the modeling errors</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >NRMSE</th></tr></thead><tr><td align="center" valign="middle" >Speed modeling error</td><td align="center" valign="middle" >0.0243</td></tr><tr><td align="center" valign="middle" >i<sub>rd</sub> current modeling error</td><td align="center" valign="middle" >0.0947</td></tr><tr><td align="center" valign="middle" >i<sub>rq</sub> current modeling error</td><td align="center" valign="middle" >0.0214</td></tr></tbody></table></table-wrap><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Multi-model, System and RLS outputs</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x23.png"/></fig><p>where E<sub>i</sub> and f identify, respectively, the impact of the unknown input on the state’s system and the unknown input vector.</p><sec id="s4_1"><title>4.1. PI Multiobserver</title><p>The PI structure is developed below in favor of achieving simultaneously estimation of both state and unknown inputs.</p><disp-formula id="scirp.57667-formula1496"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x24.png"  xlink:type="simple"/></disp-formula><p>where x<sub>i</sub> and y denote respectively the estimated state vector and output vector. V<sub>i</sub> are the validities functions calculated in the modeling phase.</p><p>This observer is known as a robust observer regarding the unknown inputs that have feebly variation. The main task of the observer design is to find out the gain matrices K<sub>I</sub> and K<sub>pi</sub>.</p><p>The PI observer uses the influence of the rebuilding output’s error with a proportional effect to estimate the system’s whereas the integral effect is used to estimate the signal of sensor or actuator’s defaults.</p><p>The sub models’ outputs y<sub>i</sub>(k), used, as modeling’s artificial signals, to represent the real system’s behavior are not exploitable to control an observer, indeed only the multi-model’s total output y(k) as it is accessible to measurement, can be designed with a system’s physical quantity. Thus an augmented system is defined below.</p><disp-formula id="scirp.57667-formula1497"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x25.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.57667-formula1498"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x26.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57667-formula1499"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x27.png"  xlink:type="simple"/></disp-formula><p>By the use of the Lyapunov’s approach defined below and taken advantage of ensuring the employed observer’s stability, and in order to guarantee the convergence of the estimation’s error, the conditions generated in forms of linear matrix inequalities (LMI) permit to compute the observer’s gains. Thus the following theorem proposed in [<xref ref-type="bibr" rid="scirp.57667-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.57667-ref15">15</xref>] suggests sufficient conditions ensuring the exponential convergence of the estimation’s error.</p><p>Theorem: The estimation error between the decoupled multi-model (12) and the PI observer (13) converges exponentially towards zero if there exists a symmetric definite positive matrix P and a matrix G verifying the LMI following:</p><disp-formula id="scirp.57667-formula1500"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x28.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.57667-formula1501"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x29.png"  xlink:type="simple"/></disp-formula><p>α is the attenuation rate which serves to quantify the convergence speed of the estimation error. Having 0 &lt; α &lt; 0.5 let to find the KI and Kp gains as:</p><disp-formula id="scirp.57667-formula1502"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x30.png"  xlink:type="simple"/></disp-formula><p>The state’s estimation e and input’s estimation errors ԑ are given by the following equation:</p><disp-formula id="scirp.57667-formula1503"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x31.png"  xlink:type="simple"/></disp-formula><p>Taking into account the fact that the unknown inputs are considered as constant or with very slow dynamics,</p><disp-formula id="scirp.57667-formula1504"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x32.png"  xlink:type="simple"/></disp-formula><p>Then</p><disp-formula id="scirp.57667-formula1505"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x33.png"  xlink:type="simple"/></disp-formula><p>Considering the definition of augmented system, the augmented error can be defined as</p><disp-formula id="scirp.57667-formula1506"><label>(23)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x34.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.57667-formula1507"><label>(24)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x35.png"  xlink:type="simple"/></disp-formula><p>The Lyapunov approach is defined as</p><disp-formula id="scirp.57667-formula1508"><label>(25)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x36.png"  xlink:type="simple"/></disp-formula><p>The exponential convergence is guaranteed if there exists a symmetric definite positive matrix P and a positive scalar α verifying the following condition:</p><disp-formula id="scirp.57667-formula1509"><label>(26)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x37.png"  xlink:type="simple"/></disp-formula><p>With</p><disp-formula id="scirp.57667-formula1510"><label>(27)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x38.png"  xlink:type="simple"/></disp-formula><p>This PI multiobserver is studied and developed by [<xref ref-type="bibr" rid="scirp.57667-ref17">17</xref>] [<xref ref-type="bibr" rid="scirp.57667-ref18">18</xref>] , however, the studied multiobserver is valid only in the case of slowly varying faults.</p></sec><sec id="s4_2"><title>4.2. Adaptive PI Multiobserver</title><p>To consider the case of variable fault, as expressed in (28), we propose to create an adaptive PI multiobserver that is based on the updating of the estimated validity functions <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-2340179x39.png" xlink:type="simple"/></inline-formula> which are computed as the Equations (29)- (31).</p><p>In fact the real system depends on the fault, so the multiobserver design should also consider this fault therefore in each instant a new validities are computed in order to guarantee always the convergence of the residue r<sub>i</sub> expressed by (29).</p><disp-formula id="scirp.57667-formula1511"><label>(28)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x40.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57667-formula1512"><label>(29)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x41.png"  xlink:type="simple"/></disp-formula><p>Then simple validities are expressed by (30).</p><disp-formula id="scirp.57667-formula1513"><label>(30)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x42.png"  xlink:type="simple"/></disp-formula><p>The final adaptive validities after normalization and reinforcement of the simple validities functions are expressed in (31).</p><disp-formula id="scirp.57667-formula1514"><label>(31)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x43.png"  xlink:type="simple"/></disp-formula><p>The Final adaptive PI multiobserver is developed in (32).</p><disp-formula id="scirp.57667-formula1515"><label>(32)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x44.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s5"><title>5. Detection and Isolation of DFIM System’s Sensors Faults</title><p>In this section, the main task to be reached is detection and isolation of the Doubly fed Induction motor’s faults by means of a designed adaptive PI multi-observer synthesized to estimate the different system outputs. Within this part the fault effect on the system state is expressed with E<sub>i</sub>, i = 1,・・・, N while M represents the faults effect on the DFIM output.</p><p>Three types of sensor faults are considered, they affect respectively the speed Ω, the i<sub>rd</sub> current and the i<sub>rq</sub> current. The faults are variable and have sinus form.</p><p>The PI multiobserver structure is equivalent to a bank of three multiobservers since</p><disp-formula id="scirp.57667-formula1516"><label>(33)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x45.png"  xlink:type="simple"/></disp-formula><p>The different faults are chosen to be occurred in the same date. In this case;</p><disp-formula id="scirp.57667-formula1517"><label>(34)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x46.png"  xlink:type="simple"/></disp-formula><p>And</p><disp-formula id="scirp.57667-formula1518"><label>(35)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x47.png"  xlink:type="simple"/></disp-formula><p>The resolution of the different LMI helps to find the matrices gains then to construct the PI multiobserver. The K<sub>Pi</sub> and K<sub>I</sub> gains are given by (36).</p><disp-formula id="scirp.57667-formula1519"><label>(36)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x48.png"  xlink:type="simple"/></disp-formula><p>The residual equations are given as follows:</p><disp-formula id="scirp.57667-formula1520"><label>(37)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x49.png"  xlink:type="simple"/></disp-formula><p>where R<sub>s</sub>, R<sub>id</sub> and R<sub>iq</sub> designate respectively the speed residual signal, i<sub>rd</sub> current residual signal and i<sub>rq</sub> current residual signal.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-2340179x50.png" xlink:type="simple"/></inline-formula>designate the i<sup>th</sup> estimated output and y<sub>i</sub> designate the i<sup>th</sup> real system output in no faulty case.</p><p>Obtained results approve the performance of the fault estimation method; the different system’s outputs y<sub>i</sub>, i = 1, ・・・, 3 follow rapidly and respectively, the real ones y<sub>i</sub> with satisfied error as shown in <xref ref-type="fig" rid="fig6">Figure 6</xref>, <xref ref-type="fig" rid="fig7">Figure 7</xref> and <xref ref-type="fig" rid="fig8">Figure 8</xref>.</p><p>The different fault evolutions approve that the estimated fault is occurred at the same time when the real fault is occurred so the detection task is verified.</p><p>The three residual signals shown in <xref ref-type="fig" rid="fig9">Figure 9</xref>, <xref ref-type="fig" rid="fig1">Figure 1</xref>0 and <xref ref-type="fig" rid="fig1">Figure 1</xref>1 follow with high precision the different sensors fault signals and approve that the faults are well identified.</p><p>The isolation task is verified since for each system output a multiobserver is synthesized, i.e. each residue is sensitive to only one sensor fault and insensitive to all other faults concerning the residual signals computed in (36).</p><p>The isolation of the different sensor faults can be summarized in the <xref ref-type="fig" rid="fig1">Figure 1</xref>2. In fact, if a speed sensor fault f1 is occurred then the residual signal Rs ≠ 0, if an i<sub>rd</sub> current sensor fault f2 is occurred then the residual signal R<sub>id</sub> ≠ 0, and if an i<sub>rq</sub> current sensor fault f3 is occurred, the residual signal R<sub>iq</sub> ≠ 0.</p><p>Considering the previous results, generated relationship can be described in a summarized table as shown in <xref ref-type="table" rid="table3">Table 3</xref>.</p></sec><sec id="s6"><title>6. Experimental Validation</title><p>In this section the proposed diagnosis approach is validated experimentally on 1 KW squirrel cage Induction motor.</p><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Speed sensor fault identification</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x51.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> ird current sensor fault identification</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x52.png"/></fig><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> irq current sensor fault identification</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x53.png"/></fig><fig id="fig9"  position="float"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> Speed residual signal Rs and speed sensor fault evolution</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x54.png"/></fig><fig id="fig10"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>0</label><caption><title> ird current residual signal Rid and current fault evolution</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x55.png"/></fig><fig id="fig11"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>1</label><caption><title> irq current residual signal Riq and current fault evolution</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x56.png"/></fig><fig id="fig12"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>2</label><caption><title> The residual sensor faults evolutions</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x57.png"/></fig><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Fault signature table</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >R<sub>s</sub></th><th align="center" valign="middle" >R<sub>id</sub></th><th align="center" valign="middle" >R<sub>iq</sub></th></tr></thead><tr><td align="center" valign="middle" >f1</td><td align="center" valign="middle" >≠0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >f2</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >≠0</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >f3</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >≠0</td></tr></tbody></table></table-wrap><p>An experimental set up drive system exposed in <xref ref-type="fig" rid="fig1">Figure 1</xref>3 is prepared to provide a set of input/output measurement with the help of MATLAB/SIMULINK and DSpace system with DS1104 controller board based on the digital signal processor (DSP) TMS320F240. The measurement of the stator current is achieved via Hall type sensors. An incremental encoder position sensor delivering 1024 pulses per revolution is used to detect The IM speed. Load torque is generated by a resistive bank fed by a DC generator.</p><p>To create parameter variation we propose to use a variable resistor that is connected in series to each phase of the motor to vary the stator resistance</p><p>The experimental set up is exploited to collect a rich database at 600 rpm with a wide range of loads and stator resistance variations.</p><p>The database is then clustered into eight clusters. Next a structural and parametric identification is performed into the obtained cluster to generate the eight local models. Finally a multi-model is created after the combination of the local models pondered by validities functions.</p><p>The obtained multi-model is next necessary to synthesize the adaptive PI multiobserver.</p><p>The resolution of the different LMI conditions helps to calculate the multiobserver gain matrices that are expressed in Equation (38).</p><disp-formula id="scirp.57667-formula1521"><label>(38)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-2340179x58.png"  xlink:type="simple"/></disp-formula><p>To test the obtained multiobserver, we inject two sensor fault, the first affect the speed sensor between t = 109 s and t = 130 s, while the second affect the current sensor between t = 14 s and 21 s.</p><fig id="fig13"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>3</label><caption><title> Experimental setup system</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x59.png"/></fig><p><xref ref-type="fig" rid="fig1">Figure 1</xref>4 and <xref ref-type="fig" rid="fig1">Figure 1</xref>5 approve that the estimated speed and current signals follow the real ones affected with the faults. <xref ref-type="fig" rid="fig1">Figure 1</xref>6 and <xref ref-type="fig" rid="fig1">Figure 1</xref>7 expose the residual signals that approve that the estimated fault follow the real ones.</p><p>According to experimental results, we can resume that the detection of the sensor fault is successfully achieved. We can notice that the value of the residual signals (speed and current residue) changes from zeros only when the fault occurs.</p><p>The identification of the faults is approved as the two residual signals follow with high precision of the two sensors fault signals.</p><fig id="fig14"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>4</label><caption><title> Speed sensor fault identification</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x60.png"/></fig><fig id="fig15"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>5</label><caption><title> Is current sensor fault identification: Real current with fault and estimated current</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x61.png"/></fig><fig id="fig16"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>6</label><caption><title> Speed residual signal Rw and speed sensor fault evolution</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x62.png"/></fig><fig id="fig17"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>7</label><caption><title> Current residual signal Ris and current sensor fault evolution</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-2340179x63.png"/></fig><p>The fault isolation is proved as for each system output an estimated output is generated. Each residue is sensitive to only one sensor fault.</p></sec><sec id="s7"><title>7. Conclusions</title><p>In this paper, a multi-model diagnosis strategy is applied to the detection and isolation of the different sensor faults that can affect the induction machine. Firstly the system’s modeling is investigated through the multi- model approach. Then considering the system’s decoupled multi-model structure an adaptive PI multi-observer is synthesized. The novel multiobserver is synthesized exploiting the classic PI multiobserver that is modified to obtain the adaptive one. The modification consists in the multiobserver validities calculation. This multiob- server is used in the fault detection and isolation of the different sensor faults that can affect the system’s outputs. The obtained experimental and simulation results performed under MATLAB/SIMULINK environment show that the applied method has an excellent capacity to describe the Induction machine under faulty case. The objectives are reached since the different computed residuals signals affirm that the detection, identification and isolation of the sensors faults are well achieved. In this paper the study is limited on simulation and without considering actuators and system faults, thus, in future work experimental study concern the induction motor, in future work we will propose to validate experimentally the diagnosis approach on a DFIM under sensor and actuator faults.</p></sec><sec id="s8"><title>Acknowledgements</title><p>This work was supported by grants from National Engineering School of Gabes, Tunisia, Photovoltaic, Wind and Geothermal Systems Research Unit, Gabes, Tunisia</p></sec></body><back><ref-list><title>References</title><ref id="scirp.57667-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Gritli, Y., Stefani, A., Rossi, C., Filippetti, F. and Chatti, A. (2011) Experimental Validation of Doubly Fed Induction Machine Electrical Faults Diagnosis under Time-Varying Conditions. Electric Power Systems Research, 81, 751-766.http://dx.doi.org/10.1016/j.epsr.2010.11.004</mixed-citation></ref><ref id="scirp.57667-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Djeghali, N., Ghanes, M., Djennoune, S. and Barbot, J.P. (2013) Sensorless Fault Tolerant Control for Induction Motors. International Journal of Control, Automation, and Systems, 11, 563-576.http://dx.doi.org/10.1007/s12555-012-9224-z</mixed-citation></ref><ref id="scirp.57667-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Saha, S., Aldeen, M. and Tan, C.P. (2011) Fault Detection in Transmission Networks of Power Systems. International Journal of Electrical Power &amp; Energy Systems, 33, 887-900. http://dx.doi.org/10.1016/j.ijepes.2010.12.026</mixed-citation></ref><ref id="scirp.57667-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Chen, J.H., Li, H.K., Sheng, D.R. and Li, W. (2015) A Hybrid Data-Driven Modeling Method on Sensor Condition Monitoring and Fault Diagnosis for Power Plants. International Journal of Electrical Power &amp; Energy Systems, 71, 274-284. http://dx.doi.org/10.1016/j.ijepes.2015.03.012</mixed-citation></ref><ref id="scirp.57667-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Darouach, M. (2012) On the Functional Observers for Linear Descriptor Systems. Systems &amp; Control Letters, 61, 427-434. http://dx.doi.org/10.1016/j.sysconle.2012.01.006</mixed-citation></ref><ref id="scirp.57667-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Manuja, S., Narasimhan, S. and Patwardhan, S.C. (2009) Unknown Input Modeling and Robust Fault Diagnosis Using Black Box Observers. Journal of Process Control, 19, 25-37. http://dx.doi.org/10.1016/j.jprocont.2008.02.004</mixed-citation></ref><ref id="scirp.57667-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Pertew, A.M., Marquez H.J. and Zhao, Q. (2005) H∞ Synthesis of Unknown Input Observers for Non-Linear Lipschitz Systems. International Journal of Control, 78, 1155-1165. http://dx.doi.org/10.1080/00207170500155488</mixed-citation></ref><ref id="scirp.57667-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Fu, Y.P., Cheng, Y.H. and Jiang, B. (2010) Robust Actuator Fault Estimation for Satellite Control System via Fuzzy Proportional Multiple-integral Observer Method. Proceedings of the 29th Chinese Control Conference, Beijing, 29-31 July 2010, 3942-3946.</mixed-citation></ref><ref id="scirp.57667-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Soren, G. and Horst, S. (2014) Diagnosis of Actuator Parameter Faults in Wind Turbines Using a Takagi-Sugeno Sliding Mode Observer. Intelligent Systems in Technical and Medical Diagnostics, Advances in Intelligent Systems and Computing, 230, 29-40.</mixed-citation></ref><ref id="scirp.57667-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Yang, X.-B., Jin, X.-Q., Du, Z.-M., Zhu, Y.-H. and Guo, Y.-B. (2013) A Hybrid Model-Based Fault Detection Strategy for Air Handling Unit Sensors. Energy and Buildings, 57, 132-143. http://dx.doi.org/10.1016/j.enbuild.2012.10.048</mixed-citation></ref><ref id="scirp.57667-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Abid, A., BenHamed, M. and Sbita, L. (2011) Induction Motor Real Time Application of Multimodel Modeling Approach. International Review of Electrical Engineering, 6, 655-660.</mixed-citation></ref><ref id="scirp.57667-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Capron, B.D.O. and Odloak, D. (2013) LMI-Based Multi-Model Predictive Control of an Industrial C3/C4 Splitter. Journal of Control, Automation and Electrical Systems, 24, 420-429. http://dx.doi.org/10.1007/s40313-013-0050-1</mixed-citation></ref><ref id="scirp.57667-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Elfelly, N., Dieulot, J.Y. and Borne, P. (2008) A Neural Approach of Multimodel Representation of Complex Processes. International Journal of Computers, Communications &amp; Control, 3, 149-160.</mixed-citation></ref><ref id="scirp.57667-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">He, S.P. (2014) Fault Estimation for T-S Fuzzy Markovian Jumping Systems Based on the Adaptive Observer. International Journal of Control, Automation and Systems, 12, 977-985. http://dx.doi.org/10.1007/s12555-013-0350-z</mixed-citation></ref><ref id="scirp.57667-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Yu, C.H. and Choi, J.W. (2014) Interacting Multiple Model Filter-Based Distributed Target Tracking Algorithm in Underwater Wireless Sensor Networks. International Journal of Control, Automation and Systems, 12, 618-627. http://dx.doi.org/10.1007/s12555-013-0238-y</mixed-citation></ref><ref id="scirp.57667-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Ichalal, D. (2009) Estimation et diagnostic de systèmes non linéaires décrits par un modèle de Takagi-Sugeno. PhD Thesis, Institut National Polytechnique de Lorraine, Nancy.</mixed-citation></ref><ref id="scirp.57667-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Orjuela, R., Marx, B., Ragot, J. and Maquin, D. (2010) Fault Diagnosis for Nonlinear Systems Represented by Heterogeneous Multiple Models. 2010 Conference on Control and Fault Tolerant Systems, Nice, 6-8 October 2010, 600-605. http://dx.doi.org/10.1109/SYSTOL.2010.5675969</mixed-citation></ref><ref id="scirp.57667-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Hamdi, H., Rodrigues, M., Mechmeche, C., Theilliol, D. and Benhadj Braiek, N. (2012) Fault Detection and Isolation in Linear Parameter-Varying Descriptor Systems via Proportional Integral Observer. International Journal of Adaptive Control and Signal Processing, 26, 224-240. http://dx.doi.org/10.1002/acs.1260</mixed-citation></ref><ref id="scirp.57667-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Zhou, D., Spée, R. and Wallace, A.K. (1993) Laboratory Control Implementations for Doubly-Fed Machines. Proceedings of the International Conference on Industrial Electronics, Control, and Instrumentation, IECON’93, Maui, 15-19 November 1993, 1181-1186. http://dx.doi.org/10.1109/iecon.1993.339234</mixed-citation></ref></ref-list></back></article>