1. 引言
目前已经有很多学者开展了关于多相流方面的研究，并提出了不同的构成相稳定性描述模型，比较常见的有：幂律模型 [1] 、Bingham模型[2] [3] 、Herschel-Bulkley模型 [4] [5] 、Casson模型 [6] [7] 和Ree- Eying模型[8] 等。但研究发现，这些模型的提出都是基于对多相流体流变特性的定量描述，相关参数往往由实验分析经验确定，缺乏理论依据，模型存在适用范围的限定；这些模型的设计上主要针对于多相流体动特性及边界处的流变分析，当其应用于多相润滑胶体时，对胶体内部的组分变化与力学描述上没有相关理论研究内容；这些模型都不能有效解释胶体内在非均匀应力源(场)下产生组分变化的原因。
此外，大多数模型构建时都把动力粘度作为一个重要的方程参数，但是对非均布应力作用下的多相润滑胶体而言：由于相漂移引起其内部相结构的变化，显然其动力粘度在应力作用方向上的粘度值是变化的而非一定值；当胶体应用于孔流动或薄膜润滑时，其实时粘度值是很难测定或不可测的；在重载或高速摩擦场合，环境温度与压力的变化将引起胶体粘度值的改变；而且对进入稳态剪切流动的胶体，没有证据表明胶体的粘度与其内部相分布相关。
在文中通过对多相润滑胶体在非均布应力下的内部相分布规律及其影响因素分析，并基于稳态系统的熵最小理论，构建了一个可用于描述稳定工作情况胶体内部构成相分布的二相楔滑流变模型。
2. 尖锋应力下的相漂移
此时如果在平板槽中的
位置对流体施加一尖峰应力扰动
，
，并标记这一特定胞元为
。设未施加尖峰应力时胞元
的初始参数为：
，
，则得到：
由于作用于胞元B相上的应力平衡被破坏，在应力差的作用下将推动胞元中的B组分沿楔形面方向进行漂移。如果不考虑邻近胞元间的影响，B组分的漂移将使得胞元内的
值减小。假设多相润滑胶体的组分配比符合最优设计，即
，漂移的发生使得胞元中
增大，
增大，并最终引起
的减小。
设经过一段时间的相漂移后，当
时胞元内实现新的应力平衡
，则新的漂移抵制角
可求得：
显然：
且
。
设由
胞元处尖锋应力引起的对槽壁处额外横向膨胀抑制力为
，参于漂移胞元达到新稳态后的横向应力差为
，
，则存在如下关系：
(1)
式中，
——为平板槽的宽度
。
事实上，胞元中组分沿楔形向的漂移还对相邻胞元产生了流动方向上的作用力，并引起槽宽方向上胶体流速的梯度变化。如果取：
，
且
，则此尖峰应力下的平板槽流动模型即成为非牛顿流体的平板流动模型。
3. 梯度面应力下的相漂移
假设在胶体流动方向上施加的是一梯度变化面力：
则可等效为作用于坐标
处胞元上的等效尖峰应力源：
(2)
则由这一等效尖峰应力引起的附加横向膨胀抑制力可计算为：
(3)
对上式变形并分类整理得到：
(4)
利用计算机工具求解上述方程即可得到发生内部组分漂移后，胞元中的B组分体积分数在槽宽方向上的分布曲线。
4. 模型的构建与验证
鉴于磁流胶体参数在外磁场下的多变性和其在孔中流变特性的实验性测量困难，就需要选用其它参数可知的胶体进行模型的有效性验证。相比之下，高温润滑脂作为混合胶体的一种，其在管道中的流变特性已经有众多的研究成果与实验数据，且润滑脂本身的数据较为完善，所以这里对仿生胞孔进行了适当的尺度放大，以高温润滑脂在细直管道中的受压流动为例对该模型的有效性分析。
在整个流动过程中，两相间的最大静摩擦系数
即为润滑脂内的最大内摩擦系数。如果不考虑应力的梯度变化，即假定作用于脂柱端面两端的作用力分别为：
和
，则因为润滑脂在稳定流动时其作用力方向与应变抗力均沿轴向方向且与流动方向相反，所以其任意半径处胞元上的受力情况应当满足下式：
(5)
得到：
即：
(6)
采用计算机工具，即可求得在半径
方向上流动发生前润滑脂内膨胀抑制力
的分布曲线，如图1所示。图中
值为负是因为其力的作用方向与坐标系的正方向相反，由图中可见其变化趋势与脂内层间剪应力
在半径方向上的分布梯度相同，且在边壁处达到最大值
。
而将求出的
，即可作出当润滑脂内发生组分漂移且达到稳态流动后，相应的润滑脂胞元中重质相在半径方向上的分布函数
，如图2所示。
Figure 1. Distribution of 
before grease flows
图1. 流动发生前润滑脂内膨胀抑制力
的分布
Figure 2. Distribution of 
when grease flows steadily
图2. 稳态流动润滑脂内胞元中重质相分数
的分布曲线
5. 结论
本章主要通过对多相润滑胶体在非均布应力下的内部相分布规律及其影响因素分析，并基于稳态系统的熵最小理论，构建了一个可用于描述稳定工作情况胶体内部构成相分布的模型。对模型的有效性进行仿真分析。研究得到的主要成果如下：
1) 构建了基于稳态系统的最小熵理论，构建了用于描述多相胶体在非均布应力场中相漂移现象与稳态组分分布的模型。
2) 提出了漂移抵制角
与横向膨胀抑制角
的概念用以表征多相胶体胞元在外载荷作用下发生相漂移的趋势，并进行了其表达式的推导。
3) 以高温润滑脂在管道中的流动为例，应用本模型进行了脂内组分漂移现象的研究。仿真结果表明，应用本模型较好地提示了润滑脂流动过程中的剪切稀化现象与组分分布的变化情况。