哈工大 流体力学 大作业 管路计算.docxVIP

Harbin Institute of Technology流体力学大作业院系：班级：姓名：学号：指导教师：?哈尔滨工业大学计算机在管路水力计算中的应用等径无分支管路的计算1.1简单管路沿程损失的计算分析：此题为第一类问题，可用如下步骤求解。1) 2)3) 调用CL子程序计算。4) 1.1.1 程序框图见图1-1。1.1.2 程序说明:在程序中，hf表示，q表示，pi表示。1.1.3 程序代码：#includemath.h#includestdio.h#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;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）取=0.032）由达西公式可求得3）4）调用CL子程序计算。5）检查和之差是否小于给定精度，如果不满足，将算出的作为重复2）3）4）步骤，直至满足后计算下一步。6）1.2.1 程序框图见图1-2。图1-2 简单管路程序框图(2)1.2.2 程序说明：在程序中，参数表示同前述程序相同，取精度=0.000005。1.2.3 程序代码：#includemath.h#includestdio.h#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）取=0.032）由代入达西公式可求得3）4）调用CL子程序计算。5）检查和之差是否小于给定精度，如果不满足，将算出的作为重复2）3）4）步骤，直至满足后计算下一步。1.3.1 程序框图见图1-3。1.3.2 程序说明：在程序中，参数表示同前述程序相同，取精度=0.000005。3.3.3 程序代码：#includemath.h#includestdio.h#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,lll;printf(please input:dlt,q,l,hf,nu\n);scanf(%f,%f,%f,%f,%f,dlt,q,l,hf,nu);lmd1=0.03;d=pow(8*l*q*q*lmd1/(pi*pi*g*hf),0.2);v=4*q/(pi*d*d);re=v*d/nu;lmd=CL(dlt,d,re,id);while(fabs(lmd-lmd1)=1e-6) {lmd1=lmd; lll=8*l*q*q*lmd1/(pi*pi*g*hf); d=pow(lll,0.2); v=4*q/(pi*d*d); re=v*d/nu; lmd=CL(dlt,d,re,id);}printf(v=%7.4f m/s,re=%7.4f,lmd=%7.4f,d=%7.4f m ,v,re,lmd,d);}图1-3 简单管路程序框图(3)串联管路的计算方法2.1串联管路总水