Адаптационные процессы в кости после остеосинтеза шейки бедра упругими резьбовыми фиксаторами
Автор: Акулич Ю.В., Акулич А.Ю., Денисов А.С., Подгаец Р.М.
Журнал: Российский журнал биомеханики @journal-biomech
Статья в выпуске: 3 (37) т.11, 2007 года.
Бесплатный доступ
Представлены математическая модель и расчётная схема решения начально-краевой задачи адаптивной пороупругости в её одномерной постановке для области шейки бедра, подвергнутой осевому сжатию упругими резьбовыми имплантатами в результате хирургической операции остеосинтеза по поводу перелома. В качестве определяющего соотношения среды принимается феноменологическое уравнение, полученное в работе (Акулич и др., Российский журнал биомеханики, 2005, №2). Губчатая костная ткань рассматриваемой области бедра считается неоднородной и изотропной, а структурная модель включает сплошной матрикс и пустые сферические поры в нём. Связь скорости изменения радиуса пор с деформационным стимулом адаптации устанавливается с помощью кинетического уравнения мезоуровня (Сотин и др., Российский журнал биомеханики, 2001, №5). Процессы в межфрагментарном пространстве не рассматриваются, а в качестве необходимых условий сращения принимается наличие сжимающего усилия, обеспечивающего постоянный контакт поверхностей перелома и обычную (не патологическую) пористость костной ткани фрагментов. Исследуется влияние жёсткости и удлинения установки пружин упругих фиксаторов на послеоперационные временные изменения параметров структуры и механических свойств костной ткани.
Шейка бедра, остеосинтез, упругие резьбовые фиксаторы, математическое моделирование, адаптация, структурные и механические свойства кости
Короткий адрес: https://sciup.org/146215898
IDR: 146215898
Текст научной статьи Адаптационные процессы в кости после остеосинтеза шейки бедра упругими резьбовыми фиксаторами
В структуре переломов у лиц пожилого и старческого возраста повреждения шейки бедра составляют около 10%. Общепринятый способ лечения переломов (остеосинтез) заключается в близком сопоставлении отломков (репозиция) и их надёжной фиксации жесткими резьбовыми фиксаторами (рис. 1а).
Однако результаты остеосинтеза переломов шейки бедра с использованием жестких резьбовых фиксаторов нельзя признать удовлетворительными.
Несращения зафиксированы у 15% пациентов с переломами шейки бедра, подвергшихся остеосинтезу жёсткими резьбовыми фиксаторами. Причиной неудачных исходов, вероятнее всего, является ослабление фиксации отломков шейки бедра до их

Рис. 1. Жёсткий (а) и упругий (б) резьбовые фиксаторы сращения, поскольку в этот период пациент передвигается на костылях без опоры возникают микроперемещения, препятствующие сращению на оперированную конечность и при ослабленной фиксации между отломками.
Так как фиксация жёсткими фиксаторами осуществляется благодаря упругому деформированию костной ткани под витками резьбы фиксаторов и деформации при этом малы (по нашим расчётам не превышают 3%), то в случае возникновения резорбции в этой области кости уровень упругих свойств костной ткани быстро снижается (модуль упругости уменьшается от 0,61 ГПа до 0,15 ГПа за 35 дней), что ведёт к существенному ослаблению фиксации.
Применение упругого имплантата (рис. 1б) позволяет распределить потенциальную энергию деформирования системы «костная ткань - фиксаторы» между костной тканью и пружинами фиксаторов. Предполагается, что рост деформации сжатия костной ткани под витками резьбы, вызванный её резорбцией, приводит к существенно более малому уменьшению потенциальной энергии пружины и, следовательно, к более медленному уменьшению величины силы сжатия отломков, то есть к более стабильной фиксации. Для обоснования этого предположения необходимо установление характера и количественных оценок послеоперационных изменений механических свойств и структурных параметров костной ткани системы «костная ткань - упругие фиксаторы».
В данной статье с помощью одномерной математической модели решается начально-краевая задача адаптивной пороупругости для проксимального отдела бедра, подвергнутого остеосинтезу тремя упругими фиксаторами (рис. 2). При этом делаются следующие допущения:
-
• процессы в межфрагментарном пространстве не рассматриваются;
-
• поверхность перелома резорбируется с постоянной скоростью;

Рис. 2. Проксимальный отдел бедра с упругими фиксаторами. 1 - ось шейки, 2 - головка, 3 - предполагаемая линия перелома, 4 - ось фиксатора, 5 - головка фиксатора
-
• активация мышц тазобедренной группы при ходьбе на костылях за время сращения равномерно уменьшается на 30%;
-
• влияние тонкого слоя кортикальной костной ткани проксимального отдела бедра на его напряженно-деформированное состояние не учитывается;
-
• губчатая костная ткань изотропна.
Структура губчатой костной ткани цилиндра моделируется сферическими пустыми порами, равномерно распределёнными в сплошном матриксе. Параметрами структуры являются радиус r и плотность пор n , имеющие различные значения в разных участках цилиндра. Изотропия губчатой костной ткани в шейке в радиальном направлении от её оси подтверждается круговой формой проекции тензора структуры на плоскость, перпендикулярную оси шейки бедра [15].
Математическая модель
Моделируются этап установки фиксаторов и процесс адаптационных изменений, в течение которого происходит сращение перелома без нагрузки на конечность (около четырёх месяцев ходьбы на костылях). Исследуются временные изменения механических свойств (модуль упругости, деформация, усилие сжатия фрагментов) и параметров структуры (радиус и плотность пор и пористость) при различных значениях жёсткости и удлинения установки пружин фиксаторов.
При моделировании используются следующие допущения:
-
1) до начала установки фиксаторов деформации и напряжения в костной ткани отломков равны нулю;
-
2) давление резьбы на кость однородно по длине резьбовой части фиксатора;
-
3) напряженно-деформированное состояние области остеосинтеза принимается одномерным (вдоль оси шейки).
Одномерное представление области остеосинтеза
Последнее допущение позволяет представить рассматриваемую область составным круговым цилиндром, вычленяемым из проксимального отдела бедра (рис. 3). Диаметр цилиндра равен наименьшему диаметру шейки, а длина распространяется от подвертельной области до купола головки, не выходя за его пределы. Анализ экспериментальных данных о распределении радиуса пор [11] позволяет выделить в цилиндре пять характерных однородных участков по его длине, показанных на рис. 3. Шестой участок включает костную ткань трёх спиралевидных лент, расположенных под витками резьбы трёх фиксаторов, и условно представляется цилиндром высотой, равной толщине ленты. Структура и механические свойства костной ткани этого участка до начала адаптации такие же, как у костной ткани пятого участка, в котором расположена резьба фиксатора. Таким образом, вычленяемый цилиндр условно достраивается дополнительным цилиндрическим участком 6. Участки считаются не связанными между собой, то есть представляют отдельные тела, контактирующие без трения.
Для определения начальных для процесса адаптации значений напряжений и деформаций ниже рассматривается равновесие системы «участки костной ткани - фиксаторы» в момент окончания затяжки фиксаторов.
Ось шейки
Линия перелома

Оси
Вычленяемый цилиндр
Рис. 3. Аппроксимация области шейки круговым цилиндром, состоящим из шести однородных участков: 1 - треугольник Варда и подвертельная область, 2 - срединная часть шейки, 3 - область
фиксаторов
сращения, 4 - субкапитальная часть головки, 5 - центральная область головки, 6 - контактный слой костной ткани с фиксатором
Установка фиксаторов
Уравнения равновесия для каждого участка без учёта объёмных сил имеют вид
^ - 0, к - 1,6. 5 x
Следовательно, присутствует однородность напряжений и деформаций по длине участков.
Напряжения в первых пяти участках равны, отличаются от напряжения под резьбой фиксатора и определяются соотношением
( 3 F 0 + К ,0 ) / s , , к - 1,5,
F 0 / sb + Rx 0 / sn , k = 6, где sb ,sn – площади спиралевидной ленты и поперечного сечения цилиндра, соответственно, Rx0 - проекция реакции тазовой кости на ось шейки в начале послеоперационного периода ( Rx0= 90 Н [12]), F0 - сила сжатия шейки одним фиксатором, определяемая законом Гука
F 0 - c A l 0 ,
где A 1 0 - удлинение пружины при установке фиксатора.
Деформации участков определяются в соответствии с определяющим соотношением адаптивной пороупругости для напряжения [2]
8 к 0 - а к 0 / ^ к 0 Ek 0, k - 1,6, где £ко, Еко - объёмное содержание костного вещества (матрикса) и модуль упругости костной ткани k-го участка, соответственно. Для сферических пор и принятой выше структуры костной ткани объёмное содержание матрикса выражается через параметры структуры (радиус пор r и объёмная плотность пор n) по формуле [10]
5 к = 1 - 4 / 3п ntr^,k = 16. (5)
Аппроксимация экспериментальных данных с применением блочной модели структуры губчатой костной ткани [4] показала, что модуль упругости костной ткани аппроксимируется кубической параболой вида
E k 0 = ^ k 0 E m , (6)
где E m - модуль упругости матрикса.
Формулы (1) - (6) используются также и для определения значений параметров структуры и механических характеристик костной ткани в равновесном состоянии при F= 0, то есть в состоянии до травмы при физиологической нагрузке. При весе тела равном 720 Н Rx =1220 Н [7]. В табл. 1 приведены значения параметров структуры, деформации и модуля упругости в каждом участке цилиндра для равновесного состояния, рассчитанные по заданному распределению величины радиуса пор вдоль оси шейки.
Таблица 1
Значения параметров структуры, деформации и модуля упругости участков в равновесном состоянии
Величина |
Номер участка |
|||||
1 |
2 |
3 |
4 |
5 |
6 |
|
Радиус пор, мм |
0,287 |
0,303 |
0,297 |
0,307 |
0,286 |
0,286 |
Плотность кости 1, г/см3 |
1,36 |
1,27 |
1,3 |
1,32 |
1,34 |
1,34 |
Плотность пор, 1/мм3 |
6,8 |
6,2 |
6,3 |
6,5 |
6,6 |
6,6 |
Объёмное содержание матрикса |
0,364 |
0,271 |
0,293 |
0,319 |
0,334 |
0,334 |
Деформация сжатия ε eq |
0,0042 |
0,0135 |
0,01 |
0,007 |
0,0058 |
0,0058 |
Модуль упругости2, ГПа |
0,79 |
0,33 |
0,42 |
0,54 |
0,61 |
0,61 |
Данные таблицы показывают, что у конкретного пациента наибольшей плотностью и жесткостью обладает губчатая костная ткань подвертельной области (участок 1), затем - головки (участки 5 и 6), а шейка (участки 2 и 3) имеет наименьшую плотность и жёсткость, и, очевидно, наименьшую прочность.
Исследуется установка четырёх типов упругих фиксаторов, различаемых величиной жёсткости пружины: 30, 35, 45 и 58 кН/м, при двух значениях удлинения пружины ∆l0 =1,5 и 2,0 мм. Выбор вариантов обусловлен требуемой прочностью пружины фиксатора.
В табл. 2 приводятся деформации сжатия участков цилиндра ε k 0 , k = 1 , 6 для рассматриваемых вариантов установки фиксаторов. Данные таблицы показывают, что во всех вариантах установки фиксаторов наибольшему сжатию подвергается костная ткань под их резьбой (участок 6), поскольку площадь опирания резьбы на кость меньше площади поперечного сечения шейки в семь раз, тогда как усилие меньше только в три раза. С ростом жёсткости и удлинения пружины деформация сжатия во всех областях шейки возрастает.
Важным этапом характеристики режима установки фиксаторов является сопоставление величин деформаций от сжатия фиксаторами (табл. 2) с уровнем деформаций равновесного состояния ε eq (табл. 1). Из диаграммы, изображенной на рис. 4, следует, что при установке всех четырёх типов фиксаторов с удлинением пружины на 1,5 мм, деформация сжатия костной ткани только в шестом участке близка к равновесной, а в остальных участках она значительно ниже этого уровня, то есть костная ткань недогружена.
При удлинении пружины на 2,0 мм (рис. 5) деформация костной ткани в шестом участке превышает равновесную (перегрузка), если применяются фиксаторы с жёсткостью пружины 58 кН/м. В остальных участках наблюдается недогрузка при установке фиксаторов всех типов.
Адаптивное изменение механических свойств
Предполагается, что после остеосинтеза в костной ткани одновременно осуществляются три главных процесса: адаптационной перестройки, предварительной резорбции поверхностей перелома и сращения отломков. С позиции биомеханики важно отметить, что если процесс адаптации охватывает все участки костной ткани, нагруженные фиксаторами, то процессы предварительной резорбции и сращения являются локальными и распространяются в пределах малой области (участок 3 имеет небольшую длину вдоль оси шейки, равную 4 мм).
Таблица 2
Распределение деформации сжатия вдоль оси шейки после установки фиксаторов
с, кН/м |
∆ l 0 ⋅ 103, м |
Номер участка |
|||||
1 |
2 |
3 |
4 |
5 |
6 |
||
30 |
1,5 |
0,0008 |
0,0025 |
0,0026 |
0,0013 |
0,0011 |
0,0046 |
2,0 |
0,0009 |
0,0030 |
0,0031 |
0,0016 |
0,0013 |
0,0051 |
|
35 |
1,5 |
0,0008 |
0,0027 |
0,0028 |
0,0014 |
0,0012 |
0,0049 |
2,0 |
0,0010 |
0,0033 |
0,0035 |
0,0017 |
0,0014 |
0,0055 |
|
45 |
1,5 |
0,0010 |
0,0032 |
0,0034 |
0,0017 |
0,0014 |
0,0054 |
2,0 |
0,0012 |
0,0040 |
0,0042 |
0,0021 |
0,0017 |
0,0062 |
|
58 |
1,5 |
0,0012 |
0,0039 |
0,0040 |
0,0020 |
0,0017 |
0,0061 |
2,0 |
0,0015 |
0,0048 |
0,0051 |
0,0025 |
0,0021 |
0,0071 |

Рис. 4. Деформации участков костной ткани от сжатия фиксаторами различной жёсткости (табл. 2) с уровнем деформаций равновесного состояния ε eq (табл. 1) при удлинении пружины на 1,5 мм

Рис. 5. Деформации участков костной ткани от сжатия фиксаторами различной жёсткости (табл. 2)
с уровнем деформаций равновесного состояния ε eq (табл. 1) при удлинении пружины на 2,0 мм
На основании выше было принято допущение о пренебрежимо малом влиянии сращения на адаптацию (однако при исследовании процесса сращения необходимо учитывать адаптацию). В то же время необходимо принять во внимание процесс предварительной резорбции, поскольку он вызывает укорочение шейки и, следовательно, влияет на величину силы сжатия отломков, то есть действует на область остеосинтеза. Учёт этого процесса осуществляется путём укорочения длины третьего участка в течение 7 - 10 дней [13] с постоянной скоростью резорбции 1 мкм/день.
Длительность τ процесса адаптации принимается равной суммарному времени процессов предварительной резорбции и сращения, то есть 125 дням. В этот период пациент передвигается на костылях, не опираясь на больную конечность, поэтому мышцы тазобедренной группы активируются только весом конечности (120 Н) и величина силы Rx существенно уменьшается по сравнению с её равновесным физиологическим значением, составляя в начале послеоперационного периода 90 Н [12] с уменьшением из-за недогрузки к его окончанию. Функция R x = R x ( t ) гипотетически предполагается линейной со скоростью убывания 0,25 Н/день, определённой из расчёта 30%-го уменьшения к концу четырех месяцев ходьбы на костылях.
В текущий момент времени процесса адаптации напряжения в участках определяются с помощью соотношений, имеющих форму (2), поэтому скорости изменения напряжений имеют вид
( 3 F + R , ) / s , k = 1,5, F / S b + R x I S n , k = 6.
Для каждого участка определяющие соотношения имеют вид [10]
Gk ( в k в k ) Rk в k + ^ kEk в k,k 1,6, где чувствительность Rk костной ткани к стимулу адаптации вk - вke определяется соотношением, полученным с учётом зависимости E = £3 Em [4],
Rk = 16пankrk Ek , где a - коэффициент, характеризующий активность и плотность костных клеток на поверхности пор [8], равный 0,055 мм/(с-ЕОД) (ЕОД - единиц относительной деформации).
Адаптивное изменение параметров структуры
Определяющее уравнение структурной перестройки имеет вид [10]
r k = - a (в k -в keq ) , k = 1,6. (10)
Предполагается, что функция rк = rк (вk) (10) является кусочно-линейной (рис. 6) и имеет три характерных точки: вcr - критическое значение деформации сжатия костной ткани, соответствующее предельному напряжению сжатия gcr = 7 МПа [14], превышение которого вызывает резорбцию; вeq - в1 , вeq + в2 - нижняя и верхняя границы зоны нечувствительности к изменению деформации [1], соответственно.
Параметры структуры (радиус пор и плотность пор) не являются независимыми, поскольку в любой момент времени процесса они связаны какой-либо физической величиной, например, плотностью сырой костной ткани ρск . В этом случае, как показано в работе [11], плотность пор определяется соотношением n (r ) =
v
РсЛr)-P^^ 3
ск ж
Р м -Р ж ) 4 п r3
где р ж , р м - плотности внутритканевой жидкости и матрикса, соответственно, а р ск ( r ) = 2 , 395 - 3 , 71 r, r е [0,23 - 0,35 мм]. (12)
График зависимости (11) (рис. 7) показывает, что в указанном интервале изменения радиуса пор плотность пор уменьшается примерно в два раза. Следует заметить, что при этом плотность сырой костной ткани находится на уровне, установленном экспериментально [11].
Изменение плотности пор при изменении их радиуса говорит о возможности описания соотношениями (10) и (11) пересечения пор в случае резорбции, то есть образование более крупных пор, что наблюдается в гистологических исследованиях, например, при остеопорозе [5].

Рис. 6. Вид функции Г- = r ( s ) . При £ < £ eq +£ 2 наблюдается резорбция костной ткани, в интервале
ε eq – ε 1 ≤ ε ≤ ε eq +ε 2 адаптация отсутствует (деформации ε 1 , ε 2 , определяющие зону нечувствительности, принимаются равными 0,005), в интервале ε cr ≤ ε ≤ ε e q –ε 1 осуществляется адаптация с производством матрикса, при £ < £ cr - резорбция, tg а = a .

0,20 0,24
0,28 0,32 0,36 0,40
Радиус пор, мм
Рис. 7. Изменение плотности пор с ростом (или убыванием) их радиуса в процессе адаптации в соответствии с зависимостями (11) и (12)
Геометрическая связь
Одномерную модель системы «костная ткань-фиксаторы» удобно представить в виде совокупности упругих и жёстких элементов (рис. 8). Участки костной ткани обозначены последовательно соединёнными упругими элементами жёсткостью Ck,k = 1,6 . Фиксаторы изображены в виде двух недеформируемых стержней, соединённых пружиной жёсткости 3с. Вертикальные плиты жёстко соединены со стержнями и внешними торцами первого и шестого участков и могут свободно скользить в горизонтальных направляющих, объединяя все упругие элементы в систему. Механическая модель показывает, что фиксаторы являются упругой связью c1
3 c 6

R x
ЛААЛ —---
Плита Костная ткань Фиксаторы
x
VZZZZZZZZZZA
Рис. 8. Механическая модель системы «костная ткань - упругие фиксаторы»:
c k - жёсткости участков костной ткани, где k =1..6; Rx - проекция реакции тазовой кости на ось шейки
для костной ткани и, следовательно, необходимо установить её аналитическое выражение. В каждый момент времени процесса адаптации рассматривается упругое равновесие между костной тканью и фиксаторами, поэтому связь является геометрической и имеет вид (см. рис. 9)
Д 1ш = (Д 1ш )о + Д10 -Д1,
где Д 1ш, Д l - абсолютные деформации костной ткани и пружин в текущий момент времени процесса адаптации, соответственно, ( Д 1ш ) 0 , Дl 0 - значения этих величин в начале процесса адаптации, то есть после установки фиксаторов. Производная по времени от соотношения (13) равна
d Д 1ш /dt = - d Д l / dt .
Поскольку
Д 1ш ХДlk X ё k -(lk )о , к=1 к=1
где Дl k - абсолютная деформация к -го участка костной ткани в текущий момент времени процесса адаптации, (lk ) о - длина к -го участка костной ткани до установки фиксаторов, и выполняется закон Гука для пружины фиксатора
F = c -Д l , (16)
то из равенства (14) следует соотношение дифференциальной интегрируемой связи
X ё к -( lk ) 0 =- F Iе . к = 1
c 2
c 4
3 с 6
(Д l ш ) о
I
II
III
iwwm/HMWNwi i
x
3 с
Д l
Д l о
ш
Рис. 9. К установлению геометрической связи (левая плита механической модели зафиксирована):
I - костная ткань до установки фиксаторов, II - равновесие после установки фиксаторов, III - равновесие в текущий момент времени адаптации
Алгоритм решения начально-краевой задачи
В уравнениях (7), (8) и (17) имеется 13 неизвестных функций: F, гk, г k,k = 1 , 6 , которые вычисляются с помощью пошаговой по времени процедуры. В качестве начальных значений деформаций принимаются величины, определяемые соотношением (4). На каждом временном шаге функции F, г k , к = 1 , 6 считаются постоянными, и формируется система пяти линейных алгебраических уравнений относительно г k , к = 1 , 5, решаемая методом Гаусса.
Результаты моделирования
Для тех же значений жёсткости и начального удлинения пружины фиксатора, с которыми осуществлялся остеосинтез (см. раздел 2.2), рассчитывалось изменение параметров структуры и механические свойства костной ткани в послеоперационном периоде. Жёсткость варьировалась на четырёх уровнях (30, 35, 45 и 58 кН/м), а начальное удлинение пружины фиксатора на двух (2,0 и 1,5 мм), то есть исследовались последствия восьми вариантов остеосинтеза. Наименьшее усилие сжатия достигается в варианте остеосинтеза при с= 58 кН/м и A l q = 2 , 0 мм, а наименьшее - при с= 30 кН/м и A l о = 1 , 5 мм.
На рис. 10 представлены процессы изменения пористости при наименьшем (А) и наибольшем (Б) усилиях сжатия отломков при остеосинтезе. Значение пористости 0,8 принимается в данной работе предельным (пористость нормальной, без патологии, губчатой костной ткани лежит в интервале 0,5 - 0,75 [3]) и обозначается на рис. 10. пунктирной линией. Сравнение кривых показывают, что в обоих случаях установки фиксаторов в участках 1 - 5 осуществляется резорбция кости в виду её недогрузки, что вызывает рост пористости. Различие этих процессов заключается в их интенсивности.

А Время, день

Б Время, день
Рис. 10. Изменение пористости участков шейки при установке фиксаторов с различной жёсткостью с и начальным удлинением A 1 0 пружины: А - с = 30 кН/м, A 1 0 =1,5 мм; Б - с = 58 кН/м, A 1 0 =2,0 мм.
Цифры у кривых обозначают номер участка костной ткани, пунктирная линия - верхнее предельное значение пористости в норме
При остеосинтезе с наименьшим усилием сжатия (рис. 10А) интенсивность резорбции выше ввиду большей недогрузки. Поэтому при данном варианте остеосинтеза предельный уровень пористости достигается на втором участке (срединная часть шейки) уже через 35 дней, на третьем (область сращения) - через 2,5 месяца. Предельная пористость костной ткани в области перелома, очевидно, не способствует процессу сращения, поэтому данный вариант остеосинтеза считается неуспешным. С увеличением усилия сжатия костной ткани при остеосинтезе интенсивность резорбции уменьшается и при наибольшем значении усилия сжатия предельный уровень пористости достигается только на втором участке, что позволяет считать этот вариант остеосинтеза более предпочтительным.
В вариантах остеосинтеза с промежуточными величинами усилий сжатия характер изменения пористости остаётся таким же, как на рис. 10А, с уменьшением наклона кривых по мере роста усилия сжатия отломков при установке фиксаторов.
В слое костной ткани под резьбой фиксаторов (участок 6) во всех вариантах остеосинтеза, кроме варианта с наибольшим усилием, пористость не изменяется, так как деформация сжатия этого участка остаётся в области нечувствительности ε eq - ε1 ≤ ε ≤ ε eq + ε2 , где адаптация не стимулируется (см. рис. 6). В варианте остеосинтеза с наибольшим усилием сжатия деформация этого участка значительно превосходит уровень равновесной деформации (рис. 5), выходит из области нечувствительности и попадает в интервал ε cr ≤ ε ≤ ε eq - ε1 (см. рис. 6), где осуществляется адаптация с производством матрикса. Поэтому пористость костной ткани этого участка уменьшается (рис. 10Б).
Процесс изменения модуля упругости показан на рис. 11. Характер изменения приведенных кривых противоположен соответствующим кривым изменения пористости (рис. 10). Этот факт объясняется тем, что поскольку пористость p и объёмное содержание матрикса ξ m костной ткани с пустыми порами связаны соотношением [11]

Рис. 11. Изменение модуля упругости костной ткани после остеосинтеза с наименьшим (А) и наибольшим (Б) усилиями сжатия отломков

p = 1 - ξ m , (18)
то производные p и ^ m имеют разные знаки и в силу соотношения (6) росту пористости соответствует уменьшение модуля упругости. Сравнение кривых 1 - 5, изображённых на рис. 11А и на рис. 11Б, показывает, что с увеличением усилия сжатия костной ткани при остеосинтезе замедляется процесс снижения упругих свойств костной ткани участков 1 - 5. Из сравнения кривых 6 следует, что при остеосинтезе с наибольшим усилием сжатия (рис. 11Б) жёсткость костной ткани под резьбой фиксатора (участок 6) даже возрастает.
На рис. 12 представлены кривые послеоперационного изменения силы осевого сжатия шейки при установке фиксаторов с различной жёсткостью с и начальным удлинением ∆ l 0 пружины. Кривые показывают, что при всех вариантах остеосинтеза происходит уменьшение сжимающей силы в послеоперационный период. Скорость и величина её уменьшения к концу периода сращения практически одинаковы и составляют в среднем 0,4 Н/день и 22%, соответственно.
Заключение
Полученные результаты моделирования послеоперационного изменения биомеханического состояния костной ткани шейки бедра, подвергнутой остеосинтезу тремя упругими резьбовыми фиксаторами, показывают, что представленная модель качественно правильно описывает как этап установки фиксаторов, так и последующий процесс адаптационных изменений структуры и механических свойств костной ткани.
Применение упругих резьбовых фиксаторов в исследованном интервале изменения их жёсткости и величины удлинения при закручивании обеспечивает осевое сжатие шейки силой не более 440 Н при установке фиксаторов и не менее 360 Н к окончанию периода сращения. Поскольку на протяжении всего периода сращения уменьшение силы сжатия происходит монотонно и достаточно медленно, то достигнутую с помощью упругих фиксаторов фиксацию можно считать стабильной.

Рис. 12. Изменение силы осевого сжатия шейки в зависимости от жёсткости пружины фиксатора с и величины её удлинения при установке ∆ l 0 . Сплошной линией обозначены процессы при ∆ l 0 =2,0 мм, пунктиром - при ∆ l 0=1,5 мм
Однако для выбора необходимой жёсткости пружины фиксатора и предельной величины её удлинения при закручивании имплантата требуется оценка наименьшего допустимого уровня усилия сжатия отломков кости, удовлетворяющего условиям сращения губчатой костной ткани при одноопорной ходьбе пациента на костылях.
Благодарности
Работа выполнена при финансовой поддержке РФФИ, проект № 07-01-00787.