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

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

..pdf
Скачиваний:
64
Добавлен:
15.11.2022
Размер:
5.47 Mб
Скачать

9.12. УПРАВЛЕНИЕ НАПРЯЖЕНИЯМИ И ДЕФОРМАЦИЯМИ

 

В СТАТИЧЕСКИ НЕОПРЕДЕЛИМОЙ АДАПТИВНОЙ») ФЕРМЕ

 

Данный пример иллюстрирует управление напряжениями и де-

формациями в статически неопределимой плоской ферме, а priori

дискретной структуре. Для управления используются собственные

деформации, свободные от напряжений, и нильпотентные собствен-

ные деформации. При этом теорема о декомпозиции дает новый ме-

тод решения подобных задач. При решении задачи показано, что

проблемы управления напряжениями и деформациями разделяются.

На рис. 9.12 показана ста-

 

 

 

 

 

F

 

тически

неопределимая

ферма,

 

 

 

 

 

 

имеющая

статически неопреде-

 

1

2

3

7

5

x

лимые реакции опор. Ферма

 

3

 

 

8

 

 

имеет n = 2 квадратных пролета

1

 

6

11

 

и составлена из Nm = 5n + 1 = 11

4

 

 

 

 

 

9

 

 

 

стержневых элементов. Предпо-

 

 

5

4

10

 

 

лагается, что модули упругости

 

2

 

 

6

 

стержней одинаковые,

Ei

= E ,

 

z

 

 

 

 

 

и площади поперечных

сече-

 

 

 

 

 

 

Рис. 9.12. Статически неопределимая

ний стержней равны,

A i = A

( i = 1, Nm = 11). Число узлов обо-

(«адаптивная») ферма, нагруженная

одной силой при наличии статически

значается

Nn = 6 , ранг статиче-

неопределимых реакций в опорах

скойнеопределенности системы

 

 

 

 

 

 

 

Обозначается как s = 3 (степень внутренней статической неопределенно-

сти системы, или число «лишних» стержней

sint

= 2 ; степень внешней

статической неопределенности системы, или число «лишних» опор

sext = 1 ). Число опорных реакций NR = 4 . Предполагается, что разруше-

ние системы происходит локально, когда какое-либо сжимающее

или какое-либо растягивающее усилие в стержне превзойдет соот-

ветствующее допустимое значение. Максимальное растягивающее

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

 

 

 

 

 

 

 

 

 

 

151

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

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

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

ε*u = −εF ,

(9.68)

поэтому определение базиса в подпространстве Hu

в данной задаче

не требуется. Базис же в подпространстве Hσ , наоборот, должен быть построен, чтобы дать решение задачи управления напряжениями.

9.12.1. Решение задачи теории упругости и построение базиса для нильпотентной собственной деформации

Решение задачи управления напряжениями делится на три этапа. На первом этапе решается изотермическая задача определения упругих напряжений для заданной силовой нагрузки, например для заданной силы F , указанной на рис. 9.12. На втором этапе строится базис в подпространстве нильпотентных собственных деформаций Hσ . Нильпотентные собственные деформации являются деформа-

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

Мы не рассматриваем здесь в деталях процедуру конструирования базиса нильпотентных собственных деформаций. Детальное описание теоремы о декомпозиции и алгоритмов построения базиса можно найти в работах [17, 23].

Здесь мы лишь упомянем, что каждый элемент подпространства Hσ (или каждая нильпотентная собственная деформация) может

152

быть представлен как

линейная комбинация базисных элементов.

Можно доказать, что

число базисных элементов подпространст-

ва Hσ (размерность Hσ ) равно в дискретной или дискретизирован-

ной системе степени

статической неопределимости (внутренней

и внешней) системы. Эта степень статической неопределимости равна сумме внутренней и внешней неопределимости.

В нашем случае мы имеем sint = 2 , sext =1 и s = sint + sext = 3 . Поэтому мы должны построить три базисных нильпотентных собственных деформации.

9.12.2. Постановка и решение задачи управления

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

Предполагается, что силы выпучивания в диагональных элементах пропорциональны силам выпучивания в недиагональных

элементах и что критические силы S Bj в коротких элементах равны друг другу.

S Bj

= (S B ) < 0 в недиагональных элементах,

(9.69)

S Bj

= (kB S B ) < 0

в диагональных элементах,

(9.70)

где 0 < kB <1 .

 

 

 

 

 

 

Предел текучести S Pj

одинаков для всех стержневых элементов

и предполагается, что он в

kP раз больше, чем силы выпучивания

в коротких элементах (например, k p

= 2 ).

 

 

 

S Pj = (kP S B ) > 0 ,

j =

 

.

 

 

1, Nm

(9.71)

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

153

 

F

= (1 + kB )S B .

(9.72)

 

max

 

 

 

2

 

 

 

 

 

 

 

 

 

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

задачи управления напряжениями может быть взята в виде следую-

щей квадратичной формы:

 

 

 

 

 

 

 

Φ(F , X ) = Nm

S

j

(F , X ) 2

min( X ) ,

(9.73)

 

 

 

 

 

 

 

j =1

 

 

 

 

 

 

учитывая следующие неравенства

 

 

 

 

S Bj S j (F , X ) S Pj , j =1, Nm .

(9.74)

Здесь

X ={X j , j =1, 2, ..., s}

вектор параметров оптимизации,

который в данном случае равен избыточным (выбранным нами) си-

лам в стержневых элементах и опорных реакциях, а именно:

X1 = S4 ,

X 2 = S9 , X 3 = R4z .

 

 

 

 

 

 

X/SB

 

 

 

 

 

 

 

0,6

 

 

 

 

 

 

 

 

 

 

 

 

 

X 2

 

0,4

 

 

 

 

 

 

 

0,2

 

 

 

 

 

 

 

0

 

 

 

 

 

Fmax /SB

 

 

 

 

 

 

X3

 

 

 

 

 

 

 

-0,2

 

 

 

 

 

X1

 

0,8

 

 

 

1

1,2

F/SB

0,6

 

 

 

Рис. 9.13. Аналитическое и численное решения для задачи

 

перераспределения напряжений для kB = 0, 5 и kP = 2 .

 

Аналитическое и численное решения представлены сплошной

 

и пунктирной линиями соответственно

 

154

Сравнение аналитического и численного решений показано на рис. 9.13, на котором хорошо видно соответствие.

Нильпотентные собственные деформации могут быть найдены через базисные элементы в Hσ , ε*σj , и значения параметров X j :

s=3

 

ε*σ = X j ε*σj .

(9.75)

j=1

9.13.ВЫВОД ОБОБЩЕННОЙ ФОРМУЛЫ МАЙЗЕЛЯ

Далее приводится вывод формулы, позволяющей вычислить перемещение любой точки тела, на которое наложена собственная деформация. Эта формула полезна для анализа как живых, так и неживых систем. В частности, с помощью этой формулы можно вычислить перемещение произвольной точки тела за счет ростовых деформаций. Первоначально эту формулу вывел советский ученый В.А. Майзель в 1940 году [24] применительно к перемещениям упругого тела при нагреве или охлаждении. Обобщение для произвольной собственной деформации сделали Ф. Циглер и Г. Иршик в 2001 году [24].

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

Рассмотрим одну единичную фиктивную силу (рис. 9.14). Принцип возможных перемещений дает

1 uk (ε) (r) = σij (k ) (ξ, r) εij (ξ) dV .

(9.76)

 

Здесь применены следующие обозначения: uk (ε) (r) перемеще-

ние в точке с радиусом-вектором r , вызванное собственной деформацией, в направлении k ; σij (k ) (ξ, r) напряжение в точке тела ξ , вы-

званное фиктивной единичной силой в точке r в направлении k .

155

Рис. 9.14. Схема данного тела с собственной деформацией и фиктивной силой, приложенной в точке r в направлении k

Полная деформация εij есть сумма упругой деформации εije ,

вызванной наложенной собственной деформацией ε*ij , и собственной деформации

εij = εije + ε*ij .

(9.77)

В соответствии с теоремой о декомпозиции собственной деформации любая собственная деформация ε*ij может быть разложена

на собственную деформацию, свободную от напряжений ε*(ij u ) , и

нильпотентную собственную деформацию ε*(ij σ)

ε* = ε*(u ) + ε*(σ)

ij ij ij . (9.78)

Затем учтем формулу для нильпотентной собственной деформации, что приводит к выражению

εij = εije + ε*ij = Cijkl σkl + ε*(ij u ) Cijkl σkl ,

εij = ε*(ij u ) .

Врезультате уравнение (9.76) приводится к виду

1 uk (ε) (r) = σij (k ) (ξ, r) ε*ij(u ) (ξ) dV .

(9.79)

V

 

156

Поскольку нильпотентная собственная деформация не создает перемещений, можно записать формулу (9.79) в форме (обобщенной) формулы Майзеля

1 uk (ε) (r) = σij (k ) (ξ, r) ε*ij (ξ) dV .

(9.80)

V

 

Если собственная деформация есть температурная деформация, то обобщенная формула Майзеля сводится к формуле, первоначально полученной Майзелем [24].

В самом деле, допустим, что

ε*ij = αT δij ,

(9.81)

где α − секущий коэффициент температурного расширения, T температура, δij символ Кронекера.

Тогда формула (9.80) принимает вид

uk (ε) (r) = σij (k ) (ξ, r) α T (ξ) δij

dV = σ(k ) (ξ, r) α T (ξ) dV ,

(9.82)

V

V

 

где

 

 

σ(k ) = σxx (k ) + σyy (k ) + σyy (k )

= σii (k ) = I1 (σ) = 3σm (k ) .

(9.83)

Здесь σm (k ) среднее нормальное лой в точке ξ от единичной силы в точке

Далее рассмотрим простой пример применения формулы Майзеля.

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

Предположим, что твердое тело является изотропным и неоднородным. Твердое тело закреплено в точке О.

напряжение, вызванное си- r в направлении k .

Рис. 9.15. Схема нагружения тела единичной фиктивной силой F в направлении k

157

В данном случае собственная деформация равна температурной деформации

ε*ij = εTij = αT δij .

Тогда из формулы Майзеля (9.82) получим

uk (ε) (r) = σ(k ) (ξ, r) αT (ξ) dV = α T σ(k ) (ξ, r) dV .

(9.84)

 

V

V

 

 

Далее докажем важную формулу

 

 

 

σdV = Fn rn .

 

(9.85)

 

V

n

 

 

где Fn внешние силы, rn

радиусы–векторы сил относительно

точки O .

 

 

 

 

Сначала запишем уравнение равновесия

 

 

 

σij , j + Fi

= 0, i =1, 2, 3 .

 

(9.86)

Из (9.86) после интегрирования по объему тела следует

 

 

σij , j xi dV + Fi xi dV = 0 .

 

(9.87)

 

V

V

 

 

Далее примем теорему Гаусса-Остроградского

 

σij , j xi dV = (σij xi ), j dV σij xi, j dV =

 

V

V

V

 

(9.88)

= σij xi n j dS σij δij dV = ti xi dS σ dV .

 

S

V

S

V

 

Поскольку ti = σij n j = 0

на граничной поверхности S , то из (9.87)

и (9.88) имеем

 

 

 

 

 

σ dV = F r dV .

 

(9.89)

 

V

V

 

 

Если объемные силы не распределены по объему, но являются сосредоточенными силами, то получим формулу (9.85).

158

В рассматриваемом случае

Fn rn = r,

σ(k ) dV = r .

 

n

V

 

В результате из соотношения (9.84) получим

 

uk (ε) (r) = αTr .

(9.90)

Это верно для любой точки твердого тела произвольной формы. Полученное перемещение направлено по радиусу r . Перпендикулярные компоненты перемещения u равны нулю, так как в этом случае r F = 0 .

9.14. СЛЕДСТВИЕ ИЗ ТЕОРЕМЫ О СОБСТВЕННОЙ ДЕФОРМАЦИИ, СВОБОДНОЙ ОТ НАПРЯЖЕНИЙ (ТЕОРЕМА 1)

9.14.1. Формулировка следствия

Пусть компоненты тензора температурной деформации имеют вторые производные по координатам, т.е. ε* (C 2 (V ))6 .

Обозначение C 2 (V ) говорит о том, что имеется функциональ-

ное пространство, куда входят функции, являющиеся непрерывными и имеющие первую и вторую производную по координатам облас-

ти V , которые тоже непрерывны. Обозначение (C 2 (V ))6 говорит о том,

что рассматриваемая величина ε* имеет шесть компонент, причем каждая из них принадлежит C 2 (V ) .

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

T (x, y, z, t) = a1 (t)x + a2 (t) y + a3 (t)z + a4 (t)

(9.91)

не вызывает появления температурных напряжений в однородном и изотропном твердом теле без опор (свободное твердое тело). Условие (9.91) является также необходимым условием для тела с опорами.

159

Другими словами, в этом случае собственная деформация (в данном случае температурная деформация) является собственной деформацией, свободной от напряжений.

9.14.2. Доказательство

Пусть температурные напряжения в каждой точке твердого тела равны нулю, т.е. σ(r) = 0, r V . Тогда температурная деформация должна удовлетворять условию ε* Hu (или температурная деформация должна принадлежать подпространству совместных собственных деформаций Hu ). Если существуют вторые производные компонент

тензора собственной деформации по координатам, ε* (C 2 (V ))6 , то

условия совместности деформаций имеют вид (условия совместности Сен-Венана)

 

 

 

 

 

 

2ε*

+

2ε*yy

= 2

 

2ε*xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xx

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

y2

 

x2

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

2ε*yy

+

2ε*

=

2

 

2ε*yz

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2

 

y2

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ε*

+

 

2ε*

=

 

 

2ε*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zz

 

 

 

xx

2

 

 

 

 

zx

,

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

z2

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

2

 

*

 

 

 

 

 

 

 

∂ε*

 

*

 

 

 

 

∂ε*

 

 

 

 

εxx =

 

 

 

 

yz

+

∂εzx +

 

 

 

 

 

 

xy

,

(9.92)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

y

 

 

 

z

 

 

y z x

 

 

 

 

 

 

 

 

 

 

2ε*

 

 

 

 

 

 

 

*

 

∂ε*

 

 

 

∂ε*

 

 

 

 

 

 

yy

=

 

 

 

∂εzx

+

 

 

xy

+

 

 

 

 

yz

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

z

 

 

 

 

 

 

x

 

 

 

x z y

 

 

 

 

 

 

 

 

 

 

 

 

2

*

 

 

 

 

 

 

 

 

 

∂ε*

 

∂ε*

 

 

 

 

 

 

 

*

 

 

 

 

 

εzz

 

 

=

 

 

 

xy

+

 

 

yz

 

 

+ ∂εzx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

x

 

 

 

 

 

 

y

 

 

 

x y z

 

 

 

 

 

 

 

 

 

 

160