Основы биомеханики
..pdfгде, как очевидно,
f (Τ) = g(1Η) . (6.6)
Эта кривая не зависит от напряженно-деформированного состояния материала и может быть получена из простейших экспериментов (например, из эксперимента по одноосному растяжению (см. рис. 4.1).
Замечание. Несколько слов более подробно о понятии девиатора тензора. Любой тензор
второго ранга можно единственным способом разложить на два тензора: шаровой тензор (у которого диагональные элементы равны друг другу, а недиагональные элементы равны нулю) и девиатор (у которого сумма диагональных элементов равна нулю).
Следовательно, длялюбоготензора Т имеетместосоотношение
Т=S+D, |
(6.7) |
где S − шаровой тензор, а D − девиатор. Или
Тij = Sij + Dij , |
(6.8) |
причем компоненты тензоров Sij и Dij имеют следующий вид:
S |
|
= |
1 |
(T |
+T |
+T |
)δ |
|
=T |
δ |
|
= |
1 |
T δ |
|
. |
(6.9) |
|
3 |
|
|
|
|
||||||||||||
|
ij |
|
11 |
22 |
33 |
|
ij |
cp |
|
ij |
|
3 kk |
ij |
|
|
||
|
|
|
|
|
Dij =Tij |
– Tcp δij . |
|
|
|
|
|
|
(6.10) |
Для наглядности тензоры Т, S и D можно представить в матричном виде.
91
|
|
|
|
|
|
|
1 |
Тkk |
|
|
|
0 |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Т11 |
|
Т13 |
|
|
3 |
|
|
|
|
|
|
|
1 0 |
0 |
|
||||||
Т12 |
|
|
|
|
1 |
|
|
|
|
|
|
|
1 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Т = |
Т21 |
Т22 |
Т23 |
, |
S = |
|
|
0 |
|
|
|
Тkk |
|
|
0 |
|
= |
|
Тkk 0 1 |
0 |
, |
|
3 |
|
|
3 |
|||||||||||||||||
|
Т31 |
Т32 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Т33 |
|
|
|
|
|
|
|
|
1 |
|
|
|
0 0 |
1 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Тkk |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
0 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
Т11 − Tcp |
|
Т12 |
|
|
|
Т13 |
|
|
|
|
|
|
||||
|
|
|
|
|
Т21 |
|
|
Т22 − Tcp |
|
|
Т23 |
|
|
|
|
|
|
||||
|
|
|
D = |
|
|
|
|
|
. |
|
|
|
|
||||||||
|
|
|
|
|
Т31 |
|
|
|
Т32 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
Т33 − Tcp |
|
|
|
|
Разделение тензора на шаровую и девиаторную составляющие во многих случаях имеет физический смысл. Так, если Т есть тензор деформации, то шаровая часть описывает изменение объема элементарного элемента, а девиатор описывает изменения формы элемента.
Разрушение элемента материала наступает тогда, когда параметр Одквиста (степень деформации) достигает предельного значения.
Параметр Одквиста (Odqvist) имеет вид
T |
T |
|
Γ = ∫ |
Ηdt = ∫ 2ξij ξij dt. |
(6.11) |
0 |
0 |
|
Условие разрушения |
|
|
|
Γ = γ*. |
(6.12) |
Здесь введен еще один инвариант: Γ − параметр Одквиста, или степень деформации.
С помощью введенных гипотез можно получить определяющие соотношения теории ползучести.
Из формулы (6.3) после проведения свертки (т.е. суммирования) по индексам i и j получаем
Sij Sij = λ2ξij ξij .
92
Далее используем формулы (6.4)
2Τ2 = λ2 Η2 |
λ = |
2Τ |
= 2g(H). |
|
Η |
||||
2 |
|
|
||
Из (6.3) получаем |
|
|
|
|
Sij |
= 2g(Η)ξij . |
|
(6.13) |
Учитывая связь g(H) и f (T) согласно (6.6), определяющие соотношения (6.13) можно записать в обращенной форме
|
|
|
|
|
ξij |
= |
1 |
f (Τ)Sij . |
|
|
|
|
|
|
(6.14) |
||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
||
Применим степенной закон Нортона–Хоффа в форме |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
Η = BΤm , |
|
|
|
|
|
|
(6.15) |
|||
тогда Η = (B Τm−1 )Τ = f (Τ)Τ. |
|
|
|
|
|
|
|
|
|
|
|
||||||
Отсюда имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
f (Τ) = B Τm−1. |
|
|
|
|
|
|
(6.16) |
||||
Или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
|
1 |
|
|
1 |
= Τ = g(Η)Η g(Η) = |
1 |
|
|
1 |
−1. |
|||||
|
Η = Τm , |
|
Η |
m |
|
|
Η |
m |
|||||||||
|
|
|
1 |
|
|
1 |
|
|
|||||||||
|
B |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
B m |
|
|
|
|
B m |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
Введем обозначения 1 =µ . Тогда m
g(H) = 1µ Hµ-1.
B
В итоге определяющие соотношения теории установившейся ползучести при использовании закона Нортона–Хоффа примут вид
Sij |
= 2 |
|
1 |
Hµ-1ξij |
, |
(6.17) |
|
Bµ |
|||||||
ξij |
= |
1 |
BTm-1Sij . |
|
(6.18) |
||
|
|
||||||
|
2 |
|
|
|
|
93
6.2. МОДЕЛЬ ТРАНСПЛАНТАТА СЕГМЕНТА АРТЕРИИ
Сегмент артерии моделируется как прямой круговой цилиндр, состоящий из однородного, несжимаемого материала, который деформируется в условиях ползучести. В целях простоты в нашем анализе мы будем игнорировать концевые эффекты вблизи опор трансплантата с естественной артерией. Начальные внутренний
и внешний радиусы – b0 |
и B0 (рис. 6.2). Длина трансплантата |
l = const(t) . Трансплантат |
составлен из искусственного биомате- |
риала и подвержен действию давления крови P = const . Исследование напряженно-деформированного состояния трансплантата в условиях линейной и нелинейной упругости было проведено ранее (см. например, работу [13]).
Рис. 6.2. Трансплантат как полый цилиндр
94
Рассмотрим изменение радиусов b(t) и B(t) в условиях ста-
ционарной ползучести, предполагая, что толщина стенки трансплантата h гораздо меньше, чем длина l ,
h = B −b << l, |
(6.19) |
и что толщина гораздо меньше радиусов трансплантата |
|
h << b. |
(6.20) |
В этой ситуации мы можем использовать мембранное уравнение Лапласа (Laplace), которое хорошо известно из курса сопротивления материалов. Это уравнение получается из рассмотрения статического равновесия элемента материала независимо от его физических свойств. В данном случае для цилиндрической тонкостенной оболочки это уравнение имеет вид
|
σϕ |
= |
|
P |
(6. 21) |
||
|
|
|
|
|
, |
||
|
R |
|
|||||
|
|
|
h |
|
|||
где σϕ − окружное (или кольцевое) нормальное напряжение, R − |
|||||||
радиус срединной поверхности оболочки. |
|
||||||
Из (6.21) с учетом условия (6.20) получим |
|
||||||
σϕ = |
PR |
>> P. |
(6.22) |
||||
|
|||||||
|
|
|
h |
|
|
|
Поскольку радиальное напряжение σr изменяется от σr = −P при r = b до σr = 0 при r = B , то, следовательно, σr P , и мы можем пренебречь σr в сравнении с σϕ . Для простоты мы также пренебрегаем продольным напряжением σz . Из-за симметрии задачи
касательные напряжения равны нулю, и в итоге мы имеем относительное напряженное состояние с единственной ненулевой компонентой σϕ .
Вычислим далее компоненты девиатора напряжений
95
Sr |
= Sz |
= − |
σϕ |
, |
|
Sϕ = σϕ − |
σϕ |
= |
2 |
|
σϕ. |
(6.23) |
|||
3 |
|
3 |
3 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Интенсивность касательных напряжений согласно (6.4) |
|
||||||||||||||
Τ = |
|
1 |
Sij Sij = |
|
1 |
(Sr2 + Sϕ2 + Sz2 ) = |
σϕ |
. |
(6.24) |
||||||
2 |
2 |
|
|||||||||||||
|
|
|
|
|
|
|
|
3 |
|
|
Далее используем определяющее соотношение (6.14) и (6.18)
ξϕ = |
1 |
f (Τ)Sϕ = |
1 |
BΤm−1Sϕ = |
1 |
BΤm−1 |
2 |
σϕ = |
1 |
BΤm−1σϕ. |
(6.25) |
|
|
|
|
|
|||||||
2 |
2 |
2 |
3 |
3 |
|
Вычислим скорость окружной деформации ξϕ (под деформацией здесь понимается логарифмическая деформация).
|
dεϕ |
|
2πR(t) |
|
ξϕ = |
|
, εϕ = log |
|
, |
d t |
|
|||
|
|
2πR0 |
ξϕ = 1 d R .
R d t
Подставляя (6.24) и (6.26) в (6.25), получим
1 dR |
|
1 |
|
σm−1 |
B |
|
|||
|
|
|
= |
|
B |
ϕ |
σϕ = |
|
σϕm . |
|
|
|
|
|
3m+1 |
||||
R d t 3 |
3m−1 |
|
Подстановка σϕ из (6.21) дает соотношение
1 dR |
= |
B |
PR m |
||||
|
|
|
|
|
|
. |
|
R d t |
|
|
|||||
|
3m+1 |
h |
(6.26)
(6.27)
(6.28)
Для решения данного дифференциального уравнения нужно найти связь между толщиной h и радиусом R . Это можно сделать с помощью условия несжимаемости, которое имеет место при деформациях пластичности и ползучести.
96
2π R h = 2π R h , |
h = |
|
R0 h0 |
. |
(6.29) |
||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
|
|
|
|
|
|
|
|
|
R |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Выражение (6.29) для h подставляем в соотношение (6.28). |
|||||||||||||||||||||||||||||
|
1 dR |
= |
|
|
B |
|
|
|
PR2 m |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
(6.30) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
R d t |
|
|
|
|
|
3m+1 |
R0 h0 |
|
||||||||||||||||||||
Из (6.30), с учетом того, что ξϕ |
= |
1 |
|
dR |
, найдем скорость ок- |
||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R d t |
|
||||
ружной деформации ξϕ0 при t = 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
m |
|
||||||||||
|
ξϕ0 = |
|
|
|
|
|
|
PR0 |
. |
(6.31) |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
3m+1 |
h0 |
|
|
|
|||||||||||||
Подстановка (6.31) в (6.30) дает |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
1 dR |
|
|
|
|
|
R |
2m |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
= ξϕ0 |
|
|
|
. |
(6.32) |
|||||||||||||
|
|
|
R d t |
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
R0 |
|
|
|
|
|
|
|
|
|
|
||||||||||||
Интегрируем дифференциальное уравнение (6.32) методом раз- |
|||||||||||||||||||||||||||||
деления переменных с учетом начального условия R = R0 |
при t = 0 . |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
R |
2m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
1 − |
|
|
|
0 |
|
= 2mξϕ0t. |
(6.33) |
|||||||||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Формула (6.33) дает закон возрастания радиуса R(t) при пол- |
|||||||||||||||||||||||||||||
зучести трансплантата артерии. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
При разрушении ( R → ∞ ) получим |
|
||||||||||||||||||||||||||||
|
|
|
|
|
T = |
1 |
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
|
(6.34) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
2mξϕ0 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где ξϕ0 определяется по формуле (6.31).
97
Рис. 6.3. Разрушение трансплантата согласно закону Нортона–Хоффа
Формула (6.34) представляет полное время вязкого разрушения рассматриваемого трансплантата в условиях установившейся ползучести, описываемой степенным законом Нортона–Хоффа (рис. 6.3).
6.3. КОНТРОЛЬНЫЕ ВОПРОСЫ
1.Что означает термин «многоосное напряженно-деформиро- ванное состояние»?
2.Что означает гипотеза единой кривой?
3.Можете ли вы привести пример имплантатов живых (твердых и мягких) тканей?
98
ГЛАВА 7. ОСНОВНЫЕ ПРИНЦИПЫ ПОСТРОЕНИЯ ОПРЕДЕЛЯЮЩИХ СООТНОШЕНИЙ ДЛЯ ЖИВЫХ ТКАНЕЙ И БИОМАТЕРИАЛОВ
По определению, определяющие соотношения для сплошной среды вводятся как соотношения между напряженным и деформированным состояниями бесконечно малого элемента среды, зависящие также от ряда других параметров (температура, химический состав, структура, внешние электромагнитные поля и т.д.). Установление определяющего соотношения является важнейшей задачей механики. В ряде случаев эти соотношения хорошо известны, например линейно упругое тело, ньютоновкая жидкость и др. Однако даже для неживых систем, например при сложных нагружениях, эти соотношения до сих пор известны не полностью и являются предметом многочисленных исследований. Еще более сложной задачей является установление определяющих соотношений для живых систем (in vitro и in vivo). Поэтому данная глава будет посвящена этому важному вопросу.
Возможные формы определяющих соотношений ограничены некоторыми фундаментальными постулатами, основанными на физических и интуитивных соображениях. Основные принципы для механических и биомеханических материалов (неживые и живые ткани) формулируются ниже [14, 13, 15].
Прежде всего, введем важные понятия «система отсчета» и «система координат». Часто эти понятия используют как синонимы. Имеются, однако, важные отличия.
По определению, система отсчета – это твердое тело, где неподвижно находится наблюдатель. Система отсчета неподвижна относительно наблюдателя, это тело может двигаться относительно других тел с другими возможными наблюдателями.
В каждой системе отсчета можно ввести одну или больше систем координат, которые можно использовать для локализации точки
99
в этом теле. По определению, система координат в данной системе отсчета фиксирована в этой системе.
На рис. 7.1 системы отсчета А и В имеют прикрепленные к ним системы координат. Система В движется относительно системы А.
Рис. 7.1. Подвижная (справа) и неподвижная (слева) системы отсчета
Рассмотрим вектор a , фиксированный в системе отсчета A , которая содержит декартову систему координат ( x , y , z ) и цилиндрическую полярную систему координат ( r , θ, z ), возможно имеющих различные начала. Относительно декартовой системы координат вектор a имеет компоненты ax , ay , az , и относительно цилиндри-
ческой системы вектор a имеет компоненты ar , aθ , az . Различные компоненты описывают один и тот же вектор a . Преобразование компонент от декартовой к цилиндрической системе есть пример не зависящего от времени преобразования координат.
Далее рассмотрим другую систему отсчета В, которая движется поступательно со скоростью u и вращается с угловой скоростью ω относительно системы отсчета А (см. рис. 7.1). Вектор a фиксирован, т.е. геометрически инвариантен в системе отсчета А, но для наблюдения в системе отсчета В он кажется изменяющимся по направлению. Относительно декартовой системы координат ( x1 , y1 , z1 ), фиксированной в системе отсчета В, компоненты изменяются со време-
100