Slide 3
VJ:冷却水在充满过程中在夹套中的体积。 VJTJ:冷却水充满过程中的模型方程所用的项； RAMP：给定值PSET随时间的变化率； PRINT：打印语句，循环之前打印出大的标题
Slide 4
Slide 5
从 100 行 开 始 进 入积分循环计算 过程 1: 计 算 速 率 常 数 K1,K2。 2. 变 送 器 方 程 ： 求出PTT； 3. 控制器方程： 由变送器来的温 度信号求出其输 出气压信号PC。 4. 阀门方程：求 出蒸汽阀和冷却
Slide 1
结束
间பைடு நூலகம்反应器
ALPHA1,ALPHA2:指前因子，分别用于求取反应 速率常数K1,K2。 AVP,BVP:组分的饱和蒸汽压方程系数（由温度求 蒸汽压) DELTA：积分时间步长；
Slide 2
KC:控制器方程系数。 TPRINT:用来控制打印时间的间 隔； PJ：夹套压力； DENS:由TJ,PJ计算出水蒸气密度； PSET:给定气动压力信号的初值。 START：指示是否加热或冷却； 初值1；在随后的程序中被修 改； FWO：冷却水流量； FULL：指示夹套是否被冷却水充 满
}
间歇反应器作业
计算框图：
赋初值和给定值 进入循环 计算反应速率常数K1,K2 列出变送器、控制器、阀门方程，最终 计算出阀门开度XW,XS 判断夹套加热或冷却,START>0加热；<0冷却
蒸气加热的方程组
冷却过程（判断是否充满FULL） 充满前方程 充满后方程
反应液和金属壁的方程组 进行积分得到新值 否，继续循环 判断时间是否停止
Slide 6
从100行开始进入积分循环计 算过程 6.如果加热，则写出加热时的 方程；求出新的密度值DENS后 ，进行一个迭代循环过程计算 TJ，然后由TJ计算PJ。
Slide 7
从100行开始进入积分循环计算 过程 7.如果冷却，则写出冷却时的 方程；首先求出由阀门开度求 出冷却水的流量，令与水蒸气 加热相关变量=0； 然后利用FULL变量值判断是否 处于充满状态；如果未充满， 则计算传热面积AO，冷却水体 积VJ新值（积分），以及传热 量QJ，导数项VJTJ的新值，并 求出TJ。
Slide 11
• do • { • //计算速率常数 • k1 = a1* pow(e,-E1/(1.99*(T+460))); • k2 = a2* pow(e,-E2/(1.99*(T+460))); • • //需要写出变送器、温度控制器、两个阀门的方程
• • • • • • • • • • • • • • • • • • // XS:蒸气阀门开度, XW：冷却水阀门开度 //判断夹套处于加热或冷却过程, 利用START变量 if(START>0)//加热过程 { //写出水蒸气加热的加套的方程组 } else//冷却过程 { //判断是否充满 if(FULL>0)//已经充满 {//写出充满后的方程 } else//未充满 { //写出未充满时的方程组 Slide 12 }
从100行开始进入积分循环计 算过程 5: 利用START的值判断处于加 热或冷却过程。。 6.如果加热，则写出加热时的 方程；首先求出蒸汽流量（由 阀门开度XS和压力差35-PJ）， 以 及 传 热 量 QJ 和 冷 凝 液 的 量 WC，然后就可以写出常微分 导数项方程DENDOT,并进行积 分求出新的密度值DENS。
从100行开始进入积分循环计算过程 9. 对CA,CB,T,TM利用导数项CADOT,CBDOT,TDOT,TMDOT进行积分，得到新的CA,CB,T,TM 值。
Slide 10
从100行开始进入积分循环计算过程 10.进行判断，是否达到结束时间，打印输出语句，并且对时间进行累加，最后返回 100行语句循环计算 。
Slide 8
从100行开始进入积分循环计算过程 8.如果已经充满，则按充满后的夹套方程写出，首先求出Q传热量（后面夹套能量衡 算方程要用到Q）。求出导数项TJDOT，并进行积分得到新的TJ值。然后就可以写出 CA,CB,T,TM的导数项CADOT,CBDOT,TDOT,TMDOT。
Slide 9