24小时热门版块排行榜    

查看: 550  |  回复: 0

shin123

新虫 (初入文坛)

[求助] 求助!!求助!!走过路过的大神帮忙看一段程序(关于精馏塔仿真软件的)

帮忙解释一下。。。最好能根据这些代码导出编写这段代码所依据的数学公式




三、模型整合

仿真模型主要精馏塔模型和调节仪表模型两部分组成。
精馏塔模型包含精馏塔的物料衡算、精馏塔塔板温度分布。
精馏塔物料衡算公式如下:
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

» 猜你喜欢

近则立身兴家,远则求知成人.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 shin123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:25 by mibUvS8DDCwf
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:12 by vZfe6xYu34yj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:57 by vZfe6xYu34yj
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 23:45 by vZfe6xYu34yj
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:33 by vZfe6xYu34yj
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:09 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 by vZfe6xYu34yj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:01 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:00 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:33 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
信息提示
请填处理意见