哈工大流体力学作业.docVIP

Harbin Institute of Technology 流体力学大作业 院 系： 班 级： 设 计 者： 学 号： 指导教师： 设计时间： 并联管路沿程损失计算 分析：此题为关联管路的第二类问题。由于只知道总流量，各分支流量及能量损失均不知道，因此计算较复杂，其具体步骤如下： 1）根据经验，先给定一个hf值 。 由此hf值求出各分支管的流量，，，（其步骤参见简单管路第二类问题）。 计算其和=++ 将总流量按如下比例分配 =./，=./，=./ 用计算出的、、去求以核对流量分配的正确性, 若求得各分支管路的损失h1、h2、h3中任意两者之间的差小于某一给定精度，则认为hf正确，否则进行下一步。 令hf= （h1+h2+h3）/3，重复（b）（c）（d）（e）步骤，直至满足所给定的精度为止。 1 程序框图 见下图 并联管路程序框图 2 程序清单： #include windows.h #includemath.h #includestdio.h #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;i3;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;i3;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;i3;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,6