Модель упругого скольжения в системе «Колесо — поверхность дороги»

УДК 629.1.036

Модель упругого скольжения в системе «колесо — поверхность дороги»1

М. В. Чувейко, С. В. Носачёв

Представлена математическая модель взаимодействия колеса и поверхности дороги. Данная модель позволяет производить качественное исследование особенностей качения колеса для случая отсутствия потери сцепления. Экспериментальные исследования, выполненные в работе, произведены с использованием методов численного моделирования. В полученной математической модели предположено, что жёсткость дороги велика и её деформациями можно пренебречь, а поверхность колеса задана совокупностью точек. Результаты исследования показали отсутствие тривиальной зависимости между угловой частотой колеса и линейной скоростью его центра. Было установлено, что линейная скорость движения колеса при фиксированной угловой частоте существенно зависит от величины внешней нагрузки. Результаты работы представляют интерес в задачах машиностроения, в частности — при построении систем точного позиционирования. Ключевые слова: система «колесо — поверхность дороги», деформация, сила трения, упругое скольжение.

Введение. С момента появления автомобильного транспорта в машиностроительной области возникла задача обоснования природы сил сцепления системы «колесо — поверхность дороги». Данная проблема была рассмотрена специалистами в области механики, триботехники и математики, которые искали ответ в виде решения контактных задач, позволяющих определить размеры пятна контакта и эпюры касательных нормальных напряжений на площадке контакта [1, 2]. В настоящее время в этой области работает ряд отечественных учёных [3-7]. Одним из основоположников нового направления изучения контактных явлений — механики фрикционного взаимодействия (трибомеханики) — является И. Г. Горячева.

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

Основная часть. Рассмотрим задачу о качении колеса по поверхности дороги. Будем полагать, что жёсткость дороги достаточно велика, в связи с чем её деформациями можно пренебречь. Поверхность колеса зададим совокупностью точек упрого-диссипативно связанными между собой и центром колеса. Положение этих точек в отсутствии деформации колеса будем называть их естественным положением. Координаты естественного положения точек могут быть заданы посредством вектора а*, определяемого углами их отклонения от некоторой фиксированной метки М на колесе. Ограничимся моделью рассматривающей динамику п равноотстоящих друг от друга точек (в своём естественном состоянии) поверхности колеса. Тогда координаты вектора а* будут иметь следующие значения:

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

ром а = |А,а2,...,ап]Т.

1 Работа выполнена в рамках инициативной НИР.

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

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

С учётом вышесказанного, моделируемое колесо будет иметь вид, приведённый на рис. 1.

Рис. 1. Контакт «колесо — поверхность дороги»

Здесь г — радиус колеса; 5 — длина контакта «колесо — поверхность дороги»; I — высота центра колеса; Y — угол контакта; ар — угол некоторой точки р поверхности колеса;

хр — расстояние точки р от центра контакта «колесо — поверхность дороги»; ус — скорость

движения центра колеса вдоль горизонтальной оси.

Как видно из рисунка, параметры модели записываются следующим образом:

Из (4) следует, что линейная скорость точки р относительно центра колеса может быть найдена по формуле:

1 - r2 - v

. s Y = arcsin —, 2r

xn = l tanan.

p p

v = dp = i dap pc dt cos2 a„ dt

Тогда абсолютная скорость точки p равна:

l dap

v = v---p

p c cos2 a dt

Сила трения, действующая на точку p со стороны поверхности дороги, может быть найдена по формуле:

Fp = -sign Vp f,

где f — некоторая константа; sign() — функция знака.

Для реализации мягкой функции знака (рис. 2), в модели использована следующая формула:

2

sign(x) = -arctan (kx),

где k — некоторый коэффициент.

-0.3 -0.2

0.2 0.3

Рис. 2. График функции sign

Введём вспомогательные матрицы G и E размера [nx л]:

С =

-2 1 0 0 •• 0 1 1 0 0 0 • 0 0

1 -2 1 0 •• 0 0 0 1 0 0 • 0 0

0 1 -2 1 •• 0 0 0 0 1 0 • 0 0

0 0 1 -2 • 0 0 ; e = 0 0 0 1 • 0 0

0 0 0 0 •• -2 1 0 0 0 0 • 1 0

1 0 0 0 •• 1 -2 0 0 0 0 • 0 1

Тогда матрицы жёсткости C и диссипации H будут иметь вид:

C = с1С - c2E, H = - /E,

где c — коэффициент жёсткости взаимодействия между соседними точками колеса, c2 — коэффициент жёсткости взаимодействия точки и центра колеса, / — коэффициент диссипации взаимодействия между соседними точками колеса, /2 — коэффициент диссипации взаимодействия точки и центра колеса.

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

7d2a „da _(~ da^

J—г = Ca + H--/Fia,— I,

dt2 dt | dt) K 1

где F — векторная функция силы, моделирующая сцепление колеса с дорогой, чьи компоненты,

с учётом (6)—(8), определяются следующим выражением:

1«  с1=

^  эдп

0,

V -1

са , + ср с сН

соб2

а, + а* + ф < Y

\а, + а, + ф > Y

где J — момент инерции точек, ф — угол поворота центра колеса; ф = шН . Динамика центра колеса будет описываться уравнением:

сСУс _ а. (_ сса

т—с = + > Г. I (а,— с вн П ' \ с

где т — масса колеса, Гвн — некоторая внешняя сила, приложенная к центру колеса (вдоль горизонтальной оси).

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

а б

Рис. 3. Угловая частота точки обода колеса

Результат моделирования наглядно подтвердил тот факт, что в случае, когда внешняя нагрузка Гнн препятствует движению, точки колеса, находящиеся в набегающей полуокружности, подвергаются сжатию, а точки с противоположной стороны — растяжению. Это наглядно продемонстрировано на рис. 4, а. В случае сонаправленных векторов скорости движения центра колеса и внешней нагрузки, имеет место обратный эффект (рис. 4, б).

Это явление приводит к появлению эффекта упругого скольжения [5].

По результатам работы математической модели стало видно, что при фиксированной угловой частоте центра колеса в отсутствии классического скольжения (движение точек колеса в зоне контакта), скорость центра колеса в установившемся режиме не является величиной фиксированной и зависит от величины внешней нагрузки Рнн. Для случая на рис. 4, а, скорость точек

колеса в зоне зацепления ниже, чем скорость вращения центра колеса. Это приводит к сгущению точек в контакте. Однако, после выхода точки из зоны зацепления, её угловая частота резко возрастает и накопленное натяжение сводится к нулю (рис. 5).

w^w/ Ш///ШШ/////Ш

а б

Рис. 4. Моделирование движения колеса при сонаправленной внешней силе и скорости

со, С

со

9.5 10 1 < Точка в зоне контакта> < .5 11 11.5 12 Точка вне зоны контакта > 12.5 13 t,

Рис. 5. Зависимость скорости отдельной точки колеса от времени

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

Заключение. При качении колеса по поверхности дороги имеет место упругое скольжение. Таким образом, классическая формула связи угловой и линейной скорости (vc = Rto) в случае неабсолютно жёсткого колеса, является приближённой. Более точное выражение для линейной скорости должно учитывать FBH и параметры колеса.

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

1. Джонсон, К. Механика контактного взаимодействия / К. Джонсон. — Москва : Мир, 1989. — 509 с.

2. Carter, F. W. On the action of a locomotive driving wheel / F. W. Carter // Proc. R. Soc. — London, 1926. — 112А. — Pp. 151-157.

3. Ляпушин, Н. Н. Модель физических процессов в пятне контакта при движении колеса по рельсу со скольжением / Н. Н. Ляпушин, А. Н. Савоськин // Наука и техника транспорта. — 2008. — № 1. — С. 33-42.

4. Максак, В. И. Предварительное смещение и жёсткость механического контакта / В. И. Максак. — Москва : Наука, 1975. — 60 с.

5. Коган, А. Я. Взаимодействие колеса и рельса при качении // Вестник ВНИИЖТ. — 2004. — № 5. — С. 33-40.

6. Лапшин, В. П. Модель связи упруго-вязкого смещения поверхности колеса, относительно рельса с тяговыми характеристиками / В. П. Лапшин, И. А. Туркин, С. В. Носачёв // Вестник Дон. гос. техн. ун-та. — 2012. — № 5. — С. 40-49.

7. Горячева, И. Г. Механика фрикционного взаимодействия / И. Г. Горячева. — Москва : Наука, 2001. — 478 с.

Материал поступил в редакцию 05.02.13.

References

1. Johnson, K. Mexanika kontaktnogo vzaimodejstviya. [Contact interaction mechanics.] Moscow, Mir, 1989, 509 p. (in Russian).

2. Carte, F. W. On the action of a locomotive driving wheel. Proc. R. Soc., London, 1926, 112А, pp. 151-157.

3. Lyapushin, N. N., Savoskin, A. N. Model" fizicheskix processov v pyatne kontakta pri dvizhenii kolesa po rel"su so skol"zheniem. [Model of physical processes in contact patch under wheel sliding motion on the track.] Nauka i texnika transporta, 2008, no. 1, pp. 33-42 (in Russian).

4. Maksak, V. I. Predvaritel"noe smeshhenie i zhestkost" mexanicheskogo kontakta. [Predis-placement and mechanical contact rigidity.] Moscow : Nauka, 1975, 60 p. (in Russian).

5. Kogan, A. Y. Vzaimodejstvie kolesa i rel"sa pri kachenii. [Wheel — rail interaction under rolling.] Vestnik VNIIZhT, 2004, no. 5, pp. 33-40 (in Russian).

6. Lapshin, V. P., Turkin, I. A., Nosachev, S. V. Model" svyazi uprugo-vyazkogo smeshheniya pov-erxnosti kolesa otnositel"no rel"sa s tyagovy"mi xarakteristikami. [Bond model of viscoelastic wheel surface shift relative to rail with pull characteristics.] Vestnik of DSTU, 2012, no. 5, pp. 40-49 (in Russian).

7. Goryacheva, I. G. Mexanika frikcionnogo vzaimodejstviya. [Frictio interaction mechanics.] Moscow : Nauka, 2001, 478 p. (in Russian).

ELASTIC SLIP MODEL IN 'WHEEL — ROAD SURFACE' SYSTEM1

M. V. Chuveyko, S. V. Nosachev

A mathematical model of the wheel — road surface interaction is presented. This model permits to carry out a qualitative investigation on the wheel swing specificity in a coupling case. The experimental investigation is implemented through the numerical modeling methods. In the obtained mathematical mode, it is assumed that the inflexibility of the road is rather great, and its deformations can be neglected. The surface of the wheel is specified by a set of points. The research results have shown a lack of the trivial dependence between the angular frequency of the wheel and the linear speed of its centre. It was established that the linear speed of the wheel movement depends considerably on the value of the external loading. The work results are of interest in the mechanical engineering problems, in particular — in the construction of the exact positioning systems. Keywords:'wheel—road surface' system, deformation, friction force, elastic slip.

1 The research is done within the frame of the independent R&D.

116