Full Text

ВВЕДЕНИЕ В последние время для оценки функциональных состояний, стабилизации вертикальной позы и диагностики неврологических нарушений человека в неврологии, ортопедии и спортивной медицине часто используют стабилоплатформы. Cтабилометрия - метод исследования функций организма человека, связанных с управлением человеком малым изменением вертикальной позы. Он основывается на оценке координат центра давления (ЦД) - точки приложения равнодействующей сил нормального давления, с которыми человек действует на опору. Одновременно показания стабилоанализатора содержат и величину равнодействующей сил нормального давления - вес. Современная стабилометрия рассматривается как простое нейрофункциональное обследование. Однако ее широкое применение затрудняется сложностью интерпретации регистрируемых сигналов в соответствии с общепринятыми в физиологии и практической медицине критериями. Существующая постурографическая литература оперирует понятиями постуральных стратегий, синергий, спектров мощности колебаний, принятыми в биомеханике и достаточно далекими от круга понятий, которые используют практикующие врачи и физиологи. В связи с этим естественной является попытка ряда авторов выделить стабилометрические показатели, наиболее информативные для медиков и биологов. Важным показателем уровня физического состояния и здоровья являются энергетические возможности человека. В связи с этим представляется интересной идея использовать энергетические оценки, высказанная в [14-16]. К сожалению, реализация этой идеи связана с существенными трудностями. Предложенный в [10] параметр, названный «индексом энергозатрат», даже в скорректированном варианте [1, 11] оценивает изменение части кинетической энергии механической системы, описывающей движение человека, в предположении, что движение медленное и ЦД совпадает с проекцией общего центра масс на опорную плоскость [5]. Таким образом, исследователь имеет дело с еще одним параметром, физическое обоснование которого имеет заметные недостатки. Отметим, что задача энергетических оценок в стабилометрии является достаточно сложной. Прежде всего следует отметить, что значимой корреляции стабилометрических показателей с оценками энергозатрат организма, получаемыми на основании метаболических оценок, полученных с помощью газоанализатора масочного типа, не обнаружено [4]. Этот результат говорит о том, что доля энергии, затрачиваемая на перемещение тела при удержании вертикальной позы, невелика по сравнению с энергозатратами, направленными на поддержание жизнедеятельности организма. Стабилометрические данные, в свою очередь, предоставляют информацию о суммарных значениях некоторых суставных моментов. Следовательно, показания стабилоанализатора позволяют оценивать механическую работу этих «суммарных суставных моментов». Разделить усилия мышц антагонистов стабилометрическая информация не позволяет в силу своего интегративного характера. Подход к подобному анализу для движений во фронтальном направлении рассмотрен ранее авторами в [9]. Цель работы - подробно рассмотреть проблемы вычисления оценок энергетических показателей по данным стабилометрии и проанализировать, опираясь на биомеханические модели движения человека, связь с этими оценками принятых стабиломерических показателей. МОДЕЛЬ ДВИЖЕНИЯ ЧЕЛОВЕКА Используем аналогичный подход к вычислению механической работы, совершаемой мышцами ног, в виде интеграла от модуля мощности, развиваемой мышцами [3]. Обсудим применение такой оценки для движения человека. Будем использовать простейшие модели перевернутого маятника, подробно описанные в [8], так как стабилографическая информация не позволяет однозначно восстанавливать состояния более сложных многозвенных систем. Для определения работы мышц нам потребуются оценки угловых скоростей изменения суставных углов и моментов в суставах. Для определения этих величин необходимо использовать математическую модель. Рассмотрим малые колебания человека в окрестности ортоградной позы, когда его тело расположено вблизи вертикали, а ноги находятся примерно на ширине плеч. В таких предположениях описание движения человека разделяется на сагиттальную и фронтальную модели. Для описания колебаний человека в сагиттальной плоскости примем традиционную модель перевернутого маятника, изображенную на рис. 1. Предполагаем, что тело человека в ходе теста допустимо моделировать недеформируемым стержнем массы m, закрепленным шарнирно в точке O, соответствующим проекции на сагиттальную плоскость оси, соединяющей голеностопные суставы человека. Центр масс стержня расположен в точке C, удаленной от точки O на расстояние Момент инерции стержня относительно фронтальной оси, проходящей через точку О, равен J. Отклонение стержня от вертикали опишем углом θ. Через обозначим момент, создаваемый мышцами в голеностопном суставе, приложенный как показано на рис. 1. Будем считать, что стопа неподвижна относительно платформы. Уравнения движения для малых значений угла θ и скорости его изменения запишем в соответствии с [8] в виде (1) Здесь g - ускорение свободного падения. Соотношения для момента получим на основании условий равновесия системы «платформа-стопа», в которых в соответствии с [8] выразим проекцию реакции опоры F на сагиттальную ось из теоремы о движении центра масс: (2) Здесь N - нормальная реакция опоры; yn - саггитальная координата ЦД, отсчитываемая от оси голеностопного сустава, а h - расстояние от точки О до плоскости, образованной чувствительными элементами сенсоров стабилоанализатора. Используем соотношения (1) и (2) для записи уравнения: (3) Здесь принято обозначение для приведенного момента инерции Рис. 1. Модель движений человека в сагиттальной плоскости. Показаны внешние силы, действующие на систему «человек - платформа» При записи соотношений для движений во фронтальной плоскости используем упрощенную модель трехзвенника с двумя параллельными звеньями, которые опираются на платформу [1, 7, 8]. Будем считать, что опора ноги на плоскость точечная и центры тазобедренных суставов и точек опоры ног образуют параллелограмм (рис. 2). За начало отсчета примем положение, которое назовем исходной позой, для которой будем далее считать, что ступни находятся непосредственно под тазобедренными суставами. Такое положение приближенно соответствует позе «ноги на ширине плеч». Обозначим через l - «длину ноги» - массу туловища с головой и руками; - массу одной ноги, т.е. Звенья и моделируют прямые ноги, третье звено - туловище с руками и головой, моделируется единым твердым телом. В соответствии с [8, 9] уравнение движения человека во фронтальной плоскости запишем в виде (4) где - угол отклонения ног от вертикали; - момент инерции ноги относительно опорной точки; - cуммарный момент мышечных усилий, в основном прилагаемых в тазобедренных суставах [7]. В случае если фронтальная координата ЦД отсчитывается от середины отрезка уравнение (4) в соответствии с [9] с учетом соотношений для равновесия платформы и основными теоремами динамики можно преобразовать к виду (5) Рис. 2. Модель движения человека во фронтальной плоскости Здесь приведенный фронтальный момент инерции a d - высота над тазобедренными суставами центра масс системы «туловище с головой и руками». Таким образом, и для сагиттали, и для фронтали уравнения математической модели представимы в одинаковом виде: (6) (7) Эти уравнения описывают фронтальные или сагиттальные движения в зависимости от значений, которые принимают параметры и переменные в соответствии с табл.1. В дальнейшем будем также пренебрегать малыми изменениями нормальной реакции N и использовать соотношение N = mg. Таблица 1 Обозначения переменных моделей Новое обозначение Фронталь Сагитталь z x y ПРОБЛЕМЫ ВЫЧИСЛЕНИЯ РАБОТЫ СУММАРНЫХ МОМЕНТОВ ОТДЕЛЬНЫХ ГРУПП МЫШЦ Механическую работу отдельных групп мышц ног человека при движении в окрестности исходной позы будем оценивать, используя интеграл от модуля мощности [3] на интервале времени [0, T], где ω - угловая скорость изменения угла наклона. Соотношения (3) и (5) используются для оценки углов наклона. Проблемы, возникающие при подобных вычислениях, будут обсуждаться в следующем разделе. Здесь же обратим внимание на то, что сопоставление начала отсчета стабилометрической платформы с проекцией середины отрезка, соединяющего центры голеностопных суставов, с миллиметровыми точностями практически невозможно. Как правило, за начало отсчета выбирают некоторое среднее положение ЦД стоящего человека. Это положение соответствует наклону туловища вперед на 1-2º [6]. Для фронтального направления проведенная центровка также может не соответствовать началу отсчета координат в рассматриваемой математической модели. В таком случае координаты zn, использованные в математической модели (6), (7), отличаются от показаний стабилоанализатора z на постоянную величину Рассмотрим подробно отмеченную ситуацию. Значения углов наклона представим в виде сумм: Постоянные составляющие значений углов связаны с координатами ЦД в соответствии c (7) формулой Аналогично значения суммарных суставных моментов также можно представить в виде сумм: Постоянные составляющие значений моментов связаны с координатами ЦД в соответствии c (6) формулой Соответствующая составляющая работы представима в виде Так как величина неизвестна, по показаниям стабилоанализатора можно оценить только второе слагаемое из правой части неравенства (8) С некоторой степенью условности можно говорить, что такая величина характеризует механическую работу, совершаемую переменной составляющей суммарного момента. Эта величина является оценкой снизу механической работы, и её можно назвать показателем механической работы для заданного направления. Общим показателем механической работы в этом случае следует считать сумму этих величин, рассчитанных для фронтали и сагиттали. АЛГОРИТМ ВЫЧИСЛЕНИЯ ПОКАЗАТЕЛЯ МЕХАНИЧЕСКОЙ РАБОТЫ Определение угловой скорости и момента M по стабилометрическим данным аналогично задаче определения проекции центра масс. Эта задача является математически некорректной, так как требует численного интегрирования уравнения (7), решения которого неустойчивы. Различные методы регуляризации этой задачи рассмотрены в [2, 17]. Модифицируем для решения поставленной задачи алгоритм, использующий метод фильтрации сигнала с помощью прямого и обратного преобразования Фурье. Этот алгоритм в отличие от использованного в [9] не требует интервалов спокойного стояния в начале и конце записи и поэтому может применяться при обработках любых обследований, в ходе которых человек не совершал значимых изменений позы. Модифицированный алгоритм предполагает, что выходные сигналы для и M не содержат экспоненциальных составляющих, порожденных решением однородного уравнения для (7). Тогда для показаний стабилоанализатора по каждой из осей алгоритм сводится к следующей последовательности операций. 1. Сигнал стабилоанализатора на интервале времени [0, T] представляем в виде где 2. Используя преобразование Фурье, строим Δ(ω) образ функции δ(t). 3. Комплекснозначные Фурье-образы W(ω) и V(ω) для составляющих угловой скорости и момента получаем в результате умножения Δ(ω) на значения передаточных функций, соответствующих уравнениям (6) и (7). (9) 4. Используя обратное преобразование Фурье, строим их праобразы w(t) и v(t). 5. Оценки угловой скорости и момента получим по формулам 6. Подставив эти значения в выражение (8), вычисляем оценку показателя механической работы Такой алгоритм очевидно требует достаточно сложных вычислений. В связи с этим интересно рассмотреть соотношения для вычисления показателя механической работы в предположениях, близких к использованным в соотношениях, обосновывающих «индекс энергозатрат» [5]. Приближенная оценка показателя механической работы для медленных движений Рассмотрим соотношения (9). Сделаем два упрощающих предположения. Это соотношение для сагиттального направления верно, если мы будем пренебрегать высотой голеностопного сустава над плоскостью чувствительных элементов силовой платформы. В этом случае и, как отмечалось ранее в [8] момент определяется изменением стабилограммы: Будем рассматривать диапазон низких частот Тогда и, следовательно, Это утверждение эквивалентно предположению о том, что скорости ЦД и проекции центра масс совпадают. Такое предположение по факту использовано в обосновании введенного в [5] «индекса энергозатрат». Тогда верно Для интервала знакопостоянства произведения парциальная составляющая На этом интервале знак не меняется, а значит, возможна грубая оценка: Здесь суммирование ведется по всем j таким, что съем измерений производится в моменты . Тогда грубая оценка показателя механической работы может быть представлена в виде где N - количество дискретных измерений на интервале [0, T]. Таким образом, для очень медленных движений для частот ниже частоты основных колебаний по [6, 12], для которых процент тактовых интервалов с переменами знаков мал, в качестве «энергетического» можно попробовать использовать показатель (10) Результаты исследования В исследовании приняли участие 10 здоровых человек в возрасте 19-22 лет с нормальным зрением, давших информированное согласие на участие в обследовании. Стабилограммы регистрировали с помощью компьютерного стабилоанализатора с биологической обратной связью «Стабилан-01-2». Обследуемый стоял на стабилографической платформе перед большим панорамным экраном с открытыми глазами. Исследование проводилось в несколько этапов, длительность которых составляла приблизительно 80 с. На каждом этапе вначале регистрировали стабилограмму обследуемого в условиях ортоградной позы, длительностью 20 с, после чего ему демонстрировали одно из изображений, вызывающих оптическую иллюзию движения: «Крутящийся силуэт», «Пульсирующий рисунок», «Train Moves Both Ways» и т.п., и на ее фоне также писали стабилограмму. В заключение у обследуемого повторно регистрировали стабилограмму в условиях ортоградной позы, длительностью 20 с. После каждого этапа обследуемому давался пятиминутный перерыв в условиях сидя. Всего в эксперименте было использовано 9 различных зрительных иллюзий движения. Для анализа корреляционных связей различных показателей были использованы результаты 20 секундных интервалов 238 стабилометрических проб. Из них в 122 интервалах обследуемые стояли с открытыми глазами в ортоградной позе и в 116 интервалах им предъявлялись различные зрительные иллюзии движения, вызывающие изменение постуральных реакций и соответственно стабилометрических показателей. В ходе анализа результатов обследований вычислялись следующие показатели: 1) показатель механической работы для фронтали и сагиттали и их сумма 2) приближенные оценки для показателя механической работы и из (10); 3) индекс энергозатрат в соответствии с [11] размахи колебаний по фронтали X и сагиттали Y, и средняя скорость центра давления V в соответствии с [13]. Для углубления анализа вычислялись и анализировались проекции средней скорости на фронтальную и сагиттальную оси, а также фронтальный и сагиттальный аналоги «индекса энергозатрат», вычисляемые по формулам Для результатов выбранных обследований характерен значительный разброс приведенных показателей. Для каждого из обследуемых среднее квадратическое отклонение скоростей и размахов составляет величины порядка 20-25 % от значений показателей, а для «энергетических» показателей более 30 %. Таким образом, изменение выбранных показателей являлось значимым. В табл. 2-4 приведены коэффициенты корреляции, вычисленные для указанных показателей по всему ансамблю результатов обследований. В приведенных таблицах опущены корреляционные коэффициенты, связывающие фронтальные и сагиттальные показатели, так как они изменяются в диапазоне 0,2-0,7 и не представляют интереса. Так же ведет себя и большая часть коэффициентов из приведенных таблиц. Это объясняется тем, что стабилометрические параметры отображают многие факторы системы удержания вертикальной позы, отражают особенности сложных движений и не являются независимыми. Однако на общем фоне следует выделить несколько пар величин, для которых коэффициент корреляции близок к единице (соответствующие ячейки таблиц выделены серым цветом). Таблица 2 Корреляционные коэффициенты фронтальных показателей Корреляция 0,39 0,38 0,86 0,51 1 0,43 0.40 0,47 - 1 0.57 0,93 - - 1 0,66 Таблица 3 Корреляционные коэффициенты сагиттальных показателей Корреляция 0,55 0,61 0,90 0,62 1 0,66 0,55 0,65 1 0,65 0,95 1 0,67 Таблица 4 Корреляционные коэффициенты суммарных показателей Корреляция 0,57 0,57 0,63 1 0,66 0,66 1 0,95 Прежде всего эти оценки показателя механической работы для движений во фронтальном и сагиттальном направлениях скоррелированы с соответствующими размахами X и Y. Вспомним, что работа связана с изменением полной механической энергии системы, которая является суммой потенциальной и кинетической энергий. При максимальном отклонении, в случае когда скорость равна нулю, полная механическая энергия равна потенциальной энергии системы. Таким образом, механическая работа на участке между двумя экстремальными отклонениями будет определяться разностью потенциальных энергий в этих положениях. Изменение потенциальной энергии в свою очередь определяется максимальной разностью значений угла наклона, который связан с размахом колебаний ЦД. Эти соображения поясняют высокую корреляцию показателей механической работы с размахами. Этот факт подтверждает применимость введенной нами оценки, учитывающей значительную часть механической работы мышц. Корреляция индекса энергозатрат E c оценками показателя механической работы как для фронтального и сагиттального направлений, так и в общем случае не превышает общего уровня значений коэффициентов корреляции для пар разнородных стабилометрических показателей. Высокие коэффициенты корреляции отмечаются для пар E и средней скорости V центра давления, а также для его фронтального и сагиттального аналогов с соответствующими проекциями средней скорости. Это значит, что изменение энергии, соответствующее индексу E, невелико и не характеризует в должной степени изменение механической энергии. Отметим, что показатель P также не коррелирован с и не характеризует в должной степени совершаемую механическую работу. Это очевидно связано с существенным вкладом быстрых движений. Тем не менее P представляет собой относительно независимый интегральный показатель, и его рассмотрение может быть интересно. Наиболее перспективным представляется доработка алгоритма расчета подобного показателя в целях его использования для оценки энергетического вклада медленных движений при рассмотрении разнотемповых составляющих движения в соответствии с [18]. Заключение Стабилоанализатор является устройством, совершающим измерения механических величин, причем стабилометрические измерения не позволяют разделять информацию об усилиях мышц антагонистов. В соответствии с этим стабилоанализатор позволяет оценивать только часть энергозатрат: механическую работу, совершаемую суммарным моментом, создаваемым отдельными группами мышц человека. Для сагиттального направления речь идет о мышцах, создающих момент в голеностопном суставе. Для фронтального - в основном в тазобедренном. Задача оценки механической работы суставных моментов является достаточно сложной и требует использования механической модели и применения специальных математических методов. Для её решения необходимо знание масс инерционных характеристик тела человека. В статье предложен алгоритм расчета оценки значимой части механической работы, основанный на простых моделях. Этот алгоритм пригоден для анализа движения обследуемых не совершающих заметных изменений углов в коленном и тазобедренном суставах, удерживающих спину прямо. При анализе движений значительных групп пациентов с двигательными нарушениями алгоритм может давать существенные погрешности. Математический анализ задачи позволяет выписать приближенную оценку P показателя механической работы для медленных движений. Этот показатель оказывается достаточно простым для вычисления и представляет собой относительно независимый интегральный показатель. Его рассмотрение может представлять определенный интерес для задач анализа энергетического вклада медленных составляющих движений. Изменение части энергии, соответствующее показателю, называемому индекс энергозатрат, невелико и не характеризует в должной степени изменение механической энергии. Таким образом, несмотря на то, что указанный индекс энергозатрат представляет собой интересный интегральный квадратичный скоростной показатель, его наименование не соответствует механическому содержанию вычисляемой величины и вводит пользователя в заблуждение.

About the authors

P. A Kruchinin

N. V Kholmogorova

References

Statistics

Views

Abstract - 70

PDF (Russian) - 38

Refbacks

• There are currently no refbacks.