JP6794823B2 - 解析装置、解析方法、及びコンピュータプログラム - Google Patents

Description

本発明は、解析装置、解析方法、及びコンピュータプログラムに関するものである。

衝突安全性と軽量化の要請から、自動車の車体等に適用されている鋼板の引張強さは、５９０ＭＰａ級を超え、ＧＰａの単位にまで達しつつある。高強度鋼板は、板厚を増加させることなく、衝突時の吸収エネルギー、強度等を高めることができる。他の軽量化素材に比べれば、高強度鋼板を用いた部品加工、組み立て等は、必ずしも設備や生産技術の大きな変革を要しない。従って、高強度鋼板を用いた部品加工、組み立て等の生産コストへの負荷は、他の軽量化素材に比べれば、比較的小さいと考えられている。

しかしながら、プレス加工の場合、鋼板の強度上昇と共に、スプリングバック、しわ等が増加し、部品の寸法精度の確保が困難となる。また、強度上昇に伴う延性の低下は、プレス成形時の破断の危険性を高める。このように、高強度鋼板による性能と生産性の両立は、従来の軟鋼板に比べれば必ずしも容易ではなく、開発工期短縮、製造コスト抑制等と相まって、量産までの加工技術開発の負荷を大きく高めている。

上述の破断、しわ、及びスプリングバックのうち、破断及びしわは、部品形状、成形条件等によっては、従来の軟鋼板を用いる場合でも発生し得る課題であり、これまでに数多くの知見、対策技術等が蓄積されている。一方、スプリングバックは、高強度鋼板を用いる場合に特有の課題であり、いまだ十分な検討がされておらず、広く普及している成形シミュレーションも、スプリングバック解析による寸法精度予測については、十分な実用信頼性を得ているとは言えない。

スプリングバックとは、成形品を成形後に金型から取り出す作業や、不要な部分をトリミングする作業等、部品の拘束を緩和する作業をすることで残留応力が駆動力となり、新たな釣り合いを満たすように部品に生じる弾性変形である。高強度鋼板ではこのスプリングバックが大きいため、最終製品として要求されている寸法精度の確保が困難になる。

スプリングバックは、現象に応じて、「角度変化」、「壁そり」、「ねじれ（ひねれ）」、「稜線そり（面そり）」、「パンチ底のスプリングバック」に分類される。いずれも、部品内での残留応力分布が曲げ又はねじりのモーメントとして働き、材料の弾性係数、板厚、部品形状等で決まる剛性に応じて部品が変形した結果として生じる。例えば、最も良く知られているスプリングバックの例は、曲げ角度変化、壁そり等である。これらのスプリングバックは、板厚方向の応力分布が駆動力となり、剛性は主に板厚で決定される。あるいは、長手方向に湾曲したハット断面のビームをドロー成形すると、壁そり及びねじれが生じるが、湾曲の曲率が小さいと部品剛性が高まり、壁そりが小さくなる。

このメカニズムに基づくと、寸法精度不良の対策方法の一つは、板厚、部品形状等を変更し、スプリングバックの弾性変形モードに応じた剛性を高め、スプリングバックに対する抵抗を増大させることである。

スプリングバックに対する抵抗の増大は、確実に寸法変化を低下させる。スプリングバックに対する抵抗は、その弾性変形モードに対する部品の剛性である。部品の剛性を高めるためには、部品の形状が重要な因子である。部品の形状は、性能、レイアウト等の要件で制約されているため、ビード、エンボス等の小さな対策が有効である。

壁そりに対するビードは、最も典型的な剛性補強対策例である。一方、複雑な部品での３次元的なスプリングバックに対して、最適な剛性補強位置を見出すのは容易でない。新たな試みとして、固有振動解析及び最適化を併用する方法が提案されている（非特許文献１を参照）。

即ち、固有振動数が剛性の平方根に比例し、且つ、密度（質量）の平方根に反比例することに着目し、スプリングバックの弾性変形モードを特定し、そのスプリングバックの弾性変形モードに対応する固有振動変形モードを選定し、その固有振動変形モードの固有振動数が上昇するように、要素の一部を同密度の高弾性材料に置き換え、その要素の最適な配置を、最適化ツール等を用いて求める。これにより、最適な剛性補強位置を容易に見出すことができる。

しかしながら、この方法では、固有振動変形モードを手動で選定しなければならないため、手間が掛かるという問題があった。また、固有振動変形モードを自動で選定するとしても、選定基準が明らかでないため、スプリングバックの弾性変形モードに対応する固有振動変形モードを選定できるとは限らず、高弾性材料の最適な配置が得られるとは限らなかった。そこで、すべての固有振動変形モードを選定し、それぞれについて高弾性材料の最適な配置を求めることも考えられるが、それではコストが掛かりすぎるという問題があった。このため、スプリングバックによって生じる部材の変形を低減するための解析を容易に行うことができなかった。

樋渡俊二、外１名、「高強度鋼板の成形課題への包括的取組み」、日本鉄鋼協会創形創質工学部会板成形フォーラム特別講演会、社団法人日本鉄鋼協会、２００５年１月１４日

本発明は、このような従来の課題を解決すべくなされたものであり、有限要素モデルデータに配置が容易な形状の矩形状ビードを配置し、当該矩形状ビードが配置された有限要素モデルデータに対して成型解析を行うことなく、スプリングバックによって生じる部材の変形を確実に低減するための解析を容易に行えるようにすることを目的とする。

本発明の解析装置は、部材の形状を示す有限要素モデルデータ、前記部材の材料物性を示す材料物性データ、及び前記部材に生じる応力分布を示す応力分布データを取得するプレス成形解析部と、前記有限要素モデルデータ、前記材料物性データ、及び前記応力分布データに基づいて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記スプリングバックによって生じる前記部材の変位分布を示す第１変位分布データを生成するスプリングバック解析部と、前記有限要素モデルデータ及び前記材料物性データに基づいて、各固有振動変形モードについて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記各固有振動変形モードでの前記部材の変位分布を示す第２変位分布データを生成する固有振動解析部と、前記第１変位分布データと前記各固有振動変形モードでの前記第２変位分布データのそれぞれとの間の一致度を求めるモード分解部と、前記一致度に基づいて、一つ又は複数の固有振動変形モードを選定するモード選定部と、前記一つ又は複数の固有振動変形モードのうち、少なくとも一部の固有振動数が上昇するように、前記有限要素モデルデータに非伏角状の矩形状ビードを配置するビード配置部とを含み、前記スプリングバック解析部は、前記矩形状ビードが配置された前記有限要素モデルデータに対して、取得された前記応力分布データに基づいてマッピングしてスプリングバック解析を行う。

本発明の解析装置において、前記部材は、フランジ部及び縦壁部を有するハット形状の鋼板であり、 前記ビード配置部は、前記部材の前記有限要素モデルデータについて、前記縦壁部のみに前記矩形状ビードを配置する。

本発明の解析方法は、部材の形状を示す有限要素モデルデータ、前記部材の材料物性を示す材料物性データ、及び前記部材に生じる応力分布を示す応力分布データを取得するプレス成形解析工程と、前記有限要素モデルデータ、前記材料物性データ、及び前記応力分布データに基づいて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記スプリングバックによって生じる前記部材の変位分布を示す第１変位分布データを生成する第１スプリングバック解析工程と、前記有限要素モデルデータ及び前記材料物性データに基づいて、各固有振動変形モードについて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記各固有振動変形モードでの前記部材の変位分布を示す第２変位分布データを生成する固有振動解析工程と、前記第１変位分布データと前記各固有振動変形モードでの前記第２変位分布データのそれぞれとの間の一致度を求めるモード分解工程と、前記一致度に基づいて、１つ又は複数の固有振動変形モードを選定するモード選定工程と、前記一つ又は複数の固有振動変形モードのうち、少なくとも一部の固有振動数が上昇するように、前記有限要素モデルデータに非伏角状の矩形状ビードを配置するビード配置工程と、前記ビード配置工程により前記矩形状ビードが配置された前記有限要素モデルデータに対して、前記プレス成形解析工程で取得された前記応力分布データに基づいてマッピングしてスプリングバック解析を行う第２スプリングバック解析工程とを少なくとも行う。

本発明の解析方法において、前記部材は、フランジ部及び縦壁部を有するハット形状の鋼板であり、前記ビード配置工程では、前記部材の前記有限要素モデルデータについて、前記縦壁部のみに前記矩形状ビードを配置する。

本発明のコンピュータプログラムは、部材の形状を示す有限要素モデルデータ、前記部材の材料物性を示す材料物性データ、及び前記部材に生じる応力分布を示す応力分布データを取得するプレス成形解析工程と、前記有限要素モデルデータ、前記材料物性データ、及び前記応力分布データに基づいて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記部材のスプリングバックによって生じる前記部材の変位分布を示す第１変位分布データを生成する第１スプリングバック解析工程と、前記有限要素モデルデータ及び前記材料物性データに基づいて、各固有振動変形モードについて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記各固有振動変形モードでの前記部材の変位分布を示す第２変位分布データを生成する固有振動解析工程と、前記第１変位分布データと前記各固有振動変形モードでの前記第２変位分布データのそれぞれとの間の一致度を求めるモード分解工程と、前記一致度に基づいて、一つ又は複数の固有振動変形モードを選定するモード選定工程と、前記一つ又は複数の固有振動変形モードのうち、少なくとも一部の固有振動数が上昇するように、前記有限要素モデルデータに非伏角状の矩形状ビードを配置するビード配置工程と、前記ビード配置工程により前記矩形状ビードが配置された前記有限要素モデルデータに対して、前記プレス成形解析工程で取得された前記応力分布データに基づいてマッピングしてスプリングバック解析を行う第２スプリングバック解析工程とをコンピュータに実行させる。

本発明のコンピュータプログラムにおいて、前記部材は、フランジ部及び縦壁部を有するハット形状の鋼板であり、前記ビード配置工程では、前記部材の前記有限要素モデルデータについて、前記縦壁部のみに前記矩形状ビードを配置する。

本発明によれば、有限要素モデルデータに配置が容易な形状の矩形状ビードを配置し、当該矩形状ビードが配置された有限要素モデルデータに対して成型解析を行うことなく、スプリングバックによって生じる部材の変形を確実に低減するための解析を容易に行うことが可能となる。

本実施形態による解析装置の動作の一例を説明するフローチャートである。 プレス成形後の応力分布が示された部材の一例を示す模式図（斜視図）である。 固有振動変形モードの次数と変位分布データ間の一致度との関係を示す特性図である。 有限要素モデルデータに矩形状ビードが配置された部材の一例を示す模式図である。 スプリングバック解析後の部材の変位分布の一例を示す模式図である。 本実施形態による解析装置の機能的な構成の一例を示すブロック図である。 パーソナルユーザ端末装置の内部構成を示す模式図である。

以下、図面を参照しつつ、本発明の様々な実施形態について説明する。但し、本発明の技術的範囲はそれらの実施形態に限定されず、特許請求の範囲に記載された発明とその均等物に及ぶ点に留意されたい。

図１は、解析装置１の動作の一例を説明するフローチャートである。なお、図１に示す処理フローは、図１１に示す解析装置１の記憶部１２に予め記憶されているプログラムに基づき、主に解析装置１の処理部１３により、解析装置１の各要素と協働して実行される。なお、解析装置１の構成については後述する。

最初に、処理部１３は、記憶部１２に記憶されている「ＣＡＤモデルデータ及び生成条件データ」に基づいて、部材の形状を、特定の形状及びサイズの複数の要素に分割することにより、有限要素モデルデータを生成する（ステップＳ１０１）。また、処理部１３は、生成した有限要素モデルデータを、記憶部１２に格納する。なお、有限要素モデルの生成処理には、市販のアプリケーションプログラムであるＨｙｐｅｒＭｅｓｈを利用する。しかしながら、他のアプリケーションプログラム（ＡＮＳＡ等）を利用することも可能である。

次に、処理部１３は、記憶部１２に記憶されている有限要素モデルデータ、材料物性データ、及びプレス成形条件データに基づいて、部材の各要素に発生する応力を有限要素法によって求めることにより、部材の応力分布データを生成する（ステップＳ１０２）。このように部材の応力分布データには、部材の各要素に発生する応力を示す情報が含まれる。また、処理部１３は、生成した応力分布データを、記憶部１２に格納する。なお、プレス成形の解析処理には、市販のアプリケーションプログラムであるＨｙｐｅｒＦｏｒｍを利用する。しかしながら、他のアプリケーションプログラム（ＬＳ−ＤＹＮＡ、ＰＡＭ−ＳＴＡＭＰ等）を利用することも可能である。

図２は、プレス成形後の応力分布が示された部材の一例を示す模式図（斜視図）である。図２では、処理部１３により求められた応力の大きさ（下死点応力分布）が色の濃淡により示されている。当該部材は、フランジ部及び縦壁部を有するハット形状の９８０ＭＰａ級の高張力鋼板であり、板厚は１．６ｍｍである。

次に、処理部１３は、記憶部１２に記憶されている有限要素モデルデータ、応力分布データ、材料物性データ、及び境界条件データに基づいて、スプリングバック後の部材の各節点の変位を有限要素法によって求めることにより、部材の第１変位分布データを生成する（ステップＳ１０３）。このように第１変位分布データには、部材の各節点の変位を示す情報が含まれる。また、処理部１３は、生成した第１変位分布データを、記憶部１２に格納する。なお、スプリングバックの解析処理には、市販のアプリケーションプログラムであるＬＳ−ＤＹＮＡを利用する。しかしながら、他のアプリケーションプログラム（ＰＡＭ−ＳＴＡＭＰ等）を利用することも可能である。

次に、処理部１３は、記憶部１２に記憶されている有限要素モデルデータ、材料物性データ、及び境界条件データに基づいて、各固有振動変形モードについて、部材の各節点の変位を有限要素法によって求めることにより、部材の第２変位分布データを生成する（ステップＳ１０４）。このように第２変位分布データには、部材の各節点の変位を示す情報が含まれる。また、処理部１３は、生成した第２変位分布データのそれぞれを、記憶部１２に格納する。なお、固有振動の解析処理には、市販のアプリケーションプログラムであるＯｐｔｉＳｔｒｕｃｔを利用する。しかしながら、他のアプリケーションプログラム（Ｎａｓｔｒａｎ等）を利用することも可能である。

次に、処理部１３は、記憶部１２に記憶されている「スプリングバックにより生じた変形に基づく第１変位分布データと、各固有振動変形モードでの第２変位分布データのそれぞれ」と、の各組み合わせについて、部材の各節点の変位の一致度を求めることにより、変位分布データ間の一致度を求める（ステップＳ１０５）。処理部１３は、部材の各節点について、各固有振動変形モードでの変位と、スプリングバックにより生じた変形に基づく変位とに、どの程度の関連性があるのかを求めることにより、変位分布データ間の一致度を求める。具体的には、処理部１３は、次の（１）式に示す連立一次方程式を解くことにより、変位分布データ間の一致度を求める。ここで、ｕ_ijはｉ次の固有振動変形モードでの節点ｊの変位を示し、ｕ_SBjはスプリングバック変形での節点ｊの変位を示し、ａ_iはｉ次の固有振動変形モードとスプリングバック変形の一致度を示す。なお、ａ_iの大きさが大きいほど変位分布データ間の一致度が高いことを示す。

あるいは、処理部１３は、部材の各節点の変位後の位置の差分を求めることにより、変位分布データ間の一致度を求めることも可能である。具体的には、処理部１３は、次の（２）式又は（３）式に示す数式を計算することにより、変位分布データ間の一致度を求めることも可能である。ここで、ｕ_ijはｉ次の固有振動変形モードでの節点ｊの変位を示し、ｕ_SBjはスプリングバック変形での節点ｊの変位を示し、及びｂ_i、ｃ_iはｉ次の固有振動変形モードとスプリングバック変形の一致度を示す。なお、ｂ_i、ｃ_iの大きさが大きいほど変位分布データ間の一致度が高いことを示す。

処理部１３は、求めた変位分布データ間の一致度のそれぞれを、記憶部１２に格納する。尚、本明細書では、ｉ次の固有振動変形モードとスプリングバック変形の一致度ａ_i、ｂ_i、ｃ_iのうち、ａ_iを用いた場合を例に挙げて説明する。

次に、処理部１３は、記憶部１２に記憶されている変位分布データ間の一致度のそれぞれに基づいて、一つ又は複数（少数）の固有振動変形モードを選定する（ステップＳ１０６）。処理部１３は、変位分布データ間の一致度の大きいものから順に、所定数（例えば、３つ）の固有振動変形モードを選定する。あるいは、処理部１３は、変位分布データ間の一致度が所定の閾値よりも大きい、又は所定の閾値以上である一つ又は複数の固有振動変形モードを選定することも可能である。このように処理部１３は、変位分布データ間の一致度に基づいて固有振動変形モードを選定する。このため、少数の固有振動変形モードを選定するのにも関わらず、スプリングバックによって生じる部材の変形に対応する固有振動変形モードが選定される可能性が高くなる。そして、処理部１３は、選定した固有振動変形モードの次数を、記憶部１２に格納する。

次に、処理部１３は、記憶部１２に記憶されている次数の固有振動変形モードから、固有振動数を同時に上昇させるものの組み合わせを生成する。処理部１３は、変位分布データ間の一致度の大きいものから順に、一つずつ固有振動変形モードを累加していくことにより、ｍ個の固有振動変形モード（例えば、８次、２次、及び１２次）からｎ個の組み合わせ（例えば、｛８次｝、｛８次、２次｝、及び｛８次、２次、１２次｝）を生成する（ｍ、ｎは、１以上の整数であり、同じ数であっても異なる数であっても良い。）。あるいは、処理部１３は、変位分布データ間の一致度の大きいものから順に、二つずつ固有振動変形モードを累加していくことも可能である。

例えば図３に示すように、固有振動変形モードの次数と変位分布データ間の一致度との関係が得られる。ここで、横軸は、固有振動変形モードの次数を示す。縦軸は、変位分布データ間（スプリングバック変形に基づく第１変位分布データと、各固有振動変形モードでの第２変位分布データとの間）の一致度を示す。この関係から、８次モードが他のモードに比べて変位分布データ間の一致度が大きいことが分かる。

次に、処理部１３は、処理部１３は、生成した固有振動変形モードの組み合わせのそれぞれについて、当該組み合わせに含まれる固有振動変形モードの固有振動数の少なくとも一つが上昇するように、記憶部１２に記憶されている材料物性データ、境界条件データ、及び配置条件データに基づいて、記憶部１２に記憶されている有限要素モデルデータにビードを配置する。例えば固有振動変形モードの次数と変位分布データ間の一致度との関係が図３のようである場合、変位分布データ間の一致度が最も大きい８次モードにより得られる部材の形状が最適化形状であることになる。

本実施形態では、当該有限要素モデルデータに、ビードとして、非伏角状の（伏角形状のない）矩形状ビードを配置する（ステップＳ１０７）。ビードの形状は、一般的に複雑となる。従って、計算通りにビードを部材に配置することは容易ではない。そこで本実施形態では、後述するスプリングバック解析を行う前に、部材の形状・取り付け位置及び他の部材との関係等とを考慮して、有限要素モデルデータに矩形状ビードを配置する。具体的には、ハット形状の鋼板である部材の有限要素モデルデータについて、縦壁部からフランジ部にかけて矩形状ビードを配置することはせず、縦壁部のみに矩形状ビードを配置するようにする。

図４は、有限要素モデルデータに矩形状ビードが配置された部材の一例を示す模式図（斜視図）であり、（ａ）が計算通りにビードを部材に配置した比較例を、（ｂ）が本実施形態をそれぞれ示している。比較例では、部材に配置されたビード１０１が極めて複雑であり、実際にビード１０１のような複雑な形状のビードを部材に付すことは困難である。これに対して本実施形態では、部材には、比較的単純な形状であるが最適化形状である矩形状ビード１０２が付されている。実際の部材に付すことができる形状の矩形状ビードを計算段階で有限要素モデルデータに付すことにより、後述するように、容易且つ確実にスプリングバックを低減することができる。

ここで、本実施形態のように部材の有限要素モデルデータに矩形状ビードを付与するのではなく、例えば図１のステップＳ１０１〜Ｓ１１１により得られた情報に基づいて、部材に対する実際の矩形状ビードを付与することも考えられる。しかしながらこの場合、計算（処理装置）によって最適化された矩形状ビードの形状に基づいて、人手により矩形状ビードを修正して付与しなければならず、十分なスプリングバックの低減を実現することができなくなるという問題がある。更にこの場合、ビードの無い形状の部材を成形した後に、追加工程で矩形状ビードの成形を行うため、工程増を招くという問題がある。

これに対して本実施形態では、部材の有限要素モデルデータに矩形状ビードを付与する構成を採るため、人手によることなく計算（処理装置）によって最適化された矩形状ビードの形状を求めることができる。これにより、スプリングバックを低減するという技術的効果を最大化して、部材のどこの位置に、どの程度の幅・高さで矩形状ビードを配置するかということを正確に求めることができる。更に本実施形態では、部材の有限要素モデルデータに矩形状ビードを付与する構成により、追加工程を行うことなく、１つの工程で矩形状ビードを備えた部材を成形することができる。

また本実施形態では、部材の有限要素モデルデータについて、縦壁部のみに矩形状ビードを配置するようにする。部材の有限要素モデルデータにおいて、フランジ部から縦壁部にかけて繋がったビードを配置すると、成形初期にしわ押えとダイスによりブランクを挟んだ際に、部材にしわが発生するという問題がある。本実施形態では、縦壁部のみに矩形状ビードを配置することにより、部材のしわの発生が抑止される。

ここで、有限要素モデルデータに矩形状ビードを配置することは、有限要素モデルデータに含まれる要素に属する部材の表面に伏角形状のない矩形状の凹凸をつけることを意味する。このように本発明におけるビードは、部材の表面に矩形状の凹凸を付けることを意味し、エンボス等も含む概念である。また、上記の組み合わせに含まれる固有振動変形モードが複数である場合、当該複数の固有振動変形モードの固有振動数の全てを上昇させるのが好ましい。しかしながら、当該組み合わせに含まれる固有振動変形モードによっては、このようにすることができない場合がある。従って、このような場合には、当該複数の固有振動変形モードの固有振動数の少なくとも１つを上昇させる。

そして、処理部１３は、矩形状ビードが配置された後の有限要素モデルデータのそれぞれを、記憶部１２に格納する。なお、矩形状ビードの配置処理には、市販のアプリケーションプログラムであるＯｐｔｉＳｔｒｕｃｔを利用する。また、矩形状ビードの配置処理に際し、部材面積に対するビード配置面積の割合を指定することができる。ここで、部材面積とは、部材の面のうち、矩形状ビードを配置することが可能な面の面積である。ビード配置面積とは、部材の面のうち、矩形状ビードが配置される部分の面積である。この割合を大きくしすぎると、部材の大部分の領域にビードを配置することになる。しかしながら、部材の形状、部材の取り付け箇所、及び他の部材との関係等によって、矩形状ビードを配置する領域に制限が課せられる場合がある。また、スプリングバックによって生じる部材の変形を抑制するのに大きく寄与する領域のみを、矩形状ビードを配置する領域とし、部材の設計を行い易くしたい場合がある。このような観点から、ユーザは、操作部１４を操作することによって、部材面積に対するビード配置面積の割合を適宜指定する。

次に、処理部１３は、記憶部１２に記憶されている「矩形状ビードが配置された後の各有限要素モデルデータ」について、当該有限要素モデルデータに、記憶部１２に記憶されているステップＳ１０２で取得した応力分布データ（残留応力データ）をマッピングし、スプリングバック解析を行う（ステップＳ１０８）。本実施形態では、「矩形状ビードが配置された後の各有限要素モデルデータ」について成形解析を行うことなく、当該スプリングバック解析を行う。処理部１３は、当該残留応力データをマッピングした有限要素モデルデータと、記憶部１２に記憶されている「材料物性データ及び境界条件データ」とに基づいて、部材の各節点の変位を有限要素法によって求めることにより、部材の第１変位分布データを生成する。また、処理部１３は、このようにして生成した「矩形状ビードが配置された後の第１変位分布データ」を、記憶部１２に格納する。なお、スプリングバックの解析処理には、市販のアプリケーションプログラムであるＬＳ−ＤＹＮＡを利用する。しかしながら、他のアプリケーションプログラム（ＰＡＭ−ＳＴＡＭＰ等）を利用することも可能である。

次に、処理部１３は、矩形状ビードが配置された後の第１変位分布データに基づいて定まる部材の形状と、当該部材の目標形状との差が所定の範囲内であるか否かを判定する（ステップＳ１０９）。この判定は、例えば、各要素の節点の変位後の複数の位置と、当該位置に対応する目標形状の位置との差の全てが所定の範囲内であるか否かを判定することにより実現できる。

この判定の結果、矩形状ビードが配置された後の第１変位分布データに基づいて定まる部材の形状と、当該部材の目標形状との差が所定の範囲内でない場合には、ステップＳ１０４に戻る。そして、矩形状ビードが配置された後の第１変位分布データに基づいて定まる部材の形状と、当該部材の目標形状との差が所定の範囲内になるまで、記憶部１２に記憶されている「矩形状ビードが配置された後の（最新の）各有限要素モデルデータ」を用いて、ステップＳ１０４〜Ｓ１０９の処理を繰り返し行う。そして、矩形状ビードが配置された後の第１変位分布データに基づいて定まる部材の形状と、当該部材の目標形状との差が所定の範囲内になると、処理部１３は、スプリングバック後の部材の形状と目標形状との差が所定の範囲内になる矩形状ビードの配置（ステップＳ１０７で得られた最新の矩形状ビードの配置）を特定する情報等を表示部１５に表示させる。そして、図１のフローチャートによる処理を終了する。矩形状ビードの配置を特定する情報としては、例えば、矩形状ビードの位置、縦横の大きさ、高さ（深さ）等の情報を含めることができる。また、スプリングバック後の部材の形状と目標形状との差が所定の範囲内になるときに選定された固有振動変形モード（ステップＳ１０６で得られた最新の固有振動変形モード）を示す情報を表示部１５に表示させても良い。

図５は、スプリングバック解析後の部材の変位分布の一例を示す模式図（斜視図）である。８次の固有振動モードの組み合わせについて処理部１３により求められた矩形状ビードの配置が行われた後の部材の変位分布を示している。

本実施形態では、複数の固有振動モードの組わせたビード配置の候補で、それぞれ固有振動数を求めて、固有振動数の最大値が最大となる固有振動モードの組み合わせのビード配置を選定すれば、部材の変位を大きく低減することができる。

解析装置１の処理部１３は、有限要素モデルデータに対して、プレス成形解析、スプリングバック解析、固有振動解析等を実行する。また、処理部１３は、スプリングバック変形との一致度に基づいて少数の固有振動変形モードを選定し、それらの各組み合わせについて、当該組み合わせに含まれる固有振動変形モードの固有振動数が上昇するように、有限要素モデルデータに矩形状ビードを配置する。そして、処理部１３は、矩形状ビードが配置された後の有限要素モデルデータに対して、プレス成形解析を行うことなく、スプリングバック解析を再度実行する。そして、スプリングバック解析後の部材の形状と目標形状との差が所定の範囲内になるまで、矩形状ビードの配置と、プレス成形解析と、スプリングバック解析とを繰り返し行う。これにより、スプリングバック変形に対応する固有振動変形モードを容易に選定することができ、容易且つ確実にスプリングバックを低減することが可能となる。

以下、本実施形態における解析装置１のハードウェア構成について説明する。図６は、解析装置１の機能的な構成の一例を示すブロック図である。

解析装置１は、すでにインストールされているプログラムを実行し、記憶部１２に記憶されているデータ及び／又は他の装置に記憶されているデータを参照し、各種の処理を実行する。また、解析装置１は、ユーザから操作部１４を介して入力された指示に従って各種の処理を実行し、その結果を、表示部１５を介してユーザに提示する。そのために、解析装置１は、通信部１１と、記憶部１２と、処理部１３と、操作部１４と、表示部１５とを有する。

通信部１１は、解析装置１を不図示のネットワークに接続するための通信インターフェース回路を有する。通信部１１は、不図示の他の装置からネットワークを介して受信したデータを、処理部１３に渡す。また、通信部１１は、処理部１３から受け取ったデータを、ネットワークを介して他の装置に送信する。

記憶部１２は、例えば、半導体メモリ、磁気ディスク装置、又は光ディスク装置のうちの少なくとも何れか一つを有する。記憶部１２は、処理部１３での処理に用いられるアプリケーションプログラム、データ等を記憶する。記憶部１２は、例えば、アプリケーションプログラムとして、有限要素モデル生成プログラム、プレス成形解析プログラム、スプリングバック解析プログラム、固有振動解析プログラム、モード分解プログラム、モード選定プログラム、ビード配置プログラム等を記憶する。

また、記憶部１２は、予め設定される初期設定データとして、以下のデータ等を記憶する。
第１の初期設定データとして、記憶部１２は、部材の位置・大きさ・形状を示すＣＡＤ（Computer Aided Design）モデルデータを記憶する。第２の初期設定データとして、記憶部１２は、部材の材料物性（寸法、板厚、材料、ヤング率、ポワソン比、質量密度、応力とひずみとの関係等）を示す材料物性データを記憶する。第３の初期設定データとして、記憶部１２は、有限要素モデルの生成条件（要素の形状・大きさ等）を示す生成条件データを記憶する。第４の初期設定データとして、記憶部１２は、プレス成形条件（摩擦係数、部材フランジ押さえ力等）を示すプレス成形条件データを記憶する。第５の初期設定データとして、記憶部１２は、境界条件（部材上の固定点等）を示す境界条件データを記憶する。第６の初期設定データとして、記憶部１２は、矩形状ビードの配置条件（縦横の大きさ、高さ、配置面積等）を示す配置条件データ等を記憶する。

また、記憶部１２は、処理部１３により生成される中間データとして、以下のデータ等を記憶する。
即ち、第１の中間データとして、記憶部１２は、ＣＡＤモデルデータに対応する有限要素モデルデータ（各要素の位置・形状・大きさ等）を記憶する。第２の中間データとして、記憶部１２は、部材の応力分布を示す応力分布データを記憶する。第３の中間データとして、記憶部１２は、スプリングバックによって生じる部材の変位分布を示す第１変位分布データを記憶する。第４の中間データとして、記憶部１２は、各固有振動変形モードでの部材の変位分布を示す第２変位分布データを記憶する。第５の中間データとして、記憶部１２は、変位分布データ間の一致度を記憶する。第６の中間データとして、記憶部１２は、固有振動変形モードの選定次数を記憶する。第７の中間データとして、記憶部１２は、矩形状ビードが配置された後の有限要素モデルデータを記憶する。第８の中間データとして、記憶部１２は、矩形状ビードが配置された後の変位分布を示す変位分布データを記憶する。更に、記憶部１２は、所定の処理に係る一時的なデータを、一時的に記憶しても良い。

処理部１３は、１個または複数個のプロセッサ及びその周辺回路を有する。処理部１３は、解析装置１の全体的な動作を統括的に制御する処理部、例えばＣＰＵ（Central Processing Unit）である。即ち、処理部１３は、解析装置１の各種の処理が操作部１４の操作、記憶部１２に記憶されているプログラム等に基づいて適切な手順で実行されるように、通信部１１、表示部１５等の動作を制御する。処理部１３は、記憶部１２に記憶されているプログラム（オペレーティングシステムプログラム、アプリケーションプログラム等）に基づいて処理を実行する。また、処理部１３は、複数のプログラム（アプリケーションプログラム等）を並列に実行することができる。

処理部１３は、図１におけるステップＳ１０１の処理を実行する有限要素モデル生成部１３１と、ステップＳ１０２の処理を実行するプレス成形解析部１３２と、ステップＳ１０３及びＳ１０８の処理を実行するスプリングバック解析部１３３と、ステップＳ１０４の処理を実行する固有振動解析部１３４と、ステップＳ１０５の処理を実行するモード分解部１３５と、ステップＳ１０６の処理を実行するモード選定部１３６と、ステップＳ１０７の処理を実行するビード配置部１３７と、ステップＳ１０９の処理を実行する形状判定部１３８と、ステップＳ１１０の処理を実行するデータ描画部１３９とを有する。処理部１３が有するこれらの各部は、処理部１３が有するプロセッサ上で実行されるプログラムによって実装される機能モジュールである。あるいは、処理部１３が有するこれらの各部は、ファームウェアとして解析装置１に実装されても良い。

データ描画部１３９は、データの描画処理を実行する。即ち、有限要素モデル生成部１３１、プレス成形解析部１３２、スプリングバック解析部１３３、固有振動解析部１３４、モード分解部１３５、モード選定部１３６、及びビード配置部１３７から与えられたデータを解析し、そのデータを所定の形式（例えば、コンター図）でレンダリングし、その描画データを生成する。そして、データ描画部１３９は、生成した描画データを表示部１５等に出力する。このようにした場合、表示部１５が、出力部として機能する。しかしながら、必ずしもこのようにする必要はない。例えば、通信部１１が、有限要素モデル生成部１３１、プレス成形解析部１３２、スプリングバック解析部１３３、固有振動解析部１３４、モード分解部１３５、モード選定部１３６、及びビード配置部１３７から与えられたデータを外部装置に送信する場合には、通信部１１が出力部として機能する。

操作部１４は、解析装置１の操作が可能であればどのようなデバイスでもよく、例えば、キーボード、タッチパネル等である。ユーザは、このデバイスを用いて、選択等の指示を入力することが可能となる。操作部１４は、ユーザにより操作されると、その操作に対応する信号を発生する。そして、発生した信号は、ユーザの指示として、処理部１３に入力される。

表示部１５も、映像、画像等の表示が可能であればどのようなデバイスでもよく、例えば、液晶ディスプレイ、有機ＥＬ（Electro−Luminescence）ディスプレイ等である。表示部１５は、処理部１３から供給される描画データに応じた映像、画像等を表示する。

以上説明してきたように、本実施形態によれば、有限要素モデルデータに配置が容易な矩形状ビードを配置し、当該矩形状ビードが配置された有限要素モデルデータに対して成型解析を行うことなく、スプリングバックによって生じる部材の変形を確実に低減するための解析を容易に行うことが可能となる。

なお、本発明は、本実施形態に限定されるものではない。例えば、本実施形態では、解析装置１は、図６に示す各部を有するとしたが、その一部については、不図示のサーバ装置が有するとしてもよい。サーバ装置は、例えば、解析装置１の記憶部１２に相当する記憶部を有し、記憶部に記憶されているプログラム、データ等を解析装置１に提供し、解析装置１に解析処理を実行させるようにしてもよい。このようにした場合、解析装置１の処理部１３は、サーバ装置から通信部１１を介してプログラム、データ等を取得する。一方、解析装置１の記憶部１２にプログラム、データ等を記憶する場合には、処理部１３が記憶部１２からプログラム、データ等を取得することになる。

また、サーバ装置は、解析装置１の記憶部１２及び処理部１３に相当する記憶部及び処理部を有し、記憶部に記憶されているプログラム、データ等を用いて解析処理を実行し、その結果のみを解析装置１に提供するようにしても良い。

また、解析装置１の処理部１３が有する各機能をコンピュータに実現させるためのプログラムは、磁気記録媒体または光記録媒体といったコンピュータによって読み取り可能な記録媒体に記録された形で提供されても良い。

即ち、上記した本実施形態は、コンピュータがプログラムを実行することによって実現することができる。本実施形態による解析装置の各構成要素（図１１の処理部１３の各部１３１〜１３９等）の機能は、コンピュータのＲＡＭやＲＯＭ等に記憶されたプログラムが動作することによって実現できる。同様に、本実施形態による解析方法の各ステップ（図１のステップＳ１０１〜Ｓ１１０等）は、コンピュータのＲＡＭやＲＯＭ等に記憶されたプログラムが動作することによって実現できる。

具体的に、上記のプログラムは、例えばＣＤ−ＲＯＭのような記録媒体に記録し、或いは各種伝送媒体を介し、コンピュータに提供される。上記のプログラムを記録する記録媒体としては、ＣＤ−ＲＯＭ以外に、フレキシブルディスク、ハードディスク、磁気テープ、光磁気ディスク、不揮発性メモリカード等を用いることができる。他方、上記のプログラムの伝送媒体としては、プログラム情報を搬送波として伝搬させて供給するためのコンピュータネットワークシステムにおける通信媒体を用いることができる。ここで、コンピュータネットワークとは、ＬＡＮ、インターネットの等のＷＡＮ、無線通信ネットワーク等であり、通信媒体とは、光ファイバ等の有線回線や無線回線等である。

また、本発明に含まれるプログラムとしては、供給されたプログラムをコンピュータが実行することにより本実施形態の機能が実現されるようなもののみではない。例えば、そのプログラムがコンピュータにおいて稼働しているＯＳ（オペレーティングシステム）或いは他のアプリケーションソフト等と共同して本実施形態の機能が実現される場合にも、かかるプログラムは本発明に含まれる。また、供給されたプログラムの処理の全て或いは一部がコンピュータの機能拡張ボードや機能拡張ユニットにより行われて本実施形態の機能が実現される場合にも、かかるプログラムは本発明に含まれる。

例えば、図７は、パーソナルユーザ端末装置の内部構成を示す模式図である。この図７において、１２００はＣＰＵ１２０１を備えたパーソナルコンピュータ（ＰＣ）である。ＰＣ１２００は、ＲＯＭ１２０２またはハードディスク（ＨＤ）１２１１に記憶された、又はフレキシブルディスクドライブ（ＦＤ）１２１２より供給されるデバイス制御ソフトウェアを実行する。このＰＣ１２００は、システムバス１２０４に接続される各デバイスを総括的に制御する。

ＰＣ１２００のＣＰＵ１２０１、ＲＯＭ１２０２またはハードディスク（ＨＤ）１２１１に記憶されたプログラムにより、本実施形態の図１におけるステップＳ１０１〜Ｓ１１０の手順等が実現される。

１２０３はＲＡＭであり、ＣＰＵ１２０１の主メモリ、ワークエリア等として機能する。１２０５はキーボードコントローラ（ＫＢＣ）であり、キーボード（ＫＢ）１２０９や不図示のデバイス等からの指示入力を制御する。

１２０６はＣＲＴコントローラ（ＣＲＴＣ）であり、ＣＲＴディスプレイ（ＣＲＴ）１２１０の表示を制御する。１２０７はディスクコントローラ（ＤＫＣ）である。ＤＫＣ１２０７は、ブートプログラム、複数のアプリケーション、編集ファイル、ユーザファイル、ネットワーク管理プログラム等を記憶するハードディスク（ＨＤ）１２１１、及びフレキシブルディスク（ＦＤ）１２１２とのアクセスを制御する。ここで、ブートプログラムとは、パソコンのハードやソフトの実行（動作）を開始する起動プログラムである。

１２０８はネットワーク・インターフェースカード（ＮＩＣ）であり、ＬＡＮ１２２０を介して、ネットワークプリンタ、他のネットワーク機器、或いは他のＰＣと双方向のデータのやり取りを行う。
なお、パーソナルユーザ端末装置を用いる代わりに、本実施形態の解析装置に特化された所定の計算機等を用いても良い。

本発明は、例えば、自動車の車体等に適用される部材の成形に利用することができる。

１１ 通信部
１２ 記憶部
１３ 処理部
１４ 操作部
１５ 表示部

Claims (6)

1. 部材の形状を示す有限要素モデルデータ、前記部材の材料物性を示す材料物性データ、及び前記部材に生じる応力分布を示す応力分布データを取得するプレス成形解析部と、
   前記有限要素モデルデータ、前記材料物性データ、及び前記応力分布データに基づいて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記スプリングバックによって生じる前記部材の変位分布を示す第１変位分布データを生成するスプリングバック解析部と、
   前記有限要素モデルデータ及び前記材料物性データに基づいて、各固有振動変形モードについて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記各固有振動変形モードでの前記部材の変位分布を示す第２変位分布データを生成する固有振動解析部と、
   前記第１変位分布データと前記各固有振動変形モードでの前記第２変位分布データのそれぞれとの間の一致度を求めるモード分解部と、
   前記一致度に基づいて、一つ又は複数の固有振動変形モードを選定するモード選定部と、
   前記一つ又は複数の固有振動変形モードのうち、少なくとも一部の固有振動数が上昇するように、前記有限要素モデルデータに非伏角状の矩形状ビードを配置するビード配置部と、
   を含み、
   前記スプリングバック解析部は、前記矩形状ビードが配置された前記有限要素モデルデータに対して、取得された前記応力分布データに基づいてマッピングしてスプリングバック解析を行うことを特徴とする解析装置。
2. 前記部材は、フランジ部及び縦壁部を有するハット形状の鋼板であり、
   前記ビード配置部は、前記部材の前記有限要素モデルデータについて、前記縦壁部のみに前記矩形状ビードを配置することを特徴とする請求項１に記載の解析装置。
3. 部材の形状を示す有限要素モデルデータ、前記部材の材料物性を示す材料物性データ、及び前記部材に生じる応力分布を示す応力分布データを取得するプレス成形解析工程と、
   前記有限要素モデルデータ、前記材料物性データ、及び前記応力分布データに基づいて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記スプリングバックによって生じる前記部材の変位分布を示す第１変位分布データを生成する第１スプリングバック解析工程と、
   前記有限要素モデルデータ及び前記材料物性データに基づいて、各固有振動変形モードについて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記各固有振動変形モードでの前記部材の変位分布を示す第２変位分布データを生成する固有振動解析工程と、
   前記第１変位分布データと前記各固有振動変形モードでの前記第２変位分布データのそれぞれとの間の一致度を求めるモード分解工程と、
   前記一致度に基づいて、１つ又は複数の固有振動変形モードを選定するモード選定工程と、
   前記一つ又は複数の固有振動変形モードのうち、少なくとも一部の固有振動数が上昇するように、前記有限要素モデルデータに非伏角状の矩形状ビードを配置するビード配置工程と、
   前記ビード配置工程により前記矩形状ビードが配置された前記有限要素モデルデータに対して、前記プレス成形解析工程で取得された前記応力分布データに基づいてマッピングしてスプリングバック解析を行う第２スプリングバック解析工程と、
   を少なくとも行うことを特徴とする解析方法。
4. 前記部材は、フランジ部及び縦壁部を有するハット形状の鋼板であり、
   前記ビード配置工程では、前記部材の前記有限要素モデルデータについて、前記縦壁部のみに前記矩形状ビードを配置することを特徴とする請求項３に記載の解析方法。
5. 部材の形状を示す有限要素モデルデータ、前記部材の材料物性を示す材料物性データ、及び前記部材に生じる応力分布を示す応力分布データを取得するプレス成形解析工程と、
   前記有限要素モデルデータ、前記材料物性データ、及び前記応力分布データに基づいて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記部材のスプリングバックによって生じる前記部材の変位分布を示す第１変位分布データを生成する第１スプリングバック解析工程と、
   前記有限要素モデルデータ及び前記材料物性データに基づいて、各固有振動変形モードについて、前記有限要素モデルデータの節点の変位を有限要素法によって求めることにより、前記各固有振動変形モードでの前記部材の変位分布を示す第２変位分布データを生成する固有振動解析工程と、
   前記第１変位分布データと前記各固有振動変形モードでの前記第２変位分布データのそれぞれとの間の一致度を求めるモード分解工程と、
   前記一致度に基づいて、一つ又は複数の固有振動変形モードを選定するモード選定工程と、
   前記一つ又は複数の固有振動変形モードのうち、少なくとも一部の固有振動数が上昇するように、前記有限要素モデルデータに非伏角状の矩形状ビードを配置するビード配置工程と、
   前記ビード配置工程により前記矩形状ビードが配置された前記有限要素モデルデータに対して、前記プレス成形解析工程で取得された前記応力分布データに基づいてマッピングしてスプリングバック解析を行う第２スプリングバック解析工程と、
   をコンピュータに実行させることを特徴とするコンピュータプログラム。
6. 前記部材は、フランジ部及び縦壁部を有するハット形状の鋼板であり、
   前記ビード配置工程では、前記部材の前記有限要素モデルデータについて、前記縦壁部のみに前記矩形状ビードを配置することを特徴とする請求項５に記載のコンピュータプログラム。