JP2010201822A - 流体流動過程の解析装置、解析方法及び解析プログラム - Google Patents
流体流動過程の解析装置、解析方法及び解析プログラム Download PDFInfo
- Publication number
- JP2010201822A JP2010201822A JP2009050927A JP2009050927A JP2010201822A JP 2010201822 A JP2010201822 A JP 2010201822A JP 2009050927 A JP2009050927 A JP 2009050927A JP 2009050927 A JP2009050927 A JP 2009050927A JP 2010201822 A JP2010201822 A JP 2010201822A
- Authority
- JP
- Japan
- Prior art keywords
- fluid
- flow
- flow conductance
- conductance
- cavity
- 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.)
- Pending
Links
Images
Abstract
【解決手段】流体が流動するキャビティの少なくとも一部を複数の微小要素に分割してなる3次元モデルを構築するモデル構築処理部32と、微小要素における流体の流動コンダクタンスを決定する流動コンダクタンス決定処理部33と、決定された流動コンダクタンスに基づいて微小要素における流体の圧力を算出する算出処理部34とを備え、算出された流体の圧力により解析された解析結果を出力する流体流動過程の解析装置1において、流動コンダクタンス決定処理部33は、キャビティの幅の大きさと、流体に含有される充填材の代表長さとの関係式を含む関数方程式を解くことにより決定するものである。
【選択図】図1
Description
C×((∂2κ/∂x2)+(∂2κ/∂y2)+(∂2κ/∂z2))=−(1/η)
但し、
C=α×exp(L/t)+β×(L/t)+γ
で示され、
α×exp(L/t)+β×(L/t)+γ<1
となるとき、C=1となる(ηは流体の材料粘度、x、y及びzは微小要素の位置、κは流動コンダクタンス、Cは流動抵抗、Lは充填材の代表長さ、tはキャビティの幅の大きさ、α、β及びγは実験的に求められた定数である。)。
C×((∂2κ/∂x2)+(∂2κ/∂y2)+(∂2κ/∂z2))=−(1/η)
但し、
C=α×exp(L/t)+β×(L/t)+γ
で示され、
α×exp(L/t)+β×(L/t)+γ<1
となるとき、C=1となる(ηは流体の材料粘度、x、y及びzは微小要素の位置、κは流動コンダクタンス、Cは流動抵抗、Lは充填材の代表長さ、tはキャビティの幅の大きさ、α、β及びγは実験的に求められた定数である。)。
C×((∂2κ/∂x2)+(∂2κ/∂y2)+(∂2κ/∂z2))=−(1/η)
但し、
C=α×exp(L/t)+β×(L/t)+γ
で示され、
α×exp(L/t)+β×(L/t)+γ<1
となるとき、C=1となる(ηは流体の材料粘度、x、y及びzは微小要素の位置、κは流動コンダクタンス、Cは流動抵抗、Lは充填材の代表長さ、tはキャビティの幅の大きさ、α、β及びγは実験的に求められた定数である。)。
C×((∂2κ/∂x2)+(∂2κ/∂y2)+(∂2κ/∂z2))=−(1/η)
ηは、溶融樹脂の材料粘度、x、y及びzは、3次元空間座標軸、κは、流動コンダクタンスを示している。またCは、流動抵抗を示すものであり、下記式(2)で示されるものである。
C=α×exp(L/t)+β×(L/t)+γ
Lは、溶融樹脂に含有されるガラス繊維(充填材)の平均長さ、tは、射出成形に用いられる金型により形成される製品の肉厚を規定するキャビティの幅の大きさ、α、β及びγは、それぞれ代表的な射出成形品の例等で実験することにより求められた定数である。つまり、上記式(2)は、キャビティの幅の大きさtと、溶融樹脂に含有されるガラス繊維の平均長さLとの関係式を有している。また、上記式(2)の右辺が下記条件式(3)を具備する場合には、「C=1」とする。尚、本実施の形態においては、充填材としてガラス繊維を用いているが、本発明においては、成形品の剛性、強度、耐熱性、絶縁性、耐蝕性等の物性を高めることができるものであればよく、例えば、炭素繊維等の強化繊維、タルク、マイカ、炭酸カルシウム、アルミナ、クレー、珪藻土、石綿、硫酸バリウム、酸化チタン、カリオン、湿式又は乾式シリカ、コロイド状シリカ、リン酸カルシウム、ジルコニア等の無機物を用いることもできる。また本実施の形態においては、代表長さの一例として平均長さを用いているが、本発明においては、平均長さ以外に、最小長さ、最大長さ、中間長さを代表長さとして用いてもよい。特に充填材としてガラス繊維、炭素繊維等の強化繊維を用いた場合に、充填材の代表長さとしては例えば30〜3000μm程度が想定される。
α×exp(L/t)+β×(L/t)+γ<1
η×((∂2U/∂x2)+(∂2U/∂y2)+(∂2U/∂z2))=∂P/∂x
Uは、x軸方向における射出成形材料の流動速度、Pは、圧力を示す。
U=−κ×(∂P/∂x)
η=a×(剪断速度)b×exp(c×温度)
ここでa、b、cは、射出成形材料に特有の係数であり、公知な粘度測定装置により実験的に求められるものである。
((∂U/∂x)+(∂V/∂y)+(∂W/∂z))=0
U、V、Wは、(Uについては上述したが)x、y、zのそれぞれの座標軸方向の射出成形材料の流動速度である。
V=−κ(∂P/∂y),W=−κ(∂P/∂z)
ここで上記式(8)は、上記式(5)と同様にダルシー流れの式と呼ばれるものであり、多孔質中の浸透流を表す式である。すなわち3次元の座標軸x、y、z方向の流動速度(U、V、W)は、それぞれの方向の圧力勾配等に比例するものと仮定したものである。
κ((∂2P/∂x2)+(∂2P/∂y2)+(∂2P/∂z2))=0
((∂2κ/∂x2)+(∂2κ/∂y2)+(∂2κ/∂z2))=−(1/η)
10 入力装置
20 表示装置
30 制御部
31 入力処理部
32 モデル構築処理部
33 流動コンダクタンス決定処理部
34 算出処理部
35 出力処理部
40 記憶装置
50 キャビティ
Claims (6)
- 流体が流動する金型のキャビティの少なくとも一部を複数の微小要素に分割してなる3次元モデルを構築するモデル構築手段と、
前記微小要素における流体の流動コンダクタンスを決定する流動コンダクタンス決定手段と、
前記流動コンダクタンス決定手段により決定された流動コンダクタンスに基づいて前記微小要素における流体の圧力、圧力変化及び流動速度の少なくとも1つを算出する算出手段と
を備え、前記算出手段により算出された流体の圧力、圧力変化及び流動速度の少なくとも1つにより解析された解析結果を出力する流体流動過程の解析装置において、
前記流動コンダクタンス決定手段は、前記金型により製作される製品の肉厚を規定する前記キャビティの幅の大きさと、前記流体に含有される充填材の代表長さとの関係式を含む関数方程式を解くことにより決定することを特徴とする流体流動過程の解析装置。 - 前記流動コンダクタンス決定手段は、下記の関数方程式を解くことにより前記微小要素における流体の流動コンダクタンスを決定することを特徴とする請求項1に記載の流体流動過程の解析装置。
C×((∂2κ/∂x2)+(∂2κ/∂y2)+(∂2κ/∂z2))=−(1/η)
但し、
C=α×exp(L/t)+β×(L/t)+γ
で示され、
α×exp(L/t)+β×(L/t)+γ<1
となるとき、C=1となる(ηは流体の材料粘度、x、y及びzは微小要素の位置、κは流動コンダクタンス、Cは流動抵抗、Lは充填材の代表長さ、tはキャビティの幅の大きさ、α、β及びγは実験的に求められた定数である。)。 - 流体が流動する金型のキャビティの少なくとも一部を複数の微小要素に分割してなる3次元モデルを構築するモデル構築工程と、
前記微小要素における流体の流動コンダクタンスを決定する流動コンダクタンス決定工程と、
前記流動コンダクタンス決定工程により決定された流動コンダクタンスに基づいて前記微小要素における流体の圧力、圧力変化及び流動速度の少なくとも1つを算出する算出工程と
を含み、前記算出工程により算出された流体の圧力、圧力変化及び流動速度の少なくとも1つにより解析された解析結果を出力する流体流動過程の解析方法において、
前記流動コンダクタンス決定工程は、前記金型により製作される製品の肉厚を規定する前記キャビティの幅の大きさと、前記流体に含有される充填材の代表長さとの関係式を含む関数方程式を解くことにより決定することを特徴とする流体流動過程の解析方法。 - 前記流動コンダクタンス決定工程は、下記の関数方程式を解くことにより前記微小要素における流体の流動コンダクタンスを決定することを特徴とする請求項3に記載の流体流動過程の解析方法。
C×((∂2κ/∂x2)+(∂2κ/∂y2)+(∂2κ/∂z2))=−(1/η)
但し、
C=α×exp(L/t)+β×(L/t)+γ
で示され、
α×exp(L/t)+β×(L/t)+γ<1
となるとき、C=1となる(ηは流体の材料粘度、x、y及びzは微小要素の位置、κは流動コンダクタンス、Cは流動抵抗、Lは充填材の代表長さ、tはキャビティの幅の大きさ、α、β及びγは実験的に求められた定数である。)。 - コンピュータを、
流体が流動する金型のキャビティの少なくとも一部を複数の微小要素に分割してなる3次元モデルを構築するモデル構築手段と、
前記微小要素における流体の流動コンダクタンスを決定する流動コンダクタンス決定手段と、
決定された流動コンダクタンスに基づいて前記微小要素における流体の圧力、圧力変化及び流動速度の少なくとも1つを算出する算出手段と、
算出された流体の圧力、圧力変化及び流動速度の少なくとも1つにより解析された解析結果を出力する出力手段と
して機能させるための流体流動過程の解析プログラムにおいて、
前記流動コンダクタンス決定手段は、前記金型により製作される製品の肉厚を規定する前記キャビティの幅の大きさと、前記流体に含有される充填材の代表長さとの関係式を含む関数方程式を解くことにより決定することを特徴とする流体流動過程の解析プログラム。 - 前記流動コンダクタンス決定手段は、下記の関数方程式を解くことにより前記微小要素における流体の流動コンダクタンスを決定することを特徴とする請求項5に記載の流体流動過程の解析プログラム。
C×((∂2κ/∂x2)+(∂2κ/∂y2)+(∂2κ/∂z2))=−(1/η)
但し、
C=α×exp(L/t)+β×(L/t)+γ
で示され、
α×exp(L/t)+β×(L/t)+γ<1
となるとき、C=1となる(ηは流体の材料粘度、x、y及びzは微小要素の位置、κは流動コンダクタンス、Cは流動抵抗、Lは充填材の代表長さ、tはキャビティの幅の大きさ、α、β及びγは実験的に求められた定数である。)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009050927A JP2010201822A (ja) | 2009-03-04 | 2009-03-04 | 流体流動過程の解析装置、解析方法及び解析プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009050927A JP2010201822A (ja) | 2009-03-04 | 2009-03-04 | 流体流動過程の解析装置、解析方法及び解析プログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2010201822A true JP2010201822A (ja) | 2010-09-16 |
Family
ID=42963789
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009050927A Pending JP2010201822A (ja) | 2009-03-04 | 2009-03-04 | 流体流動過程の解析装置、解析方法及び解析プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2010201822A (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015066873A (ja) * | 2013-09-30 | 2015-04-13 | 本田技研工業株式会社 | コンピュータ支援による金型設計装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0899341A (ja) * | 1994-08-01 | 1996-04-16 | Toray Ind Inc | 流体流動過程の解析装置、解析方法、射出成形過程の解析装置、解析方法、射出成形品および射出成形品の製造方法 |
JPH10278086A (ja) * | 1997-04-04 | 1998-10-20 | Canon Inc | 薄肉部品の充填診断装置及び方法 |
JP2002321265A (ja) * | 2001-04-25 | 2002-11-05 | Toray Ind Inc | ウェルドライン予測方法および装置 |
JP2004025796A (ja) * | 2002-06-28 | 2004-01-29 | Canon Inc | 樹脂成形品の機械物性予測方法及び機械物性予測システム |
JP2006168300A (ja) * | 2004-12-20 | 2006-06-29 | Toray Ind Inc | 成形条件設定支援方法及び装置 |
JP2006523351A (ja) * | 2003-03-03 | 2006-10-12 | モルドフロウ アイルランド リミテッド | 加工される材料の性質を予測する装置および方法 |
JP2008200859A (ja) * | 2007-02-16 | 2008-09-04 | Toray Ind Inc | そり変形解析方法およびそのプログラムならびにそり変形解析装置 |
JP2008254396A (ja) * | 2007-04-09 | 2008-10-23 | Denso Corp | 成形品の収縮による変形量を予測する方法およびその装置 |
-
2009
- 2009-03-04 JP JP2009050927A patent/JP2010201822A/ja active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0899341A (ja) * | 1994-08-01 | 1996-04-16 | Toray Ind Inc | 流体流動過程の解析装置、解析方法、射出成形過程の解析装置、解析方法、射出成形品および射出成形品の製造方法 |
JPH10278086A (ja) * | 1997-04-04 | 1998-10-20 | Canon Inc | 薄肉部品の充填診断装置及び方法 |
JP2002321265A (ja) * | 2001-04-25 | 2002-11-05 | Toray Ind Inc | ウェルドライン予測方法および装置 |
JP2004025796A (ja) * | 2002-06-28 | 2004-01-29 | Canon Inc | 樹脂成形品の機械物性予測方法及び機械物性予測システム |
JP2006523351A (ja) * | 2003-03-03 | 2006-10-12 | モルドフロウ アイルランド リミテッド | 加工される材料の性質を予測する装置および方法 |
JP2006168300A (ja) * | 2004-12-20 | 2006-06-29 | Toray Ind Inc | 成形条件設定支援方法及び装置 |
JP2008200859A (ja) * | 2007-02-16 | 2008-09-04 | Toray Ind Inc | そり変形解析方法およびそのプログラムならびにそり変形解析装置 |
JP2008254396A (ja) * | 2007-04-09 | 2008-10-23 | Denso Corp | 成形品の収縮による変形量を予測する方法およびその装置 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015066873A (ja) * | 2013-09-30 | 2015-04-13 | 本田技研工業株式会社 | コンピュータ支援による金型設計装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Heller et al. | Planar deposition flow modeling of fiber filled composites in large area additive manufacturing | |
Tseng et al. | Improved fiber orientation predictions for injection molded fiber composites | |
Chen et al. | Rheological behavior of POM polymer melt flowing through micro-channels | |
Baruffi et al. | Effects of micro-injection moulding process parameters on accuracy and precision of thermoplastic elastomer micro rings | |
Modoni et al. | Application of the Digital Twin for in process monitoring of the micro injection moulding process quality | |
Marhöfer et al. | Gate design in injection molding of microfluidic components using process simulations | |
Wu et al. | In-line steady shear flow characteristics of polymer melt in rectangular slit cavities during thin-wall/micro injection molding | |
JP2010201822A (ja) | 流体流動過程の解析装置、解析方法及び解析プログラム | |
Shen et al. | Analysis for microstructure of microlens arrays on micro-injection molding by numerical simulation | |
JP2009271781A (ja) | 成形品の変形形状の予測方法とその装置、変形形状の予測プログラムとその記憶媒体 | |
JP2007083602A (ja) | 射出成形品の成形収縮率予測方法 | |
Mani et al. | Process parameter effects on dimensional accuracy of micro-injection moulded part | |
Abdullahi et al. | Process development and product quality of micro-metal powder injection molding | |
JP4643373B2 (ja) | 押出スクリュの均一溶融シミュレーション方法およびその方法を実行するコンピュータプログラムおよび押出スクリュの均一溶融シミュレーション装置 | |
JP5134259B2 (ja) | そり変形解析方法およびそのプログラムならびにそり変形解析装置 | |
Nakhoul et al. | A multiphase Eulerian approach for modelling the polymer injection into a textured mould | |
JP2013166282A (ja) | 繊維フィラー含有熱可塑性樹脂材料の流動過程の解析方法 | |
JP2003271678A (ja) | 数値解析方法および装置 | |
Islam et al. | Injection moulding simulation and validation of thin wall components for precision applications | |
JP2998596B2 (ja) | 流体流動過程の解析装置、解析方法、射出成形過程の解析装置、解析方法、射出成形品および射出成形品の製造方法 | |
Giboz et al. | Real-time analysis of polymer flow under real processing conditions applied to microinjection molding | |
JP2017149020A (ja) | 成形品の粘弾性構造解析装置、方法、プログラム、および媒体 | |
Willems et al. | Prediction of the mechanical properties of long fiber reinforced thermoplastics | |
Hopmann et al. | Development of a methodical approach to set-up the injection velocity profile dependent on the part geometry | |
JP4844421B2 (ja) | 射出成形の流動解析方法及びその装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20100728 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20100728 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120119 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130718 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130723 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20131119 |