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

1、文档来源为：从网络收集整理.word版本可编辑.欢迎下载支持专业：机械设计及理论班级：2007学号：14姓名：方涛涛有限长径向滑动轴承挤压润滑数值分析(有限差分法)1雷诺方程此处假定润滑膜具有相同的粘度，径向挤压膜轴承，径向轴承在载荷W作用下形成挤压润滑时，轴心有一定y为轴向万向，并市入x R ,12 R2 dh(1)dt的移动速度，即挤压速度0在挤压膜轴承中认为间隙h是t的函数。dt将轴承表面沿平面展开，x为周向方向，dxRd,贝Reynolds方程变为：2一h3-ph3R2p2y若令：yYL/23文档来源为：从网络收集整理.word版本可编辑.则量纲一的Reynolds方程为:2P2P 3

2、 sin PY21 coscos d(1 cos )3 dt（2）以上各式中，d为轴承直径；L为轴承宽度；为偏心率,e / c , e为偏心距，c结束0,2为半径间隙。方程（2）中两个自变量的变化范围是：Y1,1,边界条件（Gumbel）:P1）轴向万向：在边缘Y=1和Y=-1处，P=0；在Y=0处，02）周向方向：按Gumbel边界条件：即油膜起点0处，取P=0;油膜终点在发散区，时，P=02程序框图简单说明：采用有限差分法求解Reynolds方程。由量纲一化的Reynolds方程对应得出有限差分法的计算方程，对于每个节点可写出一个方程，而在边界上的节点变量满足边界条件，他们的数值已知。这样

3、，就可以得出一组线性方程组。方程数与未知数数目一致，所以可求解。用迭代法求解代数方程组，并使计算结果满足一定的精度，最终求得整个求解域上各节点的变量值。3源程序文档来源为 :从网络收集整理.word 版本可编辑.欢迎下载支持functionq1(eps,deps,Bd,m,n)%有限差分法计算有限长径向滑动轴承挤压膜压力分布及承载量%偏心率：eps%挤压速度：deps%轴承径宽比：Bd%有限差分网格划分，轴承周向方向网格数：m%有限差分网格划分，轴承轴向方向网格数：nzspan=-1,1fispan=0,pi;%区间范围delfi=(fispan(2)-fispan(1)/m;fi=fispa

4、n(1):delfi:fispan(2);delz=(zspan(2-)zspan(1)/n;z=zspan(1):delz:zspan(2);ndeps=length(deps);%deps向量长度%pf=zeros(m+1,n+1);H=zeros(m1)*(n-1),(m-1)*(n-1);%压力系数矩阵g=zeros(m-1,n-1);%边值条件%pf(1,:)=0;pf(m+1,:)=0;pf(:,1)=0;pf(:,n+1)=0;%p=pf(2:m,2:n);%解%求有限差分法计算方程的系数fora=1:ndeps%Kndeps组deps(d/dt)值，分别求解对应压力分布和承载力

5、forj=1:n-1fori=1:m-1A1=1;#文档来源为:从网络收集整理.word 版本可编辑.文档来源为 :从网络收集整理.word 版本可编辑.欢迎下载支持 .B1=BdA2;C1=-3*eps*sin(fi(i)/(1+eps*cos(fi(i);E1=cos(fi(i)*deps(a)/(1+eps*cos(fi(i)A3;K=2*(A1/delfiA2+B1/delzA2);F=E1/K;A=(A1/delfiA2+C1/(2*delfi)/K;B=(A1/delfiA2-C1/(2*delfi)/K;C=B1/delzA2/K;D=B1/delzA2/K;E=1;%构建压力系

6、数矩阵g(i,j)=F;ii=(j-1)*(m-1)+i;H(ii,ii)=-E;ifmod(ii,(m-1)=1g(i,j尸g(i,jB*pf(1,j);%右端项组装ifi1H(ii,ii-1)=0;%MATLAEfr矩阵行列与一般写法相反endelseH(ii,ii-1)=B;%MATLA即矩阵行列与一般写法相反endifmod(ii,(m-1)=0g(i,j)=g(i,jA*pf(m-1,j)；%右端项组装ifj1H(ii,ii-m+1)=D;elseg(i,j)=g(i,j)-D*pf(i,1);end;ifjn-1H(ii,ii+m-1)=C;elseg(i,j)=g(i,j)-C*

7、pf(i,n-1);endendend%用雅克比迭代求解压力系数矩阵p(:)=in_jacobi(H,g(:),eye(length(g(:),1);pf(2:m,2:n)=p;%将负压各点置零forj=1:n+1fori=1:m+1ifpf(i,j)exps&count10000,x=A*x1+b;temp=x1;x1=x;x=temp;count=count+1;ifx0x=0;end;end;in_jacobi=x;disp(叠代次数是：)；count4压力分布图(量纲一)(图1)调用q1(0.2,1,1.25,30,10)即偏心率为0.2；挤压速度为1；径宽比1.25；m=30；n=20图1压力分布图5挤压速度d/dt与承载力Fw之间的关系曲线(图2)调用q1(0.5,1:10,1.25,30,10)即偏心率为0.2；挤压速度为1到10，步长为1；径宽比1.25；9