книги / Современные проблемы биомеханики
..pdf2.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 ,
где введены обозначения
2μ = 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.
(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) = M0e− BE 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α e− EB 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) = Ce− E1B1t + (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 |
(1− e− E1B1t ). |
|
|
|
||||
|
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