го потоку=length (Xsc);% число компонентів холодного потоку=[a1 b1 c1 d1; a2 b2 c2 d2; ...; an bn cn dn];% матриця коеф. полінома для холодного потоку=[a1 b1 c1 d1; a2 b2 c2 d2; ...; an bn cn dn];% матриця коеф. полінома для гарячого потоку=1/((1/ah) + (1/ac) + (thw/tw))% коеф.теплопередачіi=1: kc% вважаємо ізобарну теплоємність холодного потоку при Тес
Cpc=Cpc +(koefc(Xsc(1,i),1)+koefc(Xsc(1,i),2)*Tec+koefc(Xsc(1,i),3)*Tec^2+koefc(Xsc(1,i),4)*Tec^3)*Xsc(2,i);%koefc=[a b c d]
end=Gsc * Cpc * Tsc;% теплота холодного потоку на вході
for i=1: kh=Cph +(koefh(Xsh(1,i),1)+koefh(Xsh(1,i),2)*Tehs+koefh(Xsh(1,i),3)*Tehs^2+koefh(Xsh(1,i),4)*Tehs^3)*Xsh(2,i);
end=Gsc * Cpc * Tsh;% теплота гарячого потоку на вході=Tsc + 0.01;=1; Qtpf=0;% для входу в ціклabs (Qtps-Qtpf) gt;=eps=Tec + 0.1;% задаємося температурою холодного потоку на виході=0; i=1: kc% вважаємо ізобарну теплоємність холодного потоку при Тес
Cpc=Cpc +(koefc(Xsc(1,i),1)+koefc(Xsc(1,i),2)*Tec+koefc(Xsc(1,i),3)*Tec^2+koefc(Xsc(1,i),4)*Tec^3)*Xsc(2,i);%koefc=[a b c d]
end=Gsc * Cpc * Tec;% визначаємо теплоту холодного потоку на виході=Qec-Qsc;% визначаємо теплосодержание=Qsh-Qtps;% визначаємо теплоту гарячого потоку на виході=1; Tehs=Tsh; % входимо в цикл і обнуляем Tehsabs (Tehs-Tehf) gt;=0.1% визначаємо температуру гарячого потоку на виході
Tehs=Tehs - 0.1;
Cph=0;
for i=1: kh
Cph=Cph +(koefh(Xsh(1,i),1)+koefh(Xsh(1,i),2)*Tehs+koefh(Xsh(1,i),3)*Tehs^2+koefh(Xsh(1,i),4)*Tehs^3)*Xsh(2,i);
end
Tehf=Qeh/(Gsh * Cph); (tsh-tec) gt; (teh-tsc)
tb=tsh-tec;% більша різниця температур
tm=teh-tsc;% менша різниця температур
tb=teh-tsc;% більша різниця температур
tm=tsh-tec;% менша різниця температур=(tb-tm)/(log (tb/tm));% рушійна сила теплопередачі
Qtpf=F * tcp * Ktp; Tec gt; Tehf;
Перевірка адекватності.
Модель адекватна об'єкту, якщо результати моделювання підтверджуються і можуть служити основою для прогнозування процесів, що протікають в досліджуваних об'єктах.
При цьому адекватність моделі залежить від мети моделювання і прийнятих критеріїв.
Відповідність вихідних координат, що визначаються за моделі, вихідним координатам, знайдених експериментально, при однакових значеннях вхідних координат, оцінюється дисперсією апроксимації:
де порядок моделі; число точок знятих в ході експерименту.
Для оцінки точності даних знімаються з об'єкта в ході експерименту користуються дисперсією відтворюваності:
де r - число ступенів свободи.
Для оцінки моделі на адекватність використовують критерій Фішера: знаходять відношення, по відомим значенням і r з таблиць Фішера знаходять. Якщо, то можна говорити про адекватність моделі, якщо - модель неадекватна.
При неадекватності об'єкту моделювання модель доповнюють і знову перевіряють на адекватність.
Висновок
теплообмінник кожухотрубний математичний
В ході курсової роботи була розроблена стаціонарна, детермінована математична модель із зосередженими параметрами Одноходовий кожухотрубного противоточного теплообмінника - підігрівача з наступними припущеннями:
· Тип теплообмінника - одноходовой кожухотрубний;
· Режим роботи теплообмінника - безперервний;
· Режим теплообміну - стаціонарний;
· Зміна агрегатного стану речовин при теплопередачі відсутня;
· Схема руху потоків - протівоточная;
· Втрати теплоти відсутні;
· Коефіцієнти тепловіддачі «холодного» і «гарячого» потоків визначаються при початковій температурі теплоносіїв.
Дана математична модель була реалізована на М-мові системи Matlab.
Список літератури
1.Саулін Д.В. Математичне моделювання ХТС (конспект лекцій).- Перм .: ПГТУ, 2003. - 91 с.
2.Касаткін А.Г. Основні процеси та апарати хімічної технології.- М .: Госхіміздат, 1960.- 832 с.
.Технічні термодинаміка: Робоча програма, завдання на контрольні роботи, методичні вказівки до виконання контрольних і практичних робіт.- СПб .: СЗТУ, 2004. - 139 с.
4.Курс фізичної хімії, том 1. - М .: Хімія, 1964. - 624 с.
.Учебно-методичний комплекс дисципліни студента «Моделювання та ідентифікація об'єктів...
