Лечение пострадавших с наиболее сложными переломами длинных трубчатых костей, к числу которых относятся переломы костей голени, не только продолжает оставаться актуальной проблемой современной травматологии и ортопедии, но и приобретает особое значение в современных условиях нарастающей этажности городской застройки, стремительного развития высокоскоростного транспорта и техногенного давления на человека всё более усложняющейся производственной обстановки [1].
Растущие требования социума к качеству лечения переломов привели за последние несколько десятилетий к разработке современных оперативных методик, обеспечивающих высокую стабильность фиксации за счет биосовместимости использующихся конструкций, совершенствования технико-технологических характеристик процессов их установки и удаления на завершающей стадии консолидации перелома. Успехи в совершенствовании методик лечения переломов достигнуты как результат клинического опыта и применения методов биомеханики. Показательны в этом отношении исследования А.И. Сеппо [8] и Х.А. Янсона [10], выполненные до начала широкого применения компьютерных технологий. В современных условиях методы математического моделирования и информационные технологии многократно увеличивают возможности исследователей, в том числе в области совершенствования методик лечения переломов. К числу наиболее эффективных и универсальных инструментов компьютерного моделирования относится метод конечных элементов [2]. Этот метод успешно применяется в моделировании инженерных конструкций и биомеханических систем [4].
Для сравнительного исследования качества получаемых при различных вариантах остеосинтеза моделей «кость - фиксатор» и «кость - аппарат внешней фиксации» («кость -АНФ») используются, в соответствии с положениями механики, методы, применяемые при решении задач, возникающих в условиях деформации твердого тела. При этом костная ткань рассматривается как изотропный материал.
Конструкторская разработка и оптимизация современных и фьючерсных средств фиксации переломов требуют тщательной экспериментальной оценки качества всех вариантов остеосинтеза согласно требованиям биомеханики. То есть выявления мест возникновения наибольших напряжений в условиях нагрузки, испытываемых твердым телом, приводящих к разрушению костной ткани или фиксирующих устройств (чаще всего в местах их рабочего контакта) [3; 6].
Представленная работа является фрагментом исследования, посвященного проблеме выбора оптимального варианта фиксации различных оскольчатых диафизарных переломов костей голени типа, в частности В2 по классификации AO/ASIF [7], в условиях остеосинтеза блокируемым интрамедуллярным стержнем и аппаратом Г.А. Илизарова. При исследовании напряжений производилось компьютерное моделирование систем «кость - фиксирующая конструкция» и исследовалась их стабильность методом конечных элементов. При создании математической модели использовалась подготовка предварительных моделей с помощью рентгенограмм и скиаграмм переломов, а также серии томограмм, получаемых с помощью мультиспирального компьютерного томографа фирмы SIMENS. Подобное решение проблемы использовалось, по данным литературных источников, и другими исследователями [5; 9].
Учет анизотропных свойств материала кости, решение контактной задачи теории упругости позволяет представить картину распределения напряжений в элементах конструкции с высокой степень достоверности. Компьютерное моделирование с применением метода конечных элементов на основе томограмм, выполненных с шагом 0,5 мм, обеспечивает создание модели большеберцовой кости, подобной биологическому прототипу. Кроме этого, методика позволяет рассматривать напряженно-деформированное состояние фиксатора при функциональной нагрузке.
Сложные геометрические модели биологических объектов создавались в системе трехмерного проектирования по результатам данных компьютерной томографии (рис. 1).
Рис. 1. Моделирование геометрии.
Первым этапом настоящего исследования было создание компьютерной модели оскольчатого перелома типа В2 по классификации АО/ASIF большеберцовой кости на основе данных компьютерной томографии. А также построения моделей блокирующего интрамедуллярного стержня и аппарата Г.А. Илизарова на основе схем и технических характеристик конструкции (рис. 2).
Рис. 2. Построение 3D-модели фиксатора на примере блокируемого интрамедуллярного стержня фирмы Chm (Польша).
Вторым этапом являлось создание сборной компьютерной модели «кость - фиксатор» в условиях остеосинтеза блокирующим интрамедуллярным стержнем и аппаратом Г.А. Илизарова (рис. 3). Далее геометрические модели кости и фиксаторов импортировались в систему конечно-элементного моделирования, где и проводились расчеты напряженно-деформированного состояния систем.
Рис. 3. Модель фиксации перелома на примере блокируемого интрамедуллярного стержня (в разрезе продольной плоскостью).
Задача о нагружении системы «большеберцовая кость - фиксатор» решалась в трехмерной постановке. В программном комплексе ANSYS решались стационарные уравнения теории упругости в перемещениях Навье-Ламе. Нижний торец кости жестко закреплялся, а к верхнему концу (суставная площадка проксимального эпифиза большеберцовой кости) прикладывалась осевая нагрузка 450 Н (рис. 4, 5).
Рис. 4. Общая картина деформированного состояния модели «кость - фиксатор» (слева дана цветовая шкала величин перемещений в м).
Рис. 5. Общая картина деформированного состояния модели «кость - АНФ» (слева дана цветовая шкала величин перемещений в м).
Картина деформированного состояния большеберцовой кости, блокируемого интрамедуллярного стержня и аппарата Г.А. Илизарова под действием нагрузки 450 Н позволяют сделать вывод об адекватности представленных моделей «кость - фиксатор».
Между элементами фиксатора (аппарат Илизарова, интрамедуллярный гвоздь) и костными отломками создавались условия контакта с трением. Коэффициент трения составлял 0,55. Механические свойства фиксаторов задавались на основе технических данных изготовителя, а механические свойства кости на основе литературных данных. Материал кости и фиксаторов предполагался изотропным идеально-упругим.
Таким образом, алгоритм решения задачи методом конечных элементов состоял из следующих шагов.
- Выбор математической модели. Это означает, что необходимо выбрать соответствующий конечный элемент или задать дифференциальные уравнения, описывающие исследуемый процесс.
- Создание или импорт геометрической модели.
- Задание свойств материалов, граничных и начальных условий, параметров решателя.
- Разбиение модели на конечно-элементную сетку.
- Решение задачи.
- Визуализация и обработка полученных результатов.
Оценка с помощью конечно-элементной программы ANSYS напряжений в системе «кость - фиксатор» в условиях статического блокирования под действием осевой нагрузки в 450 Н выявила максимальное напряжение в 7,44 МПа (рис. 6).
Рис. 6. Общая картина напряженного состояния модели (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
В результате оценки максимальных напряжений в стержне и блокирующих элементах с помощью метода конечных элементов в условиях статического блокирования нестабильного перелома большеберцовой кости было выявлено значение 24,76 МПа (рис. 7).
Рис. 7. Общая картина напряженного состояния фиксатора (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
Максимальное напряжение, приходящееся на проксимальный отломок большеберцовой кости, равно 3,9 МПа (рис. 8).
Рис. 8. Напряженное состояние проксимального отломка модели (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
Максимальное напряжение в центральном отломке равно 4,04 МПа (рис. 9, 10).
Рис. 9. Напряженное состояние центрального отломка модели (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
Рис. 10. Напряженное состояние центрального отломка модели (повернуто).
Максимальное напряжение в дистальном отломке равно 4,1 МПа (рис. 11).
Рис. 11. Напряженное состояние дистального отломка модели (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
При оценке по методу конечных элементов максимального эквивалентного напряжения, которое возникает в стержне толщиной 10,0 мм, длиной 285 мм, а также верхних и нижних стальных блокирующих винтах диаметром 4,5 мм, в условиях статического блокирования под действием осевой нагрузки в 450 Н, максимальные эквивалентные напряжения не превышают 24,76 МПа, что делает запас прочности достаточным.
Оценка с помощью конечно-элементной программы ANSYS напряжений в системе «кость - АНФ» в условиях чрескостного остеосинтеза аппаратом Г.А. Илизарова под действием осевой нагрузки в 450 Н равна 4,52 МПа (рис. 12).
Рис. 12. Общая картина напряженного состояния модели (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
Максимальное напряжение в системе «кость - АНФ» испытывает центральный отломок в месте проведения репозиционной спицы с напайкой - 39 МПа (рис. 13).
Рис. 13. Общая картина напряженного состояния модели (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
Максимальное напряжение в проксимальном отломке равно 5,58 МПа (рис. 14).
Рис. 14. Напряженное состояние проксимального отломка модели (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
Максимальное напряжение в центральном отломке равно 7,61 МПа (рис. 15).
Рис. 15. Напряженное состояние центрального отломка модели (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
Максимальное напряжение в дистальном отломке равно 3,97 МПа (рис. 16).
Рис. 16. Напряженное состояние дистального отломка модели (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
С помощью метода конечных элементов оценивалось максимальное эквивалентное напряжение, которое возникает в кольцах, с внутренним диаметром 160 мм, толщиной 5 мм и шириной 20 мм, имеющих отверстия диаметром 7 мм через каждые 15°. Предложенная компоновка аппарата под действием осевой нагрузки в 450 Н испытывает напряжение, равное 10 МПа. В спицах диаметром 1,5 мм, в том числе с напайками, максимальное эквивалентное напряжение было равно 31 МПа, что также делает запас прочности фиксирующей системы достаточным (рис. 17).
Рис. 17. Общая картина напряженного состояния фиксатора (слева дана цветовая шкала величин эквивалентных напряжений по von Mises в Па).
Сравнение результатов компьютерного моделирования показало, что в системе «кость - блокируемый стержень» механическое напряжение в фиксаторе достигает максимальной величины 24,76 МПа (табл. 1), что меньше, чем в системе «кость - АНФ» - 31 МПа (табл. 2). Длительная концентрация напряжений на элементах фиксирующего устройства может послужить причиной разрушения последних.
Таблица 1 - Сводные результаты конечноэлементного моделирования по системе «кость - блокируемый стержень»
Локализация |
Напряжение, МПа |
Максимальное напряжение, МПа |
Система «кость - блокируемый стержень» |
7,44 |
7,44 |
Проксимальный отломок |
3,9 |
4,04 |
Центральный отломок |
4,04 |
|
Дистальный отломок |
4,1 |
|
Фиксатор |
24,76 |
24,76 |
Таблица 2 - Сводные результаты конечноэлементного моделирования по системе «кость - АНФ»
Локализация |
Напряжение, МПа |
Максимальное напряжение, МПа |
Система «кость - АНФ» |
4,52 |
4,52 |
Проксимальный отломок |
5,58 |
7,61 |
Центральный отломок |
7,61 |
|
Дистальный отломок |
3,97 |
|
Кольца |
10 |
31 |
Спицы |
31 |
|
Напряжение репозиционной спицы |
39 |
39 |
Как в системе «кость - блокируемый стержень», так и в системе «кость - АНФ» имеет место концентрация напряжений, приходящихся на кость. Так, максимальное механическое напряжение и в первом, и во втором случае имеет место в центральном отломке, однако в условиях остеосинтеза аппаратом Г.А. Илизарова концентрация напряжений выражена в большей степени. В результате компьютерного моделирования было доказано, что система «кость - блокируемый стержень» имеет преимущество, так как концентрация напряжений в отломках в этом случае меньше, чем в системе «кость - АНФ», что особенно важно в условиях неопорного перелома, так как оно может вести к смещению отломков. Максимальное напряжение в системе «кость - АНФ» имеет место в области проведения репозиционной спицы, что в условиях избыточной нагрузки может сопровождаться переломом последней и смещением отломка.
Заключение. Реализация методологии математического моделирования элементов опорно-двигательного аппарата человека как сборной конструкции на основе метода конечных элементов с решением контактной задачи теории упругости является актуальным научным направлением. Эта методология позволяет решать на более высоком уровне информативности задачи определения тактики лечения и выбора оптимального метода фиксации, а также прогнозирования хода лечебных мероприятий в области травматологии.
Рецензенты:
- Мельцер Рудольф Ионтелевич, д.м.н., профессор, заведующий курсом травматологии, ортопедии и ВПХ кафедры общей и факультетской хирургии ПетрГУ, г. Петрозаводск.
- Колесников Геннадий Николаевич, д.т.н., профессор, заведующий кафедрой механики ПетрГУ, г. Петрозаводск.