Немавхола Ф., Панделани Т., Нгвангва Г. Применение гиперупругих моделей для описания поведения разных областей овечьего сердца на основе двухосных механических испытаний. Российский журнал биомеханики, 2022, № 2, С. 19-30. БО!: 10.15593/Н7ЬВюшеЬ/2022.2.02
РОССИИСКИИ ЖУРНАЛ БИОМЕХАНИКИ № 2,2022
RUSSIAN JOURNAL OF BIOMECHANICS
https://ered.pstu. ru/index.php/rjb
Научная статья
DOI: 10.15593/RZhBiomeh/2022.2.02 УДК 531/534: [57+61]
ПРИМЕНЕНИЕ ГИПЕРУПРУГИХ МОДЕЛЕЙ ДЛЯ ОПИСАНИЯ ПОВЕДЕНИЯ РАЗНЫХ ОБЛАСТЕЙ ОВЕЧЬЕГО СЕРДЦА НА ОСНОВЕ ДВУХОСНЫХ МЕХАНИЧЕСКИХ ИСПЫТАНИЙ
Ф. Немавхола1, Т. Панделани12, Г. Нгвангва1
1 Колледж науки, техники и технологии университета Южной Африки, Претория, Южная Африка
2 Совет по научным и промышленным исследованиям, Претория, Южная Африка
О СТАТЬЕ
АННОТАЦИЯ
Получена: 21 октября 2021 Одобрена: 10 июня 2022 Принята к публикации: 14 июня 2022
Ключевые слова:
механика сердечной деятельности, экспериментальная механика, механика овечьего сердца, гиперупругие модели, подбор параметров моделей, механика мягких тканей, двухо-севое механическое испытание.
Сердечная недостаточность остается одной из самых распространенных причин смерти во всем мире, особенно среди людей старше 60 лет. Чтобы разработать и подобрать подходящие материалы для замены тканей сердца с целью эффективного лечения, необходимо понять биомеханическое поведение сердца при нагрузке. В работе исследуется механический отклик пассивного миокарда овцы, полученного из трех разных областей сердца. Так как приобретение модели сердца живого животного стоит дорого и проведение экспериментальных исследований требует прохождения строгой этической экспертизы, авторы оценивают соответствие шести различных гиперупругих моделей по механическим испытаниям ткани пассивного миокарда. Использованы образцы сердечной ткани 10 овец, которые в течение 3 часов после смерти были доставлены в лабораторию для биомеханических испытаний. Верхние области сердца над короткой осью были аккуратно изъяты. Образцы тканей были взяты из центральных областей левого и правого желудочков, межжелудочковой перегородки. Затем эпикард и эндокард осторожно срезали, чтобы выделить миокард. Были получены, обработаны и дискретизированы кривые «напряжения-деформации». Результаты показывают, что модели Чои-Вито и Фанга наилучшим образом подходят для левого желудочка, а модели Хольцапфеля (2000), Хольцапфеля (2005), полиномиальная (анизотропная) и four-fiber family - для правого желудочка.
©ПНИПУ
Введение
По структуре и функции сердечная ткань является одной из самых сложных в организме человека и животных. Она является неоднородной и анизотропной [7; 9; 14; 16]. Сердце осуществляет свою работу как под действием активных (во время систолы), так и пассивных (в диастолу) сил, обеспечивая постоянное и
непрерывное движение крови по сосудам. Венозная кровь поступает в правый желудочек и выталкивается в легочную артерию для кровоснабжения сосудов легких. Во время сокращения левого желудочка кровь через аортальный клапан попадает в аорту и разносится по всем органам. Между правым и левым желудочками находится межжелудочковая перегородка, которая подвижна во время каждого сердечного цикла.
© Немавхола Фупуфело - профессор, e-mail: masitfi@unisa.ac.za : 0000-0002-6250-5157 © Панделани Заньяни - доцент, e-mail: tpandelani@csir.co.za : 0000-0003-2107-5194 © Нгвангва Гарри - доцент, e-mail: ngwanhm@unisa.ac.za : 0000-0001-5486-8049
Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
Различия в деятельности отделов сердца показывают, что они обладают разными биомеханическими свойствами, поскольку (как показано в исследовании [10]) функция и свойства взаимосвязаны. Golob и соавт. [10] сообщают, что многие формы сердечных заболеваний сопровождаются соответствующим ухудшением биомеханических характеристик ткани миокарда желудочков.
Поэтому понимание механических свойств мягких тканей имеет решающее значение для разработки адекватных численных моделей [23; 24; 28; 30-32]. Анализ биомеханики желудочка поможет описать его сложное биомеханическое поведение в различных условиях и поможет улучшить понимание сердечной функции и недостаточности [19-22]. Эти знания также смогут стать ключом к разработке эффективных методов лечения [25-27]. Однако большинство исследований были сосредоточены на изучении поведения тканей левой части сердца [15; 18; 33; 35; 36]. Миокард играет главную роль в насосной функции сердца, поэтому основная часть исследований уделяет внимание данному слою. Недавние исследования показали, что миокард правого желудочка обладает уникальными свойствами. Sacks и Chuong [34] показали, что структура правого желудочка отличается от левого с точки зрения жесткости направления волокон и степени анизотропии. Далее выявлено, что заболевания, связанные с легкими, например, COVID-19, по-видимому, оказывает большее влияние на ткань миокарда правого желудочка. Изученные ранее авторами различные свойства миокарда крыс с физиологической точки зрения недостаточно близки к анатомии сердца человека, чем более крупные животные.
В данном исследовании изучаются различия в ткани миокарда, выделенной из трех областей сердца овцы, путем изучения и сравнения шести гиперупругих моделей, а именно модели Фанга [6], модели Чои-Вито [15], модели Хольцапфеля (2000) [11], модели Холь-цапфеля (2005) [12], модели four-fiber family [2; 8] и полиномиальной (анизотропной) модели [3]. Исследованы особенности данных моделей для описания механического поведения образцов. В миокарде крыс как экспериментальные, так и вычислительные исследования показали, что тканьправого желудочка демонстрирует более явное анизотропное поведение, чем левый и межелудочковая перегородка, что согласуется с выводами Sacks и Chuong [34]. Сложную ориентацию волокон и высокую степень анизотропии в правом желудочке трудно точно описать с помощью гиперупругих моделей, которые предполагают только ортотропию в ткани. Экспериментальные исследования тканей миокарда живых животных имеют огромные ограничения из-за относительно высоких затрат и этических требований, чем применение математических моделей. С одной стороны, процедуры проведения испытаний и подготовка образцов на реальных животных очень трудоемкие. С другой стороны, компьютерное и мате-
матическое моделирование тканей сердца не требует сложной подготовки образцов, если известны параметры модели. Кроме того, математические модели предоставляют исследователям возможности для изучения различных физиологических состояний т без угрозы здоровья и жизни для моделируемых объектов.
В работе W. Li [17] было рекомендовано проводить дополнительные испытания на одноосное и двухосное растяжения в пограничных и удаленных тканях миокарда на разных временных масштабах. Можно добавить, что испытания должны проводится при большем числе вариантов нагружения, насколько это возможно, чтобы максимально реалистично описать механическое поведение ткани миокарда. Важно для понимания процесса развития инфаркта миокарда, с целью разработки эффективных стратегий лечения, оценить про-грессирования заболеваний до полной сердечной недостаточности.
Материалы и методы
Подготовка тканей
Сердца овец (п = 10) с неизвестными заболеваниями сердца и разного возраста были собраны на местной скотобойне через один час после забоя. Далее сердца были помещены в ящик-холодильник для доставки в лабораторию биомеханики. В течение часа в лаборатории овечьи сердца были подготовлены к механическим испытаниям. Все образцы хранили в 0,9%-ном физиологическом растворе NaCl в течение 30 мин перед приготовлением. После этого из сердец вырезались образцы квадратой формы 18 х18 мм2.
Двухосные механические испытания
Методология и протокол двухосных испытаний были в основном заимствованы из ранее опубликованных работ нашей исследовательской группы [28; 29]. Используя направление от основания к вершине и изолируя сосочко-вую мышцу, квадратные образцы размером 18 х 18 мм2 были вырезаны из левого и правого желудочков, а так же межжелудочковой перегородки. В данной работе продольным считалось направление, которое проходит вдоль сосочковой мышцы; окружным - направление под углом 90° к продольному. Декартова система координат (которая применяется для микроструктурного анализа) в данном случае не использовалась, поскольку визуализация для определения направления волокон не выполнялась.
Вместо этого была рассмотрена цилиндрическая система координат. Двухосные испытания механических свойств всех образцов проводились с использованием аппарата BioTester 5000 CeUScale, Ватерлоо, Канада (рис. 1). Квадратный образец ткани миокарда был помещен в захваты BioRake, таким образом его размер уменьшался до
площади 16 х16 мм2. Затем во время испытаний к каждо- му образцу прикладывали предварительную нагрузку
Рис. 1. Экспериментальная установка для равноосного механического испытания сердца овцы для трех различных областей, включая левый желудочек, межжелудочковую перегородку и правый желудочек
Таблица 1
Список гиперупругих моделей, использованных в исследовании, для сопоставления данных о двухосевом растяжении различных областей/стенок сердца овец (левый, правый желудочки и межжелудочковая перегородка)
Номер модели Название модели Математическая формулировка Источники
1 Модель Фанга ш = с е -1), (1) где Q = ЬЕЦв + Ь2Е222 + ЬъЕ2т + 2Ь4£йв£22 + 2Ь5Е22ЕХХ + 2Ь6ЕххЕвв ; Ь - фшичтекж параметры. Модель реализована в полиномиальном формате [6]
2 Модель Чои-Вито Ш = Ь0 [ехр(ЬЕ21) + ехр (Ь2Е222) + ехр (2ЬзЕпЕ22) - з], (2) где Ь, - физический параметр. Модель реализована в экспоненциальном формате [15]
3 Модель Хольцапфеля (2000) Ш = ^ехр(с2(14 -1)2)-1], (3) где с, - физический параметр. Модель использована в экспоненциальном формате [11]
4 Модель Хольцапфеля (2005) Ш = С. {ехр [С2 ((1 - к)(/! - 3)2 + к(14 -1)2) -1], (4) где с, - физический параметр и к - параметр, который регулирует скорость сходимости. Модель указана в экспоненциальном формате [12]
5 Модель four-fiber family Ш = §(/1 -3) + 14=1^{ехр|>2,-1)2]-1)' (5) модель задана в гибридной форме (сочетание полиномиальной и экспоненциальной форм), где с, си, С21 - физические параметры [2, 8]
6 Полиномиальная (анизотропная) модель ш=ез=1«, (/1 - з)'+ЕЗ=А (/2 - з)у+Е:=2 с (/4 -1)к+Е:=2 (/б -1)т, (6) где а,, Ь, ск, и вт - физические параметры [3]
0,5 мН. Чтобы уменьшить созданные напряжения, все образцы были подвергнуты предварительной обработке, по методике описанной в предыдущем исследовании [29]. Каждый образец нагружался равномерно в продольном и окружном направлениях до 0,4% деформации. Чтобы имитировать температуру тела, образец помещали в 0,9%-ный физиологический раствор NaCl и нагревали до 37°С перед двухосными механическими испытаниями. Площадь поперечного сечения определялась путем измерения толщины образца с помощью штангенциркуля.
Гиперупругие модели
Обычно для представления анизотропии в биологических мягких тканях используются два подхода: подход, с
использованием тензора Грина - Лагранжа [13; 25-27], и подход, с использованием инвариантов деформации. Подходы, основанные на применении тензора Грина - Лагранжа, включают выражение функций плотности энергии деформации W в виде суммирования компонент тензора Eij. Ateshian и Costa [1] утверждают, что эта формулировка позволяет разделить функцию энергии деформации на шаровую и девиаторную части, что облегчает учет несжимаемости при реализации численной процедуры.
Однако в работе Chagnon и соавт. [4] сообщают, что такие модели обладают свойствами, которые влияют на их поливыпуклость и делают их менее устойчивыми. Другая трудность, связанная с этими моделями, заключается в том, что параметры материала не имеют физического
Результаты
Обнаружено, что толщины левого, правого желудочков и межжелудочковой перегородки составляют 14,5 ± 0,85 мм; 7,2 ± 0,25 мм; 5,2 ± 0,13 мм соответственно. Они были использованы для расчета достоверных напряжений и деформаций как в продольном, так и в окружном направлениях. Для оценки напряжений, площадь поперечного сечения была рассчитана с использованием толщины и ширины (16 ± 0,012 мм) образца ткани.
Табл. 2-7 содержит информацию о параметрах материала шести гиперупругих моделей, а также статистические метрики. Как и ожидалось, разброс значений параметров в параметрах материала у десяти образцов достаточно большой, что свидетельствует о вероятностном характере получившихся свойств материала миокарда у разных образцов.
Поэтому очень трудно применить набор параметров материала, полученных на одном образце, к экспериментальным данным, полученным на другом образце, даже для одной и той же модели. Тем не менее обнадеживает то, насколько малы различия в параметрах соответствия моделей (например, коэффициент детерминации Я2) для всех исследуемых моделей.
Таблица 2
Данные двухосного механического испытания, аппроксимированные гиперупругой полиномиальной (анизотропной) моделью (я1, Я2, Я3, ¿1, ¿2, ¿3, С2, С3, С4, С5, С6, ф). Средние значения получены по параметрам
материала для каждого образца
смысла, и, соответственно такие модели сложнее апрок-симировать.
С другой стороны, в подходах, основанных на инвариантах деформации 4, функция плотности энергии деформации Ш, представляется в виде комбинации изотропных и анизотропных функций. В литературе существуют три различные реализации: полиномиальная, степенная и экспоненциальная, причем последняя является наиболее популярной, поскольку она описывает эффект упрочнения при деформации в мягких тканях. В литературе существуют другие методы, с использованием инвариантов деформации как функций логарифма и тангенса. Однако эти формы подходят для моделирования небольших деформаций, особенно в период активации мягких тканей. Значимой особенностью подходов, основанных на применении инвариантов деформации, является учет анизотропии за счет 14 и/или 16 инвариантов деформации.
В исследовании рассматриваются шесть анизотропных гиперупругих моделей, две из которых основаны с использованием тензора Грина - Лагранжа, а четыре - на подходе инварианта деформации. Выражения для функций плотности энергии деформации приведены в табл. 1.
Параметры Левый желудочек Межелудочковая перегородка Правый желудочек
Коэффициент детерминации Я2 0,974 ± 0,023 0,948 ± 0,060 0,986 ± 0,009
Коэффициент корреляции г 0,117 ± 0,043 0,152 ± 0,108 0,107 ± 0,045
Стандартное отклонение МБ 0,989 ± 0,009 0,978 ± 0,026 0,994 ± 0,004
Среднеквадратичное отклонение МЯЫББ 0,097 ± 0,038 0,118 ± 0,082 0,084 ± 0,033
а\ 0,787 ± 0,397 0,516 ± 0,219 0,701 ± 0,268
а2 0,496 ± 0,681 0,187 ± 0,382 0,670 ± 0,506
аз 0,044 ± 0,296 0,143 ± 0,121 0,333 ± 0,288
Ьх 0,019 ± 0,050 -0,019 ± 0,031 0,000 ± 0,060
Ь2 0,006 ± 0,040 0,008 ± 0,027 -0,003 ± 0,043
Ьз 0,027 ± 0,090 -0,019 ± 0,060 -0,024 ± 0,034
С2 -0,212 ± 0,250 0,051 ± 0,198 -0,122 ± 0,072
С3 -0,014 ± 0,131 0,171 ± 0,205 -2,007 ± 5,971
С4 0,177 ± 0,140 0,084 ± 0,094 0,034 ± 0,063
С5 0,006 ± 0,099 -0,072 ± 0,270 -0,030 ± 0,108
С6 -0,038 ± 0,188 -0,118 ± 0,112 0,029 ± 0,061
Ф 0,027 ± 0,068 -0,016 ± 0,034 -0,017 ± 0,098
Таблица 3
Данные двухосного механического испытания, аппроксимированные гиперупругой моделью Фанга (с, ¿1, ¿2, ¿3, ¿4, ¿5, ¿6). Средние значения получены по параметрам материала для каждого образца
Параметры Левый желудочек Межелудочковая перегородка Правый желудочек
Коэффициент детерминации Я2 0,978 ± 0,012 0,921 ± 0,122 0,953 ± 0,045
Коэффициент корреляции г 0,990 ± 0,005 0,961 ± 0,070 0,984 ± 0,023
Стандартное отклонение МБ 0,126 ± 0,031 0,162 ± 0,159 0,145 ± 0,099
Среднеквадратичное отклонение МЯЫББ 0,150 ± 0,037 0,201 ± 0,177 0,184 ± 0,121
С 2,199 ± 1,380 35,178 ± 47,660 30,275 ± 63,779
Ь1 2,721 ± 0,580 0,301 ± 0,307 0,678 ± 2,536
Ь2 0,186 ± 0,352 0,126 ± 0,471 0,420 ± 0,460
Ь3 -1,098 ± 0,509 -0,070 ± 0,447 0,293 ± 1,972
Ь4 -0,531 ± 0,307 0,312 ± 0,339 0,283 ± 1,847
Ь5 1,331 ± 0,928 -0,247 ± 0,358 -1,204 ± 2,281
Ьб -1,523 ± 0,602 -0,194 ± 0,238 -0,066 ± 0,965
Таблица 4
Данные двухосного механического испытания, аппроксимированные гиперупругой моделью four-fiber family (c, cii, C2i, ci2, C22, ci34, C234, 90). Средние значения получены по параметрам материала для каждого образца
Параметры Левый желудочек Межелудочковая перегородка Правый желудочек
Коэффициент детерминации Я2 0,977 ± 0,010 0,943 ± 0,052 0,986 ± 0,052
Коэффициент корреляции г 0,132 ± 0,027 0,147 ± 0,093 0,084 ± 0,093
Стандартное отклонение МБ 0,16 ± 0,035 0,178 ± 0,108 0,105 ± 0,108
Среднеквадратичное отклонение МЯЫББ 0,99 ± 0,005 0,975 ± 0,024 0,994 ± 0,024
С 2,075 ± 0,832 0,328 ± 0,265 0,211 ± 0,265
С11 97,187 ± 304,567 0,665 ± 0,966 1,172 ± 0,966
С21 0,924 ± 0,454 1,357 ± 0,648 1,030 ± 0,648
С12 2,357 ± 1,434 1,98 ± 1,317 0,731 ± 1,317
С22 1,072 ± 0,324 0,539 ± 0,574 1,139 ± 0,574
С134 2,132 ± 0,782 0,792 ± 0,670 1,599 ± 0,670
С234 1,765 ± 0,623 0,908 ± 1,084 0,974 ± 1,084
ф0 0,85 ± 0,560 0,981 ± 0,723 0,895 ± 0,723
Таблица 5
Данные двухосного механического испытания, аппроксимированные гиперупругой моделью Чои-Вито (¿0, ¿1, ¿2, ¿3). Средние значения получены по параметрам материала для каждого образца
Параметры Левый желудочек Межелудочковая перегородка Правый желудочек
Коэффициент детерминации Я2 0,938 ± 0,021 0,772 ± 0,237 0,853 ± 0,146
Коэффициент корреляции г 0,980 ± 0,006 0,969 ± 0,030 0,976 ± 0,035
Стандартное отклонение МБ 0,174 ± 0,043 0,303 ± 0,214 0,257 ± 0,145
Среднеквадратичное отклонение МЯЫББ 0,207 ± 0,048 0,367 ± 0,241 0,315 ± 0,201
Ьс 73,594 ± 19,946 101,626 ± 96,351 117,659 ± 126,761
Ь1 0,018 ± 0,024 0,126 ± 0,190 0,731 ± 0,900
Ь2 0,040 ± 0,059 0,097 ± 0,174 1,184 ± 1,347
Ь3 0,032 ± 0,035 0,312 ± 0,446 1,880 ± 2,117
Таблица 6
Данные двухосного механического испытания, аппроксимированные гиперупругой моделью Хольцапфеля (2000) (и, к1, к2, ф ). Средние значения получены по параметрам материала для каждого образца
Параметры Левый желудочек Межелудочковая перегородка Правый желудочек
Коэффициент детерминации Я2 0,96 ± 0,024 0,068 ± 0,068 0,984 ± 0,014
Коэффициент корреляции г 0,129 ± 0,037 0,116 ± 0,116 0,089 ± 0,031
Стандартное отклонение МБ 0,156 ± 0,046 0,127 ± 0,127 0,112 ± 0,039
Среднеквадратичное отклонение МЯЫББ 0,985 ± 0,008 0,033 ± 0,033 0,994 ± 0,004
ц 0 ± 0,000 0,015 ± 0,015 0,002 ± 0,007
к1 1,501 ± 0,554 0,408 ± 0,408 1,429 ± 0,691
кг 0,397 ± 0,340 0,258 ± 0,258 0,874 ± 0,533
Ф 0,912 ± 0,434 0,389 ± 0,389 0,889 ± 0,082
Таблица 7
Данные двухосного механического испытания, аппроксимированные гиперупругой моделью Хольцапфеля (2005) (и к1, к2, ф, р). Средние значения получены по параметрам материала для каждого образца
Параметры Левый желудочек Межелудочковая перегородка Правый желудочек
Коэффициент детерминации Я2 0,966 ± 0,022 0,953 ± 0,057 0,986 ± 0,013
Коэффициент корреляции г 0,115 ± 0,040 0,121 ± 0,102 0,081 ± 0,029
Стандартное отклонение МБ 0,142 ± 0,050 0,156 ± 0,117 0,102 ± 0,038
Среднеквадратичное отклонение МЯЫББ 0,987 ± 0,008 0,978 ± 0,028 0,994 ± 0,004
Ц 0,032 ± 0,054 0,073 ± 0,131 0,109 ± 0,140
к1 1,412 ± 0,491 1,216 ± 0,455 1,295 ± 0,650
к2 0,226 ± 0,278 0,1 ± 0,206 0,552 ± 0,665
Ф 0,83 ± 0,814 0,681 ± 0,481 1,274 ± 0,351
Р 0,756 ± 0,248 0,725 ± 0,162 0,611 ± 0,348
1,2 1,2
1,0 1,0 0,8
0,8 —с^тгус- "" Л А
0,6 0,4 0,2 ЛЖ % МЖП ПЖ 0,6 0,4 0,2 ЛЖ МЖП ПЖ
0 -1-1-1-1-1-1-1-1-1-1 0 1 1 1 1 1 1 1 1 1 1
(Ч
1,2 1,0 0,8 0,6 0,4 0,2 0
1,2 1,0 0,8 % 0,6 0,4 0,2 0
1 2 3
4 5 6 Образцы а
7 8 9 10
1 2 3
4 5 6 Образцы в
7 8 9 10
ЛЖ
МЖП
ПЖ
1 2 3
4 5 6 Образцы д
7 8 9 10
1 2 3
4 5 6 Образцы б
7 8 9 10
ЛЖ
МЖП
ПЖ
(Ч
о;
1,2 1,0 0,8 0,6 0,4 0,2 0
1,2 1,0 0,8 0,6 0,4 0,2 0
ЛЖ
МЖП
ПЖ
1 2 3 4
5 6 7 Образцы
г
8 9 10
ЛЖ
МЖП
ПЖ
1 2 3 4
5 6 7 Образцы
е
8 9 10
Рис. 2. Коэффициент детерминации R2 шести гиперупругих моделей в 10 образцах, рассмотренных при двухосевом испытании ткани сердца овец в областях: ЛЖ - левый желудочек, МЖП - межжелудочковая перегородка, ПЖ - правый желудочек. а - модель Чои-Вито, б - модель Фанга, в - полиномиальная (анизотропная) модель, г - модель four-fiber family, д - модель Хольцапфеля (2000), е - модель Хольцапфеля (2005)
На рис. 2 показаны различия в коэффициенте детерминации Я2 шести гиперупругих моделей у 10 образцов.
На рис. 3 показан индекс оценки, на основе фактического коэффициента детерминации Я2 шести гиперупругих моделей у 10 образцов.
Обсуждение
В исследовании были определены характеристики подгонки шести моделей к различным областям сердца (левый и правый желудочки, межжелудочковая перегородка). Прямое сравнение гиперупругих моделей было проведено на основе коэффициента детерминации Я2 и индекса оценки для левого и правого желудочков, а также межжелудочковой перегородки.
Коэффициенты Я2 усредненных кривых левого, правого желудочков и межжелудочковой перегородки для шести различных моделей представлены в табл. 8. Таким образом, было обнаружено, что модели Чои-Вито и Фанга наилучшим образом подходят для левого желу-
дочка, а модели Хольцапфеля (2000), Хольцапфеля (2005), полиномиальная (анизотропная) и four-fiber family - для правого желудочка.
Таблица 8
R2 значения для всех стенок сердца (левый, правый желудочки и межжелудочковая перегородка)
Параметры Левый желудочек Межжелудочковая перегородка Правый желудочек
Модель Фанга 0,98 0,92 0,95
Полиномиальная модель 0,97 0,95 0,99
Модель Хольцапфеля (2000) 0,96 0,94 0,98
Модель Хольцапфеля (2005) 0,97 0,95 0,99
Модель four-fiber family 0,98 0,94 0,99
Модель Чои-Вито 0,94 0,77 0,85
100
80
I 60
Я
0
1 40
20
0
-20 120
100
80
60
40
20
0
-20 120
100
80
60
40
20
0
-20
я о
и Я о
Полиномиальная модель Модель Фанга Модель Хольцапфеля (2000) Модель Хольцапфеля (2005) Модель four-fiber family Модель Чои-Вито
5 6 Образцы а
5 6 Образцы б
■ - Полиномиальная модель п - Модель Фанга г - Модель Хольцапфеля (2000) 1 - Модель Хольцапфеля (2005) а- Модель four-fiber family Модель Чои-Вито
я - Полиномиальная модель в - Модель Фанга □ - Модель Хольцапфеля (2000) о - Модель Хольцапфеля (2005) а - Модель four-fiber family в - Модель Чои-Вито
5 6 Образцы
10
Рис. 3. Индекс оценки, определенный на основе фактического коэффициента детерминации R2 полиномиальной модели, модели Фанга, модели Хольцапфеля (2000), модели Хольцапфеля (2005), модели four-fiber family, модели Чои-Вито для левого и правого желудочков, межжелудочковой перегородки
в
1 г»г»
ЛЖ
МЖП а
ПЖ
1—=—1 ' _1_ '
1
1
1,00 0,90 0,80 0,70 0,60 0,50 0,40 0,30 0,20
1,00 0,98 0,96 0,94 0,92 0,90 0,88 0,86 0,84
ЛЖ
МЖП в
ПЖ
ЛЖ
ПЖ
ЛЖ
МЖП
г
ПЖ
1,00
0,95
0,90
0,85
0,80
0,75
МЖП
д
Рис. 4. Распределение Я2 для областей левого, правого желудо симации рассматриваемыми гиперупругими моделями: а - мод дель Хольцапфеля (2000), г - модель Хольцапфеля (2005)
Именно коэффициент детерминации Я2 использовался в качестве критерия, однако сообщается, что усреднение кривой напряжение-деформация является наилучшим вариантом [16]. Кроме того, выполнено сравнение коэффициента детерминации Я2, который получен для каждого образца (п = 10) в шести гиперупругих моделях. Этот подход может помочь сделать правильный вывод, поскольку на графике прослеживается тенденция по выбранным моделям на основе Я2.
На рис. 4 показаны диаграммы шести гиперупругих моделей при сравнении значения Я2 для областей левого и правого желудочков, а также межжелудочковой перегородки образцов овец. При использовании Я2 отмечено, что наблюдались значительные отклонения в поперечном направлении в левом желудочке - правом желудочке, незначительные в межжелудочковой
перегородке - правом желудочке и незначительные в левом желудочке - межжелудочковой перегородке, как показано в табл. 9.
Для модели Фанга различия в Я2 между левым желудочком и межжелудочковой перегородкой оказались значительными (р = 0,155), а между левым и правым желудочками были незначительными (р = 0,107). Кроме того, не было обнаружено никаких заметных различий между межжелудочковой перегородкой и правым желудочком (р = 0,436).
Для модели Чои-Вито различия в Я2 между левым желудочком и межжелудочковой перегородкой оказались значительными (р = 0,041), а между левым и правым желудочками - незначительными (р = 0,085). Кроме того, не было обнаружено никаких заметных разли-
чий между межжелудочковой перегородкой и правым желудочком (р = 0,370).
Для полиномиальной модели различия в R2 между левым желудочком и межжелудочковой перегородкой оказались незначительными (р = 0,213), а между левым желудочком и правым желудочком - незначительными (р = 0,152). Кроме того, не было обнаружено никаких заметных различий между межжелудочковой перегородкой и правым желудочком (р = 0,067).
Для модели four-fiber family различия в R2 между стенками левого желудочка и межжелудочковой перегородкой оказались менее значимыми (р = 0,063). Имелись достоверные различия между левым и правым желудочками (р = 0,014), а также между межжелудочковой перегородкой и правым желудочком (р = 0,019).
Для модели Хольцапфеля (2000) обнаружено, что различия в R2 между левым желудочком и межжелудочковой перегородкой были незначительными (p = 0,376), но значимыми были различия между левым и правым желудочками (p = 0,017). Кроме того, были обнаружены меньшие различия между межжелудочковой перегородкой и правым желудочком (р = 0,061).
В модели Хольцапфеля (2005) различия в R2 между левым желудочком и межжелудочковой перегородкой
оказались незначительными (p = 0,527), но значимыми - между левым и правым желудочками (p = 0,024). Кроме того, выявлены меньшие различия между межжелудочковой перегородкой и правым желудочком (р = 0,095).
Для модели Фанга различия в R2 между левым желудочком - межжелудочковой перегородкой, левым желудочком - правым желудочком и межжелудочковой перегородкой - правым желудочком были признаны незначительными (p = 0,155; 0,1066; 0,436 соответственно).
Для левого желудочка обнаружено, что модель Фанга имеет индекс оценки 100 (рис. 5). Значит, биомеханическое поведение этой модели наиболее близко по сравнению с другими пятью гиперупругими моделями. Выявлено, что модель Хольцапфеля (2005) имеет оценочный индекс 100 и, следовательно, считается наилучшей при установке на желудочковую стенку сердца овцы. Для правого желудочка модель four-fiber family показала наилучшие результаты с оценочным индексом 100. Модель Чои-Вито получила худшие показатели в трех областях сердца, и свидетельствует, что поведение ткани миокарда не может быть точно представлено моделями, которые не учитывают ориентацию волокон и анизотропию.
Фанга альная мо- Хольцапфеля Хольцапфеля four-fiber Чои-Вито дель (2000) (2005) family
Рис. 5. Индекс оценки вычисленный по среднему значению R2 полиномиальной модели, модели Фанга, модели Хольцапфеля (2000), модели Хольцапфеля (2005), модели four-fiber family, модели Чои-Вито (ЛЖ - левый желудочек, МЖП - межжелудочковая перегородка, ПЖ - правый желудочек)
Таблица 9
Р-значения для R2 (ЛЖ- левый желудочек, МЖП - межжелудочковая перегородка, ПЖ - правый желудочек)
Модель P-значения для R2: поперечная вариабельность всех областей сердечной стенки
ЛЖ-МЖП ЛЖ-ПЖ МЖП-ПЖ
Модель Чои-Вито 0,041 0,085 0,370
Полиномиальная модель 0,213 0,152 0,067
Модель four-fiber family 0,063 0,014 0,019
Модель Хольцапфеля (2000) 0,376 0,017 0,061
Модель Хольцапфеля (2005) 0,527 0,024 0,095
Модель Фанга 0,155 0,1066 0,436
Заключение
Настоящее исследование является важным шагом на пути к надежному структурному моделированию здорового миокарда для описания механических свойств различных тканей в условиях двухосевой нагрузки. Результаты двухосевых механических испытаний, проведенных на ткани сердца овцы из левого, правого желудочков и межжелудочковой перегородки, были сопоставлены с шестью гиперупругими моделями.
Настоящее исследование продемонстрировало важность выбора правильных моделей для различных областей сердца овцы.
Это имеет решающее значение, поскольку данные
Список литературы
1. Ateshian G.A., Costa K.D. A frame-invariant formulation of Fung elasticity // Journal of Biomechanics. - 2009. -Vol. 42(6). - P. 781-785. DOI: https://doi.org/10.1016 /j.jbiomech.2009. 01.015
2. Baek, S., Gleasona R.L., Rajagopal K.R., Humphreya J.D. Theory of small on large: potential utility in computations of fluid-solid interactions in arteries // Computer Methods in Applied Mechanics and Engineering. - 2007. - Vol. 196(31-32). - P. 3070-3078. DOI: https://doi.org/10.1016 /j.cma.2006.06.018
3. Bursa, J., Skacel P., Zemanek M., Kreuter D. Implementation of hyperelastic models for soft tissues in FE program and identification of their parameters // Conference: Proceedings of the Sixth IASTED International Conference on Biomedical Engineering. -Innsbruck: Austria. - 2008.
4. Chagnon G., Rebouah M., Favier D. Hyperelastic energy densities for soft biological tissues: a review // Journal of Elasticity. - 2015. - Vol. 120(2). - P. 129-160. DOI: https://doi.org/10.1007/s10659-014-9508-z
5. Choi H.S., Vito R. Two-dimensional stress-strain relationship for canine pericardium // Journal of Biomechanical Engineering. - 1990. - Vol. 112(2). - P. 153-159. DOI: https://doi.org/10.1115/1.2891166
6. Chuong C., Fung Y. Three-dimensional stress distribution in arteries // Journal of Biomechanical Engineering. -1983. - Vol. 105(3). - P. 268-274. DOI: https://doi.org/ 10.1115/ 1.3138417
7. Dibb R., Qi Y., Liu C. Magnetic susceptibility anisotropy of myocardium imaged by cardiovascular magnetic resonance reflects the anisotropy of myocardial filament a-helix polypeptide bonds // Journal of Cardiovascular Magnetic Resonance. - 2015. - Vol. 17(1). - P. 1-14. DOI: https://doi.org/10.1186/s12968-015-0159-4
8. Ferruzzi J., Vorp D.A., Humphrey J. On constitutive descriptors of the biaxial mechanical behaviour of human abdominal aorta and aneurysms // Journal of the Royal Society Interface. - 2011. - Vol. 8(56). - P. 435-450. DOI: https://doi.org/10.1098/rsif.2010.0299
9. Fung Y.C. Biomechanics: mechanical properties of living tissues. - Springer Science & Business Media, 2013. -568 p.
10. Golob M., Moss R.L., Chesler N.C. Cardiac tissue structure, properties, and performance: a materials science perspective // Annals of Biomedical Engineering. - 2014.
характеристики материала могут быть использованы для более точного компьютерного моделирования механической функции сердца. Эти физические параметры могут быть использованы в будущем для разработки конечно-элементных моделей для понимания того, как инфаркт миокарда влияет на глобальное функционирование сердца. Точное моделирование механических свойств в области сердца может привести к улучшению хирургических методов лечения.
Продолжение исследований в этом направлении перспективно в медицине для оценки структуры, состояния и жизнеспособности тканей, а также определения интенсивности протекающих в них патофизиологических процессов.
- Vol. 42(10). - P. 2003-2013. DOI: https://doi.org/ 10.1007 /s10439-014-1071-z
11. Holzapfel G.A., Gasser T.C., Ogden R.W. A new constitutive framework for arterial wall mechanics and a comparative study of material models // Journal of Elasticity and the Physical Science of Solids. - 2000. -Vol. 61(1). - P. 1-48. DOI: https://doi.org/10.1023 /A:1010835316564
12. Holzapfel G.A., Sommer G., Gasser C.T., Regitnig P. Determination of layer-specific mechanical properties of human coronary arteries with nonatherosclerotic intimal thickening and related constitutive modeling // American Journal of Physiology-Heart and Circulatory Physiology.
- 2005. - Vol. 289(5). - P. H2048-H2058. DOI: https://doi.org /10.1152/ajpheart.00934.2004
13. Humphrey J.D. Continuum biomechanics of soft biological tissues // Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences. - 2003. - Vol. 459(2029). - P. 346. DOI: https://doi.org/ 10.1098/rspa.2002.1060
14. Hunter P.J., McCulloch A.D., Ter Keurs H.E.D.J. Modelling the mechanical properties of cardiac muscle // Progress in Biophysics and Molecular Biology. - 1998. -Vol. 69,(2-3). - P. 289-331. DOI: https://doi.org/10.1016 /S0079-6107(98)00013-3
15. Kakaletsis S., Meador W.D., Mathur M., Sugerman G.P., Jazwiec T., Malinowski M., Lejeune E., Timek T.A., Rausch M.K. Right ventricular myocardial mechanics: Multi-modal deformation, microstructure, modeling, and comparison to the left ventricle // Acta Biomaterialia. -2021. - Vol. 123. - P. 154-166. DOI: https://doi.org/ 10.1016/ j.actbio.2020.12.006
16. Laurence D., Ross C., Jett S., Johns C., Echols A., Baumwart R., Towner R., Liao J., Bajona P., Wu Y., Lee C.-H. An investigation of regional variations in the biaxial mechanical properties and stress relaxation behaviors of porcine atrioventricular heart valve leaflets // Journal of Biomechanics. - 2019. - Vol. 83. - P. 16-27. DOI: https://doi.org/10.1016/jjbiomech.2018.11.015
17. Li W. Biomechanics of infarcted left ventricle - A review of experiments // Journal of the Mechanical Behavior of Biomedical materials. - 2020. - Vol. 103. - P. 103591. DOI: https://doi.org/10.1016/j .jmbbm.2019.103591
18. Mas P.T., Rodríguez-Palomares J.F., Antunes M.J. Secondary tricuspid valve regurgitation: a forgotten entity
// Heart. - 2015. - Vol. 101(22). - P. 1840-1848. DOI: https://doi.org/10.1136/ heartjnl-2014-307252
19. Masithulela F. Analysis of passive filling with fibrotic myocardial infarction // ASME International Mechanical Engineering Congress and Exposition. Conference proceedings. - Houston: USA, 2015. DOI: https://doi.org/ 10.1115/IMECE2015-50003
20. Masithulela F. The effect of over-loaded right ventricle during passive filling in rat heart: A biventricular finite element model // ASME International Mechanical Engineering Congress and Exposition. Conference proceedings. - Houston: USA, 2015. DOI: https://doi.org /10.1115/ IMECE2015-50004
21. Masithulela F. Bi-ventricular finite element model of right ventricle overload in the healthy rat heart // Bio-medical Materials and Engineering. - 2016. - Vol. 27(5). - P. 507-525. DOI: https://doi.org/10.3233/BME-161604
22. Masithulela F.J. Computational biomechanics in the remodelling rat heart post myocardial infarction. PhD thesis. - South Africa: Cape Town: University of Cape Town, 2016. - 233 p.
23. Ndlovu Z., Nemavhola F., Desai D. Biaxial mechanical characterization and constitutive modelling of sheep sclera soft tissue // Russian Journal of Biomechanics. - 2020. -Vol. 24(1). - P. 84-96. DOI: https://doi.org/10.15593 /RJBiomech/ 2020.1.09
24. Nemavhola F. Biaxial quantification of passive porcine myocardium elastic properties by region // Engineering Solid Mechanics. - 2017. - Vol. 5(3). - P. 155-166. DOI: https://doi.org/10.5267/j.esm.2017.6.003
25. Nemavhola F. Fibrotic infarction on the LV free wall may alter the mechanics of healthy septal wall during passive filling // Biomedical Materials and Eengineering. - 2017. - Vol. 28(6). - P. 579-599. DOI: https://doi.org/10.3233 /BME-171698
26. Nemavhola F. Detailed structural assessment of healthy interventricular septum in the presence of remodeling infarct in the free wall - A finite element model // Heliyon. - 2019. - Vol. 5(6). - P. e01841. DOI: https://doi.org/10.1016/ j.heliyon.2019.e01841
27. Nemavhola F. Mechanics of the septal wall may be affected by the presence of fibrotic infarct in the free wall at end-systole // International Journal of Medical Engineering and Informatics. - 2019. - Vol. 11(3). - P. 205-225. DOI: https://doi.org/10.1504 /IJMEI.2019.101632
28. Nemavhola F. Study of biaxial mechanical properties of the passive pig heart: material characterisation and categorisation of regional differences // International Journal of Mechanical and Materials Engineering. - 2021.
- Vol. 16(1). - P. 1-14. DOI: https://doi.org/10.1186 /s40712-021-00128-4
29. Nemavhola F., Ngwangwa H., Davies N., Franz T. Passive biaxial tensile dataset of three main rat heart myocardia: left ventricle, mid-wall and right ventricle // Preprints. - 2021. - 2021080153. DOI: https://doi.org/ 10.20944/ preprints202108.0153.v1
30. Ngwangwa H.M., Nemavhola F. Evaluating computational performances of hyperelastic models on supraspinatus tendon uniaxial tensile test data // Journal of Computational Applied Mechanics. - 2021. - Vol. 52(1).
- P. 27-43. DOI: https://doi.org/10.22059 /jcamech.2020.310491.559
31. Ngwangwa H., Nemavhola F., Pandelani T., Msibi M., Mabuda I., Davies N., Franz T. Determination of cross-directional and cross-wall variations of passive biaxial mechanical properties of rat myocardium // Preprints. -2022. - 2021090244. DOI: https://doi.org/10.3390/ pr10040629
32. Ngwangwa H.M., Pandelani T., Nemavhola F. The application of standard nonlinear solid material models in modelling the tensile behaviour of the supraspinatus tendon // Preprints. - 2021. - 2021080298. DOI: https://doi.org/10.20944/preprints202108.0298.v1
33. Rigolin V.H., Robiolio P.A., Wilson J.S., Harrison J.K., Bashore T.M. The forgotten chamber: the importance of the right ventricle // Catheterization and Cardiovascular Diagnosis. - 1995. - Vol. 35(1). - P. 18-28. DOI: https://doi.org/10.1002/ccd.1810350105
34. Sacks M., Chuong C. Biaxial mechanical properties of passive right ventricular free wall myocardium // Journal of Biomechanical Engineering. - 1993. - Vol. 115(2). - P. 202-205. DOI: https://doi.org/10.1115/1.2894122
35. Sheehan F., Redington A. The right ventricle: anatomy, physiology and clinical imaging // Heart. - 2008. -Vol. 94(11). - P. 1510-1515. DOI: https://doi.org/10.1136 /hrt.2007.132779
36. Sirry M.S., Butler J.R., Patnaik S.S., Brazile B., Bertucci R., Claude A., McLaughlin R., Davies N.H. 4 , Liao J., Franz T. Characterisation of the mechanical properties of infarcted myocardium in the rat under biaxial tension and uniaxial compression // Journal of the Mechanical Behavior of Biomedical Materials. - 2016. - Vol. 63. -P. 252-264. DOI: https://doi.org/10.1016 /j.jmbbm.2016.06.029
Благодарности. Исследование было проведено при финансовой поддержке Южно-Африканского университета. Любое мнение, выводы и заключения или рекомендации, выраженные в этой публикации, принадлежат авторам, и поэтому ЮжноАфриканский университет не несет никакой ответственности в связи с этим.
Конфликт интересов. Авторы заявляют об отсутствии конфликта интересов.
Заявления об одобрении и согласии на участие. Этика исследования на животных была одобрена Комитетом по этике Колледжа науки техники и технологий Южно-Африканского университета 4 сентября 2021 года под контрольным номером 2019/СБЕТ 80Е/ЛКШ01.
FITTING OF HYPERELASTIC CONSTITUTIVE MODELS IN DIFFERENT SHEEP HEART REGIONS BASED ON BIAXIAL MECHANICAL TESTS
F. Nemavhola1, T. Pandelani12, H. Ngwangwa1
1 College of Science Engineering and Technology, University of South Africa, Pretoria, South Africa
2 The Council for Scientific and Industrial Research, Pretoria, South Africa
ABSTRACT
Heart failure remains one of the leading causes of death especially among people over the age of 60 years worldwide. To develop effective therapy and suitable replacement materials for the heart muscle it is necessary to understand its biomechanical behaviour under load. This paper investigates the passive mechanical response of the sheep myocardia excised from three different regions of the heart. Due to the relatively higher cost and huge ethical demands in acquisition and testing of real animal heart models, this paper evaluates the fitting performances of six different constitutive models on the myocardial tissue responses. Ten sheep were sacrificed, and their hearts were excised and transported within 3 hours to the testing biomechanical laboratory. The upper sections of the hearts above the short axes were carefully dissected out. Tissues were dissected from the mid-sections of the left ventricle, mid-wall and right ventricle for each heart. The epicardia and endocardia were then carefully sliced off each tissue to leave the myocardia. Stress-strain curves were calculated, filtered and resampled. The results show that Choi-Vito model was found to provide the best fit to the left ventricle, the polynomial (anisotropic) model to right ventricle, the Four-Fiber Family model to right ventricle, Holzapfel (2000) to right ventricle, Holzapfel (2005) to right ventricle and the Fung model to left ventricle.
©PNRPU
ARTICLE INFO
Received: 21 October 2021 Approved: 10 June 2022 Accepted for publication: 14 June 2022
Key words:
cardiac mechanics, experimental mechanics, sheep heart mechanics, hy-perelastic constitutive model, model fitting, soft tissue mechanics, biaxial testing.