24小时热门版块排行榜

返回列表

shin123

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 418.5
帖子: 11
在线: 15.2小时
虫号: 2422745
注册: 2013-04-14
性别: GG
专业: 能源化工

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

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

三、模型整合

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