連続体の運動方程式
連続体を支配する運動方程式
空間表示
$ \rho\frac{\mathrm{D}\bm{v}}{\mathrm{D}t}=\bm{\nabla}\cdot\bm{\sigma}+\bm k
$ \rho:空間表示した連続体の密度
$ \bm{v}:空間速度
$ \bm{\sigma}:表面力
$ \bm k:物体力
物質表示
その1: $ \tilde\rho J\ddot{\bm\phi}=\bm P\cdot\overleftarrow{\bm\nabla}+J\bm K
$ \tilde\rho:=\left.\rho\right|_{\bm x=\bm\phi(\bm X,t)}
$ \bm\phi:運動函数
$ J:=\det\bm F:体積変化率
$ \bm F:=\bm\phi\overleftarrow{\bm\nabla}:変形勾配tensor
$ \bm P:=J\left.\bm\sigma\right|_{\bm x=\bm\phi(\bm X,t)}\cdot{\bm F^{-1}}^\top:第1Piola-Kirchhoff応力tensor
$ \bm K:=\left.\bm k\right|_{\bm x=\bm\phi(\bm X,t)}
その2: $ \tilde\rho\ddot{\bm\phi}\cdot\bm F=\left.\bm\sigma\right|_{\bm x=\bm\phi(\bm X,t)}\cdot\overleftarrow{\bm\nabla}+\bm K\cdot\bm F
cf. 物質表示した連続体の運動方程式
圧縮性、非定常含めた全ての状態を表せる
流体の運動方程式をはじめとする、構造力学や流体力学における各種の運動方程式は、この方程式の外力項$ \bm{\nabla}\cdot\bm{\sigma}+\bm kを具体的に決定することで導出される
具体的に決定すること=構成則を決定すること
この式にデファクトスタンダードとなる名前は存在しなさそう
Cauchyの運動方程式と呼ばれることがある
/takkerでは、人名を含まず法則の意味を含む名前を主に使いたいので、「連続体の運動方程式」という表現を使うことにする
この方針はGaussの発散定理を発散定理と呼ぶ理由と同じ
派生式
積分形
$ \int_V\frac{\partial\rho\pmb{v}}{\partial t}\mathrm{d}v+\int_{\partial V}\rho\pmb{v}\pmb{v}\cdot\mathrm{d}\pmb{a}=\int_V\pmb{K}\mathrm{d}v+\int_{\partial V}\pmb{\sigma}\cdot\mathrm{d}\pmb{a}
左辺第2項は、領域を通過する運動量束の総和に相当する
$ \frac{\partial}{\partial t}\int_{B_t}\frac12\rho|\bm v|^2\mathrm dv+\int_{B_t}\bm\sigma:\bm d\mathrm dv=\int_{\partial B_t}\bm v\cdot\bm\sigma\cdot\mathrm d\bm a+\int_{B_t}\bm K\cdot\bm v\mathrm dv
連続体の力学的エネルギ保存則
微分系
連続体の力学的エネルギ保存則
$ \rho\frac{\mathrm D}{\mathrm Dt}\left(\frac12|\bm v|^2\right)+\bm\sigma:\bm d=\bm\nabla\cdot(\bm\sigma\cdot\bm v)+\bm K\cdot\bm v
エネルギー流束を導出する時に使う形
渦度方程式
$ \pmb{\nabla}\timesをとる
回転 (vector解析)の各種演算の$ \pmb{a}\times(\pmb{\nabla}\times\pmb{a})=\pmb{\nabla}\left(\frac12|\pmb{a}|^2\right)-(\pmb{a}\cdot\pmb{\nabla})\pmb{a}を使って変形すると
$ \frac{\partial\pmb v}{\partial t}+\pmb\nabla\left(\frac12|\pmb v|^2\right)=\pmb v\times\pmb\nabla\times\pmb v+\frac{\pmb K}\rho+\frac{\pmb\nabla\cdot\pmb\sigma}\rho
発散をとると
$ \implies\frac{\partial\pmb\nabla\cdot\pmb v}{\partial t}+\pmb\nabla^2\left(\frac12|\pmb v|^2\right)=\pmb\nabla\cdot(\pmb v\times\pmb\nabla\times\pmb v)+\pmb\nabla\frac1\rho\cdot(\pmb K+\pmb\nabla\cdot\pmb\sigma)+\frac1\rho(\pmb\nabla\cdot\pmb K+\pmb\nabla\cdot\pmb\nabla\cdot\pmb\sigma)
うーん、一般的に何かがわかるわけじゃあないか
流体の圧力が満たすべき方程式のようなきれいなものが出てくることを期待した
座標変換
流線座標での表示
#2024-04-28 22:01:53
#2024-04-05 18:44:59
#2023-09-23 09:19:48
#2022-11-05 05:47:10