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

Курсова робота

з дисципліни: «Основи інформатики та обчислювальної техніки»

на тему «Оптимізація параметрів елементів системи теплопостачання»

Одеса 2010

Вступ

Дана курсова робота – це поєднання усіх навиків які були отримані студентом на протязі навчального курсу програми.

Курсова робота передбачає в собі розрахунок оптимального діаметру енергосистеми для найменших затрат.

В роботі використовується поєднання як текстових так і графічних редакторів, а також однієї з вивчених алгоритмічних мов програмування.

Для перевірки правильності виконаних розрахунків уся робота виконується двома способами, збіжність котрих і буде свідченням правильності виконання. Серед таких методів в курсовій роботі використаний пакет Microsoft Excel та для перевірки мова програмування Basic, алгоритм програми якої був написаний за допомогою програми QBasic.

1. Основна мета курсової роботи

Основною метою курсової роботи є визначення оптимального діаметра теплової dопт теплової мережі системи теплопостачання від джерела до споживача при змінені швидкості теплоносія у внутрішньому контурі (у діапазоні w₁=0.4÷5.0 м/с) та діаметра теплообмінника зовнішнього контуру d2 (у діапазоні 0,25÷0,4м).

При зростанні швидкості теплоносія w₁ у трубопроводі відповідно зменшується його діаметр d₁ та товщина стінки δ₁, що зменшує витрати на матеріал трубопроводу. Однак відповідно зростає гідравлічний опір у системі та, отже, витрати електроенергії на привод мережного насоса. Тому завданням КР є визначення оптимального значення швидкості теплоносія w₁, при якій щорічні витрати на енергосистему будуть мінімальні.

2. Текст завдання курсової роботи з розрахунковими формулами

Визначають діаметр внутрішнього контуру енергосистеми:

де G1 - витрати теплоносія у внутрішньому контуру, кг/с;

π- коефіцієнт, що дорівнює 3,14;

ρ₁- густина теплоносія внутрішнього контору, кг/м³;

с - питома теплоємність теплоносія, Дж/(кг*К).

Визначають товщину стінки трубопроводів внутрішнього та зовнішнього контурів енергосистеми:

,

де p_1, p₂ – відповідно тиск всередині трубопроводу внутрішнього та зовнішнього контурів, МПа;

[σ_м] - припустима міцність матеріалу трубопроводу, МПа.

d2 – діаметр зовнішнього контору, що змінюється в діапазоні 0,25÷0,4м.

Число Рейнольда для внутрішнього та зовнішнього контурів:

;

де w₂,w₂ – відповідно середня кінематична в’язкість теплоносіїв внутрішнього та зовнішнього контурів при заданій температурі, м²/с;

d_2екв=d₂-d₁-₁ - еквівалентний діаметр кільцевого каналу теплообмінника, м;

швидкість теплоносія у кільцевому каналі теплообмінника, м/с;

p₂ – густина теплоносія в зовнішнього контурі, кг/м³;

G₂- витрата теплоносія в зовнішньому контурі, кг/с;

Число Нусальта для внутрішнього та зовнішнього контурів:

; ,

Де Pr₁, Pr₂ – відповідно числа Прандтля теплоносіїв внутрішнього та зовнішнього контурів при заданій температурі.

Коефіцієнти тепловіддачі теплоносіїв внутрішнього та зовнішнього контурів, Вт/(м² К):

;

Де , - відповідно коефіцієнти теплопровідності носіїв внутрішнього та зовнішнього контурів при заданій температурі, Вт/(м К).

Коефіцієнти теплопередачі від внутрішнього до зовнішнього контурів:

,

Де - коефіцієнт теплопровідності матеріалу трубопроводів, Вт/(м К).

Розрахункова площа теплообмінника:

,

Де Q-теплова продуктивність енергосистеми, кВт;

- температурний напір, К.

Довжина теплообмінника типу «труба в трубі»:

.

Маса матеріалу енергосистеми:

Де - відповідно довжина ділянок енергосистеми, м;

- густина матеріалу трубопроводів, кг/м³

Витрата коштів на енергосистему за рік:

,

Де Е_{н -} нормативний коефіцієнт, що дорівнює 0,15:

S_м- вартість металу трубопроводу, грн/кг;

S_ел- вартість 1 кВт∙год електроенергії;

λ- коефіцієнт тертя трубопроводу;

∑- підсумковий коефіцієнт місцевого опору трубопроводу;

- ККД мережного насоса;

- ККД електродвигуна для приводу насоса;

g – прискорення вільного падіння, м/с²

3. Електронна таблиця з оптимізаційними розрахунками енергосистеми

Таблиця 1

4.Текст програми на мові програмування QBASIK

Розрахунки виконуються за допомогою програми 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

z = m * en * sm + g1 * (1 + lyamda * (l + l1 + l2 + l3) / d1 + dzeta) * kpd1 * kpd2 * h * sel * .001 * w1 ^ 2 / 2 / 9.81

PRINT "w1="; w1; TAB(15); "d1="; d1; TAB(30); "delta1="; delta1; TAB(45); "l="; l; TAB(60); "z="; z

NEXT: NEXT

END

5. Результати розрахунків у QBASIC

Висновки

1. Під час виконання даної курсової роботи був розрахований оптимальний діаметр теплової мережі системи теплопостачання від джерела до споживача при змінній швидкості теплоносію та витратах на електроенергію.

2. Розрахунки були виконані двома способами: у середовищі Microsoft Excel та за допомогою алгоритмічної мови програмування Quick Basic. При виконанні завдання були отримані навички роботи з електронними таблицями та з графічним відображенням інформації. При написанні програми до курсової роботи згадали прийоми та методи алгоритмічного програмування.

3. До курсової роботи було зроблено креслення теплової мережі.

4. Розрахунки параметрів системи теплопостачання, виконані двома способами збігаються.

5. Основною метою курсової роботи було визначення оптимального діаметра d_опт теплової мережі теплопостачання від джерела до споживача при зміні швидкості теплоносія у внутрішньому контурі (у діапазоні = 0,4÷0,5 м/с) та діаметра теплообмінника зовнішнього контуру d₂ (у діапазоні = 0,25÷0,4 м) при даних вхідних параметрах.

При зростанні швидкості теплоносія w₁ у трубопроводі. Відповідно зменшується його діаметр d₁ та товщина стінки δ₁ , що зменшує витрати на матеріал трубопроводу. Однак відповідно зростає гідравлічний опір , отже , витрати енергії на привод мережного насоса

У Курсовій роботі отримані оптимальні значення швидкості теплоносія w₁ , при якій щорічні витрати на енергосистему будуть мінімальні .

Для d₁ = 0.25 мінімальне значення швидкості витрат на енергосистему

Z= 611,6454 грн. при w₁ = 1,4 м/с , де d₁ = 0,124372858 м , δ₁ =0,00954 м

Для d₂ = 0.3 мінімальне значення швидкості витрат на енергосистему

Z= 665,8462 грн. при w₁ = 1,4 м/с , де d₁ = 0,124372858 м , δ₁ =0,00954 м

Для d₃ = 0.35 мінімальне значення швидкості витрат на енергосистему

Z= 737,4314 грн. при w₁ = 1,4 м/с , де d₁ = 0,124372858 м , δ₁ =0,0954 м

Для d₄ = 0.4 мінімальне значення швидкості витрат на енергосистему

Z= 831,278 грн. при w₁ = 1,4 м/с , де d₁ = 0,12437286 м , δ₁ =0,0095 м

Перелік використаної літератури

1. Довженко В.А., Колєсніков Ю.В. MS Excel 2003. СПБ., БВХ-Петербург , 2004.

2. Ананьев А.В.,Федоров А.М.Самоучитель Visual Basic 6.0.

3. А.В.Кузін, С.В.Левонисова Бази даних: підр. посібник для студ. вищ. навч. закладів, 2008.

4. Клима И. Оптимизация энергетических систем. - М.: Высш. шк., 1991.-302 с.

5. Гринчишин Я.Т., Алгоритмы и программы на Бейсике. - М.; “Просвещение” 1988.

6. Бундюк А.М. Програмування в середовищі QBASIC:Конспект лекцій.- Одеса: ОДПУ,1996.–124с.

7. Фаронов В.В., Delphi. Програмування на мові високого рівня, 2005г.