动压轴承的压力场分布,要求求解雷诺方程的简化形式为:!!xh3!!px!"+!!zh3!!pz!"=6U"!h!x公式中参数的含义参见流体动力学。1编程求解雷诺方程采用有限差分的五点差分格式,120120网格点计算。采用C++语言编写求解轴承压力分布的程序,可以大大提高运行速度,节省计算时间。把计算的结果存储到txt格式的文件中。存储数据文件的格式为:起始部分为标识文件属性的字符,用于说明当前文件的产生日期和功能等等,准确的区分不同时间、不同版本产生的数据。正文部分为数据,是该数据文件的主要部分。这样作的优点是:能够很清晰地知道当前文件的来源及用途,不至于混淆多个相似文件,方便查找。生成的txt文本格式如下:全圆动压轴承压力计算(ips=0.4,ratio=1)Currenttimeis:03/11/070……00.0001333860.0002632160.0003898980.0005136770.0006347110.0007531080.0008689430.0009822730.001093140.001201580.001307620.00141129……2Matlab编程绘制图形为了快捷地绘制出高质量的图形,采用Matlab图形库编写代码画图,描绘压力场分布情况。文件的起始部分是字符串,Matlab的load函数无法直接读取该文件中的数据。这样就需要编写一个可以读取这样格式文件的函数,判断出字符和数据部分。把数据和轴承系数分别读取出来,字符作为图形的名称,数据描绘成压力场分布图,这样就可以达到可视化的目的,清楚地看到压力的变化。算法:首先对文本逐行读取,存储行号到一个数组中。一旦发现该行是以数字开始的,就停止读取字符,记录字符的最大行号。以此分界,上面是字符,下面是数据。(1)首先读取第一行,判断结果,发现该行不是数据开始的,就增加行号计数器。(2)读取第二行,判断结果,发现该行也不是数据开始的,继续增加行号计数器。(3)下一行读取的数据行,发现该行是数据开始的,记录行号为数据起始行。(4)停止读取字符,读取数据直至文件末尾,存放到数据数组中去,关闭文件。(5)重新打开文件,读取文件到行号计数器的位置,把读取到的字符串赋给字符数组。(6)调整数据数组的行号和列号到真正的大小。这样就能准确地读取数据,并绘制出良好图形,很清晰地看到动压轴承的压力分布状况,下面是压力场分布的黑白图片(输出其实为彩色)。2007-3-11图的效果:从图中可以清晰地看出,动压轴承全部的压力分布状况,沿圆周方向压力递增,形成油膜的压力呈楔形分布,达到高峰后,油膜快速破裂,压力快速降低到0。对压力场进行积分运算就可以知道动压轴承的承载能力,提全圆动压轴承压力场分布(120120)等压线分布状况供轴承制造的参考。3结论(1)读取的数据最好是不要经过编辑,编辑过的数据文件经常读取不成功。用记事本打开过以后,文件的结束符改变了,读取的时候为乱码。(2)算法中打开两次文件,在一次读取多行字符的时候,很可能读取不成功。打开两次就可以保证数据的正确读取。(3)数值计算求解的压力分布情况比实验方法获得图形结果能够更好地反映轴承中的压力分布情况。数值计算绘制动压轴承压力场分布图@闫帅经$烟台职业学院机械工程系!山东烟台264205动压轴承压力场分布复杂,用实验方法求解的是二维平面图形。随着计算机技术的发展,选择合适的计算工具可以快速获得压力场的三维图。采用C++编写计算程序,生成txt文件存放计算结果。编写Matlab程序读取txt文件中的数据和字符,实现压力的分布情况的三维立