応力

応力（おうりょく、ストレス、英: stress）とは、物体の内部に生じる力の大きさや作用方向を表現するために用いられる物理量である。物体の変形や破壊などに対する負担の大きさを検討するのに用いられる。

連続体力学
法則 質量保存の法則 運動量保存の法則 エネルギー保存の法則 クラウジウス–デュエムの不等式 固体力学 固体 · 変形 · 弾性 · 弾性波 · 弾塑性 · 塑性 · フックの法則 · 応力 · ひずみ · 有限変形理論 · レオロジー · 粘弾性 · 超弾性 流体力学 流体 · 流体静力学 流体動力学 · 粘度 · ニュートン流体 非ニュートン流体 表面張力 科学者 ニュートン · ストークス · ナビエ · コーシー · フック · ベルヌーイ

応力
量記号	σ
次元	T−2 L−1 M
種類	2階テンソル
SI単位	パスカル (Pa)
FPS重力単位	重量ポンド毎平方インチ (psi)

この物理量には応力ベクトル (stress vector) と応力テンソル (stress tensor) の2つがあり、単に「応力」といえば応力テンソルのことを指すことが多い。応力テンソルは座標系などを特別に断らない限り、主に2階の混合テンソルおよび混合ベクトルとして扱われる（混合テンソルについてはテンソル積#テンソル空間を参照）。応力ベクトルと応力テンソルは、ともに連続体内部に定義した微小面積に作用する単位面積あたりの力として定義される。そのため、それらの単位は、SIではPa (N/m²)、重力単位系ではkgf/mm²で、圧力と同じである。

異なる定義

応力という物理量は、分野によって全く異なる使われ方がなされている。即ち、土木・建築分野においては連続体内部の面にかかる力（単位：ニュートン（N））のことを応力と呼び、その単位断面積当たりの力を「応力度(stress intensity)（単位：N/m² = Pa）」と呼んでいる。

応力の定義の違い

物理量	計量法、物理学、材料工学、機械工学など	土木・建築分野
力（単位：N　）	力	応力
単位断面積当たりの力（単位：N/m2 = Pa）	応力	応力度

以下では、計量法体系の定義にあるとおり、応力を「単位断面積当たりの力」の意味で用いる。

応力ベクトル

応力ベクトルとは、物体表面あるいは物体内に仮想的な微小面を考えたとき、その微小面に作用する単位面積あたりの力であり、ベクトル（1階のテンソル）で表される。後述する応力テンソルの説明にあるように、応力テンソルσの各成分の第1の下添字は「応力成分を考えている微小面の法線の向き」を、第2の下添字は「考えている微小面に作用する力の向き」をそれぞれ表している。このことから明らかなように、微小面の単位法線ベクトルを n とすると、その微小面での応力ベクトル t は次のように与えられる。

${\boldsymbol {t}}=\sigma ^{\mathrm {T} }{\boldsymbol {n}}$

この式はコーシーの式と呼ばれる。例えば、3次元デカルト座標系 (x , y , z ) において、単位法線ベクトルを ${\boldsymbol {n}}=(n_{x},n_{y},n_{z})=(\cos \alpha ,\cos \beta ,\cos \gamma )$ と表すと、応力ベクトルの成分 $t_{x},\;t_{y},\;t_{z}$ は次のようになる。

${\begin{pmatrix}t_{x}\\t_{y}\\t_{z}\end{pmatrix}}={\begin{pmatrix}\sigma _{xx}n_{x}+\sigma _{yx}n_{y}+\sigma _{zx}n_{z}\\\sigma _{xy}n_{x}+\sigma _{yy}n_{y}+\sigma _{zy}n_{z}\\\sigma _{xz}n_{x}+\sigma _{yz}n_{y}+\sigma _{zz}n_{z}\end{pmatrix}}$

応力テンソル

応力テンソルは、応力ベクトルの定め方の違いから、真応力テンソル、公称応力テンソル、第2パイオラ・キルヒホッフ応力テンソルの3種類が定義されておりいずれも2階のテンソルとなる。ただし、これらの応力テンソルに違いが生じるのは有限変形理論に基づいて物体の運動を記述した場合であり、材料力学や応用力学で多用されている微小変位・微小変形の仮定の下では、これらの応力テンソルはすべて真応力テンソルに一致する。

真応力テンソルを σ で表すものとすると、その成分は座標軸を x , y , z と定めた3次元デカルト座標の下では、

$\sigma ={\begin{pmatrix}\sigma _{xx}&\sigma _{xy}&\sigma _{xz}\\\sigma _{yx}&\sigma _{yy}&\sigma _{yz}\\\sigma _{zx}&\sigma _{zy}&\sigma _{zz}\end{pmatrix}},\ {\mbox{or}},\ \sigma ={\begin{pmatrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\end{pmatrix}}\ {\mbox{or}},\ \sigma =\sigma _{ij}({\boldsymbol {e}}_{i}\otimes {\boldsymbol {e}}_{j})$

のように表される。e_i 等は座標軸 x , y , z 方向の基底ベクトルである。このとき、各成分の第1の下添字は「応力成分を考えている微小面の法線の向き」を、第2の下添字は「考えている微小面に作用する力の向き」をそれぞれ表している。例えば、σ_xy とは、法線の方向がx 軸の向きに一致する微小面において考えている、y 軸方向の力の成分を意味する。そのため、応力テンソルの成分には、微小面の法線と力の作用方向が一致する垂直応力成分と、一致しないせん断応力成分の2種類がある。応力テンソルは、三次元空間の各面に作用する力の成分をまとめた2階のテンソルであり、一般的には3×3の行列として表される。

垂直応力とせん断応力

垂直応力は $\sigma _{xx},\;\sigma _{yy},\;\sigma _{zz}$ の3成分となる。垂直応力は、力の作用面と力の作用方向とが直交し、作用面を引っ張る方向に作用した場合には引張応力、作用面を押し込む方向に作用した場合には圧縮応力と呼ばれる。材料力学や応用力学、構造力学などにおいては、引張応力が正の垂直応力となるように応力テンソルを定義するのが一般的であるが、地盤工学においては圧縮応力が正の垂直応力となるように力の正の向きを定義することもある。

一方、せん断応力は、力の作用面の法線の向きと力の作用方向とが一致しない応力成分であり、$\sigma _{xy},\;\sigma _{yx},\;\sigma _{yz},\;\sigma _{zy},\;\sigma _{zx},\;\sigma _{xz}$の6つが該当する。微小変形の力学においては、せん断応力を記号τで表すことがある。この場合の応力テンソルの表記は以下のようになる。

$\sigma ={\begin{pmatrix}\sigma _{x}&\tau _{xy}&\tau _{xz}\\\tau _{yx}&\sigma _{y}&\tau _{yz}\\\tau _{zx}&\tau _{zy}&\sigma _{z}\\\end{pmatrix}}$

応力テンソルの対称性

応力を定義している物体内でモーメントのつりあい条件を満たすものと仮定すると、応力テンソルは対称テンソルとなる。すなわち、

$\sigma =\sigma ^{\mathrm {T} }$

が成り立つ。例えば、

$\sigma _{xy}=\sigma _{yx},\;\sigma _{yz}=\sigma _{zy},\;\sigma _{zx}=\sigma _{xz}$

が成り立ち、応力テンソルσの独立な成分は6成分となることがわかる。

この性質のため、固体物性やCAEなどの分野では、独立な6成分を並べてベクトルとする表記がしばしば用いられる。これをフォークト表記という。

$\sigma ={\begin{pmatrix}\sigma _{xx}\\\sigma _{yy}\\\sigma _{zz}\\\sigma _{xy}\\\sigma _{yz}\\\sigma _{zx}\\\end{pmatrix}}\equiv {\begin{pmatrix}\sigma _{x}\\\sigma _{y}\\\sigma _{z}\\\tau _{xy}\\\tau _{yz}\\\tau _{zx}\\\end{pmatrix}}$

任意座標系への応力の変換

応力テンソルの座標変換

真応力はテンソル量であり、座標系によってその成分は変化することとなる。以下のように座標系を変換する。$$x'=Ax$$応力テンソルの座標系変換式は以下で表される。

$\sigma '=A\sigma A^{\mathrm {T} }$

ここで、 σは変換前の座標系における応力テンソル、σ' は変換後の座標系における応力テンソル、A は回転行列 、A^T はA の転置行列である。各成分で表すと以下の通りである。

${\begin{pmatrix}\sigma '_{11}&\sigma '_{12}&\sigma '_{13}\\\sigma '_{21}&\sigma '_{22}&\sigma '_{23}\\\sigma '_{31}&\sigma '_{32}&\sigma '_{33}\\\end{pmatrix}}={\begin{pmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\\\end{pmatrix}}{\begin{pmatrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\\\end{pmatrix}}{\begin{pmatrix}a_{11}&a_{21}&a_{31}\\a_{12}&a_{22}&a_{32}\\a_{13}&a_{23}&a_{33}\\\end{pmatrix}}$

ここで、a_ij は2つの座標間の方向余弦で、各座標軸とは下記の表のような関係となる。

	$x_{1}$	$x_{2}$	$x_{3}$
$x'_{1}$	$a_{11}$	$a_{12}$	$a_{13}$
$x'_{2}$	$a_{21}$	$a_{22}$	$a_{23}$
$x'_{3}$	$a_{31}$	$a_{32}$	$a_{33}$

3次元応力状態での各応力の変換式は以下のようになる。

$\sigma _{11}'=a_{11}^{2}\sigma _{11}+a_{12}^{2}\sigma _{22}+a_{13}^{2}\sigma _{33}+2a_{11}a_{12}\sigma _{12}+2a_{11}a_{13}\sigma _{13}+2a_{12}a_{13}\sigma _{23}$

$\sigma _{22}'=a_{21}^{2}\sigma _{11}+a_{22}^{2}\sigma _{22}+a_{23}^{2}\sigma _{33}+2a_{21}a_{22}\sigma _{12}+2a_{21}a_{23}\sigma _{13}+2a_{22}a_{23}\sigma _{23}$

$\sigma _{33}'=a_{31}^{2}\sigma _{11}+a_{32}^{2}\sigma _{22}+a_{33}^{2}\sigma _{33}+2a_{31}a_{32}\sigma _{12}+2a_{31}a_{33}\sigma _{13}+2a_{32}a_{33}\sigma _{23}$

$\sigma _{12}'=a_{11}a_{21}\sigma _{11}+a_{12}a_{22}\sigma _{22}+a_{13}a_{23}\sigma _{33}+(a_{11}a_{22}+a_{12}a_{21})\sigma _{12}+(a_{12}a_{23}+a_{13}a_{22})\sigma _{23}+(a_{11}a_{23}+a_{13}a_{21})\sigma _{13}$

$\sigma _{23}'=a_{21}a_{31}\sigma _{11}+a_{22}a_{32}\sigma _{22}+a_{23}a_{33}\sigma _{33}+(a_{21}a_{32}+a_{22}a_{31})\sigma _{12}+(a_{22}a_{33}+a_{23}a_{32})\sigma _{23}+(a_{21}a_{33}+a_{23}a_{31})\sigma _{13}$

$\sigma _{13}'=a_{11}a_{31}\sigma _{11}+a_{12}a_{32}\sigma _{22}+a_{13}a_{33}\sigma _{33}+(a_{11}a_{32}+a_{12}a_{31})\sigma _{12}+(a_{12}a_{33}+a_{13}a_{32})\sigma _{23}+(a_{11}a_{33}+a_{13}a_{31})\sigma _{13}$

平面応力状態での応力変換式は以下の通りである。

$\sigma _{11}'=a_{11}^{2}\sigma _{11}+a_{12}^{2}\sigma _{22}+2a_{11}a_{12}\sigma _{12}$

$\sigma _{22}'=a_{21}^{2}\sigma _{11}+a_{22}^{2}\sigma _{22}+2a_{21}a_{22}\sigma _{12}$

$\sigma _{12}'=a_{11}a_{21}\sigma _{11}+a_{12}a_{22}\sigma _{22}+(a_{11}a_{22}+a_{12}a_{21})\sigma _{12}$

この変換を図示する方法として、モールの応力円が知られている。

主応力

せん断応力成分がゼロとなるように座標系を取ったときの垂直応力を主応力と呼ぶ。その座標系の基底ベクトルを応力テンソルの主軸あるいは主応力軸と呼ぶ。さらに主軸に垂直な面を主面あるいは主応力面と呼ぶ。各点での主軸の方向を連ねていくと、物体の中には互いに直交する曲線群を描くことができる。これを主応力線という。真応力テンソルは対称テンソルであるため、ある応力状態を表す真応力テンソルに対して、せん断応力が見掛け上現れず主応力のみが垂直応力として現れる主軸が必ず一組存在する。

せん断応力がゼロとなるときの垂直応力が主応力であるが、同時に主応力はあらゆる座標系の中で垂直応力が最大、最小となる値を示している。3つの主応力をσ₁ ≥ σ₂ ≥ σ₃ の関係となるようにとったとき、最大の主応力σ₁ を最大主応力 、最小となる主応力σ₃ を最小主応力 、これら2つに直交する主応力σ₂ を中間主応力 と呼び、ある座標系での応力状態$(\sigma _{x},\sigma _{y},\sigma _{z},\tau _{xy},\tau _{yz},\tau _{zx})$が与えられているとき、主応力は以下の関係から求められる。

$\operatorname {det} (\sigma -\lambda I)={\begin{vmatrix}(\sigma _{x}-\lambda )&\tau _{xy}&\tau _{zx}\\\tau _{xy}&(\sigma _{y}-\lambda )&\tau _{yz}\\\tau _{zx}&\tau _{yz}&(\sigma _{z}-\lambda )\end{vmatrix}}=0$

平面応力状態における主応力

2次元における一般的な応力状態

2次元における主応力面

平面応力状態では σ_z, τ_yz, τ_zx が 0 なので、主応力は以下の関係から求められる。

${\begin{vmatrix}(\sigma _{x}-\lambda )&\tau _{xy}\\\tau _{xy}&(\sigma _{y}-\lambda )\\\end{vmatrix}}=0$

これを解くと、平面応力状態での主応力 σ₁, σ₂ は次のようになる。

$\sigma _{1},\sigma _{2}={\frac {(\sigma _{x}+\sigma _{y})\pm {\sqrt {(\sigma _{x}-\sigma _{y})^{2}+4\tau _{xy}^{2}}}}{2}}$

主軸の方向は次のようになる。

$\theta _{1}=\tan {\frac {\sigma _{1}-\sigma _{x}}{\tau _{xy}}}\ ,\ \theta _{2}=\tan {\frac {\sigma _{2}-\sigma _{x}}{\tau _{xy}}}$

ここでθは、x 軸とσ₁、σ₂の主軸がなす角度である。

主せん断応力

あらゆる座標系の中で最大となるせん断応力を主せん断応力または最大せん断応力と呼ぶ。主せん断応力が働く面は、主軸に対して45°あるいは135°傾いた面となる。主せん断応力τ₁、τ₂、τ₃ は、主応力σ₁、σ₂、σ₃ より次式で求まる。

$\tau _{1}={\frac {|\sigma _{2}-\sigma _{3}|}{2}}\quad ,\quad \tau _{2}={\frac {|\sigma _{3}-\sigma _{1}|}{2}}\quad ,\quad \tau _{3}={\frac {|\sigma _{1}-\sigma _{2}|}{2}}$

一般的に、主応力とは異なり、主せん断応力が働く面にはせん断応力だけでなく垂直応力も働く。

平衡方程式

外力F を受けて静的な釣り合い状態にある物体内部の任意の点では、その応力σは次の平衡方程式あるいはつりあい方程式を満たす。

${\begin{aligned}{\frac {\partial \sigma _{xx}}{\partial x}}+{\frac {\partial \tau _{yx}}{\partial y}}+{\frac {\partial \tau _{zx}}{\partial z}}+F_{x}&=0,\\{\frac {\partial \tau _{xy}}{\partial x}}+{\frac {\partial \sigma _{yy}}{\partial y}}+{\frac {\partial \tau _{zy}}{\partial z}}+F_{y}&=0,\\{\frac {\partial \tau _{xz}}{\partial x}}+{\frac {\partial \tau _{yz}}{\partial y}}+{\frac {\partial \sigma _{zz}}{\partial z}}+F_{z}&=0\end{aligned}}$

あるいは次のような書き方もされる。

${\begin{aligned}&\sigma _{ij,j}+F_{i}=0,\quad (i=1,2,3),\\&\operatorname {div} \sigma +{\boldsymbol {F}}={\boldsymbol {0}}\end{aligned}}$

応力場σが平衡方程式と、表面力規定境界∂R_t における境界条件

${\boldsymbol {t}}=\sigma ^{\mathrm {T} }{\boldsymbol {n}},\quad {\text{at}}\;\partial R_{t}$

を満たすとき、その応力場σを静的に許容な場という。

パイオラ･キルヒホッフ応力テンソル

真応力テンソルσと変形勾配テンソルF を用いて定義される次のテンソルをパイオラ・キルヒホッフ応力テンソルという。

第1パイオラ・キルヒホッフ応力テンソル

$\Pi =\det(F)\sigma (F^{-1})^{\mathrm {T} }$

第2パイオラ・キルヒホッフ応力テンソル

$K=F^{-1}\Pi =\det(F)F^{-1}\sigma (F^{-1})^{\mathrm {T} }$

仮想仕事の原理を適用する際には、これらの応力テンソルと共役な関係にあるひずみテンソルは以下のようになる。

• 真応力 - アルマンシーひずみ
• 第1パイオラ･キルヒホッフ応力 - 変形勾配
• 第2パイオラ･キルヒホッフ応力 - グリーンひずみ

偏差応力

偏差応力(deviatoric stress) は、応力テンソルからその等方成分を差し引いたものとして定義される。物体に等方的な圧縮・引張り以外のせん断変形が生じた場合に、偏差応力が発生する。偏差応力 dev[σ] は次のように定義される。

$\operatorname {dev} [\sigma ]:=\sigma +pI$

ここでI は2階の単位テンソル、

$p:=-{\frac {1}{3}}\operatorname {tr} [\sigma ]=-{\frac {\sigma _{kk}}{3}}$

は非決定（静水圧）応力であり、平均応力（3つの主応力の平均値）のマイナスに等しい。pI は平均応力テンソルと呼ばれる。

偏差応力の固有値 s₁, s₂, s₃ は、元の応力テンソルの固有値（主応力）と次の関係がある。

$s_{i}=\sigma _{i}+p,\quad i=1,2,3$

偏差応力の主軸は元の応力テンソルの主軸と一致する。

材料の降伏と等価応力

上記にあるとおり、応力は3次元的なテンソルである。一般の応力について材料の特性値を調べるのは困難であるため、降伏に対して等価とみなせる1軸応力に対応するスカラー量である等価応力に換算すると便利である。等価応力は材料の降伏する条件に応じて以下のようなものがある。

最大主応力説

ある点で最大主応力 σ₁ が材料の降伏を決定するというのが最大主応力説である。すなわち、

$\sigma _{1}\geq \sigma _{\mathrm {Y} }$

が降伏の条件である。ここで σ_Y は材料の降伏応力である。最大主応力説はガラスなどの脆性材料で良く当てはまる。

せん断ひずみエネルギー説

単位体積あたりのせん断ひずみエネルギーが限界を越えると、材料が破壊されるという説である。ともいう。全ひずみエネルギーから静ひずみエネルギーを差し引いたせん断ひずみエネルギー U を評価基準とする。

$U={\frac {1+\nu }{6E}}((\sigma _{1}-\sigma _{2})^{2}+(\sigma _{2}-\sigma _{3})^{2}+(\sigma _{3}-\sigma _{1})^{2})$

ここで、νはポアソン比、E はヤング率である。

せん断ひずみエネルギーに比例する相当応力をMisesの相当応力σ_Misesとよび、主応力を用いて以下の式で表される。

$\sigma _{\mathrm {Mises} }^{2}={\frac {(\sigma _{1}-\sigma _{2})^{2}+(\sigma _{2}-\sigma _{3})^{2}+(\sigma _{3}-\sigma _{1})^{2}}{2}}$

降伏条件は以下の通り。

$\sigma _{\mathrm {Mises} }\geq \sigma _{\mathrm {Y} }$

せん断ひずみエネルギー説は鋼材などの延性材料に比較的良く当てはまる。

最大せん断応力説

延性材料が降伏するときすべりが観察されることに着目し、最大せん断応力が降伏を決定するという説を最大せん断応力説、またはトレスカの応力説と呼ぶ。このときに用いられる相当応力をトレスカ応力とよび、最大せん断応力を記号 τ_max 、トレスカ応力をσ_Trescaで表すと、主応力とは次式に示す関係がある。

${\begin{aligned}\tau _{\mathrm {max} }&={\frac {1}{2}}(\mathrm {max} |\sigma _{i}-\sigma _{j}|),\\\sigma _{\mathrm {Tresca} }&=2\tau _{\mathrm {max} }\end{aligned}}$

降伏条件は以下の通り。

$\sigma _{\mathrm {Tresca} }\geq \sigma _{\mathrm {Y} }$

最大せん断応力説も延性材料に当てはまることが多い。また、σ_Tresca ≥ σ₁, σ_Tresca ≥ σ_Mises であり、上記2説に対して安全側であることから評価基準として利用されることがある。

残留応力

残留応力とは、外力が作用していない物体の内部に生じている応力である。残留応力は機械的または熱的な原因で物体に不均一に弾塑性変形が生じることにより発生する。

→詳細は「残留応力」を参照

脚注

参考文献

• 日本機械学会 編『機械工学辞典』（第2版）丸善、2007年1月20日。ISBN 978-4-88898-083-8。 
• 大路清嗣、中井善一『材料強度』コロナ社、2010年10月20日。ISBN 978-4-339-04039-5。 
• 村上敬宜『弾性力学』（第14版）養賢堂、2004年3月30日。ISBN 978-4842501215。 
• 小林英男; 轟章『固体の弾塑性力学』数理工学社、2007年。ISBN 978-4-901683-51-7。 
• 渋谷陽二『塑性の物理』森北出版、2011年。ISBN 978-4-627-66761-7。 
• 野田直剛; 谷川義信; 須見尚文; 辻知章『基礎弾性力学』（8版）日新出版、1999年。ISBN 4-8173-0146-5。 

関連項目

• 力、力学、応用力学、材料力学
• 変形、ひずみ、ひずみゲージ
• 弾性係数、降伏、耐力
• 地震、アスペリティ
• 熱応力
• 応力集中、応力拡大係数
• マクスウェルの応力テンソル
• 力の流れ
• ビリアル応力

外部リンク

• 『応力』 - コトバンク