CN106908248A - The double regular empirical parameter automatic calibrating methods of weber burning of self-identifying list - Google Patents
The double regular empirical parameter automatic calibrating methods of weber burning of self-identifying list Download PDFInfo
- Publication number
- CN106908248A CN106908248A CN201710173226.5A CN201710173226A CN106908248A CN 106908248 A CN106908248 A CN 106908248A CN 201710173226 A CN201710173226 A CN 201710173226A CN 106908248 A CN106908248 A CN 106908248A
- Authority
- CN
- China
- Prior art keywords
- weber
- combustion
- burning
- fraction
- fired
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/04—Testing internal-combustion engines
- G01M15/042—Testing internal-combustion engines by monitoring a single specific parameter not covered by groups G01M15/06 - G01M15/12
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
It is an object of the invention to provide the double regular empirical parameter automatic calibrating methods of weber burning of self-identifying list.The present invention judges to need the weber equation number of selection according to fraction test data has been fired：For weber equation number be 1 when, the value according to a preliminary estimate of weber parameter is drawn using algebraic analysis, final estimate is then drawn using least-squares algorithm.For weber equation number be 2 when, determine that method draws combustion phase burble point using combustion phase burble point, fraction test data will have been fired to be divided into two parts by combustion phase burble point and processed accordingly, a weber parameter is worth according to a preliminary estimate to be shown using algebraic analysis to two parts data respectively, then final estimate is drawn using least-squares algorithm.The achievable self-identifying weber equation number of the present invention is simultaneously calibrated automatically to weber equation parameter such that it is able to is fast and accurately built the zero dimension combustion model based on weber burning rule, and be can guarantee that the stability and optimality of calibration result.
Description
Technical field
The present invention relates to a kind of combustion parameters calibration method.
Background technology
In order to solve increasingly serious problem of environmental pollution, International Emissions regulation is more and more harsher, limits having for engine
Evil emission, causes manufacturer, and the emission control to engine is seen particularly important.And in the emission performance of diesel engine and cylinder
Combustion process has a close contact, thus the emission control for realizing to the real-time control of combustion process to engine have it is important
Meaning.With the fast development of computer technology, computer simulation technique has flourishing vitality, by reality system
Abstract imitation, take out system model, people are simulated experimental study to such model on computers, have both dropped saving
Research and production costs, reduce risk, also improve scientific research efficiency.So reliability and accuracy of system model, directly
Connect the reliability and accuracy for determining simulation result.In field of internal combustion engine, the zero dimension burning based on weber (Wiebe) burning rule
Model form is simple, and modeling difficulty is small, while possessing certain simulation accuracy in certain condition range.With this weber
Wiebe) based on combustion model, successfully to direct-injection, non-direct-injection, two stroke diesel engine has carried out cylinder internal pressure to numerous researchers
The prediction of power and temperature.Weber Wiebe) empirical parameter of burning rule can directly affect the standard of weber (Wiebe) combustion model
True property, has document is studied how the empirical parameter of weber (Wiebe) combustion model is calibrated, such as algebraic analysis method
And least-squares algorithm, but algebraic analysis method and least-squares algorithm respectively have advantage and disadvantage.The former good stability, it is not necessary to give
Determine initial value, but it cannot be guaranteed that the optimality of calibration parameter；The latter can ensure the local optimality of calibration parameter, but convergence and
Calibration result depends on given initial value.It is therefore desirable to consider to be combined two methods, make both advantage and disadvantage complementary, most
The quickly and accurately regular empirical parameter of calibration weber burning is realized eventually.Research shows that single weber (Wiebe) burning rule is only suitable
For a kind of combustion phase or two kinds of combustion phases with slight mixing combustion process emulation, and for obvious two kinds combustions
Burning the combustion process of phase blending can not realize preferably emulation.Double weber (Wiebe) burning rules can be to two kinds of combustion phases
The combustion process of blending is preferably emulated, but calibration difficulty is higher, therefore the single pair of selection of weber (Wiebe) equation number
Need to weigh compromise between calibration difficulty and precision.For the combustion process for giving, how the Wei used needed for automatic identification
Primary (Wiebe) equation number and automatically calibration show that weber (Wiebe) equation empirical parameter is very crucial.
The content of the invention
It is an object of the invention to provide to have fired fraction test data to rely on rapidly and precisely automatic calibration parameter
The double regular empirical parameter automatic calibrating methods of webers burning of self-identifying list.
The object of the present invention is achieved like this：
The double regular empirical parameter automatic calibrating methods of weber burning of self-identifying list of the present invention, it is characterized in that：
(1) combustion test is carried out to diesel engine, test data sequence is collectedWhereinIt is crank angle, xbBe andIt is corresponding to have fired fraction, burning fitting initial angleIt is taken as the 1% corresponding crank angle of the fraction of combustion, burning fitting terminal
AngleIt is taken as the 99% corresponding crank angle of the fraction of combustion；
(2) by test data sequenceLinearisation：According to having fired fraction test data sequenceCalculate first
Go outWithValue according to a preliminary estimateWithWherein Turn to have fired corresponding bent axle when fraction is zero
Angle, if having fired fraction test data is consistently greater than zero, using data initial point correspondence crank angle asBy single weber of equation
Linearized, made Realize test data sequenceIt is linear
Change, the data sequence after linearisation is
(3) weber equation number is determined：It is preset asIt is rightData sequence carries out linear fit, draws fitting essence
Degree R2, the single double weber resolution E of setting, if R2>=E, weber equation number is identified as 1；If R2＜ E, the identification of weber equation number
It is 2；
(4) corresponding weber parameter automatic calibrating method is used for the weber equation number for determining, draws a weber equation
Parametric calibration result：
When weber equation number is identified as 1, according to fraction test value has been fired, draw firstWithValue according to a preliminary estimateWithWherein To have fired corresponding crank angle when fraction is zero, if having fired fraction experiment number
According to consistently greater than zero, using data initial point correspondence crank angle asThen to data sequenceLinear fit is carried out, is obtained
Go out fit slope A, by m0=A-1 calculates m0, wherein m0It is the initial value of fire burning index；Then byCalculate efficiency of combustion factor a；With
As fit equation is treated, with m0、WithRespectively as m,WithIterative initial value, using Nonlinear Least-Square Algorithm
Fitting draw m,WithFinal estimate；
When weber equation number is identified as 2, according to fraction test value has been fired, draw firstWithValue according to a preliminary estimateWithWherein To have fired corresponding crank angle when fraction is zero, if having fired fraction experiment number
According to consistently greater than zero, using data initial point correspondence crank angle asThen to linearisation afterData validation burning phase
Position burble point p, is to find a point p so that before and after this pointWithData carry out the synthesis of fitting a straight line respectively
R2Precision reaches maximum；Fraction test data sequence will have been fired according to combustion phase burble point pIt is divided into two parts, i.e.,WithWhereinx1b=[xb(1),xb(2),…,xb(p)],x2b=[xb(p+1),xb(p+2),…,xb(n)]；α0=xbP () is used as premixed combustion
Ratio initial value, and to x1bAnd x2bIt is normalized： Order
Realizing respectively willWithLinearisation, it is rightWithTwo parts data sequence enters line respectively
Property fitting, fit slope A is drawn respectively1And A2, by m10=A1-1、m20=A2- 1 draws m1 respectively0And m20, wherein m10And m20
Respectively premixed combustion fire burning index initial value, diffusion combustion fire burning index initial value；With α0、m10、m20、WithRespectively as α, m1,m2、WithIterative initial value, using non-linear least square calculate
Method fitting draw α, m1,m2、WithFinal estimate；
(5) output weber equation number and correspondence weber equation parameter collection, so as to complete the double weber burning rule of self-identifying list
Calibrate then and automatically the empirical parameter for drawing weber burning rule.
The present invention can also include：
1st, the determination method of combustion phase burble point p is as follows：
Assuming that data separating point i is by data sequenceIt is divided into two parts, the part before i-th data isPart after i-th data isIt is rightWithFitting a straight line, two are carried out respectively
The linear fit precision of partial data is respectively R2 1And R2 2, synthesis precision R2For：R2(i)=[R2 1×i+R2 2× (n-i)]/n,
Wherein n is total data amount check, and data separating point i can be made to be changed by 1~n, obtains synthesis precision R respectively successively2, then taking makes
Obtain synthesis precision R2Data separating point i when reaching maximum is used as combustion phase burble point p.
Advantage of the invention is that：The present invention according to weber (Wiebe) burn rule, with fired fraction test data be according to
Support, using double weber (Wiebe) algorithms of original self-identifying list, automatic identification weber (Wiebe) equation number, and according to identification
Weber (Wiebe) the equation number for arriving uses the corresponding automatic calibration algorithm of weber parameter, final to realize quickly and accurately certainly
The method that dynamic calibration draws the regular empirical parameter of weber burning.The double regular empirical parameters of weber (Wiebe) burning of self-identifying list are certainly
Dynamic calibration method realizes automatic identification weber (Wiebe) equation number using original single double weber (Wiebe) self-identifying algorithms,
Algebraic analysis method and least-squares algorithm are combined, makes both advantage and disadvantage complementary, realize the regular experience of weber (Wiebe) burning
The automatic calibration of parameter, the method automatic identification weber (Wiebe) equation number, during parametric calibration convergence and stability compared with
Good, accuracy is higher, is single double weber (Wiebe) the burning rules of researcher's selection in the industry and calibration weber (Wiebe) combustion
Burn regular empirical parameter and great convenience is provided.
Brief description of the drawings
Fig. 1 is flow chart of the invention.
Specific embodiment
Illustrate below in conjunction with the accompanying drawings and the present invention is described in more detail：
With reference to Fig. 1, it is usually used in single, double weber (Wiebe) the burning rule of internal combustion engine zero dimension burning modeling, respectively such as formula
(1) parametric equation and shown in formula (2).
In formula (1)：xbIt is fuel percentage；M is combustion quality index；A is the efficiency of combustion factor；- instantaneous bent
Shaft angle；- combustion continuation angle；- timing of combustion.Wherein m, a,WithBe empirical parameter to be calibrated.Formula (2)
In：xbIt is fuel percentage；m1And m2Respectively first weber (Wiebe) and second weber of burning product of (Wiebe) equation
Matter index；a1And a2Respectively first weber (Wiebe) and second weber of efficiency of combustion factor of (Wiebe) equation；- instantaneous
Crank angle；WithRespectively first weber (Wiebe) and second weber of combustion continuation angle of (Wiebe) equation；WithRespectively first weber (Wiebe) and second weber of timing of combustion of (Wiebe) equation.Wherein α, m1、m2、WithIt is empirical parameter to be calibrated.
First using beginning of ignition test value as the discreet value of timing of combustionWithAs combustion duration
Estimate linearization process is carried out according to weber (Wiebe) equation to having fired fraction test data；Then to the examination after treatment
Test data and use linear fit, with the R of linear fit2Precision uses the standard of weber (Wiebe) equation number as weighing, with
Resolutions of the value E as weber (Wiebe) equation number is manually set, works as R2>=E, weber (Wiebe) equation number elects 1 as,
Work as R2＜ E, weber (Wiebe) equation number elects 2 as；Corresponding weber (Wiebe) is used according to the weber equation number for determining
Parameter automatic calibrating method, finally draws weber (Wiebe) the equation empirical parameter calibrated.It is first when weber equation number is 1
First show that weber (Wiebe) equation parameter is worth according to a preliminary estimate using algebraic analysis method, then join in this, as weber (Wiebe)
Number iterative initial value, the final estimate of weber (Wiebe) parameter is drawn using Nonlinear Least-Square Algorithm calibration.Weber equation
When number is 2, determine that method draws combustion phase burble point using combustion phase burble point proposed by the present invention first, secondly root
Test data is divided into two parts according to combustion phase burble point, and respective handling is carried out respectively to this two parts data, it is then right
Two parts data after treatment are respectively adopted algebraic analysis method and show that weber (Wiebe) equation is worth according to a preliminary estimate, finally use
Nonlinear Least-Square Algorithm calibration draws double weber (Wiebe) final estimates of equation empirical parameter.
The double regular empirical parameter automatic calibrating methods of weber (Wiebe) burning of self-identifying list, calculation process is as follows：
Step one：Import the fraction test data sequence of combustion for determiningWhereinIt is crank angle, xbBe andIt is right
The combustion fraction answered, burning fitting initial angleThe corresponding crank angle of the fraction of combustion slightly larger than 0 (preferably 1%) is taken as,
Burning fitting terminal angleIt is taken as being slightly less than the corresponding crank angle of the fraction of combustion of 1 (preferably 99%), and extracts
Between test data.
Step 2：By test data sequenceLinearisation.According to having fired fraction test data sequenceCount first
DrawWithValue according to a preliminary estimateWithWherein To have fired corresponding bent axle when fraction is zero
Corner is (to improve the applicability and stability of the method, if having fired fraction test data is consistently greater than zero, with data initial point correspondence
Crank angle conduct).Single weber of equation is linearized, is made
Realize test data sequenceLinearisation, the data sequence after linearisation is
Step 3：It is determined that weber (Wiebe) equation number.It is preset asIt is rightData sequence carries out Linear Quasi
Close, draw fitting precision R2, (E is generally chosen for the number between 0.99~1 to single double weber (Wiebe) the resolution E of setting, preferably
For 0.995), if R2>=E, weber (Wiebe) equation number is identified as 1；If R2＜ E, weber (Wiebe) equation number is identified as
2。
Step 4：The weber equation number determined for step 3 uses corresponding weber parameter automatic calibrating method, obtains
Go out a weber equation parameter calibration result.When weber (Wiebe) equation number is identified as 1, according to fraction test value has been fired, first
Go outWithValue according to a preliminary estimateWithWherein Turn to have fired corresponding bent axle when fraction is zero
Angle is (bent with data initial point correspondence if having fired fraction test data is consistently greater than zero to improve the applicability and stability of the method
Shaft angle conduct)；Then to data sequenceLinear fit is carried out, fit slope A is drawn, by m0=A-1 is calculated
m0, wherein m0It is the initial value of fire burning index；Then by
Calculate efficiency of combustion factor a, preferably definite value 4.605, fraction combustion duration has been fired corresponding to 0%~99%；With right
It is required that the formula (1) in 1 is used as fit equation is treated, with m0、WithRespectively as m,WithIterative initial value, using non-
Linear least-squares algorithm fitting draw m,WithFinal estimate.When weber (Wiebe) equation number is identified as 2, root
According to fraction test value has been fired, draw firstWithValue according to a preliminary estimateWithWherein To have fired
Corresponding crank angle is (to improve the applicability and stability of the method, if fired the fraction test data beginning when fraction is zero
Eventually be more than zero, using data initial point correspondence crank angle as)；Then to linearisation afterData validation combustion phase
Burble point p, is to find a point p so that before and after this pointWithData carry out the comprehensive R of fitting a straight line respectively2
Precision reaches maximum；Fraction test data sequence will have been fired according to combustion phase burble point pIt is divided into two parts, i.e.,WithWhereinx1b=[xb(1),xb(2),…,xb(p)],x2b=[xb(p+1),xb(p+2),…,xb(n)]。α0=xbP () is used as premixed combustion
Ratio initial value, and to x1bAnd x2bIt is normalized： Order
Realizing respectively willWithLinearisation, it is rightWithTwo parts data sequence enters line respectively
Property fitting, fit slope A is drawn respectively1And A2, by m10=A1-1、m20=A2- 1 draws m1 respectively0And m20, wherein m10And m20
Respectively premixed combustion fire burning index initial value, diffusion combustion fire burning index initial value；A1 and a2 are both preferably definite value 4.605
(being not limited to 4.605), fraction combustion duration has been fired corresponding to 0%~99%, used as premixed combustion and diffusion combustion efficiency
Factor.With α0、m10、m20、WithRespectively as α, m1,m2、WithRepeatedly
For initial value, using Nonlinear Least-Square Algorithm fitting draw α, m1,m2、WithFinal estimate.
Step 5：Output weber equation number and correspondence weber equation parameter collection.
So far according to the test data for determining, it is possible to which double weber (Wiebe) the burning rules of self-identifying list are simultaneously calibrated automatically
Draw the empirical parameter of weber burning rule.
The determination Method And Principle of the combustion phase burble point p described in step 4 is as follows：
Combustion phase burble point p.Assuming that data separating point i is by data sequenceBe divided into two parts, i-th data it
Preceding part isPart after i-th data isIt is rightWithCarry out respectively straight
Line is fitted, and the linear fit precision of two parts data is respectively R2 1And R2 2, synthesis precision R2It is defined as follows shown in formula, wherein n is
Total data amount check, can be such that data separating point i is changed by 1~n, obtain synthesis precision R respectively successively2, then take so that comprehensive
Precision R2Data separating point i when reaching maximum is used as combustion phase burble point p.
The double regular empirical parameter automatic calibrating method concrete principles of weber (Wiebe) burning of self-identifying list are as follows：
Shown in single weber of equation such as formula (1).
ForBy xbFor 0 when corresponding crank angle determine；ForBy xbIt is 0~0.99 corresponding crank angle
Phase determines.
Generally by the corresponding crank angle of 50% thermal dischargeAs combustion centre, corresponding to formula (1), it is：
Arrangement can be obtained：
Formula (5) is converted into：
From formula (6), by calculating m by the slope for asking G and H.
Assuming that burning fitting starting point is xbs, corresponding crank angle isBurning fitting terminal is xbc, corresponding song
Shaft angle isCan be drawn by formula (1)：
For combustion duration elect as 0~99% fired the corresponding crank angle of fraction during when,
Can be seen by formula (7)
Go out, also can be different for different combustion duration and m, a.A of the present invention is preferably 4.605.
Single weber of burning rule of formula (1) is converted into：
From formula (8), the big I of m represents the speed of burning velocity.
For single combustion phase or the combustion process of slight double combustion phase blending, due to being fired in whole combustion process
Speed is burnt to be more or less the same, therefore the linear relationship of K and G is preferable, and the R of linear fit is used to K and G2Precision is higher；For obvious
For the combustion process of double combustion phases blending, because burning velocity difference is larger in whole combustion process, therefore K and G line
Sexual intercourse is poor, and the R of linear fit is used to K and G2Precision is relatively low.Analyzed more than, the R of K and G linear fits2Precision
Level can characterize the order of severity of double combustion phase blending, and E is weber equation number resolution, chooses R2>=E is used as single Wei
Primary equation mark, chooses R2＜ E are identified as double webers of equations, to realize single double weber self-recognition functions.E is preferred in the present invention
It is 0.995.The mathematic(al) representation of single double weber equation self-identifying algorithms is as follows：
In formula, Num is weber equation number.
For DID engine, always more or less there is premixed combustion in combustion process in addition to diffusion combustion mode
Pattern, double webers of equation parameter calibrations, find two kinds of burble points of combustion mode ten of premixed combustion and diffusion combustion for convenience
Divide key.This software takes a kind of original method to determine combustion mode burble point, confirms that arthmetic statement is as follows.
By formula (8) it can be seen that the size of m can reflect the speed of combustion process.For premixed combustion mode, burning speed
Degree is very fast, therefore m is smaller；For diffusion combustion mode, burning velocity is slower, therefore m is larger.Combustion mode burble point is to look for
Cause that this point G before and after and K data carry out the comprehensive R of fitting a straight line respectively to a point2Precision reaches maximum.Assuming that
The R of linear fit is used to the G and K before this point (p-th data point)2Precision is R2 1, the G and K after this point are using linear
The R of fitting2Precision is R2 2, define synthesis R2P () precision is：
Then, R is taken2Corresponding p is the combustion mode burble point to be looked for when () is for maximum p.
Fraction test data will be fired according to data separating point p and has been divided into two parts,x1b=xb(1:P),x2b=xb(p+1:End),
α0=xb(p), and to x1bAnd x2bProcessed：To treatment
Two groups of data afterwards carry out algebraic analysis respectively, draw m,WithDiscreet value, the iteration as least-squares algorithm is initial
Value.
Least-squares algorithm is usually used in nonlinear equation fitting problems, and its theory is as described below.
Given n is to independent variable and the test data (x of dependent variablei,yi), parameter set to be determined is β, the fitting side of selection
Journey is p (x, β), therefore, the quadratic sum of error is：
Least-squares algorithm is to obtain one group of β so that S (β) is minimum.
Currently preferred least-squares algorithm is Levenberg-Marquardt algorithms.Algorithm calculate start needs to
The initial value of fixed parameter set β to be calibrated.Afterwards, replaced using the value β+δ of new estimation in each iteration step β.In order to determine δ,
To equation p (xi, β+δ) and carry out Linear Estimation：
p(xi, β+δ) and=p (xi,β)+Jiδ (12)
In formulaIt is the gradient relative to β.When S (β) reaches minimum, S (β) is changed into the gradient of δ
0.According to formula (12), the single order of formula (11) is estimated as follows formula：
It is expressed as in vector form：
S(β+δ)≈||y-f(β)-Jδ||2 (14)
Formula (14) is to J derivations, and it is zero to make derived function, can be obtained：
(JTJ) δ=JT[y-f(β)] (15)
Levenberg-Marquardt algorithms are improved formula (15), are changed into following formula：
(JTJ+ λ I) δ=JT[y-f(β)] (16)
In formula, I is unit matrix, and λ is damped coefficient, the step-length for adjusting each iteration.When λ is zero, formula (16) is moved back
Formula (15) is turned to, is Gauss-Newton's algorithm；When λ is higher value, formula (16) deteriorates to gradient descent algorithm.
In order to improve the convergence rate of formula (16), J is used to formula (16)TJ replaces I, final Levenberg-Marquardt
Algorithm is shown below：
(JTJ+λdiag(JTJ)) δ=JT[y-f(β)] (17)
Because least-squares algorithm needs to give the initial value of parameter to be calibrated, and least-squares algorithm when calculating and starting
Convergence and iterative calculation time stronger to initial value dependence and final result it cannot be guaranteed that Global Optimality, can only
Ensure local optimality, therefore it is very crucial to give more rational initial value.The present invention is using algebraic analysis method according to experiment number
According to the weber empirical parameter for calculating as the initial value of least-squares algorithm, then further iterated to calculate, such as
This can ensure the optimality of calibration result, the convergence of least square method can largely be improved again, and can reduce iteration meter
Evaluation time.
Claims (2)
1. double regular empirical parameter automatic calibrating methods of webers burning of self-identifying list, it is characterized in that：
(1) combustion test is carried out to diesel engine, test data sequence is collectedWhereinIt is crank angle, xbBe andCorrespondence
Combustion fraction, burning fitting initial angleIt is taken as the 1% corresponding crank angle of the fraction of combustion, burning fitting terminal angleTake
It is the 99% corresponding crank angle of the fraction of combustion；
(2) by test data sequenceLinearisation：According to having fired fraction test data sequenceCalculate firstWithValue according to a preliminary estimateWithWherein To have fired corresponding crank angle when fraction is zero,
If having fired fraction test data is consistently greater than zero, using data initial point correspondence crank angle asSingle weber of equation is carried out
Linearisation, order Realize test data sequenceLinearisation,
Data sequence after linearisation is
(3) weber equation number is determined：It is preset asIt is rightData sequence carries out linear fit, draws fitting precision
R2, the single double weber resolution E of setting, if R2>=E, weber equation number is identified as 1；If R2＜ E, weber equation number is identified as
2；
(4) corresponding weber parameter automatic calibrating method is used for the weber equation number for determining, draws a weber equation parameter
Calibration result：
When weber equation number is identified as 1, according to fraction test value has been fired, draw firstWithValue according to a preliminary estimateWith
Wherein To have fired corresponding crank angle when fraction is zero, if having fired fraction test data is consistently greater than zero, with
The correspondence crank angle conduct of data initial pointThen to data sequenceLinear fit is carried out, fit slope A is drawn, by m0=A-1
Calculate m0, wherein m0It is the initial value of fire burning index；Then by
Calculate efficiency of combustion factor a；WithAs fit equation is treated, with m0、WithRespectively
As m,WithIterative initial value, using Nonlinear Least-Square Algorithm fitting draw m,WithFinal estimate；
When weber equation number is identified as 2, according to fraction test value has been fired, draw firstWithValue according to a preliminary estimateWith
Wherein To have fired corresponding crank angle when fraction is zero, if fired fraction test data be consistently greater than
Zero, using data initial point correspondence crank angle asThen to linearisation afterData validation combustion phase burble point p,
It is to find a point p so that before and after this pointWithData carry out the comprehensive R of fitting a straight line respectively2Precision reaches most
Greatly；Fraction test data sequence will have been fired according to combustion phase burble point pIt is divided into two parts, i.e.,With
Whereinx1b=[xb(1),xb(2),…,xb(p)],x2b=
[xb(p+1),xb(p+2),…,xb(n)]；
α0=xb(p) as premixed combustion ratio initial value, and to x1bAnd x2bIt is normalized： Order Realizing respectively willWithLinearisation, it is rightWithTwo parts data sequence carries out linear fit respectively, and fit slope A is drawn respectively1And A2, by m10=
A1-1、m20=A2- 1 draws m1 respectively0And m20, wherein m10And m20Respectively premixed combustion fire burning index initial value, spreads combustion
Burn fire burning index initial value；With α0、m10、m20、WithRespectively as α, m1,m2、WithIterative initial value, using Nonlinear Least-Square Algorithm fitting draw α, m1,m2、With
Final estimate；
(5) output weber equation number and correspondence weber equation parameter collection, so as to complete the double weber burning rules of self-identifying list simultaneously
Automatic calibration draws the empirical parameter of weber burning rule.
2. double regular empirical parameter automatic calibrating methods of webers burning of self-identifying list according to claim 1, it is characterized in that：
The determination method of combustion phase burble point p is as follows：
Assuming that data separating point i is by data sequenceIt is divided into two parts, the part before i-th data isI-th
Part after individual data isIt is rightWithFitting a straight line, the line of two parts data are carried out respectively
Property fitting precision is respectively R2 1And R2 2, synthesis precision R2For：R2(i)=[R2 1×i+R2 2× (n-i)]/n, wherein n is total number
According to number, data separating point i can be made to be changed by 1~n, obtain synthesis precision R respectively successively2, then take so that synthesis precision R2Reach
To data separating point i during maximum as combustion phase burble point p.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710173226.5A CN106908248A (en) | 2017-03-22 | 2017-03-22 | The double regular empirical parameter automatic calibrating methods of weber burning of self-identifying list |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710173226.5A CN106908248A (en) | 2017-03-22 | 2017-03-22 | The double regular empirical parameter automatic calibrating methods of weber burning of self-identifying list |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106908248A true CN106908248A (en) | 2017-06-30 |
Family
ID=59195161
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710173226.5A Pending CN106908248A (en) | 2017-03-22 | 2017-03-22 | The double regular empirical parameter automatic calibrating methods of weber burning of self-identifying list |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106908248A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107844625A (en) * | 2017-09-13 | 2018-03-27 | 上海机电工程研究所 | Spary coating type ablation jet vane, which rises, burns moment point estimation approach |
| CN108918588A (en) * | 2018-04-02 | 2018-11-30 | 广西玉柴机器股份有限公司 | The method of simulated high-pressure common rail diesel engine combustion state |
| WO2019047863A1 (en) * | 2017-09-06 | 2019-03-14 | Xuesong Li | Apparatuses and methods for optical calibration |
| CN112065584A (en) * | 2020-08-06 | 2020-12-11 | 南京瑞华动力科技有限公司 | System and method for controlling Weber correction index of gas turbine fuel |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007248119A (en) * | 2006-03-14 | 2007-09-27 | Toyota Motor Corp | Method for determining wiebe function parameter and device for presuming heat release rate of internal combustion engine |
| CN102722699A (en) * | 2012-05-22 | 2012-10-10 | 湖南大学 | Face identification method based on multiscale weber local descriptor and kernel group sparse representation |
| US20140139557A1 (en) * | 2012-11-19 | 2014-05-22 | Samsung Display Co., Ltd., | Display apparatus and control method for saving power thereof |
| CN105157055A (en) * | 2015-06-24 | 2015-12-16 | 黄红林 | Method for identifying combustion model of circulating fluidized bed on basis of least squares |
| CN105868544A (en) * | 2016-03-25 | 2016-08-17 | 华北电力大学（保定） | A method of analyzing insulating material performance by using three-parameter Weibull distribution to process flashover voltages |
-
2017
- 2017-03-22 CN CN201710173226.5A patent/CN106908248A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007248119A (en) * | 2006-03-14 | 2007-09-27 | Toyota Motor Corp | Method for determining wiebe function parameter and device for presuming heat release rate of internal combustion engine |
| CN102722699A (en) * | 2012-05-22 | 2012-10-10 | 湖南大学 | Face identification method based on multiscale weber local descriptor and kernel group sparse representation |
| US20140139557A1 (en) * | 2012-11-19 | 2014-05-22 | Samsung Display Co., Ltd., | Display apparatus and control method for saving power thereof |
| CN105157055A (en) * | 2015-06-24 | 2015-12-16 | 黄红林 | Method for identifying combustion model of circulating fluidized bed on basis of least squares |
| CN105868544A (en) * | 2016-03-25 | 2016-08-17 | 华北电力大学（保定） | A method of analyzing insulating material performance by using three-parameter Weibull distribution to process flashover voltages |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019047863A1 (en) * | 2017-09-06 | 2019-03-14 | Xuesong Li | Apparatuses and methods for optical calibration |
| CN107844625A (en) * | 2017-09-13 | 2018-03-27 | 上海机电工程研究所 | Spary coating type ablation jet vane, which rises, burns moment point estimation approach |
| CN108918588A (en) * | 2018-04-02 | 2018-11-30 | 广西玉柴机器股份有限公司 | The method of simulated high-pressure common rail diesel engine combustion state |
| CN112065584A (en) * | 2020-08-06 | 2020-12-11 | 南京瑞华动力科技有限公司 | System and method for controlling Weber correction index of gas turbine fuel |
| CN112065584B (en) * | 2020-08-06 | 2021-12-17 | 南京瑞华动力科技有限公司 | System and method for controlling Weber correction index of gas turbine fuel |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106908248A (en) | The double regular empirical parameter automatic calibrating methods of weber burning of self-identifying list | |
| Morgan et al. | Mapping surrogate gasoline compositions into RON/MON space | |
| Powell | A review of IC engine models for control system design | |
| Yeliana et al. | Estimation of double-Wiebe function parameters using least square method for burn durations of ethanol-gasoline blends in spark ignition engine over variable compression ratios and EGR levels | |
| US7181332B1 (en) | Method for controlling an operating condition of a vehicle engine | |
| Winterbone et al. | A wholly dynamic model of a turbocharged diesel engine for transfer function evaluation | |
| WO2020034630A1 (en) | Neural network-based cell delay prediction method and cell delay sensitivity calculation method | |
| JP2011106334A (en) | Method of estimating heat release rate of engine using wiebe function model | |
| CN110118130B (en) | Control device for internal combustion engine | |
| CN104615840B (en) | The modification method and system of a kind of digital simulation model | |
| JP2013100830A (en) | Combustion control in internal combustion engine | |
| Grill et al. | Quasi-dimensional and empirical modeling of compression-ignition engine combustion and emissions | |
| CN109282916A (en) | The method and apparatus for predicting dual fuel engine piston face maximum temperature | |
| WO2022237074A1 (en) | Neural network model-based intake air amount estimation method for secondary inflation model of gasoline engine | |
| CN107092723A (en) | The regular empirical parameter automatic calibrating method of double webers of burnings | |
| Sjerić | New physically based sub-models for the cycle-simulation of spark-ignition engine | |
| JP2004068729A (en) | Adapting method for engine control parameter and its system | |
| CN106548029A (en) | The determination method of the regular empirical coefficient of weber burning | |
| Linse et al. | An experimental and numerical investigation of turbulent flame propagation and flame structure in a turbo-charged direct injection gasoline engine | |
| CN105631093B (en) | A kind of Design of Mechanical Structure method based on M-BSWA multiple-objection optimizations | |
| CN106960092A (en) | The regular empirical parameter automatic calibrating method of weber burning | |
| CN116591846A (en) | Fuel injection control method and device, electronic equipment and storage medium | |
| CN112883653B (en) | Artificial intelligence-based modeling method for real-time engine model | |
| Liu et al. | A model for on-line monitoring of in-cylinder residual gas fraction (RGF) and mass flowrate in gasoline engines | |
| Sjerić et al. | Cycle-simulation turbulence modelling of IC engines |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170630 |