ВЛИЯНИЕ ПЛОТНОСТИ КОСТНЫХ ТКАНЕЙ НА НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ ВБЛИЗИ ДЕНТАЛЬНЫХ ИМПЛАНТАТОВ

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2023. Т. 27, № 1. С. 189-201 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1976

EDN: VDVFUD

УДК 539.32

Влияние плотности костных тканей на напряженно-деформированное состояние вблизи дентальных имплантатов

М. Н. Перельмутер

Институт проблем механики им. А. Ю. Ишлинского РАН,

Россия, 119526, Москва, просп. Вернадского, 101, корп. 1.

Аннотация

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

Ключевые слова: имплантат, костные ткани, напряженно-деформированное состояние, метод граничных интегральных уравнений.

Получение: 26 декабря 2022 г. / Исправление: 19 февраля 2023 г. / Принятие: 17 марта 2023 г. / Публикация онлайн: 24 марта 2023 г.

Математическое моделирование, численные методы и комплексы программ Краткое сообщение

© Коллектив авторов, 2023 © СамГТУ, 2023

3 © Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования

Перельмутер М. Н. Влияние плотности костных тканей на напряженно-деформированное состояние вблизи дентальных имплантатов // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2023. Т. 27, № 1. С. 189-201. EDN: VDVFUD. DOI: 10.14498/vsgtu1976. Сведения об авторе

Михаил Натанович Перельмутер А https://orcid.org/0000-0002-8430-5412 доктор физико-математических наук; ведущий научный сотрудник; лаб. механики прочности и разрушения материалов и конструкций; e-mail: perelm@ipmnet.ru

Введение. Эффективность дентальной имплантации помимо прочих факторов [1] зависит от плотности костной ткани, которая может изменяться в широком диапазоне [2,3]. Ввиду этого изучение влияния плотности кортикальной и губчатой костных тканей челюсти на напряженно-деформированное состояние в зоне дентальной имплантации представляет большой научный и практический интерес [4, 5]. Изменение плотности костной ткани проводит к изменению ее модуля упругости [6]. Этот факт позволяет выполнять математическое моделирование влияния плотности костной ткани на напряженно-деформированное состояние при квазистатическом приложении нагрузки посредством соответствующего изменения модуля упругости кости [3,7-9,11].

Для моделирования напряженного состояния стоматологических имплан-татов и костных тканей в последние десятилетия используются численные методы. Наибольшее распространение получил метод конечных элементов (МКЭ) [12,13]. Метод граничных элементов (МГЭ) в задачах стоматологической биомеханики используется не так широко. Одно из первых применений этого метода было выполнено в 1993 году для анализа напряженного состояния соединения цилиндрического имплантата и костной ткани [14]. МГЭ также применялся для исследования напряженного состояния костных тканей и винтовых стоматологических имплантатов, изготовленных из различных материалов (титана или керамики) [15]. Сравнительный анализ применения МКЭ и МГЭ для анализа напряженного состояния имплантатов и костной ткани был выполнен в [16]. При использовании МГЭ ввиду независимого моделирования перемещений и напряжений на поверхности тела, в зонах с высокими градиентами напряжений можно получить результаты, сопоставимые с результатами, полученными с помощью МКЭ, но без существенного измельчения сетки элементов. При решении задачи МГЭ выполняется дискретизация только поверхности (границы — в двумерном случае) расчетной модели, что является одним из основных преимуществ данного метода перед МКЭ, при использовании которого необходимо разбиение на элементы всего объема тела. В данной работе МГЭ используется для анализа влияния изменения плотности костных тканей на напряженное состояние вблизи дентальных им-плантатов.

1. Постановка задачи и методика численного решения. Расчет напряженно-деформированного состояния выполнен в постановке плоской деформации для модели, состоящей из цилиндрического имплантата с коронкой и окружающих костных тканей при допущении, что костная ткань является изотропным и однородным упругим материалом. Значения модуля упругости кортикальной и губчатой костных тканей изменяются в достаточно широком диапазоне [17]. При выполнении данного исследования полагалось, что модуль упругости губчатой костной ткани изменяется в интервале Е3 = 0.5 ^ 8 ГПа, модуль упругости кортикальной костной ткани — Ес = 6 ^ 20 ГПа.

Расчеты выполнялись при фиксированном модуле упругости костной ткани одного типа и вариации модуля упругости костной ткани другого типа: 1) для каждого из трех значений модуля упругости кортикальной кости Ес = {6; 10; 18} ГПа выполнялся расчет при различных значениях модуля упругости губчатой кости Е3 = {0.5; 1; 3; 5; 8} ГПа;

2) для каждого из трех значений модуля упругости губчатой кости Е8 = = {0.5; 1; 3} ГПа выполнялся расчет при различных значениях модуля упругости кортикальной кости Ес = {6; 8; 10; 14; 18; 20} ГПа. Полагалось также, что коронка изготовлена из керамики, а остальные детали имплантата изготовлены из титана. Базовые механические свойства материалов имплантата и костных тканей приведены в табл. 1 по данным работ [10,11,17].

Таблица 1

Механические свойства имплантата и костных тканей [Mechanical properties of the implant and bone tissues]

Material Young's modulus E, GPa Poisson's ratio Yield strength (a), Tensile strength (b), MPa

Titanium 110 0.3 880

Ceramics 70 0.19 320

Cortical bone 18 0.3 130

Spongy bone 0.5 0.35 10 ^ 20

Расчетная модель имплантата и костных тканей содержит семь подобластей (см. рис. 1). Первые три подобласти моделируют костные ткани: 1 — внешний слой кортикальной кости (средняя толщина слоя 1.7 мм); 2 — губчатая кость — основная часть костной ткани челюсти; 3 — слой кортикальной кости толщиной 1 мм, формирующий стенки нижнечелюстного канала, имеющего диаметр 3 мм. Остальные подобласти соответствуют элементам конструкции имплантата (длина цилиндрической части имплантата— 12 мм, диаметр — 4.2 мм, см. подобласть 4). Вся модель полностью с указанием вида нагрузки и граничных условий, а также увеличенное изображение верхней части модели представлены на рис. 2 (для удобства графического представления внешний слой кортикальной кости, цилиндрическая часть имплантата и коронка заштрихованы). Между всеми подобластями модели полагается выполнение условий идеального сцепления, включая зону соединения цилиндрической части имплантата и костных тканей (полная остеоинтеграция). Нагрузка величиной р = 6 МПа, соответствующая усилию жевания [18], прикладывается к окклюзионной поверхности коронки в двух вариантах (в вертикальном направлении и под углом 45° к вертикали). Полагая, что площадь контактной поверхности коронки 5 ~ 0.25 |см| , можно считать, что заданное усилие р соответствует приложенной силе Е ~ 150 |Н|. По внешней границе кортикальной кости полагаем нулевые перемещения по направлениям осей координат, тонкая заштрихованная зона на рис. 2 показывает область заданных перемещений по внешней границе модели (это обозначение только на рисунке, не является частью модели).

Интенсивность напряжений (эквивалентные напряжения, напряжения по Мизесу) для состояния плоской деформации определяется следующим образом: _

/(^1 - )2 + (^2 - ^э)2 + (оз - ^)2

= V-2-'

где и 1, 02, 0з — главные напряжения.

Рис. 1. Подобласти модели: 1, 3 — кортикальная кость, 2 — губчатая кость, 4 — имплантат, 5 — винт, 6 —абатмент, 7— керамическая коронка под действием наклонной нагрузки

[Figure 1. The model subregions: 1, 3 — cortical bone, 2 — spongy bone, 4 — implant, 5 — screw, 6 — abatment, 7 — ceramics crown under inclined load]

Рис. 2. Дискретизация границ подобластей имплантата и костных тканей; действие наклонной нагрузки; 7 подобластей; 1106 узлов во всей модели

[Figure 2. Boundary elements discretization of the implant and surrounded bone tissues under inclined load application; 7 subregions; total 1106 nodes]

Для численного анализа напряженно-деформированного состояния в работе используется прямой вариант МГЭ [19]. Двумерная составная конструкция моделируется набором однородных подобластей. Такой подход позволяет рассматривать соединения различных материалов в конструкции, а также условия неидеального соединения между подобластями. Граничное интегральное уравнение (ГИУ) теории упругости для каждой однородной подобласти конструкции без учета массовых сил и температурных воздействий имеет вид [19]

Сц

где точки р и д принадлежат границе подобласти Г; ь,г(д), ¿¿(д) — перемещения и усилия на границе подобласти; с^ (р) —функции, зависящие от локальной геометрии в точке р; для гладкой границы сц(р) = /2. Функции (р,д) и ^ (р,д) являются фундаментальными решениями Кельвина для задачи теории упругости (плоская деформация) [19].

На участках идеального контакта подобластей полагаем, что выполняются условия непрерывности для перемещений и равновесия для усилий:

иЧ

где к и п — номера соседних подобластей.

Для численного решения ГИУ (1) используются изопараметрические квадратичные граничные элементы. Дискретизация границ подобластей и аппроксимация переменных позволяет заменить вычисление интегралов по всей поверхности подобласти в ГИУ (1) суммой интегралов по совокупности граничных элементов, представляющих поверхность подобласти, и получить для

каждой точки коллокации р дискретное представление ГИУ. Сингулярное поведение фундаментальных решений в ГИУ (1) проявляется в основном при вычислении интегралов по элементам, содержащим точку коллокации р. Поэтому рассматривается два случая расположения точки р: р / е и р G е, где е — граничный элемент, по которому выполняется интегрирование. В первом случае интегрирование выполняется по квадратурным формулам Гаус-са—Лежандра с учетом сгущения точек интегрирования в направлении точки коллокации. При р G е фундаментальное решение плоской задачи теории упругости для перемещений Gij (р, q) содержит слабую логарифмическую особенность и интегралы в дискретном представлении ГИУ, содержащие эту функцию, вычисляются с использованием квадратурной формулы с логарифмической весовой функцией. Фундаментальное решение для поверхностных усилий Fij (p,q) содержит сильную особенность и вычисление суммы соответствующих интегралов и параметра Cij (р) выполняется косвенным путем исходя из смещения тела как жесткого целого (см. описание алгоритма численного интегрирования для двумерных и пространственных задач в [20]). Обходя последовательно узловые точки границ всех подобластей, получаем дискретное представление ГИУ. Учет граничных условий на внешних границах подобластей и дополнительных условий на участках соединения подобластей (2) позволяет получить из системы дискретных ГИУ для всех подобластей конструкции систему линейных алгебраических уравнений для определения узловых неизвестных задачи. Методика численного решения ГИУ реализована в комплексе программ «МЕГРЭ» (Метод ГРаничных Элементов) [20-22], который использован для получения результатов, приведенных в данной статье.

2. Результаты расчетов при нормальной нагрузке на коронку.

Первый этап расчета выполнен при сжимающей нагрузке, приложенной к ок-клюзионной поверхности коронки. Зависимости максимальных эквивалентных напряжений в костных тканях от модулей упругости кортикальной и губчатой костной ткани представлены на рис. 3-6.

При увеличении модуля упругости губчатой костной ткани Es (увеличении жесткости кости) возрастает часть нагрузки, передаваемая на эту костную ткань. Поэтому максимальные напряжения в губчатой костной ткани возрастают, см. рис. 3. Напряжения в кортикальной костной ткани при увеличении модуля упругости губчатой кости Es снижаются за счет уменьшения нагрузки, передаваемой на эту часть кости, см. рис. 4. Отметим, что напряжения в губчатой кости снижаются при увеличении модуля упругости кортикальной кости Ес, см. рис. 3 и 5. С другой стороны, уровень напряжений в кортикальной кости возрастает при увеличении модуля упругости этой кости, см. рис. 4 и 6.

На рис. 7 представлены распределения интенсивности напряжений по границе губчатой кости для двух значений модуля упругости кортикальной кости. При малом значении модуля упругости губчатой кости максимум напряжений наблюдается вблизи зоны контакта губчатой и кортикальной костных тканей с имплантатом (шейка имплантата, см. рис. 1, подобласти 1 и 2, верхняя часть), а при увеличении модуля упругости губчатой кости положение максимума смещается в нижнюю зону контакта между имплантатом и губчатой костью. Отметим, что ввиду несимметрии коронки и возникновения

4.5

4.0

св Рч

§ 3.5

3.0

2.5

0

у f V

у 1-EC = 6 GPa 2 - Ec = 10 GPa 3 - Ec = 18 GPa

6

£ 45

ь 30 15

2 4

Е„ GPa

Рис. 3. Зависимость максимальных напряжений в губчатой кости от модуля упругости губчатой кости; сжатие

[Figure 3. The maximum stresses in the spongy bone vs the modulus of elasticity of the spongy bone; compression]

0

3 1 - Ec = 6 GPa 2 — Ec = 10 GPa 3 - Ec = 18 GPa

2\

A\\

0

4

, GPa

6

Рис. 4. Зависимость максимальных напряжений в кортикальной кости от модуля упругости губчатой кости; сжатие

[Figure 4. The maximum stresses in the cortical bone vs the modulus of elasticity of the spongy bone; compression]

3

1 — E, = 0.5 GPa 2 — Es = 1.0 GPa 3 — Es = 3.0 GPa

16

80

60

40

20

18

4 8 12

Ес, GPa

Рис. 5. Зависимость максимальных напряжений в губчатой кости от модуля упругости кортикальной кости; сжатие

[Figure 5. The maximum stresses in the spongy bone vs the modulus of elasticity of the cortical bone, compression]

1 — Es = 0.5 GPa 2 — E, = 1.0 GPa 3 — Es = 3.0 GPa

3

12

Ec, GPa

16

18

Рис. 6. Зависимость максимальных напряжений в кортикальной кости от модуля упругости кортикальной кости; сжатие

[Figure 6. The maximum stresses in the cortical bone vs the modulus of elasticity of the cortical bone; compression]

Рис. 7. Интенсивность напряжений Oi вдоль контура губчатой кости; сжатие, Ес = 18 ГПа: a) Oi,max = 2.6 МПа, Es = 0.5 ГПа; b) ffj,max = 3.9 МПа, Es = 5.0 ГПа

[Figure 7. Stresses intensity Oi along the spongy bone contour; compression, Ec = 18 GPa: a) Oi,max = 2.6 MPa, Ea = 0.5 GPa; b) Oi,max = 3.9 MPa, Ea = 5.0 GPa]

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

3. Результаты расчетов при наклонной нагрузке на коронку. Второй этап расчета на модели имплантата выполнен при наклонной нагрузке, приложенной к коронке (см. рис. 2). Зависимости максимальных эквивалентных напряжений в костных тканях от модулей упругости кортикальной и губчатой кости представлены на рис. 8—11. При наклонной нагрузке тенденции изменения максимальных эквивалентных напряжений аналогичны полученным ранее при сжимающей нагрузке, но уровень максимальных напряжений выше, что связано с действием сдвиговой компоненты нагрузки рх, приводящей к возникновению изгибающих моментов М\ ~ рхк в слое кортикальной

16

12

V /

//2 / 3 1-EC = 6 GPa 2 - Ec = 10 GPa 3 - Ec = 18 GPa

0 2 4 6 8

Е„ GPa

Рис. 8. Зависимость максимальных напряжений в губчатой кости от модуля упругости губчатой кости; наклонная нагрузка

[Figure 8. The maximum stresses in the spongy bone vs the modulus of elasticity of the spongy bone; inclined load]

80

60

40

20

1 - Ec = 6 GPa 2-Ec = 10 GPa 3 - Ec = 18 GPa

3

0

4

, GPa

6

Рис. 9. Зависимость максимальных напряжений в кортикальной кости от модуля упругости губчатой кости; наклонная нагрузка

[Figure 9. The maximum stresses in the cortical bone vs the modulus of elasticity of the spongy bone; inclined load]

12

10

я Он

3

1 — E„ = 0.5 GPa 2 — Es = 1.0 GPa 3 — Es = 3.0 GPa

2X.

16

18

4 8 12

Ес, GPa

Рис. 10. Зависимость максимальных напряжений в губчатой кости от модуля упругости кортикальной кости; наклонная нагрузка [Figure 10. The maximum stresses in the spongy bone vs the modulus of elasticity of the cortical bone; inclined load]

80

70

« Oh

S 60

50

40

2/

3 1-E,= 0.5 GPa 2-Es = 1.0 GPa 3-E„= 3.0 GPa

12

E,:, GPa

16

18

Рис. 11. Зависимость максимальных напряжений в кортикальной кости от модуля упругости кортикальной кости; наклонная нагрузка [Figure 11. The maximum stresses in the cortical bone vs the modulus of elasticity of the cortical bone; inclined load]

кости и М2 ~ рх(к + Ь) — у основания имплантата, где к — высота коронки, Ь — длина имплантата.

Увеличение модуля упругости губчатой костной ткани Е3 приводит к возрастанию максимальных напряжений в губчатой кости, причем рост напряжений происходит более интенсивно, чем при нормальной нагрузке, см. рис. 3 и 8. Отметим, что при наклонной нагрузке и малых значениях модуля упругости губчатой кости (Е3 < 2 МПа) изменение напряжений в кортикальной кости слабо зависит от модуля упругости кортикальной костной ткани (см. рис. 4 и 9). Напряжения в губчатой кости снижаются при увеличении модуля упругости кортикальной кости Ес, см. рис. 8 и 10. При этом значения максимальных напряжений значительно (примерно в три раза) превышают соответствующие значения при нормальной нагрузке (см. рис. 3, 5 и рис. 8, 10), что связано с действием изгибающего момента. Уровень напряжений в кортикальной кости возрастает при увеличении модуля упругости этой кости, (см. рис. 9 и 11). Для максимальных напряжений в кортикальной кости при наклонной нагрузке отличия от случая нормальной нагрузки невелики (см. рис. 4, 6 и рис. 9, 11), что связано с меньшей величиной механического момента М\ в слое кортикальной кости.

Распределение интенсивности напряжений по границе губчатой кости при действии наклонной нагрузки представлено на рис. 12. При этом характер распределения напряжений в губчатой кости заметно изменяется: при Е8 = 0.5 ГПа максимальные напряжения наблюдаются в нижней части зоны контакта губчатой кости и основания имплантата, а при Е3 = 5.0 ГПа положение максимума напряжений смещается в зону шейки имплантата на границу с губчатой и кортикальной костными тканями.

Рис. 12. Интенсивность напряжений Oi вдоль контура губчатой кости; сжатие, Ес=18 ГПа: a) CTi,max = 4.3 МПа, Еа = 0.5 ГПа; b) aitmax = 12.2 МПа, Еа = 5.0 ГПа [Figure 12. Stresses intensity Oi along the spongy bone contour; compression, Ec = 18 GPa: a) ffi,max = 4.3 MPa, Ea = 0.5 GPa; b) aitmax = 12.2 MPa, Ea = 5.0 GPa]

Заключение. Исследование зависимости интенсивности напряжений от плотности кортикальной и губчатой костных тканей выполнено посредством параметрической вариации модулей упругости соответствующих костных тканей для модели соединения цилиндрического имплантата и костных тканей при действии нормальной и наклонной нагрузок. Полученные результаты качественно (ввиду различия геометрических параметров моделей и механических свойств) согласуются с результатами работ [3,8,11].

При увеличении модуля упругости губчатой костной ткани (увеличении ее плотности) возрастает часть нагрузки, передаваемая на эту костную ткань. Поэтому максимальные эквивалентные напряжения в губчатой костной ткани возрастают, что стимулирует процесс остеоинтеграции имплантата, сопровождающийся ростом плотности костной ткани [23]. Напряжения в кортикальной костной ткани при увеличении модуля упругости губчатой кости снижаются за счет уменьшения нагрузки, передаваемой на эту часть кости, что приводит к более однородному распределению напряжений в костных тканях при нагружении имплантата. С другой стороны, наличие слоя кортикальной кости высокой плотности (с большим модулем упругости) приводит к снижению напряжений в губчатой кости, что позволяет повысить надежность имплантации при ослабленной губчатой кости. При этом уровень максимальных напряжений в кортикальной кости возрастает и область максимальных напряжений наблюдается вблизи шейки имплантата.

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

Конкурирующие интересы. Конкурирующих интересов не имею. Авторская ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.

Финансирование. Работа выполнена в рамках государственного задания (№ госрегистрации 123021700050-1).

Библиографический список

1. Параскевич В. Л. Дентальная имплантология: Основы теории и практики. M.: МИА, 2011. 400 с.

2. Chugh T., Jain A. K., Jaiswal R. K., et al. Bone density and its importance in orthodontics // J. Oral Biol. Craniofac. Res., 2013. vol. 3, no. 2. pp. 92-97. DOI: https://doi.org/10.1016/ j.jobcr.2013.01.001.

3. Premnath K., Sridevi J., Kalavathy N., et al. Evaluation of stress distribution in bone of different densities using different implant designs: A three-dimensional finite element analysis// J. Indian Prosthodont Soc., 2013. vol. 13, no. 4. pp. 555-559. DOI: https://doi. org/10.1007/s13191-012-0189-7.

4. Wirth A. J., Muller R., van Lenthe G. H. Computational analyses of small endosseous implants in osteoporotic bone// Eur. Cell. Mater., 2010. vol.20. pp. 58-71. DOI: https:// doi.org/10.22203/ecm.v020a06.

5. Lee H., Jo M., Noh G. Biomechanical effects of dental implant diameter, connection type, and bone density on microgap formation and fatigue failure: A finite element analysis // Comput. Methods Programs Biomed., 2021. vol.200, 105863. DOI: https://doi.org/10. 1016/j.cmpb.2020.105863.

6. Zioupos P., Cook R. B., Hutchinson J. R. Some basic relationships between density values in cancellous and cortical bone// J. Biomech., 2008. vol.41, no. 9. pp. 1961-1968. DOI: https://doi.org/10.1016/j.jbiomech.2008.03.025.

7. Clift S. E., Fisher J., Watson C. J. Finite element stress and strain analysis of the bone surrounding a dental implant: Effect of variations in bone modulus // Proc. Inst. Mech. Eng. H, 1992. vol.206, no. 4. pp. 233-241. DOI: https://doi.org/10.1243/PIME_PROC_1992_206_ 295_02.

8. Рогожников Г. И., Конюхова С. Г., Няшин Ю. И. [и др.] Влияние модуля упругости губчатой и кортикальной кости на напряженное состояние в области пластинчатого

имплантата при окклюзионной нагрузке // Российский журнал биомеханики, 2004. Т. 8, №1. С. 54-60. EDN: JWSHCL.

9. Nutu E., Ahmad S., Pastrama S. Influence of bone elastic properties on the predicted stress distribution in the dental implant vicinity // Materials Today: Proceedings, 2017. vol. 4, no. 5, part 1. pp. 5904-5908. DOI: https://doi.org/10.1016/j.matpr.2017.06.067.

10. Олесова В. Н., Бронштейн Д. А., Лернер А. Я. [и др.] Напряженно-деформированное состояние в протезной конструкции на дентальном имплантате при цементной фиксации искусственной коронки// Российский журнал биомеханики, 2016. Т. 20, №4. С. 311-315. EDN: XXMNXF. DOI: https://doi.org/10.15593/RJBiomech/2016.4.02.

11. Fedorova N. V. The study of the stress-strain state of the dental ceramic implants depending on their shape and bone mineralization degree // Russian Journal of Biomechanics, 2019. vol.23, no. 3. pp. 388-394. EDN: BDBUKN. DOI: https://doi.org/10.15593/ RZhBiomeh/2019.3.10.

12. Дьяченко Д. Ю., Дьяченко С. В. Применение метода конечных элементов в компьютерной симуляции для улучшения качества лечения пациентов в стоматологии: систематический обзор// Кубанский научный медицинский вестник, 2021. Т. 28, №5. С. 98-116. EDN: KDCHLT. DOI: https://doi.org/10.25207/1608-6228-2021-28-5-98-116.

13. Buyuk F. N., Savran E., Karpat F. Review on finite element analysis of dental implants // J. Dent. Implant Res., 2022. vol.41, no. 3. pp. 50-63. DOI: https://doi.org/10.54527/ jdir.2022.41.3.50.

14. Wolfe L. A. Stress analysis of endosseous implants using the Boundary Integral Equation (BIE) method // J. Biomed. Eng., 1993. vol. 15, no. 4. pp. 319-323. DOI: https://doi.org/ 10.1016/0141-5425(93)90009-N.

15. Перельмутер М. Н. Исследование напряженно-деформированного состояния стоматологических имплантатов методом граничных интегральных уравнений // Вестник ПНИПУ. Механика, 2018. №2. С. 83-95. EDN: XUGGCL. DOI: https://doi.org/10.15593/ perm.mech/2018.2.08.

16. Citarella R., Armentani E., Caputo F., Lepore M. Stress analysis of an endosseus dental implant by BEM and FEM // The Open Mechanical Engineering Journal, 2012. vol. 6. pp. 115-124. DOI: https://doi.org/10.2174/1874155X01206010115.

17. Misch C. E., Qu Z., Bidez M W. Mechanical properties of trabecular bone in the human mandible: Implications for dental implant treatment planning and surgical placement // J. Oral Maxillofac. Surg., 1999. vol.57, no. 6. pp. 700-706. DOI: https://doi.org/10.1016/ S0278-2391(99)90437-8.

18. Король Д. М., Николов В. В., Онипко Е. Л., Ефименко А. С. Определение интенсивности окклюзионного давления у пациентов на ортопедическом приеме // Современная медицина: актуальные вопросы, 2015. №46-47. С. 40-46. EDN: UINWJH.

19. Banerjee P. K., Butterfield R. Boundary Element Methods in Engineering Science. London: McGraw-Hill, 1981. 452 pp.

20. Перельмутер М. Н. Применение метода граничных элементов при исследовании пространственного напряженного состояния составных конструкций / Проблемы прочности и динамики в авиадвигателестроении. Вып. 4 / Труды ЦИАМ, Т. 1237, 1989. С. 74-99.

21. Perelmuter M. Boundary element analysis of structures with bridged interfacial cracks // Comput. Mech., 2013. vol.51, no. 4. pp. 523-534. EDN: RFJHDN. DOI: https://doi.org/10. 1007/s00466-012-0817-4.

22. Perelmuter M. Analysis of interaction of bridged cracks and weak interfaces // Int. J. Mech. Sci., 2018. vol.149, no. 4. pp. 349-360. EDN: BLWNCV. DOI: https://doi.org/10.1016/j. ijmecsci.2018.10.011.

23. Lin D., Li Q., Li W., et al. Mandibular bone remodeling induced by dental implant // J. Biomech., 2010. vol.43, no. 2. pp. 287-293. DOI: https://doi.org/10.1016/j.jbiomech. 2009.08.024.

Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki

[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 1, pp. 189-201 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1976

MSC: 74A10, 74S15, 74M15

The effect of bone tissue density on the stress-strain state near dental implants

M. N. Perelmuter

Ishlinsky Institite for Problems in Mechanics, Russian Academy of Sciences,

101—1, pr. Vernadskogo, Moscow, 119526, Russian Federation.

Abstract

The dependence of the stress-strain state of the bone tissue on its density near the dental implant has been studied. The computations were performed by the boundary element method for the plane-deformed state of a model consisting of a cylindrical implant and surrounding bone tissues. Bone tissue is considered as an isotropic and homogeneous elastic material. Simulation the effect of bone density on the stress-strain state when applying a quasi-static load is performed by changing of elasticity modulus of the bone. It has been established that with the increasing in the spongy bone tissue elastic modulus, the maximum equivalent stresses in this bone tissue increase. Stresses in the cortical bone tissue decrease with the increasing in the spongy bone elastic modulus due to the decreasing in the load transferred to this bone part. Stresses in the spongy bone decrease with the increasing in the cortical bone layer elasticity modulus. The level of maximum stress in the cortical layer of the bone increases with the increasing of this bone tissue elastic modulus. The maximum of stresses in the cortical bone tissue are observed near the implant neck.

Keywords: implant, bone tissues, stress-strain state, boundary integral equations method.

Received: 26th December, 2022 / Revised: 19th February, 2023 / Accepted: 17th March, 2023 / First online: 24th March, 2023

Competing interests. I have no competing interests.

Author's Responsibilities. I take full responsibility for submitting the final manuscript in print. I approved the final version of the manuscript.

Mathematical Modeling, Numerical Methods and Software Complexes Short Communication

© Authors, 2023

© Samara State Technical University, 2023 (Compilation, Design, and Layout) 9 © The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:

Perelmuter M. N. The effect of bone tissue density on the stress-strain state near dental implants, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 1, pp. 189-201. EDN: VDVFUD. DOI: 10.14498/vsgtu1976 (In Russian). Author's Details:

Mikhail N. Perelmuter https://orcid.org/0000-0002-8430-5412

Dr. Phys. & Math. Sci.; Leading Researcher; Lab. of Mechanics of Strength and Fracture of Materials and Structures; e-mail: perelm@ipmnet.ru

Funding. The work was carried out within the framework of a state assignment (state registration no. 123021700050-1.

References

1. Paraskevich V. L. Dental'naia implantologiia: Osnovy teorii i praktiki [Dental Implantology: The fundamentals of Theory and Practice]. Moscow, Medical Information Agency, 2011, 400 pp.

2. Chugh T., Jain A. K., Jaiswal R. K., et al. Bone density and its importance in orthodontics, J. OralBiol. Craniofac. Res., 2013, vol. 3, no. 2, pp. 92-97. DOI: https://doi.org/10.1016/ j.jobcr.2013.01.001.

3. Premnath K., Sridevi J., Kalavathy N., et al. Evaluation of stress distribution in bone of different densities using different implant designs: A three-dimensional finite element analysis, J. Indian Prosthodont Soc., 2013, vol.13, no. 4, pp. 555-559. DOI: https://doi. org/10.1007/s13191-012-0189-7.

4. Wirth A. J., Muller R., van Lenthe G. H. Computational analyses of small endosseous implants in osteoporotic bone, Eur. Cell. Mater., 2010, vol.20, pp. 58-71. DOI:https:// doi.org/10.22203/ecm.v020a06.

5. Lee H., Jo M., Noh G. Biomechanical effects of dental implant diameter, connection type, and bone density on microgap formation and fatigue failure: A finite element analysis, Comput. Methods Programs Biomed., 2021, vol.200, 105863. DOI: https://doi.org/10. 1016/j.cmpb.2020.105863.

6. Zioupos P., Cook R. B., Hutchinson J. R. Some basic relationships between density values in cancellous and cortical bone, J. Biomech., 2008, vol. 41, no. 9, pp. 1961-1968. DOI: https:// doi.org/10.1016/j.jbiomech.2008.03.025.

7. Clift S. E., Fisher J., Watson C. J. Finite element stress and strain analysis of the bone surrounding a dental implant: Effect of variations in bone modulus, Proc. Inst. Mech. Eng. H, 1992, vol.206, no. 4, pp. 233-241. DOI: https://doi.org/10.1243/PIME_PROC_1992_206_ 295_02.

8. Rogozhnikov G. I., Konyuhova S. G., Nyashin Y. I., et al. The influence of elasticity modulus of spongy and cortical bone on stress state near planar implant under occlusal load, Russian Journal of Biomechanics, 2004, vol. 8, no. 1, pp. 49-54.

9. Nutu E., Ahmad S., Pastrama S. Influence of bone elastic properties on the predicted stress distribution in the dental implant vicinity, Materials Today: Proceedings, 2017, vol. 4, no. 5, part 1, pp. 5904-5908. DOI: https://doi.org/10.1016Zj.matpr.2017.06.067.

10. Olesova V. N., Bronshtein D. A., Lerner A. Ya., et al. Stress-strain state in prosthetic construction on dental implant with cement fixing artificial crown, Russian Journal of Biomechanics, 2016, vol.20, no. 4, pp. 266-269.

11. Fedorova N. V. The study of the stress-strain state of the dental ceramic implants depending on their shape and bone mineralization degree, Russian Journal of Biome-chanics, 2019, vol.23, no. 3, pp. 388-394. EDN: BDBUKN. DOI: https://doi.org/10.15593/ RZhBiomeh/2019.3.10.

12. Dyachenko D. Yu., Dyachenko S. V. Finite element method in computer simulation for improved patient care in dentistry: A systematic review, Kuban Scientific Medical Bulletin, 2021, vol. 28, no. 5, pp. 98-116 (InRussian). EDN: KDCHLT. DOI: https://doi.org/10.25207/ 1608-6228-2021-28-5-98-116.

13. Büyük F. N., Savran E., Karpat F. Review on finite element analysis of dental implants, J. Dent. Implant Res., 2022, vol.41, no. 3, pp. 50-63. DOI: https://doi.org/10.54527/ jdir.2022.41.3.50.

14. Wolfe L. A. Stress analysis of endosseous implants using the Boundary Integral Equation (BIE) method, J. Biomed. Eng., 1993, vol.15, no. 4, pp. 319-323. DOI: https://doi.org/ 10.1016/0141-5425(93)90009-N.

15. Perelmuter M. N. Analysis of stress-strain state of dental implants by the boundary integral equations method, PNRPU Mechanics Bulletin, 2018, no. 2, pp. 83-95 (In Russian). EDN: XUGGCL. DOI: https://doi.org/10.15593/perm.mech/2018.2.08.

16. Citarella R., Armentani E., Caputo F., Lepore M. Stress analysis of an endosseus dental implant by BEM and FEM, The Open Mechanical Engineering Journal, 2012, vol. 6, pp. 115124. DOI: https://doi.org/10.2174/1874155X01206010115.

17. Misch C. E., Qu Z., Bidez M W. Mechanical properties of trabecular bone in the human mandible: Implications for dental implant treatment planning and surgical placement, J. Oral Maxillofac. Surg., 1999, vol.57, no. 6, pp. 700-706. DOI: https://doi.org/ 10.1016/S0278-2391(99)90437-8.

18. Korol D. M., Nikolov V. V., Onipko E. L., Efimenko A. S. Determination of the intensity of occlusal pressure in patients at orthopedic examination, Modern Medicine: Current Issues, 2015, no. 46-47, pp. 40-46 (In Russian). EDN: UINWJH.

19. Banerjee P. K., Butterfield R. Boundary Element Methods in Engineering Science. London, McGraw-Hill, 1981, 452 pp.

20. Perelmuter M. N. Application of the boundary element method in the study of the spatial stress state of composite structures, In: Problems of Strength and Dynamics in Aircraft Engine. Issue 4, Proc. CIAM, 1237, 1989, pp. 74-99 (In Russian).

21. Perelmuter M. Boundary element analysis of structures with bridged interfacial cracks, Comput. Mech., 2013, vol.51, no. 4, pp. 523-534. EDN: RFJHDN. DOI: https://doi.org/10. 1007/s00466-012-0817-4.

22. Perelmuter M. Analysis of interaction of bridged cracks and weak interfaces, Int. J. Mech. Sci., 2018, vol.149, no. 4, pp. 349-360. EDN: BLWNCV. DOI: https://doi.org/10.1016/j. ijmecsci.2018.10.011.

23. Lin D., Li Q., Li W., et al. Mandibular bone remodeling induced by dental implant, J. Biomech., 2010, vol.43, no. 2, pp. 287-293. DOI: https://doi.org/10.1016/j.jbiomech. 2009.08.024.