Full Text

Введение На развитие и функционирование костной системы человека влияют постоянно изменяющиеся нагрузки, имеющие различную (физиологическую и патологическую) природу. В частности, такое изменение может происходить вследствие врачебного вмешательства. Известно [11, 17, 20, 25-26, 31-33], что подобные изменения нагрузки имеют отдаленные во времени последствия и в первую очередь отражаются на организации структуры костной ткани. Костная ткань делится на компактную (кортикальную) и губчатую (трабекулярную) в зависимости от расположения трабекул. Структурной единицей губчатой костной ткани, образующей ее архитектуру, является трабекула. Губчатая костная ткань является неоднородным пористым анизотропным материалом, механические свойства которого во многом определяются его внутренней архитектурой и подчиняются тем же принципам, что и инженерные конструкции [33, 54]. Губчатая костная ткань живого человека является динамической структурой, в которой непрерывно происходят различные адаптационные процессы. Адаптация (перестройка) структуры трабекулярной костной ткани представляет собой процесс поворота трабекул в сторону воздействия преобладающей силы на рассматриваемую область. В частности, под действием новой нагрузки в структуре трабекулярной костной ткани могут начаться адаптационные изменения, которые затем окажут влияние на кортикальную костную ткань. Механизм перестройки на клеточном уровне реализуется путем резорбции (рассасывания) отдельных старых трабекул остеобластами и образования новых трабекул (остеосинтеза) остеокластами. Известно, что адаптационные процессы, протекающие в организме человека, подчиняются закону Вольфа (Wolff’s law) [69]. Закон Вольфа применительно к живой костной ткани звучит следующим образом: кость приспосабливает свою внешнюю форму и внутреннюю структуру к тем механическим силам, которые она должна выдержать [11, 13, 14, 20, 32, 33, 54]. В частности, закон Вольфа для костной ткани говорит о том, что архитектура губчатой кости в локальной области структурно приспосабливается к местному напряженному состоянию костной ткани [13, 20, 32, 33]. При этом структурная адаптация в живой губчатой кости носит направленный характер: трабекулы располагаются закономерно, сообразно тому, какие внешние нагрузки испытывает данная кость [13, 20, 32, 33]. В частности, в областях, где приложена нагрузка и началась перестройка, происходит увеличение плотности, а трабекулы в процессе перестройки выстраиваются вдоль главных направлений тензора напряжений [13, 20, 32, 33]. Несмотря на существующие в настоящее время методы, позволяющие получать изображения губчатой костной ткани in vivo [47, 51, 52], их применение в клинической практике ограниченно возможно и практически не обладает прогностической функцией. Таким образом, возникает необходимость в прогнозировании долговременного отклика структуры на внешнее вмешательство методами математического моделирования, т.е. необходимость в разработке модели перестройки (адаптации) губчатой костной ткани с учётом взаимовлияния напряжённо-деформированного состояния и структурных изменений. Математическое моделирование перестройки губчатой костной ткани по-прежнему остается одной из актуальных задач биомеханики. Ранее автором данной статьи была описана модель перестройки губчатой костной ткани с учётом взаимовлияния напряжённо-деформированного состояния и структурных изменений, позволяющая спрогнозировать долговременный отклик структуры на внешнее вмешательство методами математического моделирования [13, 14, 64]. Данная модель может быть использована для решения прикладных задач о перестройке губчатой костной ткани в различных отделах скелета человека [3-9, 14, 30] в прогностических целях, в частности применительно к различным отделам зубочелюстной системы человека [31]. Особый интерес представляет моделирование структуры в нижней челюсти человека. Зубочелюстная система человека - объект исследования Нижняя челюсть человека является непарной и единственной подвижной костью лицевого скелета. К нижней челюсти крепится большое количество мышц, приводящих её в движение. Эти особенности определяют сложность строения нижней челюсти и её влияние на развитие лицевого скелета и покрывающих скелет мягких тканей [1, 2, 17]. Строение нижней челюсти с момента рождения и до глубокой старости изменяется в соответствии с испытываемой ею нагрузкой, имеющей различную природу [2]. При этом в процессе развития и роста нижней челюсти структура губчатой костной ткани также изменяется в соответствии с испытываемым костью нагружением, ее напряженным состоянием (рис. 1). Для объективного учета этих изменений необходимо выяснить типовые особенности нормального строения, присущие челюсти в целом и подчиняющиеся различным ее анатомо-физиологическим законам в каждом возрастном периоде. Остановимся на этих изменениях подробнее. Замечено, что тело нижней челюсти имеет ярко выраженную структуру губчатой костной ткани (рис. 2) [16]. Исследования трабекулярной структуры выявляют закономерности ее строения. Костные балки губчатого вещества располагаются по линиям действия сил при нагрузке. Всего выделяют до девяти таких траекторий [26]. Строение ветви нижней челюсти во всех возрастных периодах было подробно изучено в работах [2, 26]. Ветвь нижней челюсти человека растет в длину и ширину неравномерно [2]. Наиболее интенсивные периоды роста ветви в высоту (длину) отмечаются в возрасте от 3 до 4 лет и в 9-11 лет. Периоды значительного роста ветви в ширину наблюдаются в сроки прорезывания больших коренных зубов, т.е. в возрасте 9, 15 и 18 лет [26]. Наибольший прирост массы костного вещества восходящей ветви нижней челюсти происходит в период от 6 до 15 лет. Период от 25 до 45 лет включительно характеризуется наименьшим разбросом количества костного вещества ветви. С 50 лет и до глубокой старости количество костного вещества в восходящей ветви человека убывает. Наибольшее увеличение количества компактного вещества по отношению к губчатому происходит в четыре срока: от 9 месяцев до 1 года; от 2,5 до 4 лет; от 9 до 12 лет и от 15 до 18 лет [26]. Формирование губчатого вещества восходящей ветви по возрастам подтверждает данные о различных типах организации костной ткани в двух соприкасающихся треугольных полях (вентральном и дорсальном) восходящей ветви. Для вентрального треугольного поля типично (в зрелом возрасте) преобладание дугообразных конструкций (в виде так называемых «опрокинутых арок»), для дорсального - веерообразных, расходящихся перекладин [26]. Радиальные, расходящиеся структуры губчатого вещества ветви обнаруживаются уже к моменту рождения; количество их быстро нарастает к 3 годам и в дальнейшем увеличивается как в численном отношении, так и в отношении мощности составляющих их элементов. Дуговые перекладины вентрального отдела восходящей ветви появляются позднее радиальных. К двум годам при наличии 16-18 зубов они впервые прорисовываются и развиваются в период от 3 до 6 лет. Рис. 1. Линии напряжений при нагрузке на фронтальном участке (напряжения возникают во всей челюсти) [16] Рис. 2. Губчатая костная ткань нижней челюсти [16] а б Рис. 3. Структура (а) и траектории (б) ветви нижней челюсти [2] Затем сложившаяся система «опрокинутых арок» как бы созревает, становится более четкой, отшлифовывается в процессе жизнедеятельности в соответствии с меняющимися условиями существования (в том числе вследствие смены молочных зубов постоянными). К 15 годам губчатое вещество вентрального треугольника приближается по своему строению к кости зрелого возраста, а к 18 годам оно почти не отличается от губчатого вещества 25-летнего человека (рис. 3). Выпадение зубов и старческая инволюция проявляются в наибольших изменениях губчатого вещества в вентральном треугольном поле восходящей ветви. Регрессивные процессы старческого возраста возникают в тех отделах кости, которые в онтогенезе оформляются позднее других. Показано, что трабекулы в альвеолярном отростке и ветви нижней челюсти выстраиваются вдоль линий главных напряжений, т.е. вдоль главных направлений тензора напряжений, и образуют арочную архитектуру в ветви нижней челюсти [2, 27]. Данное утверждение согласуется с законом Вольфа для костной ткани. Отсюда видно, что необходимо иметь способ количественного описания формирующейся под воздействием изменяющегося биомеханического давления структуры костной ткани в различных отделах зубочелюстной системы. В частности, при постановке задач об определении напряжённо-деформированного состояния в нижней челюсти человека необходимо учитывать ее внутреннюю структуру [12-15, 24-26]. Материалы и методы С тех пор, как был сформулирован закон Вольфа, накопилось множество экспериментальных данных и математических моделей, описывающих этот механизм [11, 20, 30-41, 43, 44, 46, 53, 54, 56-59, 61, 62, 70]. В данной работе реализуется процесс перестройки губчатой костной ткани как упругого анизотропного тела, оснащенного структурой, согласно эволюционным соотношениям, впервые предложенным в работе [32] и примененным при практических расчетах в серии исследований [3-7, 9, 11-15, 20, 30, 31]. В качестве инструмента численного моделирования используются программный пакет конечно-элементного анализа ANSYS и пакет прикладных программ MATLAB [7]. На уровне органов и тканей (мезо- и макроуровнях) математически удобнее описывать губчатую структуру некоторой области тензорной величиной, пренебрегая свойствами отдельных трабекул. Известно, что одним из наиболее удачных способов [11-13, 20, 22, 26, 33, 38, 39, 45, 48, 49, 54, 56, 57, 60, 61, 63, 65-68] описания локальной структуры многих пористых и композиционных материалов является симметричный, положительно определенный тензор второго ранга, названный тензором структуры (fabric tensor) и обозначенный как [11-13, 20, 22, 26, 33, 60, 63]. Тензор структуры, построенный для губчатой костной ткани в соответствии с ранее описанной методикой [11, 12, 20, 26, 63], позволяет компактно в тензорной форме описать анизотропию костной структуры, причём его главные значения позволяют охарактеризовать распределение материала вдоль главных направлений [14, 15, 26, 64]. Данный тензор позволяет сформулировать соотношения, способные количественно описать свойства материала с учётом его строения и адаптационные изменения, например, костной микроструктуры. Ранее автором была рассмотрена измерительная процедура, позволяющая проводить стереологические исследования [22] для построения тензора структуры, и исследованы свойства тензора структуры и связанного с ним тензора анизотропии [12, 63]. Также была отработана методика экспериментального получения трёхмерных образцов трабекулярной костной ткани и дальнейшего построения тензора структуры [12, 63] - см. рис. 4. Установлено [14, 33, 39, 64], что в состоянии гомеостаза (физиологического равновесия) ориентация трабекул в рассматриваемой области губчатой кости совпадает с главными направлениями тензора напряжений в этой же области, в этом случае тензоры напряжений и структуры должны быть соосны. а б Рис. 4. Характерный вид костного образца после его извлечения из муфельной печи (а) и после его подготовки для построения тензора структуры с наложенным эллипсом структуры (б) [12] Совпадение главных осей тензоров напряжения и структуры имеет место в случае, когда их скалярное произведение коммутативно [20, 27, 33, 39], т.е. выполняется условие, что . (1) В предыдущих работах [13, 64] было показано, что существует определяющее соотношение, включающее в себя тензор структуры и связывающее напряжённо-деформированное состояние в губчатой костной ткани с её строением (трабекулярной микроструктурой). При этом считается, что костный матрикс в губчатой кости (пористом упругом теле) изотропен [32, 33] и вся ее неоднородность связана со строением анизотропной микроструктуры трабекулярной костной ткани. В этом случае анизотропия губчатой кости может быть описана посредством тензора структуры, при этом упругие свойства материала будут зависеть от пористости губчатой кости и от ориентации трабекул. Таким образом, на основе существующих подходов [8, 13, 14, 21, 23, 27, 28, 32-41], можно получить определяющее соотношение, позволяющее описать напряжённо-деформированное состояние губчатой кости с учётом её структуры: , (2) где - девиатор тензора структуры нормированный таким образом, что [8]; e - изменение доли твёрдого объёма кости относительно отсчётной величины ν0; g1-g6 - константы [13, 32], ГПа. Эти константы были определены в работе [62] после серии экспериментов на различных образцах губчатых костей человека и крупного рогатого скота. На основе закона Вольфа [13, 14, 20, 32, 33, 39, 64, 69] был также описан предполагаемый механизм перестройки трабекулярной микроструктуры, происходящий вследствие адаптационных процессов, перманентно протекающих в живой костной ткани. В соответствии с законом Вольфа для губчатой костной ткани [13, 14, 32, 33] были получены эволюционные соотношения, описывающие изменения в строении губчатой кости под воздействием различных нагрузок. Данные кинетические уравнения способны описать изменение девиатора тензора структуры и доли твёрдого объёма кости e. А именно: (3) и , (4) где h1-h4 и f1-f3 - константы [14, 32], сут-1, определяемые эмпирически так, чтобы перестройка костной ткани происходила за время, соответствующее природной действительности. В норме адаптация губчатой кости происходит примерно за 160 дней [14, 20, 33, 54, 64]. Для соотношений (2)-(4) была представлена постановка начально-краевой задачи о перестройке трабекулярной костной ткани. Данная постановка может быть использована для решения вопросов об исследовании напряжённо-деформированного состояния губчатой костной ткани и протекающих в ней адаптационных процессов. Результаты Для представленной в работе математической модели перестройки губчатой костной ткани был разработан и протестирован численный алгоритм ее решения. Данный алгоритм был реализован в виде комплексов проблемно-ориентированных программ [14, 15, 30, 31]. В ряде вычислительных экспериментов были осуществлены идентификация и верификация параметров модели и на ряде примеров показана эволюция губчатой костной ткани при изменении напряжённо-деформированного состояния. Далее будут кратко описаны полученные результаты. Для верификации разработанного алгоритма была решена классическая задача [32, 33], подробно разобранная в [14, 64]. Рассматривалась локальная область губчатой костной ткани (рис. 5), находящаяся в состоянии гомеостаза в течение достаточно долгого промежутка времени (t < 0). В момент времени t = 0 происходит однократное изменение условий нагружения, приводящее к перестройке трабекулярной микроструктуры. Новое напряжённое состояние в рассматриваемой области при этом не изменяется в течение достаточно долгого промежутка времени (t > 0) [14]. Данные, полученные по истечении 160 суток, показаны на рис. 6. В результате адаптационных процессов для нового напряженно-деформированного состояния была получена новая костная микроструктура [14, 64], при этом время и характер сходимости, полученные из решения, совпадают с известными [32] результатами. Заметим, что представленный пример является комплексным: действующее на трабекулярную микроструктуру давление, изначально являясь неоднородным, изменяется как по величине, так и по направлению. Данный пример сложен для детального анализа и необходим лишь для сравнения с результатами работы [32] и верификации представленного решения. Серия примеров [3-6, 15, 64], рассмотренная автором и не представленная в данной работе, показала хорошую работоспособность модели. Далее были рассмотрены классические задачи [18] о растяжении и изгибе консольно-защемленной балки, оснащенной структурой. Полученные результаты также свидетельствуют о работоспособности модели [4, 15, 30, 31]. Рис. 5. Иллюстрация процесса перестройки структуры трабекулярной костной ткани из одного состояния физиологического равновесия в другое [14] а б в г Рис. 6. Изменение компонент девиатора тензора структуры (а, б), плотности (в) и угла κ(t) между главными осями тензора напряжений и девиатора тензора структуры (г) [14] l P Рис. 7. Растяжение неоднородного стержня под действием осевой силы Рис. 8. Главные направления девиатора тензора структуры после окончания перестройки На примере растяжения ортотропного стержня [18] было показано влияние различной ориентации структуры материала на его напряженно-деформированное состояние. В задаче консольно-защемленная балка растягивалась осевой силой P (см. рис. 7), при этом в каждой точке бруса задавалась структура, соответствующая полученному в работах [12, 63] значению тензора структуры для трехмерного образца губчатой костной ткани. В результате первоначально несоосные тензор напряжений и девиатор тензора структуры стали соосными (рис. 8), а новая костная микроструктура изменилась так, как показано на рис. 9 [6, 15, 30]. Значения перемещений, напряжений и деформации также совпали со значениями, определенными в работе [18] для двумерного и трехмерного случаев (рис. 10-11). На примере изгиба консольно-защемленной балки, оснащенной структурой [18], было показано влияние различной ориентации структуры материала на его напряжённо-деформированное состояние [4]. Было также показано качественное совпадение результатов численного моделирования с результатами, полученными в работе [46], где консольно-защемленная балка при изгибе принимает характерную форму языка: утолщается в области заделки и становится тоньше в области приложения изгибающей нагрузки, - а линии главных напряжений образуют характерные траектории (рис. 12). В рассмотренной автором статьи задаче изначально однородная балка изгибалась силой F (рис. 13). а б в г Рис. 9. Изменение структуры (компонент девиатора тензора структуры K1(t) (а), K2(t) (б) и K3(t) (в)) и доли твердого объема e(t) трабекулярной микроструктуры (г) В результате начальный тензор структуры соответствующий однородной структуре, перестал быть шаровым. В процессе перестройки главные оси девиатора тензора структуры стали соосными главным осям тензора напряжений (рис. 14) [3]. Изначально одинаковая плотность балки увеличилась в области консольного защемления и уменьшилась в области приложения нагрузки (рис. 15), что можно интерпретировать как приобретение балкой характерной формы, показанной на рис. 12. Значения перемещений, напряжений и деформации [30] также совпали со значениями, определенными в работе [18] (рис. 16-17). Задача о формировании губчатой костной ткани в ветви нижней челюсти Рассмотрим задачу о формировании и дальнейшей перестройке губчатой костной ткани в ветви нижней челюсти человека. Общая постановка данной задачи может быть представлена в следующем виде (рис. 18): (5) (6) а б Рис. 10. Изолинии перемещений в неоднородном стержне после завершения перестройки для двумерного (а) и трехмерного (б) случаев (7) (8) (9) (10) (11) (12) (13) (14) В начальный момент времени структура губчатой костной ткани в ветви нижней челюсти считалась однородной и была задана шаровым тензором структуры, т.е. . а б Рис. 11. Перемещения точек средней линии (а) и относительная погрешность (б) численного решения в сравнении с аналитическим [5] Рис. 12. Адаптация формы консольно-защемленной балки, нагруженной изгибающей силой [46] Рис. 13. Изгиб консольно-защемленной балки, оснащенной в начальный момент времени однородной структурой (или ) Рис. 14. Сравнение направления главных осей тензора напряжений и девиатора тензора структуры после завершения перестройки Рис. 15. Изменение плотности (доли твердого объема) балки после завершения перестройки: - более плотные зоны; - менее плотные зоны Рис. 16. Изолинии перемещений в неоднородном стержне после завершения перестройки а б в Рис. 17. Сравнение аналитического [4] и численного решений: перемещения точек средней линии Uy (а), напряжения σz (б) и деформации εz (в) в сечении, перпендикулярном средней линии при z = l / 2 Рис. 18. Постановка задачи о перестройке губчатой костной ткани в ветви нижней челюсти человека а б в Рис. 19. Сравнение траекторий напряжения, полученных методом фотоупругости [55] (а и б), траекторий ветви нижней челюсти [2] и главных направлений тензора структуры (в) [27] Такая структура может соответствовать ситуации, когда на нижнюю челюсть не оказывали влияния никакие нагрузки, т.е. челюсти новорожденного ребенка. В дальнейшем при приложении нагрузок структура нижней челюсти должна приобрести характерный вид (см. рис. 3). Ранее было показано [20, 26], что тензор структуры с достаточной степенью точности может описывать структуру ветви нижней челюсти, при этом главные оси тензора структуры совпадали с траекториями трабекул, полученными в работе [2], и траекториями напряжений из работы [55] (рис. 19). а б Рис. 20. Интенсивность напряжений по Мизесу (а) и ее сравнение с данными из работы [19] (б) Расчёт напряженно-деформированного состояния производился с помощью метода конечных элементов в пакете ANSYS [7, 31], определение параметров структуры (доли твердого объема кости и девиатора тензора структуры) осуществлялось в MATLAB. Была построена двумерная конечно-элементная модель нижней челюсти, её размеры были определены из работы [2], свойства костного матрикса взяты из работы [50], координаты точек прикрепления жевательных мышц и нагрузка от их действия на кость - из работ [10, 19]. Для задания конечно-элементной сетки использовался двумерный конечный элемент Plane182, поддерживающий задание анизотропных свойств. Численное интегрирование кинетических уравнений осуществлялось методом Рунге - Кутта четвёртого порядка. Были рассмотрены различные варианты приложения нагрузки в начальный момент времени. Необходимо отметить, что наиболее наглядно начальное напряженное состояние позволял представить вариант области без учета зубного ряда (нагрузка прикладывается непосредственно к кости челюсти между предполагаемым нахождением второго и первого премоляра) - см. рис. 20. а б Рис. 21. Упрощенное представление распределения напряжений: главные направления, полученные в ANSYS (а) [31], и траектории, в среднем наблюдаемые в челюсти (б) [50] Видно, что полученная картина интенсивности напряжений в ветви нижней челюсти (см. рис. 20, а) качественно совпадает с траекториями напряжений, полученными методом фотоупругости [55] (см. рис. 1 и 19, а и б). Полученные результаты совпадают также с данными из работы [19] для заданной области без учета параметров структуры (см. рис. 20, б). Ещё одним наглядным и удобным для последующих вычислений способов иллюстрации напряжённого состояния является изображение главных направлений тензора напряжений, которые показывают, как ориентированы трабекулы внутри костной ткани (рис. 21). Полученная по истечении 160 суток доля твердого объема кости показана на рис. 22. Результаты распределения плотности костного вещества, полученные для геометрии челюсти с учетом зубов и без них, при рассматриваемых условиях нагружения качественно совпадают. Видно сходство в распределении пористости в сравнении с рис. 3, а. Что касается распределения трабекул в ветви нижней челюсти, то, как показали вычислительные эксперименты, отведенных 160 дней недостаточно для того, чтобы главные оси девиатора тензора структуры стали соосными с главными осями тензора напряжений. По всей видимости, структуре в ветви нижней челюсти требуется больше времени для ее адаптации к физиологическому напряженному состоянию. Это соответствует известным фактам: полностью трабекулярная структура формируется лишь к 16-18 годам. Определение параметров модели и времени перестройки является предметом нашего дальнейшего исследования. Заключение В работе рассмотрена модель перестройки губчатой костной ткани с учётом определяющего соотношения и кинетических уравнений феноменологической теории [12], описывающих изменения трабекулярной архитектуры на мезоуровне с помощью тензора структуры. Осуществлена постановка задачи о перестройке губчатой костной ткани, и на ряде модельных примеров показана ее работоспособность. Рис. 22. Распределение доли твердого объема кости по истечении 160 суток а б Рис. 23. Нормальная структура губчатого вещества нижней челюсти (препарат) мужчины 25 лет (а) и структура губчатого вещества нижней челюсти (препарат) мужчины соответственного возраста при синдроме Попова - Годона (б) [42] Рассмотрена также перестройка костной ткани в ветви нижней челюсти. Полученные результаты достаточно схожи с траекториями напряжений и структуры, известными из литературы. Таким образом, представленная модель отражает закономерности формирования костной ткани в различных отделах скелета человека и может быть использована, например, для описания отклика зубочелюстной системы на изменяющееся биомеханическое давление. В частности, представляет интерес моделирование поведения костной ткани в районе вживления имплантанта, например при синдроме Попова - Годона [25, 26] (рис. 23) или какой-либо функциональной травме [7] с последующим замещением дефектов зубного ряда. Одной из важнейших проблем современной стоматологии является атрофия челюстных костей при частичной или полной вторичной адентии. Потеря постоянных зубов и отсутствие своевременного восстановительного лечения являются основными этиологическими факторами развития «патологического» типа атрофии альвеолярного отростка. Известно, то при потере зубов губчатая костная ткань приобретает характерное диффузное строение, трабекулы располагаются хаотично, что может быть представлено в виде шарового тензора структуры (см. рис. 23 и 24, а). В дальнейшем, при замещении дефекта зубного ряда имплантатом, трабекулы выстраиваются перпендикулярно образующей корня зуба (или штифта имплантата) - см. рис. 24, б. То же можно сказать и про ортодонтическое перемещение зуба: однородная губчатая структура должна приобрести известный в ортодонтии вид (см. рис. 24). В этом случае граничные и начальные условия примут следующий вид: Решение подобных задач в дальнейшем позволит планировать лечение различных заболеваний зубочелюстной системы и глубже понять механизмы ее функционирования. Появляется возможность прогнозировать результаты отдаленного врачебного вмешательства при индивидуальном подходе к каждому пациенту. а б Рис. 24. Постановка задачи о перестройке костной ткани в случае восстановлении зубного ряда при синдроме Попова - Годона

About the authors

A. A Kichenko

References

Statistics

Views

Abstract - 82

PDF (Russian) - 18

Refbacks

• There are currently no refbacks.