Blok-sxemalar asosida oddiy gidravlik sistemalarni hisoblash uchun dastur tuzish Reja

Issiqlik almashinuvi jarayonlari modellari
Issiqlik almashinadigan apparatlarning matematik modellarinivatov ni qurishdaо, apparatda m ga ega bo'lgan bir xil issiqlik energiyasi oqimlarining barglari va turlarini aniqlash uchun strukturaviy tahlil oldindan amalgaеoshiriladi. Har bir oqim uchun matematik opvaSaniya oqim va issiqlik uzatish harakati bilan bog'liq bo'lgan sovutish suvi oqimidagi harorat o'zgarishini tavsiflovchi ifoda shaklida yoziladi. Oldindanvaoldindan taxminlar shakllantiriladi.
Agar sovutish suvi oqimining tuzilishie-aralashtirish modeliga moskeladigan bo'lsa, unda bu oqimning matematik tavsifi uchunсissiqlik uzatishni hisobga olgan holda tenglama (4.8) dan foydalanish mumkin:
, (4.8)
bu erda v ideal aralashtirish oqimining hajmi, m³;
 - sovutish suvi zichligi, kg/ m³;
– sovutgichning o'ziga xos issiqlik quvvati, J/ (kgK);
 - volumetrik oqim tezligi, m³/s;
F-issiqlik almashinuvi yuzasi, m²;
 – issiqlik uzatish koeffitsienti, W/(M²K);
(Т_T – t) – devor va sovutish suvi haroratining farqi, K;
T_VX-kirish oqim harorati, K;
t-vaqt, bilan.
Agar oqim strukturasi idealrepressiya modeliga moskeladigan bo'lsa, unda oqimning matematik tavsifi uchun issiqlik uzatishni hisobga olgan holda tenglama (4.8) ishlatilishi mumkin:
, (4.9)
bu erda s_b-oqim kesimining maydoni, m²;
L-ideal joy maydoni uzunligi, m;
l-0 dan L ga qadar o'zgaruvchan koordinatalar;
T =T (l, t) – issiqlik oqimining haroratiоvaqt o'tishi bilan mekansal koordinatalar orqali tarqatilishi funktsiyasi.
Odatda tenglamalarda (3.8) va (3.9) issiqlik uzatish koeffitsienti, sovutgichning zichligi va issiqlik постояннымиquvvati haroratning cheklangan intervalida doimiy ravishda olinadi. Bubem oqim tezligi doimiy qoladi, deb taxmin qilinadi.
Cоставим ba'zi turdagi issiqlikuzatish qurilmalarining matematik modellarini qoldiring.

Tenglama (3.2) dan ioqimdagi i - moddaning massa ulushini aniqlash mumkinго вещества в вы:

. (3.3)
Issiqlik balansi tenglamasi quyidagi shaklga ega:
, (3.4)
bu erda C_p-oqimlarning o'ziga xos issiqligi, ;
t-oqim harorati, ^оu bilan.
Keyin chiqish oqimining harorati nisbati bilan belgilanadi
. (3.5)
m oqimidagi i-moddaning o'ziga xos issiqlik sig'imining haroratga bog'liqligi harorat jfunktsiyasi sifatida ifodalanishi mumkin:
, (3.6)
bu erda a_i, b_i, C_i, D_i ampirik koeffitsientlar (ilova F, jadval). 4) har bir narsa uchun topilganеstva [10].
Теплоемкость Jго oqimining issiqlik quvvati qo'shimchalar qoidasi bo'yicha hisoblanadiаддитивности:
. (3.7)
Tenglama (3.5) ni hal qilish uchun oddiy ite radio usulidan foydalanishingiz mumkin[6-8]:
, (3.8)
bu erda k-iteratsiya raqami.
Hisobni tugatish sharti – .
Dastlabki yondashuv sifatida qabul qilinishi mumkin
. (3.9)
Shunday qilib, oqim holatining parametrlarini belgilash, yning harakatini mikserga kiritish, chiqish parametrlarini aniqlash mumkin.


Misol
Mikserga kirish uchunоdami 30 va 50 kg/s va 400 va 500 s temperaturalari bo'yicha massa rasx bilan ikkita oqim Сberiladi.
Tarkibi nоommaviy oqimlarning%. пp quyida keltirilgan.

	H2	CH4	bilan2N6	bilan3N8	n - С4N10	bilan5N12
i oqim	60	30	10	–	–	–
II oqim	–	–	–	40	50	10

Kirish oqimlarining tarkibiga qarab, mikserdan chiqadigan oqim oltita komponentni o'z ichiga oladi. Shuning составы как входных поuchun, oqimlari va chiqish oqimi haqida kiritish n tarkibi tillo, hajmi [1-6] bilan ta'riflash kerak va chiqish oqimi tarkibi hisoblash tsiklik PR zedura yordamida tashkilоqilish
Oqimlarning har birining issiqlik sig'imini hisoblash ham davriy protsedura yordamida amalga oshirilishi mumkin.
Chiqish oqimining issiqlik sig'imi bir funktsiya vaсharorat koma bilan, chunki, uning hisoblash yineleme teshik qismiеhisoblanadi va har bir yangi harorat yondashuv uchun takrorlanadiо, oldingi iterasyonda nur haqida n, va keyin bu qiymatiоtank issiqlik quyidagi harorat yondashuv bilan belgilanadiоth oqimi chiqish.
Точность Yineleme jarayonini tugatish uchun shartni belgilaydigan aniqlik haroratni o'lchash uchun ishlatiladigan qurilmaning xatosi bilan baholanadi. Qurilmalarning aniqligi sinfini hisobga olgan holda,еpromW quvurlarining haroratini o'lchash uchunленных трубопроводов, примем ishlatiladi, biz = 0,5 S ni olamiz.
Ushbu muammoni hal qilish algoritmining oqim sxemasi shakl. 3.2, Turbo Paskal tilida hisoblash dasturi ilova ilmiy-tadqiqotеinstitutida berilgan.

Shakl. 3.2. Mikserni hisoblash dasturining oqim sxemasi (boshlanishi)

1