JP6281270B2 - Method and program for creating simulation model - Google Patents

Method and program for creating simulation model Download PDF

Info

Publication number
JP6281270B2
JP6281270B2 JP2013254955A JP2013254955A JP6281270B2 JP 6281270 B2 JP6281270 B2 JP 6281270B2 JP 2013254955 A JP2013254955 A JP 2013254955A JP 2013254955 A JP2013254955 A JP 2013254955A JP 6281270 B2 JP6281270 B2 JP 6281270B2
Authority
JP
Japan
Prior art keywords
mesh
model
space region
tire
model space
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.)
Active
Application number
JP2013254955A
Other languages
Japanese (ja)
Other versions
JP2015114151A (en
Inventor
勇司 児玉
勇司 児玉
正隆 小石
正隆 小石
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yokohama Rubber Co Ltd
Original Assignee
Yokohama Rubber Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yokohama Rubber Co Ltd filed Critical Yokohama Rubber Co Ltd
Priority to JP2013254955A priority Critical patent/JP6281270B2/en
Publication of JP2015114151A publication Critical patent/JP2015114151A/en
Application granted granted Critical
Publication of JP6281270B2 publication Critical patent/JP6281270B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/80Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
    • Y02T10/82Elements for improving aerodynamics

Description

本発明は、タイヤを取り巻く空間中の流体の流れをシミュレートするためのシミュレーションモデルの作成方法、及びこの作成方法を実行するプログラムに関する。   The present invention relates to a simulation model creation method for simulating fluid flow in a space surrounding a tire, and a program for executing the creation method.

従来より、コンピュータを用いた構造体周りの空気等の流体の流れのシミュレーション計算が行われている。このシミュレーション計算では、流体の流れを再現するために、流体空間を多数の領域に分割した、すなわちメッシュ分割したモデルを作成する。現在、タイヤを含めた車両の空力抵抗を低減するためのシミュレーション計算においてもメッシュ分割したシミュレーションモデルが用いられる。   Conventionally, simulation calculation of the flow of fluid such as air around a structure using a computer has been performed. In this simulation calculation, in order to reproduce the fluid flow, a model in which the fluid space is divided into a number of regions, that is, a meshed model is created. Currently, a simulation model divided into meshes is also used in a simulation calculation for reducing the aerodynamic resistance of a vehicle including a tire.

例えば、ゴルフクラブを構造体とし、このゴルフクラブヘッドの空力特性を実物モデルによる試験評価によらずに評価して、効率的にゴルフクラブを開発するためのシミュレーション方法が提案されている(特許文献1)。   For example, a simulation method for efficiently developing a golf club by using a golf club as a structure and evaluating the aerodynamic characteristics of the golf club head without performing a test evaluation using a real model has been proposed (Patent Document). 1).

当該シミュレーション方法では、気流仮想領域を、上記ゴルフボールモデルのボール表面から離間する方向に格子区画の体積が漸次増大するように多数の格子区画で分割する。
また、ゴルフクラブヘッドモデルの表面を、三角形、四角形等の多辺形、又は略三角形、略四角形等の略多辺形の面区画を多数設定し、この各々の面区画を一面とするゴルフクラブヘッドのモデル表面に隣接した格子区画を設定する。ゴルフクラブヘッドのモデル表面に隣接する格子区画は、略四角柱形状等の略多角柱形状や略多角錐形状に設定される。そしてゴルフクラブヘッドのモデル表面に隣接する格子区画からゴルフクラブヘッドから離間する方向に格子区画の体積が漸次増大するように、気流仮想領域の残部を格子状に区画して、気流仮想領域全域を格子区画で区画する。
In the simulation method, the airflow virtual region is divided into a large number of grid sections so that the volume of the grid sections gradually increases in a direction away from the ball surface of the golf ball model.
In addition, the golf club head model has a plurality of polygonal surfaces such as triangles and rectangles, or substantially polygonal surfaces such as substantially triangles and rectangles, and a golf club having each of these surface regions as one surface. Set a grid section adjacent to the model surface of the head. The lattice section adjacent to the model surface of the golf club head is set to a substantially polygonal column shape such as a substantially quadrangular prism shape or a substantially polygonal pyramid shape. The remainder of the airflow virtual region is partitioned in a lattice shape so that the volume of the lattice partition gradually increases in the direction away from the golf club head from the lattice partition adjacent to the model surface of the golf club head. Compartment with lattice section.

また、タイヤの主溝と空力特性との関係を、シミュレーションモデルを用いて解析した研究も知られている(非特許文献1)。
当該研究では、タイヤ周りの流体空間をメッシュ分割したシミュレーションモデルが用いられている。
In addition, a study in which the relationship between the main groove of the tire and the aerodynamic characteristics is analyzed using a simulation model is also known (Non-Patent Document 1).
In this research, a simulation model in which the fluid space around the tire is divided into meshes is used.

特開2012−135344号公報JP 2012-135344 A 木村他、「タイヤ周りの流れ場における縦溝の効果に関する数値研究」,可視化情報シンポジウム講演論文集,p.12−006、2013Kimura et al., "Numerical study on the effect of longitudinal grooves on the flow field around the tire", Proceedings of the Visualization Information Symposium, p. 12-006, 2013

しかし、上述のシミュレーションモデルを、タイヤを含めた車両の空力抵抗のシミュレーション方法に用いる場合、特にタイヤの形状が空力抵抗に及ぼす影響を検討するために、形状の異なる種々のタイヤモデルを作成するが、タイヤモデルを作成するたびに、タイヤモデル周りの空間としてメッシュ分割したシミュレーションモデルを作成することになる。このようなシミュレーションモデルは、極めて多数のメッシュ分割数を含むため、メッシュ分割に要する時間がかかる。したがって、種々のタイヤモデルを用いて空力抵抗のシミュレーションを行う際に要する時間は長大になる。   However, when the simulation model described above is used in a simulation method for aerodynamic resistance of a vehicle including a tire, various tire models having different shapes are created, particularly in order to examine the effect of the tire shape on the aerodynamic resistance. Each time a tire model is created, a simulation model is created by dividing the mesh as a space around the tire model. Since such a simulation model includes a very large number of mesh divisions, it takes time for mesh division. Therefore, the time required for aerodynamic resistance simulation using various tire models becomes long.

そこで、本発明は、タイヤを取り巻く空間中の流体の流れをシミュレートするとき、シミュレーションモデルの作成に要する時間を短縮することができるシミュレーションモデルの作成方法とこの作成方法を実行するプログラムを提供することを目的とする。   Therefore, the present invention provides a simulation model creation method capable of reducing the time required to create a simulation model and a program for executing the creation method when simulating a fluid flow in a space surrounding a tire. For the purpose.

本発明の一つの態様は、タイヤを取り巻く空間中の流体の流れをシミュレートするためのシミュレーションモデルの作成方法である。当該作成方法は、
(1)少なくともタイヤの一部の形状を境界とする第1の空間領域を再現した第1のモデル空間領域と、前記第1の空間領域の周りに位置する第2の空間領域を再現した第2のモデル空間領域を別々にメッシュ分割するステップと、
(2)前記第1のモデル空間領域において作成された第1のメッシュと、前記2のモデル空間領域において作成された第2のメッシュとの間の前記流体の物理量を補間する関係を定めるステップと、含む。
前記第1のモデル空間領域と前記第2のモデル空間領域とは、互いに離間しており、前記第1のモデル空間領域と前記第2のモデル空間領域との離間距離の寸法は、前記第1のメッシュの最小メッシュサイズの寸法よりも小さく、
前記流体の物理量を補間する関係を定める前記ステップにおいて、前記第1のメッシュが前記第2のメッシュに、及び前記第2のメッシュが前記第1のメッシュに物理量を与えるとき、前記第1のメッシュ及び前記第2のメッシュの中で、物理量を受ける側のメッシュAに対向する、前記物理量を与える側のメッシュBが複数ある場合、前記メッシュAに対向する前記メッシュBの境界面の部分の大きさの比率に応じて前記メッシュBの物理量を加重平均した値を前記メッシュAに与えるように、前記関係を定める。
One aspect of the present invention is a method for creating a simulation model for simulating fluid flow in a space surrounding a tire. The creation method is
(1) A first model space region that reproduces a first space region that borders at least a part of the shape of the tire, and a second space region that is located around the first space region. Meshing the two model space regions separately;
(2) defining a relationship for interpolating the physical quantity of the fluid between the first mesh created in the first model space region and the second mesh created in the second model space region; , including the.
The first model space region and the second model space region are separated from each other, and the dimension of the separation distance between the first model space region and the second model space region is the first model space region. Smaller than the minimum mesh size dimension of the mesh,
In the step of determining a relationship for interpolating the physical quantity of the fluid, the first mesh is applied when the first mesh gives the second mesh and the second mesh gives the first mesh a physical quantity. In the second mesh, when there are a plurality of physical quantity-giving meshes B facing the physical quantity receiving mesh A, the size of the boundary surface portion of the mesh B facing the mesh A The relationship is determined so that the mesh A is given a weighted average value of the physical quantities of the mesh B according to the ratio.

このとき、前記第1のメッシュの平均サイズは、前記第2のメッシュの平均サイズに比べて小さい、ことが好ましい。   At this time, it is preferable that the average size of the first mesh is smaller than the average size of the second mesh.

また、前記第1のメッシュは、前記タイヤをモデル化したタイヤモデルと接し、
前記タイヤモデルのタイヤ幅方向における、前記タイヤモデルの幅寸法をSt[mm]とし、前記タイヤモデルの直径寸法をDt[mm]としたとき、前記第1のモデル空間領域の前記タイヤ幅方向の幅寸法S[mm]は、前記幅寸法Stの0.1倍より大きく3倍より小さい寸法であり、かつ、前記第1の空間領域の前記タイヤ径方向に沿った寸法の最大値は、前記直径寸法Dtの0.3倍より大きく2倍より小さい寸法である、ことが好ましい。
The first mesh is in contact with a tire model that models the tire,
When the width dimension of the tire model in the tire width direction of the tire model is St [mm] and the diameter dimension of the tire model is Dt [mm], the tire width direction of the first model space region is the tire width direction. The width dimension S [mm] is a dimension larger than 0.1 times and smaller than three times the width dimension St, and the maximum value of the dimension along the tire radial direction of the first space region is It is preferable that the dimension is larger than 0.3 times the diameter dimension Dt and smaller than twice.

前記第1のモデル空間領域は、前記タイヤのサイドウォールと接する空間を再現したモデル空間領域を含む、ことが好ましい。   It is preferable that the first model space region includes a model space region that reproduces a space in contact with a sidewall of the tire.

また、前記第1のモデル空間領域は、前記タイヤのトレッドパターンと接する空間を再現したモデル空間領域を含む、ことも同様に好ましい。   It is also preferable that the first model space region includes a model space region that reproduces a space in contact with the tread pattern of the tire.

前記タイヤがタイヤ回転軸の周りに回転することを再現するために、前記第1のメッシュは、タイヤ回転軸の周りに回転するメッシュである、ことが好ましい。   In order to reproduce the rotation of the tire around the tire rotation axis, the first mesh is preferably a mesh that rotates around the tire rotation axis.

本発明のさらに一つの態様は、タイヤを取り巻く空間中の流体の流れをシミュレートするためのシミュレーションモデルの作成方法をコンピュータに実行させるコンピュータが実行可能なプログラムである。当該プログラムは、
(3)少なくともタイヤの一部の形状を境界とする第1の空間領域を再現した第1のモデル空間領域と、前記第1の空間領域の周りに位置する第2の空間領域を再現した第2のモデル空間領域を別々にメッシュ分割する処理をコンピュータに行わせる手順と、
(4)前記第1のモデル空間領域において作成された第1のメッシュと、前記2のモデル空間領域において作成された第2のメッシュとの間の前記流体の物理量を補間する関係をコンピュータに定めさせる手順と、含む。
前記第1のモデル空間領域と前記第2のモデル空間領域とは、互いに離間しており、前記第1のモデル空間領域と前記第2のモデル空間領域との離間距離の寸法は、前記第1のメッシュの最小メッシュサイズの寸法よりも小さく、
前記流体の物理量を補間する関係をコンピュータに定めさせる前記手順において、前記第1のメッシュが前記第2のメッシュに、及び前記第2のメッシュが前記第1のメッシュに物理量を与えるとき、前記第1のメッシュ及び前記第2のメッシュの中で、物理量を受ける側のメッシュAに対向する、前記物理量を与える側のメッシュBが複数ある場合、前記メッシュAに対向する前記メッシュBの境界面の部分の大きさの比率に応じて前記メッシュBの物理量を加重平均した値を前記メッシュAに与えるように、前記関係を定めさせる。
Yet another embodiment of the present invention is a computer-executable program that causes a computer to execute a simulation model creation method for simulating fluid flow in a space surrounding a tire. The program is
(3) A first model space region that reproduces a first space region that borders at least a part of the shape of the tire, and a second space region that is located around the first space region. A procedure for causing a computer to perform the process of dividing the model space area of the two meshes separately;
(4) A relationship for interpolating the physical quantity of the fluid between the first mesh created in the first model space region and the second mesh created in the second model space region is determined in the computer. including a procedure for, the.
The first model space region and the second model space region are separated from each other, and the dimension of the separation distance between the first model space region and the second model space region is the first model space region. Smaller than the minimum mesh size dimension of the mesh,
In the step of causing the computer to determine a relationship for interpolating the physical quantity of the fluid, when the first mesh gives a physical quantity to the second mesh and the second mesh gives the physical quantity to the first mesh, the first mesh In the case where there are a plurality of meshes B on the physical quantity-given side facing the mesh A on the physical quantity receiving side among the one mesh and the second mesh, the boundary surface of the mesh B facing the mesh A The relationship is determined so that the mesh A is given a weighted average value of the physical quantities of the mesh B according to the ratio of the size of the parts.

上述の態様のシミュレーションモデルの作成方法及びプログラムによれば、シミュレーションモデルの作成に要する時間を短縮することができる。   According to the simulation model creation method and program of the above-described aspect, the time required to create the simulation model can be shortened.

本実施形態のシミュレーションモデルの作成方法を実行するモデル作成装置のブロック構成図である。It is a block block diagram of the model creation apparatus which performs the creation method of the simulation model of this embodiment. 本実施形態のシミュレーションモデルの作成方法の流れを説明する図である。It is a figure explaining the flow of the preparation method of the simulation model of this embodiment. 本実施形態のシミュレーションモデルの作成方法で作成された第1のモデル空間領域と第2のモデル空間領域における第1のメッシュ及び第2のメッシュの一例を示す図である。It is a figure which shows an example of the 1st mesh in the 1st model space area | region and the 2nd model space area which were created with the creation method of the simulation model of this embodiment, and a 2nd mesh. 本実施形態のシミュレーションモデルの作成方法で作成されたシミュレーションモデルの一例を示す図である。It is a figure which shows an example of the simulation model created with the creation method of the simulation model of this embodiment. (a),(b)は、本実施形態のシミュレーションモデルの作成方法で作成される第1のモデル空間領域と第2のモデル空間領域間の補間の関係の例を説明する図である。(A), (b) is a figure explaining the example of the relationship of the interpolation between the 1st model space area | region created with the creation method of the simulation model of this embodiment, and a 2nd model space area | region. 本実施形態のシミュレーションモデルの作成方法で用いる空間領域と車両モデルの一例を示す図である。It is a figure which shows an example of the space area and vehicle model which are used with the preparation method of the simulation model of this embodiment.

以下、本実施形態のシミュレーションモデルの作成方法及びプログラムについて、図面を参照して詳細に説明する。なお、本実施形態で用いる空間領域をメッシュ分割したシミュレーションモデルは、有限要素法、有限差分法、有限体積法、有限メッシュ法、格子ボルツマン法のシミュレーション計算に適用できるシミュレーションモデルである。   Hereinafter, a simulation model creation method and program according to the present embodiment will be described in detail with reference to the drawings. Note that the simulation model obtained by dividing the space region used in this embodiment into a mesh is a simulation model that can be applied to simulation calculations of the finite element method, the finite difference method, the finite volume method, the finite mesh method, and the lattice Boltzmann method.

図1は、本実施形態のシミュレーションモデルの作成方法を実行するモデル作成装置10のブロック構成図である。モデル作成装置10は、CPU12、メモリ14、ROM16、及び入出力部18を有するコンピュータにより構成されている。ROM16には、本実施形態であるシミュレーションモデルの作成方法を実施するプログラムが記憶されており、CPU12の指示により呼び出して実行することにより、領域設定部20、メッシュ分割部22、補間関係設定部24、及び制御部26を形成する。
領域設定部20は、オペレータの指示に従って第1のモデル空間領域と第2のモデル空間領域とを設定する。
メッシュ分割部22は、第1のモデル空間領域と第2のモデル空間領域を別々にメッシュ分割する。
補間関係設定部24は、第1のモデル空間領域において作成された第1のメッシュと、2のモデル空間領域において作成された第2のメッシュとの間の流体の物理量を補間する関係を設定する。この設定は、オペレータの指示に従って行われる。
制御部26は、領域設定部20、メッシュ分割部22、補間関係設定部24の動作を制御し、あるいは、判定を行う。
入出力部18には、マウス・キーボード等の入力操作デバイス27、ディスプレイ28、及びプリンタ30と接続されている。
メモリ14には、領域設定部20、メッシュ分割部22、及び補間関係設定部24で得られた処理結果や、処理前の設定条件を記憶する。
FIG. 1 is a block configuration diagram of a model creation apparatus 10 that executes the simulation model creation method of the present embodiment. The model creation device 10 is configured by a computer having a CPU 12, a memory 14, a ROM 16, and an input / output unit 18. The ROM 16 stores a program for executing the simulation model creation method according to the present embodiment, and is called and executed by an instruction from the CPU 12 to execute an area setting unit 20, a mesh dividing unit 22, and an interpolation relation setting unit 24. And the control unit 26 are formed.
The area setting unit 20 sets a first model space area and a second model space area in accordance with an instruction from the operator.
The mesh division unit 22 meshes the first model space region and the second model space region separately.
The interpolation relationship setting unit 24 sets a relationship for interpolating the physical quantity of the fluid between the first mesh created in the first model space region and the second mesh created in the second model space region. . This setting is performed in accordance with an operator instruction.
The control unit 26 controls the operation of the region setting unit 20, the mesh division unit 22, and the interpolation relation setting unit 24, or makes a determination.
The input / output unit 18 is connected to an input operation device 27 such as a mouse / keyboard, a display 28, and a printer 30.
The memory 14 stores the processing results obtained by the region setting unit 20, the mesh division unit 22, and the interpolation relation setting unit 24, and setting conditions before processing.

図2は、モデル作成装置10が実施するシミュレーションモデルの作成方法の流れを説明する図である。本実施形態のシミュレーションモデルは、タイヤを取り巻く空間中の流体の流れをシミュレートするためのモデルである。
まず、領域設定部20は、タイヤの形状を設定する(ステップS10)。タイヤの形状は3次元形状であり、例えば、トレッドパターン、サイドウォールの表面の凹凸、サイドウォールの断面形状を忠実に再現した形状である。このタイヤの形状は、メモリ14に記憶されており、メモリ14から呼び出される。あるいは、入力操作デバイス27によって形状を設定することもできる。
FIG. 2 is a diagram for explaining the flow of a simulation model creation method performed by the model creation apparatus 10. The simulation model of this embodiment is a model for simulating the flow of fluid in the space surrounding the tire.
First, the region setting unit 20 sets the tire shape (step S10). The shape of the tire is a three-dimensional shape, for example, a shape that faithfully reproduces the tread pattern, the irregularities on the surface of the sidewall, and the sectional shape of the sidewall. The shape of the tire is stored in the memory 14 and is called from the memory 14. Alternatively, the shape can be set by the input operation device 27.

次に、メッシュ分割部22は、タイヤを取り巻く空間において、少なくともタイヤの一部の形状を境界とする第1の空間領域を再現した第1のモデル空間領域30と、第1の空間領域の周りに位置する第2の空間領域を再現した第2のモデル空間領域32を作成して別々にメッシュ分割する(ステップS20)。メッシュ分割された第1のモデル空間領域30と第2のモデル空間領域32はメモリ14に記憶される。
なお、第1の空間領域には、タイヤは含まれない。したがって、第1のモデル空間領域30は、タイヤの形状を境界面によって再現した空間領域である。上述の「タイヤの少なくとも一部」には、トレッドパターンを有するトレッド部、あるいは、サイド部が含まれ、タイヤ全体が含まれてもよい。なお、第1のモデル空間領域30によって作成される第1のメッシュにおいて、メッシュの格子点が、タイヤの形状の表面上、あるいは、この表面上よりも僅かにタイヤ表面外側に位置し、これらの格子点のうちタイヤの形状に近い2つの格子点により作られる直線あるいは、少なくとも3つの格子点によ作られる平面の一部がタイヤの形状の内側を通ってもよい。メッシュ分割では、3次元の空間領域が、例えば4面体、5面体、あるいは6面体等の要素によって細かく分割される。メッシュ分割の方法については特に制限されず、公知のメッシュ分割法が用いられる。
図3は、別々に作成された第1のモデル空間領域30と第2のモデル空間領域32における第1のメッシュ及び第2のメッシュの一例を示す図である。第2のモデル空間領域32には、第1のモデル空間領域30の大きさと外周形状に概略対応した凹部38を有する。
Next, the mesh division unit 22 includes a first model space region 30 that reproduces a first space region that has at least a part of the shape of the tire as a boundary in a space surrounding the tire, and around the first space region. A second model space region 32 that reproduces the second space region positioned at is created and meshed separately (step S20). The first model space region 30 and the second model space region 32 that are divided into meshes are stored in the memory 14.
Note that the first space region does not include a tire. Therefore, the first model space region 30 is a space region in which the shape of the tire is reproduced by the boundary surface. The “at least part of the tire” described above includes a tread portion having a tread pattern or a side portion, and may include the entire tire. Note that in the first mesh created by the first model space region 30, the mesh lattice points are located on the surface of the tire shape, or slightly outside the tire surface, and these linear created by two grid points close to the shape of the tire among the lattice points or a part of the plane created Ri by the at least three grid points may be through the inner shape of the tire. In mesh division, a three-dimensional space region is finely divided by elements such as a tetrahedron, a pentahedron, or a hexahedron. The mesh dividing method is not particularly limited, and a known mesh dividing method is used.
FIG. 3 is a diagram illustrating an example of the first mesh and the second mesh in the first model space region 30 and the second model space region 32 created separately. The second model space region 32 has a recess 38 that roughly corresponds to the size and outer peripheral shape of the first model space region 30.

次に、補間関係設定部24は、第1のモデル空間領域30を第2のモデル空間領域32に組み込んで1つの空間領域を形成するとき、第1のモデル空間領域30と第2のモデル空間領域32の近接するメッシュ領域間で物理量を適切に受け渡すように補間の関係を設定する(ステップS30)。第1のモデル空間領域30あるいは第2のモデル空間領域32の互いに面する境界面を有するメッシュ領域のそれぞれについて、第2のモデル空間領域32あるいは第1のモデル空間領域30のメッシュ領域との間で、物理量の受け渡しを行うための物理量の補間関係が定められる。シミュレーション計算をするとき、物理量が第1のモデル空間領域30と第2のモデル空間領域32の接続部分で連続的に繋がることが、適切なシミュレーション結果を得る上で必要である。このため、第1のモデル空間領域30において作成された第1のメッシュと、第2のモデル空間領域32において作成された第2のメッシュとの間の流体の物理量を補間する関係が定められる。ここで、流体の物理量は、速度、圧力、密度、あるいは温度を含む。流体の物理量は、互いに近接するメッシュ領域内の平均物理量やメッシュ領域の重心位置における物理量が、第1のモデル空間領域30と第2のモデル空間領域32の近接するメッシュ領域間で受け渡される。
これにより、補間関係設定部24は、別々に作成された第1のモデル空間領域30を第2のモデル空間領域32に組み込んだ1つのモデル空間領域を有するシミュレーションモデルを作成する(ステップS40)。なお、図4は、作成された1つのシミュレーションモデルの一例を示す図である。
Next, when the interpolation relationship setting unit 24 incorporates the first model space region 30 into the second model space region 32 to form one space region, the first model space region 30 and the second model space Interpolation relations are set so that physical quantities are appropriately transferred between mesh areas adjacent to the area 32 (step S30). Between each of the mesh regions having the boundary surfaces facing each other in the first model space region 30 or the second model space region 32, between the mesh region of the second model space region 32 or the first model space region 30 Thus, the physical quantity interpolation relationship for transferring the physical quantity is defined. When performing the simulation calculation, it is necessary for the physical quantity to be continuously connected at the connecting portion between the first model space region 30 and the second model space region 32 in order to obtain an appropriate simulation result. For this reason, the relationship which interpolates the physical quantity of the fluid between the 1st mesh created in the 1st model space area | region 30 and the 2nd mesh created in the 2nd model space area | region 32 is defined. Here, the physical quantity of the fluid includes speed, pressure, density, or temperature. As for the physical quantity of the fluid, the average physical quantity in the mesh area adjacent to each other and the physical quantity at the center of gravity of the mesh area are transferred between the mesh areas adjacent to each other in the first model space area 30 and the second model space area 32.
Thereby, the interpolation relation setting unit 24 creates a simulation model having one model space area in which the separately created first model space area 30 is incorporated into the second model space area 32 (step S40). FIG. 4 is a diagram illustrating an example of one created simulation model.

図5(a),(b)は、第1のモデル空間領域30と第2のモデル空間領域32間の補間の関係の例を説明する図である。
図5(a)に示す例では、第1のモデル空間領域30と第2のモデル空間領域32の対向する境界に面する第1のメッシュ領域34及び第2のメッシュ領域36を補間対象のメッシュ領域とする。対向するメッシュ領域の大きさは同じサイズとし、メッシュ領域の境界面が互いに対向している。このとき、補間関係設定部24は、第1のメッシュ領域34に対して最も近くに位置する対向する第2のメッシュ領域36を抽出し、この第2のメッシュ領域36に、第1のメッシュ領域34における物理量を与える。あるいは、第1のメッシュ領域34に対して最も近くに位置する対向する第2のメッシュ領域36を抽出し、この第1のメッシュ領域34に、第2のメッシュ領域36における物理量を与える。すなわち、第1のメッシュ領域34と第2のメッシュ領域36における物理量を同一とする。
ここで、メッシュ領域の物理量は、メッシュ領域内の平均物理量あるいはメッシュ領域の重心位置における物理量が例示される。
FIGS. 5A and 5B are diagrams illustrating an example of an interpolation relationship between the first model space region 30 and the second model space region 32. FIG.
In the example shown in FIG. 5A, the first mesh area 34 and the second mesh area 36 facing the boundary between the first model space area 30 and the second model space area 32 are interpolated. This is an area. The sizes of the mesh regions facing each other are the same size, and the boundary surfaces of the mesh regions are opposed to each other. At this time, the interpolation relationship setting unit 24 extracts the second mesh region 36 facing the closest to the first mesh region 34, and the first mesh region 36 is extracted to the second mesh region 36. The physical quantity at 34 is given. Alternatively, the opposing second mesh region 36 located closest to the first mesh region 34 is extracted, and a physical quantity in the second mesh region 36 is given to the first mesh region 34. That is, the physical quantities in the first mesh region 34 and the second mesh region 36 are the same.
Here, the physical quantity of the mesh area is exemplified by an average physical quantity in the mesh area or a physical quantity at the center of gravity of the mesh area.

図5(b)に示す例でも、第1のモデル空間領域30と第2のモデル空間領域32の対向する境界に面する第1のメッシュ領域34及び第2のメッシュ領域36を補間対象のメッシュ領域とする。図5(b)に示す例では、対向するメッシュ領域の大きさは異なり、メッシュ領域の対向する境界面が互いにタイヤの形状の周方向において位置ずれしている。図(b)に示す例では、第1のメッシュ領域34aと第2のメッシュ領域34bは周方向に互いに隣接し、第2のメッシュ領域36aと第2のメッシュ領域36bは周方向に互いに隣接している。さらに、第1のメッシュ領域34a,34bの境界面のタイヤの形状の周方向に沿ったサイズは、第2のメッシュ領域36a,36bの境界面のタイヤの形状の周方向に沿ったサイズに比べて大きい。このような場合、第2のメッシュ領域36aの境界面全ては、第1のメッシュ領域34aの境界面に対向しているので、第2のメッシュ領域36aには、第1のメッシュ領域34aにおける物理量を与える。第2のメッシュ領域36bの境界面は、第1のメッシュ領域34aの境界面と第1のメッシュ領域34bの境界面にまたがって対向しているので、第1のメッシュ領域34aにおける物理量と第2のメッシュ領域34bにおける物理量の平均値を与える。この場合、上記平均値として、以下に記載する加重平均値を用いることが好ましい。加重平均値に用いる重み付け係数は、第2のメッシュ領域36bの境界面に対向する第1のメッシュ領域34aの境界面の部分の大きさと、第2のメッシュ領域36bの境界面に対向する第1のメッシュ領域34bの境界面の部分の大きさの比率に応じて定まる係数であり、上記比率が大きいほど重み付け係数は大きくなる。
同様に、第1のメッシュ領域34aに物理量を与えるとき、第2のメッシュ領域36aにおける物理量と第2のメッシュ領域36bにおける物理量の平均値を与える。この平均値は、上述した加重平均値を用いることが好ましい。
また、第1のモデル空間領域30と第2のモデル空間領域32間の物理量の補間について、メッシュ領域の位置情報あるいは格子点の位置情報を用いて物理量に線形外挿補間を適用してもよく、あるいは、周知のパラメトリック空間補間を適用することができる。パラメトリック空間補間は、例えば特許4548005号公報(段落0026〜0031)に記載されている。パラメトリック空間補間とは、物理量の補間をしようとする点を含むメッシュ領域を、例えば矩形形状等の単位メッシュ領域から形状関数を用いて形状変換を施したものとして定めることにより、物理量を補完しようとするメッシュ領域内の点に対応した、単位メッシュ領域の形状内における対応点の位置情報を求め、この位置情報を用いて物理量を補完する処理である。
Also in the example shown in FIG. 5B, the first mesh area 34 and the second mesh area 36 facing the boundary between the first model space area 30 and the second model space area 32 are interpolated. This is an area. In the example shown in FIG. 5 (b), the sizes of the mesh regions facing each other are different, and the opposing boundary surfaces of the mesh regions are displaced from each other in the circumferential direction of the tire shape. In the example shown in FIG. 5 (b), the first mesh region 34a and the second mesh region 34b are adjacent to each other in the circumferential direction, the second mesh region 36a and the second mesh region 36b adjacent to each other in the circumferential direction doing. Furthermore, the size along the circumferential direction of the tire shape at the boundary surface of the first mesh regions 34a and 34b is larger than the size along the circumferential direction of the tire shape at the boundary surface of the second mesh regions 36a and 36b. Big. In such a case, since all the boundary surfaces of the second mesh region 36a are opposed to the boundary surface of the first mesh region 34a, the physical quantity in the first mesh region 34a is included in the second mesh region 36a. give. Since the boundary surface of the second mesh region 36b faces the boundary surface of the first mesh region 34b across the boundary surface of the first mesh region 34b, the physical quantity in the first mesh region 34a The average value of the physical quantities in the mesh region 34b is given. In this case, it is preferable to use a weighted average value described below as the average value. The weighting coefficient used for the weighted average value is the size of the boundary surface portion of the first mesh region 34a facing the boundary surface of the second mesh region 36b and the first surface facing the boundary surface of the second mesh region 36b. This coefficient is determined according to the ratio of the size of the boundary surface portion of the mesh area 34b, and the weighting coefficient increases as the ratio increases.
Similarly, when a physical quantity is given to the first mesh area 34a, an average value of the physical quantity in the second mesh area 36a and the physical quantity in the second mesh area 36b is given. The average value is preferably the above-described weighted average value.
In addition, for the interpolation of the physical quantity between the first model space area 30 and the second model space area 32, linear extrapolation interpolation may be applied to the physical quantity using the mesh area position information or the grid point position information. Alternatively, well-known parametric spatial interpolation can be applied. Parametric spatial interpolation is described in, for example, Japanese Patent No. 4548005 (paragraphs 0026 to 0031). Parametric spatial interpolation is to try to complement the physical quantity by defining the mesh area including the point to interpolate the physical quantity as a shape converted from a unit mesh area such as a rectangular shape using a shape function. In this process, the position information of the corresponding point in the shape of the unit mesh area corresponding to the point in the mesh area to be obtained is obtained, and the physical quantity is complemented using this position information.

このようにして、物理量の補間関係を定めることにより、図4に示すような1つのシミュレーションモデルが作成される。このようなシミュレーションモデルにおいて、流体の流れが再現される。
さらに、タイヤの形状を変更して他のシミュレーションモデルを作る場合、まず、領域設定部20は、変更したタイヤの形状を設定する(ステップS50)。
変更したタイヤの形状は、メモリ14から呼び出されてもよく、また、マウス・キーボード等の入力操作デバイス27による入力により設定されてもよい。
In this way, by defining the physical quantity interpolation relationship, one simulation model as shown in FIG. 4 is created. In such a simulation model, the fluid flow is reproduced.
Furthermore, when making another simulation model by changing the shape of the tire, first, the region setting unit 20 sets the changed shape of the tire (step S50).
The changed tire shape may be called from the memory 14 or may be set by an input using an input operation device 27 such as a mouse / keyboard.

次に、メッシュ分割部22は、設定されたタイヤの形状を有するタイヤを取り巻く空間において、少なくともタイヤの一部の形状を境界とする第1の空間領域を再現した第1のモデル空間領域30を再度設定し、メッシュ分割する。メッシュ分割された第1のモデル空間領域30の外周は、ステップS20においてメッシュ分割された第2のモデル空間領域32に組み込まれるように、第2のモデル空間領域32の凹部38の大きさに対応している。ステップS60では、第2のモデル空間領域32は新たに作成されない。ステップS50で作成される第1のモデル空間領域30は、変更したタイヤの形状に合わせて作成され、外周は、第2のモデル空間領域32の凹部38の大きさに対応している。すなわち、メッシュ分割部22は、第1のモデル空間領域30のみをメッシュ分割するので、モデル作成時間を大きく低減することができる。メッシュ分割された第2のモデル空間領域32は、ステップS20で作成されたものを用いる。   Next, the mesh division unit 22 includes a first model space region 30 that reproduces a first space region that has at least a part of the shape of the tire as a boundary in a space surrounding the tire having the set shape of the tire. Set again and divide the mesh. The outer circumference of the mesh-divided first model space area 30 corresponds to the size of the concave portion 38 of the second model space area 32 so as to be incorporated into the second model space area 32 mesh-divided in step S20. doing. In step S60, the second model space region 32 is not newly created. The first model space region 30 created in step S50 is created in accordance with the changed tire shape, and the outer periphery corresponds to the size of the recess 38 of the second model space region 32. That is, since the mesh dividing unit 22 meshes only the first model space region 30, the model creation time can be greatly reduced. As the second model space region 32 divided by mesh, the one created in step S20 is used.

次に、補間関係設定部24は、ステップS60で作成された、メッシュ分割された第1のモデル空間領域30と、ステップS20で作成され、メモリ14に記憶されたメッシュ分割された第2のモデル空間領域32とを用いて、物理量の補間関係を設定する(ステップS70)。この補間関係の設定は、ステップS30における物理量の補間関係の設定と同じ方法で行う。このため、ステップS70の補間関係の設定の説明は省略する。
物理量の補間関係を境界面上に位置するメッシュ領域のそれぞれについて設定することにより、補間関係設定部24は、第1のモデル空間領域11を第2のモデル空間領域32に組み込んだシミュレーションモデルを作成する(ステップS80)。こうして、形状を変更したタイヤ周りの空間を再現したシミュレーションモデルが得られる。
Next, the interpolating relationship setting unit 24 creates the first model space area 30 divided in mesh created in step S60 and the second model divided in mesh created in step S20 and stored in the memory 14. An interpolating relationship between physical quantities is set using the space region 32 (step S70). The setting of the interpolation relationship is performed by the same method as the setting of the physical quantity interpolation relationship in step S30. For this reason, the description of the interpolation relation setting in step S70 is omitted.
By setting the physical quantity interpolation relationship for each mesh region located on the boundary surface, the interpolation relationship setting unit 24 creates a simulation model in which the first model space region 11 is incorporated into the second model space region 32. (Step S80). In this way, a simulation model that reproduces the space around the tire whose shape has been changed is obtained.

さらに、制御部26は、タイヤの形状を変更したタイヤ周りの空間を再現したシミュレーションモデルを作成するか否かを判定する(ステップS)。制御部26は、入力操作デバイス27からのオペレータの入力指示に応じて、ステップS60〜S80を繰り返すステップに進むか、あるいは、シミュレーションモデルを用いた計算を行う(ステップS100)ステップに進む。シミュレーション計算を行うとき、制御部26は、例えばROM16に記憶されているソフトウェアを呼び出して起動し、シミュレーション計算を行う機能を有するソフトウェアモジュールを立ち上げる。このソフトウェアモジュールは、ステップS40及びステップS80で作成された複数のシミュレーションモデルを用いてシミュレーション計算を行う。シミュレーション計算は、例えば、有限要素法、有限差分法、有限体積法、有限メッシュ法、あるいは格子ボルツマン法を用いて行う。   Furthermore, the control unit 26 determines whether or not to create a simulation model that reproduces the space around the tire whose shape is changed (step S). The control unit 26 proceeds to a step of repeating steps S60 to S80 or performs a calculation using a simulation model (step S100) according to an input instruction from the input operation device 27 by the operator. When performing the simulation calculation, for example, the control unit 26 calls and activates software stored in the ROM 16 and starts up a software module having a function of performing the simulation calculation. This software module performs a simulation calculation using the plurality of simulation models created in step S40 and step S80. The simulation calculation is performed using, for example, a finite element method, a finite difference method, a finite volume method, a finite mesh method, or a lattice Boltzmann method.

このように、タイヤの形状を変更して、空間領域のシミュレーションモデルを作成するとき、タイヤの形状を変更する度に、先に作成してメモリ14に記憶した第2のモデル空間領域32を用いるので、第1のモデル空間領域30だけをメッシュ分割すればよい。したがって、シミュレーションモデルの作成に要する時間を短縮することができる。   As described above, when the simulation model of the space region is created by changing the shape of the tire, the second model space region 32 previously created and stored in the memory 14 is used every time the shape of the tire is changed. Therefore, only the first model space region 30 needs to be mesh-divided. Therefore, the time required for creating the simulation model can be shortened.

なお、第1のモデル空間領域30における第1のメッシュの平均サイズは、第2のモデル空間領域32における第2のメッシュの平均サイズに比べて小さいことが、タイヤ周りの流体を正確に再現する点で、好ましい。ここで、平均サイズとは、メッシュ領域の平均体積、メッシュ領域の平均要素長さ、あるいは、各メッシュ領域の要素長さの総和の平均長さを含む。「平均要素長さ」とは、各メッシュ領域を作る線分の長さの平均を、第1のモデル空間領域30あるいは第2のモデル空間領域32の全メッシュ領域において平均したものである。「各メッシュ領域の要素長さの総和の平均長さ」とは、各メッシュ領域を作る線分の長さの総和を、第1のモデル空間領域30あるいは第2のモデル空間領域32の全メッシュ領域において平均したものである。   The average size of the first mesh in the first model space region 30 is smaller than the average size of the second mesh in the second model space region 32, so that the fluid around the tire is accurately reproduced. In terms, it is preferable. Here, the average size includes the average volume of the mesh region, the average element length of the mesh region, or the average length of the sum of the element lengths of each mesh region. The “average element length” is an average of the lengths of line segments forming each mesh region in all mesh regions of the first model space region 30 or the second model space region 32. “The average length of the sum of the element lengths of each mesh area” means the total length of the line segments that make up each mesh area, the total mesh of the first model space area 30 or the second model space area 32 It is averaged over the area.

また、第1のモデル空間領域30の第1のメッシュは、タイヤをモデル化したタイヤモデルと接し、タイヤモデルのタイヤ幅方向における、タイヤモデルの幅寸法をSt[mm]とし、タイヤモデルのタイヤ径方向の直径寸法をDt[mm]とする。このとき、第1のモデル空間領域30のタイヤ幅方向の幅寸法S[mm]は、タイヤモデルの幅寸法Stの0.1倍より大きく3倍より小さい寸法であり、かつ、第1のモデル空間領域30のタイヤ径方向の最大値は、直径寸法Dtの0.3倍より大きく2倍より小さい寸法であることが、第1のモデル空間領域30のメッシュ分割を効率よく行ってシミュレーションモデルの作成時間を短縮する点で好ましい。ここで、タイヤモデルは、メッシュ分割されたシミュレーション計算可能なメッシュモデルではなく、剛体モデルであり、第1のモデル空間領域30の境界を画するものである。したがって、シミュレーション計算時、上記剛体モデルを画する第1のモデル空間領域30の境界面において境界条件等が設定される。   The first mesh in the first model space region 30 is in contact with the tire model that models the tire, and the tire model width dimension in the tire width direction of the tire model is St [mm]. The diameter dimension in the radial direction is Dt [mm]. At this time, the width dimension S [mm] in the tire width direction of the first model space region 30 is a dimension larger than 0.1 times and smaller than three times the width dimension St of the tire model, and the first model. The maximum value in the tire radial direction of the space region 30 is a size that is larger than 0.3 times the diameter size Dt and smaller than twice, so that the mesh model of the first model space region 30 is efficiently divided and the simulation model This is preferable in terms of shortening the creation time. Here, the tire model is not a mesh model that can be calculated by simulation, but a rigid body model, and demarcates the boundary of the first model space region 30. Therefore, at the time of simulation calculation, boundary conditions and the like are set on the boundary surface of the first model space region 30 that defines the rigid body model.

第1のモデル空間領域30と第2のモデル空間領域32とは、一部分が重なっていてもよく、あるいは、互いに離間してもよい。但し、互いに離間する場合、第1のモデル空間領域30と第2のモデル空間領域32との離間距離の寸法は、第1のモデル空間領域30の第1のメッシュの最小メッシュサイズの寸法よりも小さいことが、シミュレーション計算において計算結果の精度を確保する点で好ましく、物理量の補間の精度を保持する点で好ましい。ここで、最小メッシュのサイズとは、例えば各メッシュ領域を作る線分の長さの平均の中の最小長さを意味する。第1のモデル空間領域30と第2のモデル空間領域32との離間距離は、図5(a)に示す例では図中の距離Dである。   The first model space region 30 and the second model space region 32 may partially overlap or may be separated from each other. However, when separated from each other, the dimension of the separation distance between the first model space area 30 and the second model space area 32 is larger than the minimum mesh size of the first mesh in the first model space area 30. A small value is preferable in terms of ensuring the accuracy of calculation results in simulation calculations, and is preferable in terms of maintaining the accuracy of physical quantity interpolation. Here, the size of the minimum mesh means, for example, the minimum length among the average lengths of the line segments forming each mesh region. The separation distance between the first model space region 30 and the second model space region 32 is a distance D in the drawing in the example shown in FIG.

第1のモデル空間領域30は、例えば、タイヤのサイドウォールと接する空間を再現したモデル空間領域を含むことができる。タイヤのサイドウォールの表面凹凸の形状を種々変更して空力抵抗に優れたタイヤの形状を抽出する際、シミュレーションモデルの作成時間を短縮できる点で効果は大きい。
また、第1のモデル空間領域30は、例えば、タイヤのトレッドパターンと接する空間を再現したモデル空間領域を含むことができる。タイヤのトレッドパターンやトレッドプロファイルを種々変更して空力抵抗に優れたタイヤを抽出する際、シミュレーションモデルの作成時間を短縮できる点で効果は大きい。特に、トレッドパターン近傍のメッシュ領域は、正確にトレッドパターンの形状を再現する点から複雑なメッシュ分割がなされている。このため、トレッドパターンを種々変更するとき、第1のモデル空間領域30のみを細かくメッシュ分割すればよいので、シミュレーションモデルの作成時間を大幅に短縮することができる。
ここで、第1のモデル空間領域30を作成するときに参照するタイヤの形状は、タイヤの設計図面やCAD情報から得ることができる。また、有限要素法などを用いたタイヤの動的又は静的変形解析によって算出されるタイヤの外側形状の一部から得ることもできる。タイヤの動的又は静的変形解析によって算出されるタイヤの外側形状の一部を用いることにより、タイヤが自動車に装着され接地した場合のタイヤの変形を考慮することができ、タイヤの設計図面やCAD情報を用いる場合よりもシミュレーションの精度が向上する効果がある。また、タイヤが自動車に装着され接地した場合のタイヤの変形した形状は、上記有限要素法などを用いたタイヤの動的又は静的変形解析によって算出された結果から取得する他、自動車などの車両に取り付けられた実際のタイヤの変形した形状をレーザ形状計測装置などによって計測して取得することもできる。
The first model space region 30 can include, for example, a model space region that reproduces a space in contact with a tire sidewall. When extracting the shape of a tire having excellent aerodynamic resistance by variously changing the shape of the surface unevenness of the sidewall of the tire, the effect is great in that the time for creating a simulation model can be shortened.
The first model space region 30 can include, for example, a model space region that reproduces a space in contact with the tire tread pattern. When extracting tires with excellent aerodynamic resistance by changing the tread pattern and tread profile of the tires, the effect of reducing the time required to create a simulation model is significant. In particular, the mesh region in the vicinity of the tread pattern is subjected to complicated mesh division in terms of accurately reproducing the shape of the tread pattern. For this reason, when the tread pattern is variously changed, only the first model space region 30 needs to be finely divided into meshes, so that the simulation model creation time can be greatly shortened.
Here, the shape of the tire referred to when creating the first model space region 30 can be obtained from the tire design drawing or CAD information. It can also be obtained from a part of the outer shape of the tire calculated by dynamic or static deformation analysis of the tire using a finite element method or the like. By using a part of the outer shape of the tire calculated by dynamic or static deformation analysis of the tire, it is possible to take into account the deformation of the tire when the tire is mounted on an automobile and grounded. There is an effect that the accuracy of simulation is improved as compared with the case of using CAD information. In addition, the deformed shape of the tire when the tire is mounted on the vehicle and grounded is obtained from the result calculated by the dynamic or static deformation analysis of the tire using the finite element method or the like, and the vehicle such as the vehicle. The deformed shape of the actual tire attached to the tire can also be obtained by measuring with a laser shape measuring device or the like.

さらに、タイヤがタイヤ回転軸の周りに回転することを再現するために、第1のモデル空間領域30の第1のメッシュは、タイヤ回転軸の周りに回転するメッシュとすることできる。第1のモデル空間領域30の第1のメッシュを回転可能にすることにより、シミュレーション計算結果をより正確にすることができる。この場合、回転する第1のメッシュと第2のモデル空間領域32の非回転のメッシュとの間では、上述した物理量の補間関係を設定して、物理量の受け渡しを行う。すなわち、シミュレーション計算中、第1のメッシュは回転するが、時間ステップを刻みながら時々刻々回転するので、時間ステップ毎に、図3(a),(b)に示す方法を用いて、物理量の受け渡すべきメッシュ領域を探索し、探索したメッシュ領域との間で、上述した物理量の受け渡しを行う。   Furthermore, in order to reproduce the rotation of the tire around the tire rotation axis, the first mesh in the first model space region 30 can be a mesh that rotates around the tire rotation axis. By making the first mesh in the first model space region 30 rotatable, the simulation calculation result can be made more accurate. In this case, between the rotating first mesh and the non-rotating mesh in the second model space area 32, the physical quantity is transferred by setting the above-described physical quantity interpolation relationship. That is, during the simulation calculation, the first mesh rotates, but it rotates every moment while ticking the time step. Therefore, the physical quantity is received using the method shown in FIGS. 3A and 3B for each time step. The mesh area to be transferred is searched, and the above-described physical quantity is transferred to and from the searched mesh area.

上述したタイヤを取り巻く空間中の流体の流れをシミュレートするためのシミュレーションモデルの作成方法は、ROM16に記憶されたプログラムを呼び出して起動することにより実行される。このようなプログラムは、
少なくともタイヤの一部の形状を境界とする第1の空間領域を再現した第1のモデル空間領域30と、第1の空間領域の周りに位置する第2の空間領域を再現した第2のモデル空間領域32を別々にメッシュ分割する処理をコンピュータに行わせる手順と、
第1のモデル空間領域30において作成された第1のメッシュと、2のモデル空間領域32において作成された第2のメッシュとの間の流体の物理量を補間する関係をコンピュータに定めさせる手順と、含む。
このようなプログラムは、電気通信回線を通じて入手することができ、あるいは、ROMあるいは記録媒体に記憶された形態で入手することができる。
The method for creating a simulation model for simulating the flow of fluid in the space surrounding the tire described above is executed by calling and starting a program stored in the ROM 16. Such programs are
A first model space region 30 that reproduces at least a first space region bounded by the shape of a part of the tire, and a second model that reproduces a second space region located around the first space region A procedure for causing the computer to perform a process of dividing the space region 32 into meshes separately;
A procedure for causing a computer to determine a relationship for interpolating a physical quantity of fluid between a first mesh created in the first model space region 30 and a second mesh created in the second model space region 32; Including.
Such a program can be obtained through a telecommunication line, or can be obtained in a form stored in a ROM or a recording medium.

以上のように、本実施形態では、少なくともタイヤの一部の形状を境界とする第1の空間領域を再現する第1のモデル空間領域30と、第1の空間領域の周りに位置する第2の空間領域を再現する第2のモデル空間領域32を別々にメッシュ分割し、第1のモデル空間領域30において作成された第1のメッシュと、第2のモデル空間領域32において作成された第2のメッシュとの間の流体の物理量を補間する関係を定める。このため、図4に示すように1つのシミュレーションモデルとして作成することができる。タイヤの形状を種々変更するとき、第1のモデル空間領域30を作成し直せばよいので、図4に示すようなシミュレーションモデル全体を作成し直す必要が無くなり、シミュレーションモデルの作成時間を短縮することができる。   As described above, in the present embodiment, the first model space region 30 that reproduces the first space region having at least a part of the shape of the tire as a boundary, and the second model space that is located around the first space region. The second model space region 32 that reproduces the space region is separately divided into meshes, and the first mesh created in the first model space region 30 and the second mesh created in the second model space region 32 are divided. The relationship for interpolating the physical quantity of the fluid between the meshes is defined. Therefore, it can be created as one simulation model as shown in FIG. When the tire shape is changed variously, the first model space region 30 only needs to be recreated, so that it is not necessary to recreate the entire simulation model as shown in FIG. 4, and the simulation model creation time is shortened. Can do.

上述の実施形態は、タイヤの周りに流体を満たした第1及び第2の空間領域を再現したモデル空間領域(第1のモデル空間領域、第2のモデル空間領域)を含んだシミュレーションモデルを作成したが、タイヤを装着した車両の周りに流体を流す空間を再現したモデル空間領域をシミュレーションモデルとして作成することもできる。図6は、本実施形態の作成方法で用いる空間領域と車両モデルの一例を示す図である。タイヤを含む車両の形状を再現した剛体モデルである車両モデル40を、モデル空間領域42内に配して、車両モデル40に装着されたタイヤモデルの形状を変更する場合、タイヤモデルに比べてモデル空間領域42は極めて大きい。このため、タイヤの形状を変更する度にモデル空間領域42のメッシュ分割を繰り返すことは多大な時間を要する。このため、本実施形態で説明したように、少なくともタイヤの一部の形状を境界とする第1の空間領域を再現した第1のモデル空間領域30をメッシュ分割し、メッシュ分割した第1のモデル空間領域30を、第1のモデル空間領域30の周りを囲む、既に作成されている、メッシュ分割した第2のモデル空間領域32に組み込むことにより、短時間にシミュレーションモデルを作成することができる。なお、モデル空間領域42は、第1のモデル空間領域30と第2のモデル空間領域32を含む領域である。   The above-described embodiment creates a simulation model including a model space region (first model space region, second model space region) that reproduces the first and second space regions filled with fluid around the tire. However, it is also possible to create a model space region that reproduces a space in which a fluid flows around a vehicle equipped with tires as a simulation model. FIG. 6 is a diagram illustrating an example of a space area and a vehicle model used in the creation method of the present embodiment. When the vehicle model 40, which is a rigid model that reproduces the shape of the vehicle including the tire, is arranged in the model space region 42 and the shape of the tire model mounted on the vehicle model 40 is changed, the model is compared with the tire model. The space area 42 is extremely large. For this reason, it takes a lot of time to repeat the mesh division of the model space region 42 every time the shape of the tire is changed. Therefore, as described in the present embodiment, the first model space region 30 that reproduces the first space region having at least a part of the shape of the tire as a boundary is mesh-divided, and the mesh-divided first model A simulation model can be created in a short time by incorporating the space region 30 into the already created mesh-divided second model space region 32 surrounding the first model space region 30. The model space area 42 is an area including the first model space area 30 and the second model space area 32.

(実験例)
図6に示すモデル空間領域42のサイズを、縦26m×横26m×高さ11mとし、この領域の中心部分に、実車相当サイズの車両モデル40を配し、モデル空間領域の一方の端から空気を流して車両モデル42周りの空気の流れをシミュレートするためのシミュレーションモデルの作成を行った。シミュレーションモデルのメッシュの格子点の総数は約600万個である。このとき、従来のようにモデル空間領域42全体を一度にメッシュ分割してシミュレーションモデルを作成する方法(従来の方法)と、図2に示すフローに沿って、第1のモデル空間領域30と第2のモデル空間領域32を別々に作成する本実施形態の方法と、を行った。
本実施形態のステップS40までに要するシミュレーションモデルの作成時間は、従来の方法の作成時間と同等であったが、ステップS50〜S80の間のシミュレーションモデルの作成時間は、従来の方法の作成時間対比5分の1であった。すなわち、本実施形態の方法は、第2のモデル空間領域32のメッシュ分割をしないことにより作成時間を短縮することができる。これより、本実施形態の効果は明らかである。
(Experimental example)
The size of the model space area 42 shown in FIG. 6 is 26 m long × 26 m wide × 11 m high. A vehicle model 40 having a size equivalent to an actual vehicle is arranged at the center of this area, and air is passed from one end of the model space area A simulation model for simulating the air flow around the vehicle model 42 was created. The total number of mesh points of the mesh of the simulation model is about 6 million. At this time, according to the conventional method of creating a simulation model by dividing the entire model space region 42 at a time (conventional method) and the first model space region 30 and the first model space according to the flow shown in FIG. The method of this embodiment in which the two model space regions 32 are separately created was performed.
The simulation model creation time required until step S40 of the present embodiment is equivalent to the creation time of the conventional method, but the simulation model creation time between steps S50 to S80 is compared with the creation time of the conventional method. It was 1/5. That is, the method of the present embodiment can reduce the creation time by not dividing the second model space region 32 into meshes. From this, the effect of this embodiment is clear.

以上、本発明のシミュレーションモデルの作成方法及びプログラムについて詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良や変更をしてもよいのはもちろんである。   Although the simulation model creation method and program of the present invention have been described in detail above, the present invention is not limited to the above-described embodiment, and various improvements and modifications may be made without departing from the spirit of the present invention. Of course.

10 モデル作成装置
12 CPU
14 メモリ
16 ROM
18 入出力部
20 領域設定部
22 メッシュ分割部
24 補間関係設定部
26 制御部
27 入力操作デバイス
30 第1のモデル空間領域
32 第2のモデル空間領域
34,34a,34b 第1のメッシュ領域
36,36a,36b 第2のメッシュ領域
38 凹部
40 車両モデル
42 モデル空間領域
10 Model creation device 12 CPU
14 Memory 16 ROM
18 Input / output unit 20 Region setting unit 22 Mesh division unit 24 Interpolation relation setting unit 26 Control unit 27 Input operation device 30 First model space region 32 Second model space region 34, 34a, 34b First mesh region 36, 36a, 36b Second mesh area 38 Recess 40 Vehicle model 42 Model space area

Claims (7)

タイヤを取り巻く空間中の流体の流れをシミュレートするためのシミュレーションモデルの作成方法であって、
(1)少なくともタイヤの一部の形状を境界とする第1の空間領域を再現した第1のモデル空間領域と、前記第1の空間領域の周りに位置する第2の空間領域を再現した第2のモデル空間領域を別々にメッシュ分割するステップと、
(2)前記第1のモデル空間領域において作成された第1のメッシュと、前記2のモデル空間領域において作成された第2のメッシュとの間の前記流体の物理量を補間する関係を定めるステップと、含み、
前記第1のモデル空間領域と前記第2のモデル空間領域とは、互いに離間しており、前記第1のモデル空間領域と前記第2のモデル空間領域との離間距離の寸法は、前記第1のメッシュの最小メッシュサイズの寸法よりも小さく、
前記流体の物理量を補間する関係を定める前記ステップにおいて、前記第1のメッシュが前記第2のメッシュに、及び前記第2のメッシュが前記第1のメッシュに物理量を与えるとき、前記第1のメッシュ及び前記第2のメッシュの中で、物理量を受ける側のメッシュAに対向する、前記物理量を与える側のメッシュBが複数ある場合、前記メッシュAに対向する前記メッシュBの境界面の部分の大きさの比率に応じて前記メッシュBの物理量を加重平均した値を前記メッシュAに与えるように、前記関係を定める、ことを特徴とするシミュレーションモデルの作成方法。
A method of creating a simulation model for simulating fluid flow in a space surrounding a tire,
(1) A first model space region that reproduces a first space region that borders at least a part of the shape of the tire, and a second space region that is located around the first space region. Meshing the two model space regions separately;
(2) defining a relationship for interpolating the physical quantity of the fluid between the first mesh created in the first model space region and the second mesh created in the second model space region; It includes,
The first model space region and the second model space region are separated from each other, and the dimension of the separation distance between the first model space region and the second model space region is the first model space region. Smaller than the minimum mesh size dimension of the mesh,
In the step of determining a relationship for interpolating the physical quantity of the fluid, the first mesh is applied when the first mesh gives the second mesh and the second mesh gives the first mesh a physical quantity. In the second mesh, when there are a plurality of physical quantity-giving meshes B facing the physical quantity receiving mesh A, the size of the boundary surface portion of the mesh B facing the mesh A A method for creating a simulation model , wherein the relationship is determined so that a value obtained by weighted averaging of physical quantities of the mesh B is given to the mesh A in accordance with the ratio .
前記第1のメッシュの平均サイズは、前記第2のメッシュの平均サイズに比べて小さい、請求項1に記載のシミュレーションモデルの作成方法。   The simulation model creation method according to claim 1, wherein an average size of the first mesh is smaller than an average size of the second mesh. 前記第1のメッシュは、前記タイヤをモデル化したタイヤモデルと接し、
前記タイヤモデルのタイヤ幅方向における、前記タイヤモデルの幅寸法をSt[mm]とし、前記タイヤモデルの直径寸法をDt[mm]としたとき、前記第1のモデル空間領域の前記タイヤ幅方向の幅寸法S[mm]は、前記幅寸法Stの0.1倍より大きく3倍より小さい寸法であり、かつ、前記第1の空間領域の前記タイヤ径方向に沿った寸法の最大値は、前記直径寸法Dtの0.3倍より大きく2倍より小さい寸法である、請求項1または2に記載のシミュレーションモデルの作成方法。
The first mesh is in contact with a tire model that models the tire;
When the width dimension of the tire model in the tire width direction of the tire model is St [mm] and the diameter dimension of the tire model is Dt [mm], the tire width direction of the first model space region is the tire width direction. The width dimension S [mm] is a dimension larger than 0.1 times and smaller than three times the width dimension St, and the maximum value of the dimension along the tire radial direction of the first space region is The method for creating a simulation model according to claim 1 or 2, wherein the size is larger than 0.3 times and smaller than twice the diameter dimension Dt.
前記第1のモデル空間領域は、前記タイヤのサイドウォールと接する空間を再現したモデル空間領域を含む、請求項1〜のいずれか1項に記載のシミュレーションモデルの作成方法。 The first model space region comprising said model space area that reproduces the sidewall in contact with the space of the tire, how to create a simulation model according to any one of claims 1-3. 前記第1のモデル空間領域は、前記タイヤのトレッドパターンと接する空間を再現したモデル空間領域を含む、請求項1〜のいずれか1項に記載のシミュレーションモデルの作成方法。 The simulation model creation method according to any one of claims 1 to 4 , wherein the first model space region includes a model space region that reproduces a space in contact with a tread pattern of the tire. 前記タイヤがタイヤ回転軸の周りに回転することを再現するために、前記第1のメッシュは、タイヤ回転軸の周りに回転するメッシュである、請求項1〜のいずれか1項に記載のシミュレーションモデルの作成方法。 For the tire to reproduce that rotates about the tire rotation axis, the first mesh is a mesh which rotates around the tire rotational axis, according to any one of claims 1 to 5 How to create a simulation model. タイヤを取り巻く空間中の流体の流れをシミュレートするためのシミュレーションモデルの作成方法をコンピュータに実行させるコンピュータが実行可能なプログラムであって、
(3)少なくともタイヤの一部の形状を境界とする第1の空間領域を再現した第1のモデル空間領域と、前記第1の空間領域の周りに位置する第2の空間領域を再現した第2のモデル空間領域を別々にメッシュ分割する処理をコンピュータに行わせる手順と、
(4)前記第1のモデル空間領域において作成された第1のメッシュと、前記2のモデル空間領域において作成された第2のメッシュとの間の前記流体の物理量を補間する関係をコンピュータに定めさせる手順と、含み
前記第1のモデル空間領域と前記第2のモデル空間領域とは、互いに離間しており、前記第1のモデル空間領域と前記第2のモデル空間領域との離間距離の寸法は、前記第1のメッシュの最小メッシュサイズの寸法よりも小さく、
前記流体の物理量を補間する関係をコンピュータに定めさせる前記手順において、前記第1のメッシュが前記第2のメッシュに、及び前記第2のメッシュが前記第1のメッシュに物理量を与えるとき、前記第1のメッシュ及び前記第2のメッシュの中で、物理量を受ける側のメッシュAに対向する、前記物理量を与える側のメッシュBが複数ある場合、前記メッシュAに対向する前記メッシュBの境界面の部分の大きさの比率に応じて前記メッシュBの物理量を加重平均した値を前記メッシュAに与えるように、前記関係を定めさせる、ことを特徴とするプログラム。
A computer-executable program for causing a computer to execute a method of creating a simulation model for simulating a fluid flow in a space surrounding a tire,
(3) A first model space region that reproduces a first space region that borders at least a part of the shape of the tire, and a second space region that is located around the first space region. A procedure for causing a computer to perform the process of dividing the model space area of the two meshes separately;
(4) A relationship for interpolating the physical quantity of the fluid between the first mesh created in the first model space region and the second mesh created in the second model space region is determined in the computer. includes a procedure for, the,
The first model space region and the second model space region are separated from each other, and the dimension of the separation distance between the first model space region and the second model space region is the first model space region. Smaller than the minimum mesh size dimension of the mesh,
In the step of causing the computer to determine a relationship for interpolating the physical quantity of the fluid, when the first mesh gives a physical quantity to the second mesh and the second mesh gives the physical quantity to the first mesh, the first mesh In the case where there are a plurality of meshes B on the physical quantity-given side facing the mesh A on the physical quantity receiving side among the one mesh and the second mesh, the boundary surface of the mesh B facing the mesh A A program characterized in that the relationship is determined so that a value obtained by weighted averaging of physical quantities of the mesh B is given to the mesh A according to a ratio of the sizes of the parts .
JP2013254955A 2013-12-10 2013-12-10 Method and program for creating simulation model Active JP6281270B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013254955A JP6281270B2 (en) 2013-12-10 2013-12-10 Method and program for creating simulation model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013254955A JP6281270B2 (en) 2013-12-10 2013-12-10 Method and program for creating simulation model

Publications (2)

Publication Number Publication Date
JP2015114151A JP2015114151A (en) 2015-06-22
JP6281270B2 true JP6281270B2 (en) 2018-02-21

Family

ID=53528083

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013254955A Active JP6281270B2 (en) 2013-12-10 2013-12-10 Method and program for creating simulation model

Country Status (1)

Country Link
JP (1) JP6281270B2 (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6737656B2 (en) * 2016-08-01 2020-08-12 Toyo Tire株式会社 Wind tunnel test method
JP6737657B2 (en) * 2016-08-01 2020-08-12 Toyo Tire株式会社 Simulation method and program, and simulation device
JP6861583B2 (en) * 2017-06-19 2021-04-21 Toyo Tire株式会社 Fluid analysis method, fluid analysis device, and program around the rotating body
CN109271651B (en) * 2018-07-09 2023-01-20 广东工业大学 Laser-induction-based lattice Boltzmann gas-liquid two-phase flow simulation method
CN109946096B (en) * 2018-11-26 2020-09-25 西安理工大学 High-speed train air pipe fault diagnosis method based on model space

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4466102B2 (en) * 2004-02-05 2010-05-26 横浜ゴム株式会社 Method for creating structure simulation model and structure simulation method
JP4792049B2 (en) * 2008-01-09 2011-10-12 住友ゴム工業株式会社 Tire noise performance simulation method and tire manufacturing method
JP5672797B2 (en) * 2010-06-25 2015-02-18 横浜ゴム株式会社 Simulation model creation method, simulation method, simulation model creation device, and simulation device

Also Published As

Publication number Publication date
JP2015114151A (en) 2015-06-22

Similar Documents

Publication Publication Date Title
JP6281270B2 (en) Method and program for creating simulation model
US9165093B2 (en) Simulation method for tire
JP5064669B2 (en) Method for simulating behavior of tire / wheel assembly and simulating behavior of tire
JP5329180B2 (en) Method for analyzing mechanical behavior of golf ball and golf ball
JP2008308801A (en) Method, apparatus, and program for preparing model for analyzing tire cord, and method for analyzing tire performance
JP5533185B2 (en) Tire peripheral space analysis method, tire peripheral space analysis computer program, and analysis apparatus
JP2013088871A (en) Simulation method, computer program for simulation, and simulation device
JP2002358473A (en) Gas flow simulation method
JP2006018422A (en) Tire finite element modelling method
JP2007223504A (en) Method and device for estimating tire performance, and storage medium
JP2007230403A (en) System for producing tire model, method for producing tire model, tire model, and method for simulating behavior of tire
JP2006072893A (en) Tire model generation method, and simulation method using the tire model
JP5811625B2 (en) Simulation method and simulation apparatus
JP5782684B2 (en) Simulation model creation method, simulation method, simulation model creation device, and simulation device
JP2011201390A (en) Tire performance prediction method and tire performance prediction device
JP2017067486A (en) Method for creating tire model and method for simulating tire temperature
JP2007045363A (en) Method, device and program for preparing tire model
US20200159979A1 (en) Method and device for rotation analysis of rotating body
JP2003240651A (en) Physical quantity display method for tire
JP6424543B2 (en) Tire simulation method and tire performance evaluation method
JP2012181600A (en) Tire model creation method, tire model creation device, tire model creation program, and tire performance analysis method
JP4615983B2 (en) Tire model creation method, tire model, and tire behavior simulation method
JP7439393B2 (en) Tire simulation method
JP2012056393A (en) Tire model creation method, tire model creation apparatus, and tire model creation program
JP2012192888A (en) Tire performance simulation method, tire performance simulation system, and tire performance simulation program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20161202

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170815

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170810

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20171012

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180108

R150 Certificate of patent or registration of utility model

Ref document number: 6281270

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

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