Fault Detection in Automotive Semi-active

0 downloads 0 Views 518KB Size Report
Vehicle (QoV) with suspension using a controllable shock absorber which uses ... A fault detection in an automotive semi-active ... Advanced Fault Detection and.
13AC-0123

Fault Detection in Automotive Semi-active Suspension: Experimental Results Author, co-author (Do NOT enter this information. It will be pulled from participant tab in MyTechZone) Affiliation (Do NOT enter this information. It will be pulled from participant tab in MyTechZone) Copyright ยฉ 2013 SAE International

sensors and / or actuators. The study case is a Quarter of Vehicle (QoV) with a suspension using a configurable shock absorber. The experimental testbed uses Hardware-in-theLoop (HiL) to validate the fault detection approach. The HiL is an automotive Electronic Control Unit (ECU), and the QoV model is embedded in a Field-Programmable Gate Array (FPGA). The proposed scheme is performed for discrete time systems. The LPV approach is designed to detect a sensor fault on a vertical acceleration dynamic system. A control system for a semi-active suspension can take into account that detection and be online modified to assure the best performance in comfort, as well as messages to driver of a malfunction. The main technical contribution is the capability of fault sensor detection using a differences equation of second order which uses a measurement-based scheduling parameter.

ABSTRACT A proposal of fault detection using Linear Parameter Varying (LPV) systems is experimentally validated on embedded systems. The fault detection system is oriented to faulty sensors and / or actuators. The study case is a Quarter of Vehicle (QoV) with suspension using a controllable shock absorber which uses a SkyHook controller. The experimental testbed uses Hardware-in-the-Loop (HiL) to validate the approach. The HiL is an automotive Electronic Control Unit (ECU), and the QoV model is embedded in a FieldProgrammable Gate Array (FPGA). Results exhibit the effectiveness of the approach which is consistent with numerical simulation. These results open the application of LPV approaches to commercial vehicles since it is easy of implementation for several features such as, low computation load, lumped parameter model and available for nonlinear dynamic systems.

This paper is organized as follows. Section 2 shows a literature review. Section 3 describes the automotive suspension system. Section 4 presents the proposed solution where the Hardware-in-the-Loop including the semi-active suspension control system testbed using a Controller Area Network (CAN) is described. Section 5 discusses results. Section 6 concludes the paper.

1.-INTRODUCTION A controlled automotive suspension in actual vehicles claims to improve comfort and road holding using semi-active shock absorbers. This control system uses measurements from vertical acceleration in sprung mass, unsprung mass and suspension deflection. Hence, if a fault is presented in sensors, the suspension performance will be deteriorated. The fault detection of sensors allows to reconfigurate the system to nominal performance and/or advice driver of the system fault for maintenance. Since these control systems are implemented in computation constrained embedded systems, fault detection methods with light computation are mandatory.

2.-LITERATURE REVIEW Fast fault detection of automotive suspensions in extreme applications or situations, i.e. racing (car, bike), multiple uses of vehicles (off-road to highway changes) has opportunity areas. The detection of fault in sensors in racing application will allows to always keep an optimum damping, while in vehicles will decrease the risks of loss of stability at high speeds / maneouvres after off-road suspension operation. Consumers actually expect safety, security, reliability, comfort, etc. from their vehicles as they do from their other electronic products, [1]. Advanced Fault Detection and Isolation (FDI) methods are classified into two major groups, those which do not assume any form of model information

A fault detection in an automotive semi-active suspension using Linear Parameter Varying (LPV) systems, [9], is experimental validated on automotive embedded systems. The fault detection system is oriented to faulty 1

(process history-based methods) and those which use accurate dynamic process models (model-based methods), [10]. Basically, the FDI selection depends on the availability of a reliable model or process measurements under the considered faulty modes of operation. Faults in engineering processes belong to one of the following categories, [4]: faults in sensor and actuator called soft faults or, process faults where some plant parameters change because of physical/mechanical problems.

linear viscous friction (๐น๐‘ ), a spring force (๐น๐‘˜ ) and an electrically variable force (๐น๐‘– ) which is subject to saturations. It is herby proposed to model the effects of a damper failure as an arbitrary additive force, denoted by ๐น๐‘œ . Thus, the resulting model for the damper force (๐น๐‘€๐‘… ) can be stated as:

๐น๐‘€ = ๐ผ๐‘“ ๐‘ tanh(๐‘Ž1 ๐‘งฬ‡๐‘‘๐‘’๐‘“ + ๐‘Ž2 ๐‘ง๐‘‘๐‘’๐‘“ ) + ๐‘ 1 ๐‘งฬ‡๐‘‘๐‘’๐‘“ + ๐‘ 2 ๐‘ง๐‘‘๐‘’๐‘“ + ๐น๐‘œ (2) ๏ฟฝ๏ฟฝ๏ฟฝ ๏ฟฝ๏ฟฝ๏ฟฝ ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ ๐น๐‘–

which can be rewritten:

FDI of automotive suspensions shows model based approaches have prohibitive computation [6, 12, 3, 2, 11]. In general, the approach methodology proposes the online register of data from the condition of the sensors / actuatos using mems, and then, through an online adapation of models and observers, the fault is detected. Fault detection that uses memory and online identification are prohibitive in automotive industry. Hence efficient computing and accurate methods to perform FDI in suspensions systems are an open issue.

โˆ’๐น๐‘€๐‘… โˆ’ ๐‘˜๐‘  ๐‘ง๐‘‘๐‘’๐‘“ ๐‘˜๐‘  ๐‘ง๐‘‘๐‘’๐‘“ + ๐น๐‘€๐‘… โˆ’ ๐‘˜๐‘ก (๐‘ง๐‘ข๐‘  โˆ’ ๐‘ง๐‘Ÿ )

(3)

where ๐œŒ1 = tanh(๐‘Ž1 ๐‘งฬ‡๐‘‘๐‘’๐‘“ + ๐‘Ž2 ๐‘ง๐‘‘๐‘’๐‘“ ). Therefore a state-space representation of the QoV model is given by: ๐‘งฬˆ๐‘  ๐‘งฬ‡๐‘  ๏ฟฝ ๏ฟฝ ๐‘งฬˆ๐‘ข๐‘  ๏ฟฝ ๏ฟฝ๐‘งฬ‡๏ฟฝ ๐‘ข๐‘  ๐‘ฅฬ‡

๐‘1

=

๐‘งฬˆ๐‘  ๐‘งฬˆ ๏ฟฝ ๐‘ข๐‘  ๏ฟฝ ๐‘ง๐‘‘๐‘’๐‘“ ๏ฟฝ๏ฟฝ๏ฟฝ

(1)

where ๐‘š๐‘  and ๐‘š๐‘ข๐‘  are the sprung and unsprung mass, respectively; ๐‘ง๐‘‘๐‘’๐‘“ is the suspension deflection; ๐‘ง๐‘Ÿ represents the road profile; ๐‘˜๐‘  and ๐‘˜๐‘ก are respectively the stiffness of the suspension and the tire and ๐น๐‘€๐‘… the semiactive damping force.

๐‘ฆ

๐‘1

๐‘“๐œŒ

=

๐‘ฅฬ‡ ๐‘ฆ

๐‘1

๐‘1

โŽกโˆ’ ๐‘š๐‘  โˆ’ ๐‘š๐‘  ๐‘š๐‘  โŽค ๐‘งฬ‡๐‘  ๐‘š๐‘  โŽข1 โŽฅ ๐‘ง๐‘  0 0 0 โŽข ๐‘1 ๐‘1 ๐‘1 ๐‘1 +๐‘˜๐‘ก โŽฅ ๏ฟฝ๐‘งฬ‡ ๏ฟฝ ๐‘ข๐‘  โˆ’ โˆ’ โŽข๐‘š๐‘ข๐‘  ๐‘š๐‘ข๐‘  ๐‘š๐‘ข๐‘  ๐‘š๐‘ข๐‘  โŽฅ ๐‘ง ๏ฟฝ๏ฟฝ ๏ฟฝ ๐‘ข๐‘  โŽฃ0 โŽฆ ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ 0 1 0

The dynamic behavior of a QoV with a semiactive suspension, Figure 1, is described by the following system of ordinary differential equations: = =

๐น๐‘˜

๐น๐‘€ = ๐ผ๐‘“๐‘ ๐œŒ1 + ๐‘1 ๐‘งฬ‡๐‘‘๐‘’๐‘“ + ๐‘2 ๐‘ง๐‘‘๐‘’๐‘“ + ๐น๐‘œ

3.-AUTOMOTIVE SUSPENSION SYSTEM

๐‘š๐‘  ๐‘งฬˆ๐‘  ๐‘š๐‘ข๐‘  ๐‘งฬˆ๐‘ข๐‘ 

๐น๐‘

๐‘1

๐ด

๐‘ 1 0โŽค โŽกโˆ’ ๐‘š๐‘  0 ๐ผ โŽข 0 0โŽฅ ๏ฟฝ๐‘ง ๏ฟฝ + โŽข๐‘“0๐‘๐œŒ1 โŽฅ ๐‘Ÿ ๐‘˜๐‘ก 0โŽฅ ๏ฟฝ ๐น๐‘œ โŽข ๐‘š๐‘ข๐‘  ๐‘š๐‘ข๐‘  โŽฃ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ โŽฆ 0 0 0 ๐‘ข

๐‘1

๐ต ๐‘1

๐‘ฅ

(4)

๐‘1

๐‘งฬ‡ โŽกโˆ’ ๐‘š๐‘  โˆ’ ๐‘š๐‘  ๐‘š๐‘  โŽค ๐‘  ๐‘š๐‘  ๐‘ง โŽข ๐‘1 ๐‘1 ๐‘ ๐‘ +๐‘˜ โŽฅ ๐‘  โˆ’ 1 โˆ’ 1 ๐‘ก โŽฅ ๏ฟฝ๐‘งฬ‡๐‘ข๐‘  ๏ฟฝ โŽข๐‘š๐‘ข๐‘  ๐‘š๐‘ข๐‘  ๐‘š๐‘ข๐‘  ๐‘š๐‘ข๐‘  โŽฃ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ โŽฆ ๐‘ง๐‘ข๐‘  0 1 0 โˆ’1 ๐ถ ๐‘“๐‘ ๐œŒ1

= =

1โŽค ๐ผ โŽกโˆ’ ๐‘š๐‘  0 โŽข โŽฅ ๐‘ง ๐‘“ ๐œŒ ๐‘˜ ๐‘ก + ๐‘ 1 0โŽฅ ๏ฟฝ ๐‘Ÿ ๏ฟฝ โŽข ๐‘š๐‘ข๐‘  ๐‘š๐‘ข๐‘  ๐น๐‘œ โŽฃ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ๏ฟฝ 0 0 0โŽฆ

๐ด๐‘ฅ + ๐ต๐‘ข ๐ถ๐‘ฅ + ๐ท๐‘ข

๐ท

(5)

with ๐‘1 = ๐‘˜๐‘  + ๐‘2 .

In this application, it is proposed to implement fault detection and control algorithms for a quarter of vehicle system. The QoV model is implemented in an FPGA, National Instruments Compact RIO. The fault detection and control functionalities are implemented in a MPC5517 microcontroller. Both modules communicated via a CAN network. Figure 2 illustrates the components under consideration.

Figure 1. QoV model of a semi-active suspension

The semiactive suspension is represented by an MR damper, which can be modeled as the sum of various forces: a

The QoV model is implemented in the NI-FPGA inside the real-time module of the Compact RIO โ€œcRIO-9004โ€. The 2

chassis includes a CAN module NI-9852 composed of two low-speed CAN.

Both fault detection and control algorithms are computed inside the Freescale microcontroller. The microcontroller runs at 18 ๐‘€โ„Ž๐‘ง. It is equipped with a CAN controller in order to achieve transmission with the cRIO FPGA. Thanks to the low computational load of the algorithms, real time computations are available at the constant rate specified by the system: ๐น๐‘  = 100 ๐ป๐‘ง. Computations are made in 32 bits floating point precision.

(a)

Once all data has been received, fault detection and control algorithms are computed. Then, those data are sent to cRIO. The CAN network ensures the communication between both modules. It has been chosen to use this type of network due to its great utilization in automotive industry. The CAN network is programmed in both components to work at a constant bit rate of 100 ๐‘˜๐‘๐‘–๐‘ก๐‘ /๐‘ . The National Instrument FPGA sends the 3 sensors data coded as 32 bits fixed point data. As the way of coding data is different between cRIO and Freescale, is has been chosen to multiply the decimal data to be sent by a constant ๐›ผ = 10242 and transmit directly an integer. Then, division by ๐›ผ is made in the receiving module. Each integer ๐ผ are split in 4 bytes ๐ผ = [๐ผ1 ๐ผ2 ๐ผ3 ๐ผ4 ] as CAN network work with bytes.

(b)

After experiments, it has been shown that the transmission delay was roughly 1 ms for 1 integer. Transmitting the 3 data from cRIO to Freescale took 3 ms. Figure 3 illustrates the communication.

Figure 2. Experimental implementation (a) Photo, (b) Block diagram The model runs at the sampling frequency ๐น๐‘  = 100๐ป๐‘ง in order to achieve real-time computations. The state space model is implemented as : ๏ฟฝ

๐‘‹๐‘˜+1 ๐ด ๏ฟฝ=๏ฟฝ ๐‘Œ๐‘˜ ๐ถ

๐ต ๐‘‹๐‘˜ ๏ฟฝ๏ฟฝ ๏ฟฝ ๐ท ๐‘ˆ๐‘˜

Figure 3. Communication between modules

(6)

4.-FAULT DETECTION APPROACH

where the numerical values are shown in the Appendix. It results as output the vectors of the next state ๐‘‹๐‘˜+1 and the outputs ๐‘Œ๐‘˜ of the system. The inputs of the matrix are the states ๐‘‹๐‘˜ and the control input ๐‘ˆ๐‘˜ composed by the input current ๐ผ๐‘˜ and the road profile ๐‘ง๐‘Ÿ๐‘˜ .

This section presents the fault detection approach which is used to detect a fault on the sprung mass acceleration sensor. The aim is by making use of the system inputs and outputs, recover if the actuator is faulty or healthy. The approach is facing several difficulties such as :

Once all the outputs are computed, they are sent to the fault detector and controller in order to detect a fault and compute the input current in the system.

1. The non linearity of the inputs 2. The road profile which is unmeasured The approach presented in [8] is adapted for the system under consideration as it proposes a fault detector for LPV systems with uncertainties. Nevertheless, a simpler

All decimal numbers inside the FPGA are coded in fixed point since floating point computation is not available. 3

synthesis can considerations.

be

handled

according

to

the

The parity matrix ๐’ฒ๐‘œ exists if and only if the rank of [โ„‹ | ๐’ข๐‘ข ] is degenerated. However, it is known that the rank of the concatenation is given by :

further

The fault detection scheme is based on the parityspace approach. Since the nonlinearity of the system is only in the input matrices ๐ต and ๐ท, it is proposed to consider the nonlinearity outside the matrices, directly inside the control input. For this purpose, it is considered the system : ๐‘‹๐‘˜+1 ๐‘Œ๐‘˜

= =

๐ด๐‘‹๐‘˜ + ๐ต๐ผ ๐œŒ ๏ฟฝ ๐‘˜ ๐ผ๐‘˜ + ๐ต๐‘Ÿ ๐‘ง๐‘˜ ๐ผ๐œŒ๐‘˜

๐ด๐ถ๐‘˜ + ๐ท๐ผ ๐œŒ ๏ฟฝ ๐‘˜ ๐ผ๐‘˜ + ๐ท๐‘Ÿ ๐‘ง๐‘˜ + ๐ท๐‘“ ๐น๐‘œ

rank([โ„‹

(7)

๐‘Ÿ = ๐’ฒ โ‹… (๐‘Œ๐‘  โˆ’ ๐’ข๐ผ ๐‘ˆ๐ผ๐‘  )

The parity space approach aims to express the output of the system ๐‘Œ๐‘˜ along a horizon ๐‘ . It leads to the following system :

๐ท๐‘ฅ โŽก๐ถ๐ต โŽข ๐‘ฅ ๐ถ๐ด๐ต โŽข 2๐‘ฅ where ๐’ข๐‘ฅ = โŽข๐ถ๐ด ๐ต๐‘ฅ โŽข โ‹ฎ๐‘ฅ โŽข๐ถ๐ด๐‘ โˆ’2 ๐ต๐‘ฅ โŽฃ and

๐ถ โŽก๐ถ๐ด โŽค โŽข โŽฅ โ„‹ = โŽข๐ถ๐ด2 โŽฅ โŽขโ‹ฎ โŽฅ โŽฃ๐ถ๐ด๐‘ โˆ’1 โŽฆ

0 ๐ท๐‘ฅ ๐ถ๐ต๐‘ฅ ๐ถ๐ด๐ต๐‘ฅ โ‹ฑ ๐ถ๐ด๐‘ โˆ’3 ๐ต๐‘ฅ

โ‹ฏ 0 โ‹ฑ โ‹ฑ โ‹ฑ โ‹ฏ

โ‹ฏ โ‹ฏ โ‹ฑ โ‹ฑ โ‹ฑ ๐ถ๐ด๐ต๐‘ฅ

[5] proposed a controller to minimize the vertical motion of the sprung mass by connecting a virtual damper between the body and the sky, named Sky Hook (SH). It is a practical and feasible solution to improve comfort on primary ride driving. In practice, the adjustable damper is approximated to mimic the virtual damper. The needed damping force is:

(8)

โ‹ฏ โ‹ฏ โ‹ฑ โ‹ฑ โ‹ฑ ๐ถ๐ต๐‘ฅ

0 0 โŽค โŽฅ โ‹ฎ โŽฅ โ‹ฎ โŽฅ 0 โŽฅ ๐ท๐‘ฅ โŽฅ โŽฆ

๐‘ ๐‘งฬ‡ ๐น๐‘†๐ป = ๏ฟฝ ๐‘†๐ป ๐‘  ๐‘๐‘š๐‘–๐‘› (๐‘งฬ‡๐‘  โˆ’ ๐‘งฬ‡๐‘ข๐‘  )

๐‘งฬ‡๐‘  (๐‘งฬ‡๐‘  โˆ’ ๐‘งฬ‡๐‘ข๐‘  ) โ‰ฅ 0 (12) ๐‘งฬ‡๐‘  (๐‘งฬ‡๐‘  โˆ’ ๐‘งฬ‡๐‘ข๐‘  ) < 0

The SH force was converted from Newtons to Amperes using a proportional relation of 5 A / 3000 N that is the maximum force that the identified model of the MR damper can achive in the span of 0 to 5 A, [7]. The electric current was saturated in the aforementioned span.

5.-EXPERIMENTAL RESULTS The road profile is composed of a sine wave. Table 1 presents the frequencies and corresponding normalized amplitudes. The real amplitudes has been multiplied by ๐ด๐‘œ = 15 ๐‘š๐‘š. First, the vehicle is in uncontrolled system, meaning the control input is ๐ผ = 0๐ด. Then, the SH controller is integrated. Several experiments corresponding to different excitation frequencies has been performed. The sprung mass

To fulfill the requirements of sensitivity to the fault and insensitivity to disturbances like the road profile, it can be synthesized a parity-matrix ๐’ฒ๐‘œ as : | ๐’ข๐‘ข ] = 0

๐‘–๐‘“ ๐‘–๐‘“

where ๐น๐‘†๐ป is the SH force, ๐‘๐‘†๐ป is the damping coefficient to achieve the comfort, ๐‘๐‘š๐‘–๐‘› is the minimum damping coefficient, ๐‘งฬ‡๐‘  , and ๐‘งฬ‡๐‘ข๐‘  are the vertical velocity of sprung and unsprung masses. This technique certainly improves the ride comfort. This controller is the reference in comfort oriented suspensions due to its practical implementation and design.

In order to detect the fault ๐น๐‘œ , a residual has to be synthesized. Two approaches has been carried out. First, an exact decoupling face to the state and the road profile has been handled. Then a sub-optimal optimization process face to road profile is presented, making use of lower computation time.

๐’ฒ๐‘œ โ‹… [โ„‹

(11)

This application allows to reduce the horizon ๐‘ , and so computations and computation time, [8].

๐‘œ

๐‘ˆ๐ผ๐‘ 

(10)

The parity matrix is synthesized considering ๐‘  = ๐‘› + 1, where ๐‘› is the number of states of the system. It results a residual given by :

Now, the system under consideration is an LTI system. The aim is to synthesize a residual sensitive to the fault ๐น๐‘œ , and insensitive to the road profile ๐‘ง๐‘Ÿ .

๐‘Œ๐‘ 

๐’ข๐‘ข ]) โ‰ค rank([โ„‹]) + rank([๐’ข๐‘ข ])

Finally, the condition of existence of the parity matrix leads to a long horizon time ๐‘ , so large computations. The next part makes use of an optimization approach in order to guarantee the orthogonality (decoupling face) to the states combines with an optimization approach in order to minimize the effect of the road profile.

๐ผ๐œŒ๐‘˜

๐น๐‘œ ๐‘ง๐‘Ÿ ๐‘ฆ ๐ผ โŽก ฬ‡ โŽค โŽก๐‘งฬ‡๐‘Ÿ โŽค โŽก๐‘ฆฬ‡ โŽค โŽก๐ผ ฬ‡ โŽค ๐น โŽข ๐‘œ โŽฅ โŽข๐‘งฬˆ โŽฅ โŽข๐‘ฆฬˆ โŽฅ โŽข โŽฅ โŽข โŽฅ โˆ’ ๐’ข๐ผ โŽข๐ผ ฬˆ โŽฅ = โ„‹๐‘ฅ + ๐’ข๐น โŽข๐น๐‘œฬˆ โŽฅ + ๐’ข๐‘ข โŽข ๐‘Ÿ โŽฅ โŽขโ‹ฎ โŽฅ โŽขโ‹ฎ โŽฅ โŽขร‚ โ‹ฎ โŽฅ โŽขโ‹ฎ โŽฅ (๐‘ ) โŽฃ๏ฟฝ โŽฃ๐‘ฆ (๐‘ ) โŽฆ โŽฃ๐‘ง๐‘Ÿ (๐‘ ) โŽฆ ๏ฟฝ๏ฟฝ๏ฟฝ ๐ผ (๐‘ ) โŽฆ โŽฃ๐น โŽฆ

|

(9)

4

accelerometer is considered as faulty and its fault has to be estimated. Table 1.Frequencies and relative amplitude for each sine wave tested. F(Hz) A(mm) F(Hz) A(mm)

0.5 1

0.57 1

0.65 1

0.76 1

0.85 1

0.97 1

1.10 0.98

1.26 0.95

1.44 0.91

1.64 0.87

1.87 0.83

2.14 0.79

F(Hz) A(mm)

2.44 0.76

2.78 0.72

3.17 0.67

3.62 0.57

4.13 0.47

4.72 0.36

F(Hz) A(mm)

5.38 0.29

6.14 0.26

7.01 0.23

8.00 0.20

9.13 0.14

10.42 0.10

F(Hz) A(mm)

11.89 0.09

13.53 0.08

15.48 0.07

17.66 0.06

20.16 0.05

23.00 0.04 Figure 5. Faulty outputs: left plots are for uncontrolled systems and right plots are for SH control system

Figure 4 compares the uncontrolled system and the SH control system. The SH controller has been implemented. Several experiments for different frequencies between ๐‘“๐‘– = 0.5 ๐ป๐‘ง to ๐‘“๐‘’ = 23 ๐ป๐‘ง has been carried out. A pseudoBode diagram shows the benefits of a SH control system against uncontrolled system for the resonance frequenct. At this frequency the sprung mass has a high peak when the suspension is soft (I = 0A), and low peak for the SH control system.

The fault estimation procedures estimate presented fault in Figure 6. This figure shows the effectiveness of the approach. The road profile effect is well attenuated, even in the suboptimal approach. The proposal approach is robust since the fault estimation is good no matter the changes in the road profile and the SH control compensation. This shows the exact decoupling face to the road profile and the manipulation of the semi-active damper.

The residual as proposed in (11) is implemented. However, only the most significant resonating frequency ๐‘“๐‘ โ‰ƒ 1.88 ๐ป๐‘ง is presented.

The proposed approach has been proved to be feasible for commercial implementation given its light computing time and high efficiency.

6.-CONCLUSIONS This paper presents an applicative fault detection combined with a Sky Hook Controller for a class of semiactive controllers. The proposed test-bed aims to highlight the performances of a fault detector and a controller for practical industrial applications. Two theoretical approaches has been presented for fault detection. For industrial implementation reasons, a simpler form based on the combination of a parityspace approach and optimization face to road profile is taken for study. The performances of the methodology highlight the interest of this approach. On the other hand, a SkyHook controller has been implemented. The combination of those low computational load algorithms allows for real time implementation, even for industrial applications.

Figure 4. Comparison of sprung mass acceleration for uncontrolled and SH control systems

Acknoledgement.

A fault of ๐น๐‘œ = +1๐‘š. ๐‘  โˆ’2 is present on the output all the time, except between ๐‘ก๐‘– = 2.4๐‘  and ๐‘ก๐‘’ = 6.6๐‘  in uncontrolled system and between ๐‘ก๐‘– = 6.2๐‘  and ๐‘ก๐‘’ = 12.1 in SH control system. The outputs are presented in Figure 5.

This work is partially supported by the PCP 03/10 CONACYT Mรฉxico and the Rรฉseau Universitaire pour l'Innovation Industrielle from CNRS France, the Autotronics 5

Research chair at Tecnolรณgico de Monterrey, and the French National Project INOVE / ANR 2010 BLAN 0308.

6.

Kim, J., Lee, H., 2011. โ€œSensor Fault Detection and Isolation Algorithm for a Continuous Damping Control Systemโ€. J. of Automobile Eng 225, 1347โ€“1364 . 7. Lozoya-Santos, J., Tudรณn-Martรญnez, J.C., MoralesMenendez, R., Ramirez-Mendoza, R., Sename, O., Dugard, L., 2011. โ€œControl Strategies for an Automotive Suspension with a MR Damperโ€ . 18 ๐‘กโ„Ž IFAC World Congress, Italy , pp:1820โ€“1825. 8. Varrier, S., D.Koenig, Martinez, J., 2012. โ€œRobust Fault Detection for Vehicle Lateral Dynamicsโ€, 51st IEEE Conf on Decision and Control, USA. 9. Varrier, S., Koenig, D., Martinez, J., 2012. โ€œA parity Space-Based Fault Detection On LPV Systems: Approach For Vehicle Lateral Dynamics Control Systemโ€, 8 ๐‘กโ„Ž IFAC-SAFEPROCESS, pp:1191-1196. 10. Venkatasubramanian, V., Rengaswamy, R., Kavuri, S., Yin, K., 2003. โ€œA Review of Process Fault Detection and Diagnosis Part I Quantitative Model-Based Methodsโ€. Computers and Chemical Eng 27, 293โ€“311 . 11. Vidal, Y., Acho, L., Pozo, F., Rodellar, J., 2010. โ€œFault Detection in Base-Isolation Systems Via a Restoring Force Observerโ€, Conf. on Control and Fault Tolerant Systems, France , pp. 94โ€“99 . 12. Wang, H., Song, G., 2011. โ€œFault Detection and Fault Tolerant Control of a Smart Base Isolation System with Magneto-Rheological Damperโ€. Smart Mater. Struct. 20, 1โ€“9 .

Figure 6. Fault estimation: upper plot shows the uncontrolled response, and lower plot, the SH control system.

Appendix ๐‘ง ๐‘งฬ‡ ๐‘งฬˆ ๐‘˜๐‘ ๐‘๐‘ ๐‘๐‘€๐‘… ๐ผ

REFERENCES 1.

2.

3.

4. 5.

Bertram, T., Bekes, F., Greul, R., Hanke, O., Hab, C., Hilgert, J., Hiller, M., ๐‘‚ฬˆttgen, O., Opgen-Rhein, P., Torlo, M., Ward, D., 2003. โ€œModelling and Simulation for Mechatronic Design in Automotive Systemsโ€. Control Eng. Practice 11, 179โ€“190 . Fischer, D., Isermann, R., 2004. โ€œMechatronic Semiactive and Active Vehicle Suspensionsโ€. Control Eng. Practice 12 , 1353โ€“1367 . Gรกspรกr, P, Szabรณ, Z, Bokor, J, 2010. โ€œLPV Design of Fault-Tolerant Control for Road Vehiclesโ€, Conf on Control and Fault Tolerant Systems, France, pp: 807โ€“812. Gertler, J., 1998. โ€œFault Detection and Diagnosis in Engineering Systemsโ€ . CRS Press, 1 ๐‘ ๐‘ก Ed. Karnopp, D, Crosby, MJ, Harwood, R, 1974 . โ€œVibration Control Using Semi-Active Force Generatorsโ€. Trans of ASME, J of Eng for Industry 96 , 619โ€“626 .

Variable

๐‘” ๐‘š๐‘  ๐‘š๐‘ข๐‘  ๐‘˜๐‘  ๐‘˜๐‘ก ๐‘“๐‘ ๐‘Ž1 ๐‘Ž2 ๐‘1 ๐‘2

6

deflection in meters deflection velocity in meters/sec deflection acceleration in meters/secยฒ passive stiffness constant in ๐‘/๐‘š, passive damping constant in ๐‘๐‘ /๐‘š, MR damping constant in (๐‘ โ‹… ๐ด)/๐‘š, applied current in A Value 9.81 500 100 118,000 170,000 600.95 37.85 22.15 2831 -7897

Unit ๐‘š. ๐‘  โˆ’2 ๐‘˜๐‘” ๐‘˜๐‘” ๐‘. ๐‘šโˆ’1 ๐‘. ๐‘šโˆ’1 ๐‘ ๐‘. ๐‘ . ๐‘šโˆ’1 ๐‘. ๐‘šโˆ’1 ๐‘. ๐‘ . ๐‘šโˆ’1 ๐‘. ๐‘šโˆ’1

Description gravity acceleration constant sprung mass unsprung mass spring rigidity constant tire rigidity constant Damper non-linear force Suspension non-linear damping Damper non-linear stifness Suspension passive damping Suspension passive stiffness