JP2004342010A

JP2004342010A - シミュレーション装置、シミュレーション方法、およびシミュレーションプログラム

Info

Publication number: JP2004342010A
Application number: JP2003140431A
Authority: JP
Inventors: Takashi Moriya; 岳志守屋; Rie Akazawa; 理恵赤澤
Original assignee: Nissan Motor Co Ltd
Current assignee: Nissan Motor Co Ltd
Priority date: 2003-05-19
Filing date: 2003-05-19
Publication date: 2004-12-02

Abstract

【課題】メッシュの自己交差の発生しないワンステップ解法によるシミュレーション方法を提供する。
【解決手段】ワンステップ解法によるシミュレーションを実行し（Ｓ２）、平面展開形状のアウトラインを抽出して（Ｓ３）、自己交差がある場合（Ｓ４：Ｙｅｓ）、自己交差部分の形状を変更して（Ｓ５）、変更量に応じて端部の拘束力を補正し（Ｓ６）、シミュレーションを行う。これを自己交差が無くなるまで実行する。
【選択図】図６

Description

【０００１】
【発明の属する技術分野】
本発明は、シミュレーション装置、シミュレーション方法、およびシミュレーションプログラムに関し、特にプレス形成品などのワーク形状を、ワンステッブ解法によりシミュレーションするシミュレーション装置、シミュレーション方法、およびシミュレーションプログラムに関する。
【０００２】
【従来の技術】
従来、プレス部品のドロー成形工程のシミュレーション方法として、ワンステッブ解法を用いた技術がある（たとえば、非特許文献１参照）。
【０００３】
ここでワンステップ解法とは、部品形状（ワーク形状）から計算を始めて、部品形状と展開形状の歪エネルギー差が最小となるような平面展開形状を弾塑性の収束計算により求め、部品形状と展開形状のメッシュ寸法の変更結果から応力および歪を１ステップで計算するシミュレーション方法である。したがって、増分法によるシミュレーションと比較すれば、精度的には劣るもののシミュレーション時間は大幅に少なくてすみ、さまざまな形状変化などをシミュレーションして試すのに適している。また、精度が劣るといっても部品形状によっては十分使用できる。
【０００４】
【非特許文献１】
矢野ほか、「ＯｎｅＳｔｅｐＦＥＭのブランク形状予測精度の評価」第４９回塑性加工連合講演会論文集、第１７３〜１７４頁、平成１０年１０月３０日、社団法人日本塑性加工学会発行
【０００５】
【発明が解決しようとする課題】
しかしながら、従来のワンステップ解法は、その解法上、部品形状を平面に展開して、ワークであるプレス部品のブランク形状を求めるため、部品アウトラインが小さな曲率を有する部位では、メッシュが自己交差してしまう。
【０００６】
図７は、このようなメッシュの自己交差を説明するための図面である。
【０００７】
図において図７Ａは、自動車の車体パネルの一部であり、プレスにより形成されるワークの一例である。図７Ｂは、ワンステップ解放によりシミュレーションを行った結果であり、図７Ａ中円で囲んだ部分を示す。
【０００８】
図７Ｂに示すように、シミュレーションの際に用いるメッシュのアウトライン１０１に自己交差が発生している。
【０００９】
このような自己交差が発生すると、シミュレーションの結果においてこの自己交差発生部分の板厚が極端に減少して、正しい評価ができなくなるという問題がある。また、実際のプレス部品では、このような自己交差することはありえないために、その点でもシミュレーション結果として好ましくない。
【００１０】
そこで、本発明の目的は、メッシュの自己交差の発生しないワンステップ解法によるシミュレーション装置、シミュレーション方法、およびシミュレーションプログラムを提供することである。
【００１１】
【課題を解決するための手段】
上記目的は、ワーク形状をワンステップ解法によりシミュレーションするシミュレーション実行手段と、前記シミュレーション実行手段により行われたシミュレーション結果の形状からアウトラインを抽出するアウトライン抽出手段と、前記アウトライン抽出手段が抽出したアウトラインからメッシュの自己交差の有無を検出する自己交差検出手段と、前記自己交差検出手段が自己交差を検出した部分の形状を自己交差が発生しない形状に変形する部品形状変形手段と、を有することを特徴とするシミュレーション装置により達成される。
【００１２】
上記目的は、ワーク形状をワンステップ解法によりシミュレーションする段階と、前記シミュレーションにより得られた前記ワーク形状からアウトラインを抽出する段階と、前記アウトラインからメッシュの自己交差の有無を検出する段階と、前記自己交差を検出した部分の形状を自己交差が発生しない形状に変形する段階と、を有することを特徴とするシミュレーション方法により達成される。
【００１３】
上記目的は、ワーク形状をワンステップ解法によりシミュレーションするステップと、前記シミュレーションにより得られた形状からアウトラインを抽出するステップと、前記アウトラインからメッシュの自己交差の有無を検出するステップと、前記自己交差を検出した部分の形状を自己交差が発生しない形状に変形するステップと、前記変形された部分の拘束力を変形量に応じて補正するステップと、を有することを特徴とするシミュレーションプログラムにより達成される。
【００１４】
【発明の効果】
本発明によれば、ワンステップ解法を用いたシミュレーションにおいて解析不具合となる自己交差を検出し、自己交差が発生する場合には、その部分の形状を変形して自己交差が起こらないようにしたので、シミュレーションの精度を維持した最適部品形状を得ることができる。これにより、ワンステッブ解析のやりなおし解析の工数削減や解析精度の向上を図ることができる。
【００１５】
【発明の実施の形態】
以下、図面を参照して本発明の実施の形態を説明する。
【００１６】
図１は、本発明を適用した実施の形態におけるシミュレーション装置の機能構成を説明するためのブロック図である。
【００１７】
このシミュレーション装置１の機能構成は、図１に示すように、ワーク形状のシミュレーションを行うシミュレーション部１０と、シミュレーションされたワーク形状からアウトラインを抽出するアウトライン抽出部１１と、メッシュの自己交差の発生の有無を評価する自己交差評価部１２と、自己交差部の部品形状を変更する自己交差部形状変更部１３と、部品端部の拘束力を補正する部品端拘束力補正部１４と、からなる。
【００１８】
また、シミュレーション装置１には、たとえばＣＡＤ装置２０などからシミュレーションを行うワークの部品形状データを入力するための部品形状データ入力部２１と、部品端部における拘束力（詳細後述）を入力および出力するための部品端拘束力入出力部２２と、シミュレーション結果を出力する出力装置２３が接続されている。なお、ここで部品端拘束力入出力部２２は、図示しない外部記憶装置や他のコンピュータなどに接続されており、後述する拘束力関係式などを取り入れたりまたは外部記憶装置などに対して書き出したりするためのインターフェースである。
【００１９】
出力装置２３は、シミュレーション結果を表示するためのディスプレイやプリンタ、およびシミュレーション結果を記憶するための記憶媒体などである。
【００２０】
また、このシミュレーション装置１には、図示しないキーボードやポインティングデバイスなどの入力装置が接続されている。
【００２１】
なお、このようなシミュレーション装置１は、実際にはコンピュータであり、各部の機能はこのコンピュータによって後述する手順を処理するため作成されたプログラムが実行することにより提供される。
【００２２】
以下、各部の機能について説明する。
【００２３】
シミュレーション部１０は、シミュレーション実行手段であり、すでに従来技術として説明したワンステップ解法により、ワークであるプレス形成品の形状をシミュレーションする。
【００２４】
アウトライン抽出部１１は、アウトライン抽出手段であり、シミュレーション部１０において、行われたシミュレーションの結果である部品形状データの平面展開時に、その平面展開形状からアウトラインを抽出する。このようなアウトラインの抽出は、従来からある手法を用いればよく、既存のアウトライン抽出方法としては、たとえば、フリーエッジ抽出法などが用いられる。このような抽出方法は多くのシミュレーションプログラムやＣＡＥプログラムなどに機能として組み込まれているので、実際には、それらの機能を利用することになる。このようなフリーエッジ抽出法によるアウトライン抽出機能をもつシミュレーションプログラムとしての一例をあげれば、「ＭＳＣ．ＶｉｓｕａｌＮａｓｔｒａｎｆｏｒＷｉｎｄｏｗｓ（ただしＷｉｎｄｏｗｓは登録商標）」（株式会社ＣＲＣソリューションズ製）がある。
【００２５】
自己交差評価部１２は、自己交差検出手段であり、アウトライン抽出部１１によって抽出されたアウトラインから、シミュレーションの際に用いているメッシュの線が交差してしまう部分（自己交差）を検出する。自己交差部分の検出も既存の技術を使用すればよく、既存のシミュレーションプログラムに搭載されている自己交差検出機能をそのまま用いればよい。このような自己交差の検出機能を有するプログラムの一例を挙げれば、ＧＨＬ（株式会社精密形状処理研究所製）がある。
【００２６】
自己交差部形状変更部１３は、部品形状変形手段であり、自己交差評価部１２において検出された自己交差部分を、自己交差が発生しない臨界曲率となるように、自己交差発生部分（主にＲ部分）の形状を変更する。
【００２７】
この自己交差部分の形状変更について説明する。
【００２８】
まず、ワンステップ解法を用いたシミュレーションでは、展開時の部品端形状において、メッシュが自己交差している場合、―般に展開方向から見た部品端形状は小さなＲをもった凹形状を有することとなる。
【００２９】
この自己交差の発生を防ぐために部品端形状の凹形状を大きな曲率のＲ形状へと変更することで、同等の材料における端部の移動量で比較した場合、メッシュが自己交差する可能性は低くなる。しかし、大きな曲率のＲ形状への変更は部品形状を大きく変更してしまうため、解析精度には悪い影響を与えてしまう。また、部品端の変形に対する抵抗力が変わってしまうため、Ｒ形状の変更は必要最小限にとどめる必要がある。すなわち、自己交差が発生しない臨界曲率のＲ形状となるように部品端形状に変更する必要がある。
【００３０】
そこで、本実施の形態では、自己交差部の部品形状変更に、最適化手法としてはもっともシンブルな二分法を用いた。この手法は一定増分でパラメータを増やし、対象となる目的事象がＮＧからＯＫ、またはＯＫからＮＧとなった場合に、探索方向を逆転し、増分を半減し、再度探索する、ことを繰り返す手法である。そして、増分があらかじめ指定した値に連したら部品端形状は最適化されたとみなす。
【００３１】
図２は、この部品形状の最適化を説明するための図面である。この例では、初期Ｒ増分がその部分の曲率半径として１０ｍｍ、最小増分が２．５ｍｍとした。
【００３２】
図２中、自己交差の有無を、自己交差発生無しをＯＫ、自己交差発生有りをＮＧとして表している。
【００３３】
図において、反復数＝０のときが実際の部品形状で、Ｒ５０（Ｒ＝５０ｍｍ）としたときのシミュレーション結果である。この状態では、自己交差が発生してＮＧとなる。そこで、Ｒ５０（Ｒ＝５０ｍｍ）からＲ６０（Ｒ＝６０ｍｍ）に曲率半径を１０ｍｍ増やすと、まだ、自己交差が発生してＮＧとなる。さらにＲ６０（Ｒ＝６０ｍｍ）からＲ７０（Ｒ＝７０ｍｍ）に曲率半径を１０ｍｍ増やすと自己交差が無くなり、ＯＫとなるが、これが臨界値どうかわからないので、今度は逆にＲ７０（Ｒ＝７０ｍｍ）からＲ６５（Ｒ＝６５ｍｍ）へ曲率半径を５ｍｍ減らす。そうすると、また、ＮＧとなるので、Ｒ６５（Ｒ＝６５ｍｍ）からＲ６７．５（Ｒ＝６７．５ｍｍ）と曲率半径を増加させる値を先の半分の２．５ｍｍにすると、ＯＫとなる。
【００３４】
このようにして曲率半径の値を移行してゆくことで、部品端形状の大きさが最適化されていく。このような手法を用いることで、最小限の形状変更で解析不具合を解消することができる。
【００３５】
なお、ここでは二分法を用いて臨界曲率となるような変形量を求めたが、これは、当然に他の方法でもよく、たとえば、曲率半径を少しずつ増加させて行き、自己交差が発生しなくなった値とするようにしてもよし、逆に自己交差が発生しないことが推測される大きな値から少しずつ小さな値に減少させて、自己交差が発生したところでその前の値とするようにしてもよい。
【００３６】
このように部品端形状を変えることで、自己交差の発生しないワンステップ解放によるシミュレーションを行うことができる。しかし、これだけではワンステップ解放によるシミュレーションでは、部品端形状を変えたことによる素材の伸び縮みに対する影響を無視してしまうことになる。そこで、本実施の形態では、さらに、自己交差発生部において実際のプレス時の状態を再現できるような拘束力の値に補正している。なお、この拘束力は、プレス時におけるブランク材のどこをどの程度の力で拘束するかを規定する値である。この拘束力の補正は、拘束力補正手段として機能する部品端拘束力補正部１４によって行われる。
【００３７】
部品端拘束力補正部１４による部品端拘束力の補正について説明する。
【００３８】
自己交差部形状変更部１３によって部品端形状を変更すると、部品展開時に部品形状そのものが持つ変形抵抗も変化する。そのため、設計部品形状で標準化された部品端拘束力は、形状の変更に応じて変更する必要がある。
【００３９】
図３は、部品端形状の変更と拘束力の関係を説明するための図面である。
【００４０】
図３において、部品は、自動車車体のサイドパネル部品であり、図示する部分は、そのコーナー部分拡大図である。
【００４１】
図示するように、設計値におけるＲ形状の曲率半径ｒａが、前記の手法によりｒｂに変更された場合を考える。設計形状での部品端拘束力標準値をＦａとしたとき、Ｒ形状変形後の拘束力Ｆｂを求める。実際にはコーナー角度や部品の成形深さ等多くの要因により変動するため、初等式等で表示したものを使用するのは不十分である。
【００４２】
このため、本実施の形態では基礎実験により得られたデータをもとに部品端拘束力Ｆの補正を実施する。
【００４３】
この基礎実験は、Ｒ形状の曲率半径を変えたときに、その部分の板厚が所定の許容範囲内で同じになるような拘束力を、拘束力を変えてシミュレーションすることにより求める。
【００４４】
図４および図５は、拘束力補正のための基礎実験の結果を示すグラフである。
【００４５】
ここでは、一例として、Ｒ形状を、その曲率半径がＲ５０（Ｒ＝５０ｍｍ）、Ｒ１００（Ｒ＝１００ｍｍ）、Ｒ１５０（Ｒ＝１５０ｍｍ）、Ｒ２００（Ｒ＝２００ｍｍ）、Ｒ２５０（Ｒ＝２５０ｍｍ）の、５ポイントにおいてそれぞれ試行錯誤により拘束力を変えて板厚Δｔが同じになる拘束力を求めたものである。
【００４６】
図４に示したグラフは、これらのＲ形状に対してこの板厚Δｔが所定の許容範囲内で同じ減少率となる拘束力Ｆ（シミュレーションプログラムで使用する値であり任意単位である）をプロットした結果である。なお、板厚Δｔは、プレス前のブランク材の板厚から、プレス後においてこのＲ形状部分の板厚が約３０％減少した板厚となるように設定している。
【００４７】
実際の値は、Ｒ５０のとき拘束力０で板厚減少率３０．４％、Ｒ１００のとき拘束力０．３３で板厚減少率３０．２％、Ｒ１５０のとき拘束力０．４３で板厚減少率３０．６％、Ｒ２００のとき拘束力０．５６で板厚減少率３０．０％、Ｒ２５０のとき拘束力０．６０で板厚減少率３０．５％であった。
【００４８】
ここで、目標とする板厚の減少率は、約３０％としているが、これは、自動車車体に通常使用される鋼板をプレスした場合に、経験的にこれ以上板厚が減少するような場合には、プレスによりブランク材に破断が生じる可能性があるためである。したがって、用いるブランク材の材質や厚みなどにより、基礎実験において設定する板厚は適宜変更する必要がある。
【００４９】
また、目標とする板厚の減少率に対して、約±１％程度の誤差は所定の許容範囲とした。もちろんこの許容範囲については、求めるシミュレーション精度、ブランク材の材質、および板厚などに応じて適宜決定すればよく限定されるものではない。これはたとえば、シミュレーション精度が低下してもよい場合、また、ブランク材が破断してしまう減少率を２５％としたとき、実際にはブランク材の材質や板厚によって３０％程度までは誤差として許容できる（破断しない）ような場合などには、上記所定の許容範囲を±５％程度までとしてもよいのである。
【００５０】
一方、図５は、Ｒ形状を変えたときに、拘束力を変えずにシミュレーションを行った場合の板厚Δｔ（減少率（％））の変化を示すグラフである。
【００５１】
この基礎実験の結果、図４に示すように、拘束力を変えることで、Ｒ形状の曲率半径を変化させた場合でも板厚をほぼ同じにできる拘束力を得ることができることがわかる。一方、図５から、拘束力を変えずにＲ形状を変化させた場合は、板厚も変わってしまうことがわかる。
【００５２】
基礎実験の結果から得られたＲ形状の曲率半径を変更した場合に対応して設定する拘束力は、図４に示したグラフから、曲率半径に対応して同じ板厚Δｔの減少率となる拘束力を求めるための拘束力関係式を求めて、この関係式の係数を使用することになる。
【００５３】
上記基礎実験で得られる関係式は、板厚をｙ、曲率半径をｘとすると、
ｙ＝０．３７４３Ｌｎ（ｘ）−１．４３８８ …（１）
となる。なお、この関係式によって得られる曲線を図４中に示した。
【００５４】
このように基礎実験によって、あらかじめ変形量を変えたときに板厚の減少率が決められた範囲内で一定となる変形量と拘束力との関係式を求めることで、自己交差をなくすためにＲ形状部分の曲率半径を変えた際に、それに対応した拘束力を設定することができ、Ｒ形状変更後のシミュレーション結果における板厚分布をより正確に求めることができるようになる。なお、目標とする板厚の減少率は、あらかじめ決めた範囲内に収まるように（上記の基礎実験では、ブランク材が破断しないような板厚の範囲）、ブランク材の材質や元の厚さ、また求めるシミュレーションの精度などから決定すればよく、限定されない。
【００５５】
なお、基礎実験によって求めた拘束力関係式は、データベース（不図示）としてブランク材の種類や板厚などによりあらかじめ数種類準備しておくことで、さまざまなブランク材の種類や、元の板厚などが変更になった際に、どのように部品端拘束力を変更すればいいのかを決定することができるようになる。
【００５６】
したがって、部品端拘束力補正部１４は、このようなデータベースを参照して部品端拘束力の補正を行う。
【００５７】
なお、基礎実験自体も、このシミュレーション装置１を使用して行うことができる。そして基礎実験により求められた結果である拘束力の関係式は、部品端拘束力入出力部２２を介して、外部の記憶装置にデータベースとして記憶させておく。
【００５８】
次に、このように構成されたシミュレーション装置によるシミュレーションの手順を説明する。
【００５９】
図６は、シミュレーション装置によるシミュレーションの手順を示すフローチャートである。
【００６０】
まず、シミュレーション装置１は、部品形状データ入力部２１から部品形状データを取り込むと共に、部品端拘束力入出力部２２からシミュレーションを行う際の標準拘束力を設定する（Ｓ１）。ここで、標準拘束力は、たとえばこれまでに行ったプレス成形のシミュレーションにおいて使用した値を部品形状や部品名と共にデータベースとして記憶しておき、シミュレーションする部品と形状や名称が同じものがあればそれをそのまま使用すればよい。
【００６１】
続いて、シミュレーション部１０がワンステップ解法によりシミュレーションを実行する（Ｓ２）。このシミュレーションは市販のワンステップ解法シミュレーションプログラムによって行われる。このシミュレーションは、部品の形状を展開することにより、プレスの際に必要となるブランク材の形状を求めたり、また、ブランク材からの板厚変化を求めたり、さらに、プレス時にしわや破断などが発生しないかをシミュレーションするものである。
【００６２】
続いて、アウトライン抽出部１１がシミュレーション結果から得られた平面展開形状からアウトラインを抽出する（Ｓ３）。このアウトライン抽出により得られたアウトラインがプレス前の部品の展開形状となる。
【００６３】
続いて、自己交差評価部１２が自己交差の発生を評価する（Ｓ４）。その結果、自己交差の発生があれば（Ｓ４：Ｙｅｓ）、自己交差部部品形状変更部１３が、自己交差が発生した部分の形状を変更する（Ｓ５）。
【００６４】
そして、部品端拘束力補正部１４が拘束力の補正を行って（Ｓ６）、ステップＳ２へ戻り、ワンステップ解法によるシミュレーションを実行する。
【００６５】
したがって、ステップＳ４において自己交差が検出されると、ステップＳ２からＳ６の処理を自己交差が無くなるまで実行することになる。
【００６６】
一方、ステップＳ４の自己交差評価において、自己交差の発生が検出されなければ（Ｓ４：Ｎｏ）、そのままシミュレーション結果を出力装置２３に出力（Ｓ７）および表示（Ｓ８）して（印刷や記憶動作などを含む）、処理を終了する。
【００６７】
以上のように本実施の形態により、自己交差が発生するような曲率半径を有する部分でも板厚の極端な減少や実際の部品であり得ない自己交差した形状結果などがなくなり、ワンステップ解法を用いたシミュレーションであっても、より精度の高いシミュレーション結果を得ることができるようになる。
【００６８】
また、本実施の形態においては、自己交差部分の形状を変更する際に二分法を用いたことで、形状変更にかかる時間も少なく、最小限の形状変更で解析不具合を解消することができる。
【００６９】
さらに、本実施の形態においては、自己交差が発生する部品端部におけるＲ形状を単に変更しただけでなく、変更したＲ形状に合わせて拘束力を補正することとしたので、変更後の形状でもシミュレーション精度が低下することがない。
【００７０】
そして、この拘束力は、基礎実験によりさまざまなブランク材の種類や元の板厚、目標とするプレス後の板厚などに応じてデータベース化することで、さまざまなプレスのシミュレーションに適応可能となる。
【００７１】
以上本発明の実施の形態を説明したが、本発明はこのような実施の形態に限定されるものではない。たとえば、拘束力の補正については、Ｒ形状の曲率半径を変えることによって、板厚が大きく変化しないような場合、または部品端から離れた部位の形状についてのみ求めたいような場合には、拘束力の補正は行わなくてもよい。また、拘束力を補正するための関係式についても、シミュレーションごとに入力するようにしてもよく、必ずしもデータベース化しておく必要はない。
【００７２】
さらには、当業者によるさまざまな改変が可能であることは言うまでもない。
【図面の簡単な説明】
【図１】本発明を適用した実施の形態におけるシミュレーション装置の機能構成を説明するためのブロック図である。
【図２】この部品形状の最適化を説明するための図面である。
【図３】部品端形状の変更と拘束力説明する図面である。
【図４】拘束力補正のための基礎実験の結果を示すグラフであり、曲率半径を変えたときに板厚が同じ減少率となる拘束力をプロットものである。
【図５】Ｒ形状を変えたときに、拘束力を変えずにシミュレーションを行った場合の板厚の減少率を示すグラフである。
【図６】シミュレーションの手順を示すフローチャートである。
【図７】自己交差を説明するための図面である。
【符号の説明】
１…シミュレーション装置
１０…シミュレーション部
１１…アウトライン抽出部
１２…自己交差評価部
１３…自己交差部形状変更部
１４…部品端拘束力補正部

Claims

ワーク形状をワンステップ解法によりシミュレーションするシミュレーション実行手段と、
前記シミュレーション実行手段により行われたシミュレーション結果の形状からアウトラインを抽出するアウトライン抽出手段と、
前記アウトライン抽出手段が抽出したアウトラインからメッシュの自己交差の有無を検出する自己交差検出手段と、
前記自己交差検出手段が自己交差を検出した部分の形状を自己交差が発生しない形状に変形する部品形状変形手段と、
を有することを特徴とするシミュレーション装置。
前記部品形状変形手段によって変形された部分の拘束力を変形量に応じて補正する拘束力補正手段を、さらに有することを特徴とする請求項１記載のシミュレーション装置。
前記拘束力補正手段は、あらかじめ変形量を変えたときに前記ワークの板厚の減少率が所定の許容範囲内で同じとなる変形量と拘束力との関係式を求めておき、該関係式に基づいて前記拘束力を補正することを特徴とする請求項２記載のシミュレーション装置。
ワーク形状をワンステップ解法によりシミュレーションする段階と、
前記シミュレーションにより得られた前記ワーク形状からアウトラインを抽出する段階と、
前記アウトラインからメッシュの自己交差の有無を検出する段階と、
前記自己交差を検出した部分の形状を自己交差が発生しない形状に変形する段階と、
を有することを特徴とするシミュレーション方法。
前記自己交差を検出した部分の形状を自己交差が発生しない形状に変形する段階の後に、前記変形された部分の拘束力を変形量に応じて補正する段階を、さらに有することを特徴とする請求項４記載のシミュレーション方法。
前記変形された部分の拘束力を変形量に応じて補正する段階は、あらかじめ変形量を変えたときに前記ワークの板厚の減少率が所定の許容範囲内で同じとなる変形量と拘束力との関係式を求めておき、該関係式に基づいて前記拘束力を補正することを特徴とする請求項５記載のシミュレーション装置。シミュレーション装置。
ワーク形状をワンステップ解法によりシミュレーションするステップと、
前記シミュレーションにより得られた形状からアウトラインを抽出するステップと、
前記アウトラインからメッシュの自己交差の有無を検出するステップと、
前記自己交差を検出した部分の形状を自己交差が発生しない形状に変形するステップと、
前記変形された部分の拘束力を変形量に応じて補正するステップと、
を有することを特徴とするシミュレーションプログラム。