WO2011036821A1 - 車両の旋回特性推定装置 - Google Patents

車両の旋回特性推定装置 Download PDF

Info

Publication number
WO2011036821A1
WO2011036821A1 PCT/JP2009/067128 JP2009067128W WO2011036821A1 WO 2011036821 A1 WO2011036821 A1 WO 2011036821A1 JP 2009067128 W JP2009067128 W JP 2009067128W WO 2011036821 A1 WO2011036821 A1 WO 2011036821A1
Authority
WO
WIPO (PCT)
Prior art keywords
vehicle
yaw rate
time constant
value
response time
Prior art date
Application number
PCT/JP2009/067128
Other languages
English (en)
French (fr)
Inventor
尚大 横田
Original Assignee
トヨタ自動車株式会社
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 トヨタ自動車株式会社 filed Critical トヨタ自動車株式会社
Priority to CN2009801616535A priority Critical patent/CN102574527B/zh
Priority to PCT/JP2009/067128 priority patent/WO2011036821A1/ja
Priority to DE112009005281.9T priority patent/DE112009005281B4/de
Priority to US13/389,696 priority patent/US8744686B2/en
Priority to JP2010541617A priority patent/JP5041069B2/ja
Publication of WO2011036821A1 publication Critical patent/WO2011036821A1/ja

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/11Pitch movement
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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/00Purposes 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/18Propelling the vehicle
    • B60W30/18009Propelling the vehicle related to particular drive situations
    • B60W30/18145Cornering
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/112Roll movement
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/114Yaw movement
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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
    • B60W50/00Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
    • B60W2050/0062Adapting control system settings
    • B60W2050/0075Automatic parameter input, automatic initialising or calibrating means
    • B60W2050/0083Setting, resetting, calibration
    • B60W2050/0088Adaptive recalibration
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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
    • B60W2520/00Input parameters relating to overall vehicle dynamics
    • B60W2520/14Yaw

Definitions

  • the present invention relates to a turning characteristic estimation device for a vehicle, and more particularly to turning a vehicle for estimating a steering response time constant coefficient indicating a turning characteristic of the vehicle based on a relationship between a reference yaw rate and a transient yaw rate when the vehicle turns.
  • the present invention relates to a characteristic estimation apparatus.
  • the actual yaw rate of the vehicle is in a first order lag relationship with the standard yaw rate of the vehicle, and if the coefficient for the vehicle speed of the time constant of the first order lag is the steering response time constant coefficient, the vehicle stability factor and the steering response time constant coefficient are Represents the turning characteristics of Estimating the vehicle stability factor and the steering response time constant coefficient by estimating the parameters a and b of the discrete time transfer function from the reference yaw rate to the actual yaw rate using an ARX (auto-regressive exogenous model) it can.
  • ARX auto-regressive exogenous model
  • a standard yaw rate is estimated based on travel data when a vehicle turns, and parameters a and b of a discrete time transfer function from the standard yaw rate to the actual yaw rate are estimated.
  • a vehicle turning characteristic estimation device that estimates a steering response time constant coefficient based on parameters and estimates a vehicle stability factor based on parameters a and b.
  • a large amount of data on the standard yaw rate and the actual yaw rate is necessary to estimate the parameters a and b. Therefore, there is a problem that it takes a long time to estimate the vehicle stability factor and the steering response time constant coefficient.
  • the main object of the present invention is to steer on the basis of the transient yaw rate of the vehicle and the actual yaw rate of the vehicle, which are in a first-order lag relationship with the steady reference yaw rate of the vehicle, in view of the above-described problems in the conventional turning characteristic estimation device.
  • Estimate the steering response time constant coefficient in less time than before by correcting the estimated response time constant coefficient by learning so that the estimated response time constant coefficient approaches the true steering response time constant coefficient It is to be.
  • the transient yaw rate of the vehicle that is in a first order lag relationship with respect to the steady reference yaw rate of the vehicle is regarded as a value corresponding to the actual yaw rate of the vehicle, and at the time of steering response that is a coefficient relating to the vehicle speed of the time constant of first order delay
  • a vehicle turning characteristic estimation device that corrects an estimated value of a steering response time constant coefficient based on a relationship between a transient yaw rate of a vehicle and an actual yaw rate of the vehicle so as to approach the actual yaw rate.
  • the estimated value of the steering response time constant coefficient since the estimated value of the steering response time constant coefficient is corrected so that the transient yaw rate of the vehicle approaches the actual yaw rate of the vehicle, the estimated value of the steering response time constant coefficient approaches the true steering response time constant coefficient. It can be modified as follows. Therefore, it is not necessary to estimate the parameter of the discrete time transfer function from the reference yaw rate to the actual yaw rate, and the steering response time constant coefficient can be estimated in a shorter time than in the past.
  • At least one of the magnitude of the deviation between the transient yaw rate of the vehicle and the actual yaw rate of the vehicle and the phase difference between the transient yaw rate of the vehicle and the actual yaw rate of the vehicle is reduced, so that the transient yaw rate of the vehicle is reduced.
  • the estimated value of the steering response time constant coefficient may be corrected so as to approach the actual yaw rate.
  • at least one of the magnitude of the deviation between the transient yaw rate of the vehicle and the actual yaw rate of the vehicle and the phase difference between the transient yaw rate of the vehicle and the actual yaw rate of the vehicle are reduced, whereby the transient yaw rate of the vehicle is reduced.
  • the estimated value of the steering response time constant coefficient can be corrected so as to approach the actual yaw rate.
  • the stored steering response time constant coefficient is used as the first reference value, and the first reference value, the second reference value larger than the first reference value, and the first reference value
  • the transient yaw rate of the three vehicles is calculated using the third reference value which is also smaller, and the index value of the magnitude of the deviation between the vehicle's transient yaw rate and the vehicle's actual yaw rate among the first to third reference values
  • the reference value with the smallest yaw rate deviation index value may be used as the estimated value of the steering response time constant coefficient.
  • the reference value having the smallest yaw rate deviation index value among the first to third reference values is set as the estimated value of the steering response time constant coefficient. Therefore, the value closest to the true steering response time constant coefficient among the first to third reference values can be used as the estimated value of the steering response time constant coefficient.
  • the estimated value of the steering response time constant coefficient is set as a new first reference value, a new first reference value, and a new second reference value that is larger than the new first reference value.
  • the yaw rate deviation index value is The smallest reference value may be used as a new estimated value of the steering response time constant coefficient.
  • the value closest to the true steering response time constant coefficient among the new first to third reference values is set as the estimated value of the steering response time constant coefficient. Accordingly, an estimated value of the steering response time constant coefficient closer to the true steering response time constant coefficient than the estimated value of the steering response time constant coefficient selected from the first to third reference values can be obtained. Therefore, by repeating the procedure of this configuration, the estimated value of the steering response time constant coefficient can be gradually brought closer to the true steering response time constant coefficient.
  • the difference between the new second reference value and the new first reference value is smaller than the difference between the second reference value and the first reference value.
  • the difference between the value and the new first reference value may be smaller than the difference between the third reference value and the first reference value.
  • the difference between the new second reference value and the new first reference value is greater than or equal to the difference between the second reference value and the first reference value, and the new third reference value.
  • the estimated value of the steering response time constant coefficient is increased faster during the true steering response. It can be close to a constant coefficient.
  • the yaw rate deviation index value is a deviation between the actual yaw rate of the vehicle from which the component below the first predetermined frequency is removed and the transient yaw rate of the vehicle from which the component below the second predetermined frequency is removed. It may be a value based on the size of.
  • the transient yaw rate is obtained based on the state quantity of the vehicle, and the state quantity of the vehicle and the actual yaw rate of the vehicle are detected by the detecting means. According to the above configuration, it is possible to eliminate stationary detection errors such as the zero point offset of the detection means for detecting the vehicle state quantity and the actual yaw rate of the vehicle, so that the estimation accuracy of the steering response time constant coefficient is increased. Can do.
  • the first predetermined frequency or the second predetermined frequency may be changed according to an index value of the amount of reciprocating steering performed by the driver per unit time.
  • a stationary detection error such as a zero point offset of the detection means for detecting the state quantity of the vehicle such as the actual yaw rate varies depending on the amount of reciprocating steering by the driver per unit time. According to the above configuration, it is possible to properly eliminate stationary detection errors according to the amount of reciprocating steering performed by the driver per unit time.
  • the first predetermined frequency or the second predetermined frequency may be changed according to the magnitude of the longitudinal acceleration of the vehicle.
  • the influence of the steady detection error such as the zero point offset of the detection means for detecting the state quantity of the vehicle such as the actual yaw rate on the estimation accuracy of the steering response time constant coefficient is the magnitude of the change in vehicle speed, that is, the longitudinal acceleration of the vehicle. It also varies depending on the size. According to the above configuration, it is possible to properly eliminate stationary detection errors in accordance with the magnitude of the longitudinal acceleration of the vehicle.
  • the yaw rate deviation index value may be a value obtained by converting the magnitude of the deviation between the transient yaw rate of the vehicle and the actual yaw rate of the vehicle into the magnitude of the deviation of the steering angle of the front wheels.
  • the magnitude of the deviation between the vehicle's transient yaw rate and the vehicle's actual yaw rate depends on the vehicle speed
  • the magnitude of the deviation between the vehicle's transient yaw rate and the vehicle's actual yaw rate is converted to the magnitude of the deviation of the steering angle of the front wheels.
  • the value does not depend on the vehicle speed. Therefore, according to the above configuration, the estimated value of the steering response time constant coefficient can be obtained based on the yaw rate deviation index value that does not depend on the vehicle speed, thereby estimating the steering response time constant coefficient without being influenced by the vehicle speed. Can do.
  • the deviation of the steering angle of the front wheels is a difference between the steering angle of the front wheels and the actual steering angle of the front wheels for achieving the transient yaw rate of the vehicle.
  • the estimated value of the steering response time constant coefficient is increased and corrected so that the phase of the transient yaw rate of the vehicle is the phase of the actual yaw rate of the vehicle.
  • the estimated value of the steering response time constant coefficient may be reduced and corrected.
  • the steering response time constant coefficient is estimated so that the estimated value of the steering response time constant coefficient approaches the true steering response time constant coefficient by reducing the phase difference between the transient yaw rate of the vehicle and the actual yaw rate of the vehicle.
  • the value can be modified.
  • the steering response time constant coefficient may be estimated individually for the right turn and the left turn of the vehicle. According to this configuration, even when the turning characteristic of the vehicle differs depending on the turning direction of the vehicle, for example, when the center of gravity of the vehicle is not in the center of the vehicle in the lateral direction or when the vehicle turns in the lateral direction, the right side of the vehicle The steering response time constant coefficient can be estimated for each of the turn and the left turn.
  • the steering response time constant coefficient may be estimated for each region of the magnitude of the lateral acceleration of the vehicle.
  • the magnitude of the deviation between the transient yaw rate of the vehicle and the actual yaw rate of the vehicle also varies depending on the magnitude of the lateral acceleration of the vehicle.
  • the steering response time constant coefficient is estimated for each region of the magnitude of the lateral acceleration of the vehicle, thereby eliminating the influence of the magnitude of the lateral acceleration of the vehicle and estimating the steering response time constant coefficient. Can do. In the two-wheel model of the vehicle shown in FIG.
  • the vehicle mass and yaw moment of inertia are M and I, respectively, and the distances between the center of gravity 102 of the vehicle and the front and rear axles are Lf and Lr, respectively.
  • the cornering forces of the front wheel 100f and the rear wheel 100r are respectively Ff and Fr, and the cornering powers of the front wheel and the rear wheel are respectively Kf and Kr.
  • the actual steering angle of the front wheel 100f is ⁇
  • the slip angles of the front and rear wheels are ⁇ f and ⁇ r, respectively, and the slip angle of the vehicle body is ⁇ .
  • the Laplace operator is converted to Laplace using the Laplace operator as s, and the following equations 8 to 10 are established by arranging the yaw rate ⁇ .
  • a reference yaw rate ⁇ (s) is determined.
  • Kh in the above equation 9 is a stability factor
  • Tp in the above equation 10 is a coefficient relating to the vehicle speed V of the first-order lag system having a time constant depending on the vehicle speed, that is, the “steering response time constant coefficient” in this specification. It is a coefficient to call.
  • the above equation 8 is an equation for calculating the yaw rate ⁇ of the vehicle from the actual steering angle ⁇ of the front wheels, the vehicle speed V, and the lateral acceleration Gy.
  • the yaw rate calculated from this linearization model is assumed to be a transient yaw rate ⁇ tr having a first-order lag with respect to the steady-state standard yaw rate ⁇ t expressed by the following equation 11.
  • the transient yaw rate ⁇ tr is calculated according to the above equation 8, and the estimated value of the steering response time constant coefficient Tp is based on the relationship between the transient yaw rate ⁇ tr and the actual yaw rate ⁇ so that the transient yaw rate ⁇ tr approaches the actual yaw rate ⁇ . May be modified.
  • the yaw rate deviation index value may be an integrated value of an index value of the magnitude of deviation between the vehicle transient yaw rate and the vehicle actual yaw rate.
  • the component below the first predetermined frequency may be removed from the actual yaw rate of the vehicle by the high-pass filter processing
  • the component below the second predetermined frequency may be removed from the transient yaw rate of the vehicle by the high-pass filter processing.
  • the first and second predetermined frequencies may be the same frequency.
  • the vehicle speed is V
  • the vehicle wheelbase is L
  • the magnitude of the deviation between the vehicle's transient yaw rate and the vehicle's actual yaw rate is multiplied by L / V, so that the vehicle's transient yaw rate is obtained.
  • a value obtained by converting the magnitude of the deviation between the actual yaw rate of the vehicle and the magnitude of the deviation of the steering angle of the front wheels may be calculated.
  • the estimated value of the steering response time constant coefficient may be corrected with the increase / decrease correction amount, and the magnitude of the increase / decrease correction amount may be gradually reduced.
  • the steering response time constant coefficient is estimated when the deviation between the estimated value of the current steering response time constant coefficient and the estimated value of the previous steering response time constant coefficient is less than the reference value. The correction of the value may be terminated. Further, in the above configuration, the estimated value of the steering response time constant coefficient is corrected with the increase / decrease correction amount so as not to exceed the range between the preset minimum value and maximum value of the steering response time constant coefficient.
  • the steering response time constant coefficient is not estimated. It may be.
  • FIG. 1 is a schematic configuration diagram showing a first embodiment of a turning characteristic estimation device according to the present invention applied to a vehicle motion control device.
  • FIG. 2 is a flowchart showing an estimation calculation routine of the steering response time constant coefficient Tp in the first embodiment.
  • FIG. 3 is a flow chart showing a main part of an estimation calculation routine of the steering response time constant coefficient Tp in the second embodiment of the turning characteristic estimation device according to the present invention configured as a modification of the first embodiment.
  • FIG. 4 is a graph showing the relationship between the steering frequency fs and the cutoff frequency fhc of the high-pass filter process.
  • FIG. 5 is a flowchart showing a main part of an estimation calculation routine of the steering response time constant coefficient Tp in the third embodiment of the turning characteristic estimation device according to the present invention configured as a modification of the first embodiment.
  • FIG. 6 is a graph showing the relationship between the steering frequency fs, the cutoff frequency fhc of the high-pass filter process, and the absolute value of the longitudinal acceleration Gx of the vehicle.
  • FIG. 7 is a flowchart showing a main part of a steering response time constant coefficient Tp estimation calculation routine in the fourth embodiment of the turning characteristic estimation device according to the present invention configured as a modification of the first embodiment.
  • FIG. 6 is a graph showing the relationship between the steering frequency fs, the cutoff frequency fhc of the high-pass filter process, and the absolute value of the longitudinal acceleration Gx of the vehicle.
  • FIG. 7 is a flowchart showing a main part of a steering response time constant coefficient Tp estimation calculation routine in the fourth embodiment of the turning characteristic estimation device according to the present invention configured as
  • FIG. 8 is a flowchart showing the first half of a steering response time constant coefficient Tp estimation calculation routine in the fifth embodiment of the turning characteristic estimation apparatus according to the present invention, which is configured as a modification of the first embodiment.
  • FIG. 9 is a flowchart showing the latter half of the steering response time constant coefficient Tp estimation calculation routine in the fifth embodiment of the turning characteristic estimation apparatus according to the present invention, which is configured as a modification of the first embodiment.
  • FIG. 10 is a flowchart showing the first half of the steering response time constant coefficient Tp estimation calculation routine in the sixth embodiment of the turning characteristic estimation apparatus according to the present invention, which is configured as a modification of the first embodiment.
  • FIG. 9 is a flowchart showing the latter half of the steering response time constant coefficient Tp estimation calculation routine in the fifth embodiment of the turning characteristic estimation apparatus according to the present invention, which is configured as a modification of the first embodiment.
  • FIG. 10 is a flowchart showing the first half of the steering response time constant coefficient Tp estimation calculation routine in the sixth embodiment of the turning characteristic
  • FIG. 11 is a flowchart showing the latter half of the routine for estimating the steering response time constant coefficient Tp in the sixth embodiment of the turning characteristic estimation device according to the present invention, which is configured as a modification of the first embodiment.
  • FIG. 12 is a flowchart showing an estimation calculation routine of the steering response time constant coefficient Tp in the seventh embodiment of the turning characteristic estimation apparatus according to the present invention.
  • FIG. 13 is a flowchart showing a subroutine for estimating and storing the steering response time constant coefficient Tp in step 200 of FIG.
  • FIG. 14 is a flowchart showing an estimation calculation routine of the steering response time constant coefficient Tp in the eighth embodiment of the turning characteristic estimation apparatus according to the present invention.
  • FIG. 12 is a flowchart showing an estimation calculation routine of the steering response time constant coefficient Tp in the seventh embodiment of the turning characteristic estimation apparatus according to the present invention.
  • FIG. 13 is a flowchart showing a subroutine for estimating and storing the steering response time constant coefficient Tp in step 200 of FIG.
  • FIG. 15 is a flowchart showing a subroutine for estimating and storing the steering response time constant coefficient Tp in step 200 of FIG.
  • FIG. 16 is a flowchart showing a subroutine for estimating and storing the steering response time constant coefficient Tp in step 205 of FIG.
  • FIG. 17 is a flowchart showing an estimation calculation routine of the steering response time constant coefficient Tp in the ninth embodiment of the turning characteristic estimation apparatus according to the present invention.
  • FIG. 18 is a flowchart showing a subroutine for estimating and storing the steering response time constant coefficient Tp in step 200 of FIG.
  • FIG. 19 is a flowchart showing a subroutine for estimating and storing the steering response time constant coefficient Tp in step 205 of FIG.
  • FIG. 20 is an explanatory diagram showing a two-wheel model of a vehicle for estimating a steering response time constant coefficient.
  • FIG. 1 is a schematic configuration diagram showing a first embodiment of a turning characteristic estimation device according to the present invention applied to a vehicle motion control device.
  • reference numeral 50 denotes a motion control device for the vehicle 10 as a whole, and the turning characteristic estimation device according to the present invention is a part of the motion control device 50.
  • the vehicle 10 has left and right front wheels 12FL and 12FR and left and right rear wheels 12RL and 12RR.
  • the left and right front wheels 12FL and 12FR which are steered wheels, are steered via tie rods 18L and 18R by a rack and pinion type power steering device 16 that is driven in response to turning of the steering wheel 14 by the driver.
  • the braking force of each wheel is controlled by controlling the braking pressure of the wheel cylinders 24FR, 24FL, 24RR, and 24RL by the hydraulic circuit 22 of the braking device 20.
  • the hydraulic circuit 22 includes an oil reservoir, an oil pump, various valve devices, and the like, and the braking pressure of each wheel cylinder is normally driven according to the depression operation of the brake pedal 26 by the driver. It is controlled by the master cylinder 28 and, if necessary, is controlled by the electronic control unit 30 as described later.
  • the vehicle 10 includes a yaw rate sensor 36 for detecting the actual yaw rate ⁇ of the vehicle, a longitudinal acceleration sensor 38 for detecting the longitudinal acceleration Gx of the vehicle, a lateral acceleration sensor 40 for detecting the lateral acceleration Gy of the vehicle, and a wheel speed Vwi of each wheel.
  • the steering angle sensor 34, the yaw rate sensor 36, and the lateral acceleration sensor 40 detect the steering angle, the actual yaw rate, and the lateral acceleration, respectively, with the left turning direction of the vehicle being positive.
  • a signal indicating the longitudinal acceleration Gx detected by the sensor 38, a signal indicating the lateral acceleration Gy detected by the lateral acceleration sensor 40, and a signal indicating the wheel speed Vwi detected by the wheel speed sensors 42FR to 42RL are sent to the electronic control unit 30. Entered.
  • the electronic control unit 30 includes, for example, a CPU, a ROM, an EEPROM, a RAM, a buffer memory, and an input / output port device, which are connected to each other by a bidirectional common bus.
  • a microcomputer with a general configuration.
  • the ROM stores the default values of the stability factor Kh and the steering response time constant coefficient Tp used for calculating the reference yaw rate ⁇ t. These default values are set for each vehicle when the vehicle is shipped.
  • the EEPROM stores the estimated value of the steering response time constant coefficient Tp and the like, and the estimated value of the steering response time constant coefficient Tp is calculated based on the running data of the vehicle when the vehicle is in a turning state, as will be described in detail later.
  • the estimated value of the steering response time constant coefficient Tp is calculated and stored in the buffer memory. Further, when the electronic control unit 30 calculates the estimated value of the steering response time constant coefficient Tp, based on the estimated value, the standard value Tp0 and the reference value interval ⁇ Tp for calculating the five reference values of the steering response time constant coefficient Tp. Update. Further, the electronic control unit 30 calculates the target yaw rate ⁇ tt corresponding to the transient yaw rate ⁇ tr using the estimated value of the steering response time constant coefficient Tp stored in the EEPROM, and the deviation between the yaw rate detection value ⁇ and the target yaw rate ⁇ tt. As a result, the yaw rate deviation ⁇ is calculated.
  • the electronic control unit 30 determines whether or not the turning behavior of the vehicle has deteriorated by determining whether or not the magnitude of the yaw rate deviation ⁇ exceeds the reference value ⁇ o (positive value).
  • the vehicle motion control performed by the electronic control unit 30 is an arbitrary control as long as the vehicle motion is controlled based on the target yaw rate ⁇ tt calculated using the estimated value of the steering response time constant coefficient Tp. It's okay.
  • the routine for estimating the steering response time constant coefficient Tp in the first embodiment will be described with reference to the flowchart shown in FIG. The control according to the flowchart shown in FIG.
  • step 10 control is started from step 10, and in step 20, a signal indicating the steering angle ⁇ detected by each sensor is read.
  • step 30 low-pass filter processing for removing high-frequency noise is performed on the steering angle ⁇ read in step 20 and the like.
  • the low-pass filter process in this case may be a primary low-pass filter process with a cutoff frequency of 3.4 Hz, for example.
  • step 40 the stability factor Kh is set to a default value set in advance when the vehicle is shipped. When the stability factor is estimated based on the travel data of the vehicle, the stability factor Kh may be set to the estimated value.
  • step 50 the vehicle speed V is calculated based on the wheel speed Vwi
  • the steering angle ⁇ of the front wheel is calculated based on the steering angle ⁇
  • the standard yaw rate ⁇ t is calculated according to the above equation 11.
  • step 80 high-pass filter processing for removing the influence of the zero offset of the sensor is performed on the actual yaw rate ⁇ subjected to the low-pass filter processing in step 30 and the transient yaw rate ⁇ trn calculated in step 70. .
  • the high-pass filter process in this case may be a primary high-pass filter process with a cutoff frequency of 0.2 Hz, for example.
  • step 90 it is determined whether or not the vehicle is in a turning state. If a negative determination is made, the control returns to step 20, and if an affirmative determination is made, the control proceeds to step 100.
  • whether or not the vehicle is turning is determined whether or not the absolute value of the lateral acceleration Gy of the vehicle is greater than or equal to the reference value in a situation where the vehicle is traveling at a vehicle speed greater than or equal to the reference value. The determination is made by determining whether the absolute value of the actual yaw rate ⁇ of the vehicle is greater than or equal to a reference value and whether the absolute value of the product of the yaw rate ⁇ of the vehicle and the vehicle speed V is greater than or equal to the reference value. Good.
  • the integrated value ⁇ an of the yaw rate deviation index value is cleared to 0 at step 110, and when a negative determination is made, control proceeds to step 120.
  • any of the following (A1) and (A2) it may be determined that the integrated value ⁇ an of the yaw rate deviation index value needs to be adjusted.
  • step 150 of the previous cycle the standard value Tp0 of the steering response time constant coefficient Tp is updated.
  • step 160 of the previous cycle the reference value interval ⁇ Tp is updated.
  • the stability factor Kh is estimated in step 40, the following (A3) is added to the determination condition, and when any of (A1) to (A3) is satisfied, the yaw rate deviation index value It may be determined that it is necessary to adjust the integrated value ⁇ an.
  • (A3) The absolute value of the deviation ⁇ Kh between the stability factor Kh when the integrated value ⁇ an of the yaw rate deviation index value was adjusted last time and the stability factor Kh estimated in step 40 of the current cycle is the stability factor. It is above the standard value for deviation.
  • the integrated value ⁇ an of the yaw rate deviation index value may be calculated as an integrated value of the magnitude of the deviation between the actual yaw rate ⁇ bpf after the bandpass filter processing and the transient yaw rate ⁇ trbpfn according to the following equation 19.
  • step 130 it is determined whether or not the estimation of the steering response time constant coefficient Tp is permitted. If a negative determination is made, the control returns to step 20 and an affirmative determination is made. Sometimes control proceeds to step 140. In this case, when the following (B1) or (B2) is established, it may be determined that the estimation of the steering response time constant coefficient Tp is permitted. (B1) The reference time or more has elapsed since the integrated value ⁇ an of the yaw rate deviation index value was cleared to 0 in the previous step 110.
  • step 140 The minimum value of the integrated values ⁇ a1 to ⁇ a5 of the five yaw rate deviation index values is equal to or greater than the reference value of the integrated values.
  • step 140 the smallest value ⁇ am (m is any one of 1 to 5) among the integrated values ⁇ a1 to ⁇ a5 of the five yaw rate deviation index values is determined.
  • the steering response time constant coefficient Tpm corresponding to the minimum value ⁇ am of the integrated values of the yaw rate deviation index values is set as the estimated value of the steering response time constant coefficient Tp.
  • Ka is a coefficient larger than 0 and smaller than 1
  • the guard minimum value of the reference value interval ⁇ Tp is ⁇ Tpmin (a positive constant)
  • the reference value interval ⁇ Tp is Ka ⁇ the current reference value interval ⁇ Tp and It is updated to the larger one of the guard minimum values ⁇ Tpmin and stored in the buffer memory.
  • step 170 it is determined whether or not the storage of the estimated value of the steering response time constant coefficient Tp in the EEPROM is permitted.
  • step 180 the estimated value of the steering response time constant coefficient Tp is stored in the EEPROM, whereby the estimated value of the steering response time constant coefficient Tp stored in the EEPROM is updated.
  • the steering response time constant to the EEPROM It may be determined that the storage of the estimated value of the coefficient Tp is permitted.
  • the steady-state standard yaw rate ⁇ t is calculated in step 50, and five reference values Tpn of different sizes are set in step 60.
  • the transient yaw rate ⁇ trn is calculated for the five reference values Tpn based on the steady standard yaw rate ⁇ t.
  • the integrated value ⁇ a1 ⁇ of the five yaw rate deviation index values obtained by replacing the magnitude of the deviation between the actual yaw rate ⁇ bpf and the transient yaw rate ⁇ trbpfn after the bandpass filter processing with the magnitude of the steering angle deviation of the front wheels. ⁇ a5 is calculated.
  • step 140 the smallest value ⁇ am of the integrated values ⁇ a1 to ⁇ a5 of the five yaw rate deviation index values is determined, and the steering response time constant coefficient Tpm corresponding to the minimum value ⁇ am of the integrated values of the yaw rate deviation index values is steered. This is an estimated value of the response time constant coefficient Tp.
  • the transient yaw rate ⁇ trn is calculated for the five reference values Tpn, and the reference value that minimizes the deviation between the transient yaw rate ⁇ trn and the actual yaw rate ⁇ is the steering response time constant coefficient Tp. Estimated value.
  • the steering response time constant coefficient corresponding to the transient yaw rate ⁇ trn closest to the actual yaw rate ⁇ that is, the reference value closest to the actual steering response time constant coefficient can be used as the estimated value of the steering response time constant coefficient Tp.
  • the reference value interval ⁇ Tp is updated so that the size thereof becomes smaller.
  • five reference values Tpn of different sizes are set based on the updated standard value Tp0 of the steering response time constant coefficient Tp and the reference value interval ⁇ Tp.
  • the magnitude of the deviation between the transient yaw rate ⁇ trn and the actual yaw rate ⁇ can be gradually reduced to make the transient yaw rate ⁇ trn gradually approach the actual yaw rate ⁇ . Therefore, the steering response time constant coefficient corresponding to the transient yaw rate ⁇ trn closest to the actual yaw rate ⁇ , that is, the estimated value of the steering response time constant coefficient Tp can be gradually brought closer to the actual steering response time constant coefficient.
  • the steady-state normative yaw rate ⁇ t is calculated based on the steering angle ⁇ or the like subjected to the low-pass filter processing in Step 30.
  • step 80 the high-pass filter process is performed on the transient yaw rate ⁇ trn and the actual yaw rate ⁇ calculated based on the steady reference yaw rate ⁇ t in step 70. Therefore, not only the high-frequency noise included in the detected steering angle ⁇ etc. can be removed, but also the influence of the zero point offset of the yaw rate sensor 36 etc. can be eliminated. Therefore, since the deviation between the actual yaw rate and the transient yaw rate can be calculated by eliminating the influence of the sensor zero offset, the steering response compared to the case where the high-pass filter processing is not performed on the actual yaw rate ⁇ and the transient yaw rate ⁇ trn. The time constant coefficient Tp can be accurately estimated.
  • the number of high-pass filter processes can be reduced as compared with the case where the high-pass filter process is performed on the steering angle ⁇ and the lateral acceleration Gy used for the calculation of the steady-state standard yaw rate ⁇ t.
  • the calculation load can be reduced.
  • the bandpass filter processing may be performed on the actual yaw rate ⁇ and the transient yaw rate ⁇ trn without performing the lowpass filter processing on the steering angle ⁇ or the like. In that case, it is possible to reduce the number of computations required for the filter processing as compared with the case of the first embodiment, while effectively removing high-frequency noise, and thereby the computation load of the electronic control unit 30. Can be reduced.
  • the steering response time constant coefficient Tp is estimated based on the integrated value ⁇ an of the yaw rate deviation index value based on the magnitude of the deviation between the actual yaw rate ⁇ bpf and the transient yaw rate ⁇ trbpfn after the bandpass filter processing. A value is determined. Therefore, the actual yaw rate ⁇ and the transient yaw rate ⁇ tr are instantaneous as compared with the case where the estimated value of the steering response time constant coefficient Tp is obtained based on the magnitude of the deviation between the actual yaw rate ⁇ bpf and the transient yaw rate ⁇ trbpfn after the bandpass filter processing.
  • the integrated value ⁇ an of the yaw rate deviation index value is an integrated value of a value obtained by replacing the magnitude of the deviation between the actual yaw rate ⁇ bpf and the transient yaw rate ⁇ trbpfn with the magnitude of the deviation of the steering angle of the front wheels. It is. Therefore, the steering response time constant coefficient Tp can be estimated without being affected by the vehicle speed V.
  • the steering response time constant coefficient Tp can be accurately estimated as compared with the case where the integrated value ⁇ an of the yaw rate deviation index value is, for example, an integrated value of the magnitude of the deviation between the actual yaw rate ⁇ bpf and the transient yaw rate ⁇ trbpfn. Further, the trouble of estimating the steering response time constant coefficient Tp for each vehicle speed V or changing the steering response time constant coefficient Tp used for calculating the target yaw rate ⁇ tt according to the vehicle speed V is avoided. The capacity of the storage means can be reduced.
  • step 100 it is determined in step 100 whether or not the integrated value ⁇ an of the yaw rate deviation index value needs to be adjusted, and when an affirmative determination is made, in step 110.
  • the integrated value ⁇ an of the yaw rate deviation index value is cleared to zero. Therefore, in a situation where the standard value Tp0 of the steering response time constant coefficient Tp and the reference value interval ⁇ Tp are changed, the integrated value ⁇ an of the previous yaw rate deviation index value has an adverse effect on the estimation of the steering response time constant coefficient Tp. This can be reliably prevented.
  • step 130 it is determined in step 130 whether or not the estimation of the steering response time constant coefficient Tp is permitted. Is executed.
  • step 170 it is determined whether or not the storage of the estimated value of the steering response time constant coefficient Tp is permitted, and when an affirmative determination is made.
  • step 180 the estimated value of the steering response time constant coefficient Tp is stored in the EEPROM. Therefore, the estimated value of the steering response time constant coefficient Tp can be stored in the EEPROM when the estimated value of the steering response time constant coefficient Tp substantially matches the actual steering response time constant coefficient.
  • step 90 it is determined in step 90 whether or not the vehicle is in a turning traveling state, and when an affirmative determination is made, step 100 and subsequent steps are executed.
  • FIG. 3 is a flowchart showing a main part of a routine for estimating the steering response time constant coefficient Tp in the second embodiment of the turning characteristic estimation device according to the present invention, which is configured as a modification of the first embodiment.
  • steps corresponding to the steps shown in FIG. 2 are assigned the same step numbers as those shown in FIG. The same applies to the flowchart of FIG.
  • step 72 when step 70 is completed, in step 72, the number of reciprocating steering operations by the driver per unit time is calculated as the steering frequency fs. Further, the cut-off frequency fhc is calculated from the map corresponding to the graph shown in FIG. 4 based on the steering frequency fs so that the cut-off frequency fhc of the high-pass filter processing in step 80 becomes smaller as the steering frequency fs becomes lower. . In the high-pass filter processing of the actual yaw rate ⁇ and the transient yaw rate ⁇ trn in step 80, the cutoff frequency is set to the cutoff frequency fhc calculated in step 72. In the first embodiment described above, the cut-off frequency fhc of the high-pass filter process in step 80 is constant.
  • the cutoff frequency fhc is variably set according to the steering frequency fs so that the cutoff frequency fhc becomes smaller as the steering frequency fs becomes lower.
  • FIG. 5 is a flowchart showing a main part of a steering response time constant coefficient Tp estimation calculation routine in the third embodiment of the turning characteristic estimation device according to the present invention, which is configured as a modification of the first embodiment.
  • step 74 the number of reciprocating steering operations by the driver per unit time is calculated as the steering frequency fs. Further, the cutoff frequency fhc of the high-pass filter processing decreases as the steering frequency fs decreases, and the cutoff frequency fhc of the high-pass filter processing increases as the absolute value of the longitudinal acceleration Gx of the vehicle increases. Based on the absolute value of the longitudinal acceleration Gx, the cutoff frequency fhc is calculated from a map corresponding to the graph shown in FIG.
  • the cut-off frequency is set to the cut-off frequency fhc calculated in step 74. Also, assuming that the errors of the zero offset of the sensor for the vehicle yaw rate ⁇ , lateral acceleration Gy, and front wheel steering angle ⁇ are ⁇ 0, Gy0, and ⁇ 0, respectively, the detected values of the vehicle yaw rate, lateral acceleration, and front wheel steering angle are ⁇ + ⁇ 0, Gy + Gy0 and ⁇ + ⁇ 0.
  • the design value and the true value of the stability factor are Khde and Khre, respectively, and the deviation ⁇ t between the steady reference yaw rate ⁇ t and the detected yaw rate at the time of steady turning of the vehicle is expressed by the following Expression 20.
  • the yaw rate deviation ⁇ t is converted into the steering wheel deviation ⁇ t by multiplying both sides of the above equation 20 by L / V
  • the steering angle deviation ⁇ t of the front wheel is expressed by the following equation 21. Therefore, the influence of the zero offset of the sensor is ⁇ 0 ⁇ KhdeGy0L ⁇ 0L / V.
  • ⁇ 0 ⁇ KhdeGy0L is a constant, but ⁇ 0L / V changes according to the vehicle speed V.
  • the cutoff frequency fhc is based on the absolute value of the longitudinal acceleration Gx of the vehicle so that the cutoff frequency fhc of the high-pass filter processing increases as the absolute value of the longitudinal acceleration Gx of the vehicle increases.
  • FIG. 7 is a flowchart showing a main part of an estimation calculation routine of the steering response time constant coefficient Tp in the fourth embodiment of the turning characteristic estimation device according to the present invention configured as a modification of the first embodiment.
  • step 95 is executed prior to step 100.
  • step 95 it is determined whether or not the vehicle is in a state where the steering response time constant coefficient Tp can be estimated with high reliability. If a negative determination is made, the control returns to step 20 and affirmative. When the determination is made, control proceeds to step 100.
  • C1 The absolute value of the steering angular velocity ⁇ d, which is the time differential value of the steering angle ⁇ , is greater than or equal to the reference value of the steering angular velocity.
  • the traveling road is not a bad road.
  • C3 Not braking.
  • the condition of C1 is that the steering response time constant coefficient Tp represents the characteristics of the transient turning of the vehicle. In a situation where the magnitude of the steering angular velocity is small, the steering response time constant coefficient Tp can be accurately estimated. Considering what cannot be done.
  • the condition of C2 takes into consideration that noise is superimposed on the actual yaw rate ⁇ on a rough road and that the grip state of the tire with respect to the road surface is likely to fluctuate.
  • the condition of C3 takes into consideration that the calculation of the steady-state standard yaw rate ⁇ t according to the above equation 11 is premised on that there is no influence of the braking force.
  • FIGS. 8 and 9 are respectively a first half and a second half of an estimation calculation routine for the steering response time constant coefficient Tp in the fifth embodiment of the turning characteristic estimation device according to the present invention configured as a modification of the first embodiment. It is a flowchart which shows a part. In the fifth embodiment, when step 50 is completed, it is determined in step 52 whether or not the vehicle is turning right.
  • steps 60 to 180 the same control as in steps 60 to 180 of the first embodiment is performed for the right turn of the vehicle.
  • step 80 high-pass filtering is performed on the actual yaw rate ⁇ and the transient yaw rate ⁇ trrn.
  • step 140 the smallest value ⁇ arm (m is any one of 1 to 5) among the integrated values ⁇ ar1 to ⁇ ar5 of the five yaw rate deviation index values is determined. Further, the steering response time constant coefficient Tprm corresponding to the minimum value ⁇ arm of the integrated values of the yaw rate deviation index values is set as the estimated value of the steering response time constant coefficient Tpr when turning right.
  • step 160 the reference value interval ⁇ Tpr is updated to a larger value of Ka ⁇ the current reference value interval ⁇ Tpr and the guard minimum value ⁇ Tpmin, and stored in the buffer memory. Further, when an affirmative determination is made in step 170, the estimated value of the steering response time constant coefficient Tpr at the time of right turn is stored in the EEPROM in step 180, and thereby the right turn stored in the EEPROM. The estimated value of the steering response time constant coefficient Tpr is updated. In steps 65 to 185, “r” indicating a right turn is replaced with “l” indicating a left turn, whereby the same control as in steps 60 to 180 is performed for the left turn. ⁇ Turning characteristics when turning right and turning characteristics when turning left may be different.
  • the vehicle turns The turning characteristics of the vehicle vary depending on the direction. Therefore, the steering response time constant coefficient also varies depending on the turning direction of the vehicle. According to the fifth embodiment, the same operational effects as those of the first embodiment described above are obtained, and the steering response time constant coefficient is estimated for each turning direction of the vehicle. Even when the turning characteristics are different, the steering response time constant coefficients Tpr and Tpl can be accurately estimated for each of the right turn and the left turn of the vehicle.
  • step 10 and 11 are respectively a first half and a second half of an estimation calculation routine of the steering response time constant coefficient Tp in the sixth embodiment of the turning characteristic estimation device according to the present invention configured as a modification of the first embodiment. It is a flowchart which shows a part.
  • step 50 when step 50 is completed, it is determined in step 54 whether or not the absolute value of the lateral acceleration Gy of the vehicle is larger than the first reference value Gy1 (positive constant). Is done. When a negative determination is made, the control returns to step 20, and when an affirmative determination is made, the control proceeds to step 56. In step 56, it is determined whether or not the absolute value of the lateral acceleration Gy of the vehicle is larger than a second reference value Gy2 (a positive constant larger than Gy1).
  • step 80 a high-pass filter process is performed on the actual yaw rate ⁇ and the transient yaw rate ⁇ tr1n.
  • step 140 the smallest value ⁇ a1m (m is any one of 1 to 5) among the integrated values ⁇ a11 to ⁇ a15 of the five yaw rate deviation index values is determined.
  • the steering response time constant coefficient Tp1m corresponding to the minimum value ⁇ a1m of the integrated values of the yaw rate deviation index values is the estimated value of the steering response time constant coefficient Tp1 for the first region of the lateral acceleration Gy.
  • the reference value interval ⁇ Tp1 is updated to a larger value of Ka ⁇ the current reference value interval ⁇ Tp1 and the guard minimum value ⁇ Tpmin, and stored in the buffer memory.
  • step 170 the estimated value of the steering response time constant coefficient Tp1 is stored in the EEPROM in step 180, thereby the first region of the lateral acceleration Gy stored in the EEPROM.
  • the estimated value of the steering response time constant coefficient Tp1 is updated.
  • steps 65 to 185 “1” indicating the first region of the lateral acceleration Gy is replaced with “2” indicating the second region of the lateral acceleration Gy, whereby the second acceleration of the lateral acceleration Gy.
  • the same control as in steps 60 to 180 is performed for the area.
  • the turning characteristics of a vehicle vary depending on the magnitude of the lateral acceleration Gy of the vehicle.
  • FIG. 12 is a flowchart showing an estimation calculation routine of the steering response time constant coefficient Tp in the seventh embodiment of the turning characteristic estimation device according to the present invention applied to the vehicle motion control device. In the seventh embodiment, step 60 and steps 100 to 180 in the first embodiment described above are not executed.
  • step 70 the transient yaw rate ⁇ tr is calculated according to the following equation 22 corresponding to the above equation 17. If an affirmative determination is made in step 90, the control proceeds to step 200, and the steering response time constant coefficient Tp is estimated according to the flowchart shown in FIG. As shown in FIG. 13, in step 210, whether or not the absolute value of the phase difference Dy between the transient yaw rate ⁇ trbpf and the actual yaw rate ⁇ bpf after the bandpass filter processing is larger than a reference value Dy0 (positive constant). A determination is made. When a negative determination is made, it is not necessary to correct the increase or decrease of the estimated value of the steering response time constant coefficient Tp.
  • step 220 it is determined whether or not the phase of the transient yaw rate ⁇ trbpf is earlier than the phase of the actual yaw rate ⁇ bpf. If a negative determination is made, the control proceeds to step 240, and if an affirmative determination is made, the control is performed. Proceeds to step 230. In step 230, the sum Tpf + ⁇ Tpb of the previous estimated value Tpf of the steering response time constant coefficient Tp and the increase / decrease correction amount ⁇ Tpb set in the previous step 260 and the preset maximum value of the steering response time constant coefficient.
  • Tpmax The smaller value of the steering response time constant coefficient Tp.
  • a value Tpf ⁇ Tpb obtained by subtracting the increase / decrease correction amount ⁇ Tpb from the previous estimated value Tpf of the steering response time constant coefficient Tp and a preset minimum value Tpmin (Tpmax smaller than Tpmax).
  • Tpmin a preset minimum value
  • Tpmin a preset minimum value
  • Tpe positive constant
  • step 260 the “previous estimated value Tpf” for the control of the next cycle is rewritten to the estimated value of the steering response time constant coefficient Tp, and Kb is larger than 0 and larger than 1, for example 0.5.
  • the increase / decrease correction amount ⁇ Tpb is rewritten to the product of the coefficient Kb and the current increase / decrease correction amount ⁇ Tpb.
  • step 270 the estimated value of the rudder response time constant coefficient Tp is stored in the EEPROM, and in step 280, the previous estimated value Tpf is cleared to 0 and the increase / decrease correction amount ⁇ Tpb is cleared to 0. .
  • the seventh embodiment when the phase of the transient yaw rate ⁇ trbpf is earlier than the phase of the actual yaw rate ⁇ bpf, an affirmative determination is made at step 220, and an estimated value of the steering response time constant coefficient Tp at step 230. Is corrected to increase.
  • the phase of the transient yaw rate ⁇ trbpf is later than the phase of the actual yaw rate ⁇ bpf, a negative determination is made at step 220, and the estimated value of the steering response time constant coefficient Tp is reduced and corrected at step 240.
  • the estimated value of the steering response time constant coefficient Tp approaches the true steering response time constant coefficient by reducing the phase difference between the transient yaw rate ⁇ tr of the vehicle and the actual yaw rate ⁇ of the vehicle.
  • the estimated value of the steering response time constant coefficient Tp can be corrected. Therefore, the steering response time constant coefficient Tp can be estimated to a value close to the true steering response time constant coefficient.
  • the steady-state standard yaw rate ⁇ t is calculated based on the steering angle ⁇ and the like subjected to the low-pass filter processing in step 30.
  • step 80 the high-pass filter process is performed on the transient yaw rate ⁇ tr and the actual yaw rate ⁇ calculated based on the steady reference yaw rate ⁇ t in step 70. Accordingly, as in the case of the first to sixth embodiments described above, not only the high-frequency noise included in the detected steering angle ⁇ and the like can be removed, but also the influence of the zero point offset of the yaw rate sensor 36 and the like is removed. Can do.
  • step 90 it is determined in step 90 whether or not the vehicle is in a turning traveling state, and when an affirmative determination is made, step 200 and subsequent steps are executed. Accordingly, as in the case of the first to sixth embodiments described above, step 200 and subsequent steps are unnecessarily executed in a situation where the vehicle is not in a turning state and the steering response time constant coefficient Tp cannot be accurately estimated.
  • step 260 when the absolute value of the deviation between the steering response time constant coefficient Tp and the previous estimated value Tpf is equal to or greater than the reference value Tpe for storage determination, a negative determination is made in step 250. .
  • the increase / decrease correction amount ⁇ Tpb is reduced by rewriting the product of the coefficient Kb and the current increase / decrease correction amount ⁇ Tpb. Therefore, the phase difference between the transient yaw rate ⁇ tr and the actual yaw rate ⁇ can be effectively reduced as compared with the case where the increase / decrease correction amount ⁇ Tpb is not reduced, thereby effectively reducing the estimated value of the steering response time constant coefficient Tp.
  • step 270 the estimated value of the rudder response time constant coefficient Tp is stored in the EEPROM and updated. Therefore, the estimated value of the steering response time constant coefficient Tp becomes the EEPROM at the stage where the estimated value of the steering response time constant coefficient Tp does not become a value close to the true steering response time constant coefficient, compared with the case where the determination in step 250 is not performed. And the frequency of storing and updating the estimated value of the rudder response time constant coefficient Tp can be reduced.
  • step 210 when the absolute value of the phase difference Dy between the transient yaw rate ⁇ trbpf and the actual yaw rate ⁇ bpf after the bandpass filter processing is smaller than the reference value Dy0, a negative determination is made in step 210, and the step The steps after 220 are not executed. Therefore, in the situation where the estimated value of the steering response time constant coefficient Tp is close to the true steering response time constant coefficient, it is possible to prevent the steering response time constant coefficient Tp from being estimated and updated unnecessarily. Can do. Eighth embodiment FIG.
  • step 14 is a flowchart showing an estimation calculation routine of the steering response time constant coefficient Tp in the eighth embodiment of the turning characteristic estimation device according to the present invention configured as a modification of the seventh embodiment.
  • step 52 is executed as in the case of the fifth embodiment described above. If an affirmative determination is made in step 52, the control proceeds to step 70, and if a negative determination is made, the control proceeds to step 75.
  • step 70 the transient yaw rate ⁇ trr during the right turn is calculated according to the following equation 23 corresponding to the equation 17 based on the steering response time constant coefficient Tpr during the right turn described later.
  • step 80 the high-pass filter process is performed on the actual yaw rate ⁇ and the transient yaw rate ⁇ trr, whereby the actual yaw rate ⁇ bpf and the transient yaw rate ⁇ trrbpf after the high-pass filter process are calculated. If an affirmative determination is made in step 90, the control proceeds to step 200, where the steering response time constant coefficient Tpr during right turn is estimated according to the flowchart shown in FIG. 15, and the steering response time constant coefficient Tpr is calculated. The estimated value is stored in the EEPROM.
  • step 75 based on the steering response time constant coefficient Tpl for the left turn described later, the transient yaw rate ⁇ trl for the left turn is calculated according to the following equation 24 corresponding to the above equation 17.
  • step 85 the high-pass filter process is performed on the actual yaw rate ⁇ and the transient yaw rate ⁇ trl, thereby calculating the actual yaw rate ⁇ bpf and the transient yaw rate ⁇ trlbpf after the high-pass filter process.
  • step 95 is executed in the same manner as step 90.
  • step 95 the control proceeds to step 205, where the steering response time constant coefficient Tpl during left turn is estimated according to the flowchart shown in FIG. 16, and the estimated value of the steering response time constant coefficient Tpl is determined. Is stored in the EEPROM. As shown in FIG. 15, in step 210, it is determined whether or not the absolute value of the phase difference Dyr between the transient yaw rate ⁇ trrbpf and the actual yaw rate ⁇ bpf when turning right after the bandpass filter processing is larger than the reference value Dy0. A determination is made. When a negative determination is made, it is not necessary to correct the increase or decrease of the estimated value of the steering response time constant coefficient Tpr.
  • step 220 it is determined whether or not the phase of the transient yaw rate ⁇ trrbpf is earlier than the phase of the actual yaw rate ⁇ bpf. If a negative determination is made, the control proceeds to step 240, and if an affirmative determination is made, the control is performed. Proceeds to step 230. In step 230, the sum Tprf + ⁇ Tprb of the previous estimated value Tprf of the steering response time constant coefficient Tpr for the right turn and the increase / decrease correction amount ⁇ Tprb for the right turn set in the previous step 260 is set in advance.
  • the smaller one of the maximum steering response time constant coefficients Tpmax is set as the estimated value of the steering response time constant coefficient Tpr when turning right.
  • the value Tprf ⁇ Tprb obtained by subtracting the increase / decrease correction amount ⁇ Tprb from the previous estimated value Tprf of the steering response time constant coefficient Tpr for the right turn and the preset minimum value Tpmin of the steering response time constant coefficient are set.
  • the larger value is the estimated value of the steering response time constant coefficient Tpr when turning right.
  • step 270 the control proceeds to step 270, and when a negative determination is made, the control proceeds to step 260.
  • step 260 the “previous estimated value Tprf” for the control of the next cycle is rewritten to the estimated value of the steering response time constant coefficient Tpr, and the increase / decrease correction amount ⁇ Tprb is changed to the coefficient Kb and the current increase / decrease correction amount ⁇ Tprb.
  • step 270 the estimated value of the rudder response time constant coefficient Tpr is stored in the EEPROM, and in step 280, the previous estimated value Tprf is cleared to 0 and the increase / decrease correction amount ⁇ Tprb is cleared to 0. .
  • FIG. 17 is a flowchart showing an estimation calculation routine of the steering response time constant coefficient Tp in the ninth embodiment of the turning characteristic estimation device according to the present invention configured as a modification of the seventh embodiment.
  • step 50 when step 50 is completed, steps 54 and 56 are executed as in the case of the sixth embodiment described above.
  • step 56 control proceeds to step 70, and when a negative determination is made, control proceeds to step 75.
  • step 70 the transient yaw rate ⁇ tr1 for the first region of the lateral acceleration Gy is calculated according to the following equation 25 corresponding to the above equation 17.
  • step 80 the high-pass filter process is performed on the actual yaw rate ⁇ and the transient yaw rate ⁇ tr1, thereby calculating the actual yaw rate ⁇ bpf and the transient yaw rate ⁇ tr1bpf after the high-pass filter process.
  • step 90 the control proceeds to step 200, and the steering response time constant coefficient Tp1 for the first region of the lateral acceleration Gy is estimated according to the flowchart shown in FIG.
  • the estimated value of the time constant coefficient Tp1 is stored in the EEPROM.
  • step 75 the transient yaw rate ⁇ tr2 for the second region of the lateral acceleration Gy is calculated according to the following equation 26 corresponding to the above equation 17.
  • step 85 the high-pass filter process is performed on the actual yaw rate ⁇ and the transient yaw rate ⁇ tr2, thereby calculating the actual yaw rate ⁇ bpf and the transient yaw rate ⁇ tr2bpf after the high-pass filter process.
  • step 95 is executed in the same manner as step 90. If an affirmative determination is made in step 95, the control proceeds to step 205, where the steering response time constant coefficient Tp2 for the second region of the lateral acceleration Gy is estimated according to the flowchart shown in FIG. The estimated value of the constant coefficient Tp2 is stored in the EEPROM. As shown in FIG. 18, in step 210, the absolute value of the phase difference Dy1 between the transient yaw rate ⁇ tr1bpf and the actual yaw rate ⁇ bpf after the bandpass filter processing for the first region of the lateral acceleration Gy is greater than the reference value Dy0. A determination is made as to whether it is larger.
  • step 220 it is determined whether or not the phase of the transient yaw rate ⁇ tr1bpf is earlier than the phase of the actual yaw rate ⁇ bpf. If a negative determination is made, the control proceeds to step 240, and if an affirmative determination is made, the control is performed. Proceeds to step 230.
  • step 230 the previous estimated value Tp1f of the steering response time constant coefficient Tp1 for the first region of the lateral acceleration Gy and the increase / decrease of the lateral acceleration Gy set in the previous step 260 for the first region.
  • the smaller one of the sum Tp1f + ⁇ Tp1b with the correction amount ⁇ Tp1b and the preset maximum value Tpmax of the steering response time constant coefficient is the estimated value of the steering response time constant coefficient Tp1 for the first region of the lateral acceleration Gy. .
  • step 240 a value Tp1f ⁇ Tp1b obtained by subtracting the increase / decrease correction amount ⁇ Tp1b from the previous estimated value Tp1f of the steering response time constant coefficient Tp1 for the first region of the lateral acceleration Gy and a preset steering response time constant coefficient.
  • the larger one of the minimum values Tpmin is set as the estimated value of the steering response time constant coefficient Tp1 for the first region of the lateral acceleration Gy.
  • step 270 the control proceeds to step 270, and when a negative determination is made, the control proceeds to step 260.
  • step 260 the “previous estimated value Tp1f” for the control of the next cycle is rewritten with the estimated value of the steering response time constant coefficient Tp1, and the increase / decrease correction amount ⁇ Tp1b is changed to the coefficient Kb and the current increase / decrease correction amount ⁇ Tp1b.
  • step 270 the estimated value of the rudder response time constant coefficient Tp1 is stored in the EEPROM.
  • step 280 the previous estimated value Tp1f is cleared to 0, and the increase / decrease correction amount ⁇ Tp1b is cleared to 0. .
  • steps 215 to 285 “1” indicating the first region of the lateral acceleration Gy is replaced with “2” indicating the second region of the lateral acceleration Gy, whereby the second acceleration of the lateral acceleration Gy.
  • the same control as in steps 210 to 280 is performed for the area. Therefore, according to the ninth embodiment, the same effects as those of the seventh embodiment described above can be obtained, and the magnitude of the lateral acceleration Gy can be increased even when the magnitude of the lateral acceleration Gy of the vehicle varies greatly.
  • the steering response time constant coefficient Tp can be accurately estimated for each region.
  • the transient yaw rate ⁇ trbpf and the like are calculated for each cycle as in the other embodiments. However, the steering read for each cycle is used.
  • the angle ⁇ and the like are stored in the buffer memory, and the transient yaw rate ⁇ trbpf and the like are calculated on the basis of the steering angle ⁇ and the like of each cycle after the turn is completed. It may be modified so that the process 205 is performed.
  • the present invention has been described in detail with respect to specific embodiments, the present invention is not limited to the above-described embodiments, and various other embodiments are possible within the scope of the present invention. It will be apparent to those skilled in the art.
  • the reference value may be modified to be calculated in accordance with, for example, Equations 13 to 15, and may be modified to calculate seven or more reference values. .
  • the reference value interval ⁇ Tp and the like are gradually reduced, but the reference value interval ⁇ Tp and the like may be constant values.
  • the integrated value ⁇ an or the like of the yaw rate deviation index value is cleared to zero. Is done.
  • the integrated value ⁇ an of the yaw rate deviation index value may be corrected so as to be adjusted to a positive value smaller than the current value, for example, a positive value close to 0.
  • the steering response time constant coefficient Tp is estimated for the first and second regions where the magnitude of the lateral acceleration Gy of the vehicle is different.
  • the steering response time constant coefficient Tp may be corrected so as to be estimated for three or more regions in which the lateral acceleration Gy of the vehicle is different.
  • the increase / decrease correction amount ⁇ Tpb or the like is gradually reduced, but the increase / decrease correction amount ⁇ Tpb or the like may be a constant value.
  • the fifth and sixth embodiments described above are configured as modifications of the first embodiment described above, but the configuration of the fifth or sixth embodiment is any of the second to fourth embodiments. May be applied.
  • the seventh embodiment described above is configured as a modification of the first embodiment described above, but the configuration of the seventh embodiment is applied to any of the second to fourth embodiments. Also good.

Landscapes

  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Steering Control In Accordance With Driving Conditions (AREA)
  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)

Abstract

本発明による車両の旋回特性推定装置は、車両の定常規範ヨーレートに対し一次遅れの関係にある車両の過渡ヨーレートを車両の実ヨーレートに対応する値とみなし、一次遅れの時定数の車速にかかる係数である操舵応答時定数係数を推定する。操舵応答時定数係数の推定値を使用して車両の定常規範ヨーレートに基づいて車両の過渡ヨーレートを演算し、車両の過渡ヨーレートが車両の実ヨーレートに近づくよう車両の過渡ヨーレートと車両の実ヨーレートとの関係に基づいて操舵応答時定数係数の推定値を修正する。

Description

車両の旋回特性推定装置
 本発明は、車両の旋回特性推定装置に係り、更に詳細には車両が旋回する際の規範ヨーレートと過渡ヨーレートとの関係に基づき車両の旋回特性を示す操舵応答時定数係数を推定する車両の旋回特性推定装置に係る。
 車両の規範ヨーレートに対し車両の実ヨーレートが一次遅れの関係にあり、一次遅れの時定数の車速にかかる係数を操舵応答時定数係数とすると、車両のスタビリティファクタ及び操舵応答時定数係数は車両の旋回特性を表す。ARX(auto−regressive exogenous model)を使用し、規範ヨーレートから実ヨーレートへの離散時間伝達関数のパラメータa及びbを推定することにより、車両のスタビリティファクタ及び操舵応答時定数係数を推定することができる。
 例えば特開2004−26073号公報には、車両が旋回する際の走行データに基づき規範ヨーレートを推定し、規範ヨーレートから実ヨーレートへの離散時間伝達関数のパラメータa及びbを推定し、パラメータaに基づいて操舵応答時定数係数を推定し、パラメータa及びbに基づいて、車両のスタビリティファクタを推定する車両の旋回特性推定装置が記載されている。
 しかし上記公開公報に記載されている如き従来の旋回特性推定装置に於いては、パラメータa及びbを推定するためには規範ヨーレート及び実ヨーレートについての多数のデータが必要である。そのため車両のスタビリティファクタや操舵応答時定数係数の推定に長い時間を要するという問題がある。
 本発明の主要な目的は、従来の旋回特性推定装置に於ける上述の如き問題に鑑み、車両の定常規範ヨーレートに対し一次遅れの関係にある車両の過渡ヨーレート及び車両の実ヨーレートに基づいて操舵応答時定数係数の推定値が真の操舵応答時定数係数に近付くよう操舵応答時定数係数の推定値を学習によって修正することにより、従来に比して短時間にて操舵応答時定数係数を推定することである。
 本発明によれば、車両の定常規範ヨーレートに対し一次遅れの関係にある車両の過渡ヨーレートを車両の実ヨーレートに対応する値とみなし、一次遅れの時定数の車速にかかる係数である操舵応答時定数係数を推定する車両の旋回特性推定装置であって、操舵応答時定数係数の推定値を使用して車両の定常規範ヨーレートに基づいて車両の過渡ヨーレートを演算し、車両の過渡ヨーレートが車両の実ヨーレートに近づくよう車両の過渡ヨーレートと車両の実ヨーレートとの関係に基づいて操舵応答時定数係数の推定値を修正する車両の旋回特性推定装置が提供される。
 この構成によれば、車両の過渡ヨーレートが車両の実ヨーレートに近づくよう操舵応答時定数係数の推定値が修正されるので、操舵応答時定数係数の推定値を真の操舵応答時定数係数に近づくよう修正することができる。従って規範ヨーレートから実ヨーレートへの離散時間伝達関数のパラメータを推定することは不要であり、従来に比して短時間にて操舵応答時定数係数を推定することができる。
 上記構成に於いて、車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさ及び車両の過渡ヨーレートと車両の実ヨーレートとの位相差の少なくとも一方が小さくなることによって車両の過渡ヨーレートが車両の実ヨーレートに近づくよう操舵応答時定数係数の推定値を修正するようになっていてよい。
 この構成によれば、車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさ及び車両の過渡ヨーレートと車両の実ヨーレートとの位相差の少なくとも一方が小さくなり、これにより車両の過渡ヨーレートが車両の実ヨーレートに近づくよう操舵応答時定数係数の推定値を修正することができる。
 また上記構成に於いて、記憶されている操舵応答時定数係数を第一の基準値として、第一の基準値、第一の基準値よりも大きい第二の基準値、第一の基準値よりも小さい第三の基準値を使用して三つの車両の過渡ヨーレートを演算し、第一乃至第三の基準値のうち、車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさの指標値であるヨーレート偏差指標値が最も小さくなる基準値を操舵応答時定数係数の推定値とするようになっていてよい。
 この構成によれば、第一乃至第三の基準値のうちヨーレート偏差指標値が最も小さくなる基準値が操舵応答時定数係数の推定値とされる。よって第一乃至第三の基準値のうち真の操舵応答時定数係数に最も近い値を操舵応答時定数係数の推定値とすることができる。
 また上記構成に於いて、操舵応答時定数係数の推定値を新たな第一の基準値とし、新たな第一の基準値、新たな第一の基準値よりも大きい新たな第二の基準値、新たな第一の基準値よりも小さい新たな第三の基準値を使用して三つの車両の過渡ヨーレートを演算し、新たな第一乃至第三の基準値のうち、ヨーレート偏差指標値が最も小さくなる基準値を新たな操舵応答時定数係数の推定値とするようになっていてよい。
 この構成によれば、新たな第一乃至第三の基準値のうち真の操舵応答時定数係数に最も近い値が操舵応答時定数係数の推定値とされる。従って第一乃至第三の基準値より選択される操舵応答時定数係数の推定値よりも真の操舵応答時定数係数に近い操舵応答時定数係数の推定値を求めることができる。よってこの構成の手順を繰り返すことにより、操舵応答時定数係数の推定値を真の操舵応答時定数係数に漸次近づけることができる。
 また上記構成に於いて、新たな第二の基準値と新たな第一の基準値との差は第二の基準値と第一の基準値との差よりも小さく、新たな第三の基準値と新たな第一の基準値との差は第三の基準値と第一の基準値との差よりも小さくてよい。
 この構成によれば、新たな第二の基準値と新たな第一の基準値との差が第二の基準値と第一の基準値との差以上であり、新たな第三の基準値と新たな第一の基準値との差が第三の基準値と第一の基準値との差以上である場合に比して、操舵応答時定数係数の推定値を速く真の操舵応答時定数係数に近づけることができる。
 また上記構成に於いて、ヨーレート偏差指標値は、第一の所定周波数以下の成分が除去された車両の実ヨーレートと、第二の所定周波数以下の成分が除去された車両の過渡ヨーレートとの偏差の大きさに基づく値であってよい。
 過渡ヨーレートは車両の状態量に基づいて求められ、それらの車両の状態量や車両の実ヨーレートは検出手段によって検出される。上記構成によれば、車両の状態量や車両の実ヨーレートを検出する検出手段の零点オフセットの如き定常的な検出誤差を排除することができるので、操舵応答時定数係数の推定精度を高くすることができる。
 上記構成に於いて、単位時間当たりの運転者による往復操舵の多さの指標値に応じて第一の所定周波数若しくは第二の所定周波数を変更するようになっていてよい。
 実ヨーレート等の車両の状態量を検出する検出手段の零点オフセットの如き定常的な検出誤差は、単位時間当たりの運転者による往復操舵の多さによって変動する。上記構成によれば、単位時間当たりの運転者による往復操舵の多さに応じて適正に定常的な検出誤差を排除することができる。
 また上記構成に於いて、車両の前後加速度の大きさに応じて第一の所定周波数若しくは第二の所定周波数を変更するようになっていてよい。
 実ヨーレート等の車両の状態量を検出する検出手段の零点オフセットの如き定常的な検出誤差が操舵応答時定数係数の推定精度に与える影響は、車速の変化の大きさ、即ち車両の前後加速度の大きさによっても変動する。上記構成によれば、車両の前後加速度の大きさに応じて適正に定常的な検出誤差を排除することができる。
 また上記構成に於いて、ヨーレート偏差指標値は車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさを前輪の舵角の偏差の大きさに換算した値であってよい。
 車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさは車速に依存するが、車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさを前輪の舵角の偏差の大きさに換算した値は車速に依存しない。従って上記構成によれば、車速に依存しないヨーレート偏差指標値に基づいて操舵応答時定数係数の推定値を求めることができ、これにより車速の影響を受けることなく操舵応答時定数係数を推定することができる。尚前輪の舵角の偏差とは、車両の過渡ヨーレートを達成するための前輪の舵角と前輪の実舵角との差である。
 また上記構成に於いて、車両の過渡ヨーレートの位相が車両の実ヨーレートの位相よりも早いときには操舵応答時定数係数の推定値を増大修正し、車両の過渡ヨーレートの位相が車両の実ヨーレートの位相よりも遅いときには操舵応答時定数係数の推定値を低減修正するようになっていてよい。
 この構成によれば、車両の過渡ヨーレートと車両の実ヨーレートとの位相差が小さくなることによって操舵応答時定数係数の推定値が真の操舵応答時定数係数に近づくよう操舵応答時定数係数の推定値を修正することができる。
 また上記構成に於いて、車両の右旋回及び左旋回について個別に操舵応答時定数係数を推定するようになっていてよい。
 この構成によれば、例えば車両の重心が車両の横方向の中心にない場合や車両の横方向に変動する場合の如く、車両の旋回方向によって車両の旋回特性が異なる場合にも、車両の右旋回及び左旋回の各々について操舵応答時定数係数を推定することができる。
 また上記構成に於いて、車両の横加速度の大きさの領域毎に操舵応答時定数係数を推定するようになっていてよい。
 車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさは車両の横加速度の大きさによっても変動する。上記構成によれば、車両の横加速度の大きさの領域毎に操舵応答時定数係数を推定し、これにより車両の横加速度の大きさの影響を排除して操舵応答時定数係数を推定することができる。
 図20に示された車両の二輪モデルに於いて、車両の質量及びヨー慣性モーメントをそれぞれM及びIとし、車両の重心102と前輪車軸及び後輪車軸との間の距離をそれぞれLf及びLrとし、車両のホイールベースをL(=Lf+Lr)とする。また前輪100f及び後輪100rのコーナリングフォースをそれぞれFf及びFrとし、前輪及び後輪のコーナリングパワーをそれぞれKf及びKrとする。また前輪100fの実舵角をδとし、前輪及び後輪のスリップ角をそれぞれβf及びβrとし、車体のスリップ角をβとする。更に車両の横加速度をGyとし、車両のヨーレートをγとし、車速をVとし、車両のヨー加角速度(ヨーレートγの微分値)をγdとする。車両の力及びモーメントの釣合い等により下記の式1~6が成立する。
 MGy=Ff+Fr  ……(1)
 Iγd=LfFf−LrFr  ……(2)
 Ff=−Kfβf  ……(3)
 Fr=−Krβr  ……(4)
 βf=β+(Lf/V))γ−δ  ……(5)
 βr=β−(Lr/V)γ  ……(6)
 上記式1~6より下記の式7が成立する。
Figure JPOXMLDOC01-appb-I000001
 車速Vが実質的に一定であると仮定し、ラプラス演算子をsとして上記式7をラプラス変換し、ヨーレートγについて整理することにより、下記の式8~10が成立し、よってこれらの式により規範ヨーレートγ(s)が求められる。
Figure JPOXMLDOC01-appb-I000002
 上記式9のKhはスタビリティファクタであり、上記式10のTpは車速依存の時定数をもつ一次遅れ系の車速Vにかかる係数、即ち本明細書に於いて「操舵応答時定数係数」と呼ぶ係数である。これらの値は車両のヨー運動に関する操舵応答を特徴付けるパラメータ、即ち車両の旋回特性である。また上記式8は前輪の実舵角δ、車速V、横加速度Gyより車両のヨーレートγを演算する式である。この線形化モデルより演算されるヨーレートを下記の式11にて表される定常規範ヨーレートγtに対する一次遅れの過渡ヨーレートγtrとする。
Figure JPOXMLDOC01-appb-I000003
 よって上記構成に於いて、上記式8に従って過渡ヨーレートγtrが演算され、過渡ヨーレートγtrが実ヨーレートγに近づくよう過渡ヨーレートγtrと実ヨーレートγとの関係に基づいて操舵応答時定数係数Tpの推定値が修正されてよい。
 また上記構成に於いて、ヨーレート偏差指標値は車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさの指標値の積算値であってよい。
 また上記構成に於いて、ハイパスフィルタ処理によって車両の実ヨーレートより第一の所定周波数以下の成分が除去され、ハイパスフィルタ処理によって車両の過渡ヨーレートより第二の所定周波数以下の成分が除去されてよい。
 また上記構成に於いて、第一及び第二の所定周波数は同一の周波数であってよい。
 また上記構成に於いて、車速をVとし、車両のホイールベースをLとして、車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさにL/Vが乗算されることにより、車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさが前輪の舵角の偏差の大きさに換算した値が演算されてよい。
 また上記構成に於いて、操舵応答時定数係数の推定値を増減修正量にて修正し、増減修正量の大きさを漸次低減するようになっていてよい。
 また上記構成に於いて、今回の操舵応答時定数係数の推定値と前回の操舵応答時定数係数の推定値との偏差の大きさが基準値未満になったときに操舵応答時定数係数の推定値の修正を終了するようになっていてよい。
 また上記構成に於いて、予め設定された操舵応答時定数係数の最小値と最大値との範囲を越えないよう操舵応答時定数係数の推定値を増減修正量にて修正するようになっていてよい。
 また上記構成に於いて、ハイパスフィルタ処理後の車両の過渡ヨーレートとハイパスフィルタ処理後の車両の実ヨーレートとの偏差の大きさが基準値未満であるときには操舵応答時定数係数の推定を行なわないようになっていてよい。
 図1は車両の運動制御装置に適用された本発明による旋回特性推定装置の第一の実施形態を示す概略構成図である。
 図2は第一の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンを示すフローチャートである。
 図3は第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第二の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの要部を示すフローチャートである。
 図4は操舵周波数fsとハイパスフィルタ処理のカットオフ周波数fhcとの関係を示すグラフである。
 図5は第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第三の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの要部を示すフローチャートである。
 図6は操舵周波数fsとハイパスフィルタ処理のカットオフ周波数fhcと車両の前後加速度Gxの絶対値との関係を示すグラフである。
 図7は第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第四の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの要部を示すフローチャートである。
 図8は第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第五の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの前半部を示すフローチャートである。
 図9は第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第五の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの後半部を示すフローチャートである。
 図10は第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第六の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの前半部を示すフローチャートである。
 図11は第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第六の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの後半部を示すフローチャートである。
 図12は本発明による旋回特性推定装置の第七の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンを示すフローチャートである。
 図13は図12のステップ200に於ける操舵応答時定数係数Tpの推定及び記憶のサブルーチンを示すフローチャートである。
 図14は本発明による旋回特性推定装置の第八の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンを示すフローチャートである。
 図15は図14のステップ200に於ける操舵応答時定数係数Tpの推定及び記憶のサブルーチンを示すフローチャートである。
 図16は図14のステップ205に於ける操舵応答時定数係数Tpの推定及び記憶のサブルーチンを示すフローチャートである。
 図17は本発明による旋回特性推定装置の第九の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンを示すフローチャートである。
 図18は図17のステップ200に於ける操舵応答時定数係数Tpの推定及び記憶のサブルーチンを示すフローチャートである。
 図19は図17のステップ205に於ける操舵応答時定数係数Tpの推定及び記憶のサブルーチンを示すフローチャートである。
 図20は操舵応答時定数係数を推定するための車両の二輪モデルを示す説明図である。
 以下に添付の図を参照しつつ、本発明を好ましい実施形態について詳細に説明する。
第一の実施形態
 図1は車両の運動制御装置に適用された本発明による旋回特性推定装置の第一の実施形態を示す概略構成図である。
 図1に於いて、50は車両10の運動制御装置を全体的に示しており、本発明による旋回特性推定装置は運動制御装置50の一部をなしている。車両10は左右の前輪12FL及び12FR及び左右の後輪12RL及び12RRを有している。操舵輪である左右の前輪12FL及び12FRは運転者によるステアリングホイール14の転舵に応答して駆動されるラック・アンド・ピニオン式のパワーステアリング装置16によりタイロッド18L及び18Rを介して操舵される。
 各車輪の制動力は制動装置20の油圧回路22によりホイールシリンダ24FR、24FL、24RR、24RLの制動圧が制御されることによって制御されるようになっている。図には示されていないが、油圧回路22はオイルリザーバ、オイルポンプ、種々の弁装置等を含み、各ホイールシリンダの制動圧は通常時には運転者によるブレーキペダル26の踏み込み操作に応じて駆動されるマスタシリンダ28により制御され、また必要に応じて後に説明する如く電子制御装置30により制御される。
 車輪12FR~12RLのホイールシリンダにはそれぞれ対応するホイールシリンダの圧力Pi(i=fr、fl、rr、rl)を検出する圧力センサ32FR~32RLが設けられ、ステアリングホイール14が連結されたステアリングコラムには操舵角θを検出する操舵角センサ34が設けられている。
 また車両10にはそれぞれ車両の実ヨーレートγを検出するヨーレートセンサ36、車両の前後加速度Gxを検出する前後加速度センサ38、車両の横加速度Gyを検出する横加速度センサ40、各車輪の車輪速度Vwi(i=fr、fl、rr、rl)を検出する車輪速度センサ42FR~42RLが設けられている。尚操舵角センサ34、ヨーレートセンサ36及び横加速度センサ40は車両の左旋回方向を正としてそれぞれ操舵角、実ヨーレート及び横加速度を検出する。
 図示の如く、圧力センサ32FR~32RLにより検出された圧力Piを示す信号、操舵角センサ34により検出された操舵角θを示す信号、ヨーレートセンサ36により検出された実ヨーレートγを示す信号、前後加速度センサ38により検出された前後加速度Gxを示す信号、横加速度センサ40により検出された横加速度Gyを示す信号、車輪速度センサ42FR~42RLにより検出された車輪速度Vwiを示す信号は電子制御装置30に入力される。
 尚図には詳細に示されていないが、電子制御装置30は例えばCPUとROMとEEPROMとRAMとバッファメモリと入出力ポート装置とを有し、これらが双方向性のコモンバスにより互いに接続された一般的な構成のマイクロコンピュータを含んでいる。ROMは規範ヨーレートγtの演算に使用されるスタビリティファクタKh及び操舵応答時定数係数Tpのデフォルト値を記憶している。これらのデフォルト値は車両の出荷時に車両毎に設定される。またEEPROMは操舵応答時定数係数Tpの推定値等を記憶し、操舵応答時定数係数Tpの推定値等は後に詳細に説明する如く車両が旋回状態にあるときの車両の走行データに基づいて演算されることによって適宜更新される。
 電子制御装置30は、後述の如く図2に示されたフローチャートに従い、車両が旋回を開始すると、操舵角の如き旋回走行データに基づいて定常規範ヨーレートγtを演算し、大きさが互いに異なる操舵応答時定数係数Tpの五つの基準値について、定常規範ヨーレートγtに対し一次遅れの関係をなす五つの過渡ヨーレートγtrn(n=1~5)を演算する。
 そして電子制御装置30は、車両の実ヨーレートγと過渡ヨーレートγtrnとの偏差の大きさが最も小さくなる操舵応答時定数係数Tpの基準値に基づいて、予め設定された条件が成立しているときに操舵応答時定数係数Tpの推定値を演算し、それらをバッファメモリに記憶する。また電子制御装置30は、操舵応答時定数係数Tpの推定値を演算すると、その推定値に基づいて操舵応答時定数係数Tpの五つの基準値を演算するための標準値Tp0及び基準値間隔ΔTpを更新する。
 更に電子制御装置30は、EEPROMに記憶されている操舵応答時定数係数Tpの推定値を使用して過渡ヨーレートγtrに対応する目標ヨーレートγttを演算し、ヨーレート検出値γと目標ヨーレートγttとの偏差としてヨーレート偏差Δγを演算する。そして電子制御装置30は、ヨーレート偏差Δγの大きさが上記基準値γo(正の値)を越えているか否かの判別により車両の旋回挙動が悪化しているか否かを判定し、車両の旋回挙動が悪化しているときには車両の旋回挙動が安定化するよう車両の運動を制御する。尚電子制御装置30が行う車両の運動制御は操舵応答時定数係数Tpの推定値を使用して演算される目標ヨーレートγttに基づいて車両の運動を制御するものである限り、任意の制御であってよい。
 次に図2に示されたフローチャートを参照して第一の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンについて説明する。尚図2に示されたフローチャートによる制御は図には示されていないイグニッションスイッチの閉成により開始され、所定の時間毎に繰返し実行される。このことは後述の他の実施形態についても同様である。
 まずステップ10より制御が開始され、ステップ20に於いては各センサにより検出された操舵角θを示す信号等の読み込みが行われる。
 ステップ30に於いてはステップ20に於いて読み込まれた操舵角θ等に対し高周波ノイズを除去するためのローパスフィルタ処理が行われる。この場合のローパスフィルタ処理は例えば3.4Hzをカットオフ周波数とする一次のローパスフィルタ処理であってよい。
 ステップ40に於いてはスタビリティファクタKhが車両の出荷時に予め設定されているデフォルト値に設定される。尚車両の走行データに基づいてスタビリティファクタが推定される場合には、スタビリティファクタKhはその推定された値に設定されてよい。
 ステップ50に於いては車輪速度Vwiに基づいて車速Vが演算され、操舵角θに基づいて前輪の舵角δが演算されると共に、上記式11に従って規範ヨーレートγtが演算される。
 ステップ60に於いては前サイクルのステップ150に於いて更新された操舵応答時定数係数Tpの標準値Tp0、及び前サイクルのステップ160に於いて更新された基準値間隔ΔTpに基づき、下記の式12~16に従って五つの基準値Tpn(n=1~5)が設定される。尚標準値Tp0及び基準値間隔ΔTpがまだ更新されていないときには、それぞれEEPROMに記憶されているそれらのデフォルト値に設定される。
 Tp1=Tp0−2ΔTp ……(12)
 Tp2=Tp0−ΔTp ……(13)
 Tp3=Tp0 ……(14)
 Tp4=Tp0+ΔTp ……(15)
 Tp5=Tp0+2ΔTp ……(16)
 ステップ70に於いては五つの基準値Tpnについて上記式8及び11に対応する下記の式17に従って過渡ヨーレートγtrn(n=1~5)が演算される。
Figure JPOXMLDOC01-appb-I000004
 ステップ80に於いてはステップ30に於いてローパスフィルタ処理された実ヨーレートγ及びステップ70に於いて演算された過渡ヨーレートγtrnに対しセンサの零点オフセットの影響を除去するためのハイパスフィルタ処理が行われる。この場合のハイパスフィルタ処理は例えば0.2Hzをカットオフ周波数とする一次のハイパスフィルタ処理であってよい。
 上述の如くステップ30に於いてローパスフィルタ処理が行われているので、上記ハイパスフィルタ処理が行われることにより実ヨーレートγ及び過渡ヨーレートγtrnに対しバンドパスフィルタ処理が行われることと同様の結果が得られる。よってステップ80に於いてハイパスフィルタ処理された実ヨーレートγ及び過渡ヨーレートγtrnをそれぞれバンドパスフィルタ処理後の実ヨーレートγbpf及び過渡ヨーレートγtrbpfn(n=1~5)と表記する。
 ステップ90に於いては車両が旋回走行状態にあるか否かの判別が行われ、否定判別が行われたときには制御はステップ20へ戻り、肯定判別が行われたときには制御はステップ100へ進む。この場合車両が旋回走行状態にあるか否かの判別は、車両が基準値以上の車速にて走行している状況にて、車両の横加速度Gyの絶対値が基準値以上であるか否か、車両の実ヨーレートγの絶対値が基準値以上であるか否か、車両のヨーレートγと車速Vとの積の絶対値が基準値以上であるか否かの何れかの判別により行われてよい。
 ステップ100に於いてはステップ120に於いて演算されるヨーレート偏差指標値の積算値Δγan(n=1~5)を調整する必要があるか否かの判別が行われる。肯定判別が行われたときにはステップ110に於いてヨーレート偏差指標値の積算値Δγanが0にクリアされ、否定判別が行われたときには制御はステップ120へ進む。
 この場合下記の(A1)及び(A2)の何れかが成立するときに、ヨーレート偏差指標値の積算値Δγanを調整する必要があると判定されてよい。
(A1)前サイクルのステップ150に於いて操舵応答時定数係数Tpの標準値Tp0が更新された。
(A2)前サイクルのステップ160に於いて基準値間隔ΔTpが更新された。
 またステップ40に於いてスタビリティファクタKhが推定される場合には、下記の(A3)が判断条件に追加され、(A1)~(A3)の何れかが成立するときに、ヨーレート偏差指標値の積算値Δγanを調整する必要があると判定されてよい。
(A3)ヨーレート偏差指標値の積算値Δγanが前回調整されたときのスタビリティファクタKhと現サイクルのステップ40に於いて推定されたスタビリティファクタKhとの偏差ΔKhの絶対値がスタビリティファクタの偏差についての基準値以上である。
 ステップ120に於いてはバンドパスフィルタ処理後の実ヨーレートγbpfと過渡ヨーレートγtrbpfnとの偏差の大きさが前輪の舵角の偏差の大きさに置き換えられたヨーレート偏差指標値の積算値Δγan(n=1~5)が下記の式18に従って演算される。
Figure JPOXMLDOC01-appb-I000005
 尚ヨーレート偏差指標値の積算値Δγanは下記の式19に従ってバンドパスフィルタ処理後の実ヨーレートγbpfと過渡ヨーレートγtrbpfnとの偏差の大きさの積算値として演算されてもよい。
 Δγan=現在のΔγan+|γtrbpfn−γbpf| ……(19)
 ステップ130に於いては操舵応答時定数係数Tpの推定が許可される状況であるか否かの判別が行われ、否定判別が行われたときには制御はステップ20へ戻り、肯定判別が行われたときには制御はステップ140へ進む。
 この場合下記の(B1)又は(B2)が成立するときに、操舵応答時定数係数Tpの推定が許可される状況であると判定されてよい。
(B1)前回ステップ110に於いてヨーレート偏差指標値の積算値Δγanが0にクリアされた時点より基準時間以上が経過している。
(B2)五つのヨーレート偏差指標値の積算値Δγa1~Δγa5のうちの最小値が積算値の基準値以上である。
 ステップ140に於いては五つのヨーレート偏差指標値の積算値Δγa1~Δγa5のうち最も小さい値Δγam(mは1~5の何れか)が決定される。そしてヨーレート偏差指標値の積算値の最小値Δγamに対応する操舵応答時定数係数Tpmが操舵応答時定数係数Tpの推定値とされる。
 ステップ150に於いては次サイクルのステップ60に於ける五つの基準値Tpnの設定に供される標準値Tp0が、ステップ140に於いて推定された操舵応答時定数係数Tp(=Tpm)に更新され、バッファメモリに記憶される。
 ステップ160に於いてはKaを0よりも大きく1よりも小さい係数とし、基準値間隔ΔTpのガード最小値をΔTpmin(正の定数)として、基準値間隔ΔTpがKa×現在の基準値間隔ΔTp及びガード最小値ΔTpminのうちの大きい方の値に更新され、バッファメモリに記憶される。
 ステップ170に於いてはEEPROMへの操舵応答時定数係数Tpの推定値の記憶が許可される状況であるか否かの判別が行われ、否定判別が行われたときには制御はステップ20へ戻り、肯定判別が行われたときにはステップ180に於いて操舵応答時定数係数Tpの推定値がEEPROMに記憶され、これによりEEPROMに記憶されている操舵応答時定数係数Tpの推定値が更新される。
 この場合ステップ160に於いて更新された基準値間隔ΔTpとガード最小値ΔTpminとの偏差の大きさΔTpaが偏差の基準値(0以上の値)以下であるときに、EEPROMへの操舵応答時定数係数Tpの推定値の記憶が許可される状況であると判定されてよい。
 上述の如く構成された第一の実施形態の作動に於いては、ステップ50に於いて定常規範ヨーレートγtが演算され、ステップ60に於いて大きさが異なる五つの基準値Tpnが設定される。そしてステップ70に於いて定常規範ヨーレートγtに基づき五つの基準値Tpnについて過渡ヨーレートγtrnが演算される。
 またステップ120に於いてバンドパスフィルタ処理後の実ヨーレートγbpfと過渡ヨーレートγtrbpfnとの偏差の大きさが前輪の舵角の偏差の大きさに置き換えられた五つのヨーレート偏差指標値の積算値Δγa1~Δγa5が演算される。
 そしてステップ140に於いて五つのヨーレート偏差指標値の積算値Δγa1~Δγa5のうち最も小さい値Δγamが決定され、ヨーレート偏差指標値の積算値の最小値Δγamに対応する操舵応答時定数係数Tpmが操舵応答時定数係数Tpの推定値とされる。
 かくして第一の実施形態によれば、五つの基準値Tpnについて過渡ヨーレートγtrnが演算され、過渡ヨーレートγtrnと実ヨーレートγとの偏差の大きさが最も小さくなる基準値が操舵応答時定数係数Tpの推定値とされる。従って実ヨーレートγに最も近い過渡ヨーレートγtrnに対応する操舵応答時定数係数、即ち実際の操舵応答時定数係数に最も近い基準値を操舵応答時定数係数Tpの推定値とすることができる。
 また第一の実施形態によれば、ステップ150に於いて標準値Tp0がステップ140にて推定された舵応答時定数係数Tp(=Tpm)に更新される。またステップ160に於いて基準値間隔ΔTpがその大きさが小さくなるよう更新される。更にステップ60に於いて、更新された操舵応答時定数係数Tpの標準値Tp0及び基準値間隔ΔTpに基づき、大きさが異なる五つの基準値Tpnが設定される。
 従って過渡ヨーレートγtrnと実ヨーレートγとの偏差の大きさを漸次小さくして過渡ヨーレートγtrnを漸次実ヨーレートγに近づけることができる。よって実ヨーレートγに最も近い過渡ヨーレートγtrnに対応する操舵応答時定数係数、即ち操舵応答時定数係数Tpの推定値を漸次実際の操舵応答時定数係数に近づけることができる。
 また第一の実施形態によれば、ステップ30にてローパスフィルタ処理された操舵角θ等に基づいて定常規範ヨーレートγtが演算される。またステップ80に於いてステップ70にて定常規範ヨーレートγtに基づいて演算された過渡ヨーレートγtrn及び実ヨーレートγに対しハイパスフィルタ処理が行われる。
 従って検出される操舵角θ等に含まれる高周波ノイズを除去することができるだけでなく、ヨーレートセンサ36等の零点オフセットの影響を除去することができる。よってセンサの零点オフセットの影響を排除して実ヨーレートと過渡ヨーレートとの偏差を演算することができるので、実ヨーレートγ及び過渡ヨーレートγtrnに対しハイパスフィルタ処理が行われない場合に比して操舵応答時定数係数Tpを正確に推定することができる。また定常規範ヨーレートγtの演算に供される操舵角θ及び横加速度Gyに対しハイパスフィルタ処理が行われる場合に比して、ハイパスフィルタ処理の回数を低減することができ、これにより電子制御装置30の演算負荷を低減することができる。
 尚、操舵角θ等に対しローパスフィルタ処理されることなく実ヨーレートγ及び過渡ヨーレートγtrnに対しバンドパスフィルタ処理が行われてもよい。その場合には高周波ノイズを効果的に除去しつつ、上述の第一の実施形態の場合に比してフィルタ処理に要する演算の回数を低減することができ、これにより電子制御装置30の演算負荷を低減することができる。
 また第一の実施形態によれば、バンドパスフィルタ処理後の実ヨーレートγbpfと過渡ヨーレートγtrbpfnとの偏差の大きさに基づくヨーレート偏差指標値の積算値Δγanに基づいて操舵応答時定数係数Tpの推定値が求められる。従ってバンドパスフィルタ処理後の実ヨーレートγbpfと過渡ヨーレートγtrbpfnとの偏差の大きさに基づいて操舵応答時定数係数Tpの推定値が求められる場合に比して、実ヨーレートγや過渡ヨーレートγtrの瞬間的な変動に起因して操舵応答時定数係数Tpが不正確に推定される虞れを低減することができる。
 また第一の実施形態によれば、ヨーレート偏差指標値の積算値Δγanは実ヨーレートγbpfと過渡ヨーレートγtrbpfnとの偏差の大きさが前輪の舵角の偏差の大きさに置き換えられた値の積算値である。従って車速Vの影響を受けることなく操舵応答時定数係数Tpを推定することができる。よってヨーレート偏差指標値の積算値Δγanが例えば実ヨーレートγbpfと過渡ヨーレートγtrbpfnとの偏差の大きさの積算値である場合に比して、操舵応答時定数係数Tpを正確に推定することができる。また車速V毎に操舵応答時定数係数Tpを推定したり、目標ヨーレートγttの演算に供される操舵応答時定数係数Tpを車速Vによって変更したりする煩雑さを回避し、必要な演算回数や記憶手段の容量を低減することができる。
 また第一の実施形態によれば、ステップ100に於いてヨーレート偏差指標値の積算値Δγanを調整する必要があるか否かの判別が行われ、肯定判別が行われたときにはステップ110に於いてヨーレート偏差指標値の積算値Δγanが0にクリアされる。従って操舵応答時定数係数Tpの標準値Tp0や基準値間隔ΔTpが変更されたような状況に於いて、それ以前のヨーレート偏差指標値の積算値Δγanが操舵応答時定数係数Tpの推定に悪影響を及ぼすことを確実に防止することができる。
 また第一の実施形態によれば、ステップ130に於いて操舵応答時定数係数Tpの推定が許可される状況であるか否かの判別が行われ、肯定判別が行われたときにステップ140以降が実行される。従って操舵応答時定数係数Tpの推定や標準値Tp0及び基準値間隔ΔTpの更新が不必要に頻繁に行われたり、ヨーレート偏差指標値の積算値が小さい状況に於いて行われたりすることを防止することができる。
 また第一の実施形態によれば、ステップ170に於いて操舵応答時定数係数Tpの推定値の記憶が許可される状況であるか否かの判別が行われ、肯定判別が行われたときにステップ180に於いて操舵応答時定数係数Tpの推定値がEEPROMに記憶される。従って操舵応答時定数係数Tpの推定値が実質的に実際の操舵応答時定数係数に一致した段階で操舵応答時定数係数Tpの推定値をEEPROMに記憶することができる。換言すれば、操舵応答時定数係数Tpの推定値が実際の操舵応答時定数係数に実質的に一致するまで、操舵応答時定数係数Tpの推定値をEEPROMに不必要に記憶することなく操舵応答時定数係数Tpの推定を繰り返し、操舵応答時定数係数Tpの推定値を漸次実際の操舵応答時定数係数に近づけることができる。
 また第一の実施形態によれば、ステップ90に於いて車両が旋回走行状態にあるか否かの判別が行われ、肯定判別が行われたときにステップ100以降が実行される。従って車両が旋回走行状態にはなく、操舵応答時定数係数Tpの正確な推定ができない状況に於いてステップ100以降が不必要に実行されること及び操舵応答時定数係数Tpが不正確に推定されることを防止することができる。
第二の実施形態
 図3は第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第二の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの要部を示すフローチャートである。尚図3に於いて、図2に示されたステップに対応するステップには図2に於いて付されたステップ番号と同一のステップ番号が付されており、このことは後述の他の実施形態のフローチャートについても同様である。
 この第二の実施形態に於いては、ステップ70が完了すると、ステップ72に於いて単位時間当たりの運転者による往復操舵の回数が操舵周波数fsとして演算される。また操舵周波数fsが低いほどステップ80に於けるハイパスフィルタ処理のカットオフ周波数fhcが小さくなるよう、操舵周波数fsに基づき図4に示されたグラフに対応するマップよりカットオフ周波数fhcが演算される。
 そしてステップ80に於ける実ヨーレートγ及び過渡ヨーレートγtrnのハイパスフィルタ処理に於いては、カットオフ周波数がステップ72に於いて演算されたカットオフ周波数fhcに設定される。
 上述の第一の実施形態に於いては、ステップ80に於けるハイパスフィルタ処理のカットオフ周波数fhcは一定である。従ってセンサの零点オフセットの影響が確実に除去されるようカットオフ周波数fhcが高い値に設定されると、単位時間当たりの運転者による往復操舵の回数が少ない状況に於いて操舵応答時定数係数Tpを推定することができなくなる虞れがある。逆にカットオフ周波数fhcが低い値に設定されると、単位時間当たりの運転者による往復操舵の回数が多い状況に於いてセンサの零点オフセットの影響を効果的に除去することができなくなる虞れがある。
 これに対し第二の実施形態によれば、操舵周波数fsが低いほどカットオフ周波数fhcが小さくなるよう、操舵周波数fsに応じてカットオフ周波数fhcが可変設定される。従って単位時間当たりの運転者による往復操舵の回数が多い状況に於いてセンサの零点オフセットの影響を効果的に除去しつつ、単位時間当たりの運転者による往復操舵の回数が少ない状況に於いて操舵応答時定数係数Tpを推定することができなくなることを防止することができる。
 尚カットオフ周波数fhcは操舵周波数fsに基づきマップより演算されるようになっているが、操舵周波数fsの関数として演算されてもよい。
第三の実施形態
 図5は第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第三の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの要部を示すフローチャートである。
 この第三の実施形態に於いては、ステップ70が完了すると、ステップ74に於いて単位時間当たりの運転者による往復操舵の回数が操舵周波数fsとして演算される。また操舵周波数fsが低いほどハイパスフィルタ処理のカットオフ周波数fhcが小さくなると共に、車両の前後加速度Gxの絶対値が大きいほどハイパスフィルタ処理のカットオフ周波数fhcが大きくなるよう、操舵周波数fs及び車両の前後加速度Gxの絶対値に基づき図6に示されたグラフに対応するマップよりカットオフ周波数fhcが演算される。
 そしてステップ80に於ける実ヨーレートγ及び過渡ヨーレートγtrnのハイパスフィルタ処理に於いては、カットオフ周波数がステップ74に於いて演算されたカットオフ周波数fhcに設定される。
 また車両のヨーレートγ、横加速度Gy、前輪の舵角δについてセンサの零点オフセットの誤差をそれぞれγ0、Gy0、δ0とすると、車両のヨーレート、横加速度、前輪の舵角の検出値はそれぞれγ+γ0、Gy+Gy0、δ+δ0である。よってスタビリティファクタの設計値及び真の値をそれぞれKhde及びKhreとして、車両の定常旋回時に於ける定常規範ヨーレートγtと検出ヨーレートとの偏差Δγtは下記の式20により表わされる。
Figure JPOXMLDOC01-appb-I000006
 上記式20の両辺にL/Vを掛けてヨーレート偏差Δγtを前輪の舵角の偏差Δδtに換算すると、前輪の舵角の偏差Δδtは下記の式21により表わされる。
Figure JPOXMLDOC01-appb-I000007
 よってセンサの零点オフセットの影響はδ0−KhdeGy0L−γ0L/Vである。このうちδ0−KhdeGy0Lは定数であるが、γ0L/Vは車速Vに応じて変化する。従って車速Vの変化、即ち車両の前後加速度Gxの大きさが大きいほど、ヨーレート偏差指標値に与えるセンサの零点オフセットの影響が大きくなり、逆に車両の前後加速度Gxの大きさが小さいほど、ヨーレート偏差指標値に与えるセンサの零点オフセットの影響が小さくなる。
 第三の実施形態によれば、車両の前後加速度Gxの絶対値が大きいほどハイパスフィルタ処理のカットオフ周波数fhcが大きくなるよう、車両の前後加速度Gxの絶対値にも基づいてカットオフ周波数fhcが可変設定される。従って上述の第二の実施形態と同様の作用効果が得られると共に、車速Vの変化に拘らずセンサの零点オフセットの影響を効果的に除去することができる。
 尚カットオフ周波数fhcは操舵周波数fs及び車両の前後加速度Gxの絶対値に基づきマップより演算されるようになっているが、操舵周波数fs及び車両の前後加速度Gxの絶対値の関数として演算されてもよい。
第四の実施形態
 図7は第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第四の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの要部を示すフローチャートである。
 この第四の実施形態に於いては、ステップ90に於いて車両が旋回走行状態にあると判定されると、ステップ100に先立ってステップ95が実行される。ステップ95に於いては車両が高い信頼性にて操舵応答時定数係数Tpを推定し得る状況にあるか否かの判別が行われ、否定判別が行われたときには制御はステップ20へ戻り、肯定判別が行われたときには制御はステップ100へ進む。
 この場合下記の(C1)~(C3)の全てが成立するときに、車両が高い信頼性にて操舵応答時定数係数Tpを推定し得る状況にあると判定されてよい。
(C1)操舵角θの時間微分値である操舵角速度θdの絶対値が操舵角速度の基準値以上である。
(C2)走行路が悪路ではない。
(C3)制動中ではない。
 尚C1の条件は、操舵応答時定数係数Tpが車両の過渡旋回の特性を表すものであり、操舵角速度の大きさが小さい状況に於いては操舵応答時定数係数Tpを精度よく推定することができないことを考慮したものである。またC2の条件は、悪路に於いては実ヨーレートγにノイズが畳重すること、路面に対するタイヤのグリップ状態が変動し易いことを考慮したものである。更にC3の条件は、上記式11による定常規範ヨーレートγtの演算に於いては制動力の影響がないことが前提となっていることを考慮したものである。
 従って第四の実施形態によれば、車両が高い信頼性にて操舵応答時定数係数Tpを推定し得る状況にあるか否かの判別が行われない第一乃至第三の実施形態の場合に比して、操舵応答時定数係数Tpを精度よく推定することができる。
第五の実施形態
 図8及び図9はそれぞれ第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第五の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの前半部及び後半部を示すフローチャートである。
 この第五の実施形態に於いては、ステップ50が完了すると、ステップ52に於いて車両が右旋回しているか否かの判別が行われる。肯定判別が行われたときにはステップ60乃至180に於いて車両の右旋回について第一の実施形態のステップ60乃至180と同様の制御が行われる。これに対し否定判別が行われたときにはステップ65乃至185に於いて車両の左旋回について第一の実施形態のステップ60乃至180と同様の制御が行われる。
 即ちステップ60に於いては基準値Tpnに対応する右旋回時の五つの基準値Tprn(n=1~5)が設定され、ステップ70に於いては五つの基準値Tprnについて過渡ヨーレートγtrrn(n=1~5)が演算される。
 ステップ80に於いては実ヨーレートγ及び過渡ヨーレートγtrrnに対しハイパスフィルタ処理が行われる。そしてステップ100に於いてステップ120にて演算されるヨーレート偏差指標値の積算値Δγarn(n=1~5)を調整する必要があるか否かの判別が行われる。肯定判別が行われたときにはステップ110に於いてヨーレート偏差指標値の積算値Δγarnが0にクリアされ、否定判別が行われたときには制御はステップ120へ進む。
 ステップ120に於いては右旋回時のヨーレート偏差指標値の積算値Δγarn(n=1~5)が演算される。そしてステップ140に於いて五つのヨーレート偏差指標値の積算値Δγar1~Δγar5のうち最も小さい値Δγarm(mは1~5の何れか)が決定される。またヨーレート偏差指標値の積算値の最小値Δγarmに対応する操舵応答時定数係数Tprmが右旋回時の操舵応答時定数係数Tprの推定値とされる。
 ステップ150に於いては次サイクルのステップ60に於ける五つの基準値Tprnの設定に供される標準値Tpr0が、ステップ140に於いて推定された右旋回時の操舵応答時定数係数Tpr(=Tprm)に更新され、バッファメモリに記憶される。
 ステップ160に於いては基準値間隔ΔTprがKa×現在の基準値間隔ΔTpr及びガード最小値ΔTpminのうちの大きい方の値に更新され、バッファメモリに記憶される。
 更にステップ170に於いて肯定判別が行われたときには、ステップ180に於いて右旋回時の操舵応答時定数係数Tprの推定値がEEPROMに記憶され、これによりEEPROMに記憶されている右旋回時の操舵応答時定数係数Tprの推定値が更新される。
 またステップ65乃至185に於いては、右旋回を示す「r」が左旋回を示す「l」に置き換えられることにより、左旋回について上記ステップ60乃至180と同様の制御が行われる。
 右旋回時の旋回特性及び左旋回時の旋回特性が互いに異なる場合がある。例えば乗車状況若しくは荷物の積載状況に起因して重心が車両の横方向の中心にない場合や、人の昇降若しくは積載荷物の積み下ろしにより重心位置が横方向に大きく変動する場合には、車両の旋回方向によって車両の旋回特性が異なる。そのため操舵応答時定数係数も車両の旋回方向によって異なる。
 第五の実施形態によれば、上述の第一の実施形態と同様の作用効果が得られると共に、車両の旋回方向毎に操舵応答時定数係数が推定されるので、車両の旋回方向によって車両の旋回特性が異なる場合にも、車両の右旋回及び左旋回の各々について操舵応答時定数係数Tpr及びTplを精度よく推定することができる。
第六の実施形態
 図10及び図11はそれぞれ第一の実施形態の修正例として構成された本発明による旋回特性推定装置の第六の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンの前半部及び後半部を示すフローチャートである。
 この第六の実施形態に於いては、ステップ50が完了すると、ステップ54に於いて車両の横加速度Gyの絶対値が第一の基準値Gy1(正の定数)よりも大きいか否かの判別が行われる。否定判別が行われたときには制御はステップ20へ戻り、肯定判別が行われたときには制御はステップ56へ進む。
 ステップ56に於いては車両の横加速度Gyの絶対値が第二の基準値Gy2(Gy1よりも大きい正の定数)よりも大きいか否かの判別が行われる。否定判別が行われたときにはステップ60乃至180に於いて車両の横加速度Gyの絶対値が第一の基準値Gy1よりも大きく第二の基準値Gy2以下である場合(横加速度Gyの第一の領域)について第一の実施形態のステップ60乃至180と同様の制御が行われる。これに対し肯定判別が行われたときにはステップ65乃至185に於いて車両の横加速度Gyの絶対値が第二の基準値Gy2よりも大きい場合(横加速度Gyの第二の領域)について第一の実施形態のステップ60乃至180と同様の制御が行われる。
 即ちステップ60に於いては横加速度Gyの第一の領域について基準値Tpnに対応する五つの基準値Tp1n(n=1~5)が設定され、ステップ70に於いては五つの基準値Tp1nについて過渡ヨーレートγtr1n(n=1~5)が演算される。
 ステップ80に於いては実ヨーレートγ及び過渡ヨーレートγtr1nに対しハイパスフィルタ処理が行われる。そしてステップ100に於いてステップ120にて演算されるヨーレート偏差指標値の積算値Δγa1n(n=1~5)を調整する必要があるか否かの判別が行われる。肯定判別が行われたときにはステップ110に於いてヨーレート偏差指標値の積算値Δγa1nが0にクリアされ、否定判別が行われたときには制御はステップ120へ進む。
 ステップ120に於いては横加速度Gyの第一の領域についてヨーレート偏差指標値の積算値Δγa1n(n=1~5)が演算される。そしてステップ140に於いて五つのヨーレート偏差指標値の積算値Δγa11~Δγa15のうち最も小さい値Δγa1m(mは1~5の何れか)が決定される。またヨーレート偏差指標値の積算値の最小値Δγa1mに対応する操舵応答時定数係数Tp1mが横加速度Gyの第一の領域ついての操舵応答時定数係数Tp1の推定値とされる。
 ステップ150に於いては次サイクルのステップ60に於ける五つの基準値Tp1nの設定に供される標準値Tp10が、ステップ140に於いて推定された舵応答時定数係数Tp1(=Tp1m)に更新され、バッファメモリに記憶される。
 ステップ160に於いては基準値間隔ΔTp1がKa×現在の基準値間隔ΔTp1及びガード最小値ΔTpminのうちの大きい方の値に更新され、バッファメモリに記憶される。
 更にステップ170に於いて肯定判別が行われたときには、ステップ180に於いて操舵応答時定数係数Tp1の推定値がEEPROMに記憶され、これによりEEPROMに記憶されている横加速度Gyの第一の領域ついての操舵応答時定数係数Tp1の推定値が更新される。
 またステップ65乃至185に於いては、横加速度Gyの第一の領域を示す「1」が横加速度Gyの第二の領域を示す「2」に置き換えられることにより、横加速度Gyの第二の領域について上記ステップ60乃至180と同様の制御が行われる。
 一般に車両の旋回特性は車両の横加速度Gyの大きさによって異なる。第六の実施形態によれば、車両の横加速度Gyの大きさの領域毎に操舵応答時定数係数Tpが推定されるので、上述の第一の実施形態と同様の作用効果が得られると共に、車両の横加速度Gyの大きさが大きく変動する場合にも、横加速度Gyの大きさの各領域について操舵応答時定数係数Tpを精度よく推定することができる。
第七の実施形態
 図12は車両の運動制御装置に適用された本発明による旋回特性推定装置の第七の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンを示すフローチャートである。
 この第七の実施形態に於いては、上述の第一の実施形態のステップ60及びステップ100~180は実行されない。従ってステップ70に於いては上記式17に対応する下記の式22に従って過渡ヨーレートγtrが演算される。またステップ90に於いて肯定判別が行われると、制御はステップ200へ進み、図13に示されたフローチャートに従って操舵応答時定数係数Tpが推定される。
Figure JPOXMLDOC01-appb-I000008
 図13に示されている如く、ステップ210に於いてはバンドパスフィルタ処理後の過渡ヨーレートγtrbpf及び実ヨーレートγbpfの位相差Dyの絶対値が基準値Dy0(正の定数)より大きいか否かの判別が行われる。否定判別が行われたときには操舵応答時定数係数Tpの推定値の増減修正は不要であるので、制御はステップ20へ戻り、肯定判別が行われたときには制御はステップ220へ進む。
 ステップ220に於いては過渡ヨーレートγtrbpfの位相が実ヨーレートγbpfの位相より早いか否かの判別が行われ、否定判別が行われたときには制御はステップ240へ進み、肯定判別が行われたときには制御はステップ230へ進む。
 ステップ230に於いては操舵応答時定数係数Tpの前回の推定値Tpfと前回のステップ260に於いて設定された増減修正量ΔTpbとの和Tpf+ΔTpb及び予め設定された操舵応答時定数係数の最大値Tpmaxのうち小さい方の値が操舵応答時定数係数Tpの推定値とされる。
 ステップ240に於いては操舵応答時定数係数Tpの前回の推定値Tpfより増減修正量ΔTpbを減算した値Tpf−ΔTpb及び予め設定された操舵応答時定数係数の最小値Tpmin(Tpmaxよりも小さい正の定数)のうち大きい方の値が操舵応答時定数係数Tpの推定値とされる。
 ステップ250に於いては今回の操舵応答時定数係数Tpの推定値と前回の推定値Tpfとの偏差の絶対値が記憶判定の基準値Tpe(正の定数)未満であるか否かの判別が行われ、肯定判別が行われたときには制御はステップ270へ進み、否定判別が行われたときには制御はステップ260へ進む。
 ステップ260に於いては次サイクルの制御のための「前回の推定値Tpf」が操舵応答時定数係数Tpの推定値に書き換えられると共に、Kbを例えば0.5の如く0よりも大きく1よりも小さい係数として、増減修正量ΔTpbが係数Kbと現在の増減修正量ΔTpbとの積に書き換えられる。
 ステップ270に於いては舵応答時定数係数Tpの推定値がEEPROMに記憶され、ステップ280に於いては前回の推定値Tpfが0にクリアされると共に、増減修正量ΔTpbが0にクリアされる。
 かくして第七の実施形態によれば、過渡ヨーレートγtrbpfの位相が実ヨーレートγbpfの位相より早いときには、ステップ220に於いて肯定判別が行われ、ステップ230に於いて操舵応答時定数係数Tpの推定値が増大補正される。これに対し過渡ヨーレートγtrbpfの位相が実ヨーレートγbpfの位相より遅いときにはステップ220に於いて否定判別が行われ、ステップ240に於いて操舵応答時定数係数Tpの推定値が低減補正される。
 従って第七の実施形態によれば、車両の過渡ヨーレートγtrと車両の実ヨーレートγとの位相差が小さくなることによって操舵応答時定数係数Tpの推定値が真の操舵応答時定数係数に近づくよう操舵応答時定数係数Tpの推定値を修正することができる。よって操舵応答時定数係数Tpを真の操舵応答時定数係数に近い値に推定することができる。
 また第七の実施形態によれば、ステップ50に於いてステップ30にてローパスフィルタ処理された操舵角θ等に基づいて定常規範ヨーレートγtが演算される。またステップ80に於いてステップ70にて定常規範ヨーレートγtに基づいて演算された過渡ヨーレートγtr及び実ヨーレートγに対しハイパスフィルタ処理が行われる。
 従って上述の第一乃至第六の実施形態の場合と同様、検出される操舵角θ等に含まれる高周波ノイズを除去することができるだけでなく、ヨーレートセンサ36等の零点オフセットの影響を除去することができる。よってセンサの零点オフセットの影響を排除して実ヨーレートと過渡ヨーレートとの偏差を演算することができるので、実ヨーレートγ及び過渡ヨーレートγtrに対しハイパスフィルタ処理が行われない場合に比して操舵応答時定数係数Tpを正確に推定することができる。
 また第七の実施形態によれば、ステップ90に於いて車両が旋回走行状態にあるか否かの判別が行われ、肯定判別が行われたときにステップ200以降が実行される。従って上述の第一乃至第六の実施形態の場合と同様、車両が旋回走行状態にはなく、操舵応答時定数係数Tpの正確な推定ができない状況に於いてステップ200以降が不必要に実行されること及び操舵応答時定数係数Tpが不正確に推定されることを防止することができる。
 特に第七の実施形態によれば、操舵応答時定数係数Tpと前回の推定値Tpfとの偏差の絶対値が記憶判定の基準値Tpe以上であるときには、ステップ250に於いて否定判別が行われる。そしてステップ260に於いて増減修正量ΔTpbが係数Kbと現在の増減修正量ΔTpbとの積に書き換えられることにより低減される。
 従って増減修正量ΔTpbが低減されない場合に比して、過渡ヨーレートγtrと実ヨーレートγとの位相差を効率的に小さくすることができ、これにより操舵応答時定数係数Tpの推定値を効率的に真の値に近づけることができる。
 また第七の実施形態によれば、操舵応答時定数係数Tpと前回の推定値Tpfとの偏差の絶対値が記憶判定の基準値Tpe未満になっときにステップ250に於いて肯定判別が行われ、ステップ270に於いて舵応答時定数係数Tpの推定値がEEPROMに記憶され更新される。
 従ってステップ250の判別が行われない場合に比して、操舵応答時定数係数Tpの推定値が真の操舵応答時定数係数に近い値にならない段階で舵応答時定数係数Tpの推定値がEEPROMに記憶され更新されることを防止することができ、また舵応答時定数係数Tpの推定値の記憶及び更新の頻度を低減することができる。
 更に第七の実施形態によれば、バンドパスフィルタ処理後の過渡ヨーレートγtrbpf及び実ヨーレートγbpfの位相差Dyの絶対値が基準値Dy0より小さいときには、ステップ210に於いて否定判別が行われ、ステップ220以降のステップは実行されない。
 従って操舵応答時定数係数Tpの推定値が真の操舵応答時定数係数に近い値である状況に於いて、操舵応答時定数係数Tpの推定や更新が不必要に繰り返し行われることを防止することができる。
第八の実施形態
 図14は第七の実施形態の修正例として構成された本発明による旋回特性推定装置の第八の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンを示すフローチャートである。
 この第八の実施形態に於いては、ステップ50が完了すると、上述の第五の実施形態の場合と同様にステップ52が実行される。そしてステップ52に於いて肯定判別が行われたときには制御はステップ70へ進み、否定判別が行われたときには制御はステップ75へ進む。
 ステップ70に於いては後述の右旋回時の操舵応答時定数係数Tprに基づき上記式17に対応する下記の式23に従って右旋回時の過渡ヨーレートγtrrが演算される。
Figure JPOXMLDOC01-appb-I000009
 ステップ80に於いては実ヨーレートγ及び過渡ヨーレートγtrrに対しハイパスフィルタ処理が行われ、これによりハイパスフィルタ処理後の実ヨーレートγbpf及び過渡ヨーレートγtrrbpfが演算される。
 またステップ90に於いて肯定判別が行われると、制御はステップ200へ進み、図15に示されたフローチャートに従って右旋回時の操舵応答時定数係数Tprが推定され、舵応答時定数係数Tprの推定値がEEPROMに記憶される。
 同様にステップ75に於いては後述の左旋回時の操舵応答時定数係数Tplに基づき上記式17に対応する下記の式24に従って左旋回時の過渡ヨーレートγtrlが演算される。
Figure JPOXMLDOC01-appb-I000010
 ステップ85に於いては実ヨーレートγ及び過渡ヨーレートγtrlに対しハイパスフィルタ処理が行われ、これによりハイパスフィルタ処理後の実ヨーレートγbpf及び過渡ヨーレートγtrlbpfが演算される。
 またステップ95がステップ90と同様に実行される。ステップ95に於いて肯定判別が行われると、制御はステップ205へ進み、図16に示されたフローチャートに従って左旋回時の操舵応答時定数係数Tplが推定され、舵応答時定数係数Tplの推定値がEEPROMに記憶される。
 図15に示されている如く、ステップ210に於いてはバンドパスフィルタ処理後の右旋回時の過渡ヨーレートγtrrbpf及び実ヨーレートγbpfの位相差Dyrの絶対値が基準値Dy0より大きいか否かの判別が行われる。否定判別が行われたときには操舵応答時定数係数Tprの推定値の増減修正は不要であるので、制御はステップ20へ戻り、肯定判別が行われたときには制御はステップ220へ進む。
 ステップ220に於いては過渡ヨーレートγtrrbpfの位相が実ヨーレートγbpfの位相より早いか否かの判別が行われ、否定判別が行われたときには制御はステップ240へ進み、肯定判別が行われたときには制御はステップ230へ進む。
 ステップ230に於いては右旋回時の操舵応答時定数係数Tprの前回の推定値Tprfと前回のステップ260に於いて設定された右旋回時の増減修正量ΔTprbとの和Tprf+ΔTprb及び予め設定された操舵応答時定数係数の最大値Tpmaxのうち小さい方の値が右旋回時の操舵応答時定数係数Tprの推定値とされる。
 ステップ240に於いては右旋回時の操舵応答時定数係数Tprの前回の推定値Tprfより増減修正量ΔTprbを減算した値Tprf−ΔTprb及び予め設定された操舵応答時定数係数の最小値Tpminのうち大きい方の値が右旋回時の操舵応答時定数係数Tprの推定値とされる。
 ステップ250に於いては操舵応答時定数係数Tprの推定値と前回の推定値Tprfとの偏差の絶対値が記憶判定の基準値Tpe未満であるか否かの判別が行われ、肯定判別が行われたときには制御はステップ270へ進み、否定判別が行われたときには制御はステップ260へ進む。
 ステップ260に於いては次サイクルの制御のための「前回の推定値Tprf」が操舵応答時定数係数Tprの推定値に書き換えられると共に、増減修正量ΔTprbが係数Kbと現在の増減修正量ΔTprbとの積に書き換えられる。
 ステップ270に於いては舵応答時定数係数Tprの推定値がEEPROMに記憶され、ステップ280に於いては前回の推定値Tprfが0にクリアされると共に、増減修正量ΔTprbが0にクリアされる。
 またステップ215乃至285に於いては、右旋回を示す「r」が左旋回を示す「l」に置き換えられることにより、左旋回について上記ステップ210乃至280と同様の制御が行われる。
 従って第八の実施形態によれば、上述の第七の実施形態と同様の作用効果が得られると共に、車両の旋回方向によって車両の旋回特性が異なる場合にも、車両の右旋回及び左旋回の各々について操舵応答時定数係数Tpr及びTplを精度よく推定することができる。
第九の実施形態
 図17は第七の実施形態の修正例として構成された本発明による旋回特性推定装置の第九の実施形態に於ける操舵応答時定数係数Tpの推定演算ルーチンを示すフローチャートである。
 この第九の実施形態に於いては、ステップ50が完了すると、上述の第六の実施形態の場合と同様にステップ54及び56が実行される。そしてステップ56に於いて否定判別が行われたときには制御はステップ70へ進み、否定判別が行われたときには制御はステップ75へ進む。
 ステップ70に於いては上記式17に対応する下記の式25に従って横加速度Gyの第一の領域についての過渡ヨーレートγtr1が演算される。
Figure JPOXMLDOC01-appb-I000011
 ステップ80に於いては実ヨーレートγ及び過渡ヨーレートγtr1に対しハイパスフィルタ処理が行われ、これによりハイパスフィルタ処理後の実ヨーレートγbpf及び過渡ヨーレートγtr1bpfが演算される。
 またステップ90に於いて肯定判別が行われると、制御はステップ200へ進み、図18に示されたフローチャートに従って横加速度Gyの第一の領域についての操舵応答時定数係数Tp1が推定され、舵応答時定数係数Tp1の推定値がEEPROMに記憶される。
 同様にステップ75に於いては上記式17に対応する下記の式26に従って横加速度Gyの第二の領域についての過渡ヨーレートγtr2が演算される。
Figure JPOXMLDOC01-appb-I000012
 ステップ85に於いては実ヨーレートγ及び過渡ヨーレートγtr2に対しハイパスフィルタ処理が行われ、これによりハイパスフィルタ処理後の実ヨーレートγbpf及び過渡ヨーレートγtr2bpfが演算される。
 またステップ95がステップ90と同様に実行される。ステップ95に於いて肯定判別が行われると、制御はステップ205へ進み、図19に示されたフローチャートに従って横加速度Gyの第二の領域についての操舵応答時定数係数Tp2が推定され、舵応答時定数係数Tp2の推定値がEEPROMに記憶される。
 図18に示されている如く、ステップ210に於いては横加速度Gyの第一の領域についてのバンドパスフィルタ処理後の過渡ヨーレートγtr1bpf及び実ヨーレートγbpfの位相差Dy1の絶対値が基準値Dy0より大きいか否かの判別が行われる。否定判別が行われたときには操舵応答時定数係数Tp1の推定値の増減修正は不要であるので、制御はステップ20へ戻り、肯定判別が行われたときには制御はステップ220へ進む。
 ステップ220に於いては過渡ヨーレートγtr1bpfの位相が実ヨーレートγbpfの位相より早いか否かの判別が行われ、否定判別が行われたときには制御はステップ240へ進み、肯定判別が行われたときには制御はステップ230へ進む。
 ステップ230に於いては横加速度Gyの第一の領域についての操舵応答時定数係数Tp1の前回の推定値Tp1fと前回のステップ260に於いて設定された横加速度Gyの第一の領域についての増減修正量ΔTp1bとの和Tp1f+ΔTp1b及び予め設定された操舵応答時定数係数の最大値Tpmaxのうち小さい方の値が横加速度Gyの第一の領域についての操舵応答時定数係数Tp1の推定値とされる。
 ステップ240に於いては横加速度Gyの第一の領域についての操舵応答時定数係数Tp1の前回の推定値Tp1fより増減修正量ΔTp1bを減算した値Tp1f−ΔTp1b及び予め設定された操舵応答時定数係数の最小値Tpminのうち大きい方の値が横加速度Gyの第一の領域についての操舵応答時定数係数Tp1の推定値とされる。
 ステップ250に於いては操舵応答時定数係数Tp1の推定値と前回の推定値Tp1fとの偏差の絶対値が記憶判定の基準値Tpe未満であるか否かの判別が行われ、肯定判別が行われたときには制御はステップ270へ進み、否定判別が行われたときには制御はステップ260へ進む。
 ステップ260に於いては次サイクルの制御のための「前回の推定値Tp1f」が操舵応答時定数係数Tp1の推定値に書き換えられると共に、増減修正量ΔTp1bが係数Kbと現在の増減修正量ΔTp1bとの積に書き換えられる。
 ステップ270に於いては舵応答時定数係数Tp1の推定値がEEPROMに記憶され、ステップ280に於いては前回の推定値Tp1fが0にクリアされると共に、増減修正量ΔTp1bが0にクリアされる。
 またステップ215乃至285に於いては、横加速度Gyの第一の領域を示す「1」が横加速度Gyの第二の領域を示す「2」に置き換えられることにより、横加速度Gyの第二の領域について上記ステップ210乃至280と同様の制御が行われる。
 従って第九の実施形態によれば、上述の第七の実施形態と同様の作用効果が得られると共に、車両の横加速度Gyの大きさが大きく変動する場合にも、横加速度Gyの大きさの各領域について操舵応答時定数係数Tpを精度よく推定することができる。
 尚上述の第七乃至第九の実施形態に於いては、他の実施形態と同様に各サイクル毎に過渡ヨーレートγtrbpf等が演算されるようになっているが、各サイクル毎に読み込まれた操舵角θ等がバッファメモリに記憶され、旋回終了後に各サイクルの操舵角θ等に基づいて過渡ヨーレートγtrbpf等が演算され、それらの過渡ヨーレートγtrbpf等について第七乃至第九の実施形態のステップ200又は205の処理が行われるよう修正されてもよい。
 以上に於いては本発明を特定の実施形態について詳細に説明したが、本発明は上述の実施形態に限定されるものではなく、本発明の範囲内にて他の種々の実施形態が可能であることは当業者にとって明らかであろう。
 例えば上述の第一乃至第六の実施形態に於いては、ステップ60に於いて操舵応答時定数係数Tpの標準値Tp0及び基準値間隔ΔTpに基づき、五つの基準値Tpn(n=1~5)が設定されるようになっているが、基準値は例えば式13乃至15に従って三つ演算されるよう修正されてもよく、また七つ以上の基準値が演算されるよう修正されてもよい。
 また上述の第一乃至第六の実施形態に於いては、基準値間隔ΔTp等が漸次小さくされるようになっているが、基準値間隔ΔTp等は一定の値であってもよい。
 また上述の第一乃至第六の実施形態に於いては、ヨーレート偏差指標値の積算値Δγan等を調整する必要があると判定されると、ヨーレート偏差指標値の積算値Δγan等が0にクリアされる。しかしヨーレート偏差指標値の積算値Δγan等は例えば0に近い正の値の如く、現在の値よりも小さい正の値に調整されるよう修正されてもよい。
 また上述の第六及び第九の実施形態に於いては、車両の横加速度Gyの大きさが異なる第一及び第二の領域について操舵応答時定数係数Tpが推定されるようになっているが、車両の横加速度Gyの大きさが異なる三つ以上の領域について操舵応答時定数係数Tpが推定されるよう修正されてもよい。
 また上述の第七乃至第九の実施形態に於いては、増減修正量ΔTpb等が漸次小さくされるようになっているが、増減修正量ΔTpb等は一定の値であってもよい。
 また上述の第五及び第六の実施形態は上述の第一の実施形態の修正例として構成されているが、第五又は第六の実施形態の構成が第二乃至第四の実施形態の何れかに適用されてもよい。同様に上述の第七の実施形態は上述の第一の実施形態の修正例として構成されているが、第七の実施形態の構成が第二乃至第四の実施形態の何れかに適用されてもよい。

Claims (12)

  1. 車両の規範ヨーレートに対し一次遅れの関係にある車両の過渡ヨーレートを車両の実ヨーレートに対応する値とみなし、一次遅れの時定数の車速にかかる係数である操舵応答時定数係数を推定する車両の旋回特性推定装置であって、操舵応答時定数係数の推定値を使用して車両の規範ヨーレートに基づいて車両の過渡ヨーレートを演算し、車両の過渡ヨーレートが車両の実ヨーレートに近づくよう車両の過渡ヨーレートと車両の実ヨーレートとの関係に基づいて操舵応答時定数係数の推定値を修正する車両の旋回特性推定装置。
  2. 車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさ及び車両の過渡ヨーレートと車両の実ヨーレートとの位相差の少なくとも一方が小さくなることによって車両の過渡ヨーレートが車両の実ヨーレートに近づくよう操舵応答時定数係数の推定値を修正することを特徴とする請求項1に記載の車両の旋回特性推定装置。
  3. 記憶されている操舵応答時定数係数を第一の基準値として、前記第一の基準値、前記第一の基準値よりも大きい第二の基準値、前記第一の基準値よりも小さい第三の基準値を使用して三つの車両の過渡ヨーレートを演算し、前記第一乃至第三の基準値のうち、車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさの指標値であるヨーレート偏差指標値が最も小さくなる基準値を操舵応答時定数係数の推定値とすることを特徴とする請求項2に記載の車両の旋回特性推定装置。
  4. 前記操舵応答時定数係数の推定値を新たな第一の基準値とし、前記新たな第一の基準値、前記新たな第一の基準値よりも大きい新たな第二の基準値、前記新たな第一の基準値よりも小さい新たな第三の基準値を使用して三つの車両の過渡ヨーレートを演算し、前記新たな第一乃至第三の基準値のうち、前記ヨーレート偏差指標値が最も小さくなる基準値を新たな操舵応答時定数係数の推定値とすることを特徴とする請求項3に記載の車両の旋回特性推定装置。
  5. 前記新たな第二の基準値と前記新たな第一の基準値との差は前記第二の基準値と前記第一の基準値との差よりも小さく、前記新たな第三の基準値と前記新たな第一の基準値との差は前記第三の基準値と前記第一の基準値との差よりも小さいことを特徴とする請求項4に記載の車両の旋回特性推定装置。
  6. 前記ヨーレート偏差指標値は、第一の所定周波数以下の成分が除去された車両の実ヨーレートと、第二の所定周波数以下の成分が除去された車両の過渡ヨーレートとの偏差の大きさに基づく値であることを特徴とする請求項3に記載の車両の旋回特性推定装置。
  7. 単位時間当たりの運転者による往復操舵の多さの指標値に応じて前記第一の所定周波数若しくは前記第二の所定周波数を変更することを特徴とする請求項6に記載の車両の旋回特性推定装置。
  8. 車両の前後加速度の大きさに応じて前記第一の所定周波数若しくは前記第二の所定周波数を変更することを特徴とする請求項6に記載の車両の旋回特性推定装置。
  9. 前記ヨーレート偏差指標値は車両の過渡ヨーレートと車両の実ヨーレートとの偏差の大きさを前輪の舵角の偏差の大きさに換算した値であることを特徴とする請求項3に記載の車両の旋回特性推定装置。
  10. 車両の過渡ヨーレートの位相が車両の実ヨーレートの位相よりも早いときには操舵応答時定数係数の推定値を増大修正し、車両の過渡ヨーレートの位相が車両の実ヨーレートの位相よりも遅いときには操舵応答時定数係数の推定値を低減修正することを特徴とする請求項2に記載の車両の旋回特性推定装置。
  11. 車両の右旋回及び左旋回について個別に操舵応答時定数係数を推定することを特徴とする請求項1乃至10の何れか一つに記載の車両の旋回特性推定装置。
  12. 車両の横加速度の大きさの領域毎に操舵応答時定数係数を推定することを特徴とする請求項1乃至10の何れか一つに記載の車両の旋回特性推定装置。
PCT/JP2009/067128 2009-09-24 2009-09-24 車両の旋回特性推定装置 WO2011036821A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CN2009801616535A CN102574527B (zh) 2009-09-24 2009-09-24 车辆的转弯特性估计装置
PCT/JP2009/067128 WO2011036821A1 (ja) 2009-09-24 2009-09-24 車両の旋回特性推定装置
DE112009005281.9T DE112009005281B4 (de) 2009-09-24 2009-09-24 Verfahren und Vorrichtung zum Schätzen einer Drehcharakteristik eines Fahrzeugs
US13/389,696 US8744686B2 (en) 2009-09-24 2009-09-24 Device for estimating turning characteristic of vehicle
JP2010541617A JP5041069B2 (ja) 2009-09-24 2009-09-24 車両の旋回特性推定装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2009/067128 WO2011036821A1 (ja) 2009-09-24 2009-09-24 車両の旋回特性推定装置

Publications (1)

Publication Number Publication Date
WO2011036821A1 true WO2011036821A1 (ja) 2011-03-31

Family

ID=43795591

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2009/067128 WO2011036821A1 (ja) 2009-09-24 2009-09-24 車両の旋回特性推定装置

Country Status (5)

Country Link
US (1) US8744686B2 (ja)
JP (1) JP5041069B2 (ja)
CN (1) CN102574527B (ja)
DE (1) DE112009005281B4 (ja)
WO (1) WO2011036821A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014227141A (ja) * 2013-05-27 2014-12-08 トヨタ自動車株式会社 車両挙動制御装置および車両状態推定装置
JPWO2013125047A1 (ja) * 2012-02-24 2015-07-30 トヨタ自動車株式会社 車両の旋回特性推定装置
US9261032B2 (en) 2011-01-18 2016-02-16 Toyota Jidosha Kabushiki Kaisha Air-fuel ratio control apparatus for an internal combustion engine

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5024456B2 (ja) * 2009-09-24 2012-09-12 トヨタ自動車株式会社 車両の旋回特性推定装置
ES2676170T3 (es) * 2009-12-25 2018-07-17 Yamaha Hatsudoki Kabushiki Kaisha Dispositivo de evaluación de características de conductor y vehículo de montar a horcajadas provisto del mismo
CN103303367B (zh) * 2013-06-21 2015-06-24 电子科技大学 一种四轮驱动电动车车身稳定控制方法
US9205844B2 (en) * 2014-03-05 2015-12-08 Ford Global Technologies, Llc Driver drowsiness detection
GB2536008B (en) 2015-03-03 2019-06-12 Jaguar Land Rover Ltd Vehicle state estimation apparatus and method
KR101868843B1 (ko) * 2016-11-29 2018-06-19 주식회사 만도 차량 제어 장치 및 그 제어 방법
CN109515444B (zh) * 2017-09-18 2021-02-26 百度在线网络技术(北京)有限公司 用于输出无人驾驶汽车的行驶性能指标的方法和装置
CN108749919B (zh) * 2018-04-17 2021-03-09 南京航空航天大学 一种线控四轮独立转向系统容错控制系统及其控制方法
CN109765888B (zh) * 2018-12-27 2021-12-31 东软睿驰汽车技术(沈阳)有限公司 一种车辆控制方法、装置及系统
JP7415975B2 (ja) * 2021-02-12 2024-01-17 トヨタ自動車株式会社 車両の走行支援システム、及び、車両の走行支援方法
US11724739B2 (en) * 2021-07-22 2023-08-15 GM Global Technology Operations LLC Vehicle actuation commands to affect transient handling

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10258720A (ja) * 1997-03-19 1998-09-29 Mitsubishi Motors Corp 車両の旋回制御装置
JP2000095085A (ja) * 1998-09-22 2000-04-04 Toyota Motor Corp 車輌の基準ヨーレート演算装置
JP2004026074A (ja) * 2002-06-27 2004-01-29 Toyota Motor Corp 車輌の旋回特性推定装置
JP2004026073A (ja) * 2002-06-27 2004-01-29 Toyota Motor Corp 車輌の旋回特性推定装置
JP2005008067A (ja) * 2003-06-19 2005-01-13 Toyota Motor Corp 車輌の旋回特性推定装置
JP2006027388A (ja) * 2004-07-14 2006-02-02 Advics:Kk 車両の姿勢制御装置
JP2006168630A (ja) * 2004-12-17 2006-06-29 Toyota Motor Corp 旋回制御装置
JP2009119921A (ja) * 2007-11-12 2009-06-04 Mitsubishi Electric Corp 車両用操舵装置及び車両制御装置
JP2009137512A (ja) * 2007-12-10 2009-06-25 Honda Motor Co Ltd 車両の後輪舵角制御装置および後輪舵角制御方法
JP2009274507A (ja) * 2008-05-13 2009-11-26 Mitsubishi Electric Corp 車両の旋回特性推定装置及び旋回特性推定方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4706979A (en) * 1985-07-12 1987-11-17 Nissan Motor Co., Ltd. Steering control system for wheeled vehicle
JP3357159B2 (ja) * 1993-08-10 2002-12-16 三菱自動車工業株式会社 車両運転操作状態の推定方法および車両運転特性制御方法
JP3116738B2 (ja) 1994-07-28 2000-12-11 トヨタ自動車株式会社 車輌の挙動制御装置
JP4069921B2 (ja) * 2004-10-25 2008-04-02 三菱自動車工業株式会社 車両の旋回挙動制御装置
WO2010082288A1 (ja) * 2009-01-13 2010-07-22 トヨタ自動車株式会社 車両状態推定装置

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10258720A (ja) * 1997-03-19 1998-09-29 Mitsubishi Motors Corp 車両の旋回制御装置
JP2000095085A (ja) * 1998-09-22 2000-04-04 Toyota Motor Corp 車輌の基準ヨーレート演算装置
JP2004026074A (ja) * 2002-06-27 2004-01-29 Toyota Motor Corp 車輌の旋回特性推定装置
JP2004026073A (ja) * 2002-06-27 2004-01-29 Toyota Motor Corp 車輌の旋回特性推定装置
JP2005008067A (ja) * 2003-06-19 2005-01-13 Toyota Motor Corp 車輌の旋回特性推定装置
JP2006027388A (ja) * 2004-07-14 2006-02-02 Advics:Kk 車両の姿勢制御装置
JP2006168630A (ja) * 2004-12-17 2006-06-29 Toyota Motor Corp 旋回制御装置
JP2009119921A (ja) * 2007-11-12 2009-06-04 Mitsubishi Electric Corp 車両用操舵装置及び車両制御装置
JP2009137512A (ja) * 2007-12-10 2009-06-25 Honda Motor Co Ltd 車両の後輪舵角制御装置および後輪舵角制御方法
JP2009274507A (ja) * 2008-05-13 2009-11-26 Mitsubishi Electric Corp 車両の旋回特性推定装置及び旋回特性推定方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9261032B2 (en) 2011-01-18 2016-02-16 Toyota Jidosha Kabushiki Kaisha Air-fuel ratio control apparatus for an internal combustion engine
JPWO2013125047A1 (ja) * 2012-02-24 2015-07-30 トヨタ自動車株式会社 車両の旋回特性推定装置
US9604636B2 (en) 2012-02-24 2017-03-28 Toyota Jidosha Kabushiki Kaisha Turning characteristic estimating device for vehicle
JP2014227141A (ja) * 2013-05-27 2014-12-08 トヨタ自動車株式会社 車両挙動制御装置および車両状態推定装置

Also Published As

Publication number Publication date
DE112009005281T5 (de) 2012-12-27
US8744686B2 (en) 2014-06-03
US20120173039A1 (en) 2012-07-05
JP5041069B2 (ja) 2012-10-03
DE112009005281B4 (de) 2019-07-11
JPWO2011036821A1 (ja) 2013-02-14
CN102574527A (zh) 2012-07-11
CN102574527B (zh) 2013-08-28

Similar Documents

Publication Publication Date Title
JP5041069B2 (ja) 車両の旋回特性推定装置
JP5024456B2 (ja) 車両の旋回特性推定装置
JP5858051B2 (ja) 車両の積載状態推定方法及び装置
JP3060923B2 (ja) 車両状態推定装置
JP6056954B2 (ja) 車両の走行運動制御装置
JP4151389B2 (ja) 車輌の挙動制御装置
JP3829934B2 (ja) 車輌の旋回特性推定装置
JPH0840232A (ja) 車輌の挙動制御装置
JP2000052963A (ja) 車輌の挙動制御方法
JPH11101732A (ja) 路面摩擦係数検出装置
JP2005014901A (ja) アクティブ法線力調整システムを備えたビークルダイナミクスコントロールシステムおよびその調整方法
JPH09118212A (ja) 車体の横滑り速度推定装置
KR20180068681A (ko) 차량의 거동 산출 시스템 및 방법
JP3282449B2 (ja) 車輌の横滑り状態量検出装置
JP2004026074A (ja) 車輌の旋回特性推定装置
WO2013088580A1 (ja) 車両のスタビリティファクタ推定装置
JP4228792B2 (ja) 車輌の旋回特性推定装置
JP2002173012A (ja) 車輌の挙動制御装置
JP2013129316A (ja) 車両の運動制御装置
JP2013132938A (ja) 車両のスタビリティファクタ推定装置
JP2013129351A (ja) 車両のスタビリティファクタ推定装置
JP2005206075A (ja) 車輌の制御装置
JP2005008067A (ja) 車輌の旋回特性推定装置
JP2001225733A (ja) 車輌の走行制御装置

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200980161653.5

Country of ref document: CN

WWE Wipo information: entry into national phase

Ref document number: 2010541617

Country of ref document: JP

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09849851

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 13389696

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 1120090052819

Country of ref document: DE

Ref document number: 112009005281

Country of ref document: DE

122 Ep: pct application non-entry in european phase

Ref document number: 09849851

Country of ref document: EP

Kind code of ref document: A1