Question

matlab解方程的问题 为什么算不出这个d0???

a=0.7;%小孔节流系数
%节流小孔直径mm
PS=4e6;%供油压力pa
PR=PS/1.7;%油膜压力
density=0.9e3; %油的密度(kg/m3)
h0=0.06e-3:0.005e-3:0.105e-3; %初始估算轴承间隙,单位：m
D=150e-3; R=D/2; %轴承直径、半径，单位：m
c1=15e-3; %轴向封油边
l=120e-3;
b1=25e-3; %周向封油边
theta=68.5998*pi/180/2; %2*theta是包角
u=45e-3; %粘度pa s

h0length=length(h0);
table=zeros(1,h0length);

for k=1:h0length
h0=0.06e-3+k*0.005e-3
% Q0=0.25*a*pi*d0^2*sqrt( 2*(PS-PR)/density);
QR=R*h0^3/(6*u*c1)*PR*(l*c1/R*b1+2*theta)
d0=solve('0.25*a*pi*d0^2*sqrt( 2*(PS-PR)/density)=QR','d0')
d0=vpa(d0,4)
table(1,k)=d0;
end

Answer 1

a=0.7;%小孔节流系数
D=3;%节流小孔直径mm
PS=4e6;%供油压力pa
PR=PS/1.7;%油膜压力
density=0.9e3;%油的密度(kg/m3)
h0=0.06e-3:0.005e-3:0.105e-3; %初始估算轴承间隙,单位：
%mD=150e-3;
R=D/2; %轴承直径、半径，单位：
c1=15e-3; %轴向封油边
l=120e-3;
b1=25e-3; %周向封油边
theta=68.5998*pi/180/2; %2*theta是包角
u=45e-3; %粘度pa s
h0length=length(h0);
table=zeros(1,h0length);
for k=1:h0length
h0=0.06e-3+k*0.005e-3 ;% Q0=0.25*a*pi*d0^2*sqrt( 2*(PS-PR)/density);
QR=R*h0^3/(6*u*c1)*PR*(l*c1/R*b1+2*theta);
%d0=solve('0.25*a*pi*d0^2*sqrt(2*(PS-PR)/density)=QR','d0');
d0=(QR/(0.25*a*pi*sqrt(2*(PS-PR)/density)))^0.5;
%d0=vpa(d0,4);
table(1,k)=d0;
end
table

按上面的改后看看，D值可改

追问

这里用solve，怎么都字符运算，能用SOLVE来解决这个方程么？

追答

能用一般语句解决的就不用SOLVE，节约时间，反之却出问题