第5章 反应器的优化和稳定性.ppt

第五章 反应器的优化和稳定性;第五章 反应器的优化和稳定性;5.1 目标函数和约束条件;5.1.2 约束条件 反应器所特有的及反应过程存在的特定关系式（与上述三种变量间）称为约束条件。 分为两种： 等式约束：过程方程式（物料平衡、速度式、k式等）等 不等式约束：温度极限等、不等式。;将等式约束写成 yi = ?i（x1 , x2 , x3 , … , xm）（写成操作变量的函数） 于是目标函 P = P（ x1 , x2 , x3 , … , xm）（消去 yi） 若 P 和?P /?xk 是连续的，则在满足;下面举例说明应用;间歇反应器， t 为操作时间;可解得TOPT ，因很繁，仅写出;亦即;令其为“0”，得;微分法要求从目标函数中消去状态变量 yi 当 yi 不能写成操作变量 xi的显函数时，不能用微分法，此时，可用拉格朗日乘数法 设新的目标函数为：;未知量(2n+m) 个;由 R 物料平衡;问题分为： （1） 温度T 一定，求最佳θ; （2） θ 一定，求最优温度 Topt;例如：;即搜索矢量; 例3 讨论日产10吨化工厂生产费用最小值 图e.5.3.a流程图;出发点任选，设 x1 =800，x2 =6，P0 =3133333;以后逐步计算，直到 ?P / ?xi=0为止，即为极值点。;步长很重要，P 曲面复杂时要选小些，曲面简单时可大些； 当存在两个以上极值时，一旦到达其一就不能前进，为此应选择几个出发点； 到达是极值点还是鞍点，落入鞍点就无法解脱； 变步长方法： 一次成功，下次取3倍步长； 一次失败，下次取1/2步长。;5.2.4 线性规划（Linear Programming）; 不同 K 值构成一个平行线族，线愈高，则K值愈大，但 K 必须位于“可行面积”（可行域）内 ;引入松弛变量 Zi，把不等式化为等式;写为矩阵形式;由于 fm+1，…，fm+l 不为“0”;Dantzig 程序 （计算方法）（运筹学中，迭代或旋转运算） 因方程只有 l 个，变量为 m+ l 个，故先假定 l 个变量然后求解。 （1）（设 f1 = f2 = f3 = … = fm =0）求初始可行解 根据（5.45）式， f1 … fm 为“0”，又am+1 … am+l 为1，构成单位向量;解释 ： （i）主成分（基变量） 由上可见，从原点出发可理解为选定 l 维坐标系，使符合初始条件。所以， am+1 , am+2 , … ,am+l应满足条件 ;（5.52）;举例说明： ;（a）-（b）/ 2 得 ;3 次可行解 ;适当选择 θ 值，如;下面举例说明（见 p155 ~p160）;;③第2可行解（将 f1 系数向量 a1化为单位向量）判别换入换出变量;④求第3可行解，其它同前;⑤求第4可行解;上表得，第4可行解;另外，根据（5.59）式;目标函数;泛函的极值条件是变分为0：即;例6;由三角公式;积分可得;该泛函达极值的必要条件为：;例7. 圆形管道横断面旋度平方的面积分为最小（条件），平均流速 U0，求速度分布。;旋度平方的面积分的泛函形式;积分;可求得;5.2.6 最大原理;现以数学符号引入另一个函数 J;（5.70）;说明：;对（c）Talor 级数展开;上式很复杂，可适当选取 λi 使其简化;说明：以上为问题最优解的必要条件，但不是充分条件。可能有多个解满足上述条件，只有一个为最优。也可能因约束条件限制使过程不能达最佳状态。这时可用彭策根（Pontryagain）最大原理。;将共轭函数（5.78）和极值必要条件（5.80）代入（5.83）得;解：;本题中;式中;（1）;5.2.7 动态规划;由流程图看：改变萃取剂量 w1 y1，x1变 x2 , … , xN 及 y2 , … , yN 都发生变化；同理，改变 wn 也会使其 yn 以及其后各段的 x 和 y 都发生变化。构成 N 变量问题。若用传统方法，如微分法，首先要给出总收益的目标函数，即使这样，解起来也相当复杂。;第N-1段：对N段而言，进液浓度xN-1 是唯一受上段影响的因素，因此N段的最优决策是针对xN-1而言的。 xN-1一定，N段决策不变, 也就是以N的最优决策对N-1段。 改变N-1段条件（ wN －1），可以求得2段过程的最优决策,其余类推。;用于第N段，目标函数为;如果视 m , λ为