JPH05138265A - 成形シミユレーシヨン方法及び装置 - Google Patents
成形シミユレーシヨン方法及び装置Info
- Publication number
- JPH05138265A JPH05138265A JP29925991A JP29925991A JPH05138265A JP H05138265 A JPH05138265 A JP H05138265A JP 29925991 A JP29925991 A JP 29925991A JP 29925991 A JP29925991 A JP 29925991A JP H05138265 A JPH05138265 A JP H05138265A
- Authority
- JP
- Japan
- Prior art keywords
- analysis
- particles
- region
- mold
- analysis region
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0216—Finite elements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/22—Moulding
Abstract
(57)【要約】
【目的】 有限要素法を基本として、金型と素材との接
触評価及びリメッシュを必要としない、パーティクル流
れモデルによる成形シミュレーション方法及び装置を提
供すること。 【構成】 固定された要素上をパーティクル（６）が移
動することにより、各要素自身が変形するものではな
く、要素の材料定数が変化するという記述形式を採る。
触評価及びリメッシュを必要としない、パーティクル流
れモデルによる成形シミュレーション方法及び装置を提
供すること。 【構成】 固定された要素上をパーティクル（６）が移
動することにより、各要素自身が変形するものではな
く、要素の材料定数が変化するという記述形式を採る。
Description
【０００１】
【産業上の利用分野】本発明は、塑性成形加工分野の成
形シミュレーション方法及びその装置に関し、特に製品
開発段階における事前予測としての素材の流動変形過程
を、数値シミュレーションする識別粒子（以下、パーテ
ィクルと呼ぶ）流れモデルによる成形シミュレーション
方法及びその装置に関する。
形シミュレーション方法及びその装置に関し、特に製品
開発段階における事前予測としての素材の流動変形過程
を、数値シミュレーションする識別粒子（以下、パーテ
ィクルと呼ぶ）流れモデルによる成形シミュレーション
方法及びその装置に関する。
【０００２】
【従来の技術】一般に、鍛造・圧延・押出し・射出など
の塑性成形加工分野では、適確な素材形状、金型および
工程設計がなされることにより、製品の品質向上及びコ
ストダウンが図れることは知られている。
の塑性成形加工分野では、適確な素材形状、金型および
工程設計がなされることにより、製品の品質向上及びコ
ストダウンが図れることは知られている。
【０００３】そこで、図２２、図２３及び図２４に示す
通り、従来は、通常の成形シミュレーション装置を用い
て、製品開発段階における事前予測としての素材の流動
変形過程を解析していた（平成元年度 塑性加工春季講
演会，１９８９．５．１８〜２０，「ギア素材の型鍛造
の弾塑性有限要素法解析」 後藤 学・渋谷 友次）。
簡単に言えば、図２２に示すような製品を成形するに際
して、図２３に示すような金型Ｄと素材Ｍとを想定す
る。そして、素材Ｍには図２３に示すような要素分割を
行い、分割された要素が成形過程でどのように変化する
かをシミュレートする。成形過程のシミュレーション結
果は図２４のように表示される。ここでは、図２４ｃに
おいて破線から実線で示すように要素再分割（以下、リ
メッシュと呼ぶ。）が行われている。
通り、従来は、通常の成形シミュレーション装置を用い
て、製品開発段階における事前予測としての素材の流動
変形過程を解析していた（平成元年度 塑性加工春季講
演会，１９８９．５．１８〜２０，「ギア素材の型鍛造
の弾塑性有限要素法解析」 後藤 学・渋谷 友次）。
簡単に言えば、図２２に示すような製品を成形するに際
して、図２３に示すような金型Ｄと素材Ｍとを想定す
る。そして、素材Ｍには図２３に示すような要素分割を
行い、分割された要素が成形過程でどのように変化する
かをシミュレートする。成形過程のシミュレーション結
果は図２４のように表示される。ここでは、図２４ｃに
おいて破線から実線で示すように要素再分割（以下、リ
メッシュと呼ぶ。）が行われている。
【０００４】具体的には、この従来法は、有限要素法プ
ログラムをベースとし、図２１に示すように、これに所
定の境界条件及びリメッシュに基づく処理ができるよう
に必要操作を加えて、基本的な成形過程の解析を行うも
のであった。
ログラムをベースとし、図２１に示すように、これに所
定の境界条件及びリメッシュに基づく処理ができるよう
に必要操作を加えて、基本的な成形過程の解析を行うも
のであった。
【０００５】換言すれば、有限要素法による成形シミュ
レーションプログラムをベースとした従来の解析手法
は、素材のみ、或いは素材と金型とを有限要素でモデル
化し、このモデル化した有限要素自体の変形を追跡する
ことで、素材の流動変形過程を数値シミュレーションし
ていた。
レーションプログラムをベースとした従来の解析手法
は、素材のみ、或いは素材と金型とを有限要素でモデル
化し、このモデル化した有限要素自体の変形を追跡する
ことで、素材の流動変形過程を数値シミュレーションし
ていた。
【０００６】
【発明が解決しようとする課題】しかしながら、従来の
解析手法においては、図１９に示す通り金型と素材との
接触評価が必要であり、演算エラーの原因となってい
た。
解析手法においては、図１９に示す通り金型と素材との
接触評価が必要であり、演算エラーの原因となってい
た。
【０００７】さらに、従来の解析手法においては、有限
要素自身を変形させることから、変形が進につれて、図
２０に示すように要素がいびつになり、それ以上の計算
が不可能となる。この大変形に対処するためのリメッシ
ュは、ユーザにとって極めて煩わしい作業となる欠点が
あった。
要素自身を変形させることから、変形が進につれて、図
２０に示すように要素がいびつになり、それ以上の計算
が不可能となる。この大変形に対処するためのリメッシ
ュは、ユーザにとって極めて煩わしい作業となる欠点が
あった。
【０００８】そこで、本発明の技術的課題は、上記欠点
に鑑み、有限要素法を基本として、金型と素材との接触
評価及びリメッシュを必要としない、パーティクル流れ
モデルによる成形シミュレーション方法及びその装置を
提供することである。
に鑑み、有限要素法を基本として、金型と素材との接触
評価及びリメッシュを必要としない、パーティクル流れ
モデルによる成形シミュレーション方法及びその装置を
提供することである。
【０００９】
【課題を解決するための手段】本発明によれば、金型内
に素材を含む解析空間を規定して、前記金型の移動によ
る前記素材の流動変形過程を、有限要素法により数値シ
ミュレートする成形シミュレーション方法において、前
記解析空間には、前記素材が占める第１の解析領域と、
前記素材が存在しない残りの第２の解析領域とを設定
し、前記第１の解析領域と前記第２の解析領域とを解析
対象として有限要素でモデル化して有限要素法により前
記素材の流動変形過程を数値シミュレートし、前記モデ
ルは、前記解析対象を固定もしくは準固定した複数の分
割要素とし、前記第１の解析領域には、前記素材形状を
規定する複数のパーティクルを分散して設定し、前記パ
ーティクルを前記解析空間内で移動させ、移動した領域
における材料定数を変化させるようなモデルとすること
を特徴とする成形シミュレーション方法が得られる。
に素材を含む解析空間を規定して、前記金型の移動によ
る前記素材の流動変形過程を、有限要素法により数値シ
ミュレートする成形シミュレーション方法において、前
記解析空間には、前記素材が占める第１の解析領域と、
前記素材が存在しない残りの第２の解析領域とを設定
し、前記第１の解析領域と前記第２の解析領域とを解析
対象として有限要素でモデル化して有限要素法により前
記素材の流動変形過程を数値シミュレートし、前記モデ
ルは、前記解析対象を固定もしくは準固定した複数の分
割要素とし、前記第１の解析領域には、前記素材形状を
規定する複数のパーティクルを分散して設定し、前記パ
ーティクルを前記解析空間内で移動させ、移動した領域
における材料定数を変化させるようなモデルとすること
を特徴とする成形シミュレーション方法が得られる。
【００１０】なお、前記数値シミュレートは、前記複数
のパーティクルを前記複数の分割要素上を移動させるこ
とによる前記複数の分割要素の物理量の変化に基づい
て、有限要素法により前記素材の流動変形過程をシミュ
レートするものである。
のパーティクルを前記複数の分割要素上を移動させるこ
とによる前記複数の分割要素の物理量の変化に基づい
て、有限要素法により前記素材の流動変形過程をシミュ
レートするものである。
【００１１】本発明によればまた、金型内に素材を含む
解析空間を規定して、前記金型の移動による前記素材の
流動変形過程を、有限要素法により数値シミュレートす
る成形シミュレーション装置において、前記解析空間
に、前記素材が占める第１の解析領域と、前記素材が存
在しない残りの第２の解析領域とを設定し、前記第１の
解析領域と前記第２の解析領域とを解析対象として有限
要素でモデル化する演算モデル構成部と、前記金型と前
記素材とに基づいて、予め規定される境界条件の下で、
前記解析空間全体におけるエネルギー消散率が最小とな
るように、前記複数のパーティクルが存在する前記第１
の解析領域における前記分割要素の節点速度と、前記複
数のパーティクルが存在しない前記第２の解析領域にお
ける前記分割要素の節点速度とを、前記解析空間内の速
度場として算出する速度場算出部と、前記速度場に基づ
いて、前記複数のパーティクルの移動速度を算出するパ
ーティクル移動速度算出部と、該パーティクル移動速度
算出部より算出された移動速度に基づいて、移動後にお
ける前記複数のパーティクルの位置情報を算出する移動
位置算出部と、前記速度場に基づいて、前記分割要素の
物理量を算出する物理量算出部とを有することを特徴と
する成形シミュレーション装置が得られる。
解析空間を規定して、前記金型の移動による前記素材の
流動変形過程を、有限要素法により数値シミュレートす
る成形シミュレーション装置において、前記解析空間
に、前記素材が占める第１の解析領域と、前記素材が存
在しない残りの第２の解析領域とを設定し、前記第１の
解析領域と前記第２の解析領域とを解析対象として有限
要素でモデル化する演算モデル構成部と、前記金型と前
記素材とに基づいて、予め規定される境界条件の下で、
前記解析空間全体におけるエネルギー消散率が最小とな
るように、前記複数のパーティクルが存在する前記第１
の解析領域における前記分割要素の節点速度と、前記複
数のパーティクルが存在しない前記第２の解析領域にお
ける前記分割要素の節点速度とを、前記解析空間内の速
度場として算出する速度場算出部と、前記速度場に基づ
いて、前記複数のパーティクルの移動速度を算出するパ
ーティクル移動速度算出部と、該パーティクル移動速度
算出部より算出された移動速度に基づいて、移動後にお
ける前記複数のパーティクルの位置情報を算出する移動
位置算出部と、前記速度場に基づいて、前記分割要素の
物理量を算出する物理量算出部とを有することを特徴と
する成形シミュレーション装置が得られる。
【００１２】
【作用】すなわち、本発明は、図３に示すように、成形
素材のみでなく、素材が存在しない領域も解析対象と
し、固定した分割要素上をパーティクルが移動すること
により、素材の変形挙動を表現できるようにしている。
素材のみでなく、素材が存在しない領域も解析対象と
し、固定した分割要素上をパーティクルが移動すること
により、素材の変形挙動を表現できるようにしている。
【００１３】従来のシミュレーション装置における演算
部は、素材の変形量を逐次要素自身の変形として算出す
るのに対し、本発明のシミュレーション装置における演
算部は、素材をパーティクルで表し、第１の解析領域内
に分布したパーティクルを全解析領域で移動させること
により素材の変形挙動を表現するものである。
部は、素材の変形量を逐次要素自身の変形として算出す
るのに対し、本発明のシミュレーション装置における演
算部は、素材をパーティクルで表し、第１の解析領域内
に分布したパーティクルを全解析領域で移動させること
により素材の変形挙動を表現するものである。
【００１４】換言すれば、本発明は、パーティクルを解
析領域内で移動させ、素材が存在していた第１の解析領
域から素材が存在していなかった第２の解析領域へ移動
した場合、移動した領域における材料定数が変化すると
いう記述形式を採るものであると言える。
析領域内で移動させ、素材が存在していた第１の解析領
域から素材が存在していなかった第２の解析領域へ移動
した場合、移動した領域における材料定数が変化すると
いう記述形式を採るものであると言える。
【００１５】
【実施例】次に、本発明の実施例を図面を参照して説明
する。図１において、１はデータ入力部であり、金型及
び素材に関する種々のデータが入力され、後述する演算
処理部により金型が所定量移動した場合の素材の流動変
形過程を有限要素法により数値シミュレートして、出力
部３に出力する。
する。図１において、１はデータ入力部であり、金型及
び素材に関する種々のデータが入力され、後述する演算
処理部により金型が所定量移動した場合の素材の流動変
形過程を有限要素法により数値シミュレートして、出力
部３に出力する。
【００１６】次に、図２をも参照して、図１に示された
構成要素の作動内容を説明する。
構成要素の作動内容を説明する。
【００１７】演算処理部は、演算モデル構成部７、速度
場算出部８、パーティクル移動速度算出部９、物理量算
出部１０、パーティクル移動位置算出部１１、及び解析
領域の材料定数、物理量の変更部１２とから構成され、
ＣＰＵで実現される。
場算出部８、パーティクル移動速度算出部９、物理量算
出部１０、パーティクル移動位置算出部１１、及び解析
領域の材料定数、物理量の変更部１２とから構成され、
ＣＰＵで実現される。
【００１８】演算モデル構成部７は、入力データに基づ
いて、解析空間を画定する金型４と、この解析空間に配
される素材５とを含む全解析領域を規定すると共に、解
析空間には、素材５が占める第１の解析領域Ａと、素材
５が存在しない第２の解析領域Ｂとが定義され、これら
第１、第２の解析領域ＡとＢとの両方の領域を解析対象
とする。
いて、解析空間を画定する金型４と、この解析空間に配
される素材５とを含む全解析領域を規定すると共に、解
析空間には、素材５が占める第１の解析領域Ａと、素材
５が存在しない第２の解析領域Ｂとが定義され、これら
第１、第２の解析領域ＡとＢとの両方の領域を解析対象
とする。
【００１９】さらに、演算モデル構成部７は、解析対象
を固定もしくは準固定した複数の分割要素とし、第１の
解析領域Ａに、素材５内の位置情報を規定するパーティ
クル６を分散して設定して有限要素でモデル化して、演
算モデルとする。
を固定もしくは準固定した複数の分割要素とし、第１の
解析領域Ａに、素材５内の位置情報を規定するパーティ
クル６を分散して設定して有限要素でモデル化して、演
算モデルとする。
【００２０】第１の解析領域Ａのエネルギー消散率をΦ
ａ、第２の解析領域Ｂのエネルギー消散率をΦｂとする
と、全解析領域におけるエネルギー消散率Φは、下記の
数式１で表される。
ａ、第２の解析領域Ｂのエネルギー消散率をΦｂとする
と、全解析領域におけるエネルギー消散率Φは、下記の
数式１で表される。
【００２１】第１の解析領域Ａのエネルギー消散率Φａ
は、素材の塑性変形によるエネルギー消散率Φv と、素
材の表面と金型４の面との摩擦仕事によるエネルギー消
散率ΦS とからなり、動的に可容速度場から得られるエ
ネルギー消散率が最小のとき正解と一致する。
は、素材の塑性変形によるエネルギー消散率Φv と、素
材の表面と金型４の面との摩擦仕事によるエネルギー消
散率ΦS とからなり、動的に可容速度場から得られるエ
ネルギー消散率が最小のとき正解と一致する。
【００２２】
【数１】
【００２３】ここで、数式１をさらに具体的に表せば、
次の数式２で表すことができる。
次の数式２で表すことができる。
【００２４】
【数２】
【００２５】なお、実際の計算では第２の解析領域Ｂの
エネルギー消散率Φｂには、Φｂ≦Φａとなる極小の有
限値を与える。また、相当応力として、わずかに素材の
圧縮性を考慮し、体積歪み速度を導入することにより、
ひずみ速度から応力が算出可能となる。
エネルギー消散率Φｂには、Φｂ≦Φａとなる極小の有
限値を与える。また、相当応力として、わずかに素材の
圧縮性を考慮し、体積歪み速度を導入することにより、
ひずみ速度から応力が算出可能となる。
【００２６】次に、演算処理部の具体的な作動内容を、
図４及び図５のフローチャートをも参照して説明する。
図４及び図５のフローチャートをも参照して説明する。
【００２７】演算処理部は、パーティクルを複数の分割
要素上を移動させることによる複数の分割要素の物理量
の変化に基づいて、有限要素法により素材の流動変形過
程を数値シミュレートするものである。
要素上を移動させることによる複数の分割要素の物理量
の変化に基づいて、有限要素法により素材の流動変形過
程を数値シミュレートするものである。
【００２８】速度場算出部８は、有限要素によるモデル
化を行う。すなわち、金型と素材とに基づいて、ある境
界条件の下で、解析領域全体におけるエネルギー消散率
Φが最小となるように、パーティクルが存在する第１の
解析領域Ａにおける分割要素の節点速度ｕa と、パーテ
ィクルが存在しない第２の解析領域Ｂにおける分割要素
の節点速度ｕb とを、解析領域内の速度場として算出す
るものである（フロ−１００）。
化を行う。すなわち、金型と素材とに基づいて、ある境
界条件の下で、解析領域全体におけるエネルギー消散率
Φが最小となるように、パーティクルが存在する第１の
解析領域Ａにおける分割要素の節点速度ｕa と、パーテ
ィクルが存在しない第２の解析領域Ｂにおける分割要素
の節点速度ｕb とを、解析領域内の速度場として算出す
るものである（フロ−１００）。
【００２９】すなわち、下記の数式３に示すように、解
析領域全体におけるエネルギー消散率Φは、それら節点
速度の関数として表せる。
析領域全体におけるエネルギー消散率Φは、それら節点
速度の関数として表せる。
【００３０】
【数３】
【００３１】ここで、節点数ｍ，ｎは解析中、パーティ
クルの移動にともない逐次変化することになる。
クルの移動にともない逐次変化することになる。
【００３２】パーティクル移動速度算出部９は、パーテ
ィクルの移動速度を、内挿関数Ｎｊを用いて、次の数式
４で求めるものである（フロー１１０）。
ィクルの移動速度を、内挿関数Ｎｊを用いて、次の数式
４で求めるものである（フロー１１０）。
【００３３】
【数４】
【００３４】内挿関数Ｎｊは、分割要素ｉの各節点の座
標｛Ｘｉ｝とパーティクルの座標｛Ｘｐ｝を用いる下記
の数式５より求められる。この数式５は、二次の多項式
となるため、実際には、内挿関数Ｎｊはニュートンラプ
ソン法により求められる。
標｛Ｘｉ｝とパーティクルの座標｛Ｘｐ｝を用いる下記
の数式５より求められる。この数式５は、二次の多項式
となるため、実際には、内挿関数Ｎｊはニュートンラプ
ソン法により求められる。
【００３５】
【数５】
【００３６】物理量算出部１０は、速度場算出部８によ
り得られた速度場に基づいて、分割要素の物理量を算出
する（フロー１１０）。すなわち、パーティクルは、速
度のほか、ひずみ、ひずみ速度、及び応力すべての情報
を持ち、その位置でのそれらの値は、数式４と同様に、
それぞれ、下記の数式６にて求められる。
り得られた速度場に基づいて、分割要素の物理量を算出
する（フロー１１０）。すなわち、パーティクルは、速
度のほか、ひずみ、ひずみ速度、及び応力すべての情報
を持ち、その位置でのそれらの値は、数式４と同様に、
それぞれ、下記の数式６にて求められる。
【００３７】
【数６】
【００３８】移動位置算出部１１は、パーティクル移動
速度算出部９より算出された移動速度に基づいて、移動
後におけるパーティクルの位置情報を算出する（フロー
１２０）。すなわち、移動後のパーティクルの位置座標
を、数式４で得られたパーティクルの移動速度とその時
間増分とに基づいて、次の数式７より求めるものであ
る。
速度算出部９より算出された移動速度に基づいて、移動
後におけるパーティクルの位置情報を算出する（フロー
１２０）。すなわち、移動後のパーティクルの位置座標
を、数式４で得られたパーティクルの移動速度とその時
間増分とに基づいて、次の数式７より求めるものであ
る。
【００３９】
【数７】
【００４０】解析領域の材料定数、物理量の変更部１２
は、金型４の移動方向に合わせて全解析領域を変更する
（フロー１３０）。金型移動方向のみを比例的に移動す
れば、容易に解析領域が変更される。一方、図５に示す
ように、解析領域を完全固定した場合は、フロー１３０
の操作は不要である。
は、金型４の移動方向に合わせて全解析領域を変更する
（フロー１３０）。金型移動方向のみを比例的に移動す
れば、容易に解析領域が変更される。一方、図５に示す
ように、解析領域を完全固定した場合は、フロー１３０
の操作は不要である。
【００４１】解析領域の材料定数、物理量の変更部１２
はまた、演算モデル構成部７に新たな分割要素の物理量
及びパーティクルの新座標をフィードバックして、デー
タの変更を行う（フロー１４０）。
はまた、演算モデル構成部７に新たな分割要素の物理量
及びパーティクルの新座標をフィードバックして、デー
タの変更を行う（フロー１４０）。
【００４２】以上の演算処理を繰り返すことにより、解
析が完了する。
析が完了する。
【００４３】ここで、本実施例の鍛造解析の実験例を説
明する。
明する。
【００４４】（第１の実験例）図６に示すように、直径
９０ミリ、高さ６０ミリの円柱を、形状内径１００ミリ
の金型で、４２．５ミリの据え込み中心間距離を１７．
５ミリとする場合の軸対象解析を行った。この場合の分
割要素、図７に示す。対称性を考慮し、１／２モデルと
し、全要素数を１７８１とした。パーティクル表示分が
素材に該当し、材料定数としてアルミ材相当のσ＝１５
２ε0.264 （Ｍｐａ）を用い、摩擦係数をμ＝０．２と
仮定した。
９０ミリ、高さ６０ミリの円柱を、形状内径１００ミリ
の金型で、４２．５ミリの据え込み中心間距離を１７．
５ミリとする場合の軸対象解析を行った。この場合の分
割要素、図７に示す。対称性を考慮し、１／２モデルと
し、全要素数を１７８１とした。パーティクル表示分が
素材に該当し、材料定数としてアルミ材相当のσ＝１５
２ε0.264 （Ｍｐａ）を用い、摩擦係数をμ＝０．２と
仮定した。
【００４５】図８に、解析によって得られた変形図を示
す。据え込み量が、２０ミリ，３０ミリ、４２．５ミリ
（最終）の場合について、それぞれ、パーティクル（図
ａ）と、流線（図ｂ）による表示について示している。
金型４の移動にともない、パーティクルが移動し、素材
が金型４内に充満していく状態が分かる。
す。据え込み量が、２０ミリ，３０ミリ、４２．５ミリ
（最終）の場合について、それぞれ、パーティクル（図
ａ）と、流線（図ｂ）による表示について示している。
金型４の移動にともない、パーティクルが移動し、素材
が金型４内に充満していく状態が分かる。
【００４６】（第２の実験例）次に、本発明の応用例と
して、従来、有限要素解析が不向きとされていたせん断
解析を説明する。
して、従来、有限要素解析が不向きとされていたせん断
解析を説明する。
【００４７】図９に示すように、間隙１ミリの上下の刃
で板厚１０ミリ，直径１９０ミリの円板を直径１００ミ
リの円板に打ち抜く解析を行った。
で板厚１０ミリ，直径１９０ミリの円板を直径１００ミ
リの円板に打ち抜く解析を行った。
【００４８】図１０は、分割要素の状態を示すものであ
り、要素数は１２９６とした。材料定数としては、アル
ミ材相当のσ＝１５２ε0.264 （Ｍｐａ）を用いた。
り、要素数は１２９６とした。材料定数としては、アル
ミ材相当のσ＝１５２ε0.264 （Ｍｐａ）を用いた。
【００４９】図１１に、せん断途中（図ａ）とせん断終
了時（図ｂ）の状態の解析結果を示す。これにより、完
全にせん断された状態までの過程が解析できることが分
かった。
了時（図ｂ）の状態の解析結果を示す。これにより、完
全にせん断された状態までの過程が解析できることが分
かった。
【００５０】（第３の実験例）次に、円柱据え込み解析
の他の実験例を説明する。
の他の実験例を説明する。
【００５１】図１２に示すように、直径Ｄ＝２００ミ
リ、高さ２Ｈｏ＝２００ミリの円柱を高さ１２０ミリま
で据え込みした。
リ、高さ２Ｈｏ＝２００ミリの円柱を高さ１２０ミリま
で据え込みした。
【００５２】図１３に分割要素の状態を示す。対象性を
考慮して１／４モデルとし、変形前の状態で素材（パー
ティクル）が存在する第１の解析領域の要素数を４０
０、存在しない第２の解析領域の要素数を４００、総要
素数を８００（２０×４０）とした。材料定数は、アル
ミ材相当のσ＝１５２ε0.264 （Ｍｐａ）を用い、摩擦
係数をμ＝０．２と仮定した。
考慮して１／４モデルとし、変形前の状態で素材（パー
ティクル）が存在する第１の解析領域の要素数を４０
０、存在しない第２の解析領域の要素数を４００、総要
素数を８００（２０×４０）とした。材料定数は、アル
ミ材相当のσ＝１５２ε0.264 （Ｍｐａ）を用い、摩擦
係数をμ＝０．２と仮定した。
【００５３】解析によって得られた変形形状を図１４に
示す。（図ａ）は圧力量ΔＨ＝２０ミリ（圧力率ΔＨ／
Ｈｏ）＝２０％，（図ｂ）は圧力量ΔＨ＝４０ミリ（圧
力率ΔＨ／Ｈｏ）＝４０％の場合を示す。図面右側は、
流線による表示、図面左側はパーティクルによる表示で
ある。
示す。（図ａ）は圧力量ΔＨ＝２０ミリ（圧力率ΔＨ／
Ｈｏ）＝２０％，（図ｂ）は圧力量ΔＨ＝４０ミリ（圧
力率ΔＨ／Ｈｏ）＝４０％の場合を示す。図面右側は、
流線による表示、図面左側はパーティクルによる表示で
ある。
【００５４】金型４の移動に伴い、パーテイクルが外周
方向へ移動しており、素材の変形状態が明瞭に分かっ
た。
方向へ移動しており、素材の変形状態が明瞭に分かっ
た。
【００５５】図１５は、円柱外周部の変形をプロットし
たものであり、比較のためのリップルス・フォージ（小
坂田・他２名，塑加工論，１９８４−１０，４１）によ
る結果（Ｏ印）と完全に一致している。
たものであり、比較のためのリップルス・フォージ（小
坂田・他２名，塑加工論，１９８４−１０，４１）によ
る結果（Ｏ印）と完全に一致している。
【００５６】（第４の実験例）次に、図１６に示すよう
に、軸対称形状の鍛造解析例の他の実験例を示す。
に、軸対称形状の鍛造解析例の他の実験例を示す。
【００５７】図１７は、分割要素の状態を示す。対象性
を考慮して１／２モデルとし、材料定数は、アルミ材相
当のσ＝１５２ε0.264 （Ｍｐａ）を用い、摩擦係数を
μ＝０．２と仮定した。
を考慮して１／２モデルとし、材料定数は、アルミ材相
当のσ＝１５２ε0.264 （Ｍｐａ）を用い、摩擦係数を
μ＝０．２と仮定した。
【００５８】図１８に、初期状態から、最終状態に至る
までの変形過程を示した。それぞれパーティクルと流線
による変形状態を表示しており、変形が進むにつれ、パ
ーティクルが金型４と接触し、素材が型内に充満してい
く状態が明瞭に分かる。
までの変形過程を示した。それぞれパーティクルと流線
による変形状態を表示しており、変形が進むにつれ、パ
ーティクルが金型４と接触し、素材が型内に充満してい
く状態が明瞭に分かる。
【００５９】以上、本発明を鍛造シミュレーション装置
について説明したが、本発明は射出成形機にも適用でき
ることは言うまでもない。
について説明したが、本発明は射出成形機にも適用でき
ることは言うまでもない。
【００６０】
【発明の効果】以上の説明から分かるように、本発明に
よれば、パーティクルの流れで素材の変形を記述するか
ら、接触の問題を考慮すること無く複雑な形状の鍛造解
析が可能となる。また、解析途中でのリメッシュが不要
となる効果がある。
よれば、パーティクルの流れで素材の変形を記述するか
ら、接触の問題を考慮すること無く複雑な形状の鍛造解
析が可能となる。また、解析途中でのリメッシュが不要
となる効果がある。
【図１】本発明の実施例に係わるブロック図である。
【図２】本発明の実施例に係わる解析領域の概念図であ
る。
る。
【図３】本発明の分割要素とパーティクルとの挙動を示
す概念図である。
す概念図である。
【図４】本発明の実施例に係わる解析領域を所定方向に
移動させた場合のフローチャートである。
移動させた場合のフローチャートである。
【図５】本発明の実施例に係わる解析領域を完全固定さ
せた場合のフローチャートである。
せた場合のフローチャートである。
【図６】型鍛造解析に係わる第１の実験例を示す図であ
る。
る。
【図７】第１の実験例の分割要素を示す図である。
【図８】第１の実験例の結果を、パーティクル（図ａ）
と流線（図ｂ）とで表す変形図である。
と流線（図ｂ）とで表す変形図である。
【図９】せん段解析に係わる第２の実験例を示す図であ
る。
る。
【図１０】第２の実験例に係わる分割要素を示す図であ
る。
る。
【図１１】第２の実験例に係わるせん断途中（図ａ）と
せん断終了時（図ｂ）の状態の解析結果を示す図であ
る。
せん断終了時（図ｂ）の状態の解析結果を示す図であ
る。
【図１２】円柱据え込み解析に係わる第３の実験例を示
す図である。
す図である。
【図１３】第３の実験例に係わる分割要素を示す図であ
る。
る。
【図１４】第３の実験例に係わる円柱据え込み解析結果
を示す図である。
を示す図である。
【図１５】第３の実験例に係わる円柱外周部の変形をプ
ロットした図である。
ロットした図である。
【図１６】軸対称形状の鍛造解析例に係わる第４の実験
例を示す図である。
例を示す図である。
【図１７】第４の実験例に係わる分割要素を示す図であ
る。
る。
【図１８】第４の実験例に係わる軸対称形状の鍛造解析
結果を示す図である。
結果を示す図である。
【図１９】従来技術における金型との接触評価を示す概
念図である。
念図である。
【図２０】従来技術における有限要素のつぶれを示す概
念図である。
念図である。
【図２１】従来技術の型鍛造の有限要素法解析を示す概
念図である。
念図である。
【図２２】従来の有限要素法解析に適用される成形部品
の一例を示した図である。
の一例を示した図である。
【図２３】従来の有限要素法解析における金型及び要素
分割された素材の一例を示した図である。
分割された素材の一例を示した図である。
【図２４】従来の有限要素法解析における金型及び要素
分割された素材の変形過程を示した図である。
分割された素材の変形過程を示した図である。
１ データ入力部 ３ 出力部 ４ 金型 ５ 素材 ６ パーティクル ７ 演算モデル構成部 ８ 速度場算出部 ９ パーティクル移動速度算出部 １０ 物理量算出部 １１ パーティクル移動位置算出部 １２ 解析領域の材料定数、物理量の変更部
Claims (3)
- 【請求項１】 金型内に素材を含む解析空間を規定し
て、前記金型の移動による前記素材の流動変形過程を、
有限要素法により数値シミュレートする成形シミュレー
ション方法において、 前記解析空間には、前記素材が占める第１の解析領域
と、前記素材が存在しない残りの第２の解析領域とを設
定し、前記第１の解析領域と前記第２の解析領域とを解
析対象として有限要素でモデル化して有限要素法により
前記素材の流動変形過程を数値シミュレートし、 前記モデルは、前記解析対象を固定もしくは準固定した
複数の分割要素とし、前記第１の解析領域には、前記素
材形状を規定する複数の識別粒子を分散して設定し、前
記識別粒子を前記解析空間内で移動させ、移動した領域
における材料定数を変化させるようなモデルとすること
を特徴とする成形シミュレーション方法。 - 【請求項２】 請求項１記載の成形シミュレーション方
法において、 前記数値シミュレートは、前記複数の識別粒子を前記複
数の分割要素上を移動させることによる前記複数の分割
要素の物理量の変化に基づいて、有限要素法により前記
素材の流動変形過程をシミュレートするものであること
を特徴とする成形シミュレーション方法。 - 【請求項３】 金型内に素材を含む解析空間を規定し
て、前記金型の移動による前記素材の流動変形過程を、
有限要素法により数値シミュレートする成形シミュレー
ション装置において、 前記解析空間に、前記素材が占める第１の解析領域と、
前記素材が存在しない残りの第２の解析領域とを設定
し、前記第１の解析領域と前記第２の解析領域とを解析
対象として有限要素でモデル化する演算モデル構成部
と、 前記金型と前記素材とに基づいて、予め規定される境界
条件の下で、前記解析空間全体におけるエネルギー消散
率が最小となるように、前記複数の識別粒子が存在する
前記第１の解析領域における前記分割要素の節点速度
と、前記複数の識別粒子が存在しない前記第２の解析領
域における前記分割要素の節点速度とを、前記解析空間
内の速度場として算出する速度場算出部と、 前記速度場に基づいて、前記複数の識別粒子の移動速度
を算出する識別粒子移動速度算出部と、 該識別粒子移動速度算出部より算出された移動速度に基
づいて、移動後における前記複数の識別粒子の位置情報
を算出する移動位置算出部と、 前記速度場に基づいて、前記分割要素の物理量を算出す
る物理量算出部とを有することを特徴とする成形シミュ
レーション装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP29925991A JP2775538B2 (ja) | 1991-11-14 | 1991-11-14 | 成形シミュレーション方法及び装置 |
| US07/976,505 US5402366A (en) | 1991-11-14 | 1992-11-13 | Method and apparatus for simulating a mechanical operation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP29925991A JP2775538B2 (ja) | 1991-11-14 | 1991-11-14 | 成形シミュレーション方法及び装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH05138265A true JPH05138265A (ja) | 1993-06-01 |
| JP2775538B2 JP2775538B2 (ja) | 1998-07-16 |
Family
ID=17870223
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP29925991A Expired - Lifetime JP2775538B2 (ja) | 1991-11-14 | 1991-11-14 | 成形シミュレーション方法及び装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5402366A (ja) |
| JP (1) | JP2775538B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7542817B2 (en) | 1995-02-14 | 2009-06-02 | Aztec Ip Company, L.L.C. | Method and apparatus for manufacturing objects having optimized response characteristics |
Families Citing this family (33)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3207971B2 (ja) * | 1993-06-25 | 2001-09-10 | 富士通株式会社 | 最適骨組及び板組構造の設計方法 |
| JP3501486B2 (ja) * | 1993-12-27 | 2004-03-02 | キヤノン株式会社 | 射出成形品の変形量予測方法及びその装置 |
| JP3127713B2 (ja) * | 1994-05-10 | 2001-01-29 | 株式会社村田製作所 | 粉末成形の最適条件解析方法 |
| US5572434A (en) * | 1994-06-14 | 1996-11-05 | Cornell Research Foundation, Inc. | Method for simulating mold filling of semi-solid material |
| US5581468A (en) * | 1995-01-11 | 1996-12-03 | The Dow Chemical Co. | System and method for modeling plastic molding and molding parts incorporating the same |
| US5940309A (en) * | 1996-09-06 | 1999-08-17 | White; Warren D. | System and method for modeling plastic molding and molding parts incorporating the same |
| US5754447A (en) * | 1996-10-30 | 1998-05-19 | Sandia National Laboratories | Process for predicting structural performance of mechanical systems |
| US6125333A (en) * | 1997-11-06 | 2000-09-26 | Northrop Grumman Corporation | Building block approach for fatigue spectra generation |
| JP2000197942A (ja) * | 1998-10-30 | 2000-07-18 | Toyota Motor Corp | 鍛造加工解析方法およびその実行プログラムを記録した媒体 |
| US6731996B1 (en) | 1998-11-18 | 2004-05-04 | Alcan International Limited | Method of using isotropic plasticity to model the forming of anisotropic sheet |
| BR9915486A (pt) * | 1998-11-18 | 2001-07-31 | Alcan Int Ltd | Método para modelar a formação de chapas anisotrópicas |
| US6402367B1 (en) * | 2000-06-01 | 2002-06-11 | Aemp Corporation | Method and apparatus for magnetically stirring a thixotropic metal slurry |
| US6432160B1 (en) * | 2000-06-01 | 2002-08-13 | Aemp Corporation | Method and apparatus for making a thixotropic metal slurry |
| US6399017B1 (en) * | 2000-06-01 | 2002-06-04 | Aemp Corporation | Method and apparatus for containing and ejecting a thixotropic metal slurry |
| US7024342B1 (en) | 2000-07-01 | 2006-04-04 | Mercury Marine | Thermal flow simulation for casting/molding processes |
| US6611736B1 (en) | 2000-07-01 | 2003-08-26 | Aemp Corporation | Equal order method for fluid flow simulation |
| US7010474B1 (en) * | 2000-11-03 | 2006-03-07 | General Electric Company | Performance predictor for fabricated parts |
| US20020199156A1 (en) * | 2001-05-23 | 2002-12-26 | Fuel Tech, Inc. | Hardware-adaptable data visualization tool for use in complex data analysis and engineering design |
| US6829563B2 (en) * | 2001-05-31 | 2004-12-07 | Sumitomo Rubber Industries, Ltd. | Simulation method for estimating performance of product made of viscoelastic material |
| US6742567B2 (en) * | 2001-08-17 | 2004-06-01 | Brunswick Corporation | Apparatus for and method of producing slurry material without stirring for application in semi-solid forming |
| FR2847999B1 (fr) * | 2002-12-02 | 2005-02-11 | E S I Software | Logiciel de modelisation d'emboutissage |
| US7275446B2 (en) * | 2003-07-23 | 2007-10-02 | Engis Corporation | Apparatus for measuring the crushing strength of micron superabrasives |
| DE602004023360D1 (de) * | 2003-10-17 | 2009-11-12 | Sumitomo Rubber Ind | Verfahren zur Simulation von viskoelastischem Material |
| US20040194051A1 (en) * | 2004-05-13 | 2004-09-30 | Croft Bryan L. | Finite element modeling system and method for modeling large-deformations using self-adaptive rezoning indicators derived from eigenvalue testing |
| KR101383663B1 (ko) * | 2005-10-04 | 2014-04-09 | 아즈텍 아이피 컴퍼니, 엘.엘.씨. | 가상 테스트에 근거한 파라미터화 재료 및 성능 특성 |
| US7945344B2 (en) * | 2008-06-20 | 2011-05-17 | SAKT13, Inc. | Computational method for design and manufacture of electrochemical systems |
| US9249502B2 (en) * | 2008-06-20 | 2016-02-02 | Sakti3, Inc. | Method for high volume manufacture of electrochemical cells using physical vapor deposition |
| US8357464B2 (en) | 2011-04-01 | 2013-01-22 | Sakti3, Inc. | Electric vehicle propulsion system and method utilizing solid-state rechargeable electrochemical cells |
| US10770745B2 (en) | 2011-11-09 | 2020-09-08 | Sakti3, Inc. | Monolithically integrated thin-film solid state lithium battery device having multiple layers of lithium electrochemical cells |
| US8301285B2 (en) | 2011-10-31 | 2012-10-30 | Sakti3, Inc. | Computer aided solid state battery design method and manufacture of same using selected combinations of characteristics |
| US9127344B2 (en) | 2011-11-08 | 2015-09-08 | Sakti3, Inc. | Thermal evaporation process for manufacture of solid state battery devices |
| US9627717B1 (en) | 2012-10-16 | 2017-04-18 | Sakti3, Inc. | Embedded solid-state battery |
| US9627709B2 (en) | 2014-10-15 | 2017-04-18 | Sakti3, Inc. | Amorphous cathode material for battery device |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3896666A (en) * | 1974-03-25 | 1975-07-29 | Elmer G Johnson | Induced agglomeration system for use in a system for simulating reentry flow conditions |
| US4567774A (en) * | 1983-04-28 | 1986-02-04 | Battelle Development Corporation | Determining mechanical behavior of solid materials using miniature specimens |
| JPH07104855B2 (ja) * | 1985-03-28 | 1995-11-13 | インターナショナル・ビジネス・マシーンズ・コーポレーション | 数値シミュレーション装置 |
| JPH0193474A (ja) * | 1987-02-26 | 1989-04-12 | Ulvac Corp | セラミツクス接合方法 |
| US4895027A (en) * | 1988-01-19 | 1990-01-23 | Battelle Development Corporation | Determining plane strain fracture toughness and the J-Integral for solid materials using stress field modified miniature specimens |
| US4933889A (en) * | 1988-04-29 | 1990-06-12 | International Business Machines Corporation | Method for fine decomposition in finite element mesh generation |
| US5136497A (en) * | 1990-07-12 | 1992-08-04 | Bdm International, Inc. | Material consolidation modeling and control system |
| US5202837A (en) * | 1990-07-12 | 1993-04-13 | Coe Carlos J | Material consolidation modeling and control system |
-
1991
- 1991-11-14 JP JP29925991A patent/JP2775538B2/ja not_active Expired - Lifetime
-
1992
- 1992-11-13 US US07/976,505 patent/US5402366A/en not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7542817B2 (en) | 1995-02-14 | 2009-06-02 | Aztec Ip Company, L.L.C. | Method and apparatus for manufacturing objects having optimized response characteristics |
Also Published As
| Publication number | Publication date |
|---|---|
| US5402366A (en) | 1995-03-28 |
| JP2775538B2 (ja) | 1998-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH05138265A (ja) | 成形シミユレーシヨン方法及び装置 | |
| Lee et al. | A new method of preform design in hot forging by using electric field theory | |
| CN1065175C (zh) | 分析流体流动过程的装置和方法 | |
| Koç et al. | The use of FEA and design of experiments to establish design guidelines for simple hydroformed parts | |
| Kulon et al. | Applying knowledge-based engineering to traditional manufacturing design | |
| Lin et al. | Application of abductive network and FEM to predict an acceptable product on T-shape tube hydroforming process | |
| JPH10301975A (ja) | 製品の設計パラメータ最適値決定方法、システム及びその方法プログラムを記載した記録媒体 | |
| EP4066144A1 (en) | Optimization of geometry of shaped body and manufacturing tools | |
| Choi et al. | Application of intelligent design support system for multi-step deep drawing process | |
| Vemuri et al. | BID: a knowledge-based system to automate blocker design | |
| Price et al. | Deciding parameter values with case-based reasoning | |
| Alfozan et al. | An upper bound elemental technique approach to the process design of axisymmetric forging by forward and backward simulation | |
| Lewis et al. | Alternative techniques for casting process simulation | |
| Rathod et al. | Optical CAD modelling and designing of compound die using the python scripting language | |
| Samolyk et al. | Application of the slip-line field method to the analysis of die cavity filling | |
| Jirathearanat et al. | Virtual processing–application of rapid prototyping for visualization of metal forming processes | |
| Lin et al. | Radius ratio estimation and fold situation prediction of the deformation profile in forging–extrusion process | |
| Naranje et al. | An intelligent CAD system for automatic modelling of deep drawing die | |
| CN113343523A (zh) | 一种基于空间网格的有限元数值模拟分析方法 | |
| JP2006516932A (ja) | プレスをモデル化するためのソフトウェア | |
| JP3006374B2 (ja) | 鍛造用型打方向算出装置 | |
| JP2653217B2 (ja) | 塑性加工用材料定数の設定方法 | |
| Kwan et al. | Die shape design for T-shape tube hydroforming | |
| JPH10269265A (ja) | 有限要素法による構造解析方法とその装置 | |
| JP2002052560A (ja) | 射出成形品製造パラメータ決定支援システム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19980401 |