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

й о тз

M

"¡3

-M

м (U

Ü £ л

МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ MACHINE BUILDING AND MACHINE SCIENCE

УДК 62.752, 621:534.833;888.6 10.23947/1992-5980-2017-17-3-46-59

Устройства для преобразования движения в структуре диады механической колебательной системы*

С. В. Елисеев1, А. И. Орленко2, Д. Х. Нгуен3**

1,3 Иркутский государственный университет путей сообщения, Российская Федерация

2 Красноярский институт железнодорожного транспорта Иркутского государственного университета путей сообщения, Российская Федерация

Motion translation devices in dyad structure of mechanical oscillatory system*** S. V. Eliseev1, A. I. Orlenko2, D. H. Nguyen3**

1,3 Irkutsk State Railway Transport Engineering University, Russian Federation

2 Krasnoyarsk Institute of Rail Transport, Branch of Irkutsk State Transport Engineering University, Russian Federation

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

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

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

Introduction. The possibilities of creating new approaches to the evaluation of the dynamic properties of mechanical oscillation systems are considered. The research objective is to develop a method for constructing mathematical models of mechanical oscillatory systems with several degrees of freedom based on the concepts of system-forming fundamental capabilities of dyads as certain structural formations that determine features of the states of the initial system as a whole.

Materials and Methods. Approaches and techniques of the structural mathematical modeling are used, within which a mechanically oscillating system is compared with a dynamically equivalent structural diagram of an automatic control system. Estimation of the dynamic properties and introduced additional couplings is based on the application of transfer functions, amplitude-frequency characteristics, and frequency analysis technologies. Research Results. The changeability of the dynamic system properties through introducing additional constraints that realize the transformation of states in the interactions of elements on the basis of double differentiation effects is proved. The possibilities of changes that can be initiated by devices for converting motion are shown.

Discussion and Conclusions. A technique for building mathematical models of dyads, and a technology for estimating their dynamic properties are developed. Physically realizable schemes of the design engineering solutions based on motion translation devices are proposed.

Ключевые слова: диада, устройство для преобразования движения, структурная схема, передаточная функция. Keywords: dyad, motion translation device, structural chart, transfer function.

Образец для цитирования: Елисеев, С. В. Устройства для преобразования движения в структуре диады механической колебательной системы / С. В. Елисеев, А. И. Орленко, Д. Х. Нгуен // Вестник Дон. гос. техн. ун-та. — 2017. — Т.17, № 3. — С. 46-59. For citation: S. V. Eliseev, A. I. Orlenko, D. H. Nguyen. Motion translation devices in dyad structure of mechanical oscillatory system. Vestnik of DSTU, 2017, vol. 17, no.3, pp. 46-59.

* Работа выполнена в рамках инициативной НИР (рег. номер AAAA-A16-116112350117-7). 46 ** E-mail: eliseev_s@inbox.ru, Orlenko_AI@krsk.irgups.ru, Huynhnd1987@gmail.com

*** The research is done within the frame of the independent R&D (reg. no. AAAA-A16-116112350117-7).

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

I. Общие положения. Постановка задачи исследования. Рассмотрим механическую линейную колебательную систему с двумя степенями свободы, продемонстрированную на рис. 1. На систему действуют гармонические внешние силы Q1 и приложенные непосредственно к массоинерционным элементам т1 и т2. В системе имеются упругие элементы с жесткостями к1, к2, к3, а также устройство для преобразования движения с приведенной массой Ь [1, 5].

Рис. 1. Принципиальная схема механической колебательной системы, содержащей диаду (mi, кг, mi, L)

Fig. 1. Block diagram of mechanical oscillatory system containing dyad (mi, к2, m2, L)

Используя технологии структурного математического моделирования [3, 5], основанные на применении уравнений Лагранжа 2-го рода и последующих преобразований Лапласа, получим в качестве модели структурную схему эквивалентной в динамическом отношении системы автоматического управления (рис. 2).

Q " ß

Рис. 2. Структурная схема системы с устройствами преобразования движения Fig. 2. Structural diagram of system with motion translation devices

Рассматриваемая система обладает линейными свойствами и совершает малые колебания относительно положения статического равновесия. На рис. 2 приняты следующие обозначения: p = ja — комплексная переменная; значок <-> соответствует изображению переменной по Лапласу. Понятие о приведенной массе L устройства для преобразования движения более подробно рассматривается в [7]. Отметим, что устройство для преобразования движения (УПД) может быть реализовано в различных конструктивно-технических формах, например, в виде рычажных, зубчатых или винтовых несамотормозящихся механизмов.

(U

К X <и

(U

и о X

к 3

CÖ

(U

к

X <и О

а н о о X

к 3

CÖ

Используя структурную схему (рис. 2), можно получить передаточные функции исходной системы, полагая что у и у2 являются выходными сигналами, а внешние воздействия <1, Q2 — входными. Принимая для упрощения,

что Q1 Ф 0 (Q2 = 0), запишем передаточные функции:

^(р) = | = ^2 / + k2 + kз . (1)

W2(p)=<2 ==ЬАтг. (2)

<21 А(р)

где

А(р) = (т1 р2 + Ьр2 + к1 + к2)(т2р2 + Ьр2 + к2 + к3) - (Ьр2 + к2)2 (3)

— характеристическое частотное уравнение системы.

Кроме передаточных функций (1), (2) в рассмотрение вводится передаточная функция межпарциальных связей:

^2^) = ^ = 2 ЬЬ + + к • (4)

у1 т2 р + Ьр + к2 + к3

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

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

II. Диада и ее динамические свойства. Диада рассматривается как структурное образование из типовых элементов с передаточными функциями инерционных звеньев т1, т2, упругого звена к2 и УПД. Отметим, что диада, как некоторое структурное образование, находится «в изоляции» от связей с опорными поверхностями (к1 = 0, к3 = 0). При этом не учитываются силы сопротивления со стороны опорной горизонтальной поверхности.

Структурная схема диады и некоторые формы её преобразований приведены на рис. 3, а, б, в.

й о тз

M

"¡3

-M

м <U

О, £ Л

Qi

а)

( Lp2 + k2)2 {m2 + L)p2 + k2

Г+) 1

(m + L) p2 + k2

tL

У1

Qi

У1

b) c)

Рис. 3. Структурная схема диады: а) общий вид; b) структурная схема диады при исключении координат y2 ; с) структурная схема системы с выделением объекта вибрационной защиты mi

Fig. 3. Dyad structural scheme: а) - general form; b) - dyad structure when excluding coordinates y2 ; c) - system structural scheme with allocation of vibration protection object mi

Передаточные функции диады на основе использования структурных схем (рис. 3, а, б, в) принимают вид:

у т2 р2 + Ьр2 + к2

^ ( р ) = о- =

А( р)

^ (Р) _# =

У2 _ ЬР2 + к2

01 4( р)

^2д(р) _ У- _- Ьр2 +к2

у1 т2 р + Ьр + к2

где

(5)

(6)

(7)

(8)

А1(р) _ (т1 р2 + Ьр2 + к2)(т2р2 + Ьр2 + к2) - (Ьр2 + к2)2

— характеристическое частотное уравнение диады.

Выражение (8) может быть упрощено и приведено к виду:

р4[т1т2 + Ь(т1 + т2)] + р2к2(т1 + т2) _ 0 . (9)

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

р2{[т1т2 + Ь(т1 + т2)]р2 + к2(т1 + т2)} _ 0 . (10)

Из (9), (10) очевидно, что диада имеет одну частоту собственных колебаний с»2д _ 0. Вторая частота собственных колебаний диады определяется выражением:

к2(т1 + т2)

2

Ю2д _-

т1т2 + Ь(т1 + т2) Если считать, что диада имеет без УПД приведенную массу

(11)

(12)

то с учетом (11) приведенная масса диады с УПД составит

т УПД _ т + Ь .

пр УПД пр

(13)

Введение УПД в диаду уменьшает частоту собственных колебаний и вид амплитудно-частотных характеристик (АЧХ) системы. На рис. 4, а, б приведены для сравнения АЧХ обычного вида и АЧХ с введенными УПД.

2x10" 1.45x10"

3.5x10" -2x10" -7.5x10"

- 1.3x10"

- 1.85x10*

- 24*10"

- 2.95x10" -3.5x10"

т • 7Г

*

* т.(1 -Т.С ___ , 1

ч п ) - т со «

'■в- 2 А: б ^^ 4------ й 1 т

3 1

3 /

3 •

3

3

и К X

<и

«

и ю о X

к 3

и

к

X <и О Л

н о о X

к 3

ей

а)

т =

пр

Рис. 4. Амплитудно-частотные характеристики диады: а — АЧХ по координате vi: Ь — АЧХ по координатеуг ( 1 — при L = 0, * *

Fig. 4. Dyad amplitude-frequency characteristics:

a - AFC in coordinate vi; b - AFC in coordinate vj ( ' - at L = 0,.....

— при L = 80 кг)

'•* -мЬ = 80к§)

Для модельной задачи принято, что ш1 = 10 кг, ш2 = 20 кг, к2 = 2000 Н/м; приведенная масса Ь принимает два

Л

значения — Ь = 0 и Ь = 80 кг. Графики амплитудно-частотных характеристик ^(ю) при Ь = 0 и Ь = 80 кг имеют от-

61

личия. Частоты собственных колебаний определяются положением точек (1) и (1') на оси абсцисс. При увеличении Ь

У

значение частоты собственных колебаний уменьшается. При Ь = 0 график ^(ю) пересекается с осью абсцисс в точке

й

(2). При увеличении Ь эта точка перемещается в точку (2') (рис. 4, а). Характерно, что при увеличении Ь значения частоты динамического гашения колебаний и частоты собственных колебаний сближаются.

Графики АЧХ ^(ю) (рис. 4, б) при Ь = 0 и Ь = 80 кг имеют характерные точки (1) и (1') на оси абсцисс, определяющие значения частот собственных колебаний. При Ь = 0 график ^(ю) не имеет пересечений с осью абсцисс.

й

Вместе с тем график -^(ю) при Ь = 80 кг имеет пересечение с осью абсцисс в точке (2'), что свидетельствует о воз-

й

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

Отметим также, что точка (2') на рис. 4, б находится правее точки (1'), что свидетельствует о наличии особых режимов динамического взаимодействия элементов диады. Для диады без УПД АЧХ имеют вид, характерный для систем с одной степенью свободы. При этом полагается, что циклическое движение может на данном этапе исследования не рассматриваться.

С учетом УПД движение по координате у1 имеет частоту динамического гашения:

m2 + L

При этом частота собственных колебаний определяется выражением:

к2

й о тз

"¡3

и

<и

О, £ -Й

mm L +—L-2-

(14)

(15)

При частоте возмущения ю ^ ж , система «запирается» и W1д (р) стремится к значению

^О) = +Ь+Ь + ).

р^ж ш1ш2 + Ь(ш1 + ш2) По координате у2 режим динамического гашения колебаний определится выражением:

(16)

L

соответственно на высоких частотах получим:

^д (р)=-£—г • (18)

т1т2 + Ь(ш1 + т2)

Таким образом, введение УПД привносит в диаду новые динамические эффекты. В частности, возможным становится динамическое гашение колебаний по координате у2 (17), изменяются также значения парциальных частот диады и ее собственной частоты.

Передаточная функция межпарциальных связей с учетом УПД имеет вид, определяемый выражением (7), из которого следуют особенности динамических взаимодействий в межкоординатных связях у и у2. Характерным является «обнуление» (7) или ситуация, когда у2 = 0, что не является возможным в диадах без устройств для преобразования движения. При р ^ 0, соотношения амплитуд = 1, а при р ^ да , = —£—, что можно рассматривать как

У У т2 +£

форму проявления рычажных связей [9, 10].

Одновременное действие двух внешних силовых факторов не изменяет характеристическое частотное уравнение диады, но влияет на формы передаточных функций. Примем, что Q2 = ав1, тогда

Ш (о) = У = (т2 + £) Р2 +к2 +а( £Р 2 +к2) . (19)

1 й А( р) ' 1 '

Ш ' ( р) = У*. = [(т1 + р2 +^2]а+^р2 + к2 . (20)

2д в А( р) ' 1 '

Ш ' ( р) = У2. = [(т1 + £) р2 +^2]а+^р2 +к2 (21)

"12д\^/_ — _ , 2 7 /г 2 7 ч • К^1/

у1 (т2 + £) р +к2 +а(Ьр +к2)

При действии двух силовых факторов частоты собственных колебаний не изменяются. При этом изменяются частоты динамического гашения колебаний. По координате у1 имеем:

2 к2 (1 + а)

т2 + £(1 + а)

По координате у2 соответственно получим:

2 к2 (1 + а)

= (22)

М>иН1> Т ,Л ч (23)

ат1 +£(1 + а)

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

Для случая ю2Шнд =ю2динд можно получить следующее уравнение связности:

к2 (1 + а) = к2 (1 + а)

(24)

т2 + £(1 + а) ат1 + £(1 + а) откуда следует, что

а = (25) §

т1 д

<а

т ^

Если выполняется условие а = —-, то в системе реализуется особый динамический режим, при котором и

т §

К

3

Уг = 1, что соответствует движению объекта по двум координатам как «единое целое», то есть имеет место синфаз- § У S

К

ное движение с одинаковыми амплитудами. Аналогичным образом могут быть получены и другие формы движений. щ

К

к

III. Использование энергетической функции. Если рассматривать диаду с УПД, то в системе координат y1, q

л

y2 ее кинетическая и потенциальная энергии определяются выражениями: н

т =1 my2 +1 rn2y22 +1 L(y -У2)2, (26) |

2 2 2 а

П = 2k2(y2 - У1)2. (27)

Частотная энергетическая функция может быть записана в виде:

k2 (y2 " yi)2

тогда

(mi + L)J>2 + (m2 + L)y2 - 2Lyiy2

Введем понятие коэффициента связности координат:

i = yJ y^,

2 K(i - 1)2

m1 + m2i2 + L(i -1)2

(28) (28') (29)

Рассмотрим, как влияет L на коэффициент формы связи координат i. Если УПД нет, то L = 0, а (29) примет

вид:

2 k2( i - 1) Ю =-;

Сравнивая с частотой собственных колебаний, определенной выше, получим:

к2( I -1)2 k2

m1 + m2 i2 m1m2

(30)

(31)

откуда следует, что

Если в диаду ввести УПД, то (31) трансформируется к виду:

k2( i -1)2 k

m1 + m2i2 + L(i -1)2 l + mm m1 + m2

откуда найдем, что значение i не зависит от L

m.

i = —

m2

(32)

(33)

(34)

На рис. 5 приведены результаты решения уравнения (29) по определению значений коэффициента связности i при различных значениях L (в модельной задаче приняты: m1 = 10 кг, m2 = 20 кг, к2 = 2000 Н/м, L = 0, 10, 20, 50 кг).

агщ 300

I \ ! \

! V..... 1 V 1 д.

L = й . 1 "Д f 1

1

| L28 •1---.V-- 1

120 \

1 = 20 —t-96 1 " '*. \

L = 50 1 ч , \ ч^. \ 3>2> Ч'Д.

1 ' ¡Т. (2)

й о

ТЗ

'й

и (U

Ü С -Й

-Mi - -0.5

1 -1 12

Рис. 5. Графики зависимости ю2(1) при различных значениях L Fig. 5. Dependency graphs ю2(1) at different values of L

Из графиков на рис. 5 очевидно, что динамическое состояние диады не может находиться в форме колебаний, когда i = 1, поэтому все графики ю2(1) при различных значениях L пересекаются в одной точке (1) на оси абсцисс (i = 1); при этом и2о6 = 0 при всех значениях L. В данном случае это означает, что система имеет циклическую координату, характеризующую установившееся равномерное прямолинейное движение. Вместе с тем все графики с различными значениями L имеют экстремальные свойства при i = -0.5, что следует из выражения (34). Таким образом,

m

i = -

m

2

учет Ь (приведенная масса УПД) не приводит к изменению связности координат /, но при этом изменяются соответствующим образом значения частот собственных колебаний диады. При увеличении Ь частоты собственных колебаний диады уменьшаются (в пределе до нулевых значений), что связано с представлениями о получении очень низких частот при Ь ^ да.

При I = 0, что соответствует — = 1, частоты собственных колебаний определяются вышеприведенным выра-

У1

жением (29). При / ^ +да из (29) следует, что

ш2 =

m2 + L

(35)

Отметим, что графики на рис. 5 при i

> будут также стремиться к пределу

+ L

определяемому выра-

жением (35), что соответствует парциальной частоте системы, состоящей из элементов т2, к2 и Ь.

Амплитудно-частотные характеристики (АЧХ) диады при внешнем силовом возмущении к элементу т1 приведены на рис. 6, а, б, в, г, где значение приведенной массы Ь изменяется в пределах Ь = 0 (рис. 6, а); Ь = 10 кг (рис. 6, б); Ь = 20 кг (рис. 6, в); Ь = 50 кг (рис. 6, г).

16x10" &ID"

-SxlO" -1.6*10" -2.4.10" -Si" 10"

ft"

t

w__ '0

I....J 1..... 6 IE »

{

/

1

c) d) Рис. 6. Формы амплитудно-частотных характеристик диады при различных значениях приведенной массы устройства для преобразования движения L: а) L = 0; b) L = 10 кг; c) L = 20 кг; d) L = 50 кг (сплошная линия ( ') соответствует графику зависимости по координате у1; точечная линия (......) соответствует по координате V,)

Fig. 6. Forms of dyad amplitude-frequency characteristics at different values of reduced mass of motion translation device L:

а - L = 0; b - L = 10 kg; c - L = 20 kg; d - L = 50 kg

(solid line ( ') corresponds to dependency graph in coordinate vj; dotted line (......) corresponds to the one in coordinate )

АЧХ на рис. 6, а отображают динамические свойства без УПД при значениях параметров m1 = 10 кг, m2 = 20 кг, к2 = 2000 Н/м. В рассматриваемом случае система имеет одну частоту собственных колебаний (ш2 Ф 0), которая определяется точкой (2) на оси частот. В движениях по координате y1 реализуется режим динамического гашения колебаний точки (1) на оси частот. Движение по координате y2 происходит без режима динамического гашения

y

колебаний, но АЧХ по координате y1 имеет характерную форму с наличием минимального значения

Qi

Введение УПД с приведенными массами L = 0, 10, 20, 50 кг, что отображается на АЧХ, приведенных на рис. 6, приводит к появлению ряда новых динамических эффектов. Это связано с проявлениями режимов динамического гашения колебаний по координате y1. При этом значение частоты собственных колебаний ш2 Ф 0 сдвигается в область более низких значений (точка (1) на рис. 6). На графиках АЧХ при L Ф 0, частоты динамического гашения ко-

<и К X <и

<и и о X

к

3

й

<и К X <и О Л

н о о X

к

3

й

к

2

к

2

лебаний (т. (3) на рис. 6, б, в, г) сдвигаются в сторону расположения точки (2), определяющей частоту собственных колебаний. Отметим, что сближение частот динамического гашения колебаний и собственных колебаний не является приемлемым для решения ряда задач динамики из-за возможностей проявления свойств, характерных для форм неустойчивых движений [11].

IV. Возможности физической реализации диады с использованием рычажных механизмов. Рассмотрим пример реализации диады в структуре динамического корректора колебаний в виброзащитной системе с двумя степенями свободы. На рис. 7 рассматривается механическая колебательная система с двумя степенями свободы как расчетная схема некоторого технического объекта массой да2, который опирается на УПД рычажного типа. Такое устройство является стержнем, на одном конце которого размещен груз массой m0 с возможностями изменения расстояния /2 до точки вращения (точка В). В точке А (рис. 7) рычаг имеет связь в виде кинематической пары V класса (вращательная) с объектом защиты (АВ = А). Одновременно в точке В рычаг соединяется вращательной кинематической парой с промежуточной массой ml. В системе используются четыре упругих элемента с жесткостями Ь, k2, kз и k4 (рис. 7). Движение системы описывается в системе координат у0, у1, у2, связанной с неподвижным базисом; при этом положению объекта защиты m2 соответствует координата у2, положению рычага — у0, промежуточной массы — у1. В качестве внешнего возмущения рассматриваются гармонические колебания опорной поверхности z(t). Система обладает линейными свойствами и совершает малые колебания относительно положения устойчивого статического равновесия.

Рис. 7. Принципиальная схема виброзащитной системы технического объекта с УПД Fig. 7. Block diagram of vibration protection system of technical object with motion translation device В системе координат y\ и y2 выражения для кинетической и потенциальной энергий имеют вид:

т=2 miy2+2 m2y 2+2 moy2' (36)

П = 2k1 U " Z)2 + 2k2 (у2 " yi )2 + 2k3 (o " yi )2 + 2k4 ^y2 " Уо )2 • (37)

С учетом кинематических соотношений сложного движения (36), (37) и введением коэффициента рычажной l,

связи ip = — получим:

li

T =

2 mi yi2 + 2 m2y2 + 2 m0 [ У1 (2 + ip )" y2ip ] 2

й о

ТЗ

"¡3

и (U

о, С -Й

П = i ki (i - z )2 + 2 k2 (У 2 - yi)2 + 2 k3il(yi - У2)2 + 2 k4(i + ip )2 ' (У 2 " yi)2.

В данном случае уравнения движения системы принимают вид:

.yi

+ yi

+ m0 (ip + i y2 (m2 + mfp ) + У2

ki + k2 + k3ip +

+k4 (ip + i)2 k2 + k3ip + +k4 (i + ip )2

- У2m0ip (ip + i)- У2 - yimoip (ip + i)- yi

k2 + k3ip +

+k4 (i + ip) k2 + k3i 1 + +k4 (i + ip)

= ki z(t),

= 0.

(38)

(39)

(40)

(41)

На основании дифференциальных уравнений (40), (41) после преобразований Лапласа [3, 5] может быть построена структурная схема (рис. 8), которая имеет две координаты y и y2. Что касается координаты y0, то она связана с y и y2 следующими соотношениями:

Уо = У i1 + ip )- У2ip. (42)

Коэффициенты уравнений (40) и (41) в операторной форме приводятся в таблице 1.

Таблица 1 Table 1

Коэффициенты уравнений (40), (41) в координатах y, y2 Coefficients of équations (40), (41) in coordinates y, y2

аи а12

m1 + m0 (ip +1) p2 + k1 + k2 + k3i2p + k4 (ip +1) -[mo'p (ip + 1)p2 + k2 + hip, + k4 (1 + ip )2]

а21

-[moip (ip +1) p2 + k2 + k/p + k4 (1 + ip )2] (m2 + m0ip, ) p2 + k2 + Vp + k4 (1 + ip )2

Й 02

k г (t ) 0

Примечание: р = ую — комплексная переменная (у = ); символ <-> над переменной означает её изображение по Лапласу [3, 5].

Структурная схема системы в соответствии с данными таблицы 1 имеет вид, проиллюстрированный на рис. 8. При частоте

к2 + к' + к, (1 +' )2 ю2 = --3 р . Л . р' (43)

меж . /. - \ ^ '

m0lp

(+1)

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

к2 + к3/'2 + к. (1 +' )2

= 2 3р 4 \ р} , (44)

т2 + т01р

П2 к2 + ^ + к4 (1 +'Р )2 (45)

Ю2дин =-ТТ.-Л-. (45)

т0'р ('р + 1)

Н

moip (ip +1) p2 + k2 + + k3i2 + k4(1 + i )2 3 p 4 V p у

1

[ml + m0(ip +1)2] p2 + ki + + k2 + kip + k4(ip +1)2

чи

m0ip (ip + 1)p2 + k2 +

+ k/ + k4(1 + i )2

3 p 4 V ps

1

(m2 + m0ip)p 2 + k2 +

+ kf + k4(1 + i )2

3 p 4 V p у

<D

S X

<u

«

(U

и о X

s

3

й

Рис. 8. Структурная математическая модель механической колебательной системы, изображенной на рис. 7 Fig. 8. Structural mathematical model of mechanical oscillatory system pictured in Fig. 7

Передаточные функции системы при кинематическом возмущении определяются из структурной схемы и имеют вид:

y (m2 + m„/2) p2 + k2 + kJ2 + k. (1 + i )2

W( p) =-yb = —-—-2-—---—, (46)

ki ^ A(p)

(U S X <u О Л H о о X

s 3

ce

У m0ip (i +1) p2 + k2 + k3ip + k4 (1 + ip )2

W2(p) = 0 p p-—-2----p—, (47)

k Z A(p)

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

y m0 ip ( ip +1) p2 + k2 + k3ip + k4 (1 + ip )2 W12 (p) = ^ = 0 p p 2 2-2-Ц---p-j, (48)

y1 (m2 + m0'p )p + k2 + k3lp + k4(1 + 'p )

(49)

где

A( p) = {[ m1 + m0( ip +1)2] p2 + k + k2 + k3/j; + k4( ip +1)2} x x[( m2 + moi^) p2 + k2 + кз/^ + k4(1 + ip )2] - [ moip ( ip +1) p2 + k2 + кз/p + k4(1 + ip )2]2

— характеристическое частотное уравнение системы.

Если /р = 0, то выражения (46), (47) трансформируются к виду:

^(p) = = m2pA+(k2 + k4 , (50)

k1 A ( p)

W2(p) = h = (51)

kl z A (p)

где

A'( p) = [(m! + mo) p2 + k + k2 + k4] • (m2 p2 + k2 + k4) - (k2 + k4)2 (52)

— характеристическое частотное уравнение исходной системы при ip = 0.

Парциальные частоты системы в общем случае определяются выражениями:

2 k1 + k2 + k/p + k4 (ip + 1)2 (53)

n\ =-7--, (53)

m1 + m0( ip +1)

2 k2 + ^ + k4(1 + ip )2

П2 =-+--■

m2 + m0 ip

(54)

При к\ = 0 структурная схема на рис. 8 преобразуется в структурную схему диады, приведенную на рис. 9.

m0ip (ip + 1)p2 + k2 + + кЦ + k4(1 + ip )2

1

К + m0 (ip +1)2] p2 +

+ k2 + кЛ + k4 (ip + 1)2

m0ip (ip + 1)p1 + k2 +

1

( m2 + m0ip2) p 2 + k2 +

+ k3ip2 + k4(1 + ip )2

Ö о

T3 M

"¡3

-m

M

(U

Ü С Л

Рис. 9. Структурная схема диады в общем виде

Fig. 9. Dyad structural diagram in general form

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

[ Ш\ + Шо( ¡р + \)2] р2 + к2 + к3/р + к4( /р + \)2 = Ш\р2 + Шо(1р + \) Р2 + т01р (¡р + \) р2 + ко =

= [Ш\ + Шо(1р + \)] р2 + Шо1р (1р + \) р2 + ко = [Ш\ + Шо(1р + \)] р2 + Ьр2 + ко. В свою очередь, по координате у2 получаем равенство:

(Ш2 + Шо/'р)р2 + к2 + к/р + к4 (\ + 1р )2 = (Ш2 - Шо1р)р2 + Шо^р (1р + \)р2 + ко =

= (Ш2 - Шо/'р)р2 + Ьр2 + ко.

Структурная схема (рис. 9) на основе выражений (55), (56) преобразуется к виду, изображенному на рис. \о.

(55)

(56)

+ Iii + k4 (1 + )

2

3 p 4

z

Рис. 10. Структурная схема диады с устройством для преобразования движения Fig. 10. Dyad structural diagram with motion translation device

Таким образом, в составе расчетной схемы в виде механической колебательной системы с двумя степенями свободы (рис. 7), имеющей рычажный механизм второго рода может быть выделена диада. В диаде кроме обычных массоинерционных элементов m' = m1 + m0 (ip +1) и m' = m2 - m0ip содержится УПД с передаточной функцией

L = m0ip (ip +1), (57)

а также обобщенная пружина с коэффициентом жесткости k0, который соответствует жесткости трех параллельно соединенных пружин:

k = k + kfp + k4( ip +1)2. (58)

При «обнулении» k2, k3, k4 и m0 система превращается в обычную диаду, рассмотренную в предыдущих разделах статьи.

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

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

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

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

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

Библиографический список

1. I. A. Karnovsky, E. Lebed, Theory of vibration protection, doi 10.1007/978-3-319-28020-2. Switzerland, 2016. —

669 p.

2. 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.

3. Елисеев, С. В. Динамическое гашение колебаний : концепция обратной связи и структурные методы мате-магического моделирования / С. В. Елисеев, А. П. Хоменко. — Новосибирск : Наука, 2014. — 357 с.

4. Галиев, И. И. Методы и средства виброзащиты железнодорожных экипажей / И. И. Галиев, В. А. Нехаев, В. А. Николаев. — Москва : Учебно-методический центр по образованию на ж.-д. трансп., 2010. — 340 с.

(U S X

<и

«

(U

и о X

s 3

eö

(U S X <u О

a

H о о X

s 3

eö

5. Елисеев, С. В. Динамический синтез в обобщенных задачах виброзащиты и виброизоляции технических объектов / С. В. Елисеев, Ю. И. Резник, А. П. Хоменко, А. А. Засядко. — Иркутск : ИГУ, 2008. — 523 с.

6. Генкин, М. Д. Упруго-инерционные виброизолирующие системы. Предельные возможности, оптимальные структуры / М. Д. Генкин, В. М. Рябой. — Москва : Наука, 1988. — 191 с.

7. Eliseev S. V., Lukyanov A. V., Reznik Yu. N., Khomenko A. P. Dynamics of mechanical system with additional ties. Irkutsk, 2006. — 315 p.

8. Ким, Д. П. Теория автоматического управления. Т. 1. Линейные системы / Д. П. Ким. — Москва : ФИЗ-МАТЛИТ, 2003. — 288 с.

9. Каимов, Е. В. Оценка динамических свойств в системе с рычажными связями / Е. В. Каимов, А. И. Артюнин // Информационные и математические технологии в науке и управлении : труды XIX Байкальской Всероссийской конф. — Иркутск, 2014. — С. 123-130.

10. Каимов, Е. В. Особенности рычажных связей в схемах подвесок и приборов / Е. В. Каимов, В. Б. Кашуба // Перспективы развития и эффективность функционирования транспортного комплекса юга России : труды межд. научно-практ. конф. — Ростов-на-Дону, 2015. — С. 224-226.

11. Rocard Y. General Dynamics of Vibrations. Paris: Masson, 1949. - 458 p.

References

1. Karnovsky, I.A., Lebed, E. Theory of vibration protection. Switzerland, 2016, 669 p.

2. De Silva, C. W. Vibration. Fundamentals and Practice. Boca Raton, London, New York, Washington, D.C.: CRC Press, 2000, 957 p.

3. Eliseev, S.V., Khomenko, A.P. Dinamicheskoe gashenie kolebaniy: kontseptsiya obratnoy svyazi i strukturnye metody matematicheskogo modelirovaniya. [Dynamic oscillation suppression: feedback concept and structural methods of mathematical modeling.] Novosibirsk: Nauka, 2014, 357 p. (in Russian).

4. Galiev, I.I., Nekhaev, V.A., Nikolaev, V.A. Metody i sredstva vibrozashchity zheleznodorozhnykh ekipazhey. [Methods and means of vibration protection of railway vehicles.] Moscow: Uchebno-metodicheskiy tsentr po obrazovaniyu na zh.-d. transp., 2010, 340 p. (in Russian).

5. Eliseev, S.V., Reznik, Y.I., Khomenko, A.P., Zasyadko, A.A. Dinamicheskiy sintez v obobshchennykh zadachakh vibrozashchity i vibroizolyatsii tekhnicheskikh ob''ektov. [Dynamic synthesis in generalized problems of vibration protection and vibration isolation of technical objects.] Irkutsk: IGU, 2008, 523 p. (in Russian).

6. Genkin, M.D., Ryaboy, V.M. Uprugo-inertsionnye vibroizoliruyushchie sistemy. Predel'nye vozmozhnosti, opti-mal'nye struktury. [Elastic-inertial vibration isolation systems. Limiting possibilities, optimal structures.] Moscow: Nauka, 1988, 191 p. (in Russian).

7. Eliseev S. V., Lukyanov A. V., Reznik Yu. N., Khomenko A. P. Dynamics of mechanical system with additional ties. Irkutsk, 2006, 315 p.

8. Kim, D.P. Teoriya avtomaticheskogo upravleniya. T. 1. Lineynye sistemy. [Automatic control theory. Vol. 1. Linear systems.] Moscow: FIZMATLIT, 2003, 288 p. (in Russian).

9. Kaimov, E.V., Artyunin, A.I. Otsenka dinamicheskikh svoystv v sisteme s rychazhnymi svyazyami. [Estimation of dynamic properties in linkage system.] Informatsionnye i matematicheskie tekhnologii v nauke i upravlenii: trudy XIX Baykal'skoy Vserossiyskoy konf. [Information and Mathematical Technologies in Science and Management: Proc. XIX Baikal All-Russian Conf.] Irkutsk, 2014, pp. 123-130 (in Russian).

10. Kaimov, E.V., Kashuba, V.B. Osobennosti rychazhnykh svyazey v skhemakh podvesok i priborov. [Linkage features in carriages and devices.] Perspektivy razvitiya i effektivnost' funktsionirovaniya transportnogo kompleksa yuga Rossii: trudy mezhd. nauchno-prakt. konf. [Development prospects and efficiency of transport complex in the south of Russia: Proc. Int. Sci.-Pract. Conf.] Rostov-on-Don, 2015, pp. 224-226 (in Russian).

5 11.Rocard, Y. General Dynamics of Vibrations. Paris: Masson, 1949, 458 p.

■д Поступила в редакцию 23.06.2017 "jg Сдана в редакцию 23.06.2017 .> Запланирована в номер 17.07.2017

Received 23.06.2017 Submitted 23.06.2017 Scheduled in the issue 17.07.2017

Об авторах:

Елисеев Сергей Викторович,

директор НОЦ современных технологий, системного анализа и моделирования Иркутского государственного университета путей сообщения (РФ, 664074, г. Иркутск, ул. Чернышевского, 15), доктор технических наук, профессор,

ОЯСГО: http://orcid.org/0000-0001-6876-8786 с^ссу s@inbox.ru

Орленко Алексей Иванович,

директор Красноярского института железнодорожного транспорта Иркутского государственного университета путей сообщения (РФ, 660028, г. Красноярск, ул. Ладо Кецховели, 89), кандидат технических наук, доцент, ОЯСГО: http://orcid.org/0000-0001-6639-9059 Ог1епко AI@krsk.irgups.ru

Нгуен Дык Хуинь,

аспирант Иркутского государственного университета путей сообщения (РФ, 664074, г. Иркутск, ул. Чернышевского, 15),

ОЯСГО: http://orcid.org/0000-0001-9003-5378 HuYnhnd1987@gmai1.com

Authors:

Eliseev, Sergey V.,

director, REC of High Technologies, System Analysis and Modeling, Irkutsk State Railway Transport Engineering University (RF, 664074, Irkutsk, ul. Chernyshevskogo, 15), Dr.Sci. (Eng.), professor, ORCID: http://orcid.org/0000-0001-6876-8786 eliseev s@inbox.ru

Orlenko, Alexey I.,

director, Krasnoyarsk Institute of Rail Transport, Branch of Irkutsk State Transport Engineering University (RF, 660028, Krasnoyarsk, ul. Lado Ketskhovely, 89), Cand. Sci. (Eng.), associate professor, ORCID: http://orcid.org/0000-0001-6639-9059 Orlenko AI@krsk.irgups.ru

Nguyen Duc Huynh,

postgraduate student, Irkutsk State Railway Transport Engineering University (RF, 664074, Irkutsk, ul. Chernyshevskogo, 15),

ORCID: http://orcid.org/0000-0001-9003-5378 Huynhnd1987@gmail.com