CC BY
35
5. Предложен метод построения частотных
динамических диаграмм для оценки возможности
изменения динамических свойств систем при изменениях коэффициента связности внешних сил а.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Rocard Y. General Dynamics of Vibrations. Frederick Ungar Publishing Co. New-York. 1960. pp. 245.
2. Banakh L., Kempner M. Vibrations of Mechanical Systems with Regular Structure. Berlin; Heidelberg: Springer, 2010. pp. 262.
3. Karnovsky I.A., Lebed E. Theory of Vibration Protection. Springer International Publishing, Switzerland, 2016. 708 p.
4. Елисеев C.B., Артюнин А.И. Прикладная теория колебаний в задачах динамики линейных механических систем. Новосибирск : Наука, 2016. 459 с.
5. De Silva, C. W. Vibration. Fundamentals and Practice / Clarence W. de Silva. Boca Raton, London, New York, Washington, D.C.: CRC Press, 2000. 957 p.
6. Дружинский И.А. Механические цепи. М. : Машиностроение, 1977. 238с.
7. Ленк А. Электромеханические системы. Системы с сосредоточенными параметрами. М. : Мир, 1978. 283 с.
8. Елисеев С.В., Резник Ю.И., Хоменко А.П. Ме-хатронные подходы в динамике механических колебательных систем. Новосибирск : Наука, 2011. 384 с.
9. Елисеев, С.В., Хоменко А.П. Динамическое гашение колебаний: концепция обратной связи и структурные методы математического моделирования. Новосибирск : Наука, 2014. 357 с._
тшт
10.Генкин, М.Д., Рябой В.М. Упруго-инерционные виброизолирующие системы. Предельные возможности, оптимальные структуры. М. : Наука, 1988. 191 с.
11.Eliseev S.V., Lukyanov A.V., Reznik Yu.N., Khomenko A.P. Dynamics of mechanical system with additional ties. Irkutsk, 2006. 315 p.
12.Елисеев C.B., Ермошенко Ю.В. Сочленения звеньев в динамике механических колебательных систем. Иркутск : ИрГУПС. 2012. 156 с.
13.Белокобыльский, C.B., Елисеев C.B., Ситов И.С. Динамика механических систем. Рычажные и инерционно-упругие связи. СПб. : Политехника, 2013. 319 с.
14.Елисеев C.B., Кинаш Н.Ж., Кашуба В.Б. Рычажные связи механических колебательных систем // Вестник Всерос науч.-исслед. и проект.-конструктор. ин-та электровозостроения. 2016. № 1 (69). С. 112-126.
15.Хоменко А.П., Елисеев C.B. Квазиэлементы в механических колебательных системах. особенности систем при исключении переменных динамического состояния // Современные технологии. Системный анализ. Моделирование. 2013. № 2 (38). С. 8-17.
16.0 связях между координатами движения в механических колебательных системах с рычажными устройствами / C.B. Белокобыльский и др. // Системы. Методы. Технологии. 2015. № 2(26). С. 7-13.
17.Соотношения координат движения элементов механических колебательных систем как форма проявления рычажных связей / C.B. Белокобыльский и др. // Системы. Методы. Технологии. 2015. № 3 (27). С. 7-14.
УДК 534.1 Кошелее Александр Викторович,
к. т. н., научный сотрудник, ПАО Арзамасское научно-производственное предприятие «Темп-Аеиа», тел. 8(83147)-7-83-69, 8(910)895-33-91, e-mail: [email protected]
Денцое Николай Николаевич, к. т. н., доцент,
Нижегородский государственный технический университет им. Р.Е. Алексеева, тел. 8(904)399-66-08, 8(831)-436-03-01, e-mail: [email protected]
ОПРЕДЕЛЕНИЕ ОБЛАСТЕЙ ДИНАМИЧЕСКОЙ НЕУСТОЙЧИВОСТИ В ПАРАМЕТРИЧЕСКИХ СИСТЕМАХ С МАЛЫМ КОЭФФИЦИЕНТОМ ВОЗБУЖДЕНИЯ
A. V. Koshelev, N. N. Dentsov
DYNAMIC INSTABILITY DOMAINS DETERMINATION IN PARAMETRIC SYSTEMS WITH LOW EXCITATION COEFFICIENT
Аннотация. Решается задача повышения эффективности и стабильности рабочего режима колебаний вибрационных машин. Рассматривается механическая система инерционного типа, в которой реализуются режимы многократного комбинационного параметрического резонанса. Показано, что для существования колебаний необходимо выполнение порогового условия возбуждения, которое зависит от трения в колебательной системе, коэффициентов возбуждения и собственной ча -
Механика
стоты колебаний маятников. Кроме того, использование упругой изотропной подвески для возбуждения круговых колебаний рабочего органа снижает порог возбуждения в два раза. Уменьшение коэффициента собственной частоты маятников в два раза снижает порог возбуждения в три раза.
Показано, что в механической системе реализуется синергетический эффект за счет взаимного стимулирования парциальных подсистем. Этот эффект достигается неустойчивостью положения равновесия динамической системы при выполнении условия параметрического резонанса. В зоне параметрического резонанса положение равновесия сохраняется, но оно неустойчиво. В граничной точке (точке бифуркации) положение равновесия теряет устойчивость, что приводит к самовозбуждению параметрических резонансных колебаний.
Приводятся результаты моделирования областей неустойчивости и амлитудно-частотных характеристик вибрационной машины с различным коэффициентом возбуждения и величин трения. Показано, что при малом коэффициенте возбуждения и увеличенном трении технологической нагрузки происходит расширение резонансных зон за счет малого трения на маятниках параметрического роторно-маятникового возбудителя.
Результаты исследований позволяют выработать конструкторско-технологические мероприятия для разработки и создания вибрационных машин технологического назначения со стабильным резонансным режимом работы.
Ключевые слова: резонансная зона, маятник, область неустойчивости, демпфирование, параметрический резонанс.
Abstract. The problem of improving the efficiency and stability of the operating mode of vibration machine vibrations. Is considered. The mechanical system of inertia type, with the combination of multiple modes of parametric resonance, is analgen. It is shown that the existence of vibrations is possible with to satisfy the threshold condition of excitation, which depends on the friction in the oscillating system, the driving factors and the natural frequency of oscillation of the pendulum. In addition, the use of an isotropic elastic suspension for driving the circular working body vibrations reduces the excitation threshold two times, and a reduction ratio of own pendulum frequency halves the excitation threshold three times.
It is shown that a mechanical system realizes a synergistic effect due to mutual stimulation of partial subsystems. This effect is achieved by the instability of the equilibrium position of a dynamical system when the condition of parametric resonance is satisfied. In the zone of parametric resonance, the position of the equilibrium remains, but it is unstable. At the boundary point (the bifurcation point), the equi-librium position loses stability, which leads to self-excitation of parametric resonance oscillations.
The results of the simulation of the instability regions and amplitude-frequency characteristics of vibration machine with a different coefficient of friction values and excitation are given. It is shown that at low coefficients of co-stimulation and increased friction processing load, there is an expansion of the resonance zones due to the small friction in the pendulum of the parametric rotary pendulum pathogen.
The results allow to develop engineering and design activities for the development and creation of technological purpose vibration machines with a stable resonant mode of operation.
Keywords: resonance zone, pendulum, instability region, damping, parametric resonance.
Введение
В настоящее время вибрационные машины технологического назначения пользуются большим спросом и популярностью в различных отраслях промышленности. Наибольшее распространение эти машины получили в производстве тонкодисперсных материалов и нанопорошков.
Современные вибрационные измельчители, выпускаемые как в нашей стране, так и за рубежом, являются прогрессивным видом помольного оборудования. Однако существующая вибропомольная техника не удовлетворяет современным требованиям по производительности, стабильности, технологической эффективности и энергопотреблению.
Одним из наиболее эффективных способов решения указанных проблем является использование фундаментального явления резонанса. Но резонансные режимы работы при обычном резонансе вынужденных колебаний исполнительного органа вибромашины имеют очень низкую стабильность. Такое обстоятельство ограничивает практическое использование резонансных машин. Однако оригинальные результаты [1-4], полученные в ходе работ по возбуждению и стабилизации резонансных параметрических колебаний, дают возможность создания высокопроизводительных
энергосберегающих дробильно-измельчительных машин.
Постановка задачи
Как известно, параметрические колебания -это колебательные движения, возникающие в системе за счет периодического изменения одного из энергоемких параметров системы (жесткости или массы). Характерной чертой параметрических колебаний является то, что возбуждение не происходит, если колебательная система находится в положении равновесия. Однако при определенных отношениях собственной частоты колебаний к частоте возбуждения положение равновесия может оказаться неустойчивым, и сколь угодно малые возмущения могут вызвать параметрические колебания [5].
Как известно, параметрическое возбуждение возможно при соотношении частот [6]
2шп
ш =—1, п = 1, 2, 3,... п
Здесь ю, юп - частота параметрического возбуждения и собственная частота колебаний системы (осциллятора) соответственно. При выполнении указанных соотношений возможно возникновение неограниченно возрастающих колебаний -параметрического резонанса. Однако легче всего параметрический резонанс возникает при ю = 2юп.
Необходимо отметить, что обычный резонанс вынужденных колебаний возникает только при частоте ю = юп, в то время как параметрический резонанс возбуждается вблизи частот ю = 2юп, ю = юп, ю = 2юп/3, ..., то есть при бесчисленном множестве значений частоты ю параметрического возбуждения.
Одной из важнейших задач на этапе проектирования параметрически возбуждаемой вибрационной машины является расчет областей динамической неустойчивости (резонансных зон). В данной работе задачей является определение области динамической неустойчивости вибрационной машины на многократном комбинационном параметрическом резонансе, динамика которой рассмотрена в [7].
Роторно-маятниковый возбудитель 1 этой вибромашины состоит из набора отдельных одинаковых уравновешенных дисков 2 (рис. 1). В каждом диске образована пара незамкнутых беговых дорожек 3 кругового профиля, которые расположены симметрично относительно двух взаимно перпендикулярных его диаметров, а их центры смещены от оси вращения ротора в диаметрально противоположных направлениях на одинаковые расстояния AB = I.
На беговых дорожках размещены одинаковые уравновешенные тела качения (маятники) 4 массой m каждый с возможностью обкатки. Диски соединяются между собой в единую конструкцию так, что беговые дорожки одной пары повернуты вокруг оси ротора на угол уп = тс/я относительно другой, где я - число дисков (я = 2). Ротор возбудителя содержит N = 2я тел качения (маятников), расположенных попарно в параллельных плоскостях и имеющих массу тп, в собранном виде жестко закрепляется на приводном валу, который посредством подшипников устанавливается на рабочем органе 5 массой Мп, имеющем две степени свободы: поступательное движение х, у по круговой траектории в плоскости вращения ротора в направлении координатных осей Ох, Оу. Рабочий орган 5 (например, корпус помольной камеры) связан с основанием 6 упругими элементами 7. Демпфер 8 моделирует трение от технологической нагрузки.
Система координат ААх'у'г с началом в центре масс ротора движется поступательно относительно неподвижной системы Оху1. При этом плоскость Ах 'у' расположена в плоскости вращения ротора. В положении статического равновесия оси этих координатных систем совпадают.
Рис. 1. Параметрически возбуждаемая механическая система
Методы исследования
Положение беговых дорожек определяется углами ук =ш г + 2тск/N (к = 1, 2, 3, 4; N = 4 -
число маятников), а положение маятников определяется углами фк = Акео8(ю1т + 0). Качания маятников на углы фк, а также перемещения х, у рабочего органа составляют степени свободы рассматриваемой механической системы. Эти величины принимаются за обобщенные координаты системы.
При равномерном вращении роторно-маятникового возбудителя на движения маятников накладываются нестационарные связи. При этом уравнениями
хСк = х + I Шв у к + Рс е0^У к + Фк ) ,
уСк = у +1 вт у к +Рс ^(у к +Фк)
описываются все наложенные на представленную систему голономные связи, которые явно содер
жат периодические функции времени периода 2л/ш. Здесь х , у - координаты центра тяжести к-го маятника, рс = ВСк (рис. 1).
При поступательном движении по круговой траектории в направлении координатных осей Ох, Оу рабочий орган преодолевает силы упругости упругой нелинейной изотропной подвески, характеристики которой задаются в виде
Рх = с*х + с!хх\ Ру = ^у + СХуУ\
где сх = су = с - суммарная жесткость упругих элементов в направлении осей х, у соответственно, с1х = с1у = сг - коэффициент нелинейности упругих восстанавливающих сил в соответствующих направлениях.
Кроме сил упругости рабочий орган преодолевает диссипативные силы (силы сопротивления), характеристики которых задаются в виде
Их = Фх + Ь*х2)х, Яу = фу + Ььу2)у, при этом маятники испытывают сопротивление пропорциональное относительной угловой скорости обкатки фк :
Их = (Ьх + Ьих2)х, Яу = (Ьу + Ь^у2)у,
Яфк = (а+ахф2)фк, к = 1,2,..^, где Ьх = Ьу = Ь, а - коэффициенты линейного демпфирования, Ь^х = Ь^ = Ь^ а1 - коэффициенты
нелинейного демпфирования.
Составив выражения для кинетической, потенциальной энергий, диссипативной функции Релея и воспользовавшись уравнением Лагранжа второго рода, запишем дифференциальные уравнения движения представленной вибромашины:
фк 2®2Фк =у2[—01пУк +Фк С0^к -- 2 ф/^1п у к) - у (соэу к - фк яп у к -
1 2 2 2 3 _ 2
-^ф^0^к)] + ® ф* -2(_0 + Кфк)фк,
N
N
~+_ = цоЕ[(фк-<~ фк^1пук + 2<~фкСо^к -
к =1
1 _2 2 _ _ • _ - -Ю ф2С0^к +фкф>к С0^к - 2с~ фк<фк ап Vк +
• 2 ~
+ ф 2С0^ к
1 2-- • _ 1~2 3 • ~ -ф2фк ап Vк + тЮ фк ап Vк
2 6
_ 2 _ 2 _ 3 ~ 2 • -с~ фкф^С08Ук -фкф2^п Vк ] -Р2~ - 2(~ + й~ )—,
_ + ~ = Цо Е[(-фк + <_ фк)С^ук + 2<~ фкук -к=1
1 ~2 2 • ~ ~ ~ Ю фк 81п V к +фк ф к 8Ш V к + 2® ф к ф к С08У к +
• 2 • ~ 1 2- _ + ф 2 ^1п v к + 2 ф2 (рк с0sv к
1 ~2 3 ~
-со фкк c0sVк 6
(1)
-с~фкфу^1п Vк +фкф2С^ Vк ] -Р2~ -
_ —_2 ~
- 2(п + й _ )_,
где к = 1, 2, .., N, N=4, _ = х/1, — = у /1 - безразмерные координаты, Т = А2г - безразмерное время, А2 = ^с / М0 - собственная частота рабочего органа, ц0 = трс /М/, _0 = п0 / А2, _ = п/ А2 - безразмерные коэффициенты линейного демпфирования, й0 = й0 / А2, й = й / А2 - безразмерные коэффициенты нелинейного демпфирования, р1=1/6, Р2 = у /2 - коэффициент нелинейности упругих восстанавливающих сил, с~ = ю / А2 - безразмерная частота параметрического возбуждения, V2 = трс// 1В - безразмерный параметр, определяющий собственную частоту качаний маятников во вращающейся системе координат, у = с1 / с. Здесь М = М0 + т0 + Nm - общая масса системы, п0 = а/2]В, п = Ь/2М, й0 = а1/(2/В), й = Ь1/(2М). Точка обозначает дифференцирование по т.
Примем во внимание, что множитель ц0
^0 = 2в /(^2), где
коэффициент,
пропорциональный отношению общей массы осцилляторов качения к массе всей системы (коэффициент возбуждения). Он определяет существование параметрических резонансных колебаний. Причем чем больше его значение, тем легче возбуждение. Однако с ростом коэффициента возбуждения возрастают инерционные нагрузки на подшипниковые узлы привода, поэтому задача сводится к тому, чтобы как можно меньшей массой маятников раскачать массивный рабочий орган. Слишком малая масса маятников не приведет к возбуждению колебаний, а слишком большая -снизит КПД машины. Расчет областей неустойчивости (резонансных зон) позволит подобрать оптимальные характеристики и параметры механической системы, обеспечивая 100 % выход вибрационной машины на многократный комбинационный параметрический резонанс.
Рассматриваются области динамической неустойчивости стационарного режима работы вибромашины на комбинационном параметриче-
можно записать в виде
в = V2 N1^ / 2 = V2 Nm рс /(2М) -
ском резонансе. Колебания в системе (1) возбуждаются на частотах ю1 и ю2, связанных с частотой параметрического возбуждения ш суммарным соотношением
ш = ш1 +ш2, (2)
где ш1 «А,1 = уш , ш2 «А,2 = 1 - частоты генерации маятников во вращающейся системе координат и рабочего органа соответственно. С помощью приемов метода усреднения [8, 9] получены уравнения для определения амплитуд и частот с учетом того, что (2) имеет вид ш = (1 + А) /(1 — у):
а = 3
8
Р 2|1+—у — ш1 I А2 —р^1 Ап2
— А 2Р±
± (2Яы + п0ы0ы1)'
1 + А
8|--ш]
1 — V
1 — АП 8 п
(3)
4Я0 иыпыш1
— 1,
ш | у(1 + д) + п0ы0ы1 д ш = I--1--А'
2р
1+8 Рх А |. (4)
1 — V пы Здесь А = ш — (уш + 1), расстройка частоты параметрического возбуждения относительно суммы собственных частот:
А2Р=ш2 — 1 — 3Р2ш2А2, ыП = 1 + 1к0Д} ,
ы = 1 +1 к*А2, ы1 = (1 — 3 А?)/(1 —1 А(2),
4 8 8
к* = кпп/ пп, к = к / п, к* = кп/ пп, к = к / п.
Отметим, что уравнения (3), (4), в свою очередь, получены из уравнений первого приближения [7].
Устойчивость системы определяется связанными между собой параметрами е, V, ю, яп, п . Анализ уравнения (3) показывает, что необходимым условием существования колебаний является выполнение неравенства
3 ( 1 2 ^2 пп 8ш211 — — Ад I — 4пппыпыш1 >П. (5)
Тогда максимальная амплитуда определяется
(Атах )2
= 4
1 —
1
4пппыпыш1
3
8ш 2
(6)
У
4«ппыпыш
3
8ш 2
1 < 1,
тшш
которое можно записать через коэффициент возбуждения:
4Яп пыпыш1
8 >-
(7)
Выражение (5) удовлетворяется тогда, когда при линейном демпфировании выполняется пороговое условие возбуждения многократного комбинационного параметрического резонанса. При точной настройке на комбинационный параметрический резонанс ш = уш +1, где ш1 « уш, ш2 «1 оно примет вид
4упп п
8 >-
1 — у
(8)
Из (6) следует, что максимальные амплитуды зависят от коэффициента возбуждения е, а также от линейного и нелинейного демпфирования. Положительные значения амплитуды имеют место при удовлетворении условия
Из (8) видно, что с увеличением параметра V, который может принимать значения 0 < V < 1, коэффициент возбуждения необходимо повышать, увеличивая массу маятников или уменьшая массу рабочего органа вибромашины. Так, коэффициент возбуждения уменьшается в 3 раза при переходе от V = 0,5 к V = 0,25. Поэтому в дальнейших расчетах принимаем V = 0,25. Кроме того, использование упругой изотропной подвески позволяет снизить коэффициент возбуждения в два раза по сравнению с однонаправленной упругой системой [7].
Для численного определения областей неустойчивости можно пользоваться различными методами, например 2^-алгоРитмом [10]. Он основан на численном определении собственных значений матрицы коэффициентов системы уравнений первого приближения. Этот метод неэкономичен с точки зрения производительных затрат машинного времени. В настоящей работе области неустойчивости определяются решением уравнений (3), (4) при помощи ЭВМ с применением современного программного продукта МагкЬаЪ, что позволяет выполнить расчет с высокой точностью.
Анализ результатов
На рис. 2 приведены результаты численного моделирования областей неустойчивости при изотропном закреплении рабочего органа вибромашины с учетом нелинейного демпфирования для случая малого коэффициента возбуждения е = 0,01 и оптимального параметра V = 0,25. По оси абсцисс отложены значения безразмерной частоты параметрического возбуждения ш , а по оси ординат - безразмерный коэффициент яп демпфирования маятников. Области неустойчивости построены для случая изотропной упругой системы, когда юх = юу = ю2, пх = пу = п . Области неустойчивости
заштрихованы.
Анализ областей неустойчивости показывает, что увеличение коэффициента п демпфирования рабочего органа при малых значениях ко-
3
2
эффициента ~0 демпфирования маятников приводит к расширению областей неустойчивости при малом коэффициенте возбуждения. В сравнении с [10] (рис. 3), где е = 0,04, коэффициент возбуждения в 4 раза меньше, поэтому исследование областей неустойчивости при е = 0,01 и с учетом нелинейного демпфирования представляет особый интерес, практическую ценность и восполняет пробелы предшествующих исследований.
0.22 0.180.140.1 0.06 0.02"
«=0.02
«=0.02
«=0.18
использованием высокодобротных материалов, соответствующей термической обработкой маятников и диска роторно-маятникового возбудителя. Кроме того, предпочтительнее [13] использование роторно-маятникового возбудителя в виде приводного вала, получающего вращения от вынесенного из колебательной системы электродвигателя.
"о 0.220.1! 0.14-1 0.1
0.06 0.02
«=0.02 «=0.06
«=0.02
«=0.06
0 1 1.5 ю
Рис. 2. Области неустойчивости на плоскости параметров (ю , ~о) при юх= £ = 0,01,
~0 = ~ = 0,04
Анализ результатов показывает, что ширина резонансных зон зависит от распределения конструктивного (модального) демпфирования по двум подсистемам (маятников и рабочего органа). Увеличение коэффициента демпфирования ~ рабочего органа при малых значениях коэффициента демпфирования ~0 маятников приводит к расширению области неустойчивости, увеличивая зону резонансного режима работы вибромашины. В этом случае при изменении технологической нагрузки самопроизвольного выхода машины из резонансного режима работы не произойдет, чего нельзя сказать о машинах, работающих в режиме обычных вынужденных колебаний [11, 12]. Важно отметить, что уменьшение трения на маятниках легко реализуется на практике за счет конструк-торско-технологических мероприятий, а именно снижения массы маятников и коэффициента трения. Первое мероприятие достигается использованием изотропной упругой подвески. Второе -
0 I 1.5 2 2.5 ш
Рис. 3. Области неустойчивости на плоскости параметров (ю , ~0) при юх= £ = 0,04, к00 = к = 0
Резонансные кривые представлены на рис. 4 в виде зависимостей амплитуды колебаний рабочего органа и маятников от частоты параметрического возбуждения ю .
Все кривые построены с коэффициентом частоты V = 0,25, коэффициентом возбуждения е = 0,01 и различными значениями линейного и нелинейного трения. Резонансная кривая 1 соответствует комбинационному резонансу, при настройке
~0 = ~ = 0,03, ~0 = 0,07, ~ = 0,05, Рх = 1/6, Р2 = 0,15. Кривая 2 соответствует случаю, когда ~0 = ~ = 0,03, к0 = к = 0,06, Рх = 1/6, Р2 = 0. Резонансная кривая 3 построена для величин ~0 = 0,03, ~ = 0,18, ~ = 0,03, ~ = 0,07, Рх = 1/6, Р2= 0.
Приведем анализ полученных результатов. Внутри интервала ш(1) < ю < ю(2) резонансной кри-
вой 1, где ю(1), ю(2) — граничные точки области
)(1), ю(2)
неустойчивости системы (1), состояние равновесия неустойчиво. Точка ю = ю(1) является точкой
иркутским государственный университет путей сообщения
перехода из положения равновесия в положение неравновесного состояния - область динамической неустойчивости [14] (резонансная зона параметрических колебаний). Амплитуда этих колебаний изменяется вдоль верхней ветви резонансной кривой. В точке А1 этой кривой колебания теряют устойчивость, и происходит их срыв. При движении в обратном направлении, начиная с больших значений ю , параметрические колебания можно
возбудить в диапазоне ш(2) <ш только путем жесткого ударного возбуждения. Для этого необходим толчок величиной, превосходящей значение амплитуды нижней (штриховой) ветви резонансной кривой. Граничные точки ш(1) = 1,06,
ш
(2) _
= 1,55 определяются из уравнений (3) и (4) при А = А0 = 0. Верхняя ветвь резонансной кривой 1, изображенной сплошной линией, является устойчивой, а нижняя (штриховая) - неустойчивой. Наклон резонансной кривой характеризует коэффициент нелинейности упругих восстанавливающих сил р2 (для кривой 1 р2 = 0,15). Изменение этого коэффициента приводит к изменению наклона резонансных кривых, вплоть до наклона в противоположную сторону (при р2 < 0).
1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 со
б)
Рис.4. Амплитудно-частотная характеристика рабочего органа и маятников
_ШШ^
В случаи малогистерезисной линейной упругой системы (ß2 = 0) резонансные кривые 2, 3 для амплитуд A и A0 в окрестности частоты возбуждения ш = 1,3, соответствующей точной настройке на комбинационный параметрический резонанс (Ш = vra +1, v = 0,25), представляют собой плавные кривые. Они не имеют характерного пикового значения, присущего высокодобротным колебательным системам при обычном резонансе вынужденных колебаний.
Из рис. 4, а, б (резонансные кривые 2, 3) видно, что шестикратное увеличение п линейного демпфирования рабочего органа и увеличение коэффициента h нелинейного демпфирования рабочего органа при одинаковых малых значениях коэффициента ~0 линейного демпфирования маятников приводит к сдвигу и расширению резонансной зоны. Причем амплитуда колебаний рабочего органа снижается всего лишь в 2,5 раза, тогда как при обычном резонансе вынужденных колебаний линейной системы [15] амплитуда колебаний уменьшается тоже в шесть раз, а наличие демпфирования всегда уменьшает резонансную область.
Заключение
Анализ областей динамической неустойчивости и амплитудно-частотных характеристик показывает, что на рабочий резонансный режим работы вибромашины существенное влияние оказывает трение на маятниках в роторно-маятниковом возбудителе и в упругой системе рабочего органа. При вынужденных колебаниях трение линейно уменьшает их амплитуду. В случае многократного комбинационного параметрического резонанса амплитуда колебаний рабочего органа менее зависима от трения за счет увеличения амплитуды колебаний маятников. При этом частота колебаний рабочего органа в резонансной зоне остается постоянной. Амплитуда плавно нарастает по экспоненциальному закону по мере увеличения частоты вращения роторно-маятникового возбудителя в пределах резонансной зоны, и исчезает при значении правой граничной точки.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Антипов В.И., Денцов H.H., Кошелев A.B. Энергетические соотношения в вибрационной машине на многократном комбинационном параметрическом резонансе // Вестник Нижего-род. гос. ун-та им. Н.И. Лобачевского. 2013. № 5. С. 188-194.
2. Кошелев A.B. Эффективность вибрационной измельчительной машины с параметрическим
IQBQ]
возбуждением // Вестник Машиностроения. 9. 2016. - №5. - с. 27-32.
3. Пат. № 2604005 Рос. Федерация. Вибрационная измельчительная машина / A.B. Кошелев, A.A. 10. Ермолаев. Бюл. № 34.
4. Пат. № 2532235 Рос. Федерация. Вибрационная транспортирующая машина / В.И. Антипов, Р.И. Антипова, A.B. Кошелев, H.H. Денцов. 11. Бюл. № 30.
5. Шмидт Г. Параметрические колебания. М. : 12. Мир, 1978. 336 с.
6. Ильин М.М., Колесников К.С., Саратов Ю.С. Теория колебаний. М. : Изд-во МГТУ им. Н.Э. 13. Баумана, 2003. 272с.
7. Антипов В.И., Денцов H.H., Кошелев A.B. Динамика параметрически возбуждаемой вибра- 14. ционной машины с изотропной упругой системой // Фундаментальные исследования. 2014.
№ 8. Ч. 5. С. 1037-1042. 15.
8. Боголюбов H.H., Митропольский Ю.А. Асимптотические методы в теории нелинейных колебаний. М. : Наука, 1974. 504 с._
Ланда П.С. Автоколебания в системах с конечным числом степеней свободы. М. : Наука, 1980. 360 с.
Антипов В.И. Динамика вибрационных машин с параметрическим возбуждением // Автореф. ... д-ра. техн. наук. Нижний Новгород, Изд-во НГТУ, 2001. 38 с.
Гончаревич И.Ф. Вибрация - нестандартный путь. М. : Наука, 1986. 209 с. Динамика машин и управление машинами : справочник / под ред. Крейнина Г.В. М. : Машиностроение, 1988. 239 с. Вибрация в технике : справочник. Т.4. Вибрационные процессы и машины / под ред. Э.Э. Лавендела. М. : Машиностроение, 1981. 509 с. Гукенхеймер Дж., Холмс Ф. Нелинейные колебания, динамические системы и бифуркации векторных полей. Москва-Ижевск : ИКИ, 2002. Блехман И.И. Вибрационная механика. М. : Физматлит, 1994. 400 с.
УДК 66.023
Трутаев Станислав Юрьевич,
к. т. н., заведующий отделом инновационных разработок,
АО «ИркутскНИИхиммаш», тел. 8(3952)410-336, e-mail: [email protected]
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДЕМПФИРУЮЩИХ УСТРОЙСТВ ВЯЗКОГО ТРЕНИЯ НА ОСНОВЕ РЕЗУЛЬТАТОВ СТЕНДОВЫХ ИСПЫТАНИЙ
S. Yu. Trutaev
MATHEMATICAL MODELING OF VISCOUS FRICTION DAMPER DEVICES ON THE BASIS OF THE BENCH TESTING RESULTS
Аннотация. Продемонстрированы разработки AO «ИркутскНИИхиммаш» в области создания демпфирующих устройств вязкого трения для обеспечения промышленной безопасности оборудования, эксплуатируемого в условиях интенсивных динамических воздействий. Отмечена возможность обеспечения промышленной безопасности динамически нагруженного оборудования на основе применения дискретных демпфирующих связей вязкого или сухого трения. Показана принципиальная конструкция разработанного и внедренного в промышленную эксплуатацию демпфирующего устройства вязкого трения мембранного типа, обладающего высокой чувствительностью к вибрациям малых амплитуд, а также возможностью варьирования демпфирующих свойств в широком диапазоне. Рассмотрены подходы к интерпретации работы реальных демпфирующих устройств вязкого трения с использованием обобщенных моделей Максвелла и Кельвина - Фойгта по результатам стендовых испытаний. В частности, показана возможность идентификации искомых параметров математической модели демп-фирования на основе применения методов оптимального проектирования конструкций. При этом отмечена возможность ис-пользования в качестве минимизируемой целевой функции суммы квадратов невязок между экспериментально измеренной си-лой демпфирования и силой демпфирования, рассчитанной по принятой математической модели. Продемонстрировано разра-ботанное специализированное программное обеспечение, реализующее решение задачи идентификации параметров моделей демпфирования для различного сочетания обобщенных моделей Кельвина - Фойгта и Максвелла с общим количеством цепей до 10. При этом в качестве базового алгоритма поиска оптимума принят метод прямого сканирования по сетке с по-следова-тельным уменьшением зоны поиска.
Ключевые слова: опасный производственный объект, динамическое воздействие, вибрация, сейсмика, демпфирование, демпфер вязкого трения, метод конечных элементов, идентификация, модель Максвелла, модель Кельвина - Фойгта.
Abstract. The article presents the results of work of JSC "IrkutskNIIhimmash " in the field of developing viscous friction dampers. They provide the industrial safety of equipment used under the intensive dynamic impact. The possibility of ensuring industrial safety of dynamically loaded equipment based on the use of discrete damping bonds of viscous or dry friction is noted. The principal design of a viscous friction damping device of membrane type developed and introduced into commercial operation, which has a high sensitivity to vibrations of small amplitudes, as well as the possibility of varying the damping properties in a wide range is shown. The approaches to interpretation of work of the real viscous friction damping devices based on the results of bench testing with the help of generalized Maxwell and Kelvin-Voigt models are considered in the article. In particular, the possibility of identification of the required parameters of damping mathematical model on the basis of the methods of structures optimal design is shown. Besides, it is noted that it is possible
