工程流体力学水力计算编程.docxVIP

流体力学管路水力的计算哈工大姓名：王小亮学号：1090230113 班号：0902301一.问题提出为了实现在已知参数（总流量、粘度、管长、管径、粗糙度、总作用水头等）的情况下，能直接算出已知管路系统的基本流动参数（流速、分流量、损失因数、雷诺数、沿裎损失因数等）的目的，为此特别编写了简单管路系统流动参数计算的程序。该程序能实现串联和并联管路系统流动参数的计算。需要指明的是，由于本人编程能力有限，且为了能计算书上例题的管路系统，故第二类问题的串联管路系统有且仅有两个串联管子，其余均为三个管子串联或并联。二.数学模型及算法1.算法首先，将已知参数的实际管路系统抽象简化为理想物理模型，并根据管路系统类型进行分类；然后，对其进行理论分析，计算出流动参数的计算方程；最后，通过编程实现对所求流动参数的计算。2.数学模型（1）串联管路系统的第一类问题已知流过串联管路的流量，介质参数（），管路参数（），求所需要的总水头。如下图：设为入口损失因数，对A、B两截面列伯努力方程有根据连续性方程的又由，由公式可以计算出，从而求出h。（2）串联管路系统的第二类问题已知总水头h，介质参数（），管路参数（），求通过的流量如下图：设为入口损失因数，对A、B两截面列伯努力方程有根据连续性方程的由此可得又，，由公式可以计算出。将算出的与所取得对比，若二者之差均满足所取得精度，则计算结束，否则令作为新的重新计算为止。最终可得流量（3）并联管路系统的第一类问题已知两点间的压力降（即能量损失）h，介质参数（），管路参数（），求总流量如下图：先取const，（i=1，2，3，下同）;由达西公式可求得所以由公式可以计算出，将算出的与所取得对比，若二者之差均满足所取得精度，则计算结束，否则令作为新的重新计算为止。则（4）并联管路系统的第二类问题已知总流量，介质参数（），管路参数（），求各分支管路的流量及能量损失h如下图：根据经验，先取h=const;由此h值根据并联管路第一类问题计算出各分支管路的流量（i=1，2，3，下同）;则蒋总流量按如下分配用计算出的流量，结合公式、，可以计算出，从而求出；若中任两个之差满足给定精度，则h为所求值，否则令h=，从头重新计算，直到满足精度为止。三.程序流程图四.程序变量说明及程序清单1.主要变量说明di（i=1,2,3）：各管道的直径（m）；li（i=1,2,3）：各管道的长度（m）；dlti（i=1,2,3）：各管道的粗糙度（m）nu：流体粘度（m^2/s）;Vi（i=1,2,3）：各管道中的流速（m/s）；Rei（i=1,2,3）：各管道的雷诺数；lmdi（i=1,2,3）：实际各管道的沿裎损失因数；lmdii（i=1,2,3）：假设各管道的沿裎损失因数；qvi（i=1,2,3）：实际各管道中的流量（m^3/s）；qvii（i=1,2,3）：假设各管道中的流量（m^3/s）；2.程序清单# include stdio.h# include math.h# define pi 3.1415926# define g 9.81double CL (double dlt,doubled,double Re){double lmd1,temp,lmd=0,Rec1,Rec2,Rec3; Rec1=2320; Rec2=80*d/dlt; Rec3=4160*pow(d/(2.0*dlt),0.85);if(ReRec1) {lmd=64/Re; }else if(ReRec2) {if(Re1e5) {lmd=0.3164/pow(Re,0.25); }else if(Re3e6) {lmd=0.0032+0.221*pow(Re,-0.237); } }else if(ReRec3) {lmd=0.01;while (fabs(lmd-lmd1)=1e-6) { lmd1=lmd;temp=log(dlt/(3.7*d)+2.51/Re/sqrt(lmd1))/log(10);lmd=pow(1.0/(-2.0)/temp,2); } }else {lmd=pow((1.74+2.0*log(d/(2.0*dlt))/log(10)),-2); }returnlmd;}main(){inti=0; double qv,qvt,nu,rou,h,s,ID,DK1,DK2,DK;double l1,l2,l3,d1,d2,d3,dlt1,dlt2,dlt3;double lmd1,lmd2,lmd3,lmd11=0.025,lmd22=0.03,lmd33=0.02