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

有限长径向滑动轴承挤压润滑的数值模拟专业：机械设计和理论类别：2007学号：名字：方涛有限长度径向滑动轴承挤压润滑的数值模拟(有限差分法)1雷诺方程式此处假设润滑膜的黏度相同。径向挤压膜轴承，径向轴承在负载w的作用下形成挤压润滑时，轴心具有恒定的移动速度，即挤压速度。在挤压膜轴承中，间隙h是被视为t的函数。沿平面展开轴承曲面时，x方向是主方向，y方向是轴方向，导入时，Reynolds表达式为：(1)对于命令：标注1的Reynolds方程式为：(2)在以上各种中，d是轴承直径。l是轴承宽度。偏心率，e为偏心距离，c为半径间距。在方程式(2)中，两个引数的变更范围如下：边界条件(Gumbel):1)轴方向：边Y=1和Y=-1；p=0；Y=0时，2)圆周方向：基于Gumbel边界条件：油膜起始位置，p=0；如果油膜末端位于发散区域中，则P=02方块图简要说明：雷诺方程用有限差分法求解。与有限差分法相对应的基于尺寸的Reynolds表达式可以为每个节点创建一个表达式，并且已知满足边界条件的边界节点变量的值。通过这种方法，可以得到一组线性方程。方程式数与未知数一致，所以可以求解。用迭代法求解代数方程，并确保计算结果满足一定精度，从而得到整个解域中每个节点的变量值。开始输入偏心EPS、挤压速度deps、轴承直径d、轴承宽度l、网格分割m、n初始化压力线性方程式的系数矩阵，右端常数向量设定边界条件求有限差分法计算方程的系数压力线性方程的系数矩阵和右端向量的构造迭代法求解线性方程组画油膜压力图Simpson数值积分解决方案承载力挤出速度与轴承容量关系曲线的绘制结束3源代码程序Function q1(eps、deps、Bd、m、n)%有限差分法计算有限长度径向滑动轴承挤压膜的压力分布和承载力%偏心率：EPS挤压速度%:deps轴承直径比率%:Bd%有限差分栅格分割，轴承圆周方向栅格数：m%有限差分网格，轴承轴向网格数：nZspan=-1，1；Fispan=0，pi；%部分范围delfi=(fispan(2)-fispan(1)/m；fi=fispan(1): delfi : fispan(2)；delz=(zspan(2)-zspan(1)/n；z=zspan(1): delz : z span(2)；Ndps=长度(deps)；s向量长度%Pf=zeros(m 1，n 1)；H=zeros(m-1)*(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 (23360m，23360n)；%解决方案%求有限差分法计算方程的系数For a=1:ndeps%取ndps组deps(deps/dt)的值，分别解决相应的压力分布和轴承容量For j=1:n-1For i=1:m-1a1=1；B1=BD 2；C1=-3 * EPS * sin(fi(I)/(1eps * cos(fi(I)；E1=cos(fi(I)* deps(a)/(1eps * cos(fi(I)3；k=2 *(a1/delfi 2 B1/delz 2)；f=E1/K；a=(a1/del fi 2 C1/(2 * del fi)/k；b=(a1/delfi 2-C1/(2 * delfi)/k；c=B1/delz 2/k；d=B1/delz 2/k；e=1；%配置%压力系数矩阵G(i，j)=F；ii=(j-1)*(m-1)I；H(ii，ii)=-E；If mod(ii，(m-1)=1G(i，j)=g(i，j)-B*pf(1，j)；%右结束程序集If i1H(ii，ii-1)=0；%MATLAB矩阵行与常规写相反EndElseH(ii，ii-1)=B；%MATLAB矩阵行与常规写相反EndIf mod(ii，(m-1)=0G(i，j)=g(i，j)-A*pf(m-1，j)；%右结束程序集If J1H(ii，ii-m 1)=D；ElseG(i，j)=g(i，j)-D*pf(i，1)；EndIf jexps count10000、x=a* x1 b；Temp=x1x1=x；X=tempcount=count 1；If x0x=0；EndEndin _ jacobi=x；Disp(重复次数为：)；Count4压力分布图(尺寸1)(图1)调用Q1(0.2，1，1.25，30，10)偏心率为0.2。挤压速度为1。纵横比1.25；M=30N=20图1压力分布图5挤压速度d/dt和承载力Fw