资源描述：
《中国石油大学(北京)瞬变流编程实践》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、14中国石油大学（北京）输油管道瞬变流实验报告实验课程：实验名称：班级：学号：姓名：实验台号：成绩：实验日期：年月日14题目：条件：一管道输送柴油。管道由首站和末站组成。管径：457×7.1mm站间距为130km。泵机组串联。密度842.1,e=0.01mm,粘度7.02×10-6m2/s,K=1.39×109Pa，E=207×109Pa，Kv=1.75,3.83,8.38,18.32,40.1,87.62,191.6,419.0,916.2,2003.4,（∞）1#泵转速2970RPM2#泵转速2970RPM
2、3#泵转速2970RPM流量m3/h扬程m功率kw流量m3/h扬程m功率kw流量m3/h扬程m功率kw496.66148.46241.17411.26280.39387.06350.4306.94385.32575.68146.94284.85634.45255.03458.6602.35295.16511.53855.39142.02325.77853.39232.04539.42855.11273.9644.49边界条件:首站给油泵定值30m;末站进压0.4MPa;事故工况，末站关阀;阀门从全开到全关的时间
3、为：10秒。①假设：1、正常情况下首站给油泵定值30m，油品进入各主泵及其附件产生的局部磨阻为30m；2、泵特性曲线中流量扬程曲线符合二次方程规律；3、管道处于水平，即纵断面图为直线；4、各阀之间用固定墩支撑，且管材泊松系数为0.3；5、阀门开度与时间为一次函数关系，满足阀门打开t时间后，阀门开度为(1-0.01t);6、下游边界为阀门和固定水位罐，其在固定水位罐管道末端阀门产生一恒定压力0.4MPa。简化流程图如下：②利用MATLAB软件编程如下：ticclc;clear14%求解各泵的特性曲线p1=[496
4、.66148.46575.68146.94855.39142.02];p2=[411.26280.39634.45255.03853.39232.04];p3=[350.4306.94602.35295.16855.11273.9];%用二次方程求解各泵的A、B值r1=polyfit(p1(:,1)./3600,p1(:,2),2);r2=polyfit(p2(:,1)./3600,p2(:,2),2);r3=polyfit(p3(:,1)./3600,p3(:,2),2);%常数项G=9.81;%重力加速度r
5、ou=842.1;%密度od=0.457;%外径e=0.01*0.001;%绝对粗糙度sigma=7.1*0.001;%壁厚id=od-2*sigma;%内径mu=7.02*10^(-6);%粘度mew=1.39*10^9;%水的弹性系数modulusofelasticitymei=207*10^9;%钢管弹性模量dbs=;%站间距kesi=e/id;%相对粗糙度fi=1-0.3^2;%修正系数FSP=0.4*10^6;%末站进压HFSP=FSP/(G*rou);%末站压头T=100;%计算时间%比例系数f=0
6、.0246*mu^0.25/id^4.75;%列冰宗公式，单位长度摩阻系数a1=(r1(1)+r2(1)+r3(1));%二次方系数a2=(r1(2)+r2(2)+r3(2));%一次方系数a3=(r1(3)+r2(3)+r3(3));%常数项pf=8*1.75/(pi^2*id^4*G);%末站阀摩阻%二分法求根pl=0;%左边界pr=1;%右边界delta=10^(-9);%误差maxl=1+round((log(pr-pl)-log(delta))/log(2));%最大计算步数fork=1:maxlmi
7、d=(pl+pr)/2;%求中间点ypl=-f*dbs*pl^1.75-HFSP-pf+a1*pl^2+a2*pl+a3;%左边界值ypr=-f*dbs*pr^1.75-HFSP-pf+a1*pr^2+a2*pr+a3;%右边界值ymid=-f*dbs*mid^1.75-HFSP-pf+a1*mid^2+a2*mid+a3;%中间点边界值ifymid==0pl=mid;pr=mid;%已解出根14elseifypl*ymid>0pl=mid;%左边界去除，求解区域减半elseypr*ymid>0;pr=mid;
8、%右边界去除，求解区域减半endif(pr-pl)9、
10、Re>Re1disp('流动不处于水力光滑区');e
