(57)【要約】    【課題】  有限要素法を用いた板材成形のシミュレーションシステムにおいて、シミュレーションの簡素化及び高速化を図り、スプリングバックを評価できるシステムを提供する。    【解決手段】  複数工程による成形結果をシミュレーションするシステムであって、最終設計形状及び中間設計形状を入力する入力装置と、各種処理を行う処理装置と、処理結果を出力する出力装置とを備え、処理装置は、設計形状を有限要素分割して解析モデルを作成する解析モデル作成手段と、各要素を所定平面に平行になるように回転させる要素回転手段と、各要素を所定平面に投影させて、解析モデルで同一節点を共有する節点同士が一致するように各要素を変形させて平面展開モデルを作成する展開モデル作成手段と、各要素の変形に応じた応力等を算出する演算手段と、各設計形状に対する平面展開モデルの各要素の応力等に基いて各工程による成形結果を評価する評価手段とを有する。【選択図】    図３１
2006185228.tif 000002
【特許請求の範囲】
【請求項１】
板材をブランクとする複数工程による成形結果を有限要素法を用いてシミュレーションするシステムであって、最終設計形状及び最終設計形状に至る１又は複数の中間設計形状を入力する入力装置と、入力された複数の設計形状について各種処理を行う処理装置と、該処理装置による処理結果を出力する出力装置とが備えられていると共に、上記処理装置は、上記入力装置で入力された設計形状を有限要素分割して解析モデルを作成する解析モデル作成手段と、該作成手段で作成した解析モデルの各要素を所定平面に平行になるようにその中心回りに回転させる要素回転手段と、該回転手段で回転された各要素を上記所定平面上に投影させた状態で、上記解析モデルで同一節点を共有していた節点同士が一致するように応力の釣合を保ちながら各要素を変形させて解析モデルの平面展開モデルを作成する展開モデル作成手段と、該作成手段で作成した平面展開モデルの各要素についてそれぞれの変形に応じた応力、歪、又は厚みの少なくとも１つを算出する演算手段と、各設計形状について同様に作成した平面展開モデルの各要素の応力、歪、又は厚みの少なくとも１つに基いて各工程による成形結果を評価する評価手段とを有することを特徴とする多段階成形シミュレーションシステム。
【請求項２】
要素回転手段は、要素中心の法線ベクトルと要素中心を通って所定平面に直交する所定平面直交ベクトルとを作成し、上記要素中心の法線ベクトルが上記所定平面直交ベクトルに一致するように各要素を回転させることを特徴とする請求項１に記載の多段階成形シミュレーションシステム。
【請求項３】
展開モデル作成手段は、回転手段で回転された各要素を所定平面上に投影させた状態で、要素回転前の解析モデルで同一節点を共有していた節点同士に各節点間距離に応じた引っ張り荷重を設定し、この引っ張り荷重により解析モデルで同一節点を共有していた節点同士が引っ張り荷重の釣合を保ちながら一致するように各要素を変形させることを特徴とする請求項１または請求項２に記載の多段階成形シミュレーションシステム。
【請求項４】
所定の工程で展開モデル作成手段により得られた平面展開モデルとその工程前の平面展開モデルの対応する位置にそれぞれ基準点を設定し、これらの基準点を用いて両平面展開モデルの位置合せを行う位置合せ手段が備えられ、演算手段は、該位置合せ手段によって位置合せされた当該工程前後の平面展開モデルの互いに対応する要素における歪の差を演算し、その差に基いて、当該工程後の平面展開モデルにおける各要素の応力を計算することを特徴とする請求項１から請求項３のいずれかに記載の多段階成形シミュレーションシステム。
【請求項５】
中間設計形状を得る少なくとも１つの工程にその中間設計形状として複数の設計形状が設定されている場合に、この複数の中間設計形状をそれぞれ経由して最終設計形状に至る複数の成形経路を設定する成形経路設定手段が備えられ、演算手段は、各経路ごとに経路を構成する各設計形状について平面展開モデルの各要素の応力、歪、又は厚みの少なくとも１つを算出すると共に、それらの算出結果に基いて最適な成形経路を判定する判定手段が備えられていることを特徴とする請求項１から請求項４のいずれかに記載の多段階成形シミュレーションシステム。
【請求項６】
平面展開モデルの各要素ごとに厚み方向に複数のレイヤー要素を設定するレイヤー要素設定手段と、所定の工程で展開モデル作成手段により得られた平面展開モデルとその工程前の平面展開モデルとの対応する位置にそれぞれ基準点を設定し、これらの基準点を用いて両平面展開モデルの位置合せを行う位置合せ手段と、当該工程前後の設計形状の曲率から平面展開モデルにおける各レイヤー要素の曲げ歪を算出する曲げ歪算出手段と、レイヤー要素ごとに、演算手段を用いて算出した歪と上記曲げ歪算出手段により算出した曲げ歪とを加算して合計歪を算出する合計歪算出手段と、当該工程前後の平面展開モデルの互いに対応するレイヤー要素における合計歪の差に基いて当該工程後の平面展開モデルにおける各レイヤー要素の応力を算出するレイヤー要素応力算出手段と、該算出手段により算出した応力を用いて各要素の厚み方向の応力分布を作成する応力分布作成手段と、該作成手段で作成した応力分布に基いて当該工程後のスプリングバックの評価を行うスプリングバック評価手段を行うスプリングバック評価手段とを有することを特徴とする請求項１から請求項３のいずれかに記載の多段階成形シミュレーションシステム。
【請求項７】
所定の工程で展開モデル作成手段により得られた平面展開モデルとその工程前の平面展開モデルとの対応する位置にそれぞれ基準点を設定し、これらの基準点を用いて両平面展開モデルの位置合せを行う位置合せ手段と、両平面展開モデルの輪郭形状に基いてトリミングラインを設定するトリミングライン設定手段を有することを特徴とする請求項１から請求項４のいずれかに記載の多段階成形シミュレーションシステム。
【請求項８】
最終設計形状及び最終設計形状に至る１又は複数の中間設計形状を入力する入力装置と、入力された複数の設計形状について各種処理を行う処理装置と、該処理装置による処理結果を出力する出力装置とが備えられ、上記入力装置で入力された設計形状を有限要素分割して解析モデルを作成する解析モデル作成機能を有するコンピュータに用いられ、板材をブランクとする複数工程による成形結果を有限要素法を用いてシミュレーションするためのプログラムであって、該コンピュータを、上記解析モデル作成機能により作成した解析モデルの各要素を所定平面に平行になるようにその中心回りに回転させる要素回転手段、該回転手段で回転された各要素を上記所定平面上に投影させた状態で、上記解析モデルで同一節点を共有していた節点同士が一致するように応力の釣合を保ちながら各要素を変形させて解析モデルの平面展開モデルを作成する展開モデル作成手段、該作成手段で作成した平面展開モデルの各要素についてそれぞれの変形に応じた応力、歪、又は厚みの少なくとも１つを算出する演算手段、各設計形状について同様に作成した平面展開モデルの各要素の応力、歪、又は厚みの少なくとも１つに基いて各工程による成形結果を評価する評価手段として機能させることを特徴とする多段階成形シミュレーション用プログラム。
【請求項９】
最終設計形状及び最終設計形状に至る１又は複数の中間設計形状を入力する入力装置と、入力された複数の設計形状について各種処理を行う処理装置と、該処理装置による処理結果を出力する出力装置とが備えられたコンピュータに用いられ、板材をブランクとする複数工程による成形結果を有限要素法を用いてシミュレーションするためのプログラムであって、該コンピュータを、上記入力装置で入力された設計形状を有限要素分割して解析モデルを作成する解析モデル作成手段、該作成手段により作成した解析モデルの各要素を所定平面に平行になるようにその中心回りに回転させる要素回転手段、該回転手段で回転された各要素を上記所定平面上に投影させた状態で、上記解析モデルで同一節点を共有していた節点同士が一致するように応力の釣合を保ちながら各要素を変形させて解析モデルの平面展開モデルを作成する展開モデル作成手段、該作成手段で作成した平面展開モデルの各要素についてそれぞれの変形に応じた応力、歪、又は厚みの少なくとも１つを算出する演算手段、各設計形状について同様に作成した平面展開モデルの各要素の応力、歪、又は厚みの少なくとも１つに基いて各工程による成形結果を評価する評価手段として機能させることを特徴とする多段階成形シミュレーション用プログラム。
【請求項１０】
コンピュータを、中間設計形状を得る少なくとも１つの工程にその中間設計形状として複数の設計形状が設定されている場合に、この複数の中間設計形状をそれぞれ経由して最終設計形状に至る複数の成形経路を設定する成形経路設定手段として機能させると共に、演算手段として機能させるときは、各経路ごとに経路を構成する各設計形状について平面展開モデルの各要素の応力、歪、又は厚みの少なくとも１つを算出する一方、それらの算出結果に基いて最適な成形経路を判定する判定手段として機能させることを特徴とする請求項８または請求項９に記載の多段階成形シミュレーション用プログラム。
【請求項１１】
コンピュータを、平面展開モデルの各要素ごとに厚み方向に複数のレイヤー要素を設定するレイヤー要素設定手段、所定の工程で展開モデル作成手段により得られた平面展開モデルとその工程前の平面展開モデルとの対応する位置にそれぞれ基準点を設定し、これらの基準点を用いて両平面展開モデルの位置合せを行う位置合せ手段、当該工程前後の設計形状の曲率から平面展開モデルにおける各レイヤー要素の曲げ歪を算出する曲げ歪算出手段、レイヤー要素ごとに、演算手段を用いて算出した歪と上記曲げ歪算出手段により算出した曲げ歪とを加算して合計歪を算出する合計歪算出手段、当該工程前後の平面展開モデルの互いに対応するレイヤー要素における合計歪の差に基いて当該工程後の平面展開モデルにおける各レイヤー要素の応力を算出するレイヤー要素応力算出手段、該算出手段により算出した応力を用いて各要素の厚み方向の応力分布を作成する応力分布作成手段、該作成手段で作成した応力分布に基いて当該工程後のスプリングバックの評価を行うスプリングバック評価手段を行うスプリングバック評価手段として機能させることを特徴とする請求項８または請求項９に記載の多段階成形シミュレーション用プログラム。
【請求項１２】
コンピュータを、所定の工程で展開モデル作成手段により得られた平面展開モデルとその工程前の平面展開モデルとの対応する位置にそれぞれ基準点を設定し、これらの基準点を用いて両平面展開モデルの位置合せを行う位置合せ手段、両平面展開モデルの輪郭形状に基いてトリミングラインを設定するトリミングライン設定手段として機能させることを特徴とする請求項８または請求項９に記載の多段階成形シミュレーション用プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、板状ブランクを用いた成形の結果を有限要素法を用いてシミュレーションするシステム、より詳しくは、複数工程による各中間形状を経て最終製品形状に成形する多段階成形をシミュレーションする多段階成形シミュレーションシステム及び多段階成形シミュレーション用プログラムに関し、コンピュータによる有限要素解析技術の分野に属する。
【背景技術】
【０００２】
板材をブランクとし、これをプレス成形や折り曲げ成形等により所定形状の製品に成形する場合、まず、設計された製品形状を試作、検証した後、これを量産するための金型を設計、製作し、次いでこの金型を用いて製品を試作し、ＯＫであれば量産に移行する、というプロセスを経ることになるが、このようなプロセスにおいては、金型を用いた製品の試作の段階で、しわ、われ等の欠陥の発生が判明することがあり、この場合、製品の設計形状の見直しから金型の修正、試作まで、再度上記プロセスを実行しなければならないことになり、多大なコストと時間を費やすことになる。
【０００３】
このような実情に対し、近年、製品設計の段階、或は金型設計の段階で、成形結果をシミュレーションし、金型製造前のできるだけ早い段階で欠陥の発生を予測し、事前にその対策をとることにより上記のような多大な無駄を回避することがおこなわれており、そのようなシミュレーションシステムとして各種のものが提供されている。
【０００４】
これらのシステムは、いずれも有限要素法を用いるものであって、動的陽解法、静的陽解法、静的陰解法、ワンステップ法などの手法を用いたものがある。このうち、例えば特許文献１に開示されているようなワンステップ法は、製品の設計形状を平板状に展開したときのブランクの形状をワンステップで求めると共に、その際に各要素に生じる歪や応力から製品各部の板厚の分布やしわ、われ等の欠陥の発生を評価するものであり、計算量が少ないことや、金型に関する情報が不要であるなどの長所がある。 
【特許文献１】特開２００４−１４８３８１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
ところで、板状のブランクを製品形状に成形する際に、単一工程では製品形状を物理的に成形できないことがあり、このような場合には複数工程を経て製品形状を実現する多段階成形が実施されることになる。この多段階成形では、一般に、製品の設計者により提示された製品形状に基いて、金型の設計者がこれを実現するために各工程で成形すべき中間形状を設計すると共に中間形状を成形するための金型を設計することになる。そして、各工程の金型を製造し、製品を試作した段階で、割れやしわ等の欠陥の発生が判明した場合、各中間形状の見直しや各工程の金型の修正を行うとすると、膨大なコストと時間を費やすことになる。
【０００６】
一方、このような無駄を事前に回避するために多段階成形のシミュレーションを行う場合は、まず、各中間形状及び最終形状、並びに各工程で使用される金型等の外部の部材などをモデル化し、これらのモデルを用いてブランクの初期形状から最終形状に至る加工プロセスをシミュレーションし、事前にしわや割れ等の欠陥の発生を予測することが考えられる。しかしながら、このような方法では、多くのモデル化の作業が必要な上、扱うデータ量及び計算量が多くなり、データの管理が煩雑となると共にコンピュータのパワーを要求するといった問題が生じ、さらに、中間形状の設計案を多数用意してより最適となる形状の獲得を期すほどこれらの問題が増幅されることになる。このように、従来の方法では多段階成形のシミュレーションは事実上不可能であったといえる。また、最終形状のみを対象としたワンステップ法により、最終形状で生じる歪や応力からしわや割れ等の欠陥の発生の評価を行うことも考えられるが、最終形状の評価がＯＫでも、中間工程における歪や厚さの増減を考慮すると、中間工程の評価がＯＫであるとは限らない。
【０００７】
一方、このような板材の曲げ成形においては、板材の屈曲部の外側では材料が伸びるために離型後に縮みの応力が、内側では逆に材料が縮むために伸びの応力が発生し、離型後に板厚方向の応力差を低減するようにモーメントが生じて形状にスプリングバック等の曲率変化を生じることがある。曲げ方向と反対の方向への曲率変化をスプリングバック、曲げ方向への曲率変化をスプリングゴーといい（以下、これらを単に「スプリングバック」という）、これらによって形状の寸法精度不良を招くことになる。この寸法精度不良は、部品同士の組み付け性や外観を損なう他、多段階成形にあっては、次工程の金型に対して適正な位置に配備できずに加工不良を生じる可能性がある。そのため、シミュレーションにより事前にスプリングバックの発生を予測して、シミュレーションの結果に基いて実際の成形において余分に曲げるなどの対策をとることが考えられるが、従来では、前述の多段階成形のシミュレーションと同様に、各設計形状とや周囲の部品（ダイ、パンチ等）を全てモデル化し、接触条件を付加して実物の行程に沿ってシミュレーションする方法がとられていたため、データ量及び計算量が多くなるという問題を有していた。
【０００８】
そこで、本発明は、上記のような実情に鑑み、多段階成形のシミュレーションの簡素化及び高速化を図ると共に、スプリングバックのシミュレーションをも同時に行うことができるシステムを提供することを課題とする。
【課題を解決するための手段】
【０００９】
上記課題を解決するため、本発明は次のように構成したことを特徴とする。
【００１０】
まず、請求項１に記載の発明は、板材をブランクとする複数工程による成形結果を有限要素法を用いてシミュレーションするシステムであって、最終設計形状及び最終設計形状に至る１又は複数の中間設計形状を入力する入力装置と、入力された複数の設計形状について各種処理を行う処理装置と、該処理装置による処理結果を出力する出力装置とが備えられていると共に、上記処理装置は、上記入力装置で入力された設計形状を有限要素分割して解析モデルを作成する解析モデル作成手段と、該作成手段で作成した解析モデルの各要素を所定平面に平行になるようにその中心回りに回転させる要素回転手段と、該回転手段で回転された各要素を上記所定平面上に投影させた状態で、上記解析モデルで同一節点を共有していた節点同士が一致するように応力の釣合を保ちながら各要素を変形させて解析モデルの平面展開モデルを作成する展開モデル作成手段と、該作成手段で作成した平面展開モデルの各要素についてそれぞれの変形に応じた応力、歪、又は厚みの少なくとも１つを算出する演算手段と、各設計形状について同様に作成した平面展開モデルの各要素の応力、歪、又は厚みの少なくとも１つに基いて各工程による成形結果を評価する評価手段とを有することを特徴とする。
【００１１】
また、請求項２に記載の発明は、上記請求項１に記載のシステムにおいて、要素回転手段は、要素中心の法線ベクトルと要素中心を通って所定平面に直交する所定平面直交ベクトルとを作成し、上記要素中心の法線ベクトルが上記所定平面直交ベクトルに一致するように各要素を回転させることを特徴とする。
【００１２】
そして、請求項３に記載の発明は、上記請求項１または請求項２に記載のシステムにおいて、展開モデル作成手段は、回転手段で回転された各要素を所定平面上に投影させた状態で、要素回転前の解析モデルで同一節点を共有していた節点同士に各節点間距離に応じた引っ張り荷重を設定し、この引っ張り荷重により解析モデルで同一節点を共有していた節点同士が引っ張り荷重の釣合を保ちながら一致するように各要素を変形させることを特徴とする。
【００１３】
一方、請求項４に記載の発明は、上記請求項１から請求項３のいずれかに記載のシステムにおいて、所定の工程で展開モデル作成手段により得られた平面展開モデルとその工程前の平面展開モデルの対応する位置にそれぞれ基準点を設定し、これらの基準点を用いて両平面展開モデルの位置合せを行う位置合せ手段が備えられ、演算手段は、該位置合せ手段によって位置合せされた当該工程前後の平面展開モデルの互いに対応する要素における歪の差を演算し、その差に基いて、当該工程後の平面展開モデルにおける各要素の応力を計算することを特徴とする。
【００１４】
また、請求項５に記載の発明は、上記請求項１から請求項４のいずれかに記載のシステムにおいて、中間設計形状を得る少なくとも１つの工程にその中間設計形状として複数の設計形状が設定されている場合に、この複数の中間設計形状をそれぞれ経由して最終設計形状に至る複数の成形経路を設定する成形経路設定手段が備えられ、演算手段は、各経路ごとに経路を構成する各設計形状について平面展開モデルの各要素の応力、歪、又は厚みの少なくとも１つを算出すると共に、それらの算出結果に基いて最適な成形経路を判定する判定手段が備えられていることを特徴とする。
【００１５】
一方、請求項６に記載の発明は、上記請求項１から請求項３のいずれかに記載のシステムにおいて、平面展開モデルの各要素ごとに厚み方向に複数のレイヤー要素を設定するレイヤー要素設定手段と、所定の工程で展開モデル作成手段により得られた平面展開モデルとその工程前の平面展開モデルとの対応する位置にそれぞれ基準点を設定し、これらの基準点を用いて両平面展開モデルの位置合せを行う位置合せ手段と、当該工程前後の設計形状の曲率から平面展開モデルにおける各レイヤー要素の曲げ歪を算出する曲げ歪算出手段と、レイヤー要素ごとに、演算手段を用いて算出した歪と上記曲げ歪算出手段により算出した曲げ歪とを加算して合計歪を算出する合計歪算出手段と、当該工程前後の平面展開モデルの互いに対応するレイヤー要素における合計歪の差に基いて当該工程後の平面展開モデルにおける各レイヤー要素の応力を算出するレイヤー要素応力算出手段と、該算出手段により算出した応力を用いて各要素の厚み方向の応力分布を作成する応力分布作成手段と、該作成手段で作成した応力分布に基いて当該工程後のスプリングバックの評価を行うスプリングバック評価手段を行うスプリングバック評価手段とを有することを特徴とする。
【００１６】
また、請求項７に記載の発明は、上記請求項１から請求項３のいずれかに記載のシステムにおいて、所定の工程で展開モデル作成手段により得られた平面展開モデルとその工程前の平面展開モデルとの対応する位置にそれぞれ基準点を設定し、これらの基準点を用いて両平面展開モデルの位置合せを行う位置合せ手段と、両平面展開モデルの輪郭形状に基いてトリミングラインを設定するトリミングライン設定手段を有することを特徴とする。
【００１７】
ここで、上記請求項１〜７に記載の発明の概念を図１を用いて説明する。
【００１８】
図１（ａ）は、多段階成形における中間形状と最終形状とを設計したものである。図１（ｂ）は、これらの入力データから、中間形状の解析モデルＭ１及び最終形状の解析モデルＭ２を作成したものである。さらに、図１（ｃ）は、解析モデルＭ１を展開計算して得られる平面展開モデルＯ１、及び解析モデルＭ２を展開計算して得られる平面展開モデルＯ２とを作成したものである。このとき、各要素の変形に応じた応力、歪、又は厚みを計算することができ、これらの計算結果に基いて各工程による成形結果の評価、即ち各設計形状が適切か否かを評価することができる。
【００１９】
一方、両平面展開モデルＯ１、Ｏ２に、それぞれ対応する代表的な位置に基準点Ｐ１、Ｐ２を設定し、この基準点Ｐ１、Ｐ２が一致するように両モデルＯ１、Ｏ２を重ね合わせて位置合せを行う。このとき、両モデルＯ１、Ｏ２の位置関係に基いて、各要素の位置関係を特定することができ、例えば解析モデルＭ２の点Ａ２に対応する平面展開モデルＯ２上の点Ａ２′は、平面展開モデルＯ１上においてＡ１′に位置することを特定することができる。
【００２０】
そして、図２に示すように、両平面展開モデルＯ１、Ｏ２の輪郭形状に基いて、点Ａ１′の位置を解析モデルＭ１に反映させ、同様に平面展開モデルＯ１における平面展開モデルＯ２上の全ての点について解析モデルＭ１に反映させることによって、トリミングラインを見出すことができる。
【００２１】
一方、請求項８に記載の発明は、最終設計形状及び最終設計形状に至る１又は複数の中間設計形状を入力する入力装置と、入力された複数の設計形状について各種処理を行う処理装置と、該処理装置による処理結果を出力する出力装置とが備えられ、上記入力装置で入力された設計形状を有限要素分割して解析モデルを作成する解析モデル作成機能を有するコンピュータに用いられ、板材をブランクとする複数工程による成形結果を有限要素法を用いてシミュレーションするためのプログラムであって、該コンピュータを、上記解析モデル作成機能により作成した解析モデルの各要素を所定平面に平行になるようにその中心回りに回転させる要素回転手段、該回転手段で回転された各要素を上記所定平面上に投影させた状態で、上記解析モデルで同一節点を共有していた節点同士が一致するように応力の釣合を保ちながら各要素を変形させて解析モデルの平面展開モデルを作成する展開モデル作成手段、該作成手段で作成した平面展開モデルの各要素についてそれぞれの変形に応じた応力、歪、又は厚みの少なくとも１つを算出する演算手段、各設計形状について同様に作成した平面展開モデルの各要素の応力、歪、又は厚みの少なくとも１つに基いて各工程による成形結果を評価する評価手段として機能させることを特徴とする。
【００２２】
また、請求項９に記載の発明は、最終設計形状及び最終設計形状に至る１又は複数の中間設計形状を入力する入力装置と、入力された複数の設計形状について各種処理を行う処理装置と、該処理装置による処理結果を出力する出力装置とが備えられたコンピュータに用いられ、板材をブランクとする複数工程による成形結果を有限要素法を用いてシミュレーションするためのプログラムであって、該コンピュータを、上記入力装置で入力された設計形状を有限要素分割して解析モデルを作成する解析モデル作成手段、該作成手段により作成した解析モデルの各要素を所定平面に平行になるようにその中心回りに回転させる要素回転手段、該回転手段で回転された各要素を上記所定平面上に投影させた状態で、上記解析モデルで同一節点を共有していた節点同士が一致するように応力の釣合を保ちながら各要素を変形させて解析モデルの平面展開モデルを作成する展開モデル作成手段、該作成手段で作成した平面展開モデルの各要素についてそれぞれの変形に応じた応力、歪、又は厚みの少なくとも１つを算出する演算手段、各設計形状について同様に作成した平面展開モデルの各要素の応力、歪、又は厚みの少なくとも１つに基いて各工程による成形結果を評価する評価手段として機能させることを特徴とする。
【００２３】
そして、請求項１０に記載の発明は、上記請求項８または請求項９に記載のプログラムにおいて、コンピュータを、中間設計形状を得る少なくとも１つの工程にその中間設計形状として複数の設計形状が設定されている場合に、この複数の中間設計形状をそれぞれ経由して最終設計形状に至る複数の成形経路を設定する成形経路設定手段として機能させると共に、演算手段として機能させるときは、各経路ごとに経路を構成する各設計形状について平面展開モデルの各要素の応力、歪、又は厚みの少なくとも１つを算出する一方、それらの算出結果に基いて最適な成形経路を判定する判定手段として機能させることを特徴とする。
【００２４】
また、請求項１１に記載の発明は、上記請求項８または請求項９に記載のプログラムにおいて、コンピュータを、平面展開モデルの各要素ごとに厚み方向に複数のレイヤー要素を設定するレイヤー要素設定手段、所定の工程で展開モデル作成手段により得られた平面展開モデルとその工程前の平面展開モデルとの対応する位置にそれぞれ基準点を設定し、これらの基準点を用いて両平面展開モデルの位置合せを行う位置合せ手段、当該工程前後の設計形状の曲率から平面展開モデルにおける各レイヤー要素の曲げ歪を算出する曲げ歪算出手段、レイヤー要素ごとに、演算手段を用いて算出した歪と上記曲げ歪算出手段により算出した曲げ歪とを加算して合計歪を算出する合計歪算出手段、当該工程前後の平面展開モデルの互いに対応するレイヤー要素における合計歪の差に基いて当該工程後の平面展開モデルにおける各レイヤー要素の応力を算出するレイヤー要素応力算出手段、該算出手段により算出した応力を用いて各要素の厚み方向の応力分布を作成する応力分布作成手段、該作成手段で作成した応力分布に基いて当該工程後のスプリングバックの評価を行うスプリングバック評価手段を行うスプリングバック評価手段として機能させることを特徴とする。
【００２５】
さらに、請求項１２に記載の発明は、上記請求項８または請求項９に記載のプログラムにおいて、コンピュータを、所定の工程で展開モデル作成手段により得られた平面展開モデルとその工程前の平面展開モデルとの対応する位置にそれぞれ基準点を設定し、これらの基準点を用いて両平面展開モデルの位置合せを行う位置合せ手段、両平面展開モデルの輪郭形状に基いてトリミングラインを設定するトリミングライン設定手段として機能させることを特徴とする。
【発明の効果】
【００２６】
まず、請求項１〜３に記載の発明によれば、多段階成形において、解析モデル作成手段、要素回転手段、展開モデル作成手段によって、各設計形状についての平面展開モデルを同様に繰り返し作成し、各平面展開モデルの各要素の変形に基いて応力、歪、又は厚みの少なくとも１つを求め、各要素の応力、歪、又は厚みに基いて各工程によるしわや割れ等の欠陥の発生を評価することができる。このように、各中間形状をあたかも最終形状のように扱うことによって、欠陥の発生する可能性の高い工程を特定することができると共に、得られた応力、歪、又は厚みの情報に基いて中間形状の設計変更等を行うなどの対策をとることができる。
【００２７】
このように、本発明においては、各設計形状について作成された解析モデルからワンステップで平面展開モデルを作成し、展開計算結果の各要素の変形に基づいて計算された各要素の応力、歪、又は厚みのデータにより、多段階成形のシミュレーションに使用する大方のデータが揃うことになる。従って、従来のように金型や周囲の部材等をモデル化する必要がないので扱うデータ量が少なく、各平面展開モデルの各要素についての応力、歪、厚みの値を組み合わせるだけで簡単に各工程での成形結果を評価することができ、圧倒的な計算の簡素化及び高速化を実現することができる。
【００２８】
また、各要素を回転させるときには、具体的には請求項２に記載の発明のように、まず、要素中心の法線ベクトルと要素中心を通って所定平面に直交する所定平面直交ベクトルとを作成し、上記要素中心の法線ベクトルが上記所定平面直交ベクトルに一致するように各要素を回転させることによって、各要素を所定平面に平行にすることができる。
【００２９】
そして、応力の釣り合いを保ちながら各要素を変形させるときは、請求項３に記載の発明のように、回転後の各要素を所定平面に投影させた状態で、要素回転前の解析モデルで同一節点を共有していた節点同士に各節点間距離に応じた引っ張り荷重を設定し、この引っ張り荷重により解析モデルで同一節点を共有していた節点同士を一致させるように計算を行うことによって、応力の釣合を保ちながら各要素を変形させることができる。
【００３０】
そして、平面展開モデルの各要素について応力を計算するときは、具体的には請求項４に記載の発明のように、工程前後の平面展開モデルの互いに対応する要素における歪の値を演算し、その差に基いて、当該工程後の平面展開モデルにおける各要素の応力を計算することができる。
【００３１】
また、請求項５に記載の発明によれば、金型の設計者が各工程で得られる中間形状の複数の設計案を用意する場合、各設計案における各要素の応力、歪、又は厚みに基いて、最適な成形経路、即ち最適となる中間設計形状の組合せを判定することができる。
【００３２】
このように、複数の成形経路から最適なものを選択する際も、各設計形状ごとにワンステップ法により得られた応力、歪、又は厚みのデータを組合せるだけでよく、扱うデータ量が少なく、計算の簡素化及び高速化が実現される一方、多くの中間設計形状の設計案を用意してより適当となる中間設計形状を得ることができる。
【００３３】
また、請求項６に記載の発明は、各要素を分割して求めた厚み方向の応力分布によりスプリングバックの評価を行うものである。即ち、平面展開モデル作成時に得られる各要素の歪は、要素の厚み方向に対して一様であり、この歪を用いて算出した応力は、同様に厚み方向に対して一様であるため、スプリングバックの評価には使えないことになるが、各要素の厚み方向に複数のレイヤー要素を設定して各レイヤー要素において、展開計算で得られた歪と各設計形状の曲率に基いて計算された曲げ歪とを加算して合計歪とし、工程前後の平面展開モデルの対応するレイヤー要素の合計歪の差を算出して、この差に基いて各レイヤー要素の応力が計算できると共に、各要素の厚み方向の曲げ歪を考慮した応力分布が得られ、スプリングバックを適正に評価することができる。そして、金型の設計者は、厚み方向の応力分布に基いて金型を設計することができ、例えば加工の際にスプリングバック量を見込んで余分に曲げるように金型を設計する等の対策をとることができ、工程により得られる形状の寸法精度を向上させることができる。このように、ワンステップ法により得られたデータを組み合わせてスプリングバックの評価を行うことができるので、従来のシステムに比べて計算の簡素化及び高速化を実現することができる。
【００３４】
さらに、請求項７に記載の発明によれば、工程前後の平面展開モデルの輪郭形状に基いて、各要素の対応関係を工程前の解析モデルに反映させることができ、トリミングラインを見出すことができる。これによって、金型の設計者は、各工程の中間設計形状を示すだけで、自動的にトリミングラインを設定することができる。
【００３５】
一方、請求項８、９に記載の発明は、上記請求項１に記載のシステムのプログラムに対応するもので、請求項８は、既存のプログラムとして解析モデル作成機能が予めインストールされたコンピュータ、請求項９は、解析モデル作成機能を有しないコンピュータに使用するプログラムであって、いずれのプログラムにおいても請求項１に記載のシステムと同様の作用が得られる。
【００３６】
また、請求項１０〜１２に記載の発明は、上記請求項５〜７に記載のシステムのプログラムに対応するもので、それぞれ請求項５〜７に記載のシステムと同様の作用が得られる。
【発明を実施するための最良の形態】
【００３７】
次に、本発明の実施の形態に係る多段階成形シミュレーションシステムについて説明する。
【００３８】
図３は、このシステムの中心となるコンピュータの構成を示すもので、このコンピュータ１０は、中央処理装置１１と、各種条件の設定やシステムの制御等に用いられる入力装置１２と、記録媒体２０からのデータ読み込み装置１３と、プログラム記録部１４ａ及びデータファイル記録部１４ｂが設けられた記録装置１４と、入力画面や計算結果等を表示する表示装置１５と、計算結果等を印刷する印刷装置１６とを有する。
【００３９】
上記記録装置１４のプログラム記録部１４ａには、図４に示すように、メインプログラムと複数のサブプログラムとが記録され、サブプログラムとしては、解析モデル作成サブプログラム、解析用基礎データ作成サブプログラム、解析用データ作成プログラム、解析計算サブプログラム、マッピング処理サブプログラム、描画サブプログラム、評価用データ作成サブプログラム、及び成形経路評価サブプログラムが記録されている。
【００４０】
また、上記記録装置１４は、入力するＣＡＤデータ即ち設計形状に応じた複数のデータファイル記録部１４ｂ…１４ｂを有し、各データファイル記録部１４ｂには、図５に示すように、当該ＣＡＤデータに係る解析モデルデータファイルと、解析用基礎データファイルと、材料属性データファイルと、解析用データファイルと、解析結果データファイルと、評価用データファイルとがそれぞれ記録されるようになっている。
【００４１】
これらのデータファイルのうち、解析モデルデータファイルは、一部を図６に示すような解析対象製品を有限要素分割してなる解析モデルについての各種のデータを記録したテーブルで構成され、このうち、図７に示す要素構成テーブルは、各要素Ｅ1、Ｅ２…に対してこれを構成する節点番号Ｎ１、Ｎ２…をそれぞれ記録し、例えば図６に示す要素Ｅ１の場合は、４節点Ｎ１、Ｎ２、Ｎ５、Ｎ４で構成されていることを示している。また、図８に示す節点座標テーブルは、図７のテーブルに記録された各要素Ｅ１、Ｅ２…を構成する各節点Ｎ１、Ｎ２…のｘ座標、ｙ座標、ｚ座標をそれぞれ記録するようになっている。
【００４２】
また、解析用基礎データファイルは、解析計算に先立ち、予め計算され或は設定された各種のデータを記録したテーブルで構成されている。
【００４３】
図９に示す節点構成テーブルは、前述の解析モデルデータファイルに含まれる要素構成データを図７のテーブルから読み出し、各節点を構成する要素、即ち当該節点を取り囲む要素の番号Ｅ１、Ｅ２…を記録するようになっている。例えば図６の節点Ｎ５については、これを取り囲む節点の番号Ｅ１〜Ｅ４が記録される。
【００４４】
図１０に示す要素中心座標テーブルは、各要素の中心座標を記録するようになっている。即ち、図７の要素構成テーブルにより例えば要素Ｅ１を構成する節点Ｎ１、Ｎ２…を読み出すと共に、図８の節点座標テーブルからこれらの節点Ｎ１、Ｎ２…のｘ座標、ｙ座標、ｚ座標を読み出し、これら節点のｘ座標、ｙ座標、ｚ座標のそれぞれの平均座標が当該要素Ｅ１の中心座標として記録される。
【００４５】
図１１に示す要素ベクトルテーブルは、各要素についての法線ベクトル（以下、「要素ベクトル」という）ＶＥ１、ＶＥ２…のｘ成分、ｙ成分、ｚ成分を記録するようになっている。
【００４６】
この要素ベクトルＶＥ１、ＶＥ２…の算出方法を例えば図６に示す要素Ｅ１について説明すると、まず図７の要素構成テーブルから要素Ｅ１が節点Ｎ１、Ｎ２、Ｎ５、Ｎ４で構成されることを読み出すと共に、これらの節点Ｎ１、Ｎ２、Ｎ５、Ｎ４のｘ座標、ｙ座標、ｚ座標を図８の節点座標テーブルから読み出し、これらの座標データに基づき、要素Ｅ１を平面としたときの該要素Ｅ１の傾きを求める。そして、図１０の要素中心座標テーブルから各要素の中心座標を読み出し、要素中心でその平面に直交する方向の単位ベクトルのｘ成分、ｙ成分、ｚ成分を求め、これらの成分を有するベクトルを要素Ｅ１の要素ベクトルＶＥ１とする。
【００４７】
図１２に示す要素回転軸ベクトルテーブルは、要素回転の際の回転軸となる回転軸ベクトルのｘ成分、ｙ成分、ｚ成分を各要素ごとに記録するようになっている。
【００４８】
この回転軸ベクトルＶＳＥ１、ＶＳＥ２…の算出方法を図１３に基いて説明すると、まず、図１０の要素中心座標テーブルから各要素の中心座標を読み出し、この中心座標を通りｘｙ平面に直交するｘｙ平面直交ベクトルＶｚを生成する。そして、図１１の要素ベクトルテーブルから要素ベクトルＶＥ１、ＶＥ２…を読み出し、上記ｘｙ平面直交ベクトルＶｚと要素ベクトルＶＥ１、ＶＥ２…とで決定される平面Ｒを定義する。次に、中心座標を通り上記平面Ｒに直交する方向の単位ベクトルを要素回転軸ベクトルＶＳＥ１、ＶＳＥ２…として、該ベクトルＶＳＥ１、ＶＳＥ２…のｘ成分、ｙ成分、ｚ成分が各要素ごとに図１２の要素回転軸ベクトルテーブルに記録される。
【００４９】
図１４に示す要素回転角度テーブルは、要素ベクトルＶＥ１、ＶＥ２…を要素中心回りに回転させたときに、上記ｘｙ平面直交ベクトルＶｚに一致するために必要な回転角度θを各要素ごとに記録するようになっている。
【００５０】
図１５の分離要素構成テーブルは、番号ＥＳ１、ＥＳ２…を付した回転後の各要素（分離要素）において、各分離要素を構成する各節点に番号ＮＳ１、ＮＳ２…を付し、各分離要素ごとにこれを構成する節点番号ＮＳ１、ＮＳ２…を記録するようになっている。つまり、要素回転により各要素は離反し、これに伴って節点Ｎ１、Ｎ２…は分離するので、この分離した後の要素（分離要素）を構成する節点番号ＮＳ１、ＮＳ２…が記録される。
【００５１】
図１６の元節点−分離要素の節点テーブルは、分離要素の節点ＮＳ１、ＮＳ２…が、要素回転前の節点Ｎ１、Ｎ２…（元節点）のいずれに対応するかを示したテーブルであって、要素回転前の元節点Ｎ１、Ｎ２…ごとに要素回転前に共有していた分離要素の節点番号ＮＳ１、ＮＳ２…を記録するようになっている。なお、各節点Ｎ１、Ｎ２…に対応する分離要素の節点ＮＳ１、ＮＳ２…の数は、有限要素分割の結果に応じて増減するようになっている。
【００５２】
図１７の分離要素の節点座標テーブルは、分離要素の節点ＮＳ１、ＮＳ２…ごとにｘ座標、ｙ座標を記録するようになっている。
【００５３】
図１８の節点ベクトルテーブルは、各節点ごとに、後述する曲げ歪を計算するために必要なデータである節点ベクトルＶＮ１、ＶＮ２…のｘ成分、ｙ成分、ｚ成分を記録するようになっている。
【００５４】
ここで、この節点ベクトルＶＮ１、ＶＮ２…の算出方法を例えば図１９に示す節点Ｎ５について説明すると、図９の節点構成テーブルから節点Ｎ５を構成する要素、即ち節点Ｎ５を取り囲む要素を読み出す。この場合、図１９の例では、要素Ｅ１、Ｅ２、Ｅ４、Ｅ３が読み出される。
【００５５】
そして、図１１の要素ベクトルテーブルから、これらの要素Ｅ１、Ｅ２、Ｅ４、Ｅ３についての要素ベクトルＶＥ１、ＶＥ２、ＶＥ４、ＶＥ３のｘ成分、ｙ成分、ｚ成分をそれぞれ読み出すと共に、これら４つのベクトルの各成分の平均値をｘ成分、ｙ成分、ｚ成分とする単位ベクトルを求め、これを節点Ｎ５についての節点ベクトルＶＮ５とする。
【００５６】
このようにして、各節点Ｎ１、Ｎ２…についての節点ベクトルＶＮ１、ＶＮ２…を順次求め、これを図１８の節点ベクトルテーブルに書き込む。
【００５７】
図２０の拘束条件テーブルは、後述する展開計算に際して解析モデルの位置決めをするために指定した２つ以上の節点番号ＮＡ、ＮＢ…が記録されている。そして、この拘束点ＮＡ、ＮＢは、後述する位置合せの際の基準点となる。
【００５８】
図２１の材料特性テーブルは、解析対象の製品に関するデータ、即ち当該材料のヤング率、ポアソン比、降伏応力等の解析に必要なデータを記録するようになっている。
【００５９】
図２２の解析結果テーブルは、展開計算による各分離要素変形後の各節点ＮＳ１、ＮＳ２…のｘ座標、ｙ座標をそれぞれ記録するようになっている。
【００６０】
図２３の歪テーブルは、解析計算の結果生じた各分離要素ＥＳ１、ＥＳ２…の変形に基いて、設定した評価点ごとにｘ方向の歪εｘ、ｙ方向の歪εｙ、ｚ方向の歪εｚ、ｘ方向にみたｙ方向の剪断歪εｘｙ、同様に剪断歪εｙｚ、εｚｘ、及び各要素ごとに板厚減少率を記録するようになっている。ここで、上記評価点について説明すると、要素Ｅを図２４に示すような正規化座標（ｒ，ｓ）に変換し、ガウス積分点とよばれる要素内部の座標（±１／√３，±１／√３）に対応する４点を評価点Ｄａ〜Ｄｄとし、さらに、図２５に示すように要素Ｅを７層のレイヤーＬ１〜Ｌ７に分割し、各レイヤーごとに４つの評価点ＤＬ１ａ〜ＤＬ１ｄ、…ＤＬ７ａ〜ＤＬ７ｄを設定している。
【００６１】
図２６の応力テーブルは、各評価点ごとに応力σｘｘ、σｙｙ、σｚｚ及び剪断応力σｘｙ、σｙｚ、σｚｘを記録するようになっている。ここに記録された応力の値は、図２３の歪テーブルから読み出した歪の値に基いて、図２１の材料特性テーブルのヤング率やポアソン比等のデータを適宜用いて計算されたものである。
【００６２】
また、データファイル記録部１４ｂの評価用データファイルに記録された評価用データファイルには、各工程の安全性評価に使用される各種テーブルが記録されている。
【００６３】
図２７の評価点マッピングテーブルは、後述する工程前後の平面展開モデルの位置合せを行ったときの、当該工程前後の平面展開モデルの各要素の対応関係を表すものである。このテーブルは、工程後の各評価点ごとに、該評価点の座標を含む工程前の要素番号を記録すると共に、工程前の各要素のｒｓ座標系における評価点の座標（ｒ，ｓ）を記録するようになっている。そして、これらが記録されたテーブルが、成形経路を構成する各工程ごとに用意されている。また、この評価点マッピングテーブルは、後述する成形経路ごとに工程数ｍ−１（工程２後〜最終工程ｍ後）分をワンセットとして準備されている。
【００６４】
図２８の補間歪テーブルは、図２７の評価点マッピングテーブルに記録された工程前の各要素における座標（ｒ，ｓ）での各歪成分を記録するようになっている。各歪成分は、各評価点ごとに、歪εｘ、εｙ、εｚ、剪断歪εｘｙ、εｙｚ、εｚｘ、主歪ε１、ε２、及び板厚減少率である。なお、各要素の座標（ｒ，ｓ）における歪は、展開計算の要素の変形に伴う各評価点の変位により求めることができる。また、上記評価点マッピングテーブルと同様に、これらが記録されたテーブルが、成形経路を構成する各工程ごとに用意されている。また、補間歪テーブルは、上記評価点マッピングテーブルと同様に、成形経路ごとに工程数ｍ−１分をワンセットとして準備されている。
【００６５】
図２９の補間応力テーブルは、図２５の評価点マッピングテーブルに記録された工程前の各要素における座標（ｒ，ｓ）での各応力成分を記録するようになっている。各応力成分は、各評価点ごとに応力σｘｘ、σｙｙ、σｚｚ、及び剪断応力σｘｙ、σｙｚ、σｚｘである。そして、上記補間歪テーブルと同様に、これらが記録されたテーブルが、成形経路を構成する各工程ごとに用意されている。また、補間応力テーブルは、上記補間歪テーブルと同様に、成形経路ごとに工程数ｍ−１分をワンセットとして準備されている。
【００６６】
また、図３０の評価用応力テーブルは、後述する計算により求めた工程前の形状に依存した応力が各評価点ごとに記録され、成形経路ごとに工程数ｍ（工程１後〜最終工程ｍ後）分をワンセットとして準備されることになる。
【００６７】
次に、当該システムの作用を、上記記録装置１４のプログラム記録部１４ａに記録されているメインプログラムの動作を示す図３１のフローチャートに従って説明する。ここでは、自動車の車体を構成する一部材を成形するために実施される多段階成形として、板材をブランクとし、図３２に示すようなブランクの中央部を膨出させた第１中間形状と、図３３に示すような該第１中間形状をトリミングすると共に所定の加工を施した第２中間形状と、図３４に示すような該第２中間形状の各部をさらに折り曲げて成形される最終形状とを設計している。なお、最終形状において、図３５（ｂ）に示すような折り曲げ方は単一工程では実現することができないため、図３５（ａ）に示すような第２中間形状を経て複数工程で成形する多段階成形を実施するのである。
【００６８】
まず、プロセスＰ１として、図３に示すコンピュータ１０のデータ読込み装置１３により、記録媒体２０から解析対象の多段階成形に係る各工程により成形される中間形状のＣＡＤデータを読み込む。コンピュータ１０の表示装置１５には、図３２〜図３４に示すような中間形状又は最終形状の全体を表示した画面１５ａ，１５ｂ，１５ｃが表示されることになる。
【００６９】
次に、プロセスＰ２〜Ｐ８で、上記各ＣＡＤデータごとに展開計算を行う。以下、上記図３２で示した第１中間形状のＣＡＤデータの展開計算を代表して説明する。
【００７０】
まず、プロセスＰ２で、上記プログラム記録部１４ａに記録されている解析モデル作成サブプログラムにより、上記ＣＡＤデータに基いて当該形状を有限要素分割した解析モデルが作成される。このとき、図３６の画面１５ｄに解析モデルを示すように、設計形状が各要素に分割されて表示される。そして、各要素を構成する節点番号Ｎ１、Ｎ２…を図７の要素構成テーブルに書き込み、各節点のｘ座標、ｙ座標、ｚ座標を図８の節点座標テーブルに書き込む。
【００７１】
次に、プロセスＰ３として、上記解析モデルについての解析用基礎データが上記プログラム記録部１４ａに記録されている解析用基礎データ作成サブプログラムにより作成される。
【００７２】
ここで、この解析用基礎データ作成サブプログラムの動作を図３７に示すフローチャートに従って説明する。
【００７３】
まず、プロセスＰ２１で、解析モデル作製処理により得られた節点の要素構成データを図７の要素構成テーブルから読み出し、各節点Ｎ１、Ｎ２…を取り囲む要素の番号Ｅ１、Ｅ２…を図９の節点構成テーブルに書き込む。次に、プロセスＰ２２で、節点の座標データを図８の節点座標テーブルから読み出し、要素を構成する節点の座標の平均座標を算出して、この平均座標を要素の中心座標として、これらのｘ座標、ｙ座標、ｚ座標を図１０の要素中心座標テーブルに書き込む。
【００７４】
そして、プロセスＰ２３〜Ｐ２５で、要素回転のために使用する各データを作成する。即ち、プロセスＰ２３では、前述のように各要素についての要素ベクトルＶＥ１、ＶＥ２…を生成し、図１１の要素ベクトルテーブルに書き込む。
【００７５】
次に、プロセスＰ２４で、前述のように回転軸ベクトルＶＳＥ１、ＶＳＥ２…のｘ成分、ｙ成分、ｚ成分を生成し、図１２の要素回転軸テーブルに書き込む。また、プロセス２５で図１４の要素回転角度テーブルに要素ごとの回転角度θを書き込む。
【００７６】
次に、プロセスＰ２６で、回転軸ベクトルＶＳＥ１、ＶＳＥ２…及び回転角度θのデータに基いて全要素をｘｙ平面に平行になるように回転させ、要素分離モデルを生成する。このとき、図３８に示すような解析モデルを構成する各要素Ｅ１、Ｅ２…は、図３９に示すように、分離されて分離要素ＥＳ１、ＥＳ２…として上記要素分離モデルを構成することになる。
【００７７】
一方、プロセスＰ２７で、要素分離モデルを構成する各分離要素を構成する節点の番号ＮＳ１、ＮＳ２…を図１５の分離要素構成テーブルに書き込む。
【００７８】
次に、プロセスＰ２８で、図７の要素構成テーブルと図１５の分離要素構成テーブルとから、要素回転前の節点Ｎ１、Ｎ２…（元節点）及び分離要素の節点ＮＳ１、ＮＳ２…を読み出して、元節点Ｎ１、Ｎ２…を共有していた分離要素の各節点番号ＮＳ１、ＮＳ２…を図１６の元節点−分離要素の節点テーブルに書き込む。
【００７９】
また、プロセスＰ２９で、各分離要素の節点のｘ座標、ｙ座標、ｚ座標を図１７の分離要素の節点座標テーブルに書き込む。
【００８０】
そして、以下に説明する平面展開モデル作成においては、これらの各分離要素を図４０に示すようにｘｙ平面に投影させ、つまり各分離要素のｚ座標を無視して各分離要素を同一のｘｙ平面上に並べた状態で処理を行うことになる。
【００８１】
プロセスＰ３０では、拘束条件の設定を行う。まず、解析モデルで特徴的な位置にある２つの節点ＮＡ、ＮＢを選択し、これらを図２０の拘束条件テーブルに書き込む。そして、これらの節点ＮＡ、ＮＢは、要素回転後には、各分離要素を構成する複数の節点ＮＳＡ１、ＮＳＡ２…に分割されることになる。例えば、図４１（ａ）に示すように、要素ＥＡ１〜ＥＡ４で囲まれた節点ＮＡは、要素回転後にｘｙ平面に投影すると図４１（ｂ）に示すように分離要素ＥＳＡ１〜ＥＳＡ４の節点ＮＳＡ１〜ＮＳＡ４として分離する。そして、これらの節点ＮＳＡ１〜ＮＳＡ４のｘ座標及びｙ座標を図１７の分離要素の節点座標テーブルから読み出し、これらの節点ＮＳＡ１〜ＮＳＡ４の座標の平均座標（ＮＳＡＶｘ、ＮＳＡＶｙ）を求める。次に、この平均座標の点を平均点ＮＳＡＶとすると、図４１（ｃ）に示すように、節点ＮＳＡ１〜ＮＳＡ４が平均点ＮＳＡＶに一致するように、各分離要素ＥＳＡ１〜ＥＳＡ４をｘｙ平面に平行移動させる。そして、節点ＮＢについても同様にして、図１７の分離要素の節点座標テーブルにおいて、節点ＮＳＡ１〜ＮＳＡ４、ＮＳＢ１〜ＮＳＢ４を有する分離要素を構成する各節点の座標を、前述の分離要素を平行移動させた後の各節点の座標に書き換える。
【００８２】
以上により、解析用基礎データ作成サブプログラムによる解析用基礎データの作成処理が終了し、プロセスＰ３１で、得られたデータを解析用基礎データファイルに格納する。そして、コンピュータ１０の動作は、図３１に示すメインプログラムによる処理に戻り、プロセスＰ４以下を実行する。
【００８３】
メインプログラムのプロセスＰ４では、今回の解析対象の製品に関するデータ、即ち当該材料のヤング率、ポアソン比、降伏応力、接線係数等の解析に必要なデータを図２１の材料特性テーブルから読み取る。
【００８４】
そして、プロセスＰ５で上記基礎データに基いて解析用データの作成を行う。この解析用データは、上記プログラム記録部１４ａに記録されている解析用データ作成サブプログラムにより作成される。そして、このサブプログラムの動作を図４２に示すフローチャートに従って説明すると、まず、プロセスＰ４１で、解析計算のために必要となる各要素ごとの剛性マトリクスＫｅの作成を行う。例えば４節点を有する分離要素の場合、分離要素ごとの自由度は、第１〜第４節点のｘ変位、ｙ変位の合計８自由度となり、要素剛性マトリクスＫｅは数式１のように表される。なお、行列に含まれる成分ｋ１１〜ｋ８８は、分離要素の材料に基づくヤング率に応じた値である。
【００８５】
【数１】
2006185228.tif 000003
【００８６】
ここで、各節点の未知の変位ベクトルをＵｅとすると、各要素に作用する荷重ベクトルＦｅとの間には、数式２に示す関係式が成り立つ。
【００８７】
【数２】
2006185228.tif 000004
【００８８】
但し、このとき荷重ベクトルＦｅはゼロベクトルである。
【００８９】
次に、プロセスＰ４２で、解析計算のために必要となるバネ要素ごとの剛性マトリクスＫｓの作成を行う。例えば２節点を結びつけるバネ要素の剛性マトリクスＫｓは、自由度が第１節点、第２節点のｘ変位、ｙ変位であり、合計で４自由度となって、数式３に示すように表される。
【００９０】
【数３】
2006185228.tif 000005
【００９１】
ここで、Ｅは分離要素の材料に応じたヤング率、αは十分大きな定数である。定数αは、分離要素を適正に変形させる目的でバネ要素の剛性マトリクスＫｓに乗算されるパラメータであって、例えば材料が鉄の場合、鉄のヤング率Ｅは２１０００ｋｇｆ／ｍｍ２であり、この数値のみで上記バネ要素の剛性マトリクスＫｓを作成して解析計算を行った際に、分離要素に対してバネ要素が軟らかいことになって、分離要素に十分な変形が生じていない状態で釣り合うという事態が発生しうる。そのため、ヤング率Ｅより約６桁上のオーダーの定数αを乗算して、バネ要素を十分に硬いものとして、対応する節点が一致するまで分離要素を十分に変形させるのである。なお、定数αが過度に大きいときには、解析精度が低下するおそれがあり、定数αは前述のようにヤング率Ｅの約６桁上のオーダーが妥当である。また、各分離要素が同一の材料であるときは、各ヤング率に乗算する定数αは同一にすればよいが、材料の異なる分離要素間にバネ要素を設定した場合は、分離要素ごとに材料に応じて定数αが異なるようにしてもよい。
【００９２】
一方、上記プロセスＰ４１の要素剛性マトリクスＫｅ作成時と同様に、一つのバネ要素について剛性マトリクスＫｓ、変位ベクトルＵｓ、要素に作用する荷重ベクトルＦｓの間には数式４に示す関係式が成り立つ。
【００９３】
【数４】
2006185228.tif 000006
【００９４】
次に、プロセスＰ４３で、バネ要素間に作用する荷重データを作成する。つまり、荷重データは、バネ要素が設定された節点間のｘ方向の距離ｄｘ、ｙ方向の距離ｄｙを成分とした荷重ベクトルＦｓであって数式５に示すように表される。
【００９５】
【数５】
2006185228.tif 000007
【００９６】
ここで、Ｅは分離要素の材料に応じたヤング率、αは上記数式３で示した十分大きな定数である。
【００９７】
さらに、前述の拘束条件に係る分離要素の節点については、これらの節点間にバネ要素を設定してバネ要素の剛性マトリクスＫｓを作成すると共に、荷重ベクトルｆｓにおける節点間距離ｄｘ、ｄｙはゼロになる。
【００９８】
以上のようにして要素剛性マトリクスＫｅ、バネ要素の剛性マトリクスＫｓ、バネ要素間に作用する荷重ベクトルＦｓを解析用データとし、図３１のフローチャートのプロセスＰ６に進んで解析計算を行う。この解析計算は、解析計算サブプログラムにより行われる。
【００９９】
ここで、この解析計算について詳述すると、まず、上記のようにして算出した要素剛性マトリクスＫｅ及びバネ要素の剛性マトリクスＫｓを全体剛性マトリクスＫｇへ足し込みを行うと共に、要素の荷重ベクトルＦｅ及びバネ要素の荷重ベクトルＦｓを全体荷重ベクトルＦｇへ足し込みを行う。
【０１００】
なお、以下の説明では、図４３に示すような４節点を有する２つの分離要素ＥＳ１、ＥＳ２において、節点ＮＳ２と節点ＮＳ７との間にバネ要素Ｓを設定した場合を例に説明する。即ち、図４３に示したように節点ＮＳ１〜ＮＳ４を有する分離要素ＥＳ１と節点ＮＳ５〜ＮＳ８を有する分離要素ＥＳ２との２つの分離要素を全体剛性マトリクスＫｇに足し込みを行う場合、分離要素ＥＳ１は、（×）で示す節点ＮＳ１のｘ成分、ＮＳ１のｙ成分、ＮＳ２のｘ成分、ＮＳ２のｙ成分、…、ＮＳ４のｘ成分、ＮＳ４のｙ成分の順に１〜８行を構成し、分離要素ＥＳ２は、（＋）で示すＮＳ５のｘ成分、ＮＳ５のｙ成分、ＮＳ６のｘ成分、ＮＳ６のｙ成分、…、ＮＳ８のｘ成分、ＮＳ８のｙ成分の順に９〜１６行を構成して全体剛性マトリクスＫｇは数式６に示すような１６行のマトリクスとなる。なお、全体剛性マトリクスＫｇの行列の数は総節点数の２倍となる。
【０１０１】
【数６】
2006185228.tif 000008
【０１０２】
次に、上記全体剛性マトリクスＫｇにバネ要素の足し込みを行う。この例においては、分離要素ＥＳ１の節点ＮＳ２と分離要素ＥＳ２の節点ＮＳ７とを結ぶバネ要素Ｓが設定されているので、上記バネ要素の剛性マトリクスＫｓの節点ＮＳ２、ＮＳ７に対応する成分を、全体剛性マトリクスＫｇにおける節点ＮＳ２、ＮＳ７に対応する位置（数式７で○で囲まれた箇所）にそれぞれ足し込みを行う。その結果、数式７に示すような全体剛性マトリクスＫｇが得られる。
【０１０３】
【数７】
2006185228.tif 000009
【０１０４】
続いて、上記バネ要素の荷重ベクトルＦｓの全体荷重ベクトルＦｇへの足し込みを行う。このときも同様にバネ要素の荷重ベクトルＦｓの節点ＮＳ２、ＮＳ７に応じた成分を、全体荷重ベクトルＦｇにおける節点ＮＳ２、ＮＳ７に対応する位置に足し込みを行い、数式８に示すような荷重ベクトルＦｓが得られる。
【０１０５】
【数８】
2006185228.tif 000010
【０１０６】
なお、この例は２節点間のバネ要素について示しているが、３節点以上の場合には総当りでバネ要素を使用する。つまり、３節点間にそれぞれバネ要素が設定された場合は、第１の節点−第２の節点、第１の節点−第３の節点、第２の節点−第３間の節点の３通り、４節点の場合は、第１の節点−第２の節点、第１の節点−第３の節点、第１の節点−第４の節点、第２の節点−第３の節点、第２の節点−第４の節点、第３の節点−第４の節点間の６通り、Ｎ節点の場合は、Ｎ×（Ｎ−１）／２通りのバネ要素が設定されることになる。
【０１０７】
以上の手順による足し込みを全ての分離要素及びバネ要素について行うと、全体剛性マトリクスＫｇ及び全体荷重ベクトルＦｇはそれぞれ数式９に示す式で表される。なお、ｎｅは分離要素の総数、ｎｓはバネ要素の総数を表す。
【０１０８】
【数９】
2006185228.tif 000011
【０１０９】
そして、全体剛性マトリクスＫｇの逆行列を求めて数式１０に示すように変形し、これを解くことによって変位ベクトルＵｇが求まる。
【０１１０】
【数１０】
2006185228.tif 000012
【０１１１】
このとき、節点同士が一致するように各要素が変形し、全体として荷重の釣合を保持しながら、解析モデルを平面形状に展開したことになる。そして、メインプログラムのプロセスＰ７で、展開計算の結果の各節点のｘ座標及びｙ座標を図２２の解析結果テーブルに書き込む。ここで、得られた平面形状を、図４４の画面１５ｅに平面展開モデルとして示す。なお、第１中間形状の平面展開モデルは、該形状を成形するために最適となるブランク形状である。
【０１１２】
次に、プロセスＰ８で上記展開計算の結果に基いて、つまり各分離要素の変形に基いて応力、歪、厚みの算出を行い、図２３の歪テーブル及び図２６の応力テーブルに書き込む。
【０１１３】
また、上記展開計算においては、要素のｘｙ平面における変形に基いて応力、歪（面内歪）、厚みを単に算出することができるが、これに加えて曲げ歪による影響を考慮することによって、より精度の高いシミュレーションを行うと共に、スプリングバックの発生を事前に予測することができるようになる。
【０１１４】
まず、各レイヤーの面内歪について説明すると、図２５に示したように、７層のレイヤーＬ１〜Ｌ７で構成された要素Ｅは、各レイヤーＬ１〜Ｌ７ごとに４つの評価点ＤＬ１ａ〜ＤＬ１ｄ、…ＤＬ７ａ〜ＤＬ７ｄを有している。そして、解析モデル展開時に生じた面内歪により、各評価点ＤＬ１ａ〜ＤＬ１ｄ、…ＤＬ７ａ〜ＤＬ７ｄは各レイヤーＬ１〜Ｌ７の伸び縮みによって面内方向にずれる。このときのずれが各レイヤーＬ１〜Ｌ７の面内歪に相当し、例えば図４５に示すように、評価点ＤＬ１ａ〜ＤＬ７ａにおいては同様の面内歪εｘ′となる。そして、前述の図２３の歪テーブルは、全ての評価点におけるｘ方向の面内歪εｘ′、ｙ方向の面内歪εｙ′、及びｚ方向の面内歪εｚ′、また、各レイヤーの剪断方向の歪について、ｘ方向にみたｙ方向の面内剪断歪γｘｙ′、同様に面外剪断歪γｙｚ′、γｚｘ′を記録している。
【０１１５】
次に、各評価点における曲げ歪の算出方法について説明すると、曲げ歪の計算においては、まず、図１８の節点ベクトルテーブルを作成し、特定の要素について要素ベクトルと、周囲の４つの節点ベクトルとの差を、要素に発生した全体としての曲げ変形とみなすようにしている。
【０１１６】
例えば図４６に示すように、節点Ｎｗ１、Ｎｗ２、Ｎｗ３、Ｎｗ４で構成された要素Ｅｗについては、各節点ベクトルＶＮｗ１、ＶＮｗ２、ＶＮｗ３、ＶＮｗ４と、要素ベクトルＶＥｗとのベクトル差ΔＶＮｗ１〜ΔＶＮｗ４は、それぞれ数式１１のように表される。なお、要素Ｅｗは、レイヤーＬ１〜Ｌ７で構成されている。
【０１１７】
【数１１】
2006185228.tif 000013
【０１１８】
ここで、ベクトル差ΔＶＮｗ１〜ΔＶＮｗ４をｘ成分、ｙ成分、ｚ成分に分解し、各節点ベクトルＶＮ１、ＶＮ２…の要素ベクトルＶＥ１、ＶＥ２…に対する角度差を算出する。そして、図４７に示すように、要素Ｅｗについて、その要素のなす面に最も近い平面にＸＹ座標系を設定することにより、図４８に示すように、上記ｘｙｚ座標系における３方向で表された角度差をＸ軸回りの角度差ΔθＸとＹ軸回りの角度差ΔθＹとの２方向で表すことができる。一方、曲率は角度差の変化率であるから、以下に示す数式１２により求められる。なお、ΔθＸのＸ軸方向の変化率を曲率φＸＸ、ΔθＸのＹ軸方向の変化率を曲率φＸＹ、ΔθＹのＸ軸方向の変化率を曲率φＹＸ、ΔθＹのＹ軸方向の変化率を曲率φＹＹとしている。
【０１１９】
【数１２】
2006185228.tif 000014
【０１２０】
そして、この数式により４つの評価点における曲率φＸＸ、φＸＹ、φＹＸ、φＹＹをそれぞれ求める。
【０１２１】
次に、要素Ｅｗの板厚中心からの各評価点までの厚さ方向の距離ＺＬを求める。ここでは板厚を７等分して各レイヤーＬ１〜Ｌ７を構成しているので、レイヤーＬ４における板厚中心から評価点までの厚さ方向の距離ＺＬ４はゼロ、内側のレイヤーＬ１〜Ｌ３の評価点までの距離ＺＬ１〜ＺＬ３は負の値、外側のレイヤーＬ５〜Ｌ７までの距離ＺＬ５〜ＺＬ７は正の値になる。
【０１２２】
そして、ｘ軸方向の曲げ歪εｘ″、ｙ軸方向の曲げ歪εｙ″、及びｘ軸からみたｙ軸方向の曲げ剪断歪γｘｙ″は、数式１３により求められる。
【０１２３】
【数１３】
2006185228.tif 000015
【０１２４】
図４９に示すように、例えば評価点ＤＬ５ａのｘ方向の曲げ歪εｘ″は、曲率φＸＸに距離ＺＬ５を乗算して求めた円弧長として求められる。
【０１２５】
一方、これらの曲げ歪εｘ″、εｙ″、γｘｙ″と図２３の歪テーブルに記録された面内歪εｘ′、εｙ′、γｘｙ′とを足し合わせることによりによりトータルの歪εｘ、εｙ、γｘｙが算出され、歪テーブルに評価点ごとに記録する。数式１４にｘ軸方向のトータルの歪εｘの算出過程を示す。
【０１２６】
【数１４】
2006185228.tif 000016
【０１２７】
さらに、数式εｘと同様の方法で、ｙ軸方向のトータルの歪εｙ、及びｘ軸方向からみたｙ軸方向のトータルの剪断歪γｘｙを数式１５に示す。
【０１２８】
【数１５】
2006185228.tif 000017
【０１２９】
なお、図５０に示すように、曲げ歪は、曲げの内側（Ｌ１側）が圧縮（負）で外側（Ｌ７側）が引っ張り（正）となり中央でゼロとなる直線的な分布となる。そのため、例えば面内歪が図４５に示したような一様な引張特性を有する場合には、トータルの歪の分布は、図５１に示すように、上記曲げ歪の分布が引張側にずれることになる。
【０１３０】
一方、図２３の歪テーブルに記録されたトータルの歪を用いて、応力―歪特性により、各評価点における応力σｘｘ、σｙｙ、σｚｚ、剪断応力τｘｙ、τｙｚ、τｚｘの値が計算される。そして、このように計算された応力を、評価点ごとに図２６に示す応力テーブルに書き込む。
【０１３１】
次に、プロセスＰ９で、上記プロセスＰ１で読み込んだ全ての設計形状についての展開計算が終了したか否かを判定し、終了していない場合はプロセスＰ２に戻り、終了した場合はプロセスＰ１０に進む。例えば上記説明では第１中間形状についての解析のみが終了しているので、第２中間形状及び最終形状についての解析を行うためにプロセスＰ２〜Ｐ８を繰り返すことになる。そして、このプロセスを繰り返すことで、各設計形状ごとにデータファイル記録部１４ｂ…１４ｂが作成されることになる。
【０１３２】
そして、プロセスＰ１０〜Ｐ１５で、複数の成形経路が存在する場合は、各成形経路についての安全性を評価し、また、各中間形状におけるトリミングラインの設定を行う。つまり、上記プロセスＰ２〜Ｐ８の解析においては、各工程後の中間形状の設計案が図３２の第１中間形状と図３３の第２中間形状と１つずつ示したが、ここでは工程後ごとに複数の中間形状の設計案を用意して、これらの中から最適な中間形状を選び出すと共に最適成形経路を見出すことを目的としている。例えば図５３に示すように、第１中間形状として２つの設計案、第２中間形状として３つの設計案、第３中間形状として３つの設計案がそれぞれ用意されている場合を考えると、初期形状から最終形状に至る成形経路としては１８通りあり、この中から最適成形経路を見出すことになる。
【０１３３】
まず、プロセスＰ１０として１８通りの成形経路の中から任意の成形経路を選び出す。そして、プロセスＰ１１で、当該成形経路について、上記プログラム記録部１４ａに記録されている評価用データ作成サブプログラムにより評価用データを作成する。
【０１３４】
ここで、この評価用データ作成サブプログラムの動作を図５４に示すフローチャートに従って説明する。
【０１３５】
まず、プロセスＰ５１で工程ｍ＝２を設定する。そして、プロセスＰ５２で工程ｍ前後の設計形状による平面展開モデルの位置合せを行う。例えば、図４４に示した第１中間形状の平面展開モデルと図５５の画面１５ｆに示すような第２中間形状の平面展開モデルとを、拘束点ＮＡ、ＮＢ同士を一致するように重ね合わせると、図５６の画面１５ｇに示すように両者の位置関係が決定されることになる。
【０１３６】
そして、プロセスＰ５３で、例えば図５７に示すように、第２中間形状（工程ｍ後）の平面展開モデルの評価点Ｄｍの座標（ｘ，ｙ）から、図２２の解析結果テーブルの座標データに基いて、評価点Ｄｍを含む第１中間形状（工程ｍ前）の平面展開モデルにおける要素の番号Ｅｍ−１を読み出し、同様に第２中間形状の全ての評価点についてこれを含む要素番号Ｅ１、Ｅ２…を読み出す。
【０１３７】
次に、プロセスＰ５４で、図５８に示すように要素Eｍ−１にｒｓ座標系を設定し、プロセスＰ５５で、ｒｓ座標系における評価点Ｄｍの座標（ｒｍ，ｓｍ）を算出する。
【０１３８】
そして、プロセス５６に進んで、工程ｍ前後に係る評価点マッピングテーブルを作成する。即ち、工程後の評価点ごとに、上記プロセスＰ５３で求めた当該評価点を含む工程前の要素番号、及び上記プロセスＰ５４で求めた該要素における評価点の座標のｒ値、ｓ値を図２７の評価点マッピングテーブルにそれぞれ書き込む。
【０１３９】
次に、プロセスＰ５７に進み、図２８に示す工程ｍ前の補間歪テーブルを作成する。この補間歪テーブルは、工程ｍ後の評価点を含む工程前の要素の座標（ｒ，ｓ）における歪εｘ、εｙ、εｚ、剪断歪εｘｙ、εｙｚ、εｚｘ、及び板厚減少率を記録している。このとき、前述の展開計算により各要素は変形し、各要素を構成する各節点は変位したが、要素を構成する各節点の変位量及び変位方向に基いて、当該要素各部の歪分布をｒｓ座標系で表された近似関数で定義することができ、上記評価点マッピングテーブルの座標データに基づいて、工程ｍ前の平面展開モデルにおける各評価点に対応する座標（ｒ，ｓ）の歪を補間して求めることができる。
【０１４０】
また、プロセスＰ５８で、図２７に示す工程ｍ前の補間応力テーブルを作成する。この補間応力テーブルは、各評価点ごとに座標（ｒ，ｓ）の応力σｘｘ、σｙｙ、σｚｚ、剪断応力σｘｙ、σｙｚ、σｚｘを記録している。
【０１４１】
次に、プロセスＰ５９で、評価用応力テーブルを作成する。この評価用応力テーブルは、工程１後は図２６の応力テーブルに記録された各応力と同様であり、工程２後は特別な計算を行う。この工程２後の応力の算出方法について図５９を用いて説明する。ここでは、図５３で示した各中間形状の設計案の組合せにおいて、第１中間形状を第２設計案、第２中間形状を第３設計案、第３中間形状を第１設計案とした成形経路を選択した場合を例にとる。
【０１４２】
まず、工程１後の第２設計案に係る評価用応力テーブルは、展開計算のプロセスＰ８で作成したものであり、ここで作成した各評価点の歪εＩ及び応力σＩは、後述の工程の安全性評価に採用する。
【０１４３】
ところで、応力と歪の関係が線形の場合は、数式１６に示すように、歪の増分Δεから応力の増分Δσが求められ、これを工程前の対応する評価点の応力σｎ−１に足し合わせることによって工程後の応力σｎ求めることができる。
【０１４４】
【数１６】
2006185228.tif 000018
【０１４５】
しかしながら、工程２以降の応力の計算においては、工程前における応力の分布に依存して応力と歪の関係が非線形となり、数式１２によっては正確な応力は得られない。即ち、例えば図６０の応力−歪の関係のグラフで示すように、板材を曲げたときの特性は、歪に対して応力が上に凸のカーブを描きながら点Ａから点Ｂに増加し、ここから曲げ戻しを行った場合、少しの歪の減少に対して応力は点Ｂから点Ｃに急激に減少するような特性を示し、ここから再び曲げたときには点Ｃから点Ｄに応力が急激に増加すると共に、点Ｄから点Ｅ間では緩やかに応力が増加するような特性を示す。このように、多段階成形においては、曲げ戻しも考慮しなければならず、上記数式１２では曲げ戻しの影響を加味することができない。なお、工程１後の応力については、このような非線形性を考慮する必要がないため、数式１２の式で算出した応力を安全性評価に採用することができるのである。また、歪は曲げ戻しの有無に拘らず各要素の変位に対して線形の関係にあるので、工程によらず図２３の歪テーブルに記録されたものが安全性評価に採用できる。
【０１４６】
そのため、図２６の応力テーブルにおいて、図５９の工程２後の第３設計案に係る各評価点の応力σＩＩは、上記応力と歪が線形の関係にあるという前提で算出されたものであるため正確な値でなく、以下に説明する応力σＩＩ′として別途計算し直して図３０の評価用応力テーブルに記録される。即ち、以下に示す数式１７により応力を計算する。
【０１４７】
【数１７】
2006185228.tif 000019
【０１４８】
この式においてＤはεの関数である。ここで、より具体的に説明すると、数式１８にようにＤの関数を積分することにより応力の増分Δσを算出する。
【０１４９】
【数１８】
2006185228.tif 000020
【０１５０】
そして、数式１９は、数式１８の積分表記の応力の増分Δσをコンピュータ計算用に工程ｎ−１から工程ｎまでの区間をｍ等分して足し算の和として表し、これにより計算された応力の増分Δσを工程前の応力σｎ−１に加算することになる。
【０１５１】
【数１９】
2006185228.tif 000021
【０１５２】
ここで、ｍは１０程度に設定して計算量が多くならないようにしている。
【０１５３】
一方、上記数式１７において、斜体で示した項は図２８の補間歪テーブルに記録された補間歪及び図２９の補間応力テーブルに記憶された補間応力である。このとき、図２３の歪テーブルに記録された評価点における歪から、図２７の評価点マッピングテーブルにおいて当該評価点に対応する座標（ｒ，ｓ）における補間歪を減算して当該評価点における歪の増分を算出し、この増分に基いて応力の増分を算出する。次に、工程２前の第２設計案において補間した応力に、上記応力の増分を加算することによって、工程２後の第３設計案における応力σＩＩ′が求められる。そして、求めた応力ＩＩ′を評価用応力テーブルに記録する。
【０１５４】
図３０の評価用応力テーブルの、工程３後の第１設計案に係る各評価点の応力σＩＩＩ′は、数式２０により算出される。なお、この数式におけるＤは、上記数式１７で使用したものと同じ関数である。
【０１５５】
【数２０】
2006185228.tif 000022
【０１５６】
工程３以降においては、応力の増分を加算する対象として、工程３後の第１設計案の各評価点に対応する点で、評価用応力テーブルに記録された工程２後の第３設計案に係る応力分布に基いて補間した応力σＩＩ′を用いることになる。
【０１５７】
そして、プロセスＰ６０で、工程ｍが最終工程４か否かを判定し、工程ｍが最終工程４でないときは、プロセスＰ６１でｍに１を加算して次の工程を設定し、プロセスＰ５２に戻る。
【０１５８】
そして、プロセスＰ６０で工程ｍが最終工程４である場合は、図３１のメインプログラムに戻ってプロセスＰ１２で、評価指標値算出、即ち当該成形経路の評価を行う。
【０１５９】
この評価指標値算出は、成形経路評価サブプログラムにより実行され、当該成形経路における加工の際の安全性、即ちしわや割れ等の欠陥の発生などを評価するものであって、各設計案及び最終形状における各評価点の応力、歪、板厚減少率に基いて、各成形経路ごとに評価即ちしわや割れ等の発生の指標となる評価指標値を求めることになる。なお、応力、歪、板厚減少率のデータを用いて安全性の評価を行う方法は各種のものがあるが、以下にその一例を説明する。
【０１６０】
まず、図２３の歪テーブルに記録された主歪ε１、ε２を用いた評価では、縦軸を主歪ε１、横軸をε２としたグラフに各評価点の主歪ε１、ε２をプロットして、成形限界線に対する位置関係を調べる。成形限界線は、材料に応じて決まる曲線であり、予め実験的又は理論的に導き出されたものである。そして、図６１に示すように、プロットした点が成形限界線より上方に位置する場合（●で示す）は、成形の結果欠陥が生じる可能性が高く、成形限界線より下方に位置する場合（○で示す）は、成形の結果欠陥が生じる可能性が低いと判断できる。そして、この評価を当該成形経路を構成する各設計形状について行い、評価結果に基いて当該成形経路の評価指標値を求める。また、成形限界線から大きく下方に逸脱する点が存在する場合は評価指標値をゼロに設定する。なお、主歪ε１（最大歪）、ε２（最小歪）は、数式２１により算出されたものである。
【０１６１】
【数２１】
2006185228.tif 000023
【０１６２】
また、図２３の歪テーブルに記録された板厚減少率は、成形前の評価点における板厚が成形後に減少した割合を示し、板厚減少率が８０％以上になったときにしわやわれが発生する可能性が高いなどの評価を行うことができる。そして、当該成形経路を構成する各設計形状について板厚減少率を求め、これに基いて当該成形経路の評価指標値を求める。なお、板厚減少率が８０％以上の箇所が存在する場合は評価指標値をゼロに設定する。
【０１６３】
一方、成形により一部塑性変形が生じているときには、図３０の評価用応力テーブルに記録された評価点の厚み方向の応力分布に応じて、外部からの拘束が外されたときに、板厚内部では弾性変形部分が復元しようとするのに対して、塑性変形部分がこの復元を妨げるように作用する。その結果、板厚方向の応力分布に応じた所定のモーメントが生じ、スプリングバックが発生する。そして、スプリングバックが生じた後は、図６２に示すように、厚み方向の残留応力分布が釣り合った状態で安定することになる。なお、スプリングバック量の算出に際しては、厚さ方向の応力分布から各節点に作用するモーメントを計算し、該モーメントによりたわみ量を計算することができ、該たわみ量により残留応力分布を得る方法が周知である。そして、求めたスプリングバック量により評価指標値を求める。
【０１６４】
そして、得られた各成形経路ごとの評価指標値を一旦記憶し、図３１のメインプログラムに戻ってプロセスＰ１３に進む。ここでは全成形経路、即ち図５３の例においては１８通り全ての成形経路についての評価が完了したか否かを判定し、完了していない場合はプロセスＰ１０に戻って未評価の成形経路について評価を行うことになる。そして、全成形経路について評価が完了している場合は、プロセスＰ１４に進み、記憶された成形経路の評価結果に基いて最適成形経路を決定する。具体的には、各成形経路の評価指標値を比較して、所定の閾値を越えているものから、最も数値の高いものを最適成形経路とする。また、各種の評価手法により得られた評価指標値を複合的に評価して最適成形経路を求めるようにしてもよい。例えば主歪、板厚減少率により求めた評価指標値のうち、所定の閾値を超えているものに対して、スプリングバック量による評価指標値が最も良いものを最適成形経路にする等である。
【０１６５】
なお、ここでは各成形経路について初期形状から最終形状に至る評価結果を全て揃えた上で安全性評価を行うようにしているが、各工程ごとに評価を行い、評価指標値により欠陥の発生が判明した時点でその成形経路については評価をＮＧとして、次の成形経路についての評価を開始するようにしてもよく、さらに、欠陥の発生が判明した設計案を含む成形経路すべてについて以後シミュレーションを中止するようにしてもよく、これらによって計算量の削減を図ることができる。
【０１６６】
さらに、プロセスＰ１５では、トリミングラインの決定を行う。図３２に示した第１中間形状に対してトリミングを施して図３３に示した第２中間形状に成形するときを例に説明すると、評価点マッピングテーブルで第１中間形状の平面展開モデルにおいて第２中間形状の平面展開モデルの各評価点の座標を含む要素を、第１中間形状に反映させる。つまり、図６３の表示装置１５ｈの画面に斜線の領域で示すように反映され、斜線の領域の輪郭がトリミングラインとなる。これによって、金型の設計者は、各工程の中間設計形状を示すだけで、自動的にトリミングラインを設定することができる。
【０１６７】
以上のような方法によれば、多段階成形において、各設計形状についての平面展開モデルを同様に繰り返し作成し、各平面展開モデルの各要素の変形に基いて応力、歪、又は厚みの少なくとも１つを求め、各要素の応力、歪、又は厚みに基いて各工程によるしわや割れ等の欠陥の発生などを評価することができる。このように、各中間形状をあたかも最終形状のように扱うことによって、欠陥の発生する可能性の高い工程を特定することができると共に、得られた応力、歪、又は厚みの情報（例えば、工程前後の歪の差）に基いて中間形状の設計変更等を行うなどの対策をとることができる。
【０１６８】
また、各設計形状について作成された解析モデルからワンステップで平面展開モデルを作成し、展開計算結果の各要素の変形に基づいて計算された各要素の応力、歪、又は厚みのデータにより、多段階成形のシミュレーションに使用する大方のデータが揃うことになる。従って、従来のように金型や周囲の部材等をモデル化する必要がないので扱うデータ量が少なく、各平面展開モデルの各要素についての応力、歪、厚みの値を組み合わせるだけで簡単に各工程での欠陥の発生を評価することができ、圧倒的な計算の簡素化及び高速化を実現することができる。
【０１６９】
さらに、金型の設計者が各工程で得られる中間形状の複数の設計案を用意する場合、各設計案における各要素の応力、歪、又は厚みに基いて、最適な成形経路、即ち最適となる中間設計形状の組合せを判定することができる。
【０１７０】
このように、複数の成形経路から最適なものを選択する際も、各設計形状ごとにワンステップ法により得られた応力、歪、又は厚みのデータを組合せるだけでよく、扱うデータ量が少なく、計算の簡素化及び高速化が実現される一方、多くの中間設計形状の設計案を用意してより適当となる中間設計形状を得ることができる。
【０１７１】
また、各要素の厚み方向に複数のレイヤー要素を設定して各レイヤー要素において、展開計算で得られた歪と各設計形状の曲率に基いて計算された曲げ歪とを加算して合計歪とし、工程前後の平面展開モデルの対応するレイヤー要素の合計歪の差を算出して、この差に基いて各レイヤー要素の応力が計算できると共に、各要素の厚み方向の曲げ歪を考慮した応力分布が得られ、スプリングバックを適正に評価することができる。そして、金型の設計者は、厚み方向の応力分布に基いて金型を設計することができ、例えば加工の際にスプリングバック量を見込んで余分に曲げるように金型を設計する等の対策をとることができ、工程により得られる形状の寸法精度を向上させることができる。このように、ワンステップ法により得られたデータを組み合わせてスプリングバックの評価を行うことができるので、従来のシステムに比べて計算の簡素化及び高速化を実現することができる。
【０１７２】
ところで、プログラム記録部１４ａに記録されたマッピング処理サブプログラム及び描画サブプログラムは、応力、歪、板厚の分布等を表示装置１５に出力するためのプログラムであって、マッピング処理サブプログラムにより平面展開モデル又は製品形状に各分布を反映させ、描画サブプログラムにより分布図として表示することになる。
【０１７３】
なお、上記解析モデル作成プログラム、解析計算プログラム、描画サブプログラム、及びマッピング処理サブプログラムは、コンピュータ１０に予めインストールされた既存のプログラムを用いたものであってもよい。この場合、既存の解析モデル作成サブプログラムにより作成された解析モデルは解析用基礎データ作成サブプログラム等のフォーマットに統一され、他のプログラムとの互換性が確保される。一方、解析用データ作成サブプログラムにより作成した解析用データは、既存の解析計算サブプログラムのフォーマットに統一され、解析用データは解析計算サブプログラムとの互換性が確保される。同様に、既存の解析計算サブプログラムの解析結果は、描画サブプログラム及び加工用データ作成サブプログラムのフォーマットに統一され、これらのプログラムとの互換性が確保される。また、評価用データ作成サブプログラムにより作成した評価用データは、成形経路評価サブプログラムのフォーマットに統一され、評価用データは成形経路評価サブプログラムのフォーマットに統一される。
【産業上の利用可能性】
【０１７４】
本発明は、上記のような実情に鑑み、多段階成形のシミュレーションの簡素化及び高速化を図ると共に、スプリングバックのシミュレーションをも同時に行うことができるシステムを提供する。本発明は、板状ブランクを用いた成形の結果を有限要素法を用いてシミュレーションするシステム、より詳しくは、複数工程による各中間形状を経て最終製品形状に成形する多段階成形をシミュレーションする多段階成形シミュレーションシステム及び多段階成形シミュレーション用プログラムに関し、コンピュータによる有限要素解析技術の分野に広く好適である。
【図面の簡単な説明】
【０１７５】
【図１】本発明の概念の説明に係る製品形状及び解析モデルの説明図である。
【図２】同じくトリミングラインの説明図である。
【図３】本発明の実施の形態に係るシステムの全体構成を示すブロック図である。
【図４】同システムの記録装置に記録されているプログラムの説明図である。
【図５】同じくデータファイルの説明図である。
【図６】解析モデルの一部を拡大した各要素及び節点の説明図である。
【図７】解析モデルデータファイルに含まれる要素構成テーブルの説明図である。
【図８】同じく節点座標テーブルの説明図である。
【図９】解析用基礎データファイルに含まれる節点構成テーブルの説明図である。
【図１０】同じく要素中心座標テーブルの説明図である。
【図１１】同じく要素ベクトルテーブルの説明図である。
【図１２】同じく要素回転軸ベクトルテーブルの説明図である。
【図１３】要素回転軸ベクトル生成方法の説明図である。
【図１４】解析用基礎データファイルに含まれる要素回転角度テーブルの説明図である。
【図１５】同じく分離要素構成テーブルの説明図である。
【図１６】同じく元節点−分離要素の節点テーブルの説明図である。
【図１７】同じく分離要素の節点座標テーブルの説明図である。
【図１８】同じく節点ベクトルテーブルの説明図である。
【図１９】節点ベクトルテーブル算出過程の説明図である。
【図２０】同じく拘束条件テーブルの説明図である。
【図２１】材料属性データファイルを構成する材料特性テーブルの説明図である。
【図２２】解析結果データファイルに含まれる解析結果テーブルの説明図  である。
【図２３】同じく歪テーブルの説明図である。
【図２４】評価点の説明図である。
【図２５】レイヤーの説明図である。
【図２６】解析結果データファイルに含まれる応力テーブルの説明図である。
【図２７】評価用データファイルに含まれる評価点マッピングテーブルの説明図である。
【図２８】同じく補間歪テーブルの説明図である。
【図２９】同じく補間応力テーブルの説明図である。
【図３０】同じく評価用応力テーブルの説明図である。
【図３１】システム全体を制御するメインプログラムの動作を示すフローチャートである。
【図３２】第１中間形状のＣＡＤ図面を表示する画面の説明図である。
【図３３】第２中間形状のＣＡＤ図面を表示する画面の説明図である。
【図３４】最終形状のＣＡＤ図面を表示する画面の説明図である。
【図３５】図３０中のＡ−Ａ線及び図３１中のＢ−Ｂ線による断面図である。
【図３６】第１中間形状を有限要素分割した解析モデルを表示する画面の説明図である。
【図３７】解析用基礎データ作成サブプログラムの動作を示すフローチャートである。
【図３８】図２９の領域Ｐの拡大図である。
【図３９】図３１の各要素を回転させた要素分離モデルの説明図である。
【図４０】図３２の各要素をｘｙ平面に投影したときの説明図である。
【図４１】拘束条件に係る節点の扱いの説明図である。
【図４２】解析用データ作成プログラムの動作を示すフローチャートである。
【図４３】２節点間にバネ要素を設定した例を示す説明図である。
【図４４】第１中間形状の平面展開モデルを表示する画面の説明図である。
【図４５】面内歪分布の説明図である。
【図４６】曲げ歪の計算に係る要素の説明図である。
【図４７】座標変換の説明図である。
【図４８】ＸＹ座標系における角度差の説明図である。
【図４９】曲げ歪算出過程の説明図である。
【図５０】曲げ歪分布の説明図である。
【図５１】トータルの歪分布の説明図である。
【図５２】塑性変形を含む時の応力分布の説明図である。
【図５３】工程に複数の中間形状の設計案があるときの成形経路の説明図である。
【図５４】評価用データ作成サブプログラムの動作を示すフローチャートである。
【図５５】第２中間形状の平面展開モデルを表示する画面の説明図である。
【図５６】第１、第２中間形状による平面展開モデルを重ね合わせたときの画面の説明図である。
【図５７】評価点マッピングの説明図である。
【図５８】評価点のｒｓ座標系における位置の説明図である。
【図５９】設計案の組合せ例と各設計案における応力算出の説明図である。
【図６０】応力―歪特性の非線形性を表すグラフである。
【図６１】主歪による安全性評価の例を示すグラフである。
【図６２】残留応力分布の説明図である。
【図６３】トリミングラインの説明図である。
【符号の説明】
【０１７６】
１０        コンピュータ
１１        中央処理装置
１２        入力装置
１３        データ読み込み装置
１４        記憶装置
１５        表示装置
【図１】
2006185228.tif 000024
【図２】
2006185228.tif 000025
【図３】
2006185228.tif 000026
【図４】
2006185228.tif 000027
【図５】
2006185228.tif 000028
【図６】
2006185228.tif 000029
【図７】
2006185228.tif 000030
【図８】
2006185228.tif 000031
【図９】
2006185228.tif 000032
【図１０】
2006185228.tif 000033
【図１１】
2006185228.tif 000034
【図１２】
2006185228.tif 000035
【図１３】
2006185228.tif 000036
【図１４】
2006185228.tif 000037
【図１５】
2006185228.tif 000038
【図１６】
2006185228.tif 000039
【図１７】
2006185228.tif 000040
【図１８】
2006185228.tif 000041
【図１９】
2006185228.tif 000042
【図２０】
2006185228.tif 000043
【図２１】
2006185228.tif 000044
【図２２】
2006185228.tif 000045
【図２３】
2006185228.tif 000046
【図２４】
2006185228.tif 000047
【図２５】
2006185228.tif 000048
【図２６】
2006185228.tif 000049
【図２７】
2006185228.tif 000050
【図２８】
2006185228.tif 000051
【図２９】
2006185228.tif 000052
【図３０】
2006185228.tif 000053
【図３１】
2006185228.tif 000054
【図３２】
2006185228.tif 000055
【図３３】
2006185228.tif 000056
【図３４】
2006185228.tif 000057
【図３５】
2006185228.tif 000058
【図３６】
2006185228.tif 000059
【図３７】
2006185228.tif 000060
【図３８】
2006185228.tif 000061
【図３９】
2006185228.tif 000062
【図４０】
2006185228.tif 000063
【図４１】
2006185228.tif 000064
【図４２】
2006185228.tif 000065
【図４３】
2006185228.tif 000066
【図４４】
2006185228.tif 000067
【図４５】
2006185228.tif 000068
【図４６】
2006185228.tif 000069
【図４７】
2006185228.tif 000070
【図４８】
2006185228.tif 000071
【図４９】
2006185228.tif 000072
【図５０】
2006185228.tif 000073
【図５１】
2006185228.tif 000074
【図５２】
2006185228.tif 000075
【図５３】
2006185228.tif 000076
【図５４】
2006185228.tif 000077
【図５５】
2006185228.tif 000078
【図５６】
2006185228.tif 000079
【図５７】
2006185228.tif 000080
【図５８】
2006185228.tif 000081
【図５９】
2006185228.tif 000082
【図６０】
2006185228.tif 000083
【図６１】
2006185228.tif 000084
【図６２】
2006185228.tif 000085
【図６３】
2006185228.tif 000086