変形された3Dオブジェクト内のあらゆる点における機械的応力の表現
コーシー応力テンソル 3次元における応力の構成要素
一般的な記号
σ SI単位 パスカル (Pa) その他のユニット
平方インチ当たりのポンド数 (psi)、 バール SI基本単位 では Pa = kg ⋅ m −1 ⋅ s −2 テンソル 寸法
ら
−
1
ま
T
−
2
{\displaystyle {\mathsf {L}}^{-1}{\mathsf {M}}{\mathsf {T}}^{-2}}
連続体力学 において 、 コーシー応力テンソル (記号 、 オーギュスタン=ルイ・コーシー にちなんで名付けられた)は、 真応力テンソル [1] または単に 応力テンソルとも呼ばれ、 変形 状態、配置、または構成にある材料内部の点における 応力 の状態を完全に定義します 。2 次テンソル は 9 つの要素で構成され 、単位長さ 方向ベクトル e と、 e に垂直な仮想面を横切る 牽引ベクトル T ( e ) を関連付けます。
σ
{\displaystyle {\boldsymbol {\sigma }}}
σ
私
じ
{\displaystyle \sigma_{ij}}
T
(
e
)
=
e
⋅
σ
または
T
じ
(
e
)
=
∑
私
σ
私
じ
e
私
。
{\displaystyle \mathbf {T} ^{(\mathbf {e} )}=\mathbf {e} \cdot {\boldsymbol {\sigma }}\quad {\text{または}}\quad T_{j}^{(e)}=\sum _{i}\sigma _{ij}e_{i}.}
[あ]
応力テンソルと牽引力ベクトルの SI 基本単位は、応力スカラーに対応する ニュートン毎平方メートル (N/m 2 ) または パスカル (Pa) です。単位ベクトルは 無次元 です。
コーシー応力テンソルは、 座標系の変化に応じて テンソル変換法則に従います。この変換法則をグラフィカルに表現したものが、応力の モール円 です。
コーシー応力テンソルは、 小さな変形を 受ける物体の 応力解析に使用され、 線形弾性理論 の中心的な概念です 。 有限変形とも呼ばれる大きな変形の場合、 ピオラ・キルヒホッフ応力テンソル 、 ビオ応力テンソル 、 キルヒホッフ応力テンソル など、他の応力の尺度が必要になります 。
線型運動量保存 の原理によれば 、連続体が静的平衡状態にある場合、その物体のすべての質点におけるコーシー応力テンソルの成分が平衡方程式(ゼロ加速度に対する コーシーの運動方程式 )を満たすことが証明できます。同時に、 角運動量保存 の原理によれば、平衡状態には 任意の点に対する モーメントの合計がゼロであることが必要であり、これにより 応力テンソルは対称 であり、したがって元の 9 つの独立した応力成分ではなく 6 つの独立した応力成分のみを持つという結論が導かれます。ただし、偶応力、つまり単位体積あたりのモーメントが存在する場合、応力テンソルは非対称です。これは、 クヌーセン数が 1 に近い場合、または連続体が非ニュートン流体である場合にも当てはまり、これは ポリマー などの回転不変流体につながる可能性があります 。
け
ん
→
1
{\displaystyle K_{n}\rightarrow 1}
応力テンソルには、選択した座標系や応力テンソルが作用する領域要素に依存しない不変量がいくつか存在します。これらは応力テンソルの 3 つの 固有値 であり、主応力と呼ばれます。
オイラー・コーシー応力原理 – 応力ベクトル
図2.1a 連続体の内部表面の 差にかかる接触力と偶応力の内部分布 。表面によって分離された連続体の2つの部分間の相互作用の結果として生じる。
d
ス
{\displaystyle dS}
ス
{\displaystyle S}
図2.1b 連続体の内部表面の 差にかかる接触力と偶応力の内部分布 。表面によって分離された連続体の2つの部分間の相互作用の結果として生じる。
d
ス
{\displaystyle dS}
ス
{\displaystyle S}
図 2.1c 法線ベクトル n を持つ内部表面 S 上の応力ベクトル。検討中の平面の方向によっては、応力ベクトルは必ずしもその平面に垂直、 つまり に平行であるとは限らず 、2 つの成分に分解できます。1 つは平面に垂直な成分で、 法線応力 と呼ばれ、もう 1 つはこの平面に平行な成分で、 せん断応力 と呼ばれます。
ん
{\displaystyle \mathbf {n} }
σ
ん
{\displaystyle \sigma _{\mathrm {n}}}
τ
{\displaystyle \tau}
オイラー ・コーシーの応力原理は 、 物体を分割する任意の表面(実数または虚数)上で、物体の一部が他の部分に及ぼす作用は、物体を分割する表面上の分布力と結合のシステムと等価(等分)であることを 述べています。 [2] そして、それは表面上で定義され 、表面の単位ベクトルに連続的に依存すると仮定される、 牽引ベクトル と呼ばれる フィールドによって表されます 。 [3] [4] : p.66–96
T
(
ん
)
{\displaystyle \mathbf {T} ^{(\mathbf {n} )}}
ス
{\displaystyle S}
ん
{\displaystyle \mathbf {n} }
オイラー・コーシーの応力原理を定式化するには、図 2.1a または 2.1b に示すように、連続体を 2 つのセグメントに分割する 内部の物質点を通る仮想表面を考えます (切断面図、または表面で囲まれた連続体内の任意の体積を持つ図のいずれかを使用できます )。
ス
{\displaystyle S}
ポ
{\displaystyle P}
ス
{\displaystyle S}
ニュートン と オイラー の古典力学によれば 、物体の運動は外部から加えられた 力の作用によって生み出され、その力は 表面力 と 体積力 の2種類あると想定される 。 [5] したがって、物体または物体の一部に加えられる
力の合計は次のように表される。
ふ
{\displaystyle \mathbf {F} }
b
{\displaystyle \mathbf {b} }
ふ
{\displaystyle {\mathcal {F}}}
ふ
=
b
+
ふ
{\displaystyle {\mathcal {F}}=\mathbf {b} +\mathbf {F} }
この記事では、コーシー応力テンソルに関連する表面力についてのみ説明します。
物体が外部表面力または 接触力 を受ける場合、 オイラーの運動方程式 に従って、内部接触力および接触モーメントは、物体内の点から点へ、また、連続体の一部が他の部分に機械的に接触することにより、分割面を介して 1 つのセグメントから他のセグメントへ伝達されます(図 2.1a および 2.1b)。 法線 ベクトル を持つ を含む 領域要素上では 、力の分布は、点 P で及ぼされる接触力 および表面モーメントに等分されます 。特に、接触力は次のように表されます。
ふ
{\displaystyle \mathbf {F} }
ス
{\displaystyle S}
Δ
ス
{\displaystyle \デルタS}
ポ
{\displaystyle P}
ん
{\displaystyle \mathbf {n} }
Δ
ふ
{\displaystyle \Delta \mathbf {F} }
Δ
ま
{\displaystyle \Delta \mathbf {M} }
Δ
ふ
=
T
(
ん
)
Δ
ス
{\displaystyle \Delta \mathbf {F} =\mathbf {T} ^{(\mathbf {n} )}\,\Delta S}
ここで 平均表面牽引力 は です 。
T
(
ん
)
{\displaystyle \mathbf {T} ^{(\mathbf {n} )}}
コーシーの応力原理は、 [6] : p.47–102 が 非常に小さくなり、ゼロに近づくにつれて、比率 は になり 、偶応力ベクトルは 消滅すると主張しています。連続体力学の特定の分野では、偶応力は消滅しないと想定されています。ただし、連続体力学の古典的な分野では、 偶応力やボディモーメントを考慮しない
非 極性材料を扱っています。
Δ
ス
{\displaystyle \デルタS}
Δ
ふ
/
Δ
ス
{\displaystyle \Delta \mathbf {F} /\Delta S}
d
ふ
/
d
ス
{\displaystyle d\mathbf {F} /dS}
Δ
ま
{\displaystyle \Delta \mathbf {M} }
合成ベクトル は 表面牽引力 [7] 、 応力ベクトル [ 8] 、 牽引力 [4] 、 牽引力ベクトル [6] とも呼ばれ 、
法線ベクトルを持つ平面に関連付けられた 点で 次のように定義されます 。
d
ふ
/
d
ス
{\displaystyle d\mathbf {F} /dS}
T
(
ん
)
=
T
私
(
ん
)
e
私
{\displaystyle \mathbf {T} ^{(\mathbf {n} )}=T_{i}^{(\mathbf {n} )}\mathbf {e} _{i}}
P
{\displaystyle P}
n
{\displaystyle \mathbf {n} }
T
i
(
n
)
=
lim
Δ
S
→
0
Δ
F
i
Δ
S
=
d
F
i
d
S
.
{\displaystyle T_{i}^{(\mathbf {n} )}=\lim _{\Delta S\to 0}{\frac {\Delta F_{i}}{\Delta S}}={dF_{i} \over dS}.}
この式は、応力ベクトルが物体内の位置とそれが作用する平面の方向によって決まることを意味します。
これは、内部接触力のバランスをとる作用によって 接触力密度 または コーシー牽引力場 [5] が生成されることを意味します。これは、特定の時間における 物体の 特定の構成における物体の体積全体にわたる内部接触力の分布を表します 。これはベクトル場ではありません。特定の物質点の位置だけでなく 、法線ベクトルによって定義される表面要素の局所的な方向にも依存するためです 。 [9]
T
(
n
,
x
,
t
)
{\displaystyle \mathbf {T} (\mathbf {n} ,\mathbf {x} ,t)}
t
{\displaystyle t}
x
{\displaystyle \mathbf {x} }
n
{\displaystyle \mathbf {n} }
検討中の平面の方向によっては、応力ベクトルは必ずしもその平面に垂直、 つまり 平行であるとは限らず 、2つの成分に分解できます(図2.1c)。
n
{\displaystyle \mathbf {n} }
σ
n
=
lim
Δ
S
→
0
Δ
F
n
Δ
S
=
d
F
n
d
S
,
{\displaystyle \mathbf {\sigma _{\mathrm {n} }} =\lim _{\Delta S\to 0}{\frac {\Delta F_{\mathrm {n} }}{\Delta S}}={\frac {dF_{\mathrm {n} }}{dS}},}
ここで、 差分領域に対する 力の法線成分は
d
F
n
{\displaystyle dF_{\mathrm {n} }}
d
F
{\displaystyle d\mathbf {F} }
d
S
{\displaystyle dS}
τ
=
lim
Δ
S
→
0
Δ
F
s
Δ
S
=
d
F
s
d
S
,
{\displaystyle \mathbf {\tau } =\lim _{\Delta S\to 0}{\frac {\Delta F_{\mathrm {s} }}{\Delta S}}={\frac {dF_{\mathrm {s} }}{dS}},}
ここで、 は 微分表面積に対する 力の接線成分です 。せん断応力はさらに 2 つの互いに直交するベクトルに分解できます。
d
F
s
{\displaystyle dF_{\mathrm {s} }}
d
F
{\displaystyle d\mathbf {F} }
d
S
{\displaystyle dS}
コーシーの公理
コーシーの公理 によれば 、応力ベクトルは、 点を通り、 において 同じ法線ベクトルを持つすべての表面に対して不変であり 、 [7] [10] すなわち、 において共通の 接線 を持つ。これは、応力ベクトルが法線ベクトルのみの関数であり 、内部表面の曲率の影響を受けないことを意味する。
T
(
n
)
{\displaystyle \mathbf {T} ^{(\mathbf {n} )}}
P
{\displaystyle P}
n
{\displaystyle \mathbf {n} }
P
{\displaystyle P}
P
{\displaystyle P}
n
{\displaystyle \mathbf {n} }
コーシーの基本補題
コーシーの公理の帰結として 、 コーシーの基本補題 [1] [7] [11] または コーシーの逆数定理 [12] : p.103–130 があり、これは同じ面の反対側に作用する応力ベクトルは大きさが等しく、方向が反対であると述べています。コーシーの基本補題は、 ニュートンの作用反作用の運動の第3法則 に相当し 、次のように表されます。
−
T
(
n
)
=
T
(
−
n
)
.
{\displaystyle -\mathbf {T} ^{(\mathbf {n} )}=\mathbf {T} ^{(-\mathbf {n} )}.}
コーシーの応力定理 - 応力テンソル
物体の ある点における応力の状態は、 その点を通るすべての平面(無限の数)に関連付けられたすべての応力ベクトル T ( n )によって定義されます。 [13] しかし、 コーシーの基本定理 [11] ( コーシーの応力定理 とも呼ばれる ) [1] によれば、互いに垂直な3つの平面上の応力ベクトルを知るだけで、その点を通る他の任意の平面上の応力ベクトルを座標変換方程式を通じて見つけることができます。
コーシーの応力定理は、コーシー応力テンソルと呼ばれる n に依存しない2 次 テンソル場 σ ( x , t ) が存在し、 T が n の線形関数となる ことを述べています 。
T
(
n
)
=
n
⋅
σ
or
T
j
(
n
)
=
σ
i
j
n
i
.
{\displaystyle \mathbf {T} ^{(\mathbf {n} )}=\mathbf {n} \cdot {\boldsymbol {\sigma }}\quad {\text{or}}\quad T_{j}^{(n)}=\sigma _{ij}n_{i}.}
この式は、法線単位ベクトルnを 持つ平面に関連付けられた連続体の任意 の点 P における応力ベクトル T ( n ) は、
座標軸に垂直な平面上の応力ベクトルの関数として、 つまり 応力テンソルσの 成分 σij として表すことができることを意味します 。
この式を証明するために、 座標平面に向いた 3 つの面と、 法線単位ベクトル nで指定された任意の方向に向いた微小面積 d Aを持つ 四面体を考えます (図 2.2)。四面体は、単位法線 n を持つ任意の平面に沿って微小要素をスライスすることによって形成されます 。この平面上の応力ベクトルは T ( n ) で表されます。四面体の面に作用する応力ベクトルは T ( e 1 ) 、 T ( e 2 ) 、および T ( e 3 ) で表され、定義により応力テンソル σ の成分 σ ij です。この四面体は、 コーシー四面体 と呼ばれることもあります。力の釣り合い、 つまり オイラーの運動の第 1 法則 (ニュートンの運動の第 2 法則) は、次のようになります。
T
(
n
)
d
A
−
T
(
e
1
)
d
A
1
−
T
(
e
2
)
d
A
2
−
T
(
e
3
)
d
A
3
=
ρ
(
h
3
d
A
)
a
,
{\displaystyle \mathbf {T} ^{(\mathbf {n} )}\,dA-\mathbf {T} ^{(\mathbf {e} _{1})}\,dA_{1}-\mathbf {T} ^{(\mathbf {e} _{2})}\,dA_{2}-\mathbf {T} ^{(\mathbf {e} _{3})}\,dA_{3}=\rho \left({\frac {h}{3}}dA\right)\mathbf {a} ,}
図 2.2. 法線単位ベクトルn を持つ平面に作用する応力ベクトル 。 符号規則に関する注意: 四面体は、平行六面体を任意の平面 nに沿ってスライスすることによって形成されます。したがって、平面 n に作用する力は 、平行六面体のもう半分によって及ぼされる反作用であり、反対の符号を持ちます。
ここで、右辺は四面体に囲まれた質量と加速度の積を表します。ρ は 密度、 a は加速度、 h は 平面 n を底辺とした場合の四面体の高さです。軸に垂直な四面体の面の面積は、d A を 各面に投影することで求められます (ドット積を使用)。
d
A
1
=
(
n
⋅
e
1
)
d
A
=
n
1
d
A
,
{\displaystyle dA_{1}=\left(\mathbf {n} \cdot \mathbf {e} _{1}\right)dA=n_{1}\;dA,}
d
A
2
=
(
n
⋅
e
2
)
d
A
=
n
2
d
A
,
{\displaystyle dA_{2}=\left(\mathbf {n} \cdot \mathbf {e} _{2}\right)dA=n_{2}\;dA,}
d
A
3
=
(
n
⋅
e
3
)
d
A
=
n
3
d
A
,
{\displaystyle dA_{3}=\left(\mathbf {n} \cdot \mathbf {e} _{3}\right)dA=n_{3}\;dA,}
そして、次の式に代入して d A を 打ち消します。
T
(
n
)
−
T
(
e
1
)
n
1
−
T
(
e
2
)
n
2
−
T
(
e
3
)
n
3
=
ρ
(
h
3
)
a
.
{\displaystyle \mathbf {T} ^{(\mathbf {n} )}-\mathbf {T} ^{(\mathbf {e} _{1})}n_{1}-\mathbf {T} ^{(\mathbf {e} _{2})}n_{2}-\mathbf {T} ^{(\mathbf {e} _{3})}n_{3}=\rho \left({\frac {h}{3}}\right)\mathbf {a} .}
四面体が点に縮む極限ケースを考えると、 hは 0にならなければなりません(直感的には、平面 nは nに沿って O に向かっ て平行移動します )。その結果、方程式の右辺は0に近づき、
T
(
n
)
=
T
(
e
1
)
n
1
+
T
(
e
2
)
n
2
+
T
(
e
3
)
n
3
.
{\displaystyle \mathbf {T} ^{(\mathbf {n} )}=\mathbf {T} ^{(\mathbf {e} _{1})}n_{1}+\mathbf {T} ^{(\mathbf {e} _{2})}n_{2}+\mathbf {T} ^{(\mathbf {e} _{3})}n_{3}.}
直交座標系の座標軸に垂直な平面を持つ材料要素(ページ上部の図を参照)を想定すると、各要素平面に関連付けられた応力ベクトル、 つまり T ( e 1 ) 、 T ( e 2 ) 、および T ( e 3 )は、法線成分と 2 つのせん断成分、 つまり 3 つの座標軸方向の成分に分解できます。法線 単位ベクトルが x 1 軸方向に向いている 表面の特定のケースでは 、法線応力を σ 11 で表し、2 つのせん断応力を σ 12 と σ 13 で表します。
T
(
e
1
)
=
T
1
(
e
1
)
e
1
+
T
2
(
e
1
)
e
2
+
T
3
(
e
1
)
e
3
=
σ
11
e
1
+
σ
12
e
2
+
σ
13
e
3
,
{\displaystyle \mathbf {T} ^{(\mathbf {e} _{1})}=T_{1}^{(\mathbf {e} _{1})}\mathbf {e} _{1}+T_{2}^{(\mathbf {e} _{1})}\mathbf {e} _{2}+T_{3}^{(\mathbf {e} _{1})}\mathbf {e} _{3}=\sigma _{11}\mathbf {e} _{1}+\sigma _{12}\mathbf {e} _{2}+\sigma _{13}\mathbf {e} _{3},}
T
(
e
2
)
=
T
1
(
e
2
)
e
1
+
T
2
(
e
2
)
e
2
+
T
3
(
e
2
)
e
3
=
σ
21
e
1
+
σ
22
e
2
+
σ
23
e
3
,
{\displaystyle \mathbf {T} ^{(\mathbf {e} _{2})}=T_{1}^{(\mathbf {e} _{2})}\mathbf {e} _{1}+T_{2}^{(\mathbf {e} _{2})}\mathbf {e} _{2}+T_{3}^{(\mathbf {e} _{2})}\mathbf {e} _{3}=\sigma _{21}\mathbf {e} _{1}+\sigma _{22}\mathbf {e} _{2}+\sigma _{23}\mathbf {e} _{3},}
T
(
e
3
)
=
T
1
(
e
3
)
e
1
+
T
2
(
e
3
)
e
2
+
T
3
(
e
3
)
e
3
=
σ
31
e
1
+
σ
32
e
2
+
σ
33
e
3
,
{\displaystyle \mathbf {T} ^{(\mathbf {e} _{3})}=T_{1}^{(\mathbf {e} _{3})}\mathbf {e} _{1}+T_{2}^{(\mathbf {e} _{3})}\mathbf {e} _{2}+T_{3}^{(\mathbf {e} _{3})}\mathbf {e} _{3}=\sigma _{31}\mathbf {e} _{1}+\sigma _{32}\mathbf {e} _{2}+\sigma _{33}\mathbf {e} _{3},}
指数表記ではこれは
T
(
e
i
)
=
T
j
(
e
i
)
e
j
=
σ
i
j
e
j
.
{\displaystyle \mathbf {T} ^{(\mathbf {e} _{i})}=T_{j}^{(\mathbf {e} _{i})}\mathbf {e} _{j}=\sigma _{ij}\mathbf {e} _{j}.}
応力ベクトルの9つの成分 σ ij は 、コーシー応力テンソル と呼ばれる2次直交座標テンソルの成分であり 、これはある点における応力の状態を完全に定義するのに使用でき、次のように表されます。
σ
=
σ
i
j
=
[
T
(
e
1
)
T
(
e
2
)
T
(
e
3
)
]
=
[
σ
11
σ
12
σ
13
σ
21
σ
22
σ
23
σ
31
σ
32
σ
33
]
≡
[
σ
x
x
σ
x
y
σ
x
z
σ
y
x
σ
y
y
σ
y
z
σ
z
x
σ
z
y
σ
z
z
]
≡
[
σ
x
τ
x
y
τ
x
z
τ
y
x
σ
y
τ
y
z
τ
z
x
τ
z
y
σ
z
]
,
{\displaystyle {\boldsymbol {\sigma }}=\sigma _{ij}=\left[{\begin{matrix}\mathbf {T} ^{(\mathbf {e} _{1})}\\\mathbf {T} ^{(\mathbf {e} _{2})}\\\mathbf {T} ^{(\mathbf {e} _{3})}\\\end{matrix}}\right]=\left[{\begin{matrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\\\end{matrix}}\right]\equiv \left[{\begin{matrix}\sigma _{xx}&\sigma _{xy}&\sigma _{xz}\\\sigma _{yx}&\sigma _{yy}&\sigma _{yz}\\\sigma _{zx}&\sigma _{zy}&\sigma _{zz}\\\end{matrix}}\right]\equiv \left[{\begin{matrix}\sigma _{x}&\tau _{xy}&\tau _{xz}\\\tau _{yx}&\sigma _{y}&\tau _{yz}\\\tau _{zx}&\tau _{zy}&\sigma _{z}\\\end{matrix}}\right],}
ここで、 σ 11 、 σ 22 、および σ 33 は 垂直応力、 σ 12 、 σ 13 、 σ 21 、 σ 23 、 σ 31 、および σ 32 はせん断応力です。最初のインデックス i は、応力が X i 軸 に垂直な平面に作用することを示し、2 番目のインデックス j は 応力が作用する方向を示します (たとえば、 σ 12 は、応力が 1 番目の 軸、つまり X 1 に垂直な平面に作用し、2 番目 の軸、つまり X 2 に沿って作用することを意味します )。応力成分は、座標軸の正の方向に作用し、作用する平面の外向きの法線ベクトルが正の座標方向を指している場合に正になります。
したがって、応力テンソルの成分を使用すると
T
(
n
)
=
T
(
e
1
)
n
1
+
T
(
e
2
)
n
2
+
T
(
e
3
)
n
3
=
∑
i
=
1
3
T
(
e
i
)
n
i
=
(
σ
i
j
e
j
)
n
i
=
σ
i
j
n
i
e
j
{\displaystyle {\begin{aligned}\mathbf {T} ^{(\mathbf {n} )}&=\mathbf {T} ^{(\mathbf {e} _{1})}n_{1}+\mathbf {T} ^{(\mathbf {e} _{2})}n_{2}+\mathbf {T} ^{(\mathbf {e} _{3})}n_{3}\\&=\sum _{i=1}^{3}\mathbf {T} ^{(\mathbf {e} _{i})}n_{i}\\&=\left(\sigma _{ij}\mathbf {e} _{j}\right)n_{i}\\&=\sigma _{ij}n_{i}\mathbf {e} _{j}\end{aligned}}}
あるいは、同等に、
T
j
(
n
)
=
σ
i
j
n
i
.
{\displaystyle T_{j}^{(\mathbf {n} )}=\sigma _{ij}n_{i}.}
あるいは、行列形式では、
[
T
1
(
n
)
T
2
(
n
)
T
3
(
n
)
]
=
[
n
1
n
2
n
3
]
⋅
[
σ
11
σ
12
σ
13
σ
21
σ
22
σ
23
σ
31
σ
32
σ
33
]
.
{\displaystyle \left[{\begin{matrix}T_{1}^{(\mathbf {n} )}&T_{2}^{(\mathbf {n} )}&T_{3}^{(\mathbf {n} )}\end{matrix}}\right]=\left[{\begin{matrix}n_{1}&n_{2}&n_{3}\end{matrix}}\right]\cdot \left[{\begin{matrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\\\end{matrix}}\right].}
コーシー応力テンソルの Voigt 表記法による表現では、応力テンソルの対称性を利用して 、 応力 を 次の形式の 6 次元ベクトルとして表現します。
σ
=
[
σ
1
σ
2
σ
3
σ
4
σ
5
σ
6
]
T
≡
[
σ
11
σ
22
σ
33
σ
23
σ
13
σ
12
]
T
.
{\displaystyle {\boldsymbol {\sigma }}={\begin{bmatrix}\sigma _{1}&\sigma _{2}&\sigma _{3}&\sigma _{4}&\sigma _{5}&\sigma _{6}\end{bmatrix}}^{\textsf {T}}\equiv {\begin{bmatrix}\sigma _{11}&\sigma _{22}&\sigma _{33}&\sigma _{23}&\sigma _{13}&\sigma _{12}\end{bmatrix}}^{\textsf {T}}.}
Voigt 表記法は、固体力学における応力とひずみの関係を表す際に、また数値構造力学ソフトウェアにおける計算効率を高めるために広く使用されています。
応力テンソルは反変 2 次テンソル であることが示されており、これは座標系の変更によってどのように変換されるかを示しています。 x i システムからx i ' システムへ 、初期システムのコンポーネント σ ij は 、テンソル変換規則に従って新しいシステムの
コンポーネント σ ij 'に変換されます (図 2.4)。
σ
i
j
′
=
a
i
m
a
j
n
σ
m
n
or
σ
′
=
A
σ
A
T
,
{\displaystyle \sigma '_{ij}=a_{im}a_{jn}\sigma _{mn}\quad {\text{or}}\quad {\boldsymbol {\sigma }}'=\mathbf {A} {\boldsymbol {\sigma }}\mathbf {A} ^{\textsf {T}},}
ここで Aは a ij の 成分を持つ 回転行列 である 。行列形式では、これは
[
σ
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 \left[{\begin{matrix}\sigma '_{11}&\sigma '_{12}&\sigma '_{13}\\\sigma '_{21}&\sigma '_{22}&\sigma '_{23}\\\sigma '_{31}&\sigma '_{32}&\sigma '_{33}\\\end{matrix}}\right]=\left[{\begin{matrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\\\end{matrix}}\right]\left[{\begin{matrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\\\end{matrix}}\right]\left[{\begin{matrix}a_{11}&a_{21}&a_{31}\\a_{12}&a_{22}&a_{32}\\a_{13}&a_{23}&a_{33}\\\end{matrix}}\right].}
図2.4 応力テンソルの変換
行列演算 を展開し 、 応力テンソルの対称性 を用いて項を簡略化すると、次の式が得られる。
σ
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
,
σ
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
,
σ
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
,
σ
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
,
σ
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
,
σ
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 {\begin{aligned}\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}.\end{aligned}}}
応力のモール 円は 、応力の変換をグラフィカルに表現したものです。
垂直応力とせん断応力
任意の平面上で法線単位ベクトルn を持つ任意の点に作用する任意 の 応力ベクトル T ( n ) の法線応力成分 σ n の大きさは、 応力テンソル σの成分 σ ij に関して 、 応力ベクトルと法線単位ベクトルのドット積
で表されます。
σ
n
=
T
(
n
)
⋅
n
=
T
i
(
n
)
n
i
=
σ
i
j
n
i
n
j
.
{\displaystyle {\begin{aligned}\sigma _{\mathrm {n} }&=\mathbf {T} ^{(\mathbf {n} )}\cdot \mathbf {n} \\&=T_{i}^{(\mathbf {n} )}n_{i}\\&=\sigma _{ij}n_{i}n_{j}.\end{aligned}}}
ベクトル n に直交するせん断応力成分 τ nの大きさは、 ピタゴラスの定理 を使用して求めることができます 。
τ
n
=
(
T
(
n
)
)
2
−
σ
n
2
=
T
i
(
n
)
T
i
(
n
)
−
σ
n
2
,
{\displaystyle {\begin{aligned}\tau _{\mathrm {n} }&={\sqrt {\left(T^{(\mathbf {n} )}\right)^{2}-\sigma _{\mathrm {n} }^{2}}}\\&={\sqrt {T_{i}^{(\mathbf {n} )}T_{i}^{(\mathbf {n} )}-\sigma _{\mathrm {n} }^{2}}},\end{aligned}}}
どこ
(
T
(
n
)
)
2
=
T
i
(
n
)
T
i
(
n
)
=
(
σ
i
j
n
j
)
(
σ
i
k
n
k
)
=
σ
i
j
σ
i
k
n
j
n
k
.
{\displaystyle \left(T^{(\mathbf {n} )}\right)^{2}=T_{i}^{(\mathbf {n} )}T_{i}^{(\mathbf {n} )}=\left(\sigma _{ij}n_{j}\right)\left(\sigma _{ik}n_{k}\right)=\sigma _{ij}\sigma _{ik}n_{j}n_{k}.}
バランスの法則 – コーシーの運動方程式
図4. 平衡状態にある連続体
コーシーの運動の第一法則
線形運動量保存 の原理によれば 、連続体が静的平衡状態にある場合、物体内のすべての質点におけるコーシー応力テンソルの成分が平衡方程式を満たすことが証明できます。
σ
j
i
,
j
+
F
i
=
0
{\displaystyle \sigma _{ji,j}+F_{i}=0}
、
どこ
σ
j
i
,
j
=
∑
j
∂
j
σ
j
i
{\displaystyle \sigma _{ji,j}=\sum _{j}\partial _{j}\sigma _{ji}}
たとえば、平衡状態の 静水圧流体 の場合、応力テンソルは次の形式になります。
σ
i
j
=
−
p
δ
i
j
,
{\displaystyle {\sigma _{ij}}=-p{\delta _{ij}},}
ここで は静水圧、 は クロネッカーのデルタ です 。
p
{\displaystyle p}
δ
i
j
{\displaystyle {\delta _{ij}}\ }
コーシーの運動の第二法則
角運動量保存 の原理によれば、平衡では 任意の点に対する モーメント の合計がゼロになる必要があり、応力テンソルは 対称的で あるという結論に至り、したがって、元の 9 つの独立した応力成分ではなく、6 つの独立した応力成分のみを持つことになります。
σ
i
j
=
σ
j
i
{\displaystyle \sigma _{ij}=\sigma _{ji}}
しかし、偶応力、つまり単位体積あたりのモーメントが存在する場合、応力テンソルは非対称になります。これは、 クヌーセン数が 1 に近い場合、または連続体が非ニュートン流体である場合にも当てはまり、 ポリマー などの回転不変でない流体につながる可能性があります 。
K
n
→
1
{\displaystyle K_{n}\rightarrow 1}
主応力と応力不変量
2D 回転要素の応力成分 。長方形要素の面 (エッジ) の方向の角度が変化すると、応力成分がどのように変化するかを示す例。 主応力は、 せん断応力がすべての面から同時に消えたときに発生します。これが発生する方向によって、 主方向 が決まります。この例では、長方形が水平の場合、応力は次のように表されます。
[
σ
11
σ
12
σ
21
σ
22
]
=
[
−
10
10
10
15
]
.
{\displaystyle \left[{\begin{matrix}\sigma _{11}&\sigma _{12}\\\sigma _{21}&\sigma _{22}\end{matrix}}\right]=\left[{\begin{matrix}-10&10\\10&15\end{matrix}}\right].}
応力を受ける物体の各点には、主平面 と呼ばれる少なくとも 3 つの平面と、 主方向 と呼ばれる 法線ベクトルがあり 、対応する応力ベクトルは平面に対して垂直、つまり法線ベクトルと平行または同じ方向であり 、法線せん断応力は存在しません。これらの主平面に対して垂直な 3 つの応力は、 主応力 と呼ばれます 。
n
{\displaystyle \mathbf {n} }
n
{\displaystyle \mathbf {n} }
τ
n
{\displaystyle \tau _{\mathrm {n} }}
応力テンソルの成分は 、検討中の点における座標系の方向に依存します。ただし、応力テンソル自体は物理量であるため、それを表すために選択された座標系とは無関係です。すべてのテンソルに関連付けられた 不変量 は、座標系からも独立しています。たとえば、ベクトルは階数 1 の単純なテンソルです。3 次元では、3 つの成分があります。これらの成分の値は、ベクトルを表すために選択された座標系に依存しますが、ベクトルの 大きさは 物理量 (スカラー) であり、ベクトルを表すために選択された 直交座標系とは無関係です ( 正規 である限り )。同様に、2 階数のテンソル (応力テンソルやひずみテンソルなど) には、3 つの独立した不変量が関連付けられています。このような不変量の 1 セットは、応力テンソルの主応力で、応力テンソルの固有値にすぎません。その方向ベクトルは、主方向または 固有ベクトル です。
σ
i
j
{\displaystyle \sigma _{ij}}
法線単位ベクトルに平行な応力ベクトルは 次のように表されます。
n
{\displaystyle \mathbf {n} }
T
(
n
)
=
λ
n
=
σ
n
n
{\displaystyle \mathbf {T} ^{(\mathbf {n} )}=\lambda \mathbf {n} =\mathbf {\sigma } _{\mathrm {n} }\mathbf {n} }
ここで、 は比例定数であり、この特定のケースでは、 法線応力ベクトルまたは主応力の大きさに対応します。
λ
{\displaystyle \lambda }
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
とを知っている ので 、
T
i
(
n
)
=
σ
i
j
n
j
{\displaystyle T_{i}^{(n)}=\sigma _{ij}n_{j}}
n
i
=
δ
i
j
n
j
{\displaystyle n_{i}=\delta _{ij}n_{j}}
T
i
(
n
)
=
λ
n
i
σ
i
j
n
j
=
λ
n
i
σ
i
j
n
j
−
λ
n
i
=
0
(
σ
i
j
−
λ
δ
i
j
)
n
j
=
0
{\displaystyle {\begin{aligned}T_{i}^{(n)}&=\lambda n_{i}\\\sigma _{ij}n_{j}&=\lambda n_{i}\\\sigma _{ij}n_{j}-\lambda n_{i}&=0\\\left(\sigma _{ij}-\lambda \delta _{ij}\right)n_{j}&=0\\\end{aligned}}}
これは、 未知数である 3つの線形方程式の 同次システム 、つまりゼロです。 の非自明な(ゼロでない)解を得るには、係数の行列式行列がゼロに等しくなければなりません。つまり、システムは特異です。したがって、
n
j
{\displaystyle n_{j}}
n
j
{\displaystyle n_{j}}
|
σ
i
j
−
λ
δ
i
j
|
=
|
σ
11
−
λ
σ
12
σ
13
σ
21
σ
22
−
λ
σ
23
σ
31
σ
32
σ
33
−
λ
|
=
0
{\displaystyle \left|\sigma _{ij}-\lambda \delta _{ij}\right|={\begin{vmatrix}\sigma _{11}-\lambda &\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}-\lambda &\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}-\lambda \\\end{vmatrix}}=0}
行列式を展開すると特性方程式 が得られる。
|
σ
i
j
−
λ
δ
i
j
|
=
−
λ
3
+
I
1
λ
2
−
I
2
λ
+
I
3
=
0
{\displaystyle \left|\sigma _{ij}-\lambda \delta _{ij}\right|=-\lambda ^{3}+I_{1}\lambda ^{2}-I_{2}\lambda +I_{3}=0}
どこ
I
1
=
σ
11
+
σ
22
+
σ
33
=
σ
k
k
=
tr
(
σ
)
I
2
=
|
σ
22
σ
23
σ
32
σ
33
|
+
|
σ
11
σ
13
σ
31
σ
33
|
+
|
σ
11
σ
12
σ
21
σ
22
|
=
σ
11
σ
22
+
σ
22
σ
33
+
σ
11
σ
33
−
σ
12
2
−
σ
23
2
−
σ
31
2
=
1
2
(
σ
i
i
σ
j
j
−
σ
i
j
σ
j
i
)
=
1
2
[
(
tr
(
σ
)
)
2
−
tr
(
σ
2
)
]
I
3
=
det
(
σ
i
j
)
=
det
(
σ
)
=
σ
11
σ
22
σ
33
+
2
σ
12
σ
23
σ
31
−
σ
12
2
σ
33
−
σ
23
2
σ
11
−
σ
31
2
σ
22
{\displaystyle {\begin{aligned}I_{1}&=\sigma _{11}+\sigma _{22}+\sigma _{33}\\&=\sigma _{kk}={\text{tr}}({\boldsymbol {\sigma }})\\[4pt]I_{2}&={\begin{vmatrix}\sigma _{22}&\sigma _{23}\\\sigma _{32}&\sigma _{33}\\\end{vmatrix}}+{\begin{vmatrix}\sigma _{11}&\sigma _{13}\\\sigma _{31}&\sigma _{33}\\\end{vmatrix}}+{\begin{vmatrix}\sigma _{11}&\sigma _{12}\\\sigma _{21}&\sigma _{22}\\\end{vmatrix}}\\&=\sigma _{11}\sigma _{22}+\sigma _{22}\sigma _{33}+\sigma _{11}\sigma _{33}-\sigma _{12}^{2}-\sigma _{23}^{2}-\sigma _{31}^{2}\\&={\frac {1}{2}}\left(\sigma _{ii}\sigma _{jj}-\sigma _{ij}\sigma _{ji}\right)={\frac {1}{2}}\left[\left({\text{tr}}({\boldsymbol {\sigma }})\right)^{2}-{\text{tr}}\left({\boldsymbol {\sigma }}^{2}\right)\right]\\[4pt]I_{3}&=\det(\sigma _{ij})=\det({\boldsymbol {\sigma }})\\&=\sigma _{11}\sigma _{22}\sigma _{33}+2\sigma _{12}\sigma _{23}\sigma _{31}-\sigma _{12}^{2}\sigma _{33}-\sigma _{23}^{2}\sigma _{11}-\sigma _{31}^{2}\sigma _{22}\\\end{aligned}}}
特性方程式には、 応力テンソルの対称性により虚数ではない、3 つの実根があります。 、、 は 主応力で、固有値の関数です。固有値は 特性多項式 の根です 。主応力は、特定の応力テンソルに対して一意です。したがって、特性方程式から、係数 、、 はそれぞれ 第 1、第 2、第 3 の 応力不変量 と呼ばれ、座標系の方向に関係なく常に同じ値になります。
λ
i
{\displaystyle \lambda _{i}}
σ
1
=
max
(
λ
1
,
λ
2
,
λ
3
)
{\displaystyle \sigma _{1}=\max \left(\lambda _{1},\lambda _{2},\lambda _{3}\right)}
σ
3
=
min
(
λ
1
,
λ
2
,
λ
3
)
{\displaystyle \sigma _{3}=\min \left(\lambda _{1},\lambda _{2},\lambda _{3}\right)}
σ
2
=
I
1
−
σ
1
−
σ
3
{\displaystyle \sigma _{2}=I_{1}-\sigma _{1}-\sigma _{3}}
λ
i
{\displaystyle \lambda _{i}}
I
1
{\displaystyle I_{1}}
I
2
{\displaystyle I_{2}}
I
3
{\displaystyle I_{3}}
各固有値に対して、方程式 の 非自明な解が存在します。これらの解は、主応力が作用する平面を定義する主方向または 固有ベクトル です 。主応力と主方向は、ある点における応力を特徴付けるものであり、方向とは無関係です。
n
j
{\displaystyle n_{j}}
(
σ
i
j
−
λ
δ
i
j
)
n
j
=
0
{\displaystyle \left(\sigma _{ij}-\lambda \delta _{ij}\right)n_{j}=0}
主方向に軸が向けられた座標系は、法線応力が主応力であり、応力テンソルが対角行列で表されることを意味します。
σ
i
j
=
[
σ
1
0
0
0
σ
2
0
0
0
σ
3
]
{\displaystyle \sigma _{ij}={\begin{bmatrix}\sigma _{1}&0&0\\0&\sigma _{2}&0\\0&0&\sigma _{3}\end{bmatrix}}}
主応力は、応力不変量、、、およびを形成するために組み合わせることができる 。 最初 の不変量と3番目の不変量は、それぞれ応力テンソルのトレースと行列式である。したがって、
I
1
{\displaystyle I_{1}}
I
2
{\displaystyle I_{2}}
I
3
{\displaystyle I_{3}}
I
1
=
σ
1
+
σ
2
+
σ
3
I
2
=
σ
1
σ
2
+
σ
2
σ
3
+
σ
3
σ
1
I
3
=
σ
1
σ
2
σ
3
{\displaystyle {\begin{aligned}I_{1}&=\sigma _{1}+\sigma _{2}+\sigma _{3}\\I_{2}&=\sigma _{1}\sigma _{2}+\sigma _{2}\sigma _{3}+\sigma _{3}\sigma _{1}\\I_{3}&=\sigma _{1}\sigma _{2}\sigma _{3}\\\end{aligned}}}
主座標系は単純であるため、特定の点における弾性媒体の状態を検討するときに役立ちます。主応力は、x方向とy方向の応力、または部品の軸方向と曲げ応力を評価するために、次の式で表されることがよくあります。 [14] :p.58–59 主法線応力を使用して、 フォンミーゼス応力 を計算し、最終的には安全係数と安全余裕を計算できます。
σ
1
,
σ
2
=
σ
x
+
σ
y
2
±
(
σ
x
−
σ
y
2
)
2
+
τ
x
y
2
{\displaystyle \sigma _{1},\sigma _{2}={\frac {\sigma _{x}+\sigma _{y}}{2}}\pm {\sqrt {\left({\frac {\sigma _{x}-\sigma _{y}}{2}}\right)^{2}+\tau _{xy}^{2}}}}
方程式の平方根 の下の部分だけを使用すると、 プラスとマイナスの最大せん断応力と最小せん断応力に等しくなります。これは次のように示されます。
τ
max
,
τ
min
=
±
(
σ
x
−
σ
y
2
)
2
+
τ
x
y
2
{\displaystyle \tau _{\max },\tau _{\min }=\pm {\sqrt {\left({\frac {\sigma _{x}-\sigma _{y}}{2}}\right)^{2}+\tau _{xy}^{2}}}}
最大および最小せん断応力
最大せん断応力または最大主せん断応力は、最大主応力と最小主応力の差の半分に等しく、最大主応力と最小主応力の方向の間の角度を二等分する平面に作用します。つまり、最大せん断応力の平面は 主応力平面から方向付けられます。最大せん断応力は次のように表されます。
45
∘
{\displaystyle 45^{\circ }}
τ
max
=
1
2
|
σ
max
−
σ
min
|
{\displaystyle \tau _{\max }={\frac {1}{2}}\left|\sigma _{\max }-\sigma _{\min }\right|}
そうすると
σ
1
≥
σ
2
≥
σ
3
{\displaystyle \sigma _{1}\geq \sigma _{2}\geq \sigma _{3}}
τ
max
=
1
2
|
σ
1
−
σ
3
|
{\displaystyle \tau _{\max }={\frac {1}{2}}\left|\sigma _{1}-\sigma _{3}\right|}
応力テンソルがゼロでない場合、最大せん断応力に対して平面に作用する法線応力成分はゼロでなく、次の式に等しい。
σ
n
=
1
2
(
σ
1
+
σ
3
)
{\displaystyle \sigma _{\text{n}}={\frac {1}{2}}\left(\sigma _{1}+\sigma _{3}\right)}
応力偏差テンソル
応力テンソルは、 他の 2 つの応力テンソルの合計として表すことができます。
σ
i
j
{\displaystyle \sigma _{ij}}
平均 静水圧応力 テンソル または 体積応力テンソル または 平均法線応力テンソル 、 これは応力を受ける物体の体積を変化させる傾向がある 。
π
δ
i
j
{\displaystyle \pi \delta _{ij}}
応力偏差テンソル と呼ばれる偏差成分は 、 それを歪める傾向があります。
s
i
j
{\displaystyle s_{ij}}
それで
σ
i
j
=
s
i
j
+
π
δ
i
j
,
{\displaystyle \sigma _{ij}=s_{ij}+\pi \delta _{ij},\,}
ここで 平均応力は次のように与えられる。
π
{\displaystyle \pi }
π
=
σ
k
k
3
=
σ
11
+
σ
22
+
σ
33
3
=
1
3
I
1
.
{\displaystyle \pi ={\frac {\sigma _{kk}}{3}}={\frac {\sigma _{11}+\sigma _{22}+\sigma _{33}}{3}}={\frac {1}{3}}I_{1}.\,}
圧力 ( )は、一般に、応力テンソルの トレース の負の1/3から 速度の発散が寄与する応力を引いたものとして定義されます。つまり、
p
{\displaystyle p}
p
=
λ
∇
⋅
u
→
−
π
=
λ
∂
u
k
∂
x
k
−
π
=
∑
k
λ
∂
u
k
∂
x
k
−
π
,
{\displaystyle p=\lambda \,\nabla \cdot {\vec {u}}-\pi =\lambda \,{\frac {\partial u_{k}}{\partial x_{k}}}-\pi =\sum _{k}\lambda \,{\frac {\partial u_{k}}{\partial x_{k}}}-\pi ,}
ここで、 は比例定数(つまり、 ラメパラメータ の最初のもの)、 は 発散演算子 、 は k 番目の 直交座標 、 は 流速 、は の k 番目の直交成分 です 。
λ
{\displaystyle \lambda }
∇
⋅
{\displaystyle \nabla \cdot }
x
k
{\displaystyle x_{k}}
u
→
{\displaystyle {\vec {u}}}
u
k
{\displaystyle u_{k}}
u
→
{\displaystyle {\vec {u}}}
偏差応力テンソルは、コーシー応力テンソルから静水圧応力テンソルを減算することで得られます。
s
i
j
=
σ
i
j
−
σ
k
k
3
δ
i
j
,
[
s
11
s
12
s
13
s
21
s
22
s
23
s
31
s
32
s
33
]
=
[
σ
11
σ
12
σ
13
σ
21
σ
22
σ
23
σ
31
σ
32
σ
33
]
−
[
π
0
0
0
π
0
0
0
π
]
=
[
σ
11
−
π
σ
12
σ
13
σ
21
σ
22
−
π
σ
23
σ
31
σ
32
σ
33
−
π
]
.
{\displaystyle {\begin{aligned}s_{ij}&=\sigma _{ij}-{\frac {\sigma _{kk}}{3}}\delta _{ij},\,\\\left[{\begin{matrix}s_{11}&s_{12}&s_{13}\\s_{21}&s_{22}&s_{23}\\s_{31}&s_{32}&s_{33}\end{matrix}}\right]&=\left[{\begin{matrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\end{matrix}}\right]-\left[{\begin{matrix}\pi &0&0\\0&\pi &0\\0&0&\pi \end{matrix}}\right]\\&=\left[{\begin{matrix}\sigma _{11}-\pi &\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}-\pi &\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}-\pi \end{matrix}}\right].\end{aligned}}}
応力偏差テンソルの不変量
応力偏差テンソルは2次テンソルなので、不変量 の集合も持ち、応力テンソルの不変量を計算するのと同じ手順で得ることができます。応力偏差テンソルの主方向は 、応力テンソルの主方向と同じであること が示されます 。したがって、特性方程式は次のようになります。
s
i
j
{\displaystyle s_{ij}}
σ
i
j
{\displaystyle \sigma _{ij}}
|
s
i
j
−
λ
δ
i
j
|
=
λ
3
−
J
1
λ
2
−
J
2
λ
−
J
3
=
0
,
{\displaystyle \left|s_{ij}-\lambda \delta _{ij}\right|=\lambda ^{3}-J_{1}\lambda ^{2}-J_{2}\lambda -J_{3}=0,}
ここで 、、 およびは、それぞれ 第 1、第 2、および第 3 偏差応力不変量 です。これらの値は、選択した座標系の方向に関係なく同じ (不変) です。これらの偏差応力不変量は、 の成分 またはその主値 、、 およびの関数として表すことができます。または 、 またはその主値 、、 および の関数として表すことができます 。したがって、
J
1
{\displaystyle J_{1}}
J
2
{\displaystyle J_{2}}
J
3
{\displaystyle J_{3}}
s
i
j
{\displaystyle s_{ij}}
s
1
{\displaystyle s_{1}}
s
2
{\displaystyle s_{2}}
s
3
{\displaystyle s_{3}}
σ
i
j
{\displaystyle \sigma _{ij}}
σ
1
{\displaystyle \sigma _{1}}
σ
2
{\displaystyle \sigma _{2}}
σ
3
{\displaystyle \sigma _{3}}
J
1
=
s
k
k
=
0
,
J
2
=
1
2
s
i
j
s
j
i
=
1
2
tr
(
s
2
)
=
1
2
(
s
1
2
+
s
2
2
+
s
3
2
)
=
1
6
[
(
σ
11
−
σ
22
)
2
+
(
σ
22
−
σ
33
)
2
+
(
σ
33
−
σ
11
)
2
]
+
σ
12
2
+
σ
23
2
+
σ
31
2
=
1
6
[
(
σ
1
−
σ
2
)
2
+
(
σ
2
−
σ
3
)
2
+
(
σ
3
−
σ
1
)
2
]
=
1
3
I
1
2
−
I
2
=
1
2
[
tr
(
σ
2
)
−
1
3
tr
(
σ
)
2
]
,
J
3
=
det
(
s
i
j
)
=
1
3
s
i
j
s
j
k
s
k
i
=
1
3
tr
(
s
3
)
=
1
3
(
s
1
3
+
s
2
3
+
s
3
3
)
=
s
1
s
2
s
3
=
2
27
I
1
3
−
1
3
I
1
I
2
+
I
3
=
1
3
[
tr
(
σ
3
)
−
tr
(
σ
2
)
tr
(
σ
)
+
2
9
tr
(
σ
)
3
]
.
{\displaystyle {\begin{aligned}J_{1}&=s_{kk}=0,\\[3pt]J_{2}&={\frac {1}{2}}s_{ij}s_{ji}={\frac {1}{2}}\operatorname {tr} \left({\boldsymbol {s}}^{2}\right)\\&={\frac {1}{2}}\left(s_{1}^{2}+s_{2}^{2}+s_{3}^{2}\right)\\&={\frac {1}{6}}\left[(\sigma _{11}-\sigma _{22})^{2}+(\sigma _{22}-\sigma _{33})^{2}+(\sigma _{33}-\sigma _{11})^{2}\right]+\sigma _{12}^{2}+\sigma _{23}^{2}+\sigma _{31}^{2}\\&={\frac {1}{6}}\left[(\sigma _{1}-\sigma _{2})^{2}+(\sigma _{2}-\sigma _{3})^{2}+(\sigma _{3}-\sigma _{1})^{2}\right]\\&={\frac {1}{3}}I_{1}^{2}-I_{2}={\frac {1}{2}}\left[\operatorname {tr} \left({\boldsymbol {\sigma }}^{2}\right)-{\frac {1}{3}}\operatorname {tr} ({\boldsymbol {\sigma }})^{2}\right],\\[3pt]J_{3}&=\det(s_{ij})\\&={\frac {1}{3}}s_{ij}s_{jk}s_{ki}={\frac {1}{3}}{\text{tr}}\left({\boldsymbol {s}}^{3}\right)\\&={\frac {1}{3}}\left(s_{1}^{3}+s_{2}^{3}+s_{3}^{3}\right)\\&=s_{1}s_{2}s_{3}\\&={\frac {2}{27}}I_{1}^{3}-{\frac {1}{3}}I_{1}I_{2}+I_{3}={\frac {1}{3}}\left[{\text{tr}}({\boldsymbol {\sigma }}^{3})-\operatorname {tr} \left({\boldsymbol {\sigma }}^{2}\right)\operatorname {tr} ({\boldsymbol {\sigma }})+{\frac {2}{9}}\operatorname {tr} ({\boldsymbol {\sigma }})^{3}\right].\,\end{aligned}}}
なぜなら 、応力偏差テンソルは純粋せん断状態にあるからです。
s
k
k
=
0
{\displaystyle s_{kk}=0}
等価応力またはフォン・ミーゼス応力 と呼ばれる量は 固体力学でよく使用されます。等価応力は次のように定義されます。
σ
vM
=
3
J
2
=
1
2
[
(
σ
1
−
σ
2
)
2
+
(
σ
2
−
σ
3
)
2
+
(
σ
3
−
σ
1
)
2
]
.
{\displaystyle \sigma _{\text{vM}}={\sqrt {3\,J_{2}}}={\sqrt {{\frac {1}{2}}~\left[(\sigma _{1}-\sigma _{2})^{2}+(\sigma _{2}-\sigma _{3})^{2}+(\sigma _{3}-\sigma _{1})^{2}\right]}}\,.}
八面体応力
図6. 八面体応力面
主方向を座標軸と見なして、法線ベクトルが各主軸と等しい角度をなす(つまり、方向余弦が に等しい) 平面を八面体平面 と呼びます 。八面体平面は全部で 8 個あります(図 6)。これらの平面上の応力テンソルの法線成分とせん断成分は、それぞれ 八面体法線応力 と 八面体せん断応力 と呼ばれます。原点を通る八面体平面は π 平面 と呼ばれます( π は 、上のセクションで π で示した 平均応力 と混同しないでください) 。π 平面 上では 、 です。
|
1
/
3
|
{\displaystyle |1/{\sqrt {3}}|}
σ
oct
{\displaystyle \sigma _{\text{oct}}}
τ
oct
{\displaystyle \tau _{\text{oct}}}
s
i
j
=
1
3
I
{\textstyle s_{ij}={\frac {1}{3}}I}
主軸上の点O(図6)の応力テンソルが
σ
i
j
=
[
σ
1
0
0
0
σ
2
0
0
0
σ
3
]
{\displaystyle \sigma _{ij}={\begin{bmatrix}\sigma _{1}&0&0\\0&\sigma _{2}&0\\0&0&\sigma _{3}\end{bmatrix}}}
八面体平面上の応力ベクトルは次のように表される。
T
oct
(
n
)
=
σ
i
j
n
i
e
j
=
σ
1
n
1
e
1
+
σ
2
n
2
e
2
+
σ
3
n
3
e
3
=
1
3
(
σ
1
e
1
+
σ
2
e
2
+
σ
3
e
3
)
{\displaystyle {\begin{aligned}\mathbf {T} _{\text{oct}}^{(\mathbf {n} )}&=\sigma _{ij}n_{i}\mathbf {e} _{j}\\&=\sigma _{1}n_{1}\mathbf {e} _{1}+\sigma _{2}n_{2}\mathbf {e} _{2}+\sigma _{3}n_{3}\mathbf {e} _{3}\\&={\frac {1}{\sqrt {3}}}(\sigma _{1}\mathbf {e} _{1}+\sigma _{2}\mathbf {e} _{2}+\sigma _{3}\mathbf {e} _{3})\end{aligned}}}
八面体平面に関連する点Oにおける応力ベクトルの法線成分は
σ
oct
=
T
i
(
n
)
n
i
=
σ
i
j
n
i
n
j
=
σ
1
n
1
n
1
+
σ
2
n
2
n
2
+
σ
3
n
3
n
3
=
1
3
(
σ
1
+
σ
2
+
σ
3
)
=
1
3
I
1
{\displaystyle {\begin{aligned}\sigma _{\text{oct}}&=T_{i}^{(n)}n_{i}\\&=\sigma _{ij}n_{i}n_{j}\\&=\sigma _{1}n_{1}n_{1}+\sigma _{2}n_{2}n_{2}+\sigma _{3}n_{3}n_{3}\\&={\frac {1}{3}}(\sigma _{1}+\sigma _{2}+\sigma _{3})={\frac {1}{3}}I_{1}\end{aligned}}}
これは平均法線応力または静水圧応力です。この値は8つの八面体平面すべてで同じです。八面体平面上のせん断応力は
τ
oct
=
T
i
(
n
)
T
i
(
n
)
−
σ
oct
2
=
[
1
3
(
σ
1
2
+
σ
2
2
+
σ
3
2
)
−
1
9
(
σ
1
+
σ
2
+
σ
3
)
2
]
1
2
=
1
3
[
(
σ
1
−
σ
2
)
2
+
(
σ
2
−
σ
3
)
2
+
(
σ
3
−
σ
1
)
2
]
1
2
=
1
3
2
I
1
2
−
6
I
2
=
2
3
J
2
{\displaystyle {\begin{aligned}\tau _{\text{oct}}&={\sqrt {T_{i}^{(n)}T_{i}^{(n)}-\sigma _{\text{oct}}^{2}}}\\&=\left[{\frac {1}{3}}\left(\sigma _{1}^{2}+\sigma _{2}^{2}+\sigma _{3}^{2}\right)-{\frac {1}{9}}(\sigma _{1}+\sigma _{2}+\sigma _{3})^{2}\right]^{\frac {1}{2}}\\&={\frac {1}{3}}\left[(\sigma _{1}-\sigma _{2})^{2}+(\sigma _{2}-\sigma _{3})^{2}+(\sigma _{3}-\sigma _{1})^{2}\right]^{\frac {1}{2}}={\frac {1}{3}}{\sqrt {2I_{1}^{2}-6I_{2}}}={\sqrt {{\frac {2}{3}}J_{2}}}\end{aligned}}}
参照
注記
^ 詳細:
[
T
1
(
e
)
T
2
(
e
)
T
3
(
e
)
]
=
[
e
1
e
2
e
3
]
⋅
[
σ
11
σ
12
σ
13
σ
21
σ
22
σ
23
σ
31
σ
32
σ
33
]
.
{\displaystyle \left[{\begin{matrix}T_{1}^{(\mathbf {e} )}&T_{2}^{(\mathbf {e} )}&T_{3}^{(\mathbf {e} )}\end{matrix}}\right]=\left[{\begin{matrix}e_{1}&e_{2}&e_{3}\end{matrix}}\right]\cdot \left[{\begin{matrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\\\end{matrix}}\right].}
参考文献
^ abc
Fridtjov Irgens (2008)、「連続力学」。スプリンガー。 ISBN 3-540-74297-2
^ トゥルーズデル、C . ; Toupin, RA (1960)、「The Classical Field Theories」、 Flügge, Siegfried (編)、 Principles of Classical Mechanics and Field Theory/Prinzipien der Klassischen Mechanik und Feldtheorie 、Handbuch der Physik (Encyclopedia of Physics)、vol. III/1、ベルリン – ハイデルベルク – ニューヨーク: Springer-Verlag 、pp. 226–793、 Bibcode :1960HDP....2....F、 doi :10.1007/978-3-642-45943-6 、 ISBN
978-3-540-02547-4 、 MR 0118005、 Zbl 0118.39702 。
^ Peter Chadwick (
1999)、「連続体力学: 簡潔な理論と問題」。Dover Publications、「物理学の本」シリーズ。ISBN 0-486-40180-4 。ページ
^ ab Yuan-cheng Fung および Pin Tong (2001)「古典的および計算的固体力学 」
。World Scientific。ISBN 981-02-4124-0
^ ab スミス&トゥルーズデル p.97
^ ab G. Thomas Mase および George E. Mase (1999)、「エンジニアのための連続体力学 」
(第 2 版)。CRC プレス。ISBN 0-8493-1855-6
^ abc
I-Shih Liu (2002)、「連続体力学」。Springer ISBN 3-540-43019-9
^ ab Han-Chin Wu (2005)、「連続体力学と塑性 」
。CRC Press。ISBN 1-58488-363-4
^ ルブリナー
^ バサール
^ abc
Teodor M. Atanackovic および Ardéshir Guran (2000)、「科学者とエンジニアのための弾性理論」。Springer。ISBN 0-8176-4072 -X
^
Keith D. Hjelmstad (2005)、「構造力学の基礎」(第 2 版)。Prentice-Hall。ISBN 0-387-23330 -X
^
Wai-Fah Chen および Da-Jian Han (2007)、「構造エンジニアのための可塑性」 J. Ross Publishing ISBN 1-932159-75-4
^
バーナード・ハムロック (2005)、「機械要素の基礎」、マグロウヒル、 ISBN 0-07-297682-9
^
ラビンドラナート・チャタジー (1999)、「連続体力学の数学的理論」。アルファサイエンス 。ISBN 81-7319-244-8
^ John Conrad Jaeger、NGW Cook、RW Zimmerman (2007)、「Fundamentals of Rock Mechanics 」
(第 4 版)。Wiley-Blackwell。ISBN 0-632-05759-9
^
モハメッド・アミーン (2005)、「計算弾性: 弾性理論と有限要素法および境界要素法」(書籍)。アルファサイエンス、 ISBN 1-84265-201-X
^
ウィリアム・プラガー (2004)、「連続体の力学入門」ドーバー出版 ISBN 0-486-43809-0