第二章化工系统的定常态模拟与分析2.pptVIP

阻尼因子q的取值不当可使迭代计算收敛缓慢甚至发散 Wegstein法虽然无须人为选定q值，但是也会因为q值不当导致坏的收敛行为 有界Weigstein法就是凭借经验人为地把q值限定在一定的范围内，以改善收敛行为 有界Wegstein法 FLOWTRAN流程模拟系统中取为-5、0 CHESS系统中当q0或q-10时令q=0 * * * S1 S2 S3 S1 S4 S5 S3 S3 S2 {S1 S2 S6 S7} S3 S3 {S1 S4 S5 S6 S7} S1 S4 S7 {S2 S5 S6} S5 S6 S1 S4 S7 S1 S1 S4 S7 {S2} {S3 S4 S5} S2 S3 {S4 S5 S6 S7} S5 S6 {S1 S4 S7} 非多余断裂族: 断裂组｛S1，S4，S7｝为最优断裂组。 通过断裂可以把不可分隔子系统中的回路物流打开，从而可以利用序贯模块法对该过程系统进行模拟计算。这种模拟计算的开始是首先要设定起始物流变量的猜值，计算的终点则在于该猜值与计算值的收敛。 物流 权重因子 图2-13 不可分隔子系统 6 4 2 3 7 1 5 单元1 单元2 单元3 单元4 (2) (3) (4) (9) (3) (2) (2) 权重因子 回路 流股 用Upadhye-Grens断裂法求图2-13中的最优断裂组 单元 输入流 输出流 1 5，6 1 2 1，4，7 2 3 2 3，4，5 4 3 6，7 从断裂组（S1,S2,S3）开始，反复利用替代规则 S1 S2 S3 S1 S4 S5 S3 S3 S2 {S1 S2 S6 S7} S3 S3 {S1 S4 S5 S6 S7} S1 S4 S7 {S2 S5 S6} S5 S6 S1 S4 S7 S1 S1 S4 S7 {S2} {S3 S4 S5} S2 S3 {S4 S5 S6 S7} S5 S6 {S1 S4 S7} 非多余断裂族: 断裂组｛S1，S4，S7｝为最优断裂组。 通过断裂可以把不可分隔子系统中的回路物流打开，从而可以利用序贯模块法对该过程系统进行模拟计算。这种模拟计算的开始是首先要设定起始物流变量的猜值，计算的终点则在于该猜值与计算值的收敛。 2.2.3 断裂物流变量的收敛 执行断裂物流变量收敛功能的模块称收敛单元模块 X 是断裂物流变量的猜值，y是经过程系统模型计算得出的断裂物流变量的计算值。断裂物流变量的收敛问题实际上是个迭代求解非线性方程组的问题 x = y = G(x) 图2－15不可分隔子系统的断裂物流 　 过程系统 (a) 断裂流股 (b) 当断裂物流变量猜值x与计算值y之差小于收敛容差ε时 y － x = G（x）－ x ε 则x为断裂物流变量的收敛解 收敛单元的功能总计有如下作用： (Ⅰ)获取猜值的初值x0 (Ⅱ)根据计算值y，以一定的方法确定新的猜值x (Ⅲ)比较猜值x和计算值y，若其结果满足给定精度要求则结束迭代计算，否则继续迭代计算过程 可见，收敛单元实质上就是一个数值迭代求解非线性方程组的子程序。求解非线性方程组的数值计算方法很多，适合于收敛单元的数值计算方法一般应尽可能满足下列要求： 1. 对初值的要求不高。 1)初值易得，不易引起迭代计算的发散； 2)初值的组数少。 2. 数值稳定性好 2. 数值稳定性好 2. 数值稳定性好 3. 收敛速度快 对收敛速度的影响主要有三个因素：迭代次数；函数G（x）的计算次数；矩阵求逆的次数。 每计算一次函数值就相当于做一次流程回路的模拟计算,每求一次导数就要做两次流程模拟计算 对于断裂物流的收敛，好的非线性方程组的数值迭代次数少，而且应该尽量避免导数计算和矩阵求逆 4. 占用计算机存储空间少 直接迭代法 直接迭代法是将计算值yk作为下一轮迭代的猜值xk+1而实施迭代计算 xK+1 = yK, 即xk+1＝G（xk） 非线性方程组x = y = G(x) 的另外一种形式为 F(x)＝x－G(x)= 0 → xk+1＝xk＋F(xk) 与牛顿公式相比较： 直接迭代法的雅可比矩阵为单位矩阵 优点：方法简单， 只需要一组初值，不需计算导数和逆矩阵 弱点：迭代次数多、收敛速度慢， 对初值要求较高 例: 用直接迭代法求解下列方程组 解：令猜值为X1＝2；X2＝10；X3＝5 解：令猜值为X1＝6；X2＝3.5；X3＝5 为改善直接迭代法的收敛行为，提出了阻尼直接迭