(57)【要約】【課題】残留歪みおよび残留応力を精度よく解析することができるダイカストシミュレーション方法を提供すること。【解決手段】鋳型要素、鋳物要素、及び加圧要素としての微小要素に分割する要素作成ステップＳ１０１と、熱的特性値に基づいて鋳物要素の温度と固相率を求める凝固解析ステップＳ３０１と、鋳物要素の機械的特性値を求めると共に、加圧要素によって鋳物要素に加圧力を付加し、加圧力と機械的特性値とを基にして、鋳物要素の歪み、変位及び応力を求める途中応力解析ステップＳ３０２と、鋳物要素のすべてが固相になった後に、熱的特性値に基づいて鋳物要素の温度を求める冷却解析ステップＳ４０１と、途中応力解析ステップＳ３０２において得られた鋳物要素の歪み、変位および応力を初期値として用い、鋳物要素の歪み、変位及び応力を求める熱応力解析ステップＳ４０２とを有する。【選択図】図１
2006026723.tif 000002
【特許請求の範囲】
【請求項１】
鋳型内に溶湯を加圧注入して凝固させることにより所望形状の鋳物を得るダイカストを行う際に、少なくとも上記鋳物に生じる歪み及び応力を求めるダイカストシミュレーション方法であって、
上記鋳型の形状データを基に、鋳型部分に位置する鋳型要素、鋳型内の空間に位置する鋳物要素、及び該鋳物要素に加圧力を加える部分に位置する加圧要素として、それぞれ複数の微小要素に分割する要素作成ステップと、
上記鋳物要素のすべてに溶湯が充填された状態の熱的特性値を付与し、該熱的特性値に基づいて微小時間後の上記鋳物要素の温度と固相率を経時的に求める凝固解析ステップと、
該凝固解析ステップにより得られた温度及び／又は固相率を基にして上記鋳物要素の機械的特性値を求めると共に、上記加圧要素によって上記鋳物要素に加圧力を付加し、該加圧力と上記機械的特性値とを基にして、上記鋳物要素の歪み、変位及び応力を経時的に求める途中応力解析ステップと、
上記鋳物要素のすべてが固相になった後に、熱的特性値に基づいて微小時間後の上記鋳物要素の温度を経時的に求める冷却解析ステップと、
該冷却解析ステップにより得られた温度を基にして上記鋳物要素の機械的特性値を求めると共に、該機械的特性値を基にして、かつ、上記途中応力解析ステップにおいて得られた上記鋳物要素の歪み、変位および応力を初期値として用い、上記鋳物要素の歪み、変位及び応力を経時的に求める熱応力解析ステップとを有することを特徴とするダイカストシミュレーション方法。
【請求項２】
請求項１において、上記途中応力解析ステップにおける上記加圧要素は、剛体要素とすることを特徴とするダイカストシミュレーション方法。
【請求項３】
請求項１又は２において、上記凝固解析ステップを開始する前には、溶湯が上記鋳物要素に充填される充填解析を経時的に行う充填解析ステップを行い、該充填解析ステップの解析結果に基づいて上記凝固解析ステップにおける上記熱的特性値の初期値を付与することを特徴とするダイカストシミュレーション方法。
【請求項４】
請求項１〜３のいずれか１項において、上記途中応力解析ステップにおける上記機械的特性値としては、少なくとも、固相率の値を変数とするヤング率算出式により求めたヤング率と、固相率の値を変数とするポアソン比算出式により求めたポアソン比とを用いることを特徴とするダイカストシミュレーション方法。
【請求項５】
請求項４において、上記ヤング率算出式および／または上記ポアソン比算出式としては、固相率の範囲を複数の区分に分けた場合の当該区分毎に異なる式を用いることを特徴とするダイカストシミュレーション方法。
【請求項６】
請求項１〜５のいずれか１項に記載のダイカストシミュレーション方法をコンピュータに実行させるためのダイカストシミュレーションプログラム。
【請求項７】
請求項６に記載のダイカストシミュレーションプログラムを記録したコンピュータ読み取り可能な記憶媒体。
【請求項８】
請求項１〜５のいずれか１項に記載のダイカストシミュレーション方法を実行するよう構成されたコンピュータを有することを特徴とするダイカストシミュレーション装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ダイカスト鋳造（以下、単にダイカストという）により得られた鋳物に存在する残留歪みおよび残留応力を精度よく予測できるダイカストシミュレーションに関する。
【背景技術】
【０００２】
例えば、アルミニウム合金よりなる鋳物を製造する場合には、鋳型内に溶融した材料（溶湯）を加圧充填する方法であるダイカストが多用されている。
ダイカストにより得られる鋳物における課題としては、ひけ巣欠陥、空気等の巻き込み欠陥、湯回り不良欠陥等の他に、鋳物に残る残留歪みおよび残留応力がある。
【０００３】
従来、ダイカストシミュレーションとして、ひけ巣などの鋳造欠陥を予測するものが種々考案されてきたが、残留歪みおよび残留応力を解析するものはほとんど存在していなかった。また、残留歪みおよび残留応力を解析するものであったとしても、それは、凝固後の冷却過程における線膨張係数の変化に基づいて解析するものであり、実際に鋳物に現れる残留歪みおよび残留応力と大きな差異があり、信頼性が低いものであった。
この要因は、鋳造圧力は、ひけ巣部への溶湯の補給への駆動力に使われるのみであり、固相の変化には寄与しないと考えられていたことによる。
【０００４】
例えば、下記の特許文献１に開示された金属溶湯の流動凝固解析方法では、凝固解析において、型を微小要素に分割し、その微小要素毎に熱伝導および溶質移動を解析している。この方法では、鋳造圧力を固相内の溶湯流れの駆動力として取扱い、固相の変形はなく、固相間の溶湯流れのみに鋳造圧力（溶湯圧力）が寄与しているとしてある。
しかしながら、このように、鋳造圧力が固相間の溶湯流れのみに寄与し、残留歪みは冷却時（温度３００℃以下）の線膨張差により生じるという手法では、精度よく鋳物の残留歪みおよび残留応力を求めることができなかった。
【０００５】
【特許文献１】特開平１０−３４３２０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
本発明は、かかる従来の問題点に鑑みてなされたもので、残留歪みおよび残留応力を精度よく解析することができるダイカストシミュレーション方法、その装置及びそのプログラム、並びに当該プログラムを記録した記録媒体を提供しようとするものである。
【０００７】
第１の発明は、鋳型内に溶湯を加圧注入して凝固させることにより所望形状の鋳物を得るダイカストを行う際に、少なくとも上記鋳物に生じる歪み及び応力を求めるダイカストシミュレーション方法であって、
上記鋳型の形状データを基に、鋳型部分に位置する鋳型要素、鋳型内の空間に位置する鋳物要素、及び該鋳物要素に加圧力を加える部分に位置する加圧要素として、それぞれ複数の微小要素に分割する要素作成ステップと、
上記鋳物要素のすべてに溶湯が充填された状態の熱的特性値を付与し、該熱的特性値に基づいて微小時間後の上記鋳物要素の温度と固相率を経時的に求める凝固解析ステップと、
該凝固解析ステップにより得られた温度及び／又は固相率を基にして上記鋳物要素の機械的特性値を求めると共に、上記加圧要素によって上記鋳物要素に加圧力を付加し、該加圧力と上記機械的特性値とを基にして、上記鋳物要素の歪み、変位及び応力を経時的に求める途中応力解析ステップと、
上記鋳物要素のすべてが固相になった後に、熱的特性値に基づいて微小時間後の上記鋳物要素の温度を経時的に求める冷却解析ステップと、
該冷却解析ステップにより得られた温度を基にして上記鋳物要素の機械的特性値を求めると共に、該機械的特性値を基にして、かつ、上記途中応力解析ステップにおいて得られた上記鋳物要素の歪み、変位および応力を初期値として用い、上記鋳物要素の歪み、変位及び応力を経時的に求める熱応力解析ステップとを有することを特徴とするダイカストシミュレーション方法にある（請求項１）。
【０００８】
本発明のダイカストシミュレーション方法では、少なくとも、上記要素作成ステップ、凝固解析ステップ、途中応力解析ステップ、冷却解析ステップ、および熱応力解析ステップを行い、最終的に上記熱応力解析ステップを例えば室温まで行うことによって、鋳物に生じる残留歪みおよび残留応力を精度よく求めることができる。
【０００９】
すなわち、本発明の方法では、特に上記途中応力解析ステップを有している。ここでは、上記のごとく、凝固解析ステップにより得られた温度及び／又は固相率を基にして上記鋳物要素の機械的特性値を求めると共に、上記加圧要素によって上記鋳物要素に加圧力を付加し、該加圧力と上記機械的特性値とを基にして、上記鋳物要素の歪み、変位及び応力を経時的に求める。これにより、従来のシミュレーションでは求めなかった凝固途中における各鋳物要素の歪み、変位および応力が得られる。そのため、その後の上記熱応力解析ステップにおいては、この凝固過程に生じた各鋳物要素の歪み、変位および応力を前提として、これに積み重ねるように凝固後における各鋳物要素の歪み、変位および応力を求めることができる。
それ故、従来のような凝固途中における固相の変形がないという前提の基に行うシミュレーションに比べ、鋳物に残る残留歪みおよび残留応力の解析の精度を格段に高めることができる。そして、このダイカストシミュレーション方法を用いれば、残留歪みおよび残留応力の小さい鋳物が得られる鋳物設計を支援することができる。
【００１０】
第２の発明は、上記第１の発明のダイカストシミュレーション方法をコンピュータに実行させるためのダイカストシミュレーションプログラムにある（請求項６）。
第３の発明は、上記第２の発明のダイカストシミュレーションプログラムを記録したコンピュータ読み取り可能な記憶媒体にある（請求項７）。
第４の発明は、上記第１の発明のダイカストシミュレーション方法を実行するよう構成されたコンピュータを有することを特徴とするダイカストシミュレーション装置にある（請求項８）。
【００１１】
これら第２〜第４の発明は、いずれも、これらを用いて上記ダイカストシミュレーションを実施することができ、上述したように、鋳物に残る残留歪みおよび残留応力を従来よりも精度よく解析することができる。
【発明を実施するための最良の形態】
【００１２】
第１の発明においては、上記のごとく、少なくとも、要素作成ステップと、凝固解析ステップと、途中応力解析ステップと、冷却解析ステップと、熱応力解析ステップとを実施する。以下に、これらについてさらに説明する。
【００１３】
＜要素作成ステップ＞
要素作成ステップでは、取り込む鋳型の形状データとして様々な形式のデータを用いることができるが、例えば、国際規格ＩＧＥＳ（Initial Graphics Exchange Specification）の機能を介して取り込むことができる。
なお、鋳型内の空間部分（キャビティ部分）はできる限り正確な形状が必要であるが、鋳型部分の形状としては、最小限に絞ることができる。また、上記加圧要素としては、例えば、実際の鋳造時おいて溶湯に圧力をかけるプランジャーに接する鋳物部分を想定して定義することができる。なお、この要素作成ステップにおいて、上記各微小要素を定義すると共に、それぞれの物性値等を与えることもできるが、これは後の解析ステップにおいて行ってもよい。
また、上記微小要素の大きさはできる限り小さい方が精度向上のために好ましいが、小さいほど解析時間が長くなるので、適度な大きさにすることが好ましい。また、微小要素の大きさをすべての部分で統一する必要はなく、部位によって変更することができる。
【００１４】
＜凝固解析ステップ＞
凝固解析ステップは、上記鋳物要素のすべてに溶湯が充填された状態から開始する。そして、このステップにより解析を開始する際に、各鋳物要素の熱的特性値の初期値を付与する。ここでいう熱的特性値としては、例えば、密度、比熱、熱伝導率、潜熱等がある。また、これらの特性値は、温度等によって変化させることが好ましいが、大幅な精度低下につながらない限り一定値に固定することもできる。
なお、この凝固解析ステップを実施する前に、後述する充填解析ステップを実施して、溶湯の挙動を考慮して溶湯充填完了時あるいはその後の冷却を経て凝固開始する際における各鋳物要素の熱的特性値の初期値を定めることもできる。
【００１５】
そして、この凝固解析ステップでは、上記の熱的特性値に基づいて、微小時間後の上記鋳物要素の温度と固相率を求める。ここで、固相率（ｆｓ）とは、その微小要素の体積をＡ、固相の部分の体積をＢとした場合の、Ｂ／Ａで示される値であって、その要素全体が液相の場合を０、要素全体が固相の場合を１として表す値である。
この凝固解析ステップは、鋳物要素のすべてが固相のみとなるまで、すなわち、すべての固相率が１になるまで繰り返し行われる。
【００１６】
また、凝固解析ステップにおいて鋳物要素の温度と固相率の求め方は様々な手法により行うことができる。
例えば、鋳物要素と鋳型要素のそれぞれの比熱、熱伝導率、および密度と、鋳物要素と鋳型要素の間の熱伝達率等を用いた熱力学的エネルギの式を基にして熱伝導解析を行って温度を求め、その後、溶湯の潜熱放出を考慮して固相率を求める手法がある。
また、上記の鋳物要素の比熱をみかけ比熱（潜熱放出を考慮した比熱）に代えて熱伝導解析を行って温度を求め、その後、シャイルの式を用いて固相率を求める手法もある。
【００１７】
＜途中応力解析ステップ＞
途中応力解析ステップでは、まず、上記凝固解析ステップにより温度及び固相率が得られる毎に、これらを基にして上記鋳物要素の機械的特性値を求める。ここでいう機械的特性値は、例えば、ヤング率、ポアソン比、線膨張係数、耐力、加工硬化指数等がある。なお、これらの機械的特性値は、上記温度あるいは固相率に基づいて必ずしも変化させる必要はなく、固定値を用いても誤差が少ない場合には固定値を適用することも可能である。
【００１８】
また、途中応力解析ステップでは、上記加圧要素によって上記鋳物要素に加圧力を付加する。この場合の加圧力の付加方法としては、例えば、実際の溶湯の加圧圧力に相当する加圧力を一気に加え、以後同一圧力をキープする方法、実際の加圧圧力に相当するまで徐々に加圧力を高めて行く方法等がある。
そして、上記加圧力と上記機械的特性値を基にして、上記鋳物要素の歪み、変位及び応力を求める。これを上記凝固解析ステップの結果が出る度に繰り返し行い、鋳物要素のすべてが固相、すなわちすべての固相率が１となるまで行う。
【００１９】
＜冷却解析ステップ＞
冷却解析ステップは、上記鋳物要素のすべてが固相（固相率が１）になった後、つまり、凝固が完了した後に実施し、上記熱的特性値に基づいて微小時間後の上記鋳物要素の温度を経時的に求める。この段階では、固相を求める必要がないので、少なくとも鋳物要素の温度の変化を求める。
【００２０】
＜熱応力解析ステップ＞
熱応力解析ステップでは、上記途中応力解析ステップによって得られた各鋳物要素の歪み、変位および応力を初期値として、凝固完了後の歪み、変位および応力を引き続き求めていく。
すなわち、熱応力解析ステップでは、まず上記冷却解析ステップにより得られた温度を基にして上記鋳物要素の機械的特性値を求める。ここで用いる機械的特性値の種類としては、上記途中応力解析ステップにおいて用いたものと同様である。また、ここで、上記機械的特性値は、必ずしも温度に応じて変化させる必要はなく、固定値を用いても誤差が少ない場合には固定値を適用することも可能である。
そして、熱応力解析ステップにおいては、上記のごとく、途中応力解析ステップにおいて得られた上記鋳物要素の歪み、変位および応力を初期値として用い、上記鋳物要素の歪み、変位及び応力を求める。
この熱応力解析ステップを、上記冷却解析ステップにおける温度が室温になるまで繰り返すと、最終的に室温になった鋳物における残留歪みおよび残留応力を得ることができる。
【００２１】
次に、上記第１の発明においては、上記途中応力解析ステップにおける上記加圧要素は、剛体要素とすることが好ましい（請求項２）。
すなわち、上記加圧要素としては、通常は、各要素にそれぞれ圧力が付与される状態で解析する場合が多い。しかし、この場合には、固相率が１の部分には実質的に変位が生じず、固相率が低い部分が集中して変位し、実際の加圧現象を十分に再現できない場合がある。
これに対し、上記のごとく加圧要素を剛体要素として、その全体が変位するようにして解析することにより、解析精度を高めることができる。
【００２２】
また、上記凝固解析ステップを開始する前には、溶湯が上記鋳物要素に充填される充填解析を経時的に行う充填解析ステップを行い、該充填解析ステップの解析結果に基づいて上記凝固解析ステップにおける上記熱的特性値の初期値を付与することが好ましい（請求項３）。
この場合の充填解析ステップとしては、例えば、いわゆるNavier-Stokesの式と連続の式を用いて湯流れ解析を行うと共に、いわゆる熱力学的エネルギの式を用いて充填時の温度解析を行う。さらに、充填後の鋳物要素の冷却を解析するステップを加え、上記凝固解析ステップに連続的に繋ぐこともできる。
この場合には、実際の鋳造により近いシミュレーションを行うことができ、その解析精度を更に高めることができる。なお、ステップが増加する分、解析時間の増加は否めない。
【００２３】
また、上記途中応力解析ステップにおける上記機械的特性値としては、少なくとも、固相率の値を変数とするヤング率算出式により求めたヤング率と、固相率の値を変数とするポアソン比算出式により求めたポアソン比とを用いることが好ましい（請求項４）。すなわち、上記のごとく機械的特性値としては、固定値を用いることも可能であるが、ヤング率とポアソン比とは、固相率を変数とする上記ヤング率算出式とポアソン比算出式とを用いて求めることが好ましい。固相率に応じてヤング率とポアソン比を積極的に変更することにより、各鋳物要素をあたかも固相のような状態で扱いつつ、その鋳物要素における計算上の歪み、変位及び応力を妥当な値として求めることができる。
【００２４】
また、この場合の上記ヤング率算出式および／または上記ポアソン比算出式としては、固相率の範囲を複数の区分に分けた場合の当該区分毎に異なる式を用いることが好ましい（請求項５）。例えば一例を示すと、上記溶湯がアルミニウム合金の場合には、固相率が０．３〜１の場合には、ヤング率（ＧＰａ）＝（１０＋固相率×４０）、ポアソン比＝（０．３３−０．０１×固相率）を用い、固相率が０．３未満の場合には、ヤング率（ＧＰａ）＝（７＋固相率×５０）、ポアソン比＝（０．４４−０．９９×固相率）を用いることができ、この場合にはさらに解析精度を高めることができる。
【実施例】
【００２５】
（実施例１）
本発明の実施例に係るダイカストシミュレーション方法及びその装置につき、図１〜図８を用いて説明する。
本例のダイカストシミュレーション装置１は、図２に示すごとく、各解析ステップを実行する手段を備えたコンピュータ１０を有する装置である。コンピュータ１０には、キーボード、記憶媒体読み取り装置等の入力装置１１と、ディスプレイ、プリンタ等の出力装置１２が接続されている。
【００２６】
コンピュータ１０が実行可能な手段は、後述するダイカストシミュレーション方法における各ステップを実行する手段であり、具体的には、要素作成ステップ実行手段２１と、凝固解析ステップ実行手段２３と、途中応力解析ステップ実行手段２４と、冷却解析ステップと実行手段２５と、熱応力解析ステップ実行手段２６とを有している。なお、これらの実行手段は、後述するダイカストシミュレーション方法をコンピュータに実行させるためのダイカストシミュレーションプログラムとして作成したものを読み込んで実行可能としたものである。なお、同図に示すごとく、オプションとして充填解析手段実行手段を２２を追加することもできるが、本例では省略してある。
【００２７】
本例のダイカストシミュレーション方法は、図１に示すごとく、要素作成ステップＳ１０１と、凝固解析ステップＳ３０１と、途中応力解析ステップＳ３０２と、冷却解析ステップＳ４０１と、熱応力解析ステップＳ４０２とを行う方法である。なお、この方法においては、凝固解析ステップＳ３０１の前に溶湯を充填する際の解析を行う充填解析ステップＳ２０１をオプションとして加えることができるが、本例では加えていない。
以下、これをさらに詳説する。
【００２８】
まず、本例では、図３、図４に示すごとく、アルミニウム合金よりなる直方体状の鋳物８を作ることを想定し、また、溶湯の注入加圧部８１及び第１ランナ部８２、第２ランナ部８３も鋳物の一部としたモデルを想定する。そして、これらを囲う部分が鋳型部分７であり、上記注入加圧部８１を押圧するプランジャ７０に接する部分８５を加圧部分として想定する。そして、これらを空間とした鋳型形状をＣＡＤデータとして作成し、このＣＡＤデータを用いて要素作成ステップＳ１０１を行う。
要素作成ステップＳ１０１は、上記ＣＡＤデータを基に、鋳型部分７に位置する鋳型要素、鋳型内の空間（鋳物８、注入加圧部８１及び第１ランナ部８２、第２ランナ部８３を含む部分）に位置する鋳物要素、及び該鋳物要素に加圧力を加える部分（プランジャ７０に当接する部分８５）に位置する加圧要素として、それぞれ複数の微小要素に分割するステップである。
【００２９】
次に、本例では、充填解析を行うことなく、凝固解析ステップＳ３０１を実施する。
凝固解析ステップＳ３０１では、上記鋳物要素のすべてに溶湯が充填された状態の熱的特性値を付与する。本例で鋳物要素に与える熱的特性値は、比熱、熱伝導率、密度、及び鋳型との熱伝達率である。また、鋳型要素にも比熱、熱伝導率、密度を与える。なお、これらの熱的特性値は、予め上記要素作成ステップにおいて与えておき、この凝固解析ステップＳ３０１においてあらためて置き換えることもできる。
そして、凝固解析ステップＳ３０１では、上記熱的特性値に基づいて微小時間後の上記鋳物要素の温度と固相率を求め出力する。
具体的な計算方法は、様々な方法があるが、本例では下記の式［数１］に代表される熱力学的エネルギの式及び下記の式［数２］に代表される固相率ｆｓを求める式を用いて行う。
【００３０】
【数１】
2006026723.tif 000003
【００３１】
【数２】
2006026723.tif 000004
【００３２】
次に、上記凝固解析ステップＳ３０１の解析結果が出力される度に上記途中応力解析ステップＳ３０２を行う。
途中応力解析ステップＳ３０２では、まず、凝固解析ステップＳ３０１により得られた温度及び固相率を基にして鋳物要素の機械的特性値を求める。本例では、鋳物すなわち溶湯がアルミニウム合金であるので、固相率が０．３〜１の場合には、ヤング率（ＧＰａ）＝（１０＋固相率×４０）、ポアソン比＝（０．３３−０．０１×固相率）を用い、固相率が０．３未満の場合には、ヤング率（ＧＰａ）＝（７＋固相率×５０）、ポアソン比＝（０．４４−０．９９×固相率）を用いてヤング率及びポアソン比を用いた。また、本例では、線膨張係数として固定値を用いた。
【００３３】
さらに、このステップでは、上記加圧要素によって上記鋳物要素に加圧力を付加する。本例では、加圧力としては、最終的に実際に付加する鋳造圧力を想定して加えるが、その加圧力は徐々に増やしていき、凝固完了までの間に最大圧力に到達するように設定した。また、加圧力のかけ方は、それぞれの加圧要素に与える方法で行い、剛体要素としての設定は行わなかった。
そして、この途中応力解析ステップＳ３０２では、その時の上記機械的特性値と加圧力とを基にして、鋳物要素の歪み、変位及び応力を経時的に求め、出力する。
【００３４】
次に、凝固完了判断ステップＳ３０３では、上記凝固解析ステップＳ３０１における出力値の１つである固相率が、すべての鋳物要素において１になったか否かを判定する。すべての鋳物要素の固相率が１になっていない場合には、上記凝固解析ステップＳ３０１と途中応力解析ステップＳ３０２を繰り返し行う。一方、すべての鋳物要素において固相率が１になったと判定した場合には、次の冷却解析ステップＳ４０１に移る。
【００３５】
冷却解析ステップＳ４０１では、鋳物要素における熱的特性値に基づいて微小時間後の上記鋳物要素の温度を求める。このときの熱的特性値の種類としては、上記凝固解析ステップＳ３０１の場合と同様である。なお、その値自体は、必要に応じて変更することができる。そして、温度の算出は、上記と同様に熱力学エネルギの式を用いて算出する。
【００３６】
次に、熱応力解析ステップＳ４０２では、冷却解析ステップＳ４０１により得られた温度を基にして鋳物要素の機械的特性値を求める。このときの機械的特性値としては、上記途中応力解析ステップＳ３０２で用いた特性値と同様な種類のものを用いる。なお、ヤング率とポアソン比とは、上述した算出式において固相率＝１を代入して得られる固定値を用いる。
そして、熱応力解析ステップＳ４０２では、上記の熱的特性値を基にして、かつ、上記途中応力解析ステップにおいて得られた上記鋳物要素の歪み、変位および応力を初期値として用い、上記鋳物要素の歪み、変位及び応力を求める。
【００３７】
次に、冷却完了判断ステップＳ４０３では、上記冷却解析ステップＳ４０１における出力値の温度が、室温、具体的には２０℃に到達したか否かを判定する。すべての鋳物要素の温度が室温以下の温度になっていない場合には、上記冷却解析ステップＳ４０１と熱応力解析ステップＳ４０２を繰り返し行う。一方、すべての鋳物要素において温度が室温に到達したと判定した場合には、すべての解析を終了する。
【００３８】
なお、ダイカストシミュレーション装置１は、上記各ステップにおける出力値をメモリ部分（図示略）に記憶し、様々な形態で出力できるように構成してある。
図７には、上記鋳物８のＡ部（図３、図４）における歪みの履歴を出力した例を示す。同図は、横軸に溶湯の射出を開始してから室温までの冷却が完了するまでの鋳造時間をとり、縦軸に歪みをとったものである。
そして、プロットした計算結果Ｅ１が今回の解析結果であるＡ部の歪みの経時的変化である。
【００３９】
また、本例では、図５、図６に示すごとく、シミュレーションした鋳型と同じ形状の鋳型を実際に作製し、実際にダイカストを行うと共に上記Ａ部に相当する部分の歪みを実測した。具体的には、これらの図に示すごとく、円筒状のステンレス鋼よりなるダミー材８８の内周面に歪みゲージ８９を貼設しておき、このダミー材８８を鋳ぐるむようにダイカストを行ってその過程で歪みを実測した。この実測値をＣ１として図７に合わせて示す。
同図から知られるように、本例の場合には、計算結果Ｅ１と実測値Ｃ１との差が比較的大きいが、歪みの変化傾向は忠実に再現できていることがわかる。
【００４０】
また、本例では、比較のために、上記途中応力解析ステップＳ３０２において、加圧力をまったく与えずに、その他は上記と同様にして解析を行った。その結果を比較例としての計算結果Ｃ２として図８に示す。同図にも、上記と同様の実測値Ｃ１を並記した。
この計算結果Ｃ２の結果（図８）と上記計算結果Ｅ１（図７）との比較から知られるように、凝固途中の上記途中応力解析ステップＳ３０２において鋳造圧力を十分に考慮することにより、解析結果が向上することが分かる。
【００４１】
（実施例２）
本例では、実施例１における途中応力解析ステップＳ３０２における加圧力を加える際に、その加圧要素を剛体要素に設定して解析したものである。剛体要素であるか否かによる変位の違いを図９及び図１０を用いて説明する。
図９は、加圧要素が剛体要素でない場合の説明図であり、同図（ａ）は、プランジャに当接している加圧要素８５のそれぞれに加圧力ΔＰをかけた状態を示してある。また、同図（ｂ）には、固相率の状態を示してある。固相率は、鋳型に当接している部分が１、中央が０で分布している。この状態においては、それぞれ加圧力ΔＰを加えると、変位は、変位線Ｓ１で示すように中央部が最も変位する。
【００４２】
一方、図１０は、加圧要素が剛体要素である場合の説明図であり、同図（ａ）は、プランジャに当接している加圧要素８５全体に加圧力Ｐがかけられた状態を示してある。この場合の固相率の分布は、同図（ｂ）に示すごとく上記と同様である。この状態においては、加圧力Ｐを加えると、剛体要素であるので、全体が同じように変位して、変位線Ｓ２で示される変位の位置は直線状となる。
【００４３】
このような違いに基づき、本例では加圧要素を剛体要素に設定して、その他は実施例１と同様にして解析を実施した。その結果を図１１に示す。
同図は、実施例１の場合と同様に、鋳物８のＡ部（図３、図４）における歪みの履歴を出力したものであり、本例の結果は計算結果Ｅ２としてプロットしてある。また、同図にも、実測値を実測値Ｃ１として示してある。
【００４４】
同図と、先に示した図７との比較から知られるように、加圧要素を剛体要素とすることによって、計算結果と実測値との差が格段に減少し、解析精度が向上することが分かる。
【００４５】
（実施例３）
本例は、実施例１における途中応力解析ステップＳ３０２で用いたヤング率とポアソン比とを一定値に変更して解析した場合における、解析結果の誤差の拡大度合いを調べた例である。
実施例１における計算結果Ｅ１と実測値Ｃ１との誤差（％）と、本例における計算結果と実測値Ｃ１との誤差（％）を棒グラフで示したものが図１２である。なお、これらの誤差（％）は、解析終了時における計算結果をａ、実測値をｂとした場合に、｜ａ−ｂ｜／ｂ×１００（％）の式により得られたものである。
【００４６】
同図から知られるように、実施例１の場合には、本例（実施例３）よりも誤差が小さくより高い精度となることが分かる。そして、これにより、途中応力解析ステップＳ３０２で用いるヤング率及びポアソン比を上述した固相率の値を変数とする算出式により求めること、特に、固相率の範囲を複数の区分に分けて異なる式を用いることが非常に有効であることがわかる。
【図面の簡単な説明】
【００４７】
【図１】実施例１における、ダイカストシミュレーション方法を示すフロー図。
【図２】実施例１における、ダイカストシミュレーション装置の構成を示す説明図。
【図３】実施例１における、解析対象の鋳物形状を正面から見た説明図。
【図４】実施例１における、鋳型形状を断面から見た説明図。
【図５】実施例１における、ダミー材を鋳ぐるむ位置を鋳物の正面から見た説明図。
【図６】実施例１における、ダミー材を鋳ぐるむ位置を鋳物の断面から見た説明図。
【図７】実施例１における、計算結果Ｅ１の結果を示す説明図。
【図８】実施例１における、計算結果Ｃ２の結果を示す説明図。
【図９】実施例２における、剛体要素でない加圧要素の変位の状態を示す説明図。
【図１０】実施例２における、剛体要素である加圧要素の変位の状態を示す説明図。
【図１１】実施例２における、計算結果Ｅ２の結果を示す説明図。
【図１２】実施例３における、解析誤差の比較を示す説明図。
【符号の説明】
【００４８】
１  ダイカストシミュレーション装置
１０  コンピュータ
２１  要素作成ステップ実行手段
２３  凝固解析ステップ実行手段
２４  途中応力解析ステップ実行手段
２５  冷却解析ステップと実行手段
２６  熱応力解析ステップ実行手段２６
Ｓ１０１  要素作成ステップ
Ｓ３０１  凝固解析ステップ
Ｓ３０２  途中応力解析ステップ
Ｓ４０１  冷却解析ステップ
Ｓ４０２  熱応力解析ステップ
【図１】
2006026723.tif 000005
【図２】
2006026723.tif 000006
【図３】
2006026723.tif 000007
【図４】
2006026723.tif 000008
【図５】
2006026723.tif 000009
【図６】
2006026723.tif 000010
【図７】
2006026723.tif 000011
【図８】
2006026723.tif 000012
【図９】
2006026723.tif 000013
【図１０】
2006026723.tif 000014
【図１１】
2006026723.tif 000015
【図１２】
2006026723.tif 000016