城给水管网水力计算程序及例题.DOC

给水排水管道工程 课 程 设 计 指 导 书 环境科学与工程学院第一部分 城市给水管网水力计算程序及习题 一、程序 #define M 18 #define N 6 #define ep 0.01 #include math.h int sgn(double x); main() { int k, i,ko,q,p,flag=0; double h[M]; double l[]={?}; double D[]={?}; 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[i]); r[i]=r[i]+(h[k]/Q[k]); } dq[i]=-(f[i]/(r[i]*2)); } { if (fabs(f[N])=ep) 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\n,f[i]);} printf(ep=%f\n,0.01); printf(n=%d,m=%d,ko=%d\n,N,M,ko); for(k=0;k=M-1;k++) { printf(%d),k+1); printf(k=%d, l=%f, h=%f, ,k+1,l[k],h[k]); printf(Q=%f, ,Q[k]*1000); printf(v=%f\n,4*Q[k]/(3.1416*pow(D[k],2))); } } int sgn(double x) { if(x0)return 1; else if(x==0) return 0; else return -1; } 变量说明 ep——环内水头损失闭合差允许值(m)； N ——环数； M ——管段数； cz——管道粗糙系数； k ——管段编号； k0——校正次数； L ——管段长度(m)； D ——管径(m)； Q ——管段流量(L/s)； io——管段所属环号(小环号)；初分流量为顺时针，io为正，初分流量为逆时针，io为负； jo——管段所属环号(大环号)，均为负，不相临为零； s ——管段摩阻系数； h ——管段水头损失(m)。 城市总体规划概况某市近期规划人口为12万，用水普及率预计100％，城区大部分建筑在6层，屋内有给排水卫生设备和淋浴设备，市内有工业企业甲。 城市用水情况城市生活用水量变化情况如下表： 时间 0～1 1～2 2～3 3～4 4～5 5～6 6～7 7～8 8～9 9～10 10～11 11～12 用水量 1.10 0.70 0.90 1.10 1.30 3.91 6.61 5.84 7.04 6.69 7.17 7.31 时间 12～13 13～14 14～15 15～16 16～17 17～18 18～19 19～20 20～21