応力集中のソースを表示

[[File:FEM example.svg|thumb|400px|[[FEM]]による応力分布の計算例<br>赤い部分が応力集中部（高効力部）]]
'''応力集中'''(おうりょくしゅうちゅう、{{Lang-en-short|stress concentration}})とは、物体の形状変化部で局所的に[[応力]]が増大する現象である<ref name = "機械工学辞典_151"/>。[[機械]]・[[構造物]]の[[疲労 (材料)|疲労]]破壊や[[脆性]]破壊では、この応力集中を起こす部分が[[破壊]]の起点となることが多い。

== 概要 ==
[[File:HoleForceLines.svg|thumb|150px|孔が存在する場合の応力線イメージ図]]
物体に力が負荷されると、物体内部に[[応力]]が発生する。一般に、内部の応力の分布は一様ではなく、力の[[負荷]]の仕方や物体の形状によって、応力は場所ごとに変化する。特に、[[孔]]や溝、段といった一様な形状が変化する部分では応力分布が乱れ、形状変化部の前後に比べて局所的に応力が増大する。このような現象を'''応力集中'''と呼び、応力集中を起こす箇所を'''応力集中部'''あるいは'''切欠き'''と呼ぶ。

以下に代表的な応力集中が問題になる事例を示す。
* 物体の外形が変化する場合（例：[[段]]や[[溝]]のような局所的なものから、外板の形状変化部のような骨組みの変化まで）
* 物体に空洞が存在する場合（例：貫通穴や材料中の空洞[[欠陥]]）
* 集中[[荷重]]を受ける場合（例：荷重を受ける範囲が十分に小さいと見なせる場合）
* 別の物体の接触（例：[[ヘルツの接触応力]]）
* 材料の[[弾性率]]が異なる物質が介在する場合（例：[[金属]]材料に含まれる非金属介在物）

応力集中がどの程度起こるかは、[[弾性]]力学、[[塑性]]力学といった個体[[力学]]理論による応力分布の[[解析]]により解明される。しかし、応力分布の厳密解が判明している問題は限られており、特に3次元問題の解析は2次元問題よりも非常に難しく、厳密解が得られる問題は非常に限られている<ref name = "応力集中の考え方_45"/>。そのため、実際の複雑な形状の応力分布を計算する方法としては、[[有限要素法]](FEM)による数値解析が行われている。また、実物で応力分布を計測する方法としては、[[光弾性]]応力測定、熱弾性応力測定、ひずみゲージによる応力測定がある。

応力集中部あるいは切欠きは応力が高まることから破壊の起点となり易い。疲労破壊では切欠きから発生した[[き裂]]が進展して破壊に至ることが多い。切欠きが存在する場合は、存在しない場合よりも疲労強度が低くなり、このような効果を'''切欠き効果'''と呼ぶ。脆性破壊においても、切欠きの存在により脆性破壊が起き易くなる。[[鋼|鉄鋼]]のような[[展延性|延性]]材料でも切欠きの存在により脆性的な破壊を起こすことがあり、このような現象を'''切欠脆性'''と呼ぶ<ref name = "機械工学辞典_315-316"/>。

== 応力集中係数 ==
[[File:Kt arbre epaule traction.svg|thumb|200px|引張を受けるR付段を持つ丸棒の応力集中係数のグラフ<br>公称応力は直径が小さい部分での断面積で定義されている。]]

応力集中の度合いを表すために、応力集中による最大応力を基準となる応力で除した'''応力集中係数'''(stress concentration factor)を用いる。
:<math> K_t = \frac{\sigma_{max}}{\sigma_n}</math>
:ここで
:''K''<sub>t</sub>：応力集中係数
:''σ''<sub>max</sub>：応力集中部の最大応力
:''σ''<sub>n</sub>：公称応力
応力集中係数の他に形状係数(shape factor)とも呼ぶ<ref name = "応力集中_15"/>。記号としては<math>K_t</math>や<math>\alpha</math>が用いられる。

公称応力は応力集中係数を定義するための基準の応力で任意に定義されるものである。公称応力の取り方としては大きく3つの取り方がある<ref name = "応力集中_15"/>。
* 穴などの応力集中要素がある場合、これらの要素により母体の断面そのものが減少し、応力分布の乱れによる応力集中とは別に正味断面積の平均応力が高まるが、この平均応力で公称応力を定義する場合。
* 応力集中要素による減少断面積を使わずに定義する公称応力。応力集中部手前の一様形状における遠方応力を使用する場合。
* 応力集中要素による最大応力を含む断面で定義するが、断面積の計算する際には応力集中要素は存在しない(切欠きが埋まっている)場合の断面積を使用する場合。
ハンドブックや教科書などに種々の場合の応力集中係数がまとめられているが、公称応力の取り方に注意して利用する必要がある。

== ひずみ集中 ==

[[ひずみ]]に関しても同様の係数、'''ひずみ集中係数'''(factor of strain concentration)が定義される<ref name = "機械工学辞典_1085"/>。
: <math> K_\epsilon = \frac{\epsilon_{max}}{\epsilon_n}</math>
: ここで
: ''K''<sub>ε</sub>：ひずみ集中係数
: ''ε''<sub>max</sub>：最大ひずみ
: ''ε''<sub>n</sub>：公称ひずみ

弾性範囲内ではK<sub>t</sub>=K<sub>ε</sub>だが、応力が[[降伏 (物理)|降伏条件]]を満たして弾性範囲を脱すると、ひずみ集中係数は応力集中係数と異なってくる<ref name = "疲労設計便覧_138-139"/>。弾性範囲を超えると、塑性により応力集中は緩和されるが、ひずみ集中は緩和されない<ref name = "機械工学辞典_1085"/>。部材が静的負荷を受ける場合、切欠き底で弾性範囲を超えて塑性ひずみが生じるようになると、塑性領域での応力集中係数は弾性領域での応力集中係数よりも減少し、ひずみ集中係数は弾性領域よりも増大する<ref name = "疲労設計便覧_138-139"/>。

塑性範囲でのひずみ集中係数と応力集中係数の関係の推定式としては、ノイバー則(Neuber's rule)がよく用いられる<ref name = "疲労設計便覧_140-142"/>。ノイバー(Neuber)は、深い切欠きを有する部材が面外せん断を受ける場合の計算に基づいて次式の関係を導いた<ref name = "疲労設計便覧_140-142"/>。
: <math> \sqrt{K_\sigma K_\epsilon} = K_t </math>
: ここで
: ''K''<sub>t</sub>：弾性応力集中係数
: ''K''<sub>σ</sub>：塑性応力集中係数
: ''K''<sub>ε</sub>：塑性ひずみ集中係数
ただし、有限要素法による検証によると、公称応力が材料の降伏応力を超えて切欠き底が全面降伏するような条件では、ノイバー則はひずみ集中係数にやや過大な値を与える傾向がある<ref name = "疲労設計便覧_138-139"/><ref name = "切欠きの非弾性応力・ひずみ集中の簡易推定法_32"/>。

== 2次元問題 ==
実際の物体は3次元であるが、3次元物体の応力分布を求めるのは容易ではないので、厚みあるいは高さを0とした2次元の形状の応力解析が行われてきた。3次元物体が[[平面応力]]あるいは[[平面ひずみ状態]]にあるものを2次元問題として扱える。実際の物体では完全な平面応力あるいは平面ひずみにあるものは無いが、例えば薄板や高剛性材料に挟まれた物体などを、近似的に平面応力あるいは平面ひずみ状態と見なして2次元問題から得られた解や知見を当てはめることができる。弾性率が異なる別の物体が介在する場合を除き、2次元問題の応力分布は問題の物体の弾性率によらずに形状と境界条件(荷重条件、拘束条件)のみに依存する<ref name = "応力集中_18"/>。

以下に代表的な2次元形状の切欠きの応力集中の弾性解析解を示す。

=== 円孔の応力集中 ===
[[File:Stress concentration by a hole.png|thumb|無限遠から一様引張応力を受ける無限板中の円孔による応力集中]]

遠方から一様な引張応力を受ける無限板に存在する[[円 (数学)|円孔]]について、最大応力を含む線上での垂直応力分布は次式で与えられる<ref name = "応力集中の考え方_1"/>。
:<math> \sigma_y = \sigma_0 \left( 1+\frac{a^2}{2x^2}+\frac{3a^4}{2x^4}\right) </math>
:ここで
:''σ''<sub>y</sub>：円孔中心を通り遠方応力に平行な線(x軸)上の垂直応力
:''σ''<sub>0</sub>：遠方引張応力
:''a''：円孔半径
:''x''：円孔中心を通り遠方応力に平行な線(x軸)上の円孔中心からの距離
最大応力は上式で''x'' = ''a''(円孔縁)の位置で発生し、この点で円孔の応力集中係数は円孔半径の値によらず次のようになる。
:<math> K_t = 3</math>

=== 楕円孔の応力集中 ===
[[File:Stress concentration by an ellipse (2).png|thumb|無限遠から一様引張応力を受ける無限板中の楕円孔による応力集中]]
遠方から長軸に垂直な一様引張応力を受ける無限板に存在する[[楕円|楕円孔]]について、最大応力を含む線上での応力分布は次式で与えられる<ref name = "弾性力学_56"/>。
:<math>
\begin{align}
 \sigma_y = & \sigma_0 \bigg[ \frac{1}{\xi^2-1}\left(\xi^2+\frac{a}{a-b}\right) \\
            & -\frac{1}{(\xi^2-1)^2} \left\{ \frac{\xi^2}{2} \left( \frac{a-b}{a+b}-\frac{a+3b}{a-b}\right)-\frac{b(a+b)}{(a-b)^2} \right\} \\
            & -\frac{4\xi^2}{(\xi^2-1)^3} \left( \frac{\xi^2 b}{a+b}-\frac{b}{a-b} \right)\frac{a}{a-b} \bigg] \\
\end{align}
</math>

:<math> \xi=\frac{x+\sqrt{x^2-c^2}}{c} \quad,\quad c=\sqrt{a^2-b^2} </math>
:ここで
:''σ''<sub>y</sub>：楕円孔長軸(x軸)上の垂直応力
:''σ''<sub>0</sub>：遠方引張応力
:''a''：楕円孔長辺
:''b''：楕円孔短辺
:''x''：楕円孔長軸(x軸)上の楕円孔中心からの距離
最大応力は上式で''x'' = ''a''(楕円孔長辺縁)の位置で発生し、この点で応力集中係数は次のようになる。
:<math> K_t =1+\frac{2a}{b}</math>
あるいは''x'' = ''a''の点における曲率半径を用いて次のようにも表される。
:<math> K_t =1+2\sqrt{\frac{a}{\rho}} \quad,\quad \rho=b^2/a </math>

楕円孔は、''b'' → 0とすればき裂（グリフィスき裂）の問題となり、また、等価楕円の概念を利用して任意形状の切欠きの応力集中系数を近似できる場合があるなど、他の問題への応用の広がりが大きい<ref name = "応力集中の考え方_5"/>。

=== 等価楕円の概念 ===
[[File:Equivalent ellipse.png|thumb|応力集中係数を概算するための等価楕円の概念、孔と切欠きの場合]]
理想的な円孔や楕円孔と異なる複雑な形状の孔や切欠きの応力集中係数を簡易に近似計算するために、'''等価楕円'''(equivalent ellipse)の考え方が平野により考案された<ref name = "二次元彈性体の形状係数の研究"/>。等価楕円の考え方では、板中の孔に対しては孔の曲率半径''ρ''と孔の全長2''a''と等しい楕円を、板縁切欠きに対しては切欠き底半径''ρ''と切欠き深さ''a''と等しい楕円を当てはめて応力集中係数を推定する。すなわち、これら2つのパラメータが応力集中に対しては影響が大きく、他の形状要素(例えば切欠きの開き角など)の影響は相対的に小さいと考える方法である<ref name = "応力集中_44"/>。

等価楕円による推定は万能ではなく、例えば、大きな孔縁にある非常に小さな切欠きの応力集中では、等価楕円による推定値は正確な応力集中係数値から大きく外れる<ref name = "応力集中_46"/>。ただし、上手く使用すれば実用上十分な近似値を推定できる。例えば遠方で引張を受ける無限板縁のV形切欠きの場合では、開き角が''θ'' = 90°、切欠き深さと切欠き底半径の比が''a''/''ρ'' = 4のときで、正確な数値計算結果では''K''<sub>t</sub> = 5.274、等価楕円による計算では''K''<sub>t</sub> = 5である<ref name = "応力集中の考え方_17"/>。

=== き裂の応力集中 ===
[[File:Stress concentration by a crack.png|thumb|無限遠から一様引張応力を受ける無限板中のき裂による応力集中]]
[[き裂]]の応力集中は楕円孔の短長を''b'' → 0とした極限として考えることができる。遠方からき裂に垂直な一様引張応力を受ける無限板に存在する貫通直線き裂について、き裂延長線上での応力分布は次式で与えられる<ref name = "弾性力学_59"/>。
:<math> \sigma_y = \frac{\sigma_0 x}{\sqrt{x^2-a^2}} </math>
:ここで
:''σ''<sub>y</sub>：き裂延長線(x軸)上の垂直応力
:''σ''<sub>0</sub>：遠方引張応力
:''a''：き裂半長
:''x''：き裂延長線(x軸)上のき裂中心からの距離
き裂先端の応力に注目すると、''x'' → ''a''では''σ''<sub>y</sub>→∞となり、応力集中係数も∞となる。よってき裂の問題では、材料中の最大応力のみで材料強度を論じることができない。

さらに、上式をき裂先端を原点に''x''座標を取り直し、き裂先端近傍(''x''/''a'' ≪ 1)に絞って考えると応力分布は次式となる。
:<math> \sigma_y \cong \frac{\sigma_0 \sqrt{a}}{\sqrt{2x}} </math>
:ここで
:''x''：き裂延長線(x軸)上のき裂先端からの距離
上式から、き裂先端近傍部分の応力は<math>\sqrt{x}</math>に反比例すること、応力分布はパラメータ<math>\sigma_0 \sqrt{a}</math>により決定されることが分かる。<math>\sigma_0 \sqrt{a}</math>に<math>\sqrt{\pi}</math>に乗じれば、[[破壊力学]]でき裂先端の応力状態を表すパラメータである[[応力拡大係数]]となる。

{{Main|破壊力学|応力拡大係数|グリフィス理論}}

== 3次元問題 ==
3次元問題の解析は2次元問題よりも非常に難しく、厳密解が得られる問題は非常に限られている。代わりにFEMによる数値解析が行われる。2次元問題の応力分布は、弾性係数が異なる別の物体が介在する場合を除き、問題の物体の弾性係数によらずに形状と境界条件のみにより決定されるが、3次元問題の場合は応力分布が物体の弾性係数(ポアソン比)にも依存する<ref name = "応力集中_19"/>。以下に3次元問題の代表的な切欠きの応力集中の弾性解析解を示す。

=== 球状空洞による応力集中 ===
[[File:Stress concentration by a spherical cavity (2).png|thumb|無限遠から一様引張応力を受ける無限物体中の球状空洞による応力集中]]

遠方から一様引張応力を受ける無限物体に存在する球状空洞について、最大応力を含む面上での応力分布は次式で与えられる<ref name = "応力集中の考え方_45"/>。
:<math> \sigma_z = \sigma_0 \left[ 1+\frac{4-5\nu}{2(7-5\nu)}\frac{a^3}{r^3}+\frac{9}{2(7-5\nu)}\frac{a^5}{r^5}\right] </math>
:ここで
:''σ''<sub>z</sub>：球中心を通り遠方応力に平行な面(x-y面)上の垂直応力
:''σ''<sub>0</sub>：遠方引張応力
:''a''：球半径
:''r''：球中心を通り遠方応力に平行な面(x-y面)上の球中心からの距離
:''ν''：ポアソン比
最大応力は上式でr=a(球縁)の位置で発生し、この点で応力集中係数は次のようになる。
:<math> K_t =\frac{27-15\nu}{2(7-5\nu)}</math>

== 脚注 ==
{{脚注ヘルプ}}
{{Reflist|2|refs=
<ref name = "機械工学辞典_151">[[#機械工学辞典|「機械工学辞典」p.151]]</ref>
<ref name = "機械工学辞典_315-316">[[#機械工学辞典|「機械工学辞典」pp.315-316]]</ref>
<ref name = "機械工学辞典_1085">[[#機械工学辞典|「機械工学辞典」p.1085]]</ref>
<ref name = "応力集中_18">[[#応力集中|「応力集中」p.18]]</ref>
<ref name = "応力集中_19">[[#応力集中|「応力集中」p.19]]</ref>
<ref name = "応力集中_15">[[#応力集中|「応力集中」p.15]]</ref>
<ref name = "応力集中_44">[[#応力集中|「応力集中」p.44]]</ref>
<ref name = "応力集中_46">[[#応力集中|「応力集中」p.46]]</ref>
<ref name = "応力集中の考え方_1">[[#応力集中の考え方|「応力集中の考え方」p.1]]</ref>
<ref name = "応力集中の考え方_5">[[#応力集中の考え方|「応力集中の考え方」p.5]]</ref>
<ref name = "応力集中の考え方_17">[[#応力集中の考え方|「応力集中の考え方」p.17]]</ref>
<ref name = "応力集中の考え方_45">[[#応力集中の考え方|「応力集中の考え方」p.45]]</ref>
<ref name = "弾性力学_56">[[#弾性力学|「弾性力学」p.56]]</ref>
<ref name = "弾性力学_59">[[#弾性力学|「弾性力学」p.59]]</ref>
<ref name = "疲労設計便覧_140-142">[[#疲労設計便覧|「疲労設計便覧」pp.140-142]]</ref>
<ref name = "疲労設計便覧_138-139">[[#疲労設計便覧|「疲労設計便覧」pp.138-139]]</ref>
<ref name = "切欠きの非弾性応力・ひずみ集中の簡易推定法_32">[[#切欠きの非弾性応力・ひずみ集中の簡易推定法|「切欠きの非弾性応力・ひずみ集中の簡易推定法」p.32]]</ref>
<ref name = "二次元彈性体の形状係数の研究">[[#二次元彈性体の形状係数の研究|「二次元彈性体の形状係数の研究」]]</ref>

}}

== 参考文献 ==
* {{cite book|和書|editor=日本機械学会 |title=機械工学辞典 |publisher=丸善 |date=2007-01-20 |edition=第2版 |ISBN=978-4-88898-083-8 |ref=機械工学辞典}}
* {{cite book|和書|editor=日本材料学会 |title=疲労設計便覧 |publisher=養賢堂 |date=2008-10-01 |edition=第3版 |ISBN=978-4-8425-9501-6 |ref=疲労設計便覧}}
* {{cite book|和書|author=西田正孝 |title=応力集中 |publisher=森北出版 |date=1993-12-25 |edition=増補版 |ISBN=978-4-627-94029-1 |ref=応力集中}}
* {{cite book|和書|author=村上敬宜 |title=弾性力学 |publisher=養賢堂 |date=2004-03-30 |edition=第14版 |ISBN=978-4842501215 |ref=弾性力学}}
* {{cite book|和書|author=村上敬宜 |title=応力集中の考え方 |publisher=養賢堂 |date=2005-07-01 |edition=第1版 |ISBN=978-4842503745 |ref=応力集中の考え方}}
* {{Cite journal |和書|author=左近淑郎, 伊達新吾 |title =切欠きの非弾性応力・ひずみ集中の簡易推定法 |journal=材料 |volume=48 |pages=30-37 |year=1999 |publisher=日本材料学会 |url=https://doi.org/10.2472/jsms.48.30 |issn=05145163 |naid=110002293614 |doi=10.2472/jsms.48.30 |ref=切欠きの非弾性応力・ひずみ集中の簡易推定法}}
* {{Cite journal |和書|author=平野冨士夫 |title=二次元彈性体の形状係数の研究(第2報) |journal=日本機械学會論文集 |issue=16 |volume=55 |pages=52-58 |date=1950-11-05 |publisher=日本機械学 |url=https://doi.org/10.1299/kikai1938.16.55_52 |issn=00290270 |naid=110002349142 |doi=10.1299/kikai1938.16.55_52 |ref=二次元彈性体の形状係数の研究}}

== 外部リンク==
{{Commonscat}}
*{{機械工学事典|id=07:1001454}}

{{DEFAULTSORT:おうりよくしゆうちゆう}}
[[Category:材料工学]]
[[Category:固体力学]]