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

Представим прямоугольный кирпич, лежащий на горизонтальном столе. Если мы потянем его с помощью динамометра, когда он лежит плашмя на своей широкой грани (номинальная площадь $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. Сводная таблица трибологических режимов