This invention relates to calibration techniques for dynamic systems having one or more actuators, such as engine systems.
Traditional calibration techniques for engines are based upon steady-state data collection on an operating engine, over a range of actuator settings and operating conditions. The steady state data is used with an optimization methodology to arrive at optimal or near optimal actuator settings to achieve goals for engine emissions, performance, and fuel economy.
The effectiveness of the traditional technique is high, in that optimal actuator settings can be found. However, the time required for the data collection and optimization process is too long. Considerable time is consumed by the data collection process because the engine must stabilize at each test point until steady state operation is achieved.
For situations in which true steady state is unachievable or requires too much time, transient data may be acquired, with the assumption that the persistent transient behavior of the engine will not appreciably affect the final steady state calibration. Of course, for highly transient engine operation, this assumption leads to poor calibration results.
Advanced methods for selection and use of “fractional” experimental data sets for use in the optimization process generally fall into the category of statistical design of experiments. The key to the method is in selection of experimental data points that will be most useful, while elimination of less useful data points, during acquisition, is made. Hence, the eliminated data is not available for later use. The time savings accomplished with these methods is a strong function of the numbers of experimental points not tested. A tradeoff exists in time-savings versus the ability to arrive at optimal actuator settings with a fractionated experimental data set.
The basis of the invention is a transient method for deriving steady state engine calibrations. The method is based upon non-steady state experimentally derived data, which is manipulated to arrive at steady-state optimal actuator settings.
More specifically, unlike conventional calibration methods, data acquisition is not delayed until the engine fully settles at each test point. Instead, as explained below, slow response variables are distinguished from fast ones. The engine is calibrated for the fast settling variables, and correction terms are added for slow variables.
The use of transient experimental data to derive steady state optimal actuator settings provides dramatic time-savings for engine calibration. This technique offers high value to the engine and engine-related industry by reducing the time and cost of calibrating engines.
Although the proposed calibration method does not rely on steady state data, it accounts for the transient nature of the test when optimizing the engine calibration. A key to the proposed invention is the recognition that a description of optimized actuator settings for the transient engine system will naturally include optimal calibration for the steady state system.
It is recognized that all steady state systems are a subset of a more mathematically complex transient system. In other words, a mathematical description of the dynamic system will contain a full description of steady state behavior. This is seen in mathematical descriptions of physical processes, where the solution for the current state of the system is a combination of time dependent and non-time-dependent terms.
For purposes of example, this description is in terms of a heavy duty diesel engine and typical actuators used on such an engine. However, the same concepts are applicable to other engines and their actuators, and the engines may be internal and external combustion engines. In fact, the concepts described herein could be applied to any dynamic system having actuators that are calibrated in response to system operating conditions.
(VGT) vane or other turbocharger output means Actuators associated with these parameters include fuel injectors 11, EGR valve 12, and turbocharger vane 13.
The following is a subset of the optimization variables of interest, available directly from, or derived from, measurements taken by various engine sensors:
With more than six degrees-of-freedom, one can obtain desired torque (or any other output) from the engine with many different combinations of actuator setpoints. Some of these combinations are more desirable than the others, and it is in this sense that a combination of settings is viewed as “optimal”.
Example Optimization problem
An example of an optimization problem may be stated as follows: Given an engine, find, for a specified engine speed, actuator settings to produce specified torque so as to minimize the BSFC while not exceeding specified thresholds of engine out NOx and particulate emissions, peak cylinder pressure, and VGT speed.
This is a classic constrained optimization problem, except that there are no expressions/formulae for the engine outputs as functions of the actuator settings. Such relations, however, can be “constructed” from measurements performed on the actual system.
Lumped Parameter Description of Dynamical Systems
Before proceeding, we define the “state” of the engine system. If the states are known, the system is fully defined.
More specifically, we consider a lumped parameter, state-space model of a dynamic system and how it relates to the optimization problem above. Thus,
x=f(x,u),y=g(x,u)
, where x ⊂ X ε n is the state vector and u ⊂ε m is the input vector.
Let x* denote an equilibrium point (steady state) of interest for a given input vector u. One can appeal to an implicit function theorem, and solve 0=f (x, u) for x* in terms of u at least locally around x*, say,
x*=h(u)
This solution provides a steady state value of x in terms of input u. Substituting for x in the output equation:
y=g(x*, u)=g(h(u),u)=ĝ(u).
The above equation forms the basis of “input-output” mapping for a dynamical system in steady state. It simply states that the outputs of a system in steady state depend solely upon its inputs. However, this expression gives no indication as to how long may one have to wait before a steady state is reached.
One way to address the question of how long to wait is to consider the eigensystem (eigenvalues, eigenvectors) associated with the linearized system matrix D1f|x*. This provides an insight into local behavior of the system around x*. Assume that the eigenvalues can be grouped as relatively “slow” and fast” depending on their real part.
Based on the eigenvectors, one can in turn label some states, xs, as relatively “slow” and other states, xf, as “fast”. Given this partitioning of states, one can rewrite the above as:
One way to obtain a reduced order system description is to assume that the fast states xf have settled. In that case, 0=xf−ff (x,u). Assume that one could solve this implicit equation for xf: xf=ψ(xs, u) to obtain:
xs=f(xs,u), y=ĝ(xs,u)
,where
and ĝ is similarly defined.
An important conclusion then is that in quasi-steady state where one waits just long enough for the “fast” states to settle, the system outputs can be viewed to depend solely on the slow states and the inputs. Notice that this is different from the earlier conclusion that in steady state, the outputs depend solely on the inputs.
Approximate Output Functions
We now seek to express the outputs in terms of inputs and slow states. Look-up tables are a common way to express functional relationships. A function is locally (interior to a grid box) approximated by a linear combination of its values at the corners of the box. For such a representation, the product space U×Xs must first be gridded. The system is driven to each of the gridpoints in the product space. Implicit in this step is settling of fast dynamics. Then, system outputs of interest are measured.
It must be noted that one can apply arbitrary settings (within permissible range) to actuators; however, no such direct control is necessarily possible over system states. One can manipulate a sequence of system inputs to drive the system to a certain state; such possibility depends on the controllability of the system.
Quasi-Steady State Measurements
The method proposed here employs a-priori gridding of the input space U alone, not U×Xs.
In Step 1, about each such gridpoint, the actuators are subjected to rate-limited simultaneous excursions up to the neighboring gridpoints.
In Step 2, the rates are chosen to be slower than the fast system states. Fast states are usually known to experienced test engineers and seldom need to be identified via the eigensystem analysis mentioned earlier. Note that the rates must as well be slower than the slowest sensor used.
In Step 3, a linear mathematical model of the output quantities is then derived. Note that this method does not guarantee complete coverage of the entire product space U×Xs, because we actively operate only in the product space of U. In such cases, we use the “nearest neighbor extrapolation” to fill the table.
Optimization
The optimization method proposed here is named “Sequential Linear (or Quadratic) Program with Sequential Relaxation,” and is described below:
In Step 1, for each of the hyperboxes defined by the gridding of U×Xs a linear programming (LP) problem is constructed and solved. The objective function, equality, and inequality constraints are all approximated by their respective linear representation constructed above. Note that the variables of the LP are subjected to the natural box constraint defined by the grid.
In Step 2, we begin with tightest inequality constraints (say, NOx emission less than some small number).
In Step 3, the solutions to this sequence of LP (one per hyperbox) are enumerated and the best solution is selected
In Step 4, if no solutions are found, the inequality constraints are sequentially relaxed (say, little more NOx emission allowed) and the process repeated
In Step 5, the resulting solution is then verified on the actual engine. Should this verification fail, the input space U is either refined (finer gridding) or enlarged (larger operating range during data collection) and the process repeated.
Optimization of Automated Engine Test Cell
Other hardware associated with system 20 are an analog to digital (A/D) converter 21b, an engine controller 21a, and external driver unit 22. The engine controller 21a is a card (or multiple cards, depending on the number of cylinders) is inserted into the computer 21. It contains an integrated circuitry that acts as an engine control chip and a specialized timer. The injector/spark/pulse width modulator (PWM) driver 22, controls various actuators of engine 21.
An engineer programs the computer 21 using a standard programming language, such as C, and the controller 21a transfers the instructions to the engine 10. Users of system 20 can change basic control configurations, such as number of cylinders, firing order, and other parameters, simply by entering numbers on an interface screen of computer 21.
The engine control capability of system 20 gives rise to the issue of what actuator settings are optimal for a desired operational feature of the engine. For purposes of this description, the desired operational feature is referred to as an “optimization variable”. Examples of optimization variables are the level of NOx emissions, the level of unburned hydrocarbon particles in the emissions, the brake specific fuel consumption (BSFC), the coefficient of variation for indicated mean effective pressure (COV-IMEP), and peak cylinder pressure. It should be understood that these are but a few of the many possible optimization variables that may be associated with an engine.
In accordance with the methods described above, the engine control algorithms compute optimum values for fueling, ignition, idle air, and exhaust gas recirculation from input received from the A/D converter 21b and engine controller 21a. The computed values are sent back to the engine controller 21a, which activates the driver 22, which controls various actuators of engine 10. In this manner, the above described method may be implemented with software as part of an engine control system. Once optimum values have been derived and tested, they may also be used as calibration values to calibrate the associated actuators of production engines.
Step 34 involves collection of response data. As described above, at each combination of settings, the engine is allowed to reach a “quasi-steady” state. Various combinations of settings are made and the engine is run to the quasi-steady state. In Step 35, from the response data, target values can be identified as well as the engine inputs associated with the optimal values. The target values may also be “optimal” in the sense that they are absolute optimal values for a particular response, or in the sense, that in view of all desired engine conditions, they are the “best” values.
It should be noted that the method can be used to identify inputs for steady state as well as transient engine output responses.
This application claims priority to U.S. provisional patent application No. 60/586,637, entitled “Use of Transient Data to Derive Steady State Calibrations for Dynamic Systems” filed on Jul. 9, 2004.
Number | Date | Country | |
---|---|---|---|
60586637 | Jul 2004 | US |