土木工学における幾何学的計算技術
図1. 3次元応力状態のモール円
モール円は、 コーシー応力テンソル の 変換則 を 2 次元でグラフィカルに表現したものです 。
モール円は、 機械工学 における 材料の強度 、 地盤工学 における 土壌の強度 、 構造工学 における建築物の強度に関する計算でよく使用されます。また、多くの平面における 応力を 垂直成分と水平成分に減らして計算するためにも使用されます。これらは 主応力 が計算される主平面と呼ばれます。モール円は、主平面と主応力をグラフィカルに表現するためにも使用でき、そのための最も簡単な方法の1つです。 [1]
連続体 とみなされる物質体の 応力解析 を実行すると、 座標系 に対する特定の物質点におけるコーシー応力テンソルの成分がわかります 。次に、モール円を使用して、回転した座標系、つまりその点を通過する異なる方向の平面に作用する応力成分をグラフィカルに決定します。
円上の各点の横座標 と縦座標 ( 、 ) は、回転座標系に作用する法線応力と せん断応力の 成分の大きさです。言い換えると、円は、各平面上のあらゆる方向における応力の状態を表す点の 軌跡 であり、軸は応力要素の主軸を表します。
σ
ん
{\displaystyle \sigma _{\mathrm {n} }}
τ
ん
{\displaystyle \tau _{\mathrm {n} }}
19世紀のドイツ人エンジニア、 カール・クルマンは、 曲げ加工 中の水平梁の縦応力と垂直応力を考慮した応力のグラフィカルな表現を初めて考案しました 。彼の研究は、同じドイツのエンジニアである クリスチャン・オットー・モール (応力円の由来)に影響を与え、モールはそれを2次元と3次元の両方の応力に拡張し、 応力円に基づいた 破壊基準を開発しました。 [2]
ある点における応力状態を表す別のグラフィカル手法としては、 ラメの応力楕円体 や コーシーの応力二次曲線など があります。
モール円は、 ひずみテンソル や 慣性モーメントテンソルを含む、 任意 の対称2x2 テンソル 行列に適用できます 。
モチベーション
図 2. 連続体として想定される荷重を受けた変形可能な材料本体内の応力。
内部力は、連続体 と仮定される変形可能な物体の粒子間に、加えられた 外力、すなわち 表面力 または 体積力 に対する反応として生じます。この反応は、連続体に対する オイラーの運動法則 に従います。これは、粒子に対する ニュートンの運動法則 と同等です。これらの内部 力の強さの尺度は、 応力 と呼ばれます 。物体は連続体と仮定されるため、これらの内部力は物体の体積内で連続的に分布します。
構造 工学、 機械工学 、 地質 工学 などの工学では、物体内の応力分布、たとえばトンネル周辺の岩盤、飛行機の翼、建物の柱などの応力は、 応力解析 によって判定されます。応力分布を計算するということは、物体内のすべての点(物質粒子)における応力の判定を意味します。 コーシー によれば、物体内の 任意の点(図 2)の応力は、連続体と仮定すると、 コーシー応力テンソル と呼ばれる タイプ (2,0) の 2 次 テンソル の9 つの応力成分によって完全に定義され ます
。
σ
私
じゅう
{\displaystyle \sigma_{ij}}
σ
{\displaystyle {\boldsymbol {\シグマ }}}
σ
=
[
σ
11
σ
12
σ
13
σ
21
σ
22
σ
23
σ
31
σ
32
σ
33
]
≡
[
σ
x
x
σ
x
ええ
σ
x
ず
σ
ええ
x
σ
ええ
ええ
σ
ええ
ず
σ
ず
x
σ
ず
ええ
σ
ず
ず
]
≡
[
σ
x
τ
x
ええ
τ
x
ず
τ
ええ
x
σ
ええ
τ
ええ
ず
τ
ず
x
τ
ず
ええ
σ
ず
]
{\displaystyle {\boldsymbol {\sigma }}=\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{行列}}\right]}
図 3. 平面応力条件下での連続体内の一点における応力変換。
物体内の応力分布が座標系 に関して決定された後、 回転した座標系 に関して 特定の物質点における応力テンソルの成分、すなわち 、その注目点を通過し座標系と角度を形成する異なる方向の平面に作用する応力を計算することが必要となる場合があります (図 3)。たとえば、最大垂直応力と最大せん断応力、およびそれらが作用する平面の方向を見つけることが目的です。これを実現するには、座標系を回転させながらテンソル変換を行う必要があります。 テンソル の定義から、コーシー応力テンソルは テンソル変換法則 に従います。コーシー応力テンソルのこの変換法則をグラフィカルに表現したものが、応力のモール円です。
(
x
、
ええ
)
{\displaystyle (x,y)}
ポ
{\displaystyle P}
(
x
′
、
ええ
′
)
{\displaystyle (x',y')}
(
x
、
ええ
)
{\displaystyle (x,y)}
2次元応力状態のモール円
図 4. 平面応力条件下での連続体内の点を通過する平面における応力成分。
2 次元では、任意の 2 つの垂直方向に対する特定の物質点の応力テンソルは、 3 つの応力成分のみで完全に定義されます。特定の座標系では、 これらの応力成分は、法線応力 と 、せん断応力です 。角運動量のバランスから、コーシー応力テンソルの対称性が実証されます。この対称性は、 を意味します 。したがって、コーシー応力テンソルは次のように記述できます。
ポ
{\displaystyle P}
(
x
、
ええ
)
{\displaystyle (x,y)}
σ
x
{\displaystyle \sigma_{x}}
σ
ええ
{\displaystyle \sigma_{y}}
τ
x
ええ
{\displaystyle \tau_{xy}}
τ
x
ええ
=
τ
ええ
x
{\displaystyle \tau_{xy}=\tau_{yx}}
σ
=
[
σ
x
τ
x
ええ
0
τ
x
ええ
σ
ええ
0
0
0
0
]
≡
[
σ
x
τ
x
ええ
τ
x
ええ
σ
ええ
]
{\displaystyle {\boldsymbol {\sigma }}=\left[{\begin{matrix}\sigma _{x}&\tau _{xy}&0\\\tau _{xy}&\sigma _{y}&0\\0&0&0\\\end{matrix}}\right]\equiv \left[{\begin{matrix}\sigma _{x}&\tau _{xy}\\\tau _{xy}&\sigma _{y}\\\end{matrix}}\right]}
目的は、モール円を使用して、 回転座標系、つまり、 -平面 を通り 、平面に垂直な異なる方向の平面 上の応力成分 とを見つけることです (図 4)。回転座標系は、 元の座標系 と 角度を成します 。
σ
ん
{\displaystyle \sigma _{\mathrm {n} }}
τ
ん
{\displaystyle \tau _{\mathrm {n} }}
(
x
′
、
ええ
′
)
{\displaystyle (x',y')}
ポ
{\displaystyle P}
x
{\displaystyle x}
ええ
{\displaystyle y}
(
x
′
、
ええ
′
)
{\displaystyle (x',y')}
θ
{\displaystyle \theta}
(
x
、
ええ
)
{\displaystyle (x,y)}
モール円の方程式
平面応力 と 平面ひずみ の2次元の場合のモール円の方程式を導くには、まず、平面に平行な 方向、つまりページまたは画面に垂直な
方向 に単位面積を持つ、 質点の周りの2次元微小材料要素(図4)を考えます。
ポ
{\displaystyle P}
ええ
{\displaystyle y}
ず
{\displaystyle z}
微小要素上の力の平衡から、法線応力 とせん断応力の大きさは 次のように表されます。
σ
ん
{\displaystyle \sigma _{\mathrm {n} }}
τ
ん
{\displaystyle \tau _{\mathrm {n} }}
σ
ん
=
1
2
(
σ
x
+
σ
ええ
)
+
1
2
(
σ
x
−
σ
ええ
)
コス
2
θ
+
τ
x
ええ
罪
2
θ
{\displaystyle \sigma _{\mathrm {n} }={\frac {1}{2}}(\sigma _{x}+\sigma _{y})+{\frac {1}{2}}(\sigma _{x}-\sigma _{y})\cos 2\theta +\tau _{xy}\sin 2\theta }
τ
n
=
−
1
2
(
σ
x
−
σ
y
)
sin
2
θ
+
τ
x
y
cos
2
θ
{\displaystyle \tau _{\mathrm {n} }=-{\frac {1}{2}}(\sigma _{x}-\sigma _{y})\sin 2\theta +\tau _{xy}\cos 2\theta }
両方の方程式は、既知のコーシー応力テンソルにテンソル変換則を適用することによっても得られます。これは、および の方向の力の静的平衡を実行することと同等です 。
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
これら 2 つの方程式は、 モール円の パラメトリック方程式 です。これらの方程式では、はパラメータであり、 と は座標です。つまり、横軸 と縦軸 を持つ座標系を選択し 、パラメータ に値を与えると、 得られた点が円上に配置されること
になります。
2
θ
{\displaystyle 2\theta }
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
θ
{\displaystyle \theta }
これらのパラメトリック方程式からパラメータを消去すると 、モール円の非パラメトリック方程式が得られます。これは、およびの方程式を並べ替えることで実現できます。まず 、 最初の方程式の最初の項を転置し、各方程式の両辺を2乗してから加算します。したがって、次の式が得られます。
2
θ
{\displaystyle 2\theta }
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
[
σ
n
−
1
2
(
σ
x
+
σ
y
)
]
2
+
τ
n
2
=
[
1
2
(
σ
x
−
σ
y
)
]
2
+
τ
x
y
2
(
σ
n
−
σ
a
v
g
)
2
+
τ
n
2
=
R
2
{\displaystyle {\begin{aligned}\left[\sigma _{\mathrm {n} }-{\tfrac {1}{2}}(\sigma _{x}+\sigma _{y})\right]^{2}+\tau _{\mathrm {n} }^{2}&=\left[{\tfrac {1}{2}}(\sigma _{x}-\sigma _{y})\right]^{2}+\tau _{xy}^{2}\\(\sigma _{\mathrm {n} }-\sigma _{\mathrm {avg} })^{2}+\tau _{\mathrm {n} }^{2}&=R^{2}\end{aligned}}}
どこ
R
=
[
1
2
(
σ
x
−
σ
y
)
]
2
+
τ
x
y
2
and
σ
a
v
g
=
1
2
(
σ
x
+
σ
y
)
{\displaystyle R={\sqrt {\left[{\tfrac {1}{2}}(\sigma _{x}-\sigma _{y})\right]^{2}+\tau _{xy}^{2}}}\quad {\text{and}}\quad \sigma _{\mathrm {avg} }={\tfrac {1}{2}}(\sigma _{x}+\sigma _{y})}
これは、次の形の円 (モール円)
の方程式である。
(
x
−
a
)
2
+
(
y
−
b
)
2
=
r
2
{\displaystyle (x-a)^{2}+(y-b)^{2}=r^{2}}
座標系内の 座標を持つ点を中心と する半径 。
r
=
R
{\displaystyle r=R}
(
a
,
b
)
=
(
σ
a
v
g
,
0
)
{\displaystyle (a,b)=(\sigma _{\mathrm {avg} },0)}
(
σ
n
,
τ
n
)
{\displaystyle (\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
サイン規約
モール円を使用する際に考慮する必要がある 2 つの異なる符号規則があります。1 つは「物理空間」の応力成分に対する符号規則で、もう 1 つは「モール円空間」の応力成分に対する符号規則です。さらに、2 つの符号規則のそれぞれにおいて、 工学力学 ( 構造工学 および 機械工学) の文献は 地質力学の 文献とは異なる符号規則に従います 。標準的な符号規則はなく、特定の符号規則の選択は、特定の問題に対する計算と解釈の利便性によって左右されます。これらの符号規則の詳細な説明は、以下で示します。
図 4 を使用したモール円の式の導出は、工学力学の符号規則に従っています。 この記事でも工学力学の符号規則を使用します 。
物理空間の符号規則
コーシー応力テンソルの慣例(図 3 および図 4)から、応力成分の最初の下付き文字は応力成分が作用する面を示し、2 番目の下付き文字は応力成分の方向を示します。したがって、 せん断応力は、-軸の正の方向に法線ベクトルを持つ面に作用し 、 -軸の正の方向に作用します。
τ
x
y
{\displaystyle \tau _{xy}}
x
{\displaystyle x}
y
{\displaystyle y}
物理空間の符号規則では、正の法線応力は作用面(引張)に対して外側に向かい、負の法線応力は作用面(圧縮)に対して内側に向かいます(図 5)。
物理空間の符号規則では、正のせん断応力は、軸の正の方向にある材料要素の正の面に作用します。また、正のせん断応力は、軸の負の方向にある材料要素の負の面に作用します。正の面は軸の正の方向に法線ベクトルを持ち、負の面は軸の負の方向に法線ベクトルを持ちます。たとえば、せん断応力 と は正の面に作用するため正であり、それぞれ - 軸と- 軸 の正の方向にも作用します (図 3)。同様に、負の面に作用するそれぞれの反対のせん断応力 と は、それぞれ - 軸と- 軸 の負の方向に作用するため 、負の符号を持ちます。
τ
x
y
{\displaystyle \tau _{xy}}
τ
y
x
{\displaystyle \tau _{yx}}
y
{\displaystyle y}
x
{\displaystyle x}
τ
x
y
{\displaystyle \tau _{xy}}
τ
y
x
{\displaystyle \tau _{yx}}
x
{\displaystyle x}
y
{\displaystyle y}
モール円空間記号規則
図 5. モール円を描くための工学力学の符号規則。この記事では、図に示すように符号規則 # 3 に従います。
モール円空間符号規則では、法線応力は物理空間符号規則の法線応力と同じ符号を持ちます。つまり、正の法線応力は作用面の外側に作用し、負の法線応力は作用面の内側に作用します。
ただし、モール円空間におけるせん断応力の慣習は、物理空間における慣習とは異なります。モール円空間の符号慣習では、正のせん断応力は材料要素を反時計回りに回転させ、負のせん断応力は材料を時計回りに回転させます。このように、モール円 空間ではせん断応力成分は正であり、 モール円空間ではせん断応力成分は負です。
τ
x
y
{\displaystyle \tau _{xy}}
τ
y
x
{\displaystyle \tau _{yx}}
モール円空間を描画するには、数学的に正しいモール円を生成する 2 つのオプションがあります。
正のせん断応力は上向きにプロットされます(図5、符号規則#1)
正のせん断応力は下向きにプロットされ、つまり - 軸が反転されます (図 5、符号規則 #2)。
τ
n
{\displaystyle \tau _{\mathrm {n} }}
正のせん断応力を上向きにプロットすると、モール円上の角度は 時計回りに正に回転しますが、これは物理空間の慣例とは逆です。そのため、一部の著者 [3] は、せん断応力の物理空間の慣例と同様に、モール円上の角度が 反時計回りに正に回転するように、正のせん断応力を下向きにプロットすることを好みます。
2
θ
{\displaystyle 2\theta }
2
θ
{\displaystyle 2\theta }
モール円空間でせん断応力軸が下向きになるという「問題」を克服するために、正のせん断応力は材料要素を時計回りに回転させ、負のせん断応力は材料要素を反時計回りに回転させるという 別の 符号規則があります (図 5、オプション 3)。この方法では、正のせん断応力はモール円空間で上向きにプロットされ、角度は モール円空間で反時計回りに正に回転します。この 別の 符号規則では、正のせん断応力は反時計回りのせん断応力でもあり、両方とも下向きにプロットされるため、図 5 の符号規則 #2 と同じ円が生成されます 。また、負のせん断応力は 時計回りのせん断応力であり、両方とも上向きにプロットされます。
2
θ
{\displaystyle 2\theta }
τ
n
{\displaystyle \tau _{\mathrm {n} }}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
この記事は、物理空間については工学力学の符号規則に従い、 モール円空間については 代替の符号規則に従います(図5の符号規則#3)。
モール円の描画
図 6. 平面応力および平面ひずみ条件のモール円 (2 倍角アプローチ) 。応力解析後、 物質点における応力成分 、 、 が 判明します。これらの応力成分は 、 を通る2 つの垂直な平面と に作用します。 モール円上の点 と の座標は、それぞれ物質要素の平面 とに作用する応力成分です。次にモール円を使用して、 を通る 他の平面に作用する円上の応力成分 および 、 つまり任意の応力点の座標 を見つけます 。線と の間の角度は 、 を通る平面 と の法線ベクトルの間の 角度の 2 倍です 。
σ
x
{\displaystyle \sigma _{x}}
σ
y
{\displaystyle \sigma _{y}}
τ
x
y
{\displaystyle \tau _{xy}}
P
{\displaystyle P}
A
{\displaystyle A}
B
{\displaystyle B}
P
{\displaystyle P}
A
{\displaystyle A}
B
{\displaystyle B}
A
{\displaystyle A}
B
{\displaystyle B}
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
D
{\displaystyle D}
D
{\displaystyle D}
P
{\displaystyle P}
O
B
¯
{\displaystyle {\overline {OB}}}
O
D
¯
{\displaystyle {\overline {OD}}}
θ
{\displaystyle \theta }
B
{\displaystyle B}
D
{\displaystyle D}
P
{\displaystyle P}
図 4 に示すように、研究対象の物体内の点における 応力成分 、 、 がわかっていると仮定する と、 における応力の状態に対するモール円を作成する手順は次のようになります 。
σ
x
{\displaystyle \sigma _{x}}
σ
y
{\displaystyle \sigma _{y}}
τ
x
y
{\displaystyle \tau _{xy}}
P
{\displaystyle P}
P
{\displaystyle P}
水平 軸と垂直 軸を持つ 直交座標系を描きます 。
(
σ
n
,
τ
n
)
{\displaystyle (\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
選択した符号規則に従って、両方の垂直な平面 と 上の既知の応力成分に対応する空間 に、 それぞれ2つの点 とをプロットします(図4と6)。
A
(
σ
y
,
τ
x
y
)
{\displaystyle A(\sigma _{y},\tau _{xy})}
B
(
σ
x
,
−
τ
x
y
)
{\displaystyle B(\sigma _{x},-\tau _{xy})}
(
σ
n
,
τ
n
)
{\displaystyle (\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
A
{\displaystyle A}
B
{\displaystyle B}
点と点 を 直線で結んで 円の直径を描きます 。
A
{\displaystyle A}
B
{\displaystyle B}
A
B
¯
{\displaystyle {\overline {AB}}}
モール円を描きます 。円の 中心は直径線の中点であり 、この線と軸の交点に相当します 。
O
{\displaystyle O}
A
B
¯
{\displaystyle {\overline {AB}}}
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
主法線応力の検出
2D 回転要素の応力成分 。 クリックするとアニメーションが表示されます 。長方形要素の方向の角度が変化すると、その面 (エッジ) の応力成分がどのように変化するかを示す例。 主応力は 、せん断応力がすべての面から同時に消えたときに発生します。これが発生する方向によって、 主方向 が決まります。この例では、長方形が水平の場合、応力は次のように表されます 。対応するモール円表現は下部に示されています。
[
σ
x
x
τ
x
y
τ
y
x
σ
y
y
]
=
[
−
10
10
10
15
]
.
{\displaystyle \left[{\begin{matrix}\sigma _{xx}&\tau _{xy}\\\tau _{yx}&\sigma _{yy}\end{matrix}}\right]=\left[{\begin{matrix}-10&10\\10&15\end{matrix}}\right].}
主応力 の大きさは、 円が - 軸と交差する 点 と(図 6) の横座標です 。主主応力の大きさは 、常にこれら 2 点の横座標の絶対値のうち最大のものになります。同様に、副主応力の大きさは、 常にこれら 2 点の横座標の絶対値のうち最小のものになります。予想どおり、これら 2 点の縦座標は 0 で、主平面のせん断応力成分の大きさに対応します。あるいは、主応力の値は次のように求めることもできます。
C
{\displaystyle C}
E
{\displaystyle E}
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
σ
1
{\displaystyle \sigma _{1}}
σ
2
{\displaystyle \sigma _{2}}
σ
1
=
σ
max
=
σ
avg
+
R
{\displaystyle \sigma _{1}=\sigma _{\max }=\sigma _{\text{avg}}+R}
σ
2
=
σ
min
=
σ
avg
−
R
{\displaystyle \sigma _{2}=\sigma _{\min }=\sigma _{\text{avg}}-R}
ここで、平均法線応力 の大きさは 中心の横座標であり 、次のように表される。
σ
avg
{\displaystyle \sigma _{\text{avg}}}
O
{\displaystyle O}
σ
avg
=
1
2
(
σ
x
+
σ
y
)
{\displaystyle \sigma _{\text{avg}}={\tfrac {1}{2}}(\sigma _{x}+\sigma _{y})}
円の半径 の長さ (2点を通る円の方程式に基づく)は次のように表される。
R
{\displaystyle R}
R
=
[
1
2
(
σ
x
−
σ
y
)
]
2
+
τ
x
y
2
{\displaystyle R={\sqrt {\left[{\tfrac {1}{2}}(\sigma _{x}-\sigma _{y})\right]^{2}+\tau _{xy}^{2}}}}
最大および最小せん断応力の検出
最大せん断応力と最小せん断応力は、それぞれ円の最高点と最低点の座標に対応します。これらの点は、円の中心を通る垂直線と円の交点にあります。 したがって、最大せん断応力と最小せん断応力の大きさは、円の半径の値に等しくなります。
O
{\displaystyle O}
R
{\displaystyle R}
τ
max
,
min
=
±
R
{\displaystyle \tau _{\max ,\min }=\pm R}
任意の平面上の応力成分の検出
前述のように、2 次元応力解析を実行すると、 物質点 における応力成分 、 、がわかります 。これらの応力成分は、 図 5 と 6 に示すように、2 つの垂直な平面 と に作用します。モール円は、 平面 と角度を 成し て を通過する 他の平面に作用する応力成分 と 、つまり 円上の任意の点の座標を見つけるために 使用されます。これには、2 倍角と、平面の極または原点という 2 つのアプローチを使用できます。
σ
x
{\displaystyle \sigma _{x}}
σ
y
{\displaystyle \sigma _{y}}
τ
x
y
{\displaystyle \tau _{xy}}
P
{\displaystyle P}
A
{\displaystyle A}
B
{\displaystyle B}
P
{\displaystyle P}
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
D
{\displaystyle D}
D
{\displaystyle D}
P
{\displaystyle P}
θ
{\displaystyle \theta }
B
{\displaystyle B}
ダブルアングル
図 6 に示すように、応力が作用する平面に対して反時計回り の 角度で 平面に作用する 応力成分を決定するには、既知の応力点 から点まで 円の周りを同じ反時計回り方向に 角度 移動します。 つまり、 モール円内の
線 と線の間の角度 です。
(
σ
n
,
τ
n
)
{\displaystyle (\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
D
{\displaystyle D}
θ
{\displaystyle \theta }
B
{\displaystyle B}
σ
x
{\displaystyle \sigma _{x}}
2
θ
{\displaystyle 2\theta }
B
(
σ
x
,
−
τ
x
y
)
{\displaystyle B(\sigma _{x},-\tau _{xy})}
D
(
σ
n
,
τ
n
)
{\displaystyle D(\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
2
θ
{\displaystyle 2\theta }
O
B
¯
{\displaystyle {\overline {OB}}}
O
D
¯
{\displaystyle {\overline {OD}}}
2 倍角アプローチは、 通過する任意の 2 つの物理平面に対する法線ベクトル間の角度 (図 4) が、モール円上の対応する応力点 と円の中心を結ぶ 2 本の線間の角度の半分であるという事実に基づいています。
θ
{\displaystyle \theta }
P
{\displaystyle P}
(
σ
n
,
τ
n
)
{\displaystyle (\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
この 2 倍の角度の関係は、モール円の媒介変数方程式が の関数であるという事実から生じます。また、図 5 の 周りの材料要素の平面 と平面 は角度 で分離されており 、モール円では 角度 (角度の 2 倍)
で表されている こともわかります。
2
θ
{\displaystyle 2\theta }
A
{\displaystyle A}
B
{\displaystyle B}
P
{\displaystyle P}
θ
=
90
∘
{\displaystyle \theta =90^{\circ }}
180
∘
{\displaystyle 180^{\circ }}
平面の極または原点
図 7. 平面応力および平面ひずみ条件のモール円 (極アプローチ)。極から引いた直線は、その直線と同じ方向 (平行) に傾斜した平面上の応力状態を表す点でモール円と交差します。
2 番目のアプローチでは、モール円上の点、つまり 極 または 平面の原点 を決定します。極から引いた直線はどれも、その直線と同じ方向 (平行) に傾斜した平面上の応力状態を表す点でモール円と交差します。したがって、任意の特定の平面上の応力成分とがわかっていれば 、 特定の座標 と モール円を通るその平面に平行な線を引き、その線とモール円の交点として極を見つけることができます。例として、図 7 に示すように、応力成分 、 、 を持つ応力状態があるとします。まず、点 から の 作用面に平行な 線を引くか、または別の方法を選択した場合は、点 から の 作用面に平行な 線を引くことができます 。これら 2 本の線のいずれかとモール円の交点が極です。極が決定したら、 垂直と角度をなす平面、つまり法線ベクトルが 水平面と角度をなす平面上の応力状態を調べるために、極からその平面に平行な線を引くことができます (図 7 を参照)。その平面上の法線応力とせん断応力は、線とモール円の交点の座標になります。
σ
{\displaystyle \sigma }
τ
{\displaystyle \tau }
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
σ
x
,
{\displaystyle \sigma _{x},\!}
σ
y
,
{\displaystyle \sigma _{y},\!}
τ
x
y
,
{\displaystyle \tau _{xy},\!}
B
{\displaystyle B}
σ
x
{\displaystyle \sigma _{x}}
A
{\displaystyle A}
σ
y
{\displaystyle \sigma _{y}}
θ
{\displaystyle \theta }
θ
{\displaystyle \theta }
主平面の方向を見つける
最大主応力と最小主応力が作用する平面(主平面 とも呼ばれる)の方向は 、モール円でそれぞれ角度 ∠BOC と ∠BOE を測定し、それらの各角度の半分を取ることによって決定できます。したがって、と 間の角度 ∠BOC は 、 主主平面が平面となす 角度の 2 倍です 。
O
B
¯
{\displaystyle {\overline {OB}}}
O
C
¯
{\displaystyle {\overline {OC}}}
θ
p
{\displaystyle \theta _{p}}
B
{\displaystyle B}
角度 とは 次の式からも求めることができます。
θ
p
1
{\displaystyle \theta _{p1}}
θ
p
2
{\displaystyle \theta _{p2}}
tan
2
θ
p
=
2
τ
x
y
σ
y
−
σ
x
{\displaystyle \tan 2\theta _{\mathrm {p} }={\frac {2\tau _{xy}}{\sigma _{y}-\sigma _{x}}}}
この方程式は、離れている の 2 つの値を定義します (図)。この方程式は、円の幾何学から直接導き出すことも、円の のパラメトリック方程式を 0 に等しくすることによって導き出すこともできます (主平面のせん断応力は常に 0 です)。
θ
p
{\displaystyle \theta _{\mathrm {p} }}
90
∘
{\displaystyle 90^{\circ }}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
例
図8
図9
図 8 および図 9 に示すように、応力状態にある材料要素の 1 つの側面が水平面に対して 10° の方向を向いていると仮定します。モール円を使用して、次の値を求めます。
行動面の方向。
最大せん断応力とそれらの作用面の方向。
水平面上の応力成分。
応力変換式または応力変換法則を使用して答えを確認します。
解決策:
物理空間の工学力学の符号規則 (図 5) に従うと、この例の材料要素の応力成分は次のようになります。
σ
x
′
=
−
10
MPa
{\displaystyle \sigma _{x'}=-10{\textrm {MPa}}}
σ
y
′
=
50
MPa
{\displaystyle \sigma _{y'}=50{\textrm {MPa}}}
τ
x
′
y
′
=
40
MPa
{\displaystyle \tau _{x'y'}=40{\textrm {MPa}}}
。
この特定の応力状態のモール円を描く手順に従って、まず、 -軸を上に 向けた直交座標系を描きます 。
(
σ
n
,
τ
n
)
{\displaystyle (\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
τ
n
{\displaystyle \tau _{\mathrm {n} }}
次に、図 8 と図 9 の両方に示すように、平面 A と平面 B での応力の状態を表す 2 つの点 A(50,40) と B(-10,-40) をプロットします。これらの点は、モール円空間 (図 5) の工学力学の符号規則に従っています。この規則では、材料要素から外側に向かう正の法線応力と、材料要素を時計回りに回転させる各平面での正のせん断応力を想定しています。このように、平面 B に作用するせん断応力は負で、平面 A に作用するせん断応力は正です。円の直径は、点 A と点 B を結ぶ線です。円の中心は、この線と - 軸の交点です。中心の位置と直径の長さがわかっているので、この特定の応力状態のモール円をプロットできます。
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
軸と交差する点 E と点 C (図 8 と図 9) の横座標は 、それぞれ最小および最大の法線応力の大きさです。点 E と点 C の縦座標は、それぞれ副主平面と主主平面の両方に作用するせん断応力の大きさで、主平面の場合はゼロです。
σ
n
{\displaystyle \sigma _{\mathrm {n} }}
モール円を使用する目的は、円上の異なる点の座標を実際に測定して、異なる応力成分をグラフィカルに見つけることですが、結果を解析的に確認する方が便利です。したがって、円の中心の半径と横座標は、
R
=
[
1
2
(
σ
x
−
σ
y
)
]
2
+
τ
x
y
2
=
[
1
2
(
−
10
−
50
)
]
2
+
40
2
=
50
MPa
{\displaystyle {\begin{aligned}R&={\sqrt {\left[{\tfrac {1}{2}}(\sigma _{x}-\sigma _{y})\right]^{2}+\tau _{xy}^{2}}}\\&={\sqrt {\left[{\tfrac {1}{2}}(-10-50)\right]^{2}+40^{2}}}\\&=50{\textrm {MPa}}\\\end{aligned}}}
σ
a
v
g
=
1
2
(
σ
x
+
σ
y
)
=
1
2
(
−
10
+
50
)
=
20
MPa
{\displaystyle {\begin{aligned}\sigma _{\mathrm {avg} }&={\tfrac {1}{2}}(\sigma _{x}+\sigma _{y})\\&={\tfrac {1}{2}}(-10+50)\\&=20{\textrm {MPa}}\\\end{aligned}}}
そして主応力は
σ
1
=
σ
a
v
g
+
R
=
70
MPa
{\displaystyle {\begin{aligned}\sigma _{1}&=\sigma _{\mathrm {avg} }+R\\&=70{\textrm {MPa}}\\\end{aligned}}}
σ
2
=
σ
a
v
g
−
R
=
−
30
MPa
{\displaystyle {\begin{aligned}\sigma _{2}&=\sigma _{\mathrm {avg} }-R\\&=-30{\textrm {MPa}}\\\end{aligned}}}
点Hと点G(図8と図9)の座標は、それぞれ最小せん断応力と最大せん断応力の大きさです。点Hと点Gの横軸は、それぞれ最小せん断応力と最大せん断応力が作用する同じ平面に作用する法線応力の大きさです。最小せん断応力と最大せん断応力の大きさは、次のように解析的に求めることができます。
τ
max
,
min
=
±
R
=
±
50
MPa
{\displaystyle \tau _{\max ,\min }=\pm R=\pm 50{\textrm {MPa}}}
最小せん断応力と最大せん断応力が作用する同じ平面に作用する法線応力は、
σ
a
v
g
{\displaystyle \sigma _{\mathrm {avg} }}
主法線応力と主せん断応力の方向を見つけるには、二重角度アプローチ (図 8) または極アプローチ (図 9) のいずれかを選択できます。
2倍角アプローチを使用して、モール円(図8)の∠BOCと∠BOEの角度を測定し、物理空間で主応力と副主応力が平面Bとなす角度の2倍を見つけます。これらの角度のより正確な値を取得するには、手動で角度を測定する代わりに、解析式を使用できます。
2
θ
p
=
arctan
2
τ
x
y
σ
x
−
σ
y
=
arctan
2
∗
40
(
−
10
−
50
)
=
−
arctan
4
3
{\displaystyle {\begin{aligned}2\theta _{\mathrm {p} }=\arctan {\frac {2\tau _{xy}}{\sigma _{x}-\sigma _{y}}}=\arctan {\frac {2*40}{(-10-50)}}=-\arctan {\frac {4}{3}}\end{aligned}}}
一つの解は次のようになります 。図8を見ると、この値は∠BOEの角度に対応しています。したがって、小主角は
2
θ
p
=
−
53.13
∘
{\displaystyle 2\theta _{p}=-53.13^{\circ }}
θ
p
2
=
−
26.565
∘
{\displaystyle \theta _{p2}=-26.565^{\circ }}
そして、主角は
2
θ
p
1
=
180
−
53.13
∘
=
126.87
∘
θ
p
1
=
63.435
∘
{\displaystyle {\begin{aligned}2\theta _{p1}&=180-53.13^{\circ }=126.87^{\circ }\\\theta _{p1}&=63.435^{\circ }\\\end{aligned}}}
この特定の例では 、 と は (-軸に向いている) の作用面に対する角度であり、 (-軸に向いている )
の作用面に対する角度ではないことに注意してください。
θ
p
1
{\displaystyle \theta _{p1}}
θ
p
2
{\displaystyle \theta _{p2}}
σ
x
′
{\displaystyle \sigma _{x'}}
x
′
{\displaystyle x'}
σ
x
{\displaystyle \sigma _{x}}
x
{\displaystyle x}
極アプローチを使用して、まず極または平面の原点を特定します。そのためには、モール円上の点 A を通り、水平に対して 10° 傾いた線、つまり、 作用する平面 A に平行な線を引きます。極は、この線がモール円と交差する点です (図 9)。極の位置を確認するには、作用する平面 B に平行なモール円上の点 B を通り、線を引きます 。この線は、極でもモール円と交差します (図 9)。
σ
y
′
{\displaystyle \sigma _{y'}}
σ
x
′
{\displaystyle \sigma _{x'}}
極からモール円上のさまざまな点まで線を引きます。これらの線がモール円と交差する点の座標は、線と同じ傾きを持つ物理空間の平面に作用する応力成分を示します。たとえば、極から円内の点 C までの線は、が 作用する物理空間の平面と同じ傾きを持ちます。この平面は、モール円空間でも物理空間でも、平面 B と 63.435° の角度をなします。同様に、極から点 E、D、F、G、H まで線をたどり、同じ向きの平面上の応力成分を見つけます。
σ
1
{\displaystyle \sigma _{1}}
一般的な3次元応力状態を表すモール円
図10. 3次元応力状態のモール円
一般的な 3 次元の点における応力のモール円を構築するには、まず 主応力 の値とその 主方向を 評価する必要があります。
(
σ
1
,
σ
2
,
σ
3
)
{\displaystyle \left(\sigma _{1},\sigma _{2},\sigma _{3}\right)}
(
n
1
,
n
2
,
n
3
)
{\displaystyle \left(n_{1},n_{2},n_{3}\right)}
一般的な座標系、座標系 の 代わりに主軸を座標系として考え 、 と仮定すると 、 単位ベクトル を持つ与えられた平面に対する 応力ベクトル の 法線成分とせん断成分は 、次の式を満たします。
x
1
{\displaystyle x_{1}}
x
2
{\displaystyle x_{2}}
x
3
{\displaystyle x_{3}}
σ
1
>
σ
2
>
σ
3
{\displaystyle \sigma _{1}>\sigma _{2}>\sigma _{3}}
T
(
n
)
{\displaystyle \mathbf {T} ^{(\mathbf {n} )}}
n
{\displaystyle \mathbf {n} }
(
T
(
n
)
)
2
=
σ
i
j
σ
i
k
n
j
n
k
σ
n
2
+
τ
n
2
=
σ
1
2
n
1
2
+
σ
2
2
n
2
2
+
σ
3
2
n
3
2
{\displaystyle {\begin{aligned}\left(T^{(n)}\right)^{2}&=\sigma _{ij}\sigma _{ik}n_{j}n_{k}\\\sigma _{\mathrm {n} }^{2}+\tau _{\mathrm {n} }^{2}&=\sigma _{1}^{2}n_{1}^{2}+\sigma _{2}^{2}n_{2}^{2}+\sigma _{3}^{2}n_{3}^{2}\end{aligned}}}
σ
n
=
σ
1
n
1
2
+
σ
2
n
2
2
+
σ
3
n
3
2
.
{\displaystyle \sigma _{\mathrm {n} }=\sigma _{1}n_{1}^{2}+\sigma _{2}n_{2}^{2}+\sigma _{3}n_{3}^{2}.}
がわかっているので、 ガウスの消去法 を使って 、、、 を 解くことができ 、次の式が得られます。
n
i
n
i
=
n
1
2
+
n
2
2
+
n
3
2
=
1
{\displaystyle n_{i}n_{i}=n_{1}^{2}+n_{2}^{2}+n_{3}^{2}=1}
n
1
2
{\displaystyle n_{1}^{2}}
n
2
2
{\displaystyle n_{2}^{2}}
n
3
2
{\displaystyle n_{3}^{2}}
n
1
2
=
τ
n
2
+
(
σ
n
−
σ
2
)
(
σ
n
−
σ
3
)
(
σ
1
−
σ
2
)
(
σ
1
−
σ
3
)
≥
0
n
2
2
=
τ
n
2
+
(
σ
n
−
σ
3
)
(
σ
n
−
σ
1
)
(
σ
2
−
σ
3
)
(
σ
2
−
σ
1
)
≥
0
n
3
2
=
τ
n
2
+
(
σ
n
−
σ
1
)
(
σ
n
−
σ
2
)
(
σ
3
−
σ
1
)
(
σ
3
−
σ
2
)
≥
0.
{\displaystyle {\begin{aligned}n_{1}^{2}&={\frac {\tau _{\mathrm {n} }^{2}+(\sigma _{\mathrm {n} }-\sigma _{2})(\sigma _{\mathrm {n} }-\sigma _{3})}{(\sigma _{1}-\sigma _{2})(\sigma _{1}-\sigma _{3})}}\geq 0\\n_{2}^{2}&={\frac {\tau _{\mathrm {n} }^{2}+(\sigma _{\mathrm {n} }-\sigma _{3})(\sigma _{\mathrm {n} }-\sigma _{1})}{(\sigma _{2}-\sigma _{3})(\sigma _{2}-\sigma _{1})}}\geq 0\\n_{3}^{2}&={\frac {\tau _{\mathrm {n} }^{2}+(\sigma _{\mathrm {n} }-\sigma _{1})(\sigma _{\mathrm {n} }-\sigma _{2})}{(\sigma _{3}-\sigma _{1})(\sigma _{3}-\sigma _{2})}}\geq 0.\end{aligned}}}
、およびは 非負な
ので、これらの式の分子は次式を満たす。
σ
1
>
σ
2
>
σ
3
{\displaystyle \sigma _{1}>\sigma _{2}>\sigma _{3}}
(
n
i
)
2
{\displaystyle (n_{i})^{2}}
τ
n
2
+
(
σ
n
−
σ
2
)
(
σ
n
−
σ
3
)
≥
0
{\displaystyle \tau _{\mathrm {n} }^{2}+(\sigma _{\mathrm {n} }-\sigma _{2})(\sigma _{\mathrm {n} }-\sigma _{3})\geq 0}
分母 として
σ
1
−
σ
2
>
0
{\displaystyle \sigma _{1}-\sigma _{2}>0}
σ
1
−
σ
3
>
0
{\displaystyle \sigma _{1}-\sigma _{3}>0}
τ
n
2
+
(
σ
n
−
σ
3
)
(
σ
n
−
σ
1
)
≤
0
{\displaystyle \tau _{\mathrm {n} }^{2}+(\sigma _{\mathrm {n} }-\sigma _{3})(\sigma _{\mathrm {n} }-\sigma _{1})\leq 0}
分母 として
σ
2
−
σ
3
>
0
{\displaystyle \sigma _{2}-\sigma _{3}>0}
σ
2
−
σ
1
<
0
{\displaystyle \sigma _{2}-\sigma _{1}<0}
τ
n
2
+
(
σ
n
−
σ
1
)
(
σ
n
−
σ
2
)
≥
0
{\displaystyle \tau _{\mathrm {n} }^{2}+(\sigma _{\mathrm {n} }-\sigma _{1})(\sigma _{\mathrm {n} }-\sigma _{2})\geq 0}
分母 として
σ
3
−
σ
1
<
0
{\displaystyle \sigma _{3}-\sigma _{1}<0}
σ
3
−
σ
2
<
0.
{\displaystyle \sigma _{3}-\sigma _{2}<0.}
これらの式は次のように書き直すことができる。
τ
n
2
+
[
σ
n
−
1
2
(
σ
2
+
σ
3
)
]
2
≥
(
1
2
(
σ
2
−
σ
3
)
)
2
τ
n
2
+
[
σ
n
−
1
2
(
σ
1
+
σ
3
)
]
2
≤
(
1
2
(
σ
1
−
σ
3
)
)
2
τ
n
2
+
[
σ
n
−
1
2
(
σ
1
+
σ
2
)
]
2
≥
(
1
2
(
σ
1
−
σ
2
)
)
2
{\displaystyle {\begin{aligned}\tau _{\mathrm {n} }^{2}+\left[\sigma _{\mathrm {n} }-{\tfrac {1}{2}}(\sigma _{2}+\sigma _{3})\right]^{2}\geq \left({\tfrac {1}{2}}(\sigma _{2}-\sigma _{3})\right)^{2}\\\tau _{\mathrm {n} }^{2}+\left[\sigma _{\mathrm {n} }-{\tfrac {1}{2}}(\sigma _{1}+\sigma _{3})\right]^{2}\leq \left({\tfrac {1}{2}}(\sigma _{1}-\sigma _{3})\right)^{2}\\\tau _{\mathrm {n} }^{2}+\left[\sigma _{\mathrm {n} }-{\tfrac {1}{2}}(\sigma _{1}+\sigma _{2})\right]^{2}\geq \left({\tfrac {1}{2}}(\sigma _{1}-\sigma _{2})\right)^{2}\\\end{aligned}}}
これらは、応力 、 、 に対する 3 つのモール円の方程式で、半径 は それぞれ 、 、 、 、中心の座標はそれぞれ 、 、 です。
C
1
{\displaystyle C_{1}}
C
2
{\displaystyle C_{2}}
C
3
{\displaystyle C_{3}}
R
1
=
1
2
(
σ
2
−
σ
3
)
{\displaystyle R_{1}={\tfrac {1}{2}}(\sigma _{2}-\sigma _{3})}
R
2
=
1
2
(
σ
1
−
σ
3
)
{\displaystyle R_{2}={\tfrac {1}{2}}(\sigma _{1}-\sigma _{3})}
R
3
=
1
2
(
σ
1
−
σ
2
)
{\displaystyle R_{3}={\tfrac {1}{2}}(\sigma _{1}-\sigma _{2})}
[
1
2
(
σ
2
+
σ
3
)
,
0
]
{\displaystyle \left[{\tfrac {1}{2}}(\sigma _{2}+\sigma _{3}),0\right]}
[
1
2
(
σ
1
+
σ
3
)
,
0
]
{\displaystyle \left[{\tfrac {1}{2}}(\sigma _{1}+\sigma _{3}),0\right]}
[
1
2
(
σ
1
+
σ
2
)
,
0
]
{\displaystyle \left[{\tfrac {1}{2}}(\sigma _{1}+\sigma _{2}),0\right]}
モール円に関するこれらの方程式は、すべての許容応力点が これらの円上または円で囲まれた陰影付きの領域内にあることを示しています (図 10 を参照)。 円 に関する方程式を満たす応力点は、円上または円の外側にあります。 円に関する方程式を満たす 応力点は 、円上または円の内側にあります 。そして最後に、 円に関する方程式を満たす応力点は 、円上または円の外側にあります 。
(
σ
n
,
τ
n
)
{\displaystyle (\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
(
σ
n
,
τ
n
)
{\displaystyle (\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
C
1
{\displaystyle C_{1}}
C
1
{\displaystyle C_{1}}
(
σ
n
,
τ
n
)
{\displaystyle (\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
C
2
{\displaystyle C_{2}}
C
2
{\displaystyle C_{2}}
(
σ
n
,
τ
n
)
{\displaystyle (\sigma _{\mathrm {n} },\tau _{\mathrm {n} })}
C
3
{\displaystyle C_{3}}
C
3
{\displaystyle C_{3}}
参照
参考文献
^ 「主応力と主平面」 www.engineeringapps.net . 2019年12月25日 閲覧 。
^ パリー、リチャード・ホーリー・グレイ(2004年)。モール円、応力経路、地盤工学(第2版)。テイラー&フランシス。pp. 1–30。ISBN 0-415-27297-1 。
^ Gere, James M. (2013). 材料力学 . Goodno, Barry J. (第 8 版). スタンフォード、コネチカット州: Cengage Learning. ISBN 9781111577735 。
文献
ビール、フェルディナンド・ピエール、エルウッド・ラッセル・ジョンストン、ジョン・T・デウルフ(1992)。 『材料力学』 。McGraw -Hill Professional。ISBN 0-07-112939-1 。
Brady, BHG; ET Brown (1993)。地下採掘のための岩石力学(第3版)。Kluwer Academic Publisher。pp. 17–29。ISBN 0-412-47550-2 。
デイビス、RO; セルバドゥライ。APS (1996)。弾性と地盤力学。ケンブリッジ大学出版局。pp. 16–26。ISBN 0-521-49827-9 。
ホルツ、ロバート D.、コバックス、ウィリアム D. (1981)。土質工学入門。プレンティス ホール土木工学および工学力学シリーズ。プレンティス ホール 。ISBN 0-13-484394-0 。
イェーガー、ジョン・コンラッド、クック、NGW、ジマーマン、RW (2007)。岩石力学の基礎(第4版)。ワイリー・ブラックウェル 。pp.9–41。ISBN 978-0-632-05759-7 。
ジュミキス、アルフレッド・R. (1969)。理論土質力学:土質力学と基礎工学への実際的応用。Van Nostrand Reinhold Co. ISBN 0-442-04199-3 。
パリー、リチャード・ホーリー・グレイ(2004年)。モール円、応力経路、地盤工学(第2版)。テイラー&フランシス。pp. 1–30。ISBN 0-415-27297-1 。
ティモシェンコ、スティーブン P. 、 ジェームズ・ノーマン・グーディアー (1970)。 弾性理論 (第 3 版)。マグロウヒル インターナショナル エディション 。ISBN 0-07-085805-5 。
ティモシェンコ、スティーブン P. (1983)。 材料強度の歴史:弾性理論と構造理論の歴史の簡単な説明付き 。ドーバー物理学の本。ドーバー出版 。ISBN 0-486-61187-6 。
外部リンク
ウィキバーシティにはモール円 に関する学習リソースがあります
ウィキメディア・コモンズには、モール円 に関連するメディアがあります 。
レベッカ・ブランノンによるモールの円とその他の円
DoITPoMS 教育学習パッケージ - 「応力解析とモール円」