流体力学コーシー応力の式【図でわかりやすく解説】

前回の記事で、流体力学の運動量保存則である、

運動量保存則

\begin{align*}\frac{\partial \rho v_{i}}{\partial t}+\frac{\partial\rho v_{i}v_{j}}{\partial x_{j}}=\frac{\partial \sigma_{ij}}{\partial x_{j}}+\rho K_{i}\cdot\cdot\cdot (1)\end{align*}

\(\sigma_{ij}\)（コーシー応力）が何かということを説明しました。

流体力学コーシー応力の意味は圧力と粘性抵抗【図でわかりやすく解説】流体力学の基礎方程式のひとつとして運動量保存則があります。それは、こんな感じ。運動量保存則 \begin{align*}\...

ただ上の記事で、コーシー応力は圧力と粘性抵抗であるという説明の記事を書きましたが、具体的な式を書いていませんでしたよね。

そこで・・・

本記事の内容
コーシー応力の具体な式を示したいと思います。
\begin{align*}\sigma_{ij}=\begin{pmatrix}
\sigma_{11} & \sigma_{12} & \sigma_{13}\\
\sigma_{21} & \sigma_{22} & \sigma_{23}\\
\sigma_{31} & \sigma_{32} & \sigma_{33}
\end{pmatrix}\cdot\cdot\cdot (2)\end{align*}

↑結果はこれです。

連続体近似をした時の内力をどのように式で表現するかということですので、コーシー応力は、別に流体力学だけに出てくる応力ではなくて、材料力学などでもよく出てきます。

\(x\)軸に垂直な面の内部にかかる力を考える
1. \(\boldsymbol{t}^{x}\)を\(x\)、\(y\)、\(z\)に分解する
\(x\)、\(y\)、\(z\)軸に垂直な面の内部にかかる力を考える
まとめ
しかし、立方体だけだと限定され過ぎていないか？

\(x\)軸に垂直な面の内部にかかる力を考える

内部にはたらく力を考えるわけですが、下記のように流れの領域のある1部分を切り取って2分割することを考えます。

まずは簡単のため、\(x\)軸に垂直な面で領域を2分割にしてみます。

そして、分割した片方（左側）の面にはいわゆる内力というのがはたらいています。

この内力は右側の面から受けている力から来ています。

では、その面の上の点\(P\)での\(x\)軸に垂直な単位面積あたりの力を考えましょう。

\begin{align*}\boldsymbol{t}^{x}=\frac{df^{x}}{dS^{x}}\cdot\cdot\cdot (3)\end{align*}

とこうなりますよね。

\(x\)方向の面：\(dS^{x}\)
\(x\)方向にはたらく力：\(df^{x}\)

※注意ですが、別に\(x\)軸方向の力を考えたわけではなくて、\(x\)軸に垂直な面に対する単位面積あたりの力を考えたわけです。

だから、\(x\)軸に垂直な面に対する力\(\boldsymbol{t}^{x}\)は、必ずしも面に対して垂直な方向の力とは限りません。

\(\boldsymbol{t}^{x}\)を\(x\)、\(y\)、\(z\)に分解する

左側だけの面位着目して\(x\)軸に垂直な面に対する力\(\boldsymbol{t}^{x}\)を\(x\)、\(y\)、\(z\)に分解することを考えます。

ポイント
面に垂直な方向の力と面に平行な方向の力に分解する
「どの軸に垂直な面か」と「どの方向の力か」を添え字で表現する
面に対して垂直な方向の応力：\(\sigma_{x}\)
面に対して平行な方向の応力：\(\tau_{xy}\)、\(\tau_{xz}\)

\(x\)、\(y\)、\(z\)軸に垂直な面の内部にかかる力を考える

先ほどまでは\(x\)軸に垂直な面にはたらく力しか書いていませんでしたが、\(x\)軸以外の軸についても同様の考え方ができるわけですよね。

わかりやすくするために点\(P\)のまわりを拡大して立方体を用意して見ましょう。

では、この立方体の表面に先ほど同様に応力を書いてみます。

はい、こんな感じです。

こらがコーシー応力の式です。

添え字が文字だと後で添え字が足らなくなる可能性があるので、よく添え字は数字に置き換えることが多いです。

内部にはたらく力は、こんな感じでテンソル表記になります。

テンソル表記は苦手な人が多いかもしれませんが、順をおって理解していけばだんだん慣れてくると思います。

まとめ

コーシー応力
\begin{align*}\sigma_{ij}=\begin{pmatrix}
\sigma_{11} & \sigma_{12} & \sigma_{13}\\
\sigma_{21} & \sigma_{22} & \sigma_{23}\\
\sigma_{31} & \sigma_{32} & \sigma_{33}
\end{pmatrix}\cdot\cdot\cdot (2)\end{align*}

添え字が2文字になると少々煩雑に見えますが、ポイントは

どの面に垂直な力を考えているか
どの軸の方向の力を考えているか

これを添え字で表現しているだけだと覚えていればOKです。

しかし、立方体だけだと限定され過ぎていないか？

本記事で書いた内部にはたらく力は、立方体で考えましたよね。

これって限定しすぎていないか？って思われる方もいるでしょう。

点Pを通る面なんて、無数に存在します。

だから任意の面に対する応力ベクトルとコーシー応力との関係を示しておく必要があります。

別の記事で書くとして、結果だけ載せておきます。

任意の面にはたらく応力ベクトル\(\boldsymbol{t}\)
\begin{align*}t_{i}=\begin{pmatrix}
t_{1} \\
t_{2} \\
t_{3}
\end{pmatrix}\end{align*}
↓
\begin{align*}t_{i}=\begin{pmatrix}
\sigma_{11} & \sigma_{12} & \sigma_{13}\\
\sigma_{21} & \sigma_{22} & \sigma_{23}\\
\sigma_{31} & \sigma_{32} & \sigma_{33}
\end{pmatrix}\begin{pmatrix}
n_{1} \\
n_{2} \\
n_{3}
\end{pmatrix}\cdot\cdot\cdot (4)\end{align*}

(4)式はこのように書いたりもします。
\begin{align*}t_{i}=\sigma_{ij}n_{i}\cdot\cdot\cdot (5)\end{align*}

※添え字が2回続くと和を意味しています（アインシュタインの縮約記法）。