プログラミングの記事を書きます。-korokoro前回の記事では、コーシー応力\(\sigma_{ij}\)の物理的な意味と、その式を示しました。https://takun-physics.net/?p=4966https://takun-physics.net/?p=4992このコーシー応力\(\sigma_{ij}\)は、流体力学の運動量保存則の圧力項と粘性項に関わる内力で重要です。K_{i}\end{align*}本記事を読むと以下のことがわかるようになります。(1)\end{align*}(1)式はこのように書いたりもします。(2)\end{align*}※添え字が2回続くと和を意味しています(アインシュタインの縮約記法)。やって(1)式を導くのかですが、簡単です。「任意の面\(S\)にはたらく力」と「三角錐の各面(面\(S\)以外)にはたらく力」とのつり合いの式を立てれば良いのです。こちらの記事もどうぞ任意の面\(S\)を用意する(準備段階)点\(P\)を頂点とする微小な三角錐\(PABC\)を考えます。緑の面\(S\)にはたらく応力ベクトル\(\boldsymbol{t}\)とコーシー応力\(\boldsymbol{\sigma}\)との関係を知ることが本記事の目標となります。※緑の面の面積は\(S\)法線ベクトルと応力ベクトルでは、任意の面上の法線ベクトル\(\boldsymbol{n}\)と応力ベクトル\(\boldsymbol{t}\)を下記のように置きます。ここは、文字を定義しただけですね。法線ベクトルは別の書き方では、こんな感じで、「各軸におろした角度」を使って方向余弦を定義できます。続いては、三角錐の各面の面積を定義していきます。\(x\)軸に垂直な面の面積:\(S_{x}\)\(y\)軸に垂直な面の面積:\(S_{y}\)\(z\)軸に垂直な面の面積:\(S_{z}\)このように「~軸に垂直な面の面積」を添え字を使って区別します。応力は「単位面積あたりの力」を意味していますので、各面にはたらく力を考えるときには、各面の面積を定義しておく必要があるので定義しました。各面にはたらく応力ベクトル続いて各面にはたらく応力ベクトルを書いていきます。応力ベクトルは2種類あります。添え字が多くて少々煩雑に見えますが、添え字の意味を理解しておくとたいしたことはないです。垂直応力は添え字がひとつしかないのでわかりやすいですが、せん断応力は添え字が2つあるのでどういった意味なのかを示しておきます。これで準備段階は終わりました。任意の面にはたらく応力ベクトルと、三角錐の各面にはたらく力とのつり合いの式を考えることになります。\(x\)軸方向の力のつり合いどの軸も特別な軸というわけではないので、どれかひとつだけの軸について、どのような式が導けるのかを考えれば十分ですよね。それを\(x\)軸にとることにります。\(x\)軸方向の力のつり合いを考えていきます。任意の面(緑色)にはたらく応力ベクトルの\(x\)方向の応力を\(t_{x}\)と置きます。(4)\end{align*}と、こうのようになります。これでは\(S\)や\(S_{x}\)などの面積の表記が邪魔ですよね。そこで、\(S\)や\(S_{x}\)などの関係式を示しておく必要があります。表題の通り、\(S_{x}\)は\(S\)を\(x\)軸方向に投影した面積です。(5)\end{align*}の関係式が成り立ちます。わからないっていう人は、面で考えるからわかりにくいのだと思います。法線ベクトル\(\boldsymbol{n}\)は面\(S\)に対して垂直なベクトルなのだから、面が線に見える方向から眺めてみれば良いです。\(zx\)面方向から面\(S\)を眺めてみましょう。すると、こんな感になります。(5)\end{align*}こうなることがわかりますよね。(8)\end{align*}となります。(10)\end{align*}このようになります。行列表示にするさて、これを行列での表記にするととてもすっきりします。(10)\end{align*}下記が(8)(9)(10)の行列での表記です。(11)\end{align*}ここで、\(x\)、\(y\)、\(z\)などの文字が添え字だと扱いにくい場合も多いので、添え字を数字に変えておきます。こんな感じで・・・そうすると、冒頭で示した応力ベクトルとコーシー応力の関係式が導けたことになります。(1)\end{align*}これでも十分よくまとまっていますが、もっとすっきりさせるためにテンソル表記にすることもあります。(1)式はこのように書いたりもします。(2)\end{align*}※添え字が2回続くと和を意味しています(アインシュタインの縮約記法)。順を追って行けばそれほど難しくもなかったかもしれません。(1)\end{align*}(1)式はこのように書いたりもします。(2)\end{align*}※添え字が2回続くと和を意味しています(アインシュタインの縮約記法)。(2)式は何に使うのか?よくありがちなのが、複雑で長い式を一生懸命計算していって導いた結果について、「ん?これはいったい何に使うのか?」ってなることです。流体力学の運動量保存則の圧力項と粘性項に関わる内力を表現するためです。※下記の\(\sigma_{ij}\)が何かを理解するためです。イメージがつかめてきたところで次回は、流体力学の基礎式の一般系から運動量保存則を導きたいと思います。-【理想気体の比熱の温度依存性(8)】量子論で2原子分子の比熱を計算。低温極限で調和振動子の量子化エネルギーは凍結している。【CAE解析の目的とは?】CAEは実機の代替えなのか?精度を無視してパラメータスタディすることは意味がないのか?【Udemyがお勧め理由5つ】プログラミングを勉強するなら動画学習おすすめ【まとめ】比熱の温度依存性について。古典論と量子論から理解する。22:32非常に図がわかりやすいです。00:44ありがとうございます(^^)パワポの2016バージョンを使っているので特に特別な方法は使っていません(笑)↓これを最近模索中です。15:11こんにちは。連続体力学の初学者で、いつも勉強に利用させていただいております。任意の面での応力についての解説文のなかで疑問が湧いたのですが、三角錐の各面にかかる力のつり合いを考えていますが、ここで力のつりあいが成立する物理的な要請ってあるのでしょうか?本質は力のつり合いではなく、応力がテンソル(線型写像)であることではないかと思いました。詰まる所、なぜ応力がテンソル表現なのか、という部分がうまく理解できません。単純に線型な物理モデルを扱っているだけのことなのか、それとも応力がテンソルというのは連続体力学で一般的な話なのでしょうか?。的はずれなコメントをしていたら申し訳ありません。12:49お読みいただきましてありがとうございます。少々回答が長くなり申し訳ございません。ここで力のつりあいが成立する物理的な要請ってあるのでしょうか?本質は力のつり合いではなく、応力がテンソル(線型写像)であることではないかと思いました。この2つはほぼ同じ意味であると解釈しています。ここで力のつりあいが成立する物理的な要請ってあるのでしょうか?つりあいが成り立つための物理的な要請があるのではなくて、つりあっていること自体を要請しています。これは、ニュートン力学第三法則の3つ目の作用反作用にあたります。三角錐が有限の大きさに見えますが、無限に小さくしていくと、「σx、τyx、τzx」で発生している内力は「任意の面tx」とつりあっているということを自体を要請しています。「Aさんが物体Bをある力で押すと、同じ力でAさんは物体Bから押し返されます」。これをAさんと物体Bを一体モノだと思うと今回考えている状況と同じです。内力なのでつりあっていること自体を要請しています。それでも納得がいかない場合もあるかと思います。その場合は、「最小作用の原理から仮想仕事の原理」というものを要請すると、自然と力のつりあいが導かれます。ここで考えている力は内力なので力がつりあっていると状態が最も実現可能な状態と考えるのです。本質は力のつり合いではなく、応力がテンソル(線型写像)であることではないかと思いました。これも正しいと思っています。作用反作用をFA=-FBと考えると、FBをマイナスに写像したものがFAだと考えることができるように、記事内であるように三角錐の各面をある面の応力へ移し変えただけという解釈にもなります。写像したものの式を眺めると力のつりあいの式であると後で解釈することも可能です。ぜ応力がテンソル表現なのか、という部分がうまく理解できません。単純に線型な物理モデルを扱っているだけのことなのか、それとも応力がテンソルというのは連続体力学で一般的な話なのでしょうか?力の作用面の法線の向きと力の作用方向とが一致してしていない応力成分(せん断応力)があるため、2階テンソルになっていると理解しております。2階テンソルといっても、「共変テンソル、反変テンソル、混合テンソル」などありますが、申し訳ないですがそのあたりは勉強不足なのでこれ以上は回答ができません。材料力学でも流体力学でも同じように連続体力学として応力を定義していますので、連続体力学では一般的ではないかと思っております。21:57korokoroさま返信ありがとうございます。私は社会人になってから物理の学習意欲が再燃したもので、久々に他人と物理の話ができて嬉しく思います。ところで応力がテンソル表現になる理由について、返信を読んでもすぐには腑に落ちなかったのですが、あわせて以下のリンクを読んで一応納得できました。σ(n)が成立することはすぐに納得できたのですが、ここからσが線型写像であることを結論するのは飛躍があるので、悶々としていました。リンク先の説明にあるように応力が面積力であることを加味するとたしかに応力σがテンソルでなければならないことが説明できるな、と一応わかった気にはなりました。わからないことがあったらコメントさせてください。コメントするコメントをキャンセルプロフィールカマキリ(^^)大学以上の物理を勉強してブログに書いていきます。大学の専攻は物性理論で、Fortranを使って数値計算をしていました。CAE解析(流体解析)・毎日note更新しています。内容は、「現在勉強中の内容」「日々思ったこと」を書いています。1DCAEの操作方法などを動画にしています。youtube・Qiitaではプログラミング言語の基本的な内容をまとめています。