Harbin Institute of Technology流体力学大作业院 系: 班 级: 姓 名: 学 号: 指导教师: ©哈尔滨工业大学计算机在管路水力计算中的应用1 等径无分支管路的计算1.1 简单管路沿程损失的计算分析:此题为第一类问题,可用如下环节求解。1) v= 4qVπd22) Re=vdν3) 调用 CL 子程序计算λ 。4) hf=λ ldv22g1.1.1 程序框图见图 1-1。1.1.2 程序阐明: 在程序中,hf 表达hf ,q 表达qV ,pi 表达π 。1.1.3 程序代码:#include#include#define pi 3.1415926#define g 9.81float CL(float,float,float,int);void main() NN开 始调用 CL 求 λ打印输出结束读 qV,d,l,ν,△{int id;id=0;float dlt,d,l,q,nu,hf,v,re,lmd;printf("please input:dlt,d,l,q,nu\n");scanf("%f,%f,%f,%f,%f",&dlt,&d,&l,&q,&nu);v=4*q/(pi*d*d);re=v*d/nu;lmd=CL(dlt,d,re,id);hf=lmd*l*v*v/(d*2.0*g);printf("v=%7.4fm/s,re=%7.4f,lmd=%7.4f,hf=%7.4fm",v,re,lmd,hf);} 图 1-1 简单管旅程序框图(1)1.2 简单管路流量的计算分析:此题为第二类问题,要用试算法,其环节如下:1)取λ1=0.032)由达西公式可求得v=√2 gshf√λ1lN开 始打印输出结束读 hf,d,l,ν,△λ=0.03调用 CL 求 λ|λ - λ1|<ε3)Re=vdν4)调用 CL 子程序计算λ 。5)检查 λ 和λ1之差与否不大于给定精度ε ,假如不满足,将算出的 λ 作为λ1反复2)3)4)环节,直至满足后计算下一步。6)qV= πd24 v1.2.1 程序框图见图 1-2。图 1-2 简单管旅程序框图(2)1.2.2 程序阐明:在程序中,参数表达同前述程序相似,取精度ε =0.000005。1.2.3 程序代码:#include#include#define pi 3.1415926#define g 9.81float CL(float,float,float,int);void main(){int id; id=0;float dlt,d,l,q,nu,hf,v,re,lmd,lmd1;printf("please input:dlt,d,l,hf,nu\n");scanf("%f,%f,%f,%f,%f",&dlt,&d,&l,&hf,&nu);lmd1=0.03;v=sqrt(2*g*d*hf/(l*lmd1));re=v*d/nu;lmd=CL(dlt,d,re,id);while(fabs(lmd-lmd1)>=1e-6) {lmd1=lmd; v=sqrt(2*g*d*hf/(l*lmd1)); re=v*d/nu; lmd=CL(dlt,d,re,id);}q=pi*d*d*v/4;printf("v=%7.4f m/s,re=%7.4f,lmd=%7.4f,q=%7.4f m3/s",v,re,lmd,q);}1.3 简单管路管径的计算分析:此题为第二类问题,要用试算法,其环节如下:1)取λ1=0.032)由v= 4qVπd2 代入...