CC BY
16ISSN 0868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2003, том 13, № 3, c. 77-83
ОРИГИНАЛЬНЫЕ СТАТЬИ
УДК 681.5.01: 658.5; 681.5.01: 658.512 © В. Н. Шамберов
МЕТОД АНАЛИТИЧЕСКОГО ИССЛЕДОВАНИЯ ВЛИЯНИЯ СУХОГО ТРЕНИЯ НА ПОВЕДЕНИЕ АВТОРЕГУЛИРУЕМЫХ ДИНАМИЧЕСКИХ СИСТЕМ
Рассматривается проблема влияния сухого трения в элементах приборных механизмов на возникновение автоколебательных режимов. Задача решается аналитически с использованием вырожденных динамических моделей элемента с трением.
ВВЕДЕНИЕ
Большинство технических устройств, в том числе и приборов, можно отнести к замкнутым ав-торегулируемым динамическим системам. Практический и научный интерес представляют случаи возникновения автоколебаний в таких системах. Возникновение автоколебаний часто связано с проявлением сухого трения в механических элементах этих систем.
Сухое трение не имеет до сих пор вполне достаточного физического объяснения. В связи с этим при решении вопросов анализа и синтеза технических систем трение не всегда корректно учитывается в их динамических моделях, что не позволяет понять и объяснить причины возникновения нежелательного динамического поведения.
Качественный учет физически существенных черт в изучаемом явлении часто позволяет перейти к более простой модели, тем не менее правильно отражающей основные закономерности явления, что делает возможным проведение аналитического исследования.
Учет сухого трения в элементах систем переводит систему в разряд сильно нелинейных. Исследование многомерных нелинейных систем возможно в основном лишь средствами вычислительного эксперимента. Однако зависимость поведения системы от начальных условий и возможность возникновения в системе сложных динамических режимов делает вычислительный эксперимент недостаточно эффективным.
Излагаемый ниже подход позволяет в некоторой степени преодолеть этот барьер, дав исследованию предварительную информацию (полученную на основе точных методов) о возможном поведении системы для определенной области значений ее параметров (коэффициентов) и сделать дальнейшие исследования (уже другими методами) более эффективными.
МОДЕЛИРОВАНИЕ ДИНАМИКИ ЭЛЕМЕНТА С ТРЕНИЕМ
Основными составляющими сил сухого трения (¥с.тр.) являются сила трения покоя (¥ 0) и сила
трения скольжения (¥(У) трск ). Переход от покоя
к скольжению происходит через предварительное смещение, при котором наблюдается наибольшая сила трения покоя [1]. В настоящее время определены устойчивые закономерности в законе сухого трения:
а) сила трения покоя больше силы трения скольжения;
б) переход от покоя к скольжению сопровождается скачкообразным уменьшением силы трения;
в) сила трения покоя увеличивается с увеличением времени неподвижного контакта;
г) с увеличением скорости скольжения (У) сила трения скольжения уменьшается.
Математической моделью закона сухого трения может служить описание:
рс.тр.| < ^тр.0, если Y = 0;
FC1P. = -F(Y) тр.сК. Sign(Y), если Y Ф 0,
(1)
где ¥ (У) трск — монотонно убывающая функция от У, характеризуемая максимальным значением
¥ (У) тр.ск. ^ ¥тр.0Ст. при У ^ 0;
минимальным значением
F (Y)
^ F ■ пои л ^(
тр.ск. тр.шт
максимально отрицательным наклоном — коэффициентом к' = -(а(¥(У)тР.ск.)М7 )|у|^ • 8^(7),
где ¥ 0 > ¥ > ¥ ■ > 0.
^ тр.0 тр. ост. тр.шт
У )тр.ск 1 с.тр. Щтр.0 1 тр. ост. Щ ■ / 1 тр.ШШ , /; /ч arctg(k)
/
агс^(к') 1 тр.шш / У
1 тр. ост. ^^^ У )тр.ск.
-Щтр.0
Рис. 1. Модель закона сухого трения
-0.8 1 0.8
-0 8 0.8
Рис. 2. Фазовый портрет элемента с трением
при т < 0.25к2/ с и к < к'
■1.6 1.6
Рис. 3. Фазовый портрет элемента с трением при т < 0.25к2/ с и к > к'
Долгое время практика учета трения в механических элементах ориентировалась на его простейшие идеализации, что не позволяло получить новые знания об исследуемом явлении и выработать допускающие обобщения концепции.
Динамическую модель элемента с трением с учетом принятого закона сухого трения (1) можно получить на основе принципа динамического равновесия Даламбера (масса вызывает силу инерции, пропорциональную ускорению и противоположно ему направленную):
тУ = 0,
если У = 0 и |сУ + Щ < Щтр0;
тУ + кУ + Щ(У)с тр ск. 81вд(У) + сУ + Щ = 0,
(2)
если У Ф 0;
тУ + кУ + Щтр0 81ви(сУ + Щ) + сУ + Щ = 0,
если У = 0 и |сУ + Щ > ^тр 0,
где т — масса подвижной части элемента с трением; к — коэффициент, характеризующий вязкое трение в элементе; с — коэффициент, характеризующий упругие силы сопротивления движению; Щ — внешнее воздействие. Вид фазового портрета, полученного по уравнениям (2) при
т < 0.25к2 / с, существенно зависит от соотношения коэффициентов к , к' (рис. 2, 3).
Для рассматриваемых практических случаев значение коэффициента т пренебрежимо мало.
В предельном случае (т = 0) движения разделяются на быстрые и медленные. При к < к' на фазовом портрете системы существуют участки неустойчивых медленных движений, которые на практике не реализуются. Любое быстрое движение
переводит систему либо на отрезок равновесия, либо на устойчивое медленное движение. Законность такого рассмотрения, известная ранее как "гипотеза скачка", основывается на доказательствах А.Н. Тихонова, Л.С. Понтрягина, Е.Ф. Мищенко и др. в теории релаксационных колебаний [2].
Полагая, что исходные состояния системы всегда принадлежат либо отрезку равновесия, либо устойчивому медленному движению, аппроксимируя медленные движения на фазовой плоскости линейными отрезками, можно получить математические модели элемента с трением класса кусочно-линейных неоднозначных динамических систем (рис. 4, 5).
Для некоторых элементов помимо массы можно пренебречь также и силами вязкого трения. Получить вырожденную модель (m = 0, k = 0), сохранив в ней все качественные особенности сухого трения, возможно лишь при c > 0 . При этом математическое описание модели будет следующим:
cY =
cY -
если юУ + Л <
тр.0 ■
- Л - Бщп^У -+ Л),
(3)
если \cY + Л >
тр.0
где cY- — предыстория состояния элемента с трением. Геометрическая интерпретация поведения модели (3) при гармоническом входном воздействии представлена на рис. 6.
Л
Рис. 4. Динамическая модель элемента с трением при т = 0 и к < к'
Л
Рис. 5. Динамическая модель элемента с трением при т = 0 и к > к'
У
3
-Л
тр.0
\
и1
^ 1 1 тр.тт ^ ^
-Л
агс^(-с)
Л
тр.0
Л
ж
Рис. 6. Динамическая модель элемента с трением при т = 0 и к = 0
ИССЛЕДОВАНИЕ МНОГОМЕРНЫХ ДИНАМИЧЕСКИХ СИСТЕМ С СУХИМ ТРЕНИЕМ
Применение метода рассмотрим на примере системы, динамика которой может быть представлена следующей математической моделью в век-торно-матричной форме:
тр.тт
х = А х х + В х у,
у = С х х + Б х у, (4)
у = Щу).
Здесь х, х — векторы переменных размерности п х 1; А, В, С — матрицы вещественных коэффициентов размерности п х п , п х т и т х п соответственно; Б — диагональная матрица вещественных коэффициентов размерности т х т ; у, у, у — векторы переменных размерности т х 1; ^у) — вектор нелинейных функций размерности
т х 1. В общем случае п > т > 2 . Составляющие вектора у — переменные 1>1, т>2,..., *т — определяются кусочно-линейными функциями п1(у1),
п2( У2Х пт (Ут ) вида:
п( У) =
0, при |у| < к и п (у) = 0; у, при |у| > к и п- (у) = 0 или при п - (у) Ф 0,
(5)
где g = М 1 х х — вектор-столбец новых переменных, Л = М -1 х А х М —матрица Жордана, В т = М -1 х В , С т = С х М . Преобразование не должно приводить к появлению в В т нулевых строк. Коэффициенты матриц В т и С т определяются как:
Ьт 11 = (т22Ь11 - т12 M,
ьт 12 = (т22¿12 -т^)^М,
Ьт 21 = (т11Ь21 - M,
Ьт 22 = (т11Ь22 - M,
Ст11 С11т11 + С12 т21
Ст12 С11т12 + С12т22 :
С тт С 21 + С 22 21, С тт С 21 ^^12 I С 2 2 ^^
21 12 22 22
где к > 0 — положительный коэффициент, характеризующий величину силы трения покоя; п - (у) — предыстория функции п(у) . Составляющие вектора ^у)— характеризуемые описанием (5) нелинейности пДуД п2(у2), ..., пт(ут) — определены качественным учетом сухого трения в виде силы трения покоя, существующего до момента начала движения при пренебрежении инерционными свойствами элемента с трением, и силами трения скольжения (см. рис. 5).
Для простоты рассуждений рассмотрим систему (4) при п = т = 2 , для которой применим метод, заключающийся в том, что в пространстве параметров системы определим области (сечения), для которых о поведении исходной системы высокого порядка (в данном случае четвертого) с двумя нелинейностями вида (5) можно судить по поведению нескольких (для данного рассмотрения четырех) имеющих одинаковую структуру подсистем второго порядка с одной нелинейностью (5) в каждой. Будем считать, что собственные числа \ , Я2 матрицы А — ненулевые, вещественные, различные и определены, а значения коэффициентов матриц В, С могут выбираться. Линейным преобразованием х = М х g (где М — модальная по отношению к А неособенная матрица) система (1) приводится к виду
£ = Л х g + В т х у,
у = ВД, (6)
у = С т х Ш + б х у,
Линейное преобразование приводит систему к желаемой системе координат (при этом в системе исчезает п(п -1) связей), что позволяет проще выявить ее свойства. Данный прием широко использовался при решении нелинейных задач теории автоматического регулирования [3] для одномерных систем.
Поскольку коэффициенты матриц В, С могут в какой-то степени изменяться (выбираться), возможна ситуация, когда определенная часть коэффициентов матриц В т и С т принимает значения, равные нулю (при этом структура исходной системы (4) не изменяется).
Идея аналитически точного исследования динамики нелинейной системы произвольного порядка с помощью фазовой плоскости теоретически была разработана в работах [4, 5].
Назовем сечением 1 условие Ьт = Ст = 0.
При этом система (6) может быть представлена двумя подсистемами второго порядка:
¿1 =Л & + Ьт 11 V + Ьти V2,
V=пДуД
у1 = Ст 11 ¿1 + ¿11*1 + С т12 & 2;
& 2 =^2 & 2 + Ьт 22 ,
^2 = п2(у 2),
у 2 = Ст 22 & 2 + 422 *2.
Одна из подсистем оказывает внешнее воздействие на другую подсистему в силу ненулевых значений коэффициентов Ь , с , а другая подсистема
не оказывает обратного воздействия в силу нулевых значений коэффициентов Ь , ст1 . Соответствующим образом получаем еще три сечения.
21
22
Сечение 2: (= cmíí = 0). Система (6) приво- дится к виду
дится к виду
¿1 = Я1 ¿1 + Ьт 12 ^ + Ьт22 ^
Т>2 = Щ( у 2),
У 2 = Cm 21 ¿1 + ё 22 V 2 + Cm22 Я 2;
¿ 2 = Я2 ¿ 2 + Ьт 21 V1,
V = п1 (У1),
VI = Cm 12 ¿ 2 + d11V1 .
Сечение 3: (Ьт = cm = 0). Система (6) приво-
дится к виду
¿1 = Я1 ¿1 + Ьт и ^
V = п1(У1),
VI = Cm 21 ¿1 + d11V1 ;
¿2 = Я2 ¿2 + Ьт 21V2 + Ьт21 ^
V2 = П2(У 2),
У2 = Cm 12 ¿ 2 + 422 + Cm21 ¿1.
Сечение 4: (Ьт = cm = 0). Система (6) приво-
¿1 = Я1 ¿1 + Ьт 12 V2,
V2 = П2( У 2),
У 2 = Cm 21 ¿1 + 4 22
¿2 = Я2 ¿2 + Ьт 21V1 + Ьти V2,
V = п1 (У1),
VI = Cm 12 ¿2 + + Cm11 ¿1.
Рассмотренные подсистемы характеризуются своими параметрами Я, Ьт, 4, ^, нелинейностью п(у) и имеют одинаковую структуру
^ = Я + ЬmV,
V = п( У), (7)
У=Cm¿ +
Система (7) с функцией V = п( у) вида (5) была исследована методом точечных отображений [6], при этом были получены необходимые и достаточные условия качественно различного динамического поведения. Структура разбиения пространства параметров Я, Ьт, ё, ^ на области различного динамического поведения показана на рис. 7.
а
ЬmCm
Я2
V
-ё
:й Я
Я2
Ли
б
Я
Рис. 7. Структура разбиения пространства параметров системы (7).
а — для X < 0, б — для X > 0
Допустим, что вязким трением в элементах с трением можно пренебречь (к = 0). В этом случае описание (4) принимает следующий вид:
(8)
х = А X х + В X у(у), У = С х х,
где составляющие вектора у(у) — нелинейные функции У]), У2( У 2),..., Ут (Ут ) вида
V - (у), если |у(у)-- у| < к; у, если У- (у) - У > к.
у( У) =
(9)
¿1 =Я1 ¿1 + Ьт 12 У2(У21
V 2 = Cm 21 ¿1;
^ 2 =Я2 ¿ 2 + Ьт 21 У1(У1) + Ьтп У2(У 2),
VI = '^т 12 ¿ 2 + Cm11 ¿1.
В условиях сечений о динамическом поведении системы (10) и соответственно системы (8) можно судить по динамическому поведению подсистем вида
При этом V (у) — предыстория состояния элемента с трением (см. рис. 6).
Аналогичным преобразованием (при \фЯ2) система приводится к виду
% = Л X g + В т X у(у), у = С т X Ц . (10)
Как и в предыдущем случае, динамическое поведение системы (9) в условиях рассмотренных сечений (при п = т = 2) можно представить:
в сечении 1 ( Ьт 21 = ^ 21 = 0) —
^ 2 =Я2 ¿ 2 + Ьт 22 V2(У2),
V 2 = C
т 22
¿ 2;
^ 1 =Я1 ¿1 + Ьт и V1(Уl) + Ьт12 V2( У 2 ),
VI = Cm 11 ¿1 + Cma ¿ 2;
в сечении 2 ( Ьт11 = Cm11 = 0) —
^ 2 =Я2 ¿ 2 + Ьт 21 V1(Уl),
У1 = Cm 12 ¿ 2 ;
¿1 =Я1 ¿1 + Ьт 12 V2( У 2 ) + Ьт22 V1( Уl),
У 2 = " т 21 ¿1 + "т^ ¿ 2;
в сечении 3 ( Ьт12 = cm2 = 0) —
^ 1 =Я1 ¿1 + Ьт,, V1(Уl),
У1 (;т 21 ¿1;
^ 2 =Я2 ¿ 2 + Ьт 21 V2 (У 2 ) + Ьт21 МуА У 2 = Cm 12 ¿2 + Cm21 ¿1;
в сечении 4 ( Ьт = ^ = 0) —
¿ =Я¿ + ЬтЛ У), У = Cm¿.
(11)
Для системы (11) с функцией v(у) вида (9) также получены необходимые и достаточные условия качественно различного динамического поведения. Структура разбиения пространства параметров Я, Ьт, cm представлена на рис. 8.
В отношении системы (4) можно сделать следующие утверждения:
а) система в условиях сечений устойчива, если две соответствующие подсистемы устойчивы;
б) система в условиях сечений автоколебательна, если одна из подсистем автоколебательна, а другая устойчива; автоколебания определяются автоколебательной подсистемой;
X
ЬmCn
Рис. 8. Структура разбиения пространства параметров системы (11)
в) система в условиях сечения является автоколебательной, если обе подсистемы автоколебательны;
г) судить о поведении системы в условиях сечения по поведению соответствующих подсистем нельзя, если хотя бы одна из них является неустойчивой.
Последнее утверждение связано с тем, что при рассмотрении сечений нулевые значения соответствующих коэффициентов не могут быть абсолютными.
ЗАКЛЮЧЕНИЕ
Развитие теории нелинейных динамических систем непосредственно связано с применением точных аналитических методов исследования. Основное достоинство строгого анализа это представление решения посредством математического выражения — обобщающей идеи, позволяющей перейти на более глубокий уровень понимания явления. Однако налицо и ограниченность такого метода, связанная с размерностью изучаемого динамического явления. Универсальным методом исследования можно было бы считать численные методы исследования с привлечением мощных средств вычислительной техники, если бы они (методы) обладали достаточной общностью и достоверностью представления результатов.
Современной технологией научного познания следует считать обоснованное сочетание точных аналитических методов с численными неаналитическими методами — вычислительный эксперимент. Рассматриваемый в работе метод позволяет эффективно использовать накопленный опыт изучения кусочно-линейных динамических систем невысокого порядка строгими качественными методами исследований применительно к исследованию многомерных динамических систем высокого порядка, осуществляемого средствами вычислительного эксперимента.
СПИСОК ЛИТЕРАТУРЫ
1. Гаркунов Д.Н. Триботехника. М.: Машиностроение, 1985. 424 с.
2. Мищенко Е.Ф., Розов Н.Х. Дифференциальные уравнения с малым параметром и релаксационные колебания. М.: Наука, 1975. 248 с.
3. Лурье А.И. Некоторые нелинейные задачи теории автоматического регулирования (современные проблемы механики). М.—Л.: Гостех-издат, 1951. 216 с.
4. Методы исследования нелинейных систем автоматического управления / Под ред. Р.А. Нелепина. М.: Наука, 1975. 448 с.
5. Нелепин Р.А., Камачкин А.М., Туркин И.И., Шамберов В.Н. Алгоритмический синтез нелинейных систем управления / Под ред. Р.А. Нелепина. Л.: Изд-во ЛГУ, 1990. 220 с.
6. Шамберов В.Н. Исследование типовой промышленной системы автоматического регулирования с некулоновой моделью сухого трения. Автореф. дис. ... канд. техн. наук. Л.: ЛГУ, 1988. 16 с.
7. Шамберов В.Н. Методы аналитического исследования влияния сухого трения на точность и качество авторегулируемых систем // Фундаментальные и прикладные проблемы теории точности процессов, машин, приборов и систем: Труды 5-й сессии Международной научной школы (Фридлендеровские чтения). СПб.: Институт проблем машиноведения РАН, 2002. С.93-100.
Санкт-Петербургский государственный морской технический университет
Материал поступил в редакцию 28.05.2003.
METHOD OF ANALYTICAL RESEARCH INTO THE INFLUENCE OF DRY FRICTION ON DYNAMIC BEHAVIOUR OF AUTOMATIC CONTROL SYSTEMS
V. N. Shamberov
Saint-Petersburg State Marine Technical University
The problem of dry friction impact on autooscillations appearing in adjusting devices is considered. It is solved analytically using singular dynamic models of the elements with dry friction.
