JP2000052951A - Method and apparatus for estimating vehicle sideslip angle of vehicle - Google Patents
Method and apparatus for estimating vehicle sideslip angle of vehicleInfo
- Publication number
- JP2000052951A JP2000052951A JP10227242A JP22724298A JP2000052951A JP 2000052951 A JP2000052951 A JP 2000052951A JP 10227242 A JP10227242 A JP 10227242A JP 22724298 A JP22724298 A JP 22724298A JP 2000052951 A JP2000052951 A JP 2000052951A
- Authority
- JP
- Japan
- Prior art keywords
- speed
- vehicle
- wheel
- estimated
- estimating
- 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
Landscapes
- Steering Control In Accordance With Driving Conditions (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
- Regulating Braking Force (AREA)
Abstract
(57)【要約】
【課題】車体横滑り角を推定する際に、低車速域での推
定精度を向上させると共に、ドリフトを生じることを抑
制することができる横滑り角推定方法及び装置を提供す
る。
【解決手段】車体横滑り角を推定する際に、例えば車輪
速度に基づいて前後方向速度VX を算出すると共に、ヨ
ーレートγ、前後方向加速度αX 及び横方向加速度αY
を検出し、前後方向速度VX とその推定値である推定前
後方向速度V^Xとの誤差信号をフィードバックするこ
とにより、未知量である推定横方向速度V^Y を算出
し、推定前後方向速度V^X と推定横方向速度V^Y と
から推定横滑り角β^を算出する。
(57) Abstract: A method and an apparatus for estimating a side slip angle that can improve the estimation accuracy in a low vehicle speed range and suppress the occurrence of drift when estimating a vehicle body side slip angle. When estimating A vehicle body slip angle, for example, to calculate the longitudinal direction velocity V X based on the wheel speed, yaw rate gamma, longitudinal acceleration alpha X and lateral acceleration alpha Y
Detects, by feeding back the error signal between the front-rear direction velocity V X and the estimated longitudinal direction velocity V ^ X which is the estimated value, to calculate the estimated lateral velocity V ^ Y is unknown amount, estimated longitudinal direction An estimated sideslip angle β ^ is calculated from the speed V ^ X and the estimated lateral speed V ^ Y.
Description
【0001】[0001]
【発明の属する技術分野】この発明は、車両の車体横滑
り角推定方法及び推定装置に関し、特に、低車速域での
車体横滑り角の推定精度の向上を図ったものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method and an apparatus for estimating a vehicle body sideslip angle of a vehicle, and more particularly, to an improvement in estimation accuracy of a vehicle body sideslip angle in a low vehicle speed range.
【0002】[0002]
【従来の技術】この種の従来の技術としては、例えば特
開平8−268306号公報(以下、第1従来例と称
す)及び特開平8−332934号公報(以下、第2従
来例と称す)に開示されたものが知られている。2. Description of the Related Art Japanese Patent Application Laid-Open No. 8-268306 (hereinafter referred to as a first conventional example) and Japanese Patent Application Laid-Open No. 8-332934 (hereinafter referred to as a second conventional example) are examples of such prior art. Is disclosed.
【0003】第1従来例は、オブザーバーの状態方程式
に基づいて計算した推定横滑り角及び推定ヨーレートと
検出前後速度とにより推定横滑り角を計算する一方で、
検出横加速度と前記計算した推定横加速度との比を補正
係数として計算し、その補正係数をオブザーバーの状態
方程式におけるコーナリングパワーに関する因子に乗算
するようになっていて、これにより、オブザーバーの状
態方程式に基づいて横滑り角を推定すると、コーナリン
グパワーが的確に補正されているので、横滑り角の推定
精度が良好になる、というものである。The first conventional example calculates an estimated sideslip angle based on an estimated sideslip angle and an estimated yaw rate calculated based on an observer's state equation and a detected front-rear speed,
The ratio between the detected lateral acceleration and the calculated estimated lateral acceleration is calculated as a correction coefficient, and the correction coefficient is multiplied by a factor relating to cornering power in the observer's state equation. When the sideslip angle is estimated based on the cornering power, since the cornering power is accurately corrected, the estimation accuracy of the sideslip angle is improved.
【0004】また、第2従来例は、車速、ヨーレート、
及び横加速度をパラメータとする積分演算により車両の
横滑り物理量を演算する横滑り物理量演算ブロックと、
検出された車両の状態量を入力変数として車両モデルに
より車両の横滑り量を推定する横滑り物理量推定ブロッ
クと、車両の旋回挙動が不安定なときには積分演算の積
分時定数を大きくする積分時定数設定ブロックとを有
し、横滑り物理量演算ブロックにより演算された横滑り
物理量及び横滑り物理量推定ブロックにより推定された
横滑り物理量に基づき最終演算ブロックにより最終的に
車両の横滑り物理量を演算するようにしている。[0004] The second conventional example includes a vehicle speed, a yaw rate,
And a skid physical quantity calculation block for calculating a skid physical quantity of the vehicle by an integration calculation using the lateral acceleration as a parameter,
A side slip physical quantity estimation block for estimating a side slip amount of a vehicle by a vehicle model using the detected state quantity of the vehicle as an input variable, and an integration time constant setting block for increasing an integration time constant of an integration operation when the turning behavior of the vehicle is unstable. The final calculation block finally calculates the skid physical quantity of the vehicle based on the skid physical quantity calculated by the skid physical quantity calculation block and the skid physical quantity estimated by the skid physical quantity estimation block.
【0005】[0005]
【発明が解決しようとする課題】ここで、一般的にタイ
ヤの横力やコーナリングフォースは、タイヤスリップ角
が大きな領域では非線形な特性を有するため、タイヤモ
デルに基づいた推定式を用いて車両の車体横滑り角を推
定する場合に、タイヤ特性(タイヤ横滑り角とコーナリ
ングフォースとの関係)が非線形な領域では車体横滑り
角の推定精度が低下することが判っている。Here, since the lateral force and the cornering force of the tire generally have nonlinear characteristics in a region where the tire slip angle is large, an estimation formula based on a tire model is used for the vehicle. When estimating the vehicle body side slip angle, it has been found that the estimation accuracy of the vehicle body side slip angle decreases in a region where the tire characteristics (the relationship between the tire side slip angle and the cornering force) are nonlinear.
【0006】そのため、例えば上記第1従来例のよう
に、検出横加速度と推定横加速度との比に応じた補正係
数を用いる等して推定式に含まれるコーナリングパワー
を補正することにより、車体横滑り角の推定精度が低下
しないようにすることが考えられる。しかし、このよう
な補正係数を用いる方法であっても、車体質量、車体ヨ
ーイング慣性モーメント、車両重心から前後車軸までの
水平距離、車両前後速度、コーナリングパワー等を使用
した車両の旋回運動モデルを使用するため、車速が低い
領域では近似式が有効でないために、横滑り角の推定精
度が悪化するという未解決の課題がある。特に、後輪を
駆動する自動車の発進直後の駆動力により所謂パワーオ
ーバーステアと呼ばれる現象をはじめ、低速走行時にも
車体横滑り角が増大する可能性は高くこのような領域で
横滑り角を高精度で推定することが望まれている。Therefore, as in the first prior art, for example, the cornering power included in the estimation formula is corrected by using a correction coefficient corresponding to the ratio between the detected lateral acceleration and the estimated lateral acceleration, so that the vehicle body skids. It is conceivable that the estimation accuracy of the angle is not reduced. However, even with the method using such a correction coefficient, a turning motion model of the vehicle using the vehicle mass, the vehicle yawing inertia moment, the horizontal distance from the vehicle center of gravity to the front and rear axles, the vehicle front and rear speed, the cornering power, etc. is used. Therefore, there is an unsolved problem that the approximation formula is not effective in a region where the vehicle speed is low, so that the accuracy of estimating the sideslip angle deteriorates. In particular, there is a high possibility that the vehicle body slip angle increases even at low speeds, including a phenomenon called so-called power oversteer due to the driving force immediately after the vehicle that drives the rear wheels is started. It is desired to estimate.
【0007】また、第2従来例のように、車速、ヨーレ
ート、及び横加速度をパラメータとする積分演算により
車両の横滑り物理量を演算する場合には、上記第1従来
例のような低車速域での推定精度の悪化という問題点は
発生しないものの、積分処理を使用しているためにセン
サ間の動特性の違い、積分誤差によって推定値が真値か
ら外れていくドリフトが発生する可能性があり、このた
めに車両モデルを使用する方法との併用等が必要となる
という未解決の課題がある。Further, when the vehicle side slip physical quantity is calculated by the integral calculation using the vehicle speed, the yaw rate, and the lateral acceleration as parameters, as in the second conventional example, the vehicle is operated in a low vehicle speed range as in the first conventional example. Although the problem of deterioration of the estimation accuracy does not occur, there is a possibility that the estimation value deviates from the true value due to differences in dynamic characteristics between sensors and integration errors due to the use of integration processing. However, there is an unsolved problem that it is necessary to use the method together with a method using a vehicle model.
【0008】そこで、本発明は、上記従来例の未解決の
課題に着目してなされたものであり、低車速域での車体
横滑り角の推定精度を向上させることができる横滑り角
推定方法及び推定装置を提供することを目的としてい
る。Therefore, the present invention has been made in view of the above-mentioned unsolved problems of the prior art, and a method and an apparatus for estimating a sideslip angle capable of improving the accuracy of estimating the sideslip angle of a vehicle body in a low vehicle speed range. It is intended to provide a device.
【0009】[0009]
【課題を解決するための手段】上記目的を達成するため
に、請求項1に係る車両の車体横滑り角推定方法は、車
両の前後方向及び横方向について、夫々当該車両が鉛直
軸周りの運動をしている状態を記述する運動モデルを有
し、該運動モデルは、前後方向の運動を記述する運動モ
デルから推定される推定前後運動物理量と、実際に測定
される前後運動物理量及び他の推定方法から得られる前
後運動物理量の何れかとの誤差信号をフィードバック信
号として使用し、前記運動モデルに基づいて車体の横滑
り角を推定するようにしたことを特徴としている。According to a first aspect of the present invention, there is provided a method for estimating a vehicle body skid angle of a vehicle according to the first aspect of the present invention, wherein the vehicle moves about a vertical axis in a front-rear direction and a lateral direction, respectively. A motion model that describes a state in which the motion is performed. The motion model includes an estimated front-rear motion physical quantity estimated from the motion model that describes the front-rear motion, and an actually measured front-rear motion physical quantity and other estimation methods. An error signal with any of the forward and backward movement physical quantities obtained from is used as a feedback signal, and the side slip angle of the vehicle body is estimated based on the movement model.
【0010】また、請求項2に係る車両の車体横滑り角
推定方法は、請求項1に係る発明において、前記前後方
向物理量として、前後方向速度に設定されていることを
特徴としている。According to a second aspect of the present invention, there is provided the vehicle body side slip angle estimating method according to the first aspect, wherein the longitudinal physical quantity is set to a longitudinal velocity.
【0011】さらに、請求項3に係る車両の車体横滑り
角推定方法は、請求項2に係る発明において、前記前後
方向速度が、車輪速度及び車輪速度の推定値の何れかを
使用することを特徴としている。Further, in the vehicle body side slip angle estimation method according to a third aspect, in the invention according to the second aspect, the front-rear direction speed uses one of a wheel speed and an estimated value of the wheel speed. And
【0012】さらにまた、請求項4に係る車両の車体横
滑り角推定方法は、駆動輪の車輪速の大きな方の値が、
非駆動輪のそれに対してある一定値又は一定比率以上大
きく、且つ前輪の車輪速又は車体速度がある一定値より
小さい場合に請求項1の車体横滑り角推定方法を使用
し、それ以外の場合に他の車体横滑り角推定方法を使用
するようにしたことを特徴としている。Further, in the method for estimating a vehicle body skid angle according to claim 4, the larger value of the wheel speed of the drive wheel is determined by:
The vehicle body slip angle estimation method according to claim 1 is used when the wheel speed or the vehicle body speed of the front wheels is smaller than a certain value by more than a certain value or a certain ratio with respect to that of the non-driving wheels. It is characterized in that another method of estimating a vehicle body skid angle is used.
【0013】なおさらに、請求項5に係る車両の車体横
滑り角推定装置は、車両の各車輪の車輪速度を算出する
車輪速度算出手段と、車体に発生するヨーレートを検出
するヨーレート検出手段と、前記車輪速度算出手段で検
出した車輪速度と前記ヨーレート検出手段で検出したヨ
ーレートとに基づいて前後方向速度を算出する前後方向
速度算出手段と、該前後方向速度算出手段で算出した前
後方向速度と推定前後方向速度との誤差信号をフィード
バックして前後方向及び横方向速度を推定する速度推定
手段と、該速度推定手段で推定した前後方向及び横方向
速度に基づいて横滑り角を算出する横滑り角算出手段と
を備えたことを特徴としている。Still further, according to a fifth aspect of the present invention, there is provided a vehicle body slip angle estimating apparatus for calculating a wheel speed of each wheel of a vehicle, a yaw rate detecting means for detecting a yaw rate generated in the vehicle body, Forward / backward speed calculating means for calculating a forward / backward speed based on the wheel speed detected by the wheel speed calculating means and the yaw rate detected by the yaw rate detecting means; forward / backward speed calculated by the forward / backward speed calculating means; Speed estimating means for estimating the longitudinal and lateral speeds by feeding back an error signal from the direction velocity, and a skid angle calculating means for calculating a skid angle based on the longitudinal and lateral speeds estimated by the speed estimating means; It is characterized by having.
【0014】[0014]
【発明の効果】請求項1に係る発明によれば、前後方向
の運動を記述する運動モデルから推定される推定前後運
動物理量と、実際に測定される前後運動物理量又は他の
推定方法から得られる前後運動物量との誤差信号をフィ
ードバック信号として使用するので、低車速域で高精度
の横滑り角の推定を行うことができると共に、推定値が
初期値を誤って与えた等の何らかの理由で真値を離れた
場合にも真値に収束させることができ、前述した第2従
来例のようにドリフトが生じることがないという効果が
得られる。According to the first aspect of the present invention, an estimated front-rear movement physical quantity estimated from a movement model describing a front-rear movement and obtained from an actually measured front-rear movement physical quantity or another estimation method. Since the error signal with the forward and backward moving object quantity is used as the feedback signal, it is possible to estimate the sideslip angle with high accuracy in the low vehicle speed range, and the estimated value is the true value for some reason such as incorrectly giving the initial value. Can be converged to the true value even when the distance is deviated, and the effect that drift does not occur unlike the second conventional example described above is obtained.
【0015】また、請求項2に係る発明によれば、前後
方向運動物理量として、前後方向速度が設定されている
ので、正確な実測が可能であり、より高精度の横滑り角
推定を行うことができるという効果が得られる。According to the second aspect of the present invention, since the longitudinal velocity is set as the longitudinal physical quantity, accurate actual measurement is possible, and a more accurate side slip angle estimation can be performed. The effect that it can be obtained is obtained.
【0016】さらに、請求項3に係る発明によれば、前
後方向速度を、車輪速度及び車輪速度の推定値の何れか
を使用して算出するので、実際の車両の走行状態に応じ
た値として算出することができ、より高精度の横滑り角
推定を行うことができるという効果が得られる。Further, according to the third aspect of the present invention, since the longitudinal speed is calculated using either the wheel speed or the estimated value of the wheel speed, the speed is calculated as a value corresponding to the actual running state of the vehicle. It is possible to calculate, and it is possible to obtain an effect that a more accurate side slip angle estimation can be performed.
【0017】さらにまた、請求項4に係る発明によれ
ば、発進時等の低車速域では、請求項1の発明を適用し
て、高精度の横滑り角推定を行い、中・高車速域では、
他の推定方法に切換えることにより、高精度の横滑り角
推定を継続し、すべての車速域で横滑り角を高精度に推
定することができるという効果が得られる。Furthermore, according to the invention of claim 4, in a low vehicle speed range at the time of starting or the like, the invention of claim 1 is applied to perform a highly accurate side slip angle estimation, and in a middle / high vehicle speed range. ,
By switching to another estimation method, it is possible to obtain an effect that high-accuracy sideslip angle estimation can be continued and the sideslip angle can be estimated with high accuracy in all vehicle speed ranges.
【0018】なおさらに、請求項5に係る発明によれ
ば、請求項1に係る発明と同様に、低車速域で高精度の
横滑り角の推定を行うことができると共に、推定値が初
期値を誤って与えた等の何らかの理由で真値を離れた場
合にも真値に収束させることができ、前述した第2従来
例のようにドリフトが生じることがないという効果が得
られる。Further, according to the invention of claim 5, similarly to the invention of claim 1, it is possible to perform a highly accurate estimation of the sideslip angle in a low vehicle speed range, and the estimated value is set to an initial value. The true value can be converged to the true value even if the true value is deviated for some reason, such as being given by mistake, and the effect that drift does not occur unlike the second conventional example described above is obtained.
【0019】[0019]
【発明の実施の形態】以下、本発明の実施の形態を図面
に基づいて説明する。図1は本発明の第1の実施の形態
の全体構成を示すブロック図であり、本実施の形態にお
ける車両は、車体横滑り角等の状態量を推定する車両状
態量推定装置10と、目標制動力設定装置20と、制動
力制御装置30と、を備えて構成されている。なお、目
標制動力設定装置20は、車体横滑り角推定装置10か
ら供給される各値等に基づいて適切な制動力を設定する
ための装置であって、その具体的構成は、実行される制
御(例えば、VDC,ABS,TCS等のように車両に
適用可能な制御)に応じて適宜決定されるものである
が、その内容は本発明の本質ではないため、その具体的
な説明は省略する。また、制動力制御装置30も、実行
される制御に応じて適宜公知の構成が採用可能であるか
ら、その具体的な説明も省略する。Embodiments of the present invention will be described below with reference to the drawings. FIG. 1 is a block diagram showing the overall configuration of a first embodiment of the present invention. A vehicle in this embodiment includes a vehicle state quantity estimating device 10 for estimating a state quantity such as a vehicle body side slip angle, and a target control system. A power setting device 20 and a braking force control device 30 are provided. The target braking force setting device 20 is a device for setting an appropriate braking force based on each value or the like supplied from the vehicle body side slip angle estimating device 10, and its specific configuration is the control to be executed. (For example, control applicable to a vehicle such as VDC, ABS, TCS, etc.), but the content is not the essence of the present invention, and a specific description thereof is omitted. . In addition, the braking force control device 30 can appropriately adopt a known configuration according to the control to be performed, and thus a specific description thereof is omitted.
【0020】そして、車両状態量推定装置10は、実際
にはマイクロコンピュータや必要なインタフェース回路
等によって構成され、車両に発生しているヨーレートを
検出するヨーレートセンサ11から供給される実ヨーレ
ートrと、車両に発生している前後加速度を検出する前
後加速度センサ12から供給される前後加速度αX と、
車両に発生している横加速度を検出する横加速度センサ
13から供給される横加速度αY と、各車輪の車輪速を
検出する車輪速センサ14から供給される車輪速VwFL
〜VwRRとに基づいて、所定の演算処理を実行して、車
体の横滑り角βを推定する。The vehicle state quantity estimating apparatus 10 is actually constituted by a microcomputer, necessary interface circuits and the like, and is provided with an actual yaw rate r supplied from a yaw rate sensor 11 for detecting a yaw rate generated in the vehicle. a longitudinal acceleration alpha X supplied from the longitudinal acceleration sensor 12 for detecting a longitudinal acceleration generated in the vehicle,
Lateral acceleration α Y supplied from a lateral acceleration sensor 13 detecting the lateral acceleration occurring in the vehicle, and wheel speed Vw FL supplied from a wheel speed sensor 14 detecting the wheel speed of each wheel.
A predetermined calculation process is executed based on ~ Vw RR to estimate the sideslip angle β of the vehicle body.
【0021】すなわち、車両の前後方向及び横方向につ
いて車両が鉛直軸周りの運動をしている状態を記述する
一般的な2輪運動モデルは、図2に示すように、操舵機
構を有する前輪1F及び操舵機構を有さない後輪1Rを
有し、前輪操舵角をδ、前輪横滑り角をβf 、後輪横滑
り角をβr 、重心点2での車速をV、前後方向速度をV
X 、横方向速度をVY 、前後加速度をαX 、横加速度を
αY 、ヨーレートをγ、横滑り角をβとし、重心点から
前輪側車軸までの水平距離をLf 、重心点から後輪側車
軸までの水平距離をLr としたときに、その運動方程式
は、車両の質量をM、車両のヨー方向回転慣性モーメン
トをIZ 、前後輪のタイヤコーナリングパワーをC
Pf ,CPr とすると、下記(1)式及び(2)式で表
すことができる。That is, as shown in FIG. 2, a general two-wheel motion model that describes a state in which the vehicle is moving around a vertical axis in the front-rear direction and the lateral direction of the vehicle is a front wheel 1F having a steering mechanism. And the rear wheel 1R without a steering mechanism, the front wheel steering angle is δ, the front wheel side slip angle is β f , the rear wheel side slip angle is β r , the vehicle speed at the center of gravity 2 is V, and the front-rear direction speed is V
X , lateral speed V Y , longitudinal acceleration α X , lateral acceleration α Y , yaw rate γ, sideslip angle β, horizontal distance from the center of gravity to the front wheel side axle L f , center of gravity to rear wheel Assuming that the horizontal distance to the side axle is L r , the equation of motion is that the mass of the vehicle is M, the rotational moment of inertia of the vehicle in the yaw direction is I Z , and the tire cornering power of the front and rear wheels is C.
If P f and CP r are used, they can be expressed by the following equations (1) and (2).
【0022】 MVX (β′+γ)=2CPf βf +2CPr βr …………(1) IZ γ′=2Lf CPf βf −2Lf CPf βf …………(2) そして、横滑り角が大きな状態の運動を正確に記述する
ために、図2の運動モデルから車両の前後加速度αX 横
加速度αY と前後方向速度VX 及び横方向速度VY との
関係は、下記(3)式及び(4)式で表すことができ
る。[0022] MV X (β '+ γ) = 2CP f β f + 2CP r β r ............ (1) I Z γ' = 2L f CP f β f -2L f CP f β f ............ (2 ) in order to sideslip angle accurately describe the motion of large states, the relationship between the longitudinal acceleration alpha X lateral acceleration alpha Y and the longitudinal direction velocity V X and lateral velocity V Y of the vehicle from the motion model of Figure 2 , (3) and (4) below.
【0023】 αX =VX ′+VY γ …………(3) αY =VY ′−Vx γ …………(4) ここで、VX ′は前後方向速度VX の微分値、VY ′は
横方向速度の微分値である。Α X = V X ′ + V Y γ (3) α Y = V Y ′ −V x γ (4) where V X ′ is the differential of the longitudinal velocity V X. The value, V Y ', is the derivative of the lateral velocity.
【0024】上記(3)は、車両が減速を行っている場
合に、もし車両が重心点での横方向の速度成分とヨーレ
イトとが発生していれば、前後加速度センサの検出値
と、前後方向加速度の微分値が一致しないことを示して
いる。The above (3) means that if the vehicle is decelerating, if the vehicle has a lateral velocity component and a yaw rate at the center of gravity, the detected value of the longitudinal acceleration sensor and the longitudinal This indicates that the differential values of the directional acceleration do not match.
【0025】また、上記(4)式は、車両に遠心力が働
いている場合には、車両の横方向速度VY の微分値
VY ′と横加速度センサの検出値とが一致しないことを
示している。The above equation (4) indicates that when a centrifugal force acts on the vehicle, the differential value V Y 'of the lateral speed V Y of the vehicle does not match the detection value of the lateral acceleration sensor. Is shown.
【0026】これら(3)式及び(4)式から分かるよ
うに、車両にヨーレートが発生している状態では、前後
方向及び横方向の車両の運動はヨーレートを介して密接
に関係している。As can be seen from the equations (3) and (4), when the vehicle has a yaw rate, the vehicle movement in the front-rear direction and the lateral direction is closely related via the yaw rate.
【0027】また、前後方向速度VX が横方向速度VY
に対して十分に大きくない状況においては、それを考え
て車両の運動をモデル化する必要がある。そこで、上記
(3)式及び(4)式に使われる変数のうち車両におい
てセンサによって直接測定が可能な物理量は、ヨーレー
トγ、横前後加速度α及び横方向加速度αY である。[0027] In addition, the front-rear direction speed V X is the lateral velocity V Y
In situations where the vehicle motion is not large enough, it is necessary to model the motion of the vehicle with this in mind. Therefore, the expressions (3) and (4) a physical quantity which can be measured directly by sensors in the vehicle of the variables used in the formula, yaw rate gamma, a lateral longitudinal acceleration alpha and lateral acceleration alpha Y.
【0028】また、車両の前後方向の速度VX は、車輪
の速度等によって比較的容易に測定若しくは推定するこ
とが可能である。さらに、上記(3)式及び(4)式
は、前後方向の加速度αX 及び横方向の加速度αy と前
後方向速度VX 及び横方向速度VY との関係を表す式で
あるが、同時に前後方向速度VX 及び横方向速度VY を
状態量とした下記(5)及び(6)式で表される運動方
程式と考えることが可能となる。Further, the speed V X in the longitudinal direction of the vehicle, can be relatively easily measured or estimated by the speed of the wheel. Further, the above equations (3) and (4) are equations expressing the relationship between the longitudinal acceleration α X and the lateral acceleration α y and the longitudinal velocity V X and the lateral velocity V Y. it is possible to consider longitudinal direction velocity V X and lateral velocity V Y quantity of state was the following (5) and (6) and equation of motion of the formula.
【0029】 VX ′=−VY γ+αX …………(5) VY ′=VX γ+αY …………(6) したがって、上記(5)式及び(6)式を直接使用して
車両の運動を記述することが可能となる。V X ′ = −V Y γ + α X (5) V Y ′ = V X γ + α Y (6) Therefore, the above equations (5) and (6) are directly used. Thus, the motion of the vehicle can be described.
【0030】このため、上記(5)式及び(6)式をも
とに前後方向速度VX を既知量としてのその誤差信号を
用い、未知量である横方向速VY を推定するオブザーバ
を設定すれば、横滑り角が大きな領域においても横滑り
角を高精度で求めることが可能となる。[0030] Thus, using the error signal as a known quantity in the longitudinal direction velocity V X based on the above equation (5) and (6), an observer for estimating the lateral velocity V Y is unknown amount If set, the sideslip angle can be determined with high accuracy even in a region where the sideslip angle is large.
【0031】ここで、前記(5)式及び(6)式は非線
型な運動式となっているために下記(7)式及び(8)
式に示すオブザーバを使用する。 V^X ′=-V^Y γ+αX −K11(V^X −VX )-K12sign(V^X −VX ) …………(7) V^Y ′= V^X γ+αY −K21(V^X −VX )-K22sign(V^X −VX ) …………(8) ここで、K11,K21は誤差信号に対して線型なフィード
バックを与えるゲイン、K12,K22は誤差の符号につい
てのみフィードバックを与えるゲインとなっている。こ
れらのゲインの値は実験結果などから決定することがで
きる。Here, since the above equations (5) and (6) are nonlinear motion equations, the following equations (7) and (8)
Use the observer shown in the equation. V ^ X '= -V ^ Y γ + α X -K 11 (V ^ X -V X) -K 12 sign (V ^ X -V X) ............ (7) V ^ Y' = V ^ X γ + α Y -K 21 (V ^ X -V X) -K 22 sign (V ^ X -V X) ............ (8) where, K 11, K 21 gives linear feedback with respect to the error signal The gains K 12 and K 22 are gains that provide feedback only for the sign of the error. These gain values can be determined from experimental results and the like.
【0032】そして、上記(7)式及び(8)式で推定
された前後方向速度V^X 及び横方向速度V^Y を用い
て、下記(9)式に従って最終的な横滑り角β^を算出
する。Using the longitudinal velocity V ^ X and the lateral velocity V ^ Y estimated by the above equations (7) and (8), the final side slip angle β 角 is calculated according to the following equation (9). calculate.
【0033】 β^=tan-1(V^Y /V^X ) …………(9) このように、上記第1の実施形態においては、車両状態
量推定装置10で、ヨーレイトセンサ11で検出した実
ヨーレートγ、前後加速度センサ12で検出した前後加
速度αX ,横加速度センサ13で検出した横加速度
αY 、車輪速センサ14で検出した車輪速VwFL〜Vw
RRに基づいて算出される前後方向速度VX に基づいて前
記(7)式および(8)式のオブザーバ演算を行うこと
により、横方向速度VY に推定前後方向速度V^X と実
際に検出した前後方向速度VX との誤差信号をフィード
バックして、推定横方向速度V^Y を算出し、これと推
定前後方向速度V^X とに基づいて推定横滑り角β^を
算出することができる。Β ^ = tan −1 (V ^ Y / V ^ X ) (9) As described above, in the first embodiment, the vehicle state quantity estimating apparatus 10 and the yaw rate sensor 11 actual yaw rate γ detected longitudinal acceleration detected by the acceleration sensor 12 before and after alpha X, lateral acceleration alpha Y detected by the lateral acceleration sensor 13, the wheel speed Vw FL ~Vw detected by the wheel speed sensor 14
By performing the observer computation of the equation (7) and (8) based on the longitudinal direction velocity V X calculated based on the RR, actually detected the longitudinal estimated lateral velocity V Y direction velocity V ^ X by feeding back the error signal between the front-rear direction velocity V X which is the estimated lateral velocity V ^ Y is calculated, the estimated lateral slip angle beta ^ can be calculated on the basis of the this the estimated longitudinal direction velocity V ^ X .
【0034】このとき、前記(7)式および(8)式で
推定前後方向速度V^X と実際に検出した前後方向速度
VX との誤差信号をフィードバックするようにしている
だけであり、非線型な運動を線型化してモデル化する必
要がないので、低車速の急発進時などに車両の挙動が不
安定になり車両に大きな横滑り角を発生してスピンして
いる状態においても、より正確な横滑り角の推定を行う
ことができる。[0034] In this case, the equation (7) and (8) is only so as to feed back the error signal between the front-rear direction velocity V X which actually detects the longitudinal direction velocity V ^ X estimated equation, non Since there is no need to model linear motion in a linear fashion, the behavior of the vehicle becomes unstable when suddenly starting at low vehicle speeds, and even when the vehicle is spinning due to a large skid angle, it is more accurate. It is possible to estimate an appropriate side slip angle.
【0035】しかも、推定前後方向速度V^X と実際に
検出した前後方向速度VX との誤差信号をフィードバッ
クするので、推定値が初期値を誤って与えた等の何らか
の理由で真値から離れた場合にも真値に収束することに
なり、前述した第2従来例のように積分処理を利用した
方式におけるドリフトの発生を確実に防止することがで
きる。[0035] Moreover, since the fed back error signal actually the detected longitudinal direction velocity V X and the estimated longitudinal direction velocity V ^ X, away from the true value for any reason, such an estimate is given incorrectly initial value In such a case, the value converges to the true value, and it is possible to reliably prevent the occurrence of drift in the method using the integration processing as in the second conventional example.
【0036】次に、本発明の第2の実施形態を図3につ
いて説明する。この第2の実施形態は、前後方向速度V
X を車輪速度又はその推定値から求めるようにしたもの
である。Next, a second embodiment of the present invention will be described with reference to FIG. In the second embodiment, the longitudinal velocity V
X is obtained from the wheel speed or its estimated value.
【0037】すなわち、前輪又は後輪駆動車では非駆動
輪が存在するために、駆動時において非駆動輪の車輪速
を使用して前後方向速度を算出することが可能である。
一般に、低速時に大きな横滑り角が発生するのは、後輪
駆動車などが急発進を行った場合が想定される。そのよ
うな状態では、前輪には駆動力等が発生していないため
に車輪速度は概略その位置でのタイヤの前後方向の車体
速度と考えて差し支えない。このた、下記(10)式を
用いることで、車輪の回転速度より車輪位置での前後方
向速度を得ることができる。That is, since a non-drive wheel exists in a front-wheel or rear-wheel drive vehicle, it is possible to calculate the front-rear direction speed using the wheel speed of the non-drive wheel during driving.
In general, a large skid angle occurs at low speed when a rear-wheel drive vehicle or the like suddenly starts. In such a state, since no driving force or the like is generated at the front wheels, the wheel speed may be roughly considered to be the vehicle speed in the front-rear direction of the tire at that position. By using the following equation (10), it is possible to obtain the longitudinal speed at the wheel position from the rotational speed of the wheel.
【0038】 VwX =Vw・cosδ …………(10) ここで、Vwは車輪の回転速度を並進速度に換算したも
のであり、δは舵角(後輪では零又は4輪操舵システム
登載時にはその舵角となる)、VwX はVwの車体の前
後方向の成分であり、車輪のスリップが零の状態では車
輪位置での車体の前後方向速度となる。Vw X = Vw · cos δ (10) where Vw is a value obtained by converting the rotational speed of the wheel into a translation speed, and δ is a steering angle (zero or four-wheel steering system mounted on the rear wheels). Vw X is a component of Vw in the front-rear direction of the vehicle body, and when the wheel slip is zero, it is the vehicle front-rear speed at the wheel position.
【0039】また、前記第1の実施形態における(7)
および(8)式のオブザーバ演算で使用されるのは、重
心位置での前後方向速度であるが、これは下記(11)
式を用いて重心位置での前後方向速度に変換することで
算出することができる。In the first embodiment, (7)
In the observer calculation of equations (8) and (8), the velocity in the front-rear direction at the position of the center of gravity is used.
It can be calculated by converting into the longitudinal velocity at the position of the center of gravity using the equation.
【0040】 VX =VwX +(−1)i ・γ・tr/2 …………(11) ここで、trはトレッド、iは車輪位置を表す添え字で
奇数が左輪、偶数が右輪を夫々表す。V X = Vw X + (− 1) i · γ · tr / 2 (11) where tr is a tread, i is a subscript representing a wheel position, and odd numbers are left wheels, and even numbers are right. Each represents a ring.
【0041】そして、前記(7)および(8)式のオブ
ザーバ演算で必要な前後方向速度V X として上記(1
1)式を適用することにより、後輪駆動車が加速してい
る場合の横滑り角β^を算出することができる。Then, the objects of the equations (7) and (8)
Front-rear speed V required for server calculation X(1)
By applying equation (1), the rear-wheel drive vehicle is accelerating.
Can be calculated.
【0042】一方、制動時においては、車輪速の推定を
行うために、車両の四輪の内の一輪の制動力を意図的に
下げてスリップを抑制することにより、概略スリップが
ないときと同じと見なせることになり、前後方向速度V
wX を算出することが可能となる。また、全ての車輪に
おいて制動力によって車輪のスリップが発生している状
態でも、例えば車輪の回転運動のモデルを使用し、車輪
速度を他の手段によりって推定することによって前記
(7)および(8)式のオブザーバ演算で必要な前後方
向速度VX を算出することができる。ここで、タイヤの
摩擦係数は摩擦係数推定手段によって与えられる。On the other hand, at the time of braking, in order to estimate the wheel speed, the braking force of one of the four wheels of the vehicle is intentionally reduced to suppress the slip. And the speed V in the front-rear direction
it is possible to calculate the w X. Further, even in a state where the wheels are slipping due to the braking force in all the wheels, the wheel speed is estimated by other means by using, for example, a model of the rotational motion of the wheels, and the above (7) and ( 8) can be calculated longitudinal direction velocity V X required by the observer calculation formula. Here, the friction coefficient of the tire is given by the friction coefficient estimating means.
【0043】このように加速状態および減速状態の何れ
の状態においても、車輪の回転を使用して、車輪位置若
しくは重心位置での前後方向速度を求めることができ
る。具体的には、図3に示すように、車両状態量推定装
置10に、前述したヨーレートセンサ11、前後加速度
センサ12、横加速度センサ13及び車輪速センサ14
に加えて、前輪1Fの操舵角を検出する操舵角センサ1
5で検出した操舵角δ及びマスタシリンダ圧を検出する
マスタシリンダ圧センサ16で検出したマスタシリンダ
圧MPが入力されると共に、制動力制御装置30からの
制動力制御状態であるか否かを表す制動力制御状態信号
SBが入力され、且つこの車両状態量推定装置10で、
図4に示す状態量演算処理を実行することにより、車輪
速度に基づいて推定横滑り角β^を算出する。As described above, in both the acceleration state and the deceleration state, the rotation in the forward and backward directions at the wheel position or the center of gravity can be obtained by using the rotation of the wheel. Specifically, as shown in FIG. 3, the vehicle state quantity estimating apparatus 10 includes the above-described yaw rate sensor 11, longitudinal acceleration sensor 12, lateral acceleration sensor 13, and wheel speed sensor 14 as described above.
Angle sensor 1 for detecting the steering angle of front wheel 1F
5, the master cylinder pressure MP detected by the master cylinder pressure sensor 16 for detecting the steering angle δ and the master cylinder pressure is input, and indicates whether or not the braking force control unit 30 is in the braking force control state. The braking force control state signal SB is input, and the vehicle state quantity estimation device 10
The estimated side slip angle β ^ is calculated based on the wheel speed by executing the state quantity calculation processing shown in FIG.
【0044】この状態量演算処理は、所定時間(例えば
10msec)毎のタイマ割込処理によって実行され、
先ずステップS1でマスタシリンダ圧センサ16で検出
したマスタシリンダ圧MPを読込み、これが予め設定し
た比較的小さい1MPa程度の閾値MPS 以上であるか
否かを判定することにより、運転者がブレーキペダルを
踏込んだ制動状態であるか否かを判定し、MP<MPS
であるときには運転者による制動状態でないと判断し
て、ステップS2に移行し、制動力制御装置30からの
制動力制御状態信号SBを読込んで、これがオン状態で
あるか否かを判定することにより、ブレーキペダルの踏
込以外のトラクション制御又はヨーモーメント抑制制御
による制動力制御が行われているか否かを判定し、制動
力制御が行われていないときにはステップS3に移行し
て、非駆動輪となる前輪の左右何れか一方の車輪速度V
wFL,VwFRを車輪速度Vwとして設定してからステッ
プS4に移行する。This state quantity calculation processing is executed by a timer interruption processing every predetermined time (for example, 10 msec).
First read the master cylinder pressure MP detected by the master cylinder pressure sensor 16 in step S1, by which to determine whether a relatively small 1MPa about threshold MP S above a preset, the driver on the brake pedal It is determined whether the vehicle is in the depressed braking state and MP <MP S
When it is determined that the vehicle is not in the braking state by the driver, the process proceeds to step S2, the braking force control state signal SB from the braking force control device 30 is read, and it is determined whether or not this is the on state. It is determined whether or not the braking force control by the traction control other than the depression of the brake pedal or the yaw moment suppression control is being performed, and if the braking force control is not being performed, the process proceeds to step S3 to be the non-drive wheel. Either left or right wheel speed V of the front wheel
After setting w FL and Vw FR as the wheel speed Vw, the process proceeds to step S4.
【0045】このステップS4では、車輪速度Vwをも
とに前記(10)式の演算を行って前後方向速度VwX
を算出し、次いでステップS5に移行して、ヨーレート
センサ11で検出したヨーレートγ(n) を読込んでから
ステップS6に移行する。In step S4, the equation (10) is calculated on the basis of the wheel speed Vw, and the longitudinal speed Vw X
Then, the process proceeds to step S5 to read the yaw rate γ (n) detected by the yaw rate sensor 11, and then proceeds to step S6.
【0046】このステップS6では、前後方向速度Vw
X 及びヨーレートγ(n) をもとに前記(11)式の演算
を行って、重心位置での前後方向加速度VX (n) を算出
する。In this step S6, the longitudinal speed Vw
Based on X and the yaw rate γ (n), the calculation of the equation (11) is performed to calculate the longitudinal acceleration V X (n) at the position of the center of gravity.
【0047】次いで、ステップS7に移行して、重心位
置での前後方向速度VX (n) をもとに下記(12)式及
び(13)式に従って車両の前後方向推定車体速度V^
X (n) 及び横方向推定車体速度V^Y (n) を算出する。Next, the process proceeds to step S7, and based on the longitudinal speed V X (n) at the position of the center of gravity, the estimated longitudinal vehicle speed V ^ of the vehicle in accordance with the following equations (12) and (13).
X (n) and estimated lateral vehicle speed V ^ Y (n) are calculated.
【0048】 V^X (n) =V^X (n) +Δt(−A+αX (n) −B−C) ……(12) V^Y (n) =V^Y (n) +Δt(D+αY (n) −E−F) ……(12) A=V^Y (n-1) γ(n-1) B=K11(V^X (n-1) −VX (n) ) C=K12sign(V^X (n-1) −VX (n) ) D=V^X (n-1) γ(n-1) E=K21(V^X (n-1) −VX (n) ) F=K22sign(V^X (n-1) −VX (n) ) 次いで、ステップS8に移行して、車両の前後方向推定
車体速度V^X (n) 及び横方向推定車体速度V^Y (n)
をもとに、前記(9)式に従って推定横滑り角β^(n)
を算出し、次いでステップS9に移行して、算出した車
両の前後方向推定車体速度V^X (n) 及び横方向推定車
体速度V^Y (n) を前回値V^X (n-1)及びV^Y (n-1)
として記憶してからタイマ割込処理を終了して所定の
メインプログラムに復帰する。[0048] V ^ X (n) = V ^ X (n) + Δt (-A + α X (n) -B-C) ...... (12) V ^ Y (n) = V ^ Y (n) + Δt (D + α Y (n) -E-F) ...... (12) A = V ^ Y (n-1) γ (n-1) B = K 11 (V ^ X (n-1) -V X (n)) C = K 12 sign (V ^ X (n-1) -V X (n)) D = V ^ X (n-1) γ (n-1) E = K 21 (V ^ X (n-1) −V X (n)) F = K 22 sign (V ^ X (n−1) −V X (n)) Then, the process proceeds to step S8, and the estimated vehicle speed V 前後X (n) in the front-rear direction of the vehicle is obtained. And estimated lateral vehicle speed V ^ Y (n)
And the estimated sideslip angle β ^ (n) according to the above equation (9).
Then, the process proceeds to step S9, where the calculated estimated longitudinal vehicle body speed V ( X (n) and estimated lateral vehicle speed V ^ Y (n) of the vehicle are set to the previous values V ^ X (n-1). And V ^ Y (n-1)
And then terminates the timer interrupt processing and returns to the predetermined main program.
【0049】一方、ステップS1の判定結果が、運転者
がブレーキペダルを踏込んでいる制動状態であるときに
は、ステップS10に移行して、例えばアンチロックブ
レーキ制御装置で使用する推定車体速度を参照して、こ
れを車輪速度Vwとして設定してから前記ステップS4
に移行する。On the other hand, if the result of the determination in step S1 is that the driver is depressing the brake pedal, the flow proceeds to step S10 to refer to, for example, the estimated vehicle speed used by the antilock brake control device. , This is set as the wheel speed Vw, and
Move to
【0050】また、ステップS2の判定結果が、車両が
制動力制御を行っているときには、ステップS11に移
行して、非駆動輪に制動制御力を解除可能な車輪がある
か否かを判定し、解除可能な車輪があるときにはステッ
プS12に移行して、該当車輪の制御力を解除してから
該当車輪の車輪速を車輪速度Vwとして設定してから前
記ステップS4に移行し、解除可能な車輪がないときに
は前記ステップS10に移行する。If the result of the determination in step S2 is that the vehicle is performing braking force control, the flow shifts to step S11 to determine whether or not any of the non-driven wheels has a wheel capable of releasing the braking control force. If there is a wheel that can be released, the process proceeds to step S12, in which the control force of the wheel is released, the wheel speed of the wheel is set as the wheel speed Vw, and then the process proceeds to step S4. If there is not, the process proceeds to step S10.
【0051】この図4の横滑り角演算処理において、ス
テップS1〜S6の処理前後方向速度算出手段に対応
し、ステップS7の処理が速度推定手段に対応し、ステ
ップS8の処理が横滑り角算出手段に対応している。In the side slip angle calculation processing of FIG. 4, the processing in steps S1 to S6 corresponds to the forward / backward speed calculation means, the processing in step S7 corresponds to the speed estimation means, and the processing in step S8 corresponds to the side slip angle calculation means. Yes, it is.
【0052】したがって、運転者がブレーキペダルを踏
込んでおらず、トラクション制御やヨーレート制御等の
制動力制御が行われていないときには、非駆動輪となる
前輪1Fの左右輪の何れか一方の車輪速VwFL又はVw
FRに基づいて重心位置における前後方向速度VX (n) を
算出し、これと前回の推定前後方向速度V^X (n-1)及
び推定横方向速度V^Y (n-1) と、前後加速度αX 及び
横加速度αY に基づくオブザーバ演算によって推定前後
方向速度V^X (n) 及び推定横方向速度V^Y(n) を算
出し、これらに応じて推定横滑り角β^(n) を算出す
る。Therefore, when the driver has not stepped on the brake pedal and braking force control such as traction control or yaw rate control has not been performed, the wheel speed of one of the left and right wheels of the front wheel 1F, which is a non-driven wheel, is determined. Vw FL or Vw
The longitudinal velocity V X (n) at the position of the center of gravity is calculated based on FR , and this and the previous estimated longitudinal velocity V ^ X (n-1) and the estimated lateral velocity V ^ Y (n-1) are calculated as follows: calculating the estimated longitudinal direction velocity V ^ X (n) and the estimated lateral velocity V ^ Y (n) by an observer computation based on the longitudinal acceleration alpha X and the lateral acceleration alpha Y, estimated lateral slip angle in accordance with these beta ^ (n ) Is calculated.
【0053】一方、運転者が制動操作を行っていない
が、トラクション制御やヨーレート制御が行われている
場合には、トラクション制御では非駆動輪となる前輪で
は制動力制御されていないので、これらの何れかの車輪
速が選択され、ヨーレート制御では内輪側又は外輪側の
何れかの車輪が制動力制御されないので、この車輪が選
択されて、これらに基づいて重心位置における前後方向
速度VX (n) を算出する。On the other hand, when the driver does not perform the braking operation but performs the traction control or the yaw rate control, the traction control does not perform the braking force control on the front wheels that are the non-driving wheels. Any one of the wheel speeds is selected, and in the yaw rate control, the braking force control is not performed on either the inner wheel side or the outer wheel side. Therefore, this wheel is selected, and the forward / rearward speed V X (n ) Is calculated.
【0054】また、運転者が制動操作を行っている場合
や制動操作を行っていないがトラクション制御やヨーレ
ート制御によって全ての車輪が制動力制御されている場
合には、例えばアンチロックブレーキ制御での推定車体
速度に基づいて重心位置における前後方向速度を算出す
る。When the driver is performing the braking operation or when the braking operation is not performed but all the wheels are controlled by the traction control or the yaw rate control, for example, the antilock brake control The longitudinal speed at the position of the center of gravity is calculated based on the estimated vehicle speed.
【0055】この結果、車両の制動状態に応じてそのと
きの前後方向速度を的確に表す車輪速を選択して、これ
に基づいて重心位置おける前後方向速度を算出すること
ができ、高精度で横滑り角を推定することができる。As a result, it is possible to select the wheel speed accurately representing the longitudinal speed at that time according to the braking state of the vehicle and calculate the longitudinal speed at the center of gravity based on the selected wheel speed. The sideslip angle can be estimated.
【0056】なお、上記第2の実施形態においては、運
転者の制動操作をマスタシリンダ圧センサ16でマスタ
シリンダ圧を検出することにより検出する場合について
説明したが、これに限定されるものではなく、ブレーキ
ペダルの踏込み時にオンとなるブレーキランプスイッチ
のスイッチ信号等によって検出するようにしてもよい。In the second embodiment, the case where the braking operation of the driver is detected by detecting the master cylinder pressure by the master cylinder pressure sensor 16 has been described. However, the present invention is not limited to this. Alternatively, it may be detected by a switch signal of a brake lamp switch that is turned on when the brake pedal is depressed.
【0057】次に、本発明の第3の実施形態を図5につ
いて説明する。この第3の実施形態は、前述した第1及
び第2の実施形態においては、推定横滑り角を算出する
際に、操舵角δを使用していないので、低車速域での前
後方向加速度αX に対して横方向加速度αY が大きい状
態では高精度で横滑り角を推定することができるが、そ
れ以外の高車速域では、推定精度が低下することになる
ため、高車速域では他の推定方法に切換えることによ
り、全ての車速域で高精度の横滑り角推定を行うことが
できるようにしたものである。Next, a third embodiment of the present invention will be described with reference to FIG. In the third embodiment, since the steering angle δ is not used when calculating the estimated sideslip angle in the first and second embodiments, the longitudinal acceleration α X in the low vehicle speed range is used. When the lateral acceleration α Y is large, the sideslip angle can be estimated with high accuracy, but in other high vehicle speed ranges, the estimation accuracy decreases. By switching to the method, it is possible to perform a highly accurate side slip angle estimation in all vehicle speed ranges.
【0058】すなわち、車両状態量推定装置10で、図
5に示す推定方法切換処理を実行する。この推定方法切
換処理は、所定時間(例えば10msec)毎のタイマ
割込処理として実行され、先ず、ステップS21で、前
回使用した推定方法が前述した第1従来例と同様の運動
モデルを使用した高速用横滑り角推定方法であったか否
かを判定し、高速用横滑り角推定方法であったときに
は、ステップS22に移行して、駆動輪としての後輪1
Rの左右輪のうちの高い車輪速VwRHから非駆動輪とし
ての前輪1Fの左右輪のうちの高い車輪速VwFHを減算
した値(VwRH−VwFH)が予め設定した設定値ΔVw
1 より大きく且つ非駆動輪としての前輪1Fの高い車輪
速VwFHが予め設定した設定車速VwS1より小さいか否
かを判定し、(VwRH−VwFH)>ΔVw1 且つVwFH
<VwS1であるときには、低車速域での駆動輪スピン状
態であると判断して、ステップS23に移行する。That is, the vehicle state quantity estimating apparatus 10 executes an estimation method switching process shown in FIG. The estimation method switching process is executed as a timer interrupt process at predetermined time intervals (for example, 10 msec). First, in step S21, the estimation method used last time uses the same high-speed motion model as the first conventional example described above. It is determined whether or not the method is a side slip angle estimating method for high speed. If the method is a method for estimating a side slip angle for high speed, the process proceeds to step S22, and the rear wheel 1 as a driving wheel is determined.
A value (Vw RH −Vw FH ) obtained by subtracting the high wheel speed Vw FH of the left and right wheels of the front wheel 1F as the non-drive wheel from the high wheel speed Vw RH of the left and right wheels of R is a preset value ΔVw.
It is determined whether or not the wheel speed Vw FH greater than 1 and the high speed of the front wheel 1F as the non-drive wheel is smaller than a preset vehicle speed Vw S1 , and (Vw RH −Vw FH )> ΔVw 1 and Vw FH
When <Vw S1 , it is determined that the driving wheel is spinning in the low vehicle speed range, and the process proceeds to step S23.
【0059】このステップS23では、前述した第1の
実施形態又は第2の実施形態による前後方向速度V
X (n) と推定前後方向速度V^X (n) との誤差信号をフ
ィードバックすることにより、推定前後方向速度V^X
(n) 及び推定横方向速度V^Y (n) を算出し、これらか
ら推定横滑り角β^(n) を算出する横滑り角推定方法を
使用する低速用横滑り角推定方法によって推定横滑り角
β^(n) を算出してからステップS24に移行する。In step S23, the longitudinal velocity V according to the first or second embodiment described above is used.
By feeding back an error signal between X (n) and the estimated longitudinal velocity V ^ X (n), the estimated longitudinal velocity V ^ X
(n) and the estimated lateral velocity V ^ Y (n), and the estimated sideslip angle β ^ by the low-speed sideslip angle estimation method using the sideslip angle estimation method of calculating the estimated sideslip angle β ^ (n) from these. After calculating (n), the process proceeds to step S24.
【0060】このステップS24では、後述する高車速
用横滑り角推定方法で推定ヨーレートγ^(n) を使用す
るので、ヨーレートセンサ11で検出した実ヨーレート
γ(n) を推定ヨーレートγ^(n) として記憶してからス
テップS25に移行して、今回の推定結果(γ(n) 、β
(n) 、VX (n) 、VY (n) )を前回値として記憶すると
共に、使用した横滑り推定方法に応じた種別を表すフラ
グFを設定してからタイマ割込処理を終了して所定のメ
インプログラムに復帰する。In step S24, since the estimated yaw rate γ ^ (n) is used in the high vehicle speed side slip angle estimation method described later, the actual yaw rate γ (n) detected by the yaw rate sensor 11 is used as the estimated yaw rate γ ^ (n). Then, the process proceeds to step S25, where the estimation results (γ (n), β
(n), V X (n), and V Y (n)) are stored as previous values, and a flag F indicating a type corresponding to the side slip estimation method used is set, and then the timer interrupt processing is terminated. Return to the predetermined main program.
【0061】一方、ステップS22の判定結果が設定条
件を満足しないときには、ステップS26に移行して、
高速用横滑り角推定方法によって推定横滑り角β^(n)
を算出し、次いでステップS27に移行して、ステップ
S23での低速用横滑り角推定方法で使用する推定前後
方向速度V^X として車速V(n) を設定すると共に、下
記(13)式に従って推定横方向速度V^Y (n) を算出
してから前記ステップS25に移行する。On the other hand, if the result of the determination in step S22 does not satisfy the set condition, the flow shifts to step S26,
Estimated sideslip angle β ^ (n) by sideslip angle estimation method for high speed
Then, the process proceeds to step S27, where the vehicle speed V (n) is set as the estimated front-rear direction speed V ^ X used in the low-speed side slip angle estimation method in step S23, and is estimated according to the following equation (13). After calculating the lateral velocity V ^ Y (n), the flow shifts to step S25.
【0062】 V^X (n) =V(n) *tan(β^(n) ) …………(13) また、ステップS21の判定結果が前回の推定方法が低
速用横滑り角推定方法であるときには、ステップS28
に移行して、前記ステップS22と同様に駆動輪として
の後輪1Rの左右輪のうちの高い車輪速VwRHから非駆
動輪としての前輪1Fの左右輪のうちの高い車輪速Vw
FHを減算した値(VwRH−VwFH)が予め設定した前記
設定値ΔVw1 より大きい設定値ΔVw2 より大きく且
つ非駆動輪としての前輪1Fの高い車輪速VwFHが予め
設定した前記設定車速VwS1よ大きい設定車速VwS2よ
り大きいか否かを判定し、(VwRH−VwFH)>ΔVw
1且つVwFH>VwS2であるときには、高車速域である
か駆動輪にホイールスピンを生じていないものと判断し
て前記ステップS26に移行し、そうでないときには駆
動輪にホイールスピンを生じているものと判断して前記
ステップS23に移行する。V ^ X (n) = V (n) * tan (β ^ (n)) (13) Also, the result of the determination in step S21 is that the previous estimation method is the low-speed sideslip angle estimation method. If there is, step S28
And the high wheel speed Vw of the left and right wheels of the front wheel 1F as the non-driving wheel from the high wheel speed Vw RH of the left and right wheels of the rear wheel 1R as the driving wheel, similarly to the step S22.
The value (Vw RH -Vw FH ) obtained by subtracting FH is larger than the set value ΔVw 2 which is larger than the set value ΔVw 1 set in advance, and the high wheel speed Vw FH of the front wheel 1F as the non-drive wheel is set in advance to the set vehicle speed determining whether greater than Vw S1 by large set vehicle speed Vw S2, (Vw RH -Vw FH )> ΔVw
If 1 and Vw FH > Vw S2 , it is determined that the vehicle is in the high vehicle speed range or that no wheel spin has occurred in the drive wheels, and the process proceeds to step S26. Otherwise, wheel spin has occurred in the drive wheels. Then, the process proceeds to step S23.
【0063】また、ステップS26での高速用横滑り角
推定方法は、前述した図2の2輪運動モデルから前輪の
タイヤの横滑り角βf 及び後輪のタイヤの横滑り角βr
は、下記(14)式及び(15)式で表すことができ
る。The method for estimating the side slip angle for high speed in step S26 is based on the side slip angle β f of the front wheel tire and the side slip angle β r of the rear wheel tire based on the two-wheel motion model shown in FIG.
Can be expressed by the following equations (14) and (15).
【0064】 βf =tan-1{(VY −Lf γ)/VX }+δ …………(14) βr =tan-1{(VY −Lr γ)/VX } …………(15) さらに、車速VX がある程度大きい状態では線形化が可
能であり、下記(16)式で表される車体横滑り角βと
車体前後方向速度VX 及び横方向速度VY との関係を用
い、前記(14)式及び(15)式に代えて下記(1
7)式及び(18)式の関係が得られ、以上より最終的
に下記(19)式及び(20)式の線形2輪モデルが得
られる。Β f = tan -1 {(V Y -L f γ) / V X } + δ (14) β r = tan -1 {(V Y -L r γ) / V X … ......... (15) further, in a state the vehicle speed V X is large to some extent is possible linearization, a vehicle body side slip angle β and the longitudinal direction of the vehicle body speed V X and lateral velocity V Y represented by the following formula (16) And the following equation (1) is used instead of the equations (14) and (15).
The relations of the equations (7) and (18) are obtained, and finally, the linear two-wheel model of the following equations (19) and (20) is obtained.
【0065】 β=−VY /VX …………(16) βf =−β−(Lf γ/VX )+δ …………(17) βr =−β+(Lr γ/VX ) …………(18) MVX (β′+γ)=2CPf {−β−(Lf γ/VX )+δ} +2CPr {−β+(Lr γ/VX )}……(19) IZ γ′=2Lf CPf {−β−(Lf γ/VX )+δ} −2Lr CPr {−β+(Lr γ/VX )} ……(20) この線形2輪モデルは制御設計を行うために下記(2
1)〜(28)の状態空間で表現できる。Β = −V Y / V X (16) β f = −β− (L f γ / V X ) + δ (17) β r = −β + (L r γ / V X ) (18) MV X (β ′ + γ) = 2CP f {−β− (L f γ / V X ) + δ} + 2CP r {−β + (L r γ / V X )} (19) I Z γ '= 2L f CP f {-β- (L f γ / V X) + δ} -2L r CP r {-β + (L r γ / V X)} ...... (20) the linear The two-wheel model is as follows (2
It can be represented by the state space of 1) to (28).
【0066】[0066]
【数1】 (Equation 1)
【0067】そして、上記(21)式及び(22)式の
車両モデルに対し、オブザーバは下記(30)式ように
設計される。 x^′=Ax^+Bu+Ly〜 …………(30) ここで、A,B,C,Dは車両の運動を計算するための
行列である。x^は、推定された状態量、Lは状態推定
量の補正量を決定するゲインである。ここで、誤差信号
y〜は、下記(31)式によって求められ、推定状態量
x^は下記(32)式で求められる。The observer is designed as shown in the following equation (30) for the vehicle models of the above equations (21) and (22). x ^ ′ = Ax ^ + Bu + Ly (30) where A, B, C, and D are matrices for calculating the motion of the vehicle. x ^ is an estimated state amount, and L is a gain that determines a correction amount of the state estimated amount. Here, the error signal y〜 is obtained by the following equation (31), and the estimated state quantity x ^ is obtained by the following equation (32).
【0068】 y〜=y−(Cx^+Du) …………(31)Y〜 = y− (Cx ^ + Du) (31)
【0069】[0069]
【数2】 (Equation 2)
【0070】したがって、図6におけるステップS31
で前記(25)式〜(28)式の行列A〜Dを計算し、
次いでステップS32に移行して、前記(29)式及び
(30)式を使用して推定ヨーレートγ^及び推定横滑
り量β^を算出し、次いでステップS33に移行して、
算出した推定ヨーレートγ^及び推定横滑り量β^を前
回値として記憶してから処理を終了する。Therefore, step S31 in FIG.
Calculate the matrices A to D of the above equations (25) to (28),
Then, the process proceeds to step S32 to calculate the estimated yaw rate γ ^ and the estimated side slip amount β ^ using the above equations (29) and (30), and then proceeds to step S33.
After the calculated estimated yaw rate γ ^ and the estimated side slip amount β ^ are stored as the previous values, the process ends.
【0071】この第3の実施形態によると、車両の発進
時等の低車速域で駆動輪となる後輪にホイールスピンを
生じる状態となると、前回の処理時に高速用横滑り角推
定方法が選択されているものとしても、ステップS22
からステップS23に移行して、第1又は第2の実施形
態による前後方向速度VX (n) と推定前後方向速度V^
X (n) との誤差信号をフィードバックすることにより、
推定前後方向速度V^ X (n) 及び推定横方向速度V^Y
(n) を算出し、これらから推定横滑り角β^(n) を算出
する横滑り角推定方法を使用する低速用横滑り角推定方
法によって推定横滑り角β^(n) を算出するので、高精
度で推定横滑り角β^を算出することができる。According to the third embodiment, the vehicle starts
Wheel spin on the rear wheel that is the driving wheel in low vehicle speed range such as time
If this occurs, the high-speed sideslip angle
Even if the setting method is selected, step S22
From step S23 to the first or second embodiment
Velocity V depending on the stateX(n) and estimated longitudinal velocity V 方向
XBy feeding back the error signal with (n),
Estimated velocity V 前後 X(n) and estimated lateral velocity V ^Y
(n) and the estimated sideslip angle β ^ (n)
Estimation Method for Low-Speed Side-Slip Angle Using Evolving Side-Slip Angle Estimation Method
The estimated sideslip angle β ^ (n) is calculated by the
The estimated sideslip angle β ^ can be calculated in degrees.
【0072】この状態から、駆動輪としての後輪のホイ
ールスピンが収まるか、車速が速くなると、ステップS
28での判定条件を満足しなくなると、ステップS26
に移行して、図8に示す高速用横滑り角推定処理が実行
されるので、高車速域で高精度の横滑り角推定を行うこ
とができる。From this state, if the wheel spin of the rear wheel as the driving wheel stops or the vehicle speed increases, step S
If the determination condition in step 28 is no longer satisfied, step S26
Then, the high-speed sideslip angle estimation process shown in FIG. 8 is executed, so that a highly accurate sideslip angle estimation can be performed in a high vehicle speed range.
【0073】なお、上記第3の実施形態においては、高
速用横滑り推定方法として、(25)式〜(31)式を
使用して推定ヨーレートγ^及び推定横滑り角β^を算
出する場合について説明したが、これに限定されるもの
ではなく、前述した第1従来例や第2従来例のような他
の横滑り角推定方法を適用することができる。In the third embodiment, the case where the estimated yaw rate γ ^ and the estimated side slip angle β ^ are calculated using the equations (25) to (31) as the high-speed sideslip estimation method will be described. However, the present invention is not limited to this, and other methods for estimating the sideslip angle, such as the first conventional example and the second conventional example described above, can be applied.
【図1】本発明の第1の実施形態の全体構成を示すブロ
ック図である。FIG. 1 is a block diagram showing an overall configuration of a first embodiment of the present invention.
【図2】2輪運動モデルを示す説明図である。FIG. 2 is an explanatory diagram showing a two-wheel motion model.
【図3】本発明の第2の実施形態の全体構成を示すブロ
ック図である。FIG. 3 is a block diagram showing an overall configuration of a second embodiment of the present invention.
【図4】第2の実施の形態における演算処理の概要を示
すフローチャートである。FIG. 4 is a flowchart illustrating an outline of a calculation process according to the second embodiment.
【図5】本発明の第3の実施形態における演算処理の概
要を示すフローチャートである。FIG. 5 is a flowchart illustrating an outline of a calculation process according to a third embodiment of the present invention.
【図6】第3の実施形態における高速用横滑り角演算処
理の一例を示すフローチャートである。FIG. 6 is a flowchart illustrating an example of a high-speed sideslip angle calculation process according to the third embodiment.
10 車両状態量推定装置 11 ヨーレートセンサ 12 前後加速度センサ 13 横加速度センサ 14 車輪速センサ 15 操舵角センサ 16 マスタシリンダ圧センサ 20 目標制動力設定装置 30 制動力制御装置 Reference Signs List 10 vehicle state quantity estimation device 11 yaw rate sensor 12 longitudinal acceleration sensor 13 lateral acceleration sensor 14 wheel speed sensor 15 steering angle sensor 16 master cylinder pressure sensor 20 target braking force setting device 30 braking force control device
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.7 識別記号 FI テーマコート゛(参考) B62D 113:00 137:00 ──────────────────────────────────────────────────続 き Continued on the front page (51) Int.Cl. 7 Identification symbol FI theme coat ゛ (reference) B62D 113: 00 137: 00
Claims (5)
々当該車両が鉛直軸周りの運動をしている状態を記述す
る運動モデルを有し、該運動モデルは、前後方向の運動
を記述する運動モデルから推定される推定前後運動物理
量と、実際に測定される前後運動物理量及び他の推定方
法から得られる前後運動物理量の何れかとの誤差信号を
フィードバック信号として使用し、前記運動モデルに基
づいて車体の横滑り角を推定するようにしたことを特徴
とする車両の車体横滑り角推定方法。1. A motion model that describes a state in which the vehicle is moving around a vertical axis in a front-rear direction and a lateral direction of the vehicle, wherein the motion model describes a motion in a front-rear direction. Using an error signal between the estimated longitudinal physical quantity estimated from the model and any of the actually measured longitudinal physical quantity and the longitudinal physical quantity obtained from another estimation method as a feedback signal, the vehicle body based on the motion model A vehicle side slip angle estimation method for estimating a vehicle side slip angle.
度が設定されていることを特徴とする請求項1記載の車
両の車体横滑り角推定方法。2. The method according to claim 1, wherein a longitudinal speed is set as the longitudinal physical quantity.
速度の推定値の何れかを使用することを特徴とする請求
項2記載の車両の車体横滑り角推定方法。3. The vehicle body side slip angle estimating method according to claim 2, wherein the longitudinal speed uses one of a wheel speed and an estimated value of the wheel speed.
動輪のそれに対してある一定値又は一定比率以上大き
く、且つ前輪の車輪速又は車体速度がある一定値より小
さい場合に請求項1の車体横滑り角推定方法を使用し、
それ以外の場合に他の車体横滑り角推定方法を使用する
ようにしたことを特徴とする車両の車体横滑り角推定方
法。4. The method according to claim 1, wherein the larger value of the wheel speed of the driving wheel is larger than that of the non-driving wheel by a certain value or a certain ratio or more, and the wheel speed of the front wheel or the vehicle speed is smaller than a certain value. Using the vehicle body side slip angle estimation method of item 1,
A method for estimating a vehicle body side slip angle of a vehicle, wherein another method for estimating a vehicle body side slip angle is used in other cases.
速度算出手段と、車体に発生するヨーレートを検出する
ヨーレート検出手段と、前記車輪速度算出手段で検出し
た車輪速度と前記ヨーレート検出手段で検出したヨーレ
ートとに基づいて前後方向速度を算出する前後方向速度
算出手段と、該前後方向速度算出手段で算出した前後方
向速度と推定前後方向速度との誤差信号をフィードバッ
クして前後方向及び横方向速度を推定する速度推定手段
と、該速度推定手段で推定した前後方向及び横方向速度
に基づいて横滑り角を算出する横滑り角算出手段とを備
えたことを特徴とする車両の車体横滑り角推定装置。5. A wheel speed calculating means for calculating a wheel speed of each wheel of a vehicle, a yaw rate detecting means for detecting a yaw rate generated in the vehicle body, and a wheel speed detected by the wheel speed calculating means and the yaw rate detecting means. Forward / backward speed calculating means for calculating a forward / backward speed based on the detected yaw rate; and an error signal between the forward / backward speed calculated by the forward / backward speed calculating means and the estimated forward / backward speed, fed back and forth and forward and backward. A vehicle body slip angle estimating device comprising: a speed estimating means for estimating a speed; and a skid angle calculating means for calculating a skid angle based on the longitudinal and lateral speeds estimated by the speed estimating means. .
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP22724298A JP3959860B2 (en) | 1998-08-11 | 1998-08-11 | Vehicle body slip angle estimation method and estimation apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP22724298A JP3959860B2 (en) | 1998-08-11 | 1998-08-11 | Vehicle body slip angle estimation method and estimation apparatus |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2000052951A true JP2000052951A (en) | 2000-02-22 |
| JP3959860B2 JP3959860B2 (en) | 2007-08-15 |
Family
ID=16857755
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP22724298A Expired - Fee Related JP3959860B2 (en) | 1998-08-11 | 1998-08-11 | Vehicle body slip angle estimation method and estimation apparatus |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3959860B2 (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003531066A (en) * | 2000-04-19 | 2003-10-21 | コンティネンタル・テーベス・アクチエンゲゼルシヤフト・ウント・コンパニー・オッフェネ・ハンデルスゲゼルシヤフト | How to determine the amount of driving dynamics for cars online |
| JP2007106273A (en) * | 2005-10-14 | 2007-04-26 | Nissan Motor Co Ltd | Vehicle running state estimation device |
| JP2009053039A (en) * | 2007-08-27 | 2009-03-12 | Honda Motor Co Ltd | Vehicle posture estimation device and vehicle posture estimation method |
| JP2009067300A (en) * | 2007-09-14 | 2009-04-02 | Toyota Central R&D Labs Inc | Body speed estimation device |
-
1998
- 1998-08-11 JP JP22724298A patent/JP3959860B2/en not_active Expired - Fee Related
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003531066A (en) * | 2000-04-19 | 2003-10-21 | コンティネンタル・テーベス・アクチエンゲゼルシヤフト・ウント・コンパニー・オッフェネ・ハンデルスゲゼルシヤフト | How to determine the amount of driving dynamics for cars online |
| JP2007106273A (en) * | 2005-10-14 | 2007-04-26 | Nissan Motor Co Ltd | Vehicle running state estimation device |
| JP2009053039A (en) * | 2007-08-27 | 2009-03-12 | Honda Motor Co Ltd | Vehicle posture estimation device and vehicle posture estimation method |
| JP2009067300A (en) * | 2007-09-14 | 2009-04-02 | Toyota Central R&D Labs Inc | Body speed estimation device |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3959860B2 (en) | 2007-08-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3394249B2 (en) | Vehicle motion control method | |
| US7496436B2 (en) | Method and device for assisting an operator of a vehicle in the vehicle stabilization | |
| JP3669668B2 (en) | Vehicle wheel slip angle detection device | |
| JP3539722B2 (en) | Road surface friction coefficient estimation device for vehicles | |
| EP0982206B1 (en) | Method of estimating vehicle yaw rate | |
| JP2004505855A (en) | Method for controlling yaw and lateral dynamics in road vehicles | |
| CN116513161B (en) | Vehicle stability control method based on drift balance | |
| JP3748334B2 (en) | Vehicle attitude control device | |
| JP4990384B2 (en) | Vehicle motion control method using jerk information | |
| JPH06221968A (en) | Road friction coefficient detector | |
| JP3271956B2 (en) | Road surface friction coefficient estimation device for vehicles | |
| JP3959860B2 (en) | Vehicle body slip angle estimation method and estimation apparatus | |
| US6349256B1 (en) | Turning behavior state detecting system for vehicle | |
| US6974195B2 (en) | Method for increasing the maneuverability or driving stability of a vehicle during cornering | |
| JP3039071B2 (en) | Vehicle turning limit judgment device | |
| JP2002162225A (en) | Road surface gradient estimator | |
| KR100721046B1 (en) | Vehicle stability control system | |
| KR100907868B1 (en) | Control method of vehicle stability control system | |
| JP3535358B2 (en) | Road friction coefficient estimation device | |
| JPH11115720A (en) | Estimation device of road surface friction coefficient | |
| JP3669174B2 (en) | Vehicle motion control device | |
| JP3087438B2 (en) | Vehicle turning state estimation device and vehicle motion control device | |
| JP3547956B2 (en) | Vehicle attitude control device | |
| JPH1178933A (en) | Method and apparatus for estimating vehicle sideslip angle of vehicle | |
| KR100751216B1 (en) | Vehicle stability control system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050822 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060110 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060310 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060815 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061013 |
|
| 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: 20070424 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070507 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |