JPH0415761A - 冷熱サイクル構造体温度の解析方法および金型装置系の設計装置 - Google Patents
冷熱サイクル構造体温度の解析方法および金型装置系の設計装置Info
- Publication number
- JPH0415761A JPH0415761A JP2116359A JP11635990A JPH0415761A JP H0415761 A JPH0415761 A JP H0415761A JP 2116359 A JP2116359 A JP 2116359A JP 11635990 A JP11635990 A JP 11635990A JP H0415761 A JPH0415761 A JP H0415761A
- Authority
- JP
- Japan
- Prior art keywords
- calculation
- heat flux
- cooling
- data
- mold
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims description 22
- 238000013461 design Methods 0.000 title claims description 19
- 238000001816 cooling Methods 0.000 claims abstract description 45
- 230000004907 flux Effects 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000004141 dimensional analysis Methods 0.000 claims abstract description 15
- 230000008859 change Effects 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 92
- 238000010438 heat treatment Methods 0.000 claims description 29
- 238000004904 shortening Methods 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 8
- 125000004122 cyclic group Chemical group 0.000 claims description 4
- 230000008569 process Effects 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 8
- 239000002826 coolant Substances 0.000 description 4
- 238000002347 injection Methods 0.000 description 4
- 239000007924 injection Substances 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000000465 moulding Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000005094 computer simulation Methods 0.000 description 3
- 238000012938 design process Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 239000000110 cooling liquid Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000011347 resin Substances 0.000 description 2
- 229920005989 resin Polymers 0.000 description 2
- 241000282326 Felis catus Species 0.000 description 1
- 239000003637 basic solution Substances 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Landscapes
- Moulds For Moulding Plastics Or The Like (AREA)
- Molds, Cores, And Manufacturing Methods Thereof (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔産業上の利用分野〕
本発明は、冷熱サイクル構造体温度の解析方法および金
型装置系の設計装置に関し、さらに詳しくは、例えば射
出成形用金型のように冷熱サイクルの繰り返しのもとで
使用される金型において、その金型の形状や金型に用い
る金属材料の種類、金型冷却管の配管状態や冷却液の流
量、さらには成形品の形状を規定することとなる金型キ
ャビティ形状、およびこれらの金型やその冷却系の運転
条件等、金型装置系全般の構成や運転条件の最適な設計
を行うための解析方法および設計装置に関する。
型装置系の設計装置に関し、さらに詳しくは、例えば射
出成形用金型のように冷熱サイクルの繰り返しのもとで
使用される金型において、その金型の形状や金型に用い
る金属材料の種類、金型冷却管の配管状態や冷却液の流
量、さらには成形品の形状を規定することとなる金型キ
ャビティ形状、およびこれらの金型やその冷却系の運転
条件等、金型装置系全般の構成や運転条件の最適な設計
を行うための解析方法および設計装置に関する。
一般に、金型装置系等の冷熱サイクル構造体の構成や運
転条件が最適に設計されていないと、成形効率か悪くな
ったり、成形品の不良を多く生じたりする。そしてこの
ような不具合を机上の設計段階で予測することは困難な
場合が多く、又、金型装置系を実際に試作して見てトラ
イアンドエラーで改良することも行われるか、コスト面
や効率面で有利ではない。
転条件が最適に設計されていないと、成形効率か悪くな
ったり、成形品の不良を多く生じたりする。そしてこの
ような不具合を机上の設計段階で予測することは困難な
場合が多く、又、金型装置系を実際に試作して見てトラ
イアンドエラーで改良することも行われるか、コスト面
や効率面で有利ではない。
このような点から、近時、計算機シミュレーションによ
る解析を利用して最適の設計を得ようとする手法が用い
られている(特開昭63−46511号公報、特開昭6
3−165970号公報等参照)。そしてこれらの従来
技術においては解析方法としていわゆる定常解析法や非
定常解析法が用いられ、また非定常解析法の計算手法と
して差分法、有限要素法、あるいは境界要素法が用いら
れている。
る解析を利用して最適の設計を得ようとする手法が用い
られている(特開昭63−46511号公報、特開昭6
3−165970号公報等参照)。そしてこれらの従来
技術においては解析方法としていわゆる定常解析法や非
定常解析法が用いられ、また非定常解析法の計算手法と
して差分法、有限要素法、あるいは境界要素法が用いら
れている。
ところで、例えば射出成形用金型のようにサイクリック
な冷熱プロセスの繰り返しのもとで使用される金型を、
その熱的条件面の考慮から最適に設計しようとする場合
、熱流束(単位時間(5ec)に単位面積(cn)を通
って出入りする熱量(cat))と金型温度とか刻々と
変化する冷熱サイクルを解析する為に非定常解析法が必
要となり、かつその計算手法としては、3次元解析モデ
ルの作成か著しく容易になる点から境界要素法が著しく
有利である。
な冷熱プロセスの繰り返しのもとで使用される金型を、
その熱的条件面の考慮から最適に設計しようとする場合
、熱流束(単位時間(5ec)に単位面積(cn)を通
って出入りする熱量(cat))と金型温度とか刻々と
変化する冷熱サイクルを解析する為に非定常解析法が必
要となり、かつその計算手法としては、3次元解析モデ
ルの作成か著しく容易になる点から境界要素法が著しく
有利である。
それにも拘らず、従来、上記のような設計目的で境界要
素法による非定常解析を行う技術は実用化されていなか
った。その理由は、金型温度を解析するためには、定数
として与えられる種々の条件の他に、設定時刻に従って
変化する熱流束を経時的な積分演算によって知る必要か
あり、例えばある設定時刻n+1での金型温度を知るた
めには設定時刻0〜nの全時点についての(n+1)回
の積分演算を行わなければならないという点にあった。
素法による非定常解析を行う技術は実用化されていなか
った。その理由は、金型温度を解析するためには、定数
として与えられる種々の条件の他に、設定時刻に従って
変化する熱流束を経時的な積分演算によって知る必要か
あり、例えばある設定時刻n+1での金型温度を知るた
めには設定時刻0〜nの全時点についての(n+1)回
の積分演算を行わなければならないという点にあった。
従って、ある程度以上のサイクル数に渉って金型の冷熱
サイクルを解析しようとする場合、上記設定時刻nの値
が太き(なるため、所要計算量が膨大なものとなり、通
常容量の計算機では計算時間か非常に長くなって、設計
工程の遅延か著しく、反面、大容量の計算機はコスト面
で用い難いとの理由から実用化を阻まれていたのである
。
サイクルを解析しようとする場合、上記設定時刻nの値
が太き(なるため、所要計算量が膨大なものとなり、通
常容量の計算機では計算時間か非常に長くなって、設計
工程の遅延か著しく、反面、大容量の計算機はコスト面
で用い難いとの理由から実用化を阻まれていたのである
。
C本願発明の説明〕
(発明の目的)
本願発明は、・上記の問題点に鑑み、計算機シミュレー
ションを利用した金型装置系等の冷熱サイクル構造体の
設計において、境界要素法による非定常解析を行うこと
の利点を維持したままで、しかも設計工程の短縮、コス
ト低減の面で実用性のある金型装置系の設計方法および
設計装置を提供することを目的とする。
ションを利用した金型装置系等の冷熱サイクル構造体の
設計において、境界要素法による非定常解析を行うこと
の利点を維持したままで、しかも設計工程の短縮、コス
ト低減の面で実用性のある金型装置系の設計方法および
設計装置を提供することを目的とする。
(第1発明および第2発明の構成)
第1発明の構成は、第1図のフローチャートに示すよう
に、CAE/CAD/CAMシステムを用いる解析方法
であって、冷熱サイクル構造体の3次元解析モデルに基
づき、その各冷熱サイクルの特定時刻間の熱流束の平均
値を用いて境界要素法により近似的に積分演算し、構造
体の温度変化を求めることを特徴とする冷熱サイクル構
造体温度の解析方法である。
に、CAE/CAD/CAMシステムを用いる解析方法
であって、冷熱サイクル構造体の3次元解析モデルに基
づき、その各冷熱サイクルの特定時刻間の熱流束の平均
値を用いて境界要素法により近似的に積分演算し、構造
体の温度変化を求めることを特徴とする冷熱サイクル構
造体温度の解析方法である。
第2発明の構成は、第2図のブロック図に示すように、
データの入力を行うデータ入力手段1と、前記データ入
力手段により入力された形状データに基づいて金型やそ
の冷却装置の3次元解析モデルを作成するモデリング手
段2と、前記データ入力学段1により入力された条件デ
ータを前記3次元解析モデルに付加して境界要素法によ
り金型のサイクリックな温度変化を設定時刻の経過に従
って順次積分演算する演算手段3と、前記演算手段3に
おける積分演算か各冷熱サイクルの終了時に該当する設
定時刻に至るごとにそれまでの全時刻における熱流束の
平均値を算出する平均熱流束算出手段4と、前記平均熱
流束算出手段4により算出された熱流束の平均値を前記
演算手段3に付加して、次の冷熱サイクルでの積分演算
において、それまでの全時刻に渉る積分演算の結果に1
ステップで代替する熱流束値として用いさせる演算ステ
ップ短縮手段5と、前記演算手段3、平均熱流束算出手
段4、演算ステップ短縮手段5による所要の冷熱サイク
ル数に渉る演算の結果を前記形状データあるいは出力デ
ータと対比させるべく出力する出力手段6とを有する金
型装置系の設計装置である。
データの入力を行うデータ入力手段1と、前記データ入
力手段により入力された形状データに基づいて金型やそ
の冷却装置の3次元解析モデルを作成するモデリング手
段2と、前記データ入力学段1により入力された条件デ
ータを前記3次元解析モデルに付加して境界要素法によ
り金型のサイクリックな温度変化を設定時刻の経過に従
って順次積分演算する演算手段3と、前記演算手段3に
おける積分演算か各冷熱サイクルの終了時に該当する設
定時刻に至るごとにそれまでの全時刻における熱流束の
平均値を算出する平均熱流束算出手段4と、前記平均熱
流束算出手段4により算出された熱流束の平均値を前記
演算手段3に付加して、次の冷熱サイクルでの積分演算
において、それまでの全時刻に渉る積分演算の結果に1
ステップで代替する熱流束値として用いさせる演算ステ
ップ短縮手段5と、前記演算手段3、平均熱流束算出手
段4、演算ステップ短縮手段5による所要の冷熱サイク
ル数に渉る演算の結果を前記形状データあるいは出力デ
ータと対比させるべく出力する出力手段6とを有する金
型装置系の設計装置である。
(第1発明及び第2発明の作用・効果)次に、第1発明
、第2発明の作用、効果を第1図のフローチャート、第
2図のブロック図に基づいて説明する。
、第2発明の作用、効果を第1図のフローチャート、第
2図のブロック図に基づいて説明する。
まず、冷熱サイクル構造体、例えば、金型やその冷却装
置の形状、構造を規定する形状データか、データ入力手
段1によって計算機に入力される。
置の形状、構造を規定する形状データか、データ入力手
段1によって計算機に入力される。
かかる形状データとして、金型の成形面や外側面の形状
、冷却管の内径及び外径や配管経路等のデータが挙げら
れる。
、冷却管の内径及び外径や配管経路等のデータが挙げら
れる。
次いでモデリング手段2により、上記形状データに基づ
いて金型やその冷却装置の3次元解析モデルか作成され
る。この3次元解析モデルは境界要素法に特有の境界モ
デル、サーフェイスモデル等と称されるものであり、差
分法や有限要素法において作成されるモデルに比し、そ
の作成か著しく容易である。
いて金型やその冷却装置の3次元解析モデルか作成され
る。この3次元解析モデルは境界要素法に特有の境界モ
デル、サーフェイスモデル等と称されるものであり、差
分法や有限要素法において作成されるモデルに比し、そ
の作成か著しく容易である。
一方、金型の冷熱サイクルを規定する条件データも前記
データ入力手段1により計算機に入力される。かかる条
件データとしては、金型構成材料の熱伝導度や比熱、被
成形材料の射出時温度や比熱、金型の冷却時間、冷却液
の流量等が挙げられる。
データ入力手段1により計算機に入力される。かかる条
件データとしては、金型構成材料の熱伝導度や比熱、被
成形材料の射出時温度や比熱、金型の冷却時間、冷却液
の流量等が挙げられる。
こうして、計算機に入力された条件データか前記3次元
解析モデルに付加され、演算手段3によって積分演算さ
れる。この積分演算は、初期条件として与えられる条件
データを境界要素法に従って所定の熱伝導方程式にあて
はめ、設定時刻の経過に伴う熱流束の変化と、これに基
づく金型温度の変化とを経時的に算出するものである。
解析モデルに付加され、演算手段3によって積分演算さ
れる。この積分演算は、初期条件として与えられる条件
データを境界要素法に従って所定の熱伝導方程式にあて
はめ、設定時刻の経過に伴う熱流束の変化と、これに基
づく金型温度の変化とを経時的に算出するものである。
但し、上記の通りの演算が行われるのは、第1の冷熱サ
イクルだけである。即ち、第1の冷熱サイクルの終了時
に該当する設定時刻に到ると、平均熱流束算出手段4に
よってそれまでの全時刻における熱流束の平均値(以下
、これを「平均熱流束」という。)か算出され、次いで
演算ステップ短縮手段5によってこの平均熱流束か前記
演算手段3に付加される。
イクルだけである。即ち、第1の冷熱サイクルの終了時
に該当する設定時刻に到ると、平均熱流束算出手段4に
よってそれまでの全時刻における熱流束の平均値(以下
、これを「平均熱流束」という。)か算出され、次いで
演算ステップ短縮手段5によってこの平均熱流束か前記
演算手段3に付加される。
この平均熱流束は、第2の冷熱サイクルでの積分演算の
基礎となるへき、第1の冷熱サイクルの全時刻に渉る積
分演算の結果としての熱流束の値に1ステップで代替す
るものである。従って、例えば第3図に示すように、第
1の冷熱サイクルにn点の設定時刻かあった場合、第2
の冷熱サイクルにおける設定時刻n+1での所要計算量
は、本来ならば設定時刻0〜nの全時刻にわたる(n+
1)回の累積的な積分演算を要するところが、前記平均
熱流束を用いる1回の積分演算で済むこととなる。しか
も平均熱流束はそれまでの全時刻における熱流束の平均
値として与えられるので、実際の熱流束の値と十分に近
似している。
基礎となるへき、第1の冷熱サイクルの全時刻に渉る積
分演算の結果としての熱流束の値に1ステップで代替す
るものである。従って、例えば第3図に示すように、第
1の冷熱サイクルにn点の設定時刻かあった場合、第2
の冷熱サイクルにおける設定時刻n+1での所要計算量
は、本来ならば設定時刻0〜nの全時刻にわたる(n+
1)回の累積的な積分演算を要するところが、前記平均
熱流束を用いる1回の積分演算で済むこととなる。しか
も平均熱流束はそれまでの全時刻における熱流束の平均
値として与えられるので、実際の熱流束の値と十分に近
似している。
そして、平均熱流束算出手段4と演算ステップ短縮手段
5とによる上記のような演算ステップ短縮処理か各冷熱
サイクルの終了時ごとに行われるので、特に多数の冷熱
サイクルに渉る経時的積分演算を行う場合には、解析結
果の信頼性を維持したままで、その演算ステップ短縮効
果、ひいては計算時間の短縮効果はきわめて大きなもの
となる。
5とによる上記のような演算ステップ短縮処理か各冷熱
サイクルの終了時ごとに行われるので、特に多数の冷熱
サイクルに渉る経時的積分演算を行う場合には、解析結
果の信頼性を維持したままで、その演算ステップ短縮効
果、ひいては計算時間の短縮効果はきわめて大きなもの
となる。
こうして、演算手段3、平均熱流束算出手段4、演算ス
テップ短縮手段5による所要の冷熱サイクル数に渉る演
算か行われ、金型の温度変化が算出される。この演算結
果は出力手段6によって数値化データ、図形化データ等
として出力されるが、非定常解析を行っているため、金
型の任意の部位における経時的な温度変化を示すデータ
や、任意の設定時刻における金型全体の温度分布を示す
データ等、多様なデータの出力が可能である。従って、
金型の冷熱サイクルを多面的かつ有効に解析することが
できる。
テップ短縮手段5による所要の冷熱サイクル数に渉る演
算か行われ、金型の温度変化が算出される。この演算結
果は出力手段6によって数値化データ、図形化データ等
として出力されるが、非定常解析を行っているため、金
型の任意の部位における経時的な温度変化を示すデータ
や、任意の設定時刻における金型全体の温度分布を示す
データ等、多様なデータの出力が可能である。従って、
金型の冷熱サイクルを多面的かつ有効に解析することが
できる。
そして、これらの出力データに基づき、当初に入力した
形状データや条件データが評価、修正されて、金型装置
系の最適設計が行われるのである。
形状データや条件データが評価、修正されて、金型装置
系の最適設計が行われるのである。
以上の点から、本願第1発明および第2発明によれば、
計算機シミュレーションを利用した金型装置系の設計に
おいて、境界要素法による非定常解析を行うことの利点
を維持したままで、通常容量の計算機によっても設計工
程の短縮を十分に図ることができる。
計算機シミュレーションを利用した金型装置系の設計に
おいて、境界要素法による非定常解析を行うことの利点
を維持したままで、通常容量の計算機によっても設計工
程の短縮を十分に図ることができる。
(第1発明及び第2発明のその他の発明の説明)第1発
明の「冷熱サイクル構造体」とは、サイクリックな冷熱
変化を伴う一定の構造体を言い、射出成形用金型がその
一例として挙げられるが、これに限定されない。
明の「冷熱サイクル構造体」とは、サイクリックな冷熱
変化を伴う一定の構造体を言い、射出成形用金型がその
一例として挙げられるが、これに限定されない。
前記第1発明および第2発明において、データ入力手段
の種類は限定されず、入力用キーボード、平面図形の読
み取り装置、立体図形の認識装置等の任意の手段により
データを入力することができる。
の種類は限定されず、入力用キーボード、平面図形の読
み取り装置、立体図形の認識装置等の任意の手段により
データを入力することができる。
条件データはモデリング手段2による3次元解析モデル
の作成後に計算機に入力しても良いし、あるいは3次元
解析モデルの作成前(例えば、形状データの入力と同時
あるいはそれ以前)に計算機に入力して、3次元解析モ
デルに付加されるまでの間、読み出し可能な状態で計算
機に記憶させておいても良い。
の作成後に計算機に入力しても良いし、あるいは3次元
解析モデルの作成前(例えば、形状データの入力と同時
あるいはそれ以前)に計算機に入力して、3次元解析モ
デルに付加されるまでの間、読み出し可能な状態で計算
機に記憶させておいても良い。
演算手段3での積分演算は、非定常な熱伝導現象を支配
する適切な方程式に従って行われる。このような方程式
の一例として、後に述べる実施例で記すものか挙げられ
るが、これに限定されるものではなく、設定時刻の経過
に従ってサイクリックな温度変化の経時的な積分演算を
行うための他の適切な方程式であっても用い得る。
する適切な方程式に従って行われる。このような方程式
の一例として、後に述べる実施例で記すものか挙げられ
るが、これに限定されるものではなく、設定時刻の経過
に従ってサイクリックな温度変化の経時的な積分演算を
行うための他の適切な方程式であっても用い得る。
出力手段6により出力された演算結果は、これを単に人
為的に評価して形状データや条件データを修正するとい
う方式の他、金型の温度変化や温度分布等の望ましいス
タンダードを予め計算機に記憶させておき、計算機内に
設けた評価手段によって前記演算結果をスタンダードと
の対比で評価させ、その結果から初期の形状データや条
件データか自動的に修正を受けるようにしても良い。更
に、このように修正された形状データや条件データに基
づいて3次元解析モデルやこれに付加される条件データ
が自動的に再構成され、再度積分演算が実施されるよう
に構成しても良い。この場合、上記プロセスが繰り返さ
れるごとにスタンダードとの対比で最適設計に至ったか
否かを評価手段によって判定させ、YESの判定信号に
よってのみ演算の停止と出力手段6による演算結果の出
力とか行われるように構成して、金型装置系の最適設計
に到る工程を全自動化することもできる。この場合、工
程短縮化のメリットがとりわけ大きい。
為的に評価して形状データや条件データを修正するとい
う方式の他、金型の温度変化や温度分布等の望ましいス
タンダードを予め計算機に記憶させておき、計算機内に
設けた評価手段によって前記演算結果をスタンダードと
の対比で評価させ、その結果から初期の形状データや条
件データか自動的に修正を受けるようにしても良い。更
に、このように修正された形状データや条件データに基
づいて3次元解析モデルやこれに付加される条件データ
が自動的に再構成され、再度積分演算が実施されるよう
に構成しても良い。この場合、上記プロセスが繰り返さ
れるごとにスタンダードとの対比で最適設計に至ったか
否かを評価手段によって判定させ、YESの判定信号に
よってのみ演算の停止と出力手段6による演算結果の出
力とか行われるように構成して、金型装置系の最適設計
に到る工程を全自動化することもできる。この場合、工
程短縮化のメリットがとりわけ大きい。
次に、本願第1発明および第2発明の一実施例を説明す
る。本実施例は、第4図に示すような固定型7と可動型
8とからなる箱型成形品の射出成形用金型において、両
型7.8の内部にそれぞれ一点鎖線で示す冷却管9、l
Oが配置された場合の各冷却管9.10の最適冷却液流
量を求めるために行ったものである。
る。本実施例は、第4図に示すような固定型7と可動型
8とからなる箱型成形品の射出成形用金型において、両
型7.8の内部にそれぞれ一点鎖線で示す冷却管9、l
Oが配置された場合の各冷却管9.10の最適冷却液流
量を求めるために行ったものである。
そして本実施例の設計装置は、第5図に示すように、メ
インコンピュータ11、グラフィックコンピュータ12
、モデル表示装置t13、出力表示装置14より構成さ
れる。
インコンピュータ11、グラフィックコンピュータ12
、モデル表示装置t13、出力表示装置14より構成さ
れる。
メインコンピュータ11には、前記第2図に示したデー
タ入力手段1、演算手段3、平均熱流束算出手段4、演
算ステップ短縮手段5か含まれ、グラフィックコンピュ
ータ12にはモデリング手段2が含まれると共に、3次
元解析モデルかモデル表示装fl13で表示てきるよう
になっており、出力表示装置14は前記出力手段6に相
当す、るものである。
タ入力手段1、演算手段3、平均熱流束算出手段4、演
算ステップ短縮手段5か含まれ、グラフィックコンピュ
ータ12にはモデリング手段2が含まれると共に、3次
元解析モデルかモデル表示装fl13で表示てきるよう
になっており、出力表示装置14は前記出力手段6に相
当す、るものである。
本実施例の具体的内容を第6図に基づいて説明すると、
まず、データ入力手段1の条件データ入力回路15に以
下の■〜■の条件データを入力する。
まず、データ入力手段1の条件データ入力回路15に以
下の■〜■の条件データを入力する。
■ 成形条件(冷却時間、成形時間、樹脂温度、大気温
度、冷却液温度・流量) ■ 樹脂熱特性データ (熱伝導率、密度、比熱、固化温度、 流動停止温度、潜熱) ■ 金型熱特性データ (熱伝導率、密度、比熱) ■ 冷却液熱・流動特性データ (熱伝導率、密度、比熱、粘度) これらの条件データは、条件データメモリ16に保存さ
れる。
度、冷却液温度・流量) ■ 樹脂熱特性データ (熱伝導率、密度、比熱、固化温度、 流動停止温度、潜熱) ■ 金型熱特性データ (熱伝導率、密度、比熱) ■ 冷却液熱・流動特性データ (熱伝導率、密度、比熱、粘度) これらの条件データは、条件データメモリ16に保存さ
れる。
次に、データ入力手段1の形状データ入力回路17に固
定型7、可動型8、およびそれらの冷却管9.10につ
いての3次元形状のデータか入力され、これらのデータ
は形状データメモリ18に保存される。
定型7、可動型8、およびそれらの冷却管9.10につ
いての3次元形状のデータか入力され、これらのデータ
は形状データメモリ18に保存される。
次いて形状データメモリ18に保存された形状データの
読み込みによってモデリング手段2のモデル作成回路1
9か第7図のような三次元解析モデルを作成し、前記モ
デル表示装置13にデイスプレィすると共に、この3次
元解析モデルを演算手段3のモデル入力回路20によっ
て積分演算回路21に入力する。
読み込みによってモデリング手段2のモデル作成回路1
9か第7図のような三次元解析モデルを作成し、前記モ
デル表示装置13にデイスプレィすると共に、この3次
元解析モデルを演算手段3のモデル入力回路20によっ
て積分演算回路21に入力する。
一方、前記条件データメモリ16に保存された条件デー
タも、条件データ付加回路22によって読み込まれ、積
分演算回路21に付加される。
タも、条件データ付加回路22によって読み込まれ、積
分演算回路21に付加される。
ここで、積分演算回路21が、3次元解析モデルと、こ
れに初期条件として付加された条件データとから演算を
行うにあたり、以下の熱伝導方程式を用いる。
れに初期条件として付加された条件データとから演算を
行うにあたり、以下の熱伝導方程式を用いる。
(1)式において、Tは温度、tは時間、ρは密度、C
pは比熱、λは熱伝導率、2はラプラシアンである。
pは比熱、λは熱伝導率、2はラプラシアンである。
上記の(1)式に適宜な初期条件と境界条件を与えて時
間依存性の基本解を用い、時刻nにおける積分方程式を
求め、それをマトリックス表示すると、次の(2)式か
得られる。
間依存性の基本解を用い、時刻nにおける積分方程式を
求め、それをマトリックス表示すると、次の(2)式か
得られる。
(2)式において、Sは設定時刻、Hは温度マトリック
ス、Gは熱流束マトリックス、Qは熱流束、〒。は初期
温度、Bは初期温度マトリックスである。
ス、Gは熱流束マトリックス、Qは熱流束、〒。は初期
温度、Bは初期温度マトリックスである。
上記(2)式によって熱流束の変化に基づき、積分演算
回路21か演算を行い、設定時刻の経過に伴う金型温度
のサイクリックな変化を求めるのである。
回路21か演算を行い、設定時刻の経過に伴う金型温度
のサイクリックな変化を求めるのである。
ところで、(2)式によれば、各設定時刻に於ける金型
温度を求めるには、累積的に積分演算を行う必要がある
ため、特に設定時刻の値がある程度以上に大きくなると
、積分演算の所要計算量か膨大なものとなる。そこで、
次のような演算ステップ短縮処理が行われる。
温度を求めるには、累積的に積分演算を行う必要がある
ため、特に設定時刻の値がある程度以上に大きくなると
、積分演算の所要計算量か膨大なものとなる。そこで、
次のような演算ステップ短縮処理が行われる。
即ち、第1の冷熱サイクルの終了時に該当する設定時刻
に至ると、平均熱流束算出手段4の設定時刻判定回路2
3かこのことを検出し、平均熱流束算出回路24を作動
させる。そして、平均熱流束算出回路24は平均熱流束
、即ちそれまでの全時刻における熱流束の平均値を算出
して、演算ステップ短縮手段5の平均熱流束付加回路2
5を通じ、この平均熱流束の値を積分演算回路21に付
加する。更に、初期値変更回路26か作動して、積分演
算回路21における第2の冷熱サイクルの積分演算にあ
たり、本来の初期条件に基づく設定時刻0からの積分演
算に代え、平均熱流束を用いてIステップの演算で代替
する。
に至ると、平均熱流束算出手段4の設定時刻判定回路2
3かこのことを検出し、平均熱流束算出回路24を作動
させる。そして、平均熱流束算出回路24は平均熱流束
、即ちそれまでの全時刻における熱流束の平均値を算出
して、演算ステップ短縮手段5の平均熱流束付加回路2
5を通じ、この平均熱流束の値を積分演算回路21に付
加する。更に、初期値変更回路26か作動して、積分演
算回路21における第2の冷熱サイクルの積分演算にあ
たり、本来の初期条件に基づく設定時刻0からの積分演
算に代え、平均熱流束を用いてIステップの演算で代替
する。
以後、各冷熱サイクルの終了時ごとにこのような演算ス
テップ短縮処理か行われる。
テップ短縮処理か行われる。
なお、設定時刻の進行に伴う積分演算の進行は、常に演
算終了判定回路27によりチエツクされ、予め解析の終
了時点として定められた設定時刻に至った時は、積分演
算が停止され、出力手段6のデータ出力回路28を介し
てデータ出力装置29に所望のデータ形態の解析データ
か出力される。
算終了判定回路27によりチエツクされ、予め解析の終
了時点として定められた設定時刻に至った時は、積分演
算が停止され、出力手段6のデータ出力回路28を介し
てデータ出力装置29に所望のデータ形態の解析データ
か出力される。
第8図は本実施例における固定型7および可動型8の経
時的な金型温度の変化の解析結果を示し、点線が計算値
、実線が金型試作品に熱電対を取り付けて実測した実験
値を示す。図に示すように、計算値は実測値と良好な対
応を示している。第8図の結果から、可動型8は固定型
7に比して常に高い温度域にあり、可動型8の冷却を強
化する必要のあることが判明したので、可動型8の冷却
管10における冷却液流量を増加させるように運転条件
の設計を修正した。
時的な金型温度の変化の解析結果を示し、点線が計算値
、実線が金型試作品に熱電対を取り付けて実測した実験
値を示す。図に示すように、計算値は実測値と良好な対
応を示している。第8図の結果から、可動型8は固定型
7に比して常に高い温度域にあり、可動型8の冷却を強
化する必要のあることが判明したので、可動型8の冷却
管10における冷却液流量を増加させるように運転条件
の設計を修正した。
第9図(a)、(b)はそれぞれ第5、第18の冷熱サ
イクルにおける可動型8の表面温度の分布を等混線で示
した解析結果である。これらの図によれば、可動型8の
表面温度は、特に上側の突出した角部30において高い
ことが認められ、この角部30を有効に冷却できるよう
に冷却管10の配管の設計を修正した。
イクルにおける可動型8の表面温度の分布を等混線で示
した解析結果である。これらの図によれば、可動型8の
表面温度は、特に上側の突出した角部30において高い
ことが認められ、この角部30を有効に冷却できるよう
に冷却管10の配管の設計を修正した。
なお、本実施例における第18の冷熱サイクルまでの解
析に要した計算時間は、同様の解析を演算ステップ短縮
処理なして行った場合に比し、約1/30てあった。
析に要した計算時間は、同様の解析を演算ステップ短縮
処理なして行った場合に比し、約1/30てあった。
第1図は本願第1発明を説明するフローチャート、第2
図は本願第2発明のクレーム対応ブロック図、第3図は
演算ステップ短縮処理を説明する図、第4図は実施例に
おいて解析対象とした金型とその冷却系とを示す斜視図
、第5図は本願第2発明のシステム構成図、第6図は実
施例のブロック図、第7図は実施例において作成された
3次元解析モデルを示す図、第8図は金型温度の経時的
変化の解析データを示す図、第9図(a)、(b)はそ
れぞれ、第5及び第18の冷熱サイクルにおける可動型
の表面温度の分布の解析データを示す図である。 データ入力手段 モデリング手段 演算手段 平均熱流束算出手段 演算ステップ短縮手段 出力手段 固定型 可動型 冷却管 冷却管
図は本願第2発明のクレーム対応ブロック図、第3図は
演算ステップ短縮処理を説明する図、第4図は実施例に
おいて解析対象とした金型とその冷却系とを示す斜視図
、第5図は本願第2発明のシステム構成図、第6図は実
施例のブロック図、第7図は実施例において作成された
3次元解析モデルを示す図、第8図は金型温度の経時的
変化の解析データを示す図、第9図(a)、(b)はそ
れぞれ、第5及び第18の冷熱サイクルにおける可動型
の表面温度の分布の解析データを示す図である。 データ入力手段 モデリング手段 演算手段 平均熱流束算出手段 演算ステップ短縮手段 出力手段 固定型 可動型 冷却管 冷却管
Claims (2)
- (1)CAE/CAD/CAMシステムを用いる解析方
法であって、 冷熱サイクル構造体の3次元解析モデルに基づき、その
各冷熱サイクルの特定時刻間の熱流束の平均値を用いて
境界要素法により近似的に積分演算し、 構造体の温度変化を求める、 ことを特徴とする冷熱サイクル構造体温度の解析方法。 - (2)データの入力を行うデータ入力手段と、前記デー
タ入力手段により入力された形状データに基づいて金型
やその冷却装置の3次元解析モデルを作成するモデリン
グ手段と、 前記データ入力手段により入力された条件データを前記
3次元解析モデルに付加して、境界要素法により金型の
サイクリックな温度変化を設定時刻の経過に従って順次
積分演算する演算手段と、前記演算手段における積分演
算が各冷熱サイクルの終了時に該当する設定時刻に至る
ごとにそれまでの全時刻における熱流束の平均値を算出
する平均熱流束算出手段と、 前記平均熱流束算出手段により算出された熱流束の平均
値を前記演算手段に付加して、次の冷熱サイクルでの積
分演算においてそれまでの全時刻に渉る積分演算の結果
に1ステップで代替する熱流束値として用いさせる演算
ステップ短縮手段と、前記演算手段、平均熱流束算出手
段、演算ステップ短縮手段による所要の冷熱サイクル数
に渉る演算の結果を前記形状データあるいは出力データ
と対比させるべく出力する出力手段と、 を有することを特徴とする金型装置系の設計装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2116359A JP3023969B2 (ja) | 1990-05-02 | 1990-05-02 | 冷熱サイクル構造体温度の解析方法および金型装置系の設計装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2116359A JP3023969B2 (ja) | 1990-05-02 | 1990-05-02 | 冷熱サイクル構造体温度の解析方法および金型装置系の設計装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH0415761A true JPH0415761A (ja) | 1992-01-21 |
JP3023969B2 JP3023969B2 (ja) | 2000-03-21 |
Family
ID=14685011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2116359A Expired - Fee Related JP3023969B2 (ja) | 1990-05-02 | 1990-05-02 | 冷熱サイクル構造体温度の解析方法および金型装置系の設計装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3023969B2 (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0787547A1 (en) * | 1996-01-31 | 1997-08-06 | ROLLS-ROYCE plc | A method of investment casting and a method of making an investment casting mould |
JP2001255078A (ja) * | 2000-03-13 | 2001-09-21 | Tokyo Gas Co Ltd | 熱交換体の形状評価方法 |
JP2008055488A (ja) * | 2006-09-01 | 2008-03-13 | Fuji Heavy Ind Ltd | 金型温度のシミュレーション方法及びシミュレーションプログラム |
WO2011067644A1 (en) | 2009-12-04 | 2011-06-09 | Toyota Jidosha Kabushiki Kaisha | Numerical analysis method of a vehicle drive |
-
1990
- 1990-05-02 JP JP2116359A patent/JP3023969B2/ja not_active Expired - Fee Related
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0787547A1 (en) * | 1996-01-31 | 1997-08-06 | ROLLS-ROYCE plc | A method of investment casting and a method of making an investment casting mould |
US5868194A (en) * | 1996-01-31 | 1999-02-09 | Rolls-Royce Plc | Method of investment casting and a method of making an investment casting mould |
JP2001255078A (ja) * | 2000-03-13 | 2001-09-21 | Tokyo Gas Co Ltd | 熱交換体の形状評価方法 |
JP4498527B2 (ja) * | 2000-03-13 | 2010-07-07 | 東京瓦斯株式会社 | 熱交換体における氷着形状解析方法及びそれを用いた熱交換体の形状最適化方法 |
JP2008055488A (ja) * | 2006-09-01 | 2008-03-13 | Fuji Heavy Ind Ltd | 金型温度のシミュレーション方法及びシミュレーションプログラム |
WO2011067644A1 (en) | 2009-12-04 | 2011-06-09 | Toyota Jidosha Kabushiki Kaisha | Numerical analysis method of a vehicle drive |
DE112010004655T5 (de) | 2009-12-04 | 2012-10-31 | Toyota Jidosha K.K. | Numerisches analysierverfahren einer fahrzeugantriebseinheit |
US8688315B2 (en) | 2009-12-04 | 2014-04-01 | Toyota Jidosha Kabushiki Kaisha | Numerical analysis method of a vehicle drive |
DE112010004655B4 (de) * | 2009-12-04 | 2016-09-01 | Toyota Jidosha Kabushiki Kaisha | Numerisches Analysierverfahren einer Fahrzeugantriebseinheit |
Also Published As
Publication number | Publication date |
---|---|
JP3023969B2 (ja) | 2000-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6816820B1 (en) | Method and apparatus for modeling injection of a fluid in a mold cavity | |
US5202843A (en) | CAE system for automatic heat analysis of electronic equipment | |
Shahane et al. | Finite volume simulation framework for die casting with uncertainty quantification | |
Özek et al. | Calculating molding parameters in plastic injection molds with ANN and developing software | |
CN103970928A (zh) | 模拟方法和模拟装置 | |
Hanoglu et al. | Thermo-mechanical analysis of hot shape rolling of steel by a meshless method | |
Sulaiman et al. | Modeling of the thermal history of the sand casting process | |
Pao et al. | A medial-axes-based interpolation method for solidification simulation | |
Navah et al. | Thermofluid topology optimization for cooling channel design | |
JPH0415761A (ja) | 冷熱サイクル構造体温度の解析方法および金型装置系の設計装置 | |
JP4032848B2 (ja) | 成形シミュレーション方法、成形シミュレーション装置及び成形シミュレーションプログラム並びに当該成形シミュレーションプログラムを記録したコンピュータ読みとり可能な記録媒体 | |
Zhou et al. | Computer filling simulation of injection molding based on 3D surface model | |
Babaei et al. | Improved advection algorithm of computational modeling of free surface flow using structured grids | |
JP4032755B2 (ja) | 成型シミュレーション方法、成型シミュレーション装置及び成型シミュレーションプログラム並びに当該成型シミュレーションプログラムを記録したコンピュータ読みとり可能な記録媒体 | |
Dantzig | Solidification modeling: Status and outlook | |
JP4052006B2 (ja) | 成型シミュレーション方法、成型シミュレーション装置及び成型シミュレーションプログラム並びに当該成型シミュレーションプログラムを記録したコンピュータ読みとり可能な記録媒体 | |
JP2010030259A (ja) | 成形品形状の予測方法、キャビティ形状の予測方法及びそれらの装置、成形品形状又はキャビティ形状の予測プログラムとその記憶媒体 | |
JP5889077B2 (ja) | 成形品収縮変形予測装置、成形品収縮変形予測方法及び成形品収縮変形予測プログラム | |
JP5235573B2 (ja) | 強度解析方法、強度解析装置及び強度解析プログラム | |
US20190197209A1 (en) | Anti-warping design method for resin molded article, recording medium, and anti-warping design device for resin molded article | |
JP4373353B2 (ja) | ダイカストシミュレーション方法、その装置及びそのプログラム、並びに当該プログラムを記録した記録媒体 | |
Erber et al. | Predicting the local solidification time using spherical neural networks | |
JPH10278085A (ja) | 射出成形プロセスにおける温度履歴予測装置及び方法 | |
Gebelin et al. | Modelling of the wax injection process for the investment casting process: Prediction of defects | |
Zhou et al. | Mould cooling simulation for injection moulding using a fast boundary element method approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |