Анализ деформационного поведения элементов зубочелюстной системы и однослойной защитной каппы в рамках контактного взаимодействия

• [8]    и т.д. На протяжении последних десятилетий можно отметить особый интерес отечественных и зарубежных ученых к вопросам биомеханического моделирования задач стоматологии: численный анализ влияния неоднородности зубных тканей на центр вращения [9], моделирование начальных перемещений элементов зубочелюстной системы под воздействием ортодонтических нагрузок [10], экспериментальное и численное исследование поведения материалов защитной протетической конструкции [11], численное моделирование прикуса [12], анализ контактного взаимодействия пары зубов с защитной зубной шиной [13] и др.

• 2.    ПОСТАНОВКА ЗАДАЧИ

Одним из актуальных направлений исследования является проблема профилактики травм зубочелюстной системы с использованием индивидуальных средств защиты в виде протети-ческих конструкций, которая широко отражена в [14-16 и др.]. Которая формирует новый класс задач биомеханики зубочелюстной системы. Большое количество исследований данной отрасли посвящено анализу влияния протетиче-ских конструкций на зубочелюстную систему человека и снижение нагрузки на зубной ряд при занятии спортом и тяжелым физическим трудом [8, 11, 14 др.]. Другим актуальным направлением исследования является анализ свойств материалов спортивных зубных шин [11, 15, 17-19 и др.] и их влияние на элементы зубочелюстной системы.. Например, в работе [15] большое внимание уделено анализу различных вариантов средств индивидуальной защиты в виде зубных шин, изготовленных на основе этиленвинилацетата. Авторским коллективом [11] выполнена серия натурных экспериментов по анализу свойств материалов трехслойной защитной шины: Eva и А-силикон. Можно отметить, что в связи с большой важностью использования защитных конструкций в виде зубных шин с целью профилактики травматизма зубочелюстной системы возникла необходимость в компьютерном моделировании биомеханического поведения, как самих конструкций, так и материалов из которых они изготовлены. В настоящий момент можно выделить ряд задач, на решение которых направлены исследования данной научной области: биомеханический анализ влияния свойств материалов зубных шин на деформационное поведение твердых тканей зубов; математическое моделирование контактного взаимодействия зубных шин разной геометрической конфигурации при широком диапазоне физиологических нагрузок; комплексное междисциплинарное исследование закономерностей изменения напряженно-деформированного состояния зубочелюстной системы человека при использовании защитных зубных капп (контактное давление, контактное касательное напряжение, интенсивность напряжений в твердых тканях зубов и т.д.).

Целью данной работы является выявление закономерностей изменения напряженно-деформированного состояния биомеханического контактного узла при контакте пары зубов с учетом и без учета защитной протетической конструкции от уровня физиологических воздействий, по средствам численного моделирования. А так же анализ поведения материалов протетических конструкций при контактном взаимодействии с зубочелюстной системой человека в рамках решения модельных задач в осесимметричной постановке.

В работе рассмотрена задача численного моделирования контактного взаимодействия для двух расчетных схем (рис. 1): первая – контакт зубов верхнего (1) и нижнего (2) зубного ряда без учета протетической конструкции, вторая – контакт зубной пары через протетическую конструкцию (3). Физико-механические свойства Eva, материала из которого изготовлена зубная

Рис. 1. Расчетные схемы контакта элементов зубного ряда:

• а) с учетом и б) без учета протетической конструкции шина были получены по результатам серии экспериментов [11] (натурные эксперименты проводились научным коллективом трех ВУЗов г. Перми (ПГМУ, ПНИПУ, ПГНИУ) на базе ПГНИУ). Модуль Юнга для упругого участка составил 17,1±1,58 МПа, коэффициент Пуассона ^{ν =} 0^,46 , материал Eva проявляет упругопластические свойства, поэтому по результатам серии экспериментов была построена диаграмма σ-ε [20].

• 3.    РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ ЗАДАЧИ

Математическая постановка, приведенная в [13], дополняется статическими, кинематическими и контактными граничными условиями: на границе S σ приложена постоянная функциональная нагрузка, которая варьируется от 50 до 500 Н (сила индентирования) и запрещен изгиб; на границе S U запрещены перемещения по координате y; на границах контакта S _cont рассматривается контактное взаимодействие с учетом трения по сопрягаемым поверхностям и всеми типами контактного состояния (прилипание, проскальзывание, отлипание) с постоянным коэффициентом трения 0.3 [20].

Наибольший интерес представляют характер контакта в зонах смыкания элементов зубочелюстной системы S _K для первой расчетной схемы и S _{K 1} , S _{K 2} для второй расчетной схемы.

В рамках анализа сходимости результатов численного решения задачи контакта пары зубов через однослойную протетическю конструкцию проведена серия численных расчетов. Задача реализуется в программном комплексе ANSYS с использованием осесимметричных конечных элементов с лагранжевой аппроксимацией. По результатам серии расчетов установлено, что конечно-элементная сетка с градиентным сгу- щением элементов к зонам контакта, дает оптимальное по точности и времени счета решение задачи. Максимальный размер элемента 2,50E-04 м, минимальный – 6,25E-05 м.

Для анализа влияния учета протетической конструкции при взаимодействии элементов зубочелюстной системы человека проведена серия численных экспериментов, по результатам которой получены поля распеределения интенсивности напряжений в элементах зубного ряда (рис. 2-3). Результаты представлены для уровня нагрузки 300 Н.

При контакте пары зубов через упругопластическую однослойную каппу наблюдается значительное падение максимального уровня интенсивности напряжений: для зуба верхнего зубного ряда более чем на 60 %, для зуба нижнего зубного ряда более чем на 74%. При этом зона максимальной интенсивности напряжений не носит локальный характер, как при контакте без учета протетической конструкции.

На рис. 4 показаны зависимости максимальной интенсивности напряжений от силы инден-тирования для двух вариантов расчетных схем.

Максимальная интенсивность напряжений для расчетной схемы с учетом протетической конструкции из Eva ниже, чем при контакте зубов без учета каппы: для зуба верхнего зубного ряда в среднем в 2-2,5 раза меньше; для зуба нижнего зубного ряда – в 2-3,5 раза меньше. Интенсивность напряжений в протетической конструкции минимальная, что связанно с физико-механическими свойствами материала и проявлением упругопластического деформирования каппы. На рис. 5 показан качественная картина распределения интенсивности пластической деформации в зубной шине в зависимости от уровня функциональной нагрузки.

Рис. 5 отражает характер изменения интенсивности пластических деформаций при увеличении силы индентирования. Пластические деформации проявляются при нагрузке 100 Н,

МПа 0,65 2,00 3,34 4,68 6,03 7,37 8,71 10,1 11,4 11,8

Рис. 2. Интенсивность напряжений зуба верхнего зубного ряда: модель с учетом (а) и без учета (б) однослойной каппы

Рис. 4. Зависимость максимальной интенсивности напряжения от нагрузки: 1-2 нижний и верхний зуб соответственно при контакте без каппы; 3-4 – нижний и верхний зуб соответственно при контакте через каппы; 5 – каппа

Рис. 3. Интенсивность напряжений зуба нижнего зубного ряда: модель с учетом (а) и без учета (б) однослойной каппы

Рис. 5. Качественная картина распределения

пластических деформаций в каппе

при увеличении нагрузки вблизи области смыкания зубов нижнего зубного ряда с каппой. По мере увеличения нагрузки площадь пластических деформаций увеличивается и при нагрузке 300 Н появляется вблизи области смыкания зубов верхнего зубного ряда с каппой. Как и следовало ожидать максимум интенсивности пластических деформаций в каппе наблюдается вблизи зон контактного взаимодействия.

На рис. 6 показана зависимость максимальной интенсивности пластических деформаций от силы индентирования.

По мере увеличения силы индентирования максимальный уровень пластических деформаций возрастает по нелинейному закону и при максимальной нагрузке составляет приблизительно 28 %.

Следующим этапом исследования стал анализ влияния геометрии зуба нижнего зубного ряда на напряженно-деформированное состояние и зону контакта биомеханического узла. Рассмотрено три вида геометрической конфигурации зуба (рис. 7), отличающиеся размером первоначальной площадки контакта S _{K 2} : наибольшая площадка контакта у геометрии (б); наименьшая площадка контакта (изначальный контакт точечный) у геометрии (в). Обозначим

Рис. 6. Зависимость максимального уровня интенсивности пластических деформаций от нагрузки расчетные схемы с соответствующей геометрической конфигурацией: расчетная схема (а) (рис. 7, а); расчетная схема (б) (рис. 7, б); расчетная схема (в) (рис. 7, в).

На рис. 8-9 показан характер распределения полей интенсивности напряжения для расчетных моделей биомеханического контактного взаимодействия пары зубов через однослойную протетическую конструкцию отличающиеся геометрической конфигурацией зуба нижнего зубного ряда, результаты показаны для силы индентрования 500 Н.

Установлено, что изменение геометрии нижнего зуба влияет на интенсивность напряжений в элементах зубочелюстной системы. Как и следовало ожидать, при контакте через однослойную протетическую конструкцию во всех расчетных схемах происходит снижение максимального уровня интенсивности напряжений. При этом уровень снижения max σ _И отличается у всех трех рассмотренных расчетных схем и для зуба верхнего зубного ряда составляет 2.5, 1.6 и 40 раз для расчетных схем (а), (б) и (в) соответственно.

Аналогичная ситуация с уровнем снижения максимальных интенсивностей напряжения наблюдается и в элементах нижнего зубного ряда: понижение в 3.9, 2.7 и 40 раз для расчетных схем (а), (б) и (в) соответственно. Можно отметить, что максимальное снижение интенсивностей напряжений в зубной эмали при введении в расчетную схему однослойной протетической

Рис. 8. Интенсивности напряжений зуба верхнего зубного ряда: а), б), в) расчетные схемы а, б, в соответственно

Рис. 7. Геометрические конфигурации зубов нижнего зубного ряда

Рис. 9. Интенсивности напряжений зуба верхнего зубного ряда: а), б), в) расчетные схемы а, б, в соответственно

Рис. 10. Распределение контактного давления: а) зона контакта каппа-зуб верхнего зубного ряда; б) зона контакта каппа-зуб нижнего зубного ряда; 1, 2, 3 – расчетные схемы а, б, в соответственно

Рк , МПа                                     _{б )}

--J -Ц______________;___________I—i--------------;______________;______________ji—1—1_______; г , 1'

О 0,001 0,00 2 0,003 0,004 0,00 5 0,006 0 ,007

конструкции наблюдается в случае модели с изначально точечной зоной контакта (расчетная схема (в)).

Особый интерес так же представляет влияние геометрии элементов зубочелюстной системы человека на распределение контактного давления (рис. 10).

Можно отметить, что введение в биомеханический узел однослойной протетической конструкции привело к снижению контактного давления в зоне смыкания зубов к одному порядку. Уровень контактного давления при силе индентирования 500 Н во всех зонах контакта не превышает 7 МПа. В зоне контакта каппа-зуб верхнего зубного ряда максимальным контактным давлением обладает расчетная схема (б), зоне контакта каппа-зуб нижнего зубного ряда – расчетная схема (а). Стоит отметить, что уровень интенсивности пластических деформаций в каппе для расчетных схем (а), (б) практически одинаковый и не превышает 28 %. При одинаковых нагрузках уровень пластических деформаций расчетной схемы (в) выше чем у двух других расчетных схем приблизительно в 2,5 раза и при максимальной нагрузке 500 Н составляет 78 %.

• 4.    ЗАКЛЮЧЕНИЕ

В рамках исследования построена численная модель задачи контактного взаимодействия двух зубов верхнего и нижнего зубного ряда с учетом и без учета однослойной защитной протетической конструкции из Eva и выполнен анализ деформационного поведения биомеханического контактного узла в рамках теории деформационной упругопластичности при широком диапазоне функциональных нагрузок (50-500 Н). Выполнен анализ сходимости численного решения задачи при изменении степени дискретизации системы с выбором оптимальной конечно-элементной модели разбиения конструкции с градиентным уменьшением размера элемента к области контакта. Получено распределение напряжений и деформаций в паре зубов при контактном взаимодействии при разном уровне физиологической нагрузки для разных геометрических конфигураций элементов зубочелюстной системы.

При анализе результатов серии численных экспериментов установлено:

• -    При контакте пары зубов через протетиче-скую конструкцию максимальный уровень кон-

• тактного давления в среднем ниже в 3 раза, чем при контакте элементов зубочелюстной системы человека без учета протетической конструкции. Контактное касательное напряжение в 3-4 раз меньше контактного давления. Установлено влияние упругопластического деформирования каппы на характер распределения контактного давления и контактного касательного напряжения.

• -    При контактном взаимодействии пары зубов через протетическую конструкцию наблюдается падение максимальной интенсивности напряжений в твёрдых тканях зубов. В среднем максимальная интенсивность напряжения у зуба из нижнего зубного ряда при контакте через каппу меньше в 3,5 раза, чем при контакте без каппы, у зуба из верхнего зубного ряда – в 2,5 раза (для расчетных схем 1, 2).

• -    При увеличении нагрузки увеличивается площадь пластических деформаций в каппе, при этом при нагрузке 500 Н максимальный уровень интенсивности пластических деформаций не превышает 28 %. Влияние пластической деформации на зону контакта проявляется в характере распределения контактного давления и контактного касательного напряжения.

Рассмотрено влияние геометрических характеристик элементов зубочелюстной системы на деформационное поведение биомеханического узла. В рамках анализа рассмотрено три геометрии зуба нижнего зубного ряда с первоначально разной по площади площадкой контакта. Установлено:

• -    Снижение максимального уровня интенсивности напряжений в твердых тканых зубов наблюдается во всех ранее рассмотренных расчетных схемах. Минимальное снижение уровня σ _И наблюдается у расчетной схемы (б). Величина максимальной интенсивности напряжения расчетной схемы (в) без учета защитной протети-ческой конструкции значительно превышает значение этой же величины двух других расчетных схем, практически в 20 раз. При этом стоит отметить, что уровни интенсивности напряжения в расчетных схемах с использованием защитной протетической конструкции имеют количественные отличия, в среднем около 1 МПа.

• -    C увеличением силы индентирования, уровень пластических деформаций возрастает нелинейно и достигает максимального значения ~78% для расчетной схемы (в), и ~28 % для расчетных схем (а) и (б) для разной геометрии нижнего зуба.