УДК 519.233.5 + 519.246.85
ФИЗИЧЕСКИЙ ЭКСПЕРИМЕНТ В МАТЕМАТИЧЕСКОМ МОДЕЛИРОВАНИИ ДИНАМИЧЕСКИХ СИСТЕМ С ДЕТЕРМИНИРОВАННЫМИ ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ
С.С. Торбунов
Алтайский государственный технический университет им. И.И. Ползунова, г. Барнаул E-mail: sle@agtu.secna.ru
Построение регрессионных математических моделей нелинейных и стохастических процессов во многих исследованиях является единственным методом, позволяющим получить адекватное решение поставленной задачи. Физический эксперимент в сочетании с аналитическими методами и принципом внешнего дополнения, заключающемся в эмпирическом упрощении условий задачи с целью получить конечное решение, приводит к эффективным результатам в различных отраслях науки. В настоящей статье приведены три примера его использования для исследования устойчивости двухфазной высокотемпературной струи в ограниченном пространстве, в порошковой металлургии и исследования движения твердого тела в сопротивляющейся среде.
Важность эксперимента в современной науке и технике не вызывает сомнения. Чтобы представить себе масштабы повседневной экспериментальной работы, достаточно к натурным исследованиям, проводимым в области фундаментальных наук и при проектировании новой техники, добавить испытания по образцам опытной и серийной продукции на тысячах заводов страны.
«Под экспериментом подразумевается общий процесс научного исследования и получения новых данных» [1. С. 3].
Технический прогресс приводит не только к увеличению сложности объектов испытаний - аппаратуры и оборудования, выпускаемых промышленностью; одновременно повышаются требования к глубине проникновения в сущность функционирования составных частей и взаимодействия между ними, а также к точности измерения регистрируемых величин и т. п.
Математическая модель становится неотъемлемым элементом испытаний, без построения которой невозможно осуществить планирование эксперимента, его проведение и обработку результатов. Модели, линейные по параметрам, стали в настоящее время традиционными. Гораздо реже используют модели, нелинейные по параметрам. Такие модели применяются при исследовании динамических объектов, кинетики и механизма явлений. При изучении динамических процессов зачастую на первое место ставится задача качественного исследования, а именно: выявление полного вектора определяющих параметров, установление интервалов изменения их значений для устойчивого режима или поддержания определенного режима протекания физического процесса, оптимизация процесса. Большой интерес в инженерных задачах представляют физические процессы, рассматриваемые в математической физике и механике сплошной среды, которые, как правило, являются стохастическими, что также необходимо учитывать при проведении эксперимента.
Условие детерминированности параметров позволяет упростить построение математических моделей стохастических систем в виде дифферен-
циальных уравнений динамических систем с переменными параметрами, а введение доказательных внешних дополнений, основу которых составляют уравнения и функции регрессий, не нарушая адекватности модельной системы, приводит к конечному решению; при этом аналитическая структура математической модели позволяет для качественного исследования динамической системы использовать теорию чувствительности и теорию устойчивости динамических систем по части переменных.
Конечно, даже при наличии сформулированных общих черт эксперимент не стал бы предметом отдельной науки, если бы задачи, решаемые различными исследователями, не имели ничего общего и полностью определялись спецификой той области знаний, где они работают. Однако это не так. Оказывается, можно выделить типовые задачи исследования, с которыми приходится сталкиваться каждому экспериментатору. К основным, наиболее распространенным типовым задачам исследования обычно относят:
1) получение некоторых предварительных сведений о процессе (обработка литературных данных, опросы специалистов и анализ результатов опросов, отсеивающий эксперимент);
2) получение формальных зависимостей;
3) проверку гипотез, т. е. некоторых содержательных предположений о свойствах объекта;
4) оптимизацию свойств изучаемого объекта (определение оптимальных соотношений, слежение за оптимумом и т. д.).
В построении математических моделей динамических систем физический эксперимент необходим в следующем:
1) построение статистической модели,
2) выявление кинематических зависимостей между параметрами физического процесса,
3) доказательство адекватности математической модели.
Основные принципы планирования эксперимента: отказ от полного перебора возможных входных состояний, принцип постепенного усложне-
ния математической модели (принцип последовательного планирования), принцип рандомизации (принцип приведения к случайности), принцип оптимальности планирования эксперимента.
Наиболее результативным является регрессионный эксперимент. Основное назначение регрессионного анализа - получение по экспериментальным данным регрессионных моделей объектов исследования, т. е. моделей вида у=Дх,г)+е, где х -вектор контролируемых и управляемых факторов (параметров), г - вектор только контролирующих факторов, е - аддитивная помеха. Построение регрессионной модели позволяет исследовать параметрическую устойчивость динамической системы, что решает задачу исследования устойчивости в условиях неприменимости теории устойчивости Ляпунова.
Действительно, после введения в исследование параметров Я=Я1,...,Я„ система дифференциальных уравнений первого порядка с переменными у=у1,...,у имеет вид:
у, _ Г (у, Л(0\ I), I _ 1,2,..., п,
решение которой у;=у;(А(0),/).
Возмущение параметров ДХк=рк, к=1,2,...,п, приводит к возмущенным решениям, после разложения правой части которых в ряд Тейлора и, ограничившись членами первого порядка вариаций /Зк, получим:
а ГУ ^ дГ {У ду± Ау дГ в
Введем чувствительность /-ой переменной по к-му параметру:
ду/дХк=зл
- и, учитывая независимость вариаций /Зк, приходим к уравнениям чувствительности
которые можно назвать вариационными уравнениями.
Устойчивость к возмущению параметров существует, если чувствительность ^¡к\<е, где е - малая заданная величина.
В общем случае следует из системы дифференциальных уравнений первого порядка. Однако, можно исключить влияние времени на чувствительность, полагая, согласно теории орбитальной устойчивости Пуанкаре, ¿¡¡=0. Оправданием этого служит также то, что момент времени, при котором условие я&=ю вызывает неустойчивость, абсолютно произволен, поскольку неизвестно, как учитывать время при исследовании чувствительности параметров. Поэтому для регрессионной модели
у=¥(Х1,...,Хк) получаем следующее выражение для чувствительности:
%=дк"/Д,
где Д - детерминант системы (*), а Д^ - детерминант, полученный из д заменой соответствующих столбцов векторами (-д¥/дХк).
Понятно, что s¡k=ж (т. е. неустойчивость) следует из условия Д=0. Это условие определяет Хкр, т. е. критическое значение параметра, вызывающее неустойчивость.
Одной из центральных теорем для всей теории параметрической устойчивости динамических систем является следующая теорема: «Устойчивость по Ляпунову и параметрическая устойчивость динамических систем в фазовом пространстве по чувствительности физических координат к изменению параметров эквивалентны» [2].
Достижением этой теории устойчивости динамических систем является ее применимость к исследованию устойчивости неконсервативных неавтономных динамических систем в фазовом пространстве, причем понятие чувствительности к возмущениям параметров может быть использовано в качестве способа оценки устойчивости, а устойчивость по чувствительностям всех физических координат к изменению параметров их полной совокупности является устойчивостью динамической системы с переменными параметрами в целом.
Таким образом, теория параметрической устойчивости является необходимым инструментом регрессионного анализа при качественном исследовании динамических систем с переменными параметрами в физическом эксперименте.
Построение регрессионной модели позволяет также решить задачу параметрической оптимизации динамических систем с детерминированными переменными параметрами [2].
В качестве примеров использования физического эксперимента в построении регрессионных моделей динамических систем с детерминированными переменными параметрами приведем следующие задачи.
Устойчивость двухфазной плазменной струи в ограниченном пространстве. В [2] исследовано возникновение неустойчивости и самовозбуждения колебаний, которое представляют широчайший интерес для многих разделов науки и техники. Для того, чтобы правильно моделировать техническую систему и адекватно записать её дифференциальные уравнения, необходимо иметь ясное качественное представление о механизмах и причинах возможной неустойчивости и возбуждения колебаний, позволяющее найти нужное для устранения самовозбуждения изменение конструкции или технологического режима. Такая задача возникает при исследовании движения горящей газовой (плазменной) струи в ограниченном пространстве (например, трубе) в процессе напыления легирующих элементов на подложку (например, с использованием самораспространяющегося высокотемпературного
синтеза - СВС-процесса). Причиной возникновения колебаний могут быть отрицательные трения или направленные воздействия, то есть имеют место автоколебания. Для определения условий устойчивости газовой струи необходимо знать параметры процесса, которые для данной технической системы могут быть определены только проведением специального эксперимента, а поэтому предпочтительнее такое исследование проводить по результатам построения регрессионной модели состояния газовой струи для определённого сечения в зависимости от параметров динамической системы.
При исследовании устойчивости технологического процесса, характеризуемого переносом массы и энергии (а именно это происходит при получении новых материалов наплавкой легирующих металлов в СВС-процессе) тепло- и массообмен связаны друг с другом. Состояние потока массы в некотором сечении газовой (плазменной) струи с частицами легирующего металла определяется критерием Рейнольдса
^ vd
Яе = — = Я1;
V
где V - скорость струи, V - кинематическая вязкость, й - средний диаметр частиц металла, и числом Прандля, характеризующим физические
свойства окружающей среды,
*
Рг = = к2,
а
где а* - коэффициент температуропроводности, V * - кинематическая вязкость окружающей среды. Иногда удобнее ввести критерий Нуссельта
хт а^
Ми = ,
к
где ак - коэффициент конвективного теплообмена (теплоотдачи), к - коэффициент теплопроводности окружающей среды.
Фазовой координатой движения массы является критерий Эйлера
Р =
ру2 У1'
а теплового состояния - средняя температура струи в данном сечении =у2.
По результатам экспериментальных замеров пусть построены в данном сечении струи регрессионные уравнения
у^Оъу^кА),
Уг=Рг(Ух,У2,к1,к1).
При исследовании параметрической устойчивости плоской динамической системы для двух параметров согласно (4) при Д=0 получим:
Щ = 0 Щ_дЕк = 0
5у ду1 ' ду2 ду2 '
Эти уравнения позволят определить критические значения параметров динамической системы,
при которых плазменная струя становится неустойчивой.
Следует также заметить, что фазовыми координатами могут быть другие характеристики, определяющие адекватную модельную систему.
Преимущество данного метода исследования устойчивости неконсервативных динамических систем заключается, при условии построения адекватной модели (конечно же, только при этом условии остальное имеет смысл), в определении интервалов значений параметров (по чувствительности), обеспечивающих устойчивый процесс. По заданным условиям на поверхности наплавляемого тела (подложки) данная задача может рассматриваться как обратная, для определения начальных условий технологического процесса (и даже оптимальных).
Регрессионный анализ в порошковой металлургии. Из-за сложности и даже невозможности выявления аналитическими методами влияния всех факторов на теплокинетические процессы и качественные характеристики продукта в металлургии единственным методом научного исследования является регрессионный анализ, который включает планирование, проведение и математическую обработку результатов эксперимента Теории научного планирования эксперимента и регрессионного анализа достаточно хорошо разработаны, но остаются задачи, требующие решения. В порошковой металлургии к таким задачам относятся выявление наиболее значимых факторов и их взаимодействий и определение компонентного и массового составов шихты для получения материала с требуемыми свойствами. Достижение высоких физико-механических и эксплуатационных свойств современных сплавов возможно лишь при использовании методов получения твёрдых сплавов как гетерогенных материалов, состоящих из зёрен высокотвёрдых тугоплавких соединений (карбидов, реже нитридов или боридов переходных металлов), сцементированных пластичным металлом - связкой. В качестве тугоплавкой фазы твёрдых сплавов наиболее широко используются карбиды вольфрама, титана, хрома, тантала, а в качестве металла - связки - кобальт, никель, железо и их сплавы.
Наиболее эффективным методом получения твёрдых сплавов в последнее время признан разработанный академиком А.Г. Мержановым [3] и его сотрудниками СВС, который характеризуется такой уникальной особенностью, как существование в течение короткого времени реакции высокотемпературной твердожидкой среды, допускающей различные типы воздействий (электрические, химические, термодинамические и т. д.). Кроме этих воздействий, возможно использование в СВС-си-стемах различных легирующих и инертных добавок, позволяющих изменять режим синтеза для получения заданных структуры и свойств конечного продукта. СВС из-за высоких температур, необходимых для синтеза, наиболее эффективен в хими-
ческих системах, образующих тугоплавкие материалы.
Одной из основных задач СВС материала с заданными свойствами является определение компонентного массового состава шихты для обеспечения этих свойств.
Для инструментальной стали основными свойствами являются твёрдость (У^, жаропрочность (У,) и вязкость (У3). Эти свойства зависят от того, какие легирующие компоненты представлены определёнными металлами и в каком количестве. И вот здесь представляется возможность определения компонентного состава шихты по результатам физического эксперимента, используя регрессионный анализ. К примеру, к термитной смеси добавлен карбид титана (Д), а легирующими компонентами являются вольфрам (Х2) и хром (Х3). Пусть по результатам факторного эксперимента получены регрессионные уравнения:
у=МХМ), У/ХМ), У^ХМ).
Задавая требуемые свойства У=У/*=соп81 (/=1,2,3), получим замкнутую систему алгебраических уравнений
У*=/(ХМ), /=1,2,3, решение которой Х=Х* (/=1,2,3) определяет массовый компонентный состав шихты для получения интерметаллического сплава с заданными свойствами.
Проведённые в Алтайском государственном техническом университете им. И.И. Ползунова опыты частично подтвердили правильность теоретических расчетов. Так, максимальная температура горения термитной смеси (^1300 °С) и возникновение СВС-процесса (горение Т1С) соответствовали указанным теоретическим значениям. СВС-процесс не возникал при Х1=28 мас. % термитной смеси. С увеличением количества ТЮ от 5 до 28 мас. % термитной смеси скорость фронта горения уменьшалась от 0,15 до 0,07 м/с, его максимальная температура в определенном сечении (^50 мм от верхнего конца трубки, шихта насыпалась в вертикальную трубку из кварцевого стекла диаметром 22 мм и длиной 100 мм при небольшом уплотнении) снизилась от 1290 до 1170 °С при повышении твердости получаемого слитка интерметаллида от 15 до 60 НЯС. Опыты проводили со стандартными порошками
Для исследования взаимосвязей различных воздействий на теплофизические процессы синтеза возможно применение регистраторов динамических оптических полей - пирометров, выполненных на дискретных и многоэлементных фотоприемниках, работающих в режиме прямого детектирования или в режиме накопления заряда. Главной особенностью таких устройств является бесконтактность и малая инерционность измерений, в отличие от контактных температурных датчиков, что способствует проведению тепловизионного контроля быстропротекающих высокоэнергетических процессов с повышенной точностью.
По результатам экспериментов и характеру расположения точек построены зависимости твердости сплава (в НЯС) и скорости распространения фронта горения (в м/с) от количественного содержания карбида титана в шихте (в мас. % термита) (рисунок). Эти зависимости апроксимируются следующими кубическими параболами:
7=10-4,1Х+1,25Х2-0,046Х3, К=0,15+0,0093Х+0,0008Х2-0,0000267Х3.
О 4 8 12 16 20
"ПС, мас. % термита
Рисунок. Зависимость твердости сплава стали У и скорости фронта горения V от массы карбида титана
С увеличением содержания карбида титана в шихте твердость сплава увеличивается, а скорость фронта горения уменьшается по нелинейному закону, причем изменение характера зависимостей происходит практически при одном и том же значении массы карбида титана в шихте (8...10 мас. % термита).
В качестве примера построения дифференциально-регрессионной модели движения твердого тела в сопротивляющейся среде рассмотрим движение плуга при рабочем процессе плуга - пахоте -на основе моделирования взаимодействия рабочего органа плуга (корпуса) и почвы. Теоретическое построение математических моделей сельскохозяйственных агрегатов, их рабочих процессов и систем управления является весьма сложной задачей, практически осуществимой в отдельных случаях моделирования, и то при идеализации условий работы этих объектов и самих объектов.
Исходя из того, что дифференциальные уравнения и функции регрессий наиболее полно моделируют стохастическую техническую систему, предлагается рассматривать движение плуга по дифференциальным уравнениям движения механической системы в пространстве, в которых силовые характеристики представлены уравнениями множественной регрессии на физических переменных и переменных рабочего процесса (или физических координат и переменных параметрах). Разумеется, уравнения регрессии должны быть функциональными.
При исследовании устойчивости движения плуга как динамической системы с переменными параметрами прежде всего следует определить характер движения плуга как твердого тела. Если рассматривать
его движение как поступательное, то получаем прямолинейную траекторию движения центра масс (с учетом сцепки с трактором). Очевидно, для исследования технологической устойчивости движения плуга по глубине вспашки необходимо рассмотреть плоскопараллельное движение в продольно-вертикальной плоскости. Так как одним из агротехнических требований к пахоте является прямолинейность борозды, то необходимо рассмотрение движения в горизонтальной плоскости. Боковые отклонения корпуса плуга также влияют на глубину вспашки и прямолинейность борозды, поэтому необходимо учесть эти движения в поперечно-вертикальной плоскости.
Следующим этапом является определение параметров плуга и технологических параметров, варьирование которых изменяет движение плуга. В общем случае такими параметрами являются: / -сила тяги двигателя, к - ширина плуга, а и Ь - глубина вспашки и ширина захвата корпуса плуга, е и I — угол наклона затылочной фаски лемеха ко дну борозды и ее длина, в — угол наклона лезвия лемеха ко дну борозды, V - скорость плуга, 1„ /у, - моменты инерции плуга относительно осей координат, проходящих через центр тяжести плуга, Т -твердость почвы (обобщенный показатель физико-механических свойств почвы), / - коэффициент трения корпуса плуга о дно и боковые поверхности борозды. Сюда же могут быть включены обобщенные характеристики лемешно-отвальных поверхностей, формы стойки и полевой доски, углы установки стойки корпуса плуга, характеристики установки и формы колес плуга. При этом сила тяги двигателя в общем случае не горизонтальна.
И последнее - построение уравнений движения. Обязательным элементом этих уравнений должен быть учет влияния вращательных моментов на движение плуга. Несводимость сопротивлений почвы, преодолеваемых плужным корпусом, к одной равнодействующей силе является следствием несимметричности и переменной кривизны ле-мешно-отвальной поверхности, непостоянства давления пласта почвы, а также - величины и направления сил трения в различных точках этой поверхности и затупления лезвия лемеха.
Как известно, в общем случае движение твердого тела характеризуется следующей системой уравнений:
шх = ¥х, ту = ¥у, ш2 = ¥г,
3= Мх,, 3у,а>у, = Му,, 3 а2, = М2,, (10)
где оси х и I определяют продольно-вертикальную плоскость; ось х направлена по горизонтали в сто-
СПИСОК ЛИТЕРАТУРЫ
1. Шенк Х. Теория инженерного эксперимента. - М.: Мир, 1982. - 368 с.
2. Торбунов С.С Математическое моделирование и параметрическая устойчивость динамических систем с детерминированными параметрами. - Томск: Изд-во ТПУ, 2002. - 231 с.
рону движения плуга; т - масса плуга; ¥„ /у, ^ — проекции всех сил, включая реакции почвы, на оси координат; /х, /у-, — моменты инерции плуга относительно осей координат, проходящих через центр тяжести плуга параллельно осям х,у,х, х ,у X - ускорения центра тяжести плуга вдоль осей координат; М„ Му, М1 — моменты всех сил относительно осей координат, проходящих через центр тяжести плуга, точки означают производные по времени.
В процессе эксплуатации почвообрабатывающих орудий на лемехи со стороны почвы действует пространственная система сил. Величина внешних сил и характер их распределения по поверхности лемеха зависят в общем случае от физико-механических и технологических свойств обрабатываемой почвы, от конструктивных размеров лемеха, от углов установки лемеха к дну и стенке борозды, от скорости движения агрегата и глубины обработки и от ряда других показателей.
Экспериментальные исследования проводились на пахотных почвах учебно-опытного хозяйства Новосибирского СХИ и в лаборатории СИ-БИМЭ РАСХН (г. Краснообск Новосибирской области) на почвенном канале.
В результате экспериментальной работы установлены регрессионные зависимости влияния длины и ширины лемеха, скорости движения орудия, глубины обработки почвы, ширины и угла наклона затылочной фаски лезвия на величины пространственных силовых характеристик внешних сил и исследована технологическая устойчивость движения плуга. Под технологической устойчивостью движения плуга понимается сохранение прямолинейного поступательного движения, по возможности горизонтального (то есть на склонах по агротехническим требованиям пахоту необходимо производить поперек наклона поля) при постоянной глубине вспашки. Однако в действительности глубина хода рабочих органов машины нарушается из-за непостоянства свойств почвы по длине борозды и неровностей рельефа поля.
В [2] для заданных технологических условий пахоты получены зависимости, определяющие область значений конструктивных параметров, которые обеспечивают технологическую устойчивость движения плуга, то есть качественную пахоту, что доказывает правомерность исследования движения плуга, как динамической системы с детерминированными переменными параметрами, по регрессионной модели.
Аналогично может быть исследовано движение любого твердого тела.
3. Мержанов А.Г. Твердо-пламенное горение. - Черноголовка:
ИСМАН, 2000. - 224 с.
Поступила 14.11.2006 г.