哈工大 流体力学 大作业 管路计算..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.81 float 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);} NN开 始调用CL求λ打印输出结束读qV,d,l,ν,△ N N 开 始 调用CL求λ 打印输出 结束 读qV,d,l,ν,△ 图1-1 简单管路程序框图(1) 1.2简单管路流量的计算 分析：此题为第二类问题，要用试算法，其步骤如下： 1）取=0.03 2）由达西公式可求得 3） 4）调用CL子程序计算。 5）检查和之差是否小于给定精度，如果不满足，将算出的作为重复2）3）4）步骤，直至满足后计算下一步。 6） 1.2.1 程序框图 见图1-2。N N 开 始 打印输出 结束 读hf,d,l,ν,△ λ=0.03 调用CL求λ |λ－λ1|ε 图1-2 简单管路程序框图(2) 1.2.2 程序说明： 在程序中，参数表示同前述程序相同，取精度=0.000005。 1.2.3 程序代码： #includemath.h #includestdio.h #define pi 3.1415926 #define g 9.81 float 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.03 2）由代入达西公式可求得 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.81 float 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,%