USING OF DIGITAL TECHNOLOGIES AT BIOMECHANICAL MODELLING OF ORTHODONTIC TOOTH MOVEMENT
- Authors: Dubinin A.L1, Osipenko M.A1, Nyashin Y.I1
- Affiliations:
- Issue: Vol 23, No 3 (2019)
- Pages: 359-374
- Section: Articles
- URL: https://ered.pstu.ru/index.php/rjb/article/view/2349
- DOI: https://doi.org/10.15593/RZhBiomeh/2019.3.03
- Cite item
Abstract
Dental anomalies (malocclusion, crooked teeth, etc.) are a common problem. Due to them, violations of the aesthetics of the face, blood supply to the brain, respiration, posture, facial expressions, pronunciation of words, digestion and etc. are possible. Correction of dental anomalies consists in moving the teeth to a new position corresponding to the physiological norm under the influence of the orthodontic load (bracket system, elastopositioner, etc.). The main question in orthodontic treatment is the determination of such load (value, direction, place of application). This work continues a series of publications by authors dedicated to solving this question. Earlier, the concept of “center of resistance of a tooth” was introduced and determined. Using this consept, it is possible to study the initial movement of the tooth within the framework of the bone hole. In accordance with the introduced conditions for the existence of this point, it turned out that the "center of resistance of a tooth" does not exist always. Therefore, the authors proposed a new concept of “region of resistance of a tooth” determined by sets of translational action lines and a pair of rotation axes. The purpose of this work is to study them, namely, the conditions for the existence of these lines, the relative position and, on the basis of this, to demonstrate the classification of the types of tooth resistance region.
Full Text
Введение В ранних работах авторов развит биомеханический подход к исследованию начальной стадии перемещения зубов (в рамках костной лунки, т.е. до начала процессов перестройки костной ткани) при ортодонтическом лечении [4, 6, 11, 12]. Изменение положения зуба происходит в результате действия на него ортодонтической нагрузки, развиваемой корректирующей конструкцией (брекет-система, эластопозиционер и пр.). Для исследования этого малого перемещения была построена биомеханическая модель, в основу которой легло понятие «центр сопротивления зуба». Впервые этот термин был введен Д. Фишем в 1917 г. в работе [10]. Рассматривается элемент зубочелюстной системы, включающий в себя сам зуб и окружающую его связку - периодонт, которая заполняет остальное пространство костной лунки. Зуб рассматривается как абсолютно твердое тело, периодонт принимается линейно-упругой средой. Окружающая кость учитывается в виде граничных условий на перемещении внешней границы периодонта. Нагрузка прикладывается к коронке зуба. Считается, что ортодонтические силы малы (около 100 гс) и действуют продолжительное время [1-3, 5, 7-9, 13]. В отсутствие нагрузки считается, что зуб находится в состоянии устойчивого равновесия [12]. Вследствие линейной упругости среды связь действующей системы сил и получаемых перемещений описывается следующей зависимостью: (1) где - обобщенные координаты зуба (перемещение полюса и вектор малого поворота зуба соответственно); матрицы характеризуют податливость среды (зависят от геометрической формы корня зуба, упругих свойств периодонта, положения полюса); - главный вектор системы сил; - главный момент системы сил относительно полюса. Перемещение точки зуба с радиусом-вектором находится по формуле (2) Центр сопротивления зуба и его окружения определен как точка, удовлетворяющая двум требованиям: 1) если то для любого (если приложенные силы сводятся к паре сил с моментом, равным главному моменту системы сил, то зуб повернется вокруг центра сопротивления), 2) если то для любого (если приложенные силы сводятся к равнодействующей с линией действия, проходящей через данный центр, то зуб будет двигаться поступательно (необязательно в направлении действия силы). Здесь - радиус-вектор центра сопротивления. Введены условия существования центра сопротивления: «Если матрица антисимметрична, то центр сопротивления существует, единственен и его положение определяется формулой . В противном случае центр сопротивления не существует». Установлено, что центр сопротивления существует при наличии у системы «зуб - периодонт» некоторых идеализированных элементов симметрии: плоскости симметрии, оси симметрии. В случае моделирования данной системы формой, близкой к реальной форме зуба (не имеющей элементов симметрии), данное условие нарушается и центра сопротивления быть не может. В этом случае в работах авторов [4, 6] предложено и исследовано новое понятие «область сопротивления зуба». Оно является естественным обобщением понятия «центр сопротивления» и сохраняет ряд его свойств. Для удобства дальнейшего повествования введены новые термины: 1. Прямая поступательного воздействия - линия действия силы, приводящей зуб к поступательному движению (рис. 1). 2. Ось поворота парой - прямая, вокруг которой зуб поворачивается при нагружении парой сил (рис. 2). Пользуясь новыми терминами, можно переформулировать определение известного понятия центра сопротивления. Рис. 2. Изображение оси поворота парой Рис. 1. Изображение прямой поступательного воздействия а б Рис. 3. Прямые поступательного воздействия (синие линии) и оси поворота парой (красные линии) пересекаются в едином центре (а), не пересекаются (б) «Центр сопротивления зуба - точка, имеющая следующие свойства: 1) любая прямая поступательного воздействия проходит через этот центр, и любая прямая, проходящая через этот центр, является прямой поступательного воздействия (рис. 3); 2) любая ось поворота парой проходит через этот центр, и любая прямая, проходящая через эту точку, является осью поворота парой (см. рис. 3, а)». Далее введем обозначения: - симметричная часть матрицы ; - антисимметричная часть матрицы ; - вектор, соответствующий матрице ( для любого вектора Так как условием существования центра сопротивления является антисимметричность матрицы то справедливо также следующее утверждение: центр сопротивления существует тогда и только тогда, когда , и его координаты могут быть найдены как компоненты вектора . Если не антисимметрична (т.е. ), то центра сопротивления не существует. В данном случае прямые поступательного воздействия и оси поворота парой уже не будут пересекаться в одной точке, а будут расположены каким-то другим образом (см. рис. 3, б). Цель данной работы заключается в том, чтобы исследовать наборы прямых поступательного воздействия и осей поворота парой, а именно условия их существования, взаимное расположение, и на основе этого продемонстрировать классификацию видов области сопротивления зуба. Исследование множества прямых поступательного воздействия Пусть - вектор приложенной к зубу силы, а - радиус-вектор точки приложения этой силы. Рассматривая второе уравнение из системы (1), запишем его в виде . (3) Вынося общий множитель за скобки, получим Учитывая, что - симметричная матрица, получим . Тогда . Учитывая определение матриц , - симметричная часть матрицы ; - антисимметричная часть матрицы ), а также то, что при транспонировании антисимметричной матрицы у нее меняется знак, данное выражение можно привести к виду . Введем вектор , соответствующий матрице ( для любого вектора , получим . В конечном виде данное выражение выглядит следующим образом: . Исследуем множество прямых поступательного воздействия. Для дальнейшего удобства перейдем в специальную систему координат с началом в точке, расположенной на конце вектора , и выберем за полюс эту точку. Оси этой системы координат направим вдоль главных осей матрицы . Тогда условие можно представить в виде . Учитывая, что матрица невырождена, приводя симметричную часть матрицы к диагональному виду и представляя как , получим . (4) Представим это выражение в виде системы линейных однородных уравнений: (5) Далее, при известных значениях матрицы в каждой точке поверхности приложим такую силу , чтобы она удовлетворяла условию (4). Ее компоненты могут быть найдены из решения системы (5). Линия действия вектора полученной силы будет являться прямой поступательного воздействия. Таким образом определим множество прямых поступательного воздействия. При этом перемещение полюса находится по формуле Так как приравниваем к нулю определитель матрицы (4) и получаем уравнение поверхности второго порядка (6) Запишем уравнение прямолинейных образующих поверхности второго порядка (в зависимости от знаков «+», «-» определяется первое или второе семейство соответственно), проходящих через точку (7) Значения чисел определяют вид поверхности. Таким образом, приходим к классификации множеств прямых поступательного воздействия. Множество прямых поступательного воздействия 1-го типа Все числа не равны нулю и имеют одинаковый знак: Выражение (6) в данном случае выглядит следующим образом: Приведем данное выражение к каноническому виду: В этом случае поверхность (6) является пустым множеством. Прямых поступательного воздействия не существует. При значениях выражение (6) аналогично данному. Поэтому этот набор значений не дает нового вида области сопротивления. Множество прямых поступательного воздействия 2-го типа Все числа не равны нулю и имеют различные знаки: Выражение (6) в данном случае выглядит следующим образом: Приведем данное выражение к каноническому виду: В этом случае поверхность (6) - это однополостный гиперболоид (рис. 4). В то же время из (5) определим множество прямых поступательного воздействия: Прямая поступательного воздействия Рис. 4. Множество прямых поступательного воздействия 2-го типа Линия действия вектора совпадает с прямой (7) (со знаком «+»). Следовательно, множество прямых поступательного воздействия совпадает с одним из семейств прямолинейных образующих. Набор знаков также не дает нового вида области сопротивления ввиду того, что выражение (6) не дает нового вида при смене знаков . Множество прямых поступательного воздействия 3-го типа Одно из чисел равно нулю, остальные не равны нулю и имеют одинаковый знак: Выражение (6) в данном случае выглядит следующим образом: Приведем данное выражение к каноническому виду: В этом случае поверхность (6) - мнимые пересекающиеся плоскости. Прямых поступательного воздействия не существует. Множество прямых поступательного воздействия 4-го типа Одно из чисел равно нулю, остальные отличны от нуля и имеют различные знаки: Выражение (6) в данном случае выглядит следующим образом: Приведем его к каноническому виду: В этом случае поверхность (6) - это две пересекающиеся плоскости: Определим множество прямых поступательного воздействия из (5): Если взять первую плоскость то из (7) получим: где . Эта прямая проходит через точку при любых (), принадлежащих первой плоскости. Таким образом, получаем семейство прямых в плоскости проходящих через точку (рис. 5). Если взять вторую плоскость то прямые поступательного воздействия также проходят через точку в этой плоскости (см. рис. 5). Таким образом, прямые поступательного воздействия - это два плоских пучка прямых, которые проходят через две точки и , лежащие на линии пересечения данных плоскостей (см. рис. 5). Множество прямых поступательного воздействия 5-го типа Лишь одно из чисел не равно нулю: Выражение (6) в данном случае выглядит следующим образом: При этом поверхность (6) вырождается в пару совпадающих плоскостей. С2 С1 Прямая поступательного воздействия С1 Рис. 5. Множество прямых поступательного воздействия 4-го типа Аналогично прямые поступательного воздействия собираются в плоский пучок прямых, который проходит через одну точку С3 (рис. 6). Множество прямых поступательного воздействия 6-го типа Все числа равны нулю: В данном случае (6) - все пространство. Из (5) можно показать, что Видно, что при вектор - любой. Это и означает семейство всех прямых, проходящих через точку . Таким образом, прямые поступательного воздействия - пространственный пучок прямых, проходящих через точку С4 (рис. 7). В данном случае что удовлетворяет условию существования центра сопротивления, который совпадает с точкой С4. Рис. 6. Множество прямых поступательного воздействия 5-го типа С4 Прямая поступательного воздействия Рис. 7. Множество прямых поступательного воздействия 6-го типа Исследование множества осей поворота парой Далее пусть - момент пары сил, приложенной к зубу. Тогда система (1) примет следующий вид: Выразим и подставим в (2). Тогда Учитывая определение матриц , ( - симметричная часть матрицы - антисимметричная часть матрицы ), получим Вектор представим в виде антисимметричной матрицы ( для любого вектора , найдем (8) Перейдем к уже введенной ранее системе координат с началом, расположенным в точке, координаты которой являются компонентами вектора . Из (8) и условия (при чистом повороте центр сопротивления не перемещается) получаем уравнение оси поворота (9) Представим (7) в виде системы линейных однородных уравнений (10) Далее, при известных значениях матрицы через каждую точку поверхности проведем ось вращения (линия действия вектора малого поворота удовлетворяющую условию (9). Его компоненты могут быть найдены из решения системы (10). Таким образом определим множество осей поворота парой. Так как приравниваем к нулю определитель матрицы (9) и получаем также уравнение поверхности 2-го порядка Таким образом, приходим к классификации типов области сопротивления через множество осей поворота парой. Далее рассмотрим отдельно каждый вид множеств осей поворота парой аналогично тому, как это было проведено при исследовании посредством множества прямых поступательного воздействия. Множество осей поворота парой 1-го типа Все числа не равны нулю, имеют одинаковый знак: Выражение (11) в этом случае выглядит следующим образом: (11) Приведем данное выражение к каноническому виду: При этом соотношение (11) не имеет решения и осей поворота парой нет. При значениях выражение (11) аналогично данному. Поэтому этот набор значений не дает нового вида области сопротивления. Множество осей поворота парой 2-го типа Все числа не равны нулю, имеют разные знаки: Выражение (11) в данном случае выглядит следующим образом: Приведем данное выражение к каноническому виду: В этом случае поверхность (11) - это однополостный гиперболоид (рис. 8). Из (10) определим множество осей поворота: В данном случае линия действия вектора совпадает со второй прямой (6) (со знаком «-»). Следовательно, множество осей поворота также совпадает со вторым из семейств прямолинейных образующих (рис. 8). Набор знаков также не дает нового вида области сопротивления, так как (11) не дает нового вида при смене знаков . Рис. 8. Множество осей поворота парой 2-го типа Множество осей поворота парой 3-го типа Одно из чисел равно нулю, остальные отличны от нуля и имеют одинаковый знак: Выражение (11) в данном случае выглядит следующим образом: Приведем данное выражение к каноническому виду: В этом случае поверхность (11) - мнимые пересекающиеся плоскости. Осей поворота парой не существует. Множество осей поворота парой 4-го типа Одно из чисел равно нулю, остальные отличны от нуля и имеют различные знаки: Выражение (11) в данном случае выглядит следующим образом: Приведем данное выражение к каноническому виду: Аналогично ранее рассмотренному случаю из классификации типов области сопротивления через прямые поступательного воздействия (также четвертого типа при ) получаем подобный результат: Однако в данном случае параметр . Таким образом, если взять первую плоскость тогда оси поворота парой образуют семейство прямых, которое проходит через точку , введенную ранее в этой плоскости (см. рис. 9). Если взять вторую плоскость то оси поворота парой образуют семейство прямых, которое проходит через точку в этой плоскости. Выражение (11) определяет два семейства осей поворота, которые параллельны двум плоскостям. Соответствующий набор осей поворота парой предстает в виде двух плоских пучков прямых (см. рис. 9). Они лежат в тех самых плоскостях, образованных прямыми поступательного воздействия. В то же время эти пучки прямых проходят также через те же две точки, из которых состояла область сопротивления (см. рис. 5). Однако наборы осей поворота парой и прямых поступательного воздействия, которые проходят через одну из точек области сопротивления, оказываются лежащими в разных плоскостях. Множество осей поворота парой 5-го типа Только одно из чисел отлично от нуля: С2 С1 Ось поворота парой Рис. 9. Множество осей поворота парой 4-го типа Выражение (11) в данном случае выглядит следующим образом: Выражение (11) определяет семейство направлений, параллельное некоторой плоскости. Соответствующий набор осей поворота парой предстает в виде плоского пучка прямых. Он лежит в той же плоскости, которая образовалась посредством прямых поступательного воздействия. Данный пучок прямых проходит через ту же точку, которая представляла собой область сопротивления (рис. 10). Ось поворота парой С3 Рис. 10. Множество осей поворота парой 5-го типа Множество осей поворота парой 6-го типа Все числа равны нулю: В этом случае соответствующее множество осей поворота парой (11) есть все пространство. Аналогично подобному случаю из классификации множеств прямых поступательного воздействия при соответствующих условиях видно, что оси поворота парой являются пространственным пучком прямых, который проходит через одну точку С4 (рис. 11). В данном случае и существует центр сопротивления, совпадающий с точкой С4. Рис. 11. Множество осей поворота парой 6-го типа С4 Ось поворота парой Выводы Таким образом, продемонстрирована классификация, приведены методика определения вида и положения области сопротивления зуба в зависимости от геометрических и механических параметров системы «зуб - периодонт». Пользуясь этим понятием, можно определять точку приложения и направление системы ортодонтических сил. Такие данные способны объективизировать эмпирические методы ортодонтии: они позволяют теоретически обосновать выбор оптимального варианта приложения ортодонтической нагрузки с количественной оценкой ее величины с учетом формы, положения, топографии зуба, сформулировать рекомендации по настройке и установке ортодонтического аппарата.About the authors
A. L Dubinin
M. A Osipenko
Yu. I Nyashin
References
- Астахов Н.А., Гофунг Е.М., Катц А.Я. Ортопедическая стоматология. - Л.: Медгиз, 1940.
- Величина сил ортодонтических аппаратов и морфологические изменения при их применении [Электронный ресурс]. - URL: www.mysurgeon.ru/baza_znaniy/stomatologiya_657 (дата обращения: 10.01.2016).
- Губская А.Н., Золотухин И.В., Корешков В.А. Дозирование ортодонтической нагрузки при перемещении зубов. - Киев, 1982. - 41 с.
- Дубинин А.Л. Область сопротивления зуба: экспериментальное определение // Российский журнал биомеханики. - 2015. - Т. 19, № 1. - С. 79-89. doi: 10.15593/RJBiomech/2015.1.07
- Калвелис Д.А. Ортодонтия. Зубочелюстные аномалии в клинике и эксперименте. - Л.: Медицина, 1964. - 238 с.
- Осипенко М.А., Няшин Ю.И., Няшин М.Ю., Дубинин А.Л. Область сопротивления зуба: определения и свойства // Российский журнал биомеханики. - 2013. - Т. 17, № 2. - С. 31-38.
- Персин Л.С. Ортодонтия. Диагностика и лечение зубочелюстных аномалий: руководство для врачей. - М.: Медицина, 2004. - 360 с.
- Проффит У.Р. Современная ортодонтия. - М.: МЕДпресс-информ, 2008. - 560 с.
- Хорошилкина Ф.Я. Руководство по ортодонтии. - М.: Медицина, 1999. - 797 с.
- Fish G.D. Some engineering principles of possible interest to orthodontists // Dental Cosmetics. - 1917. - Vol. 59. - P. 881-889.
- Nyashin Y.I., Nyashin M.Y., Osipenko M.A., Lokhov V.A., Dubinin A.L., Rammerstorfer F., Zhurov A. Centre of resistance and centre of rotation of a tooth: experimental determination, computer simulation and the effect of tissue nonlinearity // Computer Methods in Biomechanics and Biomedical Engineering. - 2016. - Vol. 3. - P. 229-239. doi: 10.1080/10255842.2015.1007961
- 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.
- Schwartz A.M. Tissue changes incidental to orthodontic tooth movement // International Journal of Orthodontics. - 1932. - Vol. 18. - P. 331-352.
Statistics
Views
Abstract - 81
PDF (Russian) - 39
Refbacks
- There are currently no refbacks.