METHOD OF DETECTING AND CONTROLLING DETONATION OF AN INTERNAL COMBUSTION ENGINE

Abstract: A method of detecting and controlling detonation, including the steps of: generating (10) a vibration signal (D) proportional to the intensity of vibration on the crankcase; wide-band-filtering (11) the vibration signal (D) to generate a first intermediate signal; rectifying (12) the first intermediate signal to generate a second intermediate signal; integrating (13) the second intermediate signal to generate a first numeric value; calculating (14) a logarithm of the first numeric value to obtain a second numeric value (xi); calculating (15) a mean value ( mu i) as a function of the second numeric value (xi); calculating (16) a difference between the second numeric value (xi) and the mean value ( mu i) to obtain a third numeric value ( delta i); comparing (17) the third numeric value ( delta i) with a predetermined threshold value ( delta 0); determining (18) the presence of detonation in the event the third numeric value ( delta i) has a first predetermined relationship with the predetermined threshold value ( delta 0); forming (19, 21) a detonation index (Ii) indicating the behaviour of the engine in terms of detonation; and obtaining (22) a spark lead correction value (CAi) as a function of the detonation index (Ii).

Description: METHOD OF DETECTING AND CONTROLLING DETONATION OF AN INTERNAL COMBUSTION ENGINE TECHNICAL FIELD The present invention relates to a method of detecting and controlling detonation of an internal combustion engine. BACKGROUND ART As is known, in certain operating conditions, internal combustion engines are subject to high-intensity detonation, i.e. spontaneous combustion of the mixture, which occurs in many engine cycles, and results in impaired efficiency, dangerous overheating and reduced life of the engine, as well as possibly in sudden failure of certain engine components. Various systems are therefore being studied, or have already been proposed, to detect detonation and act accordingly on combustion parameters to reduce its effects and the likelihood of it occurring. Such systems are mainly based on directly or indirectly monitoring the pressure cycle pattern in the cylinder as a function of the position of the piston inside the cylinder. The pressure cycle pattern is bell- shaped, reaching a peak close to the top dead-center position, which peak is a classic rounded shape when combustion is normal (Figure la) , and has numerous indentations in the presence of detonation (Figure lb) . Information by which to detect detonation is obtained by analyzing the indentations. Some known methods employ sensors located inside the combustion chamber to directly determine the amplitude of the indentations. While providing for highly accurate, reliable amplitude values, such a sensor arrangement calls for sophisticated high-cost technology and is therefore only suitable for laboratory or prototype application. Other methods employ sensors located, and for detecting vibration, on the crankcase. Though technically simpler and cheaper, the values obtained by such sensors are subject to greater interference than those detected directly, on account of the vibration measured on the crankcase also being the result of other phenomena besides crankcase-filtered variations in pressure inside the cylinder. Patent GB-A-2 265 006 filed by Nippondenso Co Ltd, for example, relates to a detonation control system featuring detonation sensors on the crankcase, and which detects detonation by comparing the intensity of the sensor signal with a decision threshold, and provides for eliminating detonation by accordingly controlling engine operating parameters, such as injection advance. More specifically, the system performs a logarithmic conversion of the intensity of the sensor signal, determines its distribution, calculates a value corresponding to the standard deviation of the distribution, and compares the value with a threshold calculated on the basis of the previously calculated value and a mean value of said distribution. The distribution is determined by processing numeric values representing the amplitudes of spectral components obtained by narrow-band filtering the output signal of the sensors. The output signal, in fact, comprises numerous harmonics, and narrow-band filtration is employed to only extract the highest-energy harmonic. It is important to note that, in the above patent, the detonation threshold is adapted at each cycle. In some cases, however, controlling operation of the engine on the basis of one frequency in the spectrum and one engine cycle at a time may prove restrictive. That is, considering only one frequency may result in others, different from the one considered but generated all the same by detonation, being overlooked, while control designed solely to eliminate local detonation cycle by cycle may not always be the best solution, in that, in terms of engine operation and efficiency, it sometimes pays to tolerate a given number of detonations every so many engine cycles and/or in only some of the cylinders. Trace detonation, in fact, by enabling high efficiency of the engine, is a favourable engine operating condition, and therefore one to be encouraged. Finally, adapting the decision threshold calls for numerous numeric processing operations, which further complicate implementation of the method. One proposal to eliminate the aforementioned drawbacks is the detonation detecting method described in Patent Application WO-97/21084 filed on 29.11.1996 by the same Applicant and published on 12.06.1997. The detection method in the above patent application provides for wide-band filtering and rectifying the output signal of an acceleration sensor on the crankcase; integrating the rectified signal in a time window appropriately located with respect to the combustion considered; and effecting a logarithmic conversion of the integration result. Repeating the above operations for each combustion of the engine provides for defining a distribution of values, of which are calculated, at each combustion, a mean value and a numeric value correlated to the standard deviation on the basis of previous combustion values. On the basis of the numeric value and a predetermined threshold value which is constant alongside variations in engine speed, a detonation coefficient value is then calculated at each engine combustion, and which indicates the detonation "propensity" of the engine. DISCLOSURE OF INVENTION It is an object of the present invention to provide an improved detonation detecting method as compared with that described in Patent Application WO-97/21084. According to the present invention, there is provided a method of detecting and controlling detonation of an internal combustion engine, comprising the steps of: a) acquiring a vibration signal proportional to the intensity of vibration on the crankcase; b) wide-band filtering said vibration signal to generate a first intermediate signal; c) rectifying said first intermediate signal to generate a second intermediate signal; d) integrating said second intermediate signal to generate a first numeric value; e) calculating a logarithm of said first numeric value to obtain a second numeric value; f) calculating a mean value as a function of said second numeric value; characterized by comprising the steps of: g) calculating a difference between said second numeric value and said mean value to obtain a third numeric value ; h) comparing said third numeric value with a predetermined threshold value; m) determining the presence of detonation in the event said third numeric value has a first predetermined relationship with said predetermined threshold value; n) forming a detonation index indicating the behaviour of the engine in terms of detonation; and o) calculating, from said detonation index, a correction value to be added to the spark lead. BRIEF DESCRIPTION OF DRAWINGS A preferred, non-limiting embodiment of the present invention will be described by way of example with reference to the accompanying drawings, in which: Figures la and lb show two pressure cycle graphs in the absence and presence of detonation respectively; Figure 2 shows a simplified diagram of a combustion system featuring the detection method according to the invention; Figure 3 shows a flow chart of a first embodiment of the detection method according to the invention; Figures 4 and 5 show graphs of certain quantities relative to the detection method according to the invention. BEST MODE FOR CARRYING OUT THE INVENTION Number 1 in Figure 2 indicates as a whole a detonation detecting device connected to an internal combustion engine 2, of which only the crankcase 3 and cylinders 4 are shown. Detecting device 1 comprises a known acceleration sensor 5 (not described in detail) fitted to crankcase 3 of engine 2 and generating a vibration signal D correlated to the intensity of the vibration on crankcase 3; and a central processing unit 6 for receiving vibration signal D and implementing the operations relative to the method according to the present invention and described below with reference to the Figure 3 flow chart. The operations described with reference to Figure 3 are repeated for each combustion of engine 2 to continually monitor detonation in the various engine cycles, and each of the values calculated in a generic i- th repetition is indicated hereinafter by the subscript II » J_ * In the following description, reference is made to a counter I i associated with the cylinder 4 of engine 2 in which the combustion considered takes place. More specifically, each cylinder 4 is assigned a respective counter I, and, to identify the various counters I in the following description, each is followed by an "[n]" term indicating the n-th cylinder 4 with which the counter is associated. The counter values are initialized, e.g. set to 0, at the start of detection (e.g. when engine 2 is turned on) and are subsequently modified according to whether or not detonation is established in respective cylinders 4. More specifically, and as will be made clear in the following description, the values of a counter I define a detonation index of engine 2 indicating the behaviour of engine 2 in terms of detonation. As shown in Figure 3, after each combustion of engine 2 , vibration signal D is generated by acceleration sensor 5 and acquired by central processing unit 6 (block lb) . Vibration signal D is then wide-band filtered (block 11), rectified (block 12), and integrated in an appropriate time window (block 13). More specifically, the integration time window commences at the combustion phase of the cylinder 4 considered, immediately following the top dead-center position, and is of constant angular duration with respect to the speed of engine 2, whereas the filtration band of vibration signal D is calibrated experimentally and normally ranges between 5 and 20 kHz. The logarithm x L of the integration result is then calculated (block 14), and a corresponding mean value μ ^ is determined as a function of logarithm value x 1 (block 15). More specifically, mean value μ x is calculated according to the equation: where μ x is the mean value calculated at the i-th repetition; μ _ λ is the mean value calculated at the previous repetition i-1; x ^. is the logarithm value calculated at the i-th repetition; and λ x is a predefined parameter . As can be seen from the above equation, at the first repetition, the calculated mean value μ., ^ coincides with the calculated logarithm value x 1 , whereas, from the second repetition onwards, mean value μ x is updated on the basis of both the logarithm value x calculated at the i-th repetition, and the mean value μ ^x calculated at repetition i-1, each weighted according to parameter λ x . A peculiar characteristic of the above algorithm for updating mean value μ x is that of defining a filter with a very small time constant enabling the calculated mean value to rapidly follow variations in the x values obtained following logarithmic conversion of the integration result. This characteristic is extremely useful in transient states by ensuring precise tracking of vibration signal D. Calculation of mean value μ x is followed by calculating the value δ j of the difference between value X;, and mean value μ x (block 16). Value δ x is then compared with a predetermined threshold value δ 0 (block 17) to determine the presence or absence of detonation and accordingly modify the value of the counter I x [n] associated with the cylinder in question. According to the present invention, threshold value 5 0 is constant for all the repeat detection operations, is constant alongside variations in the speed of engine 2 , and is correlated to the standard deviation value σ of the distribution of values x x in the absence of detonation. In particular, the standard deviation value σ of tne distribution of values x in the absence of detonation is conveniently calculated on the basis of laboratory engine tests, by so setting the injection advance of the engine as to eliminate detonation. Given the standard deviation value σ, the corresponding threshold value δ 0 is then calculated and conveniently memorized in central processing unit 6. Tests have shown that the best threshold value δ 0 for reliably detecting detonation is roughly 3-4 times the standard deviation value σ. Block 17 determines whether value δ 1 is greater than or equal to predetermined threshold value δ 0 . If δ x is greater than or equal to threshold value δ 0 (YES output of block 17), the presence of detonation in the cylinder 4 in question is established (block 18) and the value of counter I ± [n] associated with the cylinder is modified (block 19) according to the equation: Φ] = I J" M i δ . - δ o) where M is a proportion coefficient, which is generally a function of the speed of engine 2 and the n-th cylinder considered at the time. In other words, by depending on - li ¬ the cylinder 4 considered, proportion coefficient M takes into account the transfer function of crankcase 3 in the combustion considered. More specifically, at the i-th repetition, the value of the n-th counter Ii[n] is calculated by decreasing the value of the n-th counter I 1 [n] at repetition i-1 by an amount proportional to the difference between value δ x and threshold value δ 0 . Conversely, if δ x is less than threshold value δ 0 (NO output of block 17), the absence of detonation in the cylinder 4 in question is established (block 20) and the value of counter Iχ[n] associated with the cylinder is modified (block 21) according to the equation: /,[/.] = /,_,[#.]+1 That is, at the i-th repetition, the value of the n-th counter I 1 [n] is calculated by increasing by one unit the value of the n-th counter ∑i[n] at repetition i-1 . As stated, the current value of a counter I 1 [n] defines a detonation index of engine 2 indicating the behaviour of engine 2 in terms of detonation, i.e. the tendency of engine 2 to detonate. By increasing the value of the counter when no detonation is detected in the respective cylinder, and decreasing the counter value when detonation is detected, it therefore follows that high positive counter values indicate a tendency of engine 2 either not to detonate or to detonate in temporally widely spaced engine cycles, whereas high negative counter values indicate a tendency of engine 2 to detonate. Once the value of counter I [n] s modified according to whether the presence or absence of detonation is established in the respective cylinder (block 19 or 21), the current value I [n] of the counter is used to determine an injection advance correction value CA ± [n] relative to the n-th cylinder 4 considered, and defining the best correction to be made to the injection advance of the cylinder to control detonation in the cylinder as required (block 22). The injection advance correction value CA i [n] is determined according to the equation: CA,[n] = K F -I,[n] where I [n] is the current value of the n-th counter at the i-th repetition; and K p is a factor defining how the advance correction is to be implemented. In particular, K p may be constant alongside variations in the speed of engine 2 , or may be a function of the speed of engine 2. Purely by way of example, Figure 4 shows a graph of the value I ± [n] assumed, in successive engine cycles, by a counter associated with a cylinder 4 of engine 2; and Figure 5 shows a graph of the advance correction value CA i [n] for the same cylinder, calculated according to the above equation and as a function of value Ii[n] of the same counter. In Figure 5, advance correction value CA ± [n] varies in discrete values of 0.5° as opposed to continuously. The advantages of the method according to the present invention are as follows . In particular, to determine the presence or absence of detonation, the method according to the present invention operates on the basis of the difference between logarithm value x ± and mean value μ i r with no need, as in the method described in Patent Application WO-97/21084, to calculate a numeric value correlated to the standard deviation of the distribution of values x x . This obviously provides for greatly simplifying the processing required to determine the detonation behaviour and effectively control operation of engine 2, as well as for reducing the number of components required. Moreover, the method according to the present invention models the detonation phenomenon by defining a threshold value δ 0 valid in all operating conditions and for all types of engine 2, thus further simplifying processing and reducing the number of components required. Moreover, the method according to the present invention provides for determining advance correction value CA ± ln] using a law of correction enabling the number of detonations present at any speed of engine 2 to be "measured out" as required, and which provides for obtaining an advance correction directly proportional to the amount of detonation that has taken place. Moreover, the method according to the present invention operates on the basis of a wide-band-filtered signal to make use of the maximum possible amount of information, even though this obviously involves also detecting the effects of other phenomena not related to detonation. Finally, the method according to the present invention employs a mean value μ x updating algorithm enabling the calculated mean value to rapidly follow any variations in the x L values obtained following logarithmic conversion of the integration result. Clearly, changes may be made to the method as described and illustrated herein without, however, departing from the scope of the present invention.

Claims: CLAIMS 1 ) A method of detecting and controlling detonation of an internal combustion engine (2), comprising the steps of: a) acquiring (10) a vibration signal (D) proportional to the intensity of vibration on the crankcase; b) wide-band filtering (11) said vibration signal (D) to generate a first intermediate signal; c) rectifying (12) said first intermediate signal to generate a second intermediate signal; d) integrating (13) said second intermediate signal to generate a first numeric value; e) calculating (14) a logarithm of said first numeric value to obtain a second numeric value (x x ) ; f) calculating (15) a mean value (μ as a function of said second numeric value (x ; characterized by comprising the steps of: g) calculating (16) a difference between said second numeric value (x x ) and said mean value (μ 1 ) to obtain a third numeric value (δ x ) ; h) comparing (17) said third numeric value (δ with a predetermined threshold value (δ 0 ) ; m) determining (18) the presence of detonation in the event said third numeric value (δ has a first predetermined relationship with said predetermined threshold value (δ 0 ) ; n) forming a detonation index indicating the behaviour of the engine in terms of detonation; and o) calculating, from said detonation index, a correction value to be added to the spark lead. 2) A method as claimed in Claim 1, characterized in tήat said step i) comprises the step of determining (18) said presence of detonation in the event said third numeric value (δ- is greater than or equal to said predetermined threshold value (δ 0 ) . 3 ) A method as claimed in Claim 1 or 2 , characterized by also comprising the step of: p) determining (20) the absence of detonation in the event said third numeric value {δ t ) does not have said first predetermined relationship with said predetermined threshold value (δ 0 ) . 4 ) A method as claimed in any one of the foregoing Claims, characterized by also comprising the steps of: q) calculating (19, 21) a detonation index { I ± ln] ) indicating the behaviour of said engine (2) in terms of detonation; and r) calculating (22) an injection advance correction value (CAi nJ) as a function of said detonation index ( Ii fn ] ) . 5 ) A method as claimed in Claim 4 , characterized in that said step r) comprises the step of calculating (22) said injection advance correction value ( CA [n] ) according to the equation: Gψ] = *,■/,[/.] where CA ± [n] is said injection advance correction value, I ± [n] is said detonation index, and K p defines a law of correction for the injection advance. 6 ) A method as claimed in Claim 4 or 5 , characterized in that said step q) comprises the step of: s) calculating (19) said detonation index ( I ± [n] ] as a function of said third numeric value (δ ± ) and of said predetermined threshold value (δ 0 ) in the event said third numeric value (δ. has said first predetermined relationship with said predetermined threshold value (δ 0 ). 7 ) A method as claimed in any one of the foregoing Claims from 4 to 6, characterized by repeating said steps from a) to r) to generate a number (i) of said second numeric values (x ± ) , of said mean values ( μ ) , of said third numeric values (δ ) , of said detonation indexes { I i [n] ) r and of said injection advance correction values [ CA^n] ) . 8 ) A method as claimed in Claim 7 , characterized in that each of said detonation indexes (I ' i ' nJ) is calculated (19), in the event said third numeric value (δ ± ) has said first predetermined relationship with said predetermined threshold value (δ 0 ) , according to the equation: /,[,,] = /,_^]- -(c5,-c5 0 ) where Ii [n] is the detonation index calculated at the i- th repetition, I^^n ] is the detonation index calculated at the previous repetition i-1 , δ ± is the third numeric value calculated at the i-th repetition, δ 0 is said predetermined threshold value, and M is a proportion coefficient. 9) A method as claimed in Claim 7 or 8, characterized in that each of said detonation indexes ( ϊi[n] ) is calculated (21), in the event said third numeric value (δ L ) does not have said first predetermined relationship with said predetermined threshold value (δ 0 ) , according to the equation: /,[/!] = /,_,[!!]+1 where Ii[n] is the detonation index calculated at the i- th repetition, and I ± ^ tn] is the detonation index calculated at the previous repetition i-1. 10) A method as claimed in any one of the foregoing Claims from 7 to 9, characterized in that each of said mean values (μ is calculated (15) according to the equation: where μ ± is the mean value calculated at the i-th repetition, μ ± _ x is the mean value calculated at the previous repetition i-1 , x t is the second numeric value calculated at the i-th repetition, and λ x is a predetermined parameter . 11) A method as claimed in any one of the foregoing Claims, characterized in that said step b) comprises the step of filtering said vibration signal (D) between 5 and 25 kHz. 12 ) A method as claimed in any one of the foregoing Claims, characterized in that said step d) comprises the step of integrating said second intermediate signal in a predetermined time window during the combustion phase. 13) A method as claimed in Claim 12, characterized in that said integration time window commences, during said combustion phase, immediately following the top dead-center position, and is of constant angular duration with respect to the speed of said engine (2). 14) A method as claimed in any one of the foregoing Claims, characterized in that said threshold value (δ 0 ) used in step h) is predetermined and constant alongside a variation in the speed of said engine (2). 15) A method as claimed in Claim 7, characterized in that said threshold value (δ 0 ) used in step h) is correlated to the standard deviation value of the distribution of second numeric values ( x x ) obtained in the absence of detonation.

Inventors: MORRA GIANLUCA, RE FIORENTIN STEFANO, REINAUDI GIANPIERO, TUBETTI PAOLO, EINAUDI, GIANPIERO, MORRA, GIANLUCA, RE FIORENTIN, STEFANO, TUBETTI, PAOLO
Applicants: FIAT RICERCHE, MORRA GIANLUCA, RE FIORENTIN STEFANO, REINAUDI GIANPIERO, TUBETTI PAOLO, C.R.F. SOCIETA' CONSORTILE PER AZIONI, EINAUDI, GIANPIERO, MORRA, GIANLUCA, RE FIORENTIN, STEFANO, TUBETTI, PAOLO
Application year: 1998
Main classification: G01L23
IPC codes: F02D45, F02P5, G01L23