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

Представим прямоугольный кирпич, лежащий на горизонтальном столе. Если мы потянем его с помощью динамометра, когда он лежит плашмя на своей широкой грани (номинальная площадь \(A_1 = 200\,\text{см}^2\)), а затем поставим тот же кирпич на узкое ребро (\(A_2 = 50\,\text{см}^2\)), обыденная интуиция подсказывает, что широкая грань должна создавать гораздо большее сопротивление. Однако эксперимент показывает: сила трения \(F_f\) остается абсолютно одинаковой при четырехкратном изменении площади опоры! И наоборот: если трение действительно не зависит от площади контакта, почему болиды Формулы-1 используют сверхширокие безпротекторные шины (слики)?

Второй классический миф гласит, что конькобежец скользит по льду, поскольку высокое давление под узким лезвием конька понижает точку плавления льда и образует водную пленку. Однако строгий расчет по термодинамическому уравнению Клаузиуса — Клапейрона показывает, что конькобежец массой \(75\,\text{кг}\) снижает температуру плавления льда всего на \(0{,}12^\circ\text{C}\), что никак не объясняет скольжение при \(-10^\circ\text{C}\) или \(-20^\circ\text{C}\)!

В этой статье мы обратимся к трибологии (науке о трении, износе и смазке) и физике поверхности на наноуровне, чтобы разрешить эти парадоксы: от различия между кажущейся (номинальной) площадью контакта (\(A_{\text{apparent}}\)) и фактической площадью контакта (\(A_{\text{real}}\)) до динамики прерывистого скольжения (stick-slip), квазижидкого слоя (QLL) на поверхности льда и структурной сверхсмазываемости (superlubricity).

1. Классическая загадка: законы Амонтона — Кулона и крах гипотезы геометрического зацепления

1.1 От Леонардо да Винчи (1493) к Амонтону (1699) и Кулону (1785)

Еще около 1493 года Леонардо да Винчи в своих записных книжках (Арунделовский кодекс) экспериментально установил, что сила трения пропорциональна нагрузке и не зависит от площади соприкосновения. Переоткрытые Гийомом Амонтоном (1699), дополненные Леонардом Эйлером (1748), который ввел коэффициент \(\mu = \tan\theta_c\), и обобщенные Шарлем Огюстеном де Кулоном (1785), эти правила образуют три классических закона сухого трения:

Гипотеза жестких микрозубьев с углом наклона \(\theta\) (\(F_f = F_N \tan\theta\)) противоречит первому началу термодинамики: на абсолютно жестком или чисто упругом склоне работа \(W_{\text{up}} = F_N \Delta h\), затраченная на подъем на вершину микровыступа, полностью возвращается (\(W_{\text{down}} = -F_N \Delta h\)) при спуске с обратной стороны зубца (\(\oint \mathbf{F}\cdot d\mathbf{r} = 0\)). Истинное трение требует необратимой деформации и диссипации механической энергии в колебания кристаллической решетки — фононы.

2. Номинальная и фактическая площадь контакта: как микровыступы (асферити) решают парадокс кирпича

2.1 «Положить перевернутую Швейцарию на Австрию»: адгезионная теория Боудена и Тейбора (1950)

Как показали кембриджские физики Фрэнк Филип Боуден и Дэвид Тейбор (1950), даже зеркально отполированный стальной брусок под атомно-силовым микроскопом (АСМ) имеет среднеквадратичную шероховатость \(R_q \sim 10\,\text{нм}\text{–}1\,\text{мкм}\), что в тысячи раз превышает радиус действия межатомных сил (\(0{,}2\text{–}0{,}4\,\text{нм}\)). По знаменитому выражению Боудена: «Сложить два твердых тела вместе — все равно что перевернуть Швейцарию вверх дном и поставить ее на Австрию: соприкоснутся лишь самые высокие альпийские вершины».

Эти микроскопические пики называются микровыступами (asperities). Сумма площадей реального атомного контакта на вершинах микровыступов образует фактическую площадь контакта (\(A_{\text{real}}\)), которая ничтожно мала по сравнению с кажущейся (номинальной) площадью (\(A_{\text{apparent}}\)):

Поскольку вся нормальная нагрузка \(F_N\) сосредоточена на \(A_{\text{real}}\), локальное напряжение \(p_{\text{local}} = F_N / A_{\text{real}}\) мгновенно достигает пластической твердости материала \(H \approx 3\sigma_y\). Вершины микровыступов пластически сминаются до тех пор, пока фактическая площадь не вырастет ровно настолько, чтобы удержать нагрузку \(F_N\):

В точках атомного сближения возникают адгезионные микроспайки (холодная сварка) с пределом прочности на сдвиг \(\tau_s\). Для скольжения необходимо срезать эти спайки:

  flowchart TD
    A["Макроскопическая нагрузка F_N и номинальная площадь A_apparent"] --> B["Кирпич плашмя (A_1): Малое номинальное давление P_1 = F_N / A_1"]
    A --> C["Кирпич на ребре (A_2 = A_1 / 4): Номинальное давление в 4 раза выше P_2 = 4 P_1"]
    B --> D["Множество микровыступов деформируются слабо"]
    C --> E["Меньшее число выступов сминается в 4 раза сильнее + вступают новые пики"]
    D --> F["Фактическая площадь контакта строго одинакова: A_real = F_N / H (или интеграл Гринвуда-Уильямсона)"]
    E --> F
    F --> G["Срез атомных адгезионных микроспаек: F_f = τ_s × A_real"]
    G --> H["Закон Амонтона-Кулона: F_f = μ F_N (не зависит от A_apparent)"]

Парадокс кирпича полностью разрешен: сила трения прямо пропорциональна ФАКТИЧЕСКОЙ площади контакта \(A_{\text{real}}\)! При повороте кирпича на узкую грань номинальное давление возрастает вчетверо, сильнее сжимая микровыступы, так что сумма \(A_{\text{real}} = F_N / H\) остается неизменной.

2.2 Упругая статистическая механика Гринвуда — Уильямсона (1966) и пределы применимости закона

Для одиночной упругой сферы радиуса \(R\) теория Герца (1882) дает \(A_{\text{single}} \propto w^{2/3}\). Однако Дж. Гринвуд и Дж. Уильямсон (1966) доказали, что при интегрировании по статистическому распределению высот вершин \(\phi(z) = \frac{1}{\sigma_s} e^{-z/\sigma_s}\) вовлечение новых микровыступов оставляет среднюю площадь одного пятна постоянной (\(\pi R \sigma_s\)), тогда как число пятен растет линейно с \(F_N\):

Когда закон Амонтона нарушается? Когда \(A_{\text{real}} \to A_{\text{apparent}}\):

1. Мягкие эластомеры и гоночные слики Формулы-1 (теория JKR, 1971): Из-за малого модуля Юнга (\(E \sim 1\text{–}10\,\text{МПа}\)) и работы адгезии \(\Delta\gamma\) по уравнению Джонсона — Кендалла — Робертса (JKR) \(a_{\text{JKR}}^3 = \frac{3R}{4E^*}(w + 3\pi R \Delta\gamma + \sqrt{6\pi R \Delta\gamma w + (3\pi R \Delta\gamma)^2})\) эффективный коэффициент трения падает с ростом давления (\(\mu_{\text{eff}} \propto P_{\text{nom}}^{-1/3}\)). Широкие слики увеличивают \(A_{\text{apparent}}\), снижая давление и максимизируя сцепление.
2. Лапки геккона (setae и spatulae): Миллионы кератиновых волосков, разветвленных на нанолопаточки шириной \(100\text{–}200\,\text{нм}\), снижают эффективный модуль ниже критерия Дальквиста (\(E_{\text{eff}} < 100\,\text{кПа}\)), обеспечивая \(A_{\text{real}} \approx A_{\text{apparent}}\) за счет сил Ван-дер-Ваальса.
3. Холодная сварка в космическом сверхвысоком вакууме: Без защитной оксидной пленки (\(1\text{–}5\,\text{нм}\)) чистые металлы в вакууме (\(P < 10^{-10}\,\text{Па}\)) мгновенно свариваются (\(\mu \to \infty\), причина заклинивания антенны зонда «Галилео» в 1991 году).

3. Динамика прерывистого скольжения (Stick-Slip) и закон Дитериха — Руины

Блок массы \(m\), тянущийся со скоростью \(v_0\) через пружину жесткости \(k\), подчиняется уравнению \(m \ddot{x} = k(v_0 t - x) - F_f(\dot{x}, \theta)\). После фазы покоя (Stick) длительностью \(t_{\text{stick}} = \frac{(\mu_s - \mu_k)F_N}{k v_0}\) происходит резкий срыв (Slip) с собственной частотой \(\omega_0 = \sqrt{k/m}\):

Почему \(\mu_s > \mu_k\)? Из-за ползучести (крипа) микровыступов фактическая площадь растет логарифмически со временем покоя: \(A_{\text{real}}(t) = A_0 [1 + c \ln(1 + t/\tau_0)]\) (старение контакта / Contact Aging). Скоростной и фазовый закон трения Дитериха — Руины (Rate-and-State Friction Law, 1979–1983) связывает скорость (\(v\)) и переменную состояния (\(\theta\), критическое смещение \(D_c\)):

В стационарном режиме (\(\theta_{\text{ss}} = D_c / v\)) имеем \(\mu_{\text{ss}}(v) = \mu_0 + (a - b)\ln(v/v_*)\). В режиме скоростного разупрочнения (\(a - b < 0\)) при жесткости системы ниже критической жесткости \(k_c\) возникает автоколебательная неустойчивость Хопфа (стик-слип):

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

4. Почему лед на самом деле скользкий? Опровержение плавления под давлением и квазижидкий слой (QLL)

Рассчитаем понижение точки плавления льда (\(\text{I}_h\)) по уравнению Клаузиуса — Клапейрона (\(T_m = 273{,}15\,\text{К}\), скрытая теплота плавления \(L_f = 3{,}335 \times 10^5\,\text{Дж/кг}\), \(\rho_{\text{liq}} = 999{,}8\,\text{кг/м}^3\), \(\rho_{\text{ice}} = 916{,}7\,\text{кг/м}^3\), \(\Delta v = -9{,}07 \times 10^{-5}\,\text{м}^3/\text{кг}\)):

Для конькобежца массой \(75\,\text{кг}\) (\(F_N = 735\,\text{Н}\)) на лезвии \(30\,\text{см} \times 1{,}5\,\text{мм}\) (\(A = 4{,}5 \times 10^{-4}\,\text{м}^2\)) давление \(\Delta P = 1{,}63\,\text{МПа}\) снижает температуру плавления лишь на:

Плавление под давлением не может объяснить катание при \(-10^\circ\text{C}\)! Три истинных механизма таковы:

1. Поверхностное предплавление и квазижидкий слой Фарадея (QLL, 1850): Поверхностные молекулы \(\text{H}_2\text{O}\) имеют лишь 2 или 3 водородные связи вместо 4 (недокоординированные молекулы), свободно колеблются и образуют квазижидкий слой (QLL) толщиной \(1\text{–}20\,\text{нм}\) вплоть до \(-30^\circ\text{C}\) при нулевом давлении.
2. Фрикционный разогрев (Боуден и Хьюз, 1939): Тепловой поток \(q = \mu P v\) в тонком скин-слое \(\delta_{\text{th}} \sim \sqrt{\alpha_{\text{th}} L / v}\) мгновенно нагревает контакт до \(0^\circ\text{C}\) (поэтому полиэтиленовые лыжи скользят намного лучше металлических).
3. Вязкоупругая реология межфазной воды (Canale, Bocquet & Bonn, 2019, Phys. Rev. X): Нанопленка воды обладает вязкостью в \(50\text{–}100\) раз выше обычной воды и модулем упругости \(G'(\omega)\), создающим упругогидродинамическую подъемную силу.

5. От кривой Штрибека к структурной сверхсмазываемости

Помимо гидродинамической смазки, описываемой уравнением Рейнольдса (1886) на кривой Штрибека:

модель Прандтля — Томлинсона (1928–1929) описывает атомное сухое трение через потенциал \(V_{\text{PT}}(x, X) = -\frac{U_0}{2}\cos(2\pi x/a) + \frac{1}{2}k_{\text{eff}}(x - X)^2\) и параметр \(\eta_{\text{PT}} = \frac{2\pi^2 U_0}{k_{\text{eff}} a^2}\). При \(\eta_{\text{PT}} \le 1\) атомные скачки стик-слип и излучение фононов исчезают.

  flowchart TD
    A["Сухой контакт двух атомарно гладких кристаллических решеток (модель Френкеля-Конторовой)"] --> B["Соизмеримый контакт: Рациональное отношение b/a или угол θ = 0°, 60°"]
    A --> C["Несоизмеримый контакт: Иррациональное b/a или угол θ ≠ 0° (муаровая сверхрешетка)"]
    B --> D["Все N атомов синхронно попадают в минимумы потенциала"]
    D --> E["Коллективный барьер: F_static ∝ N ∝ A_real (Высокое трение и стик-слип)"]
    C --> F["Атомы, идущие в гору (+∇U), точно уравновешены атомами, катящимися с горы (-∇U)"]
    F --> G["Деструктивная интерференция боковых сил выше перехода Обри: Σ F_i ≈ 0"]
    G --> H["Структурная сверхсмазываемость: Исчезновение трения (μ < 0.0001) без жидкой смазки!"]

В многоатомной модели Френкеля — Конторовой Мотохиса Хирано и Кадзумаса Синдзё (1990) предсказали структурную сверхсмазываемость (подтвержденную на повернутом графите Dienwiebel et al., 2004): при несоизмеримом контакте (муаровые сверхрешетки графена или \(\text{MoS}_2\) выше перехода Обри) боковые силы на атомах взаимно уничтожаются за счет деструктивной интерференции (\(\sum F_i \approx 0\)), снижая коэффициент сухого трения до \(\mu < 0{,}0001\)!

6. Сводная таблица трибологических режимов