Заочные электронные конференции
Логин   Пароль  
Регистрация Забыли пароль?
 
     
О МЕТОДАХ РАСЧЕТА ВИБРОУДАРНЫХ СИСТЕМ С БОЛЬШИМ ЧИСЛОМ УДАРНЫХ ПАР
Крупенин В.Л.


Для чтения PDF необходима программа Adobe Reader
GET ADOBE READER

УДК 534.1

О МЕТОДАХ РАСЧЕТА ВИБРОУДАРНЫХ СИСТЕМ

С БОЛЬШИМ ЧИСЛОМ УДАРНЫХ ПАР

© Виталий Львович Крупенин

Учреждение Российской академии наук институт машиноведения РАН

им. А.А. Благонравова

Аннотация. Дается представление периодических виброударных процессов в системах с произвольным конечным числом сосредоточенных ударных пар. Это представление, названное 2n-параметрическим, используется для расчета виброударного режима движения балки, колеблющейся вблизи двустороннего решетчатого ограничителя хода. Исследованы некоторые динамические эффекты, свойственные изучаемым системам. Проведен расчет реальной конструкции.

2n-параметрическое представление. Представим периодический виброударный процесс в склерономной стационарной механической системе общего вида, содержащей n сосредоточенных ударных пар. При построении в несколько модифицированном виде используем методы частотно-временного анализа виброударных процессов [1—3], основанные на записи искомых законов движения через периодические функции Грина взаимодействующих линейных механических подсистем.

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

Рис. 1

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

(1)

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

Отметим также, что совершенно аналогично может рассматриваться случай, когда принадлежащие изучаемой системе сосредоточенные тела соударяются между собой. В этом случае соотношение типа (1) записывается для относительных координат ), где — абсолютные координаты контактирующих тел. Соответственно видоизменяются и условия удара [1—3].

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

(2)

( — неизвестные импульс и фаза удара в - й ударной паре, ()и периодическую обобщенную функцию

(3)

где - функция Дирака, которую иногда называют «последовательность -функций»; .

Известно [2], что в условиях периодичности соотношения (1) оказываются эквивалентны интегральным представлениям

(4)

где

(5)

-периодическая несимметричная функция Грина, отвечающая оператору [2].

С учетом (2) и (3) получаем

(6)

где — движение в отсутствие ударов (установившееся решение линейной задачи).

Данное представление содержит 2пнеизвестных параметров движения и называется 2n-параметрическим представлением периодического виброударного процесса. Для нахождения неизвестных параметров имеется 2псоотношений совместности, следующих из условий удара

(7)

При этом необходимо для всех и для всех k. Систему (7) можно расписать подробнее:

(8)

Первая группа уравнений — это система n линейных алгебраических уравнений относительно импульсов ударов в каждой из ударных пар. Вторую группу получим из второго соотношения (7). Для имеем

(9)

Данная группа также является системой линейных алгебраических уравнений относительно импульсов В то же время системы (8) и (9) нелинейны относительно фаз . В конкретных случаях, используя свойства ПФГ [2], уравнения (8) и (9) могут быть упрощены.

В случаях, когда ударные пары симметричны или ищутся субгармонические или комбинационные режимы, вид соотношений (6), (8), (9) остается практически тем же; меняются только конкретные значения величин и . Например, при рассмотрении симметричных ударных пар и отыскании симметричных режимов все соотношения строятся, исходя из симметричных ПФГ [2 ]. Вместо (5) здесь необходима замена и

(10)

После решения системы (8), (9) и определения представления (6) можно найти перемещение любой точки у

(11)

где смысл обозначений очевиден.

Заметим, что возможен случай, когда процесс сопровождается ударами лишь в ударных парах, так что ударных импульсов равны нулю. Если и колебания в единственной «действующей» ударной паре имеют амплитуду, существенно превосходящую амплитуды всех прочих точек системы, то это свидетельствует об установлении локализованного режима, экспериментально описанного, например, в [4, 5].

Хлопок в цепочке упругосвязанных масс.Хлопком называется режим движения, при котором все массы двигаются синфазно и одновременно ударяются об ограничитель. При хлопке виброударный процесс наиболее интенсивен. Рассмотрим модель системы с периодической структурой (рис. 1, б). Не ограничивая общности рассуждений, для простоты будем считать . Предположим, что система консервативна. Тогда уравнения движения для и граничные условия имеют вид

(12)

где и удар предполагается абсолютно упругим.

Пусть в системе осуществляется хлопок, в котором участвуют все частицы. Тогда 2n-параметрическое представление (6) и определяющая система (8), (9) при принимает вид

(13)

(14)

где удар совмещен с началом отсчета времени.

Учитывая соотношения [2 ] и решения [6 ], при найдем

(15)

где собственные частоты линейной цепочки .

Для получения периодической функции Грина (ПФГ) вне отрезка функцию (15) надо продолжить по периодичности, исходя из условия для всех .

В силу ортогональности собственных форм второе соотношение (14); при всех k является тождеством (решение автономных систем определяются с точностью до произвольных фаз). Вектор импульсов легко отыскивается из системы линейных уравнений, следующих из первого соотношения (14) при Ранее подобное решение было выписано в [1 ].

Рассмотрим трехпролетную цепочку с двумя массами (n = 2). Из (14) следует

(16)

Так как и , то . Для собственных частот находим . Соотношения (13) — (15) дают

что в точности отвечает движению обычного «ударного осциллятора» [1,2]. Из условия следует, что хлопки возможны в частотном диапазоне . Из (12) для искомых законов движения получаем

т.е. движение каждого тела аналогично движению «ударного осциллятора». Подобный вывод справедлив при любом n.

Система с решетчатым ограничителем. Рассмотрим систему, представленную на рис. 2. Шарнирно закрепленный стержень может при совершении

Рис. 2

колебаний в плоскости чертежа соударяться с nплоскими протяженными ограничителями ширины . Параметры системы: длина стержня - растягивающее усилие жесткость и погонная масса стержня Kи установочные зазоры координаты начал ограничителей - концов - Будем рассматривать симметричную систему

Предполагая, что малая диссипация учитывается только при ударах и мал уровень возбуждения, запишем уравнение движения в предположениях, что стержень моделируется балкой Бернулли

(17)

Здесь граничные условия отвечают свободному опиранию; - малый параметр - периодическое возбуждение, предполагаемое симметричным; и - упругая и малая диссипативные составляющие плотности силы ударного воздействия. Появление малого параметра в правой части говорит о том, что исследоваться будут колебания, близкие к собственным.

Плотность ударных сил, действующих при , записывается через плотность ударного импульса и распределение фазы [3, 7 ]

(18)

где — функция Хевисайда; симметричная Т— периодическая последова­тельность -функций Дирака [2];

Для искомого поля перемещений имеет место интегральное представление [3, 7], являющееся аналогом 2n-параметрического представления (6)

(19)

где — решение соответствующей линейной задачи; — плотность импульса и распределения фазы для k- пары в решетке.

Симметричная ПФГ балки при имеет вид [1, 2 ]

(20)

и отвечает оператору, действующему из точки z в точку х. Собственные частоты линейной системы . Вне отрезка фун­кция (20) продолжается на всю ось по периодичности, исходя из условий симметрии .

Воспользуемся малостью величин : волновые процессы в самом ударном элементе оказывают на режим слабое влияние и вместо (19) можно перейти к 2n-параметрическому представлению вида

(21)

где фигурирует полный импульс в k-паре

В формулу (21) входит усреднение ПФГ

(22)

Подставляя (20) в (22) , после интегрирования имеем

(23)

Представление (21) вполне аналогично (б). Оно предполагает, что масса приведенного точечного контактирующего тела , а импульс , где — коэффициент восстановления в k-й паре. Используя (21), можно получить определяющую систему уравнений, подобную (8), (9). Для получения более просто интерпретируемых решений найдем в консервативном случае одночастотный периодический режим. Затем, предположив, что на резонансе режим движения приближенно отвечает кон­сервативному решению, найдем условие существования резонансных режимов, исходя из принципа энергетического баланса, в соответствии с которым на периодическом движении балансируются работы неконсервативных сил — воз -буждения и диссипации [1 — 3].

Положив , для свободных колебаний некоторого периода Т0,ограничиваясь рассмотрением синфазных хлопков и совместив одновременные удары о правые ограничители с началом отсчета времени, найдем:

(24)

(25)

Можно показать, что с принятой точностью соотношения (25) удовлетворяются, а решения линейной алгебраической системы (24) полностью определяют представление

(26)

Следует отметить, что при получении решений консервативной задачи период усредненной ПФГ Т0не равен периоду Т. В консервативном случае период Т0есть функция величин . Далее, при изучении одночастотных резонансных режимов будем предполагать, что Т = Т0.

Расчет трехпролетной трубки. Рассмотрим трехпролетную трубку (n =2), которую корректно моделировать при посредстве балки Бернулли. Положим . Разрешая систему (24) и обозначая частоту свободных колебаний , с учетом (23) находим , где

(27)

Это вместе с (26) и (23) полностью описывает режимы в консервативной системе. Приравнивая нулю фигурную скобку (27) , можно отыскать частоты, отвечающие бесконечным импульсам.

Проведем расчеты при следующих геометрических и физических ? характеристиках трубки [7]: наружный диаметр , внутренний диаметр , модуль Юнга (отсюда изгибная жесткость , длина трубки

Частотные зоны, в которых формально могут существовать режимы движения типа хлопков, определяются условием . Для нахождения границ этих зон необходимо определить частоты, при которых . Расчеты показывают, что одна такая зона

Рис.3

существует справа от первой собственной частоты линейной системы и занимает отрезок . Расчеты были выполнены численно. Учитывали 20 членов ряда (27). Сходимость контролировали численно. Наряду с этой зоной были найдены другие, оказавшиеся весьма узкими, и их рассмотрение не имеет практического смысла. Данное замечание находится в согласии с экспериментом [4].

Ограничиваясь экспериментально наблюдаемыми хлопками — режимами, ус­танавливающимися справа от первой собственной частоты линейной системы, можно получить зависимости, показанные на рис. 3 (кривая 1 отвечает случаю ; 20,85; 31).

Пусть в (17) . Рассмотрим вынужденные резонансные колебания, для осуществления которых необходимо, чтобы частота попала в интервал собственных частот . Будем описывать движения при посредстве (26), (27), считая и . Совместив удар с началом отсчета времени, считаем фазу неизвестной. Для ее нахождения воспользуемся принципом энергетического баланса, который в случае пударных пар принимает вид

(28)

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

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

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

Работа выполнена при финансовой поддержке Российского фонда фундамен-льных исследований (проект 10-08-00500).

Литература
  1. Бабицкий В. И. Теория виброударных систем. М.: Наука, 1978. 352 с.

  2. Бабицкий В. И., Крупе нин В. Л. Колебания в сильно нелинейных системах. М.: Наука, 1985.

  3. 320 с. Веприк А. М., Вознюк 77. Д., Крупенин В. Л. и др. Широкополосные виброударные генераторымеханических колебаний. Л.: Машиностроение, 1987. 79 с. Вабицкий В. И., Веприк А. А/., Крупенин В. Л. Виброударные эффекты при ограниченных колебаниях одномерной цепочки упруго соединенных бусинок//ДАН СССР. 1989. Т. 300. № 3. С. 562—566.

  4. Крупенин В. Л. Трансформация форм колебаний струны, взаимодействующей с двумя протяженными преградами//ДАН СССР. 1990. Т. 313. № 6. С. 1390—1394.

  5. Коткин Г. Л., Сербо В. Г. Сборник задач по классической механике. М.: Наука, 1977. 320 с.

  6. Крупенин В. Л. К теории виброударных систем с распределенными ударными элементами //Изв.

  7. АН СССР. МТТ. 1986. № 1. С. 25—32. Парогенераторы АЭС вертикальные. Расчет вибраций теплообменных труб. РТМ-108. Л.: Изд-во

НПО ЦИТИ, 1986. 101 с.

Библиографическая ссылка

Крупенин В.Л. О МЕТОДАХ РАСЧЕТА ВИБРОУДАРНЫХ СИСТЕМ С БОЛЬШИМ ЧИСЛОМ УДАРНЫХ ПАР // Научный электронный архив.
URL: http://econf.rae.ru/article/6375 (дата обращения: 18.10.2019).



Сертификат Получить сертификат