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

книги / Современные проблемы биомеханики

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
2.49 Mб
Скачать

Предположим, что каждая ткань имеет ростовую деформацию ε1* и ε*2 . Можно считать, что каждая часть позвонка как

стержень имеет первоначально одинаковую длину, но после ростовой деформации длины этих стержней различны.

Соотношение σ ε для каждого стержня имеет вид

σ1 = E1 (ε1 ε1* ),

 

σ2 = E2 (ε2 ε*2 ).

(

Эти стержни соединяются пластинами и имеют одинаковую длину. Деформация в этом состоянии может быть названа остаточной деформацией εr . Напряжения в этом состоянии –

самоуравновешенные остаточные напряжения ρ1 и ρ2 . Наконец, сила Р сжимает эти стержни и создает дополнительную сжимающую деформацию, величина которой обозначена как . Напряжения в этом состоянии равны σp1 и σp2 .

Для вычисления мы должны ввести коэффициенты Ar1 и Ar 2 – коэффициенты заполнения костной ткани. Эти коэффициенты учитывают существование пор в костной ткани, 0 Ar1 1. Из экспериментов на образцах известно, что для трабекулярной кости Ar1 0,28 ; для кортикальной кости Ar 2 1 .

Также из эксперимента известно, что модуль Юнга кортикальной ткани равен E2 = 20 ГПа = 2 104 MПa.

Тогда для трабекулярной костной ткани согласно экспериментальному закону

E1 = (Ar1 )3 E2 = (0,28)3 2 104 MПa = 0,044 104 MПa.

(3.10)

Площади поперечных сечений равны

 

 

A = 2,11 104 м2 ,

A = 1,16 104

м2 .

(3.11)

1

2

 

 

71

Остаточные напряжения в позвонке можно найти из экспериментальных данных: деформация после удаления трабекулярной ткани равна ε = 34 106 . Как видно из рис. 3.9,

 

 

 

 

ε = ε*2 εr ,

 

 

(

ρ2 = − ε E2

= −34 106 2 104 MПa = −0,68MПa.

(3.13)

Из условия равновесия можно получить

 

 

 

 

 

ρ1 A1 + ρ2 A2 = 0,

 

 

ρ

= −

ρ2 A2

 

= 0,68 1,16 MПa = 0,37MПa.

(3.14)

 

 

1

 

A1

2,11

 

 

 

 

 

 

 

 

 

 

Далее найдем напряжения в позвонке после приложения

силы Р. На рис. 3.9 видно, что

 

 

 

 

σp1 = ρ1 E1 , σp2

 

= ρ2 E2

, где

= εr εp .

(3.15)

σ

 

σ1 = E1 (ε1 − ε1* )

 

σ2 = E2 (ε2 − ε*2 )

 

 

 

 

 

 

 

 

 

 

ε1*

 

 

 

 

 

σ p1

ρ1

ε*2

ε =

l

 

 

 

 

 

 

 

 

l0

0

 

 

 

 

 

σ p2

 

 

 

 

1

 

 

 

 

 

ρ2

 

 

 

 

εr

 

 

 

 

 

 

ε p

 

 

 

 

 

P = 0

 

 

 

 

 

 

P 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Рис. 3.9. Диаграмма остаточных напряжений

72

Напомним, что решаемая задача является статически неопределенной, поэтому требуется добавить еще дополнительное условие. На основании закона Вольфа сформулируем дополнительное условие: в оптимальной конструкции истинные напряжения должны быть равны, т.е. конструкция должна быть равнонапряженной.

 

σp1

 

=

 

σp2

 

,

 

0 Ar i 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ar1

 

 

Ar 2

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σp1

 

 

=

 

A

 

=

0,28

,

 

 

 

 

 

 

 

 

r1

 

 

 

 

 

σ

p

2

 

 

A

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r 2

 

 

 

 

 

 

 

 

σp2

=

 

σp1

= 3,57σp1.

(3.17)

 

 

 

0,28

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из выражений (3.15) и (3.17) имеем

 

3,57(ρ1 E1

) = −ρ2 + E2 ,

 

= 3,57ρ1 + ρ2

 

= 0,297 104.

 

E2 + 3,57E1

 

 

 

 

 

 

σp1 = ρ1 E1

= 0,36 MПa,

 

σp2 = −3,57σp1 = −1,28 MПa.

Окончательно

P = σp1 A1 + σp2 A2 = −72,5 Н.

Это аксиальная сила, сжимающая копчиковый позвонок быка.

73

Контрольные вопросы

1.Что такое остаточные напряжения?

2.Какова роль остаточных напряжений в развитии и функционировании живых систем (в частности, человека)?

3.Можете ли вы привести примеры появления остаточных напряжений в неживых тканях?

4.Могут ли остаточные напряжения быть полезными, или это невозможно?

5.Где в вашем организме существуют остаточные напря-

жения?

74

ГЛАВА 4. ПОВРЕЖДАЕМОСТЬ

ИПЕРЕСТРОЙКА КОСТИ

4.1.НАКОПЛЕНИЕ ПОВРЕЖДЕНИЙ

Понятие повреждаемости (damage) применяется при описании накопления микроповреждений, в частности микротрещин, имеющиx место при многих периодических процессах в живых и неживых системах.

Применительно к живым тканям понятие повреждаемости впервые применили P. Prendergast и J. Taylor [11]. При описании этого процесса ими были предложены следующие гипотезы.

Первая гипотеза. Существует повреждаемость (накопление внутренних микротрещин) в живой системе (в частности, кости) даже при равновесии перестройки. Эта величина обозначается как ωRE .

Например, в простейшем опыте растяжения стержня силой Р (рис. 4.1)

 

 

 

 

 

ω=

F

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F0

 

 

 

 

 

 

 

 

 

 

 

 

 

где F0 – номинальная площадь по-

 

 

 

 

 

перечного сечения; F

– площадь

 

 

 

 

 

поперечного сечения, занятая мик-

 

 

 

 

 

ротрещинами(ненесущаянагрузку).

 

P

 

В общем случае мера повре-

Рис. 4.1. Пример

ждаемости ω может

пониматься

как, например, отношение объема,

эксперимента

занятого трещинами, к номиналь-

на повреждаемость

ному объему, т.е. 0 ω1 ( ω= 0

означает, что трещины отсутству-

75

ют; ω = 1 означает, что трещины заполняют весь объем). Понятие повреждаемости ω ввел Ю.Н. Работнов (1959).

В механике известно также другое, близкое по смыслу по-

нятие «сплошность» (Л.М. Качанов, 1957):

 

ψ = 1ω,

0 ψ1.

(4.

В нашем примере

 

 

 

 

ψ = 1ω= 1

F

=

F0 F

,

F0

F0

 

 

 

т.е. это отношение несущей нагрузку площади поперечного сечения к номинальной площади.

Стимулом перестройки можно считать изменение повреждаемости, отсчитываемое от состояния равновесия перестройки. Если изменение в накопленной повреждаемости обозначить ω, то тогда имеем

ω= ωωRE ,

где ω означает актуальную повреждаемость в микроструктуре. Аргументы, доказывающие наличие повреждаемости при

равновесии перестройки:

1)в структуре имеются микротрещины, которые можно увидеть с помощью электронного микроскопа;

2)наличие природных зон концентраторов напряжений

вкости;

3)масштабный эффект в кости (зависимость прочности от размеров образца, а именно уменьшение прочности для меньших образцов кости).

Вторая гипотеза. Известно, что костная ткань является материалом, который обладает свойством восстановления после повреждаемости. Вторая гипотеза состоит в том, что скорость восстановления повреждаемости зависит от гомеостатических напряжений и не изменяется в процессе перестройки.

76

Тогда для здорового организма, а также пренебрегая эффектом старения, можно записать уравнение

d

( ω) =

d

(ωωRE ) =

d ω

 

d ω

 

 

 

 

 

(4.4)

 

 

 

 

 

 

 

= ωωRE

,

d t

d t

d t

d t

 

t =0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где ω – скорость повреждаемости при данных напряжениях, температуре и плотности кости; ωRE – скорость восстановления.

В равновесии перестройки имеем

ωt =0

= ωRE , ωt =0 = ωRE , ω= 0.

 

 

 

 

Тот факт, что скорость восстановления ωRE есть констан-

та, не зависящая от процесса перестройки, подтверждается следующим фактом.

Скорость транспорта к данному месту в кости веществ, формирующих кость, не изменяется непосредственно при изменении напряжений, т.е. эта скорость не зависит от накопленной повреждаемости.

При равновесии перестройки имеется динамическое равновесие

ω= ωRE

между образованием повреждаемости и скоростью восстановления, так что повреждаемость в равновесии перестройки не аккумулируется.

Используя Х для обозначения геометрического параметра репозиции или резорбции кости, можно написать закон перестройки

 

d X

= C ω = C (ω

ωRE ),

(4.

 

 

 

d t

 

 

 

 

 

 

 

где C – константа перестройки.

 

 

 

 

 

 

После интегрирования

(4.6) получим

 

 

0

 

t

 

 

RE )

 

) = C

(

ω

ω

(4.7)

X (t) X (t

 

 

dt.

 

 

 

 

t0

 

 

 

 

 

77

Согласно формуле (4.7) для предсказания временного хода адаптации нужно вычислить скорость накопления повреждений, но не саму повреждаемость.

В формулах (4.6), (4.7) Х есть геометрический параметр репозиции или резорбции кости, например радиус трубчатой кости. Параметр С, вообще говоря, зависит от положения в кости.

Замечание. Уравнение (4.4) вводит изменение накопленной повреждаемости (стимул перестройки) как разность между скоростью образования повреждаемости и скоростью восстановления повреждаемости. При этом нужно указать, что скорость образования повреждаемости определяется скоростью цикла нагружения, в то время как скорость восстановления определяется физиологией кости (скоростью вовлечения клеток и транспортом минералов). Многие авторы считают, что необходимо усреднять все функции времени в течение одного дня. Другими словами, мгновенная реакцияорганизма наизменениеусловий невозможна.

4.2. ПРИЛОЖЕНИЕ МОДЕЛИ

Модель перестройки при накоплении повреждений будет проиллюстрирована с использованием простого примера – перестройки диафиза кости при уменьшении вращательной нагрузки. Диафиз кости моделируется как твердый круговой цилиндр. Эта структурная модель использовалась многими исследователями ввиду ее простоты для описания диафиза длинных костей.

4.3. ОПРЕДЕЛЕНИЕ ПОВРЕЖДАЕМОСТИ

Повреждаемость может быть определена как локальная потеря механической целостности или состояние локального трещинообразования (Качанов, 1958). Механика повреждаемости используется для предсказания постепенного разрушения инженерных структур (Леметр, Шабоше, 1990). Повреждае-

78

мость определяется рядом способов в зависимости от используемого анализа. Имеется много практических мер повреждаемости и математических форм переменных повреждаемости (Крайсинович, 1981). Например, остаточная прочность или уменьшение модуля упругости при усталости, что может быть измерено (Картер, Хейес, 1977) в кортикальной кости, могут быть включены в число переменных повреждаемости.

Здесь мы используем как переменную состояния остаточную жизнь при переменных нагрузках.

В частном случае усталости, имеющей место при амплитуде напряжений, равной σi , мы можем написать, что увели-

чение повреждаемости есть функция от ni / Ni , где Ni есть число циклов до разрушения при амплитуде циклических напряже-

ний, равной σi , и ni

– число циклов, происшедших при

амплитуде напряжений

σi .

 

ω=

ni

.

 

 

 

 

Ni

Следовательно, повреждаемость увеличивается линейно с увеличением ni от нуля (неповрежденное состояние) до 1 (со-

стояние разрушения).

Таким образом, скорость увеличения повреждаемости предполагается не зависящей от уровня повреждаемости. При истории нагружения, включающей несколько значений σi , мы можем написать

ω= m ni , i=1 Ni

где m – число ступеней нагружения. Соотношение (4.9) отра-

жает так называемый «принцип линейного суммирования повреждений» (Робертсон, 1952).

79

Скорость повреждаемости можно оценить как обратную величину по отношению к времени до усталостного разрушения.

ω= d ω d ω = 1 . d t d ni Ni

Здесь вместо времени вводится пропорциональная величина – число циклов при данной амплитуде напряжений. На самом деле

t = T ni ,

где T – период изменения силы.

Эта концепция объясняется на рис. 4.2.

При больших циклических напряжениях число циклов до разрушения мало, что ведет к большой скорости накопления повреждений (это описывается на рис. 4.2 наклоном линии ω(n)). Малые циклические напряжения создают, соответственно, малую скорость накопления повреждений.

Рис. 4.2. Накопление повреждений при различных напряжениях: N1 , N2 – количество циклов

до разрушения при σ1 и σ2

Число циклов до разрушения может быть вычислено из данных, которые получены экспериментально. Например, в работе [12] получено эмпирическое соотношение для гаверсовых

80

Соседние файлы в папке книги