コンピュータを用いたプレス加工シミュレーションにおいて、3次元局所分岐理論により定めた分岐界面に対し、分岐界面に垂直な方向の伸び変形を含む変形つまり劈開変形に着目し、劈開変形の進展と劈開変形を起こすために必要な力の関係を求め、劈開変形を起こすために必要な力のうち分岐界面に垂直な方向に作用する力の増加率が0または減少になるつまり劈開変形の不安定条件を満たすことをもって破断と判定する、ことを特徴とする被加工材の破断判定方法。請求項1に記載の3次元局所分岐理論により定めた分岐界面の劈開変形の不安定条件を、分岐界面に垂直な方向の伸び変形を含む劈開変形を表す速度勾配テンソルを定め、この速度勾配テンソルを用いた公称応力速度テンソルから分岐界面に垂直な方向に作用する力の大きさの変化を求め、この力の大きさが増加しないことをもって破断と判定する、ことを特徴とする被加工材の破断判定方法。請求項1に記載の3次元局所分岐理論により定めた分岐界面の劈開変形の不安定条件を、分岐界面に垂直な方向の伸び変形を含む劈開変形を表す速度勾配テンソルを定め、この速度勾配テンソルを用いた公称応力速度テンソルから分岐界面に垂直な方向に作用する力の大きさの変化を求め、この力の大きさが増加しないことをもって破断と判定する、ことをコンピュータに実行させる事を特徴とする被加工材の破断判定プログラム。請求項3に記載の被加工材の破断判定プログラムを記録した事を特徴とするコンピュータ読み取り可能な記録媒体。本発明の課題は、前述した欠点を回避する、即ち、破断発生有無の判定を確実に行うことができるようにし、熟練者の経験や勘を加味して判断する必要のない、薄板だけでなく鍛造にも適用できる、破断判定方法を得ることである。本発明は上記課題を解決するために、三次元局所分岐理論で求めた、局所変形が起こる条件局所分岐を満たした後の変形つまり局所分岐後の変形に着目し、局所変形が起こる場所を示す分岐界面に対して垂直な方向の伸びを含む変形つまり劈開変形が起こると仮定し、この劈開変形を起こすために必要な力が増加せずに劈開変形が進展するという劈開変形の不安定条件をもって破断と判定する、ことを特徴とする被加工材の破断判定方法である。劈開変形を起こすために必要な力が増加せずに劈開変形の進展を破断と考える理由は、破断とは材料が分離する現象であることに着目すると、変形が分岐界面に集中したのち、分岐界面を引き剥がそうとする変形が起こり、この変形が大きく進展すると、最初はつながっていた材料が最終的には分離する、と考えられるためである。この劈開変形を進展させるとき、材料が変形に抵抗する力の増加により単位面積あたりの力はより多くの力を必要とする。しかし、分岐界面に垂直な方向の伸び変形が進展すると、分岐界面の面積がもとの面より小さくなるため、分岐界面全体の劈開変形に必要な力が減少に転じる場合がある。すると、劈開変形に必要な力の増加なしに劈開変形が無限に進展するため、面積が無限に小さくなりやがて0になるり、破断すると考えられる。そこで、この劈開変形に必要な力の増加率を求めることで、増加率が正の値なら劈開変形に必要な力は増加するため破断しないと判定し、増加率が0または負の値なら劈開変形に必要な力は変化しないまたは減少するため破断と判定する。具体的には三次元局所分岐理論によって分岐界面を求め、この分岐界面に垂直な方向の伸び変形を含む変形つまり劈開変形を表す速度勾配テンソルを定め、この速度勾配テンソルに従った変形の進展によって変化する分岐界面に作用する力のうち、分岐界面に垂直な方向の大きさの変化率を求め、その値が0または負の値になることをもって破断と判定する、ことを特徴とする被加工材の破断判定方法である。本発明の加工品の破断判定方法は、破断発生の判定を、ある瞬間の被加工材の状態応力テンソル、加工硬化率に基づいて行うため、背景技術の欠点である熟練者の経験や勘を加味して判断する必要のない、確実な破断の判定が可能となる。また、本発明の加工品の破断判定方法は、局所分岐後におこる劈開変形に着目した破断判定方法であるため、鍛造のような局所変形が起こる条件は満たすが破断が発生し難い加工でも、熟練者の経験や勘を加味して判断する必要のない、確実な破断の判定が可能となる。本発明は上記課題を解決するために、三次元局所分岐理論で提示されている式1によって分岐界面を求め、分岐界面に垂直な方向の伸び変形を含む変形つまり劈開変形の速度勾配テンソルを求め、速度勾配テンソルに従った変形の進展によって分岐界面に作用する力の分岐界面に垂直な方向の大きさの変化率を求め、その値が0または負の値になることをもって破断と判定する、ことを特徴とする加工品の破断判定方法である。以下、詳細に説明する。三次元局所分岐理論では、それまで一様に変形していた状態から、一部のみが変形して他に変形が広がらない局所変形へ突然切り替わる分岐現象が起こる条件分岐限界と呼ぶを、被加工材の応力σと加工硬化率hから判定出来る、ことを説明している。例えば、特許文献1参照。前記の一部は三次元局所分岐理論では、図1のハッチング面以下、分岐界面という。Cを指す。この分岐界面を引き剥がす方向は三次元局所分岐理論で提示されている式1の分岐界面の法線ベクトルnから知る事ができる。ここで、破断とは材料が分離する現象であることに着目すると、変形が分岐界面に集中したのち、分岐界面に作用する力に引き剥がす方向の力が含まれていると、分岐界面を引き剥がす方向の変形劈開変形が起こると考えられる。このとき劈開変形を起こすために必要な力が増加すると劈開変形が途中で停止し変形が止まり安定するが、劈開変形を起こすために必要な力が増加しないと劈開変形が停止せずに無限に進展するような劈開変形の不安定条件になり、最終的に材料は分離し破断が起こると考え、この劈開変形の不安定条件を用いて破断の判定を行う事に特徴がある。具体的な手順を示すと、まずはじめに分岐界面を得るため、三次元局所分岐理論で提示されている下記数式に、ある瞬間の被加工材の状態すなわち応力テンソルおよび加工硬化率を与え、分岐界面の法線ベクトルnおよび分岐界面に平行で且つ分岐変形の方向のベクトルmを求める。以下、上記式を数式1という。ここで、σ細字は、応力テンソルの大きさノルムhは、加工硬化率nは、分岐界面の法線ベクトルmは、分岐界面に平行な方向のベクトルで、nに直交するベクトルsは、応力テンソルの方向を表す単位テンソルHとΣは、変数に上記m、nを用い、パラメータにsが与えられる式minは、括弧内の最小値crは、括弧内によって定義される局所変形の条件を意味する表現である。ここで、HおよびΣの具体的な式は以下に示すようなものである。以下、上記式を数式2という。以下、上記式を数式3という。ここで、数式2と数式3のnは、分岐界面の法線ベクトルmは、分岐界面に平行な方向のベクトルで、nに直交するベクトルsは、応力テンソルの方向を表す単位テンソルIは、単位テンソルD傍線は、速度型の材料構成式の接線係数テンソルから加工硬化率hを分離したもので、具体的には下記数式に示す。以下、上記式を数式4という。ここで、Dは、速度型構成式における接線剛性テンソルhは、加工硬化率sは、応力テンソルの方向を表す単位テンソルを示している。このとき、数式1の右辺の値が最小値となるベクトルmおよびnの組を得る事ができる。このnが分岐界面の法線ベクトルを表し、mは分岐界面に平行で且つ分岐変形の方向のベクトルを表している。次いで、劈開変形の速度勾配テンソルLを作成する。図2を参照して説明すると、最も単純な劈開変形は、分岐界面である図2のハッチング面Cに繋がったA部とB部の平行を保ったまま分岐界面に垂直な方向図2のベクトルnで示す方向に引き剥がそうとする変形である。この変形では分岐界面Cの部分で、分岐界面に垂直な方向の伸び変形が起こる。ここで被加工材が変形する時に体積を一定に保つ条件を加味すると、ベクトルnに直交する2つの互いに直交するベクトル図2のベクトルmやベクトルn×mが示す方向では縮む変形が起こる。これらを考慮して、この劈開変形を表す速度勾配テンソルLは下記数式で表す事が出来る。以下、上記式を数式5という。ここで、nは、分岐界面の法線ベクトルmは、分岐界面に平行な方向のベクトルで、nに直交するベクトルn×mは、nとmに直交するベクトル○の中に×の記号は、テンソル積を表すμは、m方向とm×n方向の縮み変形の比率を表す任意の値である。前記の例では図2のA部とB部が平行に引き剥がされるとしたが、劈開変形は図2のA部とB部が傾きながら引き剥がしてもよく、数式の表現を若干変更することでどのような劈開変形の速度勾配テンソルLも作成できる。しかし、重要なのは分岐界面に垂直な方向の伸びを含む変形であることで、この条件を満たす変形の速度勾配テンソルであれば本発明の目的に利用できる。次いで、劈開変形の進展と劈開変形を起こすために必要な力の関係を求めるため、下記数式で示す公称応力速度テンソルSを求める数式に、定めた劈開変形の速度勾配テンソルLを与え、劈開変形の公称応力速度テンソルSを求める。以下、上記式を数式6という。ここで、σ太字は、加工材の応力テンソルDは、速度型構成式における接線剛性テンソルLは、作成した劈開変形を表す速度勾配テンソルεは、上記速度勾配テンソルLの対称成分ωは、上記速度勾配テンソルLの反対称成分trは、括弧内のテンソルの転置テンソルを意味する表現である。そして、得られた劈開変形の公称応力速度テンソルSと分岐界面の法線ベクトルnの内積を用いて、分岐界面に垂直な方向に作用する力の大きさの変化率Tの値を求める。具体的には下記数式で示す方法で変化率Tを計算できる。以下、上記式を数式7という。ここで、Sは、劈開変形の公称応力速度テンソルnは、分岐界面の法線ベクトルtrは、括弧内のテンソルの転置テンソルを意味する表現である。ここで、分岐界面に対して垂直な方向に作用する力の大きさを求めるのは、分岐界面を引き剥がす方向に作用する力だからであり、また、変化率を求めるのは、劈開変形の進行によって作用する力が増減のどちらの方向に変化するかに着目するからである。求めた変化率Tの値が正の値ならば、劈開変形が進行するためには今ある力では不足で増加が必要であり、途中で停止して安定することを示しているため、破断は起こらないと判定する。また求めた変化率Tの値が0または負の値になることは、劈開変形が進行するためには今ある力で十分または過剰であり、劈開変形が停止せずに無限に進展するような不安定な条件であることを示しており、材料は分離し破断が起こるため、破断と判定する。ところで、この劈開変形の不安定条件は、ある瞬間の応力テンソルσと加工硬化率hから求めるため、変形に至る過程に依存せず、前述の背景技術の項で述べたひずみ経路依存性に影響されない破断の判定を行うことができ、この点が本願発明の特徴である。また、局所変形をもって破断と判定する特許文献1の方法と異なり、局所変形後におこる劈開変形の不安定をもって破断と判定するので、鍛造のような局所変形が起こるが破断には至りにくい加工に対しても、より適切に破断の判定を行うことができ、この点が本願発明の特徴である。特開2008-105048号公報