Overview
| Comment: | Initial revision |
|---|---|
| Downloads: | Tarball | ZIP archive | SQL archive |
| Timelines: | family | ancestors | descendants | both | origin/master | trunk |
| Files: | files | file ages | folders |
| SHA3-256: |
fd97c8d5275abad1bc1811ada397b3e0 |
| User & Date: | gawthrop@users.sourceforge.net on 1998-02-25 15:03:51 |
| Other Links: | branch diff | manifest | tags |
Context
|
1998-02-25
| ||
| 15:28:34 | Added state.m rep. check-in: e0f4289844 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
| 15:03:51 | Initial revision check-in: fd97c8d527 user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
|
1998-02-24
| ||
| 22:35:33 | Added -p (preserve) switch to cp in tidy mode. check-in: ff434a775d user: gawthrop@users.sourceforge.net tags: origin/master, trunk | |
Changes
Added mttroot/mtt/lib/comp/simple/ES_eqn.m version [862ca697d4].
> > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > | 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 |
function structure = ES_eqn(name,bond_number,bonds,direction,cr,args, ...
structure,eqnfile);
% ES_eqn - equations for ES component
%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %%%%% Model Transformation Tools %%%%%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Matlab function ES_eqn
% ES_eqn(name,bond_number,bonds,direction,cr,args, ...
% structure,eqnfile);
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Copyright (c) P.J. Gawthrop, 1997.
if nargin<6
eqnfile = 'stdout';
end;
% Check that there are exactly two bonds.
ports = 2;
if check_bonds(bonds,ports,'ES')==0
return
end
% There are 2 ports; extract the information
e_bond = bonds(1,:);
s_bond = bonds(2,:);
% The efforts are the same
if e_bond(1)==1 %effort imposed on the enthalpy bond
fprintf(eqnfile,'%s := %s + t_0;\n' , ...
varname(name,bond_number(2), 1), ...
varname(name,bond_number(1), 1) );
else %effort imposed on the entropy bond
fprintf(eqnfile,'%s := %s - t_0;\n' , ...
varname(name,bond_number(1), 1), ...
varname(name,bond_number(2), 1) );
end;
% The flows need to be converted
if e_bond(2)==-1 %flow imposed on the enthalpy bond
fprintf(eqnfile,'%s := %s/(%s + t_0);\n' , ...
varname(name,bond_number(2), -1), ...
varname(name,bond_number(1), -1), ...
varname(name,bond_number(2), 1) );
else % flow imposed on the entropy bond
fprintf(eqnfile,'%s := %s*(%s - t_0);\n' , ...
varname(name,bond_number(1), -1), ...
varname(name,bond_number(2), -1), ...
varname(name,bond_number(2), 1) );
end;
|