Математическое моделирование диска височно-нижнечелюстного сустава на основе теории М.А. Био

Аннотация


Заболеваниям височно-нижнечелюстного сустава подвержено около 40 % людей в мире в возрасте 20-40 лет. В настоящее время этимология заболеваний височно-нижнечелюстного сустава изучена слабо, существует лишь их примерная медицинская классификация. Анализ литературы показывает, что, как правило, практически ни одно заболевание компонентов сустава не обходится без их механического повреждения. Диск височно-нижнечелюстного сустава играет ключевую роль в работоспособности самого сустава и зубочелюстного комплекса в целом. Основной фактор, влияющий на работоспособность диска, - величина механических напряжений, вызванных различными нагрузками. Как правило, лечение заболеваний такого рода носит неинвазивный характер, но в сложных случаях заболевания также используются хирургические методы. Во избежание операционного вмешательства и для профилактики заболеваний на начальной стадии необходимо провести биомеханическое математическое моделирование диска височно-нижнечелюстного сустава. Изучение гистологии и физиологии диска показывает, что оптимальная модель для описания механического поведения диска под нагрузкой - изотропная пороупругая модель, предложенная М.А. Био. Она позволяет оценить напряженное состояние диска и объяснить с точки зрения механики роль синовиальной жидкости внутри него. В работе на основе постановки задачи линейной изотропной теории пороупругости построена модель поведения диска височно-нижнечелюстного сустава под нагрузкой. Решение сводится к начально-краевой задаче для нестационарного параболического уравнения. Из анализа результатов видно, что распределение давления жидкости в диске при приложении постоянной силы позволяет снизить величину механических напряжений в упругом хрящевом скелете, защищая его от разрушения. Исследована зависимость величины давления синовиальной жидкости от коэффициента гидравлической проницаемости, для чего был проведен анализ результатов решения для модели с различными величинами этого коэффициента.

Полный текст

Введение Диск представляет собой сложную структуру - твердый хрящевой матрикс, насыщенный синовиальной (суставной) жидкостью - и играет ключевую роль в работоспособности самого сустава и зубочелюстного комплекса в целом. Основной фактор, влияющий на работоспособность диска, - величина механических напряжений, вызванных различными нагрузками [1]. Центральная зона диска подвержена наибольшим сжимающим нагрузкам в процессе жизнедеятельности. Более того, ввиду отсутствия в этой зоне кровеносных сосудов ее самостоятельная регенерация при механическом повреждении практически невозможна. Во избежание операционного вмешательства и для профилактики заболеваний диска на начальной стадии необходимо провести биомеханический анализ центральной зоны диска височно-нижнечелюстного сустава на предмет определения величины и распределения напряжений. 1. Материалы и методы На основании данных об анатомии, гистологии, физиологии, а также работ других авторов по биомеханическому моделированию диска височно-нижнечелюстного сустава наиболее адекватной математической моделью для описания поведения диска под нагрузкой принята изотропная пороупругая модель [2]. Такая модель позволяет учесть пористую упругую хрящевую структуру диска, а также свойства и количественное содержание синовиальной жидкости внутри него. Для моделирования поведения центральной зоны диска височно-нижнечелюстного сустава под нагрузкой и определения величины механических напряжений в диске, а также их распределения с течением времени необходимо сформулировать постановку задачи и найти зависимость напряжений от величины гидравлического давления синовиальной жидкости внутри него. Для решения такой задачи использована постановка классической задачи Манделя, которая ранее ставилась для насыщенных влагой грунтов [3, 4]. Постановка этой задачи впервые использовалась для биомеханического моделирования в данной работе. Решение задачи проводилось методом преобразования Лапласа, в отличие от первоисточника, в котором задача решалась методом Фурье, причем результаты решения совпали. Решение методом Фурье, примененным в оригинальной постановке задачи Манделя, имеет ряд особенностей, которые реализуются не во всех случаях. Универсальным способом аналитического решения таких задач является метод преобразования Лапласа. Представим центральную зону диска височно-нижнечелюстного сустава в виде плоского прямоугольного тела длиной 2а. Мыщелок нижней челюсти и чешуйчатую часть височной кости будем рассматривать как две жесткие непроницаемые пластины, окружающие диск сверху и снизу (рис. 1). Объектам анализа будет являться плоско-деформированное состояние тела. Синовиальная жидкость свободно вытекает из центра диска к лимфоузлам на периферии, закрепления на краях отсутствуют. В начальный момент диск не нагружен, далее к центру нижней и верхней пластины прикладывается постоянная сжимающая сила в направлении оси y, моделирующая жевательную нагрузку и равная F. Касательные напряжения в хрящевом скелете равны нулю на участке контакта пороупругого тела и пластин из-за того, что коэффициент трения синовиальной жидкости близок к нулю. Наличие же переработанной синовиальной жидкости (мукоида) в центральной зоне диска позволяет рассматривать его структуру как однородную. Рис. 1. Пороупругая модель диска височно-нижнечелюстного сустава Решение модели сводится к начально-краевой задаче для нестационарного параболического уравнения и производится методом преобразования Лапласа. 2. Результаты На рис. 2 представлены графические распределения гидравлического давления в пороупругом теле и напряжений в твердом упругом скелете диска височно-нижнечелюстного сустава с течением времени. Для упрощения расчетов задача решена в относительных величинах: гидравлическое давление и напряжения в скелете нормированы по их начальным значениям, координата нормирована по длине модели диска, время - по квадрату длины и значению коэффициенту консолидации. В начальный момент времени давление в центре тела близко к начальному, но превышает его на некоторую величину. Данный эффект играет положительную роль в функционировании диска - давление в порах и напряжения в упругом скелете связаны: чем выше гидравлическое давление, тем ниже напряжения хрящевого скелета, следовательно, такое возрастание давления не дает диску разрушиться под действием мгновенно приложенной нагрузки. Далее с течением времени наблюдается падение давления по всему телу по экспоненте в связи со свободным вытеканием жидкости по краям. Рис. 2. Перераспределение: а - безразмерного давления синовиальной жидкости; б - безразмерных напряжений в центральной зоне диска височно-нижнечелюстного сустава с течением времени. Кривая 1 - при , кривая 2 - при , кривая 3 - при , кривая 4 - при , кривая 5 - при Поскольку давление в жидкости падает, напряжения в диске возрастают. Такое поведение напряжений и давления в диске способствует равномерному распределению нагрузки, предотвращает его от механических повреждений, а также обеспечивает питание зон, свободных от кровеносных сосудов. Различные механические упругие константы для диска височно-нижнечелюстного сустава получены многими авторами по итогам различных экспериментов. Но применение пороупругой модели диска требует рассмотрения не только упругих параметров, но и коэффициента гидравлической проницаемости диска, входящего в закон Дарси. От величины коэффициента проницаемости непосредственно зависит возможность диска сопротивляться нагрузкам, его жесткость. Другими словами, величина коэффициента проницаемости влияет на способность диска сопротивляться сжатию и, следовательно, его разрушению. Низкая величина коэффициента говорит о том, что на перераспределение давления и напряжения в теле требуется значительно большее количество жидкости. В случае с высоким значением коэффициента проницаемости быстрый обмен жидкости в хряще приводит к существенному уменьшению жесткости и, соответственно, уменьшению напряжения в структуре. Значение коэффициента проницаемости для хрящей разных видов зависит от их функциональности и назначения. Рис. 3. Распределение безразмерного давления синовиальной жидкости в диске височно-нижнечелюстного сустава при разных значениях коэффициента гидравлической проницаемости. Кривая 1 - для диска височно-нижнечелюстного сустава; кривая 2 - для гиалинового хряща; кривая 3 - для мениска коленного сустава: а - при , б - при , в - при , г - при , д - при Для подтверждения этой гипотезы был проведен эксперимент по определению коэффициента гидравлической проницаемости диска височно-нижнечелюстного сустава на суставных дисках свиньи. По итогам эксперимента получено среднее значение коэффициента гидравлической проницаемости диска височно-нижнечелюстного [5]. Для сравнения: значения коэффициента проницаемости для суставного хряща больше в 0,3 раз, для мениска коленного сустава - примерно в 1,2 раза (рис. 3). На начальном этапе, когда течение жидкости еще развито недостаточно, разница в значениях давления незначительна. Далее с течением времени эта разница существенно возрастает. В последний рассматриваемый момент времени величина давления в гиалиновом хряще примерно в 20 раз больше, чем в диске височно-нижнечелюстного сустава и мениске коленного сустава. Из графических решений на рис. 3 видно, что величина коэффициента гидравлической проницаемости сильно влияет на значение гидравлического давления и, как следствие, величину механических напряжений в твердом хрящевом скелете. Это, в свою очередь, подтверждает, что сопротивляемость сжатию напрямую зависит от значения коэффициента проницаемости. Данные результаты служат подтверждением нашего предположения о том, что коэффициент проницаемости для различных видов хрящей, в том числе для диска височно-нижнечелюстного сустава, необходимо уточнять экспериментально.

Об авторах

В. М. Тверье

Пермский национальный исследовательский политехнический университет

А. С. Миленин

Пермский национальный исследовательский политехнический университет

Список литературы

  1. Тверье В.М., Симановская Е.Ю., Няшин Ю.И. Механический фактор развития и функционирования зубочелюстной системы человека // Российский журнал биомеханики. - 2005. - Т. 9, № 2. - С. 34-42.
  2. Экспериментальное изучение фильтрационных свойств и структурных особенностей дисков височно-нижнечелюстных суставов свиней / М.Ю. Няшин [и др.] // Российский журнал биомеханики. - 2002. - Т. 6, № 3. - С. 33-38.
  3. Coussy O. Poromechanics. - J.Wiley, 2004. - 298 p.
  4. Mandel J. Consolidation des sols // Geotechnique. - 1953. - Vol. 3, № 7. - P. 287-299.
  5. Тверье В.М., Няшин Ю.И. Коэффициент гидравлической проницаемости диска височно-нижнечелюстного сустава: экспериментальное определение // Российский журнал биомеханики. - 2010. - Т. 14, № 2. - С. 28-36.

Статистика

Просмотры

Аннотация - 11

PDF (Russian) - 6

Ссылки

  • Ссылки не определены.

Данный сайт использует cookie-файлы

Продолжая использовать наш сайт, вы даете согласие на обработку файлов cookie, которые обеспечивают правильную работу сайта.

О куки-файлах