MATLAB在化工管路设计中的应用学案.ppt

MATLAB在化工管路设计中的应用 化工管路设计 处理数据量大，应用方程复杂而无法直接求解。若用常规手工计算, 往往要对同一算式进行十几次甚至几十次反复计算, 才能得到结果。 目前有效的解决途径之一就是用计算机取代常规手工计算提高计算速度和准确性。 背景： 流体流动是众多化工过程操作单元的基础之一, 其核心就是管路设计型计算。 ｛ 连续性方程 伯努利方程 能量损失计算式 常规手工计算基本是以迭代试差法作为化工管路设计的解决方法。 Matlab计算 ? 综合运用 背景： 问题提出： 图1为一个输水管路示意图, 由液面1至截面2管路全长300m( 包括局部阻力的当量长度) , 输水管中的流速为2.36×10-3m3/s, 截面1至液面2间有一闸门阀, 其间的直管阻力可以忽略。管壁相对粗糙度ε/d=0.004, 水温为20℃, 粘度为1.0 ×10-3mPa·s 在阀门全开时, 摩擦因子采用λ=0.1(ε/d+68/Re)0.23 , 试设计管路直径大小以满足实际生产要求。 1、水流的初速度u0，通过已知条件计算出管子直径d0 2、再根据流体的密度ρ、粘度μ、速度u，求出雷诺数Re，根据经验公式求出摩擦系数λ。 3、通过能量损失公式算出此时的管径d 4、根据流速算出新的速度u，与初始的u0迭代试差，若试差10-4则退出运算，所求的d就是满足实际生产要求的管径。 问题求解： 问题求解： 采用连续性方程式、柏努利方程式与范宁公式(能量损失公式)联合求解。 结论： 运行程序以后可得到以下结果: u=1.073243m/s; d=0.052913m 因此选用与结果相符或相近规格的管路即可。 Matlab 程序： L=300; %已知管路长度, m; h=10; %液位高度差, m; den=1000; %流体密度, kg/m3; vis=1.000e-3; %流体黏度, Pa?s; V=0.00236; %设计管路的输水流量要求, m3/s; e=2.12e-4; %管壁绝对粗糙度, mm; s1=0.17;s2=1; %s1代表闸门阀全开时的 局部阻力系数, s=s1+s2; %s2代表出口突然扩大的阻力系数; g=9.81; %重力加速度, m/s2; hf=h*g; %总阻力损失,J/kg; 由柏努利 方程推导得到。 u0=input(请输入初始流速数值u0=) %初值设定范围可以较大, m/s; i=1; fprintf( 试差计算开始); while(1) d=2*sqrt(V/(pi*u0)); Re=d*den*u0/vis; if Re = 2000 Lamda=64/Re; %层流时摩擦系数求解公式; else Lamda=0.1* (e/d+68/Re)^0.23; % 湍流时摩擦系数求解公式; end u=sqrt(2*hf/(Lamda*L/d+s)); fprintf(i=%d\r u=%fm/s\n,i,u); if (abs (u-u0)1.0e-4) %迭代试差结果要求精度为10-4; break; else u0=u; i=i+1; end end d=2*sqrt(V/(pi*u)); fprintf(d=%fm,d);