شبیه سازی برج تقطیر توسط matlab

حــامد · Dec 24, 2009

Binary Distillation Column with 30 trays (cyclohexane n-heptane)

Distillation column models are generally good test models for nonlinear model reduction and identification. The concentrations at each stage or tray are highly correlated. The dynamics of the distillation process can be described by a relatively few number of underlying dynamic states. A couple papers have been published with this model as an example application. One in particular is:

Hahn, J. and T.F. Edgar, An improved method for nonlinear model reduction using balancing of empirical gramians, Computers and Chemical Engineering, 26, pp. 1379-1397, (2002)

This plot shows the system response after a step change in the reflux ratio from 3.0 to 1.5. Each trajectory represents the mole fraction of cyclohexane at each tray. The top reflux material becomes less pure (more n-heptane) due to the increased draw from the top of the column.

Distillation 1 Distillation Column (Constant Relative Volatility)

Distillation 2 Distillation Column with Wilson Eq and Psat/P constant

Distillation 3 Distillation Column with Wilson Equation

Distillation 4 Distillation Column with Enthalpy Equation

همه فایلها در یک فایل زیپ

فايل ضميمهhttp://www.www.www.iran-eng.ir/images/attach/rar.gifdistillation.rar (62.8 كيلو بايت

ghazal · Apr 18, 2010

سلام فایل زیپ رو چه جور می تونم اجرا کنم؟

حــامد · Jul 19, 2010

کد:

function [ret,x0,str,ts,xts]=colas_PI(t,x,u,flag);
%COLAS_PI	is the M-file description of the SIMULINK system named COLAS_PI.
%	The block-diagram can be displayed by typing: COLAS_PI.
%
%	SYS=COLAS_PI(T,X,U,FLAG) returns depending on FLAG certain
%	system values given time point, T, current state vector, X,
%	and input vector, U.
%	FLAG is used to indicate the type of output to be returned in SYS.
%
%	Setting FLAG=1 causes COLAS_PI to return state derivatives, FLAG=2
%	discrete states, FLAG=3 system outputs and FLAG=4 next sample
%	time. For more information and other options see SFUNC.
%
%	Calling COLAS_PI with a FLAG of zero:
%	[SIZES]=COLAS_PI([],[],[],0),  returns a vector, SIZES, which
%	contains the sizes of the state vector and other parameters.
%		SIZES(1) number of states
%		SIZES(2) number of discrete states
%		SIZES(3) number of outputs
%		SIZES(4) number of inputs
%		SIZES(5) number of roots (currently unsupported)
%		SIZES(6) direct feedthrough flag
%		SIZES(7) number of sample times
%
%	For the definition of other parameters in SIZES, see SFUNC.
%	See also, TRIM, LINMOD, LINSIM, EULER, RK23, RK45, ADAMS, GEAR.

% Note: This M-file is only used for saving graphical information;
%       after the model is loaded into memory an internal model
%       representation is used.

% the system will take on the name of this mfile:
sys = mfilename;
new_system(sys)
simver(1.3)
if (0 == (nargin + nargout))
     set_param(sys,'Location',[59,72,853,623])
     open_system(sys)
end;
set_param(sys,'algorithm',     'Gear')
set_param(sys,'Start time',    '0.0')
set_param(sys,'Stop time',     '200')
set_param(sys,'Min step size', '0.01')
set_param(sys,'Max step size', '5')
set_param(sys,'Relative error','1e-6')
set_param(sys,'Return vars',   '')
set_param(sys,'AssignWideVectorLines','on');

add_block('built-in/Clock',[sys,'/','Clock'])
set_param([sys,'/','Clock'],...
		'position',[85,10,105,30])

add_block('built-in/To Workspace',[sys,'/','time'])
set_param([sys,'/','time'],...
		'mat-name','t',...
		'position',[485,12,535,28])

add_block('built-in/Constant',[sys,'/','F'])
set_param([sys,'/','F'],...
		'orientation',1,...
		'position',[170,65,200,85])

add_block('built-in/Constant',[sys,'/','zF'])
set_param([sys,'/','zF'],...
		'orientation',1,...
		'Value','0.5',...
		'position',[116,65,144,85])

add_block('built-in/Demux',[sys,'/','Demux'])
set_param([sys,'/','Demux'],...
		'outputs','[1,1,1,1,41]',...
		'position',[380,116,435,184])

add_block('built-in/To Workspace',[sys,'/','Comp.'])
set_param([sys,'/','Comp.'],...
		'mat-name','Comp',...
		'buffer','1200',...
		'position',[495,237,545,253])

add_block('built-in/To Workspace',[sys,'/','x_B'])
set_param([sys,'/','x_B'],...
		'mat-name','y2',...
		'buffer','1200',...
		'position',[495,142,545,158])

add_block('built-in/To Workspace',[sys,'/','M_D'])
set_param([sys,'/','M_D'],...
		'mat-name','y3',...
		'buffer','1200',...
		'position',[495,172,545,188])

add_block('built-in/To Workspace',[sys,'/','M_B'])
set_param([sys,'/','M_B'],...
		'mat-name','y4',...
		'buffer','1200',...
		'position',[495,202,545,218])

add_block('built-in/Constant',[sys,'/','rMD'])
set_param([sys,'/','rMD'],...
		'orientation',2,...
		'Value','0.5',...
		'position',[490,300,550,320])

add_block('built-in/Gain',[sys,'/','Gain'])
set_param([sys,'/','Gain'],...
		'orientation',2,...
		'Gain','-10',...
		'position',[330,290,360,320])

add_block('built-in/Sum',[sys,'/','Sum'])
set_param([sys,'/','Sum'],...
		'orientation',2,...
		'inputs','-+',...
		'position',[400,295,420,315])

add_block('built-in/Constant',[sys,'/','rMB'])
set_param([sys,'/','rMB'],...
		'orientation',2,...
		'Value','0.5',...
		'position',[490,350,550,370])

add_block('built-in/Gain',[sys,'/','Gain1'])
set_param([sys,'/','Gain1'],...
		'orientation',2,...
		'Gain','-10',...
		'position',[335,340,365,370])

add_block('built-in/S-Function',[sys,'/',['Distillation',13,'column',13,'(nonlinear)',13,'',13,'LV configuration']])
set_param([sys,'/',['Distillation',13,'column',13,'(nonlinear)',13,'',13,'LV configuration']],...
		'function name','colas',...
		'position',[285,117,355,183])

add_block('built-in/Mux',[sys,'/','Mux'])
set_param([sys,'/','Mux'],...
		'inputs','7',...
		'position',[225,116,265,184])

add_block('built-in/Constant',[sys,'/','qF'])
set_param([sys,'/','qF'],...
		'orientation',1,...
		'position',[220,65,250,85])

add_block('built-in/Sum',[sys,'/','Sum1'])
set_param([sys,'/','Sum1'],...
		'orientation',2,...
		'inputs','-+',...
		'position',[405,345,425,365])

add_block('built-in/Transfer Fcn',[sys,'/',' CyD'])
set_param([sys,'/',' CyD'],...
		'orientation',2,...
		'Numerator','[3.76 1]',...
		'Denominator','[3.76 0]',...
		'position',[340,411,390,449])

add_block('built-in/Sum',[sys,'/','Sum2'])
set_param([sys,'/','Sum2'],...
		'orientation',2,...
		'inputs','-+',...
		'position',[415,420,435,440])

add_block('built-in/To Workspace',[sys,'/','Delta L'])
set_param([sys,'/','Delta L'],...
		'mat-name','u1',...
		'buffer','1200',...
		'position',[500,477,550,493])

add_block('built-in/To Workspace',[sys,'/','Delta V'])
set_param([sys,'/','Delta V'],...
		'mat-name','u2',...
		'buffer','1200',...
		'position',[505,517,555,533])

add_block('built-in/Sum',[sys,'/','Sum4'])
set_param([sys,'/','Sum4'],...
		'position',[95,110,115,130])

add_block('built-in/Sum',[sys,'/','Sum5'])
set_param([sys,'/','Sum5'],...
		'orientation',3,...
		'position',[100,200,120,220])

add_block('built-in/Sum',[sys,'/','Sum3'])
set_param([sys,'/','Sum3'],...
		'orientation',2,...
		'inputs','-+',...
		'position',[415,385,435,405])

add_block('built-in/Gain',[sys,'/','KpyD'])
set_param([sys,'/','KpyD'],...
		'orientation',2,...
		'Gain','26.1',...
		'position',[250,414,295,446])

add_block('built-in/Transfer Fcn',[sys,'/','CxB'])
set_param([sys,'/','CxB'],...
		'orientation',2,...
		'Numerator','[3.31 1]',...
		'Denominator','[3.31 0]',...
		'position',[280,377,335,413])

add_block('built-in/Gain',[sys,'/','KpxB'])
set_param([sys,'/','KpxB'],...
		'orientation',2,...
		'Gain','-37.5',...
		'position',[215,378,265,412])

add_block('built-in/To Workspace',[sys,'/','y_D'])
set_param([sys,'/','y_D'],...
		'mat-name','y1',...
		'buffer','1200',...
		'position',[495,92,545,108])

add_block('built-in/Transport Delay',[sys,'/',['Transport',13,'Delay']])
set_param([sys,'/',['Transport',13,'Delay']],...
		'orientation',1,...
		'Delay Time','1',...
		'Initial Input','0.01',...
		'position',[550,265,590,295])


%     Subsystem  'Graph1'.

new_system([sys,'/','Graph1'])
set_param([sys,'/','Graph1'],'Location',[0,59,274,252])

add_block('built-in/Inport',[sys,'/','Graph1/x'])
set_param([sys,'/','Graph1/x'],...
		'position',[65,55,85,75])

add_block('built-in/S-Function',[sys,'/',['Graph1/S-function',13,'M-file which plots',13,'lines',13,'']])
set_param([sys,'/',['Graph1/S-function',13,'M-file which plots',13,'lines',13,'']],...
		'function name','sfuny',...
		'parameters','ax, color,dt',...
		'position',[130,55,180,75])
add_line([sys,'/','Graph1'],[90,65;125,65])
set_param([sys,'/','Graph1'],...
		'Mask Display','plot(0,0,100,100,[90,10,10,10,90,90,10],[65,65,90,40,40,90,90],[90,78,69,54,40,31,25,10],[77,60,48,46,56,75,81,84])',...
		'Mask Type','Graph scope.')
set_param([sys,'/','Graph1'],...
		'Mask Dialogue','Graph scope using MATLAB graph window.\nEnter plotting ranges and line type.|Time range:|y-min:|y-max:|Line type (rgbw-:*). Seperate each plot by ''/'':')
set_param([sys,'/','Graph1'],...
		'Mask Translate','color = @4; ax = [0, @1, @2, @3]; dt = -1;')
set_param([sys,'/','Graph1'],...
		'Mask Help','This block plots to the MATLAB graph window and can be used as an improved version of the Scope block. Look at the m-file sfuny.m to see how it works. This block can take scalar or vector input signal.')
set_param([sys,'/','Graph1'],...
		'Mask Entries','20\/0\/0.02\/''y-/g--/c-./w:/m*/ro/b+''\/')


%     Finished composite block 'Graph1'.

set_param([sys,'/','Graph1'],...
		'position',[685,116,715,154])


%     Subsystem  'Graph'.

new_system([sys,'/','Graph'])
set_param([sys,'/','Graph'],'Location',[0,59,274,252])

add_block('built-in/Inport',[sys,'/','Graph/x'])
set_param([sys,'/','Graph/x'],...
		'position',[65,55,85,75])

add_block('built-in/S-Function',[sys,'/',['Graph/S-function',13,'M-file which plots',13,'lines',13,'']])
set_param([sys,'/',['Graph/S-function',13,'M-file which plots',13,'lines',13,'']],...
		'function name','sfuny',...
		'parameters','ax, color,dt',...
		'position',[130,55,180,75])
add_line([sys,'/','Graph'],[90,65;125,65])
set_param([sys,'/','Graph'],...
		'Mask Display','plot(0,0,100,100,[90,10,10,10,90,90,10],[65,65,90,40,40,90,90],[90,78,69,54,40,31,25,10],[77,60,48,46,56,75,81,84])',...
		'Mask Type','Graph scope.')
set_param([sys,'/','Graph'],...
		'Mask Dialogue','Graph scope using MATLAB graph window.\nEnter plotting ranges and line type.|Time range:|y-min:|y-max:|Line type (rgbw-:*). Seperate each plot by ''/'':')
set_param([sys,'/','Graph'],...
		'Mask Translate','color = @4; ax = [0, @1, @2, @3]; dt = -1;')
set_param([sys,'/','Graph'],...
		'Mask Help','This block plots to the MATLAB graph window and can be used as an improved version of the Scope block. Look at the m-file sfuny.m to see how it works. This block can take scalar or vector input signal.')
set_param([sys,'/','Graph'],...
		'Mask Entries','20\/0.985\/1\/''y-/g--/c-./w:/m*/ro/b+''\/')


%     Finished composite block 'Graph'.

set_param([sys,'/','Graph'],...
		'position',[685,56,715,94])

add_block('built-in/Transport Delay',[sys,'/',['Transport',13,'Delay1']])
set_param([sys,'/',['Transport',13,'Delay1']],...
		'orientation',1,...
		'Delay Time','1',...
		'Initial Input','0.99',...
		'position',[610,220,650,250])

add_block('built-in/Constant',[sys,'/','xBs'])
set_param([sys,'/','xBs'],...
		'orientation',2,...
		'Value','0.01',...
		'position',[690,390,730,410])

add_block('built-in/Constant',[sys,'/','yDs'])
set_param([sys,'/','yDs'],...
		'orientation',2,...
		'Value','0.99',...
		'position',[690,424,730,446])

add_block('built-in/Constant',[sys,'/','V0'])
set_param([sys,'/','V0'],...
		'Value','3.20629',...
		'position',[0,240,60,260])

add_block('built-in/Constant',[sys,'/','L0'])
set_param([sys,'/','L0'],...
		'Value','2.70629',...
		'position',[0,103,55,127])
add_line(sys,[440,120;480,120;490,100])
add_line(sys,[440,150;460,150;460,180;490,180])
add_line(sys,[440,165;450,165;450,210;490,210])
add_line(sys,[440,180;439,180;439,245;490,245])
add_line(sys,[185,90;185,100;165,100;165,160;220,160])
add_line(sys,[130,90;130,100;155,100;155,170;220,170])
add_line(sys,[110,20;480,20])
add_line(sys,[360,150;375,150])
add_line(sys,[440,135;475,135;475,150;490,150])
add_line(sys,[270,150;280,150])
add_line(sys,[395,305;365,305])
add_line(sys,[485,310;425,310])
add_line(sys,[325,305;175,305;175,140;220,140])
add_line(sys,[440,150;460,150;460,300;425,300])
add_line(sys,[440,165;450,165;450,350;430,350])
add_line(sys,[485,360;430,360])
add_line(sys,[400,355;370,355])
add_line(sys,[330,355;185,355;185,150;220,150])
add_line(sys,[235,90;235,104;200,104;200,180;220,180])
add_line(sys,[410,430;395,430])
add_line(sys,[440,120;460,120;460,75;680,75])
add_line(sys,[120,120;220,120])
add_line(sys,[60,115;90,115])
add_line(sys,[65,250;105,250;105,225])
add_line(sys,[110,195;110,180;125,180;125,130;220,130])
add_line(sys,[410,395;340,395])
add_line(sys,[275,395;270,395])
add_line(sys,[335,430;300,430])
add_line(sys,[245,430;80,430;90,125])
add_line(sys,[210,395;115,395;115,225])
add_line(sys,[245,430;235,430;235,485;495,485])
add_line(sys,[210,395;200,395;200,525;500,525])
add_line(sys,[440,135;680,135])
add_line(sys,[570,300;570,385;470,385;470,390;440,390])
add_line(sys,[440,135;570,135;570,260])
add_line(sys,[440,120;630,120;630,215])
add_line(sys,[630,255;630,420;470,420;470,425;440,425])
add_line(sys,[685,400;440,400])
add_line(sys,[685,435;440,435])

drawnow

% Return any arguments.
if (nargin | nargout)
	% Must use feval here to access system in memory
	if (nargin > 3)
		if (flag == 0)
			eval(['[ret,x0,str,ts,xts]=',sys,'(t,x,u,flag);'])
		else
			eval(['ret =', sys,'(t,x,u,flag);'])
		end
	else
		[ret,x0,str,ts,xts] = feval(sys);
	end
else
	drawnow % Flash up the model and execute load callback
end

Thread starter	عنوان	تالار	پاسخ ها	تاریخ
S	شبیه سازی برج تقطیر واکنشی با aspen pluse11.1	مدلسازی و شبیه سازی	1	Jun 12, 2013
	طراحی و شبیه سازی برج تقطیر	مدلسازی و شبیه سازی	91	Dec 11, 2010
	شبیه سازی برج تقطیر با matlab	مدلسازی و شبیه سازی	1	Feb 8, 2010
A	شبیه سازی برج تقطیر با matlab	مدلسازی و شبیه سازی	19	Nov 14, 2009
A	شبیه سازی با نرم افزار متریال استودیو(material studio)	مدلسازی و شبیه سازی	0	May 23, 2019

شبیه سازی برج تقطیر توسط matlab

حــامد

مدیر بازنشسته

ghazal

عضو جدید

حــامد

مدیر بازنشسته

Similar threads