化工过程过程系统的模拟优秀课件.ppt

《化工过程过程系统的模拟优秀课件.ppt》由会员分享，可在线阅读，更多相关《化工过程过程系统的模拟优秀课件.ppt（38页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、化工过程过程系统的模拟化工过程过程系统的模拟第1页，本讲稿共38页 （2）过程系统的设计（Design Problem）过程系统模型控制模型输入流股向量设备参数向量输出流股向量设计要求指标可调设备参数向量可调输入向量输出设计结果向量第2页，本讲稿共38页 （3）过程系统的优化过程系统模型输出流股向量经济分析模型优化程序给定输入给定参数经济参数输出优化结果优化变量性能指标约束条件约束条件第3页，本讲稿共38页（1）图形表示）图形表示 工艺流程图的有向图或信息流程图（Information flow diagram)基本概念：节点 设备单元 边 流股 子图 路径 循环回路或环路1271211121

2、4354678910111213683549102.1 过程系统结构的表达第4页，本讲稿共38页（2）矩阵表示）矩阵表示 （a）过程矩阵）过程矩阵(ProcessMatrix)Rp 表达过程系统单元设备与流股之间的关系，由流股将相关 设备关联起来。单元设备序号相关物流号流入该节点的流股+流出该节点的流股-12341112-1 1-2第5页，本讲稿共38页 (b)邻接矩阵（Adjacency Matrix)RA 一个由n个单元或节点组成的系统，其邻接矩阵或相邻矩阵 可表示为nn的方阵。1 2 3 4 5 6 7 8 9 10 11 12流入节点j流出节点i121112111011第6页，本讲稿共

3、38页RA=j 1 2 3 4 5 6 7 8 9 10 11 12i2143576891012110 1 0 0 0 0 0 0 0 0 0 00 0 1 0 0 0 0 0 0 0 0 00 0 0 1 0 0 0 0 0 0 0 10 0 0 0 1 0 0 0 0 0 0 00 0 0 0 0 1 0 0 0 0 0 00 0 0 0 0 0 1 0 1 0 0 00 0 0 0 0 0 0 0 1 0 0 00 0 1 0 0 1 0 0 0 0 0 00 0 0 0 0 0 0 0 0 1 0 00 0 0 0 0 1 0 0 0 0 1 00 0 0 0 0 0 0 0 0 0

4、0 00 0 0 0 0 0 0 0 0 0 0 0空的列（元素都为零）：系统中没有输入的节点；空的行（元素都为零）：系统中没有输出的节点。邻接矩阵：第7页，本讲稿共38页（c）关联矩阵（Incidence Matrix)RI元素Sij=-1，边（流股）j为节点（单元设备）i的输出流股 1，边（流股）j为节点（单元设备）i的输入流股 0，边（流股）j与节点（单元设备）i无关联RI=j 1 2 3 4 5 6 7 8 9 10 11 12 13 14i21435 7689101211-1 0 0 0 0 0 0 0 0 0 0 0 0 01-1 0 0 0 0 0 0 0 0 0 0 0 00

5、1 1 0 0 0 0 0 1 0 0 0 0 -10 0 1-1 0 0 0 0 0 0 0 0 0 00 0 0 1-1 0 0 0 0 0 0 0 0 00 0 0 0 1-1 0 1 0 -1 0 1 0 00 0 0 0 0 1-1 0 0 0 0 0 0 0 0 0 0 0 0 0 1-1-1 0 0 0 0 00 0 0 0 0 0 0 0 0 1 -1 0 0 00 0 0 0 0 1 0 0 0 0 1 -1 -1 00 0 0 0 0 0 0 0 0 0 0 0 1 00 0 0 0 0 0 0 0 0 0 0 0 0 1关联矩阵第8页，本讲稿共38页关联矩阵具有以下性质：

6、关联矩阵具有以下性质：.若有向图中有n个节点m条边，则关联矩阵为n行m列的矩阵。.每一流股(边)在矩阵中标出两次，即同一条边可是一个节点的输出又是另一节点的输入边。.列的元素之和为零。第9页，本讲稿共38页3.3过程系统的分解过程系统的分解将一个结构已定的系统分割成一些更小的次一级系统的将一个结构已定的系统分割成一些更小的次一级系统的方法。方法。将系统的总目标分解成更小的系统的目标，或者将阶将系统的总目标分解成更小的系统的目标，或者将阶数、维数数、维数 很大的系统的数学模型分解成阶数、维数较小的子很大的系统的数学模型分解成阶数、维数较小的子系统的数学模型。系统的数学模型。分解的目的：分解的目的

7、：降低计算复杂度，提高计算效率。降低计算复杂度，提高计算效率。第10页，本讲稿共38页3.2.1 问题的提出分解的必要性：所有方程联立求解困难分解的可能性：每一个方程并不含所有变量 矩阵的稀疏性 系统分解（Decomposition)步骤：（1）系统的分隔（或分割，Partitioning)识别独立的子系统 从子系统中识别循环回路或最大循环网（2）子系统（循环回路或最大循环网）的断裂第11页，本讲稿共38页3.2.2不相关子系统的识别不相关子系统的识别（了解）（了解）可分为可分为2个子系统：个子系统：第12页，本讲稿共38页 3.2.3 对不相关子系统的分隔 在不相关子系统中识别出不可再分隔的

8、子系统，即循环回路及 最大循环网，并用拟节点表示，然后按信息流方向排出有利的 计算顺序。ABCDEDFCBEA最大循环网最大循环网（包含（包含2 2个关联的循环回路）个关联的循环回路）2 2个序贯相连的循环回路个序贯相连的循环回路第13页，本讲稿共38页直观分析法：直观分析法：H单独单独1组；组；A，B，C，D，E构成构成1组；组；F，G构构成成1组；组；I单独单独1组。组。计算顺序：计算顺序：H，(A,B,C,D,E)，(F,G)，I第14页，本讲稿共38页 3.2.4 最大循环网的断裂 选择最优断裂流股的准则：I.断裂的流股数目最少；II.断裂流股包含的变量数目最少；III.对每一流股选定

9、一个权因子，该权因子数值反映了断裂该流股时迭代计算的难易程度，应当使所有的断裂流股权因子数值总和最小；IV.选择一组断裂流股，使直接代入法具有最好的收敛特性。第15页，本讲稿共38页u Lee-Rudd断裂法 该法属于第I类最优断裂准则，即断裂的流股数目最少，把一最大循环网所包含的所有回路打开。有四个回路A，B，C，D及8个流股。第16页，本讲稿共38页S1 S2 S3 S4 S5 S6 S7 S80 1 1 0 0 0 0 01 1 0 1 0 0 0 00 0 0 0 0 0 1 10 0 0 1 1 1 1 01 2 1 2 1 1 2 1ABCDf2432R回路矩阵的元素定义：Cij=

10、1 物流Si在回路i内时0 物流Si不在回路i内时f：回路频率，某一流股出现在各回路的次数。R：回路的秩，某一回路中包含的流股总数。其相应的回路矩阵（Loop matrix)为：第17页，本讲稿共38页步骤：I.除去不独立的列除去不独立的列k对于第对于第j 列与第列与第k 列，若流股频率列，若流股频率fj fk 成立，且成立，且k 列中非零列中非零值的行对应列值的行对应列j 的行也为非零值，则列的行也为非零值，则列k 不是独立的，为列不是独立的，为列j 所包含。所包含。S2 S4 S7ABCD1 0 00 0 11 1 00 1 1R1122第18页，本讲稿共38页II.选择断裂流股选择断裂流

11、股剩下的独立列构成的回路矩阵中，秩为剩下的独立列构成的回路矩阵中，秩为1的行说明该行所的行说明该行所对应的回路只剩下一股物流，为此打开该回路，必须将该行非零对应的回路只剩下一股物流，为此打开该回路，必须将该行非零元素对应的流股断裂。元素对应的流股断裂。断裂断裂S2，A、C打开；断裂打开；断裂S7，B、D打开。打开。S2 S4 S7ABCD1 0 00 0 11 1 00 1 1R1122第19页，本讲稿共38页计算顺序图示:第20页，本讲稿共38页思考题：最大循环网如何断裂？如何断裂？第21页，本讲稿共38页ABCDE ab cds1s3s5s8s4s13s7s6s1s3s4s5s6s7s8s

12、13Ra1113b1113c1113d112f12121211第22页，本讲稿共38页s1s3s4s5s6s7s8s13Ra1113b1113c1113d112f12121211第23页，本讲稿共38页s3s5s7Ra11b112c112d11断裂断裂S3，a、b打开；断裂打开；断裂S7，c、d打开。打开。第24页，本讲稿共38页是否满足是否满足第25页，本讲稿共38页 3.3 化工流程模拟计算收敛方法系统经过分隔和循环网的断裂后，给定初值，模拟计算时需要选择有效的收敛算法。当 或时，即得到收敛解 第26页，本讲稿共38页 流程模拟涉及的基本概念：（1）化工流程模拟 应用计算机作为辅助手段，利

13、用过程模拟软件对一个化工过 程进行稳态的热量和物料衡算或设备尺寸计算和费用计算，以获得过程流程系统的物性和行为。（2）隐式表达形式和显示表达形式方程形式叫做方程的显式表达形式。方程形式叫做方程的隐式表达形式。（3）局部收敛（local convergence）迭代求解不能保证收敛到真实解的特性就叫做局部收敛。第27页，本讲稿共38页 （4）全局收敛 (global convergence)对于迭代求解，如待求解的非线性方程无论只有一个解还是多 个解，算法均能保证方程的求解收敛在唯一正确的解时，则称 迭代求解具有全局收敛性。（5）收敛判据（convergence criterion)用来判定迭代

14、计算收敛精度的目标函数值称之收敛判据。或也可按相对量考虑而提出如下收敛判据：第28页，本讲稿共38页（6）收敛容差（convergence tolerance）在方程的迭代求解过程中，在收敛判据中设定的前后两次 迭代结果的差值，就叫做收敛容差，也称收敛误差。收敛 容差一般用来代表。通常为一个足够小的正数。或第29页，本讲稿共38页（7）收敛速度（convergence speed)求解方程的任何迭代法的收敛速度可用下式来衡量：lim X*是它的解，n和C都是正数。指数n愈大，收敛速度愈快。通常将n=1和n=2所对应的情况分别称为收敛速度具有线性收敛（linear convergence）和二次

15、收敛（quadratic convergence）的性质。若n比1大一些，则称为超线性收敛（superlinear convergence）。第30页，本讲稿共38页 3.3.1 直接迭代法（direct substitution method）求解显式方程式的最简单的一种迭代方法：直接迭代法比较广泛地用于流程模拟计算中，当初值选得较好时是会收敛的，但其收敛速度较慢。第31页，本讲稿共38页 3.3.2 部分迭代法（partial substitution method）其迭代公式为：或写成：w是用来调节两部分大小的一个系数，叫松弛因子。实际使用部分迭代法时，要对w的数值进行合理的估计。第32

16、页，本讲稿共38页 3.3.3 韦格施坦法（Wegstein method）其迭代公式为：其中：此法的收敛速度，具有超线性收敛的性质，比部分迭代法（包括直接迭代法）快。需设置两个初始点，但如果在第一轮迭代中采用直接迭代法，从第二轮开始再改用韦格施坦法，则只需设置一个初始点即可迭代求解。第33页，本讲稿共38页 3.3.4 牛顿拉夫森法(Newton-Raphson method)对于非线性方程组：在 处作泰勒展开，只截取一次项，则可得如下的方程形式：记作 称雅可比矩阵,则可得如下方程：上式为一线性方程组，于是，可得牛顿拉夫森法迭代公式为：第34页，本讲稿共38页牛顿拉夫森法的收敛速度很快，具有

17、二次收敛性。3.3.5 拟牛顿法（quasi Newton method）设代替雅可比矩阵逆阵的矩阵为：则拟牛顿法的迭代公式为：该法初值要求不高，收敛速度快，收敛性能大为改善。第35页，本讲稿共38页 3.3.6 各种计算收敛方法的比较 上面介绍了几种计算收敛方法。这里应特别提到另一种处上面介绍了几种计算收敛方法。这里应特别提到另一种处 理方法理方法“搜索法搜索法”。也就是说，如果把非线性方程组的求解当成一个最优化问题。也就是说，如果把非线性方程组的求解当成一个最优化问题来处理，则可以用整套非线性规划法来解非线性方程组问题，如最速下降法、来处理，则可以用整套非线性规划法来解非线性方程组问题，如

18、最速下降法、单纯形法等等单纯形法等等 这类搜索法求解特别适合于初始值较差的情况下。由于初始点距离解太远，用别的方法引起不收敛时，采用这种方法效果比较好。对于一个给定的问题，想选择一次迭代计算量最少的方法。但是，也必须考虑收敛速度。往往先试计算量最少的方法，如果失败了，再试计算量大但收敛性较好的方法。第36页，本讲稿共38页几种计算收敛方法的比较几种计算收敛方法的比较 方法预测搜索方向预测搜索距离是否要求导数是否能保证Q总是局部下降由基准点预测新方向时要算多少次函数值相对收敛速度直接迭代法 是 是 否 否 1 慢牛顿拉夫森法是 是 是 是 n 快韦格斯坦法 是 是 否 否 1 较快拟牛顿法 是 是 否 否 1相当快第37页，本讲稿共38页 3.4 过程系统模拟基本方法 序贯模块法、联立方程法和联立模块法。序贯模块法：逐个单元模块依次序贯计算求解系统模型的一种方法。联立方程法：将描述过程系统的所有方程组织起来，形成一大型非线性方程组，进行联立求解。联立模块法：取序贯模块法及联立方程法两者之长，模拟过程分两个层次：单元模块的层次，系统流程的层次。第38页，本讲稿共38页