| 查看: 532 | 回复: 0 | ||
[求助]
求助！！求助！！走过路过的大神帮忙看一段程序（关于精馏塔仿真软件的）
|
|
帮忙解释一下。。。最好能根据这些代码导出编写这段代码所依据的数学公式 三、模型整合 仿真模型主要精馏塔模型和调节仪表模型两部分组成。 精馏塔模型包含精馏塔的物料衡算、精馏塔塔板温度分布。 精馏塔物料衡算公式如下： F=W+D F*Xf=W*Xw+D*Xd 我们采用的精馏塔系统采用塔釜蒸汽间接加热的方式，所以从塔釜到塔顶温度逐步升高。在仿真模型中，将温度分布简单的归纳为一线性的变化，由塔釜到塔顶温度成线性分布。 精馏塔数学模型以4秒钟一次的运算速度在系统中反复执行，其代码如下： if(enablemodule==1){ double pi = 3.141592653579; float rr; float R1; timeFlag=timeFlag+1; R1=OldData[4][1]/F604[1]; rr=float(0.97)/float(60.16)/float(OldData[1][1]); if(rr<0.93) rr=float(0.93); if(rr>0.9999) rr=float(0.9999); for(int i=1;i<=DataRectNumBer;i++) { DataRect->calvalue(); } OldData[1][4]=float(59.90)+float(2.98)/R1*float(0.27)+(T602[1]-float(81.78))*float(0.6); OldData[2][4]=float(63.6)+float(2.98)/R1*float(0.25)+(T602[1]-float(81.78))*float(0.8); OldData[10][4]=float(87.1)+(T602[1]-float(81.78))*float(0.8)-(F604[1]-18139)/18139*float(0.2); OldData[9][1]=F604[1]*float(0.76)/rr; OldData[8][1]=F604[1]*float(0.24/0.18)/850; OldData[3][4]=OldData[3][2]/60*float(18139.14); OldData[4][4]=OldData[4][2]/60*float(42135); OldData[5][4]=OldData[5][2]/50*float(81.78)+float(170)*(float(18139.14)-F604[1])/float(18139.14); OldData[7][4]=OldData[7][2]/60*float(1.2); OldData[3][1]=OldData[3][1]+float(0.3)*(OldData[3][4]-OldData[3][1]); OldData[4][1]=OldData[4][1]+float(0.05)*(OldData[4][4]-OldData[4][1]); OldData[5][1]=OldData[5][1]+float(0.02)*(OldData[5][4]-OldData[5][1]); OldData[7][1]=OldData[7][1]+float(0.05)*(OldData[7][4]-OldData[7][1]); OldData[1][1]=OldData[1][1]+float(0.05)*(OldData[1][4]-OldData[1][1]); OldData[2][1]=OldData[2][1]+float(0.05)*(OldData[2][4]-OldData[2][1]); OldData[10][1]=OldData[10][1]+float(0.05)*(OldData[10][4]-OldData[10][1]); OldData[1][5]=OldData[1][1]+float(sin((timeFlag+1)/0.6*sin(2*timeFlag+1)))*float(0.03); OldData[2][5]=OldData[2][1]+float(sin((timeFlag+2)/0.6*sin(3*timeFlag+1)))*float(0.03); OldData[5][5]=OldData[5][1]+float(sin((timeFlag+3)/0.6*sin(4*timeFlag+1)))*float(0.03); OldData[6][5]=OldData[6][1]+float(sin((timeFlag+4)/0.6*sin(5*timeFlag+1)))*float(1.2); OldData[3][5]=OldData[3][1]+float(sin((timeFlag+5)/0.6*sin(6*timeFlag+1)))*float(40); OldData[8][5]=OldData[8][1]+float(sin((timeFlag+6)/0.6*sin(7*timeFlag+1)))*float(2); OldData[4][5]=OldData[4][1]+float(sin((timeFlag+7)/0.6*sin(8*timeFlag+1)))*float(100); OldData[7][5]=OldData[7][1]+float(sin((timeFlag+8)/0.6*sin(9*timeFlag+1)))*float(0.03); OldData[9][5]=OldData[9][1]+float(sin((timeFlag+9)/0.6*sin(10*timeFlag+1)))*float(100); OldData[10][5]=OldData[10][1]+float(sin((timeFlag+10)/0.6))*float(0.03); for(i=1;i<=10;i++) { T602=T602[i+1]; F604=F604[i+1]; } T602[10]=OldData[5][5]; F604[10]=OldData[3][5]; 调节仪表模型模拟了真实调节仪表中的比例控制，其阀门的开度与测量值和调节值的距离成比例，代码如下： void CDataRect::calvalue() { if(m_CanBeChange==1 && m_IsAuto==1) { float x,y,k,xx,yy; //x=OldData[m_Index][3]; //y=OldData[m_Index][6]; k=NewData2[m_Index]; yy=OldData[m_Index][2]; //yy=float(x/y)*60; x=OldData[m_Index][1]; y=OldData[m_Index][3]; xx=float((y-x)/k); yy=xx*100+yy; if(yy>100) yy=100; if(yy<0) yy=0; OldData[m_Index][2]=yy; } } |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 精馏塔仿真操作中的可靠性评价研究.doc
2014-04-07 16:20:50, 815.5 K
» 猜你喜欢
科研孤儿太难了
已经有19人回复
-
能否退出参与的面上项目解除限项
已经有9人回复
-
明天应该可查了！？
已经有4人回复
-
什么时候开奖？
已经有11人回复
-
范进中举一文的中心思想
已经有6人回复
-
建议基金发布提前给出明确的时间点
已经有14人回复
-
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
找到一些相关的精华帖子，希望有用哦~
Origin软件里怎么画出下面这种图形，求大神帮忙啊
已经有5人回复
- 请大神帮忙看看这是什么植物？？？？？
已经有6人回复
- 新手求助~~求大神帮忙！！
已经有3人回复
- 请大神帮忙看一下这两张图片的组织都有什么？由于不是专业的所以看起来比较困难。
已经有10人回复
- 2015年考研求助，求大神帮忙
已经有15人回复
- 积分求助 求大神帮忙！
已经有7人回复
- chemistry an Asia Journal 投稿求助！投过的大神们帮帮忙啊
已经有5人回复
- 很简单的Fortran赋值出错啦，求大神们帮忙找bug...
已经有11人回复
- 哪位大神帮忙看一下MATLAB求解微分方程的问题？
已经有3人回复
- APL这样回复还有希望吗？请各位大神帮忙看看。
已经有9人回复
- 求大神帮忙解答！
已经有8人回复
- 求助一道物理化学考研真题 还请各位路过的大神们出手相助
已经有3人回复
- 哪位大神帮忙翻译一下啊，关于混凝土自修复的，万分感激........
已经有5人回复
- 热重-差热分析的图，哪个大神帮忙解释一下呗……
已经有6人回复
- 求助大神帮忙解答如下题目
已经有5人回复
- 修改图标问题，请各位大神帮忙！
已经有5人回复
- 求各位大神帮忙下一篇文章
已经有3人回复
- 走过路过的师兄师姐师弟师妹们帮帮忙啊
已经有10人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
40