JPH10181565A - Road surface state identification device - Google Patents

Road surface state identification device

Info

Publication number
JPH10181565A
JPH10181565A JP34361296A JP34361296A JPH10181565A JP H10181565 A JPH10181565 A JP H10181565A JP 34361296 A JP34361296 A JP 34361296A JP 34361296 A JP34361296 A JP 34361296A JP H10181565 A JPH10181565 A JP H10181565A
Authority
JP
Japan
Prior art keywords
tire
road surface
rotation
signal
state
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP34361296A
Other languages
Japanese (ja)
Other versions
JP3684726B2 (en
Inventor
Nobuyoshi Onoki
伸好 小野木
Yuzo Imoto
井本  雄三
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Denso Corp
Original Assignee
Denso Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Denso Corp filed Critical Denso Corp
Priority to JP34361296A priority Critical patent/JP3684726B2/en
Publication of JPH10181565A publication Critical patent/JPH10181565A/en
Application granted granted Critical
Publication of JP3684726B2 publication Critical patent/JP3684726B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Regulating Braking Force (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a means to identify the μ on a road surface during the stable running of the wheels of a vehicle without slipping. SOLUTION: A round tire 100 is once deformed in a flat state on ground- contact surface. In this case, a longitudinal force ΔF acting on a tire tread surface varies according to μ, etc., on a road surface. Also, because the material of each part of the tire reduces a resonation during the rotation, the tire is manufactured with the stiffness of these materials varied. Therefore the parts of the tire are deformed differently at the ground-contact surface of the tire. In this way, using as a standard those variation signals (tire rotating pitch) generated when the same position of the tire is brought in contact with the road surface, the state of the road surface is identified by the magnitude of the signals and their continuity. In addition, a tire rotation pitch component is extracted by the cepstrum analysis from wheel speed signals transmitted from a vehicle speed sensor and used to identify the road surface state.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、走行中の路面状態
を判別する路面状態判別装置に関し、ドライバーへの危
険情報の通知を目的とするシステムに関するもの、アン
チスキッドブレーキ装置、トラクションコントロール装
置等の走行系の制御へ路面情報を伝達することにより車
両走行の安全性を向上するシステムに適用することもで
きる。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a road condition judging device for judging a road condition during traveling, and relates to a system for notifying a driver of danger information, an anti-skid brake device, a traction control device and the like. The present invention can also be applied to a system that improves the safety of vehicle traveling by transmitting road surface information to the control of the traveling system.

【0002】[0002]

【従来の技術】従来、路面における摩擦係数μ(以下、
μという)の判別方法として、車両の減速・加速状態に
おけるスリップ率に基づいて判別を行う方法や、旋回中
の横加速度、回転モーメントに基づいて判別を行う方法
が知られている。また、検出された摩擦係数μに基づい
て制御が行われる装置として、アンチスキッドブレーキ
装置やトラクションコントロール装置がある。
2. Description of the Related Art Conventionally, a coefficient of friction μ on a road surface (hereinafter referred to as “μ”)
As a method for determining (μ), there is known a method of making a determination based on a slip ratio in a deceleration / acceleration state of a vehicle, and a method of making a determination based on a lateral acceleration and a turning moment during turning. Further, as a device that performs control based on the detected friction coefficient μ, there are an anti-skid brake device and a traction control device.

【0003】[0003]

【発明が解決しようとする課題】しかしながら、これら
の路面のμの判別方法によると、車両に加速度が発生し
ていない一定速度状態であるとき、或いはステアリング
による操作が行われていないときには、路面のμの判別
を行うことができない。従って、上記したような状態に
おいては、路面のμが低下したという情報等を運転者に
伝達して注意を促したり、アンチスキッドブレーキ装置
等の制御装置に早いタイミングで伝達することができな
かった。
However, according to these methods for determining μ of the road surface, when the vehicle is in a constant speed state where no acceleration is generated or when the steering operation is not performed, μ cannot be determined. Therefore, in the above-described state, it is not possible to transmit information such as a decrease in the μ of the road surface to the driver to warn the driver or to promptly transmit the information to a control device such as an anti-skid brake device. .

【0004】本発明は上記点に鑑みて、車両に加速度が
発生していない一定速度状態であるとき、或いはステア
リングによる操作が行われていないときにおいても路面
のμを判別可能にすることを目的とする。
SUMMARY OF THE INVENTION In view of the foregoing, it is an object of the present invention to make it possible to determine the μ of the road surface even when the vehicle is in a constant speed state in which no acceleration is generated, or when the steering operation is not performed. And

【0005】[0005]

【課題を解決するための手段】円形のタイヤは、接地面
で一度フラットな状態に変形するが、この時タイヤトレ
ッド面に働く前後力は路面のμやタイヤ周上の各部の材
料の強弱によって変化する。図1(a)は接地面でのゴ
ムブロックが前後方向に受ける力を図示したものである
が、転動しているタイヤ100のゴムブロックは最初の
接地部で前方向に引っ張られたのち、上下に圧縮され地
面から離れようとするときに後部に引っ張られる。この
力をΔFとする。
Means for Solving the Problem A circular tire is once deformed into a flat state on a tread, but at this time, the longitudinal force acting on the tire tread surface depends on the μ of the road surface and the strength of the material of each part on the tire circumference. Change. FIG. 1A illustrates the force that the rubber block on the ground contact surface receives in the front-rear direction. After the rubber block of the rolling tire 100 is pulled forward at the first ground contact portion, It is compressed up and down and pulled to the rear as it tries to get off the ground. This force is defined as ΔF.

【0006】前後に引っ張られる力ΔFは路面のμ特性
及びタイヤ100のゴムブロックの剛性により異なる。
一般には、路面のμまたはタイヤ剛性が高いほど、力Δ
Fの変動は大きくなる。また、路面のμが高く、剛性が
高いほど或いはタイヤの剛性が高いほど力ΔFの変化率
(勾配)も大きい。ところで、一般にタイヤの各部の材
料は回転中の共振を減らすため、材料の剛性に強弱をつ
けて製作される。このため、接地部にてタイヤ各部が異
なる状態で変形してタイヤ回転毎に変動が起きる。ま
た、タイヤ100のトレッド部及びカーカス部に振動や
強弱の波を形成しながら転動する。このため、路面μが
同様である場合にはタイヤ100の同一箇所で、つまり
タイヤ100の回転と同一又はその整数倍の周波数毎に
同様な波を形成する。
[0006] The force ΔF to be pulled back and forth depends on the μ characteristics of the road surface and the rigidity of the rubber block of the tire 100.
In general, the higher the μ of the road surface or the tire stiffness, the greater the force Δ
The fluctuation of F increases. Also, as the μ of the road surface is higher and the rigidity is higher or the tire rigidity is higher, the change rate (gradient) of the force ΔF is larger. By the way, in general, the material of each part of the tire is manufactured by giving strength to the material in order to reduce resonance during rotation. For this reason, each part of the tire is deformed in a different state at the ground contacting part, and changes every time the tire rotates. Further, the tire 100 rolls while forming vibrations and strong and weak waves on the tread portion and the carcass portion of the tire 100. For this reason, when the road surface μ is the same, a similar wave is formed at the same location of the tire 100, that is, at the same frequency as the rotation of the tire 100 or every integer multiple thereof.

【0007】このため、路面のμの大小及び路面の均一
性は、上記タイヤ100の振動や剛性の強弱によってタ
イヤ100の回転やサスペンション系に発生するタイヤ
100の回転と同期する変動(以下、タイヤ回転ピッチ
という)に基づいて計測することができる。つまり、タ
イヤ100の同一箇所に発生する変動を物差しとしてそ
の大小、継続性を評価することで路面状態を判別するこ
とができる。
For this reason, the magnitude of μ on the road surface and the uniformity of the road surface vary depending on the vibration of the tire 100 and the strength of the stiffness of the tire 100 and fluctuations synchronous with the rotation of the tire 100 and the rotation of the tire 100 generated in the suspension system (hereinafter referred to as “tire”). Rotation pitch). That is, the road surface condition can be determined by evaluating the magnitude and continuity of the fluctuation occurring at the same location of the tire 100 as a ruler.

【0008】具体的には、人工的に形成されたアスファ
ルトのような路面では路面状態が安定して継続する。こ
のため、タイヤ回転ピッチの情報が乱されない。また、
路面μが大きいために力ΔFの変動が大きく、図1
(b)に示すようにタイヤ100の前後歪みが大きくな
り、タイヤ回転ピッチも大きくなる。しかし、自然的に
形成される低μ路面では路面のμが小さく、自然路面で
あるためバラツキがあり、また凹凸もアスファルトより
不均一である。このため、タイヤ回転ピッチの情報が乱
される。また路面μが小さいために力ΔFの変動が小さ
く、図1(c)に示すようにタイヤ100の前後歪みが
小さくなり、タイヤ回転ピッチ成分も小さくなる。
[0008] Specifically, on a road surface such as artificially formed asphalt, the road surface condition stably continues. Therefore, the information of the tire rotation pitch is not disturbed. Also,
Since the road surface μ is large, the fluctuation of the force ΔF is large.
As shown in (b), the front-back distortion of the tire 100 increases, and the tire rotation pitch also increases. However, on a naturally formed low μ road surface, μ of the road surface is small, the road surface is uneven because of the natural road surface, and unevenness is more uneven than asphalt. For this reason, information on the tire rotation pitch is disturbed. Further, since the road surface μ is small, the fluctuation of the force ΔF is small, and as shown in FIG. 1C, the front-back distortion of the tire 100 is reduced, and the tire rotation pitch component is also reduced.

【0009】上記理由により車両の走行中、タイヤ回転
ピッチを計測することにより路面のμの状態を判別する
ことができる。即ち、アスファルト路のようにグリップ
があり一定した路面と、雪や氷のようにμが小さくバラ
ツキがある路面をタイヤ回転ピッチを計測することによ
り区別できる。上記目的を達成するために、以下に示す
手段を採用する。
For the above reasons, the state of μ on the road surface can be determined by measuring the tire rotation pitch while the vehicle is running. That is, a road surface having a constant grip such as an asphalt road and a road surface having a small μ and a variation such as snow or ice can be distinguished by measuring the tire rotation pitch. In order to achieve the above object, the following means are adopted.

【0010】請求項1に記載の発明においては、車両の
タイヤ(1〜4)と路面の間の摩擦力に基づく変動であ
ってタイヤの回転と同期した変動成分を含む信号を出力
する手段(11〜14)を備え、この信号から変動成分
を抽出して、この変動成分に基づき路面状態を判別する
ことを特徴とする。なお、変動成分を含む信号として
は、車輪加速度、ホイールトルク、サスペンション歪み
又はホイールサスペンション振動等がある。
According to the first aspect of the present invention, a means for outputting a signal which includes a fluctuation component based on a frictional force between the tires (1 to 4) of a vehicle and a road surface and includes a fluctuation component synchronized with the rotation of the tire. 11 to 14), wherein a fluctuation component is extracted from this signal, and the road surface state is determined based on the fluctuation component. Note that the signal including the fluctuation component includes wheel acceleration, wheel torque, suspension distortion, wheel suspension vibration, and the like.

【0011】上述したように、タイヤ回転ピッチ成分か
ら路面のμを判別することができる。従って、路面判別
をタイヤ回転ピッチで行うことにより車両にスリップが
発生していない一定速度状態であるとき、或いはステア
リングによる操作が行われていないときにおいても、路
面のμを判別することができる。具体的には、請求項2
に示すように、タイヤ回転ピッチの抽出は、ケプストラ
ム解析によって行うことができる。この場合、タイヤ回
転ピッチ成分に基づき、路面のμの判別を行うことがで
きる。
As described above, the μ of the road surface can be determined from the tire rotation pitch component. Therefore, by performing the road surface determination at the tire rotation pitch, the μ of the road surface can be determined even when the vehicle is in a constant speed state in which no slip occurs or when the steering operation is not performed. Specifically, claim 2
As shown in (1), the extraction of the tire rotation pitch can be performed by cepstrum analysis. In this case, it is possible to determine μ of the road surface based on the tire rotation pitch component.

【0012】また、請求項3に記載の発明のように、車
輪速センサ(11〜14)からの信号に基づきタイヤと
路面における摩擦力に基づく変動抽出することができ
る。この場合、車輪速センサからの車輪速度信号から車
輪速度信号成分とタイヤ回転ピッチ成分とを分離して、
分離されたタイヤ回転ピッチ成分から路面のμの判別を
行うことができる。
According to the third aspect of the present invention, it is possible to extract a change based on a frictional force between a tire and a road surface based on a signal from a wheel speed sensor (11 to 14). In this case, the wheel speed signal component and the tire rotation pitch component are separated from the wheel speed signal from the wheel speed sensor,
It is possible to determine μ of the road surface from the separated tire rotation pitch components.

【0013】[0013]

【発明の実施の形態】以下、本発明を図に示す実施形態
について説明する。図2に、本実施形態における構成図
を示す。図2に示すように各車輪1〜4に対応して車輪
速センサ11〜14が取りつけられており、これらの車
輪速センサ11〜14からの信号が路面状態判別用の電
子制御装置(以下、ECUという)5に入力される。ま
た、ECU5には、ステアリング6に取りつけられたス
テアリングセンサ7により検出されるステアリング角も
入力され、これらに基づき路面状態判別が行われる。な
お、ECU5によって判別された路面状態は、ウォーニ
ングランプ8にて運転者に警告される。
BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a block diagram showing a first embodiment of the present invention. FIG. 2 shows a configuration diagram in the present embodiment. As shown in FIG. 2, wheel speed sensors 11 to 14 are attached to the wheels 1 to 4, respectively, and signals from the wheel speed sensors 11 to 14 are used to determine an electronic control unit (hereinafter, referred to as a road surface state determining device). ECU 5). The steering angle detected by the steering sensor 7 attached to the steering 6 is also input to the ECU 5, and the road surface state is determined based on the steering angle. The road condition determined by the ECU 5 is warned by a warning lamp 8 to the driver.

【0014】具体的には、上記路面状態判別は、ケプス
トラム解析を使用してタイヤ回転ピッチ成分を検出する
ことにより行う。ここで、ケプストラム解析とは、音声
認識等で普及している技術であり(参考文献:斉藤 収
三著 音声情報処理の基礎 オーム社)、音声において
ケプストラムの原理に基づき各ピッチ成分を検出してい
る。このケプストラムの原理に基づき、タイヤ回転数ピ
ッチの成分を検出する。
Specifically, the road surface condition is determined by detecting a tire rotation pitch component using cepstrum analysis. Here, cepstrum analysis is a technique that is widely used in speech recognition and the like (Reference: Shozo Saito, Basics of Speech Information Processing), and detects each pitch component in speech based on the cepstrum principle. I have. Based on the principle of the cepstrum, a component of a tire rotation speed pitch is detected.

【0015】即ち、信号s(t)を周期Tのピッチ周期
の信号とみなした場合、数式1のように、信号s(t)
は信号成分s0(t)とピッチ成分Σδの積算値の畳み
込み積分で表せる。
That is, when the signal s (t) is regarded as a signal having a pitch period T, the signal s (t)
Can be expressed by convolution of the integrated value of the signal component s0 (t) and the pitch component Σδ.

【0016】[0016]

【数1】 (Equation 1)

【0017】また、数式1をフーリエ変換(以下、FF
Tという)することにより、s(t)におけるパワース
ペクトル値:S(ω)が信号成分S0(ω)とピッチ成
分の積で表せ、数式2のようになる。
Equation (1) is converted to a Fourier transform (hereinafter referred to as FF).
T), the power spectrum value: S (ω) at s (t) can be represented by the product of the signal component S0 (ω) and the pitch component, and is expressed by Expression 2.

【0018】[0018]

【数2】 (Equation 2)

【0019】さらに、数式2について対数変換を行うと
信号成分とピッチ成分が分離でき、数式3で表される。
Further, when the logarithmic transformation is performed on Expression 2, the signal component and the pitch component can be separated, and are expressed by Expression 3.

【0020】[0020]

【数3】 (Equation 3)

【0021】但し、S(t):入力信号、t:時間、
n:整数、T:ピッチ周期、S0(t):ピッチ以外の
信号成分、*:畳み込み積分 この数式3における第1項はs0(t)のパワースペク
トル値であり、第2項の括弧内は線スペクトルを示す。
ここで、このs0(t)のパワースペクトル値はωに対
して比較的ゆるやかな変化を示し、線スペクトルはωに
対して急激な変化を示す。
Where S (t): input signal, t: time,
n: integer, T: pitch period, S0 (t): signal component other than pitch, *: convolution integral The first term in Equation 3 is the power spectrum value of s0 (t), and the parentheses in the second term are 2 shows a line spectrum.
Here, the power spectrum value of s0 (t) shows a relatively gradual change with respect to ω, and the line spectrum shows a sharp change with ω.

【0022】このため、数式3を変数ωの信号列と見な
したとき、ピッチ成分は信号成分に比べてωに対して急
激な変化を示す。また、一般には、s0(ω)はωに対
して急な変化は少ないため数式3を逆フーリエ変換する
とピッチ成分を高周波(高ケフレンシー)側に、信号成
分を低周波(低ケフレンシー)に分離することができ
る。
For this reason, when Equation 3 is regarded as a signal sequence of the variable ω, the pitch component shows a sharp change with respect to ω as compared with the signal component. Further, in general, s0 (ω) has little abrupt change with respect to ω. Therefore, when inverse Fourier transform is applied to Equation 3, the pitch component is separated into the high frequency (high quefrency) side and the signal component is separated into the low frequency (low quefrency). be able to.

【0023】つまり、S0(ω)は、低周波成分が主体
であるので高周波領域ではあまり波形として現れず、こ
の領域ではタイヤピッチ成分を主体とした波形になるた
め、逆FFTを行うことにより、低周波の車輪速度信号
成分と高周波のタイヤピッチ成分とを分離し、この波形
からタイヤ回転ピッチ成分のみを抽出することができ
る。
That is, since S0 (ω) mainly has a low frequency component, it does not appear as a waveform in a high frequency region, and has a waveform mainly composed of a tire pitch component in this region. The low-frequency wheel speed signal component and the high-frequency tire pitch component can be separated, and only the tire rotation pitch component can be extracted from this waveform.

【0024】従って、車両においてはケプストラムにて
車輪速度信号からタイヤ回転ピッチ成分及び車輪速度信
号成分に分離して、このうちのタイヤ回転ピッチ成分か
ら路面状態判別を行うことができる。また、路面がドラ
イ状態、圧雪状態、氷上状態である場合において、検出
される波形の例を図3(a)〜(c)に示す。図3
(a)に示されるように、ドライ状態においてはタイヤ
回転ピッチ成分はタイヤ1回転周期となる波形として現
れ、タイヤ回転ピッチ成分は大きく現れる。図3(b)
に示されるように、圧雪状態においては凹凸や路面μが
不均一であるため振動がランダムになり無秩序な波形と
なる。図3(c)に示されるように、氷上状態において
はタイヤ1回転周期となる波形として現れるが、タイヤ
回転ピッチ成分は小さくしか現れない。このように、各
状態で特徴ある波形となり、これらの波形に基づきケプ
ストラム解析を行いタイヤ回転ピッチ成分を検出する。
Therefore, in the vehicle, the cepstrum separates the wheel speed signal into the tire rotation pitch component and the wheel speed signal component, and the road surface state can be determined from the tire rotation pitch component. FIGS. 3A to 3C show examples of waveforms detected when the road surface is in a dry state, a snow compaction state, or an on-ice state. FIG.
As shown in (a), in the dry state, the tire rotation pitch component appears as a waveform corresponding to one rotation cycle of the tire, and the tire rotation pitch component appears largely. FIG. 3 (b)
As shown in (1), in the snow compaction state, the unevenness and the road surface μ are non-uniform, so that the vibration is random and the waveform becomes irregular. As shown in FIG. 3 (c), in the on-ice state, it appears as a waveform corresponding to one rotation cycle of the tire, but only a small tire rotation pitch component appears. Thus, characteristic waveforms are obtained in each state, and cepstrum analysis is performed based on these waveforms to detect tire rotation pitch components.

【0025】上述したケプストラム解析を使用した路面
判別検出における各処理を、図4に示す。また、図5に
これらの処理に基づき演算される波形を示す。まず、ス
テップ201では、回転角度センサからローターの歯車
毎(1パルス毎)の車輪速度パルスを検出して、この車
輪速度パルスから車輪速度パルス間隔を求めて入力す
る。なお、この車輪速度パルスの検出には汎用に普及し
ている電磁ピックアップ方式又はフォトカプラ方式等が
利用できる。
FIG. 4 shows each process in the road surface discrimination detection using the above-described cepstrum analysis. FIG. 5 shows waveforms calculated based on these processes. First, in step 201, a wheel speed pulse is detected for each gear (one pulse) of the rotor from the rotation angle sensor, and a wheel speed pulse interval is obtained and input from the wheel speed pulse. For detection of the wheel speed pulse, an electromagnetic pickup system or a photocoupler system widely used widely can be used.

【0026】ステップ202では、サンプルした車輪速
度パルス間隔から、サンプル間隔が一定時間となる車輪
速度Vw(t)を求める。具体的には、図6(a)に示
すように、ローター50の歯車50a毎に発生する車輪
速度パルス毎に、図6(b)に示すような車輪速度Vw
(n)を求め、さらに車輪速度Vw(n)から2次スプ
ライン関数で補間して、図6(c)に示すサンプリング
間隔が一定時間となる車輪速度Vw(t)を求める。
In step 202, from the sampled wheel speed pulse intervals, a wheel speed Vw (t) at which the sample interval is constant is determined. More specifically, as shown in FIG. 6A, for each wheel speed pulse generated for each gear 50a of the rotor 50, a wheel speed Vw as shown in FIG.
(N) is obtained, and further interpolated from the wheel speed Vw (n) by a quadratic spline function to obtain the wheel speed Vw (t) at which the sampling interval shown in FIG.

【0027】続くステップ203では、FFTを実行す
る際の誤差を減らすためにハイパスフィルター(以下、
HPFという)によって速度の低周波成分を除いてい
る。ステップ204〜ステップ208は、ケプストラム
演算における処理である。上記ケプストラムの原理に基
づき、これらのステップについて説明する。ステップ2
04でウィンドウ処理を行い、対象となるウィンドウエ
リアを決める。このエリアが狭いと計算精度が下がり、
広いと応答性が下がるため適切な広さのエリアに設定す
る必要がある。なお、計算範囲をタイヤ回転速度に応じ
て切り替えることにより低速から高速までカバーでき
る。
In the following step 203, a high-pass filter (hereinafter, referred to as "FFT") is used to reduce an error in performing the FFT.
HPF) removes the low frequency component of the velocity. Steps 204 to 208 are processes in the cepstrum calculation. These steps will be described based on the cepstrum principle. Step 2
At 04, window processing is performed to determine a target window area. If this area is small, the calculation accuracy will decrease,
If the area is large, the responsiveness will be reduced. By switching the calculation range according to the tire rotation speed, it is possible to cover from low speed to high speed.

【0028】本実施形態では、このウィンドウエリアを
Ts=2πr/Vb*Nrot(但し、Ts:ウィンド
ウ時間、r:タイヤ半径、Vb:車体速度、Nrot:
タイヤ回転数とする)に基づき決定している。基本的に
は計測エリアをタイヤが何回転するかで決定するもので
あり、Nrotは例えば約20回転ほどで実施する。な
お、図5(a)に示す波形は、検出された車輪速度Vw
(t)を示す波形であり、図5(b)に示す波形は、H
PFの処理後にウィンドウ処理によって決定されたウィ
ンドウエリア内の車輪速度サンプルを拡大した波形であ
る。そして、この図5(b)における波形が、数式1に
おける信号s(t)によって表された波形としており、
これにより信号s(t)が決定する。
In this embodiment, this window area is defined as Ts = 2πr / Vb * Nrot (where Ts: window time, r: tire radius, Vb: vehicle speed, Nrot:
Tire rotation speed). Basically, the measurement area is determined by the number of rotations of the tire, and Nrot is performed at about 20 rotations, for example. The waveform shown in FIG. 5A is the detected wheel speed Vw.
The waveform shown in FIG. 5 (t) and the waveform shown in FIG.
It is the waveform which expanded the wheel speed sample in the window area determined by window processing after PF processing. 5B is a waveform represented by the signal s (t) in Expression 1.
Thereby, the signal s (t) is determined.

【0029】ステップ205では、ウィンドウエリア内
において、上記した数式2に示すようなFFT演算を行
いS(ω)を算出する。ステップ206ではFFT演算
の結果からパワースペクトル値を演算する。続くステッ
プ207では、数式3に示すように、このパワースペク
トル値の対数値logVw(ω)を演算する。なお、こ
の時のlogVw(ω)をプロットすると図5(c)の
波形にて表される。そして、ステップ208では、この
対数値を逆FFT演算する。
In step 205, S (ω) is calculated in the window area by performing an FFT operation as shown in the above equation (2). In step 206, a power spectrum value is calculated from the result of the FFT calculation. In the following step 207, the logarithmic value logVw (ω) of the power spectrum value is calculated as shown in Expression 3. Note that plotting logVw (ω) at this time is represented by a waveform in FIG. In step 208, the logarithmic value is subjected to an inverse FFT operation.

【0030】ステップ209及びステップ210では、
タイヤ回転ピッチ成分と関係した成分のみを取り出すた
めに、ケフレンシー軸上で車体速度と関係したウィンド
ウをかける処理を行う。まず、ステップ209では、H
PFによってノイズ等で発生したオフセット成分を取り
除く。このとき、オフセット成分を取り除いた後の波形
は図5(d)の波形にて表される。
In steps 209 and 210,
In order to extract only the component related to the tire rotation pitch component, a process of applying a window related to the vehicle speed on the quefrency axis is performed. First, in step 209, H
The PF removes offset components generated by noise or the like. At this time, the waveform after removing the offset component is represented by the waveform in FIG.

【0031】そして、ステップ210では、車輪速度信
号成分を取り除きタイヤの回転ピッチのみを検出するた
めのウィンドウ処理を行う(図5(e)の波形に該当す
る)。ウィンドウはタイヤの回転周期にあわせて一定の
間隔でn個並べられ、その配置間隔ΔTはΔT=2πr
/Vbによって算出される。また、nはケプストラムグ
ラフ上に発生するピッチの次数成分で決められるが、通
常のタイヤでは4次程であるため、本実施形態において
はn=4で処理を行っている。また、ウィンドウにはハ
ミングウィンドウ:Hw(i,q)を使用しており、数
式4のような関数として表すことができる。
In step 210, window processing is performed to remove the wheel speed signal component and detect only the rotation pitch of the tire (corresponding to the waveform in FIG. 5 (e)). N windows are arranged at regular intervals in accordance with the rotation cycle of the tire, and the arrangement interval ΔT is ΔT = 2πr
/ Vb. Further, n is determined by the order component of the pitch generated on the cepstrum graph, but is about the fourth order in a normal tire, and therefore, in this embodiment, the processing is performed with n = 4. Further, a hamming window: Hw (i, q) is used for the window, and can be expressed as a function as shown in Expression 4.

【0032】また、一般に高μ路面の時には高次成分が
でやすいため、この場合には数式5に示すように、ウィ
ンドウゲイン:G(i)を高く設定する。なお、所定の
定数はハミングウィンドウに用いられる所定の定数であ
り、qはケフレンシー、Mはウィンドウ幅である。
In general, high-order components are likely to occur on a high μ road surface. In this case, the window gain: G (i) is set high as shown in Expression 5. Here, the predetermined constant is a predetermined constant used for the Hamming window, q is quefrency, and M is the window width.

【0033】[0033]

【数4】 (Equation 4)

【0034】[0034]

【数5】 (Equation 5)

【0035】また、上記数式4、5に基づきウィンドウ
処理後のケプストラム値:Cw(q)は数式6により求
められる(図5(f)の波形に該当する)。なお、C
(q)はオフセットをHPFで取り除いた後のケプスト
ラム値である。
The cepstrum value after window processing: Cw (q) is obtained by Expression 6 based on Expressions 4 and 5 above (corresponding to the waveform in FIG. 5F). Note that C
(Q) is a cepstrum value after the offset is removed by the HPF.

【0036】[0036]

【数6】 (Equation 6)

【0037】そして、ステップ211では、ウィンドウ
領域における2乗和を算出して、タイヤ回転ピッチの各
次数成分の総和(タイヤ回転ピッチパラメータ):Cp
を求める。
In step 211, the sum of squares in the window area is calculated, and the sum of the respective components of the tire rotation pitch (tire rotation pitch parameter): Cp
Ask for.

【0038】[0038]

【数7】 (Equation 7)

【0039】なお、実験的に設定された各種路面状態に
ついて、図7(a)に示すような速度で走行した場合に
おける各々のケプストラム解析結果Cpを図7(b)に
示す。このように、路面の種類に応じてケプストラム解
析の結果が異なる。このため、ケプストラム解析によっ
て路面状態の判別を行うことができる。ステップ212
では、先程求めたタイヤ回転ピッチ成分と車体加速度及
びステアリング角をもとに路面判別を実施する。つま
り、タイヤ回転ピッチ成分はステアリング角や車体加速
度によって変化するため、より完全に路面状態を判別す
るためである。
FIG. 7B shows the cepstrum analysis results Cp of various road surface conditions experimentally set when the vehicle travels at a speed as shown in FIG. 7A. As described above, the result of the cepstrum analysis differs depending on the type of the road surface. For this reason, the road surface state can be determined by the cepstrum analysis. Step 212
Then, the road surface is determined based on the tire rotation pitch component, the vehicle body acceleration, and the steering angle obtained earlier. That is, since the tire rotation pitch component changes according to the steering angle and the vehicle body acceleration, the road surface state is determined more completely.

【0040】具体的には、この判別は路面の状態に基づ
き設定された図8に示す路面状態選択マップに基づいて
行い、タイヤ回転ピッチパラメータCp、ステアリング
角及び車体加速度により表される座標A(タイヤ回転ピ
ッチパターン,ステアリング角,車体加速度)が図8に
示される座標軸上のどこに位置するかによって路面判別
を行う。なお、車体加速度は、4輪の車輪速度に基づい
て算出している。
More specifically, this determination is made based on the road surface condition selection map shown in FIG. 8 set based on the road surface condition, and the coordinates A (represented by the tire rotation pitch parameter Cp, the steering angle and the vehicle body acceleration). The road surface is determined based on where the tire rotation pitch pattern, steering angle, and vehicle body acceleration) are located on the coordinate axes shown in FIG. The vehicle body acceleration is calculated based on the wheel speeds of the four wheels.

【0041】具体的には、図8に示されるようにアスフ
ァルトやコンクリート等の高μ路と雪道や氷上等の低μ
路とを区別するしきい面Bがあり、座標Aがこのしきい
面Bの上側にあれば高μ路と判別され、逆に下側にあれ
ば低μ路と判別される。このとき、運転者がステアリン
グを操作していない場合にはステアリング角は零であ
り、また車速が一定である場合には車体加速度は零であ
る。このような場合には、座標Aは(Cp,0,0)と
なる。
More specifically, as shown in FIG. 8, a high μ road such as asphalt or concrete and a low μ road such as a snowy road or on ice are used.
There is a threshold surface B for distinguishing the road from the road. If the coordinates A are above the threshold surface B, the road is determined to be a high μ road, and if the coordinates A are below, the road is determined to be a low μ road. At this time, when the driver does not operate the steering, the steering angle is zero, and when the vehicle speed is constant, the vehicle body acceleration is zero. In such a case, the coordinates A are (Cp, 0, 0).

【0042】ステップ213では、判別結果に基づき運
転者への通知を行う。通知には、路面μ低下を示すウォ
ーニングランプを使用し、例えば高μ路から低μ路に移
行した場合に数秒間ウォーニングランプが点灯させる等
して通知する。また、アンチロックブレーキ装置やトラ
クションコントロール装置等に対しても予め低μ路であ
る旨の信号を送ることにより、これらの制御初期におけ
るスリップを有効に抑えることができる。
In step 213, the driver is notified based on the determination result. For the notification, a warning lamp indicating a decrease in the road surface μ is used. For example, when the road changes from a high μ road to a low μ road, the warning lamp is turned on for a few seconds to give a notification. In addition, by sending a signal indicating that the road is a low μ road in advance to an antilock brake device, a traction control device, and the like, it is possible to effectively suppress the slip in the initial stage of these controls.

【0043】このように、車両に加速度が発生していな
い一定速度状態であるとき、或いはステアリングによる
操作が行われていないときにおいても路面の摩擦係数μ
を検出することができる。 (第2の実施形態)本実施形態においては、第1実施形
態と異なるサンプル間隔でケプストラム解析を行う。
As described above, even when the vehicle is in a constant speed state in which no acceleration is generated, or when the steering operation is not performed, the friction coefficient μ of the road surface is not changed.
Can be detected. (Second Embodiment) In this embodiment, cepstrum analysis is performed at a sample interval different from that in the first embodiment.

【0044】ケプストラム解析は音声処理に使用する都
合上通常は時間軸上で解析を行うため、第1実施形態に
おいてはこの時間軸上で上記解析を行ったが、本実施形
態においては距離軸上で解析を行う。この場合、タイヤ
回転ピッチは、タイヤの回転する距離に比例してピッチ
が発生するため、距離軸上で解析を行った方が直接的な
解析を行うことができる。
Since the cepstrum analysis is usually performed on the time axis because of the use for speech processing, the above analysis was performed on this time axis in the first embodiment. Perform analysis with. In this case, since the pitch of the tire rotation pitch is generated in proportion to the rotation distance of the tire, a direct analysis can be performed by performing the analysis on the distance axis.

【0045】具体的には、図5(b)に示すように、第
1実施形態においては車輪速度におけるサンプリング間
隔を一定時間の周期毎にしているため、車輪速度Vw
(t)をローター50の歯車50a毎に取り出された車
輪速度Vw(n)から補間する必要がある。しかし、図
5(c)に示すように、車輪速度におけるサンプリング
間隔を一定距離の周期毎にすると、距離はローター50
の歯車50aの移動に比例しているため、サンプリング
間隔をローター50の歯車50a毎(1パルス毎)にす
ることができ、上記補間の必要がなくなる。
More specifically, as shown in FIG. 5 (b), in the first embodiment, the sampling interval of the wheel speed is set to a predetermined time period, so that the wheel speed Vw
It is necessary to interpolate (t) from the wheel speed Vw (n) extracted for each gear 50a of the rotor 50. However, as shown in FIG. 5 (c), if the sampling interval at the wheel speed is set to be a cycle of a fixed distance, the distance becomes equal to the rotor 50
Is proportional to the movement of the gear 50a, so that the sampling interval can be set for each gear 50a of the rotor 50 (every pulse), eliminating the need for the interpolation.

【0046】このように距離に比例した解析を行うこと
によりケプストラム値もタイヤ回転に同期した距離の位
置に並ぶ。ケプストラム解析を使用した路面判別検出に
おける各処理を図9に示す。以下、図9に基づきこれら
の処理を説明する。ステップ301にて、1個分のパル
スが伝達されるのに所要時間を1周期として、この周期
毎の車輪速度に対応した車輪速度パルスを検出する。次
に、ステップ302では、車輪速度パルスに基づき、V
=k×Δn/ΔTの関係から車輪速度Vw(n)を演算
する。
By performing the analysis in proportion to the distance in this manner, the cepstrum values are arranged at the position of the distance synchronized with the rotation of the tire. FIG. 9 shows each process in road surface discrimination detection using cepstrum analysis. Hereinafter, these processes will be described with reference to FIG. In step 301, a period required for transmitting one pulse is defined as one cycle, and a wheel speed pulse corresponding to the wheel speed in each cycle is detected. Next, in step 302, based on the wheel speed pulse, V
= K × Δn / ΔT to calculate the wheel speed Vw (n).

【0047】また、ステップ303〜ステップ309に
ついては、第1実施形態におけるステップ203〜ステ
ップ209と同様の処理を行う(図10(a)の波形に
該当する)。そして、ステップ310では、タイヤ回転
ピッチを検出するためのウィンドウ処理を行う(図10
(b)の波形に該当する)。但し、この際には算出され
たケフレンシーの単位がパルス数となるため、タイヤ1
回転の検出パルス数(歯数)をNsとすると、このパル
ス数Nsの定数倍でウィンドウ位置を決定する。
Steps 303 to 309 perform the same processing as steps 203 to 209 in the first embodiment (corresponding to the waveform in FIG. 10A). Then, in step 310, window processing for detecting the tire rotation pitch is performed (FIG. 10).
(This corresponds to the waveform of (b).) However, in this case, the unit of the calculated quefrency is the number of pulses.
Assuming that the number of rotation detection pulses (number of teeth) is Ns, the window position is determined by a constant multiple of the number of pulses Ns.

【0048】この後、ステップ311〜ステップ313
について、第1実施形態におけるステップ211〜ステ
ップ213と同様の処理を行い処理を終了する。なお、
このように距離軸上でケプストラム解析を行っているた
め回転と同期していないノイズを除くという効果が得ら
れる。つまり、距離軸上でケプストラム解析を行う場
合、サンプリング時に取り出される信号は、回転と同期
した成分を持っている。このため、その信号における最
大周波数は、回転周期に比例する周波数となる。従っ
て、サンプル周波数を回転に比例した値とすれば、サン
プリング時に取り出される信号と一定の比例関係を保つ
ことができ、車速によらずエリアスが発生しにくくな
る。
Thereafter, steps 311 to 313 are executed.
, The same processing as in steps 211 to 213 in the first embodiment is performed, and the processing ends. In addition,
As described above, since the cepstrum analysis is performed on the distance axis, an effect of removing noise that is not synchronized with rotation can be obtained. That is, when cepstrum analysis is performed on the distance axis, a signal extracted at the time of sampling has a component synchronized with rotation. Therefore, the maximum frequency of the signal is a frequency proportional to the rotation cycle. Therefore, if the sample frequency is set to a value proportional to the rotation, a constant proportional relationship can be maintained with the signal taken out at the time of sampling, and aliasing is less likely to occur regardless of the vehicle speed.

【0049】(他の実施形態)第1、第2実施形態にお
いては、タイヤ回転ピッチ成分の検出を、車輪速センサ
からの信号(車輪速度信号)に基づき行ったが、タイヤ
と路面の摩擦力に基づく変動を信号出力できるもの、例
えば上下加速度センサやホイールトルクセンサ、サスペ
ンション歪み検出センサ又はホイールサスペンション振
動検出センサ等からの信号に基づき行うこともできる。
(Other Embodiments) In the first and second embodiments, the detection of the tire rotation pitch component is performed based on the signal (wheel speed signal) from the wheel speed sensor. Can be output as a signal, for example, a signal from a vertical acceleration sensor, a wheel torque sensor, a suspension distortion detection sensor, a wheel suspension vibration detection sensor, or the like.

【0050】また、ケプストラム解析によって、タイヤ
回転ピッチ成分を求めたが、LPCケプストラムの係数
によって求めることもできる。また、スタッドレスタイ
ヤ、サマータイヤ等タイヤの種類に応じてもタイヤ回転
ピッチが変化するため、タイヤの種類別に路面状態の判
別を行うこともできる。例えば、タイヤの種類別に図6
に示すような路面状態判別用マップを設定しておき、車
両走行時には予めこのタイヤの種類を示す記号などを入
力することで、タイヤの種類に応じた路面状態判別用マ
ップに補正すればよい。
Although the rotation pitch component of the tire is obtained by the cepstrum analysis, it can be obtained by the coefficient of the LPC cepstrum. In addition, since the tire rotation pitch changes depending on the type of tire such as a studless tire and a summer tire, the road surface condition can be determined for each type of tire. For example, FIG.
A road surface condition determination map as shown in FIG. 1 is set in advance, and a symbol indicating the type of tire is input in advance when the vehicle is running, so that the road surface condition determination map according to the tire type may be corrected.

【0051】さらに、これと同様にアスファルト等の路
面が均一な高μ路を車両が走行中に、状態設定用のスイ
ッチ等を押すことにより、初期状態として入力してその
時のタイヤ回転ピッチパラメータCpを保存し、この値
を基準値として路面状態判別を行うこともできる。
Similarly, when the vehicle is running on a high μ road such as asphalt or the like having a uniform road surface, by pressing a state setting switch or the like, the state is inputted as an initial state and the tire rotation pitch parameter Cp at that time is input. Can be stored, and the road surface state can be determined using this value as a reference value.

【図面の簡単な説明】[Brief description of the drawings]

【図1】(a)はタイヤが路面からうける力を示す模式
図であり、(b)及び(c)は路面状態に対するタイヤ
回転ピッチの大きさを示す模式図である。
FIG. 1A is a schematic diagram showing a force that a tire receives from a road surface, and FIGS. 1B and 1C are schematic diagrams showing a magnitude of a tire rotation pitch with respect to a road surface state.

【図2】車両各センサにおける回路的な構成を示す模式
図である。
FIG. 2 is a schematic diagram showing a circuit configuration of each sensor of the vehicle.

【図3】路面状態に応じてサンプルとして検出される波
形を示した説明図である。
FIG. 3 is an explanatory diagram showing a waveform detected as a sample according to a road surface condition.

【図4】第1実施形態における路面状態判別処理のフロ
ーチャートである。
FIG. 4 is a flowchart of a road surface state determination process in the first embodiment.

【図5】図4における処理に対応した波形を示す特性図
である。
FIG. 5 is a characteristic diagram showing a waveform corresponding to the processing in FIG.

【図6】サンプル間隔として、時間又は距離を採用した
場合における比較図である。
FIG. 6 is a comparison diagram in a case where time or distance is adopted as a sample interval.

【図7】各種の路面状態に応じたケプストラム解析の結
果を示す説明図である。
FIG. 7 is an explanatory diagram showing results of cepstrum analysis according to various road surface conditions.

【図8】路面状態判別を行う基準とする選択マップを示
す説明図である。
FIG. 8 is an explanatory diagram showing a selection map used as a reference for performing road surface state determination.

【図9】第2実施形態における路面状態判別処理のフロ
ーチャートである。
FIG. 9 is a flowchart of a road surface state determination process according to the second embodiment.

【図10】図9における処理に対応した波形を示す特性
図である。
FIG. 10 is a characteristic diagram showing a waveform corresponding to the processing in FIG. 9;

【符号の説明】[Explanation of symbols]

1〜4…車輪、5…電子制御装置、6…ステアリングセ
ンサ、7…ウォーニングランプ、11〜14…車輪速セ
ンサ、50…ローター。
1-4: wheels, 5: electronic control unit, 6: steering sensor, 7: warning lamp, 11-14: wheel speed sensor, 50: rotor.

Claims (4)

【特許請求の範囲】[Claims] 【請求項1】 車両のタイヤ(1〜4)と路面における
摩擦力に基づく変動であって前記タイヤの回転と同期し
た変動成分を含む信号を出力する手段(11〜14)を
備え、 前記信号から前記変動成分を抽出して、この変動成分に
基づき路面状態を判別することを特徴とする路面状態判
別装置。
1. A means (11-14) for outputting a signal which includes a fluctuation component synchronized with the rotation of the tire (1-4) and the rotation of the tire based on a frictional force on a road surface, the signal comprising: A road surface state determination device that extracts the fluctuation component from the road surface and determines a road surface state based on the fluctuation component.
【請求項2】 前記信号をケプストラム解析して前記変
動成分の抽出を行うことを特徴とする請求項1に記載の
路面状態判別装置。
2. The road surface state determination device according to claim 1, wherein the signal is subjected to cepstrum analysis to extract the fluctuation component.
【請求項3】 前記信号を出力する手段は、車輪速セン
サ(11〜14)であることを特徴とする請求項1又は
2に記載の路面状態判別装置。
3. The road surface condition determination device according to claim 1, wherein the means for outputting the signal is a wheel speed sensor (11 to 14).
【請求項4】 車両のタイヤ(1〜4)と、 前記タイヤの回転に同期して発生する変動成分を含む信
号を出力する手段(11〜14)とを備え、 前記変動成分に基づき路面状態を判別することを特徴と
する路面状態判別装置。
4. A tire (1-4) for a vehicle, and means (11-14) for outputting a signal containing a fluctuation component generated in synchronization with the rotation of the tire, and a road surface condition based on the fluctuation component. A road surface condition determination device characterized in that the determination is made.
JP34361296A 1996-12-24 1996-12-24 Road surface condition discrimination device Expired - Fee Related JP3684726B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP34361296A JP3684726B2 (en) 1996-12-24 1996-12-24 Road surface condition discrimination device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP34361296A JP3684726B2 (en) 1996-12-24 1996-12-24 Road surface condition discrimination device

Publications (2)

Publication Number Publication Date
JPH10181565A true JPH10181565A (en) 1998-07-07
JP3684726B2 JP3684726B2 (en) 2005-08-17

Family

ID=18362888

Family Applications (1)

Application Number Title Priority Date Filing Date
JP34361296A Expired - Fee Related JP3684726B2 (en) 1996-12-24 1996-12-24 Road surface condition discrimination device

Country Status (1)

Country Link
JP (1) JP3684726B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000108863A (en) * 1998-10-02 2000-04-18 Toyota Central Res & Dev Lab Inc Physical quantity estimating device and abs controller
JP2005145338A (en) * 2003-11-18 2005-06-09 Mitsubishi Motors Corp Anti-lock brake device
JP2010264917A (en) * 2009-05-15 2010-11-25 Bridgestone Corp Method and device for estimating road surface state and method for controlling vehicle
CN113795415A (en) * 2019-05-07 2021-12-14 普利司通欧洲有限公司 Method and system for identifying unevenness of road surface
CN113825686A (en) * 2019-05-07 2021-12-21 普利司通欧洲有限公司 Method and system for identifying unevenness of road surface

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000108863A (en) * 1998-10-02 2000-04-18 Toyota Central Res & Dev Lab Inc Physical quantity estimating device and abs controller
JP2005145338A (en) * 2003-11-18 2005-06-09 Mitsubishi Motors Corp Anti-lock brake device
JP2010264917A (en) * 2009-05-15 2010-11-25 Bridgestone Corp Method and device for estimating road surface state and method for controlling vehicle
CN113795415A (en) * 2019-05-07 2021-12-14 普利司通欧洲有限公司 Method and system for identifying unevenness of road surface
CN113825686A (en) * 2019-05-07 2021-12-21 普利司通欧洲有限公司 Method and system for identifying unevenness of road surface

Also Published As

Publication number Publication date
JP3684726B2 (en) 2005-08-17

Similar Documents

Publication Publication Date Title
JP2836652B2 (en) Tire pressure detector
JP5837341B2 (en) Road surface condition determination method and apparatus
US6668637B2 (en) Tire air pressure estimating apparatus
JP5165603B2 (en) Tire running state estimation method, steady running state estimation device, tire wear estimation method and apparatus
JP5557569B2 (en) Road surface condition estimation method
EP1798077A2 (en) Apparatus, method and program for alarming decrease in tire air-pressure
JP2007055284A (en) Road surface state estimation method, tire for road surface state estimation, road surface state estimation device and vehicle control device
WO2006054744A1 (en) Tire slipping state detecting method and tire slipping state detecting device
EP1925516B1 (en) Brake control method and brake control device
JP4242047B2 (en) Anti-collision control device
JP2002160511A (en) Air pressure of tire estimating device
JP3684726B2 (en) Road surface condition discrimination device
JPH06138018A (en) Road surface condition measuring device
JPH06258196A (en) Detector of road face friction coefficient
JP3362522B2 (en) Tire pressure detector
JP2008039769A (en) Apparatus for evaluating traveling safety of vehicle, and method for evaluating traveling safety of vehicle
JP3136801B2 (en) Tire pressure detector
JP3391486B2 (en) Tire pressure detector
JP4027585B2 (en) Physical quantity estimation device and road surface friction state estimation device
JP3289384B2 (en) Tire pressure detector
JP3438599B2 (en) Tire pressure estimation device
JPH06328920A (en) Tire pneumatic pressure detector
JP3095095B2 (en) Tire abnormal wear detection device
JP3136772B2 (en) Tire pressure detector
JP3289312B2 (en) Tire pressure detector

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20041027

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20041102

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20041228

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: 20050510

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050523

R150 Certificate of patent (=grant) or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080610

Year of fee payment: 3

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110610

Year of fee payment: 6

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110610

Year of fee payment: 6

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120610

Year of fee payment: 7

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120610

Year of fee payment: 7

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130610

Year of fee payment: 8

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140610

Year of fee payment: 9

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees