运动微分方程推导

以应力表示的黏性流体运动微分方程的推导

1. 黏性流体的内应力

黏性流体在运动时,表面力不仅有法向应力,还有切向应力,因此黏性流体的表面力不垂直于作用面。

如在任一点取一微小的正六面体,如图所示,作用在平面ABCD上的力

有法向应力

pxx

,与切向应力xy

和xz。应力符号的第一个字母表示作

用面的外法线方向,第二个脚标表示应力方向。

流体场内任一点的应力状况,即该点流体微团在任一方向的作用面上的应力,都可以用通过该点的三个相互垂直的作用面上的九个应力分量来表示。

2. 以应力表示的运动微分方程

在黏性流体中取一边长为dx,dy,dz的长方体。各表面应力的方向如图所示。为清晰起见,其中两个面上的应力符号未标。各应力的值均为代数值,正直表示应力沿相应坐标系的正向,反之亦然。由于流体不能承受拉力,因此,

pxxpyy,pzz

必为负值。

由牛顿第二定律,x方向的运动微分方程为:

Xdxdydzpxxdydz+[-(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]

dudxdydz

等式两边同时除以

,然后分别对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

xy

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方向的运动微分方程为:

Xdxdydzpxxdydz+[-(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]

dudxdydz

等式两边同时除以

,然后分别对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

xy

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


相关文章

  • 用微积分推导简谐运动的位移时间公式
  • 用微积分推导简谐运动的位移时间公式 已知:以弹簧振子作为研究对象,弹簧振子的质量为m ,设所经历的时间为t ,且当t=0时,弹簧振子处于受力平衡位置,弹簧的拉力系数为k (单位为N/m),设弹簧振子的初速度为V 0,试推导弹簧振子位 移s ...查看


  • 爱因斯坦场方程推导过程的逻辑梳理--纪念广义相对论发表100周年
  • 爱因斯坦场方程推导过程的逻辑梳理 ---纪念广义相对论发表100周年 钟双金 (华东师范大学物理系 上海 201100) (收稿日期:20140710) 摘 要:爱因斯坦用了10年时间分5步从狭义相对论提升到广义相对论及引力场方程,本文对此 ...查看


  • 互动启发式教学在流体力学中的应用
  • 互动启发式教学法在<工程流体力学>教学中的应用 宫伟力1 王渊源 常斌2 (1.中国矿业大学(北京)力建学院 2.中国矿业大学(北京)资源学院) 1. 引言 教学方法是完成教学任务的有效手段,是教学过程中重要的组成部分,也是我国 ...查看


  • 互相垂直同频率简谐振动的合成轨迹方程的推导
  • 安徽科技学院学报,2010,24(2):40-43JournalofAnhuiScienceandTechnologyUniversity 互相垂直同频率简谐振动的合成轨迹方程的推导 熊德永,马 慧 1 2 (1.贵州师范学院物理与电子科学 ...查看


  • [椭圆及其标准方程]教学设计(第1课时)
  • 自己收藏的 觉得很有用 故上传到百度 与大家一起分享! <椭圆及其标准方程>教学设计(第1课时) 吉林省松原市第二高级中学 刘东风 邮编 138000 一.教学内容分析 本节课是普通高中新课程标准实验教科书人教A版数学选修2-1 ...查看


  • 传递过程原理读书笔记
  • 传递过程原理读书笔记 第4章 微元平衡法与通用微分方程组 1. 对矢量式作各种推导和运算  应于于正交曲线坐标系--直角坐标系.柱坐标系和球坐标系. 4.1.3连续方程的推导: 先对流动的流体作质量衡算: 质量累积4.1连续方程--质 ...查看


  • 行星齿轮机构传动
  • 行星齿轮传动速比计算方法综述 杨桂香 郭志强 王明海 杨珍 (中国一拖集团有限公司,河南 洛阳 471039) 摘要:以具有代表性的2K-H型行星齿轮传动为例,对行星齿轮传动速比常用的计算方法进行了介绍:分别用行星架固定法.力矩法.速度图解 ...查看


  • 2.1.1椭圆及其标准方程教案
  • 2.1.1椭圆及其标准方程 教学目标: 1.掌握椭圆的定义,掌握椭圆标准方程的两种形式及其推导过程: 2.能根据条件确定椭圆的标准方程,掌握运用待定系数法求椭圆的标准方程: 3.通过对椭圆概念的引入教学,培养学生的观察能力和探索能力: 4. ...查看


  • 矢量控制系统理论基础及其公式推导
  • 矢量控制系统理论基础及其公式推导 目录: 1. 坐标变换理论 2. A-B-C 静止坐标系下的感应电机数学模型 3. 任意转速旋转的d-q 坐标系下的感应电机数学模型 4. α-β坐标系下的感应电机数学模型 5. dq0坐标系下的感应电机数 ...查看


热门内容