Математические модели структуры сосудистого русла тела человека (обзор литературы)
А.В. Дмитриев, Я. А. Лысых, О.К. Зенин
Аннотация
Целью работы явился сравнительный анализ и систематизация имеющихся в литературе сведений о математическом моделировании структуры сосудистого русла человека. Был проведен анализ 51 издания среди публикаций за последние 10 лет.
Наряду с гидродинамическим моделированием сердечно-сосудистой системы существует и структурное моделирование, что исследует как строение различных тканей и органов, так и принцип построения сосудистого русла. Модели структуры тканей различаются в первую очередь точностью и производительностью, в то время как подходы к построению иерархической структуры сосудистого русла могут различаться на фундаментальном уровне. Описанные подходы структурного моделирования постоянно развиваются как путем улучшения старых моделей и разработкой новых, так и за счет увеличения доступных вычислительных мощностей. Результаты анализа свидетельствует о необходимости проведения дальнейших разработок в данной области и создания единой универсальной теории.
Ключевые слова: структурное моделирование, сердечно-сосудистая система, гемодинамика.
Abstract
Mathematical models of the vascular bed structure of the human body (literature review)
V. Dmitriev, Ya.A. Lysykh, O.K. Zenin
Background. The aim of the work is to compare and systematize information available in the literature on mathematical modeling of the structure of the human vascular bed. 51 literary sources of the last decade were analyzed.
Along with the hydrodynamic modeling of the cardiovascular system, there is also structural modeling, which explores both the structural structure of various t№ and organs and the principle of the construction of the vascular bed. Structure models differ primarily in the accuracy and performance associated with the presentation of materials, while approaches to constructing the hierarchical structure of the vascular bed may differ at a fundamental level. The described approaches of structural modeling are constantly evolving both by improving the old models and developing new ones, and by increasing the available computing power. The results of the analysis indicate the need for further development in this area and the creation of a unified universal theory.
Keywords: structural modeling, cardiovascular system, hemodynamics.
Введение
В связи с ростом количества сердечно-сосудистых заболеваний боль-шую перспективу приобрели исследования в области математического моде-лирования сердечно-сосудистой системы. Такие исследования позволят пред-сказывать развитие различных заболеваний (например атеросклероза), пере-нести эксперименты гп ъШсо при разработке лекарств, а также оценивать эф-фективность различных методов лечения, таких как протезирование клапанов сердца, ангиопластика и др. Сегодня можно выделить два направления моде-лирования системы кровотока - гемодинамическое и структурное. Первое рассматривает движение крови в сердечно-сосудистой системе и ее взаимо-действие с окружающими тканями как жидкости, второе же сосредоточено на рассмотрении принципов строения сосудистого русла в целом или его от-дельных участков. Эти два направления часто дополняют друг друга, однако могут быть самостоятельными. В данной статье рассмотрены методы и спо-собы структурного моделирования сердечно-сосудистой системы известные в настоящее время. Был проведен анализ 51 издания среди публикаций за по-следние 10 лет.
Структурное моделирование как моделирование строения различных частей сердечно-сосудистой системы можно разделить на направления - дискретное и сплошной среды [1]. Дискретное направление основано на об-работке каждой частицы отдельно, оно дает возможность моделировать твердое вещество как неоднородное для обнаружения дефектов, которые могут возникнуть лишь на ограниченной области. Очевидно, что такие воз-можности сопряжены с повышенными требованиями к вычислительным мощностям, а это приводит к использованию данного подхода только к не-большим вычислительным областям [2-4], в то время как подход сплошной среды применяется к обширным расчетным зонам [5-7]. Однако развитие цифровых технологий, увеличение их мощности и скорости обработки ин-формации постепенно убирает различие в целевых областях между данны-ми подходами.
Если говорить о предмете структурного моделирования, то в основном исследуются характеристики тканей органов и сосудов и их взаимодействие с кровью, однако существуют дискретные модели красных кровяных клеток для изучения влияния изменения их характеристик как естественных, так и паталогических на процессы в организме. Далее приведены примеры функ-ций энергии напряжения для некоторых моделей материалов, используемых в структурном моделировании и их описание.
Модель пео-Ноокеап (простая гиперэластичность)
ж=-2 (/1 - 3)-ц 1п 3 Д (1п з )2.
Эта модель и линейная эластичность дают похожие профили деформа-ции при небольших напряжениях. При больших напряжениях модель пеоНоокеап дает гораздо лучшее описание деформации и применима для напряжений вплоть до 20 % [8]. Также сердечно-сосудистые ткани считаются несжимаемыми (3 = 1), что упрощает функцию до одного слагаемого.
Анизотропная модель венечной артерии на основе экспоненци-альной функции Фанга
Данная модель реализована в [9] как расширение предыдущей модели многослойной артериальной стенки венечного сосуда сердца [10]. Модель включает в себя экспоненциальное изотропное слагаемое, предложенное Фангом [11], и анизотропное слагаемое, связанное с /4, что зависит от угла между направлением структуры волокон и периферическим направлением напряжения.
1. Модель красных кровяных клеток:
В статье [12] разработана 2D-модель сохраненного энергетического по-тенциала. Первое слагаемое определяет обычно малое напряжение деформа-ции, в то время как второе слагаемое представляет собой зависящее от изме-нения площади изотропное напряжение.
Достоинством модели, использующей приближение сплошной среды, является использование надежных уравнений механики сплошной среды, ре-шаемых не менее надежными численными методами. В первую очередь для структурного моделирования это метод конечных элементов (FE metod), поз-воляющий включать в расчет модель материала для учета особых характери-стик и поведения объекта при внешней нагрузке. Данный метод применяется в моделях как жидкостной [13, 14], так и твердой [15] части кровеносной си-стемы. Несмотря на удобство данного направления моделирования, его ис-пользование для представления локальных феноменов, таких как поврежде-ние или разрыв тканей, требует использования нестандартных уравнений. Хотя и предпринимаются попытки устранения этого недостатка, например расширенный метод конечных элементов (X-FEM) [16] или погруженный ме-тод конечных элементов [17], но дискретное моделирование продолжает быть более эффективным для локальных случаев, представляя естественное развитие таких процессов.
Дискретные методы моделирования опираются на набор точек, образующих 2D- или 3D-структуру, представляющую изучаемую поверхность или объем. Для биологических тканей данные структуры имеют треугольное расположение узлов, так как оно является наиболее репрезентативным [18]. Частицы внутри области соединяются связями, меняя соотношения нагрузки и деформации, для которых можно моделировать свойства материала объекта. Для воспроизведения сопротивления деформации могут быть созданы дополнительные связи между гранями, образованными тремя частицами. Также есть возможность указания силы контакта между несоединенными ча-стицами, что позволяет моделировать взаимодействия между раздельными твердыми телами, например такими, как створки сердечного клапана.
Благодаря моделированию отдельных частиц подобные модели могут учесть малые дефекты системы гораздо лучше, чем модели сплошной среды.
Однако из-за того, что даже малейшие объекты (например красные кровяные тельца) состоят из большого количества частиц, дискретные методы имеют высокие требования к вычислительным мощностям, что накладывает ограни-чение на их эффективное использование где-либо, кроме моделей самых ма-леньких объектов. Одна из основных моделей, используемая для построения иерархии структуры кровеносной системы, - модель бинарных деревьев. Данная модель используется в широком спектре задач - русла кровеносных и бронхиальных систем, корневых систем и крон растений, русла рек и других системах, обеспечивающих приток/отток веществ в распределенной в пространстве системе [19]. Подобное представление при определенных геометриче-ских соотношениях являет собой оптимальные трубопроводы, обеспечивающие минимальный расход энергии на транспортировку, что есть характеристика оптимальных систем [20, 21]. И так как все в природе стремится к оптимальности, данный алгоритм может достаточно точно представлять кровяное русло как различных внутренних органов, так и всего организма в целом. Данная модель иерархии структуры часто используется совместно с гидроди-намическим моделированием, позволяя рассматривать течение крови в дву-мерной или трехмерной задаче на реалистичной геометрии. Кровеносная си-стема как бинарное дерево рассмотрена в [22-37].
При построении геометрии в подобных моделях обычно используют установленные зависимости между диаметрами, углами разветвлений и дли-нами сегментов сосудов [22]. Для диаметров сегментов используется закон Мюррея:
где й?о - диаметр родительского сосуда; ^2 - диаметры дочерних сосудов; у~ 3 (у = 2,55 -- 3,02 и у = 2,76 -- 3,02 для артериальных и венозных русел, у = 2,62 -- 2,91 для бронхиального дерева) [38].
Для углов ветвления выполняются соотношения Розена [39]:
где ^ = тт {^, й2 } /тах {^, ^2 } - коэффициент асимметрии бифуркации.
Для длин сегментов артерий выполняется аллометрический принцип
где ае[2,6; 7,59], Ре [0,84; 1,16]. Для внутриорганных русел существуют аппроксимации зависимостей Ь^ (dj) полиномами и другими функциями, специфическими для разных органов [41].
Широко распространенной альтернативой модели бинарных деревьев является электрическая модель сердечно-сосудистой системы, основанная на аналогии между электрическими характеристиками и гидродинамическими свойствами сосудов. Данный подход обычно используется для рассмотрения всей системы кровообращения в целом или ее крупных частей, так как позво-ляет смоделировать, пусть и сильно упрощенно, взаимодействие как отдель-ных частей сердечно-сосудистой системы между собой, так и общую картину качества работы сердечно-сосудистой системы. Различные вариации элек-трической модели представлены в работах [42-51].
Для представления кровеносной системы как электрической ее харак-теристики и части переводятся в аналогичные элементы и характеристики электрических систем. Обычно каждая артерия представлена набором, состо-ящим из резистора (представляющим сопротивление), конденсатором (моде-лирующим расширяемость сосуда) и индуктором (реализующим инерцию крови). Сердце либо определяется как источник переменного тока, произво-дящий напряжение, которое представляет собой давление крови, либо также представляется в виде совокупности электрических элементов, что моделиру-ет работу предсердий и желудочков сердца. Характеристики системы крово-обращения, такие как давление, поток, объем, сопротивление и расширяе-мость, заменяются соответственно на напряжение, силу тока, заряд, сопро-тивление и емкость для электрической системы.
Заключение
Говоря о структурном моделировании сердечно-сосудистой системы, в зависимости от контекста имеют в виду модели либо строения и материалов тканей и органов, либо иерархической структуры сосудов. Для моделей тка-ней и органов существует много различных подходов, которые можно разде-лить на две группы - дискретные модели и модели сплошной среды. И если первая группа оперирует структурами, состоящими из отдельных частиц со своими свойствами, что дает высокую точность при расчете различных со-стояний, но и накладывает высокие требования к вычислительным мощно-стям, то вторая группа использует осредненные характеристики материалов, имея меньшую точность, чем дискретный подход, но при этом использует го-раздо меньшие вычислительные мощности. Если говорить об иерархическом моделировании, то тут также существует два основных подхода - модель би-нарных деревьев и электрическая модель. Первый подход представляет собой построение сосудистой структуры при помощи установленных математиче-ских зависимостей в виде разветвляющихся деревьев и может использоваться как самостоятельно, например для изучения оптимальности строения крове-носной системы, так и в сочетании с гидродинамическим моделированием. Электрическая же модель сочетает в себе как иерархическую структуру сер-дечно-сосудистой системы в виде электрической цепи, так и расчетную мо-дель различных параметров кровотока с использованием их электрических аналогов. Описанные подходы структурного моделирования постоянно раз-виваются как улучшением старых моделей и разработкой новых, так и за счет увеличения доступных вычислительных мощностей.
моделирование венечный артерия
Библиографический список
1. Structural modelling of the cardiovascular system / B. Owen [et al.] // Biomechanics and Modeling in Mechanobiology. - 2018. - Vol. 17, № 5. - P. 1217-1242.
2. Fedosov, D. A. A multiscale red blood cell model with accurate mechanics, rheology, and dynamics / D. A. Fedosov, B. Caswell, G. E. Karniadakis // Biophysical Journal. - 2010. - Vol. 98, № 10. - P. 2215-2225.
3. Spectrin-level modeling of the cytoskeleton and optical tweezers stretching of the eryth¬rocyte / J. Li [et al.] // Biophysical Journal. - 2005. - Vol. 88. - P. 3707-3719.
4. Nakamura, M. Spring-network-based model of a red blood cell for simulatiing mesoscopic blood flow / M. Nakamura, S. Bessho, S. Wada // International Journal for Numerical Methods in Biomedical Engineering. - 2013. - Vol. 29, № 1. - P. 114-128.
5. Blunt trauma and acute aortic syndrome: a three-layer finite-element model of the aortic wall / A. R. Zhao [et al.] // European Journal of Cardio-Thoracic Surgery. - 2008. - V. 34, № 3. - P. 623-629.
6. A coupled momentum method for modeling blood flow in three-dimensional deforma¬ble arteries / C. A. Figueroa [et al.] // Computer Methods in Applied Mechanics and Engineering. - 2006. - Vol. 195, № 41-43. - P. 5685-5706.
7. Fluidstructure interaction simulation of aortic blood flow / P. Crosetto [et al.] // Com¬puters & Fluids. - 2011. - Vol. 43, № 1. - P. 46-57.
8. Gent, A. N. Engineering with rubber: how to design rubber components / A. N. Gent. - 3rd edition. - Munich : Hanser Publishers, 2012. - 451 p.
9. Determination of layer-specific mechanical properties of human coronary arteries with nonatherosclerotic intimal thickening and related constitutive modeling / A. Gerhard [et al.] // American Journal of Physiology: Hearth and Circulatory Physiology. - 2005. - Vol. 289, № 5. - P. 2048-2058.
10. Gerhard, A. H. A New Constitutive Framework for Arterial Wall Mechanics and a Comparative Study of Material Models / A. H. Gerhard, C. G. Thomas, W. O. Ray // Journal of elasticity and the physical science of solids. - 2000. - Vol. 61, № 1-3. -
11. Fung, Y. C. Elasticity of soft tissue in simple elongation / Y. C. Fung // American Journal of Physiology. - 1967. - Vol. 213, № 6. - P. 1532-1544.
12. Strain energy function of red blood cell membranes / R. Skalak [et al.] // Biophysical Journal. - 1973. - Vol. 13, № 3. - P. 245-264.
13. Taylor, C. A. Finite element modeling of three-dimensional pulsatile flow in the abdominal aorta: relevance to atherosclerosis / C. A. Taylor, T. J. Hughes, C. K. Zarins // Annals of Biomedical Engineering. - 1998. - Vol. 26, № 6. - P. 975-987.
14. Oshima, M. Finite element simulation of blood flow in the cerebral artery / M. Oshima, R. Torii, T. Kobayashi // Computer Methods in Applied Mechanics and Engineering. - 2001. - Vol. 191, № 6-7. - P. 661-671.
15. Kim, H. J. Incorporating autoregulatory mechanisms of the cardiovascular system in three-dimensional finite element models of arterial blood flow / H. J. Kim, K. E. Jansen, C. A. Taylor // Annals of Biomedical Engineering. - 2010. - Vol. 38, № 7. - P. 2314¬2330.
16. Abdelaziz,Y. A survey of the extended finite element / Y. Abdelaziz, A. Hamouine // Computers & Structures. - 2008. - Vol 86, № 11-12. - P. 1141-1151.
17. Immersed finite element method and its applications to biological systems / W. K. Liu [et al.] // Computer Methods in Applied Mechanics and Engineering. - 2006. - Vol. 195, № 13-16. - P. 1722-1749.
18. Abraham, F. F. Diffraction from polymerized membranes: flat vs crumpled /
19. Bejan, A. Shape and Structure: from Engineering to Nature / A. Bejan. - Cambridge : Univ. Press, 2000. - 210 p.
20. Черноусько, Ф. Л. Оптимальная структура ветвящихся трубопроводов / Ф. Л. Черноусько // Прикладная математика и механика. - 1977. - T. 41, № 2. - C. 376-383.
21. La Barbera, M. Principles of design of fluid transport systems in zoology / M. La Barbera // Science. - 1990. - Vol. 249. - P. 992-1000.
22. Балабанов, В. А. Математическое моделирование и 3Б-визуализация бинар¬ных деревьев с минимумом самопересечений / В. А. Балабанов, Н. Н. Кизилова // Вісник Харківського національного університету імені В. Н. Каразіна. - 2017. -