CN108197392A - The assist characteristic curve design method of automobile electric booster steering system based on SOC - Google Patents

The assist characteristic curve design method of automobile electric booster steering system based on SOC Download PDF

Info

Publication number
CN108197392A
CN108197392A CN201810010887.0A CN201810010887A CN108197392A CN 108197392 A CN108197392 A CN 108197392A CN 201810010887 A CN201810010887 A CN 201810010887A CN 108197392 A CN108197392 A CN 108197392A
Authority
CN
China
Prior art keywords
soc
power
smoothing
assisted
value
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.)
Granted
Application number
CN201810010887.0A
Other languages
Chinese (zh)
Other versions
CN108197392B (en
Inventor
武健
徐广飞
吴海荣
王锋波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liaocheng University
Original Assignee
Liaocheng University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Liaocheng University filed Critical Liaocheng University
Priority to CN201810010887.0A priority Critical patent/CN108197392B/en
Publication of CN108197392A publication Critical patent/CN108197392A/en
Application granted granted Critical
Publication of CN108197392B publication Critical patent/CN108197392B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Power Steering Mechanism (AREA)
  • Secondary Cells (AREA)

Abstract

The invention discloses a kind of assist characteristic curve design methods of the automobile electric booster steering system based on SOC, consider the situation of battery SOC, SOC is divided into three classes, respectively SOC<40%th, 40%<SOC<70% and SOC>70%;To SOC<40% situation carries out segment processing, to 40% corresponding to shaped form assist characteristic curve<SOC<70% situation carries out double smoothing processing;To SOC>70% situation carries out single exponential smoothing processing.The present invention rationally designs assist characteristic curve by considering on the basis of electric boosted automobile batteries energy management, power-assisted steering performance can be obviously improved, it avoids since assist motor power supply is discontinuous caused by the fast or big hand torque operation of the variation of hand-power square, and then cause and turn to difficulty etc., ensure that battery can generate stable electric energy output in the various electricity stages and make to stablize by the power steering that assist characteristic curve calculates, improve power-assisted ability to the greatest extent while ensureing that power-assisted stablizes safe, it is energy saving.

Description

The assist characteristic curve design method of automobile electric booster steering system based on SOC
Technical field
The invention belongs to automobile technical fields, are related to electric boosting steering system and Energy Management System, especially one Kind considers the assist characteristic Curve Design in the case of battery SOC.
Background technology
The assist characteristic curve of EPS refers to the power-assisted electric current of motor and steering wheel input torque, Vehicle Speed, preceding Take turns the relationship of air pressure, front shaft again etc. between parameters, consider under normal circumstances steering wheel input torque, Vehicle Speed this two A major parameter can meet control requirement.The assist characteristic curve of EPS is the key technology of EPS software systems.Assist characteristic Curve determines that controller (ECU) goes the size of control power-assisted electric current according to which type of target, goes to meet different driving cycles Under requirement to booster torquemoment.Basic assist characteristic curve is roughly divided into:Linear type, broken line type, shaped form.It can expire substantially Sufficient design requirement, but the actual effect of power-assisted is poor, then, the Curve guide impeller to assist characteristic curve occurs, wherein, there is ginseng The adjustable Curve-type Assistance Characteristics Curve Design of number carries out shape by the parameter of setting to Curve-type Assistance Characteristics curve Variation, can make curve variable at any time.It, can be according to small power-assisted area and big power-assisted area not with double-gradient smooth assist curves Independent calibration is carried out to parameter of curve with feel requirement, has the advantages that power-assisted smooth-going, parameter function is independent, is easy to debugging, energy Preferably meet expectation of the driver in terms of feel is turned to.Above design does not account for battery only from curve The influence that electricity generates vehicle booster, when hand-power square changes greatly and hand-power square is larger, battery probably can not The EPS electric currents met the requirements are provided, EPS system can power off suddenly, and steering boost system can not work moment, seriously affects and drives Sail safety.
Application No. is 200810202069.7 Chinese patents, and disclosing " has the individually adjusted shaped form power-assisted of parameter The electric boosting steering system of characteristic ";There are parameter adjustments to make the Curve Design period long for it, meanwhile, the reality of vehicle can not be contacted Border situation, so as to carry out curve adjustment according to the actual conditions of vehicle.
Application No. is 201010600546.2 Chinese patents, and disclosing " has the electronic of double-gradient smooth assist curves Servo steering system, " it is segmented the reality that can not ensure power-assisted curve during carrying out diclinic counting smooth processing there are broken line Whether application effect can guarantee sustainable supply of battery capacity etc..
Invention content
In order to overcome the above-mentioned deficiencies of the prior art, the present invention provides a kind of automobile electric power-assisted steerings based on SOC The assist characteristic curve design method of system.
A kind of assist characteristic curve design method of the automobile electric booster steering system based on SOC, the assist characteristic Curve design consider battery SOC situation, during Curve Design on the basis of virgin curve type assist characteristic curve into Row smoothing processing, while constrained to electric current needed for EPS system to ensure battery even if can in the case of low battery There is steady and sustained electric current to export, ensure the realization of the entire power-assisted processes of EPS, Curve Design includes the following steps:
Step 1) estimates the SOC of battery with kalman filtering theory;SOC is divided into three classes, respectively SOC< 40%th, 40%<SOC<70% and SOC>70%;
Step 2) is to SOC<40% situation carries out segment processing corresponding to shaped form assist characteristic curve, is divided into 5 sections 5 areas, respectively:Power-assisted area 1 is power-assisted dead zone, and power-assisted area 2 is small power-assisted area, and power-assisted area 3 is middle power-assisted area, and power-assisted area 4 is Big power-assisted area, power-assisted area 5 is saturation region;Curve smoothing processing is carried out respectively to 5 sections of assist characteristic curve, and power-assisted area 2 is hand Torque T1<Td<T2, handled using single exponential smoothing;Middle Resistance Level 3 is T2<Td<T3, handled using double smoothing;It helps greatly Power area 4 is T3<Td<T4, handled using Three-exponential Smoothing, improve the change rate of shaped form power-assisted curve, so as to stablize battery Electric current exports;
Step 3) is to 40%<SOC<70% situation carries out double smoothing processing to Curve-type Assistance Characteristics curve;
Step 4) is to SOC>70% situation carries out single exponential smoothing processing to Curve-type Assistance Characteristics curve;
During step 5) is to utilizing hand Calculating Torque during Rotary target current, i.e. Ia=ζ Td, proportionate relationship ζ therein is changed Processing,
Wherein, λ is transformation coefficient, determines the corresponding power-assisted curve of speed;Wherein, IaFor target current, TdFor hand-power square, ζ For hand-power square and the conversion relation of target current;
Step 6) is under different speeds, by obtaining the complete assist characteristic under different speeds about the continuous item of speed Curve, i.e.,
I=e-av·ζTd (2)
Wherein I be each vehicle speed condition under target current, e-avFor the compensation coefficient comprising speed, a is constant.
This patent considers the remaining capacity of battery in Curve Design, meanwhile, the design of curve is helped in virgin curve type On the basis of force characteristic curve, in the case where meeting practical application, the calculating of target current can be preferably realized, ensure real Border follow current stablizes output;During power-assisted curve is designed on the one hand be on the basis of shaped form power-assisted curve into Capable smoothing processing, effect are more good;On the other hand, this patent makes the design carried out on the basis of battery capacity is considered, It can ensure the sustainable supply of battery current.
The assist characteristic curve design method of automobile electric booster steering system based on SOC uses double card in step 1) Kalman Filtering theory makes the state of dynamical system the optimal estimation in minimum variance meaning, estimates applied to battery SOC, electricity Pond is seen as dynamical system, and SOC is an internal state of system;The general mathematical form of battery model is:
State equation:xk+1=Akxk+Bkuk+wk=f (xk,uk)+wk (3)
Observational equation:yk=ckxk+vk=g (xk,uk)+vk (4)
Wherein, the input vector u of systemk, to include battery current, temperature, residual capacity and internal resistance variable;
The output y of systemk, the operating voltage for battery;
Battery SOC is included in the quantity of state x of systemkIn, f (xk,uk) and g (xk,uk) be all determined by battery model it is non- Linear equation is linearized in calculating process;SOC algorithms are estimated, including SOC estimation and reaction evaluated error, association side The recursion equation of poor matrix, covariance matrix are used for providing evaluated error range, this equation is in battery model state equation In, SOC is described as to the foundation of state vector:
Specifically double card Kalman Filtering algorithm steps are:
Step (1) gives initial value:
To enable iteration Fast Convergent, by the relatively actual value of the initial value setting of SOC, R0Initial value then lead to It crosses current SOC and tables look-up and provide, other two state can be set as 0, thus just obtain X (0) and R (0);
Step (2) SOC estimates:
The system mode of kth step is estimated using the systematic parameter of -1 step of kth, the system mode then walked again with kth is estimated Count the systematic parameter of kth step:
Step (2.1) is estimated using the system mode that current integration walks kth:
X (kk-1)=As(k)X(k-1)+Bs(k)I(k)+ωs (6)
Wherein,
Wherein, ωsFor the process noise of system, substantially determined by the noise of electric current;
The optimal estimation of step (2.2) system mode:
X (k)=X (k | k-1)+Ks(k)(V(k)-V(k)') (9)
Wherein V (k) is the battery both end voltage measured, and V (k) ' is the terminal voltage estimated using battery model;
V (k) '=Fv(S(k))-R0(k)I(k)-URC1(k)-URC2(k)+v (10)
Step (2.3) seeks Ks (k) (kalman gain), calculates variance matrix:
Wherein QsFor the covariance matrix of systematic procedure noise, rsThe variance of measurement noise for voltage;Fv(S) it is open circuit Function of the voltage about SOC;
Step (2.4) updates variance matrix P, terminates the Kalman filtering of the SOC of this step;
Ps(k)=[I-Ks(k)Cs(k)]Ps(k|k-1) (14)
It is special to choose shaped form power-assisted for the assist characteristic curve design method of automobile electric booster steering system based on SOC Linearity curve function is:
Wherein,F (v) is the corresponding power-assisted curvilinear system of friction speed;F (v) is vehicle The function of fast v, v is faster, and F (v) is smaller, takes the exponential decrease function of v, G (Td) for different directions disk torque input when power-assisted Curve is increasing function, i.e. the input torque of steering wheel is bigger, light in order to turn to, and booster torquemoment also must be bigger.
The assist characteristic curve design method of automobile electric booster steering system based on SOC, used smoothing processing Or exponential smoothing, it is a kind of Time Series Analysis Forecasting method to grow up on the basis of the method for moving average, it is to pass through meter Exponential smoothing value is calculated, cooperation regular hour sequential forecasting models predict the future of phenomenon;Its principle is any phase Exponential smoothing value is all the weighted average of current period actual observation value and previous phase index smooth value, the fundamental formular of exponential smoothing It is:
St=α yt+(1-α)St-1 (16)
In formula, St-- the smooth value of time t;yt-- the actual value of time t;St-1-- the smooth value of time t-1;α -- it is smooth Constant, value range are [0,1];
According to smooth number difference, exponential smoothing is divided into:Single Exponential Smoothing, Secondary Exponential Smoothing Method and index three times Exponential smoothing;
For small power-assisted area 2, T1<Td<T2, handled using single exponential smoothing, predictor formula is:
yt+1'=ayt+(1-a)y't (17)
Wherein, yt+1' -- the smooth value S of the predicted value of t+1 phases, i.e. current period (t phases)t (1);yt-- the actual value of t phases;yt'-- The predicted value of t phases, i.e. last smooth value St-1 (1);The formula can be write again:yt+1'=y't+a(yt-y't), it is seen then that the next period Predicted value be current period predicted value again with using a as the sum of current period actual value of discount and predicted value error;
For middle power-assisted area 3, T2<Td<T3, handled using double smoothing, double smoothing is that an index is put down Sliding is smooth again, it is suitable for the time series of tool linear trend, and predictor formula is:
yt+m=(2+am/(1-a))y't-(1+am/(1-a))yt=(2y't-yt)+m(y't-yt)a/(1-a) (18)
In formula, yt=ay't-1+(1-a)yt-1
Obviously, double smoothing is a linear equation, and intercept is:(2yt'-yt), slope is:(yt'-yt)a/(1- A), independent variable is predicts number of days, double smoothing fundamental formular:
St (2)=aSt (1)+(1-a)St-1 (2) (19)
Yt+ Y=at+btT (20)
at=2St (1)-St (2) (21)
bt=(a/1-a) (St (1)-St (2)) (22)
In formula, St (1)-- the single exponential smoothing value S of t phasest (2)-- the double smoothing value α of t phases -- smoothing factor Yt+T-- t+T phase predicted value T-- elapse issue backward by the t phases;
For big power-assisted area 4, T3<Td<T4, handled using Three-exponential Smoothing, Three-exponential Smoothing prediction is secondary smooth On the basis of it is smooth again, predictor formula is:
yt+m=(3y't-3yt)+[(6-5a)y't-(10-8a)yt+(4-3a)yt]*am/2(1-a)2+(y't-2yt+y't)* a2m2/2(1-a)2 (23)
In formula, yt=ayt-1+(1-a)yt-1
Predicted value is the weighted sum of former observation, and different data are given with different power, and new data is to larger Power, legacy data give smaller power, by trend adjustment, add trend correction value, can improve exponential smoothing to a certain extent Prediction result, the formula of the exponential smoothing after adjustment are:
Include trend prediction (YtTt)=new prediction (Yt)+trend corrects (Tt)
There are three steps for the Smoothing Prediction of progress trend adjustment:
(1) the simple index number smoothing prediction (Y of t phases is calculated using previously described methodt);
(2) trend is calculated, formula is:
Tt=(1-b) Tt-1+b(Yt-Yt-1) (24)
Wherein, Tt=the t phases smoothed trend;Tt-1=the t phases last smoothed trend;B=selections become Gesture smoothing factor;Yt=to t phase simple index number smoothing predictions;Yt-1=to t phase last issue simple index number smoothing predictions,
(3) the Smoothing Prediction value (Y after trend adjustment is calculatedtTt), calculation formula is:
YtTt=Yt+Tt。 (25)
The technical solution adopted by the present invention is when being designed to assist characteristic curve, it is contemplated that the situation of SOC, especially It is to be changed greatly in hand-power square and greatly in the case of hand-power square, change the change rate of power-assisted curve, control the electric current of EPS in electricity Within the supply district of pond.
SOC numerical value is observed in real time by Kalman Filter observer, according to SOC numerical values reciteds to virgin curve type assist characteristic Curve carries out a point situation and handles, to SOC<40 situation is segmented virgin curve type assist characteristic curve smoothly to be located respectively Reason, to SOC>40% situation with regard to carrying out level-one respectively, two level exponential smoothing handles, while the demand current of EPS is carried out about Beam ensures lasting supply of the battery to EPS electric currents, so as to be that servo steering system persistently plays a role.
Compared with prior art, the beneficial effects of the invention are as follows battery electricity is considered when assist characteristic curve is designed The situation of SOC is measured, steering boost system can be made to be run in the case of any battery capacity SOC, ensure the steering of driver Safety.
Description of the drawings
Fig. 1 is the assist characteristic curve of the prior art.
Fig. 2 is SOC of the present invention<Assist characteristic curve in the case of 40%.
Fig. 3 is the present invention 40%<SOC<Assist characteristic curve in the case of 70%.
Fig. 4 is SOC of the present invention>Assist characteristic curve in the case of 70%.
Fig. 5 is the assist characteristic Curve Design flow chart the present invention is based on SOC.
Specific embodiment
Below in conjunction with attached drawing to further instruction of the present invention.
Fig. 1 is engaged it is found that the design of the prior art does not account for battery capacity for vehicle only from curve The influence that power-assisted generates, when hand-power square changes greatly and hand-power square is larger, battery can not probably be provided and be met the requirements EPS electric currents, EPS system can power off suddenly, and steering boost system can not work moment, seriously affect driving safety.
Fig. 2-Fig. 4 is engaged to understand;A kind of assist characteristic of automobile electric booster steering system based on SOC of the present invention is bent Line design method, the design of the assist characteristic curve considers the situation of battery SOC, in virgin curve type during Curve Design It is smoothed on the basis of assist characteristic curve, even if while being constrained to ensure battery to electric current needed for EPS system Can also there be steady and sustained electric current to export in the case of low battery, ensure the realization of the entire power-assisted processes of EPS, Curve Design Include the following steps:
Step 1) estimates the SOC of battery with kalman filtering theory;SOC is divided into three classes, respectively SOC< 40%th, 40%<SOC<70% and SOC>70%;
Step 2) is to SOC<40% situation carries out segment processing corresponding to shaped form assist characteristic curve, is divided into 5 sections 5 areas, respectively:Power-assisted area 1 is power-assisted dead zone, and power-assisted area 2 is small power-assisted area, and power-assisted area 3 is middle power-assisted area, and power-assisted area 4 is Big power-assisted area, power-assisted area 5 is saturation region;Curve smoothing processing is carried out respectively to 5 sections of assist characteristic curve, and power-assisted area 2 is hand Torque T1<Td<T2, handled using single exponential smoothing;Middle Resistance Level 3 is T2<Td<T3, handled using double smoothing;It helps greatly Power area 4 is T3<Td<T4, handled using Three-exponential Smoothing, improve the change rate of shaped form power-assisted curve, so as to stablize battery Electric current exports;For small power-assisted area 2, T1<Td<T2, handled using single exponential smoothing.In small 2 sections of power-assisted area, n point is chosen a1, a2, a3..., obtain a series of point (x1, y1), (x2, y2)(x3, y3) ... it is T to take initial value1, by taken each point successively It substitutes into:yt+1'=y't+a(yt-y't).Meanwhile it is α to take smoothing factor, damped coefficient μ simultaneously changes damped coefficient, changes successively For operation until obtaining final result.
For middle power-assisted area 3, T2<Td<T3, handled using double smoothing.In middle 3 sections of power-assisted area, n point a are chosen1, a2, a3..., obtain a series of point (x1, y1), (x2, y2)(x3, y3) ... it is T to take initial value2, by taken each point successively generation Enter:yt+m=(2+am/(1-a))y't-(1+am/(1-a))yt=(2y't-yt)+m(y't-yt) a/ (1-a) simultaneously, sliding system of making even Number is α, and damped coefficient for μ and changes damped coefficient with this, is iterated operation until obtaining final result.
For big power-assisted area 4, T3<Td<T4, handled using Three-exponential Smoothing.In middle 3 sections of power-assisted area, n point a are chosen1, a2, a3..., obtain a series of point (x1, y1), (x2, y2)(x3, y3) ... it is T to take initial value2, by taken each point successively generation Enter:yt+m=(3y't-3yt)+[(6-5a)y't-(10-8a)yt+(4-3a)yt]*am/2(1-a)2+(y't-2yt+y't)*a2m2/2 (1-a)2Meanwhile it is α to choose smoothing factor, damped coefficient is μ and changes damped coefficient successively, is iterated operation until obtaining Final result.
Step 3) is to 40%<SOC<70% situation carries out double smoothing processing to Curve-type Assistance Characteristics curve;
To 40%<SOC<70% situation carries out single exponential smoothing processing to virgin curve type power-assisted curve, takes n point a1, a2, a3..., obtain a series of point (x1, y1), (x2, y2)(x3, y3) ... it is T to take initial value1, by taken each point successively It substitutes into:yt+m=(2+am/(1-a))y't-(1+am/(1-a))yt=(2y't-yt)+m(y't-yt) a/ (1-a) simultaneously, cunning of making even Coefficient is α, and damped coefficient for μ and changes damped coefficient with this, is iterated operation, until obtaining final result.
Step 4) is to SOC>70% situation carries out single exponential smoothing processing to Curve-type Assistance Characteristics curve;
To SOC>70% situation carries out single exponential smoothing processing to virgin curve type power-assisted curve, takes n point a1, a2, a3..., obtain a series of point (x1, y1), (x2, y2)(x3, y3) ... it is T to take initial value1, taken each point is substituted into successively: yt+1'=y't+a(yt-y't) simultaneously, taking smoothing factor, damped coefficient for μ and changes damped coefficient with this, is iterated fortune for α It calculates, until obtaining final result.
During step 5) is to utilizing hand Calculating Torque during Rotary target current, i.e. Ia=ζ Td, proportionate relationship ζ therein is changed Processing,
Wherein, λ is transformation coefficient, determines the corresponding power-assisted curve of speed;Wherein, IaFor target current, TdFor hand-power square, ζ For hand-power square and the conversion relation of target current;When the hand-power square size for measuring torque sensor is converted into target current, By change conversion coefficient can gentle assist characteristic curve, target current value is made to remain at the discharge current range of battery It is interior, i.e.,Curve can be made to generate the effect of horizontal elongation.
Step 6) is under different speeds, by obtaining the complete assist characteristic under different speeds about the continuous item of speed Curve, i.e.,
I=e-av·ζTd (27)
Wherein I be each vehicle speed condition under target current, e-avFor the compensation coefficient comprising speed, a is constant.
This patent considers the remaining capacity of battery in Curve Design, meanwhile, the design of curve is helped in virgin curve type On the basis of force characteristic curve, in the case where meeting practical application, the calculating of target current can be preferably realized, ensure real Border follow current stablizes output;During power-assisted curve is designed on the one hand be on the basis of shaped form power-assisted curve into Capable smoothing processing, effect are more good;On the other hand, this patent makes the design carried out on the basis of battery capacity is considered, It can ensure the sustainable supply of battery current.
The assist characteristic curve design method of automobile electric booster steering system based on SOC uses double card in step 1) Kalman Filtering theory makes the state of dynamical system the optimal estimation in minimum variance meaning, estimates applied to battery SOC, electricity Pond is seen as dynamical system, and SOC is an internal state of system;The general mathematical form of battery model is:
State equation:xk+1=Akxk+Bkuk+wk=f (xk,uk)+wk (28)
Observational equation:yk=ckxk+vk=g (xk,uk)+vk (29)
Wherein, the input vector u of systemk, to include battery current, temperature, residual capacity and internal resistance variable;
The output y of systemk, the operating voltage for battery;
Battery SOC is included in the quantity of state x of systemkIn, f (xk,uk) and g (xk,uk) be all determined by battery model it is non- Linear equation is linearized in calculating process;SOC algorithms are estimated, including SOC estimation and reaction evaluated error, association side The recursion equation of poor matrix, covariance matrix are used for providing evaluated error range, this equation is in battery model state equation In, SOC is described as to the foundation of state vector:
Specifically double card Kalman Filtering algorithm steps are:
Step (1) gives initial value:
To enable iteration Fast Convergent, by the relatively actual value of the initial value setting of SOC, R0Initial value then lead to It crosses current SOC and tables look-up and provide, other two state can be set as 0, thus just obtain X (0) and R (0);
Step (2) SOC estimates:
The system mode of kth step is estimated using the systematic parameter of -1 step of kth, the system mode then walked again with kth is estimated Count the systematic parameter of kth step:
Step (2.1) is estimated using the system mode that current integration walks kth:
X (k | k-1)=As(k)X(k-1)+Bs(k)I(k)+ωs (31)
Wherein,
Wherein, ωsFor the process noise of system, substantially determined by the noise of electric current;
The optimal estimation of step (2.2) system mode:
X (k)=X (k | k-1)+Ks(k)(V(k)-V(k)') (34)
Wherein V (k) is the battery both end voltage measured, and V (k) ' is the terminal voltage estimated using battery model;
V (k) '=Fv(S(k))-R0(k)I(k)-URC1(k)-URC2(k)+v (35)
Step (2.3) seeks Ks (k) (kalman gain), calculates variance matrix:
Wherein QsFor the covariance matrix of systematic procedure noise, rsThe variance of measurement noise for voltage;Fv(S) it is open circuit Function of the voltage about SOC;
Step (2.4) updates variance matrix P, terminates the Kalman filtering of the SOC of this step;
Ps(k)=[I-Ks(k)Cs(k)]Ps(k|k-1) (39)
It is special to choose shaped form power-assisted for the assist characteristic curve design method of automobile electric booster steering system based on SOC Linearity curve function is:
Wherein,F (v) is the corresponding power-assisted curvilinear system of friction speed;F(v) For the function of speed v, v is faster, and F (v) is smaller, takes the exponential decrease function of v, G (Td) it is when different directions disk torque inputs Power-assisted curve is increasing function, i.e. the input torque of steering wheel is bigger, light in order to turn to, and booster torquemoment also must be bigger.
The assist characteristic curve design method of automobile electric booster steering system based on SOC, used smoothing processing Or exponential smoothing, it is a kind of Time Series Analysis Forecasting method to grow up on the basis of the method for moving average, it is to pass through meter Exponential smoothing value is calculated, cooperation regular hour sequential forecasting models predict the future of phenomenon;Its principle is any phase Exponential smoothing value is all the weighted average of current period actual observation value and previous phase index smooth value, the fundamental formular of exponential smoothing It is:
St=α yt+(1-α)St-1 (41)
In formula, St-- the smooth value of time t;yt-- the actual value of time t;St-1-- the smooth value of time t-1;α -- it is smooth Constant, value range are [0,1];
According to smooth number difference, exponential smoothing is divided into:Single Exponential Smoothing, Secondary Exponential Smoothing Method and index three times Exponential smoothing;
For SOC>70% situation and SOC<Small power-assisted area 2 in the case of 40%, i.e. T1<Td<T2When.Using one Secondary exponential smoothing processing, predictor formula are:
yt+1'=ayt+(1-a)y't (42)
Wherein, yt+1' -- the smooth value S of the predicted value of t+1 phases, i.e. current period (t phases)t (1);yt-- the actual value of t phases;yt'-- The predicted value of t phases, i.e. last smooth value St-1 (1);The formula can be write again:yt+1'=yt'+a(yt-yt'), it is seen then that the next period Predicted value be current period predicted value again with using a as the sum of current period actual value of discount and predicted value error;
For 40%<SOC<70% situation and SOC<Middle power-assisted area 3 in the case of 40%, i.e. T2<Td<T3When.It should It is handled with double smoothing, double smoothing is to the smooth again of single exponential smoothing, it is suitable for tool linear trend Time series, predictor formula are:
yt+m=(2+am/(1-a))y't-(1+am/(1-a))yt=(2y't-yt)+m(y't-yt)a/(1-a) (43)
In formula, yt=ay't-1+(1-a)yt-1
Obviously, double smoothing is a linear equation, and intercept is:(2y't-yt), slope is:(yt'-yt)a/(1- A), independent variable is predicts number of days, double smoothing fundamental formular:
St (2)=aSt (1)+(1-a)St-1 (2) (44)
Yt+ Y=at+btT (45)
at=2St (1)-St (2) (46)
bt=(a/1-a) (St (1)-St (2)) (47)
In formula, St (1)-- the single exponential smoothing value S of t phasest (2)-- the double smoothing value α of t phases -- smoothing factor Yt+T-- t+T phase predicted value T-- elapse issue backward by the t phases;
For SOC<Big power-assisted area 4 in 40% situation, i.e. T3<Td<T4When, it handles using Three-exponential Smoothing, refers to three times Number smoothing predictions be it is secondary it is smooth on the basis of it is smooth again, predictor formula is:
yt+m=(3y't-3yt)+[(6-5a)y't-(10-8a)yt+(4-3a)yt]*am/2(1-a)2+(y't-2yt+y't)* a2m2/2(1-a)2 (48)
In formula, yt=ayt-1+(1-a)yt-1
Predicted value is the weighted sum of former observation, and different data are given with different power, and new data is to larger Power, legacy data give smaller power, by trend adjustment, add trend correction value, can improve exponential smoothing to a certain extent Prediction result, the formula of the exponential smoothing after adjustment are:
Include trend prediction (YtTt)=new prediction (Yt)+trend corrects (Tt)
There are three steps for the Smoothing Prediction of progress trend adjustment:
(1) the simple index number smoothing prediction (Y of t phases is calculated using previously described methodt);
(2) trend is calculated, formula is:
Tt=(1-b) Tt-1+b(Yt-Yt-1) (49)
Wherein, Tt=the t phases smoothed trend;Tt-1=the t phases last smoothed trend;B=selections become Gesture smoothing factor;Yt=to t phase simple index number smoothing predictions;Yt-1=to t phase last issue simple index number smoothing predictions,
(3) the Smoothing Prediction value (Y after trend adjustment is calculatedtTt), calculation formula is:
YtTt=Yt+Tt。 (50)
Engage Fig. 5 it is found that the assist characteristic curve design method of the automobile electric booster steering system based on SOC just It is the situation that battery SOC is considered in 1 design process of the assist characteristic curve, to the virgin curve type in the case of different SOC Assist characteristic curve carries out exponential smoothing processing, even if while being constrained to ensure battery low to electric current needed for EPS system Also there can be steady and sustained electric current to export in the case of electricity, ensure the realization of the entire power-assisted processes of EPS.
The design method process mainly includes:The prediction of battery SOC, carried out according to the value of battery SOC classification smoothing processing, Processing is changed using hand Calculating Torque during Rotary target current process comparative example relationship ζ, the curve under full speed is handled.Its In, the prediction of battery SOC ensures the precision of prediction using double Kalman filtering algorithms;It is flat that classification is carried out according to the value of battery SOC Sliding processing, considers SOC respectively<40%, 40%<SOC<70%, SOC>70% 3 kind of situation, to virgin curve type assist characteristic curve It is smoothed respectively.
The technical solution adopted by the present invention is when being designed to assist characteristic curve, it is contemplated that the situation of SOC, especially It is to be changed greatly in hand-power square and greatly in the case of hand-power square, change the change rate of power-assisted curve, control the electric current of EPS in electricity Within the supply district of pond.
SOC numerical value is observed in real time by Kalman Filter observer, according to SOC numerical values reciteds to virgin curve type assist characteristic Curve carries out a point situation and handles, to SOC<40 situation is segmented virgin curve type assist characteristic curve smoothly to be located respectively Reason, to SOC>40%
Situation with regard to carrying out level-one respectively, two level exponential smoothing is handled, while the demand current of EPS is constrained, protect Lasting supply of the battery to EPS electric currents is demonstrate,proved, so as to be that servo steering system persistently plays a role.
Compared with prior art, the beneficial effects of the invention are as follows battery electricity is considered when assist characteristic curve is designed The situation of SOC is measured, steering boost system can be made to be run in the case of any battery capacity SOC, ensure the steering of driver Safety.

Claims (4)

1. a kind of assist characteristic curve design method of the automobile electric booster steering system based on SOC, it is characterised in that:It is described The design of assist characteristic curve considers the situation of battery SOC, in virgin curve type assist characteristic curve during Curve Design On the basis of be smoothed, even if while being constrained situation to ensure battery in low battery to electric current needed for EPS system Under can also there is steady and sustained electric current to export, ensure the realization of the entire power-assisted processes of EPS, Curve Design includes the following steps:
Step 1) estimates the SOC of battery with kalman filtering theory;SOC is divided into three classes, respectively SOC<40%th, 40%<SOC<70% and SOC>70%;
Step 2) is to SOC<40% situation carries out segment processing corresponding to shaped form assist characteristic curve, is divided into 5 sections 5 Area, respectively:Power-assisted area 1 is power-assisted dead zone, and power-assisted area 2 is small power-assisted area, and power-assisted area 3 is middle power-assisted area, and power-assisted area 4 is to help greatly Power area, power-assisted area 5 are saturation regions;Curve smoothing processing is carried out respectively to 5 sections of assist characteristic curve, and power-assisted area 2 is hand-power square T1<Td<T2, handled using single exponential smoothing;Middle Resistance Level 3 is T2<Td<T3, handled using double smoothing;Big power-assisted area 4 i.e. T3<Td<T4, handled using Three-exponential Smoothing, improve the change rate of shaped form power-assisted curve, so as to stablize the electric current of battery Output;
Step 3) is to 40%<SOC<70% situation carries out double smoothing processing to Curve-type Assistance Characteristics curve;
Step 4) is to SOC>70% situation carries out single exponential smoothing processing to Curve-type Assistance Characteristics curve;
During step 5) is to utilizing hand Calculating Torque during Rotary target current, i.e. Ia=ζ Td, proportionate relationship ζ therein is changed place Reason,
Wherein, λ is transformation coefficient, determines the corresponding power-assisted curve of speed;Wherein, IaFor target current, TdFor hand-power square, ζ is hand The conversion relation of torque and target current;
Step 6) is under different speeds, by obtaining the song of the complete assist characteristic under different speeds about the continuous item of speed Line, i.e.,
I=e-av·ζTd (2)
Wherein I be each vehicle speed condition under target current, e-avFor the compensation coefficient comprising speed, a is constant.
2. the assist characteristic curve design method of the automobile electric booster steering system based on SOC according to claim 1, It is characterized in that:With pair kalman filtering theories the state of dynamical system is made in step 1) optimal in minimum variance meaning Estimation estimates that battery is seen as dynamical system, and SOC is an internal state of system applied to battery SOC;Battery model General mathematical form is:
State equation:xk+1=Akxk+Bkuk+wk=f (xk,uk)+wk (3)
Observational equation:yk=ckxk+vk=g (xk,uk)+vk (4)
Wherein, the input vector u of systemk, to include battery current, temperature, residual capacity and internal resistance variable;
The output y of systemk, the operating voltage for battery;
Battery SOC is included in the quantity of state x of systemkIn, f (xk,uk) and g (xk,uk) be all determined by battery model it is non-linear Equation is linearized in calculating process;SOC algorithms are estimated, including SOC estimation and reaction evaluated error, covariance square The recursion equation of battle array, covariance matrix are used for providing evaluated error range, this equation is in battery model state equation, is incited somebody to action SOC is described as the foundation of state vector:
Specifically double card Kalman Filtering algorithm steps are:
Step (1) gives initial value:
To enable iteration Fast Convergent, by the relatively actual value of the initial value setting of SOC, R0Initial value then by work as Preceding SOC, which tables look-up, to be provided, other two state can be set as 0, has thus just obtained X (0) and R (0);
Step (2) SOC estimates:
The system mode of kth step is estimated using the systematic parameter of -1 step of kth, the system state estimation kth then walked again with kth The systematic parameter of step:
Step (2.1) is estimated using the system mode that current integration walks kth:
X (k | k-1)=As(k)X(k-1)+Bs(k)I(k)+ωs (6)
Wherein,
Wherein, ωsFor the process noise of system, substantially determined by the noise of electric current;
The optimal estimation of step (2.2) system mode:
X (k)=X (k | k-1)+Ks(k)(V(k)-V(k)') (9)
Wherein V (k) is the battery both end voltage measured, and V (k) ' is the terminal voltage estimated using battery model;
V (k) '=Fv(S(k))-R0(k)I(k)-URC1(k)-URC2(k)+v (10)
Step (2.3) seeks Ks (k) (kalman gain), calculates variance matrix:
Wherein QsFor the covariance matrix of systematic procedure noise, rsThe variance of measurement noise for voltage;Fv(S) it is open-circuit voltage Function about SOC;
Step (2.4) updates variance matrix P, terminates the Kalman filtering of the SOC of this step;
Ps(k)=[I-Ks(k)Cs(k)]Ps(k|k-1) (14) 。
3. the assist characteristic curve design method of the automobile electric booster steering system according to claim 1 based on SOC, It is characterized in that:Choosing shaped form assist characteristic curvilinear function is:
Wherein,G(Td)=b1Td 2+b2Td+B3;F (v) is the corresponding power-assisted curvilinear system of friction speed;F (v) is speed v Function, v is faster, and F (v) is smaller, takes the exponential decrease function of v, G (Td) for the input of different directions disk torque when power-assisted it is bent Line is increasing function, i.e. the input torque of steering wheel is bigger, light in order to turn to, and booster torquemoment also must be bigger.
4. the assist characteristic curve design method of the automobile electric booster steering system based on SOC according to claim 1, It is characterized in that:
Smoothing processing or exponential smoothing are a kind of Time Series Analysis Forecastings to grow up on the basis of the method for moving average Method, it is by gauge index smooth value, and cooperation regular hour sequential forecasting models predict the future of phenomenon;It is former It is all the weighted average of current period actual observation value and previous phase index smooth value that reason, which is the exponential smoothing value of any phase, exponential smoothing The fundamental formular of method is:
St=α yt+(1-α)St-1 (16)
In formula, St-- the smooth value of time t;yt-- the actual value of time t;St-1-- the smooth value of time t-1;α -- smoothing constant, Its value range is [0,1];
According to smooth number difference, exponential smoothing is divided into:Single Exponential Smoothing, Secondary Exponential Smoothing Method and Three-exponential Smoothing Method;
For small power-assisted area 2, T1<Td<T2, handled using single exponential smoothing, predictor formula is:
yt+1'=ayt+(1-a)yt' (17)
Wherein, yt+1' -- the smooth value S of the predicted value of t+1 phases, i.e. current period (t phases)t (1);yt-- the actual value of t phases;yt' -- the t phases Predicted value, i.e., last smooth value St-1 (1);The formula can be write again:yt+1'=yt'+a(yt-yt'), it is seen then that the next period is pre- Measured value be current period predicted value again with using a as the sum of current period actual value of discount and predicted value error;
For middle power-assisted area 3, T2<Td<T3, handled using double smoothing, double smoothing is to single exponential smoothing Smooth again, it is suitable for the time series of tool linear trend, and predictor formula is:
yt+m=(2+am/(1-a))yt'-(1+am/(1-a))yt=(2yt'-yt)+m(yt'-yt)a/(1-a) (18)
In formula, yt=ayt'-1+(1-a)yt-1
Obviously, double smoothing is a linear equation, and intercept is:(2yt'-yt), slope is:(yt'-yt) a/ (1-a), from Variable is predicts number of days, double smoothing fundamental formular:
St (2)=aSt (1)+(1-a)St-1 (2) (19)
Yt+ Y=at+btT (20)
at=2St (1)-St (2) (21)
bt=(a/1-a) (St (1)-St (2)) (22)
In formula, St (1)-- the single exponential smoothing value S of t phasest (2)-- the double smoothing value α of t phases -- smoothing factor Yt+T-- t+T phase predicted value T-- elapse issue backward by the t phases;
For big power-assisted area 4, T3<Td<T4, handled using Three-exponential Smoothing, Three-exponential Smoothing prediction is secondary smooth basis On it is smooth again, predictor formula is:
yt+m=(3yt'-3yt)+[(6-5a)yt'-(10-8a)yt+(4-3a)yt]*am/2(1-a)2+(yt'-2yt+yt')*a2m2/ 2(1-a)2 (23)
In formula, yt=ayt-1+(1-a)yt-1
Predicted value is the weighted sum of former observation, and different data are given with different power, and new data gives larger power, old Data give smaller power, by trend adjustment, add trend correction value, can improve Smoothing Prediction knot to a certain extent Fruit, the formula of the exponential smoothing after adjustment are:
Include trend prediction (YtTt)=new prediction (Yt)+trend corrects (Tt)
There are three steps for the Smoothing Prediction of progress trend adjustment:
(1) the simple index number smoothing prediction (Y of t phases is calculated using previously described methodt);
(2) trend is calculated, formula is:
Tt=(1-b) Tt-1+b(Yt-Yt-1) (24)
Wherein, Tt=the t phases smoothed trend;Tt-1=the t phases last smoothed trend;The trend of b=selections is put down Sliding coefficient;Yt=to t phase simple index number smoothing predictions;Yt-1=to t phase last issue simple index number smoothing predictions,
(3) the Smoothing Prediction value (Y after trend adjustment is calculatedtTt), calculation formula is:
YtTt=Yt+Tt。 (25)。
CN201810010887.0A 2018-01-05 2018-01-05 Power-assisted characteristic curve design method of automobile electric power-assisted steering system based on SOC Active CN108197392B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810010887.0A CN108197392B (en) 2018-01-05 2018-01-05 Power-assisted characteristic curve design method of automobile electric power-assisted steering system based on SOC

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810010887.0A CN108197392B (en) 2018-01-05 2018-01-05 Power-assisted characteristic curve design method of automobile electric power-assisted steering system based on SOC

Publications (2)

Publication Number Publication Date
CN108197392A true CN108197392A (en) 2018-06-22
CN108197392B CN108197392B (en) 2021-09-17

Family

ID=62588418

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810010887.0A Active CN108197392B (en) 2018-01-05 2018-01-05 Power-assisted characteristic curve design method of automobile electric power-assisted steering system based on SOC

Country Status (1)

Country Link
CN (1) CN108197392B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110895332A (en) * 2019-12-03 2020-03-20 电子科技大学 Distributed tracking method for extended target
CN111006885A (en) * 2019-12-12 2020-04-14 北京经纬恒润科技有限公司 Method and device for calibrating EPS (electric power steering) model parameters of electric power steering system
CN117290965A (en) * 2023-11-22 2023-12-26 中汽研汽车检验中心(广州)有限公司 Vehicle model simulation test simulation method, equipment and medium of vehicle simulation software

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201230055Y (en) * 2008-07-21 2009-04-29 武汉亚格光电技术有限公司 Self-locking water, electricity path combined rapid coupler
CN101722980A (en) * 2008-10-31 2010-06-09 上海格尔汽车附件有限公司 Electric power steering system with individually adjustable-parameter and curve-type power characteristic
JP2015229359A (en) * 2014-06-03 2015-12-21 株式会社ジェイテクト Electric power steering device
CN106184351A (en) * 2016-07-11 2016-12-07 南京航空航天大学 A kind of Multipurpose Optimal Method of electric-hydraulic combined power steering system
CN106240390A (en) * 2016-08-09 2016-12-21 潍柴动力股份有限公司 A kind of method of dynamic optimization electric energy when power system for pure electric bus and low SOC thereof

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201230055Y (en) * 2008-07-21 2009-04-29 武汉亚格光电技术有限公司 Self-locking water, electricity path combined rapid coupler
CN101722980A (en) * 2008-10-31 2010-06-09 上海格尔汽车附件有限公司 Electric power steering system with individually adjustable-parameter and curve-type power characteristic
JP2015229359A (en) * 2014-06-03 2015-12-21 株式会社ジェイテクト Electric power steering device
CN106184351A (en) * 2016-07-11 2016-12-07 南京航空航天大学 A kind of Multipurpose Optimal Method of electric-hydraulic combined power steering system
CN106240390A (en) * 2016-08-09 2016-12-21 潍柴动力股份有限公司 A kind of method of dynamic optimization electric energy when power system for pure electric bus and low SOC thereof

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ZHAO XUE-PING等: "Parametric design and application of steering characteristic curve in control for electric power steering", 《MECHATRONICS》 *
王金归: "电动助力转向控制系统研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》 *
黄玉鹏: "基于dSPACE的电动助力转向系统控制策略研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110895332A (en) * 2019-12-03 2020-03-20 电子科技大学 Distributed tracking method for extended target
CN111006885A (en) * 2019-12-12 2020-04-14 北京经纬恒润科技有限公司 Method and device for calibrating EPS (electric power steering) model parameters of electric power steering system
CN111006885B (en) * 2019-12-12 2021-11-23 北京经纬恒润科技股份有限公司 Method and device for calibrating EPS (electric power steering) model parameters of electric power steering system
CN117290965A (en) * 2023-11-22 2023-12-26 中汽研汽车检验中心(广州)有限公司 Vehicle model simulation test simulation method, equipment and medium of vehicle simulation software
CN117290965B (en) * 2023-11-22 2024-04-09 中汽研汽车检验中心(广州)有限公司 Vehicle model simulation test simulation method, equipment and medium of vehicle simulation software

Also Published As

Publication number Publication date
CN108197392B (en) 2021-09-17

Similar Documents

Publication Publication Date Title
US11958365B2 (en) Method for dual-motor control on electric vehicle based on adaptive dynamic programming
CN103770832A (en) Assisted power controlling method for electric power assisted steering system
CN108197392A (en) The assist characteristic curve design method of automobile electric booster steering system based on SOC
CN101898582B (en) Method of returning steering wheel using motor
CN107207040B (en) Electric power steering apparatus
CN110739893B (en) Improved self-adaptive trackless Kalman filtering rotational inertia identification method
CN109039167B (en) Control method and system for built-in permanent magnet synchronous motor
CN107154763A (en) Permagnetic synchronous motor dead beat direct Torque Control and control method
CN111608868B (en) Maximum power tracking adaptive robust control system and method for wind power generation system
CN109507876B (en) Credibility reasoning-based PID parameter setting method for electric propulsion ship motor
CN107171612A (en) Fuzzy score rank PID switched reluctance machines method for controlling torque and system
CN104908814A (en) Fractional order PID control method for automobile steer-by-wire system
CN111290263B (en) Improved PID (proportion integration differentiation) optimization control algorithm based on RBFNN (radial basis function network) and BAS (basic automatic component analysis)
CN108563845B (en) A kind of building method of the EPS controller based on composite power source
CN104122795A (en) Novel extremal function index based intelligent self-turning PID (Proportion Integration Differentiation) indoor temperature control algorithm
CN108054972A (en) A kind of method for improving permanent magnetic linear synchronous motor dynamic control performance
CN101639665A (en) Control method of humanoid robot
CN109334757B (en) Control method of electric power steering system
CN104167968A (en) Asynchronous motor vector control method
CN103777521B (en) A kind of low speed control method of vehicle based on fuzzy control
CN108073554B (en) Power-assisted characteristic curve design method of automobile electric power-assisted steering system based on SOC
CN105292243B (en) The forecast Control Algorithm of automobile electric booster steering system
CN104270046A (en) Motor control method based on self-learning of rotating speed-current two-dimensional fuzzy model
CN114670856B (en) Parameter self-tuning longitudinal control method and system based on BP neural network
EP4166782A1 (en) Wind turbine model predictive control (mpc) with thrust and/or power correction

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
GR01 Patent grant
GR01 Patent grant