УДК 531/534: [57+61]
КОНЕЧНО-ЭЛЕМЕНТНЫЙ АНАЛИЗ ЭФФЕКТИВНОСТИ УСТАНОВКИ ВКЛАДКИ ДЛЯ ЗАМЕЩЕНИЯ ДЕФЕКТА ПРЕМОЛЯРА
В.Ю. Кирюхин*, Г.И. Рогожников**, О.В. Шулятникова**
* Кафедра теоретической механики Пермского государственного технического университета, Россия, 614990, Пермь, Комсомольский проспект, 29, e-mail: [email protected]
** Кафедра ортопедической стоматологии Пермской государственной медицинской академии имени академика Е.А. Вагнера, Россия, 614990, Пермь, ул. Куйбышева, 39
Аннотация. В данной работе на основании конечно-элементного моделирования рассматривается работа вкладки, замещающей обширный дефект эмали зуба. Предполагается, что вкладка изготовлена из циркония и укреплена в зубе с помощью цемента. Анализируются два вида формы поверхности контакта вкладки зуба: плоская и пазообразная. Также предлагаются варианты нагружения вкладки (температурное и силовое) для сравнения и выделения наиболее опасного воздействия с позиции последующих возможных осложнений. В результате показано, что температурное воздействие имеет большее влияние на состояние вкладки и окружающие ее ткани, нежели силовое. Использование дополнительного паза в теле вкладки позволяет снизить пиковые напряжения в тех областях, в которых наиболее часто возникают осложнения.
Ключевые слова: восстановление зубов, вкладка, проектирование, цирконий, метод конечных элементов.
Введение
При возникновении больших дефектов жевательной поверхности зуба (сколы, обширный и глубокий кариес) в современной стоматологии применяют технологию установки вкладки. Суть данного метода состоит в том, что сначала удаляют причины поражения, трещины и пр. Затем пораженную область обрабатывают, придавая удобную для дальнейших манипуляций форму. После этого подготавливается фрагмент (вкладка) из какого-либо искусственного материала, который с помощью цемента устанавливается в подготовленную полость, восстанавливая тем самым полностью поверхность зуба. Вкладка может быть изготовлена по форме подготовленной полости, а также иметь стандартную форму, приданную изготовителем. Тогда полость готовится под имеющуюся вкладку.
Огромный опыт установки вкладок не избавил этот тип восстановления поверхности зуба от ряда негативных последствий: выпадение вкладки
(расцементирование), откол эмали зуба, возникновение кариеса на поверхности эмали в области ее контакта с материалом вкладки и пр. Тогда кажется очевидной актуальность анализа и последующего проектирования формы и материала вкладки, способа ее крепления с целью минимизировать вероятность появления приведенных выше нежелательных осложнений.
© Кирюхин В.Ю., Рогожников Г.И., Шулятникова G.B., 2007
09806267
Рассмотрение данной проблемы с позиций биомеханики началось достаточно давно. Одной из первых работ в этой области стало исследование по применению металлических вкладок для восстановления молочных зубов [1]. Причем практика многих лет продемонстрировала эффективность использования для изготовления вкладок таких металлов, как сплавы золота, титана [2-4]. Представленные в этих исследованиях результаты показали, что, в принципе, существенных различий между металлами с точки зрения механической эффективности при работе вкладки не существует. При этом такие материалы, как титан или золото, достаточно практичны при использовании в клинике наряду с полимеризационными композиционными материалами [5].
Вместе с тем, в последнее время отмечается огромный исследовательский интерес к применению материалов неметаллической природы [6]. Однако до сих пор в научной литературе не приведено доказательств полного преимущества в использовании того или иного типа материалов (металлов или композитов). Очевидно, те и другие имеют свои положительные качества, которые принимает во внимание практикующий врач при принятии решения.
В последнее время в ортопедической стоматологии стал активно использоваться цирконий. Обладая хорошими свойствами в отношении биологической инертности, он имеет высокие механические свойства (см. далее таблицу) и удобнее в клинической практике [5]. Поэтому авторы данной работы предлагают обратиться к анализу работы вкладки, изготовленной из такого материала.
Дополнительной особенностью исследования практических задач стоматологии с позиции биомеханики является выделение фактора, воздействующего на объект анализа: силовой и/или температурный. В действительности, до сих пор остается не выясненным до конца вопрос, какой из них оказывает наибольшее влияние на запуск и протекание процессов разрушения. Так, многочисленные исследователи, к примеру [7-9], рассчитывают поведение зубочелюстной модели под воздействием силового нагружения, вызванного жеванием. Намного скромнее список работ, к примеру [10], в которых рассматривается температурное воздействие. Однако на основании результатов упомянутых исследований не следует, что тепловое или силовое воздействие играет преобладающую роль, на которую нужно ориентироваться при проектировании того или иного метода лечения какого-то поражения или заболевания в области зубов и окружающих их тканей. Правда, следует отметить, что «чистое» силовое воздействие рассматривается в таких работах чаще.
Конечно, учет обоих факторов нагружения приводит к наиболее адекватным результатам анализа состояния тканей зубочелюстной системы. Однако в этих условиях попытка спроектировать операцию или параметры применяемых при лечении средств (материалы, размеры области, способ фиксации и пр.) оборачивается в лучшем случае возрастающей сложностью. Но и не исключена ситуация, когда учет факторов температуры и силы приводит к противоречивым тенденциям при проектировании. Например, использование твердых сплавов металлов для замещения дефектов эмали зуба оправдано с точки зрения деформационной совместимости и прочности в условиях силового нагружения. Однако давно показано в теории [11] и на практике [5], что такая комбинация материалов приводит к резким негативным последствиям из-за температурной несовместности титана и эмали, выраженной в большом различии значений коэффициентов линейного температурного расширения. И наоборот, применение материалов, близких по теплопроводным свойствам к соседствующим живым тканям (стекловолоконные материалы, к примеру), имеет тенденцию к более низкому качеству по реакции на интенсивное или долговременное силовое нагружение.
Поэтому авторы данной работы считают актуальным сравнить, воздействие какого из факторов (температура или сила) соответствует картине осложнений, наиболее часто наблюдаемых на практике, в долгосрочной [12] или краткосрочной перспективе [13].
Цели работы
В данной работе на основе механического анализа предлагается исследовать основные причины, приводящие к негативным последствиям установки вкладки. А именно: какой из механических факторов (сила или температура) оказывает
наибольшее негативное воздействие на состояние эмали, и как влияет на качество установки вкладки использование паза.
Поэтому в рамках настоящего механического исследования авторы ставят целью следующие вопросы:
1) сравнить напряженное состояние целого зуба как механической конструкции под действием жевательной нагрузки и под действием температуры,
2) провести аналогичное сравнение для восстановленного вкладкой зуба без паза и с пазом.
Для решения этой проблемы предлагается следующая последовательность действий.
1. Создать трехмерную геометрическую модель расчетной схемы моляра по выбранному образцу. Задать размеры и расположение циркониевой вкладки. Задать тип цементирующего материала и его толщину между зубом и вкладкой.
2. Смоделировать на основе теории термоупругости процесс силового нагружения и нагрева системы зуб-вкладка. Использовать при этом численный метод конечных элементов как хорошо зарекомендовавший себя в решении задач биомеханики в стоматологии.
3. Провести анализ результатов расчета, а именно:
• продемонстрировать напряженно-деформированное состояние зуба и отдельных его компонентов при указанных вариантах нагрузки;
• провести качественное сравнение различных показателей напряженного состояния и выделить отличия в силовом или температурном нагружении, при использовании паза или без него;
• указать наиболее опасные зоны (по величине напряжений и деформаций).
Модель зуба как биомеханической системы
Структура расчетной схемы
Для расчета предлагается модель, полученная по реальным размерам премоляра (двухкорневого зуба). В рассматриваемую расчетную схему включены следующие элементы: зуб, сама вкладка и фиксирующий ее цемент. Зуб представлен в виде конструкции из эмали и дентина, в которых подготовлены полости для установки вкладки, см. рис. 1. На рис. 2 показана сборка модели с условно разнесенными частями. Формы цемента и вкладки проиллюстрированы на данных диаграммах с подготовленным пазом, который по предположению авторов, с одной стороны, должен улучшить качество фиксации вкладки. С другой стороны, он призван увеличить
скорость теплообмена между вкладкой и сохранившимися тканями зуба. Никакие иные компоненты (альвеолярный отросток, нижняя челюсть, зубы-антагонисты) зубочелюстной системы в данной работе не рассматриваются, поскольку они либо достаточно удалены от области интересов исследования (вкладки), либо моделируются в виде соответствующих граничных условий.
Дополнительной особенностью моделирования дентина в этой работе является введение полого канала для имитации механической работы пульпы.
Вкладка Паз во вкладке Цемент Эмаль
Канал для пульпы Дентин
Рис. 1. Сечение зуба с пазом в установленной вкладке
Пищевой комок Вкладка Цемент Эмаль Дентин
Рис. 2. Схема расчетной области с разнесенными деталями.
Здесь указана модель с пазом во вкладке
Модель в рамках теории упругости
Все расчеты по определении и анализу напряженного состояния зуба с вкладкой (и без нее) проводятся на основе линеаризованной теории упругости. Обоснование применимости данной теории к задачам такого рода подробно обсуждается в работе [14]. Ниже приводятся ее основные соотношения в классической постановке.
Пусть исследуемое тело занимает ограниченную область О трехмерного евклидова пространства Е3. Замыкание области обозначим через О, границу (которая считается достаточно гладкой) - через Г ( О = О ^ Г ).
1) Уравнение статического равновесия внутри области
Біу 5 + 0 = 0, х є О, (1)
где а - симметричный тензор напряжений, 0 - вектор объемной силы, 0 є (С(О)) . В формуле (1) и далее величины 5, , и считаются функциями координат,
представленными радиусом-вектором г є О .
2) Деформации будем считать малыми и аддитивными, т.е. тензор малой деформации представляется в виде суммы тензоров упругой е и температурной т деформации
е + т, г є О . (2)
3) Упругие деформации связаны с напряжениями законом Гука
5 = С • • е, г є О, (3)
где С - четырехвалентный тензор модулей упругости, Сум є С1(О) .
4) Соотношения деформация-перемещения записываются в рамках линеаризованной теории
5(и) = 2(У и + и V), г є О, (4)
^ ^ 2 — 3
где и - вектор перемещения, и є (С (О)) .
Заметим, что компоненты тензора деформации удовлетворяют условиям
совместности деформаций, которые эквивалентны обращению в нуль компонент
тензора второго ранга го!:(го1 ) (при условии существования вторых производных от
компонент деформации по координатам):
го!:(го1 в) = 0, г є О. (5)
5) Будем считать, что граница области Г делится на две взаимно непересекающиеся части Г = Ги + Г5 . На части границы Ги заданы кинематические граничные условия, на
части Г5 задан вектор напряжений Р є (С(Г5 ))3 :
и ^ г єГи , (6)
п •а=Р, г єГа. (7)
6) Будем считать, что температурные деформации т є (С1(О))6 будут равны
т
5т =
|ат (т) ёт , г є О, (8)
то
где ат (Т) - тензор коэффициентов линейного теплового расширения как функция температуры Т; Т - температура, при которой тело находится в естественном (ненапряженном и недеформированном) состоянии.
Г раничные условия
Прежде чем перейти к описанию граничных условий задачи, необходимо отметить предположение, которое, строго говоря, к ним не относится. При расчете напряженно-деформированного состояния авторы предлагают пренебрегать воздействием объемных сил на компоненты модели: Q = 0. Обоснование данного предположения базируется на следующих аргументах. Среди возможных источников сил Q можно указать лишь силу тяжести, порядок величины которой составляет 10 гр • 9,8 м/с2 ~ 0,098 Н. Данное значение несоизмеримо мало по сравнению с нагрузками, прикладываемыми к модели.
В механических задачах вообще и, особенно, в задачах биомеханики важной частью проведения расчетов в соответствии с выбранной моделью является конкретизация граничных условий, которые необходимы для того, чтобы передать на исследуемое тело силовую или температурную нагрузку.
Традиционно выбор граничных условий в задачах биомеханики требует анализа и предположений, поскольку приходится вычленять интересующую область из более общей биомеханической системы, обладающей в целом сложной геометрией и нетривиальным движением. В этих условиях не всегда возможно заранее, не решая задачи, определить точки (области) с известными перемещением или напряжением. Зачастую этот вопрос исследователи решают интуитивно или на основании эксперимента.
Одним из принципиальных решений данной проблемы является отдаление границ от рассматриваемой области. Тогда согласно принципу Сен-Венана способ задания граничного условия будет в меньшей степени влиять на результат в интересующей исследователя области.
Рис. 3. Закрепленные области (отмечены темным Рис. 4. Темным тоном обозначены области тоном) и жевательные нагрузки приложения температурных нагрузок: сверху -
нагрев от жидкости температурой 60?С, снизу -каналы пульпы с температурой 36,6?С
В данной работе авторы предлагают во всех моделях считать неподвижными наружные поверхности корней зуба. Температура считается постоянной и равной 36,6?С по каналу пульпы в массе дентина, рис. 3, 4.
В отношении моделирования нагрузки на коронку зуба авторы предлагают новую идею. Дело в том, что в моделях прежних исследователей нагрузки прикладывались непосредственно к коронке, то есть делались малообоснованные гипотезы от том, каким образом та или иная сила жевания распределяется по жевательной поверхности. Поэтому в данной работе предлагается ввести в рассмотрение идеализированный фрагмент пищи или употребляемой пациентом жидкости, на который и действует необходимая нагрузка, см. рис. 4.
Другой важной особенностью является постановка в соответствии друг другу силовых и температурных граничных условий, поскольку стоит задача сравнения их влияния на зуб. Конечно, приуменьшив одно воздействие и преувеличив другое, можно добиться практически любых результатов. Ведь сравниваемое температурное воздействие никак не может быть связано с предполагаемым усилием жевания. Поэтому авторы предположительно выбрали жевательное усилие, равное 200 Н на рассматриваемый зуб, и температуру 60?С. Эти данные представлены как усредненные. Но в этом случае, как показали результаты, заметно отличие температурного и силового воздействий.
Свойства
Проблема определения свойств для элементов биомеханической конструкции в данном случае возникает не из-за того, что нет соответствующих данных. Как раз наоборот, они очень разнообразны. Так, к примеру, различные марки цемента имеют существенные отличия по физико-механическим свойствам. То же самое можно сказать и про цирконий, используемый как достаточно новый для стоматологии материал.
Проблема другого рода возникает при попытке зафиксировать температурные и упругие свойства пищевого комка.
Однако авторы видят выход из возникшей неопределенности в следующем. Основным инструментом анализа данной работы является сравнение различных вариантов геометрии и нагрузки. Поэтому достаточно выбрать некоторые усредненные характеристики для любого фрагмента модели, близкие к возможным реальным значениям. Все они сведены в следующей таблице.
Таблица
Механические свойства материалов, используемых в модели________________________
Материал Модуль упругости, ГПа Предел прочности на растяжение, МПа Предел прочности на сжатие, МПа Коеффициент теплопроводности, 10-6 м2/с Коэффициент линейного температурного расширения, 10-6/?С
Цирконий 99 30 350 150,0 11,0
Эмаль 48 21 270 0,47 12,0
Дентин 18 36 260 0,18 7,5
Цемент 23 - 110 0,60 7,5
Результаты
Описанная выше расчетная схема со всеми указанными значениями материальных свойств и вариантов граничных условий была смоделирована численно с помощью метода конечных элементов как пространственная структура. И ниже авторы представляют иллюстрации с полученными данными, сопровождая их необходимыми комментариями.
Для иллюстрации и оценки результатов могут использоваться различные параметры напряженно-деформированного состояния. Авторы выбрали первый инвариант тензора напряжений:
а = а + а2 + а3, (9)
напряжения по Мизесу
и максимальное нормальное напряжение
а„ =аь (11)
где аь а 2 , а3 - главные напряжения.
Здоровый зуб
Для того, чтобы оценить эффект, создаваемый в зубе со вкладкой, температурным или силовым воздействием, авторы сначала рассматривают полноценный здоровый зуб, имеющий целые эмаль и дентин. На рис. 5-7 показаны распределения на поверхности первого инварианта тензора напряжений, напряжения по Мизесу и запаса прочности в соответствии с максимальным нормальным напряжением.
Очевидно, распределение первого инварианта тензора напряжений, напряжения по Мизесу, максимальных нормальных усилий (а также иных факторов, здесь не указанных) при температурной нагрузке менее благоприятно, поскольку большие значения в случае силовой нагрузки вызваны лишь краевыми эффектами. Можно отметить, насколько ниже запас прочности по критерию максимального нормального напряжения в случае температурного нагружения.
а , Па а , Па
4,229е+007 И 3.877е+007 . 3,525е+007 .3.172е+007 . 2,820е+007 .2.467е+007 2.115е+007 . 1 762е+007 . 1.410е+007 .1 057е+007 1.7,050е+006
^3.526е+006 1.213е+003
3.683е+007 1.3.376е+007 .3.06Эе+007 .2.763е+007 .2.45Ве+007 .2.149е+007 . 1 842е+007 . 1 535е+007 . 1 229е+007 .9.218е+00В .6.150е+006 1.3.082е+006 1.1 428е+004
Силовое нагружение Температурное нагружение
Рис. 5. Первый инвариант тензора напряжений а
а
, Па
УОМ ■
■ 3.789Є+007 К 3.473е+007 . 3.158е+007 . 2.842е+007 . 2.526Є+007 . 2.210е+007 . 1.395Є+007 і 1.579Є+007 11,263е+007 • 9.473е+006 ^(.6.316е+006 Н.3.159Є+006 Н. 1.115Є+003
а
УОЫ
, Па
I
3.198Є+007 2.931 е+007 2.665е+007 2.399е+007 2.132е+007 1 .Зббе+007 1.600е+007 1 .ЗЗЗе+007 1.0676+007 8.004Є+006 5.340Є+006 2.676Є+006 1,266е+004
Силовое нагружение
Температурное нагружение
Рис. 6. Напряжение по Мизесу а
УОЫ
Имя модели: Полный зуб
Имя упражнения: Силовое нагружение
Тип эпюры: Проверка проектирования Построениеі
Критерий: Максимальное нормальное усилие
Распределение запаса прочности: Мин. коэффициент запаса прочности = 2.6
Силовое нагружение
Имя модели: Полный зуб Имя упражнения: Температурное нагружение Тип эпюры: Проверка проектирования Построение4 Критерий: Максимальное нормальное усилие
Распределение запаса прочности: Мин. коэффициент запаса прочности = 0.45
1.9.186
■+00
|8.376е+00 . 7.564е+00 . 6.752е+00 -5.940Є+ОО 5.128е+00 14.316е+00 .3.504Є+00 . 2.692е+00 . 1.880Є+00 І1Л68Є+00' 12.563е+000
а„ , Па
Температурное нагружение Рис. 7. Оценка запаса прочности по а
1 .ОООе+002 19.170е+001 Щ 8.341 е+001 . 7.511е+001 . 6.682е+001 , 5.852Є+001 І5.023Є+001 .4.193е+001 .3.364Є+001 .2.534Є+001 і 1.704е+001 .8.749Є+000 И.4.536Є-001
Восстановленный зуб
Теперь рассматривается зуб с установленной в него вкладкой из циркония. Причем расчеты производятся для двух вариантов формы вкладки. В первом случае вкладка имеет форму замещаемой области. Тогда слой цемента между вкладкой и тканями зуба имеет форму, показанную на рис. 8а. Во втором случае во вкладке дополнительно создается полость, заполняемая при восстановлении зуба цементом. Такой технологический прием приводит к созданию объема цемента, имеющего форму, как на рис. 8б. Для того, чтобы лучше сориентироваться, читатель отсылается к рис. 1 и
2, где область цемента из рис. 8б показана в общей композиции. На рис. 9 показана расчетная область дентина и эмали для обоих типов вкладки. Видно, что при выполнении операции врачу приходится препарировать как эмаль, так и дентин в относительно больших объемах.
Дальнейшие решения искались для двух задач на одной и той же модели области: задачи теплопроводности и задачи термоупругости при сформулированных граничных условиях. Так, на рис. 10 проиллюстрировано, как распределяется температура по объему цемента при квазистатическом температурном нагружении для обоих форм вкладки. Следует отметить, что в случае использования дополнительного паза температура распределяется ровнее. На диаграмме (рис. 10б) визуально вообще невозможно выделить области с температурой, отличной от средней по объему цемента. В отношении же остальных элементов модели (рис. 11) эффект от применения паза при решении задачи теплопроводности вообще не заметен.
Рис. 8. Форма цемента для рассматриваемых случаев подготавливаемой формы вкладки
Рис. 9. Схема подготовки полости для установки вкладки с цементом в виде прослойки
Имя модели: Зуб без паза во вкладке Имя упражнения: Теппоообмен Тип эпюры: Термический Построение1 Временной шаг 1
Т
т'
, ?С
І.ОООе+ОСГ
;.805е+00
;.бЮе+оо
;.415е+00
;.220е+00
;.025е+00
-.8300+00
•.635Є+00
•.440Є+00
•.245Є+00
•.050Є+00
;.855е+00
|.660е+00
Имя модели: Зуб + паз во вкладке Имя упражнения: Теплообмен Тип эпюры: Термический Построениеі Временной шаг 1
_.6.000е+00
5.805е+00 .5.610Є+00 .5.415е+00 .5.220Є+00 . 5.025е+00 .4.830е+00 .4.635Є+00 . 4.440е+00 .4.245Є+00 4.050е+00
13.855е+00 З.ббОе+ОО
Рис. 10. Распределение температуры по объему цемента для вкладки без паза (а) и с пазом (б)
Имя модели: Зуб без паза во вкладке Имя упражнения: Теплоообмен Тип эпюры: Термический Построениеі Временной шаг 1
6.000Є+001
™ 5.805©+001 .5.610Є+001 .5.415Є+001 .5.220©+001 .5.0250+001 . 4.830е+001 .4.635Є+001 .4.440Є+001 .4.2450+001 ^.4.050©+001 1.3.8550+001 1.3.6600+001
Имя модели: Зуб + паз во вкладке
6.000Є+001
\ 5.805Є+001 . 5.610е+001 . 5.415е+001 . 5.220е+001 . 5.025е+001 . 4.830е+001 . 4.635е+001 . 4.440е+001 . 4.245Є+001 ,4.050е+001
[3.855е+001 3.6Є0Є+001
Рис. 11. Распределение температуры по объему сохранившихся тканей зуба (рис. 9) применении вкладки без паза (а) и с пазом (б)
Имя модели: Зуб без паза во вкладке Имя упражнения: Сило бое нагружение Тип эпюры: Проверка проектирования Построениеі Критерий: Максимальное усилие уоп Мібєб
Распределение запаса прочности: Мин. коэффициент запаса прочности = 8.8
У
я:
Имя модели: Зуб без паза во вкладке Имя упражнения: Температурное нагружение Тип эпюры: Проверка проектирования Построениеі Критерий : Максимальное усилие уоп Мібєб
Распределение запаса прочности: Мин. коэффициент запаса прочности = 9.1
.000е+002 ,240е+00' .480Є+00 720е+00 ,96Ое+0О .2ООЄ+0О 440е+00 ,68Ое+0О .920Є+00 .160е+00 .4ООЄ+0О .640Є+00' .804е+ООО
при
У
і
1.000е+002 9.242е+00-8.484Є+00 . 7.726е+00
-6.21ІЄ+00
. 5.453е+00 _4.695е+0О . 3.937е+00 .3.179Є+00 . 2.422е+00 1.664Є+00' Я9.060Є+000
а б
Рис. 12. Оценка запаса прочности по нормальному напряжению для вкладки без паза при силовом (а) и температурном (б) нагружениях
а
Имя модели: Зуб без паза во вкладке Имя упражнения: Силовое нагружение Тип эпюры: Проверка проектирования ПостроениеЗ Критерий : Максимальное касательное усилие
Распределение запаса прочности: Мин. коэффициент запаса прочности = 1.2
У
I
Имя модели: Зуб без паза во вкладке Имя упражнения: Температурное нагружение Тип эпюры: Проверка проектирования ПостроениеЗ Критерий: Максимальное касательное усилие
Распределение запаса прочности: Мин. коэффициент запаса прочности = 0.61
000е+002 177е+00' 354е+00 .531 е+00 708е+00 884е+00 .061 е+00 238е+00 .415е+00 592е+00 769е+00' 457е+000 226е+000
У
Я
1.000е+002 9.172е+001 |.8.343е+001 . 7.515е+001 . 6.687е+001 . 5.859е+001 .5.030е+001 .4.202е+001 . 3.374е+001 .2.545е+001 . 1.717е+001 а.8.888е+000
а б
Рис. 13. Оценка запаса прочности по максимальному касательному усилию для вкладки без паза при силовом (а) и температурном (б) нагружениях
Имя модели: Зуб без паза во вкладке Имя упражнения: Теппоообмен Тип эпюры: Термический Построение1 Временной шаг 1
6.000е+001 |.5.805е+001 . 5.610е+001 . 5.415е+001 .5.220е+001 . 5.025е+001 . 4.830е+001 .4.635е+001 .. 4.440е+001 .4.245е+001 ^.4.050е+001 1.3.855е+001 Ёз.660е+001
Рис. 14. Распределение температуры по объему цемента для вкладки без паза
Имя модели: Зуб + паз во вкладке Имя упражнения: Силовое нагружение Тип эпюры: Проверка проектирования Построение1 Критерий: Максимальное усилие уоп М1эе5
Распределение запаса прочности: Мин. коэффициент запаса прочности =
Имя модели: Зуб + паз во вкладке Имя упражнения: Температурное нагружение Тип эпюры: Проверка проектирования Построение1 Критерий : Максимальное усилие уоп М1зеэ Распределение запаса прочности: Мин. коэффициент зе
^1 .ОООе+С 9.308е+( Щ. 8.616е+( .7.924е+1 .7.233е+( , 6.541 е+( .5.349е+[ -.5.157е+С .4.465е+[ _3.773е+[ . 3.081 е+С
У
К
.000е+002 9.352е+001 18.703е+001 _8.055е+001 .7.406е+001 .6.758е+001
|.6.110е+001 5.461 е+001 . 4.81 Зе+001 .4.164 е+001 . 3.516е+001
.2.81
е+001
а б
Рис. 15. Оценка запаса прочности по нормальному напряжению для вкладки с пазом при силовом (а) и температурном (б) нагружениях
Имя модели:.Зуб + паз во вкладке
Имя упражнения: Силовое нагружение
Тип эпюры: Проверка проектирования ПостроениеЗ
Критерий: Максимальное касательное усилие
Распределение запаса прочности: Мин. коэффициент запаса прочности =
Имя модели: Зуб + паз во вкладке
Имя упражнения: Температурное нагружение
Тип эпюры: Проверка проектирования ПостроениеЗ
Критерий: Максимальное касательное усилие
Распределение запаса прочности: Мин. коэффициент запаса прочности =
I
.000е+[ 9.177е+С .355Є+С 7.532е+С 6.71 Ое+С .687Є+С 5.065е+С ,242е+С З.420е+С 2.597е+С 1,775е+С 9.522е+С 1.297Є+С
Д00е+002 9.175е+001 .349е+001 _7.524е+001 .6.699Є+001 ;л1.5.874е+001 -5.043е+001 _4.223е+001 _3.398е+001 . 2.57Зе+001 .1.7470+001
Ё- 9.222Є+000 9.697Є-001
У
У
а б
Рис. 16. Оценка запаса прочности по максимальному касательному усилию для вкладки с пазом
при силовом (а) и температурном (б) нагружениях
Рис. 17. Распределение температуры по объему цемента для вкладки без паза: а - вид снаружи, б
- вид изнутри без дентина
Рассмотрим теперь отдельно параметры напряженно-деформированного состояния вкладки без паза и с пазом и сравним ее работу при силовом и температурном нагружениях. Результаты представлены на рис. 12-17.
Выводы
На основании проделанных расчетов и полученных числовых значений можно со всей определенностью сделать следующие выводы.
1. Критерий максимального нормального напряжения дает результат, при котором эмаль зуба должна начать разрушаться уже при употреблении напитков или пищи с температурой около 60?С. Однако этого не происходит, и объяснение этому авторы видят в том, что такой критерий плохо применим для таких существенно неизотропных материалов, как эмаль.
2. В данной работе предложено новое решение проблемы приложения нагрузки к зубу с использованием пищевого комка. Такой подход избавил авторов от необходимости определять форму и распределение усилий по поверхности эмали и вкладки зуба. Авторы считают, что данная идея может быть применена и при анализе иных биомеханических конструкций, позволяя решить нетривиальную проблему задания граничных условий.
3. Еще одним значимым отличием работы от известных ранее стало использование температурной нагрузки в качестве основного воздействия на зуб. С помощью
расчетов было показано, что температурное воздействие на зуб приводит к более негативным последствиям, нежели силовое нагружение.
4. Применение паза во вкладке приводит к улучшению теплообмена и снижению напряжений в сильно нагруженных областях на величины порядка 10%. Такой вывод можно сделать на основании численных результатов, проиллюстрированных на рис. 11-17.
Список литературы
1. Ireland, R.L. Low fusing metal inlays for fillings for deciduous teeth / R.L. Ireland // International Journal of Orthodontia and Dentistry for Children. - 1935. - Vol. 21, No. 6. - P. 577-580.
2. Pascal, M. Stress distribution of inlay-anchored adhesive fixed partial dentures: A finite element analysis of the influence of restorative materials and abutment preparation design / M. Pascal, P. Nikolaos, C.B. Urs, K. Ivo // The Journal of Prosthetic Dentistry. - 2002. - Vol. 87, No. 5. - P. 516-528.
3. Wolf, B.H. Margin quality of titanium and high-gold inlays and onlays - a clinical study / B.H. Wolf, M.H. Walter, K.W. Boening, A.E. Schmidt // Dental Materials. - 1998. - Vol. 14, No. 5. - P. 370-374.
4. Walter, M.H. Plaque, gingival health and post-operative sensitivity in titanium inlays and onlays: a randomized controlled clinical trial / M.H. Walter, B.H. Wolf, A.E. Schmidt, K.W. Boening, R. Koch // Journal of Dentistry. - 2001. - Vol. 29, No. 3. - P. 181-186.
5. Рогожников, Г.И. Сплавы титана в ортопедической стоматологии / Г.И. Рогожников, В.А. Четвертных, М.Д. Кацнельсон, Н.Б. Асташина. - Пермь: Перм. гос. мед. академия, 2007.
6. Sadowsky, S.J. An overview of treatment considerations for esthetic restorations: A review of the literature/ S.J. Sadowsky // The Journal of Prosthetic Dentistry. - 2006. - Vol. 96, No. 6. - P. 433-442.
7. Hubsch, P.F. A finite element analysis of the stress at the restoration-tooth interface, comparing inlays and
bulk fillings / P.F. Hubsch, J. Middleton, J. Knox // Biomaterials. - 2000. - Vol. 21, No. 10. -
P. 1015-1019.
8. Ratledge, D.K. The effect of restorative materials on the wear of human enamel / D.K. Ratledge, B.G.N. Smith, R.F. Wilson // The Journal of Prosthetic Dentistry. - 1994. - Vol. 72, No. 2. - P. 194-203.
9. Derand, T. Stress analysis of cemented or resin-bonded loaded porcelain inlays / T. Derand // Dental Materials. - 1991. - Vol. 7, No. 1. - P. 21-24.
10. Goehring, T.N. Marginal adaptation of inlay-retained adhesive fixed partial dentures after mechanical and thermal stress: An in vitro study / T.N. Goehring, O.A. Peters, F. Lutz // The Journal of Prosthetic Dentistry.
- 2001. - Vol. 86, No. 1. - P. 81-92.
11. Pechenov, V.S. Finite Element Stress Analysis of Tooth Crowns with Titanium Inlays / V.S. Pechenov,
N.B. Astashina, Yu.I. Nyashin, G.I. Rogozhnikov // Russian Journal of Biomechanics. - 1997. - Vol. 1,
No. 1-2. - P. 108-118.
12. Malament, K.A. Survival of Dicor glass-ceramic dental restorations over 14 years. Part II: Effect of thickness of Dicor material and design of tooth preparation / K.A. Malament, S.S. Socransky // The Journal of Prosthetic Dentistry. - 1999. - Vol. 81, No. 6. - P. 662-667.
13. Gohring, T.N. Clinical and scanning electron microscopic evaluation of fiber-reinforced inlay fixed partial dentures: Preliminary results after one year / T.N. Gohring, W.H. Mormann, F. Lutz // The Journal of Prosthetic Dentistry. - 1999. - Vol. 82, No. 6. - P. 662-668.
14. Кирюхин, В.Ю. Конечно-элементный анализ механических причин возникновения вторичных деформаций / В.Ю. Кирюхин, Г.И. Рогожников, В.А. Вершинин // Российский журнал биомеханики.
- 2005. - Том 9, № 3. - С. 16-31.
FINITE ELEMENT ANALYSIS OF THE EFFECTIVENESS OF INLAY INSTALLATION FOR PREMOLAR DEFECT TREATMENT
V.Y. Kiryukhin, G.I. Rogozhnikov, O.V. Shulyatnikova (Perm, Russia)
In the present paper, the functioning of the inlay that replaces spacious defect of the tooth enamel is considered. It is assumed that inlay is made of zirconium and is set in tooth by means of cement. Two kinds of surface shapes of the tooth inlay (without and with groove) are analyzed. Also, for comparison and for searching the most dangerous action from the position of the possible complications the variants of inlay loading by temperature and by masticatory forces are considered. As a result, it is shown that thermal action influences the state of inlay and surrounded tissues more effectively than the force loading. The utilization of subsidiary groove reduces the peak stresses in regions where the complications are localized more frequently.
Key words: tooth restoration, inlay, design, zirconium, finite element method.
Получено 15 марта 2007