Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Основы биомеханики

..pdf
Скачиваний:
64
Добавлен:
15.11.2022
Размер:
5.47 Mб
Скачать
Рис. 6.1. Иллюстрация зависимости интенсивности скоростей деформаций сдвига от интенсивности касательных напряжений

где, как очевидно,

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 = 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 Τm1 )Τ = f (Τ)Τ.

 

 

 

 

 

 

 

 

 

 

 

Отсюда имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (Τ) = B Τm1.

 

 

 

 

 

 

(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Τm1Sϕ =

1

BΤm1

2

σϕ =

1

BΤm1σϕ.

(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

 

σm1

B

 

 

 

 

=

 

B

ϕ

σϕ =

 

σϕm .

 

 

 

 

 

3m+1

R d t 3

3m1

 

Подстановка σϕ из (6.21) дает соотношение

1 dR

=

B

PR m

 

 

 

 

 

 

.

R d t

 

 

 

3m+1

h

(6.26)

(6.27)

(6.28)

Для решения данного дифференциального уравнения нужно найти связь между толщиной h и радиусом R . Это можно сделать с помощью условия несжимаемости, которое имеет место при деформациях пластичности и ползучести.

96

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