Full Text

Введение На развитие и функционирование опорно-двигательной системы человека влияют постоянно изменяющиеся нагрузки, имеющие как физиологическую, так и патологическую природу (например, такое изменение может происходить из-за врачебного вмешательства). Известно [1-14], что подобные изменения нагрузки имеют отдаленные во времени последствия и в первую очередь отражаются на организации структуры костной ткани [4]. Губчатая костная ткань является неоднородным пористым анизотропным материалом, механические свойства которого во многом определяются его внутренней архитектурой и подчиняются тем же принципам, что и инженерные конструкции [4, 6, 9, 12]. В кости непрерывно происходят различные адаптационные процессы, в частности под действием новой нагрузки в структуре трабекулярной костной ткани могут начаться адаптационные изменения, которые затем окажут влияние на кортикальную костную ткань. Известно, что адаптационные процессы, протекающие в организме человека, подчиняются закону Вольфа (Wolff’s law) [9]. Применительно к костной ткани закон Вольфа гласит, что архитектура губчатой кости в локальной области структурно приспосабливается к местному напряженному состоянию костной ткани [2, 4, 6, 9, 10]. При этом трабекулы в процессе перестройки выстраиваются вдоль главных направлений тензора напряжений [4, 6, 9, 10]. Несмотря на существующие в настоящее время методы, позволяющие получать изображения губчатой костной ткани in vivo [15], их применение в клинической практике ограниченно возможно и практически не обладает прогностической функцией. Таким образом, возникает необходимость в прогнозировании долговременного отклика структуры на внешнее вмешательство методами математического моделирования. Ранее автором была предложена и рассмотрена модель перестройки губчатой костной ткани с учетом взаимовлияния напряженно-деформированного состояния и структурных изменений, позволяющая спрогнозировать долговременный отклик структуры на внешнее вмешательство методами математического моделирования [2-4, 7]. Данная модель может быть использована для решения прикладных задач о перестройке губчатой костной ткани в различных отделах скелета человека [2-4, 7] в прогностических целях, при этом особый интерес представляет моделирование структуры в нижней челюсти человека. 1. Материалы и методы Существует ряд математических моделей, описывающих адаптационные изменения в трабекулярной костной ткани, подчиняющихся закону Вольфа [2-4, 6-14]. В данной работе рассматривается процесс перестройки губчатой кости как упругого анизотропного тела, обладающего структурой согласно эволюционным соотношениям, впервые предложенным в работе [8] и примененным при практических расчетах в серии исследований [2-4, 6]. На мезо- и макроуровнях математически удобнее описывать губчатую структуру посредством тензора, пренебрегая свойствами отдельных трабекул [14]. Известно, что одним из наиболее удачных способов описания локальной структуры многих пористых и композиционных материалов является симметричный, положительно определенный тензор второго ранга, названный тензором структуры [2-4, 6-9]. Тензор структуры, построенный в соответствии с ранее рассмотренной процедурой [4, 9], позволяет компактно описать неоднородность костной структуры; данный тензор также позволяет сформулировать соотношения, способные количественно описать свойства материала с учетом его строения и адаптационные изменения, происходящие в его микроструктуре [6, 9]. В предыдущих работах [2, 4, 9, 10] было показано, что существует определяющее соотношение, включающее в себя тензор структуры и связывающее напряженно-деформированное состояние в губчатой костной ткани с ее строением. А именно: (1) где - девиатор тензора структуры, нормированный таким образом, что [8]; e - изменение доли твердого объема кости; g1-g6 - константы [2, 9], имеющие размерность [ГПа]. Эти константы были определены в работе [13] после серии экспериментов на различных образцах губчатых костей человека и крупного рогатого скота. Также на основе закона Вольфа [2, 6, 9] был описан предполагаемый механизм перестройки трабекулярной микроструктуры, происходящий вследствие адаптационных процессов и подробно рассмотренный в работах [2, 4, 8, 9]. Были получены следующие кинетические уравнения: (2) и (3) где h1-h4 и f1-f3 - константы [2, 9], имеющие размерность [сут-1], определяемые эмпирически так, чтобы перестройка костной ткани происходила за время адаптации губчатой кости к новым нагрузкам, т.е. примерно за 160 дней [3, 4, 6, 8, 9]. Для соотношений (1)-(3) была представлена постановка начально-краевой задачи о перестройке трабекулярной костной ткани [2, 4]. Данная постановка может быть использована для решения вопросов об исследовании напряженно-деформированного состояния губчатой костной ткани и протекающих в ней процессов адаптации [3, 4, 7]. 2. Результаты Для представленной математической модели перестройки губчатой костной ткани был разработан и протестирован численный алгоритм ее решения. Данный алгоритм был реализован в виде комплексов проблемно-ориентированных программ [3, 4, 7], в серии вычислительных экспериментов была осуществлена идентификация и верификация параметров модели и на ряде примеров [3, 4, 7] показана эволюция губчатой костной ткани при изменении напряженно-деформированного состояния. Для верификации разработанного алгоритма была решена классическая задача [8, 9], подробно разобранная в работе [3]. Результаты, полученные по истечении 160 сут, свидетельствуют, что в результате адаптационных процессов для нового напряженно-деформированного состояния была получена новая костная микроструктура [3], при этом время и характер сходимости, полученные из решения, совпадают с известными [8, 9] результатами. Заметим, что представленный пример является комплексным: действующее на трабекулярную микроструктуру давление, изначально являясь неоднородным, изменяется как по величине, так и по направлению. Данный пример сложен для детального анализа и необходим лишь для сравнения с результатами работы [8] и верификации представленного решения. Серия примеров [3, 4, 7], рассмотренная автором, показала хорошую работоспособность модели. Далее рассмотрим задачу о формировании губчатой костной ткани в ветви нижней челюсти человека и дальнейшей ее адаптации к испытываемым нагрузкам. Известно, что трабекулы в альвеолярном отростке и ветви нижней челюсти выстраиваются вдоль линий главных напряжений, т.е. вдоль главных направлений тензора напряжений и образуют арочную архитектуру в ветви нижней челюсти [1, 5] (рис. 1, 2). Рис. 1. Губчатая костная ткань нижней челюсти [5] а б Рис. 2. Структура ветви нижней челюсти (а) [1], траектории ветви нижней челюсти (б) [1] и главные направления тензора структуры (а и б) [4] Ранее было показано [4, 6, 7], что тензор структуры с достаточной степенью точности может описывать структуру в ветви нижней челюсти, при этом главные оси тензора структуры совпадали с траекториями трабекул, полученными в работе [1], и траекториями напряжений, полученными в работе [16] методом фотоупругости. Таким образом, представленная модель может описать особенности формирующейся под воздействием изменяющегося биомеханического давления структуры костной ткани в различных отделах зубочелюстной системы. Общая постановка данной задачи может быть представлена в следующем виде (рис. 3): (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) Рис. 3. Постановка задачи о перестройке губчатой костной ткани в ветви нижней челюсти человека В начальный момент времени структура губчатой костной ткани в ветви нижней челюсти считалась однородной и была задана нормированным шаровым тензором структуры. Такая структура может соответствовать ситуации, когда на нижнюю челюсть не оказывали влияния никакие нагрузки, т.е. челюсти новорожденного ребенка. В дальнейшем при приложении нагрузок структура нижней челюсти должна приобрести характерный вид (см. рис. 2, б). Расчет напряженно-деформированного состояния производился с помощью метода конечных элементов в пакете ANSYS, определение параметров структуры (доли твердого объема кости и девиатора тензора структуры) осуществлялось в MATLAB [4, 7]. Была построена двумерная конечно-элементная модель нижней челюсти [4, 7, 17], показанная на рис. 4. Для ее задания использовался двумерный конечный элемент Plane182. Численное интегрирование кинетических уравнений осуществлялось методом Рунге - Кутты четвертого порядка. Рис. 4. Конечно-элементная модель нижней челюсти Были рассмотрены различные варианты приложения нагрузки в начальный момент времени [4, 17]. Необходимо отметить, что наиболее наглядно начальное напряженное состояние позволял представить вариант области без учета зубного ряда. В результате получилось, что показанная на рис. 5 интенсивность напряжений в ветви нижней челюсти качественно совпадает с траекториями напряжений, полученными в работах [4, 16]. Результаты также совпадают с результатами из работы [17] для данной области без учета параметров структуры. Что касается распределения трабекул в ветви нижней челюсти, то, как показали вычислительные эксперименты, отведенных 160 дней недостаточно для того, чтобы главные оси девиатора тензора структуры стали соосными с главными осями тензора напряжений. Вероятно, структуре требуется больше времени для ее адаптации к физиологическому напряженному состоянию. Это соответствует известным фактам: полностью трабекулярная структура формируется лишь к 16-18 годам [5]. Рис. 5. Интенсивность напряжений по Мизесу Однако же по истечении 200 дней в ветви нижней челюсти можно наблюдать стремление тензора структуры ориентироваться таким образом, чтобы сформировать структуру, показанную на рис. 2, а. По истечении 300 дней можно заметить формирование характерной для данной области арочной структуры (рис. 6, а). Для удобства представления главные направления девиатора тензора структуры показаны интегрально и наложены на картину интенсивности напряжений, наблюдаемой в ветви нижней челюсти. а б Рис. 6. Главные направления девиатора тензора структуры в ветви нижней челюсти (а) и упрощенное представление распределения напряжений: главные направления, полученные в ANSYS [7] (б) Можно заметить также, что ориентация главных направлений тензора структуры в ветви нижней челюсти, которые показывают, как ориентированы трабекулы внутри костной ткани, качественно совпадает с главными направлениями тензора напряжений (рис. 6, а, б). Более точное определение параметров модели и характерного времени перестройки является предметом нашего дальнейшего исследования. Полученная по истечении 160 сут доля твердого объема кости была подробно рассмотрена в работах [4, 7] и в данной работе не приводится. Отметим, что для заданных условий нагружения результаты распределения плотности костного вещества качественно совпадают с распределением пористости, которое было показано на рис. 1. Заключение В работе рассмотрена модель перестройки губчатой костной ткани с учетом определяющего соотношения и кинетических уравнений феноменологической теории [4, 9], описывающих изменения трабекулярной архитектуры на мезоуровне с помощью тензора структуры. Осуществлена постановка задачи о перестройке губчатой костной ткани и рассмотрена перестройка костной ткани в ветви нижней челюсти. Полученные результаты достаточно схожи с траекториями напряжений и структуры, известными из литературы. Таким образом, представленная модель отражает закономерности формирования костной ткани в различных отделах скелета человека и может быть использована, например, для описания отклика зубочелюстной системы на изменяющееся биомеханическое давление. В частности, представляет интерес моделирование поведения костной ткани в районе вживления имплантата, например при синдроме Попова - Годона [5] или какой-либо функциональной травме [1, 6], с последующим замещением дефектов зубного ряда. Одной из важнейших проблем современной стоматологии является атрофия челюстных костей при частичной или полной вторичной адентии. Потеря постоянных зубов и отсутствие своевременного восстановительного лечения являются основными этиологическими факторами развития «патологического» типа атрофии альвеолярного отростка. Решение подобных задач в дальнейшем позволит планировать лечение различных заболеваний зубочелюстной системы и глубже понять механизмы ее функционирования. Появляется возможность прогнозировать результаты отдаленного врачебного вмешательства при индивидуальном подходе к каждому пациенту.

About the authors

A. A Kichenko

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 22

PDF (Russian) - 10

Refbacks

• There are currently no refbacks.