JP2018096753A - Data processing method, estimation method of tire, evaluation method of vibration ride comfort performance, and data processing device - Google Patents
Data processing method, estimation method of tire, evaluation method of vibration ride comfort performance, and data processing device Download PDFInfo
- Publication number
- JP2018096753A JP2018096753A JP2016239574A JP2016239574A JP2018096753A JP 2018096753 A JP2018096753 A JP 2018096753A JP 2016239574 A JP2016239574 A JP 2016239574A JP 2016239574 A JP2016239574 A JP 2016239574A JP 2018096753 A JP2018096753 A JP 2018096753A
- Authority
- JP
- Japan
- Prior art keywords
- tire
- value
- input data
- convergence
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims description 83
- 238000003672 processing method Methods 0.000 title claims description 33
- 238000012545 processing Methods 0.000 title claims description 26
- 238000000034 method Methods 0.000 title claims description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 53
- 238000005259 measurement Methods 0.000 claims abstract description 40
- 230000005540 biological transmission Effects 0.000 claims abstract description 14
- 238000012546 transfer Methods 0.000 claims description 51
- 238000013016 damping Methods 0.000 claims description 33
- 238000005096 rolling process Methods 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 11
- 230000008602 contraction Effects 0.000 claims description 8
- 239000011324 bead Substances 0.000 claims description 6
- 239000000725 suspension Substances 0.000 claims description 6
- 238000010521 absorption reaction Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 16
- 238000012360 testing method Methods 0.000 description 11
- 230000001953 sensory effect Effects 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 6
- 230000008859 change Effects 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Landscapes
- Tires In General (AREA)
Abstract
Description
本発明は、コンピュータが行うデータ処理方法、タイヤの評価方法、車両の振動乗り心地性能の評価方法、及びデータ処理装置に関する。 The present invention relates to a data processing method performed by a computer, a tire evaluation method, a vehicle vibration riding comfort performance evaluation method, and a data processing apparatus.
タイヤの振動乗り心地の評価は、タイヤを、突起等がドラム路面上に設けられた室内ドラム上を走行させ、この時タイヤの軸力に発生する振動を計測することによって行われる。また、タイヤの振動乗り心地の評価は、評価しようとするタイヤを車両に装着し、車両を凹凸のある路面を走行させて、走行中のタイヤ回転軸に作用する加速度や車両の室内騒音を計測することによって行われる。 The evaluation of tire vibration riding comfort is performed by running the tire on an indoor drum having protrusions or the like provided on the drum road surface, and measuring the vibration generated in the axial force of the tire at this time. In addition, the evaluation of tire vibration ride comfort is done by mounting the tire to be evaluated on the vehicle, running the vehicle on a rough road surface, and measuring the acceleration acting on the rotating tire rotation axis and the vehicle interior noise. Is done by doing.
例えば、タイヤを室内ドラム上に走行させたときの、タイヤの振動特性を簡便かつ精度よく評価する技術が知られている（特許文献１）。
具体的には、タイヤの評価方法は、転動中のタイヤに入力を与えた後にタイヤを自由減衰振動させ、入力を与えた後におけるタイヤの振動に関する情報を取得し、この振動に関する情報から得られた時刻歴の振動の波形に合うように、２自由度以上の自由減衰振動の式を定め、得られた自由減衰振動の式を用いて、転動中における前記タイヤを評価するための情報を求める。これにより、転動中におけるタイヤの固有振動数及び減衰比（減衰係数）等をタイヤ評価情報として取得する。
For example, a technique for simply and accurately evaluating the vibration characteristics of a tire when the tire is run on an indoor drum is known (Patent Document 1).
Specifically, in the tire evaluation method, after giving input to the rolling tire, the tire is subjected to free-damping vibration, and information on the vibration of the tire after giving the input is obtained and obtained from the information on this vibration. Information for evaluating the tire during rolling using a formula of a free damping vibration having two or more degrees of freedom so as to fit a vibration waveform of a given time history. Ask for. Thus, the natural frequency and damping ratio (damping coefficient) of the tire during rolling are acquired as tire evaluation information.
上記タイヤの評価方法では、タイヤの固有振動数及び減衰比がタイヤ評価情報として得られる。しかし、同じ凹凸を走行してもタイヤが路面から受ける入力は異なる場合があるため、固有振動数及び減衰比が同じであってもタイヤの評価が異なる場合がある。このため、路面の凹凸からタイヤが受ける入力を適切に取得して評価できることが好ましい。 In the tire evaluation method, the natural frequency and the damping ratio of the tire are obtained as tire evaluation information. However, since the tire may receive different inputs from the road even when traveling on the same unevenness, the evaluation of the tire may be different even if the natural frequency and the damping ratio are the same. For this reason, it is preferable that the input received by the tire from the unevenness of the road surface can be appropriately acquired and evaluated.
そこで、本発明は、路面上の段差を通過するときにタイヤ軸力に生じるタイヤ軸力の振動からタイヤ入力及びタイヤの振動伝達特性における特性パラメータを算出することができるデータ処理方法及びデータ処理装置、さらには、算出したタイヤ入力及び特性パラメータの少なくとも一方を用いたタイヤの評価方法及び車両の振動乗り心地性能の評価方法を提供することを目的とする。 Therefore, the present invention provides a data processing method and a data processing apparatus capable of calculating characteristic parameters in tire input and tire vibration transmission characteristics from vibration of tire axial force generated in tire axial force when passing through a step on a road surface. A further object of the present invention is to provide a tire evaluation method and a vehicle vibration riding comfort performance evaluation method using at least one of the calculated tire input and characteristic parameters.
本発明の一態様は、路面を転動するタイヤが前記路面上の段差を通過するときにタイヤ回転軸に生じるタイヤ軸力Ｏの振動から、前記タイヤが前記段差から受けるタイヤ入力Ｉと、前記タイヤの軸力Ｏと前記タイヤ入力との間の振動伝達特性Ｈを、コンピュータが算出するデータ処理方法である。
当該データ処理方法は、
前記コンピュータが、前記タイヤ軸力Ｏの計測データから、設定された前記振動伝達特性Ｈの特性パラメータの仮値を用いて、前記タイヤ入力Ｉの仮入力データを算出する入力データ算出ステップと、
前記コンピュータが、前記入力データ算出ステップで算出した前記仮入力データのうち、前記タイヤが前記段差を通過する前後のデータの値を０にした修正仮入力データと前記タイヤ軸力Ｏの計測データから、前記振動伝達特性Ｈの特性パラメータの仮値を修正した修正仮値を算出するパラメータ値算出ステップと、を備え、
前記仮入力データのうち前記タイヤが前記段差を通過する前後のデータの値が０に収束する条件を少なくとも含む収束条件を満足するまで、前記コンピュータは、前記修正仮値を、前記特性パラメータの前記仮値として設定して前記入力データ算出ステップ及び前記パラメータ値算出ステップを繰り返し行うことにより、前記収束条件を満足するときの前記タイヤ入力Ｉの仮入力データと前記特性パラメータの仮値のそれぞれを前記タイヤ入力Ｉの収束入力データ及び前記特性パラメータの収束値として算出する。
One aspect of the present invention is the tire input I received from the step by the tire axial force O generated on the tire rotation axis when the tire rolling on the road passes through the step on the road surface, This is a data processing method in which a computer calculates a vibration transfer characteristic H between the tire axial force O and the tire input.
The data processing method is
An input data calculating step in which the computer calculates temporary input data of the tire input I from the measurement data of the tire axial force O using the set temporary value of the characteristic parameter of the vibration transfer characteristic H;
From the temporary input data calculated by the computer in the input data calculation step, corrected temporary input data in which the value of data before and after the tire passes the step is set to 0, and measurement data of the tire axial force O A parameter value calculating step of calculating a corrected temporary value obtained by correcting the temporary value of the characteristic parameter of the vibration transfer characteristic H,
Until the computer satisfies the convergence condition including at least the condition that the value of the data before and after the tire passes through the step in the temporary input data satisfies the convergence condition, the computer sets the corrected temporary value as the characteristic parameter. By setting the provisional value and repeating the input data calculation step and the parameter value calculation step, each of the provisional input data of the tire input I and the provisional value of the characteristic parameter when the convergence condition is satisfied is obtained. The convergence input data of the tire input I and the convergence value of the characteristic parameter are calculated.
その際、前記収束条件は、さらに、前記修正仮入力データと前記振動伝達特性の特性パラメータの修正仮値から算出した前記タイヤ軸力Ｏの予測データが前記計測データに収束する条件を含む、ことが好ましい。 At this time, the convergence condition further includes a condition that the predicted data of the tire axial force O calculated from the corrected temporary input data and the corrected temporary value of the characteristic parameter of the vibration transfer characteristic converges to the measurement data. Is preferred.
前記タイヤ軸力Ｏは、タイヤ上下軸力及びタイヤ前後軸力のいずれか一方を少なくとも含むことが好ましい。 The tire axial force O preferably includes at least one of a tire vertical axial force and a tire longitudinal axial force.
前記振動伝達特性は、1自由度振動伝達特性であり、前記特性パラメータは、非減衰固有角振動数と減衰係数を含む、ことが好ましい。 Preferably, the vibration transfer characteristic is a one-degree-of-freedom vibration transfer characteristic, and the characteristic parameters include a non-damped natural angular frequency and a damping coefficient.
本発明の他の一形態であるタイヤの評価方法は、
前記データ処理方法において前記収束入力データ及び前記特性パラメータの収束値を算出するステップと、
前記収束入力データ及び前記特性パラメータの前記収束値の少なくとも一方を用いて、前記コンピュータが前記タイヤの振動乗り心地性能を評価するステップと、を備える。
A tire evaluation method according to another embodiment of the present invention is as follows.
Calculating a convergence value of the convergence input data and the characteristic parameter in the data processing method;
The computer using the convergence input data and the convergence value of the characteristic parameter to evaluate the vibration riding comfort performance of the tire.
前記タイヤ評価方法では、前記データ処理方法において前記収束入力データ及び前記特性パラメータの収束値を算出するステップと、
前記特性パラメータの前記収束値のうち、前記非減衰固有角振動数及び前記減衰係数の収束値からタイヤのサイド・ビード部の振動の吸収の程度を表すダンパ要素を前記コンピュータが評価するステップと、を備えることが好ましい。
In the tire evaluation method, calculating the convergence input data and the convergence value of the characteristic parameter in the data processing method;
The computer evaluating a damper element representing a degree of vibration absorption of a tire side bead from the convergence value of the non-damped natural angular frequency and the damping coefficient among the convergence value of the characteristic parameter; It is preferable to provide.
前記タイヤの振動乗り心地性能を評価するステップでは、
前記収束入力データの値の絶対値の積分値、あるいは、前記タイヤが前記段差を通過する期間における前記収束入力データの時間微分の絶対値の平均値を用いて、前記タイヤの振動乗り心地性能のうちの性能１の評価を行う、ことが好ましい。
In the step of evaluating the vibration riding comfort performance of the tire,
Using the integral value of the absolute value of the convergent input data or the average value of the absolute value of the time derivative of the convergent input data during the period when the tire passes the step, the vibration ride comfort performance of the tire It is preferable to evaluate performance 1 of them.
前記タイヤの振動乗り心地性能を評価するステップでは、
前記タイヤの転動する計測用転動速度と異なる評価用転動速度における前記タイヤの振動乗り心地性能のうちの性能２の評価を行うために、前記計測用転動速度における前記収束入力データの時間軸上の波形を前記評価用転動速度に応じて伸縮することで得られる伸縮波形を表した伸縮収束入力データと、前記特性パラメータの前記収束値を有する振動伝達特性を用いてコンボリューション処理を行い、該コンボリューション処理を行うことで算出された前記評価用転動速度に対応した前記タイヤ軸力Ｏの予測データの絶対値の積分値と前記伸縮収束入力データの絶対値の積分値との差分により、前記性能２の評価を行う、ことが好ましい。
In the step of evaluating the vibration riding comfort performance of the tire,
In order to evaluate
本発明のさらに他の一態様は、車両の振動乗り心地性能の評価方法は、
前記データ処理方法において前記収束入力データ及び前記特性パラメータの収束値を算出するステップと、
前記特性パラメータで表した前記タイヤの振動伝達特性を、車両のサスペンションをモデル化したサスペンションモデルに組み合わせた統合モデルに、前記収束入力データを入力することにより、車両の振動乗り心地性能の評価を行うステップと、を備える、ことが好ましい。
According to still another aspect of the present invention, there is provided a method for evaluating vibration ride comfort performance of a vehicle.
Calculating a convergence value of the convergence input data and the characteristic parameter in the data processing method;
Evaluation of the vehicle's vibration ride comfort performance is performed by inputting the convergence input data into an integrated model in which the vibration transmission characteristics of the tire represented by the characteristic parameters are combined with a suspension model obtained by modeling a vehicle suspension. Preferably comprising steps.
本発明のさらに他の一態様は、路面を転動するタイヤが前記路面上の段差を通過するときにタイヤ回転軸に生じるタイヤ軸力Ｏの振動から、前記タイヤが前記段差から受けるタイヤ入力Ｉと、前記タイヤの軸力Ｏと前記タイヤ入力との間の振動伝達特性Ｈを算出するデータ処理装置であって、
前記タイヤ軸力Ｏの計測データを取得するデータ取得部と、
前記計測データから、設定された前記振動伝達特性Ｈの特性パラメータの仮値を用いて、前記タイヤ入力Ｉの仮入力データを算出する入力データ算出部と、
前記入力データ算出部で算出した前記仮入力データのうち、前記タイヤが前記段差を通過する前後のデータの値を０にした修正仮入力データと前記タイヤ軸力Ｏの計測データから、前記振動伝達特性の特性パラメータの仮値を修正した修正仮値を算出するパラメータ値算出部と、
前記仮入力データのうち前記タイヤが前記段差を通過する前後のデータの値が０に収束する条件を少なくとも含む収束条件を満足するまで、前記修正仮値を、前記特性パラメータの前記仮値として設定して前記仮入力データの算出及び前記修正仮値の算出を繰り返し行うよう前記入力データ算出部及び前記パラメータ値算出部を制御する制御部と、
前記収束条件を満足するときの前記タイヤ入力Ｉの仮入力データと前記特性パラメータの仮値のそれぞれを前記タイヤ入力Ｉの収束入力データ及び前記特性パラメータの収束値として算出する出力部と、を備える。
Still another aspect of the present invention is a tire input I which the tire receives from the step from vibration of a tire axial force O generated on a tire rotation shaft when a tire rolling on the road passes through the step on the road surface. And a data processing device for calculating a vibration transfer characteristic H between the tire axial force O and the tire input,
A data acquisition unit for acquiring measurement data of the tire axial force O;
An input data calculation unit that calculates temporary input data of the tire input I using a temporary value of the set characteristic parameter of the vibration transfer characteristic H from the measurement data;
Of the temporary input data calculated by the input data calculation unit, the vibration transmission is based on the corrected temporary input data in which the value of the data before and after the tire passes through the step is 0 and the measurement data of the tire axial force O. A parameter value calculation unit for calculating a corrected temporary value obtained by correcting the temporary value of the characteristic parameter of the characteristic;
The corrected provisional value is set as the provisional value of the characteristic parameter until a convergence condition including at least a condition in which the value of the data before and after the tire passes through the step of the provisional input data converges to 0 is satisfied. A control unit that controls the input data calculation unit and the parameter value calculation unit so as to repeatedly calculate the temporary input data and the correction temporary value;
An output unit that calculates the temporary input data of the tire input I and the temporary value of the characteristic parameter when the convergence condition is satisfied as the convergence input data of the tire input I and the convergence value of the characteristic parameter; .
前記振動伝達特性は、１自由度振動伝達特性であり、前記特性パラメータは、非減衰固有角振動数と減衰係数を含み、
さらに、前記特性パラメータの前記収束値のうち、前記減衰係数の収束値からタイヤのサイド・ビード部の振動の吸収の程度を表すダンパ要素を前記コンピュータが評価する評価部を備える、ことが好ましい。
The vibration transfer characteristic is a one-degree-of-freedom vibration transfer characteristic, and the characteristic parameters include a non-damped natural angular frequency and a damping coefficient,
Further, it is preferable that the computer further includes an evaluation unit that evaluates a damper element that represents the degree of vibration absorption of the side bead portion of the tire from the convergence value of the damping coefficient among the convergence values of the characteristic parameters.
上述のデータ処理方法及びデータ処理装置によれば、路面上の凹凸を通過するときに生じるタイヤ軸力の振動からタイヤ入力及びタイヤの振動伝達特性における特性パラメータを算出することができる。さらに、このタイヤ入力及び特性パラメータを用いてタイヤを評価することができる。 According to the above-described data processing method and data processing apparatus, it is possible to calculate the characteristic parameters in the tire input and the tire vibration transfer characteristics from the vibration of the tire axial force generated when passing through the unevenness on the road surface. Further, the tire can be evaluated using the tire input and characteristic parameters.
以下、本発明のデータ処理方法、タイヤの評価方法、車両の振動乗り心地性能の評価方法、及びデータ処理装置について添付の図面を参照しながら詳細に説明する。 Hereinafter, a data processing method, a tire evaluation method, a vehicle vibration riding comfort performance evaluation method, and a data processing device according to the present invention will be described in detail with reference to the accompanying drawings.
図１は、本実施形態のデータ処理方法及びタイヤの評価方法を実行するデータ処理装置１０の構成を示す図である。データ処理装置１０は、ＣＰＵ１２及びメモリ１４を備えるコンピュータで構成された装置である。データ処理装置１０には、ディスプレイ１６及びキーボードやマウス等を含む入力操作系１８と接続されている。
FIG. 1 is a diagram illustrating a configuration of a
さらに、データ処理装置１０は、タイヤドラム試験機２０と接続されている。タイヤドラム試験機２０の試験ドラム２２のドラム面には、突起２２ａが設けられている。試験ドラム２２を所定の速度で転回転させ、この試験ドラム２２のドラム面にタイヤＴを所定の荷重で押し付けてタイヤＴを転動させる。このとき、タイヤＴは突起２２ａを通過する度に突起２２ａから振動の入力を受け、この入力は、タイヤＴの振動伝達特性を介してタイヤ回転軸２４に伝達される。タイヤ回転軸２４には、タイヤ軸力を計測するセンサが設けられている。このセンサからの出力がタイヤ軸力Ｏにおける振動の計測データであり、このタイヤ軸力の計測データがデータ処理装置１０に送られる。なお、タイヤ軸力Ｏは、ドラム面から受ける一定の負荷荷重の反力と、突起２２ａをタイヤＴが通過するときに生じる振動の力を含むが、センサは、一定の負荷荷重の反力（ＤＣ成分）は出力せず、振動のみを出力する。
Further, the
データ処理装置１０は、試験ドラム２２のドラム面を転動するタイヤＴがこのドラム面上の突起２２ａを通過するときにタイヤ回転軸に生じるタイヤ軸力Ｏから、タイヤＴが突起２２ａを通過するときに受けるタイヤ入力Ｉと、タイヤ軸力Ｏとタイヤ入力Ｉとの間の振動伝達特性Ｈを算出するデータ処理を行う。本実施形態では、タイヤＴを試験ドラム２２のドラム面上を転動させるが、平坦な路面上を転動させてもよい。また、本実施形態では、突起２２ａを用いるが、突起２２ａには限定されず、路面に凹部あるいは凸部の段差があればよい。
In the
データ処理装置１０は、メモリ１４に記憶されたプログラムを呼び出して起動することにより、ソフトウェアモジュールが形成される。すなわち、データ処理装置１０では、データ取得部３０、入力データ算出部３２、パラメータ値算出部３４、制御部３６、出力部３８、及び評価部４０がソフトウェアモジュールとして形成される。これらの部分の実質的な演算等は、ＣＰＵ１２が行う。
The
データ取得部３０は、タイヤドラム試験機２０から送られるタイヤ軸力Ｏの計測データを取得する。計測データは、メモリ１４に記憶される。
入力データ算出部３２は、タイヤ軸力Ｏの計測データから、設定された振動伝達特性Ｈの特性パラメータの仮値を用いて、タイヤ入力Ｉの仮入力データを算出する。ここで、振動伝達特性Ｈは、例えば１自由度の振動伝達特性の場合、すなわち、２次伝達特性の場合、特性パラメータは非減衰固有角振動数と減衰係数を含む。タイヤＴでは、特性パラメータの値は凡そ一定の範囲内にあるので、入力データ算出部３２は、この一定の範囲内のある値を仮値として予め設定しておく。この特性パラメータの仮値は、後述するように繰り返し計算により修正され、最終的に収束した値になる。設定された振動伝達特性Ｈの特性パラメータの仮値を用いて、タイヤ入力Ｉの仮入力データを算出する方法については後述する。
The
The input
パラメータ値算出部３４は、入力データ算出部３２で算出した仮入力データのうち、タイヤＴが突起２２ａを通過する前後のデータの値を０にした修正仮入力データとタイヤ軸力Ｏの計測データから、振動伝達特性Ｈの特性パラメータの仮値を修正した修正仮値を算出する。
すなわち、タイヤＴが突起２２ａを通過するときは、突起２２ａからタイヤ入力Ｉを受けるが、突起２２ａ以外の部分をタイヤＴが通過するとき、試験ドラム２２のドラム面から受ける入力データＩは０になるので（振動を生じさせるタイヤ入力Ｉはないので）、仮入力データは、本来、タイヤＴが突起２２ａを通過する前後において０にならなければならない。このため、本実施形態のパラメータ値算出部３４は、入力データ算出部３２で算出した仮入力データのうち、タイヤＴが突起２２ａを通過する前後のデータの値を０にした修正仮入力データを生成する。この修正仮入力データとタイヤ軸力Ｏの計測データとから、振動伝達特性Ｈの特性パラメータの仮値を修正した修正仮値を算出する。修正仮入力データとタイヤ軸力Ｏの計測データが既知であるので、特性パラメータの修正仮値を算出することができる。振動伝達特性Ｈの特性パラメータの仮値を修正した修正仮値を算出する方法については後述する。
The parameter
That is, when the tire T passes through the
制御部３６は、ソフトウェアモジュールの各部分の動作の管理及び制御を行う部分である。制御部３６は、入力データ算出部３２で算出される仮入力データのうちタイヤＴが突起２２ａを通過する前後のデータの値が０に収束する条件を少なくとも含む収束条件を満足するまで、パラメータ値算出部３４で算出した修正仮値を、特性パラメータの仮値として設定して入力データ算出部３２における仮入力データの算出及びパラメータ値算出部３４における修正仮値の算出を繰り返し行うよう入力データ算出部３２及びパラメータ値算出部３４を制御する。
The
出力部３８は、収束条件を満足するときのタイヤ入力Ｉの仮入力データと特性パラメータの仮値のそれぞれをタイヤ入力Ｉの収束入力データ及び特性パラメータの収束値として算出する。算出した収束入力データ及び特性パラメータの収束値は、メモリ１４に記憶され、また、ディスプレイ１６や図示されないプリンタに出力される。
The
評価部４０は、収束入力データ及び特性パラメータの収束値の少なくとも一方を用いて、タイヤＴの振動乗り心地性能を評価する。例えば、振動伝達特性Ｈが、１自由度の振動伝達特性、すなわち、２次伝達特性であって、特性パラメータは非減衰固有角振動数と減衰係数を含む場合、評価部４０は、特性パラメータから算出されるダンパＣ（Ｎ・秒／ｍｍ）を用いてタイヤＴの振動乗り心地性能を評価する。あるいは、評価部４０は、収束入力データの波形の積分値を用いて、振動乗り心地性能の各性能を評価する。この点については、後ほど詳述する。評価結果は、メモリ１２に記憶され、また、ディスプレイ１６や図示されないプリンタに出力される。
The
このようなデータ処理装置１０において、データ処理方法及びタイヤの評価方法をより詳細に説明する。
図２は、本実施形態のデータ処理方法の一例を説明する図である。図３は、本実施形態のデータ処理方法及びタイヤの評価方法の処理の流れを示す図である。図２に示す例では、振動伝達特性Ｈが、１自由度の振動伝達特性、すなわち、２次伝達特性である場合を示している。データ取得部３０は、タイヤ軸力Ｏの計測データを取得する（ステップＳ１０）。例えば、図２に示すような波形を取得する。この後、２次伝達特性のインパルス応答は、非減衰固有角振動数ωｎと減衰係数ζ（ζ＜１）の設定した値を用いて図２中に示す式で表されるので、この式を利用して、タイヤ入力Ｉの仮入力データを算出する（ステップＳ２０）。このとき、非減衰固有角振動数ωｎと減衰係数ζの値は、凡そ一定の範囲内にあるので、入力データ算出部３２は、この一定の範囲内のある値を初期の仮値として予め設定しておく。初期の仮値として、例えば、ωｎ＝２π・７０（ｒａｄ／秒）、ζ＝０.０５と定める。このとき、タイヤ軸力Ｏと、非減衰固有角振動数ωｎと減衰係数ζの値を含む図２中の式を用いて、タイヤ軸力Ｏからタイヤ入力Ｉを算出する処理方法は、デコンボリューション処理であり、公知の方法でタイヤ入力Ｉを算出することができる。具体的には、タイヤ軸力Ｏの計測データの時系列波形の離散値と、図２中の式を具体的に離散化して表現した行列（具体的には、タイヤ入力Ｉからタイヤ軸力Ｏを求めるための、Ｈ（ｔ；ωｎ，ζ）を用いた行列の逆行列）を用いて行列演算を行って、タイヤ入力Ｉの仮入力データを算出することができる。
なお、タイヤ入力Ｉからタイヤ軸力Ｏを算出する処理方法は、コンボリューション処理であり、公知の方法でタイヤ軸力Ｏを算出することができる。具体的には、タイヤ入力Ｉの時系列波形の離散値と、図２中の式を具体的に離散化して表現した行列を用いて行列演算を行って、タイヤ軸力Ｏを算出することができる。
In such a
FIG. 2 is a diagram for explaining an example of the data processing method according to the present embodiment. FIG. 3 is a diagram illustrating a processing flow of the data processing method and the tire evaluation method of the present embodiment. In the example shown in FIG. 2, the vibration transfer characteristic H is a one-degree-of-freedom vibration transfer characteristic, that is, a secondary transfer characteristic. The
The processing method for calculating the tire axial force O from the tire input I is a convolution process, and the tire axial force O can be calculated by a known method. Specifically, the tire axial force O may be calculated by performing a matrix operation using a discrete value of the time series waveform of the tire input I and a matrix obtained by specifically discretizing the expression in FIG. it can.
次に、パラメータ値算出部３４は、算出した仮入力データのうち、タイヤＴが突起２２ａを通過する前後のデータの値を０にした修正仮入力データとタイヤ軸力Ｏの計測データから、振動伝達特性Ｈの特性パラメータの仮値を修正した修正仮値を算出する（ステップＳ３０）。修正仮値の算出は、例えば、振動伝達特性Ｈの特性パラメータの仮値を種々変更した振動伝達特性Ｈと、仮入力データとを用いたコンボリューション処理の処理結果が、タイヤ軸力Ｏの計測データに最も近似するような特性パラメータの仮値を探索する。特性パラメータの仮値の探索は、ニュートン・ラフソン法等を用いて行うことができる。
次に、制御部３６は、ステップＳ２０で算出される仮入力データのうちタイヤＴが突起２２ａを通過する前後のデータの値が０に収束する条件を少なくとも含む収束条件を満足するか否かを判定する（ステップＳ４０）。判定の結果、収束条件を満足する場合、制御部３６は、ステップＳ２０で算出した仮入力データを収束入力データとし、ステップＳ３０で算出した修正仮値を特性パラメータの収束値とする。この後、タイヤ評価に移行する（ステップＳ６０）。
一方、判定の結果、収束条件を満足しない場合、ステップＳ３０で算出した修正仮値を、特性パラメータの仮値として設定して（ステップＳ５０）、ステップＳ２０に戻る。こうして、ステップＳ２０及びステップＳ３０を繰り返し行う。すなわち、ステップＳ４０において収束条件を満足するまで、ステップＳ２０及びステップＳ３０を繰り返し行う。
なお、タイヤの評価については、後述する。
Next, the parameter
Next, the
On the other hand, if the convergence condition is not satisfied as a result of the determination, the corrected provisional value calculated in step S30 is set as the provisional value of the characteristic parameter (step S50), and the process returns to step S20. Thus, step S20 and step S30 are repeated. That is, step S20 and step S30 are repeated until the convergence condition is satisfied in step S40.
The tire evaluation will be described later.
このように、本実施形態では、タイヤ軸力Ｏから、タイヤＴの振動伝達特性Ｈの特性パラメータの仮値を用いて算出した仮入力データのうち、タイヤＴが突起２２ａを通過する前後の期間の値が０になることを利用して、仮入力データにおいて、タイヤＴが突起２２ａを通過する前後の期間の値が０に収束するまで、タイヤＴの振動伝達特性Ｈの特性パラメータの仮値を変更するので、タイヤＴの振動伝達特性Ｈの特性パラメータの値及びタイヤ入力のデータを精度よく算出することができる。
Thus, in the present embodiment, the period before and after the tire T passes the
図４（ａ）〜（ｆ）は、本実施形態のデータ処理方法による演算結果の一例を示す図である。
図４（ａ）は、入力データ算出部３２が、タイヤ軸力Ｏから、非減衰固有角振動数ωｎと減衰係数ζの初期の値を仮値として用いて仮入力データを算出した結果の一例を示す。この場合、非減衰固有角振動数ωｎと減衰係数ζの初期値は、実際の非減衰固有角振動数ωｎと減衰係数ζの値とずれているので、タイヤＴが突起２２ａを通過した後において、仮入力データは０になっていない。このため、パラメータ算出部３４は、図４（ｂ）に示すように、タイヤＴが突起２２ａを通過した前後において、仮入力データの値を０にした修正仮入力データをつくる。パラメータ算出部３４は、次に、図４（ｂ）に示す修正仮入力データと、タイヤ軸力Ｏの計測データとを用いて、非減衰固有角振動数ωｎと減衰係数ζの値を初期の仮値から修正した修正仮値を算出する。この場合の修正仮値の算出では、非減衰固有角振動数ωｎと減衰係数ζの仮値を種々変更しながら、変更の度に、タイヤ軸Ｏの時系列波形を算出する。算出したタイヤ軸Ｏの時系列波形が、タイヤ軸力Ｏの計測データに最もよく近似するまで仮値を探索する。算出するタイヤ軸Ｏの時系列波形をタイヤ軸力Ｏの計測データに最もよく近似させる非減衰固有角振動数ωｎと減衰係数ζの仮値を、修正仮値とする。図４（ｃ）は、算出するタイヤ軸Ｏの時系列波形がタイヤ軸力Ｏの計測データに最もよく近似する状態の一例を示している。
FIGS. 4A to 4F are diagrams illustrating examples of calculation results obtained by the data processing method according to the present embodiment.
FIG. 4A shows the result of the temporary input data calculated by the input
図４（ｄ）は、入力データ算出部３２が、パラメータ算出部３４が算出した非減衰固有角振動数ωｎと減衰係数ζの修正仮値を用いてタイヤ軸力Ｏから、仮入力データを算出した結果の一例を示す。この場合、非減衰固有角振動数ωｎと減衰係数ζの修正仮値は、実際の非減衰固有角振動数ωｎと減衰係数ζの値とずれているので、タイヤＴが突起２２ａを通過した後において、仮入力データは０になっていない。しかし、タイヤＴが突起２２ａを通過した後の仮入力データの値は、図４（ａ）におけるタイヤＴが突起２２ａを通過した後の仮入力データの値は小さくなっている。このことは、非減衰固有角振動数ωｎと減衰係数ζの値が、実際の非減衰固有角振動数ωｎと減衰係数ζの値に近づいていることを示している。
さらに、パラメータ算出部３４は、図４（ｅ）に示すように、タイヤＴが突起２２ａを通過した前後において、仮入力データの値を０にする。パラメータ算出部３４は、次に、図４（ｅ）に示す仮入力データと、タイヤ軸力Ｏの計測データとを用いて、非減衰固有角振動数ωｎと減衰係数ζの修正仮値をさらに修正した修正仮値を、上述した方法と同じ方法で算出する。図４（ｆ）は、算出するタイヤ軸Ｏの時系列波形がタイヤ軸力Ｏの計測データに最もよく近似する状態の一例を示している。
このようにして、図４（ｄ）に示す仮入力データにおけるタイヤＴの突起２２ａの通過前後の値が許容範囲になるまで、理想的には０になるまで、繰り返し、修正仮入力データと修正仮値を求める。
In FIG. 4D, the input
Further, as shown in FIG. 4E, the
In this way, until the value before and after the passage of the
このとき、タイヤＴが突起２２ａを通過する前後の期間は、仮入力データの波形から自動的に求めることが好ましい。図５は、タイヤＴが突起２２ａを通過する前後の期間を、タイヤＴが突起２２ａを通過する期間を定めることにより定める一例を示す図である。
図５に示すように、タイヤ入力Ｉの仮入力データが０（Ｎ）の横軸とタイヤ入力Ｉの波形で囲まれた部分の面積を順番にＳ１，Ｓ２，Ｓ３，・・・・として求め、この面積の合計の５％〜１０％の面積を有する部分であって、多くても最大面積と２番目に大きな面積の部分を取り出し、この最大面積の部分と２番目の面積の部分のうち、最も早い時間から最も遅い時間の間をタイヤＴが突起２２ａを通過する期間として定める。したがって、この通過期間の前後の期間が、タイヤＴが突起２２ａを通過する前後の期間となる。図５に示す例では、最大面積Ｓ１の部分の期間はｔ１〜ｔ２であり、２番目の大きさの面積Ｓ２の部分が占める期間はｔ３〜ｔ４である。したがって、この場合のタイヤＴが突起２２ａを通過する期間は、ｔ１〜ｔ４となり、タイヤＴが突起２２ａを通過する前後の期間は、ｔ１以前、及びｔ４以降である。なお、面積Ｓ１，Ｓ２，Ｓ３，・・・・の合計の５％〜１０％の面積を有する部分が１つしかない場合、その部分の最初の時間と最後の時間の間を、タイヤＴが突起２２ａを通過するとする。
At this time, it is preferable that the period before and after the tire T passes the
As shown in FIG. 5, the area of the portion surrounded by the tire input I temporary input data of 0 (N) and the tire input I waveform is sequentially determined as S1, S2, S3,. , A portion having an area of 5% to 10% of the total of the areas, and taking out the maximum area and the second largest area portion at most, out of the maximum area portion and the second area portion The period from the earliest time to the latest time is defined as a period during which the tire T passes through the
上記実施形態では、タイヤＴのトラム面に対して直交する方向の力、いわゆる上下方向の力を対象として説明したが、ドラム面の回転方向の成分、いわゆる前後方向の力を対象とすることもできる。図６（ａ）は、突起２２ａを通過するときに生じるタイヤＴの前後方向のタイヤ軸力Ｏの振動の一例を示す図であり、図６（ｂ）は、突起２２ａを通過するときに生じるタイヤＴの前後方向のタイヤ軸力Ｏの振動から、上述のデータ処理方法を用いて算出したタイヤＴの前後方向のタイヤ入力Ｉ（収束入力データ）の振動の一例を示す図である。
このように、本実施形態のタイヤ軸力Ｏは、タイヤ上下軸力及びタイヤ前後軸力のいずれか一方を少なくとも含む、データ処理方法であることが好ましい。
In the above-described embodiment, the force in the direction orthogonal to the tram surface of the tire T, that is, the so-called vertical force has been described, but the component in the rotation direction of the drum surface, that is, the so-called front-rear direction force may be used. it can. FIG. 6A is a diagram illustrating an example of the vibration of the tire axial force O in the front-rear direction of the tire T generated when passing through the
As described above, the tire axial force O of the present embodiment is preferably a data processing method including at least one of the tire vertical axial force and the tire longitudinal axial force.
本実施形態における仮入力データのうち、タイヤＴが突起２２ａを通過する前後の期間の値が０に収束する収束条件は、タイヤＴが突起２２ａを通過する前後の期間における仮入力データの０からの残差の二乗和Ｑ１が所定の閾値以下であることを含むことが好ましい。
Among the temporary input data in the present embodiment, the convergence condition in which the value of the period before and after the tire T passes the
さらに、この収束条件は、さらに、修正仮入力データと振動伝達特性Ｈの特性パラメータの修正仮値から算出したタイヤ軸力Ｏの予測データが、取得したタイヤＴのタイヤ軸力Ｏの計測データに収束する条件を含むことが好ましい。具体的には、収束条件は、上記予測データと計測データの残差の二乗和Ｑ２と上記二乗和Ｑ１の和が予め定めた閾値以下であることを含むことが好ましい。二乗和Ｑ１と二乗和Ｑ２の和によって収束条件を判定することで、タイヤ入力Ｉのデータを精度良く求めることができる。 Further, the convergence condition is that the predicted data of the tire axial force O calculated from the corrected temporary input data and the corrected temporary value of the characteristic parameter of the vibration transfer characteristic H is converted into the acquired measurement data of the tire axial force O of the tire T. It is preferable to include a condition for convergence. Specifically, it is preferable that the convergence condition includes that the sum of the square sum Q2 of the residual of the prediction data and the measurement data and the sum of the square sum Q1 is equal to or less than a predetermined threshold value. By determining the convergence condition based on the sum of the square sum Q1 and the square sum Q2, the tire input I data can be obtained with high accuracy.
本実施形態における振動伝達特性Ｈは、1自由度振動伝達特性あるいは２自由度振動伝達特性、さらには、高次自由度振動伝達特性でもよく、特に制限されないが、簡単な伝達特性で、後述する１００Ｈｚ以下の振動による性能評価を行うことができる点で、特性パラメータは、非減衰固有角振動数ωｎと減衰係数ζを含む、１自由度振動伝達特性であることが好ましい。 The vibration transmission characteristic H in the present embodiment may be a one-degree-of-freedom vibration transmission characteristic, a two-degree-of-freedom vibration transmission characteristic, or a higher-order degree-of-freedom vibration transmission characteristic, and is not particularly limited. The characteristic parameter is preferably a one-degree-of-freedom vibration transfer characteristic including a non-damped natural angular frequency ω n and a damping coefficient ζ in that performance evaluation can be performed with vibration of 100 Hz or less.
図３に示すステップＳ６０におけるタイヤの評価では、上述したデータ処理方法において算出した収束入力データ及び特性パラメータの収束値の少なくとも一方を用いて、評価部４０がタイヤＴの振動乗り心地性能を評価する。
例えば、１自由度振動伝達特性の場合、特性パラメータは、非減衰固有角振動数ωｎと減衰係数ζを含むので、評価部４０は、減衰係数ζの値からタイヤのサイド・ビード部の振動の吸収の程度を表すダンパ要素Ｃ（Ｎ・秒／ｍｍ）を評価することができる。
１自由度振動伝達特性の場合、非減衰固有角振動数ωｎと減衰係数ζは、１つの質量ｍ（ｋｇ）と１つのばね要素Ｋ（Ｎ／ｍｍ）と、１つのダンパ要素（Ｎ・秒／ｍｍ）とを用いて、以下のように表すことできる。
ωｎ ＝ （Ｋ／ｍ）（１／２）
ζ ＝ Ｃ／（２・（ｍ・Ｋ）（１／２））
したがって、ｍ＝Ｋ／ωｎ ２ 、Ｃ＝２Ｋζ／ωｎ と表される。
ここで、Ｋは、予めサイド・ビード部のばね要素Ｋであり、実験により算出される公知のばね定数に対応するので、実験から既知の値として求めておくことができる。したがって、Ｃ＝２Ｋζ／ωｎ を用いてダンパ要素Ｃを算出することができる。
なお、ばね要素Ｋの算出方法は、特開平０１−１５６６３４号公報あるいは特開平０１−１５６６３５号公報に開示されている。
In the tire evaluation in step S60 shown in FIG. 3, the
For example, in the case of a one-degree-of-freedom vibration transfer characteristic, since the characteristic parameters include the non-damped natural angular frequency ω n and the damping coefficient ζ, the
In the case of the one-degree-of-freedom vibration transfer characteristic, the undamped natural angular frequency ω n and the damping coefficient ζ are one mass m (kg), one spring element K (N / mm), and one damper element (N · Second / mm) and can be expressed as follows.
ω n = (K / m) (1/2)
ζ = C / (2 · (m · K) (1/2) )
Therefore, m = K / ω n 2 and C = 2Kζ / ω n are expressed.
Here, K is a spring element K of the side bead portion, and corresponds to a known spring constant calculated by experiment, and can be obtained as a known value from the experiment. Therefore, it is possible to calculate the damper element C using C = 2Kζ / ω n.
The method for calculating the spring element K is disclosed in Japanese Patent Laid-Open No. 01-156634 or Japanese Patent Laid-Open No. 01-156635.
図７（ａ）は、Ｃ＝２Ｋζ／ωｎ を用いて算出した、構造の異なるタイヤＴＡ〜ＴＤのダンパ要素Ｃの算出結果の一例を示す図である。図７（ａ）からわかるように、タイヤＴＤのダンパ要素Ｃが最も高く、大きい順に、タイヤＴＡ、タイヤＴＢ、タイヤＴＣである。ダンパ要素Ｃの値は大きいほど、振動伝達特性において振動を減衰させる程度が高く、振動を素早く減衰させるので、タイヤ入力Ｉがない状態ではダンパ要素Ｃの値が大きいほど好ましい。したがって、突起２２ａ等の段差を通過した後のタイヤＴの振動乗り心地の性の評価では、ダンパ要素Ｃの値が大きいことが好ましい。
また、タイヤＴが突起２２ａを通過するとき、タイヤ入力Ｉによってコツゴツとした振動を受けるので、この振動の評価のためには、タイヤ入力Ｉの入力データを用いることが好ましい。図７（ｂ）は、タイヤＴＡ〜ＴＤのタイヤ入力Ｉの入力波形の一例を示す図である。図７（ｂ）によると、タイヤ構造が異なることにより、タイヤ入力Ｉの入力波形も変化することがわかる。したがって、この入力波形の大きさや時間変化の大きさによって、コツゴツとした振動の大小を評価することができる。
7 (a) is calculated using the C = 2Kζ / ω n, is a diagram illustrating an example of a calculation result of the damper element C of different tire TA~TD structurally. As can be seen from FIG. 7A, the damper element C of the tire TD is the highest, and the tire TA, the tire TB, and the tire TC are in descending order. The larger the value of the damper element C, the higher the degree of damping of the vibration in the vibration transfer characteristics and the quick damping of the vibration. Therefore, the larger the value of the damper element C is preferable in the absence of the tire input I. Therefore, it is preferable that the value of the damper element C is large in the evaluation of the vibration ride comfort of the tire T after passing through the steps such as the
Further, when the tire T passes through the
図８は、本実施形態の評価部４０が行う評価の一例を説明する図である。評価部４０は、２つの性能（性能１、性能２）の評価を行う。評価部４０は、性能１について、タイヤ入力Ｉの収束入力データの積分値、すなわち、タイヤＴが突起２２ａを通過する期間の積分値（積分１）によって評価する。評価部４０は、性能２について、タイヤ軸力Ｏの計測データの積分値（積分２）から、タイヤ入力Ｉの収束入力データの積分値（積分１）を差し引いた値によって評価する。なお、積分２の積分範囲については、計測データの値が所定のレベル以上の期間を範囲とする。
すなわち、評価部４０は、積分１によって性能１を評価し、積分２−積分１によって性能２を評価する。
FIG. 8 is a diagram illustrating an example of evaluation performed by the
That is, the
図９（ａ）は、図７（ｂ）に示すタイヤＴＡ〜ＴＤの収束入力データを重ね書きした図である。図９（ｂ）は、タイヤＴＡ〜ＴＤの積分１の値と、性能１のドライバによる官能評価結果を示す図である。図９（ｂ）によると、積分１の大きさに関して、タイヤＴＡ＞タイヤＴＤ＞タイヤＴＢ＞タイヤＴＣの順番になっている。ドライバによる官能評価結果では、タイヤＴＢ＝タイヤＴＣ＞タイヤＴＤ＞タイヤＴＡとなっている。これより、性能１は、積分１の値が小さいほど性能１の評価が高いことがわかる。
したがって、評価部４０は、積分１の値が小さいものほど、性能１が優れているタイヤであると評価する。
FIG. 9A is a diagram in which the convergence input data of the tires TA to TD shown in FIG. FIG. 9B is a diagram illustrating a value of integral 1 of tires TA to TD and a sensory evaluation result by a driver of performance 1. According to FIG. 9B, regarding the magnitude of the integral 1, the order is tire TA> tire TD> tire TB> tire TC. In the sensory evaluation result by the driver, tire TB = tire TC> tire TD> tire TA. From this, it can be understood that the performance 1 has a higher evaluation of the performance 1 as the value of the integral 1 is smaller.
Therefore, the
本実施形態では、積分１を用いて性能１の評価を行うが、タイヤＴが突起２２を通過する期間におけるタイヤ入力Ｉの収束入力データの時間微分の絶対値の平均値を用いて、性能１の評価を行うことも好ましい。図１０（ａ）は、図９（ａ）に示す収束入力データを時間微分した波形を示す図である。図１０（ｂ）は、上記時間微分した波形の、タイヤＴが突起２２を通過する期間における絶対値の平均値（“傾斜平均値”）と、性能１のドライバによる官能評価結果を示す図である。図１０（ｂ）によると、平均値の大きさに関して、タイヤＴＡ＞タイヤＴＤ＞タイヤＴＢ＞タイヤＴＣの順番になっている。ドライバによる官能評価結果では、タイヤＴＢ＝タイヤＴＣ＞タイヤＴＤ＞タイヤＴＡとなっている。これより、性能１は、上記平均値（“傾斜平均値”）の値が小さいほど性能１の評価が高いことがわかる。
したがって、評価部４０は、タイヤＴが突起２２を通過する期間におけるタイヤ入力Ｉの収束入力データの時間微分の絶対値の平均値が小さいものほど、性能１が優れているタイヤであると評価することも好ましい。
In this embodiment, the performance 1 is evaluated using the integral 1, but the performance 1 is calculated using the average value of the absolute values of the time differentiation of the convergence input data of the tire input I during the period in which the tire T passes the
Therefore, the
また、評価部４０は、性能２に関して、以下の処理を行ってタイヤＴを評価することが好ましい。
タイヤＴのタイヤ軸力Ｏの計測データの取得に用いた計測転動速度の条件と異なる評価用転動速度の条件で性能２の評価を行うために、まず、評価部４０は、計測に用いた計測転動速度におけるタイヤ入力Ｉの収束入力データの時間軸上の波形を、図１１に示すように、評価用転動速度に応じて伸縮することで得られる伸縮波形を表した伸縮収束入力データを作成する。図１１は、性能２の評価を行うために行う処理の一例を説明する図である。
Further, with respect to the
In order to evaluate the
次に、評価部４０は、特性パラメータの収束値を有する振動伝達特性Ｈを用いてコンボリューション処理を行い、このコンボリューション処理を行うことで算出された評価用転動速度に対応したタイヤ軸力Ｏの予測データの絶対値の積分値（積分２）と伸縮収束入力データの絶対値の積分値（積分１）と差分により、性能２の評価を行う。
図１２は、評価用転動速度を種々変化したときの積分値（積分２−積分１）の変化の一例を示す図である。評価用転動速度が約６０ｋｍ／時において、積分値（積分２−積分１）は低下し、評価用転動速度が約６０ｋｍより高くなると、積分値（積分２−積分１）は上昇する挙動を示す。
Next, the
FIG. 12 is a diagram illustrating an example of a change in the integral value (integral 2−integral 1) when the rolling speed for evaluation is variously changed. When the evaluation rolling speed is about 60 km / hour, the integral value (integral 2 -integral 1) decreases, and when the evaluation rolling speed becomes higher than about 60 km, the integral value (integral 2 -integral 1) increases. Indicates.
図１３は、タイヤＴＡ〜タイヤＴＤの積分値（積分２−積分１）と、性能２のドライバによる官能評価結果の一例を示す図である。性能２は、特に６０〜８０ｋｍ／時の条件で発生し易い振動として知られている。したがって、性能２に関して、６０〜８０ｋｍ／時におけるタイヤＴＡ〜タイヤＴＤの積分値（積分２−積分１）を注目すると、積分値（積分２−積分１）の大きさに関して、タイヤＴＣ＞タイヤＴＢ＞タイヤＴＡ＞タイヤＴＤの順番になっている。ドライバによる性能２の官能評価結果では、タイヤＴＡ＝タイヤＴＤ＞タイヤＴＢ＞タイヤＴＣとなっている。これより、性能２は、上記積分値（積分２−積分１）が小さいほど性能２の評価が高いことがわかる。
したがって、評価部４０は、上記積分値（積分２−積分１）が小さいものほど、性能２が優れているタイヤであると評価することが好ましい。
FIG. 13 is a diagram illustrating an example of an integrated value (integral 2−integral 1) of tire TA to tire TD and a sensory evaluation result by a driver of
Therefore, it is preferable that the
性能２は、図８に示すように、積分１の期間を除去しているので、突起２２ａを通過後の振動の大小に依存する。この場合の振動の大小は、上述したようにダンパ要素Ｃの大きさに依存する。このため、性能２に関しては、ダンパ要素Ｃが大きいほど、性能２が優れているともいえる。この点で、ダンパ要素Ｃを用いて性能２の評価を行うこともできる。
As shown in FIG. 8, the
また、本実施形態では、図３に示すステップ４０を満足するタイヤ入力Ｉの収束入力データ及び特性パラメータの収束値を、特性パラメータで表したタイヤの振動伝達特性Ｈを、車両のサスペンションをモデル化したサスペンションモデルに組み合わせた統合モデルに入力することにより、車両の振動乗り心地性能の評価を行うことも好ましい。これにより、タイヤＣのダンパ要素Ｃを含んだ形態で車両の振動乗り心地性能の評価を行うことができる。
Further, in the present embodiment, the tire input I satisfying
以上、本発明のデータ処理方法、タイヤの評価方法、車両の振動乗り心地性能の評価方法、及びデータ処理装置について詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良や変更をしてもよいのはもちろんである。 The data processing method, the tire evaluation method, the vehicle vibration ride comfort performance evaluation method, and the data processing device of the present invention have been described in detail above. However, the present invention is not limited to the above-described embodiment, and the gist of the present invention. It goes without saying that various improvements and changes may be made without departing from the scope of the invention.
１０ データ処理装置
１２ ＣＰＵ
１４ メモリ
１６ ディスプレイ
１８ 入力操作系
２０ タイヤドラム試験機
２２ 試験ドラム
２２ａ 突起
２４ タイヤ回転軸
３０ データ取得部
３２ 入力データ算出部
３４ パラメータ算出部
３６ 制御部
３８ 出力部
４０ 評価部
10
14
Claims (11)
前記コンピュータが、前記タイヤ軸力Ｏの計測データから、設定された前記振動伝達特性Ｈの特性パラメータの仮値を用いて、前記タイヤ入力Ｉの仮入力データを算出する入力データ算出ステップと、
前記コンピュータが、前記入力データ算出ステップで算出した前記仮入力データのうち、前記タイヤが前記段差を通過する前後のデータの値を０にした修正仮入力データと前記タイヤ軸力Ｏの計測データから、前記振動伝達特性Ｈの特性パラメータの仮値を修正した修正仮値を算出するパラメータ値算出ステップと、を備え、
前記仮入力データのうち前記タイヤが前記段差を通過する前後のデータの値が０に収束する条件を少なくとも含む収束条件を満足するまで、前記コンピュータは、前記修正仮値を、前記特性パラメータの前記仮値として設定して前記入力データ算出ステップ及び前記パラメータ値算出ステップを繰り返し行うことにより、前記収束条件を満足するときの前記タイヤ入力Ｉの仮入力データと前記特性パラメータの仮値のそれぞれを前記タイヤ入力Ｉの収束入力データ及び前記特性パラメータの収束値として算出する、ことを特徴とするデータ処理方法。 The tire input I received by the tire from the step, the tire axial force O, and the tire from the vibration of the tire axial force O generated on the tire rotation shaft when the tire rolling on the road surface passes through the step on the road surface. A data processing method in which a computer calculates a vibration transfer characteristic H between tire inputs,
An input data calculating step in which the computer calculates temporary input data of the tire input I from the measurement data of the tire axial force O using the set temporary value of the characteristic parameter of the vibration transfer characteristic H;
From the temporary input data calculated by the computer in the input data calculation step, corrected temporary input data in which the value of data before and after the tire passes the step is set to 0, and measurement data of the tire axial force O A parameter value calculating step of calculating a corrected temporary value obtained by correcting the temporary value of the characteristic parameter of the vibration transfer characteristic H,
Until the computer satisfies the convergence condition including at least the condition that the value of the data before and after the tire passes through the step in the temporary input data satisfies the convergence condition, the computer sets the corrected temporary value as the characteristic parameter. By setting the provisional value and repeating the input data calculation step and the parameter value calculation step, each of the provisional input data of the tire input I and the provisional value of the characteristic parameter when the convergence condition is satisfied is obtained. A data processing method characterized by calculating the convergence input data of the tire input I and the convergence value of the characteristic parameter.
前記収束入力データ及び前記特性パラメータの前記収束値の少なくとも一方を用いて、前記コンピュータが前記タイヤの振動乗り心地性能を評価するステップと、を備えることを特徴とするタイヤの評価方法。 A step of calculating a convergence value of the convergence input data and the characteristic parameter in the data processing method according to any one of claims 1 to 4,
And a step of evaluating the vibration ride performance of the tire using at least one of the convergence input data and the convergence value of the characteristic parameter.
前記特性パラメータの前記収束値のうち、前記非減衰固有角振動数及び前記減衰係数の収束値からタイヤのサイド・ビード部の振動の吸収の程度を表すダンパ要素を前記コンピュータが評価するステップと、を備えることを特徴とするタイヤの評価方法。 A step of calculating a convergence value of the convergence input data and the characteristic parameter in the data processing method according to claim 4,
The computer evaluating a damper element representing a degree of vibration absorption of a tire side bead from the convergence value of the non-damped natural angular frequency and the damping coefficient among the convergence value of the characteristic parameter; A method for evaluating a tire, comprising:
前記収束入力データの値の絶対値の積分値、あるいは、前記タイヤが前記段差を通過する期間における前記収束入力データの時間微分の絶対値の平均値を用いて、前記タイヤの振動乗り心地性能のうちの性能１の評価を行う、請求項５に記載のタイヤの評価方法。 In the step of evaluating the vibration riding comfort performance of the tire,
Using the integral value of the absolute value of the convergent input data or the average value of the absolute value of the time derivative of the convergent input data during the period when the tire passes the step, the vibration ride comfort performance of the tire The tire evaluation method according to claim 5, wherein performance 1 is evaluated.
前記タイヤの転動する計測用転動速度と異なる評価用転動速度における前記タイヤの振動乗り心地性能のうちの性能２の評価を行うために、前記計測用転動速度における前記収束入力データの時間軸上の波形を前記評価用転動速度に応じて伸縮することで得られる伸縮波形を表した伸縮収束入力データと、前記特性パラメータの前記収束値を有する振動伝達特性を用いてコンボリューション処理を行い、該コンボリューション処理を行うことで算出された前記評価用転動速度に対応した前記タイヤ軸力Ｏの予測データの絶対値の積分値と前記伸縮収束入力データの絶対値の積分値との差分により、前記性能２の評価を行う、請求項５に記載のタイヤの評価方法。 In the step of evaluating the vibration riding comfort performance of the tire,
In order to evaluate performance 2 out of the vibration riding comfort performance of the tire at an evaluation rolling speed different from the rolling speed for measurement at which the tire rolls, the convergence input data at the measurement rolling speed is Convolution processing using expansion / contraction convergence input data representing the expansion / contraction waveform obtained by expanding / contracting the waveform on the time axis according to the evaluation rolling speed, and the vibration transfer characteristic having the convergence value of the characteristic parameter And the integral value of the absolute value of the prediction data of the tire axial force O corresponding to the rolling speed for evaluation calculated by performing the convolution process and the integral value of the absolute value of the expansion / contraction convergence input data, The tire evaluation method according to claim 5, wherein the performance 2 is evaluated based on the difference between the two.
前記特性パラメータで表した前記タイヤの振動伝達特性を、車両のサスペンションをモデル化したサスペンションモデルに組み合わせた統合モデルに、前記収束入力データを入力することにより、車両の振動乗り心地性能の評価を行うステップと、を備えることを特徴とする車両の振動乗り心地性能の評価方法。 A step of calculating a convergence value of the convergence input data and the characteristic parameter in the data processing method according to any one of claims 1 to 4,
Evaluation of the vehicle's vibration ride comfort performance is performed by inputting the convergence input data into an integrated model in which the vibration transmission characteristics of the tire represented by the characteristic parameters are combined with a suspension model obtained by modeling a vehicle suspension. And a method for evaluating vibration ride comfort performance of a vehicle.
前記タイヤ軸力Ｏの計測データを取得するデータ取得部と、
前記計測データから、設定された前記振動伝達特性Ｈの特性パラメータの仮値を用いて、前記タイヤ入力Ｉの仮入力データを算出する入力データ算出部と、
前記入力データ算出部で算出した前記仮入力データのうち、前記タイヤが前記段差を通過する前後のデータの値を０にした修正仮入力データと前記タイヤ軸力Ｏの計測データから、前記振動伝達特性の特性パラメータの仮値を修正した修正仮値を算出するパラメータ値算出部と、
前記仮入力データのうち前記タイヤが前記段差を通過する前後のデータの値が０に収束する条件を少なくとも含む収束条件を満足するまで、前記修正仮値を、前記特性パラメータの前記仮値として設定して前記仮入力データの算出及び前記修正仮値の算出を繰り返し行うよう前記入力データ算出部及び前記パラメータ値算出部を制御する制御部と、
前記収束条件を満足するときの前記タイヤ入力Ｉの仮入力データと前記特性パラメータの仮値のそれぞれを前記タイヤ入力Ｉの収束入力データ及び前記特性パラメータの収束値として算出する出力部と、を備えることを特徴とするデータ処理装置。 The tire input I received by the tire from the step, the tire axial force O, and the tire from the vibration of the tire axial force O generated on the tire rotation shaft when the tire rolling on the road surface passes through the step on the road surface. A data processing device for calculating a vibration transfer characteristic H between tire inputs,
A data acquisition unit for acquiring measurement data of the tire axial force O;
An input data calculation unit that calculates temporary input data of the tire input I using a temporary value of the set characteristic parameter of the vibration transfer characteristic H from the measurement data;
Of the temporary input data calculated by the input data calculation unit, the vibration transmission is based on the corrected temporary input data in which the value of the data before and after the tire passes through the step is 0 and the measurement data of the tire axial force O. A parameter value calculation unit for calculating a corrected temporary value obtained by correcting the temporary value of the characteristic parameter of the characteristic;
The corrected provisional value is set as the provisional value of the characteristic parameter until a convergence condition including at least a condition in which the value of the data before and after the tire passes through the step of the provisional input data converges to 0 is satisfied. A control unit that controls the input data calculation unit and the parameter value calculation unit so as to repeatedly calculate the temporary input data and the correction temporary value;
An output unit that calculates the temporary input data of the tire input I and the temporary value of the characteristic parameter when the convergence condition is satisfied as the convergence input data of the tire input I and the convergence value of the characteristic parameter; A data processing apparatus.
さらに、前記特性パラメータの前記収束値のうち、前記減衰係数の収束値からタイヤのサイド・ビード部の振動の吸収の程度を表すダンパ要素を前記コンピュータが評価する評価部を備える、請求項１０に記載のデータ処理装置。 The vibration transfer characteristic is a one-degree-of-freedom vibration transfer characteristic, and the characteristic parameters include a non-damped natural angular frequency and a damping coefficient,
Furthermore, the evaluation part which the said computer evaluates the damper element showing the degree of absorption of the vibration of the side bead part of a tire from the convergence value of the above-mentioned attenuation coefficient among the convergence values of the above-mentioned characteristic parameter is provided. The data processing apparatus described.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016239574A JP6805777B2 (en) | 2016-12-09 | 2016-12-09 | Data processing method, tire evaluation method, vehicle vibration ride comfort evaluation method, and data processing device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016239574A JP6805777B2 (en) | 2016-12-09 | 2016-12-09 | Data processing method, tire evaluation method, vehicle vibration ride comfort evaluation method, and data processing device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2018096753A true JP2018096753A (en) | 2018-06-21 |
| JP6805777B2 JP6805777B2 (en) | 2020-12-23 |
Family
ID=62631367
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016239574A Active JP6805777B2 (en) | 2016-12-09 | 2016-12-09 | Data processing method, tire evaluation method, vehicle vibration ride comfort evaluation method, and data processing device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6805777B2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2021085824A (en) * | 2019-11-29 | 2021-06-03 | 富士通株式会社 | Characteristic calculation program, characteristic calculation method, and information processing system |
| JP7465180B2 (en) | 2020-08-28 | 2024-04-10 | Ｔｏｙｏ Ｔｉｒｅ株式会社 | Tire characteristic evaluation method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11352024A (en) * | 1998-06-04 | 1999-12-24 | Bridgestone Corp | Method and device for predicting higher-order component of high-speed uniformity of tire and method for manufacturing the tire |
| JP2000241309A (en) * | 1999-02-23 | 2000-09-08 | Bridgestone Corp | Method for simulating vibration and noise of tire |
| JP2004198219A (en) * | 2002-12-18 | 2004-07-15 | Yokohama Rubber Co Ltd:The | Testing device for tire vibrational characteristic |
| JP2006242581A (en) * | 2005-02-28 | 2006-09-14 | Bridgestone Corp | Tire uniformity correction method |
| JP2006292547A (en) * | 2005-04-11 | 2006-10-26 | Sumitomo Rubber Ind Ltd | Measuring technique of vibration characteristic of tire and measuring device of vibration characteristic of tire |
-
2016
- 2016-12-09 JP JP2016239574A patent/JP6805777B2/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11352024A (en) * | 1998-06-04 | 1999-12-24 | Bridgestone Corp | Method and device for predicting higher-order component of high-speed uniformity of tire and method for manufacturing the tire |
| JP2000241309A (en) * | 1999-02-23 | 2000-09-08 | Bridgestone Corp | Method for simulating vibration and noise of tire |
| JP2004198219A (en) * | 2002-12-18 | 2004-07-15 | Yokohama Rubber Co Ltd:The | Testing device for tire vibrational characteristic |
| JP2006242581A (en) * | 2005-02-28 | 2006-09-14 | Bridgestone Corp | Tire uniformity correction method |
| JP2006292547A (en) * | 2005-04-11 | 2006-10-26 | Sumitomo Rubber Ind Ltd | Measuring technique of vibration characteristic of tire and measuring device of vibration characteristic of tire |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2021085824A (en) * | 2019-11-29 | 2021-06-03 | 富士通株式会社 | Characteristic calculation program, characteristic calculation method, and information processing system |
| JP7465180B2 (en) | 2020-08-28 | 2024-04-10 | Ｔｏｙｏ Ｔｉｒｅ株式会社 | Tire characteristic evaluation method |
Also Published As
| Publication number | Publication date |
|---|---|
| JP6805777B2 (en) | 2020-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4199273B2 (en) | Vehicle state estimation device | |
| ITMI952399A1 (en) | CHECK METHOD FOR DETERMINING THE DEGREE OF COMFORT OF A TIRE OF A WHEEL PER VEHICLE | |
| JP2018096753A (en) | Data processing method, estimation method of tire, evaluation method of vibration ride comfort performance, and data processing device | |
| JP7390731B2 (en) | Road surface profile estimation device, road surface profile estimation system, road surface profile estimation method, and road surface profile estimation program | |
| JP2011076240A (en) | Planning support system for vehicle | |
| Kido et al. | Tire and road input modeling for low-frequency road noise prediction | |
| JP4861897B2 (en) | Vehicle state estimation device | |
| Papadrakakis et al. | Tmeasy 6.0-a handling tire model that incorporates the first two belt eigenmodes | |
| Yun et al. | Modal identification of building structures under unknown input conditions using extended Kalman filter and long-short term memory | |
| JP6952531B2 (en) | Railway vehicle abnormality diagnosis system | |
| JPH11264785A (en) | Method and device for inspecting installed shock absorber | |
| Patten et al. | Validation of a nonlinear automotive seat cushion vibration model | |
| Kindt et al. | Measurement of the tire dynamic transfer stiffness at operational excitation levels | |
| JP2004340849A (en) | Tire model generating method, and tire performance predicting method and program | |
| JP2010261770A (en) | Accumulated fatigue calculation method, and vibration test method using the same | |
| JP6346526B2 (en) | Deformation simulation method for tire and spring assembly | |
| Baecker et al. | From road excitation to spindle forces in frequency domain: Linearization of the rolling tire | |
| Ortega Almirón et al. | Structure-Borne Prediction on a Tire-Suspension Assembly Using Experimental Invariant Spindle Forces | |
| JP4201821B2 (en) | Tire model determination method, tire transient response data calculation method, tire evaluation method, and tire design method | |
| JP6138631B2 (en) | Apparatus, method and computer program for predicting tire vibration | |
| JPH11264784A (en) | Method and device for inspecting installed shock absorber | |
| Saraf et al. | Integration of real and virtual tools for suspension development | |
| JP2022037666A (en) | Tire simulation method | |
| Boere | Prediction of road texture influence on rolling resistance | |
| JP7074579B2 (en) | Tire performance simulation method and equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20191204 |
|
| TRDD | Decision of grant or rejection written | ||
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20201023 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20201104 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20201117 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6805777 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |