CN113511211A - Torsional vibration control method based on electric driving system of electric vehicle - Google Patents
Torsional vibration control method based on electric driving system of electric vehicle Download PDFInfo
- Publication number
- CN113511211A CN113511211A CN202110604573.5A CN202110604573A CN113511211A CN 113511211 A CN113511211 A CN 113511211A CN 202110604573 A CN202110604573 A CN 202110604573A CN 113511211 A CN113511211 A CN 113511211A
- Authority
- CN
- China
- Prior art keywords
- motor
- actual
- torque
- rotating speed
- vehicle
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000012937 correction Methods 0.000 claims abstract description 17
- 238000011156 evaluation Methods 0.000 claims description 27
- 230000001133 acceleration Effects 0.000 claims description 12
- UONOETXJSWQNOL-UHFFFAOYSA-N tungsten carbide Chemical compound [W+]#[C-] UONOETXJSWQNOL-UHFFFAOYSA-N 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 9
- 238000012887 quadratic function Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims 1
- 230000010355 oscillation Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 239000003638 chemical reducing agent Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W30/00—Purposes of road vehicle drive control systems not related to the control of a particular sub-unit, e.g. of systems using conjoint control of vehicle sub-units
- B60W30/18—Propelling the vehicle
- B60W30/20—Reducing vibrations in the driveline
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W2510/00—Input parameters relating to a particular sub-units
- B60W2510/08—Electric propulsion units
- B60W2510/081—Speed
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W2510/00—Input parameters relating to a particular sub-units
- B60W2510/08—Electric propulsion units
- B60W2510/083—Torque
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W2552/00—Input parameters relating to infrastructure
- B60W2552/15—Road slope, i.e. the inclination of a road segment in the longitudinal direction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/72—Electric energy management in electromobility
Landscapes
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Control Of Electric Motors In General (AREA)
Abstract
The invention discloses a torsional vibration control method based on an electric automobile electric driving system, which obtains torsional vibration compensation torque by multiplying the difference value between the model rotating speed and the actual rotating speed of a motor by a torsional vibration compensation coefficient, and further corrects the target torque of the motor. The motor model rotating speed is the theoretical rotating speed of the motor without torsional vibration influence, the actual rotating speed of the motor is compared with the rotating speed of the motor model, and then the actual rotating speed of the motor is used as a correction basis to be close to the rotating speed of the motor model as far as possible, the correction basis is accurate, the torque correction is timely, and the driving smoothness of the whole vehicle can be improved.
Description
Technical Field
The invention belongs to the field of electric vehicle control, and particularly relates to a torsional vibration control method based on an electric drive system of an electric vehicle.
Background
The electric driving system of the electric automobile can be simplified into a single-mass spring system when the automobile runs dynamically. When the running condition of the vehicle is changed violently, the rotating speed oscillation impact can be caused, the vehicle shakes forwards and backwards, and torsional vibration is generated. Torsional vibration control is the reduction or attenuation of such oscillations by torque intervention to the point where they are imperceptible to, or within acceptable limits for, the driver.
CN106740272A discloses a pure electric vehicle and a method and a device for controlling low-speed shaking of the pure electric vehicle, wherein whether the pure electric vehicle shakes at a low speed is judged according to the rotating speed of a motor output shaft; when the pure electric vehicle shakes at a low speed, the rotating speed of the output shaft of the motor is controlled to be maintained within a preset range, and therefore the low-speed shaking of the pure electric vehicle is better inhibited; the method specifically comprises the following steps: and determining the output torque change amount corresponding to the absolute value of the difference between the currently acquired rotating speed of the motor output shaft and the upper limit value or the lower limit value of the preset range according to the prestored corresponding relation between the rotating speed change amount of the motor output shaft and the output torque change amount of the motor, and reducing or increasing the output torque of the motor by the output torque change amount so as to effectively adjust the rotating speed of the motor output shaft within the preset range. The method is characterized in that torque correction is carried out when the pure electric vehicle shakes at a low speed, so that the rotating speed of the output shaft of the motor is maintained in a preset range, delay exists in torque correction, the rotating speed of the output shaft of the motor is maintained in the preset range, and a driver still can feel irregularity.
Disclosure of Invention
The invention aims to provide a torsional vibration control method based on an electric driving system of an electric vehicle, so that torsional vibration is effectively controlled, and the driving smoothness of the whole vehicle is improved.
The invention discloses a torsional vibration control method based on an electric vehicle electric driving system, which comprises the following steps:
read motor model rotational speed MotSpdMod of moment t(t)Actual rotation speed MotSpd of motor(t)Actual torque M of motor(t)Pressure P of brake master cylinder(t)And the slope of the road on which the vehicle is located(t);
According to said slope(t)And the actual rotational speed MotSpd of the motor(t)Determining the running resistance at t moment and converting the running resistance into the running resistance moment WR at the motor end(t);
According to the pressure P of the brake master cylinder(t)And the actual rotational speed MotSpd of the motor(t)Determining the brake master cylinder pressure at the time t and converting the pressure into the brake torque WB at the motor end(t);
Using the integral formula of the rotation speed: MotSpdMod(t+dt)=MotSpdMod(t)+(M(t)-WR(t)-WB(t)) XKf x dt, calculating the motor model rotating speed MotSpdMod at the moment of t + dt(t+dt)And storing; kf is a preset power assembly characteristic correction coefficient, the vehicle is static at the initial moment, the initial value of the rotating speed of the motor model is equal to zero, the initial value of the actual rotating speed of the motor is equal to zero, the initial value of the actual torque of the motor is equal to zero, and the sum of the initial value of the driving resistance torque converted to the motor end by the driving resistance and the initial value of the braking torque converted to the motor end by the pressure of the brake master cylinder is equal to zero;
obtaining the actual rotating speed MotSpd of the motor at the moment of t + dt(t+dt)Actual torque M of motor(t+dt)Pressure P of brake master cylinder(t+dt)And the slope of the road on which the vehicle is located(t+dt)And storing;
according to the actual rotating speed MotSpd of the motor(t+dt)Actual torque M of motor(t+dt)Inquiring a rotating speed-torque-torsional vibration resistance compensation coefficient table to obtain a corresponding torsional vibration resistance compensation coefficient k(t+dt)(ii) a The rotating speed-torque-torsional vibration compensation coefficient table is a pre-stored corresponding relation table of the actual rotating speed of the motor, the actual torque of the motor and the torsional vibration compensation coefficient;
using the formula for torsional vibration compensation torque: TC (tungsten carbide)(t+dt)=k(t+dt)×(MotSpdMod(t+dt)-MotSpd(t+dt)) And calculating to obtain the torsional vibration compensation torque TC at the moment of t + dt(t+dt);
Compensating torque TC for torsional vibration at the moment of t + dt(t+dt)And superposing the target torque with the original target torque of the motor to obtain the corrected target torque of the motor at the moment of t + dt.
Preferably, according to said slope(t)And the actual rotational speed MotSpd of the motor(t)Determining the driving resistance at time t and converting into the driving resistance moment WR at the motor end(t)The specific mode is as follows:
according to said slope(t)Inquiring a gradient-fitting parameter table to obtain a specific numerical value of the corresponding fitting parameter A, B, C; wherein, the gradient-fitting parameter table is a pre-stored corresponding relation table of the gradient of the road where the vehicle is located and the fitting parameter A, B, C;
the corresponding specific value of the fitting parameter A, B, C and the actual rotating speed MotSpd of the motor(t)And carrying in a running resistance moment calculation formula: WR (pulse Width modulation)(t)=A×MotSpd(t) 2+B×MotSpd(t)+ C, calculating the running resistance at t moment and converting the running resistance into the running resistance moment WR at the motor end(t)。
Preferably, the master cylinder pressure P is based on(t)And the actual rotational speed MotSpd of the motor(t)Determining the brake master cylinder pressure at the time t and converting the pressure into the brake torque WB at the motor end(t)The specific mode is as follows:
if the actual rotating speed of the motor MoySpd(t)Is less than or equal to the preset rotation speed threshold Mthre, the vehicle is stationary, and WB is enabled(t)=-WR(t);
If the actual speed MotSpd of the motor(t)If the absolute value of (b) is greater than the preset rotation speed threshold value Mthre, it indicates that the vehicle is moving, and WB is made(t)=sign×Kb×P(t);
Wherein, when MotSpd(t)When > 0, sign is 1, when MotSpd(t)When < 0, sign is-1, KbThe method is a preset conversion coefficient for converting the pressure of the brake master cylinder to the braking torque of the motor end.
Preferably, the slope-fitting parameter table is obtained by the following test:
firstly, enabling the vehicle to run at a constant speed at different speeds on the same slope, and recording the actual motor torque M at t moment corresponding to different speeds on the slope after the actual motor torque and the actual motor rotating speed are stable(t)And the actual rotating speed MotSpd of the motor(t)(ii) a Wherein, under the condition of constant-speed running of the vehicle, the running resistance at the time t is converted into the running resistance moment WR of the motor end(t)Actual torque M of motor equal to t moment(t);
Secondly, the actual rotating speed MotSpd of the motor at the t moment corresponding to different vehicle speeds under the slope(t)Converted to a driving resistance torque WR at the motor end from the driving resistance(t)And fitting a quadratic function to obtain a running resisting moment calculation formula: WR (pulse Width modulation)(t)=A×MoySpd(t) 2+B×MotSpd(t)+ C; then extracting fitting parameters A, B, C to obtain a group of fitting parameters A, B, C;
step three, changing the gradient for n-1 times, and then repeatedly executing the step one and the step two until n groups of fitting parameters A, B, C are obtained;
and fourthly, corresponding n groups of fitting parameters A, B, C to n gradients one by one to form the gradient-fitting parameter table.
Preferably, the preset powertrain characteristic correction coefficient Kf is obtained by the following test:
the vehicle is enabled to run on a flat road at uniform acceleration, and the actual torque of the motor at m moments in the process of running at uniform acceleration is recordedActual rotating speed of motor at m moments Wherein, under the condition of uniform acceleration running of the vehicle, the motor model rotating speed at m moments Respectively equal to the actual rotational speed of the motor at the corresponding moment The pressure of the brake master cylinder at m moments is converted into the braking torque of the motor endAre all equal to zero;
fitting parameters A, B, C corresponding to the flat road and the actual rotating speed of the motor at m moments Respectively carrying into a driving resistance moment calculation formula: WR (pulse Width modulation)(t)=A×MotSpd(t) 2+B×MotSpd(t)+ C, calculating the running resistance at m times and converting the running resistance into the running resistance torque at the motor end
Respectively and correspondingly substituting the actual torque of the motor at m moments, the rotating speed of the motor model at m moments, the pressure of the brake master cylinder at m moments into the braking torque at the motor end and the driving resistance at m moments into the driving resistance torque at the motor end into formulas: MotSpdMod(t+dt)=MotSpdMod(t)+(M(t)-WR(t)-WB(t))×Kf1X dt, m Kf are calculated1;
For m Kf1Averaging, m Kf1Is used as the preset power assembly characteristic correction coefficientKf。
Preferably, the table of the rotational speed-torque-torsional vibration compensation coefficients is obtained by the following test:
firstly, presetting a torsional vibration compensation coefficient as an empirical value or zero;
secondly, enabling the vehicle to run at the actual rotating speed and the actual torque of a certain motor; if the vehicle meets the driving smoothness evaluation standard (for example, no torsional vibration occurs or the torsional vibration is within an acceptable range), recording the actual rotating speed of the motor, the actual torque of the motor and the torsional vibration resistance compensation coefficient at the moment, if the vehicle does not meet the driving smoothness evaluation standard, changing the torsional vibration resistance compensation coefficient, enabling the vehicle to continuously run at the actual rotating speed of the motor and the actual torque of the motor until the vehicle meets the driving smoothness evaluation standard, and recording the actual rotating speed of the motor, the actual torque of the motor and the torsional vibration resistance compensation coefficient when the vehicle meets the driving smoothness evaluation standard;
thirdly, changing the actual rotating speed and/or the actual torque of the motor, enabling the vehicle to run at the changed actual rotating speed and/or the changed actual torque of the motor, recording the actual rotating speed, the actual torque and the torsional vibration resistance compensation coefficient of the motor at the moment if the vehicle meets the driving smoothness evaluation standard, changing the torsional vibration resistance compensation coefficient if the vehicle does not meet the driving smoothness evaluation standard, enabling the vehicle to continue running at the actual rotating speed and the actual torque of the motor until the vehicle meets the driving smoothness evaluation standard, and recording the actual rotating speed, the actual torque and the torsional vibration resistance compensation coefficient of the motor when the vehicle meets the driving smoothness evaluation standard;
and fourthly, combining all the actual rotating speeds of the motor, the actual torque of the motor and the torsional vibration resistance compensation coefficient which meet the driving smoothness evaluation standard to form the rotating speed-torque-torsional vibration resistance compensation coefficient table.
According to the invention, the torsional vibration compensation torque is obtained by multiplying the difference value between the model rotating speed and the actual rotating speed of the motor by the torsional vibration compensation coefficient, and the target torque of the motor is corrected. The motor model rotating speed is the theoretical rotating speed of the motor without torsional vibration influence, the actual rotating speed of the motor is compared with the rotating speed of the motor model, and then the actual rotating speed of the motor is used as a correction basis to be close to the rotating speed of the motor model as far as possible, the correction basis is accurate, and the torque correction is timely, so that the torsional vibration is more effectively controlled, and the driving smoothness of the whole vehicle is improved.
Drawings
Fig. 1 is a schematic diagram of torsional vibration control based on an electric drive system of an electric vehicle according to the present embodiment.
Fig. 2 is a flowchart of a torsional vibration control method based on an electric drive system of an electric vehicle according to the present embodiment.
Detailed Description
As shown in fig. 1 and fig. 2, the torsional vibration control method based on the electric drive system of the electric vehicle in the present embodiment includes:
step one, reading the stored motor model rotating speed MotSpdMod at the moment t(t)Actual rotation speed MotSpd of motor(t)Actual torque M of motor(t)Pressure P of brake master cylinder(t)And the slope of the road on which the vehicle is located(t)。
Step two, slope according to t moment(t)And the actual motor speed MotSpd at time t(t)Determining the running resistance at t moment and converting the running resistance into the running resistance moment WR at the motor end(t)。
First, the slope according to time t(t)And inquiring the gradient-fitting parameter table to obtain the specific numerical value of the corresponding fitting parameter A, B, C. Wherein, the gradient-fitting parameter table is a pre-stored corresponding relation table of the gradient of the road where the vehicle is located and the fitting parameter A, B, C;
then, the corresponding specific value of the fitting parameter A, B, C and the actual motor rotation speed MotSpd at the time t are compared(t)And carrying in a running resistance moment calculation formula: WR (pulse Width modulation)(t)=A×MotSpd(t) 2+B×MotSpd(t)+ C, calculating the running resistance at t moment and converting the running resistance into the running resistance moment WR at the motor end(t)。
As an example, the slope-fitting parameter table described above can be obtained by the following tests:
firstly, enabling the vehicle to run at a constant speed at different speeds under the same gradient, and waiting for the actual torque of the motorAfter the actual rotating speed of the motor is stable, recording the actual torque M of the motor at the t moment corresponding to different vehicle speeds under the slope(t)And the actual rotating speed MotSpd of the motor(t). Wherein, under the condition of constant-speed running of the vehicle, the running resistance at the time t is converted into the running resistance moment WR of the motor end(t)Actual torque M of motor equal to t moment(t)。
Secondly, the actual rotating speed MotSpd of the motor at the t moment corresponding to different vehicle speeds under the slope(t)Converted to a driving resistance torque WR at the motor end from the driving resistance(t)And fitting a quadratic function to obtain a running resisting moment calculation formula: WR (pulse Width modulation)(t)=A×MotSpd(t) 2+B×MotSpd(t)+ C; the fitting parameters A, B, C are then extracted to obtain a set of fitting parameters A, B, C. The vehicle running resistance includes rolling resistance, air resistance, and gradient resistance. The difference in driving force and running resistance causes the vehicle to accelerate. On the premise that the calculation period dt is small enough, the output torque of the motor is approximately considered to be unchanged, and the theoretical rotating speed (namely the rotating speed of the motor model) of the next sampling period of the motor can be estimated if the running resistance is known. The air resistance of the vehicle is proportional to the square of the vehicle speed, and the running resistance is converted into the resistance moment of the motor end by fitting a quadratic function of the actual rotating speed of the motor.
And thirdly, changing the gradient n-1 times, and then repeatedly executing the first step and the second step until n groups of fitting parameters A, B, C are obtained.
And fourthly, corresponding n groups of fitting parameters A, B, C to n gradients one by one to form the gradient-fitting parameter table.
Step three, according to the pressure P of the brake master cylinder at the time t(t)And the actual motor speed MotSpd at time t(t)Determining the brake master cylinder pressure at the time t and converting the pressure into the brake torque WB at the motor end(t)。
Judging whether the actual rotating speed MotSpd of the motor exists or not(t)The absolute value of (a) is less than or equal to a preset rotation speed threshold value Mthre;
if yes, the vehicle is shown to be stationary, and WB is driven(t)=-WR(t)。
Otherwise (i.e. actual speed of the motor)MotSpd(t)Is greater than a preset rotational speed threshold value Mthre), indicating that the vehicle is moving, and making WB(t)=sign×Kb×P(t)。
Wherein Mthre is more than or equal to 0; when MotSpd(t)When > 0, sign is 1, when MotSpd(t)When < 0, sign is-1, KbFor a predetermined conversion factor, K, for converting the pressure of the brake master cylinder to the braking torque at the motor endbRelated to braking system, tire radius, speed ratio of the retarder, etc.
Step four, utilizing a rotating speed integral formula: MotSpdMod(t+dt)=MotSpdMod(t)+(M(t)-WR(t)-WB(t)) XKf x dt, calculating the motor model rotating speed MotSpdMod at the moment of t + dt(t+dt)And stored. Kf is a preset power assembly characteristic correction coefficient; the vehicle is static at the initial moment, and the initial value of the rotating speed of the motor model (namely the rotating speed MotSpdMod of the motor model at the initial moment)(0)) Equal to zero, the initial value of the actual rotation speed of the motor (i.e. the actual rotation speed MotSpd of the motor at the initial moment)(0)) Equal to zero, the initial value of the actual torque of the motor (i.e. the actual torque M of the motor at the initial moment)(0)) Is equal to zero; initial value of master cylinder pressure (i.e. master cylinder pressure P at initial time)(0)) Satisfies the following conditions: the running resistance is converted to an initial value of the running resistance torque of the motor end (i.e., the running resistance at the initial time is converted to the running resistance torque WR of the motor end(0)) The initial value of the braking torque converted to the motor end from the pressure of the brake master cylinder (namely the braking torque WB converted to the motor end from the pressure of the brake master cylinder at the initial moment)(0)) The sum is equal to zero; i.e. WR if the vehicle is stationary on a flat road with a gradient of 0 at the initial moment(0)=0,WB(0)0; WR if the vehicle is stationary on a grade other than 0 at the initial time(0)≠0,WB(0)=-WR(0)The braking torque maintains the vehicle on the grade.
The vehicle dynamics model simplified equation is:wherein M is(t)At time tActual torque (i.e. motor-side output torque), WR, of the motor(t)Converting the driving resistance at time t into the driving resistance torque at the motor end, WB(t)The pressure of the brake master cylinder at the time t is converted into the brake torque at the end of the motor, J is the rotational inertia of the whole driving and transmission chain,is the angular acceleration of the motor model.
According to the following steps:(i.e. theIs proportional to) Andand (4) obtaining a rotation speed integral formula by sorting: MotSpdMod(t+dt)=MotSpdMod(t)+(M(t)-WR(t)-WB(t))×Kf×dt。
Wherein the content of the first and second substances,representing the motor model rotational acceleration, MotSpdMod(t+dt)Motor model speed, MotSpdMod, representing time t + dt(t)Expressing the rotating speed of the motor model at the moment t, Kf is a preset power assembly characteristic correction coefficient,(i.e., Kf is proportional to). Kf depends on the characteristics of the power assembly, such as the gear of a speed reducer, the vehicle servicing mass and the rotational inertia of a shafting of a motor system.
By way of example, Kf may be obtained by the following test:
firstly, making the vehicle carry out uniform acceleration running on a flat road (the gradient is 0), and recording the actual torque of the motor at m moments in the uniform acceleration running processActual rotating speed of motor at m momentsWherein, under the condition of uniform acceleration running of the vehicle, the motor model rotating speed at m moments Respectively equal to the actual rotational speed of the motor at the corresponding moment The pressure of the brake master cylinder at m moments is converted into the braking torque of the motor endAre all equal to zero (since the brake pedal is not stepped on when the vehicle is in uniform acceleration running on a flat road, the brake master cylinder pressure at each moment is equal to zero, and therefore the brake master cylinder pressure at each moment is converted into the brake torque at the motor end and is also equal to zero).
Secondly, fitting parameters A, B, C corresponding to the flat road and the actual rotating speed of the motor at m moments Respectively carrying into a driving resistance moment calculation formula: WR (pulse Width modulation)(t)=A×MotSpd(t) 2+B×MotSpd(t)+ C, calculating the running resistance at m timesForce conversion to motor-end driving moment
Thirdly, the actual torque of the motor at m momentsmotor model rotation speed at m moments The pressure of the brake master cylinder at m times is converted into the braking torque (all equal to zero) of the motor end, and the running resistance at m times is converted into the running resistance torque of the motor end Respectively correspondingly substituting into the formulas: MotSpdMod(t+dt)=MotSpdMod(t)+(M(t)-WR(t)-WB(t))×Kf1X dt, m Kf are calculated1。
For m Kf1Averaging, m Kf1The average value of (d) is used as a preset powertrain characteristic correction coefficient Kf.
Step five, acquiring the actual rotating speed MotSpd of the motor at the moment of t + dt(t+dt)Actual torque M of the motor at time t + dt(t+dt)Brake master cylinder pressure P at time t + dt(t+dt)And the slope of the road on which the vehicle is located at time t + dt(t+dt)And stored.
Step six, according to the actual rotating speed MotSpd of the motor at the moment of t + dt(t+dt)Actual torque M of motor(t+dt)Inquiring a rotating speed-torque-torsional vibration resistance compensation coefficient table to obtain a corresponding torsional vibration resistance compensation coefficient k(t+dt). Wherein, the rotating speed-torque-anti-torsional vibration compensation coefficient table is a prestored motor actual rotation tableAnd the corresponding relation table of the speed, the actual torque of the motor and the torsional vibration compensation coefficient.
As an example, the table of the rotational speed-torque-torsional vibration compensation coefficients can be obtained by the following tests:
in the first step, the torsional vibration compensation coefficient is preset to an empirical value or zero.
Secondly, enabling the vehicle to run at the actual rotating speed and the actual torque of a certain motor; if the vehicle meets the driving smoothness evaluation standard (for example, no torsional vibration occurs or the torsional vibration is within an acceptable range), recording the actual rotating speed of the motor, the actual torque of the motor and the torsional vibration compensation coefficient at the moment, if the vehicle does not meet the driving smoothness evaluation standard, changing the torsional vibration compensation coefficient, enabling the vehicle to continuously run at the actual rotating speed of the motor and the actual torque of the motor until the vehicle meets the driving smoothness evaluation standard, and recording the actual rotating speed of the motor, the actual torque of the motor and the torsional vibration compensation coefficient when the vehicle meets the driving smoothness evaluation standard.
And thirdly, changing the actual rotating speed and/or the actual torque of the motor, enabling the vehicle to run at the changed actual rotating speed and/or the changed actual torque of the motor, recording the actual rotating speed, the actual torque and the torsional vibration resistance compensation coefficient of the motor at the moment if the vehicle meets the driving smoothness evaluation standard, changing the torsional vibration resistance compensation coefficient if the vehicle does not meet the driving smoothness evaluation standard, enabling the vehicle to continue running at the actual rotating speed and the actual torque of the motor until the vehicle meets the driving smoothness evaluation standard, and recording the actual rotating speed, the actual torque and the torsional vibration resistance compensation coefficient of the motor when the vehicle meets the driving smoothness evaluation standard.
And fourthly, combining all the actual rotating speeds of the motor, the actual torque of the motor and the torsional vibration resistance compensation coefficient which meet the driving smoothness evaluation standard to form a rotating speed-torque-torsional vibration resistance compensation coefficient table.
Step seven, utilizing a torsional vibration compensation torque formula: TC (tungsten carbide)(t+dt)=k(t+dt)×(MotSpdMod(t+dt)-MoySpd(t+dt)) And calculating to obtain the torsional vibration compensation torque TC at the moment of t + dt(t+dt)。
Step eight, compensating the torsional vibration resistance torque TC at the time of t + dt(t+dt)And superposing the target torque with the original target torque of the motor to obtain the corrected target torque of the motor at the moment of t + dt.
And outputting the corrected motor target torque to a downstream control module for torque control, so that the driving smoothness of the whole vehicle is improved.
Claims (6)
1. A torsional vibration control method based on an electric drive system of an electric vehicle is characterized by comprising the following steps:
read motor model rotational speed MotSpdMod of moment t(t)Actual rotation speed MotSpd of motor(t)Actual torque M of motor(t)Pressure P of brake master cylinder(t)And the slope of the road on which the vehicle is located(t);
According to said slope(t)And the actual rotational speed MotSpd of the motor(t)Determining the running resistance at t moment and converting the running resistance into the running resistance moment WR at the motor end(t);
According to the pressure P of the brake master cylinder(t)And the actual rotational speed MotSpd of the motor(t)Determining the brake master cylinder pressure at the time t and converting the pressure into the brake torque WB at the motor end(t);
Using the integral formula of the rotation speed: MotSpdMod(t+dt)=MotSpdMod(t)+(M(t)-WR(t)-WB(t)) XKf x dt, calculating the motor model rotating speed MotSpdMod at the moment of t + dt(t+dt)And storing; kf is a preset power assembly characteristic correction coefficient, the vehicle is static at the initial moment, the initial value of the rotating speed of the motor model is equal to zero, the initial value of the actual rotating speed of the motor is equal to zero, the initial value of the actual torque of the motor is equal to zero, and the sum of the initial value of the driving resistance torque converted to the motor end by the driving resistance and the initial value of the braking torque converted to the motor end by the pressure of the brake master cylinder is equal to zero;
obtaining the actual rotating speed MotSpd of the motor at the moment of t + dt(t+dt)Actual torque M of motor(t+dt)Pressure P of brake master cylinder(t+dt)And slope of road on which vehicle is located(t+dt)And storing;
according to the actual rotating speed MotSpd of the motor(t+dt)Actual torque M of motor(t+dt)Inquiring a rotating speed-torque-torsional vibration resistance compensation coefficient table to obtain a corresponding torsional vibration resistance compensation coefficient k(t+dt)(ii) a The rotating speed-torque-torsional vibration compensation coefficient table is a pre-stored corresponding relation table of the actual rotating speed of the motor, the actual torque of the motor and the torsional vibration compensation coefficient;
using the formula for torsional vibration compensation torque: TC (tungsten carbide)(t+dt)=k(t+dt)×(MDtSpdMod(t+dt)-MotSpd(t+dt)) And calculating to obtain the torsional vibration compensation torque TC at the moment of t + dt(t+dt);
Compensating torque TC for torsional vibration at the moment of t + dt(t+dt)And superposing the target torque with the original target torque of the motor to obtain the corrected target torque of the motor at the moment of t + dt.
2. Method for torsional vibration control based on an electric vehicle electric drive system according to claim 1, characterized in that the gradient slope is dependent on(t)And the actual rotational speed MotSpd of the motor(t)Determining the running resistance at t moment and converting the running resistance into the running resistance moment WR at the motor end(t)The specific mode is as follows:
according to said slope(t)Inquiring a gradient-fitting parameter table to obtain a specific numerical value of the corresponding fitting parameter A, B, C; wherein, the gradient-fitting parameter table is a pre-stored corresponding relation table of the gradient of the road where the vehicle is located and the fitting parameter A, B, C;
the corresponding specific value of the fitting parameter A, B, C and the actual rotating speed MotSpd of the motor(t)And carrying in a running resistance moment calculation formula: WR (pulse Width modulation)(t)=A×MotSpd(t) 2+B×MotSpd(t)+ C, calculating the running resistance at t moment and converting the running resistance into the running resistance moment WR at the motor end(t)。
3. Method for controlling torsional vibrations based on an electric drive system of an electric vehicle according to claim 1 or 2, characterized in that it is based on the brake master cylinder pressure P(t)And the actual rotational speed MotSpd of the motor(t)Determining the brake master cylinder pressure at the time t and converting the pressure into the brake torque WB at the motor end(t)The specific mode is as follows:
if the actual speed MotSpd of the motor(t)Is less than or equal to the preset rotation speed threshold Mthre, the vehicle is stationary, and WB is enabled(t)=-WR(t);
If the actual speed MotSpd of the motor(t)If the absolute value of (b) is greater than the preset rotation speed threshold value Mthre, it indicates that the vehicle is moving, and WB is made(t)=sign×Kb×P(t);
Wherein, when MotSpd(t)When > 0, sign is 1, when MotSpd(t)When < 0, sign is-1, KbThe method is a preset conversion coefficient for converting the pressure of the brake master cylinder to the braking torque of the motor end.
4. The method of claim 2 or 3, wherein the slope-fit parameter table is obtained by testing:
firstly, enabling the vehicle to run at a constant speed at different speeds on the same slope, and recording the actual motor torque M at t moment corresponding to different speeds on the slope after the actual motor torque and the actual motor rotating speed are stable(t)And the actual rotating speed MotSpd of the motor(t)(ii) a Wherein, under the condition of constant-speed running of the vehicle, the running resistance at the time t is converted into the running resistance moment WR of the motor end(t)Actual torque M of motor equal to t moment(t);
Secondly, the actual rotating speed MotSpd of the motor at the t moment corresponding to different vehicle speeds under the slope(t)Converted to a driving resistance torque WR at the motor end from the driving resistance(t)And fitting a quadratic function to obtain a running resisting moment calculation formula: WR (pulse Width modulation)(t)=A×MotSpd(t) 2+B×MotSpd(t)+ C; then extracting fitting parameters A, B, C to obtain a group of fitting parameters A, B, C;
step three, changing the gradient for n-1 times, and then repeatedly executing the step one and the step two until n groups of fitting parameters A, B, C are obtained;
and fourthly, corresponding n groups of fitting parameters A, B, C to n gradients one by one to form the gradient-fitting parameter table.
5. The method of claim 4, wherein the predetermined powertrain characteristic correction factor Kf is obtained by testing:
the vehicle is enabled to run on a flat road at uniform acceleration, and the actual torque of the motor at m moments in the process of running at uniform acceleration is recordedActual rotating speed of motor at m moments Wherein, under the condition of uniform acceleration running of the vehicle, the motor model rotating speed at m moments Respectively equal to the actual rotational speed of the motor at the corresponding moment The pressure of the brake master cylinder at m moments is converted into the braking torque of the motor endAre all equal to zero;
fitting parameters A, B, C corresponding to the flat road and the actual rotating speed of the motor at m moments Respectively carrying into a driving resistance moment calculation formula: WR (pulse Width modulation)(t)=A×MotSpd(t) 2+B×MotSpd(t)+ C, calculating the running resistance at m times and converting the running resistance into the running resistance torque at the motor end
Respectively and correspondingly substituting the actual torque of the motor at m moments, the rotating speed of the motor model at m moments, the pressure of the brake master cylinder at m moments into the braking torque at the motor end and the driving resistance at m moments into the driving resistance torque at the motor end into formulas: MotSpdMod(t+dt)=MotSpdMod(t)+(M(t)-WR(t)-WB(t))×Kf1X dt, m Kf are calculated1;
For m Kf1Averaging, m Kf1Is used as the preset powertrain characteristic correction coefficient Kf.
6. The method of claim 5, wherein the table of speed-torque-torsional compensation coefficients is obtained by testing:
firstly, presetting a torsional vibration compensation coefficient as an empirical value or zero;
secondly, enabling the vehicle to run at the actual rotating speed and the actual torque of a certain motor; if the vehicle meets the driving smoothness evaluation standard, recording the actual rotating speed of the motor, the actual torque of the motor and the anti-torsional vibration compensation coefficient at the moment, if the vehicle does not meet the driving smoothness evaluation standard, changing the anti-torsional vibration compensation coefficient, enabling the vehicle to continuously run at the actual rotating speed of the motor and the actual torque of the motor until the vehicle meets the driving smoothness evaluation standard, and recording the actual rotating speed of the motor, the actual torque of the motor and the anti-torsional vibration compensation coefficient when the vehicle meets the driving smoothness evaluation standard;
thirdly, changing the actual rotating speed and/or the actual torque of the motor, enabling the vehicle to run at the changed actual rotating speed and/or the changed actual torque of the motor, recording the actual rotating speed, the actual torque and the torsional vibration resistance compensation coefficient of the motor at the moment if the vehicle meets the driving smoothness evaluation standard, changing the torsional vibration resistance compensation coefficient if the vehicle does not meet the driving smoothness evaluation standard, enabling the vehicle to continue running at the actual rotating speed and the actual torque of the motor until the vehicle meets the driving smoothness evaluation standard, and recording the actual rotating speed, the actual torque and the torsional vibration resistance compensation coefficient of the motor when the vehicle meets the driving smoothness evaluation standard;
and fourthly, combining all the actual rotating speeds of the motor, the actual torque of the motor and the torsional vibration resistance compensation coefficient which meet the driving smoothness evaluation standard to form the rotating speed-torque-torsional vibration resistance compensation coefficient table.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110604573.5A CN113511211B (en) | 2021-05-31 | 2021-05-31 | Torsional vibration control method based on electric driving system of electric vehicle |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110604573.5A CN113511211B (en) | 2021-05-31 | 2021-05-31 | Torsional vibration control method based on electric driving system of electric vehicle |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113511211A true CN113511211A (en) | 2021-10-19 |
CN113511211B CN113511211B (en) | 2022-09-06 |
Family
ID=78065325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110604573.5A Active CN113511211B (en) | 2021-05-31 | 2021-05-31 | Torsional vibration control method based on electric driving system of electric vehicle |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113511211B (en) |
Citations (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07245806A (en) * | 1994-03-01 | 1995-09-19 | Fuji Heavy Ind Ltd | Motor torque controller for electric car |
US20040084231A1 (en) * | 2001-11-24 | 2004-05-06 | Dieter Hoetzer | Method for controlling the operating response of a hybrid drive of a vehicle |
US20050080535A1 (en) * | 2003-10-14 | 2005-04-14 | Steinmetz Todd M. | Speed control for an electrically variable transmission |
US20110029172A1 (en) * | 2009-07-30 | 2011-02-03 | Hyundai Motor Company | Apparatus and method for controlling motor position and creep of electric vehicle |
JP2012090442A (en) * | 2010-10-20 | 2012-05-10 | Nissan Motor Co Ltd | Controller for electric vehicle |
US20130024057A1 (en) * | 2011-07-21 | 2013-01-24 | Hyundai Motor Company | Vehicle travel torque control system and control method thereof |
CN102933421A (en) * | 2010-06-11 | 2013-02-13 | 日立建机株式会社 | Pitching control device for electric vehicle |
US20130080012A1 (en) * | 2010-05-27 | 2013-03-28 | Nissan Motor Co., Ltd | Vibration-damping control device for vehicle |
US20130124024A1 (en) * | 2010-05-31 | 2013-05-16 | Nissan Motor Co., Ltd. | Torque response control apparatus for electric motor of vehicle |
US20140107877A1 (en) * | 2012-10-17 | 2014-04-17 | Kia Motors Corporation | Method and system for controlling anti-jerk of electric vehicle |
JP2014108732A (en) * | 2012-12-03 | 2014-06-12 | Nissan Motor Co Ltd | Braking/driving force control device for vehicle |
CN104755339A (en) * | 2012-08-16 | 2015-07-01 | 捷豹路虎有限公司 | System and method for controlling vehicle speed |
WO2016151661A1 (en) * | 2015-03-20 | 2016-09-29 | 日産自動車株式会社 | Damping control device for electric vehicle |
US20160375892A1 (en) * | 2015-06-26 | 2016-12-29 | Hyundai Motor Company | System and method for engine stop control of hybrid vehicle |
US20170197629A1 (en) * | 2014-05-30 | 2017-07-13 | Scania Cv Ab | Torque control of a vehicle powertrain based on a time derivative for a dynamic torque |
US20170225673A1 (en) * | 2015-08-26 | 2017-08-10 | Hyundai Motor Company | Apparatus and method for controlling vehicle having motor |
CN107839688A (en) * | 2017-10-19 | 2018-03-27 | 吉林大学 | A kind of electric automobile crawling method for controlling driving speed |
CN108422852A (en) * | 2018-03-21 | 2018-08-21 | 重庆长安汽车股份有限公司 | A kind of hybrid vehicle and its dynamic coupling system integrating device |
WO2018155625A1 (en) * | 2017-02-23 | 2018-08-30 | マツダ株式会社 | Drive control method and drive control device for hybrid vehicle |
JP2019093917A (en) * | 2017-11-22 | 2019-06-20 | トヨタ自動車株式会社 | Vehicular control device |
US20190381895A1 (en) * | 2017-01-24 | 2019-12-19 | Nissan Motor Co., Ltd. | Vehicle control device and control method |
CN111376907A (en) * | 2018-12-29 | 2020-07-07 | 北京宝沃汽车有限公司 | Engine torque compensation value obtaining method and device, controller and automobile |
US20200292075A1 (en) * | 2019-03-15 | 2020-09-17 | Hyundai Motor Company | Method for reducing torsional shock of driving system of electric vehicle |
CN112172811A (en) * | 2019-07-03 | 2021-01-05 | 丰田自动车株式会社 | Vehicle with a steering wheel |
CN112440971A (en) * | 2019-09-05 | 2021-03-05 | 长城汽车股份有限公司 | Torque control method and device and vehicle |
-
2021
- 2021-05-31 CN CN202110604573.5A patent/CN113511211B/en active Active
Patent Citations (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07245806A (en) * | 1994-03-01 | 1995-09-19 | Fuji Heavy Ind Ltd | Motor torque controller for electric car |
US20040084231A1 (en) * | 2001-11-24 | 2004-05-06 | Dieter Hoetzer | Method for controlling the operating response of a hybrid drive of a vehicle |
US20050080535A1 (en) * | 2003-10-14 | 2005-04-14 | Steinmetz Todd M. | Speed control for an electrically variable transmission |
US20110029172A1 (en) * | 2009-07-30 | 2011-02-03 | Hyundai Motor Company | Apparatus and method for controlling motor position and creep of electric vehicle |
US20130080012A1 (en) * | 2010-05-27 | 2013-03-28 | Nissan Motor Co., Ltd | Vibration-damping control device for vehicle |
US20130124024A1 (en) * | 2010-05-31 | 2013-05-16 | Nissan Motor Co., Ltd. | Torque response control apparatus for electric motor of vehicle |
CN102933421A (en) * | 2010-06-11 | 2013-02-13 | 日立建机株式会社 | Pitching control device for electric vehicle |
JP2012090442A (en) * | 2010-10-20 | 2012-05-10 | Nissan Motor Co Ltd | Controller for electric vehicle |
US20130024057A1 (en) * | 2011-07-21 | 2013-01-24 | Hyundai Motor Company | Vehicle travel torque control system and control method thereof |
CN104755339A (en) * | 2012-08-16 | 2015-07-01 | 捷豹路虎有限公司 | System and method for controlling vehicle speed |
US20140107877A1 (en) * | 2012-10-17 | 2014-04-17 | Kia Motors Corporation | Method and system for controlling anti-jerk of electric vehicle |
JP2014108732A (en) * | 2012-12-03 | 2014-06-12 | Nissan Motor Co Ltd | Braking/driving force control device for vehicle |
US20170197629A1 (en) * | 2014-05-30 | 2017-07-13 | Scania Cv Ab | Torque control of a vehicle powertrain based on a time derivative for a dynamic torque |
WO2016151661A1 (en) * | 2015-03-20 | 2016-09-29 | 日産自動車株式会社 | Damping control device for electric vehicle |
US20160375892A1 (en) * | 2015-06-26 | 2016-12-29 | Hyundai Motor Company | System and method for engine stop control of hybrid vehicle |
US20170225673A1 (en) * | 2015-08-26 | 2017-08-10 | Hyundai Motor Company | Apparatus and method for controlling vehicle having motor |
US20190381895A1 (en) * | 2017-01-24 | 2019-12-19 | Nissan Motor Co., Ltd. | Vehicle control device and control method |
WO2018155625A1 (en) * | 2017-02-23 | 2018-08-30 | マツダ株式会社 | Drive control method and drive control device for hybrid vehicle |
CN107839688A (en) * | 2017-10-19 | 2018-03-27 | 吉林大学 | A kind of electric automobile crawling method for controlling driving speed |
JP2019093917A (en) * | 2017-11-22 | 2019-06-20 | トヨタ自動車株式会社 | Vehicular control device |
CN108422852A (en) * | 2018-03-21 | 2018-08-21 | 重庆长安汽车股份有限公司 | A kind of hybrid vehicle and its dynamic coupling system integrating device |
CN111376907A (en) * | 2018-12-29 | 2020-07-07 | 北京宝沃汽车有限公司 | Engine torque compensation value obtaining method and device, controller and automobile |
US20200292075A1 (en) * | 2019-03-15 | 2020-09-17 | Hyundai Motor Company | Method for reducing torsional shock of driving system of electric vehicle |
CN112172811A (en) * | 2019-07-03 | 2021-01-05 | 丰田自动车株式会社 | Vehicle with a steering wheel |
CN112440971A (en) * | 2019-09-05 | 2021-03-05 | 长城汽车股份有限公司 | Torque control method and device and vehicle |
Non-Patent Citations (4)
Title |
---|
ANONYMOUS: "《IP.COM非专利全文库》", 5 May 2011, IP.COM * |
张农等: "双电机混合动力汽车起步与换挡过程协调控制", 《长安大学学报(自然科学版)》 * |
胡坚耀等: "基于模糊逻辑算法的纯电动客车起步加速能力优化与再生制动控制", 《集成技术》 * |
莫旭辉等: "电动汽车坡道起步电机转速控制研究", 《计算机仿真》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113511211B (en) | 2022-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6416441B1 (en) | process and a device for determining the driving torque of a vehicle | |
US8214096B2 (en) | Apparatus and method for controlling an accelerator for electric vehicles | |
KR960010678B1 (en) | Simulated engine characteristic control system | |
CN113267351B (en) | Automobile multidimensional working condition testing system and testing method | |
US6911791B2 (en) | Method and device for the reduction of load cycle oscillations in the drive train of a motor vehicle | |
EP1186462A3 (en) | Engine output and gear transmission ratio controller | |
EP1781484B1 (en) | Methods and apparatus for reducing vehicle tire wear | |
KR20160084426A (en) | Method and system for controlling the regenerative braking of an electric or hybrid motor vehicle | |
US8457823B2 (en) | Method and device for operating a hybrid vehicle | |
US6154702A (en) | Method and apparatus for estimating applied wheel torque in a motor vehicle | |
CN107575570A (en) | A kind of double-clutch automatic gearbox is creeped vibration control method | |
US20220234498A1 (en) | Sound signal generation device, sound signal generation method, and sound signal generation program | |
CN108844745B (en) | Kinetic parameter identification method of hybrid electric vehicle transmission system | |
CN101890956A (en) | Measure the apparatus and method of the quality of motor vehicle driven by mixed power | |
JP4665790B2 (en) | Vehicle vibration reduction control device | |
CN113650615A (en) | Gear shifting control method and device and storage medium | |
CN113511211B (en) | Torsional vibration control method based on electric driving system of electric vehicle | |
CN113200034A (en) | Torque control method and device, storage medium and vehicle | |
CN107428240B (en) | Method for operating a drive train of a motor vehicle and corresponding drive train | |
US9533601B2 (en) | Apparatus for controlling rotary machine | |
CN113386578B (en) | Torque output control method and device, vehicle and equipment | |
CN112503170A (en) | Starting gear decision method, device and system for mechanical automatic transmission | |
CN105046015B (en) | The load transfer method of analogy material hesitation | |
CN1940270A (en) | Method and device for running internal-combustion engine | |
CN113464641A (en) | Automatic gearbox clutch torque control system and control method thereof |
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 |