角联网路通风系统风阻的分析.docxVIP

角联网路通风系统风阻的分析 通风系统的重要参数之一是通风系统的重要参数之一。角网络道路的总风阻对分析和计算矿山通风系统具有重要意义。由于角联网路的复杂性, 其总风阻不能用求串、并联网路总风阻的公式计算。本文以矿井通风基本原理为基础, 通过任意给定角联网路的总风量, 借助于计算机编程运算, 求得了角联网路的总风阻, 且误差很小。 1 总风阻r计算 图1为角联网路的示意图。假定巷道AC、CB、AD、CD、DB的风阻均已知, 且分别为R1、R2、R3、R4、R5。首先任意给定通过角联网路的总风量Q, 通过所编计算机程序运算可求得在已知各风阻及总风量条件下通过巷道AC、CB、AD、CD、DB的风量Q1、Q2、Q3、Q4和Q5, 然后可求得总风压HAB。再利用总风量Q、总风压HAB与总风阻R之间的关系, 求得总风阻R。 设风流方向如图1中箭头所示, 按风压平衡定律, 有以下两方程成立: HAC+HCB=HAD+HDB(1) HAC+HCD=HAD(2) 按照通风阻力定律, 由式 (1) 、式 (2) 可得: R1Q2112+R2Q2222=R3Q2332+R5Q2552(3) R1Q2112+R4Q2442=R3Q2332(4) 由风量平衡定律, 有下面三式成立: Q2=Q1-Q4(5) Q3=Q-Q1(6) Q5=Q3+Q4(7) 为便于求解, 设Q1=x,Q4=y, 由式 (3) ～ (7) 可得以下二元方程组: R1x2+R2(x-y)2=R3(Q-x)2+R5(Q-x+y)2(8) R1x2+R4y2=R3(Q-x)2(9) 由方程组 (8) 和 (9) 整理, 可得方程: Ux4+Vx3+Wx2+Hx+I=0 (10) 其中,U、V、W、H和I均为已知数R1、R2、R3、R4、R5和Q的函数。由于目前不能用一般解方程的方法求出x的解, 必须借助于计算机编程运算才能求出方程 (10) 中x的近似根。 2 wx+h算法 为了建立计算机程序, 必须先用数值计算方法获得方程 (10) 中的x的迭代计算公式。为此, 先将方程 (10) 的左边用函数F(x) 表示, 即: F(x) =Ux4+Vx3+Wx2+Hx+I F(x) 的导数F′ (x) =4Ux3+3Vx2+2Wx+H 按数值计算法中的牛顿迭代法有: xn+1=xn-F(xn) /F′ (xn) (n=0, 1, 2, …) (11) 为了减少迭代次数, 在n=0时, 给定初值x0(x0=Q/2) , 然后按照迭代公式 (11) 逐步求出x1、x2、x3……xn, 直到前后两次x值之差 (即xn+1-xn) 小于0.01时, 所求xn即为最终解x(即Q1) 。然后将Q1值代入式 (9) 、 (5) 、 (6) 和 (7) , 分别求出Q4、Q2、Q3、Q5之值, 最后求出总风阻R。 R= (R1Q2112+R2Q2222) /Q2 为了验证所求各风量及风阻R值的精度是否满足要求, 最后需求出网孔A—C—B—D—A的压差C1和网孔A—C—D—A的压差C2。如果C1和C2值均很小, 则所求结果精度满足要求。计算机程序框图如图2所示。 3 网孔a—程序运行结果及分析 任意给定角联网路风阻R1、R2、R3、R4、R5的4组数值, 且任意给定每组数值对应的角联网路总风量Q的5个数值 (表1) 。在每一总风量Q条件下分别求出角联网路总风阻R, 通过巷道AC和CD的风量Q1和Q4, 网孔A—C—B—D—A和A—C—D—A的压差C1和C2, 风量Q2、Q3、Q5值表中未打印。 以表1中的第一组数值为例对运行结果进行分析:给定角联网路风阻R1、R2、R3、R4和R5及总风量Q值, 由表1可以看出, 在不同风量下求得总风阻R的前6位小数均相同。因此, 它们之间的差小于10-6N·s2/m8, 并且压差C1值均小于10-3Pa,C2值均小于10-5Pa, 即网孔A—C—B—D—A和网孔A—C—D—A的压差基本接近零, 所求结果误差很小。分析表1中第二组、第三组和第四组数值所运行的结果, 不难发现有与上述基本相同的规律。因此, 可以认为, 网路的总风阻R与给定总风量Q的大小无关, 这说明了风阻R反映巷道本身的特性。此外, 给定总风量Q较小时, 所求C1和C2值就越小。因此, 任意给定总风量Q一个较小的值后, 求得的总风阻R即为真值。 4 迭代算法的收敛 利用上述方法求角联网路总风阻, 虽然数学方程建立较复杂, 但所用迭代计算方法收敛快, 一般迭代3～4次即可达到表中结果的精度。此外, 程序运算时间短, 输入已知数据即可得到结果。