CC BY
25
ПРОГНОЗИРОВАНИЕ ПОВРЕЖДЕНИЙ В МЯГКИХ ТКАНЯХ СТОПЫ ЧЕЛОВЕКА ПРИ ДЕЙСТВИИ ЭКСТРЕМАЛЬНЫХ НАГРУЗОК
Р.О. Сарачев, Ю.В. Калашников, Ю.И. Няшин
Россия, 614600, Пермь, Комсомольский проспект, 29а, Пермский государственный технический университет, кафедра теоретической механики
Аннотация: В работе рассматриваются методика прогнозирования накопления повреждений в мягких тканях стопы спортсмена при приземлении на жесткую поверхность. Гимнаст, в том числе его стопа, моделируется многозвенной стержневой системой с учетом упругих межзвенных связей. Задача динамики движения позволяет найти внешние и внутренние усилия. По найденным нагрузкам, действующим на стопу, определяются напряжения в составляющих ее элементах и исследуется повреждаемость биоструктуры связок и сухожилий.
Ключевые слова: накопление повреждений, напряжения, связки, сухожилия, стопа
Введение
Глобальной целью исследований является рассмотрение процесса накопления повреждений в мягких тканях стопы спортсмена (мышцах, связках и сухожилиях) при приземлении на жесткую поверхность. Важную роль при изучении этих процессов играют напряжения, которые имеют место во всех структурах стопы. Знание напряжений имеет фундаментальное значение при понимании физиологии и патофизиологии стопы, и может представить очень ценную информацию для предотвращения и лечения повреждений. Нагрузка, действующая на структуры локомоторной системы человека при спортивных упражнениях, является одним из возможных стимулов сохранения или увеличения сопротивления биологических тканей повреждаемости. Однако, чрезмерная нагрузка может быть причиной микро- или макротравм анатомических структур.
Одно из первых рассмотрений вопроса о перенапряжении тканей у спортсменов состоялось на XII Международном конгрессе по спортивной медицине в 1958 г. в Москве. Большой доклад по этой проблеме сделал известный итальянский ученый Д. Ла Кава (D. La Cava) [1]. В нем он показал, что перенапряжение как процесс является патологическим фактором и его нельзя смешивать с физиологическим изнашиванием тканей, вызванным самой жизнью. Ла Кава высказал мысль о том, что в развитии микротравматического патологического процесса большую роль играют состояния, при которых оказывается нарушенным обмен веществ как местных тканей, так и всего организма. Наличие в организме очагов воспалительных процессов является предрасполагающим условием для перенапряжения тканей.
В некоторых исследованиях считается, что нагрузка является чрезмерной во многих спортивных упражнениях. Известно, что в течение одного года получают травмы от 25% до 50% бегунов [2], от 21% до 52 % теннисистов получают травмы опорно-двигательного аппарата в течение сезона [3]. Травмам подвергаются и спортсмены профессионалы и спортсмены любители. Так, например, известно, что среди людей, занимающихся спортивной аэробикой и получивших травму, 76% -инструкторы и 43% - их подопечные [4]. Эти цифры говорят о том, что проблема
существует и необходимо предпринимать определенные шаги к исправлению ситуации.
При анализе спортивных повреждений должны быть рассмотрены самые различные вопросы, такие как: непосредственное исследование поврежденной
анатомической структуры, тип повреждений, внешние силы, действующие на атлета при упражнениях, которые могут привести к травмам, действительные внутренние силы или напряжения, действующие в анатомической структуре до и во время повреждения, материальные свойства поврежденных структур и их критические пределы и другие. Однако одной из главных целей исследований, относящихся к спортивным повреждениям, должно быть уменьшение частоты этих повреждений.
По данным отделения спортивной травмы Центрального института травматологии и ортопедии [5, 6] частота повреждений мягких тканей по локализации распределяется следующим образом: нижние конечности - 62%; верхние конечности -22%; прочие - 16%. При этом наблюдается выраженная преимущественная зависимость между определенной группой биологических тканей и занятием определенным видом спорта.
В результате чрезмерных и однообразных нагрузок перенапряжение испытывает место перехода мышц в сухожилие и область прикрепления сухожилий к костной ткани [7]. Это нередко вызывает боли у спортсменов и невозможность продолжать тренировки.
На наш взгляд, анализ процесса накопления повреждений в мягких тканях стопы должен состоять из нескольких этапов. В качестве примера такого рода анализа спортивных движений предполагается рассмотреть приземление гимнаста на жесткую поверхность и возникающие при этом повреждения в мягких тканях стопы (связках, сухожилиях, мышцах).
Этапы анализа
1. Определение начальных условий для последующего анализа, то есть применительно к нашей задаче необходимо определить кинематические характеристики движения человека до контакта с поверхностью приземления.
2. Построение биомеханической модели движения человека и определение по заданному движению реакций, действующих на стопу.
3. Построение структурной биомеханической модели стопы, которая отражала бы все биомеханические свойства объекта исследования (взаимодействие структурных элементов, их подвижность и так далее) и определение относительных перемещений структурных элементов и напряжений в мягких тканях стопы.
4. Определение по известным напряжениям, возникающим в связках и сухожилиях стопы, возможного уровня накопления повреждений, учитывая их залечивание.
5. Идентификация факторов, влияющих на возникновение повреждений в подвергнутых опасности структурах, и разработка стратегии обеспечения уменьшения возможных повреждений.
Модель приземления
Существует несколько подходов к построению биомеханических моделей, описывающих движение человека. Использование того или иного подхода зависит от предпочтений исследователя и тех целей, которые он перед собой ставит. Тем не менее, основной принцип построения остается неизменным - моделирование частей тела
человека в виде твердых тел, соединенных между собой шарнирами и описание движения полученной системы при помощи уравнений Лагранжа II рода. Отличия могут заключаться в моделировании действия мышц и усложненном представлении некоторых элементов модели. Так, например, в работах [8, 9] действие мышц заменяется эффективными моментами в суставах, а в работе [10] вводится в рассмотрение модель мышцы с эффективным плечом мышцы, а также дополнительный элемент, описывающий рессорные свойства стопы.
Мы будем придерживаться подхода использованного в работах [8, 9]. Рассмотрим приземление гимнаста при вертикальном полете с высоты к, когда вектор скорости его центра масс перпендикулярен горизонтальной поверхности приземления, механические свойства которой известны. Тело спортсмена моделируется шестизвенной шарнирно соединенной стержневой системой: ступни, голени, бедра, торс и голова, плечи, предплечья. Г еометрические и инерционные характеристики всех звеньев известны, а правые и левые части тела действуют синхронно, образуя единый элемент.
Характер полета и относительно небольшая скорость приземления позволяют пренебречь сопротивлением воздуха. Общие теоремы динамики, эквивалентные для случая вертикального полета вытянувшегося гимнаста задаче динамики точки, определяют скорость центра масс (ступней) в момент касания поверхности по известной формуле:
которая является начальным условием для дальнейших исследований.
Движение гимнаста с момента касания матов описывается системой дифференциальных уравнений второго порядка, число которых определяется числом степеней свободы стержневой модели. Сегменты тела предполагаются недеформируемыми, суставы заменяются шарнирами без трения (идеальными). Положение спортсмена определяется восемью обобщенными координатами: хт, ут, ф;, ф2, ф3, ф4, ф5, ф6, которые представлены на рисунке 1 (хт, ут - координаты тазобедренного сустава).
(1)
D MD
Рис. 1. Стержневая модель гимнаста.
В нашей работе система уравнений получена с помощью уравнений Лагранжа II рода может быть записана в виде:
Jk (q;, qj, q}, xd ,yd , md , Mi)=° j,k =1,S’i=i,5, (2)
где qj - обобщенные координаты, XD, YD - горизонтальное и вертикальное усилия, действующие на стопу (реакция поверхности приземления), MD - приведенный к точке касания момент реакции опоры, Mt - приведенные к суставам моменты мышечных
усилий, S - число степеней свободы механической системы.
Конкретный вид уравнений (2) зависит от упрощений в шестизвенной стержневой модели спортсмена в связи с теми или иными задачами исследования движения гимнаста. Если известны все силовые воздействия (в том числе приведенные к суставам моменты мышечных усилий), то к системе уравнений (2) добавляются соответствующие начальные условия, и она решается аналитически для простых случаев или интегрируется численно. Обратная задача заключается в нахождении силовых воздействий по известной кинематике процесса приземления, которая определяется из эксперимента.
Модель стопы
Построению биомеханической модели стопы и исследованию ее механических свойств посвящены работы [11, 12]. Использование результатов этих исследований применительно к решаемой нами задаче напрямую невозможно, так как нас интересует поведение стопы в процессе движения человека, что не учитывается в предлагаемой [11, 12] модели. Тем не менее, некоторые идеи из этих работ будут использованы при построении модели стопы в нашей работе.
Стопа имеет аркообразную форму, воспринимающую нагрузку, приложенную к таранной кости со стороны большеберцовой и малоберцовой костей. Такая аркообразная форма сохраняется при действии нагрузки, что достигается за счет plantar aponeurosis - мышцы, которая простирается от выступа пяточной кости до плюснефаланговых суставов [13]. Основные опорные области стопы сконцентрированы под окончаниями (головками) пяти плюсневых костей и выступом пяточной кости. Определение распределения нагрузки для этих шести областей является объектом отдельного исследования. Отметим только один факт, что при наличии большого числа научных работ, посвященных этому вопросу, полученные результаты сильно отличаются не только у разных исследователей, но и в рамках одного исследования [14,
15, 16, 17].
Введем следующие обозначения:
O - центр вращения голеностопного сустава.
XDY - неподвижная система координат, связанная со стопой.
EOA - сегмент, при помощи которого моделируются сразу две кости - пяточная и таранная. Эти кости представлены в виде одного тела в силу их малой подвижности друг относительно друга.
AB, BC, CD - сегменты, отвечающие за ладьевидную, клиновидную и плюсневую кости стопы, соответственно.
Fi (i=1..2), T - мышечные силы, действующие в данной рассматриваемой части стопы, F¡ - силы имеющие одну точку крепления на стопе, а другую вне ее, T -отвечает за мышцу (plantar aponeurosis), оба конца которой крепятся на стопе.
YO, XO, M0 - реакции и момент в голеностопном суставе
YD, Xd , MD - реакции и момент, приложенные к головке плюсневой кости.
хг, Уг ( = 1,4) - координаты центров масс каждого из сегментов стопы.
Будем рассматривать стопу как шарнирно-соединенную систему твердых тел. Пять аналогичных друг другу продольных сводов стопы заменим одним. Схема стержневой модели стопы представлена на рисунке 2. Чтобы в дальнейшем перейти к трехмерной модели, необходимо будет рассмотреть все пять продольных сводов стопы, учесть их взаимное расположение и влияние друг на друга.
Предположим, что нам известны условия нагружения свода стопы, тогда необходимо определить:
1. какие при заданных внешних усилиях должны быть усилия в мышцах и связках,
2. как при данных условиях нагружения будет менять свою форму свод стопы.
Иными словами, мы хотим при заданных усилиях У0, Х0, Уа, ХП и моменте МП
найти углы Д (г = 1,4), определяющие конфигурацию свода стопы, усилия ^, ^2, Т и моменты М0, МА, Мв, Мс , ЯА, Яв, Яс .
Уравнения движения свода стопы можно получить аналогично уравнениям (2), но тогда будут исключены межзвенные реакции, которые необходимы для нахождения напряжений в связках. Поэтому для решения поставленной задачи рассмотрим сегменты стопы ЕА, АВ, ВС и СП по отдельности и для каждого выделенного сегмента запишем уравнения движения.
*4
>
X
Рис.3. Стержневая схема сегмента ЕА стопы.
На рисунке 3 изображен сегмент ЕА, который включает в себя пяточную и таранную кости стопы. Для этого сегмента получим следующие три уравнения: сегмент ЕА
т1Х 1 - F2x + XO + XA + Tx ,
т1 У1 - ~P\ + F2 y + YO + YA + Ty ,
I i'$i - Ma + м о + M„ + мЯа + + мт.
(3)
(4)
(5)
Здесь т
и
Mf , м
F2 RA
суммарные масса и вес пяточной и таранной костей, ^ , ±у.1Т ,±у.1К моменты, создаваемые силами ¥2, Ял , Т и Р0 относительно
центра масс сегмента ЕА.
Величины моментов МА,М0 и силы Т зависят от свойств соединительной ткани. Вопрос построения определяющих соотношений для мягких тканей (мышц, связок и сухожилий) в данной работе не рассматривается. На сегодняшний день имеется огромное количество определяющих соотношений от простых соотношений теории упругости до соотношений, в которых учитываются происходящие в мягких тканях химические процессы. На наш взгляд, выбор вида определяющих соотношений должен обуславливаться целью исследования. В первом приближении мы будем использовать определяющие соотношения теории упругости для моментов МА,М0, а
также Мв, Мс .
Величину силы Т также будем считать подчиняющейся упругому закону:
fL - О
, где Д Ь0 - текущая и начальная (в недеформированном состоянии)
'0 У
длина отрезка, соединяющего точки Е и П (рис.2), ЕТ - модуль упругости. Впоследствии предполагается использовать более сложные определяющие
соотношения, учитывающие вязкоупругие свойства тканей, например, такие как предложены в работе Фанга [18].
Аналогично тому, как было сделано для сегмента ЕА, запишем уравнения движения для оставшихся частей стопы: сегмент АВ
сегмент ВС
m 2 x: 2 II X B 1 X A (6)
m2y 2 P2 ^ YB YA , (7)
I202 — —MA + MB + MRa + MRb , (8)
m3x3 - XC — XB, (9)
тзУ3 - — P3 ^ YC — YB , (10)
Із0з - —Mb + Mc + MRb ^ Mrc , (11)
ID4 04 - -Mc + MRc + MD + MP4 . (12)
сегмент СD
Данная система уравнений (3)-(12) получена в предположении, что в данный момент времени нет контакта пяточной кости с поверхностью, в результате чего в последних трех уравнениях этой системы отсутствует сила ¥], которая отвечает за поднятие носка стопы.
Из системы уравнений (3)-(12) можно временно исключить все внутренние реакции в суставах, проекции силы ¥2 на оси X, У и ускорения центров масс с помощью
кинематических соотношений. Замкнутая система уравнений получается после добавления соотношений для координат точки О:
хо (і )= Ргх (Рі } и І = 1,4, (13)
Уо (і)= Ргу (рі } г, І = 1,4 • (14)
После чего мы получим следующую систему дифференциальных уравнений для определения изменения неизвестных углов Р г (г = 1,4):
а11Р1 ^ а12р2 ^ а13р3 + а14р4 = /1 , (15)
а22р2 ^ а23р3 ^ а24р4 = /2 , (16)
а33р3 ^ а34р4 = /з , (17)
а44р4 = /4 • (18)
где коэффициенты а у являются функциями масс-инерционных характеристик стопы и углов Р і, / - функциями заданных усилий Уо, Хо, Ув, Хв, момента Мв, углов Р і и
хо (і), уо (і). Полученная система (15)-(18) интегрируется численно при
соответствующих начальных условиях.
Математическая модель повреждаемости
Одной из причин травм в мягких тканях в медицинской литературе называется асимметричная быстроменяющаяся нагрузка. Любые мышцы состоят из пучков поперечно полосатых волокон [19]. В типичной мышце с параллельно расположенными волокнами каждое из волокон идет по всей ее длине. На обоих концах волокна прикрепляются к элементам скелета непосредственно или с помощью сухожилия.
Сухожилия построены из пучков волокон, которые состоят из коллагена. Волокна эти расположены параллельно продольной оси сухожилия, благодаря чему механические свойства коллагена приобретают решающую роль в формировании механических свойств сухожилия. Основные из этих свойств - большая прочность на разрыв, гибкость, незначительная растяжимость и эластичность.
Введем параметр, отвечающий за целостность материала / - сплошность. Пусть, сначала происходит “хрупкое разрушение” [20] одного из волокон при /=0, после чего нагрузка перераспределяется на остальные волокна, для которых 0 <у ^ < 1. Далее
происходит последовательное “разрушение” всех остальных волокон мышцы.
Примем, что скорость роста поврежденности определяется эффективным
напряжением актах/у (для сокращения записи <5к =aЫax), т.е.
5=1 N. <19>
ш V у У
Рассмотрим более простой случай степенной зависимости
dy /а ^
—= - A dt
к
(20)
V '
A>0, n>0 - некоторые постоянные. После интегрирования это уравнение дает при у/=0
tk =-----1-----, (21)
* A(n + 1)а П V 7
tk - время “разрушения” волокна при заданной нагрузке. Последовательный разрыв
волокон будет определять ступенчатый характер нагрузки (напряжения) для каждого
неразорвавшегося еще волокна. Тогда справедлив закон линейного суммирования повреждений, и можно записать (21) в виде [20]:
i At ------
Z—* = 1, i = 1, N 0, (22)
k=1 tik
где N0 - число волокон в мышце; tik - время “разрушения” i-ого волокна при действии k-ой нагрузки, то есть отвечающее текущему значению напряжения; At* - время действия k-ой нагрузки. Последнее выражение есть ни что иное, как система
алгебраических уравнений порядка N0 относительно времени последовательного
разрушения Atk действия N0 нагрузок. При этом в каждом последующем уравнении на одно неизвестное больше чем в предыдущем, а первое есть уравнение вида
At1 = t11 = Л/ 1 И n > (23)
A(n + 1)оп
здесь g.* - напряжение в волокне с номером i ( i = 1,k ) при количестве
неразорвавшихся слоев равном N0 - k ( k = 1,N0 ), k - номер последнего разорвавшегося волокна.
Поэтому возможно получение рекуррентной формулы для вычисления Ati и тогда время разрушения слоя толщиной H определяется следующим выражением
N0
t = At, A^., (24)
i=2
в котором Ati можно определить из системы (22), а для At1 будем использовать (23).
Исследование накоплений повреждений в мышцах в данной работе проводилось на основе анализа состояния этого вопроса в литературе. Все литературные источники, известные авторам работы, которые в той или иной степени касаются рассмотрения этого процесса, носят описательный характер, так как этой проблемой не занимались ни математики, ни механики. В литературе можно встретить анализ причин и механизмов возникновения, накопления повреждений, основанный на наблюдениях лечащих врачей, а также общие рекомендации более рационального проведения тренировок и выполнения различных упражнений на основе практического опыта самих спортсменов.
Предложенная в работе методика оценки повреждений в мышцах позволяет провести сравнительный анализ спортивного упражнения с разными уровнями нагрузок по времени условного разрыва волокон мягких тканей и определение наиболее безопасного выполнения данного упражнения. Количественная оценка времени накопления повреждений возможна после получения из эксперимента констант уравнения повреждаемости (20), что выходит за рамки данной работы.
Численная реализация и обсуждение
Численная реализация предложенной методики проводилась для спортсмена с данными, приведенными в работе [8]. Не имея литературных данных или собственного эксперимента по динамике приземления, авторы провели численный эксперимент при отталкивании гимнаста от жесткой поверхности от полного касания ступни до отрыва пятки на некоторую высоту. Задав простую кинематику, близкую к реальной, получили значение усилий Уа, Х0, Ув, Хв и момента Мв. Координаты х0, у0 точки О входили в заданную кинематику.
Модель нагружения стопы и определение напряжений в мягких тканях рассчитывалась как в статическом, так и динамическом вариантах при начальных условиях полученных выше. Были использованы значения материальных констант приведенные в [11]. Геометрические размеры стопы снимались с натурной модели скелета стопы. Полученные усилия и напряжения близки к литературным экспериментальным [17] и расчетным [11] данным.
Литература
1. ЛА КАВА Д. Хронические повреждения в спорте как результат повторных микротравм. В книге: Спортивная медицина. Труды XII Международного конгресса. Москва, 400-408, 1959.
2. NOVACHECK T.F. Running injuries: A biomechanical approach. The Journal of Bone and Joint Surgery: 1220-1233, 1999.
3. NIGG B.M., BOBBERT M. On the potential of various approaches in load analysis to reduce the frequency of sports injuries. J Biomechanics, 25: 3-12, 1990.
4. RICHIE D.H., KELSON S.F., BELLUCCI P.A. Aerobic dance injuries: a retrospective study of instructors and participants. Phisician Sports Med, 13: 130-140, 1985.
5. Труды Центрального института травматологии и ортопедии. Выпуск 30: Повреждения и заболевания опорно-двигательного аппарата, 1986.
6. Труды Центрального института травматологии и ортопедии. Выпуск 31: Диагностика и лечение повреждений и заболеваний опорно-двигательного аппарата, 1987.
7. МИРОНОВА З.С., МЕРКУЛОВА Р.И., БОГУЦКАЯ Е.В., БАДНИН И.А. Перенапряжения опорно-двигательного аппарата у спортсменов. Москва, Физкультура и спорт, 1982.
8. ВОРОНОВ А.В., ЛАВРОВСКИЙ Э.К. О моделировании рациональных вариантов техники бега на коньках. Современные проблемы биомеханики. Выпуск 7: 144-163, 1993.
9. KALACHNIKOV Y.V., NYACHIN Y.I., SARATCHEV R.O. Dynamics of a gymnast touchdown on elasto-plastic mats. Russian Journal of Biomechanics, 2(1-2): 81-88, 1998.
10. WANG Y., MATSUHISA H., HONDA Y. Loads on Lower Joints and Optimal Action of Muscles for Shock Reduction. JSME International Journal, 42(3): 574-582, 1999.
11. SALATHE E.P. JR, ARANGIO G.A., SALATHE E.P. A biomechanical model of the foot. J Biomechanics, 19(12): 989-1001,1986.
12. SALATHE E.P. JR, ARANGIO G.A., SALATHE E.P. The foot as shock absorber. J Biomechanics, 23(7): 655-659, 1990.
13. HICKS J.H. The foot as a support. Acta Anat, 25: 34-45, 1955.
14. BETTS R.P., FRANKS C.I., DUCKWORTH T., Burke J. Static and dynamic foot-pressure measurements in clinical orthopaedics. Medical & Biological Engineering & Computing, 18: 674-684, 1980.
15. POLLARD J.P., Le QUESHE L.P., TAPPIN J.W. Forces under the foot. J Biomedical Engineering, 5: 37-40, 1983.
16. SIMKIN A., F. STOKES I.A. Characterisation of the dynamic vertical force distribution under the foot. Medical & Biological Engineering & Computing, 20: 12-18, 1982.
17. SOAMES R.W. Foot pressure patterns during gait. J Biomedical Engineering, 7: 120-126, 1985.
18. FUNG Y.-C. Stress-strain-history relations of soft tissue in simple elongation. In Biomechanics: Its Foundations and Objectives, Prentice Hall, Englewood Cliffs, New Jersey, 181-208, 1972. (Фанг Я.-Ч. Математические модели зависимости напряжение - деформация для живых мягких тканей. Механика полимеров, 5: 850 - 867, 1975).
19. АЛЕКСАНДЕР Р. Биомеханика. Москва: Мир, 1970.
20. КАЧАНОВ Л.М. Основы механики разрушения. Москва: Наука,1974.
INJURIES PREDICTION IN THE SOFT TISSUES OF THE MAN FOOT DURING THE EXTREME LOADS ACTION
R.O. Saratchev, Y.V. Kalachnikov, Y.I. Nyashin (Perm, Russia)
In this paper, the method has been proposed for the injuries accumulation prediction in the soft tissues of the sportsman foot during a vertical landing on a hard surface. The gymnast, including his foot, was simulated by the rigid bodies system, the elastic relationships between rigid bodies were considered. The external and internal forces on the foot were determined from the movement dynamics problem. These forces applied on the foot were used for the stress determination in elements included in the foot biomechanical model to diagnose injuries in the ligaments and tendons.
Key words: accumulation of injuries, stress, ligament, tendon, foot
Получено 24 мая 2000
