METHOD FOR CALCULATING THE ADVANCE SPEED OF A RAILWAY VEHICLE

Abstract

Disclosed is a method for calculating or estimating the speed of a railway vehicle. The method comprises the steps of generating speed signals indicating the angular speed () of the wheels of said at least one controlled axle; estimating, as a function of such angular speed (), the value of the adhesion () in the contact area of the wheels of such axle and the rails, using an adhesion observer and computing the value of the speed slip () of the wheels of such controlled axle, generating signals representative of the derivative

[00001]$\left(\frac{d.Math..Math.}{d.Math..Math.}\right)$

of said adhesion () as a function of the slip (); generating a driving signal (C(T.sub.j+1)) for torque control devices controlling the torque applied to the wheels applying the driving signal (C(T.sub.j+1)) to said torque control devices and therefore computing the vehicle speed as the linear advance speed of such at least one controlled axle

Claims

1. A method for calculating or estimating the speed of a railway vehicle, at least one axle of which has an associated control system of the adhesion of the wheels to the rails; the method comprising the steps of: generating speed signals indicating the angular speed () of the wheels of said at least one axle; estimating, as a function of said angular speed (), the value of the adhesion () in the contact area of the wheels of said axle and the rails, using an adhesion observer and calculating the value of the speed slip () of the wheels of said controlled axle, generating signals representative of the derivative $\left(\frac{d.Math..Math.}{d.Math..Math.}\right)$ of said adhesion () as a function of the slip () of the wheels of said axle; generating a driving signal (C(T.sub.j+1)) for torque control means controlling the torque applied to the wheels of said axle, by means of an adaptive control of said derivative signals $\left(\frac{d.Math..Math.}{d.Math..Math.}\right)$ as a function of an error signal (e((T.sub.j+1)) indicative of the difference between the value of said derivative $\left(\frac{d.Math..Math.}{d.Math..Math.}\right)$ and a predetermined reference value such as to reduce and keep said difference substantially at zero; and applying said driving signal (C(T.sub.j+1)) to said torque control means and then computing the vehicle speed as the linear advance speed of said at least one controlled axle.

2. A method according to claim 1, wherein said driving signal (C(Tj+1)) is generated by means of an adaptive filtering of LMS type.

3. A method according to claim 1, wherein said driving signal (C(Tj+1)) is generated by the time integration of the derivative $\left(\frac{d.Math..Math.}{d.Math..Math.}\right)$ of the adhesion () as a function of the slip ().

4. A method according to claim 1, wherein said driving signal (C(Tj+1)) is generated by time integration of the sign of the derivative $\left(\frac{d.Math..Math.}{d.Math..Math.}\right)$ of the adhesion () as a function of the slip ()

5. A method according to claim 1, wherein said reference value $\left(\frac{d.Math..Math.}{d.Math..Math.}\right)$ is equal or greater than zero.

6. A method for computing or estimating the speed of a railway vehicle of claim 1 in a braking condition.

7. A method for computing or estimating the speed of a railway vehicle of claim 1 in a traction condition.

Description

[0009] Further features and advantages of the invention will become apparent from the following detailed description, provided purely by way of non-limiting example with reference to the accompanying drawings, in which:

[0010] FIG. 1 is a block diagram of an anti-skid control system of the wheels of a railway vehicle;

[0011] FIG. 2 is a block diagram of a closed loop control system of an axle's rotation speed;

[0012] FIG. 3 is a graph showing qualitatively the trend of the adhesion coefficient of the wheels of an axle, shown on the y-axis, as a function of the slip , shown on the x-axis;

[0013] FIG. 4 is a diagram illustrating the forces applied to an axle's wheel;

[0014] FIGS. 5A and 5B are graphs used to illustrating in detail the control criterion to which the present invention refers, FIG. 5B showing in enlarged scale a part of the graph of FIG. 5A;

[0015] FIGS. 6, 7 and 8 are block diagrams relating to different systems for the implementation of the method according to the present invention, and

[0016] FIG. 9A, already described, is a train with two independent cars, and FIG. 9B, likewise already described, is a train with two cars bound with a Jacobs bogie.

[0017] Electronic systems are installed on board most modern rail vehicles, which typically include wheel slide control subsystems, intended to intervene both when the vehicle is in the traction phase and when it is in the braking phase. These subsystems are known as anti-skid or anti-slide systems, or also WSP (Wheel Slide Protection) systems.

[0018] A system for controlling the adhesion of the wheels, as an anti-skid function, according to the prior art, is schematically represented in FIG. 1 of the accompanying drawings, which refers to a vehicle with n controlled axles A1, A2, . . . , An. The axles A1, A2, . . . , An comprise a respective shaft S1, S2, . . . , Sn and a respective wheelset W1, W2, . . . , Wn integral in rotation to it.

[0019] In the drawings, only one wheel of each axle is generally illustrated.

[0020] The WSP system of FIG. 1 comprises an electronic control unit ECU, typically based on microprocessor architecture, that receives tachometer signals relating to the angular speed of each axle A1, A2, . . . An from sensors SS1, SS2, . . . , SSn respectively associated to these axles. The electronic control unit ECU is also connected to the torque control apparatuses TC1, TC2, . . . , TCn, each associated to a respective axle A1, A2, . . . , An.

[0021] The electronic control unit ECU is provided to operate a modulation of the torque applied to each axle according to a predetermined algorithm if, in the case of applying torque during traction or braking in a degraded adhesion condition, the wheels of one or more axles end up in a possible incipient skidding condition. Torque modulation is implemented in such a way as to prevent a total locking of the axles, possibly so as to bring each axle into a condition of controlled slipping with the intention of recovering adhesion and in any case for the entire duration of the degraded adhesion condition.

[0022] FIG. 2 shows a block diagram illustrating an adhesion control/recovery system for a generic axle: the error or difference E(t) between the reference speed value V.sub.R(t) at which one wishes to slide the controlled axle A and the measured speed V.sub.M(t) detected by the associated sensor SS and conditioned by an acquisition and processing module APM is applied as an input signal to a control module CM, which outputs a drive signal Y(t) to the torque control apparatus TC associated with the axle A.

[0023] The reference speed V.sub.R(t) is obtained as a fraction of the instantaneous speed of the vehicle, for example, according to the expression:

V.sub.s(t)=V.sub.s(t)(1) (1)

[0024] where V.sub.V(t) is the instantaneous (computed) speed of the vehicle, and represents the related slip of the axle A to be obtained during the skidding phase.

[0025] It is evident how the knowledge of the vehicle's instantaneous speed V.sub.V(t) is essential for properly controlling skidding.

[0026] In the event of braking, the algorithm most used for the estimation of the vehicle's actual speed V.sub.V(t) normally uses a function of the type:

V.sub.v(T.sub.j)=max [S.sub.1(T.sub.j), . . . , S.sub.n(T.sub.j), (V.sub.v(T.sub.j1)+a.sub.max.Math.T)](2)

[0027] while in the event of traction, the following function is used:

V.sub.v(T.sub.j)=min[S.sub.1(T.sub.j), . . . , S.sub.n(T.sub.j), (V.sub.v(T.sub.j1)+a.sub.max.Math.T) (3)

[0028] where a.sub.max is the maximum acceleration permitted for the vehicle in operation, this acceleration having a positive sign in the case of a traction condition and a negative sign in the case of a braking condition.

[0029] The contribution (V.sub.v(T.sub.j1)+a.sub.max.Math.T) in the relationships (2), (3) serves to contain the variation of the V.sub.v(t) within physical limits allowed by the train, when excessive instantaneous and concurrent variations of the speeds of the axles due to particularly degraded adhesion conditions, in traction or braking conditions, could lead to a loss of significance of the speed V.sub.v(t) computed with these relationships (2), (3).

[0030] More accurate variants of the relationships (2), (3) are known but still based on the instantaneous measurement of the individual speed of the axles. It becomes evident here how the availability of a dead axle would make (2), (3) extremely accurate if all the axles were subjected to torque during skidding phase.

[0031] By way of non-limiting example, one possible embodiment of a torque control apparatus TC of the torque applied to an axle is described and illustrated in the previous Italian patent application No. 102015000086465 filed Dec. 22, 2015.

[0032] A torque control apparatus may be made according to many variants known to persons skilled in the art.

[0033] As is well known, the adhesion coefficient () between wheels and rails varies according to the slip substantially in the way illustrated in FIG. 3. Based on the expression (1) above, may be expressed as with 0V.sub.rV.sub.v and 01.

[00002]$\begin{array}{cc}=\frac{\mathrm{Vv}-\mathrm{Vr}}{\mathrm{Vv}}& \left(4\right)\end{array}$

[0034] In FIG. 3, the curves 1, 2 and 3 qualitatively represent the trend of the adhesion according to the environmental conditions: curve 1 corresponds to an adhesion condition in dry contact conditions between the wheels and rails, curve 2 corresponds to an adhesion condition in the presence of moisture between the wheels and rails, and curve 3 represents an adhesion condition in the presence of viscous material between the wheels and rails, such as oil or rotten leaves (typical condition in the autumn period), or even rust mixed with moisture (typical condition in railway depots).

[0035] It has been found experimentally that the values of at the adhesion peaks a.sub.1, a.sub.2, a.sub.3 vary with the change in the adhesion conditions, which move along a curve as indicated at A in FIG. 3.

[0036] Experimental measurements demonstrate how the curve A lies in an area corresponding to values 00.02 even in very degraded adhesion conditions.

[0037] If one or more axles, for example the one previously defined as the dead axle, can be maintained on the curve A during traction or braking, it is achieved the dual effect of using, for said axles, the maximum available adhesion and at the same time tracking the actual speed of the train, corresponding to =0, with a maximum error of 2%.

[0038] FIG. 4 is a diagram illustrating forces applied to an axle's wheel A. From this figure, it is clear that:

F.sub.m.Math.R=F.sub.A.Math.RJ.Math.{dot over ()}(5)

[0039] where:

F.sub.A=.Math.m.Math.g (6)

[0040] for which:

F.sub.m=.Math.m.Math.gJ/R.Math.{dot over ()}(7)

[0041] where F.sub.m is the tangential force applied to a wheel by the traction and/or braking system, R is the radius of the wheel, J is the moment of inertia of the axle, m is the mass applied to the wheel-rail contact area, and {dot over ()} is the instantaneous angular acceleration of the axle.

[0042] It is clear that at the same instantaneous angular acceleration, the maximum applicable force F.sub.m is obtained in correspondence with the maximum value of adhesion , i.e. at the points lying on the curve A of FIG. 3.

[0043] The method according to the present invention uses an adhesion observer to evaluate in real time the adhesion value at the contact area between the wheels and rails for one or more axles during a skidding phase and, by processing these values in real time, identifies continuously over time the value to be assigned to a slip control system.

[0044] An adhesion observer adapted to dynamically identify the instantaneous value (T.sub.j) of the adhesion in a generic sampling period T.sub.j of a predetermined duration T at the wheel-rail contact area during skidding is definable using the equation (7) from which with some simple steps the following relationship is obtained:

[0045] where

[00003]$\begin{array}{cc}\left(T_j\right)=\frac{1}{m.Math.g}.Math.\left[F_m\left(T_j\right)+J/R.Math.\stackrel{.}{}\left(T_j\right)\right]& \left(8\right)\end{array}$ [0046] {dot over ()} is the angular acceleration of the axle, i.e. the time derivative of the angular speed co of the axle; the value of this acceleration is already available in real time within a control and adhesion recovery system, because angular acceleration is one of the variables on which the control function implemented by the block CM of FIG. 2 is normally based for controlling the slip of the axle; the sign of {dot over ()} depends on the instantaneous acceleration or deceleration condition of the axle; [0047] m is the mass on the wheel-rail contact area; in the latest generation trains, the m value is known in real time, as it is commonly available to the system that computes the accelerating/braking force to apply to the axle to obtain the desired accelerations/decelerations; [0048] J is the moment of inertia of the axle and is a parameter whose value is always known, being supplied by the manufacturer of the carriages, as it is fundamental for the computation of stopping distances; [0049] F.sub.m, already defined above in relation to FIG. 4, can be obtained by multiplying the pressure applied to the brake cylinder, known to the braking system, for pressure/force conversion coefficients typical of the brake cylinder, as well as the transmission and efficiency coefficients of the levers and of the friction coefficient between the brake linings and discs (in the case of disc brakes); in the case of electrodynamic type traction or braking, the value of the force F.sub.m may be obtained from the electric current value supplied/regenerated by the motor in traction or, respectively, in braking; in the case of so-called blended braking, the intensity of the force F.sub.m may be determined as the sum of the respective contributions of the pneumatic brake and of the electrodynamic brake, appropriately weighed with respective coefficients; and [0050] T.sub.j is the generic j-th sampling period of the system with which the adhesion observer and more generally the method according to the invention is carried out; in the description that follows, T.sub.j will replace the use of the variable t representing time.

[0051] Downstream of the adhesion observer, a low-pass type filter may appropriately be provided, to remove or at least mitigate instantaneous and noise variations present outside of the frequency band useful for a correct observation of the adhesion values.

[0052] An embodiment of a system for implementing a method according to the present invention is illustrated in FIG. 6.

[0053] The method provides for identifying and tracking the slip value of at least one axle, such that the curve () illustrated in FIG. 5 shows the maximum value, i.e. the value for which

[00004]$\frac{d.Math..Math.\left(T\right)}{d.Math..Math.\left(T\right)}=0.$

[0054] For this purpose, a system implementing an LMS algorithm (Least Mean Square) may be used. For an accurate description of the general features of the convergence criteria and the implementation variants of LMS algorithms, please refer to the available literature and in particular to the text: B. Widrow, S. D. Stearns, Adaptive Signal Processing, New Jersey, Prentice-Hall, Inc., 1985.

[0055] With reference to FIG. 6, an adhesion observer 1201 receives input signals representative of the value of the speed of the wheel of the controlled axle An that is to be maintained on the adhesion peak, together with an X vector containing the values of the magnitudes m(T.sub.j), J, R and F.sub.m(T.sub.j) previously described, for the estimation of the instantaneous value of adhesion (T.sub.j) relating to the controlled axle.

[0056] The output of the adhesion observer 1201 is connected to the input of a module 1202 which computes the value of the derivative

[00005]$\frac{d.Math..Math.}{d.Math..Math.},$

e.g. according to the equation:

[00006]$\begin{array}{cc}\frac{d.Math..Math.\left(T_j\right)}{d.Math..Math.\left(T_j\right)}=\frac{\left(T_j\right)-\left(T_{j-1}\right)}{\left(T_j\right)-\left(T_{j-1}\right)}& \left(9\right)\end{array}$

[0057] where the value of is obtained in real time in accordance with the equation (4).

[0058] An adder 1203 outputs an error signal e(T.sub.j) as the difference between the desired value of said derivative (i.e. the value 0) and its instantaneous value computed by the module 1202, and such error is used to adapt the LMS algorithm implemented in a block 1204.

[0059] The latter provides in output a torque request C(T.sub.j+1) for said axle, which is transmitted to a torque control module 1205 of a per se known type, having for example the architecture described in the previous aforementioned Italian patent application with reference to FIG. 3.

[0060] In a manner known per se, the module 1204 continuously corrects the output C(T.sub.j+1) in order to minimize or nullify the error e(T), i.e. in order to obtain a nullification of the aforementioned derivative, i.e. in order to bring and maintain said controlled axle to the adhesion peak value.

[0061] By applying, therefore, the solution according to FIG. 6 to at least one axle, said axle will always advance at a linear speed equal to that of the vehicle (less than a maximum error that can be estimated within 2%), even in degraded adhesion conditions, at the same time providing the maximum force value, in traction or braking, made possible by the available adhesion.

[0062] A simplified implementation of the group of modules included in the dashed line block 1206 of FIG. 6 is illustrated in FIG. 7, where the block 1204, which implements the LMS algorithm, is replaced with a simple integrator 805, the output of which, amplified with a gain K, generates the torque value C(T.sub.j+1) to be assigned to the adhesion control and recovery system 1205. In such case, when

[00007]$\frac{d.Math..Math.}{d.Math..Math.}>0$

the integrator 805 increases the torque value C(T.sub.j+1), when

[00008]$\frac{d.Math..Math.}{d.Math..Math.}<0$

the integrator 805 decreases the torque value C(T.sub.j+1), and when

[00009]$\frac{d.Math..Math.}{d.Math..Math.}=0$

the integrator 805 keeps the torque value C(T.sub.j+1) stable.

[0063] In this way, the system brings and maintains said controlled axle to the peak adhesion value.

[0064] The gain K regulates the identification speed of the average adhesion peak value and simultaneously ensures the stability of the closed loop system.

[0065] A further simplified variant of embodiment of the dashed block 1206 of FIG. 6 is shown in FIG. 8: the module 903 determines the sign of the derivative

[00010]$\frac{d.Math..Math.}{d.Math..Math.}.$

The output of the block 903 being equal to +1 or 1 (positive and, respectively, negative sign), a subsequent integrator 805 performs simple unitary sums.

[0066] The integrator 805 may be replaced with an up/down type counter updated with period T=T.sub.j+1T.sub.j.

[0067] The diagrams according to FIGS. 7 and 8 perform a continuous tracking of the average adhesion peak , continuously adapting to the change in adhesion conditions, similarly to what was achieved with the diagram according to FIG. 6. The latter allows rapid and accurate tracking of the condition

[00011]$\frac{d.Math..Math.}{d.Math..Math.}=0,$

but requires me use or a certain number of computations in real time.

[0068] The diagram according to FIG. 8 greatly reduces the number of computations necessary, but also reduces the tracking speed of the condition

[00012]$\frac{d.Math..Math.\mathrm{.Math.}}{d.Math..Math.}=0.$

[0069] The diagram according to FIG. 7 has features intermediate between those of the diagrams according to FIGS. 6 and 8.

[0070] Therefore, the two expressions (2), (3) provided above always allow a very reliable value of the vehicle's speed V.sub.v to be provided, even in very degraded adhesion conditions.

[0071] If it is desired to further increase the accuracy of the tracking of the train's speed, it is sufficient to compute the error with respect to values of

[00013]$\frac{d.Math..Math.\mathrm{.Math.}}{d.Math..Math.}>0,$

i.e. on the left side of the curve illustrated in FIG. 5B, at the expense of the applied torque value, which will prove to be lower than the maximum peak as a function of the increase of the applied reference value

[00014]$\frac{d.Math..Math.\mathrm{.Math.}}{d.Math..Math.}.$

[0072] Naturally, without altering the principle of the invention, the embodiments and the details of implementation may vary widely with respect to those described and illustrated purely by way of non-limiting example, without thereby departing from the scope of the invention as defined in the appended claims.