IGNITION TIMING CONTROL FOR AN INTERNAL COMBUSTION ENGINE

Abstract: An ignition timing control for an internal combustion engine by pressure ratio management. The ratio of combustion chamber pressure to motored pressure is determined for several predetermined crankshaft rotational angles. These pressure ratios are used to control ignition timing for MBT, diluent gas flow rate (EGR) and fuel balance among combustion chambers. The control requires from three to five pressure samples for each combustion event, depending on sensor and control options. The pressure sensors need not sense absolute pressure if they are linear and do not change gain or bias significantly during one cycle.

Description: This invention relates to ignition timing control in an internal combustion engine, and particularly to the control of ignition timing in response to the ratio of combustion chamber pressure to motored pressure as a function of crankshaft rotational angle. Currently, control of internal combustion engines is based on sensing a set of variables such as coolant temperature, manifold pressure (MAP), engine speed and exhaust oxygen concentration and adjusting ignition timing, exhaust gas recirculation (EGR) rate and fuel flow to a prescribed calibration for a baseline engine. The problem with this approach is that, because of manufacturing differences and aging, the engine being controlled is not necessarily the same as the baseline engine used for reference. One approach to the solution of this problem requires extensive effort in the design of the engine to ensure both that all combustion chambers behave the same and that components will not significantly alter characteristics with aging. But even with such design and tight quality control in the manufacturing process, the required uniformity and stability over time of such engines cannot always be achieved. A second approach to the solution of the problem is to implement a control with the capability to adjust for these differences and changes. Such a control is possible using pressure sensors for the combustion chamber and applying feedback control to ignition timing, dilution gas rate and fuel rate. In a typical current engine control, the three controlled combustion parameters are spark timing (or fuel injection timing in a diesel engine), EGR rate and air/fuel ratio. The first affects the timing of the initiation of the combustion process while the latter two affect the speed and duration of the combustion process. In current practice, air/fuel ratio is generally controlled in closed loop by an exhaust oxygen sensor to produce a constant stoichiometric ratio for emission control by oxidizing and reducing catalysts in the exhaust system. Since the efficiency of one or the other catalyst falls rapidly as the air/fuel ratio strays even slightly from stoichiometric in either direction, this parameter must be strictly controlled and is not available for maximizing power or fuel efficiency. EGR is generally controlled by a combination of different parameters such as exhaust backpressure, engine coolant temperature, engine speed, throttle position or manifold pressure and has proven difficult to control accurately. Spark or injection timing is generally determined from a stored table addressed with engine speed and load parameters with additional retard, in some cases, in response to a knock, MAP or throttle movement detector. Closed loop combustion control has been suggested in the prior art in various forms. LPP spark (location of peak pressure) timing controls have been proposed, in which spark timing is controlled to maintain a predetermined location of peak combustion pressure. This stability of the timing of peak combustion pressure has been found to produce MBT (minimum ignition advance for best torque) operation for many engine operating conditions. The location of peak pressure may be sensed using headbolts with an embedded piezoelectric material which responds to the stresses created in reaction to the pressure on the engine cylinder head. However, this approach has difficulties when the combustion charge is highly dilute or the engine is under light load. Some arrangements have been suggested which adjust ignition timing to control maximum absolute combustion pressure with respect to some predetermined reference level. There is, in addition, a suggestion in the US patent No. 4,449,50l that ignition timing be controlled to maintain the ratio of maximum combustion chamber pressure to maximum motored pressure in accordance with a stored table addressed by engine speed and load factors. These approaches, however, involve stored references determined for a baseline engine which may not be the desired reference for any particular engine at any given time. Other approaches have involved sensing flame position using ion gaps or engine speed variations. These approaches have drawbacks related, respectively, to cyclic variability in flame position and to difficulties with vehicle transients. If a measure of combustion pressure is available, several other approaches have been suggested. With unlimited computing capability, one might be able to perform a detailed heat release analysis or to obtain combustion duration to estimate air/fuel ratio. However, the computational requirements of these schemes prohibits their application for the foreseeable future. The closest prior art is believed to be the above mentioned US patent. In this publication, a pressure ratio is determined and used to determine ignition timing. However, the pressure ratio is not the ratio of combustion chamber pressure to motored pressure at a single crankshaft rotational angle, but is the ratio of maximum combustion chamber pressure, which occurs at MBT at about l5 degrees atdc (defined below), to maximum motored pressure, which occurs at zero degrees atdc. This ratio is computed and compared with a stored, desired value, with ignition timing adjusted in closed loop to maintain this value. This pressure ratio must be compared to a reference which is obtained in a "ratio map" of values stored in memory as a function of engine speed and load factor such as fuel quantity. This ratio map is necessary because the absolute value of the pressure ratio at one crankshaft rotational angle will vary with such factors. In the present invention, there is no need to reference such factors since the pressure ratio is related to final pressure ratio at complete combustion. In other words the reference is computed during each combustion event from the same pressure sensor rather than having to look it up in a table addressed by the readings of several other sensors. The present invention therefore has advantages in cost, since memories and sensors are minimised, as well as accuracy, since the reference is determined in the same combustion event under the same engine age and environmental status. Further, US patent no. 4,449,50l requires a calibrated pressure sensor as fixed stored reference values are used. In a preferred embodiment of the present invention, an extra pressure measurement in each combustion event eliminates the need for a calibrated pressure sensor by, in effect, recalibrating the pressure sensor at each combustion event, as long as the pressure sensor is linear and does not vary its gain appreciably during one combustion event. US patent no. 4,449,50l does not determine or recognise the usefulness of the motored pressure and a useful combustion pressure at substantially complete combustion, which values are key to the present invention in establishing a reference for determining the fraction of fuel burned at the earlier crankshaft rotational angle. Further, it does not compute the final pressure ratio which provides an accurate reference for a particular combustion event relating to computations to the actual total fuel burned. The primary object of this invention is to control the operation of an internal combustion engine from combustion chamber pressures measured at a few predetermined crankshaft rotational angles during an engine cycle. The main aim being to maintain ignition timing at an optimum value (MBT) by fixing the point at which a predetermined fraction of the fuel is burnt at a predetermined crankshaft rotational angle. To this end, an ignition timing control in accordance with the present invention is characterised by the features specified in the characterising portion of Claim l. In particular, the ignition timing for MBT may be determined in response to the computed ratio of first and second pressure ratios, where the latter are the ratios of sensed combustion chamber pressure to motored pressure at predetermined crankshaft rotational angles placed, respectively, during combustion and substantially at the end of combustion. The motored pressures may be derived from the combustion chamber pressure sensed at a predetermined crankshaft angle prior to combustion. If combustion chamber pressure is sensed at yet another crankshaft rotational angle prior to combustion, the pressure sensor for the combustion chamber need not sense absolute pressure, as long as it is linear and does not change its gain or bias significantly during one engine cycle. Thus four sampled combustion chamber pressure values, and only three if an absolute pressure sensor is used, suffice to provide ignition timing control. In addition, control of dilution and, in some cases, fuel charge matching are possible from the same plus, at most, one additional sampled combustion chamber pressure. This invention is further described, by way of example, with reference to the following description of a preferred embodiment, and the accompanying drawings, in which:- Figure l is a schematic diagram of an engine with a control according to this invention; Figure 2 is a graphical representation of combustion pressure and motored pressure as a function of crankshaft rotational position (crank angle) during combustion in the engine of Figure l for several possible combustion durations; Figure 3 is a graphical representation of the integral A2 of equation (l6) below as a function of crankshaft rotational position during combustion in the engine of Figure l for selected combustion initiation timings; Figure 4 is a graphical representation of the ratio of combustion chamber pressure to motored pressure after completion of combustion as a function of relative dilution of the fuel charge in the engine of Figure l; Figure 5 is a graphical representation of the ratio of combustion chamber pressure to motored pressure as a function of crankshaft rotational position during combustion in the engine of Figure l for an average MBT cycle and a late burn cycle; and Figure 6 is a graphical representation of a large number of combustion events on a plot of final pressure ratio vs. fractional pressure ratio at l0 degrees atdc (defined below). Mathematical Background First the mathematical equations upon which this invention is based will be derived. The goal of this derivation is to produce a first order approximation based on combustion chamber pressure measurement which is computationally fast enough to be usable in engine operation. The contents of the combustion chamber of an engine are assumed to behave as an ideal gas with constant properties. The derivation starts with the First Law of Thermodynamics for the mixture, with terms defined at the end of the specification before the claims. dU = dQ - dW (l) The work done by the piston may be expressed. dW = P dV (2) Other relations for an ideal gas are listed. dU = m cv dT (3) P V = m R T (4) cp - cv = R (5) cp/cv = Gm (6) The heat addition may be viewed as a balance between energy supplied to the combustion chamber contents by the combustion process and heat exchange with the combustion chamber walls. dQ = hc d(mb) - d(qw) (7) These relations may be combined to yield a differential equation in terms of the combustion chamber pressure and volume and the net heat addition. dP + (Gm P/V) dV = (R/cv V) [hc d(mb) - d(qw)] (8) When the differentials are expressed as derivatives with respect to crankshaft rotational position or angle, this equation represents a first order linear differential equation with the following solution. P(ϑ) = Pmot(ϑ) [l.0 + (cv To m)⁻¹ INTl] (9) In the preceding equation, INTl is the integral from ϑspk to ϑ of the following expression. [hc dmb(ϑ)/dϑ - dqw(ϑ)/dϑ][V(ϑ)/Vo] (Gm-1) dϑ By defining a term for the total heat of combustion, hc mf, and factoring it from within the integral, Equation 9 may be rewritten. P(ϑ) = Pmot(ϑ) [l.0 + (hc mf)/­ (cv To m)]INT2 (l0) In the preceding equation INT2 is the integral from ϑspk to ϑ of the following expression. [dmb′(ϑ)/dϑ - dqw′(ϑ)/dϑ][V(ϑ)/Vo] (Gm-1) dϑ In integral INT2, the following terms apply. mb′(ϑ) = mb(ϑ)/mf (ll) gw′(ϑ) = qw(ϑ)/(hc mf) (l2) The term in equation ll represents the mass fraction of fuel burned at any crank angle, while the term in equation l2 represents the wall heat transfer as a fraction of total fuel energy. Recognizing that the mass of the combustion chamber contents is the sum of the intake air, intake fuel and residual gases (the gases remaining in the combustion chamber from the previous cycle and recirculated exhaust gases), the following equation applies. m = ma + mr + mf (l3) This allows equation ll to be rewritten. P(ϑ) = Pmot(ϑ) [l.0 + Al A2(ϑ)] (l4) Al = (hc/cv To)[mf/(ma + mr + mf)] (l5) A2 = INT2 (l6) The factor Al depends primarily on the mass of fuel relative to the mass of diluent, decreasing with greater dilution. The term A2 depends on the time for combustion and the timing of combustion in the cycle along with the wall-heat-loss rate. Closed Loop Control of Ignition Timing Losses due to less than complete expansion of the combustion chamber charge increase as combustion is shifted away from being centred at the point of minimum combustion chamber volume (tdc). Heat transfer losses, however, decrease as combustion timing is retarded. Thus the ignition timing for best efficiency (MBT) is that timing that places the combustion event in the cycle to minimize the sum of these two losses. This generally corresponds to placing the middle of the combustion event in the range of 8 to l5 degrees atdc, depending on the distribution of heat release and factors such as engine speed and compression ratio which affect heat transfer. The optimization of ignition timing in feedback control requires that we infer from the combustion pressure trace where the combustion event is placed in the cycle. From equation l4, A2(ϑ) is the term dependent on the timing of the combustion event. It is plotted as a function of crankshaft rotational angle (crank angle) in Figure 3 as calculated for a compression ratio of 8 and with the wall heat loss neglected. The duration of combustion is specified l20 degrees. The mass burn rate is assumed to increase linearly with crankshaft rotational angle up to the 50% mass burned point and then to decrease linearly to the end of combustion, which approximates empirically observed behaviour. Curves 40-45 correspond to 50% mass burned fraction timings (CA50%) of -l0, 0, l0, 20, 30 and 40 degrees atdc, respectively. According to the relation defined in equation l4, for the same combustion chamber contents (Al constant) and the same motored pressure, differences in the value of A2 are related to differences in the fired pressure. The most retarded burn (CA50% = 40 degrees atdc) yields a higher value of A2 at l00 degrees atdc than more advanced burns. This occurs because there is less energy removed in expansion work relative to the more advanced cycles, which results in a higher internal energy and thus a higher pressure. It has been found that the location of the 50% mass burned fraction for best efficiency generally occurs in vicinity of l0 degrees atdc for conventional spark ignition engine compression ratios and for moderate to high engine speeds. Higher compression ratios and lower engine speeds favour somewhat later burning for MBT. For the purposes of this mathematical development, the l0 degree placement of CA50% for MBT will be assumed, although it will turn out that constancy of this figure with varied operating conditions is not a requirement for the invention. From Figure 3, if the value of the integral A2(ϑ) is sampled at l0 degrees atdc and divided by the value at a crank angle when combustion is over, one observes that the cycle with the CA50% at l0 degrees yields a value of this ratio of approximately 0.5. Cycles which were more advanced would yield higher values, while cycles which were retarded would yield lower values. The inclusion of heat transfer does little to modify this result. This behaviour forms the basis for controlling the ignition timing by pressure ratio management. The fractional excursion of A2 between its value before combustion and the value when combustion is complete approximates the fraction of fuel burned at any crank angle. To develop this example further, it is assumed that combustion is completed before 90 degrees atdc. From equation l4 the ratio A2(l0)/A2(90), defined as B below, can be computed. B = ([P(l0)/Pmot(l0)] - l)/([P(90)/­ Pmot(90)] - l) (l7) The determination of B requires values of the combustion chamber pressure at four crankshaft rotational angles, including two fired angles which can be sampled directly - namely P(l0) and P(90) - and two motored values, which must be approximated. It is a fair assumption to assume that the pressure at 90 degrees atdc, had combustion not occurred (Pmot(90)), is close to the pressure at 90 degrees btdc. This pressure can also be sampled directly. From the pressure measured at 90 degrees btdc, the motored pressure at l0 degrees atdc can be estimated assuming the compression process is polytropic. Pmot(l0) = P(-90)[V(-90)/V(l0)] n (l8) The polytropic exponent n is nominally l.3; and the volumes at these two crankshaft rotational angles are known for a particular engine geometry. The value of B is therefore calculable from pressures that can be measured with appropriate transducers. Generally the sensing of absolute pressure requires a transducer of known gain and bias. This can lead to impractical costs which are possible to circumvent. If a pressure sensing device produces a voltage which is linearly related to the pressure, the following relation applies. E(ϑ) = G P(ϑ) + Ebias (l9) The gain G of this transducer is assumed constant for a given cycle but allowed to vary on a longer time scale. A similar assumption is applied to the bias voltage, which is calculable from voltages sampled during the compression stroke according to the following expressions, in which ϑr is some crankshaft rotational angle in the compression stroke prior to combustion. Ebias = [E(-90) - K2 E(ϑr)]/(l.0 - K2) (20) K2 = [V(ϑr)/V(-90)] n (2l) These equations lead to the following expression for B. B = Kl [(E(l0) - Ebias) - (E(-90) - Ebias)/Kl]/[E(90) - E(-90)] (22) Kl = [V(l0)/V(90)] n (23) Ignition timing can be controlled by sampling the voltage from a pressure sensor at the four specified crankshaft rotational angles and computing B from equations 20-23. The ignition timing is then adjusted from the current value either in proportion to the error signals generated by comparing B to the targeted value of mass fraction burned, Bo, or in constant increments depending on the sign of the error: if (B - Bo) > 0 then retard combustion; if (B - Bo) < 0 then advance combustion. Because combustion events vary somewhat stochastically on a cycle to cycle basis, averaging the value of B from several consecutive cycles is probably desirable for stability. Since a signal directly related to the combustion pressure is available, knock may be detected using proper filtering of the signal for knock control. The preferred embodiment would be one wherein each combustion chamber is controlled separately. Although the polytropic exponent may be assumed as previously stated, it may also be calculated by sampling one more point along the compression curve. Certain extreme operating conditions might create problems with the suggested approach to ignition timing control. For conditions of high dilution (idle, high EGR or very lean), it is possible that advancing the ignition could yield cycles that burn later than any for a more retarded timing. Excessive advance under these conditions could produce misfire. Ignition timing under these conditions is best supplemented with dilution control, which may also be accomplished with the pressure ratio management techniques, as described below. Closed Loop Control of Dilution In equations l4-l6, the total energy liberated in the combustion process is separated from its temporal behaviour. The total energy is incorporated in the factor Al with its dependence on the amount of fuel in the combustion chamber, whereas A2 is affected by the timing and rate of energy release. An increase in the amount of diluent decreases the value of Al. Increased dilution generally reduces flame speed, thereby increasing the time for combustion and altering the integration for computing A2. The influence of combustion duration on the value of the term A2 is examined next. Figure 2 shows a plot of motored pressure (curve 50) and the pressure calculated from equation l4, neglecting the heat loss term, for combustion durations of 40, 80 and l20 degrees (curves 5l, 52 and 53) and a constant value of Al. The combustion events are timed to place the CA50% at l0 degrees atdc. The pressure traces are observed to be significantly different. They have certain things in common, however. First, at l0 degrees atdc, the pressures are nearly identical. Furthermore, after combustion is completed, the pressures again are nearly identical. This implies that, for MBT timing, the pressure at any point in a cycle depends primarily on the amount of fuel burned to that point and not significantly on the previous history in the cycle. It does depend on the ratio of fuel to diluent, as dictated by term Al. The ratio of fired to motored pressure after combustion has been completed is a measure of the fuel to diluent ratio. The level of dilution depends on the excess air and/or EGR entering the combustion chamber during intake and on the amount of residual gas remaining in the combustion chamber from the previous cycle, which is a function of load, valve timing and exhaust pressure. In general, it is known that increasing dilution increases efficiency and reduces oxides of nitrogen (NOx) emissions up to the point where combustion rates slow sufficiently to affect efficiency adversely. In Figure 4, the ratio of fired pressure to the motored pressure evaluated at the crankshaft rotational angle at which combustion is completed is plotted as a function of diluent to fuel ratio. Increasing dilution is observed to decrease the pressure ratio, approaching l.0 for high dilution. In an actual engine, poor combustion limits the allowable level of dilution and the corresponding pressure ratio to values significantly higher than l.0. Poor combustion manifests itself in several forms depending on the design of the engine. The various options for controlling the level of dilution with pressure ratio management can be divided into three primary categories. The first is based on controlling the amount of diluent to a specified value as a function of speed and load. The second option is to use diluent level to adjust the burn rate or combustion time to a specified value. The third is based on using the variability in the pressure ratio among cycles as a dilution level criterion. Each of these options will be discussed. In the testing of a given baseline engine, it is possible to define a dilution limit based on one of a number of criteria. Whatever the criterion, this dilution level provides a target value for the ratio P(90)/Pmot(90). Since signals for P(90) and Pmot(90) are available from the scheme for ignition timing control, comparison of their ratio to the target value indicates the change required in EGR level. Such control must account for the response time of the EGR valve and the delay associated with the transit time of the flow of exhaust gases from the EGR valve to the combustion chamber. This can be accomplished by averaging a number of cycles to evaluate the dilution level. Depending on the distribution of EGR among combustion chambers, it might be desirable to adjust the dilution to yield the targeted ratio for either a specific combustion chamber, probably the most dilute, or the average of all combustion chambers. The targeted level of dilution might be a function of both engine speed and load. The dependence of residual gas fraction on load is compensated in this approach. However, turbulence, and therefore flame speed, are lower at light load, which suggests the need for less diluent and a higher pressure ratio target than used at higher loads. Engine speed has an influence on the pressure ratio target because of the longer time available for heat transfer at low engine speed. For the same total dilution, the pressure ratio will be lower at low engine speed because of this effect, which suggests that the target ratio should decrease as speed decreases. The algorithm for controlling dilution with this option is developed from equation l4, in which the pressure ratio after combustion can be related to the dilution level. Al A2 = [P(90)/Pmot(90)] - l.0 (24) In this expression A2 has assumed its final value because combustion is complete. The relations specified in equations 2l-23 allow evaluation of P(90)/Pmot(90) in terms of sampled voltage. [P(90)/Pmot(90)] - l.0 = [(K2 - l)/K2] ([E(90) - E(-90)]/E(-90) - E(ϑr)]) (25) The relative dilution is thus calculable from three sampled voltages. The voltage relation may be defined as follows. C = [E(90) - E(-90)]/[E(-90) - E(ϑr)] (26) The targeted value of Al A2 may be determined by engine testing. This value will be specified as Co, which might be a function of speed and load. The condition for control of EGR is expressed in the following conditional statements. C > Co[K2/(K2-l)] increase EGR C < Co[K2/(K2-l)] decrease EGR A possible problem condition for the ignition timing control comes when a cycle is both over dilute and retarded. For such a condition it is possible for the combustion timing control to satisfy the control criterion B=0.5 at l0 degrees atdc but not provide true MBT because combustion is retarded sufficiently to be incomplete at the second sampling point at 90 degrees atdc. Even if this is the case, the dilution algorithm will properly direct a decrease in the level of dilution. The ignition timing will eventually stabilize to a true MBT as the level of dilution is decreased. This situation would most commonly be encountered after a rapid reduction in throttle, for which a quick return to MBT is probably not important. Although the algorithm as described above uses the average final pressure ratio over a plurality of engine cycles, it may also be used for a single cycle for better transient performance with a slight modification. Figure 6 shows the final pressure ratio and fractional pressure ratio at l0 degrees atdc for a large number of combustion events at three different air/fuel ratios: l8/l (curve 70), 2l/l (curve 7l) and 24/l (curve 72). Air/fuel ratio was varied by changing the airflow with fuel flow held constant. For each air/fuel ratio the separate symbol types correspond to cycles from three different ignition timings. The fractional increase in the pressure ratio at l0 degrees atdc is the timing parameter previously designated as B. For the engine speed corresponding to the data shown, the best efficiency spark timing produced a cycle average value of 0.6 for B. Cycles with higher values are overadvanced, while cycles with lower values are retarded. As Figure 6 shows, overadvanced cycles produce lower final pressure ratios than more retarded cycles with the same mixture strength. Two factors are responsible for this behaviour. The more advanced a cycle is, the greater are the heat losses to the chamber walls because of higher temperatures and pressures during the flame process. This results in a lower pressure when combustion is complete. As the combustion process is moved away from being centred at tdc there is a reduction in the amount of work removed from the combustion gases, which increases the pressure after combustion. This, combined with the reduced heat loss, increases the final pressure ratio as the timing of a cycle is more retarded. The influence of timing on the final pressure ratio must be considered when estimating relative dilution. This is accomplished by estimating the value of the pressure ratio had the cycle occurred with an optimum value of the timing parameter. Since the slopes of the lines in Figure 6 are essentially equal, the sensitivity of the final pressure ratio to the timing parameter is independent of the dilution level. The final pressure ratio for optimum timing can be expressed in terms of the slope of the line, d(PRF)/dB, the final pressure ratio for the cycle, PRFi, and the difference between the timing parameter of the cycle, Bi, and the value for optimum timing, Bo. PRF(Bo) = PRFi + (d(PRF)/dB)(Bi - Bo) (27) Subtracting the value of l.0 from PRF(Bo) produces the parameter proportional to the relative dilution, mf/(ma + mr + mf). With the single cycle approach, it is apparent that the combustion events at each air/fuel ratio cluster about an essentially straight line and that the lines appear to be parallel. The individual combustion events may vary greatly in fractional pressure ratio at l0 degrees atdc, although the average is being maintained at a target of 0.6, for example. For a given air/fuel ratio, therefore, the observed final pressure ratio of each individual combustion event may be converted to the equivalent final pressure ratio at 0.6 ten degree atdc pressure ratio as long as the fractional pressure ratio at l0 degrees atdc is known for the individual combustion event, which it is. As an example of this conversion, assume a calculated final pressure ratio of 3.l and a ten degree atdc pressure ratio of 0.3 for a single combustion event. Referring to Figure 6, this corresponds to a point on curve 72. Moving down curve 72 until a ten degree atdc pressure ratio of 0.6 is reached, one may read out a corrected final pressure ratio of 3.0. If the point had not been on one of curves 70-72, another curve could have been drawn parallel to them through the point. With the single cycle approach, the control may be used in a much faster mode with individual combustion events rather than being based solely on the average of many combustion events, with improved transient performance. As described to this point, this option for dilution control is straightforward. However, the pressure ratio is being used to control the amount of diluent to a specified level, or percentage of the intake level, based on the performance of a baseline engine. This limit is not necessarily the same for the engines eventually manufactured or for aged engines. It would be preferred that the dilution level be related to a more fundamental characteristic of combustion at the dilute limit, which the pressure signal can sense. One possibility is to update the stored target final pressure ratio during engine operation; and an improvement using this approach will be described along with the third option of dilution control. Before this, however, a second option of dilution control will be described. In the previous option of pressure dilution level control by pressure ratio management, the magnitude of the final pressure ratio was used to control the dilution level. The time for combustion also depends on dilution level. Alternative dilution control algorithms can be formulated based on this dependence. Since the combustion timing is to be controlled by placing the 50% mass fraction burned at a given crankshaft rotational angle, it is possible to adjust the mass fraction burned at a different crankshaft rotational angle by varying the level of dilution. This is equivalent to specifying the time to burn a given mass fraction or a combustion duration. The mass burn fraction profile is thus tailored by using the level of dilution. Presume that, for optimum dilution, combustion is 90% complete at a given crankshaft rotational angle: 60 degrees atdc, for example. By sampling the voltage at this crankshaft rotational angle and computing a parameter equivalent to parameter B in the combustion timing application, dilution should be increased if the value is greater than 0.9 or decreased if it is less than 0.9. This application requires the sampling of two more voltages, E(60) and E(-60), the latter being assumed to be equivalent to Emot(60). If E(-60) is also selected as E(ϑr) in the combustion timing application, this implementation requires the sampling of only one extra voltage, E(60). In this form the mass fraction burned at 60 degrees atdc can be specified by D, as follows. D = [Pmot(90)/Pmot(60)][P(60) - Pmot(60)]/[P(90) - Pmot(90)] (28) This may be expressed in terms of the sampled voltages. D = [V(60)/V(90)] n [E(60) - E(-60)]/­ [E(90) - E(-90)] (29) This approach has an attractive feature if the value of D is chosen as l.0. In those cases where extreme dilution and retard might produce a burn which satisfies the combustion timing criterion, the value of D will be less than l.0 whenever burning is not complete at 60 degrees atdc. This will yield a request for less dilution and stabilize the ignition timing. In specifying pressure ratio management control based on sampling at a crankshaft rotational angle after l0 degrees atdc, characteristics of combustion late in the burn are emphasized. Similar control options exist to emphasize the time for the first portion of the burn by specifying the mass fraction burned at a crank angle before l0 degrees atdc. Alternatively, the time from combustion initiation to the 50% mass burned fraction point could be specified by adjusting dilution level until MBT timing occurs at a targeted advance. The previous two options of dilution control are based on the average behaviour of collections of cycles. As previously mentioned, it has been observed that, as the dilution limit is approached, a greater number of abnormally late burning cycles are produced. This behaviour suggest a third possible criterion for the control of dilution. Figure 5 shows typical traces of the pressure ratio as a function of crankshaft rotational angle. Curve 60 corresponds to a nominal cycle with the burn placed optimally so that the 50% mass fraction burned point, and therefore 50% of the final pressure ratio, occurs at l0 degrees atdc. Curve 6l corresponds to a late burn cycle. A late burning cycle will produce a low value of mass fraction burned at the l0 degree atdc sampling point, whereas the average of the collection of cycles will yield a value of 0.5 because the combustion timing is adjusted to yield this value. A current cycle average of B at the l0 degree atdc point could be used as a late burn indicator. If late burns are detected over a sampling period of some number of cycles, dilution should be decreased. If no late burns are sensed, the dilution should be increased. This technique places a bound on the extremes of the variation of mass burned at a given point in the cycle. Similarly, the variability in B could be computed and adjusted to a targeted value using dilution level. The late burn technique also shows the possibility of being used to update the target final pressure ratio of the first option in dilution control. This final pressure ratio was described earlier as being determined for a baseline engine and stored as a reference in a read only memory ROM, perhaps in a table as a function of one or more other engine parameters such as speed or load. However, the normal differences among engines and changes due to engine age and wear may be compensated if some parameter exists which allows correction of the stored reference during engine operation. The late burn technique determines the lean limit for the engine, since it identifies the operating condition on the verge of misfire, where efficiency begins to suffer. The first option may be run on a continuous basis, with dilution level controlled to maintain a predetermined final pressure ratio. However, the number of late burns may be counted over a comparatively long series of cycles with the result used to adjust the target final pressure ratio. Thus the stored reference could be regularly updated in engine operation to produce engine operation at the lean limit. It is usually desirable to use the various dilution control options discussed above in combination with the combustion pressure control of ignition timing as previously described, since many of the equations are common and efficiency in calculation and computer time is promoted. However, it is also possible to use any of the dilution control methods by itself, with a slightly different method of ignition timing control. When a dilution level is computed according to any of the methods described, ignition timing is derived from a stored table in a read only memory ROM or random access memory RAM addressed by engine speed and dilution level. Since the dilution level input depends on combustion pressure, this factor also affects ignition timing through the table lookup. However, this method may have some advantage in that it is thought that it may be used without the averaging (lag filter function) normally required for ignition timing control as described earlier and thus enable the ignition timing to respond more quickly to transients while still allowing dilution flow control by combustion chamber pressure. Control of Fuel Flow and the Lean Burn Engine EGR to control NOx emissions has been widely implemented because difficulties were encountered when the alternative of lean operation was attempted. Generally, for the same masses of EGR and excess air, flame temperatures are somewhat higher and oxygen more plentiful with excess air, which favours higher NOx production. With sufficient dilution with excess air (air/fuel ratios above l8/l), flame temperatures are sufficiently low to inhibit NOx. At this level of dilution, combustion rates are slowed. Because most engines do not distribute fuel equally to all combustion chambers, some will operate leaner than the overall mixture ratio, some richer. Richer mixture combustion chambers will be operating over advanced and producing higher NOx than would be indicated by the overall air/fuel ratio, while leaner mixture combustion chambers will be retarded and yield less power than if they were timed optimally. The overall dilution level of the engine is additionally limited by poor combustion in the combustion chamber with the leanest mixture. For these reasons it is desirable for the optimal implementation of a lean calibrated engine to adjust combustion timing and air/fuel ratio on an individual combustion chamber basis. The potential benefits arise from NOx control without the reducing catalyst and difficult to control EGR system. It will be assumed that fuel is supplied to each combustion chamber of the engine through a separate fuel injector. The amount of fuel injected is presumed to be adjusted by specifying the duration of an electrical pulse to the injector driver. Because of manufacturing tolerances and differences in temperature, the amount of fuel delivered to each combustion chamber can vary. It is assumed that the intake manifold and intake ports are designed to produce minimal differences in airflow among the combustion chambers. This being the case, differences among the combustion chambers in air/fuel ratio can be attributed to differences in fuel flow rate. Once it has been specified that dilution is controlled by adjusting excess air, it is reasonable to make fuel rate the variable that responds to the driver's load request rather than air as regulated by the traditional intake throttle. In the proposed scheme, load demand is specified by the accelerator pedal, which directly adjusts the mean pulse duration of the fuel injectors. The air throttle is adjusted to correspond to the mean injector pulse but trimmed to adjust the diluent flow rate. As shown in equations l4 and l5, the pressure ratio after combustion is directly related to the ratio of fuel mass to total mass in the combustion chamber. Differences in individual combustion chamber values of the pressure ratio can be attributed to fueling differences. In terms of sampled voltages, the term proportional to the fuel rate is represented by equation 26, which is repeated here with the letter "i" to indicate that the value is computed for each combustion chamber. Ci = [Ei(90) - Ei(-90)]/­ [Ei(-90) -Ei(ϑr)] (30) It is assumed that the total fuel rate to the engine will remain approximately constant if the average pulse duration remains constant. The pulse duration for each fuel injector can then be adjusted by multiplying the mean pulse duration by the ratio of the average value of C among all combustion chambers to the individual combustion chamber value. Ti = Tmean (Cmean/Ci) (3l) The airflow to the engine is controlled to maintain the dilution limit as sensed by one of the three techniques suggested previously. Naturally, this method of fuel flow control fits nicely with the ignition timing control based on combustion chamber pressure as described earlier, since the pressure sensor and some of the equations are in common; however, this fuel control method has utility of its own and may be used with other ignition timing systems and with engines using dilution gas rather than lean air/fuel ratios for dilution. In operation the engine behaves as follows. At idle the residual fraction is high, resulting in an air/fuel ratio approximating stoichiometric. As the load demand is increased, approximately constant dilution is maintained and air/fuel ratio thus increases. Eventually a point is reached in the vicinity of 60% load where the air throttle is wide open and the engine operates at its leanest allowable air/fuel ratio (nominally 22/l). At this point the airflow is maximum and a request for more power results in a higher fuel flow and a richer mixture. At approximately 75% load the air/fuel ratio has been decreased to about l8/l. Further enrichment may be achieved without greatly increasing NOx production by "switching" the air/fuel ratio to the other side of stoichiometry one combustion chamber at a time and thus avoiding the region of high NOx production. A lean burn engine control based on pressure ratio management should operate with higher driveability than an engine calibrated with high levels of dilution. In the latter engine, increased load demand is initiated by a change in the air throttle. The airflow thus leads the fuel flow, which temporarily produces even leaner cycles and a possible "sag" in engine response. If the control structure is reversed, as in the lean implementation suggested here, fuel flow will lead the airflow and produce richer cycles in the initial stages of acceleration to provide more rapid combustion. This approach might result in a greater tendency to knock; but actual knock may be contained with a knock control based on the same pressure signals. Hardware and Software Description Apparatus capable of controlling an internal combustion engine according to the principles discussed above is shown schematically in Figure l. Engine l0 is a spark ignited internal combustion engine having four combustion chambers or cylinders ll, an intake manifold l2 with a throttle l3 , an exhaust manifold l5 with an oxidizing catalyst l6 and a rotating crankshaft, not shown, through which power is transferred from combustion chambers ll to the vehicle drivetrain. Each combustion chamber ll is provided with an intake valve l7 communicating with intake manifold l2, an exhaust valve l8 communicating with exhaust manifold l5, a spark plug 20, a fuel injector 2l in intake manifold l2 adjacent intake valve l7 and a pressure sensor 22 for each combustion chamber. Computer 23 is a digital computer including the standard elements of central processing unit CPU, random access memory RAM, read only memory ROM, input/output circuitry and clock circuitry. Appropriate inputs and outputs are indicated for sending signals from computer 23 to fuel injectors 2l (Fl-F4) and spark plugs 20 (Sl-S4) and to computer 23 from pressure sensors 22 (Pl-P4). The computer 23 may be a variant of those used for computerized engine control on prior art vehicles based on the Motorola 68000 (R) series microprocessor. Spark plugs 20 and fuel injectors 2l are of the type normally used in the prior art. Pressure sensors 22 may be any known type of pressure sensor adapted or adaptable to sensing the pressure inside combustion chambers ll, providing they generate a linear voltage output signal with a gain and bias that do not change significantly during a cycle of engine l0. One pressure sensor 22 which might be appropriate is a piezoelectric deflection type of transducer which includes a finger projecting through the water jacket of the combusion chamber ll to sense pressure caused deflections of the wall of the combustion chamber. Another pressure sensor 22 which might be appropriate is a piezoelectric sensor attached directly to the wall of the combustion chamber ll. Pressure sensors 22 with membranes of piezoelectric or other materials exposed directly to the combustion chamber interior or combined with the spark plug 20 are further possibilities. There is another promising pressure sensor in which a fibre optic conductor is routed through a tortuous course in the spark plug 20 body to be deflected by combustion chamber pressure variations in such a way as to affect the light transmissibility thereof. Further research is required to determine the best pressure sensor. An ignition circuit 25 is effective when signaled by computer 23 to generate a high voltage pulse to activate a spark plug 20 and to direct the pulse to the proper spark plug. Ignition circuit 25 and spark plugs 20 define ignition apparatus. It may include the standard distributor and take any appropriate form from the prior art. The invention is not concerned with how the pulses are generated or directed to the spark plugs 20 but with when the signal pulses are generated within computer 23 for ignition timing. Signals indicative of crankshaft rotational position are generated in crankshaft rotational position sensor 26 which may include a toothed wheel 27 and a Hall effect, magnetic or optic pickup 28. Any reasonably accurate rotational position sensing apparatus from the prior art may be used; however, the more accurate the apparatus and its mounting on engine l0, the more the benefit gained from this invention. An EGR valve 30 controls the flow of diluent gases from oxidizing catalyst l6 back to intake manifold l2 and includes a positioning motor controlled by signals from computer 23 over the line marked EGR. Pressure sensors 22 and computer 23 define pressure responsive means responsive to combustion chamber pressure to generate signals for a combustion event indicating combustion pressure and motored pressure at a first predetermined crankshaft rotational angle during combustion and at a second predetermined crankshaft rotational angle at substantially complete combustion. Computer 23 also defines computing means effective to compute, from the signals, the first and second pressure ratios of combustion pressure to motored pressure at the first and second crankshaft rotational angles respectively. Computer 23, ignition circuit 25, crankshaft rotational position sensor 26, toothed wheel 27 and pick-up 28 define adjusting means responsive to the computing means to adjust the ignition timing to maintain a predetermined ratio between the first and second pressure ratios each decreased by one, whereby combustion timing is controlled in closed loop to an MBT value. The implementation of the algorithms already discussed should be clear to any programmer skilled in the art of engine control. The signal pulses from crankshaft rotational position sensor 26 may be used to generate interrupts in the main operating program of computer 23. If these interrupts occur with a frequency as high as the required accuracy of crankshaft angle determination, the interrupts may simply trigger a counter or increment a memory location used to keep track of the angle. If, as appears more likely, some interpolation between interrupts is required, higher frequency pulses from a real time clock may be used to measure the time between interrupts with the time being computationally subdivided to generate the counter triggering pulses. Numbers representing the required crankshaft rotational angles for pressure sensing are stored in memory and continually compared to the counter for the generation of pressure sensing signals at the proper time. Each such pressure sensing signal causes the digitized pressure voltage signal to be captured and stored in the appropriate memory location or register. The pressure sensing takes place for each engine cycle during the crankshaft rotational angle between 90 degrees btdc and 90 degrees atdc. Following the last required sensed pressure value, the computer 23 is signaled to commence computation of the required parameters, comparison with appropriate stored or computed references and output of control signals to output circuitry effective to generate ignition signals, fuel injection signals and control signals for the EGR valve 30. The particulars of the program will depend on choice of apparatus and of algorithms. An example of the control algorithm will be given. Assume that the combustion chamber pressure is sampled at 90 and 60 degrees btdc (-90, -60) and at l0 and 90 degrees atdc. For ignition timing, equation (20) is first evaluated, with E(ϑr) being the reading at -60 degrees. Next, the resulting value of Ebias is used in the solution of equation (22) for the parameter B. In these equations, Kl and K2 are predetermined from the known combustion chamber geometry from equations (23) and (2l), respectively, and stored in ROM for use in the calculations. The parameter B is averaged over a plurality of engine cycles with a first order lag function and, with each new update the average is compared with a predetermined constant from a table or single location in memory to determine whether ignition timing is advanced or retarded. Continuing with the first dilution control option, parameter C is calculated from equation (26) and either averaged over a plurality of engine cycles or converted as described earlier to an equivalent value for MBT timing and compared with a predetermined target value. This target value may be a reference stored in ROM or an updated reference based on the late burn criterion of the third option as described earlier and below. EGR is increased or decreased depending on the sign of the comparison. If the second dilution control option is chosen, the parameter D, rather than C, is calculated from equation (29). This, of course, requires one more sampled sensor voltage, at 60 degrees atdc. The value of D is compared to a predetermined target value stored in ROM. If the third dilution control option is chosen by itself, neither C nor D need be calculated. An average value of B is maintained as a reference to which the individual values of B are compared. The individual values of B which fall short of the average value by an amount greater than a predetermined, stored reference value are summed during a sample period of predetermined duration and compared with another predetermined, stored reference value. The result of this comparison determines whether EGR is increased or decreased, with decrease resulting from too many low values of B signaling too many late burns, as previously described. If the third option is used to update the target value of C for the first option, the result of the comparison may be used to adjust the targeted value of C in one direction or the other. If the engine uses one of the dilution control options but does not use the ignition timing option, dilution is determined first by any of the algorithms described above and the ignition timing is then derived from a lookup table in ROM or RAM addressed by dilution level and engine speed, the latter being determined in a manner well known in the art from the pulses generated by crank angle signal apparatus 26. It is expected that individual values of dilution level will be used without averaging for good transient response. If the engine is a lean burn engine with individually controlled fuel injectors and balanced intake airflow for each combustion chamber, fuel charge balancing may also be included. The parameter Ci is calculated from equation (30) for each combustion chamber. From these values the average or mean value of C, Cmean, may be calculated and the fraction Cmean/Ci determined for each combustion chamber. The injector pulse width Ti for each fuel injector, determined by whatever fuel control algorithm is selected for a desired average fuel flow rate, is multiplied by the fraction Cmean/Ci for the fuel injector to produce a balanced fuel charge. Naturally, all the normal techniques of feedback control may be used as desired. For example, the value of a parameter such as B, C or D, which is to be compared to a stored reference, may be averaged with a lag filter function; a dead band may be established; and the correction value may be of constant absolute value or some variant of proportional control, all at the choice of the designer to fit the particular engine. When calculations are complete, the output correction values for ignition timing and EGR, if the engine is not a lean burn engine, are provided to the normal control routines or circuits for controlling these parameters. The correction factor for fuel injection timing is provided to the injection timing routine. The actual control of the ignition, EGR valve and fuel injectors from this point is well known in the prior art. In the present invention, the pressure sensors 22 provide a signal which is sampled at several crankshaft rotational angles to obtain or derive the actual combustion pressure (5l,52 or 53 in Figure 2) and the motored pressure (50 in Figure 2; the pressure with no ignition) at two crankshaft rotational angles; a first angle (for example, l0 degrees atdc) during combustion at which a desired mass fraction of fuel (for example, 50%) is to be burnt, and a second angle (for example, 90 degrees atdc) after combustion at which l00% of the fuel is burnt. The ratio of combustion pressure to motored pressure is calculated for each angle. The pressure ratio minus one at the first should be the given fraction (50%) of the pressure ratio minus one at the second, since the increase in pressure ratio varies with the mass fraction of fuel burnt. Ignition timing is varied in closed loop to maintain the desired ratio. Attention is drawn to our other applications, nos. ........(MJD/l836) and ........(MJD/l837) filed the same day as this application. Nomenclature Table atdc - after top dead centre btdc - before top dead centre B - timing parameter, fractional pressure ratio at l0 degrees atdc Bo - timing parameter corresponding to optimum timing Ca50% - crank angle of 50% mass burn (atdc) Ci - individual combustion chamber fuel/diluent parameter Cmean - average of individual Ci's cp - specific heat at constant pressure (J/kg-K) cv - specific heat at constant volume (J/kg-K) EGR - exhaust gas recirculation E - voltage (V) Ebias - bias voltage (V) Ei - transducer voltage for individual combustion chamber G - transducer gain (V/kPa) Gm - ratio of specific heats (cp/cv) hc - heat of combustion (J/kg) m - mass (kg) ma - mass of intake air (kg) mb - mass of fuel burned (kg) MBT - minimum ignition (spark) advance for best torque (btdc) mf - mass of fuel (kg) mr - mass of residual gas (kg) n - polytropic exponent P - pressure (kPa) Pmot - motored pressure, i.e., pressure in combustion chamber with no combustion and with crankshaft turned by outside agent (kPa) PRF - Final Pressure Ratio, pressure ratio after the completion of combustion PRF(Bo) - Final Pressure Ratio with optimum timing qw - wall heat transfer (J) Q - heat addition (J) R - gas constant (J/kg-K) T - temperature Ti - fuel pulse for individual fuel injector Tmean - average fuel pulse among all injectors To - temperature at intake valve closing (K) tdc - top dead centre U - internal energy (J) V - volume (m³) Vo - volume at intake valve closing (m³) ϑ - crankshaft rotational angle (tdc=0, atdc>0, btdc<0) ϑr - reference crank angle during compression ϑspk - crank angle of ignition

Claims: 1. An ignition timing control for an internal combustion engine of the type having a combustion chamber(ll); a rotating crankshaft; and ignition apparatus (20,25) for the combustion chamber; the ignition timing control being characterised by pressure responsive means (22,23) responsive to combustion chamber pressure to generate signals for a combustion event indicating combustion pressure and motored pressure at a first predetermined crankshaft rotational angle during combustion and at a second predetermined crankshaft rotational angle at substantially complete combustion; computing means (23) effective to compute, from the signals, first and second pressure ratios of combustion pressure to motored pressure at the first and second predetermined crankshaft rotational angles, respectively; and adjusting means (23,28) responsive to the computing means to adjust the ignition timing to maintain a predetermined ratio between the first and second pressure ratios each decreased by one, whereby combustion timing is controlled in closed loop to an MBT value. 2. An ignition timing control according to claim l in which the second predetermined crankshaft rotational angle is 90 degrees after tdc. 3. An ignition timing control according to claim l or in claim 2 in which the signals of combustion pressure are derived directly from the sensed combustion chamber pressures at the first and second predetermined crankshaft rotational angles and the motored pressures are derived from combustion chamber pressure at a third predetermined crankshaft rotational angle prior to combustion in which combustion chamber volume is equal to the combustion chamber volume at one of the first and second predetermined crankshaft rotational angles. 4. An ignition timing control according to claim 3 in which the third predetermined crankshaft rotational angle is 90 degrees before tdc. 5. An ignition timing control according to claim 4 in which the motored pressure at 90 degrees after tdc is derived directly from the sensed combustion chamber pressure at 90 degrees before tdc and the motored pressure at the first predetermined crankshaft rotational angle is derived from the combustion chamber pressure at 90 degrees before tdc as the inverse ratio of the volumes of the combustion chamber at the two predetermined crankshaft angles raised to the power of the polytropic exponent. 6. An ignition timing control as claimed in Claim l, in which the pressure responsive means comprises a pressure sensor (23) responsive to combustion chamber pressure to generate an output voltage, the pressure sensor being linear with a gain and bias which do not significantly vary during a single engine cycle; and means (23) effective to sample the output voltage of the pressure sensor at third and fourth predetermined crankshaft rotational angles in the compression stroke prior to combustion, at the first predetermined crankshaft rotational angle during combustion and at the second predetermined crankshaft rotational angle at substantially complete combustion, the third and second predetermined crankshaft angles having identical combustion chamber volumes. 7. An ignition timing control according to claim 6 in which the ratio of the first and second pressure ratios is derived according to the following equations: B = Kl[(E3 - Ebias) - (El - Ebias)/Kl /(E4 - El) wherein Ebias = (El - K2 E2)/(l - K2), Kl = (V3/V4)n, K2 = (V2/Vl)n, n is the polytropic exponent and El, E2, E3, E4 and Vl, V2, V3, V4 are sensor output voltages and combustion chamber volumes, respectively, at the third, fourth, first and second predetermined crankshaft rotational angles. 8. An ignition timing control according to claim 6 or claim 7 in which the second predetermined crankshaft rotational angle is 90 degrees after tdc and the third predetermined crankshaft rotational angle is 90 degrees before tdc. 9. An ignition timing control according to claim 8 in which the fourth predetermined crankshaft rotational angle is 60 degrees before tdc and the first predetermined crankshaft rotational angle is l0 degrees after tdc. l0. An ignition timing control according to any one of the preceding claims, comprising a diluent gas flow control apparatus (30) and gas adjusting means (23,30) effective to adjust the diluent gas flow to maintain the second pressure ratio at a predetermined value, whereby combustion timing and diluent flow are controlled optimally in closed loop. 11. An ignition timing control according to Claim l0, in which the pressure responsive means (22,23) generates a signal for a combustion event indicating combustion pressure and motored pressure at a fifth predetermined crankshaft rotational angle during combustion; the computing means (23) computing, from the signal, third pressure ratios of combustion pressure to motored pressure at the fifth predetermined crankshaft rotational angle; and the gas adjusting means (23,30) and adjusting diluent gas flow to maintain a predetermined ratio between the second and third pressure ratios each decreased by one. 12. An ignition timing control according to Claim l0 in which the gas adjusting means (23,30) is effective to maintain an average of consecutive ratios of the first and second pressure ratios each decreased by one and adjust diluent gas flow in response to the degree of variation of individual ones of the consecutive ratios from the average. 13. An ignition timing control according to claim l2 in which the diluent gas flow is decreased when a predetermined fraction of a number of the consecutive ratios are smaller by a predetermined amount than the average of the consecutive ratio and otherwise to increase the diluent gas flow. 14. An ignition timing control as claimed in any one of the preceding claims for an internal combustion engine (l0) having a plurality of combustion chambers (ll), air supply apparatus (l2,l3) for the combustion chambers, and fuel supply apparatus (2l) having individual control for each combustion chamber, the ignition timing control comprising means (23) effective to derive from values of the second pressure ratio an average second pressure ratio for all combustion chambers and adjust the fuel supplied to each combustion chamber in accordance with the average second pressure ratio divided by the second pressure ratio for that combustion chamber to even fuel distribution among the combustion chambers, whereby the ignition timing is controlled to MBT and the fuel flow is optimized for lean burn operation.

Inventors: MATEKUNAS, FREDERIC ANTON
Applicants: GENERAL MOTORS CORPORATION
Application year: 1986
Main classification: F02P5
IPC codes: F02P5