《基于MATLAB的油膜轴承载荷特性研究.doc》由会员分享，可在线阅读，更多相关《基于MATLAB的油膜轴承载荷特性研究.doc（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油膜轴承流体润滑理论的形成14
2、 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油膜轴承弹性变形数学模型的建立及求解464
4、.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 其他参数对油膜轴承性能的影响655.
5、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不考虑轴瓦弹性变形情况下油膜轴承性能数值计算程序2考虑轴瓦弹性变
6、形情况下油膜轴承性能的数值计算程序外文翻译116基于MATLAB的油膜轴承载荷特性研究Study Basing on the MATLAB about Loading Property of Film Bearing专业：机械设计制造及其自动化（机电）班级：机自0227班 学号：02424112作者：指导教师： 职称：副教授摘 要以往人们大多是在假设轴瓦为刚性的前提下进行油膜轴承的研究和设计。而现代油膜轴承常在高速、重载条件下运行，产生很高的油膜压力，加上由于薄壁瓦的兴起和轴瓦材料向低弹性模量方向的发展，使得轴瓦的弹性变形变得不容忽视。通过对一系列油膜轴承计算研究，分析了在考虑弹性变形和不考虑
7、弹性变形两种情况下润滑性能的变化。利用MATLAB软件对油膜轴承载荷特性进行研究，掌握油膜轴承工作原理、大型软件MATLAB编程及程序开发、数值计算方法有限差分法和有限元方法的原理及相关算法的编写，学习模型数据的前后处理，掌握有关润滑及摩擦等相关知识的应用。此次设计了解了油膜轴承的工作原理类型及结构，设计出油膜轴承载荷特性的计算模型并无量纲化，计算刚性情况下的油膜轴承载荷情况特性和考虑弹性变形的油膜轴承相关特性。此外，还考虑了黏度，温度，速度，相对间隙等参数对油膜特性的影响。研究结果表明，在大偏心率的情况下，弹性变形对油膜轴承润滑性能的影响不容忽视。关键词:油膜轴承、MATLAB、数值计算、有
8、限元法、耦合算法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 well as
9、 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 in the co
10、ndition 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 numeric
11、al 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 property calcu
12、lation 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 analysis resul
13、ts 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，Fortran或C语言
23、结合的教学试验，但是除计算机专业学生外，其他专业的学生还必须花大力气去学习计算机语言的编程，两者结合的效果并不理想。从“使用”算法考虑，易学好用的MATLAB软件比C或Fortran等语言平台更适合大部分学习者。同时，要想在科学计算上与国际接轨，就必须学会MATLAB。1.3.2 MATLAB的特点当今世界上流行的30多个数学类软件中，MATLAB语言处于数值计算型软件的主导地位，适用范围几乎涵盖了工程数学的各个方面。自从1984年美国MathWorks公司把它推向市场以来，由于功能强大，语法简单，使用方便，界面友好，很快便受到科技界的普遍欢迎。MATLAB是英文“MATRIX LABORAT
24、ORY”（矩阵实验室）的缩写，该软件集计算、绘图和仿真于一身，其特点概括如下：（1）数值计算（Numeric）它具有庞大的数学函数库，可以进行工程数学中几乎所有的数值计算。对科学和工程中的各种数值计算、统计等各类数学问题都能提供方便快捷的检索方法。 （2）符号运算（Symbolic） 它具有对数学公式进行符号处理的功能，像因式分解、公式推导、求函数的导数、幂级数展开及解微分方程组等这类数学问题，只要用一个指令，调用一个M-函数文件，就可以轻而易举地解决。（3）数据可视化（Graphic） 它具有很强的数据可视化能力，可以把计算所得的数据根据不同的情况和以前轻松自如地绘制出二维、三维图形，并可对
25、图形的线条、色彩、视角等加以变换。还可以用不同坐标系进行绘图，绘出各种特殊的几何图形（如直方、柱状、网面等），把数据关系的特征形象化的淋漓尽致。（4）建模仿真（Simulink） 它具有动态系统建模、仿真和分析的集成环境。用户只要在视窗里通过简单的鼠标操作，就可以建立起直观的方框图式系统模型，比用传统软件包从微分方程、差分方程出发进行的建模更直观、方便和灵活。通过改变模块中的参数和模块间的连接方式，观察系统中发生的变化，结合仿真的结果，可对系统进行全面或局部的调试和修改，以求达到理想效果。1.4有限元与有限差分耦合算法的建立1.4.1有限元法的概述有限元法亦称有限单元法或有限元素法。它是将弹性
26、连续体离散成为有限个单元的一种近似数值解法。由于离散后的单元与单元之间只通过节点相互连接，且离散后的单元数目和节点数目都是有限的，所以称这种方法为有限元法。将连续体离散成为有限个单元，并建立描述其形态的各种公式，此项工作称为连续体的离散化。有限单元法的基本思想是将连续的求解区域离散为一组有限个、并且按一定方式相互联结在一起的单元的组合体。由于每个单元能按不同的联结方式进行组合，且单元本身又可以有不同的形状，因此可以模型化几何形状复杂的求解域。有限单元法作为数值分析方法的另一个重要特点是利用在每一个单元内假设的近似函数来分片地表示全求解域上待求的未知场函数。单元内的近似函数通常由未知场函数或及其
27、导数在单元的各个结点的数值和其插值函数来表达。这样一来，在一个问题的有限元分析中，未知场函数或及其导数在各个结点上的数值就成为新的未知量(也即自由度)，从而使一个连续的无限自由度问题变成离散的有限自由度问题。一经求解出这些未知量，就可以通过插值函数计算出各个单元内场函数的近似值，从而得到整个求解域上的近似解。显然随着单元数目的增加，也即单元尺寸的缩小，或者随着单元自由度的增加及插值函数精度的提高，解的近似程度将不断改进。如果单元是满足收敛要求的，近似解最后将收敛于精确解。有限元法的优点是显而易见的，因为离散后的单元都是形态容易了解的标准单元，可为每个单元单独建立方程，并可用有限个参数加以描述。而整个结构是由有限
