资源描述：
《哈工大-流体力学-大作业-管路计算.docx》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、HarbinInstituteofTechnology流体力学大作业院系：班级：姓名：学号：指导教师：©哈尔滨工业大学计算机在管路水力计算中的应用1等径无分支管路的计算1.1简单管路沿程损失的计算分析：此题为第一类问题，可用如下步骤求解。1)2)3)调用CL子程序计算。4)1.1.1程序框图见图1-1。1.1.2程序说明:在程序中，hf表示，q表示，pi表示。1.1.3程序代码：#include#include#definepi3.#defineg9.81floatCL(float,float,float,int);voidmain(
2、){intid;id=0;floatdlt,d,l,q,nu,hf,v,re,lmd;printf("pleaseinput:dlt,d,l,q,nu");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);}NN开始调用CL求λ打印输出结束读qV,d,l,ν,△图1-1简单管
3、路程序框图(1)1.2简单管路流量的计算分析：此题为第二类问题，要用试算法，其步骤如下：1）取=0.032）由达西公式可求得3）4）调用CL子程序计算。5）检查和之差是否小于给定精度，如果不满足，将算出的作为重复2）3）4）步骤，直至满足后计算下一步。6）1.2.1程序框图见图1-2。N开始打印输出结束读hf,d,l,ν,△λ=0.03调用CL求λ
4、λ－λ1
5、<ε图1-2简单管路程序框图(2)1.2.2程序说明：在程序中，参数表示同前述程序相同，取精度=0.。1.2.3程序代码：#include#include#definepi3.#
6、defineg9.81floatCL(float,float,float,int);voidmain(){intid;id=0;floatdlt,d,l,q,nu,hf,v,re,lmd,lmd1;printf("pleaseinput:dlt,d,l,hf,nu");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=sq
7、rt(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.4fm/s,re=%7.4f,lmd=%7.4f,q=%7.4fm3/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程序说明：在程序中，参数表示同前述
8、程序相同，取精度=0.。3.3.3程序代码：#include#include#definepi3.#defineg9.81floatCL(float,float,float,int);voidmain(){intid;id=0;floatdlt,d,l,q,nu,hf,v,re,lmd,lmd1,lll;printf("pleaseinput:dlt,q,l,hf,nu");scanf("%f,%f,%f,%f,%f",&dlt,&q,&l,&hf,&nu);lmd1=0.03;d=pow(8*l*q*q*lmd1/(pi*pi*
9、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.4fm/s,re=%7.4f,lmd=%7.4f,d=%7.4fm",v,re,lmd,d);}YNN开始打印输出结束读hf,qV,l,ν,△λ=0.03调用CL
