THE POISSON'S RATIO OF DENTIN AS ANISOTROPIC MEDIum WITH HEXAGONAL SYMMETRY
- Authors: Muslov S.A1, Lisovenko D.S1, Arutyunov A.S1, Pivovarov A.A1, Manin A.I1, Kirakosyan L.G1, Kharah Y.N1, Arutyunov S.D1
- Affiliations:
- Issue: Vol 22, No 4 (2018)
- Pages: 537-545
- Section: Articles
- URL: https://ered.pstu.ru/index.php/rjb/article/view/2391
- DOI: https://doi.org/10.15593/RZhBiomeh/2018.4.09
- Cite item
Abstract
Full Text
Введение Упругие свойства твердых тканей зуба и гидроксиапатита как минеральной их составляющей достаточно часто вызывают интерес у исследователей [4-7, 12-17]. Однако остаются актуальными и до конца не изученными некоторые аспекты этой тематики, в частности величина коэффициента Пуассона поперечной деформации эмали и дентина как анизотропных сред с гексагональной симметрией. Этот вопрос особенно интересен с позиций прочностных свойств тканей зуба, представляющих в контексте теории разрушения микронеоднородную среду со сложным характером напряженно-деформированного состояния в микрообъемах. Анизотропия и неоднородность твердых тканей зуба сред повышают степень концентрации микронапряжений, ответственных за возникновение микротрещин Гриффитса и пластических деформаций, в зависимости от степени анизотропии и неоднородности. В данном сообщении авторы представляют достаточно подробное исследование величины коэффициента Пуассона дентина как анизотропной среды гексагональной сингонии - три элементарных ячейки образуют правильную призму на шестигранном основании. В связи с тем что в литературе до сих пор не сложились единые обозначения и расчетные формулы для определения упругих характеристик анизотропных сред, данное исследование проведено на основании двух различных методик расчета [2, 16]. Расчеты и результаты 1. Коэффициент Пуассона (коэффициент поперечной деформации) µ является весьма информативной характеристикой при анализе деформационного поведения твердых тел. Для реально упругих тел, увеличивающих свой объем при растяжении и уменьшающих при сжатии, величина µ может лежать только в пределах от 0 до 0,5 [1]. Практически для большинства материалов коэффициент Пуассона находится еще в более узком интервале: от 0,2 до 0,4 [3]. Однако, по современным представлениям, диапазон возможных значений коэффициента Пуассона существенно расширен за оба предела: может быть отрицательным µ < 0 (ауксетики) и больше 0,5. Известно, что коэффициент Пуассона характеризует отношение поперечной деформации к продольной и является одной из упругих констант материала. Для дентина коэффициент, по литературным данным, лежит в интервале 0,29-0,33 [9]. По результатам ультразвуковых измерений [15] он равен 0,32. 2. Считается, что дентин представляет собой биокомпозит и состоит приблизительно на 45-70% из неорганического материала в форме кристаллов гидроксиапатита [10, 14]. Кристаллы гидроксиапатита расположены между коллагеновыми волокнами и по классу симметрии относятся к гексагональной сингонии. Поэтому есть все основания рассматривать дентин как анизотропную среду, а это значит, что его упругие свойства описываются матрицей упругих постоянных или коэффициентов податливости Коэффициенты Пуассона такой среды в генерализованной общей форме могут быть определены как где - компоненты тензора деформации. В результате гексагональные системы и дентин тоже описываются не одним, а двумя коэффициентами Пуассона и [16]. Штрихи показывают, что данные компоненты µ отнесены не к кристаллографической, а к специальной (лабораторной) системе координат, связанной с зубом с помощью углов Эйлера и отличной, вообще говоря, от кристаллографической. По литературным данным, значения упругих постоянных дентина достаточно хорошо известны [12, 13], поэтому, зная их, можно вычислить коэффициенты податливости и Пуассона дентина по формулам [8, 16] для гексагональных кристаллов: и где и - углы Эйлера (описывающие поворот абсолютно твердого тела в трехмерном евклидовом пространстве). Значения упругих постоянных , ГПа: с11 = 37; с12 = 16,6; с13 = 8,7; с33 = 39; с44 = 5,7. Значения коэффициентов податливостей дентина , ГПа-1: s11 = 0,0346; s12 = -0,0145; s13 = -0,0045; s33 = 0,0276; s44 = 0,1754. Результаты вычисления коэффициентов Пуассона дентина как кристаллической системы с гексагональной структурой представлены на рис. 2 (в полярных координатах (а) и в декартовых координатах (б). Вследствие симметрии и для простоты изложения графики приведены только для Из рис. 2 отчетливо следует выраженное анизотропное поведение коэффициента Пуассона кристаллов дентина. Его значения в зависимости от направления в пространстве изменяются в весьма широких пределах (в 4 с лишним раза). Рис. 1. Гексагональная ячейка и ортогональная система координат Xi, где ОX3 совпадает с направлением оси кристалла. Вращение ортогональной системы Xi' относительно системы Xi описывается в терминах трех углов Эйлера θ, φ и ψ a Min µ31(θ) = 0,132 Max µ31(θ) = 0,534 б Рис. 2. Полярные координаты (а) и декартовы координаты (б). Коэффициент Пуассона дентина для различных направлений и , а-f: радиан соответственно Отметим аномально высокую величину максимального значения коэффициента Пуассона (0,534) вдоль ряда направлений, что крайне необычно для материалов. Это означает, что при сжатии локальных областей дентина вдоль этих направлений их объем будет увеличиваться, а при растяжении - наоборот, уменьшаться. Упругое поведение под нагрузкой реставрационного материала принципиально другое (μ пломбировочного материала существенно меньше 0,5), поэтому при сжатии его объем уменьшается, а при растяжении - увеличивается. Это несоответствие в деформационном поведении на границе разнородных сред пломбировочного материала и дентина может привести к образованию доменов перенапряжений на данной границе, ослаблению сцепления реставрационного материала с дентином и, как негативный итог, деградации фиксации и часто встречаемой ситуации разрушения прямой, а порой и непрямой реставрации, особенно композиционными материалами. Кроме того, известно, что величина коэффициента Пуассона тесно связана и прямо влияет на деформационную прочность кристаллов, а именно на трещиностойкость. Как известно, очагами разрушения материалов служат небольшие зародышевые микротрещины различного происхождения. Считается, что микротрещины в момент зарождения имеют длину ~ 10-6-10-4 мм. Распространение трещин может привести к разрушению зуба, т.е. разделению его на части. Разрушение носит хрупкий характер и происходит, как правило, при малых деформациях, без заметных пластических. В теории Гриффитса для объемного напряженного состояния показано, что величина критического напряжения σ, начиная с которого хрупкая трещина стартует и неограниченно растет, зависит от значения µ [9]. Это следует из нижеприведенной формулы и графика на рис. 3: где γs - удельная поверхностная энергия дентина; l - длина микротрещины; Е - модуль Юнга; E* - приведенный модуль Юнга, E*= E/1 - µ2. Из рис. 3 видно, что чем выше коэффициент Пуассона, тем сильнее критическое напряжение зависит от него. С учетом того, что максимальные напряжения действуют в периферийной части ткани зуба, вероятность старта трещины там также максимальна. В средней части зуба, где располагается пульпа, внутренние силы минимальны, соответственно, вероятность старта будет меньше. Мы можем сделать некоторые оценки. Если положить и l = 1 мкм, то получим при µ = 0,13 σ = 2,64 ГПа, а при µ = 0,52 σ = 3,12 ГПа. Далее, если гипотетически предположить, что данные напряжения приложены к площади 1 мм2, то сила, вызвавшая это напряжение, будет равна в первом случае 2,64 кН, во втором - 3,12 кН. По порядку величины это почти силы, действующие в зубочелюстной системе человека. Напомним, что жевательная мышца, по данным [17], развивает усилия до 1 кН (например, при раскалывании ореха, косточек слив или абрикосов). 3. Заслуживает внимания подход, основанный на расчете экстремальных значений коэффициента Пуассона кристаллических сред и построении изоповерхностей коэффициента в пространстве. Результаты вычисления коэффициентов Пуассона дентина как кристаллической гексагональной структуры: (рис. 5). Последние четыре цифры в квадратных скобках указывают направление растяжения, а первые четыре - направление поперечной деформации. Четырехиндексовая система Миллера-Браве используется для определения направлений в гексагональной решетке (рис. 4). Рис. 3. Зависимость напряжения старта трещины, ГПа, от коэффициента Пуассона для объемного напряженного состояния дентина, [5] и l = 10-6 м = 1 мкм Рис. 4. Примеры кристаллографических направлений в гексагональной решетке (базисная плоскость, ось z перпендикулярна плоскости чертежа) В гексагональной решетке начало координат помещают в центр основания элементарной ячейки. Кристаллографические оси х и у проходят из этого центра через вершины шестиугольного основания элементарной ячейки, располагаясь под углом 120° одна к другой, а ось z является вертикальной осью гексагональной призмы. За единицу измерения вдоль осей х и у принимают период решетки, а, а вдоль оси z - период, с. В плоскости базиса проводят дополнительную ось u, расположенную под углом 120° к осям х и у. Направление -u находится между направлениями +x и +y. Дополнительный индекс i определяют точно так же, как и индексы Миллера, и ставят на третьем месте [hkil]. a б Рис. 5. Изменчивость коэффициента Пуассона µ дентина в пространстве На рис. 5. представлена поверхность коэффициента Пуассона для дентина. Он подтверждает выраженный анизотропный характер коэффициента Пуассона. Его значения в зависимости от направления (углов Эйлера θ и ψ) изменяются в широких пределах (до четырех с лишним раз). Отметим также аномально высокую величину максимального значения коэффициента Пуассона (0,54) вдоль ряда направлений, что оказывается выше значения верхнего предела для коэффициента Пуассона изотропных материалов. Кроме того, максимальное и минимальное значения коэффициента Пуассона хорошо согласуются с данными, полученными с помощью первого подхода для расчета коэффициента Пуассона как анизотропной среды (см. рис. 1, 2, a). В заключение отметим, что аналогичный анализ может быть выполнен и для эмали зубов (см. рис. 5, б), состоящей из неорганических минералов, в основном кристаллов апатитов с гексагональной структурой, на 90-97% [10, 14]. В частности, минеральная основа эмали представлена гексагональными кристаллами гидрокси-, карбонат-, хлор-, фторапатитов, а эмалевые призмы являются основным структурным образованием. • Значения упругих постоянных для эмали, ГПа: с11 = 115; с12 = 42,4; с13 = 30; с33 = 125; с44 = 22,8. • Значения коэффициентов податливостей для эмали, ГПа -1: s11 = 0,0104; s12 = -0,0034; s13 = -0,0017; s33 = 0,0088; s44 = 0,0439, • Экстремальные значения коэффициента Пуассона для эмали, а также значения коэффициентов Пуассона в частных ориентациях: . Выводы 1. Дентин зубов не является изотропной средой вследствие симметрии его минеральной составляющей - кристаллов гидроксиапатита. 2. Установлена выраженная анизотропия коэффициента Пуассона дентина на основе расчетов по формулам упругих постоянных и коэффициентов податливости для гексагональной сингонии. Максимальное значение коэффициента Пуассона (0,53-0,54) более чем в 4 раза (4,15) превышает минимальное (0,13). 3. Максимальное значение коэффициента Пуассона оказывается выше значения верхнего предела для коэффициента Пуассона изотропных, в том числе реставрационных, материалов, используемых в стоматологии, что может локально сказываться на качестве реставраций. В данном контексте высказано предположение, что установленная упругая анизотропия модели дентина с кристаллической гексагональной симметрией является клинически нежелательным фактором. 4. Более тщательный анализ упругой анизотропии дентина как минерально-органического комплекса и микронеоднородной гетерофазной системы возможен с привлечением достижений теории анизотропных сред с кристаллографической текстурой (несмотря на то что все призмы минеральных веществ в дентине обладают одинаковым или схожим кристаллическим строением, они отличаются взаимной ориентацией кристаллографических осей). Дальнейший анализ должен быть основан на знании характера пространственного распределения кристаллографических осей отдельных минеральных призм дентина. 5. Изучение показателей анизотропии дентина как анизотропной неоднородной среды приобретает практическое значение при исследовании проблем прочности тканей зуба и качества реставраций. Значения коэффициента Пуассона эмали принципиально не отличаются от таковых у дентина. Они по-прежнему весьма велики и свидетельствуют об анизотропии упругих свойств эмали зуба (хотя и в меньшей степени:About the authors
S. A Muslov
D. S Lisovenko
A. S Arutyunov
A. A Pivovarov
A. I Manin
L. G Kirakosyan
Y. N Kharah
S. D Arutyunov
References
- Беломестных В.Н., Теслева Е.П. Коэффициент Пуассона и параметр Грюнайзена твердых тел // Известия Томского политехнического университета. - 2003. - Т. 306, № 5. - С. 8-12.
- Гольдштейн Р.В., Городцов В.А., Лисовенко Д.С. Изменчивость упругих свойств гексагональных ауксетиков // Доклады Академии наук. - 2011. - Т. 441, № 4. - С. 468-471.
- Иванов Г.П., Лебедев Т.А. О физическом смысле коэффициента Пуассона // Труды Ленингр. политехн. ин-та им. М.И. Калинина. - 1964. - № 236. - С. 38-46.
- Лебеденко И.Ю., Арутюнов С.Д., Муслов С.А., Усеинов А.С. Нанотвердость и модуль Юнга зубной эмали // Вестник Российского университета дружбы народов. Серия: Медицина. - 2009. - № 4. - С. 637-638.
- Лебеденко И.Ю., Арутюнов С.Д., Муслов С.А., Усеинов А.С. Исследование наномеханических свойств зубной эмали // Кафедра. - № 32. - С. 24-28.
- Лебеденко И.Ю., Арутюнов С.Д., Муслов С.А., Болатаев З.Б. Упругая анизотропия твердых тканей зуба // Медицинская физика - 2010: сб. материалов III Евразийского конгресса по медицинской физике и инженерии. - 2010. - М.: Изд-во МГУ, 2010. - Т. 1. Медицинская биофизика. - С. 279-280.
- Лебеденко И.Ю., Арутюнов С.Д., Усеинов А.С., Муслов С.А., Брандт Н.Н., Апресян С.В. Исследование воздействия сверхнизких (4,2 K) температур на наномеханические свойства зубной эмали // Медицинская физика - 2010: сб. материалов III Евразийского конгресса по медицинской физике и инженерии. - 2010. - М.: Изд-во МГУ, 2010. - Т. 4. Инновационные технологии в медицине и здравоохранении. - Ч. I. - С. 323-324.
- Най Д. Физические свойства кристаллов. - М.: ИЛ, 1960. - 385 с.
- Перцов С.С., Стюрева Г.М., Муслов С.А., Синицын А.А., Корнеев А.А., Зайцева Н.В. Основы биомеханики для стоматологов. - М.: МГМСУ, 2017. - 115 с.
- Хышиктуев Б.С., Хышиктуева Н.А. Биохимия полости рта: учеб. пособие. - Чита: ИИЦ ЧГМА, 2004. - 84 с.
- Goldstein R.V., Gorodtsov V.A., Komarova M.A., Lisovenko D.S. Extreme values of the shear modulus for hexagonal crystals // Scripta Mater. - 2017. - Vol. 140. - P. 55-58. doi: 10.1016/j.scriptamat.2017.07.002
- Katz J.L., Ukraincik K. On the anisotropic elastic properties of hydroxyapatite // J. Biomech. - 1971. - Vol. 4, iss. 3. - P. 221-227. doi: 10.1016/0021-9290(71)90007-8
- Katz J.L. Hard tissue as a composite material - I. Bounds on the elastic behavior // J. Biomech. - 1971. - Vol. 4, iss. 5. - P. 455-473. doi: 10.1016/0021-9290(71)90064-9
- Lees S., Davidson C.L. The role of collagen in the elastic properties of calcified tissues // J. Biomech. - 1977. - Vol. 10, № 8. - P. 473-486.
- Lees S., Rollins F.R. Jr. Anisotropy in hard dental tissues // J. Biomech. - 1972. - Vol. 5, № 6. - P. 557-664.
- Povolo F. Poisson's ratio in zirconium single crystals // J. of Nuclear Materials. - 1983. - Vol. 118, iss. 1. -P. 78-82.
- Prium G.J., De Jongh H.J., Ten Bosch J.J. Forces acting on the mandible during bilateral static bite at different bite force levels // J. Biomech. - 1980. - Vol. 13. - P. 735-763.
Statistics
Views
Abstract - 58
PDF (Russian) - 219
Refbacks
- There are currently no refbacks.