Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Современные проблемы теплогазоснабжения населенных мест и предприятий. Учебное пособие (курс лекций)

.pdf
Скачиваний:
0
Добавлен:
08.09.2026
Размер:
2 Мб
Скачать
☆
21
6 + + + – + + – + – – – 7 – + + + – + + – + – – 8 – – + + + – + + – + –
9 – – – + + + – + + – + 10 + – – – + + + – + + – 11 – + – – – + + + – + + 12 – – – – – – – – – – –
Сверхнасыщенные планы
Планы данного типа применяют в тех случаях, когда фактиче­ски значимых эффектов значительно меньше общего их теорети­ческого числа. Под термином «эффект» понимается влияние от­дельного (j-го) фактора и парного взаимодействия 2-х факторов.
Сверхнасыщенным является план, в котором число опытов меньше числа оцениваемых эффектов. Используя эти планы, нель­зя дать качественную оценку влияния всех факторов и их взаимо­действий на изучаемую величину. Однако основная цель поста­новки отсеивающих экспериментов заключается не в количествен­ной оценке влияния факторов и их взаимодействий на изучаемую величину, а в выделении значимых эффектов из их общего числа.
Сущность и основное достоинство методик предварительного экспериментального отсеивания факторов состоит в том, что в ре­зультате небольшого числа опытов удаётся отсеять группу (ряд) незначимых факторов, которые в противном случае (попав в об­щую матрицу планирования основной рабочей серии опытов) вы­звали бы лавину дополнительных опытов:
n
)1m(n22N
-++=
. (1.19)
В итоге мы бы их всё равно отсеяли, но гораздо более доро­гой ценой. Отрицательную роль несущественных факторов можно сравнить с вирусной инфекцией, которую гораздо вы­годнее погасить в зародыше. В качестве матрицы сверхнасы­щенного плана используют случайные выборки из систематиче­ских дробных реплик.
Число опытов в матрице планирования принимают таким, что­бы оно оказалось не меньше числа ожидаемых значимых эффек­тов. Обычно число опытов в матрице принимают кратным 8 и пре­вышающим (n+1), где n – число факторов.
22
Например, в задаче с десятью факторами требуется выделить значимые эффекты из числа 55 потенциально возможных (10 ли­нейных эффектов от каждого фактора и 45 эффектов парных взаи­модействий).
Предполагается провести 16 опытов. Для построения матрицы планирования факторы разбиваются на две группы по пять факто­ров в каждой. Для первой группы факторов в качестве матрицы планирования используем полуреплику
15
-
Матрицу планирования для остальных пяти факторов
)ххххх1(2
´´´´=
.
54321
х......х
106
получаем случайным выбором строк из той же полуреплики (мат­рицы
хх ¸ ).
51
Результаты проведенной выборки представляем в виде таблицы 1.8.
Таблица 1.8
Полуреплика 2
5–1
(
)ххххх1
´´´´=
54321
опыта
Факторы Факторы Номер
х1 х2 х3 х4 х
5
Номер
опыта
х1 х2 х3 х4 х
5
1 + + + + + 9 + + + – –
2 – + + + – 10 – + + – +
3 + – + + – 11 + – + – +
4 – – + + + 12 – – + – –
5 + + – + – 13 + + – – +
6 – + – + + 14 – + – – –
7 + – – + + 15 + – – – –
8 – – – + – 16 – – – – +
23
1.2.4. Принципы переноса результатов научных
лабораторных исследований на промышленные условия
В современном мире усложняющейся техники ее дальнейшее совершенствование невозможно без экспериментальных исследо­ваний, подтверждающих те или иные теоретические разработки и научные положения.
Физические модели – это реальные лабораторные стенды, кото­рыми в ходе эксперимента заменяют оригинал. Но после заверше­ния серии экспериментов необходимо по определенным методи­кам и правилам перенести полученные результаты на промышлен­ные условия. Эти задачи успешно решаются с помощью теории подобия.
Подобие – это условие, при котором возможен количественный перенос результатов опыта с модели на оригинал. Более того, по­добие модели и оригинала дает нам правило такого переноса.
По первой теории подобия одноименные критерии подобия в сходственных точках подобных объектов (модели и оригинал) равны, то есть количественный перенос сводится к тривиальному приравниванию критериев подобия.
Таким образом, при моделировании на основе метода подобия критерии подобия играют двоякую роль:
– во-первых, на основе этих критериев определяют, какой должна быть модель, чтобы обеспечить условия моделирования;
– во-вторых, значения этих критериев и есть та количественная мера, которая переносится с модели на оригинал.
Рассмотрим характерный пример: канал оригинал и модели.
Рис. 1.4. Сравнительная схема оригинала
и модели: 1, 2 – сходственные точки
24
Необходимо найти распределение скоростей жидкости в гори-
w
l
×
l
l
l
зонтальном канале сложной формы и скорости в разных его точ­ках, например, для предсказания возможности отложения осадков.
Прежде чем строить канал, необходимо знать это распределе­ние, и если оно окажется неудачным, то форму канала потребуется изменить. С этой целью распределение изучают на модели.
Каковы должны быть форма и размеры модели канала?
Форма модели определяется условиями геометрического подо­бия, а размеры – соображениями целесообразности: большая мо­дель дорога, а на маленькой резко возрастает роль мелких погреш­ностей и побочных факторов.
Какова должна быть скорость жидкости в модели?
В данном случае подобие определяется критерием Рейнольдса:
w
орор
где w – средняя скорость в канале.
Удобно использовать на модели ту же жидкость, что и в ориги­нале (особенно, если это вода). Тогда
модмодорор
n
ор
=
ww
n
×
мод
×=×
l
модмод
, (1.20)
n
=
n
ор
, следовательно:
мод
; (1.21)
w
мод
Получили, что скорость воды в модели должна быть обратно пропорциональна ее размеру: чем меньше модель, тем больше скорость. Если условие (1.22) соблюдено, то распределение скоро­стей в модели подобно их распределению в оригинале, и значения критерия Рейнольдса в сходственных точках равны.
В рассмотренном примере сравнение модели и объекта осуще­ствлялось по одному критерию. В случаях, когда сравнение прово­дят по двум и более критериям, построить модель оказывается крайне трудно.
Например, в ходе анализа критериев может выясниться обстоя­тельство, что для моделирования ускорения надо создать 1000·gор – величину, далеко превышающую технические возможности. Дру­гим вариантом нереальной задачи может служить вывод о необхо­димости развития скорости объекта до 1000 м/с [1].
w
ор
=
l
ор
. (1.22)
мод
25
Метод аналогии и его использование при моделировании
l
Аналогичными называются объекты и процессы, описываемые аналогичными уравнениями.
Уравнения аналогичны, если они содержат разные физические величины, но все операторы (знаки алгебраических действий, сим­волы функции, символы дифференцирования и интегрирования и т. д.) в них совпадают и следуют в одном и том же порядке.
В математике мы по существу не имеем дела с физическими ве­личинами, а оперируем лишь символами их обозначающими. По­этому аналогичность уравнений, описывающих различные физи­ческие объекты, приводит к одинаковому распределению числен­ных значений этих величин в пространстве и времени, хотя эти значения отображают разные физические величины [1].
Процесс Прохождение Тепловой
фильтрации электрического тока поток
r
-=
w
Закон Закон Закон
Дарси-Вейсбаха Ома Ньютона – Рихмана
Pqrad
R
, м/с
.фильтр
r
i -=
uqrad
, а
R
.эл
r
q -=
tqrad
R
, Вт/м∙К
.терм
Рассмотрим пример. Необходимо исследовать процесс неста- ционарного распространения тепла вдоль длинного стержня (тем­пература стержня меняется и по длине, и по времени).
Распределение температуры по длине в условиях нестационар­ного процесса описывается уравнением Фурье:
2
а
Т
Т
×=
, (1.23)
2
Т
l
=
а
где
t
Т
r
×
с
р
– время; l – длина; Т – температура.
t
, м2/с – коэффициент температуропроводности;
Уравнение, описывающее нестационарную электропроводность проводника, обладающего сопротивлением (R) и емкостью (С), имеет вид:
26
2
l
l
t
×
С
R
l
l
l
×
t
Аналогичность уравнений (1.4) и (1.5) очевидна, следовательно, тепловой оригинал можно моделировать электрической моделью.
Критерий подобия для оригинала – число Фурье:
t
1U
=
×
U
×
CR
. (1.24)
2
a
F
Соответствующий критерий для модели:
Условие подобия – это равенство критериев:
а
×
2
ор
Если оригинал не очень велик, то может оказаться удобным сделать модель приблизительно в натуральную величину.
Тогда
=
o
П
=
э
a
ор
2
мод
t
орор
=
. (1.25)
2
1
×
t
. (1.26)
2
××
t
=
мод
CR
××
t
мод
СR
2
мод
,
. (1.27)
2
ll ××
модмодмод
t
=×
t
t
орор
мод
ор
или
t
Для общеупотребительных значений а, R, С их произведение составляет величину порядка 10-5÷10-9.
Следовательно, предложенная электрическая модель чрезвы­чайно удобна, так как многократно сокращает время проведения исследования.
мод
СRа×
(1.28)
модмод
СRа ××=
модмодор
. (1.29)
27
Моделирование с помощью ЭВМ
Применение аналогии расширяет возможности моделирования, но не снимает трудностей, связанных с обеспечением подобия фи­зической модели и оригинала. Аналоговую модель оригинала большой сложности построить не удается.
Преодолеть эти препятствия и разработать методы моделирова­ния сложных тепло-массообменных процессов возможно, если при исследовании отказаться не только от физической идентичности модели и оригинала, но и от их подобия во времени и пространст­ве. Но тогда возникают вопросы: как строить модель? Будет ли это модель вообще? Что ее будет объединять с оригиналом?
Обратимся к сформулированным выше требованиям, предъяв­ляемым к процессу моделирования: это экономичность и возмож­ность переноса данных с модели на оригинал, что подтверждается или отвергается практикой.
Для практической технологии совершенно безразлично, с по­мощью каких технических средств реализуются эти требования. Более того, если можно решить задачу «в уме», не прибегая к мо­делированию, тем лучше.
В достаточно простых случаях расчетный метод, несомненно, предпочтительнее моделирования, потому что дешевле сосчитать, чем строить модель и экспериментировать. Например, для испыта­ния противоударного устройства наручных часов, необходимо уз­нать: какую скорость будут иметь часы, бросаемые с различной высоты? Для этого достаточно воспользоваться уравнением сво­бодного падения Ньютона:
В данном случае есть математическая модель, но нет процесса моделирования. Это не будет моделированием, так как моделиро­вание это экспериментальный метод, связанный с работой модели, а при решении уравнения (1.30) эксперимент отсутствовал.
По мере развития науки и совершенствования технологических процессов усложняются и математические модели. Часто склады­вается такая ситуация, когда математическая модель есть, но ре­шить это сложнее уравнение или систему уравнений человеку крайне сложно. Например, прочностные модели в авиастроении в 50–60-е годы описывались системами из 8÷10 дифференциальных уравнений, для решения которых требовалось более месяца.
w
Hg2 ××=
. (1.30)
28
В середине ХХ века число не решаемых задач возрастало лави­нообразно. Создавалось критическое положение, когда огромный потенциал, накопленный наукой и практикой, невозможно было использовать. Выход из этого кризиса, как известно, был найден с появлением быстродействующих вычислительных машин.
Рассмотрим процедуру решения технологической задачи:
а) расчетным путем;
б) на основе моделирования.
В первом случае главная часть работы проводится в уме (бума­га и калькулятор играют лишь вспомогательную роль): выбирается физический подход, метод решения, подбираются системы урав­нений и т. д.
Во втором случае решение получают, включая в работу некото­рое техническое устройство (стенд) – материальную модель, соз­данную на основе мысленной модели процесса. ЭВМ осуществля­ет перебор функций, систем уравнений в поисках оптимального описания. Результаты работы этого устройства позволяют полу­чить решение. С одной стороны, это не физический эксперимент, а расчет. Но в то же время для нас не важно, как работает та или иная установка. Важнее, кто или что работает – приходится ли всю задачу решать человеку, или значительная часть работы перело­жена на техническое устройство (модель или ЭВМ).
В этом смысле применение ЭВМ имеет существенное сходство с моделированием: здесь мы на основе мысленной модели процес­са – его математического описания – формулируем некоторую ус­тановку для ЭВМ, то есть, вводим в нее программу. Результатом работы этой программы является итоговое решение.
1.2.5. Проведение научных исследований
в условиях промышленного производства
При проведении научных исследований в условиях промыш­ленного производства наибольшее распространение получили ме­тоды адаптационной оптимизации.
Оптимизация процессов в производственных условиях имеет ряд особенностей, которые обусловлены следующими причинами.
1. В условиях действующего производства экстремальные зада­чи оптимизации приходится решать при непрекращающемся вы­пуске продукции.
29
2. В условиях действующего производства не представляется
возможным варьировать факторы в широких интервалах, так как на области изменения факторов налагаются ограничения, связан­ные с риском получения бракованной продукции.
3. В производственных условиях большую роль могут играть неконтролируемые переменные, изменение которых приводит к смещению (дрейфу) оптимума функции по отношению к контро­лируемым переменным. Так, например, в условиях массового про­изводства станки поточных линий работают в оптимальных режи­мах по контролируемым параметрам (скорость, подача, глубина резания и т. д.), обеспечивающих минимальный процентный вы­ход бракованной продукции.
Однако резкое повышение температуры наружного воздуха и относительной влажности летом или соответствующее их пониже­ние зимой могут
вызвать резкое увеличение выхода бракованной продукции.
4. В условиях действующего производства дисперсия воспро­изводимости эксперимента значительно меньше дисперсии вос­производимости лабораторного опыта.
В 1962 г. Спиндлей, Хетд и Химсуорз предложили метод адап­тационной оптимизации процессов. Этот метод планирования по­строен таким образом, что все решения в ходе экспериментирова­ния принимаются на основании определённых правил, то есть они как бы заранее регламентируются.
Формализация стратегии поиска оптимума позволяет использо­вать для управления технологическим процессом управляющие ЭВМ или персонал невысокой квалификации. Применение машин позволяет полностью автоматизировать процесс управления. Ли­ния с числовым программным управлением может работать как саморегулирующаяся система.
В поисках технологического оптимума очень часто используют шаговый метод, который предусматривает выполнение следующих этапов работ:
а) изучение некоторой области поверхности отклика с целью определения направления движения к оптимуму;
б) движение к оптимуму по установленному направлению.
Классическим примером реализации такой методологии поиска является метод крутого восхождения Бокса – Уилсона.
30
Рассмотрим функцию отклика
=
)x,x(fy21=
, которая пред-
ставляет собой уравнение поверхности. Эту поверхность можно представить линиями равного отклика на факторной плоскости
xx
: линии
21
consty1i=
-
consty
;
i
;
consty1i=
+
. Вблизи мак-
симального значения величин (у) такая поверхность имеет вид холма, вершина которого представляет собой точку, соответст­вующую оптимальным условиям протекания процесса.
Для поиска оптимума вначале реализуют полный факторный эксперимент 22 (опыты 1, 2, 3, 4), что отвечает шагу «а» – опреде­ление направления движения к оптимуму. По итогам выполнен­ных опытов определяют коэффициенты линейного уравнения рег­рессии и направление градиента:
ygrad
+
=
x
. (1.31)
x
21
y
y
После этого производят крутое восхождение (опыты 5, 6, 7, 8). Приняв точку с наибольшим значением у (опыт 7) за центр плана, вновь реализуют полный факторный эксперимент 22 (опыты 9, 10, 11, 12) и определяют новое направление градиента.
Последним этапом поиска оптимума является движение по но­вому градиенту 22 (опыты 13, 14, 15, 16), где и находят максимум функции (опыт 15). Схема движения к оптимуму функции по ме­тоду «крутого восхождения» представлена на рисунке 1.5.
Симплексный метод решения задачи базируется на построе­нии геометрических фигур – симплексов.
Симплекс – простейшая геометрическая фигура, образованная множеством n+1-независимых точек в n-мерном факторном про­странстве. Точки, образующие симплекс, называют его вершина­ми. Так, в двухфакторном эксперименте симплекс – это треуголь­ник, а в трехфакторном эксперименте – треугольная пирамида.
Основной особенностью симплексного планирования является совмещение процесса изучения поверхности отклика с процессом движения по самой поверхности отклика.
Указанное совмещение достигается в результате того, что усло­вие проведения каждого очередного опыта устанавливается на ос­нове оценки результатов предшествующих опытов, поставленных в точках факторного пространства, соответствующих вершинам симплекса.
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]