Уравнение Хилла для мышцы: сила и скорость

Уравнение Хилла связывает силу, которую развивает мышца, со скоростью её укорочения: чем тяжелее нагрузка, тем медленнее мышца сокращается. Это эмпирическое соотношение Арчибалд Хилл вывел в 1938 году из опытов с теплопродукцией портняжной мышцы лягушки, и оно стало основой количественной физиологии мышечного сокращения. В этой статье разберём само характеристическое уравнение, выведем из него формулу скорости укорочения, поясним смысл констант $a$ и $b$ и посчитаем, при какой нагрузке мышца отдаёт максимум мощности. Калькулятор ниже строит кривую сила-скорость по вашим параметрам.

Характеристическое уравнение Хилла

В классической форме уравнение Хилла записывают как соотношение между силой $F$ (нагрузкой) и скоростью укорочения $v$:

$(F + a)(v + b) = (F_0 + a)\, b$

Здесь $F_0$ - максимальная изометрическая сила, то есть нагрузка, при которой мышца уже не укорачивается ($v = 0$); $a$ и $b$ - две постоянные Хилла. Правая часть уравнения - константа, зависящая только от свойств мышцы. Произведение $(F + a)(v + b)$ остаётся неизменным при любой нагрузке, поэтому график зависимости $v$ от $F$ - это равнобочная гипербола, сдвинутая в новые оси на $-a$ по силе и $-b$ по скорости.

Физический смысл прост: при увеличении нагрузки $F$ скорость $v$ должна падать так, чтобы произведение сохранялось. Когда нагрузка достигает $F_0$, скорость обращается в ноль - это изометрический режим. Если же нагрузки нет вовсе ($F = 0$), скорость максимальна.

Вывод формулы скорости укорочения

Из характеристического уравнения легко выразить скорость как функцию нагрузки. Раскроем и решим относительно $v$:

$v = \frac{b\,(F_0 - F)}{F + a}$

Это рабочая форма уравнения Хилла. Подставив $F = 0$, получаем максимальную скорость укорочения:

$V_{\max} = \frac{b\,F_0}{a}$

а при $F = F_0$ числитель обнуляется и $v = 0$. Две точки - $(0, V_{\max})$ и $(F_0, 0)$ - это концы кривой, а между ними тянется гипербола.

Рассчитать для своей задачи

Кривизна гиперболы задаётся отношением $a/F_0$: чем оно меньше, тем сильнее кривая прогнута и тем резче падает скорость при росте нагрузки. Для большинства скелетных мышц $a/F_0 \approx 0{,}25$, причём из размерностей следует, что одновременно $b/V_{\max} = a/F_0$, поэтому одна безразмерная константа задаёт форму всей кривой.

Физиологический смысл констант a и b

Хилл показал, что его константы - не просто подгоночные коэффициенты. Постоянная $a$ имеет размерность силы и связана с теплотой укорочения: при сокращении мышца выделяет дополнительное тепло, пропорциональное пройденному пути, и коэффициент пропорциональности численно равен $a$. Постоянная $b$ имеет размерность скорости и определяет, как быстро высвобождается эта энергия.

Через эти же константы выражается максимальная скорость $V_{\max}$, которая отражает скорость скольжения актиновых и миозиновых нитей друг относительно друга. Молекулярную основу этого скольжения описывает теория скользящих нитей: мостики миозина циклически прикрепляются к актину, и чем выше нагрузка, тем дольше держится каждый мостик, тем медленнее идёт скольжение - именно это и даёт гиперболическую форму кривой Хилла.

Мощность мышцы и оптимальная нагрузка

Механическая мощность мышцы - это произведение силы на скорость: $P = F\,v$. Подставив уравнение Хилла, получим

$P = \frac{b\,F\,(F_0 - F)}{F + a}$

При $F = 0$ скорость максимальна, но сила нулевая, поэтому мощность равна нулю. При $F = F_0$ скорость нулевая - снова $P = 0$. Между этими крайностями мощность достигает максимума.

Рассчитать для своей задачи

Дифференцируя $P(F)$ и приравнивая производную к нулю, можно показать, что максимум мощности приходится примерно на $F \approx 0{,}3\,F_0$ (для типичного $a/F_0 = 0{,}25$ - точнее, около $0{,}31\,F_0$). Эта закономерность объясняет, почему в спорте максимальную мощность развивают при средних, а не предельных нагрузках - связанные эффекты разбираются в материале про коэффициент полезного действия мышцы. Калькулятор выше отмечает точку пика на кривой и показывает значение оптимальной нагрузки.

Связь с режимами сокращения

Уравнение Хилла описывает концентрический режим - мышца укорачивается под нагрузкой меньше $F_0$. При $F = F_0$ мышца работает изометрически: длина постоянна, скорость нулевая. Эти крайние точки кривой соответствуют двум каноническим режимам, подробно разобранным в статье про изотонический и изометрический режимы сокращения.

Если нагрузка превышает $F_0$, мышца удлиняется под нагрузкой (эксцентрический режим), и классическое уравнение Хилла туда не распространяется - там действуют другие соотношения. Поэтому при подстановке всегда следите, чтобы $F$ лежала в диапазоне от $0$ до $F_0$.

Пример расчёта по уравнению Хилла

Пусть $F_0 = 100$ Н, $a = 25$ Н (то есть $a/F_0 = 0{,}25$), $b = 10$ мм/с. Сначала найдём максимальную скорость:

$V_{\max} = \frac{b\,F_0}{a} = \frac{10 \cdot 100}{25} = 40 \text{ мм/с}$

Теперь скорость при нагрузке $F = 40$ Н:

$v = \frac{b\,(F_0 - F)}{F + a} = \frac{10\,(100 - 40)}{40 + 25} = \frac{600}{65} \approx 9{,}2 \text{ мм/с}$

Мощность в этой точке: $P = F\,v = 40 \cdot 9{,}2 \approx 369$ единиц (Н·мм/с). Перебрав нагрузки, убедимся, что пик мощности приходится на $F \approx 31$ Н - близко к теоретическому $0{,}3\,F_0$.

Как определяют константы из эксперимента

На практике постоянные Хилла находят из серии опытов с постнагрузочными сокращениями: мышце задают разные грузы $F$ и измеряют начальную скорость укорочения $v$. Получив набор точек, их линеаризуют. Удобный приём - преобразовать уравнение Хилла к линейному виду. Если обозначить $P_0 = F_0$ и переписать соотношение через $(F_0 - F)/v$, то зависимость становится прямой линией, по наклону и пересечению которой считывают $a$ и $b$.

Изометрическую силу $F_0$ измеряют отдельно - это максимальное усилие при фиксированной длине, когда укорочения нет. Максимальную скорость $V_{\max}$ либо экстраполируют к нулевой нагрузке, либо измеряют методом быстрого освобождения, когда мышце мгновенно снимают нагрузку. Полученные $F_0$, $V_{\max}$ и $a/F_0$ полностью задают кривую, и дальше уравнение Хилла предсказывает скорость при любой промежуточной нагрузке без новых измерений. Именно эта предсказательная сила сделала модель стандартом в учебниках физиологии и биомеханики, а также в моделях мышечной работы для спорта и протезирования.

Частые ошибки

• Путают $F_0$ с текущей нагрузкой. $F_0$ - это предельная изометрическая сила (постоянная мышцы), а $F$ - переменная нагрузка в конкретном опыте. В формулу скорости подставляют именно текущую $F$.
• Берут $F > F_0$. При нагрузке выше максимальной силы числитель $(F_0 - F)$ становится отрицательным, и формула даёт бессмысленную отрицательную скорость. Уравнение Хилла работает только для $0 \le F \le F_0$.
• Считают $a$ и $b$ независимыми. Из размерностей следует жёсткая связь $b = a\,V_{\max}/F_0$, поэтому форму кривой задаёт одно отношение $a/F_0$, а не два свободных параметра.
• Ищут максимум мощности на краях. На концах кривой ($F = 0$ и $F = F_0$) мощность нулевая. Максимум - внутри, около трети от $F_0$.
• Смешивают размерности. Скорость в мм/с, сила в ньютонах - следите, чтобы $a$ было в единицах силы, а $b$ в единицах скорости.

FAQ

Чем характеристическое уравнение Хилла отличается от формулы скорости? Это одно и то же соотношение в разных формах. Характеристическое уравнение $(F + a)(v + b) = (F_0 + a)b$ записано симметрично, а формула $v = b(F_0 - F)/(F + a)$ - то же уравнение, решённое относительно скорости для подстановки чисел.

Почему кривая сила-скорость именно гипербола, а не прямая? Потому что произведение $(F + a)(v + b)$ постоянно. Геометрически условие «произведение двух смещённых координат равно константе» задаёт равнобочную гиперболу. Молекулярно это следствие того, что под нагрузкой миозиновые мостики отцепляются медленнее.

Применимо ли уравнение Хилла к одиночному волокну и к целой мышце? Да, форма гиперболы сохраняется на обоих уровнях, меняются лишь численные значения $F_0$, $a$ и $b$. Для целой мышцы это интегральные характеристики, для волокна - локальные.

Коротко

Уравнение Хилла $(F + a)(v + b) = (F_0 + a)b$ описывает мышцу как гиперболическую зависимость скорости укорочения от нагрузки. Из него выводится рабочая формула $v = b(F_0 - F)/(F + a)$, максимальная скорость $V_{\max} = bF_0/a$ и кривая мощности с пиком около $0{,}3\,F_0$. Константы $a$ и $b$ связаны через $a/F_0 = b/V_{\max}$ и несут энергетический смысл, а сама модель остаётся базовой для биомеханики и спортивной физиологии.