The determination of spatial stress-strain state of human proximal femur

Abstract


Work is dedicated to the problem of modeling the structural adaptation of spongy bone tissue of the proximal femur after hip surgery osteosynthesis. In this paper, the elastic problem of determining of equilibrium strain field in the proximal femur during single support phase step in the norm. Equilibrium strain is strain at which there is no adaptation of bone tissue. Individual model of the proximal femur with an inhomogeneous distribution of mechanical properties in each finite element of model is constructed. Fields of equilibrium strain, which do not contradict to known experimental data, both qualitatively and quantitatively, are calculated with finite element method. The equilibrium strain are necessary for the formation of local incentives to adapt in the future.

Full Text

Перелом шейки бедренной кости - серьезная травма, которая при неправильном или несвоевременном лечении может привести к тяжелым последствиям, связанным с полным протезированием сустава, или летальному исходу. В связи с ростом числа дорожно-транспортных происшествий увеличилось число переломов шейки бедра у людей в молодом возрасте, что ставит проблему лечения в разряд социальных. Традиционный способ лечения переломов - остеосинтез (соединение костей) с помощью резьбовых фиксаторов, устанавливаемых хирургами с максимальным моментом закручивания (рис. 1). Рис. 1. Фиксация отломков проксимального отдела бедра резьбовыми фиксаторами Современная технология предполагает учет индивидуальных свойств костной ткани пациента путем предоперационного рентгенологического исследования и компьютерного моделирования адаптации в процессе сращения. Технология является контролируемым остеосинтезом шейки бедра [1]. Однако в этой технологии используемое одноосное напряженно-деформированное состояние не позволяет учесть влияние таких факторов, как особенность индивидуальной формы и размеров, пространственной неоднородности распределения структурных и механических свойств костной ткани проксимального отдела бедра пациента, а также вид перелома и ориентация фиксаторов в кости на величину индивидуального момента закручивания фиксаторов. В качестве исследуемой области берется фрагмент проксимального отдела бедра (рис. 2). На границах области задаются перемещения и нагрузки. Все тензорные величины и их инварианты, входящие в приведенные ниже уравнения, являются функциями радиуса-вектора частицы М. Рис. 2. Рассматриваемая область кости - проксимальный отдел бедра Система уравнений начально-краевой задачи содержит следующие соотношения: - уравнение равновесия (1) предполагает неоднородность напряженно-деформированного состояния фрагмента кости; - определяющее соотношение пороупругости [2] , (2) где - тензор жесткости губчатой костной ткани является функцией от - плотность кости, зависящая от радиуса-вектора [3]; - кинематическое соотношение Коши для тензора малых деформаций ; (3) - граничные условия: а) силовые: (4) где - вектор внешней нормали; - давление ацетабулярной впадины, , (5) где - вектор распределенной мышечной нагрузки б) кинематические: граница фиксируется так, что одна из ее точек (С) закреплена в глобальной системе координат , (6) а остальные точки закреплены только в направлении нормали , (7) где - вектор перемещения частицы. Значения давлений ацетабулярной впадины и нагрузок от мышечных сил представлены в работах Miyanaga и Ueo [4, 5]. С помощью конечно-элементного программного пакета Ansys 12.1 (Ansys, USA) был произведен расчет линейного напряженно-деформированного состояния. На рис. 3 представлено рассчитанное поле равновесных деформаций. Отсюда следует, что наибольшие деформации, равные 0,0012, возникают в менее жестких областях внутри зоны большого вертела, равные 0,0017 - в области серповидного пятна контакта, а также равные 0,0012 - в наиболее нагруженных областях в шейке бедра. Полученное поле деформаций находится в области значений, не противоречащих Рис. 3. Поле равновесных деформации проксимального отдела бедра приведенным экспериментальным данным [6], согласно которым физиологические равновесные деформации изменяются в переделах ±0,003 напряженно-деформированного состояния.

About the authors

A. F. Shuliat'ev

Perm National Research Polytechnic University

Iu. V. Akulich

Perm National Research Polytechnic University

References

  1. Акулич Ю.В., Акулич А.Ю., Денисов А.С. Предоперационное определение индивидуальных физических характеристик губчатой костной ткани проксимального отдела бедра человека // Российский журнал биомеханики. - 2011. - Т. 15, № 1 (51). - С. 33-41.
  2. Определяющие соотношения структурной адаптации костной ткани / Ю.В. Акулич, П.А. Брюханов, М.В. Мерзляков, А.В. Сотин // Известия Саратовского государственного университета. Новая серия. Серия: Математика. Механика. Информатика. - 2011. - Т. 11, № 2. - С. 54-61.
  3. Relations of mechanical properties to density and CT numbers in human bone / J.Y. Rho [et al.] // Medical Engineering and Physics. - 1995. - Vol. 17, no. 5. - P. 347-355.
  4. Miyanaga Y., Fukubayashi H., Kurosava N. Contact study of the hip joint: load deformation pattern, contact area and contact pressure // Archives of Orthopaedic and Traumatic Surgery. - 1984. - Vol. 103, № 13. - P. 13-17.
  5. Biomechanical aspects of the development of aseptic necrosis of the femoral head / T. Ueo [et al.] // Archives of Orthopaedic and Traumatic Surgery. - 1985. - Vol. 104. - P. 145-149.
  6. Регирер С.А., Штейн А.А., Логвенков С.А. Свойства и функции костных клеток: биомеханические аспекты // Современные проблемы биомеханики. Механика роста и морфогенеза. - М.: Изд-во Моск. ун-та, 2000. - № 10. - С. 174-224.

Statistics

Views

Abstract - 7

PDF (Russian) - 2

Refbacks

  • There are currently no refbacks.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies