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

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

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

2.3. ОПРЕДЕЛЯЮЩИЕ СООТНОШЕНИЯ ДЛЯ ИЗОТРОПНОГО РАСТУЩЕГО УПРУГОГО ТЕЛА

Предположим, что ростовые деформации развиваются в рассматриваемой области изотропно, т.е. компоненты тензоров A и B не зависят от вращения (и отражения) декартовой ортогональной системы координат.

В этом случае тензоры называются изотропными, их компоненты по определению не зависят от вращения (и отражения) осей координат. Другими словами, эти компоненты не зависят от ортогонального преобразования базисных векторов.

Известно [5], что среди тензоров второго ранга имеется только один линейно независимый изотропный тензор

Aij = δij ,

(2.7)

где δij – символы Кронекера,

 

 

1,

если i = j,

 

т.е. δij

 

 

0, если i j.

 

Среди тензоров четвертого ранга имеется три линейно не-

зависимых тензора [5]

 

 

Bijkl = δij δkl ,

Bijkl = δik δjl ± δil δjk .

 

В связи с этим предположим, что тензор Bijkl

есть линей-

ная комбинация линейно независимых тензоров (2.8)

Bijkl = λ δij δkl + μ (δik δjl

+ δil δjk )+ μ1 (δik δjl δil δjk ). (2.9)

Простое вычисление позволяет получить

 

(B σ)ij = Bijkl σkl = λ δij σkk

+ μ (σij + σji )+ μ1 (σij

σji ).

Здесь и далее применяется правило суммирования А. Эйнштейна: повторяющийся индекс означает суммирование при изменении данного индекса от 1 до 3.

21

Далее, предполагая

симметрию тензора

напряжений

( σij = σji ) и учитывая, что σkk = σ11 + σ22 + σ33 = I1 (σ)

есть первый

инвариант тензора напряжений σ , получим

 

Bijkl σkl

= λδij σkk + 2μσij .

 

Окончательно скорость деформации роста можно представить в следующем виде:

ξijg = Aδij + λδij σkk + 2μσij .

(А)

Во многих случаях удобно разложить тензор напряжений на сферический и тензор-девиатор

σij = 13 σkk δij + sij .

Вэтом случае соотношение (А) принимает вид

ξijg = Aδij + 1χ σkk δij + 21ηsij ,

где введены обозначения

= 21η, λ+ 23μ = 1χ.

Для сравнения запишем компоненты упругой деформации

εije =

1

σkk δij +

1

sij ,

 

 

 

3K

2G

а также компоненты скорости упругой деформации

e

=

d εije

=

1

 

 

+

1

s

,

ξij

dt

3K

δij

 

σkk

2G ij

где введены упругие модули

G =

E

,

K =

E

, σkk

= Kεkk .

2(1+ ν)

1

(2.

(2.12

(2.

(2.

(2.16)

22

Очевидно, что соотношения (2.12) и (2.15) довольно похожи друг на друга. Кроме члена с тензором A и заменой напряжений на скорости изменения напряжений можно отметить следующее различие формул (2.12) и (2.15).

Из второго начала термодинамики следует, что в формуле (2.15) G 0 и K 0 , что ведет к ограничениям

1ν 12.

В соотношении (2.12) для ростовых деформаций такие ог-

раничения отсутствуют.

Соотношение(2.12) можетбытьнаписановтензорнойформе

ξg =

A +

I1 (σ)

g +

1

s,

(2.18

χ

2η

 

 

 

 

 

где g – метрический тензор; s – тензор-девиатор напряжений.

Далее приведем определяющее соотношение для ростовой деформации в другой форме.

Допустим, что в соотношении (А) заложены дополнительные условия λ = 0, 2μ = B (новое обозначение).

Тогда в соответствии с формулой (2.10)

(B σ)ij = Bσij .

Соотношение (А) примет вид

ξijg = Aδij + Bσij

или в тензорной форме

ξg = Ag + Bσ.

Для одноосного напряженного состояния получим

ξg = A + Bσ.

Жизнь каждого человека проходит стадии эволюции (усложнения структуры и функций органов и тканей) и инволюции

23

(обратный процесс упрощения структуры и функций органов и тканей). На стадии эволюции нужно предположить, что A > 0 ; напротив, на стадии инволюции предполагается, что A < 0.

Предположим также, что в соотношении (2.21) B > 0. При таком предположении растяжение будет стимулировать рост ткани, а сжатие, наоборот, подавлять рост (закон Хютера– Фолькмана, Германия, XIX в.).

В настоящее время всемирно известен компрессионнодистракционный метод для лечения различных патологий костной ткани, в частности при переломах и недоразвитии конечностей (акад. Г.А. Илизаров, Курган, Россия). В данном случае управление ростовыми деформациями позволяет достичь удлинения конечностей, доходящего до 50 см. Например, на этом основано лечение ахондроплазии (нарушение роста кости, приводящее к недостаточному росту конечностей при нормальных размерах туловища).

Заметим также, что если в соотношении (2.21) положить A = 0 , то получим определяющее соотношение линейно-вязкой (или ньютоновской) жидкости.

В заключение этого параграфа отметим, что соотношение (2.22) помогает анализировать сколиоз – заболевание позвоночника, заключающеесяв егоискривлении во фронтальной плоскости.

 

Пусть имеется начальная кри-

 

визна позвоночника во фронтальной

 

плоскости. Причиной этого искрив-

 

ления может быть какая-либо асим-

 

метрия в организме (различная дли-

 

на нижних конечностей, например

 

вследствие ахондроплазии; асиммет-

 

рия в зубочелюстной системе; асим-

 

метрия в мышечной системе и т.д.)

Рис. 2.2. Искривление

(рис. 2.2).

позвоночникавофрон-

Тогда на выпуклой стороне по-

тальнойплоскости

звоночника возникают растягивающие

24

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

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

Вопрос к читателю: что можно предложить для лечения искривления позвоночника (сколиоза) с помощью биомеханических методов?

2.4. МОДЕЛЬ РАЗВИТИЯ СКОЛИОЗА

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

Дадим следующую постановку рассматриваемой проблемы. Пусть ребенок имеет ахондроплазию, приведшее к различной длине ног. В результате позвоночник наклонен на угол α к горизонту. Усилия мышц представляются парой сил с моментом М, удерживающей голову в вертикальном по-

ложении (рис. 2.3).

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

Рис. 2.3. Изгибпозвоночника Далее рассмотрим решение сле- при ахондроплазии дующеймодельной задачи.

25

 

 

В момент t = 0 вследствие действия момента М в позво-

ночнике имеет место только упругая деформация. Затем начина-

ется также ростовая деформация.

 

 

необходим для сохране-

 

 

Какой мышечный момент

M (t)

ния имеющейся конфигурации позвоночника?

 

 

 

На рис. 2.3 величина

w означает расстояние от элемента

балки до оси Ox.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Применяя классическую теорию изгиба балок Бернулли–

Эйлера, можно записать следующее:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

g

 

σ

 

 

 

 

 

ξ =

ξ0 + yk

=

ξ

 

 

+ ξ

 

=

 

 

+ A + Bσ,

(2.23

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

ξ

– продольная (вдоль оси балки)

скорость деформации;

ξ0

– продольная скорость деформации нейтральной оси позво-

ночника (как балки); у – расстояние волокна балки от нейтраль-

ной оси; k – кривизна нейтральной оси;

k – скорость измене-

ния кривизны нейтральной оси.

 

 

 

 

 

 

 

 

 

 

 

Далее вводятся силовые факторы:

 

 

 

 

M = σ y ds – изгибающий момент,

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M = σ y ds – скорость изменения изгибающего момента.

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из выражения (2.23) можно получить

 

 

 

 

 

 

 

=

 

σ

 

+ A + Bσξ0 ,

 

 

 

 

 

 

yk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

κ y2ds =

1

σ y ds + B σ y ds + (A ξ0 ) y ds ,

 

 

E

 

 

s

 

s

 

 

 

 

 

 

 

s

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

1

 

 

+ B M

,

(2.

 

 

 

 

 

 

kI

 

 

 

 

M

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

I

– момент инерции поперечного сечения; M

– изгибаю-

щий момент, равный моменту мышц.

 

 

 

 

 

26

Внашем случае имеем k = 0 .

Витоге получаем дифференциальное уравнение относительно искомой функции M (t)

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M + B M = 0.

 

(2.

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение уравнения (2.25) имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

M (t) = M0eBE t ,

 

(2.

где M0 = M0 (x) = const(x).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку k = 0 , то k =

d 2w

= M0 = const(t).

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

При t = 0

имеем

EI

d 2w

 

= const(x).

 

 

 

 

dx2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

Тогда

 

 

 

 

 

= C

w(x) = C

 

 

+ Kx + D.

dx

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее для определения констант С, K, D можно учесть

следующие граничные условия:

 

 

 

 

 

 

 

 

 

dw

 

 

 

= − tgα,

 

 

dw

 

 

= 0,

w

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

x=0

 

 

 

 

 

 

 

dx

 

x= L

 

 

 

 

 

x= L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L tgα

 

 

 

 

C = tgα, K = − tgα,

 

D =

.

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

В результате прогиб балки описывается выражением

 

 

 

 

 

 

 

w(x) =

1 tgα(x L)2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 L

 

 

 

 

 

 

 

При t = 0

 

 

 

M0

=

EI tgα

и, наконец,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M (t) = EI tgα eEB t .

L

27

В данной задаче кривизна позвоночника постоянна и равна tg(α) / L. При t = 0 она определяется упругой деформацией по-

звоночника, а далее упругая деформация уменьшается ( M (t) 0

при t → ∞ ) и неупругая часть кривизны (из-за ростовой деформации) возрастает.

Можно также сформулировать другую задачу. Пусть при t = 0 известны начальная конфигурация позвоночника и момент мышц. Как должен изменяться угол α(t), чтобы сколио-

тическая деформация не росла совсем или изменялась как можно медленнее?

В некоторых методах лечения сколиоза управление углом α(t) достигается путем ношения обуви с различной толщиной

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

2.5. РАСЧЕТ РОСТОВЫХ ДЕФОРМАЦИЙ

Имеются два независимых стержня: 1 и 2 (рис. 2.4). Каждый стержень имеет упругие и ростовые свойства: Ei – модуль Юнга ( i = 1,2 ); Si – площадь поперечного сечения; Ai , Bi

реологические коэффициенты определяющего соотношения для ростовых деформаций. При t = 0 длины стержней равны.

 

l

 

1

F1

x

 

 

2

F2

 

 

 

Рис. 2.4. Два растущих стержня

28

Каковы должны быть допустимые силы F1(t) и F2 (t), создающие одинаковые длины стержней?

l1 (t) = l2 (t), или ξ1 (t) = ξ2 (t), t.

В данной задаче скорость деформации имеет вид

 

 

1 Fi

 

Fi

 

, i = 1,2 .

(2.2

ξi

=

 

 

 

+ Bi

 

+ Ai

Ei Si

Si

 

 

 

 

 

 

Считаем, что S1, S2 – const.

Из определяющего соотношения (2.27) для обоих стержней имеем

1

 

F

 

F

 

 

 

1

 

F

 

F

 

 

 

 

1

+ B1

1

+ A1

 

=

 

2

+ B2

2

+ A2

.

(2.28)

E1 S1

S1

E2 S2

S2

 

 

 

 

 

 

 

 

Уравнение (2.28) дает соотношение между силами F1 и F2 ,

но оно недостаточно для нахождения этих сил. Для их определения нужна дополнительная информация. Например, пусть известно, что F2 (t) = 0.

Тогда из (2.28) можно получить

1

 

F

 

F

 

 

= A2 ,

 

1

+ B1

1

+ A1

 

E1

S1

S1

 

 

 

 

 

d F1 + E1 B1 F1 = ( A2 A1 )E1S1 . dt

Общее решение дифференциального уравнения (2.29) имеет вид

F1 (t) = CeE1B1t + (A2 A1 ) S1 . B1

Постоянная интегрирования C может быть найдена из условия, что начальнаядеформация равнанулю, т.е. при t = 0 F1 = 0.

29

Тогда

0

= C + (A A )

S1

 

C = −(A

A )

S1

.

 

 

 

 

2 1 B

 

2

1 B

 

 

 

 

1

 

 

 

1

 

 

 

В результате имеем

 

 

 

 

 

 

 

 

 

F1 (t) = (A2 A1 )

S1

(1eE1B1t ).

 

 

 

 

B

 

 

 

 

 

 

 

1

 

 

 

 

 

 

При A2

> A1 решение для F1 (t) показано на рис. 2.5.

F1

 

 

 

 

 

(A2 A1 )

 

S

 

 

 

 

 

 

1

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

t

 

Рис. 2.5. Зависимость силы F1 от времени при F2 = 0

Асимптотическое значение F1 = (A2

A1 )

S1

 

может быть

B

 

 

 

 

1

 

(2.29) при

сразу найдено из дифференциального уравнения

F1 (t) = const .

2.6. ИСПОЛЬЗОВАНИЕ ТЕОРИИ РОСТОВЫХ ДЕФОРМАЦИЙ ДЛЯ УЛУЧШЕНИЯ МЕТОДОВ ЛЕЧЕНИЯ ВРОЖДЕННОЙ РАСЩЕЛИНЫ ТВЕРДОГО НЁБА (ВОЛЧЬЕЙ ПАСТИ)

Расщелины твердого неба (волчья пасть) и губ (заячья губа) являются весьма частым проявлением челюстной врожденной патологии, вызывающей нарушение таких жизненно важных функций, как сосание, глотание, дыхание, речь, слух. Они имеют выраженное отрицательное косметическое проявление.

30

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