JP5533277B2 - Evaluation method of tire inclusion jumping characteristic, tire inclusion jumping characteristic evaluation apparatus and program - Google Patents
Evaluation method of tire inclusion jumping characteristic, tire inclusion jumping characteristic evaluation apparatus and program Download PDFInfo
- Publication number
- JP5533277B2 JP5533277B2 JP2010125706A JP2010125706A JP5533277B2 JP 5533277 B2 JP5533277 B2 JP 5533277B2 JP 2010125706 A JP2010125706 A JP 2010125706A JP 2010125706 A JP2010125706 A JP 2010125706A JP 5533277 B2 JP5533277 B2 JP 5533277B2
- Authority
- JP
- Japan
- Prior art keywords
- model
- tire
- road surface
- inclusion
- particle
- 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.)
- Expired - Fee Related
Links
Landscapes
- Tires In General (AREA)
Description
本発明は、路面に設けられる介在物上をタイヤが通過するときのタイヤの介在物飛び跳ね特性を評価する方法および装置、さらには、この方法を実行するプログラムに関する。 The present invention relates to a method and an apparatus for evaluating the inclusion jumping characteristic of a tire when the tire passes over an inclusion provided on a road surface, and a program for executing the method.
近年、コンピュータの処理速度の向上に伴って、車両に装着される空気入りタイヤ(以降、単にタイヤという)のウェット性能を、有限要素法や有限体積法を用いて予測する方法が種々提案されている。
タイヤのウェット性能とは、例えば、ハイドロプレーニング現象に代表されるように、タイヤと路面間に水が介在することによってタイヤ特性が低下することをいう。従って、タイヤのウェット性能が良好であるとは、タイヤと路面間に水が介在することなく、たとえ水が介在しても、その介在を可能な限り抑えられ、また、この水の介在による影響が小さいことをいう。
In recent years, various methods for predicting the wet performance of a pneumatic tire (hereinafter simply referred to as a tire) mounted on a vehicle by using a finite element method or a finite volume method have been proposed along with an improvement in computer processing speed. Yes.
The wet performance of a tire means that tire characteristics are deteriorated due to water intervening between the tire and a road surface, for example, as represented by a hydroplaning phenomenon. Therefore, good wet performance of the tire means that water does not intervene between the tire and the road surface, even if water intervenes, the interposition is suppressed as much as possible, and the influence of this water interposition Is small.
一方、路面上に形成された水膜をタイヤが通過するとき、水膜の一部が路面から跳ね上がり、水が飛散する場合がある。このような水跳ねは、濡れた路面を車両が通過したとき引き起こされ、道路の側方にいる人の通行を妨げる場合がある。あるいは、高速道路を走行するトラックやバスに対して隣接する車線を走行中の車両がトラックやバスの水跳ねを被って走行の障害になる場合もある。 On the other hand, when the tire passes through the water film formed on the road surface, a part of the water film may jump from the road surface and the water may be scattered. Such water splashes are caused when the vehicle passes on a wet road surface, and may hinder traffic on the side of the road. Alternatively, a vehicle traveling on a lane adjacent to a truck or bus traveling on a highway may be subject to traveling obstacles due to the water splash of the truck or bus.
このような状況下、タイヤのウェット性能等の路面状態に応じたタイヤ特性を予測、評価する際、膨大なモデルとならず、従来に比べて精度の高い演算結果を得て評価することのできるタイヤ特性予測方法が知られている(特許文献1)。 Under such circumstances, when predicting and evaluating the tire characteristics according to the road surface condition such as the wet performance of the tire, it is not an enormous model, and it is possible to obtain and evaluate a calculation result with higher accuracy than the conventional model. A tire characteristic prediction method is known (Patent Document 1).
当該文献1では、有限要素モデルのタイヤモデルと路面モデルと、この路面モデル上の少なくとも一部分に複数の粒子モデルを配した路面状態再現モデルとが作成される。このタイヤモデルが粒子モデルの配された路面モデル上の領域を踏み込む際、タイヤモデルおよび路面状態再現モデルの変形計算を行う時間ステップ毎にタイヤモデルと粒子モデルとの相互作用の演算を行いながらタイヤモデルおよび路面状態再現モデルの変形計算を行う。この計算結果から、タイヤモデルまたは路面状態再現モデルの所定の特性物理量を算出してタイヤ特性を予測し評価する。
当該文献1では、特性物理量として、路面状態再現モデルがタイヤモデルに作用する揚力、路面モデルがタイヤモデルに作用する路面反力、および路面状態再現モデルにおける粒子モデルの移動位置や流速を算出できるが、タイヤの水跳ね特性を評価することはできない。
In Document 1, a tire model and a road surface model of a finite element model, and a road surface state reproduction model in which a plurality of particle models are arranged on at least a part of the road surface model are created. When this tire model steps into the area on the road surface model where the particle model is placed, the tire model and the particle surface model are calculated while calculating the interaction between the tire model and the particle model at each time step when the tire model and the road surface state reproduction model are calculated. Deformation calculation of model and road surface state reproduction model is performed. From this calculation result, a predetermined characteristic physical quantity of the tire model or the road surface state reproduction model is calculated to predict and evaluate the tire characteristics.
In the document 1, as the characteristic physical quantity, the lift force that the road surface state reproduction model acts on the tire model, the road surface reaction force that the road surface model acts on the tire model, and the movement position and flow velocity of the particle model in the road surface state reproduction model can be calculated. The water splash characteristics of a tire cannot be evaluated.
一般に、タイヤを装着した車両を用いて実際の水跳ねを行って、タイヤの水はね特性を評価することは可能であるが、定量的な評価は難しい。また、実際のタイヤを作製して車両による試験を行う煩雑な作業を要する。このため、タイヤの水跳ね特性を効率よく、定量的に評価することはできない。このような問題は、水跳ねに限られず、雪、砂、泥、および砂利等の跳ね上げについても同様に生じる。 In general, it is possible to evaluate water splash characteristics of a tire by performing actual water splash using a vehicle equipped with a tire, but quantitative evaluation is difficult. Moreover, the complicated work which produces an actual tire and performs the test by a vehicle is required. For this reason, the water splash characteristic of a tire cannot be evaluated efficiently and quantitatively. Such a problem is not limited to water splashing, but also occurs for snow, sand, mud, gravel and the like.
そこで、本発明は、路面に設けられる介在物上をタイヤが通過するときのタイヤの介在物飛び跳ね特性を効率よく、定量的に評価することができる、タイヤの介在物飛び跳ね特性の評価方法、タイヤの介在物飛び跳ね特性の評価装置およびプログラムを提供することを目的とする。 Accordingly, the present invention provides a method for evaluating inclusion jumping characteristics of a tire, a tire that can efficiently and quantitatively evaluate the inclusion jumping characteristics of a tire when the tire passes over inclusions provided on a road surface. An object of the present invention is to provide an evaluation device and program for the jumping property of inclusions.
本発明の一態様は、路面に設けられる介在物上をタイヤが通過するときのタイヤの介在物の飛び跳ね特性を評価する方法であって、
タイヤを再現したタイヤモデルと、路面を再現した路面モデルと、路面上に設けられる介在物を再現した、複数の微小モデルがお互いに分離可能に含まれる介在物モデルとを作成するステップと、
前記介在物モデルが設けられた前記路面モデル上で前記タイヤモデルを転動させることにより、前記微小モデルを飛散させるステップと、
前記飛散した前記微小モデルを、設定された時間間隔毎に、前記路面モデルの面に対して垂直な平面に投影するステップと、
前記平面に投影された前記微小モデルの時系列の投影像を用いて、前記投影像毎に、飛散する前記微小モデルの前記投影像の包絡線を求め、前記投影像毎に求めた前記包絡線の中から、前記包絡線の形状が最も大きい包絡線を用いてタイヤの介在物の飛び跳ね特性を評価するステップと、を有することを特徴とする。
One aspect of the present invention is a method for evaluating the jumping characteristics of tire inclusions when the tire passes over inclusions provided on a road surface,
Creating a tire model that reproduces a tire, a road surface model that reproduces a road surface, and an inclusion model that reproduces inclusions provided on the road surface and includes a plurality of minute models that are separable from each other;
Scattering the micro model by rolling the tire model on the road surface model provided with the inclusion model;
Projecting the scattered micro model on a plane perpendicular to the plane of the road surface model at set time intervals ;
Using the time-series projection image of the micro model projected on the plane, an envelope of the projection image of the micro model scattered is obtained for each projection image, and the envelope obtained for each projection image And the step of evaluating the jumping characteristics of the inclusions of the tire using the envelope having the largest shape of the envelope .
本発明の別の一態様は、路面に設けられる介在物上をタイヤが通過するときのタイヤの介在物の飛び跳ね特性を評価する装置であって、
タイヤを再現したタイヤモデルと、路面を再現した路面モデルと、路面上に設けられる介在物を再現した、複数の微小モデルがお互いに分離可能に含まれる介在物モデルとを作成するモデル作成部と、
前記介在物モデルが設けられた前記路面モデル上で前記タイヤモデルを転動させることにより、前記微小モデルを飛散させるシミュレーション演算部と、
前記飛散した前記微小モデルを、設定された時間間隔毎に、前記路面モデルの面に対して垂直な平面に投影するシミュレーション結果処理部と、
前記平面に投影された前記微小モデルの時系列の投影像を用いて、前記投影像毎に、飛散する前記微小モデルの前記投影像の包絡線を求め、前記投影像毎に求めた前記包絡線の中から、前記包絡線の形状が最も大きい包絡線を用いて、タイヤの介在物の飛び跳ね特性を評価する評価部と、を有することを特徴とする。
Another aspect of the present invention is an apparatus for evaluating the jumping characteristics of tire inclusions when the tire passes over inclusions provided on a road surface,
A model creation unit that creates a tire model that reproduces a tire, a road surface model that reproduces a road surface, and an inclusion model that reproduces inclusions provided on the road surface and includes a plurality of minute models that are separable from each other; ,
A simulation calculation unit that scatters the minute model by rolling the tire model on the road surface model provided with the inclusion model;
A simulation result processing unit that projects the scattered micro model on a plane perpendicular to the plane of the road surface model for each set time interval ;
Using the time-series projection image of the micro model projected on the plane, an envelope of the projection image of the micro model scattered is obtained for each projection image, and the envelope obtained for each projection image And an evaluation unit that evaluates the jumping characteristic of the inclusions of the tire using the envelope having the largest shape of the envelope .
本発明のさらに別の一態様は、路面に設けられる介在物上をタイヤが通過するときの、タイヤの介在物飛び跳ね特性をコンピュータに評価させるプログラムであって、
タイヤを再現したタイヤモデルと、路面を再現した路面モデルと、路面上に設けられる介在物を再現した、複数の微小モデルがお互いに分離可能に含まれる介在物モデルとをコンピュータが作成する手順と、
前記介在物モデルが設けられた前記路面モデル上で前記タイヤモデルを転動させることにより、コンピュータが前記介在物モデルの前記微小モデルを飛散させる手順と、
前記飛散した前記微小モデルを、コンピュータが、設定された時間間隔毎に、前記路面モデルの面に対して垂直な平面に投影する手順と、
前記平面に投影された前記微小モデルの時系列の投影像を用いて、前記投影像毎に、飛散する前記微小モデルの前記投影像の包絡線を求め、前記投影像毎に求めた前記包絡線の中から、前記包絡線の形状が最も大きい包絡線を用いて、コンピュータが、タイヤの介在物の飛び跳ね特性を評価する手順と、を有する。
Yet another aspect of the present invention is a program that causes a computer to evaluate inclusion jumping characteristics of a tire when the tire passes over inclusions provided on a road surface.
A procedure for a computer to create a tire model that reproduces a tire, a road surface model that reproduces a road surface, and an inclusion model that reproduces inclusions provided on the road surface and includes a plurality of minute models that can be separated from each other. ,
A step of causing the computer to scatter the minute model of the inclusion model by rolling the tire model on the road surface model provided with the inclusion model;
The computer projects the scattered micro model on a plane perpendicular to the plane of the road surface model at set time intervals , and
Using the time-series projection image of the micro model projected on the plane, an envelope of the projection image of the micro model scattered is obtained for each projection image, and the envelope obtained for each projection image The computer has a procedure for evaluating the jumping characteristics of the inclusions of the tire using the envelope having the largest shape of the envelope .
上記タイヤの介在物飛び跳ね特性の評価方法、タイヤの介在物飛び跳ね特性の評価装置およびプログラムは、路面に設けられる介在物上をタイヤが通過するときのタイヤの介在物飛び跳ね特性を効率よく、定量的に評価することができる。 The method for evaluating the inclusion jumping characteristic of the tire, the apparatus for evaluating the inclusion jumping characteristic of the tire, and the program efficiently and quantitatively determine the inclusion jumping characteristic of the tire when the tire passes over the inclusion provided on the road surface. Can be evaluated.
以下、本発明のタイヤの介在物飛び跳ね特性の評価方法、タイヤの介在物飛び跳ね特性の評価装置およびプログラムについて詳細に説明する。 Hereinafter, the evaluation method of the inclusion jumping characteristic of the tire, the evaluation device for the inclusion jumping characteristic of the tire, and the program according to the present invention will be described in detail.
(タイヤ評価装置)
図1は、本実施形態のタイヤの水跳ね特性を評価する評価装置10の概略の構成を示す図である。水跳ね特性の評価とは、タイヤが路面上の水膜を通過するとき、タイヤ周囲に水が跳ね上がるときの、到達距離および到達高さ、あるいは水跳ね密度等を評価することをいう。
(Tire evaluation device)
FIG. 1 is a diagram illustrating a schematic configuration of an evaluation apparatus 10 that evaluates the water splash characteristics of a tire according to the present embodiment. Evaluation of water splash characteristics refers to evaluation of reach distance and reach height, or water splash density when water splashes around the tire when the tire passes a water film on the road surface.
評価装置10は、タイヤの水跳ね特性を、タイヤモデルを用いた転動シミュレーションを用いて評価する。評価装置10は、作成されたタイヤモデルを、別途作成された水膜モデルが設けられた路面モデル上に転動させる転動シミュレーションを行う。水膜モデルには、複数の粒子モデルがお互いに分離可能に含まれる。タイヤモデルが水膜モデルを通過するとき、水膜モデル中の粒子モデルがお互いに分離して飛散するので、評価装置10は、この粒子モデルの飛散状態を、1つの平面に投影して、粒子モデルの投影像を用いて、タイヤの水跳ね特性を評価する。
評価装置10が水跳ね特性を評価するタイヤは、例えば乗用車用タイヤあるいは重荷重用(バス、トラック用)タイヤ等が含まれ、タイヤサイズに限定されない。
The evaluation device 10 evaluates the water splash characteristics of the tire using a rolling simulation using a tire model. The evaluation apparatus 10 performs a rolling simulation in which the created tire model is rolled on a road surface model provided with a separately created water film model. The water film model includes a plurality of particle models so as to be separable from each other. When the tire model passes through the water film model, since the particle models in the water film model are separated from each other and scattered, the evaluation apparatus 10 projects the scattering state of the particle model onto one plane to generate particles. The model's projected image is used to evaluate the water splash characteristics of the tire.
The tire for which the evaluation device 10 evaluates the water splash characteristic includes, for example, a passenger car tire or a heavy load (bus or truck) tire, and is not limited to the tire size.
評価装置10は、CPU12、バス14、メモリ16、および、入出力インターフェース部18を主に有するコンピュータを含む。評価装置10は、入出力インターフェース部18を通して、マウスやキーボード等の入力操作系20と、ディスプレイやプリンタ等の出力装置22と接続されている。評価装置10は、メモリ16に記憶されたプログラムを呼び出して、プログラムを実行することにより、モデル作成部24、シミュレーション演算部26、シミュレーション結果処理部28、評価部30が形成されて、処理モジュール群36が形成される。 The evaluation apparatus 10 includes a computer that mainly includes a CPU 12, a bus 14, a memory 16, and an input / output interface unit 18. The evaluation device 10 is connected to an input operation system 20 such as a mouse and a keyboard and an output device 22 such as a display and a printer through an input / output interface unit 18. The evaluation device 10 calls a program stored in the memory 16 and executes the program, thereby forming a model creation unit 24, a simulation calculation unit 26, a simulation result processing unit 28, and an evaluation unit 30, and a processing module group. 36 is formed.
モデル作成部24は、複数の要素によりタイヤを再現したタイヤモデルと、タイヤが接地する路面を再現した路面モデルと、路面上に設けられる水膜を再現した水膜モデルと、を作成する。水膜モデルは、複数の粒子モデルが分離可能に含まれる。本実施形態は介在物として水を用いるが、水の代わりに、雪、砂、泥、および砂利等が再現されてもよい。本実施形態では、水、すなわち水膜を介在物の代表として用いる。水膜モデルは、後述するように、粒子法を用いてモデル化された複数の粒子モデルが用いられる。この複数の粒子モデルは、タイヤモデルから受ける外力によってお互いに分離可能な程度の弱い拘束力を有する。 The model creation unit 24 creates a tire model that reproduces a tire by a plurality of elements, a road surface model that reproduces a road surface on which the tire contacts the ground, and a water film model that reproduces a water film provided on the road surface. The water film model includes a plurality of particle models in a separable manner. Although this embodiment uses water as inclusions, snow, sand, mud, gravel, or the like may be reproduced instead of water. In this embodiment, water, that is, a water film is used as a representative inclusion. As described later, a plurality of particle models modeled using the particle method are used as the water film model. The plurality of particle models have weak binding forces that can be separated from each other by an external force received from the tire model.
図2(a)は、タイヤモデルTM、路面モデルRMおよび水膜モデルSMを示す図である。
作成されるタイヤモデルTMは、図2(a)に示すような3次元形状をなした有限要素法に基づいたタイヤモデルである。図2に示すタイヤモデルTMでは要素分割するメッシュは表示されていない。
路面モデルRMは、図2(a)に示すような剛体の路面モデルRMあるいは、複数の要素で構成された有限要素モデルである。
水膜モデルSMは、図2(a)に示すように、路面モデルRMのタイヤモデルTMの走行部分の一部に形成されている。
FIG. 2A is a diagram showing a tire model T M , a road surface model R M, and a water film model S M.
The created tire model T M is a tire model based on the finite element method having a three-dimensional shape as shown in FIG. In the tire model T M shown in FIG. 2, the mesh for element division is not displayed.
The road surface model R M is a rigid road surface model R M as shown in FIG. 2A or a finite element model composed of a plurality of elements.
As shown in FIG. 2A, the water film model S M is formed in a part of the traveling portion of the tire model T M of the road surface model R M.
図2(b)は、水膜モデルSMを詳細に説明する図である。
路面モデルRM上の一部分に設けられる水膜モデルSMは、複数の粒子モデルPMを一定の間隔で配したモデルである。モデル作成部24は、水膜モデルSMの作成時、具体的には、水膜を再現した流体モデルとして、図2(b)に示すように、複数の粒子モデルPMをX方向、Y方向、Z方向に一定の間隔で等方状に配列する。モデル作成部24は、例えば、粒子モデルPMの間隔を2.5mm、粒子モデルの総数Nを数万個として、各粒子モデルPMの質量が同じになる様に、等方状に規則的に配置して、例えば、横幅400mm、厚さ10mm、長さ120mmとする水膜モデルSMを作成する。
水膜モデルSMの粒子モデルPMに適用される粒子法については、後述する。
FIG. 2B illustrates the water film model SM in detail.
The water film model S M provided in a part on the road surface model R M is a model in which a plurality of particle models P M are arranged at a constant interval. When creating the water film model S M , specifically, the model creation unit 24 converts a plurality of particle models P M in the X direction and Y as a fluid model reproducing the water film, as shown in FIG. Isotropically arranged at regular intervals in the Z direction. For example, the model creation unit 24 sets the interval between the particle models P M to 2.5 mm and the total number N of the particle models N to several tens of thousands, and is regular in an isotropic manner so that the mass of each particle model P M becomes the same. For example, a water film model S M having a lateral width of 400 mm, a thickness of 10 mm, and a length of 120 mm is created.
The particle method applied to the particle model P M of the water film model S M will be described later.
シミュレーション演算部26は、水膜モデルSMが設けられた路面モデルルRM上でタイヤモデルルTMを転動させることにより、水膜モデルルSMの粒子モデルルPMを飛散させる。粒子モデルルPMとタイヤモデル粒子モデルルTMとの相互作用については、後述する粒子法の説明の中で述べる。 Simulation unit 26 by the water film model S M is to roll the tire model Le T M on the road Moderuru R M provided, to scatter the particles Moderuru P M of the water film Moderuru S M. The interaction between the particle model P M and the tire model particle model T M will be described in the description of the particle method described later.
シミュレーション結果処理部28は、シミュレーション中、飛散した粒子モデルルPMを、路面モデルRMの面に対して垂直な平面に投影する処理を行う。
具体的には、シミュレーション結果処理部28は、タイヤモデルTMが、水膜モデルSM上の通過開始から通過終了までの所定の時間範囲を一定の時間間隔毎に、飛散する粒子モデルPMの像を上記平面に投影する。投影に用いる上記平面は、タイヤモデルTMの中心を通る路面モデルRMの面に対して垂直な面である。転動タイヤモデルTMは、路面モデルRMに対して移動するため、設定された時間毎に上記平面が設けられて、飛散する粒子モデルPMが上記平面に投影される。上記平面は、この平面の法線方向がタイヤモデルTMの移動方向となるように設けられる。この他に、図3に示されるように、上記平面は、路面モデルRMの面に対して垂直な平面であれば、タイヤモデルTMの移動方向に対してどのような方向に傾斜角θをもって傾斜して設けられてもよい。図3は、タイヤモデルTMを路面モデルRMの上方から見た図である。
シミュレーション結果処理部28は、飛散する粒子モデルルPMを平面に投影する処理を、設定された時間間隔毎に行い、粒子モデルルPMの投影画像を時系列に取得する。
Simulation result processing section 28, during the simulation, scattered particles Moderuru P M, performs a process of projecting the plane perpendicular to the plane of the road surface model R M.
Specifically, the simulation result processing section 28, the tire model T M is a predetermined time range to the passage end from passing the start of the water film model S M every predetermined time interval, the particle model scattering P M Are projected onto the plane. The plane used for the projection is a plane perpendicular to the plane of the road surface model R M passing through the center of the tire model T M. Since the rolling tire model T M moves with respect to the road surface model R M , the plane is provided every set time, and the scattered particle model P M is projected onto the plane. The plane is normal to the direction of this plane is provided such that the moving direction of the tire model T M. In addition, as shown in FIG. 3, if the plane is a plane perpendicular to the plane of the road surface model R M , the inclination angle θ in any direction with respect to the moving direction of the tire model T M May be provided with an inclination. FIG. 3 is a view of the tire model T M as viewed from above the road surface model R M.
Simulation result processing section 28, a process of projecting the particles Moderuru P M that scattered plane, performed every time the set interval, acquires the projected image of the particle Moderuru P M in a time series.
評価部30は、上記平面に投影された粒子モデルルPMの投影像を用いて、タイヤの水跳ね特性を評価する。
図4は、(a)〜(c)は、粒子モデルPMの投影像を用いた水跳ねの評価を行う流れを示す図である。図中、Y方向は、図1に示すように、タイヤモデルTMの移動方向と直交するタイヤ幅方向に該当し、タイヤモデルTMの回転軸に平行な方向である。Z方向は、路面モデルRMの面に垂直な方向である。
まず、図4(a)に示されるように、評価部30は、設定された時間における粒子モデルルPMの投影像毎に、Y方向の予め定められた各位置における粒子モデルルPMのZ方向の最高位置をプロットする。
次に、評価部30は、図4(b)に示されるように、プロットされた点を用いて曲線によるカーブフィットを行う。すなわち、評価部30は、粒子モデルPMの投影像の包絡線の形状を作成する。包絡線は、例えば、多項式によって表される曲線であることが好ましく、多項式は、a・(x−x0) b+c(a,b,c,x0は定数であり、bは自然数)であることが好ましい。多項式a・(x−x0) b+cは、粒子モデルPMの投影像のプロットされた各点を精度良くカーブフィットすることができる。
次に、評価部30は、図4(c)に示されるように、作成された包絡線において、路面モデルRMの面からの最高高さH、包絡線の路面モデルRMの面への交差位置L、および、包絡線のタイヤモデルTMからの出射角度Φを算出する。算出された最高高さH、交差位置Lおよび出射角度Φは、メモリ16に記憶される。
Evaluation unit 30 uses the projection image of the projected particles Moderuru P M in the plane, to evaluate the water splash property of the tire.
Figure 4, (a) ~ (c) is a diagram showing the flow of the evaluation of water splash with projected image of the particle model P M. In the figure, Y-direction, as shown in FIG. 1, corresponds to the tire width direction orthogonal to the moving direction of the tire model T M, is a direction parallel to the rotation axis of the tire model T M. Z direction is a direction perpendicular to the plane of the road surface model R M.
First, as shown in FIG. 4 (a), the evaluation unit 30, for each projection image of the particle Moderuru P M at the time set, Z direction of the particle Moderuru P M at predetermined respective positions in the Y-direction Plot the highest position of.
Next, as shown in FIG. 4B, the evaluation unit 30 performs curve fitting using a curve using the plotted points. That is, the evaluation unit 30 creates an envelope shape of the projection image of the particle model P M. The envelope is preferably a curve represented by a polynomial, for example, and the polynomial is a · (x−x 0 ) b + c (a, b, c, x 0 are constants, and b is a natural number). Preferably there is. The polynomial a · (x−x 0 ) b + c can curve-fit each plotted point of the projection image of the particle model P M with high accuracy.
Next, the evaluation unit 30, as shown in FIG. 4 (c), the envelope is created, the maximum height H from the surface of the road surface model R M, to the surface of the road surface model R M of the envelope intersection L, and calculates the output angle Φ from the tire model T M envelope. The calculated maximum height H, intersection position L, and emission angle Φ are stored in the memory 16.
評価部30は、時系列の粒子モデルルPMの投影画像毎に、図4(a)〜(c)に示す処理を行い、最も大きな包絡線を用いて、水跳ね特性を評価する。
評価部30は、最も大きな包絡線の最高高さH、交差位置Lおよび出射角度Φを用いて、他のタイヤにおける対応する最高高さH、交差位置Lおよび出射角度Φと比較することにより、水跳ね特性に関するタイヤの順位付けの評価を行う。また、評価部30は、カーブフィットにより得られた包絡線の形状を描画するために、ディスプレイ、プリンタ等の出力装置22に出力する。
なお、評価部30は最高高さH、交差位置Lおよび出射角度Φの3つを算出して評価に用いるが、最高高さH、交差位置Lおよび出射角度Φの少なくとも1つを算出して評価に用いることもできる。
Evaluation unit 30, for each projection image of the particle Moderuru P M of the time series, performs the processing shown in FIG. 4 (a) ~ (c) , using the largest envelope, to evaluate the water splash property.
The evaluation unit 30 uses the maximum height H, the intersection position L, and the emission angle Φ of the largest envelope to compare with the corresponding maximum height H, intersection position L, and emission angle Φ in other tires, Evaluate the ranking of tires with respect to water splash characteristics. Further, the evaluation unit 30 outputs the envelope shape obtained by curve fitting to an output device 22 such as a display or a printer.
The evaluation unit 30 calculates and uses the maximum height H, the intersection position L, and the exit angle Φ for evaluation, but calculates at least one of the maximum height H, the intersection position L, and the exit angle Φ. It can also be used for evaluation.
図5は、評価部30においてカーブフィットにより得られる包絡線の例を示す図である。図5では、転動速度が40km/時と60km/時である2つの走行速度における、タイヤA,Bの包絡線を示している。タイヤA,Bは、タイヤ形状が異なっている。タイヤBはタイヤAに対して、いずれの走行速度においても最高高さHが低いことがわかる。このように、カーブフィットにより得られる包絡線を用いて、水跳ね特性を効率よく、定量的に評価することができる。 FIG. 5 is a diagram illustrating an example of an envelope obtained by curve fitting in the evaluation unit 30. FIG. 5 shows the envelopes of the tires A and B at two traveling speeds where the rolling speed is 40 km / hour and 60 km / hour. The tires A and B have different tire shapes. It can be seen that the maximum height H of the tire B is lower than that of the tire A at any travel speed. Thus, the water splash characteristic can be efficiently and quantitatively evaluated using the envelope obtained by curve fitting.
また、評価装置10は、粒子モデルPMを用いるので、水跳ね特性の評価は、複数の粒子モデルPMが路面モデルRMの面へ落下する位置における粒子モデルPMの密度分布を用いて行われてもよい。例えば、Y方向の各位置において路面モデルRM上に落下する粒子モデルPMの密度分布を求め、この密度分布により水跳ね特性を評価することもできる。例えば、密度が所定値を超えるY方向位置を評価指標とする。 The evaluation device 10, since use particles model P M, the evaluation of water splash characteristics, using the density distribution of the particle model P M at the position where the plurality of particles model P M falls into the plane of the road model R M It may be done. For example, determine the density distribution of the particle model P M falling on the road surface model R M at each position in the Y direction, it is also possible to evaluate the water splash property by the density distribution. For example, the Y direction position where the density exceeds a predetermined value is used as the evaluation index.
なお、評価装置10は、コンピュータを実行させることで上記機能を発揮する装置である。このようなコンピュータのメモリ16には、以下のようなプログラムが記憶され、このプログラムが随時呼び出されて評価装置10として機能する。
すなわち、メモリ16には、路面に設けられる水膜等の介在物上をタイヤが通過するときの、タイヤの介在物飛び跳ね特性をコンピュータに評価させるプログラムが記憶される。このプログラムは、
タイヤを再現したタイヤモデルTMと、路面を再現した路面モデルRMと、路面上に設けられる水膜を再現した複数の粒子モデルPMをお互いに分離可能に含む水膜モデルSMとをコンピュータが作成する手順と、
水膜モデルSMが設けられた路面モデルRM上でタイヤモデルTMを転動させることにより、コンピュータが粒子モデルPMを飛散させる手順と、
飛散した粒子モデルPMを、コンピュータが、路面モデルRMの面に対して垂直な平面に投影する手順と、
平面に投影された粒子モデルPMの投影像を用いて、コンピュータが、タイヤの水跳ね特性を評価する手順と、を有する。
The evaluation device 10 is a device that exhibits the above functions by causing a computer to execute. The following program is stored in the memory 16 of such a computer, and this program is called as needed to function as the evaluation device 10.
That is, the memory 16 stores a program that causes the computer to evaluate the inclusion jumping characteristic of the tire when the tire passes over the inclusion such as a water film provided on the road surface. This program
A tire model T M that reproduces a tire, a road surface model R M that reproduces a road surface, and a water film model S M that includes a plurality of particle models P M that reproduce a water film provided on the road surface in a separable manner. The steps the computer creates,
A procedure for causing the computer to scatter the particle model P M by rolling the tire model T M on the road surface model R M provided with the water film model S M ;
A procedure in which the computer projects the scattered particle model P M onto a plane perpendicular to the plane of the road surface model R M ;
The computer has a procedure for evaluating the water splash characteristics of the tire using the projection image of the particle model P M projected on the plane.
(タイヤの水跳ね特性の評価方法)
図6は、タイヤの水跳ね特性の評価方法のフローの一例を示すフローチャートである。
まず、オペレータからの入力操作系20を通して入力を受け、この入力に基づいて、各モデルの作成条件やタイヤの転動シミュレーションの走行条件等を含んだ各種条件が設定される。設定された条件はメモリ16に記憶される。
(Evaluation method for water splash characteristics of tires)
FIG. 6 is a flowchart illustrating an example of a flow of a method for evaluating the water splash characteristics of a tire.
First, an input is received from an operator through an input operation system 20, and various conditions including conditions for creating each model and running conditions for tire rolling simulation are set based on this input. The set condition is stored in the memory 16.
次に、モデル作成部24は、メモリ16から呼び出されたモデルの作成条件に応じて、タイヤモデルTM、路面モデルRM、および、粒子モデルPMが一定間隔で配列された水膜モデルSMを作成する(ステップS10)。
タイヤモデルTMは、例えば3次元有限要素モデルであり、各要素に材料定数が付与されている。
路面モデルRMは、剛体モデル、あるいは有限要素により作られた有限要素モデルである。この有限要素モデルには材料定数が付与されている。
水膜モデルSMは粒子モデルPMが図2(b)に示すように層状に形成される。水膜モデルSMは、少なくとも2層以上の粒子モデルPMを有するのが、精度の良い水跳ね特性の評価を得る点で好ましい。なお、粒子モデルPMの挙動を支配する後述する定式化された式中のパラメータの値を変更することにより、水の代わりに、雪、砂、泥、および砂利のいずれか1つを含む介在物を再現することができる。
Next, the model creation unit 24 sets the water film model S in which the tire model T M , the road surface model R M , and the particle model P M are arranged at regular intervals according to the model creation conditions called from the memory 16. M is created (step S10).
The tire model TM is, for example, a three-dimensional finite element model, and a material constant is given to each element.
The road surface model RM is a rigid body model or a finite element model made of finite elements. A material constant is given to this finite element model.
In the water film model S M, the particle model P M is formed in layers as shown in FIG. The water film model S M preferably has at least two particle models P M from the viewpoint of obtaining an accurate evaluation of water splash characteristics. In addition, by changing the value of a parameter in a formula that will be described later that governs the behavior of the particle model P M, the inclusion includes any one of snow, sand, mud, and gravel instead of water. The thing can be reproduced.
次に、シミュレーション演算部26は、タイヤモデルTMに対して走行条件を設定して、タイヤの転動シミュレーションを実行し、粒子モデルPMを飛散させる(ステップS20)。
走行条件は、空気圧、負荷荷重、転動速度、横力、キャンバ角、制駆動力を含む。これらの条件は、予めメモリ16に記憶されており、シミュレーション演算部26から呼び出されて設定される。
タイヤモデルTMの転動シミュレーションでは、タイヤモデルTMに対して空気圧充填処理、接地処理、および転動処理が実行される。
空気圧充填処理では、タイヤモデルTMのタイヤ空洞領域に面する内周面の各節点に、設定された圧力に相当する力が与えられる。
Next, the simulation calculation unit 26 sets traveling conditions for the tire model T M , executes a tire rolling simulation, and scatters the particle model P M (step S20).
The traveling conditions include air pressure, load load, rolling speed, lateral force, camber angle, braking / driving force. These conditions are stored in advance in the memory 16 and are called from the simulation calculation unit 26 and set.
The rolling simulation of the tire model T M, tire inflation process to the tire model T M, grounding, and the rolling process is performed.
Pneumatic filling process, each node of the inner peripheral surface facing the tire cavity region of the tire model T M, given the force corresponding to the set pressure.
接地処理では、タイヤモデルTMが路面モデルRMに、所定の刻み幅で近づけられて、路面モデルRMとの間での接触の有無が判定される。さらに、接触するときの路面モデルRMからタイヤモデルTMの各節点が受ける力が算出され、この算出した力の合計値が設定された負荷荷重になるまで、タイヤモデルTMが路面モデルRMに近づけられる。 In grounding, the tire model T M is the road surface model R M, is brought close at a certain interval, the presence or absence of contact between the road surface model R M is determined. Furthermore, each node receives the force of the tire model T M from the road surface model R M when the contact is calculated, until the sum of the calculated force is set applied load, the tire model T M is the road surface model R Can approach M.
転動処理では、タイヤモデルTMに転動速度0から設定された転動速度まで徐々に速度が付与されることで、設定された転動速度で転動するタイヤモデルTMを作成する。このとき、タイヤモデルTMの各節点に、所定の時間刻み幅と転動速度とによって定まる変位を強制的に与えて回転させることで、路面モデルRM上を転がる転動状態が再現される。この転動処理では、タイヤモデルTMの各節点に路面モデルRMとの間の摩擦係数が与えられ、タイヤモデルTMの路面モデルRMに対する粘着、滑りの状態が再現される。
なお、タイヤモデルTMが水膜モデルSM上を通過するときの、粒子モデルPMとタイヤモデルTMとの相互作用については、後述する。粒子モデルPMは、粒子モデルPMとタイヤモデルTMとの相互作用により、飛散する。
In the rolling process, by gradually speed up the rolling speed set by the rolling speed 0 to the tire model T M is applied, creating a tire model T M rolling at a rolling speed which is set. In this case, each node of the tire model T M, is rotated by applying forced displacement determined by the predetermined time step size and the rolling speed, the rolling condition is reproduced rolling on the road surface model R M . In rolling process, the friction coefficient between the road model R M is applied to each node of the tire model T M, adhesion to the road surface model R M of the tire model T M, the state of slip is reproduced.
Incidentally, when the tire model T M passes over the water film model S M, the interaction between the particle model P M and the tire model T M will be described later. The particle model P M is scattered by the interaction between the particle model P M and the tire model T M.
次に、シミュレーション結果処理部28は、飛散した粒子モデルPMを、設定された時間間隔毎に、路面モデルRMの面に対して垂直な平面上に投影する(ステップS30)。上記平面は、タイヤモデルTMの中心を通り、タイヤ幅方向に平行な面であることが、水跳ね特性を精度良く評価する点で好ましい。 Next, the simulation result processing unit 28 projects the scattered particle model P M on a plane perpendicular to the plane of the road surface model R M for each set time interval (step S30). The plane through the center of the tire model T M, it is a plane parallel to the tire width direction is preferable in terms of accurately evaluating the water splash property.
次に、評価部30は、飛散する粒子モデルPMの上記平面上の投影像から、図4(a),(b)に示すように、粒子モデルPMの投影像の包絡線を形成する(ステップS40)。具体的には、評価部30は、図4(a)に示されるように、Y方向の各位置における粒子モデルPMの投影像の最高位置をプロットし、このプロットした点を用いて、図4(b)に示すように、多項式で表される曲線、例えば、a・(x−x0) b+c(a,b,c,x0は定数であり、bは自然数)によるカーブフィットを行う。こうして、評価部30は粒子モデルPMの投影像の包絡線を作成する。なお、飛散する粒子モデルPMの上記平面上の投影像は、設定された時間間隔毎に時系列で得られるので、評価部30は、粒子モデルPMの投影像の包絡線のうち最も大きな包絡線を水跳ね特性の評価用として選択する。 Next, the evaluation unit 30 from the projection image on the plane of the particle model P M scattered, as shown in FIG. 4 (a), (b) , to form the envelope of the projection image of the particle model P M (Step S40). Specifically, the evaluation unit 30, as shown in FIG. 4 (a), by plotting the highest position of the projected image of the particle model P M at each position in the Y direction, using the point that the plot, FIG. 4 (b), a curve fit by a curve represented by a polynomial, for example, a · (x−x 0 ) b + c (a, b, c, x 0 are constants and b is a natural number) is performed. Do. Thus, the evaluation unit 30 creates an envelope of the projection image of the particle model P M. In addition, since the projected image on the plane of the scattering particle model P M is obtained in time series for each set time interval, the evaluation unit 30 is the largest of the envelopes of the projected image of the particle model P M. An envelope is selected for evaluation of water splash characteristics.
次に、評価部30は、選択された包絡線のパラメータ、すなわち、最高高さH、交差位置Lおよび出射角度Φの少なくともいずれか1つを、水跳ね特性の評価パラメータとして抽出する(ステップS50)。
水跳ね特性の評価は、抽出された包絡線のパラメータの値を用いて、異なる仕様のタイヤモデル間の優劣あるいは順位付けにより行われる。あるいは、水跳ね特性の評価は、抽出された包絡線のパラメータの値が、予め定められた評価パラメータの許容範囲に含まれるか否かの判定により行われてもよい。
Next, the evaluation unit 30 extracts at least one of the selected envelope parameters, that is, the maximum height H, the intersection position L, and the emission angle Φ as an evaluation parameter for water splash characteristics (step S50). ).
The evaluation of water splash characteristics is performed by using superiority or inferiority or ranking between tire models having different specifications, using the extracted envelope parameter values. Alternatively, the evaluation of the water splash characteristic may be performed by determining whether or not the value of the extracted envelope parameter is included in the predetermined allowable range of the evaluation parameter.
なお、タイヤモデルTM、路面モデルRMおよび水膜モデルSMを作成する際、評価装置10は、タイヤを装着する車両のホイールハウスを少なくとも再現した車両モデルを作成し、粒子モデルPMを飛散させる際、タイヤモデルTMを車両モデルに装着してもよい。この場合、飛散した粒子モデルPMが、車両モデルに衝突する状態を調べることもできる。 When creating the tire model T M , the road surface model R M, and the water film model S M , the evaluation apparatus 10 creates a vehicle model that reproduces at least the wheel house of the vehicle on which the tire is mounted, and generates the particle model P M. when to scatter, it is also possible to attach a tire model T M on the vehicle model. In this case, scattered particles model P M is also possible to examine the state of collision to the vehicle model.
このようなタイヤモデルTMの水跳ね特性を精度良く、効率よく、定量的に評価を行うためには、以下のような粒子モデルPMを用いることが好ましい。
以下、粒子モデルPM、および粒子モデルPMとタイヤモデルTMの相互作用について説明する。
In order to evaluate the water splash characteristics of the tire model T M accurately, efficiently, and quantitatively, it is preferable to use the following particle model P M.
Hereinafter, the particle model P M and the interaction between the particle model P M and the tire model T M will be described.
(粒子モデル)
一般に、有限要素法によって作成されるモデルは、図7(a)に示すように、有限要素の応力や歪み等の物理量の分布を、有限要素を構成する節点の物理量を用いて表現する。一方、図7(b)に示されるように、空間に固定された空間格子内で流体要素(図7(b)中の斜線部分)を移動させ、空間格子で仕切られた領域における物理量を用いて水膜モデルを表現する方法もある。しかし、本実施形態の水膜モデルSMは、図7(c)に示すように、複数の粒子モデルPMで構成し、この粒子モデルPMを一定の間隔で等方状に配列した状態から、後述する移動規定条件の下に粒子モデルPMを移動させて、粒子モデルPMの持つ速度や密度や全エネルギー量を用いて表現する。そのため、従来の方法のように、空間格子や有限要素が設けられない。しかも、粒子モデルPMは、粒子モデル移動規定条件の下に自由に移動するので、大きな移動や飛散を伴う挙動を流体に固定した座標系で記述するLagrangianに基づいて表現しても、数値計算上の適切な解を得ることができる。
(Particle model)
In general, a model created by the finite element method represents the distribution of physical quantities such as stress and strain of a finite element using physical quantities of nodes constituting the finite element, as shown in FIG. On the other hand, as shown in FIG. 7 (b), the fluid element (the hatched portion in FIG. 7 (b)) is moved in the space lattice fixed in the space, and the physical quantity in the region partitioned by the space lattice is used. There is also a way to express a water film model. However, as shown in FIG. 7C, the water film model S M of the present embodiment is composed of a plurality of particle models P M , and the particle models P M are arranged in an isotropic manner at regular intervals. From this, the particle model P M is moved under the movement regulation condition described later, and is expressed using the speed, density, and total energy amount of the particle model P M. Therefore, unlike the conventional method, no spatial grid or finite element is provided. Moreover, since the particle model P M moves freely under the particle model movement regulation conditions, even if the behavior accompanied by large movement and scattering is expressed based on Lagrangian described in the coordinate system fixed to the fluid, numerical calculation You can get the appropriate solution above.
このような粒子モデルPMを用いて物理的挙動を演算し予測することは、一般に天文学の分野における圧縮性流れの数値解析法におけるSPH(Smoothed Particle Hydrodynamics)法として公知である(Monaghan,J.J.,Smoothed Particle Hydrodynamics ,Annu.Rev.Astron.Astrophys.vol.30,1992,pp543-574) 。本実施形態は、一定の間隔で等方状に配列した複数の粒子モデルPMからなる水膜モデルSMを用いて、タイヤモデルTMと粒子モデルPMとの相互作用を演算して、粒子モデルPMの飛散を含む移動を計算する。 The calculation and prediction of physical behavior using such a particle model P M is generally known as the SPH (Smoothed Particle Hydrodynamics) method in the numerical analysis method of compressible flow in the field of astronomy (Monaghan, JJ, Smoothed Particle Hydrodynamics, Annu. Rev. Astron. Astrophys. Vol. 30, 1992, pp 543-574). In the present embodiment, the interaction between the tire model T M and the particle model P M is calculated using a water film model S M composed of a plurality of particle models P M arranged in an isotropic manner at regular intervals. The movement including the scattering of the particle model P M is calculated.
粒子モデルPMを用いて表す流体の挙動は、公知のように、下記式(1)〜(3)に示す流体の連続の式、運動方程式およびエネルギー方程式の各式を支配方程式として用いて表される。
ここで、ρは流体の密度、Uは流体の速度、σは流体の応力テンソルであり、流体の圧力pとσ=−pI(Iは単位テンソル)の関係にある量、eは単位体積当たりの全エネルギー量である。
Table behavior of fluids expressed using a particle model P M, using as known, continuity equation of fluid represented by the following formula (1) to (3), as the governing equation of the equations of the motion equation and the energy equation Is done.
Here, ρ is the density of the fluid, U is the velocity of the fluid, σ is the stress tensor of the fluid, the quantity is in the relationship between the pressure p of the fluid and σ = −pI (I is the unit tensor), e is per unit volume The total amount of energy.
一方、流体の持つ速度(流速)等の注目物理量の分布f(x)(xは空間の3次元位置座標を表す)は、下記式(4)に示すように、所定の位置座標x’における注目物理量f(x’)とカーネル関数W(|x−x’|,h)(hはパラメータ)とを用いて求められる近似物理量の分布<f(x)>によって近似することができる。
そこで、このときの粒子モデルPMの位置座標x’における注目物理量f(x’)を、粒子モデルPj の持つ注目物理量fj (jは1以上N以下の自然数で、Nは粒子モデルPの総数)とし、粒子モデルPMの密度ρj と、粒子モデルPMの質量mj とを用いて下記式(5)のように離散化して表すことができる。例えば、流体の速度U(x)について下記式(6)のように表すことができる。
同様に流体の密度ρ(x)および全エネルギー量e(x)も上記式(5)の注目物理量f(x)の替わりに用いて表すことができる。すなわち、上記式(1)〜(3)における流体の速度Uや密度ρや全エネルギー量eを、粒子モデルPMの有する速度Uj や密度ρj を用いた近似速度<U(x)>、近似密度<ρ(x)>、近似全エネルギー量<e(x)>を用いて表すことができる。そこで、上記式(1)〜(3)のU、ρ、eに近似速度<U(x)>、近似密度<ρ(x)>、近似全エネルギー量<e(x)>を代入して下記式(7)〜(9)を求めることができる。iは、複数の粒子モデルPMのうち注目する粒子モデルの番号を指す。以降注目する粒子モデルをPi とする。
ここで、▽Wijは、下記式(10)を表す。
なお、▽は空間微分を表す。また、πijは人工粘性を表し、後述する変形計算において数値上の振動解の発生を抑制し、また異なる粒子モデルPMが衝突する際に互いに相手をすり抜けないように、運動量や全エネルギー量を交換するもので、予め設定されるものである。
Similarly, the fluid density ρ (x) and the total energy amount e (x) can also be expressed using the physical quantity f (x) of interest in the above equation (5). That is, the fluid velocity U, density ρ, and total energy amount e in the above formulas (1) to (3) are approximated using the velocity U j and density ρ j of the particle model P M <U (x)>. , Approximate density <ρ (x)>, approximate total energy <e (x)>. Therefore, substituting approximate velocity <U (x)>, approximate density <ρ (x)>, and approximate total energy <e (x)> into U, ρ, and e in the above equations (1) to (3). The following formulas (7) to (9) can be obtained. i indicates the number of the particle model of interest among the plurality of particle models P M. Hereinafter, the particle model to be noted is P i .
Here, ▽ W ij represents the following formula (10).
In addition, ▽ represents spatial differentiation. In addition, π ij represents artificial viscosity, suppresses the generation of numerical vibration solutions in the deformation calculation described later, and does not slip through each other when different particle models P M collide with each other. Are set in advance.
このように、速度Uや密度ρや全エネルギー量e等の注目物理量の分布f(x)を粒子モデルPの持つ注目物理量fj を用いて近似することによって定められる式(7)〜(9)が作成される。なお、本実施形態では、粒子モデルPMが式(7)〜(9)をすべて備える必要はなく、式(9)を、粒子モデルPMの移動を規定する粒子モデル移動規定条件として備えた水膜モデルSMのであっても良い。以降では、式(8)を粒子モデル移動規定条件として備えた水膜モデルSMに基づいて説明する。 Thus, the equations (7) to (9) defined by approximating the distribution f (x) of the physical quantity of interest such as the velocity U, the density ρ, and the total energy amount e using the physical quantity f j of the particle model P. ) Is created. In the present embodiment, the particle model P M does not have to include all of the equations (7) to (9), and the equation (9) is included as a particle model movement defining condition that defines the movement of the particle model P M. it is in the water film model S M may be. In the following, description will be made based on a water film model S M provided with Equation (8) as a particle model movement regulation condition.
ここで、カーネル関数W(|x−x’|,h)は、特に限定されないが、|x−x’|/h=νとして、例えば、下記式(11)に示すように、3次のB−スプライン関数を用いるとよい。
式(11)に示す3次のB−スプライン関数をカーネル関数W(|x−x’|,h)として用いると、νが2以上の場合、カーネル関数W(|x−x’|,h)は0となり、νが2未満の場合正の値となるため、式(7)〜(9)に示すように、注目する粒子モデルPi の挙動は、νが2未満の範囲内に位置する粒子モデルPj に依存することになる。すなわち、図8(a)に示すように、粒子モデルPi を中心とする半径2hの球の範囲(この範囲を粒子モデルPi に対応した近傍領域Ri という)内に含まれる粒子モデルPj によって、粒子モデルPiの移動が規定される。この半径2hはスムージングレングスともいう。また、水膜モデルSMでは、図8(b)に示すように、水膜モデルSMの全領域を覆うように、粒子モデルPi のそれぞれの周りに近傍領域Ri が粒子モデルPi のそれぞれに対応して定められる。しかも、この近傍領域Ri 内の速度U等の注目物理量の分布をこの近傍領域Ri に含まれる粒子モデルPj の持つ速度Uj 等の注目物理量を用いて近似することによって定められる、粒子モデルPi の移動を規定する速度Ui 等に関する方程式(式(7)〜(9))を備える。このような複数の粒子モデルPMから構成される水膜モデルSMによって水膜が再現される。 When the cubic B-spline function shown in Expression (11) is used as the kernel function W (| x−x ′ |, h), when ν is 2 or more, the kernel function W (| x−x ′ |, h ) Becomes 0 and becomes a positive value when ν is less than 2, as shown in equations (7) to (9), the behavior of the particle model P i of interest is located within the range where ν is less than 2. It depends on the particle model P j to be used. That is, as shown in FIG. 8 (a), the particle model P included in a range of a sphere having a radius 2h centered on the particle model P i (this range is referred to as a neighboring region R i corresponding to the particle model P i ). j defines the movement of the particle model P i . This radius 2h is also called a smoothing length. Further, the water film model S M, as shown in FIG. 8 (b), so as to cover the entire area of the water film model S M, particle model P region near around each of the i R i particle model P i It is determined corresponding to each of. Moreover, it determined by approximating with attention physical quantity of speed U j or the like having a distribution of attention physical quantity of speed U and the like of the neighboring region R i particle model P j included in the neighboring region R i, particles Equations (Expressions (7) to (9)) relating to the speed U i and the like that define the movement of the model P i are provided. A water film is reproduced by such a water film model S M composed of a plurality of particle models P M.
(タイヤモデルと粒子モデルとの相互作用の計算)
次に、タイヤモデルTMと複数の粒子モデルPMとの相互作用の演算を説明する。
具体的には、タイヤモデルTMの変形計算と複数の粒子モデルPMの移動の計算が、互いのモデルが及ぼす相互作用を考慮して行われる。これらの計算、すなわち、タイヤモデルTMの変形計算、粒子モデルPMの変形計算、およびタイヤモデルTMと複数の粒子モデルPとの相互作用の演算は、所定の時間間隔毎の時間ステップ毎に逐次行われる。具体的には、タイヤモデルTMが水膜モデルSM上の通過を開始する直前のタイヤモデルTMの各要素の応力や速度や加速度等の物理量が取り出され、この物理量を初期条件として、粒子モデルPMに速度が与えられる。
(Calculation of interaction between tire model and particle model)
Next, the calculation of the interaction between the tire model T M and the plurality of particle models P M will be described.
Specifically, the deformation calculation of the tire model T M and the movement calculation of the plurality of particle models P M are performed in consideration of the interaction between the models. These calculations, that is, the deformation calculation of the tire model T M , the deformation calculation of the particle model P M , and the calculation of the interaction between the tire model T M and the plurality of particle models P are performed at each time step at predetermined time intervals. It is performed sequentially. Specifically, the physical quantity such as stress or speed and acceleration of each element of the tire model T M immediately before the tire model T M starts passing over the water film model S M is taken out, the physical quantity as an initial condition, A speed is given to the particle model P M.
図8は、転動シミュレーションの要部の処理フローを示すフローチャートである。
まず、それぞれの粒子モデルPMの位置が既知の状態にある、時刻T(n) の時間ステップにおいて、粒子モデルPM各々のスムージングレングスの設定、すなわち、上記半径2hの球の範囲(近傍領域Ri ) を粒子モデルPi のそれぞれについて設定し(ステップS100)、近傍領域Ri に含まれる近接する粒子モデルPj を探索し求める(ステップS102)。この後、近接する粒子モデルPj の数に基づいて粒子モデルPi における密度ρi を算出するとともに、近傍領域Ri における歪みや歪み速度を粒子モデルPj の位置および速度を用いて算出する(ステップS104)。この後、粒子モデルPi における密度ρi から水を支配する状態方程式と等温変化の条件とを用いて粒子モデルPi にかかる圧力を求める。具体的には、密度に比例する形で粒子モデルPi にかかる圧力を算出する。さらに、粒子モデルPi のそれぞれの有する運動エネルギーおよび歪みエネルギーを算出する(ステップS106)。算出される粒子モデルPi における圧力をpとすると式(8)中のσi とσi =−pI(Iは単位テンソル)の関係があり、この関係を用いて、後述する粒子モデルPMの移動の計算において用いられる。さらに、粒子モデルPMの移動の計算において粘性を考慮する場合、近傍領域Ri における歪み速度から粒子モデルPi における粘性応力が求められ、式(8)中のσi に加えられる。これにより、雪、砂、泥、および砂利の介在物の挙動を再現することができる。なお、近傍領域Ri が粒子モデルPj を全く含まない場合、式(8)中のmjは0となり右辺は0となる。
FIG. 8 is a flowchart showing a processing flow of a main part of the rolling simulation.
First, in the time step of time T (n) when the position of each particle model P M is in a known state, the setting of the smoothing length of each particle model P M , that is, the range of the sphere having the radius 2h (neighboring region) R i ) is set for each of the particle models P i (step S100), and the adjacent particle models P j included in the neighboring region R i are searched for (step S102). Thereafter, to calculate the density [rho i in the particle model P i on the basis of the number of particles model P j adjacent, is calculated by using the position and velocity of the strain and strain rate of particles model P j in the neighboring region R i (Step S104). Thereafter, determining the pressure applied to the particle model P i by using the condition of the state equation and the isothermal change governing the water from the density [rho i in the particle model P i. Specifically, the pressure applied to the particle model P i is calculated in proportion to the density. Further, the kinetic energy and strain energy of each of the particle models P i are calculated (step S106). If the pressure in the calculated particle model P i is p, there is a relationship between σ i in equation (8) and σ i = −pI (I is a unit tensor), and a particle model P M to be described later using this relationship. Used in the calculation of movement. Furthermore, when considering the viscosity in the calculation of the movement of the particles model P M, viscosity stress is obtained from the strain rate in the neighboring region R i in the particle model P i, is added to the sigma i in equation (8). Thereby, the behavior of inclusions of snow, sand, mud, and gravel can be reproduced. Note that when the neighboring region R i does not include the particle model P j at all, m j in the equation (8) is 0 and the right side is 0.
一方、粒子モデルPMと同時刻T(n) におけるタイヤモデルTMの各有限要素における歪みが算出され(ステップS108)、さらに、タイヤモデルTMの各有限要素における応力、運動エネルギーや歪みエネルギーが算出される(ステップS110)。次に、時刻T(n) が予め定められた経過時刻を経過しているか否かが判断され(ステップS112)、予め定められた経過時刻を経過していないと判断した場合、時刻T(n) を時刻T(n+1) とし、タイヤモデルTMと粒子モデルPMの相互作用の演算(ステップS114)と、粒子モデルPMの移動の計算(水膜モデルSMの変形計算)(ステップS116)と、タイヤモデルTMの変形計算(ステップS118)とが行われる。ステップS112において時刻T(n) が予め定められた経過時刻を経過している場合、相互作用の計算は終了する。このようにして、予め定められた経過時刻までステップS100〜110、ステップS114〜118が繰り返される。 On the other hand, the strain in each finite element of the tire model T M at the same time T (n) as the particle model P M is calculated (step S108), and further the stress, kinetic energy and strain energy in each finite element of the tire model T M. Is calculated (step S110). Next, it is determined whether or not the time T (n) has passed a predetermined elapsed time (step S112). If it is determined that the predetermined elapsed time has not elapsed, the time T (n) ) At time T (n + 1) , the calculation of the interaction between the tire model T M and the particle model P M (step S114) and the movement of the particle model P M (deformation calculation of the water film model S M ) ( and step S116), deformation calculation of the tire model T M (step S118) is performed. If the time T (n) has passed a predetermined elapsed time in step S112, the calculation of the interaction ends. In this way, steps S100 to 110 and steps S114 to 118 are repeated until a predetermined elapsed time.
図10(a)〜(c)は、粒子モデルP2 における近傍領域Ri に含まれる粒子モデルPMが粒子モデルP1 、P3 であり、粒子モデルP2 がタイヤモデルTMと接触して、速度U2 で移動する例を示している。この例に沿って説明すると、時刻T(1) の時間ステップにおいて速度uを有して移動するタイヤモデルTMが粒子モデルP2 と接触する直前の状態(図10(a))から、時刻T(2) の時間ステップにおいて粒子モデルP2 が一旦タイヤモデルTMの内部に位置する状態(図10(b))に変化し、さらに、時刻T(3) の時間ステップにおいて粒子モデルP2 がタイヤモデルTMの外部に位置して移動する状態(図10(c))に変化するように、相互作用の演算(ステップS114)と、粒子モデルPの移動計算(ステップS116)と、タイヤモデルの変形計算(ステップS118)と、が行われる。 10A to 10C, the particle model P M included in the neighboring region R i in the particle model P 2 is the particle models P 1 and P 3 , and the particle model P 2 is in contact with the tire model T M. Thus, an example of moving at a speed U 2 is shown. Explaining along this example, from the state (FIG. 10A) immediately before the tire model T M that moves at a speed u at the time step of time T (1) contacts the particle model P 2 , the time At a time step of T (2) , the particle model P 2 is temporarily changed to a state (FIG. 10B) positioned inside the tire model T M , and further, at a time step of time T (3) , the particle model P 2 is changed. so it changed to a state (FIG. 10 (c)) which moves in a position outside the tire model T M, the calculation of the interaction (step S114), the mobile computing particle model P (step S116), the tire Model deformation calculation (step S118) is performed.
まず、時刻T(2) の時間ステップにおける粒子モデルPMに対してステップS100〜ステップS106の処理が行われるとともに、タイヤモデルTMに対してステップS108、ステップS110の処理が行われる。時刻T(2) が所定の経過時刻を経過しない時刻であるとステップS112において判断された場合、時刻T(2) の時間ステップの状態から時刻T(3) の時間ステップの状態へ移行する際のタイヤモデルTMと粒子モデルPMの相互作用の演算が行われる(ステップS114)。相互作用の演算では、まず、粒子モデルPMがタイヤモデルTMの内部に進入しているか否かが、粒子モデルPMの各々について判別される。粒子モデルP2 がタイヤモデルTM内部に進入したと判別されると、この状態の粒子モデルP2 が時刻T(3) の時間ステップにおいてタイヤモデルTMの外部に位置するように、図10(b)のようにタイヤモデルTMの内部に進入した粒子モデルP2 の進入量に応じた相互作用エネルギー項を算出し、この項が粒子モデルP2 に関する式(8)に付加される。 First, steps S100 to S106 are performed on the particle model P M at the time step of time T (2) , and steps S108 and S110 are performed on the tire model T M. If the time T (2) is determined in step S112 When it is time has not yet passed a predetermined time elapsed, when shifting from the state of the time steps in the time T (2) to the time step of the state at time T (3) The interaction between the tire model T M and the particle model P M is calculated (step S114). In the calculation of the interaction, first, whether the particle model P M has entered the inside of the tire model T M is determined for each particle model P M. If the particle model P 2 is determined to have entered the inside of the tire model T M, so as to be positioned outside of the tire model T M at time step of the particle model P 2 in this state the time T (3), 10 As shown in (b), an interaction energy term corresponding to the amount of the particle model P 2 entering the tire model T M is calculated, and this term is added to the equation (8) relating to the particle model P 2 .
この後、相互作用エネルギー項が付加された粒子モデルP2 に関する式(8)および相互作用エネルギー項が付加されていない粒子モデルP1 、P3 に関する式(8)を後述する差分スキームを用いて解き、粒子モデルP1 〜P3 の加速度を求め、速度および位置を算出する。これによって、時刻T (2)の時間ステップにおける粒子モデルPの移動の計算が行われ(ステップS116)、時刻T(3) の時間ステップにおける粒子モデルPMの移動の状態が算出される。一方、タイヤモデルTMには、粒子モデルP2 に付与される相互作用エネルギー項の反作用として、相互作用エネルギー項が付与されるため、この相互作用エネルギー項の付与された、時刻T(2) の時間ステップにおけるタイヤモデルTMの変形計算が行われ(ステップS118)、時刻T(3) の時間ステップにおけるタイヤモデルTMの変形状態が算出される。相互作用エネルギー項は粒子モデルP2がタイヤモデルTMに与える衝突圧力である。こうして算出された時刻T(3) の時間ステップにおける粒子モデルPMに対してステップS100〜S106の処理が行われるとともに、タイヤモデルTMに対してステップS108、ステップS110の処理が行われた後、ステップ112の判断が行われる。このようにして、所定の経過時刻を経過するまで、時間ステップ毎にタイヤモデルTMと粒子モデルPMの相互作用の演算と粒子モデルPMの移動の計算とタイヤモデルTMの変形計算が繰り返され、ステップS112で判断結果が肯定されると相互作用の計算は終了する。
なお、本実施形態では、粒子モデルPMの飛散の挙動を計算するので、相互作用の計算が終了した後も、粒子モデルPMの移動の計算は続行される。
Thereafter, Expression (8) relating to the particle model P 2 to which the interaction energy term is added and Expression (8) relating to the particle models P 1 and P 3 to which the interaction energy term is not added are used by a differential scheme described later. Solve, find the acceleration of the particle models P 1 to P 3 , and calculate the velocity and position. Thereby, the movement of the particle model P at the time step at time T (2) is calculated (step S116), and the movement state of the particle model P M at the time step at time T (3) is calculated. On the other hand, since an interaction energy term is given to the tire model T M as a reaction of the interaction energy term given to the particle model P 2 , the time T (2) at which this interaction energy term was given. The deformation calculation of the tire model T M at the time step is performed (step S118), and the deformation state of the tire model T M at the time step of time T (3) is calculated. The interaction energy term is the collision pressure that the particle model P 2 applies to the tire model T M. After the processing in steps S100 to S106 is performed on the particle model P M at the time step of time T (3) calculated in this way, and the processing in steps S108 and S110 is performed on the tire model T M. In step 112, a determination is made. In this manner, the calculation of the interaction between the tire model T M and the particle model P M , the calculation of the movement of the particle model P M , and the deformation calculation of the tire model T M are performed for each time step until a predetermined elapsed time elapses. Repeatedly, when the determination result in step S112 is affirmative, the calculation of the interaction ends.
In the present embodiment, since calculating the behavior of scattering of particles model P M, even after the calculation of the interaction is completed, the calculation of the movement of the particles model P M is continued.
なお、式(8)における時間微分の方程式は、中心差分による時間積分スキームを用いた解法を用いて解かれ、速度Ui が計算される。すなわち、下記式(12)のような差分スキームをとり、Ui n+1/2 (第n時間ステップと第(n+1)時間ステップとの間の中間の速度)から、下記式(13)を用いて第(n+1)時間ステップにおける粒子モデルPi の変位量di n+1 が求められる。さらに、変位量di n+1 が下記式(14)を用いて初期状態における位置座標xi 0 に加算されて、第(n+1)時間ステップにおける粒子モデルPi の位置座標xi n+1 が求められる。Δtn は、第n時間ステップにおける時間ステップの時間間隔である。時間ステップが一定の時間間隔毎のステップの場合、Δtn はnに関わらず、すべて一定である。こうして、各粒子モデルPMの位置座標が求められる。 Note that the time differential equation in equation (8) is solved using a solution using a time integration scheme based on the center difference, and the velocity U i is calculated. That is, a differential scheme such as the following equation (12) is taken, and the following equation (13) is obtained from U i n + 1/2 (an intermediate speed between the n-th time step and the (n + 1) -th time step). Using this, the displacement amount d i n + 1 of the particle model P i at the (n + 1) th time step is obtained. Further, the displacement d i n + 1 is added to the position coordinates x i 0 in the initial state using the following equation (14), and the position coordinates x i n + 1 of the particle model P i in the (n + 1) th time step. Is required. Δt n is a time interval of time steps in the nth time step. When the time steps are steps at regular time intervals, Δt n is all constant regardless of n. Thus, the position coordinates of each particle model P M is obtained.
なお、次の時間ステップにおいてタイヤモデルTM内に進入する粒子モデルPi が有るか否かを判別し、進入した粒子モデルPi に対して相互作用エネルギー項を式(8)に与えて、次の時間ステップで粒子モデルPi のタイヤモデルTM外に移動させる上述の方法の他に、粒子モデルPi のタイヤモデルTM内への進入が発生しない様に、粒子モデルPi の移動を計算する際、タイヤモデルTMの表面形状をラグランジェの未定定数法における拘束条件として付加し、この未定定数を求める方法を利用してもよい。この場合、数値として求まる未定定数が粒子モデルPi の進入を阻止する粒子モデルPMのタイヤモデルTMへの衝突圧力となる。 Incidentally, giving to determine whether particles model P i that enters the tire model T M at the next time step is present, the interaction energy terms against ingress particles model P i in equation (8), in addition to the above-described method of moving the tire model T M outside the particle model P i in the next time step, as entry into the tire model T M of the particle model P i does not occur, the movement of the particle model P i when calculating, adding the surface shape of the tire model T M as a constraint condition in undetermined constant method of Lagrange, it may be used a method of determining the undetermined constants. In this case, the undetermined constant obtained as a numerical value is the collision pressure of the particle model P M that prevents the particle model P i from entering the tire model T M.
以上のようにして、タイヤモデルTMの変形計算および粒子モデルPMの移動の計算が行われる。なお、本実施形態では、タイヤモデルTMに並進速度と並進速度に応じた回転速度を与えて、タイヤモデルTMが路面モデルRMに対して滑りがない(滑り率0)転動状態を演算するが、タイヤモデルTMの滑り率が0でない状態を演算してもよい。 As described above, the deformation calculation of the tire model T M and the movement of the particle model P M are calculated. In the present embodiment, it is given a rotational speed corresponding to the translational velocity and the translational speed of the tire model T M, the tire model T M will not slip relative to the road surface model R M (the slip ratio 0) tumbling condition Although the calculation is performed, a state where the slip rate of the tire model T M is not 0 may be calculated.
粒子法は、上記SPH法の他に、RKPM法(Reproducing Kernel Particle Method) や有限粒子法(Finite Particle Method)およびEFGM(Element Free Galerkin Method)法を含む。
あるいは、粒子法による粒子モデルを用いる代わりに、粒状体や水等の流体を複数の粒子モデルPMでモデル化し、複数の粒子モデルPM間を単純なばねやダッシュポット等で結合しあるいは粒子モデルPM間同士の接触を摩擦力を用いて表す方法(DEM:Discrete Element Method )を粒子モデルPMに適用することもできる。
In addition to the SPH method, the particle method includes an RKPM method (Reproducing Kernel Particle Method), a finite particle method (Finite Particle Method), and an EFGM (Element Free Galerkin Method) method.
Alternatively, instead of using the particle model by particle method, a fluid such as granules or water modeled by a plurality of particles model P M, and coupling between a plurality of particles model P M simple spring and dashpot or the like or particles A method (DEM: Discrete Element Method) that expresses contact between the models P M using frictional force can also be applied to the particle model P M.
このように、本実施形態の評価装置10および評価方法では、タイヤモデルTMが水を再現した粒子モデルPMを飛散させる際、路面モデルRMの面に対して垂直な平面に粒子モデルPMを投影し、投影された粒子モデルPMの投影像を用いるので、タイヤの水跳ね特性を効率よく、かつ定量的に評価することができる。 Thus, in the evaluation device 10 and the evaluation method of this embodiment, a tire model T M when a to scatter the particles model P M reproduces the water, road model R M particle model P in a plane perpendicular to the plane of the Since M is projected and the projected image of the projected particle model P M is used, the water splash characteristics of the tire can be evaluated efficiently and quantitatively.
以上、本発明のタイヤの介在物飛び跳ね特性の評価方法、タイヤの介在物飛び跳ね特性の評価装置およびプログラムについて詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良や変更をしてもよいのはもちろんである。 As described above, the method for evaluating the inclusion jumping characteristic of the tire, the evaluation device for the inclusion jumping characteristic of the tire, and the program according to the present invention have been described in detail. Of course, various improvements and changes may be made within the range not to be performed.
10 評価装置
12 CPU
14 バス
16 メモリ
18 入出力インターフェース部
20 入力操作系
22 出力装置
24 モデル作成部
26 シミュレーション演算部
28 シミュレーション結果処理部
30 評価部
36 処理モジュール群
10 Evaluation device 12 CPU
14 Bus 16 Memory 18 Input / output interface unit 20 Input operation system 22 Output device 24 Model creation unit 26 Simulation operation unit 28 Simulation result processing unit 30 Evaluation unit 36 Processing module group
Claims (9)
タイヤを再現したタイヤモデルと、路面を再現した路面モデルと、路面上に設けられる介在物を再現した、複数の微小モデルがお互いに分離可能に含まれる介在物モデルとを作成するステップと、
前記介在物モデルが設けられた前記路面モデル上で前記タイヤモデルを転動させることにより、前記微小モデルを飛散させるステップと、
前記飛散した前記微小モデルを、設定された時間間隔毎に、前記路面モデルの面に対して垂直な平面に投影するステップと、
前記平面に投影された前記微小モデルの時系列の投影像を用いて、前記投影像毎に、飛散する前記微小モデルの前記投影像の包絡線を求め、前記投影像毎に求めた前記包絡線の中から、前記包絡線の形状が最も大きい包絡線を用いてタイヤの介在物の飛び跳ね特性を評価するステップと、を有することを特徴とするタイヤの介在物飛び跳ね特性の評価方法。 A method for evaluating the jumping characteristics of a tire inclusion when the tire passes over an inclusion provided on a road surface,
Creating a tire model that reproduces a tire, a road surface model that reproduces a road surface, and an inclusion model that reproduces inclusions provided on the road surface and includes a plurality of minute models that are separable from each other;
Scattering the micro model by rolling the tire model on the road surface model provided with the inclusion model;
Projecting the scattered micro model on a plane perpendicular to the plane of the road surface model at set time intervals ;
Using the time-series projection image of the micro model projected on the plane, an envelope of the projection image of the micro model scattered is obtained for each projection image, and the envelope obtained for each projection image A method for evaluating the jumping characteristic of a tire inclusion using the envelope having the largest shape of the envelope , and the step of evaluating the jumping characteristic of the tire inclusion.
前記微小モデルを飛散させる際、前記タイヤモデルを前記車両モデルに装着して、前記微小モデルを飛散させる、請求項1〜5のいずれか1項に記載の評価方法。 Furthermore, when creating the tire model, the road surface model and the inclusion model, create a vehicle model that reproduces at least the wheel house of the vehicle on which the tire is mounted,
Time of scattering the said micro model, by mounting the tire model on the vehicle model, is scattered the micro model, the evaluation method according to any one of claims 1-5.
タイヤを再現したタイヤモデルと、路面を再現した路面モデルと、路面上に設けられる介在物を再現した、複数の微小モデルがお互いに分離可能に含まれる介在物モデルとを作成するモデル作成部と、
前記介在物モデルが設けられた前記路面モデル上で前記タイヤモデルを転動させることにより、前記微小モデルを飛散させるシミュレーション演算部と、
前記飛散した前記微小モデルを、設定された時間間隔毎に、前記路面モデルの面に対して垂直な平面に投影するシミュレーション結果処理部と、
前記平面に投影された前記微小モデルの時系列の投影像を用いて、前記投影像毎に、飛散する前記微小モデルの前記投影像の包絡線を求め、前記投影像毎に求めた前記包絡線の中から、前記包絡線の形状が最も大きい包絡線を用いて、タイヤの介在物の飛び跳ね特性を評価する評価部と、を有することを特徴とするタイヤの介在物飛び跳ね特性の評価装置。 An apparatus for evaluating the jumping characteristics of tire inclusions when the tire passes over inclusions provided on a road surface,
A model creation unit that creates a tire model that reproduces a tire, a road surface model that reproduces a road surface, and an inclusion model that reproduces inclusions provided on the road surface and includes a plurality of minute models that are separable from each other; ,
A simulation calculation unit that scatters the minute model by rolling the tire model on the road surface model provided with the inclusion model;
A simulation result processing unit that projects the scattered micro model on a plane perpendicular to the plane of the road surface model for each set time interval ;
Using the time-series projection image of the micro model projected on the plane, an envelope of the projection image of the micro model scattered is obtained for each projection image, and the envelope obtained for each projection image And an evaluation unit that evaluates the jumping characteristics of the inclusions of the tire using the envelope having the largest shape of the envelope .
タイヤを再現したタイヤモデルと、路面を再現した路面モデルと、路面上に設けられる介在物を再現した、複数の微小モデルがお互いに分離可能に含まれる介在物モデルとをコンピュータが作成する手順と、
前記介在物モデルが設けられた前記路面モデル上で前記タイヤモデルを転動させることにより、コンピュータが前記介在物モデルの前記微小モデルを飛散させる手順と、
前記飛散した前記微小モデルを、コンピュータが、設定された時間間隔毎に、前記路面モデルの面に対して垂直な平面に投影する手順と、
前記平面に投影された前記微小モデルの時系列の投影像を用いて、前記投影像毎に、飛散する前記微小モデルの前記投影像の包絡線を求め、前記投影像毎に求めた前記包絡線の中から、前記包絡線の形状が最も大きい包絡線を用いて、コンピュータが、タイヤの介在物の飛び跳ね特性を評価する手順と、を有することを特徴とするプログラム。 A program for causing a computer to evaluate the inclusion jumping characteristic of a tire when the tire passes over inclusions provided on a road surface,
A procedure for a computer to create a tire model that reproduces a tire, a road surface model that reproduces a road surface, and an inclusion model that reproduces inclusions provided on the road surface and includes a plurality of minute models that can be separated from each other. ,
A step of causing the computer to scatter the minute model of the inclusion model by rolling the tire model on the road surface model provided with the inclusion model;
The computer projects the scattered micro model on a plane perpendicular to the plane of the road surface model at set time intervals , and
Using the time-series projection image of the micro model projected on the plane, an envelope of the projection image of the micro model scattered is obtained for each projection image, and the envelope obtained for each projection image The computer has a procedure for evaluating the jump characteristics of the inclusions of the tire using the envelope having the largest shape of the envelope .
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010125706A JP5533277B2 (en) | 2010-06-01 | 2010-06-01 | Evaluation method of tire inclusion jumping characteristic, tire inclusion jumping characteristic evaluation apparatus and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010125706A JP5533277B2 (en) | 2010-06-01 | 2010-06-01 | Evaluation method of tire inclusion jumping characteristic, tire inclusion jumping characteristic evaluation apparatus and program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011252748A JP2011252748A (en) | 2011-12-15 |
JP5533277B2 true JP5533277B2 (en) | 2014-06-25 |
Family
ID=45416803
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010125706A Expired - Fee Related JP5533277B2 (en) | 2010-06-01 | 2010-06-01 | Evaluation method of tire inclusion jumping characteristic, tire inclusion jumping characteristic evaluation apparatus and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5533277B2 (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6045898B2 (en) * | 2012-12-12 | 2016-12-14 | 東洋ゴム工業株式会社 | Tire simulation apparatus, method and program thereof |
JP6093717B2 (en) * | 2014-01-10 | 2017-03-08 | 株式会社トヨタ車体研究所 | Ice fracture model and ice fracture analysis simulation method |
JP6637459B2 (en) * | 2017-03-07 | 2020-01-29 | 公益財団法人鉄道総合技術研究所 | Snow accretion simulation method and snow accretion simulation program |
JP7395456B2 (en) * | 2020-11-18 | 2023-12-11 | 住友重機械工業株式会社 | Simulation device and program |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3643199B2 (en) * | 1997-01-13 | 2005-04-27 | 株式会社 日立インダストリイズ | Vehicle splash test method and apparatus |
JP3133738B2 (en) * | 1998-09-07 | 2001-02-13 | 株式会社ブリヂストン | Tire performance prediction method, fluid simulation method, tire design method, tire vulcanization mold design method, tire vulcanization mold production method, pneumatic tire production method, recording medium recording tire performance prediction program |
JP3655531B2 (en) * | 2000-04-21 | 2005-06-02 | 株式会社ブリヂストン | Tire performance prediction method, fluid simulation method, tire design method, tire vulcanization mold design method, tire vulcanization mold manufacturing method, pneumatic tire manufacturing method, recording medium on which tire performance prediction program is recorded |
JP4152109B2 (en) * | 2002-01-31 | 2008-09-17 | 横浜ゴム株式会社 | Tire characteristic predicting method, tire manufacturing method, and tire characteristic predicting method |
JP4080931B2 (en) * | 2003-03-28 | 2008-04-23 | 住友ゴム工業株式会社 | Tire simulation method |
JP4496716B2 (en) * | 2003-05-19 | 2010-07-07 | 横浜ゴム株式会社 | Tire simulation method, tire performance prediction method, tire manufacturing method, tire and program |
JP2005207775A (en) * | 2004-01-20 | 2005-08-04 | Toyota Motor Corp | Apparatus for analyzing three-dimensional behavior of road surface condition factor |
JP4573018B2 (en) * | 2004-01-22 | 2010-11-04 | トヨタ自動車株式会社 | Vehicle wet road running environment reproduction method and apparatus |
JP5494246B2 (en) * | 2010-06-01 | 2014-05-14 | 横浜ゴム株式会社 | Tire design method, tire design apparatus and program |
-
2010
- 2010-06-01 JP JP2010125706A patent/JP5533277B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2011252748A (en) | 2011-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101513759B1 (en) | Method for simulating tire noise performance and method for manufacturing tire | |
US6564625B1 (en) | Method of designing a tire, method of designing a vulcanizing mold for a tire, method of making a vulcanizing mold for a tire, method of manufacturing a pneumatic tire, and recording medium with tire designing program recorded thereon | |
JP4782392B2 (en) | Tire wear simulation method, apparatus, program, and medium | |
JP4496716B2 (en) | Tire simulation method, tire performance prediction method, tire manufacturing method, tire and program | |
EP1777634A1 (en) | Method for tire rolling simulation on sand | |
JP6293478B2 (en) | Tire wear simulation apparatus, method and program thereof | |
JP2005186900A (en) | Simulation method for tire | |
JP5533277B2 (en) | Evaluation method of tire inclusion jumping characteristic, tire inclusion jumping characteristic evaluation apparatus and program | |
JP3133738B2 (en) | Tire performance prediction method, fluid simulation method, tire design method, tire vulcanization mold design method, tire vulcanization mold production method, pneumatic tire production method, recording medium recording tire performance prediction program | |
EP2650806A2 (en) | Simulation method for a tire | |
JP6045898B2 (en) | Tire simulation apparatus, method and program thereof | |
JPH11201875A (en) | Method for simulating tire performance | |
Hermange et al. | In-depth analysis of hydroplaning phenomenon accounting for tire wear on smooth ground | |
JP4437884B2 (en) | Tire performance prediction method, fluid simulation method, tire design method, recording medium, and tire performance prediction program | |
JP5494246B2 (en) | Tire design method, tire design apparatus and program | |
JP3655531B2 (en) | Tire performance prediction method, fluid simulation method, tire design method, tire vulcanization mold design method, tire vulcanization mold manufacturing method, pneumatic tire manufacturing method, recording medium on which tire performance prediction program is recorded | |
JP4152109B2 (en) | Tire characteristic predicting method, tire manufacturing method, and tire characteristic predicting method | |
JP4594030B2 (en) | Tire performance prediction method, tire simulation method, tire performance prediction program, and recording medium | |
JP2013121800A (en) | Tire performance prediction method, tire simulation method, tire performance prediction program, and recording medium | |
JP2012063949A (en) | Simulation method, system, and program of tire wear | |
JP2006076404A (en) | Tire model, tire behavior simulation method, program, and recording medium | |
JP3452863B2 (en) | Tire drainage simulation method | |
JP5122790B2 (en) | Vehicle running body performance prediction method, apparatus and program | |
JP6025304B2 (en) | Tire performance simulation method, tire performance simulation apparatus, and tire performance simulation program | |
JP4080931B2 (en) | Tire simulation method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130527 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131224 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140107 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140307 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20140401 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5533277 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140414 |
|
LAPS | Cancellation because of no payment of annual fees |