过长系统过程-2.5化工过程系统分解.pptVIP

分隔只解决了框块——即单元模块组合之间的计算顺序问题，而框块内 部（如ABCDE内部）还包括若干单元和循环回路，这些单元组合的求解 顺序，就要用“断裂”方法解决。 也就是说，分隔问题中的组合节点或回路存在一个“断裂”问题。 （1）断裂的基本概念 先考察方程组的断裂。 f1 ( x2 , x3 ) = 0 f2 ( x2 , x3 , x4 ) = 0 f3 ( x1 , x2 , x3 , x4 ) = 0 f4 ( x1 , x2 , x3 ) = 0 × × ? × f3 × × f4 × ? × f2 × ? ? f1 x3 x1 x4 x2 ? × × × × 2.5.2 断裂（Tearing） 输出变量的指定为： 该方程组已无法分隔，必须联立求解。可以按一般非线性方程组 F（X）= 0 的数值计算方法求解，而通过“断裂”，可以使该四维方程组降解为四个一维方程来求解： 若以 x3 为断裂变量，也称试差变量，则： （a）设 x3 初值 x30 ，由 f1 求 x2：x2 = G1（x30） （b）由 f2 求 x4：x4 = G2（x2，x30） （c）由 f3 求 x1：x1 = G3（x2，x30，x4） （d）由 f4 求 x3 计算值：x3 = G4（x1，x2），亦即 x3 = G（x30） （e）检验 ？若不满足，返回（a）；若满足，则得方程组的解。 在返回（a）时，若以 x30 = x3 产生新的初值，即为直接迭代。 若以x3初值 = qx3原初值 +（1-q）x3 产生，则为 Wegstein 迭代。 这种通过数值迭代将高维方程组降维依次求解的方法，即称为“断裂”。 一般，断裂的方法不只一个。例如 ，这里，流股权代表流股变量数。 该闭环中，S10 有 5 个变量，S11 有 4 个变量，比较少，选择 S11 断裂，会使计算量少易于收敛。也就是说，存在一个最优断裂问题。 S10（5） F G S11（4） （2）最优断裂准则 选择哪些流股为断裂流股？以什么准则来判断断裂流股的选择是否最优？ 迄今为止，已提出的最佳断裂准则基本上有四类： ① 使断裂流股的数目最少（Barkley Motard 1972）； ② 使断裂流股的变量数目最少（Rubin 1962，Sargent Westerberg 1964），即需收 敛的变量数最少； ③ 使断裂流股的总权和最小（Christensen Rudd 1969, Pho,T.K Lapidus 1973）； 给每一流股分配一个“权”，反映断裂该流股时迭代计算的困难程度。“权”是一 个经验数据。一个简单的办法：“权”可以等于每个流股中断裂变量的数目； ④ 使被断裂的回路总数最少（Upadhye Grens 1975）。或者说，找出对直接迭代法 （信息不作加工而反馈回来）具有最佳收敛性能的断裂，即每一回路只断裂一次的 非多余断裂。 上述四类最优断裂准则可用一个统一的数学式表达： s.t． ≥ 1 i=1，2，…，M（回路号） j=1，2，…，N（流股号） ? pj —— 权 有效断裂：{xj┃xj = 1} 0 j 流股未断裂 0 j 流股不属于 i 回路 1 j 流股断裂 1 j 流股属于 i 回路 约束的含义：每个回路至少要断裂一次 xj = aij = 对于准则①： 令 pj = 1 ，即目标是使断裂流股数目最少； 对于准则②： 令 pj = 变量数。 对于准则③： pj 根据流股性质赋值。例如依据流股变量对计算过程灵敏度大小的估计 赋值：灵敏度大，则流股变量的微小变化便会引起整个模拟过程较大的 波动，使计算不易收敛，这时该流股的权就应取较大的值。 对于准则④