# Online Estimation of the Adhesion Coefficient and Its Derivative Based on the Cascading SMC Observer

Authors: Changfan Zhang; Jian Sun; Jing He; Linfan Liu
Journal: Journal of Sensors
 (2017)

---

## Abstract

The adhesion coefficient of wheel-rail surface cannot be directly measured, so a cascaded sliding-mode observer is proposed to observe adhesion coefficient and its derivative. A kinetic model for running heavy-duty locomotive is also established. The state equation of wheel adhesion control system is derived from the equation of traction motor torque balance, and adhesion coefficient is proposed to be calculated by load torque. Then, the cascaded sliding-mode observer is designed, and its stability is justified by Lyapunov stability. Based on the equivalence control principle for sliding-mode variable structure, an algorithm to estimate adhesion coefficient and its derivative is established. The simulation and experimental results are used to verify the effectiveness of the observer with load variations or wheel-rail status changes.


---

## Body

## 1. Introduction

Traction power of heavy-duty locomotives has improved along with the rapid development of railway transportation industry. However, the traction power of locomotives is limited by wheel-rail adhesion. The actual value of locomotive wheel-rail adhesion depends on the utilization of locomotive adhesion, which is generally expressed by adhesion coefficient. The wheel-rail adhesion is highly nonlinear and varies rapidly with time because of changes in rails and creep speed. Moreover, the adhesion coefficient cannot be measured directly [1–3]. Smart sensors are a crucial part of automatic detection and control owing to their functions of supervision and diagnosis. High measuring accuracy, quick response, working stability, and data analysis and state judgment functions are all required for sensors to guarantee safety and reliability of locomotives in motion [4]. Therefore, a smart sensor is necessary to supervise and diagnose the adhesion coefficient to maintain optimal adhesion.Adhesive control algorithms can be divided into two types according to the maximum adhesion point tracking: readhesion control and optimized adhesion control [5]. The readhesion control type [6] adjusts the motor output torque through the driver’s handle to restore adhesion after idling, and obtaining a higher adhesion utilization is difficult. Meanwhile, the optimal adhesion control type is mainly the creep speed and the adhesion characteristic curve slope algorithms [7–9]. The creep speed algorithm adjusts the adhesion coefficient repeatedly to maximize the adhesion coefficient, but it requires the measurement of the locomotive speed. Thus, its implementation is difficult. The adhesion characteristic curve slope algorithm is used to let the locomotive work at the peak after obtaining the adhesion characteristic curve in real time. Obtaining the derivative of the adhesion coefficient to the creep speed in this algorithm is difficult. To solve this problem, the adhesive derivative algorithm is proposed to estimate the actual adhesion coefficient and its derivative based on the load torque and its derivative obtained by designing the observer. Takaoka et al. [10, 11] and Teramoto et al. [12] proposed observing the load torque with a reduced-dimension observer. The algorithm is simple and can satisfy the basic requirements of the system. However, it is susceptible to noise-signal interference because of its direct differential processing when obtaining the adhesion coefficient derivative. Wenli et al. [13] and Liao et al. [14] introduced an omnidirectional observer for the motor load torque and its derivative, thereby improving the system’s immunity. The algorithm is based on the assumption that the traction motor load torque does not change. However, ensuring that the motor load torque remains constant with changes in track condition, creep speed, and so on during the actual running of the locomotive is difficult.The variable structure sliding-mode control (SMC) is robust in the presence of model uncertainties, parameter variations, or external disturbances. The SMC can effectively reduce the dependence on any accurate system model, and its application in the field of state observers has widened [15–17]. This study aims to develop a smart sensor based on the cascaded sliding-mode observer to handle the complex working conditions of locomotive. The sensor determines the load torque and its derivative and then estimates the adhesion coefficient and its derivative.This paper is structured as follows: the second part describes the wheel-rail kinetic model in the locomotive running process. The third part presents the cascading SMC observer designed to observe the adhesion coefficient and its derivative. The fourth part introduces the simulation analysis. The fifth part provides the conclusions.


## 2. Description of the Wheel-Rail System

The wheels roll forward under the action of driving torqueT when the train is running. The wheels relative to the tracks produce creep as a result of elastic deformation on the wheel-rail contact surface. They also form an adhesive force on the contact surface, which is the only power for driving the wheels forward, when rolling. Figure 1 shows the simplified uniaxial wheel-rail model.Figure 1
Simplified wheel-rail model.The creep speed is defined as follows [18]:(1)vs=vd-vt.Obtaining the locomotive speed vt more accurately is hard.The adhesion coefficient is defined as follows [19]:(2)μvs=FsWg.The pulling motor motion equation is as follows [20]:(3)Jmdωmdt=Tm-TL,(4)TL=FsrRg.The kinetic equation of the wheel set is(5)T=TmRg,Ft=Tr,Mjdvddt=Ft-Fs,Mj=Jr2.The train motion equation is presented as follows:(6)Mdvtdt=Fs-Fdvt,where vd is the wheel speed, Fs is the adhesion traction, Ft is the wheel-rail driving force, M is the locomotive mass, Wg is the axle load, Jm is the equivalent moment of the motor inertia, J is the moment of inertia of the wheel set, Tm is the electromagnetic torque of the motor, TL is the motor load torque, T is the driving torque, r is the wheel set radius, Rg is the gearbox ratio, and Fd(vt) is the drag.Figure2 shows an adhesion characteristic curve, which indicates that the creep speed increases with the traction when the train is running in the creep region. The adhesion coefficient also increases until a peak is reached. The train traction enters a nonstable section when the train traction required exceeds the maximum adhesion force of the wheel-rail. The adhesion coefficient within the section sharply decreases with higher creep speed. Moreover, the adhesion force of the wheel also reduces, hence the train idling. Locomotive is prone to wheel set idling at two transitions: (1) locomotive suddenly travels from dry high-adhesion rail-wheel surface to low-adhesion rail-wheel surface and (2) driver pulls up the control handle of locomotive, causing the traction beyond maximum available traction.Figure 2
Characteristic curve of the adhesion.The peak point of the adhesion in the locomotive traction process isdμ/dvs=0, with the stable section dμ/dvs>0 on the left side and the idling section dμ/dvs<0 on the right side (Figure 2).Equation (8) is derived from (7):(7)dμdvs=dμ/dtdvs/dt,(8)dμdvs=0⟹dμdt=0.Using (8), dμ/dt and dμ/dvs are equivalent.Based on the following:(9)μ=RgWgrTLthe load torque TL can be estimated as follows to determine the peak point of the adhesion coefficient.Using (9),(10)μ˙=RgWgrT˙L.By selecting the state variables,x1=θ, x2=ω, and x3=TL, where θ is the displacement angle of the motor and ω is the motor speed, T˙L is assumed as a bounded variable. Among them, only x1 is measurable. The following state equation can be obtained as follows using (3):(11)x˙1=x2,x˙2=-1Jmx3+1JmTm,x˙3=T˙L,(12)y=x1.In literature [9, 10], the following state equation is often obtained based on the hypothesis that the rate of change in the motor load torque is 0:(13)ω˙m=-1JmTL+1JmTm,T˙L=0.However, in actual locomotive run, change rate in motor load torque being zero is difficult to ensure.


## 3. Observer Design

### 3.1. Design of Load Torque Observer

As for objectT˙L≠0 described in formula (11), only x1 is measurable, and load torque is observed using a cascaded sliding-mode observer as follows:(14)x^˙1=x^2+v1,x^˙2=1JmTm+v2,(15)y^=x^1.The SMC item is as follows:(16)v1=k1sgn⁡x1-x^1,v2=k2sgn⁡x2-x^2,where x^1, x^2 are the observed values of x1 and x2, k1, k2 are constants to be determined, and sgn⁡· is a sign function.One defines the observer bias:(17)e1=x1-x^1,e2=x2-x^2,where e1 can be obtained from x1 and x^1. And because x2 is unknown, therefore, e2 can not be directly obtained.The following bias equation can be obtained using (11) and (14) as follows:(18)e˙1=x˙1-x^˙1=x2-x^2+v1=e2-k1sgn⁡e1,(19)e˙2=x˙2-x^˙2=-1Jmx3+1JmTm-1JmTm+v2=-1Jmx3-k2sgn⁡e2.Proposition 1.
Consider system (11) and its observer (14). For any initial x0∈Rn, if ki  (i=1,2,3) are sufficiently large, the estimated load torque can be calculated by T^L=-Jmk2sgn⁡e2.Proof.
Step 1. Selecting the Lyapunov function V1=1/2e12 and considering (18), the derivative of time is(20)V˙1=e1e˙1=e1e2-k1sgn⁡e1=e1e2-k1e1≤e1e2-k1e1=e1e2-k1.
Whenk1 is large enough,(21)k1≥e2+η1
andη1 is an arbitrary positive constant. Then,(22)V˙1=e1e˙1≤-η1e1.
The system meets the conditions for the SMC. Accordingly, the system will reach the SMC plane within a limited time. After reaching the SMC plane and based on the principle of equivalence of the SMC [21],(23)e1=e˙1=0.
Using (18),(24)e2=k1sgn⁡e1.
Step 2. In selecting the Lyapunov function V2=1/2e22+V1, the derivative of time is obtained as follows:(25)V˙2=V˙1+e2e˙2.
Using (19) and (22),(26)V˙2≤e2e˙2=e2-1Jmx3-k2sgn⁡e2=-1Jme2x3-k2e2≤1Jme2x3-k2e2=e21Jmx3-k2.
Whenk2 is large enough,(27)k2≥1Jmx3+η2
andη2 is an arbitrary positive constant; then,(28)V˙2≤e2e˙2≤-η2e2.
Likewise, the following is derived using (19):(29)x^3=-Jmk2sgn⁡e2,
wheree2 is derived from (24).
This completes the proof.


### 3.2. Design of Load Torque Derivative Observer

As for objectT˙L≠0 described in formula (11), only x1 is measurable, and load torque derivative is observed using a cascaded sliding-mode observer as follows:(30)z^˙1=z^2+v3,z^˙2=-1Jmz^3+1JmTm+v4,z^˙3=v5,(31)y^=z^1.The SMC item is as follows:(32)v3=κ1sgn⁡x1-z^1,v4=κ2sgn⁡x2-z^2,v5=κ3sgn⁡x3-z^3.wherez^1, z^2, and z^3 are the observed values of x1, x2, and x3. Meanwhile, κ1, κ2, and κ3 are constants to be determined, and sgn⁡(·) is a sign function.The state estimation error is defined as(33)E1=x1-z^1,(34)E2=x2-z^2,(35)E3=x3-z^3,where E1 can be obtained from x1 and z^1. And because x2 and x3 are unknown, therefore, E2 and E3 can not be directly obtained.The following bias equation can be obtained using (11) and (30) as follows:(36)E˙1=x˙1-z^˙1=x2-z^2+v3=E2-κ1sgn⁡E1,(37)E˙2=x˙2-z^˙2=-1Jmx3+1JmTm--1Jmz^3+1JmTm+v4=-1JmE3-κ2sgn⁡E2,(38)E˙3=x˙3-z^˙3=T˙L-v5=T˙L-κ3sgn⁡E3.Proposition 2.
Consider system (11) and its observer (30) or any initial x0∈Rn, and if κi  (i=1,2,3) are sufficiently large, the estimated load torque derivative can be calculated by T˙^L=κ3sgn⁡E3.Proof.
Step 1. Selecting the Lyapunov function V3=1/2E12 and considering (36), the derivative of time is(39)V˙3=E1E˙1=E1E2-κ1sgn⁡E1=E1E2-κ1E1≤E1E2-κ1E1=E1E2-κ1.
Whenκ1 is large enough,(40)κ1≥E2+δ1
andδ1 is an arbitrary positive constant. Then,(41)V˙3=E1E˙1≤-δ1E1.
The system meets the conditions for the SMC. Accordingly, the system will reach the SMC plane within a limited time. After reaching the SMC plane and based on the principle of equivalence of the SMC [21],(42)E1=E˙1=0.
Using (36),(43)E2=κ1sgn⁡E1.
Step 2. In selecting the Lyapunov function V4=1/2E22+V3, the derivative of time is obtained as follows:(44)V˙4=V˙3+E2E˙2.
Using (37) and (41),(45)V˙4≤E2E˙2=E2-1JmE3-κ2sgn⁡E2=-1JmE2E3-κ2E2≤1JmE2E3-κ2E2=E21JmE3-κ2.
Whenκ2 is large enough,(46)κ2≥1JmE3+δ2
andδ2 is an arbitrary positive constant; then,(47)V˙4≤E2E˙2≤-δ2E2.
Likewise, the following is derived using (37):(48)E3=-Jmκ2sgn⁡E2,
whereE2 is derived from (43).
Step 3. In selecting the Lyapunov function V5=1/2E32+V4, the derivative of time is(49)V˙5=V˙4+E3E˙3.
Using (38) and (47),(50)V˙5≤E3E˙3=E3T˙L-κ3sgn⁡E3=E3T˙L-κ3E3≤E3T˙L-κ3E3=E3T˙L-κ3.
Whenκ3 is large enough,(51)κ3≥T˙L+δ3
andδ3 is an arbitrary positive constant; then(52)V˙5≤E3E˙3≤-δ3E3.
Likewise, the following is obtained using (38):(53)T˙^L=κ3sgn⁡E3,
whereE3 is derived from (48).
This completes the proof.The SMC-based observer algorithm of the adhesion coefficient and its derivative is summarized as follows:(1)
Load torque observers and its derivative are established by formulas (14) and (30), respectively.(2)
By (24), e2 can be estimated.(3)
By (29), T^L can be estimated.(4)
By (43), E2 can be estimated.(5)
By (48), E3 can be estimated.(6)
By (53), T˙^L can be estimated.The sign functionsgn⁡ei is replaced with the relay method, sgn⁡ei=ei/ei+δi, to solve the chattering problem of the SMC observer. δ1 is a very small positive constant.(7)
By (9) and (10), u and its derivative u˙ can be estimated.Hence, the adhesion coefficientu and its derivative u˙ can be obtained, with the corresponding SMC observer structured as Figure 3.Figure 3
System diagram of the derivative observer of the adhesion coefficient.


## 3.1. Design of Load Torque Observer

As for objectT˙L≠0 described in formula (11), only x1 is measurable, and load torque is observed using a cascaded sliding-mode observer as follows:(14)x^˙1=x^2+v1,x^˙2=1JmTm+v2,(15)y^=x^1.The SMC item is as follows:(16)v1=k1sgn⁡x1-x^1,v2=k2sgn⁡x2-x^2,where x^1, x^2 are the observed values of x1 and x2, k1, k2 are constants to be determined, and sgn⁡· is a sign function.One defines the observer bias:(17)e1=x1-x^1,e2=x2-x^2,where e1 can be obtained from x1 and x^1. And because x2 is unknown, therefore, e2 can not be directly obtained.The following bias equation can be obtained using (11) and (14) as follows:(18)e˙1=x˙1-x^˙1=x2-x^2+v1=e2-k1sgn⁡e1,(19)e˙2=x˙2-x^˙2=-1Jmx3+1JmTm-1JmTm+v2=-1Jmx3-k2sgn⁡e2.Proposition 1.
Consider system (11) and its observer (14). For any initial x0∈Rn, if ki  (i=1,2,3) are sufficiently large, the estimated load torque can be calculated by T^L=-Jmk2sgn⁡e2.Proof.
Step 1. Selecting the Lyapunov function V1=1/2e12 and considering (18), the derivative of time is(20)V˙1=e1e˙1=e1e2-k1sgn⁡e1=e1e2-k1e1≤e1e2-k1e1=e1e2-k1.
Whenk1 is large enough,(21)k1≥e2+η1
andη1 is an arbitrary positive constant. Then,(22)V˙1=e1e˙1≤-η1e1.
The system meets the conditions for the SMC. Accordingly, the system will reach the SMC plane within a limited time. After reaching the SMC plane and based on the principle of equivalence of the SMC [21],(23)e1=e˙1=0.
Using (18),(24)e2=k1sgn⁡e1.
Step 2. In selecting the Lyapunov function V2=1/2e22+V1, the derivative of time is obtained as follows:(25)V˙2=V˙1+e2e˙2.
Using (19) and (22),(26)V˙2≤e2e˙2=e2-1Jmx3-k2sgn⁡e2=-1Jme2x3-k2e2≤1Jme2x3-k2e2=e21Jmx3-k2.
Whenk2 is large enough,(27)k2≥1Jmx3+η2
andη2 is an arbitrary positive constant; then,(28)V˙2≤e2e˙2≤-η2e2.
Likewise, the following is derived using (19):(29)x^3=-Jmk2sgn⁡e2,
wheree2 is derived from (24).
This completes the proof.


## 3.2. Design of Load Torque Derivative Observer

As for objectT˙L≠0 described in formula (11), only x1 is measurable, and load torque derivative is observed using a cascaded sliding-mode observer as follows:(30)z^˙1=z^2+v3,z^˙2=-1Jmz^3+1JmTm+v4,z^˙3=v5,(31)y^=z^1.The SMC item is as follows:(32)v3=κ1sgn⁡x1-z^1,v4=κ2sgn⁡x2-z^2,v5=κ3sgn⁡x3-z^3.wherez^1, z^2, and z^3 are the observed values of x1, x2, and x3. Meanwhile, κ1, κ2, and κ3 are constants to be determined, and sgn⁡(·) is a sign function.The state estimation error is defined as(33)E1=x1-z^1,(34)E2=x2-z^2,(35)E3=x3-z^3,where E1 can be obtained from x1 and z^1. And because x2 and x3 are unknown, therefore, E2 and E3 can not be directly obtained.The following bias equation can be obtained using (11) and (30) as follows:(36)E˙1=x˙1-z^˙1=x2-z^2+v3=E2-κ1sgn⁡E1,(37)E˙2=x˙2-z^˙2=-1Jmx3+1JmTm--1Jmz^3+1JmTm+v4=-1JmE3-κ2sgn⁡E2,(38)E˙3=x˙3-z^˙3=T˙L-v5=T˙L-κ3sgn⁡E3.Proposition 2.
Consider system (11) and its observer (30) or any initial x0∈Rn, and if κi  (i=1,2,3) are sufficiently large, the estimated load torque derivative can be calculated by T˙^L=κ3sgn⁡E3.Proof.
Step 1. Selecting the Lyapunov function V3=1/2E12 and considering (36), the derivative of time is(39)V˙3=E1E˙1=E1E2-κ1sgn⁡E1=E1E2-κ1E1≤E1E2-κ1E1=E1E2-κ1.
Whenκ1 is large enough,(40)κ1≥E2+δ1
andδ1 is an arbitrary positive constant. Then,(41)V˙3=E1E˙1≤-δ1E1.
The system meets the conditions for the SMC. Accordingly, the system will reach the SMC plane within a limited time. After reaching the SMC plane and based on the principle of equivalence of the SMC [21],(42)E1=E˙1=0.
Using (36),(43)E2=κ1sgn⁡E1.
Step 2. In selecting the Lyapunov function V4=1/2E22+V3, the derivative of time is obtained as follows:(44)V˙4=V˙3+E2E˙2.
Using (37) and (41),(45)V˙4≤E2E˙2=E2-1JmE3-κ2sgn⁡E2=-1JmE2E3-κ2E2≤1JmE2E3-κ2E2=E21JmE3-κ2.
Whenκ2 is large enough,(46)κ2≥1JmE3+δ2
andδ2 is an arbitrary positive constant; then,(47)V˙4≤E2E˙2≤-δ2E2.
Likewise, the following is derived using (37):(48)E3=-Jmκ2sgn⁡E2,
whereE2 is derived from (43).
Step 3. In selecting the Lyapunov function V5=1/2E32+V4, the derivative of time is(49)V˙5=V˙4+E3E˙3.
Using (38) and (47),(50)V˙5≤E3E˙3=E3T˙L-κ3sgn⁡E3=E3T˙L-κ3E3≤E3T˙L-κ3E3=E3T˙L-κ3.
Whenκ3 is large enough,(51)κ3≥T˙L+δ3
andδ3 is an arbitrary positive constant; then(52)V˙5≤E3E˙3≤-δ3E3.
Likewise, the following is obtained using (38):(53)T˙^L=κ3sgn⁡E3,
whereE3 is derived from (48).
This completes the proof.The SMC-based observer algorithm of the adhesion coefficient and its derivative is summarized as follows:(1)
Load torque observers and its derivative are established by formulas (14) and (30), respectively.(2)
By (24), e2 can be estimated.(3)
By (29), T^L can be estimated.(4)
By (43), E2 can be estimated.(5)
By (48), E3 can be estimated.(6)
By (53), T˙^L can be estimated.The sign functionsgn⁡ei is replaced with the relay method, sgn⁡ei=ei/ei+δi, to solve the chattering problem of the SMC observer. δ1 is a very small positive constant.(7)
By (9) and (10), u and its derivative u˙ can be estimated.Hence, the adhesion coefficientu and its derivative u˙ can be obtained, with the corresponding SMC observer structured as Figure 3.Figure 3
System diagram of the derivative observer of the adhesion coefficient.


## 4. Simulation and Experiment

### 4.1. Simulation Results and Analysis

The adhesion coefficient and its derivative are testified by a uniaxial locomotive simulation model established using MATLAB/Simulink.Case 1.
Simulation environment for locomotive starts with dry track surface, andTm remains 6500 N·m. Figures 4 and 5 show wheel-rail adhesion coefficient and adhesion coefficient derivative, respectively. In the figures, solid red and black-dotted lines represent actual and observed values, respectively. Figures 6 and 7 show motor angular displacement and speed, respectively. In the figures, solid red, blue-dotted, and black-dotted lines represent actual value, observed value from load torque observer, and observed value from derivative observer, respectively. As shown in Figures 4 and 5, adhesion coefficient immediately increases from zero to 0.11, and adhesion coefficient derivative is zero. At 2.5 s, locomotive travels from dry to wet rail surface, and wheel-rail adhesion does not meet locomotive traction needs. When the locomotive runs to idling section, adhesion coefficient declines to zero. A majority of solid red line tallies with black-dotted line; and the observed value matches the actual value well. As shown in Figures 6 and 7, when initial observer value is different, values of angular displacement and speed observed from two observers are in line with actual values.Figure 4
Wheel-rail adhesion coefficient.
(a)
Original image (b)
Partially enlarged viewFigure 5
Wheel-rail adhesion coefficient derivative.
(a)
Original image (b)
Partially enlarged viewFigure 6
Motor angular displacement.
(a)
Original image (b)
Partially enlarged viewFigure 7
Motor speed.
(a)
Original image (b)
Partially enlarged viewCase 2.
Simulation environment for locomotive is always in a dry state, and initial value ofTm is zero, with an increase by 2000 N·m/s. Figures 8 and 9 show wheel-rail adhesion coefficient and adhesion coefficient derivative, respectively. In the figures, red solid, blue-dotted, and black-dotted lines represent actual value, observed value from load torque observer, and observed value from derivative observer, respectively. As shown in Figures 8 and 9, adhesion coefficient increases proportionally; adhesion coefficient derivative rises rapidly to a constant; and locomotive runs in creep section. Observed value is in line with actual value. As shown in Figures 10 and 11, observed values of motor angular displacement and speed from observers agree with actual values well, considering the varying initial observer value.Figure 8
Wheel-rail adhesion coefficient.
(a)
Original image (b)
Partially enlarged viewFigure 9
Wheel-rail adhesion coefficient derivative.
(a)
Original image (b)
Partially enlarged viewFigure 10
Motor angular displacement.
(a)
Original image (b)
Partially enlarged viewFigure 11
Motor speed.
(a)
Original image (b)
Partially enlarged view


### 4.2. Experimental Results and Analysis

RT-LAB is a powerful, modular, distributed, real-time platform from OPAL-RT Technologies, Inc. It supports model-based design using rapid control prototyping (RCP) and hardware-in-the-loop simulation (HILS) of complex dynamic systems. HILS differs from pure real-time simulation and RCP by the use of the real controller in the loop; this controller is connected to the rest of the system simulated by RT-LAB. Therefore, in HILS, the plant is simulated, but the controller is real [22, 23]. To implement the proposed approach, experiments have been carried out on an OP5600 RT-LAB platform using HILS. The RT-LAB platform is shown in Figure 12, and the configuration is shown in Figure 13. The observer is a TMS320F2812 digital signal processor, which implements high-performance control and computation. The motor and wheel-rail system are simulated by RT-LAB. The sampling period is chosen as 50 μs, which is the same as the SMO. The semiphysical simulation results are shown in Figures 14–17.Figure 12
RT-LAB experiment setup.Figure 13
The structure of RT-LAB hardware-in-the-loop system.Figure 14
Wheel-rail adhesion coefficient. Adhesion coefficient: 0.01/grid; time: 0.5 s/grid.Figure 15
Wheel-rail adhesion coefficient derivative. Adhesion coefficient derivative: 0.1/grid; time: 0.5 s/grid.Figure 16
Wheel-rail adhesion coefficient. Adhesion coefficient: 0.01/grid; time: 0.5 s/grid.Figure 17
Wheel-rail adhesion coefficient derivative. Adhesion coefficient derivative: 0.005/grid; time: 0.5 s/grid.Case 1.
Experiment environment for locomotive starts with dry track surface; at 2.5 s, locomotive travels from dry to wet rail surface.Case 2.
Experiment environment for locomotive is always in a dry state, and initial value ofTm is zero, with an increase by 2000 N·m/s.Results Analysis. Figures 14 and 15 are the experimental waveforms under the change of the rail surface. In Figure 14, the oscilloscope shows the adhesion coefficient; in Figure 15, the oscilloscope shows the adhesion coefficient derivative, CH3 channel shows the actual value, and CH1 channel shows the observed value. The rail surface changes from dry to wet state in 2.5 s, the locomotive running into the skid region. Figures 16 and 17 are experimental waveforms under the change of traction force. In Figure 16, the oscilloscope shows the adhesion coefficient; in Figure 17, the oscilloscope shows the adhesion coefficient derivative, CH3 channel shows the actual value, and CH2 channel shows the observed value. It can be seen that the adhesion coefficient is increased by the positive proportion, and the locomotive runs in the microslip area. The experimental results are in agreement with the simulation results.


## 4.1. Simulation Results and Analysis

The adhesion coefficient and its derivative are testified by a uniaxial locomotive simulation model established using MATLAB/Simulink.Case 1.
Simulation environment for locomotive starts with dry track surface, andTm remains 6500 N·m. Figures 4 and 5 show wheel-rail adhesion coefficient and adhesion coefficient derivative, respectively. In the figures, solid red and black-dotted lines represent actual and observed values, respectively. Figures 6 and 7 show motor angular displacement and speed, respectively. In the figures, solid red, blue-dotted, and black-dotted lines represent actual value, observed value from load torque observer, and observed value from derivative observer, respectively. As shown in Figures 4 and 5, adhesion coefficient immediately increases from zero to 0.11, and adhesion coefficient derivative is zero. At 2.5 s, locomotive travels from dry to wet rail surface, and wheel-rail adhesion does not meet locomotive traction needs. When the locomotive runs to idling section, adhesion coefficient declines to zero. A majority of solid red line tallies with black-dotted line; and the observed value matches the actual value well. As shown in Figures 6 and 7, when initial observer value is different, values of angular displacement and speed observed from two observers are in line with actual values.Figure 4
Wheel-rail adhesion coefficient.
(a)
Original image (b)
Partially enlarged viewFigure 5
Wheel-rail adhesion coefficient derivative.
(a)
Original image (b)
Partially enlarged viewFigure 6
Motor angular displacement.
(a)
Original image (b)
Partially enlarged viewFigure 7
Motor speed.
(a)
Original image (b)
Partially enlarged viewCase 2.
Simulation environment for locomotive is always in a dry state, and initial value ofTm is zero, with an increase by 2000 N·m/s. Figures 8 and 9 show wheel-rail adhesion coefficient and adhesion coefficient derivative, respectively. In the figures, red solid, blue-dotted, and black-dotted lines represent actual value, observed value from load torque observer, and observed value from derivative observer, respectively. As shown in Figures 8 and 9, adhesion coefficient increases proportionally; adhesion coefficient derivative rises rapidly to a constant; and locomotive runs in creep section. Observed value is in line with actual value. As shown in Figures 10 and 11, observed values of motor angular displacement and speed from observers agree with actual values well, considering the varying initial observer value.Figure 8
Wheel-rail adhesion coefficient.
(a)
Original image (b)
Partially enlarged viewFigure 9
Wheel-rail adhesion coefficient derivative.
(a)
Original image (b)
Partially enlarged viewFigure 10
Motor angular displacement.
(a)
Original image (b)
Partially enlarged viewFigure 11
Motor speed.
(a)
Original image (b)
Partially enlarged view


## 4.2. Experimental Results and Analysis

RT-LAB is a powerful, modular, distributed, real-time platform from OPAL-RT Technologies, Inc. It supports model-based design using rapid control prototyping (RCP) and hardware-in-the-loop simulation (HILS) of complex dynamic systems. HILS differs from pure real-time simulation and RCP by the use of the real controller in the loop; this controller is connected to the rest of the system simulated by RT-LAB. Therefore, in HILS, the plant is simulated, but the controller is real [22, 23]. To implement the proposed approach, experiments have been carried out on an OP5600 RT-LAB platform using HILS. The RT-LAB platform is shown in Figure 12, and the configuration is shown in Figure 13. The observer is a TMS320F2812 digital signal processor, which implements high-performance control and computation. The motor and wheel-rail system are simulated by RT-LAB. The sampling period is chosen as 50 μs, which is the same as the SMO. The semiphysical simulation results are shown in Figures 14–17.Figure 12
RT-LAB experiment setup.Figure 13
The structure of RT-LAB hardware-in-the-loop system.Figure 14
Wheel-rail adhesion coefficient. Adhesion coefficient: 0.01/grid; time: 0.5 s/grid.Figure 15
Wheel-rail adhesion coefficient derivative. Adhesion coefficient derivative: 0.1/grid; time: 0.5 s/grid.Figure 16
Wheel-rail adhesion coefficient. Adhesion coefficient: 0.01/grid; time: 0.5 s/grid.Figure 17
Wheel-rail adhesion coefficient derivative. Adhesion coefficient derivative: 0.005/grid; time: 0.5 s/grid.Case 1.
Experiment environment for locomotive starts with dry track surface; at 2.5 s, locomotive travels from dry to wet rail surface.Case 2.
Experiment environment for locomotive is always in a dry state, and initial value ofTm is zero, with an increase by 2000 N·m/s.Results Analysis. Figures 14 and 15 are the experimental waveforms under the change of the rail surface. In Figure 14, the oscilloscope shows the adhesion coefficient; in Figure 15, the oscilloscope shows the adhesion coefficient derivative, CH3 channel shows the actual value, and CH1 channel shows the observed value. The rail surface changes from dry to wet state in 2.5 s, the locomotive running into the skid region. Figures 16 and 17 are experimental waveforms under the change of traction force. In Figure 16, the oscilloscope shows the adhesion coefficient; in Figure 17, the oscilloscope shows the adhesion coefficient derivative, CH3 channel shows the actual value, and CH2 channel shows the observed value. It can be seen that the adhesion coefficient is increased by the positive proportion, and the locomotive runs in the microslip area. The experimental results are in agreement with the simulation results.


## 5. Conclusions

Obtaining the wheel-rail adhesion state is extremely hard. This study designs a method to estimate the adhesion coefficient and its derivative online based on the cascading SMC observer. Mathematical model of wheel-rail is obtained by analyzing wheel-rail adhesion mechanism, which is proposed to estimate adhesion coefficient by load torque. The derivative of the torque load in the existing algorithms is required to be 0; that is,T˙L=0. However, in the practical operation of locomotives, this is difficult to guarantee because of the changing track conditions and creep speed. In the present study, this assumption is relaxed; that is, it should only be assumed as a bounded quantity. The derivative of the load torque is estimated using the principle of equivalence in SMC. The chattering problem of the SMC observer is then solved. The algorithm proposed in this study has a simple structure with high robustness and is easy to implement. By means of simulation and experiment, the observed values are consistent with the actual values, which confirms the algorithm validity



---

*Source: 8419295-2017-01-03.xml*