(57)【要約】鋳造／成形過程１０２にて低圧縮性粘性流体の振舞いをシミュレートする方法である。この方法は、流体の運動量及びエネルギを解くため古典的流体力学方程式及び制御容積有限要素技術を使用する。この方法は、シミュレーションする間、ショットスリーブ１０９内の流体の流れ、収縮巣の形成、ダイと熱伝導流体との間の熱伝導、ダイにおける噴霧潤滑剤の冷却効果及びメンドラインの形成を選択的に考慮に入れる。
2004506515.tif 000002
【特許請求の範囲】
【請求項１】
鋳造／成形過程に対する熱流れをシミュレーションする方法において、
（ａ）シミュレーション過程に組み込まれる少なくとも１つのモジュールを提供し、該モジュールが、ショットスリーブ内の流れモジュール、収縮巣予測モジュール、ダイ内の熱伝導流体管モジュール、ダイ潤滑剤冷却モジュールを含む群から選ばれることと、
（ｂ）流れ領域を有限要素に離散化することと、
（ｃ）前記流れ領域内の質量の釣合い状態を更新することと、
（ｄ）速度及び圧力の場を更新することと、
（ｅ）温度の場を更新することとを備える、鋳造／成形過程に対する熱流れをシミュレーションする方法。
【請求項２】
鋳造／成形過程にて熱流れをシミュレーションする方法において、
シミュレートするため、鋳造／成形過程の少なくとも１つの成分を選び、該成分が、ショットスリーブ内の熱流れと、収縮巣の形成と、熱伝導流体の加熱する熱流と、該潤滑剤の冷却効果と、メンドラインの位置とから成る群から選ばれるようにすることと、
計算領域を要素のメッシュに離散化することと、
前記計算領域内の材料の釣合い状態を更新することと、
圧力及び速度の場を更新することと、
温度の場を更新することと、
シミュレーション停止基準に到達するまで、材料の釣合い状態、圧力、速度及び温度の前記更新を繰り返すこととを備える、鋳造／成形過程にて熱流れをシミュレーションする方法。
【請求項３】
請求項２の方法において、前記圧力及び速度の場を求めることが、有限要素方程式に従って離散化した運動量方程式及び質量方程式を誘導することと、等順位法に従って前記方程式を解くこととを備える、方法。
【請求項４】
請求項３の方法において、前記有限要素方程式が制御容積に基づく有限要素方程式である、方法。
【請求項５】
請求項２の方法において、ショットスリーブ構成要素内の前記熱流れがシミュレートされ、ラムが速度制御状態にあるとき、前記ショットスリーブ内のラムの速度から前記速度及び圧力が決定され、前記ラムが圧力制御状態にあるとき、ラムの力と前記流体からの抵抗力との釣合いから前記速度及び圧力が決定される、方法。
【請求項６】
請求項２の方法において、前記計算領域の前記温度の場を求めることが、潤滑剤の噴霧の冷却効果を決定することを更に備え、経験的に求められた冷却係数を自動的に識別し且つ前記冷却係数から熱の損失を計算することにより、潤滑剤噴霧の前記冷却効果が決定される、方法。
【請求項７】
鋳造シミュレーションにショットスリーブ内の熱伝導及び流体の流れを組み込む方法において、
前記ショットスリーブ及びラムを流れ領域内に組み込み、前記ラムの動作が速度又は圧力の何れかによって制御されるようにすることと、
前記ラムに隣接する要素に亙る質量熱流を決定し、前記ラムが速度制御状態にあるとき、該質量熱流の決定が該ラムの速度に基づき、前記ラムが圧力制御状態にあるとき。前記質量熱流の決定が、該ラムの押し力と前記流体の抵抗力との間の釣合いに基づくようにすることと、
圧力及び速度を更新することと、
前記ラムと、前記ショットスリーブと前記流れとの間の不一致的な熱伝導を想定することにより、温度の場を更新することとを備える、方法。
【請求項８】
請求項７の方法において、前記ラムの前記押し力が付与された液圧圧力及び液圧シリンダの断面積から計算される、方法。
【請求項９】
請求項８の方法において、前記圧力を前記更新することが、前記ラムの力をラムの表面積で割ることにより、前記ラムに隣接する要素の結節点に対する圧力境界条件を計算することを更に備える、方法。
【請求項１０】
請求項９の方法において、前記圧力を前記更新することが、第１の位置に対して前記圧力境界条件を使用することと、流体の排除量を計算することと、前記排除量をラム速度の境界条件として使用することと、前記ラムの前記押し力が前記流体の前記抵抗力に等しくなるまで、繰り返す工程を反復することとを更に備える、方法。
【請求項１１】
請求項７の方法において、前記不一致的な熱伝導が記録法により補正される、方法。
【請求項１２】
請求項１１の方法において、前記温度を更新することが、前記ラムに隣接した前記ショットスリーブの第１組の要素を識別し且つ該第１組の要素を熱伝導の計算に含めることを備え、前記温度を更新することが、前記ラムにより完全に占められた前記ショットスリーブの第２組の要素を識別し且つ該第２組の要素を前記計算領域から排除することを更に備える、方法。
【請求項１３】
請求項１２の方法において、前記第１組及び第２組の要素が記録技術によって識別される、方法。
【請求項１４】
請求項７の方法において、速度プロフィールを提供することを更に備え、前記ラムに対する前記休止圧力プロフィール、前記速度プロフィール及び前記プロフィールがユーザによって提供される、方法。
【請求項１５】
請求項１４の方法において、前記ラムの前記速度制御状態から圧力制御状態への切り換えが開始休止状態により決定される、方法。
【請求項１６】
請求項１５の方法において、前記開始休止状態が金型の充填割合であり、該割合が９９％充填である、方法。
【請求項１７】
請求項１５の方法において、前記開始休止条件が時間である、方法。
【請求項１８】
請求項１５の方法において、前記開始休止条件が前記ラムの位置である、方法。
【請求項１９】
溶融体の鋳造金型充填シミュレーションに巣を組み込む方法において、
巣の効果を連続性制約内に組み込むことと、
圧力休止を開始することと、
巣の形成圧力よりも低い圧力を有する要素を識別することと、
前記低下した圧力勾配に基づいて溶融体の速度及び正味流入量並びに巣の増分的容積割合を決定することと、
巣を更新することとを備える、溶融体の鋳造金型充填シミュレーションに巣を組み込む方法。
【請求項２０】
請求項１９の方法において、
前記改変することが、見掛けの密度βを連続性制約内に組込み、該改変した連続性制約が次式で表現されるようにする、方法。
【数１】
2004506515.tif 000037
ここで、
β＝（１−φ）ρ合金の見掛け密度
ｘ、ｙ、ｚ＝座標系の３つの成分
τ＝時間
ｕ、ｖ、ｗ＝３つの速度ベクトル成分
φ＝密度の容積分（０−巣無し、１−１００％多孔度）
ρ＝温度の関数としての合金密度である。
【請求項２２】
請求項２１の方法において、前記巣を形成する圧力が、相応する温度にて前記溶融体中に溶融した気体が核形成するときの圧力である、方法。
【請求項２３】
請求項２１の方法において、前記巣を形成する圧力が雰囲気圧力である、方法。
【請求項２４】
請求項２３の方法において、前記更新することが、次式によって表現される増分的巣を計算することを備える、方法。
【数２】
2004506515.tif 000038
ここで、
ρｓは、固相線温度における理論密度
ρは、所定の温度における理論密度
Ｖは、制御容積の容積
ν’（実際は、「ν」の上に「・」が付されている）は、制御容積に入ってくる質量の正味流量
ｄτは、増分時間ステップである。
【請求項２５】
請求項２４の方法において、前記巣を更新することが、液体溶融体を保持するが、所定の温度以下の温度を有する溶融体を保持する要素によって取り巻かれた要素を識別することと、前記液体溶融体の凝固に起因する巣を計算することと、該計算された巣を全体の巣に追加することとを更に備える、方法。
【請求項２６】
請求項２５の方法において、前記所定の温度が前記溶融体が凝固するときの温度である、方法。
【請求項２７】
請求項２６の方法において、前記巣を更新することが、前記休止段階の後、前記金型内非零液体溶融体の割合を点検することと、前記非零の液体溶融体の凝固に起因する巣を計算することと、該計算した巣を前記全体的な巣に追加することとを更に備える、方法。
【請求項２８】
請求項２７の方法において、前記巣を更新することが、前記負の圧力を有する前記制御容積の各々に対する正味流入量を計算することを更に備える、方法。
【請求項２９】
請求項２８の方法において、増分時間だけ時間を進めるステップと、前記更新された見掛けの密度βを使用して速度及び圧力の場を計算し、これにより前記制御容積の各々における前記収縮巣の時間経歴を確立するステップとを更に備える、方法。
【請求項３０】
ダイと熱伝導流体（ＨＴＦ）との間の熱伝導を熱流れ中に組み込む方法において、
ほぼ円形の断面幾何学的形態を有するものとしてＨＴＦ管をモデル化し、速度及び温度並びにＲｅ、Ｐｒ、Ｎｕ数の間、解析／経験的な相関関係に対する十分に確立された解析解法が存在するようにすることと、
前記ＨＴＦとの間の熱伝導を１次元的問題としてモデル化することと、
前記ＨＴＦの熱伝導係数を決定することと、
温度を更新することとを備える、方法。
【請求項３１】
請求項３０の方法において、前記モデル化することが、前記ＨＴＦ管を複数の部分に分割することを備え、該部分の各々が上流の結節点及び下流の結節点により境が設定され、前記ＨＴＦの温度が前記結節点にて評価され、各部分内の前記ＨＴＦの温度分布を前記２つの結節点から線形補間することができ、
１組のダイ要素の表面と前記ＨＴＦ管の部分の各々とを関係付けることを備える、方法。
【請求項３２】
請求項３１の方法において、前記モデル化することが、次のようにして前記ＨＴＦと前記ダイとの間の熱伝導を制御する１次元的エネルギ方程式を規定することを更に備える、方法。
【数３】
2004506515.tif 000039
ここで、
Ａｆｌｏｗ＝ＨＴＦ管の断面積
Ｌ＝流れ部分の長さ
Ｃｐ＝比熱
ρ＝ＨＴＦの密度
Δτ＝時間ステップ
ｈｋ＝ＨＴＦとダイ部分ｋとの間の対流熱伝導係数
Ａｋ＝部分ｋにおけるＨＴＦ管の円形面積
Ｔｍｋ＝部分ｋにおける平均ダイ温度
Ｔ０ｊ＝時間ステップの終了時の結節点ｊにおけるＨＴＦの温度
Ｔ０ｊ＝時間ステップ開始時結節点ｊのＨＴＦの温度である。
【請求項３３】
請求項３２の方法において、ＨＴＦの前記密度、粘度及び熱伝導率が前記ＨＴＦの温度から決定される、方法。
【請求項３４】
請求項３２の方法において、前記ＨＴＦの前記粘度、密度及び熱伝導率が検索テーブルから決定される、方法。
【請求項３５】
請求項３３の方法において、前記熱伝導係数が、ＨＴＦのナッセル数から決定される一方、ＨＴＦナッセル数が．前記ＨＴＦのプラントル数及びレイノルズ数から決定され、前記プラントル数及び前記レイノルズ数が前記ＨＴＦ管の前記粘度及び直径から決定される、方法。
【請求項３６】
請求項３５の方法において、前記ナッセル数が解析方程式に基づいて前記Ｒｅ及びＰｒ数から決定される、方法。
【請求項３７】
請求項３５の方法において、前記ナッセル数が経験的相関関係に基づいて前記Ｒｅ及びＰｒ数から決定される、方法。
【請求項３８】
請求項３５の方法において、前記熱流が、前記部分の表面積、前記ダイ要素と前記ＨＴＦの熱係数及びその平均温度差から決定される、方法。
【請求項３９】
請求項３５の方法において、前記ＨＴＦ管の前記結節点における温度が、次のものを備える三角形マトリックスを使用して前記エネルギ方程式を解くことにより同時に決定される、方法。
（ａ）前記輸送方程式の項を次のように規定し、
【数４】
2004506515.tif 000040
【数５】
2004506515.tif 000041
【数６】
2004506515.tif 000042
（ｂ）冷却通路の温度を次のように表現し、
【数７】
2004506515.tif 000043
ここで、
【数８】
2004506515.tif 000044
【数９】
2004506515.tif 000045
【数１０】
2004506515.tif 000046
である。
【請求項４０】
請求項３９の方法において、前記組みのダイ要素の表面を亙る前記熱伝導の全てを積分し且つ吸収された前記エネルギ量を時間の関数としてプロットすることにより、前記ＨＴＦにより吸収された熱エネルギの全体量を決定することを更に備える、方法。
【請求項４１】
請求項３０の方法において、前記全体的な熱流体の流れシミュレーションの前記温度の場を更新することを更に備える、方法。
【請求項４２】
ダイにおける潤滑剤の噴霧の冷却効果を溶融体の鋳造シミュレーションに組み込む方法において、
噴霧パラメータと、噴霧ノズルの動作プロフィールとを備える過程パラメータを提供することと、
前記動作プロフィールに基づいて前記ノズルの位置を更新することと、
前記ノズルの噴霧円錐体及び該噴霧円錐体によりカバーされるダイ要素の面を決定することと、
形成されたモデルに基づいて冷却係数を自動的に識別することと、
前記冷却係数から前記ダイ表面から失われた熱を決定することと、
前記ダイの温度の場を決定すべく前記失われた熱を前記ダイのエネルギ方程式に組み込むこととを備える、方法。
【請求項４３】
請求項４２の方法において、前記噴霧の前記動作が局所的座標系内で規定され、前記ダイの有限要素メッシュが前記局所的座標内で変換される、請求項６ａの方法において、前記局所的噴霧座標系が、前記ダイを構成する３次元的メッシュの要素を前記局所的噴霧座標内に配置する変換マトリックスを更に備える、方法。
【請求項４４】
請求項４３の方法において、前記局所的噴霧座標系が、原点と、ｘ軸と、ｙ軸と、ｚ軸とを備え、前記原点が前記噴霧ノズルにて画成され、前記ｘ軸が前記噴霧ノズルの長軸により画成され、前記ｙ軸が前記噴霧ノズルの短軸により画成され、前記ｚ軸が前記ノズルから前記ダイまで伸びる軸により画成される、方法。
【請求項４５】
請求項４４の方法において、前記噴霧ノズルの前記動作が、その当初の位置と、方向ベクトルと、その後の偏位ベクトル及び回転ベクトルにより画成される、方法。
【請求項４６】
請求項４５の方法において、前記噴霧ノズルの前記動作が、色々な時間間隔における連続的な位置ベクトル及び方向ベクトルと、方向ベクトル及び後続の偏位ベクトル並びに回転ベクトルにより画成される、方法。
【請求項４７】
請求項４６の方法において、前記噴霧円錐体が、ノズルの形状及び向きを幾何学的に解析することにより決定される、方法。
【請求項４８】
請求項４２の方法において、噴霧円錐体によりカバーされる前記ダイ表面を決定することと、
目の位置に前記ノズルを配置し且つ前記局所的座標系の前記ｚ軸に沿って前記ダイの方向を見て且つ前記噴霧の前方にあり且つ前記噴霧に面する直角ベクトルを有するダイ要素の面を識別することと、
その他の面の後方に隠れた要素を排除することとを備える、方法。
【請求項４９】
請求項４８の方法において、前記要素を前記排除することが、
【数１１】
2004506515.tif 000047
であるならば、影響試験の１つの範囲を実行することを備える、方法。
【請求項５０】
請求項４９の方法において、前記冷却係数が検索テーブルから決定される、方法。
【請求項５１】
請求項５０の方法において、前記検索テーブルが次の手順に従って経験的に決定される、方法。
（ａ）次式により、前記ダイ要素の表面の前記投影面積に達する噴霧量を計算することと、
【数１２】
2004506515.tif 000048
Ｑｆ＝金型の表面に到達する噴霧量
Ｑｓ＝噴霧容積流れ
ここで、
ηｓ＝噴霧実施効率
ｆｈｐ＝視認空間の割合
（ｂ）前記既知の潤滑剤の量を加熱された平坦な板の上に噴霧し、該加熱された平坦な板の全ての熱物理的パラメータが既知であるようにすることと、
（ｃ）前記既知潤滑剤量の噴霧に起因する温度効果を測定することと、
（ｄ）前記使用された既知潤滑剤噴霧量に相応する冷却係数を計算することと、
（ｅ）前記使用された噴霧量に関する前記冷却係数を検索テーブルに記憶させる。
【請求項５２】
溶融体の流れメンドラインの位置を予測する方法において、
グラフを表示し得るようにされたプロセッサを提供することと、
溶融体の前面に粒子を形成することと、
前記溶融体が流れるとき、前記粒子を追跡することと、
粒子の位置を時間の関数として表示することと、
異なる記号を有する粒子が衝突する箇所として前記メンドラインの位置を視覚的に識別することとを備える、方法。
【請求項５３】
請求項５３の方法において、前記粒子が記号により表現され、前記記号が、異なる時点にて且つ異なる溶融体の前面にて形成された粒子について相違する、方法。
【請求項５４】
請求項５２の方法において、熱流れシミュレーションからのデータを提供することを更に備える、方法。
【請求項５５】
請求項５３の方法において、異なる記号にて表現された粒子の混合レベルから修復の質及び強さを評価することを更に備える、方法。
【請求項５６】
請求項５５の方法において、前記停止条件がユーザが決定したものである、方法。
【請求項５７】
請求項５５の方法において、前記増分時間が、ユーザが決定したものである、方法。
【請求項５８】
請求項５７の方法において、前記増分時間が前記熱流れシミュレーションにて使用される増分時間と同一である、方法。
【請求項５９】
請求項５２の方法において、前記自由面が前記溶融体の容積割合から規定される、方法。
【請求項６０】
請求項５２の方法において、前記粒子の位置を更新することが、次のような数学的表現から粒子の新たな位置を計算する手順を備える、方法。
【請求項６１】
請求項５２の方法において、前記溶融体の前面における前記粒子が十分であるかどうかを決定することを更に備える、方法。
【請求項６２】
請求項６１の方法において、前記粒子の十分性を決定することが、前記溶融体の表面に付着した粒子数を決定することを備える、方法。
【請求項６３】
請求項６２の方法において、前記粒子と溶融体の前面との間の距離が所定の許容公差範囲内にあるならば、粒子が溶融体の前面に付着したものと見なされる、方法。
【請求項６４】
請求項５２の方法において、前記材料の粒子の位置を前記更新することが、前記粒子の各々を新たな位置まで前進させるステップを備え、前記新たな位置が前記溶融体の局所的速度ベクトルにより決定される、方法。
【発明の詳細な説明】
【０００１】
【関連出願】
本発明は、その内容の全体を参考として引用し、本明細書に含めた、２０００年７月１日付けで出願された、「鋳造／成形法過程用の熱流れシミュレーション（ＴＨＥＲＭＡＬ　ＦＬＯＷ　ＳＩＭＵＬＡＴＩＯＮ　ＦＯＲ　ＣＡＳＴＩＮＧ／ＭＯＬＤＩＮＧ　ＰＲＯＣＥＳＳＥＳ）」という名称の同時所有する出願係属中の米国特許出願第０９／６１０，３９０号の一部継続出願である。引用した出願は、その内容の全体を参考として引用し本明細書に含めてある。
【０００２】
【発明の分野】
本発明は、ダイキャスティング過程、又は金型充填過程の間、流体又は溶融体の振舞いをシミュレーションする方法、特に、非圧縮性、粘性な溶融体の流れ及びダイ内での熱伝導をシミュレートする数値的方法に関する。
【０００３】
【発明の背景】
製造のための鋳造／成形過程において、成形物品の製品の質について、金型の設計及び加工条件を最適化することが重要である。従来、かかる最適化は、専門家の大雑把なやり方で繰り返して行われる多数回の試行錯誤後に実現される。コンピュータ技術の進歩に伴い、コンピュータシミュレーションによって鋳造／成形過程の間、ある材料の流れ及び熱伝導の振舞いを解析することが可能となっている。現実的にモデル化したとき、かかるシミュレーションの結果を使用して、鋳造振舞いを予測し、これにより製品及び金型の効率及び加工条件の最適化を改良することができる。
【０００４】
鋳造／成形過程における流体の流れをシミュレートするための従来の方法は、全体として、運動、連続性及びエネルギの方程式を成形部品の離散化したデルに適用することと、加工ラインに沿った速度、圧力及び温度の分布状態を得るためその数値計算を行うこととを含む。これらのシミュレーションにおいて、材料の熱物理的性質、材料の温度、金型の温度、及び射出速度等のような作動上の特徴がシミュレーションのために必要とされる入力である。これらの入力は、全体として、解析結果と実際のモデリングとを繰り返して比較することにより得られた実験上のノウハウに基づいて行われる概その推測値である。従来のモデルは、全体として、入力条件が適当であるか否かを決定するための手段は何も提供しない。このため、正確なシミュレーションを実現する前に、現実的な境界条件を供給することが必要である。従来の流体の流れモデルは、物理的境界条件を一層良く推定することを可能にするであろう多くの改善を欠如しており、従って、正確なシミュレーションを実現し得ないことがしばしばである。
【０００５】
望ましい１つの改良点は、ショットスリーブ内の流れをモデルに組み込むことである。従来のダイキャスティング法において、典型的に、プランジャ（又はラム）と、ショットスリーブと、ゲートと、ランナーとを有する射出装置から金型キャビティ内に溶融材料が供給される。溶融材料をショットスリーブからゲートを通じてランナー内に射出するため、圧力がプランジャを介して加えられる。プランジャ及び溶融体が移動すると、それらの間の境界面は移動する不連続部となり、このことは、数値モデルに不安定さをしばしば生じさせる。このため、多くの従来技術のモデルは、ゲートの入口にてシミュレーションを開始し、ラムによって加えられる温度及び圧力又は速度がゲートの表面に亙って均一で且つ一定であると推測する。実際には、ショットスリーブ内の溶融体の温度は、材料がショットスリーブ及びプランジャに対して熱を失うため、均一でも一定でもない。同様に、ショットスリーブ内の速度及び圧力の場は極めて複雑なものとなる可能性がある。より重要なことは、プランジャにより加えられる圧力、及びキャビティが充填された後、ショットスリーブ内に残るビスケットの厚さ及び温度は、成形し／鋳造した部品の巣の形成を抑制する加圧効果に大きな影響を与える。本発明の目的の１つは、鋳造／成形のシミュレーションに短いスリーブ内の熱の流れを組み込むため、効果的で高信頼性のアルゴリズムを提供することである。
【０００６】
別の改良点は、ダイと熱伝導流体（ＨＴＦ）との間の熱交換をモデル内に組み込むことである。金型の温度は、全ての金属の鋳造又はプラスチックの射出成形過程に共通する工程パラメータである。射出サイクル間の安定的なダイ温度を保つため、ＨＴＦ管は、ダイの本体内に埋め込まれることがしばしばである。加熱した油又は水を使用してダイを予熱し、始動時間を短くし且つそれに伴うスクラップを減少させることができる。連続運転中、ＨＴＦの温度は、ダイ内で循環的な熱勘定を保ち得るように調節することができる。過程のシミュレーションにダイとＨＴＦとの間の熱伝導を組み込むための多数の試みが為されている。最も厳しく且つリソースを要求する方策は、完全な３Ｄモデルを使用して管路内でのＨＴＦの流れ及びその対流熱伝導をシミュレートすることである。計算時間を短くするため、既存のモデルの殆どは、ダイ内部のＨＴＦの流れ又は温度変化を無視する。それに代えて、ユーザは、熱伝導管の要素の表面を選び且つ熱的境界条件として、ＨＴＦ温度及び熱伝導係数を割り当てることが必要とされる。ＨＴＦの粘度は、その温度と共に変化し、また、ＨＴＦの温度は、全ループに亙って変化するため、レイノルズ数及びプラントル数は、ＨＴＦ管に沿って一定ではない。このため、ユーザが適切な熱伝導係数を決定することは難しい。このため、複雑な計算を行わずに、ＨＴＦの温度を正確に更新し、また、ダイからＦへの熱流が正確に予測され得るように、ＨＴＦの熱伝導係数を決定し得ることが望ましい。
【０００７】
収縮巣の形成を予測する方法は、また、金属の鋳造のシミュレーションモデルに組み込むことが望ましい特徴である。収縮巣の形成は、金属鋳造過程、例えば、半固体形成法に共通の問題点である。溶融金属が溶融金属よりも低温度であるダイのキャビティを充填するとき、溶融金属は、ダイにその熱エネルギを失い且つ凝固し始める。液相が固体に変化するとき、その密度は増大する。キャビティ内に供給された追加的な材料が無いならば、金属の正味量は、キャビティの正味容積よりも小さくなり、収縮巣の形成を誘発する。その寸法及び分布に依存して、収縮巣は、特に、その伸び及び疲労強度の点にて製品の機械的性質を著しく損なう可能性がある。圧力密の製品の場合、収縮巣は、多数の鎖状通路を無秩序に形成して漏洩を生じさせる可能性がある。特に、例えば、制御アーム及びナックル、空調機のハウジング及び燃料噴射レールのような、構造用又は漏洩無しの製品を設計するとき、収縮巣を正確に予測することは有利な点である。
【０００８】
従来、現在技術のダイキャスティングのシミュレーションプログラムは、鋳造品の容積内の高温箇所（凝固した低温の金属により取り巻かれたより高温の非凝固金属）に基づき収縮巣の潜在的な位置を表示できるのみである。温度分布に加えて、プログラムの幾つかは、また、収縮巣を予測するための金属の冷却率も組み込んでいる。しかし、収縮巣の形成は、実際には、単に熱的因子にのみ関係するのみならず、流れ及び加圧状態にも関連している。ダイキャビティ内の溶融金属が例えば、ショットスリーブ内に残る溶融金属のような、加圧された溶融材料と直かに接触しているならば、容積収縮を補償し且つ収縮巣の形成を抑制するため追加的な材料をダイキャビティ内に押し込むことができる。必要とされる圧力は、凝固する材料の粘度、溶融材料の温度、及びダイキャビティの幾何学的形態に依存する。収縮巣を正確に予測し、従って、収縮巣の形成を防止するため、巣が質量保存方程式に組み込まれた加圧された熱の流れをシミュレートするモデルが必要とされるであろう。
【０００９】
任意のダイキャスティング過程又は金型充填過程にとって重要な別の特徴は、潤滑剤が噴霧されるとき、ダイ又は金属の冷却状態を監視することである。一般に、圧力鋳造過程において、金属がキャビティ内で凝固するとき、ダイは開いて部品を突き出す。ダイを解放し且つ引き抜き傷又ははんだを防止するため、ダイの壁に水系潤滑剤が噴霧されることがしばしばである。潤滑剤は、典型的に、９８％の水を含むから、この潤滑剤は、合金からダイに伝達された熱エネルギのかなりの量を除去する。ダイの温度及び合金とのその熱伝導を予測するため、モデルに潤滑剤の噴霧による冷却効果を組み込む必要がある。全体として、異なる角度又は方向から潤滑剤を噴霧する多数の可動又は静止型ノズルかせあるようにすることができる。場合によっては、これらのノズルは、潤滑剤をより均一に噴霧するため前後に動くことができる。更に、ノズルの型式、噴霧圧力及びダイの温度に依存して、形成される冷却効果の分布は均一でない。その結果、特に、複雑なキャビティの形状を有するダイにとって、ユーザがダイ又は金属表面における相応する冷却係数を規定することは極めて難しい。既存のプログラムの幾つかにおいて、ユーザは、噴霧によりカバーすることのできるこれらの要素の表面における冷却係数を手操作で選び且つ割り当てなければならない。ノズルが動く場合、冷却係数を設定するには長い時間がかかる。噴霧ノズルの幾何学的形態及びその動作から冷却係数を計算するモデルが望ましい。
【００１０】
メンドラインの位置及び取り込んだ空気を予測するモデルもまた、鋳造／成形模擬に組み込むのに望ましい特徴である。溶融材料がキャビティを充填するとき、溶融体の前面は、割れて、例えば、コア又は穴の周りを進み且つ再結合する。溶融体の前面は、ダイのスチール及び空気と接触しているため、その前面は、残留する潤滑剤又は水分によって過度に冷却し且つ（又は）汚れる可能性がある。更に、ダイの設計が不適切であるならば、キャビティは、漸進的に充填されず、溶融材料中に空気が取り込まれる可能性がある。その結果、溶融体の前面が合う何処の箇所にても、過早の凝固、取り込まれた空気又は不純物のため、境界面の結合程度が弱体化する可能性がある。
【００１１】
ユーザは、メンドラインを探知すべく色々な角度から充填パターンを示すことができるが、特に、部品が大型で且つ複雑である場合、このことは時間がかかり且つ信頼し得ない。更に、溶融体の前面が合った後、これらの前面は、そのとき、キャビティが完全に充填されていなかったならば、動きを続けよう。その結果、弱体化した結合境界面の最終位置が最初のメンド位置と異なる可能性がある。更に、材料が境界面を亙って強力に混合するならば、結合強度を向上させることができる。これらの潜在的な鋳造の欠点を定量的に識別し、防止可能な設計を設計段階にて早期に組み込むことができることは、設計段階にて有益なことであろう。さもなければ、加工範囲は極めて狭くなり、均一な質を保つことは難しくなるであろう。
【００１２】
このため、鋳造／成形過程にて非圧縮性、粘性な流体の流れの振舞いをコンピュータシミュレーションによって正確にモデル化することができるように、上述した改良点を具体化する新規な方法が必要とされている。
【００１３】
【発明の概要】
従って、本発明の１つの目的は、鋳造／成形過程にて非圧縮性、粘性な溶融体の流れをシミュレートする改良された方法を提供することである。
本発明の関連する目的及び有利な点は、以下の説明から明らかになるであろう。
【００１４】
本発明は、ナビエル−ストークス（Ｎａｖｉｅｒ−Ｓｔｏｋｅｓ）の方程式及びエネルギ方程式を満足させる、ダイキャスティング又は金型充填過程中の熱の流れをシミュレートする方法を提供することである。この方法は、ショットスリーブ内の流れ、熱伝導流体管内の熱交換、潤滑剤によるダイの冷却、収縮巣の形成、及びメンドラインの形成という５つの流れ及び熱モジュールを提供し、これらの各々が別個に又は共に鋳造／成形過程に組み込むことのできる選択可能な構成要素を表現するようにすることである。これらのモジュールを組み込むことは、境界条件を一層良く予測し、従って、より正確にシミュレートすることを可能にする。モジュールは、任意の組み合わせにて鋳造モデルの全体に組み込むことができる。
【００１５】
本発明の１つの面は、鋳造／成形過程に対する熱流れをシミュレートする方法であって、
（ａ）全体的なシミュレーション過程に組み込まれるシミュレーションモジュールであって、ショットスリーブ内の流れモジュール、収縮巣予測モジュール、熱伝導流体管モジュール、ダイ潤滑剤冷却モジュールを含み、結節点及び熱流面を有する要素に流れ領域を離散化し、最初のパラメータ値及び過程調節限定符号を入力するシミュレーションモジュールを選ぶことと、（ｂ）増分的な時間ステップだけ時間を増分させることと、（ｃ）材料の釣合い状態を更新することと、（ｄ）運動量方程式を形成し且つ解き、速度及び圧力の場を更新することと、（ｅ）エネルギ方程式を形成し且つ解き、温度の場を更新することと、（ｆ）パラメータが収束されているかどうかを点検し、否であるならばステップ（ｄ）に戻ることと、（ｇ）最後のシミュレーション事象が生じているかどうかを決定し、否であるならば、ステップ（ｃ）に戻ることと、（ｈ）メンドライン予測モジュールを含むプロセッサ後シミュレーションを実行すべきかどうか決定することとを含む方法である。
【００１６】
【好ましい実施の形態の説明】
本発明の原理の理解を促進する目的のため、図面に図示した実施の形態に関して以下に説明し、その説明のために特定の用語を使用する。しかし、これにより本発明の範囲を何ら限定することを意図するものではないことが理解されよう。図示した装置の任意の変更例及び更なる改変例、また、本発明が関係する技術分野の当業者が通常、案出されると考えられるものとして図示した本発明の原理の更なる適用例も包含されるものである。
【００１７】
本発明の図示した実施の形態は、鋳造／成形過程中、非圧縮性、粘性な流体又は溶融体の振舞いをコンピュータシミュレーションする数値的方法である。図示した方法は、古典的な流体力学モデルに基づくものであり、ユーザが規定した時間間隔にて制御方程式に対する答えを数的に得るものである（一般に、コンピュータである計算装置を使用して）。図示した方法は、金属のダイキャスティング過程及びポリマの成形過程における全ての熱流れに対して一般的に適用される。
【００１８】
本発明は、ショットスリーブ内の流れモジュール、ダイ内の熱伝導管モジュール、ダイの潤滑剤冷却モジュール、収縮巣予測モジュール、メンドライン予測モジュールという５つの新規なモジュールを提供し、これらモジュールの各々が鋳造／成形過程の１つの面をシミュレートするようにするもので、従来技術の方法を改良するものである。これらのモジュールは、制御する流体力学方程式を解くための現実的な境界条件を提供する。
【００１９】
熱流れシミュレーションの目的に依存して、ユーザは、選択的に、全て又は任意の数のモジュールをシミュレーションに組み込むことを選ぶ。各モジュールの順次的な作動は、全体的なシミュレーション工程と別個に個別のフローチャートに記載してある。これらのフローチャートにおいて、モジュールは、主として説明の目的のため、完全なサブルーチンとして記述されている。しかし、これらのモジュールが組み込まれる実際のシミュレーションにおいて、モジュールの作動を異なる程度にて積分し、全体的なシミュレーション工程とし且つ識別不能となるようにすることも可能であることは当該技術分野の当業者に明らかであろう。
【００２０】
本発明の図示した実施の形態は、プロセッサ利用の計算装置内で作動し得る設計とされている。典型的に、本発明の方法の論理ルーチンは、例えば、ハードディスク、光学ディスク、磁気テープ又はプロセッサがアクセス可能な任意の現在技術の記憶媒体のような記憶媒体にある。次に、図１を参照すると、本発明の方法を実行することのできる典型的なプロセッサシステムが図示されている。読み取り専用記憶装置（ＲＯＭ）及びランダムアクセス記憶装置（ＲＡＭ）を有するプロセッサ４を含む従来のコンピュータシステム２は、ケーブル及びその他の手段によって入力装置及び出力装置に作用可能に接続されている。図示するように、これらの入力装置及び出力装置は、ディスプレイ５、キーボード６、マウス７、プリンタ８、入力／出力記憶装置９、又はその任意の組み合わせを含むことができる。プロセッサ４は、作動システムと協働して、本発明の方法の論理ルーチンを実行し、入力装置からデータを受け取り且つ記憶し、計算を行い、データを出力装置に伝送する。これと代替的に、実行可能なプログラム、入力データ及び出力データはネットワークを通じてアクセス可能であるローカルシステム又は遠方のシステムに保存してもよい。従来のコンピュータシステムが図示されているが、本発明から逸脱せずに、その他のシステムにて置換してもよい。具体的には、本発明の図示した実施の形態は、ユニックス演算システムを使用してシリコングラフィックシステム内にて作動する。しかし、同様の計算能力、記憶容量及びグラフィックディスプレイを有するその他のシステムも使用可能である。
【００２１】
本発明方法の一実施形態は、流れ領域が例えば、速度、圧力及び溶融体の前面位置のパラメータのような可変値が画成される結節点（又はグリッド点）を有する有限要素グリッド又はメッシュに分割される有限微分又は有限要素理法に基づいて開発されたものである。これらの可変値は、例えば、圧力及び温度のようなスカラー、例えば、２次元／３次元的問題における速度のようなベクトルを含む。従来の有限微分又は有限要素理論に従って、可変値は、その一次結節点値又は形状関数を使用して補間することができる。
【００２２】
制御方程式は、いくつかの従来の有限微分法（Ｓ．Ｖ．パタンカー、ヘミェスフェアパブリッシングコーポレーション（Ｓ．Ｖ．Ｐａｔａｎｋａｒ，Ｈｅｍｉｓｐｈｅｒｅ　Ｐｕｂｌｉｓｈｉｎｇ　Ｃｏｒｐｏｒａｔｉｏｎ））により数値熱伝導及び流体の流れ）又は、例えば、リッツ、ガラーキンの最小二乗法、採取法又は一般的に、重み付けした残留有限要素法のような有限要素法に従って一次数学等式を形成し得るように離散化させる（Ｊ．Ｎ．リディ、マックグロウ−ヒルインターナショナルブックカンパニー（Ｊ．Ｎ．Ｒｅｄｄｙ，ＭｃＧｒａｗ−Ｈｉｌｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｂｏｏｋ　Ｃｏｍｐａｎｙ）による有限要素法への導入）又はブログレスインアストロノミカル及びエアロノティカルエンジニアリングのＶｏｌ．８６ｐｐ３０５−３２７、１９８３の「熱伝導方程式の容積制御式の有限要素方程式Ｃｏｎｔｒｏｌ−Ｖｏｌｕｍｅ　Ｂａｓｅｄ　Ｆｉｎｉｔｅ−Ｅｌｅｍｅｎｔ　Ｆｏｒｍｕｌａｔｉｏｎ　ｏｆ　ｔｈｅ　Ｈｅａｔ　Ｃｏｎｄｕｃｔｉｏｎ　Ｅｑｕａｔｉｏｎ」」。
【００２３】
形成される分割した運動量方程式及び質量保存方程式を解くための多数の方法がある。許容可能な誤りの範囲内で収束した答えを得るための任意の解法が提案された方法にて適用可能である。好ましい実施の形態において、離散化方程式は、等順次法にて解かれる。等順次法は、ガウス点速度の場を導入し、各結節点にて圧力及び速度を評価することを可能にする。この方法は、流体の速度及び圧力の正確な解答を提供し、流体の質量保存を計算装置のトランケーション限界（機械の丸めレベル）まで満足させることができる。等順位法については、以下に詳細に説明する。等順位法は、好ましいが、運動量方程式を解き且つ質量を機械の丸めレベルまで保存することができる任意の方法も使用可能である。
【００２４】
この方法において、制御方程式（１）から（３）は次のように表現することができる。
質量の保存
【００２５】
【数１３】
2004506515.tif 000003
【００２６】
運動量の保存
【００２７】
【数１４】
2004506515.tif 000004
【００２８】
エネルギ方程式
【００２９】
【数１５】
2004506515.tif 000005
【００３０】
ここで、τ＝時間
ρ＝密度
Ｐ＝圧力
Ｖ＝速度ベクトル
μ＝粘度
ｇ＝重力成分
Ｔ＝温度
ＣＰ＝熱容量
Φ＝粘性加熱源の項
ｑ´´＝熱流
【００３１】
【全体的なシミュレーション手順】
本発明の一実施形態による全体的なシミュレーション手順１００のフローチャートを示す図２に関して説明する。
【００３２】
段階１０１にて、流れ領域は、要素のメッシュに離散化させる。メッシュ及び初期パラメータ値及びその他の境界条件を表現するマトリックスを入力する。
【００３３】
流れ領域をモデル化するため、有限微分グリッド又は有限要素メッシュを使用することができる。有限微分グリッドは、座標に対して平行であり、従って複雑な領域をモデル化する自由度の無い正規線にて１つの空間を離散化る。このため、不規則的な有限要素メッシュが好ましい。有限要素メッシュは、従来の手順により生じさせることができる。かかる適宜な手順は、例えば、アラカワ（Ａｒａｋａｗａ）に対し１９９１年４月２４日付けで発行された、「３次元的幾何学的処理方法及び装置（Ｔｈｒｅｅ−Ｄｉｍｅｎｓｉｏｎａｌ　Ｇｅｏｍｅｔｒｙ　Ｐｒｏｃｅｓｓｉｎｇ　Ｍｅｔｈｏｄ　ａｎｄ　Ａｐｐａｒａｔｕｓ　Ｔｈｅｒｅｆｏｒ）」という名称の米国特許第５，０１０，５０１号及びメッシュカット（Ｍｅｓｈｋａｔ）に対し１９９６年９月３日付けで発行された、「物体のメッシュ項を生産させる方法及びシステム（Ｍｅｔｈｏｄ　ａｎｄ　Ｓｙｓｔｅｍ　ｆｏｒ　Ｐｒｏｄｕｃｉｎｇ　Ｍｅｓｈ　Ｒｅｐｒｅｓｅｎｔａｔｉｏｎｓ　ｏｆ　Ｏｂｊｅｃｔｓ）」という名称の米国特許第５，５５３，２０６号に開示されている。
【００３４】
有限微分及び有限要素の双方は、モデル化する流れ領域に依存して一次元的、２次元的又は３次元的とすることができる。２次元的流れの問題をモデル化するためには、図３及び図４に図示したものと同様の三角形有限要素のメッシュがその幾何学的自由さの点で好ましい。しかし、四辺形要素もまた適用可能である。次に、図３を参照すると、三角形有限要素１０は、３つの頂点結節点１２を含んでいる。流体の可変パラメータ（すなわち温度、速度及び圧力）は、要素１０の結節点１２にて評価される。等順位方程式に特定的であるように、ガウス点速度Ｖｇ及びガウス点圧力勾配▽Ｐｇが各要素内にて規定され、これらが同一の順位を有し且つ結節点１２の可変値から計算されるようにする。
【００３５】
中間側部点１４を中心点１６に連結することにより有限要素１０を３つの副制御容積（ＳＣＶ）１８に分割することができる。副制御容積１８の各々は、材料が１つの副制御容積から別の容積まで流れることができる４つの熱流面３０によりその境が設定される。頂点結節点１２を取り巻くこれらの副制御容積１８は、その結節点に対する制御容積（ＣＶ）を形成する。例えば、図４に図示するように、頂点結節点２２は、５つの副制御容積２３から２７によって取り巻かれ、影領域は、頂点結節点２２の制御容積２０を表わす。制御容積２０は、５つの副制御容積２３から２７から成るものとして図示されているが、副制御容積の一部のみを流体で充填し又は部分的に充填することが可能である。このため、計算領域内にはこれらの充填した制御容積及び部分的に充填した制御容積のみが含まれる。１つの所定のメッシュの場合、ＳＣＶの各々にて充填状態を画成することは、制御容積の全体における充填状況のみを画成する場合よりも予想される移動自在面に対するよりも大きい自由度及びより優れた解法を提供することができる。しかし、全制御容積の充填状況を画成することは、プログラムの構造を簡略化することになる。
【００３６】
３次元的問題をモデル化するため、利用可能な自動の四面体メッシュの自由さがある点にて、四面体要素の方が六面体要素よりも好ましい。図４に図示するように、線形の四面体要素４０は、４つの頂点結節点４２及び４つの面４３を有している。有限要素４０は、表面中心点４４（四面体要素４０の各面における中心点）を容積中心点４６（六面体の中心点）に連結することにより有限要素４０を４つの副制御容積（ＳＣＶ）４８に分割することができる。副制御容積４８の各々は、６つの熱流面を有する六面体である。
【００３７】
三角形及び四面体有限要素が図示され且つ好ましいが、その他の幾何学的形状の有限要素、その組み合わせ又は有限識別法を使用することもできることも理解すべきである。更に、簡略さのため線形有限要素を図示した。しかし、有限要素を任意の順位のものとすることが考えられる。
【００３８】
好ましくは、不必要な可変値を排除し且つ計算時間を短縮するため計算領域が画成されるようにする。計算流れ領域の境界は、米国特許第５，５７２，４３４号に開示されたように、ワング（Ｗａｎｇ）らに従って描くことができる。この米国特許の開示内容は、その全体が明示的に本明細書に含まれている。ワングは、流れ領域内の要素の充填状況に従って計算領域の境界を描く。ワングは、図６に図示した三角形要素に対し、空の要素６２はその副制御容積の全てが空要素であると規定し、充填した要素６４は、その副制御容積の全てが充填され、部分的に充填された要素６６は、その副制御容積の少なくとも１つ、但し全てはなく充填された要素である。更に、ワングは、流れ領域の結節点を内部又は外部結節点と特定する。１つの結節点が充填した副制御容積にのみ取り巻かれるならば、この結節点は「内部結節点」６７として特定される。その他の全ての結節点が「外部結節点」である。これらの外部結節点は、いくつかの群に更に分割される、すなわち、充填され且つ部分的に充填された副制御容積と関係した外部容積は、「自由面結節点」６８として把握され、自由面結節点に近接するが、それ自体充填した副制御容積と関係しない結節点は、「被駆動結節点」６９と呼ばれる。計算領域の境界７０は、空の副制御容積と関係した外部結節点を連結し且つ被駆動結節点を接続する線として規定される。計算流れ領域を画成するワングの方法が図示されているが、計算流れ領域を画成し且つ有意義な解を得ることを許容する任意のその他の方法が本発明の範囲から逸脱せずに使用可能である。
【００３９】
段階１０２において、ユーザはシミュレーションに組み込むべき熱／流れモジュールを選び且つシミュレーションに組み込むべき過程条件を選ぶ。ユーザは、シミュレートすべきサイクル数及び可変値パラメータに対する収束基準を更に選ぶ。
【００４０】
ユーザは、次の熱及び流れモジュールの１つ又は任意の組み合わせをシミュレーションに含むことを選ぶことができる。すなわち、ショットスリーブ２００内の流れ、収縮巣の予測３００、ダイ４００内の熱伝導流体管、ダイ潤滑剤の冷却５００、メンドライン予測６００である。次に、この方法は、鋳造／成形過程にて実際の手順に基づいてその選んだモジュールを全体的な解答手順に積分する。メンドライン予測モジュール６００を除いて、このモジュールは独立的に実行されることを理解すべきである。モジュールのステップ手順は、全体のシミュレーションモジュールのステップ手順間にて分割し、２つのステップ手順が相互に依存し且つ識別不能とならないようにする。
【００４１】
選ばれたモジュールの各々に対し、特定の過程パラメータ及び境界条件が必要とされる。例えば、ショットスリーブ内の流れジュール２００にラム速度及び圧力プロフィールが入力することが必要とされる。この方法は、ユーザが必要な条件を供給するように指令する。しかし、ユーザがそれを拒絶するならば、この方法は、その記憶バンクから不作動過程条件を提供する。
【００４２】
ユーザは、シミュレートすべきサイクル数を選ぶことができる。１つのサイクルは、ショットスリーブに対し新たな流体を供給することからダイキャビティを充填し且つ鋳造部品を最終的に突き出すまでの完全な充填事象である。この方法は、経過したサイクル数を決定するカウンタを提供する。ユーザは、対象とする過程の現象に基づいてシミュレートするためサイクル数を設定する必要がある。全体として、過程作動の始動時、構成要素が安定状態となるのを許容し得るようにいくつかのサイクルを通じて作動させる必要がある。
【００４３】
ユーザは、また可変パラメータの収束を決定する基準を選ぶこともできる。更に繰り返しても独立的な可変値の正確さが向上しないとき、繰り返し過程は収束したものとされる。実際上、この繰り返し過程は、所望の精度が得られたときに終了する。運動量方程式及びエネルギ方程式は、非線形であり、流体の性質は速度の場及び温度と共に変化するから、増分的な変化が無視し得る程度に小さくなるまで、繰り返しが必要とされる。このため、連続的な２つの繰り返し間の可変パラメータの増分的な差を収束基準として選ぶことが好ましい。しかし、その他の適宜な基準を選ぶこともできる。更に、ユーザはこの収束基準のための許容公差を選ぶことができる。
【００４４】
段階１０３において、可変パラメータの全てが初期化される。この方法は、サイクルカウンタ及びマスタークロックをシミュレーションの開始を表示する零にリセットする。
【００４５】
段階１０４において、この方法は、ユーザがショットスリーブ内の流れモジュール及びダイ内のＨＴＦ管モジュールを選んだならば、ショットスリーブ及びＨＴＦ管を含むキャビティの可変値を初期化する。
【００４６】
段階１０５において、その後の運転サイクルの各々にて、サイクルカウンタは１だけ進む。
【００４７】
段階１０６において、その後の各時間ステップに対し、時間τを増分間隔Δτだけ進める。ユーザは、Δτを固定し又はΔτが変化することを許容することも選ぶことができる。好ましくは、Δτは、変化し得るものとする。Δτを変化させることは、ユーザが各時間ステップにて可変パラメータの増分的変化を規定し、シミュレーションの全体に亙る可変パラメータの収束を僅かな回数、繰り返すことにより実現することができるようにする。
【００４８】
段階１０７において、ユーザがショットスリーブ内の流れモジュール２００を組み込むことを選んだならば、この方法は、プランジャ（又はラム）が圧力制御又は速度制御状態にあるかどうかを決定する。プランジャが速度制御状態にあるならば、この方法は、プランジャによって排除された流体の量を決定する。この流体の量は、ダイキャビティにおける質量の熱流を決定する。次に、この方法は、後続の時間ステップの準備としてプランジャの位置を更新する。プランジャが圧力制御状態にあるならば、この方法は、プランジャが流体に加える正味圧力（プランジャの押す力と流体の抵抗力との差）を決定し且つ境界面における速度を決定する。ユーザがショットスリーブ内の流れモジュール２００を選ばなかったならば、この方法は段階１０８に進む。
【００４９】
段階１０８にて、副制御容積の充填状態、溶融体の前面の位置、取り込まれた空気ポケット及びその他の新たな温度及び速度の場に基づく依存可能な変数が更新される。
【００５０】
副制御容積内の流体の瞬間的な量は、その初期の容積と最初から時間ステップが解析されるまでの溶融体の正味流入量の合計値である。質量熱流は、流体の密度、速度、熱流面積及び時間から計算される。結節点制御容積内の溶融体の量は、供給された全ての副制御容積内溶融体合計量である。溶融体の充填部分は、制御容積の各々を追跡することで計算することができる。この方法の１つの好ましい実施の形態において、結節点の制御容積ｊの周囲質量熱流の寄与分はガウス点速度の場（ｕｇ、νｇ、ｗｇ）を使用して計算されることが好ましい。
【００５１】
【数１６】
2004506515.tif 000006
【００５２】
ここでＪｊは、結節点制御容積ｊの周りの質量熱流、ρは密度、ωは収束表面積である。ガウス点速度の場を使用して質量熱流を計算すれば、大局的質量保存を計算装置の截頭限界値まで計算することが可能となる（機械の丸めレベル）。
【００５３】
更に、この方法は、大局的質量の釣合い状態が保存されることを確実にすべく質量熱流を計算するため重み付けスキームを使用する。結節点制御容積を見出す質量熱流の場合、上流の充填部分を使用して充填束に重み付けする。ある質量熱流が結節点制御容積を排出するならば、次に、下流の充填部分を使用してその熱流に重み付けする。この重み付けスキームは、ソーラバッフ（ＳＯＬＡＶＯＦ）が使用するドーナーアクセプターアルゴリズムと全く同様である。
【００５４】
この段階にて、自由表面の位置も決定される。溶融体の容積部分、ｆｍを使用して次のように規定される自由表面の位置；ｆｍを推定する。
【００５５】
【数１７】
2004506515.tif 000007
【００５６】
ここで、Ｖ及びＶｍは、制御容積の容積及び溶融体が占める容積をそれぞれ表わす。溶融体の前面が存在する要素において、いくつかの結節点における溶融体の容積部分ｆｍは１であり、それ以外は１以下である。自由表面は、溶融体の容積部分を結節点にて補間した値が０．５である表面であると規定される。
【００５７】
この段階にて、流れ領域内に取り込まれた空気ポケットを識別することもできる。取り込まれた空気ポケットは、充填し又は部分的に充填した副制御容積によって取り巻かれた空の副制御容積の単一又は１つの群である。これらの取り込まれた空気ポケットは、その初期圧力と称され且つ雰囲気圧力が割り当てられる。これらの取り込まれた空気ポケットの位置及びその溶融体分は、後続の時間ステップの各々にて追跡し且つ更新される。取り込まれた空気ポケット内の圧力は、特定時間間隔間のポケット内への全体的な質量熱流を使用することにより理想気体法則に従って計算される。更新された圧力は、次の繰り返し時における任意の充填間事象の境界条件となる。
【００５８】
結節点制御容積の充填状況を更新し且つ溶融体の前面を決定した後、プログラムは、前に説明した規定に従って計算領域の境界を更新する。
段階１０９にて、この方法は、圧力及び速度の値を求めるためナビエルストローク方程式の近似解を求める。この方法は、独立した方程式を形成し且つその他の方程式を等順位法に基づいて解く。等順位法の展開例がその開示内容の全体を参考として引用し本明細書に含めた「流体流れシミュレーションのための等順位法（Ｅｑｕａｌ　Ｏｒｄｅｒ　Ｍｅｔｈｏｄ　ｆｏｒ　Ｆｌｕｉｄ　Ｆｌｏｗ　Ｓｉｍｕｌａｔｉｏｎ）」という名称の出願に開示されている。次に、この方法について簡単に説明する。
【００５９】
簡略化のため２次元的ケースについて考えると、組み立てた運動量等式は、等式（６）、（７）によりモデル化することができる。
【００６０】
【数１８】
2004506515.tif 000008
【００６１】
及び
【００６２】
【数１９】
2004506515.tif 000009
【００６３】
ここで、ａｉｊは、影響力係数マトリックスを構成する項、ｎは隣接する副制御容積の収束面の合計数である。この領域の整数外の圧力勾配を引けば、等式（６）、（７）は以下の等式（８）から（１２）に書き直される。
【００６４】
【数２０】
2004506515.tif 000010
【００６５】
【数２１】
2004506515.tif 000011
【００６６】
ここで、
【００６７】
【数２２】
2004506515.tif 000012
【００６８】
【数２３】
2004506515.tif 000013
【００６９】
【数２４】
2004506515.tif 000014
【００７０】
結節点のハット速度Ｕｉ、νｉは、隣接する副制御容積の結節点速度ｕｊ、νｊから成り、圧力項は含まない。ガウス点は、速度Ｕ＾ｉ、ν＾ｉ（実際は、「Ｕ」、「ν」の上にそれぞれ「＾」が付されている。）を有し、ガウス点の係数項Ｋｇは次のように規定される。
【００７１】
【数２５】
2004506515.tif 000015
【００７２】
【数２６】
2004506515.tif 000016
【００７３】
【数２７】
2004506515.tif 000017
【００７４】
ここで、ｎｅは、要素内の収束面の数である。ガウス点速度は、次式で表わされる。
【００７５】
【数２８】
2004506515.tif 000018
【００７６】
【数２９】
2004506515.tif 000019
【００７７】
ガウス点速度ｕｇ、νｇに対する上記の２次元的表現において、圧力勾配∂Ｐ／∂ｘ、∂Ｐ／∂ｙは、収束面にて一定である。整合性を図るため、項ｕｇ、νｇ、Ｋｇは、同様に収束面にて一定である。
【００７８】
ガウス点収束速度ｕｇ、νｇの項を連続等式（１）に代入することにより、圧力のポアソン方程式が得られる。圧力方程式は、次式で表示される。
【００７９】
【数３０】
2004506515.tif 000020
【００８０】
この方程式を要素の熱流面領域に亙って積分することにより、圧力の場は次のように表現することができる。
【００８１】
【数３１】
2004506515.tif 000021
【００８２】
グリーン公式を使用すれば、これら面積の積分値は、線積分値に変換される。
【００８３】
【数３２】
2004506515.tif 000022
【００８４】
上記に規定されたガウス点のハット速度Ｕ＾ｇ、ν＾ｇ、Ｋｇを使用することにより、ポアソン圧力方程式を組み立て且つ解いて圧力の場を決定することができる。
結節点速度ｕｉ、νｉ及びガウス点速度Ｕ＾ｇ、ν＾ｇは、新たに誘導された圧力の場及びＵ＾ｉ、ν＾ｉ、Ｋｉ、Ｕ＾ｇ、ν＾ｇ、Ｋｇを運動量方程式に代入することにより計算される。
【００８５】
等順位法が図示され且つこれが好ましいが、質量及び運動量方程式の独立した形態を誘導し且つ解くその他の方法も本発明の精神及び範囲から逸脱せずに、使用可能であることを理解すべきである。例えば、運動量方程式、質量保存方程式を組み合わせて１組の一次行列方程式とし且つ例えば、逆行列法のような直接的な解法又は例えば、共役勾配法のような繰り返し的な解法を使用して解くこともできる。
【００８６】
更に、巣予測モジュール３００が選ばれるならば、連続方程式は、合金の見かけ密度にβ＝（１−φ）ρ（ここでφは、巣の容積部分（０−巣無し、１から１００％多孔質）であるその容積部分の項にて表現された巣の効果を示すように改変される。
【００８７】
段階１１０、１１０Ａにて、巣予測モジュール３００が組み込まれるならば、この方法は、繰り返す間、結節点圧力が巣形成圧力（ＰＦＰ）よりも低いかどうかを点検する。是であるならば、圧力はＰＦＰにリセットされ、収縮容積と正味流入容積との差から同一の制御容積内の増分的制御容積が計算される。巣は更新される。ＰＦＰは、合金の各々に対する独特の値であり、ユーザが画成するものである。
【００８８】
巣予測モジュールが組み込まれない場合、又はＰＦＰ以下のモジュール圧力が存在しないならば、この方法は段階１１１へと進む。段階１１１において、この方法はプランジャが圧力制御状態にあるかどうかを点検する。是であるならば、この方法は、収縮を補償すべくダイキャビティ内に搾り出された材料の容積に基づいてプランジャの新たな位置を決定する。
【００８９】
段階１１２において、この方法は、適宜なエネルギ方程式（３）を解くことによりダイ及び溶融体双方温度の場を解くべく進む。ユーザが収縮巣モジュール３００及びショットスリーブ内の流れモジュール２００を組み込むことを選んだならば、ショットスリーブ及びプランジャと溶融体との間の熱伝導もこの温度の計算に含まれることになる。更新されたダイ温度を使用してダイとＨＴＦとの間の熱伝導係数を決定する。
【００９０】
エネルギ方程式は、有限要素方程式に従って離散化されることが好ましい。好ましくは、溶融体のエネルギ方程式は、慣性力、移流、拡散、容積熱源、粘性加熱、同時的境界面加熱の項を含むようにする。同様に追加的なディリヒレット及びニューマン境界条件を含む必要がある。ダイに対するエネルギ方程式は、伝導及び輻射に対する項を含むことが好ましい。組み立て過程の後、大局的エネルギ方程式の左側は、温度解ベクトルに基づいて作動する係数行列によって表現することができる。右側は、ソース項のベクトルから成っている。このエネルギ方程式の離散化した形態は、各結節点における温度を求めるため、解くことができる。図示した実施の形態において、離散化したエネルギ方程式は、繰り返して解かれる。溶融体及びダイ双方温度の場が更新される。各結節点における固体の割合は、温度対固体割合の曲線を使用して決定することができる。このような曲線は、経験的に求められ、当該技術分野にて全体として既知である。
【００９１】
段階１１３において、ユーザがシミュレーションにダイ内のＨＴＦ管モジュール４００を組み込むことを選んだならば、エネルギ方程式を解く前に、モジュール４００のステップ手順が組み込まれよう。流れ領域は、またダイ内の熱伝導管を含む。入口温度、流量及び流体の型式が規定される。モジュール４００は、ＡＴＦ管が円形であると想定し、また、確立された流れ及び熱モジュールを使用して熱流を計算する。モジュール４００の結果としてダイから熱伝導流体への熱流速は、ダイの温度を予測するために使用される。ＨＴＦ管モジュール４００の順序的ステップについて以下に説明する。
【００９２】
段階１１５において、この方法は材料の全ての性質を更新する。
段階１１６において、この方法は、速度及び温度の場の収束を点検する。この方法は、速度及び温度値の全てが、ユーザが供給した許容公差に収束するまで、流れ方程式及びエネルギ方程式の解の近似値を繰り返して求める。結節点又はガウス点速度又は温度Ｔｉの何れかが収束しないならば、この方法は、段階１０７に戻り、この段階にて、質量方程式、運動量方程式及びエネルギ方程式が形成され且つ更新された材料の性質を使用して再度、繰り返される。
【００９３】
段階１１７において、この方法は、充填停止又は休止停止状態にあるかどうかを点検する。ユーザは、時間、充填割合又は圧力とすることのできる充填停止又は休止停止を規定することができる。充填割合は、充填停止基準として好ましく、時間は好ましい休止停止基準である。停止時間に達しないならば、この方法は段階１０６にて別の時間間隔に対する増分時間に戻り、別の繰り返しと進む。流れシミュレーション停止時間に到達したならば、次に、シミュレーションの金型充填部分が完了する。しかし、ユーザがシミュレーションにダイ循環剤冷却モジュール５００を組み込むことを選んだならば、シミュレーションの熱的部分が続行する。
【００９４】
ステップ１１８において、部品は拒絶され、この方法はダイから噴霧潤滑剤への熱の損失を決定することを続ける。この熱流は、ダイの温度を計算する境界条件として使用される。モジュール５００の作動手順が完了したとき、シミュレーションが完了する。ダイ循環剤冷却モジュール５００の順次的な作動について以下に説明する。
【００９５】
段階１１９において、この方法は、サイクル状態の終了を点検する。ユーザが規定したサイクル数に到達しないならば、次に、更新された材料の特性がモデルに加えられ、新たなサイクルが開始される。
【００９６】
段階１２０及び段階１２１において、ユーザは、メンドラインの位置を決定するため、メンドライン予測モジュール６００を作動させ且つそれ以前のシミュレーションの修復の質を評価する。メンドラインの予測は、処理後の措置であり、全体的なシミュレーション手順１００の計算サイクルと相互作用しない。メンドラインの予測モデルの結果は、図で表示することができる。メンドラインの予測モジュール６００の順次的作動について以下に説明する。
【００９７】
【ショットスリーブ内の流れモジュール】
ショットスリーブ内の流れのシミュレーションは、キャビティの充填及び休止のためのシミュレーションと完全に結合されている。ショットスリーブ内の流れを組み込むことにより、我々は、初期温度の変化及びショットスリーブ内の供給材料の熱損失、不均一な速度、ゲートにおける温度の場、ビスケットの寸法及び凝固並びに休止圧力の降下に起因する予測の正確さを向上させることができる。
【００９８】
実際の過程制御の場合と同様に、ユーザは、ラム速度プロフィール又は射出圧力プロフィールを制御することを選ぶことができる。速度制御を使用するとき、ラムの速度は、隣接する結節点における速度境界条件として供給され、各時間間隔におけるラムの位置が監視される。ラムに近接する各制御容積に対し、ラムによって排除された流体の量が決定され且つ質量保存方程式における正味流入量の一部として組み込まれる。圧力制御が選ばれたならば、ラムの圧力は、最初の繰り返しの間、これら隣接する結節点における境界条件として供給される。次に、これらの隣接する結節点から流れ出る流体の全体量をラムの断面積で割ってその時間間隔におけるラムの有効変位を求める。その後の繰り返しにおいて、射出した材料からの全ラムの力が液圧圧力の力に等しくなるように排出量が調節される点を除いて、ラム速度制御の場合と同一のアルゴリズムを使用することにより、この排出量が境界条件として供給される。熱勘定のため、ラムの要素は、ラム自体と共に、移動する。容器ラムと接触したこれらのスリーブ要素は、入念な記録を通じて識別され、これらの間の熱伝導をエネルギ方程式に含むことができるようにする。ラムは、スリーブチャンバを通じて移動し、チャンバの多くの要素はラムによって占められ且つラムの後に残ることを理解すべきである。これらの要素は、計算領域から除外される。
【００９９】
次に、本発明のショットスリーブ内の流れモジュール２００の一実施形態のフローチャートを示す図７について説明する。モジュール２００は、完全な独立的なルーチンとして図７に図示されているが、このモジュールは、概念的理解及び明確化のために図示したものである。当該技術分野の当業者には、ショットスリーブ内の流れモジュール２００の順次的作動は、全体的なシミュレーション手順の間に散在されており、この手順が作動することによりその機能が依存することが容易に明らかであることが考えられ且つ理解される。
【０１００】
モジュール２００は、ショットスリーブ、全てのコア、挿入体、キャビティ、ゲート、ランナー及びラムを含む、流れ領域の幾何学的モデルを入力する段階２０１にて開始する。次に、流れ領域及びラムを複数の有限要素又は複数の積分グリッド内に組み込み、メッシュ／グリッドを表現する行列を入力する。段階２０２において、ユーザは、ラムの当初の位置及び速度又は圧力プロフィールを含む１組みの過程パラメータ及び境界条件を規定することができる。速度及び圧力プロフィール又はその双方の組み合わせは、時間、位置又はその双方の組み合わせとして規定することができる。ユーザが規定可能であるその他のパラメータは、ビレットの当初の容積及び温度分布を含む。好ましくは、ビレットの当初の位置及び温度分布は、ショットスリーブ及びキャビティ内の容積内で規定されるようにする。最後に、材料の熱物理的性質も入力される。
【０１０１】
段階２０３において、マスタークロックは、増分時間Δτだけ進む。時間クロックは、全体的な解手順１００の段階１０６のものと同一であり、増分時間Δτは、また全体的な解手順１００と同一であることも理解すべきである。
【０１０２】
各時間ステップにおいて、速度又は圧力の何れかがラムの動作を制御するが、その双方がラムの動作を制御することはない。このため、各時間ステップの開始時、この方法は、ラムの動作が速度により又は圧力によって制御されているかを最初に点検する。ラムが速度によって制御されるならば、ラムの速度は、境界条件としてラムに最も近い流体領域内の結節点に割り当てられる。ラムの動作が圧力によって制御されるならば、それ以前の時間ステップからのラムの速度が初期条件として使用されよう。その後の繰り返しにおいて、付与された液圧圧力に起因する押し力が溶融材料の抵抗力よりも大きいならば、ラムの速度は増し、またその逆にその抵抗力の方が大きいならば、ラムの速度は低下する。溶融材料に加わる平均圧力は、液圧圧力及び液圧シリンダの断面積から計算されたラムの押し力をラムの断面積で割ることにより計算することができる。この圧力は、ラムに最も近い熱流れ領域内の結節点における境界条件として付与することができる。
【０１０３】
段階２０４において、ラムの位置が更新される。ラムは、常にその軸方向に移動するため、ラムの位置は、ユーザが規定したラムの速度プロフィール又は所定のラム圧力により押し出された流体の容積及び時間増分Δτの何れかから記帳技術を使用して決定することができる。
【０１０４】
段階２０５において、この方法は、ラムの更新された位置及びショットスリーブの記録された行列を比較することにより、ラムの要素と直接接触したスリーブの要素を識別する。ラムに隣接したショットスリーブ表面における要素は、その後の熱伝導計算にて使用される。ラムによって完全に占められたチャンバの要素は、計算領域から除外されることが好ましい。
【０１０５】
段階２０６において、この方法は、全体的なシミュレーション手順１００の段階１０９と同一に作用する圧力及び速度の場を更新する。しかし、ショットスリーブ内の流れモジュール２００を組み込んであるため、速度及び圧力場計算のための境界条件は、改変されている。速度制御において、ラムの速度は、隣接する結節点における境界条件として使用される。前進するラムによって排除される溶融金属の量は、ショットスリーブの幾何学的形態及び増分的ラムの位置から計算される。この溶融金属の量は、制御容積の収束面における質量熱流に等しいようにされる。質量熱流は、圧力に対するポアソン方程式を形成するため、質量保存方程式に代入される。圧力方程式を解いて圧力の場を求める。圧力制御において、液圧圧力と液圧シリンダの面積との倍数であるラムの力は、ラムの面積で割り、また、最初の繰り返し位置の隣接する結節点における圧力境界条件として使用される。圧力及び速度が補正された後、隣接する制御容積の間でラムによって排除された流体の全体量を分割してその時間間隔におけるラムの有効排除量を求める。次の繰り返し時、この排除量は上述したように、ラム速度の制御の境界条件として使用されるが、各繰返しの後、この排除量が調節されて、射出した材料からの全圧力が与えられた液圧圧力により発生されたラムの全体力に等しいようにする点を除く。
【０１０６】
更新された圧力の場を使用して、速度の場を求めるため、大局的運動量方程式を解く。運動量方程式は、等順位法によって形成され且つ解かれる。等順位法であることが好ましいが、本発明の精神及び範囲から逸脱せずに、質量方程式及び運動量方程式を形成し且つ解くその他の方法も使用可能である。
【０１０７】
段階２０７にて温度の場を更新することは、全体的なシミュレーション手順１００の段階１１２と同一の工程である。この場合にも、モジュール２００を組み込むことは、温度の場を計算する境界条件を改変する。要素が移動するラムの要素と接触する箇所であるショットスリーブ内でこの改変が行われる。これらの非整合状態要素間の熱伝導は、境界面の境界を正確に記帳することで補償される。
【０１０８】
この方法は、段階２０８にて速度及び温度の場が収束しているかどうかを点検する（全体的なシミュレーション手順１００の段階１１６と同一）。速度の場又は温度の場が収束していないならば、この方法は、更新された圧力、速度及び温度値を使用して更に繰り返すべく、段階２０６に戻る。速度及び温度が収束しているならば、段階２０９にて、この方法は、速度から休止に切り換える基準にラムの制御が到達したかどうかを点検する。ユーザは、切り換えのための基準値を選ぶことができる。好ましくは、金型の充填割合が切り換えの基準として選ばれるようにし、より好ましくは、金型の９９％の充填割合が切り換えの基準として選ばれるようにする。これと代替的に、ユーザは、ラムの速度制御を圧力制御に切り換えるため、時間又はラムの位置を選ぶようにしてもよい。
【０１０９】
段階２１０にて圧力滞在状態にあるとき、上述したものと同一のアルゴリズムを使用して、圧力及び速度を解くことができる。ラムの動作は、典型的に、休止段階にて小さいため、ラムに隣接する結節点の間の圧力勾配は小さい。ラムの力及びラムの断面積から計算された圧力を、繰り返しの全体に対する境界条件として使用することができる。圧力の場が更新された後、ラムに隣接する制御容積から押し出された流体の全容積をラムの断面積で割ることによりラムの動作が計算される。
【０１１０】
段階２１１にて、この方法は、休止終了時に到達したかどうかを点検する。休止終了時間はユーザが選ぶ。休止終了時間に到達していないならば、この方法は段階２０３に戻って、時間増分Δτだけ進め、圧力及び速度の場が更新される。休止時、溶融体の流れに対する駆動力は、質量熱流ではなくて、圧力である。段階２１２において、休止終了時間に到達したことを点検した後、この方法は、全体的な解手順１００の段階１１８へと進む。
【０１１１】
【収縮巣予測モジュール】
収縮巣予測モジュール３００は、結合された熱―流れモデルを使用して収縮巣の形成を予測する。図示した実施の形態は、半固体形成法又はその他の鋳造法に適用可能である。最初に、熱―流れモデルは、金属及びダイが占める領域内の温度、圧力及び速度の予測値を提供する。質量保存の制約の下、合金が時間と共により高密度となる領域の結果、その周囲のマトリックス内の流れ抵抗が材料の充填を遅らせるのに十分に大きいならば、運動量方程式における圧力値は小さくなる。局所的圧力は、典型的に、極めて低く又は零に近い、巣を形成する圧力よりも低圧となるとき、更なる収縮は、その領域内に巣を生じさせることになろう。このため、収縮巣の形成及び移流は、質量保存方程式及び運動量方程式に巣及び密度の変化を組み込むことにより直接、計算することができる。この過程は、モデル内の各位置にて巣の時間経歴を形成するため各時間レベルにて行うことができる。
【０１１２】
次に、本発明の収縮巣予測モジュール３００による一実施形態のフローチャートを示す、図８について説明する。
段階３０１にて、ユーザは、収縮巣予測モジュール３００を選ぶ。この方法は、予測モジュール３００のステップ手順を全体的な流れモジュール１００のシミュレーションに組み込む。ユーザは、キャビティを充填する間のラムの速度／圧力に加えて、休止段階中の時間の関数として圧力プロフィールを規定することができる。
【０１１３】
運動量方程式を離散化して、速度成分の項を全体的なシミュレーション手順１００の段階１０９にて説明したように形成する。質量保存に対する巣の影響は、見掛けの密度β＝（１−φ）ρを規定することにより組み込み、ここでφは、巣の容積割合、及びρは、所定の温度における合金の理論密度である。巣の効果を含み得るように、連続性の制約を改変する。この改変された連続性の制約は、次式で表わされる。
【０１１４】
【数３３】
2004506515.tif 000023
【０１１５】
ここで、
β＝（１−φ）ρ合金の見掛け密度
ｘ、ｙ、ｚ＝座標系の３つの成分
τ＝時間
ｕ、ｖ、ｗ＝３つの速度ベクトル成分
φ＝密度の容積分（０−巣無し、１−１００％多孔度）
ρ＝温度の関数としての合金密度
【０１１６】
図示した実施の形態において、この方法は、ガウス点速度の場を使用することにより離散化した運動量方程式及び質量保存方程式を解くため、等順位方程式に従う。ガウス点速度ベクトルの項を改変した連続性の制約内に入れて圧力に対するポアソン方程式を形成する。この方程式は、計算領域内の各結節点の圧力を決定するために解く。等順位方程式であることが好ましいが、結合した運動量質量方程式及び運動量方程式を誘導し且つ解くためのその他の方法が本発明の範囲から逸脱せずに使用することも可能である。
【０１１７】
収縮巣の解析は、圧力休止の開始後に始まる。段階３０２において、圧力の場が解析される。制御容積の圧力が、各合金に独特の値であり、ユーザによって入力データの一部として規定される、巣の形成圧力（ＰＦＰ）と称される基準圧力以下に降下したとき、収縮巣が形成され又は増大する。ＰＦＰは、相応する温度にて合金内に融解した気体の各形成圧力として規定することができる。これと代替的に、ＰＦＰは、雰囲気圧力として規定することもできる。この段階において、ＰＦＰ以下の圧力を有するこれらの制御容積が識別される。
【０１１８】
例えば、部分液体を保持する領域が加圧状態から隔離されたならば、すなわちその領域を取り巻く材料が完全に凝固したならば、残る液体割合の相応する収縮が計算され且つ巣に加えられる。次に、この領域は、計算領域から排除される。材料は、固体と見なすことができ、所定の温度、例えば、ユーザが規定した固相線温度よりもその温度が低いとき、流れることはない。コンピュータの計算限界を上廻る材料の極めて大きい（無限）粘度を回避するため、本発明にてこの方策が好ましい。
【０１１９】
段階３０３にて、この方法は、ＰＦＰ以下の圧力を有するこれらの制御容積の圧力をＰＦＰにリセットし、その結果、いくつかの要素における圧力勾配を減少させる。
【０１２０】
段階３０４にて、減少した圧力勾配を使用して、この方法は、結節点速度を決定し且つ速度の場を更新する。
【０１２１】
段階３０５において、この方法は、減少した圧力勾配を使用して、ＰＦＰ以下の圧力を有する制御容積の各々にて正味容積流入量νを決定する。この減少した圧力勾配は、特定の制御容積内の圧力がより低圧であるときよりも少量の正味流量を生じさせる。このより少ない正味流量は、空隙が占める空間を考慮するものとして妥当なものとすることができる。
【０１２２】
段階３０６において、制御容積中の巣の増分量は、高密度化率とより少ない正味流れとの差を求めることにより計算される。制御容積内の巣の増分的容積分ｄφを推定するため、次の等式を使用することができ、
【０１２３】
【数３４】
2004506515.tif 000024
【０１２４】
ここで、
ρは、所定の温度における理論密度
Ｖは、制御容積の容積
ν・は、制御容積に入ってくる質量の正味流量
ｄτは、増分時間ステップ
【０１２５】
休止が停止したとき（ゲート又はラムにおける圧力が解放されるとき）、この方法は、流れ領域内に残る、零以外の任意流体割合の存在を点検する。
零以外の液体割合が残るならば、この方法は、相応する収縮を計算し、最終的な巣を更新する。
【０１２６】
段階３０７にて、各結節点の見掛けの密度βを更新した巣の値φを使用して計算する。更新した見掛けの密度βは、次の時間ステップにて速度及び圧力の場を計算するときに使用される。このため、制御時の各々における収縮巣の時間経歴を確立することができる。
【０１２７】
段階３０８にて、ショットスリーブ内の流れモジュール２００が同一の全体的なシミュレーション手順１００にて選ばれるならば、この方法は、最初に収縮を補償するため、キャビティ内に搾り出された材料の量に基づいて、ラムの新たな位置を更新し且つ全体的なシミュレーション手順１００の段階１０９へに進む。
【０１２８】
［ダイ内の熱伝導流体管モデル］
ダイ内のＨＴＦ管モジュール４００は、ダイの表面と熱伝導流体との間の熱伝導をモデル化し、温度の場が正確に決定され、また、処理条件が適正に設定されるようにする。２つの３次元物体間の熱伝導をモデル化することは、複雑であり、また、解が常に得られるとは限らない。このモデルは、ＨＴＦ管がほぼ円形の断面を有し、ＨＴＦの軸対称の速度及び温度分布を十分に確立された解析法に基づいて計算することができるものと想定する。簡単な計算技術により、ダイ表面とＨＴＦ管との間にて熱交換が生ずる箇所である境界面を正確に決定することができる。このように、このモデルは、３次元的問題を解法が制御可能である１次元的問題に簡略化する。ＨＴＦ管は、典型的に、円形の断面にて穿孔されるため、このモデルは、流れ領域のエネルギ勘定に対する現実の境界条件を提供する。　【０１２９】
図９には、ダイ要素９に対する関係にて熱伝導管６０の概略図６１が図示されている。熱伝導管６０の全体は、多数の線形部分ｋｊ（ここで、ｊ＝１、２、……Ｎ）に分割されている。部分ｋｊの各々は、ダイから冷却通路への熱伝導が計算される箇所である既知の組みのダイ要素表面６２と関係している。部分ｋｊの各々は、流れ方向ｍに対して名称が付された、上流結節点ｌｊ−１及び下流結節点ｌｊにより境が設定されている。部分ｋｊの各々における平均温度分布は、その両端の結節点温度から線形に補間される。
【０１３０】
この熱伝導管の制御方程式は、１次元エネルギ方程式の暗黙的な形態である。
【０１３１】
【数３５】
2004506515.tif 000025
【０１３２】
ここで、
Ａｆｌｏｗ＝ＨＴＦ管の断面積
Ｌ＝流れ部分の長さ
Ｃｐ＝比熱
ｐ＝ＨＴＦの密度
Δτ＝時間ステップ
ｈｋ＝ＨＴＦとダイ部分ｋとの間の対流熱伝導係数
Ａｋ＝部分ｋにおけるＨＴＦ管の円形面積
Ｔｍｋ＝部分ｋにおける平均ダイ温度
Ｔｊ＝時間ステップの終了時の結節点ｊにおけるＨＴＦの温度
Ｔ０ｊ＝時間ステップ開始時結節点ｊのＨＴＦの温度
次の項を規定することにより、
【０１３３】
【数３６】
2004506515.tif 000026
【０１３４】
【数３７】
2004506515.tif 000027
【０１３５】
【数３８】
2004506515.tif 000028
【０１３６】
冷却通路の温度Ｔｊを次式にて表すことができる。
【０１３７】
【数３９】
2004506515.tif 000029
【０１３８】
ここで、
【０１３９】
【数４０】
2004506515.tif 000030
【０１４０】
【数４１】
2004506515.tif 000031
【０１４１】
【数４２】
2004506515.tif 000032
【０１４２】
本発明の方法は、３次元解法を使用して上記のエネルギ方程式を数値的に解く。時間ステップの各々にて、レイノルズ数（Ｒｅ）、プラントル数（Ｐｒ）、ヌセルト数（Ｎｕ）及び熱伝導管の要素の各表面における相応する熱伝導係数ｈｋは、ユーザが規定した容積流量、ＨＴＦの熱物理的性質、及び更新したＨＴＦ温度に基づいて計算される。ＨＴＦとダイとの間で交換された熱流は、熱交換係数ｈｋ、ＨＴＦとダイとの間の局所的温度差及び時間ステップからダイ要素の表面６２の各々にて計算される。管部分ｋｊの各々内におけるダイ要素の表面の熱流にて、上流結節点ｌｊ−１におけるＨＴＦの温度、ＨＴＦの容積流量及び熱物理的性質、ＨＴＦの熱エネルギ、下流の結節点Ｔ１の温度である第１結節点温度Ｔ０＝Ｔｉｎｌｅｔの境界条件の設定値を計算することができる。この方策により、Ｒｅ、Ｐｒ、ＮｕとＨＴＦの流量、温度及び粘度の変化との間の解析／経験的相関関係を組み込んで、計算ステップを著しく増すことなく、熱伝導を正確に予測することができる。
【０１４３】
次に、本発明のダイ内熱伝導モジュール４００の一実施形態のフローチャートを示す図１０について説明する。
段階４０１にて、ユーザは、モジュール４００を全体的なシミュレーション手順１００内に組み込むことを選ぶ。ユーザは、ダイの有限要素モデルにおけるＨＴＦ管６０の寸法及び位置、流体の型式、その入口温度、容積流量及び初期温度を規定することができる。
【０１４４】
この方法は、ＨＴＦ管の入口及び出口のユーザが特定した要素の端縁とダイとの有限要素メッシュから、ＨＴＦ管６０と関係したダイ要素の面６２の全てを識別する。次に、これらのＨＴＦ要素の面は、その入口からの流れ距離に基づいて、多数の部分ｋｊ（ここでｊ＝１、Ｎ）に順次に群分けされる。この部分の数Ｎは、ＨＴＦ管の全長に基づいてユーザが決定する。
【０１４５】
段階４０３にて、この方法は、時間ステップ開始時熱伝導流体温度Ｔｋに基づいてＨＴＦの粘度μを更新する。最初に、Τｋは、ユーザ又はその方法の何れかが提供した境界値である。流れが進むに伴い、Τｋは、ＨＴＦがダイヤから吸収した蓄積熱から更新される。ＨＴＦの温度に基づいて、この方法は、部分ｋｊにおけるＨＴＦの密度、粘度μ及び熱伝導率を決定する。これらのパラメータは、これらのパラメータを一般的なＨＴＦの温度と相関させる検索テーブルから直接、読み取ることができる。このデータテーブルは、一般的に利用可能であり且つ当該技術分野の当業者に既知である。レイノルズ数、Ｒｅ＝ρνｄ／μ、及びプラントル数、Ｐｒ＝Ｃｐμ／ｋは、ＨＴＦ管の既知の性質及び直径から計算される。次に、相応するヌセルト数をレイノルズ数及びプラントル数から決定することができる。ナッセル数は、経験的な相関関数に基づいてＲｅ及びＰｒ数から決定することができる。これと代替的に、ナッセル数は、解析方程式に基づいて決定してもよい。Ｒｅ、Ｐｒ及びＮｕ数の間の関係は、モデルと独立的である。当該技術分野の当業者には、どの特定のモデルがその過程を最も良く表現するかが既知である。
【０１４６】
段階４０４において、この方法は、次の関係に従って、ナッセル数から部分ｋｊの熱伝導係数ｈｋを決定する。
ｈｋ＝Ｎｕｋ／ｄ
ここで、ｈｋ＝部分ｋにおけるＨＴＦの熱伝導率
ｋ＝熱伝導率
ｄ＝ＨＴＦ管の直径
【０１４７】
段階４０５において、この方法は、ＨＴＦ管の部分ｋと関係した組みのダイ要素６１の平均結節点温度Τｍｋを方程式（２３）に入力する。ダイ要素６１の結節点温度０は、その前の繰返し時のダイ温度場から読み出される。次に、熱伝導係数ｈｋ、ダイとＨＴＦとの間の温度差、ダイ要素の表面積及び時間から、各ダイ要素の表面６２におけるＨＴＦとダイとの間の熱流を決定することができる。
【０１４８】
この方法は、熱流の全てをダイとの境界面にて積分することによりＨＴＦ管内でＨＴＦにより吸収された熱エネルギの全体量を更に計算する。吸収された熱エネルギの量は、各時間ステップにて保存され且つ要求されたとき、時間の関数としてプロットすることができる。
【０１４９】
段階４０６において、この方法は、輸送方程式を構成し且つＨＴＦの全ての結節点温度を同時に解くため三角形マトリックスを形成する。１つの部分の下流結節点におけるＨＴＦの平均熱エネルギは、流入するＨＴＦの熱エネルギ及びダイから吸収されたエネルギを追加し且つ合計値をＨＴＦの質積流量で割ることにより、計算される。Τｊ＝０＝Τｉｎｌｅｔである境界条件を与えることにより、次に、各下流の平均温度は、その平均熱エネルギ及び熱性質、すなわち比熱から部分ごとに、順次に決定することができる。
【０１５０】
段階４０７にて、ＨＴＦ管の更新した結節点温度はダイとＨＴＦとの間の熱流を計算するための境界条件として使用される。新たに計算した熱流は、ダイの新たな温度を解くための熱境界条件としてダイのエネルギ方程式に組み込まれる。
【０１５１】
［ダイ潤滑剤冷却モジュール］
ダイ潤滑剤冷却モジュール５００は、部品を突き出した後、潤滑剤を噴霧することによりダイから除去された熱の量を計算する。この熱流は、ダイの新たな温度を解くための熱境界条件としてダイのエネルギ方程式に組み込まれる。このモジュールにおいて、潤滑剤の噴霧によりカバーされる要素の表面は、ノズルの位置から金型表面に向けて噴霧円錐体を放出することにより数学的に識別される。金型キャビティにおける要素の表面の相対的位置を噴霧円錐体と比較することにより、噴霧によりカバーされるこれらの要素の表面を識別することができる。噴霧ノズルが動くならば、噴霧円錐体及び相応する噴霧した要素の表面の位置は、各時間ステップに対し更新されよう。要素の表面における冷却係数は、経験的に決定される。実際の始動時、ユーザは、噴霧ノズル、ノズルの型式、噴霧圧力及び流量の位置、角度及び動作経路及び流量のみを規定すればよい。プログラムは、噴霧円錐体によって捕捉されたこれらの表面において、各時間ステップ時に、時間に依存する対流冷却係数を自動的に識別し且つ割り当てる。この関数により、ユーザは、ダイの温度及び製品の質に対するその影響を評価するため噴霧ノズル、噴霧パターン及び噴霧パラメータを迅速に変化させることができる。
【０１５２】
次に、本発明ダイ潤滑剤冷却モジュール５００の一実施形態のフローチャートを示す、図１１について説明する。モジュール５００は、シミュレーションの流れ部分に対する流れ停止時間に達し、ダイが開き且つ部品が突き出された後、開始する。
【０１５３】
ユーザが全体のシミュレーションにモジュール５００を組み込むことを選ぶならば、ユーザは、段階５０１にて噴霧パラメータを規定することができる。例えば、ユーザは、例えば、ノズルの形状及び寸法のような、噴霧ノズルのスタイル、ダイ表面に対するその最初の位置及びその動作プロフィールを特定することができる。ノズルは、任意の形状でよいが、円形又は楕円形の何れかであることが好ましい。更に、ユーザは、潤滑流体、その熱性質、流量、噴霧圧力及び噴霧時間のような、噴霧パラメータを特定することもできる。
【０１５４】
噴霧ノズルの動作は、その最初の位置及び方向ベクトル、並びにその後の偏位及び回転ベクトルにより、又は可変の時間間隔における連続的な位置ベクトル及び方向ベクトルにより規定することができる。その何れの場合でも、ベクトルは、全有限要素モデルの大局的座標又はダイ上に画成された局所的座標に対して規定することができる。
【０１５５】
局所的座標系が使用されるならば、ノズルの長軸及び短軸がＸ及びＹ軸を画成し、ノズルからダイの表面まで伸びる線がＺ軸を画成するように、ノズルがこの座標系の原点に配置されるようにすることが好ましい。次に、ダイの有限要素モデルを整合のため同一の局所的座標系に変換する。噴霧円錐体をモデル化するため上述の座標系も使用することができる。
【０１５６】
ダイ潤滑剤冷却モジュール５００を熱流れシミュレーションの大局的エネルギ勘定に一体化し、その作動順序を全体的なシミュレーション手順１００のマスタークロックによって制御する。段階５０２においてマスタークロックが時間間隔Δτだけ進むと、モジュール５００に対する新たな時間ステップも開始される。
【０１５７】
各時間ステップの開始時、段階５０３にて、各潤滑剤噴霧ノズルの位置が所定の動作プロフィールに基づいて更新される。段階５０４にて、ノズルの３次元噴霧円錐体の数学的表現をノズルの形状及びその向きから計算することができる。
【０１５８】
段階５０５において、噴霧円錐体中の要素が識別され、その他の面の後方に隠れた要素は排除される。上述の識別及び排除は、視覚的試験によって実現され、この場合、ノズルは眼の位置に置かれ、局所的座標系のＺ軸に沿ってダイの方向を見る。噴霧の前方にあり、また、噴霧に向けた正規ベクトルを有するダイ要素の面が選ばれる。選ばれた要素の面の上にて、第２試験が行なわれて、何れの要素の面が実際に噴霧した曝露され（すなわち、影響の円錐体の範囲内にあり）、その他のダイの作用部分から隠されないのかを決定する。
【０１５９】
この試験は次式となる。
【０１６０】
【数４３】
2004506515.tif 000033
【０１６１】
次に、要素の表面は、影響の円錐体の内側にあると決定される。否であるならば、その要素の表面は、その後のモデル化から排除される。下付き文字ｌは、噴霧により画成された局所的座標系を表す。
【０１６２】
段階５０６にて、各噴霧円錐体内部ダイ要素表面における視認領域（投影領域）が半球状の投影により決定される。半球状の投影技術は、輻射熱伝導の文献として豊富な文献に記載されており、当該技術分野の当業者に周知である。例えば、１９８５年にマグローヒルから出版された、Ｍ．Ｎ．オジシック（Ｍ．Ｎ．　Ｏｚｉｓｉｋ）、「熱伝導−基本的アプローチ（Ｈｅａｔ　Ｔｒａｎｓｆｅｒ−Ａ　Ｂａｓｉｃ　Ａｐｐｒｏａｃｈ）」には、半球状投影技術が記載されている。
【０１６３】
段階５０７にて、影響を受ける円錐体内のダイ要素の表面からの熱流が計算される。この熱流は、次のステップにダイの新しい温度を解くための熱境界条件としてダイのエネルギ方程式に組み込まれる。
【０１６４】
キャビティの表面からの熱流又は損失は、ダイ要素の表面の投射面積、ダイの温度、潤滑剤の温度、周囲温度及び噴霧パターンから決定することができる。一般に、熱流は、ダイの温度、噴霧量、噴霧流量、噴霧圧力等のような、噴霧パラメータの関数である熱伝導係数を使用して計算される。好ましくは、熱伝導係数は予め決め且つ検索テーブルに記憶させる。
【０１６５】
熱伝導係数テーブルは、経験的に構成することができる。次に、かかるテーブルを構成する一例を示す。最初に、要素の表面の投影面積に到達する噴霧量を次式で計算する。
【０１６６】
【数４４】
2004506515.tif 000034
【０１６７】
ここで、
Ｑｆ＝金型表面に到達する噴霧量
Ｑｓ＝噴霧の容積流れ
ここで、
ηｓ＝噴霧の実施効率
ｆｈｐ＝視認空間の割合
【０１６８】
次に、噴霧実験を行うことができ、この場合、既知量の潤滑剤が加熱された平坦な板に噴霧される。平坦な板の全ての熱物理的パラメータは既知である。好ましくは、この板は、ダイと同様の材料とする。噴霧に起因する温度降下を測定し、熱伝導係数を計算する。次に、ダイの要素の表面に到達する噴霧が既知であるとき、テーブルを使用して熱伝導係数を検索する。熱伝導係数用の検索テーブルを確立する方法について示唆したが、ダイ要素の表面の投影面積からの熱伝導係数を決定する、経験的又は数値的なその他の方法も使用することができる。
【０１６９】
熱流が計算され、ダイの温度が更新された後、段階５０８において、この方法は、ユーザが規定した噴霧停止時に達したかどうかを点検する。噴霧停止時に達したとき、その手順は、ある間隔だけ時間を増分する段階５０２に戻り、噴霧手順を繰り返す。噴霧停止時に達したならば、ダイ潤滑剤の冷却モジュール５００が計算され、シミュレーションは、全体的な手順１００の段階１１９まで続く。
【０１７０】
［メンドライン予測モジュール］
メンドライン予測モジュール６００は、流体流れ中のメンドラインの位置をシミュレートする。モジュール６００は、金型の充填シミュレーションの間に採取された流れ及び熱のデータからメンドラインの位置を予測する。このモジュールは、シミュレーションが完了した後に実行される。これは工程後の措置であり、全体的なシミュレーション手順１００の計算サイクルと相互作用しない。しかし、メンドラインを、流れシミュレーションが進行するに伴って追跡することが可能であると考えられる。
【０１７１】
溶融体の前面の自由面に沿って、材料粒子を表現する数学点が形成される。視覚的に、これらの粒子は、色及び記号の組み合わせで表現される。各粒子の位置が追跡される。各時間ステップにて、粒子はそれ以前の位置における計算された溶融体の速度に基づいて新たな位置まで進む。シミュレーションが完了した後、ユーザは発生時（追加されたとき）又は局所的な温度を表わす異なる色又は記号が割り当てられたこれら粒子の分布状態を表現することができる。このようにして、メンドラインを識別することができる。
【０１７２】
次に、本発明メンドライン予測モジュール６００の一実施形態のフローチャートを示す図１２を参照する。メンドライン予測モジュール６００が作動されたとき、ユーザは、例えば、単位面積当たりの粒子数のような目標粒子の密度を規定することができる。溶融体前面の当初位置にて、粒子は選ばれた目標粒子の密度に基づいて表面に形成される。更に、ユーザは、全体的なシミュレーションの範囲内で特定の時間間隔においてのみメンドラインの形成を予測すべくシミュレーションの開始時間及び終了時間を規定することができる。
【０１７３】
シミュレーションの開始時、段階６０１にて、タイムクロックが初期化される。プロセッサ後のシミュレーションであるため、モジュール６００の時間ステップは、熱流れの全体的なシミュレーションのマスタークロックによって制御されない。
【０１７４】
段階６０２において、タイムクロックは、時間間隔Δτまで進む。この時間間隔は、ユーザが規定することができる。ユーザは、流れデータが採取されるシミュレーションの間に使用されたものと同一の時間間隔Δτに従うことができる。これと代替的に、ユーザは、溶融体の動作に関する望ましい解法に従って時間増分を規定してもよい。
【０１７５】
段階６０３にて、部分的に充填された制御容積の各々における溶融体の前面の位置及び形状は、次式で規定される溶融体の容積部分ｆｍに基づいて数学的に識別される。
【０１７６】
【数４５】
2004506515.tif 000035
【０１７７】
ここで、Ｖｍは溶融体が占める制御容積の容積
Ｖは制御容積の全容積
【０１７８】
本発明の一実施形態において、溶融体の前面の自由面は、結節点点から補間された溶融体の容積部分が０．５すなわちｆｍ＝０．５である表面として規定される。溶融体の前面が存在する要素は、最初に、その結節点を容積部分から識別する。溶融体の前面を保持する１つの要素において、溶融体の容積割合は、いくつかの結節点に対し０．５以上であり、その他の結節点に対し０．５以下である。溶融体の容積割合が０．５に等しい点は、その後に計算される要素の端縁と交差する。それにより集められた点を使用してその中心心が自由面の位置ベクトルＸｆｓにより示される多面体を画成する。溶融体の容積割合が溶融体の前面を規定するために使用されるが、本発明の範囲から逸脱せずに、多数のその他の定義を使用することができる。
【０１７９】
シミュレーションの開始時、この方法は、溶融体の前面の自由面に沿って材料の粒子を形成する。これら粒子は、標識により視覚的に表現する数学的点である。標識は色又は記号或いはその組み合わせとすることができる。好ましくは、同時に且つ同一の溶融体の前面に形成される粒子には、同一の標識が割り当てられるようにする。最初に、粒子の各々には位置ベクトルＸｐ゜、温度Ｔｐ゜、速度ベクトルＶｐ゜も割り当てられる。粒子の初期速度及び温度は、結節点から１つの要素内で補間された関係する溶融体の全面の速度及び温度に等しい。粒子の各々の位置、温度及び速度が追跡される。
【０１８０】
流体の流れが進むのに伴い、溶融体の前面は拡張し又は収縮し、溶融体の前面に沿った粒子の密度を変化させる。段階６０４において、この方法は、溶融体の前面における粒子数が十分であるかどうかを決定する。この方法は、溶融体の前面に付着したものと考えられる粒子数をカウントする。
【０１８１】
この方法は、各粒子の位置ベクトルＸｐを自由面の平均位置Ｘｆｓと比較する。粒子がユーザが規定した自由面の点の許容公差範囲内にあるならば、その粒子は附着した粒子として規定される。粒子が規定された自由面の点の許容公差外にあるならば、その粒子は分離した粒子である。
【０１８２】
図示した実施の形態において、この方法は、判定基準として粒子と溶融体の前面との間の距離を使用する。溶融体の前面の自由面からの粒子間の距離がユーザが規定した値以内であるならば、粒子は付着している、すなわち、溶融体の前面にあるものと見なされる。距離の基準は、粒子が依然として溶融体の前面に付着しているかどうかを決定するために使用されるが、液体の割合の範囲、粒子密度等のようなその他の基準を本発明の精神及び範囲から逸脱せずに使用することができる。
【０１８３】
溶融体の前面における粒子数が、ユーザが特定した粒子数密度よりも少ないならば、モジュールは、段階６０５にて新たな粒子を自由面に追加する。この方法は、これら新たに追加された粒子を追跡するため異なる色又は記号を割り当てる。粒子の各々が形成される時点も記録される。この場合にも、粒子が一度び形成されたならば、それら粒子の各々の位置、温度及び速度が保存され且つ連続的に追加される。自由面における粒子数が十分であると判断されたならば、モジュールは段階６０６へと進む。
【０１８４】
段階６０６にて、全ての粒子の位置が更新される。時間ステップの各々にて、粒子は、それ以前の位置にて計算された溶融体の速度に基づいて新たな位置まで前進する。この新たな位置は、局所的な速度ベクトルを使用して更新される。
【０１８５】
【数４６】
2004506515.tif 000036
【０１８６】
ここで、Ｘは粒子の位置ベクトル
Ｖは粒子の位置における流体の速度ベクトル
Δτは時間ステップ
ｉ、ｉ＋１は時間ステップである。
【０１８７】
これらの移流過程の間、粒子がモデルの境界に進入しないように粒子は拘束される。その新たな位置が記憶される。
【０１８８】
段階６０７にて、要素の各々における自由面の平均位置が更新される。各時間ステップにて、溶融体の増分的正味流入量のため、溶融体の前面の位置をこれに応じて調整しなければならないから、各制御容積中、溶融体の容積割合は変化する。溶融体の前面の位置は、各時間ステップにて更新される。
【０１８９】
段階６０８にて、この方法は上述した手順に従って各粒子の付着／分離状況を更新する。粒子が最早、溶融体の前面に付着していないならば、最初に分離した粒子に割り当てられた記号を変更し、視覚的識別が容易であるように異なる色／記号に変える。
【０１９０】
段階６０９において、視覚的解析のため、粒子の位置がコンピュータスクリーンに表現される。付着した粒子及び分離した粒子は、色／記号で識別される。視覚的に、異なる溶融体の前面を表わす異なる記号及び色の粒子が衝突する箇所がメンドラインの位置である。更に修復の質及び強度を粒子の温度及び混合程度によって視覚的に評価することもできる。メンドラインにおける温度が高く且つ混合が完全であればある程、修復は益々、良好となる。
【０１９１】
段階６１０にて、この方法は停止時間を点検する。停止時間に達していないならば、この方法は段階６０２に戻り、異なる時間ステップにてシミュレーションを繰り返す。停止時間に達しているならば、シミュレーションは完了する。ユーザはそのシミュレーションの結果を出力することができる。
【図面の簡単な説明】
【図１】本発明を具体化する機械の概略図である。
【図２】本発明の方法の全体的な解答手順の一実施形態のフローチャートである。
【図３】本発明方法の２次元流体の流れをモデル化するのに適した三角形有限要素の概略図である。
【図４】２次元結節点制御容積を示す概略図である。
【図５】本発明方法の３次元流体の流れをモデル化する四面体有限要素の概略図である。
【図６】計算領域の境界を示す２次元的流れ領域の一例である。
【図７】本発明方法のショットスリーブ内の流れモジュールの一実施形態を示すフローチャートである。
【図８】本発明方法の収縮巣予測モジュールの一実施形態を示すフローチャートである。
【図９】金型表面に対する伝導線部分を示す熱伝導流体供給管の概略図である。
【図１０】本発明方法の熱伝導流体供給モジュールの一実施形態を示すフローチャートである。
【図１１】本発明方法のダイ潤滑剤冷却モジュールの一実施形態を示すフローチャート図である。
【図１２】本発明方法のメンドライン予測モジュールの一実施形態を示すフローチャート図である。
【図１】
2004506515.tif 000049
【図２】
2004506515.tif 000050
2004506515.tif 000051
【図６】
2004506515.tif 000052
【図７】
2004506515.tif 000053
【図８】
2004506515.tif 000054
【図１０】
2004506515.tif 000055
【図１１】
2004506515.tif 000056
【図１２】
2004506515.tif 000057
【国際公開パンフレット】
2004506515.tif 000058
2004506515.tif 000059
2004506515.tif 000060
2004506515.tif 000061
2004506515.tif 000062
2004506515.tif 000063
2004506515.tif 000064
2004506515.tif 000065
2004506515.tif 000066
2004506515.tif 000067
2004506515.tif 000068
2004506515.tif 000069
2004506515.tif 000070
2004506515.tif 000071
2004506515.tif 000072
2004506515.tif 000073
2004506515.tif 000074
2004506515.tif 000075
2004506515.tif 000076
2004506515.tif 000077
2004506515.tif 000078
2004506515.tif 000079
2004506515.tif 000080
2004506515.tif 000081
2004506515.tif 000082
2004506515.tif 000083
2004506515.tif 000084
2004506515.tif 000085
2004506515.tif 000086
2004506515.tif 000087
2004506515.tif 000088
2004506515.tif 000089
2004506515.tif 000090
2004506515.tif 000091
2004506515.tif 000092
2004506515.tif 000093
2004506515.tif 000094
2004506515.tif 000095
2004506515.tif 000096
2004506515.tif 000097
2004506515.tif 000098
2004506515.tif 000099
2004506515.tif 000100
2004506515.tif 000101
2004506515.tif 000102
2004506515.tif 000103
2004506515.tif 000104
2004506515.tif 000105
2004506515.tif 000106
2004506515.tif 000107
2004506515.tif 000108
2004506515.tif 000109
2004506515.tif 000110
2004506515.tif 000111
2004506515.tif 000112
2004506515.tif 000113
2004506515.tif 000114
2004506515.tif 000115
2004506515.tif 000116
2004506515.tif 000117
2004506515.tif 000118
2004506515.tif 000119
2004506515.tif 000120
2004506515.tif 000121
2004506515.tif 000122
2004506515.tif 000123
2004506515.tif 000124
2004506515.tif 000125
2004506515.tif 000126
2004506515.tif 000127
2004506515.tif 000128
2004506515.tif 000129
2004506515.tif 000130
2004506515.tif 000131
2004506515.tif 000132
【国際調査報告】
2004506515.tif 000133
2004506515.tif 000134
2004506515.tif 000135