CC BY
27
УДК 531/534: [57+61]
Российский Журнал
www.biomech.ru
АКТУАТОРНАЯ ФУНКЦИЯ МЫШЦЫ: МОДЕЛЬ ГЕНЕРАЦИИ СИЛЫ ПРИ ИЗОМЕТРИЧЕСКОМ ВОЗБУЖДЕНИИ
Д.А. Черноус, С.В. Шилько
Институт механики металлополимерных систем им. В.А. Белого Национальной академии наук Беларуси, Беларусь, 246050, Гомель, ул. Кирова, 32а, e-mail: Shilko_mpri@mail.ru
Аннотация. Предложено биомеханическое описание скелетной мышцы как системы множества волокон различной длины и ориентации. В рамках этого описания модифицирована (упрощена) ранее разработанная модель сократительной функции мышечного моноволокна. В упрощенной модели не учитывается жесткость моноволокна на сжатие и вводится допущение о равномерности распределения деформаций и напряжений по длине моноволокна. Сопоставлены результаты расчета скорости свободного изотонического сокращения мышечного моноволокна по исходной и упрощенной моделям. Определены временные зависимости силы, генерируемой при изометрическом возбуждении, и распределение осевого напряжения по площади поперечного сечения мышцы в виде крученой нити заданной длины, радиуса и угла крутки.
Ключевые слова: скелетная мышца, мышечное моноволокно, крученая нить, актиномиозиновое взаимодействие, осевое напряжение, продольная деформация.
Введение
Ранее авторами [3] была разработана методика описания механического поведения мышечного моноволокна при реализации сократительной функции. Использование гипотезы «скользящих нитей» [1] и соотношений теории упругости позволяет определить распределение напряжений и деформаций по моноволокну при изотоническом и изометрическом возбуждении. Для этого составляется и решается система нелинейных дифференциальных уравнений в частных производных. Прямое и непосредственное использование данной методики при переходе от моноволокна к мышце представляется достаточно трудоемкой задачей. Поэтому при описании механического поведения мышцы приходится использовать ряд упрощений. Целью настоящей работы является адекватное описание деформирования всей скелетной мышцы на основе упрощенной модели мышечного моноволокна.
Упрощенная модель мышечного моноволокна
Известно [7], что моноволокна в скелетной мышце непрямолинейны. Данное обстоятельство обусловливает практически нулевую жесткость мышцы на сжатие и растяжение при малых деформациях. При этом градиенты напряжений в активном состоянии прямолинейного моноволокна для мышцы в целом будут пренебрежимо
© Черноус Д.А., Шилько С.В., 2008
09806267
малы. Данное обстоятельство позволяет считать распределение напряжений и деформаций в мышце однородным по длине и исследовать только временные зависимости параметров напряженно-деформированного состояния.
Нулевая жесткость на сжатие мышцы, находящейся в пассивном состоянии, позволяет при описании реализации сократительной функции учитывать только актиномиозиновое взаимодействие между нитями саркомера [6]. Данное взаимодействие описывается на основе гипотезы «скользящих нитей», которую можно сформулировать в виде двух дифференциальных уравнений
Здесь п, т - количество «тянущих» и «тормозящих» поперечных актиномиозиновых мостиков, соответственно; к1, к2 - временные константы
образования и разрыва мостика, соответственно; а - общее число мостиков на половине длины саркомера; 8 - расстояние между положениями равновесия тянущего и тормозящего мостика; ё ~ 2,2 мкм - длина саркомера; е - продольная деформация мышечного моноволокна. Для адекватного описания функционирования мышечной ткани необходимо учитывать зависимость общего числа актиномиозиновых мостиков а от текущей длины саркомера при деформировании. Данная зависимость имеет сложный вид [10]. Не нарушая общности дальнейших выкладок, будем использовать следующую аппроксимацию данной зависимости
Актиномиозиновое взаимодействие создает в моноволокне осевое напряжение
сила, генерируемая одним поперечным мостиком. Так как относительное скольжение нитей саркомера сопровождается перемещением массы, то при описании сократительной активности необходимо учитывать инерционные силы. Для осевого напряжения можно записать
Здесь к - коэффициент, характеризующий инерционность подвижных элементов саркомера. Значение к определяется экспериментально. Но приближенную оценку коэффициента к можно получить, зная длину моноволокна Ь и его плотность р. Для этого определим осевое напряжение аин, обусловленное только инертностью нитей саркомера:
Здесь т = р£Ь - масса моноволокна; ііс - ускорение центра масс. Используя допущение об однородном распределении деформации и напряжения по длине моноволокна, для ускорения центра масс можно записать:
Подставив выражения для массы и ускорения центра масс, получим приближенную оценку коэффициента к.
. , ч пё . пё . ,
п = к1(а-п - т) +-є; ті =-----------є - к2т.
28 ’ 28
(1)
а(є) = атах(1 - 4Є2).
(2)
/
а = — (п - т) . Здесь £ - площадь поперечного сечения мышечного моноволокна; /
/
а = — (п - т) + кє.
£
(3)
ь
ис = є—. с 2
к =1 Ь2р.
2
Используя формулы (1) и (2), для описания процесса сокращения мышечного моноволокна составим два дифференциальных уравнения
— ( - к ё) = к1а(е) - п(2к1 - к2) + — (а-к ё:) (к1 - к2) +—е; / / 8
П = к,
—
а(е) - 2п + (а - к ё)
пё . +—е. 28
(4)
При решении системы (4) следует учитывать, что входящее в него напряжение а обусловлено только актиномиозиновым взаимодействием. До активизации сократительной функции оно равно нулю. Если на моноволокно действует внешняя нагрузка, то процесс сокращения начнется только в тот момент, когда напряжение а уравновесит эту нагрузку.
Сопоставим результаты использования системы (4) с данными, полученными на основе более сложной модели моноволокна. Подробная модель деформирования мышечной ткани, разработанная в [3], позволяет учесть сопротивление моноволокна на сжатие и описать неоднородное распределение деформаций и напряжений по длине. В соответствии с данной моделью при описании процесса активизации сократительной функции мышечного моноволокна составляется и решается система четырех дифференциальных уравнений в частных производных
• 1 / \ — ё а'
п = к,(а- п - т) + п—и , т =-к2т - п—и , и = —,
14 28 2 28 р
1 Н , /
а = — а + Ей +------------и +—
т т —
к1 (а - п - т) + к2т + п—и + — (п - т)
8т
(5)
В системе (5) Е - мгновенный модуль упругости материала моноволокна; Н -длительный модуль упругости; т - время релаксации. Штрихами при перемещении и напряжении обозначена производная по координате х, отсчитываемой вдоль длины моноволокна. Математическую модель (5) в дальнейшем будем называть моделью вязкоупругого стержня».
В качестве примера использования систем (4) и (5) рассмотрим свободное изотоническое сокращение моноволокна при следующих значениях параметров: длина
9 2
Ь = 0,07 м; площадь поперечного сечения — = 5-10 м ; условная плотность материала р = 1000 кг / м3; время релаксации т = 1 мс. Параметры актиномиозинового взаимодействия примем следующими: к1 = 160 с-1; к2 = 88,9 с-1; 8 = 8-10-9 м; / = 3-10-12 Н; атах = 2-106. Данные значения параметров, как было показано в работе [4], соответствуют портняжной мышце лягушки при температуре 0°С.
На рис. 1 представлены временные зависимости скорости свободного
сокращения, рассчитанные на основе модели работы [3] и упрощенной системы (4). Видно, что через некоторый промежуток времени tyст ~ 50 мс после начала активации мышечного моноволокна свободное сокращение является стационарным, т.е. при t > tyст скорость сокращения изменяется со временем незначительно. При этом расчетные оценки значения ^ст, полученные на основе точной (5) и упрощенной (4) моделей, практически совпадают.
V, см/с
Рис. 1. Временная зависимость скорости свободного сокращения. Сплошные кривые -использование математической модели (5); пунктирная - использование упрощенной модели (4). Кривая 1 построена для Е = 105 Па, Н = 104 Па; кривая 2 - Е = 104 Па,
Н = 103 Па
Следует отметить, что при использовании подробной модели (5) скорость сжатия моноволокна для I > 1уст несколько уменьшается со временем. Нарушение стационарности обусловлено жесткостью прямолинейного моноволокна на сжатие. В соответствии с моделью (5) характер зависимости v(t) скорости при ^ > ^уст определяется вязкоупругими характеристиками (Е, Н, т) материала моноволокна. При использовании предложенной упрощенной модели (4) стационарность свободного сокращения не нарушается. В соответствии с системой (4) скорость свободного сокращения моноволокна в стационарном режиме определяется следующим образом:
V уст = - Ьє| = 28к
1г>густ
ь
уст и Муст 2 а •
Если сокращение моноволокна осуществляется при наличии ненулевого постоянного напряжения а0, то установившаяся скорость определяется по формуле
= 28Ьк,к 2 (/а- Ла р)
уст а(к1 + к2)(0 + /ак1 /(к1 + к2))' )
Введя обозначения:
/ак1 2§ьк1к 2 „ х ті о
а =------—, Ь =-------Р0 = /а, Р = Ла0,
к +к2 а(к1 + к2) 0 0
выражение (6) можно переписать в форме, аналогичной уравнению Хилла [8].
ЬР0 -Р
V = Ь—0-----.
Р + а
Известные экспериментальные данные [8, 9] о механическом поведении мышц свидетельствуют о неизменности скорости свободного изотонического сокращения в стационарном режиме вплоть до прекращения процесса активации. Значение скорости изотонического сокращения определяется по уравнению Хилла и зависит только от параметров актиномиозинового взаимодействия. Таким образом, предложенная упрощенная методика (4), являясь менее подробной по сравнению с моделью (4), позволяет, однако, более адекватно отразить известные особенности деформирования мышечной ткани. Важным преимуществом системы (4) является относительная простота расчетов, что существенно упрощает переход от отдельного моноволокна к скелетной мышце в целом.
При переходе от отдельного мышечного моноволокна к мышце необходимо учитывать сложный характер взаимного расположения и ориентации волокон. Для описания структуры скелетной мышцы используем функцию Ь(г), которая определяет зависимость длины моноволокна Ь от расстояния до оси симметрии мышцы. При рассмотрении мышцы как крученой нити [5] длиной Ь0, радиусом Я и углом крутки 9 функция Ь(г) будет иметь вид:
Если перед активацией сократительной функции мышце длиной Ь0 было придано удлинение и, то функция Ь(г) позволяет определить начальную продольную деформацию каждого моноволокна
При описании сокращения первоначально растянутой мышцы все моноволокна можно подразделить на две группы: 1) моноволокна, длина которых совпадает с заданной длиной мышцы Ь0 + и; 2) моноволокна, длина которые превышает длину мышцы. Волокна, относящиеся к первой группе, находится в изометрических условиях. При описании сократительной активности этих волокон в системе (4) следует положить:
В этом случае система (3) позволяет определить две функции времени: п({), а(^). При описании процесса генерации напряжения в изометрических условиях система (4) может быть решена аналитически
Описание активного состояния мышцы
(7)
0, Ь(г) > Ь0 + и.
) = е|,=0 = ( + ь0 - ь(г)), є0) = є(ґ) = Щ) = а п[=0 = а{=0
ь(г)
=0.
а(ї) = —п(Г) = а — + Се к1. Л Л
Для этих моноволокон система (4) позволяет вычислить п({) и е(^). В ходе свободного сокращения длина моноволокна уменьшается. Как только величина Ь(г )[ + еО)] сравняется с длиной мышцы Ь0 + и, моноволокно, первоначально относящееся ко второй группе, перейдет в первую.
Численная реализация методики позволяет при заданном удлинении и определить временную зависимость осевого напряжения а во всех моноволокнах мышцы. Проинтегрировав напряжение а по площади поперечного сечения мышцы, вычислим общую сжимающую силу Р, которая генерируется мышцей при активации сократительной функции.
Если длина мышцы Ь0 + и поддерживается постоянной, сила Р, обусловленная актиномиозиновым взаимодействием, возрастает от нуля до максимально возможного при данном удлинении и значения Ртах и в дальнейшем остается неизменной до завершения процесса активации. Данный режим функционирования мышцы называют изометрическим возбуждением [3]. Если поддерживается постоянной внешняя растягивающая сила Q, вызвавшая удлинение мышцы в пассивном состоянии, то генерируемая мышцей сила Р перестанет изменяться, достигнув значения Q. С этого момента начнется уменьшение общей длины мышцы. Данный процесс называют изотоническим сокращением.
Пример расчета
В качестве примера использования разработанной методики рассмотрим изометрическое возбуждение мышцы, моделируемой крученой нитью длины Ь0 = 0,1 м и радиусом Я = 1 см. Площадь поперечного сечения и параметры актиномиозинового взаимодействия моноволокон, образующих данную мышцу, те же, что и для ранее рассмотренного моноволокна (см. рис. 1). Моделируемую мышцу удлиним на и и активируем сократительную функцию в изометрических условиях.
На рис. 2 представлены расчетные временные зависимости силы Р,
генерируемой мышцей в рассматриваемых условиях. Зависимости на рис. 2 позволяют сделать вывод о том, что увеличение угла крутки 9 при малых начальных удлинениях приводит к росту промежутка времени tP, необходимого для достижения силой Р максимального значения Ртах. Само значение Ртах практически не зависит от угла крутки. При больших удлинениях (относительное удлинение более 10%) влияние угла крутки на промежуток времени tP весьма незначительно, но увеличение 9 приводит к некоторому увеличению максимальной силы Ртах.
Предложенная методика позволяет определить на только интегральную силу Р, но и осевое напряжение а в каждом моноволокне. Так, для рассматриваемой мышцы (9 = 30°, и = 10 мм) были получены распределения напряжения по площади поперечного сечения (рис. 3) на различных стадиях изометрического возбуждения. Момент времени 5 мс соответствует достижению 1/3 от максимального значения силы Р. В момент времени 100 мс увеличение силы Р практически прекращается (^ = 100 мс).
Кроме того, на рис. 3 представлено распределение по площади поперечного сечения мышцы напряжения ап, обусловленного вязкоупругим сопротивлением материала моноволокон на сжатие. Это напряжение определяется начальной растягивающей деформацией моноволокон (7) а п = Н е| . Использование в последнем
равенстве длительного модуля упругости (Н = 104 Па) связано с тем, что активизация
Рис. 2. Временные зависимости генерируемой мышцей силы Р в изометрических условиях. Числа у кривых соответствуют значению угла крутки в градусах. Сплошные кривые построены при удлинении мышцы на и = 20 мм; пунктирные - и = 1 мм
Рис. 3. Зависимость осевого напряжения а в моноволокнах от расстояния г до оси симметрии мышцы. Числа у кривых соответствуют моменту времени в мс. Пунктиром показано напряжение ап до начала активации
сократительной функции осуществляется после прекращения процесса релаксации напряжений.
Графики на рис. 3 свидетельствуют о том, что для всех моноволокон, находящихся в изометрических условиях (первая группа моноволокон), напряжение а практически мгновенно выравнивается. Процесс завершения роста силы Р
приближенно соответствует распространению одинакового напряжения а по всей площади поперечного сечения мышцы.
Заключение
Установлено, что разработанная упрощенная математическая модель
сокращения мышечного моноволокна точнее по сравнению с ранее предложенной моделью «вязкоупругого стержня», описывает известные особенности деформирования мышечных тканей. Эта модель позволила выявить влияние параметров внутренней структуры (в частности, угла крутки) на процесс генерации сжимающей силы в условиях изометрического возбуждения скелетной мышцы.
Список литературы
1. Дещерский, В.И. Математические модели мышечного сокращения / В.И. Дещерский. - М.: Наука, 1977.
2. Морецкий, А. Исследование и моделирование биомеханических свойств скелетных мышц / А. Морецкий, К. Федюлис, К. Кендзиор // Механика полимеров. - 1975. - № 4. - С. 629-633.
3. Черноус, Д.А. Вязкоупругая модель мышечного моноволокна / Д.А. Черноус // Механика. Научные исследования и учебно-методические разработки. Сб. научн. трудов. Вып. 1, Гомель: УО БелГУТ, 2007. - С. 58-64.
4. Черноус, Д.А. Моделирование сократительной активности мышечной ткани / Д.А. Черноус, С.В. Шилько // Российский журнал биомеханики. - 2006. - Том 10, № 3. - С. 53-62.
5. Черноус, Д.А. Описание диаграммы растяжения крученой нити / Д.А. Черноус, С.В. Шилько // Пластические массы. - 2006. - № 7. - С. 41-44.
6. Шестаков, Д.А. Математическая модель механических свойств волокон скелетной мышцы с учетом растяжимости актиновых нитей / Д.А. Шестаков, А.К. Цатурян // Биофизика. - 1998. - Том 43, № 2. -С. 329-334.
7. Шмерлинг, М.Д. Скелетная мышца: структурно-функциональные аспекты адаптации / М. Д. Шмерлинг, Е. Е. Филюшина, И.И. Бузуева и др. - Новосибирск: Наука, Сиб. отделение РАН, 1991.
8. Hill, A.V. First and last experiments in muscle mechanics / A.V. Hill. - Cambridge, 1970.
9. Lee, Eun-J. New insights into the passive force enhancement in skeletal muscles / Eun-J. Lee, V. Joumaa, W. Herzog // Journal of Biomechanics. - 2007. - Vol. 40. - P. 719-727.
10. Rassier, D.E. Stretch-induced, steady-state force enhancement in single skeletal muscle fibers exceeds the isometric force at optimum fiber length / D.E. Rassier, W. Herzog, J. Wakeling, D.A. Syme // Journal of Biomechanics. - 2003. - Vol. 36, No. 9. - P. 1309-1316.
ACTUATING MUSCULAR FUNCTION: A MODEL OF FORCE GENERATION AT ISOMETRIC EXCITATION
D.A. Chernous, S.V. Shilko (Gomel, Belarus)
A biomechanical description of a skeletal muscle is proposed in the form of a set of fibers of different length and orientation. This description was used to modify (simplify) the model developed previously for the contractive function of a muscle monofiber. The compressive stiffness of the monofiber is disregarded in the simplified model. An assumption of a stress and strain homogeneous distribution along the monofiber length is made. The calculation results of a free isotonic contraction velocity of a muscle monofiber of the original model are compared with the simplified one. Time dependences of the force generated at
isometric excitation have been calculated along with axial stress distribution over a crosssectional area of a muscle taken in a twisted thread form of a given length and torsion angle.
Key words: skeletal muscle, muscle monofiber, twist fiber, actinomyosin interaction, axial stress, longitudinal strain.
Получено 15 января 2008
