CC BY
86
УДК: 517.956.32
Ф. В. Голованёва, С. А. Шабров, Меач Мон
АДАПТАЦИЯ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ
ДЛЯ ОДНОЙ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ВТОРОГО ПОРЯДКА
С НЕГЛАДКИМИ РЕШЕНИЯМИ
В работе метод конечных элементов распространяется на математическую модель второго порядка с негладкими решениями.
Ключевые слова: математическая модель, метод конечных элементов, негладкие решения.
Метод конечных элементов мы применим для нахождения приближённого решения математической модели, задаваемой следующей системой
м'Лх>"„= „~uQст=
rlu(0,t)-u'x(0,t) = 0, (1)
• y2u(l,t) + ux(l,t) = 0, и(х, 0) = <р0(х\ ut(x, 0) = ЙО),
которая возникает при описании малых поперечных собственных колебаний растянутой с силой p(x) вдоль отрезка [0;/] струны, оба конца
которой имеют упругое закрепление с помощью пружин жесткости у1 и у2 соответственно. Кроме того, вся система помещена во внешнюю среду с локальным коэффициентом упругости dQ, а на струне произвольным образом распределены массы. Отметим, что скачки функции Q(x) соответствуют локализованным особенностям (типа пружины). Функции р0(х) и <ру(х) — это начальные отклонения от положения равновесия и скорость объекта соответственно.
Решение задачи (1) мы ищем в классе E — непрерывных по сово-
Голованёва Фаина Валентиновна — кандидат физико-математических наук, доцент (Воронежский государственный университет, Воронеж); e-mail: shaspoteha@mail.ru.
Шабров Сергей Александрович — кандидат физико-математических наук, доцент (Воронежский государственный университет, Воронеж); e-mail: shabrov_s_a@math.usu.ru.
Меач Мон — кандидат физико-математических наук, доцент (Университет Khmerak, Пномпень, Королевство Камбоджа); e-mail: meach_mon@yahoo.com.
© Голованёва Ф. В., Шабров С. А., Меач Мон, 2016
89
купносги переменных на [0;/]х[0;Г] функций, имеющих непрерывные частные производные и и иа при каждом фиксированном хе[0;/]. При этом: и(х,?) — абсолютно непрерывна на [0;/] при фиксированном /; ри:(х,1) — <т — абсолютно непрерывна на [0;/] для каждого постоянного / е [0; 7 ]; смешанные частные производные и 1х и Их1 равны почти всюду. Функция <т(х), порождающая на [0;/] меру, содержит все особенности системы, а именно: у ег(х) есть скачки в тех и только тех точках, в которых имеется скачок хотя бы у одной из функций р(х), Q(х) или М(х), входящих в дифференциальное уравнение системы (1).
Для нахождения приближённого решения изучаемой модели мы адаптируем к ней метод конечных элементов. Решение (х, ^) ищем в
N
виде ^ак(1><Рк(х>г где <рк(х) — классические базисные функции, ак(1) —
к=0
неизвестные функции. Построение <рк(х) произведём следующим образом: отрезок [0;/] разобьём на N равных частей упорядоченным набором точек хкг к = 0,1,...,7^ (х0=0, хы=1 и добавим точки х_1=-1/Ы и =/ + |/л ); (рк(х) равна единице в точке хк, нулю на |О;хк , | и[хм;/] и продолжена линейной функцией на | хк ,: хк | и | хк: х, , |.
Применим стандартную схему метода конечных элементов: умножим уравнение в (1) на базисную функцию, проинтегрируем по мере сг (именно в этом и состоит отличие от классического метода). В результате для нахождения ак (^) получим линейную систему линейных дифференциальных уравнений второго порядка, для численного решения которой можно применять различные разностные схемы. Можно показать, что получающаяся система, дополненная начальными данными, которые определяются из начальных условий модели (1), имеет единственное решение.
Доказано, что скорость сходимости приближённого решения к точному имеет порядок Ы~112.
Метод конечных элементов удалось применить к изучаемой модели ввиду того, что дифференциальное уравнение в (1) является поточечно-заданным. Последнее стало возможно после применения производных по мере в рамках концепции Ю. В. Покорного [5], которая получила развитие в работах [1—4], [6—8] его учеников для различных математических моделей с особенностями.
Список литературы
1. Голованёва Ф. В. О функции Грина некоторых негладких задач: дисс. ... канд. физ.-мат. наук; Воронежский государственный университет. Воронеж, 2007. 101 с.
2. Зверева М. Б. О некоторых вопросах качественной теории дифференциаль-
90
ных уравнений с производными Стилтьеса: дисс. ... канд. физ.-мат. наук; Воронежский государственный университет. Воронеж, 2005. 120 с.
3. Зверева М. Б., Шабров С. А., Лилов Е. В. Об адаптации метода конечных элементов для решения граничной задачи с дифференциалами Стилтьеса на геометрическом графе / / Вестник Воронежского государственного университета. Серия: Физика. Математика. 2014. № 1. С. 97 — 105.
4. Покорный Ю. В. Дифференциал Стилтьеса в импульсных задачах с разрывными решениями / Ю. В. Покорный, М. Б. Зверева, С. А. Шабров, М. Б. Давыдова / / Доклады Академии Наук. 2009. Т. 428. № 5. С. 595 — 597.
5. Покорный Ю. В. Интеграл Стилтьеса и производные по мере в обыкновенных дифференциальных уравнениях / / Доклады Академии НАУК. 1999. Т. 364. № 2. С. 167—169.
6. Покорный Ю. В., Зверева М. Б., Шабров С. А. Осцилляционная теория Штурма-Лиувилля для импульсных задач / / Успехи математических наук. 2008. Т. 63. № 1. С. 111 — 154.
7. Шабров С. А. Об одной математической модели малых деформаций стержневой системы с внутренними особенностями // Вестник Воронежского государственного университета. Серия: Физика. Математика. 2013. № 1. С. 232 — 250.
8. Pokornyi Yu. V. An Irregular Extension of the Oscillation Theory of the Sturm-Liouville Spectral Problem / Yu. V. Pokornyi, M. B. Zvereva, S. A. Shabrov, A. S. Ishchenko // Mathematical Notes. 2007. Т. 82. № 3—4. С. 518 — 521.
•Jc -Jc -Jc
Golovaneva Faina V.1, Shabrov Sergey A.2, Meach Mon3
DAPTATION FINITE ELEMENT METHODS FOR ONE MATHEMATICAL MODEL
OF THE SECOND ORDER WITH NONSMOOTH SOLUTIONS
i1, 2 Voronezh State University, Voronezh, Russia;
University Khmerak, Phnom Penh, Kingdom of Cambodia)
In the work the finite element method covers the mathematical model of the second order with
nonsmooth solutions.
Keywords: mathematical model, the finite element method, nonsmooth solutions.
References
1. Golovaneva F. V. O funktsii Grina nekotorykh negladkikh zadach (On the green function of some nonsmooth problems): the dissertation on competition of a scientific degree of candidate of physical and mathematical Sciences, Voronezh, 2007. 101 p.
2. Zvereva, M. B. O nekotorykh voprosakh kachestvennoy teorii differentsial'-nykh uravneniy s proizvodnymi Stilt'esa (On some problems of qualitative theory of differential equations with derivatives of the Stieltjes): the dissertation on competition of a scientific degree of candidate of physical and mathematical Sciences, Voronezh, 2005. 120 p.
3. Zvereva M. B., Shabrov S. A., Lilov E. V. Adaptation of the Finite Elements Method for Solution of a Boundary Value Problem with Stieltjes Differentials on a Graph [Ob adaptatsii metoda konechnykh elemen-tov dlya resheniya granichnoy zadachi s differentsialami Stilt'esa na geo-metricheskom grafe], Vestnik Voronezhskogo gosudarstvennogo universiteta, Series: Fizika. Matematika, 2014, no. 1, pp. 97—105.
4. Pokornyi Yu. V., Zvereva M. B., Shabrov S. A., Davydova M. B. Stieltjes differential pulsed problems with discontinuous solutions [Differentsial Stilt'esa v
91
impul'snykh zadachakh s razryv-nymi resheniyami], Doklady Akademii Nauk, 2009, vol. 428, no. 5, pp. 595-597.
5. Pokornyi Yu. V. The Stieltjes Integral and Derivatives with Respect to the Measure in Ordinary Differential Equations [Integral Stilt'esa i proizvodnye po mere v obyknovennykh differentsial'nykh uravneniyakh], Doklady Akademii Nauk, 1999, vol. 364, no. 2, pp. 167-169.
6. Pokorny Yu. V., Zvereva M. B., Shabrov S. A. Oscillation theory of the Sturm-Liouville problem for impulsive problems [Ostsillyatsionnaya teoriya Shturma-Liuvillya dlya impul'snykh zadach], Uspekhi matematicheskikh nauk, 2008, vol. 63, no. 1, pp. 111-154.
7. Shabrov S. A. Mathematical model for small deformations of the rod system with internal features [Ob odnoy matematicheskoy modeli malykh deformatsiy sterzhne-voy sistemy s vnutrennimi osobennostyami], Vestnik Voronezhskogo gosudarstvennogo universiteta, Series: Fizika. Matematika, 2013, no. 1, pp. 232 — 250.
8. Pokornyi Yu. V., Zvereva M. B., Shabrov S. A., Ishchenko A. S. An Irregular Extension of the Oscillation Theory of the Sturm-Liouville Spectral Problem, Mathematical Notes, 2007, vol. 82, no. 3—4, pp. 518 — 521.
•Jc -Jc -Jc
92
