Курсовая работа: Оптимізація параметрів елементів системи теплопостачання

Розрахунки виконуються за допомогою програми QBASIK

Програма оптимізації діаметра системи теплопостачання

Склав студент групи ТЕ –0901 Незгодюк Л.И.

Вхідні дані для розрахунку

rom = 8000: ro1 = 1000: ro2 = 1000: en = .15: dzeta = 30: lyamda = .02

kpd2 = .9: sel = .4: c = 4180: lyamdam = 40: pr1 = 2.2: pr2 = 3: lyamda1 = .67

lyamda2 = .66: mu1 = 3.6E-07: mu2 = 4.8E-07: g1 = 17: g2 = 32: p1 = 1.2: p2 = .7

sigma = 22: l1 = 22: l2 = 32: l3 = 42: q = 120: sm = 1.7: dt = 17: h = 5700:

kpd1 = .7

Тіло програми

DIM d1, delta1, l, z, w1, delta2

FOR d2 = .25 TO .4 STEP .05

PRINT USING "#.##"; d2

delta2 = .1 * (d2 * p2 / (2 * sigma / 2.5 - p2)) ^ .5

FOR w1 = .4 TO 1 STEP .2

d1 = (4 * g1 / 3.14 / ro1 / w1) ^ .5

delta1 = .1 * (d1 * p1 / (2 * sigma / 2.5 - p1)) ^ .5

re1 = w1 * d1 / mu1

nu1 = .021 * re1 ^ .8 * pr1 ^ .43

alfa1 = nu1 * lyamda1 / d1

d2ekv = d2 - d1 - delta1

w2 = 4 * g2 / (3.14 * ro2 * (d2 ^ 2 - (d1 + delta1) ^ 2))

re2 = w2 * d2ekv / mu2

nu2 = .021 * re2 ^ .8 * pr2 ^ .43

alfa2 = nu2 * lyamda2 / d2ekv

k = 1 / (1 / alfa1 + delta1 / lyamdam + 1 / alfa2)

fr = q * 1000 / k / dt

l = fr / 3.14 / (d1 + delta1)

m = 3.14 / 4 * (((d1 + delta1) ^ 2 - d1 ^ 2) * (l + l1 + l2 + l3) + ((d2 + delta2) ^ 2 - d2 ^ 2) * l) * rom