CC BY
40
DOI: 10.15593/RZhBiomeh/2018.3.04 УДК [53:57:61+004](082)
МОДЕЛИРОВАНИЕ ДИНАМИЧЕСКОГО ПОВЕДЕНИЯ МОСТОВИДНОГО ЗУБНОГО ПРОТЕЗА С ПОМОЩЬЮ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ
А.Е. Крупнин1,2, Я.Н. Харах3, Л.Г Киракосян3, С.Д. Арутюнов3
1 Национальный исследовательский центр «Курчатовский институт», Россия, 123098, Москва, пл. Академика Курчатова, 1, e-mail: [email protected]
2 Московский государственный технический университет имени Н.Э. Баумана, Россия, 105005, Москва, ул. 2-я Бауманская, 5
Кафедра пропедевтики стоматологических заболеваний Московского государственного медико-стоматологического университета им. А.И. Евдокимова, Россия, 127473, Москва, ул. Делегатская, 20, стр. 1, e-mail: [email protected]
Аннотация. Цель исследования - применение методов динамического анализа в задачах биомеханики зубочелюстного аппарата на примере мостовидного зубного протеза из полиметилметакрилата. Задача решена в плоско-деформированной постановке. Проведена экспериментальная верификация на основании данных эксперимента для керамического протеза идентичной конфигурации. На основании результатов модального анализа определены собственные частоты и константы для модели пропорционального демпфирования Рэлея. Исследовано влияние длительности действия силы на напряженно-деформированное состояние биомеханической системы. Определены области наибольших эквивалентных (по Мизесу) напряжений. Установлено, что при стремлении длительности действия возбуждающей силы к периоду собственных колебаний системы амплитудные значения напряжений на окклюзионной поверхности протеза могут превышать предел текучести материала. При увеличении длительности действия возбуждающей силы значения динамических и статических напряжений не различаются. Показано, что влияние пульпы на результаты статического расчета незначительно, в то время как учет пульпы приводит к более интенсивному поглощению энергии, затуханию колебаний и снижению амплитудных значений перемещений и напряжений.
Ключевые слова: биомеханика зубочелюстного аппарата, мостовидный зубной протез, метод конечных элементов, динамика, модальный анализ, демпфирование.
Введение
Исследования последних лет, основанные на данных Всемирной организации здравоохранения и отечественных авторов, свидетельствуют о высоком уровне распространенности частичного отсутствия зубов и динамике его роста [3, 5, 17-19, 25, 27, 39, 43]. Основной причиной потери зубов являются кариес и его осложнения. Но замещение дефектов протезами ведет к деформации зубных рядов, нарушению их окклюзионных взаимоотношений, перепрограммированию функций мышц и височно-
© Крупнин А.Е., Харах Я.Н., Киракосян Л.Г, Арутюнов С. Д., 2018
Крупнин Артур Евгеньевич, инженер-исследователь лаборатории полимерных материалов; ассистент кафедры прикладной механики, Москва
Харах Ясер Насерович, ассистент кафедры пропедевтики стоматологических заболеваний, Москва Киракосян Левон Гамлетович, старший лаборант кафедры пропедевтики стоматологических заболеваний, Москва
Арутюнов Сергей Дарчоевич, д.мед.н., профессор, заведующий кафедрой пропедевтики стоматологических заболеваний, Москва
нижнечелюстного сустава, формированию трудноизлечимых, а порой неизлечимых заболеваний зубочелюстного аппарата.
В стоматологической практике несъемные зубные протезы являются самыми распространенными ортопедическими конструкциями для замещения включенных дефектов зубных рядов [13]. Однако при их использовании достаточно часто возникают поломки [33, 36], обусловленные врачебными ошибками прогнозирования результатов лечения из-за необоснованного выбора опорных зубов, конструкции протеза, технологии изготовления и конструкционного материала. Для снижения риска разрушения и прогнозирования влияния установленного протеза на твердые и мягкие ткани широко применяются численные методы, в частности, метод конечных элементов [6, 22, 23, 30, 35, 36, 38, 40, 46, 48].
Метод конечных элементов позволяет оценить максимальные нагрузки, которые выдерживает изделие, напряжения и деформации, возникающие в протезе и окружающих тканях. Такой подход дает возможность еще на стадии проектирования выбрать наилучший вариант конструктивного исполнения, соответствующий биомеханическим требованиям, сократив тем самым ресурсы и предотвратив возможные усложнения после имплантации.
Результат численного исследования напряженно-деформированного состояния конструкций зависит от многих аспектов. Важными шагами являются редукция исследуемого нативного объекта до адекватного модельного представления и интерпретация нагрузок, действующих на биомеханическую систему. В первом случае речь идет о выборе геометрии исследуемой области и представлении механических свойств материалов, из которых эта система состоит. Многие задачи биомеханики зубочелюстного аппарата успешно решаются в двумерной (плоской) постановке [8-10, 21, 22, 32, 35, 38, 40].
Несмотря на определенные допущения и ограничения, накладываемые сведением трехмерных систем к двумерным, решение плоских задач имеет ряд преимуществ [30]: качественная густая сетка, уменьшение времени решения задачи, возможность получения быстрых результатов на нескольких моделях. Это особенно ценно в задачах, где требуется оценить влияние одного или нескольких параметров на другие или на отклик системы в целом. Выбор двумерной схемы хорошо согласуется с экспериментальными результатами [22], а полученные решения плоских задач соответствуют решениям, аналогичным трехмерным в исследуемом сечении. Решение задач в трехмерной постановке дает более полную картину деформирования системы, но априори имеет большую размерность и требует значительных вычислительных мощностей для получения адекватного интерпретируемого результата [30].
Сложившаяся на данный момент культура и методология решения задач биомеханики зубочелюстного аппарата с использованием численных методов позволяет говорить о применимости линейно-упругих моделей мягких (пульпа, периодонт) и твердых (эмаль, дентин, кортикальная, альвеолярная и губчатая кости) тканей. При численной оценке прочности протезов доминирует использование статического анализа, применяющегося в случаях не зависящей от времени, медленно меняющейся нагрузки, действующей на биомеханическую систему [30, 32, 35, 40]. Однако исследования показывают, что жевательные нагрузки, воспринимаемые конструкцией протеза и окружающими тканями, изменяются с течением времени [46], возбуждая в системе колебательные процессы, а отклик зависит не только от значения действующей силы, но и от демпфирующих характеристик [16]. Актуальной становится задача исследования напряженно-деформированного состояния мостовидного зубного протеза методом конечных элементов с учетом временной зависимости силы, сопоставление с результатами статического анализа. Расширение возможностей математического инструментария позволит оценить динамический отклик и
напряженно-деформированное состояние системы, выявить закономерности, присущие исследуемым процессам и в дальнейшем сформулировать требования, предъявляемые к проектируемому изделию.
Несмотря на большое число работ, посвященных биомеханике зубочелюстной системы, очень малая часть из них посвящена исследованию динамических процессов. Так, в работе [45] на примере титанового дентального имплантата был численно проведен модальный анализ и показано, что характеристики динамического отклика системы зависят от времени действия импульса. В исследовании [16] численно изучался отклик системы, состоящей из резца верхней челюсти и окружающих тканей, при действии сосредоточенной силы, изменяющейся по синусоидальному закону в течение 8 мс. Показано, что демпфирующие свойства системы влияют на напряжения, возникающие в зубе, а максимальные динамические напряжения меньше статических напряжений. В то же время в работе [47] на примере титанового дентального имплантата показано, что динамические напряжения, возникающие в имплантате, выше статических. Определение собственных частот и форм колебаний сегментов верхней челюсти, а также зависимость собственных частот от модуля упругости периодонта и введения демпфирования подробно изложено в [24]. Подход к идентификации динамических нагрузок, действующих на мостовидные протезы двух исполнений при жевании, представлен в [46]. Влияние толщины окружающей зуб кости на резонансную частоту описан в [15]. Анализ отечественных и зарубежных работ обнаруживает недостаток освещения методов динамического анализа применительно к расчету мостовидных зубных протезов и дентальных протезов в целом. В то же время выстроенный аппарат теории механических колебаний [1, 4] позволяет предположить, что в мостовидных зубных протезах при определенных условиях нагружения могут возникать колебательные, близкие к резонансным, режимы. Игнорирование динамической природы сил, действующих на биомеханическую систему, может повлечь разрушение протеза раньше прогнозируемого срока. Следствием этого являются не только повторная операция, сопряженная с финансовыми затратами, но и травмы мягких тканей зуба, стресс и дискомфорт пациента.
В ряде работ при оценке напряженно-деформированного состояния зубных протезов и окружающих тканей влияние пульпы на распределение деформаций, напряжений и энергии не учитывается [6, 10, 24, 46, 48]. Несмотря на известные в медицине и биологии функции пульпы, ее механические функции остаются малоизученными. Например, в работе [14] утверждается, что пульпа играет роль естественного амортизатора зуба без приведения количественных оценок.
Настоящая работа является первым этапом на пути создания методологии проектирования и расчета зубных протезов, в которой устанавливается связь между результатами статического и динамического анализа, формулируются условия применимости каждого из видов анализа. Количественно характеризуется вклад пульпы в демпфирование вынужденных колебаний биомеханической системы. Задача решается в плоско-деформированной постановке на примере мостовидного зубного протеза, сделанного из полиметилметакрилата, широко использующегося в качестве материала для дентального протезирования благодаря хорошей биосовместимости и адгезии, отсутствию вкусовых свойств и нерастворимости в биологических жидкостях (кровь, слюна) [41]. На основании результатов экспериментального исследования керамического мостовидного зубного протеза схожей конструкции проводится верификация численной модели. Сопоставляются результаты статического и динамического анализа. Количественно характеризуется влияние пульпы на напряженно-деформированное состояние протеза и демпфирующие свойства биомеханической системы в целом.
Методы и материалы
Теоретические сведения
Разрешающее уравнение метода конечных элементов при решении статических задач имеет вид [2]
Ku = f, (1)
где K - матрица жесткости, u - вектор узловых перемещений, f - вектор узловых сил. Это классическое уравнение применяется при решении задач, когда силы, действующие на систему, не зависят от времени, т.е. постоянны, или являются медленно меняющимися функциями времени. Из теории колебаний [2, 4] известно, что при стремлении частоты возбуждающей силы к собственной частоте колебаний системы амплитуда колебаний увеличивается. Если сила действует в течение короткого интервала времени, сопоставимого с периодом свободных колебаний системы, то в системе могут возникать колебательные процессы. Уравнение движения системы с конечным числом n степеней свободы в векторно-матричной форме записи имеет вид
[1,46]
MÜ (t ) + Bü (t ) + Ku (t ) = f (t ),
<u(0) = u°, (2)
«(0) = ù°,
где u(V), ù(i), ü(i) - векторы узловых перемещений, скоростей и ускорений, М, В и К - матрицы масс, диссипации и жесткостей соответственно, u(0) = u°, ü(0) = ü°-
векторы узловых перемещений и скоростей в начальный момент времени t = 0.
Установлено [16, 24], что демпфирование оказывает значительное влияние на динамическое поведение зубочелюстных систем. Если распределение сил трения в системе заведомо не известно, то общепринятым [1, 12, 24, 46] является использование модели Рэлея (пропорциональное демпфирование), компоненты в которой определяются следующим образом:
B = aM + ßK. (3)
Коэффициенты a и ß отвечают за демпфирование низкочастотных и высокочастотных форм соответственно. Для i-й формы колебаний параметры a и ß связаны через безразмерный коэффициент демпфирования ç; следующим соотношением:
çi=■ (4)
где шг - круговая собственная частота i -й формы колебаний в системе без учета демпфирования. Известно [16], что коэффициенты ç; в задачах биомеханики зубочелюстного аппарата могут принимать значения от 0,091 до 0,24. По итогам модального анализа определяются n собственных форм колебаний, суммарно задействующих более 90% общей массы системы. Для задания ç; используется методика, предложенная в работе [12]. Выражения для коэффициентов матрицы диссипации в модели Рэлея выражаются следующим образом:
а = ■
^ n n r n \
2 S ©2 S 3L - nlL ©г Ь
V i=1 i=1 ©i i=1
n n 1
S и' S 7? -
i=1 i=1 ©г
(5)
n
f
2
ß =
1
V i=1 i=1 ©i
i=1 ©i J
n n 1
T ©2 TЛ - n
(6)
i=1
i=1
Для энергетической оценки демпфирующих свойств системы введем параметр, равный отношению энергии, рассеиваемой за один период колебаний, к максимальной упругой энергии - коэффициент поглощения у [1]:
W U - U
у = — = ■ 1 2
U
U
(7)
В этом выражении и, и - значения потенциальной энергии упругой деформации для первого и второго периода колебаний соответственно. Коэффициент поглощения позволяет оценить, насколько существенно демпфирующие параметры системы влияют на процесс затухания колебаний.
Геометрия модели
Объект исследования представлен мезиодистальным сечением участка нижней челюсти, зубного ряда и мостовидного протеза, включающего опорные искусственные коронки на первый премоляр и первый моляр, промежуточную часть - фасетку (недостающий второй премоляр), кортикальную и губчатую кости, периодонт и пульпу (рис. 1). Геометрия зубов получена по результатам компьютерной томографии с последующим воссозданием трехмерных моделей.
^=50Н
0,5г* г*
Губчатая кость Дентин
Кортикальная кость Периодонт
Полиметилметакрилат Пульпа
t, c
Рис. 1. Геометрия модели, нагрузка и граничные условия
Высота в сечении премоляра составляет 21 мм (6,19 мм - коронка, 14,81 мм -корень), моляра - 20,38 мм (6,19 мм - коронка, 14,19 - корень). Длина мостовидного протеза - 31 мм, толщина шейки - 5,4 мм, радиусы скруглений шейки - 0,4 мм. Толщина периодонта равномерна по длине корня зуба и составляет 0,3 мм [28, 29]. Толщина кортикального слоя также постоянна по длине корня зуба и равна 1,5 мм для премоляра и 2 мм для моляра [11, 20, 37].
Нагрузки и граничные условия
Для решения задачи определения напряженно-деформированного состояния протеза в статической и динамической постановках к бугорку на окклюзионной поверхности фасетки приложена сосредоточенная сила [40, 46] Г = 50 Н (см. рис. 1). При статическом расчете сила линейно возрастает от 0 до 50 Н. При решении динамической задачи сила линейно возрастает от 0 до 50 Н в течение времени t* / 2, после чего линейно убывает со значения 50 Н до нуля за то же время {* /2. Так как длительность действия нагрузки t* на зубы и протез может измеряться сотыми [46], тысячными [16] и десятитысячными [45] долями секунды, исследовался отклик модели при { = 2 • 10"5 с, ^ = 2 • 10"4 с, { = 2 • 10"3 с, ^ = 2-10""2 с, ^ = 2-10"1 с. На границы расчетной области наложены кинематические граничные условия [22, 42].
Задача экспериментальной верификации решалась в статической постановке, исходные данные приведены [48].
Конечно-элементная модель
Задача решалась в плоско-деформированной постановке [8-10, 21, 22, 32, 35, 38, 40]. Разбиение модели на конечные элементы и численное решение проводилось в программном пакете ЛЫ8У8. Для определения числа конечных элементов, в полной мере отражающих процесс распределения нагрузок в системе и динамический отклик при действии возбуждающей силы, проведен тест на сходимость значений эквивалентных (по Мизесу) напряжений от числа элементов [23]. После проведения теста на сходимость окончательное число 6-узловых [32] элементов РЬЛЫЕ183 составило 63 289 (рис. 2).
Рис. 2. Конечно-элементная модель (а, б, в - области сгущения сетки)
Размер элементов на окклюзионной поверхности и в областях скруглений шейки - 0,005 мм. Качество сетки контролировалось средними значениями параметров
_о
Element quality = 0,94, Skewness = 9,58-10 и Orthogonal quality = 0,94. Толщина слоя цемента между мостовидным протезом и поверхностями препарированных зубов считалась незначительной и в расчетах не учитывалась [6].
Материалы протезов, твердые и мягкие ткани нижней челюсти в процессе моделирования предполагались однородными, изотропными и линейно-упругими [47, 48]. Механические свойства материалов взяты из [7, 26, 34, 48] и представлены в таблице.
Механические свойства материалов
Материал Модуль Юнга Е, МПа Коэффициент Пуассона ц Плотность р, кг/м3
Кортикальная кость 12 200 0,26 1100
Губчатая кость 1220 0,31 270
Дентин 18 600 0,31 203
Пульпа 2 0,45 95
Периодонт 70 0,45 135
Полиметилметакрилат 3000 0,38 1190
Керамика 15 600 0,28 -
Результаты
Тест на сходимость решения
Результаты теста на сходимость решения представлены на рис. 3. Из графика видно, что выбранная конфигурация сетки обеспечивает незначительное расхождение в результатах (менее 1%) по сравнению с результатами, полученными на сетке с большей густотой.
60
3 Я
50
<u d 50
§ С
§ ^ 40
«
s «
m S
£ 5 30
i * л
i 1 о и л
« & 20
л
10
50 000 100 000 150 000 200 000 Количество элементов Окклюзионная поверхность Скругление шейки
Рис. 3. Результаты теста на сходимость напряжений в локациях протеза от числа элементов сетки
0
Экспериментальная верификация
Результаты экспериментальной верификации представлены на рис. 4.
Рис. 4. Результаты экспериментальной верификации двухмерной модели: а - зависимость «перемещение - сила», воспроизведенная по результатам работы
[48]; б - расхождение в результатах
Видно, что максимальное расхождение перемещений составляет менее 8%. Меньшие значения перемещений объясняются тем, что рассматриваемый протез имеет увеличенную толщину шейки (5,4 мм против 4,08 мм) по сравнению с аналогом из работы [48]. Таким образом, можно говорить о правомерности применения двумерной постановки для рассматриваемого класса задач.
Статический анализ
При решении задачи как в статической, так и в динамической постановке максимальные эквивалентные напряжения в протезе возникают на окклюзионной поверхности и в области скругления шейки ближе к моляру (рис. 5).
58,524 52,022 а 45,519
! 39,016 32,514 26,011 19,509 13,006 6,5034 0,0008
Рис. 5. Области максимальных эквивалентных напряжений в протезе, МПа (а, б - области наибольших эквивалентных напряжений)
В модели с учетом пульпы максимальные эквивалентные напряжения на окклюзионной поверхности составили 58,5 МПа, в шейке - 44,1 МПа. В модели без учета пульпы - 58,6 и 44,2 МПа соответственно. Перемещения точки приложения силы в модели с учетом пульпы - 0,0943 мм, без учета пульпы - 0,0949 мм. Различие между величинами составляет менее 1%, что свидетельствует о применимости упрощенной модели (без учета пульпы) при выполнении статических расчетов в задачах проектирования зубных протезов.
Динамический анализ
По результатам модального анализа установлено, что влияние демпфирования и учета пульпы на первую собственную частоту несущественно. Значение первой собственной частоты равно 12 000 Гц. Таким образом, период свободных
колебаний, соответствующий ^, равен Т = 8,3 • 105 с. Поэтому случай { = 8,3 • 105 с также был проанализирован. Для модели с пульпой число первых собственных форм колебаний, суммарно задействующих более 90% общей массы системы, равно 28. Для модели без учета пульпы - 23. Для модели с пульпой коэффициенты а = 4,2440 4, Р = 1,84 -10"6, для модели без учета пульпы - а = 6,13 -10"2, Р = 9,45 -10"7. Графики перемещений точки приложения силы представлены на рис. 6 (здесь и далее: на всех графиках синий цвет соответствует модели с пульпой, оранжевый - без пульпы).
Из графиков видно, что ярко выраженный затухающий колебательный процесс имеет место при кратковременных воздействиях на систему. В случае {* > 2 -10"3 с нагружение можно считать статическим, без возникновения в системе колебательных процессов, поэтому соответствующие графики не приводятся. Графики максимальных эквивалентных напряжений, возникающих на окклюзионной поверхности и в шейке, представлены на рис. 7 и 8 соответственно.
В модели с учетом пульпы перемещения и максимальные напряжения ниже, чем в модели без учета пульпы, а скорость затухания колебаний в первом случае значительно выше. Анализ зависимостей перемещений и напряжений позволяет сделать вывод о том, что пульпа обладает энергопоглощающей функцией и выступает природным демпфером в зубочелюстном аппарате. Амплитудные значения перемещений и напряжений для первого колебания (являющиеся максимальными) в моделях с учетом и без учета пульпы различаются не более чем на 3%. В случае, когда {* = 8,3 • 105 с, максимальные напряжения на окклюзионной поверхности достигают предела текучести (для полиметилметакрилата предел текучести при растяжении составляет 65 МПа). Именно в этой области будет происходить стирание материала. Важно отметить, что при {* = 8,3 • 105 с в системе начинают накапливаться пластические деформации, а на окклюзионной поверхности возникает чередование областей наибольших напряжений в зависимости от амплитудных положений фасетки зубного протеза. Этим обусловлен вид зависимости на рис. 7, б.
Коэффициент поглощения энергии
На рис. 9 показаны значения коэффициента поглощения в зависимости от длительности действия силы F. Из графика следует, что пульпа может поглощать до 30% накопленной за первый цикл колебаний энергии. Для случаев {* > 2 40"3 с предполагалось, что происходит полное рассеивание накопленной энергии, поэтому у = 1.
0,05
S 0
я
<D
I -0,05 а
<D
с
-0,10
0,05
к я
<D
3
S а
с:
-0,05
-0,10
s я
<D
3
<D
'S
<D Л
<D
с
0,05
0,00
-0,05
-0,10
А IL л. _
0,0002 0,0004 0,0006 0,0008 0,001 Время, с а
0,0002
0,0004 0,0006 Время, с
б
0,0008
0,001
J
\ ryW
V
0,0002 0,0004 0,0006 Время, с
0,0008
0,001
Рис. 6. Перемещения точки приложения силы при: а - £ = 2 • 10 с;
* —5 * —4
б - £ = 8,3 • 10 с; в - £ = 2 • 10 с. Знак «минус» соответствует перемещению
точки вниз
0
0
0
0
в
а > s £ °
<D
« S
S и S и со о
й §3
Л ® Я л
Я Я
о
NH га
Я Е
2 S
«
<D %
« Л
S «
s
H о
о «
а «
о с
60
50
40
30
20
10
0,0002
0,0004
Время, с
ЛЛЛУУУЧ^.
0,0006 0,0008
0,001
а > s
л о
« я
е
п
а «
s
S и
S и
со о
й §3
л ® Я л
S
с
И й
S
« «
о
NH га
я с: s S
S н
с
о я
X
р
е
«
о с
я
е
*
я р
s «
60
50
40
30
20
10
60
50
40
30
s
H
с
о «
р
ё 20
о
с
10
0,0002 0,0004 0,0006 Время, с б
0,0008
0,001
/
0,0002 0,0004 0,0006 Время, с в
0,0008
0,001
Рис. 7. Максимальные эквивалентные напряжения на окклюзионной поверхности
протеза при: а - t = 2 • 10 с; б - t = !
,3 • 10"5 с; в - t* = 2• 10"4 с
0
а
0
0
о ® §
Ё и
о о
ч ч
^у cd
х g
о
S
60
50
40
30
20
10
0
0,0000 0,0002 0,0004 0,0006
Время, с
0,0008
0,0010
щ
е
4
^ iS
* i о
а К
и
5 ^
S cj
д В
е
&
cd Щ
60
50
40
30
20
10
0,0000 0,0002
0,0004
0,0006
0,0008
0,0010
Время, с
б
е
и
Ё и ее ч ч л 1-н — оэ ^ J2 s л С са и ^ ьй о ^
Щ К
ЬЙ
S
S cj
д В
е
60
50
40
30
20
10
0
0,0000 0,0002 0,0004 0,0006
Время, с
0,0008
0,0010
Рис. 8. Максимальные эквивалентные напряжения на скруглениях шейки протеза
при: а - t = 2-10 5 с; б - t =
,3 -10-5 с; в - t = 2-10-4 с
а
0
в
«
<D
3
о «
1-
о с
£
<D
S
я
S
со
о «
1,0 0,8 0,6 0,4 0,2 0,0
10-4 0,001 0,010 0,100 Время действия силы, с
Рис. 9. Значения коэффициента поглощения энергии при разных режимах нагружения
Сравнение результатов статического и динамического анализа
На рис. 10 представлены зависимости амплитудных значений перемещений и напряжений для исследуемых режимов нагружения в моделях с учетом и без учета пульпы.
В случаях, когда время действия силы { > 103 с, динамические перемещения и напряжения стремятся к статическим. При значении {* = 105 с динамические перемещения и напряжения значительно меньше статических. При стремлении длительности действия силы к периоду свободных колебаний, соответствующих первой собственной частоте, наблюдается возрастание амплитудных значений динамических перемещений до 0,12 мм, напряжений - до 65 и 55 МПа соответственно. Таким образом, динамические напряжения могут превышать статические на 11%, достигая предела текучести материала.
Выводы
При решении задач биомеханики зубочелюстного аппарата с применением численных методов очень важно грамотно интерпретировать условия, в которых функционирует протез in vivo. Это определяет методику, инструментарий решения сформулированной модельной задачи и, как следствие, результат. В работе показано, что использование динамического анализа при оценке напряженно-деформированного состояния мостовидного зубного протеза в случаях действия кратковременных сил помогает глубже понять природу и назначение отдельных звеньев биомеханической системы. Так, установлено, что пульпа играет роль естественного демпфера, гасителя вынужденных колебаний в системе, что позволяет избежать дислокации корней зубов, которые могут привести к травме периодонта. Кроме того, показано, что важным этапом при проектировании мостовидных зубных протезов и имплантатов в целом является проведение модального анализа с целью определения собственных частот и резонансных режимов. Изменяющиеся во времени жевательные нагрузки, действующие на зубной протез, могут привести биомеханическую систему в колебательное движение с амплитудными значениями напряжений, близкими и достигающими предела текучести. Таким образом, оценка напряженно-деформированного состояния протеза под действием статических сил может быть заниженной, так как такой подход не учитывает всего спектра физиологических нагрузок.
а
4
ет
s с 2 d
« S а я
к Й л ю
R <D
I §
5 « О
и л
S
0,12
0,11
0,10
0,09
0,08
св
й с
и 5
^'S -
ha g н® н s 8 о
S «
Я
я о к
ййД
ч s 2
Й И со а
10-4 0,001 0,010 0,100 Время действия силы, с а
а *
с л
М М
о Й И ый S i
65
60
55
50
45
40
35
10
0,001
0,010
0,100
Время действия силы, с б
55
50
45
40
35
30
//
; /
ff §
•
- 9
10
,-4
0,001 0,010 0,100
Время действия силы, с в
Рис. 10. Максимальные амплитудные значения при разных режимах нагружения: а - перемещения; б - напряжения на окклюзионной поверхности; в - напряжения
в скруглениях шейки
Список литературы
1. Бидерман В.Л. Теория механических колебаний. - М.: Высшая школа, 1980. - 480 с.
2. Зенкевич О. Метод конечных элементов в технике. - М.: Мир, 1975. - 543 с.
3. Лабунец В.А., Диева Т.В., Лабунец О.В. Повозрастной характер распространенности дефектов зубных рядов и дефектов коронковой части зубов, требующих ортопедического лечения у лиц молодого возраста // Cardiovascular risk factors. - 2010. - Vol. 51, № 2. - С. 373-375.
4. Пановко Я.Г. Основы прикладной теории колебаний и удара. - Л.: Машиностроение, 1976. - 320 с.
5. Фелькер Е.В., Ячменева Л.А., Евдокимова Е.И. Распространенность и локализация дефектов зубных рядов среди населения г. Курска // Международный журнал экспериментального образования. -2015. - № 5. - С. 1.
6. Чумаченко Е.Н. Математическое моделирование напряженно-деформированного состояния зубных протезов: учеб. пособие / Е.Н. Чумаченко, С.Д. Арутюнов, И.Ю. Лебеденко. - М.: МГМСУ, 2003. -271 с.
7. Abdel-Wahab A.A., Ataya S., Silberschmidt V.V. Temperature-dependent mechanical behaviour of PMMA: Experimental analysis and modelling // Polymer Testing. - 2017. - Vol. 58. - P. 86-95. DOI: 10.1016/j.polymertesting.2016.12.016
8. Aggarwal S., Garg V. Finite element analysis of stress concentration in three popular brands of fiber posts systems used for maxillary central incisor teeth // Journal of Conservative Dentistry: JCD. - 2011. - Vol. 14, № 3. - С. 293-296. DOI: 10.4103/0972-0707.85819
9. Al-Zahawi A.R., Saeed K.M.M., Saeed H.M.M. The influence of shifting the class I cavity position prepared in posterior teeth buccally and lingualy on stress distribution (finite element analysis study) // Journal of Baghdad College of Dentistry. - 2015. - Vol. 27, № 1. - P. 11-17.
10. Amarante M.V., Pereirab M.V.S., Darwishc F.A.I., Camaraod A.F. Virtual analysis of stresses in human teeth restored with esthetic posts // Materials Research. - 2008. - Vol. 11, № 4. - P. 459-463.
11. Chen Y.C. Nature-inspired shape optimization of dental restoration for replacing single missing tooth. Ph.D. thesis [Электронный ресурс], 2014, 246 p. - URL: conservancy.umn.edu/bitstream/handle/11299/ 167010/Chen_umn_013 0E_15347.pdf%3B sequence=1 (дата обращения: 03.09.2018).
12. Chowdhury I., Dasgupta S.P. Computation of Rayleigh damping coefficients for large systems // The Electronic Journal of Geotechnical Engineering. - 2003. - Vol. 8. - P. 1-11.
13. Faggion C.M., Giannakopoulos N.N., Listl S. How strong is the evidence for the need to restore posterior bounded edentulous spaces in adults? Grading the quality of evidence and the strength of recommendations // Journal of Dentistry. - 2011. - Vol. 39, № 2. - P. 108-116. DOI: 10.1016/j.jdent.2010.11.002
14. Huang G.T.J. Dental pulp and dentin tissue engineering and regeneration-advancement and challenge // Frontiers in Bioscience (Elite Edition). - 2011. - Vol. 3. - P. 788-800.
15. Huang H.M., Chiu C.L., Yeh C.Y., Lee S.Y. Factors influencing the resonance frequency of dental implants1 // Journal of Oral and Maxillofacial Surgery. - 2003. - Vol. 61, № 10. - P. 1184-1188.
16. Huang H.M., Tsai C.Y., Lee H.F., Lin C.T., Yao W.C., Chiu W.T., Lee S.Y. Damping effects on the response of maxillary incisor subjected to a traumatic impact force: a nonlinear finite element analysis // Journal of Dentistry. - 2006. - Vol. 34, № 4. - P. 261-268.
17. Hyde S., Dupuis V., Mariri B.P., Dartevelle S. Prevention of tooth loss and dental pain for reducing the global burden of oral diseases // International Dental Journal. - 2017. - Vol. 67. - P. 19-25. DOI: 10.1111/idj.12328
18. Jin L.J., Lamster I.B., Greenspan J.S., Pitts N.B., Scully C., Warnakulasuriya S. Global burden of oral diseases: emerging concepts, management and interplay with systemic health // Oral Diseases. - 2016. -Vol. 22, № 7. - P. 609-619. DOI: 10.1111/odi.12428
19. Kassebaum N.J., Smith A.G.C., Bernabé E., Fleming T.D., Reynolds A.E., Vos T., Murray C.J.L., Marcenes W., GBD 2015 Oral Health Collaborators. Global, regional, and national prevalence, incidence, and disability-adjusted life years for oral conditions for 195 countries, 1990-2015: a systematic analysis for the global burden of diseases, injuries, and risk factors // Journal of Dental Research. - 2017. - Vol. 96, № 4. - P. 380-387. DOI: 10.1177/0022034517693566
20. Kim J.H., Park Y.C. Evaluation of mandibular cortical bone thickness for placement of temporary anchorage devices (TADs) // The Korean Journal of Orthodontics. - 2012. - Vol. 42, № 3. - P. 110-117. DOI: 10.1186/s13005-015-0058-9
21. Lemos C.A.A., Almeida D.A.F., Batista Victor E.S., Mello C.C., Verri F.R., Pellizzer E.P., Mazaro J.V.Q. Influence of diameter and intraradicular post in the stress distribution. Finite element analysis // Revista de Odontologia da UNESP. - 2016. - Vol. 45, №. 3. - P. 171-176. DOI: 10.1590/1807-2577.21215
22. Li Q., Ichim I., Loughran J., Li W., Swain M., Kieser J. Numerical simulation of crack formation in all ceramic dental bridge // Key Engineering Materials. - 2006. - Vol. 312. - P. 293-298.
23. Li W., Swain M.V., Li Q., Steven G.P. Towards automated 3D finite element modeling of direct fiber reinforced composite dental bridge // Journal of Biomedical Materials Research Part B: Applied Biomaterials. - 2005. - Vol. 74, № 1. - С. 520-528. DOI: 10.1002/jbm.b.30233
24. Liao Z., Chen J., Zhang Z., Li W., Swain M., Li Q. Computational modeling of dynamic behaviors of human teeth // Journal of Biomechanics. - 2015. - Vol. 48, № 16. - P. 4214-4220. DOI: 10.1016/j.jbiomech.2015.10.019
25. Listl S., Galloway J., Mossey P.A., Marcenes W. Global economic impact of dental diseases // Journal of Dental Research. - 2015. - Vol. 94, № 10. - С. 1355-1361. DOI: 10.1177/0022034515602879
26. Merdji A., Della N., Benaissa A., Bachir Bouiadjra B.-A., Serier B., Mootanah R., Muslih I., Mukdadi O.M. Numerical analysis of dental caries effect on the biomechanical behavior of the periodontal system // Journal of Nanotechnology in Engineering and Medicine. - 2015. - Vol. 6, № 3. - P. 031004. DOI: 10.1115/1.4032689
27. Misch C.E. Dental implant prosthetics. - Amsterdam: Elsevier Health Sciences, 2014.
28. Nanci A. Ten cate's oral histology development, structure, and function, 8/e. - Elsevier India, 2012.
29. Nanci A., Bosshardt D.D. Structure of periodontal tissues in health and disease // Periodontology. - 2000. -2006. - Vol. 40, № 1. - P. 11-28. DOI: 10.1111/j.1600-0757.2005.00141.x
30. Piccioni M.A.R.V., Campos E.A., Saad J.R.C., Andrade M.F., Galvao M.R., Rached A.A. Application of the finite element method in dentistry // RSBO Revista Sul-Brasileira de Odontologia. - 2013. - Vol. 10, № 4. - P. 369-377.
31. Pilliar R.M., Deporter D.A., Ватсон P.A., Валикетт N. Dental implant design-effect on bone remodeling // Journal of Biomedical Materials Research. - 1991. - Vol. 25, № 4. - P. 467-483. DOI: 10.1002/jbm.820250405
32. Poiate I.A.V.P., Vasconcellos A.B., Mori M., Poiate E.Jr. 2D and 3D finite element analysis of central incisor generated by computerized tomography // Computer Methods and Programs in Biomedicine. -2011. - Vol. 104, № 2. - P. 292-299. DOI: 10.1016/j.cmpb.2011.03.017
33. Pol C.W.P., Raghoebar G.M., Kerdijk W., Boven G.C., Cune M.S., Meijer H.J.A. A systematic review and meta-analysis of 3-unit fixed dental prostheses: Are the results of 2 abutment implants comparable to the results of 2 abutment teeth? // Journal of Oral Rehabilitation. - 2018. - Vol. 45, № 2. - P. 147-160. DOI: 10.1111/joor. 12575
34. Ridwan-Pramana A., Marcian P., Borak L., Narra N., Forouzanfar T., Wolff J. Finite element analysis of 6 large PMMA skull reconstructions: A multi-criteria evaluation approach // PloS one. - 2017. - Т. 12, № 6. -P. e0179325. DOI: 10.1371/journal.pone.0179325
35. Romeed S.A., Fok S.L., Wilson N.H.F. A comparison of 2D and 3D finite element analysis of a restored tooth // Journal of Oral Rehabilitation. - 2006. - Vol. 33, № 3. - P. 209-215. DOI: 10.1111/j.1365-2842.2005.01552.x
36. Sailer I., Makarov N.A., Thoma D.S., Zwahlen M., Pjetursson B.E. All-ceramic or metal-ceramic tooth-supported fixed dental prostheses (FDPs). A systematic review of the survival and complication rates. Part I: Single crowns (SCs) // Dental Materials. - 2015. - Vol. 31, № 6. - С. 603-623. DOI: 10.1016/j.dental.2015.02.011
37. Sathapana S., Forrest A., Monsour P., Naser-ud-Din S. Age-related changes in maxillary and mandibular cortical bone thickness in relation to temporary anchorage device placement // Australian Dental Journal. -2013. - Vol. 58, № 1. - P. 67-74. DOI: 10.1111/adj.12018
38. Shinya A., Yokoyama D. Finite element analysis for dental prosthetic design [Электронный ресурс] // Finite element analysis. Ed. by David Moratal. - InTech, 2010. - URL: www.intechopen.com/books/finite-element-analysis/finite-element-analysis-for-dental-prosthetic-design (дата обращения: 03.09.2018). DOI: 10.5772/10006
39. Silva-Junior M.F., Batista M.J., de Sousa M.L.R. Incidence of tooth loss in adults: a 4-year population-based prospective cohort study // International Journal of Dentistry. - 2017. - Vol. 2017. - Article ID 6074703. -7 p. DOI: 10.1155/2017/6074703
40. Smielak B., Swiniarski J., Wolowiec-Korecka E., Klimek L. 2D-Finite element analysis of inlay-, onlay bridges with using various materials // Archives of Materials Science and Engineering. - 2016. - Vol. 79, Iss. 2. - P. 71-78.
41. Spasojevic P., Zrilic M., Panic V., Stamenkovic D., Seslija S., Velickovic S. The mechanical properties of a poly (methyl methacrylate) denture base material modified with dimethyl itaconate and di-n-butyl itaconate // International Journal of Polymer Science. - 2015. - Vol. 2015. - Article ID 561012. - 9 p. DOI: 10.1155/2015/561012
42. Teixeira E.R., Sato Y., Akagawa Y., Shindoi N. A comparative evaluation of mandibular finite element models with different lengths and elements for implant biomechanics // Journal of Oral Rehabilitation. -1998. - Vol. 25, № 4. - P. 299-303. DOI: 10.1111/j.1365-2842.1998.00244.x
43. Tyrovolas S., Koyanagi A., Panagiotakos D.B., Haro J.M., Kassebaum N.J., Chrepa V., Kotsakis G.A. Population prevalence of edentulism and its association with depression and self-rated health // Scientific Reports. - 2016. - Vol. 6. - P. 37083. D01:10.1038/srep37083
44. Vaillancourt H., Pilliar R.M., McCammond D. Factors affecting crestal bone loss with dental implants partially covered with a porous coating: a finite element analysis // International Journal of Oral & Maxillofacial Implants. - 1996. - Vol. 11, № 3. - P. 351-359.
45. Williams K.R., Williams A.D.C. Impulse response of a dental implant in bone by numerical analysis // Biomaterials. - 1997. - Vol. 18, № 10. - P. 715-719.
46. Zhang D., Zhang Z., Liu J., Jiang C., Yoda N., Meng X., Li Q. Identification of dynamic load for prosthetic structures // International Journal for Numerical Methods in Biomedical Engineering. - 2017. - Vol. 33, № 12. - P. e2889. DOI: 10.1002/cnm.2889
47. Zhang J.K., Chen Z.Q. The study of effects of changes of the elastic modulus of the materials substitute to human hard tissues on the mechanical state in the implant-bone interface by three-dimensional anisotropic finite element analysis // West China J. Stomatol. - 1998. - Vol. 16. - P. 274-278.
48. Zhang Z., Thompson M., Field C., Li W., Li Q., Swain M.V. Fracture behavior of inlay and onlay fixed partial dentures - an in-vitro experimental and XFEM modeling study // Journal of the Mechanical Behavior of Biomedical Materials. - 2016. - Vol. 59. - P. 279-290. DOI: 10.1016/j.jmbbm.2016.01.035
MODELLING OF DYNAMIC BEHAVIOUR OF DENTAL BRIDGE USING
FINITE ELEMENT METHOD
A.E. Krupnin, Y.N. Kharakh, L.G. Kirakosyan, S.D. Arutyunov (Moscow, Russia)
The goal of research is application of dynamic analysis methods in dentofacial biomechanics using polymethylmethacrylate dental bridge as a pilot example. Plane strain formulation is used. Experimental verification is carried out based on results of ceramic dental bridge with identical design. Natural frequencies and Rayleigh damping ratios are determined. Influence of force duration on stresses occurring in dental bridge is investigated. Local areas of maximal equivalent stresses are determined using von Mises criterion. Moreover, the study shows high increase in deformation and stresses (excessing the yield stress) in the case when force duration equal or close to natural period. The results of dynamic analysis tend to results of static analysis in the case of force duration increase. The influence of pulp can be neglected during solving static analysis problems. But, the research approved pulp's role as a natural shock absorber in the dentofacial system. As a result, consideration of pulp leads to decrease of deformation, stresses, intense energy absorption and oscillation damping.
Key words: dentofacial biomechanics, dental bridge, finite element method, dynamics, modal analysis, damping.
Получено 4 сентября 2018
