CC BY
39
МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ MACHINE BUILDING AND MACHINE SCIENCE
УДК 681.5.017+ 681.5.075
DOI 10.12737/ 11606
Моделирование влияния динамики углового движения колеса на вертикальные колебания
*
транспортного средства
В. П. Лапшин1, В. С. Кравец2, С. А. Конев3**
1 2 3 Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация Simulation of wheel angular movement dynamic impact on vehicle bouncing *** V. P. Lapshin1, V. S. Kravets2, S. A. Konev3**
1 2 3 Don State Technical University, Rostov-on-Don, Russian Federation
Целью работы является повышение эффективности управления транспортными средствами, поэтому помимо синтеза модели приведен пример анализа динамики современного пассажирского электровоза ЭП10. В качестве предмета работы выступают математические модели, отражающие динамику вертикального движения транспортного средства с учетом нелинейного преобразования углового движения колеса в поступательное движение центра масс колеса. Рассматривается контакт колеса и дороги. Представлен синтез модели, отражающей динамику вертикальных колебаний транспортного средства. От известных моделей приведенная в данной статье отличается тем, что в ней учтена динамика углового движения. Таким образом, полученная модель связана с координатами состояния углового движения колеса. Результаты анализа представлены фазовыми портретами, отражающими динамику вертикального движения центра масс колеса, а также динамику вертикальных колебаний кузова. Полученные результаты позволяют сделать следующий вывод. Существует связь между динамикой вертикального и углового движения в транспортном средстве. Как следствие, скорость углового движения влияет на длину площадки контакта колеса с дорогой (рельсом). Результаты могут быть применены при синтезе систем управления транспортными средствами, учитывающих выявленную связанность.
Ключевые слова: динамика транспортного средства, площадка контакта, вертикальные колебания транспортного средства.
As the subject of work, are mathematical mod-els that reflect the dynamics of vertical move-ment of the vehicle, taking into account the nonlinear transformation of the angular motion of the wheel in the forward movement of the center of mass of the wheel, the contact wheel and the road. The article provides a synthesis of models, reflecting the dynamics of vertical oscillations of the vehicle. The difference of the present model from existing consists in taking into account the dynamics of angular motion, resulting in a model associated with the state coordinates the angular motion of the wheel. The aim of this work was to increase the effec-tiveness of the management of modern vehi-cles, so in addition to the synthesis model, the example analysis of the dynamics of a specific vehicle which was organized by the contempo-rary pas-passenger electric locomotive EP. The output of the analysis phase portraits reflecting the dynamics of the vertical motion of the cen-ter of mass of the wheel and the phase portrait, reflecting the dynamics of vertical oscillations of the body. In General, the obtained results allow to make infer-ence about the connectivity of the dynamics of vertical and angular move-ment in the TRANS-tailor the tool and, as a consequence, the effect of the rate of angular movement in the area of wheel contact with the road (rail). As the field of application of the results, see the synthesis of control systems of vehicles, taking into account the revealed relat-edness.
Keywords: vehicle dynamics, ground contact, vertical vibrations of the vehicle.
Введение. В настоящее время довольно остро стоит проблема удешевления железнодорожных и автотранспортных перевозок. Существенно сократить издержки в этой сфере позволит проектирование современных систем автоматического управления транспортными средствами, в контуре управления которыми нет человека. Для современных сложных динамических транспортных систем принципиальное значение приобретают вопросы построения моделей, описывающих контактное взаимодействие колеса и дороги в координатах состояния системы управления [1 -5]. Для синтеза систем управления транспортным средством необходимо решить ряд проблем — в частности, проблему учета
*Работа выполнена в рамках инициативной НИР.
**E-mail: [email protected], [email protected], [email protected]
***The research is done within the frame of the independent R&D.
fa
ш
s
X о о т
43
0
01
X
к
CD
fa
ш
s
X о
Cd
0 Й
01
X S о
влияния динамики вертикальных колебаний центров масс транспортного средства на тяговые и скоростные характеристики процесса управления.
Задача синтеза модели, описывающей вертикальные колебания, не нова. Ее решению посвящено значительное количество работ — как о железнодорожной технике [6, 7], так и об автомобильной [8-11]. Отметим, что в этих моделях не учитывается взаимовлияние координат вращательного движения колеса и координат, описывающих вертикальное движение центра масс колеса и кузова транспортного средства. Вместе с тем, такое влияние подтверждено экспериментально [10, 11]. Это позволяет сделать вывод о необходимости математического описания данного явления. Итак, синтез качественной модели, отражающей указанную связь, является актуальной научной и практической задачей, чему и посвящен представленный материал.
Модель динамики вертикальных колебаний. При построении моделей, описывающих сложную динамику систем управления, часто возникает противоречие между требованием к точности модели и ее сложности. Однако для качественного анализа основных форм колебаний в транспортном средстве достаточно принять двухмассовую модель, которая в упрощенном виде представлена на рис. 1.
а)
7(о)
7(о)
б)
Рис. 1. Модель динамики транспортного средства: в координатах состояния (а); упрощенная двухмассовая модель (б)
й О ТЗ
М
'й
и
и
О, £ -Й
66
На рис. 1 т1, т2 — массы колесной пары и кузова соответственно; с1, с2 — характеристики жесткости колеса и подвески соответственно; ц1, ц2 —коэффициенты, характеризующие диссипацию в колесе и в подвеске соответственно. Горизонтальной пунктирной линией отмечены уровни, соответствующие положению статического равновесия центра масс кузова и колеса. Начало отсчета поместим на этих уровнях, а ось 02 будет направлена вверх. Если центр масс кузова или колеса находится выше положения статического равновесия, то их координаты 21, 22 > 0, если ниже, то 21, < 0. С точки зрения анализа вертикальных колебаний транспортного средства, возникающих в нормальной плоскости, удобно представлять колесо в виде упруго-жесткой связи. Так, к примеру, делают при решении задач моделирования в автомобильной технике [8-11].
Значение деформации подвески с достаточной степенью точности можно определить как:
где 8
т2 Я
А12 = 8стат.2 ~ 22 + =
статическая осадка пружины подвески под действием силы тяжести.
X
2
Для построения модели, отражающей динамику вертикальных колебаний центров масс колеса и кузова автомобиля, используем широко известный подход, предполагающий использование уравнения Лагранжа второго рода [4]. Однако при этом примем в расчет влияние времени существования контакта. Для этого с учетом направления осей х\, введем дополнительную функцию, отражающую эффект «всплывания колеса», который проявляется при увеличении скорости вращения. Подобный эффект достаточно хорошо изучен в динамке металлорежущих станков. Например, в [12, 13] этот эффект объясняется гидродинамической составляющей сил, возникающих в зоне контактирующих поверхностей. В рассматриваемом случае учитывается существование контакта не только в режиме относительного скольжения колеса, но и в режиме качения. Поэтому «всплывание колеса» можно объяснить тем, что контакт колеса и дороги формируется не мгновенно, а в течение некоторого времени. Под временем формирования контакта колеса и дороги мы понимаем время, которое будет потрачено на деформацию в нормальной плоскости колеса. Обозначим время формирования контакта, которое также является и временем существования контакта, величиной А/к. Не вдаваясь в глубокий анализ процессов, происходящих в зоне контакта при сближении поверхностей колеса и дороги, отметим, что существует некоторое эффективное значение А&-, при котором сцепление колеса и дороги достигает некоторого максимального значения и при котором деформация колеса в нормальной плоскости останавливается.
С учетом вышеизложенного и при условии, что Аt < Ак, значение деформационных изменений колеса в нормальной плоскости можно определить:
Д. =5„л - ^ - / (ДО, (1)
где / (А/) — некоторая неизвестная функция, учитывающая влияние времени существования контакта на величину
(ш + ш ) р
деформационных изменений колеса в нормальной плоскости, а б 1 = —1-— статическая деформация (осадка) колеса под действием веса кузова и колеса.
Если А/ > А/к, то
Д, =б , -(2)
1 стат.1 1 V /
При А/ < А/к будут справедливы два следующих предела:
Иш /(Д) ^б„л,
Д/^0 (3)
Иш / (Д/) ^ 0. ()
С достаточной степенью точности такое изменение ^ (Д^) можно описать при помощи следующей функции:
Дк
1 - е^7
(4)
/ (Д) = б
Подставим (4) в (1) получим:
Д,=б ,еД - 2,. (5)
1 стат.1 1 V '
Для дальнейшего анализа уравнения (5) выразим А/ и А/к в координатах состояния. Для этого рассмотрим зависимость времени существования контакта от длины площадки контакта и скорости углового движения колеса. В
К
Пусть Я — радиус колеса; 9 — угол между двумя прямыми, соединяющими центр колеса с двумя противопо- 0 ложными краями площадки контакта; А/. — величина, характеризующая степень деформации колеса, под действием ^ нормально приложенной нагрузки (5); Ь — формируемая длина площадки контакта. 0
С учетом принятых обозначений можно определить угловую величину сектора колеса, принимающего уча- д стие в контакте — 9. Для этого найдем Ь:
Ь = 2^Д/,(2Я-Д.) , (6) ^
и
а затем, используя известные тригонометрические соотношения, определим значение 9:
К
нн
Для железнодорожной техники угол 9/2 достаточно мал. С учетом этого, используя известную теорему о пре- К дельном значении бесконечно малых величин, будет справедливо:
2
ю
. 0 Л/Д/,(2Я -Д.)
¡¡ш- = ---. (7)
ю Я
0= Д/, (2Я -Д.) . (8)
о
Л £ Л
Время формирования площадки контакта можно связать с ее угловой характеристикой (8) и скоростью вращения колеса следующим образом:
0
А/ = -, (9)
ю
где ю — угловая скорость колеса.
Отметим, что контактируют и колесо с дорогой (в железнодорожной технике — с рельсом), и дорога с коле-
П
сом. Время формирования контакта колеса с дорогой будет определяться как А/ = — , а время формирования контакта
ю
Л Ь
дороги с колесом — как А/ = — .
Ь 0
Для режима торможения будет справедливо А/ = —, а для остальных режимов (тяги и боксования) — А/ = — .
V ю
Нелинейная зависимость 9 от Д11 несколько усложняет наши рассуждения, но (8) легко раскладывается в ряд Тейлора. Это позволяет сделать вывод о допустимости линейной интерпретации функции / (А/1) в виде:
/(А/,) = КкА11, (10)
где коэффициент Кк будет определяться углом наклона характеристики /(Д/1) в зоне возможных изменений Д1г и радиусом колеса (для К = 0,5 Кк ~ 14):
а/ = 2КЖ. (П)
Кю
Если известно время формирования контакта Д/к (а мы считаем, что это эмпирически определимая величина), то также известно и значение скоростей юк и V:, то есть это заранее заданные константы. Преобразуем выражение (5) к виду
ю
А/1 =8„ле^ - V (12)
Как видно из (12), от скорости вращения колеса зависит степень его деформации под действием нормально приложенной нагрузки. Это позволяет говорить о взаимосвязанности динамики углового движения колеса (колесной пары в железнодорожной технике) и динамики вертикального колебательного движения в нормальной плоскости центров масс колеса и кузова.
С учетом вышеизложенного и направления осей хьх2 величины деформации колеса и подвески:
ю
А/1 = 5сшашЛе "" - ^ (13)
А/2 = 5сшат.2 - х2 + хР
Уравнение Лагранжа, описывающее динамику взаимосвязанного движения центров масс кузова и колеса в нормальной плоскости, примет следующий вид:
й2 х2 йх2 йх.
т —^ + Ц2 —- + с2 х2 = ц2 —L + с2 х1,
й й/ й/
й2 х. йх, йх2
т1-~2Г + + Ц2 ) + (с1 +С2) х1 = С2 х2 +^2-7" - (14)
ш ш ш
^(т + т2)(1 - в'*) - Ц'5".' в'* ^ .
юк й
Система (14) описывает динамику вертикального движения центров масс кузова и колеса. Отметим, что, если § угловое движение отсутствует, то последние два члена второго уравнения системы будут равны нулю, и система (14) 2 преобразуется к линейному виду.
й Рассматриваемую систему от предложенных ранее [8-11] отличает связь координат состояния, описывающих:
и
и — вертикальное движение центра масс колеса,
— динамику углового движения.
Кроме того, с учетом связи через проскальзывание в данном случае принимается во внимание и поступательное движение.
Результаты численного моделирования. Проанализируем динамику вертикальных колебаний, описываемых системой (14). Для этого в качестве параметров уравнений, входящих в систему, подберем величины, соответствую-
68
щие современному пассажирскому электровозу ЭП10 [1]. В примере, представленном в этой работе, т2 = 50000 кг;
С2 = 1107 Н/м; Ц2 = 40000 Н/(м/с); т = 300 кг; а = 1,6-107 Н/м; щ = 4000 Н/(м/с). Тогда значение 8стат1 = 0,0314 м; ^стат2 = 0,05 м; Кстаб. примем равным 0,99; Юк = 40 рад/с (72 км/ч) при Я = 0,5.
Для проведения эксперимента по оценке динамики транспортного средства в качестве возмущающего воздействия использовались изменения начальных условий по координате х\. То есть моделировалась неровность дороги (рельса), при этом скорость вращения колеса совершала колебания относительно некоторого стационарного значения, что может соответствовать порывам встречного ветра. Результаты моделирования вертикальной реакции колеса представлены рис. 2, 3.
_41_I-1-1--I_I-1-1-■ к
-0 03 -0 02 -0 01 0 0 01 002 0 03 0 04 0 05
21
Рис. 2. Фазовый портрет колебательного движения центра масс колеса при V = 72 км/ч
Реакция на возникающие неровности, представленная фазовым портретом на рис. 2, носит колебательный характер. Взяты значения встречающихся неровностей, не превышающие степень статической деформации колеса (^ стат 1 = 0,0314 м). Однако, как видно из рис. 2, при увеличении скорости транспортного средства на неровной поверхности возможна потеря контакта колеса с дорогой. Из опыта эксплуатации транспортных систем известно, что на неровной дороге транспортное средство должно уменьшать свою скорость. В противном случае возможна потеря управляемости. Для проверки адекватности модели мы увеличили скорость. Результаты моделирования системы представлены на рис. 3.
dz1
¿г
Рис. 3. Фазовый портрет колебательного движения центра масс колеса при V = 135 км/ч
Как видно из рис. 3, с увеличением скорости движения транспортного средства растет вероятность потери управления (контакта колеса с дорогой) на неровной поверхности. Таким образом, адекватность принятой модели подтверждается.
ш
х
X о о т
43
о о X
X сг>
ш
к
X о со о Й о X
к
о
Кузов транспортного средства тоже реагирует на возникающие неровности. Результатом такой реакции должны быть колебания центра масс кузова. Эффект «всплывания» колеса при увеличении скорости движения определяет сложный характер этих колебаний. С точки зрения проверки адекватности модели, правильно давать приращения координат состояния системы относительно стационарного состояния центра масс колеса и контролировать реакцию на такие изменения центра масс кузова (г2). Иными словами, речь идет о моделировании зависимости колебаний кузова от величины возникающих на пути неровностей. Результаты моделирования вертикальных колебаний кузова транспортного средства представлены на рис. 4.
Как видно из рис. 4, кузов транспортного средства совершает более сложное колебательное движение, чем колесо. Однако в отличие от динамики вертикальных колебаний центра масс колеса, центр масс кузова меньше реагирует на колебательные изменения скорости движения транспортного средства, то есть колебания затухают.
Заключение. Представленная в статье модель позволяет численно определять степень влияния скорости транспортного средства на динамику вертикальных колебаний центров масс кузова и колеса. Это позволит проектировать сложные системы управления движением транспорта. Результаты численного моделирования позволяют сделать вывод об адекватности принятой модели. Колебательный характер нарастания скорости транспортного средства обусловливает колебательный характер движения центра масс колеса вблизи равновесного состояния. При увеличении скорости движения увеличивается «всплывание» колеса и отмечается некоторый подъем кузова. При синтезе модели была допущена возможность полного «всплытия» колеса — как в системе (3). Однако проведенные эксперименты (см. рис. 2-4) это не подтвердили, что связано с влиянием подвески. Таким образом, «всплывание» в модели ограничено не только весом, но и упругой силой, связанной с деформацией пружины подвески. Отметим, что рассматриваемый подход к моделированию динамики связанного движения транспортных систем уже известен [14-17], но предлагаемая в работе модель приводиться впервые.
В качестве дальнейшего развития предлагаемого подхода авторы рассматривают синтез связанной модели колебательного движения всего транспортного средства с учетом взаимодействия углового и поступательного движения колеса при проскальзывании.
2 Библиографический список
3 1. Лапшин, В. П. Модель связи вертикальных деформаций с возникновением циркуляционных сил в системах |2 «колесо — рельс» / В. П. Лапшин // Вестник Дон. гос. техн. ун-та. — 2011. — Т. 2, № 8 (59). — С. 1424-1431.
-§ 2. Туркин, И. А. Модель связи упруго-вязкого смещения поверхности колеса относительно рельса с тяговыми
характеристиками / И. А. Туркин, С. В. Носачев, В. П. Лапшин // Вестник Дон. гос. техн. ун -та. — 2012. — Т. 12, № 5 (66). — С. 40-49.
_> 3. Заковоротный, В. Л. Проблемы динамики транспортных трибосистем / В. Л. Заковоротный, В. А. Шапова-
^ лов // Сборка в машиностроении, приборостроении. — 2005. — № 12. — С. 19-24.
£ 4. Заковоротный, В. Л. Исследование комплексного коэффициента трения / В. Л. Заковоротный, В. А. Шапо-
валов // Трение и износ. — 1987. — Т. 9. — № 12. — С. 722-728.
5. Заковоротный, В. Л. Исследование динамических характеристик машины с учетом нелинейности процессов 70 трения / В. Л. Заковоротный, В. А. Шаповалов // Трение и износ. — 1986. — Т. 7, № 4. — С. 681-685.
6. Обобщение передового опыта тяжеловесного движения: вопросы взаимодействия колеса и рельса / У. Дж. Харрис [и др.]. — Москва : Интекст, 2002. — 408 с.
7. Теория электрической тяги / В. Е. Розенфельд [и др.]. — Москва : Транспорт, 1995. — 294 с.
8. Яценко, H. H. Поглощающая и сглаживающая способность шин / H. H. Яценко. — Москва : Машиностроение, 1978. — 131 с.
9. Жуков, A. B. Исследование колебаний автомобилей при переезде неровностей на режимах торможения или разгона / A. B. Жуков // Автомобильная промышленность. — 1973. — № 2. — С. 21-25.
10. Поросятковский, В. А. Исследование параметров вертикальных колебаний автомобильной шины и их влияние на сопротивление качению колеса по грунту : дис. ... канд. техн. наук / В. А. Поросятковский. — Киев, 1974. — 217 с.
11. Смирнов, Г. А. Теория движения колесных машин / Г. А. Смирнов. — Москва : Машиностроение, 1990. —
352 с.
12. Кудинов, В. А. Динамика станков / В. А. Кудинов. — Москва : Машиностроение, 1967. — 359 с.
13. Крагельский, И. В. Фрикционные автоколебания / И. В. Крагельский, Н. В. Гитис. — Москва: Наука, 1987.
— 182 с.
14. Interrelation between tribosystem evolution and parameters of dynamic friction system // V. L. Zakovorotny [et al.] // Journal of Friction and Wear. — 1998. — Т. 19. — № 6. — P. 54-64.
15. Zakovorotny, V. L. Mathematical modelling and control of the evolution of dynamic systems interacting with medium / V. L. Zakovorotny, Nguyen Dong Anh, Pham Dinh Tung // Vietnam Journal of Mechanics. — 2011. — Vol. 33, № 3.
— P. 148-161.
16. Kuznetsov, E. A. Effect of Fluid Lubricant on the Contact Characteristics of Rough Elastic Bodies in Compression / E. A. Kuznetsov // Wear. — 1985. — Vol. 157. — P. 177-194.
17. Ben-David, O. Slip-Stick and the Evolution of Frictional Strength / O. Ben-David, S. M. Rubinstein, J. Fineberg // Nature. — 2010. — V. 463, № 7277. — P. 76-79.
References
1. Lapshin, V. P. Model' svyazi vertikal'nykh deformatsiy s vozniknoveniem tsirkulyatsionnykh sil v sistemakh «koleso — rel's». [Buckling communication model with circulation forces in wheel-rail systems.] Vestnik of DSTU, 2011, vol. 2, no. 8 (59), pp. 1424-1431 (in Russian).
2. Turkin, I. A., Kosachev, S. V., Lapshin, V. P. Model' svyazi uprugo-vyazkogo smeshcheniya poverkhnosti kolesa otnositel'no rel'sa s tyagovymi kharakteristikami. [Bond model of viscoelastic wheel surface shift relative to rail with pull characteristics.] Vestnik of DSTU, 2012, vol. 12, no. 5 (66), pp. 40-49 (in Russian).
3. Zakovorotny, V. L., Shapovalov, V. A. Problemy dinamiki transportnykh tribosystem. [Problems of the transport tribosystems dynamics.] Assembling in Mechanical Engineering and Instrument-Making, 2005, no. 12, pp. 19-24 (in Russian).
4. Zakovorotny, V. L., Shapovalov, V. A. Issledovanie kompleksnogo koeffitsienta treniya. [Study of the complex coefficient of friction. ] Friction and Wear, 1987, vol. 9, no. 12, pp. 722-728 (in Russian).
5. Zakovorotny, V. L., Shapovalov, V. A. Issledovanie dinamicheskikh kharakteristik mashiny s uchetom nelineynosti protsessov treniya. [Investigation of machine dynamic characteristics with account for friction nonlinearity.] Friction and g Wear, 1986, vol. 7, no. 4, pp. 681-685 (in Russian). |
6. Harris W. J., et al. Obobshchenie peredovogo opyta tyazhelovesnogo dvizheniya: voprosy vzaimodeystviya kolesa i H rel'sa. [Best practice analysis of heavy haul traffic: wheel - rail interaction.] Moscow: Intext, 2002, 408 p. (in Russian). с
7. Rosenfeld, V. E., et al. Teoriya elektricheskoy tyagi. [Electric propulsion theory.] Moscow:Transport, 1995, 294 p. тз (in Russian). H
8. Yatsenko, N. N. Pogloshchayushchaya i sglazhivayushchaya sposobnost' shin. [Absorbing and smoothing capacity е of tires.] Moscow: Mashinostroenie, 1978, 131 p. (in Russian). К
9. Zhukov, A. V. Issledovanie kolebaniy avtomobiley pri pereezde nerovnostey na rezhimakh tormozheniya ili ^ razgona. [Study of vehicle vibrations when crossing irregularities in the modes of braking or acceleration.] Avtomobil'naya U promyshlennost', 1973, no. 2, pp. 21-25 (in Russian). о
10. Porosyatkovskiy, V. A. Issledovanie parametrov vertikal'nykh kolebaniy avtomobil'noy shiny
i ikh vliyanie na so- е
protivlenie kacheniyu kolesa po gruntu : dis. ... kand. tekhn. nauk. [Investigation of the parameters of motor tire vertical vibra- е tions and their effect on the ground resistance: Cand.Sci. (Eng.) diss.] Kiev, 1974, 217 p. (in Russian). е
11. Smirnov, G. A. Teoriya dvizheniya kolesnykh mashin. [Theory of wheel carrier movement.] Moscow: Mashinostroenie, 1990, 352 p. (in Russian).
12. Kudinov, V. A. Dinamika stankov. [Dynamics of machines.] Moscow: Mashinostroenie, 1967, 359 p. (in Russian). 71
13. Kragelskiy, I. V., Gitis, N. V. Friktsionnye avtokolebaniya. [Frictional self-oscillations.] Moscow: Nauka, 1987, 182 p. (in Russian).
14. Zakovorotny, V. L., Marchak, M., Usikov, I. V., Lukyanov, A. D. Interrelation between tribosystem evolution and parameters of dynamic friction system. Journal of Friction and Wear, 1998, vol. 19, no. 6, pp. 54-64.
15. Zakovorotny, V. L., Nguyen Dong Anh, Pham Dinh Tung. Mathematical modeling and control of the evolution of dynamic systems interacting with medium. Vietnam Journal of Mechanics, VAST, vol. 33, no. 3 (2011), pp. 148-161.
16. Kuznetsov, E. A. Effect of fluid lubricant on the contact characteristics of rough elastic bodies in compression. Wear, 1985, vol. 157, pp. 177-194.
17. Ben-David, O., Rubinstein, S. M., Fineberg, J. Slip-stick and the evolution of frictional strength. Nature, 2010, no. 7277, vol. 463, pp. 76-79.
Поступила в редакцию 03.03.2015 Сдана в редакцию 04.03.2015 Запланирована в номер 10.04.2015
й о тз
M
"¡3
-M
м
(U >
Л £ -Й
