Решение начально-краевых задач с подвижной границей Текст научной статьи по специальности «Математика»

Фундаментальные проблемы теоретической и прикладной механики Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (5), с. 2423-2424

2423

УДК 539.3,534.1

РЕШЕНИЕ НАЧАЛЬНО-КРАЕВЫХ ЗАДАЧ С ПОДВИЖНОЙ ГРАНИЦЕЙ © 2011 г. Н.А. Пинчук, А.М. Столяр

Южный федеральный университет, Ростов-на-Дону

ajoiner@mail.ru

Поступила в редакцию 24.08.2011

Исследуются начально-краевые задачи для уравнений гиперболического типа - линейных и нелинейных - с подвижной границей. Они моделируют продольные и продольно-поперечные колебания упругого троса переменной длины. Рассмотренные постановки восходят к работам А.Ю. Ишлинского [1], Г.Н. Савина и О.А. Горошко [2], А.И. Весницкого [3] и др. Предлагаются асимптотический и численные методы решения. При асимптотическом анализе в качестве малого параметра используется отношение скорости изменения длины троса к скорости распространения волны. Численные методы являются модификациями методов ко -нечных разностей и Рунге-Кутта. Проводится сравнение результатов асимптотического и численного интегрирования.

Ключевые слова: начально-краевые задачи, подвижные границы, асимптотическое и численное интегрирование.

1. Продольные колебания троса переменной длины. Построение асимптотики

Моделируемая система состоит из троса, ко -торый может разматываться с катушки или наматываться на нее с заданной скоростью, и висящего на ней груза [1]. Поведение данной системы можно описать следующими уравнениями [4]:

д 2и ~д^Г

2 д2и

= а —- + £,

Е

дх

2

= _, и| () = 0,

р 1х=_(

и(х,ґ) = ^(ґ) - х - и(х,ґ),

д 2и

х=0

ЕЕ

т

ди

дх

+1

х=0 у

&__

&ґ

1+

ди(_, ґ)

д_

ди

Ч=о = фі(хх дґ _ «>1,=0 = _0.

= 8±у(ґ),

= Ф 2 (X),

ґ=0

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

д 2у

+ Р(е±, ґ) = а

+ Р(в±, ґ) =

2 _д^ дх2 ’ ЕЕ ду

х=0

т дх

х=0

Іґ= 0

= Фц( х),

(2)

ду

д"

= Ф 22 (х), _ (ґ) = _0 +8±_і(ґ),

ґ = 0

(1)

_і (ґ)=1 ^(т)&т, у

х=_ (ґ)

= 0.

Здесь и(х, ґ) - продольное смещение троса в подвижной системе координат, связанной с грузом; ^(ґ) - расстояние от точки подвеса троса до груза, _(ґ) - длина троса в недеформированном состоянии; є±у (ґ) — скорость изменения длины троса, е± = е у 0, когда длина троса увеличивается, е± = -е Р 0, когда длина троса уменьшается, е = V* / а; V* — максимальная скорость изменения длины троса за время рассматриваемого процесса; предполагается, что е РР 1; остальные обозначения понятны.

Здесь функции Р(8±,X),Ф11(х),Ф22(х) выражаются через исходные параметры. Задача (2) содержит малый параметр £ как в уравнении, так и в краевых условиях, причем на краю х = /(X) он входит в выражение для самой границы. Применяя метод асимптотического интегрирования, решение задачи (2) строим в виде

(3)

V = Х V, (х, X)8 *

к=0

и для функций Vк (х, X) получаем начально-краевые задачи уже на постоянном отрезке:

д'2vk . 2 д2v

—2- + а к Рк(х) = а ----------

2

д Ч

+ а к рк(ґ) =

дх ЕЕ дук

х=0

т дх

х=0

2

а

k-1

vkL=/0 =-Z

[ll(t)]

k-i

д k-ivk

~0 (k -i)! dx

k-i

Х=І0

vk\t=0 =pk Ф11( x)

dvk

dt

= Yk Ф 22(x). (4)

t=0

2. Продольно-поперечные колебания троса переменной длины. Построение асимптотики

Для описания продольно-поперечных колебаний в плоской постановке добавим к уравнениям задачи (1) уравнения:

д 2 w

dt

w

2

2 д (du dw ^ . .

=a q(xt)’

dx І dx dx I

(5)

x= 0

= fl(t),

w

w

= *.( x), £

l (t) = f2(t), = %( x).

t=0

Іґ=° дґ

Здесь через д(х, ґ) обозначена приведенная внешняя нагрузка, действующая на трос в поперечном направлении; функции/^),/2(ґ) и Т1 (х), Т2 (х) задают краевые и начальные условия соответственно. В выражение одной из границ задачи (5) входит малый параметр. Разлагая функцию ^(х, ґ) в ряд Тейлора в окрестности х = /0, строим ее попутно с разложением (3) функции у(х, ґ) в виде

w

= Z wk(x 0є k •

(б)

k=0

Для функций (х, ґ) аналогично п. 1 получаем

начально-краевые задачи на постоянном отрезке.

3. Численное интегрирование

Проведено численное решение исходных задач (1) и (1), (5) на переменном промежутке [0, /(0] и решение с использованием разложений (3) и (3), (6) соответственно. Для решения исходных задач разработаны модификации метода ко -

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

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

Работа выполнена при частичной поддержке РФФИ, грант № 06-01-00287.

Список литературы

1. Ишлинский А.Ю. // Укр. мат. журн. 1953. Т. 5, №4. С. 370—374.

2. Савин Г.Н., Горошко О. А. Динамика нити переменной длины. Киев, 1962.

3. Весницкий А.И. Волны в системах с движущимися границами и нагрузками. М.: Физматлит, 2001. С. 320.

4. Кечеджиян Л.О., Пинчук Н.А., Столяр А.М. // Изв. вузов. Сев.-Кав. рег. Естеств. науки. 2008. № 1. С. 22—27.

x

SOLUTION OF INITIAL BOUNDARY-VALUE PROBLEMS WITH MOVING BOUNDARIES

N.A. Pinchuk, A.M. Stolyar

Initial boundary-value problems with moving boundaries for linear and nonlinear hyperbolic type equations are investigated. They model longitudinal and transverse-longitudinal oscillations of elastic cable of a variable length. Formulations of the problem considered go back to the papers by A.J. Ishlinsky [1], G.N. Savin and O.A. Goroshko [2], A.I. Vesnitsky [3]. Asymptotic and numerical methods are developed. Asymptotic integration is carried out when the small parameter is equal to the relation of velocities of the cable length variation and wave propagation. Numerical integration is done using the modified finite-difference and Runge-Kutta methods. Comparison of asymptotic and numerical solutions is carried out.

Keywords: initial boundary-value problems, moving boundaries, asymptotic and numerical integration.