Область сопротивления зуба: экспериментальное определение Текст научной статьи по специальности «Математика»

DOI: 10.15593/RZhBiomeh/2015.1.07 УДК 531/534:[57+61]

Российский

Журнал

Биомеханики

www.biomech.ru

ОБЛАСТЬ СОПРОТИВЛЕНИЯ ЗУБА: ЭКСПЕРИМЕНТАЛЬНОЕ ОПРЕДЕЛЕНИЕ

А.Л. Дубинин

Кафедра теоретической механики и биомеханики Пермского национального исследовательского политехнического университета, Россия, 614990, Пермь, Комсомольский проспект, 29, e-mail: [email protected]

Аннотация. Для изучения начального перемещения зуба было введено понятие центра сопротивления. Для данной точки определены следующие свойства: если приложенная система сил сводится к равнодействующей с линией действия, проходящей через центр сопротивления, то зуб начнет поступательное перемещение; если приложенная система сил сводится к паре сил с моментом, равным главному моменту, то зуб повернется вокруг данного центра. В дальнейшем было показано, что центр сопротивления существует при выполнении исключительных условий (наличие оси симметрии у корня зуба). Поэтому для случаев, когда данный центр не существует, было предложено новое понятие - «область сопротивления зуба». Такой областью называется область минимального диаметра, имеющая следующие свойства: всякая прямая поступательного воздействия проходит через эту область и через всякую точку этой области проходит прямая поступательного воздействия. Показано, что она может быть эллипсом, двумя точками и одной точкой. В настоящей работе создана методика нахождения области сопротивления для частного случая двумя способами: 1) по совокупности всех прямых поступательного воздействия; 2) аналитически. Экспериментально установлена и показана связь между видом области сопротивления и формой корня зуба, а также рассмотрен переход от одного вида области к другому при варьировании геометрических параметров тела. Исследованы некоторые свойства отдельных видов области. Приведен случай существования центра сопротивления, не рассмотренный ранее.

Ключевые слова: центр сопротивления зуба, область сопротивления зуба, зубочелюстная система, ортодонтическое перемещение зубов, экспериментальное исследование.

Куда и как нужно приложить силу, чтобы переместить зуб в нужное положение? Данный вопрос является определяющим при ортодонтическом лечении и довольно непросто дать на него ответ ввиду необходимости учета большого количества факторов, влияющих на данный процесс: реакция тканей, окружающих зуб, их индивидуальные механические свойства, сложные физиологические процессы, характеристики аппарата, оказывающего силовое воздействие, особенности его установки и прочее [1, 5, 11].

В настоящей статье изучается вопрос о начальном перемещении зуба. Принимается модель (характерная для большинства работ, в которых также изучается данная проблема [2-4, 6, 7, 9, 10, 12, 13]): зуб, окруженный периодонтальной связкой,

© Дубинин А.Л., 2014

Дубинин Алексей Лаврентьевич, аспирант кафедры теоретической механики и биомеханики, Пермь

Введение

рассматривается как абсолютно твердое тело, погруженное в линейно-упругую среду и находящееся в состоянии устойчивого равновесия. При приложении системы сил тело испытывает малые перемещения в пределах костной лунки. Для изучения начального перемещения зуба часто пользуются понятием центра сопротивления, введенным в 1917 г. [8]. Данный центр является точкой, такой, что через нее проходят все прямые поступательного воздействия, и любая прямая, проходящая через данную точку, является прямой поступательного воздействия [2]. На его положение не влияют ни величина, ни направление приложенной к зубу силы, однако оно зависит от формы корня и упругих свойств окружающей среды. Известно также, что на условия существования центра сопротивления наложены довольно сильные ограничения, такие как наличие у корня зуба оси симметрии [1]. То есть данным понятием можно пользоваться, если рассматривать зуб в виде идеализированного тела вращения. Однако для более глубокого изучения вопроса необходима новая идея.

В работе [2] введено соответствующее понятие - «область сопротивления зуба», которое является обобщающим для центра сопротивления. Такой областью называется область минимального диаметра, имеющая следующие свойства: всякая прямая поступательного воздействия проходит через эту область и через всякую точку этой области проходит прямая поступательного воздействия. Она существует всегда и сохраняет основные свойства центра сопротивления. Показано, что она может быть набором точек, образующих эллипс, двумя точками или одной точкой (рис. 1).

Известно, что вид области сопротивления, представленный на рис. 1, а, характерен для общего случая задания формы корня зуба. Каждая из возможных прямых поступательного воздействия обязательно проходит через какую-либо точку области сопротивления (эллипс минимального диаметра). Случаи, представленные на рис. 1, б, в, соответствуют наличию элементов симметрии у тела. Видно, что на рис. 1, в область сопротивления состоит из одной точки, через которую проходят все прямые поступательного воздействия, что соответствует центру сопротивления.

Стоит отметить, что область сопротивления - не единственное понятие, появившееся в научной литературе за период 2013-2014 гг., предложенное для изучения начальной подвижности зуба в общем случае при отсутствии центра сопротивления [5, 6, 8].

а

б

в

Рис. 1. Виды областей сопротивления: а - эллипс; б - две точки; в - одна точка; черная линия - прямая поступательного воздействия, серая точка -

точка области сопротивления

fP л fa Y ] Г R1

1Ф J ~T IY ß J V ^ j

Определение положения области сопротивления

Настоящая работа является прямым продолжением идей, развитых в статьях [10] и [4], где теория центра и области сопротивления строится на следующей модели. Зуб, окруженный периодонтальной связкой, рассматривается как абсолютно твердое тело, погруженное в линейно-упругую среду, и находится в состоянии устойчивого равновесия. При приложении системы сил тело испытывает малые перемещения в пределах костной лунки. Пусть р - вектор-столбец компонент перемещения полюса;

Ф - вектор-столбец компонент малого поворота зуба вокруг полюса; Я - главный

вектор системы сил; М - главный момент системы сил относительно полюса. Тогда вследствие линейной упругости среды

(1)

где а, в, у - матрицы, которые определяются формой корня зуба, упругими свойствами периодонта и положением полюса. Матрица системы (1) является симметричной и положительно определенной, следовательно, а, в также симметричные и положительно определенные, у - несимметрична в общем случае.

Пусть 8 - симметричная часть матрицы ув; в - антисимметричная часть матрицы ув; В - вектор, соответствующий матрице в. Можно показать, что центр сопротивления существует, если и только если 8 = 0, тогда компоненты вектора в

являются его координатами. Если же 8 Ф 0, то центр сопротивления не существует, но прямые поступательного воздействия могут существовать, и их множество будет определять область сопротивления. Также данную область можно найти аналитически. В этом случае компоненты вектора в являются координатами центра области сопротивления (для случая на рис. 1, а - центр эллипса; рис. 1, б - середина между двумя точками; рис. 1, в - центр сопротивления). Направив оси системы координат по

главным осям матрицы 8, найдем положение области сопротивления, а ее размеры

определятся через собственные значения матрицы 8 .

а б

Рис. 2. Схема задачи: а - пространственная модель зуба, окруженного пружинами; б - положение 7-й пружины в пространстве

Однако для определения положения области сопротивления аналитическим методом необходимо знать коэффициенты матрицы системы (1). Для их вычисления представим упругую среду (периодонтальную связку) в виде набора пружин с заданной жесткостью ci, ориентация в пространстве каждой из которых задается координатами

точки крепления к зубу (xi, yi, zi) и углами ni, ^, 9i с осями системы координат (рис. 2). Далее применим принцип возможных перемещений. Сумма виртуальных работ внешних сил равна сумме виртуальных работ напряжений

£4 = £84,.

j 1

Сообщим точкам виртуальные перемещения и запишем виртуальную работу внешних

сил £ Pi8ri и виртуальную работу напряжений £ ciAli8ri

i i

£ Рг 8?i =£ Ci A?, (2)

i i

где 8ri - виртуальное перемещение точки прикрепления к зубу i-й пружины, Ali - вектор удлинения i-й пружины.

Виртуальное перемещение точки складывается из поступательного виртуального перемещения вместе с полюсом (началом системы координат) 8r00 и поворота на бесконечно малый угол вокруг полюса 8ф

8ri = 8r0 + 8p х r . Аналогично для действительного малого перемещения A ri

Ar =Ar0 + Ap х r ,

где Ap - вектор малого угла поворота вокруг полюса под действием приложенных сил.

Представим 8ri в виде 8ri = 8rxi + 8ryj + 8rzk, a удлинение i-й пружины Ali как

скалярное произведение вектора действительного перемещения Ari и единичного

вектора ni вдоль оси пружины Ali = (A rini) • ni = (A rix cos ni + A riy cos 2, i + A riz cos 9 i) • ni.

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

£Ci Ai Щ = £Ci (Axoa + Ay0b + AZod + AФxe + APyf + APzg)

i i (3)

(a8 Xo+b8 Уо+d8 Zo+е8фх+/8фу + g8 pZ), где a = cos n, b = cos ^, d = cos 9i, e = (zicos ^ - ycos 9i), / = (xt cos 9 i - ztcos n), g = (xicos ^i - yicos ni).

Сумму виртуальных работ внешних сил запишем через обобщенные силы

£ p 8rt = Qx8x0 + Qy8yo + Qz8zo + Mx8px + My8py + Mz8pz. (4)

i

Тогда, сравнивая выражения (3) и (4), обобщенные силы находим как коэффициенты при одинаковых вариациях обобщенных координат. Объединив их в систему уравнений, получим

Ох = X С (АХ0а + АУоЬ + 4 + АФхе + АФу/ + Афг£К

г

Оу = X С (^Хо а + АУоЬ + Аго й + Афх^ + Афу/ + Афг^)),

г

0г = X С (Ахо а + АУоЬ + Аго й + Афх^ + Афу/ + Афг^ )й,

г

< _

мх = X С (Ахоа + АуоЬ + Ай + Афх^ + Афу/ + Афг<?)

г

Му =Х С (Ахоа + АуоЬ + Ай + Афх^ + Афу,/ + Афг<?)/,

г

Мг = X С (Ахоа + АуоЬ + Ай + Афх^ + Афу/ + Афг^))

г

Запишем данную систему уравнений в матричном виде

fа' ' W2 ciab ciad cae caf cagg fAxo ^

б, ciba cb cibd cfie cbf cbg A,o

Q, cida cidb cd2 cde cdf cdg Azo

мх ciea cieb ced ce2 cef ceg АФ*

My cfa cf cfd cfe cß2 cfg Аф ,

lMz J v c,ga cigb ciSd cge cg J

Матрица (5) является симметричной, положительно определенной, имеет размер 6^6, включает 21 независимую компоненту, которые зависят от формы корня зуба, упругих свойств периодонта, положения полюса. Видно, что искомая матрица системы (1) может быть получена как обратная от найденной матрицы (5).

Данный алгоритм запрограммирован в пакете Ма^аЬ и позволяет определять положение центра/области сопротивления в частных случаях.

Связь между видом области сопротивления и формой корня зуба

Для исследования вида области сопротивления направим оси системы координат

по главным осям матрицы 8 , а полюс О* поместим в точку, координаты которой суть компоненты вектора в. Для поиска прямых поступательного воздействия в работе [4] получено уравнение поверхности второго порядка

51 х2 + 5 2 у2 + 5 3 г2 + 515 25 3 = о, (6)

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

Общий случай

Корень зуба задается близким к физиологической форме верхнего резца (не имеющей никаких элементов симметрии) (рис. 3, а). Пользуясь программой, в данном случае можно определить числа 5г: 5Х = -о, об53; 82 = -о, оо7; 83 = о, о344.

б

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

красная линия - х-полуось)

Уравнение поверхности (6) принимает вид канонического уравнения однополостного гиперболоида

х

У

- + -

0,015 0,047 0,021

= 1.

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

Случаи наличия элементов симметрии

а) Корень зуба задается формой, имеющей одну плоскость симметрии (рис. 4). Собственные значения равны 81 =-0,41; б2 = 0; 83 = 0,41. Уравнение поверхности (6) принимает вид канонического уравнения пары пересекающихся плоскостей.

0,41г2 - 0,41х2 = 0.

а

в

z

Рис. 4. Иллюстрация к случаю с одной плоскостью симметрии. Синяя точка - центр сопротивления в плоскости, красные линии - прямые поступательного воздействия

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

б) Корень зуба задается формой прямоугольного параллелепипеда, имеющего две плоскости симметрии (рис. 5). Собственные значения равны 8Х =-0,33; 82 = 0;

83 = 0,33. Уравнение поверхности (6) также принимает вид канонического уравнения пары пересекающихся плоскостей. Стоит отметить, что эти плоскости взаимно перпендикулярны (так как |53| = ).

0,33г2 - 0,33л2 = 0.

Область сопротивления - две точки, лежащие на прямой пересечения двух плоскостей.

У

-с

Ш///М.

111 ¡¡¡р.

4W//li ISP^

ч

x

z 4

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

У

iÜ ш

x

41 ж

щ /Ш\\ш

У

Рис. 6. Иллюстрация к случаю с осью симметрии: а - ось симметрии 4-го порядка; б - ось симметрии 6-го порядка; красные линии - прямые поступательного

воздействия

Z

а

б

в) Корень зуба задается формой тела, имеющей ось симметрии и-порядка (и > 2) (рис. 6). Собственные значения получились равны 81 = 0; 82 = 0; 83 = 0. Следовательно, область сопротивления совпадает с центром сопротивления и все прямые поступательного воздействия пересекаются в одной точке.

Остальные виды областей сопротивления, приведенные в работе [4], здесь не рассматриваются, так как, вероятно, не существует зубов такой формы, для которой бы 8г. приняли соответствующие комбинации знаков. Данное утверждение

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

Связь между видами областей сопротивления

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

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

затем они были объединены на графике (рис. 8) для отражения зависимости от изменяемых параметров. Участок АВ характеризует переход от произвольной формы А к появлению двух плоскостей симметрии В (см. рис. 7, 8). Со стремлением 82 к нулю

область сопротивления в виде эллипса претерпевает следующие изменения в размерах и ориентации: 1) значения полуосей эллипса уменьшаются, причем х-полуось равна 0 в момент В (см. рис. 7, а); 2) происходит совмещение ненулевого собственного вектора

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

симметрии и-порядка (в данном случае и = 4), все 8г. и полуоси равны нулю, означая вырождение области сопротивления в единственную точку - центр сопротивления.

Рис. 7. Иллюстрация к эксперименту: а - экспериментальная начальная модель; б - процесс деформирования формы корня зуба

Z, мм

-0,4

0

5,.

0,4

z, мм

2

-0,4

B

/ C

0

0,4

Полуоси, мм

а б

Рис. 8. Зависимости изменения чисел 8; (а) и полуосей эллипса (б) от изменения формы корня зуба по оси г: а) синяя линия - 81; зеленая линия - 82; красная линия - 83; б) синяя линия - х-полуось; зеленая линия - _у-полуось

4

3

Z, мм

15

-15

-2

Z, мм

15

-15

Рис. 9. Зависимости изменения чисел 8i (а) и полуосей

эллипса (б) от нефизиологичных изменений формы корня зуба по оси г. Обозначения такие же, как на рис. 8

Рис. 10. Изображение стадии деформации тела при г = 11

(т.е. 82 = 0)

0

0

а

Также представляют интерес выход за рамки физиологичности формы корня зуба и наблюдение более широкой картины изменения значений 5,.. Для этого выбрано

тело в форме прямоугольного параллелепипеда. Деформируется оно таким образом, чтобы нарушились все элементы симметрии тела, при этом вытягивается одна из вертикальных граней. Анализируя график зависимости чисел 5,. от изменяемых

параметров (см. рис. 9), можно увидеть, что 52 принимает нулевые значения в трех

точках Zj = -2,4; z2 = 4; z3 = 11. В данных точках вид области сопротивления

соответствует 5j < 0; 52 = 0; 53 > 0. В точке z2 данный результат ожидаем, однако в

точках zx, z3 тело не имеет какой-либо наблюдаемой внешней симметрии, но данный

вид области ему все равно присущ (рис. 10). Уравнение поверхности (6) в данных точках принимает вид пары пересекающихся плоскостей, но они не будут взаимно перпендикулярны, так как угловой коэффициент к = ±^/53/5j Ф1.

Заключение

Данная работа является продолжением идей и теории, развитой в статьях [4, 10], в которых показано, что центр сопротивления зуба существует не всегда, и для тех случаев, когда его нет, введено новое понятие области сопротивления. Результатами данной работы являются создание методики для определения положения области сопротивления аналитически или нахождением всего множества прямых поступательного воздействия, установление связи между видом области сопротивления и геометрическими параметрами системы «зуб-периодонт», рассмотрен процесс перехода из одного вида области в другой. Также экспериментально показано, что центр сопротивления может существовать у тел с осью симметрии n-порядка (n > 2).

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

Автор выражает благодарность Ю.И. Няшину, М.А. Осипенко, В. С. Туктамышеву за помощь, ценные замечания и постоянный интерес к работе.

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта №15-01-04932

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

1. Дубинин А. Л., Няшин Ю.И., Осипенко М.А. Анализ развития понятия «центр сопротивления зуба» // Российский журнал биомеханики. - 2014. - Т. 18, № 4. - С. 452-470.

2. Дубинин А.Л., Няшин Ю.И., Осипенко М.А., Туктамышев В.С. Ортодонтическое перемещение зубов. Понятие центра сопротивления и области сопротивления зуба // Тезисы докладов XI Всероссийской конференции с международным участием и школы-семинара для молодых ученых «Биомеханика-2014», 1-4 декабря 2014. - Пермь, 2014. - C. 53-54.

3. Дубинин А.Л., Осипенко М.А., Няшин Ю.И., Туктамышев В.С. Понятие центра сопротивления и области сопротивления зуба // Математическое моделирование в естественных средах. - Пермь, 2014. - C. 91-92.

4. Осипенко М.А., Няшин Ю.И., Няшин М.Ю., Дубинин А.Л. Область сопротивления зуба: определения и свойства // Российский журнал биомеханики. - 2013. - Т. 17, № 2. - С. 31-38.

5. Профит У.Р. Современная ортодонтия. - М.: МЕДпресс-информ, 2008. - 560 с.

6. Burstone C.J. Biomechanics of Tooth Movement. - Philadelphia: Lea & Febiger, 1962. - 213 p.

7. Dathe H., Nägerl H., Kubein-Meesenburg D. A caveat concerning center of resistance // Journal of Dental Biomechanics. - 2013. - Vol. 4. - P. 1-7.

8. Fish G.D. Some engineering principles of possible interest to orthodontists // Dental Cosmetics. - 1917. -Vol. 59. - P. 881-889.

9. Geiger M.E., Schmidt F., Lapatki B.G. Centre of resistance of constructed tooth models of various morphologies // 18th Symposium on Computational Biomechanics in Ulm, 13-14 May 2013. - Ulm, 2013.

10. Osipenko M.A., Nyashin M.Y., Nyashin Y.I. Center of resistance and center of rotation of a tooth: the definitions, conditions of existence, properties // Russian Journal of Biomechanics. - 1999. - Vol. 3, № 1. -P. 1-11.

11. Ravinda N. Biomechanics and Estetic Strategies in Clinical Orthodontics. - St. Louis: Elsevier, 2009. - 396 p.

12. Viecilli R.F., Budiman A., Burstone C.J. Axes of resistance for tooth movement: does the center of resistance exist in 3-dimensional space? American Journal of Orthodontics and Dentofacial Orthopaedics. -2013. - Vol. 143, № 2. - P. 163-172.

13. Vollmer D., Bourauel C., Maier K., Jäger A. Determination of the centre of resistance in an upper human canine and idealized tooth model // European Journal of Orthodontics. - 1999. - Vol. 21, № 6. -P. 633-648.

REGION OF RESISTANCE OF TOOTH: EXPERIMENTAL DETERMINATION

A.L. Dubinin (Perm, Russia)

In order to study the initial tooth movement, the concept of the center of resistance was introduced. Properties of this point are: if the applied force system is reduced to the single force with the line of action passing through the center of resistance, then the tooth will start to move translationally; if the applied force system is reduced to the force couple, then the tooth will rotate around this center. Later, it was shown that the center of resistance exists at the exceptional conditions (existing of the axis of symmetry for the root of the tooth). Therefore, in cases where there is no center of resistance, in work a new concept "region of resistance of the tooth" was introduced. The region is called the region of resistance of a tooth with the minimal diameter while possessing the following properties: any line of translational action passes through this region and through any point of this region passes some line of translational action. It has been shown that region of resistance can be an ellipse, two points, and one point. In presented study, the method of finding the region of resistance for the particular case was developed by two approaches: 1) finding the set of all lines of translational action; 2) analytically. The relation between the region of resistance and tooth root shape was established experimentally. Also, the transition from one type of region of resistance to other was shown by varying geometrical parameters of the body. Some properties of certain types of region investigated. Another case of the existence of the center of resistance was shown, which has not been considered earlier.

Key words: center of resistance of tooth, region of resistance of tooth, dentition, orthodontic tooth movement, experimental determination.

Получено 20 февраля 2014