给水排水管道工程课程设计指导书.docVIP

给水排水管道工程 课 程 设 计 指 导 书 环境科学与工程学院 第一部分 城市给水管网水力计算程序及例题 1 程序框图 开始 输入原始数据 计算管段摩阻 计算各环水头损失闭合差及校正流量 计算管段流量 闭合差 EP 打印输出原始数据及计算结果 结束 程序 #define M 22 #define N 8 #define ep 0.01 #includestdio.h #includestring.h #include math.h int sgn(double x); main() { FILE *fp; int k, i,ko,q,p,flag=0; double h[M]; double l[]={750,885,715,850,610,610,610,610,610,750,885,715, 850,800,680,680,680,680,290,885,715,850}; double D[]={0.45,0.45,0.3,0.25,0.5,0.25,0.25,0.25,0.4,0.7,0.6,0.45,0.2,0.5,0.25,0.25,0.15,0.2,0.45,0.45,0.45,0.35}; double Q[]={………}; int io[]={………}; int jo[]={………}; double f[N+1],r[N+1],dq[N+1]; for(k=0;k=M-1;k++) { Q[k]=Q[k]*0.001; } for(k=0;k=M-1;k++) { Q[k]=Q[k]*sgn(io[k]); } ko=0; loop: for(k=0;k=M-1;k++) { h[k]=10.67*pow(fabs(Q[k]),1.852)*l[k]; h[k]=h[k]/(pow(100,1.852)*pow(D[k],4.87))*sgn(Q[k]); } for(i=1;i=N;i++) { f[i]=0; r[i]=0; dq[i]=0; for(k=0;k=M-1;k++) { if(abs(io[k])!=i) goto map; f[i]=f[i]+h[k]; r[i]=r[i]+(h[k]/Q[k]); map: if( abs(jo[k])!=i) continue; f[i]=f[i]+h[k]*sgn(jo[k]); r[i]=r[i]+(h[k]/Q[k]); } dq[i]=-(f[i]/(r[i]*2)); } for(i=1;i=N;i++) { if(fabs(f[i])0.01) break; } if(i==N+1) flag=1; if(flag==1) goto like; for(k=0;k=M-1;k++) { p=abs(io[k]); q=abs(jo[k]); Q[k]=Q[k]+dq[p]+(dq[q]*sgn(jo[k])); } ko=ko+1; if(flag==0) goto loop; like: printf(\n\n); for(i=1;i=N;i++) { printf(%f,%f\n,f[i],dq[i]); } printf(ep=%f\n,0.01); printf(n=%d,m=%d,ko=%d\n,N,M,ko); for(k=0;kM;k++) { printf(%d),k+1); printf(k=%d, l=%f, h=%f,,k+1,l[k],h[k]/(io[k]/fabs(io[k]))); pri