基于MATLAB的油膜轴承载荷特性研究

《基于MATLAB的油膜轴承载荷特性研究》由会员分享，可在线阅读，更多相关《基于MATLAB的油膜轴承载荷特性研究（130页珍藏版）》请在装配图网上搜索。

1、 本科生毕业论文目 录摘要3Abstract4第一章 绪论51.1 现代油膜轴承的发展概况及应用51.2 研究油膜轴承轴瓦弹性变形的必要性61.3 MATLAB程序开发工具简介71.3.1 MATLAB的概述7 1.3.2 MATLAB的特点71.4有限元与有限差分耦合算法的建立81.4.1有限元法的概述81.4.2耦合算法的建立101.4.3差分计算的基本原理111.4.4数值方法121.5 本文的研究内容及研究方法131.5.1 研究内容131.5.2 研究方法和技术路线13第二章 油膜轴承流体动力润滑原理及基本方程142.1 油膜轴承流体润滑理论及油膜形成原理142.1.1油膜轴承流体润

2、滑理论的形成14 2.1.2润滑油膜的形成原理142.2 流体动压润滑的基本方程流体润滑152.2.1研究的假设152.2.2根据流体力学原理推导雷诺方程162.3用于向心滑动轴承中的雷诺方程192.4雷诺方程的求解和边界条件的确定方法202.4.1雷诺方程的求解202.4.2边界条件的确定方法20第三章不考虑弹性变形时油膜轴承的润滑性能研究243.1二维雷诺方程的无量纲化243.2差分计算的基本原理253.3差分法求解二维雷诺方程的原理253.3.1求解区域离散化263.3.2偏微分方程的离散化263.3.3逐点松弛迭代法283.3.4收敛准则293.3.5雷诺边界条件的引入293.4 油膜

3、轴承无量纲雷诺方程的求解293.4.1 无量纲雷诺方程的求解过程293.4.2 油膜轴承无量纲油膜压力分布303.5 不考虑弹性变形时油膜轴承静特性分析323.5.1无量纲油膜合力的计算与分析323.5.2偏位角的计算与分析333.5.3无量纲润滑油量的计算与分析343.5.4无量纲摩擦力的计算与分析37第四章 有限元法及其在考虑弹性变形时油膜轴承润滑性能研究中的应用414.1有限元法的一般原理和表达格式414.1.1单元位移模式及插值函数414.1.2应变矩阵和应力矩阵434.1.3利用最小位能原理建立有限元方程444.1.4刚度矩阵的性质及边界条件的引入454.2油膜轴承弹性变形数学模型的

4、建立及求解464.2.1模型的建立及边界条件464.2.2单元刚度矩阵的建立484.2.3模型的求解过程534.3 单元网格的划分和计算成果的整理554.3.1 单元网格的划分554.3.2 计算成果的整理55第五章 油膜轴承性能的研究结果分析565.1研究参数的确定565.2弹性变形对油膜轴承润滑性能的影响585.2.1弹性变形对油膜轴承油膜压力分布的影响585.2.2弹性变形对油膜轴承油膜压力峰值和承载能力的影响615.2.3弹性变形对油膜轴承最小油膜厚度的影响625.2.4弹性变形对油膜轴承摩擦力及耗油量的影响635.2.5弹性变形对油膜轴承的偏位角的影响645.3 其他参数对油膜轴承性

5、能的影响655.3.1 有关承载系数的曲线655.3.2 不同速度, 不同间隙对润滑性能的影响665.3.3 轴承温度场对性能参数的影响675.3.4 黏度对轴承性能的影响685.3.5 等直径不同宽径比对轴承润滑性能的影响69第六章 MATLAB软件的应用716.1引言 716.2基础知识71 6.2.1数据和函数的可视化74 6.2.2计算结果的可视化756.3 程序设计766.3.1 M文件简介及其程序结构776.3.2 程序流控制796.3.3 程序的调试和评述79结论83致谢84参考文献85附录86油膜轴承性能数值计算MATLAB程序1不考虑轴瓦弹性变形情况下油膜轴承性能数值计算程序

6、2考虑轴瓦弹性变形情况下油膜轴承性能的数值计算程序外文翻译116基于MATLAB的油膜轴承载荷特性研究Study Basing on the MATLAB about Loading Property of Film Bearing专业：机械设计制造及其自动化（机电）班级：机自0227班 学号：02424112作者：指导教师： 职称：副教授摘 要以往人们大多是在假设轴瓦为刚性的前提下进行油膜轴承的研究和设计。而现代油膜轴承常在高速、重载条件下运行，产生很高的油膜压力，加上由于薄壁瓦的兴起和轴瓦材料向低弹性模量方向的发展，使得轴瓦的弹性变形变得不容忽视。通过对一系列油膜轴承计算研究，分析了在考虑

7、弹性变形和不考虑弹性变形两种情况下润滑性能的变化。利用MATLAB软件对油膜轴承载荷特性进行研究，掌握油膜轴承工作原理、大型软件MATLAB编程及程序开发、数值计算方法有限差分法和有限元方法的原理及相关算法的编写，学习模型数据的前后处理，掌握有关润滑及摩擦等相关知识的应用。此次设计了解了油膜轴承的工作原理类型及结构，设计出油膜轴承载荷特性的计算模型并无量纲化，计算刚性情况下的油膜轴承载荷情况特性和考虑弹性变形的油膜轴承相关特性。此外，还考虑了黏度，温度，速度，相对间隙等参数对油膜特性的影响。研究结果表明，在大偏心率的情况下，弹性变形对油膜轴承润滑性能的影响不容忽视。关键词:油膜轴承、MATLA

8、B、数值计算、有限元法、耦合算法AbstractIn the past, the research and design of film bearing were taken on the premise of assuming the bearing is rigid. However, due to more and more film bearings nowadays operating under conditions of high speed and heavy load would produce higher oil film pressure than before, as

9、 well as the prevalence of thin bearing shell and the developing tendency of lower elasticity modulus of bearing material, the elastic deformation of film bearing can no longer be neglected.Though a bank of film bearings calculation research, the paper analysis lubricating property of film bearing i

10、n the condition of considering and without considering elastic deformation. Using MATLAB researches loading property of film bearing, in order to marster film bearing operating principle, program and program development of large software MATLAB, learn principle and compiling correlation algorithm of

11、 numerical calculation method - the finite difference method and finite element method, study model data fore treatment and after treatment, apply knowledge about lubrication and friction.The paper introduces operating principle type and structure of film bearing, designs film bearing loading proper

12、ty calculation model and non-dimensionalize, calculates character of film bearing in the condition of considering and without considering elastic deformation. otherwise, the paper bids considering the parameters of viscosity, temperature, speed, relative space etc. effect on film property.The analys

13、is results show that the influences of elastic deformation of bearing on the performance of film bearing should no longer be neglected, especially in the case of high eccentricity ratio.Key words: film bearing; MATLAB; numerical calculation; finite element method; coupling algorithm第一章 绪论1.1现代油膜轴承的发

14、展概况及应用油膜轴承，也叫液体摩擦轴承和理想滑动轴承，它是滑动轴承的一种，是一种主要零件的加工精度、表面粗糙度以及各种相关参数（包括润滑油及载荷等）的匹配都非常理想的滑动轴承，因此，在工作时其轴与轴承的工作区域形成一个完整的压力油膜，使金属脱离接触，造成纯液体摩擦，故而称为油膜轴承。油膜轴承作为回转轴支承元件在机械领域的应用十分广泛。与滚动轴承相比，油膜轴承工作平稳、可靠、噪音较低。如果能够保证充分的液体润滑，使得滑动表面被润滑油分开而不发生直接接触，则还可以大大减少摩擦损失和表面磨损，甚至消除磨损。另外，润滑油膜同时还具有一定的吸振能力，这对提高轴承运转的稳定性和运转精度都是十分有益的。流体

15、动压润滑，就是依靠被润滑的一对固体摩擦面间的相对运动，使介于固体摩擦面间的流体润滑膜内产生压力，以承受外载荷而免除固体相互接触，从而起到减少摩擦阻力和保护固体摩擦表面的作用。自从B. Tower在其著名实验中发现了动压现象，继由O.Reynolds分析了动压润滑的机理并导出了描述润滑膜压力分布的微分方程，即著名的雷诺方程，遂奠定了流体动力润滑理论的原始基础。此外，流体动力润滑理论中还有其它的一些方程，如：流动的连续性方程、润滑剂的状态方程(粘度和密度方程)、表面的弹性方程、以及能量方程等。当然，应用最广的还是雷诺方程及其在各种具体条件下的变形形式，以及它们的求解。由于这些理论的建立，使得滑动轴

16、承的研究取得了很多成果，促进了其在实际中的应用。国外研究的比较早，大约在60年代H. D. Conway和H. C. Lee引用Winkler假定的弹性位移方程对轴瓦的弹性变形进行简化，研究了装在刚性轴承座中很薄的圆筒形径向轴瓦表面弹性位移的影响。上世纪80年代，美国学者孙大成对表面变形的影响给出了专门的论述。我国开始对这一问题的研究是始于80年代初的，国内学者中较系统的论述这一联立方程组解法的是张冠坤和钟洪，他们对这一问题的提出开始是基于对动静压轴承的研究，因为动静压轴承的静压供油压力有的已达100MPa，平均压力也达到17MPa，甚至达到23MPa，在这种情况下必须考虑轴承的弹性变形。西安

17、交通大学的梭志文等讨论了弹性流体动力润滑理论中应用权余值法的有限元分析，对Boussinesq积分推荐了一个适用性强、使用方便的子域等参变换，并以有限长的径向轴承为例，采用矩形八节点等参有限元，获得了刚性解和重载下的弹性解。随着新型轴承材料的出现，如弹性金属塑料瓦(EMP)，这些材料的弹性模量较金属低一个数量级，因此在对于由这些材料所形成的径向滑动轴承的流体动力润滑研究时，必须考虑轴瓦的弹性变形对油膜厚度的影响作用。上海大学的张国贤等对这方面的问题进行了专门的研究。近年来，尽管也有一些学者对这方面的问题作了研究，但仍然存在很多的不足。油膜轴承，作为滑动轴承的一种，是个庞大的家族，其应用范围很广

18、，被应用在轧机上作为工作辊轴承（两辊轧机）、支承辊轴承（四辊轧机）的称作轧机油膜轴承。轧机油膜轴承基本上属于低速重载、中速重载轴承，但油膜轴承不只适合应用在轧机上，还适合应用在低速度、大载荷、多重集、工作环境差、能源消耗低及安全程度高的其他重型机械上，如冶金机械、矿山机械、起重运输机械、建筑机械等，尤以固定设备为宜。油膜轴承是现代化轧机关键核心部件之一。 主要用在大型钢板轧机轧锟轴承和高速线材轧机的粗、中连轧机和预精轧、精轧机组的轧锟轴承上 。随着科学技术的不断发展，我们对轧制产品的质量要求越来越高，轧制速度急剧增加，如高速线材已超过100m/s级。油膜轴承由于具有滚动轴承所不具备的承载能力大

19、(比滚动轴承一般大3倍以上)、使用寿命长(理论寿命为1015年)、速度范围宽、抗冲击能力强等特点，因此在轧制行业的应用越来越广。1.2研究油膜轴承轴瓦弹性变形的必要性 通过资料查询，笔者注意到，人们在早期的油膜轴承研究中，在对雷诺方程进行求解时(包括解析解法和数值解法)，往往都是基于组成轴承间隙的轴颈和轴瓦都是刚性的这一假设为前提的。对于一般的轻载、低速油膜轴承，由于运动过程中产生的油膜压力较小，采用这种理论计算的结果和实际实验结果能够较好的吻合。但是，对于高速，重载油膜轴承，尤其近年来薄壁瓦的流行和轴瓦材料向低弹性模量发展的趋势，使得它在运动过程中产生的高油膜压力对轴瓦产生了较大的弹性变形，

20、如果对这种条件下的雷诺方程求解，仍然忽略弹性变形，则理论计算结果和实际实验结果有较大的出入。笔者还注意到，国内外虽有许多学者开始对考虑轴瓦弹性变形这一问题进行研究，也得出了一些结论。但该问题的研究仍处于不完善阶段。对其研究的方法也有很大差别，导致所得结论也往往不一致，还不能很好地指导实际应用，尤其对实际中的重载、高速轴承的研究更是如此。多数研究者对这一问题的研究不够深入、片面，所得的结论也只停留在理论分析水平上。至今尚未形成能够被广泛承认的理论体系和研究结果。综上所述可见，在对油膜轴承进行理论研究的数值分析中引入轴瓦弹性变形因素是十分必要的。以往在该方面进行的研究较为薄弱，而生产实际对油膜轴承

21、提出的高速、重载问题又对这一问题的研究有着十分迫切的要求，为此，本课题拟订以研究轴瓦弹性变形对轴承静特性性能的影响为主要研究目标，探讨更为科学合理的计算方法，利用计算机对考虑轴瓦弹性变形条件下的油膜轴承动压润滑理论进行更符合实际工况的模拟和求解。争取通过这些数值模拟计算，弄清轴瓦弹性变形等因素对动压润滑状态下工作的油膜轴承性能影响的过程和规律，所得出的有关理论结论将会对该领域的进一步研究和实际应用起指导和推动作用。1.3 MATLAB程序开发工具简介1.3.1 MATLAB的概述实验研究、理论分析和科学计算已经成为当代科技研究中不可或缺的三种主要手段。处于计算机时代的今天，科学计算是以数学模型

22、为基础、以计算机和数学软件为工具进行的模拟研究，它是数学通向其他学科的桥梁，是当今盛行的计算机仿真技术的重要基石。当今科学计算的研究方法是计算机和数学的有机结合，计算机已经成为科学计算比不可少的物质基础。数值计算确实需要理论上的指导，但落脚点必须是计算，学习数值计算就要自己会算，更应该会使用计算机算。在计算机高度发展和普及的今天，科技工作者科技算法意识的建立和计算能力的培养，必须在计算机环境下，在适当的平台上通过实际操作进行。“工欲善其事，必先利其器”。在明确了数值计算必须与计算机相结合的原则后，选择实现“数值计算”的工具计算机软件就显得非常重要。以前，很多人进行过数值计算和Basic，For

23、tran或C语言结合的教学试验，但是除计算机专业学生外，其他专业的学生还必须花大力气去学习计算机语言的编程，两者结合的效果并不理想。从“使用”算法考虑，易学好用的MATLAB软件比C或Fortran等语言平台更适合大部分学习者。同时，要想在科学计算上与国际接轨，就必须学会MATLAB。1.3.2 MATLAB的特点当今世界上流行的30多个数学类软件中，MATLAB语言处于数值计算型软件的主导地位，适用范围几乎涵盖了工程数学的各个方面。自从1984年美国MathWorks公司把它推向市场以来，由于功能强大，语法简单，使用方便，界面友好，很快便受到科技界的普遍欢迎。MATLAB是英文“MATRIX

24、 LABORATORY”（矩阵实验室）的缩写，该软件集计算、绘图和仿真于一身，其特点概括如下：（1）数值计算（Numeric）它具有庞大的数学函数库，可以进行工程数学中几乎所有的数值计算。对科学和工程中的各种数值计算、统计等各类数学问题都能提供方便快捷的检索方法。 （2）符号运算（Symbolic） 它具有对数学公式进行符号处理的功能，像因式分解、公式推导、求函数的导数、幂级数展开及解微分方程组等这类数学问题，只要用一个指令，调用一个M-函数文件，就可以轻而易举地解决。（3）数据可视化（Graphic） 它具有很强的数据可视化能力，可以把计算所得的数据根据不同的情况和以前轻松自如地绘制出二维、

25、三维图形，并可对图形的线条、色彩、视角等加以变换。还可以用不同坐标系进行绘图，绘出各种特殊的几何图形（如直方、柱状、网面等），把数据关系的特征形象化的淋漓尽致。（4）建模仿真（Simulink） 它具有动态系统建模、仿真和分析的集成环境。用户只要在视窗里通过简单的鼠标操作，就可以建立起直观的方框图式系统模型，比用传统软件包从微分方程、差分方程出发进行的建模更直观、方便和灵活。通过改变模块中的参数和模块间的连接方式，观察系统中发生的变化，结合仿真的结果，可对系统进行全面或局部的调试和修改，以求达到理想效果。1.4有限元与有限差分耦合算法的建立1.4.1有限元法的概述有限元法亦称有限单元法或有限元

26、素法。它是将弹性连续体离散成为有限个单元的一种近似数值解法。由于离散后的单元与单元之间只通过节点相互连接，且离散后的单元数目和节点数目都是有限的，所以称这种方法为有限元法。将连续体离散成为有限个单元，并建立描述其形态的各种公式，此项工作称为连续体的离散化。有限单元法的基本思想是将连续的求解区域离散为一组有限个、并且按一定方式相互联结在一起的单元的组合体。由于每个单元能按不同的联结方式进行组合，且单元本身又可以有不同的形状，因此可以模型化几何形状复杂的求解域。有限单元法作为数值分析方法的另一个重要特点是利用在每一个单元内假设的近似函数来分片地表示全求解域上待求的未知场函数。单元内的近似函数通常由

27、未知场函数或及其导数在单元的各个结点的数值和其插值函数来表达。这样一来，在一个问题的有限元分析中，未知场函数或及其导数在各个结点上的数值就成为新的未知量(也即自由度)，从而使一个连续的无限自由度问题变成离散的有限自由度问题。一经求解出这些未知量，就可以通过插值函数计算出各个单元内场函数的近似值，从而得到整个求解域上的近似解。显然随着单元数目的增加，也即单元尺寸的缩小，或者随着单元自由度的增加及插值函数精度的提高，解的近似程度将不断改进。如果单元是满足收敛要求的，近似解最后将收敛于精确解。有限元法的优点是显而易见的，因为离散后的单元都是形态容易了解的标准单元，可为每个单元单独建立方程，并可用有限

28、个参数加以描述。而整个结构是由有限个数目的单元所组成的。将有限个单元的方程汇集起来，称作单元的组集或集合，也可以用有限个参数来描述，其基本方程是一组代数方程组。在数学上讲，就是把微分方程的连续形式转化为描述等效集合体形态的代数方程组，以便于进行数值解。从数学角度来看:有限元法是求解数学物理方程的一种数值方法。它是各种经典数值方法，如瑞利一里兹法、伽辽金法等的新形式。有限元法与经典方法都能把一个连续体的偏微分方程组离散化为等效的代数方程组。从这一点讲，它们的性质是相同的。但有限元法与上述经典方法也有重大差别，其基本差别在于插值函数(或试探函数)的选取方式不同。在经典方法中，是在整个求解域上选取统

29、一的插值函数，并要求该函数在域内和域的边界上均满足一定的条件。而在有限元法中，插值函数要分片地分别在子域上或单元上选取，并要求插值函数在各个子域内部、子域之间的分界面上(称内部边界)以及子域与外界分界面上(称外部边界)均满足一定条件。由于有这种差别，使得有限元法的实用价值远远超过了经典方法。简言之，有限元方法的基础是变分原理和加权余量法，将邻域分割成有限个的部分邻域（有限元素），在代表这些部分邻域值的点（节点，多取元素的顶点）中，导出有关未知函数的连立一次方程，并解析这个方程。其根本原理是，在弹性体是用最小势能原理、虚功原理等，或在更为一般的意义上为变分原理。因为能量积分是在整个邻域进行的，所

30、以需要考虑对象物体或问题的整个邻域。将邻域划分为部分邻域，将有关每个部分邻域的关系加起来作为整个邻域。在每个部分邻域函数值的变动，除满足最小限条件外可自由变化。我们知道这一部分邻域划分得越小，便越接近原来问题的正确解答。主要思路做法归纳如下：1）、建立积分方程，根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值问题等价的积分表达式，这是有限元法的出发点。2）、区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作，这部分工作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外，还要表示节点的

31、位置坐标，同时还需要列出自然边界和本质边界的节点序号和相应的边界值。3）、确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函数作为单元基函数。有限元方法中的基函数是在单元中选取的，由于各单元具有规则的几何形状，在选取基函数时可遵循一定的法则。4）、单元分析：将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近；再将近似函数代入积分方程，并对单元区域进行积分，可获得含有待定系数(即单元中各节点的参数值)的代数方程组，称为单元有限元方程。5）、总体合成：在得出单元有限元方程之后，将区域中所有单元有限元方程按一定法则进行累加，形成总体有限元方程。6）、边界条件

32、的处理：一般边界条件有三种形式，分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。对于自然边界条件，一般在积分表达式中可自动得到满足。对于本质边界条件和混合边界条件，需按一定法则对总体有限元方程进行修正满足。7）、解有限元方程：根据边界条件修正的总体有限元方程组，是含所有待定未知量的封闭方程组，采用适当的数值计算方法求解，可求得各节点的函数值。1.4.2耦合算法的建立本论文所提到的耦合算法是指两种不同的数值计算方法相互嵌套，相互影响，不断的迭代直至给定的收敛精度。由于不同的数值计算方法有其适用的计算范围和计算对象，当采用某种数值计算方法处理给定问

33、题时，如果求解出来的结果反过来会影响求解时所用的已知参数，并且在反求时所采用的是另外一种数值计算方法，在这种情况下就应当考虑采用耦合算法。在具体的计算过程中，耦合算法并不是简单的先采用一种数值方法求出结果，然后采用另外一种数值方法去求解结果对已知参数的影响，这没有达到真正意义上的耦合。真正的耦合是指在采用一种数值方法根据已知参数计算结果的每一步迭代过程中，都要考虑采用另外一种数值方法计算此步迭代结果对已知参数的影响，两种计算几乎同步交错进行，最终达到收敛精度。这种计算方法能够合理的反映数学模型的求解过程，且所得到的结果也比较准确。本章在计算轴瓦弹性变形对轴承性能的影响时，建立了耦合算法。其实质

34、是在有限差分法的基础上，将有限元法融入有限差分法的每一迭代过程中。采用有限差分法计算油膜轴承油膜压力的分布的过程，是通过不断的迭代，最终收敛到给定的精度，从而得到所需的结果。但在考虑弹性变形以后，由于压力分布导致弹性变形，而弹性变形反过来又会影响压力分布。对于这种问题，采用上述的耦合算法具有明显的优点，即:在用有限差分法求解油膜压力分布的每一迭代过程中，都计入了采用有限元法计算出的本次迭代压力结果对弹性变形(油膜厚度)的影响，这样两种数值计算方法不断的循环迭代，最终达到收敛，从而求得油膜压力值和轴瓦弹性变形值。1.4.3差分计算概述所谓差分法，就是将弦的两个端点的坐标代入曲线方程，再将所得的两

35、个方程相差、分解因式，便可与直线的斜率以及弦的中点联系起来了，我们称这种方法为差分法。差分法有它的局限性，仅适用于与弦中点、弦斜率同时有关的问题，不及韦达定理法的一般性。差分法较韦达定理的优势是减少了运算量，同时也就增加了准确度。仅凭这一点就足以体现我们研究它的价值，这将为提高解题能力开辟了一条途径。在差分法中，将连续函数用一些结点上的函数值来代替，并从而将微分方程及其边界条件交换为差分（代数）方程，使问题易于求解。在这种方法中，采用了将函数离散的手段。差分法的优点是：l差分法是解微分方程边值问题和弹性力学问题的有效方法。我们总可以将微分方程化为差分方程并得出其数值解答； 2差分法简便易行，且

36、借助于计算机以取较密的网格进行分析，以求得足够精确的解答。例如，对于矩形薄板的弯曲、稳定和振动等问题，可以用差分法很方便地得出解答； 3对于某些结构，为了更精确地分析局部的应力状态，可以用差分法进行分析。例如，若用结构力学方法计算出刚架结构的整体内力分布后，则可用差分法进一步分析刚架结点附近的局部应力状态。差分法的缺点是：1对于由曲线边界和斜边界等产生的不等间距网格，虽然可以得出相应的不等间距的差分公式，或者可改造成为等间距的网格进行分析，但都比较麻烦和易于出错；2差分法比较适用于解二维问题或平面问题，这时的网格较为直观，易于图示；3差分法比较适用于等间距网格，对于应力等变化较为剧烈时，需采用

37、二次网格进行计算。4凡是近似解，在进行求导运算时会降低精度。由于计算机技术的飞速发展，对连续的的数学模型，数值计算的解需要离散化，即变成差分方程求解。差分方程就是含有取离散值得单变量的函数及其差分的方程，也可以将其视为迭代关系： （1）或者记为 （2）这里假定x与f均为一元函数，这种最简单的情形称为一阶自治的或定常的差分方程，有时为了数学上讨论的需要，将其描述为如下的函数迭代的离散过程： （3）简称为一维离散的动力学系统。 1.4.4数值方法数值方法是指数值计算面向数学问题适合于计算机使用的方法，它决定了计算量的大小及计算精度。在基本的“计算方法”中已经知道了一些常用的数值算法，可概括为以下几

38、种：1、迭代法迭代法即逐次逼近法，它指按同一公式重复计算的一个数值过程。在同一问题中可产生各种不同的迭代法，通常在这些方案中选取收敛较快的方法。2、 直代曲法将非线性问题线性化，也就是在一个局部范围中用直线近似代替曲线，即用多项式逼近复杂函数。其中又有牛顿法和梯形法。 3、 整为零法将整体问题分割为若干小部分处理，类似于积分化为和的形式又称为复化梯形法，只要划分的步长足够小，就可以达到预期的精度。4、外推法对于依赖于步长的数值计算公式，为提高其计算精度可用外推法，即把步长增大数倍于原来函数作加减，忽略高次项，可达到计算要求。1.5 本文的研究内容及研究方法1.5.1 研究内容1、 了解了油膜轴

39、承的工作原理类型及结构；2、 设计出油膜轴承载荷特性的计算模型并进行无量纲化；3、 计算刚性情况下的油膜轴承载荷情况特性，采用有限差分法求解刚性轴承的各项性能，并将所得结果与有关经典结论进行对比，以验证所采用算法的可靠性和准确性；4、 以刚性油膜轴承的数值计算作为研究的基础进行外延，研究考虑弹性变形情况下的数学模型，以较快的收敛速度和较高的收敛精度为目标研究构造新的数值算法来求解轴瓦的弹性变形，计算考虑弹性变形的油膜轴承相关特性；5、 计算了速度，相对间隙，黏度，温度等参数对油膜特性参数的影响； 6、 对己经确定的算法，选择合适的语言进行编程，进行计算机数值求解；7、 对计算机求解出的数值结果

40、进行整理、分析，对比在考虑和不考虑弹性变形两种情况下轴承性能，得出有关结论。1.5.2 研究方法和技术路线1、通过对考虑弹性变形对轴承性能的影响进行专门、系统的研究，可以进一步丰富油膜轴承的润滑理论体系，为该领域的深入研究提供理论帮助和指导；2、通过对考虑弹性变形对轴承性能影响进行定量的分析，从数值计算结果中可以得出一系列与实际吻合的结论，对实际应用中关于油膜轴承的结构设计、材料选择以及制造土艺具有理论上的指导意义；3、通过对油膜轴承及流体动力润滑理论的探讨，认为研究轴瓦弹性变形对油膜轴承性能的影响具有较大的理论和实际意义。通过对国内外研究概况的分析，认为目前虽然对这一领域育一定的研究，但研究

41、方法和研究思路仍有一定的局限性。以此为出发点，提出并建立一套比较合理的理论、算法，进而系统计算、分析和研究轴瓦弹性变形对油膜轴承性能的影响。第二章 油膜轴承流体动力润滑原理及基本方程2.1 油膜轴承流体润滑理论及油膜形成原理2.1.1油膜轴承流体润滑理论的形成润滑的目的是用某种润滑材料构成的一层润滑膜把作相对运动的两个表面分开，以不损坏这些表面，并使滑动过程中表面间的摩擦阻力尽可能地减小。对于油膜轴承来说润滑过程和润滑的机理十分关键，决定了其工作性能，是生产实践中必须要加以研究的问题。如第一章所述，工程实际中的润滑一般可以分为固体润滑、边界润滑、流体润滑、混合润滑，在这些所有的润滑方式中以流体

42、润滑最为常用、润滑的效果也最好，油膜轴承的润滑方式以流体润滑为主。流体润滑方式又包括流体动压润滑、流体静压润滑和弹性流体动压润滑三种基本形式，本论文研究的重点在于考虑轴瓦弹性变形的情况下对流体动压润滑状态下工作的油膜轴承的静特性进行研究。关于动压流体润滑的研究已有很长的历史，粘滞流动的基本定律由牛顿1687年提出，但人们对于润滑机理的研究直到十九世纪末才有所突破。1883年英国博钱托尔(Beauchamp Tower)对火车轮轴的轴承进行实验研究，发现了轴承中有流体动压现象存在，同一年，俄国的彼得洛夫发表了同心圆柱体间的摩擦及润滑问题的论文，他的论文的研究结论至今仍常在润滑问题研究中被引用。1

43、886年雷诺(O.Reynolds)应用流体力学原理提出了流体动压润滑理论，为现代轴承研究建立了理论基础。此后，索默菲尔德(A.Sommerfeld )在1904年对无限宽圆柱轴承求得了雷诺方程的解析解，米歇尔(A.GM.Michell )在1929年对无限窄轴承提出了近似解，沃克维克(F. W Ocvirk)在1953年更进一步发展了窄轴承理论，使之在工程上获得了实际应用。特别是从上世纪60年代以来，摩擦学作为一门独立的学科出现，加上相关的理论体系不断完善和计算机技术的发展和应用，使得油膜轴承的研究快速发展，形成了较为系统的理论体系。本章将根据后续章节研究内容的需要，对有关油膜轴承流体动压润滑的理论进行归纳、分析，为下一步研究建立理论基础。2.1.2润滑油膜的形成原理