daiweisc 发表于 2008-3-6 21:43

大家帮我看看这个程序有什么问题

请大家帮我看看这个程序有什么问题？见附件
谢谢

[ 本帖最后由 eight 于 2008-3-6 22:16 编辑 ]

eight 发表于 2008-3-6 22:17

原帖由 daiweisc 于 2008-3-6 21:43 发表 http://www.chinavib.com/forum/images/common/back.gif
请大家帮我看看这个程序有什么问题？见附件
谢谢 最好不要让大伙通过下载运行你的程序才知道问题所在

sigma665 发表于 2008-3-6 22:30

function swirlnozzle
% xFsolve.m
% 喷嘴计算，采用《简式旋流喷嘴的设计》文章中的方法，作者“姚增权．马进”

clear all
clc


Q = 0.6689/3600; % 流量，单位NM3/H
alpha = 60; % alpha单位为度
m = 1; % 旋流腔半径与排放口直径之比
a = 2*pi*alpha/360; % a将alpha单位转为弧度
Pl = 1e5; % 压力为100kPa
density = 1000; % 水的密度
miu = 1.0034e-6; % 水的运动粘度
x0 = 0.5;
theta = fsolve(@NonlinEqs,x0,[],a);
Cd = theta*sqrt(theta/(2-theta)); % 流量系数
Dout = sqrt(4*Q/(pi*Cd*sqrt(2*Pl/density))); % 排放口直径
R = m*Dout; % 旋流腔半径
K = sqrt(2)*(1-theta)/(theta*sqrt(theta));
Dp = sqrt(2*R*Dout/K);
Re = 4*Q/(pi*miu*Dp); % 雷诺数
s = 10^(25.8/(log10(Re))^2.58-2); % 摩擦系数
Ky = R*(Dout/2)/((Dp/2)^2+s/2*R*(R-Dout/2)); % 几何常数
theta1 = fsolve(@NonlinEqs1,x0,[],Ky);
Cd1 = theta1*sqrt(theta1/(2-theta1));
Dout1 = sqrt(4*Q/(pi*Cd1*sqrt(2*Pl/density)));
R1 = m*Dout1;
K1 = sqrt(2)*(1-theta1)/(theta1*sqrt(theta1));
Dp1 = sqrt(2*R1*Dout1/K1);
Dp2 = Dp1/sqrt(0.9); % 修正后的入口直径
disp(' 单位：mm')
fprintf('排放口直径=%.1f,修正后的入口直径=%.1f',Dout1*1000,Dp2*1000)


% ------------------------------------------------------------------
function f = NonlinEqs(x,a)
f(1) = tan(a/2)-(2*sqrt(2)*(1-x)/(sqrt(x)*(1+sqrt(1-x))));
% -------------------------------------------------------------------
function f = NonlinEqs1(x,Ky)
f(2) = Ky-sqrt(2)*(1-x)/(x^1.5)

程序见上

dexy 发表于 2008-3-7 09:13

不是能运行出来吗

eight 发表于 2008-3-7 10:55

这种提问方式论坛不太提倡，啥都不说就摆个程序上来，以为大家是神？除非有同行路过，否则，祝楼主好运
页: [1]
查看完整版本: 大家帮我看看这个程序有什么问题