УДК 531/534:57+612.7
А.В. Александров*, А.А. Фролов*, Ж. Масьон**
* Институт высшей нервной деятельности и нейрофизиологии РАН, ул. Бутлерова, 5А, 117865, Москва, Россия, e-mail: labdo@mail.ru
** Университет Прованса, 13621, Экс ан Прованс, Франция
СТРАТЕГИЯ ПОДДЕРЖАНИЯ РАВНОВЕСИЯ ЧЕЛОВЕКА ПРИ НАКЛОНЕ КОРПУСА ВПЕРЕД НА УЗКОЙ ОПОРЕ
Аннотация. Две поведенческие цели одновременно решаются при наклоне корпуса человеком: наклон как таковой и поддержание равновесия. С помощью
биомеханической модели исследуется, как эти цели достигаются. Исследование основано на разложении движения по собственным векторам динамического уравнения. Теоретически найдена координация между позой и движением, при которой наклон корпуса может быть осуществлен на бесконечно узкой опоре («оптимальный наклон» без перемещения центра давления). Оптимальное движение сравнивается с экспериментально наблюдаемым. Предлагается гипотеза, что движения вдоль собственных векторов являются независимо управляемыми, но координируемыми основной и позной компонентами движения. Основная компонента обеспечивает наклон как таковой, а позная - упреждающим образом компенсирует возмущение равновесия, вызванное основной компонентой. Результаты моделирования позволяют объяснить модель электромиографической активности, наблюдаемую при наклонах корпуса вперед.
Ключевые слова: двигательная координация, поддержание равновесия у человека, синергия, моторный контроль, поза и движение, многосуставное движение, компьютерное моделирование.
Введение
Наклон корпуса у вертикально стоящего человека относится к классу движений, при которых кроме решения основной кинематической задачи (наклон как таковой) важное значение имеет и одновременное решение второй, функционально отличной задачи сохранения равновесия. Г ипотеза о том, что в любом двигательном акте можно выделить две составляющих, основную, направленную на достижение цели, и позную, связанную с необходимостью сохранения равновесия, была выдвинута в работах [1, 2]. Дальнейшие исследования показали, что позная компонента часто начинается раньше основной, упреждающим образом подготавливая условия, при которых возмущение равновесия, вызванное ожидаемым основным движением, будет минимизировано [3, 4].
Координация между позой и движением исследовалась у человека в основном на примерах, в которых целью движения было перемещение рук при стабилизации корпуса. Однако проблема поддержания равновесия значительно усложняется, когда целью движения становится не стабилизация, а перемещение корпуса. При взаимном
© А.В. Александров, А.А. Фролов, Ж. Масьон, 2002
перемещении массивных звеньев корпуса разделение движения на основную и позную компоненты не очевидно, как при движении рук. Известно, что наклон корпуса сопровождается противоположным наклоном ног [5]. В этом случае было предложено считать, что изменение угла в тазобедренном суставе есть целенаправленное движение, управляемое «основным» контроллером, а противоположное «непроизвольное» вращение в голеностопном суставе есть сопутствующая позная компонента [6]. Однако такому разделению движения на основную и позную компоненты противоречат известные данные электромиографического (ЭМГ) анализа активации аксиальной мускулатуры при наклонах корпуса у человека. При наклоне корпуса вперед ноги отклоняются назад. Если предположить, что такое отклонение ног есть упреждающая позная компонента, то движение должно начинаться с активации экстензоров голеностопного сустава. Однако реально оно начинается с активации его флексоров [7, 8]. Таким образом, остается неясным, можно ли выделить позную компоненту при организации двигательного управления во время наклона корпуса у человека, и если можно, то что именно относится к основной, а что к позной компонентам, и как осуществляется их координация?
Основной целью данной работы является теоретическое исследование двигательной стратегии (координации между позой и движением), позволяющей решить проблему сохранения равновесия при выполнении человеком наклонов корпуса вперед. Для этого теоретически проанализирована возможность выполнения движения при усложнении условий сохранения равновесия путем неограниченного уменьшения размеров опоры. Полученные результаты численного эксперимента позволяют объяснить с биомеханической точки зрения наблюдаемую экспериментально кинематическую модель движения во время наклона корпуса у человека, а также объяснить наблюдаемую при этом модель ЭМГ.
Трехзвенная биомеханическая модель
В работах, посвященных исследованию биомеханики тела человека в сагиттальной плоскости, достаточно широко используется модель тела в виде перевернутого маятника из трех жестких звеньев, вращающихся в трех цилиндрических шарнирных суставах [9-12]. При этом верхнее звено объединяет корпус, руки и голову, средним звеном являются бедра, а нижним звеном - голени. Тремя суставами являются тазобедренный, коленный и голеностопный суставы. В настоящей работе использовалась такая же трехзвенная модель.
Три суставных угла были выбраны в качестве независимых переменных, описывающих кинематику движения тела человека в рамках модели трехзвенника: угол в голеностопном (pA (A - Ankle), коленном (pK (K - Knee) и в тазобедренном (pH (H - Hip) суставах (рис. 1).
Перевернутый трехзвенный маятник предполагается стоящим на стопе, которая в сагиттальном сечении имеет треугольную форму [13]. Проскальзывание между стопой и опорой исключается. Шарниры, соединяющие звенья модели, считаются идеальными, т.е. безфрикционными, а звенья рассматриваются как твердые тела. Движение модели рассматривается в сагиттальной плоскости без отрыва точек опоры. Кроме того, принято допущение, что центры масс звеньев расположены на линиях, соединяющих соответствующие шарниры.
Тазобедренный сустав (Нір)
Голеностопны
Коленный сустав (Knee)
сустав (Ankle)
X
>
Рис. 1. Биомеханическая трехзвенная модель тела человека в сагиттальной плоскости.
ФА, ф к и ф н - суставные углы
Уравнение движения в поле сил тяжести для системы из N звеньев, соединенных последовательно в цепочку с помощью N цилиндрических шарниров, приведено в работе [9] в терминах углов наклона звеньев относительно вертикали. В случае системы из трех звеньев с тремя цилиндрическими шарнирами уравнение движения имеет вид:
где 0 является вектором углов наклона звеньев относительно вертикали, Т есть вектор суставных моментов, С' и Б' являются соответственно инерционной и гравитационной матрицами, А' является вектором, определяющим центробежные и кориолисовы силы,
Для перехода от записи уравнения движения в терминах углов наклона 0 к записи в терминах суставных углов ф необходимо сделать замену переменных 0 = Цф,
где С(ф) = ЦтС'(0)и, Б(ф) = ЦтБ'(0)Ц и А(ф, ф) = ЦтА'(0,0). Здесь и ниже первая, вторая и третья компоненты каждого из векторов относятся соответственно к голеностопному, коленному и тазобедренному суставам.
В уравнение движения входят антропометрические параметры: длины звеньев, их массы, расстояния между их центрами масс относительно нижележащего сустава, а также их моменты инерции относительно центров масс. Для оценки данных параметров использовались усредненные «стандартные» антропометрические коэффициенты, приводимые в антропометрических таблицах [14].
Выражения для элементов матриц С и Б и вектора А, а также для перемещений центра тяжести и центра давления приведены в Приложении.
С'(0)0 - D'(0)0 + A'(0,0) = NT ,
' І - І 0^ а N= 0 І - І
v0 0 І,
(І)
Линейное приближение: собственные векторы уравнения движения
Для детального исследования особенностей двигательной координации был проведен анализ уравнения движения в линейном приближении вблизи вертикального положения звеньев тела |^|<<1.
В этом приближении уравнение движения (1) принимает вид
Сер - Бф = Т, (2)
где, в отличие от уравнения движения (1), элементы инерционной матрицы С и гравитационной матрицы Б не зависят от суставных углов и определяются как С = С(0), Б = Б(0), а также принято допущение о пренебрежимой малости центробежных и кориолисовых сил по сравнению с инерционными и гравитационными, т.е. в уравнении (1) положено А = 0. В работе [15] показано, что ошибка линейного приближения при расчете суставных углов по заданным суставным
моментам (т.е. при решении прямой задачи динамики) приводит к ошибкам не более
4% от амплитуд углов вплоть до амплитуды наклона 55°.
Линеаризованное уравнение движения (2) позволяет провести разложение любого трехсуставного движения на три компоненты (три «собственных движения»), каждая из которых в пространстве трех суставных углов представляет собой движение вдоль одного из трех собственных векторов wг■ (, = А, Н, К) уравнения движения (2):
ф = , (3)
где ’ есть матрица 3 х 3, столбцами которой являются собственные векторы wA, wн и wк, при этом собственный вектор wг■ (, = А, Н, К) по определению удовлетворяет уравнению
С^,- = х, Dw,, (4)
где X, есть соответствующее собственное значение. Вектор £={Ъа, £к, Ън}Т назовем вектором кинематических амплитуд собственных движений. Его компоненты ^л(0, £,н(^ и находятся по известным суставным углам ф путем обратного решения уравнения (3):
Ц1) = ’-1 ф(Г). (5)
Значение нижних индексов А, н и к в применении их к собственным векторам поясняется ниже.
Векторное уравнение движения (2), записанное в терминах суставных углов ф ,
в терминах кинематических амплитуд £ распадается на три независимые скалярные уравнения
-х&+£, =л, (, = А, Н, К), (6)
каждое для отдельной кинематической амплитуды .
Назовем скалярные функции (^) в уравнении (6) «динамическими
амплитудами» собственных движений (по аналогии с кинематическими амплитудами ), каждая из которых определяет временной ход моментов сил в суставах в соответствующем собственном движении.
Вектор суставных моментов Т представляется в виде:
Т = -Б’п = и а П а 0) + и н Пн 0) + ик Пк 0), (7)
где векторы
и, = -Dw, (8)
являются столбцами матрицы -Б’.
А-синергия Н-синергия К-синергия
Рис. 2. Кинематическая и динамическая модели собственных движений (синергий). Размеры кружков, обозначающих суставы, пропорциональны абсолютным величинам моментов сил в суставах в каждой из синергий
Как видно из уравнений (3) и (7), если движение является собственным для данной биомеханической системы, т.е. происходит вдоль какого-либо собственного вектора wг■, то компоненты этого вектора определяют отношения между изменениями
суставных углов, а компоненты соответствующего вектора и, определяют отношения между изменениями суставных моментов в данном собственном движении.
Три собственных вектора w,, соответствующие им собственные значения X, и
векторы и, были вычислены согласно уравнениям (4) и (8), с использованием
антропометрических коэффициентов [14] для «стандартного испытуемого» массой 70 кг и ростом 1,7 м. Для такой «стандартной» антропометрии они составили
X А = 0,108 с2, X н = 0,020 с2, X к = 0,0021 с2,
'- 0,94 ^ Г- 0,24N
w A = - 0,09 , w H = 0,06
V - 0,32 j V 0^7,
Г668,4Л Г 32,1 ^
u A = 382,5 , u H = - 41,4
V143,1 J v- 83,4 J
W к
u
к
Г- 0,33 ^ 0,80
V- 0,52 J
Г 6,3 ^ - 96,0 7,7
(9)
J
В уравнении (9) собственные векторы W; и соответствующие векторы Ui (i = A, H, K) упорядочены по возрастанию собственных значений X;. Схематические диаграммы, представленные на рис. 2, иллюстрируют кинематические (соотношения между изменениями суставных углов) и динамические (соотношения между моментами сил в суставах) модели в каждом из трех собственных движений (собственных синергий), параметры которых приведены в (9).
В каждом из векторов wа и ua первая компонента, соответствующая голеностопному суставу, является доминирующей. Поэтому трехсуставное движение вдоль собственного вектора wа будет именоваться А-синергией (А - Ankle), и все обозначения для параметров этого движения отмечаются нижним индексом «А». В векторе wA компоненты для изменений углов во всех трех суставах имеют один и тот же знак. Таким образом, A-синергия соответствует наклону, при котором отклонения нижнего, среднего и верхнего звеньев происходят в одну сторону, т.е. в целом происходит однонаправленное вращение всего тела относительно голеностопного сустава (рис. 2 слева).
В каждом из векторов wh и uh доминирующей является третья компонента, соответствующая тазобедренному суставу. Поэтому трехсуставное движение вдоль собственного вектора wн именуется ниже H-синергией (H - Hip), и все обозначения для параметров этого движения отмечаются нижним индексом «Н». В векторе wн компоненты для углов в голеностопном и тазобедренном суставах имеют противоположный знак, в то время как компонента для коленного сустава (вторая) близка к нулю. Таким образом, Н-синергия соответствует наклону, при котором флексия в тазобедренном суставе сопровождается экстензией в голеностопном при незначительной экстензии в коленном (рис. 2 в центре). Подобная кинематическая модель близка к экспериментально наблюдаемому наклону корпуса вперед (рис. 1).
Наконец, в каждом из векторов wk и uk доминирующей является компонента для коленного сустава. Поэтому трехсуставное движение вдоль данного собственного вектора wк именуется ниже К-синергией (K - Knee) и для этого движения используется нижний индекс «К». К-синергия напоминает приседание за счет сгибания коленей при близком к поступательному вертикальному движению корпуса (рис. 2 справа).
Каждое из трех уравнений (6) представляет уравнение движения для перевернутого маятника. Три маятника (i=A, H, K) имеют разные инерционные свойства, определяемые собственными значениями Xа , Xн и Xк. Как видно из (9), отношения собственных значений равны приблизительно 50:10:1, т.е. инерция движения по А-синергии приблизительно в 5 раз выше, чем инерция движения по Н-синергии. В свою очередь, инерция движения по Н-синергии приблизительно в 10 раз выше, чем по К-синергии.
Перемещения центра тяжести и центра давлений в собственных движениях
Согласно уравнениям (3), (7), (П1-П3) координаты общего центра тяжести СО и
CG CP
центра давлений СР в переднезаднем направлении, X и X , представляют собой алгебраическую сумму соответствующих координат СО и СР в каждом из собственных движений, которые определяются через кинематические и динамические амплитуды собственных движений:
XCG = Ь, , Хср = Ь, (л,- - h ^/ш), (10)
где Ь, = иц/(Mg), а Ыц является первой компонентой вектора и,, сответствующей
голеностопному суставу. Следует отметить, что согласно уравнению (4) знак любого из собственных векторов произволен. При выбранных в (9) знаках компонент
CG
собственных векторов wг■ перемещение X, в положительном направлении (вперед)
соответствует положительному изменению кинематической амплитуды .
Для «стандартного испытуемого» [14] массой 70 кг и ростом 1,7 м
коэффициенты Ь, составляют
ЬА = 9,3-10-3 м/рад, Ьн = 4,5 -10-4 м/рад, Ьк = 8,5 -10-5 м/рад. (11)
Учитывая уравнения (6) и (10), получаем следующее уравнение, связывающее перемещения центра тяжести и центра давлений в каждом из собственных движений:
- (X, + Щш )Хр + хр = хСР. (12)
250 мс
Рис. 3 Перемещение центра тяжести вперед (жирная линия) на 1 см за 500 мс и обеспечивающие данное перемещение изменения центра давлений (тонкие линии) в
каждой (А-, Н- и К-) из синергий
Общее решение уравнения (12) имеет вид:
ХС0 (I) = С^ + С2/х + ХСр (I) - 0,5(й1(;) г-Ч^ + gi г(г) гЧ'<), (13)
где С и С2 - произвольные постоянные, X2 = X; + И/g,
I I
gn(t) = / ХСр (р) е-р/хф, g,l (I) = /ХСр (р) ер'>ф. (14)
0 0
Из всего класса решений (13) рассмотрим далее только те, которые начинаются из статического состояния и заканчиваются также статическим состоянием, как это происходит при экспериментально выполняемом наклоне корпуса с остановкой в наклонном положении.
Тогда перемещение центра давления происходит на некотором ограниченном отрезке 0 < I < Т, и перемещение центра тяжести задается уравнением
ХС0(I) = ХСр(I)-0,5gл(t)е-'/х + 0,5(^2(Т)-gi2(0)е‘/х . (15)
Рис. 3 демонстрирует соотношения между перемещениями центра тяжести СО (жирные линии) и центра давлений СР (тонкие линии) для каждого из собственных движений. При этом для иллюстрации было взято полное смещение центра тяжести величиной в 1 см, а переход из начального положения в конечное происходил за
СО
500 мс. Для движения центра тяжести был выбран «типичный» временной ход Х; (I)
с колоколообразным профилем скорости и с нулевыми значениями первой и второй производной по времени в начале и в конце движения [16]. Затем по уравнению (12) вычислялись перемещения центра давления.
В каждом из собственных движений начальное смещение центра давления СР инициирует «толчок» центра тяжести СО, в результате которого СО приобретает ускорение в противоположном направлении (рис. 3). Затем СР «возвращается» и «обгоняет» СО, что приводит к торможению и последующей остановке центра тяжести. Поскольку инерция движения по А-синергии самая высокая среди трех
СО
собственных движений, то один и тот же сдвиг центра тяжести АХ за одно и то же время х может быть осуществлен по А-синергии за счет наибольшей экскурсии центра давления (т.е. наибольшим необходимым размером опоры, 6,5 см во взятом примере). И, наоборот, благодаря малой инерции движения по К-синергии экскурсия СР при таком движении самая малая среди трех собственных движений (1 ,2 см), а траектории центра давления СР и центра тяжести СО отличаются в К-синергии незначительно.
Таким образом, чтобы обеспечить смещение СО на 1 см за 500 мс, ширина опоры при движении по А-синергии должна быть не менее 6,5 см, а при движении по К-синергии
- не менее лишь 1 ,2 см.
Казалось бы, что К-синергия является наиболее эффективной среди трех собственных движений для контроля за движением центра тяжести в переднезаднем направлении, поскольку она позволяет выполнять движение на самой узкой опоре. Однако, в соответствии с (10) и (11), для того чтобы обеспечить смещение центра тяжести заданной величины, К-синергия требует самой значительной амплитуды движения. Так, например, для смещения СО на 1 см по К-синергии «стандартному испытуемому» потребовалось бы осуществить движение с размахом кинематической амплитуды в 65°, что соответствует изменению угла в коленном суставе на 51°. В то же время аналогичное смещение центра тяжести по А-синергии достигается при кинематической амплитуде собственного движения лишь в 0,61 °, что соответствует изменению угла в голеностопном суставе на 0,58°.
Моделирование наклона без перемещения центра тяжести или центра давления
В настоящем разделе теоретически исследуются две рабочие гипотезы о том, какой параметр, СО или СР, стабилизируется биомеханической системой при выполнении наклона корпуса для того, чтобы сохранить равновесие. С этой целью моделируются две крайние «идеальные» двигательные стратегии, при одной из которых наклон корпуса происходит при отсутствии каких-либо перемещений центра тяжести СО в переднезаднем направлении, а при второй - отсутствует перемещение центра давления СР. Экспериментальные данные показали, что К-синергия при наклоне корпуса вперед практически не дает вклада в кинематику движения и в перемещения СО и СР [17]. Поэтому ниже мы везде ограничиваем анализ рассмотрением лишь двух, Н- и А- синергий.
В соответствии с уравнением (10) движение по Н-синергии приводит к
СО
перемещению центра тяжести Хн = Ън £ н. Для полной компенсации этого
перемещения собственное движение по А-синергии должно вызывать
СО СО
противоположное перемещение СО ( Х А = - Хн ). Тогда
Чл(>) = -(Ън/ЪлЪн (О- (16)
Учитывая (11), видим, что кинематическая амплитуда А-синергии составляет лишь 5% от амплитуды Н-синергии, т.е. полные угловые изменения на 95% описываются Н-синергией, в связи с чем мы будем называть ее основной, в то время как А-синергия, несмотря на ее низкую «представленность» в кинематике, полностью компенсирует перемещение СО по основной синергии, в связи с чем ее логично назвать сопутствующей позной компонентой движения.
Для численного эксперимента временной ход перемещения центра тяжести в
СО
движении по основной Н-синергии Хн (I) был взят с «типичным» колоколообразным
профилем скорости (рис. 3, жирная линия). Амплитуда перемещения центра тяжести
СО
по Н-синергии АХн и его длительность хн были взяты равными соответственно
СО
АХн = 4,5 см и х н =510 мс. Такие параметры движения по Н-синергии соответствуют средним величинам, наблюдаемым экспериментально [16].
А Н Сумма
500 мс
Рис. 4. Перемещения центра тяжести CG (жирные линии) и центра давлений CP (тонкие линии) в трех моделируемых типах наклона корпуса. Полные перемещения (правые графики) представлены в виде суперпозиции соответствующих перемещений в Н- и А-синергиях. Диаграммы сверху иллюстрируют кинематическую модель синергий. Движение моделируется с неподвижным полным (суммарным) CG (а) или с неподвижным полным CР (б, в).
Вертикальные линии на всех графиках отмечают начало и конец движения в Н-синергии
Рассматриваемое движение представлено на рис. 4а. Перемещение центра тяжести в Н- и А-синергиях показано жирными линиями. Перемещения центра давления, рассчитанные по формуле (12), показаны тонкими линиями. Вертикальными отрезками на всех графиках на рис. 4 отмечены начало ( = 0) и конец ( = Т) составляющей наклона по основной Н-синергии.
На рис. 4а (справа) показаны полные перемещения СР и СО, которые являются суммами соответствующих перемещений в Н- и А-синергиях. В то время как центр тяжести СО неподвижен, полный размах центра давлений СР составил 17,6 см. Таким образом, рассматриваемый наклон может быть выполнен только при условии, что размеры опоры в переднезаднем направлении превышают 17,6 см, несмотря на то, что центр тяжести не перемещается совсем.
Наблюдаемый экспериментально размах СР составляет 5-8 см [17], что существенно меньше полученного для движения без перемещения СО. Это позволяет заключить, что стратегия поддержания равновесия при реальном наклоне корпуса человека не направлена на полную стабилизацию перемещения центра тяжести в переднезаднем направлении.
Возникает естественный вопрос, возможно ли с биомеханической точки зрения в принципе выполнить наклон корпуса без какого-либо перемещения центра давления, сохранив при этом равновесие? Такое гипотетическое движение можно было бы совершить на бесконечно узкой планке («на лезвии ножа»), и в этом смысле его можно назвать «оптимальным наклоном» с точки зрения сохранения равновесия.
СС СР
Пусть перемещение центра тяжести и давления Хн и Хн по основной Н-синергии, обеспечивающей наклон как таковой, будут те же, как на рис. 4а (слева). Для наглядности кривые продублированы на рис. 4б (слева).
СР
Для компенсации перемещения центра давлений Хн , вызываемого Н-синергией, А-синергия должна быть скоординирована таким образом, чтобы вызывать
перемещение СР, равное по абсолютной величине, но противоположное по знаку:
СР СР
X А = - Хн (рис. 4б, тонкая линия на среднем графике). Тогда согласно уравнению (15) для перемещения центра тяжести по А-синергии справедливо уравнение
ХА0(') = -ХНР(')-0,58а1(')е~'/тЛ + 0,5(8а2(Т)-8Л2«)е'1'Л ,
где
' '
8 Л1(') = -{ХН (р)е~р/тАЛр, 8а2(') = -{ХН (р)е?/тАЛр. (17)
0 0
Перемещение CG по А-синергии показано на рис. 4б (центр) жирной линией.
Суммарное перемещение центра тяжести показано на графике рис. 4б (справа) жирной
линией. Экскурсия суммарного перемещения центра тяжести CG составила 1,55 см, в то время как суммарный центр давления СР неподвижен. Сравнительно малая экскурсия CG свидетельствует о том, что кинематическая амплитуда движения по Асинергии так же мала, как в случае без перемещения CG, а основная Н-синергия и в оптимальном наклоне описывает около 95% кинематики.
Как видно из рис. 4б, оптимальный наклон разделяется на три фазы. В начальной фазе ('<0) присутствует только движение по А-синергии. В этой фазе
перемещение CG задается уравнением Хл°(') = 0,58а2(Т)е'ХА . В средней фазе (0 < ' < Т) присутствуют составляющие движения как по А-синергии, так и по основной Н-синергии. И, наконец, в заключительной фазе (' > Т) вновь движение совершается
только по А-синергии ХС° (') = - ХН (') - 0,5 8л1(' ) е ~'' А .
В начальной фазе гипотетический испытуемый пассивно падает вперед под действием силы тяжести, как перевернутый маятник, выведенный из положения равновесия. Такое начальное падение вдоль собственного вектора wа предполагает, что движение инициируется тремя бесконечно малыми моментами сил в суставах, причем соотношения между величинами запускающих моментов определяются вектором и а , приведенном в уравнении (9). Как видно на рис. 4б (справа), на протяжении начальной фазы центр тяжести CG расположен спереди от центра давления СР. При этом он постепенно увеличивает свою скорость под действием направленной вперед силы реакции опоры. Характеристическое время пассивного падения составляет т л = 0,34 с.
В течение средней фазы смещение центра тяжеста CG вперед, вызванное высокоинерционной А-синергией, останавливается, и CG даже смещается назад по отношению к начальному положению благодаря включению основного собственного движения по менее инерционной Н-синергии. Такая реверсия движения общего центра тяжести связана с тем, что Н-синергия при наклоне корпуса вперед вызывает смещение СG назад. С динамической точки зрения это происходит под действием силы реакции опоры, которая в данном случае при инициации движения по Н-синергии направлена назад. К концу средней фазы движение по основной Н-синергии прекращается, а центр тяжести, находясь позади точки опоры, вызывает силу реакции опоры, направленную назад. При этом CG имеет скорость по А-синергии, направленную вперед.
Наконец, в заключительной фазе движение СG вперед по А-синергии тормозится направленной назад силой реакции опоры. Движение вдоль вектора w л в заключительной фазе уменьшается до нуля под пассивным действием силы тяжести, аналогично перевернутому маятнику, который возвращается в положение неустойчивого равновесия.
Возникает вопрос, возможно ли совершить наклон на бесконечно узкой опоре без потери равновесия при условии, чтобы избежать движения в первой и/или третьей
со
фазах. Иными словами, возможно ли найти такое перемещение Хн ('), чтобы
8Л1(Т) = 0, и/или 8Л2(Т) = 0 . Подстановка в формулу (17) выражения для Х^0, полученного по формуле (12), дает
Т
8а1(Т) = а\ХН0 (')е-чтлл,
0
Т
8а2(Т ) = А\Х Н0 (') еЧт ЛЛ,
0
где А = (XН + И/8)/(Ха + V8) -1 ^ 0 . С учетом предположения о колоколообразности
С0
профиля скорости собственного движения Хн оба интеграла отличны от нуля, т.е.
С0
искомого собственного движения ХА для наклона корпуса человека на бесконечно узкой опоре не существует. Следует, однако, заметить, что такое движение возможно, если наклон осуществляется немонотонно.
Рис. 4в иллюстрирует невозможность выполнения наклона со стандартным колоколообразным профилем скорости без предваряющего движения по А-синергии. Как видно из рисунка, отрицательное начальное перемещение СР в высокоинерционной А-синергии (тонкая линия в центре) слишком незначительно, чтобы вызвать достаточно большое перемещение СG вперед в этой синергии (жирная линия в центре). В результате перемещение СG назад в основной Н-синергии остается нескомпенсированным и гипотетический испытуемый начинает падать назад со все нарастающей скоростью ( рис. 4в, справа).
Обсуждение результатов
Две поведенческие задачи решаются одновременно при выполнении наклона корпуса: наклон как таковой и поддержание равновесия. Основная цель работы состояла в исследовании двигательной стратегии, которую использует человек для решения этих двух задач.
В линейном приближении, которое оказалось справедливым при наклоне корпуса вплоть до 50-60°, движение было представлено в виде линейной суперпозиции трех синергичных «собственных» движений, соответствующих трем собственным векторам и собственным значениям линеаризованного уравнения движения. В каждое собственное движение вовлечены все три суставных угла пропорционально значению соответствующей компоненты собственного вектора. По доминирующему вкладу от того или иного сустава в изменение суставных углов и моментов сил в каждом из собственных движений (рис. 2) они были названы «A»-, «H»- и «К»-синергиями (Ankle
- голеностопный сустав, Hip - тазобедренный сустав, Knee - коленный сустав). Основным преимуществом перехода от описания движения в терминах суставных углов к его описанию в терминах собственных движений является распад системы связанных динамических уравнений на три независимых уравнения для каждого собственного движения в отдельности. Каждое из трех независимых уравнений есть уравнение движения для перевернутого маятника со значительно отличающимися
инерционными характеристиками. Собственные движения по А- и К-синергиям оказались соответственно наиболее и наименее инерционными (собственное значение для А-синергии приблизительно в 50 раз больше, чем для К-синергии), а собственное движение по Н-синергии - «среднеинерционным» (собственное значение для Н-синергии приблизительно в 5 раз меньше, чем для А-синергии). Дополнительным преимуществом такого методического подхода явилась возможность исключить из рассмотрения К-синергию, поскольку при реальном выполнении наклона корпуса вперед ее вклад в кинематику и динамику оказался пренебрежимо малым [17]. Таким образом, анализ биомеханической системы с тремя степенями свободы свелся к анализу системы с двумя степенями свободы (А- и Н-синергиям).
Было смоделировано «оптимальное» движение, в котором выполняются максимально жесткие требования к возможности сохранения равновесия. Поскольку равновесие не нарушается лишь в том случае, если центр давления находится внутри границ опоры, было принято по определению считать наклон оптимальным, если равновесие не нарушается даже тогда, когда эти границы можно было бы предельно сузить. Иными словами, при оптимальном наклоне равновесие сохраняется при неподвижном центре давления. Показано, что выполнение наклона корпуса без перемещения центра давления требует очень сложной координации двух собственных движений: моменты запуска и остановки движений и величины усилий должны быть тщательно подобраны. Координация собственных движений состоит в тонко настроенной балансировке между моделью внутренних мышечных усилий и внешними гравитационными силами с целью обеспечить в результате необходимый временной ход горизонтальной компоненты силы реакции опоры, которая управляет движением СО в переднезаднем направлении. Запуск движения по Н-синергии (левые вертикальные отрезки на рис. 4б) должен быть выбран точно в момент, когда скорость движения СО вперед по А-синергии уже достаточно большая, чтобы скомпенсировать предстоящее движение СО назад по Н-синергии. Аналогично остановка движения по Н-синергии (правые вертикальные отрезки на рис. 4б) должна произойти точно в момент, когда СО уже находится сзади от точки опоры, а величина остаточной скорости движения СО вперед по А-синергии в точности такова, чтобы СО вернулся и остановился в начальном положении.
В экспериментально наблюдаемом движении [17], как и в гипотетическом оптимальном движении, движение по Н-синергии описывает более 95% угловых изменений. Движение по А-синергии также сопровождается смещением центра тяжести вперед, компенсируя его смещение назад, провоцируемое Н-синергией. При этом разгон СО вперед по А-синергии также начинается раньше движения по Н-синергии.
500 мс
Рис.5. Перемещения центра давлений в Н-синергии (тонкая линия), в А-синергии (пунктир) и его полное (суммарное) перемещение (жирная линия) при типичном экспериментальном наклоне корпуса вперед.
Различие между экспериментальным и гипотетическим оптимальным наклонами состоит в различном способе инициации упреждающего движения центра тяжести вперед по А-синергии. На рис. 5 показана типичная запись перемещений центра давления при экспериментальном наклоне корпуса вперед (описание методики эксперимента приведено в [17]. В начальной фазе экспериментально наблюдаемого движения присутствует вклад флексорного момента сил в голеностопном суставе. Такой момент сил, разгоняющий центр тяжести вперед, вызывает небольшое перемещение центра давления назад, соответствующий начальному участку кривых на рис. 5. Таким образом становится понятным функциональное значение наблюдаемой обычно ранней активации флексоров голеностопного сустава [6-8, 18-20]. Смысл этой ранней активации - упреждающее ускорение центра тяжести вперед по А-синергии с целью снижения возмущения равновесия при последующем запуске Н-синергии.
Результаты исследования позволяют выдвинуть гипотезу, что собственные движения не является лишь формальным инструментом для удобного анализа уравнения движения, а представляют для центральной нервной системы две независимые, но координируемые единицы двигательного управления, которые служат для решения двух функционально разных поведенческих задач: Н-синергия
обеспечивает выполнение основной кинематической задачи (наклон как таковой), а Асинергия есть сопутствующая позная компонента движения, обеспечивающая сохранение равновесия. В пользу данного предположения свидетельствуют следующие аргументы.
1. В каждом из собственных движений суставные углы и моменты сил в суставах изменяются синхронно во всех трех суставах. Таким образом, собственное движение, как по кинематике, так и по динамике, имеет относительно простую организацию в виде одномерной временной развертки с одной определенной инерцией. Движение тела в сагиттальной плоскости с любой другой моделью изменения суставных углов является «составным» по своим инерционным характеристикам (представляет суперпозицию разных по инерции собственных движений) и, следовательно, требует более сложного управления. Сложность управления «составным» движением связана с тем, что, как оказалось, собственные движения имеют существенно разные инерции и, таким образом, требуют управляющих сигналов с существенно различающимися временными характеристиками.
2. Экспериментально наблюдаемая кинематическая модель наклона корпуса вперед в основном обеспечивается движением по Н-синергии. Следовательно, логично рассматривать Н-синергию как основное движение, обеспечивающее наклон как таковой. А-синергия незначительно представлена в наблюдаемой кинематике. Тем не менее, ее вклад в наблюдаемые перемещения центров тяжести и давления мог быть значительным. У некоторых испытуемых она практически не наблюдалась, у других она перекомпенсировала смещение центра тяжести, провоцируемого Н-синергией [17]. Следовательно, логично рассматривать А-синергию как сопутствующую (полезную, но не обязательную) позную компоненту движения, обеспечивающую компенсацию динамических возмущений равновесия, которые вызывает основная Н-синергия. При уменьшении размеров опоры амплитуда А-синергии возрастала у тех испытуемых, у которых в норме она была мала, и снижалась в противоположном случае [17]. При этом амплитуда Н-синергии не изменялась. Таким образом, рассматриваемые синергии могут управляться независимо.
3. В пользу того, что А-синергия является позной составляющей движения, служит тот факт, что она в среднем опережала и длилась дольше, чем основная Н-синергия. Как показали результаты моделирования оптимального наклона, такая координация основной и позной компонент движения объясняется разницей их
инерционных характеристик. Опережающая активация позной компоненты является типичной при организации многосуставных движений, в которые вовлечены сегменты со значительно различающимися инерционными характеристиками [3, 21].
Работа поддерживалась грантами РФФИ (проекты 02-04-48410а и 01-04-48924а) и РГНФ (проект 00-06-00242а).
Список литературы
1. Hess W.R. Teleokinetisches und ereismatisches Kraftesystem in der Biomotorik? // Helv. Physiology. Pharmacology Acta. 1943. V. 1. P. 62-C63.
2. Jung R., Hassler R. The extrapyramidal motor system. In: Field J., Magoun H.W., Hall V.E. (Editors). Handbook of Physiology, sect. 1, Neurophysiology, 2: Chap. 35. Bethesda: American Physiological Society. 1960.
3. Беленький В.Е., Гурфинкель В.С., Пальцев Е.И. Об элементах управления целенаправленными движениями // Биофизика. 1967. Т. 12. С. 135-141.
4. Cordo P.J., Nashner L.M. Properties of postural adjustments associated with rapid arm movements // Journal of Neurophysiology. 1982. Т. 47. С. 287-302.
5. Babinski J. De l'asynergie cerebelleuse // Revew Neurology. 1899. V. 7. P. 806-816.
6. Oddsson L., Thorstensson A. Fast voluntary trunk flexion movements in standing: primary movements and associated postural adjustments // Acta Physiologica Scandinavica. 1986. V. 128. P. 341-349.
7. Crenna P., Frigo C., Massion J., Pedotti A. Forward and backward axial synergies in man // Experimental Brain Research. 1987. V. 65. P. 538-548.
8. Pedotti A., Crenna P., Deat A., Frigo C., Massion J. Postural synergies in axial movements: short and longterm adaptation // Experimental Brain Research. 1989. V. 74. P. 3-10.
9. Barin K. Evaluation of a generalized model of human postural dynamics and control in the sagittal plane // Biological cybernetics. 1989. V. 61. P. 37-50.
10. Yung J.F., Winter D.A., Wells R.P. Postural dynamics in the standing human // Biological Cybernetics. 1990. V. 62. P. 309-320.
11. Kuo A.D., Zajac F.E. A biomechanical analysis of muscle strength as a limiting factor in standing posture // Journal of Biomechanics.1993. V. 26. Suppl 1. P. 137-150.
12. Kuo A., Zajac F. Human standing posture: Multijoint movement strategies based on biomechanical constraints // Progress in Brain Research. 1993. V. 97. P. 349-358.
13. Gurfinkel V.S., ShikM.L. The control of posture and locomotion. In: Gydikov, A.A., Tankov N.T., Kosarov D.S. (Editors). Motor Control, New York: Plenum Press, 1973. P. 217-234.
14. Winter D.A. Biomechanics and motor control in human movement. 2nd Ed, New York: John Wiley and Sons, 1990.
15. Alexandrov A.V., Frolov A.A., Massion J. Biomechanical analysis of movement strategies in human forward trunk bending. I. Modeling // Biological Cybernetics. 2001. V. 84. P. 425-434.
16. Alexandrov A.V., Frolov A.A., Massion J. Axial synergies during human upper trunk bending // Experimental Brain Research. 1998. V. 118. P. 210-220.
17. Alexandrov A.V., Frolov A.A., Massion J. Biomechanical analysis of movement strategies in human forward trunk bending. II. Experimental study // Biological Cybernetics. 2001. V. 84. P. 435-443.
18. Oddsson L., Thorstensson A. Fast voluntary trunk flexion movements in standing: motor patterns // Acta Physiologica Scandinavica. 1987. V. 129. P. 93-106.
19. Crenna P., Frigo C., Massion J., Pedotti A., Deat A. Forward and backward axial movements: Two modes of central control. In: Gurfinkel V.S., Ioffe M.E., Massion J., Roll J.P. (Editors). Stance and motion, facts and concepts, New York and London: Plenum Press, 1988. P. 195-201.
20. Oddsson L., Thorstensson A. Control of voluntary trunk movements in man - Mechanisms for postural equilibrium during standing // Acta Physiologica Scandinavica. 1990. V. 110. Supp. 1595.
21. Bouisset S., Zattara M. Biomechanical study of the programming of anticipatory postural adjustments associated with voluntary movements // Journal of Biomechanics. 1987. V. 20. P. 735-742.
Приложение
Элементы матрицы C'(0), матрицы D'(0) и компоненты вектора центробежных и кориолисовых сил Л'(0,0) имеют вид:
CUU = mUcU + IU ; CUM = CMU = mUcUlM cos(0U - 0m ) ;
CUL = CLU = mUcUlL cos(0U -0L ) ; CMM = mMcM + mUlM + IM ;
CML = C LM = (mMcMlL + mU^LlM )cos(0M -0 L ) ; CLL = mLcL + (mM + mU )lL + IL ;
DUu = mU%£ sin(0U V0U ; DMm = (mMcM + mUlM )g sin(0M V0M ;
DLL = (mLcL + (mM + mU )lL )g sin(0L V0L ; DUM = DUL = DMU = DUM = DLU = DLM = 0 ;
AU7 = mUcU (lL sin(0U -0L ) 0 L + lM sin(0U -0M ) 0M );
AM = (mMcMlL + mUlLlM )sin(0M -0 L )0L - mUcUlM )sin(0U - 0M )0U ;
AjL = “(mMcMlL + mUlLlM )sin(0M -0L )0M - mUcUlL sin(0U -0L )0U ,
где индексы U, M, L для звеньев соответствуют верхнему (U - Upper), среднему (М -Middle) и нижнему (L - Lower) звеньям, а m, l, c и I обозначают массу соответствующего звена, его длину, расстояние от нижележащего сустава до его центра масс и момент инерции соответствующего звена относительно его центра масс.
При заданных суставных углах координаты общего центра тяжести тела в
переднезаднем XCG и вертикальном YCG направлениях вычисляются по формулам:
xCG =-((mLcL + mMlM + mUlU )sin(9A ) + (mMcM + mUlU )sin(9 A + Фк ) +
+ mUcU sin(9A +ФK +ФH ))/M , (п1)
xCG =-((mLcL + mMlM + mUlU ) cos(фa ) + (mMcM + mUlU )cosOA + Фк ) +
+ mUcU cos(ФA +ФK +ФH ))/M , где M = mL + mM + mU есть суммарная масса тела без массы стоп.
Горизонтальная FX и вертикальная FY компоненты силы реакции опоры определяются из (п1) по ускорению центра тяжести CG:
FX = MXCG , Fy = M(g + YCG) . (п2)
CP
Перемещение центра давления CP в переднезаднем направлении, X , определяется через момент сил в голеностопном суставе TA и сил реакции опоры
XCP = (Ta - hFx )/Fj , (п3)
где h есть высота стопы (рис. 1 ).
A. Alexandrov, A. Frolov, J. Massion (Moscow, Russia, Aix-en-Provence, France)
STRATEGY OF EQUILIBRIUM MAINTENANCE DURING HUMAN FORWARD UPPER TRUNK BENDING ON THE NARROW SUPPORT
Two behavioral goals are achieved simultaneously during upper trunk bending, the bending per se and equilibrium maintenance. A biomechanical model is used for the analysis of how these goals are achieved. The investigation is based on the movement decomposition into the movements along the eigenvectors of the motion equation. The coordination between posture and movement that provides the trunk bending on the infinitely narrow support (an «optimal bending» without any displacement of the center of pressure) is found theoretically. The optimal bending is compared with the experimental one. A hypothesis that the movements along the eigenvectors represent the prime and postural movement components under independent but coordinated controls is suggested. The prime component is responsible for the bending per se, whereas the postural one compensates in the anticipatory manner an equilibrium disturbance caused by the prime component. The results of simulations clarify from biomechanical point of view the pattern of electromyographic activity usually observed in human during forward upper trunk bending.
Key words: coordination, equilibrium control, computer simulation, synergy, motor control, posture and movement.
Получено 25 декабря 2002 года