О задаче управления смещением одного из концов тонкого стержня Текст научной статьи по специальности «Математика»

12

Вестник Самарского университета. Естественнонаучная серия. Том 23 № 3 2017

УДК 517.95, 624.07

Б01: 10.18287/2541-7525-2017-23-3-12-17

А.Б. Бейлин1

О ЗАДАЧЕ УПРАВЛЕНИЯ СМЕЩЕНИЕМ ОДНОГО ИЗ КОНЦОВ

ТОНКОГО СТЕРЖНЯ

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

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

Введение

Исследованиям математических моделей колебательных процессов различной природы посвящено большое количество работ. Отметим здесь лишь некоторые из них, наиболее близкие к содержанию настоящей статьи [1—4]. Интерес к этой тематике вызван прежде всего тем, что колебательные процессы возникают в различных механических системах [5-8]. Для того чтобы обеспечить надежную работу сложной современной техники в процессе ее эксплуатации необходим контроль, который производят с помощью диагностических процедур. Однако приборы диагностики не всегда можно установить в непосредственном контакте с исследуемым объектом. Следовательно, возникает необходимость осуществлять диагностику на основе косвенных данных. С другой стороны, во избежание непредвиденных ситуаций нужно уметь управлять процессом и до проведения испытаний реального объекта знать допустимые интервалы значений входных данных. Например, как нужно закрепить концы стержня, чтобы его колебания не выводили систему из заданного режима. В связи с этим и возникают задачи об управлении колебаниями, которые также называют обратными задачами математической физики [9-12].

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

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

1. Постановка задачи.

В области Qт = (0,I) х (0, Т) рассмотрим задачу о колебании стержня, один конец которого, х = = 0, жестко закреплен, а закон движения второго конца, х = I, не задан и подлежит определению. Таким образом, мы приходим к обратной задаче, которая заключается в следующем: найти пару функций (и(х,Ь),Н(Ь)), обладающих свойствами

и е с1(^т) п с2^т), н е с1 [о,т] п с2(0,т),

!© Бейлин А.Б., 2017

Бейлин Александр Борисович ([email protected]), кафедра АСиИС, Самарский государственный технический университет, 443010, Российская Федерация, г. Самара, ул. Молодогвардейская, 133.

удовлетворяющих уравнению

Lu = utt - a2uxx = f (x,t), (1.1)

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

u(x, 0)=0, ut(x, 0) = 0, (1.2)

краевым условиям

u(0,t)=0, ux(l,t) = h(t) (1.3)

и условию переопределения, в качестве которого часто задается интегральное среднее искомого решения [10; 11]

i

¡K (x)u(x-t)dx =E (t)■ (14)

0

Заметим, что однородность начальных данных не ограничивает общности, но упрощает многие преобразования.

Задача (1.1)—(1.4) может быть трактована как задача управления, а именно, как задача определения таких условий на входные данные, а именно, функцию h(t), позволяющую сохранить заданную энергию E(t).

2. Разрешимость задачи.

Будем предполагать, что выполнены следующие условия:

a > 0, f,ftt € C(Qt), K € C[0, /], K(l)=0, K(0) = 0,E € C2 [0,T], а также условия согласования E(0) = E'(0) = 0. Очевидно, что искомая функция h(t) должна удовлетворять условиям h(0) = h'(0) = 0.

Прежде всего получим соотношение между h(t) и u(x,t), которое в дальнейшем и приведет к желаемому результату.

Пусть (u,h) — решение задачи (1.1)—(1.4), K(l) = 0 Проинтегрируем равенство (1.1), умноженное на K(x), по промежутку (0,1). Учитывая условия (1.3) и (1.4), получим после элементарных преобразований

E''(t) K'(l) „ .

h(t) = aét) + m u(l,t)-i i

- aKt)j K ''(x)u(x,t)dx — а22ЩТ)1 K (x)f (x,t)dx. (2.5)

0 0 Теперь видно, что найдя u(x,t) как решение прямой задачи (1.1)—(1.3), где h(t) считаем временно известной, мы можем рассчитывать, что найдем h(t) из соотношения (2.5). Нашей ближайшей целью является обоснование возможности реализовать описанную схему и получить таким образом решение задачи (1.1)—(1.4). В качестве первого шага мы рассмотрим прямую задачу (1.1)—(1.3). Так как второе из краевых условий (2.5) неоднородно, то для того чтобы воспользоваться методом разделения переменных перейдем к задаче с однородными условиями, введя новую неизвестную функцию

v(x, t) = u(x, t) — xh(t).

Действительно, теперь нам нужно найти функцию v(x,t) как решение следующей задачи

vtt — a2vxx = f (x,t) — xh''(t), (2.6)

v(x, 0) = 0, vt(x, 0) = 0, (2.7)

v(0,t) = 0, vx(0,t)=0. (2.8)

Применив метод разделения переменных, полагая v(x,t) = X(x)T(t), легко находим собственные функции соответствующей задачи Штурма-Лиувилля:

• n(l + 2n)x Xn(x) = sin---.

Так как начальные и граничные условия нулевые, то единственное решение однородного уравнения тривиально. Поэтому сразу переходим к нахождению решения неоднородного уравнения (2.6). Введем еще некоторые упрощения исключительно для наглядности и не ограничивая общности, а именно, пусть f (x,t) = = 0, K(x) = x2. Итак, решение задачи (2.6)—(2.8) ищем в виде

v(x,t) = ¿ Vn(t)sin n(1+2^x, (2.9)

n=0

где Уп(г) подлежат определению и должны удовлетворять условиям

Уп(0) = УП(0) = 0

для того чтобы начальные условия (2.7) были выполнены для у(х,г). Заметим, что выполнение граничных условий гарантируют собственные функции Хп(х). Подставив (2.9) в (2.6), приходим к задаче Коши:

УП (г) + и>П Уп(г) = Нп(г), Уп(0) = УП(0) = 0, (2.10)

где шп = п(1+^п)а, нп(г)- коэффициенты Фурье функции -хН''(г). Решение этой задачи, как известно, выражается формулой

Vn(t) = — í hn(r)sin wn(t — т)dr.

^n J

(2.11)

Найдем коэффициенты Фурье функции —xh''(t):

, .. ,,, 2 Г . n(1 + 2n)x ( — 1)n+l4l2h''(t)

hn(t) = —h''(t) lj x Sin^2^dx = J(1 + 2n)2 ■

Теперь решение задачи Коши можно записать так:

Vn(t)

l (—l)n+4l2

n2(1 + 2n)2

t

J h''(т) sinZn(t — т)dT.

(2.12)

Тогда, подставив в полученную формулу шп = п(1+2п)а, получим решение задачи (2.6)—(2.8)

, . 813 ^ ( — lYf,,,,^. , ч, • n(1 + 2n)x

,(x>t) = (1 + 2n)3 h (T)Sin Zn(t — T)dT Sin

n=0 ^ '

2l

(2.13)

Сходимость этого ряда не вызывает сомнений. Однако мы по-прежнему рассматриваем функцию h(t) как известную. Приступим теперь к нахождению этой функции, применив к решению (2.13) условие переопределения, представленное в виде (2.5). Предварительно сделаем некоторые преобразования. Проинтегрировав дважды по частям правую часть формулы (2.12), получим

t t

j h''(т) sinшп(t — т)dT = znh(t) — шП j h(T) sinzn(t — т)dT.

Тогда (2.13) примет вид

16l2 ~ ( — 1)n+1 >(x>t) = шnh(t)—

—u>2J h(r) sinzn(t — т)dr\ sin

. n(1 + 2n)

2l

(2.14)

Обратимся к формуле (2.5) и подставим в нее найденное решение и(х,г) = у(х,г) + хН(г) где у(х,г) определяется равенством (2.14). После элементарных преобразований получим

1 — a2,,, E''(t) 2

—— h(t) = К) '

2

a2

a

2l2

i

22

+ -jv(l,t) — aip v(x,t)dx.

(2.15)

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

Пусть Т < ^. Тогда можно эффективно воспользоваться следующими формулами [14]:

Е-

1

^ (—1)к+1 П

= —, (0.234(2, 4)),

(2k — 1)2 8' (2k — 1)3 32

sin(2k — 1)x n . ..

2k_i =4, 0 <x<n (^m))

k=i

t

2

n

,1, sin(2k + 1)x п п . . ..

i-1)k ,+J = -jx, 0 < x < -, (1.444(5)).

^ (2k + 1)2 4 ' ^ ^ 2

Учтя эти формулы и вычислив /0 sin 7T^1+22n)x dx = п(1+2п) , получим

t

v(l,t) = a J h(т d - lh(t),

-ait — т) п a .

sin —--- H--(t - т)

2l + 8l ( )

dт.

j v(x,t)dx = -—h(t) + ~—r j ^т) 00 Подставив полученные выражения в (2.14), приходим к уравнению Вольтерра второго рода

t

h(t) - J H (t, т )h(т d = g(t), (2.16)

0

где обозначено

rr, . 16 . -a(t — т) 2 . .

H) =a+a-sin ~2l +12(t-т),

(t) E''(t)

g(t) = aw ■

Уравнение (2.16) в силу ограниченности ядра однозначно разрешимо [15]. Существование производных функции h(t) вытекает из свойств E(t) и гладкости ядра.

Подводя итоги проведенным исследованиям, сформулируем полученный результат. Если выполнены условия

a> 0, f,ftt G С(Qt),

K G С[0,l], K(l) = 0, K(0)=0,E G С3[0,Т] П С4(0,Т),

а также условия согласования E(0) = E'(0) =0, то в области Qt существует единственное решение задачи (1.1)—(1.4). В частном случае, f = 0, K (x) = x2, T ^ a, интегральное уравнение, определяющее режим смещений правого конца стержня, принимает простой вид, пригодный для применения приближенных методов его решения.

Литература

[1] Rao J.S. Advanced Theory of Vibration. N.Y.: Wiley, 1992. 431 с.

[2] Федотов И.А., Полянин А.Д., Шаталов М.Ю. Теория свободных и вынужденных колебаний твердого стержня, основанная на модели Рэлея // Доклады РАН. 2007. T. 417. № 1. С. 56-61.

[3] Бейлин А.Б., Пулькина Л.С. Задача о продольных колебаниях стержня с динамическими граничными условиями // Вестник СамГУ. Естественнонаучная серия. 2014. № 3(114). С. 9-19.

[4] Бейлин А.Б. Задача о продольных колебаниях упруго закреплённого нагруженного стержня // Вестник Самарского гос. Тех. Ун-та. Серия: Физ.-мат. науки. 2016. Т. 20. № 2. С. 249-258.

[5] Бабаков И.М. Теория колебаний. М.: Наука, 1968. 560 с.

[6] Хазанов Х.С. Механические колебания систем с распределёнными параметрами: Учеб. пособие. Самара: Са-мар. Госуд. Аэрокосмич. Ун-т, 2002. 80 с.

[7] Вейц В.Л., Дондошанский В.К., Чиряев В.И. Вынужденные колебания в металлорежущих станках. М-Л.: Машгиз, 1959 288 с.

[8] Кумабэ Д. Вибрационное резание. М.: Машиностроение, 1985. 424 с.

[9] Ильин В.А., Тихомиров В.В. Волновое уравнение с граничным управлением на двух концах и задача о полном успокоении колебательного процесса // Дифференц. уравнения. 1990. Т. 35. № 5. С. 692-704.

[10] Камынин В.Л. Обратная задача определения младшего коэффициента в параболическом уравнении при условии интегрального наблюдения // Матем. заметки. 2013. 94(2). C. 207-217.

[11] Cannon J.R., Lin Y. Determination of a parameter p(t) in some quasi-linear parabolic differential equations. Inverse Problems. 1988. № 4. P. 35-45.

[12] Денисов А.М. Обратная задача для гиперболического уравнения с нелокальным краевым условием, содержащим запаздывающий аргумент // Труды института математики и механики УрО РАН. 2012. Т. 18. № 1.

[13] Бейлин А.Б., Пулькина Л.С. Задача о колебаниях стержня с неизвестным режимом на части границы // Вестник Самарского университета. Естественнонаучная серия. 2017. Т. 2.

[14] Градштейн B.C., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: Физматгиз. 1963. 1100 с.

[15] Михлин С.Г. Лекции по линейным интегральным уравнениям. М.: Физматгиз, 1959. 232 с.

References

[1] Rao J.S. Advanced Theory of Vibration. N.Y.: Wiley, 1992, 431 p. [in Russian].

[2] Fedotov I.A., Polyanin A.D., Shatalov M.Yu. Teoriia svobodnykh i vynuzhdennykh kolebanii tverdogo sterzhnia, osnovannaia na modeli Releia [Theory of free vibration of rigid rod based on Rayleigh model]. Doklady RAN [Proceedings of the Russian Academy of Sciences], 2007, Vol. 417, no. 1, pp. 56-61 [in Russian].

[3] Beylin A.B., Pulkina L.S. Zadacha o prodol'nykh kolebaniiakh sterzhnia s dinamicheskimi granichnymi usloviiami [A Problem on Longitudinal Vibration in a Short Bar with Dynamical Boundary Conditions]. Vestnik SamGU. Estestvennonauchnaia seriia [Vestnik of Samara State University. Natural Science Series], 2014, no. 3(114), pp. 9-19 [in Russian].

[4] Beylin A.B. Zadacha o prodol'nykh kolebaniiakh uprugo zakreplennogo nagruzhennogo sterzhnia [Task on longitudinal vibration of an elastic attached loaded rod]. Vestnik Samarskogo gos. Tekh. Un-ta. Seriia: Fiz.-mat. nauki [Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences], 2016, Vol. 20, no. 2, pp. 249-258 [in Russian].

[5] Babakov I.M. Teoriia kolebanii [Theory of vibrations]. M.: Nauka, 1968, 560 p. [in Russian].

[6] Khazanov Kh.S. Mekhanicheskie kolebaniia sistem s raspredelennymi parametrami: Ucheb. posobie [Mechanical oscillations of systems with distributed parameters: textbook]. Samara: Samar. Gosud. Aerokosmich. Un-t, 2002, 80 p. [in Russian].

[7] Veiz V.L., Dondoshansky V.K., Chiriaev V.I. Vynuzhdennye kolebaniia v metallorezhushchikh stankakh [Forced oscillations in metal-cutting machines]. M-L.: Mashgiz, 1959, 288 p. [in Russian].

[8] Kumabe D. Vibratsionnoe rezanie [Vibration Cutting]. M.: Mashinostroenie, 1985, 424 p. [in Russian].

[9] Ilin V.A., Tikhomirov V.V. Volnovoe uravnenie s granichnym upravleniem na dvukh kontsakh i zadacha o polnom uspokoenii kolebatel'nogo protsessa [The wave equation with boundary control at two ends and the problem of complete damping of a vibration process]. Differents. uravneniia [Differential Equations], 1990, Vol. 35, no. 5, pp. 692-704 [in Russian].

[10] Kamynin V.L. Obratnaia zadacha opredeleniia mladshego koeffitsienta v parabolicheskom uravnenii pri uslovii integral'nogo nabliudeniia [The inverse problem of determining the lower-order coefficient in parabolic equation with integral observation]. Matem. zametki [Mathematical notes], 2013, 94(2), pp. 207-217 [in Russian].

[11] Cannon J.R., Lin Y. Determination of a parameter p(t) in some quasi-linear parabolic differential equations. Inverse Problems, 1988, no. 4, pp. 35-45 [in English].

[12] Denisov A.M. Obratnaia zadacha dlia giperbolicheskogo uravneniia s nelokal'nym kraevym usloviem, soderzhashchim zapazdyvaiushchii argument [The inverse problem for a hyperbolic equation with nonlocal boundary condition involving retarding argument]. Trudy instituta matematiki i mekhaniki UrO RAN [Proceedings of Institute of Mathematics and Mechanics. Ural Branch of Russian Academy of Sciences], 2012, Vol. 18, no. 1 [in Russian].

[13] Beylin A.B., Pulkina L.S. Zadacha o kolebaniiakh sterzhnia s neizvestnym rezhimom na chasti granitsy [On Certain Problem for Hyperbolic Equation with unknown behavior on a part of the boundary]. Vestnik SamGU. Estestvennonauchnaia seriia [Vestnik of Samara State University. Natural Science Series], 2017, no. 2, pp. 9-19 [in Russian].

[14] Gradstein V.S., Ryzhik I.M. Tablitsy integralov, summ, riadov i proizvedenii [Tables of integrals, sums, series, and products]. M: Fizmatgiz, 1963, 1100 p. [in Russian].

[15] Mikhlin S.G. Lektsii po lineinym integral'nym uravneniiam [Lectures on linear integral equations]. М.: Fizmatgiz, 1959, 232 p. [in Russian].

A.B. Beylin2

ON CERTAIN CONTROL PROBLEM OF DISPLACEMENT AT ONE

ENDPOINT OF A THIN BAR

In this paper, we study an inverse problem for hyperbolic equation. This problem arises when we consider vibration of a thin bar if one endpoint is fixed but behavior of the other is unknown and is the subject to find. Overdetermination is given in the form of integral with respect to spacial variable. The problem is reduced to the second kind Volterra integral equation. Special case is considered.

Key words: hyperbolic equation, vibration of a thin bar, inverse problem, integral overdetermination.

Статья поступила в редакцию 28/V/2017. The article received 28/V/2017.

2 Beylin Alexander Borisovich ([email protected]), Department of Automated Machining and Tool Systems, Samara State Technical University, 133, Molodogvardeiskaya str., Samara, 443010, Russian Federation.