CC BY
4Электронный научный журнал "Математическое моделирование, компьютерный и натурный эксперимент в естественных науках" http://mathmod.esrae.ru/ URL статьи: mathmod.esrae.ru/37-134 Ссылка для цитирования этой статьи:
Сопенко А.А., Петрунин И.А. Численное моделирование колебаний пологой нелинейной оболочки под действием подвижной нагрузки // Математическое моделирование, компьютерный и натурный эксперимент в естественных науках. 2022. №1_
УДК 539.3; 541.1 DOI: 10.24412/2541-9269-10.24412/2541-9269-2022-1-22-33
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ КОЛЕБАНИЙ ПОЛОГОЙ НЕЛИНЕЙНОЙ ОБОЛОЧКИ ПОД ДЕЙСТВИЕМ ПОДВИЖНОЙ
НАГРУЗКИ
Сопенко А. А.1, Петрунин И.А.2
1 Саратовский государственный технический университет имени Гагарина Ю.А., Россия, Саратов, saasar@mail.ru
2 Саратовский государственный технический университет имени Гагарина Ю.А., Россия, Саратов, ivan.petrunin2000@yandex.ru
NUMERICAL SIMULATION OF OSCILLATIONS OF A SHALLOW NONLINEAR SHELL UNDER THE ACTION OF A MOVING LOAD
Sopenko A.A.1, Petrunin I.A.2
1 Yuri Gagarin State Technical University of Saratov,
Russia, Saratov, saasar@mail.ru
2 Yuri Gagarin State Technical University of Saratov,
Russia, Saratov, ivan.petrunin2000@yandex.ru
Аннотация. В работе рассматриваются сложные колебания геометрически нелинейной пологой оболочки под действием равномерно распределённой по некоторой части плана оболочки нагрузки постоянной интенсивности, движущейся с постоянной скоростью вдоль одной из координатных осей. Рассматривается смена форм колебаний оболочки в зависимости от скорости движения и величины приложенной нагрузки.
Ключевые слова: сложные колебания, пологие оболочки, подвижная нагрузка.
Abstract. The paper considers complex oscillations of a geometrically nonlinear shallow shell under the action of a load of constant intensity uniformly distributed over some part of the shell plan, moving at a constant speed along one of the coordinate axes. The change of shell oscillation forms is considered depending on the speed of movement and the magnitude of the applied load.
Keywords: complex oscillations, shallow shells, moving load.
Задачи изучения колебаний балок, пластин, оболочек под действием подвижных нагрузок ведут свою историю еще с XIX века. В исследованиях академика А.Н. Крылова, в частности, в начале XX века проводится наиболее полное на то время исследование реакции балки на подвижную нагрузку. Задачи о поведении конструкций под действием подвижных нагрузок возникают при изучении монорельсовых конструкций [1], движений колеса [2] и огромного количества других задач.
Укажем на работу [3], в которой, в частности, содержится исторический и научный обзоры данного вопроса, приводятся математические модели различных задач, связанных с подвижными нагрузками, рассматриваются приближенные аналитические и численные методы решения задач, связанных с балками, арочными системами, цилиндрическими и другими оболочками.
Дальнейшее развитие вычислительное техники расширяет возможности численного решения задач. Данная работа посвящена изучению колебаний прямоугольной в плане сферической геометрически нелинейной пологой оболочки под действием равномерно распределённой на части плана оболочки нагрузки интенсивности q. Эта часть плана представляет собой полосу размером 0 < х1 < а, где а - размер оболочки вдоль ^. Относительно оси х2 ширина этой полосы составляет Ь /8, где Ь - размер оболочки относительно х2. Эта полоса нагрузки движется вдоль х2 с постоянной скоростью V. Как только полоса нагружения уходит с плана оболочки на линии х2 = Ь, на линии х2 = 0 начинает свое движение новая такая же полоса, обеспечивая периодическое нагружение оболочки.
При этом не учитываются силы инерции подвижной нагрузки. Ширина нагруженной части плана оболочки в 1/8 размера относительно оси х2 в данной задаче не менялась. Рассматривались колебания оболочки в зависимости от интенсивности нагружения q и скорости движения нагрузки V.
Система уравнений, описывающих указанную задачу, широко известна в механике тонкостенных конструкций. В рамках модели Кирхгофа - Лява в смешанной форме в безразмерном виде она имеет следующий вид [4]:
-V4 я - Ця, Р) - V? Р - q + к (>+ ея) = 0
12(1 - V?) к
(1)
V4 Р+V? я+1 Ь(я, я) = 0 к 2
Использованы обозначения: Р и я - функция напряжений и прогиб в направлении, перпендикулярном срединной поверхности соответственно. я считается положительным, если направлен к центру кривизны. Е и V - модуль Юнга и коэффициент Пуассона, р- плотность, е - коэффициент демпфирования. Значение к приведено в (2).
Использована декартова система координат х1, х2, х3, прогибы w откладываются вдоль оси х3, размеры оболочки вдоль указанных осей обозначатся а, Ь, И. Безразмерные переменные для (1) вводились следующим образом, здесь черточки стоят над безразмерными величинами, во всех формулах, кроме (2), черточки для удобства и краткости опущены
X
XI = —, Х2 = а
х„ -
Ь ' Л =
Хз =
хз
И
— w w = —. И
г = F
ЕИ
а
г =
г-а
И2
(2)
27 2 2 - а Ь — 7 а а = а—г, к = к— 4 4 ЕИ4 1 1 И
Т- . Ь2 -
к2 = к2
И2
И
8 = 8-
а
к =
а Ь ра
ЕИ6
В (2) к, К - безразмерные кривизны оболочки, а - коэффициент температуропроводности. Под обозначением в (1) понимается
д\-)д\о) д\-)д\о) =--—---1---——
2-
3х15х2 3х15x2
(3)
Система (1) неоднократно приводилась во многих работах, поэтому остальные использованные обозначения более подробно описаны в [5].
В начальный момент времени оболочка находится в покое. По всему контуру оболочка считается шарнирно опертой на гибкие, нерастяжимые в касательной плоскости рёбра. Так, для краёв х1 = 0; а использованы граничные условия
w =—Т = 0, 822
0, 5822
58
12
"2 5х^ 5х 2
0.
(4)
Алгоритм численного решения (1) также подробно изложен в [5]. Укажем только, что использован метод конечных разностей для дискретизации производных по х, х2, х3, в результате чего первое уравнение из (1) сводится к системе обыкновенных дифференциальных уравнений, для решения которых достаточно использовать метод Рунге - Кутты 4 порядка точности. Второе уравнение (1) сводится к системе линейных уравнений, которая решается на каждом шаге по времени методом Гаусса.
Использование для расчётов подвижной нагрузки интенсивности а ,
движущейся со скоростью V относительно оси х2, приводит к тому, что
решается задача, не симметричная относительно этой оси. Если обычно задача решается на 1 / 4 плана оболочки в предложении, что решение симметрично
относительно осей симметрии оболочки, то здесь приходится решать задачу на 1 / 2 плана, что приводит к удвоению затрат машинного времени на решение. Была рассмотрена оболочка с параметрами а = Ь = 0,1м, а / И = 100,
Е = 69ГПа, V = 0,3, р = 2800кг / м3, а = 6,4 • 10-5м2 / сек. Физические
параметры материала соответствуют сплаву АМц. к 1 = к 2 = 24.
Отметим также, что безразмерная собственная частота колебаний данной оболочки, найденная численно [5], составила со = 30,22.
Далее все полученные результаты приводятся в безразмерном виде (кроме скорости V) для центральной точки плана оболочки.
V = 940м / сек
Таблица 1
Продолжение таблицы 1
Максимально возможная скорость V при численном решении задачи фактически определяется шагом по времени, используемым решения системы обыкновенных дифференциальных уравнений из (1). В табл. 1 приводятся результаты, имеющие скорее научный интерес, чем практический. Там показана зависимость прогиба w от времени (далее будем называть это сигнал) для указанной выше точки плана оболочки для V = 940м / сек = 3384км / ч. Колебания оболочки после завершения начального периода возмущений прекращаются, прогибы становятся постоянными. Видно, что между значениями q = 416,5 и 416,6 происходит прохлопывание оболочки.
Если уменьшить значение V, переходя к более «земным» скоростям, то после начального периода мы будем наблюдать колебания в расширяющемся диапазоне значений. При этом, если не рассматривать начальный возмущённый период, при малых значениях q всегда будут наблюдаться гармонические колебания на различных частотах - в зависимости от значения V. Поэтому нам показалось интересным рассмотреть значения V, вызывающие колебания оболочки на её собственной частоте колебаний, а также при близких к этому значениях V.
Далее приводятся результаты, полученные для значений следующих значений V: 15,78м / сек = 56,82км / ч; 11,69 м / сек = 42,09км / ч; 10,52м / сек = 37,88км / ч
V = 15,78м / сек
Таблица 2
Продолжение таблицы 2
Сигнал
о
5
Спектр мощности
;_{_ч V____
о
1Г)
Спектр мощности
В табл. 2 приводятся полученные значения сигналов w(t) для нескольких значений а при значении V = 15,78м / сек. Даже если рассматривать только установившиеся колебания оболочки, не удаётся точно указать значение а, при которых происходит прохлопывание оболочки. Гармонические колебания имеют место при малых значениях а . С ростом интенсивности нагружения колебания очень быстро переходят в квазипериодические, а затем в хаотические, что показано на последних графиках табл. 2. Следует отметить, что максимальный пик частоты возбуждения, примерно равный 15 при а = 50 =,
исчезает, при а = 150 максимум достигается при значении, примерно равном
21. В табл. 3 приведены результаты решения поставленной задачи для значения V = 10,52м / сек.
V = 10,52м / сек_Таблица 3
Сигнал
0 5 10 15 20
I
3.5
О 5 10 15 20
I
Продолжение таблицы 3
о in
Сигнал
Спектр мощности
о in
Спектр мощности
Отметим, что здесь также не удаётся указать конкретное значение а или достаточно малый диапазон изменения а , при котором происходит резкое прохлопывание оболочки. С увеличением значения а прогибы увеличиваются достаточно монотонно. Спектр мощности, приведённый в табл. 3, в отличии от результатов приведённых в табл. 2, показывает наличие постоянного пика возбуждающей частоты, на которой начинаются колебания оболочки при малых а .
V = 11,69 м / сек
Таблица 4
Сигнал
Фазовый портрет
Спектр мощности
о
о 10
А ' А ' А ' I \ I \ I \ \ I ! I ! I Л ' Л \ 1 I1 ! 1
! I I 1 1 1 | ' 1 ! 1 11 / 1/ и \/ 1 1 1 1 / . V/
50 50.2 50.4
50.8 51
О О
сч
т сч
Продолжение таблицы 4
Сигнал
Фазовый портрет
Спектр мощности
.2
.3 4 2
|Г \1ж и И 1 1 1 1 ! Г, 1 Л м 1 1 1 ! 1 I / и
| 1 11 и 1/ 1 М! II ! и 1 V | * || V И ! I 11 1 || V ! ' || 1 || ' У V
50 50.2 50.4 50.6 50.8 51
О
3 3
А 11 1 ' г . Л |1 , I , А А
Л п Л1 1 1
' 1 1 '' 1 ] 1/ 1 ' / ! '
1 1 и. У у V V
50 50.2 50.4 50.6 50.8 51
О 00
3
И, наконец, в табл. 4 приведены результаты, полученные для нагрузки, движущейся со скоростью V = 11,69м / сек и вызывающей колебания оболочки на наиболее опасной частоте, близкой к частоте собственных колебаний конструкции. Последняя была найдена в [5] и оказалась равна <» = 30,22. При движении ц с указанной скоростью колебания начинались на частоте со = 29,9.
Колебания с частотой, близкой к собственной частоте колебаний оболочки, оказались самыми опасными для конструкции в смысле достижения максимально амплитуды колебаний при сравнении с нагрузками той же интенсивности, но перемещающихся с другими скоростями.
Сценарий смены форм колебаний заключался в переходе от гармонических колебаний к квазипериодическим и далее к хаотическим без
наблюдения эффектов перемещаемости форм колебаний. В указанном случае оказалось возможным определить нагрузку, при которой осуществляется прохлопывание оболочки. Это происходит между значениями q = 243,l и q = 243,2. После прохлопывания оболочки с ростом q вскорости наступает безвозвратный переход к хаотическим колебаниям.
Таким образом, расчёт поведения конструкций при подвижной нагрузке позволяет выявить некоторые необычные эффекты динамического поведения оболочки, определить наиболее опасные для конструкций режимы движения нагрузки.
Литература
1. Чиркин В.В., Петренко О.С., Михайлов А.С., Галонен Ю.М. Пассажирские монорельсовые дороги. М.: Машиностроение, 1969. 240 с.
2. Zeng-Xin Yu, Hui-Feng Tan, Xing-Wen Du, Li Sun. A simple analysis method for contact deformation of rolling tire // Vehicle System Dynamics. 2001. Vol. 36. P. 435-443.
3. Коноплев Ю.Г., Якушев Р. С. Лекции по динамике сооружений с подвижными нагрузками. Казань: Отечество, 2003. 208 с.
4. Вольмир А.С. Устойчивость деформируемых систем. М.: Наука, 1967. 984 с.
5. Сопенко А.А., Майорова О.А., Черепанов М.Д. Сложные колебания геометрически и физически нелинейных пологих оболочек // Математическое моделирование, компьютерный и натурный эксперимент в естественных науках. 2016. № 3. URL: mathmod.esrae.ru/3-16.
