Оптимізація параметрів елементів системи теплопостачання
Оптимізація параметрів елементів системи теплопостачання
Курсова робота з дисципліни: «Основи інформатики та обчислювальної техніки» на тему «Оптимізація параметрів елементів системи теплопостачання» Одеса 2010 Вступ Дана курсова робота - це поєднання усіх навиків які були отримані студентом на протязі навчального курсу програми. Курсова робота передбачає в собі розрахунок оптимального діаметру енергосистеми для найменших затрат. В роботі використовується поєднання як текстових так і графічних редакторів, а також однієї з вивчених алгоритмічних мов програмування. Для перевірки правильності виконаних розрахунків уся робота виконується двома способами, збіжність котрих і буде свідченням правильності виконання. Серед таких методів в курсовій роботі використаний пакет Microsoft Excel та для перевірки мова програмування Basic, алгоритм програми якої був написаний за допомогою програми QBasic. 1. Основна мета курсової роботи Основною метою курсової роботи є визначення оптимального діаметра теплової dопт теплової мережі системи теплопостачання від джерела до споживача при змінені швидкості теплоносія у внутрішньому контурі (у діапазоні w1=0.4ч5.0 м/с) та діаметра теплообмінника зовнішнього контуру d2 (у діапазоні 0,25ч0,4м). При зростанні швидкості теплоносія w1 у трубопроводі відповідно зменшується його діаметр d1 та товщина стінки д1, що зменшує витрати на матеріал трубопроводу. Однак відповідно зростає гідравлічний опір у системі та, отже, витрати електроенергії на привод мережного насоса. Тому завданням КР є визначення оптимального значення швидкості теплоносія w1, при якій щорічні витрати на енергосистему будуть мінімальні. 2. Текст завдання курсової роботи з розрахунковими формулами Визначають діаметр внутрішнього контуру енергосистеми:
де G1 - витрати теплоносія у внутрішньому контуру, кг/с; р- коефіцієнт, що дорівнює 3,14; с1- густина теплоносія внутрішнього контору, кг/м3; с - питома теплоємність теплоносія, Дж/(кг*К). Визначають товщину стінки трубопроводів внутрішнього та зовнішнього контурів енергосистеми: , де p1, p2 - відповідно тиск всередині трубопроводу внутрішнього та зовнішнього контурів, МПа; [ум] - припустима міцність матеріалу трубопроводу, МПа. d2 - діаметр зовнішнього контору, що змінюється в діапазоні 0,25ч0,4м. Число Рейнольда для внутрішнього та зовнішнього контурів: ; де w2,w2 - відповідно середня кінематична в'язкість теплоносіїв внутрішнього та зовнішнього контурів при заданій температурі, м2/с; d2екв=d2-d1-1 - еквівалентний діаметр кільцевого каналу теплообмінника, м;
швидкість теплоносія у кільцевому каналі теплообмінника, м/с; p2 - густина теплоносія в зовнішнього контурі, кг/м3; G2- витрата теплоносія в зовнішньому контурі, кг/с; Число Нусальта для внутрішнього та зовнішнього контурів: ; , Де Pr1, Pr2 - відповідно числа Прандтля теплоносіїв внутрішнього та зовнішнього контурів при заданій температурі. Коефіцієнти тепловіддачі теплоносіїв внутрішнього та зовнішнього контурів, Вт/(м2 К): ; Де , - відповідно коефіцієнти теплопровідності носіїв внутрішнього та зовнішнього контурів при заданій температурі, Вт/(м К). Коефіцієнти теплопередачі від внутрішнього до зовнішнього контурів: , Де - коефіцієнт теплопровідності матеріалу трубопроводів, Вт/(м К). Розрахункова площа теплообмінника: , Де Q-теплова продуктивність енергосистеми, кВт; - температурний напір, К. Довжина теплообмінника типу «труба в трубі»: . Маса матеріалу енергосистеми:
Де - відповідно довжина ділянок енергосистеми, м; - густина матеріалу трубопроводів, кг/м3 Витрата коштів на енергосистему за рік: , Де Ен - нормативний коефіцієнт, що дорівнює 0,15: Sм- вартість металу трубопроводу, грн/кг; Sел- вартість 1 кВт•год електроенергії; л- коефіцієнт тертя трубопроводу; ?- підсумковий коефіцієнт місцевого опору трубопроводу; - ККД мережного насоса; - ККД електродвигуна для приводу насоса; g - прискорення вільного падіння, м/с2 3. Електронна таблиця з оптимізаційними розрахунками енергосистеми Таблиця 1 |
Початкові дані для розрахунків | | | Витрата теплоносія у внутр. контурі, G1 кг/с | 29 | | Витрата теплоносія у зовн. контурі, G2 кг/с | 44 | | Густина теплоносія внутр. контуру, с1, кг/м3 | 1000 | | Густина теплоносія зовн. контуру, с2, кг/м3 | 1000 | | Густина матеріалу трубопроводів, см, кг/м3 | 8000 | | Нормативний коефіцієнт, Ен | 0,15 | | Тиск теплоносія внутрішнього контуру, Р1, Мпа | 2,4 | | Тиск теплоносія зовнішнього контуру, Р2. Мпа | 1,9 | | Міцність матер. трубопроводів, у, Мпа | 34 | | Довжина ділянки трубопроводу, L1, м | 34 | | Довжина ділянки трубопроводу, L2, м | 44 | | Довжина ділянки трубопроводу, L3, м | 54 | | Число годин роботи системи за рік, h, год | 6800 | | ККД мережного насоса, зн | 0,65 | | ККД електродвигуна, зэл | 0,9 | | Вартість електроенергії, Sел, грн/(кВт) | 0,4 | | Вартість металу трубопроводу, Sм, грн/кг | 2,8 | | Підсумковий коеф. місц. опору трубопроводу, Уом | 30 | | Коеф. Тертя трубопроводу, л | 0,02 | | Теплова продуктивність системи, Q, кВт | 240 | | Питома теплоємність теплоносіїв, с, Дж/(кг•К) | 4180 | | Температура на вих. теплообмінника, Т2вих, гр. С | 60 | | Температурний напір, ?t, гр. С | 29 | | Коеф. Теплопровідності матер. трубопроводу, лм, Вт/(м•К) | 40 | | Число Прандтля для внутр.. контуру, Pr1 | 2,2 | | Число Прандтля для зовн. Контуру,Pr2 | 3 | | Коеф. Теплопровідності теплонос. Внутр.. контуру, л1, Вт/(м•К) | 0,67 | | Коеф. Теплопровідності теплонос.зовн.. контуру, л2, Вт/(м•К) | 0,66 | | Кінемат. В'язкість теплонос. Внутр.. контуру, v1, м.кв/с | 3,60Е-07 | | Кінемат. В'язкість теплонос.зовн.. контуру, v2, м.кв/с | 4,80Е-07 | | |
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 |
| б2 | d2 | | | | | 0.0073 | 0.25 | | | | w1, м/с | d1 , м | б1 , м | L , м | Z,грн/год | | .4 | .23268031 | .01305 | 8.14001 | 1097.26 | | .6 | .18998268 | .01179 | 9.00539 | 848.947 | | .8 | .16452983 | .01097 | 9.76593 | 724.761 | | 1 | .14715995 | .01038 | 10.4481 | 657.525 | | 1.2 | .13433804 | .00991 | 11.0702 | 623.369 | | 1.4 | .12437286 | .00954 | 11.6447 | 611.645 | | 1.6 | .11634016 | .00923 | 12.1802 | 616.918 | | 1.8 | .10968655 | .00896 | 12.6834 | 636.157 | | 2 | .1040578 | .00873 | 13.1591 | 667.572 | | 2.2 | .09921522 | .00852 | 13.611 | 710.07 | | 2.4 | .09499134 | .00834 | 14.0422 | 762.967 | | 2.6 | .09126473 | .00817 | 14.455 | 825.841 | | 2.8 | .08794489 | .00802 | 14.8516 | 898.436 | | 3 | .08496284 | .00788 | 15.2335 | 980.609 | | 3.2 | .08226491 | .00776 | 15.6021 | 1072.29 | | 3.4 | .07980869 | .00764 | 15.9587 | 1173.48 | | 3.6 | .0775601 | .00753 | 16.3042 | 1284.19 | | 3.8 | .07549146 | .00743 | 16.6397 | 1404.49 | | 4 | .07357998 | .00734 | 16.9657 | 1534.46 | | 4.2 | .0718067 | .00725 | 17.2831 | 1674.18 | | 4.4 | .07015575 | .00716 | 17.5924 | 1823.76 | | 4.6 | .06861368 | .00709 | 17.8942 | 1983.32 | | 4.8 | .06716902 | .00701 | 18.189 | 2152.98 | | 5 | .06581193 | .00694 | 18.4771 | 2332.86 | | | б2 | d2 | | | | | 0.008 | 0.3 | | | | w1, м/с | d1 , м | б1 , м | L , м | Z,грн/год | | 0.4 | .32680312 | .013048 | 9.593718 | 1138.274 | | 0.6 | .189982679 | .01179 | 10.77816 | 892.6906 | | 0.8 | .164529826 | .010972 | 11.80706 | 771.2051 | | 1 | .14715995 | .010377 | 12.72505 | 706.5952 | | 1.2 | .134338041 | .009914 | 13.55975 | 675.0131 | | 1.4 | .124372858 | 0.00954 | 14.3292 | 665.8462 | | 1.6 | .16340156 | .009226 | 15.04585 | 673.6911 | | 1.8 | .09686551 | .008959 | 15.71869 | 695,5492 | | 2 | .104057799 | .008726 | 16.35445 | 729.6609 | | 2.2 | .099215218 | .00852 | 16.95832 | 774.9588 | | 2.4 | .094991339 | .008337 | 17.53438 | 830.7873 | | 2.6 | .091264727 | .008172 | 18.08594 | 896.7474 | | 2.8 | .087944891 | .008022 | 18.6157 | 972.6068 | | 3 | .084962837 | .007885 | 19.12589 | 1058.245 | | 3.2 | .082264913 | .007758 | 19.6184 | 1153.619 | | 3.4 | .079808688 | .007642 | 20.09485 | 1258.738 | | 3.6 | .077560104 | .007533 | 20.5566 | 1373.652 | | 3.8 | .075491462 | .007432 | 21.00486 | 1498.439 | | 4 | .073579975 | .007338 | 21.44066 | 1633.195 | | 4.2 | .071806703 | .007249 | 21.86493 | 1778.035 | | 4.4 | .070155754 | .007165 | 22.27846 | 1933.082 | | 4.6 | .068613681 | .007086 | 22.682 | 2098.471 | | 4.8 | .06716902 | .007011 | 23.07617 | 2274.341 | | 5 | .06581193 | .006939 | 23.46155 | 2460.836 | | | б2 | d2 | | | | | 0.0086 | 0.35 | | | | w1, м/с | d1 , м | б1 , м | L , м | Z,грн/год | | .4 | .232680312 | .013048 | 11.25716 | 1190.611 | | .6 | .189982679 | .01179 | 12.80816 | 948.84146 | | .8 | .164529826 | .010972 | 14.14511 | 831.47296 | | 1 | .14715995 | .010377 | 15.33363 | 777.6505 | | 1.2 | .134338041 | .009914 | 16.41215 | 742.9351 | | 1.4 | .124372858 | .00954 | 17.40517 | 737.43137 | | 1.6 | .116340156 | .009226 | 18.32934 | 748.90376 | | 1.8 | .109686551 | .008959 | 19.19659 | 774.39992 | | 2 | .104057799 | .008726 | 2.01578 | 812.20129 | | 2.2 | .099215218 | .00852 | 2.79372 | 861.27755 | | 2.4 | .094991339 | .008337 | 21.53574 | 921.0069 | | 2.6 | .091264727 | .008172 | 22.24615 | 991.02195 | | 2.8 | .087944891 | .008022 | 22.92845 | 1071.1202 | | 3 | .084962837 | .007885 | 23.58555 | 1161.2094 | | 3.2 | .082264913 | .007758 | 24.21991 | 1261.2731 | | 3.4 | .079808688 | .007642 | 24.83358 | 1371.3476 | | 3.6 | .077560104 | .007533 | 25.42836 | 1491.5074 | | 3.8 | .075491462 | .007432 | 26.00579 | 1621.8544 | | 4 | .073579975 | .007338 | 26.56722 | 1762.5105 | | 4.2 | .071806703 | .007249 | 27.11382 | 1913.6125 | | 4.4 | .070155754 | .007165 | 27.64663 | 2075.3082 | | 4.6 | .068613681 | .007086 | 28.1666 | 2247.7536 | | 4.8 | .06716902 | .007011 | 28.67455 | 2431.111 | | 5 | .06581193 | .006939 | 29.17121 | 2625.547 | | | б2 | d2 | | | | | 0.0092 | 0.4 | | | | w1, м/с | d1 , м | б1 , м | L , м | Z,грн/год | | .4 | .23268031 | .01305 | 13.1251 | 1255.3841 | | .6 | .18998268 | .01179 | 15.0886 | 1020.7086 | | .8 | .16452983 | .01097 | 16.7721 | 909.3363 | | 1 | .14715995 | .01038 | 18.2647 | 854.21246 | | 1.2 | .13433804 | .00991 | 19.6173 | 831.67772 | | 1.4 | .12437286 | .00954 | 20.8615 | 831.27846 | | 1.6 | .11634016 | .00923 | 22.0188 | 847.74343 | | 1.8 | .10968655 | .00896 | 23.1045 | 878.18555 | | 2 | .1040578 | .00873 | 24.1297 | 920.94196 | | 2.2 | .09921522 | .00852 | 25.1031 | 975.03117 | | 2.4 | .09499134 | .00834 | 26.0315 | 1039.8751 | | 2.6 | .09126473 | .00817 | 26.9203 | 1115.1463 | | 2.8 | .08794489 | .00802 | 27.7738 | 1200.6792 | | 3 | .08496284 | .00788 | 28.5959 | 1296.4163 | | 3.2 | .08226491 | .00776 | 29.3894 | 1402.3737 | | 3.4 | .07980869 | .00764 | 30.1572 | 1518.6195 | | 3.6 | .0775601 | .00753 | 30.9013 | 1645.2582 | | 3.8 | .07549146 | .00743 | 31.6237 | 1782.4207 | | 4 | .07357998 | .00734 | 32.3261 | 1930.2573 | | 4.2 | .0718067 | .00725 | 33.01 | 2088.9325 | | 4.4 | .07015575 | .00716 | 33.6767 | 2258.6211 | | 4.6 | .06861368 | .00709 | 34.3273 | 2439.5056 | | 4.8 | .06716902 | .00701 | 34.963 | 2631.7743 | | 5 | .06581193 | .00694 | 35.5845 | 2835.6195 | | |
Висновки 1. Під час виконання даної курсової роботи був розрахований оптимальний діаметр теплової мережі системи теплопостачання від джерела до споживача при змінній швидкості теплоносію та витратах на електроенергію. 2. Розрахунки були виконані двома способами: у середовищі Microsoft Excel та за допомогою алгоритмічної мови програмування Quick Basic. При виконанні завдання були отримані навички роботи з електронними таблицями та з графічним відображенням інформації. При написанні програми до курсової роботи згадали прийоми та методи алгоритмічного програмування. 3. До курсової роботи було зроблено креслення теплової мережі. 4. Розрахунки параметрів системи теплопостачання, виконані двома способами збігаються. 5. Основною метою курсової роботи було визначення оптимального діаметра dопт теплової мережі теплопостачання від джерела до споживача при зміні швидкості теплоносія у внутрішньому контурі (у діапазоні = 0,4ч0,5 м/с) та діаметра теплообмінника зовнішнього контуру d2 (у діапазоні = 0,25ч0,4 м) при даних вхідних параметрах. При зростанні швидкості теплоносія w1 у трубопроводі. Відповідно зменшується його діаметр d1 та товщина стінки д1 , що зменшує витрати на матеріал трубопроводу. Однак відповідно зростає гідравлічний опір , отже , витрати енергії на привод мережного насоса У Курсовій роботі отримані оптимальні значення швидкості теплоносія w1 , при якій щорічні витрати на енергосистему будуть мінімальні . Для d1 = 0.25 мінімальне значення швидкості витрат на енергосистему Z= 611,6454 грн. при w1 = 1,4 м/с , де d1 = 0,124372858 м , д1 = 0,00954 м Для d2 = 0.3 мінімальне значення швидкості витрат на енергосистему Z= 665,8462 грн. при w1 = 1,4 м/с , де d1 = 0,124372858 м , д1 = 0,00954 м Для d3 = 0.35 мінімальне значення швидкості витрат на енергосистему Z= 737,4314 грн. при w1 = 1,4 м/с , де d1 = 0,124372858 м , д1 = 0,0954 м Для d4 = 0.4 мінімальне значення швидкості витрат на енергосистему Z= 831,278 грн. при w1 = 1,4 м/с , де d1 = 0,12437286 м , д1 = 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 г.
|