润滑理论的matlab求解程序_包括弹流和刚流润滑_matlab

clear %bearingfeature.m % This program is used to calculate the feature of journal bearing without taking the elastic deformation into account % e is the eccentricity % ratio is the ratio of L/D % m is the number of columns % n is the number of rays % S is the sum of difference between P (k) and P(k-1) % T is the sum of P(k) % r is the radius of the bearing journal % W is the total of pressure e=input('enter the eccentricity e=') ratio=input('enter the ratio of L/D ratio=') m=input('enter the number of columns m=') n=input('enter the number of rays n=' ) [l,k]=size(e); %求输入偏心率向量的维数 [l,s]=size(ratio); %求输入宽径比向量的维数 H1=figure; H2=figure; H3=figure; H4=figure; H5=figure; H6=figure; %循环计算在不同的宽径比和不同的偏心率的情况下的无量纲承载能力和偏位角 for j=1:s for i=1:k P=dispressure(e(i),m,n,ratio(j)); %disp(‘Begin to test,press return to ….’); %keyboard; hmin=thickness(m,n,e(i),P); Pmax=pressness(P); [Hr,Hi]=filmrupture(P); [W(i),angle(i)]=totalpressure(m,n,P); [Q1(i),Q2(i),Q3(i)]=flow_quantity(m,n,e(i),ratio(j),P); Fe(i,j)=friction(m,n,P,e(i)); end W; angle; %disp(‘检查W和angle的值。。。’); %keyboard; Pmax; hmin; Hr; Hi; figure(H1); hold on; plot(e,W); title('不同宽径比偏心率和无量纲承载力的关系曲线'); xlabel('偏心率\epsilon'); ylabel('无量纲承载力W'); figure(H2); hold on; plot(e,angle); title('不同宽径比偏心率和偏位角的关系曲线'); xlabel ('偏心率\epsilon'); ylabel ('偏位角\psi'); figure(H3); hold on; plot(e,Fe(:,j)); title('不同宽径比偏心率和无量纲摩擦力的关系曲线'); xlabel ('偏心率\epsilon'); ylabel ('无量纲摩擦力F'); figure (H4); hold on; plot(e,Q3); title('不同宽径比偏心率和润滑油端泄流量的关系曲线'); xlabel('偏心率\epsilon'); ylabel('润滑油端泄流量Q'); figure (H5); hold on; plot(e,Q1); title('不同宽径比偏心率和轴承油膜起始边进油量的关系曲线'); xlabel('偏心率\epsilon'); ylabel('轴承油膜起始边进油量Q1'); figure (H6); hold on; plot(e,Q2); title('不同宽径比偏心率和轴承油膜终止边出油量的关系曲线'); xlabel('偏心率\epsilon'); ylabel('轴承油膜终止边出油量Q2'); end