JP4232707B2 - ダイカストシミュレーション方法、その装置及びそのプログラム、並びに当該プログラムを記録した記録媒体 - Google Patents
ダイカストシミュレーション方法、その装置及びそのプログラム、並びに当該プログラムを記録した記録媒体 Download PDFInfo
- Publication number
- JP4232707B2 JP4232707B2 JP2004213356A JP2004213356A JP4232707B2 JP 4232707 B2 JP4232707 B2 JP 4232707B2 JP 2004213356 A JP2004213356 A JP 2004213356A JP 2004213356 A JP2004213356 A JP 2004213356A JP 4232707 B2 JP4232707 B2 JP 4232707B2
- Authority
- JP
- Japan
- Prior art keywords
- casting
- analysis step
- stress
- solid phase
- die
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Description
ダイカストにより得られる鋳物における課題としては、ひけ巣欠陥、空気等の巻き込み欠陥、湯回り不良欠陥等の他に、鋳物に残る残留歪みおよび残留応力がある。
この要因は、鋳造圧力は、ひけ巣部への溶湯の補給への駆動力に使われるのみであり、固相の変化には寄与しないと考えられていたことによる。
しかしながら、このように、鋳造圧力が固相間の溶湯流れのみに寄与し、残留歪みは冷却時(温度300℃以下)の線膨張差により生じるという手法では、精度よく鋳物の残留歪みおよび残留応力を求めることができなかった。
上記鋳型の形状データを基に、鋳型部分に位置する鋳型要素、鋳型内の空間に位置する鋳物要素、及び該鋳物要素に加圧力を加える部分に位置する加圧要素として、それぞれ複数の微小要素に分割する要素作成ステップと、
上記鋳物要素のすべてに溶湯が充填された状態の熱的特性値を付与し、該熱的特性値に基づいて微小時間後の上記鋳物要素の温度と固相率を経時的に求める凝固解析ステップと、
該凝固解析ステップにより得られた温度及び/又は固相率を基にして上記鋳物要素の機械的特性値を求めると共に、上記加圧要素によって上記鋳物要素に加圧力を付加し、該加圧力と上記機械的特性値とを基にして、上記鋳物要素の歪み、変位及び応力を経時的に求める途中応力解析ステップと、
上記鋳物要素のすべてが固相になった後に、熱的特性値に基づいて微小時間後の上記鋳物要素の温度を経時的に求める冷却解析ステップと、
該冷却解析ステップにより得られた温度を基にして上記鋳物要素の機械的特性値を求めると共に、該機械的特性値を基にして、かつ、上記途中応力解析ステップにおいて得られた上記鋳物要素の歪み、変位および応力を初期値として用い、上記鋳物要素の歪み、変位及び応力を経時的に求める熱応力解析ステップとを有することを特徴とするダイカストシミュレーション方法にある(請求項1)。
それ故、従来のような凝固途中における固相の変形がないという前提の基に行うシミュレーションに比べ、鋳物に残る残留歪みおよび残留応力の解析の精度を格段に高めることができる。そして、このダイカストシミュレーション方法を用いれば、残留歪みおよび残留応力の小さい鋳物が得られる鋳物設計を支援することができる。
第3の発明は、上記第2の発明のダイカストシミュレーションプログラムを記録したコンピュータ読み取り可能な記憶媒体にある(請求項7)。
第4の発明は、上記第1の発明のダイカストシミュレーション方法を実行するよう構成されたコンピュータを有することを特徴とするダイカストシミュレーション装置にある(請求項8)。
要素作成ステップでは、取り込む鋳型の形状データとして様々な形式のデータを用いることができるが、例えば、国際規格IGES(Initial Graphics Exchange Specification)の機能を介して取り込むことができる。
なお、鋳型内の空間部分(キャビティ部分)はできる限り正確な形状が必要であるが、鋳型部分の形状としては、最小限に絞ることができる。また、上記加圧要素としては、例えば、実際の鋳造時おいて溶湯に圧力をかけるプランジャーに接する鋳物部分を想定して定義することができる。なお、この要素作成ステップにおいて、上記各微小要素を定義すると共に、それぞれの物性値等を与えることもできるが、これは後の解析ステップにおいて行ってもよい。
また、上記微小要素の大きさはできる限り小さい方が精度向上のために好ましいが、小さいほど解析時間が長くなるので、適度な大きさにすることが好ましい。また、微小要素の大きさをすべての部分で統一する必要はなく、部位によって変更することができる。
凝固解析ステップは、上記鋳物要素のすべてに溶湯が充填された状態から開始する。そして、このステップにより解析を開始する際に、各鋳物要素の熱的特性値の初期値を付与する。ここでいう熱的特性値としては、例えば、密度、比熱、熱伝導率、潜熱等がある。また、これらの特性値は、温度等によって変化させることが好ましいが、大幅な精度低下につながらない限り一定値に固定することもできる。
なお、この凝固解析ステップを実施する前に、後述する充填解析ステップを実施して、溶湯の挙動を考慮して溶湯充填完了時あるいはその後の冷却を経て凝固開始する際における各鋳物要素の熱的特性値の初期値を定めることもできる。
この凝固解析ステップは、鋳物要素のすべてが固相のみとなるまで、すなわち、すべての固相率が1になるまで繰り返し行われる。
例えば、鋳物要素と鋳型要素のそれぞれの比熱、熱伝導率、および密度と、鋳物要素と鋳型要素の間の熱伝達率等を用いた熱力学的エネルギの式を基にして熱伝導解析を行って温度を求め、その後、溶湯の潜熱放出を考慮して固相率を求める手法がある。
また、上記の鋳物要素の比熱をみかけ比熱(潜熱放出を考慮した比熱)に代えて熱伝導解析を行って温度を求め、その後、シャイルの式を用いて固相率を求める手法もある。
途中応力解析ステップでは、まず、上記凝固解析ステップにより温度及び固相率が得られる毎に、これらを基にして上記鋳物要素の機械的特性値を求める。ここでいう機械的特性値は、例えば、ヤング率、ポアソン比、線膨張係数、耐力、加工硬化指数等がある。なお、これらの機械的特性値は、上記温度あるいは固相率に基づいて必ずしも変化させる必要はなく、固定値を用いても誤差が少ない場合には固定値を適用することも可能である。
そして、上記加圧力と上記機械的特性値を基にして、上記鋳物要素の歪み、変位及び応力を求める。これを上記凝固解析ステップの結果が出る度に繰り返し行い、鋳物要素のすべてが固相、すなわちすべての固相率が1となるまで行う。
冷却解析ステップは、上記鋳物要素のすべてが固相(固相率が1)になった後、つまり、凝固が完了した後に実施し、上記熱的特性値に基づいて微小時間後の上記鋳物要素の温度を経時的に求める。この段階では、固相を求める必要がないので、少なくとも鋳物要素の温度の変化を求める。
熱応力解析ステップでは、上記途中応力解析ステップによって得られた各鋳物要素の歪み、変位および応力を初期値として、凝固完了後の歪み、変位および応力を引き続き求めていく。
すなわち、熱応力解析ステップでは、まず上記冷却解析ステップにより得られた温度を基にして上記鋳物要素の機械的特性値を求める。ここで用いる機械的特性値の種類としては、上記途中応力解析ステップにおいて用いたものと同様である。また、ここで、上記機械的特性値は、必ずしも温度に応じて変化させる必要はなく、固定値を用いても誤差が少ない場合には固定値を適用することも可能である。
そして、熱応力解析ステップにおいては、上記のごとく、途中応力解析ステップにおいて得られた上記鋳物要素の歪み、変位および応力を初期値として用い、上記鋳物要素の歪み、変位及び応力を求める。
この熱応力解析ステップを、上記冷却解析ステップにおける温度が室温になるまで繰り返すと、最終的に室温になった鋳物における残留歪みおよび残留応力を得ることができる。
すなわち、上記加圧要素としては、通常は、各要素にそれぞれ圧力が付与される状態で解析する場合が多い。しかし、この場合には、固相率が1の部分には実質的に変位が生じず、固相率が低い部分が集中して変位し、実際の加圧現象を十分に再現できない場合がある。
これに対し、上記のごとく加圧要素を剛体要素として、その全体が変位するようにして解析することにより、解析精度を高めることができる。
この場合の充填解析ステップとしては、例えば、いわゆるNavier-Stokesの式と連続の式を用いて湯流れ解析を行うと共に、いわゆる熱力学的エネルギの式を用いて充填時の温度解析を行う。さらに、充填後の鋳物要素の冷却を解析するステップを加え、上記凝固解析ステップに連続的に繋ぐこともできる。
この場合には、実際の鋳造により近いシミュレーションを行うことができ、その解析精度を更に高めることができる。なお、ステップが増加する分、解析時間の増加は否めない。
本発明の実施例に係るダイカストシミュレーション方法及びその装置につき、図1〜図8を用いて説明する。
本例のダイカストシミュレーション装置1は、図2に示すごとく、各解析ステップを実行する手段を備えたコンピュータ10を有する装置である。コンピュータ10には、キーボード、記憶媒体読み取り装置等の入力装置11と、ディスプレイ、プリンタ等の出力装置12が接続されている。
以下、これをさらに詳説する。
要素作成ステップS101は、上記CADデータを基に、鋳型部分7に位置する鋳型要素、鋳型内の空間(鋳物8、注入加圧部81及び第1ランナ部82、第2ランナ部83を含む部分)に位置する鋳物要素、及び該鋳物要素に加圧力を加える部分(プランジャ70に当接する部分85)に位置する加圧要素として、それぞれ複数の微小要素に分割するステップである。
凝固解析ステップS301では、上記鋳物要素のすべてに溶湯が充填された状態の熱的特性値を付与する。本例で鋳物要素に与える熱的特性値は、比熱、熱伝導率、密度、及び鋳型との熱伝達率である。また、鋳型要素にも比熱、熱伝導率、密度を与える。なお、これらの熱的特性値は、予め上記要素作成ステップにおいて与えておき、この凝固解析ステップS301においてあらためて置き換えることもできる。
そして、凝固解析ステップS301では、上記熱的特性値に基づいて微小時間後の上記鋳物要素の温度と固相率を求め出力する。
具体的な計算方法は、様々な方法があるが、本例では下記の式[数1]に代表される熱力学的エネルギの式及び下記の式[数2]に代表される固相率fsを求める式を用いて行う。
途中応力解析ステップS302では、まず、凝固解析ステップS301により得られた温度及び固相率を基にして鋳物要素の機械的特性値を求める。本例では、鋳物すなわち溶湯がアルミニウム合金であるので、固相率が0.3〜1の場合には、ヤング率(GPa)=(10+固相率×40)、ポアソン比=(0.33−0.01×固相率)を用い、固相率が0.3未満の場合には、ヤング率(GPa)=(7+固相率×50)、ポアソン比=(0.44−0.99×固相率)を用いてヤング率及びポアソン比を用いた。また、本例では、線膨張係数として固定値を用いた。
そして、この途中応力解析ステップS302では、その時の上記機械的特性値と加圧力とを基にして、鋳物要素の歪み、変位及び応力を経時的に求め、出力する。
そして、熱応力解析ステップS402では、上記の熱的特性値を基にして、かつ、上記途中応力解析ステップにおいて得られた上記鋳物要素の歪み、変位および応力を初期値として用い、上記鋳物要素の歪み、変位及び応力を求める。
図7には、上記鋳物8のA部(図3、図4)における歪みの履歴を出力した例を示す。同図は、横軸に溶湯の射出を開始してから室温までの冷却が完了するまでの鋳造時間をとり、縦軸に歪みをとったものである。
そして、プロットした計算結果E1が今回の解析結果であるA部の歪みの経時的変化である。
同図から知られるように、本例の場合には、計算結果E1と実測値C1との差が比較的大きいが、歪みの変化傾向は忠実に再現できていることがわかる。
この計算結果C2の結果(図8)と上記計算結果E1(図7)との比較から知られるように、凝固途中の上記途中応力解析ステップS302において鋳造圧力を十分に考慮することにより、解析結果が向上することが分かる。
本例では、実施例1における途中応力解析ステップS302における加圧力を加える際に、その加圧要素を剛体要素に設定して解析したものである。剛体要素であるか否かによる変位の違いを図9及び図10を用いて説明する。
図9は、加圧要素が剛体要素でない場合の説明図であり、同図(a)は、プランジャに当接している加圧要素85のそれぞれに加圧力ΔPをかけた状態を示してある。また、同図(b)には、固相率の状態を示してある。固相率は、鋳型に当接している部分が1、中央が0で分布している。この状態においては、それぞれ加圧力ΔPを加えると、変位は、変位線S1で示すように中央部が最も変位する。
同図は、実施例1の場合と同様に、鋳物8のA部(図3、図4)における歪みの履歴を出力したものであり、本例の結果は計算結果E2としてプロットしてある。また、同図にも、実測値を実測値C1として示してある。
本例は、実施例1における途中応力解析ステップS302で用いたヤング率とポアソン比とを一定値に変更して解析した場合における、解析結果の誤差の拡大度合いを調べた例である。
実施例1における計算結果E1と実測値C1との誤差(%)と、本例における計算結果と実測値C1との誤差(%)を棒グラフで示したものが図12である。なお、これらの誤差(%)は、解析終了時における計算結果をa、実測値をbとした場合に、|a−b|/b×100(%)の式により得られたものである。
10 コンピュータ
21 要素作成ステップ実行手段
23 凝固解析ステップ実行手段
24 途中応力解析ステップ実行手段
25 冷却解析ステップと実行手段
26 熱応力解析ステップ実行手段26
S101 要素作成ステップ
S301 凝固解析ステップ
S302 途中応力解析ステップ
S401 冷却解析ステップ
S402 熱応力解析ステップ
Claims (8)
- 鋳型内に溶湯を加圧注入して凝固させることにより所望形状の鋳物を得るダイカストを行う際に、少なくとも上記鋳物に生じる歪み及び応力を求めるダイカストシミュレーション方法であって、
上記鋳型の形状データを基に、鋳型部分に位置する鋳型要素、鋳型内の空間に位置する鋳物要素、及び該鋳物要素に加圧力を加える部分に位置する加圧要素として、それぞれ複数の微小要素に分割する要素作成ステップと、
上記鋳物要素のすべてに溶湯が充填された状態の熱的特性値を付与し、該熱的特性値に基づいて微小時間後の上記鋳物要素の温度と固相率を経時的に求める凝固解析ステップと、
該凝固解析ステップにより得られた温度及び/又は固相率を基にして上記鋳物要素の機械的特性値を求めると共に、上記加圧要素によって上記鋳物要素に加圧力を付加し、該加圧力と上記機械的特性値とを基にして、上記鋳物要素の歪み、変位及び応力を経時的に求める途中応力解析ステップと、
上記鋳物要素のすべてが固相になった後に、熱的特性値に基づいて微小時間後の上記鋳物要素の温度を経時的に求める冷却解析ステップと、
該冷却解析ステップにより得られた温度を基にして上記鋳物要素の機械的特性値を求めると共に、該機械的特性値を基にして、かつ、上記途中応力解析ステップにおいて得られた上記鋳物要素の歪み、変位および応力を初期値として用い、上記鋳物要素の歪み、変位及び応力を経時的に求める熱応力解析ステップとを有することを特徴とするダイカストシミュレーション方法。 - 請求項1において、上記途中応力解析ステップにおける上記加圧要素は、剛体要素とすることを特徴とするダイカストシミュレーション方法。
- 請求項1又は2において、上記凝固解析ステップを開始する前には、溶湯が上記鋳物要素に充填される充填解析を経時的に行う充填解析ステップを行い、該充填解析ステップの解析結果に基づいて上記凝固解析ステップにおける上記熱的特性値の初期値を付与することを特徴とするダイカストシミュレーション方法。
- 請求項1〜3のいずれか1項において、上記途中応力解析ステップにおける上記機械的特性値としては、少なくとも、固相率の値を変数とするヤング率算出式により求めたヤング率と、固相率の値を変数とするポアソン比算出式により求めたポアソン比とを用いることを特徴とするダイカストシミュレーション方法。
- 請求項4において、上記ヤング率算出式および/または上記ポアソン比算出式としては、固相率の範囲を複数の区分に分けた場合の当該区分毎に異なる式を用いることを特徴とするダイカストシミュレーション方法。
- 請求項1〜5のいずれか1項に記載のダイカストシミュレーション方法をコンピュータに実行させるためのダイカストシミュレーションプログラム。
- 請求項6に記載のダイカストシミュレーションプログラムを記録したコンピュータ読み取り可能な記憶媒体。
- 請求項1〜5のいずれか1項に記載のダイカストシミュレーション方法を実行するよう構成されたコンピュータを有することを特徴とするダイカストシミュレーション装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004213356A JP4232707B2 (ja) | 2004-07-21 | 2004-07-21 | ダイカストシミュレーション方法、その装置及びそのプログラム、並びに当該プログラムを記録した記録媒体 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004213356A JP4232707B2 (ja) | 2004-07-21 | 2004-07-21 | ダイカストシミュレーション方法、その装置及びそのプログラム、並びに当該プログラムを記録した記録媒体 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2006026723A JP2006026723A (ja) | 2006-02-02 |
JP4232707B2 true JP4232707B2 (ja) | 2009-03-04 |
Family
ID=35893673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004213356A Expired - Fee Related JP4232707B2 (ja) | 2004-07-21 | 2004-07-21 | ダイカストシミュレーション方法、その装置及びそのプログラム、並びに当該プログラムを記録した記録媒体 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4232707B2 (ja) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4289367B2 (ja) | 2006-05-22 | 2009-07-01 | トヨタ自動車株式会社 | 鋳造部品特性推定装置 |
JP4692402B2 (ja) * | 2006-06-12 | 2011-06-01 | 株式会社豊田中央研究所 | 鋳造シミュレーション方法、その装置、そのプログラム、及び当該プログラムを記録した記録媒体、並びに鋳造方法 |
JP4582067B2 (ja) * | 2006-08-02 | 2010-11-17 | トヨタ自動車株式会社 | 形状予測装置と形状予測方法とそのためのコンピュータプログラム |
JP4914429B2 (ja) * | 2008-12-05 | 2012-04-11 | 株式会社豊田中央研究所 | 合金溶湯の凝固解析方法およびその凝固解析プログラム |
JP6268584B2 (ja) * | 2014-01-14 | 2018-01-31 | 日産自動車株式会社 | 熱変形解析方法、熱変形解析プログラム、および熱変形解析装置 |
JP7279607B2 (ja) * | 2019-10-02 | 2023-05-23 | マツダ株式会社 | 鋳造品の弾塑性応力解析方法、解析システム、解析プログラム、及び記録媒体 |
-
2004
- 2004-07-21 JP JP2004213356A patent/JP4232707B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2006026723A (ja) | 2006-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102395972B (zh) | 模具填充之后起模的仿真 | |
JP4692402B2 (ja) | 鋳造シミュレーション方法、その装置、そのプログラム、及び当該プログラムを記録した記録媒体、並びに鋳造方法 | |
JP4232707B2 (ja) | ダイカストシミュレーション方法、その装置及びそのプログラム、並びに当該プログラムを記録した記録媒体 | |
KR100611689B1 (ko) | 수지 성형품의 설계 지원 장치 및 방법 | |
US20230105651A1 (en) | Injection molding warp prediction | |
JP2010052019A (ja) | 砂型鋳物のシミュレーション方法 | |
JP4855866B2 (ja) | 射出成形解析方法、そり変形解析方法およびその装置 | |
JP4501844B2 (ja) | 鋳造解析装置及び鋳造解析方法 | |
Shinde et al. | Optimization of mold yield in multicavity sand castings | |
Anglada et al. | Prediction and validation of shape distortions in the simulation of high pressure die casting | |
JP7125906B2 (ja) | 抱付力評価方法及び収縮量評価方法 | |
WO2010064131A2 (en) | Molten alloy solidification analyzing method and solidification analyzing program for performing the same | |
JPH0622840B2 (ja) | 成形プロセスシミユレ−シヨンシステム | |
JP4032848B2 (ja) | 成形シミュレーション方法、成形シミュレーション装置及び成形シミュレーションプログラム並びに当該成形シミュレーションプログラムを記録したコンピュータ読みとり可能な記録媒体 | |
JP4018712B2 (ja) | 射出成形品の変形量算出方法、プログラム、記憶媒体、情報処理装置 | |
JP4935594B2 (ja) | 凝固解析方法および凝固解析装置 | |
JP2001293748A (ja) | 射出成形プロセスシミュレーション装置および形状精度予測方法 | |
JP4052006B2 (ja) | 成型シミュレーション方法、成型シミュレーション装置及び成型シミュレーションプログラム並びに当該成型シミュレーションプログラムを記録したコンピュータ読みとり可能な記録媒体 | |
JP5442242B2 (ja) | ダイカスト鋳造鋳物の鋳巣解析方法およびその鋳巣解析プログラム | |
JP5349859B2 (ja) | 成形品形状の予測方法、成形品の製造方法、成形品形状の予測プログラムとその記憶媒体 | |
JP2007125593A (ja) | 鋳造解析装置及び鋳造解析方法 | |
JP2009298035A (ja) | 金型設計方法 | |
JP4373353B2 (ja) | ダイカストシミュレーション方法、その装置及びそのプログラム、並びに当該プログラムを記録した記録媒体 | |
JP5889077B2 (ja) | 成形品収縮変形予測装置、成形品収縮変形予測方法及び成形品収縮変形予測プログラム | |
JP2003112349A (ja) | 射出成形解析方法、射出成形解析装置および射出成形解析プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070125 |
|
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: 20081118 |
|
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: 20081201 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111219 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111219 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121219 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121219 Year of fee payment: 4 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313532 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121219 Year of fee payment: 4 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121219 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131219 Year of fee payment: 5 |
|
LAPS | Cancellation because of no payment of annual fees |