Method of Estimating the Instantaneous Engine Speed Produced by Each Cylinder of an Internal-Combustion Engine

Abstract
The invention is a method for real-time estimation of the instantaneous engine speed produced by each cylinder of an internal-combustion engine, from an instantaneous engine speed measurement at the end of the engine transmission system. A physical model, representing in real time the dynamics of the transmission system according to the crankshaft angle and to coefficients of a Fourier series decomposition of the instantaneous speed produced by each cylinder, is constructed. These coefficients are determined in real time from coupling between the model and an adaptive type non-linear estimator. The instantaneous speed produced by each cylinder is then deduced from these coefficients. The mean torque produced by each cylinder can also be deduced therefrom. An application is: engine controls.
Description
BACKGROUND OF THE INVENTION

1. Field of the Invention


The present invention relates to a method intended for real-time estimation of the instantaneous engine speed produced by each cylinder of an internal-combustion engine from the instantaneous speed detector located at the end of the transmission system.


2. Description of the Prior Art


Knowledge of the instantaneous speed for each cylinder allows estimation of the mean torque produced by each cylinder.


Estimation of the mean torque produced by each cylinder is important for all vehicles, whether equipped with gasoline or diesel engines. In the first case, it conditions good combustion of the mixture when the fuel/air ratio is close to 1, and therefore sensitive to cylinder to cylinder difference problems. In the second case, the knowing the torque allows readjustment so as to obtain optimum running conditions. Catalysts using a NOx trap lose efficiency in the course of time. In order to recover optimum efficiency, the torque of each cylinder has to be kept identical for some seconds, prior to returning to normal running conditions with a lean mixture. Removing pollution with DeNox catalysis therefore requires precise control of the torque cylinder by cylinder.


An instantaneous engine speed detector is therefore arranged at the end of the transmission system. This measurement is greatly distorted by the transmission and it affected by noise.


In order to control more precisely, and in particular individually, injection of the fuel masses into the cylinders, reconstruction of the torque cylinder to cylinder is necessary. Installing a digital torquemeter below each cylinder of a vehicle cannot be done considering the cost price thereof.


The method according to the invention provides an estimator, working from the measurement performed at the end of the transmission chain, to estimate the instantaneous engine speed below each cylinder.


SUMMARY OF THE INVENTION

The invention relates to a method for real-time estimation of the instantaneous engine speed produced by each cylinder of an internal-combustion engine comprising at least one transmission system connected to the cylinders and a detector performing real-time measurement (x1) of the instantaneous engine speed at the end of said transmission system.


The method comprises:


a) constructing a physical model representing in real time the dynamics of the transmission system according to: the measurement (x1), coefficients of a Fourier series representing decomposition of the instantaneous engine speed produced by each cylinder, and according to a damping and to a natural frequency of the transmission system;


b) determining, in real time, the coefficients of the Fourier series representing decomposition by coupling the model with an adaptive type non-linear estimator; and


c) carrying out real-time estimation of the instantaneous engine speed produced by each cylinder from the Fourier coefficients.


The mean torque of each cylinder can also be estimated in real time from the estimation of these coefficients.


The method according to the invention can be applied to an engine control to control the fuel masses injected into each cylinder so as to adjust the mean torque produced by each cylinder.





BRIEF DESCRIPTION OF THE DRAWINGS

Other features and advantages of the method according to the invention will be clear from reading the description hereafter of embodiments given by way of non limitative example, with reference to the accompanying figures wherein:



FIG. 1 illustrates the estimation of the instantaneous engine speed below the cylinders by means of the method according to the invention, on a working point of 1250 rpm at medium load; and



FIG. 2 illustrates the estimation of the mean torque cylinder to cylinder by means of the method according to the invention, on a working point of 1500 rpm.





DETAILED DESCRIPTION

The method according to the invention allows estimation of the instantaneous engine speed produced by each cylinder of an internal-combustion engine comprising at least one transmission system connected to the cylinders. At the end of this transmission system, a detector performs real-time measurement of the instantaneous engine speed. This signal is denoted by x1. Measurement of the instantaneous engine speed below the cylinders, distorted by the drive shaft, is thus performed. The first stage of the invention thus is “reversing” the effects of the transmission to obtain the relevant information, that is the instantaneous engine speed produced by each cylinder. This relevant information is a periodic signal denoted by x0.


The method mainly comprises:


1—establishing, in an angular scale (that depending on the crankshaft angle and not on time), a physical model representing in real time the dynamics of the transmission system;


2—describing the instantaneous engine speed produced by each cylinder by quasi time-invariant parameters such as the coefficients of the Fourier analysis of the instantaneous engine speed;


3—coupling the physical model with an adaptive type non-linear estimator; and


4—carrying out real-time estimation of the instantaneous engine speed produced by each cylinder from the adaptive type non-linear estimator.


1—Physical Model of the Transmission System Dynamics


To estimate signal x0, that is the instantaneous engine speed below the cylinders, a physical model of the transmission system dynamics is first defined. One therefore considers that this system behaves like a second-order system made up of two parameters:



ω: the natural frequency of the transmission in the rotating reference frame

ζ: transmission damping.


Thus, considering the angular scale, the dynamics of the drive shaft is written as follows:








{








2



(


x
1

-

x
0


)





α
2



=



-

ξ
_




ω
_






(


x
1

-

x
0


)




α



-



ω
_

2



(


x
1

-

x
0


)









y
=

x
1










with:


x1: instantaneous engine speed at the end of the transmission chain: the measurement


x0: instantaneous engine speed below the cylinders which is the unknown



ω: natural frequency of the transmission system in the rotating reference frame



ζ: damping of the transmission system


α: crankshaft angle of the transmission system.


A variable change can be performed by putting:










w
0

=





2



x
0









α
2



+


ξ
_



ω
_






x
0




α



+



ω
_

2



x
0







(
1
)







The instantaneous engine speed below the cylinders x0 is periodic, therefore w0 is also periodic. The dynamics can therefore be rewritten in the form as follows:









{







x



α


=


A
·
x

+


A
0

·

w
0









y
=

C
·
x









(
2
)







with:






x
=

[




x
1







x
1




α





]







A
=

[



0


1





-


ω
_

2






-

ξ
_




ω
_





]








A
0

=

[



0




1



]







C
=

[



1


0



]





This equation (2) is the physical model representing in real time the transmission system dynamics. An estimation of signal w0 allows determination of an estimation of signal x0 from equation (1).


2—Description of Signal x0 by Quasi Time-Invariant Parameters


It is attempted to estimate, from this physical model and from measurement y (equal to x1), signal x0, that is the instantaneous engine speed produced by each cylinder. To perform this real-time estimation, the method according to the invention describes this signal x0 with quasi time-invariant parameters. In other words, signal x0 is defined by means of parameters which, at a given time, are constants. Therefore the fact is exploited that signal x0 is mechanically periodic. Thus, instead of performing a highly variable signal estimation x0, the Fourier coefficients of this signal can be estimated. It is also possible to use any parameter allowing description of signal x0 in connection with the periodic character thereof. The Fourier coefficient analysis of signal x0, developed into complex numbers for clarity reasons, is written as follows:











x
0



(
α
)


=




j
=

-
n


n




d
j





(







ωα

)








(
3
)







The dj represent the 2n+1 Fourier coefficients of the decomposition of signal x0.


We thus define a signal expressing the instantaneous engine speed x0 according to the time-invariant parameters dj.


To estimate parameters dj, it is possible to use again variable change w0 and to use the physical model described by system (2). Signal w0 is also mechanically periodic, and its Fourier coefficient analysis, developed into complex numbers for clarity reasons, is written as follows:








w
0



(
α
)


=




j
=

-
n


n




c
j





(







ωα

)








The cj represent the 2n+1 Fourier coefficients.


Estimation of these coefficients cj thus allows estimation of the Fourier coefficient decomposition of signal x0 and therefore signal x0 itself.


Using only a finite number of harmonics ([−n; +n]), the physical model representing in real time the transmission system dynamics is then written as follows:









{











x



α


=


A
·
x

+


A
0

·

(




j
=

-
n


n




c
j





(







ω





α

)




)













c
j




α


=
0










y
=

C
·
x


,












j


[


-
n

,
n

]








(
4
)







3—Coupling with an Adaptive Type Non-Linear Estimator


From the physical model described by system (4), an adaptive type non-linear estimator is defined comprising, on the one hand, a term linked with the dynamics and, on the other hand, a correction term:









{









x
^




α


=


A
·

x
^


+


A
0

·




j
=

-
n


n





c
^

j

·



(







ω





α

)





-

L
·

(


C
·

x
^


-
y

)















c
^

j




α


=


-



(


-







ω





α

)



·

L
j

·

(


C
·

x
^


-
y

)



,












j


[


-
n

,
n

]








(
5
)







with:


{circumflex over (x)}: estimator of x


ĉj: estimator of cj

L: a matrix to be calibrated


Lj: matrices to be calibrated.


A selection of matrices L and Lj providing convergence of the estimator is:






L
=


[




2






ξ
_



ω
_







2







ω
_

2





]






and








j


[


-
n

,
n

]











L
j

=

1


j
2

+
1






The system of equations (5) represents an adaptive type non-linear estimator allowing estimation of coefficients cj of the Fourier coefficient analysis of the signal w0.


This estimator (5) is constructed from variable change w0, but it is clear that it is possible to construct in the same manner an adaptive type non-linear estimator directly from x0.


4—Real-Time Estimation of the Instantaneous Engine Speed Produced by Each Cylinder


It is then estimated, from estimation ĉj of coefficients cj, the instantaneous engine speed produced by each cylinder x0.


Estimator (5) allows reconstruction of w0 through its Fourier coefficients cj. The goal is to reconstruct x0. By means of the expression of w0 given by equation (1), coefficients dj are expressed as a function of coefficients cj:











d
j

=





ω
_

2

-


(

j
·
ω

)

2

-

i
·
j
·
ω
·

ξ
_

·

ω
_






(



ω
_

2

-


(

j
·
ω

)

2


)

2

+


(


j
·
ω
·

ξ
_




·

ω
_



)

2



·

c
j











j


[


-
n

,
n

]







(
6
)







thus the expression of the instantaneous engine speed produced by each cylinder, by means of equations (3) and (6), and the coefficients of its Fourier decomposition by means of equation (6) is obtained.


Estimation of the Mean Torque Produced by Each Cylinder

According to the invention, it is possible to provide an estimation of the mean torque produced by each cylinder from the estimation of the instantaneous engine speed produced by each cylinder (x0) and more precisely from the estimation of its Fourier analysis into coefficients dj.


Knowledge of the mean torque produced by each cylinder is fundamental and relevant information for combustion estimation; it is the image of the combustion that takes place in the engine.


The previous estimator (5) allows estimation of the signal of the engine speed below the cylinders as well as the Fourier analysis thereof. Now, the higher the torque, the higher the excitation on the shaft. It is thus possible to correlate the torque produced by the cylinder and the Fourier coefficients of the analysis of the instantaneous engine speed signal (x0).


In general terms, it is thus possible to identify a function φ that allows determination of the MIP (Mean Indicated Pressure) or, in an equivalent manner, the mean torque from coefficients dj:






ϕ


:










R


2





n

+
1







R





{

d
j

}






RMI







This function φ can be a polynomial function. It can be determined empirically from tests. The following function φ can be selected for example:










ϕ


(

d
j

)


=





j
=

-
n


,

j

0


n




d
j
2


ϕ
0







(
7
)







with φ0 being a constant to be calibrated according to the engine speed used, by means of correlations with engine test bench measurements. This calibration can be carried out from a tabulation obtained from a linear optimization consisting in adjusting the value of φ0 so that the estimations are as close as possible to the engine parameters (parameters allowing engine calibration and provided by the manufacturer).


Results


FIG. 1 illustrates the estimation (Rest) of the instantaneous engine speed x0 below the cylinders from the estimator according to the invention (5) described above on a working point of 1250 rpm at medium load. FIG. 1 also illustrates the reference instantaneous engine speed Rref (calculated from the cylinder pressure measurements on the engine test bench). A very good signal estimation is observed.



FIG. 2 illustrates the estimation (PMIest) of the torque cylinder to cylinder with a working point at 1500 rpm, from the estimator according to the invention (5) and a function φ defined by equation (7). FIG. 2 also illustrates the reference mean torque (PMIref) (calculated from the cylinder pressure measurements on the engine test bench). A very good signal estimation is observed.


The adaptive filter thus achieved is efficient and, in particular, it requires no additional adjustment in case of working point change. No identification stage is required, only a measurement noise and model adjustment has to be performed once.


An engine control can thus, from the reconstructed torques, adjust the fuel masses injected into each cylinder so that the torques are balanced in all the cylinders.


An estimation of the instantaneous engine speed produced by each cylinder and the estimation of the mean torque cylinder to cylinder have many advantages:

    • emissions reduction,
    • improved driveability (delivered torque regulation),
    • fuel consumption reduction,
    • injection system diagnosis (detection of the drift of an injection nozzle or of the failure of the injection system).

Claims
  • 1-3. (canceled)
  • 4. A method of real-time estimation of instantaneous engine speed produced by each cylinder of an internal-combustion engine including at least one transmission system connected to the cylinders and a detector performing real-time measurement of the instantaneous engine speed coupled to the transmission comprising: a) constructing a physical model, representing in real time, dynamics of the transmission system according to the crankshaft angle, the measurement, coefficients of a Fourier series analysis of the instantaneous engine speed produced by each cylinder, and a damping and a natural frequency of the transmission system;b) determining in real time the coefficients of the Fourier series analysis by coupling the model with an adaptive type non-linear estimator; andc) carrying out real-time estimation of the instantaneous engine speed produced by each cylinder from the Fourier coefficients.
  • 5. A method as claimed in claim 4, wherein mean torque of each cylinder is estimated in real time from an estimation of the coefficients.
  • 6. The method as claimed in claim 4 comprising using the real-time estimation of the instantaneous engine speed to provide engine control to control fuel masses injected into each cylinder in order to adjust mean torque produced by each cylinder.
  • 7. The method as claimed in claim 5 comprising using the real-time estimation of the instantaneous engine speed to provide engine control to control fuel masses injected into each cylinder in order to adjust mean torque produced by each cylinder.
Priority Claims (1)
Number Date Country Kind
0509624 Sep 2005 FR national
PCT Information
Filing Document Filing Date Country Kind 371c Date
PCT/FR2006/002127 9/18/2006 WO 00 9/2/2008