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

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

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

mar = f + (mae ) + (mac ).

Полученное соотношение (6.5) индифферентно и верно

влюбой системе отсчета.

6.8.ОБЪЕКТИВНАЯ ПРОИЗВОДНАЯ

ВМОДЕЛИ МАКСВЕЛЛА

Данные пример сформулирован впервые в книге [16]. Стержень OC вращается с постоянной угловой скоростью ω вплоскости (рис. 6.3) вокруг точки О. Стержень растягивается постоянной силой P , которая вращается вместе со стержнем. Для простоты центробежные силы не рассматриваются. Определяющее состояние для материала стержня имеет видсоотношенияМаксвелла

y

x '

 

 

C

 

 

P

 

O

ϕ = ωt

x

 

Рис. 6.3. Вращающийся стержень

ξ = Aσ+ Bσ,

где A и B – константы.

 

 

 

 

 

 

Наша цель – найти компоненту

ξxx тензора скорости де-

формации ξ , где

 

 

 

 

 

 

ξij = 1

 

vi

+

v

j

 

 

 

.

 

 

 

2

 

xj

 

 

 

 

xi

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

101

Рис. 6.4. Модель Максвелла

Тогда

ε = ε + ε

 

,

ε =

σ

,

ε

 

=

σ

,

 

 

 

η

1

2

 

1

E

 

2

 

 

где E – модель Юнга пружины; η – коэффициент вязкости демпфера.

 

 

ξ = ε = ε

+ ε

 

=

σ

+

σ

,

 

 

 

2

 

 

 

 

1

 

 

 

E

 

η

 

 

 

 

 

 

 

 

 

 

 

 

ξ = Aσ+ Bσ,

A =

1 ,

B =

1

.

(6.7)

 

 

 

 

 

 

 

η

 

 

 

 

 

E

 

Предполагая, что соотношение (6.7) может быть обобще-

 

но в тензорном виде, мы можем написать соотношение (6.6).

 

Сначала решим эту задачу с точки зрения наблюдателя,

 

который вращается вместе со стержнем.

 

Этот наблюдатель находит следующие величины:

 

σx ' x ' =

P

= const ( S = const

 

представляет собой площадь

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

поперечного сечения),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σx ' x ' = 0 ,

 

 

 

 

 

 

 

 

и в соответствии с определяющим соотношением (6.6)

 

 

 

ξx ' x ' = Aσx ' x ' = A

P

.

(6

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

Компонента ξxx может быть найдена с помощью преобразования компонент тензора.

102

ξxx = ξx ' x ' cos2 ωt = A P cos2 ωt.

S

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

 

 

 

 

σxx = σx ' x ' cos2 ωt =

P

cos2 ωt,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

=

 

dσxx

= −

2Pω

 

 

 

t

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σxx

 

 

dt

 

S

 

 

sin(ω ) cos(ω ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вновь применим определяющее соотношение (6.6)

ξxx =

A

+

B

 

=

 

AP

cos

2

ω

t

2BPω

sin(ω

t

t

). (6.10)

 

 

 

 

 

 

 

 

σxx

σxx

 

S

 

 

S

 

) cos(ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

В чем же заключается ошибка наших рассуждений? Естественно предположить, что определяющее соотноше-

ние (6.6) верно в подвижной системе отсчета, вращающейся вместе со стержнем, и производная по времени вычисляется подвижным наблюдателем. Тогда решение (6.9) верно. В нашем решении неподвижный наблюдатель «забыл», что σ в (6.6) не есть производная по времени относительно неподвижной системы отсчета, это производная по времени относительно подвижной системы. Соответствующие производные по времени относительно подвижной системы являются коротационными, или объективными, производными. Можно заметить, что тензор σ всегда является индифферентным, и тензор ξ также является индифферентным,

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

Чаще всего предполагается, что определяющее соотношение верно в системе отсчета, вращающейся вместе с представи-

103

тельным объемом с угловой скоростью ω = 0,5rot v (в соответствии с теоремой Коши–Гельмгольца о разложении движения бесконечно малого элемента сплошной среды). Тогда, решая краевую задачу в неподвижной системе отсчета, мы должны ввести соответствующую коротационную производную.

В механике существует ряд таких производных в зависимости от того, в какой подвижной системе отсчета мы считаем верным определяющее соотношение. Если такой системой является подвижная система, вращающаяся с угловой скоростью ω = 0,5rot v , то мы получим производную Яуманна–Нолла (JaumannNoll). Эта производная была первоначально введена австрийским ученым Яуманном в 1911 г. и опубликована в Вене. Затем она была введена в механику сплошной среды американским ученым Ноллом (1951).

Сначала найдем производную Яуманна–Нолла для вектора. Напомним, что эта производная есть производная по времени для наблюдателя, вращающегося с угловой скоростью ω = 0,5rot v относительно неподвижной системы отсчета. Доказательство вполне аналогично теореме сложения скоростей в классическом курсе теоретической механики.

Пусть мы имеем произвольный вектор a = ak ek (напомним, что по повторяющемуся индексу k проводится суммирование от 1 до 3), где ek – базисные векторы декартовой (для про-

стоты доказательства) системы координат в движущейся вращающейся системе отсчета.

Тогда a = ak ek + ak ek .

Здесь слагаемое ak ek есть относительная скорость вектора a (производная Яуманна–Нолла для вектора, обозначаемая как a ),

иek можетбытьнайденопоформулеПуассона ek = ω× ek .

Врезультате получим

a

=

a

e

k

=

a

a

k ω×

e

k =

a

a

(6.11

 

k

 

 

 

 

 

ω× .

104

Из тензорного анализа известно, что для каждого псевдовектора ω имеется ассоциированный тензор W второго ранга в соответствии с формулой

ω× a = W a,

a.

(6.1

Несколько слов о различии понятий вектора и псевдовектора.

Оба могут быть представлены в виде

a = ak ek

= a1e1 + a2e2 + a3e3 ,

где e1,e2 ,e3 – орты декартовой ортогональной системы коорди-

нат. Понятия вектора и псевдовектора возникают при переходе к новой декартовой ортогональной системе, т.е. при операциях поворота и отражения исходной системы координат. Компоненты вектора при переходе от системы Ох1х2х3 к системе Ох1х2х3 преобразуются согласно формуле ai′ = αik ak (сумма по k от 1 до 3),

где αik – косинусугла между «новой» осью xiи «старой» осью xk

(i, k = 1, 2, 3). Компоненты псевдовектора при повороте системы координат преобразуются так же, а при отражении они приобретают обратный знак.

 

Например,

пусть

 

мы

имеем

 

 

 

x1

 

исходную систему Ох1х2х3,

азатем

 

 

 

 

 

 

 

 

 

переходим к новой системе Ох1х2х3

 

 

 

x

 

путем отражения в плоскости Ох1х2,

 

 

 

 

 

 

 

3

 

т.е. х1=х1, х2= х2, х3= – х3 (рис. 6.5).

 

 

 

 

 

O

 

 

В

 

этом

случае

α11 = α22 = 1,

 

x

3

x

2

α33 = −1 ,

остальные αik

= 0.

Прини-

 

 

 

 

 

 

 

 

мая

формулу

ai′ = αik ak

,

получим

 

Рис. 6.5. Отражение оси

 

a′ = a , a

= a , a

= −a .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

2

2

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Именно так преобразуются, например, компоненты векто-

ра скорости точки v , т.е. v

= v , v′ = v , v

= −v . Компоненты же

 

 

 

 

 

 

 

 

 

1

1

2

2

 

3

3

 

псевдовектора ω = 0,5rot v

преобразуются иначе, они приобре-

тают обратный знак.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

v

 

 

 

 

 

 

 

 

Например, ω3 = 1

 

 

2

1

 

.

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

x2

 

 

 

 

 

105

Для компонент вектора должно быть ω

= −ω

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

3

 

Здесь же (поскольку х1 = х1, х2 = х2, v1 = v1, v2 = v2) полу-

чим ω

=

1

v

v

=

1

v

v

 

= ω

 

.

 

 

2

2

1

 

 

 

2

1

 

 

 

 

x

2

x

x

 

 

 

3

 

x

 

 

 

 

3

 

 

 

 

 

 

1

 

2

 

 

 

 

1

 

2

 

 

 

 

 

 

Таким образом, компоненты псевдовектора при вращении системы координат преобразуются так же, как компоненты вектора, а при отражении меняют знак.

Вычисление показывает, что для вектора ω = 0,5rot v ассоциированный тензор W есть тензор вихря (антисимметричная часть тензораградиента скорости v ) с компонентами.

Wij = 1

 

vi

v

j

 

, i = 1, 2, 3,

 

 

 

 

 

 

2

 

xj

 

 

 

 

 

xi

 

или в тензорных обозначениях

W = 12 (v v).

Следовательно, из формул (6.11) и (6.12) найдем производную Яуманна–Нолла для вектора a.

a = a W a.

Аналогичное вычисление может быть сделано для тензора второго ранга T = Tijeie j .

d T

=

dTij

eie j

+ Tij

de

i

e j + Tijei

de j

=

dt

dt

 

 

dt

 

 

 

 

dt

 

= Tij

eie j + Tij

ω× ei e j + Tijei ω× e j =

= Tij

eie j + Tij

W ei e j + Tijei W e j =

= Tij eie j + W Tij ei e j Tijei e j W,

(6

(6.1

(6

106

T = Tij eie j + W T T W,

где было использовано соотношение

W a = a W = −a W,

(6.16)

здесь W – тензор транспонированный к W , а так как в данном случае W тензор антисимметричный, то W = −W .

Окончательный результат для производной Яуманна–Нолла тензоравторогоранга

T

 

 

(6.17)

 

= T W T + T W.

В заключение еще раз укажем, что существуют различные возможности выбрать подвижную систему отсчета (причем сис-

тема отсчета может быть даже деформируемой).

Это ведет

к различным коротационным производным. Эти производные

носят имена Зарембы, Олдройда, Коттера, Ривлина, Трусделла

и других ученых.

 

Вернемся к примеру (рис. 6.3).

 

Мы имеем определяющее соотношение (6.6)

 

ξ = Aσ+ Bσ,

 

которое, по предположению, верно в системе отсчета, вращаю-

щейся с угловой скоростью ω = 0,5rot v .

 

Это дает результат

 

ξxx =

AP

cos2 ωt.

(6.

 

 

S

 

Получим этот результат с точки зрения неподвижного на-

 

блюдателя. Вместо выражения (6.6) мы должны записать опре-

 

деляющее соотношение в виде

 

ξ = Aσ+ Bσ ,

(6

где

 

σ = σW σ+ σ W

(6.20

107

есть производная Яуманна–Нолла с тензором W .

W = 12 (v v) тензорвихря.

Подставляя (6.20) в (6.19), мы получаем правильное определяющее соотношение

ξ = Aσ+ B (σW σ+ σ W).

В данной задаче мы имеем

 

 

 

 

 

 

σx ' x ' =

P

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

σxx

=

P

 

cos2 ωt,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

σyy

=

P

sin2 ωt,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

σxy

=

P

sin(ωt )cos(ωt ).

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

ξxx = Aσxx + Bσxx B(Wxxσxx + Wxy σyx )+

B

(σxx xx

+ σxy

yx )

= σxx

 

+ σxx 2

xy σxy .

 

W

W

 

A

 

 

B

 

BW

 

 

 

 

 

 

 

 

v

 

v

y

 

 

 

 

 

Wxy

= 1

 

x

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

y

 

x

 

Вектор скорости точки

v = ω× r,

vx = ωy z ωz y = −ωy, vy = ωz x ωx z = ωx,

Wxy = 12 (ωω) = −ω.

108

Окончательно получим

 

 

 

 

 

ξxx

= A

P

cos2 ωt + B

P

(2ωcos(ωt ) sin(ωt ))+

 

 

 

 

S

S

 

 

 

 

+2Bω

P

cos(ωt ) sin(ωt )

,

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

ξxx

= A

P

cos2 ωt.

(6.

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

Мы вновь получили тот же результат.

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

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

1.Что такое определяющее соотношение?

2.Приведите примеры определяющих соотношений.

3.Почему важно знать определяющее соотношение?

4.Почему трудно найти определяющее соотношение в живых системах?

5.Перечислите основные принципы для определяющих соотношений.

109

ГЛАВА 7. ОПРЕДЕЛЯЮЩИЕ СООТНОШЕНИЯ ДЛЯ ЖИВЫХ ТКАНЕЙ И БИОМАТЕРИАЛОВ

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

Пусть мы имеем некоторую систему А (рис. 7.1).

«Черный ящик»

X

A

Y

 

 

Рис. 7.1. Система А с входными

ивыходными параметрами

Х– входные параметры, описывающие внешнее воздействие на систему; Y – выходные параметры системы. Х и Y – тензоры некоторого ранга.

Вобщем случае определяющее соотношение для системы А есть соотношение между Х и Y.

Y = Φ(X),

где функциональная зависимость Φ определяется процессами, проходящими в системе.

В простейшем случае пусть A – частица (материальная точка). Тогда X = f есть сила, действующая на эту частицу, величина Y должна описывать движение этой частицы как результат действия силы. В соответствии со вторым законом Ньютона

110

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