xwzcsu
2009-03-10, 16:32
程序如下:
v=0.01:0.01:10;
delta_p=[8000000 15000000 25000000 35000000];
Cs=0.8e-9;
u=0.051;
n=2200;
Cv=0.2e+6;
Cf=0.01;
Tc=5;
Vmax=45e-6;
nc=0.9;
rd=0.3223;
for i=1:1:length(delta_p)
for j=1:1:length(v)
npv(i,j)=1-(1.07882*Cs)*(60*delta_p(i))/((u*n)*v(j));
npt(i,j)=1/(1+(Cv*(u*n)/(60*delta_p(i)+Cf+(2*3.141592*Tc)/(delta_p(i)*Vmax))*(1.07882/v(j))));
np(i,j)=npv(i,j)*npt(i,j);
nmv(i,j)=1/(1+(Cs*delta_p(i))/(u*(v(j)/(0.377*rd))));
nmt(i,j)=1-Cf-(Cv*(u*(v(j)/(0.377*rd))/delta_p(i)));
nm(i,j)=nmv(i,j)*nmt(i,j);
n(i,j)=np(i,j)*nm(i,j)*nc;
end
end
figure(1)
plot(v,n(1,:),'-b')
hold on
plot(v,n(2,:),':r')
plot(v,n(3,:),'-.b')
plot(v,n(4,:),'--k')
title('泵—马达系统的理论效率分析')
xlabel('钻机行走速度(Km/h)')
ylabel('总效率')
legend('8Mpa','15Mpa','25Mpa','35Mpa')
grid on
hold off
运行时出现
??? Error using ==> mrdivide
Matrix dimensions must agree.
Error in ==> DTH at 14
npv(i,j)=1-(1.07882*Cs)*(60*delta_p(i))/((u*n)*v(j));
还有15,17,18都有同样的问题,如何解决
要计算的公式如下:
v=0.01:0.01:10;
delta_p=[8000000 15000000 25000000 35000000];
Cs=0.8e-9;
u=0.051;
n=2200;
Cv=0.2e+6;
Cf=0.01;
Tc=5;
Vmax=45e-6;
nc=0.9;
rd=0.3223;
for i=1:1:length(delta_p)
for j=1:1:length(v)
npv(i,j)=1-(1.07882*Cs)*(60*delta_p(i))/((u*n)*v(j));
npt(i,j)=1/(1+(Cv*(u*n)/(60*delta_p(i)+Cf+(2*3.141592*Tc)/(delta_p(i)*Vmax))*(1.07882/v(j))));
np(i,j)=npv(i,j)*npt(i,j);
nmv(i,j)=1/(1+(Cs*delta_p(i))/(u*(v(j)/(0.377*rd))));
nmt(i,j)=1-Cf-(Cv*(u*(v(j)/(0.377*rd))/delta_p(i)));
nm(i,j)=nmv(i,j)*nmt(i,j);
n(i,j)=np(i,j)*nm(i,j)*nc;
end
end
figure(1)
plot(v,n(1,:),'-b')
hold on
plot(v,n(2,:),':r')
plot(v,n(3,:),'-.b')
plot(v,n(4,:),'--k')
title('泵—马达系统的理论效率分析')
xlabel('钻机行走速度(Km/h)')
ylabel('总效率')
legend('8Mpa','15Mpa','25Mpa','35Mpa')
grid on
hold off
运行时出现
??? Error using ==> mrdivide
Matrix dimensions must agree.
Error in ==> DTH at 14
npv(i,j)=1-(1.07882*Cs)*(60*delta_p(i))/((u*n)*v(j));
还有15,17,18都有同样的问题,如何解决
要计算的公式如下: