Intelligent Contact Force Regulation of Pantograph–Catenary Based on Novel Type-Reduction Technology

: In this paper, an intelligent control scheme is proposed to suppress vibrations between the pantograph and the catenary by regulating the contact force to a reference value, thereby achieving stable current collection. In order to reduce the computational cost, an interval Type-2 adaptive fuzzy logic control with the Moradi–Zirhohi–Lin type reduction method is applied to deal with model uncertainties and exterior interference. Based on a simpliﬁed pantograph–catenary system model, the comparative simulation results show that variation of the contact force can be attenuated and variation disturbances can be repressed simultaneously. Furthermore, in terms of computational burden, the proposed type reduction method outperforms other type reduction methods.


Introduction
Railway electrification systems that transmit electrical energy from the traction substation to the moving high-speed trains are an efficient alternative to highway and air transportation.The vital subsystem of high-speed railways is the pantograph-catenary system (PCS).When the train moves, the pantograph on its roof slides along and sends electrical current from contact wires, the so-called catenary, to the train motors.The catenary lines are fixed to the irregular-interval support point along the railway line.Moreover, an appropriate contact force between the pantograph and the catenary makes the current collection process effective.However, some main factors affect the contact force between the pantograph and the catenary.To keep the contact force stable, the system output feedback considers the stiffness of the contact wires along each span, the speed of the train, and the distance of adjacent catenary lines.Unfortunately, this may result in loss of contact, wear, and arc, which can severely damage the pantograph and the catenary, whether or not the contact force is too large or too small.Generally, the quality of electric transmission is directly dependent on the interaction performance between the pantograph and the catenary.
Some research results focused on pantograph-catenary structure analysis and experimental aspects have been published.The vertical effect of the vehicle-track vibration on the pantograph-catenary interaction was presented in [1], and both pantograph-catenary and the vehicle-track models were conducted in [2].The contact wire irregulation stochastics and effect on high-speed railway pantograph-catenary interaction was analyzed in [1].The dynamic characteristics of the pantograph and catenary components, the systems' dynamic properties, and the environmental influences on the pantograph-catenary interaction were demonstrated in the review article [3].The simulation of pantograph-catenary interaction benchmarks proposed by researchers from nine countries across Europe and Asia was accomplished in [4].Meanwhile, [5] describes line testing methods for the certification of pantographs and catenaries and covers the reliability and availability of the pantograph-overhead line system and life-cycle costs.The effect of locomotive vibrations on pantograph-catenary interactions is investigated in [6].A pantograph test rig developed for fault detection and diagnosis of high-speed pantograph dynamic behavior is provided in [7].In this paper, we focus attention on the content of pantograph-catenary dynamics and the main objective to provide a novel type reduction (TR) technology for intelligent contact force regulation of the PCS.
Therefore, to effectively decrease the uncertain fluctuation of the contact force, the dynamic uplift force which is applied to pantograph frame by actuator plays an important part in the PCS.In the previous research, three methods achieve uniformity of the contact force.One is to consider the PCS in detail, including uncertain factors.In [8], performance to evaluate the contact quality adopted the dynamic behavior of the pantograph based on a multibody dynamics theory.In [9], contact resistance temperature distribution and microstructure were analyzed to address the evolutionary process of the sliding electric contact of the PCS.Meeting maintenance requirements by using automatic intelligent detection was also analyzed.In [10], online intelligent perception of the adaptive parameter PCS status was proposed using robust positioning of the PCS image contact point.In [11], automatically extracting the contact point combined with the deep pantograph network under complex backgrounds was proposed.Another proposal was the simple and effective concept of active control of PCS.The control strategy robustly overcame external disturbances by using a linear time-varying uncertain system [12].A combination of a backstepping robust controller [13], a model predictive controller [14], and an optimal controller [15] made the fluctuation of the contact force achieve minimization in the amplitude of the oscillations with good tracking performance.To regulate the fluctuation of the contact force, it is very important to apply adequate uplift force.
By using the advantages of fuzzy logic control (FLC) [16], intelligent management systems have been used for active control of PCS [17][18][19].Although Type-1 fuzzy logic controllers (T1FLCs) have been successfully applied in many application, they cannot fully handle linguistic uncertainties and noisy training data [20].Fortunately, the interval Type-2 fuzzy logic system (IT2FLS), bound from above and below by two Type-1 fuzzy sets (T1FSs), also called the footprint of uncertainty (FOU), is an extension of the concept of T1FSs.Hence, IT2FLSs can deal with uncertain parameters, modeling vagueness, linguistic uncertainties, and unreliability of information.However, Type-2 fuzzy logic system require a large computational cost (CC).Up to now, many iterative approaches for implementing TR for interval Type-2 fuzzy sets (IT2FSs) have been proposed [21,22].Much research has been dedicated to providing non-iterative methods to make the effective implementation of the IT2FLS easier, such as the uncertainty bound (UB) method [23], the Nagar-Bardini (NB) method [24], and Begian-Melek-Mendel (BMM) algorithms [25].A novel and noniterative TR method using the Moradi-Zirkohi-Lin (MZL) algorithms [26] was created to simultaneously reduce the computational burden and ensure accuracy.
In this paper, an intelligent contact force regulator incorporating MZL algorithms is presented to deal with oscillations produced from uncertain stiffness and to govern the contact force to a pre-specified reference value.Meanwhile, comparative simulations were used to confirm that the proposed TR method outperforms other TR methods.
The rest of this paper is organized as follows: Section 2 describes the simplified PCS model and problem formulation.Section 3 develops the intelligent contact force regulation with the MZL TR method.Section 4 provides the comparative simulation results.Section 5 gives the conclusion is given.

A Brief Description of the PCS Mathematical Model and Problem Formulation
The dynamics of a PCS are expressed in two common type models: multi-body and lumped-mass.The former is complex and high-computation-time; the latter is relatively simple, standard, and composed of a kinematical chain with elastic joints.Moreover, the lumped-parameters model is discussed in much research to test the performance of the controller [13,18,19,27].Therefore, this paper focuses on the linear lumped-mass model with two degrees of freedom; the motion equation of the PCS model is shown in Figure 1 [18,19].
where F up (t) and F c (t) are the dynamic uplift force (input) and the contact force (output) of the PCS, respectively.The PCS parameters are listed in Table 1.
The dynamics of a PCS are expressed in two common type models: multi-body an lumped-mass.The former is complex and high-computation-time; the latter is relative simple, standard, and composed of a kinematical chain with elastic joints.Moreover, th lumped-parameters model is discussed in much research to test the performance of t controller [13,18,19,27].Therefore, this paper focuses on the linear lumped-mass mod with two degrees of freedom; the motion equation of the PCS model is shown in Figure [18,19].

A Brief Description of the PCS Mathematical Model and Problem Formulation
The dynamics of a PCS are expressed in two common type models: multi-body and lumped-mass.The former is complex and high-computation-time; the latter is relatively simple, standard, and composed of a kinematical chain with elastic joints.Moreover, the lumped-parameters model is discussed in much research to test the performance of the controller [13,18,19,27].Therefore, this paper focuses on the linear lumped-mass model with two degrees of freedom; the motion equation of the PCS model is shown in Figure 1 [18,19].

A Brief Description of the PCS Mathematical Model and Problem Formulation
The dynamics of a PCS are expressed in two common type models: multi-body and lumped-mass.The former is complex and high-computation-time; the latter is relatively simple, standard, and composed of a kinematical chain with elastic joints.Moreover, the lumped-parameters model is discussed in much research to test the performance of the controller [13,18,19,27].Therefore, this paper focuses on the linear lumped-mass model with two degrees of freedom; the motion equation of the PCS model is shown in Figure 1 [18,19].If state variables are defined as , and x x =  , the state-space representation of the PCS model can be given as k pan

Catenary
x cat PEER REVIEW 3 of 16

A Brief Description of the PCS Mathematical Model and Problem Formulation
The dynamics of a PCS are expressed in two common type models: multi-body and lumped-mass.The former is complex and high-computation-time; the latter is relatively simple, standard, and composed of a kinematical chain with elastic joints.Moreover, the lumped-parameters model is discussed in much research to test the performance of the controller [13,18,19,27].Therefore, this paper focuses on the linear lumped-mass model with two degrees of freedom; the motion equation of the PCS model is shown in Figure 1 [18,19].
where ( ) are the dynamic uplift force (input) and the contact force (output) of the PCS, respectively.The PCS parameters are listed in Table 1.If state variables are defined as , and If state variables are defined as x 1 = x h , x 2 = .
x 3 , the state-space representation of the PCS model can be given as . where The main objective is to design the dynamic uplift force F up (t) to regulate the contact force F c (t) to the prescribed constant value F cd (t) = 100 N. If the time-varying stiffness of the catenary influenced by the stiffnessvariation coefficient α, the operational speed of train V, and the length in a span L are considered, the contact force (5) will be rewritten as In this paper, using the novel MZL TR algorithm [26] to lower CC and to fully handle or accommodate linguistic and numerical uncertainties, the interval Type-2 adaptive fuzzy logic control (IT2AFLC) incorporating with the identifier is developed to suppress the vibration resulting from the stiffness-variation coefficient α, the operational speed of train V, and the length in a span L.
Remark: It is well known that stiffness variation and wave propagation in catenaries are the main sources of contact loss between the pantograph and catenary [28].Particularly, wave reflection becomes a major cause for the loss of contact when the speed of the rail vehicle is above 250 km/h.However, the maximum speed of the rail vehicle in our model is 70 km/h and the wave propagation effect can be neglected.

Intelligent Controller Design
The objective of the controller is to design an adequate uplift force to effectively regulate the uncertain fluctuation of the contact force F c (t) to the desired constant value F cd (t) = 100 N, thereby achieving stable current collection.Based on the IT2AFLC incorporating the interval Type-2 adaptive fuzzy neural network identifier (IT2AFNNI) to predict the contact force used to tune the design parameters of the IT2AFLC by a back propagation (BP) algorithm, the intelligent control scheme with output feedback including uncertain stiffness is as shown in Figure 2. To begin with, the proposed control configuration includes two structures, IT2AFNNI and IT2AFLC; the details of the procedures are illustrated as follows.

Training Algorithm for IT2AFNNI
The structure of IT2AFNNI is as shown in Figure 2   To begin with, the proposed control configuration includes two structures, IT2AFNNI and IT2AFLC; the details of the procedures are illustrated as follows.

Training Algorithm for IT2AFNNI
The structure of IT2AFNNI is as shown in Figure 2 and the two inputs are F up (k) and F c (k − 1).The output is the predicted contact force F cI (k) used to replace F c (k) in the regulation phase.In order to train the IT2AFNNI, the performance index (PI) is defined as The BP algorithm used to adjust the weights of the IT2AFNNI by minimizing the PI is briefly described as By using the chain rule, ( 8) can be rewritten as where e I (k) = F c (k) − F cI (k) is the training error, η I is the learning rate, and W I is the output weight.Subscript I symbolizes the identifier.

Training Algorithm for IT2AFLC
e C (k) and .
e C (k) are two inputs of the IT2AFLC, as shown in Figure 2. The output is contact force F c (k).The PI used to train the IT2AFLC to regulate the contact force is defined as To update the weights of the IT2AFLC by attenuating the PI, the BP algorithm is briefly described as Based on the chain rule, (11) can be rewritten as where e C (k) = F cd (k) − F c (k) is the regulation error, η C is the learning rate, and W C is the output weight.Subscript C denotes the controller.The system sensitivity ∂F up (k) must be known in order to update the output weights in (12) to guarantee the convergence of the IT2FLC.The IT2AFNNI can achieve this goal by replacing (k) and ( 12) can be rewritten as

Main Results
The IT2AFLC with the MZL TR method is presented to suppress the vibration between the pantograph and the catenary and to deal with the uncertainties resulted from the noisy training data, noisy measurements, and linguistic uncertainties.It is a non-iterative algorithm which can reduce CC and ensure accuracy without computing the center of IT2FLS by sorting weights.
To start with, consider a Type-2 FLS having n inputs x 1 ∈ X 1 , • • • , x n ∈ X n and one output y ∈ Y.The Type-2 fuzzy rule base consists of a collection of IF-THEN rules.The rule of a Type-2 relation between the input space X 1 × X 2 × • • • × X n and the output space Y can be expressed as where F l j are antecedent Type-2 sets (j = 1, 2, • • • , n) and W l are consequent Type-2 sets.Let f i ∈ F i ≡ f i , f i , f i and f i are the lower and upper bounds of the firing strength of the ith rule of the IT2FLSs, and y i ∈ Y i ≡ y i , y i , y i and y i are the lower and upper bounds of the interval sets considered for the consequent part of the ith rule of the IT2FLSs.
The MZL TR algorithm is described in Table 2.

Computing y l
Computing y r Match the firing strength f i with their respective weights y i Match the firing strength f i with their respective weights To summarize the above analysis, intelligent contact force regulation of PCS based on the MZL TR algorithm is proposed as follows. Procedure: Step (1): Specify the membership functions µ F up and µ F c , and adjust the output weights of the IT2AFNNI by BP algorithm (9) to predict the contact force F cI (k) of the PCS.
Step (2): Specify the membership functions µ e , µ de , and set up the rule table of the IT2AFLC.
Step (3) Incorporate the IT2AFLC with IT2AFNNI based on the MZL algorithm by replacing (k) and using the BP algorithm (13) to tune the output weights of the IT2AFLC to regulate the contact force to the prescribed level, 100N.

Simulation Results
In order to verify the regulation performance, the control algorithm proposed in this paper was and compared with T1AFLC and passive control.Furthermore, a comparative study with other TR methods is provided.Simulation results show that the advocated control schemes can not only robustly attenuate the uncertain fluctuation of the contact force, but that the CC can also be efficiently reduced.
The parameters of the PCS used in the following simulation [18] are as specified: 23 × 10 4 N/m.According to the design procedure, the design is given in the following steps.
Step (1): The membership functions µ F up and µ F c for IT2AFNNI are selected as follows and as shown in Figure 3.
upper bound The identification performance is shown in Figure 4 and we can see that the IT2AFNNI can predict the contact force very quickly by (9).The identification performance is shown in Figure 4 and we can see that the IT2AFNNI can predict the contact force very quickly by (9).The identification performance is shown in Figure 4 and we can IT2AFNNI can predict the contact force very quickly by (9).Step (2) The membership functions µ e and µ de for IT2AFLC are selected as shown in Figure 5.  Step (2) The membership functions e μ and de μ for IT2AFLC are selected as shown in Figure 5.
, with the interval weights corresponding rule consequents presented.Table 3 gives the rule table of the IT2AFLC with M = 49 fuzzy IF-THEN linguistic rules based on the input membership functions µ e and µ de , and output membership function with the interval weights corresponding rule consequents presented.Step (3) In order to obtain the appropriate uplift force F up (t), which is the output of the IT2AFLC to force the contact force to the reference level 100 N, the non-iterative MZL algorithm was used to calculate the center of the IT2FLS to reduce the CC and to ensure accuracy.In accordance with the fuzzy IF-THEN linguistic rule table, the MZL employed to achieve the regulation objective is demonstrated as follows.The most-left uplift force F up l and the most-right uplift force F up r will be derived from the following process.
(1) Set average W avg by (2) Determine weights and then match firing strengths of rules (3) Compute all possible combinations Compute Compute crisp output, uplift force F up The following two case simulations will be analyzed and discussed: Case 1: To check robustness, the comparison between the T1AFLC with interference from a white Gaussian noise with SNR 20 dB, the traditional passive control (F up (t) = 100 N), and the proposed IT2AFLC will be considered.Due to the inclusion of linguistic uncertainty and data noise, the IT2AFLC can more efficiently regulate the contact force than Type-1 adaptive fuzzy logic control (T1AFLC) and the traditional passive control, as shown in Figure 6.The dynamic uplift force is severely unsteady, i.e., much more control effort must be consumed for T1AFLC, resulting in implementation difficultly as shown in Figure 7.In contrast, the control force is stable for IT2AFLC.
Electronics 2022, 10, x FOR PEER REVIEW 10 of 16 (2) Determine weights and then match firing strengths of rules (3) Compute all possible combinations The following two case simulations will be analyzed and discussed: Case 1: To check robustness, the comparison between the T1AFLC with interference from a white Gaussian noise with SNR 20 dB, the traditional passive control ( ( ) 100 and the proposed IT2AFLC will be considered.Due to the inclusion of linguistic uncertainty and data noise, the IT2AFLC can more efficiently regulate the contact force than Type-1 adaptive fuzzy logic control (T1AFLC) and the traditional passive control, as shown in Figure 6.The dynamic uplift force is severely unsteady, i.e., much more control effort must be consumed for T1AFLC, resulting in implementation difficultly as shown in Figure 7.In contrast, the control force is stable for IT2AFLC.Case 2: To consider the influence resulting from the time-varying stiffness of the catenary, the feedback signal of contact force ( ) c F t blended with the stiffness-variation coefficient α , the operational speed of train V , and the length in a span L as ( 6), the vibration suppression efficiency (VSE) [29] will be analyzed and compared between the proposed MZL algorithm, the Karnik-Mendel (KM) algorithm [30][31][32][33], and simplified interval Type-2 Takagi-Sugeno-Kang fuzzy systems (IT2 TSK FS) [34][35][36].
First of all, the VSE is defined as where Φ is the amplitude of the steady contact force value and 0 2 ( ) is the amplitude of the contact force by traditional passive control, as shown in Figure 8.  Case 2: To consider the influence resulting from the time-varying stiffness of the catenary, the feedback signal of contact force F c (t) blended with the stiffness-variation coefficient α, the operational speed of train V, and the length in a span L as ( 6), the vibration suppression efficiency (VSE) [29] will be analyzed and compared between the proposed MZL algorithm, the Karnik-Mendel (KM) algorithm [30][31][32][33], and simplified interval Type-2 Takagi-Sugeno-Kang fuzzy systems (IT2 TSK FS) [34][35][36].
First of all, the VSE is defined as where Φ is the amplitude of the steady contact force value and Φ 0 = 2αF cd (t) is the amplitude of the contact force by traditional passive control, as shown in Figure 8. Case 2: To consider the influence resulting from the time-varying stiffness of the catenary, the feedback signal of contact force ( ) c F t blended with the stiffness-variation coefficient α , the operational speed of train V , and the length in a span L as ( 6), the vibration suppression efficiency (VSE) [29] will be analyzed and compared between the proposed MZL algorithm, the Karnik-Mendel (KM) algorithm [30][31][32][33], and simplified interval Type-2 Takagi-Sugeno-Kang fuzzy systems (IT2 TSK FS) [34][35][36].
First of all, the VSE is defined as where Φ is the amplitude of the steady contact force value and 0 2 ( ) is the amplitude of the contact force by traditional passive control, as shown in Figure 8.            Furthermore, comparative results of the VSE and the CC for α = 0.3, α = 0.4, α = 0.5, and α = 0.6, with speed V = 60 and V = 70 km/hr and L, are reported in Tables 4 and 5, respectively.To provide a better comparison, the bar-line chart of the data reported in Tables 4 and 5 is also shown in Figure 13.It can be clearly seen that the VSE are very similar for the three different TR algorithms.However, for real-time applications, the computation time is a major factor; therefore, we can notice that the proposed MZL TR algorithm is more operant than the KM algorithm and the simplified IT2 algorithm in terms of contact force regulation accuracy and computational cost.In order to test the effect of parameter variations on the performance of vibration suppression and contact force regulation, the stiffness of the pantograph shoe suspension to 120% of their nominal values given in the above simulation.Table 6 summarizes the performance changes when there were 20% ± parameter perturbations in the spring and damping coefficients.The results show excellent contact force regulation and vibration suppression performance robust with our proposed control scheme.

Conclusions
In this paper, based on the MZL TR method, an intelligent control scheme is presented to suppress the vibration between the pantograph and the catenary by regulating the contact force to a reference value.The simple MZL algorithm is used to reduce the computational time cost and to deal with model uncertainties, linguistic uncertainties, and exterior interference; stable current collection can then be achieved.Comparative simulation results show that the vibration of the contact force resulting from the time-varying stiffness of the catenary can be effectively suppressed and the proposed MZL TR algorithm is more efficient than the simplified IT2 structure algorithm, a hybrid structure of four Type-1 fuzzy systems, and the KM TR algorithm, all common iterative algorithms, in regards to accuracy and computational burden.Furthermore, the proposed control scheme will be applied in real life in future experimental implementations.

Conclusions
In this paper, based on the MZL TR method, an intelligent control scheme is presented to suppress the vibration between the pantograph and the catenary by regulating the contact force to a reference value.The simple MZL algorithm is used to reduce the computational time cost and to deal with model uncertainties, linguistic uncertainties, and exterior interference; stable current collection can then be achieved.Comparative simulation results show that the vibration of the contact force resulting from the time-varying stiffness of the catenary can be effectively suppressed and the proposed MZL TR algorithm is more efficient than the simplified IT2 structure algorithm, a hybrid structure of four Type-1 fuzzy systems, and the KM TR algorithm, all common iterative algorithms, in regards to accuracy and computational burden.Furthermore, the proposed control scheme will be applied in real life in future experimental implementations.
 , the state-spa representation of the PCS model can be given as
the dynamic uplift force (input) and the contact force (out- put) of the PCS, respectively.The PCS parameters are listed in Table1 .

Figure 1 .Table 1 .
Figure 1.The prototype of a light rail vehicle.(a) Pantograph-catenary system components; (b) lumped-mass model.Reprinted with permission from ref. [18].Copyright 2016 IEEE.Table 1. Pantograph-catenary system parameter values.Parameters Displacement Speed Mass Viscous Damping Coefficient Stiffness of Suspension the state-space representation of the PCS model can be given as k cat

16 Figure 2 .
Figure 2. Block diagram of the proposed intelligent IT2AFLC scheme.
and the two inputs are ( ) .The output is the predicted contact force ( ) cI F k used to replace ( ) c F k in

Figure 2 .
Figure 2. Block diagram of the proposed intelligent IT2AFLC scheme.

Figure 3 .F
Figure 3.The membership functions for up F μ and

Figure 3 .
Figure 3.The membership functions for µ F up and µ F c .(a) The membership function µ F up .(b) The membership function µ F c .

Figure 3 .
Figure 3.The membership functions for

Figure 5 .
Figure 5.The membership functions for e μ and de μ .(a) The membership function e μ .(b) The

Figure 5 .
Figure 5.The membership functions for µ e and µ de .(a) The membership function µ e .(b) The membership function µ de .

Figure 6 .
Figure 6.The contact force regulation performance.(a) The contact force for T1AFLC and the proposed IT2AFLC.(b) The contact force for passive control and the proposed IT2AFLC.

Figure 6 .
Figure 6.The contact force regulation performance.(a) The contact force for T1AFLC and the proposed IT2AFLC.(b) The contact force for passive control and the proposed IT2AFLC.

Figure 7 .
Figure 7.The dynamic uplift force.(a) The dynamic uplift force for T1AFLC.(b) The dynamic uplift force for proposed IT2AFLC.

Figure 8 .
Figure 8.The definition of the Φ and 0 Φ for VSE.Two different sets of the different ( 0.3, 70, 40) V L α =

Figure 7 .
Figure 7.The dynamic uplift force.(a) The dynamic uplift force for T1AFLC.(b) The dynamic uplift force for proposed IT2AFLC.

Figure 7 .
Figure 7.The dynamic uplift force.(a) The dynamic uplift force for T1AFLC.(b) The dynamic uplift force for proposed IT2AFLC.

Figure 8 .
Figure 8.The definition of the Φ and 0 Φ for VSE.Two different sets of the different ( 0.3, 70, 40) V L α =

Figure 8 .
Figure 8.The definition of the Φ and Φ 0 for VSE.Two different sets of the different (α = 0.3, V = 70, L = 40) and (α = 0.6, V = 70, L = 80) are simulated and analyzed as follows.The contact force regulation performance for IT2AFLC with different three TR methods are almost the same, as shown in Figure 9.Similarly, the vertical position, speed, and acceleration of the pantograph head and frame are definitely the same as those demonstrated in Figures10 and 11, respectively.Finally, the simplified IT2 algorithm uses larger dynamic uplift force than the proposed MZL algorithm and the KM algorithm, as illustrated in Figure12.

Figure 9 .
Figure 9.The contact force regulation performance for IT2AFLC with three different type reduction methods.

Figure 10 .
Figure 10.The resulting vertical position, velocity, and acceleration of the pantograph head.

Figure 11 .
Figure 11.The resulting vertical position, velocity, and acceleration of the pantograph frame.

Figure 9 . 16 Figure 9 .
Figure 9.The contact force regulation performance for IT2AFLC with three different type reduction methods.

Figure 10 .
Figure 10.The resulting vertical position, velocity, and acceleration of the pantograph head.

Figure 11 .
Figure 11.The resulting vertical position, velocity, and acceleration of the pantograph frame.

Figure 10 . 16 Figure 9 .
Figure 10.The resulting vertical position, velocity, and acceleration of the pantograph head.

Figure 10 .
Figure 10.The resulting vertical position, velocity, and acceleration of the pantograph head.

Figure 11 .
Figure 11.The resulting vertical position, velocity, and acceleration of the pantograph frame.

Figure 11 .
Figure 11.The resulting vertical position, velocity, and acceleration of the pantograph frame.
pan k , the pan-head suspension h k , the catenary suspension cat k , the viscous damping of the pan-head suspension h c , and the frame suspension f c were changed from 80%
are the dynamic uplift force (input) and the contact force (ou put) of the PCS, respectively.The PCS parameters are listed in Table1.
c F t

Displace- ment Speed Mass Viscous Damp- ing Coefficient Stiffness of Suspe sion Pantograph head
h x h x 

Table 3
gives the rule table of the IT2AFLC with 49 M = fuzzy IF-THEN linguistic rules based on the input membership functions e

Table 6 .
Test of performance robustness.

Table 6 .
Test of performance robustness.