Лекция №8 Построение математических моделей технологических объектов и систем аналитическим методом
Вид материала | Лекция |
- Утверждаю, 89.56kb.
- Программа по курсу " Моделирование систем управления, 30.71kb.
- План изучения дисциплины № п/п, 155.57kb.
- Моделирование и оптимизация систем с распределенными параметрами, 14.61kb.
- Лекция 1 Виды математических моделей сложных систем, 201.52kb.
- Лекции по дисциплине «Математическое моделирование» для студентов и магистрантов специальности, 21.92kb.
- Учение студентов основам математического моделирования, необходимых при проектировании,, 28.67kb.
- Название Лекция-семинар: Построение математических моделей целочисленного линейного, 64.42kb.
- К. А. Тимирязева Экономический факультет Кафедра экономической кибернетики Светлов, 241.83kb.
- 1. Введение Основы анализа данных. Методология построения моделей сложных систем. Модель, 399.94kb.
Лекция № 8
Построение математических моделей технологических объектов и систем аналитическим методом
В основу аналитического метода составления математических моделей положен теоретический анализ конструкции исследуемого объекта (системы) и происходящих в нем физических процессов. Основу данного анализа составляют фундаментальные физические законы, к которым, прежде всего, относятся законы сохранения (массы, энергии, количества движения).
Общая формулировка закона сохранения следующая: изменение во времени некоторой субстанции в элементарном объеме равно сумме притока-стока этой субстанции через его поверхность с учетом скорости генерации или уничтожения субстанции в этом объеме.
Уравнение, соответствующее данной формулировке, имеет вид –

где


G – скорость генерации или уничтожения субстанции.
Дивергенция вектора плотности потока характеризует сумму притока-стока субстанции через поверхность элементарного объема и представляет собой скалярную величину, образуемую суммой частных производных вектора


На основе общей трактовки закона сохранения могут быть получены уравнения для различных субстанций. В частности, уравнение закона сохранения массы имеет вид –

где



В одномерном случае, когда скорость направлена вдоль одной оси декартовой системы координат, плотность потока массы измеряется в кг/(м2с). Примером использования данного закона является описание свойства непрерывности потока жидкости в трубопроводе –

Уравнение закона сохранения энергии может быть представлено в виде –

где




В одномерном случае плотность потока энергии измеряется в Дж/(м2с). Примером использования данного закона является уравнение теплопроводности, связывающее изменение температуры во времени и пространстве со свойствами среды –

где Q – количество тепловой энергии в единице объема, Дж/м3,


Уравнение закона сохранения количества движения для потока идеальной жидкости (без учета сил трения, обусловленных вязкостью) имеет вид –


где


Примером использования данного закона является уравнение Навье-Стокса, приближенная форма которого для анализа движения жидкости в трубопроводе имеет вид –

где



Практика показывает, что при разработке математических моделей чаще всего используются закономерности, носящие характер уравнений баланса. К ним относятся уравнения материального и энергетического баланса, уравнения экономического баланса, а также уравнения равновесия сил.
Уравнение материального баланса основано на законе сохранения массы вещества и может быть записано в следующей форме:
Приход вещества – Расход вещества = Накопление вещества
Разность между приходом и расходом вещества равна изменению его количества в рассматриваемом объекте. В стационарном (установившемся) режиме не может происходить ни убыль, ни накопление. В этом случае уравнение материального баланса для n-го количества веществ может быть записано в виде:

где


С учетом изменения количества вещества в объекте уравнение материального баланса может быть записано в виде:

где


При определении статистических зависимостей можно воспользоваться интегральной записью уравнения материального баланса в следующей форме:

В качестве примера использования уравнения материального баланса рассмотрим получение модели, описывающей смешивающее устройство, в которое поступает два потока веществ, а выходит один смешанный поток.
Материальный баланс в переходном процессе может быть выражен как изменение количества вещества, находящегося в растворе емкости, равное разности между количеством притекающих веществ и количеством вытекающего вещества за одно и то же время:

где Q1, Q2 – потоки поступления в смешивающее устройство 1 и 2 вещества соответственно;
C1, C2 – концентрация веществ 1 и 2 в приходящих потоках;
Qc, Cc – расход вещества из емкости и его концентрация;
V – объем смешивающего устройства.
В установившемся режиме количество вносимого вещества должно быть равно количеству выносимого раствора, то есть:

Уравнения (12) и (13) могут быть записаны и в стандартном виде через «вход – выход»:

где


Не менее редко, чем материальный баланс применяется энергетический баланс, основанный на законе сохранения энергии. В установившемся режиме количество энергии, притекающей в объект, равно количеству энергии, уходящей из него. По аналогии с выражением (9) уравнение энергетического баланса может быть записано в виде:

где

Примером использования данного типа уравнений является уравнение описания процессов в электрическом колебательном контуре, которое получается на основе 2 закона Кирхгофа (Для любого контура электрической цепи с сосредоточенными параметрами алгебраическая сумма напряжений на ветвях контура равна нулю), отражающего по сути энергетический баланс в электрической цепи:

где q(t) – заряд конденсатора в момент времени t.
Аналитические соотношения, основанные на экономическом балансе, описывают показатели эффективности процессов управления производством. Эти уравнения отражают прибыль, себестоимость продукции, приведенные затраты, производительность и другие экономические характеристики производства.
Кроме указанных уравнений баланса для математического описания объектов также используют:
- уравнения элементарных процессов для локальных элементов объектов (уравнения теплообмена, уравнения химических реакций, уравнения напряжений и токов элементов электрических цепей и т.п.);
- ограничения на параметры процесса (например, при моделировании технологических процессов химического производства на концентрации компонентов в многокомпонентных смесях накладывают ограничения по их значениям от 0 до 1).
Все уравнения, описывающих физические процессы в моделируемой системе можно условно разделить на компонентные и топологические. Компонентные уравнения отражают свойства отдельных элементов моделируемой системы и основываются на физических законах протекания элементарных процессов. Примерами компонентных уравнений для различных систем являются:
- уравнение поступательного движения твердого тела, полученное на основе второго закона Ньютона -


- уравнение Эйлера для трубопровода постоянного сечения -



- уравнение резистора, полученное на основе закона Ома -

Топологические уравнения базируются на физических законах, выражающих условия равновесия и непрерывности фазовых переменных и описывают взаимодействия между элементами моделируемой системы через соотношения между однотипными фазовыми переменными. В частности, к топологическим относятся приведенные выше уравнения баланса. Примерами топологических уравнений для различных систем являются:
- уравнение равновесия сил механической системы, основанное на принципе Даламбера (геометрическая сумма всех сил, приложенных к твердому телу, включая силу инерции, равна нулю) -

- балансное уравнение (13), отражающее условие непрерывности потоков жидкости;
- уравнение (17), базирующееся на втором законе Кирхгофа.
Одним из важнейших свойств аналитического метода построения математических моделей является идентичность структуры уравнений, описывающих совершенно разные процессы. В этом проявляется единство физических законов, несмотря на многообразие форм существования материи. Один и тот же тип элемента в системах различных видов описывается аналогичными по структуре компонентными уравнениями.
Моделирование систем автоматического управления осуществляется, как правило, с использованием алгебраических и дифференциальных уравнений. К алгебраическим уравнениям обычно сводится математическое описание стационарных режимов работы объектов. Для описания динамики исследуемых объектов используют дифференциальные уравнения.
Рассмотрим решение уравнения динамики САУ операционным методом при ненулевых начальных условиях в общем виде. Пусть уравнение имеет вид:

где



Уточним понятие начальных условий для уравнения (22). Под начальными условиями понимают совокупность начальных значений искомой функции y(t) и ее производных до (n-1)-й включительно. Различают левые (y(-0), y’(-0) …) и правые (y(+0), y’(+0) …) начальные условия, различие между которыми заключается в направлении подхода к точке t=0. Левые начальные условия (называемые также предначальными) характеризуют систему до начала динамического процесса, создаваемого воздействием x(t). Если в момент времени t=0 входное воздействие и его производные отличны от нуля и происходит мгновенное изменение выходной величины, то начальные условия будут отличны от предначальных. ( При этом количество начальных условий, отличных от предначальных определяется видом передаточной функции объекта W - чем меньше степень числителя W отличается от степени знаменателя, тем большее число начальных условий отличается от предначальных ).
Если заданы предначальные условия, то преобразование уравнения (22) по Лапласу в соответствии с известной теоремой дифференцирования оригинала (


где

Если выражение (24) разложить на суммы и сгруппировать слагаемые относительно множителей «s в степени», то получим следующее выражение для М:

Из уравнения (23) легко найти изображение искомой функции:

Если передаточная функция W имеет полюсы (корни полинома D(s))






Если заданы начальные условия, то осуществляя преобразование по Лапласу уравнения (22) необходимо учитывать как начальные значения искомой переменной и ее производных, так и начальные значения внешнего воздействия и его производных. Тогда после преобразования по Лапласу уравнение (22) примет вид:

где


Начальные значения

Заключительным этапом решения дифференциального уравнения операционным методом является отыскание оригинала (искомой функции времени) по изображению. Для отыскания оригинала по известному изображению в операционном исчислении используются так называемые теоремы разложения и теорема свертывания оригиналов, с которыми вы можете познакомиться в специальной литературе по операционному исчислению.
Исходя из вышеизложенного, методика построения математических моделей объектов аналитическим методом сводится к следующей последовательности действий:
- Физическое описание объекта моделирования. При этом выделяют «элементарные» процессы, протекающие в объекте моделирования, которые подлежат математическому описанию, и формулируют основные допущения, принимаемые при их описании, а также определяются фундаментальные законы, которые должны быть положены в основу составления математических уравнений исследуемого объекта (системы).
- Составление уравнений статики или динамики объекта моделирования. Данный этап включает составление компонентных и топологических уравнений на базе выделенных на первом этапе физических законов в алгебраической или дифференциальной форме.
- Определение начальных и граничных условий моделирования. Данные условия выбираются исходя из особенностей функционирования моделируемого объекта, обусловленных технологическим процессом, в котором он задействован или для которого он разрабатывается.
- Выбор метода решения уравнений математического описания объекта, разработка алгоритма и составление программы. При выборе метода решения системы уравнений обычно руководствуются требованиями обеспечения максимальной быстроты получения решения.
- Проверка соответствия (адекватности) модели объекту. Математическая модель объекта является лишь его определенным в рамках принятых допущений аналогом. Адекватность модели проверяется путем сравнения значений переменных, получаемых на модели и на реальных объектах, и проверки выполнимости критериев адекватности, которые базируются на методах дисперсионного анализа и анализа остатков, существо которых мы с вами рассматривали на прошлых лекциях.
- Изучение свойств объекта моделирования на математической модели. Данное изучение осуществляется в целях определения оптимальных условий протекания процесса, оптимизации управления процессом, а также выработки решений на создание новых объектов.
При проектировании принципиально новых процессов и объектов аналитический метод является единственным приемлемым методом математического описания исследуемых объектов.