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

matmod

.pdf
Скачиваний:
11
Добавлен:
21.05.2015
Размер:
898.18 Кб
Скачать

Л.2. Типы математических моделей

Пример 1 (продолжение).

Математическая (структурная) модель рассматриваемой задачи:

x1 = 0, x˙1 = v,

 

 

x2 = 0,

 

2 = 0 (t = 0),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m11 + k1x1 = 0,

 

 

 

m2(¨x1 + x¨2) + k2x2 = 0 (0 ≤ t ≤ T ).

Решение системы ОДУ:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1 = ω1

 

sin ω1t

 

 

 

 

ω1 = r

 

m1

 

,

T =

ω1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

k1

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

sin ω2t

 

 

sin ω1t

 

 

 

 

 

r

k2

 

 

 

 

 

 

 

ω12

 

− ω22

 

 

ω2

 

 

 

ω1

 

 

 

 

 

 

m2

 

 

 

 

x2 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω2 =

 

 

 

 

 

.

 

 

 

 

При t = T :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω12

 

 

v

 

 

 

ω2

π, x˙

2(T ) =

 

 

ω12

 

 

1 + sin

 

 

ω2

 

 

 

ω12 − ω22

 

ω2

 

 

ω1

 

ω12 − ω22

 

 

ω1

 

x2(T ) =

 

 

 

sin

 

 

 

 

 

 

 

 

v

 

 

 

 

π .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2014

21 / 74

Л.2. Типы математических моделей

Пример 1 (продолжение).

Из сохранения энергии колебательной системы:

 

 

k2

2

 

 

k2

2

 

 

m2

2

 

 

 

 

v2

 

 

 

ω12

 

 

2

 

 

2 πω2

 

 

 

 

 

A

 

=

 

 

x2

(T ) +

 

(x˙2) (T ) = 2k2

 

 

 

cos

 

 

 

.

 

 

2

 

2

2

ω22

ω12 − ω22

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

ω12

 

 

 

 

 

πω2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω2

 

ω12

− ω22

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A = hv,

h =

 

 

 

 

 

 

cos

 

.

 

 

 

 

 

 

 

 

 

 

Для m1 = 10

4

кг, m2 = 200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

Н/м:

 

 

кг, k1 = 9.8 10

 

Н/м, k2 = 1.96

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×

 

5

 

 

 

 

 

 

 

 

 

×

 

5

 

 

h

 

 

 

 

 

 

−3

с,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1.79 × 10

 

 

3v (м/с).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A (м)= 1.79 × 10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последнюю формулу можно рассматривать как функциональную

 

 

модель.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ее можно получить методом идентификации математической

 

 

 

модели: определением параметров после выбора схемы модели.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2014

 

22 / 74

Л.2. Типы математических моделей

Дискретные и непрерывные модели

Модели – как содержательные, так и математические – могут быть

дискретными,

непрерывными,

смешанными.

Между этими типами нет принципиального барьера и при уточнении или видоизменении модели дискретная картина может стать непрерывной и обратно; то же может произойти в процессе решения математической задачи

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

2014

23 / 74

Л.2. Типы математических моделей

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

Замена распределенного груза сосредоточенным

Пусть q(x) – плотность распределения груза на малом интервале (l), расположенном вблизи x = a.

R

Q = (l) q(x)dx – груз, сосредоточенный в точке.

q(x) ↔ q˜(x) = Qδ(x − a)

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

целесообразен, если конкретный вид функции q(x) нам неизвестен, но суммарное значение Q мы знаем.

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

2014

24 / 74

Л.2. Типы математических моделей

Пример 2 (продолжение).

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

Эта плотность получается с помощью осреднения исходного распределения:

q(x) = q(x; x) = 1x Q(x; x).

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

2014 25 / 74

Л.2. Типы математических моделей

Аналогичная процедура осреднения происходит при переходе

от среды, состоящей из отдельных частиц к сплошной среде, параметры которой (плотность, температура и т. п.) распределены по пространству;

от поезда с дискретными вагонами к его непрерывной модели;

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

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

2014

26 / 74

Л.2. Типы математических моделей

Внесение возможных неадекватностей при переходе дискретная ↔ непрерывная модель (Пример 3).

Пример 3. Продольные свободные упругие колебания прямолинейного однородного стержня.

u = u(x, t) – смещение в момент t сечения с координатой , отсчитываемой вдоль стержня.

2u

= a

2

2u

,

a

2

= E /ρ

 

 

∂t2

 

∂x2

 

 

 

 

 

 

 

Заменим непрерывную модель на дискретную: последовательность с шагом h материальных точек массы m = ρSh, соединенных пружинами нулевой массы и жесткости k = SE /h.

2014

27 / 74

Л.2. Типы математических моделей

Пример 3 (продолжение).

ui (t) – смещение i-й точки.

Из второго закона Ньютона получаем систему ОДУ:

m

d2ui

= k(ui+1 − ui ) − k(ui

− ui−1)

dt2

или

 

d2ui

= a2

ui+1 − 2ui + ui−1

 

 

dt2

 

 

 

 

 

h2

 

 

Замечание: ∂2u/∂x2 (непр. мод.) → симметричная разделенная разность второго порядка.

2014

28 / 74

Л.2. Типы математических моделей

Пример 3 (продолжение).

Обратный переход от дискретной модели к непрерывной:

ui+1 = u(xi , t),

 

 

 

 

 

 

 

 

h2 ′′

ui±1 = u(xi ± h, t) = u(xi , t) ± hux (xi , t) +

 

 

uxx (xi , t) + . . .

 

2

 

ui+1 − 2ui + ui−1

u′′

(x

, t) +

O

(h) (h

0).

 

h2

 

xx

i

 

 

 

Решение системы ОДУ (дискретная модель) при малом h хорошо имитирует решение уравнения частных производных.

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

2014 29 / 74

Л.2. Типы математических моделей

Пример 3 (продолжение).

Проверка закона сохранения энергии E для стержня конечной длины ℓ, жестко закрепленного на концах (u|x=0,ℓ = 0, ut|x=0,ℓ = 0).

Непрерывная модель.

1

 

 

 

1

 

 

 

 

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z0

 

Z0

 

 

 

Z0

 

 

 

 

 

 

 

 

 

 

 

 

E =

 

 

ρ(ut )2dx +

 

 

 

E (ux )2dx

=

 

 

h(ut )2

+ a2

(ux )2idx

2

2

2

 

Дифференцируем по параметру (t):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂t

 

 

 

 

h

t tt

 

 

 

x xt i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 Z0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂E

=

ρ

 

 

2uu′′ + 2a2uu′′

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя уравнение utt′′ = a2uxx′′ , получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂t

t xx

x xt

 

 

 

 

x t

 

 

 

 

 

 

 

x t

 

 

0

 

 

Z0

 

Z0 ∂x

 

 

 

 

 

 

 

 

 

 

 

 

∂E = ρa2

 

h

 

 

 

 

i

 

 

 

 

uu

 

 

 

 

 

 

 

 

u

 

 

 

uu′′ + uu′′ dx = E

 

 

dx = E u

 

 

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2014

30 / 74

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]