Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2012. № 3 (28). С. 145—151
УДК 517.956.3:534.11
ОБ ОДНОМ МЕТОДЕ ПОЛУЧЕНИЯ АНАЛИТИЧЕСКОГО РЕШЕНИЯ ВОЛНОВОГО УРАВНЕНИЯ, ОПИСЫВАЮЩЕГО КОЛЕБАНИЯ СИСТЕМ С ДВИЖУЩИМИСЯ ГРАНИЦАМИ
В. Н. Анисимов, В. Л. Литвинов, И. В. Корпен
Сызранский филиал Самарского государственного технического университета,
446001, Россия, Самарская обл., Сызрань, ул. Советская, 45.
E-mails: [email protected], [email protected]
Описан аналитический метод решения волнового уравнения с условиями, заданными на движущихся границах. С помощью замены переменных в системе функциональных уравнений исходная краевая задача сведена к системе разностных уравнений с одним постоянным смещением, которая может быть решена с помощью интегрального преобразования Лапласа. Получено выражение для амплитуды колебаний, соответствующих n-ной динамической моде в случае граничных условий первого рода. Данный метод позволяет рассмотреть более широкий класс граничных условий по сравнению с другими аналитическими методами решения краевых задач с движущимися границами.
Ключевые слова: волновое уравнение, колебания систем с движущимися границами, законы движения границ, амплитуда колебаний.
Одномерные системы с движущимися границами широко распространены в технике (канаты грузоподъемных установок [1], гибкие звенья передач [2] и т. д.). Наличие движущихся границ вызывает значительные затруднения при описании таких систем, поэтому здесь в основном используются приближенные методы решения [3,4]. Из аналитических методов наиболее эффективным является метод, предложенный в [5], который заключается в подборе новых переменных, останавливающих границы и сохраняющих инвариантность уравнения. В [6] решение ищется в виде суперпозиции двух волн, бегущих навстречу друг другу. В результате этого автору удалось решить волновое уравнение с граничными условиями первого рода, заданными на одной движущейся и одной неподвижной границах. Эффективен также метод, используемый в [7] и заключающийся в замене геометрической переменной на чисто мнимую, что позволяет свести волновое уравнение к уравнению Лапласа и применить для решения методику теории функций комплексного переменного.
Решения, полученные с помощью перечисленных методов, ограничены граничными условиями первого рода. К недостаткам методов относится также то, что в случае двух движущихся границ начальные условия, заданные при t = 0, не могут быть учтены.
Перечисленных недостатков лишен развиваемый в данной статье метод решения таких задач, в котором сочетается методика, используемая в [5,6]. Пусть движение системы описывается волновым уравнением
UTT (£,т) - % (£,Т )=° (1)
Валерий Николаевич Анисимов (к.ф.-м.н., доц.), доцент, каф. общетеоретических дисциплин. Владислав Львович Литвинов, преподаватель, каф. общетеоретических дисциплин. Инна Владимировна Корпен (к.п.н.), доцент, каф. общетеоретических дисциплин.
145
Анисимов В.Н., Литвинов В.Л., Корпен И.В.
при граничных условиях первого рода
U (li (т ),т ) = Fi(t ), li (0) = 0;
U(I2(т),т) = F2(т), 12(0) = 1, I2(т) > Ii(t)
(2)
и начальных условиях
и(£,о) = Фо(о, ur(e,o) = $i(e). (3)
Здесь т ^ 0 — безразмерное время; £ € [0,1] —безразмерная пространственная координата (li(т) ^ ^ 12(т)); 11(т), l2(т) —законы движения границ;
Фо(£), Ф1(£), F1 (т), F2(t) —заданные функции, допускающие разрывы первого рода.
Для решения задачи используем представление Даламбера. Общее решение уравнения (1) имеет вид
u(£,т)= 9(т + £) + G(t - £) (4)
где g(z) и G(z) — произвольные функции, которые необходимо определить из начальных и граничных условий, z — независимая переменная.
Подставляя решение (4) в граничные (2) и начальные (3) условия, получим
д(т + Ii(t)) + G(t - Ii(t)) = Fi(t), д(т + I2(т)) + G(t - I2(т)) = F2(т);
д(£) + G(-£) = Фо(£), 0 < £ < 1,
д'(£) + G'(-£) = Ф1(£), 0 < £ < 1.
Из системы (6) найдём функции д(£) и G(£):
д(£)
G(£)
1
2
1
2
Фо(£)+ /V(CR ,
J о J
Фо(-£) + / Ф1(-СК ,
о
0 < £ < 1; 1 < £ < 0.
(5)
(6)
(7)
В отличие от метода А. И. Весницкого [5], где в дифференциальное уравнение вводятся новые переменные, останавливающие границы и сохраняющие инвариантность уравнения, для упрощения задачи введём в систему (5) новые функции
д(^ = R(z)); G(z) = R(^(z))- (8)
Тогда система уравнений (5) примет вид I
I г(^(т + I1 (т))) + Я(^(т - I1 (т))) = F1(t); |R(t + I2(т))) + R(^(t - I2(т))) = F2(т).
Введём обозначения в первом уравнении системы (9):
Р(т + I1 (т)) = z; ^(т - I1 (т)) = z,
(9)
146
Об одном методе получения аналитического решения волнового уравнения ...
а во втором уравнении этой системы —
р(т + Ь(т)) = z; ф(т - l2(т)) = z - 1.
Другими словами, если функции <^(z) и ^(z) удовлетворяют системе
<^(т + li(т)) = ф(т - 1\(т));
^(т + I2(т)) = ф(т - I2(т)) + 1,
то система (5) принимает вид
r(z) + R(z) = di (z); r(z) + R(z - 1) = 02(z),
(10)
(11)
где di(z) = Fi(0,5<£(z) + 0,5'0(z)); 6>2(z) = F2(0,5<£(z)+ 0,5^(z - 1)). Здесь <£(z), •0(z) — функции, обратные к <^(z) и ^(z).
Заметим, что из системы (10) функции <^(z) и ^(z) определяются с точностью до константы в том смысле, что если <^(z) и ^(z) — решение системы (10), то <^(z) + C и ^(z) + C также являются решением (здесь C — произвольная постоянная). Поэтому для определённости можно выбрать такую функцию ^(z), что •0(-1) = -1. При этом из второго уравнения системы (10) при т = 0 следует, что <^(1) = 0. Из первого уравнения системы (10) при т = 0 получим ^>(0) = ^(0).
С учётом замены (8) начальные условия (6) примут вид
r(z) = g(<p(z)), р(0) < z < 0;
R(z) = G(ip(z)), -1 < z < ^(0),
где функции g(z) и G(z) определяются выражениями (7).
Таким образом, начальная задача (1)—(3) сведена к системе разностных уравнений (11) с одним постоянным смещением при начальных условиях (12). Из первого уравнения системы (11) получим
R(z) = 0i(z) - r(z). (13)
После подстановки (13) во второе уравнение (11) будем иметь
r(z) - r(z - 1) = 0(z), (14)
где 0(z) = 62(z) - 0i(z - 1). Для решения задачи (14), (12) применим интегральное преобразование Лапласа:
r(p)
6(р) _р f° r(z)e pzdz
1 -е~Р+е J_x (1 - е-Р) ’
где (9(p) —изображение функции 0(z). Оригинал данного изображения имеет вид
r(z)
/V)
J 0
го го
^2 e2nni(z-c)d( +
n=-ro П= — ГО
e2 nniz
Г
r(Z )e-2nniZ dZ.
147
Анисимов В.Н., Литвинов В.Л., Корпен И.В.
Объединяя члены при положительных и отрицательных n, получим
rz ГО rz
r(z) = 0(СЖ + 2V / )cos[2nn(z — ()]dZ+
Jo n=1Jо
/о ro r 0
r(()d( + 2V / r(()cos[2nn(z - Z)\d(. (15)
1 n=lJ-1
Рассмотрим свободные колебания системы (9(z) = 0). В этом случае из
(4) с учётом (8), (13) и (15) следует, что
и (z,t ) = £ кж,т),
(16)
n=1
где
КТ(К т) = sin{nnK<T + О — ф(т — Щх
х (АП cos{nnf^(T + О + ф(т — Щ — В„ sin{nn[^(T + О + Ф(т — Z)]}); (17)
/о г о
^ r(Z) sin(2nnZ)dZ; ВП = 4 у ^ r(Z) cos(2nn()d(.
Функцию КГ(С,т) будем называть n-ным собственным колебанием системы, а функцию w(Z,t) = (д/дт){пп[^(т + Z) + ф(т — Z)]} — мгновенной собственной частотой n-ного собственного колебания. Функция sin {пп[^(т+Z) — ф(т—Z)]} характеризует форму колебаний, и её принято называть n-ной динамической модой системы.
Рассмотрим теперь вынужденные колебания системы. При нулевых начальных условиях из (4) с учётом (8), (13) и (15) получим
ГО
U(фт) = £ Vn(Z,T) + D(Z,t), (18)
n=1
где
V.n(Z,T) = 4sin{nn[^(T + Z) — ф(т — Z)]}x
Г'Ф(т-£)
x (cos{nn[^(T + Z) + ф(т — Z)]} 0(Z)sin(2nnZ)d(—
о
f Ф(т-С)
— 4sin{nn[^(T + Z) + ф(т — Z)}} e(Z)cos(2nnZ)dZ; (19)
о
D(Z,t )
гВут+0 ф(т-С)
9(Z )dZ + 9фф(т
Z)) +
ГО
+ 2£
n= 1
г<р(т+0 Ф(т-С)
d(Z) cos{2nn[^(T + Z) — Z ]}dZ-
(20)
148
Об одном методе получения аналитического решения волнового уравнения ...
Выражение (20) представляет собой ряд Фурье для функции
вх{ф{т - £)), С = к(т);
0(<р(т + С)) + ^0(т - С)), 1i(t) < С < 12(т).
При решении задач на резонансные свойства рассматриваются, главным образом, резонансные явления в механических объектах с движущимися границами, когда амплитуда колебаний во много раз превосходит амплитуду возмущающего воздействия. Поэтому в равенстве (18) функцией D(C,t) можно пренебречь как функцией одного порядка малости с функциями F1(t) и F(т), характеризующими возмущающие воздействия.
Сравнивая выражения (16) и (18), (17) и (19), нетрудно заметить, что вынужденные колебания представляют собой суперпозицию собственных колебаний с изменяющимися во времени амплитудами:
гФ(т-£) гФ(т-О
Ап = 4 / в(() sin(2nnZ)d(; Bn = 4 / в(()cos(2nn()d(. (21)
J 0 J 0
Заметим, что из всех видов внешних воздействий наиболее распространенными являются гармонические нагрузки. Ограничимся рассмотрением случая, когда
d(z) = B cos W (z),
где функция W (z) класса C1, а B — константа, характеризующая интенсивность нагрузки. В этом случае равенства (21) можно переписать следующим образом:
Г Ф(т-?) Г Ф(т-?)
Ап = 2B / sin Фп1(С )d( + / sin Фп2(С )dZ
0 0
Г Ф(т-?) Г Ф(т-0
Bn = 2B / cos Фп1 (Z )d( + / cos Фп2(С )dZ
0 0
где Фп1(С) = 2nn( - W(Z); Фп2(С) = 2nn( + W(Z).
Так как функции 2пп( и W(Z) монотонно возрастают, фаза Фп2(() быстро изменяется, что приводит к осциллированию с небольшой амплитудой соответствующих интегралов. Фаза же ФП1(С) может изменяться очень медленно. При этом наблюдается резонансное явление, которое характеризуется ростом интегралов, содержащих фазу Фп1(С). Из изложенного следует, что при возникновении резонанса рост амплитуды связан с возрастанием интегралов с фазой Фп1((), интегралами же с фазой Фп2(С) можно пренебречь. Тогда полная амплитуда, определяемая по формуле Ап(т) = А2п + Bп, в точке С = Со(т), соответствующей максимальному размаху колебаний, будет иметь вид
Ап(т)
4B2
г ь(т) 2 г ь(т)
/ cos Фп(С )d( + / sin Фп(С )dZ
0 0
где Ь(т) = ф(т - Со(т)), Фп(С) = 2пп( - W(Z).
149
Анисимов В.Н., Литвинов В.Л., Корпен И.В.
Возможность решения задачи (1)—(3) зависит от степени сложности граничных условий, а также возможности решения системы (10). Для решения таких систем в [5] использован обратный метод: по заданным <^(z) и ф(г) из получающейся системы уравнений находятся законы движения границ 1\(т) и I2(т). При задании функций <^(z) и ф(г) в них вводится несколько произвольных постоянных. Зависимость найденных законов движения 11 (т) и ^(т) от величин этих констант позволяет аппроксимировать достаточно разнообразные законы движения границ законами, полученными из решения обратной задачи.
Совокупность обратных решений достаточно широка [8]. Для примера приведём два таких решения:
1) для функций p(z) = B(e-az - 1) - C(e-a - 1) - 1, B = C + (e-a - 1)-1; ф(z) = C(eaz - 1) - C(e-a - 1) - 1 из системы (10) находятся следующие законы движения границ: 11(т) = a-1 ln[(Be-aT - CeaT)/(B - C)]; 12(т) = 1 + 11(т); здесь a, B, C, v — постоянные величины;
2) для функций
^(z) = ^(z)
Ln[{vz + 1)/(1 - v)} Ln[(l +v)/(l -v)]
находятся следующие законы движения границ: 11(т)=0; 12(т)=vт+1.
Первое из решений может быть использовано при изучении колебаний гибких звеньев передач [2]. Второе — при изучении колебаний канатов грузоподъемных установок при равномерном подъёме (спуске).
В заключение отметим, что приведённая в настоящей работе методика позволяет получить выражение для амплитуды колебаний, соответствующих n-ной динамической моде в случае граничных условий первого рода. Кроме того, описанный выше метод позволяет также рассмотреть следующий класс краевых условий [8]:
Ut (I1 (т), т) = F1 (т), Щ (12 (т), т) = ^2(т);
U (I1 (т),т)= F1 (т), Щ (I2 (т), т) = F2 (т); (22)
U (I1 (т), т) = F1 (т), Ut (I2 (т), т) = F2^).
Насколько известно, решения волнового уравнения при граничных условиях (22) другими известными аналитическими методами получить крайне сложно.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Савин Г.Н., Горошко О. А. Динамика нити переменной длины: Киев, 1962. 332 с. [Savin G. N. Goroshko O. A. Dynamics of variable length strands: Kiev, 1962. 332 pp.]
2. Самарин Ю.П., Анисимов В.Н. Вынужденные поперечные колебания гибкого звена при разгоне// Изв. вузов. Машиностроение, 1986. №12. С. 17-21. [Samarin Yu.P., Anisimov V. N. Forced transverse vibrations of the flexible link at dispersal // Izv. Vuzov. Mashinostroenie, 1986. no. 12. Pp. 17-21].
3. Анисимов В. Н., Литвинов В. Л. Исследование резонансных свойств механических объектов с движущимися границами при помощи метода Канторовича-Галёркина // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2009. №1(18). С. 149-158. [Anisimov V.N., Litvinov V. L. Investigation of resonance characteristics of mechanical objects with moving borders by application of the Kantorovich-Galyorkin method // Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2009. no. 1(18). Pp. 149-158].
150
Об одном методе получения аналитического решения волнового уравнения ...
4. Горошко О. А., Савин Г.Н. Введение в механику деформируемых одномерных тел переменной длины. Киев: Наукова думка, 1971. 270 с. [Goroshko O.A., Savin G. N. Introduction in mechanics of one dimensional deformable bodies of variable length. Kiev: Naukova Dumka, 1971. 270 pp.]
5. Весницкий А. И. Волны в системах с движущимися границами и нагрузками. М.: Физматлит, 2001. 320 с. [Vesnitskii A.I. Waves in systems with moving boundaries and loads. Fizmatlit: Moscow, 2001. 320 pp.]
6. Весницкий А. И. Обратная задача для одномерного резонатора изменяющего во времени свои размеры // Изв. вузов. Радиофизика, 1538-1542. №10; англ. пер.: Vesnitskii A. I. The inverse problem for a one-dimensional resonator the dimensions of which vary with time // Radiophys. Quantum Electronics, 1971. Vol. 14, no. 10. Pp. 1209-1215.
7. Барсуков К. А., Григорян Г. А. К теории волновода с подвижными границами// Изв.вузов. Радиофизика, 1976. №2. С. 280-285; англ. пер.: Barsukov K.A., Grigoryan G. A. Theory of a waveguide with moving boundaries // Radiophys. Quantum, Electronics, 1976. Vol. 19, no. 2. Pp. 194-200.
8. Анисимов В.Н., Литвинов В. Л. Резонансные свойства механических объектов с движущимися границами. Самара: СамГТУ, 2009. 131 с. [Anisimov V. N., Litvinov V. L. The resonance properties of mechanical objects with moving boundaries. Samara: SamGTU, 2009. 131 pp.]
Поступила в редакцию 10/V/2012; в окончательном варианте — 25/VII/2012.
MSC: 35L05; 35L04, 35Q70
ON A METHOD OF ANALYTICAL SOLUTION OF WAVE EQUATION DESCRIBING THE OSCILLATIONS SISTEM WITH MOVING BOUNDARIES
V. N. Anisimov, V. L. Litvinov, I. V. Korpen
Syzran’ Branch of Samara State Technical University,
45, Sovetskaya str., Syzran’, Samara region, 446001, Russia.
E-mails: [email protected], [email protected]
The method of analytical solution of wave equation with the conditions, assigned on the moving boundaries, is described. With the aid of the change of variables in the system of functional equations the original boundary-value problem is brought to the system of difference equations with one fixed bias, which can be solved using the Laplace integral transform. The expression for amplitude of oscillation corresponding to n-th dynamic mode is obtained for the first kind boundary conditions. This method makes it possible to examine the broader class of boundary conditions in comparison with other exact methods of solving the boundary-value problems with the moving boundaries.
Key words: wave equation, variations of systems with moving boundaries, laws of boundary moving, amplitude of oscillation.
Original article submitted 10/V/2012; revision submitted 25/VII/2012.
Valeriy N. Anisimov (Ph. D. (Phys. & Math.)), Associate Professor, Dept. of General-Theoretical Disciplines. Vladislav L. Litvinov, Teacher, Dept. of General-Theoretical Disciplines. Inna V. Korpen (Ph. D. (Educat.)), Associate Professor, Dept. of General-Theoretical Disciplines.
151