JP2006224757A - Wheel abnormality detecting device - Google Patents

Wheel abnormality detecting device Download PDF

Info

Publication number
JP2006224757A
JP2006224757A JP2005039215A JP2005039215A JP2006224757A JP 2006224757 A JP2006224757 A JP 2006224757A JP 2005039215 A JP2005039215 A JP 2005039215A JP 2005039215 A JP2005039215 A JP 2005039215A JP 2006224757 A JP2006224757 A JP 2006224757A
Authority
JP
Japan
Prior art keywords
wheel
waveform
frequency
data
synchronization
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.)
Withdrawn
Application number
JP2005039215A
Other languages
Japanese (ja)
Inventor
Katsuhiro Kobayashi
克宏 小林
Hiroyuki Ueda
寛之 上田
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.)
Bridgestone Corp
Original Assignee
Bridgestone 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 Bridgestone Corp filed Critical Bridgestone Corp
Priority to JP2005039215A priority Critical patent/JP2006224757A/en
Publication of JP2006224757A publication Critical patent/JP2006224757A/en
Withdrawn legal-status Critical Current

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To provide a wheel abnormality detecting device capable of detecting abnormality about all wheels even while a vehicle travels at high speed. <P>SOLUTION: The wheel abnormality detecting device 10 comprises a vibration measuring means 1 continuously measuring the vibration strength of each wheel, a wheel speed measuring means 2 continuously measuring the rotation speed of each wheel, an extracting means 3 extracting a waveform synchronizing with rotations of the wheel from data measured by the vibration measuring means and the wheel speed measuring means, and the warning means 4 judging a state of the wheel from the waveform extracted by the extracting means 3 and warning abnormality when the judged state is out of a predetermined normal range, provided in the plurality of wheels equipped to a vehicle during traveling. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

本発明は、走行する車両に装着された車輪の異常を検知する装置に関し、特に、検知に要する処理速度を向上させたものに関する。   The present invention relates to an apparatus for detecting an abnormality of a wheel mounted on a traveling vehicle, and more particularly to an apparatus for improving a processing speed required for detection.

車輪、例えば、これを構成するタイヤの異常を検知する装置として、タイヤ内圧が一定値以下になるとアラームを発してドライバーに警告をする内圧異常警報装置が知られている。この装置では、タイヤのホイールに取り付けられた内圧センサによって内圧を測定し、ある一定の内圧以下になった場合に警告を発するよう構成されている。   As a device for detecting an abnormality of a wheel, for example, a tire constituting the wheel, an internal pressure abnormality alarm device is known that issues an alarm and warns a driver when the tire internal pressure becomes a predetermined value or less. In this apparatus, an internal pressure is measured by an internal pressure sensor attached to a wheel of a tire, and a warning is issued when the pressure falls below a certain internal pressure.

しかしながら、この内圧異常警報装置では、内圧低下以外の原因で発生する故障は検知できない。例えば、そのような故障の一例として、トレッドとベルト間、ベルトを構成するコード間、及び、サイドゴムとカーカスプライ間などでの剥離、プライコードやベルトコードの破断、ならびに、トレッドゴムのチャンクアウト(例えば、トレッドに設けられたブロック陸部がもぎ取られた状態)などがある。これらのタイヤ故障が生じた状態で走行を続ければ、突然のタイヤバーストが生じて走行不能となる可能性があり、加えて、大きな事故を招く恐れがある。   However, this internal pressure abnormality alarm device cannot detect a failure caused by a cause other than a decrease in internal pressure. For example, as an example of such a failure, separation between the tread and the belt, between the cords constituting the belt, between the side rubber and the carcass ply, breakage of the ply cord or the belt cord, and chunk-out of the tread rubber ( For example, there is a state in which the block land portion provided in the tread is torn off. If the vehicle continues to run in a state where these tire failures occur, a sudden tire burst may occur and the vehicle may become unable to run, and in addition, a serious accident may occur.

このような内圧以外の故障を検知するシステムとして、タイヤの振動や音のデータを計測し、予め求めた正常時のデータと比較することで、タイヤの異常等、回転体の異常を検知システムが知られていて(例えば、特許文献1参照)、この異常検知システムでも、タイヤの異常を検知することはできるが、より簡単な構成で、さらに、性能の良い車両の異常検知システムが近年求められていた。
特開2003−80912号公報
As a system for detecting failures other than such internal pressure, the system for detecting abnormalities in the rotating body such as tire abnormalities is measured by measuring tire vibration and sound data and comparing them with normal data obtained in advance. Although known (see, for example, Patent Document 1), even with this abnormality detection system, it is possible to detect a tire abnormality. However, in recent years, there is a demand for a vehicle abnormality detection system with a simpler configuration and higher performance. It was.
JP 2003-80912 A

このような背景のもと、本出願人は、例えば車輪の振動を連続的に計測し、この振動データから異常状態を検出するに際し、適応デジタルフィルタを用いて、車輪の回転に同期する波形だけを取り出し、その波形を正常なものと比較することにより、車輪の異常を検知する装置を提案した。   Under such a background, the present applicant, for example, continuously measures the vibration of the wheel, and uses an adaptive digital filter to detect only the waveform synchronized with the rotation of the wheel when detecting an abnormal state from the vibration data. We have proposed a device that detects wheel abnormalities by taking out the motor and comparing the waveforms with normal ones.

しかしながら、この提案は、原理的には、高精度に車輪の異常を検出することができるものの、高速で走行する車両において、短時間で車両に装着された複数の車輪の異常検出を、検出に十分な頻繁で行うためには、更なる工夫を必要とするものであった。   However, although this proposal can detect the abnormality of the wheels with high accuracy in principle, the detection of the abnormality of a plurality of wheels attached to the vehicle in a short time in a vehicle traveling at a high speed can be detected. In order to do it frequently enough, further contrivances were required.

本発明は、このような問題点に鑑みてなされたものであり、複数の車輪を装着する車両において、この車両が高速走行をしている状態であっても、これらの車輪のすべてについて、頻繁に、その異常をチェックし、異常があればそれを検出することのできる車輪異常検知装置を提供することを目的とする。   The present invention has been made in view of such problems, and in a vehicle in which a plurality of wheels are mounted, even if the vehicle is traveling at a high speed, the wheel is frequently Another object of the present invention is to provide a wheel abnormality detection device capable of checking the abnormality and detecting the abnormality.

<1>は、走行中の車両に装着された複数の車輪に対して、各車輪の振動強度を連続的に計測する振動計測手段と、各車輪の回転速度を、直接的に、もしくは、車輪速と相関する車両特性を測定することによって間接的に、連続的して計測する車輪速計測手段と、前記振動計測手段および車輪速計測手段で計測されたデータから車輪の回転に同期した波形を抽出する抽出手段と、抽出手段で抽出された波形から車輪の状態を判定しこの状態が予め定めた正常範囲から外れた場合に異常を警告する異常警告手段とを具えた車輪異常検知装置において、
前記抽出手段は、
振動計測手段によって連続的に計測された各車輪についての振動強度データを、予め定められた固定の時間間隔δfでサンプリングして、時間順に並んだデータ列よりなる「計測波形」を生成する固定サンプリング処理部と、
各車輪について、前記計測波形の中から、その車輪の回転に同期した成分だけを、時間順に並んだデータ列として抽出し、これを、「同期成分波形」として前記異常警告手段に出力する適応デジタルフィルタ部と、
各車輪の前記抽出のモデルとなる「モデル波形」を、時間順に並んだデータ列として生成するモデル波形生成部とよりなり、
前記適応デジタルフィルタ部は、前記モデル波形を入力し、前記計測波形と前記同期成分波形との差分、および、入力したモデル波形の両方に基づいて前記同期成分波形を生成するとともに、前記同期成分波形を生成する際のフィルタリングパラメータを、前記差分を用いて適応制御するよう構成され、
各車輪の一回転分の周期に対応する基本周波数に自然数nを乗じた周波数を、その車輪に対する「n次の同期周波数」としたとき、
各車輪に対する前記モデル波形は、このモデル波形の各データを前記δfの時間間隔で並べた時間領域波形をフーリエ変換してできる周波数領域波形が、その車輪の、前記抽出に係る次数の同期周波数の位置で山となり、隣接する次数の同期周波数同士の間では谷となる形状をなし、前記山となる部分の頂部は無限小もしくは有限の幅の周波数範囲に点もしくは平坦面として延在し、頂部の高さは、前記抽出に係るすべての次数の同期周波数に対して同じ値となる特性を有するよう構成され、
モデル波形生成部は、前記車輪速計測手段からのデータにもとづいて、各車輪の基本周波数、もしくは、この基本周波数が含まれる範囲を設定し、設定した基本周波数もしくはその範囲から前記周波数領域波形特性を満足するモデル波形を生成するよう構成されてなる車輪異常検知装置である。
<1> is the vibration measuring means for continuously measuring the vibration intensity of each wheel with respect to a plurality of wheels mounted on the running vehicle, and the rotational speed of each wheel directly or A wheel speed measuring means for continuously measuring the vehicle characteristics correlated with the speed, and a waveform synchronized with the rotation of the wheel from the data measured by the vibration measuring means and the wheel speed measuring means. In the wheel abnormality detection device comprising the extracting means for extracting and the abnormality warning means for determining the state of the wheel from the waveform extracted by the extracting means and warning the abnormality when this state is out of the predetermined normal range,
The extraction means includes
Fixed to generate a “measurement waveform” consisting of a data sequence arranged in time order by sampling vibration intensity data for each wheel continuously measured by the vibration measuring means at a predetermined fixed time interval δ f A sampling processing unit;
For each wheel, only the component synchronized with the rotation of the wheel is extracted from the measurement waveform as a data sequence arranged in time order, and this is output as the “synchronous component waveform” to the abnormality warning means. A filter section;
It consists of a model waveform generation unit that generates a “model waveform” as a model of the extraction of each wheel as a data sequence arranged in time order,
The adaptive digital filter unit receives the model waveform, generates the synchronization component waveform based on both the difference between the measurement waveform and the synchronization component waveform, and the input model waveform, and the synchronization component waveform Is configured to adaptively control the filtering parameter in generating the difference using the difference,
When the frequency obtained by multiplying the fundamental frequency corresponding to the period of one rotation of each wheel by a natural number n is the `` nth synchronization frequency '' for that wheel,
The model waveform for each wheel, frequency domain waveform of time domain waveform obtained by arranging the data in the time interval of the [delta] f can be the Fourier transform of this model waveform, the wheel, the next number of the synchronization frequency according to the extracted It becomes a mountain at the position of, and forms a valley between adjacent synchronization frequencies of the order, the top of the portion that becomes the mountain extends to a frequency range of infinitesimal or finite width as a point or flat surface, The height of the top is configured to have the same value for all orders of synchronization frequency for the extraction,
The model waveform generation unit sets a basic frequency of each wheel or a range including the basic frequency based on data from the wheel speed measuring unit, and sets the frequency domain waveform characteristics from the set basic frequency or the range. Is a wheel abnormality detection device configured to generate a model waveform that satisfies

ここで、「各車輪の回転速度を、車輪速と相関する車両特性を測定することによって間接的に、連続的して計測する車輪速計測手段」とは、例えば、車輪速と相関する車両特性である車速を測定することによって各車輪の車輪速の範囲を算出する場合の車速センサを言う。   Here, “wheel speed measuring means for continuously and continuously measuring the rotational speed of each wheel by measuring the vehicle characteristics correlated with the wheel speed” means, for example, the vehicle characteristics correlated with the wheel speed. The vehicle speed sensor in the case of calculating the range of the wheel speed of each wheel by measuring the vehicle speed.

<2>は、<1>において、前記適応デジタルフィルタ部は、複数の車輪に対するそれぞれの同期成分波形を、これらに共通の、一つのモデル波形に基づいて生成するよう構成され、
この共通のモデル波形に対応する前記周波数領域波形において、n次の同期周波数に対応する山の頂部の周波数範囲は、前記共通のモデル波形に基づく同期成分波形の抽出に係る全ての車輪に対するn次の同期周波数を含むことを特徴とする車輪異常検知装置である。
<2> in <1>, wherein the adaptive digital filter unit is configured to generate respective synchronous component waveforms for a plurality of wheels based on a single model waveform common to these,
In the frequency domain waveform corresponding to the common model waveform, the frequency range of the peak corresponding to the n-th synchronization frequency is the n-th order for all wheels related to the extraction of the synchronization component waveform based on the common model waveform. It is a wheel abnormality detection apparatus characterized by including these synchronous frequencies.

<3>は、<2>において、Aを予め定められた百分率として、前記、n次の同期周波数に対応する山の頂部の周波数範囲を、前記共通のモデル波形に基づく同期成分波形の抽出に係る複数の車輪のうち、基本周波数が最小となる車輪のn次同期周波数の0〜+A%として設定してなる車輪異常検知装置である。   <3> is the extraction of the synchronization component waveform based on the common model waveform, with the frequency range at the peak corresponding to the n-th synchronization frequency being defined as a predetermined percentage in <2>. Among the plurality of wheels, the wheel abnormality detection device is set as 0 to + A% of the n-th order synchronization frequency of the wheel having the lowest fundamental frequency.

<4>は、<2>において、前記、n次の同期周波数に対応する山の頂部の周波数範囲を、前記共通のモデル波形に基づく同期成分波形の抽出に係る複数の車輪のうち、基本周波数が最小となる車輪と、基本周波数が最大となる車輪との、両方の車輪に対するn次同期周波数に基づいて設定してなる車輪異常検知装置である。   <4> is the fundamental frequency among the plurality of wheels according to <2>, in which the frequency range of the peak corresponding to the n-th synchronization frequency is extracted from the synchronous component waveform based on the common model waveform Is a wheel abnormality detection device that is set on the basis of the n-th synchronization frequency for both wheels, the wheel having the smallest and the wheel having the largest fundamental frequency.

<5>は、<2>〜<4>のいずれかにおいて、前記前記適応デジタルフィルタ部は、前記共通のモデル波形を構成するデータを、一つ入力するごとに、この共通のモデル波形に基づく同期成分波形の抽出に係る全車輪のそれぞれについて、同期成分波形のデータを一つだけ出力する処理を行うよう構成されてなる車輪異常検知装置である。   <5> is any one of <2> to <4>, wherein the adaptive digital filter unit is based on the common model waveform every time one piece of data constituting the common model waveform is input. The wheel abnormality detection device is configured to perform a process of outputting only one piece of data of the synchronous component waveform for each of all the wheels related to the extraction of the synchronous component waveform.

<6>は、走行中の車両に装着された複数の車輪に対して、各車輪の振動強度を連続的に計測する振動計測手段と、各車輪の回転速度を、直接的に、もしくは、車輪速と相関する車両特性を測定することによって間接的に、連続的して計測する車輪速計測手段と、前記振動計測手段および車輪速計測手段で計測されたデータから車輪の回転に同期した波形を抽出する抽出手段と、抽出手段で抽出された波形から車輪の状態を判定しこの状態が予め定めた正常範囲から外れた場合に異常を警告する異常警告手段とを具えた車輪異常検知装置において、
前記抽出手段は、
各車輪について、振動計測手段により時間軸に沿って連続的に計測された振動強度データを、車輪速計測手段からのデータに応じて変化する時間間隔δvでサンプリングして、時間順に並んだデータ列よりなる「計測波形」を生成する可変サンプリング処理部と、
各車輪について、前記計測波形の中から、その車輪の回転に同期した成分だけを、時間順に並んだデータ列として抽出し、これを、「同期成分波形」として前記異常警告手段に出力する適応デジタルフィルタ部と、
各車輪の前記抽出のモデルとなる「モデル波形」を、時間順に並んだデータ列として生成するモデル波形生成部とよりなり、
前記可変サンプリング処理部において各車輪に対する前記計測波形を生成するのに用いられる前記時間間隔δvは、前記車輪速計測手段からのデータにもとづいて算出されるこの車輪の1回転分の周期に対応する基本周波数、もしくは、前記車輪速計測手段からのデータにもとづいて予め定められた手順で算出される仮想車輪の1回転分の周期に対応する基本周波数faと、これらの車輪の1回転を等分する際の分割数となる予め定められた値kとの積の逆数として求められ、
前記適応デジタルフィルタ部は、前記モデル波形を入力し、前記計測波形と前記同期成分波形との差分、および、入力したモデル波形の両方に基づいて前記同期成分波形を生成するとともに、前記同期成分波形を生成する際のフィルタリングパラメータを、前記差分を用いて適応制御し、複数の車輪に対するそれぞれの同期成分波形は、これらの車輪に共通な一つのモデル波形に基づいて生成されるよう構成され、
予め定められた固定の周波数を固定基本周波数f0として、この固定基本周波数に自然数nを乗じた周波数を、「n次の固定同期周波数」とし、前記固定基本周波数f0と前記分割数kとの積の逆数を時間間隔δ0としたとき、
前記共通のモデル波形は、このモデル波形の各データを前記時間間隔δ0おきに並べた時間領域波形をフーリエ変換してできる周波数領域波形が、前記抽出に係る次数の固定同期周波数の位置で山となり、隣接する次数の固定同期周波数同士の間では谷となる形状をなし、前記山となる部分の頂部は無限小もしくは有限の幅の周波数範囲で点もしくは平坦面として延在し、頂部の高さは、前記抽出に係るすべての次数の固定同期周波数に対して同じ値となる特性を有するよう構成されてなる車輪異常検知装置である。
<6> is the vibration measuring means for continuously measuring the vibration intensity of each wheel with respect to a plurality of wheels mounted on the running vehicle, and the rotational speed of each wheel directly or the wheel. A wheel speed measuring means for continuously measuring the vehicle characteristics correlated with the speed, and a waveform synchronized with the rotation of the wheel from the data measured by the vibration measuring means and the wheel speed measuring means. In the wheel abnormality detection device comprising the extracting means for extracting and the abnormality warning means for determining the state of the wheel from the waveform extracted by the extracting means and warning the abnormality when this state is out of the predetermined normal range,
The extraction means includes
For each wheel, the vibration intensity data continuously measured along the time axis by the vibration measuring means is sampled at a time interval δ v that changes according to the data from the wheel speed measuring means, and the data is arranged in time order. A variable sampling processing unit that generates a "measurement waveform" consisting of columns;
For each wheel, only the component synchronized with the rotation of the wheel is extracted from the measurement waveform as a data sequence arranged in time order, and this is output as the “synchronous component waveform” to the abnormality warning means. A filter section;
It consists of a model waveform generation unit that generates a “model waveform” as a model of the extraction of each wheel as a data sequence arranged in time order,
The time interval δ v used to generate the measurement waveform for each wheel in the variable sampling processing unit corresponds to a cycle of one rotation of the wheel calculated based on data from the wheel speed measuring means. Or a fundamental frequency f a corresponding to a cycle of one rotation of the virtual wheel calculated by a predetermined procedure based on data from the wheel speed measuring means and one rotation of these wheels. It is obtained as the reciprocal of the product with a predetermined value k which is the number of divisions when equally dividing,
The adaptive digital filter unit receives the model waveform, generates the synchronization component waveform based on both the difference between the measurement waveform and the synchronization component waveform, and the input model waveform, and the synchronization component waveform The filtering parameter when generating is adaptively controlled using the difference, and each synchronous component waveform for a plurality of wheels is configured to be generated based on one model waveform common to these wheels,
A predetermined fixed frequency is set as a fixed basic frequency f 0 , and a frequency obtained by multiplying the fixed basic frequency by a natural number n is referred to as an “n-th fixed synchronization frequency”, and the fixed basic frequency f 0 and the division number k are When the reciprocal of the product of the time interval δ 0 ,
The common model waveform is a frequency domain waveform obtained by Fourier transforming a time domain waveform in which each data of the model waveform is arranged at every time interval δ 0 , and is a peak at the position of the fixed synchronization frequency of the order related to the extraction. It forms a valley between the fixed synchronization frequencies of adjacent orders, and the top of the peak portion extends as a point or flat surface in a frequency range of infinitesimal or finite width, and the height of the top The wheel abnormality detection device is configured to have the same value with respect to the fixed synchronization frequencies of all orders related to the extraction.

ここで、「各車輪の回転速度を、車輪速と相関する車両特性を測定することによって間接的に、連続的して計測する車輪速計測手段」とは、例えば、車輪速と相関する車両特性である車速を測定することによって各車輪の車輪速の範囲を算出する場合の車速センサを言う。   Here, “wheel speed measuring means for continuously and continuously measuring the rotational speed of each wheel by measuring the vehicle characteristics correlated with the wheel speed” means, for example, the vehicle characteristics correlated with the wheel speed. The vehicle speed sensor in the case of calculating the range of the wheel speed of each wheel by measuring the vehicle speed.

また、「前記車輪速計測手段からのデータにもとづいて予め定められた手順で算出される仮想車輪の1回転分の周期に対応する基本周波数」における「仮想車輪」とは、例えば、計測波形の生成の対象となる車輪とは異なる車輪を「仮想車輪」とする場合や、スリップ率がゼロの状態で路面を転動しながら車速と同じ速度で進行する所定の直径の車輪を「仮想車輪」とする場合の「仮想車輪」をいう。   In addition, “virtual wheel” in “basic frequency corresponding to a cycle of one rotation of a virtual wheel calculated by a predetermined procedure based on data from the wheel speed measuring means” is, for example, a measurement waveform When a wheel that is different from the wheel to be generated is a “virtual wheel”, or a wheel with a predetermined diameter that travels at the same speed as the vehicle speed while rolling on the road surface with a slip rate of zero is a “virtual wheel” The “virtual wheel” in this case.

<7>は、<6>において、各車輪に対応するサンプリング処理に用いられる前記時間間隔δvは、前記同期成分波形の抽出に係る複数の車輪のすべてに対して共通のものが用いられ、これらの車輪の基本周波数のうち、最小のものをfmin、最大のものをfmax、前記時間間隔δvの算出に用いたものをfaとしたとき、前記モデル波形に対応する周波数領域波形において、n次の固定同期周波数に対応する山の頂部の周波数範囲は、少なくとも、f0・(fmin/fa) のn倍以上でかつf0・(fmax/fa)のn倍以下の範囲を含むよう構成されてなる車輪異常検知装置である。 <7> in <6>, the time interval δ v used for the sampling process corresponding to each wheel is the same for all of the plurality of wheels related to the extraction of the synchronous component waveform, Of the fundamental frequencies of these wheels, the minimum frequency is f min , the maximum frequency is f max , and the frequency used for calculating the time interval δ v is f a, and the frequency domain waveform corresponding to the model waveform , The peak frequency range corresponding to the n-th fixed synchronization frequency is at least n times f 0 · (f min / f a ) and n times f 0 · (f max / f a ) The wheel abnormality detection device is configured to include the following ranges.

<8>は、<7>において、前記共通のモデル波形に基づく同期成分波形の抽出に係る全ての車輪に対する基本周波数のうち、最小のものを、前記、共通の時間間隔δvの算出に用いるとともに、Aを予め定められた百分率として、n次の固定同期周波数に対応する山の頂部の周波数範囲を、n次の固定同期周波数の、0〜+A%の範囲として設定してなる車輪異常検知装置である。 <8>, in <7>, of the basic frequency for all of the wheels according to the extraction of the synchronizing component waveform based on said common model waveform, the smallest ones, the, used to calculate the common time interval [delta] v In addition, the wheel abnormality is defined by setting the frequency range of the peak corresponding to the n-th fixed synchronization frequency as a range of 0 to + A% of the n-th fixed synchronization frequency, where A is a predetermined percentage. It is a detection device.

<9>は、<7>において、n次の固定同期周波数に対応する山の頂部の周波数範囲を、前記共通のモデル波形に基づく同期成分波形の抽出に係る複数の車輪のうち、基本周波数が最小となる車輪と、基本周波数が最大となる車輪との、両方のn次の同期周波数に基づいて設定してなる車輪異常検知装置である。   <9> is the peak frequency range of the plurality of wheels related to extraction of the synchronous component waveform based on the common model waveform in the frequency range of the peak corresponding to the n-th fixed synchronous frequency in <7>. This is a wheel abnormality detection device that is set based on both n-th order synchronization frequencies of the smallest wheel and the wheel having the maximum fundamental frequency.

<10>は、<7>〜<9>のいずれかにおいて、前記適応デジタルフィルタ部は、前記共通のモデル波形を構成するデータを、一つ入力するごとに、この共通のモデル波形に基づく同期成分波形の抽出に係る全車輪のそれぞれについて、同期成分波形のデータを一つだけ出力する処理を行うよう構成されてなる車輪異常検知装置である。   <10> is any one of <7> to <9>, wherein the adaptive digital filter unit synchronizes based on the common model waveform every time one piece of data constituting the common model waveform is input. This is a wheel abnormality detection device configured to perform processing for outputting only one piece of synchronous component waveform data for each of all the wheels related to component waveform extraction.

<1>によれば、振動強度データを所定時間間隔でサンプリングして生成され、時間順に並んだデータ列よりなる「計測波形」と、時間順に並んだデータ列となる「モデル波形」とから、車輪の回転に同期する同期成分波形を抽出するので、フーリエ変換等の時間のかかる処理を行うことなく、所定の時間間隔δfでサンプリングしたものをそのまま、入力順に処理することができ、このことにより、1回の抽出工程を短いものにして、車両が高速走行状態であっても、車輪のすべてについて、異常のチェックを十分な頻度で行うことができる。 According to <1>, vibration intensity data is generated by sampling at predetermined time intervals, and a “measurement waveform” composed of a data sequence arranged in time order and a “model waveform” composed of a data sequence arranged in time order, Since the synchronous component waveform that is synchronized with the rotation of the wheel is extracted, it is possible to process the samples sampled at the predetermined time interval δ f as they are, without performing time-consuming processing such as Fourier transformation. Thus, one extraction process can be shortened, and even if the vehicle is in a high-speed running state, abnormality check can be performed with sufficient frequency for all of the wheels.

<2>によれば、適応デジタルフィルタ部を、複数の車輪に対するそれぞれの同期成分波形が、これらに共通の、一つのモデル波形に基づいて生成されるよう構成したので、モデル波形の生成にかかる時間を低減して、同期成分波形の抽出工程を一層短縮化することができる。   According to <2>, the adaptive digital filter unit is configured so that the respective synchronous component waveforms for the plurality of wheels are generated based on a single model waveform that is common to the plurality of wheels. Time can be reduced and the synchronous component waveform extraction process can be further shortened.

<3>によれば、モデル波形に対応する周波数領域波形における、n次の同期周波数に対応する山の頂部の周波数範囲を、最高速の車輪の基本周波数だけに基づいて設定するので、詳細を後述するように、装置をより安価で簡易なものにすることができる。   According to <3>, the frequency range of the peak corresponding to the n-th synchronization frequency in the frequency domain waveform corresponding to the model waveform is set based on only the fundamental frequency of the fastest wheel. As will be described later, the apparatus can be made cheaper and simpler.

<4>によれば、モデル波形に対応する周波数領域波形における、n次の同期周波数に対応する山の頂部の周波数範囲を、最高速の車輪の基本周波数と、最低速の車輪の基本周波数との両方に基づいて設定するので、詳細を後述するように、車輪に同期する周波数成分の抽出を確実に行わせることができる。   According to <4>, in the frequency domain waveform corresponding to the model waveform, the frequency range of the top of the mountain corresponding to the n-th synchronization frequency is defined as the basic frequency of the fastest wheel and the basic frequency of the slowest wheel. Therefore, the frequency components synchronized with the wheels can be reliably extracted as will be described in detail later.

<5>によれば、適応デジタルフィルタ部は、前記共通のモデル波形を構成するデータを、一つ入力するごとに、この共通のモデル波形に基づく同期成分波形の抽出に係る全車輪のそれぞれについて、同期成分波形のデータを一つだけ出力する処理を行うよう構成されているので、i番目の同期成分波形のデータを生成するのに、モデル波形のメモリもしくはバッファからの呼び出しを、車輪ごとに行うのではなく、全車輪に対して1回だけで済ませることができ、同期成分波形の抽出工程をなお一層短縮化することができる。   According to <5>, the adaptive digital filter unit inputs each piece of data constituting the common model waveform for each of the wheels related to extraction of the synchronous component waveform based on the common model waveform. Since it is configured to perform processing to output only one piece of synchronous component waveform data, in order to generate the i-th synchronous component waveform data, a call from the model waveform memory or buffer is made for each wheel. This can be done only once for all wheels, and the synchronization component waveform extraction process can be further shortened.

<6>によれば、時間順に並んだデータ列よりなる「計測波形」と、時間順に並んだデータ列となる「モデル波形」とから、車輪の回転に同期する同期成分波形を抽出するので、<1>と同様に、1回の抽出工程を短いものにして、車両が高速走行状態であっても、異常のチェックを頻繁に行うことができる。   According to <6>, since a “measurement waveform” composed of a data sequence arranged in time order and a “model waveform” composed of a data sequence arranged in time order, a synchronous component waveform synchronized with the rotation of the wheel is extracted. Similar to <1>, it is possible to frequently check for abnormalities even if the vehicle is in a high-speed traveling state by shortening one extraction step.

<7>によれば、可変サンプリング処理部におけるサンプリングを、これらに共通の、一つの時間間隔に基づいて行われるよう構成したので、サンプリン時間間隔の算出にかかる時間を低減して、同期成分波形の抽出工程を一層短縮化することができる。   According to <7>, since the sampling in the variable sampling processing unit is configured to be performed based on one common time interval, the time taken to calculate the sampling time interval is reduced, and the synchronous component waveform The extraction process can be further shortened.

<8>によれば、モデル波形に対応する周波数領域波形における、n次の同期周波数に対応する山の頂部の周波数範囲を、最高速の車輪の基本周波数だけに基づいて設定するので、<3>と同様に、装置をより安価で簡易なものにすることができる。   According to <8>, since the frequency range of the peak corresponding to the n-th synchronization frequency in the frequency domain waveform corresponding to the model waveform is set based only on the fundamental frequency of the fastest wheel, <3 As with>, the device can be made cheaper and simpler.

<9>によれば、モデル波形に対応する周波数領域波形における、n次の同期周波数に対応する山の頂部の周波数範囲を、最高速の車輪の基本周波数と、最低速の車輪の基本周波数との両方に基づいて設定するので、<4>と同様に、車輪に同期する周波数成分の抽出を確実に実施させることができる。   <9> According to the frequency domain waveform corresponding to the model waveform, the frequency range of the top of the mountain corresponding to the n-th synchronization frequency is expressed by the basic frequency of the fastest wheel and the basic frequency of the slowest wheel. Therefore, it is possible to reliably extract the frequency component synchronized with the wheel as in <4>.

<10>によれば、適応デジタルフィルタ部は、前記共通のモデル波形を構成するデータを、一つ入力するごとに、この共通のモデル波形に基づく同期成分波形の抽出に係る全車輪のそれぞれについて、同期成分波形のデータを一つだけ出力する処理を行うよう構成されているので、<5>と同様に、モデル波形のメモリもしくはバッファからの呼び出しを、車輪ごとに行うのではなく、全車輪に対して1回だけで済ませることができ、同期成分波形の抽出工程をなお一層短縮化することができる。   According to <10>, the adaptive digital filter unit inputs each piece of data constituting the common model waveform for each of the wheels related to extraction of the synchronous component waveform based on the common model waveform. Since it is configured to output only one piece of synchronous component waveform data, as in <5>, the model waveform memory or buffer is not called for each wheel. However, the process of extracting the synchronous component waveform can be further shortened.

本発明の実施形態として、以下に、第一〜第四の実施形態を図に基づいて説明する。これらの説明に用いる主たる記号の意味は下記の通りである。
<記号の説明>
r:車輪番号(r=1〜R、Rは車輪数)
i:時間軸に沿うデータ列におけるデータ順番(i=1〜M、Mはサンプリング数)
p:実施形態の番号(p=1〜4)
ar:第r輪の振動強度データ
vr:第r輪の車輪速強度データ
fr (n):第r輪のn次の同期周波数(fr (1)は基本周波数)
Xpr、Xpr (i):第p実施形態における第r輪の計測波形、および、そのi番目のデータ
Ypr、Ypr (i):第p実施形態における第r輪の同期成分波形、および、そのi番目のデータ
Rpr、Rpr (i):第p実施形態における第r輪のモデル波形、および、そのi番目のデータ
Rp、Rp (i):第p実施形態における各輪に共通なモデル波形、および、そのi番目のデータ
Epr、Epr (i):第p実施形態における第r輪の誤差、および、そのi番目のデータ
δf:固定サンプリングにおける時間間隔
δvr:第r輪の車輪速データに基づく可変サンプリングにおける時間間隔
As embodiments of the present invention, first to fourth embodiments will be described below with reference to the drawings. The meanings of the main symbols used in these explanations are as follows.
<Explanation of symbols>
r: Wheel number (r = 1 to R, R is the number of wheels)
i: Data order in the data sequence along the time axis (i = 1 to M, M is the number of samplings)
p: Embodiment number (p = 1 to 4)
a r : R wheel vibration intensity data
v r : Wheel speed strength data of the r-th wheel
f r (n) : n-th synchronization frequency of the r-th wheel (f r (1) is the fundamental frequency)
X pr , X pr (i): measurement waveform of the r-th wheel and its i-th data in the p-th embodiment
Y pr , Y pr (i): the synchronous component waveform of the r-th wheel in the p-th embodiment and its i-th data
R pr , R pr (i): model waveform of the r-th wheel in the p-th embodiment and its i-th data
R p , R p (i): model waveform common to each wheel in the p-th embodiment and its i-th data
E pr , E pr (i): error of the r-th wheel in the p-th embodiment and its i-th data δ f : time interval in fixed sampling δ vr : in variable sampling based on wheel speed data of the r-th wheel Time interval

図1(a)は、車両における車輪異常検知装置の配置を示す平面模式図であり、図1(b)は、車輪の一部に生じた異常部分を示す車輪の側面模式図であり、また、図2は、車輪異常検知装置の構成を示すブロック線図である。   Fig.1 (a) is a plane schematic diagram which shows arrangement | positioning of the wheel abnormality detection apparatus in a vehicle, FIG.1 (b) is a side surface schematic diagram of the wheel which shows the abnormal part which arose in a part of wheel, FIG. 2 is a block diagram showing the configuration of the wheel abnormality detection device.

車輪異常検知装置10は、走行中の車両に装着された複数の車輪90A、90B、90C、90Dに対して、各車輪の振動強度を連続的に計測する振動計測手段1と、各車輪の回転速度を連続的に計測する車輪速計測手段2と、前記振動計測手段1および車輪速計測手段2で計測されたデータから車輪の回転に同期した波形を抽出する抽出手段3と、抽出手段3で抽出された波形から車輪の状態を判定しこの状態が予め定めた正常範囲から外れた場合に異常を警告する異常警告手段4とを具えて構成され、車両が例えば4輪車の場合、振動計測手段1は、それぞれの車輪90A〜90Dの近く、例えば、車輪を支持するナックルに取り付けられた第1輪〜第4輪の振動センサ1a〜4dで構成され、これらの振動センサ1a〜4aは、それぞれ、連続的に、振動強度データa1〜a4を抽出手段3に出力する。 The wheel abnormality detection apparatus 10 includes a vibration measuring unit 1 that continuously measures the vibration intensity of each wheel with respect to the plurality of wheels 90A, 90B, 90C, and 90D mounted on the traveling vehicle, and the rotation of each wheel. A wheel speed measuring means 2 for continuously measuring the speed; an extracting means 3 for extracting a waveform synchronized with the rotation of the wheel from the data measured by the vibration measuring means 1 and the wheel speed measuring means 2; An abnormality warning means 4 for determining the state of the wheel from the extracted waveform and warning the abnormality when the state is out of a predetermined normal range is provided. When the vehicle is, for example, a four-wheeled vehicle, vibration measurement is performed. Means 1 is composed of first to fourth vibration sensors 1a to 4d attached to the knuckles that support the wheels 90A to 90D, for example, and these vibration sensors 1a to 4a are That , Continuously outputs the vibration intensity data a 1 ~a 4 to the extraction unit 3.

車輪速計測手段2は、車両が4輪車の場合、それぞれの車輪の回転速度を検出する第1輪〜第4輪の車輪速センサ2a〜2dで構成され、これらの車輪速センサ2a〜2dは、それぞれ、車輪速データv1〜v4を抽出手段3に連続的に出力する。振動センサ1a〜4dの例としては、車両の前後方向の加速度を検出する加速度センサを挙げることができ、また、車輪速センサ2a〜2dの例としては、それぞれの車軸に取り付けられABSの制御に用いられるギヤパルサを用いたものを挙げることができる。 When the vehicle is a four-wheeled vehicle, the wheel speed measuring means 2 includes wheel speed sensors 2a to 2d for the first to fourth wheels that detect the rotational speed of each wheel, and these wheel speed sensors 2a to 2d. Respectively output the wheel speed data v 1 to v 4 to the extraction means 3 continuously. Examples of the vibration sensors 1a to 4d include acceleration sensors that detect the acceleration in the longitudinal direction of the vehicle. Examples of the wheel speed sensors 2a to 2d are attached to the respective axles to control the ABS. The thing using the used gear pulsar can be mentioned.

抽出手段3は、振動強度データa1〜a4よりなるそれぞれの波形から、各車輪の回転に同期した波形b1〜b4だけを抽出し、これらを異常警告手段4に出力する。抽出手段3および異常警告手段4は車体側に設けられ、異常警告手段4は異常判定を行う異常制御部4aと、異常の場合に運転者にその旨を警告する異常表示部4bとよりなり、抽出手段3および異常制御部4aは中央制御装置5に配置される。 Extracting means 3, from each of the waveform consisting of vibration intensity data a 1 ~a 4, and extracts only a waveform b 1 ~b 4 synchronized with the rotation of the wheels, and outputs them to the abnormality warning means 4. The extraction unit 3 and the abnormality warning unit 4 are provided on the vehicle body side, and the abnormality warning unit 4 includes an abnormality control unit 4a that performs abnormality determination and an abnormality display unit 4b that warns the driver in the event of an abnormality. The extraction means 3 and the abnormality control unit 4a are arranged in the central control device 5.

中央制御装置5のハードウェア構成は、一般の制御装置に用いられるもので構成することができ、例えば、CPU、メモリ、AD変換を含む入出力インターフェース等で構成される。   The hardware configuration of the central control device 5 can be configured by what is used for a general control device, and includes, for example, a CPU, a memory, an input / output interface including AD conversion, and the like.

図1(b)に示すように、走行中の車両において、いずれかの車輪90Aに何らかの異常、例えば、タイヤのトレッドとベルト間で層間剥離91が起こると、車輪の回転に同期する振動波形が発生する。   As shown in FIG. 1B, in a running vehicle, if any abnormality occurs in any of the wheels 90A, for example, delamination 91 occurs between the tire tread and the belt, a vibration waveform synchronized with the rotation of the wheels is generated. appear.

本発明に係る車輪異常検知装置10は、路面からの振動等が混在する振動強度データa1から、車輪の回転に同期する振動波形だけを抽出して、この波形が所定の範囲を超えた場合、車両になんらかの異常か発生したとして、これを運転者に警告する機能を有するよう構成されている。 The wheel abnormality detection device 10 according to the present invention extracts only the vibration waveform synchronized with the rotation of the wheel from the vibration intensity data a 1 in which vibrations from the road surface and the like are mixed, and the waveform exceeds a predetermined range. The vehicle is configured to have a function of warning the driver of any abnormality that has occurred in the vehicle.

まず、本発明の車両異常検出装置10の第一実施形態について説明する。図3は、この形態の車両異常検出装置の一部を構成する抽出手段3Aを、4輪車に設けられたものを例として示すブロック線図であり、抽出手段3Aは、振動計測手段1によって連続的に計測された各車輪についての振動強度データa1〜a4を、予め定められた固定の時間間隔δfでサンプリングして、時間順に並んだデータ列よりなる計測波形X11〜X14を生成する固定サンプリング処理部11と、各車輪について、計測波形X11〜X14の中から、その車輪の回転に同期した成分だけを、時間順に並んだデータ列として抽出し、これを、同期成分波形b1〜b4(もしくはY11〜Y14)として前記異常警告手段4に出力する適応デジタルフィルタ部12と、各車輪の前記抽出のモデルとなるモデル波形R11〜R14を、時間順に並んだデータ列として生成するモデル波形生成部13とよりなる。 First, a first embodiment of the vehicle abnormality detection device 10 of the present invention will be described. FIG. 3 is a block diagram showing, as an example, the extraction means 3A that constitutes a part of the vehicle abnormality detection device of this embodiment, which is provided in a four-wheeled vehicle. The vibration intensity data a 1 to a 4 for each wheel continuously measured are sampled at a predetermined fixed time interval δ f and measured waveforms X 11 to X 14 consisting of a data sequence arranged in time order. a fixed sampling processing unit 11 for generating, for each wheel, from among the measured waveform X 11 to X 14, only components synchronized with the rotation of the wheels, and extracted as a data string arranged in time order, this synchronization The adaptive digital filter unit 12 that outputs the component waveforms b 1 to b 4 (or Y 11 to Y 14 ) to the abnormality warning means 4 and the model waveforms R 11 to R 14 that serve as the models for the extraction of each wheel, Generated as a sequence of data The more the model waveform generating unit 13 that.

固定サンプリング処理部11は、車輪90A〜90Dのそれぞれに対応する処理部(1)〜(4)よりなり、例えば処理部(1)では、図4に横軸に時間をとって例示する振動強度データ波形a1を、車輪速v1とは無関係な、予め定められた時間間隔δfでサンプリングし、図5に、横軸を時間軸に沿うデータ順としたグラフで示すように、サンプリングされたデータを、時間順に並べた計測波形X11を出力する。図5において、X11(i-1)、X11(i)、X11(i+1)は、それぞれ、時間順に連続して並んでデータ列を構成する、i-1、i、i+1番目のデータを示す。 The fixed sampling processing unit 11 includes processing units (1) to (4) corresponding to the wheels 90A to 90D. For example, in the processing unit (1), the vibration intensity illustrated by taking time on the horizontal axis in FIG. The data waveform a 1 is sampled at a predetermined time interval δ f that is unrelated to the wheel speed v 1, and is sampled as shown in the graph of FIG. 5 with the horizontal axis representing the data order along the time axis. the data and outputs measurement waveform X 11 which time-sequenced. In FIG. 5, X 11 (i-1), X 11 (i), and X 11 (i + 1) are respectively arranged in a time sequence continuously to form a data string, i-1, i, i + The first data is shown.

一回の抽出工程に対するサンプリング数Mは、予め定められた値を用いるが、適応デジタルフィルタを用いた抽出処理を十分な精度を有するものにするためには、車輪の2〜5回転分を含むとともに、タイヤ一周に対するサンプリング数を表す分割数kを、360〜3600とするのが好ましく、したがって、サンプリング数Mは、720〜18000とするのが好ましい。そして、このとき、固定サンプリングの時間間隔δfは、設計速度100km/h、車輪周長を2m程度として、20〜200μsとするのが好ましい。 The sampling number M for one extraction step uses a predetermined value, but includes 2 to 5 rotations of the wheel in order to make the extraction process using the adaptive digital filter have sufficient accuracy. At the same time, the division number k representing the sampling number for one round of the tire is preferably 360 to 3600, and therefore the sampling number M is preferably 720 to 18000. At this time, the fixed sampling time interval δ f is preferably set to 20 to 200 μs with a design speed of 100 km / h and a wheel circumference of about 2 m.

適応デジタルフィルタ部12は、車輪ごとに設けられたに適応デジタルフィルタ、ADF(1)〜(4)よりなる。図6は、例として、ADF(1)より出力される同期成分波形Y11の一部分を、横軸に時間軸に沿うデータ順をとって示すグラフであり、図6において、Y11(i-1)、Y11(i)、Y11(i+1)は、それぞれ、データ列のうち時間順に連続して並んだ、i-1、i、i+1 番目のデータを示す。 The adaptive digital filter unit 12 includes an adaptive digital filter, ADF (1) to (4), provided for each wheel. Figure 6 shows, as an example, is a graph showing a portion of a synchronous component waveform Y 11 output from ADF (1), taking data order along the horizontal axis a time axis, in FIG. 6, Y 11 (i- 1), Y 11 (i), and Y 11 (i + 1) respectively indicate i−1, i, and i + 1th data arranged sequentially in time order in the data string.

各ADF、例えば第1輪のADF(1)は、対応するモデル波形R11を入力し、計測波形X11と同期成分波形Y11との差分をE11として、差分E11と入力したモデル波形R11とに基づいて同期成分波形Y11を生成するとともに、同期成分波形Y11を生成する際のフィルタリングパラメータを、差分E11を用いて適応制御するよう構成される。 Each ADF, for example, the first wheel of the ADF (1) receives the corresponding model waveform R 11, a difference between the measured waveform X 11 synchronous component waveform Y 11 as E 11, model waveform type differential E 11 to generate a synchronizing component waveform Y 11 on the basis of the R 11, the filtering parameters in generating the synchronizing component waveform Y 11, configured to adaptively controlled using the difference E 11.

ここで、差分E1rをフィルタリングパラメータにフィードバックして最適化を図る方法としては、従来からフィルタリングパラメータ更新アルゴリズムとして知られているLMS(Least Mean Square)法、ニュートン法あるいは最急下法などの適応アルゴリズムを用いることができる、また、その他好適に利用できる適応アルゴリズムとして、複素LMSアルゴリズム(Complex Least Mean Square Algorithm)、Normalized LMSアルゴリズム(Normalized Least Mean Square Algorithm)、射影アルゴリズム(Projection Algorithm)、SHARFアルゴリズム(Simple Hyperstable Adaptive Recursive Filter Algorithm)、RLSアルゴリズム(Recursive Least Square Algorithm)、FLMSアルゴリズム(Fast Least Mean Square Algorithm)、DCTを用いた適応フィルタ(Adaptive Filter using Discrete Cosine Transform)、SANフィルタ(Single Frequency Adaptive Notch Filter)、ニューラルネットワーク(Neural Network)、遺伝的アルゴリズム(Genetic Algorithm)を用いることができる。 Here, as a method of optimizing by feeding back the difference E 1r to the filtering parameter, adaptation such as the LMS (Least Mean Square) method, Newton method, or the steepest bottom method, which is conventionally known as a filtering parameter update algorithm, is used. Algorithms that can be used, and other adaptive algorithms that can be suitably used include Complex LMS Algorithm (Complex Least Mean Square Algorithm), Normalized LMS Algorithm (Normalized Least Mean Square Algorithm), Projection Algorithm, SHARF Algorithm ( Simple Hyperstable Adaptive Recursive Filter Algorithm), RLS algorithm (Recursive Least Square Algorithm), FLMS algorithm (Fast Least Mean Square Algorithm), adaptive filter using DCT (Adaptive Filter using Discrete Cosine Transform), SAN filter (Single F requency Adaptive Notch Filter), neural network, and genetic algorithm can be used.

また、モデル波形R11〜R14は、計測波形X11〜X14の中から、その車輪の回転に同期した成分だけを抽出する際、抽出対象とする部分を打ち抜くのに用いる打ち抜き型として機能するものであり、以下にこの波形について詳細を説明するが、以下の説明において、車輪一回転分の周期に対応する周波数を基本周波数と呼び、この基本周波数に1を含む自然数nを乗じた周波数を、その車輪に対するn次の同期周波数と呼ぶこととする。 Further, the model waveform R 11 to R 14, from the measured waveform X 11 to X 14, when extracting only a component synchronized with the rotation of the wheel, functions as a cutting die used to punch the portion to be extracted The details of this waveform will be described below. In the following description, the frequency corresponding to the period of one wheel rotation is called a fundamental frequency, and the fundamental frequency is multiplied by a natural number n including 1. Is referred to as the nth order synchronization frequency for that wheel.

モデル波形R11〜R14のそれぞれは、モデル波形生成部13の、各車輪に対応する生成部(1)〜(4)で生成される。図7は、例として第1輪に対する生成部(1)で生成される同期成分波形R11の一部分を、横軸に時間軸に沿うデータ順をとって示すグラフであり、図7において、R11(i-1)、R11(i)、R11(i+1)は、それぞれ、データ列のうち、時間順に連続して並ぶ、i-1、i、i+1番目のデータを示す。 Each of the model waveforms R 11 to R 14 is generated by the generation units (1) to (4) corresponding to the wheels of the model waveform generation unit 13. FIG. 7 is a graph showing, as an example, a part of the synchronous component waveform R 11 generated by the generation unit (1) for the first wheel, with the horizontal axis representing the data order along the time axis. 11 (i-1), R 11 (i), and R 11 (i + 1) respectively indicate i-1, i, and i + 1th data arranged sequentially in time order in the data string. .

そして、各車輪に対するモデル波形R11〜R14の各データR11(i)〜R14(i)を、振動強度データa1〜a4のサンプリングに用いた固定の時間間隔δfで並べてできる時間領域波形をフーリエ変換した、周波数領域波形は、図8に、モデル波形R11に対応するものを、周波数を横軸にとって例示するように、同期成分波形Y11の抽出に係る次数の同期周波数の位置で山となり、隣接する次数の同期周波数同士の間では谷となる形状をなし、前記山となる部分の頂部は無限小もしくは有限の幅の周波数範囲に点もしくは平坦面として延在し、頂部の高さは、前記抽出に係るすべての次数の同期周波数に対して同じ値となる特性を有する。 The data R 11 (i) to R 14 (i) of the model waveforms R 11 to R 14 for each wheel can be arranged at a fixed time interval δ f used for sampling the vibration intensity data a 1 to a 4. Fourier transform the time domain waveform, the frequency domain waveform, Figure 8, those corresponding to the model waveform R 11, as illustrated frequency abscissa, the following number of synchronization frequency according to the extracted synchronization component waveform Y 11 It becomes a mountain at the position of, and forms a valley between adjacent synchronization frequencies of the order, the top of the portion that becomes the mountain extends to a frequency range of infinitesimal or finite width as a point or flat surface, The height of the top has the same value with respect to the synchronization frequencies of all orders related to the extraction.

すなわち、図8において、車輪90Aのモデル波形R11に対応する周波数波形は、車輪90Aの車輪速センサからのデータに基づいて算出された基本周波数f1 (1)の整数倍n(n=1〜20)の周波数f1 (1)〜f1 (20)(ただし、f1 (n)=n・f1 (1))の位置で山となり、これら山の間では谷となる形状を有し、図8に示したものの場合、山となる部分の頂部は無限小の点として配置され、頂部の高さは、f1 (1)〜f1 (20)の周波数のすべての位置において、同じ値となる形状をなしている。 That is, in FIG. 8, the frequency waveform corresponding to the model waveform R 11 of the wheel 90A is an integer multiple n (n = 1 ) of the basic frequency f 1 (1) calculated based on the data from the wheel speed sensor of the wheel 90A. ~ 20) with a frequency f 1 (1) ~ f 1 (20) (where f 1 (n) = n · f 1 (1) ) In the case of what is shown in FIG. 8, the top of the peak portion is arranged as an infinitely small point, and the height of the top is at all the positions of the frequencies f 1 (1) to f 1 (20) . The shape has the same value.

一方、算出された基本周波数f1 (1)から、時間順にならべられたデータ列となるモデル波形R11を生成するには、例えば、式(1)によって時間領域波形R1(t)を算出し、算出された波形を時間間隔δfでサンプリングして時間軸に沿うデータ列とすればよい。ただし、Nは、前記抽出に係る最大次数を表し、Aは、予め定めた定数である。

Figure 2006224757
On the other hand, in order to generate the model waveform R 11 that is a data sequence arranged in time order from the calculated fundamental frequency f 1 (1) , for example, the time domain waveform R 1 (t) is calculated by Equation (1). Then, the calculated waveform may be sampled at a time interval δ f to form a data string along the time axis. Here, N represents the maximum order related to the extraction, and A 1 is a predetermined constant.
Figure 2006224757

上記のように構成された抽出手段3Aの、車輪の回転に同期する成分の抽出の成否を確認するため、走行中の車輪90Aの振動強度データa1から固定の時間間隔δfでサンプリングされた計測波形X11、および、デジタルフィルタ部から出力される同期成分波形Y11のそれぞれを、これらの波形を構成する各データX11(i)、およびY11(i)を、時間間隔δfで並べて時間領域波形を生成し、さらにこれらをフーリエ変換して周波数領域波形を求めた。図9は、計測波形X11に対応する周波数領域波形であり、図10は、同期成分波形Y11に対応する周波数領域波形を示すグラフである。 In order to confirm the success or failure of extraction of the component synchronized with the rotation of the wheel, the extraction means 3A configured as described above was sampled at a fixed time interval δ f from the vibration intensity data a 1 of the traveling wheel 90A. Each of the measurement waveform X 11 and the synchronous component waveform Y 11 output from the digital filter unit, and the data X 11 (i) and Y 11 (i) constituting these waveforms at the time interval δ f A time domain waveform was generated side by side, and these were further Fourier transformed to obtain a frequency domain waveform. Figure 9 is a frequency domain waveform corresponding to the measured waveform X 11, FIG. 10 is a graph showing a frequency domain waveform corresponding to the synchronizing component waveform Y 11.

図9と、図10とを対比すると、図10は、同期周波数f1 (1)〜f1 (20)において、図9における振動強度と同じ値を有する高さの山よりなり、これらの間には谷部が形成された形の波形をなしており、計測波形X11から車輪の回転に同期する波形、すなわち、この車輪の同期周波数よりなる波形を抽出することができたことを示しており、異常警告装置4では、同期周波数f1 (1)〜f1 (20)における山の高さを判定する処理、もしくは、これと同等の処理を行うことにより、この高さが予め定めた限界値を超えた場合には、車輪に何かの異常が発生したと判定するよう構成することができる。 When FIG. 9 is compared with FIG. 10, FIG. 10 is composed of peaks having a height having the same value as the vibration intensity in FIG. 9 at the synchronization frequencies f 1 (1) to f 1 (20) . and forms a waveform in the form of troughs are formed in the waveform to synchronize the measurement waveform X 11 to the rotation of the wheel, i.e., it indicates that it is possible to extract the waveform consisting synchronizing frequency of the wheel In the abnormality warning device 4, this height is determined in advance by performing processing for determining the height of the peaks at the synchronization frequencies f 1 (1) to f 1 (20) or processing equivalent thereto. When the limit value is exceeded, it can be configured to determine that some abnormality has occurred in the wheel.

なお、図9の計測波形は、サンプリングの際の時間間隔δfを140μsとし、サンプリング数Mを2048としたサンプリング処理に基づくものであり、また、測定時の基本周波数は13.5Hzであった。 The measurement waveform in FIG. 9 is based on a sampling process in which the time interval δ f at the time of sampling is 140 μs and the sampling number M is 2048, and the fundamental frequency at the time of measurement is 13.5 Hz.

図11は、第一実施形態の抽出手段3Aが実行する処理を示すフローチャートであり、抽出手段3は、車輪数を4として、まず、車輪番号rを1に設定し(ステップs101)、第1車輪について、車輪速データv1を入力し(ステップs102)、次いで、車輪速データv1に基づいて第1車輪の基本周波数を算出し、モデル波形R11(i) (i=1〜M、ただし、Mは1回の抽出に用いるサンプリング数)を生成する(ステップs103)。この生成にあったっては、先に説明した手順のとおり、式(1)を用いて、時間領域波形を求め、これを時間間隔δfで離散化してデータ列とすればよい。 FIG. 11 is a flowchart showing the processing executed by the extraction means 3A of the first embodiment. The extraction means 3 first sets the wheel number r to 1 with the number of wheels set to 4 (step s101). For the wheels, wheel speed data v 1 is input (step s102), then the basic frequency of the first wheel is calculated based on the wheel speed data v 1 and the model waveform R 11 (i) (i = 1 to M, However, M generates the sampling number used for one extraction) (step s103). For this generation, as in the procedure described above, the time domain waveform is obtained by using the equation (1), and this is discretized at the time interval δ f to form a data string.

次いで、モデル波形R11のi番目のデータR11(i)と計測波形X11のi番目のデータX11 (i)とから同期成分波形Y11のi番目のデータY11 (i)を算出するが、これは次のような手順に従う。すなわち、まず、同期成分波形を初期化するとともに(ステップs104)、データ番号を初期化してi=1としたあと(ステップs105)、生成部(1)で生成されたモデル波形のデータR11(1)をADF(1)に入力し(ステップs106)、一方、計測波形に関し、まず、振動強度データa1を取り込みながら(ステップs107)、これを時間間隔δfでサンプリングして、計測波形X11の1番目のデータX11 (1)を生成し(ステップs108)、X11 (1)と初期化済みのY11 (0)とから差分E11(1)を求め(ステップs109)これをADF(1)に入力する。 Then, calculate the i-th data Y 11 synchronization component waveform Y 11 from the i-th data X 11 of the i-th data R 11 (i) and the measured waveform X 11 of model waveform R 11 (i) (i) This follows the following procedure. That is, first, the synchronous component waveform is initialized (step s104), the data number is initialized to i = 1 (step s105), and then the model waveform data R 11 ( 1) is input to the ADF (1) (step s106). On the other hand, the measurement waveform X is first sampled at the time interval δ f while acquiring the vibration intensity data a 1 (step s107). 11 first data X 11 (1) is generated (step s108), and difference E 11 (1) is obtained from X 11 (1) and initialized Y 11 (0) (step s109). Input to ADF (1).

ADF(1)では、差分E11(1)によってフィルタリングパラメータを最適化処理するとともに、先に入力したR11(1)と、差分E11(1)とから、同期成分波形Y11の1番目のデータY11 (1)を求め(ステップs110)、これを出力して、次の差分E11(2)の算出に用いるとともに、異常警告装置4に出力する(ステップs111)。 In ADF (1), together with processes optimized filtering parameters by the difference E 11 (1), and R 11 (1) entered previously, because the difference E 11 (1), the first synchronization component waveform Y 11 Data Y 11 (1) is obtained (step s110), and is output for use in calculation of the next difference E 11 (2) and is output to the abnormality warning device 4 (step s111).

以上で、Y11 (1)を求める一巡のステップを終了して、i=i+1として(ステップs112)、Y11 (2)を求めるステップにはいるが、その前にiがサンプリング数Mに達しているかを判定し(ステップs113)、もし達している場合には、第1輪についての同期成分波形Y11の抽出を完了し、r=r+1として(ステップs114)、第2輪について、同様の手順で同期成分波形Y11の抽出を行い、これを最終輪まで(ステップs115)繰り返えして一回の抽出工程を完了する。そして、走行中、この抽出工程を所定周期で繰り返し行うが、この繰り返しを短周期で行えば行うほど、よりリアルタイムな異常の検出を行うことができる。 With the above, one round of steps for obtaining Y 11 (1) is completed, i = i + 1 is set (step s112), and Y 11 (2) is entered, but before that, i is the sampling number M determines whether the reached (step s113), if they reached if completes the extraction of the sync component waveform Y 11 for the first wheel, as r = r + 1 (step s114), the second wheel for, it was extracted synchronization component waveform Y 11 in the same procedure, which to a final wheel (step s115) repeatedly Kaee to complete the single extraction step. And while driving | running | working, this extraction process is repeatedly performed with a predetermined period, The more real-time abnormality detection can be performed, so that this repetition is performed with a short period.

なお、各車輪についての処理を、説明の都合上、時間的に直列したものとしたが、装置の処理時間が許容する限りにおいて、並列処理とすることができ、この場合、1回の抽出工程に要する時間を短くすることができる。   In addition, although the processing for each wheel is assumed to be in series in time for convenience of explanation, it can be parallel processing as long as the processing time of the apparatus is allowed. In this case, one extraction step Can be shortened.

以上説明したように、同期成分波形Y1rの算出に際して、入力した振動強度データarを、フーリエ変換等、時間のかかる処理を行うことなく、所定の時間間隔δfでサンプリングしたものをそのまま、入力順に処理して同期成分波形Y11を生成することができ、このことによって、前記抽出の精度を犠牲にすることなく、1回の抽出工程に要する時間を短いものとすることができる。 As described above, when calculating the synchronous component waveform Y 1r , the input vibration intensity data a r is sampled at a predetermined time interval δ f without performing time-consuming processing such as Fourier transform, processed input order to be able to generate a synchronizing component waveform Y 11, by this, without sacrificing accuracy of the extraction, may be shorter the time required for one extraction step.

次に、以上に説明した第一の実施形態の車輪異常検知装置を基本形として、抽出工程をさらに短縮化した第二の実施形態について説明する。この実施形態は、第一実施形態の抽出手段3Aを、抽出手段3Bに置換した構成になり、図3にブロック線図で示すように、抽出手段3Bは、振動計測手段1によって連続的に計測された各車輪についての振動強度データa1〜a4を、予め定められた固定の時間間隔δfでサンプリングして、時間順に並んだデータ列よりなる計測波形X21〜X24を生成する固定サンプリング処理部21と、各車輪について、計測波形X21〜X24の中から、その車輪の回転に同期した成分だけを、時間順に並んだデータ列として抽出し、これを、同期成分波形b1〜b4(もしくはY21〜Y24)として前記異常警告手段4に出力する適応デジタルフィルタ部22と、各車輪の前記抽出のモデルとなるモデル波形R2を、時間順に並んだデータ列として生成するモデル波形生成部23とよりなる。 Next, a second embodiment in which the extraction process is further shortened will be described using the wheel abnormality detection device of the first embodiment described above as a basic form. In this embodiment, the extraction means 3A of the first embodiment is replaced with an extraction means 3B. The extraction means 3B is continuously measured by the vibration measurement means 1 as shown in a block diagram in FIG. the vibration intensity data a 1 ~a 4 for each wheel, which is, by sampling at a predetermined fixed time interval [delta] f, fixed to generate a measurement waveform X 21 to X 24 consisting of a data string arranged in time order For the sampling processing unit 21 and each wheel, only the component synchronized with the rotation of the wheel is extracted from the measurement waveforms X 21 to X 24 as a data sequence arranged in time order, and this is extracted as the synchronization component waveform b 1. The adaptive digital filter unit 22 that outputs to the abnormality warning means 4 as ˜b 4 (or Y 21 ˜Y 24 ) and the model waveform R 2 that becomes the model of the extraction of each wheel are generated as a data sequence arranged in time order Model waveform raw It becomes more a part 23.

図3に示した、第一実施形態の抽出手段3Aに対するブロック線図と、図12に示した、本実施形態の抽出手段3Bに対するブロック線図とを比較すると明らかなように、これらの相違点は、第一実施形態における適応デジタルフィルタ部12の各ADF(r)は、車輪ごとに異なるモデル波形R1rを入力するのに対して、第二の実施形態においては、これらの車輪に共通の一つのモデル波形R2を入力してこれを各ADF(r)における同期成分波形Y2rの生成に用いるところにあり、このことによって、後述する通り、抽出に掛かる処理時間の更なる短縮化を実現しようとするものである。 As is apparent from a comparison between the block diagram for the extracting unit 3A of the first embodiment shown in FIG. 3 and the block diagram for the extracting unit 3B of the present embodiment shown in FIG. The ADF (r) of the adaptive digital filter unit 12 in the first embodiment inputs a different model waveform R 1r for each wheel, whereas in the second embodiment, it is common to these wheels. One model waveform R 2 is input and used to generate a synchronous component waveform Y 2r in each ADF (r). This further reduces the processing time required for extraction, as will be described later. It is something that is going to be realized.

このためのモデル波形R2が具備すべき要件としては、これが各車輪に対して共通化したものであっても、すべての車輪に対して同期成分波形Y2rを高精度に抽出できるものでなくてはならということであり、以下に、この要件を満足させることのできるモデル波形R2の構成について説明する。 As a requirement that the model waveform R 2 for this purpose should have, the synchronous component waveform Y 2r cannot be extracted with high accuracy for all the wheels even if it is common to each wheel. In the following, the configuration of the model waveform R 2 that can satisfy this requirement will be described.

図13は、モデル波形生成部23で生成されるモデル波形R2の一部分を、横軸に時間軸に沿うデータ順をとって示すグラフであり、図7において、R2(i-1)、R2(i)、R2(i+1)は、それぞれ、データ列のうち、時間順に連続して並ぶ、i-1、i、i+1番目のデータを示す。 FIG. 13 is a graph showing a part of the model waveform R 2 generated by the model waveform generation unit 23 with the horizontal axis representing the data order along the time axis, and in FIG. 7, R 2 (i−1), R 2 (i) and R 2 (i + 1) respectively indicate i−1, i, and i + 1th data that are continuously arranged in time order in the data string.

そして、この実施形態において、各車輪に対して共通なモデル波形R2の各データR2(i)を、前記固定の時間間隔δfで並べた時間領域波形をフーリエ変換してできる周波数領域波形は、図14に示すように、どの車輪に対しても、抽出に係る次数の同期周波数の位置で山となり、隣接する次数の同期周波数同士の間では谷となる形状をなし、前記山となる部分の頂部は無限小もしくは有限の幅の周波数範囲に点もしくは平坦面として延在し、頂部の高さは、前記抽出に係るすべての次数の同期周波数に対して同じ値となる特性を有し、この点においては、第一の実施形態に用いられるそれぞれのモデル波形R1r(r=1〜4)と同様であるが、モデル波形R2がこれらの波形R1rと異なる点は、頂部が点ではなく、有限の幅を有する平坦面として延在する点であり、これは、モデル波形R2を、車輪速の異なるどの車輪に対しても同期周波数成分を抽出することができる打ち抜き型として機能させなければならないからであり、そのため、モデル波形R2のn次の頂部を、すべての車輪に対する同期周波数fr (n)を含むよう延在させたものである。 In this embodiment, each data R 2 (i) of the model waveform R 2 common to each wheel is a frequency domain waveform obtained by Fourier transforming a time domain waveform arranged at the fixed time interval δ f. As shown in FIG. 14, for any wheel, it becomes a mountain at the position of the synchronization frequency of the order of extraction, and forms a valley between the synchronization frequencies of the adjacent orders, and becomes the mountain The top of the part extends as a point or flat surface in a frequency range of infinitesimal or finite width, and the height of the top has the same value for all the synchronization frequencies of the extraction. This point is the same as each model waveform R 1r (r = 1 to 4) used in the first embodiment, except that the model waveform R 2 is different from these waveforms R 1r at the top. Extends as a flat surface with a finite width, not a point , And the this is a model waveform R 2, is because must function as a punching die which can be extracted synchronization frequency component for any wheel of different wheel speeds, therefore, the model waveform R 2 The top of the nth order is extended to include the synchronization frequency f r (n) for all wheels.

すなわち、全車輪のうち、最高車速の車輪に対する基本周波数(周波数は最小)をfminとし、最低車速の車輪に対する基本周波数(周波数は最大)をfmaxとしたとき、モデル波形R2のn次の頂部の延在領域の上限fu (n)および下限fw (n)は、それぞれ、式(2)および(3)の条件を満足することが必要である。

Figure 2006224757
That is, of all the wheels, when the fundamental frequency (frequency is minimum) for the wheel at the highest vehicle speed is f min and the fundamental frequency (frequency is maximum) for the wheel at the lowest vehicle speed is f max , the nth order of the model waveform R 2 It is necessary that the upper limit f u (n) and the lower limit f w (n) of the extension region at the top of the above satisfy the conditions of the expressions (2) and (3), respectively.
Figure 2006224757

一方、車輪速データv1〜v4を基に算出された基本周波数f1 (1)〜f1 (4)から、時間順にならべられたデータ列となるモデル波形R2を生成するには、これらの基本周波数f1 (1)〜f1 (4)から、上記式(2)および(3)に基づいて、n次の頂部の延在領域をそれぞれ設定し、この設定領域を適当な個数J、例えば20個の周波数の集まりからなるものと近似して、これらの周波数に対する正弦波を足し合わせて、式(4)に示すように、時間領域波形R2(t)を生成し、この時間領域波形を固定サンプリングタイムδfでサンプリングして離散データを作り、モデル波形とすればよい。なお、式(5)における、fw(1)、およびfu(1)は、1次の頂部の延在領域に対する下限周波数、および上限周波数を表し、A2は予め定められた定数である。

Figure 2006224757
On the other hand, in order to generate a model waveform R 2 that is a data sequence arranged in time order from the basic frequencies f 1 (1) to f 1 (4) calculated based on the wheel speed data v 1 to v 4 , From these fundamental frequencies f 1 (1) to f 1 (4) , an extension region of the n-th apex is set based on the above equations (2) and (3), and an appropriate number of these setting regions are set. J, for example, approximating that it consists of a collection of 20 frequencies, and adding the sine waves for these frequencies to generate a time domain waveform R 2 (t) as shown in equation (4), A time domain waveform may be sampled at a fixed sampling time δ f to create discrete data, which is used as a model waveform. In Equation (5), fw (1) and fu (1) represent a lower limit frequency and an upper limit frequency for the extension region of the primary apex, and A 2 is a predetermined constant.
Figure 2006224757

以上のように、モデル波形R2を周波数領域に変換してできる波形は、有限の幅の頂部を有し、この頂部の延在領域は、同期成分波形の抽出に係る複数の車輪のうち、最高速の車輪の基本周波数fminと、最低速の車輪の基本周波数fmaxとの両方に基づいて決めることができるが、より簡便な手段として、この代わりに、Aを予め定められた数値として、fminの0〜+A%、もしくは、fmaxの-A〜0%とするもこともできる。 As described above, the waveform can convert the model waveform R 2 in the frequency domain, it has a top portion of finite width, extending region of the top portion of the plurality of wheels according to the extraction of the synchronizing component waveform, It can be determined based on both the fundamental frequency f min of the fastest wheel and the fundamental frequency f max of the slowest wheel, but as a simpler means, instead of this, A is a predetermined numerical value. , F min of 0 to + A%, or f max of −A to 0%.

このように設定した場合、fmaxもしくはfminのいずれか一つの車輪速だけに基づいてより簡単に頂部の延在周波数領域を決めることができので好ましく、特に、車両の制御装置として、最高速の車輪の車輪側に基づいて車両の制御を行うものが車両に設けられている場合には、最高速の車輪の基本周波数fminだけに基づいて頂部の延在周波数領域を決めることにより、本装置以外の装置で選択された車輪番号、もしくは、本装置以外の装置で収集された、最高速車輪の基本周波数fminやこれに相関するデータを流用することにより、本装置を極めて簡便なものにすることができる。 When set in this way, it is preferable because the extended frequency region of the top can be determined more easily based on only one wheel speed of f max or f min. If the vehicle is equipped with a vehicle that controls the wheel based on the wheel side of the wheel, the extension frequency region of the top is determined based only on the fundamental frequency f min of the fastest wheel. By using the wheel number selected by a device other than the device, or the fundamental frequency f min of the highest speed wheel collected by a device other than this device and the data correlated therewith, this device is extremely simple. Can be.

また、車輪速の変わりに、車輪速に代わってこれを表す特性、例えば、車体ヨーレート等で、他の目的のためにリアルタイムに収集される特性値があれば、これらを転用することができる。   Also, if there are characteristic values collected in real time for other purposes such as a characteristic representing the wheel speed instead of the wheel speed, such as a vehicle body yaw rate, these can be diverted.

ここで、前記数値Aは、車両のどんな走行条件下でも、どの車輪に対しても、n次の同期周波数が頂部の延在範囲内に位置するように設定すればよく、例えば、5〜15(%)が好ましい範囲である。   Here, the numerical value A may be set so that the n-th order synchronization frequency is located within the extension range of the top for any wheel under any driving condition of the vehicle. (%) Is a preferred range.

上記のように構成された抽出手段3Bの、車輪の回転に同期する成分の抽出の成否を確認するため、走行中の第1車輪90Aの振動強度データa1から固定の時間間隔δfでサンプリングされた計測波形X21、および、デジタルフィルタ部から出力される同期成分波形Y21のそれぞれを、これらの波形を構成する各データX21(i)、およびY21(i)を、時間間隔δfで並べて時間領域波形を生成し、さらにこれらをフーリエ変換して周波数領域波形を求めた。計測波形X21に対応する周波数領域波形は、先に図9に示した通りであり、同期成分波形Y21に対応する周波数領域波形は、これを図15に示した。 In order to confirm the success or failure of extraction of the component synchronized with the rotation of the wheel by the extraction means 3B configured as described above, sampling is performed at a fixed time interval δ f from the vibration intensity data a 1 of the traveling first wheel 90A. The measured measurement waveform X 21 and the synchronization component waveform Y 21 output from the digital filter unit, the data X 21 (i) and Y 21 (i) constituting these waveforms, and the time interval δ A time domain waveform was generated by arranging them at f , and these were further Fourier transformed to obtain a frequency domain waveform. The frequency domain waveform corresponding to the measurement waveform X 21 is as shown in FIG. 9, and the frequency domain waveform corresponding to the synchronization component waveform Y 21 is shown in FIG.

また、すべての車輪に共通なモデル波形の比較例として、周波数領域波形における各山の頂部が平坦面となるものの代わりに、第2車輪を代表車輪として、周波数領域波形が第2車輪の同期周波数の位置で点状の頂部を有するモデル波形R22を用いて抽出した場合の同期成分波形Y21に対応する周波数領域波形を図16に示した。 In addition, as a comparative example of the model waveform common to all wheels, the frequency domain waveform is the synchronous frequency of the second wheel with the second wheel as a representative wheel instead of the peak of each mountain in the frequency domain waveform being a flat surface. FIG. 16 shows a frequency domain waveform corresponding to the synchronous component waveform Y 21 when extracted using the model waveform R 22 having a point-like apex at the position.

図9、図15、および図16を比較すると、図15は、同期周波数f1 (1)〜f1 (20)において、図9における振動強度と同じ値を有する高さの山をもち、これらの間には谷部が形成された形の波形をなしており、各車輪について、計測波形X11から、車輪の同期周波数よりなる同期成分波形を抽出することができたことを示しているが、図16に示した比較例の周波数領域波形の場合、当然ながら、図9における計測波形X21に対応する周波数領域波形における基本周波数の位置と、抽出に用いたモデル波形R22に対応するそれとがずれていることにより、図15に示したものに対比して、山の高さ、例えばn=1における山の高さが低く、同期成分波形の抽出が不満足なものであることがわかる。 Comparing FIG. 9, FIG. 15, and FIG. 16, FIG. 15 has peaks of height having the same value as the vibration intensity in FIG. 9 at the synchronization frequencies f 1 (1) to f 1 (20) . It has a waveform with a valley formed between them, and it shows that a synchronous component waveform consisting of the synchronous frequency of the wheel could be extracted from the measured waveform X 11 for each wheel. If the frequency domain waveforms of the comparative example shown in FIG. 16, of course, the position of the fundamental frequency in the frequency domain waveform corresponding to the measured waveform X 21 in FIG. 9, it and corresponding to the model waveform R 22 used in the extraction As a result, the height of the peak, for example, the height of the peak at n = 1, is lower than that shown in FIG. 15 and the extraction of the synchronous component waveform is unsatisfactory.

なお、図9の波形は、サンプリングの際の時間間隔δfを140μsとし、サンプリング数Mを2048としたサンプリング処理に基づくものであり、測定時の基本周波数は13.5Hzであった。 The waveform in FIG. 9 is based on a sampling process in which the time interval δ f at the time of sampling is 140 μs and the sampling number M is 2048, and the fundamental frequency at the time of measurement is 13.5 Hz.

図17は、第二実施形態の抽出手段3Bが実行する処理を示すフローチャートであり、抽出手段3Bは、車輪数を4として、まず、モデル波形生成部23に各車輪の車輪速データv1〜v4を入力し(ステップs201)、これに基づいて先に説明した手順に従ってモデル波形R2(i) (i=1〜M、ただし、Mは1回の抽出に用いるサンプリング数)を生成する(ステップs202)。 Figure 17 is a flowchart illustrating a process of extracting means 3B of the second embodiment is executed, the extraction unit 3B as 4 wheel number, first, the model waveform generator 23 on the wheel speed data v 1 ~ of the wheels v 4 is input (step s201), and a model waveform R 2 (i) (i = 1 to M, where M is the number of samplings used for one extraction) is generated based on the procedure described above. (Step s202).

次に、同期成分波形Y2rのデータを一つずつデータ順に求める手順にはいるが、これに先立って、各車輪に対する同期成分波形を初期化しておく(ステップs203)。次いで、i=1としたあと(ステップs204)、モデル波形生成部23で生成したモデル波形R2の1番目のデータR2(i)をADF(1)に入力する(ステップs205)。 Next, although the procedure for obtaining the data of the synchronous component waveform Y 2r one by one in the order of data is performed, prior to this, the synchronous component waveform for each wheel is initialized (step s203). Next, after i = 1 (step s204), the first data R 2 (i) of the model waveform R 2 generated by the model waveform generation unit 23 is input to the ADF (1) (step s205).

次に、このモデル波形データR2(1)を用いて、1〜4輪のすべての車輪に対する同期成分波形の1番目のデータY2r(1)(r=1〜4)を算出するが、これを第1輪から始めるためr=1として(ステップs206)、固定サンプリング処理部21の処理部(1)に、振動強度データa1を入力し(ステップs207)、次いで、これを時間間隔δfでサンプリングして計測波形データX21(i)を生成し(ステップs208)、続いて差分E21(1)を求め(ステップs109)これをADF(1)に入力する。 Next, using this model waveform data R 2 (1), the first data Y 2r (1) (r = 1 to 4) of the synchronous component waveform for all of the 1 to 4 wheels is calculated. as r = 1 to start it from the first wheel (step s206), the processing unit of the fixed sampling processing section 21 (1), enter the vibration intensity data a 1 (step s 207), then this time interval δ Sampling at f generates measurement waveform data X 21 (i) (step s208), and then a difference E 21 (1) is obtained (step s109), which is input to ADF (1).

ADF(1)では、差分E21(1)によってフィルタリングパラメータを最適化処理するとともに、先に入力したR21(1)と、差分E21(1)とから、同期成分波形Y21の1番目のデータY21(1)を求め(ステップs210)、次の差分E21(2)の算出に用いるとともに、異常警告装置4に出力する(ステップs211)。 In ADF (1), the filtering parameter is optimized by the difference E 21 (1), and the first input of the synchronous component waveform Y 21 from the previously input R 21 (1) and the difference E 21 (1). Data Y 21 (1) is obtained (step s210), used to calculate the next difference E 21 (2), and output to the abnormality warning device 4 (step s211).

以上で、Y21 (1)を求める一巡のステップを終了し、次に同様のステップを第2輪以降についても実施するため、r=r+1として(ステップs212)、r=5になるまでこれを繰り返す(ステップs213)。i=1に対する上記の処理を完了すると、i=i+1として(ステップs214)、i=M(ステップs215)になるまで以上を繰り返して一回の抽出工程を完了する。 This completes the round of steps for obtaining Y 21 (1), and then executes the same steps for the second and subsequent wheels, so r = r + 1 (step s212) until r = 5 This is repeated (step s213). When the above processing for i = 1 is completed, i = i + 1 is set (step s214), and the above is repeated until i = M (step s215) to complete one extraction step.

以上説明したように、同期成分波形Y11の算出に際して、入力した振動強度データarを、フーリエ変換等の時間のかかる処理を行うことなく、所定の時間間隔δfでサンプリングしたものをそのまま、入力順に処理して同期成分波形Y1rを生成して、短時間で1回の抽出工程を終えることができる点においては、第一の実施形態と同様であるが、この第二の実施形態の装置では、以下の点において、第一の実施形態のものよりも、処理の時間を短くすることができる。 As described above, when calculating the synchronizing component waveform Y 11, the vibration intensity data a r entered, without performing the time-consuming processing such as Fourier transform, as it is sampled at predetermined time intervals [delta] f, It is the same as the first embodiment in that it can be processed in the order of input to generate the synchronous component waveform Y 1r , and one extraction process can be completed in a short time. In the apparatus, the processing time can be made shorter than that of the first embodiment in the following points.

これは、ステップs205における、モデル波形R2を全車輪について共通のものとしたことによるものであり、このことによる第一の効果は、第一実施形態においては、各車輪の同期成分波形の抽出にそれぞれ別個のモデル波形データR1rを用いているので、各車輪について各同期成分波形データY1r(i)を計算する都度、モデル波形データR1r(i)をメモリあるいはバッファから呼び出さなければならなかったのに対して、本実施形態においては、4輪のY2r(i)(r=1〜4)のすべてに対して1回だけR2(i)をメモリもしくはバッファから呼び出すだけで済ませることができ点である。すなわち、Mをサンプリング数、Rを車輪数として、第一の実施形態のおいては、全部で、M・R回、モデル波形を呼び出す必要があったのに対し、第二の実施形態のおいては、これをM回で済ませることができる。 This is in step s205, is due to the model waveform R 2 were of common to all the wheels, the first effect of this is, in the first embodiment, the extraction of the wheels of the synchronization component waveform since respectively using separate model waveform data R 1r in each time to calculate the respective synchronous component waveform data Y 1r (i) for each wheel, must call model waveform data R 1r and (i) from the memory or buffer In contrast, in this embodiment, R 2 (i) can be called from the memory or buffer only once for all four wheels Y 2r (i) (r = 1 to 4). Can be a point. That is, with M as the number of samplings and R as the number of wheels, in the first embodiment, it was necessary to call the model waveform M · R times in total, whereas in the second embodiment, And you can do this M times.

モデル波形R2を全車輪について共通のものとしたことによる第二の効果は、ステップs202における、モデル波形R2の生成に関し、これに要する処理回数を減らせる点であり、第一の実施形態においては、各車輪に対して別個のモデル波形R1rを生成しなければならなかったところ、第二の実施形態においては、これを共通のモデル波形R2を一回生成するだけで済ませることができ、このことのよっても、抽出工程を短縮することができる。 The second effect obtained by making the model waveform R 2 common to all the wheels is that the number of processes required for generating the model waveform R 2 in step s202 can be reduced. First embodiment In the second embodiment, a separate model waveform R 1r had to be generated for each wheel. However, in the second embodiment, it is only necessary to generate a common model waveform R 2 once. This also makes it possible to shorten the extraction process.

その他の処理部分については、処理のステップ数は同じであり、したがって、上記第一および第二の点において処理数を減らした分だけ、全体の処理時間を短縮化することができる。   For the other processing parts, the number of processing steps is the same, and accordingly, the entire processing time can be shortened by the amount of the processing number reduced in the first and second points.

次に、第三の実施形態について以下に説明する。第三の実施形態は、第一実施形態の抽出手段3Aを、抽出手段3Cに置換した構成になり、図18にブロック線図で示すように、抽出手段3Cは、振動計測手段1によって連続的に計測された各車輪についての振動強度データa1〜a4を、それぞれ、車輪速に応じて変化する時間間隔δvrでサンプリングして、時間順に並んだデータ列よりなる計測波形X31〜X34を生成する可変サンプリング処理部31と、各車輪について、計測波形X31〜X34の中から、その車輪の回転に同期した成分だけを、時間順に並んだデータ列として抽出し、これを、同期成分波形b1〜b4(もしくはY31〜Y34)として前記異常警告手段4に出力する適応デジタルフィルタ部32と、各車輪の前記抽出のモデルとなるモデル波形R3を、時間順に並んだデータ列として生成するモデル波形生成部33とよりなる。 Next, a third embodiment will be described below. The third embodiment has a configuration in which the extraction unit 3A of the first embodiment is replaced with an extraction unit 3C. As shown in a block diagram in FIG. The vibration intensity data a 1 to a 4 measured for each wheel are sampled at a time interval δ vr that changes according to the wheel speed, respectively, and measurement waveforms X 31 to X consisting of a data sequence arranged in time order are sampled. For the variable sampling processing unit 31 that generates 34 and each wheel, only the components synchronized with the rotation of the wheel are extracted from the measurement waveforms X 31 to X 34 as a data sequence arranged in time order, An adaptive digital filter unit 32 that outputs the synchronous component waveforms b 1 to b 4 (or Y 31 to Y 34 ) to the abnormality warning means 4 and a model waveform R 3 that is a model of the extraction of each wheel are arranged in time order. Generate as a data column The more the model waveform generator 33 that.

図3に示した、第一実施形態の抽出手段3Aに対するブロック線図と、図18に示した、本実施形態の抽出手段3Cに対するブロック線図とを比較すると明らかなように、これらの相違点は、第一実施形態においては、振動波形データarを固定の時間間隔δfでサンプリングするとともに、適応デジタルフィルタ部12の各車輪に対応するADFは、車輪ごとに異なるモデル波形R1rを入力する構成としたのに対して、第三の実施形態においては、車輪速に応じて時間とともに変化する可変の時間間隔δvrでサンプリングし、その代わり、ADFには、時間に対して不変の、かつ、各車輪に共通のモデル波形R3を入力する構成とした点にある。 As is apparent from a comparison between the block diagram for the extraction unit 3A of the first embodiment shown in FIG. 3 and the block diagram for the extraction unit 3C of the present embodiment shown in FIG. In the first embodiment, the vibration waveform data a r is sampled at a fixed time interval δ f and the ADF corresponding to each wheel of the adaptive digital filter unit 12 inputs a different model waveform R 1r for each wheel. On the other hand, in the third embodiment, sampling is performed at a variable time interval δ vr that changes with time according to the wheel speed. Instead, the ADF has a time-invariant, In addition, the configuration is such that a common model waveform R 3 is input to each wheel.

可変サンプリング処理部31は、車輪90A〜90Dのそれぞれに対応する処理部(1)〜(4)よりなり、例えば処理部(1)では、図19に横軸に時間をとって例示する振動強度データ波形a1を、時間とともに変化する車輪速v1に応じて可変の時間間隔δv1でサンプリングし、サンプリングされたデータを、図20に、横軸を時間軸に沿うデータ順としたグラフで示すように、時間順に並べたデータ列よりなる計測波形X11を出力する。図5において、X11(i-1)、X11(i)、X11(i+1)は、それぞれ、時間順に連続して並んでデータ列を構成する、i-1、i、i+1番目のデータを示す。 The variable sampling processing unit 31 includes processing units (1) to (4) corresponding to the wheels 90A to 90D. For example, in the processing unit (1), the vibration intensity illustrated with time on the horizontal axis in FIG. The data waveform a 1 is sampled at a variable time interval δ v1 according to the wheel speed v 1 that changes with time, and the sampled data is shown in a graph in FIG. 20, with the horizontal axis representing the data order along the time axis. as shown, it outputs a measured waveform X 11 consisting of data string time-sequenced. In FIG. 5, X 11 (i-1), X 11 (i), and X 11 (i + 1) are respectively arranged in a time sequence continuously to form a data string, i-1, i, i + The first data is shown.

r番目の車輪に対する前記時間間隔δvrは、その車輪の1回転分の周期に対応する基本周波数fr(1)と、車輪1回転を等分する際の分割数となる予め定められた値kとの積の逆数とて設定される。すなわち、車輪が1/k回転するごとにサンプリングを行うよう設定され、これは、車輪1周内のデータ数を一定となるよう設定することによって始めてモデル関数R3を時間に不変のものにすることができるからである。 The time interval δ vr for the r-th wheel is a predetermined value that is a basic frequency f r (1) corresponding to the cycle of one rotation of the wheel and the number of divisions when equally dividing the rotation of the wheel. Set as the inverse of the product of k. That is, it is set to sample every 1 / k rotation of the wheel, and this is the first time that the model function R 3 is made invariant with time by setting the number of data in one wheel revolution to be constant. Because it can.

一回の抽出工程に対するサンプリング数Mは、予め定められた値を用いるが、適応デジタルフィルタを用いた抽出処理を十分な精度を有するものにするためには、車輪の2〜5回転分を含むとともに、タイヤ一周に対するサンプリング数を表す分割数kを、360〜3600とするのが好ましく、したがって、サンプリング数Mは、720〜18000とするのが好ましい。   The sampling number M for one extraction step uses a predetermined value, but includes 2 to 5 rotations of the wheel to make the extraction process using the adaptive digital filter have sufficient accuracy. At the same time, the division number k representing the sampling number for one round of the tire is preferably 360 to 3600, and therefore the sampling number M is preferably 720 to 18000.

なお、可変時間間隔δvrを算出するため、各車輪に対して、入力した車輪速vrから時間間隔δvrを算出するδv算出部(r)が、可変サンプリング処理部31の処理部(r)の前に設けられる。 In order to calculate the variable time interval δ vr , for each wheel, a δ v calculation unit (r) for calculating the time interval δ vr from the input wheel speed v r is a processing unit of the variable sampling processing unit 31 ( provided before r).

適応デジタルフィルタ部12は、それぞれの車輪に対応する適応デジタルフィルタであるADF(1)〜(4)よりなり、各ADF、例えば第1輪のADF(1)は、固定のモデル波形R3を入力し、計測波形X31と同期成分波形Y31との差分をE31として、差分E31と入力したモデル波形R31とに基づいて同期成分波形Y31を生成するとともに、同期成分波形Y31を生成する際のフィルタリングパラメータを、差分E31を用いて適応制御するよう構成される点において、第一の実施形態について説明したのと同様であり、詳細の説明を省略する。 Adaptive digital filter unit 12 is made of an adaptive digital filter corresponding to the respective wheels ADF (1) ~ (4) , each ADF, for example, the first wheel of the ADF (1) is a model waveform R 3 fixed type, the E 31 the difference between the measured waveform X 31 and synchronizing component waveform Y 31, and generates a synchronizing component waveform Y 31 on the basis of the model waveform R 31 input the difference E 31, synchronization component waveform Y 31 Is the same as that described in the first embodiment in that it is configured to adaptively control the filtering parameter when generating the difference E 31 , and detailed description thereof is omitted.

そして、予め定められた固定の周波数を固定基本周波数f0として、この固定基本周波数に自然数nを乗じた周波数を、「n次の固定同期周波数」とし、前記固定基本周波数f0と前記分割数kとの積の逆数を時間間隔δ0としたとき、モデル波形生成部33で生成される固定のモデル波形R3は、このモデル波形R3の各データを前記時間間隔δ0おきに並べた時間領域波形をフーリエ変換してできる周波数領域波形が、前記抽出に係る次数nの固定同期周波数f0 (n)=n・f0の位置で山となり、隣接する次数の固定同期周波数同士の間では谷となる形状をなし、前記山となる部分の頂部は無限小もしくは有限の幅の周波数範囲で点もしくは平坦面として延在し、頂部の高さは、前記抽出に係るすべての次数の固定同期周波数に対して同じ値となる特性を有するよう構成される。 A predetermined fixed frequency is defined as a fixed fundamental frequency f 0 , and a frequency obtained by multiplying the fixed fundamental frequency by a natural number n is referred to as an “nth-order fixed synchronization frequency”, and the fixed fundamental frequency f 0 and the number of divisions When the reciprocal of the product of k is the time interval δ 0 , the fixed model waveform R 3 generated by the model waveform generation unit 33 is obtained by arranging the data of the model waveform R 3 every time interval δ 0 . The frequency domain waveform formed by Fourier transform of the time domain waveform is a peak at the position of the fixed synchronization frequency f 0 (n) = n · f 0 of the extraction order n, and between the fixed synchronization frequencies of the adjacent orders. In the shape of a valley, the peak of the peak extends as a point or flat surface in a frequency range of infinitesimal or finite width, and the height of the peak is fixed for all orders related to the extraction. Has the same value for the synchronization frequency As constructed.

すなわち、図21において、モデル波形R3に対応する周波数波形は、固定基本周波数f0の整数倍n(n=1〜20)の周波数f0 (1)〜f0 (20)(ただし、f0 (n)=n・f0 (1))の位置で山となり、これら山の間では谷となる形状を有し、図8に示したものの場合、山となる部分の頂部は無限小の点として配置され、頂部の高さは、f0 (1)〜f0 (20)の周波数のすべてにおいて、同じ値となる形状をなしている。 That is, in FIG. 21, the frequency waveform corresponding to the model waveform R 3 has frequencies f 0 (1) to f 0 (20) (where f is an integer multiple n (n = 1 to 20) of the fixed fundamental frequency f 0 ). 0 (n) = n · f 0 (1) ), and has a shape that becomes a valley between these peaks. In the case shown in FIG. 8, the top of the portion that becomes the peak is infinitely small. It is arranged as a point, and the height of the top portion has a shape having the same value at all frequencies of f 0 (1) to f 0 (20) .

一方、固定基本周波数f0から、時間順にならべられたデータ列となるモデル波形R3を生成するには、例えば、式(6)によって時間領域波形R3(t)を算出し、算出された波形を時間間隔δfでサンプリングして時間軸に沿うデータ列とすればよい。ただし、Nは、前記抽出に係る最大次数を表し、A3は、予め定めた定数である。

Figure 2006224757
On the other hand, in order to generate the model waveform R 3 that is a data sequence arranged in time order from the fixed fundamental frequency f 0 , for example, the time domain waveform R 3 (t) is calculated by the equation (6). The waveform may be sampled at a time interval δ f to form a data string along the time axis. Here, N represents the maximum order related to the extraction, and A 3 is a predetermined constant.
Figure 2006224757

上記のように構成された抽出手段3Cの、車輪の回転に同期する成分の抽出の成否を確認するため、走行中の車輪90Aの振動強度データa1から可変時間間隔δv1でサンプリングされた計測波形X31、および、デジタルフィルタ部から出力される同期成分波形Y31のそれぞれを、これらの波形を構成する各データX31(i)、およびY31(i)を、固定の時間間隔δ0で並べて時間領域波形を生成し、さらにこれらをフーリエ変換して周波数領域波形を求めた。図22は、計測波形X31に対応する周波数領域波形であり、図23は、同期成分波形Y31に対応する周波数領域波形を示すグラフである。 In order to confirm the success or failure of extraction of the component synchronized with the rotation of the wheel by the extraction means 3C configured as described above, the measurement sampled at the variable time interval δ v1 from the vibration intensity data a 1 of the traveling wheel 90A. Each of the waveform X 31 and the synchronous component waveform Y 31 output from the digital filter unit is replaced with each data X 31 (i) and Y 31 (i) constituting these waveforms at a fixed time interval δ 0. Were arranged side by side to generate a time domain waveform, and these were further Fourier transformed to obtain a frequency domain waveform. Figure 22 is a frequency domain waveform corresponding to the measured waveform X 31, FIG. 23 is a graph showing a frequency domain waveform corresponding to the synchronizing component waveform Y 31.

図22と、図23とを対比すると、図23は、同期周波数f0 (1)〜f0 (20)において、図22における振動強度と同じ値を有する高さの山よりなり、これらの間には谷部が形成された形の波形をなしており、計測波形X31から車輪の回転に同期する波形、すなわち、車輪の同期周波数よりなる波形を抽出することができたことを示しており、異常警告装置4では、同期周波数f1 (1)〜f1 (20)における山の高さを判定する処理、もしくは、これと同等の処理を行うことにより、この高さが予め定めた限界値を超えた場合には、車輪に何かの異常が発生したと判定することができる。 FIG. 23 is compared with FIG. 23. FIG. 23 is composed of peaks having a height having the same value as the vibration intensity in FIG. 22 at the synchronization frequencies f 0 (1) to f 0 (20) . Shows a waveform with a valley formed, and shows that the waveform synchronized with the rotation of the wheel, that is, the waveform consisting of the synchronization frequency of the wheel was extracted from the measured waveform X 31. In the abnormality warning device 4, this height is set to a predetermined limit by performing processing for determining the height of the peaks at the synchronization frequencies f 1 (1) to f 1 (20) or processing equivalent thereto. When the value is exceeded, it can be determined that some abnormality has occurred in the wheel.

なお、上記算出に用いた、振動強度データa1のサンプリングにおける時間間隔δv1は、118μsとし、固定基本周波数f0は16.5Hz、分割数kは512であった。 Incidentally, used in the above calculation, the time interval [delta] v1 in the sampling of the vibration intensity data a 1, a 118Myuesu, fixed fundamental frequency f 0 is 16.5 Hz, the division number k was 512.

図24は、第三実施形態の抽出手段3Cが実行する処理を示すフローチャートであり、抽出手段3は、車輪数を4として、まず、車輪番号rを1に設定し(ステップs301)、第1車輪について、車輪速データv1を入力し(ステップs102)、次いで、車輪速データv1に基づいて第1車輪の基本周波数を求めてサンプリングにおける可変時間間隔δv1を算出し(ステップs303)、そのあと、同期成分波形を初期化する(ステップs304)。 FIG. 24 is a flowchart showing the processing executed by the extraction means 3C of the third embodiment. The extraction means 3 first sets the wheel number r to 1 with the number of wheels set to 4 (step s301). For the wheels, wheel speed data v 1 is input (step s102), then the basic frequency of the first wheel is obtained based on the wheel speed data v 1 to calculate a variable time interval δ v1 in sampling (step s303), Thereafter, the synchronous component waveform is initialized (step s304).

次いで、モデル波形R3のi番目のデータR3(i)と計測波形X31のi番目のデータX31 (i)とから同期成分波形Y31のi番目のデータY31 (i)を算出するが、これは次のような手順に従う。すなわち、データ番号を初期化してi=1としたあと(ステップs305)、モデル波形のデータR3(1)をADF(1)に入力し(ステップs306)、一方、計測波形に関し、まず、振動強度データa1を取り込みながら(ステップs307)、これを時間間隔δv1でサンプリングして、計測波形X31の1番目のデータX31 (1)を生成し(ステップs308)、差分E31(1)を求め(ステップs309)これをADF(1)に入力する。 Then, calculate the i-th data Y 31 synchronization component waveform Y 31 from the i-th data X 31 of the i-th data R 3 of model waveform R 3 (i) and measuring the waveform X 31 (i) (i) This follows the following procedure. That is, after initializing the data number and setting i = 1 (step s305), the model waveform data R 3 (1) is input to the ADF (1) (step s306). While taking the intensity data a 1 (step s307), this is sampled at the time interval δ v1 to generate the first data X 31 (1) of the measurement waveform X 31 (step s308), and the difference E 31 (1 ) (Step s309) and this is input to ADF (1).

ADF(1)では、差分E31(1)によってフィルタリングパラメータを最適化処理するとともに、先に入力したR3(1)と、差分E31(1)とから、同期成分波形Y31の1番目のデータY31 (1)を求め(ステップs310)、これを出力して、次の差分E31(2)の算出に用いるとともに、異常警告装置4に出力する(ステップs311)。 In ADF (1), the filtering parameter is optimized by the difference E 31 (1), and the first input of the synchronous component waveform Y 31 from the previously input R 3 (1) and the difference E 31 (1). Data Y 31 (1) is obtained (step s310), and this is output and used to calculate the next difference E 31 (2) and is output to the abnormality warning device 4 (step s311).

以上で、Y31 (1)を求める一巡のステップを終了して、i= i +1として(ステップs312)、Y11 (2)を求めるステップにはいるが、その前にiがサンプリング数Mに達しているかを判定し(ステップs313)、もし達している場合には、第1輪についての同期成分波形Y31の抽出を完了し、r=r+1として(ステップs314)、第2輪について、同様の手順で同期成分波形Y32の抽出を行い、これを最終輪まで(ステップs315)繰り返し、一回の抽出工程を完了する。そして、走行中、この抽出工程を所定周期で繰り返し行う。 With the above, one round of steps for obtaining Y 31 (1) is completed, i = i + 1 is set (step s312), and Y 11 (2) is entered, but before that, i is the sampling number M. determines whether the reached (step s313), if they reached if completes the extraction of the sync component waveform Y 31 for the first wheel, as r = r + 1 (step s314), the second wheel for, was extracted synchronization component waveform Y 32 in the same procedure, this is repeated until the last wheel (step s315), to complete the single extraction step. During the traveling, this extraction process is repeated at a predetermined cycle.

なお、各車輪についての処理を、説明の都合上、時間的に直列したものとしたが、装置の処理時間が許容する限りにおいて、並列処理とすることができ、この場合、1回の抽出工程を短縮化することができる。   In addition, the processing for each wheel is assumed to be serially connected in time for convenience of explanation, but can be parallel processing as long as the processing time of the apparatus is allowed. Can be shortened.

以上説明したように、第一の実施形態と同様に、同期成分波形Y3rの抽出に際して、入力した振動強度データarを、フーリエ変換等の時間のかかる処理を行うことなく、時間間隔δvrでサンプリングしたものをそのまま、入力順に処理して同期成分波形Y3rを生成することができ、車輪に同期する成分の抽出における精度を犠牲にすることなく、抽出工程を短いものとすることができる。 As described above, as in the first embodiment, when the synchronous component waveform Y 3r is extracted, the input vibration intensity data a r is not subjected to time-consuming processing such as Fourier transformation, and the time interval δ vr Can be processed in the order of input to generate the synchronized component waveform Y 3r , and the extraction process can be shortened without sacrificing the accuracy in extracting the components synchronized with the wheels. .

次に、以上に説明した第三の実施形態の車輪異常検知装置を基本形として、抽出工程をさらに短縮化した第四の実施形態について説明する。この実施形態は、第三実施形態の抽出手段3Cを、抽出手段3Dに置換した構成になり、図25にブロック線図で示すように、抽出手段3Dは、振動計測手段1によって連続的に計測された各車輪についての振動強度データa1〜a4のそれぞれを、車輪速に応じて時間とともに変化する可変の時間間隔δvaでサンプリングして、時間順に並んだデータ列よりなる計測波形X41〜X44を生成する可変サンプリング処理部41と、各車輪について、計測波形X41〜X44の中から、その車輪の回転に同期した成分だけを、時間順に並んだデータ列として抽出し、これを、同期成分波形b1〜b4(もしくはY41〜Y44)として前記異常警告手段4に出力する適応デジタルフィルタ部42と、各車輪の前記抽出のモデルとなるモデル波形R4を、時間順に並んだデータ列として生成するモデル波形生成部43とよりなる。 Next, a fourth embodiment in which the extraction process is further shortened will be described based on the wheel abnormality detection device of the third embodiment described above as a basic form. In this embodiment, the extraction means 3C of the third embodiment is replaced with an extraction means 3D. The extraction means 3D is continuously measured by the vibration measurement means 1 as shown in a block diagram in FIG. Each of the measured vibration intensity data a 1 to a 4 for each wheel is sampled at a variable time interval δ va that changes with time according to the wheel speed, and a measurement waveform X 41 comprising a data sequence arranged in time order For the variable sampling processing unit 41 that generates ~ X 44 and for each wheel, only the components synchronized with the rotation of the wheel are extracted from the measurement waveforms X 41 to X 44 as a data sequence arranged in time order. Are output to the abnormality warning means 4 as synchronization component waveforms b 1 to b 4 (or Y 41 to Y 44 ), and a model waveform R 4 serving as a model for the extraction of each wheel is converted into time. De in order The more the model waveform generating unit 43 for generating a data sequence.

図18に示した、第三の実施形態の抽出手段3Aに対するブロック線図と、図25に示した、第四の実施形態の抽出手段3Dに対するブロック線図とを比較すると明らかなように、第三実施形態における可変サンプリング処理部31におけるサンプリングの時間間隔δvrが車輪ごとに異なるのに対して、第四の実施形態におけるサンプリングの時間間隔δaは、全車輪に対して共通の値が用いられる点が相違しており、(もちろん、δaは時間とともに変化する)このことによって、第四の実施形態のものは、詳細を後述するように、抽出処理に掛かる時間を一層の短縮化することができる。 As is apparent from a comparison between the block diagram for the extracting unit 3A of the third embodiment shown in FIG. 18 and the block diagram for the extracting unit 3D of the fourth embodiment shown in FIG. The sampling time interval δ vr in the variable sampling processing unit 31 in the third embodiment differs for each wheel, whereas the sampling time interval δ a in the fourth embodiment uses a common value for all wheels. that is has different, by (of course, [delta] a change with time) this, those of the fourth embodiment, as will be described in greater detail below, to further shorten the time required for the extraction process be able to.

可変サンプリング処理部31は、車輪90A〜90Dのそれぞれに対応する処理部(1)〜(4)よりなり、例えば処理部(1)では、図26に横軸に時間をとって示した、振動強度データ波形a1を、時間とともに変化する車輪速vaに応じて可変の時間間隔δvaでサンプリングし、サンプリングされたデータを、図27に、横軸を時間軸に沿うデータ順としたグラフで示すように、時間順に並べたデータ列よりなる計測波形X41を出力する。図27において、X41(i-1)、X41(i)、X41(i+1)は、それぞれ、時間順に連続して並んでデータ列を構成する、i-1、i、i+1番目のデータを示す。 The variable sampling processing unit 31 includes processing units (1) to (4) corresponding to the wheels 90A to 90D. For example, in the processing unit (1), the vibration shown in FIG. The intensity data waveform a 1 is sampled at a variable time interval δ va according to the wheel speed v a that changes with time, and the sampled data is shown in FIG. 27, with the horizontal axis representing the data order along the time axis. As shown, the measurement waveform X 41 composed of a data sequence arranged in time order is output. In FIG. 27, X 41 (i−1), X 41 (i), and X 41 (i + 1) are respectively arranged in a time sequence continuously to form a data string, i−1, i, i + The first data is shown.

上記時間間隔δvaの算出の根拠となる車輪速va、あるいはこれに対応する基本周波数faは、他の車輪の車輪速、あるいは、これ以外の、それぞれの車輪の車輪速に応じて変化する一つの時間の関数とすることができ、例えば、これを車の速度から計算してもいいし、全車輪の車輪速の平均値に対する基本周波数としてもよい。そして、このときの時間間隔δvaは、このようにして設定された基本周波数faと、車輪1回転を等分する際の分割数となる予め定められた値kとの積の逆数として求められる。 The wheel speed v a that is the basis for calculating the time interval δ va , or the fundamental frequency f a corresponding thereto, changes according to the wheel speeds of other wheels or other wheel speeds. For example, this may be calculated from the speed of the car, or may be a fundamental frequency with respect to the average value of the wheel speeds of all the wheels. Then, the time interval [delta] va at this time is determined as the reciprocal of the product of the fundamental frequency f a which is set this way, the predetermined value k as the number of divisions when equally dividing the wheel 1 rotating It is done.

ここで、予め定められた固定の周波数を固定基本周波数f0として、この固定基本周波数に自然数nを乗じた周波数を、「n次の固定同期周波数」とし、前記固定基本周波数f0と前記分割数kとの積の逆数を時間間隔δ0とする。 Here, a predetermined fixed frequency is defined as a fixed fundamental frequency f 0 , and a frequency obtained by multiplying the fixed fundamental frequency by a natural number n is referred to as an “nth-order fixed synchronization frequency”, and the fixed fundamental frequency f 0 and the division Let the reciprocal of the product of the number k be the time interval δ 0 .

図28は、同期成分波形Y4rの抽出に用いるモデル波形R4の各データR4(i)を、予め定められた前記時間間隔δ0おきに並べた時間領域波形をフーリエ変換してできる周波数領域波形を示し、この周波数領域波形は、前記抽出に係る次数nの固定同期周波数f0 (n)=n・f0の位置で山となり、隣接する次数の固定同期周波数同士の間では谷となる形状をなし、前記山となる部分の頂部は無限小もしくは有限の幅の周波数範囲で点もしくは平坦面として延在し、頂部の高さは、前記抽出に係るすべての次数の固定同期周波数に対して同じ値となる特性を有するよう構成され、この点においては、第三の実施形態に用いられるモデル波形R3と同様である。 FIG. 28 shows a frequency obtained by Fourier-transforming a time domain waveform in which each data R 4 (i) of the model waveform R 4 used for extraction of the synchronous component waveform Y 4r is arranged every predetermined time interval δ 0. This frequency domain waveform has a peak at the position of the fixed synchronization frequency f 0 (n) = n · f 0 of the order n related to the extraction, and a trough between the fixed synchronization frequencies of the adjacent orders. The top of the peak portion extends as a point or flat surface in a frequency range of infinitesimal or finite width, and the height of the top corresponds to the fixed synchronization frequency of all orders related to the extraction. In contrast, this is the same as the model waveform R 3 used in the third embodiment.

しかし、モデル波形R4が波形R3と異なる点は、頂部が点ではなく、有限の幅を有する平坦面として延在する点であり、モデル波形R4は、車輪速の異なるどの車輪に対しても同期周波数成分を抽出することができる打ち抜き型として機能させるための要件を満足する必要があったからである。 However, that model waveform R 4 is different from the waveform R 3 is not a top a point, a point extending as a flat surface having a finite width, the model waveform R 4, compared which wheels of different wheel speeds This is because it is necessary to satisfy the requirements for functioning as a punching die that can extract the synchronization frequency component.

すなわち、全車輪のうち、最高車速の車輪に対する基本周波数(周波数は最小)をfminとし、最低車速の車輪に対する基本周波数(周波数は最大)をfmaxとしたとき、このような要件を満足するモデル波形Rのn次の頂部の延在領域の上限fu (n)および下限fw (n)は、それぞれ、式(2)および(3)の条件を満足することが必要である。

Figure 2006224757
That is, among all the wheels, such a requirement is satisfied when the basic frequency (frequency is minimum) for the wheel at the highest vehicle speed is f min and the basic frequency (frequency is maximum) for the wheel at the lowest vehicle speed is f max. model waveform upper limit of the extension region of the n-th top of the R 4 f u (n) and the lower limit f w (n), respectively, it is necessary to satisfy the condition of formula (2) and (3).
Figure 2006224757

一方、周波数領域波形が、上記のような要件を満足するように、モデル波形R4を生成するには、式(7)および(8)に基づいて、既知のfo、fmin、fmax、faから、n次の頂部の延在領域をそれぞれ設定し、この設定領域を適当な個数J、例えば20個の周波数の集まりからなるものと近似して、これらの周波数に対する正弦波を足し合わせて、式(9)に示す時間領域波形R4(t)を生成し、この時間領域波形を前記固定の時間間隔δ0でサンプリングして離散データを作り、モデル波形とすればよい。なお、式(5)における、fw(1)、およびfu(1)は、1次の頂部の延在領域に対する下限、および上限を表し、A4は予め定められた定数である。

Figure 2006224757
On the other hand, in order to generate the model waveform R 4 so that the frequency domain waveform satisfies the above-described requirements, a known f o , f min , f max is generated based on the equations (7) and (8). , F a , an extension region of the top of the n-th order is set, and this setting region is approximated to an appropriate number J, for example, a set of 20 frequencies, and sine waves for these frequencies are added. In addition, a time domain waveform R 4 (t) shown in Equation (9) is generated, and this time domain waveform is sampled at the fixed time interval δ 0 to create discrete data, which is used as a model waveform. In Equation (5), fw (1) and fu (1) represent a lower limit and an upper limit for the extension region of the primary apex, and A 4 is a predetermined constant.
Figure 2006224757

モデル波形R4を周波数領域に変換してできる波形の頂部の延在領域は、第二の実施形態について説明したのと同様に、最高速の車輪の基本周波数fminと、最低速の車輪の基本周波数fmaxとの両方に基づいて決めることができるほか、より簡便な手段として、この代わりに、Aを定数として、fminの0〜+A%、もしくは、fmaxの-A〜0%とするもこともでき、さらに、車体ヨーレート等の車輪速に相関する特性も用いることができる。 The extension region at the top of the waveform formed by converting the model waveform R 4 into the frequency domain is similar to that described in the second embodiment, and the basic frequency f min of the highest speed wheel and the lowest speed wheel. It can be determined based on both the fundamental frequency f max and, as a simpler means, instead of A as a constant, 0 to + A% of f min or -A to 0% of f max Furthermore, characteristics correlating with the wheel speed such as the vehicle body yaw rate can also be used.

上記のように構成された抽出手段3Dの、車輪の回転に同期する成分の抽出の成否を確認するため、走行中の第一車輪90Aの振動強度データa1を、第二車輪90Bに対応する基本周波数f2(1)に基づく時間間隔δvaでサンプリングして生成された計測波形X41、および、デジタルフィルタ部42から出力された同期成分波形Y41について、これらを構成する各データX41(i)、およびY41(i)を、時間間隔δ0で並べて時間領域波形を生成し、さらにこれらの時間領域波形をフーリエ変換して周波数領域波形を算出した。図29は、計測波形X41に対応する周波数領域波形を示し、図30は、同期成分波形Y21に対応する周波数領域波形を示す。 The configured extractor 3D as described above, to confirm the success or failure of extraction of components synchronizing the rotation of the wheel, the vibration intensity data a 1 of the first wheel 90A during running, corresponding to the second wheel 90B With respect to the measurement waveform X 41 generated by sampling at the time interval δ va based on the fundamental frequency f2 (1) and the synchronous component waveform Y 41 output from the digital filter unit 42, each data X 41 ( A time domain waveform is generated by arranging i) and Y 41 (i) at a time interval δ 0 , and further, a frequency domain waveform is calculated by Fourier transforming these time domain waveforms. Figure 29 shows the frequency domain waveform corresponding to the measured waveform X 41, Fig. 30 shows a frequency domain waveform corresponding to the synchronizing component waveform Y 21.

図29と図30とを比較すると、図30は、同期周波数f1 (1)〜f1 (20)において、図29における振動強度と同じ値を有する高さの山を有するとともに、これらの間には谷部が形成された形の波形をなしており、第一車輪90Aの計測波形X41を生成するのに、第二車輪90Bの車輪速を基にした時間間隔を用いてサンプリングしたにもかかわらず、所期の抽出が行えたことを示している。 When FIG. 29 is compared with FIG. 30, FIG. 30 has a height crest having the same value as the vibration intensity in FIG. 29 at the synchronization frequencies f 1 (1) to f 1 (20) . and forms a waveform in the form of troughs are formed on, in to generate a measurement waveform X 41 of the first wheel 90A, and sampled using a time interval based on a wheel speed of the second wheel 90B Nevertheless, it shows that the expected extraction was performed.

なお、上記算出に用いた、振動強度データa1のサンプリングにおける時間間隔δv2は、
118μsとし、固定基本周波数f0は16.5Hz、分割数kは512であった。
The time interval δ v2 in the sampling of the vibration intensity data a 1 used for the above calculation is
The fixed fundamental frequency f 0 was 16.5 Hz and the division number k was 512.

図11は、第四実施形態の抽出手段3Dが実行する処理の流れを示すフローチャートであり、抽出手段3Dは、車輪数を4として、まず、モデル波形生成部43に各車輪の車輪速データv1〜v4を入力し(ステップs401)、これに基づいて先に説明した手順に従ってモデル波形R4(i) (i=1〜M、ただし、Mは1回の抽出に用いるサンプリング数)を生成する(ステップs402)。 FIG. 11 is a flowchart showing the flow of processing executed by the extraction means 3D of the fourth embodiment. The extraction means 3D first sets the number of wheels to 4, and first sends the wheel speed data v of each wheel to the model waveform generation unit 43. 1 to v 4 are input (step s401), and the model waveform R 4 (i) (i = 1 to M, where M is the number of samplings used for one extraction) is input in accordance with the procedure described above based on this. Generate (step s402).

そして、同期成分波形Y4rのデータを一つずつデータ順に求める手順にはいるが、これに先立って、各車輪に対する同期成分波形を初期化する(ステップs403)とともに、前記車輪速v1〜v4(あるいは、この代用特性)に基づいて可変サンプリングの時間間隔δvaを算出しておく(ステップs404)。次いで、i=1としたあと(ステップs405)、モデル波形生成部43で生成したモデル波形R4の1番目のデータR4(i)をADF(1)に入力する(ステップs406)。 Then, although the procedure for obtaining the data of the synchronous component waveform Y 4r one by one in the order of data is performed, prior to this, the synchronous component waveform for each wheel is initialized (step s403) and the wheel speeds v 1 to v Based on 4 (or this substitute characteristic), the variable sampling time interval δ va is calculated (step s404). Next, after setting i = 1 (step s405), the first data R 4 (i) of the model waveform R 4 generated by the model waveform generation unit 43 is input to the ADF (1) (step s406).

次に、このモデル波形データR4(1)を用いて、1〜4輪のすべての車輪に対する同期成分波形の1番目のデータY4r(1)(r=1〜4)を算出するが、これを第1輪から始めるためr=1として(ステップs407)、固定サンプリング処理部41の処理部(1)に、振動強度データa1を入力し(ステップs408)、次いで、これを時間間隔δvaでサンプリングして計測波形データX41(i)を生成し(ステップs409)、差分E41(1)を求め(ステップs410)これをADF(1)に入力する。 Next, using this model waveform data R 4 (1), the first data Y 4r (1) (r = 1 to 4) of the synchronous component waveform for all of the 1 to 4 wheels is calculated. as r = 1 to start it from the first wheel (step s407), the processing unit of the fixed sampling processing section 41 (1), enter the vibration intensity data a 1 (step s408), then this time interval δ Sampling with va generates measurement waveform data X 41 (i) (step s409), and obtains a difference E 41 (1) (step s410), which is input to ADF (1).

ADF(1)では、差分E41(1)によってフィルタリングパラメータを最適化処理するとともに、先に入力したR41(1)と、差分E41(1)とから、同期成分波形Y41の1番目のデータY41(1)を求め(ステップs411)、次の差分E41(2)の算出に用いるとともに、異常警告装置4に出力する(ステップs412)。 In ADF (1), the filtering parameter is optimized by the difference E 41 (1), and the first input of the synchronous component waveform Y 41 from the previously input R 41 (1) and the difference E 41 (1). Data Y 41 (1) is obtained (step s411), is used to calculate the next difference E 41 (2), and is output to the abnormality warning device 4 (step s412).

以上で、Y41 (1)を求める一巡のステップを終了し、続いて同様のステップを第2輪以降についても実施するため、r=r+1として(ステップs413)、r=5になるまでこれを繰り返す(ステップs414)。i=1に対する上記の処理を完了すると、i=i+1として(ステップs415)、i=M(ステップs416)になるまでこれを繰り返し、一回の抽出工程を完了する。 In this way, the round of steps for obtaining Y 41 (1) is completed, and then the same steps are performed for the second and subsequent wheels, so r = r + 1 (step s413) until r = 5 This is repeated (step s414). When the above processing for i = 1 is completed, i = i + 1 is set (step s415), and this is repeated until i = M (step s416), thereby completing one extraction step.

以上説明したように、同期成分波形Y4rの算出に際して、入力した振動強度データarを、フーリエ変換等の時間のかかる処理を行うことなく、所定の時間間隔δvaでサンプリングしたものをそのまま、入力順に処理して同期成分波形Y4rを生成して、短時間で抽出工程を終えることができる点においては、第三の実施形態と同様であるが、この第四の実施形態の装置では、以下の点において、第三の実施形態のものよりも、処理の時間を短くすることができる。 As described above, when calculating the synchronous component waveform Y 4r , the input vibration intensity data a r is sampled at a predetermined time interval δ va without performing time-consuming processing such as Fourier transform, It is the same as the third embodiment in that it can be processed in the input order to generate the synchronous component waveform Y 4r and the extraction process can be completed in a short time, but in the apparatus of the fourth embodiment, In the following points, the processing time can be made shorter than that of the third embodiment.

これは、ステップs409における、計測波形データX4r(i)を生成するに際し、各車輪に共通な時間間隔δvaを用いることによるものであり、このことによる効果の第一は、第三実施形態においては、サンプリングの時間間隔δvrが車輪によって異なるので、i番目の同期成分波形データY4r(i)を求めるのに、これを車輪同士で同期させることができず、それぞれの車輪についての処理を独立して行わざるを得ず、そのため、各車輪について各同期成分波形データY41(i)を計算する都度、モデル波形データR1r(i)をメモリあるいはバッファから呼び出さなければならないのに対して、本実施形態においては、1回だけR4(i)をメモリもしくはバッファから呼び出すだけで、全車輪のi番目のデータの処理をすべて済ませることができる点である。すなわち、Mをサンプリング数、Rを車輪数として、第三の実施形態のおいては、全部で、M・R回モデル波形を呼び出す必要があったのに対し、第四の実施形態のおいては、これをM回で済ませることができる。 This is due to the use of the time interval δ va common to each wheel when generating the measurement waveform data X 4r (i) in step s409. The first effect of this is the third embodiment. Since the sampling time interval δ vr differs depending on the wheel, the i-th synchronization component waveform data Y4r (i) cannot be synchronized between the wheels to obtain the i-th synchronization component waveform data Y4r (i). For this reason, the model waveform data R 1r (i) must be called from the memory or buffer each time the synchronous component waveform data Y 41 (i) is calculated for each wheel. In this embodiment, the processing of the i-th data of all the wheels can be completed by calling R 4 (i) from the memory or buffer only once. That is, with M as the number of samplings and R as the number of wheels, in the third embodiment, it was necessary to call the M / R model waveform in all, whereas in the fourth embodiment Can do this M times.

時間間隔δvaを全車輪について共通のものとしたことによる第二の効果は、ステップs404における、時間間隔δvaの生成に要する処理回数を減らせる点であり、第三の実施形態においては、各車輪に対して別個の時間間隔δvrを生成しなければならなかったところ、第四の実施形態においては、共通の時間間隔δvaを一個算出するだけで済ませることができ、このことのよっても、抽出工程の短縮に寄与させることができる。 The second effect of having the time interval δ va common to all the wheels is that the number of processes required to generate the time interval δ va in step s404 can be reduced. In the third embodiment, Since a separate time interval δ vr had to be generated for each wheel, in the fourth embodiment, it is only necessary to calculate one common time interval δ va , Can also contribute to shortening the extraction process.

本発明の車輪異常検知装置は、種々の車両における車輪の異常を検出するのに用いることができる。   The wheel abnormality detection device of the present invention can be used to detect wheel abnormality in various vehicles.

車輪異常検知装置の配置を示す平面模式図、および、車輪の一部に生じた異常部分を示す車輪の側面模式図である。It is the plane schematic diagram which shows arrangement | positioning of a wheel abnormality detection apparatus, and the side surface schematic diagram which shows the abnormal part which arose in a part of wheel. 車輪異常検知装置の構成を示すブロック線図である。It is a block diagram which shows the structure of a wheel abnormality detection apparatus. 第一の実施形態の車両異常検出装置の抽出手段を示すブロック線図である。It is a block diagram which shows the extraction means of the vehicle abnormality detection apparatus of 1st embodiment. 第一の実施形態に係る振動強度データ波形a1を示すグラフである。It is a graph which shows vibration intensity data waveform a1 concerning a first embodiment. 第一の実施形態に係る計測波形X11の一部分を示すグラフである。It is a graph showing a portion of a measurement waveform X 11 according to the first embodiment. 第一の実施形態に係る同期成分波形Y11の一部分を示すグラフである。It is a graph showing a portion of a synchronous component waveform Y 11 according to the first embodiment. 第一の実施形態に係るモデル波形R11の一部分を示すグラフである。It is a graph showing a portion of the model waveform R 11 according to the first embodiment. 第一の実施形態に係るモデル波形R11に対応する周波数領域波形を示すグラフである。It is a graph showing a frequency domain waveform corresponding to the model waveform R 11 according to the first embodiment. 第一の実施形態に係る計測波形X11に対応する周波数波形を示すグラフである。It is a graph showing a frequency waveform corresponding to the measured waveform X 11 according to the first embodiment. 第一の実施形態に係る同期成分波形Y11に対応する周波数波形を示すグラフである。It is a graph showing a frequency waveform corresponding to the synchronizing component waveform Y 11 according to the first embodiment. 第一の実施形態に係る抽出手段の処理を示すフローチャートである。It is a flowchart which shows the process of the extraction means which concerns on 1st embodiment. 第二の実施形態の車両異常検出装置の抽出手段を示すブロック線図である。It is a block diagram which shows the extraction means of the vehicle abnormality detection apparatus of 2nd embodiment. 第二の実施形態に係るモデル波形R21の一部分を示すグラフである。It is a graph showing a portion of the model waveform R 21 according to the second embodiment. 第二の実施形態に係るモデル波形R21に対応する周波数領域波形を示すグラフである。It is a graph showing a frequency domain waveform corresponding to the model waveform R 21 according to the second embodiment. 第二の実施形態に係る計測波形X21に対応する周波数波形を示すグラフである。It is a graph showing a frequency waveform corresponding to the measured waveform X 21 according to the second embodiment. 第二の実施形態に係る同期成分波形Y21に対応する周波数波形を示すグラフである。It is a graph showing a frequency waveform corresponding to the second according to the embodiment synchronizing component waveform Y 21. 第二の実施形態に係る抽出手段の処理を示すフローチャートである。It is a flowchart which shows the process of the extraction means which concerns on 2nd embodiment. 第三の実施形態の車両異常検出装置の抽出手段を示すブロック線図である。It is a block diagram which shows the extraction means of the vehicle abnormality detection apparatus of 3rd embodiment. 第三の実施形態に係る振動強度データ波形a1を示すグラフである。Is a graph showing the vibration intensity data waveform a 1 according to the third embodiment. 第三の実施形態に係る計測波形X31の一部分を示すグラフである。It is a graph showing a portion of a measurement waveform X 31 according to the third embodiment. 第三の実施形態に係るモデル波形R31に対応する周波数領域波形を示すグラフである。It is a graph showing a frequency domain waveform corresponding to the model waveform R 31 according to the third embodiment. 第三の実施形態に係る計測波形X31に対応する周波数波形を示すグラフである。It is a graph showing a frequency waveform corresponding to the measured waveform X 31 according to the third embodiment. 第三の実施形態に係る同期成分波形Y31に対応する周波数波形を示すグラフである。It is a graph showing a frequency waveform corresponding to the synchronizing component waveform Y 31 according to the third embodiment. 第三の実施形態に係る抽出手段の処理を示すフローチャートである。It is a flowchart which shows the process of the extraction means which concerns on 3rd embodiment. 第四の実施形態の車両異常検出装置の抽出手段を示すブロック線図である。It is a block diagram which shows the extraction means of the vehicle abnormality detection apparatus of 4th embodiment. 第四の実施形態に係る振動強度データ波形a1を示すグラフである。Is a graph showing the vibration intensity data waveform a 1 according to the fourth embodiment. 第四の実施形態に係る計測波形X41の一部分を示すグラフである。It is a graph showing a portion of a measurement waveform X 41 according to the fourth embodiment. 第四の実施形態に係るモデル波形R41に対応する周波数領域波形を示すグラフである。It is a graph showing a frequency domain waveform corresponding to the model waveform R 41 according to the fourth embodiment. 第四の実施形態に係る計測波形X41に対応する周波数波形を示すグラフである。It is a graph showing a frequency waveform corresponding to the measured waveform X 41 according to the fourth embodiment. 第四の実施形態に係る同期成分波形Y41に対応する周波数波形を示すグラフである。It is a graph showing a frequency waveform corresponding to the synchronizing component waveform Y 41 according to the fourth embodiment. 第四の実施形態に係る抽出手段の処理を示すフローチャートである。It is a flowchart which shows the process of the extraction means which concerns on 4th embodiment.

符号の説明Explanation of symbols

1 振動計測手段
1a〜1d 振動センサ
2 車輪速計測手段
2a〜2d 車輪速センサ
3、3A、3B、3C、3D 抽出手段
4 異常警告手段
4a 異常制御部
4b 異常表示部
5 中央制御装置
10 車輪異常検知装置
11 固定サンプリング処理部
12 適応デジタルフィルタ部
13 モデル波形生成部
21 固定サンプリング処理部
22 適応デジタルフィルタ部
23 モデル波形生成部
31 可変サンプリング処理部
32 適応デジタルフィルタ部
33 モデル波形生成部
41 可変サンプリング処理部
42 適応デジタルフィルタ部
43 モデル波形生成部
90a〜90d 車輪
91 層間剥離
DESCRIPTION OF SYMBOLS 1 Vibration measuring means 1a-1d Vibration sensor 2 Wheel speed measuring means 2a-2d Wheel speed sensor 3, 3A, 3B, 3C, 3D Extraction means 4 Abnormal warning means 4a Abnormal control part 4b Abnormal display part 5 Central controller 10 Wheel abnormality Detection apparatus 11 Fixed sampling processing unit 12 Adaptive digital filter unit 13 Model waveform generation unit 21 Fixed sampling processing unit 22 Adaptive digital filter unit 23 Model waveform generation unit 31 Variable sampling processing unit 32 Adaptive digital filter unit 33 Model waveform generation unit 41 Variable sampling Processing unit 42 Adaptive digital filter unit 43 Model waveform generation unit 90a to 90d Wheel 91 Delamination

Claims (10)

走行中の車両に装着された複数の車輪に対して、各車輪の振動強度を連続的に計測する振動計測手段と、各車輪の回転速度を、直接的に、もしくは、車輪速と相関する車両特性を測定することによって間接的に、連続的して計測する車輪速計測手段と、前記振動計測手段および車輪速計測手段で計測されたデータから車輪の回転に同期した波形を抽出する抽出手段と、抽出手段で抽出された波形から車輪の状態を判定しこの状態が予め定めた正常範囲から外れた場合に異常を警告する異常警告手段とを具えた車輪異常検知装置において、
前記抽出手段は、
振動計測手段によって連続的に計測された各車輪についての振動強度データを、予め定められた固定の時間間隔δfでサンプリングして、時間順に並んだデータ列よりなる「計測波形」を生成する固定サンプリング処理部と、
各車輪について、前記計測波形の中から、その車輪の回転に同期した成分だけを、時間順に並んだデータ列として抽出し、これを、「同期成分波形」として前記異常警告手段に出力する適応デジタルフィルタ部と、
各車輪の前記抽出のモデルとなる「モデル波形」を、時間順に並んだデータ列として生成するモデル波形生成部とよりなり、
前記適応デジタルフィルタ部は、前記モデル波形を入力し、前記計測波形と前記同期成分波形との差分、および、入力したモデル波形の両方に基づいて前記同期成分波形を生成するとともに、前記同期成分波形を生成する際のフィルタリングパラメータを、前記差分を用いて適応制御するよう構成され、
各車輪の一回転分の周期に対応する基本周波数に自然数nを乗じた周波数を、その車輪に対する「n次の同期周波数」としたとき、
各車輪に対する前記モデル波形は、このモデル波形の各データを前記δfの時間間隔で並べた時間領域波形をフーリエ変換してできる周波数領域波形が、その車輪の、前記抽出に係る次数の同期周波数の位置で山となり、隣接する次数の同期周波数同士の間では谷となる形状をなし、前記山となる部分の頂部は無限小もしくは有限の幅の周波数範囲に点もしくは平坦面として延在し、頂部の高さは、前記抽出に係るすべての次数の同期周波数に対して同じ値となる特性を有するよう構成され、
モデル波形生成部は、前記車輪速計測手段からのデータにもとづいて、各車輪の基本周波数、もしくは、この基本周波数が含まれる範囲を設定し、設定した基本周波数もしくはその範囲から前記周波数領域波形特性を満足するモデル波形を生成するよう構成されてなる車輪異常検知装置。
For a plurality of wheels mounted on a running vehicle, a vibration measuring means for continuously measuring the vibration intensity of each wheel, and a vehicle that correlates the rotational speed of each wheel directly or with the wheel speed. Wheel speed measuring means for measuring continuously and indirectly by measuring characteristics; and extracting means for extracting a waveform synchronized with the rotation of the wheel from the data measured by the vibration measuring means and the wheel speed measuring means; In the wheel abnormality detection device comprising the abnormality warning means for determining the state of the wheel from the waveform extracted by the extraction means and warning the abnormality when this state is outside the predetermined normal range,
The extraction means includes
Fixed to generate a “measurement waveform” consisting of a data sequence arranged in time order by sampling vibration intensity data for each wheel continuously measured by the vibration measuring means at a predetermined fixed time interval δ f A sampling processing unit;
For each wheel, only the component synchronized with the rotation of the wheel is extracted from the measurement waveform as a data sequence arranged in time order, and this is output as the “synchronous component waveform” to the abnormality warning means. A filter section;
It consists of a model waveform generation unit that generates a “model waveform” as a model of the extraction of each wheel as a data sequence arranged in time order,
The adaptive digital filter unit receives the model waveform, generates the synchronization component waveform based on both the difference between the measurement waveform and the synchronization component waveform, and the input model waveform, and the synchronization component waveform Is configured to adaptively control the filtering parameter in generating the difference using the difference,
When the frequency obtained by multiplying the fundamental frequency corresponding to the period of one rotation of each wheel by a natural number n is the `` nth synchronization frequency '' for that wheel,
The model waveform for each wheel, frequency domain waveform of time domain waveform obtained by arranging the data in the time interval of the [delta] f can be the Fourier transform of this model waveform, the wheel, the next number of the synchronization frequency according to the extracted It becomes a mountain at the position of, and forms a valley between adjacent synchronization frequencies of the order, the top of the portion that becomes the mountain extends to a frequency range of infinitesimal or finite width as a point or flat surface, The height of the top is configured to have the same value for all orders of synchronization frequency for the extraction,
The model waveform generation unit sets a basic frequency of each wheel or a range including the basic frequency based on data from the wheel speed measuring unit, and sets the frequency domain waveform characteristics from the set basic frequency or the range. A wheel abnormality detection device configured to generate a model waveform that satisfies the above.
前記適応デジタルフィルタ部は、複数の車輪に対するそれぞれの同期成分波形を、これらに共通の、一つのモデル波形に基づいて生成するよう構成され、
この共通のモデル波形に対応する前記周波数領域波形において、n次の同期周波数に対応する山の頂部の周波数範囲は、前記共通のモデル波形に基づく同期成分波形の抽出に係る全ての車輪に対するn次の同期周波数を含むことを特徴とする請求項1に記載の車輪異常検知装置。
The adaptive digital filter unit is configured to generate respective synchronous component waveforms for a plurality of wheels based on a single model waveform common to them,
In the frequency domain waveform corresponding to the common model waveform, the frequency range of the peak corresponding to the n-th synchronization frequency is the n-th order for all wheels related to the extraction of the synchronization component waveform based on the common model waveform. The wheel abnormality detection device according to claim 1, further comprising:
Aを予め定められた百分率として、前記、n次の同期周波数に対応する山の頂部の周波数範囲を、前記共通のモデル波形に基づく同期成分波形の抽出に係る複数の車輪のうち、基本周波数が最小となる車輪のn次同期周波数の0〜+A%として設定してなる請求項2に記載の車輪異常検知装置。   Assuming that A is a predetermined percentage, the frequency range at the top of the peak corresponding to the n-th synchronization frequency is the basic frequency among the plurality of wheels related to the extraction of the synchronization component waveform based on the common model waveform. The wheel abnormality detection device according to claim 2, wherein the wheel abnormality detection device is set as 0 to + A% of the n-th order synchronization frequency of the minimum wheel. 前記、n次の同期周波数に対応する山の頂部の周波数範囲を、前記共通のモデル波形に基づく同期成分波形の抽出に係る複数の車輪のうち、基本周波数が最小となる車輪と、基本周波数が最大となる車輪との、両方の車輪に対するn次同期周波数に基づいて設定してなる請求項2に記載の車輪異常検知装置。   The frequency range at the top of the mountain corresponding to the n-th synchronization frequency is a wheel having a minimum fundamental frequency among a plurality of wheels related to extraction of a synchronization component waveform based on the common model waveform, and the fundamental frequency is The wheel abnormality detection device according to claim 2, wherein the wheel abnormality detection device is set based on the n-th synchronization frequency for both wheels with the largest wheel. 前記前記適応デジタルフィルタ部は、前記共通のモデル波形を構成するデータを、一つ入力するごとに、この共通のモデル波形に基づく同期成分波形の抽出に係る全車輪のそれぞれについて、同期成分波形のデータを一つだけ出力する処理を行うよう構成されてなる請求項2〜4のいずれかに記載の車輪異常検知装置。   Each time the adaptive digital filter unit inputs data constituting the common model waveform, for each of all wheels related to extraction of the synchronous component waveform based on the common model waveform, The wheel abnormality detection device according to any one of claims 2 to 4, wherein the wheel abnormality detection device is configured to perform a process of outputting only one data. 走行中の車両に装着された複数の車輪に対して、各車輪の振動強度を連続的に計測する振動計測手段と、各車輪の回転速度を、直接的に、もしくは、車輪速と相関する車両特性を測定することによって間接的に、連続的して計測する車輪速計測手段と、前記振動計測手段および車輪速計測手段で計測されたデータから車輪の回転に同期した波形を抽出する抽出手段と、抽出手段で抽出された波形から車輪の状態を判定しこの状態が予め定めた正常範囲から外れた場合に異常を警告する異常警告手段とを具えた車輪異常検知装置において、
前記抽出手段は、
各車輪について、振動計測手段により時間軸に沿って連続的に計測された振動強度データを、車輪速計測手段からのデータに応じて変化する時間間隔δvでサンプリングして、時間順に並んだデータ列よりなる「計測波形」を生成する可変サンプリング処理部と、
各車輪について、前記計測波形の中から、その車輪の回転に同期した成分だけを、時間順に並んだデータ列として抽出し、これを、「同期成分波形」として前記異常警告手段に出力する適応デジタルフィルタ部と、
各車輪の前記抽出のモデルとなる「モデル波形」を、時間順に並んだデータ列として生成するモデル波形生成部とよりなり、
前記可変サンプリング処理部において各車輪に対する前記計測波形を生成するのに用いられる前記時間間隔δvは、前記車輪速計測手段からのデータにもとづいて算出されるこの車輪の1回転分の周期に対応する基本周波数、もしくは、前記車輪速計測手段からのデータにもとづいて予め定められた手順で算出される仮想車輪の1回転分の周期に対応する基本周波数faと、これらの車輪の1回転を等分する際の分割数となる予め定められた値kとの積の逆数として求められ、
前記適応デジタルフィルタ部は、前記モデル波形を入力し、前記計測波形と前記同期成分波形との差分、および、入力したモデル波形の両方に基づいて前記同期成分波形を生成するとともに、前記同期成分波形を生成する際のフィルタリングパラメータを、前記差分を用いて適応制御し、複数の車輪に対するそれぞれの同期成分波形は、これらの車輪に共通な一つのモデル波形に基づいて生成されるよう構成され、
予め定められた固定の周波数を固定基本周波数f0として、この固定基本周波数に自然数nを乗じた周波数を、「n次の固定同期周波数」とし、前記固定基本周波数f0と前記分割数kとの積の逆数を時間間隔δ0としたとき、
前記共通のモデル波形は、このモデル波形の各データを前記時間間隔δ0おきに並べた時間領域波形をフーリエ変換してできる周波数領域波形が、前記抽出に係る次数の固定同期周波数の位置で山となり、隣接する次数の固定同期周波数同士の間では谷となる形状をなし、前記山となる部分の頂部は無限小もしくは有限の幅の周波数範囲で点もしくは平坦面として延在し、頂部の高さは、前記抽出に係るすべての次数の固定同期周波数に対して同じ値となる特性を有するよう構成されてなる車輪異常検知装置。
For a plurality of wheels mounted on a running vehicle, a vibration measuring means for continuously measuring the vibration intensity of each wheel, and a vehicle that correlates the rotational speed of each wheel directly or with the wheel speed. Wheel speed measuring means for measuring continuously and indirectly by measuring characteristics; and extracting means for extracting a waveform synchronized with the rotation of the wheel from the data measured by the vibration measuring means and the wheel speed measuring means; In the wheel abnormality detection device comprising the abnormality warning means for determining the state of the wheel from the waveform extracted by the extraction means and warning the abnormality when this state is outside the predetermined normal range,
The extraction means includes
For each wheel, the vibration intensity data continuously measured along the time axis by the vibration measuring means is sampled at a time interval δ v that changes according to the data from the wheel speed measuring means, and the data is arranged in time order. A variable sampling processing unit that generates a "measurement waveform" consisting of columns;
For each wheel, only the component synchronized with the rotation of the wheel is extracted from the measurement waveform as a data sequence arranged in time order, and this is output as the “synchronous component waveform” to the abnormality warning means. A filter section;
It consists of a model waveform generation unit that generates a “model waveform” as a model of the extraction of each wheel as a data sequence arranged in time order,
The time interval δ v used to generate the measurement waveform for each wheel in the variable sampling processing unit corresponds to a cycle of one rotation of the wheel calculated based on data from the wheel speed measuring means. Or a fundamental frequency f a corresponding to a cycle of one rotation of the virtual wheel calculated by a predetermined procedure based on data from the wheel speed measuring means, and one rotation of these wheels. It is obtained as the reciprocal of the product with a predetermined value k which is the number of divisions when equally dividing,
The adaptive digital filter unit receives the model waveform, generates the synchronization component waveform based on both the difference between the measurement waveform and the synchronization component waveform, and the input model waveform, and the synchronization component waveform The filtering parameter when generating is adaptively controlled using the difference, and each synchronous component waveform for a plurality of wheels is configured to be generated based on one model waveform common to these wheels,
A predetermined fixed frequency is set as a fixed basic frequency f 0 , and a frequency obtained by multiplying the fixed basic frequency by a natural number n is referred to as an “n-th fixed synchronization frequency”, and the fixed basic frequency f 0 and the division number k are When the reciprocal of the product of the time interval δ 0 ,
The common model waveform is a frequency domain waveform obtained by performing Fourier transform on a time domain waveform in which each data of the model waveform is arranged at every time interval δ 0 , and is a peak at the position of the fixed synchronization frequency of the order related to the extraction. Between the fixed synchronization frequencies of adjacent orders, forming a valley, and the top of the peak is extended as a point or flat surface in a frequency range of infinitesimal or finite width, and the height of the top The wheel abnormality detection device is configured to have the same value with respect to the fixed synchronization frequencies of all orders related to the extraction.
各車輪に対応するサンプリング処理に用いられる前記時間間隔δvは、前記同期成分波形の抽出に係る複数の車輪のすべてに対して共通のものが用いられ、これらの車輪の基本周波数のうち、最小のものをfmin、最大のものをfmax、前記時間間隔δvの算出に用いたものをfaとしたとき、前記モデル波形に対応する周波数領域波形において、n次の固定同期周波数に対応する山の頂部の周波数範囲は、少なくとも、f0・(fmin/fa) のn倍以上でかつf0・(fmax/fa)のn倍以下の範囲を含むよう構成されてなる請求項6に記載の車輪異常検知装置。 The time interval δ v used for the sampling process corresponding to each wheel is the same for all of the plurality of wheels related to the extraction of the synchronous component waveform, and the minimum of the fundamental frequencies of these wheels is used. F min , the largest one is f max , and the one used to calculate the time interval δ v is f a, and the frequency domain waveform corresponding to the model waveform corresponds to the nth-order fixed synchronization frequency The frequency range of the top of the peak is configured to include at least a range not less than n times f 0 · (f min / f a ) and not more than n times f 0 · (f max / f a ) The wheel abnormality detection device according to claim 6. 前記共通のモデル波形に基づく同期成分波形の抽出に係る全ての車輪に対する基本周波数のうち、最小のものを、前記、共通の時間間隔δvの算出に用いるとともに、Aを予め定められた百分率として、n次の固定同期周波数に対応する山の頂部の周波数範囲を、n次の固定同期周波数の、0〜+A%の範囲として設定してなる請求項7に記載の車輪異常検知装置。 Of fundamental frequency for all of the wheels according to the extraction of the synchronizing component waveform based on said common model waveform, the smallest ones, said, with use in the calculation of the common time interval [delta] v, as a percentage which is determined a A pre The wheel abnormality detection device according to claim 7, wherein a frequency range of the peak corresponding to the n-th fixed synchronization frequency is set as a range of 0 to + A% of the n-th fixed synchronization frequency. n次の固定同期周波数に対応する山の頂部の周波数範囲を、前記共通のモデル波形に基づく同期成分波形の抽出に係る複数の車輪のうち、基本周波数が最小となる車輪と、基本周波数が最大となる車輪との、両方のn次の同期周波数に基づいて設定してなる請求項7に記載の車輪異常検知装置。   The frequency range at the top of the mountain corresponding to the n-th fixed synchronization frequency is set to the wheel having the lowest fundamental frequency and the largest to the fundamental frequency among the plurality of wheels related to the extraction of the synchronization component waveform based on the common model waveform. The wheel abnormality detection device according to claim 7, wherein the wheel abnormality detection device is set on the basis of both n-order synchronization frequencies with the wheel. 前記適応デジタルフィルタ部は、前記共通のモデル波形を構成するデータを、一つ入力するごとに、この共通のモデル波形に基づく同期成分波形の抽出に係る全車輪のそれぞれについて、同期成分波形のデータを一つだけ出力する処理を行うよう構成されてなる請求項7〜9のいずれかに記載の車輪異常検知装置。   Each time the adaptive digital filter unit inputs data constituting the common model waveform, the data of the synchronous component waveform for each of all the wheels related to the extraction of the synchronous component waveform based on the common model waveform. The wheel abnormality detection device according to claim 7, wherein the wheel abnormality detection device is configured to perform a process of outputting only one of the two.
JP2005039215A 2005-02-16 2005-02-16 Wheel abnormality detecting device Withdrawn JP2006224757A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005039215A JP2006224757A (en) 2005-02-16 2005-02-16 Wheel abnormality detecting device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005039215A JP2006224757A (en) 2005-02-16 2005-02-16 Wheel abnormality detecting device

Publications (1)

Publication Number Publication Date
JP2006224757A true JP2006224757A (en) 2006-08-31

Family

ID=36986487

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005039215A Withdrawn JP2006224757A (en) 2005-02-16 2005-02-16 Wheel abnormality detecting device

Country Status (1)

Country Link
JP (1) JP2006224757A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010181319A (en) * 2009-02-06 2010-08-19 Bridgestone Corp Method and device for decision of tire internal fault
JP2015101122A (en) * 2013-11-21 2015-06-04 Ntn株式会社 Abrasion loss detection device of automobile tire
US9463674B1 (en) 2016-01-15 2016-10-11 International Business Machines Corporation Real-time detection of wheel anomalies

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010181319A (en) * 2009-02-06 2010-08-19 Bridgestone Corp Method and device for decision of tire internal fault
JP2015101122A (en) * 2013-11-21 2015-06-04 Ntn株式会社 Abrasion loss detection device of automobile tire
US10112444B2 (en) 2013-11-21 2018-10-30 Ntn Corporation Wear amount detection device for automobile tires
US9463674B1 (en) 2016-01-15 2016-10-11 International Business Machines Corporation Real-time detection of wheel anomalies

Similar Documents

Publication Publication Date Title
CN110307994B (en) Abnormal sound detection device and abnormal sound detection method
EP1813930B1 (en) Tire deformation calculating method and tire deformation calculating apparatus
US20210011457A1 (en) Online monitoring device and system for a 3D printing device
EP2012106B1 (en) Deflection calculating method at tire wheeling time, data storing method at tire wheeling time, and grounding length calculating method at tire wheeling time
EP1906165B1 (en) Method and device for calculating magnitude of cornering force generated in wheel
JP2010165242A (en) Method and system for detecting abnormality of mobile body
EP1813493A1 (en) Tire slipping state detecting method and tire slipping state detecting device
KR20090057417A (en) Method and device for monitoring the noise from a sensor
CN110307895A (en) Allophone decision maker and allophone determination method
CN103842189A (en) Improved method of detecting wear of a tyre
JP3657876B2 (en) Tire pressure estimation device
CN109655267A (en) A kind of automobile hub bearing fault signature extracting method
JP2006224757A (en) Wheel abnormality detecting device
JP4307591B2 (en) Method for monitoring vehicles with acceleration sensors, in particular helicopter planetary gearboxes
CN114734762A (en) Tire pressure monitoring method, controller, system and electronic equipment
CN109649399B (en) Driver distraction judgment method based on transverse driving speed
CN109086508B (en) Threshold value dynamic matching method for automobile emergency call system
CN105940608A (en) Evaluation method for sensor signals
TW201633025A (en) Diagnostic method for malfunction mode of machine tool main shaft and system thereof
JP4279431B2 (en) Operating state determination device
CN115610403A (en) Torque compensation method and device, computer readable storage medium and vehicle
CN208101702U (en) A kind of non-security driving detection device of automobile based on image recognition
CN111400648B (en) Method for quantitatively evaluating waveform intensity of frontal collision of automobile
JP5187850B2 (en) Road surface state detection method and road surface state detection device
FR3116900B1 (en) Method for determining the occupancy state of a seat and corresponding determination system

Legal Events

Date Code Title Description
RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20060609

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080109

A761 Written withdrawal of application

Free format text: JAPANESE INTERMEDIATE CODE: A761

Effective date: 20091028