"注目バー要素に最も近いバー要素までの距離を抽出して、スポット間距離L1の情報を取得する工程と、前記注目バー要素の端点に最も近い距離にある母材のエッジを検索し、エッジ間距離L2の情報を取得する工程と、前記注目バー要素の端点に最も近い距離にある母材の稜線を検索し、稜線間距離L3の情報を取得する工程と、前記スポット間距離L1、エッジ間距離L2を2倍した距離、稜線間距離L3を2倍した距離のうち、最も短い距離をスポット溶接の有効幅Bに決定する工程と、前記決定した有効幅Bを破断解析に用いる工程と、を含むことを特徴とするスポット溶接部の破断解析方法。注目シェル要素の線分に隣接するシェル要素が存在するか否かを検索する工程と、隣接するシェル要素が存在した場合、この注目シェル要素の法線ベクトルと隣接するシェル要素の法線ベクトルのなす角θが閾値以上であるか否かを判定する工程と、前記なす角θが閾値以上であった場合、隣接するシェル要素が存在する注目シェル要素の線分を、母材の稜線とする工程と、をさらに含み、共通の母材に属するすべてのシェル要素に対して、隣接するシェル要素が存在するか否かの検索を実行することを特徴とする請求項1に記載のスポット溶接部の破断解析方法。注目シェル要素の線分に隣接するシェル要素が存在するか否かを検索する工程と、隣接するシェル要素が存在しなかった場合、この注目シェル要素の線分を母材のエッジとする工程と、をさらに含み、共通の母材に属するすべてのシェル要素に対して、隣接するシェル要素が存在するか否かの検索を実行することを特徴とする請求項1又は2に記載のスポット溶接部の破断解析方法。前記なす角θの閾値が10度であることを特徴とする請求項2に記載のスポット溶接部の破断解析方法。注目バー要素に最も近いバー要素までの距離を抽出して、スポット間距離L1の情報を取得する工程と、前記注目バー要素の端点に最も近い距離にあるシェル要素のエッジを検索し、エッジ間距離L2の情報を取得する工程と、前記注目バー要素の端点に最も近い距離にあるシェル要素の稜線を検索し、稜線間距離L3の情報を取得する工程と、前記スポット間距離L1、エッジ間距離L2を2倍した距離、稜線間距離L3を2倍した距離のうち、最も短い距離を有効幅Bに決定する工程と、前記決定した有効幅Bを破断解析に用いる工程と、を含む破断解析処理をコンピューターに実行させることを特徴とするスポット溶接部の破断解析コンピュータープログラム。注目バー要素に最も近いバー要素までの距離を抽出して、スポット間距離L1の情報を取得する手段と、前記注目バー要素の端点に最も近い距離にあるシェル要素のエッジを検索し、エッジ間距離L2の情報を取得する手段と、前記注目バー要素の端点に最も近い距離にあるシェル要素の稜線を検索し、稜線間距離L3の情報を取得する手段と、前記スポット間距離L1、エッジ間距離L2を2倍した距離、稜線間距離L3を2倍した距離のうち、最も短い距離を有効幅Bに決定する手段と、前記決定した有効幅Bを破断解析に用いる手段と、を備えたことを特徴とするスポット溶接部の破断解析装置。しかしながら、特許文献1-3に開示されている破断解析方法では、スポット溶接部の性質や負荷状態に応じて、荷重型破断、モーメント型破断、ナゲット内破断のいずれの破断モードが発生するかを判別できないので、衝突シミュレーションにおける解析精度が必ずしも高いとは言えなかった。そのため、本願の出願人によって、荷重型破断、モーメント型破断、ナゲット内破断のいずれの破断モードが発生するかを判別可能な、改善された破断解析方法が検討されている特願2010-088271。この破断解析方法は、コンピューターを用い、スポット溶接される鋼板のそれぞれの板厚t、引張強さTS、伸びEl、化学成分、溶接部のナゲット径d、隣接する溶接部、エッジ又は稜線との距離で決まるスポット溶接の有効幅B、断面高さHに基づいて、荷重型破断、モーメント型破断、ナゲット内破断のそれぞれの破断モードの破断限界を求め、スポット溶接部の状態量が前記いずれかの破断モードの破断限界に達したときに該破断モードで破断したと評価する。上述の改善された破断解析方法は、衝突シミュレーションにおける解析精度が高いという長所がある反面、全てのスポット溶接部に対して有効幅Bを決定しなければならず、準備設定に時間がかかるという問題点があった。自動車の車体には、通常、数千点にも及ぶスポット溶接部が存在するため、これら全てのスポット溶接部に対して隣接する溶接部との距離,エッジとの距離,稜線との距離を算出し、その中から有効幅Bを決定し、コンピューターに入力するには長時間を要する。更に、スポット溶接部によっては、隣接する溶接部が存在しなかったり、稜線が存在しなかったりする等、必ずしも前述の3つの距離が揃わない場合もあり、数千点にも及ぶスポット溶接部のそれぞれについて隣接する溶接部との距離,エッジとの距離,稜線との距離のいずれを有効幅Bとすれば良いのかも明確でなかった。本発明は、前述のような問題点に鑑みてなされたものであり、スポット溶接部の破断解析を行うにおいて、多数存在するスポット溶接部のそれぞれについて、有効幅Bを短時間で決定することができるアルゴリズムを備えたスポット溶接部の破断解析方法、破断解析コンピュータープログラム、破断解析装置を提供することにある。また本発明の他の目的は、スポット溶接部の有効幅Bを適切に設定して、より精度の高い破断解析を行うことにある。即ち、本発明の要旨とするところは以下のとおりである。1注目バー要素に最も近いバー要素までの距離を抽出して、スポット間距離L1の情報を取得する工程と、前記注目バー要素の端点に最も近い距離にある母材のエッジを検索し、エッジ間距離L2の情報を取得する工程と、前記注目バー要素の端点に最も近い距離にある母材の稜線を検索し、稜線間距離L3の情報を取得する工程と、前記スポット間距離L1、エッジ間距離L2を2倍した距離、稜線間距離L3を2倍した距離のうち、最も短い距離をスポット溶接の有効幅Bに決定する工程と、前記決定した有効幅Bを破断解析に用いる工程と、を含むことを特徴とするスポット溶接部の破断解析方法。2注目シェル要素の線分に隣接するシェル要素が存在するか否かを検索する工程と、隣接するシェル要素が存在した場合、この注目シェル要素の法線ベクトルと隣接するシェル要素の法線ベクトルのなす角θが閾値以上であるか否かを判定する工程と、前記なす角θが閾値以上であった場合、隣接するシェル要素が存在する注目シェル要素の線分を、母材の稜線とする工程と、をさらに含み、共通の母材に属するすべてのシェル要素に対して、隣接するシェル要素が存在するか否かの検索を実行することを特徴とする前記1に記載のスポット溶接部の破断解析方法。3注目シェル要素の線分に隣接するシェル要素が存在するか否かを検索する工程と、隣接するシェル要素が存在しなかった場合、この注目シェル要素の線分を母材のエッジとする工程と、をさらに含み、共通の母材に属するすべてのシェル要素に対して、隣接するシェル要素が存在するか否かの検索を実行することを特徴とする前記1又は2に記載のスポット溶接部の破断解析方法。4前記なす角θの閾値が10度であることを特徴とする請求項2に記載のスポット溶接部の破断解析方法。5注目バー要素に最も近いバー要素までの距離を抽出して、スポット間距離L1の情報を取得する工程と、前記注目バー要素の端点に最も近い距離にあるシェル要素のエッジを検索し、エッジ間距離L2の情報を取得する工程と、前記注目バー要素の端点に最も近い距離にあるシェル要素の稜線を検索し、稜線間距離L3の情報を取得する工程と、前記スポット間距離L1、エッジ間距離L2を2倍した距離、稜線間距離L3を2倍した距離のうち、最も短い距離を有効幅Bに決定する工程と、前記決定した有効幅Bを破断解析に用いる工程と、を含む破断解析処理をコンピューターに実行させることを特徴とするスポット溶接部の破断解析コンピュータープログラム。6注目バー要素に最も近いバー要素までの距離を抽出して、スポット間距離L1の情報を取得する手段と、前記注目バー要素の端点に最も近い距離にあるシェル要素のエッジを検索し、エッジ間距離L2の情報を取得する手段と、前記注目バー要素の端点に最も近い距離にあるシェル要素の稜線を検索し、稜線間距離L3の情報を取得する手段と、前記スポット間距離L1、エッジ間距離L2を2倍した距離、稜線間距離L3を2倍した距離のうち、最も短い距離を有効幅Bに決定する手段と、前記決定した有効幅Bを破断解析に用いる手段と、を備えたことを特徴とするスポット溶接部の破断解析装置。本発明によれば、注目バー要素に最も近いバー要素までの距離を抽出して、スポット間距離L1の情報を取得し、注目バー要素の端点に最も近い距離にある母材のエッジを検索してエッジ間距離L2の情報を取得し、注目バー要素の端点に最も近い距離にある母材の稜線を検索して稜線間距離L3の情報を取得し、スポット間距離L1、エッジ間距離L2を2倍した距離、稜線間距離L3を2倍した距離のうち、最も短い距離を有効幅Bに決定するアルゴリズムを採用したことによって、多数存在するスポット溶接部のそれぞれについて、有効幅Bの決定を短時間で行うことが可能となる。さらに、本発明によれば、スポット間距離L1、エッジ間距離L2を2倍した距離、稜線間距離L3を2倍した距離のうち、最も短い距離を有効幅Bとすることにより、スポット溶接の破断に影響を及ぼす周辺要因のうち、最も強い影響を及ぼすものの存在を破断解析において考慮することができ、その結果、精度の高い破断解析を行うことが可能となる。以下、本発明の好ましい実施形態に従うスポット溶接部の破断解析方法について、添付図面を参照しながら詳しく説明する。但し、以下に説明する実施形態によって本発明の技術的範囲は何ら限定解釈されることはない。まず、スポット溶接部の破断解析モデルについて説明する。破断解析は、通常、有限要素法を用いてコンピューター上で行われる。そのため、スポット溶接部は、バー要素ビーム要素と称されることもある、シェル要素、ソリッド要素等を用いてモデル化され、ベクトル状態量として要素に負荷される荷重やモーメントの値を取得して破断解析を行う。例えば母材として鋼材Aと鋼材Bの2枚の板材を重ねて溶接したスポット溶接部の場合、図1に模式的に示すように、鋼材Aと鋼材Bの接続部位がバー要素aと端点A,Bで構成されているとしてモデル化する。一方、鋼材Aと鋼材Bは、シェル要素によってモデル化する。そして破断解析を実行する際には、コンピューターを用いてモデル化されたバー要素aと端点A,Bに働く応力やせん断力などを演算し、破断基準値を超えるか否かを判定する。そのための必要情報として、上述の改善された破断解析方法特願2010-088271では、スポット溶接される鋼板それぞれの板厚t、引張強さTS、伸びEl、化学成分、溶接部のナゲット径d、隣接する溶接部、エッジ又は稜線との距離で決まる有効幅B、断面高さHをコンピューターに入力している。一方、鋼材A、鋼材B及び鋼材Cの3枚の板材を重ねて溶接したスポット溶接部の場合、鋼材Aと鋼材B、鋼材Bと鋼材Cを、2枚重ねのときと同様にバー要素と端点でモデル化する。その結果、図1に模式的に示すように、鋼材Aと鋼材Bを接続するバー要素aと、鋼材Bと鋼材Cを接続するバー要素bを有し、鋼材B上にある端点Bをバー要素aとバー要素bが共有した構成にモデル化される。しかし、上述の改善された破断解析方法では、この3枚重ねのモデルをそのまま解析することができないため、見掛け上、2枚重ねのモデルに調整して解析を行う。具体的には、鋼材Aと鋼材Bを接続するバー要素aに注目した場合、鋼材Bと鋼材Cを1つの要素とみなし、鋼材Aの板厚の値と、鋼材Bと鋼材Cの板厚を足した値をそれぞれ入力する。一方、鋼材Cと鋼材Bを接続するバー要素bに注目した場合には、鋼材Aと鋼材Bを1つの要素とみなし、鋼材Cの板厚の値と、鋼材Aと鋼材Bの板厚を足した値をそれぞれ入力する。さらに、シェル要素によるモデル化について説明する。例えば解析対象物が自動車の車体の場合、図2に模式的に示すように、設計CADにおいて車体がシェル要素にモデル化され、モデル化された各シェル要素の情報がメモリ等の記憶部に格納されている。シェル要素の情報には、各シェルが属する母材の識別番号母材番号、及び、各シェルの節点の座標x,y,zが含まれる。また、各スポット溶接部の情報として、モデル化された各バー要素の情報がメモリ等の記憶部に格納されている。バー要素の情報には、バー要素が属する母材の識別番号母材番号、及び、各バー要素に付された識別番号要素番号に対応付けた端点の座標x,y,z、バー要素の代表点が含まれる。図3~図6は、モデル化された各要素情報を格納するデータベースの一例を示す。さらに、モデル化された対象物一部の模式図を図7に示す。まず、図3は、母材のデータベースである。このデータベースには、各母材に割り付けた識別番号母材番号、母材の断面情報番号の情報が含まれる。更に、各断面情報番号に対応つけて断面情報のデータベースが形成されている。断面情報には、母材の板厚tが含まれる。なお、本実施形態では特に使用しないが、板厚t以外の情報として、ELFORM:要素タイプ、SHRF:せん断面積係数、NIP:断面内板厚方向積分点数、PROPT:プリント出力オプション、QRIRID:要素積分方法、ICOMP:材料軸角度定義フラグ、SETYP:2次元ソリッド要素タイプ、NLOC:3次元シェル要素参照面位置、MAREA:非構造部の単位面積あたりの質量、IDOF:板厚フィールドの連続/不連続フラグ、EDGSET:節点セットの情報も、このデータベースに格納することができる。図4は、シェル要素のデータベースである。このデータベースには、各シェル要素に割り付けた識別番号要素番号、各シェル要素が属する母材番号、各シェル要素の節点の識別番号節点番号の情報が含まれる。図4は、母材を四角形のメッシュにしたシェル要素のデータベースの例であり、従って、一つのシェル要素が4つの節点n1~n4を有している。母材を三角形のメッシュにした場合には、一つのシェル要素が3つの節点n1~n3を有することになる。節点n1~n4の座標x,y,zの情報は、図5に示すように、節点番号に対応付けて、節点端点のデータベースに格納されている。図6は、バー要素のデータベースである。このデータベースには、各バー要素に割り付けた識別番号要素番号、各バー要素が属する母材番号、各バー要素の両端点n1,n2の識別番号端点番号の情報が含まれる。両端点n1,n2の座標x,y,zの情報は、図5に示すように、端点番号に対応付けて節点端点のデータベースに格納されている。バー要素のデータベースには、さらに代表点の情報が含まれている。バー要素の代表点は、例えばバー要素の長さ方向の中央点であり、両端点の座標x,y,zの情報から演算により算出される座標とすることができる。図6では、バー要素の識別番号要素番号に対応付けた代表点の識別番号代表点番号を格納しており、代表点番号に対応する座標x,y,zの情報は、節点端点のデータベースのような別のデータベース不図示に格納する。スポット溶接の有効幅Bとなり得る要素としては、図8a~cの3つがある。第1に、同一の母材上に複数のスポット溶接部が存在する場合には、図8aに模式的に示すように、注目しているスポット溶接部に最も近いスポット溶接部までの距離スポット間距離L1が有効幅Bの候補となり得る。シェル要素にモデル化した場合は、同一シェル要素上に存在する最も近いバー要素までの距離が、スポット間距離L1となる。第2に、図8bに模式的に示すように、注目しているスポット溶接部によって溶接された部品鋼材AとBの溶接部品の最も近いエッジまでの距離エッジ間距離L2を2倍にした値=L2×2が、有効幅Bの候補となり得る。スポット溶接部から最短のエッジ間距離L2を2倍にすることにより、スポット溶接部が部品上の偏った位置にあった場合にも、有効幅Bを適切な値に設定することができる。なお、エッジとは、1つのシェル要素だけに属した2つの節点で構成される線分を意味する。第3に、図8cに模式的に示すように、注目しているスポット溶接部と、このスポット溶接部が溶接している部品上で最も近い稜線までの距離稜線間距離L3を2倍にした値=L3×2が、有効幅Bの候補となり得る。ここでも、スポット溶接部から稜線間距離L3を2倍にすることにより、スポット溶接部が部品上の偏った位置にあった場合にも、有効幅Bを適切な値に設定することができる。なお、ここでいう稜線とは、2つのシェル要素に共通で、且つシェルの法線ベクトル同士のなす角θが閾値10度以上ある辺2つの節点で構成される線分を意味する。本実施形態に従う破断解析では、有効幅Bとなり得る上述の3つ候補値を抽出し、それらの候補の一つを有効幅Bに決定する。その手順について、図9のフローチャートを参照しながら説明する。有効幅Bの設定手順まず、図9のステップS100に示すように、注目シェル要素の線分について、隣接するシェル要素が存在するかを検索する。そして、図9のステップS101において隣接するシェル要素の存在の有無を判定する。隣接するシェル要素が存在しなかった場合、ステップS102において、この注目シェル要素の線分を母材のエッジとしてコンピューターのデータベースDBに格納する。前述の図3~図6のデータベースの例では、シェル要素のデータベースから情報を読み出し、注目したシェル要素が有する4つの節点番号のうち、節点番号が2つ同じシェル要素が有るか否かを判定する。節点番号が2つ同じシェル要素が有った場合には、隣接するシェル要素が存在すると判定する。節点番号が2つ同じシェル要素が見つからない場合は、隣接するシェル要素は存在しないと判定する。また節点番号が1つだけ同じ場合、両者は対角位置にあるので隣接するシェル要素は存在しないと判定する。一方、ステップS101の判定において隣接するシェル要素が存在していた場合には、図9のステップS103に示すように、注目シェル要素と隣接シェル要素の法線ベクトルが成す角度θを算出し、さらにステップS104において角度θが閾値10度以上であるか否かを判定する。判定の結果が閾値10度以上であった場合にのみ、ステップS105に示すように、注目シェル要素と隣接シェル要素をつないでいる線分を母材の稜線正確には母材の稜線の一部としてデータベースDBに格納する。なお、閾値は10度が好ましいが、必ずしも10度でなくともよく、10度~15度の範囲内で任意に設定することができる。10度未満の場合、母材の曲面を稜線と誤認する場合があり、抽出の精度が低下する場合がある。一方、15度を超えると稜線の抽出漏れがある場合がある。上述のステップS100~S105までの処理は、解析対象物の全てのシェルについて行い、図9のステップS106において全てのシェルについての処理が完了したと判定した場合に、次のステップS107に進行する。続いて、図9のステップS107に示すように、注目バー要素について隣接バー要素を検索し、共通する母材上に隣接バー要素が存在している場合には、その中の最短距離にあるバー要素までの距離を、スポット間距離L1とする。バー要素間の距離は、例えば各バー要素の、共通する母材上にある端点同士の距離である。この距離は、母材の面上に沿った行程距離とするのが望ましいが、3次元空間での絶対距離としてもよい。当該ステップにおいては、一つのバー要素についてスポット間距離L11の情報を取得すると、次の注目バー要素についてスポット間距離L12の情報を取得し、これを繰り返すことによって全てのバー要素についてスポット間距離L11,L12,,L1nの情報を取得する。次に、図9のステップS108に示すように、注目バー要素について、上述のデータベースDBに登録されているエッジまでの距離情報を読み出し、その中から最短のエッジまでの距離を、エッジ間距離L2とする。エッジ間距離L2は、例えば注目バー要素の端点とエッジと判定されたシェル要素の線分までの距離である。この距離は、シェル要素の面上に沿った行程距離とするのが望ましいが、3次元空間での絶対距離としてもよい。当該ステップにおいても、一つのバー要素についてエッジ間距離L21の情報を取得すると、次の注目バー要素についてエッジ間距離L22の情報を取得し、これを繰り返すことによって全てのバー要素についてエッジ間距離L21,L22,,L2nの情報を取得する。次に、図9のステップS109に示すように、注目バー要素について、上述のデータベースDBに登録されている稜線までの距離情報を読み出し、その中から最短の稜線までの距離を、稜線間距離L3とする。稜線間距離L3は、例えば注目バー要素の端点と稜線と判定されたシェル要素の線分までの距離である。この距離は、シェル要素の面上に沿った行程距離とするのが望ましいが、3次元空間での絶対距離としてもよい。当該ステップにおいても、一つのバー要素について稜線間距離L31の情報を取得すると、次の注目バー要素について稜線間距離L32の情報を取得し、これを繰り返すことによって全てのバー要素について稜線間距離L31,L32,,L3nの情報を取得する。次に、図9のステップS110に示すように、各バー要素について、ステップS107にて取得したスポット間距離L1の値と、ステップS108にて取得したエッジ間距離L2を2倍した値と、ステップS109にて取得した稜線間距離L3を2倍した値とを比較し、その中で最も短い値距離を選択する。そして、この選択した距離を、当該バー要素でモデル化されたスポット溶接部の有効幅Bに設定する。なお、スポット間距離L1の値、エッジ間距離L2を2倍した値、稜線間距離L3を2倍した値の全部又は2つが同じ値であっても、最も小さい値であればその値を有効幅Bに設定すればよい。続いて、破断解析の手順について説明する。以下の説明は、本願の出願人によって改善された破断解析方法特願2010-088271を一例に挙げているが、必ずしもこの破断解析方法に限定はされない。破断解析は、準備設定作業としてコンピューターに入力情報を読み込ませる。この入力情報には、図9の手順に従ってコンピューターが設定した有効幅Bの情報が含まれる。鋼板A,Bをスポット溶接する場合の入力項目を表1に示す。破断モードによって判定に用いられる入力項目は異なるが、ここに挙げたすべての入力項目を用いることで、すべての破断モードに対する評価を行うことができるため、最も早く限界値に達した破断モードを知ることが可能となる。次に、荷重型破断、モーメント型破断、ナゲット内破断のそれぞれの破断モードに応じて表1の○印に示す入力項目を使用して破断限界クライテリアを算出する。なお、それぞれの破断モードにおける破断限界クライテリアの算定方法は、特に限定されないが、例えば下記の方法を用いることが好ましい。まず、荷重型破断の場合は、スポット溶接部を有する試験片のせん断引張試験又は十字形引張試験を行い、ナゲット径dmmと前記試験片の幅Wmmとの比d/Wと、1式による応力集中係数αの関係を予め求め、任意の引張強さを有する材料を対象として2式によりせん断引張試験によるスポット溶接部の破断限界荷重FsNを算定する方法が好ましい。α=TSWt/F1ここで、TS:引張強さMPa、t:試験片の厚さmm、F:破断限界張力NFs=TSWt/α2また、モーメント型破断の場合は、スポット溶接部を有するフランジ引張試験を行い、スポット溶接部の端部に加えた曲げモーメントMNmと、試験材の板厚、板幅、強度特性から理論的に求まる全塑性モーメントMpNmから、3式によるモーメント効率γを予め求め、このモーメント効率γと、任意の板厚、板幅、強度特性を有する材料に対する全塑性モーメントMp ́から4式によるフランジ引張試験によるスポット溶接部の破断限界モーメントMlimNmを算定する方法が好ましい。γ=Mp/M3Mlim=Mp ́/γ4また、ナゲット内破断の場合は、例えば、下記5式により、スポット溶接部の破断限界荷重FsNを算定する方法が好ましい。Fs=e×Πd22×f×Ceqg5ここで、d:ナゲット径mm、Ceqナゲット部炭素当量の厚み方向の重み付き平均、efg:係数そして、各時刻ごとにスポット溶接要素の荷重モーメント出力に基づくモードごとの状態変数を計算する。次に、破断モードごとに前述の限界値と状態変数を比較する。そして、いずれかのモードの状態変数が限界値に達している場合、以後、破断発生済みと判定し、その後の溶接要素の相対変位ひずみに応じて許容荷重値を低下させる。全プロセスの力学計算終了後、破断詳細情報を出力する。上述の実施形態によれば、注目バー要素に最も近いバー要素までの距離を抽出して、スポット間距離L1の情報を取得し、注目バー要素の端点に最も近い距離にある母材のエッジを検索してエッジ間距離L2の情報を取得し、注目バー要素の端点に最も近い距離にある母材の稜線を検索して稜線間距離L3の情報を取得し、スポット間距離L1、エッジ間距離L2を2倍した距離、稜線間距離L3を2倍した距離のうち、最も短い距離を有効幅Bに決定するアルゴリズムを採用したことによって、多数存在するスポット溶接部のそれぞれについて、有効幅Bの決定を短時間で行うことが可能となる。さらに、上述の実施形態によれば、スポット間距離L1、エッジ間距離L2を2倍した距離、稜線間距離L3を2倍した距離のうち、最も短い距離を有効幅Bとすることにより、スポット溶接の破断に影響を及ぼす外乱要因のうち、最も強い影響を及ぼすものの存在を破断解析で考慮することができ、その結果、精度の高い破断解析を行うことが可能となる。最後に、上述の破断解析を実行可能なコンピューターシステムの一例について、図10のブロック図を参照しながら説明する。同図において、1200はコンピューターPCである。PC1200は、CPU1201を備え、ROM1202又はハードディスクHD1211に記憶された、或いはフレキシブルディスクドライブFD1212により供給されるデバイス制御ソフトウェアを実行し、システムバス1204に接続される各デバイスを総括的に制御する。前記PC1200のCPU1201、ROM1202又はハードディスクHD1211に記憶されたプログラムにより、本実施形態の処理を行う各機能手段が構成される。1203はRAMで、CPU1201の主メモリ、ワークエリア等として機能する。1205はキーボードコントローラKBCであり、キーボードKB1209から入力される信号をシステム本体内に入力する制御を行う。1206は表示コントローラCRTCであり、表示装置CRT1210上の表示制御を行う。1207はディスクコントローラDKCで、ブートプログラム起動プログラム:パソコンのハードやソフトの実行動作を開始するプログラム、複数のアプリケーション、編集ファイル、ユーザファイルそしてネットワーク管理プログラム等を記憶するハードディスクHD1211、及びフレキシブルディスクFD1212とのアクセスを制御する。1208はネットワークインターフェースカードNICで、LAN1220を介して、ネットワークプリンタ、他のネットワーク機器、或いは他のPCと双方向のデータのやり取りを行う。上述した実施形態の機能は、コンピューターがコンピュータープログラムを実行することによっても実現される。また、コンピュータープログラムをコンピューターに供給するための手段、例えばかかるプログラムを記録したCD-ROM等のコンピューター読み取り可能な記録媒体又はかかるプログラムを伝送するインターネット等の伝送媒体も本発明の実施形態として適用することができる。また、前記のプログラムを記録したコンピューター読み取り可能な記録媒体等のコンピュータプログラムプロダクトも本発明の実施形態として適用することができる。前記のコンピュータープログラム、記録媒体、伝送媒体及びコンピュータプログラムプロダクトは、本発明の範疇に含まれる。記録媒体としては、例えばフレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、CD-ROM、磁気テープ、不揮発性メモリ、ROM等を用いることができる。実施例1続いて、有効幅Bを自動算定して破断解析を行った実施例1について説明する。本実施例では、図11に示すように、2枚重ね溶接によって形成された構造部材の破断解析を行った。詳しい解析条件については同図に示す。破断解析を実行するにあたり、図9のフローチャートに従う処理を実行可能なプログラムがインストールされたコンピューターを用いて有効幅Bの自動算定を行い、各スポット溶接の有効幅Bを破断解析プログラムに認識させた。図12aは、破断個所を確認し易いように、図11aから図11bへの変形が20進行したときの解析結果を示している。同図に示される通り、本実施例では10点在るスポット溶接部のうち8点で破断が発生している。一方、すべてのスポット溶接の有効幅Bを破断解析プログラムのデフォルト値40mmとした場合、同図bに示されるように、実際には8点で破断するのに6点で破断する結果となっている。よって、図9に示すフローチャートに従って有効幅Bを自動算定すれば、作業時間を大幅に短縮できるのに加えて、衝突シミュレーションにおける解析精度を向上できる効果がある。そして、溶接部の数が数百点、数千点に増えるに伴いその効果がより顕在化する。"