JP7541241B2 - プレス成形シミュレーション方法 - Google Patents
プレス成形シミュレーション方法 Download PDFInfo
- Publication number
- JP7541241B2 JP7541241B2 JP2021018149A JP2021018149A JP7541241B2 JP 7541241 B2 JP7541241 B2 JP 7541241B2 JP 2021018149 A JP2021018149 A JP 2021018149A JP 2021018149 A JP2021018149 A JP 2021018149A JP 7541241 B2 JP7541241 B2 JP 7541241B2
- Authority
- JP
- Japan
- Prior art keywords
- press
- distribution
- model
- linear spring
- spring constant
- 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.)
- Active
Links
Images
Landscapes
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Description
本発明は、プレス成形シミュレーション方法に関する。
従来、プレス成形時におけるプレス金型に見込まれる変形をプレス金型の設計に反映させるため、プレス金型の局所の変形を実測する方法があった。
しかしながら、従来の方法は、特定のプレス金型を使用した場合の、その特定のプレス金型の局所の変形を実測する方法であり、プレス金型の設計にプレス機の剛性を見込んだ方法ではなかった。
しかしながら、従来の方法は、特定のプレス金型を使用した場合の、その特定のプレス金型の局所の変形を実測する方法であり、プレス金型の設計にプレス機の剛性を見込んだ方法ではなかった。
本発明は、上記背景技術の問題点に鑑み、計算負荷を抑えて精度良くプレス成形を再現できるプレス成形シミュレーション方法を提供することを課題とする。
本発明の要旨は以下の通りである。
（１）本発明の一態様に係るプレス成形シミュレーション方法は、プレス金型を備えたプレス機によってワークを所定の荷重でプレス加工してプレス成形品を成形する際におけるワークの物理量の変化を予測するプレス成形シミュレーション方法であって、前記プレス機を、前記プレス金型を弾性支持する弾性支持部と、前記弾性支持部を剛体支持する剛体支持部に分類する分類工程と、前記剛体支持部を剛体としてモデル化し、前記弾性支持部を前記剛体支持部で支持された弾性体としてモデル化した解析モデルを準備する解析モデル準備工程と、前記弾性支持部の平面領域を有限個に分割した各分割領域における変位の分布である離散化変位分布を、前記各分割領域における初期線形ばね定数の分布である初期線形ばね定数分布と前記荷重との積で表した離散化たわみモデルを準備する離散化たわみモデル準備工程と、前記解析モデルによるプレス成形シミュレーション結果又は前記プレス機によるプレス成形試験における実測結果に基づく、前記各分割領域における変位の分布である参照変位分布を準備する参照モデル準備工程と、前記参照変位分布と前記離散化変位分布との差分に基づいて、前記初期線形ばね定数分布を補正して、有限個の線形ばね定数の分布である線形ばね定数分布を計算する線形ばね定数分布計算工程と、を含む。
（２）上記（１）において、前記弾性支持部は、前記プレス金型に直接接するスライド又はボルスタを含んでよい。
（３）上記（１）又は（２）において、前記剛体支持部は、クラウンを含んでよい。
（４）上記（１）から（３）のいずれかにおいて、前記線形ばね定数分布と前記荷重との積で表された修正離散化たわみモデルによって前記ワークをモデル化したワークモデルに前記荷重を作用させた際の前記ワークモデルの弾塑性解析を行う弾塑性解析工程と、を含んでよい。
（２）上記（１）において、前記弾性支持部は、前記プレス金型に直接接するスライド又はボルスタを含んでよい。
（３）上記（１）又は（２）において、前記剛体支持部は、クラウンを含んでよい。
（４）上記（１）から（３）のいずれかにおいて、前記線形ばね定数分布と前記荷重との積で表された修正離散化たわみモデルによって前記ワークをモデル化したワークモデルに前記荷重を作用させた際の前記ワークモデルの弾塑性解析を行う弾塑性解析工程と、を含んでよい。
本発明によれば、計算負荷を抑えて精度良くプレス成形を再現できるプレス成形シミュレーション方法を提供できる。
図１は、プレス機１０とプレス金型２０とを備えたプレス成形装置１の説明図である。
図１に示すように、プレス成形装置１は、被加工体となるワークＷ（例えば、金属材料である。ブランクともいう。）に強い圧力を加え、ワークＷをプレス金型の形状に変形させてプレス成形品を製造する機械である。一般的に、プレス成形装置１は、プレス機１０と、プレス機１０に取り付けられたプレス金型２０とを備えている。
図１に示すように、プレス成形装置１は、被加工体となるワークＷ（例えば、金属材料である。ブランクともいう。）に強い圧力を加え、ワークＷをプレス金型の形状に変形させてプレス成形品を製造する機械である。一般的に、プレス成形装置１は、プレス機１０と、プレス機１０に取り付けられたプレス金型２０とを備えている。
プレス機１０は、下金型２１が上部に取り付けられるボルスタ１１と、上金型２２が下部に取り付けられるスライド１２と、ボルスタ１１とスライド１２とを近づけたり遠ざけたりするように駆動する駆動装置１３と、スライド１２を支持するクラウン１４と、ボルスタ１１を支持するベッド１５と、を含んでいる。
プレス金型２０は、下金型２１及び上金型２２を含んでいる。
プレス機１０の構造は、そのプレス機１０のプレス能力、仕様、設計コンセプト等によって異なる。プレス能力、仕様、設計コンセプト等が同じであっても、プレス機１０の構造には、製造誤差等に起因する、機体差がある。また、プレス機１０に取り付けられるプレス金型２０の構造は、製造するプレス成形品の仕様により異なる。したがって、プレス機１０の構造、プレス金型２０の構造の相違により、プレス荷重が同じであっても、プレス機１０の剛性が異なるので、プレス機１０の変形の状態は異なる。
ところで、プレス成形によってプレス成形品を精度良く製造するため、プレス成形品の被加工体となるワークＷに接する成形面を含むプレス金型２０の形状に、あらかじめプレス機１０のたわみω（鉛直方向の変位の平面分布を意味する。以下、変位又は変形という場合もある。）を見込んでおくことが必要となる場合がある。この場合は、荷重Ｆ（プレス荷重）とプレス機１０のたわみωとの関係であるプレス機１０の剛性を把握しておくことが好ましい。
しかしながら、プレス機１０の剛性を把握するため、プレス機１０のたわみωを精度良く実測するには、測定点が多くなるので、相応の設備等のコストを要する。また、プレス機１０の剛性を、ＦＥＭ等により計算しても、実際の剛性との乖離があったり、計算の負荷が過大となったりする。
しかしながら、プレス機１０の剛性を把握するため、プレス機１０のたわみωを精度良く実測するには、測定点が多くなるので、相応の設備等のコストを要する。また、プレス機１０の剛性を、ＦＥＭ等により計算しても、実際の剛性との乖離があったり、計算の負荷が過大となったりする。
そこで、本発明は、プレス金型を備えたプレス機によってワークを所定の荷重でプレス加工してプレス成形品を成形する際におけるワークの物理量の変化を予測するプレス成形シミュレーション方法である。そして、本発明のプレス成形シミュレーション方法は、プレス機を、プレス金型を弾性支持する弾性支持部と、弾性支持部を剛体支持する剛体支持部に分類する分類工程を含む。また、本発明のプレス成形シミュレーション方法は、剛体支持部を剛体としてモデル化し、プレス金型及び弾性支持部を剛体支持部で支持された弾性体としてモデル化した解析モデルを準備する解析モデル準備工程を含む。また、本発明のプレス成形シミュレーション方法は、弾性支持部の平面領域を有限個に分割した各分割領域における変位の分布である離散化変位分布を、各分割領域における初期線形ばね定数と荷重との積で表した離散化たわみモデルを準備する離散化たわみモデル準備工程を含む。また、本発明のプレス成形シミュレーション方法は、解析モデルによるプレス成形シミュレーション結果又はプレス機によるプレス成形試験における実測結果に基づく、各分割領域における変位の分布である参照変位分布を準備する参照モデル準備工程を含む。また、本発明のプレス成形シミュレーション方法は、参照変位分布と離散化変位分布との差分に基づいて、初期線形ばね定数分布を補正して、有限個の線形ばね定数の分布である線形ばね定数分布を計算する線形ばね定数分布計算工程を含む。このように、本発明のプレス成形シミュレーションによれば、有限個の線形ばね定数の分布で表された、簡略化された離散化たわみモデルで実際のプレス成形装置を模擬するので、計算負荷を小さくでき、しかも、実際のプレス成形装置１の剛性に十分近いモデルを作成することができる。よって、計算負荷を抑えて精度良くプレス成形をシミュレーションできる。
以下、本発明の実施形態を詳細に説明する。
以下、本発明の実施形態を詳細に説明する。
（実施形態）
図２は、解析モデル１００Ｍの説明図である。図３は、離散化たわみモデルＭを示す説明図である。図３（Ａ）は平面視における分割領域Ａｎを示す説明図である。図３（Ｂ）は正面視における離散化たわみモデルＭを示す説明図である。図４は、プレス成形シミュレーション方法のフロー図である。図５は、弾性支持部モデル１０ＡＭの説明図である。図６は、線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓの概念図である。図７は、線形ばね定数分布計算工程Ｓ５のフロー図である。
以下、本実施形態に係るプレス成形シミュレーション方法を、主に、図４に示すフロー図に沿って説明する。
図２は、解析モデル１００Ｍの説明図である。図３は、離散化たわみモデルＭを示す説明図である。図３（Ａ）は平面視における分割領域Ａｎを示す説明図である。図３（Ｂ）は正面視における離散化たわみモデルＭを示す説明図である。図４は、プレス成形シミュレーション方法のフロー図である。図５は、弾性支持部モデル１０ＡＭの説明図である。図６は、線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓの概念図である。図７は、線形ばね定数分布計算工程Ｓ５のフロー図である。
以下、本実施形態に係るプレス成形シミュレーション方法を、主に、図４に示すフロー図に沿って説明する。
（プレス成形シミュレーション方法）
本実施形態に係るプレス成形シミュレーション方法は、プレス金型２０を備えたプレス機１０によってワークＷを所定の荷重でプレス加工してプレス成形品を成形する際におけるワークＷの物理量の変化を予測するプレス成形シミュレーション方法である。
図４に示すように、本実施形態に係るプレス成形シミュレーション方法は、分類工程Ｓ１と、解析モデル準備工程Ｓ２と、離散化たわみモデル準備工程Ｓ３と、参照モデル準備工程Ｓ４と、線形ばね定数分布計算工程Ｓ５と、を含んでいる。
また、線形ばね定数分布計算工程Ｓ５で計算した線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓに基づいて離散化たわみモデル準備工程Ｓ３で準備した離散化たわみモデルＭを修正して線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表された修正離散化たわみモデルＭＭを作成する修正離散化たわみモデル作成工程Ｓ６を含んでよい。
また、線形ばね定数分布と荷重Ｆとの積で表された修正離散化たわみモデルＭＭによってワークＷに荷重Ｆを作用させた際のワークＷの弾塑性解析を行う弾塑性解析工程Ｓ７を含んでよい。
以下、プレス成形シミュレーション方法の各工程を個別に説明する。以下の各工程において、特に説明のない限り、分類工程Ｓ１において、プレス機１０のスライド１２及びボルスタ１１を弾性支持部１０Ａとし、その他の部分を剛体支持部１０Ｂとして分類した場合で説明する。
本実施形態に係るプレス成形シミュレーション方法は、プレス金型２０を備えたプレス機１０によってワークＷを所定の荷重でプレス加工してプレス成形品を成形する際におけるワークＷの物理量の変化を予測するプレス成形シミュレーション方法である。
図４に示すように、本実施形態に係るプレス成形シミュレーション方法は、分類工程Ｓ１と、解析モデル準備工程Ｓ２と、離散化たわみモデル準備工程Ｓ３と、参照モデル準備工程Ｓ４と、線形ばね定数分布計算工程Ｓ５と、を含んでいる。
また、線形ばね定数分布計算工程Ｓ５で計算した線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓに基づいて離散化たわみモデル準備工程Ｓ３で準備した離散化たわみモデルＭを修正して線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表された修正離散化たわみモデルＭＭを作成する修正離散化たわみモデル作成工程Ｓ６を含んでよい。
また、線形ばね定数分布と荷重Ｆとの積で表された修正離散化たわみモデルＭＭによってワークＷに荷重Ｆを作用させた際のワークＷの弾塑性解析を行う弾塑性解析工程Ｓ７を含んでよい。
以下、プレス成形シミュレーション方法の各工程を個別に説明する。以下の各工程において、特に説明のない限り、分類工程Ｓ１において、プレス機１０のスライド１２及びボルスタ１１を弾性支持部１０Ａとし、その他の部分を剛体支持部１０Ｂとして分類した場合で説明する。
（分類工程）
まず、プレス機１０を、プレス金型２０を弾性支持する弾性支持部１０Ａと、弾性支持部１０Ａを剛体支持する剛体支持部１０Ｂに分類する（分類工程Ｓ１）。
まず、プレス機１０を、プレス金型２０を弾性支持する弾性支持部１０Ａと、弾性支持部１０Ａを剛体支持する剛体支持部１０Ｂに分類する（分類工程Ｓ１）。
詳細には、例えば、弾性支持部１０Ａは、プレス金型２０に直接接するスライド１２又はボルスタ１１を含んでよい。これにより、プレス金型２０の変形に影響を及ぼすプレス機１０の範囲を部分的に弾性体としてモデル化できる。よって、プレス成形シミュレーション方法の計算負荷を軽減できる。なお、弾性支持部１０Ａに、ボルスタ１１を支持するベッド１５を含めてもよい。
具体的には、例えば、プレス機１０のスライド１２及びボルスタ１１を、弾性支持部１０Ａに分類し、その他の部分を、剛体支持部１０Ｂに分類する。
具体的には、例えば、プレス機１０のスライド１２及びボルスタ１１を、弾性支持部１０Ａに分類し、その他の部分を、剛体支持部１０Ｂに分類する。
また、例えば、剛体支持部１０Ｂは、クラウン１４を含んでよい。これにより、プレス金型２０の変形にほとんど影響を及ぼさないことが見込まれるプレス機１０の範囲を部分的に剛体としてモデル化できる。よって、プレス成形シミュレーション方法の計算負荷を軽減できる。
（解析モデル準備工程）
次に、図２に示すように、剛体支持部１０Ｂを剛体としてモデル化し、弾性支持部１０Ａを剛体支持部１０Ｂで支持された弾性体としてモデル化した解析モデル１００Ｍを準備する（解析モデル準備工程Ｓ２）。
詳細には、弾性支持部１０Ａを、例えば、ソリッド等の３次元要素、シェル等の２次元要素又はばね等の１次元要素の弾性体としてモデル化し、弾性支持部モデル１０ＡＭとしてよい。
また、剛体支持部１０Ｂを、弾性支持部モデル１０ＡＭの境界における変位を拘束する（例えば、ヤング率を無限大とした）境界条件としてモデル化し、剛体支持部モデル１０ＢＭとしてよい。剛体支持部モデル１０ＢＭは、弾性支持部モデル１０ＡＭに力が作用する部分において、弾性支持部モデル１０ＡＭの変形を拘束しつつ、弾性支持部モデル１０ＡＭのその部分の移動（変位）に伴って移動（変位）する剛体としてモデル化されたものであってよい。剛体支持部モデル１０ＢＭは、弾性支持部モデル１０ＡＭとの境界に接し、荷重が作用しても変形しない剛体として設定されたソリッド等の３次元要素としてモデル化されたものであってよい。
なお、図２に示すように、解析モデル１００Ｍは、弾性支持部モデル１０ＡＭを、プレス金型２０を、例えば、弾性体としてモデル化したプレス金型モデル２０Ｍで支持されたものとしてよい。解析モデル１００Ｍは、プレス金型モデル２０Ｍの上下両側を解析モデル１００Ｍで支持されたものとして、同時に、上側の解析モデル１００Ｍでスライド１２を再現し、下側の解析モデル１００Ｍでボルスタ１１を再現してよい。
次に、図２に示すように、剛体支持部１０Ｂを剛体としてモデル化し、弾性支持部１０Ａを剛体支持部１０Ｂで支持された弾性体としてモデル化した解析モデル１００Ｍを準備する（解析モデル準備工程Ｓ２）。
詳細には、弾性支持部１０Ａを、例えば、ソリッド等の３次元要素、シェル等の２次元要素又はばね等の１次元要素の弾性体としてモデル化し、弾性支持部モデル１０ＡＭとしてよい。
また、剛体支持部１０Ｂを、弾性支持部モデル１０ＡＭの境界における変位を拘束する（例えば、ヤング率を無限大とした）境界条件としてモデル化し、剛体支持部モデル１０ＢＭとしてよい。剛体支持部モデル１０ＢＭは、弾性支持部モデル１０ＡＭに力が作用する部分において、弾性支持部モデル１０ＡＭの変形を拘束しつつ、弾性支持部モデル１０ＡＭのその部分の移動（変位）に伴って移動（変位）する剛体としてモデル化されたものであってよい。剛体支持部モデル１０ＢＭは、弾性支持部モデル１０ＡＭとの境界に接し、荷重が作用しても変形しない剛体として設定されたソリッド等の３次元要素としてモデル化されたものであってよい。
なお、図２に示すように、解析モデル１００Ｍは、弾性支持部モデル１０ＡＭを、プレス金型２０を、例えば、弾性体としてモデル化したプレス金型モデル２０Ｍで支持されたものとしてよい。解析モデル１００Ｍは、プレス金型モデル２０Ｍの上下両側を解析モデル１００Ｍで支持されたものとして、同時に、上側の解析モデル１００Ｍでスライド１２を再現し、下側の解析モデル１００Ｍでボルスタ１１を再現してよい。
（離散化たわみモデル準備工程）
次に、弾性支持部１０Ａの平面領域Ａを有限個ｎに分割した各分割領域Ａｎにおける変位ω（ｉ，ｊ）の平面的な分布である離散化変位分布ω（ｉ，ｊ）ｄｉｓを、分割領域Ａｎごとの初期線形ばね定数Ｋｉｎｉｔ（ｉ，ｊ）の平面的な分布である初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表した離散化たわみモデルＭを準備する（離散化たわみモデル準備工程Ｓ３）。
次に、弾性支持部１０Ａの平面領域Ａを有限個ｎに分割した各分割領域Ａｎにおける変位ω（ｉ，ｊ）の平面的な分布である離散化変位分布ω（ｉ，ｊ）ｄｉｓを、分割領域Ａｎごとの初期線形ばね定数Ｋｉｎｉｔ（ｉ，ｊ）の平面的な分布である初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表した離散化たわみモデルＭを準備する（離散化たわみモデル準備工程Ｓ３）。
図３は、離散化たわみモデルＭを示す説明図である。図３（Ａ）は、平面視における分割領域を示す説明図である。図３（Ｂ）は正面視における離散化たわみモデルＭを示す説明図である。
詳細には、図３（Ａ）に示すように、弾性支持部１０Ａの平面領域Ａを有限個ｎに分割した各分割領域Ａｎにおける変位ω（ｉ，ｊ）の分布である離散化変位分布ω（ｉ，ｊ）ｄｉｓを、図３（Ｂ）に示すように、分割領域Ａｎごとの初期線形ばね定数Ｋｉｎｉｔ（ｉ，ｊ）の分布である初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表した離散化たわみモデルＭを準備する。なお、図３（Ａ）において、交差する直線同士で囲まれる矩形の領域が分割領域Ａｎを表している。
例えば、離散化たわみモデルＭは、次のような関数で表すことができる。
Ｆ／ｎ＝Ｋｉｎｉｔ（ｉ，ｊ）＊ω（ｉ，ｊ）
Ｆ＝Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓ＊ω（ｉ，ｊ）ｄｉｓ
詳細には、図３（Ａ）に示すように、弾性支持部１０Ａの平面領域Ａを有限個ｎに分割した各分割領域Ａｎにおける変位ω（ｉ，ｊ）の分布である離散化変位分布ω（ｉ，ｊ）ｄｉｓを、図３（Ｂ）に示すように、分割領域Ａｎごとの初期線形ばね定数Ｋｉｎｉｔ（ｉ，ｊ）の分布である初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表した離散化たわみモデルＭを準備する。なお、図３（Ａ）において、交差する直線同士で囲まれる矩形の領域が分割領域Ａｎを表している。
例えば、離散化たわみモデルＭは、次のような関数で表すことができる。
Ｆ／ｎ＝Ｋｉｎｉｔ（ｉ，ｊ）＊ω（ｉ，ｊ）
Ｆ＝Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓ＊ω（ｉ，ｊ）ｄｉｓ
（参照モデル準備工程）
次に、解析モデル１００Ｍによるプレス成形シミュレーション結果又はプレス機１０によるプレス成形試験における実測結果に基づく、弾性支持部１０Ａの平面領域Ａを有限個に分割した各分割領域Ａｎにおける変位ωｒｅｆ（ｉ，ｊ）の分布である参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを準備する（参照モデル準備工程Ｓ４）。以下、荷重Ｆと参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓとの関係を表す関数を、参照モデルＭｒｅｆという場合がある。
参照モデル準備工程Ｓ４は、離散化たわみモデル準備工程Ｓ３と並行して実施してよい。
次に、解析モデル１００Ｍによるプレス成形シミュレーション結果又はプレス機１０によるプレス成形試験における実測結果に基づく、弾性支持部１０Ａの平面領域Ａを有限個に分割した各分割領域Ａｎにおける変位ωｒｅｆ（ｉ，ｊ）の分布である参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを準備する（参照モデル準備工程Ｓ４）。以下、荷重Ｆと参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓとの関係を表す関数を、参照モデルＭｒｅｆという場合がある。
参照モデル準備工程Ｓ４は、離散化たわみモデル準備工程Ｓ３と並行して実施してよい。
詳細には、解析モデル準備工程Ｓ２で準備した解析モデル１００Ｍを用いて、プレス成形シミュレーションを行い、その結果に基づいて、弾性支持部１０Ａの平面領域Ａを有限個に分割した各分割領域Ａｎにおける変位ωｒｅｆ（ｉ，ｊ）の分布である参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを計算する。
例えば、解析モデル１００Ｍに基づいて、図５に示すように、プレス機１０のスライド１２及びボルスタ１１を弾性体としてモデル化して、スライドモデル１２ＡＭ及びボルスタモデル１１ＡＭで構成された弾性支持部モデル１０ＡＭとした参照モデルＭｒｅｆを作成する。なお、この際、弾性支持部モデル１０ＡＭを、例えば、実際のスライド１２及びボルスタ１１の形状を模擬した、３次元要素で構成されたソリッドモデルにしてよい。この解析モデル１００Ｍに基づいて作成した参照モデルＭｒｅｆにおいて、スライド１２の下面にプレス金型モデル２０Ｍの上金型モデル２２Ｍを取り付けた状態にする。また、この解析モデル１００Ｍに基づいて作成した参照モデルＭｒｅｆにおいて、ボルスタ１１の下面にプレス金型モデル２０Ｍの下金型モデル２１Ｍを取り付けた状態にする。この解析モデル１００Ｍに基づいて作成した参照モデルＭｒｅｆにおいて、上金型モデル２２Ｍと下金型モデル２１Ｍとの間に、ワークモデルＷＭを挟んだ状態にする。
そして、実際のプレス機１０を模擬できるように、上方から下方に向けて、スライドモデル１２ＡＭの上面に、荷重Ｆを作用させる。そして、スライドモデル１２ＡＭとボルスタモデル１１ＡＭとの相対的な距離を徐々に短縮して、ワークモデルＷＭに荷重Ｆを作用させて、実際のプレス成形を模擬する。この際、荷重Ｆと、スライドモデル１２ＡＭ又はボルスタモデル１１ＡＭの変位との関係から、各分割領域Ａｎにおける変位ωｒｅｆ（ｉ，ｊ）の分布である参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを計算する。
例えば、解析モデル１００Ｍに基づいて、図５に示すように、プレス機１０のスライド１２及びボルスタ１１を弾性体としてモデル化して、スライドモデル１２ＡＭ及びボルスタモデル１１ＡＭで構成された弾性支持部モデル１０ＡＭとした参照モデルＭｒｅｆを作成する。なお、この際、弾性支持部モデル１０ＡＭを、例えば、実際のスライド１２及びボルスタ１１の形状を模擬した、３次元要素で構成されたソリッドモデルにしてよい。この解析モデル１００Ｍに基づいて作成した参照モデルＭｒｅｆにおいて、スライド１２の下面にプレス金型モデル２０Ｍの上金型モデル２２Ｍを取り付けた状態にする。また、この解析モデル１００Ｍに基づいて作成した参照モデルＭｒｅｆにおいて、ボルスタ１１の下面にプレス金型モデル２０Ｍの下金型モデル２１Ｍを取り付けた状態にする。この解析モデル１００Ｍに基づいて作成した参照モデルＭｒｅｆにおいて、上金型モデル２２Ｍと下金型モデル２１Ｍとの間に、ワークモデルＷＭを挟んだ状態にする。
そして、実際のプレス機１０を模擬できるように、上方から下方に向けて、スライドモデル１２ＡＭの上面に、荷重Ｆを作用させる。そして、スライドモデル１２ＡＭとボルスタモデル１１ＡＭとの相対的な距離を徐々に短縮して、ワークモデルＷＭに荷重Ｆを作用させて、実際のプレス成形を模擬する。この際、荷重Ｆと、スライドモデル１２ＡＭ又はボルスタモデル１１ＡＭの変位との関係から、各分割領域Ａｎにおける変位ωｒｅｆ（ｉ，ｊ）の分布である参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを計算する。
また、詳細には、上述の解析モデル１００Ｍに基づいて作成した参照モデルＭｒｅｆによるプレス成形シミュレーション結果に基づく方法に換えて、実際のプレス機１０によるプレス成形試験において、プレス機１０の変位を実測し、その実測結果に基づいて、弾性支持部１０Ａの平面領域Ａを有限個に分割した各分割領域Ａｎにおける変位ωｒｅｆ（ｉ，ｊ）の分布である参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを計算してもよい。例えば、ワークＷを挟むプレス金型２０を取り付けた実際のプレス機１０に荷重Ｆを作用させるプレス成形試験において、荷重Ｆが作用した際の、プレス機１０におけるスライド１２又はボルスタ１１の鉛直方向の変位を、平面的に散らされた複数の点で、例えば、変位計を用いて実測する。そして、荷重Ｆとそれらの複数点での変位との関係から、各分割領域Ａｎにおける変位ωｒｅｆ（ｉ，ｊ）の分布である参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを計算してもよい。
（線形ばね定数分布計算工程）
そして、参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓと離散化変位分布ω（ｉ，ｊ）ｄｉｓとの差分に基づいて、初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓを補正して、有限個の線形ばね定数の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算する（線形ばね定数分布計算工程Ｓ５）。
詳細には、例えば、離散化変位分布ω（ｉ，ｊ）ｄｉｓから参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを差し引いた値がゼロに近づくように、離散化たわみモデルＭの初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓを修正する。具体的には、例えば、後述する勾配法による収束計算により線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算してよい。これにより、図６に示すような、分割領域Ａｎごとにおける線形ばね定数（Ｚ軸）のＸＹ平面的な分布となる線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算できる。
そして、参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓと離散化変位分布ω（ｉ，ｊ）ｄｉｓとの差分に基づいて、初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓを補正して、有限個の線形ばね定数の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算する（線形ばね定数分布計算工程Ｓ５）。
詳細には、例えば、離散化変位分布ω（ｉ，ｊ）ｄｉｓから参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを差し引いた値がゼロに近づくように、離散化たわみモデルＭの初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓを修正する。具体的には、例えば、後述する勾配法による収束計算により線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算してよい。これにより、図６に示すような、分割領域Ａｎごとにおける線形ばね定数（Ｚ軸）のＸＹ平面的な分布となる線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算できる。
（勾配法による収束計算）
線形ばね定数分布計算工程Ｓ５における線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓは、以下に説明するような、勾配法による収束計算によって求めることができる。
（１）図７に示すように、まず、分割領域Ａｎごとの初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓを仮設定する（初期線形ばね定数仮設定工程Ｓ５１）。
（２）次に、参照モデル準備工程Ｓ４で準備した参照モデルＭｒｅｆにより参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを計算する（参照変位分布計算工程Ｓ５２）。
（３）参照変位分布計算工程Ｓ５２と並行して、離散化たわみモデル準備工程Ｓ３で準備した離散化たわみモデルＭにより離散化変位分布ω（ｉ，ｊ）ｄｉｓを計算する（離散化変位分布計算工程Ｓ５３）。
（４）離散化変位分布ω（ｉ，ｊ）ｄｉｓと参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓとの差分の絶対値が、閾値ｔｏｌを超えるか否かを判定する（閾値判定工程Ｓ５４）。
離散化変位分布ω（ｉ，ｊ）ｄｉｓと参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓとの差分の絶対値である誤差分布ｅ（ｉ，ｊ）ｄｉｓが、閾値ｔｏｌを超える場合（例えば、分割領域Ａｎごとの変位ω（ｉ，ｊ）と分割領域Ａｎごとの参照変位ωｒｅｆ（ｉ，ｊ）との差分の絶対値である、分割領域Ａｎごとの誤差ｅ（ｉ，ｊ）が、閾値ｔｏｌを超える場合）、分割領域Ａｎごとの誤差ｅ（ｉ，ｊ）に任意の補正係数を乗じた値を、その分割領域Ａｎでの初期線形ばね定数Ｋｉｎｉｔ（ｉ，ｊ）に対して（２回目以降は修正された線形ばね定数Ｋ（ｉ，ｊ）に対して）加減して、分割領域Ａｎでの修正された線形ばね定数Ｋ（ｉ，ｊ）を得る。そして、分割領域Ａｎごとの修正された線形ばね定数Ｋ（ｉ，ｊ）の集合となる、修正された線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを得る。続いて、参照変位分布計算工程Ｓ５２及び離散化変位分布計算工程Ｓ５３を再び実施する。
離散化変位分布ω（ｉ，ｊ）ｄｉｓと参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓとの差分の絶対値である誤差分布ｅ（ｉ，ｊ）ｄｉｓが、閾値ｔｏｌ以下となった場合、後述の修正離散化たわみモデル作成工程Ｓ６に移行する。
（５）これらの工程を、離散化変位分布ω（ｉ，ｊ）ｄｉｓと参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓとの差分の絶対値である誤差分布ｅ（ｉ，ｊ）ｄｉｓが、閾値ｔｏｌ以下となるまで繰り返す。そして、閾値ｔｏｌ以下となった際の修正された線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓが、求める線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓとなる。よって、実際のプレス機１０の剛性を、有限個の線形ばね定数の平面的な分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを含む離散化たわみモデルＭで表すことができる。
線形ばね定数分布計算工程Ｓ５における線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓは、以下に説明するような、勾配法による収束計算によって求めることができる。
（１）図７に示すように、まず、分割領域Ａｎごとの初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓを仮設定する（初期線形ばね定数仮設定工程Ｓ５１）。
（２）次に、参照モデル準備工程Ｓ４で準備した参照モデルＭｒｅｆにより参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを計算する（参照変位分布計算工程Ｓ５２）。
（３）参照変位分布計算工程Ｓ５２と並行して、離散化たわみモデル準備工程Ｓ３で準備した離散化たわみモデルＭにより離散化変位分布ω（ｉ，ｊ）ｄｉｓを計算する（離散化変位分布計算工程Ｓ５３）。
（４）離散化変位分布ω（ｉ，ｊ）ｄｉｓと参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓとの差分の絶対値が、閾値ｔｏｌを超えるか否かを判定する（閾値判定工程Ｓ５４）。
離散化変位分布ω（ｉ，ｊ）ｄｉｓと参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓとの差分の絶対値である誤差分布ｅ（ｉ，ｊ）ｄｉｓが、閾値ｔｏｌを超える場合（例えば、分割領域Ａｎごとの変位ω（ｉ，ｊ）と分割領域Ａｎごとの参照変位ωｒｅｆ（ｉ，ｊ）との差分の絶対値である、分割領域Ａｎごとの誤差ｅ（ｉ，ｊ）が、閾値ｔｏｌを超える場合）、分割領域Ａｎごとの誤差ｅ（ｉ，ｊ）に任意の補正係数を乗じた値を、その分割領域Ａｎでの初期線形ばね定数Ｋｉｎｉｔ（ｉ，ｊ）に対して（２回目以降は修正された線形ばね定数Ｋ（ｉ，ｊ）に対して）加減して、分割領域Ａｎでの修正された線形ばね定数Ｋ（ｉ，ｊ）を得る。そして、分割領域Ａｎごとの修正された線形ばね定数Ｋ（ｉ，ｊ）の集合となる、修正された線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを得る。続いて、参照変位分布計算工程Ｓ５２及び離散化変位分布計算工程Ｓ５３を再び実施する。
離散化変位分布ω（ｉ，ｊ）ｄｉｓと参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓとの差分の絶対値である誤差分布ｅ（ｉ，ｊ）ｄｉｓが、閾値ｔｏｌ以下となった場合、後述の修正離散化たわみモデル作成工程Ｓ６に移行する。
（５）これらの工程を、離散化変位分布ω（ｉ，ｊ）ｄｉｓと参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓとの差分の絶対値である誤差分布ｅ（ｉ，ｊ）ｄｉｓが、閾値ｔｏｌ以下となるまで繰り返す。そして、閾値ｔｏｌ以下となった際の修正された線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓが、求める線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓとなる。よって、実際のプレス機１０の剛性を、有限個の線形ばね定数の平面的な分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを含む離散化たわみモデルＭで表すことができる。
（修正離散化たわみモデル作成工程）
次に、線形ばね定数分布計算工程Ｓ５で計算された線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表された修正離散化たわみモデルＭＭを作成する（修正離散化たわみモデル作成工程Ｓ６）。
このようにして、修正離散化たわみモデルＭＭを作成できる。したがって、プレス機１０を、離散化たわみモデルＭで表すことができるので、計算負荷を低減できるとともに、実際のプレス機１０の変形状態をより正確に再現できる。よって、プレス金型２０の設計、プレス成形品の設計等に、計算負荷を抑制しつつ、プレス機１０の剛性の影響を見込むことができる。
次に、線形ばね定数分布計算工程Ｓ５で計算された線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表された修正離散化たわみモデルＭＭを作成する（修正離散化たわみモデル作成工程Ｓ６）。
このようにして、修正離散化たわみモデルＭＭを作成できる。したがって、プレス機１０を、離散化たわみモデルＭで表すことができるので、計算負荷を低減できるとともに、実際のプレス機１０の変形状態をより正確に再現できる。よって、プレス金型２０の設計、プレス成形品の設計等に、計算負荷を抑制しつつ、プレス機１０の剛性の影響を見込むことができる。
（弾塑性解析工程）
修正離散化たわみモデル作成工程Ｓ６で作成した、修正離散化たわみモデルＭＭにより、ワークＷをモデル化したワークモデルＷＭに対して弾塑性解析を行う（弾塑性解析工程Ｓ７）。
詳細には、線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表された修正離散化たわみモデルＭＭによって、弾塑性体としてモデル化されたワークモデルＷＭに荷重Ｆを作用させた際のワークモデルＷＭの弾塑性解析を行う。
これにより、ワークモデルＷＭの変形に、計算負荷を抑制しつつ、プレス機１０の剛性の影響を見込むことができる。よって、精度良くプレス成形品を成形できるプレス金型２０を設計でき、精度良くプレス成形品を設計できる。
修正離散化たわみモデル作成工程Ｓ６で作成した、修正離散化たわみモデルＭＭにより、ワークＷをモデル化したワークモデルＷＭに対して弾塑性解析を行う（弾塑性解析工程Ｓ７）。
詳細には、線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表された修正離散化たわみモデルＭＭによって、弾塑性体としてモデル化されたワークモデルＷＭに荷重Ｆを作用させた際のワークモデルＷＭの弾塑性解析を行う。
これにより、ワークモデルＷＭの変形に、計算負荷を抑制しつつ、プレス機１０の剛性の影響を見込むことができる。よって、精度良くプレス成形品を成形できるプレス金型２０を設計でき、精度良くプレス成形品を設計できる。
（プレス成形シミュレーションプログラム）
次に、上述のプレス成形シミュレーション方法を、コンピュータに実行させるプレス成形シミュレーションプログラムについて説明する。
プレス成形シミュレーションプログラムは、コンピュータに、参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓと離散化変位分布ω（ｉ，ｊ）ｄｉｓとの差分に基づいて、有限個ｎの線形ばね定数Ｋ（ｉ，ｊ）の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算する線形ばね定数分布計算ステップを実行させるものである。なお、参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓ及び離散化変位分布ω（ｉ，ｊ）ｄｉｓのデータは、プレス成形シミュレーションプログラムの実行時に呼び出せるように、事前にメモリに記憶しておく。
このように、プレス成形シミュレーションプログラムは、有限個ｎの線形ばね定数Ｋ（ｉ，ｊ）の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算できるので、ＣＡＥに利用でき、計算負荷を抑制しつつ、プレス機１０の変形を見込んでプレス金型２０を設計できる。よって、プレス機１０の変形を見込んで設計されたプレス金型２０により、プレス成形品を精度良く製造できる。
次に、上述のプレス成形シミュレーション方法を、コンピュータに実行させるプレス成形シミュレーションプログラムについて説明する。
プレス成形シミュレーションプログラムは、コンピュータに、参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓと離散化変位分布ω（ｉ，ｊ）ｄｉｓとの差分に基づいて、有限個ｎの線形ばね定数Ｋ（ｉ，ｊ）の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算する線形ばね定数分布計算ステップを実行させるものである。なお、参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓ及び離散化変位分布ω（ｉ，ｊ）ｄｉｓのデータは、プレス成形シミュレーションプログラムの実行時に呼び出せるように、事前にメモリに記憶しておく。
このように、プレス成形シミュレーションプログラムは、有限個ｎの線形ばね定数Ｋ（ｉ，ｊ）の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算できるので、ＣＡＥに利用でき、計算負荷を抑制しつつ、プレス機１０の変形を見込んでプレス金型２０を設計できる。よって、プレス機１０の変形を見込んで設計されたプレス金型２０により、プレス成形品を精度良く製造できる。
（プレス成形シミュレーション装置）
次に、上述のプレス成形シミュレーションプログラムを搭載した、プレス成形シミュレーション装置について説明する。
プレス成形シミュレーション装置は、一般的なコンピュータと同様に、ストレージ、プロセッサ、メモリ、インタフェース及び出力装置を備えている。
詳細には、プレス成形シミュレーション装置は、有限個ｎの線形ばね定数Ｋ（ｉ，ｊ）の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを表示する出力装置を備えている。
このように、プレス成形シミュレーション装置は、有限個ｎの線形ばね定数Ｋ（ｉ，ｊ）の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを表示する出力装置を備えているので、ＣＡＥに利用でき、プレス機１０の変形を見込んでプレス金型２０を設計できる。よって、プレス機１０の変形を見込んで設計されたプレス金型２０により、プレス成形品を精度良く製造できる。
次に、上述のプレス成形シミュレーションプログラムを搭載した、プレス成形シミュレーション装置について説明する。
プレス成形シミュレーション装置は、一般的なコンピュータと同様に、ストレージ、プロセッサ、メモリ、インタフェース及び出力装置を備えている。
詳細には、プレス成形シミュレーション装置は、有限個ｎの線形ばね定数Ｋ（ｉ，ｊ）の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを表示する出力装置を備えている。
このように、プレス成形シミュレーション装置は、有限個ｎの線形ばね定数Ｋ（ｉ，ｊ）の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを表示する出力装置を備えているので、ＣＡＥに利用でき、プレス機１０の変形を見込んでプレス金型２０を設計できる。よって、プレス機１０の変形を見込んで設計されたプレス金型２０により、プレス成形品を精度良く製造できる。
（その他の実施形態）
以上、図面を参照して実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、様々な設計変更等をすることが可能である。
以上、図面を参照して実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、様々な設計変更等をすることが可能である。
本実施形態に係るプレス成形シミュレーション方法は、プレス金型２０を備えたプレス機１０によってワークＷを所定の荷重Ｆでプレス加工してプレス成形品を成形する際におけるワークＷの物理量の変化を予測するプレス成形シミュレーション方法である。プレス成形シミュレーション方法は、プレス機１０を、プレス金型２０を弾性支持する弾性支持部１０Ａと、弾性支持部１０Ａを剛体支持する剛体支持部１０Ｂに分類する分類工程Ｓ１と、剛体支持部１０Ｂを剛体としてモデル化し、弾性支持部１０Ａを剛体支持部１０Ｂで支持された弾性体としてモデル化した解析モデル１００Ｍを準備する解析モデル準備工程Ｓ２と、弾性支持部１０Ａの平面領域Ａを有限個ｎに分割した各分割領域Ａｎにおける変位の分布である離散化変位分布ω（ｉ，ｊ）ｄｉｓを、分割領域Ａｎごとの初期線形ばね定数Ｋｉｎｉｔ（ｉ，ｊ）の分布である初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓと荷重Ｆとの積で表した離散化たわみモデルＭを準備する離散化たわみモデル準備工程Ｓ３と、解析モデル１００Ｍによるプレス成形シミュレーション結果又はプレス機１０によるプレス成形試験における実測結果に基づく、各分割領域Ａｎにおける変位の分布である参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓを準備する参照モデル準備工程Ｓ４と、参照変位分布ωｒｅｆ（ｉ，ｊ）ｄｉｓと離散化変位分布ω（ｉ，ｊ）ｄｉｓとの差分に基づいて、初期線形ばね定数分布Ｋｉｎｉｔ（ｉ，ｊ）ｄｉｓを補正して、有限個ｎの線形ばね定数の分布である線形ばね定数分布Ｋ（ｉ，ｊ）ｄｉｓを計算する線形ばね定数分布計算工程Ｓ５と、を含む。これにより、有限個ｎの線形ばね定数Ｋ（ｉ，ｊ）の分布で表された、簡略化された離散化たわみモデルＭで実際のプレス成形装置１を模擬するので、計算負荷を小さくでき、しかも、実際のプレス成形装置１の剛性に十分近いモデルを作成することができる。よって、計算負荷を抑えて精度良くプレス成形を再現できるプレス成形シミュレーション方法を提供できる。
１ プレス成形装置
１０ プレス機
１０Ａ 弾性支持部
１０ＡＭ 弾性支持部モデル
１０ＢＭ 剛体支持部モデル
１１ ボルスタ
１１ＡＭ ボルスタモデル
１２ スライド
１２ＡＭ スライドモデル
１３ 駆動装置
１４ クラウン
１５ ベッド
２０ プレス金型
２０Ｍ プレス金型モデル
２１ 下金型
２１Ｍ 下金型モデル
２２ 上金型
２２Ｍ 上金型モデル
１００Ｍ 解析モデル
Ａ 平面領域
Ａｎ 分割領域
Ｆ 荷重
Ｍ 離散化たわみモデル
Ｍｒｅｆ 参照モデル
Ｓ１ 分類工程
Ｓ２ 解析モデル準備工程
Ｓ３ モデル準備工程
Ｓ４ 参照モデル準備工程
Ｓ５ 定数分布計算工程
Ｓ６ モデル作成工程
Ｓ７ 弾塑性解析工程
Ｓ５１ 初期線形ばね定数仮設定工程
Ｓ５２ 参照変位分布計算工程
Ｓ５３ 離散化変位分布計算工程
Ｓ５４ 閾値判定工程
Ｗ ワーク
ＷＭ ワークモデル
１０ プレス機
１０Ａ 弾性支持部
１０ＡＭ 弾性支持部モデル
１０ＢＭ 剛体支持部モデル
１１ ボルスタ
１１ＡＭ ボルスタモデル
１２ スライド
１２ＡＭ スライドモデル
１３ 駆動装置
１４ クラウン
１５ ベッド
２０ プレス金型
２０Ｍ プレス金型モデル
２１ 下金型
２１Ｍ 下金型モデル
２２ 上金型
２２Ｍ 上金型モデル
１００Ｍ 解析モデル
Ａ 平面領域
Ａｎ 分割領域
Ｆ 荷重
Ｍ 離散化たわみモデル
Ｍｒｅｆ 参照モデル
Ｓ１ 分類工程
Ｓ２ 解析モデル準備工程
Ｓ３ モデル準備工程
Ｓ４ 参照モデル準備工程
Ｓ５ 定数分布計算工程
Ｓ６ モデル作成工程
Ｓ７ 弾塑性解析工程
Ｓ５１ 初期線形ばね定数仮設定工程
Ｓ５２ 参照変位分布計算工程
Ｓ５３ 離散化変位分布計算工程
Ｓ５４ 閾値判定工程
Ｗ ワーク
ＷＭ ワークモデル
Claims (4)
- プレス金型を備えたプレス機によってワークを所定の荷重でプレス加工してプレス成形品を成形する際におけるワークの物理量の変化を予測するプレス成形シミュレーション方法であって、
前記プレス機を、前記プレス金型を弾性支持する弾性支持部と、前記弾性支持部を剛体支持する剛体支持部に分類する分類工程と、
前記剛体支持部を剛体としてモデル化し、前記弾性支持部を前記剛体支持部で支持された弾性体としてモデル化した解析モデルを準備する解析モデル準備工程と、
前記弾性支持部の平面領域を有限個に分割した各分割領域における変位の分布である離散化変位分布を、前記各分割領域における初期線形ばね定数の分布である初期線形ばね定数分布と前記荷重との積で表した離散化たわみモデルを準備する離散化たわみモデル準備工程と、
前記解析モデルによるプレス成形シミュレーション結果又は前記プレス機によるプレス成形試験における実測結果に基づく、前記各分割領域における変位の分布である参照変位分布を準備する参照モデル準備工程と、
前記参照変位分布と前記離散化変位分布との差分に基づいて、前記初期線形ばね定数分布を補正して、有限個の線形ばね定数の分布である線形ばね定数分布を計算する線形ばね定数分布計算工程と、
を含むことを特徴とするプレス成形シミュレーション方法。 - 前記弾性支持部は、前記プレス金型に直接接するスライド又はボルスタを含む
ことを特徴とする請求項１に記載のプレス成形シミュレーション方法。 - 前記剛体支持部は、クラウンを含む
ことを特徴とする請求項１又は請求項２に記載のプレス成形シミュレーション方法。 - 前記線形ばね定数分布と前記荷重との積で表された修正離散化たわみモデルによって前記ワークをモデル化したワークモデルに前記荷重を作用させた際の前記ワークモデルの弾塑性解析を行う弾塑性解析工程と、を含む
ことを特徴とする請求項１から請求項３のいずれか１項に記載のプレス成形シミュレーション方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021018149A JP7541241B2 (ja) | 2021-02-08 | 2021-02-08 | プレス成形シミュレーション方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021018149A JP7541241B2 (ja) | 2021-02-08 | 2021-02-08 | プレス成形シミュレーション方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2022121024A JP2022121024A (ja) | 2022-08-19 |
| JP7541241B2 true JP7541241B2 (ja) | 2024-08-28 |
Family
ID=82849566
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021018149A Active JP7541241B2 (ja) | 2021-02-08 | 2021-02-08 | プレス成形シミュレーション方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7541241B2 (ja) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005138120A (ja) | 2003-11-04 | 2005-06-02 | Toyota Motor Corp | プレス成形中の金型たわみ分布のシミュレーション方法 |
| JP2008263788A (ja) | 2007-04-16 | 2008-11-06 | Sasaki Corporation | 畦成形機 |
| JP4590377B2 (ja) | 2006-07-14 | 2010-12-01 | 関東自動車工業株式会社 | 樹脂成形品の剛性最適化解析プログラム及びその記録媒体 |
| JP6380536B2 (ja) | 2014-07-30 | 2018-08-29 | 新日鐵住金株式会社 | モデル設定方法、成形シミュレーション方法、プログラム、及びプログラムを記録したコンピュータ読み取り可能な記録媒体 |
| JP2018156127A (ja) | 2017-03-15 | 2018-10-04 | 株式会社Ｊｓｏｌ | プレス曲げ成形解析モデル生成システム及びプログラム |
| JP6760546B1 (ja) | 2019-01-25 | 2020-09-23 | 日本製鉄株式会社 | 剛性影響度分析方法、剛性影響度分析装置およびプログラム |
| JP6774583B1 (ja) | 2020-03-26 | 2020-10-28 | 株式会社Ｊｓｏｌ | 高精度成形解析用メッシュ作成システム、高精度成形解析用メッシュ作成方法、および高精度成形解析用メッシュ作成プログラム |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102395973B (zh) * | 2009-04-20 | 2013-08-28 | 新日铁住金株式会社 | 成形模拟方法、成形模拟装置、成形模拟程序及其记录介质 |
| DE102017215395B4 (de) * | 2017-09-04 | 2022-12-15 | Bayerische Motoren Werke Aktiengesellschaft | Verfahren zum Betreiben einer Umformpresse |
-
2021
- 2021-02-08 JP JP2021018149A patent/JP7541241B2/ja active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005138120A (ja) | 2003-11-04 | 2005-06-02 | Toyota Motor Corp | プレス成形中の金型たわみ分布のシミュレーション方法 |
| JP4590377B2 (ja) | 2006-07-14 | 2010-12-01 | 関東自動車工業株式会社 | 樹脂成形品の剛性最適化解析プログラム及びその記録媒体 |
| JP2008263788A (ja) | 2007-04-16 | 2008-11-06 | Sasaki Corporation | 畦成形機 |
| JP6380536B2 (ja) | 2014-07-30 | 2018-08-29 | 新日鐵住金株式会社 | モデル設定方法、成形シミュレーション方法、プログラム、及びプログラムを記録したコンピュータ読み取り可能な記録媒体 |
| JP2018156127A (ja) | 2017-03-15 | 2018-10-04 | 株式会社Ｊｓｏｌ | プレス曲げ成形解析モデル生成システム及びプログラム |
| JP6760546B1 (ja) | 2019-01-25 | 2020-09-23 | 日本製鉄株式会社 | 剛性影響度分析方法、剛性影響度分析装置およびプログラム |
| JP6774583B1 (ja) | 2020-03-26 | 2020-10-28 | 株式会社Ｊｓｏｌ | 高精度成形解析用メッシュ作成システム、高精度成形解析用メッシュ作成方法、および高精度成形解析用メッシュ作成プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2022121024A (ja) | 2022-08-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6947809B2 (en) | Method of modifying stamping tools for spring back compensation based on tryout measurements | |
| KR101893312B1 (ko) | 모델 설정 방법, 성형 시뮬레이션 방법, 성형용 공구의 제조 방법, 프로그램, 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체 및 유한 요소 모델 | |
| US20100036646A1 (en) | Analytical model preparation method, and simulation system method for predicting molding failure | |
| JP2010009574A (ja) | 金型設計装置およびその方法 | |
| EP3751344A1 (en) | Simulation method, simulation apparatus, and program | |
| JP2013071120A (ja) | プレス成形におけるスプリングバック要因分析方法及び装置 | |
| JP5969886B2 (ja) | 対応点算出システム及びプログラム、並びに金型形状生成システム及びプログラム | |
| JP2024062444A (ja) | 金型形状データの作成方法、金型の製造方法、プレス成形品の製造方法、プログラム、及び金型形状データ作成システム | |
| JP7541241B2 (ja) | プレス成形シミュレーション方法 | |
| JP4221589B2 (ja) | プレス成形中の金型たわみ分布のシミュレーション方法 | |
| JP7347715B1 (ja) | 積層造形体の寸法品質評価方法及びその評価システム | |
| JP4305645B2 (ja) | 板成形のシミュレーション方法 | |
| JP6397149B1 (ja) | 金型たわみモデル作成システム、および金型たわみモデル作成プログラム | |
| JP5107595B2 (ja) | シミュレーション解析方法および金型の設計方法 | |
| JP7492139B2 (ja) | 金型設計システム、金型設計プログラム、金型設計方法、プレス装置及びプレス成形品の製造方法 | |
| JP4884840B2 (ja) | プレス金型修正形状データの作成方法 | |
| JP2010176573A (ja) | 金型設計装置およびその方法 | |
| CN120180802A (zh) | 一种用于冲压模具模面变形补偿值的获取方法 | |
| JP7590646B2 (ja) | プレス機の剛性予測方法 | |
| JP3966075B2 (ja) | ワーク形状評価装置およびその方法ならびにそのプログラム | |
| JP2002126834A (ja) | 金型およびその設計装置、設計方法ないし製造装置、並びに、成形品 | |
| JP6809941B2 (ja) | プレス曲げ成形解析モデル生成システム及びプログラム | |
| JP2024162294A (ja) | プレス成形シミュレーション方法 | |
| CN114282412A (zh) | 一种导轨滑块系统有限元分析的等效刚度参数设置方法 | |
| CN115438535A (zh) | 一种利用有限元综合分析结果驱动模具加工数据型面变形的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20231019 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20240716 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20240717 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240729 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7541241 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |