cooleroptdes.m

matlab在化学工程中的应用实例程序,对大家研究学习有指导作用。

function CoolerOptDes
% 冷却器的最优化设计(Optimal Design of a Cooler )
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2002/05/6 $

clear all
clc

global T1 T2 G t1 JA beta K theta Jw Cw Cc Q
T1 = 140;       % ℃
T2 = 40;        % ℃
G = 2e4;        % kg/h
t1 = 30;        % ℃
JA = 200;       % yuan/m^2
beta = 0.15;
K = 836.8;      % kJ/(m^2 h ℃)
theta = 8000;   % h
Jw = 0.04;      % yuan/ton
Cw = 4.184;     % kJ/(kg ℃)
Cc = 2.092;     % kJ/(kg ℃)

Q = G*Cc*(T1-T2);
t0 = 50;        % 给定初值

t2 = fminsearch(@TotalFee,t0); 
fprintf('   优化结果：\n\n')
fprintf('冷却器最优出口温度为：%.2f %s\n',t2,'℃' )
allFee = TotalFee(t2);
fprintf('最小年费用为：%.3f 元\n',allFee)
[A w] = Area_Water(t2);
fprintf('冷却器传热面积为：%.3f m^2\n',A)
fprintf('每小时冷却水用量为：%.1f kg/h\n',w)
fee1 = JA*A*beta;
fee2 = Jw*theta*w/1000;
fprintf('年度投资费为：%.1f 元， 占总费用： %.2f %s\n', fee1,fee1/allFee*100,'%')
fprintf('年度操作费为：%.1f 元， 占总费用： %.2f %s\n', fee2,fee2/allFee*100,'%')

% ------------------------------------------------------------------
function J = TotalFee(t2)
global T1 T2 G t1 JA beta K theta Jw Cw Cc Q
[A w] = Area_Water(t2);
J = JA*A*beta + Jw*theta*w/1000;

% ------------------------------------------------------------------
function [A,w] = Area_Water(t2)
global T1 T2 G t1 JA beta K theta Jw Cw Cc Q
var1 = T1-t2;
var2 = T2-t1;
dtm = (var1-var2)/log(var1/var2);
A = Q/(K*dtm);
w = Q/Cw/(t2-t1);