The feces flow in colon: Biomechanics, rheology, simulation

Abstract


Modeling can allow us to realize the functionality of a living organ and analyze the influence of various factors on it. In particular, modeling the movement of fecal masses through the intestine allows us to investigate the relationship between motility and the processing of intestinal contents, the time of its progress, and the role played by peristaltic waves. The intestine is a complicated system of complex geometric shape, wave movements of several configurations, and multi-level processing of internal contents. At the same time, a medical examination of the individual features of its structure and function is quite time-consuming. Building a personalized model that takes into account all these factors is a non-trivial task. However, such a model (if it can be developed) will minimize the use of invasive assessments of a patient's condition. It can also be used as the basis for planning therapeutic measures for a particular patient. This article presents an overview of current mathematical models developed in the past couple of decades for simulating the movement of fecal masses through the large intestine.

Full Text

Более 40 % населения земного шара страдают функциональными желудочно-кишечными расстройствами [78]. Согласно данным [2], 2017 г. во всем мире было зарегистрировано 6,8 млн случаев воспалительных заболеваний кишечника (рис. 1). В то же время функциональные нарушения кишечника являются причиной постоянного дискомфорта и боли, которые приводят к снижению качества жизни человека. Согласно исследованию [6], существует определенная связь между синдромом раздраженного кишечника, курением и употреблением алкоголя. При оценке смертности от дивертикулярной болезни [26] было установлено, что более 40 % приходится на развитые страны. При исследовании хронического колостаза у детей [55] были установлены структурные изменения стенки толстой кишки, ее тонуса и характера перистальтики. Ученые выявили характерные признаки дегенеративных расстройств, таких как: неравномерное утолщение стенки толстой кишки (более чем на 5 мм), потеря ее слоистой структуры и неравномерная хаустрация, уплощение дна хаустры. Другим распространенным кишечным заболеванием является опухоль. По разным данным, доброкачественные опухоли кишечника выявляются у 16 - 40 % населения. Во многих случаях доброкачественный процесс со временем превращается в злокачественный. Во всем мире колоректальный рак является третьей, наиболее часто диагностируемой, злокачественной опухолью и второй по значимости причиной смерти от рака [31]. Таким образом, можно утверждать, что заболевания толстой кишки являются важным объектом внимания исследователей. Медицинское сообщество сигнализирует о росте распространенности кишечных заболеваний, в то время как полное излечение от воспалительных процессов кишечника маловероятно и приводит к хроническому течению или к обострению. Наблюдение и лечение кишечных заболеваний является острой проблемой как для пациентов (боль и дискомфорт), так и для врачей (длительное лечение и диагностика осложнений). Поэтому изучение кишечника в норме и патологических состояниях, выявление причинно-следственных связей нарушений его здорового функционирования становится все более актуальным. В современном мире изучение физиологических процессов живых организмов в нормальных и патологических состояниях все чаще описывается с использованием экспериментальных или математических моделей. Такой подход позволяет нам в некоторой степени объективизировать эти процессы и использовать их в прогностических целях. Целью этого исследования является анализ основных исследований толстой кишки, чтобы описать ее механическое функционирование и выделить основные факторы, влияющие на ее двигательную функцию и возникновение воспалительных процессов. Результаты анализа будут полезны для Рис. 1. Глобальный, региональный и национальный показатель распространенности воспалительных заболеваний кишечника (ВЗК) в 195 странах и территориях, 1990 - 2017 г.: a - стандартизированный по возрасту показатель распространенности (на 100 000 населения) ВЗК у обоих полов для 195 стран и территорий, 2017 г.: б - распространенность в зависимости от возраста с учетом половой принадлежности в 2017 г. от общего числа распространенных случаев и показатель распространенности ВЗК в зависимости от возраста на глобальном уровне; в - тенденции с 1990 по 2017 г в стандартизированном по возрасту уровне смертности от ВЗК в семи регионах [2] биомеха а нического моделирования толстой кишки и в б дальнейшего использования при оценке параметров кишечника. Моделирование нижних отделов пищеварительного тракта Моделирование процессов в пищеварительном тракте осуществлялось на разных уровнях, в том числе и в толстом кишечнике [34-39; 71; 77]. Создание модели для представления процесса пищеварения в толстой кишке является одной из важных проблем биомеханики пищеварительного тракта. Заключительной стадией пищеварения является перемещение фекалий по толстой кишке и их выведение через прямую кишку. Моделирование процесса пищеварения в толстой кишке позволит понять и описать здоровую физиологию этого процесса. Кроме того, модель может быть использована для описания патологических процессов (запор, диарея, опухоль), что может быть полезно для их диагностики и лечения. Моделирование всасывания веществ в кишечнике также может быть применено к фармацевтическим препаратам при медикаментозном лечении. Основными функциями толстой кишки являются: всасывание воды, минералов, двигательная функция (образование и выведение кала за счет перистальтики) и остаточное расщепление веществ за счет микрофлоры и кишечного сока тонкой кишки. Пищеварительный тракт - это сложная биомеханическая система. Смоделировать процессы, происходящие в нем, - непростая задача. Существует ряд исследований в этом направлении на разных уровнях. Обзор описания биомеханического, функционального и сенсорного моделирования желудочно-кишечного тракта представлен в статье [43]. Много исследований посвящено движению химуса по тонкому кишечнику, есть исследования, посвященные перевариванию в двенадцатиперстной кишке. В статье [12] был представлен обзор моделирования тонкой кишки in silico с учетом моторики, иннервации и патологических состояний. В обзоре [9] представлены параметры физиологически значимых моделей кишечника in vitro с акцентом на важность механоструктурных сигналов, таких как: топография, подвижность и поток, присутствующие в кишечном эпителии, которыми часто пренебрегают в классических системах культивирования клеток. Моделирование тонкого кишечника Рис. 2. Анатомия толстого кишечника [93] В работе [22] было представлено CFD-моделирование двенадцатиперстной кишки человека. Авторы рассматривали перистальтические движения стенки двенадцатиперстной кишки вместе с потоком химуса. Движение химуса описывалось уравнениями Навье - Стокса, а сокращения двенадцатиперстной кишки аппроксимировались сглаженной ступенчатой функцией Хевисайда. С помощью уравнения конвекции - диффузии описывается гидролиз крахмала с образованием глюкозы для имитации потока глюкозы через стенку двенадцатиперстной кишки. Лентл и др. [41] провели обзор перемешивания и потока химуса в тонком кишечнике. В работе говорится о том, что твердые субстраты и растворимые питательные вещества распределены по всей длине просвета кишечника, чтобы обеспечить оптимальное перемешивание и ферментативное переваривание. По мнению авторов, перемешивание в просвете происходит из-за локализованного складывания относительно неэластичной слизистой оболочки и перемешивания содержимого в результате повторяющихся радиальных и продольных мышечных сокращений, и усиливается нестабильным пространственным расположением сокращений. Математическая модель перистальтики в тонком кишечнике рассмотрена в статье [42]. Исследователи изучали течение ньютоновской жидкости с малым числом Рейнольдса в круглой цилиндрической трубе с серией перемещающихся сужений. Общее решение было получено линейной комбинацией двух фундаментальных решений: (1) перистальтический перенос с перемещением жидкости (при котором движение жидкости создается перистальтическими волнами без градиента давления); (2) перистальтическое сжатие без переноса жидкости. Моксоном и соавторами была разработана модель in silico в тонком кишечнике [52] для исследования переваривания и всасывания крахмала и глюкозы в тонком кишечнике. Авторы попытались определить влияние двух факторов: времени опорожнения желудка и вязкости слизи слизистой оболочки кишечника на скорость всасывания глюкозы. В статье использованы кинетические уравнения изменения концентрации веществ и массопереноса. При низкой вязкости слизи (менее 0,1 Па·с) скорость всасывания регулируется временем опорожнения желудка. При вязкости выше 0,1 Па·с 10-кратное увеличение вязкости может привести к 4-кратному снижению абсорбции глюкозы. Позже те же авторы [51] представили математическую модель механизма обратной связи между всасыванием питательных веществ и контролем опорожнения желудка. В статье [63] рассмотрена математическая модель перистальтического течения химуса в тонкой кишке вместе с помещенным эндоскопом. Химус аппроксимируется как жидкость Уильямсона. Поток рассматривается в кольцевой области, образованной между двумя концентрическими трубками, причем внешняя трубка представляет собой тонкую кишку, а внутренняя трубка представляет собой эндоскоп. Поток моделировался двумя синусоидальными перистальтическими волнами разной длины, проходящими вдоль стенки кишечника с одинаковой скоростью. Авторами получены выражения для осевой скорости, силы трения, давления и расхода. В статье [63] была создана модель для отдельной клетки гладкой мускулатуры тощей кишки человека на основе мембранного потенциала. Авторы получили описание клеточного поведения, которое способно воспроизводить экспериментально зарегистрированную медленноволновую активность. Экспериментальное исследование [64] с использованием МРТ метки было проведено для визуализации и оценки движения содержимого восходящей ободочной кишки. Авторы показали различия такого движения у здоровых испытуемых и испытуемых с запорами. В статье [73] представлена численная модель перистальтики двенадцатиперстной кишки с использованием суспензии твердых частиц в вязкой ньютоновской жидкости. Стенка кишечника моделировалась вязкоупругой мембраной, а движение стенки было в виде перистальтической волны. В исследовании моделировалось пропульсивное движение. Авторы пришли к выводу о ламинарном характере потока и неоднородном распределение твердых частиц по длине двенадцатиперстной кишки. Исследователями Таракан и др. [83] была разработана экспериментальная установка, моделирующая д) в) движение сегментации в тонкой кишке. Целью исследования было изучение явлений массопереноса, происходящих в просвете, и их потенциального влияния на концентрацию веществ, доступных для всасывания. При моделировании сегментационного движения фи а цент массопереноса увеличился в 7 раз. Авторы установили, что увеличение вязкости кишечного содержимого существенно влияло на коэффициент б) массопереноса и на содержание глюкозы, образующейся в результате гидролиза крахмала. д г е в б Рис. 3. Численное моделирование толстого кишечника: а - развитие перистальтической волны [73]; б - распределение максимальных главных напряжений [60]; в - влияние рациона питания на скорость транспортировки и на структуру каловых масс [40]; г - распределение слизи в толстой кишке [40]; д - распределение напряжений в момент максимального сжатия зажима [19]; е - распределение скорости с перистальтической волной в ректальном отделе через 0,5; 2,0 и 6,0 с [80] В исследовании [59] применялось вычислительное гидродинамическое моделирование потока и перемешивания кишечного содержимого в двенадцатиперстной кишке человека с анатомически реалистичной геометрией и паттернами сокращения. Исследовались особенности потока, такие как точки застоя и обратный поток в зависимости от реологии дигеста. Авторы обнаружили, что увеличение амплитуды сокращения приводит к увеличению скорости продвижения пищи по кишечнику и повышению степени перемешивания. Ряд работ [14; 23; 27; 28; 33; 57; 58] был выполнен исследователями на основе численного моделирования гидродинамики с использованием анатомически реалистичной геометрии желудочно-кишечного тракта. Результаты послужили основой для количественной оценки взаимосвязи между перистальтикой желудка, перемешиванием внутри и опорожнением желудка, уделяя особое внимание нарушению координации закрытия привратника с антральными сокращениями. Моделирование процесса пищеварения активно ведется для всех отделов пищеварительного тракта. Большое количество исследований направлено на исследование тонкой кишки и ее отделов (тощей кишки, двенадцатиперстной кишки, подвздошной кишки). Для построения модели толстой кишки необходим ряд параметров: волновые характеристики, параметры кала и стенки кишечника, а также граничные условия. Ученые пытаются найти ответы на все эти проблемы в своих исследованиях. Моделирование толстого кишечника В отличие от тонкого кишечника, моделированию толстого кишечника уделялось мало внимания (рис. 2). Возможно, сложное геометрическое строение (наличие гаустр) или недостаточно изученная сократительная способность (пропульсивные, медленные волны) являются причинами этого. Однако некоторые математические подходы, реализуемые для тонкого кишечника, могут быть использованы и для толстого кишечника. В частности, авторами в работе [54] рассмотрена теоретическая модель перистальтического потока вязкой жидкости через изогнутый канал под воздействием магнитного поля (число Рейнольдса предполагались малыми величинами, а длина волны - большой). Результаты исследования - это распределения давления, КПД, эффективность потока. Авторы сделали вывод, что величина давления возрастает с увеличением кривизны канала и магнитного параметра. Исследователи установили, что рефлюкс в изогнутом канале выше, чем в прямом. По мнению авторов, эффективность потока ниже в изогнутом, чем в прямом канале, а также при высоком значении магнитного параметра. В исследовании каловые массы моделировались как магнитная вязкая несжимаемая жидкость. Подобные модели магнитной гидродинамики были рассмотрены и в других исследованиях [1; 4; 24; 25; 61; 67; 69; 81]. В исследовании [84] применялся принцип Бернулли для описания движения содержимого толстого кишечника при запоре и диарее. Исследователи анализировали изменение времени прохождения толстой кишки и высокоамплитудные перистальтические сокращения дефекации у здоровых людей. Такой подход был реализован в рамках модели здорового человека, модели запора и модели диареи. Для каждой модели авторы вычисляли скорость течения каловых масс на основе объемной скорости притока из слепой кишки и данных о длине каждой части толстой кишки. Исследователи установили, что время прохождения толстой кишки обратно пропорционально объемной скорости притока из слепой кишки и зависит от размеров кишечного тракта. Йео и соавт. [90] была разработана модель описания поведения клетки гладкомышечной ткани толстой кишки человека. В модели авторы рассмотрели физиологические механизмы во время базальной электрической активности, с описанием основных ионных каналов и внутриклеточного гомеостаза кальция. При построении модели использовался подход Ходжкина - Хаксли для клеточной мембраны. Для моделирования ионных каналов были использованы экспериментальные данные для изолированных миоцитов толстого кишечника человека. В работе [79] была экспериментально исследована физическая модель дефекационной системы для количественного описания влияния аноректального угла и давления на сфинктер при воздержании. Моделирование течения каловых масс (вязкая ньютоновская жидкость) в гибкой вязкоупругой цилиндрической трубке, сокращающейся в радиальном направлении, осуществлялось исследователями в работе [72]. Исследовалось внутрипросветное давление, скорость течения, вихри, время течения. Похожие исследования проводились в работе Алексиадисом и др. [3] и Шуттом и др. [70]. Авторы сравнили три разные модели: полностью заполненную кишку, частично заполненную кишку и частично заполненную кишку с присутствием газовой фазы (газожидкостная модель). В исследованиях [20; 21] используется конечно-элементное проектирование для моделирования взаимодействия между зажимом и гиперэластичным материалом толстого кишечника. Результаты нужны для оценки контактной силы во время фиксации зажимом. Анализ проводился на двух геометриях (цилиндрическая область и реальная геометрия) с использованием гиперупругой модели Муни - Ривлина. На толстую кишку оказывается давление, имитирующее хирургический зажим, который постепенно соприкасается с поверхностью ткани. В статье [60] рассматривалась модель дивертикула толстой кишки в программе Abaqus. Расчетная площадь представляла собой осесимметричную цилиндрическую область со сферическим включением. Результаты расчетов показали, что высокие напряжения концентрируются вокруг шейки мешочка, и их значения и распространение увеличиваются с размером шейки мешочка, а не с площадью поверхности мешка. Моделирование толстого кишечника может осуществляться для разных целей [16; 21; 40; 45; 60; 72] (рис. 3): описание движения каловых масс, всасывание жидкости и питательных веществ, химические реакции (табл. 1). Для построения моделей требуется учитывать целый ряд факторов. В частности, для построения модели движения нужны граничные условия, значения давления на входе и выходе из толстого кишечника, учет анизотропии и проницаемости стенки, параметры волновых сокращений, реологические свойства каловых масс и слоя слизи. Таблица 1 Модели толстого кишечника Модель Автор Магнитная гидродинамика Narla et al., 2017[54], Tanveer et al., 2018 [81], Pavan et al., 2012 [61], Hina et al., 2012 [24], Sankad et al., 2016 [69], Ramanamurthy et al., 2013 [67], Abdulhadi et al., 2018 [1], Asghar et al., 2019 [4], Tanveer et al., 2017 [81], Hina et al., 2015 [25] Принцип Бернулли Yoshiharu Uno et al., 2018 [84] Модель, описывающая поведение клетки гладкой мускулатуры кишечника человека Yeoh et al., 2017 [90] Экспериментальная модель дефекации Stokes et al., 2018 [79] Моделирование взаимодействия зажима и гиперпластичного материала толстого кишечника Guachi et al., 2019 [20], Guachi et al., 2020 [19] Модель дивертикула толстого кишечника Patel et al., 2020 [60] Двухслойные конечно-элементные модели толстого кишечника Zhao et al., 2022 [92] Модель течения вязкой жидкости по гибкому вязкоупругому каналу Sinnott et al., 2012 [72], Schutt et al., 2020 [70], Alexiadis et al., 2017 [3] Интерес к физиологическим процессам в толстом кишечнике в последние годы большой. Применение математического аппарата реализуется на разных уровнях (клетка, поток). Однако целый спектр вопросов остаются открытыми для исследования. Реология каловых масс Важным параметром в моделировании движения пищи по пищеварительному тракту является реологическое поведение химуса, каловых масс, и влияние этого параметра на пищеварение. В соответствии с экспериментальными исследованиями [87; 88] (рис. 4), реологическое поведение каловых масс можно описать степенным законом, в котором учитывается содержание влаги: (1) где τ - кажущаяся вязкость, γ - скорость сдвига, S - содержание влаги. Кажущаяся вязкость человеческих фекалий варьировалась от 27 до 2014 Па·с для диапазона содержания влаги от 58,5 до 88,7 % при скорости сдвига 1 с-1. Исследователи обнаружили, что человеческие фекалии являются тиксотропными. В частности, их вязкость постоянно снижается до 0,5 % от исходного значения после 20-секундного воздействия скорости сдвига 10 с-1. Авторы получили линейную зависимость между вязкостью и температурой, при этом было зарегистрировано снижение вязкости на 30,6 % при повышении температуры до 35,6 °C с 13,4°C. Лубенс и др. [44] рассматривали каловые массы как пластичный материал. Авторы определяли предел текучести ряда фекалий различной консистенции и оценивали его влияние на патофизиологию дефекации. Предел текучести фекалий возрастал экспоненциально с увеличением содержания в них твердого вещества от 20 до 8000 Па. Исследователи показали, что диаметр ректоанального соединения увеличивается, чтобы обеспечить прохождение кала при высоких значениях предела текучести фекалий. Янг и др. [21] исследовали дефекацию млекопитающих от кошек до слонов, используя размеры толстого кишечника и кала, реологические измерений кала и слизи и математическую модель дефекации. Исследователи утверждают, что не только прямая, но и толстая кишка является резервуаром для фекалий. Кроме того, время испражнения млекопитающих составляет 12 с (N = 23). Авторы объясняют этот эффект с помощью математической модели, которая показывает, что фекалии скользят по толстой кишке по слою слизи. У более крупных животных не только больше кала, но и толще слой слизи, что облегчает их выход. Авторы описывают реологическое поведение каловых масс с помощью степенного закона (2) где η - кажущаяся вязкость фекалий, Kf - индекс консистенции потока, а nf = 0,21 - степенной показатель. Эти параметры исследователи определили следующим образом: 1) для каловых масс: Кf - составляет от 20 до 3000, nf = 0,21; 2) для слизи: Кs - составляет от 5.06, ns = 0,25. В обзоре [88] выделены основные механические параметры каловых масс и их состав. В частности, авторы ссылаются на величину вязкости, равную 3,5 - 5,5 Па·с, плотность - 1,06 г/мл3. В модели тонкого кишечника [53] применялась модель жидкости Вильямса. В работе [11] каловые массы описывают реологическим законом жидкости Джеффри, в статье [17] применялся закон для жидкости Сиско. При использовании подхода магнитной электродинамики в моделировании толстого кишечника жидкость принимается как жидкость Каро [82], вязкая жидкость [1; 3; 4; 24; 25; 54; 67; 69; 70; 72; 79; 84; 90] жидкость Джеффри [1], жидкость Сиско [81] [4], псевдопластическая жидкость [25], жидкость, реология которой описывается степенным законом [61]. Обобщение свойств текучести для каловых масс, рассмотренных в работах ученых, приведено в табл. 2. а г в б Рис. 4. Реологические свойства каловых масс: а - зависимость напряжения - скорости сдвига [88]; б - вязкость при различном содержании воды [87]; в - вязкость каловых масс и слизи [89]; г - предел текучести [44] Таблица 2 Реологические модели каловых масс Реологический закон Автор Степенной закон Yang et al., 2017 [89] Woolley et al., 2014 [88], Pavan et al., 2012 [61] Линейно-вязкая жидкость Hina et al., 2012 [24], Sankad et al., 2016 [69], Narla et al., 2017 [54], Ramanamurthy et al., 2013 [67] Жидкость Вильямса Nadeem et al., 2012 [53], Hina et al., 2015 [25] Жидкость Джеффри Disu et al., 2017 [11], Abdulhadi et al., 2018 [1] Жидкость Сиско German et al., 2018 [17], Asghar et al., 2019 [4], Tanveer et al., 2017 [81] Жидкость Каро Tanveer et al., 2018 [82] Механические свойства стенок толстого кишечника При моделировании пищеварения необходимы параметры, характеризующие механическое поведение стенки толстого кишечника под нагрузкой. Такие экспериментальные исследования проводились для тканей животных и человека. В работе [18] исследовались свойства и микроструктура стенки толстого кишечника при экспериментальном колите у мышей. Исследователи определили предельные напряжения при растяжении 232 ± 33 и 183 ± 25 кПа и конечные растяжения при разрыве 1,43 ± 0,04 и 1,51 ± 0,06 для контрольной группы и группы с колитом соответственно. В результате гистологического анализа авторы установили нарушение микроструктуры стенки толстой кишки из-за воспаления и отека подслизистой оболочки. Кроме того, большое количество коллагеновых волокон в мышечном слое было разрушено, а пучки волокон были тоньше при колите по сравнению с контрольной группой. Для анализа поведения колоректальной ткани под нагрузкой исследователи в работе [19] применяли конечно-элементную модель. Поведение материала описывалось тремя механическими моделями: 1) линейно-упругая модель; 2) гиперупругая модель Муни - Ривлина и 3) модель Йоха. Авторы пришли к выводу, что гиперупругие модели являются более точными по сравнению с линейно-упругими для прогнозирования полей напряжений и деформаций. В работе [47] были проведены механические экспериментальные исследования на тканях толстого кишечника мышей, расположенных в разных ее местах (восходящая, поперечная, нисходящая и сигмовидная кишка). Испытания проводились в продольном и окружном направлениях со скоростью 1 м/с. Механический отклик на динамическое одноосное нагружение показал двухслойный и двухфазный процесс, зависящий от направления воздействия и локализации. Подобное исследование этой же группой ученых [49] было проведено для толстого кишечника человека in vitro в экспериментах на растяжение при различных скоростях нагружения. Образцы были подвергнуты одноосным испытаниям на растяжение под воздействием динамической нагрузки (1 м/с), промежуточной нагрузки (10 см/с) и квазистатической нагрузки (1 см/с). Авторы определили, что в окружном направлении образцы выдерживают более высокие уровни нагрузки до повреждения, чем в продольном направлении. Исследователи также обнаружили, что присутствие taenia coli влияет на механическую реакцию при низкоскоростной нагрузке. Авторы расширили свое исследование в работе [48], оценив влияние пола и возраста с использованием трех различных протоколов: 1) охлажденные трупы без бальзамирования, 2) забальзамированные трупы и 3) свежая ткань толстой кишки на механические свойства толстой кишки. В первом приближении данные, полученные от животных, могут быть использованы для моделирования толстой кишки, но свойства тканей кишечника мыши и человека отличаются друг от друга. Исследователи [65] провели ряд экспериментов на тканях толстого кишечника свиней. Авторы обнаружили, что в зависимости от места изъятия образца вдоль кишечника по мере приближения к прямой кишке свойства материала менялись с более изотропного на более анизотропный материал, а податливость в окружном направлении увеличивалась. В работе [68] были экспериментально исследованы семь органов брюшной полости свиней, в том числе толстый кишечник, для оценки механических характеристик тканей. Для описания нелинейного поведения материала использовались экспоненциальные законы и закон Блац - Ко. Авторы в работе [74] экспериментально исследовали механические свойства толстого кишечника у мышей. Полученные результаты свидетельствуют, что эластичность тканей толстого кишечника постепенно увеличивается по направлению к прямой кишке. Исследователи установили, что материал прямого кишечника обладает более вязкоупругими свойствами в окружном направлении, чем в осевом. Также была выявлена богатая коллагеновая сеть в подслизистой оболочке с ориентацией волокон примерно на 30° по направлению к осевому направлению. Этот факт подтверждает результаты двухосного растяжения (где жесткость почти вдвое больше в осевом направлении по сравнению с окружным направлением). В исследовании [45] была проведена количественная оценка толщины каждого отдельного слоя прямой кишки, а также основные ориентации и распределения коллагеновых волокон внутри каждого из этих слоев и их диаметр. В работе [10] в модели магнитной гидродинамики течения каловых масс по кишечнику стенка кишечника моделировалась, как вязкоупругая ткань. Экспериментальные исследования, проведенные для исследования механических свойств стенки толстого кишечника, приведены в табл. 3. Таким образом, можно сделать вывод, что толстый кишечник является анизотропным (более эластичным в радиальном направлении) гиперупругим материалом, для описания поведения которого могут быть использованы разные гиперупругие модели. Волновые характеристики потока каловых масс Движение каловых масс по толстому кишечнику осуществляется под действием волновых движений. В большинстве исследований волны моделируются как синусоидальные. В частности, в работе [86] в модели магнитной гидродинамики была рассмотрена синусодальная волна и предложены треугольная и квадратная волна. Таблица 3 Свойствах стенки толстого кишечника. Экспериментальные данные Экспериментальные данные Автор Образцы Предельные растягивающие напряжения составляли 232 кПа Предельные растяжения при разрыве 1,43 Gong et al., 2017 [18] Мышь Окружное напряжение при разрыве - 1,02 МПа, а продольное - 0,70 МПа. Окружная деформация при разрыве - 69 %, а продольная - 55 % МПа Massalou et al., 2016 [47] Мышь Для продольных образцов модуль Юнга - 3,17 МПа при динамическом нагружении (1 м/с), 1,74 МПа при промежуточном нагружении (10 см/с) и 1,76 МПа при квазистатическом нагружении (1 см/с). Для окружного образца модуль Юнга - 3,15 МПа при динамическом нагружении (1 м/с), 2,14 МПа при промежуточном нагружении (10 см/с) и 0,63 МПа при квазистатическом нагружении (1 см/с) Massalou et al., 2019 [49] Человек Модуль Юнга для различных протоколов хранения: бальзамированный - 3,08; свежий - 2,97; и охлажденный - 3,17 Massalou et al., 2021 [48] Человек Двухосные тесты проводились на различных участках толстой кишки. Были определены параметры гиперупругих моделей (модель Хольцапфеля - Гассера - Огдена, модель Гассера - Огдена - Хольцапфеля, модель для четырех волокон, микроструктурные модели, модель микроволокон Фон Мизеса, модель микроволокон Бингама) Puértolas et al., 2020 [65] Свинья Получены кривые напряжение - деформация (in vivo). Эластичная модель: Модель стресс-релаксации: Rosen et al., 2008 [68] Свинья Получены средние осевые (5,010) и окружные (9,041) деформации образца и графики зависимости напряжений Коши от растяжения при двухосном растяжении Siri et al., 2019, [74; 75] Свинья Для восходящей ободочной кишки, поперечной ободочной кишки, нисходящей ободочной кишки, сигмовидной кишки определены прочность на разрыв (г), площадь поперечного сечения (мм2), предел прочности при растяжении (г/мм2), внутренний диаметр (мм), ширина при разрыве (мм), удлинение, релаксация напряжения Watters et al., 1985 [85] Человек Получены параметры напряжения - деформации толстой кишки: для поперечных образцов: уровень деформации, при котором начинается резкое увеличение напряжения 33,29 %, уровень деформации разрушения - 176,66 %, первое максимальное напряжение - 0,645 МПа, уровень деформации, соответствующий первому максимальному напряжению 87,85 %, второе максимальное напряжение 0,826 МПа, уровень деформации, соответствующий второму максимуму 135,98 %; для продольных образцов: уровень деформации, при котором начинается резкое увеличение напряжения 18,42 % (Tenia) и 27,12 % (Haustra), уровень деформации разрушения - 205,49 % (Tenia) и 205,65 % (Haustra), первое максимальное напряжение 1,188 МПа (Tenia), уровень деформации, соответствующий первому максимальному напряжению, - 40,94 % (Tenia), минимум напряжения 0,168 МПа (Tenia), уровень деформации, соответствующий минимуму напряжения 52,43 % (Tenia), второму максимуму напряжения - 0,652 МПа (Tenia) и 0,533 МПа (Haustra), уровень деформации, соответствующий второму максимуму - 173,45 % (Tenia) и 162,78 % (Haustra). Egorov et al., 2002 [13] Человек Касательные модули вдоль разных направлений: E1=59,9 кПа, E2=147 кПа и E3=37,3 кПа Qiao et al., 2005 [66] Свинья Предел прочности при растяжении (0,58; 0,87 МПа), относительное удлинение при разрушении (113,19; 62,81 %) и модуль упругости (1,83; 5,18 МПа) для образцов свиньи и человека, соответственно. Christensen et al., 2015 [8] Свинья и человек Испытания на одноосное растяжение. Предельные растягивающие напряжения составляли 233 кПа. Конечные растяжения при разрыве составили 1,43 . Gong et al., 2017 [18] Мышь Динамические испытания на равноосное растяжение. Средняя деформация разрушения Грина - Лагранжа - 0,158 и 0,139 в окружном и продольном направлениях соответственно. Среднее 2-е напряжение Пиолы - Кирхгофа при разрушении - 2,35 МПа в окружном направлении и 3,20 МПа в продольном направлении Howes et al., 2013 [46] Человек В обзорной статье [30] представлен общий анализ аномалий медленных волн ЖКТ и их классификация. В работе [75] описаны методы, которые используются для измерения транзита и сократимости желудочно-кишечного тракта или для оценки морфологических изменений нервно-мышечных структур желудочно-кишечного тракта. В работе [22] сокращения в модели двенадцатиперстной кишки аппроксимируются функцией Хевисайда. В экспериментальных исследованиях [29; 91] были зарегистрированы сверхмедленные анальные волны частотой 1,14 сокращений в мин, которые соответствуют результатам исследования [91]. Давление в толстом кишечнике Для построения модели толстого кишечника требуется ряд характеристик, таких как давление, массовый или объемный расход, объем кишечника, его геометрические размеры и др. Первые попытки измерить давление в толстой кишке при дивертикулезе были предприняты в [56]. В работе [32] было установлено максимальное давление, измеренное при колоноскопии в пределах 34 - 138 мм рт. ст., а среднее давление - 22 мм рт. ст. Исследования давления внутри толстого кишечника [76] при накачивании воздухом и корреляция этого параметра с болевыми ощущениями показало величину 38,3 мм рт. ст. В работе [90] определена величина минимального ректального давления для животных (Pmin = 0,64 кПа). В работе [29] исследовались сверхмедленные анальные волны (частотой 1,14 сокращений/мин) и анальное давление. Анальное давление в группе со сверхмедленными волнами (106,0 ± 37,0 см вод. ст.) было значительно выше, чем в группе без проблем с дефекацией (56,0 ± 27,0 см вод. ст.) или у пациентов с запорами без сверхмедленных волн (55,0 ± 26,0 см вод. ст.). Давление в слепой кишке находилось в диапазоне 60 - 80 мм рт. ст., в центральной части толстой кишки - в диапазоне 20 - 40 мм рт. ст., а в непосредственной близости от прямой кишки оно находилось в диапазоне 60 - 100 мм рт. ст. согласно [84], определяемому с помощью беспроводной капсулы (SmartPill). Анализ колоректальных повреждений в результате трансанальной баротравмы высокого давления был проведен в работе [5]. Авторы утверждают, что давление анального сфинктера в состоянии покоя составляет 40 - 80 мм рт. ст., а безопасный уровень воздушного потока при колоноскопическом исследовании должен составлять 1,46 л/мин (при давлении 80 мм рт. ст.). В работе [15] также утверждается, что давление в прямой кишке является физиологичным. Амплитуда изолированных переходных процессов давления варьировалась от 5 до 230 мм рт. ст. в моделях давления с манометрией в толстом кишечнике в исследовании [7]. В табл. 4 приведены основные результаты исследования изменения давления в толстом кишечнике. В работе [50] определены значения средней амплитуды давления высокоамплитудных волн в разных отделах толстого кишечника. В проксимальном отделе толстой кишки средняя амплитуда составила 75,3 мм рт.ст. Для волн, которые начинались в проксимальном отделе толстой кишки и заканчивались в поперечной или нисходящей ободочной кишке, средняя амплитуда составила 87,9 мм рт. ст., а для волн, которые начинались и заканчивались в поперечной или нисходящей ободочной кишке, средняя амплитуда составила 109,3 мм рт. ст. Исследование [62] установило разницу значений давления между мужчинами и женщинами, оцениваемых при аноректальной манометрии. Давление в состоянии покоя составило для мужчин 78,28 мм рт. ст., для женщин - 6,51 мм рт. ст и давление полного сжатия для мужчин - 231,27 мм рт. ст., для женщин - 145,63 мм.рт.ст. Также авторы установили, что средняя длина функционального анального канала была 2,85 см у мужчин и 2,45 см у женщин. Средние значения параметров аноректальной манометрии приведены в табл. 5. Таблица 4 Изменение давления в толстом кишечнике Параметр Одновременные волны давления Антеградные волны давления Ретроградные волны давления Возникновение (/30 мин) 14,7 0,9 0,6 Максимальная частота (/мин) 1,4 2,0 0,4 Средняя длительность (с) 15,5 16,1 21,9 Расстояние распространения (см) 32,3 29,3 26,0 Средняя амплитуда (мм рт. ст.) 14,6 20,5 15,8 Таблица 5 Средние значения параметров аноректальной манометрии Параметр Давление Давление в покое, мм рт. ст. 70,89 Общее давление сжатия, мм рт. ст. 188,45 Сокращение внешнего сфинктера, мм рт. ст. 116,83 Функциональный анальный канал, см 2,65 Заключение Давление внутри толстого кишечника зависит от многих параметров и условий эксперимента. На этот параметр влияют пол испытуемого, место измерения и момент активности толстой кишки. Моделирование потока кала по толстому кишечнику является одной из наиболее важных задач в гастроэнтерологии. Математическая модель толстого кишечника позволит сформулировать патологию физиологических процессов, участвующих в ней. Для реализации такого подхода требуется хорошая экспериментальная база для изучения механических характеристик потока. Исследования в этом направлении ведутся, но многие аспекты остаются неисследованными. Биомеханическое моделирование толстого кишечника является одной из перспективных областей исследований, которые будут изучены в будущих работах.

About the authors

M. I Shmurak

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 195

PDF (Russian) - 660

Refbacks

  • There are currently no refbacks.

Copyright (c) 2022 Шмурак М.И.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies