序贯二次规划西安交通大学材料.pptVIP

5.化工中的应用实例—求解结果 第一次迭代:取初始点x0=(54.0,38.0,3628.0,1.01,0.51)T, 表5-9为模拟结果(略去量纲)。 Company name 5.化工中的应用实例—求解结果 Company name 5.化工中的应用实例—求解结果 Company name 5.化工中的应用实例—求解结果 Company name 5.化工中的应用实例—求解结果 Company name 5.化工中的应用实例—求解结果 Company name 5.化工中的应用实例—求解结果 Company name 化工学院 * * * * * * * * * * * * * * * 步长α的选取应使得精确罚函数(5-66)的值下降, 即 这就可以用以前那种精确的一维搜索对罚函数 求极小来确定步长α。 Company name 3.一维搜索—确定步长 Company name 4.矩阵Qk的校正 Company name 4.矩阵Qk的校正 Company name 4.矩阵Qk的校正 Company name 5.化工中的应用实例 有一苯、甲苯闪蒸过程, 苯、甲苯混合物通过一单级绝热闪蒸器进行分离, 如图5-9所示。为提高苯的回收率,闪蒸器底部物流④一部分作为产品⑤, 另一部分作为循环流⑥经泵与新鲜料液混合后再次进入闪蒸器。试调整闪蒸压力p及分流系数s( =物流⑤的流量/ 物流④的流量),使气相产物苯的产量最大。 max Company name 5.化工中的应用实例 已知进料状态见表5-8。 且可认为混合物为一理想体系。 Company name 5.化工中的应用实例—数学模型 Company name 5.化工中的应用实例—数学模型 Company name 5.化工中的应用实例—数学模型 Company name 5.化工中的应用实例—数学模型 Company name 5.化工中的应用实例—不可行路径法的求解思想 上述原模型原则上可用序贯二次规划法求解。然而,由于化工生产的复杂性, 一般流程优化问题的数学模型 min(max) f(x) s .t . g(x) ≥0 (5-95) h(x) =0 其规模相当庞大，常有数以万计的过程变量x和非线性很强的过程方程式h(x) =0以及技术要求不等式g(x)≥0, 从以上简单的闪蒸问题也可以看出这一点。 Company name 5.化工中的应用实例—不可行路径法的求解思想 Biegler 和Hughes 提出了一种将优化与序贯模块法模拟结合的方法。参与调优的变量或决策变量的数目相对来说是比较少的,一般只占2‰～3 ‰, 剩下的变量大多只是模拟过程中的内部变量。 为此,可将目标函数及约束函数中的变量分为两部分, 一部分叫调优变量,仍记作x, 另一部分叫从属变量,记作r, 即 x=(x, r ) Company name 5.化工中的应用实例—不可行路径法的求解思想 一部分由原模型中的目标函数、技术要求不等式约束g(x, r)≥0和切割方程(若有回路)C(x, r)≥0组成的调优模型,即 min(max) f(x, r) s .t . g(x, r) ≥0 C(x, r) ≥0 (5-96) 另一部分则由原模型中的过程方程式构成的过程模拟模型,仍记作: h(x, r) =0 (5-97) Company name 5.化工中的应用实例—调优模型和过程模拟模型 (1) 调优模型 考虑到数学上叙述的方便和与前面序贯二次规划的连贯性,我们不妨设: x1 ———切割变量,物流⑦中苯流量的估算值; x2 ———切割变量,物流⑦中甲苯流量的估算值; x3 ———切割变量,物流⑦焓值的估算值; x4 ———调优变量,闪蒸压力; x5 ———调优变量, x5 =F5/ F4 ,分流系数; x6 ———物流⑦中苯流量的计算值( x6 =F2 C71 ); x7 ———物流⑦中甲苯流量的计算值( x7 =F7 C72 ); x8 ———物流⑦的焓值计算值( x8 =I7 )。 选定闪蒸压力x4 和分流系数x5 为调优变量, 则以气相物流③中苯的流量F3 C31为最大的并与式(5-96)相对应的调优模型为 Company name 5.化工中的应用实例—调优模型和过程模拟模型 min f