Model-Based FDI Applied to a Piezoelectric Active Vibration Suppression System for Smart Flexible Spacecraft

Fault Detection and Isolation (FDI) techniques have captured extensive interest and attention in modern autonomous systems; in particular, they are of foremost importance in space applications, due to their scientific relevance, cost and current inability of doing on-orbit maintenance of space systems. In this scenario, FDI strategies are required to counteract possible failure events that, if not properly handled, can reduce system performance or compromise the realization of the mission objectives. In this paper, a model-based FDI strategy is implemented onboard a satellite equipped with a very large mesh reflector on which a distributed network of smart actuators/sensors is mounted to actively counteract undesired elastic vibrations. In particular, the detection and isolation of a possible piezo-actuator failure occurring in the Active Vibration Control (AVC) system of the antenna is addressed by a bank of Unknown Input Observers (UIOs). The design of the proposed UIOs is derived by solving a Linear Matrix Inequality (LMI) problem, which provides the conditions for their existence, and it is based on the linearized 3D state-space model of the controlled spacecraft, under the assumption that all the uncertainties, exogenous disturbances and measurement noises are neglected. Furthermore, pole assignment in the sense of D-stability is integrated in the standard formulation of the UIO to guarantee an adequate transient behaviour of the observers. Finally, an extensive Monte Carlo simulation campaign is conducted to assess the effectiveness of the proposed FDI architecture and its robustness against modelling uncertainties and measurement noise.


Introduction
Modern Earth Observation (EO) satellites are required to fulfil high demanding stability and pointing performance to benefit by the high-resolution capacity of their new instruments. In addition, large flexible appendages (as antennas and solar panels) are currently equipped to such spacecraft, thus demanding for the capacity of performing fast and precise manoeuvres without provoking excessive elastic displacement of such scientific instruments. These large systems-due to their scientific relevance, cost and extremely limited possibilities of maintenance and repair-are good candidate to gain significant advantage from Fault Detection and Isolation (FDI) techniques.
A detailed knowledge of the failure process and an early detection of malfunctions or faults is a mandatory requirement for such safety critical systems a prompt awareness of failures appearance can effectively help in re-designing operations to avoid any loss of system performance or potential instability phenomena. In most applications, FDI functionalities are provided to the system via hardware redundancy [1]; however, although the use of these methods is common, the main drawbacks associated to them are the added complexity and costs resulting from the weight and volume of the redundant units. To overcome these critical disadvantages a wide variety of model-based FDI strategies, exploiting the mathematical model of the monitored system to make a fault/no-fault decision, has been extensively studied [2,3].
FDI methods have received considerable attention in the past two decades for several applications in the aeronautic 1 3 and aerospace fields; in particular, many studies and FDI techniques have been developed for the Attitude and Orbital Control Subsystem (AOCS) of spacecrafts. In [4], a robust fault detector composed by a H − ∕H ∞ dynamical filter, associated to a GRL test for decision making, and a bank of Unknown Input Observers (UIOs) was proposed for thruster fault isolation in the context of Mars Sample Return Mission. The work reported in [5] addresses too the terminal rendezvous sequence of the Mars Sample Return mission and proposes a comparative study between position and attitude model-based FDI via Eigenstructure Assignment (EA) for thrusters' FDI. In [6] a bank of interacting multiple model unscented Kalman filters is, instead, developed to detect and isolate reaction wheel failures based on a highfidelity nonlinear model of the device.
Additionally, also the application of learning-based algorithms as a viable FDI solution has received particular attention. In this regard, one can mention: the Iterative Learning observer proposed in [7] to achieve estimation of time-varying thruster faults; the study of Patton et al. [8] where a robust approach to FDI for non-linear systems using neuro-fuzzy observers as a combination of the UIO and neuro-fuzzy structures, called the Optimal Robust Disturbance Decoupling Observer (ORDDO), is developed and [9], where a hierarchical methodology consisting of three combined Neural Network FDI schemes is used for pulsed plasma thrusters failure in formation flying satellites.
However, most of such studies have been performed with regard to the problem of satellite AOCS faults and little work has been conducted in FDI for active vibration control applications in space system. This is because the flexible dynamics of continuous systems, such as beams and plates or more complex structural systems, require high-order models to ensure an accurate dynamical representation and many FDI techniques become intractable for high-order system models [10]. The objective of this paper is to address the design of a model-based FDI architecture, based on a bank of UIOs, able to detect and isolate, with adequate responsiveness, a piezo-actuator failure occurring in the Active Vibration Control (AVC) system during classical attitude manoeuvres. Section 2 presents a brief description of the basic structure and capabilities of the Unknown Input Observer (UIO), the mathematical formulation to synthetize this particular observer and its application in FDI problems. In Sect. 3, first, a Large Mesh Reflector Model (LMRM) is introduced as study case for the proposed FDI strategy. The considered Active Vibration Control (AVC) system is realized by a distributed network of five collocated and optimally placed piezoelectric sensors and actuators. Secondly the equations describing the tri-dimensional dynamics of a spacecraft equipped with flexible appendages are reported and their linearization around an operative point introduced. In Sect. 4, the structure of the bank of UIOs is specialized in the context of possible piezo-actuator failure occurring in the AVC system of the antenna. The effectiveness of the proposed FDI architecture and its robustness against modelling uncertainties and measurement noise is, then, assessed and discussed based on a Monte Carlo campaign by simulating a typical spacecraft attitude manoeuvre profile. Finally, concluding remarks on the presented method are discussed and future developments proposed in the last section.

Model-Based FDI
A model-based FDI architecture is generally composed of two-stages structure as shown in Fig. 1: (1) a residual generator (i.e. a detection filter structured in such a way to produce non-zero residual signals whenever a failure is occurring) and (2) a decision making process applied on the residual signals to discriminate the failure from false alarms caused by the presence of noise or other disturbances. This step could consist of simple threshold tests [11] or it may be a statistical decision method such as the Generalized Likelihood Ratio (GLR) test [12] or the Sequential Probability Ration Test (SPRT) [13]. Generally, the fault information extracted from this supervision module is then used to actively change the controller so that stability and acceptable performance of the entire system can be still verified.
In this section, we will firstly present the structure of the UIO and describe the conditions in terms of Linear Matrix Inequality (LMI) that provide the existence of the UIO for linear systems. Furthermore, the D-stability concept is integrated in the standard formulation of the UIO to adjust the dynamics of the observer. Afterwards, it is shown how multiple UIOs can be integrated together into a bank of observers to provide the system with FDI functionalities.

Formulation of the Linear Unknown Input Observer (UIO)
Consider the following general description of the continuous Linear Time Invariant (LTI) system: where besides x(t) , u(t) , y(t) , that are the state vector, the known input vector, and the output vector of the plant respectively as in Eq. 1, an unknown input vector v(t) , which will be the unknown fault vector later on, is present in the dynamics of the system. Without loss of generality, it is assumed that is a full column rank matrix. An observer is defined as UIO for the system in Eq. 1 if its state estimation error vector approaches zero asymptotically, regardless of the presence of the unknown input vector in the system (e.g. a certain disturbance input from which the observer's output is desired to be independent). By assuming no a priori knowledge of the unknown inputs the UIO structure, in contrast to traditional Luenberger observers, is more suitable for use in fault detection applications. Although the actual unknown input itself does not need to be known, the primary requirement for a UIO or other disturbance decouplingbased observer is to have knowledge of the input distribution matrix .
Let us consider an observer, shown in Fig. 2, with the following structure: where x is the state estimate vector,z is the state vector of the UIO and the matrices , , 1 are defined as follows [14]: (1) where denotes the identity matrix with compatible dimensions with the expressions in Eq. 3 and the matrices , 2 which are design matrices used to guarantee the desired state estimation error dynamics of the observer. The following theorem, extracted from [14], is a sufficient condition of the existence of the UIO and gives us a systematic procedure to generate the observer based on the solution of a LMI problem in Eq. 4. The previous formulation of the UIO was derived as the special (linear) case of the work done by Chen and Saif [14] for which a full-order nonlinear observer was constructed for a class of Lipschitz nonlinear systems with unknown inputs. However, direct application of the LMI condition may result in unsatisfactory dynamical behaviour of the state estimation error as stated in [15]. To overcome this problem the D-stability concept proposed in [16] can be used jointly with the previous theorem to constrain the poles of UIO (poles of matrix ), as reported in Theorem 2.

Theorem 2 The eigenvalues of the observer state space matrix
can be assigned into a prescribed region D = ⋂ n s k=1 D k if there exists a common Lyapunov matrix = T > 0 and matrices and such that the set of n s LMIs reported in Eq. 6 is simultaneously satisfied.

Fig. 2 Schematics of an UIO applied to a LTI system
where k and k are matrices of appropriate dimension defining each region D k . The symbol ⊗ represents the Kronecker product operator.
In the following study, the assigned LMI region for each UIO has been selected as the intersection of two elementary LMI region defined as follows: , that define the LMI region assignment approach, must be adjusted in such a way that the dynamics of the estimation error react quickly enough to any type of failure considered, allowing a first distinction between healthy/faulty devices.

Fault Detection and Isolation: Bank of UIOs
Similarly to [17], by exploiting the capability of the UIO to generate a state estimation error decoupled from one or more unknown input signals, it is possible to design a bank of observers so that the output of each of them is sensitive to every possible failure except one. Therefore, the UIO with the minimum estimation error signal in the sense of the L2-norm reveals that a fault has occurred in the associated actuator. Considering s possible failure happening in the system, if for each one of them in possible to verify the rank condition in Theorem 1 then a bank of s UIOs can be assembled as in Fig. 3.
The bank of UIOs takes as input the measured output of the system y(t) and the output of the controller u(t) (i.e. the undamaged input of the system) and produces as output the estimation errors e yi (t) with i = 1, 2, ..., s . As described in Sect. 1, the output of the residual generator is used for the decision-making process to determine the presence and location of the faulty component. In this work, this process starts by computing the L2-norms of e yi (t).
In particular, both detection and isolation were tested using two similar procedures, one involving the L2-norm of the signals ||e yi (t)||, i = 1, 2, ..., s and one involving the derivative of the L2-norm of the same signals ||̇e yi (t)||, i = 1, 2, ..., s . In the detection phase, one of the following checks is implemented to verify the presence of an actuator failure: o Procedure 1: each UIO is associated to a binary label H Ni , i = 1, 2, ..., s , which can be 1 if the related ||e yi || is over a predefined threshold T n for an assigned time interval (defined as confirmation time t c ) or 0 in the other cases, as described in Eq. 7.
If at least two out of the s labels are equal to one, then a failure is detected by the system. o Procedure 2: similarly to the previous procedure, if at least 2 out of all the ||̇e yi || are over a predefined threshold T n for the confirmation time t c , then a failure is detected, as described in Eq. 8.
The threshold and the confirmation time are parameters that must be selected and tuned to balance the trade-off between fast confirmation of the failure and robustness to false alarms. Their selected values usually differ between Procedure 1 and Procedure 2.
In addition, to ensure the FDI architecture has also isolation capabilities, two alternative approaches can be used: • Approach 1: Search for the minimum between all ||e yi (t)|| (or all ||̇e yi (t)|| ). However, it should be noticed that such a method can lead to a higher rate of failed detection and false alarms in presence of noise (which cause the L2-norm of the signal to randomly change at different time instants); • Approach 2: Verify which one of the H Ni (or H Di ) is not over the preselected threshold. In this case, if just one out of the s labels is under the threshold (ideally equal to zero), the failed actuators is the one related to such UIO. This method can guarantee a higher rate of isolation in case of noise and uncertainties.

Test Case: Spacecraft Equipped with a Large Mesh Reflector
To test the proposed FDI architecture, the case of a spacecraft equipped with a large mesh reflector model (LMRM) is considered. In detail, the appendage is attached to a central platform-considered to be rigid with respect to the flexible structure-in correspondence of an attachment point P 1 defined in the S/C reference frame, whose origin O is located at the centre of the launch vehicle payload adapter (see Fig. 4). The inertial properties of the platform, without considering the appendage, are listed in Table 1.

Large Mesh Reflector Structural Model
The flexible appendage has been realized so to be representative (i.e. Model Replica) of the payload of a realistic space mission. Indeed, a structural model based on available information about the AstroMesh® antenna [18,19] has been realized in a finite element commercial environment [20]. To properly benefit from a distributed active vibration control system, the properties of the developed model have been modified to achieve a mass reduction of up to 30% with respect to the original AstroMesh® reflector. The main properties of the test case structure compared with AstroMesh® data and its replica the are reported in Table 2.
To allow the in-orbit deployment of the reflector and to support the antenna during operations, an extendable boom is generally implemented on those spacecrafts equipped with such a type of scientific instrument. Therefore, an 8-m long, 20 kg boom has been included in the structural model of the LMRM, which has been chosen as test case. The first frequencies and modal shapes of the structure-compared with AstroMesh®'s ones-can be found in Table 3.
Concerning the available data, the order of the modes for the AstroMesh® clamped system is: yaw, pitch and roll, while the related frequencies are not available. Nevertheless, having compared the behaviour of the mesh reflector with the original AstroMesh® model in Table 3, it is reasonable to assume their dynamic behaviour is quite similar.

Dynamics of the Assembled Spacecraft
To simplify the design of both controller and observers, the nonlinear dynamic model has been linearized starting from the non-linear equations of motion for spacecraft attitude control problems. For brevity's sake, only the final equations are reported in this paper, as the steps leading to such formulation have already been derived in authors' previous works [21,22]. In general, a state vector X = [X G , Θ, ] can be considered as including the position of the platform center of gravity G with respect to an ECI inertial frame X G , the attitude of the body reference frame with respect to the inertial system Θ and the modal amplitudes of the equipped flexible appendage . Hence, the non-linear dynamics of the system can be written as indicated in Eq. 9.
In detail, the system total mass matrix is reported in Eq. 10 as where m b is the mass of the platform, m A i the mass of the appendage, B̃ A i × G is the skew matrix of the static moment vector of the system with respect to the platform center of gravity G, defined in the body reference frame, B̃ A i G,TOT the moment of inertia of the system with respect to G in the body frame, B A i k and B̃ A i k translation and rotation modal participation factors (coupling with the rigid motion), is the identity matrix. In addition, the matrices T and T are defined as where is a diagonal matrix listing all angular frequencies of the appendages as cantilevered to satellite and is a diagonal matrix containing the k-th damping factor k of the corresponding elastic mode. The term F G T indicates the generalized forces (forces f G , torques c G and projection of forces on the modal base f e ) applied to the spacecraft in G.
In this work, a network of smart actuators is implemented on the flexible structure. The forces exerted by the active devices on the modal base are included in the following expression: where ̃ g UΨ is the electro-mechanical coupling matrix of the piezoelectric actuators projected on the modal basis and p the input voltage vector [23]. The non-linear terms of Eq. 9 are reported below: where is the angular velocity of the system with respect of the inertial system written in the body reference frame, is the variation of the inertia tensor due to the flexibility [22].

Distributed Active Vibration Control (AVC) System
In the proposed study, the LMRM is equipped with a distributed AVC system composed of five actuators and (10)

rd
Roll mode: torsion about X axis Frequency: n.a.
Roll mode: torsion about X axis Frequency: 0.85 Hz sensors, whose properties are listed in Table 4. Their locations have been selected based on a placement procedure that made use of the controllability and observability gramians of the system to attempt to increase-from a control perspective-the relative controllability and observability of the most important system modes. More details about the placement process may be found in Authors' previous work [20]. The final configuration-position and related number of actuators and sensors-is reported in Fig. 5.
In particular a (pi)-shaped piezoelectric stack actuator (P-PSA) [24] was proposed to exert the strong actuating power of piezoelectric stack; when a driving voltage is applied to a piezoelectric stack, expansion or contraction along the axial direction of the stack will be induced because of the inverse piezoelectric effect. If the stack is constrained at its two ends to a host support, actuating moments will be generated at the constrained edges (see Fig. 6).
As reported in literature [24,25], the bending moment generated by the actuator about the local Z-axis in Fig. 6 can be described as where h a is the distance between the longitudinal axis of the actuator and the neutral plane of the passive structure, p is the voltage needed to actuate the piezo device,d 33 is the piezoelectric material coefficient, A PZ , E PZ , l PZ and n l are the area, the Young module the length and the layer of the device respectively. If we consider a beam element with two nodes m and n (three translations and three rotations each), the degrees of freedom to be considered will be where u, v, z indicate the translational degrees of freedom along the local X, Y and Z-axis in Fig. 6 respectively, while , , the rotational ones. Consequently, the electromechanical coupling matrix for one smart actuator can be written in the element local 3-D reference frame as:

Linearized Equations of Motion (LEoM)
In this study, the orbital dynamics are assumed to have a negligible effect when the orbital period is very low with respect to the attitude controller bandwidth [26,27] (i.e. when the duration of the manoeuvre is short with respect to the orbital period). Indeed, such a situation occurs in a major part of space applications, where the angular rates are not (15) V mn = u m , v m , z m, m , m , m , u n , v n , z n, n , n , n , 1 3 high enough to retain the gyroscopic stiffness term. In such conditions, the velocity, and eventually the operating point, can be assumed equal to zero. Therefore, by operating a Taylor expansion around the operating point, one can obtain where the overline sign refers to the system written in linearized form and ( ) is the matrix function of the rotation angles depending on the type of assumed rotation sequence. It can be demonstrated [28] that the mass and stiffness matrices correspond to the ones reported in Eq. 9, while the damping matrix, also multiplying the rate terms, may be expressed as depending on the Jacobian matrix of the gyroscopic term jac Therefore, according to the previous considerations, it is reasonable to assume Finally, the linearized equations can be presented in the state-space standard form as: with: Assuming the flexible structure is not subjected to forces other than the actuators actions, the input forces u will include u = [f T g , c T g , T p ] T and the matrix m =̃ g UΨ . Given that a colocated active control system is considered in this study, and the output matrix is = ov . A voltage model amplifier is used to amplify the current signal from the piezoelectric sensor to measurable output voltage signal ps . As reported in [29], in this work, a high pass filter with a lower cut-off frequency of 0.1 Hz is used to model the voltage model amplifier circuit while the cut-off at high frequency has been neglected since it has no influence in this particular sensing application. Figure 7 shows the complete closed-loop configuration of the spacecraft and the failure management unit made of a bank of five UIOs, each one of them prescribed to the isolation of the failure of one piezoelectric actuator.
To test the proposed FDI architecture, the spacecraft equipped with the flexible LMRM is controlled using a classical PID (Proportional-Integral-Derivative) attitude control system in parallel to a DVF (Direct velocity Feedback) vibration suppression system [30]. Delays equal to 0.1 s and 0.05 s are introduced respectively for the attitude command and piezo command with the transfer function matrices G a (s) and G p (s) using a 1/5 degrees Padé approximation [31].
Since the design of the UIOs is based on the state-space model of the assembled system under the assumption that all the uncertainties (as inertial properties, natural frequencies, damping and chain delays), exogenous disturbances (as disturbance forces and torques applied on the spacecraft) and measurements noise are neglected, in Sect. 4.4 a Monte Carlo simulations campaign has been carried out to evaluate the off-nominal performance of the proposed FDI approach.
To tackle the abovementioned uncertainties in an efficient manner, the uncertain state space model is put into the linear fractional representation (LFT) framework [32] as expressed in Eq. 22 and the uncertainty values are reported in Table 5. where [n , n , n ] T is the output measurement noise vector and F u (P, ) is the upper LFT of system matrix P(s) with block-diagonal matrix containing the uncertainties extracted from the uncertain state space model of Eq. 20.
Besides the plant's uncertainties, also some random white noise is introduced in the measurements, with values reported in Table 6, to obtain more realistic output data.

Simulation and Results
In this study, the algorithms used to design the FDI system and to determine its overall performances are implemented in MATLAB/SIMULINK; in detail, the parameters of the bank are initialized in the YALMIP environment [33], a MATLAB toolbox for modelling and optimization extremely suited for applications in system and control theory. The parameters of each UIO, defined in Theorem 1 of Sect. 2.1 are computed using the semidefinite programming solver SDPT3 [34] called as external optimization routine by YALMIP.

Flight Control Law and Fault Modelling
Let us consider a manoeuvre that excites the most significant elastic modes of the mesh reflector. The gains of the PID for platform attitude control have been selected to reach desired attitude state in about 150 s, as shown in Fig. 8. Figure 9 shows a comparison between the angular rates obtained while considering the structure as passive (red curve) with respect to the case in which the distributed network of smart actuators/sensors is activated (blue curve) in the first 5 s of the manueuver.
In this study, the following piezo-actuator fault scenarios have been considered: -Scenario A: a failure in the piezoelectric device is considered so that the actuator is completely switched off at a certain time instant (in this study, the failure is supposed to occur at t = 1 s), thus providing zero control on the structural system. -Scenario B: a failure in the control system is assumed; therefore, one actuator is fed with a non-zero constant control voltage (instant of the failure t = 2 s).
The two scenarios are illustrated in Figs. 10 and 11, where the nominal actuators voltages are reported in black, and the voltage of the failed actuator is represented in red.

Results: Nominal Model
First, the performance of the UIOs for both piezo-actuator fault scenarios are tested on the nominal linear plant of the system. Figures 12 and 13 show the behaviour of the L2-norm of the estimation error respectively for Scenario A and Scenario B in case the actuator #5 fails. Clearly, as described in Sect. 2.2, it possible to define both the threshold value and the confirmation time to promptly detect and isolate the failure. Furthermore, we observe that, as expected, in the nominal case the failure of actuator #5 can be perfectly decoupled from the output of the fifth observer. Similar results, not reported here for the sake of brevity, are obtained when using the L2-norm of the derivative of the estimation error.

Results (Scenario A): Uncertain Model with the Presence of Measurement Noise
As shown in Fig. 14, as soon as uncertainties on the model are considered, the perfect decoupling of the UIO from its corresponding failure is not valid anymore. Furthermore, since it is meaningful to study the failure of the Scenario A only in the first seconds of the manoeuvre (since reasonable value of voltages are given in the time interval [0, 2] seconds as shown in Fig. 10 the detection and isolation is compromised by the transient dynamics of the observer in presence of uncertainties. Similar conclusion can be drawn, although not reported, if the FDI procedure is done via the derivative of the L2-norm of the estimation error.

Results (Scenario B): Uncertain Model with the Presence of Measurement Noise
In this section, a comparison for the failure of the Scenario B is performed between the two FDI procedures (the L2-norm of the signals ||e yi ||, i = 1, 2, ..., 5 and the one involving the L2-norm of the derivative). In this case, since a failure in Scenario B can happen at any time during the operative life of the system, it has been assumed that a fault on the actuator #5 occurs at t = 2 s. It is worth noting that also in this configuration it is not possible to detect and isolate a failure if that happens in the first instants of the manoeuvre. Compared to Scenario A, it is possible to appreciate that, at the instant of the fault, each observer of the bank responds correctly to the failure. However, due to the uncertainties, the measurement noise and the very low magnitude of the response of some observers, as shown in Fig. 15, it is practically impossible to define an optimal It should be noted that doing this operation means exchanging an increase in performance in terms of sensitivity to observer failure (higher is better) with the need to lower the confirmation time thus increasing the possibility of false alarms in the process of identifying the action.
Finally, the robustness against modelling uncertainties and measurement noise has been assessed based on an extensive Monte Carlo simulation campaign. The individual failure of each actuator (failure Scenario B) is investigated by randomly sampling the uncertain state space system and the results of the analysis, in terms of percentage of successful failure isolations, successful detections (but no isolation)  Table 7. The campaign is run after an a priori manual tuning of the threshold value and confirmation time, which have been set as T n = [0.1, 0.15, 0.1, 0.1, 0.15] and t c = 0.05 s. The analysis shows that in case of a failure in the actuators #2, #3 and #5 the proposed FDI scheme manages almost all the time to isolate or at least detect (meaning that the failure can be circumscribed to a small set of actuators) the failure. In the case of individual failures on actuator #1 and #4 the results of the campaign are not optimal.

Conclusions
This paper presents the implementation of an FDI architecture capable of detecting and isolating the possible failure of a piezoelectric actuator used for active vibration control of a flexible spacecraft in orbit. The interesting aspect of this study lies in the fact that-while the literature contains a large body of work on various fault diagnostic systems for common attitude control subsystem actuators (such as reaction wheels, control moment gyroscopes and magnetorquers)-the development of fault diagnostic systems for distributed active vibration control in space structures has not been-to the authors' knowledge-yet extensively addressed.
Furthermore, while in attitude control the failure of an actuator (i.e. thruster stuck open) greatly influence the state dynamics and so the output of the observer, in the case of piezo-actuator, although a failure could bring important reduction in the performance of the spacecraft, the effect on the overall dynamics could be rather subtle and could be potentially concealed by other effect such as uncertainties and noise. It is also worth noting that, since each piezoelectric actuator has a contribution on each modal amplitude, the failure of one of them does not directly imply a loss of controllability of the flexible dynamics, whereas the elastic displacements will be still damped, with less efficiency, by the remanent set of actuators.
The simulations of the FDI architecture have shown promising results; however, it is still necessary to verify the proposed strategy in the framework of the non-linear system dynamics. Further developments could also concern the definition of a proper and systematic procedure for the tuning of the identification parameters (thresholds and confirmation time) to both improve the detection and isolation rate and decrease false alarms/no detection percentage. A possible further solution is in the investigation of adaptive thresholds method in which a time-varying threshold [35] is determined for the decision-making process. Moreover, additional studies may be carried out to re-design the FDI unit with increased sensitivity to failure and improved robustness to uncertainties, thus potentially allowing also to tackle the problem of loss of actuator's efficiency.
Funding Open access funding provided by Università degli Studi di Roma La Sapienza within the CRUI-CARE Agreement.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.