Механика контакта

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

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

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

Оригинальная работа по контактной механике восходит к 1881 году, когда была опубликована статья «О контакте упругих тел» («Ueber die Berührung fester elastischer Körper» ) Генрих Герц. Герц пытался понять, как оптические свойства нескольких уложенных друг на друга линз могут измениться под силы, удерживая их вместе. Контактное напряжение Герца относится к локализованным напряжениям, которые развиваются, когда две криволинейные поверхности входят в контакт и слегка деформируются под воздействием приложенных нагрузок. Эта величина деформации зависит от модуля упругости материала, находящегося в контакте. Он дает контактное напряжение как функция нормальной контактной силы, радиусов кривизны обоих тел и модуля упругости обоих тел. Контактное напряжение Герца формирует основу для определения несущей способности и усталости долговечности подшипников, шестерен и других тел, где две поверхности соприкасаются.

Содержание
  • 1 История
  • 2 Классические решения для неклейкого эластичного контакта
    • 2.1 Контакт между сферой и полупространством
    • 2.2 Контакт между двумя сферами
    • 2.3 Контакт между двумя пересеченными цилиндры равного радиуса R {\ displaystyle R}R
    • 2.4 Контакт между жестким цилиндром и упругостью полупространством
    • 2.5 Контакт между жестким коническим индором и упругим полупространством
    • 2.6 Контакт между двумя цилиндрами с параллельными осями
    • 2.7 Контакт с подшипником
    • 2.8 Метод уменьшения размерности
  • 3 Теория Герца неклейкого упругого контакта
    • 3.1 Допущения теории Герца
    • 3.2 Методы аналитического решения
      • 3.2.1 Точечный контакт на (2D) полуплоскости
      • 3.2.2 Линейный контакт на (2D) полуплоскости
        • 3.2.2.1 Нормальная нагрузка в области (a, b) {\ displaystyle (a, b)}(a,b)
        • 3.2.2.2 Сдвиговая нагрузка по области (a, b) {\ displaystyle (a, b)}(a,b)
      • 3.2.3 Точечный контакт на a (3D) полупрос транство
    • 3.3 Техника численного решения iques
  • 4 Контакт между шероховатыми поверхностями
    • 4.1 Обзор модели GW
    • 4.2 Обзор модели GT
  • 5 Адгезионный контакт между упругими телами
    • 5.1 Модель жесткого контакта Брэдли
    • 5.2 Модель упругого контакта Джонсона-Кендалла-Робертса (JKR)
    • 5.3 Модель упругого контакта Дерягина-Мюллера-Топорова (DMT)
    • 5.4 Параметр Табора
    • 5.5 Модель упругого контакта Моугиса-Дагдейла
    • 5.6 Модель Carpick -Ogletree-Salmeron (COS)
    • 5.7 Влияние формы контакта
  • 6 См. Также
  • 7 Ссылки
  • 8 Внешние ссылки
История
Когда сфера прижимается к упругому материалу, площадь контакта увеличивается.

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

Только почти почти лет спустя Джонсон, Кендалл и Робертс нашли аналогичное решение для случая контакта клея. Эту теорию отвергли Борис Дерягин и его сотрудники, предложившие другую теорию адгезии в 1970-х годах. Модель Дерягина стала известна как модель ДМТ (после Дерягина, Мюллера и Топорова), а модель Джонсона и др. Модель стала известна как модель JKR (в честь Джонсона, Кендалла и Робертса) для адгезионного эластичного контакта. Этот отказ инструментом для разработки параметров Табора и более поздних параметров Моугиса, которые количественно определяют, какая модель контакта (из моделей JKR и DMT) лучше представляет адгезивный контакт для конкретных материалов.

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

Вклад Арчарда (1957) также следует представить при обсуждении новаторских работ в этой области. Арчард пришел к выводу, что даже для шероховатых упругих поверхностей площадь контакта пропорциональна нормальная силе. Дальнейшие важные идеи в этом направлении были предоставлены Гринвудом и Уильямсоном (1966 г.), Бушем (1975 г.) и Перссоном (2002 г.). Основные результаты этих работ заключаются в том, что истинная поверхность контакта в шероховатых материалах обычно пропорциональна нормальным результатам, в то время как параметры отдельных микроконтактов (т.е. давление, размер микроконтактов) слабо зависят от нагрузки..

Классические решения для неклейкого эластичного контакта

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

Контакт между сферой и полупространством

Контакт упругой сферы с упругим полупространством

Эластичная сфера с радиусом R {\ displaystyle R }Rвдавливает упругое полупространство, где общая деформация равна d {\ displaystyle d}d, вызывая площадь контакта радиусом

a = R d {\ displaystyle a = {\ sqrt {Rd}}}{\displaystyle a={\sqrt {Rd}}}

Приложенная сила F {\ displaystyle F}F связано со смещением d {\ displaystyle d}dна

F = 4 3 E ∗ R 1 2 d 3 2 {\ displaystyle F = {\ frac {4} {3}} E ^ {*} R ^ {\ frac {1} {2}} d ^ { \ frac {3} {2}}}{\ displaystyle F = {\ frac {4} {3}} E ^ {*} R ^ {\ frac {1} {2}} d ^ {\ frac { 3} {2}}}

где

1 E ∗ = 1 - ν 1 2 E 1 + 1 - ν 2 2 E 2 {\ displaystyle {\ frac {1} {E ^ {* }}} = {\ frac {1- \ nu _ {1} ^ {2}} {E_ {1}}} + {\ frac {1- \ nu _ {2} ^ {2}} {E_ {2 }}}}{\displaystyle {\frac {1}{E^{*}}}={\frac {1-\nu _{1}^{2}}{E_{1}}}+{\frac {1-\nu _{2}^{2}}{E_{2}}}}

и E 1 {\ displaystyle E_ {1}}E_{1},E 2 {\ displaystyle E_ {2}}E_ {2 } - это модули упругости и ν 1 {\ Displaystyle \ ню _ {1}}\ nu_1 ,ν 2 {\ Displaystyle \ п u _ {2}}\nu_2коэффициенты Пуассона, связанные с каждым телом.

Распределение нормального давления в области контакта как функция расстояния от центра круга:

p (r) = p 0 (1 - r 2 a 2) 1 2 {\ displaystyle p (r) = p_ {0 } \ left (1 - {\ frac {r ^ {2}} {a ^ {2}}} \ right) ^ {\ frac {1} {2}}}{\ displaystyle p (r) = p_ {0} \ left (1 - {\ frac {r ^ {2}} {a ^ {2}}} \ right) ^ {\ frac {1} {2}}}

где p 0 {\ displaystyle p_ {0}}p_{0}- максимальное контактное давление, определяемое по формуле

p 0 = 3 F 2 π a 2 = 1 π (6 FE ∗ 2 R 2) 1 3 {\ displaystyle p_ { 0} = {\ frac {3F} {2 \ pi a ^ {2}}} = {\ frac {1} {\ pi}} \ left ({\ frac {6F {E ^ {*}} ^ {2 }} {R ^ {2}}} \ right) ^ {\ frac {1} {3}}}{\displaystyle p_{0}={\frac {3F}{2\pi a^{2}}}={\frac {1}{\pi }}\left({\frac {6F{E^{*}}^{2}}{R^{2}}}\right)^{\frac {1}{3}}}

Радиус круга зависит от приложенной нагрузки F {\ displaystyle F}F по уравнению

a 3 = 3 FR 4 E ∗ {\ displaystyle a ^ {3} = {\ cfrac {3FR} {4E ^ {*}}}}{\displaystyle a^{3}={\cfrac {3FR}{4E^{*}}}}

Общая деформация d {\ displaystyle d}dсвязана с максимальным контактным давлением воздуха

d = a 2 R = (9 F 2 16 E ∗ 2 R) 1 3 {\ displaystyle d = {\ frac {a ^ {2}} {R }} = \ left ({\ frac {9F ^ {2}} {16 {E ^ {*}} ^ {2} R}} \ right) ^ {\ fr ac {1} {3}}}{\displaystyle d={\frac {a^{2}}{R}}=\left({\frac {9F^{2}}{16{E^{*}}^{2}R}}\right)^{\frac {1}{3}}}

Максимальное напряжение сдвига возникает внутри при z ≈ 0,49 a {\ displaystyle z \ приблизительно 0,49a}{\displaystyle z\approx 0.49a}для ν = 0, 33 {\ displaystyle \ nu = 0,33}\nu =0.33.

Контакт между двумя сферами

Контакт между двумя сферами. Контакт между двумя скрещенными цилиндрами равного радиуса.

Для контакта между двумя сферами радиуса R 1 {\ displaystyle R_ {1}}R_{1}и R 2 {\ displaystyle R_ {2}}R_{2}, область контакта представляет собой круг радиуса a {\ displaystyle a}a. Уравнения такие же, как для сферы, контактирующей с полуплоскостью, за исключением того, что эффективный радиус R {\ displaystyle R}Rопределяет как

1 R = 1 R 1 + 1 R 2 { \ displaystyle {\ frac {1} {R}} = {\ frac {1} {R_ {1}}} + {\ frac {1} {R_ {2}}}}{\displaystyle {\frac {1}{R}}={\frac {1}{R_{1}}}+{\frac {1}{R_{2}}}}

Контакт между двумя пересеченными цилиндры равного радиуса R {\ displaystyle R}R

Это эквивалентно контакту между сферой радиуса R {\ displaystyle R}Rи плоскостью.

Контакт между жестким цилиндром с концом концом и упругим полупространством

Контакт между жестким цилиндрическим индентором и упругим полупространством.

Если жесткий цилиндр вдавлен в упругую полупространство. пространство, он создает давление, описываемое формулой

p (r) = p 0 (1 - r 2 R 2) - 1 2 {\ displaystyle p (r) = p_ {0} \ left (1 - {\ frac {r ^ {2}} {R ^ {2}}} \ right) ^ {- {\ frac {1} {2}}}}{\displaystyle p(r)=p_{0}\left(1-{\frac {r^{2}}{R^{2}}}\right)^{-{\frac {1}{2}}}}

где R {\ displaystyle R}R- радиус цилиндра, а

p 0 = 1 π E ∗ d R {\ displaystyle p_ {0} = {\ frac {1} {\ pi}} E ^ {*} {\ frac {d} { R}}}{\ displaystyle p_ {0} = { \ frac {1} {\ pi}} E ^ {*} {\ frac {d} {R}}}

Связь между глубиной вдавливания и нормальной эксплуатации определяется выражением

F = 2 RE ∗ d {\ displaystyle F = 2RE ^ {*} d}{\ Displaystyle F = 2RE ^ {*} d}

Контакт между жестким коническим индентором и упругое полупространство

Контакт между жестким коническим индентором и упругим полупространством.

В случае вдавливания упругого полупространства модуля Юнга E {\ displaystyle E}Eс использованием жесткого конического индентора, глубины контакта ϵ {\ displaystyle \ epsilon}\ epsilon и радиуса контакта a {\ displaystyle a}aсвязаны с использованием

ϵ = загар ⁡ (θ) {\ displaystyle \ epsilon = a \ tan (\ theta)}{\displaystyle \epsilon =a\tan(\theta)}

с θ {\ displaystyle \ theta}\ theta определяет угол между плоскостью и боковой поверхностью конуса. Общая глубина отступа d {\ displaystyle d}dопределяется как:

d = π 2 ϵ {\ displaystyle d = {\ frac {\ pi} {2}} \ epsilon}{\displaystyle d={\frac {\pi }{2}}\epsilon }

Общая сила равна

F = π E 2 (1 - ν 2) a 2 загар ⁡ (θ) = 2 E π (1 - ν 2) d 2 загар ⁡ (θ) {\ displaystyle F = {\ frac {\ pi E} {2 \ left (1- \ nu ^ {2} \ right)}} a ^ {2} \ tan (\ theta) = {\ frac {2E} {\ pi \ left (1- \ nu ^ {2} \ right)}} {\ frac {d ^ {2}} {\ tan (\ theta)}}}{\displaystyle F={\frac {\pi E}{2\left(1-\nu ^{2}\right)}}a^{2}\tan(\theta)={\frac {2E}{\pi \left(1-\nu ^{2}\right)}}{\frac {d^{2}}{\tan(\theta)}}}

Распределение давления определяется как

p (r) = E d π a ( 1 - ν 2) пер ⁡ (ar + (ar) 2-1) = E d π a (1 - ν 2) cosh - 1 ⁡ (ar) {\ displaystyle p \ left (r \ right) = {\ frac {Ed} {\ pi a \ left (1- \ nu ^ {2} \ right)}} \ ln \ left ({\ frac {a} {r}} + {\ sqrt {\ left ({\ frac { a} {r}} \ right) ^ {2} -1}} \ right) = {\ frac {Ed} {\ pi a \ left (1- \ nu ^ {2} \ right)}} \ cosh ^ {- 1} \ left ({\ frac {a} {r}} \ right)}{\displaystyle p\left(r\right)={\frac {Ed}{\pi a\left(1-\nu ^{2}\right)}}\ln \left({\frac {a}{r}}+{\sqrt {\left({\frac {a}{r}}\right)^{2}-1}}\right)={\frac {Ed}{\pi a\left(1-\nu ^{2}\right)}}\cosh ^{-1}\left({\frac {a}{r}}\right)}

Напряжение имеет логарифмическую сингулярность на конце конуса.

Контакт между двумя цилиндрами с параллельными осями

Контакт между двумя цилиндрами с параллельными осями

Контакт между двумя цилиндрами с параллельными осями сила линейно пропорциональна длине цилиндров L и углублению глубина d:

F ≈ π 4 E ∗ L d {\ displaystyle F \ приблизительно {\ frac {\ pi} {4}} E ^ {*} Ld}{\displaystyle F\approx {\frac {\pi }{4}}E^{*}Ld}

Радиусы кривизны полностью отсутствуют в этом отношении. Радиус контакта описывается обычным воздействием

a = R d {\ displaystyle a = {\ sqrt {Rd}}}{\displaystyle a={\sqrt {Rd}}}

с

1 R = 1 R 1 + 1 R 2 {\ displaystyle {\ displaystyle {\ frac {1} {R}} = {\ frac {1} {R_ {1}}} + {\ frac {1} {R_ {2}}}}{\displaystyle {\frac {1}{R}}={\frac {1}{R_{1}}}+{\frac {1}{R_{2}}}}

как при контакте двух сфер. Максимальное давление равно

p 0 = (E ∗ F π LR) 1 2 {\ displaystyle p_ {0} = \ left ({\ frac {E ^ {*} F} {\ pi LR}} \ справа) ^ {\ frac {1} {2}}}{\ displaystyle p_ {0} = \ left ({\ frac {E ^ {*} F} {\ pi LR}} \ right) ^ {\ frac {1 }{2}}}

Контакт подшипника

Контакт в случае подшипников часто представляет собой контакт между выпуклой поверхностью (охватываемый цилиндр или сфера) и вогнутой поверхностью (охватывающий цилиндр или сфера: отверстие или полусферическая чашка ).

Метод уменьшения размерности

Контакт между сферой и упругим полупространством и одномерная замененная модель.

Некоторые контактные задачи могут быть решены с помощью метода уменьшения размераности (MDR). В этом методе исходной трехмерной системы заменяется контактом тела с линейно-упругим или вязкоупругим основанием (см. Рис.). Свойства одномерных систем в точности совпадают со свойствами исходной трехмерной системы. MDR основан на решении осесимметричных контактных задач, впервые полученных Людвигом Фёпплом (1941) и Герхардом Шубертом (1942)

Однако требуется точных аналитических результатов, чтобы контактная задача была осесимметричной, контакты были компактным.

Теория Герца наличие неадгезионного упругого контакта

Классическая теория контакта фокусировалась в первую очередь на неклейком контакте, когда в зоне контакта не допускается силы натяжения, т. Е. Контактирующие тела могут отделяться без сил сцепления. Несколько аналитических и численных подходов были использованы для решения контактных задач, удовлетворяющих условию отсутствия прилипания. Сложные силы и моменты передаются между телами в местах их соприкосновения, поэтому проблемы в механике контакта могут стать довольно сложными. Кроме того, контактные напряжения обычно нелинейно зависит от деформации. Для упрощения процедур решения обычно используется система отсчета, в которой используются объекты (возможно, движущиеся относительно друг друга) статичны. Они взаимодействуют посредством поверхностного натяжения (или давления / напряжения) на их границе раздела.

В качестве примера рассмотрим два объекта, которые встречаются на некоторой поверхности S {\ displaystyle S}Sв (x {\ displaystyle x}x,y {\ displaystyle y}y) с осью z {\ displaystyle z}z, перпендикулярной поверхности. Одно из тел выполняет нормально направленное давление распределение pz = p (x, y) = qz (x, y) {\ displaystyle p_ {z} = p (x, y) = q_ {z } (x, y)}p_{z}=p(x,y)=q_{z}(x,y)и распределение сцепления в плоскости на поверхности qx = qx (x, y) {\ displaystyle q_ {x} = q_ {x} ( x, y)}q_{x}=q_{x}(x,y)и qy = qy (x, y) {\ displaystyle q_ {y} = q_ {y} (x, y)}q_{y}=q_{y}(x,y)по региону S {\ Displaystyle S}S. В терминах баланса ньютоновских сил, силы:

P z = ∫ S p (x, y) d A; Q x = ∫ S q x (x, y) d A; Q Y знак равно ∫ S qy (x, y) d A {\ displaystyle P_ {z} = \ int _ {S} p (x, y) ~ \ mathrm {d} A ~; ~~ Q_ {x} = \ int _ {S} q_ {x} (x, y) ~ \ mathrm {d} A ~; ~~ Q_ {y} = \ int _ {S} q_ {y} (x, y) ~ \ mathrm {d} A}P_{z}=\int _{S}p(x,y)~{\mathrm {d}}A~;~~Q_{x}=\int _{S}q_{x}(x,y)~{\mathrm {d}}A~;~~Q_{y}=\int _{S}q_{y}(x,y)~{\mathrm {d}}A

должно быть равным и противоположным силойм, установленным в другом теле. Моменты, соответствующие этой силе:

M x = ∫ S y q z (x, y) d A; M y = ∫ S - x q z (x, y) d A; Знак MZ равно ∫ S [xqy (x, y) - yqx (x, y)] d A {\ displaystyle M_ {x} = \ int _ {S} y ~ q_ {z} (x, y) ~ \ mathrm {d} A ~; ~~ M_ {y} = \ int _ {S} -x ~ q_ {z} (x, y) ~ \ mathrm {d} A ~; ~~ M_ {z} = \ int _ {S} [x ~ q_ {y} (x, y) -y ~ q_ {x} (x, y)] ~ \ mathrm {d} A}{\ displaystyle M_ {x} = \ int _ {S} y ~ q_ {z} (x, y) ~ \ mathrm {d} A ~; ~~ M_ {y } = \ int _ {S} -x ~ q_ {z} (x, y) ~ \ mathrm {d} A ~; ~~ M_ {z} = \ int _ {S} [x ~ q_ {y} ( x, y) -y ~ q_ {x} (x, y)] ~ \ mathrm {d} A}

также требуется для отмены между телами, чтобы они кинематически неподвижен.

Допущения теории Герца

При формировании решений контактных задач Герца делаются следующие допущения:

  • образования малы и находятся в пределах упругости.
  • Поверхности сплошные и не соответствуют друг другу (это означает, что площадь контакта намного меньше характерных размеров контактирующих тел).
  • Каждое тело можно рассматривать как упругое полупространство.
  • Поверхности не имеют трения.

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

методами аналитического решения

Контакт между двумя сферами.

Аналитические методы решения проблемы неклейкого контакта можно разделить на два типа в зависимости от геометрии области контакта. Соответствующий контакт - это контакт, при котором два тела соприкасаются в нескольких точках до того, как происходит какая-либо деформация (т.е. они просто «подходят друг к другу»). несоответствующий контакт - это контакт, в котором достаточно различны, поэтому при нулевой нагрузке они касаются только в точке (или, возможно, вдоль линии). В случае несоответствия площади контакта по сравнению с размерами объектов, и напряжения сконцентрированы в этой области. Такой контакт называется концентрированным, иначе - разноплановым.

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

Точечный контакт на (2D) полуплоскости

Схема нагружения на плоскости силой P в точке (0, 0).

Отправной точкой для решения проблем контакта является понять эффект «точечной нагрузки», приложенной к изотропной, однородной и линейной упругой полуплоскости, показанной на рисунке справа. Проблема может заключаться в плоском напряжении или плоской деформации. Это краевая задача линейной упругости с учетом тяги граничных условий :

σ x z (x, 0) = 0; σ Z (Икс, Z) знак равно - п δ (Икс, Z) {\ Displaystyle \ sigma _ {xz} (х, 0) = 0 ~; ~~ \ sigma _ {z} (x, z) = - P \ delta (x, z)}{\displaystyle \sigma _{xz}(x,0)=0~;~~\sigma _{z}(x,z)=-P\delta (x,z)}

где δ (x, z) {\ displaystyle \ delta (x, z)}{\displaystyle \delta (x,z)}- это дельта-функция Дирака. Граничные условия утверждают, что на поверхности нет касательных напряжений и дает в точке (0, 0) приложена сингулярная нормальная сила P. Применение этих условий к основным уравнениям упругости результат

σ xx = - 2 P π x 2 z (x 2 + z 2) 2 σ zz знак равно - 2 п π z 3 (x 2 + z 2) 2 σ xz = - 2 п π xz 2 (x 2 + z 2) 2 {\ displaystyle {\ begin {align} \ sigma _ {xx} = - {\ frac {2P} {\ pi}} {\ frac {x ^ {2} z} {\ left (x ^ {2} + z ^ {2} \ right) ^ {2} }} \\\ sigma _ {zz} = - {\ frac {2P} {\ pi}} {\ frac {z ^ {3}} {\ left (x ^ {2} + z ^ {2} \ справа) ^ {2}}} \\\ sigma _ {xz} = - {\ frac {2P} {\ pi}} {\ frac {xz ^ {2}} {\ left (x ^ {2} + z ^ {2} \ right) ^ {2}}} \ end {выравнивается}}}{\displaystyle {\begin{aligned}\sigma _{xx}=-{\frac {2P}{\pi }}{\frac {x^{2}z}{\left(x^{2}+z^{2}\right)^{2}}}\\\sigma _{zz}=-{\frac {2P}{\pi }}{\frac { z^{3}}{\left(x^{2}+z^{2}\right)^{2}}}\\\sigma _{xz}=-{\frac {2P}{\pi }}{\frac {xz^{2}}{\left(x^{2}+z^{2}\right)^{2}}}\end{aligned}}}

для некоторой точки, (x, y) {\ displaystyle (x, y)}(x,y), в полуплоскости. Круг, показанный на рисунке, указывает на поверхность, на максимальное напряжение сдвига постоянно. Из этого поля напряжений и таким образом быть связаны компоненты деформации, с ущербом всех материальных точек.

Линейный контакт на (2D) полуплоскости

Нормальная нагрузка в области (a, b) {\ displaystyle (a, b)}(a,b)

Предположим, скорее чем точечная нагрузка P {\ displaystyle P}P, вместо этого к применяемой распределенной нагрузке p (x) {\ displaystyle p (x)}p(x), в диапазоне a < x < b {\displaystyle aa<x<b. Принцип линейной суперпозиции может использоваться для определения результирующего поля напряжений как решения интегральных уравнений:

σ xx = - 2 z π ∫ abp (x ′) (x - x ′) 2 dx ′ [(x - х ') 2 + z 2] 2; σ zz = - 2 z 3 π ∫ abp (x ′) dx ′ [(x - x ′) 2 + z 2] 2 σ xz = - 2 z 2 π ∫ abp (x ′) (x - x ′) dx ′ [(Икс - Икс) 2 + Z 2] 2 {\ Displaystyle {\ begin {Выровнено} \ sigma _ {xx} = - {\ frac {2z} {\ pi}} \ int _ {a} ^ { b} {\ frac {p \ left (x '\ right) \ left (xx' \ right) ^ {2} \, dx '} {\ left [\ left (xx' \ right) ^ {2} + z ^ {2} \ right] ^ {2}}} ~; ~~ \ sigma _ {zz} = - {\ frac {2z ^ {3}} {\ pi}} \ int _ {a} ^ {b} {\ frac {p \ left (x '\ right) \, dx '} {\ left [\ left (xx' \ right) ^ {2} + z ^ {2} \ right] ^ {2}}} \\ [3pt] \ sigma _ {xz} = - {\ frac {2z ^ {2}} {\ pi}} \ int _ {a} ^ {b} {\ frac {p \ left (x '\ right) \ left (xx' \ right) \, dx '} { \ left [\ left (xx '\ right) ^ {2} + z ^ {2} \ right] ^ {2}}} \ end {align}}}{\displaystyle {\begin{aligned}\sigma _{xx}=-{\frac {2z}{\pi }}\int _{a}^{b}{\frac {p\left(x'\right)\left(x-x'\right)^{2}\,dx'}{\left[\left(x-x'\right)^{2}+z^{2}\right]^{2}}}~;~~\sigma _{zz}=-{\frac {2z^{3}}{\pi }}\int _{a}^{b}{\frac {p\left(x'\right)\,dx'}{\left[\left(x-x'\right)^{2}+z^{2}\right]^{2}}}\\[3pt]\sigma _{xz}=-{\frac {2z^{2}}{\pi }}\int _{a}^{b}{\frac {p\left(x'\right)\left(x-x'\right)\,dx'}{\left[\left(x-x'\right)^{2}+z^{2}\right]^{2}}}\end{aligned}}}
Сдвиговая нагрузка по области (a, б) {\ displaystyle (a, b)}(a,b)

Тот же принцип работы для нагрузки на поверхность в плоскости поверхности. Эти виды сцепления могут возникнуть в результате трения. Решение аналогично приведенному выше (как для единичных нагрузок Q {\ displaystyle Q}Q, так и для распределенных нагрузок q (x) {\ displaystyle q (x)}q(x)), но с небольшими изменениями:

σ xx = - 2 π ∫ abq (x ′) (x - x ′) 3 dx ′ [(x - x ′) 2 + z 2] 2; σ zz = - 2 z 2 π ∫ abq (x ′) (x - x ′) dx ′ [(x - x ′) 2 + z 2] 2 σ xz = - 2 z π ∫ abq (x ′) (x - x ′) 2 dx ′ [(x - x ′) 2 + z 2] 2 {\ displaystyle {\ begin {align} \ sigma _ {xx} = - {\ frac {2} {\ pi}} \ int _ {a} ^ {b} {\ frac {q \ left (x '\ right) \ left (xx' \ right) ^ {3} \, dx '} {\ left [\ left (x- x') \ right) ^ {2} + z ^ {2} \ right] ^ {2}}} ~; ~~ \ sigma _ {zz} = - {\ frac {2z ^ {2}} {\ pi}} \ int _ {a} ^ {b} {\ frac {q \ left (x '\ right) \ left (xx '\ right) \, dx'} {\ left [\ left (xx '\ right) ^ {2} + z ^ {2} \ right] ^ {2}}} \\ [3pt] \ sigma _ {xz} = - {\ frac {2z} {\ pi}} \ int _ {a} ^ {b} {\ frac {q \ left (x '\ right) \ left (xx' \ right) ^ { 2} \, dx '} {\ left [\ left (x-x') \ right) ^ {2} + z ^ {2} \ right] ^ {2}}} \ end {align}}}{\displaystyle {\begin{aligned}\sigma _{xx}=-{\frac {2}{\pi }}\int _{a}^{b}{\frac {q\left(x'\right)\left(x-x'\right)^{3}\,dx'}{\left[\left(x-x'\right)^{2}+z^{2}\right]^{2}}}~;~~\sigma _{zz}=-{\frac {2z^{2}}{\pi }}\int _{a}^{b}{\frac {q\left(x'\right)\left(x-x'\right)\,dx'}{\left[\left(x-x'\right)^{2}+z^{2}\right]^{2}}}\\[3pt]\sigma _{xz}=-{\frac {2z}{\pi }}\int _{a}^{b}{\frac {q\left(x'\right)\left(x-x'\right)^{2}\,dx'}{\left[\left(x-x'\right)^{2}+z^{2}\right]^{2}}}\end{aligned}}}

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

Точечный контакт в (3D) полупространстве

Аналогично решение Flamant для 2D полуплоскости, фундаментальные решения известны и для линейно упругого 3D полупространства. Они были обнаружены Буссинеском для концентрированной нормальной нагрузки и Черрути для тангенциальной нагрузки. См. Раздел по этому вопросу в Линейная упругость.

Методы численного

При использовании схем численного решения для контактных проблем не нужно делать различия между и несоответствующим контактом. Эти методы не основаны на предположениях в процессе решения, поскольку они основаны исключительно на общих основаниях. Помимо стандартных правил, описывающих деформацию и движение тел, можно указать два дополнительных неравенства. Первый просто ограничивает движение и деформацию, предполагаемая, что проникновение. Следовательно, зазор h {\ displaystyle h}hмежду двумя телами может быть только положительным или нулевым

h ≥ 0 {\ displaystyle h \ geq 0}{\displaystyle h\geq 0}

где h = 0 {\ displaystyle h = 0}h=0обозначает контакт. Второе допущение в механике связано с тем, что в зоне контакта не должно возникать силы натяжения (контактирующие тела можно поднимать без сил сцепления). Это приводит к неравенству. Он разработан для нормального напряжения σ n = t ⋅ n {\ displaystyle \ sigma _ {n} = \ mathbf {t} \ cdot \ mathbf {n}}{\displaystyle \sigma _{n}=\mathbf {t} \cdot \mathbf {n} }.

в местах, где есть контакт между на поверхности зазор равенство нулю, то есть h = 0 {\ displaystyle h = 0}h=0, и там нормальное напряжение отличается от нуля, действительно, σ n < 0 {\displaystyle \sigma _{n}<0}{\displaystyle \sigma _{n}<0}. В местах, где поверхности не соприкасаются, нормальное напряжение равно нулю; σ N = 0 {\ Displaystyle \ sigma _ {n} = 0}{\displaystyle \sigma _{n}=0}, разрыв положительный; т.е. h>0 {\ displaystyle h>0}{\displaystyle h>0} . Этот тип формулировки комплементарности может быть выражен в так называемой форме Куна - Такера, а именно

h ≥ 0, σ n ≤ 0, σ NH знак равно 0, {\ Displaystyle ч \ GEQ 0 \,, \ quad \ sigma _ {n} \ leq 0 \,, \ quad \ sigma _ {n} \, h = 0 \,.}{\displaystyle h\geq 0\,,\quad \sigma _{n}\leq 0\,,\quad \sigma _{n}\,h=0\,.}

Эти условия Математическая формулировка зазора зависит от кинематики лежащей в основе теории твердого тела (например, линейное или нелинейное твердое тело в двух- или трех измерениях, балка или оболочка модель). Пересчитав нормальное напряжение σ n {\ displaystyle \ sigma _ {n}}\sigma_nс точки зрения контактного давления, p {\ displaystyle p}p; т. Е. p = - σ n {\ displaystyle p = - \ sigma _ {n}}{\displaystyle p=-\sigma _{n}}проблема Куна-Такера может быть переформулиро вана как в стандартной комплементарности форма, т.е.

h ≥ 0, p ≥ 0, ph = 0. {\ di splaystyle h \ geq 0 \,, \ quad p \ geq 0 \,, \ quad p \, h = 0 \,.}{\displaystyle h\geq 0\,,\quad p\geq 0\,,\quad p\,h=0\,.}В линейно-упругом случае зазор можно сформулировать как h = h 0 + g + u, {\ displaystyle {h} = h_ {0} + {g} + u,}{\ displaystyle {h} = h_ {0} + {g} + u,} где h 0 {\ displaystyle h_ {0}}{\displaystyle h_{0}}- жесткий разделение тела, g {\ displaystyle g}g- это геометрия / топография контакта (цилиндр и шероховатость), а u {\ displaystyle u}{\displaystyle u}- эластичная деформация / прогиб. Если контактирующие тела аппроксимированы линейными упругими полупространствами, решение интегрального уравнения Буссинеска-Черрути может быть применено для выражения деформации (u {\ displaystyle u}u ) как функции контактного давления (п {\ displaystyle p }p); т. е. U знак равно ∫ ∞ ∞ К (х - s) п (s) ds, {\ displaystyle u = \ int _ {\ infty} ^ {\ infty} K (xs) p (s) ds,}{\displaystyle u=\int_{\infty }^{\infty }K(x-s)p(s)ds,}где K (x - s) = 2 π E ∗ ln ⁡ | х - с | {\ displaystyle K (xs) = {\ frac {2} {\ pi E ^ {*}}} \ ln | xs |}{\displaystyle K(x-s)={\frac {2}{\pi E^{*}}}\ln |x-s|}для линейной загрузки эластичного полупространства и K (x - s) Знак равно 1 π E ∗ 1 (Икс 1 - s 1) 2 + (Икс 2 - s 2) 2 {\ Displaystyle K ( xs) = {\ frac {1} {\ pi E ^ {*}}} {\ гидроразрыв {1} {\ sqrt {\ left (x_ {1} -s_ {1} \ right) ^ {2} + \ left (x_ {2} -s_ {2} \ right) ^ {2}}}}}{\displaystyle K(x-s)={\frac {1}{\pi E^{*}}}{\frac {1}{\sqrt {\left(x_{1}-s_{1}\right)^{2}+\left(x_{2}-s_{2}\right)^{2}}}}}для точечного нагружения упругого полупространства.

После дискретизации задачи линейной упругой контактной механики может быть сформулирована в стандартной форме задачи линейной дополнительности (LCP).

час = час 0 + g + C p час ⋅ п знак равно 0, п ≥ 0, час ≥ 0, {\ Displaystyle {\ begin {выровнено} \ mathbf {h} ​​= \ mathbf {h} _ {0} + \ mathbf {g} + \ mathbf {Cp}, \\\ mathbf {h} \ cdot \ mathbf {p} = 0, \, \, \, \ mathbf {p} \ geq 0, \, \, \, \ mathbf {h} \ geq 0, \\\ end {align}}}{\displaystyle {\begin{aligned}\mathbf {h} =\mathbf {h} _{0}+\mathbf {g} +\mathbf {Cp},\\\mathbf {h} \cdot \mathbf {p} =0,\,\,\,\mathbf {p} \geq 0,\,\,\,\mathbf {h} \geq 0,\\\end{aligned}}}

где C {\ displaystyle \ mathbf {C}}\ mathbf {C} - матрица, элементы, которые связывают контактное давление и деформация. Представленная выше строгая формулировка задачи LCP CM позволяет напрямую применять хорошо зарекомендовавшие себя численные методы решения, такие как алгоритм поворота Лемке. Преимущество алгоритма Лемке в том, что он находит численно точное решение за конечное число итераций. Реализация MATLAB, представленная Алмквистом и др., является одним из примеров, который можно использовать для численного решения проблемы. Кроме того, пример кода для LCP-решения задачи механики линейного упругого контакта 2D также опубликован при обмене файлов MATLAB Альмквистом и др.

Контакт между шероховатыми поверхностями

Когда два тела с шероховатой поверхностью прижимаются друг к другу, истинная площадь контакта, образованная между двумя телами, A {\ displaystyle A}A , намного меньше, чем кажущаяся или номинальная площадь контакта A 0 {\ displaystyle A_ {0}}A_{0}. Механика контакта с шероховатыми поверхностями обсуждается с точки зрения механики контакта и статических фрикционных взаимодействий. Природные и инженерные поверхности обычно демонстрируют особенности шероховатости, известные как неровности, в широком диапазоне масштабов длины до молекулярного уровня, при этом поверхностные структуры демонстрируют само сродство, также известное как фрактальность поверхности. Признано, что самоаффинная поверхность является средством линейного масштабирования истинной площади контакта с приложенным давлением. Допуская модель сдвига сварных контактов в трибологических взаимодействий, эту повсеместно наблюдаемую линейность между площадью контакта и давлением также можно рассматривать как источник линейности зависимости между статическим трением и приложенной нормальной силой.

При контакте между «случайной шероховатой» поверхностью и упругим полупространством истинная площадь взаимодействия с нормальной силой F {\ displaystyle F}F by

A = κ E ∗ h ′ F {\ displaystyle A = { \ frac {\ kappa} {E ^ {*} h '}} F}A={\frac {\kappa }{E^{*}h'}}F

с h ′ {\ displaystyle h'}h'равным среднеквадратичному (известное также как среднее квадратичное) уклона поверхности и κ ≈ 2 {\ displaystyle \ kappa \ приблизительно 2}\kappa \approx 2. Среднее давление на истинной контактной поверхности

pav = FA ≈ 1 2 E ∗ h ′ {\ displaystyle p _ {\ mathrm {av}} = {\ frac {F} {A}} \ приблизительно {\ frac {1} {2}} E ^ {*} h '}p_{{{\mathrm {av}}}}={\frac {F}{A}}\approx {\frac {1}{2}}E^{*}h'

можно разумно оценить как половину эффективного модуля упругости E ∗ {\ displaystyle E ^ {*}}E^{*}, умнож на средний квадрат уклона поверхности h ′ {\ displaystyle h '}h'.

Обзор модели GW

Гринвуд и Уильямсон в 1966 г. (GW) предложили теорию механики упругого контакта грубых поверхностей, которые сегодня составляют теорий трибологии (трение, адгезия, тепловая и электрическая проводимость, износ и т. Д.). Они рассмотрели контакт между гладкой жесткой плоскостью и номинально плоской деформируемой шероховатой поверхностью, покрытой круглыми выступами на вершине того же радиуса R. Их теория предполагает, что деформация каждой выступающей части не зависит от деформации ее соседей и описывается моделью Герца.. Высота выступов распределена случайным образом. Вероятность того, что высота выступов находится между z {\ displaystyle z}zи z + dz {\ displaystyle z + dz}{\displaystyle z+dz}, равна ϕ (z) dz {\ displaystyle \ phi (z) dz}{\displaystyle \phi (z)dz}. Авторы рассчитали количество пятен контакта n, общую площадь контакта A r {\ displaystyle A_ {r}}A_{r}и общую нагрузку P в общем случае. Они представили эти формулы в двух формах: в базовой и с использованием стандартизованных переменных. Если предположить, что N выступов покрывают шероховатую поверхность, то ожидаемое количество контактов будет

n = N ∫ d ∞ ϕ (z) dz {\ displaystyle n = N \ int _ {d} ^ {\ infty} \ phi (z) dz}{\displaystyle n=N\int _{d}^{\infty }\phi (z)dz}

Ожидаемую общую площадь контакта можно рассчитать по формуле

A a = N π R ∫ d ∞ (z - d) ϕ (z) dz {\ displaystyle A_ {a} = N \ pi R \ int _ {d} ^ {\ infty} (zd) \ phi (z) dz}{\displaystyle A_{a}=N\pi R\int _{d}^{\infty }(z-d)\phi (z)dz}

, а ожидаемая общая сила определяется как

P = 4 3 NE r R ∫ d ∞ (z - d) 3 2 ϕ (z) dz {\ displaystyle P = {\ frac {4} {3}} NE_ {r} {\ sqrt {R}} \ int _ {d} ^ {\ infty} ( zd) ^ {\ frac {3} {2}} \ phi (z) dz}{\displaystyle P={\frac {4}{3}}NE_{r}{\sqrt {R}}\int _{d}^{\infty }(z-d)^{\frac {3}{2}}\phi (z)dz}

где:

R, радиус кривизны микровыступа,
z, высота измеренного микровыступа от линии профиля,
d, закрыть поверхность,
E r = (1 - ν 1 2 E 1 + 1 - ν 2 2 E 2) - 1 {\ displaystyle E_ {r} = \ left ({\ frac {1- \ nu _ {1} ^ {2}} {E_ {1}}} + {\ frac {1- \ nu _ {2} ^ {2}} {E_ {2}) }} \ right) ^ {- 1}}{\displaystyle E_{r}=\left({\frac {1-\nu _{1}^{2}}{E_{1}}}+{\frac {1-\nu _{2}^{2}}{E_{2}}}\right)^{-1}}, модуль упругости составного Юнга,
E i {\ displaystyle E_ {i}}E_{i}, модуль упругости поверхности лицо,
ν i {\ displaystyle \ nu _ {i}}\nu _{i}, поверхностные коэффициенты Пуассона.

Они ввели стандартизированное разделение h = d / σ {\ displaystyle h = d / \ sigma}{\displaystyle h=d/\sigma }и стандартизированное распределение высот ϕ ∗ (s) {\ displaystyle \ phi ^ {*} (s)}{\displaystyle \phi ^{*}(s)}, стандартное отклонение которого равно единице. Ниже представлены формулы в стандартизованном виде.

F n (h) = ∫ h ∞ (s - h) n ϕ ∗ (s) dsn = η A n F 0 (h) A a = π η AR σ F 1 (h) P = 4 3 η AE р р σ 3 2 F 3 2 (час) {\ Displaystyle {\ begin {выровнено} F_ {n} (h) = \ int _ {h} ^ {\ infty} (sh) ^ {n} \ phi ^ {*} (s) ds \\ n = \ eta A_ {n} F_ {0} (h) \\ A_ {a} = \ pi \ eta AR \ sigma F_ {1} (h) \\ P = {\ frac {4} {3}} \ eta AE_ {r} {\ sqrt {R}} \ sigma ^ {\ frac {3} {2}} F _ {\ frac {3} {2}} (h) \ end {align}}}{\displaystyle {\begin{aligned}F_{n}(h)=\int _{h}^{\infty }(s-h)^{n}\phi ^{*}(s)ds\\n=\eta A_{n}F_{0}(h)\\A_{a}=\pi \eta AR\sigma F_{1}(h)\\P={\frac {4}{3}}\eta AE_{r}{\sqrt {R}}\sigma ^{\frac {3}{2}}F_{\frac {3}{2}}(h)\end{aligned}}}

где:

d - расстояние,
A {\ displaystyle A}A - номинальная площадь контакта,
η {\ displaystyle \ eta}\eta - поверхностная плотность неровностей,
E ∗ {\ displaystyle E ^ {*}}E^{*}- эффективный модуль Юнга.

В последнее время точное приближение к A r {\ displaystyle A_ {r}}A_{r}и P {\ displaystyle P}Pбыли опубликованы Jedynak. Они задаются следующими рациональными формулами, которые являются очень точными приближениями к интегралам F n (h) {\ displaystyle F_ {n} (h)}{\ displaystyle F_ {n} (h)} . Они рассчитаны для гауссова распределения неровностей

F n (h) = a 0 + a 1 h + a 2 h 2 + a 3 h 3 1 + b 1 h + b 2 h 2 + b 3 h 3 + б 4 час 4 + b 5 час 5 + b 6 час 6 ехр ⁡ (- час 2 2) {\ displaystyle F_ {n} (h) = {\ frac {a_ {0} + a_ {1} h + a_ { 2} h ^ {2} + a_ {3} h ^ {3}} {1 + b_ {1} h + b_ {2} h ^ {2} + b_ {3} h ^ {3} + b_ {4 } h ^ {4} + b_ {5} h ^ {5} + b_ {6} h ^ {6}}} \ exp \ left (- {\ frac {h ^ {2}} {2}} \ right)}{\displaystyle F_{n}(h)={\frac {a_{0}+a_{1}h+a_{2}h^{2}+a_{3}h^{3}}{1+b_{1}h+b_{2}h^{2}+b_{3}h^{3}+b_{4}h^{4}+b_{5}h^{5}+b_{6}h^{6}}}\exp \left(-{\frac {h^{2}}{2}}\right)}

Для F 1 (h) {\ displaystyle F_ {1} (h)}{\displaystyle F_{1}(h)}коэффициенты

[a 0, a 1, a 2, a 3] = [0,398942280401, 0,159773702775, 0,0389687688311, 0,00364356495452] [b 1, b 2, b 3, b 4, b 5, b 6] = [1,653807476138, 1,170419428529, 0,448892964428, 0,0951971709160, 0,003692964428, 0,0951971709160, 0,00369389 - 6,34280 \ displaystyle {\ begin {align} [a_ {0}, a_ {1}, a_ {2}, a_ {3}] = [0.398942280401,0.159773702775,0.0389687688311,0.00364356495452] \\ [b_ {1}, b_ {2 }, b_ {3}, b_ {4}, b_ {5}, b_ {6}] = \ left [1.653807476138,1.170419428529,0.448892964428,0.0951971709160,0.00931642803836, -6.383774657279 \ times 10 ^ {- 6} \ right] \ end {ali gned}}}{\displaystyle {\begin{aligned} [a_{0},a_{1},a_{2},a_{3}]=[0.398942280401,0.159773702775,0.0389687688311,0.00364356495452]\\ [b_{1},b_{2},b_{3},b_{4},b_{5},b_{6}]=\left[1.653807476138,1.170419428529,0.448892964428,0.0951971709160,0.00931642803836,-6.383774657279\times 10^{-6}\right]\end{aligned}}}

Максимальная относительная ошибка составляет 9.93 × 10-8% {\ displaystyle 9.93 \ times 10 ^ {- 8} \%}{\displaystyle 9.93\times 10^{-8}\%}.

Для F 3 2 (h) { \ displaystyle F _ {\ frac {3} {2}} (h)}{\displaystyle F_{\frac {3}{2}}(h)}коэффициенты:

[a 0, a 1, a 2, a 3] = [0.430019993662, 0.101979509447, 0.0229040629580, 0.000688602924 ] [b 1, b 2, b 3, b 4, b 5, b 6] = [1,671117125984, 1,199586555505, 0,46936532151, 0,102632881122, 0,010686348714, 0,0000517200271] {\ displaystyle { \ begin {align} [a_ {0}, a_ {1}, a_ {2}, a_ {3}] = [0.430019993662,0.101979509447,0.0229040629580,0.000688602924] \\ [b_ {1}, b_ {2}, b_ {3}, b_ {4}, b_ {5}, b_ {6}] = [1.671117125984,1.199586555505, 0.46936532151,0.102632881122,0.010686348714,0.0000517200271] \ end {align}}{\displaystyle {\begin{aligned} [a_{0},a_{1},a_{2},a_{3}]=[0.430019993662,0.101979509447,0.0229040629580,0.000688602924]\\ [b_{1},b_{2},b_{3},b_{4},b_{5},b_{6}]=[1.671117125984,1.199586555505,0.46936532151,0.102632881122,0.010686348714,0.0000517200271]\end{aligned}}}

Максимальная относительная ошибка 1.91 × 10-7% {\ displaystyle 1.91 \ times 10 ^ {- 7} \%}{\displaystyle 1.91\times 10^{-7}\%}. В документе также используются точные выражения для F n (h) {\ displaystyle F_ {n} (h)}{\ displaystyle F_ {n} (h)}

F 1 (h) = 1 2 π exp ⁡ (- 1 2 h 2) - 1 2 h erfc ⁡ (h 2) F 3 2 (h) = 1 4 π exp ⁡ (- h 2 4) h ((h 2 + 1) K 1 4 (h 2 4) - h 2 K 3 4 (h 2 4)) {\ displaystyle {\ begin {align} F_ {1} (h) = {\ frac {1} {\ sqrt {2 \ pi}}} \ exp \ left (- {\ frac {1} { 2}} h ^ {2} \ right) - {\ frac {1} {2}} h \, \ operatorname {erfc} \ left ({\ frac {h} {\ sqrt {2}}} \ right) \\ F _ {\ frac {3} {2}} (h) = {\ frac {1} {4 {\ sqrt {\ pi}}}} \ exp \ left (- {\ frac {h ^ { 2}} {4}} \ right) {\ sqrt {h}} \ left (\ left (h ^ {2} +1 \ right) K _ {\ frac {1} {4}} \ left ({\ frac {h ^ {2}} {4}} \ right) -h ^ {2} K _ {\ frac {3} {4}} \ left ({\ frac {h ^ {2}} {4}} \ right) \ right) \ end {align}}}{\displaystyle {\begin{aligned}F_{1}(h)={\frac {1}{\sqrt {2\pi }}}\exp \left(-{\frac {1}{2}}h^{2}\right)-{\frac {1}{2}}h\,\operatorname {erfc} \left({\frac {h}{\sqrt {2}}}\right)\\F_{\frac {3}{2}}(h)={\frac {1}{4{\sqrt {\pi }}}}\exp \left(-{\frac {h^{2}}{4}}\right){\sqrt {h}}\left(\left(h^{2}+1\right)K_{\frac {1}{4}}\left({\frac {h^{2}}{4}}\right)-h^{2}K_{\frac {3}{4}}\left({\frac {h^{2}}{4}}\right)\right)\end{aligned}}}

где erfc (z) означает дополнительную функцию ошибок, а K ν (z) {\ displaystyle K _ {\ nu} (z)}{\ displaystyle K _ {\ nu} (z)} - модифицированная функция Бесселя второго рода.

Для ситуаций, когда неровности на двух поверхностях имеют распределение по высоте, а пики можно считать сферическими, средним контактным давлением достаточно, чтобы вызвать текучесть, когда p av = 1,1 σ y ≈ 0, 39 σ 0 {\ displaystyle p _ {\ text {av}} = 1.1 \ sigma _ {y} \ приблизительно 0,39 \ sigma _ {0}}{\displaystyle p_{\text{av}}=1.1\sigma _{y}\approx 0.39\sigma _{0}}где σ y {\ displaystyle \ sigma _ {y}}\sigma_y- одноосный предел текучести и σ 0 {\ displaystyle \ sigma _ {0}}\sigma _{0}- твердость при вдавливании. Гринвуд и Уильямсон определили безразмерный параметр Ψ {\ displaystyle \ Psi}\Psi , названный индексом пластичности, который можно использовать для определения того, будет ли контакт эластичным или пластичным.

Модель Гринвуда-Уильямсона требует знания двух статистически зависимых величин; стандартное отклонение шероховатости поверхности и кривизны выступов неровностей. Альтернативное определение пластичности было дано Микичем. Податливость возникает, когда давление превышает одноосный предел текучести. Предел текучести пропорционален твердости при вдании σ 0 {\ displaystyle \ sigma _ {0}}\sigma _{0}, Микич определил индекс пластичности для упруго-пластического контакта как

Ψ = E ∗ h ′ σ 0>2 3. {\ displaystyle \ Psi = {\ frac {E ^ {*} h '} {\ sigma _ {0}}}>{\ frac {2} {3}} ~.}{\displaystyle \Psi ={\frac {E^{*}h'}{\sigma _{0}}}>{\ frac {2} {3}} ~.}

В этом определении Ψ {\ displaystyle \ Psi}\Psi представляет микрошероховатость в состоянии полной пластичности, требуется только одна статистическая, из, среднеквадратичная, статистическая, среднеквадратичная. Ψ < 2 3 {\displaystyle \Psi <{\frac {2}{3}}}{\displaystyle \Psi <{\frac {2}{3}}}система упруго во время контакта.

в соответствии с нагрузкой на нагрузку деформированной площади. от приложенной нормальной силы.

Обзор модели GT

Модель, предложенная Гринвудом и Триппом (GT), расширила модель GW на контакт между двумя шероховатыми поверхностями. Модель GT широко используется в области эластогидродинамического анализа.

Наиболее часто цитируемые уравнения, приводимые в модели GT, относятся к области контакта неровностей

A a = π 2 (η β σ) 2 AF 2 (λ), {\ displaystyle A_ {a} = \ pi ^ {2} (\ eta \ beta \ sigma) ^ {2} AF_ {2} (\ lambda),}{\displaystyle A_{a}=\pi ^{2}(\eta \beta \ sigma)^{2}AF_{2}(\lambda),}

и нагрузка на выступах

P = 8 2 15 π (η β σ) 2 σ β E ′ AF 5 2 (λ), {\ Displaystyle P = {\ frac {8 {\ sqrt {2}}} {15}} \ pi (\ eta \ beta \ sigma) ^ {2} {\ sqrt {\ frac {\ sigma} {\ beta}}} E'AF _ {\ frac {5} {2}} (\ lambda),}{\displaystyle P={\frac {8{\sqrt {2}}}{15}}\pi (\eta \beta \sigma)^{2}{\sqrt {\frac {\sigma }{\beta }}}E'AF_{\frac {5}{2}}(\lambda),}

где:

η β σ {\ displaystyle \ eta \ beta \ sigma }{\displaystyle \eta \beta \sigma }, параметр шероховатости,
A {\ displaystyle A}A , номинальная площадь контакта,
λ {\ displaystyle \ lambda}\lambda , масло Stribeck параметр фильма, впервые определенный Стрибеком \ cite {gt} как λ = h / σ {\ displaystyle \ lambda = h / \ sigma}{\displaystyle \lambda =h/\sigma },
E ′ {\ displaystyle E '}E', эффективный модуль упругости,
F 2, F 5 2 (λ) {\ displaystyle F_ {2}, F _ {\ frac {5} {2}} (\ lambda)}{\displaystyle F_{2},F_{\frac {5}{2}}(\lambda)}, введены статистические функции для соответствия предполагаемому гауссовскому распределению неровностей.

Точное решения для A a {\ displaystyle A_ {a}}A_{a}и P {\ displaystyle P}Pвпервые представлены Jedynak. Они выражаются как F n {\ displaystyle F_ {n}}F_{n}следующим образом:

F 2 = 1 2 (h 2 + 1) erfc ⁡ (h 2) - h 2 π exp ⁡ (- h 2 2) F 5 2 = 1 8 π exp ⁡ (- h 2 4) h 3 2 ((2 h 2 + 3) K 3 4 (h 2 4) - (2 h 2 + 5) K 1 4 (час 2 4)) {\ displaystyle {\ begin {align} F_ {2} = {\ frac {1} {2}} \ left (h ^ {2} +1 \ right) \ operatorname {erfc } \ left ({\ frac {h} {\ sqrt {2}}} \ right) - {\ frac {h} {\ sqrt {2 \ pi}}} \ exp \ left (- {\ frac {h ^ {2}} {2}} \ right) \\ F _ {\ frac {5} {2}} = {\ frac {1} {8 {\ sqrt {\ pi}}}} \ exp \ left (- {\ frac {h ^ {2}} {4}} \ right) h ^ {\ frac {3} {2}} \ left (\ left (2h ^ {2} +3 \ right) K _ {\ frac { 3} {4}} \ left ({\ frac {h ^ {2}} {4}} \ right) - \ left (2h ^ {2} +5 \ right) K _ {\ frac {1} {4} } \ left ({\ frac {h ^ {2}} {4}} \ right) \ right) \ end {align}}}{\displaystyle {\begin{aligned}F_{2}={\frac {1}{2}}\left(h^{2}+1\right)\operatorname {erfc} \left({\frac {h}{\sqrt {2}}}\right)-{\frac {h}{\sqrt {2\pi }}}\exp \left(-{\frac {h^{2}}{2}}\right)\\F_{\frac {5}{2}}={\frac {1}{8{\sqrt {\pi }}}}\exp \left(-{\frac {h^{2}}{4}}\right)h^{\frac {3}{2}}\left(\left(2h^{2}+3\right)K_{\frac {3}{4}}\left({\frac {h^{2}}{4}}\right)-\left(2h^{2}+5\right)K_{\frac {1}{4}}\left({\frac {h^{2}}{4}}\right)\rig ht)\end{aligned}}}

где erfc (z) означает дополнительную функцию ошибок, а K ν (z) {\ displaystyle K _ {\ nu} (z)}{\ displaystyle K _ {\ nu} (z)} - модифицированная функция Бесселя второго рода.

В статье можно найти исчерпывающий обзор существующих приближений к F 5 2 {\ displaystyle F _ {\ frac {5} {2}}}{\displaystyle F_{\frac {5}{2}}}. Новые предложения дают наиболее точные аппроксимации для F 5 2 {\ displaystyle F _ {\ frac {5} {2}}}{\displaystyle F_{\frac {5}{2}}}и F 2 {\ displaystyle F_ {2}}F_{2}, о которых сообщается в литературе. Они задаются следующими рациональными формулами, которые являются очень точными приближениями к интегралам F n (h) {\ displaystyle F_ {n} (h)}{\ displaystyle F_ {n} (h)} . Они рассчитаны для гауссова распределения неровностей

F n (h) = a 0 + a 1 h + a 2 h 2 + a 3 h 3 1 + b 1 h + b 2 h 2 + b 3 h 3 + б 4 час 4 + b 5 час 5 + b 6 час 6 ехр ⁡ (- час 2 2) {\ displaystyle F_ {n} (h) = {\ frac {a_ {0} + a_ {1} h + a_ { 2} h ^ {2} + a_ {3} h ^ {3}} {1 + b_ {1} h + b_ {2} h ^ {2} + b_ {3} h ^ {3} + b_ {4 } h ^ {4} + b_ {5} h ^ {5} + b_ {6} h ^ {6}}} \ exp \ left (- {\ frac {h ^ {2}} {2}} \ right)}{\displaystyle F_{n}(h)={\frac {a_{0}+a_{1}h+a_{2}h^{2}+a_{3}h^{3}}{1+b_{1}h+b_{2}h^{2}+b_{3}h^{3}+b_{4}h^{4}+b_{5}h^{5}+b_{6}h^{6}}}\exp \left(-{\frac {h^{2}}{2}}\right)}

Для F 2 (h) {\ displaystyle F_ {2} (h)}{\displaystyle F_{2}(h)}коэффициенты

[a 0, a 1, a 2, a 3] = [0,5, 0,182536384941, 0,039812283118, 0,003684879001] [b 1, b 2, b 3, b 4, b 5, b 6] = [1,960841785003, 1,708677456715, 0,856592986083, 0,264996791567, 0,04925784389313yle] \ begin выровнено} [a_ {0}, a_ {1}, a_ {2}, a_ {3}] = [0.5,0.182536384941,0.039812283118,0.003684879001] \\ [b_ {1}, b_ {2}, b_ {3 }, b_ {4}, b_ {5}, b_ {6}] = [1.960841785003,1.708677456715,0.856592986083,0.264996791567,0.049257843893,0.004640740133] \ end {align}}}{\displaystyle {\begin{aligned} [a_{0},a_{1},a_{2},a_{3}]=[0.5,0.182536384941,0.039812283118,0.003684879001]\\ [b_{1},b_{2},b_{3},b_{4},b_{5},b_{6}]=[1.960841785003,1.708677456715,0.856592986083,0.264996791567,0.049257843893,0.004640740133]\end{aligned}}}

Максимальная относительная ошибка: 1,68 × 10 - 7% {\ displ aystyle 1,68 \ times 10 ^ {- 7} \%}{\displaystyle 1.68\times 10^{-7}\%}.

Для F 5 2 (h) {\ displaystyle F _ {\ frac {5} {2}} (h)}{\displaystyle F_{\frac {5}{2}}(h)}коэффициенты:

[a 0, a 1, a 2, a 3] = [0,616634218997, 0,108855827811, 0,023453835635, 0,000449332509] [b 1, b 2, b 3, b 4, b 5, b 6] = [ 1.919948267476, 1.635304362591, 0.799392556572, 0.240278859212, 0.043178653945, 0.003863334276] {\ displaystyle {\ begin {выровнено} [a_ {0}, a_ {1}, a_ {2}, a_ {3}] = [0,616.10883421899738, 0.000449332509] \\ [b_ {1}, b_ {2}, b_ {3}, b_ {4}, b_ {5}, b_ {6}] = [1.919948267476,1.635304362591,0.799392556572,0.240278859212,0.043178653945276.00386 ] \ end {align}}}{\displaystyle {\begin{aligned} [a_{0},a_{1},a_{2},a_{3}]=[0.616634218997,0.108855827811,0.023453835635,0.000449332509]\\ [b_{1},b_{2},b_{3},b_{4},b_{5},b_{6}]=[1.919948267476,1.635304362591,0.799392556572,0.240278859212,0.043178653945,0.003863334276]\end{aligned}}}

Максимальная относительная погрешность 4,98 × 10-8% {\ displaystyle 4.98 \ times 10 ^ {- 8} \%}{\displaystyle 4.98\times 10^{-8}\%}.

Адгезионный контакт между упругими телами

Когда две твердые поверхности находятся в непосредственной близости, они испытывают притягивающие силы Ван-дер-Ваальса. Модель Ван-дер-Ваальса Брэдли позволяет рассчитать силу растяжения между двумя жесткими сферами с идеально гладкими поверхностями. Модель контакта Герца не считает адгезию возможной. Однако в конце 1960-х годов при сравнении теории Герца с экспериментами с контактом резиновых и стеклянных сфер было обнаружено несколько противоречий.

Было замечено, что, хотя теория Герца применялась при больших нагрузках, при низких нагрузках

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

Это указывало на то, что действуют силы сцепления. Модель Джонсона-Кендалла-Робертса (JKR) и модели Дерягина-Мюллера-Топорова (DMT) были первыми, в которых адгезия была включена в контакт Герца.

Модель жесткого контакта Брэдли

Обычно предполагается, что поверхностная сила между двумя атомными плоскостями на расстоянии z {\ displaystyle z}zдруг от друга может быть получено из потенциала Леннарда-Джонса. При таком предположении

F (z) = 16 γ 3 z ​​0 [(zz 0) - 9 - (zz 0) - 3] {\ displaystyle F (z) = {\ cfrac {16 \ gamma} {3z_ { 0}}} \ left [\ left ({\ cfrac {z} {z_ {0}}} \ right) ^ {- 9} - \ left ({\ cfrac {z} {z_ {0}}} \ right) ^ {- 3} \ right]}F (z) = {\ cfrac {16 \ gamma} {3z_ {0}}} \ left [\ left ( {\ cfrac {z} {z_ {0}}} \ right) ^ {{- 9}} - \ left ({\ cfrac {z} {z_ {0}}} \ right) ^ {{- 3}} \ справа]

где F {\ displaystyle F}F - сила (положительная при сжатии), 2 γ {\ displaystyle 2 \ gamma}2\gamma - общая поверхностная энергия обеих поверхностей на единицу площади, а z 0 {\ displaystyle z_ {0}}z_{0}- равновесное разделение двух атомных плоскостей.

В модели Брэдли был использован потенциал Леннарда-Джонса, чтобы найти силу сцепления между двумя жесткими сферами. Общая сила между сферами составляет

F a (z) = 16 γ π R 3 [1 4 (z z 0) - 8 - (z z 0) - 2]; 1 р = 1 р 1 + 1 р 2 {\ displaystyle F_ {a} (z) = {\ cfrac {16 \ gamma \ pi R} {3}} \ left [{\ cfrac {1} {4}} \ left ({\ cfrac {z} {z_ {0}}} \ right) ^ {- 8} - \ left ({\ cfrac {z} {z_ {0}}} \ right) ^ {- 2} \ right ] ~; ~~ {\ frac {1} {R}} = {\ frac {1} {R_ {1}}} + {\ frac {1} {R_ {2}}}}F_{a}(z)={\cfrac {16\gamma \pi R}{3}}\left[{\cfrac {1}{4}}\left({\cfrac {z}{z_{0}}}\right)^{{-8}}-\left({\cfrac {z}{z_{0}}}\right)^{{-2}}\right]~;~~{\frac {1}{R}}={\frac {1}{R_{1}}}+{\frac {1}{R_{2}}}

где R 1, R 2 {\ displaystyle R_ {1}, R_ {2}}R_{1},R_{2}- радиусы двух сфер.

Две сферы полностью разделяются, когда сила отрыва достигается в z = z 0 {\ displaystyle z = z_ {0}}z=z_{0}в этой точке

F а = F c = - 4 γ π R. {\ displaystyle F_ {a} = F_ {c} = - 4 \ gamma \ pi R.}F_{a}=F_{c}=-4\gamma \pi R.

Модель упругого контакта Джонсона-Кендалла-Робертса (JKR)

Схема площади контакта для модели JKR. Тест JKR с жестким валиком на деформируемом плоском материале: полный цикл

Чтобы учесть эффект адгезии при контакте Герца, Джонсон, Кендалл и Робертс сформулировали теорию адгезионного контакта JKR, используя баланс между накопленной упругой энергией и потерями в поверхностной энергии. Модель JKR учитывает влияние контактного давления и адгезии только внутри области контакта. Общее решение для распределения давления в зоне контакта в модели JKR:

p (r) = p 0 (1 - r 2 a 2) 1 2 + p 0 ′ (1 - r 2 a 2) - 1 2 {\ displaystyle p (r) = p_ {0} \ left (1 - {\ frac {r ^ {2}} {a ^ {2}}} \ right) ^ {\ frac {1} {2}} + p_ {0} '\ left (1 - {\ frac {r ^ {2}} {a ^ {2}}} \ right) ^ {- {\ frac {1} {2}}}}{\displaystyle p(r)=p_{0}\left(1-{\frac {r^{2}}{a^{2}}}\right)^{\frac {1}{2}}+p_{0}'\left(1-{\frac {r^{2}}{a^{2}}}\right)^{-{\frac {1}{2}}}}

Обратите внимание, что в первоначальной теории Герца термин, содержащий p 0 '{\ displaystyle p_ {0}'}p_{0}', не учитывался на том основании, что напряжение не могло поддерживаться в зоне контакта. Для контакта двух сфер

p 0 = 2 a E ∗ π R; п 0 ′ = - (4 γ E ∗ π a) 1 2 {\ displaystyle p_ {0} = {\ frac {2aE ^ {*}} {\ pi R}}; \ quad p_ {0} '= - \ left ({\ frac {4 \ gamma E ^ {*}} {\ pi a}} \ right) ^ {\ frac {1} {2}}}{\displaystyle p_{0}={\frac {2aE^{*}}{\pi R}};\quad p_{0}'=-\left({\frac {4\gamma E^{*}}{\pi a}}\right)^{\frac {1}{2}}}

где a {\ displaystyle a \, }a \, - радиус области контакта, F {\ displaystyle F}F - приложенная сила, 2 γ {\ displaystyle 2 \ gamma}2\gamma - полная поверхностная энергия обеих поверхностей на единицу площади контакта, R i, E i, ν i, i = 1, 2 {\ displaystyle R_ {i}, \, E_ {i}, \, \ nu _ {i}, ~~ i = 1,2}{\displaystyle R_{i},\,E_{i},\,\nu _{i},~~i=1,2}- радиусы, модули Юнга и коэффициенты Пуассона двух сфер, а

1 R = 1 р 1 + 1 р 2; 1 E ∗ знак равно 1 - ν 1 2 E 1 + 1 - ν 2 2 E 2 {\ displaystyle {\ frac {1} {R}} = {\ frac {1} {R_ {1}}} + {\ frac {1} {R_ {2}}}; \ quad {\ frac {1} {E ^ {*}}} = {\ frac {1- \ nu _ {1} ^ {2}} {E_ {1} }} + {\ frac {1- \ nu _ {2} ^ {2}} {E_ {2}}}}{\displaystyl e {\frac {1}{R}}={\frac {1}{R_{1}}}+{\frac {1}{R_{2}}};\quad {\frac {1}{E^{*}}}={\frac {1-\nu _{1}^{2}}{E_{1}}}+{\frac {1-\nu _{2}^{2}}{E_{2}}}}

Расстояние приближения между двумя сферами равно

d = π a 2 E * (P 0 + 2 p 0 ′) = a 2 R {\ displaystyle d = {\ frac {\ pi a} {2E ^ {*}}} \ left (p_ {0} + 2p_ {0} '\ right) = {\ frac {a ^ {2}} {R}}}{\displaystyle d={\frac {\pi a}{2E^{*}}}\left(p_{0}+2p_{0}'\right)={\frac {a^{2}}{R}}}

Уравнение Герца для площади контакта между двумя сферами, модифицированное с учетом поверхностной энергии, имеет вид

a 3 = 3 R 4 E ∗ (F + 6 γ π R + 12 γ π RF + (6 γ π R) 2) {\ displaystyle a ^ {3} = {\ frac {3R} {4E ^ {*}}} \ left (F + 6 \ gamma \ pi R + {\ sqrt {12 \ gamma \ pi RF + (6 \ gamma \ pi R) ^ {2}}} \ right)}{\displaystyle a^{3}={\frac {3R}{4E^{*}}}\left(F+6\gamma \pi R+{\sqrt {12\gamma \pi RF+(6\gamma \pi R)^{2}}}\right)}

Когда поверхностная энергия равна нулю, γ = 0 {\ displaystyle \ gamma = 0}\gamma = 0восстанавливается уравнение Герца для контакта между двумя сферами. Когда приложенная нагрузка равна нулю, радиус контакта равен

a 3 = 9 R 2 γ π E ∗ {\ displaystyle a ^ {3} = {\ frac {9R ^ {2} \ gamma \ pi} {E ^ {*}}}}{\displaystyle a^{3}={\frac {9R^{2}\gamma \pi }{E^{*}}}}

Растягивающая нагрузка, при которой сферы разделяются (т. Е. a = 0 {\ displaystyle a = 0}a=0), по прогнозам, будет

F c = - 3 γ π R {\ displaystyle F _ {\ text {c}} = - 3 \ gamma \ pi R \,}{\displaystyle F_{\text{c}}=-3\gamma \pi R\,}

Эта сила также называется силой отрыва . Обратите внимание, что эта сила не зависит от модулей двух сфер. Однако есть другое возможное решение для значения a {\ displaystyle a}aпри этой нагрузке. Это критическая область контакта ac {\ displaystyle a _ {\ text {c}}}{\displaystyle a_{\text{c}}}, задаваемая

ac 3 = 9 R 2 γ π 4 E ∗ {\ displaystyle a_ { \ text {c}} ^ {3} = {\ frac {9R ^ {2} \ gamma \ pi} {4E ^ {*}}}}{\displaystyl e a_{\text{c}}^{3}={\frac {9R^{2}\gamma \pi }{4E^{*}}}}

Если мы определим работу адгезии как

Δ γ = γ 1 + γ 2 - γ 12 {\ displaystyle \ Delta \ gamma = \ gamma _ {1} + \ gamma _ {2} - \ gamma _ {12}}\Delta \gamma =\gamma _{1}+\gamma _{2}-\gamma _{{12}}

где γ 1, γ 2 {\ displaystyle \ gamma _ {1}, \ gamma _ {2}}\gamma _{1},\gamma _{2}- энергии сцепления двух поверхностей, а γ 12 {\ displaystyle \ gamma _ { 12}}\gamma _{{12}}- член взаимодействия, мы можем записать радиус контакта JKR как

a 3 = 3 R 4 E ∗ (F + 3 Δ γ π R + 6 Δ γ π RF + (3 Δ γ π R) 2) {\ displaystyle a ^ {3} = {\ frac {3R} {4E ^ {*}}} \ left (F + 3 \ Delta \ gamma \ pi R + {\ sqrt {6 \ Delta) \ gamma \ pi RF + (3 \ Delta \ gamma \ pi R) ^ {2}}} \ right)}{\displaystyle a^{3}={\frac {3R}{4E^{*}}}\left(F+3\Delta \gamma \pi R+{\sqrt {6\Delta \gamma \pi RF+(3\Delta \gamma \pi R)^{2}}}\right)}

Растягивающая нагрузка при отрыве

F = - 3 2 Δ γ π R {\ displaystyle F = - {\ frac {3} {2}} \ Delta \ gamma \ pi R \,}{\ displaystyle F = - {\ frac {3} {2}} \ Delta \ gamma \ pi R \,}

, а критический радиус контакта определяется как

ac 3 = 9 R 2 Δ γ π 8 E ∗ {\ displaystyle a _ {\ текст {c}} ^ {3} = {\ frac {9R ^ {2} \ Delta \ gamma \ pi} {8E ^ {*}}}}{\displaystyle a_{\text{c}}^{3}={\frac {9R^{2}\Delta \gamma \pi }{8E^{*}}}}

Критическая глубина проникновения

dc = ac 2 R знак равно (R 1 2 9 Δ γ π 4 E ∗) 2 3 {\ displaystyle d _ {\ text {c}} = {\ frac {a_ {c} ^ {2}} {R}} = \ left ( R ^ {\ frac {1} {2}} {\ frac {9 \ Delta \ gamma \ pi} {4E ^ {*}}} \ right) ^ {\ frac {2} {3}}}{\displaystyle d_{\text{c}}={\frac {a_{c}^{2}}{R}}=\left(R^{\frac {1}{2}}{\frac {9\Delta \gamma \pi }{4E^{*}}}\right)^{\frac {2}{3}}}

Модель Дерягина-Мюллера-Топорова (DMT) упругого контакта

Модель Дерягина-Мюллера-Топорова (DMT) является альтернативной моделью для адгезионного контакта, которая предполагает, что профиль контакта остается таким же, как в контакте Герца, но с дополнительные привлекательные взаимодействия за пределами области контакта.

Радиус контакта между двумя сферами по теории ДМТ равен

a 3 = 3 R 4 E ∗ (F + 4 γ π R) {\ displaystyle a ^ {3} = {\ cfrac {3R } {4E ^ {*}}} \ left (F + 4 \ gamma \ pi R \ right)}a^{3}={\cfrac {3R}{4E^{*}}}\left(F+4\gamma \pi R\right)

и сила отрыва равна

F c = - 4 γ π R {\ displaystyle F_ { c} = - 4 \ gamma \ pi R \,}F_{c}=-4\gamma \pi R\,

Когда достигается сила отрыва, площадь контакта становится равной нулю, и нет сингулярности в контактных напряжениях на краю области контакта.

С точки зрения работы адгезии Δ γ {\ displaystyle \ Delta \ gamma}\Delta \gamma

a 3 = 3 R 4 E ∗ (F + 2 Δ γ π R) {\ displaystyle a ^ {3} = {\ cfrac {3R} {4E ^ {*}}} \ left (F + 2 \ Delta \ gamma \ pi R \ right)}a^{3}={\cfrac {3R}{4E^{*}}}\left(F+2\Delta \gamma \pi R\right)

и

F c = - 2 Δ γ π R {\ displaystyle F_ {c} = - 2 \ Delta \ gamma \ pi R \,}F_{c}=-2\Delta \gamma \pi R\,

параметр Табора

В 1977 году Табор показал, что очевидное противоречие между теориями JKR и ДМТ может быть решено, отметив, что две теории были крайними пределами единой теории, параметризованной параметром Табора (μ {\ displaystyle \ mu}\mu ), определенным как

μ : = dcz 0 ​​≈ [R (Δ γ) 2 E ∗ 2 z 0 3] 1 3 {\ displaystyle \ mu: = {\ frac {d_ {c}} {z_ {0}}} \ приблизительно \ left [{ \ frac {R (\ Delta \ gamma) ^ {2}} {{E ^ {*}} ^ {2} z_ {0} ^ {3}}} \ right] ^ {\ frac {1} {3} }}{\displaystyle \mu :={\frac {d_{c}}{z_{0}}}\approx \left[{\frac {R(\Delta \gamma)^{2}}{{E^{*}}^{2}z_{0}^{3}}}\right]^{\frac {1}{3}}}

где z 0 {\ displaystyle z_ {0}}z_{0}- равновесное расстояние между двумя контактирующими поверхностями. Теория JKR применима к большим, податливым сферам, для которых μ {\ displaystyle \ mu}\mu является большим. Теория ДМТ применима к маленьким жестким сферам с малыми значениями μ {\ displaystyle \ mu}\mu .

. Впоследствии Дерягин и его сотрудники, применив закон силы Брэдли к упругому полупространству, подтвердили, что параметр Табора увеличивается, сила отрыва падает от значения Брэдли 2 π R Δ γ {\ displaystyle 2 \ pi R \ Delta \ gamma}{\displaystyle 2\pi R\Delta \gamma }до значения JKR (3/2) π R Δ γ {\ displaystyle (3/2) \ pi R \ Delta \ gamma}{\ displaystyle (3/2) \ pi R \ Delta \ gamma} . Позже Гринвуд провел более подробные расчеты, показавшие S-образную кривую нагрузки / подхода, которая объясняет эффект подпрыгивания. Более эффективный метод выполнения расчетов и дополнительные результаты были даны Фэном

моделью упругого контакта Моугиса-Дагдейла

Схема зоны контакта для модели Моугиса-Дагдейла.

Дальнейшее улучшение Табора. Идея была предложена Моугисом, который представил поверхностную силу в терминах приближения Дагдейла, так что работа адгезии определяется как

Δ γ = σ 0 h 0 {\ displaystyle \ Delta \ gamma = \ sigma _ {0} ~ h_ {0}}\Delta \gamma =\sigma _{0}~h_{0}

где σ 0 {\ displaystyle \ sigma _ {0}}\sigma _{0}- максимальная сила, предсказанная потенциалом Леннарда-Джонса, а h 0 {\ displaystyle h_ {0}}h_{0}- максимальное разделение, полученное путем сопоставления областей под кривыми Дагдейла и Леннарда-Джонса (см. рисунок рядом). Это означает, что сила притяжения постоянна для z 0 ≤ z ≤ z 0 + h 0 {\ displaystyle z_ {0} \ leq z \ leq z_ {0} + h_ {0}}z_{0}\leq z\leq z_{0}+h_{0}. Дальнейшего проникновения при сжатии нет. Идеальный контакт происходит в области радиусом a {\ displaystyle a}aи силами сцепления величиной σ 0 {\ displaystyle \ sigma _ {0}}\sigma _{0}extend в область радиусом c>a {\ displaystyle c>a}c>a . В области a < r < c {\displaystyle aa<r<cдве поверхности разделены расстоянием h (r) {\ displaystyle h (r)}h(r)с h (a) = 0 {\ displaystyle h (a) = 0}h(a)=0и h (c) = h 0 {\ displaystyle h (c) = h_ { 0}}h(c)=h_{0}. Отношение m {\ displaystyle m}mопределяется как

m: = ca {\ displaystyle m: = {\ frac {c} { a}}}{\displaystyle m:={\frac {c}{a}}}.

В теории Моугиса-Дагдейла распределение силы тяги на поверхности делится на две части: одна из-за контактного давления Герца, а другая из-за адгезионного напряжения Дагдейла. Контакт Герца предполагается в области - a < r < a {\displaystyle -a. Вклад в поверхностную тягу от Her tz давление определяется как

p H (r) = (3 FH 2 π a 2) (1 - r 2 a 2) 1 2 {\ displaystyle p ^ {H} (r) = \ left ({\ frac {3F ^ {H}} {2 \ pi a ^ {2}}} \ right) \ left (1 - {\ frac {r ^ {2}} {a ^ {2}}} \ right) ^ {\ frac {1} {2}}}{\displaystyle p^{H}(r)=\left({\frac {3F^{H}}{2\pi a^{2}}}\right)\left(1-{\frac {r^{2}}{a^{2}}}\right)^{\frac {1}{2}}}

, где контактная сила Герца FH {\ displaystyle F ^ {H}}F^{H}определяется как

FH = 4 E ∗ a 3 3 R {\ displaystyle F ^ {H} = {\ frac {4E ^ {*} a ^ {3}} {3R}}}{\displaystyle F^{H}={\frac {4E^{*}a^{3}}{3R}}}

Проникновение из-за упругого сжатия составляет

d H = a 2 R { \ displaystyle d ^ {H} = {\ frac {a ^ {2}} {R}}}{\displaystyle d^{H}={\frac {a^{2}}{R}}}

Вертикальное смещение в r = c {\ displaystyle r = c}r = c равно

U ЧАС (с) знак равно 1 π R [a 2 (2 - м 2) грех - 1 ⁡ (1 м) + a 2 м 2 - 1] {\ displaystyle u ^ {H} (c) = { \ cfrac {1} {\ pi R}} \ left [a ^ {2} \ left (2-m ^ {2} \ right) \ sin ^ {- 1} \ left ({\ frac {1} {m }} \ right) + a ^ {2} {\ sqrt {m ^ {2} -1}} \ right]}{\displaystyle u^{H}(c)={\cfrac {1}{\pi R}}\left[a^{2}\ left(2-m^{2}\right)\sin ^{-1}\left({\frac {1}{m}}\right)+a^{2}{\sqrt {m^{2}-1}}\right]}

и разделение между двумя поверхностями в r = c {\ displaystyle r = c}r = c равно

час H (c) = c 2 2 R - d H + u H (c) {\ displaystyle h ^ {H} (c) = {\ frac {c ^ { 2}} {2R}} - d ^ {H} + u ^ {H} (c)}{\displaystyle h^{H}(c)={\frac {c^{2}}{2R}}-d^{H}+u^{H}(c)}

Распределение силы сцепления на поверхности n из-за клеевого напряжения Дагдейла

p D (r) = {- σ 0 π cos - 1 ⁡ [2 - m 2 - r 2 a 2 m 2 (1 - r 2 m 2 a 2)] для r ≤ a - σ 0 для a ≤ r ≤ c {\ displaystyle p ^ {D} (r) = {\ begin {cases} - {\ frac {\ sigma _ {0}} {\ pi}} \ cos ^ {-1} \ left [{\ frac {2-m ^ {2} - {\ frac {r ^ {2}} {a ^ {2}}}} {m ^ {2} \ left (1- { \ frac {r ^ {2}} {m ^ {2} a ^ {2}}} \ right)}} \ right] \ quad {\ text {for}} \ quad r \ leq a \\ - \ sigma _ {0} \ quad {\ text {for}} \ quad a \ leq r \ leq c \ end {cases}}}{\displaystyle p^{D}(r)={\begin{cases}-{\frac {\sigma _{0}}{\pi }}\cos ^{-1}\left[{\frac {2-m^{2}-{\frac {r^{2}}{a^{2}}}}{m^{2}\left(1-{\frac {r^{2}}{m^{2}a^{2}}}\right)}}\right]\quad {\text{for}}\quad r\leq a\\-\sigma _{0}\quad {\text{for}}\quad a\leq r\leq c\end{cases}}}

Тогда общая сила сцепления определяется как

FD = - 2 σ 0 м 2 a 2 [соз - 1 ⁡ (1 м) + 1 м 2 м 2 - 1] {\ displaystyle F ^ {D} = - 2 \ sigma _ {0} m ^ {2} a ^ {2} \ left [\ cos ^ {- 1} \ left ({\ frac {1} {m}} \ right) + {\ frac {1} {m ^ {2}}} {\ sqrt {m ^ {2} -1}} \ right]}{\ displaystyle F ^ {D} = - 2 \ sigma _ {0} m ^ {2} a ^ {2} \ left [\ cos ^ {- 1} \ left ({\ frac {1 } {m}} \ right) + {\ frac {1} {m ^ {2}}} {\ sqrt {m ^ {2} -1}} \ right]}

Сжатие из-за адгезии Дагдейла составляет

d D = - (2 σ 0 a E ∗) m 2 - 1 {\ displaystyle d ^ {D} = - \ left ( {\ frac {2 \ sigma _ {0} a} {E ^ {*}}} \ right) {\ sqrt {m ^ {2} -1}}}{\displaystyle d^{D}=-\left({\frac {2\sigma _{0}a}{E^{*}}}\right){\sqrt {m^{2}-1}}}

и разрыв в r = c {\ displaystyle r = c}r = c равно

h D (c) = (4 σ 0 a π E ∗) [m 2 - 1 cos - 1 ⁡ (1 m) + 1 - m ] {\ displaystyle h ^ {D} ( c) = \ left ({\ frac {4 \ sigma _ {0} a} {\ pi E ^ {*}}} \ right) \ left [{\ sqrt {m ^ {2} -1}} \ cos ^ {- 1} \ left ({\ frac {1} {m}} \ right) + 1-m \ right]}{\displaystyle h^{D}(c)=\left({\frac {4\sigma _{0}a}{\pi E^{*}}}\right)\left[{\sqrt {m ^{2}-1}}\cos ^{-1}\left({\frac {1}{m}}\right)+1-m\right]}

Тогда чистое сцепление с поверхностью контакта определяется как p (r) знак равно p H (r) + p D (r) {\ displaystyle p (r) = p ^ {H} (r) + p ^ {D} (r)}p(r)=p^{H}(r)+p^{D}(r)и чистая контактная сила равна F = FH + FD {\ displaystyle F = F ^ {H} + F ^ {D}}F=F^{H}+F^{D}. Когда h (c) = h H (c) + h D (c) = h 0 {\ displaystyle h (c) = h ^ {H} (c) + h ^ {D} (c) = h_ {0}}h(c)=h^{H}(c)+h^{D}(c)=h_{0}сцепление клея падает до нуля.

Безразмерные значения a, c, F, d {\ displaystyle a, c, F, d}a,c,F,dвводятся на этом этапе, которые определяются как

a ¯ = α a; c ¯: = α c; d ¯: = α 2 R d; α: = (4 E ∗ 3 π Δ γ R 2) 1 3; A ¯: = π c 2; F ¯ знак равно F π Δ γ R {\ displaystyle {\ bar {a}} = \ alpha a ~; ~~ {\ bar {c}}: = \ alpha c ~; ~~ {\ bar {d}}: = \ alpha ^ {2} Rd ~; ~~ \ alpha: = \ left ({\ frac {4E ^ {*}} {3 \ pi \ Delta \ gamma R ^ {2}}} \ right) ^ {\ frac {1} {3}} ~; ~~ {\ bar {A}}: = \ pi c ^ {2} ~; ~~ {\ bar {F}} = {\ frac {F} {\ pi \ Delta \ gamma R}}}{\displaystyle {\bar {a}}=\alpha a~;~~{\bar {c}}:=\alpha c~;~~{\bar {d}}:=\alpha ^{2}Rd~;~~\alpha :=\left({\frac {4E^{*}}{3\pi \Delta \gamma R^{2}}}\right)^{\frac {1}{3}}~;~~{\bar {A}}:=\pi c^{2}~;~~{\bar {F}}={\frac {F}{\pi \Delta \gamma R}}}

Кроме того, Моугис предложил параметр λ {\ displaystyle \ lambda}\lambda , который эквивалентен параметру Табора μ {\ displaystyle \ mu}\mu . Этот параметр определяется как

λ: = σ 0 (9 R 2 π Δ γ E ∗ 2) 1 3 ≈ 1,16 μ {\ displaystyle \ lambda: = \ sigma _ {0} \ left ({\ frac {9R } {2 \ pi \ Delta \ gamma {E ^ {*}} ^ {2}}} \ right) ^ {\ frac {1} {3}} \ приблизительно 1,16 \ mu}{\displaystyle \lambda :=\sigma _{0}\left({\frac {9R}{2\pi \Delta \gamma {E^{*}}^{2}}}\right)^{\frac {1}{3}}\approx 1.16\mu }

где ступенчатое когезионное напряжение σ 0 {\ displaystyle \ sigma _ {0}}\sigma _{0}равно теоретическому напряжению потенциала Леннарда-Джонса

σ th = 16 Δ γ 9 3 z 0 {\ displaystyle \ sigma _ {\ text {th}} = {\ frac {16 \ Delta \ gamma} {9 {\ sqrt {3}} z_ {0}}}}{\displaystyle \sigma _{\text{th}}={\frac {16\Delta \gamma }{9{\sqrt {3}}z_{0}}}}

Чжэн и Ю предложили другое значение для ступенчатого когезионного напряжения

σ 0 знак равно ехр ⁡ (- 223 420) ⋅ Δ γ z 0 ≈ 0,588 Δ γ z 0 {\ displaystyle \ sigma _ {0} = \ exp \ left (- {\ frac {223} {420}} \ справа) \ cdot {\ frac {\ Delta \ gamma} {z_ {0}}} \ приблизительно 0,588 {\ frac {\ Delta \ gamma} {z_ {0}}}}{\displaystyle \sigma _{0}=\exp \left(-{\frac {223}{420}}\right)\cdot {\frac {\Delta \gamma }{z_{0}}}\approx 0.588{\frac {\Delta \gamma }{z_{0}}}}

для соответствия потенциалу Леннарда-Джонса, что приводит к

λ ≈ 0,663 μ {\ displaystyle \ lambda \ приблизительно 0,663 \ mu}\lambda \approx 0.663\mu

Тогда чистая контактная сила может быть выражена как

F ¯ = a ¯ 3 - λ a ¯ 2 [m 2-1 + м 2 сек - 1 ⁡ м] {\ displaystyle {\ bar {F}} = {\ bar {a}} ^ {3 } - \ lambda {\ bar {a}} ^ {2} \ left [{\ sqrt {m ^ {2} -1}} + m ^ {2} \ sec ^ {- 1} m \ right]}{\displaystyle {\bar {F}}={\bar {a}}^{3}-\lambda {\bar {a}}^{2}\left[{\sqrt {m^{2}-1}}+m^{2}\sec ^{-1}m\right]}

и упругое сжатие как

d ¯ = a ¯ 2 - 4 3 λ a ¯ m 2 - 1 {\ displaystyle {\ bar {d}} = {\ bar {a}} ^ {2} - { \ frac {4} {3}} ~ \ lambda {\ bar {a}} {\ sqrt {m ^ {2} -1}}}{\displaystyle {\bar {d}}={\bar {a}}^{2}-{\frac {4}{3}}~\lambda {\bar {a}}{\sqrt {m^{2}-1}}}

Уравнение когезионного зазора между двумя телами принимает вид

λ a ¯ 2 2 [(m 2 - 2) sec - 1 ⁡ m + m 2 - 1] + 4 λ a ¯ 3 [m 2 - 1 sec - 1 ⁡ m - m + 1] = 1 {\ displaystyle {\ frac {\ lambda {\ bar {a}} ^ {2}} {2}} \ left [\ left (m ^ {2} -2 \ right) \ sec ^ {- 1} m + {\ sqrt {m ^ {2} -1}} \ right] + {\ frac {4 \ lambda {\ bar {a}}} {3}} \ left [{\ sqrt {m ^ {2} -1}} \ sec ^ {- 1} m-m + 1 \ right] = 1}{\displaystyle {\frac {\lambda {\bar {a}}^{2}}{2}}\left[\left(m^{2}-2\right)\sec ^{-1}m+{\sqrt {m^{2}-1}}\right]+{\frac {4\lambda {\bar {a}}}{3}}\left[{\sqrt {m^{2}-1}}\sec ^{-1}m-m+1\right]=1}

Это уравнение может быть решено для получения значений c {\ displaystyle c}cдля различных значений a {\ displaystyle a}aи λ {\ displaystyle \ lambda}\lambda . Для больших значений λ {\ displaystyle \ lambda}\lambda , m → 1 {\ displaystyle m \ rightarrow 1}m\rightarrow 1и получается модель JKR. Для малых значений λ {\ displaystyle \ lambda}\lambda извлекается модель DMT.

Модель Карпика-Оглетри-Салмерона (COS)

Модель Моугиса-Дагдейла может быть решена итеративно, только если значение λ {\ displaystyle \ lambda}\lambda Не известно априори. Приближенное решение Карпика-Оглетри-Салмерона упрощает процесс за счет использования следующего соотношения для определения радиуса контакта a {\ displaystyle a}a:

a = a 0 (β) (β + 1 - F / F c ( β) 1 + β) 2 3 {\ displaystyle a = a_ {0} (\ beta) \ left ({\ frac {\ beta + {\ sqrt {1-F / F_ ​​{c} (\ beta)}}}} {1+ \ beta}} \ right) ^ {\ frac {2} {3}}}{\displaystyle a=a_{0}(\beta)\left({\frac {\beta +{\sqrt {1-F/F_{c}(\beta)}}}{1+\beta }}\right)^{\frac {2}{3}}}

где a 0 {\ displaystyle a_ {0}}a_{0}- это контактная область в нулевая нагрузка, а β {\ displaystyle \ beta}\beta - параметр перехода, который связан с λ {\ displaystyle \ lambda}\lambda с помощью

λ ≈ - 0,924 ln ⁡ (1 - 1,02 β) {\ displaystyle \ lambda \ приблизительно -0,924 \ ln (1-1,02 \ beta)}{\ displaystyle \ lambda \ приблизительно -0,924 \ ln (1–1,02 \ бета)}

Случай β = 1 {\ displaystyle \ beta = 1}\ beta = 1 точно соответствует теории JKR, а β = 0 {\ displaystyle \ beta = 0}\beta =0соответствует теории DMT. Для промежуточных случаев 0 < β < 1 {\displaystyle 0<\beta <1}0<\beta <1модель COS близко соответствует решению Моугиса-Дагдейла для 0,1 < λ < 5 {\displaystyle 0.1<\lambda <5}0.1<\lambda <5.

Влияние формы контакта

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

См. Также
Ссылки
Внешние ссылки
  • [2] : Подробнее о контактных напряжениях и эволюции подшипниковых напряжений уравнений можно найти в этой публикации исследовательского центра NASA Glenn возглавит Bearing НАСА, трансмиссии и раздел передач, Эрвин Зарецкий
  • [3] :. MATLAB процедура для решения линейной упругой контактной задачи механики под названием ; «Решение LCP задачи линейной упругой контактной механики» предоставляется при обмене файлами в MATLAB Central.
  • [4] : Калькулятор контактной механики.
  • [5] : подробные расчеты и формулы теории JKR для двух сфер.
  • [5] : код Matlab для анализа контактов Герца (включает линейные, точечные и эллиптические случаи).
  • [6] : JKR, MD и DMT модели адгезии (процедуры Matlab).
Последняя правка сделана 2021-05-15 10:46:35
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте