以应力表示的黏性流体运动微分方程的推导
1. 黏性流体的内应力
黏性流体在运动时,表面力不仅有法向应力,还有切向应力,因此黏性流体的表面力不垂直于作用面。
如在任一点取一微小的正六面体,如图所示,作用在平面ABCD上的力
有法向应力
pxx
,与切向应力xy
和xz。应力符号的第一个字母表示作
用面的外法线方向,第二个脚标表示应力方向。
流体场内任一点的应力状况,即该点流体微团在任一方向的作用面上的应力,都可以用通过该点的三个相互垂直的作用面上的九个应力分量来表示。
2. 以应力表示的运动微分方程
在黏性流体中取一边长为dx,dy,dz的长方体。各表面应力的方向如图所示。为清晰起见,其中两个面上的应力符号未标。各应力的值均为代数值,正直表示应力沿相应坐标系的正向,反之亦然。由于流体不能承受拉力,因此,
pxxpyy,pzz
必为负值。
由牛顿第二定律,x方向的运动微分方程为:
Xdxdydzpxxdydz+[-(pxx
yx
pxx
x
dy)dydz]+
yxdxdz+[-(yx
y
zx
dy)dxdz]+
zxdxdy+[-(zx
等式两边分别除以
z
dz)dxdy]dxdydz
dux
dt
,然后分别对x,y,z求偏导,得到:
X
(1)
1
PXX
X
1
(
yx
y
duxzx
)zdt
同理,在y方向,由牛顿第三定律得:
Ydxdydz
p
yy
dxdz[(p
yy
p
y
yy
dy)dxdz]
dydz[xy
xy
zy
xy
x
dx)dydz]
zy
dxdy[(zy
dt
y
z
dz)dxdy]
dudxdydz
等式两边同时除以
,然后分别对x,y,z求偏导得:
1
Y
1
p
yy
y
(
zy
z
xy
x
)
dt
y
(2)
同理,在z方向,由牛顿第二定律得:
zdxdydz
p
zz
dxdy[(
p
zz
p
z
zz
dz)dxdy]
dydz[(xz
xz
yz
xy
yz
dx)dydz]dy)dxdz]
dxdz[(yz
yz
dxdydz
等式两边分别除以
dt
z
,并分别对想x,y,z求偏导得:
zz
Z
1
p
z
1
(xz
x
yz
y
)
dt
z
以应力表示的黏性流体运动微分方程的推导
1. 黏性流体的内应力
黏性流体在运动时,表面力不仅有法向应力,还有切向应力,因此黏性流体的表面力不垂直于作用面。
如在任一点取一微小的正六面体,如图所示,作用在平面ABCD上的力
有法向应力
pxx
,与切向应力xy
和xz。应力符号的第一个字母表示作
用面的外法线方向,第二个脚标表示应力方向。
流体场内任一点的应力状况,即该点流体微团在任一方向的作用面上的应力,都可以用通过该点的三个相互垂直的作用面上的九个应力分量来表示。
2. 以应力表示的运动微分方程
在黏性流体中取一边长为dx,dy,dz的长方体。各表面应力的方向如图所示。为清晰起见,其中两个面上的应力符号未标。各应力的值均为代数值,正直表示应力沿相应坐标系的正向,反之亦然。由于流体不能承受拉力,因此,
pxxpyy,pzz
必为负值。
由牛顿第二定律,x方向的运动微分方程为:
Xdxdydzpxxdydz+[-(pxx
yx
pxx
x
dy)dydz]+
yxdxdz+[-(yx
y
zx
dy)dxdz]+
zxdxdy+[-(zx
等式两边分别除以
z
dz)dxdy]dxdydz
dux
dt
,然后分别对x,y,z求偏导,得到:
X
(1)
1
PXX
X
1
(
yx
y
duxzx
)zdt
同理,在y方向,由牛顿第三定律得:
Ydxdydz
p
yy
dxdz[(p
yy
p
y
yy
dy)dxdz]
dydz[xy
xy
zy
xy
x
dx)dydz]
zy
dxdy[(zy
dt
y
z
dz)dxdy]
dudxdydz
等式两边同时除以
,然后分别对x,y,z求偏导得:
1
Y
1
p
yy
y
(
zy
z
xy
x
)
dt
y
(2)
同理,在z方向,由牛顿第二定律得:
zdxdydz
p
zz
dxdy[(
p
zz
p
z
zz
dz)dxdy]
dydz[(xz
xz
yz
xy
yz
dx)dydz]dy)dxdz]
dxdz[(yz
yz
dxdydz
等式两边分别除以
dt
z
,并分别对想x,y,z求偏导得:
zz
Z
1
p
z
1
(xz
x
yz
y
)
dt
z