哈工大流体力学大作业.doc

Harbin Institute of Technology 流体力学大作业 院 系： 班 级： 设 计 者： 学 号： 指导教师： 设计时间： 并联管路沿程损失计算 分析：此题为关联管路的第二类问题。由于只知道总流量，各分支流量及能量损失均不知道，因此计算较复杂，其具体步骤如下： 1）根据经验，先给定一个hf值 。 由此hf值求出各分支管的流量，，，（其步骤参见简单管路第二类问题）。 计算其和 ++ 将总流量按如下比例分配 ./， ./， ./ 用计算出的、、去求以核对流量分配的正确性, 若求得各分支管路的损失h1、h2、h3中任意两者之间的差小于某一给定精度，则认为hf正确，否则进行下一步。 令hf （h1+h2+h3）/3，重复（b）（c）（d）（e）步骤，直至满足所给定的精度为止。 1 程序框图 见下图 并联管路程序框图 2 程序清单： #include #include #include #define pi 3.1415926 #define g 9.81 float CL float,float,float ; void main int i; float dlt[3],lmd[3],re[3],d[3],l[3],v[3],h[3],q[3]; float q1[3],lmd1[3]; float qt,hf,hff,nu,qq; for i 0;i 3;i++ printf input:d,l,dlt ; scanf %f,%f,%f,d[i],l[i],dlt[i] ; printf input:qq,nu ; scanf %f,%f,qq,nu ; qt 0; hff 100; loop: for i 0;i 3;i++ hf hff; lmd1[i] 0.03; v[i] sqrt 2*g*d[i]*hf/ l[i]*lmd1[i] ; re[i] v[i]*d[i]/nu; lmd[i] CL dlt[i],d[i],re[i] ; if fabs lmd1[i]-lmd[i] 1e-6 lmd1[i] lmd[i]; v[i] sqrt 2*g*d[i]*hf/ l[i]*lmd1[i] ; re[i] v[i]*d[i]/nu; lmd[i] CL dlt[i],d[i],re[i] ; q1[i] pi*d[i]*d[i]*v[i]/4.0; qt qt+q1[i]; for i 0;i 3;i++ q[i] qq*q1[i]/qt; v[i] 4*q[i]/ pi*d[i]*d[i] ; re[i] v[i]*d[i]/nu; lmd[i] CL dlt[i],d[i],re[i] ; h[i] lmd[i]*l[i]*v[i]*v[i]/ d[i]*2.0*g ; hff h[0]+h[1]+h[2] /3.0; if fabs h[0]-h[1] 0.01||fabs h[1]-h[2] 0.01||fabs h[2]-h[1] 0.01 goto loop; printf q1 %7.4fm3/s,v1 %7.2fm,re1 %7.1f,lmd1 %7.2f,h1 %7.3fm\n,q[0],v[0],re[0],lmd[0],h[0] ; printf q2 %7.4fm3/s,v2 %7.2fm,re2 %7.1f,lmd2 %7.2f,h2 %7.3fm\n,q[1],v[1],re[1],lmd[1],h[1] ; printf q3 %7.4fm3/s,v3 %7.2fm,re3 %7.1f,lmd3 %7.2f,h3 %7.3fm\n,q[2],v[2],re[2],lmd[2],h[2] ; 3 举例 已知：图5-1表示三管路关联，l1 900m，l2 600m，l3 1200m，d1 0.3m，d2 0.2m，d3 0.4m，1 0.0003，2 0.00003，3 0.000024， 1×10-6 m2/s， 0.4m3/s. 求：各分支流量及能量损失hf。 运行程序： Input d 0 ,l 0 ,dlt 0 ,d 1 , l 1 ,dlt 1 ,d 2 ,l 2 ,dlt 2 0.3,900,0.0003,0.2,600,0.00003,0.4,1200,0.000024 Input qq,nu 0.4,1e-6 程序结果如下： q1 0.1056m/s,v1 1.49m/s,re1 448221.4,lmd1 0.02,h1 6.929m q2 0.0552m/s,v2 1.76m/s,re2 351634.8,lmd2 0.01,h2 6.577m q3 0.2392m/s,v3