Output feedback control of pantograph–catenary system with adaptive estimation of catenary parameters

With the increasing railway vehicle speed, pantograph–catenary (PAC) system has become an important part as its incidents still stand among the principal causes of railway traffic interruption. Indeed, when a rail vehicle moves, the pantograph should constantly press against the underside of the catenary. Nonetheless, it is difficult to get around the complexity of the physical interaction between the pantograph and the contact wire, which could deteriorate the quality of the electricity transfer. Thus, PAC system performances could dramatically be reduced because of bad current collection. Therefore, in this paper, we present an output feedback solution in order to design an active control of PAC system. The proposed solution is based on the backstepping control and an adaptive observer that estimates both the (unknown) catenary parameters and the system state. All synthesis steps are given and the closed-loop analysis shows asymptotic tracking behavior regardless of the time-varying catenary stiffness. Furthermore, a numerical example shows that the PAC contact can be regulated with desired effect.


Introduction
In high-speed rail vehicle systems, the main problem is related to the interaction between the train pantograph and the catenary. Indeed, when the train runs, the pantograph deforms the catenary and oscillatory motions are induced. Moreover, when the train speed increases, these oscillations may become larger and the loss of contact between the pantograph head and the collector wire could occur. Thus, when the pantograph moves along the catenary, it is fluctuated due to propagation and reflection of the wave on the catenary, which leads to a modification in its dynamics depending on the position [1].
Actually, both the pantograph and the catenary could be damaged if the contact force is too large as this could cause a contact wire breaking and stop the current collection. On the other hand, the pantograph and the catenary could lose their contact if the contact force is too small, which could cause an electric arc and accelerate the degradation of the contact wire. Therefore, considering at least these effects, many researchers have tried to design active controllers in order to ensure that the contact force remains as constant as possible [2]. Thus, to deal with this problem, the catenary device was first modeled with a constant stiffness. Then, optimal control strategies have been proposed [3]. This solution was sufficiently effective to reject external disturbances, but finding the optimal control gains while satisfying simultaneous objectives and hard constraints was beyond the adopted strategy. Thereafter, in order to improve the active control system performances, approximate models of the pantograph-catenary (PAC) system were considered [4,5]. Accordingly, many researchers assumed that the complex dynamics of the catenary might be well approximated by a linear mechanical system with space-varying lumped parameters. In this sense, the catenary parameters could be considered time-varying with a rate determined by the train speed [6]. Then, as uncertainties are present in almost all designed models, use of robust control techniques was proposed [7], while sometimes the problem was solved by tuning standard PID controllers [6]. More recently, a second-order sliding mode-based control scheme [8,9] that estimates the contact force using the measured displacements of the upper and lower pantograph frames was formally presented using the algebraic observability theory [10]. Nonetheless, this approach turns out to require more knowledge, as it requires the value of the control force applied to the upper frame, the velocity and the acceleration of the upper and lower frames, which could render the control system complex and expensive. Alternately, in other works, to perform output feedback, the contact force is evaluated by means of load cells whose measurements are compensated by accelerometers [11], which could also be very expensive mainly for their quick deterioration during the runs of the train.
In this work, our aim is to attenuate the time-varying stiffness fluctuation between the pantograph head and the contact wire. To this aim, the pantograph frames are modelled in terms of lumped masses, springs and dampers [4]. Without loss of generality, the mechanical parameters of the pantograph mechanism model (Fig. 1) are usually supposed to be constant and known. Then, to tackle this design task, we begin by constructing an adaptive observer that allows estimating both the system state and the unknown catenary equivalent stiffness. Thereafter, taking benefit of the asymptotic behavior of an adaptive observer, we synthesize a backstepping controller to ensure the output tracking of a nominal reference trajectory. In this sense, the proposed control strategy is peculiar as it joins control action with parameter estimation. Finally, analyzing the control system, we show that the proposed solution is appropriate to deal with the considered problem. Indeed, first rudiments of this work have been presented in [12] with a simplifying hypothesis that the catenary stiffness is known, which not the case is the present work.
The organization of the paper is as follows: Sect. 2 presents the PAC system modelling and the problem statement. Section 3 describes the adaptive observer. Then, Sect. 4 presents the controller synthesis, Sect. 5 gives the main results of the paper, and Sect. 6 presents a numerical example. Finally, Sect. 7 outlines some concluding remarks. To lighten the paper reading, all proofs are appended.

PAC system model and problem statement 2.1 PAC model
The pantograph and the catenary form a dynamically coupled vibrating system interacting through the contact force. Roughly, when the pantograph runs along the catenary, the variation of the catenary stiffness produces a periodic excitation that causes the pantograph vibration and leads to contact force fluctuation. As a result, the contact force is composed of a lift force that is static, and a dynamic force that depends on the vibration of the pantograph-catenary system and the vehicle speed.
Historically, a variety of catenary models are proposed in the literature from simpler models that consider only static variation of stiffness along a span to a complete finite element model (FEM) which describes the nonlinear dynamic interaction between the pantograph and the catenary system [13]. Especially, high accuracy models are required in the high-speed range, as the wave reflection becomes a major cause of contact loss [14]. Moreover, sometimes, the effect of other elements such as brackets, registration arms, and droppers is also taken into consideration. Nonetheless, without loss of generality, the simplified models with lumped and possibly time-varying parameters have been shown to be sufficiently accurate for control analysis and design purposes [6]. Thus, it turns out that a time-varying linear system can approximate the PAC dynamics with sufficient accuracy around a working configuration profile [4,7]. For this reason, to present the main idea of this paper, we consider a 2-DOF representation of the PAC system where, on one hand, the pantograph is modelled in terms of lumped masses, springs and dampers and, on the other hand, the catenary is modeled by a timevarying equivalent stiffness kðtÞ [15]. Indeed, a simplified  railway overhead contact line model, which neglects the wave propagation effects, may consider only the stiffness variation along a span. In this sense, if k max and k min are respectively the largest and smallest stiffness values in a span, the catenary average stiffness k 0 and the stiffness variation coefficient a may be approximated by Then, if we omit the stiffness variation between the vertical droppers, it is commonly assumed that the catenary equivalent stiffness can be written as where k 0 is the average stiffness, a is the stiffness variation coefficient in a span, k is the vehicle speed, and L is the span length. In addition, in this paper, we denote by b a possible phase shift at time t = 0, in the time-varying expression of the stiffness kðtÞ.
On the other hand, although the pantographs present several differences between each other, the two degrees of freedom (2-DOF) lumped-parameters (Figs. 1 and 2) may be considered as a reference model.
Roughly speaking, when the railway vehicle is running, this vibration system is in contact and its dynamics could be described by the following mathematical equations: where the subscript i = 1, 2 stand for the upper and the lower frame respectively, z i is the vertical coordinate of the pantograph frames, _ z i and € z i denote the first and the second derivative of displacement z i , m i is the mass, k i is the frame stiffness, b i is the damping coefficient, k is the time-varying catenary stiffness, F is the contact force applied by the panhead on the catenary lower wire, and u is the external actuator control action applied on the PAC system. Now, considering the displacements (z 1 ; z 2 ) and their respective derivatives ( _ z 1 ; _ z 2 ) as the state system, it follows that: Using this choice, it follows that the pantographcatenary system could be described by the following linear time-varying representation [1]: where the time-variable t is here omitted in order to alleviate the text. Then, k has the approximated form mentioned in Eq. (2) and recalled here for clarity: with w ¼ 2pk=L: Now, using Eq. (7), the (unknown) parameters of the stiffness k can be extracted by writing it as follows: where v T :¼ 1 cos(wtÞ sin(wtÞ ½ ; h :¼ k 0 k 0 a cosðbÞ k 0 asin(bÞ

Problem statement
In this paper, our main objective is to regulate the output contact force F around its nominal reference F r despite the time-varying catenary parameters. Then, as the catenary stiffness k, the output contact force F, and the state z are assumed to be unknown, our strategy is based on an output feedback control. Thus, assuming that only the frame positions (namely, z 1 and z 2 ) are measurable, it becomes necessary to recover all the unavailable variables and parameters using an adaptive observer. Thus, in the next Section, the aim is to compute an estimate of the state z, the catenary equivalent stiffness k, and the actual contact force F.

Observer synthesis
In this section, we propose to design an observer that estimates both the system state z and the stiffness k. First, let us define the measured output vector: Now, using Eqs. (5) and (9), one can easily verify that the considered PAC system can be described by the following state-affine representation: where A is a constant matrix; the components of the vector u u ð Þ and the matrix w y ð Þ are known, uniformly bounded, continuous functions that depend on the input u (that is here assumed to be bounded) and the measured output y; z denotes the unknown system state vector and h is the unknown catenary parameters vector.
Comparing Eq. (5) with Eq. (12), it follows that: wðyÞ ¼ where I n and 0 n are respectively the ðn; nÞ identity and null matrices.
As mentioned above, the representation Eq. (11) shows that we are considering a state affine system, where the state z and the parameters vector h are both involved in affine relationships. Moreover, let us notice that the matrix w y ð Þ and the vector u u ð Þ depend on measured signals u; y ð Þ.
Remark 1 Actually, it can easily be checked that the pair (A, C) is observable. Then, there exists a bounded K such that A-KC is Hurwitz, which means that any system _ xðtÞ ¼ ðA À KCÞxðtÞ f gis exponentially stable. Consequently, in order to estimate both z and h, we propose to use an adaptive observer whose estimation state error does vanish asymptotically.

Assumption
1 The solution is persistently exciting in the sense that, for some t ! t 0 and some bounded positive definite matrix R, there exist a 1 , b 1 and T 1 such that: Then, as candidate observer for system Eq. (11), we propose to use the following one [16]: where q z and q h are sufficiently large positive constants, and R is a bounded positive definite matrix, where R is the vector of real numbers. Now, let us consider the state and the parameters estimation errors, which are defined by Then, if Assumption 1 holds, the above system Eq. (17) is an asymptotic observer for system Eq. (11), in the sense that for any set of initial conditions z(0) and hð0Þ, both e z and e h do exponentially decay to zero. Specifically, the analysis of this observer dynamics, described by Eq. (17), states the following lemma [16].
Lemma 1 Consider the system described by Eq. (11), where the parameters h and the state z are both estimated using the observer dynamics Eq. (17). Then, 9 q [ 0; 8ẑð0Þ 2 R 4 ; 8ĥð0Þ 2 R 3 , the state estimation error e z and the parameters vector estimation error e h , exponentially go to zero with a rate driven by q ¼ minðq z ; q h Þ, where the estimation errors e z and e h , defined in Eq. (19), involve any trajectoryẑ and any parameters vectorĥ associated to the input u and the measured output vector y h Output feedback control of pantograph-catenary system with adaptive estimation of catenary… 255 For sake of clarity, the proof of Lemma 1 is placed in Appendix 1.

Remark 2
(i) The above global convergence result is obtained thanks to the fact that wðyÞ is globally Lipchitz in y, and u u ð Þ is locally Lipchitz in u. This property is a direct consequence of the fact that u u ð Þ is linear in z and the signal k is bounded. (ii) Using Eqs. (11), (17) and (19), we can easily verify that:: g 21 g 22 g 31 g 32 g 41 g 42 In this sense, bearing in mind Eq. (20), the observer dynamics Eq. (17) can be rewritten as follows: where f 1 ðẑ 1 ;ẑ 2 ;ẑ 3 Þ ¼ Àðk 1 =m 1 Þẑ 1 þ ðk 1 =m 1 Þẑ 2 À ðb 1 =m 1 Þẑ 3 ;

Controller synthesis
Recall that our aim is to ensure the output regulation of the PAC system. Namely, we wish to ensure that the contact force FðtÞ remains constant, despite the fluctuation of the catenary equivalent stiffness kðtÞ. Thereafter, as we are dealing with a linear time-varying system, we use backstepping techniques [17] since this control tool could ensure a robust regulation.
Hence, as FðtÞ ¼ kðtÞz 1 ðtÞ (where kðtÞ [ 0), aiming to force the contact force FðtÞ to track constant reference F r ðtÞ turns out to make z 1 ðtÞ track the following reference signal: Now, getting benefit of the fact that the state estimation errorzðtÞ vanishes exponentially, the proposed regulator may be performed based on Eq. (21). To this end, let us assume that the output reference F r ðtÞ and the catenary equivalent stiffness kðtÞ are differentiable as many times as necessary. Roughly, this is possible as F r ðtÞ could be made continuous using an appropriate low-pass filter. In addition, the stiffness kðtÞ is continuous as it is a physical signal.
Then, using the backstepping approach, it follows from the representation Eq. (17) that the control design will include three steps.
Step 1 Introducing the output tracking error: Using (21), it follows that whereẑ 3 stands for a virtual control. Let the corresponding stabilizing function be denoted by a 1 . Then, to stabilize Eq. (25) around e 1 ¼ 0, let us consider the Lyapunov function: Then, using Eq. (25), deriving V 1 with respect to time yields This suggests that the virtual controlẑ 3 is chosen equal to a 1 , with where c 1 is positive real. Asẑ 3 is not the actual control action, it cannot be forced toẑ 3 ¼ a 1 . Then, let us retain the expression of a 1 and introduce the following error: Now, from Eqs. (25), (27) and (29), _ e 1 and _ V 1 can be rewritten as follows: Step 2 Let us notice that, as a 1 depends on measurable signals, it follows that _ a 1 does exist. Then, using Eqs. (21) and (29), _ e 2 can be computed as whereẑ 4 stands for a virtual control. Let the corresponding stabilizing function be denoted by a 2 . To stabilize Eq. (32) around e 2 ¼ 0, let us consider the Lyapunov function: Now, using Eqs. (31) and (32), deriving V 2 with respect to time yields Equation (34) suggests thatẑ 4 should be chosen equal to a 2 with where c 2 is positive real constant number. However, asẑ 4 is not the actual control action, it cannot be forced toẑ 4 ¼ a 2 .
Then, let us retain the expression of a 2 and introduce the following error Afterwards, from Eqs. (32), (34) and (36), _ e 2 and _ V 2 can be rewritten in terms of e 1 ,e 2 and e 3 as follows: Step 3 From Eqs. (21) and (36), the computation of _ e 3 yields To stabilize Eq. (39) around e 3 ¼ 0, let us consider the Lyapunov function: Using Eqs. (38) and (39), we can derive V c with respect to time as: Then, letting with c 3 [ 0, it follows that: Now, from Eqs. (39), (41), and (43), _ V c and _ e 3 can be rewritten as follows: The output feedback controller established consists of the observer dynamics Eq. (17) and the control law Eq. (43). Then, for sake of clarity, before analyzing the closed loop control system, let us define the following error vector: From Eqs. (30), (37) and (45), one gets where H c :¼ G y :¼ 5 Control system analysis Based on Eqs. (17) and (43), the performance of the output feedback controller will now be formally analyzed. First, let us use the following transformation [18]: Afterwards, the closed loop control system could be described by the following error vector: Now, the main result is summarized in the following theorem: Theorem 1 Consider the control system consisting of the PAC system described by Eq. (11) and the output feedback controller consisting of the state observer dynamics Eq. (17) with the control law Eq. (43). Then, (1) The closed-loop control system is described by the following representation: Output feedback control of pantograph-catenary system with adaptive estimation of catenary… 257 In order to show the effectiveness of the proposed control scheme and the derived results, the control system including the PAC system, the adaptive observer and the backstepping regulator, have been simulated considering mechanical parameters with common values (Table 1). Of course, these parameters are only used in the implementation of the PAC system model. For controller synthesis, all the catenary parameters (k 0 ,a, b) are obviously assumed to be unknown. The (continuous) contact force reference F r ðtÞ is obtained by filtering a square sequence F c ðtÞ switching between 0 and F N ¼ 100 N (that is a common nominal value in real context). Afterwards, using a constant low acceleration, the vehicle velocity is assumed to have reached its high speed nominal value k ¼ 360 km/h , which means that the control system operates in steady state. Then, using a simulation tool, it turns out that good behaviour of the closed loop system can be ensured with the following controller and observer parameters: The simulation results are illustrated in Figs. 3, 4 and 5. First, in Fig. 3, we can see that the contact force FðtÞ is almost confounded with the filtered reference F r ðtÞ. The corresponding closed-loop control action uðt) is given in Fig. 4 and the approximated catenary equivalent stiffness is shown in Fig. 5. Therefore, as the contact force FðtÞ is almost constant (negligible fluctuation), it turns out that the proposed control strategy could provide almost good robustness despite the stiffness variation and good tracking performances even in high speed context.

Conclusion
In this paper, we deal with the problem of contact force regulation and tracking control design with respect to the time-varying catenary parameters in active train pantographs. The proposed control strategy is based on the backstepping approach and the adaptive observer. The main result of this work ensures that the contact force applied on the catenary could be maintained almost constant despite the time-varying catenary parameters. Interestingly, a Kalman-like adaptive observer provides the estimation of the catenary parameters, which means that the use of a contact force sensor is here avoided. A numerical example and a formal detailed analysis of the closed-loop control system show that both output regulation and tracking objectives of step-like contact force references can be asymptotically guaranteed. Additional perspectives of this study include the extension to the case of parameter uncertainties and the consideration of a more complex model of the PAC system.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Table 1 Typical mechanical parameters of a 2-DOF PAC system [15] Phase shift b ¼ p=4 Output feedback control of pantograph-catenary system with adaptive estimation of catenary… 259