Файл: Математическое моделирование в экологии.doc

Определим с помощью П -теоремы критерий подобия для электрической цепи, которая была рассмотрена в примере 2.3. Цепь содержит активное сопротивление R и индуктивность L. При включении цепи на постоянное напряжение U в ней протекает процесс, который можно записать в виде функциональной зависимости

U = f (t,i,R,L).

Допустим, что определяющие параметры t,i,R имеют независимые размерности

[t] = T, [i] = I , [R] = L²MT^-3I^-2.

Тогда параметры U и L будут размерно-зависимыми.

[U] = L²MT^-3I^-1, [L] = L²MT^-2 I^-2 .

В данном случае k = 3, n - k =1. Анализируя размерности параметров, находим критерии подобия.

Заметим, что в примерах 2.3 и 2.4 критерии подобия найдены на основе общих функциональных зависимостей между определяющими и определяемыми параметрами. Точные уравнения физических зависимостей нами не рассматривались. Независимость физической закономерности (2.14) от выбора

единиц измерения означает, что эту зависимость можно представить

в виде уравнения

Уменьшение числа аргументов упрощает исследование. Пусть в уравнении (2.14) для выяснения зависимости величины a от некоторого определяющего параметра a_iнадо измерить эту величину при десяти значениях данного аргумента. Тогда для экспериментального нахождения величины a - как функции n определяющих параметров a₁,…, a_n - следует произвести 10ⁿ экспериментов.

Согласно П - теореме, если все величины a₁,a₂ ,…,a_n выражаются через k независимых размерностей, задача сводится к определению функции n - k безразмерных аргументов П₁, …, П_n_-_k, для нахождения которых достаточно 10ⁿ^-^k опытов, т.е. в 10^k раз меньше. Трудоемкость установления искомой функции сокращается на столько порядков, сколько среди определяющих параметров величин с независимыми размерностями.

В примере 2.4 исходное уравнение можно заменить зависимостью

П = Ф (П₁).

Следовательно, в координатах П, П₁все опытные точки должны располагаться на единой кривой. Таким образом, заранее проведенный анализ размерностей сокращает объем экспериментальной работы во много раз.

Пример 2.5

Пусть известно, что период T колебаний маятника не зависит от его начального отклонения и скорости, а определяется лишь его длиной l, массой m и ускорением свободного падения g . Функциональная связь

T = f (l,m, g )

содержит четыре размерных величины. Допустим, что определяющие параметры l, m и g имеют независимые размерности.

[l] = L, [m] = M, [g] = LT ^-2.

Тогда параметр T следует принять размерно-зависимым.

[T ] = T.

Анализируя размерности, находим критерий подобия

Здесь, k = n = 3. Следовательно, выражение f (l,m, g ) = const и

С точностью до безразмерного множителя данная формула дает решение колебаний маятника. Причем выясняется, что их период не зависит от m. Таким образом, анализ размерностей дает возможность получить ценную предварительную информацию о процессе, не находя полностью решения задачи.

Пример 2.6

При атомном взрыве в его центре мгновенно выделяется значительная энергия E . От центра взрыва распространяется мощная ударная волна. Радиус фронта ударной волны r через промежуток времени τ после взрыва зависит от E, τ и плотности воздуха ρ [11].

r = f (E,τ,ρ).

Таким образом, n = 3, а размерности определяющих параметров в классе MLT есть:

В данном случае k = n = 3, так что функция f (E,τ,ρ) = const. Критерий подобия

Эта формула показывает, что если измерить тем или иным способом радиус ударной волны в разные моменты времени, то в логарифмических координатах (5/2)lg(r), lg(τ) экспериментальные точки должны располагаться на прямой

имеющей наклон, равный единице. Более детальный анализ показывает, что эта const =1. В 1950 г. английский физик Джеффри Ингрем Тейлор обработал

кинофильм о распространении огненного шара, снятый во время американских испытаний ядерного взрыва в Нью-Мехико в 1945 г., и определил, что энергия взрыва равнялась примерно 10¹⁷Дж.

Публикация Тейлором этой величины вызвала в свое время смущение в американских официальных кругах. В различных уравнениях математической физики (задачах гидродинамики, переноса тепла, диффузии, химической кинетики) появляются своеобразные частные решения, которые при изменении времени преобразуются одно из другого по правилу подобия. Такие решения называются автомодельными. Наиболее важное свойство автомодельного решения состоит в том, что зависимость от аргументов входит через единственный комплекс. Благодаря этому уравнение в частных производных можно привести к обыкновенному дифференциальному уравнению. Его интегрирование существенно проще, чем нахождение решения начально-граничной задачи. В сложных нелинейных задачах получение автомодельного решения зачастую остается единственно возможным средством найти аналитическую зависимость и понять качественные особенности явления. Более подробные сведения об автомодельных решениях можно найти в литературе [2], [9].

11.3.1. Постановка задачи

Дифференциальными называются уравнения, содержащие одну или несколько производных. Инженеру очень часто приходится сталкиваться с ними при разработке новых изделий или технологических процессов, так как большая часть законов физики формулируется именно в виде дифференциальных равнений. Любая задача проектирования, связанная с расчетом потоков энергии или движением тел, в конечном счете, сводится к решению дифференциальных уравнений. Лишь очень немногие из них удается решить без помощи вычислительных машин. Поэтому численные методы решения дифференциальных уравнений играют такую важную роль в практике инженерных расчетов и в моделировании.

Так при реализации цифровых систем управления инженеру-системотех-нику приходится решать дифференциальные уравнения в реальном масштабе времени, т.е. непосредственно в процессе управления объектом. Примером могут служить цифровые регуляторы в системах управления электроприводами металлорежущих станков и промышленных роботов, а также цифровые системы управления автомобильными и авиационными двигателями, летательными аппаратами, морскими судами и т.д.

Рис. 11.3.1. Решения дифференциального уравнения

Известные математические программы – MathCAD, Matlab, Mathematica и др. непригодны для решения таких задач. Эти программы занимают в ЭВМ много памяти и, кроме того, они не могут работать в реальном масштабе времени. Поэтому для построения компактных, работающих в реальном времени цифровых моделей и систем инженеру приходится самостоятельно разрабатывать алгоритмы и программы для решения дифференциальных уравнений тем или иным численным методом.

В зависимости от числа независимых переменных и, следовательно, типа входящих в них производных дифференциальные уравнения делятся на две различные категории:

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

Рассмотрим методы решения обыкновенных дифференциальных уравнений (ОДУ). Дифференциальное уравнение первого порядка можно записать в виде

y' = f(y,t).

Это уравнение имеет семейство решений y(t) . Например, если f(y,t) = y , то для произвольной константы С функция y(t) = Ce^t является решением (рис.3.1). Выбор начального значения, скажем y(0) , служит для выделения одной кривой из семейства кривых. Зачастую имеется более чем одна зависимая переменная, и тогда задача заключается в решении системы уравнений первого порядка, например,

Решение этой системы содержит две постоянные интегрирования, и, следовательно, нужны два начальных условия, чтобы определить эти константы.

Если значения y и z указаны при одном и том же значении независимой переменной t₀, то система будет иметь единственное решение. Задача определения y и z для будущих значений t > t₀ называется задачей с начальными условиями или задачей Коши. Если же условия задаются при двух или более значениях независимой переменной, то задача называется краевой. В задаче Коши дополнительные условия называют начальными, а в краевой задаче – граничными.

Любое обыкновенное дифференциальное уравнение порядка n, которое можно записать так, что его левая часть есть производная наивысшего порядка, а в правой части эта производная не встречается, может быть записана из n уравнений первого порядка путем введения n-1 новых переменных.

Например, уравнение

u'' = g(u,u' t)

можно записать как систему

где z'(t) = u''(t).

При обсуждении методов для задачи Коши удобно представлять

себе единственное уравнение

y '= f(y,t)

с начальным условием y(t₀ ) = y₀. Однако методы с равным успехом применимы и к системам уравнений. Часто в задаче Коши в роли независимой переменной выступает время.