Например, Бобцов

АВТОНОМИЗАЦИЯ НЕЛИНЕЙНЫХ ДИНАМИЧЕСКИХ СИСТЕМ

С.Е. Иванов, Г.И. Мельников

УДК 531
АВТОНОМИЗАЦИЯ НЕЛИНЕЙНЫХ ДИНАМИЧЕСКИХ СИСТЕМ
С.Е. Иванова, Г.И. Мельникова a Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, Санкт-Петербург, Россия, sivanov@mail.ifmo.ru
Математической моделью многих механических систем является система динамических уравнений полиномиальной структуры с периодическими или постоянными параметрами. Такие механические системы широко применяют в динамике виброзащиты приборов и устройств. Исследование нелинейных систем с конечным числом степеней свободы представляет сложную актуальную проблему по сравнению с линейными системами. Исследование нелинейных систем не сводится к определению конечного числа частных решений, поскольку нелинейные системы не обладают свойством суперпозиции решений. Рассматривается математическая модель нелинейной динамической системы с тремя степенями свободы, которая содержит многочлены до четвертой степени от фазовых переменных. Для исследования данной модели представлены алгоритмические формулы метода автономизации нелинейных динамических систем. Нелинейная математическая модель динамической системы преобразуется к автономной форме, и определяются главные параметры динамической системы. Представлен алгоритм метода автономизации для исследования виброзащитной нелинейной системы с тремя степенями свободы. При решении задач виброзащиты широко применяются линейные системы, хотя линейность функций недостаточно точно аппроксимирует характеристики системы, внося погрешности при анализе. В работе решена задача получения и исследования более точной нелинейной модели виброзащитной системы. Рассмотрена нелинейная виброзащитная система с тремя степенями свободы с нелинейными правыми частями в виде многочлена третьей степени от фазовых переменных с постоянными и периодическими параметрами. Система состоит из объекта виброзащиты, установленного на две платформы, находящиеся одна под другой, нижняя из которых поставлена на вибрирующее основание. Внешнее гармоническое возмущение воздействует на основание. Предполагается, что упругие элементы системы описываются многочленами третьей степени, демпфирующие элементы имеют нелинейную кубическую характеристику. В результате применения метода нелинейная система преобразуется к более простому автономному виду, и существенно сокращается количество параметров нелинейной динамической системы без ухудшения качества решения. Применение метода существенно упрощает исследование переходных и установившихся процессов нелинейных динамических систем. Решаемая задача исследования виброзащитной системы в нелинейной постановке является новой, имеет теоретическое и практическое значение. Ключевые слова: автономизация нелинейных систем, методы исследования, нелинейные системы с тремя степенями свободы.
OFF-LINE INTERACTION OF THE NONLINEAR DYNAMIC SYSTEMS
S. Ivanovb, G. Melnikovb b Saint Petersburg National Research University of Information Technologies, Mechanics and Optics, Saint Petersburg, Russia, sivanov@mail.ifmo.ru
Mathematical model of many mechanical systems is the system of dynamic equations with polynomial structure and periodic or constant parameters. Such mechanical systems are widely applied in the dynamics of a vibration guard of devices. Research of nonlinear systems with a finite number of freedom degrees represents a complicated actual problem in comparison with linear systems. Research of nonlinear systems does not come to definition of a finite number of private solutions because nonlinear systems do not possess superposition property of solutions. Mathematical model of nonlinear dynamic system with three degrees of freedom which contains polynomials to the fourth degree from phase variables is considered. Algorithmic formulas of an off-line interaction method of nonlinear dynamic systems are presented for the given model research. The nonlinear mathematical model of dynamic system is transformed to the autonomous form, and principal parameters of dynamic system are defined. The algorithm of an off-line interaction for research of nonlinear system protected from vibrations with three degrees of freedom is presented. In case of vibration guard solution of problems linear systems are widely used though linearity of functions does not approximate precisely enough the system performance, causing analysis errors. The problem of deriving and research of more exact nonlinear system model protected from vibrations is solved. The nonlinear system protected from vibrations with three degrees of freedom and nonlinear right members in the form of a polynomial of the third degree from phase variables with constant and periodic parameters is considered. The system consists of a vibration guard plant established on two platforms, one under another, the lower one is put on the vibrating foundation. Exterior harmonious perturbation exerts influence upon the foundation. Elastic system elements are supposed to be described by polynomials of the third degree and damping elements have nonlinear cubic performance. As a result of method application the nonlinear system will be transformed to the more simple autonomous form and the amount of parameters of nonlinear dynamic system is essentially reduced. The autonomous system contains essentially less parameters, than initial system, without making worse solution quality. Method application simplifies essentially the research of transitive and settled processes of nonlinear dynamic systems. The solved research problem for the system protected from vibrations in nonlinear statement is new and has theoretical and practical value. Keywords: off-line interaction of nonlinear systems, research methods, nonlinear systems with three degrees of freedom.
Введение
Математической моделью многих механических систем является система динамических уравнений полиномиальной структуры с периодическими или постоянными параметрами. Такие механические системы широко применяют в динамике виброзащиты приборов и устройств [1]. Исследование нелиней-

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics 2014, №1 (89)

151

АВТОНОМИЗАЦИЯ НЕЛИНЕЙНЫХ ДИНАМИЧЕСКИХ СИСТЕМ .

ных систем с конечным числом степеней свободы представляет сложную актуальную проблему по сравнению с линейными системами. Исследование нелинейных систем не сводится к определению конечного числа частных решений, поскольку нелинейные системы не обладают свойством суперпозиции решений.
Предметом рассмотрения является математическая модель нелинейной динамической системы с тремя степенями свободы [2]. Целью работы является приведение математической модели нелинейной динамической системы к автономному виду. Рассматривается математическая модель в виде системы дифференциальных уравнений шестого порядка с нелинейными многочленами четвертой степени относительно фазовых переменных [3]. В результате применения метода автономизации математическая модель динамической системы преобразуется к автономному виду, и определяются главные параметры динамической системы. Метод применен для исследования приборной динамической системы с малыми нелинейными частями [4]. Проведенная проверка полученных результатов и сравнение с решением, полученным численным методом Рунге–Кутта [5], показали точность метода автономизации. В результате применения метода нелинейная система преобразуется к более простому автономному виду, и существенно сокращается количество параметров нелинейной динамической системы без ухудшения качества решения. Применение метода существенно упрощает исследование переходных и установившихся процессов нелинейных динамических систем.

Математическая модель динамической системы

Объектом исследования является математическая модель динамической системы, в виде обыкно-

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

периодическими коэффициентами. Рассматривается система трех дифференциальных уравнений второго

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

Запишем рассматриваемую систему в матричном виде:

4

Aq  Bq  Cq 



cos(ωt)ν1

sin(ωt ) ν2

q qν3 ν4 12

q ν5 3

q1ν6

q2ν7

q3ν8

,

ν 1

(1)

где q  q1, q2 , q3 T – вектор координат, A,B,C – матрицы третьего порядка, ν  (ν1ν2 ν3ν4 ν5ν6 ν7 ν8 ) – век-

торный индекс,   1  2  ... 8 . Предполагается, что система содержит малые нелинейные части с постоянными и периодическими

параметрами hS ν < ε .

Метод автономизации

Приведем алгоритм метода автономизации, который применим для программирования на языке

высокого уровня.

После записи математической модели системы в виде (1) выполним следующие преобразования.

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

q0  exp(iωt) и q0  exp(iωt) , λ1  iω . Запишем в новых переменных (2) периодические функции:

(2)

cos(ωt) 

1 2 (q0

 q0 )

и

sin(ωt) 

1 2i (q0

 q0 ) .

Систему дифференциальных уравнений (1) можно представить в нормальной форме Коши.

W

X



 

0

 A 1H

0 0  A 1C

0  0

I

 

X



 

0

 

,

A1B A1G 

(3)

где фазовый вектор имеет вид X   q0q0q1q2q3q1q2q3 T . Вначале применим линейное преобразование к

системе (3):

Y  DX .
Посредством линейного преобразования (4) система приводится к диагональному виду: Y  ΛY  R XD1Y ,

(4) (5)

где   diag[λ1,λ1,...,λ4 ,λ4 ] , R – преобразованная правая часть. Выполним многочленное преобразование системы (5), которое учитывает многочлены до четвертых степеней:

4

yS  zS 

aνS Ζν

, (s



3,...,8) , Zν



z ν1 1

z ν2 2

...z8ν8

.

ν 2

(6)

В результате преобразования (6) система (5) приводится к автономному виду:

152

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics
2014, №1 (89)

С.Е. Иванов, Г.И. Мельников

4
zS  λS zS  qνS Ζν , (s  3,...,8) . ν 2

(7)

Для нахождения коэффициентов преобразования aνs и преобразованной системы qνs получены формулы (8), удобные для программирования [6]:

4 48

484

     qνS Ζν  (aνS Ζν ( λk νk  λS ))  (aνS Ζν νk zk1 qμk Ζμ )  RS (Z) (s  3,...,8) .

ν 2 ν 2 k 1

ν 2 k 3 μ 2

(8)

При неособых значениях индексов коэффициенты преобразованной системы qνs приравнивают к нулю и при таких значениях вычисляют коэффициенты преобразования aνs . При особых значениях индексов коэффициенты преобразования aνs приравнивают к нулю и вычисляют коэффициенты преобразованной системы qνs . Для нахождения особых значений векторного индекса при фиксированном s необходимо найти решение двух уравнений:

88
 λk νk  λS  0, νk  2, 3, 4 .
k 1 k 1
Произведем следующую замену переменных:

zS  ρS expi(t Im λS  θS ) ; zS1  ρS exp i(t Im λS1  θS ); s  3,5, 7; z1,2  exp(itω).

(9)

В результате замены переменных (9) получим аналитическое решение автономной системы (7):
ρ S  Re(λS )ρS  Re ψS ; θ S  ρS1 Im ψS ; s  3,5, 7;

4

ψS 

q ρ ρ ρS ν3 +ν4 ν5 +ν6 ν7 +ν8 ν3 5 7

exp(i(θ3 (ν3



ν4 )



θ5 (ν5ν6 )



θ7

(ν7



ν8 )



θS )).

ν 2

(10)

Для нахождения стационарных решений приравниваем к нулю правые части системы (10). Для нахождения фазовых координат системы (1) выразим вектор Y по формулам, обратным линейной замене (4):
X  D1Y .

Исследование нелинейной виброзащитной системы

При решении задач виброзащиты широко применяются линейные системы, хотя линейность функций не достаточно точно аппроксимирует характеристики системы, внося погрешности при анализе [7]. В работе решена задача получения и исследования более точной нелинейной модели виброзащитной системы. Решаемая задача виброзащиты в нелинейной постановке является новой, имеет теоретическое и практическое значение.
Применим представленный метод к нелинейной динамической системе с тремя степенями свободы. Рассмотрим виброзащитную систему (рис. 1), состоящую из объекта виброзащиты массой m1 , уста-
новленного на платформы массой m2 и m3 , нижняя из которых закреплена на вибрирующем основании [8]. Предполагается, что упругие элементы системы описываются полиномом третьей степени kx  lx2  px3 , демпфирующие элементы имеют нелинейную кубическую характеристику cx  dx3 .

m1

m2

m3

Рис. 1. Схема нелинейной виброзащитной системы
Математическая модель рассматриваемой системы представима в виде системы дифференциальных уравнений:

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics 2014, №1 (89)

153

АВТОНОМИЗАЦИЯ НЕЛИНЕЙНЫХ ДИНАМИЧЕСКИХ СИСТЕМ .

m1x1  c1(x1  x2 )  d1(x1  x2 )3  k1(x1  x2 )  l1(x1  x2 )2  p1(x1  x2 )3  0, m2 x2  c1 (x2  x1)  d1 (x2  x1)3  k1(x2  x1)  l1(x2  x1)2  p1(x2  x1 )3  c2 (x2  x3 )  d2 (x2  x3 )3  k2 (x2  x3 )  l2 (x2  x3 )2  p2 (x2  x3 )3  0, m2 x3  c2 (x3  x2 )  d2 (x3  x2 )3  k2 (x3  x2 )  l2 (x3  x2 )2  p2 (x3  x2 )3  c3 (x3  f )  d3 (x3  f )3  k3 (x3  f )  l3 (x3  f )2  p3 (x3  f )3  0,
где x1, x2 , x3 – абсолютные перемещения по отношению к положению равновесия системы. На основание действуют вертикальные колебания

(11)

f (t)  asin (ωt).

Введем относительное перемещение: x1  x1  f , x2  x2  f , x3  x3  f .
Запишем уравнения движения (11) в новых переменных (12): m1x1  c1(x1  x2 )  d1(x1  x2 )3  k1 (x1  x2 )  l1(x1  x2 )2  p1 (x1  x2 )3  m1 f1, m2 x2  c1(x2  x1)  d1 (x2  x1)3  k1(x2  x1)  l1(x2  x1 )2  p1(x2  x1)3  c2 (x2  x3 )  d2 (x2  x3 )3  k2 (x2  x3 )  l2 (x2  x3 )2  p2 (x2  x3 )3  m2 f, m2 x3  c2 (x3  x2 )  d2 (x3  x2 )3  k2 (x3  x2 )  l2 (x3  x2 )2  p2 (x3  x2 )3  c3 x3  d3 x33  k3 x3  l3 x32  p3 x33  m3 f.
Здесь f(t)  aω2 sin(ωt).

(12) (13)

Выполним многочленное преобразование системы (13) согласно приведенному алгоритму [9]. C точностью до членов четвертого порядка получаем автономную дифференциальную систему


 

z3 z5

 

(λ3 (λ5



q3 11100000

)

z3

,



q5 11001000

)

z5

,

z4  z3 , z6  z5 ,

 

z7

 (λ7



q7 11000010

)

z5

,

z8  z7

(14)

Решение системы (14) записывается в виде

z3,4

 ρ01 exp(t Re(λ3



q3 11100000

)



i

(θ01

 t Im(λ3



q3 11100000

))),

z5,6

 ρ02

exp(t Re(λ5



q5 11001000

)



i(θ02

 t Im(λ5



q5 11001000

))),

(15)

z7,8

 ρ03

exp(t Re(λ7



q7 11000010

)



i(θ03

 t Im(λ7



q7 11000010

))).

С помощью пакета программ рассчитан переходный и установившийся режимы движения виброзащитной системы при следующих параметрах:
m1  1,13, c1  0, 23, d1  0, 01, k1  1, 23, l1  0, 04, p1  0, 02,
m2  3,17, c2  0, 61, d2  0, 03, k2  3,13, l2  0,13, p2  0, 06,
m3  8, 71, c3  1, 73, d3  0, 09, k3  9,11, l3  0,37, p3  0,17. На систему действует внешнее возмущение: ω  2, a  0,5 .

Определены коэффициенты преобразованной автономной системы (14)

q3 11100000



0, 086  0, 011i,

q5 11001000



0,154  0, 012i,

q7 11000010



0, 044  0, 001i.

Решение системы (15) записывается в виде

z3,4  ρ01 exp(0, 297t  i(θ01 1, 480t)),

z5,6  ρ02 exp(0, 273t  i(θ02 1,130t)),

z7,8  ρ03 exp(0, 083t  i(θ03  0, 635t)).

Установившийся полигармонический режим движения системы имеет вид:

x1  0, 011 0, 035 cos(ωt)  0, 470 sin(ωt)  0, 00001cos(2ωt)  0, 00001sin(2ωt), x2  0, 011 0,113cos(ωt)  0,526 sin(ωt)  0, 00016 cos(2ωt)  0, 00007 sin(2ωt), x3  0, 008  0, 261cos(ωt)  0, 612 sin(ωt)  0, 00045cos(2ωt)  0, 00056sin(2ωt).

154

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics
2014, №1 (89)

С.Е. Иванов, Г.И. Мельников

Колебания системы происходят с частотами, кратными частоте внешней силы, воздействующей на основание. На рис. 2 показано абсолютное перемещение объекта виброзащиты при установлении колебаний нелинейной динамической системы. Представленный график показывает, что вынужденные колебания в начале установления являются квазипериодическими с частотами, кратными внешней силе.
x, мм

0,4
0,2 t, c
0

–0,2

–0,4
0 4 8 12 16 20 24 Рис. 2. Абсолютное перемещение объекта виброзащиты при установлении колебаний
В результате применения метода для анализа виброзащитной системы с тремя степенями свободы получены графики переходных и установившихся режимов виброзащитной системы. Уравнения движения системы приведены к автономному виду. Найдены существенные константы, характеризующие переходные процессы и установившиеся режимы колебаний. Приведена схема метода для систем с тремя степенями свободы, описывается алгоритм программной реализации метода.
Метод позволяет получить достаточно подробные качественные и количественные характеристики изучаемых движений [10], исследовать установившиеся режимы колебаний для систем, находящихся в условиях периодического внешнего воздействия, а также изучать переходные процессы.
Заключение
Рассмотрена нелинейная динамическая система с тремя степенями свободы и нелинейными правыми частями в виде многочлена третьей степени от фазовых переменных с постоянными и периодическими параметрами. Система состоит из объекта виброзащиты, установленного на две платформы, находящиеся одна под другой, нижняя из которых поставлена на вибрирующее основание [11]. Внешнее гармоническое возмущение воздействует на основание. Предполагается, что упругие элементы системы описываются многочленами третьей степени, демпфирующие элементы имеют нелинейную кубическую характеристику [12]. Система уравнений движения приведена к автономному виду. Определены существенные константы, характеризующие переходные процессы и установившиеся режимы колебаний.
В работе представлен алгоритм метода автономизации нелинейных динамических систем с тремя степенями свободы [13]. Разработан алгоритм метода для исследования нелинейных динамических задач. В результате применения метода нелинейная периодическая система преобразуется к автономному виду, и определяются параметры установившихся и переходных режимов движения. Для повышения точности метода в работе [14] предложено применение аппроксимации Чебышева для членов высоких степеней. Приведенные алгоритмические формулы метода позволяют приводить к автономному виду нелинейные динамические системы [15], математическая модель которых представима в виде системы дифференциальных уравнений шестого порядка с нелинейными многочленами до четвертых степеней относительно фазовых переменных.
Литература
1. Бабаков И.М. Теория колебаний. М.: Дрофа, 2004. 591 с. 2. Иванов С.Е. Исследование нелинейных динамических систем с тремя степенями свободы // Научно-
технический вестник СПбГУ ИТМО. 2011. № 4 (74). С. 62–64. 3. Бибиков Ю.Н. Курс обыкновенных дифференциальных уравнений. СПб: Лань, 2011. 304 с. 4. Мельников В.Г. Многочленные преобразования нелинейных систем управления // Изв. вузов. Прибо-
ростроение. 2007. Т. 50. № 5. С. 20–25.

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics 2014, №1 (89)

155

СОВРЕМЕННОЕ СОСТОЯНИЕ И ПЕРСПЕКТИВЫ РАЗВИТИЯ ОСНОВНЫХ ПОНЯТИЙ… .

5. Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа. Приближение функций, дифференциальные и интегральные уравнения. СПб: Лань, 2010. 400 с.
6. Иванов С.Е. Алгоритмическая реализация метода исследования нелинейных динамических систем // Научно-технический вестник информационных технологий, механики и оптики. 2012. № 4 (80). С. 90–92.
7. Матросов В.М., Румянцев В.В., Карапетян А.В. Нелинейная механика. М.: ФИЗМАТЛИТ, 2001. 432 c. 8. Мельников В.Г., Мельников Г.И., Иванов С.Е. Компьютерные технологии в механике приборных
систем: Учеб. пособие / Под ред. В.Г. Мельникова. СПб: СПбГУ ИТМО, 2006. 127 с. 9. Мельников Г.И. Динамика нелинейных механических и электромеханических систем. Л.: Машино-
строение, 1975. 198 с. 10. Иванов С.Е. Определение установившихся режимов работы виброзащитной системы с двумя степе-
нями свободы // Научно-технический вестник СПбГУ ИТМО. 2010. № 4 (68). С. 44–46. 11. Каток А.Б., Хасселблат Б. Введение в теорию динамических систем с обзором последних достижений.
М.: МЦНМО, 2005. 464 с. 12. Hamming R.W. Numerical methods for scientists and engineers. NY: Dover, 1986. 721 p. 13. Мельников В.Г. Энергетический метод параметрической идентификации тензоров инерции тел // На-
учно-технический вестник СПбГУ ИТМО. 2010. № 1 (65). С. 59–63. 14. Мельников В.Г. Преобразование динамических многочленных систем с применением аппроксимации
Чебышева // Научно-технический вестник информационных технологий, механики и оптики. 2012. № 4 (80). С. 85–89. 15. Melnikov V.G. Chebyshev economization in Poincare-Dulac transformations of nonlinear systems // Nonlinear Analysis. 2005. V. 63. N 5–7. P. e1351–e1355.

Иванов Сергей Евгеньевич

– кандидат физ.-мат. наук, доцент, Санкт-Петербургский национальный ис-

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

оптики, Санкт-Петербург, Россия, sivanov@mail.ifmo.ru

Мельников Геннадий Иванович – доктор физ.-мат. наук, профессор, Санкт-Петербургский национальный

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

оптики, Санкт-Петербург, Россия, melnikov@mail.ifmo.ru

Sergei Ivanov Gennady Melnikov

PhD, Associate professor, Saint Petersburg National Research University of Information Technologies, Mechanics and Optics, Saint Petersburg, Russia, sivanov@mail.ifmo.ru D.Sc., Professor, Saint Petersburg National Research University of Information Technologies, Mechanics and Optics, Saint Petersburg, Russia, melnikov@mail.ifmo.ru

УДК 629.78
СОВРЕМЕННОЕ СОСТОЯНИЕ И ПЕРСПЕКТИВЫ РАЗВИТИЯ ОСНОВНЫХ
ПОНЯТИЙ В ОБЛАСТИ МЕХАТРОНИКИ
Е.В. Шалобаевa, Р.-Т. Толочкаb
a Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, Санкт-Петербург, Россия, shalobaev47@mail.ru b Каунасский технологический университет, Каунас, Литва, tadas.tolocka@ktu.lt.
Проанализировано современное состояние базовых терминов в области мехатроники, которое зафиксировано в документах комиссии по стандартизации и терминологии Международной федерации по теории машин и механизмов. Указаны новые термины, которые были введены в научный оборот с начала XXI века. К таким терминам относятся «мехатронизированный объект», «мехатронный класс», объединяющый «мехатронные» и «мехатронизированные» объекты. Дано понятие об «уровневом подходе» к мехатронике, который позволяет увязать между собой собственно мехатронику, микросистемную технику и наноиндустрию. Изложены некоторые соображения о применении термина «мехатронный комплекс», которые стали особенно актуальным в связи с появлением орбитальных группировок малых космических аппаратов. Внесены коррективы по таким критериям оценки мехатронных объектов, как ремонтопригодность и долговечность, а также по использованию термина «мехатронные модули». Поставлены вопросы о расширении понимания мехатроники от компьютерного управления движением до управления состоянием объекта, о терминах «авионика» и «автоника», о взаимосвязи мехатроники и логистики, которую стали относить к предметной области мехатроники. Ключевые слова: мехатроника, мехатронные комплексы, модули и узлы, орбитальная группировка, малые космические аппараты, ремонтопригодность, долговечность, уровневый подход, мехатронизированные объекты.

156

Научно-технический вестник информационных технологий, механики и оптики Scientific and Technical Journal of Information Technologies, Mechanics and Optics
2014, №1 (89)