Home > Mathematical Modelling, Theoretical Mechanics > Phương trình chuyển động Cauchy

Phương trình chuyển động Cauchy


Ngoại lực tác động lên một thể vật chất gồm có:
– Lực khối (body force) tác động vào mỗi đơn vị thể tích, ký hiệu bởi trường vector \mathbf{f}=({{f}_{i}}).
– Lực bề mặt (surface force) hay ứng suất (stress) là lực tác động mỗi đơn vị diện tích bề mặt.
Xét thể vật chất chiếm thể tích trong miền \Omega \subset {{\mathbb{R}}^{d}} với biên bị chặn, trơn từng miếng \partial \Omega với trường vector pháp tuyến tương ứng là \mathbf{n}. Giả sử {{\mathbf{T}}^{(\mathbf{n})}}=(T_{i}^{(\mathbf{n})}) là một trường vector ứng suất tác động lên mỗi đơn vị diện tích của bề mặt \partial \Omega . khi đó định lý ứng suất Cauchy (Cauchy’s stress theorem) phát biểu rằng tồn tại trường tensor bậc hai \mathbf{\sigma }=({{\sigma }_{ij}}) sao cho

\displaystyle{{\mathbf{T}}^{(\mathbf{n})}}=\mathbf{\sigma }\,.\,\mathbf{n} hay T_{j}^{(\mathbf{n})}={{\sigma }_{ij}}{{n}_{i}}.             (1)

https://i0.wp.com/upload.wikimedia.org/wikipedia/commons/thumb/b/b3/Components_stress_tensor_cartesian.svg/500px-Components_stress_tensor_cartesian.svg.png

Ta gọi \mathbf{\sigma }=({{\sigma }_{ij}}) là trường tensor ứng suất Cauchy (Cauchy stress tensor), và nó mô tả trường vector ứng suất như là một hàm phụ thuộc tuyến tính theo trường vector pháp tuyến.
Theo định luật thứ hai của Newton, ta có

\displaystyle \int_{\partial \Omega }{{{\mathbf{T}}^{(\mathbf{n})}}}ds+\int_{\Omega }{\mathbf{f}}dV=\frac{d}{dt}\int_{\Omega }{\rho }\mathbf{v}dV,      (2)

trong đó \mathbf{v} là trường vector vận tốc của mỗi phần tử của thể vật chất, \rho là mật độ vật chất.
Ta lại có

\displaystyle \frac{d}{dt}\int_{\Omega }{\rho \,\mathbf{v}dV} =\int_{\Omega }{\frac{\partial (\rho \,\mathbf{v})}{\partial t}dV}+\int_{\partial\Omega }{(\rho \,\mathbf{v}.\,\mathbf{n})\,\mathbf{v}\,ds}

\displaystyle =\int\limits_{\Omega }{\left[ \left( \rho \frac{\partial \mathbf{v}}{\partial t}+\frac{\partial \rho }{\partial t}\mathbf{v} \right)+\,\left( \nabla .\,\rho \,\mathbf{v}\, \right)\,\mathbf{v}+\left( \rho \mathbf{v}.\nabla \right)\mathbf{v} \right]dV}.

\displaystyle =\int\limits_{\Omega }{\rho \left[ \frac{\partial }{\partial t}+\mathbf{v}.\nabla \right]\mathbf{v}dV}+\int\limits_{\Omega }{\left[ \frac{\partial \rho }{\partial t}+\nabla .(\rho \,\mathbf{v}) \right]\mathbf{v}dV}.      (3)
trong đó chú ý áp dụng định lý Gauss–Ostrogradsky (Gauss–Ostrogradsky’s divergence theorem). Giả sử khối lượng được bảo tồn, khi đó số hạng thứ hai trong biểu thức cuối ứng với sự thay đổi của luồng khối lượng theo thời gian sẽ bị triệt tiêu (chúng ta chú ý phương trình của thể vật chất liên tục

\displaystyle \frac{\partial \rho }{\partial t}+\nabla .(\rho \,\mathbf{v}) =0

)

Thay thế (1)-(3) vào (2) và ta có

\displaystyle \int\limits_{\Omega }{\left( \nabla .\,\mathbf{\sigma }+\mathbf{f}-\rho \frac{D\,\mathbf{v}}{Dt} \right)}\,\,dV=\mathbf{0}.

Trong đó D/Dt là ký hiệu của đạo hàm vật chất (material derivative), xác định với mỗi trường vector \mathbf{\varphi } như sau

\displaystyle\frac{D\mathbf{\varphi }}{Dt}=\frac{\partial \mathbf{\varphi }}{\partial t}+(\mathbf{v}\,.\nabla )\mathbf{\varphi }.

Từ đó ta thu được

\displaystyle \nabla .\,\mathbf{\sigma }+\mathbf{f}=\rho \frac{D\mathbf{v}}{Dt}          (4)

và được gọi là phương trình chuyển động Cauchy (Cauchy’s equation of motion). Có thể viết (4) dưới dạng tọa độ Cartesian như sau

\displaystyle\sum\limits_{j=1}^{n}{\frac{\partial }{\partial {{x}_{j}}}{{\sigma }_{ij}}}+{{f}_{i}}=\rho \left[ \frac{\partial }{\partial t}+\sum\limits_{j=1}^{d}{{{v}_{j}}\frac{\partial }{\partial {{x}_{j}}}} \right]{{v}_{i}}

  1. November 19, 2011 at 11:42 pm

    Anh cho em quote của cái này đi.

  2. November 27, 2011 at 3:20 am

    Tuan Minh :
    okie. you’re welcome

    Đâu anh, quote đâu anh ?

  3. November 27, 2011 at 8:15 am

    ý em là gì?
    em có thể trích lại, nhưng nhớ ghi tên người viết🙂

  4. November 30, 2011 at 1:51 pm

    Vâng, em cảm ơn anh.
    Anh trích đi, bài viết của anh hay thiệt.😀

  1. No trackbacks yet.

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: