Full Text

Введение Закон Вольфа применительно к живой костной ткани звучит следующим образом: кость приспосабливает свою внешнюю форму и внутреннюю структуру к тем механическим силам, которые она должна выдержать [21-23]. В 1872 г. Вольф провозгласил, что «когда внешняя нагрузка в какой-либо кости изменяется вследствие травмы, патологических изменений или образа жизни, функциональная перестройка со временем переориентирует трабекулы таким образом, что они выстраиваются в линии по новым траекториям напряжения» [21-23]. В работах [14, 16, 17] обсуждался этот закон, и было установлено, что трабекулы выстраиваются по линиям главных напряжений, а математической формулировки закона Вольфа до настоящего времени нет. В работах [14-16, 18] было предложено определяющее соотношение перестройки костной ткани. Структура описывалась с помощью тензора структуры (fabric tensor), входящего в определяющее соотношение как тензорная переменная наряду с тензором деформации. В работе [18] экспериментально определены материальные константы для определяющего соотношения трабекулярной костной ткани. Предложенные в работе [16] кинетические уравнения для тензора структуры были построены по тем же законам механики, что использовались для получения определяющего соотношения. Рассмотренный аналитический пример, носящий методический характер, для перестройки кости при всестороннем сжатии показал достаточную эффективность предложенных уравнений. В дальнейшем этот подход был развит в работах [7, 8, 12, 19, 20], в которых на основе разработанных постановок и алгоритмов численно был рассмотрен ряд как классических задач механики, так и задач, моделирующих перестройку костной ткани в зубочелюстной системе человека. Однако опыт биомеханического моделирования показал ряд недостатков этого подхода: 1. Закон Вольфа выполняется в определяющем соотношении при сходимости в кинетических уравнениях к новому гомеостатическому состоянию. 2. Имеется система кинетических уравнений для каждой компоненты тензора. В каждом уравнении как минимум две материальные константы, которые необходимо уметь определять. 3. В настоящее время невозможно in vivo следить за историей формирования новой структуры. Единственно, что известно экспериментально: при изменении нагрузки новая структура в среднем формируется за 160 дней [16-18]. Таким образом, необходимо построить кинетическое уравнение, в которое бы входило минимальное число материальных констант. Построение кинетического уравнения для перестройки костной ткани Рассмотрим геометрическую интерпретацию закона Вольфа. Трабекулярная структура костной ткани описывается симметричным положительно определенным тензором структуры , имеющим три положительных действительных собственных значения. Напряженное состояние определяется симметричным тензором , имеющим действительные произвольные собственные значения. В произвольной декартовой системе координат поверхность Коши тензора представляет собой эллипсоид, полуоси которого (главные направления) наклонены к осям координат. У тензора в зависимости от знаков собственных значений поверхность Коши может быть эллипсоидом, однополостным или двуполостным гиперболоидом [2, 11]. Поэтому рассмотрим пространство главных напряжений [2, 11], т.е. пространство, в котором оси координат совпадают с главными осями тензора напряжений. Аналогично определяется пространство главных направлений тензора структуры Если совместить эти два пространства, то в процессе нагружения для выполнения закона Вольфа необходимо совпадение главных осей тензоров напряжений и структуры. В теории упругопластических процессов используется геометрическое представление процессов деформирования и нагружения в специальных пятимерных векторных пространствах деформации и напряжений Ильюшина [1, 3-6, 24]. Компоненты вектора деформации и компоненты вектора напряжении линейно и взаимно однозначно связаны соответственно с компонентами девиатора деформации и девиатора напряжений . При этом длины векторов и равны соответственно интенсивностям деформации и напряжений. Таким образом, разделены векторные свойства (ориентация главных осей и угла вида соответствующих тензоров определяется направлением векторов и в построенных пространствах) и скалярные свойства напряжений и деформации. Годограф вектора образует в непрерывную кривую - траекторию деформации. Траектория деформации определяет (с точностью до объемной деформации, закон изменения которой должен быть задан отдельно) всю историю деформирования. При совмещении геометрических пространств и вектор в каждой точке траектории определенным образом ориентирован по отношению к траектории деформации и тем самым составляет определенный угол с вектором . Полученные в этих совмещенных пространствах определяющие соотношения теории упругопластических процессов геометрически означают, что направление вектора напряжений в процессе деформирования меняется. Вектор стремится занять положение касательной к траектории деформирования. На рис. 1 показана эта ситуация при простом деформировании материала по траектории OAC с изломом в точке A. 0 A Рис. 1. Процесс деформирования материала В процессе нагружения для выполнения закона Вольфа угол между векторами напряжений и структуры должен стремиться к нулю, так как при совпадении направлений векторов совпадут направления главных осей тензоров напряжений и структуры. Все три инварианта тензоров сохраняют свои значения при преобразованиях вращения координатных осей в точке тела физического пространства [4, 6]. Третий инвариант тензора, т.е. угол вида тензора, не сохраняет свою величину при ортогональных преобразованиях (сохраняются значения первого и второго инвариантов), соответствующих вращению и отражению траекторий в пространстве. Поэтому воспользоваться аналогами определяющих соотношений Ильюшина и его учеников не представляется возможным. В плоских задачах третий инвариант тензора не играет никакой роли (выражается через первый и второй инварианты). Но во всех вариантах определяющих соотношений вектор напряжений ложится на касательную к траектории деформации и угол между векторами определяется процессом нагружения. Поэтому воспользуемся геометрическими основами теории упругопластических процессов для определения кинетических уравнений перестройки трабекулярной костной ткани [1, 3-6, 24]. При этом будем учитывать, что имеется определяющее соотношение для такого материала, опубликованное в работах Cowin [14-16]. Материальные константы для определяющего соотношения получены экспериментально [18]. Определяющее соотношение имеет вид (1) где , - тензоры напряжений и деформации, - плотность костного материала. Для описания истории деформирования вводятся кинетические уравнения для и . Формулируемые кинетические уравнения для тензора структуры должны быть согласованы с уравнением (1). Пусть построены девиаторные пространство напряжений и пространство структуры Совместим эти пространства (рис. 2). Рассмотрим траектории нагружения средней кривизны (для двухзвенных траекторий будут выполняться аналогичные соотношения). В работах [9-10] показано, что для траекторий нагружения с большой степенью точности вектор напряжений находится в соприкасающейся плоскости траектории деформирования. Будем предполагать, что пятимерный вектор структуры находится в соприкасающейся плоскости траектории напряжений (рис. 2). M 0 Рис. 2. Траектория напряжений и вектор структуры : - угол между векторами напряжений и структуры ; - единичный вектор напряжений; - угол между векторами структуры и единичным вектором касательной ; - единичный вектор главной нормали к траектории напряжений; - угол между векторами и По аналогии с постулатом изотропии Ильюшина для траекторий деформации и напряжений [1, 3-6] сформулируем следующее утверждение: ориентация вектора структуры относительно траектории напряжений зависит только от внутренней геометрии этой траектории и не зависит от ее положения в пространстве Второе предположение: векторы напряжений и структуры находятся в соприкасающейся плоскости в каждой точке траектории нагружения. В соответствии с законом Вольфа угол между векторами напряжений и структуры должен стремиться к нулю с течением времени t (рис. 2), т.е. (2) В общем случае правая часть в кинетическом уравнении (2) является функционалом процесса нагружения. До настоящего времени экспериментально известно, что процесс перестройки заканчивается примерно за 160 дней. Это позволяет определить единственную константу в кинетическом уравнении. Поэтому в простейшем случае, имеющем ясный физический смысл, будем считать, что Уравнение (2) определяет только изменение векторных свойств тензора структуры. Необходимо построить кинетическое уравнение, описывающее изменение со временем скалярных свойств тензора структуры (длины вектора структуры ), согласованное как с определяющим соотношением (1), так и с кинетическим уравнением (2). Для этого поступим следующим образом. Определим в каждой точке траектории напряжений М естественный репер [1, 3-6] векторы которого определяются по длине (т.е. временем нагружения t) дуги ОМ траектории напряжений и кривизнами по обобщенным формулам Френе. Длина дуги связана со временем t следующим соотношением: Здесь - компоненты вектора напряжений, Вектор структуры на основе постулата изотропии будет иметь вид где - функционалы процесса. Так как вектор структуры находится в соприкасающейся плоскости траектории напряжений (см. рис. 2), то получаем (3) где и - единичные векторы касательной и главной нормали к траектории напряжений (см. рис. 2). Определим угол между вектором структуры и вектором касательной Пусть - единичный вектор напряжений. Функционалы и представим в виде . (4) Пусть - угол вектора с вектором Очевидно, что (5) (6) Подставим (4), (6) в (3) и после преобразований получим (7) где Угол определяется равенством Теперь получим дифференциальное уравнение для изменения модуля вектора структуры Воспользуемся тем, что и формулами Френе , Так как и выполняется (5), то (8) Продифференцировав (7) по , после преобразований получим (9) где и Умножим (9) скалярно на а затем на Учитывая ортогональность единичных векторов, получим т.е. (10) Здесь и Для плоской траектории напряжений формулы Френе дают выражения Подставим эти выражения в (10), тогда получим Умножим первое равенство на а второе - на вычтем из первого второе, тогда можно исключить скорость изменения угла наклона вектора к касательной траектории напряжений: т.е. получено кинетическое уравнение, описывающее изменение модуля вектора структуры и тем самым изменение интенсивности тензора структуры Окончательно (11) Таким образом, для решения эволюционной задачи перестройки трабекулярной костной ткани имеются определяющее соотношение (1) и кинетические уравнения (2), (11), позволяющие проследить историю формирования как структуры, так и напряженно-деформированного состояния костной ткани. Для верификации предложенных кинетических уравнений в задаче о перестройке при всестороннем сжатии было использовано аналитическое решение Угол γ 0 20 40 60 80 100 120 140 160 180 176 172 168 164 160 156 152 148 144 140 136 t, день Рис. 3. Зависимость угла между векторами напряжений и структуры от времени перестройки костной ткани из работы [16]. В этом решении нагрузка сжатия плоского образца, оснащенного начальной структурой при заданной деформации, поворачивалась на 35°. Прослеживалось изменение тензора структуры. Численное решение этой задачи подробно представлено в работах [7, 8, 19, 20]. Для этого решения введены векторы напряжений и структуры и построена зависимость угла от времени перестройки (рис. 3). Таким образом, подтверждается возможность применения предложенного кинетического уравнения (2), которое совместно с уравнением (11) может быть использовано для описания эволюции трабекулярной костной ткани. Очевидно, что предположения, сделанные при формулировке кинетического уравнения, выполняются с хорошей точностью. Заключение Таким образом, построенные кинетические уравнения могут использоваться в процессе моделирования перестройки костной ткани при использовании одной материальной константы. Это дает возможность изучить историю формирования тензора структуры при различных видах нагружения.

About the authors

V. M Tverier

References

Statistics

Views

Abstract - 85

PDF (Russian) - 36

Refbacks

• There are currently no refbacks.