Методы максимальной формализации машинных расчетов
Рассматриваемые методы основаны на системном анализе принципов расчета теплообменных аппаратов, направлены на создание иерархической системы различных видов расчетов. При разработке таких методов предлагаются четкие функциональные классификации, обобщенные структуры расчетов, построение специфических модулей на основе ограниченного их числа [28J.
Для разработки обобщенных методов расчета теплообменников предлагается следующая структура: элементарным звеном теплообменника является элемент — характерная часть аппарата (или весь аппарат) с элементарной схемой тока теплоносителей, аппарат — следующее звено иерархической структуры, состоящее из теплопередающнх элементов и полостей для дннженпя теплоносителей (число, состав, схема соединения элементов в аппарате могут быть любые); совокупность аппаратов, соединенных в любой последовательности, образует следующую иерархическую ступень — теплообменник (его основное назначение—изменение температуры, агрегатного состояния и физико-химиче-
С пуск ) |
/-Я/-1— /Вбод и 5ы5од Исходных Данных |
Вычисление Не *и1 -1й1 |
Т |
® @) — Л2 -------------------------- г - /и-1— |
■Л2 Вычисление Сг£ - *21, /л, Г - В2 Во/числение ап ' °‘21 ' Г-^ Вычисление & £1 |
Вычисление Й) |
Г-112 |
Т |
I |
Вычисление ^г! Г *-221 ~^г |
Г—2/--1- - Вычисление |
Г-£Г± |
Вычисление |
|
||
|
||
|
||
5.3. Блок-схема поимтервально-нтерационного расчета теплообменника
Ского состава теплоносителей); высшим звеном иерархии является система теплообменников, представляющая собой совокупность аппаратов, расположенных в любой последовательности.
Иерархическая структура такого типа позволяет на каждом уровне разрабатывать комплекс программ, формализующих методику расчета звена данного уровня.
В основу проектирования положен элемент («теплопередающий элемент»). Следовательно, при расчете коэффициента
|
|
|
Теплоотдачи на поверхностях элемента и коэффициента теплопередачи для выбранного типа элемента необходима четкая классификация его по форме основной поверхности, наружной поверхности наружного и внутреннего оребренйя. Поскольку на особенности механизма теплопереноса оказывает определенное влияние и взаимное расположение элементов (например, одиночная труба, шахматный пучок и т. д.), то способ компоновки, поверхности должен быть классифицирован.
Следующим звеном внутренней иерархии элемента должна быть классификация схем тока теплоносителей в элементе, так как она определяет универсализацию и формализацию расчета среднего температурного напора в рамках элемента («теплопередающего элемента»), В качестве примера классификации элементарных схем тока [28] предложен рис. 5.4.
Разработка методов максимальной формализации машинных расчетов требует создания иерархии вычислительных структур. Для этого необходимо проанализировать состав всех основных видов расчета теплообменников, чтобы выделить минимум элементарных структур (модулей), из которых можно построить алгоритм любого расчета.
Рассмотрим структуру проектных расчетоп теплообменников (рис. 5.5):
А1. Ввод исходных данных. В зависимости от состава проектного расчета состав исходных данных может быть различным. В практике построения алгоритмов расчет свойств теплоносителей рассматривают как самостоятельную задачу (внешнюю задачу).
51. Решение задачи теплового баланса (БС — ТБ) сводится к определению одной из неизвестных величин из набора? ц, ^12,
(22, Си С2. Пример алгоритма решения уравнения теплового баланса приведен в работе [37].
1, П. Расчет скорости теплоносителя, отдающего (БС — х>]) и воспринимающего тепло (БС — ш2), необходим для определения гидравлических сопротивлений.
£1, Н1. Проверка ограничений скоростей (БС — Ш1ДОп, БС — ^2доп) ДЛЯ обеспечения УСЛОВИЙ
Ш[ € (ЬУ|ш1п> ^1тах)» ^2 С (К^гтпп» ^2тах)-
Если не удовлетворяется одно из ограничений, вырабатывается признак П>0.
21. Расчет коэффициента теплопередачи (БС — к). Точность расчета к обеспечивается итерационным определением коэффициентов теплоотдачи на поверхностях.
К1. Расчет требуемой поверхности теплообмена (БС — И).
А2, С2. Расчет гидравлических сопротивлений теплоносителей— отдающего (БС — ДР|) и воспринимающего (БС — ДР2)-
В2, 02. Проверка ограничений гидравлических сопротивлений служит ДЛЯ обеспечения условий ДЯ1 < ДР|доп, Ь. Р-2 < ДРгдоп - Если не удовлетворяется любое условие, вырабатывается признак П > 0.
Е2, Расчет мощности нагнетателей теплоносителя — отдающего (БС — N1) и воспринимающего (БС — Л'г).
12. Печать результатов — перечень величин, выводимых на печать.
К2. Проверка наличия выхода за ограничения. Если результаты расчетов не выходят за ограничения, все признаки П равны пулю. Соответственно 1П = 0 и проектный расчет заканчивается.
ЛЗ. Изменение типоразмера аппарата (БС — ИТА). Обращение к этой структуре следует, если хотя бы один из признаков выхода за ограничения не равен нулю.
Аналогична блочная структура и других типов расчетов. Среди типовых структур можно выделить такие: условно-постоянные, основанные на использовании нормативных документов, методик, в которые вносятся изменения и коррективы; оспопан - ные на использовании фундаментальных закономерностей, методов, надежных логических связей, рассчитанных па длительное использование. Поэтому в обобщенных методах максимальной
Рис. 5.5. Блок-схема просктпого расчета рекуперативной теплообменника |
Формализации существенна разработка постоянных структур:
БС — к. Расчет коэффициента теплопередачи.
БС — Р. Расчет площади теплопередающей поверхности теплообменника.
БС — Фээ. Расчет функции эффективности элемента.
БС — Расчет распределения температур теплоносителей в теплообменнике.
БС — ПЭЦФ. Поиск экстремума целевой функции.
Основой метода максимальной формализации служит универсальная модель переноса тепла в элементе:
Ф’=тЬтГ=^ (5.2) |
*11 ~*12 Р *11 — (21
Где Р — безразмерный температурный комплекс (1.46);
2 ехр (52) ^
Здесь КТГП)Г=Г4Д(б.4);5 =-1^ (б-5); Р€[0,11-
Индекс противоточности элемента. При 52 > 7 с погрешностью + + 0,1 %
Фэ = г + (А+ о-
Из уравнений (5.1) — (5.6) следует выражение для площади теплопередающей поверхности;
<У, . 2А --Р [г - (А + I)] кг 2А — Р1Х + (А + 1)] •
Блок-схема (рис. 5.6) расчета поверхности теплообмена элемента, вытекающая из уравнения (5.7), достаточно проста.
Для вычисления параметра 1 необходимо знать индекс противоточности. Индекс противоточности р — единственная характеристика, однозначно определяющая схему тока в элементе. Вместе с тем она является косвенным показателем теплопередающего совершенства схемы тока. При противотоке р= 1, при прямотоке р = 0. Во всех остальных случаях рЄ [0,1]. Чем ближе р к единице, тем эффективнее теплопередача в элементе (табл. 5.1).
Расчет функции тепловой эффективности Фэ можно осуществить по схеме рис. 5.7. Эта структура достаточно проста, может использоваться в машинном и ручном счете.
И & |
|
-Ьг •£*г |
И-об- |
Несимметричный разный элемент |
: (кР) |
(№) |
Прям |
Про Г |
(кР) Т |
Прот Прям |
|
|||
|
|||
|
|
||
|
|
||
|
|||
|
|||
|
|
||
|
|||
|
|||
|
|||
|
|||
|
|||
|
|||
|
|||
|
|||
|
|||
|
|||
|
|||
|
|||
|
|||
|
|||
— |
|||
№ Схемы |
Схема |
Вид схемы ток? |
Р |
9 |
Параллельно смешанный ток с двумя противоточ - нымн ходами |
І Г-+ |
0,70 |
10 |
Перекрестный ток с перемешиванием обеих сред однократно |
1" , -~1гг Ип |
0,56 |
||
11 |
Перекрестный ток с перемешиванием обеих сред, двухкратный общий противоток |
22 |
^гг |
0,88 |
|
12 |
Перекрестный ток с разделом одной среды на струи, однократный |
Г12 |
0,.і7 |
||
13 |
Перекрестный ток с разделами обеих сред на отдельные струн, однократный |
! * Vг |
0,76 |
||
И Перекрестный ток с. разделами одной среды на струи, двухкратный общий протшюток
© Г |
||
Вычисление Р |
Вычисленил? |
|
Формула а45) |
Ірормула (5.7) |
|
Д* 1 |
I |
|
2а |
||
1 +(А+1} |
© Г-сз^Х---------- Изменение услади а расчета |
|
Ртах'^>—(Ну |
Г” Л2- БС - К |
С п*ск ) |
■ А1 |
/ІВВод Исходных Данных ZZJ~ |
© |
Вычисление А формула. |
І 5.1/ |
(—01 -1------ Вычисление I /рормула (5, 4) |
Печать
Ос тан об
SHAPE \* MERGEFORMAT
Рис. 5.6. Блок-схема расчета поверхности теплообмена
Хая расчета температурной поправки гд, при схемах течения, отличных от противоточной [28], предложено следующее уравнение:
1 - А’Р
1п-
Ел/ = у----------------- ; (5.8)
Л'Р
При Н = 1 ал, = у. --- • (5.9)
1 ІУІтах
2
Здесь Ртах = (д+-1) + г (5.10); г =]/(/?+ 1)2 _4р/? (5.11); ^ і г /і ур (5*12), где N = Р/Ртах (5.13).
ІП I 'г тах)
1 — N
( Пуск} ( 5£сЗ Заьпс!*. |
__ В1 - I--- Вычисление А формула (5-1) 1Г77~Х= |
Вычисление I формула (5 4) |
^0! Вычисление 5 формула (5 5) Г |
Рис. 5.7. Блок-схема расчета тепловой эффективности Ф, |
Зачисление Ф} формула (5-6) |
Вычисление &} (рсрмула (53) |
^’ Ос та »об ^ |
Блок-схема расчета температурной поправки достаточно проста (рис. 5.8).
Более детальные сведения о методе определения поверхности ряда элементов даны в работе [28].