微纳米间隙流动的边界滑移及其流体动力学研究

摘要

经典流体力学和经典润滑力学均认为流体在固体表面上流动时表面流体分子与固体分子的速度绝对相等，即传统的无滑移边界假设。近年来，随着微纳米测试技术及其相关领域科学技术的飞速发展，人们发现微纳米尺度下的流体间隙流动与宏观尺度下流体的流动问题有着本质区别，尺度效应突现，“边界滑移”就是其中最有代表性的一类问题。所谓边界滑移，指的是固体表面上的流体与固体表面存在相对运动速度。研究表明，在微纳米尺度下的边界滑移对间隙流体的流动特性有着重要甚至是决定性的影响。目前有关边界滑移问题的研究主要集中于发现边界滑移现象和探索滑移的内在机理以及各因素对边界滑移的影响规律，包括试验研究和分子动力学模拟，而有关边界滑移对流体系统流体动力学行为影响的数值分析则相对较少。微纳米间隙流动边界滑移及其流体动力学研究中主要存在两方面的问题：一是目前最常用的线性滑移长度模型认为边界滑移总是存在的，并且边界滑移速度与剪切率成正比，但是许多实验发现在低剪切率下没有滑移，而在高剪切率下滑移呈很强的非线性行为，即线性滑移长度模型常常不能准确描述流体(液体)流动的边界滑移；二是采用极限剪应力滑移模型时，由于流体在边界上的滑移速度和方向都是未知的，导致数值求解上的困难，传统的基于有限差分方法的迭代求解技术计算量太大而且收敛性较差，尤其在二维流动问题中该方法基本不可行。针对上述问题，本文展开了以下几个方面的研究。 基于流变学中的极限剪应力模型，结合参变量变分原理及其有限元参数二次规划方法，本文首次给出了求解球形挤压膜流动边界滑移问题的数值方法，为此类边界滑移问题的计算与分析提供了新方法。此外，通过对数值计算结果的分析，本文提出了当上下表面具有相同滑移性质时球形挤压膜流体动压承载力的拟合公式，为边界滑移的间接实验测量提供了有力的理论分析工具。通过数值计算结果与相关实验数据之间的定量比较和分析发现，极限剪应力滑移模型的理论预报值与实验测量值吻合很好，尤其在高剪切率时。一方面说明极限剪应力滑移模型可以用来准确描述间隙流动的边界滑移，另一方面也说明基于极限剪应力滑移模型的参变量变分原理及其有限元参数二次规划方法可以有效求解边界滑移问题，该方法计算效率高、收敛性好。此外，数值计算还表明边界滑移使挤压膜的流体动压力减小，而且固体表面的滑移性质对挤压膜的流体动压特性有着重要影响。当上下表面同时发生大面积的边界滑移时，挤压膜可能完全失去流体动压承载力，但若其中一个表面不发生边界滑移，则即使另一表面发生理想滑移(表面极限剪应力等于零)，挤压膜也仍能保持一定的承载力，且恰为无滑移时挤压膜承载力的四分之一。 通过对Troian等的分子动力学模拟结果的考察，本文提出了一种包含线性滑移长度模型和极限剪应力滑移模型的非线性滑移模型，从而弥补了前者在高剪切率时的不足和后者在低剪切率时对微小边界滑移的忽略。在小剪切率时，该模型认为滑移长度为常数，与线性滑移模型一致。而在高剪切率时，该模型则和极限剪应力模型一致。基于该非线性模型，对平行板剪切流和挤压膜流动的边界滑移问题进行了数值分析并与他人实验数据进行了比较，发现理论预报与实验吻合很好。 利用极限剪应力滑移模型和相应的参变量变分原理及其有限元参数二次规划技术，本文对各种一维间隙流动(与各种液体滑动轴承相对应，其典型间隙厚度在微米量级)的边界滑移及其流体动力学进行了数值分析。研究发现，当同一固体表面滑移性质处处相同时，边界滑移使流体系统的流体动压效应减小，而且与哪个表面发生运动有关。在运动的固体表面发生大面积边界滑移时，或运动和静止固体表面(两表面具有相同滑移性质)同时发生大面积滑移时，流体系统的流体动压效应剧减或完全消失。而间隙厚度减小、剪切率增加、流体粘性增加以及表面极限剪应力(其大小与表面的材料和液体的种类等有关)的减小均会使边界滑移容易发生或者使边界滑移加剧，这就从边界滑移的角度解释了为何在高速、窄隙(重载)以及轴承两表面为同种材料时轴承更容易失效，也解释了为何不同材料的轴承会有不同的极限转速(高于该转速轴承会失效)。计算分析还发现，边界滑移在大滑滚比时更容易使滑滚间隙流动系统(滚柱轴承、金属轧制润滑系统等属于此类间隙)的流体动压承载力减小，而且在纯滑动时系统可能完全丧失承载力，但纯滚动时则不会完全丧失承载力，这解释了实际当中轴承为何在纯滑动时比纯滚动时更容易发生破坏。 当静止固体表面具有不同滑移性质(复合滑移表面)时，边界滑移对系统的流体动力学影响较为复杂，与表面滑移区域的位置、尺寸、极限剪应力大小以及间隙几何等参数有关，但滑移区域内边界滑移越严重，系统的流体动力学效应反而越好。以滑块轴承和轴颈轴承为例，通过简单的优化计算，不仅从理论上发现了“复合滑移轴承系统”(具有复合滑移表面的轴承系统)比传统无滑移轴承具有更高的流体动压承载能力和更低的表面摩擦系数，而且发现平行或发散间隙时系统仍能获得相当高的流体动压承载力，这打破了经典无滑移理论认为只有收敛间隙才能获得流体动压承载力的观点，为新型轴承系统的设计与制造提供了新思路。 采用极限剪应力滑移模型求解二维间隙流动边界滑移问题时，由于边界滑移的大小和方向都事先未知，使求解变得更加困难。本文提出了分段线性化方法来逼近系统的非线性滑移控制函数，然后结合相应的参变量变分原理和有限元参数二次规划方法，提出了新的二维间隙流动边界滑移问题的数值计算方法。通过数值分析，讨论了在控制方程的线性化技术中所采用的正多边形的边数、多边形种类以及有限单元数目对求解精度的影响，数值算例表明，该方法计算效率高、精度好。使用该方法，对二维间隙流动边界滑移问题及其流体动力学进行了数值分析，发现当静止固体表面具有复合滑移性质时系统的流体动力学效应较好。进一步的分析还发现，二维复合滑移滑块轴承在静止固体表面滑移区为梯形且间隙大致平行时获得最大流体动压承载力，最大流体动压承载力是传统无滑移轴承最大承载力的约 2.5 倍，而表面摩擦系数则降低了50％以上。通过对复合滑移转子-轴承系统的动力学分析发现，与无滑移转子系统相比，复合边界滑移转子系统不仅可以使系统具有更高的承载力和更低的摩擦系数，而且能够提高转子系统的运转稳定性，并且从理论上发现系统会由于复合边界滑移的存在具有一种自我稳定运转的功能。这个发现或许有助于我们探索高性能、高运转稳定性的新型结构。但是对于静止表面滑移性质处处相同的转子-轴承系统而言，边界滑移则会使系统的动力学稳定性降低。 最后在第九章介绍了作者在博士论文前期工作中开展的关于高温(900℃)强磁场处理技术改善镍铝合金力学性能的工作。尽管这部分工作相对独立，但由于工作比较新颖，挑战性和开拓性并存，我们还是简单介绍给读者，以期获得有益的讨论并在以后继续深入研究。本文发现利用强磁场处理技术使得镍铝合金的室温力学性能获得大幅度改善，与未经磁场处理的同种合金相比，抗弯强度提高～75％，拉伸和压缩强度都提高一倍以上，而且试件的断口形貌分析表明，强磁场处理改变了镍铝合金的室温断裂方式，使合金呈现韧性断裂特征。 总之，本文首先提出和发展了极限剪应力滑移模型和非线性滑移模型用于描述边界滑移问题，并与目前较常用的线性滑移长度模型和相关实验数据进行了比较，发现前两者的适用范围远大于后者；其次本文成功解决了二维间隙流动问题的系统方程非线性求解困难；最后本文对间隙流动的边界滑移及其流体动力学分析揭示了许多以前未曾知晓的物理现象和规律，尤其是从理论上发现了复合边界滑移能使微纳米间隙流动系统的流体动力学等综合性能提高。这些对今后边界滑移相关问题的计算、分析及工程应用都具有重要的指导意义。

目录

文摘

英文文摘

声明

1绪论

1.1引言

1.2工程研究背景

1.2.1微纳米尺度下边界滑移研究的重要性(边界滑移影响的尺度效应)

1.2.2实际工程中的边界滑移问题及其潜在应用

1.3流体流动边界滑移问题的研究历史和现状

1.3.1边界滑移问题的早期研究历史简介

1.3.2近期关于流体流动边界滑移问题的实验研究和分子动力学模拟

1.3.3边界滑移的物理模型及边界滑移问题的流体动力学研究

1.4边界滑移相关研究目前主要存在的一些问题

1.5本文的主要研究内容

2计算理论与方法基础

2.1标准二次规划问题

2.2参变量变分原理的特点

2.3非线性杆系结构分析及其参变量变分原理的建立

2.3.1三杆结构的数学模型

2.3.2参变量变分原理的建立过程

2.4一维边界滑移问题中的参变量变分原理及其有限元二次规划求解

2.4.1滑移本构状态方程

2.4.2变分泛函表达式

2.4.3有限元法及其二次规划求解的实施过程

2.5本章小结

3极限剪应力滑移模型和非线性滑移模型

3.1基于极限剪应力滑移模型的平行板Couette流边界滑移分析

3.2基于极限剪应力滑移模型的球形挤压膜系统的边界滑移分析

3.2.1球形挤压膜边界滑移问题求解过程的理论推导

3.2.2上下表面具有相同滑移性质时边界滑移对挤压膜的影响

3.2.2上下表面具有不同滑移性质时边界滑移对挤压膜的影响

3.2.3计算结果与实验结果的比较与分析

3.2.4本节讨论

3.3非线性边界滑移模型

3.3.1平行板Couette流的非线性边界滑移

3.3.2球形挤压膜的非线性边界滑移

3.3.3本节讨论

3.4本章小结

4流体一维间隙流动的边界滑移及其流体动力学分析

4.1滑滚间隙内流体膜边界滑移和其流体动力学

4.1.1模型描述和计算方法

4.1.2滑滚表面具有不同滑移性质

4.1.3滑滚表面具有相同滑移性质

4.1.4本节讨论

4.2金属轧制入口处粘塑性润滑分析

4.2.1金属轧制的润滑模型

4.2.2数值计算结果与分析

4.2.3与实验的比较

4.2.4本节讨论

4.3边界滑移对轴颈轴承流体动力学性能的影响

4.3.1模型描述及计算方法

4.3.2轴颈、轴套表面具有相同滑移性质

4.3.3轴颈、轴套表面具有不同滑移性质

4.3.4本节讨论

4.4边界滑移对滑块轴承流体动力学性能的影响

4.4.1模型描述

4.4.2数值计算结果与分析

4.4.2本节讨论

4.5本章小结

5复合边界滑移对一维间隙流动的流体动力学影响

5.1复合边界滑移对无限宽滑块轴承的流体动力学影响

5.1.1复合滑移滑块轴承的模型描述

5.1.2数值计算结果与分析

5.2复合边界滑移对无限宽圆柱型轴颈轴承的流体动力学影响

5.2.1复合滑移圆柱型轴颈轴承的模型描述

5.2.2数值计算结果与分析

5.3本章小结

6二维间隙流动边界滑移问题的计算方法及其流体动力学分析

6.1二维间隙流动边界滑移问题的计算方法

6.1.1边界滑移控制方程

6.1.2问题的有限元及其二次规划表述

6.1.3数值算例一——平行圆盘挤压膜

6.1.4数值算例二——滑块轴承单面滑移问题

6.2边界滑移对二维间隙流动的流体动力学影响

6.2.1有限宽滑块轴承的单面滑移流体动力学分析

6.2.2有限宽滑块轴承双面滑移时的流体动力学分析

6.2.3初始极限剪应力转换点

6.3本章小结

7复合边界滑移对二维间隙流动的流体动力学影响

7.1复合边界滑移对有限宽滑块轴承的流体动力学影响

7.1.1有限宽复合边界滑移滑块轴承的模型描述

7.1.2滑移区域的优化及计算结果

7.2边界滑移对有限宽轴颈轴承的流体动力学影响

7.2.1有限宽轴颈轴承的模型描述

7.2.2均一边界滑移对轴承流体动力学性能的影响

7.2.3复合边界滑移对轴承流体动力学性能的影响

7.3本章小结

8边界滑移对刚性对称转子—轴承系统的动力特性影响

8.1有边界滑移时流体膜刚度和阻尼的计算方法

8.1.1流体膜压力的计算

8.1.2流体膜刚度和阻尼系数的计算

8.2刚性对称转子—轴承系统的动力稳定性

8.3数值计算结果与分析

8.3.1无限宽轴承

8.3.2有限宽轴承

8.4本章小结

9高温强磁场处理使NiAl-Cr(Mo)-Hf合金改性的初步探索

9.1工程研究背景

9.2实验材料和方法

9.3实验结果与分析

9.3.1力学性能试验

9.3.2弯曲试样的破坏断裂特征

9.3.3弯曲试样的显微组织分析及强磁场作用机理的探讨

9.4本章小结

结论

参考文献

创新点摘要

攻读博士学位期间发表学术论文情况

致 谢