DEPENDENCE OF PRIMARY STABILITY OF DENTAL IMPLANT ON THE DIRECTION OF THE OCCLUSAL LOAD
- Authors: Dashevskiy I.N1, Shushpannikov P.S1
- Affiliations:
- Issue: Vol 23, No 3 (2019)
- Pages: 391-399
- Section: Articles
- URL: https://ered.pstu.ru/index.php/rjb/article/view/2351
- DOI: https://doi.org/10.15593/RZhBiomeh/2019.3.05
- Cite item
Abstract
Dental implantology practice shows that oblique load on the implant is more dangerous than the vertical one. Of particular importance is the study of the implant-jaw system under oblique load studying the primary stability of dental implants, when osseointegration has not yet occurred and there is no adhesion at the bone-implant interface. The damaging effect of excessive loads in this case is associated with the danger of superfluous micromotions on the bone-implant interface (mutual displacements of corresponding points) arising under their action, which leads to disruption of the osseointegration process. On the simplest model of a dental implant with a square thread profile, the influence of the load tilt angle on the primary stability of the implant in the jaw is examined. Calculations are carried out by the finite element method in the ANSYS package. Graphs of changes in micromotion along the bone-implant interface are presented, tables of magnitudes and localizations of micromotion maximums and graphs of their dependence on the loading angle are given. It is shown that for non-integrated implants, transition from vertical to horizontal loading leads to a sharp decrease in implant stability, which at sufficiently high occlusal loads impairs osseointegration.
Full Text
Введение В дентальной имплантологии хорошо известно, что косая нагрузка на имплантат опаснее вертикальной [4] (то же справедливо и для нативных зубов [12]). Имеется и ряд публикаций, касающихся расчета напряженно-деформированного состояния челюсти либо ее фрагментов с имплантатами под действием наклонной нагрузки [6-9, 11, 13, 15, 17, 22]. Большинство из них посвящено остеоинтегрированным (сросшимся с костью) имплантам. Однако особое значение исследование системы «имплантат - челюсть» под действием косой нагрузки имеет при изучении первичной стабильности дентальных имплантатов, когда остеоинтеграция еще не произошла и на интерфейсе «кость -имплант» нет сцепления [21, 14, 16-20, 23]. Повреждающее воздействие чрезмерных нагрузок при этом связано с опасностью возникновения под их влиянием избыточных микродвижений (взаимных смещений ответных точек) на интерфейсе «кость - имплант», что ведет к нарушению процесса остеоинтеграции. В работах [1, 2] на минимальной модели имплантатов в челюсти было рассмотрено влияние на первичную стабильность имплантата и микродвижения на интерфейсе характеристик резьбы при вертикальной нагрузке. В настоящей работе на той же модели изучается влияние на микроподвижность угла наклона нагрузки. Постановка задачи В качестве базового расчетного варианта был принят тот же, что и в [1, 2], но здесь сила P той же абсолютной величины образует с вертикалью угол α (α = 30°, 45° и 90°). Поэтому, в отличие от [1, 2], рассматриваемая задача не является осесимметричной (но все еще обладает одной плоскостью симметрии). Имея в виду возможную экспериментальную проверку, краевые условия поставили аналогично работе [3]: образец (цилиндр из кости с ввинченным имплантатом) как бы вставлен в жёсткий гладкий стакан (обойму), дно образца приклеено к стакану. Длина и диаметр обоймы приняты соответственно L = 30 мм, D = 20 мм [3]. Сводка всех использованных значений параметров базового расчетного варианта приведена ниже: E, ГПа ν l, мм d, мм p, мм h, мм w, мм L, мм D, мм P, Н α, град титан кость титан кость 110 1 0,3 0,3 8 4 1 0,2 0,2 30 20 700 0, 90, 45, 30 Здесь E, ν - модуль Юнга и коэффициент Пуассона, l, d - габариты импланта (соответственно длина и диаметр); p, h, w - характеристики резьбы: соответственно шаг (период) p (pitch), глубина (высота) h и ширина основания зубца (нитки резьбы) w. Профиль резьбы был принят квадратным, поскольку в [1, 2] было показано, что в случае вертикальной нагрузки именно он давал минимальные микродвижения на интерфейсе. На интерфейсе «имплантат - кость» задавалось условие скольжения. а б Рис. 1. Нагрузка на имплантат и точка максимальной микроподвижности A (а); общий вид импланта с четвертным вырезом (б) Расчёты выполнялись с использованием метода конечных элементов в программном комплексе ANSYS (версия 15.0). При построении конечно-элементной сетки использовались билинейные восьмиузловые конечные элементы в форме прямоугольного параллелепипеда. Вследствие наличия в задаче плоскости симметрии разбиению подлежала только часть (половина) модели, ограниченная этой плоскостью. Общее число элементов, использованных при разбиении, приблизительно равнялось 200 000. Для моделирования контакта между костью и имплантатом использовался метод штрафа. Жесткость (модуль Юнга) имплантата на два порядка больше жесткости кости. Поэтому при действии на верхний торец имплантата горизонтальной силы имплантат практически как твердое тело наклоняется и поворачивается вокруг некоторой точки вблизи своего нижнего торца, называемой центром сопротивления [5]. При этом с одной (лицевой) стороны имплантат будет вдавливаться в кость (и здесь величины микродвижений будут минимальны - точка B), а с другой (тыльной) - отстанет от кости - и здесь микродвижения будут максимальны (точка A). Наибольший интерес представляют места и величины максимальных микродвижений, поэтому в качестве фокусной точки были выбраны точка A (рис. 1), отвечающая максимуму микродвижений в определяемом ею поперечном сечении имплантата, и соответствующие этой точке вертикальное сечение и путь AD (рис. 1). Были проведены расчеты нормальной и тангенциальной компонент, а также модуля полного вектора микродвижений |u*| на интерфейсе вдоль пути AD и аналогичного пути BD. Результаты расчетов и их обсуждение На рис. 2, а-г представлены графики |u*| вдоль AD соответственно для α = 0°, 90°, 45°, 30°. На рис. 2, д для сравнения приведен график |u*| вдоль BD для α = 90°. Здесь первый штриховой столбик - верхняя грань зубца (нитки резьбы), второй столбик - боковая грань, третий - нижняя грань. Прежде всего из сравнения рис. 2, б и 2, д видно, что микродвижения в точке B на порядок меньше, чем в точке A. Далее, из практики хорошо известно, что главные проблемы (вплоть до резорбции костной ткани) чаще всего возникают наверху имплантата, в «пришеечной» области. Полученные расчетные результаты могут служить возможным объяснением такого явления: из графиков 3, а-г видно, что при вертикальной нагрузке максимальные перемещения возникают в апексной зоне, а при горизонтальной и косой нагрузке - наверху имплантата, при этом микродвижения при горизонтальной нагрузке на порядок больше, чем при вертикальном нагружении. Тот факт, что горизонтальные микродвижения больше вертикальных, отражает меньшую жесткость конструкции при горизонтальном нагружении в сравнении с вертикальным, что, в свою очередь, обусловлено геометрически. Далее известно [14, 21], что лишь чрезмерная микроподвижность прямо связана с нарушением остеоинтеграции и формированием волокнистой инкапсуляции. Согласно [21] допустимый порог микродвижений лежит между 50 и 150 мкм: при |u*| < 50 мкм обеспечена остеоинтеграция, при |u*| > 150 мкм всегда происходит фиброинтеграция, при 50 мкм < |u*| < 150 мкм результат зависит от других (дополнительных) факторов. В наших расчетах значения микродвижений при изменении угла нагрузки от 0° (вертикальная) до 90° (горизонтальная) менялись больше чем на порядок (от 10-12 до 180 мкм) и таким образом переходили из диапазона остеоинтеграции через промежуточную зону в область фиброинтеграции. Впрочем, надо оговориться, что принятые значения нагрузки близки к экстремальным [4] и, следовательно, в норме возникают редко. Считается, что типичные значения нагрузки P ~ 200 Н [4, 10], а при таких усилиях микродвижения будут примерно в 3,5 раза меньше, т.е. около 50 мкм, что, тем не менее, лежит на границе гарантированной остеоинтеграции. Косая нагрузка представляет собой комбинацию вертикальной и горизонтальной нагрузок. Поскольку перемещения от вертикальных нагрузок на порядок меньше, чем от горизонтальных, при не слишком малых углах α ими можно пренебречь и учитывать лишь перемещения от горизонтальных усилий, которые равны P→ = Psin α. Дистанция вдоль пути A, мм Дистанция вдоль пути B, мм Дистанция вдоль пути A, мм а б Дистанция вдоль пути B, мм в г Дистанция вдоль пути A, мм д Рис. 2. Изменение абсолютной величины микродвижений на интерфейсе «кость - имплантат» вдоль пути AD при α = 0° (а), 90° (б), 45° (в), 30° (г), то же для BD α = 90° (д), w = 0,2 мм; h = 0,2 мм; p = 1,0 мм При изменении силы P площадки контакта будут меняться, и таким образом задача нелинейна по P. Если эта нелинейность невелика, то график зависимости u (sin α) должен быть близок к прямой. В таблице приведены значения максимальных относительных перемещений (микродвижений) на интерфейсе в зависимости от угла приложения нагрузки, а на рис. 3 - те же данные показаны графически. Видно, что обе зависимости (u (α), u (sin α)) существенно нелинейны, однако, в отличие от первой, график второй не имеет точек перегиба (сохраняет знак кривизны). Максимальные микродвижения на интерфейсе «кость - имплант» при разных углах нагрузки Направление нагрузки, град (угол α к вертикали) sin α Максимум перемещений, мкм Положение максимума, (виток резьбы) 00 0 10,71 На последнем 30 0,5 67,36 На первом (шейка) 45 108,65 На первом (шейка) 90 1 184,34 На первом (шейка) α, град sin α Рис. 3. Зависимость максимальных значений микродвижений на интерфейсе «кость - имплантат» от угла нагружения Для того чтобы лучше понять полученный результат, дадим для смещений имплантата оценку на основе простейшей качественной плоской модели (рис. 4), в которой будем считать, что имплантат - жесткий гладкий стержень, а кость будем описывать как винклеровский слой. Как уже было отмечено выше, при действии на верхний торец имплантата горизонтальной силы имплантат практически как твердое тело наклоняется и поворачивается вокруг некоторой точки вблизи своего нижнего торца, называемой центром сопротивления. Такой поворот описывается соотношением u(y) = k(y - y0), где y - ордината, отсчитываемая от нижнего торца имплантата; u - смещение точек стержня по горизонтали (оси x); y0 - центр сопротивления; k - тангенс угла поворота. Две неизвестные y0 и k определяются из двух уравнений баланса: сил и моментов. Для расстояния между имплантатом и жесткой стенкой введем обозначение Рис. 4. Поворот имплантата под действием на его верхний торец горизонтальной силы P Тогда уравнение баланса сил даст Из уравнения баланса моментов (1) В рассмотренной элементарной модели максимальные микродвижения (равные зазору-отставанию имплантата от кости в точке A) как раз равны абсолютным перемещениям имплантата, и таким образом разница между точным и оценочным расчетом - в 4 раза. Эта разница связана с неучетом в элементарной плоской модели сопротивления областей кости, расположенных снизу и по бокам от имплантата, тем, что в точной модели зазор будет уменьшаться за счет выдавливания в него кости и др. Выведенные формулы с точностью до числового множителя могут быть получены безо всяких вычислений из соображений размерности и линейности. Параметры элементарной модели u || P, E, l, d, δ (в конечно-элементной модели есть еще ν, габариты обоймы, параметры резьбы). В силу линейности и из качественных соображений Далее по размерности d может войти единственным образом - в знаменатель, и с точностью до числового множителя получим, как и в (1), Конечно, все эти качественные формулы малопригодны для оценки реальных величин интересующих параметров. Их ценность в том, что они приближенно (в тенденции), зато в явном виде показывают характер влияния некоторых основных параметров на нужные величины. Здесь оценки были проведены для перемещений и микродвижений, но то же можно сделать и для напряжений. Отметим, что принятое для δ значение (δ = 8 мм) соответствует стабильности (перемещениям) в щечно-язычном направлении либо перемещениям одиночного импланта в медиально-латеральном направлении в отсутствие смежных зубов; при проведении подобных расчетов и оценок на практике принимаемые значения параметров должны отвечать конкретной ситуации. Выводы Изменение направления действия силы на имплантат от вертикального к горизонтальному в отсутствие сцепления на интерфейсе «кость - имплантат» (при немедленном нагружении, при неполной остеоинтеграции) приводит к драматическому - больше чем на порядок - росту микродвижений (взаимных смещений ответных точек) на интерфейсе и перемещению их максимума с апекса на шейку имплантата. Если действующие при этом окклюзионные нагрузки достаточно велики, существует опасность возникновения на интерфейсе чрезмерных микродвижений, что может привести к нарушению остеоинтеграции, в особенности в пришеечной зоне импланта.About the authors
I. N Dashevskiy
P. S Shushpannikov
References
- Дашевский И.Н., Шушпанников П.С. Влияние профиля резьбы на первичную стабильность дентальных имплантатов // Методы компьютерной диагностики в биологии и медицине - 2017. Материалы Всерос. школы-семинара. - Саратов: Саратовский источник, 2017. - С. 133-136.
- Дашевский И.Н., Шушпанников П.С. Влияние характеристик резьбы на первичную стабильность дентальных имплантатов // Российский журнал биомеханики. - 2018. - Т. 22, № 3. - С. 361-377.
- Ерошин В.A., Джалалова М.В. Напряженно-деформированное состояние биомеханической системы имплантат - упругое основание // Российский журнал биомеханики. - 2012. - Т. 16, № 3. - С. 82-95.
- Миш К.Е. Ортопедическое лечение с опорой на дентальные имплантаты. - М.: Рид Элсивер, 2010.
- Наумович С.А., Ивашенко С.В., Босяков С.М., Крушевский А.Е. Биомеханика системы зуб -периодонт - костная ткань. - Минск, 2009.
- Олесова В.Н., Бронштейн Д.А., Берсанов Р.У., Лернер А.Я. Результаты трехмерного математического моделирования функциональных напряжений в имплантате и несъемной протезной конструкции // Международный журнал прикладных и фундаментальных исследований. - 2014. - № 2. - С. 135-139.
- Олесова В.Н., Бронштейн Д.А., Лернер А.Я., Олесов Е.Е., Бобер С.А., Узунян Н.А. Напряженно-деформированное состояние в протезной конструкции на дентальном имплантате при цементной фиксации искусственной коронки // Российский журнал биомеханики. - 2016. - Т. 20, № 4. - С. 311-315.
- Олесова В.Н., Бронштейн Д.А., Магамедханов Ю.М., Кузнецов А.В., Каирбеков Р.Д., Заславский С.А. Экспериментально-математическое изучение функциональных параметров нижнего зубного ряда // Стоматология для всех. - 2011. - № 2. - С. 45-47.
- Олесова В.Н., Дубинский С.И., Бронштейн Д.А., Магамедханов Ю.М., Кащенко П. В., Юффа Е.П. Сравнительное математическое моделирование прочностных и деформационных параметров металлокерамических коронок с винтовой и цементной фиксацией к имплантатам // Кубанский научный медицинский вестник. - 2013. - № 6 (141). - С. 140-142.
- Параскевич В.Л. Дентальная имплантология. - М., 2006.
- Перельмутер М.Н. Исследование напряженно-деформированного состояния стоматологических имплантатов методом граничных интегральных уравнений // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2018. - № 2. - С. 83-95. doi: 10.15593/perm.mech/2018.2.08
- Шварц А.Д. Биомеханика и окклюзия зубов. - М.: Медицина, 1994.
- Шумаков Ф.Г., Олесова В.Н., Захаров П.А., Печенихина В.С., Гришков М.С. Сравнительная биомеханика керамического и титанового внутрикостных дентальных имплантатов // Российский вестник дентальной имплантологии. - 2017. - № 3-4 (37-38). - С. 4-7.
- Davarpanah M., Szmukler-Moncler S. Immediate loading of dental implants: theory and clinical practice. - Quintessence Publishing, 2009. - 356 p.
- Desai S.R., Desai M.S., Katti G., Karthikeyan I. Evaluation of design parameters of eight dental implant designs: A two-dimensional fnite element analysis // Niger J. Clin. Pract. - 2012. - Vol. 15. - P. 176-181.
- Hong D.G.K., Oh J.H. Recent advances in dental implants // Maxillofac Plast Reconstr Surg. - 2017. - Vol. 5, № 39 (1). - P. 33. doi: 10.1186/s40902-017-0132-2
- Horita S., Sugiura T., Yamamoto K., Murakami K., Imai Y., Kirita T. Biomechanical analysis of immediately loaded implants according to the “All-on-Four” concept // Journal of Prosthodontic Research. - 2017. - Vol. 61. - P. 123-132.
- Huang H.L., Hsu J.T., Fuh L.J., Tu M.G., Ko C.C., Shen Y.W. Bone stress and interfacial sliding analysis of implant designs on an immediately loaded maxillary implant: a non-linear finite element study // J. Dent. - 2008. - Vol. 36. - P. 409-417.
- Ma P., Xiong W., Tan B. Influence of thread pitch, helix angle, and compactness on micromotion of immediately loaded implants in three types of bone quality: a three-dimensional finite element analysis // BioMed Research International. 2014. - Vol. 2014. - Article ID 983103. - P. 13. doi: 10.1155/2014/983103
- Ryu H.S., Namgung C., Lee J.H., Lim Y.J. The influence of thread geometry on implant osseointegration under immediate loading: a literature review // J. Adv. Prosthodont. - 2014. - Vol. 6. - P. 547-554.
- Szmukler-Moncler S., Salama H., Reingewirtz Y., Dubruille J.H. Timing of loading and effect of micromotion on bone-dental implant interface: review of experimental literature // Journal of Biomedical Materials Research. Part A. - 1998. - Vol. 43, № 2. - P. 192-203.
- Winter W., Klein D., Karl M. Effect of model parameters on finite element analysis of micromotions in implant dentistry // Journal of Oral Implantology. - 2013. - Vol. XXXIX, № 1. - P. 23-29.
- Wu J.C., Chen C.S., Yip S.W., Hsu M.L. Stress distribution and micromotion analyses of immediately loaded implants of varying lengths in the mandible and fibular bone grafts: a three-dimensional finite element analysis // Int. J. Oral Maxillofac. Implants. - 2012. - Vol. 27. - P. 77-84.
Statistics
Views
Abstract - 61
PDF (Russian) - 48
Refbacks
- There are currently no refbacks.