有限长径向滑动轴承挤压润滑数值分析matlab程序.docx

2022年有限长径向滑动轴承挤压润滑数值分析matlab程序专业：机械设计及理论班级：2007 学号：3112058014 姓名：方涛涛 有限长径向滑动轴承挤压润滑数值分析(有限差分法）1 雷诺方程此处假定润滑膜具有相似的粘度,径向挤压膜轴承，径向轴承在载荷W作用下形成挤压润滑时，轴心有一定的移动速度,即挤压速度在挤压膜轴承中觉得间隙h是t的函数将轴承表面沿平面展开，ｘ为周向方向,y为轴向方向，并带入，，则Ｒeynoｌds方程变为: 　　 　 　　 　　 　　 　（1)若令： 　 　 　 　　 则量纲一的Rｅｙｎoｌds方程为：（2）以上各式中,d为轴承直径;L为轴承宽度；为偏心率，,ｅ为偏心距,c为半径间隙方程（2）中两个自变量的变化范畴是:，边界条件（Gumbel):1)轴向方向:在边沿Y=１和Y＝-1处，P=0;在Ｙ=0处,2)周向方向：按Gumｂｅl边界条件:即油膜起点处,取Ｐ＝０;油膜终点在发散区,时，P=02 程序框图简朴阐明:采用有限差分法求解Reｙnｏldｓ方程由量纲一化的Reynolds方程相应得出有限差分法的计算方程，对于每个节点可写出一种方程，而在边界上的节点变量满足边界条件,她们的数值已知。

这样,就可以得出一组线性方程组方程数与未知数数目一致，因此可求解用迭代法求解代数方程组,并使计算成果满足一定的精度,最后求得整个求解域上各节点的变量值开始 输入偏心率eps，挤压速度deps，轴承直径d，轴承宽度l，网格划分m，n初始化压力线性方程组的系数矩阵，右端常数向量设定边界条件求有限差分法计算方程的系数ﻩ构建压力线性方程组的系数矩阵及右端向量迭代法解线性方程组绘制油膜压力图Simpson数值积分求解承载力 　绘制挤压速度与承载力的关系曲线结束 3 源程序fｕnction q１（epｓ，dｅps,Bd,m，n)%有限差分法计算有限长径向滑动轴承挤压膜压力分布及承载量%偏心率：ｅps%挤压速度:deps%轴承径宽比:Bd%有限差分网格划分,轴承周向方向网格数:m%有限差分网格划分,轴承轴向方向网格数:nzspan=[-1,1];fispaｎ=[0,pｉ］;％区间范畴dｅlｆｉ=（fisｐan(２)－fisｐan（1))/m;ｆi=fispan(１):deｌfｉ:ｆiｓｐａn(２);ｄｅlz=（zspａｎ(2）-zspaｎ(1）)/n;z=zsｐaｎ（１)：delｚ:ｚspan(2);ndｅps＝length（deps)；%dｅps向量长度%％%%%%%％％%%%%％%%%pf＝zeros(m+1,ｎ＋1);H=zeｒos(（m-1)*（n-１),(ｍ-1）＊(ｎ-1));％压力系数矩阵g=zｅros(m-1,ｎ－１）；%%%%％%%%%%%％％边值条件ｐf(1,:)=0;ｐf(ｍ+1,：）=0;pｆ(:,1)=0;ｐｆ(:，n+１)=0;%％%%％％%%%%％%%%%％％p=pf（[2:m]，[２：n］);%解％%%%％%%%%％%%%%%%求有限差分法计算方程的系数foｒ a=1:ndeps　％取nｄeｐｓ组depｓ（dε/ｄt）值，分别求解相应压力分布和承载力fｏr　ｊ=1:n－1 　　 for ｉ＝1:m－１　　 　 Ａ１＝1;　 　　 B1=Ｂｄ^2； 　　 　 Ｃ1=-3*eps*ｓin(fi(i))／(1+ｅps＊coｓ(fi(i))）; 　　 E1=cｏｓ（fi(i）)*deps（a)/(1+ｅps*coｓ(fｉ(i)))＾３; 　K=2*（Ａ1／ｄeｌfi^２+B1/dｅlz^2);　 　 　 F=E1/K； 　 Ａ＝（A１/delfｉ^２＋C1/（2*delfｉ))/Ｋ； B=(Ａ1／dｅlfｉ^2-C1/(2*delfi))/K； 　 　 Ｃ＝B1/delz^2/Ｋ; 　 　 D＝B1/ｄelz＾2/Ｋ; 　 　 E=１;%％%%％%%%%%%%％%%％％%构建压力系数矩阵 　 ｇ（ｉ,j)=F; ｉi＝(j－1)*(m-1)＋i;　 　Ｈ(ｉｉ，ii)=-E； 　 　 iｆ mｏd(ｉi,（m-1))==1 　　　 　　　 ｇ(i,j）=ｇ(i,j）－Ｂ*pf（１,j);%右端项组装　　 　 　 ｉｆ ｉ＞1 　 　　 　 　 　 Ｈ(iｉ,iｉ－1)＝0;％MAＴLAB中矩阵行列与一般写法相反 　 　 　 　 end　 　 　　 else 　 　　 　 　 H（ii,iｉ-1)=B;%ＭAＴLＡB中矩阵行列与一般写法相反　 　 　 　ｅnd 　 　 　 if　moｄ(ii,(m－1)）＝=0 g（i，ｊ）=g(i,j)-Ａ*pf(m－1，j);％右端项组装 　 if　ｊ1 　　 　 　H(ｉi,ii－ｍ＋1)=D; 　 　 eｌse 　 　　 g(i,ｊ)=g(i,j)-D*pf(i，1）; 　　　 end; 　 　 　　 　 　ｉｆ jexps & count<1000０, x＝A*x１+ｂ; 　 　tｅｍp=x1; x１=x；　　 x=temp;　 　ｃｏunｔ=cｏunｔ+１； 　　　iｆ　x＜０ 　 　x=0;　 ｅnd;end;iｎ_ｊaｃobi=x;diｓp('叠代次数是:＇);coｕｎt4 压力分布图（量纲一）(图1）调用 ｑ１（0.2，1,1．25,３0，10)即偏心率为0.2；挤压速度为1;径宽比1．25；m=30；n=２０图1 压力分布图5 挤压速度dε/dｔ与承载力Fw之间的关系曲线(图2）调用q１(０.5,[1:１０],1．25,3０,10）即偏心率为0．２;挤压速度为1到10,步长为１；径宽。