Реальная система, физическое тело, или процесс, динамические свойства которой(ого) изучаются с помощью моделирования

Страницы работы

Фрагмент текста работы

Основные понятия и определения

Объект

Реальная система, физическое тело, или процесс, динамические свойства которой(ого) изучаются с помощью моделирования.

Физическая модель

Физический макет, динамические свойства которого соответствуют динамическим свойствам изучаемого объекта.

Математическая модель

Система дифференциальных уравнений или структурная схема, адекватно описывающая объект.

Моделирование

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

Структура программы моделирования

Существует два подхода к моделированию:

1.  Моделирование по системе дифференциальных уравнений.

2.  Моделирование по структурной схеме.

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

1.  Основная часть программы выполняет следующие действия:

a.  вызов подпрограммы ввода исходных данных и расчета коэффициентов дифференциальных уравнений;

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

2.  Подпрограмма ввода исходных данных и расчета коэффициентов дифференциальных уравнений.

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

4.  Подпрограмма определения значений правых частей уравнений. Данная подпрограмма составляется на основе конкретной системы дифференциальных уравнений в форме Коши. Для определения значений входных воздействий модели вызывается соответствующая подпрограмма.

5.  Подпрограмма определения значений входных воздействий модели.

При моделировании по структурной схеме программа должна содержать:

1.  Основную часть программы, выполняющую следующие действия:

a.  вызов подпрограммы ввода исходных данных и расчета параметров звеньев;

b.  организация цикла по времени, в котором вызывается подпрограмма реализации модели и выводятся результаты.

2.  Подпрограмма ввода исходных данных и расчета параметров звеньев.

3.  Подпрограмма реализации модели. Данная подпрограмма составляется на основе конкретной модели. В ней осуществляется вызов подпрограммы определения значений входных воздействий модели и подпрограмм для элементарных звеньев в порядке следования сигналов в модели.

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

5.  Подпрограмма определения значений входных воздействий модели.

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

АЭПНачальная

<Предыдущая

>Следующая

На уровень вышеВыше

Посещений - 570. , 2002.

Численные методы решения систем дифференциальных уравнений

В табл. 1 приведены расчетные соотношения для нескольких методов численного решения дифференциальных уравнений и их систем.

Таблица 1

Расчетные соотношения для численных методов решения систем дифференциальных уравнений

Название метода

Расчетное соотношение

Эйлера

Yn+1 = Yn + Δt · ƒn

Усовершенствованный Эйлера

Yn+1(1) = Yn + Δt · ƒn; Yn+1(2) = Yn + 0.5 · Δt · (ƒn + ƒn+1(1))

Неявный Эйлера

Yn+1(ν+1) = Yn + Δt · ƒn(ν)

Рунге-Кутта

K1 = ƒ(Yn , Xn); K2 = ƒ(Yn + 0.5 · Δt · K1 , Xn+0.5); K3 = ƒ(Yn + 0.5 · Δt · K2 , Xn+0.5); K4 = ƒ(Yn + Δt · K3 , Xn+1);

Yn+1 = Yn

Δt

 · (K1 + 2 · K2 + 2 · K3 + K4)

6

Адамса-Бэшфорта

Yn+1 = Yn + 0.5 · Δt · (3 · ƒn - ƒn-1)

Адамса-Мултона

Yn+1(1) = Yn

Δt

 · (55 · ƒn - 59 · ƒn-1 + 37 · ƒn-2 - 9 · ƒn-3)

;

24

Yn+1(2) = Yn

Δt

 · (9 · ƒn+1(1) + 19 · ƒn - 5 · ƒn-1 + ƒn-2)

24

Хэмминга

Yn+1(1) = Yn-3

4

 · Δt · (2 · ƒn - ƒn-1 + 2 · ƒn-2)

;

3

Yn+1(2) = Yn+1(1)

112

 · (Yn - Yn(1))

;

121

Yn+1(ν+1)

1

 · (9 · Yn - Yn-2) + 

3

 · Δt · (ƒn+1(ν) + 2 · ƒn - ƒn-1)

8

8

Башарина

Yn+1(ν+1) = Yn + Δt · ƒ(

Yn+1(ν) + Yn

 , Xn+0.5)

2

В расчетных соотношениях обозначено:

X - 

входные воздействия модели;

Y - 

внутренние переменные модели;

ƒ - 

значения правых частей дифференциальных уравнений;

Δt - 

шаг расчета;

верхний индекс - 

номер итерации (ν + 1 - текущая итерация, ν - предыдущая и т.д.);

нижний индекс - 

номер шага расчета (n + 1 - текущий шаг расчета, n - предыдущий

Похожие материалы

Информация о работе