第九章 化学二元反应蒸馏塔仿真mathematica实现程序-马耀元.docxVIP

(*已知参数*)F = 100;xF = 0.5;xD = 0.95;xW = 0.05;q = 1;han = 6944;af = 2; DD = 50;W = 50;erro = 10^-8;time = 20;cur = 15;fun = NSolve[{yq*(q - 1) == q*xq - xF, yq*((af - 1)*xq + 1) == af*xq}, {xq, yq}];Rmin = FromDigits[(xD - Evaluate[yq /. fun])/(Evaluate[yq /. fun] - Evaluate[xq /. fun])];(*************************************控制部分*******************************************)(*可知F=D +W,所以w=0,newton法在此建立，即w-50=0,而w-time=LN-VN=L1+qF-Q/han-time=0=\R*D+q*F-Q/han-50=0,即： R*D+q*F-Q/han-W=0 \,其中，R,Q的值影响整个反应体系，具体分离目标和它俩的取值有关，利用newton-rapphson迭代法xn+1=xn-f \(xn)/f(xn)求取合适的R和Q 。因为R的取值范围一般为R=1.1~2.0Rmin,而Rmin=(xD-Yq)/(Yq-Xq) ,因此尽量满足R的值在此范围内*)pc = 1/3;ic = 1;dc = 1;R = 2;erropre1 = 5;erropre2 = 5;errold = 0;flag = 1;caltime = 0;failtime = 0;While[flag != 0, ++caltime; L1 = R*DD; LN = L1 + F*q; Q = (LN - W)*han; VN = Q/han; V1 = VN + F*(1 - q); (*物料平衡微分方程36块板，一个冷凝器，一个再沸器，共38个微分方程*)function = NDSolve[ {(*初始化方程，初始值可以任意给定*)x1[0] == 0, x2[0] == 0, x3[0] == 0, x4[0] == 0, x5[0] == 0, x6[0] == 0,x7[0] == 0, x8[0] == 0, x9[0] == 0, x10[0] == 0, x11[0] == 0, x12[0] == 0,x13[0] == 0, x14[0] == 0, x15[0] == 0, x16[0] == 0, x17[0] == 0, x18[0] == 0,x19[0] == 0, x20[0] == 0, x21[0] == 0, x22[0] == 0, x23[0] == 0, x24[0] == 0,x25[0] == 0, x26[0] == 0, x27[0] == 0, x28[0] == 0, x29[0] == 0, x30[0] == 0,x31[0] == 0, x32[0] == 0, x33[0] == 0, x34[0] == 0, x35[0] == 0, x36[0] == 0,x37[0] == 0, x38[0] == 0, (*冷凝器方程*)x1[t] == V1*(af*x2[t])/((af - 1)*x2[t] + 1) - L1*x1[t] - DD*x1[t], (*1-36块板的物料平衡方程*) (*精馏段方程*)x2[t] == V1*(af*x3[t])/((af - 1)*x3[t] + 1) + L1*x1[t] - L1*x2[t] - V1*(af*x2[t])/((af - 1)*x2[t] + 1),x3[t] == V1*(af*x4[t])/((af - 1)*x4[t] + 1) + L1*x2[t] - L1*x3[t] - V1*(af*x3[t])/((af - 1)*x3[t] + 1),x4[t] == V1*(af*x5[t])/((af - 1)*x5[t] + 1) + L1*x3[t] - L1*x4[t] - V1*(af*x4[t])/((af - 1)*x4[t] + 1),x5[t] == V1*(af*x6[t])/((af - 1)*x6[t] + 1) + L1*x4[t] - L1*x5[t] - V1*(