- •Министерство образования Российской Федерации тамбовский государственный технический университет
- •Утверждаю Зав. Кафедрой
- •Задание на курсовой проект (работу) Студент Манякина в.С. Группа cзc-41
- •Тгту 280202.Пз.Кр
- •Манякина в.С., группа cзc-41 Тамбов 2010
- •Содержание
- •Введение
- •Описание процесса
- •Математическое описание модели
- •Этапы имитационного исследования модели
- •Исследование модели
- •Результаты анализа
- •Приложение
- •Список литературы
Приложение
Исходная математическая модель:
#include <conio.h>
#include <stdio.h>
#include <math.h>
#define n 5
int main()
{ int i;
double L=0.0001,G=0.001,k=5e-5,F=1,dz=0.1,nu=0.142,Tk=3600,h,dt,V,Vb,Vv,m,xr,J;
double x1[n+2],x0[n+2]={0,0,0,0,0,0,0};
double y1[n+2],y0[n+2]={0,0,0,0,0,0,0};
clrscr();
x0[n+1]=0;
y0[0]=0.3;
V=F*dz;
Vb=V*nu;
Vv=V*(1-nu);
dt=10;
h=Tk/dt;
for(m=0;m<h;m++)
{
for(i=1;i<n+1;i++)
{ xr= y0[i]*18;
J=k*(xr-x0[i]);
x1[i]=x0[i]+dt*(L*(x0[i+1]-x0[i])+J)/Vb;
y1[i]=y0[i]+dt*(G*(y0[i-1]-y0[i])-J)/Vv;
}
for(i=1;i<n+1;i++)
{
x0[i]=x1[i];
y0[i]=y1[i];
printf("\n%d %lf %lf",i,x0[i],y0[i]);
} printf("\n\n");
getch();
}return 0;
}
Имитационная модель:
#include <conio.h>
#include <stdio.h>
#include <math.h>
# include <stdlib.h>
#define n 5
double xn(double m, double s, double min, double max)
{
int i;
double xr,x,sum;
m: for(i=0, sum=0;i<6;i++)
sum+=(double)rand()/RAND_MAX;
x=sqrt(2)*s*(sum-3)+m;
if(x<min||x>max) goto m;
return x;
}
int main()
{ int i,N=10000,l,n1,n2;
double L=0.0001,G,k=5e-5,F=1,dz=0.1,nu=0.142,Tk=3600,h,dt,V,Vb,Vv,m,xr,J,P1;
double x1[n+2],x0[n+2]={0,0,0,0,0,0,0};
double y1[n+2],y0[n+2]={0,0,0,0,0,0,0};
clrscr(); randomize ();
n1=0; n2=0;
for (l=0;l<N;l++)
{
G=xn(0.001,0.0001,0.0007,0.0013);
x0[n+1]=0;
y0[0]=0.3;
V=F*dz;
Vb=V*nu;
Vv=V*(1-nu);
dt=10;
h=Tk/dt;
for(m=0;m<h;m++)
{
for(i=1;i<n+1;i++)
{ xr= y0[i]*18;
J=k*(xr-x0[i]);
x1[i]=x0[i]+dt*(L*(x0[i+1]-x0[i])+J)/Vb;
y1[i]=y0[i]+dt*(G*(y0[i-1]-y0[i])-J)/Vv;
}
for(i=1;i<n+1;i++)
{
x0[i]=x1[i];
y0[i]=y1[i];
}
}
//printf("\n g=%g x1=%g yn=%g ",G,x0[1],y0[n]);getch();
if(y0[n]>0.05) n1++;
}
P1=(double)n1/N;
printf ("\n P1=%lf", P1);
getch();
return 0;
}
Список литературы
1. Кафаров В.В., Глебов М.Б. Математическое моделирование основных процессов химических производств: Учеб. пособие для вузов. – М.: Высш. шк., 1991. – 400с.: ил.
2. Александров И.А. Ректификационные и абсорбционные аппараты. Методы расчета и основы конструирования. – М.: Химия, 1978.
3. Дытнерский Ю.И. Основные процессы и аппараты химической технологии.
4. Плановский А.Н., Николаев П.И. Процессы и аппараты химической и нефтехимической технологии. – М.: Химия, 1972. – 496с.