Алгоритм решения прикладных задач для обыкновенных дифференциальных уравнений четвертого порядка с методом дифференциальной прогонки | Статья в журнале «Молодой ученый»

Отправьте статью сегодня! Журнал выйдет 28 декабря, печатный экземпляр отправим 1 января.

Опубликовать статью в журнале

Авторы: , ,

Рубрика: Технические науки

Опубликовано в Молодой учёный №6 (86) март-2 2015 г.

Дата публикации: 19.03.2015

Статья просмотрена: 2205 раз

Библиографическое описание:

Олимов, Муродилла. Алгоритм решения прикладных задач для обыкновенных дифференциальных уравнений четвертого порядка с методом дифференциальной прогонки / Муродилла Олимов, О. О. Жакбаров, Ф. С. Ирискулов. — Текст : непосредственный // Молодой ученый. — 2015. — № 6 (86). — С. 193-196. — URL: https://moluch.ru/archive/86/16292/ (дата обращения: 19.12.2024).

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

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

 

Приближенные решения уравнений математической физики, дифференциальных и интегральных чаще всего строятся на основе так называемых прямых методов, позволяющих свести построение приближенного решения данной задачи к решению систем линейных (если данная задача также линейная) алгебраических уравнений. Как хорошо известно, прямые методы прекрасно зарекомендовали себя на практике и получили широкое распространение; чаще других используются различные варианты метода сеток. Однако наряду со многими достоинствами, прямые методы имеют и один важный недостаток: чтобы получить более высокую точность приближения, приходится повышать порядок алгебраической системы, что часто приводит к неустойчивости вычислительного процесса. Для метода сеток увеличение порядка алгебраической системы ведёт к увеличению числа ее обусловленностей, что влечет за собой накопление ошибок при решении системы и, как следствие понижение точности. Добавим к этому, что машинное решение алгебраических систем высокого порядка наталкивается еще на специфические трудности, связанные с ограниченностью памяти ЭВМ.

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

Будем рассматривать краевые задачи только для линейных обыкновенных дифференциальных уравнений четвертого порядка с переменными коэффициентами. Распространение предлагаемого метода в сочетании с методом Канторовича — Власова или методом прямых на уравнения с частными производными в ряде случаев очевидно.

Рассмотрим дифференциальное уравнение

                                             (1)

где -известные функции при различных граничных условиях.

I.          Пусть граничных условиях для уравнения (1) будут следующие:

                                                                                   (2)

Введя обозначения

                                                                                                     (3)

Напишем уравнение (1) в виде

                                                                      (4)

Где

                                                             (5)

Ищем решение уравнение (4)

                                                                                             (6)

С начальным условием

                                                                                                                       (7)

Где - 2-мерная матрица, - 2-мерный вектор. Тогда для нахождения  и получим уравнения

                                                                       (8)

С начальными условиями

                                                                                                            (9)

Таким образом, краевая задача (1)-(2) свелась к задачам Коши (6)-(7), (8)-(9).

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

1. Изгиб балки переменного сечения, лежащей на упругом основании. Основное разрешаемое уравнение-

                                                                   (10)

где Е — модуль упругости, - момент инерции сечения, - переменный коэффициент.

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

Следовательно, полагая ,  и , из решений уравнения (1) найдем решения уравнения (3).

Если один край балки жестко защемлен, а другой — свободный, то граничные условия для уравнения (10) будут такие:

                               (11)

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

Основное разрешаемое уравнение имеет вид [3]

                                                                   (12)

где - переменная толщина, - радиус оболочки, - коэффициент Пуассона.

Уравнение (12) полностью идентично уравнению (16), поэтому решения о соответствующими граничными условиями получаем аналогичным образом.

3. Изгиб прямоугольной пластинки переменной жесткости.

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

                                       (13)

Граничные условия для пластинок следующие:

a)        жестко защемленной

 при ;                                                                                   (14)

b)        свободно опертой

 при ;                                            (15)

c)         свободной

 при ;                                                                 (16)

Решение уравнения (13) ищем в виде

                                                                                                 (17)

где - известные координатные функции, удовлетворяющие граничным условиям пластинки при y=0, y=b. Подставляя предполагаемое решение (17) и (13) и выполняя процедуры метода Канторовича — Власова, придем к уравнению

,                                          (18)

где

,

Таким образом, уравнения (10), (12) и (18) соответствующими граничными условиями с успехом могут быть решены предлагаемым методом.

 

Литература:

1.         С. С. Ирискулов, М. Олимов, К. Исманова, А. Имомов. Численные методы и алгоритмы. Учебное пособие., Наманган, Изд-во «Наманган», 2013, 278с.

2.         Олимов М., Каримов П., Исмоилов Ш.. К решению краевых задач пространственных стержней при переменных упруго-пластических нагружений с учетом разгрузки // Научно- технический журнал ФерПИ, Фергана — 2014. — № 3 — с. 113–116.

3.         Тимошенко С. П., Войновский-Кригер С. Пластинки и оболочки. М., «Наука», 1966.

Основные термины (генерируются автоматически): уравнение, задача, переменная толщина, предлагаемый метод, Решение уравнения, алгебраическая система, основное разрешаемое уравнение, переменная жесткость, переменное сечение, прямоугольная пластинка.


Похожие статьи

Решение уравнения колебаний балки при шарнирном закреплении на границах

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

Об одном методе решения линейных интегральных уравнений

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

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

Рассматривается задача разработки и использования методов параллельного программирования при численном решении задач оптимизации методами координатного и градиентного спусков. Задача оптимизации рассматривается в контексте решения задачи гашения коле...

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

Рассматривается задача исследования устойчивости разностных схем для численного решения уравнений колебаний прямоугольной мембраны и прямоугольной пластины. Исследование проводится методом Неймана. Выводятся соотношения зависимости шага по времени от...

Математическое моделирование динамики вязкоупругих трубопроводов с протекающей жидкостью

На примере вязкоупругой оболочки рассмотрены задачи о колебаниях вязкоупругих трубопроводов с протекающей жидкостью. С помощью метода Бубнова — Галеркина математическая модель задачи сведена к исследованию системы обыкновенных интегро-дифференциальны...

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

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

Многомерная интерполяция сеточной вектор-функции

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

Применение метода вариационных итераций к приближенному решению нелинейных обыкновенных дифференциальных уравнений

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

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

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

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

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

Похожие статьи

Решение уравнения колебаний балки при шарнирном закреплении на границах

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

Об одном методе решения линейных интегральных уравнений

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

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

Рассматривается задача разработки и использования методов параллельного программирования при численном решении задач оптимизации методами координатного и градиентного спусков. Задача оптимизации рассматривается в контексте решения задачи гашения коле...

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

Рассматривается задача исследования устойчивости разностных схем для численного решения уравнений колебаний прямоугольной мембраны и прямоугольной пластины. Исследование проводится методом Неймана. Выводятся соотношения зависимости шага по времени от...

Математическое моделирование динамики вязкоупругих трубопроводов с протекающей жидкостью

На примере вязкоупругой оболочки рассмотрены задачи о колебаниях вязкоупругих трубопроводов с протекающей жидкостью. С помощью метода Бубнова — Галеркина математическая модель задачи сведена к исследованию системы обыкновенных интегро-дифференциальны...

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

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

Многомерная интерполяция сеточной вектор-функции

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

Применение метода вариационных итераций к приближенному решению нелинейных обыкновенных дифференциальных уравнений

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

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

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

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

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

Задать вопрос