CC BY
76
© Группа авторов, 2008
Механизм жесткости1 и прочности2 кости в норме и при старении организма. Наноуровневая модель
А.С. Аврунин*, Р.М. Тихилов*, Л.К. Паршин**, И.И. Шубняков*, Б.Е. Мельников**
The mechanism of bone rigidity1 and strength2 in normal conditions
and at ageing. A nanolevel model
A.S. Avrunin*, R.M. Tikhilov*, L.K. Parshin**, I.I. Shubniakov*, B.E. Melnikov**
*ФГУ «Российский ордена Трудового Красного Знамени научно-исследовательский институт травматологии и ортопедии им. Р.Р. Вредена Росмедтехнологий», Россия, Санкт-Петербург (директор - д. м. н., профессор Р.М. Тихилов) **Кафедра сопротивления материалов СПб государственного политехнического университета (заведующий кафедрой - д.т.н., профессор Б.Е. Мельников), Россия, Санкт-Петербург
Цель: разработать концепцию наноуровневой модели механизма, определяющего, с одной стороны, жесткость и прочность кости на макроуровне, а с другой - увеличение хрупкости внеклеточного матрикса при старении. В основу теории положено представление о том, что каждый кристаллит жёстко связан со структурированной водой, обеспечивающей их соединение между собой. Предложенная трехфазная модель элементарного локуса костной ткани на наноуровне удовлетворительно описывает связь между модулем упругости гидроксиапатита и модулем упругости костной ткани на макроуровне. Модель обеспечивает понимание причин разрушения костной ткани на молекулярном уровне путём среза по гидратному слою. Рассматриваемая теория объясняет необходимость существования винтовых (спиральных) структур в костях оптимизацией функционирования конвекционного механизма движения жидкости в канальцевой системе и дает теоретические предпосылки создания новых конструкционных материалов с особенностями строения на наноуровне, например, искусственной кости. Ключевые слова: механизм прочности кости, механизм жесткости кости, старение.
Aim: To develop the conception of the mechanism of nanolevel model, that determines, on the one hand, the rigidity and strength of bone at macrolevel, and on the other hand - the increasing fragility of extracellular matrix during ageing. This theory is based on the idea that every crystallite is rigidly bound with the structurized water, which provides their connection with each other. The three-phase model of elementary bone tissue locus proposed satisfactorily describes the connection between hydroxyapatite modulus of elasticity and bone tissue modulus of elasticity at macrolevel. The model supports the insight into the causes of bone destruction at molecular level by cutting of the hydrated layer. The theory considered explains the necessity of existing screw (spiral) structures in bones by optimizing the functioning of the convectional mechanism of fluid flow in the system of tubules and gives the theoretical premises for creation of new constructional materials with some structural features at nanolevel, for example, for that of an artificial bone.
Keywords: mechanism of bone strength, mechanism of bone rigidity, ageing.
ВВЕДЕНИЕ
Старение организма сопровождается развитием остеопороза - «безмолвной болезни», протекающей асимптоматично до момента нарушения целостности костей [29]. Только в США по этой причине ежегодно происходит 700000 переломов позвоночника, по 250000 - шейки бедра и дис-тального отдела лучевой кости и 300000 - в других частях скелета. Здравоохранению США стоимость лечения последствий остеопороза обходится в 10 млрд. долларов (без учета стоимости длительного ухода на дому). В мире количество больных с переломами шейки бедра постоянно увеличивается и согласно расчетам возрастет с 1,7
Жесткость - способность материала сопротивляться изменению размеров и формы при механической нагрузке. Чем выше жесткость, тем меньше деформация при одной и той же нагрузке [5, 13].
2 Прочность - способность материала выдерживать механическую нагрузку без разрушения [5, 13].
млн. в 1990 г. до 6,3 млн. в 2050 г. При этом около 20 % пострадавших умирают в течение первого года [12]. Вышеизложенное свидетельствует о том, что остеопороз является социально-экономической проблемой высокой значимости.
Одним из элементов, обеспечивающих решение этой проблемы, является понимание первичных механизмов изменения механических характеристик внеклеточного матрикса при старении. Именно эти изменения приводят к последующим сдвигам архитектуры на более высоком уровне структурной организации скелета и к потере костной массы [3, 4]. До настоящего времени основной и господствующей моделью, объясняющей механизм обеспечения жесткости и прочности костной ткани, является гипотеза J.D. Currey, разработанная в отделе зоологии Оксфордского университета [23].
Концепция Сиггеу. По мнению этого ав-
тора, кость представляет собой двухфазный материал, аналогичный фиберглассу, где смола - аналог коллагена, а стекловолокно - апатита. J.D. Cur-rey подчеркивает, что двухфазные материалы могут функционировать эффективно при наличии устойчивой связи между волокнами и окружающим их веществом, добавляя при этом, что на момент создания модели характер связей между коллагеном и апатитом не ясен. Важным, с нашей точки зрения, является утверждение о том, что любая теория требует очень сильных связей между этими двумя компонентами, иначе модуль упругости кости не будет сколь либо существенно отличаться от его значений у фибрилл коллагена.
В последние десятилетия детализированы не только характер и пространственная организация связей между органическими и минеральными структурами, но и особенности их иерархической организации [1, 2, 7]. Эти данные, с одной стороны, позволяют говорить о несоответствии модели J.D. Currey сегодняшним представлениям, а с другой - разработать теорию, определяющую на современной основе механические свойства костной ткани в норме и при старении организма.
Цель: предложить общую концепцию базисной наноуровневой модели механизма, определяющего, с одной стороны, жесткость и прочность кости на макроуровне, а с другой - увеличение хрупкости внеклеточного матрикса при старении.
Теоретико-экспериментальная основа физико-механической модели. Выполнение поставленной цели (создание модели) регламентируется следующими условиями:
1) отношение модуля упругости Ект костной ткани к модулю упругости Еап гидроксиапатита должно соответствовать экспериментальным данным, представленным в литературе;
2) механизм, определяющий жёсткость кости, должен одновременно являться элементом механизма локального повреждения костных структур на наноуровне для случая, когда действующая сила превышает уровень допустимых физиологических нагрузок;
3) механизм, определяющий прочность костной ткани в норме, должен обеспечивать увеличение ее хрупкости при старении.
Результаты исследования модуля упругости. Для костной ткани, образующей компактное вещество у человека, модуль упругости Ект колеблется, по разным данным, в пределах 4-20 ГПа (в среднем 14 ГПа) [24] или 10-22 ГПа (в среднем 16 ГПа) [26], что практически одно и то же. У гидроксиапатита этот показатель значительно выше: от 40 до 90 ГПа (в среднем 65 ГПа) [26] или 114,2 ГПа [16], а у «зелёно-синего» апатита от 54 до 62 ГПа [14]. Таким образом, приведенные значения модулей упругости для каждого объекта являются величинами одного порядка и их отношение, если ориентироваться на данные A.C. Lawson и J.T. Czernuszka [26], составляет:
Ещ/Еап = 0,18-0,30 (в среднем 0,24). При старении это отношение меняется, так как происходит снижение модуля упругости Ект от его уровня 15,2 ГПа в 35 лет на 2,3 % за каждое последующее десятилетие [32]. Путем простых расчетов получаем, что к 75 годам отношение Ект/Еап составит 0,16-0,27 (в среднем 0,215).
Модуль упругости коллагена Ексл колеблется от 0,01 до 0,57 ГПа [26], то есть на два порядка ниже Еап. Это позволяет утверждать, что при одинаковых деформациях минерального и органического компонента в последнем напряжения составляют 1 % (и менее) от уровня в минеральном матриксе. Следовательно, коллаген даже теоретически не может обеспечить наблюдаемую жёсткость кости (не только в условиях сжатия, но и в условиях растяжения), которая определяется лишь минеральным компонентом с учетом особенностей структуры его иерархической организации.
Для дальнейших рассуждений выделим следующие основополагающие для разработки модели факты.
1. Кристаллиты гидроксиапатита имеют крайне малый размер, и их величины близки не только у разных организмов, но и при формировании in vitro: ширина - от 7,0 до 15,8 нм, длина - от 13,0 до 37,3 нм и толщина - от 1,5 до 4,0 нм [1, 2].
2. Кристаллиты соединяются между собой конец в конец или бок в бок, формируя копла-нарные структуры. При этом крайне важен тот факт, что полного слияния кристаллических решеток расположенных рядом кристаллитов не происходит. Между ними имеется водная пленка (гидратный слой) [1, 6, 11].
3. Гидратный слой жестко сцеплен с кристаллитами и представляет собой ультратонкий нанослой адсорбата на поверхности адсорбента, что обеспечивает весьма высокий уровень межмолекулярного взаимного притяжения [11, 18].
4. В результате межмолекулярного притяжения водная плёнка связывает кристаллиты между собой.
5. Все структуры начиная с надмолекулярного уровня и до уровня остеонов имеют иерархически построенную винтовую (спиральную) организацию [1, 6, 22], в соответствии с которой выстраиваются и элементы минерального мат-рикса.
6. При старении минерализация внеклеточного матрикса увеличивается на 12 % и более [19] за счет вытеснения водной [28, 30, 31] и органической [30, 31] составляющих. Доля ос-теонов, минерализованных на 75 % и выше, достигает 40 % и более от их общего числа, в то время как в молодом возрасте не превышает 10 %. Повышается также удельный вес центральных каналов и остеоцитарных лакун, заполненных минерализованной тканью [25].
7. Размеры копланарных объединений кристаллитов увеличиваются с возрастом (длина на
17-20 % и ширина на 5-9 % при неизменной толщине) [6, 7]. Подобные сдвиги вызваны соответствующими изменениями размеров собственно кристаллитов, их длины и ширины [21], причем размеры по А-оси меняются быстрее и сильнее, чем по С-оси [20].
8. В старших возрастных группах происходит увеличение содержания ионов карбоната (СО3--), замещающих в кристаллической решетке ионы фосфата (НРО4--), а также повышается количество фторидов и таких микроэлементов, как стронций, кремний и медь [9, 19].
9. Возникающие при старении замены в кристаллической решетке возможны только при условии возникновения аналогичных изменений в гидратном слое. Сложный процесс гетероион-ного обмена - основной фактор, определяющий состав минерального компонента кости [11].
В заключение теоретико-экспериментального обоснования модели следует указать, что она должна:
1) удовлетворительно описывать снижение модуля упругости Ект по сравнению с модулем упругости Еап;
2) учитывая наличие структурированной воды, дать теоретически обоснованную оценку ее модуля упругости Евод;
3) обеспечить понимание условий возникновения первичных микроповреждений;
4) указать причины уменьшения прочности и увеличения хрупкости костной ткани при старении организма;
5) дать теоретически обоснованную оценку снижения модуля упругости костной ткани Ект при старении.
Общие характеристики элементов физико-механической модели
Гидратный слой. Обозначим толщину этого слоя t (рис. 1). Каждый из двух расположенных рядом кристаллитов окружён своим гидратным слоем. Поэтому толщина разделяющей их водной плёнки составляет 2t или 1,3 нм, что соответствует 6-10 размерам молекул воды [11].
Рис. 1. Схематическая организация минерального мат-рикса на уровне отдельных кристаллитов с учетом структурированной воды. Обозначения: t - толщина слоя жестко связанной структурированной воды вокруг одного кристаллита; 2t - расстояние между двумя кристаллитами; а, Ь - толщина и ширина кристаллита соответственно; ЛБСЭ - фрагмент, являющийся структурной основой механической модели
Гидратная пленка (адсорбат) жестко сцепле-
на с кристаллитами (адсорбент) [11]. В результате адсорбции минимум один слой молекул адсорбата остается неподвижным при внешнем механическом воздействии [18]. В рассматриваемом случае пленка имеет несколько подобных слоев, так как молекулы воды, обладая существенным дипольным моментом, выстраиваются «в цепочки» под действием электрического поля поверхностно расположенных ионов кристаллита. Такое структурирование дополнительно стимулируется тем, что водный раствор насыщен ионами фосфата и кальция, одноименными с ионами кристаллической решетки адсорбента. В итоге структура водной пленки имеет много общего с кристаллитом и по своим механическим свойствам приближается к твердому телу, что проявляется, в том числе, увеличением плотности этого раствора. В некоторых случаях она достигает значений в 10 % и более на расстоянии 1 -2 молекулярных диаметров ад-сорбата от кристаллической поверхности адсорбента [18]. В результате по сравнению с обычным раствором гидратированный слой, содержащий ионы фосфата и кальция, ведет себя как твердое тело и демонстрирует высокое сопротивление не только отрыву от кристаллита, но и отрыву составляющих его частиц друг от друга.
По мере старения, как уже отмечалось выше, происходят изменения структуры гидратного слоя и соответствующие сдвиги кристаллической решетки минералов (нестехиометрические замены ионов фосфата на ионы карбоната и так далее). В результате снижается прочность гидрат-ного слоя, так как изменяется характер межмолекулярных связей как внутри него, так и между элементами гидратного слоя и кристаллита.
«Пористость» костных структур. Этот термин в настоящей работе носит отличный от общепринятого смысл, так как рассматривается с точки зрения механических свойств. Для количественной оценки механической пористости весь объём кости разделяется на две части. Первая вбирает в себя минеральную составляющую, лимитирующую жесткость и прочность костной ткани, вторая объединяет остальные компоненты, которые не оказывают существенного влияния на эти механические свойства и поэтому определены как «поры».
При первоначальной оценке степени пористости исходим из того, что доля минерального компонента (золы) составляет в костной ткани около 50 % по массе, то есть в 1 кг кости имеем 0,5 кг минерального компонента и 0,5 кг - органического компонента и воды. Плотность р гид-роксиапатита находится в пределах 30003200 кг/м3 [10, 15], плотность коллагена примем равной плотности воды, то есть 1000 кг/м3. В этом случае 0,5 кг минерала имеет объём 0,5/3,2 = 0,156 л, 0,5 кг органики и воды = 0,5 л, а суммарно органика, вода и минерал = 0,656 л.
Таким образом, объёмная доля минерала (0,156/0,656)х100 = 24 %.
Отношение объёма гидратного слоя к объёму кристаллита составляет 1,9 [11]. Так как относительный объём минерального компонента равен 24 %, то относительный объём гидратных слоёв 24 % х 1,9 = 46 %. Соответственно минерал и связанная с ним вода занимают 24 % + 46 % = 70 % объёма костной ткани. Оставшиеся 30 % объёма приходятся на свободную воду, органические структуры и связанную с ними воду, которые с механической точки зрения в настоящем исследовании рассматриваются как «пористость» (см. выше).
Если исходить из двухфазной механической модели I Б. Сштеу [23] и того факта, что объёмное содержание минерала и жестко связанной с ним воды в костной ткани составляет 70 %, то, следовательно, величина Eкт должна быть на уровне 70 % от значения Eап. Согласно же представленным выше экспериментальным данным отношение упомянутых модулей упругости находится в пределах от 18 % до 30 %. Этот факт легко объяснить, учитывая дополнительную податливость водных прослоек, входящих в состав минерального «монолита».
Что же в этой связи происходит при старении? Данных об изменении толщины гидратно-го слоя мы не нашли, в то время как множество авторов, как было сказано ранее, отмечают рост размера кристаллитов. В настоящей работе исходим из того, что толщина гидратного слоя остается неизменной.
Описание модели
Для дальнейших расчетов, основываясь на описанных выше параметрах, примем за норму следующие исходные размеры кристаллитов, имеющих форму прямоугольного параллелепипеда: толщина (а) - 4 нм, ширина (Ь) - 8 нм, длина (I) - 14 нм. Первоначально обратимся варианту, при котором все кристаллиты располагаются правильными рядами (рис. 1), при толщине водной плёнки между ними 2t = 1,3 нм. Локус ЛБСБ, включающий в себя один кристаллит с окружающей его водной плёнкой, является по существу наноуровневой структурой.
Продольная сила N (рис 2) сжимает или растягивает рассматриваемый локус. Исходя из гипотезы, что в этом элементарном локусе водная пленка ведёт себя подобно твёрдому телу со специфическими механическими свойствами, рассмотрим механическую стержневую модель, состоящую из двух параллельных стержней (рис. 2, а). Первый сформирован водной пленкой и представляет собой трубку прямоугольного сечения с площадью Авод и длиной I + 2t. Внутри этой трубки расположен второй стержень (с площадью сечения Аап), который состоит из двух последовательных участков, первый длиной I - кристаллит гидроксиапатита с моду-
лем упругости Еап и второй - водная прослойка длиной 2t с модулем упругости Евод. Здесь две водные плёнки, относящиеся к торцевым поверхностям, условно сведены в один слой удвоенной толщины 2t (рис. 2, б).
Продольная сила N деформирует всю рассматриваемую модель и имеет две составляющие Nвод и Nап, (рис. 2, б):
N = Neод + Naп. (1)
Сила ^^ деформирует водный столбик (водную трубку), сила Nап - внутренний составной стержень.
Под действием силы N первый и второй стержни деформируются совместно, одинаково изменяя свою длину (рис. 2, б):
А1, = А2 (2)
Для предварительных расчетов предположим, что модуль упругости водной плёнки Евод соответствует модулю объёмного сжатия воды К при температуре 40 °С и атмосферном давлении. По справочным данным имеем Евод= К = 2,3 ГПа [17]. В этом случае Евод меньше модуля Еап в 20-40 раз. Исходя из формулы (2) путём несложных выкладок получаем, что сила Nп превосходит силу Nод более чем на порядок. Это позволяет пренебречь величиной ^од и для упрощения дальнейших рассуждений принять ее равной нулю.
Рис. 2. Фрагмент ЛБСБ (рис. 1), представляющий единичный кристаллит, окруженный структурированной водой: а) - аксонометрическая проекция; б) - продольное сечение. Обозначения: I - длина кристаллита; АЪ -деформация продольной водной пленки; А12 - суммарная деформация кристаллита и поперечной водной пленки; - продольная сила, деформирующая воду, расположенную в продольном направлении; Ыап -продольная сила, деформирующая минерал и воду, расположенную в поперечном направлении; Аж - площадь поперечного сечения кристаллита; Ао - площадь поперечного сечения продольной водной пленки
Величина А12 представляет собой сумму деформаций стержня, состоящего из минерала и поперечной прослойки воды. Так как они соединены последовательно, то силы, действующие на них, одинаковы, как одинаковы и их площади сечения. Следовательно, исходя из закона Гука [5, 13] получаем:
А = КпЩаЛап) + Кп2У(ЕвоАап). (3)
В связи с тем, что мы пренебрегли силой
Ывод, имеем Ыап = N. В этом случае формула (3) принимает вид:
Л12 = М/(ЕаЛап) + Ыт/ЕоЛ^ (4) При экспериментальной оценке Ект костная ткань рассматривается в виде гомогенной структуры. В этом случае закон Гука для элементарного локуса записывается следующим образом:
Л = (Ы/Ект) ((¡+2^/(Л вод +Аап)). (5) Здесь символом Л1 обозначена продольная деформация этого локуса при условии:
Л = Л12. (6)
Исходя из этого, приравняем выражения (4) и (5) и, проведя преобразования, найдем выражение для расчета модуля упругости водной плёнки:
Евод = Ет(21/(1 + 20)((Лвод + Лап)/Лап)(1/С), (7) где:
С = 1 - (Ект/Еап)(1/(1 + 20)((Лвод + Ла1)/Лап). (8) Вычислим значения площадей поперечного сечения центрального стержня (Лап) и окружающей его трубки (Лвод):
Лап = аЬ = 4x8 = 32 нм2, Лвод = ^+Ь)^+а)-Лап = (1,3+8)(1,3+4)-32 = = 17,3 нм2.
Подставляя полученные значения площадей в формулы (7) и (8), находим:
С = 1 - 1,54(Ект/Еап)(1/(1 + 2) (9) Евод =1,54Ет2/(1 + 2ф(1/С). (10) Далее с учетом значений 2t = 1,3 нм и I = 14 нм формулы (9) и (10) принимают вид:
С = 1 - 1,409(ЕКт/Еап), (11)
Евод = 0,131Ект(1/С). (12)
Как отмечено выше, отношение Ект/Еап составляет в среднем 0,24. Подставляя это значение в формулы (11) и (12), получаем: С = 1 - 1,409 х 0,24= 0,662, Евод = 0,131Ект(1/0,662)= 0,198Ект. В дальнейших расчетах исходим из того, что Ект = 15 ГПа [26]. Поэтому, подставляя данное значение в последнюю формулу, определяем величину модуля упругости ультратонкой водной плёнки:
Евод = 0,198x15 = 2,97 ГПа. Полученное значение модуля упругости гид-ратной плёнки Евод всего на 29 % отличается от модуля объёмного сжатия свободной воды К, что вполне укладывается в интервал случайного разброса тех экспериментальных данных, которые мы положили в основу расчёта.
Подводя итоги проведённого здесь анализа, можно утверждать, что предлагаемая трёхфазная модель удовлетворительно объясняет снижение модуля упругости Ект костной ткани в 35 раз по сравнению с модулем упругости Еап гидроксиапатита. Следовательно, первое регламентирующее условие выполнено. Изменение параметров модели при старении Вернёмся к уравнениям (4) и (5), вновь приравняем друг к другу их правые части согласно
условию (2), а затем полученное уравнение разрешим относительно модуля упругости Ект костной ткани. Результат можно представить в следующем виде:
Erm=Ean((l+2t))/(l+2t Еап/ЕоЖАс/Ап+Ао)). (13) Выше было указано, что с изменением возраста от 35 до 75 лет средние значения длины l и ширины b увеличиваются примерно на 20 % и 9 % соответственно. Толщина a кристаллита и толщина 2t водной плёнки не изменяются. В итоге имеем новые значения параметров уравнения (13):
l = 14x1,2 = 16,8 нм, b = 8x1,09 = 8,72 нм,
Аап = ba = 8,72x4 = 34,9 нм2, AB0à=(b+2t)(a+2t)-Aan=(8,72+1,3)(4+1,3)-34,9 = = 18,2 нм2,
Эти новые значения параметров подставляем в формулу (13). Кроме них в формулу также вводим значение Еап = 65 ГПа [26] и значение Евод = 2,97 ГПа, полученное в результате вышеприведенных расчетов. Введение этих численных значений в формулу (13) позволяет найти теоретически рассчитанное на основании предлагаемой модели соотношение модулей упругости старческой костной ткани к гидроксиапатиту:
Ект = Еап((16,8+1,3)/ /(16,8+1,3 x65/2,97))(34,9/(34,9 + 18,2)) = 0,26Еап.
Как уже отмечалось выше, согласно экспериментальным данным при старении (после 35 лет) это отношение снижается от 0,24 до 0,215 [32]. В наших же теоретических расчетах установлено, наоборот, увеличение отношения модуля упругости костной ткани к модулю упругости гидроксиапатита с 0,24 до 0,26. Это, казалось бы на первый взгляд, несоответствие расчетных данных экспериментальным становится понятным, если учесть тот факт, что в эксперименте определение модуля упругости оценивается по достаточно большому костному фрагменту. А согласно данным Martin, R. B. [27] при старении увеличивается пористость костных структур, которая, по нашему мнению, и является причиной снижения модуля упругости на фоне увеличения хрупкости костной ткани. Жёсткость костной ткани и винтовые
структуры В предыдущем разделе рассмотрен вариант, при котором продольная ось кристаллита совпадает с направлением господствующих в кости растягивающих или сжимающих сил. Однако в естественных условиях это совпадение встречается редко в связи с наличием винтообразной иерархической организации минерального матрикса. Именно это создает условия, при которых продольное направление рассматриваемого выше элементарного локуса костной ткани не всегда совпадает с направлением господствующих растягивающих (сжимающих) усилий (рис. 3 и 4). В результате возникают касательные (сдвиговые) напряжения г,
направленные вдоль поперечной водной плёнки. В зависимости от угла р касательные напряжения могут варьировать (см. ниже).
При наличии касательных напряжений г прямоугольный элемент деформируется путём искажения прямых углов (рис. 3). Соответствующий угол искажения у принято называть углом сдвига. Величины ти усвязаны законом Гука:
у= т/О. (14)
Здесь О - модуль сдвига, который всегда существенно меньше модуля упругости Е, и поэтому проекция абсолютного сдвига на направление господствующей растягивающей (сжимающей) силы становится значимой. Эта проекция на рис. 3 обозначена символом А.
На рис. 3 изображена лишь деформация сдвига водной прослойки, хотя в гидроксиапа-тите возникают точно такие же касательные напряжения. Дело в том, что модуль сдвига упомянутого слоя весьма невелик даже в условиях значительных сил адсорбции. Именно ввиду малости модуля сдвига водного слоя имеем заметную сдвиговую деформацию этого слоя.
Рис. 3. Схематическое изображение возникновения дополнительных продольных деформаций в результате сдвига по водной прослойке. Обозначения: Ыап - внешняя нагрузка; ф - угол несовпадения продольной оси (х) кристаллита и направления действия силы; у - угол сдвига; 5 - абсолютный сдвиг; А - проекция абсолютного сдвига по направлению внешней нагрузки; у - поперечная ось
Упомянутую проекцию деформации А на направлении N следует складывать с аналогичной проекцией А' деформации поперечной водной пленки, которая была учтена выше в качестве фрагмента уравнения (4): А' = Де-
формация А' определяется нормальными напряжениями а, деформация А - касательными т. Возрастающую суммарную деформацию (А + А) следует считать причиной еще большего увеличения податливости модели (т. е. костных структур) в направлении господствующих сил N.
В этих обстоятельствах по-новому представляется роль винтовых (спиральных) структур в кости. Именно подобная структурная организация дополнительно увеличивает податливость костной ткани за счет деформации сдвига водных плёнок. В этой связи важен тот факт, что это сдвиговое деформирование структурированной воды на наноуровне возникает при лю-
бом типе деформации кости (сжатие, растяжение, скручивание, изгиб), происходящей на макроуровне. И в этих условиях, по нашему мнению, доля сдвиговых деформаций имеет один порядок с деформацией растяжения-сжатия. Повышение деформируемости костной ткани, в свою очередь, увеличивает производительность системы, обеспечивающей движение жидкости в канальцевой системе, то есть улучшает процесс поступления питательных веществ к остеоцитам и удаления от них шлаков.
Таким образом, можно говорить о создании усовершенствованной, более сложной модели элементарного локуса. Эта модель, помимо деформации растяжения-сжатия гидратной плёнки, должна учитывать её деформацию сдвига. Однако в данном случае возникают существенные технические проблемы, а именно: в настоящее время не известны не только значения модуля упругости и модуля сдвига гидратной плёнки, но даже их соотношение. В изотропном твёрдом теле отношение к = Е/О находится в пределах от 2 до 3 [5, 13]. Параметры структуры и физических свойств гидратной плёнки близки к аналогичным параметрам твердого тела, но не достигают их. Поэтому отношение к должно быть существенно больше 3. В этих условиях податливость предложенной здесь трехфазной модели будет соответствовать податливости реальной костной ткани, если принять модель упругости водной прослойки Евод превосходящим значение, которое было использовано в наших вышеприведенных расчетах Евод ~ К ~ 2ГПа. Высказанные здесь предположения не противоречат современным представлениям физики наноструктур. Однако все это требует экспериментальной проверки.
Инициация первичных локальных разрушений структуры минерального матрикса на наноуровне при изменении условий нагрузки. Для понимания механизмов инициации повреждений на наноуровне предлагаемую трёхфазную механическую модель удобнее рассматривать с помощью схемы, представленной на рис. 3 и 4. Здесь нет продольного водного слоя, потому что он практически не оказывает влияния на ход процессов деформирования в поперечной водной плёнке.
На рис. 4, а внешняя сила N приложена к модели строго по её продольной оси х. В подобной ситуации возникающая по поперечному сечению (в гидратном слое) внутренняя сила равна и параллельна внешней силе и направлена по нормали (перпендикуляру) к этому сечению. Степень воздействия внутренней силы на материал характеризуется напряжением (сила, приходящаяся на единицу площади). В данном случае - это нормальное напряжение ст, которое рассчитывается по формуле [5, 13]:
ст= ШАап. (15)
С учетом свойств гидратного слоя, располагающегося между двумя кристаллитами, необходимо еще раз подчеркнуть, что он демонстрирует высокое сопротивление отрыву или, другими словами, для его разрушения необходимо создание относительно больших нормальных напряжений. Несколько другая ситуация складывается при действии внешних сил под углом р к продольной оси кристаллита (рис. 4б). Здесь в поперечном сечении стержня возникают, как сказано выше, помимо нормальных напряжений ст еще и касательные т. Между ними в данном случае имеется следующая зависимость [5, 13]: т= 0,5ст57п2р. (16)
Рис. 4. Схематическое изображение двух вариантов направления действия нагрузки по отношению к продольной оси кристаллита: а) направление действующей силы совпадает с продольной осью кристаллита; б) направление действующей силы не совпадает с продольной осью кристаллита. Обозначения: I - длина кристаллита; 2t - толщина водной плёнки; N - внешняя нагрузка; ф - угол несовпадения продольной оси (х) кристаллита и направления действия силы; у - поперечная ось
Это крайне важно, так как в результате винтовой иерархической организации костных структур при любом направлении нагрузки значительная часть кристаллитов находится в условиях, иллюстрируемых схемой 4б. Параметры многоуровневых винтовых структур в реальной кости таковы, что угол ( может достигать 10 -20°, а в отдельных случаях 40° и более. В этих обстоятельствах из формулы (16) получаем, что касательные напряжения т могут достигать значительных величин (до половины нормальных). А в связи с тем, что водная прослойка сопротивляется срезу под действием касательных напряжений значительно слабее, чем отрыву под действием нормальных напряжений, высока вероятность локального возникновения критических ситуаций, при которых происходят эле-
ментарные акты повреждения в виде среза по водным прослойкам.
По мере роста нагрузки эти повреждения сначала возникают лишь в тех локусах, где угол наклона ( максимален. В результате нарушается локальная целостность единого минерального массива и происходит уменьшение эффективной с механической точки зрения площади поперечного сечения кости. Последнее немедленно ведёт к росту среднего уровня напряжений в костной ткани, что повышает опасность дальнейших повреждений. Кроме этого вблизи повреждений искривляются траектории господствующих растягивающих и сжимающих напряжений. Это, в свою очередь, в отдельных локусах приводит к дальнейшему увеличению угла ( (рис. 4), то есть критичность ситуации обостряется. Таким образом, можно утверждать, что один и тот же механизм, определяющий жесткость костных структур, при определенных условиях лежит в основе инициации повреждений этих же структур на наноуровне.
Снижение прочности и жесткости костной ткани при её старении. Выше было сказано, что с возрастом в кристаллической решётке гидроксиапатита происходят всякого рода нестехиометрические замены одних ионов на другие. К числу наиболее существенных факторов относится замещение ионов фосфата на ионы карбоната. Аналогичные изменения происходят в гидратном слое, представляющем собой насыщенный водный раствор, который по химическому составу близок к минеральному компоненту. Указывалось также, что гидратный слой находится в особых условиях, при которых его физические свойства напоминают твёрдое тело, более того - монокристалл. В этих обстоятельствах нестехиометрические замены вызывают дефекты его структуры, приводящие к уменьшению силы межмолекулярных связей и соответственно к уменьшению сопротивления сдвигу. В результате отмечается снижение всех прочностных и жесткостных характеристик. Именно изменение химического состава гид-ратного слоя, а не гиперминерализация костной ткани является первопричиной ухудшения прочности и уменьшения жесткости кости.
ВЫВОДЫ
1. Разработана трёхфазная наноуровневая модель элементарного локуса костной ткани, согласно которой гидратный слой воды, разделяющий кристаллиты гидроксиапатита, лежит в основе механизма, определяющего жесткость костных структур в физиологических условиях.
2. Данная модель предполагает, что этот же механизм лежит в основе инициации разрушения костных структур на наноуровне при изменении характера нагрузки.
3. Разработанная модель элементарного ло-куса костной ткани на наноуровне удовлетворительно описывает связь между модулем упругости гидроксиапатита и модулем упругости костной ткани на макроуровне.
4. Гидратный слой воды способен воспринимать не только сжимающие и растягивающие усилия, но и усилия сдвига. В этих обстоятельствах трёхфазная модель элементарного локуса обеспечивает понимание условий разрушения костной
ткани на самом низком иерархическом (молекулярном) уровне путём среза по гидратному слою.
5. Предлагаемая модель объясняет необходимость существования винтовых (спиральных) структур в костях оптимизацией функционирования конвекционного механизма движения жидкости в канальцевой системе.
6. Снижение прочности и жесткости кост-
ной ткани согласно анализу модели вызвано гетероионными заменами, происходящими в гидратном слое при старении организма, а не гиперминерализацией костной ткани.
7. Трёхфазная модель указывает на теоретические предпосылки создания новых конструкционных материалов с особенностями строения на наноуровне.
1. Лекция по остеологии. Многоуровневый характер структуры минерального матрикса и механизмы его формирования / А. С. Аврунин [и др.] // Гений ортопедии. - 2005. - № 2. - C. 89-94.
2. Уровни организации минерального матрикса костной ткани и механизмы, определяющие параметры их формирования (аналитический обзор) / А. С. Аврунин [и др.] // Морфология. - 2005. - N° 2. - C. 78-82.
3. Аврунин, А. С. Взаимосвязь морфофункциональных сдвигов на разных уровнях иерархической организации кортикальной кости при старении / А. С. Аврунин, Л. К. Паршин, А. Б. Аболин // Морфология. - 2006. - № 3. - С. 22-29.
4. Аврунин, А. С. Адаптационная модель потери губчатой кости при старении / А. С. Аврунин [и др.] // Гений ортопедии. -2007. - № 1. - С. 100-111.
5. Беляев, Н. М. Сопротивление материалов / Н. М. Беляев. - М. : Наука, 1976. - 608 с.
6. Ультраструктурная организация минерального компонента пластинчатой костной ткани у людей зрелого и старческого возраста / Ю. И. Денисов-Никольский [и др.] // Морфология. - 2002. - № 5. - С. 79-83.
7. Актуальные проблемы теоретической и клинической остеоартрологии / Ю. И. Денисов-Никольский [и др.]. - М. : Новости, 2005. - 336 с.
8. Денисов-Никольский, Ю. И. Минеральная фаза костной ткани / Ю. И. Денисов-Никольский, И. В. Матвейчук, Б А. Жилкин // Актуальные проблемы теоретической и клинической остеоартрологии. - М., 2005. - С. 72-84.
9. Докторов, А. А. Структурная организация минеральной фазы костной ткани / А. А. Докторов // Биомедицинские технологии. -
1999. - Вып. 12. - С. 42-52.
10. Матвейчук, И. В. Изучение особенностей морфогенеза в условиях изменённой функциональной нагрузки и его прикладное значение / И. В. Матвейчук, Ю. И. Денисов-Никольский А. А. Докторов // Биомедицинские технологии. - 1999. - Вып. 12. -С. 56-70.
11. Ньюман, У. Минеральный обмен кости / У. Ньюман, М. Ньюман. - М. : Иностранная литература, 1961. - 270 с.
12. Риггз, Б. Л. Предисловие / Б. Л. Риггз, Л. Д. Мелтон // Остеопороз : Этиология, диагностика, лечение : пер. с англ. - М. ; СПб.,
2000. - C. 3-5.
13. Сопротивление материалов / под ред. Б. Е. Мельникова. - 2-е изд. - СПб. : Лань, 2007. - 560 с.
14. Справочник (кадастр) физических свойств горных пород / под ред. Н. В. Мельникова. - М. : Недра, 1975. - 276 с.
15. Стрюбель, Г. Кристаллография / Г. Стрюбель, Д. Х. Юыммер. - М. : Недра, 1987. - 268 с.
16. Турчанинов, И. А. Комплексное исследование физических свойств горных пород / И. А. Турчанинов, Р. В. Медведев. - Л. : Недра, 1973. - 124 с.
17. Физические величины : справочник / под ред. И. С. Григорьева, Е. З. Мейлихова. - М. : Энергоатомиздат, 1991. - 1232 с.
18. Фридрихсберг, Д. С. Курс коллоидной химии / Д. С. Фридрихсберг. - СПб. : Химия, 1995. - 400 с.
19. Akkus, O. Age-related changes in physicochemical properties of mineral crystals are related to impaired mechanical function of cortical bone / O. Akkus, F. Adar, M. B. Schaffler // Bone. - 2004. - Vol. 34. - P. 443-453.
20. Burnell, G. Normal maturational changes in bone matrix, mineral and crystal size in the rat / G. Burnell, L. Teubner, A. Miller // Calciff Tiss. Int. - 1980. - Vol. 31, No 1. - P. 13-19.
21. Chatterji, S. Changes in structure ofhuman bone with age / S. Chatterji, J. W. Jeffery // Nature. - 1968. - Vol. 219, No 5153. - P. 482484.
22. Cooper, R. R. Morphollogy of the osteon / R. R. Cooper, J. W. Milgram, R. A. Robinson // J. Bone Joint Surg. - 1966. - Vol. 48-А, No 7. - P. 1239-1271.
23. Currey, J. D. Three analogies to explain the mechanical properties of bone / J. D. Currey // Biorheologi. - 1964. - Vol. 2. - P. 1-10.
24. Currey, J. D. Mechanical propertis of vertebrate hard tissues / J. D. Currey // Proc. Inst. Mech. Engrs. - 1998. - Vol. 212, Part H. -P. 399-411.
25. Jowsey, J. Age changes in human bone / J. Jowsey // Clin. Orthop. - 1960. - No 17. - P. 210-217.
26. Lawson, A. C. Collagen-calcium-phosphate composites / A. C. Lawson, J. T. Czernuszka // Proc. Inst. Mech. Engrs. - 1998. - Vol. 212, Part H. - P. 413-425.
27. Martin, R. B. On the significance of remodeling space and activation rate changes in bone remodeling / R. B. Martin // Bone. - 1991. -Vol. 12. - P. 391-400.
28. Robinson, R. A. The water content of bone / R. A. Robinson, S. R. Elliott // J. Bone Joint Surg. - 1957. - Vol. 39-A, No 1. - P. 167-188.
29. Sinaki, M. Реабилитация скелетно-мышечной системы / M. Sinaki // Остеопороз. Этиология, диагностика, лечение. - СПб., 2000. - C. 471-470.
30. Smith, J. W. Age Ganges in the organic fraction of bone / J. W. Smith // J. Bone Joint Surg. - 1963. - Vol. 45-B, No 4. - P. 761-769.
31. Smith, J. W. Observation on the water content of bone / J. W. Smith // J. Bone Joint Surg. - 1964. - Vol. 46-B, No 3. - P. 553-562.
32. Zioupos, P. Changes in the stiffness, strength, and toughness of human cortical bone with age / P. Zioupos, J. D. Currey // Bone. -1998. - Vol. 22, No. 1. - P. 57-66.
Рукопись поступила 02.06.2008.
