АНАЛИЗ МЕХАНИЧЕСКИХ СВОЙСТВ ВОЛОС ЧЕЛОВЕКА С ПОМОЩЬЮ ГИПЕРУПРУГИХ МОДЕЛЕЙ МУНИ-РИВЛИНА Текст научной статьи по специальности «Медицинские науки и общественное здравоохранение»

Дата публикации: 01.06.2023 Publication date: 01.06.2023

DOI: 10.24412/2588-0500-2023_07_02_37 DOI: 10.24412/2588-0500-2023_07_02_37

УДК 612.799.1 UDC 612.799.1

АНАЛИЗ МЕХАНИЧЕСКИХ СВОЙСТВ ВОЛОС ЧЕЛОВЕКА С ПОМОЩЬЮ ГИПЕРУПРУГИХ МОДЕЛЕЙ МУНИ-РИВЛИНА С.А. Муслов1, С.Д. Арутюнов1, С.С. Перцов12, К.Г. Караков3

'Московский государственный медико-стоматологический университет имени А.И. Евдокимова, г. Москва, Россия

2ФГБНУ «НИИ нормальной физиологии им. П.К. Анохина», г. Москва, Россия

3ФГБОУ ВО «Ставропольский государственный медицинский университет» Минздрава России,

г. Ставрополь, Россия

Аннотация. Рассмотрены 2-, 3-, 5- и 9-параметрические гиперупругие модели Муни-Ривлина волос человека и определены их параметры. Наименьший коэффициент корреляции с опытными данными продемонстрировала 2-параметрическая модель (0,974), наибольший - 9-параметрическая (0,999). Начальный модуль Юнга 9-параметрической модели составил E0=2,84 ГПа, что по порядку величины совпадает с эмпирическими литературными данными. Результаты расчетов, полученные в системе компьютерной алгебры Mathcad, существенно отличаются от полученных с помощью комплекса ANSYS, что, очевидно, связано с различными алгоритмами вычислений, используемыми программами. При вводе данных в ПК ANSYS применяли как инженерные напряжения и деформации, так и истинные значения. Устойчивость моделей оценивали с помощью критерия Хилла и Друкера. Было установлено, что определенные ограничения на параметры моделей Cij Муни-Ривлина не выполняются, что было подтверждено анализом кривых a-L

Ключевые слова: волос, напряжения, деформации, гиперупругие модели, модель Муни-Ривлина, модуль упругости.

ANALYSIS OF MECHANICAL PROPERTIES OF HUMAN HAIR USING

HYPERELASTIC MOONEY-RIVLIN MODELS

S.A. Muslov1, S.D. Arutyunov1, S.S. Pertsov12, K.G. Karakov3

'Moscow State University of Medicine and Dentistry, Moscow, Russia

2Research Institute of Normal Physiology named after P.K. Anokhin, Moscow, Russia

3Stavropol State Medical University, Stavropol, Russia

Annotation. We have examined 2, 3, 5 and 9 parameter hyperelastic Muni-Rivlin models of human hair and identified their parameters. The lowest correlation coefficient with experimental data was demonstrated by the 2 parameter model (0.974), the highest - by the 9 parameter model (0.999). The initial Young's modulus ofthe 9 parameter model was Eo=2.84 GPa, which coincides with the empirical literature data in order of magnitude. The calculation results obtained in the Mathcad system differ significantly from those obtained using the ANSYS software, which is obviously due to the different calculation algorithms used by the programs. When entering data into the ANSYS PC, both engineering stresses, deformations and true values were applied. The stability of the models was evaluated using the Hill and Drucker stability criterion. It was found that certain restrictions on the parameters of the Mooney-Rivlin Cij models are not met, which was confirmed by the analysis of the o-A, curves.

Keywords: hair, stresses, deformations, hyperelastic models, Mooney-Rivlin model, modulus of elasticity.

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

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

эпидермиса, представляющие собой многослойные нитеобразные структуры из эпителиальных клеток (кератиноциты, меланоциты), ороговевающих по мере удаления из точки роста в волосяных фолликулах [1]. На голове человека около 150 тысяч волос, основной компонент волос -кератин. Здоровый волос весьма эластичен, в сухом состоянии может вытянуться на 20-30% от своей длины, а в мокром - более чем на 50% и затем вернуться в исходное состояние. Волокна человеческого волоса испытывают силы растяжения, когда за ними ухаживают и укладывают. Следовательно, поведение волос при натяжении представляет интерес для исследователей. Кроме того, эти знания могут быть полезны при алопеции для трансплантации и наращивании волос, изготовителям искусственных волос, в криминалистике, анимации и т.д. Но в отличие от механических

свойств кожи лица и тела человека, которые достаточно хорошо изучены (в зависимости от анатомического расположения, ориентации образцов относительно линий Лангера, возраста и пола, а также здоровья пациентов и даже косметических процедур), упругие и прочностные свойства волос исследованы недостаточно, и совсем нет данных о гиперупругих моделях данной биологической ткани. Как известно, к уникальным свойствам гиперупругих (резиноподобных) материалов относятся: способность подвергаться большим деформациям под нагрузкой, выдерживая деформации до 500% в технических приложениях; выраженное нелинейное поведение при нагрузке-растяжении; ответ на нагрузку почти без изменения объема. В связи с этим весьма актуально исследование гиперупругих характеристик волос человека с помощью существующих моделей.

Таблица 1

Тип кривой нагружения и модели Муни-Ривлина [10]

№

Тип кривой «напряжение-деформация»_

Рекомендуемый тип модели Муни-Ривлина

1.

2-параметрическая или 3-параметрическая модель

2.

3-параметрическая или 5-параметрическая модель

3.

5-параметрическая или 9-параметрическая модель; недостаток: модели требуют большего количества констант, описывающих материал

2-параметрическая модель

3-параметрическая модель

а

АА

aMR2(t)

Коэффициент деформации

А) SD=15,742, 5=8,233%, R=0,974, E0= 599,89 МПа

I АЛ

Н aMR3(t)

Коэффициент деформации

Б) SD=4,994, 5=2,915%, R=0,997, Ec=1,39 ГПа

250

250

0

0

Л, t

Л, t

5-параметрическая модель

9-параметрическая модель

S а

I АЛ

N aMR5(t)

I аа

й aMR9(t)

1 Л, t 1.44 1 Л, t 1.44

Коэффициент деформации Коэффициент деформации

В) SD = 3,37, 5 = 2,005%, R = 0,9985, Г) SD = 2,857, 5 =2,065%, R = 0,99889,

E0 = 1,92 ГПа E0 = 2,84 ГПа

Рис. 1. Сравнение эмпирических и расчетных данных 2- (А), 3- (Б), 5- (В) и 9- (Г) параметрических моделей Муни-Ривлина Примечание: SD - стандартное отклонение (standard deviation); 5 - максимальная приведенная погрешность; R - коэффициент корреляции; E0 - начальный модуль Юнга

250

250

200

а

100

0

Мы располагали экспериментальными кривыми «напряжение-деформация», полученными при одноосных испытаниях волос пациентов в возрасте от 20 до 29 лет [2]. Их вид существенно отличается от «обычных» J-кривых, присущих для большинства мягких тканей [3-5] в том числе и кожи [6]. Исследовали волосы головы, подмышечные и лобковые. Измеряли максимальное напряжение (кГ/мм2) и удлинение (%) образцов волос мужчин и женщин, из США и Японии, подвергнутых радиации и нет, сразу после удаления, через 6 и 12 месяцев. Из других

данных можно отметить работу [7], но в ней авторы, проведя детальное исследование оптическими методами, ограничились только определением начального упругого модуля и коэффициента Пуассона волос (3,727 ГПа и 0,377 соответственно) в проксимальной, средней и дистальной их части. Весьма информативны обзор исследований структуры и механики шерсти и волосяных волокон животных [8], а также экспериментально-расчетное исследование волос человека [9], где установлено, что начальный модуль Юнга равен 4,2 ГПа, а

устанавливающее уравнение деформации растяжения кератинового волокна имеет вид:

идентичны. Ошибку аппроксимации вычисляли с помощью опции Normalized error ANSYS. Использовали расчетные формулы,

g(s) = Eis + E2(8-8ci)H(8ci) + E3(s-8c2)H(sc2), приведенные в [10].

где Ei, E2, E3 - упругие модули диапазона упругой деформации, диапазона структурного a-ß перехода и постпереходного периода соответственно, 8ci = 0,02, 8c2 = ~0,25, H =1/2(1 + sgn x) - функция Хэвисайда.

Чтобы выбрать подходящую модель, мы использовали несколько моделей Муни-Ривлина: 2-, 3-, 5- и 9-параметрическую, получивших распространение в механике больших деформаций и доступных в литературе. Действительно, выбор определяющих соотношений и адекватных гиперупругих моделей является актуальной задачей современной механики. Несмотря на то, что в нашем случае исходная эмпирическая кривая о-8 имела одну точку перегиба (табл. 1, рис. 1), мы подвергли тестированию все существующие модели Муни-Ривлина.

Методы и организация исследования. Определение и исследование параметров гиперупругих моделей производили в системе компьютерной алгебры Mathcad 13.0 и многоцелевом программном комплексе (ПК) ANSYS 2022 R2. В пакете Mathcad применяли встроенные функции linfit, при расчете корреляции моделей -функцию corr. При вводе данных в ПК ANSYS использовали инженерные (условные) и истинные (Коши) напряжения и деформации. Как известно,

°ист = °и„А ^ист = ln(1 + ^инж )> ^= е ■+ 1 , где X - коэффициент деформации, 8 -относительная деформация соответственно. При малых деформациях значения инженерных и истинных параметров почти В 3 -параметрической модели:

Модель Муни-Ривлина является дальнейшим развитием неогуковской модели и одной из самых популярных среди всех гиперупругих моделей. Она была предложена М. Муни в 1940 году и выражена в терминах инвариантов Р. Ривлином в 1948 году. Чтобы описать поведение материала, модель Муни-Ривлина как функция инвариантов деформации с различными 2, 3, 5 или 9 параметрами, доступна в литературе, и подходящую модель можно выбрать в зависимости от типа кривой напряжения-деформации, как показано в таблице 1.

В 2-параметрической модели Муни-Ривлина функция плотности энергии деформированного несжимаемого материала (1з=1), каким предполагаются мягкие биологические ткани, определяется в виде линейной комбинации двух инвариантов правого тензора деформации Коши-Грина:

w2p = ад - 3)+c01(i2

3) =

(1)

где II и 12 - первый и второй инварианты тензора, равные I! = А + + , т - — + — + — соответственно, А* = 81 + 1

Т2 - 2 + л 2 + л 2 '

Л2 Л3

- главные компоненты, 81 - главные относительные деформации (1 = 1, 2, 3), а две материальные константы С10 и С01 имеют размерность напряжений (упругих модулей) и могут оцениваться по результатам испытаний. При этом начальный модуль сдвига для моделей Муни-Ривлина

Цо = 2(Сю + С01).

W3p = Сю(1х - 3) + CoX(I2 - 3) + Сп (Ix -1)(I2 -1) в 5-параметрической:

W5p = СД - 3) + CoX(l2 -3) + C2o(Ix -3)2 + Co2(I2 -3)2 + Сп& - 1)(I2 - 1),

(2)

БИОМЕДИЦИНЫ BIOMEDICINE 2023, T. 7 (2)_2023, Vol. 7 (2)

- для 5-параметрической:

= 2С10 (л--1-) + 2С01 (1 --1-) + 6Cn (л2-л-1+ -v + ±--1

л л л л л

12 11 +4Смл(1 - -)(л2 + - - 3) + 4C02 (2л + —2 " 3)(1 -

- для 9-параметрической:

а9р = 2Cio (Л--1-) + 2С(1 (1 - ) + 6C11 (Л2-Л-1 + + ) +

Л Л Л Л Л

(4)

в 9-параметрической:

= ад -3)+Со1(12 -3)+020(1! -3)2+с02(12 -3)2+сп& -1)(12 -1)+

+С21(11 -1)2(12 -1) + С12(11 -1)(12 -1)2 + С3о(11 -3)3 + Со3(12 - 3)3. Соответствующие выражения для одноосных условных напряжений имеют вид:

- для 2-параметрической модели:

с 1 1 1

^2р = 2(сш -^х^-^) = 2сю(х-^0 + 2со1 (1 (5)

- для 3-параметрической:

^3р = 2с10(х -^) + 2со1 (1 --1) + 6с11 (x2 -х-1 + ^ + ^-) (6)

(7)

+4С2(Л(1 - р)(Л2 + 2 - 3) + 4С(2 (2Л + ^ - 3)(1 - +

11,1 +2С21 (1 - -Г)(2Л + -у - 3)(2Л3 - 4Л + -2 +1) + (8)

Л Л Л

1 1 , 5 , 2 , 1

+2С12(1 - —)(2Л + — -3)(4Л2 - - - 3Л- 6) + 6С30(Л2 + 2- 3)2(Л-—) +

Л Л Л Л Л

+6С(3(2Л + ^2 - 3)2(1 - J,).

Для простоты анализа чаще всего применяется двухпараметрическая модель.

Результаты исследования и их обсуждение. Эмпирическая (треугольниками) и

модельные кривые, построенные на основании формул (5-8), представлены в рисунке 1, численные значения параметров моделей -в таблице 2.

Таблица 2

Численные значения параметров моделей Муни-Ривлина, МПа

Порядок модели C10 C01 C11 C20 C02 C21 C12 C30 C03

2-парам. 117 -17,1

3-парам. -9,87102 1,22103 3,78102

5-парам. -3.26103 3,58103 -3,40104 1,24104 2,52104

9-парам. -1,12104 1,17104 4,30108 -7,81107 -5,8010s 4,54.107 -4,56107 -6,94.106 -2,84107

Примечание: C10, C01, C11, C20, C02, C21, C12, C30, C03 ■ лей Муни-Ривлина

Из полученных данных следует, что 2-параметрической модели Муни-Ривлина соответствует наибольшее стандартное отклонение от опытных данных (15,74) и максимальная приведенная погрешность

материальные константы гиперупуругих моде-

(8,23%), при этом коэффициент корреляции - наименьший среди всех рассмотренных моделей (0,97). Наименьшее стандартное отклонение от экспериментальных данных (2,85) и минимальную приведенную

погрешность (2,06%) продемонстрировала 9-параметрическая модель, при этом коэффициент корреляции расчетных и эмпирических данных был наибольшим среди всех моделей и составил 0,999. Отметим, что уже 3-параметрическая модель имела высокие показатели стандартного отклонения, погрешности и корреляции данных в сравнении с 2-параметриче-ской. Таким образом, 9-параметрическая модель наилучшим образом описывает механические свойства волос человека.

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

Инженерные напряжения и деформации - 2-параметрическая модель: Cio = -47,96, Coi = 198,61; 3-параметрическая модель: Cio = -i,i8i03, Coi = i,43 i03, Cii = 4,45 i02; 5-параметрическая модель: Ci0 = -4,64i03, C0i = 5,0i i03, Cii = -5,6i i04; C20 = 2,03 i04, C02 = 4,i3 i04; 9-параметрическая модель: Ci0 = -i,i9 i04, C0i = i,24 i04, Cii = -4,52i08; C20 = 2,26 i08, C02 = 2,26 i08; C2i = -5,96 i06, Ci2 = -4,2i i07, C30 = 8,88i05, C03 = i,4ii07 МПа. «Невязка» (normalized error) составила для 2-, 3- и 5-параметрических моделей 2,64, 1,33 и 1,08 соответственно.

Истинные напряжения и деформации -2-параметрическая модель: Ci0 = 67,99, C0i = i02,90; 3-параметрическая модель: Ci0 = -i,80 i03, C0i = 2,ii i03, Cii = 8,92i02; 5-параметрическая модель: Ci0 = -7,49i03, C0i = 7,95 i03, Cii = -i,77 i05; C20 = 6,79 i04, C02 = i,2i i05; 9-параметрическая модель: Ci0 = -i,48 i04, C0i = i,54 i04, Cii = i,5ii07; C20=-7,8i i06, C02 = -7,28 i06; C2i = i,i9 i07, Ci2 = -i,52i07, C30=-2,83 i06, C03=7,80 i06 МПа. «Невязка» (normalized error) составила для 2-, 3- и 5-параметричесских моделей 2,05, 0,46 и 0,11 соответственно.

В связи со значительным отличием результатов моделирования в программах Mathcad и ПК ANSYS, а также при

использовании инженерных и истинных напряжений и деформаций в ПК ANSYS мы выполнили специальное исследование этого вопроса. Результаты представлены в рисунке 2. Как видно, все кривые, практически идентичные в исследованном интервале деформаций, сильно различаются вне его. Очевидно, что Mathcad и ANSYS используют разные алгоритмы при нахождении решения и «движутся» к нему по разным траекториям. Поэтому данные, полученные с помощью этих приложений, значительно различаются. Отметим, что переход в ANSYS на ввод истинных (Коши) напряжений и деформаций вместо инженерных (условных), рекомендованный в [11], не дал положительных результатов в плане уменьшения расхождений с данными Mathcad. На наш взгляд, различия в результатах расчетов Mathcad и ANSYS должны стать предметом дальнейшего отдельного изучения. Действительно, с одной стороны существуют вышеупомянутые рекомендации [11] и справочная информация ANSYS с другой [12], согласно которой при вводе экспериментальных данных в гиперупругие модели нужно использовать инженерные значения. По утверждению Т.В. Муриной и соавт. «... все напряжения в ANSYS должны быть заданы в условном виде, за исключением объемных, которые задаются в истинном виде» [13].

Модель Муни-Ривлина должна удовлетворять критерий устойчивости, чтобы воспроизводить реальное поведение материала. Критерий устойчивости может быть записан как:

TT * °> (9)

где Gij и 8ij - компоненты тензоров напряжений и деформаций соответственно. В литературной форме критерий означает, что деформации должны расти при увеличении соответствующих приложенных напряжений. В контексте энергии деформации, связанной с увеличением напряжения, она должна быть больше нуля согласно данному критерию.

Муни-Ривлина 2-параметр. (Mathcad)

А) Mathcad-ANSYS (инженерные напряжения и деформации)

Б) Mathcad-ANSYS (истинные напряжения и деформации)

В) ANSYS (инженерные напряжения и деформации) -ANSYS (истинные напряжения и деформации)

я С

к

g CTMR2_Mathcad(t) Й -

OMR2_ANSYS_hhx

а

К АА

g GMR2_Mathcad( t)

S -

aMR2_ANSYS_iicKt>

§ К

К Я

I

В! Оч

С

а

К

а

▲ А

aMR2_ANS YSnHaÄ t) aMR2_ANSYS_ramtl

-10 А ,t 10

Коэффициент деформации Муни-Ривлина 2-параметр. (ANSYS)

Коэффициент деформации Муни-Ривлина 2-параметр. (ANSYS)

Коэффициент деформации

Рис. 2.Сравнение расчетных кривых, полученных в системе Mathcad и ПК ANSYS (2-параметрическая модель Муни-Ривлина)

Таблица 3

Критерии стабильности моделей Муни-Ривлина_

Порядок модели Критерии Выполнение

2-параметрическая C10 + C01 > 0 +

С01 > 0 -

3 -параметрическая С10 + С01 > 0 +

С11 > 0 +

5-параметрическая С10 + С01 > 0 +

С20 > 0 +

С02 < 0 -

С20 + С02 + С11 > 0 +

9-параметрическая C10 + C01 > 0 +

С30 > 0 -

С03 < 0 +

С20 + С02 + С11 > 0 -

С30 + С03 + С12 + С21 > 0 -

Примечание: Cío, Coi, C11, C20, C02, C21, C12, C30, C03 - материальные константы гиперупуру-гих моделей Муни-Ривлина

Это условие, разработанное Hill R. при одноосной нагрузке, как указано

(1958) [14] и Drucker D.C. (1959) [15], ниже (табл. 3). Например, для 2-параметри-

предусматривает определенные ограниче- ческой модели при малых одноосных

ния на параметры моделей Муни-Ривлина деформациях:

а = 2C10 (1 + е - (1 - 2s)) + 2C01 (1 - (1 - 3s)) = 2C10 (3s) + 2C01 (3s) = 6s(C10 + C01)

da

откуда ds

— = 6(C10 + C01)

Здесь учтено, что Х=8+1 и разложение в ряд: (1 + 8)п =1+п8 при 8<<1.

Видно, что в ряде случаев некоторая часть неравенств не выполняется. Впрочем, это видно и из графиков (рис. 2), построенных для расширенного диапазона деформаций, хотя он отсутствует в исследованном деформационном интервале (рис. 1 ).

Отметим, что модельный начальный модуль Юнга Ео, определенный по графикам (рис. 1) как производная ¿оМХ (А=1), полностью совпал со значениями, найденными по формуле Ео~3цо=6(Сю+Со1), и составил 599,89 МПа, 1,39, 1,92 и 2,84 ГПа. Для 9-параметрической модели эти данные весьма близки к результатам, полученным опытным путем.

Заключение:

1. Рассмотрены 2-, 3-, 5- и 9-параметрические гиперупругие модели

Муни-Ривлина волос человека и определены их параметры.

2. Наименьший коэффициент корреляции с опытными данными и, соответственно, наибольшее стандартное отклонение продемонстрировала 2-параметрическая модель ^=15,742, R=0,974), наибольший коэффициент и наименьшее стандартное отклонение - 9-параметрическая модель ^=2,857, R=0,99889). 9-параметрическая модель наилучшим образом описывает механические свойства волос человека, однако подтверждено, что и 3-параметрической модели уже достаточно, чтобы описать кривую нагружения с одной точкой перегиба.

3. Начальный модуль Юнга 9-параметрической модели составил Ео=2,84 ГПа, что весьма близко к литературным данным, полученным опытным путем.

4. Данные, полученные в системе компьютерной алгебры Mathcad, существенно отличаются от полученных с помощью ПК ANSYS, что, очевидно, связано с различными алгоритмами вычислений, используемыми программами. При вводе данных в ПК ANSYS оперировали как с инженерными напряжениями и деформациями, так и с их истинными значениями.

5. Оценку устойчивости моделей Муни-Ривлина при одноосном нагружении

СПИСОК ЛИТЕРАТУРЫ

1. Дегтярев, В. П. Нормальная физиология /

B. П. Дегтярев, Н. Д. Сорокина. - Москва: ГЭОТАР-Медиа, 2019. - 480 с.

2. Yamada, H. Strength of Biological Materials / H. Yamada. - Baltimore, 1973. - 297 p.

3. Fung, Y. C. Biomechanics: Mechanical Properties of Living Tissues, 2nd Edition / Y. C. Fung. - Springer, June 18, 1993. - 586 p.

4. Муслов, С. А. Упругость и гиперупругость урогенитальных тканей человека и животных /

C. А. Муслов, Е. А. Лапшихина, Д. С. Кобзев // Эффективная фармакотерапия. Урология и нефрология. - 2021. - Т. 17. - № 25. - С. 6-24.

5. Измерение и расчет характеристик упругости стенки общего желчного протока человека / С. А. Муслов, Н. В. Зайцева, А. А. Корнеев, А. А. Синицын // Актуальные вопросы биологической физики и химии. - 2022. - Т. 7. - № 1. -С. 92-98.

6. Grebeniuk, L. A. The mechanical properties of human skin / L. A. Grebeniuk, A. A. Utenkin // Fiziologiia Cheloveka. - 1994. - № 20. - P. 157.

7. Measurement of Young's modulus and Poisson's ratio of Human Hair using Optical techniques / Hu Z., Li G., Xie H. [et al] // Article in Proceedings of SPIE - The International Society for Optical Engineering. - December, 2009. - 10 p. DOI: 10.1117/12.851415.

8. Hearle, J. W. S. A critical review of the structural mechanics of wool and hair fibres / J. W. S. Hearle // International Journal of Biological Macromo-lecules. - 2000. - № 27(2). - pp. 123-138.

9. Structure and mechanical behavior of human hair / Y. Yu, W. Yang, B. Wang, M. Meyers // Materials Science and Engineering: C. - 2017. -№ 73. - pp. 152-163.

10. Kumar, N. Hyperelastic Mooney-Rivlin model: Determination and physical interpretation of material constants / N. Kumar, V. V. Rao //

производили с помощью критерия Хилла и Друкера. Установлено, что определенные ограничения на параметры моделей не выполняются, что было подтверждено анализом кривых о-Х.

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

MIT International Journal of Mechanical Engineering. - 2016. - № 6(1). - pp. 43-46.

11. Browell, R. The Power of Nonlinear Materials Capabilities / R. Browell, G. Lin. - In: Ray Browell, ANSYS Product Line Manager. - ANSYS Solutions, 2000. - Vol. 2. - № 1. - 7 p.

12.Интернет-ресурс. Официальный сайт ANSYS. Раздел HELP - Mechanical APDL 2022 R2. Input the Hyperelastic Experimental Data. URL: https://ansyshelp.ansys.com/account/secured ?returnurl=/Views/Secured/corp/v222/en/ans_mat/ Hlp_G_STRmcfhyp.html ?q=hyperelastic. (дата обращения: 20.12.22)

13. Мурина, Т. А. Определение параметров гиперупругих моделей материалов в пакете ANSYS / Т. А. Мурина, В. В. Кузнецов, В. Н. Водяков // Сборник научных трудов международной научно-практической конференции «Энергоэффективные и ресурсосберегающие технологии и системы». Саранск, 2016. - С. 450-456.

14. Hill, R. General theory of uniqueness and stability in elastic-plastic solids / R. Hill // Journal of the Mechanics and Physics of Solids. - 1958. -№ 6(3). - pp. 236-249. DOI: 10.1016/0022-5096(58)90029-2.

15. Drucker, D. C. A definition of a stable inelastic material / D. C. Drucker // Journal of Applied Mechanics. - 1959. - № 26 (1). - pp. 101-195. DOI: 10.1115/1.4011929.

REFERENCES

1. Degtyarev V.P., Sorokina N.D. Normal Physiology. Moscow: GEOTAR-Media, 2019. 480 p. (in Russ.)

2. Yamada H. Strength of Biological Materials. Baltimore, 1973. 297 p.

3. Fung Y.C. Biomechanics: Mechanical Properties of Living Tissues, 2nd Edition. Springer, June 18, 1993. 586 p.

4. Muslov S.A., Lapshikhina E.A., Kobzev D.S. Elasticity and hyperelasticity of urogenital tissues of human and animals. Effective Pharmacotherapy. Urology and Nephrology, 2021, vol. 17, no. 25, pp. 6-24. (in Russ.)

5. Muslov S.A., Zajtseva N.V., Korneev A.A., Sinitsyn A.A. Measurement and calculation of characteristics describing elasticity of the common bile duct's wall in human. Relevant Issues of Biological Physics and Chemistry, 2022, vol. 7, no. 1, pp. 92-98. (in Russ.)

6. Grebeniuk L.A., Uten'kin A.A. The mechanical properties of human skin. I. Fiziologiia Cheloveka, 1994, no. 20, p. 157.

7. Hu Z., Li G., Xie H., Hua T., Chen P., Huang F. Measurement of Young's modulus and Poisson's ratio of Human Hair using Optical techniques. Article in Proceedings of SPIE - The International Society for Optical Engineering. December, 2009. 10 p. DOI: 10.1117/12.851415.

8. Hearle J.W.S. A critical review of the structural mechanics of wool and hair fibres. International Journal of Biological Macromolecules, 2000, no. 27(2), pp. 123-138.

9. Yu Y., Yang W., Wang B., Meyers M.A. Structure and mechanical behavior of human hair. Materials Science and Engineering: C, 2017. no. 73, pp.152-163.

10. Kumar N., Rao V.V. Hyperelastic Mooney-Rivlin model: Determination and physical interpretation of material constants. MIT

International Journal of Mechanical Engineering, 2016, no. 6(1), pp. 43-46.

11. Browell R., Lin G. The Power of Nonlinear Materials Capabilities. In: Ray Browell, ANSYS Product Line Manager. ANSYS Solutions, 2000, vol. 2, no. 1, 7 p.

12. Official ANSYS website. HELP section -Mechanical APDL 2022 R2. Input the Hyperelastic Experimental Data. Available at: https://ansyshelp. ansys.com/account/secured?returnurl=/Views/Secu red/corp/v222/en/ans_mat/Hlp_G_STRmcfhyp.ht ml?q=hyperelastic (accessed 20.12.22)

13. Murina V.A., Kuznetsov V.V., Vodikov V.N. Identifying parameters of hyperelastic material models with the ANSYS package. Collection of works of the International Scientific and Practical Conference "Energy Efficient and Resource Preserving Technologies and Systems. Saransk, 2016. pp. 450-456. (in Russ.)

14. Hill R. General theory of uniqueness and stability in elastic-plastic solids. Journal of the Mechanics and Physics of Solids, 1958, no. 6(3), pp. 236-249. DOI: 10.1016/0022-5096(58)90029-2.

15. Drucker D.C. A definition of a stable inelastic material. Journal of Applied Mechanics, 1959, no. 26 (1), pp. 101-195. DOI: 10.1115/1.4011929.

СВЕДЕНИЯ ОБ АВТОРАХ:

Сергей Александрович Муслов - кандидат физико-математических наук, доктор биологических наук, профессор кафедры нормальной физиологии и медицинской физики МГМСУ им. А.И. Евдокимова, Москва, e-mail: muslov@mail.ru, elibrary AuthorID: 185513, ORCID ID: 0000-0002-97526804.

Сергей Дарчоевич Арутюнов - Заслуженный врач РФ, Заслуженный деятель науки РФ, доктор медицинских наук, профессор, заведующий кафедрой цифровой стоматологии, МГМСУ им. А.И. Евдокимова, Москва, e-mail: sd_arutjunov@mail.ru.

Сергей Сергеевич Перцов - доктор медицинских наук, профессор, директор ФГБНУ «НИИ нормальной физиологии им. П.К. Анохина», Москва, e-mail: s_pertsov@mail.ru. Карен Григорьевич Караков - Заслуженный врач РФ, доктор медицинских наук, профессор, заведующий кафедрой терапевтической стоматологии ФГБОУ ВО «Ставропольский государственный медицинский университет» Минздрава России, Ставрополь, e-mail: terstomsgma@yandex.ru.

INFORMATION ABOUT THE AUTHORS:

Sergej Aleksandrovich Muslov - Candidate of Physical and Mathematical Sciences, Doctor of Biological Sciences, Professor of the Department of Normal Physiology and Medical Physics, Moscow State University of Medicine and Dentistry, Moscow, e-mail: muslov@mail.ru, elibrary AuthorID: 185513, ORCID ID: 0000-0002-9752-6804.

Sergej Darchoevich Arutyunov - Honored Doctor of Russia, Honored Scientist of Russia, Doctor of Medical Sciences, Professor, Head of the Department of Digital Dentistry, Moscow State University of Medicine and Dentistry, Moscow, e-mail: sd_arutjunov@mail.ru.

Sergej Sergeevich Pertsov - Doctor of Medical Sciences, Professor, Director, Research Institute of Normal Physiology named after P.K. Anokhin, Moscow, e-mail: s_pertsov@mail.ru.

Karen Grigor'evich Karakov - Honored Doctor of Russia, Doctor of Medical Sciences, Professor, Head of the Department of Therapeutic Dentistry, Stavropol State Medical University, Stavropol, e-mail: ter-stomsgma@yandex.ru.

Для цитирования: Анализ механических свойств волос человека с помощью гиперупругих моделей Муни-Ривлина / С. А. Муслов, С. Д. Арутюнов, С. С. Перцов, К. Г. Караков // Современные вопросы биомедицины. - 2023. - Т. 7. - № 2. DOI: 10.24412/2588-0500-2023_07_02_37

For citation: Muslov S.A., Arutyunov S.D., Pertsov S.S., Karakov K.G. Analysis of mechanical properties of human hair using hyperelastic Mooney-Rivlin models. Modern Issues of Biomedicine, 2023, vol. 7, no. 2. DOI: 10.24412/2588-0500-2023 07 02 37