JP2007001330A - Braking and driving force distribution device of four wheel independent driving vehicle - Google Patents
Braking and driving force distribution device of four wheel independent driving vehicle Download PDFInfo
- Publication number
- JP2007001330A JP2007001330A JP2005180371A JP2005180371A JP2007001330A JP 2007001330 A JP2007001330 A JP 2007001330A JP 2005180371 A JP2005180371 A JP 2005180371A JP 2005180371 A JP2005180371 A JP 2005180371A JP 2007001330 A JP2007001330 A JP 2007001330A
- Authority
- JP
- Japan
- Prior art keywords
- braking
- driving force
- wheel
- force distribution
- δfx
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 319
- 230000008859 change Effects 0.000 claims abstract description 78
- 238000012937 correction Methods 0.000 claims description 152
- 230000035945 sensitivity Effects 0.000 claims description 29
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000001052 transient effect Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims 3
- 238000001514 detection method Methods 0.000 claims 1
- 238000013021 overheating Methods 0.000 abstract description 21
- 238000000034 method Methods 0.000 description 74
- 230000006399 behavior Effects 0.000 description 73
- 241001125929 Trisopterus luscus Species 0.000 description 70
- 239000013589 supplement Substances 0.000 description 65
- 230000006870 function Effects 0.000 description 34
- 230000000295 complement effect Effects 0.000 description 33
- 230000014509 gene expression Effects 0.000 description 26
- 230000008569 process Effects 0.000 description 26
- 230000005484 gravity Effects 0.000 description 24
- 230000001133 acceleration Effects 0.000 description 20
- 238000011156 evaluation Methods 0.000 description 15
- 239000000446 fuel Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 230000009467 reduction Effects 0.000 description 9
- 238000006557 surface reaction Methods 0.000 description 9
- 238000010187 selection method Methods 0.000 description 6
- 239000000725 suspension Substances 0.000 description 5
- 238000009499 grossing Methods 0.000 description 4
- 238000000137 annealing Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 108010077840 Complement C3a Proteins 0.000 description 2
- 108010078018 Complement C3d Proteins 0.000 description 2
- 108010028777 Complement C8 Proteins 0.000 description 2
- 102000016916 Complement C8 Human genes 0.000 description 2
- 102000020897 Formins Human genes 0.000 description 2
- 108091022623 Formins Proteins 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 230000035939 shock Effects 0.000 description 2
- 102000016917 Complement C1 Human genes 0.000 description 1
- 108010028774 Complement C1 Proteins 0.000 description 1
- 108090000955 Complement C2 Proteins 0.000 description 1
- 102000004381 Complement C2 Human genes 0.000 description 1
- 108010028773 Complement C5 Proteins 0.000 description 1
- 108010028771 Complement C6 Proteins 0.000 description 1
- 108010028776 Complement C7 Proteins 0.000 description 1
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 229910001416 lithium ion Inorganic materials 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 229910052987 metal hydride Inorganic materials 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- PXHVJJICTQNCMI-UHFFFAOYSA-N nickel Substances [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 1
- -1 nickel metal hydride Chemical class 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/72—Electric energy management in electromobility
Landscapes
- Regulating Braking Force (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Steering-Linkage Mechanisms And Four-Wheel Steering (AREA)
- Arrangement And Driving Of Transmission Devices (AREA)
- Steering Control In Accordance With Driving Conditions (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
Abstract
Description
本発明は、四輪独立駆動車の駆動力配分装置に関するものである。 The present invention relates to a driving force distribution device for a four-wheel independent drive vehicle.
車輪を駆動または制動するアクチュエータとしてのモータを各輪毎に有する四輪独立駆動車の駆動力配分装置が知られている(特許文献1参照)。 A driving force distribution device for a four-wheel independent drive vehicle having a motor as an actuator for driving or braking a wheel for each wheel is known (see Patent Document 1).
この装置では、4個の駆動輪のうち一輪のみがスリップしているときに、左側及び右側のうちスリップ輪と同じ側にある非スリップ輪に、スリップが発生していなければスリップ輪に配分されるはずであった出力トルクを配分する。また、スリップ輪が二輪あり、それらが左側及び右側に1個ずつあるときには、スリップしていなければそのスリップ輪に配分されるはずであったトルク出力を、同じ側にある非スリップ輪に配分するようにしている。以上により、各輪の制駆動力を補正する前後で、前後方向の車両加速度の変化と、各輪の制駆動力によって発生する車両重心周りのヨーモーメントの変化とを抑えるようにしている。
しかしながら、上記特許文献1の技術は、車両の左右各々の車輪において前後輪に配分している制駆動力を変化させないよう補正して車両重心周りのヨーモーメントの変化を抑えるものであり、制駆動力補正の前後における各輪と路面との間で発生している横方向力の変化を考慮していない。このため、前輪および後輪でそれぞれ発生している横方向力が制駆動力補正前後において大きく変化し、横方向の加速度と、それに基づく車両重心周りのヨーモーメントとが変化してしまう場合がある。こうした横方向の加速度の変化やその加速度変化に基づく車両重心周りのヨーモーメント変化はドライバーが意図しない変化であり、運転性を損なう恐れがある。
However, the technique disclosed in
そこで本発明は、各輪のスリップやモータの過熱等で特に二輪以上の制駆動力が変化した場合あるいは特に二輪以上の制駆動力を変化させる場合においても、車両挙動の変化をゼロとするような制駆動力の再配分を行う装置を提供することを目的とする。 Therefore, the present invention makes the change in vehicle behavior zero even when the braking / driving force of two or more wheels changes due to slip of each wheel, overheating of the motor or the like, or particularly when the braking / driving force of two or more wheels is changed. An object of the present invention is to provide a device that redistributes braking / driving force.
本発明は、四輪を独立に駆動可能な車両において、車両の平面運動での目標となる車両挙動を決定し、この目標となる車両挙動を実現する四輪の制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を演算し、この演算された制駆動力配分の集合の中から一つの制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択するように構成する。 The present invention determines a target vehicle behavior in a plane motion of a vehicle in a vehicle capable of independently driving the four wheels, and a set of four-wheel braking / driving force distributions {Fx for realizing the target vehicle behavior. 1 (j), Fx 2 (j), Fx 3 (j), Fx 4 (j)} are calculated, and one braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } are selected.
本発明によれば、車両の平面運動での目標となる車両挙動を決定し、この目標となる車両挙動を実現する四輪の制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を演算し、この演算された制駆動力配分の集合の中から一つの制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択するように構成したので、各輪のスリップ(各輪を駆動または制動する各アクチュエータがモータである場合には各モータの過熱)等によって特に二輪以上の制駆動力が変化した場合あるいは二輪以上の制駆動力を変化させる場合においても、従来技術より車両挙動の変化をゼロまたはそれに近い値に抑えることができる。 According to the present invention, a target vehicle behavior in the planar movement of the vehicle is determined, and a set of four-wheel braking / driving force distributions {Fx 1 (j), Fx 2 (j ), Fx 3 (j), Fx 4 (j)}, and one braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } from the set of calculated braking / driving force distributions Therefore, when the braking / driving force of two or more wheels changes due to slip of each wheel (overheating of each motor when each actuator that drives or brakes each motor is a motor) or the like, or two wheels Even when the above braking / driving force is changed, the change in the vehicle behavior can be suppressed to zero or a value close thereto as compared with the prior art.
以下、本発明の実施の形態を説明する。 Embodiments of the present invention will be described below.
図1は、本実施形態の四輪独立駆動車の概略構成図である。 FIG. 1 is a schematic configuration diagram of a four-wheel independent drive vehicle of the present embodiment.
バッテリ9からは、インバータ31〜34を介してモータ11〜14に電力が供給され、このバッテリ9より供給される電力によりモータ11が左前輪1を、モータ12が右前輪2を、モータ13が左後輪3を、モータ14が右後輪4を独立に駆動しあるいは独立に制動する。電力を供給するバッテリ9としてはニッケル水素電池あるいはリチウムイオン電池が望ましい。
Electric power is supplied from the
上記各モータ11〜14(各輪1〜4を駆動または制動するアクチュエータ)は三相同期電動機や三相誘導電動機等の力行運転及び回生運転ができる交流機である。
The
上記インバータ31〜34はモータ1〜4で発電された交流電流を直流電流に変換しバッテリ9に充電する、あるいはバッテリ9が放電した直流電流を交流電流に変換しモータ1〜4に供給する。
The
各輪1〜4の回転半径は所定値Rで全て等しく、各モータ11〜14と各車輪1〜4間は減速比1つまり直接連結されている。またさらに、各輪1〜4の輪荷重と横滑り角と路面摩擦係数が四輪で等しい場合には、制駆動力とタイヤ横力との関係は四輪で同一となる。つまり、四輪とも同じタイヤ特性を有する。
The rotation radii of the
前輪1,2は、ステアリングギヤ15を介してステアリング5の操舵により転舵可能であり、その舵角は運転者によるステアリング5の操舵によりステアリングギヤ15を介して機械的に調整される。ここで、前輪1,2の舵角変化量はステアリング5の操舵角変化量に対して1/16になるように設定されている。一方、後輪3,4の舵角は、コントローラ8から送信された指令値に追従するようステアリングアクチュエータ16によって調整される。
The
CPU、ROM、RAM、インターフェース回路及びインバータ回路等からなるコントローラ8には、車輪速センサ21〜24によって検出される各輪1〜4の速度、ステアリング角センサ25によって検出される運転者によるステアリング5の回転角、アクセルストロークセンサ26によって検出されるアクセルペダル6の踏込量、ブレーキストロークセンサ27によって検出されるブレーキペダル7の踏込量、舵角センサ41〜44によって検出される各輪1〜4の舵角、車両重心位置に取り付けられた加速度センサ100によって検出される車両の横方向加速度、ヨーレートセンサ101によって検出される車両のヨーレート等の信号が、各モータ11〜14に内蔵されている温度センサ(図示しない)によって検出される各モータ11〜14の温度の信号と共に入力され、これらの信号に基づいて各モータ11〜14にトルク配分を行う等の制御を行う。
A
このような四輪独立駆動車において、各輪1〜4のスリップや各モータ11〜14の過熱等によって、いずれか一輪の制駆動力が変化した場合あるいはいずれか一輪の制駆動力を任意に変化させる場合に、車両前後方向及び横方向の各加速度と車両重心周りのヨーモーメントとが変化しないように残り三輪の制駆動力補正量を求める技術(先願装置)を本発明者が既に提案している。この先願装置における制駆動力補正量の求め方について図2を用いて先に説明する。
In such a four-wheel independent drive vehicle, when the braking / driving force of any one wheel changes due to the slip of each
図2は図1に示した四輪独立駆動車をモデル化したものである。すなわち、左前輪1、右前輪2、左後輪3、右後輪4をそれぞれ独立に駆動できる車両において、各輪1〜4の制駆動力と、タイヤ横力と、舵角、そして車両重心周りのヨーモーメントを表した図である。当該車両ではサスペンションに加わる力によって発生するホイールアライメント変化や舵角変化等は十分小さいものと仮定し、これらサスペンション特性は無視している。
FIG. 2 is a model of the four-wheel independent drive vehicle shown in FIG. That is, in a vehicle that can independently drive the
図2において、各輪1〜4それぞれの舵角δ1,δ2,δ3,δ4[rad]は車両を鉛直上方から見た場合に時計回りを正とし、各舵角δ1,δ2,δ3,δ4は、各輪1〜4の回転方向が車両前後方向と一致している状態をゼロとする。各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4[N]は各舵角δ1,δ2,δ3,δ4が全てゼロのときに車両を前方に加速させる方向を正とし、各輪1〜4のタイヤ横力Fy1,Fy2,Fy3,Fy4[N]は各舵角δ1,δ2,δ3,δ4が全てゼロのときに車両を左方向に加速させる方向を正とする。
In FIG. 2, the steering angles δ 1 , δ 2 , δ 3 , and δ 4 [rad] of the
なお、各輪1〜4の舵角をδi(i=1〜4)、各輪1〜4の制駆動力をFxi(i=1〜4)、各輪1〜4のタイヤ横力をFyi(i=1〜4)と表すことがあり、記号δ,Fx,Fyに添える整数iについては左前輪を1で、右前輪2で、左後輪を3で、右後輪を4で表すものとする。記号に付す整数iは他の記号についても同じ意味で使用する。
The steering angle of each
また、車両のヨー回転方向の車両重心位置から前輪車軸までの距離をLf[m]、ヨー回転方向の車両重心位置から後輪車軸までの距離をLr[m]、前後輪のトレッド長さを共にLt[m]とする。従って、車両のホイールベース長さLl[m]はLl=Lf+Lrである。 Further, the distance from the vehicle center of gravity position in the yaw rotation direction of the vehicle to the front wheel axle is L f [m], the distance from the vehicle gravity center position in the yaw rotation direction to the rear wheel axle is L r [m], and the tread length of the front and rear wheels Let both be L t [m]. Therefore, the wheel base length L l [m] of the vehicle is L l = L f + L r .
ここで、まず各輪1〜4で発生する制駆動力とタイヤ横力の合力(つまりタイヤ力)の車体前後方向成分Fx’i(i=1〜4)及び車体横方向成分Fy’i(i=1〜4)を考える。すると、図3のように各輪1〜4の舵角をδi(i=1〜4)だけ切った場合におけるタイヤ力の車体前後方向成分Fx’iとタイヤ力の車体横方向成分Fy’iは次の式(1)及び式(2)で表される。ただし、タイヤ力の車体前後方向成分Fx’iは車両を前方に加速させる方向を、タイヤ力の車体横方向成分Fy’iは車両を左方向に加速させる方向をそれぞれ正とする。
Here, first, the vehicle longitudinal component Fx ′ i (i = 1 to 4) and the vehicle lateral component Fy ′ i (the tire force) of the braking / driving force and tire lateral force generated in each
Fx’i=Fxicosδi−Fyisinδi …(1)
Fy’i=Fxisinδi+Fyicosδi …(2)
従って、各輪1〜4の制駆動力が微小な量であるΔFxi(i=1〜4)だけ変化したときのタイヤ横力の変化量をΔFyi(i=1〜4)とすると、車体前後方向成分Fx’iの変化量ΔFx’i(i=1〜4),車体横方向成分Fy’iの変化量ΔFy’i(i=1〜4)は次の式(3)及び式(4)で表される。
Fx ′ i = Fx i cosδ i −Fy i sinδ i (1)
Fy ′ i = Fx i sin δ i + Fy i cos δ i (2)
Accordingly, when the amount of change in the tire lateral force when the braking / driving force of each
ΔFx’i=ΔFxicosδi−ΔFyisinδi …(3)
ΔFy’i=ΔFxisinδi+ΔFyicosδi …(4)
ここで、さらに各輪1〜4について制駆動力とタイヤ横力の関係は図4に示す関係にある。図4は輪荷重と路面摩擦係数に変化が無いとしたときの制駆動力(図では「駆動力」で略記。)とタイヤ横力の関係を表した図で、制駆動力を横軸に、タイヤ横力を縦軸にとっている。この図4の関係を利用して、各輪1〜4の現在の制駆動力Fxi(i=1〜4)とタイヤ横力Fyi(i=1〜4)における、制駆動力変化ΔFxi(i=1〜4)に対するタイヤ横力の感度をki(i=1〜4)とおく。つまり、制駆動力変化に対するタイヤ横力の感度kiは図4に示すように車体前後方向成分の変化量ΔFxi及び車体横方向成分の変化量ΔFyi(i=1〜4)が微小のときの次の式(5)により定義される値である。
ΔFx ′ i = ΔFx i cosδ i −ΔFy i sinδ i (3)
ΔFy ′ i = ΔFx i sin δ i + ΔFy i cos δ i (4)
Here, the relationship between the braking / driving force and the tire lateral force for each of the
ki=ΔFyi/ΔFxi …(5)
すると、車体前後方向成分の変化量ΔFxi及び車体横方向成分の変化量ΔFyiが微小で、この式(5)の近似が十分成り立つとすると、ΔFyi=ki・ΔFxiとおけるので、各輪1〜4の制駆動力Fxiが十分微小なΔFxiだけ変化したときの車体前後方向成分Fx’iの変化量ΔFx’iと車体横方向成分Fy’iの変化量ΔFy’iとは次の式(6)及び式(7)へと変形される。
k i = ΔFy i / ΔFx i (5)
Then, a variation DerutaFy i variation DerutaFx i and the vehicle body lateral component of the vehicle longitudinal direction component is small, the approximation of the equation (5) is sufficiently satisfied, Okeru the ΔFy i = k i · ΔFx i , variation ΔFy'i and 'variation ΔFx'i the vehicle body lateral component Fy of i' i vehicle longitudinal direction component Fx when the longitudinal force Fx i of each wheel 1-4 changes by enough small DerutaFx i Is transformed into the following equations (6) and (7).
ΔFx’i=(cosδi−kisinδi)ΔFxi=piΔFxi …(6)
ΔFy’i=(sinδi+kicosδi)ΔFxi=qiΔFxi …(7)
ここで、式(6)、式(7)の係数pi,qi(i=1〜4)は、舵角δiと、制駆動力変化に対するタイヤ横力の感度kiにより表される次の値である。
ΔFx ′ i = (cosδ i −k i sinδ i ) ΔFx i = p i ΔFx i (6)
ΔFy ′ i = (sin δ i + k i cos δ i ) ΔFx i = q i ΔFx i (7)
Here, the coefficients p i and q i (i = 1 to 4) in the equations (6) and (7) are expressed by the steering angle δ i and the sensitivity k i of the tire lateral force with respect to the braking / driving force change. It is the next value.
pi=cosδi−ki・sinδi、 …(補1)
qi=sinδi+ki・cosδi、 …(補2)
さて、上記の各輪1〜4の車体前後方向成分Fx’i及び車体横方向成分Fy’iを用いると、図2の状態において、タイヤ力の総和の車体前後方向成分Fxと、タイヤ力の総和の車体横方向成分Fyと、各輪のタイヤ力によって発生する車両重心周りのヨーモーメントの総和Mとは、次の式(8)〜(10)により表すことができる。ただし、各輪1〜4のタイヤ力によって発生する車両重心周りのヨーモーメントの総和Mは、図2のように車両を鉛直上方からみたときに反時計回りを正とする。
p i = cos δ i −k i · sin δ i , (Supplement 1)
q i = sin δ i + k i · cos δ i , (Supplement 2)
Now, when the vehicle body longitudinal component Fx ′ i and the vehicle lateral component Fy ′ i of each of the
Fx=Fx’1+Fx’2+Fx’3+Fx’4 …(8)
Fy=Fy’1+Fy’2+Fy’3+Fy’4 …(9)
M={(Fx’2+Fx’4)−(Fx’1+Fx’3)}×Lt/2
+{(Fy’1+Fy’2)×Lf−(Fy’3+Fy’4)×Lr}
…(10)
従って、各輪1〜4の制駆動力Fxiがそれぞれ微小な量であるΔFxiだけ変化したときのタイヤ力の総和の車体前後方向成分Fx,タイヤ力の総和の車体横方向成分Fy,各輪のタイヤ力によって発生する車両重心周りのヨーモーメントの総和Mの各変化量ΔFx,ΔFy,ΔMは上記の式(補1)及び式(補2)の係数pi,qiを用いて、次の式(11)〜(13)により表される。
Fx = Fx ′ 1 + Fx ′ 2 + Fx ′ 3 + Fx ′ 4 (8)
Fy = Fy ′ 1 + Fy ′ 2 + Fy ′ 3 + Fy ′ 4 (9)
M = {(Fx ′ 2 + Fx ′ 4 ) − (Fx ′ 1 + Fx ′ 3 )} × L t / 2.
+ {(Fy ′ 1 + Fy ′ 2 ) × L f − (Fy ′ 3 + Fy ′ 4 ) × L r }
(10)
Therefore, the vehicle body lateral component Fy in the vehicle longitudinal direction component Fx, the sum of the tire forces of the sum of the tire forces when longitudinal force Fx i of each wheel 1-4 has changed by DerutaFx i is a small amount, respectively, each The change amounts ΔFx, ΔFy, ΔM of the total sum M of yaw moments around the center of gravity of the vehicle generated by the wheel tire force are calculated using the coefficients p i , q i in the above equations (complement 1) and equation (complement 2), It is represented by the following formulas (11) to (13).
ΔFx=ΔFx’1+ΔFx’2+ΔFx’3+ΔFx’4
=p1ΔFx1+p2ΔFx2+p3ΔFx3+p4ΔFx4
…(11)
ΔFy=ΔFy’1+ΔFy’2+ΔFy’3+ΔFy’4
=q1ΔFx1+q2ΔFx2+q3ΔFx3+q4ΔFx4
…(12)
ΔM={(ΔFx’2+ΔFx’4)−(ΔFx’1+ΔFx’3)}×Lt/2
+{(ΔFy’1+ΔFy’2)×Lf−(ΔFy’3+ΔFy’4)×Lr}
=(−p1Lt/2+q1Lf)ΔFx1
+(p2Lt/2+q2Lf)ΔFx2
+(−p3Lt/2+q3Lf)ΔFx3
+(p4Lt/2+q4Lf)ΔFx4 …(13)
これら3つの式(11)〜(13)をまとめると次の行列式(14)で表すことができる。
ΔFx = ΔFx ′ 1 + ΔFx ′ 2 + ΔFx ′ 3 + ΔFx ′ 4
= P 1 ΔFx 1 + p 2 ΔFx 2 + p 3 ΔFx 3 + p 4 ΔFx 4
... (11)
ΔFy = ΔFy ′ 1 + ΔFy ′ 2 + ΔFy ′ 3 + ΔFy ′ 4
= Q 1 ΔFx 1 + q 2 ΔFx 2 + q 3 ΔFx 3 + q 4 ΔFx 4
(12)
ΔM = {(ΔFx ′ 2 + ΔFx ′ 4 ) − (ΔFx ′ 1 + ΔFx ′ 3 )} × L t / 2.
+ {(ΔFy ′ 1 + ΔFy ′ 2 ) × L f − (ΔFy ′ 3 + ΔFy ′ 4 ) × L r }
= (- p 1 L t / 2 + q 1 L f)
+ (P 2 L t / 2 + q 2 L f ) ΔFx 2
+ (− P 3 L t / 2 + q 3 L f ) ΔFx 3
+ (P 4 L t / 2 + q 4 L f ) ΔFx 4 (13)
These three formulas (11) to (13) can be summarized by the following determinant (14).
上記の行列式(14)の左辺、つまりタイヤ力の総和の車体前後方向成分の変化量ΔFx、タイヤ力の総和の車体横方向成分の変化量ΔFy、各輪のタイヤ力によって発生する車両重心周りのヨーモーメントの総和の変化量ΔMをすべてゼロとした次の行列式(15)を満たす制駆動力の変化量ΔFx2,ΔFx3,ΔFx4は、行列式(15)をΔFx2,ΔFx3,ΔFx4に関する連立方程式と見立てて解くと、左前輪1の制駆動力の変化量ΔFx1を用いて次の式(16)〜式(18)により表される。
The left side of the determinant (14), that is, the change amount ΔFx of the vehicle front-rear direction component of the sum of tire forces, the change amount ΔFy of the vehicle side direction component of the sum of tire forces, and the vehicle center of gravity generated by the tire force of each wheel The change amounts ΔFx 2 , ΔFx 3 , and ΔFx 4 of the braking / driving force that satisfy the following determinant (15) with the total change amount ΔM of the yaw moment of zero being the following determinants (15) are ΔFx 2 , ΔFx 3 , ΔFx 4 is expressed as the following equations (16) to (18) using the change amount ΔFx 1 of the braking / driving force of the left
ΔFx2=[{q1(p4q3−p3q4)L1+p4(p3q1−p1q3)Lt}
/{q2(p3q4−p4q3)Ll+p3(p2q4−p4q2)Lt}]
・ΔFx1
…(16)
ΔFx3=[{q4(p2q1−p1q3)Ll+p1(p4q2−p2q4)Lt}
/{q2(p3q4−p4q3)Ll+p3(p2q4−p4q2)Lt}]
・ΔFx1
…(17)
ΔFx4=[{q3(p1q2−p2q1)Ll+p2(p1q3−p3q1)Lt}
/{q2(p3q4−p4q3)Ll+p3(p2q4−p4q2)Lt}]
・ΔFx1
…(18)
従って、これら3つの式(16)〜式(18)から明らかなように、各輪1〜4の制駆動力の変化量ΔFx1,ΔFx2,ΔFx3,ΔFx4が次の式(19)の比となる場合には、タイヤ力の総和の車体前後方向成分の変化量ΔFx、タイヤ力の総和の車体横方向成分の変化量ΔFy、各輪のタイヤ力によって発生する車両重心周りのヨーモーメントの総和の変化量ΔMが全てゼロとなる。
ΔFx 2 = [{q 1 (p 4 q 3 −p 3 q 4 ) L 1 + p 4 (p 3 q 1 −p 1 q 3 ) L t }
/ {Q 2 (p 3 q 4 -p 4 q 3) L l + p 3 (p 2 q 4 -p 4 q 2) L t}]
・ ΔFx 1
... (16)
ΔFx 3 = [{q 4 (p 2 q 1 −p 1 q 3 ) L 1 + p 1 (p 4 q 2 −p 2 q 4 ) L t }
/ {Q 2 (p 3 q 4 -p 4 q 3) L l + p 3 (p 2 q 4 -p 4 q 2) L t}]
・ ΔFx 1
... (17)
ΔFx 4 = [{q 3 (p 1 q 2 −p 2 q 1 ) L 1 + p 2 (p 1 q 3 −p 3 q 1 ) L t }
/ {Q 2 (p 3 q 4 -p 4 q 3) L l + p 3 (p 2 q 4 -p 4 q 2) L t}]
・ ΔFx 1
... (18)
Accordingly, as is apparent from these three equations (16) to (18), the changes ΔFx 1 , ΔFx 2 , ΔFx 3 , and
ΔFx1:ΔFx2:ΔFx3:ΔFx4
=q2(p3q4−p4q3)Ll+p3(p2q4−p4q2)Lt
:q1(p4q3−p3q4)Ll+p4(p3q1−p1q3)Lt
:q4(p2q1−p1q2)Ll+p1(p4q2−p2q4)Lt
:q3(p1q2−p2q1)Ll+p2(p1q3−p3q1)Lt
={(Lt/Ll)(h4−h2)+h2(h4−h3)}
/(cosδ1−k1sinδ1)
:{−(Lt/Ll)(h3−h1)−h1(h4−h3)}
/(cosδ2−k2sinδ2)
:{−(Lt/Ll)(h4−h2)−h4(h2−h1)}
/(cosδ3−k3sinδ3)
:{(Lt/Ll)(h3−h1)+h3(h2−h1)}
/(cosδ4−k4sinδ4)
…(19)
ここで、式(19)の係数h1,h2,h3,h4は、各輪1〜4の舵角δ1,δ2,δ3,δ4と、制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4とから次の式(補3a)〜式(補3d)により表される値である。
ΔFx 1 : ΔFx 2 : ΔFx 3 : ΔFx 4
= Q 2 (p 3 q 4 -p 4 q 3) L l + p 3 (p 2 q 4 -p 4 q 2) L t
: Q 1 (p 4 q 3 -p 3 q 4) L l + p 4 (p 3 q 1 -p 1 q 3) L t
: Q 4 (p 2 q 1 -p 1 q 2) L l + p 1 (p 4 q 2 -p 2 q 4) L t
: Q 3 (p 1 q 2 -p 2 q 1) L l + p 2 (p 1 q 3 -p 3 q 1) L t
= {(L t / L l ) (h 4 −h 2 ) + h 2 (h 4 −h 3 )}
/ (Cosδ 1 −k 1 sinδ 1 )
: {- (L t / L l) (h 3 -h 1) -h 1 (h 4 -h 3)}
/ (Cosδ 2 −k 2 sinδ 2 )
: {- (L t / L l) (h 4 -h 2) -h 4 (h 2 -h 1)}
/ (Cosδ 3 −k 3 sinδ 3 )
: {(L t / L l ) (h 3 −h 1 ) + h 3 (h 2 −h 1 )}
/ (Cosδ 4 −k 4 sinδ 4 )
... (19)
Here, the coefficients h 1 , h 2 , h 3 , h 4 in the equation (19) are the steering angles δ 1 , δ 2 , δ 3 , δ 4 of the wheels 1 to 4 and the side of the tire against the change in braking / driving force. These are values expressed by the following equations (complement 3a) to (complement 3d) from the force sensitivities k 1 , k 2 , k 3 , and k 4 .
h1=(sinδ1+k1cosδ1)/(cosδ1−k1sinδ1) …(補3a)
h2=(sinδ2+k2cosδ2)/(cosδ2−k2sinδ2) …(補3b)
h3=(sinδ3+k3cosδ3)/(cosδ3−k3sinδ3) …(補3c)
h4=(sinδ4+k4cosδ4)/(cosδ4−k4sinδ4) …(補3d)
また、各輪1〜4の舵角δ1,δ2,δ3,δ4が十分小さい場合には、式(19)において各輪の舵角δ1≒0,δ2≒0,δ3≒0,δ4≒0とすることにより、各輪1〜4の制駆動力変化量ΔFx1,ΔFx2,ΔFx3,ΔFx4の比を式(19)に代えて次の式(20)のように求めることができる。
h 1 = (sin δ 1 + k 1 cos δ 1 ) / (cos δ 1 −k 1 sin δ 1 ) (Supplement 3a)
h 2 = (sin δ 2 + k 2 cos δ 2 ) / (cos δ 2 −k 2 sin δ 2 ) (Supplement 3b)
h 3 = (sin δ 3 + k 3 cos δ 3 ) / (cos δ 3 −k 3 sin δ 3 ) (Supplement 3c)
h 4 = (sin δ 4 + k 4 cos δ 4 ) / (cos δ 4 −k 4 sin δ 4 ) (Supplement 3d)
When the steering angles δ 1 , δ 2 , δ 3 , δ 4 of the wheels 1 to 4 are sufficiently small, the steering angles δ 1 ≈0, δ 2 ≈0, δ 3 of each wheel in the equation (19). ≒ 0, [delta] 4 with ≒ 0, longitudinal force variation DerutaFx 1 of each
ΔFx1:ΔFx2:ΔFx3:ΔFx4
=(Lt/Ll)(k4−k2)+k2(k4−k3)
:−(Lt/Ll)(k3−k1)−k1(k4−k3)
:−(Lt/Ll)(k4−k2)−k4(k2−k1)
:(Lt/Ll)(k3−k1)+k3(k2−k1)
…(20)
このように先願装置では、四輪1〜4のうちいずれか一輪の制駆動力が各輪1〜4のスリップや各モータ11〜14の故障等で変化した場合あるいは四輪1〜4のうちいずれか一輪の制駆動力を任意に変化させた場合に、制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4がこれら式(19)あるいは式(20)に表した比となるように制御する装置を提案している。即ち、サスペンション特性を無視した車両の平面運動において、一輪例えば左前輪1の制駆動力補正量ΔFx1を用い、上記の式(19)あるいは式(20)に従って残り三輪の制駆動力補正量ΔFx2,ΔFx3,ΔFx4を算出し、これら制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4で制駆動力Fx1,Fx2,Fx3,Fx4を補正することにより、タイヤ力の総和の車体前後方向成分Fxと、タイヤ力の総和の車体横方向成分Fyと、各輪のタイヤ力によって発生する車体重心周りのヨーモーメントの総和Mの変化をすべてゼロにすることができる装置を創出している。
ΔFx 1 : ΔFx 2 : ΔFx 3 : ΔFx 4
= (L t / L l ) (k 4 −k 2 ) + k 2 (k 4 −k 3 )
:-( L t / L l) ( k 3 -k 1) -k 1 (k 4 -k 3)
:-( L t / L l) ( k 4 -k 2) -k 4 (k 2 -k 1)
: (L t / L l ) (k 3 −k 1 ) + k 3 (k 2 −k 1 )
... (20)
In this way, in the prior application device, when the braking / driving force of any one of the four
さて、先願装置では、四輪1〜4のうちいずれか一輪の制駆動力が各輪1〜4のスリップや各モータ11〜14の過熱等で変化した場合あるいは四輪1〜4のうちいずれか一輪の制駆動力を任意に変化させた場合のみ考えており、四輪1〜4のうち二輪以上の制駆動力がともに変化した場合あるいは四輪1〜4のうち二輪以上の制駆動力を任意に変化させた場合には対応できない。
In the prior application device, when the braking / driving force of any one of the four
これについて説明すると、先願装置の適用対象としているサスペンション特性を無視した四輪独立駆動車の平面運動について再考してみる。 To explain this, let us reconsider the planar motion of a four-wheel independent drive vehicle ignoring the suspension characteristics to which the prior application device is applied.
四輪独立駆動車は4つの独立な駆動輪1〜4を有するので、自由度は4つである。これに対して、平面運動では車両の運動の方向は3つ(車両前後方向,車両横方向,ヨー回転方向)であるから、自由度は3つである。これは上記の式(8)〜式(10)において、制駆動力配分{Fx1,Fx2,Fx3,Fx4}によってタイヤ力の総和の車体前後方向成分Fx、タイヤ力の総和の車体横方向成分Fy、各輪1〜4のタイヤ力によって発生する車両重心周りのヨーモーメントの総和Mが決定されることと同意である。ここで、車両挙動は、タイヤ力の総和の車体前後方向成分Fx、タイヤ力の総和の車体横方向成分Fy及び各輪1〜4のタイヤ力によって発生する車両重心周りのヨーモーメントの総和Mの3つによって一義的に決まるので、これらFx、Fy及びMの3つを以下「車両挙動」で総称する。
Since the four-wheel independent drive vehicle has four
従って、上記の式(8)〜式(10)を、4つの制駆動力Fx1,Fx2,Fx3,Fx4を未知数とした連立方程式と見立てると、その拡大係数行列のサイズは3行5列になるので、階数は必ず未知数Fx1,Fx2,Fx3,Fx4の数である4より小さい。これは、つまり、上記の式(8)〜式(10)によって一つの車両挙動Fx,Fy,Mが定まるとき、その定まる車両挙動と同じ車両挙動Fx,Fy,Mを実現する制駆動力配分{Fx1,Fx2,Fx3,Fx4}は無数に存在することを意味する。 Therefore, if the above equations (8) to (10) are regarded as simultaneous equations with the four braking / driving forces Fx 1 , Fx 2 , Fx 3 , and Fx 4 as unknowns, the size of the expansion coefficient matrix is 3 rows. Since there are 5 columns, the rank is always smaller than 4 , which is the number of unknowns Fx 1 , Fx 2 , Fx 3 , Fx 4 . In other words, when one vehicle behavior Fx, Fy, M is determined by the above equations (8) to (10), the braking / driving force distribution that realizes the same vehicle behavior Fx, Fy, M as the determined vehicle behavior is determined. {Fx 1 , Fx 2 , Fx 3 , Fx 4 } means that there are countless numbers.
一方、先願装置で対象としている各輪1〜4のスリップや各モータ11〜14の過熱等によって各輪1〜4の制駆動力が変化した場合あるいは各輪1〜4の制駆動力を変化させる場合を考える。
On the other hand, when the braking / driving force of each
例えば左前輪1にスリップが発生して、路面に伝達できる制駆動力が所定値a[N](a>0)減少した場合に、先願装置では左前輪1のスリップを収束させるためにこの減少代aを左前輪1の制駆動力補正量ΔFx1として(ΔFx1=−a)、残り三輪の制駆動力補正量ΔFx2,ΔFx3,ΔFx4を上記の式(19)を用いて算出し、これら各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4で対応する各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4を補正することとなる。
For example, when slip occurs in the left
しかしながら、左前輪1のスリップを収束させるためには、この減少代a以上に左前輪1の制駆動力を減少させてもよい。つまり、左前輪1のスリップを収束させるためには左前輪1の制駆動力Fx1を先願装置のようにFx1+ΔFx1=Fx1−aに設定するのではなく、Fx1−a以下にすればよいので、例えば左前輪1の制駆動力補正量ΔFx1をΔFx1=−b(b>a)のようにおいて残り三輪2〜4の制駆動力補正量ΔFx2,ΔFx3,ΔFx4を上記の式(19)を用いて求めても本来の目的である左前輪1のスリップを収束させながら車両挙動の変化をゼロとすることが可能である。
However, in order to cause the slip of the left
ここで、左前輪1の制駆動力補正量ΔFx1をΔFx1=−aとしたときと、左前輪1の制駆動力補正量ΔFx1をΔFx1=−bとしたときとでは、残り三輪2〜4の制駆動力補正量ΔFx2,ΔFx3,ΔFx4が異なる。つまり、制駆動力配分{Fx1,Fx2,Fx3,Fx4}が、左前輪1の制駆動力補正量ΔFx1をΔFx1=−aとしたときと、左前輪1の制駆動力補正量ΔFx1をΔFx1=−bとしたときとで異なる。
Here, when the braking / driving force correction amount ΔFx 1 of the left
すると、左前輪1がスリップする上記のケースにおいて、さらに左前輪1以外のもう一輪つまり右前輪2、左後輪3、右後輪4のうちのいずれか一輪がスリップしてそのスリップしたもう一輪の制駆動力を減少させる必要がある場合には、左前輪1の制駆動力補正量ΔFx1をΔFx1=−aとしたときの制駆動力配分{Fx1,Fx2,Fx3,Fx4}では、このもう一輪のスリップを収束させることができないが、左前輪1の制駆動力補正量ΔFx1をΔFx1=−bとしたときの制駆動力配分{Fx1,Fx2,Fx3,Fx4}では、このもう一輪のスリップをも収束させ得るケースが存在する可能性がある。即ち、二輪のスリップや二つのモータの過熱等を抑えるためには、スリップやモータ過熱等が発生した二輪以上の車輪の制駆動力を適当な範囲に制限すればよいことから、四輪1〜4のうち二輪以上の制駆動力が変化した場合あるいは四輪1〜4のうち二輪以上の制駆動力を変化させる場合においても、無数に存在する同じ車両挙動を実現する制駆動力配分{Fx1,Fx2,Fx3,Fx4}の中に、この左前輪1ともう1輪との複数輪の制駆動力制限を満たす制駆動力配分{Fx1,Fx2,Fx3,Fx4}が存在する可能性があるのである。
Then, in the case where the left
しかしながら、このような複数輪の制駆動力制限を満たす制駆動力配分{Fx1,Fx2,Fx3,Fx4}が存在する場合においても、先願装置によれば車両挙動の変化をゼロにすることができない。 However, even when there is a braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } satisfying such a braking / driving force limitation of multiple wheels, according to the prior application device, the change in vehicle behavior is zero. I can't.
そこで本発明は、四輪1〜4のうちスリップやモータの過熱等で二輪以上の制駆動力が変化した場合あるいは四輪1〜4のうち二輪以上の制駆動力を変化させる場合においても、車両挙動の変化をゼロにすることが可能である場合には、車両挙動の変化をゼロとするように制駆動力の再配分を行う装置を提案するものである。
Therefore, in the case where the braking / driving force of two or more wheels changes due to slip, motor overheating or the like among the four
本発明によるこの制駆動力再配分を行う方装置の概要を次に説明する。 An outline of a device for performing this braking / driving force redistribution according to the present invention will be described below.
本発明によって制御される四輪独立駆動車は、まず要求された車両挙動(目標となる車両挙動)を実現する四輪の制駆動力配分の集合を求めるかあるいは予め求めておく。そして、スリップやモータの過熱等が発生し、四輪1〜4のうち二輪以上の制駆動力が変化した場合あるいは四輪1〜4のうち二輪以上の制駆動力を変化させる場合に、変化したあるいは変化させる複数の車輪の各輪が出力可能な制駆動力の範囲をそれぞれ求め、先に求めておいた制駆動力配分{Fx1,Fx2,Fx3,Fx4}の集合の中で、各輪の出力可能な制駆動力範囲を満たす一つの制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択することにより、制駆動力の再配分後の車両挙動が制駆動力の再配分前の車両挙動と同じになるように制駆動力を再配分する。
In the four-wheel independent drive vehicle controlled by the present invention, a set of four-wheel braking / driving force distributions for realizing the requested vehicle behavior (target vehicle behavior) is first obtained or obtained in advance. And when slip, overheating of the motor, etc. occur, the braking / driving force of two or more wheels among the four
この要求された車両挙動を実現する制駆動力配分{Fx1,Fx2,Fx3,Fx4}の集合は、次のようにして求める(演算する)ことができる。 A set of braking / driving force distributions {Fx 1 , Fx 2 , Fx 3 , Fx 4 } that realizes the requested vehicle behavior can be obtained (calculated) as follows.
まず、要求された車両挙動を実現する制駆動力配分{Fx1,Fx2,Fx3,Fx4}を少なくとも一つみつけるかあるいは設定する。次に、この一つみつけたか設定した制駆動力配分{Fx1,Fx2,Fx3,Fx4}に対し、いずれか一輪例えば左前輪1の制駆動力Fx1を十分微小な量であるΔFx1だけ変化させ、各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4から残り三輪2〜4の制駆動力補正量Fx2,Fx3,Fx4を上記の式(19)あるいは式(20)を用いて算出し、現在と同じ車両挙動Fx,Fy,Mを実現する別の制駆動力配分{Fx1,Fx2,Fx3,Fx4}を求める。この左前輪1の制駆動力をさらに微少に変化させて上記の式(19)あるいは式(20)による制駆動力配分{Fx1,Fx2,Fx3,Fx4}の補正処理を繰り返すことによって、要求された車両挙動を実現する制駆動力配分{Fx1,Fx2,Fx3,Fx4}の集合を求めることができる。
First, at least one braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } that realizes the requested vehicle behavior is found or set. Next, with respect to the braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } that has been found or set, the braking / driving force Fx 1 of any one wheel, for example, the left front wheel 1 is sufficiently small. By changing ΔFx 1, the tire lateral force sensitivity k 1 , k 2 , k 3 , k 4 with respect to changes in the braking / driving force of each
この方法に従って求めた現在と同じ車両挙動を実現する制駆動力配分{Fx1,Fx2,Fx3,Fx4}の集合の一例を図5に示すと、図5は前輪1,2のみを転舵する四輪独立駆動車において、前輪1,2を左方向に転舵し、走行抵抗に釣り合う制駆動力を輪荷重比で制駆動力配分を行うことによって左方向に定速旋回している走行状態での制駆動力配分として、丸印で表された制駆動力配分{Fx1(0),Fx2(0),Fx3(0),Fx4(0)}を設定し、この定速旋回を実現する制駆動力配分の集合を本発明の上記方法に従って求めた結果である。
FIG. 5 shows an example of a set of braking / driving force distributions {Fx 1 , Fx 2 , Fx 3 , Fx 4 } that realizes the same vehicle behavior obtained in accordance with this method. FIG. 5 shows only the
図5では左前輪1の制駆動力を横軸に採り、この左前輪1の制駆動力に対する四輪1〜4それぞれの制駆動力を縦軸に採って表している(図では制駆動力を「駆動力」と略記。)。図示のように、左前輪1の制駆動力の集合の軌跡は直線であるのに対して、残り三輪である右前輪2、左後輪3、右後輪4の制駆動力の集合の軌跡はいずれも閉じた輪になっている。
In FIG. 5, the braking / driving force of the left
ここで、図5では、丸印の制駆動力配分{Fx1(0),Fx2(0),Fx3(0),Fx4(0)}から、各輪1〜4の制駆動力を上記の式(19)の比が維持されるように微小に変化させながら求めている。従って、左前輪1の制駆動力Fx1がFx1=Fx1(0)の場合に、この定速旋回を実現する制駆動力配分としては、丸印で表された制駆動力配分{Fx1(0),Fx2(0),Fx3(0),Fx4(0)}のほかに、三角印で表された制駆動力配分{Fx1(0),Fx2(1),Fx3(1),Fx4(1)}が存在することがわかる。左前輪1の制駆動力Fx1がFx1=Fx1(0)の場合に、この定速旋回を実現する制駆動力配分が、このように2つの制駆動力配分のみとなるのは、Fx1(0)における一本の縦軸と、残り三輪の制駆動力配分の集合の軌跡を表す各閉じた輪との交点が最大2つになるためである。
Here, in FIG. 5, the braking / driving forces of the
ただし、例えば{Fx1(0),Fx2(1),Fx3(0),Fx4(1)}といった制駆動力配分ではこの定速旋回を実現することはできないことに注意すべきである。 However, it should be noted that this constant speed turning cannot be realized by braking / driving force distribution such as {Fx 1 (0), Fx 2 (1), Fx 3 (0), Fx 4 (1)}. is there.
そして、スリップやモータの過熱等で四輪のうち二輪以上の制駆動力が変化した場合あるいは四輪のうち二輪以上の制駆動力を変化させる場合には、この現在と同じ車両挙動を実現する制駆動力配分の集合の中から、各輪の出力可能な制駆動力範囲を満たす一つの制駆動力配分を選択すればよい。 And when the braking / driving force of two or more of the four wheels changes due to slip, motor overheating, etc., or when changing the braking / driving force of two or more of the four wheels, the same vehicle behavior as this present is realized. From the set of braking / driving force distributions, one braking / driving force distribution satisfying the range of braking / driving force that can be output from each wheel may be selected.
例えば、図5において四角印の制駆動力配分{Fx1(2),Fx2(2),Fx3(2),Fx4(2)}で定速旋回中に、左前輪1と右前輪2がともにスリップした場合を考える。このとき、左前輪1のスリップを収束させるためには左前輪1の制駆動力Fx1を所定の範囲、例えば−η≦Fx1≦η(η[N]>0)の範囲で制限し、かつ右前輪2のスリップを収束させるためには右前輪2の制駆動力Fx2を所定の範囲、例えば−ζ≦Fx2≦ζ(ζ[N]>0)の範囲で制限しなければならない。ここで、所定値ηは左前輪1の制駆動力上限を、−ηは左前輪1の制駆動力下限を表す。また、所定値ζは右前輪2の制駆動力上限を、−ζは右前輪2の制駆動力下限を表している。
For example, in FIG. 5, the left
この場合に、次のようにして一つの制駆動力配分を選択する。 In this case, one braking / driving force distribution is selected as follows.
まず、図6に示すように、左前輪1の制駆動力制限−η≦Fx1≦ηを満たす制駆動力配分の集合を求めると、この集合の軌跡は図6において太線で示した範囲(A)の線分及び弧となる。同様にして、図7に示すように右前輪2の制駆動力制限−ζ≦Fx2≦ζを満たす制駆動力配分の集合を求めると、この集合の軌跡は図7において太線で示した範囲(B)の線分及び弧となる。従って、左前輪1の制駆動力制限−η≦Fx1≦ηと、右前輪2の制駆動力制限−ζ≦Fx2≦ζとを同時に満たす制駆動力配分の集合の軌跡は図8において太線で示した範囲(C)の線分及び弧となるので、この太線で示した範囲(C)の中から任意の制駆動力配分を一つ選択すればよい。選択方法としては、例えば四輪をすべてモータで駆動しているような場合には、モータの総消費電力を最小化するような制駆動力配分を選択すればよい。
First, as shown in FIG. 6, when a set of braking / driving force distributions satisfying the braking / driving force limit -η ≦ Fx 1 ≦ η of the left
このように、本発明を用いることによってスリップやモータの過熱等で四輪のうち二輪以上の制駆動力が変化した場合あるいは四輪のうち二輪以上の制駆動力を変化させる場合においても、車両挙動の変化をゼロにすることが可能である場合には、車両挙動の変化をゼロとするような制駆動力の再配分を行うことが可能となることがわかる。 As described above, even when the braking / driving force of two or more of the four wheels changes due to slip, motor overheating, or the like by using the present invention, or when the braking / driving force of two or more of the four wheels is changed, the vehicle It can be seen that when the change in behavior can be made zero, the braking / driving force can be redistributed so that the change in vehicle behavior becomes zero.
コントローラ8で実行されるこの制御を以下のフローチャートを参照して詳述する。
This control executed by the
図9(A)及び図9(B)は各モータ11〜14へのトルク配分制御を実行するためのもので、操作の手順を記載している。一定時間毎に実行するものではない。
FIG. 9A and FIG. 9B are for executing torque distribution control to the
図9(A)においてステップ10では、車輪速センサ21〜24で各輪1〜4の回転速度ω1,ω2,ω3,ω4[rad/s]をそれぞれ検出し、各輪1〜4の半径Rを乗じて各車輪速度V1,V2,V3,V4[m/s]を得る。
In FIG. 9A, in
また、アクセルストロークセンサ26によってアクセルペダル6の踏込量AP[%]を、ブレーキストロークセンサ27によってブレーキペダル7の踏込量BP[%]を、ステアリング角センサ25によってステアリング5の回転角θ[rad]を、加速度センサ100によって車両の前後方向加速度Xg[m/s2]と横方向加速度Yg[m/s2]を、ヨーレートセンサ101によってヨーレートγ[rad/s]を、各モータ内蔵の温度センサによって各輪1〜4のモータ11〜14の各温度T1,T2,T3,T4[℃]をそれぞれ検出する。
Further, the
ここで、各車輪速度V1,V2,V3,V4は車両前進方向を正とし、ステアリングの回転角θは反時計回りを正とし、車両前後方向加速度Xgは車両が前方に加速する方向を正とし、車両横方向加速度Ygは車両が左旋回時に車両重心位置から旋回中心に向かう方向を正とし、ヨーレートγは車両を鉛直上方からみたときに反時計回りを正とする。 Here, the wheel speeds V 1 , V 2 , V 3 and V 4 are positive in the vehicle forward direction, the steering rotation angle θ is positive in the counterclockwise direction, and the vehicle longitudinal acceleration Xg is accelerated forward by the vehicle. The direction is positive, the vehicle lateral acceleration Yg is positive in the direction from the center of gravity of the vehicle toward the turning center when the vehicle is turning left, and the yaw rate γ is positive in the counterclockwise direction when the vehicle is viewed from vertically above.
ステップ20では、車体速度V[m/s]を次の式(21)により算出する。車体速度Vは車両前進方向を正とする。
In
V=(V1+V2+V3+V4)÷4 …(21)
ステップ30では、各輪1〜4の輪荷重W1,W2,W3,W4[N]を次の式(22)〜式(25)により算出する。
V = (V 1 + V 2 + V 3 + V 4 ) ÷ 4 (21)
In
W1=mLrg/2Ll−mhXg/2Ll−mhYg/2Lt
…(22)
W2=mLrg/2Ll−mhXg/2Ll+mhYg/2Lt
…(23)
W3=mLrg/2Ll+mhXg/2Ll−mhYg/2Lt
…(24)
W4=mLrg/2Ll+mhXg/2Ll+mhYg/2Lt
…(25)
ここで、上記の式(22)〜式(25)においてmは車両の質量[kg]、hは車両の重心高さ[m]、gは重力加速度[m/s2]である。各輪1〜4の輪荷重W1,W2,W3,W4の算出方法は特開2002−211378公報に記載されており、ここでは当該公報に従って求めている。
W 1 = mL r g / 2L l -mhXg / 2L l -mhYg / 2L t
... (22)
W 2 = mL r g / 2L l -mhXg / 2L l + mhYg / 2L t
... (23)
W 3 = mL r g / 2L l + mhXg / 2L l -mhYg / 2L t
... (24)
W 4 = mL r g / 2L l + mhXg / 2L l + mhYg / 2L t
... (25)
Here, in the above equations (22) to (25), m is the mass of the vehicle [kg], h is the height of the center of gravity of the vehicle [m], and g is the acceleration of gravity [m / s 2 ]. A method for calculating the wheel loads W 1 , W 2 , W 3 , and W 4 of each of the
ステップ40では、後輪3,4の舵角δ3,δ4がステアリング5の回転角θに対して次の式(26)の応答となるようにステアリングアクチュエータ16を制御する。
In
δ3=δ4=(1/16)[k0/(1+Tes)
−(Kf/Kr)Tes/(1+Tes)]θ
…(26)
ここで、式(26)のTe、k0は次の値である。
δ 3 = δ 4 = (1/16) [k 0 / (1 + T e s)
- (K f / K r) T e s / (1 + T e s)] θ
... (26)
Here, T e, k 0 of the formula (26) is the next value.
Te=IV/(2LlLfKf+mLrV2) …(26a)
k0=−{Lr+(mLf/2LlKr)KfV2}
/{Lf+(mLr/2LlKf)KrV2}
…(26b)
なお、上記の式(26)のmは車両の質量[kg]、Iは車体重心周りのヨー慣性モーメント[kgm2]であり、Kf,Krは前輪1,2及び後輪3,4の横滑り角が十分小さいときの単位横滑り角あたりのコーナーリングフォース[N/rad]である。また、式(26)右辺の左端の1/16はステアリング5の回転角θの変化に対する前輪舵角の感度であることは既に述べた通りである。
T e = IV / (2L l L f K f + mL r V 2) ... (26a)
k 0 = - {L r + (mL f / 2L l K r) K f V 2}
/ {L f + (mL r / 2L l K f ) K r V 2 }
... (26b)
In the above equation (26), m is the mass of the vehicle [kg], I is the yaw inertia moment [kgm 2 ] around the center of gravity of the vehicle body, and K f and K r are the
またこのとき、舵角センサ41〜44で各輪1〜4の舵角δ1,δ2,δ3,δ4[rad]を検出する。各輪1〜4の舵角δ1,δ2,δ3,δ4は、前述したように各輪1〜4の向きが車体前後方向と一致している状態をゼロとし、車両を鉛直上方からみたときに反時計回りを正とする。
At this time, the
このように後輪3,4の舵角δ3,δ4の目標応答をステアリング5の回転角θに対して決定することによって、車体速度Vの変化や各輪1〜4の舵角δ1,δ2,δ3,δ4が十分小さく、左右輪の制駆動力差がゼロの場合に、車体横滑り角βをゼロにできることが知られている(安部正人著,「自動車の運動と制御」,第二版,株式会社山海堂,平成15年4月,第8章8.3.1節)。
In this way, by determining the target responses of the steering angles δ 3 and δ 4 of the
ステップ50では、車両に対する運転者の要求制駆動力tF’[N]を次の式(27)により算出する。
In
tF’=tFa+tFb …(27)
式(27)右辺第1項のtFa[N]はアクセルペダル6の踏込量AP及び車体速度Vに対応した要求駆動力を要求駆動力マップに基づいて、また右辺第2項のtFb[N]はブレーキペダル7の踏込量BPに対応した要求制動力(ゼロ及び負の値)を要求制動力マップに基づいてそれぞれ算出したものである。ここで、要求駆動力マップはアクセルペダル6の踏込量AP及び車体速度Vに対応する要求駆動力を、要求制動力マップはブレーキペダルの踏込量BPに対応する要求制動力をそれぞれコントローラ8のROMに予め記憶させておいたものであり、要求駆動力tFa、要求制動力tFbはそれぞれ例えば図10及び図11のように設定している。
tF ′ = tFa + tFb (27)
Expression (27) tFa [N] in the first term on the right side is the required driving force corresponding to the depression amount AP of the
このようにして設定した要求制駆動力tF’に対して、次の式(27a)によりなまし処理を行い、なまし処理後の値を目標制駆動力tF[N]とする。 The required braking / driving force tF ′ set in this way is subjected to an annealing process according to the following equation (27a), and the value after the annealing process is set as a target braking / driving force tF [N].
tF={Ns×tF’(k)+(1−Ns)×tF’(k-1)}÷2 …(27a)
実際には、式(27a)におけるなまし処理の係数Ns(ただし、0<Ns<1)を用いて、式(27)で求めた要求制駆動力の現在値であるtF’(k)と、同じく式(27)で求めた要求制駆動力の前回値であるtF’(k-1)の平均値を演算する。なまし処理の係数Nsは制駆動力変動による加減速が運転者の不快感を招かないように本実施形態ではNsとして0.5を設定している。
tF = {Ns × tF ′ (k) + (1−Ns) × tF ′ (k−1)} / 2 (27a)
Actually, tF ′ (k), which is the current value of the requested braking / driving force obtained by the equation (27), using the smoothing coefficient Ns (where 0 <Ns <1) in the equation (27a), Similarly, the average value of tF ′ (k−1), which is the previous value of the required braking / driving force obtained by the equation (27), is calculated. In this embodiment, the smoothing coefficient Ns is set to 0.5 so that acceleration / deceleration due to fluctuations in braking / driving force does not cause discomfort to the driver.
上記の要求制駆動力tF’,目標制駆動力tF,要求駆動力tFa,要求制動力tFbはいずれも車両を前方に加速させる向きを正とする。要求制動力tFbはゼロまたは負の値である(図11参照)。 The required braking / driving force tF ′, target braking / driving force tF, required driving force tFa, and required braking force tFb are all positive in the direction in which the vehicle is accelerated forward. The required braking force tFb is zero or a negative value (see FIG. 11).
ステップ60では、まずステアリング5の回転角θと車両速度Vから、基本左右輪制駆動力差ΔF’[N]を基本左右輪制駆動力差マップに基づいて算出する。基本左右輪制駆動力差マップは、ステアリング回転角θと車両速度Vに対応する左右輪制駆動力差ΔF’を予めコントローラ8のROMに記憶しておいたものであり、基本左右輪制駆動力差ΔF’は例えば図12のように設定している。
In
そして、この基本左右輪制駆動力差ΔF’に対して次の式(27b)によりなまし処理を行うことによって目標左右輪制駆動力差ΔF[N]を求める。 Then, the target left / right wheel braking / driving force difference ΔF [N] is obtained by performing an annealing process on the basic left / right wheel braking / driving force difference ΔF ′ by the following equation (27b).
ΔF={Nss×ΔF’(k)+(1−Nss)×ΔF’(k-1)}÷2
…(27b)
実際には上記の式(27a)と同様に、式(27b)におけるなまし処理の係数Nss(但し、0<Nss<1)を用いて、基本左右輪制駆動力差の現在値であるΔF’(k)と、基本左右輪駆動力差の前回値であるΔF’(k-1)の平均値を演算する。なまし処理の係数Nssは左右輪制駆動力差の変動によるヨーレートや横Gの変化が運転者の不快感を招かないように本実施形態ではNssとして0.5を設定している。
ΔF = {Nss × ΔF ′ (k) + (1−Nss) × ΔF ′ (k−1)} / 2
... (27b)
Actually, like the above equation (27a), ΔF which is the current value of the basic left and right wheel braking / driving force difference using the smoothing coefficient Nss in equation (27b) (where 0 <Nss <1). The average value of '(k) and ΔF' (k-1) which is the previous value of the basic left and right wheel driving force difference is calculated. In this embodiment, the smoothing coefficient Nss is set to 0.5 so that the change in the yaw rate or the lateral G caused by the change in the left / right wheel braking / driving force difference does not cause the driver to feel uncomfortable.
ステップ70では、目標制駆動力tFと目標左右輪制駆動力差ΔFとから各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4を次の式(28)及び式(29)により算出する。各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4は、前述したように車両を前進させる方向に働く力を正とする。
In step 70, the braking / driving forces Fx 1 , Fx 2 , Fx 3 , and Fx 4 of the
Fx1=Fx3=tF/4−ΔF/4 …(28)
Fx2=Fx4=tF/4+ΔF/4 …(29)
例えば、ステアリング5の回転角θが正、つまりステアリング5を中立位置より図2のように反時計方向に回したとき、基本左右輪駆動力差ΔF’(従って目標左右輪駆動力差ΔF)が正の値となり(図12参照)、式(28)あるいは式(29)より左車輪1,3の制駆動力Fx1,Fx3が右車輪2,4の制駆動力Fx2,Fx4より小さくなる。この状態が図2に示されている。
Fx 1 = Fx 3 = tF / 4−ΔF / 4 (28)
Fx 2 = Fx 4 = tF / 4 + ΔF / 4 (29)
For example, when the rotation angle θ of the
ここで、式(28)及び式(29)の4つの制駆動力からなる制駆動力配分{Fx1,Fx2,Fx3,Fx4}が、要求された車両挙動(目標となる車両挙動)を実現する四輪の制駆動力配分である。 Here, the braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } composed of the four braking / driving forces of Expression (28) and Expression (29) is the required vehicle behavior (target vehicle behavior). ) Is a four-wheel braking / driving force distribution.
ステップ80では、各輪1〜4の横滑り角β1,β2,β3,β4[rad]を推定する。この各輪1〜4の横滑り角(スリップ角ともいう)とは、車両の進行方向とタイヤの前後方向のなす現時点における車輪のスリップ角(現状スリップ角)のことである。推定方法については種々あるが、ここでは一例として次の方法を用いる(特開平10-329689号公報参照)。
In
即ち、ステップ10で読み込んだ横方向加速度Yg及びヨーレートγと、ステップ20で求めている車体速度Vとから車体横滑り角βをまず推定する。その上で、この車体横滑り角β,ヨーレートγ,車体速度V,操舵角θから各輪1〜4の横滑り角β1,β2,β3,β4を推定する。具体的には先ず、車体横滑り角β[rad]を次の式(補4)により算出する。
That is, the vehicle body side slip angle β is first estimated from the lateral acceleration Yg and yaw rate γ read in
β=∫(Yg/V−γ)dt …(補4)
次いで、前輪1,2の横滑り角β1,β2を次の式(補5)により及び後輪3,4の横滑り角β3,β4を次の式(補6)によりそれぞれ算出する。
β = ∫ (Yg / V−γ) dt (Supplement 4)
Next, the sideslip angles β 1 and β 2 of the
β1=β2=β+θ/Gs−γ×Lf/V …(補5)
β3=β4=β+γ×Lr/V …(補6)
ここで、Gsはステアリングギヤ15のギヤ比である。
β 1 = β 2 = β + θ / Gs−γ × L f / V (Supplement 5)
β 3 = β 4 = β + γ × L r / V (Supplement 6)
Here, Gs is a gear ratio of the
各輪1〜4の横滑り角β1,β2,β3,β4の符号は、各輪1〜4の前後方向から車輪速度の方向までの角度が鉛直上方から見て反時計回りになっている場合を正とする。車体横滑り角βの符号も、各輪1〜4の横滑り角β1,β2,β3,β4の符号と同様である。
The signs of the sideslip angles β 1 , β 2 , β 3 , β 4 of the
ステップ90では、各車輪1〜4の路面摩擦係数μ1,μ2,μ3,μ4を推定する。推定方法については種々あるが、ここでは一例として次の方法(特開平6−98418号公報参照)を用いる。この方法は、各車輪1〜4が路面から受ける路面反力F1,F2,F3,F4を推定し、この路面反力F1,F2,F3,F4とステップ20で求めた輪荷重W1,W2,W3,W4とから各車輪1〜4の路面摩擦係数μ1,μ2,μ3,μ4を推定するものである。
In
簡単に概説すると、一輪についてモータに電磁トルクTmが加えられ、車輪には路面反力Fに車輪半径Rを乗じた路面反力トルクF×Rがモータによるトルクと逆方向に加えられているとする。そして、モータと車輪とは直結状態であり、車軸のねじり剛性κが十分に大きいと仮定でき(車軸のねじり変形を無視する)、モータの回転速度と車輪の回転速度とは同一の回転速度ω[rad/s]であるとの関係が成り立つとすると、モータと車輪との回転系の運動方程式は次の式(補7)にまとめることができる。 Briefly, the electromagnetic torque Tm is applied to the motor for one wheel, and the road surface reaction force torque F × R, which is the road surface reaction force F multiplied by the wheel radius R, is applied to the wheel in the opposite direction to the motor torque. To do. And it can be assumed that the motor and the wheel are directly connected, and the torsional rigidity κ of the axle is sufficiently large (ignoring torsional deformation of the axle), and the rotational speed of the motor and the rotational speed of the wheel are the same rotational speed ω. Assuming that the relationship of [rad / s] holds, the equation of motion of the rotating system of the motor and wheels can be summarized as the following equation (Supplement 7).
(Jm+Jw)dω/dt=Tm−Cm・ω−Rm・ω−F・R …(補7)
ここで、Jmはモータの慣性モーメント、Jwは車輪の慣性モーメント、Cmはモータの回転系の粘性減衰定数、Cwは車輪の回転系の粘性減衰定数、Rmはモータの回転系の個体摩擦、Rwは車輪の回転系の個体摩擦である。
(Jm + Jw) dω / dt = Tm−Cm · ω−Rm · ω−F · R (Supplement 7)
Here, Jm is the moment of inertia of the motor, Jw is the moment of inertia of the wheel, Cm is the viscous damping constant of the rotating system of the motor, Cw is the viscous damping constant of the rotating system of the wheel, Rm is the individual friction of the rotating system of the motor, Rw Is the individual friction of the rotating system of the wheel.
その結果、路面反力Fを、上記の式(補7)を変形した次の式(補8)より算出することができる。 As a result, the road surface reaction force F can be calculated from the following equation (complement 8) obtained by modifying the above equation (complement 7).
F={Tm−(Jm+Jw)dω/dt−Cm・ω−Rm・ω}/R
…(補8)
従って、この式(補8)を用いて、各輪1〜4についての路面反力F1,F2,F3,F4を次の式(補9a)〜式(補9d)により算出する。
F = {Tm− (Jm + Jw) dω / dt−Cm · ω−Rm · ω} / R
... (Supplement 8)
Therefore, the road surface reaction forces F 1 , F 2 , F 3 , and F 4 for the
F1={Tm−(Jm+Jw)dω1/dt−Cm・ω1−Rm・ω1}/R
…(補9a)
F2={Tm−(Jm+Jw)dω2/dt−Cm・ω2−Rm・ω2}/R
…(補9b)
F3={Tm−(Jm+Jw)dω3/dt−Cm・ω3−Rm・ω3}/R
…(補9c)
F4={Tm−(Jm+Jw)dω4/dt−Cm・ω4−Rm・ω4}/R
…(補9d)
このようにして算出した各輪1〜4の路面反力F1,F2,F3,F4と、各輪1〜4の輪荷重W1,W2,W3,W4とにより、次の式(補10a)〜式(補10d)を用いて路面摩擦係数μ1、μ2、μ3、μ4を算出する。
F 1 = {Tm− (Jm + Jw) dω 1 / dt−Cm · ω 1 −Rm · ω 1 } / R
... (Supplement 9a)
F 2 = {Tm− (Jm + Jw) dω 2 / dt−Cm · ω 2 −Rm · ω 2 } / R
... (Supplement 9b)
F 3 = {Tm− (Jm + Jw) dω 3 / dt−Cm · ω 3 −Rm · ω 3 } / R
... (Supplement 9c)
F 4 = {Tm− (Jm + Jw) dω 4 / dt−Cm · ω 4 −Rm · ω 4 } / R
... (Supplement 9d)
The road surface reaction forces F 1 , F 2 , F 3 , F 4 of the wheels 1 to 4 calculated in this way and the wheel loads W 1 , W 2 , W 3 , W 4 of the wheels 1 to 4 are The road surface friction coefficients μ 1 , μ 2 , μ 3 , and μ 4 are calculated using the following equations (complement 10a) to equation (complement 10d).
μ1=F1/W1 …(補10a)
μ2=F2/W2 …(補10b)
μ3=F3/W3 …(補10c)
μ4=F4/W4 …(補10d)
ステップ100では、ステップ70で求められた制駆動力配分{Fx1,Fx2,Fx3,Fx4}により実現される車両挙動(タイヤ力の総和の車体前後方向成分Fxと、タイヤ力の総和の車体横方向成分Fyと、各輪のタイヤ力によって発生する車体重心周りのヨーモーメントMの総和との3つ)を求める。
μ 1 = F 1 / W 1 (Supplement 10a)
μ 2 = F 2 / W 2 (Supplement 10b)
μ 3 = F 3 / W 3 (Supplement 10c)
μ 4 = F 4 / W 4 (Supplement 10d)
In
具体的にはステップ70で算出されている各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4と、各輪1〜4のタイヤ横力Fy1,Fy2,Fy3,Fy4と、ステップ40で算出されている各輪1〜4の舵角δ1,δ2,δ3,δ4とから、次式により各輪1〜4で発生するタイヤ力の車体前後方向成分Fx’1,Fx’2,Fx’3,Fx’4と、各輪1〜4で発生するタイヤ力の車体横方向成分Fy’1,Fy’2,Fy’3,Fy’4とを算出する。
Specifically, the braking / driving forces Fx 1 , Fx 2 , Fx 3 , Fx 4 of the
Fx’1=Fx1cosδ1−Fy1sinδ1 …(補11a)
Fx’2=Fx2cosδ2−Fy2sinδ2 …(補11b)
Fx’3=Fx3cosδ3−Fy3sinδ3 …(補11c)
Fx’4=Fx4cosδ4−Fy4sinδ4 …(補11d)
Fy’1=Fx1sinδ1+Fy1cosδ1 …(補12a)
Fy’2=Fx2sinδ2+Fy2cosδ2 …(補12b)
Fy’3=Fx3sinδ3+Fy3cosδ3 …(補12c)
Fy’4=Fx4sinδ4+Fy4cosδ4 …(補12d)
ここで、式(補11a)〜式(補12d)は前述の式(1),式(2)と同じ式である。
Fx ′ 1 = Fx 1 cos δ 1 −Fy 1 sin δ 1 (Supplement 11a)
Fx ′ 2 = Fx 2 cos δ 2 −Fy 2 sin δ 2 (Supplement 11b)
Fx ′ 3 = Fx 3 cosδ 3 −Fy 3 sinδ 3 (Supplement 11c)
Fx ′ 4 = Fx 4 cosδ 4 −Fy 4 sinδ 4 (Supplement 11d)
Fy ′ 1 = Fx 1 sin δ 1 + Fy 1 cos δ 1 (Supplement 12a)
Fy ′ 2 = Fx 2 sin δ 2 + Fy 2 cos δ 2 (Supplement 12b)
Fy ′ 3 = Fx 3 sinδ 3 + Fy 3 cosδ 3 (Supplement 12c)
Fy ′ 4 = Fx 4 sinδ 4 + Fy 4 cosδ 4 (Supplement 12d)
Here, the expressions (complement 11a) to (complement 12d) are the same expressions as the expressions (1) and (2) described above.
そして、式(補11a)〜式(補12d)により算出したタイヤ力の車体前後方向成分Fx’1,Fx’2,Fx’3,Fx’4,タイヤ力の車体横方向成分Fy’1,Fy’2,Fy’3,Fy’4を次の式(補13)〜式(補15)に代入して、タイヤ力の総和の車体前後方向成分Fx、タイヤ力の総和の車体横方向成分Fy、各輪1〜4のタイヤ力によって発生する車両重心周りのヨーモーメントの総和Mをそれぞれ算出する。
Then, the vehicle body longitudinal components Fx ′ 1 , Fx ′ 2 , Fx ′ 3 , Fx ′ 4 of the tire force calculated by the equations (complement 11a) to (complement 12d), the vehicle lateral component Fy ′ 1 of the tire force, Substituting Fy ′ 2 , Fy ′ 3 , and Fy ′ 4 into the following formulas (complement 13) to formula (complement 15), the vehicle body longitudinal component Fx of the total tire force, and the vehicle lateral component of the total tire force Fy, the sum M of yaw moments around the center of gravity of the vehicle generated by the tire force of each
Fx=Fx’1+Fx’2+Fx’3+Fx’4 …(補13)
Fy=Fy’1+Fy’2+Fy’3+Fy’4 …(補14)
M={(Fx’2+Fx’4)−(Fx’1+Fx’3)}×Lt/2
+{(Fy’1+Fy’2)×Lf−(Fy’3+Fy’4)×Lr}
…(補15)
ここで、式(補13)〜式(補15)は上記の式(8)〜式(10)と同じ式である。
Fx = Fx ′ 1 + Fx ′ 2 + Fx ′ 3 + Fx ′ 4 (Supplement 13)
Fy = Fy ′ 1 + Fy ′ 2 + Fy ′ 3 + Fy ′ 4 (Supplement 14)
M = {(Fx ′ 2 + Fx ′ 4 ) − (Fx ′ 1 + Fx ′ 3 )} × L t / 2.
+ {(Fy ′ 1 + Fy ′ 2 ) × L f − (Fy ′ 3 + Fy ′ 4 ) × L r }
... (Supplement 15)
Here, the expressions (complement 13) to the expression (complement 15) are the same as the expressions (8) to (10).
上記の式(補11a)〜式(補12d)のタイヤ横力Fy1,Fy2,Fy3,Fy4は次のようにして求める。例えば、左前輪1の場合を例にとって説明すると、コントローラ8のROMに、左前輪1の輪荷重W1,横滑り角β1毎に、左前輪1の制駆動力とタイヤ横力との関係を表すタイヤ特性マップを複数用意しておく。この複数のタイヤ特性は、図14のように設定されている。すなわち、図14において上段に左前輪1の輪荷重W1が大きい場合、中段に左前輪1の輪荷重W1が中位の場合、下段に左前輪1の輪荷重W1が小さい場合の各タイヤ特性の内容を示しており、左前輪1の輪荷重W1が大きくなるほど、つまり下段、中段、上段の順にタイヤ特性を表す曲線が全体的に膨らんでいる。また、上段、中段、下段の3つの各場合ではさらに左前輪1の代表的な6つの横滑り角α1,α2,α3,α4,α5,α6(|α1|<|α2|<|α3|<|α4|<|α5|<|α6|)のときのタイヤ特性を重ねて示しており、同じ輪荷重W1でも横滑り角の絶対値が|α1|より|α4|へと大きくなるほど上方への膨らみが増した後に、|α5|、|α6|となれば今度は全体的に縮んでいる。従って、今かりに、左前輪1の輪荷重W1が大きく、左前輪1の横滑り角β1がα2であったとすれば、図14に示す複数のタイヤ特性のうち上段のα2のタイヤ特性を選択し、その選択したタイヤ特性から左前輪1の制駆動力Fx1に対応するタイヤ横力Fy1を求めることができる。
The tire lateral forces Fy 1 , Fy 2 , Fy 3 and Fy 4 in the above equations (complement 11a) to equation (complement 12d) are obtained as follows. For example, the case of the left
残りの三輪2〜4についても同様に、各輪2〜4の輪荷重Wi,横滑り角βi(i=2〜4)に対応した複数のタイヤ特性を図14と同様に用意しておき、各輪2〜4の制駆動力Fxiに対応するタイヤ横力Fyiを、各輪2〜4毎の複数のタイヤ特性のうち輪荷重Wi、横滑り角βiのときのタイヤ特性を選択して求める。
Similarly, for the remaining three
ステップ101では、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動(Fx,Fy,M)を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}(j=1,2,…)を求める。ここで、{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}は、{Fx1,Fx2,Fx3,Fx4}={Fx1(1),Fx2(1),Fx3(1),Fx4(1)}、{Fx1,Fx2,Fx3,Fx4}={Fx1(2),Fx2(2),Fx3(2),Fx4(2)}、{Fx1,Fx2,Fx3,Fx4}={Fx1(3),Fx2(3),Fx3(3),Fx4(3)}、{Fx1,Fx2,Fx3,Fx4}={Fx1(4),Fx2(4),Fx3(4),Fx4(4)}、…といった制駆動力配分{Fx1,Fx2,Fx3,Fx4}の集合を表している。
In
ここで、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動(Fx,Fy,M)を実現する制駆動力配分の集合の演算方法として次に2つを挙げる。
<1>制駆動力配分の集合の演算方法1(請求項2に記載の発明)
現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を、ステップ50で求めている目標制駆動力tFと、ステップ60で求めている目標左右輪制駆動力差ΔFと、ステップ30で求めている各輪1〜4の輪荷重W1,W2,W3,W4と、ステップ80で求めている各輪1〜4の横滑り角β1,β2,β3,β4とから、図13に示すフローで求めたマップに基づいて求める。
Here, as a calculation method of a set of braking / driving force distributions that realizes the same vehicle behavior (Fx, Fy, M) as obtained in
<1>
A set {Fx 1 (j), Fx 2 (j), Fx 3 (Fx 1 (j), Fx 3 (j), Fx 3 (j), Fx 3 (j), Fx 3 (j), Fx 3 (j)] j), Fx 4 (j)}, the target braking / driving force tF obtained in
ここで、図13に示すフローで求めたマップとは、目標制駆動力tFと、目標左右輪制駆動力差ΔFと、輪荷重W1,W2,W3,W4と、横滑り角β1,β2,β3,β4とから、図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合を求めるマップである。なお、図13のフローによる制駆動力配分の集合の求め方については後述する。
Here, the map obtained by the flow shown in FIG. 13 includes the target braking / driving force tF, the target left / right wheel braking / driving force difference ΔF, the wheel loads W 1 , W 2 , W 3 , W 4, and the
ここではさらに、図14に示すような各輪1〜4の制駆動力とタイヤ横力との関係を表すタイヤ特性を、ステップ90で推定した各輪1〜4の路面摩擦係数μ1,μ2,μ3,μ4毎に求めてマップにしておくことにすれば、制駆動力配分の集合をより高精度に求めることができる。
<2>制駆動力配分の集合の算出方法2(請求項3、4に記載の発明)
現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を、上記の式(20)または式(19)を用いた図15に示すフローに従って求める。図15に示すフローを用いた制駆動力配分の集合の求め方については後述する。
Here, the tire characteristics representing the relationship between the braking / driving force and the tire lateral force of each of the
<2>
A set {Fx 1 (j), Fx 2 (j), Fx 3 (Fx 1 (j), Fx 3 (j), Fx 3 (j), Fx 3 (j), Fx 3 (j), Fx 3 (j)] j), Fx 4 (j)} is obtained according to the flow shown in FIG. 15 using the above equation (20) or equation (19). A method of obtaining a set of braking / driving force distributions using the flow shown in FIG. 15 will be described later.
これで制駆動力配分の集合の算出方法の説明を終える。 This concludes the description of the method for calculating the braking / driving force distribution set.
ステップ110では、各輪1〜4の制駆動力上限Fmaxi[N](i=1〜4)及び各輪1〜4の制駆動力下限Fmini[N](i=1〜4)を求める。
In step 110, the braking / driving force upper limit Fmax i [N] (i = 1 to 4) of each
この各輪1〜4の制駆動力上限Fmaxi及び制駆動力下限Fminiの算出方法として次に3つを挙げる。
<3>制駆動力上限、制駆動力下限の算出方法1(請求項5、6に記載の発明)
各輪1〜4において、スリップあるいは車輪ロックを起こさない各輪1〜4の制駆動力上限Fsmaxi[N](i=1〜4)及び制駆動力下限Fsmini[N](i=1〜4)を求め、このうちスリップあるいは車輪ロックを起こさない各輪1〜4の制駆動力上限Fsmaxiを制駆動力上限Fmaxiに、スリップあるいは車輪ロックを起こさない各輪1〜4の制駆動力下限Fsminiを制駆動力下限Fminiにする。
Three methods for calculating the braking / driving force upper limit Fmax i and the braking / driving force lower limit Fmin i of each of the
<3>
In each of the
ここで、スリップあるいは車輪ロックを起こさない各輪1〜4の制駆動力上限Fsmaxi及び制駆動力下限Fsminiの求め方としては、例えば、上記の式(補9a)〜式(補9d)により算出している各輪1〜4の路面反力Fi(i=1〜4)を、スリップあるいは車輪ロックを起こさない各輪1〜4の制駆動力上限Fsmaxiと、また−Fiを、スリップあるいは車輪ロックを起こさない各輪1〜4の制駆動力下限Fsminiとする。
<4>制駆動力上限、制駆動力下限の算出方法2(請求項7、8に記載の発明)
各モータ11〜14が過熱しまたは破損しないようにする各輪1〜4の制駆動力上限Fdmaxi[N](i=1〜4)及び制駆動力下限Fdmini[N](i=1〜4)を求め、このうち各モータ11〜14が過熱しまたは破損しないようにする各輪1〜4の制駆動力上限Fdmaxiを制駆動力上限Fmaxiに、各モータ11〜14が過熱しまたは破損しないようにする各輪1〜4の制駆動力下限Fdminiを制駆動力下限Fminiにする。
Here, as a method of obtaining the braking / driving force upper limit Fsmax i and the braking / driving force lower limit Fsmin i of each of the
<4>
Longitudinal force upper limit fdmax i of each wheel 1-4
ここで、各モータ11〜14が過熱しまたは破損しないようにする各輪1〜4の制駆動力上限Fdmaxi及び制駆動力下限Fdminiの求め方としては、例えば、各モータ11〜14の温度Ti(i=1〜4)から、現在のモータ温度とモータ過熱を抑えることができる最大出力Ptmax[W]との関係を求めたマップに基づいて、各モータ11〜14の最大出力Ptmaxi(i=1〜4)[W]を求め、この各モータ11〜14の最大出力Ptmaxiと、ステップ10で検出している各輪1〜4の回転速度ωiとから次の式(30)及び式(31)により、各モータ11〜14が過熱しまたは破損しないようにする各輪1〜4の制駆動力上限Fdmaxi及び制駆動力下限Fdminiを各輪1〜4毎に求める。
Here, as a method of obtaining the braking / driving force upper limit Fdmax i and the braking / driving force lower limit Fdmin i of each of the
Fdmaxi=Ptmaxi÷ωi …(30)
Fdmini=−Ptmaxi÷ωi …(31)
現在のモータ温度とモータ過熱を抑えることができる最大出力Ptmaxとの関係を求めたマップでは例えば図16のようにモータ温度とPtmaxとの関係を設定している。
Fdmax i = Ptmax i ÷ ω i (30)
Fdmin i = −Ptmax i ÷ ω i (31)
In the map for determining the relationship between the current motor temperature and the maximum output Ptmax that can suppress the motor overheating, the relationship between the motor temperature and Ptmax is set as shown in FIG. 16, for example.
ここで、各輪1〜4のメカブレーキによる制動力と各モータ11〜14の制駆動力を協調制御できる四輪独立駆動車であれば、各モータ11〜14が過熱しまたは破損しないようにする各輪1〜4の制駆動力下限Fdminiに、各輪1〜4のメカブレーキの最大制動力を加算することによって各輪1〜4の制動性能の限界を向上させることができる。
<5>制駆動力上限、制駆動力下限の算出方法3
上記<4>で算出した制駆動力上限Fdmaxiと、上記<3>で算出した制駆動力上限Fsmaxiとを比較して小さい方の値を駆動力上限Fmaxiに設定し、同様に上記<4>で算出した制駆動力下限Fdminiと、上記<3>で算出した制駆動力下限Fsminiとを比較して大きい方の値を制駆動力下限Fminiに設定する。この方法によれば、路面に伝達可能な駆動力制限と、モータ過熱による出力低下とを同時に考慮できるので、なおよい。
Here, if the four-wheel independent drive vehicle is capable of cooperatively controlling the braking force by the mechanical brakes of the
<5>
The braking / driving force upper limit Fdmax i calculated in <4> is compared with the braking / driving force upper limit Fsmax i calculated in <3>, and the smaller value is set as the driving force upper limit Fmax i. a longitudinal force limit Fdmin i calculated in <4>, it sets the larger value by comparing the longitudinal force limit Fsmin i calculated for the longitudinal force limit Fmin i above <3>. According to this method, it is even better because it is possible to simultaneously consider the driving force limitation that can be transmitted to the road surface and the output reduction due to motor overheating.
これで制駆動力上限、制駆動力下限の算出方法の説明を終える。 This concludes the description of the method for calculating the braking / driving force upper limit and the braking / driving force lower limit.
図9(B)に移り、ステップ120では、必要制駆動力補正量ΔFsiの絶対値|ΔFsi|が、閾値Fthより大きい車輪が一輪以上あるか否かをみる。必要制駆動力補正量ΔFsiの絶対値|ΔFsi|が、閾値Fthより大きい車輪が一輪以上ある場合にはステップ130に進む。そうでない場合にはステップ300に進む。 Turning to FIG. 9B, in step 120, it is determined whether or not there are one or more wheels whose absolute value | ΔFs i | of the required braking / driving force correction amount ΔFs i is larger than the threshold value Fth. If one or more wheels have an absolute value | ΔFs i | of the required braking / driving force correction amount ΔFs i greater than the threshold value Fth, the routine proceeds to step 130. Otherwise, go to step 300.
ここで、必要制駆動力補正量ΔFsi(i=1〜4)は各輪1〜4毎に次のように設定している。すなわち、各輪1〜4について制駆動力Fxiが制駆動力下限Fminiを超えかつ制駆動力上限Fmaxi未満(Fmini<Fxi<Fmaxi)の領域にあれば必要制駆動力補正量ΔFsi[N]=0である。各輪1〜4について制駆動力Fxiが制駆動力上限Fmaxi以上(Fxi≧Fmaxi)のとき必要制駆動力補正量ΔFsi[N]=Fmaxi−Fxiであり、各輪1〜4について制駆動力Fxiが制駆動力下限Fmini以下(Fxi≦Fmini)のとき必要制駆動力補正量ΔFsi=Fmini−Fxiである。これより、各輪1〜4の必要制駆動力補正値ΔFsiとは、各輪1〜4の制駆動力Fxiを各輪1〜4が出力可能な制駆動力にするために最低限どれだけ補正しなければならないかを表す値である。
Here, the necessary braking / driving force correction amount ΔFs i (i = 1 to 4) is set for each
また、上記の閾値Fthは、例えば路面反力Fiと制駆動力Fxiとの差が大きくなる、つまりスリップの傾向が強くなっていることや、モータ11〜14の過熱が進行している等の現象を判断するための閾値である。本実施形態では車両重量W[N]の1%(つまり0.01W)を閾値Fthとしている。
In addition, the threshold value Fth is, for example, that the difference between the road surface reaction force F i and the braking / driving force Fx i is large, that is, the tendency of slip is strong, or the
ステップ130では、ステップ101で求めた制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}をステップ110で求めた制駆動力上限Fmaxi及び制駆動力下限Fminiで各輪1〜4毎に制限する。制限する方法について次に説明する。
In step 130, the braking / driving force distribution set {Fx 1 (j), Fx 2 (j), Fx 3 (j), Fx 4 (j)} obtained in
まず、左前輪1の制駆動力Fx1(j)がFmin1≦Fx1(j)≦Fmax1の条件を満たす整数jを全て抽出し、配列Stに記憶する。次に、配列Stに記憶された整数jの中から右前輪2の制駆動力Fx2(j)がFmin2≦Fx2(j)≦Fmax2の条件を満たす整数jをさらに抽出する。この作業をさらに、左後輪3、右後輪4の順番で行い、最終的に実現可能な制駆動力配分を指し示す整数jの集合を配列Stという形で得る。また、実現可能な制駆動力配分を指し示す整数jの数を配列数Stnumに設定する。
First, all integers j in which the braking / driving force Fx 1 (j) of the left
このように制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を制駆動力上限Fmaxi及び制駆動力下限Fminiで制限することによって、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合より、スリップやモータ過熱等の制駆動力制限によって実現不可能な制駆動力配分{Fx1,Fx2,Fx3,Fx4}を除外することができる。
As described above, the set {Fx 1 (j), Fx 2 (j), Fx 3 (j), Fx 4 (j)} of the braking / driving force distribution is limited by the braking / driving force upper limit Fmax i and the braking / driving force lower limit Fmin i . As a result, the braking / driving force limit such as slip and motor overheating is determined from the set of braking / driving force distributions that realizes the same vehicle behavior as that obtained in
ステップ140では、配列数Stnumがゼロでない場合、つまりステップ100で求めたと同一の車両挙動を実現する制駆動力配分{Fx1,Fx2,Fx3,Fx4}が存在する場合にはステップ150に進み、ステップ130で得られた実現可能な制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}の中から、一つの制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択する(制駆動力の再配分を行う)。
In step 140, if the sequence number Stnum is not zero, i.e. the braking-driving force distribution to realize the same vehicle behavior and calculated in step 100 {Fx 1, Fx 2, Fx 3, Fx 4}
ステップ150におけるこの一つの制駆動力配分の選択方法として次に3つを挙げる。
<6>一つの制駆動力配分の選択1
配列Stに記録された各整数jに対応する制駆動力配分と、ステップ70で設定している各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4とから次の式(32)の評価関数J(j)を配列Stに記録された整数j毎に計算し、それら計算した評価関数J(j)のうち、評価関数J(j)の値を最小とする制駆動力配分を選ぶ。
Three methods for selecting this one braking / driving force distribution in
<6> Selection of one braking / driving
From the braking / driving force distribution corresponding to each integer j recorded in the array St and the braking / driving forces Fx 1 , Fx 2 , Fx 3 , Fx 4 of the wheels 1 to 4 set in Step 70, The evaluation function J (j) of (32) is calculated for each integer j recorded in the array St, and the braking / driving that minimizes the value of the evaluation function J (j) among the calculated evaluation functions J (j). Choose power distribution.
J(j)=(Fx1−Fx1(j))+(Fx2−Fx2(j))
+(Fx3−Fx3(j))+(Fx4−Fx4(j))
…(32)
例えば、配列Stに記録された整数jが1,3,5であったとすると、次の式(補16a)〜式(補16c)により評価関数J(1),J(3),J(5)を計算する。
J (j) = (Fx 1 −Fx 1 (j)) + (Fx 2 −Fx 2 (j))
+ (Fx 3 −Fx 3 (j)) + (Fx 4 −Fx 4 (j))
... (32)
For example, if the integer j recorded in the array St is 1, 3, and 5, the evaluation functions J (1), J (3), J (5) are obtained by the following equations (complement 16a) to (complement 16c). ).
J(1)=(Fx1−Fx1(1))+(Fx2−Fx2(1))
+(Fx3−Fx3(1))+(Fx4−Fx4(1))
…(補16a)
J(3)=(Fx1−Fx1(3))+(Fx2−Fx2(3))
+(Fx3−Fx3(3))+(Fx4−Fx4(3))
…(補16b)
J(5)=(Fx1−Fx1(5))+(Fx2−Fx2(5))
+(Fx3−Fx3(5))+(Fx4−Fx4(5))
…(補16c)
そして、J(1)>J(3)>J(5)であれば、{Fx1(5),Fx2(5),Fx3(5),Fx4(5)}を求めるべき制駆動力配分{Fx1,Fx2,Fx3,Fx4}として選択する。
J (1) = (Fx 1 −Fx 1 (1)) + (Fx 2 −Fx 2 (1))
+ (Fx 3 −Fx 3 (1)) + (Fx 4 −Fx 4 (1))
... (Supplement 16a)
J (3) = (Fx 1 −Fx 1 (3)) + (Fx 2 −Fx 2 (3))
+ (Fx 3 −Fx 3 (3)) + (Fx 4 −Fx 4 (3))
... (Supplement 16b)
J (5) = (Fx 1 −Fx 1 (5)) + (Fx 2 −Fx 2 (5))
+ (Fx 3 −Fx 3 (5)) + (Fx 4 −Fx 4 (5))
... (Supplement 16c)
And if J (1)> J (3)> J (5), the braking / driving to obtain {Fx 1 (5), Fx 2 (5), Fx 3 (5), Fx 4 (5)} Select force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 }.
このように、評価関数J(j)の値を最小とする制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択することは、スリップやモータ過熱等による制駆動力制限によって制駆動力を再配分する場合に、再配分しない場合との各輪1〜4の制駆動力差を最小にする方法の一つであり、制駆動力の急変によるショックを抑制する効果がある。
<7>一つの制駆動力配分の選択2(請求項9に記載の発明)
配列Stに記録された各整数jに対応する制駆動力配分の中でモータ11〜14の総消費電力が最小となる制駆動力配分を選択する。選択の方法としては以下の通り。
Thus, selecting the value minimum to braking-driving force distribution of the evaluation function J (j) {Fx 1, Fx 2,
<7> One braking / driving force distribution selection 2 (Invention of Claim 9)
The braking / driving force distribution that minimizes the total power consumption of the
モータ11〜14が消費する電力と、実際に制駆動力として得られる出力との差、つまり制駆動トルクと回転速度とをパラメータとして予め求めておいたモータロス(モータ11〜14の損失エネルギ)のマップから、配列Stに記録された各整数jに対応する制駆動力配分毎に、各輪1〜4のモータロスPlossi(j)[W](i=1〜4)を求め、これらモータロスPlossi(j)の和Ploss(j)[W]を次式(補17)により算出し、この算出した総モータロスPloss(j)が最小となる制駆動力配分を選択する。
The difference between the electric power consumed by the
Ploss(j)=Ploss1(j)+Ploss2(j)+Ploss3(j)+Ploss4(j)
…(補17)
各輪1〜4のモータロスは、図17のように制駆動トルク[Nm]と回転速度[rad/s]とに応じ、回転速度が同じであれば制駆動トルク(=制駆動力×R)が正の値であるいは負の値で増えるほどモータロスが大きくなるように、また制駆動トルクが同じであれば回転速度が高くなるほどモータロスが大きくなるように予め設定している。図17において、図示の5つのα1〜α5が値の異なるモータロスであり、α1よりα2のほうが、α2よりα3のほうが、α3よりα4のほうが、α4よりα5のほうがモータロスが大きくなっている(α5>α4>α3>α2>α1)。なお、図17でもαを用いているが、図14のαとは異なる。例えば、左前輪1の制駆動力Fx1と回転速度ω1から、左前輪1のモータロスPloss1(j)を求めるには次のようにする。すなわち、図17において横軸がω1、縦軸がFx1×RのときのモータロスPlossは所定値α5であるので、この所定値α5を左前輪1のモータロスPloss1(j)とする。
Ploss (j) = Ploss 1 (j) + Ploss 2 (j) + Ploss 3 (j) + Ploss 4 (j)
... (Supplement 17)
As shown in FIG. 17, the motor loss of each of the
具体的に考えてみる。配列Stに記録された整数jが1,3,5であったとすると、次の式(補18a)〜式(補18c)により総モータロスPloss(1),Ploss(3),Ploss(5)を計算する。 Think concretely. If the integer j recorded in the array St is 1, 3, and 5, the total motor loss Ploss (1), Ploss (3), Ploss (5) is calculated by the following equations (complement 18a) to (complement 18c). calculate.
Ploss(1)=Ploss1(1)+Ploss2(1)+Ploss3(1)+Ploss4(j)
…(補18a)
Ploss(3)=Ploss1(3)+Ploss2(3)+Ploss3(3)+Ploss4(3)
…(補18b)
Ploss(5)=Ploss1(5)+Ploss2(5)+Ploss3(5)+Ploss4(5)
…(補18c)
そして、Ploss(1)>Ploss(3)<Ploss(5)であれば、{Fx1(3),Fx2(3),Fx3(3),Fx4(3)}を求めるべき制駆動力配分{Fx1,Fx2,Fx3,Fx4}として選択する。
<8>一つの制駆動力配分の選択3(請求項10に記載の発明)
配列Stに記録された各整数jに対応する制駆動力配分の中で、各輪1〜4のタイヤ横力と制駆動力との合力であるタイヤ力を輪荷重W1,W2,W3,W4で除した値であるタイヤ摩擦円使用率H1(j),H2(j),H3(j),H4(j)の最大値と最小値の差が最も小さくなる制駆動力配分を選択する。選択の方法は次の通り。
Ploss (1) = Ploss 1 (1) + Ploss 2 (1) + Ploss 3 (1) + Ploss 4 (j)
... (Supplement 18a)
Ploss (3) = Ploss 1 (3) + Ploss 2 (3) + Ploss 3 (3) + Ploss 4 (3)
... (Supplement 18b)
Ploss (5) = Ploss 1 (5) + Ploss 2 (5) + Ploss 3 (5) + Ploss 4 (5)
... (Supplement 18c)
And if Ploss (1)> Ploss (3) <Ploss (5), {Fx 1 (3), Fx 2 (3), Fx 3 (3), Fx 4 (3)} should be determined to be controlled. Select force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 }.
<8> Selection of one braking / driving force distribution 3 (Invention of Claim 10)
In the braking / driving force distribution corresponding to each integer j recorded in the array St, the tire force that is the resultant of the tire lateral force and the braking / driving force of each
配列Stに記録された各整数jに対応する制駆動力配分毎に、各輪1〜4のタイヤ横力Fy1(j),Fy2(j),Fy3(j),Fy4(j)を推定(算出)し、このタイヤ横力Fy1(j),Fy2(j),Fy3(j),Fy4(j)と、ステップ30で求めている輪荷重W1,W2,W3,W4とから各輪1〜4のタイヤ摩擦円使用率H1(j),H2(j),H3(j),H4(j)を次の式(33a)〜式(33d)により算出する。
For each braking / driving force distribution corresponding to each integer j recorded in the array St, tire lateral forces Fy 1 (j), Fy 2 (j), Fy 3 (j), Fy 4 (j ) Is estimated (calculated), and the tire lateral forces Fy 1 (j), Fy 2 (j), Fy 3 (j), Fy 4 (j) and the wheel loads W 1 , W 2 obtained in
H1(j)={(Fx1(j)2+Fy1(j)2)}1/2/W1 …(33a)
H2(j)={(Fx2(j)2+Fy2(j)2)}1/2/W2 …(33b)
H3(j)={(Fx3(j)2+Fy3(j)2)}1/2/W3 …(33c)
H4(j)={(Fx4(j)2+Fy4(j)2)}1/2/W4 …(33d)
ここで、各輪1〜4のタイヤ横力Fy1(j),Fy2(j),Fy3(j),Fy4(j)の算出方法の一例はステップ100で説明済みであるので、ここでの説明は省略する。
H 1 (j) = {(Fx 1 (j) 2 + Fy 1 (j) 2 )} 1/2 / W 1 (33a)
H 2 (j) = {(Fx 2 (j) 2 + Fy 2 (j) 2 )} 1/2 / W 2 (33b)
H 3 (j) = {(Fx 3 (j) 2 + Fy 3 (j) 2 )} 1/2 / W 3 (33c)
H 4 (j) = {(Fx 4 (j) 2 + Fy 4 (j) 2 )} 1/2 / W 4 (33d)
Here, since an example of the calculation method of the tire lateral forces Fy 1 (j), Fy 2 (j), Fy 3 (j), and Fy 4 (j) of each
そして配列Stに記録された整数jに対応する各制駆動力配分毎に、四輪1〜4のタイヤ摩擦円使用率H1(j),H2(j),H3(j),H4(j)のうちの最大値Hmax(j)と、タイヤ摩擦円使用率H1(j),H2(j),H3(j),H4(j)のうちの最小値Hmin(j)とを求め、その差であるΔH(j)=Hmax(j)−Hmin(j)が最も小さくなる整数jを採る制駆動力配分を選択する。
For each braking / driving force distribution corresponding to the integer j recorded in the array St, the tire friction circle usage rates H 1 (j), H 2 (j), H 3 (j), H for the four
ここでも、配列Stに記録された整数jが1,3,5であったとする場合で具体的に考えてみると、次の式(補19a)〜式(補21d)により各輪のタイヤ摩擦円使用率H1(1),H2(1),H3(1),H4(1)、H1(3),H2(3),H3(3),H4(3)、H1(5),H2(5),H3(5),H4(5)を計算する。 Here again, when the integer j recorded in the array St is 1, 3, and 5, specifically, tire friction of each wheel is calculated by the following equations (complement 19a) to equation (complement 21d). Yen usage rate H 1 (1), H 2 (1), H 3 (1), H 4 (1), H 1 (3), H 2 (3), H 3 (3), H 4 (3) , H 1 (5), H 2 (5), H 3 (5), and H 4 (5) are calculated.
H1(1)={(Fx1(1)2+Fy1(1)2)}1/2/W1 …(補19a)
H2(1)={(Fx2(1)2+Fy2(1)2)}1/2/W2 …(補19b)
H3(1)={(Fx3(1)2+Fy3(1)2)}1/2/W3 …(補19c)
H4(1)={(Fx4(1)2+Fy4(1)2)}1/2/W4 …(補19d)
H1(3)={(Fx1(3)2+Fy1(3)2)}1/2/W1 …(補20a)
H2(3)={(Fx2(3)2+Fy2(3)2)}1/2/W2 …(補20b)
H3(3)={(Fx3(3)2+Fy3(3)2)}1/2/W3 …(補20c)
H4(5)={(Fx4(3)2+Fy4(5)2)}1/2/W4 …(補20d)
H1(5)={(Fx1(5)2+Fy1(5)2)}1/2/W1 …(補21a)
H2(5)={(Fx2(5)2+Fy2(5)2)}1/2/W2 …(補21b)
H3(5)={(Fx3(5)2+Fy3(5)2)}1/2/W3 …(補21c)
H4(5)={(Fx4(5)2+Fy4(5)2)}1/2/W4 …(補21d)
そして、H1(1),H2(1),H3(1),H4(1)のうちの最大値と最小値の差ΔH(1)、H1(3),H2(3),H3(3),H4(3)のうちの最大値と最小値の差ΔH(3)、H1(5),H2(5),H3(5),H4(5)のうちの最大値と最小値の差ΔH(5)を求めたところ、ΔH(1)<ΔH(3)<ΔH(5)であったとすれば、{Fx1(1),Fx2(1),Fx3(1),Fx4(1)}を求めるべき制駆動力配分{Fx1,Fx2,Fx3,Fx4}として選択する。
H 1 (1) = {(Fx 1 (1) 2 + Fy 1 (1) 2 )} 1/2 / W 1 (Supplement 19a)
H 2 (1) = {(Fx 2 (1) 2 + Fy 2 (1) 2 )} 1/2 / W 2 (Supplement 19b)
H 3 (1) = {(Fx 3 (1) 2 + Fy 3 (1) 2 )} 1/2 / W 3 (Supplement 19c)
H 4 (1) = {(Fx 4 (1) 2 + Fy 4 (1) 2 )} 1/2 / W 4 (Supplement 19d)
H 1 (3) = {(Fx 1 (3) 2 + Fy 1 (3) 2 )} 1/2 / W 1 (Supplement 20a)
H 2 (3) = {(Fx 2 (3) 2 + Fy 2 (3) 2 )} 1/2 / W 2 (Supplement 20b)
H 3 (3) = {(Fx 3 (3) 2 + Fy 3 (3) 2 )} 1/2 / W 3 (Supplement 20c)
H 4 (5) = {(Fx 4 (3) 2 + Fy 4 (5) 2 )} 1/2 / W 4 (Supplement 20d)
H 1 (5) = {(Fx 1 (5) 2 + Fy 1 (5) 2 )} 1/2 / W 1 (Supplement 21a)
H 2 (5) = {(Fx 2 (5) 2 + Fy 2 (5) 2 )} 1/2 / W 2 (Supplement 21b)
H 3 (5) = {(Fx 3 (5) 2 + Fy 3 (5) 2 )} 1/2 / W 3 (Supplement 21c)
H 4 (5) = {(Fx 4 (5) 2 + Fy 4 (5) 2 )} 1/2 / W 4 (Supplement 21d)
Then, the difference ΔH (1), H 1 (3), H 2 (3) between the maximum value and the minimum value of H 1 (1), H 2 (1), H 3 (1), H 4 (1). ), H 3 (3), H 4 (3), the difference between the maximum and minimum values ΔH (3), H 1 (5), H 2 (5), H 3 (5), H 4 (5 ), The difference ΔH (5) between the maximum value and the minimum value is obtained. If ΔH (1) <ΔH (3) <ΔH (5), {Fx 1 (1), Fx 2 ( 1), Fx 3 (1), Fx 4 (1)} are selected as braking / driving force distributions {Fx 1 , Fx 2 , Fx 3 , Fx 4 } to be obtained.
このように一つの制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択することによって、各輪1〜4の路面摩擦係数が急激に低下しても、各輪1〜4のタイヤ摩擦円が飽和しにくくなり車両の操縦安定性をより向上させることができる。
<9>一つの制駆動力配分の選択4
車両の操縦安定性が重視される急旋回中には上記<8>で述べた選択方法を、これに対して燃費を重視する直進・緩旋回中は上記<6>や<7>で述べた選択方法を採ることもできる。
By selecting one braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } in this way, even if the road surface friction coefficient of each of the
<9> Selection of one braking / driving
The selection method described in <8> above during a sharp turn where vehicle handling stability is important, while the above <6> and <7> are described during straight and slow turns that emphasize fuel efficiency. A selection method can also be adopted.
選択方法の例としては、例えばステップ10で検出した横方向加速度Ygの絶対値が閾値Ygth(本実施形態では3.0 m/s2)以下のときに燃費を重視する直進・緩旋回中と判断し、上記の式(32)の評価関数J(j)を最小化する制駆動力配分を、これに対してステップ10で検出した横方向加速度Ygの絶対値が閾値Ygthより大きければ車両の操縦安定性が重視される急旋回中と判断し、配列Stに記録された整数jに対応する各制駆動力配分毎に、四輪1〜4のタイヤ摩擦円使用率H1(j),H2(j),H3(j),H4(j)の最大値Hmax(j)と最小値Hmin(j)を求め、その最大値と最小値の差であるΔH(j)=Hmax(j)−Hmin(j)が最も小さくなる制駆動力配分をそれぞれ選択する。
As an example of the selection method, for example, when the absolute value of the lateral acceleration Yg detected in
これで制駆動力の再配分の説明を終える。 This completes the explanation of the redistribution of braking / driving force.
一方、ステップ140で配列数Stnumがゼロである場合、つまり、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分が各輪1〜4の制駆動力制限により存在しない場合には、ステップ210〜270に進む。
On the other hand, if the arrangement number Stnum is zero in step 140, that is, the braking / driving force distribution for realizing the same vehicle behavior as that obtained in
ステップ210〜270は、必要制駆動力補正量の絶対値|ΔFsi|(i=1〜4)が最も大きくなっている車輪の制駆動力を出力可能な範囲に再設定しながら、車両挙動(前後方向の加速度Xg、横方向の加速度Yg、車両重心周りのヨーモーメントM)をなるべく乱さない各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を用いて、各輪の制駆動力Fx1,Fx2,Fx3,Fx4を求める部分である。言い換えると、ステップ210,220,230,240,250,260,270は、制駆動力上限Fmaxi(FdmaxiまたはFsmaxi)、制駆動力下限Fmini(FdminiまたはFsmini)によって制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}が制限され、ステップ70で定まっている制駆動力配分{Fx1,Fx2,Fx3,Fx4}を連続的に変化させることができない場合に、各輪1〜4の制駆動力変化が不連続となることを許可し、ステップ70で定まっている制駆動力配分とは連続しない制駆動力配分を求める部分である。
In
なお、ステップ210よりステップ270までの間で、前述の式(19)を用いて制駆動力の再配分を行っているが、これは車両挙動変化をゼロとする制駆動力再配分が行えないときの対策の一例であり、本発明が必ずステップ210〜270の制御を用いなければならないということではない。
Note that the braking / driving force redistribution is performed from
まずステップ210では、必要制駆動力補正量の絶対値|ΔFsi|が最も大きくなっている車輪の必要制駆動力補正量ΔFsiを制駆動力補正量最大値ΔFkとする。
First, in
ステップ220では、この制駆動力補正量最大値の絶対値|ΔFk|が閾値Fthb以下の場合にフラグflg=1に設定し、仮の制駆動力補正量ΔFkrにこの制駆動力補正量最大値ΔFkを設定する。また、制駆動力補正量最大値の絶対値|ΔFk|が閾値Fthbを超えている場合にフラグflg=0に設定し、かつこのとき制駆動力補正量最大値ΔFkが正の値(ゼロである場合を含む)であれば仮の制駆動力補正量ΔFkr=Fthbとし、これに対して制駆動力補正量最大値ΔFkが負の値であれば仮の制駆動力補正量ΔFkr=−Fthbとする。 In step 220, when the absolute value | ΔFk | of the braking / driving force correction amount maximum value is equal to or less than the threshold value Fthb, the flag flg = 1 is set, and the braking / driving force correction amount maximum value is set to the temporary braking / driving force correction amount ΔFkr. ΔFk is set. When the absolute value | ΔFk | of the braking / driving force correction amount maximum value exceeds the threshold value Fthb, the flag flg = 0 is set. At this time, the braking / driving force correction amount maximum value ΔFk is a positive value (zero). In some cases, the provisional braking / driving force correction amount ΔFkr = Fthb. On the other hand, if the braking / driving force correction amount maximum value ΔFk is a negative value, the provisional braking / driving force correction amount ΔFkr = −Fthb. And
上記のフラグflg及び閾値Fthbについて次に説明する。 Next, the flag flg and the threshold value Fthb will be described.
四輪1〜4のうちいずれか一輪の制駆動力が変化した場合あるいは四輪1〜4のうち任意に一輪の制駆動力を変化させた場合に、車両挙動(Fx,Fy,M)を乱さない残り三輪の制駆動力補正量ΔFxiを求めるため上記の式(19)を用いる際には、各輪1〜4の制駆動力変化量が微小であるということが前提条件になっている。従って、この制駆動力補正量最大値ΔFkが十分微小であると仮定できないほど大きな場合には上記の式(19)を用いて残り三輪の制駆動力補正量ΔFxiを正確に求めることが難しくなる。これを判断するフラグが上記のフラグflgであり、制駆動力補正量最大値ΔFkが十分微小であると仮定できないほど大きな変化である場合にフラグflg=0に、またそうでない(つまり制駆動力補正量最大値ΔFkが十分微小であると仮定できるほど小さな変化である)場合にフラグflg=1に設定している。
When the braking / driving force of any one of the four
また、この十分微小であると仮定できる制駆動力変化量最大値の絶対値が閾値Fthbであり、制駆動力補正量最大値ΔFkがこの閾値Fthb以上の場合には、必要制駆動力補正量の絶対値|ΔFsi|が最も大きくなっている車輪の制駆動力が閾値Fthbだけ変化したと仮置きして、後述するステップ250で各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4を補正すると共に、ステップ270で制駆動力補正量最大値ΔFkが正の値(ゼロである場合を含む)であるならΔFk←ΔFk−Fthbとし、制駆動力補正量最大値ΔFkが負の値であるならΔFk←ΔFk+Fthbとすることによって制駆動力補正量最大値ΔFkを小さくする。
Further, when the absolute value of the braking / driving force variation maximum value that can be assumed to be sufficiently small is the threshold value Fthb, and the braking / driving force correction amount maximum value ΔFk is equal to or larger than the threshold value Fthb, the necessary braking / driving force correction amount is obtained. Assuming that the braking / driving force of the wheel having the largest absolute value | ΔFsi | of the wheel has changed by the threshold value Fthb, the braking / driving force Fx 1 , Fx 2 , Fx of each of the
ステップ220〜270でのこれらの操作を制駆動力補正量最大値ΔFkが十分小さくなる(つまりステップ220で|ΔFk|≦Fthbとなる)まで繰り返すことによって、制駆動力補正量最大値ΔFkが十分微小であると仮定できないほど大きい場合でも、残り三輪の制駆動力補正量ΔFxiを得ることができる。本実施形態では上記の閾値Fthbを車両重量W[N]の4%(つまり0.04W)としている。 By repeating these operations in steps 220 to 270 until the braking / driving force correction amount maximum value ΔFk becomes sufficiently small (that is, | ΔFk | ≦ Fthb in step 220), the braking / driving force correction amount maximum value ΔFk is sufficient. even if that can not be assumed to be small large, it is possible to obtain the braking and driving force correction amount DerutaFx i rest Miwa. In the present embodiment, the threshold value Fthb is 4% of the vehicle weight W [N] (that is, 0.04 W).
ステップ230では、ステップ30,80で推定(算出)した各輪1〜4の輪荷重W1,W2,W3,W4と各輪1〜4の横滑り角β1,β2,β3,β4から、各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4を求める。
In step 230, the wheel loads W 1 , W 2 , W 3 , and W 4 estimated (calculated) in
ここで、各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4の求め方を左前輪1の場合を例にとって説明する。
Here, how to obtain the tire lateral force sensitivities k 1 , k 2 , k 3 , and k 4 with respect to the change in braking / driving force of each of the
まず、前提としてコントローラ8には左前輪1の輪荷重W1,横滑り角β1毎に、左前輪1の制駆動力とタイヤ横力との関係を表す複数のタイヤ特性のマップを用意しておく。この複数のタイヤ特性は、ステップ100のところで既に述べた通り、図14のように設定されている。
First, as a premise, the
左前輪1の制駆動力Fx1に対応するタイヤ横力Fy1と、左前輪1の制駆動力Fx1+dFx1に対応するタイヤ横力Fy1+dFy1とを、複数のタイヤ特性のうち輪荷重W1、横滑り角β1のときのタイヤ特性を選択して求め、左前輪1の制駆動力変化に対するタイヤ横力の感度k1を次の式(34a)により求める。
A tire lateral force Fy 1 corresponding to the braking / driving force Fx 1 of the left
k1={(Fy1+dFy1)−Fy1}/{(Fx1+dFx1)−Fx1}
=dFy1/dFx1 …(34a)
式(34a)は上記の式(5)と基本的に同じ式である。式(34a)右辺の分母のdFx1は左前輪1の輪荷重W1と比較して十分微小な制駆動力変化量[N](dFx1>0)であり、予め与えておく。つまり、左前輪1の制駆動力Fx1+dFx1に対応するタイヤ横力Fy1+dFy1のマップ値から、左前輪1の制駆動力Fx1に対応するタイヤ横力Fy1のマップ値を差し引いて、左前輪1の制駆動力Fx1が微小な変化量dFx1だけ変化したときのタイヤ横力の変化量dFy1を求めることによって、左前輪1の制駆動力Fx1の変化に対するタイヤ横力Fy1の感度k1を上記の式(34a)により求めることができる。
k 1 = {(Fy 1 + dFy 1 ) −Fy 1 } / {(Fx 1 + dFx 1 ) −Fx 1 }
= DFy 1 / dFx 1 (34a)
Expression (34a) is basically the same as Expression (5) above. An equation (34a) dFx 1 in the denominator of the right-hand side as compared to the wheel load W 1 of the left
残り3つの車輪2〜4についても同様に、各輪2〜4の輪荷重Wi,横滑り角βiに対応した複数のタイヤ特性のマップを用意しておき、各輪2〜4の輪荷重W2,W3,W4と比較して十分微小な変化量dFx2,dFx3,dFx4を定義して、各輪2〜4の制駆動力Fx2,Fx3,Fx4が微小な変化量dFx2,dFx3,dFx4だけ変化したときのタイヤ横力の感度k2,k3,k4を次の式により求める。
Similarly, for the remaining three
k2=dFy2/dFx2 …(34b)
k3=dFy3/dFx3 …(34c)
k4=dFy4/dFx4 …(34d)
このときさらにタイヤ特性を、各輪1〜4の輪荷重W1,W2,W3,W4、各輪1〜4の横滑り角β1,β2,β3,β4に加えて、各輪1〜4の路面摩擦係数μ1,μ2,μ3,μ4毎に用意することによって各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4をさらに高精度に求めることができる。
k 2 = dFy 2 / dFx 2 (34b)
k 3 = dFy 3 / dFx 3 (34c)
k 4 = dFy 4 / dFx 4 (34d)
At this time, the tire characteristics are further added to the wheel loads W 1 , W 2 , W 3 , W 4 of the
ステップ240では、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めた車両挙動を乱さない各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を求める。
In step 240, the braking / driving force correction amounts ΔFx 1 , ΔFx 2 , ΔFx 3 of the wheels 1 to 4 that do not disturb the vehicle behavior obtained in
具体的には各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4の比を次の式(35)の通りにすれば、車両挙動の変化をゼロに抑えることができるので、例えば、左前輪1の必要制駆動力補正量の絶対値|ΔFs1|が4つの必要制駆動力補正量の絶対値|ΔFs1|,|ΔFs2|,|ΔFs3|,|ΔFs4|の中で最も大きい場合には、各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を、ステップ230で求めた各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4と、ステップ40で求めている各輪1〜4の舵角δ1,δ2,δ3,δ4とから、次の4つの式(36)〜式(39)により求める。
Specifically, if the ratio of the braking / driving force correction amounts ΔFx 1 , ΔFx 2 , ΔFx 3 , and ΔFx 4 of each
ΔFx1:ΔFx2:ΔFx3:ΔFx4
={(Lt/Ll)(h4−h2)+h2(h4−h3)}
/(cosδ1−k1sinδ1)
:{(−Lt/Ll)(h3−h1)−h1(h4−h3)}
/(cosδ2−k2sinδ2)
:{(−Lt/Ll)(h4−h2)−h4(h2−h1)}
/(cosδ3−k3sinδ3)
:{(Lt/Ll)(h3−h1)+h3(h2−h1)}
/(cosδ4−k4sinδ4)
…(35)
ΔFx1=ΔFkr …(36)
ΔFx2=[{−(Lt/Ll)(h3−h1)−h1(h4−h3)}
/{(Lt/Ll)(h4−h2)+h2(h4−h3)}]
・{(cosδ1−k1sinδ1)/(cosδ2−k2sinδ2)}ΔFkr
…(37)
ΔFx3=[{−(Lt/Ll)(h4−h2)−h4(h2−h1)}
/{(Lt/Ll)(h4−h2)+h2(h4−h3)}]
・{(cosδ1−k1sinδ1)/(cosδ3−k3sinδ3)}ΔFkr
…(38)
ΔFx4=[{(Lt/Ll)(h3−h1)+h3(h2−h1)}
/{(Lt/Ll)(h4−h2)+h2(h4−h3)}]
・{(cosδ1−k1sinδ1)/(cosδ4−k4sinδ4)}ΔFkr
…(39)
ここで、式(35)、式(37)〜式(39)において係数h1,h2,h3,h4は次の値であるので、これらの係数h1,h2,h3,h4も、ステップ230で求めた各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4と、ステップ40で求めている各輪1〜4の舵角δ1,δ2,δ3,δ4とから算出する。
ΔFx 1 : ΔFx 2 : ΔFx 3 : ΔFx 4
= {(L t / L l ) (h 4 −h 2 ) + h 2 (h 4 −h 3 )}
/ (Cosδ 1 −k 1 sinδ 1 )
: {(- L t / L l) (h 3 -h 1) -h 1 (h 4 -h 3)}
/ (Cosδ 2 −k 2 sinδ 2 )
: {(- L t / L l) (h 4 -h 2) -h 4 (h 2 -h 1)}
/ (Cosδ 3 −k 3 sinδ 3 )
: {(L t / L l ) (h 3 −h 1 ) + h 3 (h 2 −h 1 )}
/ (Cosδ 4 −k 4 sinδ 4 )
... (35)
ΔFx 1 = ΔFkr (36)
ΔFx 2 = [{− (L t / L l ) (h 3 −h 1 ) −h 1 (h 4 −h 3 )}
/ {(L t / L l ) (h 4 −h 2 ) + h 2 (h 4 −h 3 )}]
{(Cosδ 1 −k 1 sinδ 1 ) / (cosδ 2 −k 2 sinδ 2 )} ΔFkr
... (37)
ΔFx 3 = [{− (L t / L l ) (h 4 −h 2 ) −h 4 (h 2 −h 1 )}
/ {(L t / L l ) (h 4 −h 2 ) + h 2 (h 4 −h 3 )}]
{(Cosδ 1 −k 1 sinδ 1 ) / (cosδ 3 −k 3 sinδ 3 )} ΔFkr
... (38)
ΔFx 4 = [{(L t / L l ) (h 3 −h 1 ) + h 3 (h 2 −h 1 )}
/ {(L t / L l ) (h 4 −h 2 ) + h 2 (h 4 −h 3 )}]
· {(Cosδ 1 −k 1 sinδ 1 ) / (cosδ 4 −k 4 sinδ 4 )} ΔFkr
... (39)
Here, in the equations (35) and (37) to (39), the coefficients h 1 , h 2 , h 3 , h 4 are the following values, so these coefficients h 1 , h 2 , h 3 , h 4 is also the sensitivity of the tire lateral force k 1 , k 2 , k 3 , k 4 to the change in braking / driving force of each
h1=(sinδ1+k1cosδ1)/(cosδ1−k1sinδ1) …(補22a)
h2=(sinδ2+k2cosδ2)/(cosδ2−k2sinδ2) …(補22b)
h3=(sinδ3+k3cosδ3)/(cosδ3−k3sinδ3) …(補22c)
h4=(sinδ4+k4cosδ4)/(cosδ4−k4sinδ4) …(補22d)
式(35)は前述の式(19)と同じである。また、3つの式(37)〜式(39)は前述の3つの式(16)〜式(18)と基本的に同じものである。式(補22a)〜式(補22d)は前述の式(補3a)〜式(補3d)と同じである。
h 1 = (sin δ 1 + k 1 cos δ 1 ) / (cos δ 1 −k 1 sin δ 1 ) (Supplement 22a)
h 2 = (sin δ 2 + k 2 cos δ 2 ) / (cos δ 2 −k 2 sin δ 2 ) (Supplement 22b)
h 3 = (sin δ 3 + k 3 cos δ 3 ) / (cos δ 3 −k 3 sin δ 3 ) (Supplement 22c)
h 4 = (sin δ 4 + k 4 cos δ 4 ) / (cos δ 4 −k 4 sin δ 4 ) (Supplement 22d)
Expression (35) is the same as Expression (19) described above. Further, the three expressions (37) to (39) are basically the same as the above-described three expressions (16) to (18). Expressions (complement 22a) to (complement 22d) are the same as the aforementioned expressions (complement 3a) to (complement 3d).
残り三輪2〜4のうちのいずれか一輪の必要制駆動力補正量の絶対値|ΔFsi|(i=2〜4)が4つの必要制駆動力補正量の絶対値|ΔFs1|,|ΔFs2|,|ΔFs3|,|ΔFs4|の中で最も大きい場合についても、同様にして各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を求める。
The absolute value | ΔFs i | (i = 2 to 4) of the required braking / driving force correction amount for any one of the remaining three
なお、本実施形態では各輪1〜4の舵角δ1,δ2 ,δ3 ,δ4をステップ40で直接計測し、4つの制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4の比を与える式(35)及び各制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を与える式(36)〜式(39)に基づいて各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を求めている。これは即ち、各輪1〜4の制駆動力やタイヤ横力等がサスペンションに加わることによって発生する舵角変化(コンプライアンスステア,ロールステア等)を考慮して各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を求めていることに他ならない。従って、各輪1〜4の舵角を、直接検出する手段を有さないためにステアリング5の回転角θとステアリングギア比等から推定するような車両では、このサスペンション特性から舵角変化量を推定し、各輪1〜4の舵角の推定値を補正することによって各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4をより高精度に求めることができる。
In the present embodiment, the steering angles δ 1 , δ 2 , δ 3 , δ 4 of the
ステップ250では、ステップ240で求めた各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4で各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4を補正する、つまり次の式(40)〜式(43)により各輪の制駆動力Fx1,Fx2,Fx3,Fx4を算出する。
In step 250, the braking / driving force correction amounts ΔFx 1 , ΔFx 2 , ΔFx 3 , ΔFx 4 of the
Fx1←Fx1+ΔFx1 …(40)
Fx2←Fx2+ΔFx2 …(41)
Fx3←Fx3+ΔFx3 …(42)
Fx4←Fx4+ΔFx4 …(43)
ステップ260では、フラグflgをみる。フラグflg=1(制駆動力補正量最大値ΔFkが十分微小であると仮定できるほど小さな変化である)ならばステップ300へ進む。
Fx 1 <-Fx 1 + ΔFx 1 (40)
Fx 2 <-Fx 2 + ΔFx 2 (41)
Fx 3 <-Fx 3 + ΔFx 3 (42)
Fx 4 <-Fx 4 + ΔFx 4 (43)
In
フラグflg=0(制駆動力補正量最大値ΔFkが十分微小であると仮定できないほど大きな変化である)ときには制駆動力補正量最大値ΔFkを小さくするため、ステップ270に進み、そのときの制駆動力補正量最大値ΔFkが正の値(ゼロである場合を含む)であれば、制駆動力補正量最大値ΔFkより閾値Fthbだけ差し引いた値を改めて制駆動力補正量最大値ΔFkとした後、ステップ220の操作を実行する。これに対して、そのときの制駆動力補正量最大値ΔFkが負の値であれば、制駆動力補正量最大値ΔFkに閾値Fthbを加算した値を改めて制駆動力補正量最大値ΔFkとした後、ステップ220の操作を実行する。 When the flag flg = 0 (the change is so large that the braking / driving force correction amount maximum value ΔFk cannot be assumed to be sufficiently small), the process proceeds to step 270 to reduce the braking / driving force correction amount maximum value ΔFk. If the driving force correction amount maximum value ΔFk is a positive value (including a case where the driving force correction amount maximum value ΔFk is zero), the value obtained by subtracting the threshold Fthb from the braking / driving force correction amount maximum value ΔFk is again set as the braking / driving force correction amount maximum value ΔFk. Thereafter, the operation of step 220 is executed. On the other hand, if the braking / driving force correction amount maximum value ΔFk at that time is a negative value, a value obtained by adding the threshold Fthb to the braking / driving force correction amount maximum value ΔFk is changed to the braking / driving force correction amount maximum value ΔFk. After that, the operation of step 220 is executed.
2度目のステップ220において、制駆動力補正量最大値の絶対値|ΔFk|が閾値Fthb以下に収まりフラグflg=0となれば制駆動力補正量最大値ΔFkが十分微小であると仮定できるほど小さな変化になったと判断し、ステップ230、240、250で各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を求めて各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4を補正した後、ステップ260よりステップ300に進む。
In the second step 220, if the absolute value | ΔFk | of the braking / driving force correction amount maximum value falls below the threshold value Fthb and the flag flg = 0, it can be assumed that the braking / driving force correction amount maximum value ΔFk is sufficiently small. It is determined that the change is small, and the braking / driving force correction amounts ΔFx 1 , ΔFx 2 , ΔFx 3 , and ΔFx 4 of each
一方、2度目のステップ220においても制駆動力補正量最大値の絶対値|ΔFk|が閾値Fthb以下に収まらずフラグflg=1となったときには、再びステップ230、240、250で制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を求めて制駆動力Fx1,Fx2,Fx3,Fx4を補正した後、ステップ260よりステップ270へと進み、制駆動力補正量最大値ΔFkを閾値Fthbの分だけさらに小さくしてステップ220に戻る。こうしたステップ270、220の操作の繰り返しによりやがて制駆動力補正量最大値の絶対値|ΔFk|が閾値Fthb以下に収まりフラグflg=0となるので、そのときにはステップ230、240、250で制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を求めて制駆動力Fx1,Fx2,Fx3,Fx4を補正した後、ステップ260よりステップ300に進む。
On the other hand, when the absolute value | ΔFk | of the braking / driving force correction amount maximum value does not fall below the threshold value Fthb and becomes flag flg = 1 in the second step 220, the braking / driving force correction is performed again in steps 230, 240, and 250. After obtaining the amounts ΔFx 1 , ΔFx 2 , ΔFx 3 , ΔFx 4 and correcting the braking / driving forces Fx 1 , Fx 2 , Fx 3 , Fx 4 , the process proceeds from
ステップ300では、このようにして求めた制駆動力Fx1,Fx2,Fx3,Fx4にそれぞれ各輪1〜4の半径Rを乗じた値、即ちトルク指令値(Fx1×R,Fx2×R,Fx3×R,Fx4×R)が得られるようにモータ11〜14を制御する。
In step 300, the longitudinal force Fx 1 obtained in this way, Fx 2, Fx 3, respectively Fx 4 value obtained by multiplying the radius R of the
次に、図13のフローチャートに基づき、ステップ50で求めている要求制駆動力tFと、ステップ60で求めている目標左右輪制駆動力差ΔFと、ステップ30で求めている各輪1〜4の輪荷重W1,W2,W3,W4と、ステップ80で求めている各輪1〜4の横滑り角β1,β2,β3,β4とから、図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を決定するマップを求める方法について説明する。
Next, based on the flowchart of FIG. 13, the required braking / driving force tF obtained in
ステップ510では、10個の変数ht,hΔ,hw1,hw2,hw3,hw4,hβ1,hβ2,hβ3,hβ4にそれぞれ最小値の1を設定する。
In
ステップ520では目標制駆動力の関数tF(ht)を次の式(51)により設定する。
In
tF(ht)=(tFmax−tFmin)÷100×ht+tFmin …(51)
ここで、式(51)右辺のtFmaxは図9(A)のステップ50において目標制駆動力tFが採り得る最大値[N]、式(51)右辺のtFminは図9(A)のステップ50において要求制駆動力tFが採り得る最小値[N]で、図10、図11より予め分かっている。従って、目標制駆動力の関数tF(ht)はhtを変数として最小値tFminから最大値tFmaxまでを動く直線となる。
tF (h t ) = (tF max −tF min ) ÷ 100 × h t + tF min (51)
Here, tF max on the right side of Equation (51) is the maximum value [N] that can be taken by the target braking / driving force tF in
ステップ530では目標左右輪制駆動力差の関数tF(hΔ)を次の式(52)により設定する。
In
ΔF(hΔ)=(ΔFmax−ΔFmin)÷100×hΔ+ΔFmin…(52)
ここで、式(52)右辺のΔFmaxは図9(A)のステップ60において目標左右輪制駆動力差ΔFが採り得る最大値[N]、式(52)右辺のΔFminは図9(A)のステップ60において目標左右輪制駆動力差ΔFが採り得る最小値[N]で、図12より予め分かっている。従って、目標左右輪制駆動力差の関数tF(hΔ)はhΔを変数として最小値ΔFminから最大値ΔFmaxまでを動く直線となる。
ΔF (hΔ) = (ΔF max −ΔF min ) ÷ 100 × hΔ + ΔF min (52)
Here, ΔF max on the right side of equation (52) is the maximum value [N] that can be taken by the target left and right wheel braking / driving force difference ΔF in
ステップ540では左前輪1の輪荷重の関数W1(hw1)を次の式(53)により設定する。
In step 540, the wheel load function W 1 (h w1 ) of the left
W1(hw1)=(W1max−W1min)÷100×hw1+W1min …(53)
ここで、式(53)右辺のW1maxは本車両の左前輪1が採り得る輪荷重の最大値[N]、式(53)右辺のW1minは本車両の左前輪1が採り得る輪荷重の最小値[N]で、予め定めている。従って、左前輪1の関数W1(hw1)はhw1を変数として最小値W1minから最大値W1maxまでを動く直線となる。
W 1 (h w1 ) = (W 1max −W 1 min ) ÷ 100 × h w1 + W 1 min (53)
Here, W 1max on the right side of Expression (53) is the maximum value [N] of the wheel load that can be taken by the left
ステップ550〜570では残り三輪2,3,4の輪荷重の関数W2(hw2),W3(hw3),W4(hw4)を左前輪1と同じく次の式(54)〜式(56)により設定する。
In
W2(hw2)=(W2max−W2min)÷100×hw2+W2min …(54)
W3(hw3)=(W3max−W3min)÷100×hw3+W3min …(55)
W4(hw4)=(W4max−W4min)÷100×hw4+W4min …(56)
ここで、式(54)〜式(56)右辺のW2max,W3max,W4maxは各輪2,3,4が採り得る輪荷重の最大値[N]、式(54)〜式(56)右辺のW2min,W3min,W4minは各輪2,3,4が採り得る輪荷重の最小値[N]で、これらも予め定めている。従って、残り三輪2〜4の関数W2(hw2),W3(hw3),W4(hw4)はhw2,hw3,hw4を変数として最小値W2min,W3min,W4minから最大値W2max,W3max,W4maxまでを動く直線となる。
W 2 (h w2) = ( W 2max -W 2min) ÷ 100 × h w2 + W 2min ... (54)
W 3 (h w3 ) = (W 3max −W 3min ) ÷ 100 × h w3 + W 3min (55)
W 4 (h w4 ) = (W 4max −W 4min ) ÷ 100 × h w4 + W 4min (56)
Here, W 2max , W 3max , W 4max on the right side of the equations (54) to (56) are the maximum wheel loads [N] that can be taken by the
ステップ580では左前輪1の横滑り角の関数β1(hβ1)を次の式(57)により設定する。
In
β1(hβ1)=(β1max−β1min)÷100×hβ1+β1min …(57)
ここで、式(57)右辺のβ1maxは本車両の左前輪1が採り得る横滑り角の最大値[rad]、式(57)右辺のβ1minは本車両の左前輪1が採り得る横滑り角の最小値[rad]で、予め定めている。従って、左前輪1の横滑り角の関数β1(hβ1)はhβ1を変数として最小値β1minから最大値β1maxまでを動く直線となる。
β 1 (hβ 1 ) = (β 1max −β 1 min ) ÷ 100 × hβ 1 + β 1min (57)
Here, β 1max on the right side of equation (57) is the maximum value [rad] of the skid angle that can be taken by the left
ステップ590〜610では残り三輪2,3,4の横滑り角の関数β2(hβ2),β3(hβ3),β4(hβ4)を左前輪1と同じく次の式(58)〜式(60)により設定する。
In Steps 590 to 610, the functions β 2 (hβ 2 ), β 3 (hβ 3 ), and β 4 (hβ 4 ) of the remaining slip angles of the remaining three
β2(hβ2)=(β2max−β2min)÷100×hβ2+β2min …(58)
β3(hβ3)=(β3max−β3min)÷100×hβ3+β3min …(59)
β4(hβ4)=(β4max−β4min)÷100×hβ4+β4min …(60)
ここで、式(58)〜式(60)右辺のβ2max,β3max,β4maxは残り三輪2,3,4が採り得る輪荷重の最大値[N]、式(58)〜式(60)右辺のβ2min,β3min,β4minは、残り三輪2,3,4が採り得る輪荷重の最小値[N]で、これらも予め定めている。従って、残り三輪の関数β2(hβ2),β3(hβ3),β4(hβ4)はhβ2,hβ3,hβ4を変数として最小値β2min,β3min,β4minから最大値β2max,β3max,β4maxまでを動く直線となる。
β 2 (hβ 2) = ( β 2max -β 2min) ÷ 100 ×
β 3 (hβ 3 ) = (β 3max −β 3min ) ÷ 100 × hβ 3 + β 3min (59)
β 4 (hβ 4 ) = (β 4max −β 4min ) ÷ 100 × hβ 4 + β 4min (60)
Here, β 2max , β 3max , β 4max on the right side of the expressions (58) to (60) are the maximum wheel loads [N] that the remaining three
ステップ620では、目標制駆動力の関数tF(ht)と、目標左右輪制駆動力差の関数ΔF(hΔ)とから、各輪1〜4の制駆動力配分{Fx1,Fx2,Fx3,Fx4}を、図9(A)のステップ70と同じく、上記の式(28)及び式(29)と同様の式である次式(補23)、(補24)を用いて求める。
In step 620, the braking / driving force distribution {Fx 1 , Fx 2 , Fx 2 for each
Fx1=Fx3=tF(ht)/4−ΔF(hΔ)/4 …(補23)
Fx2=Fx4=tF(ht)/4+ΔF(hΔ)/4 …(補24)
ステップ630では、後述する図15のフローチャートに従って、ステップ540〜620で設定した各輪1〜4の輪荷重の関数W1(hw1),W2(hw2),W3(hw3),W4(hw4)と、各輪1〜4の横滑り角の関数β1(hβ1),β2(hβ2),β3(hβ3),β4(hβ4)において、図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}(j=1,2,…)を求める。
Fx 1 = Fx 3 = tF (h t ) / 4−ΔF (hΔ) / 4 (Supplement 23)
Fx 2 = Fx 4 = tF (h t ) / 4 + ΔF (hΔ) / 4 (Supplement 24)
In step 630, wheel load functions W 1 (h w1 ), W 2 (h w2 ), W 3 (h w3 ), wheel loads set in steps 540 to 620 according to the flowchart of FIG. In W 4 (h w4 ) and the functions β 1 (hβ 1 ), β 2 (hβ 2 ), β 3 (hβ 3 ), β 4 (hβ 4 ) of the sideslip angles of the
ステップ640〜830は10個の変数ht,hΔ,hw1,hw2,hw3,hw4,hβ1,hβ2,hβ3,hβ4の一つ一つを1ずつ101まで増やす部分であり、ステップ650において変数hβ4が101となったときにステップ660へ、ステップ670において変数hβ3が101となったときにステップ680へ、ステップ690において変数hβ2が101となったときにステップ700へ、ステップ710において変数hβ1が101となったときにステップ720へ、ステップ730において変数hw4が101となったときにステップ740へ、ステップ750において変数hw3が101となったときにステップ760へ、ステップ770において変数hw2が101となったときにステップ780へ、ステップ790において変数hw1が101となったときにステップ800へ、ステップ810において変数hΔが101となったときにステップ820へとそれぞれ進み、ステップ830において変数htが101となったときに図13の処理を終了する。これを逆に言えば、ステップ640〜830を経ることによって、ステップ630を必ず通る。ということは、各変数が1より100までの数を採りうる10個の変数ht,hΔ,hw1,hw2,hw3,hw4,hβ1,hβ2,hβ3,hβ4の全ての組み合わせに対してステップ630の操作が実行されることを意味する。
Steps 640 to 830 are steps for increasing each of the ten variables h t , hΔ, h w1 , h w2 , h w3 , h w4 , hβ 1 , hβ 2 , hβ 3 , hβ 4 to 101 one by one. Yes, when the variable hβ 4 becomes 101 at
このように、図13のフローによって、本車両が採り得る、値の異なる100個の目標制駆動力tFと、値の異なる100個の目標左右輪制駆動力差ΔFと、値の異なる100個の輪荷重W1,W2,W3,W4と、値の異なる100個の横滑り角β1,β2,β3,β4の全ての組み合わせに対して、図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を全て求めてマップ化することができる。なお、100という数値は一つの例であり、これに限定されるものでない。 As described above, according to the flow of FIG. 13, 100 target braking / driving forces tF having different values, 100 target left / right wheel braking / driving force differences ΔF having different values, and 100 having different values can be adopted by the vehicle. 9A for all combinations of the wheel loads W 1 , W 2 , W 3 , W 4 and 100 sideslip angles β 1 , β 2 , β 3 , β 4 having different values. It is possible to obtain and map all the sets {Fx 1 (j), Fx 2 (j), Fx 3 (j), Fx 4 (j)} of braking / driving force distributions that realize the same vehicle behavior as obtained in 100. it can. In addition, the numerical value of 100 is an example, and is not limited to this.
次に、図15のフローチャートに基づき、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動(Fx,Fy,M)を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}(j=1,2,…)を求める方法について説明する。
Next, based on the flowchart of FIG. 15, the braking / driving force distribution for realizing the same vehicle behavior (Fx, Fy, M) as obtained in
ステップ1010では、まず現在の制駆動力配分{Fx1,Fx2,Fx3,Fx4}(図13のステップ620で求めている)を最初の探索ステップ数lの制駆動力配分{Fx1(l),Fx2(l),Fx3(l),Fx4(l)}に設定すると共に、探索ステップ数lに初期値の1を入れる。これで、探索ステップ数lが初期値の1のときの制駆動力配分、つまり最初の制駆動力配分{Fx1(1),Fx2(1),Fx3(1),Fx4(1)}が定まる。この最初の制駆動力配分{Fx1(1),Fx2(1),Fx3(1),Fx4(1)}は後述するステップ1048において必要になるので、コントーラ8のメモリに記憶しておく。ステップ1010ではまた、制駆動力補正基準量Δ[N]に0.2Nを設定する。
In step 1010, the current braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } (obtained in step 620 in FIG. 13) is first determined as the braking / driving force distribution {Fx 1 (l), Fx 2 (l), Fx 3 (l), Fx 4 (l)} and an initial value of 1 is entered in the number of search steps l. Thus, the braking / driving force distribution when the number of search steps l is the
ここで、探索ステップ数lは、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分をいくつ求めたかを表す数である。制駆動力補正基準量Δは、後述するステップ1022,1034において上記の式(19)あるいは式(20)を用いて車両挙動を変化させない制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4を求める際に使用する。
Here, the number of search steps l is a number representing how many braking / driving force distributions for realizing the same vehicle behavior as those obtained in
ステップ1012では、各輪1〜4の制駆動力変化量が十分微小であると仮定したときにおける車両挙動を変化させない制駆動力補正量の比ΔFx1:ΔFx2:ΔFx3:ΔFx4を求め、左前輪1の制駆動力補正量ΔFx1の絶対値と閾値thを比較する。これは、左前輪1の制駆動力補正量ΔFx1(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができるか否かをみる部分である。
In step 1012, the ratio of the braking / driving force correction amount ΔFx 1 : ΔFx 2 : ΔFx 3 : ΔFx 4 that does not change the vehicle behavior when it is assumed that the braking / driving force change amount of each of the
左前輪1の制駆動力補正量ΔFx1の絶対値が閾値th以上の場合には左前輪1の制駆動力補正量ΔFx1(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができると判断し、ステップ1014に進んでフラグFl=1とし、これに対して左前輪1の制駆動力補正量ΔFx1の絶対値が閾値th未満の場合には左前輪1の制駆動力補正量ΔFx1(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができないと判断し、ステップ1012よりステップ1016に進んでフラグFl=4とする。
When the absolute value of the braking / driving force correction amount ΔFx 1 for the left
ここで、上記の制駆動力補正量の比ΔFx1:ΔFx2:ΔFx3:ΔFx4は次のようにして求める。 Here, the ratio ΔFx 1 : ΔFx 2 : ΔFx 3 : ΔFx 4 of the braking / driving force correction amount is obtained as follows.
まず、ステップ1010で設定している制駆動力配分{Fx1(1),Fx2(1),Fx3(1),Fx4(1)}における各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4を求める。求め方は図9(B)のステップ230で既に述べているので、ここでの説明は省略する。
First, the braking / driving force distributions {Fx 1 (1), Fx 2 (1), Fx 3 (1), Fx 4 (1)} set in step 1010 to the braking / driving force changes of the
そして、求めた各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4から上記の式(20)を用いて、あるいは求めた各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4と、各輪1〜4の舵角δ1,δ2,δ3,δ4(図9(A)のステップ40で算出されている)とから上記の式(19)を用いて制駆動力補正量の比ΔFx1:ΔFx2:ΔFx3:ΔFx4を求めることができる。
The tire lateral force sensitivities k 1 , k 2 , k 3 , and k 4 with respect to the obtained braking / driving force changes of the
上記の閾値thを導入した理由は次の通りである。上記の式(20)あるいは上記の式(19)を用いて得られるのは制駆動力補正量の比ΔFx1:ΔFx2:ΔFx3:ΔFx4であるので、いずれか一輪の制駆動力補正量を基準にして、残り三輪の制駆動力補正量を求めなくてはならない。従って、この制駆動力補正量の比ΔFx1:ΔFx2:ΔFx3:ΔFx4がゼロに近い値となる車輪の制駆動力補正量を基準とすると、ゼロ割(ゼロで割ること)等が発生し、コントローラ8での演算が不正確になる。そこで、閾値thを導入し、制駆動力補正量の比ΔFx1:ΔFx2:ΔFx3:ΔFx4がゼロに近い値となる車輪の制駆動力補正量を基準としないように制御を行うのである。具体的には本実施形態では閾値thとして0.01を設定している。
The reason why the threshold th is introduced is as follows. Since the ratio ΔFx 1 : ΔFx 2 : ΔFx 3 : ΔFx 4 of the braking / driving force correction amount is obtained using the above equation (20) or the above equation (19), the braking / driving force correction of any one wheel is performed. Based on the amount, the braking / driving force correction amount for the remaining three wheels must be determined. Accordingly, when the braking / driving force correction amount of the wheel at which the ratio ΔFx 1 : ΔFx 2 : ΔFx 3 : ΔFx 4 of the braking / driving force correction amount is close to zero is used as a reference, zero division (divide by zero) or the like Occurs, and the calculation in the
ステップ1018では、ステップ1010で求めている制駆動力配分{Fx1(l),Fx2(l),Fx3(l),Fx4(l)}における各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4を求める。求め方はステップ1012でも説明したように、図9(B)のステップ230と同様である。
In
ステップ1020ではフラグFlをみる。フラグFl=4の場合(つまり左前輪1の制駆動力補正量ΔFx1(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができないと判断している場合)には、ステップ1022に進む。
In step 1020, the flag Fl is checked. In the case of the flag Fl = 4 (that is, the braking / driving force distribution {Fx 1 (l + 1) which is a value when the number of
ステップ1022では、右後輪4の駆動力Fx4(1)に制駆動力補正基準量Δを加算したときに、図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求める。
In
この探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}の求め方を次に説明する。
Braking / driving force distribution {Fx 1 (l + 1), Fx 2 (l + 1), Fx 3 (l + 1), Fx 4 (l + 1), which is a value when the number of
まず、ステップ1018で求めているタイヤ横力の感度k1,k2,k3,k4から上記の式(20)を用いて、あるいはステップ1018で求めているタイヤ横力の感度k1,k2,k3,k4と、各輪1〜4の舵角δ1,δ2,δ3,δ4(図9(A)のステップ40で算出されている)とから上記の式(19)を用いて制駆動力補正量の比ΔFx1:ΔFx2:ΔFx3:ΔFx4を求める。
First, the sensitivity k 1 of the tire lateral force is calculated in
その上で、ステップ1010で設定している制駆動力配分{Fx1(l),Fx2(l),Fx3(l),Fx4(l)}、この制駆動力補正量の比ΔFx1:ΔFx2:ΔFx3:ΔFx4及び制駆動力補正基準値Δを用いて、探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を次の式(101)〜式(104)により求め、求めた制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}をコントローラ8のメモリに記憶する。
Then, the braking / driving force distribution {Fx 1 (l), Fx 2 (l), Fx 3 (l), Fx 4 (l)} set in step 1010, and the ratio ΔFx of the braking / driving force correction amount are set. 1 : ΔFx 2 : ΔFx 3 : Braking / driving force distribution {Fx 1 (l + 1), Fx 2 that is a value when the number of
Fx1(l+1)=Fx1(l)+(ΔFx1/ΔFx4)×Δ …(101)
Fx2(l+1)=Fx2(l)+(ΔFx2/ΔFx4)×Δ …(102)
Fx3(l+1)=Fx3(l)+(ΔFx3/ΔFx4)×Δ …(103)
Fx4(l+1)=Fx4(l)+Δ …(104)
ステップ1024では、ステップ1022で求めた探索ステップ数lを1増やしたときの値である右後輪4の制駆動力補正量ΔFx4(l+1)の絶対値と閾値thを比較する。これは、ステップ1022では右後輪4の制駆動力補正量ΔFx4(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めたが、次回にも右後輪4の制駆動力補正量ΔFx4(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができるか否かをみる部分である。
Fx 1 (l + 1) = Fx 1 (l) + (ΔFx 1 / ΔFx 4 ) × Δ (101)
Fx 2 (l + 1) = Fx 2 (l) + (ΔFx 2 / ΔFx 4 ) × Δ (102)
Fx 3 (l + 1) = Fx 3 (l) + (ΔFx 3 / ΔFx 4 ) × Δ (103)
Fx 4 (l + 1) = Fx 4 (l) + Δ (104)
In
探索ステップ数lを1増やしたときの右後輪4の制駆動力補正量ΔFx4(l+1)の絶対値が閾値th以上の場合には、次回にも右後輪4の制駆動力補正量ΔFx4(l)を基準として探索ステップ数を1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができると判断し、ステップ1026、1028、1030、1032を飛ばしてステップ1046に進み、探索ステップ数lを1増やす。
When the absolute value of the braking / driving force correction amount ΔFx 4 (l + 1) of the right
これに対して、ステップ1024において、ステップ1022で求めた探索ステップ数lを1増やしたときの値である右後輪4の制駆動力補正量ΔFx4(l+1)の絶対値が閾値th未満となった場合には、次回に右後輪4の制駆動力補正量ΔFx4(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができないと判断し、今度は左前輪1の制駆動力補正量ΔFx1(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めさせるため、ステップ1026、1028、1030、1032に進む。
On the other hand, in
ステップ1026では、探索ステップ数lを1増やしたときの値である左前輪1の制駆動力Fx1(l+1)と探索ステップ数lを1増やす前の値である左前輪1の制駆動力Fx1(l)の差であるFx1(l+1)−Fx1(l)とゼロを比較する。Fx1(l+1)とFx1(l)の差が正の値(ゼロの場合を含む)であればステップ1030に進み、制駆動力補正基準値Δ=|Δ|とする(制駆動力補正基準値Δの符号を正とする)。これに対して、Fx1(l+1)とFx1(l)の差が負の値であればステップ1028に進み、制駆動力補正基準値Δ=−|Δ|とする(制駆動力補正基準値Δの符号を負とする)。
In step 1026, the braking / driving force Fx 1 (l + 1) of the left
ステップ1032では次回に左前輪1の制駆動力補正量ΔFx1(l)を基準として、探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めるためフラグFl=1とした後、ステップ1046に進んで探索ステップ数lを1増やす。
In step 1032, the braking / driving force distribution {Fx 1 (l + 1), Fx which is a value when the number of
一方、ステップ1020でフラグFl=1の場合(つまり左前輪1の制駆動力補正量ΔFx1(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができると判断している場合)には、ステップ1034に進む。
On the other hand, if the flag Fl = 1 in step 1020 (that is, the braking / driving force distribution {Fx 1 ) that is a value when the number of
ステップ1034〜1044の操作はステップ1022〜1032の操作と対をなす操作であるが、操作の内容そのものは同じである。
The operations in steps 1034 to 1044 are operations that make a pair with the operations in
ステップ1034では、左前輪1の制駆動力Fx1(l)に制駆動力補正基準量Δを加算したときに図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求める。
In step 1034, the braking / driving force distribution that realizes the same vehicle behavior as that obtained in
この探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}の求め方を次に説明する。
Braking / driving force distribution {Fx 1 (l + 1), Fx 2 (l + 1), Fx 3 (l + 1), Fx 4 (l + 1), which is a value when the number of
まず、ステップ1018で求めているタイヤ横力の感度k1,k2,k3,k4から上記の式(20)を用いて、あるいはステップ1018で求めているタイヤ横力の感度k1,k2,k3,k4と、各輪1〜4の舵角δ1,δ2,δ3,δ4(図9(A)のステップ40で算出されている)とから上記の式(19)を用いて制駆動力補正量の比ΔFx1:ΔFx2:ΔFx3:ΔFx4を求める。
First, the sensitivity k 1 of the tire lateral force is calculated in
その上で、ステップ1010で設定している制駆動力配分{Fx1(l),Fx2(l),Fx3(l),Fx4(l)}、この制駆動力補正量の比ΔFx1:ΔFx2:ΔFx3:ΔFx4及び制駆動力補正基準値Δを用いて、探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を次の式(105)〜式(108)により求め、求めた制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}をコントローラ8のメモリに記憶する。
Then, the braking / driving force distribution {Fx 1 (l), Fx 2 (l), Fx 3 (l), Fx 4 (l)} set in step 1010, and the ratio ΔFx of the braking / driving force correction amount are set. 1 : ΔFx 2 : ΔFx 3 : Braking / driving force distribution {Fx 1 (l + 1), Fx 2 that is a value when the number of
Fx1(l+1)=Fx1(l)+Δ …(105)
Fx2(l+1)=Fx2(l)+(ΔFx2/Fx1)×Δ …(106)
Fx3(l+1)=Fx3(l)+(ΔFx3/Fx1)×Δ …(107)
Fx4(l+1)=Fx4(l)+(ΔFx4/Fx1)×Δ …(108)
ステップ1036では、ステップ1034で求めた探索ステップ数lを1増やしたときの値である左前輪1の制駆動力補正量ΔFx1(l+1)の絶対値と閾値thを比較する。これは、ステップ1034では左前輪1の制駆動力補正量ΔFx1(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めたが、次回にも左前輪1の制駆動力補正量ΔFx1(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができるか否かをみる部分である。
Fx 1 (l + 1) = Fx 1 (l) + Δ (105)
Fx 2 (l + 1) = Fx 2 (l) + (ΔFx 2 / Fx 1 ) × Δ (106)
Fx 3 (l + 1) = Fx 3 (l) + (ΔFx 3 / Fx 1 ) × Δ (107)
Fx 4 (l + 1) = Fx 4 (l) + (ΔFx 4 / Fx 1 ) × Δ (108)
In step 1036, the absolute value of the braking / driving force correction amount ΔFx 1 (l + 1) of the left
探索ステップ数lを1増やしたときの左前輪1の制駆動力補正量ΔFx1(l+1)の絶対値が閾値th未満の場合には、次回にも左前輪1の制駆動力補正量ΔFx1(l)を基準として探索ステップ数を1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができると判断し、ステップ1038、1040、1042、1044を飛ばしてステップ1046に進み、探索ステップ数lを1増やす。
If the absolute value of the braking / driving force correction amount ΔFx 1 (l + 1) of the left
これに対して、ステップ1036において、ステップ1034で求めた探索ステップ数lを1増やしたときの値である左前輪1の制駆動力補正量ΔFx1(l+1)の絶対値が閾値th未満となった場合には、次回に左前輪1の制駆動力補正量ΔFx1(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めることができないと判断し、今度は右後輪4の制駆動力補正量ΔFx4(l)を基準として探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めさせるため、ステップ1038、1040、1042、1044に進む。
On the other hand, in step 1036, the absolute value of the braking / driving force correction amount ΔFx 1 (l + 1) of the left
ステップ1038では、探索ステップ数lを1増やしたときの値である右後輪4の制駆動力Fx4(l+1)と探索ステップ数lを1増やす前の値である右後輪4の制駆動力Fx4(l)の差であるFx4(l+1)−Fx4(l)とゼロを比較する。Fx4(l+1)とFx4(l)の差が正の値(ゼロの場合を含む)であればステップ1042に進み、制駆動力補正基準値Δ=|Δ|とする(制駆動力補正基準値Δの符号を正とする)。これに対して、Fx4(l+1)とFx4(l)の差が負の値であればステップ1040に進み、制駆動力補正基準値Δ=−|Δ|とする(制駆動力補正基準値Δの符号を負とする)。
In step 1038, the braking / driving force Fx 4 (l + 1) of the right
ステップ1044では次回に右後輪4の制駆動力補正量ΔFx4(l)を基準として、探索ステップ数lを1増やしたときの値である制駆動力配分{Fx1(l+1),Fx2(l+1),Fx3(l+1),Fx4(l+1)}を求めるためフラグFl=4とし、ステップ1046に進んで探索ステップ数lを1増やす。
In step 1044, the braking / driving force distribution {Fx 1 (l + 1), which is a value when the number of
ステップ1048では、探索ステップ数lと100を比較する。探索ステップ数lが100未満であるときにはステップ1018に戻ってステップ1018〜1046の操作を繰り返す。すなわち、ステップ1022とステップ1034のいずれか一方の操作が100回繰り返されるので、探索ステップ数lが100に達したときには、100組の制駆動力配分{Fx1(l),Fx2(l),Fx3(l),Fx4(l)}(l=1〜100)が得られる。
In step 1048, the number of search steps l is compared with 100. When the number of search steps l is less than 100, the process returns to step 1018 and the operations of
また、探索ステップ数lが100に達したときには、最初の制駆動力配分{Fx1(1),Fx2(1),Fx3(1),Fx4(1)}と、探索ステップ数lが100であるときの制駆動力配分{Fx1(l),Fx2(l),Fx3(l),Fx4(l)}とを用いて次の式(109)の評価関数Pを計算し、この評価関数Pと1を比較する。 When the search step number l reaches 100, the first braking / driving force distribution {Fx 1 (1), Fx 2 (1), Fx 3 (1), Fx 4 (1)} and the search step number l The evaluation function P of the following equation (109) is calculated using the braking / driving force distribution {Fx 1 (l), Fx 2 (l), Fx 3 (l), Fx 4 (l)} when the value is 100 Calculate and compare this evaluation function P with 1.
P=(Fx1(l)−Fx1(1))2+(Fx2(l)−Fx2(1))2
+(Fx3(l)−Fx3(1))2+(Fx4(l)−Fx4(1))2
…(109)
式(109)の評価関数Pが1以下になっていれば、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合が一通り求められたと判断し図15のフローを終了する。このとき、制駆動力配分{Fx1(l),Fx2(l),Fx3(l),Fx4(l)}の集合はちょうど100組である。
P = (Fx 1 (l) −Fx 1 (1)) 2 + (Fx 2 (l) −Fx 2 (1)) 2
+ (Fx 3 (l) −Fx 3 (1)) 2 + (Fx 4 (l) −Fx 4 (1)) 2
... (109)
If the evaluation function P of the equation (109) is 1 or less, the braking / driving force distribution for realizing the same vehicle behavior as that obtained in
一方、探索ステップ数lが100に達しても上記の式(109)の評価関数Pが1を超えていればなおもステップ1018〜1048の操作を繰り返す。
On the other hand, even if the number of search steps l reaches 100, if the evaluation function P of the above equation (109) exceeds 1, the operations of
やがて、ステップ1048で探索ステップ数lが100を超えておりかつ上記の式(109)の評価関数Pが1以下となれば、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合が一通り求められたと判断し図15のフローを終了する。この場合には、制駆動力配分{Fx1(l),Fx2(l),Fx3(l),Fx4(l)}の集合は100組を超えている。
Eventually, if the number of search steps l exceeds 100 in step 1048 and the evaluation function P of the above equation (109) becomes 1 or less, the current tire characteristics of each
ここで、ステップ1048において上記の式(109)の評価関数Pが1以下となったとき、つまり最初の制駆動力配分{Fx1(1),Fx2(1),Fx3(1),Fx4(1)}と、探索ステップ数lが100以上のときの制駆動力配分{Fx1(l),Fx2(l),Fx3(l),Fx4(l)}とがほぼ同じになったときに、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合が一通り求められたとしたのは次の理由による。
Here, when the evaluation function P of the above equation (109) becomes 1 or less in step 1048, that is, the first braking / driving force distribution {Fx 1 (1), Fx 2 (1), Fx 3 (1), Fx 4 (1)} and braking / driving force distribution {Fx 1 (l), Fx 2 (l), Fx 3 (l), Fx 4 (l)} when the number of search steps l is 100 or more When it became the same, it was said that a set of braking / driving force distributions for realizing the same vehicle behavior as that obtained in
図14に示すようなタイヤ特性を有する車輪を4つ備える四輪独立駆動車において、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分は、ある一輪(例えば左前輪1)の制駆動力を横軸に採ると、図5のように残り三輪(右前輪2、左後輪3、右後輪4)の制駆動力の軌跡は閉じた輪となる、つまり、図15のフローチャートのように各輪1〜4の制駆動力、例えば左前輪1または右後輪4の制駆動力を制駆動力補正基準値Δずつ変化させて、図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分を求めていくと、制駆動力補正基準値Δが十分微小であれば最初の制駆動力配分に戻ることがシミュレーション上で確認されている。
In a four-wheel independent drive vehicle having four wheels having tire characteristics as shown in FIG. 14, the same vehicle behavior as that obtained in
そこで本実施形態では、図15のフローチャートによって現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分を求めていき、最初の制駆動力配分{Fx1(1),Fx2(1),Fx3(1),Fx4(1)}と、探索ステップ数lが100以上のときの制駆動力配分{Fx1(l),Fx2(l),Fx3(l),Fx4(l)}とがほぼ同じになったときを、現在の各輪1〜4のタイヤ特性において図9(A)のステップ100で求めたと同じ車両挙動を実現する制駆動力配分の集合が一通り求められたものと判断しているのである。
Accordingly, in the present embodiment, the braking / driving force distribution for realizing the same vehicle behavior as that obtained in
ここで、本実施形態の作用を説明する。 Here, the operation of the present embodiment will be described.
本実施形態(請求項1に記載の発明)によれば、図9(A)のステップ100において要求される車両挙動(車両の平面運動での目標となる車両挙動)を決定し、この要求される車両挙動を実現する四輪の制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を図9(A)のステップ101で演算し、図9(B)のステップ120、150において、この演算された制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}の中から、一つの制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択するように構成したので、スリップやモータ過熱等によって四輪のうち二輪以上の制駆動力が変化した場合あるいは四輪のうち二輪以上の制駆動力を変化させる場合においても、従来技術より車両挙動の変化を抑えることができる。
According to the present embodiment (the invention described in claim 1), the vehicle behavior required in
各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4を求めるために必要な、制駆動力Fx1,Fx2,Fx3,Fx4とタイヤ横力Fy1,Fy2,Fy3,Fy4の関係を表すタイヤ特性は、図14に示したように主に輪荷重と横すべり角によって変化することに対応して、本実施形態(請求項2に記載の発明)によれば、目標となる車両挙動を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を、図9(A)のステップ101において各輪1〜4の輪荷重W1,W2,W3,W4と、横すべり角β1,β2,β3,β4とに基づいて演算するので(上記の<1>参照)、各輪1〜4のタイヤ特性の変化がより精密に得られる。これにより、目標となる車両挙動を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}をより精密に求めることができ、スリップやモータ過熱等によって四輪のうち二輪以上の制駆動力が変化した場合あるいは四輪のうち二輪以上の制駆動力を変化させる場合における車両挙動変化を抑えることができる。
Braking / driving forces Fx 1 , Fx 2 , Fx 3 , Fx 4 and tires necessary for obtaining the sensitivity k 1 , k 2 , k 3 , k 4 of the tire lateral force with respect to changes in the braking / driving force of the
本実施形態(請求項3に記載の発明)によれば、目標となる車両挙動を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を、図9(A)のステップ101において各輪1〜4の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4に基づいて演算するので(上記の<2>参照)、各輪1〜4のタイヤ特性の変化がより直接的に得られる。これにより、目標となる車両挙動を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}をより直接的に求めることができ、各輪1〜4のスリップや各モータ11〜14の過熱等によって特に二輪以上の制駆動力が変化した場合あるいは特に二輪以上の制駆動力を変化させる場合における車両挙動変化を抑えることができる。
According to the present embodiment (the invention according to claim 3), a set of braking / driving force distributions {Fx 1 (j), Fx 2 (j), Fx 3 (j), Fx that realizes a target vehicle behavior 4 (j)} is calculated based on the tire lateral force sensitivity k 1 , k 2 , k 3 , k 4 with respect to the change in braking / driving force of each wheel 1-4 in
本実施形態(請求項4に記載の発明)によれば、目標となる車両挙動を実現する制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を、車両の平面運動での車両挙動(Fx,Fy,M)の変化をゼロとする各輪1〜4の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4の比を求める上記の式(35)、式(補22a)〜式(補22d)に基づいて演算するので、不要な車両挙動の変化を抑えて車体横滑り角βをゼロにすることができる。
According to the present embodiment (the invention described in claim 4), the set of braking / driving force distributions {Fx 1 (j), Fx 2 (j), Fx 3 (j), Fx that realize the target vehicle behavior 4 (j)} is a braking / driving force correction amount ΔFx 1 , ΔFx 2 , ΔFx 3 , ΔFx 4 for each of the
本実施形態(請求項5、6に記載の発明)によれば、図9(A)のステップ110で各輪1〜4が路面に伝達可能な制駆動力上限Fsmaxiを推定し(上記<3>参照)、制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を、図9(B)のステップ130においてこの各輪1〜4が路面に伝達可能な制制駆動力上限Fsmaxiに制限するか、または図9(A)のステップ110で各輪1〜4が路面に伝達可能な制駆動力下限Fsminiを推定し(上記<3>参照)、制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を、図9(B)のステップ130においてこの各輪1〜4が路面に伝達可能な制制駆動力下限Fsminiに制限するので、各輪1〜4がスリップやロック等によって路面に伝達できる制駆動力Fx1,Fx2,Fx3,Fx4が変化した場合においても、実現可能な制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択することができる。
According to the present embodiment (the invention described in
本実施形態(請求項7、8に記載の発明)によれば、図9(A)のステップ110で各モータ11〜14(各輪1〜4を駆動または制動するアクチュエータ)が出力可能な制駆動力上限Fdmaxiを推定し(上記<4>参照)、制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を、図9(B)のステップ130においてこの各モータ11〜14が出力可能な制駆動力上限Fdmaxiに制限するか、または図9(A)のステップ110で各モータ11〜14が出力可能な制駆動力下限Fdminiを推定し(上記<4>参照)、制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を、図9(B)のステップ130においてこの各モータ11〜14が出力可能な制駆動力下限Fdminiに制限するので、各モータ11〜14の異常等によって各輪1〜4が出力できる制駆動力が変化した場合においても、実現可能な制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択することができる。
According to the present embodiment (the inventions described in
本実施形態(請求項9に記載の発明)によれば、図9(B)のステップ150において、配列Stに記録された各整数jに対応する制駆動力配分毎に、各輪1〜4のモータロス(各輪1〜4を駆動するときに発生する各エネルギー損失)Ploss1(j),Ploss2(j),Ploss3(j),Ploss4(j)を求め(推定)し、各輪1〜4のエネルギーロスの和Ploss(j)が最小となる制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択するので(上記<7>参照)、各モータ11〜14で各輪1〜4を独立に駆動する車両では、バッテリ9の消費電力を最小化することができ、車両の航続距離の増大等の効果を得ることができる。
According to the present embodiment (the invention described in claim 9), in
本実施形態(請求項10に記載の発明)によれば、図9(B)のステップ150において、配列Stに記録された各整数jに対応する制駆動力配分の中で、タイヤ横力と制駆動力との合力であるタイヤ力を輪荷重で除したタイヤ摩擦円使用率H1(j),H2(j),H3(j),H4(j)を各輪1〜4毎に求め、この求めた各輪1〜4のタイヤ摩擦円使用率H1(j),H2(j),H3(j),H4(j)のうち最大値と最小値の差が最も小さくなる制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択するので(上記<8>参照)、各輪1〜4の路面摩擦係数が急激に低下しても、各輪1〜4のタイヤ摩擦円が飽和しにくくなり車両の操縦安定性をより向上させることができる。
According to the present embodiment (the invention described in claim 10), in
本実施形態(請求項11に記載の発明)によれば、制駆動力配分が変化する過渡状態で各輪1〜4の制駆動力変化が連続となる制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択するので、過渡状態において各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4が不連続となることによって発生する不要な衝撃をなくすことができ、車両の操縦性を高めることができる。
According to the present embodiment (the invention described in claim 11), the braking / driving force distribution {Fx 1 , Fx 2 , in which the braking / driving force change of each of the
本実施形態(請求項12、13に記載の発明)によれば、図9(B)のステップ210〜270において、各輪1〜4において、スリップあるいは車輪ロックを起こさない各輪1〜4の制駆動力上限Fsmaxiまたはスリップあるいは車輪ロックを起こさない各輪1〜4の制駆動力下限Fsminiにより制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}が制限され、制駆動力配分{Fx1,Fx2,Fx3,Fx4}を連続的に変化させることができない場合に、各輪1〜4の制駆動力変化が不連続となることを許可し、ステップ70で定まっている制駆動力配分とは連続しない制駆動力配分を求める(つまり各輪1〜4の制駆動力変化が不連続とできるようにする)か、または各モータ11〜14が過熱しまたは破損しないようにする各輪1〜4の制駆動力上限Fdmaxiまたは各モータ11〜14が過熱しまたは破損しないようにする各輪1〜4の制駆動力下限Fdminiにより制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}が制限され、制駆動力配分{Fx1,Fx2,Fx3,Fx4}を連続的に変化させることができない場合に、各輪1〜4の制駆動力変化が不連続となることを許可し、ステップ70で定まっている制駆動力配分とは連続しない制駆動力配分を求める(つまり各輪1〜4の制駆動力変化が不連続とできるようにする)ので、スリップやモータ過熱などによって各輪1〜4の制駆動力配分の集合が制限され、四輪1〜4の制駆動力配分を連続的に変化させることができない場合においても、車両挙動の変化を抑えることができる。
According to the present embodiment (the inventions described in
図18は第2実施形態の四輪独立駆動車の概略構成図で、第1実施形態の図1と置き換わるものである。図1と同一部分には同一の番号をつけている。 FIG. 18 is a schematic configuration diagram of the four-wheel independent drive vehicle of the second embodiment, which replaces FIG. 1 of the first embodiment. The same parts as those in FIG.
図18において図1と相違するのは、ジェネレータ51、エンジン52、コンバータ53からなる発電装置を追加した点にある。ジェネレータ51とエンジン52とは直接的または間接的に連結され、ジェネレータ51はエンジン52により駆動されて発電し、発電した電力は直接的に各モータ11〜14に供給されるようになっている。
18 differs from FIG. 1 in that a power generation device including a
図19(A)及び図19(B)は第2実施形態の各モータ11〜14へのトルク配分を実行するためのもので、第1実施形態の図9(A)及び図9(B)と置き換わるものである。図19(A)及び図19(B)も操作の手順を示し、一定時間毎に実行するものではない。図19(A)、図19(B)において第1実施形態の図9(A)、図9(B)と同一部分には同一のステップ番号を付している。
19 (A) and 19 (B) are for executing torque distribution to the
第2実施形態は、ステップ101で得られている実現可能な制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}の中から、各制駆動力配分を実現するための総消費電力が車両の出力可能な電力を上回ることない一つの制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択するようにしたものである。
In the second embodiment, from the set of feasible braking / driving force distributions {Fx 1 (j), Fx 2 (j), Fx 3 (j), Fx 4 (j)} obtained in
第1実施形態の図19(A)及び図19(B)との違いを主に説明すると、図19(A)においては第1実施形態の図9(A)に対してステップ11、91、102、103を新たに加えている。また、図19(B)においては第1実施形態の図9(B)のステップ120に代えてステップ111、112、113、114、115、116、117で、また第1実施形態の図9(B)のステップ210、220、230、240、250、260、270に代えてステップ211、221、231、241、251で置き換えている。
19A and 19B of the first embodiment will be mainly described. In FIG. 19A, steps 11 and 91 are compared with FIG. 9A of the first embodiment. 102 and 103 are newly added. In FIG. 19B, steps 111, 112, 113, 114, 115, 116, 117 are used instead of step 120 of FIG. 9B of the first embodiment, and FIG. B)
具体的に説明すると、図19(A)においてステップ11ではジェネレータ51の発電電力Pg[W]を検出し、電流センサ(図示しない)により検出されるバッテリ9の電流値を積算することによりバッテリ9の蓄電量Bc[%]を算出する。
More specifically, in
ステップ91ではエンジン52の目標運転点(目標回転速度と目標トルク)を設定し、その目標運転点でエンジン52を運転する。エンジン52の目標運転点の設定方法を次に述べる。
In step 91, a target operating point (target rotational speed and target torque) of the
まず、エンジン52の目標出力パワーtPe[W]を、ステップ50で求めている目標制駆動力tF(車両前後方向力の目標値)と、ステップ10で求めている車体速度Vとを用いて、次の(110)により算出する。
First, the target output power tPe [W] of the
tPe=tF×V−100×eSOC …(110)
式(110)右辺のeSOC[%]はステップ11で検出しているバッテリ9の蓄電量Bcから50%を引いた値である。すなわち、蓄電量eSOC(=Bc−50)が50%より小さくなると充電、大きくなると放電するように目標出力パワーtPeを定めている。
tPe = tF × V−100 × e SOC (110)
E SOC [%] on the right side of Expression (110) is a value obtained by subtracting 50% from the charged amount Bc of the
そして、この目標出力パワーtPeから図20を内容とするマップを参照してエンジン52の目標回転速度と目標トルクを定める。図20においてEf1〜Ef6[cm3/min]はエンジンの燃料消費率(Ef6>Ef5>Ef4>Ef3>Ef2>Ef1)、Pe eq1〜Pe eq3[Nm/s]はエンジンの出力(Pe eq1>Pe eq2>Pe eq3)で、太実線はエンジン52が同じ出力パワーを得るのに最も燃料消費率が少ない点を結んだ線(最良燃費線)である。この最良燃費線上で目標出力パワーtPeの得られる運転点を探索し、その探索した運転点に対するエンジン回転速度及びエンジントルクを目標回転速度及び目標トルク(つまり目標運転点)として設定する。
Then, the target rotational speed and the target torque of the
ステップ102では、ステップ101で求めている制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}(j=1,2,…,n)のうちの各制駆動力配分を各モータ11〜14が出力したときの総消費電力Pout(j)[W](j=1,2,…,n)を推定する。
In step 102, the set of braking / driving force distributions determined in step 101 {Fx 1 (j), Fx 2 (j), Fx 3 (j), Fx 4 (j)} (j = 1, 2,..., n), the total power consumption Pout (j) [W] (j = 1, 2,..., n) when the
各制駆動力配分を各モータ11〜14が出力したときの総消費電力Pout(j)(各制駆動力配分を実現するための総エネルギー必要量)の推定方法として次に4つを挙げる。
<10>総消費電力Pout(j)の推定方法1(請求項14に記載の発明)
各制駆動力配分を各モータ11〜14が出力したときの総消費電力Pout(j)を、ステップ20で求めている車体速度Vを用いて次の式(111)により算出する。
Next, there are four methods for estimating the total power consumption Pout (j) (the total energy requirement for realizing each braking / driving force distribution) when each motor 11-14 outputs each braking / driving force distribution.
<10>
The total power consumption Pout (j) when each motor 11-14 outputs each braking / driving force distribution is calculated by the following equation (111) using the vehicle body speed V obtained in
Pout(j)=V×(Fx1(j)+Fx2(j)+Fx3(j)+Fx4(j))
…(111)
<11>総消費電力Pout(j)の推定方法2(請求項15に記載の発明)
各制駆動力配分を各モータ11〜14が出力したときの総消費電力Pout(j)を各車輪速度V1,V2,V3,V4を用いて次の式(112)により算出する。
Pout (j) = V × (Fx 1 (j) + Fx 2 (j) + Fx 3 (j) + Fx 4 (j))
... (111)
<11>
Each control total power consumption Pout of the drive force distribution when the
Pout(j)=Fx1(j)×V1+Fx2(j)×V2+Fx3(j)×V3+Fx4(j)×V4
…(112)
<12>総消費電力Pout(j)の推定方法3(請求項16に記載の発明)
各制駆動力配分を各モータ11〜14が出力したときの総消費電力Pout(j)を、各輪1〜4のモータロスPlossi(j)[W](i=1〜4)を考慮した次の式(113)により算出する。
Pout (j) = Fx 1 (j) × V 1 + Fx 2 (j) × V 2 + Fx 3 (j) × V 3 + Fx 4 (j) × V 4
... (112)
<12>
The total power consumption Pout (j) when each motor 11-14 outputs each braking / driving force distribution, considering the motor loss Ploss i (j) [W] (i = 1-4) of each wheel 1-4. It calculates with the following formula | equation (113).
Pout(j)=(Fx1(j)×V1+Ploss1(j))
+(Fx2(j)×V2+Ploss2(j))
+(Fx3(j)×V3+Ploss3(j))
+(Fx4(j)×V4+Ploss4(j)) …(113)
各輪1〜4のモータロスPlossi(j)(各輪1〜4を駆動する時に発生する各エネルギー損失)は、第1実施形態で前述したように各輪1〜4毎に図17のマップを参照して求めればよい。図17のマップは各モータ駆動時の電気的,機械的なロスを制駆動トルクと車輪の回転速度毎に予め求めておいたマップである。
<13>総消費電力Pout(j)の推定方法4
上記<10>〜<12>において総消費電力Pout(j)を求める際に、電力を使用する他の車載機器(エアコン,カーオーディオ,ヘッドライト等),エンジン補機,モータ冷却装置等の消費電力を上乗せする。
Pout (j) = (Fx 1 (j) × V 1 + Ploss 1 (j))
+ (Fx 2 (j) × V 2 + Ploss 2 (j))
+ (Fx 3 (j) × V 3 + Ploss 3 (j))
+ (Fx 4 (j) × V 4 + Ploss 4 (j)) (113)
The motor loss Ploss i (j) (each energy loss generated when driving each
<13>
When calculating the total power consumption Pout (j) in the above <10> to <12>, consumption of other on-vehicle equipment (air conditioner, car audio, headlight, etc.), engine auxiliary equipment, motor cooling device, etc. that uses power Add more power.
これで各制駆動力配分を各モータ11〜14が出力したときの総消費電力Pout(j)の推定方法の説明を終える。 This concludes the description of the method for estimating the total power consumption Pout (j) when each motor 11-14 outputs each braking / driving force distribution.
ステップ103では、ステップ11で検出しているジェネレータ51の発電電力Pgとバッテリ9の放電可能最大電力Pb max[W]とを用いてバッテリ9とジェネレータ51が供給できる電力の上限Ps max[W]を、同じくステップ11で検出しているジェネレータ51の発電電力Pgとバッテリ9の充電可能最大電力Pb min[W]とを用いてバッテリ9とジェネレータ51が供給できる電力の下限Ps min[W]を、つまり次の式(114)及び式(115)により供給電力上限Ps max[W]と供給電力下限Ps min[W]を算出する。
In step 103, the generated power Pg of the
Ps max=Pg+Pb max …(114)
Ps min=Pg+Pb min …(115)
ここで、式(114)の放電可能最大電力Pb max,式(115)の充電可能最大電力Pb minは全て放電側を正の値とする。
Ps max = Pg + Pb max (114)
Ps min = Pg + Pb min (115)
Here, the maximum dischargeable power Pb in the equation (114) max, the maximum chargeable power Pb of the equation (115) All the min values are positive on the discharge side.
式(114)の放電可能最大電力Pb max,式(115)の充電可能最大電力Pb minは、バッテリ9の蓄電量Bcから図21の上段、下段を内容とするテーブルを参照して求める。図21上段のPb m max[W]及び図21下段のPb m min[W]は、バッテリ9が充分に充放電できる状態(本実施形態では蓄電量BSが50%の状態)で、バッテリ9が破損や急速な劣化を起こさない最大放電可能電力及び最大充電可能電力である。
The maximum dischargeable power Pb of the formula (114) max, the maximum chargeable power Pb of the equation (115) The min is obtained from the charged amount Bc of the
図19(B)に移りステップ111では、ステップ102で求めている総消費電力Pout(j)のうち1番目の総消費電力Pout(1)とステップ103で求めている供給電力上限Ps max及び供給電力下限Ps minを比較する。なお、ここでは、1番目の総消費電力Pout(1)を採用しているが、これに限定されるものでなく、2番目,3番目,…の総消費電力Pout(2),Pout(3),…でもかまわない。 19B, in step 111, the first total power consumption Pout (1) out of the total power consumption Pout (j) obtained in step 102 and the supply power upper limit Ps obtained in step 103 are obtained. max and lower limit Ps of power supply Compare min. Here, the first total power consumption Pout (1) is adopted, but the present invention is not limited to this, and the second, third,... Total power consumption Pout (2), Pout (3) ), ... But it doesn't matter.
1番目の総消費電力Pout(1)が供給電力上限Ps max以下でかつ供給電力下限Ps min以上(Ps min≦Pout(1)≦Ps max)であればステップ112に進んでフラグfp=1とする。一方、1番目の総消費電力Pout(1)が供給電力上限Ps maxを超えているかまたは供給電力下限Ps min未満であるときにはステップ111よりステップ113に進んでフラグfp=0とする。 The first total power consumption Pout (1) is the supply power upper limit Ps. Less than max and supply power lower limit Ps min or more (Ps min ≦ Pout (1) ≦ Ps max), the routine proceeds to step 112 where the flag fp = 1 is set. On the other hand, the first total power consumption Pout (1) is the supply power upper limit Ps. exceeds max or supply power lower limit Ps When it is less than min, the routine proceeds from step 111 to step 113, where the flag fp = 0 is set.
ここで、1番目の総消費電力Pout(1)は、ステップ70で設定している駆動力配分{Fx1,Fx2,Fx3,Fx4}を各モータ11〜14で出力するのに必要な総消費電力である。すなわち、ステップ111はステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}が、ジェネレータ51とバッテリ9が供給可能な電力で実現可能か否かでフラグfpの値を変えるものである。
Here, the first total power consumption Pout (1) is necessary to output the driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } set in step 70 by the
言い換えると、ステップ111は、1番目の総消費電力Pout(1)(制駆動力配分の集合のうちの各制駆動力配分を実現するための総エネルギー必要量)が車両の出力可能な電力(車両の出力可能なエネルギー)の範囲内にあるか否かを判定する部分である。すなわち、1番目の総消費電力Pout(1)が供給電力上限Ps max以下でかつ供給電力下限Ps min以上であれば車両の出力可能な電力の範囲内にあると、これに対して、1番目の総消費電力Pout(1)が供給電力上限Ps maxを超えているかまたは供給電力下限Ps min未満であるときには車両の出力可能な電力の範囲内にないと判定している。 In other words, in step 111, the first total power consumption Pout (1) (total energy requirement for realizing each braking / driving force distribution in the braking / driving force distribution set) is the power that can be output by the vehicle ( This is a part for determining whether or not the energy is within the range of the energy that can be output from the vehicle. That is, the first total power consumption Pout (1) is the supply power upper limit Ps. Less than max and supply power lower limit Ps If it is equal to or greater than min, and within the range of power that can be output by the vehicle, the first total power consumption Pout (1) is the supply power upper limit Ps. exceeds max or supply power lower limit Ps When it is less than min, it is determined that it is not within the range of electric power that the vehicle can output.
ステップ114では、必要制駆動力補正量ΔFsi[N]の絶対値|ΔFsi|とゼロを比較する。必要制駆動力補正量の絶対値|ΔFsi|がゼロより大きい車輪が一輪もない場合にはステップ115に進んでフラグft=1とする。一方、必要制駆動力補正量の絶対値|ΔFsi|がゼロより大きい車輪が一輪以上ある場合にはステップ114よりステップ116に進んでフラグft=0とする。 In step 114, the absolute value | ΔFs i | of the required braking / driving force correction amount ΔFs i [N] is compared with zero. If no wheel has an absolute value | ΔFs i | of the required braking / driving force correction amount greater than zero, the routine proceeds to step 115 where the flag ft = 1 is set. On the other hand, if there is one or more wheels whose absolute value | ΔFs i | of the required braking / driving force correction amount is greater than zero, the routine proceeds from step 114 to step 116, where the flag ft = 0 is set.
必要制駆動力補正量ΔFsi(i=1〜4)は、第1実施形態で前述したように制駆動力Fx1,Fx2,Fx3,Fx4から各輪1〜4毎に次のように設定している。すなわち、各輪1〜4について制駆動力Fxiが制駆動力下限Fminiを超えかつ制駆動力上限Fmaxi未満(Fmini<Fxi<Fmaxi)の領域にあれば必要制駆動力補正量ΔFsi[N]=0である。各輪1〜4について制駆動力Fxiが制駆動力上限Fmaxi以上(Fxi≧Fmaxi)のとき必要制駆動力補正量ΔFsi[N]=Fmaxi−Fxiであり、各輪1〜4について制駆動力Fxiが制駆動力下限Fmini以下(Fxi≦Fmini)のとき必要制駆動力補正量ΔFsi[N]=Fmini−Fxiである。
The necessary braking / driving force correction amount ΔFs i (i = 1 to 4) is determined from the braking / driving forces Fx 1 , Fx 2 , Fx 3 , and Fx 4 for each
従って、必要制駆動力補正量ΔFsiがゼロでない車輪がある(フラグft=0)ということは、ステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}を出力できない車輪が存在することを意味している。 Therefore, there is a wheel where the required braking / driving force correction amount ΔFs i is not zero (flag ft = 0), which means that the braking / driving force distribution set in step 70 {Fx 1 , Fx 2 , Fx 3 , Fx 4 } This means that there are wheels that cannot output.
ステップ117では、フラグfpとフラグftが共に1であるか否かをみる。フラグfp=1かつフラグft=1の場合、つまりステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}がジェネレータ51とバッテリ9が供給可能な電力で実現可能でありかつステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}を出力できない車輪が存在しない場合にはステップ130以降に進むことなく、ステップ300に進み、ステップ101で求めている制駆動力配分の集合のうちの一つの制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択し、その選択した制駆動力Fx1,Fx2,Fx3,Fx4にそれぞれ各輪1〜4の半径Rを乗じた値、即ちトルク指令値(Fx1×R,Fx2×R,Fx3×R,Fx4×R)が得られるようにモータ11〜14を制御する。
In step 117, it is determined whether or not both the flag fp and the flag ft are 1. When the flag fp = 1 and the flag ft = 1, that is, the braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } set in step 70 is realized with the power that can be supplied by the
一方、フラグfp=0またはフラグft=0の場合、つまりステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}がジェネレータ51とバッテリ9が供給可能な電力で実現不可能である場合またはステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}を出力できない車輪が存在する場合にはステップ117よりステップ130に進み、ステップ101で求めている制駆動力配分の集合{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}の中で、各輪1〜4の制駆動力制限を満たし、かつジェネレータ51とバッテリ9が供給可能な電力で実現可能なものを選び出す。
On the other hand, when the flag fp = 0 or the flag ft = 0, that is, the braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } set in step 70 is the power that can be supplied by the
ここではまず、ステップ101で求めている制駆動力配分の集合のうちの各制駆動力配分を、ステップ103で求めている供給電力上限Ps max,供給電力下限Ps min及びステップ110で求めている制駆動力上限Fmaxi及び制駆動力下限Fminiを用いて各輪1〜4毎に制限する方法について次に説明する。
Here, first, each braking / driving force distribution in the set of braking / driving force distributions obtained in
まず、図22に示すような配列St2を考える。配列St2の列数はステップ101で求めている制駆動力配分の集合の数nと同じ数で、行数は7である。
First, consider an array St2 as shown in FIG. The number of columns in the array St2 is the same as the number n of sets of braking / driving force distributions obtained in
まず整数j(j=1,2,…,n)を1から順にnまで変化させつつ、j番目の総消費電力Pout(j)が供給電力下限Ps min以上かつ供給電力上限Ps max以下(Ps min≦Pout(j)≦Ps max)にある場合に、対応するjの列の第2行目(j×2)に1を、これに対してj番目の総消費電力Pout(j)が供給電力下限Ps min未満または供給電力上限Ps maxを超えている場合に、対応するjの列の2行目にゼロを設定する。例えば図22では整数j=1,2,3,4のとき1が、整数j=5,6のときゼロが設定されているので、1,2,3,4番目の総消費電力Pout(1),Pout(2),Pout(3),Pout(4)は供給電力下限Ps min以上かつ供給電力上限Ps max以下にあるのに対して、5,6番目の総消費電力Pout(5),Pout(6)は供給電力下限Ps min以上かつ供給電力上限Ps max以下にないことがわかる。 First, while changing the integer j (j = 1, 2,..., N) from 1 to n, the j-th total power consumption Pout (j) min or more and upper limit of power supply Ps max or less (Ps min ≦ Pout (j) ≦ Ps 1) in the second row (j × 2) of the corresponding column j, and the j-th total power consumption Pout (j) is the supply power lower limit Ps Less than min or upper limit of power supply Ps If it exceeds max, zero is set in the second row of the corresponding j column. For example, in FIG. 22, 1 is set when the integer j = 1, 2, 3, 4 and 0 is set when the integer j = 5, 6, so the first, second, third, and fourth total power consumption Pout (1 ), Pout (2), Pout (3), and Pout (4) are lower limits of supply power Ps. min or more and upper limit of power supply Ps The fifth and sixth total power consumptions Pout (5) and Pout (6) are lower than the power supply lower limit Ps while being below max. min or more and upper limit of power supply Ps It turns out that it is not below max.
次にまた整数jを1から順にnまで変化させつつ、j番目の左前輪1の制駆動力であるFx1(j)が制駆動力下限Fmin1以上かつ制駆動力上限Fmax1以下(Fmin1≦Fx1(j)≦Fmax1)にある場合に、対応するjの列の第3行目(j×3)に1を、これに対してj番目の左前輪1の制駆動力であるFx1(j)が制駆動力下限Fmin1未満である場合または制駆動力上限Fmax1を超えている場合に、対応するjの列の第3行目にゼロを設定する。これは、各駆動力配分がそれぞれ左前輪1の駆動力制限を満たすか否かをチェックするものである。例えば図22では整数j=1,2,3のとき1が、整数j=4,5,6のときゼロが設定されているので、1,2,3番目の左前輪1の制駆動力Fx1(1),Fx1(2),Fx1(3)は制駆動力下限Fmin1以上かつ制駆動力上限Fmax1以下にあるのに対して、4,5,6番目の左前輪1の制駆動力Fx1(4),Fx1(5),Fx1(6)は制駆動力下限Fmin1以上かつ制駆動力上限Fmax1以下にないことがわかる。
Next, while changing the integer j sequentially from 1 to n, the braking / driving force Fx 1 (j) of the j-th left
同様にして、右前輪2,左後輪3,右後輪4についても制駆動力制限を満たす場合には1を、これに対して制駆動力制限を満たさない場合にはゼロを上記の配列St2の第4、第5、第6行目に設定する(右前輪2の場合はj×4に、左後輪3の場合はj×5に、右後輪4の場合はj×6に設定する)。
Similarly, for the right
そして最後に、また整数jを1から順にnまで変化させつつ、各j列の第2行目〜第6行目までの要素のANDを採り第7行目にその結果を記録する。この配列St2の第7行目が1である列のjは、対応する制駆動力配分{Fx1,Fx2,Fx3,Fx4}が各輪1〜4の制駆動力制限を満たし、かつジェネレータ51とバッテリ9が供給可能な電力で実現可能であることを示している。例えば図22では整数j=1〜6の中では2,3番目の制駆動力配分{Fx1,Fx2,Fx3,Fx4}のみが各輪1〜4の制駆動力制限を満たし、かつジェネレータ51とバッテリ9が供給可能な電力で実現可能であることを示している。
Finally, while changing the integer j from 1 to n in order, AND of the elements from the second row to the sixth row of each j column is taken and the result is recorded in the seventh row. In the column j in which the seventh row of this array St2 is 1, the corresponding braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } satisfies the braking / driving force limitation of each
また、ステップ130では配列St2の第7行目の和を配列数Stnum2に設定する。従って、配列数Stnum2は実現可能な制駆動力配分の数を示すこととなる。 In step 130, the sum of the seventh row of the array St2 is set to the array number Stnum2. Therefore, the arrangement number Stnum2 indicates the number of braking / driving force distributions that can be realized.
ステップ140では、この配列数Stnum2とゼロを比較する。配列数Stnum2≠0の場合(つまりステップ101で求めている制駆動力配分の集合の中に実現可能な制駆動力配分が存在する)にはステップ150に進み、ステップ130で求めている実現可能な制駆動力配分の集合の中から、一つの制駆動力配分{Fx1,Fx2,Fx3,Fx4}を選択し、ステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}にこの選択した制駆動力配分{Fx1,Fx2,Fx3,Fx4}を上書きする。
In step 140, this array number Stnum2 is compared with zero. If the number of arrays Stnum2 ≠ 0 (that is, there is a realizable braking / driving force distribution in the set of braking / driving force distributions determined in step 101), the process proceeds to step 150, and the realization determined in step 130 is possible. From one set of braking / driving force distributions, one braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } is selected, and the braking / driving force distribution {Fx 1 , Fx set in step 70 is selected. 2, Fx 3, Fx 4} to the selected braking-driving
ステップ150におけるこの一つの制駆動力配分の選択方法として次に3つを挙げる。
<14>一つの制駆動力配分の選択1
配列St2の第7行目が1となるjに対応する制駆動力配分{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}と、ステップ70で設定している各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4とから次の式(116)の評価関数J(j)をそれぞれ計算し、それら計算した評価関数J(j)のうち、評価関数J(j)を最小とする制駆動力配分を選択する。
Three methods for selecting this one braking / driving force distribution in
<14> Selection of one braking / driving
Braking / driving force distribution {Fx 1 (j), Fx 2 (j), Fx 3 (j), Fx 4 (j)} corresponding to j in which the 7th row of array St2 is 1 is set in step 70 From the braking / driving forces Fx 1 , Fx 2 , Fx 3 , and Fx 4 of the wheels 1 to 4, the evaluation function J (j) of the following equation (116) is calculated, and the calculated evaluation function J ( Among j), the braking / driving force distribution that minimizes the evaluation function J (j) is selected.
J(j)=(Fx1−Fx1(j))2+(Fx2−Fx2(j))2
+(Fx3−Fx3(j))2+(Fx4−Fx4(j))2
…(116)
<15>一つの制駆動力配分の選択2
配列St2の第7行目が1となるjに対応する総消費電力Pout(j)のうち、総消費電力Pout(j)が最小となる制駆動力配分を選ぶ。
<16>一つの制駆動力配分の選択3
配列St2の第7行目が1となるjに対応する制駆動力配分の中で、各輪1〜4のタイヤ横力と駆動力の合力であるタイヤ力を輪荷重W1,W2,W3,W4で除した値であるタイヤ摩擦円使用率H1(j),H2(j),H3(j),H4(j)の最大値と最小値の差が最も小さくなる制駆動力配分を選択する。選択の方法は次の通り。
J (j) = (Fx 1 −Fx 1 (j)) 2 + (Fx 2 −Fx 2 (j)) 2
+ (Fx 3 −Fx 3 (j)) 2 + (Fx 4 −Fx 4 (j)) 2
... (116)
<15> Selection of one braking / driving
Of the total power consumption Pout (j) corresponding to j in which the seventh row of the array St2 is 1, the braking / driving force distribution that minimizes the total power consumption Pout (j) is selected.
<16> Selection of one braking / driving
In the braking / driving force distribution corresponding to j in which the 7th row of the array St2 is 1, the tire force which is the resultant force of the tire lateral force and the driving force of each of the
まず、配列St2の第7行目が1となるjに対応する制駆動力配分{Fx1(j),Fx2(j),Fx3(j),Fx4(j)}を用いて、各輪1〜4のタイヤ摩擦円使用率H1(j),H2(j),H3(j),H4(j)を次の式(117a)〜式(117d)により算出する。
First, using the braking / driving force distribution {Fx 1 (j), Fx 2 (j), Fx 3 (j), Fx 4 (j)} corresponding to j in which the seventh row of the array St2 is 1, The tire friction circle usage rates H 1 (j), H 2 (j), H 3 (j), and H 4 (j) of each
H1(j)=(Fx1(j)2+Fy1(j)2)1/2/W1 …(117a)
H2(j)=(Fx2(j)2+Fy2(j)2)1/2/W2 …(117b)
H3(j)=(Fx3(j)2+Fy3(j)2)1/2/W3 …(117c)
H4(j)=(Fx4(j)2+Fy4(j)2)1/2/W4 …(117d)
そして次に、各j毎にH1(j),H2(j),H3(j),H4(j)のうちの最大値Hmax(j)と最小値Hmin(j)とを求め、その差であるΔH(j)=Hmax(j)−Hmin(j)を計算し、この差ΔH(j)が最小となるjに対応する制駆動力配分を選択する。
H 1 (j) = (Fx 1 (j) 2 + Fy 1 (j) 2 ) 1/2 / W 1 (117a)
H 2 (j) = (Fx 2 (j) 2 + Fy 2 (j) 2 ) 1/2 / W 2 (117b)
H 3 (j) = (Fx 3 (j) 2 + Fy 3 (j) 2 ) 1/2 / W 3 (117c)
H 4 (j) = (Fx 4 (j) 2 + Fy 4 (j) 2 ) 1/2 / W 4 (117d)
Then, for each j, the maximum value Hmax (j) and the minimum value Hmin (j) of H 1 (j), H 2 (j), H 3 (j), and H 4 (j) are obtained. The difference ΔH (j) = Hmax (j) −Hmin (j) is calculated, and the braking / driving force distribution corresponding to j that minimizes the difference ΔH (j) is selected.
式(117a)〜式(117d)の各輪1〜4のタイヤ横力Fy1(j),Fy2(j),Fy3(j),Fy4(j)の算出方法の一例は第1実施形態においてステップ100で説明済みであるので、ここでの説明は省略する。
An example of a method for calculating the tire lateral forces Fy 1 (j), Fy 2 (j), Fy 3 (j), and Fy 4 (j) of the
これでステップ150における一つの制駆動力配分の選択方法の説明を終える。
This completes the description of one braking / driving force distribution selection method in
一方、ステップ140で配列数Stnum2=0の場合(つまりステップ101で求めている制駆動力配分の集合の中に実現可能な制駆動力配分が存在しない)には制駆動力配分を再設定するためステップ211以降に進む。 On the other hand, if the number of arrays Stnum2 = 0 in step 140 (that is, there is no realizable braking / driving force distribution in the set of braking / driving force distributions determined in step 101), the braking / driving force distribution is reset. Therefore, the process proceeds to step 211 and thereafter.
ステップ211では、フラグftとゼロを比較する。フラグft=0の場合(つまりステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}を出力できない車輪が存在する)にはステップ221に進み、ステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}を、ステップ110で求めている制駆動力上限Fmaxiと制駆動力下限Fminiとの間に収まるように再設定する。この制駆動力配分の再設定方法について次に述べる。
In
まず、各車輪1〜4についてステップ110で求めている制駆動力上限Fmaxi、制駆動力下限Fminiを用い、制駆動力Fxiの値(正、負、ゼロ)に応じて次の式(118)〜式(120)により減少比率Qi[無名数]を算出する。
First, using the braking / driving force upper limit Fmax i and the braking / driving force lower limit Fmin i obtained in step 110 for each of the
制駆動力Fxi>0の場合:Qi=Fmaxi/Fxi …(118)
制駆動力Fxi<0の場合:Qi=Fmini/Fxi …(119)
制駆動力Fxi=0の場合:Qi=1 …(120)
ここで、制駆動力上限Fmaxi≧0、制駆動力下限Fmini≦0であるので、制駆動力Fxi>0の場合、制駆動力Fxi<0の場合とも減少比率Qiは正またはゼロの値となる。
When braking / driving force Fx i > 0: Q i = Fmax i / Fx i (118)
When braking / driving force Fx i <0: Q i = Fmin i / Fx i (119)
When braking / driving force Fx i = 0: Q i = 1 (120)
Here, since the braking / driving force upper limit Fmax i ≧ 0 and the braking / driving force lower limit Fmin i ≦ 0, when the braking / driving force Fx i > 0, the reduction ratio Q i is positive even when the braking / driving force Fx i <0. Or a zero value.
次に、各輪1〜4の減少比率Qiの中で最も小さい値を改めて減少比率Qとして採用し、ステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}を次の式(121)〜式(124)により再設定する。
Next, the smallest value among the reduction ratios Q i of the
Fx1←Q×Fx1 …(121)
Fx2←Q×Fx2 …(122)
Fx3←Q×Fx3 …(123)
Fx4←Q×Fx4 …(124)
式(121)〜式(124)は、各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4の絶対値を同じ割合(減少比率Q)で減少させることにより、各輪1〜4の制駆動力制限を満たすようにしたものである。
Fx 1 ← Q × Fx 1 (121)
Fx 2 ← Q × Fx 2 (122)
Fx 3 ← Q × Fx 3 (123)
Fx 4 ← Q × Fx 4 (124)
Expressions (121) to (124) are obtained by reducing the absolute values of the braking / driving forces Fx 1 , Fx 2 , Fx 3 , and Fx 4 of the
ステップ231では、ステップ221で再設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}を各モータ11〜14で出力する際に必要とされる総消費電力Pout’を、上記のステップ102と同じ手順で推定する。
In step 231, the total power consumption Pout ′ required when the
ステップ241では、この推定した総消費電力Pout’が供給電力上限Ps max以下にありかつ供給電力下限Ps min以上にある(Ps min≦Pout’≦Ps max)か否か、つまりステップ221で再設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}がバッテリ9とジェネレータ51で供給できる電力で実現可能か否かをチェックする。総消費電力Pout’が供給電力上限Ps max以下にありかつ供給電力下限Ps min以上にある場合には、ステップ221で再設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}がバッテリ9とジェネレータ51で供給できる電力で実現可能であると判断しステップ300に進んでステップ300の操作を実行する。
In
一方、総消費電力Pout’が供給電力上限Ps maxを超えている場合(Pout’>Ps max)または供給電力下限Ps min未満である場合(Pout’<Ps min)には、ステップ221で再設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}はバッテリ9とジェネレータ51で供給できる電力では実現不可能であると判断し、上記の総消費電力Pout(1)にこの総消費電力Pout’を再設定してステップ251に進む。
On the other hand, the total power consumption Pout ′ is the supply power upper limit Ps. When max is exceeded (Pout ′> Ps max) or supply power lower limit Ps If less than min (Pout ′ <Ps min), it is determined that the braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } reset at step 221 cannot be realized with the power that can be supplied by the
また、ステップ211でフラグft=1の場合(つまりステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}を出力できない車輪が存在しない)にもステップ221、231、241を飛ばしてステップ251に進む。
Further, when the flag ft is 1 in step 211 (that is, there is no wheel that cannot output the braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } set in step 70), step 221,
ステップ251では、ステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}がバッテリ9とジェネレータ51が供給可能な電力制限内に収まるように、ステップ103で求めている供給電力上限Ps max,供給電力下限Ps minを用いて、ステップ70で設定している制駆動力配分{Fx1,Fx2,Fx3,Fx4}を次の式(125a)〜式(125d)、式(126a)〜式(126d)により再設定する。
In step 251, the braking / driving force distribution {Fx 1 , Fx 2 , Fx 3 , Fx 4 } set in step 70 is obtained in step 103 so that it falls within the power limit that can be supplied by the
総消費電力Pout(1)>供給電力上限Ps maxの場合:
Fx1←(Ps max/Pout(1))×Fx1×τ …(125a)
Fx2←(Ps max/Pout(1))×Fx2×τ …(125b)
Fx3←(Ps max/Pout(1))×Fx3×τ …(125c)
Fx4←(Ps max/Pout(1))×Fx4×τ …(125d)
総消費電力Pout(1)<供給電力下限Ps minの場合:
Fx1←(Ps min/Pout(1))×Fx1×τ …(126a)
Fx2←(Ps min/Pout(1))×Fx2×τ …(126b)
Fx3←(Ps min/Pout(1))×Fx3×τ …(126c)
Fx4←(Ps min/Pout(1))×Fx4×τ …(126d)
式(125a)〜式(125d)、式(126a)〜式(126d)において、係数τを省略した次の式は、バッテリ9とジェネレータ51のパワー不足を解消するように、各輪1〜4の制駆動力Fx1,Fx2,Fx3,Fx4を、(不足したパワー)÷(回転速度)だけ単純に減少させたものである。
Total power consumption Pout (1)> Supply power upper limit Ps For max:
Fx 1 ← (Ps max / Pout (1)) × Fx 1 × τ (125a)
Fx 2 ← (Ps max / Pout (1)) × Fx 2 × τ (125b)
Fx 3 ← (Ps max / Pout (1)) × Fx 3 × τ (125c)
Fx 4 ← (Ps max / Pout (1)) × Fx 4 × τ (125d)
Total power consumption Pout (1) <Supply power lower limit Ps For min:
Fx 1 ← (Ps min / Pout (1)) × Fx 1 × τ (126a)
Fx 2 ← (Ps min / Pout (1)) × Fx 2 × τ (126b)
Fx 3 ← (Ps min / Pout (1)) × Fx 3 × τ (126c)
Fx 4 ← (Ps min / Pout (1)) × Fx 4 × τ (126d)
In the equations (125a) to (125d) and (126a) to (126d), the following equation in which the coefficient τ is omitted is as follows. The braking / driving forces Fx 1 , Fx 2 , Fx 3 , and Fx 4 are simply reduced by (insufficient power) / (rotational speed).
総消費電力Pout(1)>供給電力上限Ps maxの場合:
Fx1←(Ps max/Pout(1))×Fx1 …(補25a)
Fx2←(Ps max/Pout(1))×Fx2 …(補25b)
Fx3←(Ps max/Pout(1))×Fx3 …(補25c)
Fx4←(Ps max/Pout(1))×Fx4 …(補25d)
総消費電力Pout(1)<供給電力下限Ps minの場合:
Fx1←(Ps min/Pout(1))×Fx1 …(補26a)
Fx2←(Ps min/Pout(1))×Fx2 …(補26b)
Fx3←(Ps min/Pout(1))×Fx3 …(補26c)
Fx4←(Ps min/Pout(1))×Fx4 …(補26d)
式(補25a)〜式(補25d)、式(補26a)〜式(補26d)によれば、各モータ11〜14の運転点(トルク,回転速度)の変化に伴って各モータ11〜14の損失が増加し、パワー不足が解消されない場合がある。そこでこのような現象を回避するために、係数τ(0<τ<1)を用いた上記の式(125a)〜式(125d)、式(126a)〜式(126d)により制駆動力配分{Fx1,Fx2,Fx3,Fx4}を再設定するものである。つまり、係数τは、制駆動力変化に伴い各モータ損失が増加する場合等を考慮したもので、本実施形態では係数τとして0.8を設定している。
Total power consumption Pout (1)> Supply power upper limit Ps For max:
Fx 1 ← (Ps max / Pout (1)) × Fx 1 (Supplement 25a)
Fx 2 ← (Ps max / Pout (1)) × Fx 2 (Supplement 25b)
Fx 3 ← (Ps max / Pout (1)) × Fx 3 (Supplement 25c)
Fx 4 ← (Ps max / Pout (1)) × Fx 4 (Supplement 25d)
Total power consumption Pout (1) <Supply power lower limit Ps For min:
Fx 1 ← (Ps min / Pout (1)) × Fx 1 (Supplement 26a)
Fx 2 ← (Ps min / Pout (1)) × Fx 2 (Supplement 26b)
Fx 3 ← (Ps min / Pout (1)) × Fx 3 (Supplement 26c)
Fx 4 ← (Ps min / Pout (1)) × Fx 4 (Supplement 26d)
According to the equations (complement 25a) to the equation (complement 25d) and the equations (complement 26a) to the equation (complement 26d), the
ステップ300では、このようにして再設定した制駆動力Fx1,Fx2,Fx3,Fx4にそれぞれ各輪1〜4の半径Rを乗じた値、即ちトルク指令値(Fx1×R,Fx2×R,Fx3×R,Fx4×R)が得られるようにモータ11〜14を制御する。
In step 300, this way the longitudinal force Fx 1 in which and reconfiguring, Fx 2, Fx 3, Fx 4 to a value obtained by respectively multiplying the radius R of the
第2実施形態(請求項14、15、16に記載の発明)によれば、ステップ101で求めている制駆動力配分の集合のうちの各制駆動力配分を各モータ11〜14で出力するのに必要な総消費電力Pout(j)(総エネルギー必要量)を推定し、この総消費電力Pout(j)がジェネレータ51とバッテリ9が供給可能な電力(車両の出力可能なエネルギー)の範囲内(Ps min≦Pout(1)≦Ps max)にある制駆動力配分を選択するので(図19(B)のステップ111、112、114、115、117、300参照)、ステップ101で求めている制駆動力配分の集合のうちの各制駆動力配分{Fx1,Fx2,Fx3,Fx4}を各モータ11〜14で出力するのに必要な総消費電力(Pout(j))が、ジェネレータ51とバッテリ9が供給可能な電力(車両の発電能力)を超えることがなく、車両の操縦安定性をより向上させることができる。
According to the second embodiment (the invention described in
実施形態では、前輪転舵機構に加えて後輪転舵機構を有する四輪独立駆動車で説明したが、後輪転舵機構を有さず前輪のみ転舵する四輪独立駆動車に対しても本発明を適用することができる。この四輪独立駆動車に本発明を適用する場合には、例えば図9(A),図9(B)において後輪3,4の舵角δ3,δ4の値を定数(通常はδ3,δ4=0)とすればよい。
In the embodiment, the four-wheel independent drive vehicle having the rear wheel steering mechanism in addition to the front wheel steering mechanism has been described, but the present invention is also applied to a four-wheel independent drive vehicle that does not have the rear wheel steering mechanism and steers only the front wheels. The invention can be applied. When the present invention is applied to this four-wheel independent drive vehicle, for example, in FIGS. 9A and 9B, the values of the steering angles δ 3 and δ 4 of the
また、ジェネレータ51とエンジン52の代りに、燃料電池を備える四輪独立駆動車への適用も可能である。この燃料電池を備える四輪独立駆動車に本発明を適用する場合には、図19(A)、図19(B)のステップ103においてジェネレータ51の発電電力Pgに代えて、この燃料電池の発電電力を用いるようにすればよい。
Further, instead of the
実施形態では、各輪1〜4を駆動または制動するアクチュエータがモータ11〜14である場合で説明したが、各輪1〜4を制動する際に、各モータ11〜14に加えてメカブレーキを併用するものにも本発明を適用することができる。
In the embodiment, the description has been given of the case where the actuators that drive or brake the
請求項1に記載の目標車両挙動決定手段の機能は図9(A)のステップ100により、制駆動力配分集合演算手段の機能は図9(A)のステップ101、図13、図15により、制駆動力配分選択手段の機能は図9(B)のステップ150によりそれぞれ果たされている。
The function of the target vehicle behavior determining means according to
1〜4 車輪
8 コントローラ
11〜14 モータ(アクチュエータ)
1 to 4
Claims (16)
車両の平面運動での目標となる車両挙動を決定する目標車両挙動決定手段と、
この目標となる車両挙動を実現する四輪の制駆動力配分の集合を演算する制駆動力配分集合演算手段と、
この演算された制駆動力配分の集合の中から一つの制駆動力配分を選択する制駆動力配分選択手段と
を備えることを特徴とする四輪独立駆動車の制駆動力配分装置。 In vehicles that can drive all four wheels independently,
A target vehicle behavior determining means for determining a target vehicle behavior in the planar movement of the vehicle;
Braking / driving force distribution set calculating means for calculating a set of four-wheel braking / driving force distribution for realizing the target vehicle behavior;
A braking / driving force distribution device for a four-wheel independent drive vehicle, comprising: a braking / driving force distribution selecting unit that selects one braking / driving force distribution from the calculated set of braking / driving force distributions.
各輪の横すべり角を推定する横すべり角推定手段と
を備え、
前記制駆動力配分集合演算手段は、前記目標となる車両挙動を実現する四輪の制駆動力配分の集合を、各輪の輪荷重と、各輪の横すべり角とに基づいて演算することを特徴とする請求項1に記載の四輪独立駆動車の制駆動力配分装置。 Wheel load estimating means for estimating the wheel load of each wheel;
A side slip angle estimating means for estimating a side slip angle of each wheel, and
The braking / driving force distribution set calculating means calculates a set of four wheel braking / driving force distributions that realize the target vehicle behavior based on the wheel load of each wheel and the side slip angle of each wheel. The braking / driving force distribution device for a four-wheel independent drive vehicle according to claim 1,
を備え、
前記制駆動力配分集合演算手段は、前記目標となる車両挙動を実現する制駆動力配分の集合を、各輪の制駆動力変化に対するタイヤ横力の感度に基づいて演算することを特徴とする請求項1または2に記載の四輪独立駆動車の制駆動力配分装置。 Tire lateral force sensitivity estimation means for estimating the sensitivity of the tire lateral force with respect to the braking / driving force change of each wheel,
The braking / driving force distribution set calculating means calculates a set of braking / driving force distributions that realizes the target vehicle behavior based on the sensitivity of the tire lateral force with respect to a change in braking / driving force of each wheel. The braking / driving force distribution device for a four-wheel independent drive vehicle according to claim 1 or 2.
各輪の制駆動力変化に対するタイヤ横力の感度k1,k2,k3,k4を推定するタイヤ横力感度推定手段と
を備え、
前記制駆動力配分集合演算手段は、前記目標となる車両挙動を実現する制駆動力配分の集合を、車両の平面運動での車両挙動の変化をゼロとする各輪の制駆動力補正量ΔFx1,ΔFx2,ΔFx3,ΔFx4の比を求める次の式(A)〜式(E)に基づいて演算することを特徴とする請求項1または2に記載の四輪独立駆動車の制駆動力配分装置。
ΔFx1:ΔFx2:ΔFx3:ΔFx4
={(Lt/Ll)(h4−h2)+h2(h4−h3)}
/(cosδ1−k1sinδ1)
:{−(Lt/Ll)(h3−h1)−h1(h4−h3)}
/(cosδ2−k2sinδ2)
:{−(Lt/Ll)(h4−h2)−h4(h2−h1)}
/(cosδ3−k3sinδ3)
:{(Lt/Ll)(h3−h1)+h3(h2−h1)}
/(cosδ4−k4sinδ4)…(A)
ただし、h1=(sinδ1+k1cosδ1)/(cosδ1−k1sinδ1)…(B)
h2=(sinδ2+k2cosδ2)/(cosδ2−k2sinδ2)…(C)
h3=(sinδ3+k3cosδ3)/(cosδ3−k3sinδ3)…(D)
h4=(sinδ4+k4cosδ4)/(cosδ4−k4sinδ4)…(E)
h1,h2,h3,h4:係数、
Lt:トレッド長さ、
Ll:ホイールベース長さ、 Rudder angle detecting means for detecting the rudder angles δ 1 , δ 2 , δ 3 , δ 4 of each wheel;
Tire lateral force sensitivity estimating means for estimating tire lateral force sensitivity k 1 , k 2 , k 3 , k 4 with respect to the braking / driving force change of each wheel, and
The braking / driving force distribution set calculating means sets the braking / driving force distribution set for realizing the target vehicle behavior as a braking / driving force correction amount ΔFx for each wheel with zero change in the vehicle behavior due to the planar movement of the vehicle. 3. The control of the four-wheel independent drive vehicle according to claim 1, wherein the calculation is performed based on the following formulas (A) to (E) for obtaining a ratio of 1 , ΔFx 2 , ΔFx 3 , and ΔFx 4. Driving force distribution device.
ΔFx 1 : ΔFx 2 : ΔFx 3 : ΔFx 4
= {(L t / L l ) (h 4 −h 2 ) + h 2 (h 4 −h 3 )}
/ (Cosδ 1 −k 1 sinδ 1 )
: {- (L t / L l) (h 3 -h 1) -h 1 (h 4 -h 3)}
/ (Cosδ 2 −k 2 sinδ 2 )
: {- (L t / L l) (h 4 -h 2) -h 4 (h 2 -h 1)}
/ (Cosδ 3 −k 3 sinδ 3 )
: {(L t / L l ) (h 3 −h 1 ) + h 3 (h 2 −h 1 )}
/ (Cosδ 4 −k 4 sinδ 4 ) (A)
However, h 1 = (sin δ 1 + k 1 cos δ 1 ) / (cos δ 1 −k 1 sin δ 1 ) (B)
h 2 = (sin δ 2 + k 2 cos δ 2 ) / (cos δ 2 −k 2 sin δ 2 ) (C)
h 3 = (sin δ 3 + k 3 cos δ 3 ) / (cos δ 3 −k 3 sin δ 3 ) (D)
h 4 = (sin δ 4 + k 4 cos δ 4 ) / (cos δ 4 −k 4 sin δ 4 ) (E)
h 1 , h 2 , h 3 , h 4 : coefficients,
L t : tread length,
L l : wheelbase length,
前記制駆動力配分集合演算手段により演算された制駆動力配分の集合を、この各輪が路面に伝達可能な制制駆動力上限に制限することを特徴とする請求項1から4までのいずれか一つに記載の四輪独立駆動車の制駆動力配分装置。 Road surface transmission braking / driving force upper limit estimating means for estimating the braking / driving force that each wheel can transmit to the road surface,
5. A braking / driving force distribution set calculated by the braking / driving force distribution set calculating means is limited to an upper limit of braking / driving force that each wheel can transmit to a road surface. The braking / driving force distribution device for a four-wheel independent drive vehicle according to claim 1.
前記制駆動力配分集合演算手段により演算された制駆動力配分の集合を、この各輪が路面に伝達可能な制制駆動力下限に制限することを特徴とする請求項1から4までのいずれか一つに記載の四輪独立駆動車の制駆動力配分装置。 A braking / driving force lower limit estimating means for estimating a braking / driving force lower limit that each wheel can transmit to the road surface,
5. The braking / driving force distribution set calculated by the braking / driving force distribution set calculating means is limited to a lower limit of braking / driving force that each wheel can transmit to the road surface. The braking / driving force distribution device for a four-wheel independent drive vehicle according to claim 1.
前記制駆動力配分集合演算手段により演算された制駆動力配分の集合を、この各輪を駆動または制動するアクチュエータが出力可能な制駆動力上限に制限することを特徴とする請求項1から4までのいずれか一つに記載の四輪独立駆動車の制駆動力配分装置。 A braking / driving force upper limit estimating means for estimating an upper limit of braking / driving force that can be output by an actuator that drives or brakes each wheel;
5. The braking / driving force distribution set calculated by the braking / driving force distribution set calculating means is limited to an upper limit of braking / driving force that can be output by an actuator that drives or brakes each wheel. The braking / driving force distribution device for a four-wheel independent drive vehicle according to any one of the above.
前記制駆動力配分集合演算手段により演算された制駆動力配分の集合を、この各輪を駆動または制動するアクチュエータが出力可能な制駆動力下限に制限することを特徴とする請求項1から4までのいずれか一つに記載の四輪独立駆動車の制駆動力配分装置。 A braking / driving force lower limit estimating means for estimating a braking / driving force lower limit that can be output by an actuator that drives or brakes each wheel;
5. The braking / driving force distribution set calculated by the braking / driving force distribution set calculating means is limited to a braking / driving force lower limit that can be output by an actuator that drives or brakes each wheel. The braking / driving force distribution device for a four-wheel independent drive vehicle according to any one of the above.
前記制駆動力配分選択手段は、各輪のエネルギー損失の和が最小となる制駆動力配分を選択することを特徴とする請求項1から8までのいずれか一つに記載の四輪独立駆動車の制駆動力配分装置。 Energy loss estimating means for estimating each energy loss generated when each wheel is driven,
The four-wheel independent drive according to any one of claims 1 to 8, wherein the braking / driving force distribution selection unit selects a braking / driving force distribution that minimizes a sum of energy losses of the wheels. Vehicle braking / driving force distribution device.
各輪の再配分後の制駆動力におけるタイヤ横力を推定するタイヤ横力推定手段と
を備え、
前記制駆動力配分選択手段は、タイヤ横力と制駆動力との合力であるタイヤ力を輪荷重で除したタイヤ摩擦円使用率を各輪毎に求め、この求めた各輪1〜4のタイヤ摩擦円使用率の最大値と最小値の差が最も小さくなる制駆動力配分を選択することを特徴とする請求項1から8までのいずれか一つに記載の四輪独立駆動車の制駆動力配分装置。 Wheel load estimating means for estimating the wheel load of each wheel;
Tire lateral force estimating means for estimating tire lateral force in braking / driving force after redistribution of each wheel,
The braking / driving force distribution selection means obtains a tire friction circle usage rate for each wheel obtained by dividing the tire force, which is the resultant force of the tire lateral force and braking / driving force, by the wheel load, for each of the obtained wheels 1 to 4. 9. The braking / driving force distribution according to any one of claims 1 to 8, wherein a braking / driving force distribution that minimizes a difference between a maximum value and a minimum value of a tire friction circle usage rate is selected. Driving force distribution device.
前記制駆動力配分選択手段は、前記総エネルギー必要量が車両の出力可能なエネルギーの範囲内にある制駆動力配分を選択することを特徴とする請求項1から8までのいずれか一つに記載の四輪独立駆動車の制駆動力配分装置。 A total energy requirement estimating means for estimating a total energy requirement for realizing each braking / driving force distribution among the set of braking / driving force distributions of the four wheels;
9. The braking / driving force distribution selecting unit selects a braking / driving force distribution in which the total energy requirement is within a range of energy that can be output by the vehicle. The braking / driving force distribution device for the described four-wheel independent drive vehicle.
前記総エネルギー必要量推定手段は、前記総エネルギー必要量をこの各輪の回転速度に基づいて推定することを特徴とする請求項14に記載の四輪独立駆動車の制駆動力配分装置。 Each wheel rotation speed detection means for detecting the rotation speed of each wheel,
15. The braking / driving force distribution device for a four-wheel independent drive vehicle according to claim 14, wherein the total energy requirement estimation means estimates the total energy requirement based on the rotational speed of each wheel.
前記総エネルギー必要量推定手段は、前記総エネルギー必要量をこの各エネルギー損失に基づいて推定することを特徴とする請求項14に記載の四輪独立駆動車の制駆動力配分装置。 Energy loss estimating means for estimating each energy loss generated when each wheel is driven,
15. The braking / driving force distribution device for a four-wheel independent drive vehicle according to claim 14, wherein the total energy requirement estimation means estimates the total energy requirement based on each energy loss.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005180371A JP4956772B2 (en) | 2005-06-21 | 2005-06-21 | Braking / driving force distribution device for four-wheel independent drive vehicle |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005180371A JP4956772B2 (en) | 2005-06-21 | 2005-06-21 | Braking / driving force distribution device for four-wheel independent drive vehicle |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2007001330A true JP2007001330A (en) | 2007-01-11 |
JP4956772B2 JP4956772B2 (en) | 2012-06-20 |
Family
ID=37687326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005180371A Expired - Lifetime JP4956772B2 (en) | 2005-06-21 | 2005-06-21 | Braking / driving force distribution device for four-wheel independent drive vehicle |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4956772B2 (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008074191A (en) * | 2006-09-20 | 2008-04-03 | Nissan Motor Co Ltd | Driving force distribution controller for vehicle |
JP2008239115A (en) * | 2007-03-29 | 2008-10-09 | Honda Motor Co Ltd | Vehicle motion control device |
JP2009035047A (en) * | 2007-07-31 | 2009-02-19 | Nissan Motor Co Ltd | Turning behavior control method and device for vehicle |
KR100911554B1 (en) | 2008-04-23 | 2009-08-10 | 현대자동차주식회사 | Brake piping device |
JP2009208492A (en) * | 2008-02-29 | 2009-09-17 | Nissan Motor Co Ltd | Travelling posture control device and travelling posture control method of vehicle |
JP2009214742A (en) * | 2008-03-11 | 2009-09-24 | Nissan Motor Co Ltd | Lane deviation preventive device and its method |
JP2010120600A (en) * | 2008-11-21 | 2010-06-03 | Fuji Heavy Ind Ltd | Estimation device of friction coefficient of road surface |
JP2010149740A (en) * | 2008-12-25 | 2010-07-08 | Fuji Heavy Ind Ltd | Vehicle motion control apparatus |
US8195361B2 (en) | 2007-06-27 | 2012-06-05 | Toyota Jidosha Kabushiki Kaisha | Vehicle suspension system |
KR20140133276A (en) * | 2013-05-10 | 2014-11-19 | 현대위아 주식회사 | Control method of 4 wheel drive vehicle |
JP2015023691A (en) * | 2013-07-19 | 2015-02-02 | 富士重工業株式会社 | Vehicle drive control device |
CN115139815A (en) * | 2022-06-27 | 2022-10-04 | 重庆金康赛力斯新能源汽车设计院有限公司 | Torque distribution method, device, equipment and storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005073457A (en) * | 2003-08-27 | 2005-03-17 | Toyota Motor Corp | Vehicle control device |
JP2005112007A (en) * | 2003-10-02 | 2005-04-28 | Toyoda Mach Works Ltd | Vehicular integrated control device |
JP2005119647A (en) * | 2003-09-25 | 2005-05-12 | Nissan Motor Co Ltd | Driving force control device for wheel independent drive system vehicle |
-
2005
- 2005-06-21 JP JP2005180371A patent/JP4956772B2/en not_active Expired - Lifetime
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005073457A (en) * | 2003-08-27 | 2005-03-17 | Toyota Motor Corp | Vehicle control device |
JP2005119647A (en) * | 2003-09-25 | 2005-05-12 | Nissan Motor Co Ltd | Driving force control device for wheel independent drive system vehicle |
JP2005112007A (en) * | 2003-10-02 | 2005-04-28 | Toyoda Mach Works Ltd | Vehicular integrated control device |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008074191A (en) * | 2006-09-20 | 2008-04-03 | Nissan Motor Co Ltd | Driving force distribution controller for vehicle |
US9132731B2 (en) | 2006-09-20 | 2015-09-15 | Nissan Motor Co., Ltd. | Vehicle drive force distribution control system |
JP2008239115A (en) * | 2007-03-29 | 2008-10-09 | Honda Motor Co Ltd | Vehicle motion control device |
US8195361B2 (en) | 2007-06-27 | 2012-06-05 | Toyota Jidosha Kabushiki Kaisha | Vehicle suspension system |
JP2009035047A (en) * | 2007-07-31 | 2009-02-19 | Nissan Motor Co Ltd | Turning behavior control method and device for vehicle |
JP2009208492A (en) * | 2008-02-29 | 2009-09-17 | Nissan Motor Co Ltd | Travelling posture control device and travelling posture control method of vehicle |
JP2009214742A (en) * | 2008-03-11 | 2009-09-24 | Nissan Motor Co Ltd | Lane deviation preventive device and its method |
KR100911554B1 (en) | 2008-04-23 | 2009-08-10 | 현대자동차주식회사 | Brake piping device |
JP2010120600A (en) * | 2008-11-21 | 2010-06-03 | Fuji Heavy Ind Ltd | Estimation device of friction coefficient of road surface |
JP2010149740A (en) * | 2008-12-25 | 2010-07-08 | Fuji Heavy Ind Ltd | Vehicle motion control apparatus |
KR20140133276A (en) * | 2013-05-10 | 2014-11-19 | 현대위아 주식회사 | Control method of 4 wheel drive vehicle |
KR102052683B1 (en) | 2013-05-10 | 2019-12-05 | 현대위아 주식회사 | Control method of 4 wheel drive vehicle |
JP2015023691A (en) * | 2013-07-19 | 2015-02-02 | 富士重工業株式会社 | Vehicle drive control device |
CN115139815A (en) * | 2022-06-27 | 2022-10-04 | 重庆金康赛力斯新能源汽车设计院有限公司 | Torque distribution method, device, equipment and storage medium |
CN115139815B (en) * | 2022-06-27 | 2024-04-09 | 重庆金康赛力斯新能源汽车设计院有限公司 | Torque distribution method, device, equipment and storage medium |
Also Published As
Publication number | Publication date |
---|---|
JP4956772B2 (en) | 2012-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7216943B2 (en) | Hybrid vehicle | |
US8924055B2 (en) | Vehicle control apparatus | |
JP4929944B2 (en) | Vehicle driving force distribution control device | |
US20040186647A1 (en) | Vehicle control method and vehicle control apparatus | |
CN105848946B (en) | Vehicle and vehicle control method | |
JP4956772B2 (en) | Braking / driving force distribution device for four-wheel independent drive vehicle | |
JP6379148B2 (en) | Vehicle driving force control device | |
US10518775B1 (en) | Regenerative braking of vehicle with primary and secondary drive axles | |
CN107249967A (en) | Vehicle | |
JP5055992B2 (en) | Hybrid vehicle driving force control device | |
JP4779943B2 (en) | Hybrid vehicle | |
JP2008141875A (en) | Traveling device and drive control device | |
JP4428162B2 (en) | Brake control device | |
JP4892983B2 (en) | Driving force distribution device for four-wheel independent drive vehicle | |
JP4961751B2 (en) | Vehicle driving force distribution device | |
JP7176360B2 (en) | electric vehicle | |
JP2008024204A (en) | Control device of hybrid vehicle | |
JP2010151205A (en) | Driving force distribution control device for vehicle | |
JP5163500B2 (en) | Hybrid vehicle driving force distribution control device | |
JP2012171616A (en) | Torque distributor and method of distributing torque | |
JP5104102B2 (en) | Vehicle driving force distribution control device | |
KR20150018010A (en) | Apparatus and method for controlling posture of hybrid electric vehicle using regenerative breaking | |
JP4918787B2 (en) | Driving force distribution device for four-wheel independent drive vehicle | |
JP5096623B2 (en) | Torque distribution device and torque distribution method | |
US10173664B2 (en) | Expanding regenerative capacity up to vehicle dynamic limits through integration with mitigative subsystems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080422 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20110117 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110621 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110802 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20120214 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120227 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150330 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4956772 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |