応力テンソルは、応力ベクトルの定め方の違いから、真応力テンソル・コーシー応力テンソル、公称応力テンソル・第1パイオラ・キルヒホッフ応力テンソル、第2パイオラ・キルヒホッフ応力テンソルの3種類が定義されておりいずれも(行列の形式で記述できる)2階のテンソルとなる。ただし、これらの応力テンソルに違いが生じるのは有限変形理論に基づいて物体の運動を記述した場合であり、材料力学 や応用力学 で多用されている微小変位・微小変形の仮定の下では、これらの応力テンソルはすべて真応力テンソルに一致する。
真応力テンソル(微小変形理論における応力テンソル)を σ で表すものとすると、その成分は座標軸を x , y , z と定めた3次元デカルト座標の下では、
σ
=
(
σ
x
x
σ
x
y
σ
x
z
σ
y
x
σ
y
y
σ
y
z
σ
z
x
σ
z
y
σ
z
z
)
,
or
,
σ
=
(
σ
11
σ
12
σ
13
σ
21
σ
22
σ
23
σ
31
σ
32
σ
33
)
or
,
σ
=
σ
i
j
(
e
i
⊗
e
j
)
{\displaystyle \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 軸方向の力の成分を意味する。そのため、応力テンソルの成分には、微小面の法線と力の作用方向が一致する垂直応力 (normal stress ) 成分と、一致しない(異なっている)せん断応力 (shear stress ) 成分の2種類に分類することができる。
上に示した3次元デカルト座標系における応力テンソルの成分について考えた場合、垂直応力 は
σ
x
x
,
σ
y
y
,
σ
z
z
{\displaystyle \sigma _{xx},\;\sigma _{yy},\;\sigma _{zz}}
の3成分となる。垂直応力は、力の作用面と力の作用方向とが直交し、作用面を引っ張る方向に作用した場合には引張応力 (tensile stress )、作用面を押し込む方向に作用した場合には圧縮応力 (compressive stress ) と呼ばれる。材料力学 や応用力学 、構造力学 などにおいては、引張応力が正の垂直応力となるように応力テンソルを定義するのが一般的であるが、地盤工学 (土質力学 )においては圧縮応力が正の垂直応力となるように力の正の向きを定義することもある。
一方、せん断応力 は、力の作用面の法線の向きと力の作用方向とが一致しない応力成分であり、
σ
x
y
,
σ
y
x
,
σ
y
z
,
σ
z
y
,
σ
z
x
,
σ
x
z
{\displaystyle \sigma _{xy},\;\sigma _{yx},\;\sigma _{yz},\;\sigma _{zy},\;\sigma _{zx},\;\sigma _{xz}}
の6つが該当する。なお、微小変形の力学においては、せん断応力を記号τで表すことがある。この場合の応力テンソルの表記は以下のようになる。
σ
=
(
σ
x
τ
x
y
τ
x
z
τ
y
x
σ
y
τ
y
z
τ
z
x
τ
z
y
σ
z
)
{\displaystyle \sigma ={\begin{pmatrix}\sigma _{x}&\tau _{xy}&\tau _{xz}\\\tau _{yx}&\sigma _{y}&\tau _{yz}\\\tau _{zx}&\tau _{zy}&\sigma _{z}\\\end{pmatrix}}}
応力を定義している物体内でモーメント のつりあい条件(角運動量保存則 )を満たすものと仮定する[ 注 3] と、応力テンソル(真応力テンソル)は対称テンソルとなる[ 注 4] 。すなわち、
σ
=
σ
T
{\displaystyle \sigma =\sigma ^{\mathrm {T} }}
が成り立つ。例えば、上に示した3次元デカルト座標系での成分については、
σ
x
y
=
σ
y
x
,
σ
y
z
=
σ
z
y
,
σ
z
x
=
σ
x
z
{\displaystyle \sigma _{xy}=\sigma _{yx},\;\sigma _{yz}=\sigma _{zy},\;\sigma _{zx}=\sigma _{xz}}
が成り立ち、応力テンソルσの独立な成分は6成分となることがわかる。
この性質のため、固体物性 やCAE などの分野では、独立な6成分を並べてベクトルとする表記がしばしば用いられる。これをフォークト表記 (Voigt notation )という。
σ
=
(
σ
x
x
σ
y
y
σ
z
z
σ
x
y
σ
y
z
σ
z
x
)
≡
(
σ
x
σ
y
σ
z
τ
x
y
τ
y
z
τ
z
x
)
{\displaystyle \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
′
=
A
x
{\displaystyle x'=Ax}
応力テンソルの座標系変換式は以下で表される。
σ
′
=
A
σ
A
T
{\displaystyle \sigma '=A\sigma A^{\mathrm {T} }}
ここで、 σは変換前の座標系における応力テンソル、σ' は変換後の座標系における応力テンソル、A は回転行列 、A T はA の転置行列 である。各成分で表すと以下の通りである。
(
σ
11
′
σ
12
′
σ
13
′
σ
21
′
σ
22
′
σ
23
′
σ
31
′
σ
32
′
σ
33
′
)
=
(
a
11
a
12
a
13
a
21
a
22
a
23
a
31
a
32
a
33
)
(
σ
11
σ
12
σ
13
σ
21
σ
22
σ
23
σ
31
σ
32
σ
33
)
(
a
11
a
21
a
31
a
12
a
22
a
32
a
13
a
23
a
33
)
{\displaystyle {\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}}}
ここで、aij は2つの座標間の方向余弦で、各座標軸とは下記の表のような関係となる。
x
1
{\displaystyle x_{1}}
x
2
{\displaystyle x_{2}}
x
3
{\displaystyle x_{3}}
x
1
′
{\displaystyle x'_{1}}
a
11
{\displaystyle a_{11}}
a
12
{\displaystyle a_{12}}
a
13
{\displaystyle a_{13}}
x
2
′
{\displaystyle x'_{2}}
a
21
{\displaystyle a_{21}}
a
22
{\displaystyle a_{22}}
a
23
{\displaystyle a_{23}}
x
3
′
{\displaystyle x'_{3}}
a
31
{\displaystyle a_{31}}
a
32
{\displaystyle a_{32}}
a
33
{\displaystyle a_{33}}
上式を展開すると、3次元応力状態での各応力の変換式は以下のようになる[ 7] 。
σ
11
′
=
a
11
2
σ
11
+
a
12
2
σ
22
+
a
13
2
σ
33
+
2
a
11
a
12
σ
12
+
2
a
11
a
13
σ
13
+
2
a
12
a
13
σ
23
{\displaystyle \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}}
σ
22
′
=
a
21
2
σ
11
+
a
22
2
σ
22
+
a
23
2
σ
33
+
2
a
21
a
22
σ
12
+
2
a
21
a
23
σ
13
+
2
a
22
a
23
σ
23
{\displaystyle \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}}
σ
33
′
=
a
31
2
σ
11
+
a
32
2
σ
22
+
a
33
2
σ
33
+
2
a
31
a
32
σ
12
+
2
a
31
a
33
σ
13
+
2
a
32
a
33
σ
23
{\displaystyle \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}}
σ
12
′
=
a
11
a
21
σ
11
+
a
12
a
22
σ
22
+
a
13
a
23
σ
33
+
(
a
11
a
22
+
a
12
a
21
)
σ
12
+
(
a
12
a
23
+
a
13
a
22
)
σ
23
+
(
a
11
a
23
+
a
13
a
21
)
σ
13
{\displaystyle \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}}
σ
23
′
=
a
21
a
31
σ
11
+
a
22
a
32
σ
22
+
a
23
a
33
σ
33
+
(
a
21
a
32
+
a
22
a
31
)
σ
12
+
(
a
22
a
33
+
a
23
a
32
)
σ
23
+
(
a
21
a
33
+
a
23
a
31
)
σ
13
{\displaystyle \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}}
σ
13
′
=
a
11
a
31
σ
11
+
a
12
a
32
σ
22
+
a
13
a
33
σ
33
+
(
a
11
a
32
+
a
12
a
31
)
σ
12
+
(
a
12
a
33
+
a
13
a
32
)
σ
23
+
(
a
11
a
33
+
a
13
a
31
)
σ
13
{\displaystyle \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}}
平面応力状態 での応力変換式は以下の通りである[ 8] 。
σ
11
′
=
a
11
2
σ
11
+
a
12
2
σ
22
+
2
a
11
a
12
σ
12
{\displaystyle \sigma _{11}'=a_{11}^{2}\sigma _{11}+a_{12}^{2}\sigma _{22}+2a_{11}a_{12}\sigma _{12}}
σ
22
′
=
a
21
2
σ
11
+
a
22
2
σ
22
+
2
a
21
a
22
σ
12
{\displaystyle \sigma _{22}'=a_{21}^{2}\sigma _{11}+a_{22}^{2}\sigma _{22}+2a_{21}a_{22}\sigma _{12}}
σ
12
′
=
a
11
a
21
σ
11
+
a
12
a
22
σ
22
+
(
a
11
a
22
+
a
12
a
21
)
σ
12
{\displaystyle \sigma _{12}'=a_{11}a_{21}\sigma _{11}+a_{12}a_{22}\sigma _{22}+(a_{11}a_{22}+a_{12}a_{21})\sigma _{12}}
ここで座標軸間の角度θを用いて上式を書き直した場合は以下の通りである。
σ
11
′
=
σ
11
cos
2
θ
+
σ
22
sin
2
θ
+
2
σ
12
cos
θ
sin
θ
{\displaystyle \sigma _{11}'=\sigma _{11}\cos ^{2}\theta +\sigma _{22}\sin ^{2}\theta +2\sigma _{12}\cos \theta \sin \theta }
σ
22
′
=
σ
11
sin
2
θ
+
σ
22
cos
2
θ
−
2
σ
12
cos
θ
sin
θ
{\displaystyle \sigma _{22}'=\sigma _{11}\sin ^{2}\theta +\sigma _{22}\cos ^{2}\theta -2\sigma _{12}\cos \theta \sin \theta }
σ
12
′
=
(
σ
11
−
σ
22
)
cos
θ
sin
θ
+
σ
12
(
cos
2
θ
−
sin
2
θ
)
{\displaystyle \sigma _{12}'=(\sigma _{11}-\sigma _{22})\cos \theta \sin \theta +\sigma _{12}(\cos ^{2}\theta -\sin ^{2}\theta )}
x
1
{\displaystyle x_{1}}
x
2
{\displaystyle x_{2}}
x
1
′
{\displaystyle x'_{1}}
a
11
(
=
cos
θ
)
{\displaystyle a_{11}(=\cos \theta )}
a
12
(
=
sin
θ
)
{\displaystyle a_{12}(=\sin \theta )}
x
2
′
{\displaystyle x'_{2}}
a
21
(
=
−
sin
θ
)
{\displaystyle a_{21}(=-\sin \theta )}
a
22
(
=
cos
θ
)
{\displaystyle a_{22}(=\cos \theta )}
この変換を図示する方法として、モールの応力円 が知られている。
せん断応力成分がゼロとなるように座標系を取ったときの垂直応力を主応力 (principal stress ) と呼ぶ。その座標系の基底ベクトルを応力テンソルの主軸 あるいは主応力軸 と呼ぶ。さらに主軸に垂直な面を主面 あるいは主応力面 と呼ぶ[ 9] 。各点での主軸の方向(主方向)を連ねていくと、物体の中には互いに直交する曲線群を描くことができる。これを主応力線 という[ 10] 。なお、真応力テンソル(コーシー応力テンソル)は対称テンソルであるため、ある応力状態を表す真応力テンソルに対して、せん断応力が見掛け上現れず主応力のみが垂直応力として現れる主軸が必ず一組存在する。
せん断応力がゼロとなるときの垂直応力が主応力であるが、同時に主応力はあらゆる座標系の中で垂直応力が最大、最小となる値を示している[ 11] 。3つの主応力をσ1 ≥ σ2 ≥ σ3 の関係となるようにとったとき、最大の主応力σ1 を最大主応力 、最小となる主応力σ3 を最小主応力 、これら2つに直交する主応力σ2 を中間主応力 と呼び、ある座標系での応力状態
(
σ
x
,
σ
y
,
σ
z
,
τ
x
y
,
τ
y
z
,
τ
z
x
)
{\displaystyle (\sigma _{x},\sigma _{y},\sigma _{z},\tau _{xy},\tau _{yz},\tau _{zx})}
が与えられているとき、主応力は以下の関係から求められる[ 11] 。
det
(
σ
−
λ
I
)
=
|
(
σ
x
−
λ
)
τ
x
y
τ
z
x
τ
x
y
(
σ
y
−
λ
)
τ
y
z
τ
z
x
τ
y
z
(
σ
z
−
λ
)
|
=
0
{\displaystyle \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}
上式を展開したλに関する3次方程式 の根が主応力となる。実際に上式を展開すると、
λ
3
−
J
1
λ
2
+
J
2
λ
−
J
3
=
0
,
J
1
=
σ
x
+
σ
y
+
σ
z
,
J
2
=
σ
x
σ
y
+
σ
y
σ
z
+
σ
z
σ
x
−
τ
x
y
2
−
τ
y
z
2
−
τ
z
x
2
,
J
3
=
σ
x
σ
y
σ
z
+
2
τ
x
y
τ
y
z
τ
z
x
−
σ
x
τ
y
z
2
−
σ
y
τ
z
x
2
−
σ
z
τ
x
y
2
{\displaystyle {\begin{aligned}&\lambda ^{3}-J_{1}\lambda ^{2}+J_{2}\lambda -J_{3}=0,\\&J_{1}=\sigma _{x}+\sigma _{y}+\sigma _{z},\\&J_{2}=\sigma _{x}\sigma _{y}+\sigma _{y}\sigma _{z}+\sigma _{z}\sigma _{x}-\tau _{xy}^{2}-\tau _{yz}^{2}-\tau _{zx}^{2},\\&J_{3}=\sigma _{x}\sigma _{y}\sigma _{z}+2\tau _{xy}\tau _{yz}\tau _{zx}-\sigma _{x}\tau _{yz}^{2}-\sigma _{y}\tau _{zx}^{2}-\sigma _{z}\tau _{xy}^{2}\end{aligned}}}
となる。一方、上式の根はσ1 、σ2 、σ3 となるので、上式は以下のようも書き表せる。
0
=
(
λ
−
σ
1
)
(
λ
−
σ
2
)
(
λ
−
σ
3
)
=
λ
3
−
(
σ
1
+
σ
2
+
σ
3
)
λ
2
+
(
σ
1
σ
2
+
σ
2
σ
3
+
σ
3
σ
1
)
λ
−
(
σ
1
σ
2
σ
3
)
{\displaystyle {\begin{aligned}0&=(\lambda -\sigma _{1})(\lambda -\sigma _{2})(\lambda -\sigma _{3})\\&=\lambda ^{3}-(\sigma _{1}+\sigma _{2}+\sigma _{3})\lambda ^{2}+(\sigma _{1}\sigma _{2}+\sigma _{2}\sigma _{3}+\sigma _{3}\sigma _{1})\lambda -(\sigma _{1}\sigma _{2}\sigma _{3})\end{aligned}}}
以上の2式を等値すれば、
J
1
=
σ
1
+
σ
2
+
σ
3
=
σ
x
+
σ
y
+
σ
z
=
tr
(
σ
)
,
J
2
=
σ
1
σ
2
+
σ
2
σ
3
+
σ
3
σ
1
=
σ
x
σ
y
+
σ
y
σ
z
+
σ
z
σ
x
−
τ
x
y
2
−
τ
y
z
2
−
τ
z
x
2
,
J
3
=
σ
1
σ
2
σ
3
=
σ
x
σ
y
σ
z
+
2
τ
x
y
τ
y
z
τ
z
x
−
σ
x
τ
y
z
2
−
σ
y
τ
z
x
2
−
σ
z
τ
x
y
2
=
det
(
σ
)
{\displaystyle {\begin{aligned}J_{1}&=\sigma _{1}+\sigma _{2}+\sigma _{3}=\sigma _{x}+\sigma _{y}+\sigma _{z}=\operatorname {tr} (\sigma ),\\J_{2}&=\sigma _{1}\sigma _{2}+\sigma _{2}\sigma _{3}+\sigma _{3}\sigma _{1}=\sigma _{x}\sigma _{y}+\sigma _{y}\sigma _{z}+\sigma _{z}\sigma _{x}-\tau _{xy}^{2}-\tau _{yz}^{2}-\tau _{zx}^{2},\\J_{3}&=\sigma _{1}\sigma _{2}\sigma _{3}=\sigma _{x}\sigma _{y}\sigma _{z}+2\tau _{xy}\tau _{yz}\tau _{zx}-\sigma _{x}\tau _{yz}^{2}-\sigma _{y}\tau _{zx}^{2}-\sigma _{z}\tau _{xy}^{2}=\operatorname {det} (\sigma )\end{aligned}}}
を得る。J 1 、J 2 、J 3 は、ある応力状態において座標系に関わらず常に一定値となるので応力不変量 (stress invariant )と総称される。それぞれ第一次応力不変量 、第二次応力不変量 、第三次応力不変量 と呼ぶ[ 11] 。第一次応力普遍量、第三次応力不変量は、それぞれ応力テンソルの跡 、行列式 に等しい。応力不変量は以下のように表されることもある[ 12] 。
I = J 1 , II = σ1 2 + σ2 2 + σ3 2 = tr(σ2 ), III = J 3
2次元における一般的な応力状態
2次元における主応力面
平面応力状態では σz , τyz , τzx が 0 なので、主応力は以下の関係から求められる[ 11] 。
|
(
σ
x
−
λ
)
τ
x
y
τ
x
y
(
σ
y
−
λ
)
|
=
0
{\displaystyle {\begin{vmatrix}(\sigma _{x}-\lambda )&\tau _{xy}\\\tau _{xy}&(\sigma _{y}-\lambda )\\\end{vmatrix}}=0}
上式を展開するとλに関する2次方程式 が得られ、これを解くと、平面応力状態での主応力 σ1 , σ2 は次のようになる。
σ
1
,
σ
2
=
(
σ
x
+
σ
y
)
±
(
σ
x
−
σ
y
)
2
+
4
τ
x
y
2
2
{\displaystyle \sigma _{1},\sigma _{2}={\frac {(\sigma _{x}+\sigma _{y})\pm {\sqrt {(\sigma _{x}-\sigma _{y})^{2}+4\tau _{xy}^{2}}}}{2}}}
主軸の方向は次のようになる。
θ
1
=
tan
σ
1
−
σ
x
τ
x
y
,
θ
2
=
tan
σ
2
−
σ
x
τ
x
y
{\displaystyle \theta _{1}=\tan {\frac {\sigma _{1}-\sigma _{x}}{\tau _{xy}}}\ ,\ \theta _{2}=\tan {\frac {\sigma _{2}-\sigma _{x}}{\tau _{xy}}}}
ここでθは、x 軸とσ1 、σ2 の主軸がなす角度である。
あらゆる座標系の中で最大となるせん断応力を主せん断応力 または最大せん断応力 と呼ぶ。主せん断応力が働く面は、主軸に対して45°あるいは135°傾いた面となる。主せん断応力τ1 、τ2 、τ3 は、主応力σ1 、σ2 、σ3 より次式で求まる[ 11] 。
τ
1
=
|
σ
2
−
σ
3
|
2
,
τ
2
=
|
σ
3
−
σ
1
|
2
,
τ
3
=
|
σ
1
−
σ
2
|
2
{\displaystyle \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 を受けて静的な釣り合い状態にある物体内部の任意の点では、その応力σは次の平衡方程式 あるいはつりあい方程式 を満たす。
∂
σ
x
x
∂
x
+
∂
τ
y
x
∂
y
+
∂
τ
z
x
∂
z
+
F
x
=
0
,
∂
τ
x
y
∂
x
+
∂
σ
y
y
∂
y
+
∂
τ
z
y
∂
z
+
F
y
=
0
,
∂
τ
x
z
∂
x
+
∂
τ
y
z
∂
y
+
∂
σ
z
z
∂
z
+
F
z
=
0
{\displaystyle {\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}}}
あるいは次のような書き方もされる。
σ
i
j
,
j
+
F
i
=
0
,
(
i
=
1
,
2
,
3
)
,
div
σ
+
F
=
0
{\displaystyle {\begin{aligned}&\sigma _{ij,j}+F_{i}=0,\quad (i=1,2,3),\\&\operatorname {div} \sigma +{\boldsymbol {F}}={\boldsymbol {0}}\end{aligned}}}
応力場σが平衡方程式と、表面力規定境界∂Rt における境界条件(コーシーの式)
t
=
σ
T
n
,
at
∂
R
t
{\displaystyle {\boldsymbol {t}}=\sigma ^{\mathrm {T} }{\boldsymbol {n}},\quad {\text{at}}\;\partial R_{t}}
を満たすとき、その応力場σを静的に許容な場 という。
真応力(コーシー応力)テンソルσと変形勾配 テンソルF を用いて定義される次のテンソルをパイオラ・キルヒホッフ応力テンソル (Piola-Kirchhoff stress tensor)という。
第1パイオラ・キルヒホッフ応力テンソル
Π
=
det
(
F
)
σ
(
F
−
1
)
T
{\displaystyle \Pi =\det(F)\sigma (F^{-1})^{\mathrm {T} }}
第2パイオラ・キルヒホッフ応力テンソル
K
=
F
−
1
Π
=
det
(
F
)
F
−
1
σ
(
F
−
1
)
T
{\displaystyle K=F^{-1}\Pi =\det(F)F^{-1}\sigma (F^{-1})^{\mathrm {T} }}
真応力に関するコーシーの式は上述のとおり現配置での応力ベクトルt と法線ベクトルn で表されるが、パイオラ・キルヒホッフ応力テンソルを用いても類似の関係式が成り立つ。
t
0
=
Π
n
0
,
t
^
=
K
n
0
{\displaystyle {\begin{aligned}{\boldsymbol {t}}_{0}&=\Pi {\boldsymbol {n}}_{0},\\{\hat {\boldsymbol {t}}}&=K{\boldsymbol {n}}_{0}\end{aligned}}}
ここで、
n
0
{\displaystyle {\boldsymbol {n}}_{0}}
:基準配置の微小面の法線ベクトル
t
0
{\displaystyle {\boldsymbol {t}}_{0}}
:現配置の微小面に作用している力を、基準配置の微小面の面積で割って定義される応力ベクトル
t
^
{\displaystyle {\hat {\boldsymbol {t}}}}
:現配置の微小面に作用している力を基準配置で求めなおし、それを基準配置の微小面の面積で割って定義される応力ベクトル
である。
仮想仕事の原理 を適用する際には、これらの応力テンソルと共役な関係にあるひずみテンソルは以下のようになる。
コーシー応力 - アルマンシーひずみ
第1パイオラ・キルヒホッフ応力 - 変形勾配
第2パイオラ・キルヒホッフ応力 - グリーンひずみ