JPWO2004099761A1 - スポット溶接構造の疲労寿命予測方法 - Google Patents
スポット溶接構造の疲労寿命予測方法 Download PDFInfo
- Publication number
- JPWO2004099761A1 JPWO2004099761A1 JP2004571574A JP2004571574A JPWO2004099761A1 JP WO2004099761 A1 JPWO2004099761 A1 JP WO2004099761A1 JP 2004571574 A JP2004571574 A JP 2004571574A JP 2004571574 A JP2004571574 A JP 2004571574A JP WO2004099761 A1 JPWO2004099761 A1 JP WO2004099761A1
- Authority
- JP
- Japan
- Prior art keywords
- nugget
- welded structure
- spot welded
- spot
- load
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N19/00—Investigating materials by mechanical methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0073—Fatigue
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0214—Calculations a priori without experimental data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0216—Finite elements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/24—Sheet material
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Immunology (AREA)
- Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Resistance Welding (AREA)
Abstract
Description
自動車車体は薄板で構成され、その板は一般的にスポット溶接で結合される。車体に作用する荷重は、この溶接部を通じて各部材に伝達されるため、溶接部に応力が集中し強度上の弱点となる場合が多い。しかもその数は膨大で、更にきまざまな形態の複合荷重が作用することから、スポット溶接部の疲労寿命を簡便に、かつ精度良く予測できる手法の開発要求が高まっている。
この要求に対して、ディ.ラダイ(Dieter Radaj;Design and Analysis of Fatigue Resistant Welded Structures,Abington Publishing,1990 378p)らは自動車車体をスポット溶接構造から切り出したナゲットを中心とした直径(D)の円板を想定してFEMシェル解析より得られたナゲット部の分担荷重とデータベースに格納した直径Dの値とを用いて公称構造応力σnsを求め、これを用いてスポット溶接部の疲労寿命を予測する手法を提案している。自動車技術会疲労信頼性部門委員会においても、新たにねじり成分による影響を取り入れて同様な手法を提案している。ここで公称構造応力σnsとはスポット溶接部(ナゲット部)端に生する最大主応力である。
これらの手法は、図22に示す如く、ステップS1でスポット溶接構造を形成し、ステップS2で有限要素法解析用シェルモデル(FEMモデル)を形成したFEMモデルに荷重データD1を与えることにより、ステップS3でスポット溶接部(ナゲット部)の6分力を求め、次いで、ステップS4でD値を格納したデータベースD2を参照して板理論より公称構造応力σnsを求め、ステップS5で公称構造応力σns−破断繰り返し数Nの関係を表すデータベースD3からなるマップを公称構造応力σnsをもとに参照して疲労寿命を予測するものである。
これら従来の手法では、分担荷重のうちはく離荷重と曲げモーメントに対しては、ナゲット部を剛体とし、これを中心とした外周直径(D)の円板を考え、弾性学の円板曲げ理論を適用して公称構造応力を求めている。この円板外周条件は変位と放射方向傾斜をゼロとする固定条件、すなわち全自由度拘束条件を用いている。これを実際のスポット溶嬢構造に適用して公称構造応力を求めようとすると、円板直径(D)をいくらにしたらよいかという非常に難しい問題が生じる。現状では図22にフローチャートで示す如く、D値データベースD2を作成することで対処している。このD値データベースD2の作成には非常に労力を要する。
図23には、スポット溶接構造から切り出されたナゲット部(1)(Nugget,スポット溶接部)を中心とした直径(D)の円板(2)が示されており、ナゲット部(1)には、一般に、
▲1▼はく離荷重 Fz
▲2▼曲げモーメント Mx, My
▲3▼せん断力 Fx, Fy
▲4▼ねじりモーメント Mz
の6分力が作用する。
ディ.ラダイらは,ナゲット分担荷重のうち曲げモーメントとはく離荷重に対して、疲労強度パラメータである公称構造応力σnsを次のように求めている。
〔数1〕、〔数2〕は図23の円板2でナゲット部(1)は剛体とし、円板(2)外周では変位完全拘束として導かれている。
従って、〔数1〕、〔数2〕を用いて、実際のスポット溶接構造の公称構造応力σnsを求めようとするとき、円板外周の変位完全拘束条件が満たされる場合にはD値を必然的に決めることができ、公称構造応力σnsを精度良く推定することができる。しかし一般には実構造のナゲット部(1)周辺の円周(D)上の変位が完全に拘束されている場合は少なく、このような場合上式に含まれる円板(2)の外周直径(D)をいくらにすればよいのかという問題が生ずる。
せん断荷重Fx,Fyに対する応力は、ナゲットに相当する直径dの円形剛体介在物を有する無限板に、その中心にx軸方向のせん断荷重Fxが作用する問題として取り扱い、剛体介在物を除く領域におけるx軸上の応力成分σxは、
図24は引張せん断を受ける等厚三枚重ね一点スポット溶接継手の中央板外表面の負荷方向垂直応力分布である。σ0は一様負荷引張応力である。式(3)による理論解はFEM三次元弾性解析の結果と比較すると、公称構造応力に相当するナゲット端の応力は近い値を示すが、ナゲット端から離れるにつれて両者の差異が大きくなる。これは、ナゲット端から遠方で零に漸近する式(3)では負荷応力σ0がナゲット近傍に及ぼす影響が計算上現れてこないためである。
例えば、内外圧を受ける中空円板の応力分布を考えてみる。この問題の厳密解は式(5)となる。
図27はブラケット付きスポット溶接継手を想定した大小2枚の平板を一点スポット溶接した継手に、その両端に一様引張応力σ0が負荷したときの負荷方向垂直応力分布である。負荷の作用している平板はそれ自身で平衡状態にあるので、これをシェル要素でモデル化し有限要素法解析を行うと、ナゲットの分担荷重は零となる。したがって式(3)ではこの場合の応力を計算することができない。
本発明の目的は、D値決定問題を解消するために、円板のたわみ、放射方向の傾斜、曲げモーメント、はく離荷重、せん断力及びねじりモーメントに基づいて公称構造応力を求め、この公称構造応力からスポット溶接構造の疲労寿命を簡易迅速に予測することができる経済性に優れたスポット溶接構造の疲労寿命予測方法を提案するものである。
また、従来の一様負荷引張応力の値がナゲット端から離れるにつれて誤差が多くなるという問題を解決する新しい方法を提供することにある。
この請求項1に係る発明では、スポット溶接構造から切り出したナゲット部を中心とする直径Dの円板を想定し、この円板外周部にスポット溶接構造における変位を与え、その円板中央にあるナゲット部にははく離荷重及び曲げモーメント、せん断力及びねじりモーメントを与えることによりこれらに基づく公称構造応力を求めることができるので、D値データベースを作成する必要がなく、D値決定問題を解消することができ、スポット溶接構造の疲労寿命を簡易迅速に予測することができる。
また、本発明の請求項2に係るスポット溶接構造の疲労寿命予測方法は、請求項1に係る発明において、前記2枚の板が平板とL型板とからなり、前記有限要素法解析用シェルモデルがスポット溶接構造のナゲット部のL型板のフランジ幅を一辺とする正方形内を細かく分割し、放射方向にナゲット内を2分割、外側を4分割し、周方向に8分割したことを特徴とする。
更に、本発明の請求項3に係るスポット溶接構造の疲労寿命予測方法は、請求項1または請求項2に係る発明において、前記スポット溶接構造の疲労寿命の予測は、予めスポット溶接構造品について引張せん断疲労試験、はく離疲労試験及び複合荷重疲労試験を行って、公称構造応力と破断繰り返し数との関係を表すマップを形成し、該マップを公称構造応力をもとに参照して破断繰り返し数を算出することにより行うことを特徴とする。
図1に示すような簡単なスポット溶接構造であるLPモデル(11)を用いて解説する。このLPモデル(11)は平板(12)とL型板(13)とをあわせて三点でスポット溶接した構造である。平板(12)の一対の辺(12a)、(12a)を変位完全拘束し、L型板(13)の垂直片(13a)の上端に荷重(P)を種々の方向に加え、中央のナゲット部(15)に注目しこの公称構造応力σnsを求めることを考える。
先ず、有限要素法解析用シェルモデル(21)を作成する。LPモデル(11)に対して有限要素法解析を実施するために通常自動車のボディ解析に用いられるシェル要素を用いて解析モデル(21)を作成する。図2にLPシェル解析モデルの中央のナゲット部(15)周りの拡大図を示す。一般部は粗い格子状に要素分割をし、ナゲット部(15)付近はフランジ幅(L1)を一辺とする正方形内を細かく分割した。図2に示す実施の形態では放射方向にナゲット内を2分割し、外側を4分割し、周方向を8分割した。
前記有限要素法解析用シェルモデル(21)を用いて有限要素法線形弾性解析を実施し、中央のナゲット部(15)に作用する分担力を算出する。更に、ナゲット部(15)中心から放射線上(図2の場合4本)にある節点のたわみ(z方向の変位)を用いて、八角形の頂点(8点)でのたわみと放射方向の傾斜を求める。このたわみと放射方向の傾斜を求めた八角形は、図1の破線で示す直径(D)の円に内接し、このたわみと傾きは次に示す直径(D)の円板の外周部の支持条件として用いられる。
図1に破線で示すように、スポット溶接構造であるLPモデル(11)からナゲット部(15)を中心とする直径Dの円板(31)を切り出して考える。この直径(D)の円板(31)を図3に示す。円板(31)の外周は前記八角形に外接する。この円板(31)の応力と変形を、弾性学における円板曲げ理論を用いて求める。この際、この円板(31)に作用する外力が前記ナゲット部(15)に作用する分担力である。この場合z方向のはく離荷重Fzと曲げモーメントM(MxとMyの合モーメント)である。
円板(31)外周の支持条件は前記八角形の頂点のたわみと傾斜である。ただし、これらは離散値であるので、その周方向分布をフーリエ級数で近似的に表し、これを境界条件、即ち外周支持条件として与える。
ディ.ラダイらが提案した公称構造応力式〔数1〕、〔数2〕は、図3に示す円板(31)外周支持条件が変位完全固定のもとに導かれたのものであるので、これを用いて図1の中央ナゲット部(15)の公称構造応力を求めようとするとき、〔数1〕、〔数2〕に含まれるD値をいま図1の破線で示す円の外径(D)と同じにすると、図1のスポット溶接構造における直径(D)の円周上の変位は完全固定でなく、ディ.ラダイらの円板外周支持条件(変位完全固定)と異なることから、ディ・ラダイらの式〔数1〕、〔数2〕から図1の中央ナゲット部(15)の公称構造応力σnsが精度良く得られない。これを〔数1〕、〔数2〕のD値を補正することで対処している。この補正はスポット溶接構造や荷重条件などによって異なるため、いろいろな場合についてD値をいくらにしたら良いかを別途求め、これをD値データベースD2として公称構造応力σnsを求めている。
これに対して、本発明のスポット溶接構造の疲労寿命予測方法では図3の円板外周支持条件として、スポット溶接構造における直径(D)の円の外周上の実際の変位を与えているため、より精度の高い公称構造応力σnsが得られる。また、この方法では前記D値補正が不要で、D値をいくらにしたら良いかという非常に困難な問題も解消され、したがってD値データベースD2を作成する必要がなくなる。
ナゲット分担荷重のうち,せん断成分(Fx,Fy)とねじり成分、また、6分力が同時に作用する複合荷重に対する公称構造応力は、後述する理論式より求められる。簡単なスポット溶接継手試験片を用いて複合荷重下疲労試験を行い作成した、公称構造応力σnと破断繰り返し数Nfデータベースから、前述した方法で得られた公称構造応力σnsに対応する破断寿命Nfを予測することができる。以上の流れをフローチャートにしたのが図4である。
図4に示す如く、ステップ(S11)でスポット溶接構造を形成し、ステップ(S12)で有限要素法解析用シェルモデル(FEMモデル)を形成し、ステップ(S13)で荷重データD1をFEMモデルに与えて有限要素法線形弾性解析を行うことにより、スポット溶接部(ナゲット部)の6分力と節点変位とを求める。次いで、ステップ(S14)で、算出した6分力と、節点変位に基づいて所定の演算を行って公称構造応力σnsを算出し、次いでステップ(S15)で算出した公称構造応力σnsをもとに予め実験によって作成した公称構造応力σn sに対応する破断寿命NfのデータベースD3を参照して、破断寿命Nfを算出し、これに基づいて疲労寿命を予測する。
図21と比較すると公称構造応力を算出する際にD値データベースD2が不要となっていることが分かる。
前記のように、本発明においては、D値決定問題を解消するために、曲げに対する公称構造応力を求める新しい手法を提供する。
図3のように、外径(D)、厚さ(t)の円板を考え、その中心にあるナゲット(直径d)にはく離荷重Fo,曲げモーメントMが作用するとする。曲げモーメントは図1のMx,Myの合モーメントである。また、ナゲットは剛体とする。
この新たな方法は、図3に示す如く、直径(D)、厚さ(t)の円板(31)を考える。前述したように、この円板(31)に作用する外力ははく離荷重Fz,曲げモーメントM(図1のMx,Myの合モーメント)で、これらは円形のナゲット部(15)(直径d)に作用し、円板(31)の外周の支持条件は実際のスポット溶接構造のものと合わせる。ここの解析ではナゲット部(15)を剛体として取り扱う。円板外周の境界条件は式(1),(2)のように変位完全拘束とせず、後述のようにスポット溶接構造の注目するナゲット周辺の変形状態に合わせるようにする。さらにナゲット部(15)の中心に原点を持つ極座標系(rθz)を設ける。
円板(31)は外力によって曲げられて内部にひずみと応力を生ずるので、これを求める。まず、弾性論における次の円板(31)の曲げ問題の支配方程式は次の、
前記係数を決める条件は次の2種類である。
▲1▼ナゲットが剛体である条件
図3の円形なナゲット部(15)が変形せず剛体変位をすることから、円板のたわみwと半径方向の傾斜∂w/∂rはナゲット端r=d/2で、〔数9〕となる。
▲2▼円板外周支持条件
円板外周では実際のスポット溶接構造の支持条件と合わせることが必要である。したがって、円板外周(r=D/2〉ではたわみWr=D/2と傾斜∂w/∂rは〔数10〕となる。
(i)注目するナゲットを中心に直径(D)の同心円を描き、FEMシェル解析により円周上の節点変位w,θrの離散値を求める。wは面外変位、θrは半径方向のたわみwの傾斜角
(ii)この離散値から変位w,θrの周方向分布を3次の周期スプライン補間関数で補間する。例えば補間関数を構成する節点qiの数を9とすると、外円周上の変位wのスプライン補間関数は次式となる。
変位θr(θ)r=D/2についても同様にしてスプライン補間関数で表す。
(iii)これらの補間式をフーリエ級数で表し、式(10)右辺のように円板外周の変位境界条件として与える。式(10)に含まれるフーリエ級数の係数は次式より計算される。
たわみ関数が決まると、これを用いて円板の断面力、断面モーメントは次式のように得られる。
式(15)の断面モーメントを使うと板の応力成分は次式で与えられる。
したがって、スポット溶接構造のFEMシェル解析を実施し、その結果、注目するナゲットの分担荷重とそれを中心とした直径(D)の円周上の節点変位が与えられると、式(8)の未定係数はすべて定まる。
本論文で提唱したはく離荷重と曲げモーメントに対する公称構造応力は、式(1)、(2)のように簡単な形に表わせないが、実構造に適用する際に円板外径Dをいくつにしたら良いかという問題がなくなる。さらにスポット溶接構造のナゲット周辺の変形状態を境界条件として与えることができるので、公称構造応力の精度向上を計ることができる。
本理論の解析プログラムは簡単であり、これをCAEに組み込むことによりスポット溶接構造の寿命予測が容易になる。
次に、せん断荷重とねじりモーメントによる公称構造応力の理論解について、図23に基づいてせん断力とねじりモーメントによる公称構造応力を求める方策について述べる。この問題は弾性論の平面応力問題として取り扱える。基礎方程式は、
応力関数式(20)を式(21)に代入すると応力成分は次のように表せる。
境界条件は、ナゲット端(r=d/2)で、
円板外周(r=D/2)では、
式(27)右辺は円板外周に課さられる変位であるが、これはスポット溶接構造をFEMシェル解析し、得られた節点変位から与えられる。具体的には以下のようである。
▲1▼注目するナゲットを中心に直径(D)の同心円を描き、FEMシェル解析により得られた直径(D)の円周上の節点変位u,νの離散値を求める。
▲2▼この離散値から変位u,νの周方向分布を3次の周期スプライン補間関数で補間する。例えば補間関数を構成する節点qiの数を9とすると、外円周上の変位uのスプライン補間関数は次式となる。
変位ν(θ)r=D/2についても同様にしてスプライン補間関数で表す。
(iii)これらの補間式をフーリエ級数で表し、境界条件として与える。
式(27)のuxc,uyc,θcは、ナゲットに作用する分担荷重とナゲット端周辺の応力分布との釣合式を解くことによって得ることができる。そのつり合い式は、
したがって、スポット溶接構造のFEMシェル解析を実施し、その結果、注目するナゲット部の分担荷重と、ナゲット部を中心とした直径(D)の円周上の節点変位が与えられると、式(20)の未定係数はすべて定まり、ナゲット周りの変位場と応力場が求まることになる。
次に、複合荷重下の公称構造応力について説明する。
図22のように,はく離荷重、曲げモーメント、せん断荷重、ねじりモーメントが同時に作用するとき、次式で表されるナゲット端での最大主応力を公称構造応力とし、疲労強度パラメータとして用いる。
ここで本発明の有効性を確認するために、図1のLPモデル(11)の中央ナゲット部(15)の公称構造応力を求める。解の精度を検証するためにLPモデル(11)の詳細な三次元ソリッドモデルを作成し有限要素法解析でも公称構造応力を求めた。通常自動車のボディの解析ではこのような詳細なソリッドモデルは使わず、計算負荷が軽いシェルモデルを用いる。
スポット溶接構造へ適用する前に本手法による解を検証するために、図3の円板で外周の変位を完全拘束とし、はく離荷重及び曲げモーメントが個々に作用したときのたわみ分布の厳密解との比較を行ったのが図5である。
本手法の解は、円板の外周条件(式(6)右辺)として、厳密解の円板中間円周上(はく離荷重:r=3.82(d/2)、曲げモーメント:r=3.82(d/2))のたわみと傾斜を用いたもので、厳密解とよく一致している。
図6(a),(b)は本発明によって求めたナゲット部(15)の中心を通るx軸上の応力σr分布をFEMソリッド解と比較したものである。図1でθx,θy方向に負荷した場合をLP#θx θyで表している.図6(a)はLP#90#90、図6(b)はLP#45#90であり、いずれの場合も両者は良く一致している。
図7は本発明で求めた中央ナゲットの公称構造応力を示す。いずれの荷重の場合も、本発明で求めた公称構造応力はソリッドモデルによるFEM解とよく一致している。以上から本発明によって寿命予測パラメータである公称構造応力が精度良く得られることが分かる。尚、本実施の形態では、スポット溶接構造にLPモデル11を用いたが、これに限定されるものではなく、2枚以上の板のスポット溶接構造であれば、どのような構成であっても適用することができる。
〔例2〕
上記例1は分担荷重のうちはく離荷重または曲げモーメントに対する公称構造応力を求めるためのものであるので、この例2では分担荷重としてせん断力が支配的となる場合について考える。この典型的な例として、等厚な3枚の短冊形平板を一点でスポット溶接した継手に引張せん断が作用する場合を取り上げる(図8)。
境界条件として用いる節点変位データを得るために、有限要素法シェル解析を実施する。図9は三枚重ね継手のFEMシェル解析モデルである。一般部は粗い格子状に要素分割をし、ナゲット付近はラップ長さを一辺とする正方形内を、放射方向にナゲット内を2分割、外側を4分割し、周方向に8分割した。さらにナゲット内では放射方向にシェル要素の辺に沿ってバー要素を設け、ナゲット相当の剛性を持つビーム要素を用いて上下の板をナゲット中心で結合した。FEMソルバーにCOSMOS/Mを用いて、線形弾性解析を行った。
ここで、理論を適用してナゲット周辺の応力分布を求める手順について、図9のシェル解析モデルを用いて説明する。
▲1▼図9のようなシェルモデルを作成し、FEM解析を実施する。
▲2▼▲1▼の解析から、ナゲット部の分担荷重と、直径(D)の同心円周上の8節点の半径方向および周方向変位u,νを求める。
▲3▼▲1▼で求めた節点変位については変位u,νの周方向分布を3次の周期スプライン補間関数で補間する。さらに、これをフーリエ級数で表し、これを図23の円板外周の境界条件とする。
▲4▼▲3▼の境界条件の下で、上記の理論よりナゲット周辺の応力分布を求める。
以上の手順で求めた解析結果を図10に示す。これは三枚重ね継手中央板外表面の中心軸に沿った軸応力分布で、FEMソリッド解析の結果と比較した。図24に示したように従来の理論ではナゲットから離れるにつれて応力は減衰してゼロとなり、FEMソリッド解と合わなかったが、本手法ではFEMソリッド解と良い一致を示している。
〔例3〕
図27に示した大小2枚の平板を一点スポット溶接した継手に、その両端に一様引張応力σ0が負荷した場合について負荷方向垂直応力分を求める。まず、図11に示すようにシェル要素で解析モデルを作成し、ナゲット周りの八角形の頂点に位置する節点変位を用いて上に述べた理論で応力分布を求めた結果を図12である。本理論を用いると従来法では得られなかった応力分布が求められており、FEMソリッド解とも良く一致していることが分かる。
〔例4〕
スポット溶接構造例として図13に示すLPモデルを取り上げる。平板とL形板とをあわせて三点でスポット溶接し、平板の一対の辺を全自由度拘束とし、L形板の上端に荷重Pを種々の方向に負荷する。荷重Pがθx,θy方向に作用する場合をLP_θx_θyと表すことにする。L形板のフランジ幅15mm、長さ135mm、平板の幅は45mmである。スポット溶接ピッチ45mmで、板厚は両方とも0.8mmである。
FEMシェル解析モデルとしてL形板の一部を図14に示す。一般部は粗い格子状に要素分割をし、ナゲット付近はフランジ幅w∫を一辺とする正方形内を細かく分割し、放射方向にナゲット内を2分割、外側を4分割し、周方向に4分割したShell_04と8分割したShell_08モデルである。さらにナゲット内では放射方向にシェル要素の辺に沿ってバー要素を設け、ナゲット相当の剛性を持つビーム要素を用いて上下の板をナゲット中心で結合した。
図15は本手法の解の精度を検証するために行った3次元弾性解析用FEMソリッドモデルで、L形板の一部を示す。板厚方向に4分割し、ナゲット周辺は詳細に分割した。平板も同様である。荷重条件(LP_θx_θy)と変位拘束条件は図13に示すとおりである。
上記理論を用いてナゲット周辺の応力分布を求める手順について、図14(a)のシェル解析モデルShell_04を用いて説明する。
▲1▼図14(a)のようなシェルモデルを作成し、FEM解析を実施する。
▲2▼▲1▼の解析から、ナゲット中心から放射線上(図14(a)の場合AC,BDの2本)にある節点のたわみ値を用いて、多角形と正方形とが接する節点(A,B,C,Dの4点)におけるたわみとその放射方向の傾斜を求める。
▲3▼これら4点のたわみと傾斜から周方向の分布を補間し(式(12))、これをフーリエ級数で表し、図3の円板の外周境界条件(式(10))とする。
▲4▼▲1▼のFEMシェル解析から注目するナゲットの分担力を求め、このうちはく離成分Fxと曲げ成分MxとMyを式(9)で用いる。
▲5▼▲3▼▲4▼で与えられた円板外周境界条件と荷重条件のもとで、円板外径D=w∫とし、上記理論よりはく離成分と曲げ成分に対してナゲット周辺の応力分布を求める。
▲6▼ナゲット分担荷重のせん断成分(Fx,Fy)とねじり成分Mzに対しては曲げの場合と同じように▲2▼〜▲5▼の操作を行い、ナゲット周辺の応力分布を求める。
▲7▼▲5▼▲6▼の解析結果を重ね合わせることで複合荷重下の応力が得られ、上記公称構造応力(式(32))を求める。
以上の手順で求めた解析結果を図16に示す。これは荷重LP_90_90の場合で、L形板中央ナゲットのx軸上の曲げ応力σr分布で、要素分割の影響をみるためにシェル解析モデルShell_04とShell_08の結果を比較したものである。両者の結果には大きな差異は見られず、粗い要素分割であるShell_04モデルでもFEMソリッド解と良く一致し、良好な結果が得られている。以後Shell_04モデルを用いて解析した結果を示す。図中のσ0は荷重をL形板上面の作用面積で除した平均応力である。
図17は本手法によって求めた公称構造応力をFEMソリッド解析の結果と比較したもので、いずれの荷重の場合も良好な解が得られていることが分かる。
また、同一形状の試験片にて、種々の複合荷重疲労試験が可能となるGieke,Hahnにより考案されたDC試験片(単点スポット溶接した対向カップ形試験片)を図18、代表的なスポット溶接疲労試験片である引張せん断(TS)と十字形引張試験片(CT)を図19に示す。
これらの試験片に対して本手法で公称構造応力を求め、疲労試験データを整理したのが図20で、狭いバンド幅内に整理出来ていることが分かる。
本手法を用いたスポット溶接構造の疲労寿命予測法をフローチャートで示すと図21となる。
▲1▼スポット溶接構造に対してシェルモデルを作成し、有限要素法弾性解析を実施する。
▲2▼注目するナゲットの分担荷重とナゲットを中心とした同心円周上の節点変位値を求める。
▲3▼上記で述べた理論に基づき公称構造応力σnsを算出する。
別途用意したDC試験片によるσns−N∫線図より疲労寿命を予測する。
また、本発明の請求項2に係るスポット溶接構造の疲労寿命予測方法によれば、複数枚、例えば前記2枚の板が平板とL型板とからなり、前記有限要素法解析用シェルモデルがスポット溶接構造のナゲット部のL型板のフランジ幅を一辺とする正方形内を細かく分割し、放射方向にナゲット内を2分割、外側を4分割し、周方向に8分割したので、簡易迅速にスポット溶接構造の疲労寿命を予測することができる。
更に、本発明の請求項3に係るスポット溶接構造の疲労寿命予測方法によれば、スポット溶接構造の疲労寿命の予測は、予めスポット溶接構造品について引張せん断疲労試験、はく離疲労試験及び複合荷重疲労試験を行って、公称構造応力と破断繰り返し数との関係を表すマップを形成し、該マップを公称構造応力をもとに参照して破断繰り返し数を算出することにより行うので、疲労寿命を正確に推定することができる。
【数16】
で、Eはヤング率、νはポアソン比である。
式(15)の断面モーメントを使うと板の応力成分は次式で与えられる。
【数17】
式(17)に式(15)を代入して整理すると(z=t/2)、
【数17−1】
さらに、図3の円板についてz軸方向の内外力の釣り合いと、x、y軸に関するモーメントの釣り合いを考えると、
【数18】
が得られる。境界条件によって決定した係数An〜Dnを式(8)のたわみ関数に代入した後、これを式(18)に代入すると、式(9)に含まれる定数wcとθcx、θcyはナゲットに作用するはく離荷重Fzと曲げモーメントMx、Myで表すことができる。
したがって、スポット溶接構造のFEMシェル解析を実施し、その結果、注目するナゲットの分担荷重とそれを中心とした直径(D)の円周上の節点変位が与えられると、式(8)の未定係数はすべて定まる。
Claims (3)
- 複数の板をあわせてスポット溶接構造を形成し、該スポット溶接構造に対して有限要素法解析用シェルモデルを作成し、作成した有限要素法解析用シェルモデルを用いて有限要素法線形弾性解析を行ってスポット溶接部中央のナゲット部のはく離荷重、曲げモーメント、せん断力及びねじりモーメントの分担荷重、そのナゲット部を中心に描いた直径(D)の円周上の変位を算出し、算出した該分担荷重、円周上の変位とに基づいて前記ナゲット部におけるはく離荷重、曲げモーメント、せん断力及びねじりモーメントの公称構造応力を弾性学の円板曲げ理論および2次元弾性論を用いて求め、該公称構造応力よりスポット溶接構造の疲労寿命を予測することを特徴とするスポット溶接構造の疲労寿命予測方法。
- 前記複数の板が平板とL型板とからなり、前記有限要素法解析用シェルモデルがスポット溶接構造のナゲット部のL型板のフランジ幅を一辺とする正方形内を細かく分割し、放射方向にナゲット内を2分割、外側を4分割し、周方向に8分割したことを特徴とする請求項1に記載のスポット溶接構造の疲労寿命予測方法。
- 前記スポット溶接構造の疲労寿命の予測は、予めスポット溶接構造品について引張せん断疲労試験、はく離疲労試験及び複合荷重疲労試験を行って、公称構造応力と破断繰り返し数との関係を表すマップを形成し、該マップを公称構造応力をもとに参照して破断繰り返し数を算出することにより行うことを特徴とする請求項1または請求項2に記載のスポット溶接構造の疲労寿命予測方法。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2003/005857 WO2004099761A1 (ja) | 2003-05-12 | 2003-05-12 | スポット溶接構造の疲労寿命予測方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2004099761A1 true JPWO2004099761A1 (ja) | 2006-07-27 |
JP4135946B2 JP4135946B2 (ja) | 2008-08-20 |
Family
ID=33428603
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004571574A Expired - Fee Related JP4135946B2 (ja) | 2003-05-12 | 2003-05-12 | スポット溶接構造の疲労寿命予測方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US7752917B2 (ja) |
EP (1) | EP1640702A4 (ja) |
JP (1) | JP4135946B2 (ja) |
KR (1) | KR101011844B1 (ja) |
CN (1) | CN100480674C (ja) |
AU (1) | AU2003304096A1 (ja) |
WO (1) | WO2004099761A1 (ja) |
Families Citing this family (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4165551B2 (ja) * | 2005-10-20 | 2008-10-15 | トヨタ自動車株式会社 | スポット溶接破断解析方法 |
WO2008034498A1 (en) * | 2006-09-21 | 2008-03-27 | Kompetenzzentrum - Das Virtuelle Fahrzeug Forschungsgesellschaft Mbh | Device for and method of simulating spotweld failure |
JP4743112B2 (ja) * | 2006-12-26 | 2011-08-10 | Jfeスチール株式会社 | スポット溶接構造体の疲労寿命予測方法 |
JP4599439B2 (ja) * | 2008-08-07 | 2010-12-15 | 株式会社日立製作所 | 制御装置の劣化診断システム |
US8657605B2 (en) * | 2009-07-10 | 2014-02-25 | Lincoln Global, Inc. | Virtual testing and inspection of a virtual weldment |
US7945432B2 (en) * | 2008-11-26 | 2011-05-17 | Livermore Software Technology Corporation | Spot weld failure determination method in a finite element analysis |
JP5403256B2 (ja) * | 2009-11-18 | 2014-01-29 | 株式会社Ihi | 多結晶金属材の疲労強度評価装置と方法 |
CN102192858B (zh) * | 2010-02-26 | 2013-04-03 | 株式会社神户制钢所 | T型焊接接头结构体的t接头部的疲劳特性评价方法 |
JP5573633B2 (ja) * | 2010-11-26 | 2014-08-20 | Jfeスチール株式会社 | 溶接構造体の疲労寿命予測方法 |
CN102759504B (zh) * | 2012-07-06 | 2014-04-23 | 西安交通大学 | 一种基于焊接微区性能测试的残余应力检测方法 |
CN102955881B (zh) * | 2012-10-30 | 2015-06-24 | 温州大学 | 一种计算集成电路芯片焊点热疲劳失效概率的方法 |
KR101624092B1 (ko) * | 2012-12-17 | 2016-05-24 | 신닛테츠스미킨 카부시키카이샤 | 스폿 용접부의 해석 방법, 기록 매체 및 스폿 용접부의 해석 장치 |
TWI495873B (zh) * | 2012-12-17 | 2015-08-11 | Nippon Steel & Sumitomo Metal Corp | The analytical method of the spot welding portion, the resolver of the spot welding portion, the recording apparatus, and the analyzing apparatus of the spot welding portion |
CN103076116B (zh) * | 2012-12-29 | 2015-02-25 | 上海捷锦电力新材料有限公司 | 预制式电缆附件界面压力检测方法 |
CN103995919A (zh) * | 2014-04-17 | 2014-08-20 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机结构焊点疲劳寿命计算分析方法 |
KR101683979B1 (ko) * | 2014-09-17 | 2016-12-07 | 현대자동차주식회사 | 접합부 내구성 시험 결과 데이터 분석 방법 |
US20160196378A1 (en) * | 2015-01-06 | 2016-07-07 | Livermore Software Technology Corporation | Methods and Systems For Numerically Simulating Bi-Phase Material That Changes Phase After Crossing A Directional Spatial Boundary |
CN105528480B (zh) * | 2015-11-30 | 2019-05-07 | 奇瑞汽车股份有限公司 | 焊点分析方法和装置 |
CN105445328B (zh) * | 2015-12-11 | 2019-04-05 | 工业和信息化部电子第五研究所 | 综合应力下微互连焊点的疲劳寿命评价方法、装置和系统 |
CN105447268B (zh) * | 2015-12-11 | 2018-09-21 | 奇瑞汽车股份有限公司 | 一种帽形截面的汽车前纵梁设计方法 |
US10801932B2 (en) * | 2016-02-05 | 2020-10-13 | Nippon Steel Corporation | Fracture prediction method and device, program, and recording medium |
CN105891210B (zh) * | 2016-05-04 | 2019-01-15 | 中车株洲电力机车有限公司 | 电阻焊点疲劳寿命的检测方法 |
US10746641B2 (en) * | 2017-03-24 | 2020-08-18 | Honda Motor Co., Ltd. | Peel bending moment calculation |
US11471982B2 (en) * | 2017-08-18 | 2022-10-18 | The Regents Of The University Of Michigan | Unified fatigue life evaluation method for welded structures |
CN108593462A (zh) * | 2018-06-21 | 2018-09-28 | 河海大学 | 一种结合干湿膨胀法和完全软化法测试膨胀土的强度 |
CN109175759B (zh) * | 2018-10-25 | 2024-03-19 | 程力专用汽车股份有限公司 | L型焊接接头应力值的评定方法 |
US20200134112A1 (en) * | 2018-10-26 | 2020-04-30 | Universal City Studios Llc | Systems and methods to estimate nominal stress using finite element calculated stress |
CN109766633A (zh) * | 2019-01-11 | 2019-05-17 | 江西省科学院应用物理研究所 | 有限元算法中一种用于点焊连接处的结构疲劳计算处理方法 |
RU2711082C1 (ru) * | 2019-04-05 | 2020-01-15 | Акционерное Общество "Российский Концерн По Производству Электрической И Тепловой Энергии На Атомных Станциях" (Ао "Концерн Росэнергоатом") | Способ определения остаточных напряжений в металле шва сварных соединений трубопроводов (варианты) |
JP7328516B2 (ja) * | 2019-06-28 | 2023-08-17 | 日本製鉄株式会社 | スポット溶接による溶接継手の破断予測方法、スポット溶接による溶接継手の破断予測プログラム、及び、スポット溶接による溶接継手の破断予測装置 |
CN110991111B (zh) * | 2019-11-27 | 2024-01-30 | 南京安维士传动技术股份有限公司 | 一种基于摩擦接触的风电齿轮箱行星架的疲劳计算方法 |
CN111079329B (zh) * | 2019-12-04 | 2022-10-18 | 中国直升机设计研究所 | 一种基于相似结构试验的疲劳寿命评估方法 |
CN112129660B (zh) * | 2020-08-10 | 2024-02-23 | 北京工业大学 | 一种喷丸处理后的搅拌摩擦焊接头疲劳寿命预测方法 |
CN113237746B (zh) * | 2020-12-29 | 2024-04-09 | 中国航空工业集团公司西安飞机设计研究所 | 舵面操纵试验台架及铁鸟台架舵面模拟装置强度分析方法 |
KR102653004B1 (ko) * | 2021-08-18 | 2024-04-01 | 한국과학기술원 | 비정규 상태기반 페리다이나믹스와 디지털 이미지 상관법을 이용한 불연속 영역을 포함하는 구조물의 변위장 및 응력장 측정방법 |
CN113887108A (zh) * | 2021-10-20 | 2022-01-04 | 西安热工研究院有限公司 | 一种大型风电机组轮毂疲劳寿命预测方法 |
JP7363964B1 (ja) * | 2022-05-25 | 2023-10-18 | Jfeスチール株式会社 | 車体の接合位置の最適化解析方法、装置及びプログラム |
CN117420034B (zh) * | 2023-12-18 | 2024-02-23 | 哈尔滨学院 | 一种锻炼器材焊接强度检测装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1337942B1 (en) * | 2000-11-17 | 2016-10-12 | Battelle Memorial Institute | Method and system for structural stress analysis |
JP2003149130A (ja) * | 2001-11-12 | 2003-05-21 | Univ Nihon | スポット溶接構造の疲労寿命予測方法 |
JP4165551B2 (ja) * | 2005-10-20 | 2008-10-15 | トヨタ自動車株式会社 | スポット溶接破断解析方法 |
US7640146B2 (en) * | 2006-06-27 | 2009-12-29 | Livermore Software Technology Corporation | Method and system for modeling spot welds in a finite element analysis |
FR2904577B1 (fr) * | 2006-08-03 | 2009-06-05 | Snecma Sa | Methode pour evaluer la resistance en fatigue de joints soudes |
-
2003
- 2003-05-12 US US10/556,388 patent/US7752917B2/en not_active Expired - Fee Related
- 2003-05-12 KR KR1020057021436A patent/KR101011844B1/ko not_active IP Right Cessation
- 2003-05-12 AU AU2003304096A patent/AU2003304096A1/en not_active Abandoned
- 2003-05-12 EP EP03816876A patent/EP1640702A4/en not_active Withdrawn
- 2003-05-12 JP JP2004571574A patent/JP4135946B2/ja not_active Expired - Fee Related
- 2003-05-12 WO PCT/JP2003/005857 patent/WO2004099761A1/ja active Application Filing
- 2003-05-12 CN CNB038264544A patent/CN100480674C/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN100480674C (zh) | 2009-04-22 |
AU2003304096A1 (en) | 2004-11-26 |
CN1771434A (zh) | 2006-05-10 |
US20090211366A1 (en) | 2009-08-27 |
WO2004099761A8 (ja) | 2006-05-04 |
EP1640702A4 (en) | 2009-11-11 |
JP4135946B2 (ja) | 2008-08-20 |
EP1640702A1 (en) | 2006-03-29 |
WO2004099761A1 (ja) | 2004-11-18 |
KR20060080123A (ko) | 2006-07-07 |
AU2003304096A8 (en) | 2004-11-26 |
KR101011844B1 (ko) | 2011-01-31 |
US7752917B2 (en) | 2010-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPWO2004099761A1 (ja) | スポット溶接構造の疲労寿命予測方法 | |
JP7257055B2 (ja) | 溶接構造物の統一疲労寿命評価方法 | |
Palmonella et al. | Finite element models of spot welds in structural dynamics: review and updating | |
US20030154451A1 (en) | Method, system and computer program product for multidisciplinary design analysis of structural components | |
Fermer et al. | Industrial experiences of FE‐based fatigue life predictions of welded automotive structures | |
JP6319451B2 (ja) | 構造体の部材における変形モード分析方法、変形モード分析方法を用いた構造体の補強方法、構造体の変形モード分析用コンピュータプログラム、及び記録媒体 | |
Kang et al. | Evaluation of multiaxial spot weld fatigue parameters for proportional loading | |
JP6119451B2 (ja) | 薄板レーザ溶接部の疲労寿命予測方法 | |
JP2003149130A (ja) | スポット溶接構造の疲労寿命予測方法 | |
Abambres et al. | Finite element analysis of steel structures–a review of useful guidelines | |
Rahman et al. | Fatigue life prediction of spot-welded structures: a finite element analysis approach | |
US7543732B2 (en) | Evaluation method and evaluation apparatus for spot welded portion | |
JP4862408B2 (ja) | スポット溶接構造体の疲労寿命予測方法 | |
Vivio et al. | A new analytical model for the elastic–plastic behaviour of spot welded joints subjected to orthogonal load | |
Hamidian et al. | Fatigue life enhancement of automobile wheel disc considering multi-axial stresses based on critical plane approach | |
Abdullah et al. | Strategies of Finite Element Modeling for Spot Welded Joints and its Modal Correlation with Experimental Data | |
Shahhosseini et al. | Major compliance joint modelling survey for automotive body structures | |
Osborne et al. | Finite element concept modelling methodologies for pickup truck boxes | |
US20230281353A1 (en) | User element technique for enabling coarse-mode/high-fidelity computer-aided engineering durability evaluation of spot-joined structures | |
Rhew | Strain-gage balance axial section design optimization using design of experiments | |
Wei | Development of Fatigue Life Prediction Method for Spot-Welded Joints of Mixed Materials | |
Coleman et al. | Test verification of the Cassini spacecraft dynamic model | |
Vasekar et al. | Failure analysis and optimization of universal joint yoke subjected by torsion and shear | |
Shahhosseini et al. | Beam-like major compliant joint methodology for automotive body structures | |
Yunus et al. | The improvement of bolted joints model via finite element model updating method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A529 | Written submission of copy of amendment under article 34 pct |
Free format text: JAPANESE INTERMEDIATE CODE: A5211 Effective date: 20051025 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060501 |
|
A529 | Written submission of copy of amendment under article 34 pct |
Free format text: JAPANESE INTERMEDIATE CODE: A5211 Effective date: 20051025 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060501 |
|
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: 20080508 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080602 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110613 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |