КОМПЬЮТЕРНОЕ ПРОГНОЗИРОВАНИЕ ТЕПЛОВЛАЖНОСТНЫХ И ВОЗДУШНЫХ РЕЖИМОВ В ПРОИЗВОДСТВЕННЫХ ПОМЕЩЕНИЯХ
Так как почти все факторы, формирующие микроклимат, непрерывно изменяются во времени, тепловое равновесие даже в помещениях с автоматически регулируемыми отопительно-вентиляционными приборами никогда не достигается. Поэтому стабильность внутренней температуры может рассматриваться только как результат взаимно уравновешивающегося воздействия очень большого числа динамических элементарных процессов.
Регрессионная модель микроклимата.
При планировании эксперимента по исследованию микроклимата в качестве параметров оптимизации были назначены три главные хаpaктеристики тепло-влажностного и воздушного режимов помещения: температура T, относительная влажность j и подвижность воздуха V.
Функциональная связь между рассматриваемым климатическим параметром y и прострaнcтвенно-временными координатами определялась в виде полинома второй степени:
. (1)
Регрессионные коэффициенты
при при при вычисляются по методу наименьших квадратов из системы нормальных уравнений
(2)
где матрица неизвестных матрица выходных параметров испытаний Общее число опытов N зависит от типа факторного плана, по которому проводится эксперимент. Информационная матрица Фишера строится с помощью матрицы условий
Условия проведения опытов организуются по одному из симметричных центральных планов второго порядка - рототабельному композиционному плану РКП или D-оптимальному четырехфакторному плану типа Бокса B4 |2, 3]. РКП строится униформным с величиной безразмерного момента матрицы планирования l3 = 0,8705. Управляемые факторы кодируются на пяти уровнях, звездное плечо a == 2, общее число опытов N = 31. ядро плана состоит из 16 точек, количество дублирующих опытов в центре n0 = 7.
При исследованиях по плану B4 входные параметры кодируются на трех уровнях, N = 24, число звездных точек с плечом a = 1 составляет пa = 8.
Элементы стандартной матрицы С = находят в каталоге планов эксперимента для факторных и полиномиальных моделей.
После проведения натурных экспериментов на основе полученного статистического материала во всех помещениях здания по формулам (2) на ЭВМ рассчитывают значения регрессионных коэффициентов функции отклика (1) для каждого из параметров оптимизации t, j и v.
Проверка статистической значимости полученных коэффициентов осуществляется с помощью доверительных интервалов. Исключаются из модели (1) без пересчета статистически незначимые оценки при исключении одного из значений требуется пересчитывать остальные коэффициенты из этой группы по формулам (2).
Гипотезы об информационной способности и адекватности построенных полиномиальных моделей проверяются методами регрессионного анализа.
Пример использования композиционных симметрических планов Бокса В4 и РКП при оптимизации микроклимата.
Приведем результаты исследований по выбранным факторным планам тепло-влажностного и воздушного режимов в помещениях, расположенных в двухэтажном кирпичном здании [3].
Эксперимент проводился в переходный период отопительного сезона (март-апрель). В здании работала двухтрубная тупиковая система отопления с нижней разводкой, для прогрева высоких помещений были установлены напольные конвекторы, в специальном помещении первого этажа функционировал элеваторный узел.
Значительный налет ржавчины и солей на внутренних стенках отопительных приборов и трубопроводов, возникший из-за сложной конфигурации системы, отсутствия гибкой регулировки между отдельными узлами, невыполнения в течение длительного периода эксплуатации должной очистки и профилактики; отсутствие пылеуловителей, кондиционеров и других воздухоочистительных приборов необходимой мощности, позволяющих отводить избытки тепла, влаги, углекислого газа и других нежелательных веществ, а также многие другие вредные факторы вызвали нeблагоприятные изменения метеорологических условий воздуха в обслуживаемой зоне помещений (табл. 1).
Таблица 1. Максимальные абсолютные величины отклоненийэкспериментальных значений от оптимальных норм
Параметры |
Номер помещения |
Примечание |
|||||
1 |
9 |
3 |
4 |
5 |
6 |
||
t, °C |
7,3 |
3,5 |
7,2 |
7,7 |
7,7 |
7,6 |
Перетоп |
j, % |
42,4 |
45,4 |
46,6 |
44,4 |
44.7 |
45,2 |
Сухость |
v×10-1, м/с |
0,75 |
0,76 |
0,77 |
0,70 |
0,55 |
0,77 |
Недостаточная вентиляция |
Инструментальные обследования метеорологических условий воздуха в помещениях здания показали существенные отклонения экспериментальных значений параметров микроклимата от оптимальных норм, рекомендуемых СНИП.
Для разработки мероприятий по нормализации микроклимата в каждом из помещений при одинаковых (по возможности) условиях были проведены полные факторные эксперименты по планам Бокса B4 и РКП.
Процесс организации вычислительных процедур продемонстрируем на примере моделирования температурного поля в помещении размерами 3,00 x 9,04 x 8,12 м, расположенном на 2-м этаже здания.
Таблица 2. Кодирование факторов в РКП
Кодовое обозначение |
Уровни варьирования |
||||
-2 |
-1 |
0 |
1 |
2 |
|
х1 |
10 |
12 |
14 |
16 |
18 |
х2 |
0,50 |
1,00 |
1,50 |
2,00 |
2,50 |
х3 |
0,50 |
2,51 |
4,52 |
6,53 |
8,54 |
х4 |
0,50 |
2,28 |
4,06 |
5.84 |
7,62 |
Для того чтобы воспользоваться стандартными матрицами планирования X натурные значения факторов времени х1 и прострaнcтвенных координат х2, х3, х4 переводятся в безразмерные. В табл. 2 показаны 5 уровней варьирования независимых переменных в рототабельном композиционном плане. По плану Бокса В4 управляемые факторы кодируются на трех уровнях с интервалами варьирования Dх1 = 4, Dх2 = 1. Dх3 = 4.02 и Dх4 = 3,56.
Опытные данные, полученные при исследовании температурного поля в выбранном помещении, приведены в табл. 3. Регрессионные коэффициенты уравнения температур (1) рассчитываются по формулам (2), табл. 3. Проверка информационной способности и адекватности полученных полиномиальных моделей показала: по РКП Fинф = 2,39, Fад =1,19, Fтаб = 3,91; по плану Бокса Fинф = 1,93, Fад = 2,10.
Так как вычисленные критерии Фишера меньше табличных, можно предположить. что обе модели адекватны реальному процессу изменения температуры в помещении.
Анализ результатов.
В результате анализа рассчитанных полей температур, влажности и подвижности воздуха был сделан вывод, что наблюдается перегрев помещений и здания в целом, который приводит к всплытию воздушных потоков в верхнюю часть помещений и возникновению тепловых подушек с одновременным понижением и без того низкой влажности, с недостаточностью воздухообмена и образованием застойных зон. На основе полученных данных был разработан комплекс мероприятий, который позволил привести метеорологические условия и чистоту воздуха в обслуживаемой зоне помещений в соответствие с нормами.
Применение статистических методов дает возможность в оптимальные сроки с минимальным количеством опытов определить и исследовать математические зависимости основных параметров тепло-влажностного и воздушного режимов от времени суток и размеров помещения при сформировавшихся в установившемся эксплуатационном режиме здания метеорологических условиях. Планирование экспериментов по исследованию микроклимата в помещениях одновременно по двум симметричным композиционным факторным планам 2-го порядка показало, что количество опытов, требуемых для получения одинаковой максимальной погрешности результатов эксперимента при расчетах на моделях РКП в 2-3 раза меньше, чем при использовании планов типа В4 . Однако в случаях, когда хотя бы одна из стен помещения является наружной, РКП может приводить к неадекватным математическим моделям. Такой эффект можно объяснить асимметричностью расположения в помещении оконных и дверных проемов, через которые происходит неорганизованное поступление наружных воздушных масс, тогда как униформность РКП предполагает пpaктически постоянную дисперсию рассеяния прогнозируемых значений функции отклика вокруг центра эксперимента.
Таблица 3. Коэффициенты уравнения регрессии для температуры
По РКП |
По плану В4 |
||||
Рассчитанные |
Интервал значимости |
Исправленные |
Рассчитанные |
Интервал значимости |
Исправленные |
b0 = 25,4910 |
Db0 = 0,1578 |
b0 = 25,3615 |
b0 = 24,5083 |
Db0 = 0,5238 |
b0 = 24,5083 |
B1 = 0,2584 |
Db1 = 0,08523 |
b1 = 0,2584 |
b1 = - 0,09444 |
Db1 = 0,2829 |
b1= 0 |
b2 = 0,9417 |
b2 = 0,9417 |
b2 = 0,2722 |
b2 = 0 |
||
b3 = - 0,01667 |
b3 = 0 |
b3 = - 0,1167 |
b3 = 0 |
||
b4 = 0,08333 |
b4 = 0 |
b4 = 0,005556 |
b4 = 0 |
||
b11 = -0,3568 |
Db2 = 0,07809 |
b11 = - 0,3481 |
b11 = - 0,8583 |
Db2 = 0,2592 |
b11 = - 0,8583 |
b22 = - 0,9318 |
b22 = - 0,9231 |
b22 = - 0,7583 |
b22 = - 0,7583 |
||
b33 = - 0,05683 |
b33 = 0 |
b33 = 1,0417 |
b33 = 1,0417 |
||
b44 = - 0,05683 |
b44 = 0 |
b44 = 1,1417 |
b44 = 1,1417 |
||
b12 = 0,1000 |
Db3 = - 0,1044 |
b12 = 0 |
b12 = - 0,3375 |
Db3 = 0,3465 |
b12 = 0 |
b13 = 0,01250 |
b13 = 0 |
b13 = - 0,02500 |
b13 = 0 |
||
b14 = - 0,0875 |
b14 = 0 |
b14 = - 0,03750 |
b14 = 0 |
||
b23 = - 0,1875 |
b23 = - 0,1875 |
b23 = 0,08750 |
b23 = 0 |
||
b24 = 0,1625 |
b24 = 0,1625 |
b24 = 0,02500 |
b24 = 0 |
||
b34 = 0,2250 |
b34 = 0,2250 |
b34 = 0,03750 |
b34 = 0 |
СПИСОК ЛИТЕРАТУРЫ
- Адлер Ю.П., Маркова Е.В., Грановский Ю.В. Планирование эксперимента при поиске оптимальных условий - М.: Наука, 1976.
- Сергиенко Л.С., Житов В.Г. Исследование метеорологических условий в помещениях жилых и общественных зданий с применением математических методов планирования эксперимента //Известия высших учебных заведений Министерства образования РФ / Ежемес. науч.-теорет. журн. "Строительство" - № 6 (534). - Новосибирск: Изд-во Новосиб. гос. архитектурно-строительного ун-та, 2003. - С. 63-67.
- Сергиенко Л.С., Житов В.Г. О компьютерном моделировании микроклимата в здании //Матем. модели и методы их исследования / Тр. междунар. конф., Т. 2. - Красноярск: Изд-во ИВМ СО РАН, 2001. - С. 191-195.
Статья в формате PDF 113 KB...
19 04 2024 7:22:49
Статья в формате PDF 258 KB...
18 04 2024 0:35:44
Статья в формате PDF 311 KB...
16 04 2024 5:32:59
Статья в формате PDF 236 KB...
15 04 2024 5:32:18
Статья в формате PDF 119 KB...
14 04 2024 7:27:15
Статья в формате PDF 104 KB...
13 04 2024 1:57:17
Статья в формате PDF 246 KB...
11 04 2024 7:19:29
Статья в формате PDF 102 KB...
10 04 2024 15:50:20
Статья в формате PDF 100 KB...
09 04 2024 23:19:59
Статья в формате PDF 309 KB...
08 04 2024 21:12:42
Статья в формате PDF 218 KB...
07 04 2024 5:15:26
Существующие методы атомной эмиссионной спектроскопии для исследования состава металлов и сплавов используются во всех отраслях машиностроения. По мнению авторов, современные методы уже не обеспечивают необходимых точностей измерений. В данной работе авторами проведены исследования влияния внешних факторов на точность измерений прибора атомно-эмиссионной спектроскопии. ...
06 04 2024 17:40:31
Статья в формате PDF 161 KB...
05 04 2024 12:57:31
Статья в формате PDF 267 KB...
04 04 2024 6:57:21
Статья в формате PDF 284 KB...
02 04 2024 21:38:19
Статья в формате PDF 139 KB...
01 04 2024 19:24:11
Статья в формате PDF 126 KB...
31 03 2024 14:15:24
Статья в формате PDF 161 KB...
30 03 2024 2:33:41
Статья в формате PDF 129 KB...
28 03 2024 6:13:25
Статья в формате PDF 116 KB...
26 03 2024 16:52:25
Статья в формате PDF 125 KB...
25 03 2024 10:27:30
Применение хитинсодержащих препаратов оказывает положительное влияние на мясную продуктивность бычков, а превосходство по хаpaктеристикам химического состава и энергетической ценности мякоти имеют бычки, получавшие сукцинат хитозана. ...
24 03 2024 8:51:53
Статья в формате PDF 108 KB...
23 03 2024 14:56:19
Статья в формате PDF 115 KB...
22 03 2024 17:41:39
Статья в формате PDF 262 KB...
21 03 2024 0:14:27
Статья в формате PDF 134 KB...
20 03 2024 5:59:34
Статья в формате PDF 396 KB...
19 03 2024 23:54:20
Статья в формате PDF 102 KB...
18 03 2024 3:40:38
Статья в формате PDF 104 KB...
17 03 2024 7:16:22
Статья в формате PDF 117 KB...
16 03 2024 8:21:51
Статья в формате PDF 100 KB...
15 03 2024 7:28:41
Статья в формате PDF 254 KB...
14 03 2024 0:55:35
Статья в формате PDF 127 KB...
13 03 2024 19:43:16
Статья в формате PDF 149 KB...
12 03 2024 4:53:32
Статья в формате PDF 114 KB...
11 03 2024 18:31:20
Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::