Solution of the contact problem on indentation of rectangular punch in an elastic roughnesses half-space in the presence of coulomb friction

Abstract


The numerical solution of the static three-dimensional contact problem of the indentation of a rectangular stamp with a flat base in an elastic rough half-space in the presence of Coulomb friction and previously unknown adhesion and slip zones is obtained. Accounting for surface roughness in this problem is carried out based on the spherical model of microroughnesses by introducing the nonlinear terms describing surface microroughnesses crushing and shearing to the expression of relative displacement of the interacting bodies. The influence of the values of the friction coefficient and the parameters of the microscopic irregularities on the size and shape of the zone of adhesion and the distribution of the tangential contact stresses are analyzed. It is shown that the inclusion of surface microroughness shear forming roughness can lead to a substantial increase in the size of the zone of adhesion.

Full Text

Введение. При контактировании упругих шероховатых тел возникает трение, вызывающее сдвиг поверхностных микронеровностей, образующих шероховатость. Естественно ожидать, что этот сдвиг может оказать определенное влияние на качественные характеристики контактного взаимодействия © 2014 Самарский государственный технический университет. Образец для цитирования А л е к с а н д р о в А. А., Г р а б к о Е. В. Решение контактной задачи о вдавливании прямоугольного штампа в упругое шероховатое полупространство при наличии кулонова трения // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2014. № 4 (37). С. 42-52. doi: 10.14498/vsgtu1367. 42 Решение контактной задачи о вдавливании прямоугольного штампа . . . тел, на размеры и форму зон сцепления, на распределение касательных контактных напряжений. Однако в большинстве работ, посвященных решению задач о контакте шероховатых упругих тел при наличии трения Кулона [1-9], сдвиг микронеровностей либо не учитывался совсем [1, 3, 4, 7-9] либо учитывался в упрощенной линейной форме [2,5,6]. Кроме этого, условия решаемых задач часто упрощаются за счет игнорирования влияния нормальных контактных напряжений на относительные касательные упругие смещения тел [3, 5-8], либо за счет рассмотрения контакта тел лишь при их полном проскальзывании [4]. Использование нелинейных интегральных уравнений для моделирования контактного взаимодействия упругих тел [3-6, 10, 11] позволяет рассматривать различные типы граничных условий контактных задач и разрабатывать эффективные алгоритмы для получения приближенных решений таких задач. За счет незначительной модификации уравнений [11], описывающих контактное взаимодействие упругих тел при наличии кулонова трения, можно обеспечить учет шероховатости поверхностей этих тел и разработать эффективный алгоритм для решения дискретного аналога модифицированных уравнений. Целью данной статьи является разработка такого алгоритма и выявление эффектов, вызванных учетом сдвига поверхностных микронеровностей, в задаче о контакте штампа с упругим шероховатым полупространством, постановка которой характеризуется отсутствием отмеченных выше упрощений. 1. Постановка контактной задачи. Рассмотрим трехмерную статическую задачу о контакте упругих тел, имеющих шероховатые поверхности. В этой задаче не исключается случай, когда одно из контактирующих тел является абсолютно жестким (в таком случае абсолютно жесткое тело условно считается упругим, имеющим бесконечное значение модуля Юнга). Будем полагать, что взаимодействие тел сопровождается кулоновым трением, и считать поверхность контакта, а также зоны проскальзывания и сцепления на этой поверхности заранее неизвестными. При определенных допущениях такая задача сводится к отысканию неизвестных функций p1 (s), p2 (s), p3 (s), непре¯ рывных на замыкании Ω ограниченной плоской области Ω (охватывающей неизвестную заранее площадку контакта тел) и удовлетворяющих в каждой ¯ точке s множества Ω следующим соотношениям [12]: v1 (s) 0, p1 (s) 0, p2 (s) + p2 (s) 2 3 v1 (s)p1 (s) = 0; µp1 (s); 2 2 v2 (s) + v3 (s)p2 (s) + µp1 (s)v2 (s) = 0; (1) 2 2 v2 (s) + v3 (s)p3 (s) + µp1 (s)v3 (s) = 0. В этих соотношениях через v1 (s), v2 (s), v3 (s) обозначены нормальная и касательные составляющие относительных смещений противолежащих поверхностных точек взаимодействующих тел, функции p1 (s), p2 (s), p3 (s) представляют собой нормальную и касательные составляющие удельной контактной нагрузки; µ - коэффициент трения. Функции v1 (s), v2 (s), v3 (s) имеют следующий вид [13]: 43 А л е к с а н д р о в А. А., Г р а б к о Е. В. 3 A1j (pj ) (s) - ∆1 (s), v1 (s) = f1 p1 (s) + j=1 3 A2j (pj ) (s) - ∆2 (s), v2 (s) = f2 p2 (s), p3 (s), p1 (s) + (2) j=1 3 A3j (pj ) (s) - ∆3 (s). v3 (s) = f2 p3 (s), p2 (s), p1 (s) + j=1 В правых частях равенств (1) слагаемые f1 p1 (s) , f2 p2 (s), p3 (s), p1 (s) , f2 p3 (s), p2 (s), p1 (s) задают сжатие и сдвиг поверхностных микронеровностей полупространства, образующих шероховатость; функции ∆1 (s), ∆2 (s), ¯ ∆3 (s) ∈ C Ω [13] задают конфигурацию взаимодействующих тел и условия ¯ их нагружения; линейные интегральные операторы влияния Aij : C Ω → ¯ C Ω определены соотношениями Aij (pj ) (s) = Kij (s, t)pj (t)dt, i, j = 1, 2, 3. Ω В этих соотношениях ядра Kij (s, t) заданы в соответствии с формулами Буссинеска-Черутти [14], что означает принятие гипотезы о возможности аппроксимации взаимодействующих тел упругими полупространствами. Первое из условий (1) означает отсутствие взаимного проникания тел, второе выражает знакопостоянство контактного давления, и третье означает, что за пределами площадки контакта давление отсутствует. Последние три соотношения системы (1) выражают закон трения Кулона [12]. При формулировании этого закона трения принята гипотеза [15], позволяющая заменить в последних двух соотношениях системы (1) скорости относительного проскальзывания тел их относительными проскальзываниями. ¯ Для неизвестных функций p1 (s), p2 (s), p3 (s) ∈ C Ω система соотношений (1), (1) эквивалентна следующей системе уравнений [12]: p1 (s) = h p1 (s) - Ev1 (s) ; p2 (s) = q p2 (s) - Ev2 (s), p3 (s) - Ev3 (s), µh (p1 (s)) ; p3 (s) = q p3 (s) - Ev3 (s), p2 (s) - Ev2 (s), µh (p1 (s)) , (3) ¯ где s ∈ Ω; E - произвольное положительное число; выражения v1 (s), v2 (s), v3 (s) имеют вид (1), а функции h(x), q(x, y, z) заданы следующим образом: 1 h(x) = (x + |x|); 2 q(x, y, z) =      x, xz x2 + y 2 if , if x2 + y 2 z; x2 + y 2 > z. 2. Учет шероховатости поверхности упругого полупространства. Используя сферическую модель микровыступов [16], образующих шероховатость, и 44 Решение контактной задачи о вдавливании прямоугольного штампа . . . пренебрегая влиянием касательной нагрузки, приложенной к этим выступам, на их сжатие, можно получить для f1 , f2 следующие соотношения [17-19]: 2/3 f1 (u) = α h(u) ,   0,   f2 (x, y, z)= 1 - 1 - ( x2 + y 2 /˜) ˜ ˜ z βx ˜    (˜2 + y 2 )1/6 x ˜ ( x2 + y 2 /˜)2/3 ˜ ˜ z if x2 + y 2 = 0; ˜ ˜ (4) 2/3 , if x2 ˜ + y2 ˜ > 0; где x = q x, y, µh(z) , y = q y, x, µh(z) , z = µh(z). В этих соотношениях па˜ ˜ ˜ раметры α и β имеют одинаковую размерность и определяются следующими выражениями: α = 0.8255ρ β = 0.4127ρ 3 3 ˜ kπ 2 µ 2 1-ν1 E1 + 2 1-ν2 E2 2 2 1 - ν1 1 - ν2 ˜ kπ 2 + E1 E2 2 ; (5) (1 + ν1 )(2 - ν1 ) (1 + ν2 )(2 - ν2 ) + . E1 E2 В равенствах (1) параметры ν1 , ν2 , E1 , E2 представляют собой коэффициенты Пуассона и модули Юнга взаимодействующих тел, параметр ρ есть радиус ˜ сферического выступа, параметр k принимает значения 2 или 1 (в зависимости от того, оба ли взаимодействующих тела имеют шероховатые поверхности или лишь одно из них). Если одно из взаимодействующих тел является абсолютно жестким, то в равенствах (1) нужно положить равными нулю все дроби, в знаменателях которых содержится модуль Юнга абсолютно жесткого тела. 3. Алгоритм численного решения контактной задачи. Для получения численного решения системы интегральных уравнений (1) зададим область Ω в виде открытого прямоугольника, ограниченного отрезками прямых, параллельных координатным осям x и y декартовой системы координат, введенной на общей для взаимодействующих тел касательной плоскости. Разобьем Ω на k непересекающихся одинаковых прямоугольных областей ω1 , ω2 , . . . , ωk , ориентированных подобно прямоугольнику Ω (размеры области Ω выберем так, чтобы эта область включала в себя прямоугольную подошву штампа). Полагая, что искомые функции p1 (s), p2 (s), p3 (s) принимают на каждом граничном элементе ωi постоянные значения x3i-2 , x3i-1 , x3i (i = 1, 2, . . . , k), можно для определения этих значений получить из (1) следующую систему 3k скалярных уравнений [19]: 3k ; j=1 a3i-2j xj - b3i-2 x3i-1 = q x3i-1 - E f2 (x3i-1 , x3i , x3i-2 ) + 3k a3i-1 j xj - b3i-1 j=1 3k x3i - E f2 (x3i , x3i-1 , x3i-2 ) + j=1 a3i j xj - b3i , µx3i-2 ; x3i = q x3i - E f2 (x3i , x3i-1 , x3i-2 ) + 3k a3i j xj - b3i , j=1 x3i-2 = h x3i-2 - E f1 (x3i-2 ) + x3i-1 - E f2 (x3i-1 , x3i , x3i-2 ) + i = 1, 2, . . . , k. 3k j=1 a3i-1 j xj , (6) - b3i-1 , µx3i-2 , 45 А л е к с а н д р о в А. А., Г р а б к о Е. В. Входящие в эту систему числовые параметры aij и bj задаются соотношениями [19]: k d ωi 3j-l = b3i-l = a3i-r ∀l = 0, 1, 2; i = 1, 2, . . . k; ∆3-l (s)ds (7) k d K3-r 3-l (s, t)dt ds ∀r, l = 0, 1, 2; i, j = 1, 2, . . . k, ωi ωj где d есть площадь прямоугольника Ω. Для получения приближенного решения системы уравнений (1) можно использовать итерационный процесс (0) xi = 0 ∀i = 1, 2, . . . , 3k; (m+1) (m) , (m+1) (m) , βi (m) , αi x3i-2 = h γi x3i-1 = q αi (m+1) x3i = q βi (m) где величины αi (8) (m) , µh x3i-2 , (m) , µ · h x3i-2 , (m) , βi (m) (m) (m) , γi i = 1, 2, . . . , k; m = 0, 1, 2, . . . , , i = 1, 2, . . . , k, определяются равенствами 3k (m) γi = (m) x3i-2 -E (m) f1 x3i-2 (m) + - b3i-2 , a3i-2 j xj j=1 3k (m) αi = (m) x3i-1 -E (m) (m) (m) f2 x3i-1 , x3i , x3i-2 (m) + a3i-1 j xj - b3i-1 , (9) j=1 3k (m) βi (m) (m) (m) (m) (m) = x3i - E f2 x3i , x3i-1 , x3i-2 + a3i j xj - b3i . j=1 Будем полагать, что входящая в равенства (1) константа E удовлетворяет неравенству 3k Полученные результаты численного решения рассматриваемой пространственной контактной задачи представлены на рис. 1 и 2. Здесь приведены распределения нормальных σz и касательных τzx напряжений на поверхности полупространства вдоль поперечной оси симметрии контактной площадки, а также граничные контуры зон сцепления на этой площадке при μ = 0.2 и μ = 0.3. На рис. 1, а и 2, а изображено распределение нормальных контактных напряжений, на рис. 1, б и 2, б - распределение касательных контактных напряжений и на рис. 1, в и 2, в - граничные контуры зон сцепления (зоны сцепления представляют собой внутренность этих контуров). Здесь сплошной линией изображены распределения напряжений и граница зоны сцепления при отсутствии шероховатости, квадратиками - при наличии незначительной шероховатости и кольцами - при наличии существенной шероховатости. Приведенные на рис. 1 и 2 результаты свидетельствуют, что значение напряжений σz и τzx при отсутствии шероховатости несущественно отличаются от их значений при наличии шероховатости почти на всей площадке контакта (кроме участков, расположенных близ границ площадки контакта и близ границы между зоной сцепления и зоной проскальзывания). Следует отметить, что расхождение сравниваемых величин напряжений возрастает с ростом значения коэффициента трения μ. Из рис. 1, в и 2, в следует, что размеры зон сцепления, полученные при отсутствии шероховатости, могут очень сильно отличаться от этих размеров для случая наличия существенной шероховатости. Это означает, что качественные показатели контактного взаимодействия штампа и полупространства могут существенно зависеть от того, учитывается ли шероховатость поверхности полупространства или нет. Выводы. Полученные результаты свидетельствуют о том, что предложенный алгоритм решения статических задач о контакте упругих шероховатых тел при наличии кулонова трения между ними позволяет находить распределение контактных напряжений, а также конфигурацию зон проскальзывания и сцепления в задачах такого класса. Анализ полученных результатов решения конкретной контактной задачи свидетельствует о том, что учет шероховатости поверхности полупространства может приводить к существенному росту размеров зоны сцепления и к заметным изменениям в распределениях нормальных и касательных контактных напряжений по сравнению со случаем отсутствия шероховатости.

About the authors

Alexandr I Alexandrov

Zaporizhzhya National University

Email: heepper@gmail.com
66, Zhukovskogo st., Zaporizhzhya, 69600, Ukraine
(Cand. Tech. Sci.; heepper@gmail.com), Associate Professor, Dept. of Mathematical Analysis.

Elena V Grabko

Zaporizhzhya National University

Email: elenagrabko@rambler.ru
66, Zhukovskogo st., Zaporizhzhya, 69600, Ukraine
(elenagrabko@rambler.ru; Corresponding Author), Postgraduate Student, Dept. of Mathematical Analysis

References

  1. Александров В. М. Контактные задачи в трибологии / Механика и научно-технический прогресс. Т. 3, Механика деформируемого твердого тела. М.: Наука, 1988. С. 170- 180.
  2. Вовкушевский А. В. Вариационная постановка и методы решения контактной задачи с трением при учете шероховатости поверхностей // Изв. АН СССР. МТТ, 1991. № 3. С. 56-62.
  3. Горячева И. Г. Механика фрикционного взаимодействия. М.: Наука, 2001. 478 с.
  4. Александров В. М., Пожарский Д. А. Трехмерные контактные задачи при учете трения и нелинейной шероховатости // ПММ, 2004. Т. 68, № 3. С. 516-527.
  5. Pauk V., Zastrau B. W. Plane contact problems with partial slip for rough halfspace // J. Theor. Appl. Mech., 2004. vol. 42, no. 1. pp. 107-124
  6. Pauk V. Plane elastic contact involving friction and boundary roughness // EJPAU, 2006. vol. 9, no. 1, #31
  7. Дьяченко Н. Н., Шашкова Е. В. Контакт параболоидного штампа с упругим шероховатым полупространством в условиях частичного проскальзывания // Вiсник Запорiзького нацiонального унiверситету. Фiз.-мат. науки, 2010. No 2. С. 29-37.
  8. Дьяченко Н. Н., Манько Н. И.-В., Шашкова Е. В. Задача контакта квадратного в плане штампа с шероховатым полупространством в условиях частичного проскальзывания / Методи розв'язавння прикладних задач механiки деформiвного твердого тiла: Збiрник наукових праць Днiпропетровського нацiонального унiверситету, Вип. 13. Днiпропетровськ, 2012. С. 159-168.
  9. Грабко Е. В. Численное решение статической задачи о контакте упругих шероховатых тел при наличии кулонова трения / Проблеми обчислювальної механiки i мiцностi конструкцiй, Вип. 18. Днiпропетровськ: Лiра, 2012. С. 39-47.
  10. Галанов Б. А. Метод граничных уравнений типа Гаммерштейна для контактных задач теории упругости в случае неизвестных областей контакта // ПММ, 1985. Т. 49, № 5. С. 827-835.
  11. Александров А. И. Метод нелинейных граничных интегральных уравнений для решения пространственных контактных задач о взаимодействии упругих тел при наличии трения // Вiсник Днiпропетровського унiверситету. Сер. Механiка, 2010. Т. 18, № 14(1). С. 26-38
  12. Kalker J. J. A survey of the mechanics of contact between solid bodies // ZAAM, 1977. vol. 57, no. 5. pp. T3-T17. doi: 10.1002/zamm.19770570503.
  13. Александров А. И., Грабко Е. В. Теоремы существования решения для контактной задачи о взаимодействии упругих тел, имеющих шероховатые поверхности // Вiсник Запорiзького нацiонального унiверситету. Фiз.-мат. науки, 2010. No 2. С. 5-11.
  14. Love A. E. H. A treatise on the mathematical theory of elasticity. New York: Dover Publ., 1944. xviii+643 pp.
  15. Кравчук А. С. К постановке краевых задач теории упругости с трением на границе / Механика деформируемого твердого тела. Куйбышев: Куйбыш. ун-т, 1976. С. 102-105.
  16. Демкин Н. Б. Контактирование шероховатых тел. М.: Наука, 1970. 227 с.
  17. Писаренко Г. С., Яковлев А. П., Матвеев В. В. Справочник по сопротивлению материалов. Киев: Наукова думка, 1988. 736 с.
  18. Johnson K. L. Contact mechanics. Cambridge: Cambridge University Press. xi+452 pp. doi: 10.1017/cbo9781139171731
  19. Александров А. И., Грабко Е. В. Решение задач о контакте упругих шероховатых тел с использованием нелинейных интегральных уравнений / Методи розв'язавння прикладних задач механiки деформiвного твердого тiла: Збiрник наукових праць Днiпропетровського нацiонального унiверситету, Вип. 13. Днiпропетровськ, 2012. С. 14-21.
  20. Галин Л. А. Контактные задачи теории упругости. М.: Гостехиздат, 1953. 250 с.

Statistics

Views

Abstract - 29

PDF (Russian) - 8

Cited-By


Article Metrics

Metrics Loading ...

PlumX

Dimensions

Refbacks

  • There are currently no refbacks.

Copyright (c) 2014 Samara State Technical University

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies