JP4548005B2 - 不均質材料のシミュレーションモデル作成用処理装置の動作方法 - Google Patents
不均質材料のシミュレーションモデル作成用処理装置の動作方法 Download PDFInfo
- Publication number
- JP4548005B2 JP4548005B2 JP2004173619A JP2004173619A JP4548005B2 JP 4548005 B2 JP4548005 B2 JP 4548005B2 JP 2004173619 A JP2004173619 A JP 2004173619A JP 2004173619 A JP2004173619 A JP 2004173619A JP 4548005 B2 JP4548005 B2 JP 4548005B2
- Authority
- JP
- Japan
- Prior art keywords
- boundary
- heterogeneous material
- model
- shape
- finite element
- 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
Images
Description
そして、この有限要素モデルがあたかも上下、左右、奥行き方向に連続して無限に連なってバルクを成すように、有限要素モデルの対向する面に周期境界条件を付与してシミュレーションモデルを作成してシミュレーション演算を行う。
しかし、自動メッシュ分割では各有限要素の形状が同一でないため、代表領域の対向する2つの境界面上に位置する各有限要素の面の形状も異なり、したがって境界面上に位置する有限要素の節点の位置が互いに異なる。このため、自動メッシュ分割されて生成された代表領域の境界面上の節点の位置は、マニュアル修正されなければならない。
しかし、このマニュアル修正では、例えば縦、横、奥行き方向にそれぞれ256個有限要素を配列して有限要素モデルを構成した場合、修正すべき節点の数は極めて多数になる。このため、マニュアル修正に要する時間は極めて長くなり、また、その修正処理も煩雑である。
また、前記境界離散点の投影先が、前記投影面を形成する複数の有限要素の境界上に位置する場合、前記複数の有限要素のうちいずれか1つの有限要素の離散点を用いて前記境界離散点の挙動を規制するのが好ましい。
処理装置10は、入力操作系12、コンピュータ14及びディスプレイ16を有する。
入力操作系12は、マウスやキーボードであり、各種情報をオペレータの指示により入力するデバイスである。
ディスプレイ16は、入力操作系12を用いてオペレータが指示できるように入力画面を表示し、又後述する有限要素モデル及びシミュレーション演算結果を表示する部分である。
本実施形態で作成される有限要素モデルは3次元モデルであるが、2次元モデルであってもよい。
具体的には、後述するように拘束される境界面上に位置する節点の自由度を、他方の境界面上に位置する拘束する節点の自由度によって規制する。このときの拘束される節点の自由度は、拘束する各節点の自由度に所定の重み係数を乗算して加算した関係式によって表される。この関係式は、有限要素モデル作成部22にて作成された全体行列に適用されて、全体行列から拘束される節点の自由度が消去されて周期境界条件が付与されたシミュレーション可能な縮小化行列が作成される。
上記関係式及び縮小化行列の作成については後述する。
図2は不均質材料のシミュレーションモデルの作成方法の一例の流れを示すフローチャートである。図3は不均質材料の代表領域の一例の断面図である。この不均質材料は、第1のポリマー相A(図3中白色領域)と、第2のポリマー相B(図3中黒領域)と、粒系フィラー相C(図3中の円形状の灰色領域)と、粒系フィラー相の周りを取り巻くフィラー・ポリマー境界相D(図3中の灰色領域)とが不均質に分散配置されている。すなわち、材料特性の異なる複数のエラストマーが材料相として分散配置され、さらに、粒状の補強材が材料相として分散配置されている。この不均質材料は、複数の材料相における弾性率のうち、最大の弾性率は最小の弾性率の100倍以上である。不均質材料を構成する材料相は、エラストマーやフィラー等の固体相であるばかりでなく、気体や液体が満たされた空隙相であってもよい。
このような不均質材料について以下説明するようにシミュレーションモデルが作成される。
代表領域は直方体形状であり、互いに対向する3対の境界面を有する。各対の境界面を、同一の代表領域の対応する境界面と隣接して無限に配置すると、代表領域が周期的に接続されたバルク状の不均質材料が構成されるように代表領域は設定されている。
このような代表領域には異なる材料相が分散配置されて構成され、材料相の形状に従って、公知の自動メッシュ分割アルゴリズムによりメッシュ分割されて有限要素モデルが作成される。有限要素モデルは、設定された代表領域について複数の節点によって形状が特徴付けられた複数の有限要素を配置して1つの代表領域を再現した離散化モデルである。
境界節点は、例えば図4では代表領域の隣接配置方向における境界面FAを投影面として、境界面FB上に位置する節点を投影したとき、境界面FA上の投影面(境界面FA)上に投影されて、境界面FB上に位置する拘束すべき節点をいう。
このような拘束すべき節点が検出される。
図6は、境界節点についてわかり易く説明した図である。
図6では、投影面となる境界面FAが境界線として表され、境界面FBが境界線として表され、拘束されるべき節点50の境界面FAへの投影が示されている。その際、節点50の投影面上の投影先が、境界面FA上のどの有限要素の境界面上に位置するか、検出し、そのときの境界面を形成する節点の位置情報が同時に求められる。
ここで、境界面FB上に位置する節点50を境界面FAに投影したとき、図7に示すように節点50が節点51,52,54,55によって形成される有限要素の境界面上に投影されたとする。この場合、ステップS20において、節点51,52,54,55の位置情報が求められる。そして、この節点51,52,54,55の位置情報と、節点50の投影先の位置情報とを用いて、以下に示す方法により重み係数が求められる。
節点51,52,54,55で形成される有限要素の面上の物理空間上の形状を、所定のパラメトリック空間上の単位セル形状、例えば矩形形状から形状関数を用いて形状変換を施したものとして定めることにより、境界節点の単位セル形状内における対応点の位置情報を求め、この位置情報および形状関数を用いて周期境界条件に用いる重み係数を求める。
したがって、有限要素モデルの境界面に位置する有限要素毎にこの有限要素の面に投影される、あるいは有限要素の辺上に投影される節点に対して、変換Tによりパラメトリック空間の正方形上における対応点60を求めることができる。
なお、上記パラメトリック空間上の基準形状として正方形を用いているが、3角形の形状であってもよく、本発明において基準形状は正方形に限定されない。また、本発明では、3次元形状で表されるモデルの境界面の形状を変換Tにて変換するものであるが、2次元形状で表されるモデルの場合前記境界面は境界線となり、この場合パラメトリック空間上の基準形状は線分となる。したがって、この場合、変換Tによりモデルの境界線はパラメトリック空間上の基準形状である線分へ変換される。
このようにして、重み係数Wが算出される。
節点51,52,55,54の所定の物理量、例えばX方向の変位をu(51),u(52),u(55),u(54)とすると、節点50におけるX方向の変位u(50)は下記式(2)のように定める。
こうして定められた関係式が周期境界条件として有限要素モデルへ付与される。(ステップS40)。
例えば、ステップS10において作成された有限要素モデルの全体行列から、上記関係式を用いて拘束される境界節点の自由度を消去した縮小化行列を作成する。この縮小化行列の作成は不均質材料の挙動を静的に再現するシミュレーション演算に好適であり、境界節点の自由度を消去して行列の自由度を少なくするので解析に要する演算時間を短縮することができる。例えば、下記式(3)のような有限要素の行列において、u3=u2の関係式が定められた場合、下記式(4)のように、u3のマトリクス成分を削除して、u1,u2,u4のマトリクス成分で表す。
次に、この縮小化行列を用いてシミュレーション演算が成される(ステップS50)。
図9(a),(b)は、縮小化剛性行列を作成してシミュレーション演算を行った結果を示す図である。図9(a)は、所定の外力が付与されたときの有限要素モデルの変形状態を示した変形図であり、図9(b)はその時の最大主歪み分布である。
図9(a)中の白色、灰色、黒色の領域は、それぞれ図3における第2のポリマー相B、第1のポリマー相A及びフイラー・ポリマー境界相D、粒系フィラー相Cを示している。有限要素モデルの境界面は外力により変形している。
例えば、力学挙動を再現するシミュレーション演算において、有限要素モデルの境界面上に位置する節点のうち、拘束される境界節点に作用する力とこの境界節点の質量を、拘束する節点に、重み係数Wに応じた分配を行なって、拘束する節点の質量および作用する力を修正する。この場合、拘束される境界節点すべてに対して上記分配を行なった後、拘束される境界節点を除いた有限要素モデルの節点の加速度を計算する。この後、拘束する節点の加速度は、重み係数Wを用いて拘束される境界節点に振り分けられる。こうして、拘束される境界節点の加速度を求めることができる。求められた境界節点における加速度とこの境界節点における質量とを用いてこの境界節点に加わる力を求める。こうして、次の時間ステップにおいて作用する力を求める。勿論外力として次の時間ステップにおいて境界節点に力が付加された場合、この力も加算される。このような陽解法は、不均質材料の動解析のシミュレーション演算に好適である。
以上が、不均質材料のシミュレーションモデルの作成方法の説明である。
メッシュフリー法とは、図10に示すように、複数の要素点(離散点)が連続して配置される状態において、注目する要素点50a(図10中、黒丸)を中心として距離ρ0 内の範囲をサポートとし、この範囲において図10に示すようなスプライン関数、例えば4次のスプライン関数で重み関数w(x,y)が設定されている。すなわち、各要素点において距離ρ0内の範囲をサポートとして単位セルが設定されている。
12 入力操作系
14 コンピュータ
16 ディスプレイ
18 CPU
20 メモリ
22 有限要素モデル作成部
24 境界条件設定部
26 シミュレーション演算部
Claims (8)
- CPUと、
記憶手段と、
ディスプレイと、
を備えた、材料特性の異なる複数の材料相が分散配置された不均質材料のシミュレーションモデル作成用処理装置の動作方法であって、
前記CPUが、
不均質材料における代表領域を少なくとも2方向に連続的に隣接配置することによって不均質材料の塊を構成する際の前記代表領域を定め、この代表領域について、複数の離散点によって形状が特徴付けられた複数の有限要素を用いて1つの代表領域を再現した離散化モデルを作成するステップと、
前記代表領域の隣接配置方向における前記離散化モデルの一方の境界線又は境界面を投影線又は投影面として、他方の境界線又は境界面に位置する複数の離散点を投影したとき、前記投影線又は投影面上に投影される前記他方の境界線又は境界面に位置する複数の離散点を境界離散点としてそれぞれ検出するステップと、
前記境界離散点がそれぞれ投影される各有限要素の前記投影線又は投影面の形状を、所定のパラメトリック空間上の単位セル形状から形状関数を用いて形状変換したものとして定めることにより、前記境界離散点の前記単位セル形状内における対応点の位置情報を求め、この位置情報および前記形状関数を用いて前記境界離散点のそれぞれの挙動を各離散点が投影された有限要素の前記投影線又は投影面を形成する離散点の挙動によって拘束し、前記代表領域が連続して連なるように前記離散化モデルの対向する境界線又は境界面に付与する周期境界条件を設定すると共に、前記離散化モデルのそれぞれの対向する境界線又は境界面について周期境界条件を順次前記記憶手段に記憶保持するステップと、
この周期境界条件を前記離散化モデルに付与してシミュレーションモデルを作成するステップと、
作成されたシミュレーションモデルを前記ディスプレイにより出力するステップと
を実行することを特徴とする不均質材料のシミュレーションモデル作成用処理装置の動作方法。 - 前記周期境界条件は、前記有限要素の前記投影面の形状を規定する離散点の挙動に重み付け係数を乗算して前記境界離散点の挙動を規制した多項式で表されたものである請求項1に記載の不均質材料のシミュレーションモデル作成用処理装置の動作方法。
- 前記境界離散点の投影先が、前記投影面を形成する複数の有限要素の境界上に位置する場合、前記複数の有限要素のうちいずれか1つの有限要素の離散点を用いて前記境界離散点の挙動を規制する請求項2に記載の不均質材料のシミュレーションモデル作成用処理装置の動作方法。
- 前記代表領域は、複数の境界線又は境界面が角度を成して接続された角部を有する多角形または多面体形状を有し、
前記離散化モデルは、前記角部に対応する位置に前記有限要素の離散点を有し、この離散点の挙動は固定される請求項1〜3のいずれか1項に記載の不均質材料のシミュレーションモデル作成用処理装置の動作方法。 - 前記シミュレーションモデルは、前記不均質材料の歪み・応力解析のシミュレーションに用いられる請求項1〜4のいずれか1項に記載の不均質材料のシミュレーションモデル作成用処理装置の動作方法。
- 前記離散化モデルの対象となる不均質材料の前記複数の材料相における弾性率のうち、最大の弾性率は、最小の弾性率の100倍以上である請求項5に記載の不均質材料のシミュレーションモデル作成用処理装置の動作方法。
- 前記不均質材料には、粒状の補強材が材料相として分散配置されている請求項1〜6のいずれか1項に記載の不均質材料のシミュレーションモデル作成用処理装置の動作方法。
- 前記不均質材料には、材料特性の異なる複数のエラストマーが材料相として分散配置されている請求項1〜7のいずれか1項に記載の不均質材料のシミュレーションモデル作成用処理装置の動作方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004173619A JP4548005B2 (ja) | 2004-06-11 | 2004-06-11 | 不均質材料のシミュレーションモデル作成用処理装置の動作方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004173619A JP4548005B2 (ja) | 2004-06-11 | 2004-06-11 | 不均質材料のシミュレーションモデル作成用処理装置の動作方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2005352819A JP2005352819A (ja) | 2005-12-22 |
JP4548005B2 true JP4548005B2 (ja) | 2010-09-22 |
Family
ID=35587248
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004173619A Expired - Fee Related JP4548005B2 (ja) | 2004-06-11 | 2004-06-11 | 不均質材料のシミュレーションモデル作成用処理装置の動作方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4548005B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103017167A (zh) * | 2012-11-21 | 2013-04-03 | 广东电网公司电力科学研究院 | 多孔介质燃烧器的制造方法 |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4931465B2 (ja) * | 2006-03-31 | 2012-05-16 | 株式会社日本総合研究所 | メッシュ生成装置、メッシュ生成方法、及びコンピュータプログラム |
JP4931511B2 (ja) * | 2006-08-08 | 2012-05-16 | 中国電力株式会社 | 応力緩和を考慮した逆解析による高温部材のクリープ余寿命の評価方法 |
JP4685747B2 (ja) * | 2006-11-09 | 2011-05-18 | 住友ゴム工業株式会社 | ゴム材料解析モデルの作成方法 |
JP5169279B2 (ja) * | 2008-02-14 | 2013-03-27 | 横浜ゴム株式会社 | 不均質材料モデルの作成装置の動作方法、この方法を用いた不均質材料のシミュレーション装置の動作方法、不均質材料モデルの作成装置及び不均質材料のシミュレーション装置 |
JP5169278B2 (ja) * | 2008-02-14 | 2013-03-27 | 横浜ゴム株式会社 | 複合材料モデルの作成装置の動作方法、複合材料の変形挙動のシミュレーション装置の動作方法及び複合材料シミュレーション装置 |
JP5180735B2 (ja) * | 2008-08-20 | 2013-04-10 | 住友ゴム工業株式会社 | ゴム材料のシミュレーション方法 |
EP2189918A1 (en) * | 2008-11-07 | 2010-05-26 | Dassault Systèmes | Computer-implemented method of computing, in a computer aided design system, of a boundary of a modeled object. |
JP5771935B2 (ja) * | 2010-10-04 | 2015-09-02 | 横浜ゴム株式会社 | シミュレーションモデルの作成方法及びそのコンピュータプログラム、シミュレーション方法及びそのコンピュータプログラム、並びにシミュレーションモデルの作成装置 |
JP5503618B2 (ja) | 2011-10-03 | 2014-05-28 | 住友ゴム工業株式会社 | ゴム材料のシミュレーション方法 |
JP5854067B2 (ja) * | 2014-02-28 | 2016-02-09 | 横浜ゴム株式会社 | 不均質材料のシミュレーションモデルの作成方法、不均質材料のシミュレーション方法、及びプログラム |
CN105631078B (zh) * | 2014-11-07 | 2019-01-18 | 中国石油化工股份有限公司 | 天然裂缝性油藏自适应介质的数值模拟方法 |
CN107357946B (zh) * | 2016-05-10 | 2022-10-21 | 无锡时代天使医疗器械科技有限公司 | 基于热压膜成型技术的壳状牙科器械制作工艺的验证方法 |
CN114818197B (zh) * | 2022-05-10 | 2024-04-12 | 西安交通大学 | 基于边界元模型的高速电主轴热弹性变形模拟方法及系统 |
JP2024039497A (ja) * | 2022-09-09 | 2024-03-22 | 株式会社ブリヂストン | タイヤモデル作成方法、タイヤモデル作成装置、及びタイヤモデル作成プログラム |
CN115310336B (zh) * | 2022-09-29 | 2023-01-17 | 深圳市威兆半导体股份有限公司 | 一种基于igbt高温特性的有限元建模方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11328157A (ja) * | 1998-05-13 | 1999-11-30 | Mitsubishi Electric Corp | データマッピング装置及びデータマッピング方法 |
JP2003347301A (ja) * | 2002-05-29 | 2003-12-05 | Fujitsu Ltd | 均質化法を用いた構造解析方法及び構造設計方法 |
-
2004
- 2004-06-11 JP JP2004173619A patent/JP4548005B2/ja not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11328157A (ja) * | 1998-05-13 | 1999-11-30 | Mitsubishi Electric Corp | データマッピング装置及びデータマッピング方法 |
JP2003347301A (ja) * | 2002-05-29 | 2003-12-05 | Fujitsu Ltd | 均質化法を用いた構造解析方法及び構造設計方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103017167A (zh) * | 2012-11-21 | 2013-04-03 | 广东电网公司电力科学研究院 | 多孔介质燃烧器的制造方法 |
CN103017167B (zh) * | 2012-11-21 | 2016-06-08 | 广东电网公司电力科学研究院 | 多孔介质燃烧器的制造方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2005352819A (ja) | 2005-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4548005B2 (ja) | 不均質材料のシミュレーションモデル作成用処理装置の動作方法 | |
CN102682152B (zh) | 使用扩展的无网格有限元法对结构性能的数值仿真 | |
JP5169279B2 (ja) | 不均質材料モデルの作成装置の動作方法、この方法を用いた不均質材料のシミュレーション装置の動作方法、不均質材料モデルの作成装置及び不均質材料のシミュレーション装置 | |
JP5003303B2 (ja) | 複合材料のシミュレーションモデルの作成方法およびシミュレーション方法 | |
JP5854067B2 (ja) | 不均質材料のシミュレーションモデルの作成方法、不均質材料のシミュレーション方法、及びプログラム | |
Jansson et al. | A discrete mechanics model for deformable bodies | |
JP4093994B2 (ja) | 不均質材料のシミュレーション方法 | |
US6259453B1 (en) | Meshing method and apparatus | |
JP5324820B2 (ja) | 解析モデルの作成方法 | |
Tschisgale et al. | A constraint-based collision model for Cosserat rods | |
JP4631319B2 (ja) | 不均質材料のシミュレーションモデル作成方法 | |
JP4914004B2 (ja) | 変形した素材に関連する機構をモデリングするためのコンピュータ化された方法、ソフトウェアプログラムおよびシステム | |
KR20160065120A (ko) | 필러 배합 고무의 유한 요소 모델의 작성 방법 | |
JP2010205253A (ja) | 複合材料の力学的材料定数を算出する方法、複合材料中の材料の体積分率を算出する方法および記録メディア | |
JP4466101B2 (ja) | 構造体のシミュレーションモデルの作成方法およびシミュレーション方法 | |
JP5774404B2 (ja) | 解析装置、その方法及びそのプログラム | |
JP4466118B2 (ja) | シミュレーション装置の動作方法 | |
JP5404516B2 (ja) | コンピュータ支援工学解析における近傍判定方法およびシステム | |
JP4487582B2 (ja) | タイヤのシミュレーションモデル作成方法 | |
JP5771935B2 (ja) | シミュレーションモデルの作成方法及びそのコンピュータプログラム、シミュレーション方法及びそのコンピュータプログラム、並びにシミュレーションモデルの作成装置 | |
JP2005078416A (ja) | 解析モデル生成方法および装置ならびにプログラムおよびその記憶媒体 | |
WO2021205967A1 (ja) | 階層的縮退行列生成装置 | |
Hellen | How to use elements effectively | |
JP3872469B2 (ja) | データ処理方法、データ処理装置 | |
Zou et al. | A geometrically exact isogeometric blended shell: formulation, benchmarking, and automotive application |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070322 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100330 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100526 |
|
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: 20100615 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100628 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130716 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130716 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130716 Year of fee payment: 3 |
|
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 |
|
LAPS | Cancellation because of no payment of annual fees |