编制程序流图如下:
程序如下:
em=1;tc=6;H0=92.12;a=1100;D=0.9;Pair=101325;g=9.806;f=0.02;
dens_water=1000;L=1300;n=10;Cd0Ag0=0.01817; tc=6;Tmax=10;
A=pi*(D^2)/4;B=a/(g*A);
delta_t=L/(a*n);delta_x=L/n;
R=f*delta_x/(2*g*D*A^2);
t=0;
Q0=Cd0Ag0*sqrt(2*g*H0);
V0=Q0/A;
Hp1=H0+Pair/(dens_water*g)-V0^2/(2*g);
Qp=Q0+zeros(1,n+1);Hp=Hp1+zeros(1,n+1);
Hns=zeros(1,85);Qns=zeros(1,85);
Hns (1)=Hp(n+1);Qns(1)=Qp(n+1);
for k=1:84
t=k*delta_t;
if t<=Tmax
Hpp=zeros(1,n+1);Qpp=zeros(1,n+1);
for i=1:n+1
if i<2
Cm=Hp(i+1)-B*Qp(i+1)+R*Qp(i+1)*abs(Qp(i+1));
Hpp(i)=Hp(i);Qpp(i)=(Hpp(i)-Cm)/B;
通过 k, 增加时间 delt_t
计算管道各节点初值
t<Tmax
各节点值更新
Y
阀门节点初值保存
赋值时间
t=0
开 始
输入数据
计算常量 B,R , delt_t,delt_x,
计算管道入口节点、内节点、阀门节点的值
结束
N
置存阀门节点值
评论6
最新资源