УДК 531, 51-7 Вестник СПбГУ. Прикладная математика. Информатика... 2020. Т. 16. Вып. 2 МБС 92С10, 70Е55
Модельное изучение влияния многосуставных мышц на частотные характеристики тела человека
В. П. Трегубое, Н. К. Егорова
Санкт-Петербургский государственный университет, Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Трегубов В. П., Егорова Н. К. Модельное изучение влияния многосуставных мышц на частотные характеристики тела человека // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2020. Т. 16. Вып. 2. С. 150-164. https://doi.org/10.21638/11701/spbu10.2020.207
Проведен анализ работ, посвященных экспериментальным исследованиям частотных характеристик тела сидящего человека, подверженного вибрации, и построению его механических моделей. Отмечается, что в большинстве случаев измерения вибрации производятся на сиденье и голове человека, что позволяет построить передаточную функцию от сиденья к голове. Но предполагаемая внутренняя структура механической модели человека остается не подтвержденной. В этих условиях была подчеркнута необходимость проведения измерений вибрации на всех моделируемых частях тела человека. Кроме того, с анатомической точки зрения, основной вклад в механические свойства тела сидящего человека вносят многосуставные мышцы позвоночника. Однако данное важное обстоятельство до настоящего времени не учитывалось. В связи с этим была поставлена задача исследования влияния многосуставных мышц на частотные свойства тела сидящего человека. Был построен ряд механических моделей, в которых многосуставные мышцы моделировались многозвенными вязкоупругими соединениями. Так, на простейшей модели с двумя степенями свободы было показано, как наложение двухзвен-ных соединений в дополнение к однозвенным приводит к появлению антирезонансной частоты на верхней массе, что невозможно при их отсутствии. Предлагается механическая модель с произвольным числом степеней свободы с многозвенными соединениями, для которой получены формулы передаточной функции и амлитудно-частотной характеристики (АЧХ) для верхней массы. Кроме того, рассматривается механическая модель с произвольным набором многозвенных соединений. В качестве примера приводятся результаты численного расчета АЧХ для механической модели с восемью степенями свободы при наличии и отсутствии многозвенных соединений.
Ключевые слова: механическая модель, тело человека, вибрация, передаточная функция, входной механический импеданс, амплитудно-частотная характеристика, многосуставные мышцы.
Введение. Механические модели тела человека разрабатываются на основе экспериментальных исследований частотных характеристик тела человека. Измерения, как правило, производятся на сиденье и голове сидящего человека. Их результаты можно рассматривать как входной и выходной сигналы «черного ящика» (рис. 1).
По ним можно построить входной механический импеданс (ВМИ) [2], амплитудно-частотную характеристику (ЛЧХ) [3] или кажущуюся массу как функцию частоты [4]. При этом была установлена схожесть полученных характеристик с характеристиками механических систем, состоящих из абсолютно твердых тел, соединенных между собой последовательно пружинами и демпферами, что послужило основанием для со© Санкт-Петербургский государственный университет, 2020
Рис. 1. Последовательность решения задачи моделирования тела человека, подверженного вибрации [1]
здания механических моделей с сосредоточенными параметрами цепной структуры. В частности, следует отметить механические модели тела человека для трех сидячих поз, разработанные в Институте машиноведения АН СССР [5] (рис. 2).
Рис. 2. Механические модели цепной структуры тела сидящего человека, подверженного вибрации [5]
Построение таких моделей основывается на предположении о внутренней структуре черного ящика, которая не может быть определена единственным образом. Это достаточно наглядно продемонстрировано на моделях, представленных в [6-10] (см., например, рис. 3 и 4). В данных моделях абсолютно твердые тела с массами ш\,.Ш4 никак не связаны с частями тела человека. По существу это не модели человека, а модели АЧХ, для которой построена механическая модель. Модель должна содержать моделируемые части тела человека, для каждой из которых экспериментально получены передаточные функции (ПФ). Для однозначного определения параметров модели понадобятся и средние значения масс моделируемых частей тела человека [11].
Механические модели человека, в которых авторы старались приблизить их структуру к структуре тела человека, как правило, состояли из абсолютно твердых
Рис. 3. Механические модели с четырьмя степенями свободы, описанные в работах [6] (а, б), [7, 8] (в), [9, 10] (г)
Рис. 4. Механическая модель, представленная в [6]
тел, соединенных шарнирами, которые снабжены спиральными пружинами и демпферами. Вместе с тем такие модели помогали понять, как вращения частей тела человека формируют вертикальные перемещения головы относительно сиденья. В рабо-
те [12] это было продемонстрировано на моделях с пятью степенями свободы. Однако параметры данных моделей (коэффициенты упругости и демпфирования) определялись так же, как и в моделях в [6-10], по единственной ПФ от сиденья к голове.
В работе [13] голова, верхняя часть туловища и крестец моделируются абсолютно твердыми телами, соединенными шарнирно. Кроме того, поясничный отдел позвоночника представлен пятью абсолютно твердыми телами, также соединенными плоскими шарнирами (рис. 5).
В целом механическая модель обладает восемью вращательными степенями свободы, а также двумя степенями свободы прямолинейного перемещения основания модели по горизонтали и вертикали. Все шарнирные соединения снабжены спиральными пружинами и демпферами, а крестец соединен с сиденьем вертикальной и горизонтальной продольной пружинами. Для определения параметров модели были экспериментально найдены ПФ для вертикальной вибрации от сиденья ко всем моделируемым частям тела сидящего человека. С этой точки зрения, были удовлетворены все требования к моделям тела человека, сформулированные выше. В то же время в этой модели, так же как и во всех вышеприведенных, не моделируются многосуставные мышцы позвоночника, которые в значительной степени обусловливают механические свойства позвоночника. Кроме того, датчики крепились на коже, а не к остистым отросткам позвонков, как сделано в работе [14], а размеры, массы и моменты инерции взяты из [15].
В частности, хорошо известно, что в отсутствие мышц позвоночный столб не в состоянии сохранять свою конфигурацию. Однако влияние многосуставных мышц на механические характеристики тела человека в условиях вибрации не изучалось, хотя в случае ударных воздействий такие попытки предпринимались [16].
Математическое моделирование работы скелетных мышц. Математическое моделирование работы скелетных мышц было основано на механической модели мышцы Хилла, состоящей из сократительной компоненты и двух упругих компонент: параллельной и последовательной по отношению к сократительной [17]. Им же было получено уравнение, связывающее силу и скорость изотонического сокращения мышцы. Экспериментальные исследования микроструктуры мышцы позволили Хаксли сформулировать теорию сокращения мышц на микроскопическом уровне с использованием теории скользящих нитей [18]. На ее основе и гипотезе о свободном, тянущем и тормозящем состояниях актин-миозиновых мостиков Дещеревским были выведены уравнения, описывающие укорочение сократительной компоненты скелетной мышцы [19]. Они послужили основой для развития кинетической теории мышечного сокращения [20], включая переходные режимы сокращения [21] и циклические движения мышцами-антагонистами [22].
Специфика работы скелетных мышц в условиях вибрации заключается в следующем. Во-первых, при воздействии с частотой до 1 Гц нервная система реагирует на каждое колебание, но при более высоких частотах напряжение мышц становится
Рис. 5. Модель сидящего человека [13]
практически постоянным. Было установлено, что при длительном действии вибраций на тело сидящего человека ускорения на голове уменьшаются, в то время как смещения увеличиваются [23]. Во-вторых, поскольку уровень воздействия вибрации на организм человека определяется ускорением [24], то допустимые значения амплитуды колебаний при экспериментальных исследованиях на испытуемых составляют в области первого резонансного пика (от 4 до 6 Гц) несколько миллиметров с последующим снижением в каждой октавной полосе частот. С кинетической точки зрения столь малые амплитуды вибрации не приводят к разрыву актин-миозиновых мостиков. В связи с этим силу, развиваемую мышцей, можно считать пассивной и вязко-упругой, а многосуставную мышцу моделировать многозвенным вязкоупругим соединением.
Система с двумя степенями свободы и двухзвенными соединениями.
Для начала ограничимся рассмотрением простой двухмассовой модели, снабженной двухзвенными соединениями (рис. 6). Уравнения движения механической системы,
Рис. 6. Механическая система цепной структуры с двумя степенями свободы с двухзвенными соединениями
представленной на рис. 6, имеют вид
шх¿1 + Ъх(¿1 - х2) + а (хх - Х2) + Ьо(¿1 - у) + со(хх - у) = 0, (1)
Ш2¿2 - Ъх (Хх - Х2) - Сх (¿1 - ¿2) + Ъ2 (¿2 - У)+ С2 (¿2 - У) = 0,
где — перемещение тела массы шг (г = 1, 2) относительно неподвижного основания, отсчитываемое от положения его статического равновесия; у — перемещение массы то относительно неподвижного основания; с^ и bj = 0,2) — жесткости пружин и коэффициенты вязкого трения демпферов соответственно.
Для построения ПФ рассматриваемой механической системы применим к системе уравнений (1) преобразование Лапласа (с нулевыми начальными данными)
сю
Р(р) = / /М,
о
здесь р — комплексная переменная.
Заменяя функции (£) и у(£) по правилу преобразования Лапласа их изображениями Хг (р) и У (р), получаем систему линейных неоднородных алгебраических уравнений относительно Хг (р) и У (р)
(mip2 + (bi + bo)p + (ci + co))X1 - (bip + c1)X2 = (bop + co)Y, - (bip + ci)Xi + (m2p2 + (bi + b2)p + (ci + c2))X2 = (b2P + c2)Y.
Как известно, ПФ, обозначаемая через Hi(p) и представляющая собой отношение изображений входной и выходной координат объекта, вычисляется по формуле
Xi(p)
Hi(p)
1, 2.
Y(p) '
Поделим каждое из уравнений системы (2) на Y (p):
(mip2 + (bi + bo)p + (ci + co))Hi - (bip + ci)H2 = (bop + co), - (bip + ci)Hi + (m2p2 + (bi + b2)p + (ci + c2))H2 = (b2p + c2). Разрешая систему (3) относительно Hi(p), получим равенства
Hi(p) = H2(p) =
3 2
а3рл + + <УлР + «о S^p4 + S3p3 + Ô2p2 + Sip + •
h p3 + & p2 + Pip + A) S4p4 + S3p3 + S2p2 + Sip + •
(3)
(4)
(5)
где
50 = ao = во = ci c2 + co(ci + c2 ),
51 = ai = fîi = ci b2 + c2bi + bo(ci + c2 ) + co(bi + b2),
52 = (mi + m2)ci + c2mi + bib2 + com2 + bo(bi + b2),
53 = bi(mi + m2) + b2mi + bom2,
54 = mim2,
a2 = bib2 + bo(bi + b2) + com2, a3 = bom2,
(32 = Ь1Ь2 + е2т\ + Ъо(Ь1 + Ъ2), вз = Ь2т1.
АЧХ, соответствующая ПФ И\(р) (см. (4)), как известно, получается в результате подстановки в нее вместо комплексной переменной р мнимой переменной ш:
(«0 — а2Ш2) + — аз и3)
Hi(iu) =
(So - S2и"2 + S4U4) + i(Siu - S3U3)
с последующим определением ее модуля, который равен
\Hi(u)\ =
(ao - a2U2)2 + (aiи - a3U3)2
(S0 - S2LV2 + S^LV4)2 + (SiLV - S3LV3)2 '
(6)
Аналогично рассчитывается АЧХ, соответствующая ПФ H2(p) (см. (5)):
\H2(u)
(eo - в2и2)2 + (eiU - в3и3)2
V (So - S2U2 + S4U4)2 + (Siu - S3и3)2 '
Значения коэффициентов аг и 5г в (6) определяются таким образом, чтобы максимально приблизить теоретическую АЧХ к экспериментальной. Наличие этих коэффициентов всегда позволяет построить и ПФ в изображениях.
Следует заметить, что числитель АЧХ \Нх(и)| при bj = 0 равен
ао — а2 и
2
и обращается в 0 при
и =
/С1С2 + С0(С1 + с2)
с0га2
Это означает, что при наложении двухзвенных соединений на верхней массе появляется антирезонансная частота, что в отсутствие данных соединений невозможно.
Система с произвольным числом степеней свободы. Рассмотрим механическую систему аналогичной структуры, имеющую п степеней свободы (рис. 7).
Рис. 7. Многомассовая механическая система с многозвенными соединениями
Уравнения движения такой системы в абсолютной системе координат будут
иметь вид
шг XI + Ь1 (¿1 — х 2) + сх (хх — Х2) + Ьо (хх — у) + со (хх — у) = 0,
тг¿г + Ьг (±г — хг+х) + а (¿г — хг+1) — Ь—1 (х—1 — Xг) — с—х (х—х — ¿г) = 0, (7)
Шпхп + Ьп (хп — У) + Сп (хп — У) — Ьп-1 (хп-1 — хп) — СП-1 (хта-1 — хп) = 0.
Преобразование по Лапласу п линейных дифференциальных уравнений относительно функций хг (Ь) и у(Ь) (7) приводит к системе п линейных неоднородных алгебраических уравнений относительно изображений Хг (р) и У (р). Разделив каждое из уравнений полученной системы на У (р), имеем систему уравнений
( aii ai2 0 ... 0 (Hi(p)\ fDA
a2i a22 a-23 ... 0 H2(p) D2
0 аз2 азз ... 0 Нз(р) = D3
0 0 0 .. 0 . . . ann/ \Hn(P)) \DJ
A • H = D,
где матрица A является симметричной матрицей Якоби, в которой элементы, стоящие на главной диагонали, можно записать как
_ ( rriip2 + (6i + b0)p + (ci + c0) при г = 1, ["Чр2 + (6/_i + bi)p + (cj_i + с») при i = 2,n,
а элементы, стоящие слева и справа от нее, — как
= аФ = ~(hp + Ci), i = 1, n, j = г + 1,
столбец H представляет собой столбец ПФ, а элементы столбца D имеют следующий вид:
{(bop + со) при к = 1, О при i = 2, п — 1,
(bnp + Cn) при k = n.
Разрешая систему (8) относительно Hi(p), получим уравнения
Нг(Р) =
в которых Д — определитель матрицы A, а Ai — определитель, построенный из матрицы A заменой j-го столбца на столбец D.
АЧХ, соответствующая ПФ Н\(р) при bj = 0 (j = 0,п), будет равна
HiH\ =
' (а0 - а2и;2 + <У-4^4 ~ • • • + а2п-2^2" 2)2
(¿0 — S<2L02 + S4LV4 — . . . + 5<2nL02n)2 '
(9)
В соответствии с формулой (9) число антирезонансных частот на верхней массе системы, представленной на рис. 7, равно п — 1.
Следует отметить, что многозвенные соединения могут связывать различные массы системы в зависимости от того, какие многосуставные мышцы включаются в разрабатываемую модель. Рассмотрим систему, приведенную на рис. 8.
Уравнение движения массы в общем случае имеет вид
^miXi + ^b—! + bi + ^bi^ xH + ^Ci-1 + Ci + ^C^ - (bi+iii+1 + Ci+ixi+1) — (bikxk + CikхЦ.) = 0.
(10)
В (10) — суммирование по к £ [1 < i — 2^ [i + 2,n+1](n +1 = 0), а xxn+1 = y. k
или
x
Рис. 8. Различные способы наложения многозвенных соединений на механическую систему с п степенями свободы
Для жесткостей пружин и коэффициентов вязкого трения у соединений соседних масс введено обозначение с одним индексом, а для многозвенных соединений — с двумя. В том случае, когда какие-то массы тг и тк не связаны многозвенными соединениями, соответствующие коэффициенты сгк и Ьгк будут равны нулю.
Переходя от системы дифференциальных уравнений относительно функций хг (Ь) и У (Ь) в (9) с помощью преобразования Лапласа к системе алгебраических уравнений относительно изображений Хг (р) и У (р) и разрешая полученную систему относительно Нг (р), находим ПФ для каждой массы рассматриваемой механической системы. Для верхней массы ПФ будет иметь вид
Нх (р) =
ат р'
+ ат- х рт 1 + ... + ахр + ао
52п р2п + 52п-1 р2п-1 + ... + 5х р + 5о
(11)
здесь
т < 2п — 1, 5о = ао, 5х = ах,
п — число масс изучаемой механической системы.
Подставляя в (11) вместо комплексной переменной р мнимую переменную ги, получим АЧХ и ее модуль. Значения коэффициентов аг и 5г определяются аналогично системе с двумя степенями свободы, рассмотренной выше.
Определение числа антирезонансных частот. Для линейной механической системы цепной структуры без трения в условиях кинематического возбуждения число антирезонансных частот на ¿-й массе (г = 1,п) равно числу собственных частот системы той же структуры с (г — 1) степенями свободы. Данное утверждение представим как формулу
Бг = г — 1, г = 1, п, (12)
где Бг
число антирезонансных частот на г-й массе.
Исходя из (12), число антирезонансных частот на каждой массе системы цепной структуры зависит от того, какое число масс находится выше нее. Нетрудно заметить, что на верхней массе такой системы антирезонансных частот нет. Введение в систему традиционной структуры многозвенных соединений в виде пружин приводит к появлению на верхней массе антирезонансных частот. Их число $1 будет принимать значения из интервала [1,п — 1] в зависимости от того, каким образом наложены эти упругие соединения.
Остановимся на случае, когда на систему цепной структуры наложено одно многозвенное упругое соединение. Рассмотрим следующие возможные варианты введения в систему такого соединения:
1) верхняя масса соединена с нижней массой то;
2) верхняя масса соединена с массой mj (з ^ 3);
3) масса т^ = 2, п) соединена с массой то;
4) масса т^ = 2, п) соединена с массой (к = 4, п, причем к ^ з + 2).
Для вычисления числа антирезонансных частот на верхней массе системы с многозвенным упругим соединением для каждого из перечисленных выше случаев применяем формулы
1) $1 = п — 1;
2) $1 = з — 2;
3) $1 = п — з;
4) $1 = (к — 1) — з.
(13)
(14)
(15)
(16)
Для примера рассмотрим вопрос наличия антирезонансных частот на верхней массе модели с тремя степенями свободы в зависимости от того, как наложено многозвенное соединение. На рис. 9 представлены возможные варианты.
Рис. 9. Механические системы с тремя степенями свободы с многозвенными соединениями
(объяснение в тексте)
вид
В первом случае (рис. 9, а) условие определения антирезонансных частот имеет
Со ((С1 + С2 ) — т2 и2 )((С2 + Сз ) — тз и2 ) — Со С2 + С1С2 Сз = 0. Решая данное уравнение при
(Со т2 (С2 + Сз) — тз (С1 + С2 ))2 > 4со т2 тз С2 С Сз — Со С2),
получаем 2 вещественных корня. Таким образом, число антирезонансных частот на верхней массе $1 = 2, что подтверждает справедливость формулы (13). Для случая, показанного на рис. 9, б, условие
СзСо((С1 + С2) - Ш2и2) + С1С2С3 = 0
определяет одну частоту антирезонанса, что соответствует формуле (14).
В третьем случае (рис. 9, в) условие определения антирезонансных частот имеет
вид
С1 Со ((С2 + Сз ) - тзи2) + С1С2Сз = 0
и $1 = 1, что подтверждает справедливость формулы (15).
Вариант механической системы, подтверждающей формулу (16), здесь не рассматривается.
Система с восемью степенями свободы. Для иллюстрации существенного влияния многозвенных соединений на АЧХ механической модели были проведены численные расчеты АЧХ механической системы с восемью степенями свободы. Эта система может служить механической моделью тела сидящего человека, если т1 — это масса головы, т0 — масса туловища, а массы т2,.т8 — массы семи шейных позвонков вместе с окружающими их мягкими тканями. На рис. 10 проиллюстрировано различие двух АЧХ такой системы при отсутствии и наличии многозвенных соединений, моделирующих многосуставные мышцы. Важно подчеркнуть, что основная цель этого рисунка была в том, чтобы обратить внимание на качественное, а не количественное отличие двух характеристик. В частности, видно, что наличие многозвенных соединений (кривая 2) может привести к появлению частот, при которых АЧХ близка к нулю, что в отсутствии этих соединений принципиально невозможно.
10.00
■ ' ' ...................................... ......П
0.00 5.00 10.00 15.00 20.00 /Гц
Рис. 10. АЧХ для модели цепной структуры с восемью степенями свободы в отсутствие дополнительных многозвенных соединений (1) и при их наличии (2)
Заключение. Модельное изучение влияния многосуставных мышц на частотные характеристики тела человека показало, что введение многозвенных соединений в структуру модели в отсутствии трения приводит к увеличению числа антирезонансных частот на всех массах системы. В частности, антирезонансные частоты
появляются и на верхней массе, что в отсутствии многозвенных соединений невозможно. Их число зависит от того, сколько многозвенных соединений наложено на многомассовую систему и какие именно массы ими соединены. В качестве иллюстрации влияния многосуставных мышц на частотные свойства тела сидящего человека приведены графики АЧХ механических систем с восемью степенями свободы при наличии и отсутствии многозвенных соединений.
Литература
1. Tregubov V. P. Problems of mechanical model identification for human body under vibration // Mechanism and Machine Theory. 2000. Vol. 35. P. 491-504.
2. Coermann R. R. The mechanical impedance of the human body in sitting and standing position at low frequencies // Human factors. 1962. Vol. 4. P. 227-253.
3. Griffin M. J. Vertical vibration of seated subjects: effects of posture, vibration level, and frequency // Aviation Space and Environmental Medicine. 1975. Vol. 46. P. 269-276.
4. Fairley T. E., Griffin M. J. The apparent mass of the seated human body: Vertical vibration // Journal of Biomechanics. 1989. Vol. 22. N 2. P. 81-94.
5. Потёмкин Б. А., Фролов К. В. Построение динамической модели тела человека-оператора, подверженного действию широкополосных случайных вибраций // Виброизоляция машин и виброзащита человека-оператора / отв. ред. К. В. Фролов. М.: Наука, 1973. C. 17-30.
6. Bai Xian-Xu, Xu Shi-Xu, Cheng Wei, Qian Li-Jun. On 4-degree-of-freedom biodynamic models of seated occupants: Lumped-parameter modeling // Journal of Sound and Vibration. 2017. Vol. 402. N 18. P. 122-141.
7. Wan Y., Schimmels J. M. A simple model that captures the essential dynamics of a seated human exposed to whole body vibration // Advances in Bioengineering. 1995. Vol. 31. ASME publ.-BED (Bioengineering Division). P. 333-334.
8. Abbas W., AbouelaMa, O. B., El-Azab M., Elsaidy M., Megahed A. A. Optimization of biodynamic seated human models using genetic algorithms // Engineering. 2010. Vol. 2. P. 710-719.
9. Boileau P. E., Rakheja S. Whole-body vertical biodynamic response characteristics of the seated vehicle driver: measurement and model development // International Journal of Ind. Ergonomics. 1998. Vol. 22. P. 449-472.
10. Zhang E., Xu L. A., Liu Z. H., Li X. L. Dynamic modeling and vibration characteristics of multi-DOF upper part system of seated human body // Chine Journal of Engineering Design. 2008. Vol. 15. P. 244-249.
11. Трегубое В. П., Егорова Н. К. О единственности решения задачи определения параметров механических моделей тела человека, подверженного вибрационному воздействию // Вестн. С.-Петерб. ун-та. Прикладная математика. Информатика. Процессы управления. 2019. Т. 15. Вып. 4. С. 565-577. https://doi.org/10.21638/11702/spbu10.2019.412
12. Bae J.-J., Kang N. Development of a five-degree-of-freedom seated human model and parametric studies for its vibrational characteristics // Shock and Vibration. 2018. Article ID 1649180. 15 p. https://doi.org/10.1155/2018/1649180
13. Yoshimura T., Nakai K., Tamaoki G. Multi-body dynamics modelling of seated human body under exposure to whole-body vibration // Industrial Health. 2005. Vol. 43. P. 441-447.
14. Po,njo,bi M. M., Andersson G. B., Jorneus L., Hult E., Mattsson L. In vivo measurements of spinal column vibrations // The Journal of Bone and Joint Surg. 1986. Vol. 68-A. P. 695-702.
15. Kitazaki S., Griffin M. A modal analysis of whole body vertical vibration using a finite element model of the human body // Journal of Sound Vibration. 1997. Vol. 200. P. 83-105.
16. Трегубое В. П., Селезнева Н. А. Математическое моделирование динамики шейного отдела позвоночника при импульсных воздействиях // Вестн. С.-Петерб. ун-та. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2016. Вып. 1. С. 53-65.
17. Хилл А. В. Механика мышечного сокращения. Старые и новые опыты / пер. с англ. Ю. А. Шаронова; под ред. Г. М. Франка. М.: Мир, 1972. 183 с. (Hill A. V. First and last experiments in muscle mechanics.)
18. Huxley A. F. Muscle structure and theories of contraction // Progr. Biophys. and Biophys. Chem. 1957. Vol. 7. P. 255-318.
19. Дещеревский В. И. Математические модели мышечного сокращения. М.: Наука, 1977. 160 с.
20. Tregubov V. P. Development of the muscle kinetic theory and cyclic movements // Acta of Bioengineering and Biomechanics. 2003. Vol. 5. Suppl. 1. P. 512-519.
21. Кликунова К. А., Трегубое В. П. Математическое моделирование переходных режимов
мышечного сокращения // Вестн. С.-Петерб. ун-та. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2008. Вып. 3. С. 56—62.
22. ТтвупЬву V. P., Azanchevsky V. V., Zarin A., Klikunova K. Mechanical models of the intervertebral disk and cervical spine // Lecture Notes of the ICB Seminars. Warsaw, 2007. Vol. 79. P. 41-54.
23. Глухарев К. К., Потёмкин Б. А., Фролов К. В. Особенности биодинамики тела человека при вибрациях // Виброзащита человека-оператора и вопросы моделирования / отв. ред. К. В. Фролов. М.: Наука, 1973. С. 22-28.
24. Санитарные нормы 2.2.4/2.1.8.566-96. М.: Министерство здравоохранения России, 1997. 35 с.
Статья поступила в редакцию 7 марта 2020 г. Статья принята к печати 28 мая 2020 г.
Контактная информация:
Трегубов Владимир Петрович — д-р физ.-мат. наук, проф.; v.tregubov@spbu.ru Егорова Надежда Константиновна — аспирант; nadezhda_ego@mail.ru
Model study of the influence of multi-joint muscles on the frequency characteristics of the human body
V. P. Tregubov, N. K. Egorova
St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Tregubov V. P., Egorova N. K. Model study of the influence of multi-joint muscles on the frequency characteristics of the human body. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2020, vol. 16, iss. 2, pp. 150-164. https://doi.org/10.21638/11701/spbu10.2020.207 (In Russian)
The analysis of previous works devoted to experimental studies of the frequency characteristics of the body of a sitting human body subjected to vibration and the construction of its mechanical models is conducted. It is noted that in most cases measurements of vibration are made on the seat and the head of the person that allows us to build a transmission function from the seat to the head. In doing so, the supposed internal structure of the human mechanical model remained unconfirmed. Under these conditions, it was emphasized that vibration measurements should be performed on all modeled parts of a human body. In addition, from an anatomical point of view, the main contribution to the mechanical properties of the body of a sitting person is made by the multi-jointed muscles of the spine. However, this important fact has not been taken into account yet. In this regard, the task was set to study the influence of multi-articular muscles on the frequency properties of the body of a sitting person. For this purpose, a number of mechanical models were constructed in which multi-jointed muscles were modeled by multi-link viscoelastic connections. In particular, on the simplest model with two degrees of freedom it was shown how the imposition of two-link connections in addition to a single-link one leads to the appearance of an anti-resonant frequency on the upper mass, that is impossible in their absence. A mechanical model with an arbitrary number of degrees of freedom with multi-link connections is given, for which the formulas of the transfer function and the amplitude-frequency response for the upper mass are obtained. In addition, we consider a mechanical model with an arbitrary set of multi-link connections. As an example, the results of numerical calculation of the frequency response for a mechanical model with eight degrees of freedom in the presence and absence of multi-link connections are given.
Keywords: mechanical model, human body, vibration, transfer function, input mechanical impedance, amplitude-frequency response, multi-joint muscles.
References
1. Tregubov V. P. Problems of mechanical model identification for human body under vibration. Mechanism and Machine Theory, 2000, vol. 35, pp. 491—504.
2. Coermann R. R. The mechanical impedance of the human body in sitting and standing position at low frequencies. Human factors, 1962, vol. 4, pp. 227—253.
3. Griffin M. J. Vertical vibration of seated subjects: effects of posture, vibration level, and frequency. Aviation Space and Environmental Medicine, 1975, vol. 46, pp. 269—276.
4. Fairley T. E., Griffin M. J. The apparent mass of the seated human body: Vertical vibration. Journal of Biomechanics, 1989, vol. 22, N 2, pp. 81—94.
5. Potemkin B. A., Frolov K. V. Postroenie dinamitheskoy modeli tela theloveka-operatora, podverzhennogo deistviju shirokopolosnih sluthainih vibratsiy [Construction of dynamic model of human body for man-operator exposed to broadband random vibrations]. Vibration-isolation of machine and vibration-protection of man-operator. By otv. red. K. V. Frolov. Moscow, Nauka Publ., 1973, pp. 17—30. (In Russian)
6. Bai Xian-Xu, Xu Shi-Xu, Cheng Wei, Qian Li-Jun. On 4-degree-of-freedom biodynamic models of seated occupants: Lumped-parameter modeling. Journal of Sound and Vibration, 2017, vol. 402, N 18, pp. 122-141.
7. Wan Y., Schimmels J. M. A simple model that captures the essential dynamics of a seated human exposed to whole body vibration. Advances in Bioengineering, 1995, vol. 31, ASME-publ. BED (Bioengineering Division), pp. 333-334.
8. Abbas W., Abouelatta O. B., El-Azab M., Elsaidy M., Megahed A. A. Optimization of biodynamic seated human models using genetic algorithms. Engineering, 2010, vol. 2, pp. 710-719.
9. Boileau P. E., Rakheja S. Whole-body vertical biodynamic response characteristics of the seated vehicle driver: measurement and model development. International Journal of Ind. Ergonomics, 1998, vol. 22, pp. 449-472.
10. Zhang E., Xu L. A., Liu Z. H., Li X. L. Dynamic modeling and vibration characteristics of multi-DOF upper part system of seated human body. Chine Journal of Engineering Design, 2008, vol. 15, pp. 244-249.
11. Tregubov V. P., Egorova N. K. O edinstvennosti resheniya zadachi opredeleniya parametrov mekhanicheskih modelej tela cheloveka, podverzhennogo vibracionnomu vozdejstviyu [On the uniqueness of the solution to the problem of determining the parameters of mechanical models of the human body exposed to vibration]. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2019, vol. 15, iss. 4, pp. 565-577. https://doi.org/10.21638/11701/spbu10.2019.412 (In Russian)
12. Bae J.-J., Kang N. Development of a five-degree-of-freedom seated human model and parametric studies for its vibrational characteristics. Shock and Vibration, 2018, Article ID 1649180, 15 p. https://doi.org/10.1155/2018/1649180
13. Yoshimura T., Nakai K., Tamaoki G. Multi-body dynamics modelling of seated human body under exposure to whole-body vibration. Industrial Health, 2005, vol. 43, pp. 441-447.
14. Panjabi M. M., Andersson G. B., Jorneus L., Hult E., Mattsson L. In vivo measurements of spinal column vibrations. The Journal of Bone and Joint Surgery, 1986, vol. 68-A, pp. 695-702.
15. Kitazaki S., Griffin M. A modal analysis of whole body vertical vibration using a final element model of the human body. Journal of Sound Vibration, 1997, vol. 200, pp. 83-105.
16. Tregubov V. P., Selezneva N. A. Matematicheskoe modelirovanie dinamiki shejnogo otdela pozvonochnika pri impul'snyh vozdejstviyah [Mathematical modeling of the dynamics of the cervical spine under impulse effects]. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2016, iss. 1, pp. 53-66. (In Russian)
17. Hill A. V. First and last experiments in muscle mechanics. Cambridge, Cambridge University Press, 1970, 180 p. (Russ. ed.: Hill A. V. Mehanika mishechnogo sokraschenia. Starie i novie opyty. Moscow, Mir Publ., 1972, 183 p.)
18. Huxley A. F. Muscle structure and theories of contraction. Progr. Biophys. and Biophys. Chem., 1957, vol. 7, pp. 255-318.
19. Descherevsky V. I. Matematicheskie modeli myshechnogo sokrashcheniya [Mathematical models of muscle contraction]. Moscow, Nauka Publ., 1977, 160 p. (In Russian)
20. Tregoubov V. P. Development of the muscle kinetic theory and cyclic movements. Acta of Bioengineering and Biomechanics, 2003, vol. 5, suppl. 1, pp. 512-519.
21. Klikunova K. A., Tregubov V. P. Matematicheskoe modelirovanie dinamiki shejnogo otdela pozvonochnika pri impul'snyh vozdejstviyah [Mathematical modeling of transient modes of muscle contraction]. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2008, iss. 3, pp. 56-62. (In Russian)
22. Tregubov V. P., Azanchevsky V. V., Zarin A., Klikunova K. Mechanical models of the intervertebral disk and cervical spine. Lecture Notes of the ICB Seminars. Warsaw, 2007, vol. 79, pp. 41—54.
23. Glukharev K. K., Potemkin B. A., Frolov K. V. Osobennosti biodinamiki tela cheloveka pri vibraciyah [Features of human body biodynamics under vibrations]. Vibration-protection of man-operator and the issues of modeling. By otv. red. K. V. Frolov. Moscow, Nauka Publ., 1973, pp. 22—28. (In Russian)
24. Sanitarnye normy 2.2.4/2.1.8.566-96 [Sanitary standards 2.2.4/2.1.8.566-96]. Moscow, Ministry of Health of Russia, 1997, 35 p. (In Russian)
Received: March 07, 2020.
Accepted: May 28, 2020.
Authors' information:
Vladimir P. Tregubov — Dr. Sci. in Physics and Mathematics, Professor; v.tregubov@spbu.ru
Nadezhda K. Egorova — Postgraduate Student; nadezhda_ego@mail.ru