Сетевое издание

Современные проблемы науки и образования

ISSN 2070-7428

"Перечень" ВАК

ИФ РИНЦ = 1,006

• Авторы
• Резюме
• Файлы
• Ключевые слова
• Литература

Галкин А.В. 1 Дятчина Д.В. 1

---

1 Липецкий государственный технический университет

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

Статья в формате PDF

164 KB

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

инерционный трансформатор вращающего момента

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

численное интегрирование

метод Фельберга

метод Рунге-Кутты

1. Арушанян О. Б., Залеткин С. Ф. Решение задачи Коши для обыкновенных дифференциальных уравнений одношаговыми разностными методами: практикум на ЭВМ по вычислительным методам. - М.: МГУ, 2002. - 51 с.

2. Баженов С.П. Бесступенчатые передачи тяговых и транспортных машин: Учеб. пособие. - Липецк: Изд-во ЛГТУ, 2003. - 81 с.

3. Баженов С.П., Галкин А.В., Дедяев М.И. Характеристики выпрямителя момента инерционного трансформатора городского автобуса // 16-ая Всероссийская научно-техническая конференция «Проблемы и достижения автотранспортного комплекса». Материалы. - Екатеринбург, 2008. - С. 32-34.

4. Бахвалов Н. С., Жидков Н. П., Кобельков Г. М. Численные методы. - М.: Наука, 1987. - 630 с.

5. Леонов А.И. Инерционные автоматические трансформаторы вращающего момента. - М.: Машиностроение, 1978. - 224 с.

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

В качестве примера модели, заданной составными системами дифференциальных уравнений, рассмотрим инерционный трансформатор вращающего момента (ИТВМ) с учетом упругих свойств механизма свободного хода [2]. ИТВМ является бесступенчатой передачей механического типа, обладающей внутренним автоматизмом, то есть способностью автоматически изменять передаточное отношение в зависимости от угловой скорости выходного вала и величины нагрузки внешнего сопротивления [5].

Математической моделью являются системы обыкновенных дифференциальных уравнений, описывающих движение обобщенной физической модели инерционного трансформатора [3]. Работа инерционного трансформатора имеет циклический характер. В течение цикла выделяются четыре участка, на которых дифференциальные уравнения движения не меняются:

• разгон реактора до угловой скорости ведомого маховика;
• совместное движение реактора и ведомого маховика;
• торможения реактора до остановки;
• выстой реактора.

Система уравнений, описывающая первый участок:

 (1)

Условием перехода является достижение угловой скорости реактора угловой скорости ведомого звена, т.е.   В момент времени t₁ происходит переход рабочего процесса во второй такт.

Система уравнений, описывающая второй участок:

 (2)

Условием перехода является поворот сателлита в относительном движении на π радиан, т.е.  В момент времени t₂ происходит переход рабочего процесса в третий такт.

Система уравнений, описывающая третий участок, имеет такие же уравнения, что и на первом.

Условием перехода является достижение реактором угловой скорости равной нулю, т.е.  В момент времени t₃ происходит переход рабочего процесса в четвертый такт.

Система уравнений, описывающая четвертый участок:

 (3)

Условием перехода является поворот сателлита в относительном движении на 2π радиан, т.е. ,

где      - обобщенные координаты и обобщенные скорости,

a - внутреннее передаточное отношение; n - число грузовых звеньев;

m - масса грузового звена; d - расстояние от оси вращения грузового звена до его центра тяжести; e - расстояние от оси вращения ИТВМ до оси вращения грузового звена; J₂₁ - приведенный момент инерции ведущих элементов; J₂₂ - приведенный момент инерции ведущей части реактора; J_П - приведенный момент инерции ведомых элементов; J_Г- приведенный момент инерции грузового звена.

Выражения  и  являются упругими моментами и входят в уравнения в линейном виде. Коэффициенты  определяют жесткость пружины и материалов между ведомым валом и реактором.           

Начальные условия для первого участка ,   ,   , ,,

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

Решение каждой из систем дифференциальных уравнений (1)-(3) находится численно. Если использовать метод Рунге - Кутты четвертого порядка с шагом 10^-5, то при условиях (4) требуется 1634 шага с четырьмя вычислениями функции правой части на каждом шаге.

   (4)

Применим для решения данной задачи модификацию метода Рунге - Кутты, предложенную Фельбергом [1]:

  (5)

Формула имеет пятый порядок точности.

На каждом шаге требуется производить шесть вычислений функции правой части, что больше, чем при классическом методе Рунге - Кутты четвертого порядка.

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

Формула Фельберга четвертого порядка имеет следующий вид [1]:

Коэффициенты определяются по тем же формулам (5), что и для формулы Фельберга пятого порядка точности.

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

и имеет порядок точности .

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

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

• 1) если , то шаг уменьшается вдвое , где ε- заданная точность;
• 2) если , то шаг удваивается ;
• 3) если , то шаг не изменяется.

Кроме погрешности самого метода, существенное влияние на точность полученного решения имеет и погрешность начальных условий. Учитывая, что для рассматриваемой модели рабочего процесса ИТВМ в качестве начальных значений для последующих участков используются конечные значения предыдущего участка, вычисленные приближенно, то возникает погрешность начальных условий. С целью ее минимизации используем алгоритм деления шага пополам.

Обозначим C=0 условия перехода, где C определяется при численном решении систем дифференциальных уравнений. Если , то уменьшаем шаг вдвое .

Решение модели рабочего процесса методом Фельберга с ,  и начальным шагом  с представлено на рис. 1.

Рис. 1. Графики изменения угловых скоростей

Сравним результаты, получаемые методом Рунге - Кутты четвертого порядка с постоянным шагом и методом Фельберга с переменным шагом. Сравнения представлены в таблице 1.

Таблица 1. Сравнение численных методов интегрирования

Из таблицы видно, что количество вычислений при использовании метода Рунге - Кутты значительно выше, чем при использовании метода Фельберга. Так, общее количество вычислений функции правой части методом Рунге - Кутты 6536, что в 2.73 раза больше, чем методом Фельберга.

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

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

Работа поддержана грантом РФФИ № 11-07-00580-а.

Рецензенты:

• Шмырин А.М., д.т.н., профессор кафедры прикладной информатики в экономике, Липецкий эколого-гуманитарный институт, г. Липецк.
• Малыш В.Н., д.т.н., профессор, заведующий кафедрой электроники, телекоммуникаций и компьютерных технологий, Липецкий государственный педагогический университет, г. Липецк.

Библиографическая ссылка