EP0236207B1 - Procédé et système d'injection électronique à régulation par sonde lambda pour moteur à combustion interne - Google Patents

Procédé et système d'injection électronique à régulation par sonde lambda pour moteur à combustion interne Download PDF

Info

Publication number
EP0236207B1
EP0236207B1 EP87400342A EP87400342A EP0236207B1 EP 0236207 B1 EP0236207 B1 EP 0236207B1 EP 87400342 A EP87400342 A EP 87400342A EP 87400342 A EP87400342 A EP 87400342A EP 0236207 B1 EP0236207 B1 EP 0236207B1
Authority
EP
European Patent Office
Prior art keywords
value
sensor signal
richness
engine
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
EP87400342A
Other languages
German (de)
English (en)
Other versions
EP0236207A1 (fr
Inventor
Rémi Lefebre
Francis Prampolini
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Regie Nationale des Usines Renault
Original Assignee
Regie Nationale des Usines Renault
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Regie Nationale des Usines Renault filed Critical Regie Nationale des Usines Renault
Priority to AT87400342T priority Critical patent/ATE51681T1/de
Publication of EP0236207A1 publication Critical patent/EP0236207A1/fr
Application granted granted Critical
Publication of EP0236207B1 publication Critical patent/EP0236207B1/fr
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02D—CONTROLLING COMBUSTION ENGINES
    • F02D41/00—Electrical control of supply of combustible mixture or its constituents
    • F02D41/02—Circuit arrangements for generating control signals
    • F02D41/14—Introducing closed-loop corrections
    • F02D41/1438—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
    • F02D41/1477—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the regulation circuit or part of it,(e.g. comparator, PI regulator, output)
    • F02D41/1481—Using a delaying circuit
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02D—CONTROLLING COMBUSTION ENGINES
    • F02D41/00—Electrical control of supply of combustible mixture or its constituents
    • F02D41/02—Circuit arrangements for generating control signals
    • F02D41/14—Introducing closed-loop corrections
    • F02D41/1438—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
    • F02D41/1444—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases
    • F02D41/1454—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases the characteristics being an oxygen content or concentration or the air-fuel ratio
    • F02D41/1456—Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases the characteristics being an oxygen content or concentration or the air-fuel ratio with sensor output signal being linear or quasi-linear with the concentration of oxygen

Definitions

  • the invention relates to a probe-controlled electronic injection system and system for internal combustion engine, of the type comprising at least one injector whose opening time is controlled by an electronic computer as a function of engine operating parameters. and the state of the ⁇ probe.
  • the ⁇ probe is a sensor whose output voltage can switch between a high level (rich mixture) and a low level (lean mixture) located on either side of a threshold corresponding to the stoichiometric ratio (richness "1") .
  • the output signal from the X probe is formatted in the electronic injection computer and the resulting logical information is a rectangular signal to which, by convention, the value "+1" is assigned when it is at the level high and steam "-1" when it is low.
  • the regulation adapted to this kind of information is the superposition of a proportional type regulation and an integral type regulation.
  • Proportional correction makes it possible to increase the frequency of regulation, while integral correction makes it possible more particularly to pass from one operating point to another adapted differently in richness, that is to say to face a mismatch of wealth.
  • the invention provides a method and a system for injection with regulation by probe which are fundamentally different from traditional solutions while allowing the regulation frequency to be substantially increased.
  • the invention can be combined with some of these traditional solutions to further increase the efficiency of regulation.
  • the aims of the invention are achieved by means of a method of dosing the fuel supplied to an internal combustion engine by at least one injector controlled by an electronic computer associated with a probe delivering a signal capable of taking one or more other of two function states of the composition of the exhaust gases, according to which the electronic computer determines the opening time of the injector from a nominal time function of engine operating parameters and a term of proportional and integral correction as a function of the probe signal, and according to which a predictive estimate of the richness of the exhaust gases is carried out from engine operating parameters and from the pure delay, determined experimentally, between the injector and the probe, characterized in that at least one simulated probe signal is developed from said predictive wealth estimate, said correction term is developed at from the simulated probe signal and said correction term is periodically modified in response to the detection of a difference between the state of the measured probe signal and the state of a delayed simulated probe signal, the latter being obtained by a time shift of the first signal equal to said pure delay.
  • a first simulated probe signal is produced by comparison of the predictive wealth estimation with first high and low thresholds equal respectively to the high and low thresholds for change of state of the probe, the probe signal is obtained simulated delayed by a time shift of said first signal equal to said pure delay, a second simulated probe signal is produced by comparison of the predictive wealth estimate with second high and low thresholds respectively higher than the first high and low thresholds, and we develop said correction term from the second simulated probe signal.
  • a reference term representative of the correction to be made at said nominal time is developed to obtain a state of the probe representative of the richness "1" and a raw estimated value of richness is calculated as a function of the difference between the correction term and the reference term.
  • the invention also relates to an electronic injection system for implementing the method defined above, comprising at least one fuel injector on the engine intake side, a probe sensitive to the composition of the exhaust gases, sensors for measuring engine operating parameters and an electronic computer which controls the opening time of the injector as a function of said parameters and of the output signal of said probe, characterized in that said system comprises a digital value read-only memory pure delay addressable by the computer as a function of the air pressure at the engine intake.
  • FIG. 1 shows an internal combustion engine with positive ignition 1 equipped with an injector 2 on the intake side 3 and a catalyst 4 for purifying exhaust gas on the exhaust side 5.
  • the injector 2 is controlled by a programmed microcomputer 6 via a power circuit 7.
  • the microcomputer 6 determines the nominal time Tin of opening of the injector 2 as a function of the air pressure measured by a pressure sensor 8 placed on the intake side 3 and of the engine speed. This latter information is delivered by a sensor 9 in front of which the teeth of a target 10 pass in rotation with the crankshaft of the engine.
  • the target 10 can also be provided with one or more accidents placed in a predetermined angular position to provide angular position information via the sensor 9, or a second target associated with an additional sensor can be provided for this purpose.
  • the nominal time Tin can be corrected by the microcomputer 6 as a function of other information such as the temperature of the atmospheric air, the temperature of the engine cooling water, etc., which it possibly receives. on auxiliary inputs 11.
  • the nominal time Tin is also corrected on the basis of the information delivered by a probe ⁇ 12 disposed on the exhaust side 5, between the engine 1 and the catalyst 4.
  • the output signal from the probe ⁇ is shaped in the microcomputer 6 and it then presents the shape of the signal SA of FIG. 5.
  • This signal S ⁇ contains information on the residual oxygen content of the exhaust gases, and also on the momentary ratio of air and fuel of the mixture drawn in by the engine. .
  • the high and low levels of this signal S ⁇ , to which numerical values, "+1" and "-1" are respectively assigned, correspond to riches respectively higher and lower than the stoichiometric ratio (wealth "1").
  • the state of the probe ⁇ is not the instantaneous image of the richness of the mixture admitted into the engine because there is a pure delay between the injector 2 and the probe ⁇ 12.
  • This pure delay determined experimentally, is stored in the form of digital values in a read only memory 13 addressable by the computer 6 as a function of the air pressure at the intake of the engine.
  • the ROM 13 can be internal or external to the computer 6.
  • FIG. 1 relates to a four-cylinder engine comprising a single injector opening for a time Ti at each engine half-turn.
  • the invention applies to any type of positive-ignition engine, regardless of the number of injectors and cylinders which equip it.
  • the parameters for calculating the nominal time Tin of opening of the injector 2 are given only by way of examples and it is possible, among other things, to use an air flow sensor instead of the pressure sensor. 8 intake side 3. In this case, the memory 13 containing the pure delay digital values is addressed as a function of the air flow instead of the pressure.
  • the gain of the integral correction may be an increasing function of the time elapsed since the last tilting of the probe 12, for example a parabolic function if acl is of the following form: or
  • G is then a predetermined fixed coefficient involved in determining the gain of the integral correction.
  • the invention differs from the solutions recalled above for determining the correction term aci by the fact that for this purpose it does not call directly on the measured probe signal S ⁇ , but on a simulated probe signal Ss ⁇ .
  • the method according to the invention is in fact based on the fact that the injection time Ti is adjusted on the basis of an estimate of richness and that the observation of the measured probe signal SA serves to periodically readjust this estimate. This makes it possible to get rid of the pure delay between the injector and the probe and, therefore, of this step wait for its tilting to effect the proportional correction, which has the consequence of increasing the frequency of richness oscillation.
  • the term is representative of the difference in richness with respect to richness 1 at the level of the probe ⁇ at the instant n + m + 1, m representing the pure delay between the injector and the probe.
  • Blocks 13 and 14 represent the initial values of acl and a and the difference acl- is made at 15 ⁇ .
  • Block 17 represents the entire system of FIG. 1 and includes in particular as a quantity the injection time Ti and the angular position ⁇ m of the engine from which is deduced the number n of engine U-turns that have elapsed since the last tilting of the probe.
  • the output quantities of block 17 are the measured pressure P and the measured probe signal S ⁇ .
  • Block 18 represents a low-pass filtering of the pressure P and the difference PP is performed at 19. This difference is multiplied by the coefficient K 'in 20, the term K' (PP) being positive in acceleration and negative in deceleration, and making it possible to take into account the problems of wetting of the walls of the intake manifold by the fuel.
  • Block 23 represents the hysteresis of the probe ⁇ and reconstitutes at time n (half-turn n) a simulated probe signal S's ⁇ which is a predictive estimate of what the measured probe signal S ⁇ will be at instant n + m + 1.
  • block 24 represents the determination of the pure delay m as a function of the air pressure P measured at the intake of the engine.
  • Block 25 represents a pure delay m brought to the signal S's ⁇ (n), corresponding to a transfer function eK "/ NP, and the difference between the simulated probe signal S" s ⁇ (n is made at 26) ) and the measured probe signal S ⁇ (n). This difference is multiplied by the coefficient K2 at 27 to be reinjected at 22 as explained above. Furthermore, this same difference S "s ⁇ (n) -S ⁇ (n) is multiplied by K1 at 28 to be reinjected at 15.
  • the proportional and integral correction on the coefficient acl is carried out on the basis of a second simulated probe signal Ss ⁇ (n) produced by block 29 from the non-delayed richness estimate Re.
  • This block 29 has a higher hysteresis than block 23, which allows more direct tilting of the probe 12 since the wealth excursions are amplified.
  • the blocks 30 and 31 respectively represent the integral and proportional corrections and the difference obtained at 32 represents the term ⁇ cl which is subtracted at 33 from the initial ⁇ cl. We therefore obtain at the exit of 33 the term acl which is injected at 15 with the term a resulting from the difference made at 34 between the initial a (block 14) and the ⁇ ⁇ calculated (block 28).
  • Figure 4 is a flowchart of operation of the computer 6 which allows to implement the servo diagram of Figure 2.
  • Figure 3 is a flowchart of an initialization program which takes place when the engine starts.
  • step 40 Upon receipt of an initialization instruction (step 40), the computer 6 loads a predetermined value XX into a CPT counter which counts the engine U-turns (step 41).
  • Step 45 represents the end of the initialization program.
  • step 50 The sequence of the main program in FIG. 4 will now be described with reference also to the timing diagram in FIG. 5.
  • This program takes place on each detection of the passage of the motor through a predetermined angular position, for example the passage of a piston through the top dead center (step 50).
  • step 51 is a test to determine whether or not the engine is still in its starting phase. If this is the case, the counter CPT has not yet been counted down and the flag set in step 42 of the initialization program is always 0.
  • step 52 the correction term acl is conventionally calculated from the measured probe signal SA (n):
  • the reference term a is calculated: where x is a coefficient of predetermined fixed value.
  • the step 54 which follows consists in giving the simulated probe signal S "s ⁇ (n) the value presented by the measured probe signal S ⁇ (n) at the nth half-turn.
  • step 55 the delayed simulated probe signal S "s ⁇ (n) is then given the value presented by the simulated non-delayed probe signal S's ⁇ (nm) m half-turns earlier. This can be viewed on the FIG. 5 where it can be seen that the signal S "s ⁇ is offset with respect to the signal S's ⁇ of the pure delay m as a function of the pressure P.
  • step 56 the reference term a is calculated:
  • Steps 54 and 56 both lead to step 57 of calculating the term ⁇ : where k is a coefficient of predetermined fixed value.
  • the next step 59 is the calculation of the predictive estimate of non-delayed wealth Re:
  • Step 59 is followed by a series of tests to compare the wealth estimate Re with thresholds D + and D- on the one hand and D '+ and D'- on the other hand.
  • the thresholds D + and D- are symmetrical with respect to richness 1, as are the thresholds D '+ and D'- which are higher than the thresholds D + and D- respectively.
  • the thresholds D + and D '+ have been shown in FIG. 5, which corresponds to operation in a rich mixture, but we can immediately deduce the various signals which would be obtained in the event of operation in a lean mixture by comparison of the estimated wealth Re with the thresholds D- and D'-.
  • the first test 60 which follows step 59 consists in comparing Re to the threshold D +. If Re is greater than or equal to D +, the value Ss ⁇ (n) is set to +1 (step 61). Otherwise, we go to test 62 where we compare Re to the threshold D-. If Re is less than or equal to D-, the value Ss ⁇ (n) is set to the value -1 (step 63). Steps 61 and 63 or a negative response to test 62 lead to test 64 where Re is compared to the threshold D '+. If the response to this test is positive, the value S's ⁇ (n) is set to the value +1 (step 65), while in the opposite case we go to test 66 where we compare Re to the seil D'- .
  • Step 67 the value S's ⁇ (n) is set to the value -1 (step 67). Steps 65 and 67, as well as a negative response to test 66, lead to test 68. In the case of a negative response to tests 60 and 62, Ss ⁇ (n) retains the value it had at the time n-1 and, similarly, in the case of a negative response to tests 64 and 66, S's ⁇ (n) retains the value it had at time n-1.
  • Test 68 relates to the value of the flag. If the engine is still starting, the flag always has the value 0 fixed in step 42 of the initialization program and the response to test 68 is negative and leads to a test 69 concerning the content of the counter CPT initialized to the value XX in step 41 of the initialization program. In the start-up phase, the content of the CPT counter has still not been brought to 0 and the negative response to test 69 leads to step 70 where the CPT counter is decremented by one.
  • the next step 71 consists in calculating the injection time using the correction term ⁇ cl calculated in step 52:
  • the next step 72 marks the end of the program which remains pending the next interruption due to the passage of the motor through a predetermined angular position.
  • step 73 When the engine has made XX U-turns, the content of the CPT counter has been counted down to 0 and the response to test 69 is positive. The flag is then set to 1 (step 73), after which we go to step 71 for calculating the injection time Ti.
  • step 74 the correction term acl is calculated depending on the simulated probe signal Ss ⁇ :
  • the term ⁇ cl is always calculated from the simulated probe signal Ss ⁇ , although one can envisage a return to a traditional calculation based on the measured probe signal SA under certain specific operating conditions. of the motor.
  • FIG. 5 clearly shows the pure delay which exists between the injector 2 and the probe 12: it can be seen that the real richness at the level of the probe Rr, supposed to be initially at a level to facilitate the understanding of the phenomenon described, does not begins to increase as m half-turns after the increase in the richness of the mixture on admission due to the proportional correction introduced by the presence of a jump in the term ⁇ cl at the supposedly initial time I has appeared. against that the predictive wealth estimate Re begins to increase from the U-turn 1 until the U-turn p where it reaches the threshold D +. This causes a change in state of the simulated probe signal Ss ⁇ used for the calculation of aci which, therefore, immediately provides a proportional correction followed by an integral correction.
  • the simulated probe signal S's ⁇ has already changed state at half-turn n when the predictive richness estimate Re has reached the threshold D '+ below the threshold D +, so that the simulated probe signal delayed S "s ⁇ will also change state m half-turns later, namely at the half-turn n + m.
  • the real probe had switched some time earlier, at half turn q. This means that the predictive wealth estimate Re does not exactly coincide with the evolution of the real wealth Rr at the probe ⁇ and this results in a turn around q correction, both in the predictive estimate of wealth Re (coefficient K2), than of coefficient a (coefficient K1).
  • acl of proportional and integral correction is assumed to be calculated in a conventional manner as indicated by the above formula (2). It should however be noted that this term can be calculated just as well by formula (3) ensuring an integral correction of the parabolic type or by any other appropriate formula. Indeed, the invention does not lie in the actual calculation formula of this term, but in the use, for this purpose, of a simulated probe signal based on a prediction of the evolution of the richness of the gases. exhaust at the actual probe. In other words, the described method and system ensure richness looping from a signal simulated by an internal model and a recalibration of this internal model is carried out periodically from the observation of the state of the ⁇ probe.
  • FIG. 6 represents at different excitation frequencies of the term ⁇ cl the efficiency of a trifunctional catalyst of 54,000 miles as a function of the peak-to-peak amplitude of the richness oscillations at the entry of the catalyst pot. Efficiency is calculated as follows, expressed as a percentage:
  • FIGS. 7 and 8 relate to an internal combustion engine controlled respectively by a conventional fuel injection method and the method according to the invention.
  • These curves represent the spectral analysis of the engine period T expressed in milliseconds at idle speed. It appears that in the first case the fundamental line is located around 0.9 Hz, while it is close to 2 Hz with the method of the invention.
  • This increase in frequency results not only in a gain in the efficiency of the catalytic converter, but also in a reduction in low-frequency pumping of the engine speed at idle, from which it results an improvement in the vibrational comfort on the vehicle perceptible by a driver .

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)
  • Testing Of Engines (AREA)
  • Fuel-Injection Apparatus (AREA)

Description

  • L'invention concerne un procédé et un système d'injection électronique à régulation par sonde pour moteur à combustion interne, du type comportant au moins un injecteur dont te temps d'ouverture est commandé par un calculateur électronqiue en fonction de paramètres de fonctionnement du moteur et de l'état de la sonde λ.
  • La sonde λ est un capteur dont la tension de sortie peut basculer entre un niveau haut (mélange riche) et un niveau bas (mélange pauvre) situés de part et d'autre d'un seuil correspondant au rapport stoechiométrique (richesse "1"). Le signal de sortie de la sonde X est mis en forme dans le calculateur électronique d'injection et l'information logique qui en résulte est un signal rectangulaire auquel, par convention, on attribue la valeur "+1" lorsqu'il est au niveau haut et la vapeur "-1" lorsqu'il est au niveau bas.
  • La régulation adaptée à ce genre d'information est la superposition d'une régulation du type proportionnel et d'une régulation du type intégral. La correction proportionnelle permet d'accroître la fréquence de régulation, tandis que la correction intégrale permet plus particulièrement de passer d'un point de fonctionnement à un autre adapté différemment en richesse, c'est-à-dire de faire face à une désadaptation de richesse.
  • La régulation en boucle fermée de l'injection au moyen d'une sonde À est essentiellement utilisée lorsque le système d'échappement du moteur est équipé d'un catalyseur destiné à réduire les émissions de composants indésirables (polluants) des gaz d'échappement. En effet, cette régulation permet d'asservir la richesse du mélange air/carburant admis dans le moteur autour de la valeur À=1, ce qui est une condition indispensable à la combustion satisfaisante des composants toxiques par le catalyseur.
  • Le bon fonctionnement de ce dernier nécessite donc une régulation du mélange qui soit précise et présente la plus faible inertie possible. En d'autres termes, une fréquence de régulation ou de basculement de la sonde élevée favorise un fonctionnement efficace du catalyseur.
  • A cet effet, différentes solutions ont été proposées pour optimiser les termes de correction proportionnelle et intégrale et, malgré leur efficacité certaine, ces solutions se heurtent à certaines limites.
  • Les documents US―A―4.282.842 et US―A―4.111.171 décrivent un système de dosage de carburant comprenant un calculateur électronique associé à une sonde λ produisant un signal à deux états en fonction de la richesse de la composition des gaz d'échappement.
  • L'invention propose un procédé et un système d'injection à régulation par sonde qui se différencient fondamentalement des solutions traditionnelles tout en permettant d'accroître sensiblement la fréquence de régulation. En outre, l'invention peut être combinée avec certaines de ces solutions traditionnelles pour accroître encore l'efficacité de la régulation.
  • Les buts de l'invetnion sont atteints au moyen d'un procédé de dosage du carburant fourni à un moteur à combustion interne par au moins un injecteur commandé par un calculateur électronique associé à une sonde délivrant un signal susceptible de prendre l'un ou l'autre de deux états fonctions de la composition des gaz d'échappement, suivant lequel le calculateur électronique détermine le temps d'ouverture de l'injecteur à partir d'un temps nominal fonction de paramètres de fonctionnement du moteur et d'un terme de correction proportionnelle et intégrale fonction du signal de la sonde, et suivant lequel on effectue une estimation prédictive de richesse des gaz d'échappement à partir de paramètres de fonctionnement du moteur et du retard pur, déterminé expérimentalement, entre l'injecteur et la sonde, caractérisé en ce qu'on élabore au moins un signal de sonde simulé à partir de ladite estimation prédictive de richesse, on élabore ledit terme de correction à partir du signal de sonde simulé et on modifié périodiquement ledit terme de correction en réponse à la détection d'une différence entre l'état du signal de sonde mesuré et l'état d'un signal de sonde simulé retardé, ce dernier étant obtenu par un décalage temporel du premier signal égal audit retard pur.
  • Suivant une caractéristique, on élabore un premier signal de sonde simulé par comparaison de l'estimation prédictive de richesse avec des premiers seuils haut et bas égaux respectivement aux seuils haut et bas de changement d'état de la sonde, on obtient le signal de sonde simulé retardé par un décalage temporel dudit premier signal égal audit retard pur, on élabore un second signal de sonde simulé par comparaison de l'estimation prédictive de richesse avec des seconds seuils haut et bas supérieurs respectivement aux premiers seuils haut et bas, et on élabore ledit terme de correction à partir du second signal de sonde simulé.
  • Suivant une autre caractéristique, on élabore un terme de référence représentatif de la correction à apporter audit temps nominal pour obtenir un état de la sonde représentatif de la richesse "1" et on calcule une valeur brute estimée de richesse en fonction de la différence entre le terme de correction et le terme de référence.
  • L'invention à également pour objet un système d'injection électronique pour la mise en oeuvre du procédé défini ci-dessus, comprenant au moins un injecteur de carburant côté admission du moteur, une sonde sensible à la composition des gaz d'échappement, des capteurs de mesure de paramètres de fonctionnement du moteur et un calculateur électronique qui commande le temps d'ouverture de l'injecteur en fonction desdits paramètres et du signal de sortie de ladite sonde, caractérisé en ce que ledit système comprend une mémoire morte de valeur numérique de retard pur adressable par le calculateur en fonction de la pression d'air à l'admission du moteur.
  • D'autres avantages de l'invention ressortiront de la description qui va suivre d'un mode de sa réalisation donné uniquement à titre d'exemple et illustré par les dessins annexés sur lesquels:
    • La figure 1 est un schéma-bloc d'un système d'injection pour la mise en oeuvre du procédé suivant l'invention.
    • La figure 2 est un schéma fonctionnel illustrant la régulation prédictive en boucle fermée du procédé suivant l'invention.
    • La figure 3 est un organigramme illustrant le programme d'initialisation du calculateur d'injection.
    • La figure 4 est un organigramme de fonctionnement'du calculateur d'injection pour la mise en oeuvre du procédé de l'invention.
    • La figure 5 est un chronogramme montrant, en fonction du nombre N de demi-tours moteur, l'évolution d'un certain nombre de signaux représentatifs du fonctionnement du système d'injection suivant l'invention.
    • Les figures 6 à 8 sont des graphiques montrant respectivement l'efficacité d'un pot catalyseur et l'analyse spectrale de la période moteur sans et avec le procédé de l'invention.
  • Le schéma-bloc de la figure 1 montre un moteur à combustion interne à allumage commandé 1 équipé d'un injecteur 2 côté admission 3 et d'un catalyseur 4 d'épuration de gaz d'échappement côté échappement 5.
  • L'injecteur 2 est commandé par un microcalculateur programmé 6 par l'intermédiaire d'un circuit de puissance 7. Le microcalculateur 6 détermine le temps nominal Tin d'ouverture de l'injecteur 2 en fonction de la pression d'air mesurée par un capteur de pression 8 placé côté admission 3 et de la vitesse de rotation du moteur. Cette dernière information est délivrée par un capteur 9 devant lequel défilent les dents d'une cible 10 solidaire en rotation du vilebrequin du moteur. La cible 10 peut également être pourvue d'un ou plusieurs accidents placés dans une position angulaire prédéterminée pour fournir une information de position angulaire par l'intermédiaire du capteur 9, ou une deuxième cible associée à un capteur supplémentaire peuvent être prévus à cet effet.
  • Le temps nominal Tin peut être corrigé par le microcalculateur 6 en fonction d'autres informations telles que la température de l'air atmosphérique, la température de l'eau de refroidissement du moteur, etc., ..., qu'il reçoit éventuellement sur des entrées auxiliaires 11.
  • Le temps nominal Tin est également corrigé à partir de l'information délivrée par une sonde À 12 disposée côté échappement 5, entre le moteur 1 et le catalyseur 4. Le signal de sortie de la sonde À est mis en forme dans le microcalculateur 6 et il présente alors l'allure du signal SA de la figure 5. Ce signal Sλ contient une information sur la teneur en oxygène résiduel des gaz d'échappement, et également sur le rapport momentané d'air et de carburant du mélange aspiré par le moteur. Les niveaux haut et bas de ce signal Sλ, auxquels sont affectées respectivement des valeurs numériques, "+1" et "-1", correspondent à des richesses respectivement supérieure et inférieure au rapport stoechiométrique (richesse "1").
  • Comme on le verra dans la suite, l'état de la sonde À n'est pas l'image instantanée de la richesse du mélange admis dans le moteur car il existe un retard pur entre l'injecteur 2 et la sonde À 12. Ce retard pur, déterminé expérimentalement, est stocké sous forme de valeurs numériques dans une mémoire morte 13 adressable par le calculateur 6 en fonction de la pression d'air à l'admission du moteur. La mémoire morte 13 peut être interne ou externe au calculateur 6.
  • On considérera dans la suite que l'ensemble décrit à la figure 1, concerne un moteur à quatre cylindres comportant un seul injecteur s'ouvrant pendant un temps Ti à chaque demi-tour du moteur. Cependant, il doit être compris que l'invention s'applique à tout type de moteur à allumage commandé, quel que soit le nombre d'injecteurs et de cylindres qui l'équipent. De même, les paramètres de calcul du temps nominal Tin d'ouverture de l'injecteur 2 sont donnés uniquement à titre d'exemples et l'on peut, entre autres, utiliser un capteur de débit d'air au lieu du capteur de pression 8 côté admission 3. Dans ce cas, la mémoire 13 contenant les valeurs numériques de retard pur est adressée en fonction du débit d'air au lieu de la pression.
  • Le temps Ti d'ouverture de l'injecteur 2 calculé par le microcalculateur 6 est donné par la formule suivante:
    Figure imgb0001
    où
    • Tin représente le temps nominal d'ouverture calculé de manière classique en fonction des paramètres principaux et auxiliaires de fonctionnement du moteur mentionnés plus haut;
    • TiO est le temps nécessaire pour que l'injecteur commence à débiter après son excitation par le circuit de puissance 7;
    • acl est le terme de correction ou coefficient de régulation par sonde λ;
    • K est un coefficient de valeur-prédéterminée en fonction de la dynamique désirée sur la correction de Tin.
  • Dans le cas d'un microcalculateur 8 bits, on peut prendre K=256 et acl peut varier entre 0 et 256, sa valeur nominale étant 128. Par conséquent, la formule (1) peut s'écrire:
    Figure imgb0002
    le terme 1+acl/256 pouvant varier entre la valeur 1 (acl=0) et la valeur 2 (acl=256), la valeur nominale étant de 1,5 (acl=128).
  • Le coefficient acl détermine les corrections proportionnellç et intégrale et est généralement exprimé par la formule:
    Figure imgb0003
    où:
    • n représente le nombre de demi-tours effectués par le moteur depuis le dernier basculement de la sonde λ12;
    • H est un coefficient fixe ou variable déterminant l'amplitude de la correction proportionnelle;
    • G est un coefficient fixe ou variable déterminant le gain de la correction intégrale.
  • En variante, le gain de la correction intégrale peut être une fonction croissante du délai écoule depuis le dernier basculement de la sonde 12, par exemple une fonction parabolique si acl est de la forme suivante:
    Figure imgb0004
    où
  • G est alors un coefficient fixe prédéterminé intervenant dans la détermination du gain de la correction intégrale.
  • L'invention se distingue des solutions rappelées ci-dessus de détermination du terme de correction aci par le fait qu'à cet effet elle ne fait pas appel directement au signal de sonde mesuré Sλ, mais à un signal de sonde simulé Ssλ.
  • Le procédé suivant l'invention est en effet basé sur le fait que le temps d'injection Ti est réglé à partir d'une estimation de richesse et que l'observation du signal de sonde mesuré SA sert à réajuster périodiquement cette estimation. Ceci permet de s'affranchir du retard pur entre l'injecteur et la sonde et, donc, de ce pas attendre son basculement pour effectuer la correction proportionnelle, ce qui a pour conséquence d'augmenter la fréquence d'oscillation de richesse.
  • Si l'on pose en effet que a est la valeur de acl pour laquelle le temps d'injection conduit à un mélange de richesse 1, soit Ti(a)=Ti richesse 1, on peut écrire:
    Figure imgb0005
  • Compte tenu de la fonction de transfert entre l'injecteur et la sonde λ, on a:
    Figure imgb0006
  • Par conséquent, le terme
    Figure imgb0007
    est représentatif de l'écart de richesse par rapport à la richesse 1 au niveau de la sonde λ à l'instant n+m+1, m représentant le retard pur entre l'injecteur et la sonde.
  • Le procédé de régulation prédictive en boucle fermée mettant en oeuvre ce concept sera maintenant décrit en se reportant au schéma fonctionnel de la figure 2.
  • Les blocs 13 et 14 représentent les valeurs initiales de acl et a et l'on effectue en 15 la différence acl-α. Le bloc 16 représente un filtrage passe-bas à la sortie duquel on obtient le terme [β=(αcl-α).
  • Le bloc 17 représente l'ensemble du système de la figure 1 et comporte notamment comme grandeur le temps d'injection Ti et la position angulaire θm du moteur dont est déduit le nombre n de demi-tours moteur qui se sont écoulés depuis le dernier basculement de la sonde. Les grandeurs de sortie du bloc 17 sont la pression mesurée P et le signal de sonde mesuré Sλ. Le bloc 18 représente un filtrage passe-bas de la pression P et.l'on effectue en 19 la différence P-P. Cette différence est multipliée par le coefficient K' en 20, le terme K' (P-P) étant positif en accélération et négatif en décélération, et permettant de prendre en compte les problèmes de mouillage des parois du collecteur d'admission par le carburant. On effectue en 21 la différence [β-K'(P-P) et l'on retranche en 22 de cette différence le term K2[S"sλ(n)-Sλ(n)]. On obtient donc en 22 une grandeur Re qui constitue une estimation de richesse non retardée au niveau des gaz d'échappement.
  • Le bloc 23 représente l'hystérésis de la sonde λ et reconstitue à l'instant n (demi-tour n) un signal de sonde simulé S'sλ qui est une estimation prédictive de ce que sera le signal de sonde mesuré Sλ à l'instant n+m+1.
  • Par ailleurs, le bloc 24 représente la détermination du retard pur m en fonction de la pression d'air P mesurée à l'admission du moteur. Le bloc 25 représente un retard pur m apporté sur le signal S'sλ(n), correspondant à une fonction de transfert e-K"/N P, et l'on effectue en 26 la différence entre le signal de sonde simulé S"sÀ(n) et le signal de sonde mesuré SÀ(n). Cette différence est multipliée par le coefficient K2 en 27 pour être réinjectée en 22 comme expliqué précédemment. Par ailleurs, cette même différence S"sλ(n)-Sλ(n) est multipliée par K1 en 28 pour être réinjectée en 15.
  • Comme le montre la figure 5, un écart entre le signal de sonde simulé S"sλ(n) et le signal de sonde mesuré Sλ(n) à pour cause une mauvaise estimation, donc une erreur, sur le terme a (acl est connu), qui est donc corrigée par l'intermédiaire du coefficient K1 (correction intégrale), et a pour conséquence une erreur sur le calcul d'estimation de richesse en cours qui est corrigée par l'intermédiaire du coefficient K2 (correction intégrale).
  • La correction proportionnelle et intégrale sur le coefficient acl est effectuée à partir d'un deuxième signal de sonde simulé SsÀ(n) élaboré par le bloc 29 à partir de l'estimation de richesse non retardée Re. Ce bloc 29 possède une hystérésis plus importante que le bloc 23, ce qui permet des basculements plus francs de la sonde 12 puisque les excursions de richesse se trouvent amplifiées. Les blocs 30 et 31 représentent respectivement les corrections intégrale et proportionnelle et la différence obtenue en 32 représente le terme Δαcl qui est retranché en 33 du αcl initial. On obtient donc à la sortie de 33 le terme acl qui est injecté en 15 avec le terme a résultant de la différence effectuée en 34 entre le a initial (bloc 14) et le Δα calculé (bloc 28).
  • La figure 4 est un organigramme de fonctionnement du calculateur 6 qui permet de mettre en oeuvre le schéma d'asservissement de la figure 2. La figure 3 est un organigramme d'un programme d'initialisation qui se déroule lors du démarrage du moteur.
  • Dès réception d'une instruction d'initialisation (étape 40), le calculateur 6 charge une valeur prédéterminée XX dans un compteur CPT qui décompte les demi-tours moteur (étape 41 ). L'étape suivante 42 consiste à émettre un drapeau de valeur 0, puis on fixe les valeurs initiales de a et acl à l'étape 43 (a=acl=128) et de β et Re à l'étape 44 (β=Re=0). L'étape 45 représente la fin du programme d'initialisation.
  • Le déroulement du programme principal de la figure 4 sera maintenant décrit en se référant également au chronogramme de la figure 5. Ce programme se déroule à chaque détection du passage du moteur par une position angulaire prédéterminée, par exemple le passage d'un piston par le point mort haut (étape 50). L'étape suivante 51 est un test pour déterminer si le moteur est ou non encore dans sa phase de démarrage. Si tel est le cas, le compteur CPT n'a pas encore été décompté et le drapeau fixé à l'étape 42 du programme d'initialisation est toujours à 0. La réponse négative au test 51 conduit alors à l'étape 52 où le terme de correction acl est calculé de manière classique à partir du signal de sonde mesuré SA(n):
    Figure imgb0008
    A l'étape suivante 53, on calcule le terme de référence a:
    Figure imgb0009
    où x est un coefficient de valeur fixe prédéterminée. L'étape 54 qui suit consiste à donner au signal de sonde simulé S"sÀ(n) la valeur que présente le signal de sonde mesuré Sλ(n) au nième demi-tour.
  • Si, comme cela sera expliqué par la suite, la phase de démarrage du moteur est terminée et celui-ci a effectué un nombre de demi-tours au moins égal à celui XX fixé dans le compteur CPT à l'étape 41 du programme d'initialisation, la valeur du drapeau s'est trouvée fixée à 1 et la réponse au test 52 est positive. A l'étape 55 on donne alors au signal de sonde simulé retardé S"sÀ(n) la valeur que présentait le signal de sonde simulé non retardé S'sλ(n-m) m demi-tours plus tôt. Ceci peut être visualisé sur la figure 5 où l'on constate que le signal S"sλ est décalé par rapport au signal S'sλ du retard pur m fonction de la pression P. On calcule ensuite à l'étape 56 le terme de référence a:
    Figure imgb0010
  • Les étapes 54 et 56 conduisent toutes deux à l'étape 57 de calcul du terme β:
    Figure imgb0011
    où k est un coefficient de valeur fixe prédéterminée. On effectue ensuite en 58 un filtrage passe-bas sur la pression mesurée: P=P+x' (P-P) où x' est un coefficient de valeur fixe prédéterminée.
  • L'étape suivante 59 est le calcul de l'estimation prédictive de richesse non retardée Re:
    Figure imgb0012
  • L'étape 59 est suivie d'une série de tests pour comparer l'estimation de richesse Re à des seuils D+ et D- d'une part et D'+ et D'- d'autre part. Les seuils D+ et D- sont symétriques par rapport à la richesse 1, de même que les seuils D'+ et D'- qui sont supérieurs aux seuils D+ et D- respectivement. Pour des raisons de commodité, seuls les seuils D+ et D'+ ont été représentés sur la figure 5, ce qui correspond à un fonctionnement en mélange riche, mais on peut en déduire immédiatement les différentes signaux qui seraient obtenus en cas de fonctionnement en mélange pauvre par comparaison de la richesse estimée Re avec les seuils D- et D'-.
  • Le premier test 60 qui suit l'étape 59 consiste à comparer Re au seuil D+. Si Re est supérieur ou égal à D+, on fixe au signal Ssλ(n) la valeur +1 (étape 61). Dans le cas contraire, on passe au test 62 où l'on compare Re au seuil D-. Si Re est inférieur ou égal à D-, on fixe au signal Ssλ(n) la valeur -1 (étape 63). Les étapes 61 et 63 ou une réponse négative au test 62 conduisent au test 64 où l'on compare Re au seuil D'+. Si la réponse à ce test est positive, on fixe au signal S'sλ(n) la valeur +1 (étape 65), tandis que dans le cas contraire on passe au test 66 où l'on compare Re au seil D'-. Si la réponse à ce test est positive, on fixe au signal S'sλ(n) la valeur -1 (étape 67). Les étapes 65 et 67, ainsi qu'une réponse négative au test 66, conduisent au test 68. Dans le cas d'une réponse négative aux tests 60 et 62, Ssλ(n) conserve la valeur qu'il avait à l'instant n-1 et, de même, dans le cas d'une réponse négative aux tests 64 et 66, S'sλ(n) conserve la valeur qu'il avait à l'instant n-1.
  • Le test 68 porte sur la valeur du drapeau. Si l'on est toujours en phase de démarrage du moteur, le drapeau a toujours la valeur 0 fixée à l'étape 42 du programme d'initialisation et la réponse au test 68 est négative et conduit à un test 69 concernant le contenu du compteur CPT initialisé à la valeur XX à l'étape 41 du programme d'initialisation. Dans la phase de démarrage, le contenu du compteur CPT n'a toujours pas été ramené à 0 et la réponse négative au test 69 conduit à l'étape 70 où l'on décrémente d'une unité le compteur CPT.
  • L'étape suivante 71 consiste à calculer le temps d'injection en utilisant le terme de correction αcl calculé à l'étape 52:
    Figure imgb0013
  • . L'étape suivant 72 marque la fin du déroulement du programme qui reste dans l'attente de la prochaine interruption due au passage du moteur par une position angulaire prédéterminée.
  • Lorsque le moteur a effectué XX demi-tours, le contenu du compteur CPT à été décompté jusqu'à 0 et la réponse au test 69 est positive. Le drapeau est alors fixé à 1 (étape 73), après quoi l'on passe à l'étape 71 de calcul du temps d'injection Ti.
  • La phase de démarrage du moteur est alors achevée et au demi-tour suivant le test 51 conduit à dérouler les étapes 55 et 56. De même, la réponse positive au test 68 conduit à l'étape 74 où le terme de correction acl est calculé en fonction du signal de sonde simulé Ssλ:
    Figure imgb0014
  • Ensuite, tant que le moteur fonctionne, le terme αcl est toujours calculé à partir du signal de sonde simulé Ssλ, bien que l'on puisse envisager un retour à un calcul traditionnel fondé sur le signal de sonde mesuré SA dans certaines conditions particulières de fonctionnement du moteur.
  • La figure 5 montre clairement le retard pur qui existe entre l'injecteur 2 et la sonde 12: on constate en effet que la richesse réelle au niveau de la sonde Rr, supposée être initialement à un palier pour faciliter la compréhension du phénmène décrit, ne commence à augmenter que m demi-tours après que soit apparue l'augmentation de richesse du mélange à l'admission due à la correction proportionnelle introduite pas la présence d'un saut du terme αcl à l'instant supposé initial I. On voit par contre que l'estimation de richesse prédictive Re commence à croître dès le demi-tour 1 jusqu'au demi-tour p où elle atteint le seuil D+. Ceci provoque un changement d'état du signal de sonde simulé Ssλ utilisé pour le calcul de aci qui, de ce fait, apporte immédiatement une correction proportionnelle suivie d'une correction intégrale. Cependant, le signal de sonde simulé S'sλ a déjà changé d'état au demi-tour n lorsque l'estimation de richesse prédictive Re a atteint le seuil D'+ inférieur au seuil D+, de sorte que le signal de sonde simulé retardé S"sλ changera lui aussi d'état m demi-tours plus tard, à savoir au demi-tour n+m. Cependant, on a supposé dans l'exemple représenté que la sonde réelle avait basculé quelque temps plus tôt, au demi-tour q. Ceci signifie que l'estimation prédictive de richesse Re ne coïncide pas exactement avec l'évolution de la richesse réelle Rr au niveau de la sonde λ et il en résulte au demi-tour q une correction, tant de l'estimation prédictive de richesse Re (coefficient K2), que du coefficient a (coefficient K1).
  • Dans l'exemple de mise en oeuvre de l'invention par l'organigramme de la figure 4, le terme acl de correction proportionnelle et intégrale est supposé être calculé de manière classique comme indiqué par la formule (2) précitée. Il doit être cependant noté que ce terme peut être calculé tout aussi bien par la formule (3) assurant une correction intégrale de type parabolique ou par toute autre formule appropriée. En effet, l'invention ne réside pas dans la formule de calcul proprement dite de ce terme, mais dans l'utilisation, à cet effet, d'un signal de sonde simulé basé sur une prédiction de l'évolution de la richesse des gaz d'échappement au niveau de la sonde réelle. En d'autres termes, le procédé et le système décrits assurent un bouclage de richesse à partir d'un signal simulé par un modèle interne et un recalibrage de ce modèle interne est effectué périodiquement à partir de l'observation de l'état de la sonde λ. Dans l'exemple décrit, on doit donc distinguer la régulation du temps d'injection qui est base sur le signal de sonde simulé Ssλ de la régulation du modèle interne qui fait appel à l'autre signal de sonde simulé S'sλ et au signal de sonde simulé retardé S"sÀ. Bien entendu, il est possible d'effectuer la régulation du temps d'injection Ti directement sur le signal de sonde simulé S'sλ mais, comme indiqué précédemment, la solution décrite permet d'assurer des basculements plus francs de la sonde réelle 12 du fait que les seuils de richesse D+ et D- utilisés pour l'élaboration du signal de sonde simulé Ssλ sont supérieurs aux seuils réels D'+ et D'- de basculement de la sonde. D'autres modifications peuvent naturellement être apportées à l'exemple de réalisation décrit sans pour autant sortir du cadre et de l'objet de l'invention.
  • La figure 6 représente à différentes fréquences d'excitation du terme αcl l'efficacité d'un catalyseur trifonctionnel de 54.000 miles en fonction de l'amplitude crête à crête des oscillations de richesse à l'entrée du pot catalyseur. L'efficacité est calculée de la manière suivante, exprimée en pourcentage:
    Figure imgb0015
  • Ce graphique montre que l'efficacité est d'autant plus importante que la fréquence est élevée et que l'amplitude des oscillations de richesse est faible. Or, précisément, le procédé suivant l'invention permet d'accroître cette fréquence et, pour un gain H1 donné, de diminuer l'amplitude des oscillations de richesse.
  • Cette augmentation de fréquence est mise en évidence sur les figures 7 et 8 qui concernent un moteur à combustion interne commandé respectivement par un procédé d'injection de carburant classique et le procédé suviant l'invention. Ces courbes représentent l'analyse spectrale de la période moteur T exprimée en millisecondes au régime de ralenti. Il en ressort que dans le premier cas la raie fondamentale est située autour de 0,9 Hz, alors qu'elle est voisine de 2 Hz avec le procédé de l'invention.
  • Cette augmentation de fréquence se traduit non seulement par un gain sur l'efficacité du pot catalyseur, mais également par une réduction du pompage basse fréquence du régime moteur au ralenti, d'où il résulte une amélioration du confort vibratoire sur véhicule perceptible par un conducteur.

Claims (15)

  1. . 1. Procédé de dosage du carburant fourni à un moteur à combustion interne (1) par au moins un injecteur (2) commandé par un calculateur électronique (6) associé à une sonde (12) délivrant un signal (Sλ) suscpetibie de prendre l'un ou l'autre de deux états fonctions de la composition des gaz d'échappement, suivant lequel le calculateur électronique (6) détermine le temps d'ouverture (Ti) de l'injecteur (2) à partir d'un temps nominal (Tin) fonction de paramètres de fonctionnement du moteur (1) et d'un terme (acl) de correction proportionnelle et intégrale fonction de l'état du signal de la sonde, et suivant lequel on effectue une estimation prédictive (Re) de richesse des gaz d'échappement à partir de paramètres de fonctionnement du moteur et du retard pur (m), déterminée expérimentalement entre l'injecteur (2) et la sonde (12), caractérisé en ce qu'on élabore au moins un signal de sonde simulé (Ssλ) à partir de ladite estimation prédictive de richesse, on élabore ledit terme de correction (αcl) à partir du signal de sonde simulé (Ssλ) et on modifie périodiquement ledit terme de correction (acl) en réponse à la détection d'une différence entre l'état du signal de sonde mesuré (SA) et l'état d'un signal de sonde simulé retardé (S"sÀ), ce dernier étant obtenu par un décalage temporel du premier signal (S'sA) égal audit retard pur (m).
  2. 2. Procédé selon la revendication 1, caractérisé en ce qu'on élabore un premier signal de sonde simulé (S'sλ) par comparaison de l'estimation prédictive de richesse (Re) avec des premiers seuils haut et bas (D'+, D'-) égaux respectivement aux seuils haut et bas de changement d'état de la sonde, on obtien le signal de sonde simulé retardé (S"sλ) par un décalage temporel dudit premier signal (S'sλ) égal audit retard pur (m), on élabore un second signal de sonde simulé (Ssλ) par comparaison de l'estimation prédictive de richesse (Re) avec des seconds seuils haut et nas (D+, D-) supérieurs respectivement aux premiers seuils haut et bas, et on élabore ledit terme de correction (acl) à partir du second signal de sonde simulé (SsÀ).
  3. 3. Procédé selon la revendication 2, caractérisé en ce qu'on élabore un terme de référence (a) représentatif de la correction à apporter audit temps nominal (Tin) pour obtenir un état de la sonde représentatif de la richesse "1" et on calcule une valeur brute estimée de richesse (β) en fonction de la différence entre le terme de correction (acl) et le terme de référence (a).
  4. 4. Procédé selon la revendication 3, caractérisé en ce qu'on élabore ledit terme de référence (a) en fonction de la différence entre l'état du signal de sonde simulé retardé (S"sλ) et l'état du signal de sonde mesuré (Sλ).
  5. 5. Procédé selon l'une quelconque des revendications 3 et 4, caractérisé en ce que l'estimation prédictive de richesse (Re) est égale à la somme de ladite valeur brute (β), d'un premier terme (K'(P-P)) fonction de la pression d'air à l'admission du moteur et d'un second terme (K2(S"sλ-Sλ)) fonction de la différence entre l'état du signal de sonde simulé retardé (S"sλ) et l'état du signal de sonde mesuré (Ssλ).
  6. 6. Procédé selon la revendication 5, suivant lequel on détecte l'état de la sonde et on calcule le temps (Ti) d'ouverture de l'injecteur cycliquement à chaque tour ou fraction de tour du moteur, caractérisé en ce que l'estimation prédictive de richesse Re est exprimée par:
    Figure imgb0016
    où:
    β est la valeur brute estimée de richesse,
    K' est un coefficient de valeur prédéterminée,
    P est la pression d'air à l'admission du moteur,
    P représente la pression P vue à travers un filtre passe-bas,
    K2 est un coefficient de valeur prédéterminée,
    Sλ(n) est la valeur que présente le signal de sonde mesuré n cycles après le dernier basculement de la sonde,
    S"λ(n) est la valeur que présente le signal de sonde simulé retardé au cycle n,
    les signaux SÀ(n) et S"sλ(n) étant susceptibles de prendre les valeurs +1 ou -1.
  7. 7. Procédé selon la revendication 6, caractérisé en ce que la valeur brute estimée de richesse au cycle n est exprimée par:
    Figure imgb0017
    où:
    β(n-1) est la valeur brute estimée de richesse au cycle n-1,
    k est un coefficient de valeur prédéterminée,
    acl(n-1) est le terme de correction calculé au cycle n-1,
    a(n) est le terme de référence calculé au cycle n.
  8. 8. Procédé selon la revendication 7, caractérisé en ce que la valeur du terme de référence au cycle n est exprimée par:
    Figure imgb0018
    où:
    a(n-1) est la valeur du terme de référence au cycle n-1,
    K1 est un coefficient de valeur prédéterminée.
  9. 9. Procédé selon l'une quelconque des revendications 2 à 8, caractérisé en ce que la valeur du terme de correction n cycles après le dernier basculement de la sonde est exprimée par:
    Figure imgb0019
    où
    αcl(n-1) est la valeur du terme de correction au cycle n-1,
    G1 est un coefficient intervenant dans la détermination du gain de la correction intégrale,
    H1 est un coefficient déterminant l'amplitude de la correction proportionnelle,
    Ssλ(n) est la valeur que présente le second signal de sonde simulé au cycle n,
    Ssλ(n-1) est la valeur que présente le second signal de sonde simulé au cycle n-1,
    le signal Ssλ étant susceptible de prendre les valeurs +1 ou -1.
  10. 10. Procédé selon l'une quelconque des revendications 2 à 8, caractérisé en ce que la valeur du terme de correction n cycles après le dernier basculement de la sonde est exprimée par:
    Figure imgb0020
    où:
    acl(n-1) est la valeur du terme de correction au cycle n-1,
    H2 est un coefficient déterminant l'amplitude de la correction proportionnelle,
    G2 est un coefficient intervenant dans la détermination de la correction intégrale,
    Ssλ(n) est la valeur que présente le second signal de sonde simulé au cycle n,
    Ssλ(n-1) est la valeur que présente le second signal de sonde simulé au cycle n-1,
    le signal Ssλ étant susceptible de prendre les valeurs +1 ou -1.
  11. 11. Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce que ledit retard pur (m) est une fonction de la pression d'air (P) à l'admission du moteur.
  12. 12. Procédé selon la revendication 11, caractérisé en ce que ledit retard pur (m) est une valeur numérique exprimée en nombre de tours ou fractions de tour moteur.
  13. 13. Système de dosage du carburant fourni à un moteur à combustion interne (1), selon l'une quelconque des revendications 1 à 12 par au moins un injecteur (2) commandé par un calculateur électronique (6) associé à une sonde (12) délivrant un signal (Sλ) susceptible de prendre l'un ou l'autre de deux états fonctions de la composition des gaz d'échappement, suivant lequel le calculateur électronique (6) comprenant un premier moyen qui détermine le temps d'ouverture (Ti) de l'injecteur (2) à partir d'un temps nominal (Tin) fonction de paramètres de fonctionnement du moteur (1) et d'un terme (acl) de correction proportionnelle et intégrale fonction de l'état du signal de la sonde, et suivant lequel un deuxième moyen effectue une estimation prédictive (Re) de richesse des gaz d'échappement à partir de paramètres de fonctionnement du moteur et du retard pur (m), déterminée expérimentalement entre l'injecteur (2) et la sonde (12), caractérisé en ce qu un troisième moyen élabore au moins un signal de sonde simulé (Ssλ) à partir de ladite estimation prédictive de richesse, un quatrième moyen élabore ledit terme de correction (αcl) à partir du signal de sonde simulé (Ssλ) et modifie périodiquement ledit terme de correction (acl) en réponse à la détection d'une différence entre l'état du signal de sonde mesuré (Sλ) et l'état d'un signal de sonde simulé retardé (S"sÀ), ce dernier étant obtenu par un moyen de décalage temporel du premier signal (S'sÀ) égal audit retard pur (m).
  14. 14. Système d'injection électronique pour la mise en oeuvre du procédé selon la revendication 13, comprenant au moins un injecteur de carburant côté admission du moteur, une sonde sensible à la composition des gaz d'échappement, des capteurs de mesure de paramètres de fonctionnement du moteur et un calculateur électronique qui commande le temps d'ouverture de l'injecteur en fonction desdits paramètres et du signal de sortie de ladite sonde, caractérisé en ce que ledit système comprend une mémoire morte (13) de valeurs numériques de retard pur (m) adressable par ledit calculateur (6) en fonction de la pression d'air (P) à l'admission du moteur (1).
EP87400342A 1986-02-25 1987-02-17 Procédé et système d'injection électronique à régulation par sonde lambda pour moteur à combustion interne Expired - Lifetime EP0236207B1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
AT87400342T ATE51681T1 (de) 1986-02-25 1987-02-17 Elektronisches einspritzverfahren und -system mit lambdasondenregelung fuer brennkraftmaschinen.

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR8602557 1986-02-25
FR8602557A FR2594890B1 (fr) 1986-02-25 1986-02-25 Procede et systeme d'injection electronique a regulation par sonde l pour moteur a combustion interne

Publications (2)

Publication Number Publication Date
EP0236207A1 EP0236207A1 (fr) 1987-09-09
EP0236207B1 true EP0236207B1 (fr) 1990-04-04

Family

ID=9332492

Family Applications (1)

Application Number Title Priority Date Filing Date
EP87400342A Expired - Lifetime EP0236207B1 (fr) 1986-02-25 1987-02-17 Procédé et système d'injection électronique à régulation par sonde lambda pour moteur à combustion interne

Country Status (5)

Country Link
US (1) US4766871A (fr)
EP (1) EP0236207B1 (fr)
AT (1) ATE51681T1 (fr)
DE (1) DE3762145D1 (fr)
FR (1) FR2594890B1 (fr)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3810829A1 (de) * 1988-03-30 1989-10-12 Bosch Gmbh Robert Verfahren und vorrichtung zur lambdaregelung
KR960016086B1 (ko) * 1991-03-28 1996-11-27 미쯔비시지도오샤고오교오 가부시기가이샤 내연기관의 제어장치
JP2541051B2 (ja) * 1991-09-27 1996-10-09 株式会社島津製作所 限界電流式酸素濃度測定装置
IT1250530B (it) * 1991-12-13 1995-04-08 Weber Srl Sistema di controllo della quantita' di carburante iniettato per un sistema di iniezione elettronica.
US5305727A (en) * 1992-06-01 1994-04-26 Ford Motor Company Oxygen sensor monitoring
DE69326217T3 (de) * 1992-06-12 2009-11-12 Toyota Jidosha Kabushiki Kaisha, Toyota-shi Abgasemissionssteuerungssystem für verbrennungsmotoren
US5600056A (en) * 1994-06-20 1997-02-04 Honda Giken Kogyo Kabushiki Kaisha Air/fuel ratio detection system for multicylinder internal combustion engine
WO1996035048A1 (fr) * 1995-05-03 1996-11-07 Siemens Aktiengesellschaft Procede de regulation lambda d'un cylindre individuel d'un moteur a combustion interne multi-cylindre
FR2749350B1 (fr) * 1996-06-03 1998-07-10 Renault Systeme de regulation de la richesse par mode de glissement
FR2749613B1 (fr) * 1996-06-11 1998-07-31 Renault Systeme de regulation de la richesse dans un moteur a combustion interne
FR2849112B1 (fr) * 2002-12-18 2005-02-04 Renault Sa Procede de commande d'elements d'execution de fonctions elementaires de moteur a combustion interne

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5854253B2 (ja) * 1975-05-12 1983-12-03 日産自動車株式会社 クウネンピセイギヨソウチ
JPS6045297B2 (ja) * 1977-07-22 1985-10-08 株式会社日立製作所 内燃機関の燃料制御装置
JPS562548A (en) * 1979-06-22 1981-01-12 Nissan Motor Co Ltd Controller for air fuel ratio of internal combustion engine
DE3036107C3 (de) * 1980-09-25 1996-08-14 Bosch Gmbh Robert Regeleinrichtung für ein Kraftstoffzumeßsystem
US4359993A (en) * 1981-01-26 1982-11-23 General Motors Corporation Internal combustion engine transient fuel control apparatus
US4397278A (en) * 1981-04-03 1983-08-09 Ford Motor Company Air fuel ratio control using time-averaged error signal
JPS57193743A (en) * 1981-05-25 1982-11-29 Mitsubishi Electric Corp Oxygen sensor receiving device
US4548185A (en) * 1984-09-10 1985-10-22 General Motors Corporation Engine control method and apparatus
DE3441390A1 (de) * 1984-11-13 1986-05-15 M.A.N. Maschinenfabrik Augsburg-Nürnberg AG, 8000 München Verfahren zur regelung der schadstoffreduzierung bei gasmotoren

Also Published As

Publication number Publication date
DE3762145D1 (de) 1990-05-10
FR2594890B1 (fr) 1990-03-09
US4766871A (en) 1988-08-30
EP0236207A1 (fr) 1987-09-09
ATE51681T1 (de) 1990-04-15
FR2594890A1 (fr) 1987-08-28

Similar Documents

Publication Publication Date Title
EP0236207B1 (fr) Procédé et système d'injection électronique à régulation par sonde lambda pour moteur à combustion interne
FR2712388A1 (fr) Procédé et dispositif pour former un signal simulé de la température des gaz d'échappement de la sonde des gaz d'échappement ou du catalyseur.
FR2551798A1 (fr) Procede de commande d'alimentation en combustible d'un moteur a combustion interne immediatement apres le demarrage
FR2768181A1 (fr) Moteur a combustion interne notamment pour vehicule automobile
FR2466633A1 (fr) Systeme de commande de synchronisation d'etincelles pour moteurs a combustion interne
FR2486161A1 (fr) Procede et systeme pour regler l'avance a l'allumage d'un moteur en fonction de la presence ou de l'absence du cliquetis ou cognement dans ce moteur
FR2877695A1 (fr) Procede de gestion d'un moteur a combustion interne et dispositif pour sa mise en oeuvre
FR2688545A1 (fr) Systeme pour commander un moteur a combustion interne.
FR2780447A1 (fr) Procede de correction de la duree d'injection dans un moteur a combustion interne a systeme d'air secondaire
FR2785555A1 (fr) Procede de regeneration de sulfates d'un catalyseur a accumulation nox pour moteur a combustion interne a regime pauvre
EP0175596A1 (fr) Procédé adaptatif de régulation de l'injection d'un moteur à injection
FR2529255A1 (fr) Procede de commande du rapport air-combustible pour un moteur a combustion interne dans des regions de fonctionnement a faible charge
FR2548272A1 (fr) Procede de commande d'injection de combustible a l'acceleration d'un moteur a combustion interne a plusieurs cylindres du type a injection sequentielle
FR2541728A1 (fr) Procede de reglage du ralenti par reaction dans un moteur a combustion interne
EP0961875B1 (fr) Procede de commande de purge d'un pot catalytique de traitement des gaz d'echappement d'un moteur a combustion interne
FR2790282A1 (fr) Procede et appareil de mise en oeuvre d'un moteur a combustion interne notamment d'un vehicule avec un melange carburant/air pauvre
EP0636778B1 (fr) Procédé et dispositif de correction de la durée d'injection en fonction du débit de purge d'un circuit de purge à canister, pour moteur à injection
FR3101673A1 (fr) Procédé de réglage de la richesse d’un moteur à combustion interne à allumage commandé
FR2799234A1 (fr) Procede de surveillance d'un moteur a combustion interne
EP0146426B1 (fr) Procédé de commande d'un moteur à injection de carburant régulée par sonde lambda et à allumage commandé
FR2851298A1 (fr) Procede de commande d'un moteur thermique, appareil de commande et/ou de regulation pour un moteur thermique, programme d'ordinateur et support de memoire electrique d'un moteur thermique
FR2746851A1 (fr) Procede et dispositif de regulation en boucle fermee de la richesse d'un melange air/carburant destine a l'alimentation d'un moteur a combustion interne
FR2527691A1 (fr) Procede permettant de commander les dispositifs de commande de moteurs a combustion interne immediatement apres la fin d'une coupure de carburant
FR2868129A1 (fr) Procede de gestion d'un moteur a combustion interne
FR2634824A1 (fr) Procede de regulation de la richesse d'un melange air-carburant d'alimentation d'un moteur a combustion interne, du type boucle fermee, sans exploitation de mesure physique de richesse

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 19870220

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT CH DE GB IT LI SE

17Q First examination report despatched

Effective date: 19880229

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AT CH DE GB IT LI SE

REF Corresponds to:

Ref document number: 51681

Country of ref document: AT

Date of ref document: 19900415

Kind code of ref document: T

ITF It: translation for a ep patent filed

Owner name: JACOBACCI & PERANI S.P.A.

REF Corresponds to:

Ref document number: 3762145

Country of ref document: DE

Date of ref document: 19900510

GBT Gb: translation of ep patent filed (gb section 77(6)(a)/1977)
PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DE

Payment date: 19901213

Year of fee payment: 5

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 19910117

Year of fee payment: 5

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: CH

Payment date: 19910118

Year of fee payment: 5

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: AT

Payment date: 19910124

Year of fee payment: 5

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: SE

Payment date: 19910128

Year of fee payment: 5

ITTA It: last paid annual fee
26N No opposition filed
PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GB

Effective date: 19920217

Ref country code: AT

Effective date: 19920217

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Effective date: 19920218

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Effective date: 19920229

Ref country code: CH

Effective date: 19920229

GBPC Gb: european patent ceased through non-payment of renewal fee
REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DE

Effective date: 19921103

EUG Se: european patent has lapsed

Ref document number: 87400342.9

Effective date: 19920904

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES;WARNING: LAPSES OF ITALIAN PATENTS WITH EFFECTIVE DATE BEFORE 2007 MAY HAVE OCCURRED AT ANY TIME BEFORE 2007. THE CORRECT EFFECTIVE DATE MAY BE DIFFERENT FROM THE ONE RECORDED.

Effective date: 20050217