Contact problem for two-layer elastic coating of solid cylinder

Full Text

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

About the authors

L. V. Bozhkova

Moscow State University of Mechanical Engineering (MAMI)

Email: tm@mami.ru
Dr.Eng., Prof.; +7 (495) 223-05-23, ext. 1318

V. G. Ryabov

Moscow State University of Mechanical Engineering (MAMI)

Email: tm@mami.ru
Ph.D., Prof.; +7 (495) 223-05-23, ext. 1318

G. I. Noricina

Moscow State University of Mechanical Engineering (MAMI)

Email: tm@mami.ru
Ph.D.; +7 (495) 223-05-23, ext. 1318

References

Supplementary files