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

Полный текст

При уточненном методе решения контактной задачи для кольцевого слоя с учетом сил трения в зоне контакта радиальные контактные давления представлены в виде бесконечного ряда : (1) где: ( - модуль упругости слоя), - величина зоны контакта, и - безразмерные неизвестные константы. Для определения неизвестных констант и получена система двух функциональных уравнений : (2) где: , - радиусы кольцевого слоя , - радиус жесткого цилиндра, с которым контактирует упругий слой (3) - коэффициент трения, , и - некоторая последовательность чисел . Анализ последовательностей чисел и показал, что последовательность чисел стремится к некоторому числу ( , - коэффициент Пуассона упругого слоя), а последовательность чисел стремится к числу . Это позволило улучшить сходимость рядов, входящих в (3), и придать им к следующий вид: (4) где: - некоторая затабулированная функция, - интегральный синус . Как следует из (1) для определения закона распределения радиальных контактных давлений при любом угле необходимо знать две группы констант и , которые входят в систему двух функциональных уравнений (2). Названные константы можно определить двумя способами: либо путем замены системы функциональных уравнений парной бесконечной системой линейных алгебраических уравнений, либо методом ортогонализации систем функций. Вначале охарактеризуем метод замены системы двух функциональных уравнений (2) парной бесконечной системой линейных алгебраических уравнений. С этой целью необходимо выбрать полную систему функций на отрезке , затем умножить обе части функциональных уравнений (2) на функции этой системы и интегрировать их в указанных пределах. В результате получим: (5) где: (6) В качестве системы функций была выбрана система , являющаяся полной на отрезке . Парную систему линейных алгебраических уравнений (5) можно решать либо методом последовательных приближений, либо методом редукции. Рассмотрим теперь способ определения констант и , основанный на построении ортогональных систем функций по заданным системам функций. С этой целью построим прежде всего ортогональную систему функций на базе системы функций : , (7) где: есть « ый» элемент некоторой обратной матрицы . При этом элементы матрицы определяются по формулам: (8) Из первого уравнения системы (2) следует: . (9) На основании (7) и (9) получаем: (10) (11) где: . (12) Подставим (11) в (10), предварительно поменяв в (10) индекс на , а индекс на . В результате получим: (13) где: . (14) Из второго уравнения системы (2) следует: . (15) Подставив (13) в (15), предварительно поменяв в (13) индекс на и наоборот, получим: (16) где: . (17) Для определения постоянных нужно построить ортогональную систему функций на базе системы функций : , (18) тогда: , (19) где: . (20) При этом есть « ый» элемент обратной матрицы , а элементы матрицы определяются по формулам: (21) Таким образом, получены соотношения (7) - (21), с помощью которых можно вычислить все неизвестные константы и при заданном угле и, следовательно, определить закон изменения радиальных контактных давлений по формуле (1). На рисунке 1 показан закон изменения радиальных контактных давлений , полученный при приближенном (пунктирная линия) и уточненном (сплошная линия) методах решения контактной задачи в случае несжимаемого материала и при ; ; . Рисунок 1. Изменение радиальных контактных давлений Анализ численного решения задачи показал, что расхождение между законами изменения радиальных контактных давлений , полученных на основании приближенного и уточненного методов решения, увеличивается с ростом коэффициента трения .

Об авторах

Л. В Божкова

Университет машиностроения

Email: tm@mami.ru
д.т.н. проф.; 8(495) 223-05-23

Г. И Норицина

Университет машиностроения

Email: tm@mami.ru
к.т.н. доц.; 8(495) 223-05-23

В. Г Рябов

Университет машиностроения

Email: tm@mami.ru
к.т.н. проф.; 8(495) 223-05-23

Список литературы

Дополнительные файлы