- 在线时间
- 7 小时
- 最后登录
- 2015-7-21
- 注册时间
- 2012-11-30
- 听众数
- 9
- 收听数
- 0
- 能力
- 0 分
- 体力
- 2 点
- 威望
- 0 点
- 阅读权限
- 10
- 积分
- 4
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 8
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 2
升级   80% TA的每日心情 | 郁闷 2014-12-20 22:54 |
---|
签到天数: 1 天 [LV.1]初来乍到
- 自我介绍
- 数模爱好者
 |
原来的模型:
![]() ;
![]() ;
![]()
Model:
SETS:
Stream/1..8/:F; !the number of streams;
Com/1..3/:r;!the number of components,1,2,3=B,E,W;
Dis/1..2/:N,Nmin,RR,Rmin,QC,QR,O;!number of units operation or the sets of parameters;
Con(Stream,Com):x;!the component concentration of each stream;
Vol(Dis,Com):a;! relative volatility;
ENDSETS
DATA:
a=2 2.5 1
2.1 1 0.4;
r=30.783 41.5 43.94;!latent heat of three components;
ENDDATA
min=N(1)+N(2)+QC(1)+QC(2)+QR(1)+QR(2);
F(1)=100000;
F(7)=20000;
F(8)=80000;
x(1,1)=0;x(1,2)=0.8;x(1,3)=0.2;x(2,1)=1;x(2,2)=0;x(2,3)=0;
x(8,2)=0.98;
x(7,3)=0.90;
@FOR(Com(I):F(1)*x(1,I)+F(3)*x(3,I)=F(4)*x(4,I));
@FOR(Com(I):F(5)*x(5,I)+F(2)*x(2,I)=F(6)*x(6,I));
@FOR(Com(I):F(4)*x(4,I)=F(5)*x(5,I)+F(7)*x(7,I));
@FOR(Com(I):F(6)*x(6,I)=F(3)*x(3,I)+F(8)*x(8,I));
F(1)+F(3)=F(4);
F(4)=F(5)+F(7);
F(5)+F(2)=F(6);
F(6)=F(3)+F(8);
@FOR(Stream(I) SUM(Com(J):x(I,J))=1);
Nmin(1)=@log(x(5,1)*x(7,3)/(x(5,3)*x(7,1)))/@log(a(1,1));
Nmin(2)=@log(x(3,1)*x(8,2)/(x(3,2)*x(8,1)))/@log(a(2,1));
@SUM(Com(I):a(1,I)*x(4,I)/(a(1,I)-O(1)))=0;
@SUM(Com(I):a(2,I)*x(6,I)/(a(2,I)-O(2)))=0;
Rmin(1)=@SUM(Com(I):a(1,I)*x(5,I)/(a(1,I)-O(1)))-1;
Rmin(2)=@SUM(Com(I):a(2,I)*x(3,I)/(a(2,I)-O(2)))-1;
RR(1)=1.3*Rmin(1);
RR(2)=1.3*Rmin(2);
@FOR(Dis(K) N(K)-Nmin(K))/(N(K)+1)=0.75*(1-@POW(((R(K)-Rmin(K))/(R(K)+1)),0.566)));
QC(1)=@SUM(Com(I):x(5,I)*r(I))*F(5)*(RR(1)+1);
QR(1)=@SUM(Com(I):x(7,I)*r(I))*F(5)*(RR(1)+1);
QC(2)=@SUM(Com(I):x(3,I)*r(I))*F(3)*(RR(2)+1);
QR(2)=@SUM(Com(I):x(8,I)*r(I))*F(3)*(RR(2)+1);
x(3,2)=0.8933-1.0724*x(3,1);
x(5,2)=0.8933-1.0724*x(5,1);
x(4,2)-0.8933+1.0724*x(4,1)<0;
x(4,2)-(17/12)*x(4,1)>0;
x(6,2)-0.8933+1.0724*x(6,1)>0;
x(6,2)-1+(30*x(6,1))/23<0;
x(3,1)>=0.36;
x(3,1)<=0.46;
x(5,1)<0.36;
x(5,1)>0;
END
|
zan
|