Quantity in fracture mechanics; predicts stress intensity near a crack's tip
亀裂先端の極座標。 破壊力学 において 、 応力拡大係数 ( K )は、 遠隔荷重または 残留応力によって引き起こされる 亀裂 または 切欠きの 先端近傍の 応力 状態(「応力拡大」)を予測するために使用されます。 [1]これは、均質な線形 弾性 材料 に通常適用される理論的構成であり、 脆性 材料の 破壊基準 を提供するのに役立ち、 損傷許容 の分野において重要な手法です。この概念は 、亀裂先端で 小規模な 降伏を 示す材料にも適用できます。
K の大きさは 、試験片の形状、亀裂または切欠きの大きさと位置、そして材料にかかる荷重の大きさと分布に依存します。Kは次のように表されます。 [2] [3]
K = σ π a f ( a / W ) {\displaystyle K=\sigma {\sqrt {\pi a}}\,f(a/W)} ここで 、 は亀裂長さ a と試験片幅 W に依存する試験片形状の関数であり、 σ は印加応力です。 f ( a / W ) {\displaystyle f(a/W)}
線形弾性 理論によれば、き裂先端近傍の 応力分布( )は、き裂先端を原点とする 極座標 ( )において、次の形をとると予測される [4]。 σ i j {\displaystyle \sigma _{ij}} r , θ {\displaystyle r,\theta }
σ i j ( r , θ ) = K 2 π r f i j ( θ ) + h i g h e r o r d e r t e r m s {\displaystyle \sigma _{ij}(r,\theta )={\frac {K}{\sqrt {2\pi r}}}\,f_{ij}(\theta )+\,\,{\rm {higher\,order\,terms}}} ここで、 K は応力拡大係数(単位は応力×長さの 1/2 )であり、 荷重と形状によって変化する 無次元量 です。理論的には、 rが 0に近づくと応力はに近づき、応力特異点が生じ ます 。 [5] しかし実際には、この関係は先端( rが 小さい)に非常に近いところでは成り立ちません。これは、 塑性は 通常、材料の 降伏強度 を超える応力で発生し、線形弾性解はもはや適用できないためです。それでも、亀裂先端の塑性領域が亀裂長さに比べて小さい場合は、亀裂先端近くの漸近応力分布が依然として適用可能です。 f i j {\displaystyle f_{ij}} σ i j {\displaystyle \sigma _{ij}} ∞ {\displaystyle \infty }
さまざまなモードの応力拡大係数 モード I、モード II、およびモード III の亀裂荷重。 1957 年、 G. アーウィンは、亀裂周辺の応力は 応力拡大係数 と呼ばれる尺度係数で表せることを発見しました 。彼は、任意の荷重を受ける亀裂は、3 種類の線形独立亀裂モードに分解できることを発見しました。 [6] これらの荷重タイプは、図に示すようにモード I、II、または III に分類されます。モード I は、 亀裂表面が直接離れる開口 ( 引張 ) モードです。モード II は、亀裂表面が亀裂の 前縁 に垂直な方向に互いの上を滑り合う滑り (面内 せん断) モードです。モード III は、亀裂表面が互いに対して、かつ亀裂の前縁に平行に移動する引き裂き ( 面外せん断 ) モードです。モード I は、工学設計で遭遇する最も一般的な荷重タイプです。
3つの異なるモードの応力拡大係数を表すために、異なる添え字が用いられる。モードIの応力拡大係数は、 き裂開口モードに指定され適用される。モードIIの応力拡大係数は 、き裂すべりモードに適用され、モードIIIの応力拡大係数は、 き裂引張モードに適用される。これらの係数は、正式には次のように定義される。 [7] K I {\displaystyle K_{\rm {I}}} K I I {\displaystyle K_{\rm {II}}} K I I I {\displaystyle K_{\rm {III}}}
K I = lim r → 0 2 π r σ y y ( r , 0 ) K I I = lim r → 0 2 π r σ y x ( r , 0 ) K I I I = lim r → 0 2 π r σ y z ( r , 0 ) . {\displaystyle {\begin{aligned}K_{\rm {I}}&=\lim _{r\rightarrow 0}{\sqrt {2\pi r}}\,\sigma _{yy}(r,0)\\K_{\rm {II}}&=\lim _{r\rightarrow 0}{\sqrt {2\pi r}}\,\sigma _{yx}(r,0)\\K_{\rm {III}}&=\lim _{r\rightarrow 0}{\sqrt {2\pi r}}\,\sigma _{yz}(r,0)\,.\end{aligned}}}
エネルギー放出率とJ積分との関係 平面応力 条件下 では、 純粋モードIまたは純粋モードII荷重下の亀裂の ひずみエネルギー解放率 ( )は、応力集中係数と次の関係があります。 G {\displaystyle G}
G I = K I 2 ( 1 E ) {\displaystyle G_{\rm {I}}=K_{\rm {I}}^{2}\left({\frac {1}{E}}\right)} G I I = K I I 2 ( 1 E ) {\displaystyle G_{\rm {II}}=K_{\rm {II}}^{2}\left({\frac {1}{E}}\right)} ここで 、は ヤング率 、は材料の ポアソン比 である 。材料は等方性、均質性、線形弾性体であると仮定する。亀裂は初期亀裂の方向に沿って伸びると仮定する。 E {\displaystyle E} ν {\displaystyle \nu }
平面ひずみ 条件の場合 、同等の関係はもう少し複雑になります。
G I = K I 2 ( 1 − ν 2 E ) {\displaystyle G_{\rm {I}}=K_{\rm {I}}^{2}\left({\frac {1-\nu ^{2}}{E}}\right)\,} G I I = K I I 2 ( 1 − ν 2 E ) . {\displaystyle G_{\rm {II}}=K_{\rm {II}}^{2}\left({\frac {1-\nu ^{2}}{E}}\right)\,.} 純粋なモードIII負荷の場合、
G I I I = K I I I 2 ( 1 2 μ ) = K I I I 2 ( 1 + ν E ) {\displaystyle G_{\rm {III}}=K_{\rm {III}}^{2}\left({\frac {1}{2\mu }}\right)=K_{\rm {III}}^{2}\left({\frac {1+\nu }{E}}\right)} ここで、 せん断弾性係数 です 。平面ひずみにおける一般的な荷重の場合、 線形結合は次 式で表されます。 μ {\displaystyle \mu }
G = G I + G I I + G I I I . {\displaystyle G=G_{\rm {I}}+G_{\rm {II}}+G_{\rm {III}}\,.} 平面応力についても、3 つのモードの寄与を加えることで同様の関係が得られます。
上記の関係は、J積分 と応力拡大係数 を結び付けるためにも使用できる。
G = J = ∫ Γ ( W d x 2 − t ⋅ ∂ u ∂ x 1 d s ) . {\displaystyle G=J=\int _{\Gamma }\left(W~dx_{2}-\mathbf {t} \cdot {\cfrac {\partial \mathbf {u} }{\partial x_{1}}}~ds\right)\,.}
臨界応力強度係数 応力拡大係数 は 、幾何学的パラメータ(荷重タイプ)を含む適用応力の大きさを増幅するパラメータです 。どのモード状況でも、応力拡大係数は材料に加えられる荷重に正比例します。材料に非常に鋭い亀裂や V ノッチ ができると、 の最小値を 経験的に決定できます。これが亀裂を伝播させるために必要な応力拡大係数の臨界値です。 平面ひずみでのモード I 荷重 に対して決定されたこの臨界値は、材料の 臨界破壊靭性( )と呼ばれます。 の 単位は、応力に距離の平方根を乗じた値です(例:MN/m 3/2 )。 の単位は、 に達して亀裂の伝播が発生する ためには、ある臨界距離をかけて材料の破壊応力に達しなければならないことを意味します 。モード I 臨界応力拡大係数 は、 破壊力学で最も頻繁に使用される工学設計パラメータであるため、橋梁、建物、航空機、さらには鐘などに使用される破壊耐性材料を設計する場合には理解しておく必要があります。 K {\displaystyle K} Y {\displaystyle Y} K I {\displaystyle K_{\mathrm {I} }} K I c {\displaystyle K_{\mathrm {Ic} }} K I c {\displaystyle K_{\mathrm {Ic} }} K I c {\displaystyle K_{\mathrm {Ic} }} K I c {\displaystyle K_{\mathrm {Ic} }} K I c {\displaystyle K_{\mathrm {Ic} }}
研磨ではひび割れを検出できません。通常、ひび割れが見られる場合、それは 応力拡大係数 [ 要出典 ]によって予測される 臨界応力状態 に非常に近い状態です。
G基準 G 基準は 、臨界応力拡大係数(または破壊靭性)と3つのモードの応力拡大係数を関連付ける 破壊基準 である。この破壊基準は [8]のように表される。
K c 2 = K I 2 + K I I 2 + E ′ 2 μ K I I I 2 {\displaystyle K_{\rm {c}}^{2}=K_{\rm {I}}^{2}+K_{\rm {II}}^{2}+{\frac {E'}{2\mu }}\,K_{\rm {III}}^{2}} ここで、 は 平面ひずみ に対する 破壊靭性、 は 平面応力 に対する破壊靭性です。 平面応力 の臨界応力拡大係数は、 と表記されることが多いです 。 K c {\displaystyle K_{\rm {c}}} E ′ = E / ( 1 − ν 2 ) {\displaystyle E'=E/(1-\nu ^{2})} E ′ = E {\displaystyle E'=E} K c {\displaystyle K_{\rm {c}}}
例
無限領域におけるペニー型の亀裂
幅と高さ の有限板の中央に亀裂がある場合 、応力集中係数のおおよその関係は [7]である。 2 b {\displaystyle 2b} 2 h {\displaystyle 2h}
K I = σ π a [ 1 − a 2 b + 0.326 ( a b ) 2 1 − a b ] . {\displaystyle K_{\rm {I}}=\sigma {\sqrt {\pi a}}\left[{\cfrac {1-{\frac {a}{2b}}+0.326\left({\frac {a}{b}}\right)^{2}}{\sqrt {1-{\frac {a}{b}}}}}\right]\,.} 亀裂が幅の中央に位置していない場合、すなわち、位置 A における応力拡大係数は、 級数展開によって近似することができる [7] [9] d ≠ b {\displaystyle d\neq b}
K I A = σ π a [ 1 + ∑ n = 2 M C n ( a b ) n ] {\displaystyle K_{\rm {IA}}=\sigma {\sqrt {\pi a}}\left[1+\sum _{n=2}^{M}C_{n}\left({\frac {a}{b}}\right)^{n}\right]} ここで、係数は応力拡大係数曲線 [7]の フィッティングから、様々な値に対して 6で 得られる。同様の(ただし同一ではない)式が亀裂先端 B にも存在する。A と B における 応力拡大係数の代替式は [10] : 175 である。 C n {\displaystyle C_{n}} d {\displaystyle d}
K I A = σ π a Φ A , K I B = σ π a Φ B {\displaystyle K_{\rm {IA}}=\sigma {\sqrt {\pi a}}\,\Phi _{A}\,\,,K_{\rm {IB}}=\sigma {\sqrt {\pi a}}\,\Phi _{B}} どこ
Φ A := [ β + ( 1 − β 4 ) ( 1 + 1 4 sec α A ) 2 ] sec α A Φ B := 1 + [ sec α A B − 1 1 + 0.21 sin { 8 tan − 1 [ ( α A − α B α A + α B ) 0.9 ] } ] {\displaystyle {\begin{aligned}\Phi _{A}&:=\left[\beta +\left({\frac {1-\beta }{4}}\right)\left(1+{\frac {1}{4{\sqrt {\sec \alpha _{A}}}}}\right)^{2}\right]{\sqrt {\sec \alpha _{A}}}\\\Phi _{B}&:=1+\left[{\frac {{\sqrt {\sec \alpha _{AB}}}-1}{1+0.21\sin \left\{8\,\tan ^{-1}\left[\left({\frac {\alpha _{A}-\alpha _{B}}{\alpha _{A}+\alpha _{B}}}\right)^{0.9}\right]\right\}}}\right]\end{aligned}}} と
β := sin ( π α B α A + α B ) , α A := π a 2 d , α B := π a 4 b − 2 d ; α A B := 4 7 α A + 3 7 α B . {\displaystyle \beta :=\sin \left({\frac {\pi \alpha _{B}}{\alpha _{A}+\alpha _{B}}}\right)~,~~\alpha _{A}:={\frac {\pi a}{2d}}~,~~\alpha _{B}:={\frac {\pi a}{4b-2d}}~;~~\alpha _{AB}:={\frac {4}{7}}\,\alpha _{A}+{\frac {3}{7}}\,\alpha _{B}\,.} 上記の式では 、 は亀裂の中心から点 A に最も近い境界までの距離です。 上記の式は 、中心亀裂の近似式に簡略化され ないことに注意してください。 d {\displaystyle d} d = b {\displaystyle d=b}
モード I 荷重下の有限プレートの亀裂。
一軸応力を受ける板のエッジクラック
無限板:二軸応力場における斜め亀裂
平面内点力を受ける板の亀裂 長さ の亀裂を含む 寸法 の板を考えます 。 と の成分を持つ点力 が 板の 点 ( ) に作用します。 2 h × 2 b {\displaystyle 2h\times 2b} 2 a {\displaystyle 2a} F x {\displaystyle F_{x}} F y {\displaystyle F_{y}} x , y {\displaystyle x,y}
板が亀裂の大きさに比べて大きく、かつ、力の位置が亀裂に比較的近い場合、すなわち 、、、、 の場合 、 板は無限大とみなすことができる。この場合、 亀裂先端 B ( )における応力拡大係数は [11] [12] h ≫ a {\displaystyle h\gg a} b ≫ a {\displaystyle b\gg a} x ≪ b {\displaystyle x\ll b} y ≪ h {\displaystyle y\ll h} F x {\displaystyle F_{x}} x = a {\displaystyle x=a}
K I = F x 2 π a ( κ − 1 κ + 1 ) [ G 1 + 1 κ − 1 H 1 ] K I I = F x 2 π a [ G 2 + 1 κ + 1 H 2 ] {\displaystyle {\begin{aligned}K_{\rm {I}}&={\frac {F_{x}}{2{\sqrt {\pi a}}}}\left({\frac {\kappa -1}{\kappa +1}}\right)\left[G_{1}+{\frac {1}{\kappa -1}}H_{1}\right]\\K_{\rm {II}}&={\frac {F_{x}}{2{\sqrt {\pi a}}}}\left[G_{2}+{\frac {1}{\kappa +1}}H_{2}\right]\end{aligned}}} どこ
G 1 = 1 − Re [ a + z z 2 − a 2 ] , G 2 = − Im [ a + z z 2 − a 2 ] H 1 = Re [ a ( z ¯ − z ) ( z ¯ − a ) z ¯ 2 − a 2 ] , H 2 = − Im [ a ( z ¯ − z ) ( z ¯ − a ) z ¯ 2 − a 2 ] {\displaystyle {\begin{aligned}G_{1}&=1-{\text{Re}}\left[{\frac {a+z}{\sqrt {z^{2}-a^{2}}}}\right]\,,\,\,G_{2}=-{\text{Im}}\left[{\frac {a+z}{\sqrt {z^{2}-a^{2}}}}\right]\\H_{1}&={\text{Re}}\left[{\frac {a({\bar {z}}-z)}{({\bar {z}}-a){\sqrt {{\bar {z}}^{2}-a^{2}}}}}\right]\,,\,\,H_{2}=-{\text{Im}}\left[{\frac {a({\bar {z}}-z)}{({\bar {z}}-a){\sqrt {{\bar {z}}^{2}-a^{2}}}}}\right]\end{aligned}}} ここで 、 は 平面ひずみ 、 は 平面 応力 、 は ポアソン比 である。 先端 B における応力拡大係数 は、 z = x + i y {\displaystyle z=x+iy} z ¯ = x − i y {\displaystyle {\bar {z}}=x-iy} κ = 3 − 4 ν {\displaystyle \kappa =3-4\nu } κ = ( 3 − ν ) / ( 1 + ν ) {\displaystyle \kappa =(3-\nu )/(1+\nu )} ν {\displaystyle \nu } F y {\displaystyle F_{y}}
K I = F y 2 π a [ G 2 − 1 κ + 1 H 2 ] K I I = − F y 2 π a ( κ − 1 κ + 1 ) [ G 1 − 1 κ − 1 H 1 ] . {\displaystyle {\begin{aligned}K_{\rm {I}}&={\frac {F_{y}}{2{\sqrt {\pi a}}}}\left[G_{2}-{\frac {1}{\kappa +1}}H_{2}\right]\\K_{\rm {II}}&=-{\frac {F_{y}}{2{\sqrt {\pi a}}}}\left({\frac {\kappa -1}{\kappa +1}}\right)\left[G_{1}-{\frac {1}{\kappa -1}}H_{1}\right]\,.\end{aligned}}} 先端 A における応力拡大係数( )は上記の関係から決定できる。 位置 における 荷重については 、 x = − a {\displaystyle x=-a} F x {\displaystyle F_{x}} ( x , y ) {\displaystyle (x,y)}
K I ( − a ; x , y ) = − K I ( a ; − x , y ) , K I I ( − a ; x , y ) = K I I ( a ; − x , y ) . {\displaystyle K_{\rm {I}}(-a;x,y)=-K_{\rm {I}}(a;-x,y)\,,\,\,K_{\rm {II}}(-a;x,y)=K_{\rm {II}}(a;-x,y)\,.} 同様に負荷についても 、 F y {\displaystyle F_{y}}
K I ( − a ; x , y ) = K I ( a ; − x , y ) , K I I ( − a ; x , y ) = − K I I ( a ; − x , y ) . {\displaystyle K_{\rm {I}}(-a;x,y)=K_{\rm {I}}(a;-x,y)\,,\,\,K_{\rm {II}}(-a;x,y)=-K_{\rm {II}}(a;-x,y)\,.} およびの 成分を持つ局所的な力の作用下にあるプレートの亀裂 。 F x {\displaystyle F_{x}} F y {\displaystyle F_{y}}
プレートに負荷のかかった亀裂 亀裂が点力によって荷重を受けている場合 、 点 B における応力集中係数 は [7] F y {\displaystyle F_{y}} y = 0 {\displaystyle y=0} − a < x < a {\displaystyle -a<x<a}
K I = F y 2 π a a + x a − x , K I I = − F x 2 π a ( κ − 1 κ + 1 ) . {\displaystyle K_{\rm {I}}={\frac {F_{y}}{2{\sqrt {\pi a}}}}{\sqrt {\frac {a+x}{a-x}}}\,,\,\,K_{\rm {II}}=-{\frac {F_{x}}{2{\sqrt {\pi a}}}}\left({\frac {\kappa -1}{\kappa +1}}\right)\,.} 力が の間に均一に分布している場合、先端 B における応力拡大係数 は − a < x < a {\displaystyle -a<x<a}
K I = 1 2 π a ∫ − a a F y ( x ) a + x a − x d x , K I I = − 1 2 π a ( κ − 1 κ + 1 ) ∫ − a a F y ( x ) d x , . {\displaystyle K_{\rm {I}}={\frac {1}{2{\sqrt {\pi a}}}}\int _{-a}^{a}F_{y}(x)\,{\sqrt {\frac {a+x}{a-x}}}\,{\rm {d}}x\,,\,\,K_{\rm {II}}=-{\frac {1}{2{\sqrt {\pi a}}}}\left({\frac {\kappa -1}{\kappa +1}}\right)\int _{-a}^{a}F_{y}(x)\,{\rm {d}}x,\,.} プレートに生じたひび割れ。
無限板における平行亀裂の積み重ね 出典: [13]
亀裂間隔が亀裂長さよりもはるかに大きい場合(h >> a)、隣接する亀裂間の相互作用効果は無視でき、応力集中係数は長さ 2a の単一亀裂の応力集中係数と等しくなります。
すると、亀裂先端の応力拡大係数は
K I = σ π a {\displaystyle {\begin{aligned}K_{\rm {I}}&=\sigma {\sqrt {\pi a}}\end{aligned}}}
亀裂の長さが間隔よりもはるかに大きい場合 (a >> h )、亀裂は半無限亀裂のスタックと見なすことができます。
すると、亀裂先端の応力拡大係数は
K I = σ h {\displaystyle {\begin{aligned}K_{\rm {I}}&=\sigma {\sqrt {h}}\end{aligned}}}
コンパクト張力試験片
片刃ノッチ曲げ試験片
参照
参考文献 ^ ab Anderson, TL (2005). 破壊力学:基礎と応用 . CRC Press. ^ Soboyejo, WO (2003). 「11.6.2 亀裂の駆動力と相似の概念」. 工学材料の機械的性質 . マルセル・デッカー. ISBN 0-8247-8900-8 . OCLC 300921090。 ^ Janssen, M. (Michael) (2004). 破壊力学 . Zuidema, J. (Jan), Wanhill, RJH (第2版). ロンドン: Spon Press. p. 41. ISBN 0-203-59686-2 . OCLC 57491375。 ^ 多田 宏 、 PC Paris 、 George R. Irwin (2000年2月)。 『ひび割れの応力解析ハンドブック (第3版)』アメリカ機械学会。 ^ abc Liu, M.; et al. (2015). 「丸い先端のノッチにおける応力に対する改良された半解析解」 (PDF) . 工学破壊力学 . 149 : 134– 143. doi :10.1016/j.engfracmech.2015.10.004. S2CID 51902898. ^ ab Suresh, S. (2004). 材料の疲労 . ケンブリッジ大学出版局. ISBN 978-0-521-57046-6 。 ^ abcdefg Rooke, DP; Cartwright, DJ (1976). 応力集中係数の概要 . HMSO 国防省. 調達担当部署. ^ Sih, GC; Macdonald, B. (1974)、「工学問題への破壊力学の応用 - ひずみエネルギー密度破壊基準」、 工学破壊力学 、 6 (2): 361– 386、 doi :10.1016/0013-7944(74)90033-2 ^ Isida, M.、1966、 「偏心亀裂ストリップの張力に対する応力集中係数」 、ASME応用力学部門の取引、v. 88、p.94。 ^ Kathiresan, K.; Brussat, TR; Hsu, TM (1984). 高度な寿命解析手法. アタッチメントラグのき裂進展解析手法 . 飛行力学研究所, 空軍ライト航空研究所, AFSC WP空軍基地, オハイオ州. ^ ab Sih, GC; Paris, PC & Erdogan, F. (1962), "平面伸長と板曲げ問題におけるき裂先端応力集中係数", Journal of Applied Mechanics , 29 (2): 306– 312, Bibcode :1962JAM....29..306S, doi :10.1115/1.3640546 ^ エルドガン、F. ( 1962)「任意荷重下における共線切断板の応力分布について」、 第4回米国応用力学会議論文集 、 1 : 547-574 ^ クンドゥ、トリビクラム (2008年1月30日). 破壊力学の基礎. doi :10.1201/9781439878057. ISBN 9780429110054 。 ^ ab Bower, AF (2009). 固体の応用力学 . CRC Press.
外部リンク Kathiresan, K.; Hsu, TM; Brussat, TR, 1984, 「Advanced Life Analysis Methods. Volume 2. Attachment Lugのき裂進展解析法」 www.fracturemechanics.org の応力拡大係数(Bob McGinty 著)