动态模拟及解析总结计划.pptVIP

* 第六章 化工过程动态模拟与分析 主要内容 一、化工过程的动态模型 二、CSTR的动态模型 三、状态空间分析法 四、常微分方程（组）求解 6.1 化工过程的动态模型 1、动态特性是化工过程系统最基本的特性之一 间歇过程、开（停）工过程、系统的老化… 2、化工过程动态模型分类 模型类型 数学形式 应用实例 集中参数模型 代数－常微分方程组 理想搅拌罐反应器动态模型 分布参数模型 代数－偏微分方程组 管式反应器动态模型 多级集中参数模型 代数－常微分方程组 串联CSTR动态模型 混合模型 上述方程组的混合形式 多个单元过程组合而成的系统 3、确定性动态模型的数学处理 （1）正问题－模型方程组的求解； （2）逆问题－模型参数的估计； （3）过程系统的定性分析。 定常态的多重性和稳定性、参数敏感性… 6.2 常微分方程的数值求解 一、一阶常微分方程初值问题 计算出解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 节点间距 为步长，通常采用等距节点，即取 hi = h (常数)。 x0 x1 I 欧拉法 ? 欧拉公式： 向前差商近似导数 记为 x0 x1 x2 y0 y1 y2 II 隐式欧拉法 向后差商近似导数 )) ( , ( ) ( 1 1 0 1 x y x f h y x y + ? x0 x1 ) 1 , ... , 0 ( ) , ( 1 1 1 - = + = + + + n i y x f h y y i i i i 一般先用显式计算一个初值，再迭代求解。 III 改进欧拉法 迭代求解 IV 龙格－库塔法 将每两点之间分为若干点，以各点的处的斜率平均值作为中值，并根据误差分析，可以给出下列常微分初值问题的数值计算公式： 例 6.1 求解如下常微分初值问题： t = [0:10:100]; h = 10; T0 = 93.3; Ts = zeros(1,11); Ts(1) = T0; for i = 1:10 T_ = Ts(i) + h*f(t(i), Ts(i)) ; Ts(i+1) = Ts(i) + 0.5*h*( f( t(i), Ts(i) ) + f( t(i+1), T_ )); end 改进欧拉法程序： t = [0:10:100]; h = 10; T0 = 93.3; Ts = zeros(1,11); Ts(1) = T0; for i = 1:10 xi = t(i); yi = Ts(i); K1 = f(xi,yi); K2 = f(xi+h/2, yi+h*K1/2); K3 = f(xi+h/2, yi+h*K2/2); K4 = f(xi+h, yi+h*K3); Ts(i+1) = yi + h*(K1+2*K2+2*K3+K4)/6; end 龙格-库塔法程序： 二、一阶常微分方程组初值问题 一阶常微分方程组： 一阶常微分方程初值问题的改进欧拉法： 一阶常微分方程组初值问题的改进欧拉法 例 6.2 某一温度为T=360 K的间歇式反应器中进行反应 A→ B→ C，设xA、xB分别为A、B的浓度，反应速率方程如下，试求时间段[0，1]反应器中A、B浓度变化情况。 xA0 = 1; xB0 = 0; % 初始值 t = 0:0.1:1; % 时间划分 h = 0.1; % 步长 A = zeros(1,11); A(1) = xA0; B = zeros(1,11); B(1) = xB0; for i = 1:10 [dA,dB] = fs(A(i),B(i)); A_ = A(i) + h*dA; B_ = B(i) + h*dB; [dA_,dB_] = fs(A_, B_); A(i+1) = A(i) + 0.5*h*(dA + dA_); B(i+1) = B(i) + 0.5*h*(dB + dB_); end function [dA,dB] = fs(xA,xB) T = 300; % 反应温度 dA = -4e3*exp( -2500/T )*xA^2; dB = 4e3*exp( -2500/T )*xA^2 - 620000*exp( -5000/T )*xB; 三、高阶常微分方程组初值问题 思路 转化为一阶常微分方程组求解 可用改进欧拉法或龙格－库塔法求解。 例 6.4 求解常微分初值问题：如图所示储水罐，底面积为A＝2.5，初始高度H0=1，进水量Fi=10，k=0.04，求液面