城给水管网水力计算程序及例题
给水排水管道工程
课程设计指导书
环境科学与工程学院
第一部分城市给水管网水力计算程序及习题一、程序
#define M 18
#define N 6
#define ep 0.01
#include
int sgn(doublex);
main()
{ int k, i,ko,q,p,flag=0;
double h[M];
double
l[]={?};
double
D[]={?};
double
Q[]={?};
int io[]={?};
int jo[]={?};
doublef[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(doublex)
{ if(x>0)return 1;
elseif(x==0) return 0; elsereturn -1;
}
变量说明
ep --- 环内水头损失闭合差允许值(m);
N ——环数；
M ——管段数；cz ——管道粗糙系数；k ——管段编号；kO ――校正次数；
L ——管段长度(m)；
D ---- 管径(m);
Q ――管段流量(L/s)；
io ――管段所属环号(小环号)；初分流量为顺时针，io 为正，初分流量为逆时针，io 为负；
jo ――管段所属环号(大环号)，均为负，不相临为零；
s ――管段摩阻系数；h 管段水头损失(m)。
1基础资料
(1)城市总体规划概况:
某市近期规划人口为12万，用水普及率预计100%，城区大部分建筑在6 层，屋内有给排水卫生设备和淋浴设备，市内有工业企业甲。
(2)城市用水情况：城市生活用水量变化情况如下表：
(3)工业企业基本情况
甲企业职工人数为1200人，分三班制(0、8 16时)，每班8小时，每班
400人，无高温车间，每班淋浴人数为250人；生产用水量为3000立方米/ 日，均匀使用，工业用水要求水压不小于24米，水质同生活饮用水：工厂房屋最大体积为5000立方米(厂房)，房屋耐火等级为三，生产品危险等级为乙。
2、要求
进行该市给水管网和流量调节构筑物的设计计算，具体包括：
(1)计算最咼日用水量、最咼日最咼时用水量；