Абстрактное моделирование

Содержание

Слайд 2

Цель лекции. Изложение основных понятий и методов моделирования технических процессов. Этапы

Цель лекции.

Изложение основных понятий и методов моделирования технических процессов.
Этапы абстрактного моделирования.
Классификация

моделей по признаку наличия предварительной информации.
Методы построения моделей детерминированного и стохастического типа.
Модели с использованием марковских цепей и иммитационные модели.
Слайд 3

Физическая модель – это копия прибора, приспособления или машины, называемой натуральной,

Физическая модель – это копия прибора, приспособления или машины, называемой натуральной,

которая подчиняется определенным правилам.
Цель моделирования – изучение определенной ситуации для выбора оптимального решения, удовлетворяющего какому-то критерию.
Слайд 4

Этапы абстрактного моделирования Информационная модель Логико-математическая модель Алгоритм и программа Исследование

Этапы абстрактного моделирования

Информационная модель
Логико-математическая модель
Алгоритм и программа
Исследование на ЭВМ
Интерпретация результатов, проверка

и, если необходимо, корректировка моделей и программ
Слайд 5

– это конкретное словесное описание ситуации, изучаемого явления, где отвечается на

– это конкретное словесное описание ситуации, изучаемого явления, где отвечается на

вопросы: что происходит, почему происходит и при каких условиях происходит (при этом очень важно уточнить начальные и граничные условия ).
Информационная модель переводится на логико-математический язык .

Информационная модель

Слайд 6

Перед построением математических моделей необходимо составить информационную модель. Классификация моделей по

Перед построением математических моделей необходимо составить информационную модель.
Классификация моделей по признаку

наличия информации:
Первая группа - модели, в которых информация отсутствует .
Вторая группа - модели, в которых информация достаточна.
Слайд 7

Классификация моделей: детерминированные, к которым относятся статические (алгебраические) и динамические, представляемые

Классификация моделей:

детерминированные, к которым относятся статические (алгебраические) и динамические, представляемые в виде

систем дифференциальных уравнений;
стохастические (вероятностные), когда имеем дело с массовыми явлениями.
Слайд 8

Если расчет на модели дает результаты такого же порядка, как и

Если расчет на модели дает результаты такого же порядка, как и

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

Детерминированные процессы: Закон сохранения массы: dV=(Q1-Q2).dt при плотности Q1=Q2=const, dV –

Детерминированные процессы:

Закон сохранения массы:
dV=(Q1-Q2).dt
при плотности Q1=Q2=const, dV – изменение объема жидкости

в резервуаре, Q1 и Q2 – входящий и исходящий расход жидкости за время dt.
Слайд 10

Детерминированные процессы: 2. Закон сохранения импульса силы или количества движения: d(mV)=(ΣF).dt,

Детерминированные процессы:

2. Закон сохранения импульса силы или количества движения:
d(mV)=(ΣF).dt,
m –

масса тела, V – его скорость, ΣF – внешние силы и dt – элементарный промежуток времени.
Слайд 11

Детерминированные процессы: 3. Закон сохранения момента количества движения. Этот закон, как

Детерминированные процессы:

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

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

Детерминированные процессы: 4. Закон сохранения энергии: гидромеханика: P – давление, γ

Детерминированные процессы:

4. Закон сохранения энергии:
гидромеханика:
P – давление, γ - объемный вес,

z - высота, V - скорость потока, g - ускорение силы тяжести.
Слайд 13

Детерминированные процессы: Термодинамика Эп+Эк=U=const сумма потенциальной, кинетической и внутренней энергии есть величина постоянная.

Детерминированные процессы:

Термодинамика
Эп+Эк=U=const
сумма потенциальной, кинетической и внутренней энергии есть величина постоянная.

Слайд 14

Детерминированные процессы: явления переноса: диффузия; внутреннее трение или вязкость; теплопроводность.

Детерминированные процессы:

явления переноса:
диффузия;
внутреннее трение или вязкость;
теплопроводность.

Слайд 15

Распространение молекул примеси в жидкости или газе при отсутствии макроскопических перемещений

Распространение молекул примеси в жидкости или газе при отсутствии макроскопических перемещений

подчиняется закону Фука
где q — поток диффундирующего вещества; D — коэффициент диффузии;
ΔC/ΔY — градиент концентрации вещества.
Слайд 16

Сила внутреннего трения равна где η - коэффициент вязкости Через газ,

Сила внутреннего трения равна
где η - коэффициент вязкости
Через газ, заключенный

между параллельными стенкам имеющими разные температуры Т1 и Т2 в направлении, нормальном к стенке, будет распространятся поток тепла, описываемый законом Фурье:
где λ - коэффициент теплопроводности, ΔТ — приращение температуры.
Слайд 17

Многие модели в электротехнике строятся на основании закона Ома для постоянного

Многие модели в электротехнике строятся на основании закона Ома для постоянного

тока:
где i — сила тока, Е — электродвижущая сила, г — внутреннее сопротивление источника, R — нагрузка.
Для переменного тока при наличии индуктивности L и емкости С закон Ома представим в виде:
Слайд 18

Так как Q=V.C и dQ=i.dt, имеем dV/dt = i/C. Подставляя dV/dt

Так как Q=V.C и dQ=i.dt, имеем dV/dt = i/C. Подставляя dV/dt

в уравнение и дифференцируя его, получим:
Слайд 19

На основании первого закона Кирхгофа для сходящихся в точке токов имеем

На основании первого закона Кирхгофа для сходящихся в точке токов имеем

 Σi=0 - как выражение закона сохранения электрического заряда.
По второму закону Кирхгофа для замкнутого контура имеем
ΣiR+ Σir= ΣE,
т.е. алгебраическая сумма произведений токов на сопротивление равна алгебраической сумме электродвижущих сил, действующих в замкнутом контуре.
Слайд 20

В физике имеется группа уравнений в частных производных, известных под названием

В физике имеется группа уравнений в частных производных, известных под названием

уравнений математической физики. Уравнение теплопроводности или диффузии представим в общем виде:
Слайд 21

Явление выделения Пусть за время dt выделяется масса dm, определяемая выражением

Явление выделения
Пусть за время dt выделяется масса dm, определяемая выражением
-dm= p·m·dt,
где

р - коэффициент пропорциональности, (1/с). Тогда
dm/m=-p·dt
Интегрируя это выражение, получаем:
lnm = -pt + C
при t0=O, C·ln m0. Потенцируя выражение получим m=m0 ·e-pt . Если обозначить р через 1/t0, то m=m0exp(-t/t0). В безразмерном виде m=m0ехр (- τ ), где τ=t/t0 или μ=ехр(- τ).
Слайд 22

Уравнение движения поплавка Будем исходить из уравнения баланса массы. Изменение объема

Уравнение движения поплавка

Будем исходить из уравнения баланса массы. Изменение объема

воды в резервуаре за время dt равно разнице между входом и выходом, что можно представить в виде уравнения
S·dH = (Q - q) · dt,
где Q - приток воды, q - расход воды, S - площадь сечения резервуара, Н - высота.
Слайд 23

Для короткого трубопровода справедливо соотношение: где f(y), [м2] – функция высоты

Для короткого трубопровода справедливо соотношение:

где f(y), [м2] – функция высоты Н.

Можно предположить, что f(y)=B·y в случае прямоугольного сечения отводящей трубы. При μ=const:

В установившемся режиме:

Слайд 24

Динамические модели на основе модели поплавка Учет колебаний воды в трубопроводе

Динамические модели на основе модели поплавка

Учет колебаний воды в трубопроводе

приводит к сложному явлению «гидравлического удара».
На попловок действует сила веса G, подъемная сила
Fn=γвSn(H-y),
А также сила вязкости сопротивления
Слайд 25

где m – масса поплавка, где a=γВS/m и 1/T=B/m.

где m – масса поплавка,

где a=γВS/m и 1/T=B/m.

Слайд 26

при dE/dt ≠ 0 получаем типичный вынужденный колебательный процесс Решение этого

при dE/dt ≠ 0 получаем типичный вынужденный колебательный процесс

Решение этого

типа уравнений имеет вид затухающего "скачка" или затухающих колебаний при большом сопротивлении и dE/dt=0

при dE/dt ≠ 0 получаем типичный вынужденный колебательный процесс

где
φ - угол поворота вала; Т - постоянная, учитыва­ющая момент инерции масс (с2);
U - напряжение.

Слайд 27

Вероятностные (стохастические) модели В стохастических моделях участвуют случайные величины.

Вероятностные (стохастические) модели

В стохастических моделях участвуют случайные величины.

Слайд 28

рассмотрим процесс осаждения группы тяжелых частиц в канале при распределенной скорости

рассмотрим процесс осаждения группы тяжелых частиц в канале при распределенной скорости

воды, являющейся функцией глубины канала.

частицы в канале имеют нулевую скорость в точке Х=0 и h = h0 =2м
скорость воды в канале соответствует v=(1-0,25h) · h
при h=h max=2m максимальная скорость составляет величину vmax =1м/с

Слайд 29

частицы испытывают воздействия турбулентных пульсаций скорости потока fx и fy частицы

частицы испытывают воздействия турбулентных пульсаций скорости потока fx и fy

частицы

имеют скорость опускания U= f (d)

пульсации имеют нормальное распределение со средним значением
f х= f y= 0 и средним квадратическим отклонением Se=qv, где q - коэффициент пропорциональности.

Таким образом, величина пульсации пропорциональна скорости потока.

Слайд 30

вычисляют два случайных числа, равномерно распределенных в пределах 0≤ р ≤1

вычисляют два случайных числа, равномерно распределенных в пределах 0≤ р ≤1


вычисляют уже нормально распределенные числа при Рср=0 и σ2=1, где
σ2- дисперсия,
Рср -среднеквадратичное отклонение

Слайд 31

скорости перемещения частиц будут следующими: по высоте (вниз) — (u+fv) по

скорости перемещения частиц будут следующими:
по высоте (вниз) — (u+fv)
по

длине — (v+fx).
За время d t имеем перемещения частиц по высоте и длине соответственно:
h = h- (u + fv) ·dt,
Х = Х+ (V+ fx ) ·dt.
Опускание частицы заканчивается при h=0. При этом отмечается расстояние х опускания частиц.
Слайд 32

Численный эксперимент состоит в последовательном опускании частиц и регистрации расстояния опускания

Численный эксперимент состоит в последовательном опускании частиц и регистрации расстояния опускания

Lx. Целью эксперимента может быть, например, изучение распределения частиц по дну канала в зависимости от глубины, скорости воды и ее градиента по высоте, а также анализ влияния крупности частиц или пульсаций скорости потока воды на распределение частиц по высоте канала до момента опускания их на дно.
Слайд 33

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

Пример построения математической модели с последующим усложнением

Необходимо построить математическую

модель смесителя в виде резервуара, в который поступает тяжелая жидкость с плотностью -
ρ1 (кг/м3) и расходом Q1 ( м3/ с). Растворителем служит вода с плотностью -
ρ2 =1000 (кг/м3) и расходом Q2 (м3/ с) соответственно. Из резервуара вытекает смесь с расходом Q3 (м3/ с) при заданной плотности ρ3 (кг/м3).
Слайд 34

При равновесии систем имеем уравнение баланса количества жидкости Q1+Q2=Q3 а также

При равновесии систем имеем уравнение баланса количества жидкости Q1+Q2=Q3

а также

уравнение баланса массы
ρ1Q1+ ρ 2Q2= ρ 3Q3
откуда
ρ 3=(Q1 ρ 1+ Q2 ρ 2)/(Q1+Q2)
Слайд 35

Например, ς3=1500 кг/м3 при ς2= -1800 кг/м3, откуда Q2=Q1·(ρ1- ρ3)/(ρ3- ρ2)=0,6·Q1

Например, ς3=1500 кг/м3 при ς2= -1800 кг/м3, откуда

Q2=Q1·(ρ1- ρ3)/(ρ3- ρ2)=0,6·Q1

Слайд 36

В динамике также будем исходить из уравнения сохранения массы - массового

В динамике также будем исходить из уравнения сохранения массы - массового

баланса за время d t.

Изменение массы в системе:
dm = (Q1 ρ1+ Q2 ρ2- Q3 ρ)·dt
однако
dm = d (V ρ) = ρ ·dV+V·d ρ
dV = d(Fh) = h·dF+F·dh

Слайд 37

Введем первое ограничение, что резервуар цилиндрический, и тогда dV=F·dh и dm=

Введем первое ограничение, что резервуар цилиндрический, и тогда dV=F·dh и dm=

ρ F · dh + V ·d ρ.

Пусть резервуар заполнен водой до уровня, при котором
Q3=Q1+Q2
следовательно

dV=0, так как V=const

Слайд 38

Тогда Vd ρ = (Q1 ρ1+Q2 ρ2-Q3 ρ) dt получаем Решение

Тогда Vd ρ = (Q1 ρ1+Q2 ρ2-Q3 ρ) dt

получаем

Решение -1/b ·ln(a-by)=x+c или

–ln(ρ∞-ρ)=t/T+C, где T=V/Q3

И окончательно ρ = (ρ∞-ρ2)·e-t/T(кг/м3)

при t=0 C= - ln(ρ∞-ρ2)

Слайд 39

Рассмотрим более общую модель dV=F·dh+h·dF=(Q1+Q2-Q3) · dt, где При F=const, dV=F·dh

Рассмотрим более общую модель

dV=F·dh+h·dF=(Q1+Q2-Q3) · dt, где

При F=const, dV=F·dh

Слайд 40

Баланс массы

Баланс массы

Слайд 41

Подставим вместо dv/dt

Подставим вместо dv/dt

Слайд 42

Марковские цепи Для описания многих явлений, которые можно представить как совокупность

Марковские цепи

Для описания многих явлений, которые можно представить как совокупность

ряда состояний, в современной вычислительной математике находит применение теория цепей Маркова. Очень эффективно эта теория используется в случае явлений, описываемых моделями с распределенными параметрами.
Слайд 43

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

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

Частица,

находящаяся в i -м слое, в результате случайных блужданий может перейти за время d t в прилегающие соседние слои, один из которых расположен выше, а другой – ниже данного слоя с вероятностью соответственно Pi,i-1 и Pi,i+1 или оставаться в данном слое с вероятностью Pi,i. Первый индекс вероятностей указывает, из какого слоя переходит данная частица, а второй — в какой слой.
Слайд 44

Частицы могут переходить и на более удалённые слои, то есть Pi,i-j , …, Pi,i+1.

Частицы могут переходить и на более удалённые слои, то есть Pi,i-j

, …, Pi,i+1.
Слайд 45

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

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

в

совокупность состояний для момента времени

при помощи умножения вектора аt на матрицу Р в виде

Слайд 46

Например, при размерности матрицы, равной трем, имеем: Тогда

Например, при размерности матрицы, равной трем, имеем:

Тогда

Слайд 47

Если речь идет о количестве частиц в замкнутом резервуаре, то соблюдается закон сохранения

Если речь идет о количестве частиц в замкнутом резервуаре, то соблюдается

закон сохранения
Слайд 48

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

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

Маркова

уравнения Фоккера-Планка:

где в общем случае A,D = const.

Слайд 49

производные в виде конечных разностей

производные в виде конечных разностей

Слайд 50

Конечно-разностный аналог уравнения

Конечно-разностный аналог уравнения

Слайд 51

Обозначая имеем

Обозначая

имеем

Слайд 52

Обозначим также Рi,i=1-2c, Pi-1,i= - b+c, Pi+1,i= - b+c Так как

Обозначим также Рi,i=1-2c, Pi-1,i= - b+c, Pi+1,i= - b+c Так как

b>0 и с>0, то для устойчивости решения требуется лишь, чтобы Рi,j ≥ 0. Отметим также, что здесь имеет место закон сохранения величины (массы, импульса и т.д.)

Расчет марковской цепи

Слайд 53

Формулы вычислений имеют вид Формы для вычисления шага интегрирования или, точнее, дискретного шага расчета марковской цепи.

Формулы вычислений имеют вид

Формы для вычисления шага интегрирования или, точнее, дискретного

шага расчета марковской цепи.
Слайд 54

Имитационное моделирование Имитационное моделирование является перспективным направлением моделирования явлений и процессов

Имитационное моделирование

Имитационное моделирование является перспективным направлением моделирования явлений и процессов в

при­роде и технике. Оно возникло с появлением ЭВМ и получит еще более широкое применение по мере развития вы­числительной техники.
Слайд 55

рассмотрим модель, в которой волк преследует зайца. Пусть в момент времени

рассмотрим модель, в которой волк преследует зайца. Пусть в момент времени

t=t0 заяц находится в точке x=xz, у = 0, а волк — в точке х=хv, у = уy. Заяц может перемещаться лишь вдоль оси х с постоянной скоростью Vz.
Слайд 56

Составим - дифференциальное уравнение, описывающее движение волка. Пусть за время dt

Составим - дифференциальное уравнение, описывающее движение волка. Пусть за время dt

заяц переместился на dxz=Vz dt, а волк — на Vvdt. В проекции на оси координат имеем dx = Vvdt·sinΘ, -dy = Vvdt·cosΘ,

где

Отсюда

Слайд 57

Так как при t=0 заяц находится в начале координат, то xz=Vz··

Так как при t=0 заяц находится в начале координат, то xz=Vz··

t. Таким образом, в качестве модели преследования имеем систему нелинейных дифференциальных уравнений: