STUDY OF THE STRENGTH OF A HIP ENDOPROSTHESIS MADE OF POLYMERIC MATERIAL

Abstract


Developing individualized implant designs in the shortest possible time is a challenging problem in personalized medicine; advanced digital technologies are increasingly integrated in biomedicine to tackle this challenge. While acceptable results are obtained in hip arthroplasty, there is still a number of serious issues that require further research into the design and materials of personalized implants. In particular, studies on the behaviour of biomechanical systems under functional loads have raised a number of questions concerning correct application of loads accounting for the muscle-joint interactions and the implant material selected that could ensure long-term biocompatibility with the human body without replacement. As a way to solving the problem, we considered personalized finite element models prepared based on patient computer tomography data obtained from NMRC for Traumatology and Orthopaedics named after R.R. Vreden (St. Petersburg). We constructed a realistic finite element model of the pelvis, investigating the stress-strain state of the pelvic bones and implants during simulation of primary and revision hip arthroplasty. A series of simulations was performed for different shapes of the endoprosthesis and fixation methods, analysing the potential applications for the personalized endoprosthesis made of a polymer material.

Full Text

Эндопротезирование тазобедренного сустава является одной из самых распространенных хирургических операций и одной из фундаментальных задач ортопедической биомеханики. Несмотря на значительный опыт, накопленный во врачебной практике, и большое количество теоретических и экспериментальных исследований искусственного тазобедренного сустава, вопросы прочности и длительной стабильности эндопротеза остаются открытыми. Успешная клиническая эффективность тотального эндопротезирования суставов может определяться многими факторами, в том числе биомедицинскими, биомеханическими и трибологическими соображениями. Установлено, что износ и биологическая реакция на него играют важную роль в долговечности имплантатов. Выбор материала имплантата делается на основе его механических, технологических, физических и химических свойств. Как правило, конструкции стандартного и персонализированного эндопротезов, а также винты изготавливают из титанового сплава Ti-6Al-4V [1], который является биосовместимым материалам с высоким потенциалом к остеоинтеграции и имеет высокие прочностные характеристики [5]. Для обеспечения биологической фиксации имплантата поверхность имплантата, контактирующая с костью, подвергается дополнительной обработке с целью формирования пористой структуры. Несмотря на ряд преимуществ титана как весьма популярного в медицине материала и значительный опыт применения титановых сплавов, нельзя однозначно утверждать, что характеристики сплошного материала полностью наследуют материалы и конструкции, получаемые с помощью аддитивных технологий. Современные исследования показывают значительный разброс механических свойств и снижение предела усталости для титановых сплавов, получаемых с помощью селективного лазерного плавления или аналогичными методами [17]. Другой существенной проблемой титановых имплантатов, обладающих не только высоким пределом текучести, но и большими значениями модулей упругости по сравнению с костной тканью [27], что приводит к так называемому эффекту экранирования напряжений, который впервые был замечен вокруг ножки эндопротеза в интрамедулярном канале бедренной кости [19]. За счет биомеханической адаптации костной ткани к более низкому уровню напряжений вокруг ножки эндопротеза начинается процесс резорбции живой ткани, что и приводит в результате к снижению ее прочности и разрушению вблизи имплантата. Данный эффект рассасывания костной ткани представляется также существенным для некоторых типов индивидуальных конструкций ацетабулярных компонентов эндопротезов из титанового сплава. В частности, для эндопротезов рожкового типа, имеющих в качестве основного крепежного элемента стержень, входящий в специально подготовленное отверстие в гребне тазовой кости [2; 3]. Для уменьшения эффекта экранирования и поддержания напряжений в костной ткани вокруг эндопротеза на физиологическом уровне, соответствующем напряжениям в неповрежденной бедренной кости, предлагается использовать биосовместимые композитные материалы для производства новых типов бедренных компонентов искусственного сустава [6; 32]. При этом существенное преимущество композитных материалов состоит в том, что путем вариации механических составляющих композита можно регулировать его жесткость в широких пределах при достаточно высокой постоянной прочности материала [12]. В случае стандартного решения при первичном эндопротезировании как титановые, так и полимерные вертлужные компоненты, как правило, крепятся с использованием костного цемента. Помимо этого, установка вертлужного компонента может осуществляться методом бесцементной механической фиксации чашки за счет разницы диаметров костного ложа и чашки (технология press-fit). При этом пористое покрытие имплантата стимулирует эффективную остеоинтеграцию [28] и обеспечивает долговременную фиксацию эндопротеза в костных структурах таза. Можно также отметить, что выбор, как правило, остается за лечащим врачом и основывается на собственном опыте и медицинских рекомендациях. Отметим, что сложность лечения обширных дефектов вертлужной впадины в ревизионном эндопротезировании тазобедренного сустава подтверждается и тем, что существует множество вариантов реконструкции, ни один из которых не обладает явным преимуществом перед остальными [2; 41]. Следовательно, существует необходимость в оценке стабильности и работоспособности разрабатываемых конструкций, которая определяется напряженно-деформированным состоянием биомеханической системы при функциональных нагрузках. Здесь необходимо понимать, что стабильность суставных имплантатов будет зависеть от всех заранее учтенных условий нагружения, в том числе и от мышечных усилий, которые вносят свой вклад в напряженное состояние системы «кость - имплантат». Компьютерное моделирование биомеханического поведения различных отделов человеческого тела получило широкое развитие именно благодаря современным программным комплексам, позволяющим проектировать трехмерные модели костей конкретного пациента на основании данных томографии. В свою очередь 3D-печать дает возможность производить медицинские изделия под конкретного пациента. Данное направление персонализированной медицины является крайне важным в настоящее время, что приводит к необходимости поиска новых методов проектирования и производства. Таким образом, в работе будут представлены результаты исследования напряженно-деформированного состояния биомеханической системы «таз - полимерный эндопротез с пористым покрытием» при использовании стандартного и индивидуального медицинских решений. Материалы и методы Постановка задачи В работе рассматривается модель биомеханической системы, образованной костными структурами таза и а б Рис. 1. Поверхностная цифровая модель таза пациента, построенная на основе компьютерной томографии: а - послеоперационный период для стандартной операции первичного эндопротезирования; б - послеоперационный период при установке индивидуального имплантата персонализированной конструкцией имплантата, установленного в результате операции повторного эндопротезирования тазобедренного сустава пациентке (рост 164 см, вес 94 кг), проходившей лечение в ФГБУ «НМИЦ травматологии и ортопедии имени Р.Р. Вредена» (Санкт-Петербург). С помощью компьютерной томографии (КТ) и стандартного программного обеспечения обработки снимков КТ получены пространственные геометрические модели в формате STL здоровой половины таза и, соответственно, тазобедренного сустава. Виртуальная модель целого таза была получена путем симметричного отображения здоровой стороны и состояла из соединенных между собой тазовых костей и крестца. Это было необходимо для сравнительного анализа двух видов операции: первичной и ревизионной, требующей замену вертлужной впадины и окружающих костных тканей на искусственную структуру (рис. 1). Отметим, что анатомически таз человека состоит из тазовых костей и крестца, соединенных между собой суставами пояса нижних конечностей. Поэтому с целью более корректного моделирования системы на этапе создания конечно-элементной сетки поверхностные геометрические модели таза были разделены на отдельные анатомические костные структуры. При исследовании поведения системы при функциональных нагрузках возникает ряд вопросов, касающихся, корректного задания условий нагружения, которые учитывали бы взаимосвязь мышц и суставов. Основываясь на опыте российских и зарубежных исследователей [1; 3; 10; 30], было принято решение об использовании экспериментальных значений сил, возникающих в костно-мышечной системе при различных видах движений человека, представленных на открытом ресурсе OrthoLoad (https://orthoload.com/) и в программе HIP98 [4]. Задача формулируется в квазистатической постановке на основе уравнений инфинитезимальной теории упругости, предполагающими малые перемещения и деформации. Принято, что все входящие в систему «таз - эндопротез» ткани и элементы конструкции описываются моделью линейно упругой сплошной среды. Тем не менее постановка задачи подразумевает нелинейные эффекты вследствие наличия контактного взаимодействия с трением между костными структурами и компонентами эндопротеза. Более подробно контактные области, типы контактов и значение коэффициента трения приведены далее в разделе, посвященном описанию конечно-элементной модели системы. Физико-механические свойства костных тканей Особое внимание при анализе напряженного состояния любых механических и, особенно, биомеханических систем, следует уделить физико-механическим свойствам используемых материалов. Тазовая кость образована губчатым веществом низкой плотности, формирующим спонгиозную (трабекулярную) биологическую ткань, и компактным веществом высокой плотности, формирующим кортикальную кость. Большая часть нагрузки передается через внешний достаточно тонкий кортикальный слой, а спонгиозная ткань является связующим звеном, удерживающим компактное вещество от потери устойчивости и разрушения. Благодаря этому тазовая структура сочетает в себе малый вес и высокую прочность конструкции. Проведенный литературный обзор демонстрирует, что представленные в публикациях данные имеют значительный разброс. В силу возраста и других причин, которые могут вызвать деградацию костной ткани, механические свойства могут изменяться в широком диапазоне [1; 8; 23; 30]. Губчатое и компактное вещество тазовых костей человека демонстрирует анизотропию упругих свойств [3; 8; 32; 40], хотя и в меньшей степени, чем аналогичные вещества длинных трубчатых костей скелета, а также пространственное распределение плотности и модуля Юнга. Однако изменение модулей упругости как кортикальной, так и спонгиозной тканей таза незначительно как было показано ранее [8; 9]. Во многих работах, посвященных конечно-элементному моделированию тазовых костных структур, кортикальный слой и спонгиозная ткань описываются моделью линейно-упругой изотропной однородной сплошной среды [10; 14; 15; 17; 19; 21; 36; 39]. После анализа первоисточников выбраны наиболее репрезентативные значения механических свойств (табл. 1). Таблица 1 Физико-механические свойства биологических тканей Материал Модуль упругости, МПа Коэффициент Пуассона Плотность, кг/м3 Предел текучести, МПа Кортикальная кость [9; 32] 17000 0,3 1700 150 [32] Спонгиозная кость [8; 22] 70 0,2 200 2,1 [22] Лобковый симфиз [5; 13; 34] 5 0,495 2 - Хрящ крестцово-подвздошных суставов [5; 13; 34] 350 0,495 1,05 - Важно отметить, что условный предел текучести, характеризующий начало необратимых деформаций в цикле «нагружение - разгрузка», для губчатого вещества, о котором сообщалось в публикациях [8], кажется чрезмерно высоким. В работе [22] экспериментально показано, что с клинической точки зрения спонгиозная кость может быть нагружена до значений деформации в 2-3 % без видимой пластической деформации, повреждений или переломов. Поэтому в данном исследовании в качестве предела прочности спонгиозной кости использован диапазон значений 1,4-2,1 МПа. Данные [22] были подтверждены исследованиями авторов на цилиндрических образцах из губчатого вещества, выделенных из головки бедренной кости после операции [23]. Физико-механические свойства материала имплантата В работе с помощью разработанной цифровой модели исследована возможность применения органического биосовместимого полимера полиэфирэфиркетона (PEEK) [28] в качестве материала имплантата с точки зрения обеспечения необходимой прочности всей биомеханической конструкции. Полимерный материал PEEK и композиты на его основе начали использоваться в биомедицине в 90-х гг. прошлого века. Благодаря своей устойчивости к усталостным нагрузкам полиэфирэфиркетон в настоящее время широко используется при операциях на позвоночнике [37]. Основным достоинством данного материала является то, что за счет контроля на этапе армирования другими высокопрочными биосовместимыми материалами прочность и дизайн имплантата могут быть изменены и адаптированы к механическому и физиологическому состоянию тканей пациента. Еще одним преимуществом использования PEEK и его композиционных производных является переменная настраиваемая жесткость, которой обычно способствует добавка в виде углеродных волокон [6; 24]. При этом в зависимости от концентрации углеродных волокон физико-механические свойства композита на основе PEEK могут быть существенно изменены [37]. Модуль упругости PEEK-композитов, армированных углеродным волокнами (PEEK/CF), может быть адаптирован для имитации свойств компактного костного вещества (18 ГПа) или титанового сплава (110 ГПа), а также может воспроизводить модули кортикальной и спонгиозной ткани. Поэтому предполагается, что композитные материалы могут дать больший потенциал структурной биосовместимости, чем однородные монолитные материалы. Помимо этого, к достоинствам стоит отнести отсутствие коррозионно-усталостного разрушения полиэфирэфиркетона, что в перспективе исключает расшатывание имплантата и иммунные реакции у пациента. Механические свойства композитов PEEK, армированных ориентированными углеродными волокнами, очевидно, наследуют анизотропию их пространственного распределения. Модуль продольного растяжения и прочность однонаправленного армированного длинными волокнами композита, в котором непрерывные волокна выровнены параллельно направлению нагрузки, значительно выше, чем модуль поперечного растяжения и соответствующая прочность (волокна перпендикулярны направлению растягивающей нагрузки). Напротив, механические свойства композитов PEEK, армированных случайно ориентированными короткими углеродными волокнами (SCF), являются изотропными. Короткие волокна SCF намного дешевле, чем длинные непрерывные углеродные волокна, а композит PEEK/SCF можно легко изготовить с использованием обычных методов экструзии и литья под давлением [25]. Композит PEEK с 30 %-ным содержанием SCF обладает жесткостью 24 ГПа и пределом прочности 214 МПа [16] и может применяться для изготовления компонентов эндопротезов тазобедренного сустава [6]. В ходе проведенного анализа мировых первоисточников сделан вывод, что «чистый» PEEK неоднократно применялся при операциях на позвоночнике, а также есть упоминания о возможности использования неармированного полиэфирэфиркетона в качестве чашки ацетабулярного компонента искусственного тазобедренного сустава. Однако отсутствуют публикации, в которых было бы рассмотрено применение неармированного PEEK для изготовления тазового компонента индивидуального имплантата тазобедренного сустава и приведен анализ его напряженно-деформированного состояния при функциональных нагрузках. Поэтому в рамках данной работы анализ будет вестись именно с его использованием. Физико-механические свойства используемых в расчетах материалов приведены в табл. 2. Отметим, что поверхность имплантата, контактирующая с костной тканью, покрыта пористой структурой, обеспечивающей со временем эффективную остеоинтеграцию [28]. Прямое конечно-элементное моделирование такой структуры сильно усложнит численную Таблица 2 Физико-механические свойства компонентов эндопротеза Материал Модуль упругос-ти, МПа Коэффициент Пуассона Предел текучести, МПа Неармированный PEEK [37] 3600 0,42 92 Пористый PEEK (68 %) 370 0,25 63 модель системы и потребует больших вычислительных мощностей при решении задачи. Во избежание этого пористый полимер представлен в виде однородной изотропной среды, имеющей упругие характеристики, эквивалентные макроскопическим упругим характеристикам пористой структуры (табл. 2). В биомеханическую систему «таз - эндопротез» включены также дополнительные компоненты для связи костей и имплантата. Такими связующими компонентами являются винты для крепления имплантата для индивидуального решения, а также фиксирующая ножка имплантата рожкового типа [2]. Материал винтов подразумевается тем же, что и основная часть эндопротеза, т.е. неармированный PEEK (табл. 2). Конечно-элементные модели Тазовые кости, крестец и эндопротез имеют сложную пространственную геометрию (рис. 2), что приводит к значительным трудностям в построении конечно-элементных моделей биологических и искусственных структур. Геометрические пространственные модели таза были построены в специализированной системе для сегментирования данных томографии на основе серии рентгенологических снимков в формате DICOM, полученных с помощью компьютерной томографии нового поколения с высокой разрешающей способностью. Затем первичные геометрические модели подвержены постобработке и позиционированию и сохранены в формате STL в виде поверхностной треугольной сетки (рис. 3). Для доработки геометрических моделей и создания объемной конечно-элементной сетки были использованы специализированные программные комплексы для построения конечно-элементных сеток и препроцессинга расчетных моделей, которые предлагают широкий выбор инструментов для подготовки и улучшения трехмерных поверхностных сеток. В ходе обработки данных КТ оказалось, что тазовые кости и крестец представляют собой одно тело. Поскольку анатомически это не так, то модель таза была разделена на отдельные анатомические кости с введением в рассмотрение хрящевой прослойки, моделирующей связывающие их сочленения. На основе рассмотренных геометрических моделей разработаны математические численные модели компонентов системы с использованием программы конечно-элементного анализа SIMULIA Abaqus (рис. 4, 5) для анализа стандартной и индивидуальной операций. В качестве основного конечного элемента в соответствии с описанной выше постановкой задачи применен трехмерный тетраэдральный линейный четырехузловой элемент (C3D4 из библиотеки Abaqus). а б Рис. 2. Геометрические модели оперируемой правой половины тазовых костных структур пациента в случае: а - стандартного первичного эндопротезирования решения; б - индивидуального решения Рис. 3. Поверхностная треугольная сетка исследуемой модели а б Рис. 4. Конечно-элементные модели: а - тазовых костей; б - стандартной конструкции имплантата из полимерного материала а б Рис. 5. Конечно-элементные модели: а - индивидуальной конструкции имплантата из полимерного материала; б - ножки имплантата и фиксирующего винта Суммарное количество элементов в системе «таз - стандартный имплантат» составляет 1 539 271, узлов - 338 339. Общее число степеней свободы - 1 027 431. Суммарное количество элементов в системе «таз - индивидуальный имплантат» составляет 2 810 322, узлов - 636 724. Общее число степеней свободы - 1 910 862. При разделении геометрии костных структур на кортикальную и спонгиозную ткани толщина кортикального слоя была выбрана фиксированной на уровне 1,5 мм (рис. 4, а), что соответствует результатам анализа КТ таза и реферативным данным [9]. В качестве конструкций имплантата были рассмотрены стандартная полусферическая вертлужная чашка (рис. 4, б) и индивидуальная конструкция (рис. 5), изготовленные из сплошного полимера с пористым внешним слоем толщиной 1,5 мм. В задаче возникают нелинейности вследствие наличия контактных взаимодействий между костными структурами и компонентами эндопротеза. Контактные области с трением находятся в местах соединения имплантата и костей. Коэффициент трения принят равным 0,35 [3]. Нагрузки и кинематические ограничения В соответствии с постановкой задачи предполагается, что при ходьбе на систему «таз - эндопротез» действуют реактивные силы в центре тазобедренного сустава и усилия от мышц, связанных с тазом, также имеет место полный запрет трансляционных перемещений проксимальной поверхности крестца (рис. 6). Кроме того, в случае индивидуальной конструкции имплантата (рис. 5, а), изготовленного из сплошного полимера, при наличии пористого слоя его фиксация осуществляется медицинскими винтами (рис. 5, б). Сила затяга винтов равна 250 Н в силу настройки медицинского инструмента, а ее численное моделирование описано в [3]. Анализ стандартного эндопротеза проводится для следующих комбинаций «тип фиксации - материал вертлужной чашки»: 1) плотная бесцементная механическая фиксация (press-fit) чашки из сплошного полимера внутри при наличии внешнего пористого слоя; 2) фиксация press-fit чашки из полностью пористого полимера (рис. 4, б). В перечисленных случаях винты для крепления чашки не используются. Рис. 6. Схема нагружения модели таза в глобальной системе координат при движении человека: FL и FR - векторы сил реакций, возникающие при движении человека в левом и правом тазобедренных суставах; ML и MR - векторы сил, развиваемых левой и правой отводящими мышцами при движении человека Рис. 7. Графики проекций сил реакций тазобедренного сустава при одном шаге в глобальной системе координат для левой и правой ноги: FLx, FLy, FLz - силы реакций, возникающие в левом тазобедренном суставе; FRx, FRy, FRz - силы реакций, возникающие в правом тазобедренном суставе Рис. 8. Графики проекций мышечных сил при одном шаге в трех проекциях для левой и правой ноги, MLx, MLy, MLz, MRx, MRy, MRz - мышечные силы, возникающие в левой и правой ногах соответственно, в трех проекциях глобальной системы координат Также в работе [3] подробно представлено задание реактивных сил при ходьбе человека на основе данных ресурса OrthoLoad (https://orthoload.com/) и программы HIP98 [4]. Для рассматриваемого в настоящей работе пациента в соответствии с его весом 940 Н рассчитаны силы реакций в трех проекциях глобальной системы координат (рис. 7). Скелетно-мышечные нагрузки в тазобедренном суставе определяются контактной силой в суставе и силами более 20 различных мышц, охватывающих сустав [31]. Существует ряд публикаций, где доказывается, что мышцы вносят значительный вклад в общую нагрузку на скелет человека [11; 25; 35]. Для уменьшения количества точек приложения и исключения мышц, не вносящих существенного вклада, была использована концепция «мышц одного сустава», то есть мышц, которые охватывают исключительно тазобедренный сустав. В этом упрощении все волокна ягодичных мышц с аналогичной функцией (большая, средняя и малая ягодичная мышцы) сгруппированы в одну «отводящую мышцу» с одним местом прикрепления. Анализ экспериментальных данных, представленных на ресурсе OrthoLoad (https://orthoload.com/) и в программе HIP98 [4], показал, что наибольший силовой вклад в нагрузку на кости таза вносит «отводящая» мышца, поэтому именно она включена в расчет напряженного состояния рассматриваемой системы. Аналогично заданию сил в тазобедренном суставе для отводящей мышцы использованы данные репрезентативного осредненного субъекта. Для рассматриваемого в работе конкретного пациента в соответствии с его весом 940 Н рассчитаны мышечные силы в осях глобальной системы координат (рис. 8). Результаты и обсуждение Полученные результаты описывают решение задачи с учетом взаимосвязи усилий мышц и суставов в двух конструктивных вариантах - стандартный имплантат при первичном эндопротезировании и индивидуальная конструкция при повторной операции и сильно разрушенной тазовой кости. Результаты расчетов представлены в материале имплантата и костных тканях. Рассматривались эквивалентные по Мизесу напряжения в контактных парах: имплантат - кортикальная костная ткань, имплантат - спонгиозная костная ткань, а также в парах: имплантат - винт, кость - винты. Изменение напряжений при ходьбе представлено полями и графиками в характерных точках на границе раздела материалов, описывающих кортикальный слой и спонгиозную ткань тазовых костей, сплошной и пористый слои вертлужной чашки. Такими точками являются: точка А - в правой нижней части костного ложа и чашки, точка Б - в левой нижней части костного ложа и чашки, точка В - в верхней части на кромке костного ложа и чашки. Поле максимальных напряжений, возникающих в узлах модели за цикл нагружения, при стандартном эндопротезировании показано на примере второго расчетного случая (рис. 9). Нижними индексами «к» и «г» обозначены компактное и губчатое костное вещество кортикальной и спонгиозной тканей. На рис. 10 приведены графики распределения эквивалентных напряжений по Мизесу в кортикальной и спонгиозной костной ткани при ходьбе в случае стандартного решения (первичное эндопротезирование). Заметим, что максимальные эффективные напряжения по Мизесу в кортикальном слое тазовых костей не превышают условный предел текучести (рис. 10, а). Однако имеет место превышение напряжений условного предела текучести в рассматриваемых зонах для спонгиозной ткани (рис. 10, б). Можно предположить, что в области контакта внешней поверхности чашки со спонгиозной тканью костного ложа возможны локальные разрушения губчатого вещества при ходьбе пациента. Поскольку проверка качества конечно-элементной сетки показала хорошие результаты, то такое поведение материала нельзя считать численной ошибкой. Для более адекватной интерпретации полученных результатов необходимо оценить, сколько циклов нагрузки «выдержит» материал до непосредственного своего разрушения. Поскольку зона повышенных напряжений является относительно небольшой, то процессы начальной регенерации и постоянного ремоделирования спонгиозной ткани могут предотвратить глобальное разрушение кости в области вертлужной впадины. Особый интерес для анализа представляет имплантат. Напомним, что в качестве материала стандартной чашки ацетабулярного компонента с бесцементной фиксацией в вертлужной впадине рассматривается органический полимер полиэфирэфиркетон (PEEK) без армирования углеродными волокнами в виде сплошной пористой или полностью пористой структуры. Распределение эквивалентных по Мизесу напряжений, возникающих в пористом и сплошном слоях вертлужной чашки, продемонстрировано на рис. 11 для варианта чашки из двух слоев: сплошного и пористого (случай 1). Расположение точек в области чашки соответствует ранее показанным точкам на поверхности костных тканей (рис. 9). При этом нижними индексами «с» и «п» обозначены сплошной и пористый материал чашки (рис. 11). а б Рис. 9. Распределение максимальных эквивалентных по Мизесу напряжений за цикл ходьбы в случае стандартного решения для второго расчетного случая, МПа: а - в кортикальном слое кости; б - в спонгиозной костной ткани а б Рис. 10. Зависимости эквивалентных по Мизесу напряжений от фазы ходьбы в случае стандартного решения: а - в кортикальном слое кости; б - в спонгиозной костной ткани а б Рис. 11. Распределение максимальных эквивалентных по Мизесу напряжений за цикл ходьбы в случае стандартного решения для первого расчетного случая, МПа: а - в пористом слое стандартной чашки; б - в сплошном материале стандартной чашки На рис. 12, 13 показаны графики эквивалентных по Мизесу напряжений при ходьбе в характерных точках вертлужной чашки для рассматриваемых случаев. На протяжении всего цикла ходьбы максимальные напряжения не превышают предельных показателей как в сплошном, так и пористом слоях вертлужной чашки (см. рис. 12, 13). Графики демонстрируют, что имплантат обеспечивает адекватное напряженно-деформированное состояние системы в целом, при этом серьезных повреждений как чашки, так и костных тканей таза не ожидается. Таким образом, результаты анализа напряженного состояния системы «таз - эндопротез» показывают, что рассмотренные варианты имплантата из полимерного материала стандартной конструкции с фиксацией press-fit являются приемлемыми. Хотя локальные разрушения в зоне выреза кости под вертлужную чашку не исключены, но сам протез при этом остается практически ненагруженным, и, как следствие, в данном случае нет необходимости дополнительно совершенствовать эндопротез в плане повышения его прочностных характеристик, добавляя в материал углеродные волокна. При этом метод фиксации press-fit показал хорошие результаты. На рис. 14-17 приведены аналогичные результаты для индивидуальной конструкции эндопротеза тазобедренного сустава. Используются номера отверстий под медицинские винты, а также прежние обозначениянижних индексов: «к» и «г» - компактное и губчатое вещество кортикальной и спонгиозной тканей. Кроме того, буквами обозначены области: П - подвздошной кости, Л - лонной кости, С - область отверстия в подвздошной кости под ножку (фиксирующий стержень) имплантата. Рис. 12. Зависимости эквивалентных по Мизесу напряжений от фазы ходьбы в пористом слое эндопротеза: второй и третий расчетные случаи Рис. 13. Зависимости эквивалентных по Мизесу напряжений от фазы ходьбы в сплошном слое эндопротеза: первый и второй расчетные случаи Рис. 14. Распределение максимальных эквивалентных по Мизесу напряжений за цикл ходьбы в случае индивидуального решения в кортикальном слое кости, МПа Рис. 15. Распределение максимальных эквивалентных по Мизесу напряжений за цикл ходьбы в случае индивидуального решения в спонгиозной костной ткани, МПа Рис. 16. Зависимости эквивалентных по Мизесу напряжений от фазы ходьбы в кортикальном слое кости пациента в случае индивидуального решения Рис. 17. Зависимости эквивалентных по Мизесу напряжений от фазы ходьбы в спонгиозной костной ткани пациента в случае индивидуального решения Рис. 18. Распределение максимальных эквивалентных по Мизесу напряжений за цикл ходьбы в пористом слое индивидуальной конструкции имплантата, МПа Рис. 19. Распределение максимальных эквивалентных по Мизесу напряжений за цикл ходьбы в сплошном материале индивидуальной конструкции имплантата, МПа Оценивая поведение костной ткани (рис. 16, 17), заметим, что в как в кортикальном слое, так и в спонгиозной ткани достигается «условный» предел текучести. Стоит отметить, что превышения пределов в момент затяга винтов носят локальный характер и определенно связаны с численными эффектами. Однако при максимальной нагрузке в момент ходьбы видно, что отверстия под винты оказываются нагруженными. Чтобы не загромождать графики (рис. 16), для оценки изменения напряжений во время ходьбы были исключены ненагруженные винты. Отметим, что все рассматриваемые зоны за исключением контакта нижней части тазовой кости с протезом имеют высокие значения напряжений, близкие к предельным показателям или превышающие их. Анализ напряжений, возникающих в губчатом веществе кости, проведен в зонах контакта элементов конструкции эндопротеза (винта, дополнительного Рис. 20. Зависимости эквивалентных по Мизесу напряжений от фазы ходьбы в пористом слое индивидуальной конструкции имплантата Рис. 21. Зависимости эквивалентных по Мизесу напряжений от фазы ходьбы в сплошном материале индивидуальной конструкции имплантата элемента имплантата, собственно имплантата) со спонгиозной костной тканью. Чтобы не загромождать график (рис. 17), показаны результаты только в тех зонах, которые оказались самыми нагруженными еще на этапе затяга винтов. Аналогично распределению напряжений в кортикальном слое тазовой кости напряжения во всех рассмотренных зонах спонгиозной ткани либо близки к предельному значению, либо превышают его. Поскольку в данной задаче рассматриваются кратковременные нагрузки, то, возможно, наблюдаемая картина напряжений не будет критична для кости. Отметим также, что спонгиозная костная ткань как биологическая структура имеет способность к восстановлению и такое распределение напряжений не будет критичным. Тем не менее можно сделать предварительный вывод, что согласно проведенному анализу «опасной» зоной будет контакт верхней части костной ткани с протезом. Индивидуальная конструкция выполнена из органического полимера, поверхности которого, контактирующие с костью, покрыты пористой структурой. Ниже приведен анализ напряжений по отдельности в пористом покрытии и в сплошном материале протеза (рис. 18-21). Рассматриваются зоны отверстий под винты и фиксирующий стержень, или ножка, эндопротеза, а также область контакта имплантата с костью. Используются аналогичные предыдущему номера отверстий под винты, обозначения нижних индексов («с» и «п» - сплошной и пористый материал имплантата), обозначения областей имплантата (С - область отверстия под ножку имплантата, К - специфическая область контакта сферической поверхности имплантата с кортикальной тканью). Анализируя графики напряжений при ходьбе в пористом слое индивидуальной конструкции (рис. 20), заметим, что предельный показатель 63 МПа в рассматриваемых зонах имплантата не достигается. Тем не менее область контакта с костью дает самые высокие значения напряжений. Причиной высоких напряжений в зоне контакта может служить дисконгруэнтность контактирующих поверхностей, т.е. их неполное геометрическое соответствие. В таком случае нарушение топологической совместимости двух поверхностей ведет к деформации соединения, что препятствует нормальному функционированию и, как следствие, приводит к высоким напряжениям. Анализируя графики напряжений в сплошном материале индивидуальной конструкции (рис. 21), отметим, что напряжения, возникающие в зонах отверстий области имплантата под винты и ножку эндопротеза, превышают установленный предел прочности на протяжении значительного интервала цикла ходьбы. Возможно, как и в рассмотренных ранее случаях, это связано с совместимостью поверхностей, а также с «локальными» всплесками напряжений на острых краях отверстий, приводящих к сингулярности напряжений. Данная теоретическая проблема известна в линейной механике сплошных сред, и подобные эффекты отмечались в предыдущих работах авторов [3; 27]. Принято, в том числе в ходе консультаций с медицинскими специалистами, что если зона концентрации напряжений не слишком обширная, то они могут считаться допустимыми с точки зрения инфинитезимальной теории упругости. Заключение В работе рассмотрена одна из важных проблем клинической биомеханики опорно-двигательного аппарата: исследование напряженно-деформированного состояния системы «скелет - полимерный имплантат» тазобедренного сустава при ходьбе. Несмотря на приемлемые результаты при эндопротезировании тазобедренного сустава, все еще существует ряд серьезных проблем, требующих дальнейшего изучения конструкции и материалов кастомизированных имплантатов. В частности, при исследовании поведения системы при функциональных нагрузках возникает ряд вопросов, касающихся корректного задания нагрузок, которые учитывали бы взаимосвязь мышц и суставов, и выбора материала имплантата, который мог бы обеспечивать долгосрочную биосовместимость с организмом человека. На основе данных компьютерной томографии пациента разработана реалистичная цифровая модель таза и смоделированы две ситуации, имитирующие первичную и повторную операции эндопротезирования тазобедренного сустава. В результате компьютерного анализа с использованием универсального конечно-элементного комплекса инженерного анализа SIMULIA Abaqus рассчитано распределение напряжений в полимерных компонентах эндопротеза тазобедренного сустава и костных тканях таза пациента. На основе полученных результатов и литературных данных по допустимым значениям напряжений проанализирована прочность костных структур таза и эндопротеза из полимерного материала. Система фиксации при использовании полимерного имплантата в случае первичного эндопротезирования показала удовлетворительные результаты. Однако анализ индивидуального решения показал превышение предельных напряжений во всех рассматриваемых контактных парах. Можно предположить, что при больших дефектах тазобедренного сустава и высоких уровнях нагрузки более предпочтительным будет полимер, армированный углеродными волокнами. Таким образом, актуальным остается вопрос разработки таких конструкций и материалов для эндопротезирования тазобедренных суставов, которые могли бы обеспечивать надежную остеоинтеграцию и функционирование имплантата в организме без замены в течение длительного времени.

About the authors

L. B Maslov

Ivanovo State Power Engineering University; Peter the Great St. Petersburg Polytechnic University

A. Y Dmitryuk

Center for Technical Projects

M. A Zhmaylo

Peter the Great St. Petersburg Polytechnic University

A. N Kovalenko

National Medical Research Center for Traumatology and Orhopaedics named after R.R. Vreden

References

Statistics

Views

Abstract - 106

PDF (Russian) - 107

Refbacks

  • There are currently no refbacks.

Copyright (c) 2022 Маслов Л.Б., Дмитрюк А.Ю., Жмайло М.А., Коваленко А.Н.

This website uses cookies

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

About Cookies