【摩擦学学习笔记1】雷诺方程的推导

雷诺方程的推导

雷诺方程是Navier-Stokes方程的特殊形式。

1 基本假设

  1. 忽略体积力(如重力或磁力等)的作用。
  2. 流体在界面上无滑动,即贴于表面的流体流速与表面速度相同。
  3. 在沿润滑膜厚度方向不计压力的变化。由于膜厚仅几十微米或更小,压力不会发生明显的变化。
  4. 与油膜厚度相比较,轴承表面的曲率半径很大,因而忽略油膜曲率的影响,并用平移速度代替转动速度。
  5. 润滑剂是牛顿流体。这对于一般工况条件下使用的矿物油而言是合理的。
  6. 流动为层流,油膜中不存在涡流和湍流。对于高速大型轴承,可能处于湍流润滑。
  7. 与黏性力比较,可忽略惯性力的影响,包括流体加速的力和油膜弯曲的离心力。然而,对于高速大型轴承需要考虑惯性力的影响。
  8. 沿润滑膜厚度方向黏度数值不变。这个假设没有实际根据,只是为了数学运算方便所做的简化。

对于一般流体润滑问题,以上假设(1)~假设(4)基本上是正确的;而假设(5)~假设(8)是为简化而引入的,只能有条件地使用,在某些工况下必须加以修正。

2 方程推导

采用微元体分析方法。

主要步骤:

  1. 由微元体受力平衡条件,求出流体沿膜厚方向的流速分布;
  2. 将流速沿润滑膜厚度方向积分,求得流量;
  3. 应用流量连续条件,最后推导出雷诺方程的普遍形式。

2.1 微元体平衡

x方向的受力:

微元体X方向的受力

这个微元体只受其他流体微元给它的压力p和由流体黏性造成的切应力tau作用。

三个方向的速度分量u,v,wZ方向的变化率,dfrac{partial u}{partial z}dfrac{partial v}{partial z}比较大,而dfrac{partial w}{partial z}则非常小(可以认为在膜厚方向,即Z方向的速度沿着Z方向基本没有变化)。

X方向的受力平衡,可得:

p mathrm{d} y mathrm{d} z+left(tau+dfrac{partial tau}{partial z} mathrm{d} zright) mathrm{d} x mathrm{d} y=left(p+dfrac{partial p}{partial x} mathrm{d} xright) mathrm{d} y mathrm{d} z+tau mathrm{d} x mathrm{d} y /

将括号里的项展开,并两边同除mathrm{d} xmathrm{d} ymathrm{d} z即可得:

dfrac{partial p}{partial x} = dfrac{partial tau}{partial z} /

由于假设(5)和假设(6),流体满足牛顿黏性定律,tau=eta dfrac{partial u}{partial z},故

dfrac{partial p}{partial x}=dfrac{partial}{partial z}left( eta dfrac{partial u}{partial z}right) tag{1}

再由Y方向的力平衡,同理也可得:

dfrac{partial p}{partial y}=dfrac{partial}{partial z}left( eta dfrac{partial v}{partial z}right) /

由假设(3),有

dfrac{partial p}{partial z}=0 /

将式(1)直接对z进行积分两次,可得:

begin{aligned} eta dfrac{partial u}{partial z}&=int dfrac{partial p}{partial x} mathrm{~d} z=dfrac{partial p}{partial x} z+C_{1} / eta u&=intleft(dfrac{partial p}{partial x} z+C_{1}right) mathrm{d} z=dfrac{partial p}{partial x} dfrac{z^{2}}{2}+C_{1} z+C_{2} end{aligned} /

再用边界条件来确定C_1C_2。由假设(2),界面上的流体速度等于表面速度,设固体表面的速度为U_0U_h,即当z=0时,u=U_0,当z=h时,u=U_h,可得:

C_{2}=eta U_{0}, quad C_{1}=left(U_{h}-U_{0}right) dfrac{eta}{h}-dfrac{partial p}{partial x} dfrac{h}{2} /

因此,润滑膜中任意点沿X方向的流速为:

u=dfrac{1}{2 eta} dfrac{partial p}{partial x}left(z^{2}-z hright)+left(U_{h}-U_{0}right) dfrac{z}{h}+U_{0} /

同理,可以得到润滑膜中任意点沿Y方向的流速为:

v=dfrac{1}{2 eta} dfrac{partial p}{partial y}left(z^{2}-z hright)+left(V_{h}-V_{0}right) dfrac{z}{h}+V_{0} /

流速组成

上表示流速u沿Z向的分布,它由3部分组成:第1项按抛物线分布,表示由dfrac{partial p}{partial x}引起的流动,故称“压力流动”;第2项按线性(三角形)分布,代表由于两表面的相对滑动速度(U_h-U_0)引起的流动,称为“速度流动”;第3项是常数,表示整个润滑膜以速度U_0运动,沿膜厚方向即Z向各点的速度相同。

2.2 连续方程的积分

流体力学连续方程(质量守恒定律):

dfrac{partial rho}{partial t}+left[dfrac{partial(rho u)}{partial x}+dfrac{partial(rho v)}{partial y}+dfrac{partial(rho w)}{partial z}right]=0 tag{2}

将式(2)沿膜厚Z方向进行积分,有:

int_{0}^{h(x, y)} dfrac{partial rho}{partial t} mathrm{~d} z+int_{0}^{h(x, y)} dfrac{partial(rho u)}{partial x} mathrm{~d} z+int_{0}^{h(x, y)} dfrac{partial(rho v)}{partial y} mathrm{~d} z+int_{0}^{h(x, y)} dfrac{partial(rho w)}{partial z} mathrm{~d} z=0 tag{3}

由于被积函数关于各个变量均连续可微,所以可以将式(3)的积分、微分次序交换:

dfrac{partial }{partial t}int_{0}^{h(x, y)} rho mathrm{~d} z+ dfrac{partial}{partial x}int_{0}^{h(x, y)} (rho u) mathrm{~d}z+dfrac{partial}{partial y}int_{0}^{h(x, y)} (rho v) mathrm{~d} z+dfrac{partial}{partial z} int_{0}^{h(x, y)} (rho w)mathrm{~d} z = 0 /

让我们一项一项来查看,第一项积分后:

dfrac{partial }{partial t}int_{0}^{h(x, y)} rho mathrm{~d} z=dfrac{partial }{partial t}(rho h) /

将速度u的表达式代入第二项,进行积分:

begin{aligned} dfrac{partial}{partial x}int_{0}^{h(x, y)} (rho u) mathrm{~d}z &= dfrac{partial}{partial x}int_{0}^{h(x, y)} left(rho left[ dfrac{1}{2 eta} dfrac{partial p}{partial x}left(z^{2}-z hright)+left(U_{h}-U_{0}right) dfrac{z}{h}+U_{0}right] right) mathrm{~d}z / &= dfrac{partial}{partial x} rho left[ dfrac{1}{2 eta} dfrac{partial p}{partial x}left(dfrac{z^{3}}{3}-dfrac{z^2}{2} hright)+left(U_{h}-U_{0}right) dfrac{z^2}{2h}+U_{0}z right]bigg|_{0}^{h(x, y)} / &=dfrac{partial}{partial x} rho left[ dfrac{1}{2 eta} dfrac{partial p}{partial x}left(dfrac{h ^{3}}{3}-dfrac{h ^2}{2} hright)+left(U_{h}-U_{0}right) dfrac{h ^2}{2h}+U_{0}h right] / &=-dfrac{partial}{partial x} left(dfrac{rho h^3}{12eta}dfrac{partial p}{partial x} right)+dfrac{rho h}{2}dfrac{partialleft(U_h-U_0 right)}{partial x}+dfrac{ (U_h-U_0)}{2}dfrac{partial (rho h)}{partial x} + rho hdfrac{partial U_0}{partial x}+U_0dfrac{partial (rho h)}{partial x} / &=-dfrac{partial}{partial x} left(dfrac{rho h^3}{12eta}dfrac{partial p}{partial x} right)+dfrac{rho h}{2}dfrac{partialleft(U_h+U_0 right)}{partial x}+dfrac{ (U_h+U_0)}{2}dfrac{partial (rho h)}{partial x} end{aligned} /

第三项同理,积分后可得:

begin{aligned} dfrac{partial}{partial y}int_{0}^{h(y, y)} (rho v) mathrm{~d} z &= dfrac{partial}{partial y}int_{0}^{h(y, y)} left(rho left[ dfrac{1}{2 eta} dfrac{partial p}{partial y}left(z^{2}-z hright)+left(V_{h}-V_{0}right) dfrac{z}{h}+V_{0}right] right) mathrm{~d}z / &= dfrac{partial}{partial y} rho left[ dfrac{1}{2 eta} dfrac{partial p}{partial y}left(dfrac{z^{3}}{3}-dfrac{z^2}{2} hright)+left(V_{h}-V_{0}right) dfrac{z^2}{2h}+V_{0}z right]bigg|_{0}^{h(y, y)} / &=dfrac{partial}{partial y} rho left[ dfrac{1}{2 eta} dfrac{partial p}{partial y}left(dfrac{h ^{3}}{3}-dfrac{h ^2}{2} hright)+left(V_{h}-V_{0}right) dfrac{h ^2}{2h}+V_{0}h right] / &=-dfrac{partial}{partial y} left(dfrac{rho h^3}{12eta}dfrac{partial p}{partial y} right)+dfrac{rho h}{2}dfrac{partialleft(V_h-V_0 right)}{partial y}+dfrac{ (V_h-V_0)}{2}dfrac{partial (rho h)}{partial y} + rho hdfrac{partial V_0}{partial y}+V_0dfrac{partial (rho h)}{partial y}/ &=-dfrac{partial}{partial y} left(dfrac{rho h^3}{12eta}dfrac{partial p}{partial y} right)+dfrac{rho h}{2}dfrac{partialleft(V_h+V_0 right)}{partial y}+dfrac{ (V_h+U_0)}{2}dfrac{partial (rho h)}{partial y} end{aligned} /

第四项积分后可得:

dfrac{partial}{partial z} int_{0}^{h(x, y)} (rho w)mathrm{~d} z= dfrac{partial}{partial z} (rho wh) approx 0 /

将每项积分代入式(3)并两边同乘12,整理后可得:

begin{aligned} dfrac{partial}{partial x}left(dfrac{rho h^{3}}{eta} dfrac{partial p}{partial x}right)+dfrac{partial}{partial y}left(dfrac{rho h^{3}}{eta} dfrac{partial p}{partial y}right)=& 6left(U_{0}+U_{h}right) dfrac{partial(rho h)}{partial x}+6left(V_{0}+V_{h}right) dfrac{partial(rho h)}{partial y}+/ & 6 rho h dfrac{partialleft(U_{0}+U_{h}right)}{partial x}+6 rho h dfrac{partialleft(V_{0}+V_{h}right)}{partial y}+12 dfrac{partial(rho h)}{partial t} end{aligned} tag{4}

(4)即为一般形式的雷诺方程。

若令U=U_0+U_hV=V_0+V_h,并认为流体密度不随时间变化(但可以随着位置而变化),则雷诺方程可以写成

dfrac{partial}{partial x}left(dfrac{rho h^{3}}{eta} dfrac{partial p}{partial x}right)+dfrac{partial}{partial y}left(dfrac{rho h^{3}}{eta} dfrac{partial p}{partial y}right)=6left[ dfrac{partial}{partial x}(Urho h)+dfrac{partial}{partial y}(Vrho h)+2rho(w_h-w_0)right] /

3 雷诺方程的矢量形式

用于不可压缩流体润滑计算的雷诺方程普遍形式为:

dfrac{partial}{partial x}left(dfrac{h^{3}}{12 eta} dfrac{partial p}{partial x}right)+dfrac{partial}{partial y}left(dfrac{h^{3}}{12 eta} dfrac{partial p}{partial y}right)=dfrac{1}{2} dfrac{partial(h U)}{partial x}+dfrac{1}{2} dfrac{partial(h V)}{partial y}+dfrac{partial h}{partial t} /

写作矢量形式:

nablacdot left( dfrac{h^3}{12eta} nabla pright)=dfrac{1}{2}nabla cdot (h mathbf{U})+dot{h} /

其中,nabla = mathbf{i}dfrac{partial}{partial x}+mathbf{j}dfrac{partial}{partial y}mathbf{U}为速度矢量;dot{h}=dfrac{partial h}{partial t}

3.1 详细的变换过程

nabla = mathbf{i}dfrac{partial}{partial x}+mathbf{j}dfrac{partial}{partial y} = begin{pmatrix} 1 / 0 end{pmatrix}dfrac{partial}{partial x}+ begin{pmatrix} 0 / 1 end{pmatrix}dfrac{partial}{partial y} = begin{pmatrix} dfrac{partial}{partial x} / dfrac{partial}{partial y} end{pmatrix} /mathbf{U}=begin{pmatrix} U/ V end{pmatrix} /nabla p = begin{pmatrix} dfrac{partial p}{partial x}/ dfrac{partial p}{partial y} end{pmatrix} /

则由向量点乘可得:

nablacdot left( dfrac{h^3}{12eta} nabla pright)= begin{pmatrix} dfrac{partial}{partial x} / dfrac{partial}{partial y} end{pmatrix} cdot left(dfrac{h^3}{12eta} begin{pmatrix} dfrac{partial p}{partial x}/ dfrac{partial p}{partial y} end{pmatrix} right) = dfrac{partial}{partial x}left(dfrac{h^{3}}{12 eta} dfrac{partial p}{partial x}right)+dfrac{partial}{partial y}left(dfrac{h^{3}}{12 eta} dfrac{partial p}{partial y}right) /

同理

dfrac{1}{2}nabla cdot (h mathbf{U})=dfrac{1}{2} begin{pmatrix} dfrac{partial}{partial x} / dfrac{partial}{partial y} end{pmatrix} cdot left(h begin{pmatrix} U/ V end{pmatrix} right) = dfrac{1}{2} dfrac{partial(h U)}{partial x}+dfrac{1}{2} dfrac{partial(h V)}{partial y} /

发表回复

相关推荐

香港公司必备知识 | 3分钟了解商业登记证

根据香港《商业登记条例》(第310章)的规定,任何在香港经营业务的人士,都需要在公司开业之后的一个月内,以书面通知向税局 ...

· 3分钟前

被聯合抵制!四川大學走到今天這一步,純屬自找的

“落實立德樹人根本任務有差距,思想政治工作不夠到位,師德師風建設不夠紮實,加強學科建設不夠有力。”“落實全面從嚴治黨兩個...

· 11分钟前

小户型福利,7款迷你衣帽间设计方案,超实用!

其实在心底, 每个女生都希望有一间属于自己的衣帽间(包括丽酱在内)。 我们在很多影视剧里都看过这样梦幻敞亮的衣帽间,简 ...

· 11分钟前

直流電機和交流電機有什麼不同

1、首要兩者的外部供電不同,直流電機運用直流電做為電源;而交流電機則是運用交流電做為電源。2、從結構上說,前者的原理相...

· 13分钟前

转帖 欲加之罪,何患无辞–袁崇焕十二宗罪辨析

这是《袁崇焕之孤城》第二部《孤忠》后记之二(已在起点首发),看到凤凰网“袁崇焕是怎样走上神坛的?”等几组专稿以及起点中 ...

· 14分钟前