(57)【要約】スポット溶接構造に対して有限要素法解析用シェルモデルを作成し、作成した有限要素法解析用シェルモデルを用いて有限要素法線形弾性解析を行ってスポット溶接部中央のナゲット部の分担荷重、そのナゲット部を中心に描いた直径（Ｄ）の円周上の変位とに基づいて前記ナゲット部における公称構造応力を弾性学の円板曲げ理論および２次元弾性論を用いて求め、該公称構造応力よりスポット溶接構造の疲労寿命を予測する溶接構造の疲労寿命予測方法。  本法によれば、スポット溶接構造の疲労寿命を簡易迅速正確に予測することができる。
2004099761.tif 000002
【特許請求の範囲】
【請求項１】
複数の板をあわせてスポット溶接構造を形成し、該スポット溶接構造に対して有限要素法解析用シェルモデルを作成し、作成した有限要素法解析用シェルモデルを用いて有限要素法線形弾性解析を行ってスポット溶接部中央のナゲット部のはく離荷重、曲げモーメント、せん断力及びねじりモーメントの分担荷重、そのナゲット部を中心に描いた直径（Ｄ）の円周上の変位を算出し、算出した該分担荷重、円周上の変位とに基づいて前記ナゲット部におけるはく離荷重、曲げモーメント、せん断力及びねじりモーメントの公称構造応力を弾性学の円板曲げ理論および２次元弾性論を用いて求め、該公称構造応力よりスポット溶接構造の疲労寿命を予測することを特徴とするスポット溶接構造の疲労寿命予測方法。
【請求項２】
前記複数の板が平板とＬ型板とからなり、前記有限要素法解析用シェルモデルがスポット溶接構造のナゲット部のＬ型板のフランジ幅を一辺とする正方形内を細かく分割し、放射方向にナゲット内を２分割、外側を４分割し、周方向に８分割したことを特徴とする請求項１に記載のスポット溶接構造の疲労寿命予測方法。
【請求項３】
前記スポット溶接構造の疲労寿命の予測は、予めスポット溶接構造品について引張せん断疲労試験、はく離疲労試験及び複合荷重疲労試験を行って、公称構造応力と破断繰り返し数との関係を表すマップを形成し、該マップを公称構造応力をもとに参照して破断繰り返し数を算出することにより行うことを特徴とする請求項１または請求項２に記載のスポット溶接構造の疲労寿命予測方法。
【発明の詳細な説明】
【技術分野】  本発明は、スポット溶接構造の疲労寿命予測方法に関し、更に詳細に説明すると、２枚以上の板をあわせてスポット溶接構造を形成し、該スポット溶接構造に対して有限要素法解析用シェルモデルを作成し、前記スポット溶接構造のナゲット部の分担荷重を用いて公称構造応力を求め、該公称構造応力よりスポット溶接構造の疲労寿命を予測するスポット溶接構造の疲労寿命予測方法に関する。
【背景技術】  自動車の開発において、軽量化、開発期間短縮、試作車両の削減の課題に応えるため、近年コンピュータを用いたＣＡＥによる各性能の予測が積極的に推進されている。車体の強度評価についても例外ではなく、その技術も加速的に進歩している。
自動車車体は薄板で構成され、その板は一般的にスポット溶接で結合される。車体に作用する荷重は、この溶接部を通じて各部材に伝達されるため、溶接部に応力が集中し強度上の弱点となる場合が多い。しかもその数は膨大で、更にきまざまな形態の複合荷重が作用することから、スポット溶接部の疲労寿命を簡便に、かつ精度良く予測できる手法の開発要求が高まっている。
この要求に対して、ディ．ラダイ（Ｄｉｅｔｅｒ  Ｒａｄａｊ；Ｄｅｓｉｇｎ  ａｎｄ  Ａｎａｌｙｓｉｓ  ｏｆ  Ｆａｔｉｇｕｅ  Ｒｅｓｉｓｔａｎｔ  Ｗｅｌｄｅｄ  Ｓｔｒｕｃｔｕｒｅｓ，Ａｂｉｎｇｔｏｎ  Ｐｕｂｌｉｓｈｉｎｇ，１９９０  ３７８ｐ）らは自動車車体をスポット溶接構造から切り出したナゲットを中心とした直径（Ｄ）の円板を想定してＦＥＭシェル解析より得られたナゲット部の分担荷重とデータベースに格納した直径Ｄの値とを用いて公称構造応力σｎｓを求め、これを用いてスポット溶接部の疲労寿命を予測する手法を提案している。自動車技術会疲労信頼性部門委員会においても、新たにねじり成分による影響を取り入れて同様な手法を提案している。ここで公称構造応力σｎｓとはスポット溶接部（ナゲット部）端に生する最大主応力である。
これらの手法は、図２２に示す如く、ステップＳ１でスポット溶接構造を形成し、ステップＳ２で有限要素法解析用シェルモデル（ＦＥＭモデル）を形成したＦＥＭモデルに荷重データＤ１を与えることにより、ステップＳ３でスポット溶接部（ナゲット部）の６分力を求め、次いで、ステップＳ４でＤ値を格納したデータベースＤ２を参照して板理論より公称構造応力σｎｓを求め、ステップＳ５で公称構造応力σｎｓ−破断繰り返し数Ｎの関係を表すデータベースＤ３からなるマップを公称構造応力σｎｓをもとに参照して疲労寿命を予測するものである。
これら従来の手法では、分担荷重のうちはく離荷重と曲げモーメントに対しては、ナゲット部を剛体とし、これを中心とした外周直径（Ｄ）の円板を考え、弾性学の円板曲げ理論を適用して公称構造応力を求めている。この円板外周条件は変位と放射方向傾斜をゼロとする固定条件、すなわち全自由度拘束条件を用いている。これを実際のスポット溶嬢構造に適用して公称構造応力を求めようとすると、円板直径（Ｄ）をいくらにしたらよいかという非常に難しい問題が生じる。現状では図２２にフローチャートで示す如く、Ｄ値データベースＤ２を作成することで対処している。このＤ値データベースＤ２の作成には非常に労力を要する。
図２３には、スポット溶接構造から切り出されたナゲット部（１）（Ｎｕｇｇｅｔ，スポット溶接部）を中心とした直径（Ｄ）の円板（２）が示されており、ナゲット部（１）には、一般に、
▲１▼はく離荷重      Ｆｚ
▲２▼曲げモーメント    Ｍｘ，  Ｍｙ
▲３▼せん断力          Ｆｘ，  Ｆｙ
▲４▼ねじりモーメント  Ｍｚ
の６分力が作用する。
ディ．ラダイらは，ナゲット分担荷重のうち曲げモーメントとはく離荷重に対して、疲労強度パラメータである公称構造応力σｎｓを次のように求めている。
【数１】
2004099761.tif 000003
【数２】
2004099761.tif 000004
ここで、Ｄ，ｄ，ｔは図２３のナゲット部（１）を中心に持つ円板（２）の外周直径、ナゲット径、板厚である。ＦｚとＭはナゲット部（１）に作用するはく離荷重、曲げモーメントである。
〔数１〕、〔数２〕は図２３の円板２でナゲット部（１）は剛体とし、円板（２）外周では変位完全拘束として導かれている。
従って、〔数１〕、〔数２〕を用いて、実際のスポット溶接構造の公称構造応力σｎｓを求めようとするとき、円板外周の変位完全拘束条件が満たされる場合にはＤ値を必然的に決めることができ、公称構造応力σｎｓを精度良く推定することができる。しかし一般には実構造のナゲット部（１）周辺の円周（Ｄ）上の変位が完全に拘束されている場合は少なく、このような場合上式に含まれる円板（２）の外周直径（Ｄ）をいくらにすればよいのかという問題が生ずる。
せん断荷重Ｆｘ，Ｆｙに対する応力は、ナゲットに相当する直径ｄの円形剛体介在物を有する無限板に、その中心にｘ軸方向のせん断荷重Ｆｘが作用する問題として取り扱い、剛体介在物を除く領域におけるｘ軸上の応力成分σｘは、
【数３】
2004099761.tif 000005
として与えられる。νはポアソン比である。公称構造応力はナゲット端の最大主応力であるとして、この場合、
【数４】
2004099761.tif 000006
となる。
図２４は引張せん断を受ける等厚三枚重ね一点スポット溶接継手の中央板外表面の負荷方向垂直応力分布である。σ０は一様負荷引張応力である。式（３）による理論解はＦＥＭ三次元弾性解析の結果と比較すると、公称構造応力に相当するナゲット端の応力は近い値を示すが、ナゲット端から離れるにつれて両者の差異が大きくなる。これは、ナゲット端から遠方で零に漸近する式（３）では負荷応力σ０がナゲット近傍に及ぼす影響が計算上現れてこないためである。
例えば、内外圧を受ける中空円板の応力分布を考えてみる。この問題の厳密解は式（５）となる。
【数５】
2004099761.tif 000007
これを無限平板に内径と同径の穴がある問題の解である式（６）
【数６】
2004099761.tif 000008
と比較したのが図２６である。外圧が内圧に比べて小さいときには無限平板の解は厳密解と良く一致するが、内圧と同程度の大きさになると、両者の解には大きな差異が見られる。これと同じ現象が三枚重ね継手の応力分布である図２４に現れていると考えられる。したがってこのような場合無限平板の解では精度の良い応力を得ることが難しい。
図２７はブラケット付きスポット溶接継手を想定した大小２枚の平板を一点スポット溶接した継手に、その両端に一様引張応力σ０が負荷したときの負荷方向垂直応力分布である。負荷の作用している平板はそれ自身で平衡状態にあるので、これをシェル要素でモデル化し有限要素法解析を行うと、ナゲットの分担荷重は零となる。したがって式（３）ではこの場合の応力を計算することができない。
本発明の目的は、Ｄ値決定問題を解消するために、円板のたわみ、放射方向の傾斜、曲げモーメント、はく離荷重、せん断力及びねじりモーメントに基づいて公称構造応力を求め、この公称構造応力からスポット溶接構造の疲労寿命を簡易迅速に予測することができる経済性に優れたスポット溶接構造の疲労寿命予測方法を提案するものである。
また、従来の一様負荷引張応力の値がナゲット端から離れるにつれて誤差が多くなるという問題を解決する新しい方法を提供することにある。
【発明の開示】
上記目的を達成するために、本発明の請求項１に係るスポット溶接構造の疲労寿命予測方法は、複数、例えば２枚の板をあわせてスポット溶接構造を形成し、該スポット溶接構造に対して有限要素法解析用シェルモデルを作成し、作成した有限要素法解析用シェルモデルを用いて有限要素法線形弾性解析を行ってスポット溶接部中央のナゲット部の分担荷重、そのナゲット部を中心に描いた直径（Ｄ）の円周上の変位とを算出し、算出した分担荷重、円周上の変位とに基づいて前記ナゲット部における公称構造応力を弾性学の円板曲げ理論および２次元弾性論を用いて求め、該公称構造応力よりスポット溶接構造の疲労寿命を予測することを特徴とする。
この請求項１に係る発明では、スポット溶接構造から切り出したナゲット部を中心とする直径Ｄの円板を想定し、この円板外周部にスポット溶接構造における変位を与え、その円板中央にあるナゲット部にははく離荷重及び曲げモーメント、せん断力及びねじりモーメントを与えることによりこれらに基づく公称構造応力を求めることができるので、Ｄ値データベースを作成する必要がなく、Ｄ値決定問題を解消することができ、スポット溶接構造の疲労寿命を簡易迅速に予測することができる。
また、本発明の請求項２に係るスポット溶接構造の疲労寿命予測方法は、請求項１に係る発明において、前記２枚の板が平板とＬ型板とからなり、前記有限要素法解析用シェルモデルがスポット溶接構造のナゲット部のＬ型板のフランジ幅を一辺とする正方形内を細かく分割し、放射方向にナゲット内を２分割、外側を４分割し、周方向に８分割したことを特徴とする。
更に、本発明の請求項３に係るスポット溶接構造の疲労寿命予測方法は、請求項１または請求項２に係る発明において、前記スポット溶接構造の疲労寿命の予測は、予めスポット溶接構造品について引張せん断疲労試験、はく離疲労試験及び複合荷重疲労試験を行って、公称構造応力と破断繰り返し数との関係を表すマップを形成し、該マップを公称構造応力をもとに参照して破断繰り返し数を算出することにより行うことを特徴とする。
【図面の簡単な説明】  第１図は、本発明に係わるスポット溶接構造の疲労寿命予測方法のＬＰモデルを示す斜視図であり、第２図は、本発明に係わるスポット溶接構造の疲労寿命予測方法の有限要素法解析用シェルモデルのナゲット周辺の拡大斜視図であり、第３図は、本発明に係わるスポット溶接構造の疲労寿命予測方法の直径Ｄの円板の説明図であり、第４図は、本発明に係わるスポット溶接構造の疲労寿命予測方法のフローチャートであり、第５図は、本発明に係わるはく離荷重及び曲げモーメントが個々に作用したときのたわみ分布の厳密解との比較を行った図でああり、第６図は、本発明に係わるスポット溶接構造の疲労革命予測方法によって求めたナゲット中心を通るＸ軸上の応力σｒ分布をＦＥＭソリッド解と比較したもので、図１でθｘ，θｙ方向に負荷した場合を夫々示し、（ａ）はＬＰ＃９０＃９０、（ｂ）はＬＰ＃４５＃９０の応力σｒ分布図であり、第７図は、本発明に係わるスポット溶接構造の疲労寿命予湘方法によって求めた中央ナゲットの公称構造応力を示すグラフ、第８図は、本発明における等厚な３枚の短冊形平板を一点でスポット溶接して継手に引張せん断が作用する場合を取り上げた図であり、第９図は、本発明における三枚重ね継手のＦＥＭシェル解析モデルを示した図であり、第１０図は、本発明に係る等厚三枚重ね一点スポット溶接継手の応力分布を示した図であり、第１１図は、本発明に係るシェル要素で解析モデルを作成した図であり、第１２図は、本発明に係るブラケット付き平板の応力分布を示した図であり、第１３図は、本発明に係わるＦＥＭシェル解析モデルとしてＬ形板の一部を示した図であり、第１４図は、本発明に係わるシェル解析モデルを示した図の一例であり、第１５図は、本発明に係わる解の精度を検証するために行った３次元弾性解析用ＦＥＭソリッドモデルで、Ｌ形板の一部を示した図であり、第１６図は、本発明に係わる解析結果を示した図であり、第１７図は、本発明に係わる公称構造応力を示した図であり、第１８図は、本発明に係わる種々の複合荷重疲労試験が可能となるＤＣ試験片を示した図であり、第１９図は、本発明に係わる代表的なスポット溶接疲労試験片である引張せん断（ＴＳ）と十字形引張試験片（ＣＴ）を示した図であり、第２０図は、本発明に係わる各試験片に対して本手法で公称構造応力を求め、疲労試験データを整理した図であり、第２１図は、本発明に係わるスポット溶接構造の疲労寿命予測法をフローチャートで示した図であり、第２２図は、従来のスポット溶接構造の疲労寿命予測方法のフローチャートであり、第２３図は、スポット溶接構造の疲労寿命予測方法のナゲットに作用する６分力を示す説明図であり、第２４図は、本発明に係る等厚三枚重ね一点スポット溶接継手の応力分布を示した図であり、第２５図は、本発明に係わる内外圧を受ける中空円板の応力分布を示した図であり、第２６図は、本発明に係わる無限平板に内径と同径の穴がある問題の解と比較した図であり、第２７図は、本発明に係るブラケット付き平板の応力分布を示した図である。
【発明を実施するための最良の形態】  以下、本発明のスポット溶接構造の疲労寿命予測方法の実施の形態を図面に基づいて説明する。
図１に示すような簡単なスポット溶接構造であるＬＰモデル（１１）を用いて解説する。このＬＰモデル（１１）は平板（１２）とＬ型板（１３）とをあわせて三点でスポット溶接した構造である。平板（１２）の一対の辺（１２ａ）、（１２ａ）を変位完全拘束し、Ｌ型板（１３）の垂直片（１３ａ）の上端に荷重（Ｐ）を種々の方向に加え、中央のナゲット部（１５）に注目しこの公称構造応力σｎｓを求めることを考える。
先ず、有限要素法解析用シェルモデル（２１）を作成する。ＬＰモデル（１１）に対して有限要素法解析を実施するために通常自動車のボディ解析に用いられるシェル要素を用いて解析モデル（２１）を作成する。図２にＬＰシェル解析モデルの中央のナゲット部（１５）周りの拡大図を示す。一般部は粗い格子状に要素分割をし、ナゲット部（１５）付近はフランジ幅（Ｌ１）を一辺とする正方形内を細かく分割した。図２に示す実施の形態では放射方向にナゲット内を２分割し、外側を４分割し、周方向を８分割した。
前記有限要素法解析用シェルモデル（２１）を用いて有限要素法線形弾性解析を実施し、中央のナゲット部（１５）に作用する分担力を算出する。更に、ナゲット部（１５）中心から放射線上（図２の場合４本）にある節点のたわみ（ｚ方向の変位）を用いて、八角形の頂点（８点）でのたわみと放射方向の傾斜を求める。このたわみと放射方向の傾斜を求めた八角形は、図１の破線で示す直径（Ｄ）の円に内接し、このたわみと傾きは次に示す直径（Ｄ）の円板の外周部の支持条件として用いられる。
図１に破線で示すように、スポット溶接構造であるＬＰモデル（１１）からナゲット部（１５）を中心とする直径Ｄの円板（３１）を切り出して考える。この直径（Ｄ）の円板（３１）を図３に示す。円板（３１）の外周は前記八角形に外接する。この円板（３１）の応力と変形を、弾性学における円板曲げ理論を用いて求める。この際、この円板（３１）に作用する外力が前記ナゲット部（１５）に作用する分担力である。この場合ｚ方向のはく離荷重Ｆｚと曲げモーメントＭ（ＭｘとＭｙの合モーメント）である。
円板（３１）外周の支持条件は前記八角形の頂点のたわみと傾斜である。ただし、これらは離散値であるので、その周方向分布をフーリエ級数で近似的に表し、これを境界条件、即ち外周支持条件として与える。
ディ．ラダイらが提案した公称構造応力式〔数１〕、〔数２〕は、図３に示す円板（３１）外周支持条件が変位完全固定のもとに導かれたのものであるので、これを用いて図１の中央ナゲット部（１５）の公称構造応力を求めようとするとき、〔数１〕、〔数２〕に含まれるＤ値をいま図１の破線で示す円の外径（Ｄ）と同じにすると、図１のスポット溶接構造における直径（Ｄ）の円周上の変位は完全固定でなく、ディ．ラダイらの円板外周支持条件（変位完全固定）と異なることから、ディ・ラダイらの式〔数１〕、〔数２〕から図１の中央ナゲット部（１５）の公称構造応力σｎｓが精度良く得られない。これを〔数１〕、〔数２〕のＤ値を補正することで対処している。この補正はスポット溶接構造や荷重条件などによって異なるため、いろいろな場合についてＤ値をいくらにしたら良いかを別途求め、これをＤ値データベースＤ２として公称構造応力σｎｓを求めている。
これに対して、本発明のスポット溶接構造の疲労寿命予測方法では図３の円板外周支持条件として、スポット溶接構造における直径（Ｄ）の円の外周上の実際の変位を与えているため、より精度の高い公称構造応力σｎｓが得られる。また、この方法では前記Ｄ値補正が不要で、Ｄ値をいくらにしたら良いかという非常に困難な問題も解消され、したがってＤ値データベースＤ２を作成する必要がなくなる。
ナゲット分担荷重のうち，せん断成分（Ｆｘ，Ｆｙ）とねじり成分、また、６分力が同時に作用する複合荷重に対する公称構造応力は、後述する理論式より求められる。簡単なスポット溶接継手試験片を用いて複合荷重下疲労試験を行い作成した、公称構造応力σｎと破断繰り返し数Ｎｆデータベースから、前述した方法で得られた公称構造応力σｎｓに対応する破断寿命Ｎｆを予測することができる。以上の流れをフローチャートにしたのが図４である。
図４に示す如く、ステップ（Ｓ１１）でスポット溶接構造を形成し、ステップ（Ｓ１２）で有限要素法解析用シェルモデル（ＦＥＭモデル）を形成し、ステップ（Ｓ１３）で荷重データＤ１をＦＥＭモデルに与えて有限要素法線形弾性解析を行うことにより、スポット溶接部（ナゲット部）の６分力と節点変位とを求める。次いで、ステップ（Ｓ１４）で、算出した６分力と、節点変位に基づいて所定の演算を行って公称構造応力σｎｓを算出し、次いでステップ（Ｓ１５）で算出した公称構造応力σｎｓをもとに予め実験によって作成した公称構造応力σｎｓに対応する破断寿命ＮｆのデータベースＤ３を参照して、破断寿命Ｎｆを算出し、これに基づいて疲労寿命を予測する。
図２１と比較すると公称構造応力を算出する際にＤ値データベースＤ２が不要となっていることが分かる。
前記のように、本発明においては、Ｄ値決定問題を解消するために、曲げに対する公称構造応力を求める新しい手法を提供する。
図３のように、外径（Ｄ）、厚さ（ｔ）の円板を考え、その中心にあるナゲット（直径ｄ）にはく離荷重Ｆｏ，曲げモーメントＭが作用するとする。曲げモーメントは図１のＭｘ，Ｍｙの合モーメントである。また、ナゲットは剛体とする。
この新たな方法は、図３に示す如く、直径（Ｄ）、厚さ（ｔ）の円板（３１）を考える。前述したように、この円板（３１）に作用する外力ははく離荷重Ｆｚ，曲げモーメントＭ（図１のＭｘ，Ｍｙの合モーメント）で、これらは円形のナゲット部（１５）（直径ｄ）に作用し、円板（３１）の外周の支持条件は実際のスポット溶接構造のものと合わせる。ここの解析ではナゲット部（１５）を剛体として取り扱う。円板外周の境界条件は式（１），（２）のように変位完全拘束とせず、後述のようにスポット溶接構造の注目するナゲット周辺の変形状態に合わせるようにする。さらにナゲット部（１５）の中心に原点を持つ極座標系（ｒθｚ）を設ける。
円板（３１）は外力によって曲げられて内部にひずみと応力を生ずるので、これを求める。まず、弾性論における次の円板（３１）の曲げ問題の支配方程式は次の、
【数７】
2004099761.tif 000009
で与えられ、この〔数７〕を満たすようなたわみｗ（ｚ方向の変位）の関数を見つけることが必要である。一般に〔数７〕を満足するたわみ関数は次の〔数８〕で示すように三角関数を用いたｒとθの関数として表すことができる。
【数８】
2004099761.tif 000010
〔数８〕に含まれるＡｎ〜Ｄｎは境界条件により定まる未知係数で４（ｎ＋１）個ある。この係数が確定ればたわみ関数〔数８〕が確定したことになる。
前記係数を決める条件は次の２種類である。
▲１▼ナゲットが剛体である条件
図３の円形なナゲット部（１５）が変形せず剛体変位をすることから、円板のたわみｗと半径方向の傾斜∂ｗ／∂ｒはナゲット端ｒ＝ｄ／２で、〔数９〕となる。
【数９】
2004099761.tif 000011
Ｗｃとθｘｃ、θｙｃは後述の条件より定まる未定係数である。
▲２▼円板外周支持条件
円板外周では実際のスポット溶接構造の支持条件と合わせることが必要である。したがって、円板外周（ｒ＝Ｄ／２〉ではたわみＷｒ＝Ｄ／２と傾斜∂ｗ／∂ｒは〔数１０〕となる。
【数１０】
2004099761.tif 000012
式（１０）の右辺は図３の円板外周に課される変位であるが、これはスポット溶接構造をＦＥＭシェル解析し、得られた節点変位から与えられる。具体的には以下のようである。
（ｉ）注目するナゲットを中心に直径（Ｄ）の同心円を描き、ＦＥＭシェル解析により円周上の節点変位ｗ，θｒの離散値を求める。ｗは面外変位、θｒは半径方向のたわみｗの傾斜角
【数１１】
2004099761.tif 000013
である。
（ｉｉ）この離散値から変位ｗ，θｒの周方向分布を３次の周期スプライン補間関数で補間する。例えば補間関数を構成する節点ｑｉの数を９とすると、外円周上の変位ｗのスプライン補間関数は次式となる。
【数１２】
2004099761.tif 000014
ここでＢｉ，４は３次のＢスプラインで次の漸化式より得られる．
【数１３】
2004099761.tif 000015
式（１２）の係数αｉはθ＝ｑｉでｗ（θ）ｒ＝Ｄ／２がその節点値になるように決める。
変位θｒ（θ）ｒ＝Ｄ／２についても同様にしてスプライン補間関数で表す。
（ｉｉｉ）これらの補間式をフーリエ級数で表し、式（１０）右辺のように円板外周の変位境界条件として与える。式（１０）に含まれるフーリエ級数の係数は次式より計算される。
【数１４】
2004099761.tif 000016
境界条件を用いると，たわみ関数式（８）に含まれる未定定数はすべて決定する。
たわみ関数が決まると、これを用いて円板の断面力、断面モーメントは次式のように得られる。
【数１５】
2004099761.tif 000017
ここで、Ｄｐは板の曲げ剛性
【数１６】
2004099761.tif 000018
で、Ｅはヤング率、νはポアソン比である。
式（１５）の断面モーメントを使うと板の応力成分は次式で与えられる。
【数１７】
2004099761.tif 000019
さらに、図３の円板についてｚ軸方向の内外力の釣り合いと、ｘ、ｙ軸に関するモーメントの釣り合いを考えると、
【数１８】
2004099761.tif 000020
が得られる。境界条件によって決定した係数Ａｎ〜Ｄｎを式（８）のたわみ関数に代入した後、これを式（１８）に代入すると、式（９）に含まれる定数ｗｃとθｃｘ、θｃｙはナゲットに作用するはく離荷重Ｆｚと曲げモーメントＭｘ、Ｍｙで表すことができる。
したがって、スポット溶接構造のＦＥＭシェル解析を実施し、その結果、注目するナゲットの分担荷重とそれを中心とした直径（Ｄ）の円周上の節点変位が与えられると、式（８）の未定係数はすべて定まる。
本論文で提唱したはく離荷重と曲げモーメントに対する公称構造応力は、式（１）、（２）のように簡単な形に表わせないが、実構造に適用する際に円板外径Ｄをいくつにしたら良いかという問題がなくなる。さらにスポット溶接構造のナゲット周辺の変形状態を境界条件として与えることができるので、公称構造応力の精度向上を計ることができる。
本理論の解析プログラムは簡単であり、これをＣＡＥに組み込むことによりスポット溶接構造の寿命予測が容易になる。
次に、せん断荷重とねじりモーメントによる公称構造応力の理論解について、図２３に基づいてせん断力とねじりモーメントによる公称構造応力を求める方策について述べる。この問題は弾性論の平面応力問題として取り扱える。基礎方程式は、
【数１９】
2004099761.tif 000021
与えられ，これを満足する応力関数は一般に次のようになる。
【数２０】
2004099761.tif 000022
式（２０）に含まれる係数ａ０〜ｄ’ｎは未定係数で境界条件から決まる。応力成分は応力関数φを用いると、
【数２１】
2004099761.tif 000023
と表される。ひずみ成分と変位成分との関係は、
【数２２】
2004099761.tif 000024
である。ｕ，νはそれぞれγ，θ方向の変位である。平面応力問題における応力とひずみの関係は、
【数２３】
2004099761.tif 000025
となる。Ｅはヤング率、Ｇはせん断弾性係数である。
応力関数式（２０）を式（２１）に代入すると応力成分は次のように表せる。
【数２４】
2004099761.tif 000026
これらを式（２３）の第一、二式に代入してひずみεｘ，εｙを求め、このひずみを式（２２）の第一、二式に代入し積分をすることにより変位ｕ，νが得られる。
【数２５】
2004099761.tif 000027
ｆ（θ）はθのみの未知関数、ｇ（ｒ）はｒのみの未知関数である。変位式（２５）を式（２２）の第三式に代入してせん断ひずみを求め、式（２３）の第三式を用いるとせん断応力が得られる。これは式（２４）の第三式と等しくなければならないから、両式を比較することにより未知関数ｆ（θ）、ｇ（ｒ）は次のように求まる。
【数２６】
2004099761.tif 000028
Ｈ，Ｋ，Ｆは未知定数である。
境界条件は、ナゲット端（ｒ＝ｄ／２）で、
【数２７】
2004099761.tif 000029
である。式（２７）のｕｘｃ，ｕｙｃ，θｃは後述の条件より定まる定数で、ナゲットを剛体と仮定していることから、その外周上における変位ｕ、ｖの周方向分布は式（２７）のように与えられる。
円板外周（ｒ＝Ｄ／２）では、
【数２８】
2004099761.tif 000030
となる。
式（２７）右辺は円板外周に課さられる変位であるが、これはスポット溶接構造をＦＥＭシェル解析し、得られた節点変位から与えられる。具体的には以下のようである。
▲１▼注目するナゲットを中心に直径（Ｄ）の同心円を描き、ＦＥＭシェル解析により得られた直径（Ｄ）の円周上の節点変位ｕ，νの離散値を求める。
▲２▼この離散値から変位ｕ，νの周方向分布を３次の周期スプライン補間関数で補間する。例えば補間関数を構成する節点ｑｉの数を９とすると、外円周上の変位ｕのスプライン補間関数は次式となる。
【数２９】
2004099761.tif 000031
ここでＢｉ，４（θ）は３次のＢスプラインで次の漸化式より得られる。
【数３０】
2004099761.tif 000032
式（２９）の係数αｉはθ＝ｑｉでｕ（θ）ｒ＝Ｄ／２がその節点値になるように決める。
変位ν（θ）ｒ＝Ｄ／２についても同様にしてスプライン補間関数で表す。
（ｉｉｉ）これらの補間式をフーリエ級数で表し、境界条件として与える。
式（２７）のｕｘｃ，ｕｙｃ，θｃは、ナゲットに作用する分担荷重とナゲット端周辺の応力分布との釣合式を解くことによって得ることができる。そのつり合い式は、
【数３１】
2004099761.tif 000033
となる。
したがって、スポット溶接構造のＦＥＭシェル解析を実施し、その結果、注目するナゲット部の分担荷重と、ナゲット部を中心とした直径（Ｄ）の円周上の節点変位が与えられると、式（２０）の未定係数はすべて定まり、ナゲット周りの変位場と応力場が求まることになる。
次に、複合荷重下の公称構造応力について説明する。
図２２のように，はく離荷重、曲げモーメント、せん断荷重、ねじりモーメントが同時に作用するとき、次式で表されるナゲット端での最大主応力を公称構造応力とし、疲労強度パラメータとして用いる。
【数３２】
2004099761.tif 000034
ここで、σｒｓｕｍ，σｓｕｍ，τｒ  ｓｕｍはそれぞれの分担荷重の解析結果を重ね合わせたものである。
ここで本発明の有効性を確認するために、図１のＬＰモデル（１１）の中央ナゲット部（１５）の公称構造応力を求める。解の精度を検証するためにＬＰモデル（１１）の詳細な三次元ソリッドモデルを作成し有限要素法解析でも公称構造応力を求めた。通常自動車のボディの解析ではこのような詳細なソリッドモデルは使わず、計算負荷が軽いシェルモデルを用いる。
スポット溶接構造へ適用する前に本手法による解を検証するために、図３の円板で外周の変位を完全拘束とし、はく離荷重及び曲げモーメントが個々に作用したときのたわみ分布の厳密解との比較を行ったのが図５である。
本手法の解は、円板の外周条件（式（６）右辺）として、厳密解の円板中間円周上（はく離荷重：ｒ＝３．８２（ｄ／２）、曲げモーメント：ｒ＝３．８２（ｄ／２））のたわみと傾斜を用いたもので、厳密解とよく一致している。
図６（ａ），（ｂ）は本発明によって求めたナゲット部（１５）の中心を通るｘ軸上の応力σｒ分布をＦＥＭソリッド解と比較したものである。図１でθｘ，θｙ方向に負荷した場合をＬＰ＃θｘ  θｙで表している．図６（ａ）はＬＰ＃９０＃９０、図６（ｂ）はＬＰ＃４５＃９０であり、いずれの場合も両者は良く一致している。
図７は本発明で求めた中央ナゲットの公称構造応力を示す。いずれの荷重の場合も、本発明で求めた公称構造応力はソリッドモデルによるＦＥＭ解とよく一致している。以上から本発明によって寿命予測パラメータである公称構造応力が精度良く得られることが分かる。尚、本実施の形態では、スポット溶接構造にＬＰモデル１１を用いたが、これに限定されるものではなく、２枚以上の板のスポット溶接構造であれば、どのような構成であっても適用することができる。
〔例２〕
上記例１は分担荷重のうちはく離荷重または曲げモーメントに対する公称構造応力を求めるためのものであるので、この例２では分担荷重としてせん断力が支配的となる場合について考える。この典型的な例として、等厚な３枚の短冊形平板を一点でスポット溶接した継手に引張せん断が作用する場合を取り上げる（図８）。
境界条件として用いる節点変位データを得るために、有限要素法シェル解析を実施する。図９は三枚重ね継手のＦＥＭシェル解析モデルである。一般部は粗い格子状に要素分割をし、ナゲット付近はラップ長さを一辺とする正方形内を、放射方向にナゲット内を２分割、外側を４分割し、周方向に８分割した。さらにナゲット内では放射方向にシェル要素の辺に沿ってバー要素を設け、ナゲット相当の剛性を持つビーム要素を用いて上下の板をナゲット中心で結合した。ＦＥＭソルバーにＣＯＳＭＯＳ／Ｍを用いて、線形弾性解析を行った。
ここで、理論を適用してナゲット周辺の応力分布を求める手順について、図９のシェル解析モデルを用いて説明する。
▲１▼図９のようなシェルモデルを作成し、ＦＥＭ解析を実施する。
▲２▼▲１▼の解析から、ナゲット部の分担荷重と、直径（Ｄ）の同心円周上の８節点の半径方向および周方向変位ｕ，νを求める。
▲３▼▲１▼で求めた節点変位については変位ｕ，νの周方向分布を３次の周期スプライン補間関数で補間する。さらに、これをフーリエ級数で表し、これを図２３の円板外周の境界条件とする。
▲４▼▲３▼の境界条件の下で、上記の理論よりナゲット周辺の応力分布を求める。
以上の手順で求めた解析結果を図１０に示す。これは三枚重ね継手中央板外表面の中心軸に沿った軸応力分布で、ＦＥＭソリッド解析の結果と比較した。図２４に示したように従来の理論ではナゲットから離れるにつれて応力は減衰してゼロとなり、ＦＥＭソリッド解と合わなかったが、本手法ではＦＥＭソリッド解と良い一致を示している。
〔例３〕
図２７に示した大小２枚の平板を一点スポット溶接した継手に、その両端に一様引張応力σ０が負荷した場合について負荷方向垂直応力分を求める。まず、図１１に示すようにシェル要素で解析モデルを作成し、ナゲット周りの八角形の頂点に位置する節点変位を用いて上に述べた理論で応力分布を求めた結果を図１２である。本理論を用いると従来法では得られなかった応力分布が求められており、ＦＥＭソリッド解とも良く一致していることが分かる。
〔例４〕
スポット溶接構造例として図１３に示すＬＰモデルを取り上げる。平板とＬ形板とをあわせて三点でスポット溶接し、平板の一対の辺を全自由度拘束とし、Ｌ形板の上端に荷重Ｐを種々の方向に負荷する。荷重Ｐがθｘ，θｙ方向に作用する場合をＬＰ＿θｘ＿θｙと表すことにする。Ｌ形板のフランジ幅１５ｍｍ、長さ１３５ｍｍ、平板の幅は４５ｍｍである。スポット溶接ピッチ４５ｍｍで、板厚は両方とも０．８ｍｍである。
ＦＥＭシェル解析モデルとしてＬ形板の一部を図１４に示す。一般部は粗い格子状に要素分割をし、ナゲット付近はフランジ幅ｗ∫を一辺とする正方形内を細かく分割し、放射方向にナゲット内を２分割、外側を４分割し、周方向に４分割したＳｈｅｌｌ＿０４と８分割したＳｈｅｌｌ＿０８モデルである。さらにナゲット内では放射方向にシェル要素の辺に沿ってバー要素を設け、ナゲット相当の剛性を持つビーム要素を用いて上下の板をナゲット中心で結合した。
図１５は本手法の解の精度を検証するために行った３次元弾性解析用ＦＥＭソリッドモデルで、Ｌ形板の一部を示す。板厚方向に４分割し、ナゲット周辺は詳細に分割した。平板も同様である。荷重条件（ＬＰ＿θｘ＿θｙ）と変位拘束条件は図１３に示すとおりである。
上記理論を用いてナゲット周辺の応力分布を求める手順について、図１４（ａ）のシェル解析モデルＳｈｅｌｌ＿０４を用いて説明する。
▲１▼図１４（ａ）のようなシェルモデルを作成し、ＦＥＭ解析を実施する。
▲２▼▲１▼の解析から、ナゲット中心から放射線上（図１４（ａ）の場合ＡＣ，ＢＤの２本）にある節点のたわみ値を用いて、多角形と正方形とが接する節点（Ａ，Ｂ，Ｃ，Ｄの４点）におけるたわみとその放射方向の傾斜を求める。
▲３▼これら４点のたわみと傾斜から周方向の分布を補間し（式（１２））、これをフーリエ級数で表し、図３の円板の外周境界条件（式（１０））とする。
▲４▼▲１▼のＦＥＭシェル解析から注目するナゲットの分担力を求め、このうちはく離成分Ｆｘと曲げ成分ＭｘとＭｙを式（９）で用いる。
▲５▼▲３▼▲４▼で与えられた円板外周境界条件と荷重条件のもとで、円板外径Ｄ＝ｗ∫とし、上記理論よりはく離成分と曲げ成分に対してナゲット周辺の応力分布を求める。
▲６▼ナゲット分担荷重のせん断成分（Ｆｘ，Ｆｙ）とねじり成分Ｍｚに対しては曲げの場合と同じように▲２▼〜▲５▼の操作を行い、ナゲット周辺の応力分布を求める。
▲７▼▲５▼▲６▼の解析結果を重ね合わせることで複合荷重下の応力が得られ、上記公称構造応力（式（３２））を求める。
以上の手順で求めた解析結果を図１６に示す。これは荷重ＬＰ＿９０＿９０の場合で、Ｌ形板中央ナゲットのｘ軸上の曲げ応力σｒ分布で、要素分割の影響をみるためにシェル解析モデルＳｈｅｌｌ＿０４とＳｈｅｌｌ＿０８の結果を比較したものである。両者の結果には大きな差異は見られず、粗い要素分割であるＳｈｅｌｌ＿０４モデルでもＦＥＭソリッド解と良く一致し、良好な結果が得られている。以後Ｓｈｅｌｌ＿０４モデルを用いて解析した結果を示す。図中のσ０は荷重をＬ形板上面の作用面積で除した平均応力である。
図１７は本手法によって求めた公称構造応力をＦＥＭソリッド解析の結果と比較したもので、いずれの荷重の場合も良好な解が得られていることが分かる。
また、同一形状の試験片にて、種々の複合荷重疲労試験が可能となるＧｉｅｋｅ，Ｈａｈｎにより考案されたＤＣ試験片（単点スポット溶接した対向カップ形試験片）を図１８、代表的なスポット溶接疲労試験片である引張せん断（ＴＳ）と十字形引張試験片（ＣＴ）を図１９に示す。
これらの試験片に対して本手法で公称構造応力を求め、疲労試験データを整理したのが図２０で、狭いバンド幅内に整理出来ていることが分かる。
本手法を用いたスポット溶接構造の疲労寿命予測法をフローチャートで示すと図２１となる。
▲１▼スポット溶接構造に対してシェルモデルを作成し、有限要素法弾性解析を実施する。
▲２▼注目するナゲットの分担荷重とナゲットを中心とした同心円周上の節点変位値を求める。
▲３▼上記で述べた理論に基づき公称構造応力σｎｓを算出する。
別途用意したＤＣ試験片によるσｎｓ−Ｎ∫線図より疲労寿命を予測する。
【産業上の利用可能性】  以上説明したように、本発明の請求項１に係るスポット溶接構造の疲労寿命予測方法によれば、スポット溶接構造から切り出したナゲット部を中心とする直径Ｄの円板を想定し、この円板外周部にスポット溶接構造における変位を与え、その円板中央にあるナゲット部にははく離荷重、曲げモーメント、せん断力及びねじりモーメントの分担荷重を与えることにより公称構造応力を求めることができるので、Ｄ値データベースを作成する必要がなく、Ｄ値決定問題を解消することができ、スポット溶接構造の疲労寿命を簡易迅速に予測することができる。
また、本発明の請求項２に係るスポット溶接構造の疲労寿命予測方法によれば、複数枚、例えば前記２枚の板が平板とＬ型板とからなり、前記有限要素法解析用シェルモデルがスポット溶接構造のナゲット部のＬ型板のフランジ幅を一辺とする正方形内を細かく分割し、放射方向にナゲット内を２分割、外側を４分割し、周方向に８分割したので、簡易迅速にスポット溶接構造の疲労寿命を予測することができる。
更に、本発明の請求項３に係るスポット溶接構造の疲労寿命予測方法によれば、スポット溶接構造の疲労寿命の予測は、予めスポット溶接構造品について引張せん断疲労試験、はく離疲労試験及び複合荷重疲労試験を行って、公称構造応力と破断繰り返し数との関係を表すマップを形成し、該マップを公称構造応力をもとに参照して破断繰り返し数を算出することにより行うので、疲労寿命を正確に推定することができる。
【図１】
2004099761.tif 000035
【図２】
2004099761.tif 000036
【図３】
2004099761.tif 000037
【図４】
2004099761.tif 000038
【図５】
2004099761.tif 000039
【図６】
2004099761.tif 000040
【図７】
2004099761.tif 000041
【図８】
2004099761.tif 000042
【図９】
2004099761.tif 000043
【図１０】
2004099761.tif 000044
【図１１】
2004099761.tif 000045
【図１２】
2004099761.tif 000046
【図１３】
2004099761.tif 000047
【図１４】
2004099761.tif 000048
【図１５】
2004099761.tif 000049
【図１６】
2004099761.tif 000050
【図１７】
2004099761.tif 000051
【図１８】
2004099761.tif 000052
【図１９】
2004099761.tif 000053
【図２０】
2004099761.tif 000054
【図２１】
2004099761.tif 000055
【図２２】
2004099761.tif 000056
【図２３】
2004099761.tif 000057
【図２４】
2004099761.tif 000058
【図２５】
2004099761.tif 000059
【図２６】
2004099761.tif 000060
【図２７】
2004099761.tif 000061
【手続補正書】
【提出日】平成15年10月17日(2003.10.17)
【手続補正００１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項１
【補正方法】変更
【補正の内容】
【請求項１】複数の板をあわせてスポット溶接構造を形成し、該スポット溶接構造に対して有限要素法解析用シェルモデルを作成し、作成した有限要素法解析用シェルモデルを用いて有限要素法線形弾性解析を行ってスポット溶接部中央のナゲット部のはく離荷重、曲げモーメント、せん断力及びねじりモーメントの分担荷重、そのナゲット部を中心に描いた直径（Ｄ）の円周上の変位を算出し、算出した該分担荷重、円周上の変位とに基づいて前記ナゲット部におけるはく離荷重、曲げモーメント、せん断力及びねじりモーメントの公称構造応力を弾性学の円板曲げ理論および２次元弾性論を用いて、せん断力及びねじりモーメントによる公称構造応力を、
２次元弾性論の応力成分σｒ、σθ、τｒθの式
【数１】
2004099761.tif 000062
（ここで、ａｏ〜ｄ’ｎは荷重条件と境界条件より定まる未定係数）により求め、また、曲げモーメント及びはく離荷重による公称構造応力を、
曲げの応力成分σｒ、σθ、τｒθの式
【数２】
2004099761.tif 000063
ここで、
【数３】
2004099761.tif 000064
【数４】
2004099761.tif 000065
（ｒ，θは極座標、Ｅはヤング率、νはポアソン比、ｔは板厚、Ａｏ〜Ｄ’ｎは荷重条件と境界条件より定まる未定係数）により求め、該公称構造応力よりスポット溶接構造の疲労寿命を予測することを特徴とするスポット溶接構造の疲労寿命予測方法。


【手続補正００２】
【補正対象書類名】明細書
【補正対象項目名】００１６
【補正方法】変更
【補正の内容】
【００１６】
【数１６】
2004099761.tif 000066
で、Ｅはヤング率、νはポアソン比である。
式（１５）の断面モーメントを使うと板の応力成分は次式で与えられる。
【数１７】
2004099761.tif 000067
式（１７）に式（１５）を代入して整理すると（ｚ＝ｔ／２）、
【数１７−１】
2004099761.tif 000068
さらに、図３の円板についてｚ軸方向の内外力の釣り合いと、ｘ、ｙ軸に関するモーメントの釣り合いを考えると、
【数１８】
2004099761.tif 000069
が得られる。境界条件によって決定した係数Ａｎ〜Ｄｎを式（８）のたわみ関数に代入した後、これを式（１８）に代入すると、式（９）に含まれる定数ｗｃとθｃｘ、θｃｙはナゲットに作用するはく離荷重Ｆｚと曲げモーメントＭｘ、Ｍｙで表すことができる。
したがって、スポット溶接構造のＦＥＭシェル解析を実施し、その結果、注目するナゲットの分担荷重とそれを中心とした直径（Ｄ）の円周上の節点変位が与えられると、式（８）の未定係数はすべて定まる。
【国際調査報告】
2004099761.tif 000070
2004099761.tif 000071