Механические модели в стабилометрии Текст научной статьи по специальности «Медицинские технологии»

ISSN 1812-5123. Российский журнал биомеханики. 2014. Т. 18, № 2: 184-193

УДК 531.391:612.76

Российский

Журнал

Биомеханики

www.biomech.ru

МЕХАНИЧЕСКИЕ МОДЕЛИ В СТАБИЛОМЕТРИИ

П.А.Кручинин

Механико-математический факультет Московского государственного университета имени М.В. Ломоносова, Россия, 119992, Москва, Ленинские горы, 1, е-mail: [email protected]

Аннотация. Приводится описание упрощенных моделей стабилометрических измерений. Для описания движений человека приняты традиционные модели. Движение в сагиттальной плоскости описывается моделью перевернутого маятника. Движение во фронтальной плоскости описывается упрощенной моделью трехзвенника, два звена которого параллельны и опираются на пол. Эти звенья моделируют ноги, а третье звено моделирует туловище, руки и голову, которые не меняют взаимного расположения. Звенья соединены шарнирно в тазобедренных суставах. Предполагается, что плита стабилоанализатора (стабилографа) представляет собой недеформируемое твердое тело. Эти модели используются для выражения показаний стабилоанализатора через углы наклона тела и моменты в голеностопном и тазобедренном суставах. Их приближенный анализ показывает, что изменения сагиттальной и фронтальной координат центра давления в наибольшей степени характеризуют изменение моментов в голеностопном и тазобедренном суставах соответственно. Полученная оценка справедлива для произвольных по темпу движений человека и позволяет объяснить «ухудшение устойчивости» у некоторых больных в процессе реабилитации. Предполагается, что такой эффект связан с тем, что в начале реабилитационного периода у указанных пациентов ослаблена мышечный система, и нервная система осуществляет стабилизацию вертикальной позы меньшими усилиями. В ходе реабилитации мышечная система восстанавливается, и пациент может затрачивать большие усилия на стабилизацию вертикальной позы.

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

Одним из основных инструментов постурологии служит стабилометрическая платформа. Такая платформа имеет датчики силы, которые измеряют только вертикальные компоненты усилий. Вместе с компьютерным комплексом преобразования и обработки сигналов такая платформа образует стабилоанализатор (стабилограф).

Механические модели удержания человеком вертикальной позы строились уже на самых ранних стадиях развития стабилометрии. В литературе известна интерпретация результатов исследования удержания человеком вертикальной позы на основе модели перевернутого маятника [3, 18]. Для описания фронтальных движений корректнее использовать более сложные модели. Обзор различных задач удержания человеком вертикальной позы на основании механических моделей приведен в диссертациях [4] и [14]. Подобные модели используются для обсуждения результатов стабилометрических исследований. Изложение вопроса о связи этих моделей

© Кручинин П. А., 2014

Кручинин Павел Анатольевич, к.ф.-м.н., доцент кафедры прикладной механики и управления, Москва

Введение

с показаниями стабилоанализатора часто бывает недостаточно аккуратным и содержит неточности. Необходимую аккуратность выводов позволяет сделать механический анализ, основанный на исследовании физической модели. Подобный анализ частично проведен в статье [16], однако указанная работа практически недоступна и известна в основном по ссылкам. В то же время популярные издания, посвященные анализу стабилометрических исследований [2, 13, 15], не рассматривают этот вопрос, что в ряде случаев затрудняет понимание и интерпретацию результатов стабилометрических исследований. В связи с вышеизложенным представляется полезным подробное изложение и анализ соотношений указанной механической модели.

Не описывая подробно принципов работы стабилометрической платформы, укажем, что основной информацией, которая выдается программным обеспечением стабилоанализатора, являются фронтальная х^ и сагиттальная уц координаты точки Оп приложения равнодействующей N сил вертикальных реакций, действующих в опорах платформы, а также величина этой равнодействующей.

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

Модель движений в сагиттальной плоскости

Задача описания движений человека в сагиттальной плоскости рассматривалась многократно. Однако при ее анализе часто опускают некоторые особенности, связанные с пренебрежением при моделировании малыми параметрами. В связи с этим обсудим здесь аккуратное обоснование этой традиционной модели.

Для описания колебаний человека в сагиттальной плоскости примем традиционную модель перевернутого маятника, изображенную на рис. 1. Традиционно предполагаем, что тело человека в ходе теста допустимо моделировать недеформируемым стержнем массой т, закрепленным шарнирно в точке О, соответствующей голеностопному суставу. В более сложном случае антропоморфного многозвенника такая модель будет описывать наиболее значимую - первую - форму колебаний [9, 14].

Центр масс стержня расположен в точке С, удаленной от точки О на расстояние 1с. Момент инерции стержня относительно фронтальной оси, проходящей через точку О, равен Л. Отклонение стержня от вертикали опишем углом 0. Будем считать, что обследуемый ориентирован так, что сагиттальная плоскость человека параллельна оси чувствительности платформы, а его стопа неподвижна относительно платформы.

Уравнения движения для малых значений угла 0 и скорости его изменения запишем из уравнений плоскопараллельного движения тела [1] в соответствии со схемой рис. 1, б:

\ml~00 = -Яу,

0 = К-mg, (1)

Л0 = т1 с^0 + Мх.

Третье уравнение этой системы полностью описывает движение тела человека и традиционно используется для описания движения:

0=тЬк 0(2)

Л Л

Рис. 1. Модель движений человека в сагиттальной плоскости. Внешние силы, действующие: а - на систему человек - платформа; б - на модель стержня, имитирующего тело человека; в - на систему стопы ног - платформа

стабилоанализатора

Для записи соотношений для сагиттальной стабилограммы рассмотрим равновесие системы платформа - стопа под действием силы тяжести mpg этой системы, составляющих N и ¥у равнодействующей сил в опорах, реакций Ку и Я2 связи в голеностопном суставе и момента Мх мышц в голеностопном суставе, приложенных как показано на рис. 1, в.

Уравнения равновесия этой системы имеют вид

' р> - Ку = 0

< N - Яг - т^ = 0, (3)

,Nп + ^И - Мх - mpgУp = 0,

где уп - сагиттальная координата центра давления; И - расстояние от точки О до плоскости, образованной чувствительными элементами сенсоров; ур - координата центра масс объединенной системы.

Выражение для измерений стабилоанализатора получим, преобразовав уравнения систем (1) и (3):

^ , , у ^ т1 И Л т2 я12И ^

^п = Мх - ¥уИ + mpgУp = Мх 11 + I--у- 9 + mpgУp.

Для грубой оценки предполагаем, что момент инерции J можно приближенно считать равным моменту инерции стержня длиной 2lc, что соответствует известным антропометрическим данным [6] :

Nyn = Mx

1 +

3h 1 mgh 4L ~

6 + mvgyv •

С J

Таким образом, переменная составляющая сагиттальной стабилограммы совпадает с переменной составляющей величины

У, =■

M,

f

mg

Г1+mл

V

m

1 +

3h

Л

41

h

c J

fi+mл

-L6 •

(4)

J

V

m

/

J

В этом соотношении /с0 - расстояние вдоль оси у от оси голеностопного сустава

до проекции центра масс человека.

Чтобы сравнить слагаемые в соответствии с теорией размерностей [8, 12],

приведем это соотношение к безразмерному виду, перейдя к переменным 5М = ,

Мх*

у 0

5у = — и 50 = —. В качестве характерных значений 0*, у* и Мх* примем величину

у* 0*

амплитуды изменения угла 0, отклонение у* центра давления и модуль момента в голеностопном суставе при равновесии системы 0 = 0*:

у* = / с0* и М* = mg/c0*.

Для этих переменных выражение (4) примет вид

И_

Т

л 3h

1 + — 41

Sy =-c_ SM,--

i+mp

m

fi+mp1

se •

V

m

J

Учитывая соотношения, бесспорно выполняющиеся для взрослых обследуемых, т И

—— << 1 и — << 1 и отбросив постоянные слагаемые, получим основное соотношение

т К

5у = 5Мх, которое для размерной переменной составляющей сагиттальной стабилограммы эквивалентно выражению

АМх

Ау,

mg

(5)

где AMx - изменение момента в голеностопном суставе.

Э I т— И I

Это соотношение выполняется с погрешностью порядка тах <——, — > для

Iт К)

любых по темпу движений. Таким образом, сагиттальная стабилограмма в большей степени описывает нормированное изменение момента в голеностопном суставе. Погрешность этого представления для человека среднего телосложения не превышает 10%.

Модель движений во фронтальной плоскости

Как отмечалось еще в работе [3], простая модель перевернутого маятника не описывает полностью удержание человеком вертикальной позы во фронтальной плоскости. Наличие двух опор во фронтальной плоскости и незначительные величины фронтальных моментов, развиваемых мышцами в голеностопном суставе, существенно меняют механическую схему, которую необходимо использовать при моделировании движения. Эта схема зависит от постановки ног [4, 5]. При постановке единственной ноги на платформу стабилоанализатора изменение фронтальной стабилограммы, аналогично модели измерений в сагиттальной плоскости, описывает фронтальный момент в голеностопном суставе [4].

При установке обеих ног на платформу стабилоанализатора задача существенно усложняется. Для упрощения рассмотрим только один случай движения человека, когда ступни располагаются под тазобедренными суставами. Будем считать, что опора ноги о плоскость точечная и центры тазобедренных суставов А1, А2 и точки опоры ног О1, О2 образуют параллелограмм, как показано на рис. 2. Такое положение приближенно соответствует позе «ноги на ширине плеч».

Пусть М1 и М2 - моменты, создаваемые мышцами в правом и левом голеностопных суставах соответственно; I - «длина ноги» ОА; mъ - масса туловища с головой и руками; mf - масса одной ноги, т.е. масса человека

m = mь + 2mf. (6)

Примем гипотезу из работы [4] о том, что удержание позы во фронтальной плоскости осуществляется за счет мышц тазобедренного сустава. В качестве обобщенной координаты выберем угол ф, образованный ногой ОАг с вертикалью.

R,

А

о.

х„

*JV R

* О,

m,,S

F, tff,

VR.

Рис. 2. Модель движений человека во фронтальной плоскости. Активные силы, действующие в системе человек - платформа (а), и внешние силы, действующие

на платформу стабилоанализатора (б)

б

а

Уравнения движения такой системы запишем с помощью уравнения Лагранжа [1]. Кинетическая энергия системы имеет вид

т=2 ^+ 2 2

где ^ - момент инерции ноги относительно сагиттальной оси, проходящей через точку О а - скорость центра масс «туловища».

Используем для вычисления моментов инерции ^ соотношения для моментов инерции стержня. Тогда для малых значений угла ф

T =

( + mл

3 2

l (p .

Обобщенную силу в этих же предположениях запишем в виде

0 = (шъ + шг ) glф + Мх + М2.

Тогда уравнение движения системы с использованием уравнения Лагранжа

запишем в виде

mb + mf g . M + M 2 Ф = b „f T Ф + 7-\ \ . (7)

+ 2 l ( 2 .. mh + — mf | mh + —mf 11

2

3 I 3

Это уравнение приближенно описывает движение стоящего человека во фронтальной плоскости и имеет тот же вид, что и уравнение (2), описывающее движение в сагиттальной плоскости.

Для выражения фронтальной составляющей стабилограммы потребуются данные о горизонтальных Ях1, Ях2 и вертикальных Яг1, Я22 составляющих реакции

опорной поверхности. Для их отыскания воспользуемся теоремой о движении центра масс системы, показанной на рис. 2, а, в проекциях на горизонтальную и вертикальную оси:

/ + / + ШЪ*Ъ = - Кх1 - Кх 2 ,

о = Rzl + R2 -(mh + 2m/)g,

и теоремой об изменении кинетического момента [1] этой системы относительно точки 01:

| (( + К/ 1 + Кг 2 ) Мт8 + £ М,

где у1 и у2 - горизонтальные составляющие скорости центров масс ног; \ъ - скорость поступательного движения, совершаемого туловищем; Къ, К/1, К,2 - кинетические

моменты каждого из тел системы относительно фронтальной оси, проходящей через точку 01; ^ М, ^ Мя - суммы моментов сил тяжести и реакций связей

относительно фронтальной оси, проходящей через точку 01.

С учетом малости угла ф, скорости и ускорения его изменения запишем эти соотношения в линеаризованном виде:

( + mf ) = -Rxi -Rx2, 0 = Rzl + Rz2 -(ть + 2тf )g,

(8)

т

(/ + hb) + -3mf/ I/ф = -(ть + 2mf )b + (ть + mf ))ф + 2Rz2Ь,

где Иъ - высота центра масс «туловища» над тазобедренными суставами; Ь -расстояние 0О2.

Из двух последних уравнений системы (8) и уравнения (7) выразим вертикальные реакции опорных поверхностей в виде

(ть + 2mf )g (р f 2 ,2 ^ / Xg/ф

Rz1 = --—---—I ть/(/ + h ) + -тЛ \ + (mh+mf))-,

z1 2 2ь i 3 f I [ ь f' 2ь

(ть + 2mf )g ф f 2 ,2 Xg/ф

Rz 2 = "-+ i ть/(/ + hь ) + 3 т/ 1-(ь+тг )

(9)

20

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

Уравнения равновесия имеют вид

F - R 1 - R 2 = 0,

x x1 x2 "

N - R1 - R2 - m g = 0,

z1 z 2 p<J ?

Nx + Fh - RA-R2 (A + 2Ъ)- m gx = 0,

n x p z1 z 2\ / pO p 5

(10)

где хп - фронтальная координата центра давления; ¥х - проекция реакций опор платформы на горизонтальную ось фронтальной плоскости; Ьр - расстояние от точки Ог до плоскости, образованной чувствительными элементами сенсоров платформы; А - фронтальная координата точки О1 в осях стабилоанализатора (на рис. 2, б изображено отрицательное значение параметра А).

Используем эти соотношения и выражения (8)-(9) для того, чтобы получить выражение для измерений стабилоанализатора в виде

2

Nxn =(ть + 2mf )g(ь + A) + mpgxp -(ть + mf )ф + ть ( + hp ) + 3 mf/ + (ть + mf )p

/ф .

Выразим угловое ускорение из (7) и нормальную реакцию с помощью вторых уравнений систем (8) и (10), а также используем соотношение для масс (6). Отбросив постоянные слагаемые, получим соотношение для величины хя, изменение которой описывает переменную составляющую фронтальной стабилограммы:

т - — 11 + (т - 2т/) + (т -тг ")Ир

AM,.

xs =

m/g

4

т — т.

4 mf | + (m - 2mf )) + (m - mf )hp

m -

4

m — m.

(11)

1 -

m.

m

/ф,

где ДМ- - изменение суммарного момента в тазобедренных суставах, /ф - фронтальное перемещение центра масс. Выражение (11) для

С mf ^ 1--f

V m J

фронтальной составляющей содержит «хороший» малый параметр -у << 1 и шг

отношение —— « 0,2, которое также в первом приближении можно считать малым. т

-

Пренебрегая этими величинами и учитывая, что —Ь < 0,5, получим грубую оценку для величины х, в виде

X, Г1 + -ф. (12)

Чтобы сравнить слагаемые в соответствии с теорией размерностей [8, 12],

приведем это соотношение к безразмерному виду, перейдя к переменным ЪЫ = ,

М4

Ъх, = и 5ф = —. В качестве характерных значений ф* и Ы* примем величину /ф* ф*

амплитуды изменения угла ф и модуль момента в тазобедренном суставе, которое обеспечивает равновесие системы при ф = ф*.

Ы* = (т - тг ) g/ф* « mg/ф*.

Для этих переменных выражение (12) примет вид

Ъх, « ЪЫ [1 + —■^ + 5ф .

-

Поскольку -у- < 0,5, вклад изменения момента в тазобедренном суставе по

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

Обсуждение

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

Несмотря на значительные различия в структуре механической схемы для сагиттальных и фронтальных движений, приближенные уравнения (2) и (7), описывающие эти движения, однотипны. Постоянные времени для этих уравнений

mb + 3 mf l

имеют вид Т, = и Тг =.

]|3g ^ т + тг g

В соответствии с данными [6] для человека ростом 1,7 м эти постоянные времени близки по величине и приближенно составляют 0,16 и 0,2 с. Таким образом, уравнения, описывающие удержание человеком вертикальной позы, совпадают с уравнениями

колебаний активно стабилизируемого перевернутого маятника с небольшой асимметрией свойств во фронтальном и сагиттальном направлениях.

Соотношения (4) и (11) для изменения сагиттальной и фронтальной стабилограмм совпадают по форме. Однако весовой вклад слагаемых, определяемых положением центра масс и моментом в суставе, в них различен. Если для сагиттальных колебаний показания стабилоанализатора в соответствии с (5) хорошо отображают изменение момента в голеностопном суставе, то для фронтальной плоскости подобное утверждение носит сильно приближенный характер. Тем не менее в первом приближении стабилограмму целесообразно рассматривать как измерение нормированной величины моментов. Погрешности такого представления вряд ли сильно превышают погрешности определения центра масс человека в результате приближенных громоздких процедур, описанных в [17]. Отметим, что приведенная оценка справедлива для произвольных движений человека, в то время как оценки для движения центра масс часто справедливы только для медленных составляющих, характерных для спокойного стояния.

Такой подход позволяет объяснить, например, «ухудшение устойчивости» у некоторых больных в процессе реабилитации и части космонавтов в ходе послеполетной реабилитации (см., например, [10, 11]). С этим же фактором связаны проблемы использования стабилоанализатора при анализе движений [7].

Можно предположить, что такой эффект связан с тем, что в начале реабилитационного периода у указанных пациентов ослаблена мышечная система, а нервная система, хотя тоже возможно «разбалансирована», но адаптируется быстрее и в итоге осуществляет стабилизацию вертикальной позы меньшими усилиями. В ходе реабилитации мышечная система восстанавливается, и пациент может позволить затрачивать большие усилия на стабилизацию вертикальной позы, так как минимизация суммарных моментов в суставах, в соответствии с гипотезой из работы [3], не является главной целью системы управления вертикальной позой человека.

Благодарности

Работа выполнена при поддержке Российского фонда фундаментальных исследований (грант 12-01-00839).

Список литературы

1. Бухгольц Н.Н. Основной курс теоретической механики. - М.: Наука, 1972. - Ч. II. - 332 с.

2. Гаже П.М., Вебер Б. Постурология. Регуляция и нарушения равновесия тела человека / Санкт-Петербург. мед. акад. последипломного образования. - СПб., 2008. - 314 с.

3. Гурфинкель В.С., Коц Я.М., Шик М.Л. Регуляция позы человека. - М.: Наука, 1965. - 256 с.

4. Денискина Н.В. Фронтальная устойчивость вертикальной позы человека: дис. ... канд. биол. наук. - М., 2009. - 122 с.

5. Денискина Н.В., Левик Ю.С., Гурфинкель В.С. Сравнительная роль мышц голеностопного и бедренного суставов в регуляции позы человека во фронтальной плоскости при стоянии // Физиология человека. - 2001. - Т. 27, № 3. - С. 66-70.

6. Зациорский В.М., Аруин А. С., Селуянов В.Н. Биомеханика двигательного аппарата человека. - М.: Физкультура и спорт, 1981. - 143 с.

7. Карими М.Т. Оценка достоверности временных и основанных на колебаниях центра давления характеристик при анализе устойчивости стояния во время выполнения различных упражнений // Российский журнал биомеханики. - 2010. - Т. 14, № 3 (49). - С. 79-84.

8. Новожилов И.В. Фракционный анализ. - М.: Изд-во МГУ, 1995. - 224 с.

9. Новожилов И.В., Терехов А.В., Забелин А.В., Левик Ю.С., Шлыков В.Ю., Казенников О.В. Трехзвенная математическая модель для задачи стабилизации вертикальной позы человека // Математическое моделирование движений человека в норме и при некоторых видах патологии. -М.: Изд-во МГУ, 2004. - С. 7-20.

10. Пурахин Ю.Н., Какурин Л.И., Георгиевский B.C., Петухов Б.Н., Михайлов В.М. Регуляция вертикальной позы после полета на кораблях «Союз-6 - Союз-8» и 120-суточной гипокинезии // Космическая биология и медицина. - 1972. - Т. 6, № 6. - С. 47-53.

11. Савельев М.Ю. Физиологическое обоснование стабилометрии в оценке статического равновесия у детей младшего школьного возраста в норме и при нарушениях двигательной функции: автореф. дис. ... д-ра мед. наук. - Архангельск, 2005. - 19 с.

12. Седов Л.И. Методы подобия и размерности в механике. - М.: Наука, 1977. - 440 с.

13. Скворцов Д.В. Клинический анализ движений. Стабилометрия. - М.: Антидор, 2000. - 188 с.

14. Терехов А.В. Математическое моделирование регуляции позы человека: дис. ... канд. физ.-мат. наук. - М., 2007. - 147 с.

15. Усачев В.И. Стабилометрические параметры. - Таганрог: Ритм, 2011. - 29 с.

16. Gurfinkel E.V. Physical foundations of stabilography // Agressologie. - 1973. - Vol. 14 (Spec. No. C). -P. 9-13.

17. Lafond D., Duarte M., Prince F. Comparison of three methods to estimate the center of mass during balance assessment // Journal of Biomechanics. - 2004. - Vol. 37. - P. 1421-1426.

18. Nashner L.M. Sensory feedback in human posture control: science degree thesis. - Cambrige, 1970. -198 p.

MECHANICAL MODELS OF STABILOMETRY

P.A. Kruchinin (Moscow, Russia)

The paper is devoted to the simplified models of force plate measurements. The traditional models are used to describe the motion of a person. Human motion in the sagittal plane is described by the inverted pendulum model. Simplified three link model is used to frontal plane motion description. Two links of this model are parallel and support on the floor. These units model the legs, and the third link simulates the body, arms and head, which do not change the mutual arrangement. Links pivotally connected at the hip joints. It is assumed that the force plate (stabilograph) is a rigid body. These models are used to express force plate measurements through body model angles and ankle and hip torques. The approximate analysis shows that changes in the sagittal and frontal coordinates of the centre of pressure mostly characterize the change of ankle and hip torques, respectively. The obtained estimation is valid for arbitrary tempo of human movements, and it allows us to explain "deterioration of stability" for some patients during the rehabilitation. It is assumed that this effect is associated with the fact that these patients have weakened the muscular system in the beginning of the rehabilitation period, and the nervous system performs stabilization of vertical posture with fewer efforts. The muscular system is restored during the rehabilitation, and the patient can expend greater efforts for stabilization of the vertical posture.

Key words: description of human movements, postural regulation model, stabilometry, ankle and hip joints, leg shape correction.

Получено 17 марта 2014