matmod
.pdfЛ.2. Типы математических моделей
Пример 3 (продолжение). Проверка закона сохранения энергии E.
Дискретная модель. 0 6 i 6 N, u0(t) = uN (t) = 0.
|
|
|
|
|
|
|
|
|
|
N−1 m |
|
du |
|
2 |
|
N−1 k |
|
|
− ui )2 |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
E = i=1 |
2 dti |
+ i=0 |
|
2 (ui+1 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Дифференцируем по t и применяем систему уравнений для ui (t): |
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
∂E |
|
|
ma2 |
N−1 |
dui |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N−1 |
|
|
|
|
|
|
|
|
|
dui+1 |
|
dui |
|
|||||||
|
|
= |
|
|
(u |
|
|
|
2u |
+ u |
− |
) + k |
|
|
(u |
|
|
|
u |
) |
|
|
|||||||||||||||||||
|
h2 |
i=1 |
|
|
|
i+1 − |
|
i=0 |
i+1 − |
|
|
dt − dt |
|||||||||||||||||||||||||||||
|
∂t |
|
|
dt |
|
|
i |
|
|
i |
1 |
|
|
|
|
|
i |
|
|
||||||||||||||||||||||
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Преобразуем вторую сумму: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
N−1 |
|
|
|
dui+1 |
|
dui |
= |
N−1 dui |
|
|
|
|
|
|
|
N−1 |
|
|
|
|
|
|
dui |
||||||||||||||||||
|
i=0 |
(ui+1 − ui ) |
dt |
|
− dt |
i=1 |
dt (ui − ui−1) − i=1 (ui+1 − ui ) dt |
||||||||||||||||||||||||||||||||||
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|||||||||
Подставляя m = ρSh и k = SE /h, получим ∂E/∂t = 0. |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2014 |
31 / 74 |
Л.2. Типы математических моделей
Пример 3 (продолжение).
Отличия между решениями уравнений в непрерывной и дискретной моделях.
Скорость распространения любых деформаций согласно непрерывной модели равна a
Для дискретной модели влияние внешних воздействий распространяется с формально бесконечной скоростью
Адекватность можно восстановить, если пренебречь слишком малыми значениями смещений: u(xi ) 6 ǫ → u(xi ) ≡ 0.
2014 32 / 74
Л.3. Типы математических моделей
Лекция 3. Типы математических моделей (продолжение)
2014 |
33 / 74 |
Л.3. Типы математических моделей
Линейные и нелинейные модели
Линейная зависимость одной величины от другой – это пропорциональность их приращений, т.е.
z = ax + by + c z = a x + b y
Типичные линейные зависимости между физическими величинами – закон Гука,
закон Ома, закон теплопрового расширения и т.д.
Понятие линейной модели.
Применяется для моделей объектов, рассматриваемых как преобразователи (операторы), для которых каждому входу соответствует некоторый выход.
Пусть начала отсчета входа и выхода выбраны так, что нулевому входу отвечает нулевой выход.
Модель называется линейной, если в ней выполнен принцип суперпозиции (наложения), т. е. при сложении входов складываются и выходы, а при умножении входа на любое число выход умножается на то же число. Если этот принцип не выполнен, модель называется нелинейной.
2014 |
34 / 74 |
Л.3. Типы математических моделей
Линейные модели обычно описываются линейными неоднородными уравнениями – алгебраическими, дифференциальными и т.д., в которых неоднородный член отвечает входу, а решение выходу. Примеры:
1. Краевая задача о малом поперечном прогибе стержня с шарнирно закрепленными концами x = a, x = b:
|
d2 |
EI |
d2y |
= q(x) (a 6 x 6 b), |
|
dx2 |
dx2 |
||
y(a) = 0, |
y′′(a) = 0, y(b) = 0, y′′(b) = 0. |
q(x) – вход, y(x) – выход.
2. Задача Коши для осциллятора под действием внешней силы F (t) (вход). Координата x – выход.
mx¨ + kx = F (t) (0 6 t < ∞), x(0) = 0, x˙ (0) = 0.
2014 |
35 / 74 |
Л.3. Типы математических моделей
Линеаризация
Линеаризация – приближенная замена нелинейной модели на линейную.
Применяется в случаях:
1)эксперимент показывает, что отклонение от линейности в рассматриваемых диапазонах изменения переменных невелико и несущественно,
2)эти диапазоны малы и мы заменяем приращения переменных на их дифференциалы, отбрасывая члены высшего порядка малости.
(Во втором случае применяется также линейное интерполирование.)
2014 |
36 / 74 |
Л.3. Типы математических моделей
Линеаризация. Примеры.
1.Пусть величины х, у, z связаны уравнением
x3 + xy + e2z + 5 = 0.
y3 − xz
Это уравнение удовлетворяется при x = 2, y = −1, z = 0. Пусть ξ, η, ζ – малые изменения x, y, z соответствено:
x = 2 + ξ, y = −1 + η, z = ζ.
Требуется найти линейное соотношение между ξ, η, ζ, справедливое с точностью до членов высшего порядка малости.
Продифференцируем обе части исходного уравнения:
3x2dx + ydx + xdy |
− |
(x3 + xy)(3y2dy − zdx − xdz) |
+ 2e2z dz = 0, |
|
y3 − xz |
(y3 − xz)2 |
|||
|
и подставим x = 2, y = −1, z = 0 и dx = ξ, dy = η, dz = ζ (пренебрегаем величинами высшего порядка малости).
2014 37 / 74
Л.3. Типы математических моделей
|
|
x3 + xy |
|
||
|
|
|
+ e2z + 5 = 0, |
|
|
|
|
|
|
||
3x2dx + ydx + xdy |
|
y3 − xz |
|
||
− |
(x3 + xy)(3y2dy − zdx − xdz) |
+ 2e2z dz = 0, |
|||
y3 − xz |
|||||
|
(y3 − xz)2 |
|
(11ξ + 2η)(−1) − 6(3η − 2ζ) + 2ζ = 0 11ξ + 20η − 14ζ = 0.
Или в переменных x, y, z:
11x + 20y − 14z = 0.
Геометрический смысл: получили уравнение касательной плоскости к заданной поверхности в точке (x, y, z) = (2, −1, 0).
2014 |
38 / 74 |
Л.3. Типы математических моделей
2. Рассмотрим дифференциальное уравнение
y′′ + (1 + y2) sin y′ + y = 1,
имеющее частное решение y0(x) = 1.
Задача получить линеаризованное уравнение для функции η(x), определеящей близкое решение y(x) = 1 + η(x).
Варьируем исходное уравнение и подставляем y = y0 = 1, δy = η:
(δy)′′ + 2y(δy) sin y′ + (1 + y2)(cos y′)(δy)′ + δy = 0,
η′′ + 2η′ + η = 0.
2014 |
39 / 74 |
Л.3. Типы математических моделей
На последовательной линеаризации основан один из самых эффективных методов приближенного решения нелинейных уравнений различных типов – метод Ньютона.
2014 |
40 / 74 |