有限長徑向滑動軸承擠壓潤滑數(shù)值分析matlab程序

1、文檔來源為：從網(wǎng)絡(luò)收集整理.word版本可編輯.歡迎下載支持專業(yè)：機械設(shè)計及理論班級：2007學號：14姓名：方濤濤有限長徑向滑動軸承擠壓潤滑數(shù)值分析(有限差分法)1雷諾方程此處假定潤滑膜具有相同的粘度，徑向擠壓膜軸承，徑向軸承在載荷W作用下形成擠壓潤滑時，軸心有一定y為軸向萬向，并市入x R ,12 R2 dh(1)dt的移動速度，即擠壓速度0在擠壓膜軸承中認為間隙h是t的函數(shù)。dt將軸承表面沿平面展開，x為周向方向，dxRd,貝Reynolds方程變?yōu)椋?一h3-ph3R2p2y若令：yYL/23文檔來源為：從網(wǎng)絡(luò)收集整理.word版本可編輯.則量綱一的Reynolds方程為:2P2P 3

2、 sin PY21 coscos d(1 cos )3 dt（2）以上各式中，d為軸承直徑；L為軸承寬度；為偏心率,e / c , e為偏心距，c結(jié)束0,2為半徑間隙。方程（2）中兩個自變量的變化范圍是：Y1,1,邊界條件（Gumbel）:P1）軸向萬向：在邊緣Y=1和Y=-1處，P=0；在Y=0處，02）周向方向：按Gumbel邊界條件：即油膜起點0處，取P=0;油膜終點在發(fā)散區(qū)，時，P=02程序框圖簡單說明：采用有限差分法求解Reynolds方程。由量綱一化的Reynolds方程對應(yīng)得出有限差分法的計算方程，對于每個節(jié)點可寫出一個方程，而在邊界上的節(jié)點變量滿足邊界條件，他們的數(shù)值已知。這樣

3、，就可以得出一組線性方程組。方程數(shù)與未知數(shù)數(shù)目一致，所以可求解。用迭代法求解代數(shù)方程組，并使計算結(jié)果滿足一定的精度，最終求得整個求解域上各節(jié)點的變量值。3源程序文檔來源為 :從網(wǎng)絡(luò)收集整理.word 版本可編輯.歡迎下載支持functionq1(eps,deps,Bd,m,n)%有限差分法計算有限長徑向滑動軸承擠壓膜壓力分布及承載量%偏心率：eps%擠壓速度：deps%軸承徑寬比：Bd%有限差分網(wǎng)格劃分，軸承周向方向網(wǎng)格數(shù)：m%有限差分網(wǎng)格劃分，軸承軸向方向網(wǎng)格數(shù)：nzspan=-1,1fispan=0,pi;%區(qū)間范圍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);%壓力系數(shù)矩陣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);%解%求有限差分法計算方程的系數(shù)fora=1:ndeps%Kndeps組deps(d/dt)值，分別求解對應(yīng)壓力分布和承載力

5、forj=1:n-1fori=1:m-1A1=1;#文檔來源為:從網(wǎng)絡(luò)收集整理.word 版本可編輯.文檔來源為 :從網(wǎng)絡(luò)收集整理.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;%構(gòu)建壓力系

6、數(shù)矩陣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%用雅克比迭代求解壓力系數(shù)矩陣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(疊代次數(shù)是：)；count4壓力分布圖(量綱一)(圖1)調(diào)用q1(0.2,1,1.25,30,10)即偏心率為0.2；擠壓速度為1；徑寬比1.25；m=30；n=20圖1壓力分布圖5擠壓速度d/dt與承載力Fw之間的關(guān)系曲線(圖2)調(diào)用q1(0.5,1:10,1.25,30,10)即偏心率為0.2；擠壓速度為1到10，步長為1；徑寬比1.25；9