Определение силовых факторов, развиваемых ортодонтическим аппаратом для расширения верхней челюсти Текст научной статьи по специальности «Физика»

УДК 539.3+(616.314-089.23)

Российский Журнал

www.biomech.ru

ОПРЕДЕЛЕНИЕ СИЛОВЫХ ФАКТОРОВ, РАЗВИВАЕМЫХ ОРТОДОНТИЧЕСКИМ АППАРАТОМ ДЛЯ РАСШИРЕНИЯ ВЕРХНЕЙ ЧЕЛЮСТИ

Г.И. Михасев, С.М. Босяков, К.С. Юркевич

Белорусский государственный университет, Республика Беларусь, 220030, Минск, проспект

Независимости, 4, e-mail: bosiakov@bsu.by

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

Ключевые слова: ортодонтический аппарат, расширение верхней челюсти, математическое моделирование, конечно-элементный расчет, реакции и реактивные моменты.

Введение

Одним из факторов успешного лечения аномалий зубочелюстной системы является исследование биомеханического воздействия ортодонтических аппаратов на перемещаемые зубы [4]. Обзор некоторых теоретических и экспериментальных результатов, полученных в этом направлении, представлен, в частности, в монографии [2]. Здесь указано на важность корректного дозирования нагрузки, создаваемой ортодонтическими аппаратами, поскольку превышение развиваемых ими усилий или их некорректное приложение может привести к патологическим процессам, возникающим в периодонте и пародонте при перемещении зубов. Кроме этого, при дозировании нагрузки внимание в основном уделяется нагрузке, действующей со стороны стержня аппарата на корень зуба и направленной вдоль осей координат, поскольку предполагается корпусное перемещение зуба. Поэтому важно корректно выбрать углы наклона стержней аппарата к определенным осям координат, а также точку закрепления стержня на коронке ортодонтического аппарата с тем, чтобы по возможности избежать появления моментов сил, поворачивающих зуб. В настоящей работе представлены результаты анализа сосредоточенных сил и моментов сил, создаваемых ортодонтическим аппаратом с четырьмя стержнями и винтом, применяемым для лечения аномалий верхней челюсти.

© Михасев Г.И., Босяков С.М., Юркевич К.С., 2010

Михасев Геннадий Иванович, д.ф.-м.н., проф., завкафедрой био- и наномеханики БГУ, Минск Босяков Сергей Михайлович, д.ф.-м.н., доцент кафедры теоретической и прикладной механики БГУ, Минск

Юркевич Кирилл Сергеевич, аспирант кафедры теоретической и прикладной механики БГУ, Минск

09806267

Реактивные силы и моменты сил

Верхняя челюсть с установленным ортодонтическим аппаратом представлена на рис. 1. Аппарат развивает усилия, приводящие к перемещению зубов, за счет раскручивания винта, соединяющего пластинки, к которым прикреплены стержни с коронками, жестко закрепленные на зубах. Верхний конец стержня, соединенный с пластинкой ортодонтического аппарата, получает перемещения, равные смещению пластинки аппарата. Далее рассмотрим статическую модель, согласно которой нижний конец стержня соединенный с коронкой, жестко заделан и при перемещении верхнего конца стержня создает в точке закрепления реактивные силы и моменты, под действием которых по истечении определенного промежутка времени начинается перемещение зубов. Такое допущение возможно, поскольку одновременно с раскручиванием винта перемещения зубов не происходит. Расчетная схема для стержня АВ ортодонтического аппарата представлена на рис. 2 (точка А соответствует верхнему концу стержня, точка В - нижнему концу стержня).

На рис. 2 точка А, совпадающая с началом координат системы 0X72, в которой верхний конец стержня прикреплен к пластинке ортодонтического аппарата, при раскручивании винта получает перемещение и0, направленное вдоль координатной оси 07. Для определения реактивных силы и момента силы перейдем в координатную плоскость х 0 у, проходящую через ось 07 и сам стержень.

Перемещение иу вдоль оси 0у найдем из уравнения поперечного изгиба

стержня

3 4uy 3 Ч

EI-У - Rx—2у = 0, (1)

4 л -ч 2

3х 3х

р у 1 4

где E - модуль упругости; I = ~рг - осевой момент инерции поперечного сечения

стержня; r - радиус поперечного сечения стержня; Rx = -1 Esu0х - реакция коронки,

направленная вдоль оси 0х; u0х = u0 cos (a7 ) - полное перемещение точки A вдоль оси 0х; s - площадь поперечного сечения стержня; l - длина стержня.

Рис. 1. Верхняя челюсть с установленным ортодонтическим аппаратом: 1 - пластинки; 2 - винт; 3 - стержень; 4 - коронки; 5, 6 - точки, совпадающие с положением точек А и В на расчетной схеме для стержня ортодонтического аппарата

С учетом того, что точка B жестко закреплена на зубе пациента, граничные условия представим в следующем виде:

( bMy ( 0 ) = duy ( l ) = 0

() Эх Эх '

(2)

Здесь u0 = u0 sin(aY) - полное перемещение точки A вдоль оси 0y. Из уравнения (1) при граничных условиях (2) получим

uy =

*0 y

2 (Ф cos (Ф)- sin (Ф))

х

х

2Ф

Ф

(3)

— (l - x) cos (Ф)-sin (Ф)-sin I — (l - 2x)

1 К s где Ф = — < .

2\ Л

Из выражения (3) с учетом дифференциальных зависимостей между кривизной оси стержня, сосредоточенной силой и изгибающим моментом будем иметь следующие выражения для реактивных силы и момента сил, возникающих в зубе пациента и действующих в плоскости x0у :

R

Esu,

0 x

1 +

u0xu0ys COs

(Ф)

M

4 Il (Ф cos (Ф) - sin (Ф))

Eu0xu0ys sin (Ф)

l (ф cos (ф) - sin (ф))2

2 ’

(4)

(5)

l

Для определения усилий, развиваемых ортодонтическим аппаратом в медио-дистальном (касательном к зубному ряду), вестибуло-оральном (перпендикулярном касательной к зубному ряду) и вертикальном направлениях, а также моментов сил, возникающих относительно этих направлений, найдем проекции реакции Я на координатные оси системы 0ХТ2 и реактивные моменты сил, действующие относительно осей этой системы:

Здесь аХ, а7 и а2 - углы между стержнем и осями координат 0Х, 0У и 02 соответственно. Тогда для реакций ЯтЛ (для медио-дистального направления), Ят (для вестибуло-орального направления), Яу (для вертикального направления) и моментов сил , Ыт, Ыу, действующих относительно соответствующих направлений, будем

иметь

Проведем расчет сил и моментов сил, развиваемых ортодонтическим аппаратом при раскручивании винта в медио-дистальном, вестибуло-оральном и вертикальном направлениях. Будем считать, что расстояние, на которое перемещаются пластинки аппарата и0 = 0,1 мм. Модуль упругости материала, из которого изготовлены стержни, Е = 200 ГПа, радиус стержня г = 1 мм, длина стержня I = 2,5 см. На рис. 3 и 4 представлены зависимости реакций, возникающих в медио-дистальном, вестибуло-оральном и вертикальном направлениях при расположении стержня в плоскости 702 ( аХ = 90°) и в плоскости, составляющей угол я/4 с плоскостями Х02 и У02, от угла а7 и угла а2 соответственно. При построении принимаем угол 82 = 30°. Значения реакций по оси ординат приведены в кН.

R cos (aX)

^cos2 (aX) + cos2 (aZ) + sin2 (aY) cot2 (aY + )

Ry = -R cos (aY + 5j),

Rz

Rcos(aZ)

(6)

, MY = 0,

(7)

Rmd = RX cos (d2 )-RY sin (d2 ), Rvo = RX sin (d2 ) + RY cos (d2 ), Rv = RZ ,

Mmd = MX C0S (d2 ), Mvo = MX Sin (d2 ) > Mv = MZ >

где 82 - угол между касательной к зубному ряду и координатной осью 0X.

Влияние геометрических параметров

R, кН

О 0,25 0,5 0,75 1 1,25 ау,рад

Рис. 3. Зависимости реакций от угла aY при расположении стержня в плоскости Y0Z :

0 0,25 0,5 0,75 1 1,25 а7,рад

Рис. 4. Зависимости реакций от угла а2 при расположении стержня в плоскости, составляющей угол к 4 с плоскостями X02 и У02 : 1 - Я ,; 2 - Я ; 3 - Я

I та ^ уо 5 V

Из рис. 3 видно, что при расположении стержня в плоскости У02 усилия, развиваемые аппаратом в медио-дистальном и вестибуло-оральном направлениях, максимальны при аУ = 0, причем Яуо > Ята. С увеличением угла аУ реакции Ята и

Яуо уменьшаются до нулевого значения. Реакция в вертикальном направлении при

аУ = 0 и аУ = 90° равна нулю; максимальное значение реакции Яу наблюдается при

аУ = 45°.

Как следует из рис. 4, при расположении стержня в плоскости, составляющей угол к/4 с плоскостями Х02 и У02, с увеличением угла а2 усилия, развиваемые аппаратом в медио-дистальном и вестибуло-оральном направлениях возрастают по абсолютной величине от нуля до максимальных значений. Зависимость реакции Яу от

угла а2 аналогична предыдущему случаю.

При увеличении радиуса поперечного сечения стержня усилия, развиваемые ортодонтическим аппаратом по всем направлениям, резко возрастают.

М, Н м

О

-0,01 -0,02 -0,03 -0,04 -0,05 -0,06

0 0,25 0,5 0,75 1 1,25 ау, рад

Рис. 5. Зависимости реактивных моментов от угла аУ при расположении стержня в плоскости У 0 2 : 1 - М ,; 2 - М

та 5 уо

М, Н м

0

-0,01 -0,02 -0,03 -0,04 -0,05 -0,06

0 0,25 0,5 0,75 1 1,25 ау, рад

Рис. 6. Зависимости реактивных моментов от угла а2 при расположении стержня в плоскости, составляющей угол к 4 с плоскостями X02 и У02 : 1 - М ,; 2 - М ;

-1 / та уо

3 - М

у

Визуализация зависимостей для реактивных моментов сил, возникающих относительно медио-дистального и вестибуло-орального направлений, а также относительно оси 02, в двух рассматриваемых плоскостях выполнена на рис. 5 и 6 (значения моментов сил по оси ординат приведены в Нм). При построении приняты прежние числовые данные.

Из рис. 5 видно, что при расположении стержня в плоскости У02 при аУ = 0 все реактивные моменты равны нулю. С увеличением аУ моменты сил Мта и Муо возрастают по абсолютной величине. Момент сил Му равен нулю. В случае, если стержень располагается в плоскости, составляющей угол к/ 4 с плоскостями Х02 и У02, отличными от нуля являются все три реактивных момента.

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

Конечно-элементный анализ реакций и реактивных моментов

Выполним конечно-элементный расчет реакций и реактивных моментов, возникающих в точке закрепления стержня на коронке аппарата, с применением функциональных возможностей пакета ANSYS. Г еометрическая модель ортодонтического аппарата, построенная в CAD-пакете CATLA V5 для случая, если углы наклона стержней к координатным осям равны между собой, длина стержня составляет 25 мм, радиус стержня равен 1 мм, представлена на рис. 7.

Один из концов стержня закреплен на насадке аппарата, которая надевается на зуб. Другой конец каждого из стержней закреплен на одной из граней двух пластинок, которые моделировались в форме параллелепипедов. Далее рассмотрим один из стержней модели ортодонтического аппарата, считая корень зуба жестко заделанным в костной ткани в момент раскручивания винта ортодонтического аппарата. Конец стержня, закрепленный на пластинке, получает перемещение, равное 0,1 мм, в направлении оси абсцисс. Расчетная модель с визуализацией положения стержня до и после перемещения одного из его концов показана на рис. 8. Упругие свойства стержня описываются константами E = 200 ГПа, V = 0,3; радиус поперечного сечения стержня равен 1 мм, длина стержня составляет 25 мм.

В ходе проведения вычислительного эксперимента стержень моделировался линейными пространственными балочными конечными элементами Beam188 и Beam4. Элемент Beam 188 построен на основе балки Тимошенко, имеет шесть степеней свободы и пригоден для линейных, а также нелинейных задач с большими поворотами и деформациями. В качестве входных данных элемента используются геометрические параметры его поперечного сечения. Балочный конечный элемент Beam4 является элементом с одной осью и имеет возможность воспринимать растяжение, сжатие кручение и изгиб. Так же как и элемент Beam 188, Beam4 имеет шесть степеней свободы в каждом узле. Элемент определяется двумя узлами, площадью поперечного сечения, двумя моментами инерции и двумя толщинами [1]. В нашем случае использовался балочный конечный элемент Beam4.

Рис. 8. Положение стержня ортодонтического аппарата до и после перемещения

пластинки

В результате конечно-элементного расчета получены значения реакций и реактивных моментов, возникающих в точке закрепления стержня аппарата на коронке, для различных углов наклона стержня по отношению к осям декартовой системы координат. На рис. 9 и 10 выполнена одновременная визуализация значений реакций, найденных в пакете ЛЫ8У8 для определенных углов наклона стержня а7 и а2 в плоскости 702 и в плоскости, составляющей угол я /4 с плоскостями Х02 и У02 (обозначены точками), и зависимости (6) от угла аг и а2 в тех же плоскостях. При построении зависимостей (6) для стержня принимались числовые данные, используемые при конечно-элементном расчете.

Из рис. 9 и 10 видно, что результаты конечно-элементного расчета и вычисления реакций на основании решения уравнения (4) практически совпадают. Расхождения в некоторых случаях наблюдаются при предельных углах наклона стержня аппарата к осям координат, равных нулю и я / 2. Значения реакций в этих случаях могут отличаться в 4-7 раз. Для остальных углов наклона стержня отличие результатов конечно-элементного расчета при использовании элемента Веат4 по сравнению с вычислением реакций на основании аналитических выражений не превышало 0,5 %. Также отметим, что при проведении аналогичных расчетов с применением балочного элемента Веат 188 расхождение результатов не превышало 1,5%. В связи с этим

можно сделать вывод о том, что для получения более корректных результатов целесообразнее использовать балочный конечный элемент Веат4.

Я, кН

О 0,25 0,5 0,75 1 1,25 ау,рад

Рис. 9. Зависимости реакций от угла ау при расположении стержня в плоскости У02 (сплошные линии) и значения реакций, полученные в ходе конечно-элементного расчета (обозначены точками): 1 - Ку; 2 - Кг

К, кН

0 0,25 0,5 0,75 1 1,25 а2,рад

Рис. 10. Зависимости реактивных моментов от угла а2 при расположении стержня

в плоскости, составляющей угол я/4 с плоскостями X02 и У02 (сплошные линии),

и значения реакций, полученные в ходе конечно-элементного расчета (обозначены

точками): 1 - К = К ; 2 - Я

' X у I

Значения реактивных моментов для определенных углов наклона стержня аУ и а2 в плоскости У02 и в плоскости, составляющей угол я /4 с плоскостями Х02 и У02, и зависимости (7) от угла аУ и а2 в тех же плоскостях приведены на рис. 11 и 12. При построении зависимостей (7) для стержня используются те же числовые данные, что и при конечно-элементном расчете.

Максимальное отличие результатов конечно-элементного расчета реактивных моментов по сравнению со значениями, найденными в соответствии с выражениями (7), составляет в первом случае 9,8%, во втором - 12,5%. Приемлемая величина погрешности конечно-элементного расчета указывает на пригодность использования

Мх, Н-м О

-0,01 -0,02 -0,03 -0,04 -0,05 -0,06 -0,07

0 0,25 0,5 0,75 1 1,25 ау, рад

Рис. 11. Зависимость реактивного момента Мх от угла ау при расположении стержня

в плоскости У02 (сплошная линия) и значения реактивного момента, полученные в ходе конечно-элементного расчета (обозначены точками)

М, Н-м

о

-0,01 -0,02 -0,03 -0,04 -0,05 -0,06 -0,07

0 0,25 0,5 0,75 1 1,25 а2,рад

Рис. 12. Зависимости реактивных моментов от угла а2 при расположении стержня

в плоскости, составляющей угол к/4 с плоскостями X02 и У02 (сплошные линии),

и значения реакций, полученные в ходе конечно-элементного расчета (обозначены

точками): 1 - М , 2 - М

' х ’ I

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

Заключение

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

2. Также следует отметить, что не следует пренебрегать моментами сил, действующих на зуб со стороны стержня ортодонтического аппарата, поскольку из-за малой жесткости периодонта при поворотах корня [5] действие даже незначительных моментов сил может привести к появлению достаточно высоких углов поворота зуба относительно первоначального положения. Это тем более следует при атрофии костной ткани, приводящей к существенному снижению жесткостей периодонта [3].

3. Также отметим, что приемлемое совпадение результатов конечноэлементного расчета реакций и результатов, полученных на основании соотношений (6), (7), показывает, что полученные формулы для нахождения усилий, развиваемых ортодонтическим аппаратом для расширения верхней челюсти, могут быть непосредственно использованы на практике стоматологами-ортодонтами для корректного прогнозирования нагрузки, приходящейся на перемещаемые зубы.

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

Благодарности

Работа выполнена при поддержке Государственного комитета по науке и технологиям Республики Беларусь (проект «Разработать методику

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

Список литературы

1. Басов К.А. ANSYS: справочник пользователя. - М.: ДМК Пресс, 2005. - 640 с.

2. Маилян П.Д. Новые средства ортодонтического лечения. - Коломна: Коломенская типография, 1998. - 124 с.

3. Наумович С.А., Крушевский А.Е. Биомеханика системы зуб-периодонт. - Минск: Экономические технологии, 2000. - 132 с.

4. Хорошилкина Ф.Я., Малыгин Ю.М. Основы конструирования и технология изготовления ортодонтических аппаратов. - М.: Медицина, 1977.

5. Чуйко А.Н., Вовк В.Е. Особенности биомеханики в стоматологии. - Харьков: Прапор, 2006. - 304 с.

DETERMINATION OF THE POWER FACTORS DEVELOPED BY THE ORTHODONTIC APPLIANCE FOR THE MAXILLA EXTENSION

G.I. Mikhasev, S.M. Bosiakov, K.S. Yurkevich (Minsk, Belarus)

In present paper, results of calculation of the force factors acting on a tooth on the side of a rod of the orthodontic appliance, intended for expansion of the upper jaw are submitted. Numerical calculation of reactive forces and moments of forces for arrangements of the rod of the appliance in two coordinate planes are executed. The obtained results can be used at

designing constructions of orthodontic appliances in view of specific features of patients, and also for calculations of displacements of the fang during orthodontic treatment and the pressure arising in the periodontium and the bone tissue.

Key words: orthodontic appliance, extension of the upper jaw, mathematical modelling, finite element calculation, reactions and reaction moments.

Получено 23 декабря 2009