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

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2015. T. 19, № 2. С. 382—397

ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1330

Математическое моделирование, численные методы и комплексы программ

УДК 517.958:531.12; 534.11

МАТЕМАТИЧЕСКИЕ МОДЕЛИ НЕЛИНЕЙНЫХ ПРОДОЛЬНО-ПОПЕРЕЧНЫХ КОЛЕБАНИЙ ОБЪЕКТОВ С ДВИЖУЩИМИСЯ ГРАНИЦАМИ

В. Н. Анисимов, В. Л. Литвинов

Сызранский филиал Самарского государственного технического университета,

Россия, 446001, Самарская обл., Сызрань, ул. Советская, 45.

Аннотация

Произведены нелинейные постановки задач, описывающих продольнопоперечные колебания объектов с движущимися границами. Полученные математические модели состоят из системы двух нелинейных дифференциальных уравнений в частных производных с наибольшей производной по времени второго порядка и по пространственной переменной — четвёртого порядка. Нелинейные условия на движущейся границе имеют максимальную производную по времени второго порядка и по пространственной переменной третьего порядка. Учтены геометрическая нелинейность, вязкоупругость, изгибная жёсткость колеблющегося объекта, а также упругость подложки, на которой расположен объект. Получены граничные условия в случае наличия энергетического обмена между частями объекта слева и справа от движущейся границы. Движущаяся граница имеет присоединённую массу. Учтён упругий характер присоединения границы. С помощью полученной математической модели описываются продольно-поперечные колебания большой интенсивности объектов с движущимися границами. При получении математических моделей использован вариационный принцип Гамильтона.

Ключевые слова: продольно-поперечные колебания, движущиеся границы, краевые задачи, математические модели, граничные условия, нелинейная система в частных производных, вариационные принципы. doi: http://dx.doi.org/10.14498/vsgtu1330

© 2015 Самарский государственный технический университет.

Образец для цитирования

Анисимов В. Н., Литвинов В. Л. Математические модели нелинейных продольнопоперечных колебаний объектов с движущимися границами // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2015. T. 19, № 2. С. 382-397. doi: 10.14498/vsgtu1330. Сведения об авторах

Валерий Николаевич Анисимов (к.ф.-м.н., доц.; anisimov170159@mail.ru; автор, ведущий переписку), заведующий кафедрой, каф. общетеоретических дисциплин.

Владислав Львович Литвинов (vladlitvinov@rambler.ru), старший преподаватель, каф. общетеоретических дисциплин.

382

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

Колеблющиеся механические объекты с движущимися границами широко распространены в технике. Это канаты грузоподъёмных установок [1,2], гибкие звенья передач [3], лентопротяжные механизмы [4], балки [5,6], рельсовые пути [7,8], бурильные колонны [9], струны музыкальных инструментов [10], нити [11], тесёмочные передачи и т. д. До настоящего времени задачи о колебаниях таких систем решались в основном при линейной постановке и жёстком закреплении границ, когда отсутствует энергетический обмен через границу [1,3,5,12-17]. В редких случаях учитывалось действие демпфирующих сил. Как показано в работе [12], при описании колебаний волновым уравнением в случаях, когда скорость движения границы стремится к скорости распространения волн и при стремлении длины объекта к нулю энергия колебаний стремится к бесконечности, чего в реальных объектах быть не может. Возникла необходимость в построении более сложных математических моделей продольно-поперечных колебаний объектов с движущимися границами, учитывающих большое число факторов, влияющих на колебательный процесс. В данной статье построен ряд нелинейных математических моделей такого рода задач. В отличие от существующих моделей получены граничные условия в случае наличия энергетического обмена между частями объекта слева и справа от движущейся границы, учтено взаимодействие между продольными и поперечными колебаниями. При получении математических моделей использован вариационный принцип Гамильтона. Полученные математические модели позволяют описывать колебания большой интенсивности систем с движущимися границами.

1. Описание объекта. В настоящей работе рассматриваются объекты, обобщённая схема которых схематически изображена на рис. 1. Объекты могут быть струной, стержнем, балкой, канатом и т. д.

Для объекта введены следующие обозначения: р — объёмная плотность массы; S — площадь поперечного сечения; I — осевой момент инерции поперечного сечения объекта; E — модуль упругости материала объекта; р — коэффициент, характеризующий свойство вязкоупругости объекта на основе структурной модели Фойгта; е0 — начальная продольная деформация объекта, создающая натяжение T = ESeо; x — расстояние от левой границы до точки объекта, находящегося в недеформированном состоянии; L(t) —длина недеформированного в продольном направлении объекта слева от движу-

Рис. 1. Обобщённая схема объекта [Figure 1. Generalized sketch of object]

383

Анисимов В. Н., Литвинов В. Л.

щейся границы; t <Е [0,ti] —время; L0 — общая длина недеформированного объекта.

Для характеристики окружения объекта введены следующие параметры: ко — жёсткость подложки, на которой лежит объект (сила, действующая на единицу длины при единичном поперечном смещении); V(t) —окружная скорость роликов; Л — коэффициент, характеризующий действие сил сопротивления внешней среды (силы сопротивления пропорциональны скорости поперечного движения и длине). На объект в направлении вектора ё2 действует распределённая нагрузка f (x,t). На движущуюся границу действует сила F(t) = Fi(t)ei + F2(t)e2.

Движущаяся граница состоит из жёстко соединённых роликов массой m2. Масса системы роликов и каркаса равна mi. Пружина жёсткости к2 реагирует на поперечное смещение системы роликов. В продольном направлении имеет место жёсткое соединение между системой роликов и каркасом. Между роликами и объектом проскальзывание отсутствует, поэтому при предположении малости продольных деформаций имеет место равенство L'(t) = V(t).

Для характеристики продольно-поперечных колебаний объекта введём функции: u1(x,t), u2(x,t) —смещения точек объекта с координатой x в момент времени t в направлении базисных векторов ei, в2 соответственно.

Для сокращения записей введём новую функцию

u1(x,t)= u1(x,t)+ x + e0x. (1)

Функция ui(x,t) представляет собой расстояние от левой границы до точки с координатой x деформированного объекта.

2. Геометрические и дифференциальные характеристики объекта. Рассмотрим геометрию объекта, то есть найдем изменение длины элемента dx и его кривизну.

Элементарный вектор dl, в который перейдет элемент dx (рис. 1), в момент времени t равен

dl = (ui,x + 1 + £o)dxei + u2,x(x,t)dxe2.

Здесь и далее под u| и u2 понимаются функции u\(x,t) и u2(x,t). После запятой указана переменная, по которой берётся частная производная. Длина вектора находится по формуле

\Ш\ = ^(1+ £о + ui,x)2 + u2,x dx.

Найдём относительную деформацию элемента:

, \dl\ — dx г — - , s

£(x, t) = --- = \/(1 + £0 + u*l ,x)2 + u2,x — 1. (2)

С учётом (1) формула (2) будет иметь вид

£(x,t) — ^/vj7xvjTx 1. (3)

Здесь и далее по повторяющимся индексам предусматривается свёртывание.

384

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

В предположении малости изгибных деформаций кривизна элемента dI (см. рис. 1) находится по формуле

K (x,t) = U2,xx. (4)

3. Описание областей колебаний и соотношений на границе. Для учёта энергетического обмена, происходящего на границе, разобьём область колебаний в координатах x ~ t на две открытые области Wi, W2 (рис. 2), которые соответствуют частям объекта слева и справа от движущейся границы. Через W обозначена объединённая область, состоящая из Wi и W2. Через 71 и 72 обозначены замкнутые контуры, окружающие области (положительное направление обхода против часовой стрелки). Обозначим через u(L(t) — 0,t), Ui(L(t) + 0,t) значения функций слева и справа от движущейся границы.

Рис. 2. Области колебаний [Figure 2. Oscillations domains]

Из условия неразрывности деформаций объекта следует

Ui(L(t) — 0,t) = Ui(L(t) + 0,t), i = 1, 2. (5)

Характер взаимодействия роликов и объекта (см. рис. 1) исключает угловое смещение в точке их контакта, поэтому при наличии изгибной жёсткости имеет место равенство

U2,x(L(t) — 0, t) = U2,x(L(t) + 0, t) = 0. (6)

Продифференцируем соотношение (5) по t:

Ui,x(L(t) — 0, t)L'(t) + Ui,t(L(t) — 0, t) = Ui,x(L(t) + 0, t)L'(t) + Ui,t(L(t) + 0, t).

С учётом (6) полученное выражение преобразуется к виду

U2,t(L(t) — 0, t) = U2,t(L(t) + 0 ,t) = 0. (7)

4. Нахождение составляющих вариации интеграла действия. Для получения модели, описывающей продольно-поперечные колебания объекта, предлагается использовать вариационный принцип Гамильтона [19]: если заданы начальное и конечное состояния системы, то из всех возможных законов движения реализуется такой, для которого действие

Г11

I = (T — U) dt

Jto

385

Анисимов В. Н., Литвинов В. Л.

принимает стационарное значение (T — кинетическая энергия системы, U — потенциальная энергия). Для применения принципа понадобятся вариации интеграла действия от кинетической и потенциальной энергии системы.

4.1. Вариация от кинетической энергии объекта. Интеграл действия от кинетической энергии объекта определяется формулой

1 rti rLo

It 1 = т dt pSuj,tUj,t dx.

2 J 0 J 0

Представим двукратный интеграл в виде двойного:

Lt i

Uj,tUj,t dW

и найдём вариацию:

6Iti = pS ([[ Uj,tSuj,t dW +[[ Uj,tSuj,t dW^ . (8)

\J JWx J JW2 J

Здесь интеграл по области W представлен как сумма интегралов по областям Wi и W2.

Представим подынтегральное выражение в виде

Uk,t&Uk

д_

dt

{Uk,t^Uk) Uk,tt&Uk.

С помощью формулы Грина преобразуем выражение (8):

5It 1 = pS I Ф Uk,t^Uk dx + Ф Uk,t^Uk dx - Uk,tt^Uk dW

\Jyi «/ 72 «/ Jw

(9)

Здесь yi и Y2 — замкнутые контуры, окружающие области Wi и W2.

4.2. Вариация интеграла действия от кинетической энергии системы роликов и каркаса. Выразим скорости движения границы в продольном и поперечном направлениях:

Vi(t) = —Ui(L(t),t), i = 1, 2.

Здесь и далее (d/dt) — полная производная. Кинетическая энергия движущейся границы определяется выражением

2 (mi Vf(t) + m,2V>2(t)) .

Интеграл действия от кинетической энергии системы роликов и каркаса

d

1 ftl d d

It2 = 2mi J dtUi(L(t),t)—Ui(L(t),t)dt.

386

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

Для вариации действия получим

ftl d d / | \

5It2 = w* J —Ui(L(t),t)dt: vUilx=m)dt• (10)

Представим подынтегральное выражение из (10) в виде

Ли‘(ад'^Н'-Ч«)

А (d

dt V dt

Ui(L(t),t)Sui

x-L(t)

dt2

Ui(L(t),t)Sui

x-L(t)'

В результате для вариации (10) получим

SIT2 = wi(dl^AWAHx-Lw) L-t - m{dtUi(L(t),t)^Ui|x-L(t))

dt

tl d2 |

J dtp Ui(L(t) > t)6Ui |x-L(t) dt (11)

t-0

4.3. Вариация интеграла действия от упругой потенциальной энергии и вязкоупругих сил при продольной деформации объекта. Для учёта зависимости между деформациями и напряжениями используется модель Фойгта

a(t) = Ee(t) + ^e(t), (12)

где a(t) — напряжение; E — модуль упругости; e(t) — относительная деформация; у — коэффициент вязкоупругости; e(t) — производная по времени от e(t). Упругая потенциальная энергия находится по формуле

U1 = 1 ES f ( e2(x,t)dW,

2 J Jw

где e(x,t) — определяется выражением (3).

Вариация интеграла действия представляет собой изменение энергии при изменении функций ui, U2. Вариация интеграла действия от упругих и демпфирующих сил задаётся соотношением

61ц 1 = S (Ee(x, t) + yet(x, t)) 6e(x, t) dW.

J Jw

После преобразований (13) примет вид

61ц i = S

E (VUj5xUj,x l) + Ц-

Uj,xUj,xt

Ui,x

W \ 4 y/Uj,xUj,x' y/Uj,xUj,x

Представим подынтегральное выражение из (14) в виде

Uj,xUj,xt \ U

E (VUj7^Uj7x !-) + у

jjx' дщ

ix

Uj,xU

-6Ui,x —

j x j x

(13)

6Ui,x dW.

(14)

387

Анисимов В. Н., Литвинов В. Л.

д (/ , ,_______ ч Uj,xUj,xt \ ui,x5ui

= q^((E {Vjj - U jX ^ХМ ^X i

л/и j,XUj,X' у/u j,xuj,x

д

dx\\E (Vuj5xuj5^ l) + U

uj,xuj,xt \ ui,x

yjuj,xuj,x' yfuj,xuj,x

Sui. (15)

После подстановки (15) в (14) и использовании формулы Грина для областей W\ и W2 выражение (14) примет вид

д

SIU 1— S I lw dx \ (E iVuj,xuj,x l) +u

uj,xuj,xt \ u

ix

y/uj7xuj7x' y/uj,xuj,x

Sui dW—

- Si ((E (VWj - l) + u Uj’‘Uj”‘ \ '“'x Sui I d t-

Y1

yjuj,xuj,x ' yfuj,xuj,x

- Si llE (jjx - l) + u "jj \ u:xSu I dt. (16)

Y2

y/uj,xuj,x ' yfuj,xuj,x

4.4. Вариация интеграла действия от потенциальной энергии и вязкоупругих сил для деформаций изгиба объекта. При описании состояния объекта с помощью модели Фойгта (12) изгибающий момент в сечении записывается следующим образом:

M — I (EK(x, t) + uKt(x, t)).

Здесь: I — момент инерции поперечного сечения объекта, K(x,t) — кривизна объекта, которая определяется выражением (4). Интеграл действия энергии упругой деформации определяется следующим выражением:

IU2 — lEI JJwu2xx dW.

Вариация интеграла действия упругой деформации и вязкоупругих сил имеет вид

IU2 — 1 I I (Eu2,xxSu2,xx + uu2,xxtSu2,xx) dW. (17)

W

Выполним следующие преобразования подынтегральных составляющих в (17):

дд

u2,xxSu2,xx — "Ц (u2,xxSu2,x) Ц (u2,xxxSu2) + u2,xxxxSu2; (18)

dx dx

дд

u2,xxtSu2,xx — "Ц (u2,xxtSu2,x) Ц (u2,xxxtSu2) + u2,xxxxtSu2. (19)

dx dx

С помощью формулы Грина, учитывая (18), (19), для (17) получим

IU2 — I I (Eu2 ,xxxx + uu2,xxxxt) Su2 dW+

W

+ I Ф i.(Eu2,xxx + uu2,xxxt) Su2 (Eu2,xx + uu2,xxt) Su2,x) dt+

J71

388

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

+ I Ф ((Eu2,xxx + Eu2,xxxt) Su2 (Eu2,xx + ^u2,xxt) Su2,x) dt (20)

"'72

4.5. Вариации интегралов действия от факторов внешнего окружения объекта. Вариация интеграла действия от упругой деформации подложки задаётся соотношением

SIu з = ko

u2 5u2 d W.

W

(21)

Вариация интеграла действия от упругого присоединения системы роликов выражается соотношением

SIU 4

k2

ft1

I U2(L(t),t)SU2

o

Tt,, dt.

x=L(t)

Вариация интеграла действия внешних сил сопротивления

(22)

SIu 5 = А

W

U2,t SU2 dW,

(23)

Вариация интеграла действия от внешней распределенной нагрузки задаётся соотношением

SIu 6

I f (x, t) Su2 dW,

W

(24)

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

ill

SIu7 = — / Fk Suk

o

I x=L(t)

dt.

(25)

5. Получение системы дифференциальных уравнений, граничных и начальных условий, описывающих колебания объекта. Согласно вариационному принципу Гамильтона, вариация интеграла действия должна быть равна нулю:

SI = —SIti — SIt2 + SIui + SIu2 + SIu3 + SIu4 + SIu5 + SIu6 + SIu7 = 0. (26)

Составляющие вариации определены выражениями (9), (11), (16), (20)—(25).

5.1. Система дифференциальных уравнений в частных производных.

Приравняем нулю коэффициенты перед Sui, Su2 составляющих вариации (26), содержащих двойные интегралы. В результате получим систему дифференциальных уравнений, описывающую продольно-поперечные колебания объекта:

д

pSui,tt S дх ' 'E (/uj?xu j,x l) + E

ui,x

uj,xuj,xt

sjuj,xuj,x' y/uj,xuj,x

= 0; (27)

n n д ((ти ( I---------- 1^1 uj,xuj,xt \ u2,x \

pSu2,tt — S— ( [E — ^ + Ц u, J u J

д x uj xuj x uj xuj x

+

389

Анисимов В. Н., Литвинов В. Л.

+ I (Eu2,xxxx + M-U2,xxxxt) + ^u2,t + K0u2 — /(x — 0. (28)

5.2. Граничные условия при x — 0, x — Lo. Граничные условия на неподвижных границах (см. рис. 1) имеют классический вид:

U^(0, t) — 0, U2(0, t) — 0, U2,xx(0, t) — 0, u);(Lo, t) — 0, U2(Lo, t) — 0, U2,x(Lo, t) — 0.

(29)

5.3. Граничные условия при x — L(t). Выпишем части вариаций по замкнутым контурам Yi и 72 выражений (9), (11), (16), (20), (22), (25), соответствующих линии x — L(t):

J ^oS(ufc,t(L(t) - 0,t) - Ufc,t(L(t) + 0,t))L'(t)6ufc|x=L(t) +

d2 i

+ mi dt2 Ui (L(t), t)^Uj|x=L(t) +

+ ((ES(^Uj,x(L(t) - 0,t)Uj,x(L(t) - 0,t) - 1) +

Ufc,x(L(t) - 0, t)Ufc,xt(L(t) - 0, t) /uj,x(L(t) - 0,t)Uj,x(L(t) - 0,t) Ui,x(L(t) - 0,t)

+ ^S

y/Uj,x(L(t) - 0,t)Uj,x(L(t) - 0,t)

- (^E^0\/Ui5^0yt) - ^ +

+ ^S

Ufc,x(L(t) + 0, t)Ufc,xt(L(t) + 0, t) / Uj,x(L(t) + 0,t)Uj,x(L(t) + 0, t) Ui,x(L(t) +0,t)

^Ui

/uj,x(L(t) +0,t)Uj,x(L(t) + 0,t)' i x L(t)

- I (EU2 , xxx (L(t) - 0, t) + ^U2 , xxxt (L(t) - 0,t) -

- Eu2,xxx(L(t) + 0, t) - ^U2,xxxt(L(t) + 0, 'O)^U2 |x=L(t) +

+ I(EU2,xx(L(t) - 0, t) + ^U2,xxt(L(t) - 0, t) -

- EU2,xx(L(t) +0,t) - ^U2,xxt(L(t) + 0,t))^U2,x|x=L(t) +

+ K2U2(L(t),t)^U2|x=L(t) - Efc(t)^Ufc|x=L(t^ dt — 0. (30)

Приравнивая коэффициенты перед вариациями к нулю и учитывая соотношения (5)—(7), получим

pS(ui,t(L(t) - 0, t) - ui,t(L(t) + 0, t))L'(t) + miui(L(t), t) +

+ (ES(^Uj,x(L(t) - 0, t)uj,x(L(t) - 0,t) - 1) +

+ Ufc,x(L(t) - 0, t)ufc,xt(L(t) - 0, t) \

^ /Uj,x(L(t) - 0, t)uj,x(L(t) - 0, t)'

390

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

v _________Ui,x(L(t) - 0,t)________

/uj,x(L(t) - 0,t)uj,x(L(t) - 0,t)

— H —1)+

+ pS

Uk,x(L(t) + 0, t)ufc,xt(L(t) + 0, t) \ ^ VUj,X(L(tyTo7t)Uj,X^^(ty+o7t^ /

V ________ui,x(L(t) + 0,t)_______

\/ uj,x(L(t) +0,t)uj,x (L(t) +0,t)

Fi(t) = 0;

(31)

d2

m2dt2u2(L(t),t) + EIu2,xxx(L(t) + 0, t) + ju/u2,xxxt(L(t) + 0, t)-EIu2,xxx(L(t) 0, t) PIu2,xxxt(L(t) 0, t) +

+ K2u2(L(t),t) - F2(t) = 0. (32)

К полученным граничным условиям на границе необходимо добавить соотношения

u1(L(t) - 0, t) = u1(L(t) + 0, t);

u,2(L(t) - 0, t) = u,2(L(t) + 0,t); (33)

u2,x(L(t) - 0, t) = u2,x(L(t) + 0,t) = 0.

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

m2dt2u2(L(t),t) + \ ES(^uj,x(L(t) - 0,t)uj,x(L(t) - 0,t) - l) +

+ S uk,x(L(t) - 0, t)uk,xt(L(t) - 0,t) \

^ /uj,x(L(t) - 0, t)uj,x(L(t) - 0,ty

v ________u2,x(L(t) - 0,t)_________

/uj,x(L(t) + 0, t)uj,x(L(t) + 0,t)

- (E^0\Zuj,x(L^y"+oyt)uj,x^L(ty"+oyt) -1)+

+ rS ufc,x(L(t) + 0, t)uk,xt(L(t) + 0, t) \

^ /uj,x(L(t) + 0, t)uj,x(L(t) + 0, t)'

u2,x(L(t) +0,t) ,

X ---------------------------------+

y/uj,x(L(t) + 0, t)uj,x(L(t) + 0,t)

+ pS(u2,t(L(t) - 0, t) - u2,t(L(t) + 0, t))L'(t) +

+ Ku2(L(t),t) - F2(t) = 0.

Полученные граничные условия называются естественными [19].

В выражении (30) присутствуют вариации

^uilx=L(t), ^u2lx=L(t), $u2,x lx=L(t).

391

Анисимов В. Н., Литвинов В. Л.

Если функции ui(x,t), u2(x,t) заданы, то их вариации равны нулю. Граничные условия, кроме естественных, могут быть также следующими:

Ui(L(t) - 0,t) = Ui(L(t) + 0, t) = /1 (t);

U2(L(t) - 0,t) = U2(L(t) + 0,t) = /2(t);

U2,x(L(t) - 0, t) = U2,x(L(t) + 0, t) = /3(t),

где /i(t), /2(t), /3(t) —заданные функции.

5.4. Начальные условия. Выпишем части вариаций выражений (9),(11), соответствующих интервалу от t = 0 до t = ti:

/ d I \ ti fL(0) 0 .

m4dtU(L(t),t)^Ui|x=L(tJ - Uk,t(x, 0)6ufc|t=Q dx-

/•Lo rL(0)-0

- Uk,t(x, 0)6ufc . 0 dx +/ Uk,t(x, 0)6ufc. dx+

-/l(o)+o ло 1

f Lo

+ / Ufc,t(x, 0)6ufc . dx = 0. (34)

Jl( o)+o t=ti

Начальные условия обычно задаются при t = 0 в виде

Ui(x, 0) = ^i(x);

U2 (x, 0) = ^2(x);

Ui,t(x, 0) = ^3(x);

U2,t(x, 0) = ^4(x),

(35)

где ^i(x), ^>2(x), ^3(x), ^>4(x) — заданные функции. Вариации заданных функций равны нулю, что обеспечивает равенство нулю частей выражения (34) при t = 0.

Если решение задачи однозначно, то функции Ui(x,t) и U2(x,t) при t = ti однозначно определены и их вариации равны нулю. Это обеспечивает равенство нулю частей выражения (34) при t = ti.

Таким образом, модель, описывающая колебания объекта, изображенного на рис. 1, определяется системой дифференциальных уравнений (28); граничными условиями (29), (31), (32); соотношениями (33) и начальными условиями (35).

6. Линеаризация задачи. Предположим, что смещения Ui, U2 и их производные малы. При этом квадратичными членами, а также значениями £о и ul,x, U2,x по сравнению с единицей можно пренебречь. При линеаризации задача получается расщеплённой, то есть отсутствует взаимодействие между продольными и поперечными колебаниями.

6.1. Задача для продольных колебаний.

pUi,tt Eul,xx ^ui,xxt 0;

Ui(0,t)=0; ui (Lo, t) = 0;

392

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

d2

mid^u(L(t),t) - Sp(u*}t(L(t) + 0,t) - - 0,t))L'(t)-

- SE(u1,x(L(t) +0,t) - ul,x(L(t) - 0,t))-

- SMui,xt(L(t)+ 0,t) - ui,xt(L(t) - 0,t)) - Fi(t) =0;

u\ (L(t) + 0, t) = u\ (L(t) - 0, t); u[(x, 0) = ^i(x); ul,t(x, 0) = ^2(x).

6.2. Задача для поперечных колебаний.

pSu2,tt ES£qu2,xx + EIu2,xxxx + PIu2,xxxxt + ^u2,t + K0u2 f(x, t) = 0;

u2(0, t) = 0; u2,xx(0, t) = 0; u2(Lo, t) = 0; u2,x(Lo, t) = 0; d2

m2dt2u2(L(t),t) + EIu2,xxx(L(t) + 0, t) + jUIu2,xxxt(L(t) + 0, t)-

EIu2,xxx(L(t) 0, t) PIu2,xxxt(L(t) 0, t) +

+ K2u2(L(t),t) - F2(t) = 0;

u2(L(t) - 0, t) = u2(L(t) + 0, t); u2,x(L(t) - 0,t) = 0; u2,x(L(t) + 0,t) = 0; u2 (ж, 0) = ^э(ж); u2,t(x, 0) = ^4(ж).

Полученные линеаризованные задачи соответствуют классическим постановкам таких задач [1,5,12,19,20].

Заключение. Таким образом, в работе поставлены новые нелинейные математические модели продольно-поперечных колебаний объектов с движущимися границами, при этом соблюдается принцип однородности: в частном случае малых колебаний предложенные модели совпадают с классическими моделями, что свидетельствует о корректности полученных результатов. Решение конкретных краевых задач на основе полученных моделей представляет собой самостоятельную задачу и не является целью данного исследования. ORCIDs

Валерий Николаевич Анисимов: http://orcid.org/0000-0002-1346-167X Владислав Львович Литвинов: http://orcid.org/0000-0002-6108-803X

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Горошко О. А., Савин Г. Н. Введение в механику деформируемых одномерных тел переменной длины. Киев: Наукова думка, 1971. 270 с.

2. Zhu W. D., Chen Y. Theoretical and experimental investigation of elevator cable dynamics and control// J. Vib. Acoust., 2006. vol. 128, no. 1. pp. 66-78. doi: 10.1115/1.2128640.

3. Самарин Ю. П., Анисимов В. Н. Вынужденные поперечные колебания гибкого звена при разгоне// Изв. вузов. Машиностроение, 1986. №12. С. 17-21.

4. Boyle (Jr) J. M., Bhushan B. Vibration modeling of magnetic tape with vibro-impact of tape-guide contact // J. Sound Vibr., 2006. vol. 289, no. 3. pp. 632-655. doi: 10.1016/j. jsv.2005.02.033.

5. Лежнева А. А. Изгибные колебания балки переменной длины // Изв. АН СССР. МТТ, 1970. №1. С. 159-161.

393

Анисимов В. Н., Литвинов В. Л.

6. Ding Hu, Chen Li-Qun Galerkin methods for natural frequencies of high-speed axially moving beams// J. Sound Vibr., 2010. vol. 329, no. 17. pp. 3484-3494. doi: 10.1016/j. jsv.2010.03.005.

7. Guo Y., Yang S., Guo W. Analysis of dynamic characteristics of steel spring supported floating track bed // Zhendong Ceshi Yu Zhenduan = Journal of Vibration, Measurement and Diagnosis, 2006. vol. 26, no. 2. pp. 146-150 (In Chinese).

8. Lei X.-Y Effects of abrupt changes in track foundation stiffness on track vibration under moving loads // Zhendong Gongcheng Xuebao = Journal of Vibration Engineering, 2006. vol. 19, no. 2. pp. 195-199 (In Chinese).

9. Sahebkar S. M., Ghazavi M. R., Khadem S. E.,Ghayesh M. H. Nonlinear vibration analysis of an axially moving drillstring system with time dependent axial load and axial velocity in inclined well // Mechanism and Machine Theory, 2011. vol. 46, no. 5. pp. 743-760. doi: 10. 1016/j.mechmachtheory.2010.12.003.

10. Inacio O., Antunes J., Wright M. C. M. Computational modelling of string-body interaction for the violin family and simulation of wolf notes // J. Sound Vibr., 2008. vol. 310, no. 1-2. pp. 260-286. doi: 10.1016/j.jsv.2007.07.079.

11. Тихонов В. С., Абрамов А. А. Поперечные колебания гибкой нити переменной длины в потоке // Вестник МГУ. Сер. 1, 1993. №5. С. 45-48.

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

13. Анисимов В. Н., Литвинов В. Л. Исследование резонансных свойств механических объектов при помощи метода Канторовича-Галёркина // Вестн. Сам. гос. техн. унта. Сер. Физ.-мат. науки, 2009. №1(18). С. 149-158. doi: 10.14498/vsgtu658.

14. Анисимов В. Н., Литвинов В. Л., Корпен И. В. Об одном методе получения аналитического решения волнового уравнения, описывающего колебания систем с движущимися границами// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2012. №3(28). С. 145-151. doi: 10.14498/vsgtu1079.

15. Kotera T. Vibration of a string with time-varying length / Memoirs of the Faculty of Engineering, Kobe University, 24, 1978. pp. 45-54 (in Japanese).

16. Zhu W. D., Zheng N. A. Exact response of a translating string with arbitrarily varying length under general excitation / ASME 2007 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. vol. 1 (September 4-7, 2007), 21st Biennial Conference on Mechanical Vibration and Noise, Parts A, B, and C. Las Vegas, Nevada, USA, 2008. pp. 1995-2013. doi: 10.1115/detc2007-34590.

17. Zhu W. D., Zheng N. A. Exact response of a translating string with arbitrarily varying length under general excitation// J. Appl. Mech., 2008. vol. 75, no. 3, 031003. 14pp.. doi: 10.1115/ 1.2839903

18. Brake M. R., Wickert J. A. Frictional vibration transmission from a laterally moving surface to a traveling beam// J. Sound Vibr., 2008. vol. 310, no. 3. pp. 663-675. doi: 10.1016/j. jsv.2007.04.029.

19. Мышкис А. Д. Математика для технических вузов. СПб.: Лань, 2002. 640 с.

20. Анисимов В. Н., Литвинов В. Л., Корпен И. В. Постановка задачи о колебаниях балки с движущейся подпружиненной опорой // Вестн. Сам. гос. техн. ун-та. Сер. Техн. науки, 2013. №1(37). С. 93-98.

Поступила в редакцию 05/IX/2014; в окончательном варианте — 18/II/2015; принята в печать — 08/IV/2015.

394

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

Vestn. Samar. Gos. Techn. Un-ta. Ser. Fiz.-mat. nauki

[J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2015, vol. 19, no. 2, pp.382—397

ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1330

MSC: 35R37; 35G30, 35Q70

MATHEMATICAL MODELS OF NONLINEAR LONGITUDINAL-CROSS OSCILLATIONS OF OBJECT WITH MOVING BORDERS

V. N. Anisimov, V. L. Litvinov

Syzran’ Branch of Samara State Technical University,

45, Sovetskaya str., Syzran’, Samara region, 446001, Russian Federation.

Abstract

The nonlinear formulation of problems for describing longitudinal-cross oscillations of objects with moving borders is noted. These mathematical models consist of a system of two nonlinear partial differential equations with the higher time derivative of the second order and the fourth-order by the spatial variable. The nonlinear boundary conditions on moving boundary have a higher time derivative of the second order and the third-order by the spatial variable. The geometric nonlinearity, visco-elasticity, the flexural stiffness of the oscillating object and the elasticity of the substrate of object are taken into account. Boundary conditions in the case of energy exchange between the parts of the object on the left and right of the moving boundary are given. The moving boundary has got a joined mass. The elastic nature of borders joining is considered. The longitudinal-cross oscillations of objects with moving borders of high intensity can be described by the resulting differential model. The Hamilton’s variational principle is used in the formulation of the problem.

Keywords: longitudinal-cross oscillations, moving borders, boundary value problems, mathematical models, boundary conditions, nonlinear system of partial differential equations, variational principles.

doi: http://dx.doi.org/10.14498/vsgtu1330

ORCIDs

Valeriy N. Anisimov: http://orcid.org/0000-0002-1346-167X Vladislav L. Litvinov: http://orcid.org/0000-0002-6108-803X

© 2015 Samara State Technical University.

Please cite this article in press as:

Anisimov V. N., Litvinov V. L. Mathematical models of nonlinear longitudinal-cross oscillations of object with moving borders, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2015, vol. 19, no. 2, pp. 382-397. doi: 10.14498/vsgtu1330. (In Russian)

Authors Details:

Valeriy N. Anisimov (Cand. Phys. & Math. Sci.; anisimov170159@mail.ru; Corresponding Author), Head of Dept., Dept. of General-Theoretical Disciplines.

Vladislav L. Litvinov (vladlitvinov@rambler.ru), Teacher, Dept. of General-Theoretical Disciplines.

395

Анисимов В. Н., Литвинов В. Л.

REFERENCES

1. Goroshko O. A., Savin G. N. Vvedenie v mekhaniku deformiruemykh odnomernykh tel peremennoi dliny [Introduction to the Mechanics of Deformable One-Dimensional Bodies of Variable Length]. Kiev, Naukova dumka, 1971, 270 pp. (In Russian)

2. Zhu W. D., Chen Y. Theoretical and experimental investigation of elevator cable dynamics and control, J. Vib. Acoust., 2006, vol. 128, no. 1, pp. 66-78. doi: 10.1115/1.2128640.

3. Samarin Yu. P., Anisimov V. N. Forced transverse vibrations of the flexible link at dispersal, Izv. vuzov. Mashinostroenie, 1986, no. 12, pp. 17-21 (In Russian).

4. Boyle (Jr) J. M., Bhushan B. Vibration modeling of magnetic tape with vibro-impact of tape-guide contact, J. Sound Vibr., 2006, vol. 289, no. 3, pp. 632-655. doi: 10.1016/ j.jsv.2005.02.033.

5. Lezhneva A. A. Bending vibration of beam of variable length, Izv. Akad. Nauk USSR. MTT, 1970, no. 1, pp. 159-161 (In Russian).

6. Ding Hu, Chen Li-Qun Galerkin methods for natural frequencies of high-speed axially moving beams, J. Sound Vibr., 2010, vol. 329, no. 17, pp. 3484-3494. doi: 10.1016/j.jsv. 2010.03.005.

7. Guo Y., Yang S., Guo W. Analysis of dynamic characteristics of steel spring supported floating track bed, Zhendong Ceshi Yu Zhenduan = Journal of Vibration, Measurement and Diagnosis, 2006, vol. 26, no. 2, pp. 146-150 (In Chinese).

8. Lei X.-Y Effects of abrupt changes in track foundation stiffness on track vibration under moving loads, Zhendong Gongcheng Xuebao = Journal of Vibration Engineering, 2006, vol. 19, no. 2, pp. 195-199 (In Chinese).

9. Sahebkar S. M., Ghazavi M. R., Khadem S. E.,Ghayesh M. H. Nonlinear vibration analysis of an axially moving drillstring system with time dependent axial load and axial velocity in inclined well, Mechanism, and Machine Theory, 2011, vol. 46, no. 5, pp. 743-760. doi: 10. 1016/j.mechmachtheory.2010.12.003.

10. Inacio O., Antunes J., Wright M. C. M. Computational modelling of string-body interaction for the violin family and simulation of wolf notes, J. Sound Vibr., 2008, vol. 310, no. 1-2, pp. 260-286. doi: 10.1016/j.jsv.2007.07.079.

11. Tikhonov V. S., Abramov A. A. Transverse Vibrations of a Flexible String with Time-Varying Length in Flow, Vest. Mosk. Univ. Ser 1. Matematika, Mekhanika, 1993, no. 5, pp. 45-48 (In Russian).

12. Vesnitskii A. I. Volny v sistemakh s dvizhushchimisia granitsami [Waves in systems with moving boundaries and loads]. Moscow, Fizmatlit, 2001, 320 pp. (In Russian)

13. 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 [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2009, no. 1(18), pp. 149-158 (In Russian). doi: 10.14498/vsgtu658.

14. Anisimov V. N., Litvinov V. L., Korpen I. V. On a method of analytical solution of wave equation describing the oscillations sistem with moving boundaries, Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2012, no. 3(28), pp. 145-151 (In Russian). doi: 10.14498/vsgtu1079.

15. Kotera T. Vibration of a string with time-varying length, Memoirs of the Faculty of Engineering, Kobe University, 24, 1978, pp. 45-54 (in Japanese).

16. Zhu W. D., Zheng N. A. Exact response of a translating string with arbitrarily varying length under general excitation, ASME 2007 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, vol. 1 (September 4-7, 2007), 21st Biennial Conference on Mechanical Vibration and Noise, Parts A, B, and C. Las Vegas, Nevada, USA, 2008, pp. 1995-2013. doi: 10.1115/detc2007-34590.

17. Zhu W. D., Zheng N. A. Exact response of a translating string with arbitrarily varying length under general excitation, J. Appl. Mech., 2008, vol. 75, no. 3, 031003, 14 pp.. doi: 10.1115/ 1.2839903

396

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

18. Brake M. R., Wickert J. A. Frictional vibration transmission from a laterally moving surface to a traveling beam, J. Sound Vibr., 2008, vol. 310, no. 3, pp. 663-675. doi: 10.1016/j.jsv. 2007.04.029.

19. Myshkis A. D. Matematika dlia tekhnicheskikh vuzov [Mathematics for Technical Colleges]. St. Petersburg, Lan’, 2002, 640 pp. (In Russian)

20. Anisimov V. N., Litvinov V. L., Korpen I. V. The formulation of the problem of the beam fluctuations with moving spring-loaded support, Vestn. Samar. Gos. Tekhn. Univ. Ser. Tekhn. Nauki, 2013, no. 1(37), pp. 93-98 (In Russian).

Received 05/IX/2014;

received in revised form 18/II/2015;

accepted 08/IV/2015.

397