标准节流装置的设计与计算

课程设计报告(2013—2014年度第一学期)课程：过程参数检测及仪表题目：标准节流装置的设计与计算院系：自动化系班级：学号：学生姓名：Acceler指导教师：田沛设计周数：一周成绩：日期：2014年1月15日一、课程设计目的与要求本课程设计为检测技术与仪器、自动化专业《过程参数检测及仪表》专业课的综合实践环节。通过本课程设计，使学生加深过程参数检测基本概念的理解，掌握仪表的基本设计方法和设计步骤。二、设计正文第一类命题：已知条件：流体为水，工作压力MPap7.14，工作温度215t℃；管道mmD23320，材料为20号钢新无缝钢管；节流件为法兰取压标准孔板，材料为1Cr18Ni9Ti；mmd34.11720；差压kPap91.156，求（1）给定差压值p下的水流量mq；（2）计算测量时的压力损失。解：（1）辅助计算：查表得到水和水蒸气密度1=856.85kg/3m，动力粘度=127.3610Pa·s，管道线膨胀系数D=12.78610/℃，节流件线膨胀系数d=17.2610/℃，可膨胀性系数=1。mmtDDDt58.233)]20(1[20mmtddtd73.117)]20(1[20（2）查表可知，新无缝钢管的绝对粗糙度K=0.05~0.1mm，(410K/D)max=4.294.9，所以直管段粗糙度符合要求。（3）迭代计算水流量mq：由Stolz方程，得：令式中0ReD=，此时流出系数初始值为0C=0.60274。取精密度判据6101z，利用Matlab进行迭代计算，程序代码如下：A=7912885.84;yita=127.3e-6;b=0.504;Dt=233.58;c0=0.5959+0.0312*b^2.1-0.184*b^8+2.286*b^4/Dt/(1-b^4)-0.856*b^3/Dt;c=c0;z=1;%初值预设5040.0ttDd84.79128851004.0354.0412ttDpdA)62.58,:(856.0)1(286.2)Re10(0029.0184.00312.05959.02034475.065.281.2mmDmmDDDCDwhilez1e-6%迭代计算c1=c;Re=A*c;c=c0+0.0029*b^2.5*(1e6/Re)^0.75;z=abs((c-c1)/c);endcqm=A*Dt*yita*c/0.354迭代2次后，8104.5z，精密度符合要求。此时C=0.602902，mq=400721.5kg/h。对于标准孔板，其压力损失的计算式为pp)1(9.1，将与p值代入，得p=114.23kPa。答：给定差压值p下的水流量mq约为400721.5kg/h；测量时的压力损失约为114.23kPa。第二类命题：已知条件：流体为锅炉给水，最大流量hkgqm/102303max，常用流量hkgqmch/101903，最小流量hkgqm/101003min；工作压力MPap15，工作温度215t℃；管道mmD19920，材料为20号钢新无缝钢管，节流件为角接取压标准孔板，材料为1Cr18Ni9Ti；允许压力损失kPay40，上游第一阻力件为全孔球阀，上游第二阻力件为同一平面上两个90°弯头。（1）选定差压计；（2）计算20dC、、、；（3）确定最小直管段长度210lll、、；（4）计算压力损失并进行核算；（5）计算基本误差。本题设计任务书如下：序号项目符号单位数值1被测流体名称锅炉给水2流量测量范围：常用最大最小qmchqmaxmqminmkg/hkg/hkg/h1900002300001000003工作温度t°C2154工作压力pMPa155管道内径（20°C实测值）D20mm1996管道材料20号钢新无缝管7正常流量下允许的压力损失ypkPa401、准备计算根据锅炉给水工作状态下的压力1p=15MPa和温度t=215°C，查表得流体密度ρ1=857.17kg/m3,动力粘度η=127.4106Pa·s，管道线膨胀系数D=12.78610/℃，节流件线膨胀系数d=17.2610/℃。可膨胀系数=1。计算Dt：mmtDDDt496.199)]20(1[20计算正常流量和最小流量下雷诺数：ReDch=tD4mchq=3600104.12710199.496π190000463-=2.644×106ReminD=tminD4mq=3600104.12710496.991π100000463-=1.3921062、选定差压计由于采用角接取压标准孔板，最大差压的计算为：ΔPmax≤（2~2.5）ypΔPmax=2yp=2×40=80kPaΔpch=Δpmax（qmch/qmmax）2=80×103×（331023010190）2=54.594kPa流量指示仪表的流量刻度上限是有规定的，即为1，1.2，1.6，2，2.5，3，4，5，6，8×10n其中n为任意整数。题中最大流量为230t/h，流量计流量刻度上限为mmaxq=250t/h。故选用1151HP型差压变送器，其量程范围为0～186.8kPa，耐静压为31.5MPa，仪表基本误差为±0.25%。变送器差压量程调整在0～80kPa。3.计算，C，，20d流量公式：qmch=41Cε4β2Dt2chp12，其中qmch、4、Dt2chp12已知，所以不变量A=chtchpDq12m24得迭代方程为：A=421C，令X=421=CA。精密度判据为：nE=AXCA1106。β=25.0211X，dt=βDt=Dt25.0221XXqmch=190×103kg／h；Dt=199.496mm；=857.17kg/m3；Δpch=54.594kPa计算A值：A=chtchpDq242m=32310594.5417.8572)10496.199(43600/190000=0.1745经以上初步计算后，用Matlab进行迭代计算，程序代码如下：A=0.1745;yita=127.3e-6;Re=2.644e6;i=2;Dt=199.496;xn(1:30)=0;b(1:30)=0;c(1:30)=0;dta(1:30)=0;en(1:30)=1;c(1)=0.5959;%初值预设while(abs(en(i-1))1e-6)&&(i4)%开始迭代计算xn(i)=A/c(i-1);b(i)=(xn(i)*xn(i)/(1+xn(i)*xn(i)))^0.25;c(i)=0.5959+0.0312*b(i)^2.1-0.184*b(i)^8+0.0029*b(i)^2.5*(1e6/Re)^0.75;dta(i)=A-xn(i)*c(i);en(i)=dta(i)/A;i=i+1;endwhileabs(en(i-1))1e-6%快速弦截法xn(i)=xn(i-1)-dta(i-1)*(xn(i-1)-xn(i-2))/(dta(i-1)-dta(i-2));b(i)=(xn(i)*xn(i)/(1+xn(i)*xn(i)))^0.25;c(i)=0.5959+0.0312*b(i)^2.1-0.184*b(i)^8+0.0029*b(i)^2.5*(1e6/Re)^0.75;dta(i)=A-xn(i)*c(i);en(i)=dta(i)/A;i=i+1;enddt=b(i-1)*Dt;得：β=0.5271；C=0.6032；td=105.1605mm。又由)]20(1[20tddtd得，20d=104.809mm。综上，β=0.5271；C=0.6032；=1；20d=104.809mm。4.确定最小直管段长度和对管道粗糙度的要求管道布置如下图所示：上游第一阻力件为全孔球阀，上游第二阻力件为同一平面上两个90°弯头。最小直管段L由三段长度（L1，L2，L0）组成，L1是节流件与上游第一阻流件之间的直管段最小长度，由第一阻流件形式和β值决定。L0是上游第一阻流件与第二阻流件之间的直管段，由第二阻流件和β=0.7决定，按表中给出值的一半计算。L2是节流件到下游阻流件直管段的最小长度，与下游节流件的形式无关，只与β值有关。根据β约等于0.55，查表并计算可得210、、lll：0l=1/23620D=3582mm1l=1420D=2786mm2l=620D=1194mm标准孔板对管道相对粗糙度的限值可查表得到，管道相对粗糙度上限值为4.9×104，所以要求节流件上游侧10D以内管道的绝对粗糙度Ks为Ks199.496×4.9×104=0.0978mm新无缝钢管的绝对粗糙度Ks=0.05~0.1mm，符合粗糙度要求。5.计算压力损失并进行核算chp=(1-9.1）Δpch=（1-0.52719.1）×54.594=38.42kPa＜40kPa即chp＜yp，压力损失验算合格。6.计算基本误差及测量结果的不确定度求取流量值相对不确定度的计算公式为：21212224224422])(41)(41)()12()()12()()[(1PdDCqPdDcmqm因为β=0.5271，所以有Ccδ=0.6%；锅炉给水为不可压缩流体，则ε=1，0；DD与dd是与量具的误差，λD、λd值误差和工作温度t值有关。在节流件工作在设计工作温度时，主要取决于量具的误差。按置信概率95%估算，dd=±0.07％，DD=±0.4％；11=±0.06％；由于的选择差压变送器为0.25级，选择显示仪表为1级，则pp=594.54125.0803222%=1.01％，综上：%643.0]%06.04101.14107.0)5271.012(4.0)5271.015271.02(06.0[2222422442mqmq7.结论（1）差压计选用1151HP电容式差压变送器；（2）β=0.5271，C=0.6032，ε=1，20d=104.809mm；（3）最小直管段l0=3582mm，l1=2786mm，l2=1194mm；（4）压力损失chp=38.42kPa＜40kPa，压力损失验算合格；（5）不确定度mqqm=0.643%。三、课程设计总结通过此次课程设计，我得以加深对过程参数检测基本概念的理解，掌握了仪表的基本设计方法和设计步骤，学会了第一类命题的解决方法和第二类命题的设计方法，知道了如何选择节流件形式和差压计以及确定量程，理解了迭代方法的原理和计算方法。本次课程设计也使我分析问题和解决问题的能力得到了提高。四、参考文献[1]苏杰，常太华.过程参数检测及仪表.中国电力出版社,2009.6[2]周建兴，岂兴明，矫津毅，张延伟.MATLAB从入门到精通.人民邮电出版社，2012.6