книги / Устойчивость упругих тел при конечных деформациях
..pdfк фиксируемой форме связи между напряжениями и дефор мациями. Для иллюстрации рассмотрим соотношение
о,-/ = б,уХе,ш |
2цец, |
К, р, = const, |
(1.134) |
где 0,7 — компоненты |
тензора |
истинных |
напряжений |
(1.47); в/, — компоненты тензора |
деформаций |
Грина (1.8). |
Хотя соотношения (1.134) напоминают соотношения за кона Гука в линейной теории упругости и переходят в них при переходе к теории малых деформаций (линейный ва риант), однако они не могут быть соотношениями между напряжениями и деформациями в теории упругости ко нечных деформаций не только для гиперупругих, но и для общего упругого изотропного тела. Дело в том, что в ли нейных тензорных изотропных соотношениях (1.134) слева стоят компоненты несимметричного тензора истинных на пряжений 0,7 [36, 391, справа— симметричного тензора деформаций Грина (1.8). Таким образом, результаты, по лученные для аппарата теории упругости конечных дефор маций на базе соотношений (1.134), лишены смысла.
Простейшие обобщения закона Гука для теории упругос ти конечных деформаций можно получить, если выбрать упругий потенциал в виде (1.96). В случае линейной теории упругости потенциал (1.96) согласно упрощениям§8 гл. I переходит в упругий потенциал, соответствующий закону Гука для изотропного тела. Для теории упругости конеч ных деформаций из выражений (1.30), (1.70) и (1.96) находим
оц = 6,/Яеп„ + 2ре,/. |
(1.135) |
Заметим, что в выражении (1.135) в отличие от выраже |
|
ния (1.134) слева стоят компоненты симметричного |
тензора |
обобщенных напряжений, а не несимметричного тензора истинных напряжений.
Существует еще одно обобщение закона Гука для тео рии упругости конечных деформаций. Это так называемый потенциал гармонического типа [651. В работе [351 мате риал с потенциалом гармонического типа назван полули нейным. В этом случае упругий потенциал выбирают в виде
ф = -i-X s? + p52, |
(1.136) |
где Я, р — параметры Ляме.
41
Если перейти к главным осям тензора деформаций Грииа, то инварианты
S |
i + |
( * i - l ) + ( * з - 1); |
|
|
|
4 = |
(К - 1)г + |
(К - 1)г + |
{К - 1)а, |
( |
' |
где — 1 — относительные главные |
удлинения, |
которые |
согласно выражениям (1.12) и (1.20) представляются через
главные значения деформаций |
Грина: |
|
|
Я„ - 1 |
= V l + 2en- l . |
(1.138) |
|
В случае линейной |
теории |
упругости |
из выражений |
(1.138) согласно (1.119), если пренебречь нелинейными чле нами в последнем выражении, получаем
Л.п— 1 ® 1 -)- £„ 1 = |
= Un.ru |
($*)• (1.139) |
Из выражений (1.137) и (1.139) находим, |
что величины |
|
Si и Sg представляют собой первый и второй инварианты |
тензора деформаций линейной теории упругости, отнесен ные к главным осям. Следовательно, упругий потенциал (1.136) для линейной теории упругости совпадает с упругим потенциалом, соответствующим закону Гука для изотроп ного тела.
Хотя, насколько известно автору, и не существует экспе риментальных исследований, подтверждающих существо вание. материала, свойства которого при конечных дефор мациях описываются потенциалом (1.136), все же много численные публикации выполнены в рамках потенциала (1.136) . В связи с этим вычислим потенциал (1.136) через алгебраические инварианты тензора деформаций Грина, так как при дальнейшем изложении будем считать, что упругий потенциал для изотропного тела зависит от алге браических инвариантов тензора деформаций Грийа. В хо де вычислений проиллюстрируем то положение, что соотно шения настоящей главы позволяют выразить произвольные инварианты тензора деформаций через алгебраические ин варианты тензора деформаций Грина. Согласно выраже ниям (1.35) и (1.138) получаем
%п— 1 = |
\ ^ \ + 2е + |
2 |
I |
cos ф„ |
— 1; |
(1 И0) |
|
„ |
Г |
Т Т |
|
Vn |
|
||
ф1 ==ф; |
ф, = ф + -§-я; |
фз = ф ---- —я. |
|
42
Инварианты е, et и ф выражаются через алгебраические инварианты тензора деформаций Грина по формулам (1.34). С учетом изложенного, выразим инварианты Sj и s2 через алгебраические инварианты тензора деформаций Грина:
, = |
[ / 1 + - т А 1 + -* .у Ъ У З А я- А * х |
|
|||||
X cos 4 |
- arccos У 2 |
9А3— 9 A^A3+ 2А] |
+ |
||||
У (3J4S—4j)* |
|
||||||
|
|
|
|
|
|
||
+ [У |
+ |
|
|
|
X |
|
|
Xcos ^ |
arccos 1 /2 |
9i43 — 9A ^ |
+ 2A] |
2 |
+ |
||
|
|
|
|
V(3A3- A ^ 3 |
|
|
|
+ | д / 1 + -т а 1 + ^ - У ^ * а л- а \ X ^ |
|
||||||
* |
T\ |
arccos 1Л> |
943-9Л Л + 2Л? |
|
Л |
||
х |
с ф |
|
у М |
| -------- |
|
||
|
|
|
|
|
|
|
(1.141) |
«в = [ ^ / 1 + -§ -i4 1 + - f l / 2 ^ , - i l f x |
|
||||||
|
|
|
|
)As - 9 A IAs + |
2Af _ |
i ] 8 ( |
|
|
xcos -5- arccos V 2 |
|
|
1 + |
|||
|
|
|
|
1/(3A^~Aff |
|
||
+ W |
I + 4 ^Лг + ^ - У й У з Л '- А ? |
x |
|
||||
|
Г» |
„Й |
9i4s— 9i4,i4t + 2J4| |
|
T* |
||
х с Ц у arccos K 2 |
|
^ |
|
+ T » ) ~ » J + |
43
>»•*
+[V 1+4 А*+-fv2 V*A*- A'x
, r 5 9 Л -9 /1 А + 2/1? Xcosa arccos V 2
V ( 3 A , - A f r
(1.142)
Таким образом, согласно соотношениям (1.136), (1.141) и (1.142) выражаем потенциал гармонического типа через алгебраические инварианты тензора деформаций Грина. Рассмотрим случай плоской деформации в плоскости х^ох^. При этом Яд = 1 и инварианты st и Sj принимают следующий вид:
si = (*i — !) + (* * -1 ); ss — (^i — I)2+ |
(^-2— I)2- 0143) |
|
Относительные главные удлинения |
определяются |
по |
формулам (1.138). |
|
пло |
Главные значения тензора деформаций Грина для |
||
ской деформации можно записать в виде |
|
|
ё|,2 = 4 -(eu + еа)± -j- Y |
(еи—е2г)2 + 4е?2. (1.144) |
|
Учитывая выражения (1.30), запишем правую часть |
||
(1.144) через алгебраические инварианты |
|
|
eifi = - Y(Al ± V 2 A t — Ab. |
(1.145) |
|
Из формул (1.136), (1.138), |
(1.144) и (1.145) |
получим |
для плоской деформации потенциал гармонического типа, выраженный через алгебраические инварианты
<D= ± - { [ Y l + A1 + V 2 Ai - A 2i - l ] +
+ [ Y \ + A 1- V 2 A %- A * l - l \ } ' +
+ Р |[ |/^ 1 + Ai + |
~V |
— Ai — 1 1 + |
|
+ [ Y l + A 1 - V |
2 A |
t - A al - l J } . |
(1-146) |
44
Преобразовав выражение (1.146), найдем |
|
||||
Ф - (А, + 2р)(2 + А,) + |
к | / l |
+ 2A1+ 2A2l — 2At ~ |
|||
— 2 (Я |
р) ^ | / 1 |
+ |
Ar |
2Аг — А? -f- |
|
+ / l |
+ АХ- У |
2Аг - А \ ) . |
(1.147) |
Хотя формулы (1.136), (1.141), (1.142) и (1.147) имеют громоздкий вид, однако при вычислениях коэффициентов линеаризированных уравнений состояния для однородных начальных состояний они значительно упрощаются.
Глава II
ЛИНЕАРИЗИРОВАННЫЕ ЗАДАЧИ ТЕОРИИ УПРУГОСТИ КОНЕЧНЫХ ДЕФОРМАЦИЙ
В настоящей главе сформулируем основные линеари зированные задачи теории упругости при произвольных (конечных) начальных деформациях. Рассмотрим три со стояния тела. Первое — естественное недеформированное состояние, когда в теле отсутствуют напряжения и деформа ции, второе — начальное деформированное состояние (все величины, относящиеся к нему, будут отмечены индексом нуль), третье — возмущенное деформированное состояние, все величины которого равны сумме соответствующих ве личин начального состояния и малых величин, называе мых возмущениями. Величины возмущенного состояния отметим штрихом, возмущения не будем отмечать никаким индексом. Рассмотрим постановку задач только для обоб
щенных напряжений а*ц (1.50) и компонент тензора дефор маций Грина (1.8), поскольку к ней согласно результатам гл. I сводятся все постановки задач. В отдельных случаях будем использовать компоненты несимметричного тензора напряжений Кирхгофа. Таким образом, можно записать
оц = в$* + |
Ъ) = е°ц+ е,,-; . . . ; ит = и°т + ит. (11.1) |
Предполагая, что все величины с индексом нуль зна чительно больше возмущений, линеаризируем соотноше ния и уравнения гл. I. Во всем последующем изложении под линеаризированными уравнениями и соотношениями
45
будем понимать уравнения и соотношения для возмущений. Заметим, что если задана величина у в виде у = / (лг), то для возмущений получаем формулу
(П-2)
Формулой (11.2) будем пользоваться при вычислении возмущений. Заметим, что все величины с индексом нуль вычисляются по формулам гл. I, если везде поставить ин декс нуль.
§ 1. Геометрические соотношения
Линеаризируем формулу для определения компонент тензора деформаций Грина (1.8):
2tfj = llij -f- Ujj + U°njUnJ + Un,itlnj = (finl + Hn./) Un,i +
+ («*/• + «"./)«»./• |
(П.З) |
Учитывая выражения (II.2) и (11.3), из (1.12) — (1.15) получаем формулы для нахождения возмущений всех гео метрических объектов. Из выражений (1.30), (II.2) и (Н.З) выводим формулы для определения возмущений ал гебраических инвариантов
А1 — (6 ni + Чп.{) unj", |
Л2 = |
2 6// (6„/ -f- tin,/) Un.h |
|
я |
о о |
0 /5 |
. о \ |
**3 — «Эб/двд/(On/ "р UfiJ)
Учитывая результаты § 3 гл. I и формулы (II.4), можно, следуя (И. 2), вычислить возмущения инвариантов любой системы инвариантов. В качестве примера запишем форму лы для /х и / 2 (I. 33):
Л — 2 (6«/ + U n j) Un,i\
/а = 4 [(1 + Е°и) &U— е</] (6„/ + ’
Рассмотрим линеаризированное условие несжимаемости (I. 73). С учетом формул (1.74) и (П.З), после вычислений получаем
Go (6п/ + «»./) unj = 0. |
(И.6) |
Величины (JO можно найти по формулам (1.75) или (1.75а), если во всех величинах, входящих в эти формулы, поставить индекс нуль.
46
Рассмотрим упрощения, которые возникают в выраже ниях (II.3) — (II.6), когда начальное состояние является однородным и определяется соотношениями
и°п = 8/п (Я, — 1) х(, Я, = const. |
(II.7) |
Заметим, что для несжимаемого материала из условия (1.73) и выражений (11.7) получаем зависимость между удлинениями Я/
ЯдЯ^Яз = 1. |
(П.8) |
Из выражений (II.3) и (II.7) получаем |
|
2ец = ЯfUit, + ’kg'-j.i, |
(П.9) |
Из выражений (II.4) и (II.7) выводим |
|
= ЯпнП>п5 Ая = (Я„ |
1) ЯпИп.п\ |
|
(11. 10) |
(Яп — I)2ЯлыП,„.
При выводе выражений (11.10) было учтено, что в случае (II.7)
2 4 = 4 (Я) - |
1); |
|
А° = 4 - (ЯА |
- 3); |
|
= 4 |
(А -*-1) ( * * - !); |
(II. 11) |
|||
Л°3= - L (Я2 - |
1) (Я2 - |
1) (Я2 - |
1). |
||
Из выражений (II.5) |
и |
(II.7) |
находим |
||
|
Л = |
2Я„нп>„; |
(И. 12) |
||
/ 2= 4 (l + |
|
Ялн„.„ - |
2 (Я*- |
||
|
1) Ялнп.„. |
||||
Для несжимаемого тела из (II.7), (II.8) и (1.75а) полу |
|||||
чаем |
|
|
|
|
|
|
Go = 8//ЯГ2- |
(И.13) |
Из выражений (II.6), (II.7) и (11.13) при однородной начальной деформации определяем линеаризированное условие несжимаемости
Я Л м = 0. |
(11.14) |
Перейдем к линеаризации уравнений движения, граничных и начальных условий.
47
§ 2. Уравнения движения. Граничные и начальные условия
Учитывая выражения (II.1) и (II.2), получаем линеари зированные уравнения движения, выраженные через ком поненты несимметричного тензора напряжений Кирхгофа и соответствующие (1.57)
tim,l 4* Х т— Рит— 0, |
(11.15) |
где Хт — возмущения объемных сил, отнесенных к едини це объема в естественном состоянии; р — плотность среды
в естественном состоянии. |
|
перемещениях |
Граничные условия в напряжениях и |
||
(линеаризированные условия) имеют вид |
|
|
hmN^s, = |
Р'т; |
(II-16) |
и « к = |
°* |
(И-17) |
где Р'т — возмущения поверхностных сил, |
действующих |
в возмущенном состоянии, но отнесенных к единице площа ди поверхности в естественном состоянии; N( — проекции орта нормали к поверхности тела в естественном состоянии
на орты у(. |
|
динамических |
гра |
Линеаризированные условия для |
|||
ничных задач согласно (1.60) запишем в форме |
|
||
ит|t=0= 0; |
ит|t=r = |
0. |
(11.18) |
Линеаризированные начальные условия для динами |
|||
ческих смешанных задач принимают вид |
|
||
ит|т=о = 0; |
йт|т=о = |
°- |
(И. 19) |
Линеаризированные смешанные условия на части по верхности S a запишем в форме
(tnNi |s,—Л ) (1—6/i)—0; |
(taN t |s3—Л2) (1—6«)=0; |
|
(1/ Л к - /> !)( 1 - 6 0 = |
0; |
( 11.20) |
&пи1 |s, — 0; 6/2MS|s, — 0; |
6/3ы3|s3= 0. |
Необходимо отметить, что под объемом тела и поверх ностью S = S t + S2 + S a при использовании лагранжевых координат следует понимать объем и поверхность тела в естественном состоянии.
Учитывая (IIЛ), получаем линеаризированные уравне ния движения, выраженные через компоненты тензора обобщенных напряжений и соответствующие (1.64)
[С|„Фтп ”Ь Чт.п) Л |
|
Н" Хт |
Р&т— |
|
(П.21) |
|||
Линеаризированные граничные условия |
на S lt |
соответ |
||||||
ствующие (1.65), запишем в виде |
|
|
|
|
|
|||
• [Oln фтп + Mm.n) + |
GlnUm,n] |
|s, = |
Pm. |
|
(П.22) |
|||
Условия (11.17) — (11.19) |
остаются |
в |
силе, |
а |
условия |
|||
на S a (11.20) принимают вид |
|
|
|
|
|
|
||
( l ° m (6in + u?,n) 4* ст/лП1,п] N { \ s , — -Pi} (1 — |
6 д ) = |
0 ; |
|
|||||
{[Ощ (8jn -f- i<2,n) -f |
NI Is, — Р2} (1 |
|
fie) = |
0; |
(П .2 3 ) |
|||
{[огг„ (63л -f- Из,п) -f* с/пОз.л1 Nf |s, *— P3} (1 - |
6/3) — 0; |
|||||||
|
||||||||
6««j Is, — 0; |
6/2«aIs, = 0; 6ю«3|s, = |
0. |
|
Линеаризировав выражение (1.53), получим связь меж ду возмущениями компонент тензора напряжений Кирх гофа и тензора обобщенных напряжений
t l m = *Т/л (fi/лл "4 Ч т , п ) 4 & 1 п И т . П ’ |
(П .2 4 ) |
Рассмотрим упрощения, которые возникают для одно родных начальных состояний в виде (*11.7). В этом случае выражения (II. 15) — (11.20) остаются без изменений. В вы» ражениях же (11.21) — (11.24) следует принять
б т л 4* Ч т , п — f > m n k t o (&■)* |
(И .2 5 ) |
Перейдем к вычислению линеаризированных уравнений состояния для сжимаемых и несжимаемых тел.
§ 3. Уравнения состояния
Для линеаризированных задач теории упругости урав нения состояния составляются на основе результатов, приведенных в § 7 и 8 гл. I и § 1 гл. II. Будем требовать, что бы для гиперупругого тела упругий потенциал представ лял собой дважды непрерывно дифференцируемую функ цию, а для общего упругого изотропного тела — чтобы вхо дящие в уравнения состояния функции (см. § 8 гл. 1) были непрерывно дифференцируемыми функциями»
49
В общем случае линеаризированные уравнения состоя* ния для сжимаемого тела можно представить в виде
О/л = ^fnapWa.pI |
(11.26) |
tim — ®/mapHa,p* |
(11.27) |
Величины А*„ар и ю;тар. как нетрудно убедиться в ре зультате непосредственной проверки, представляют собой компоненты тензоров четвертого ранга. Согласно (11.24) между ними существует связь
W/meP = Фтп Ч* Чт.п) ^/лаР 4“ 6ат°/р- |
(11.28) |
Для несжимаемого упругого тела в общем случае ли неаризированные уравнения состоянияпредставим в виде
°in — р/ларИа.р -f- GQр\ |
(11.29) |
tim — K;mapWa,p Ч” Go Фтп Ч" Пт,л) р, |
(11.30) |
где р/„ар и x;map — компоненты тензоров четвертого ран га. Согласно выражениям (11.24) между ними существует связь
ЭДглосР ” Фтп + Чт.п) Р/лар Ч* балДгр- |
(Н-3 1) |
Соотношения (11.26), (11.27), (11.29) и (11.30) записаны в общем случае. Ниже вычислим величины р«„<хр, x/map, A./nap и toimap для частных случаев упругих тел.
§ 4. Определение коэффициентов уравнений состояния для сжимаемого тела
Для гиперупругого тела, учитывая формулы (1.70), (11.2), (11.26) и (11.27), получаем
"т <<w+ |
(«£■+4 |
г)(4 |
г+4 |
г) |
|
|
|
|
|
|
(П.32) |
©<maP = |
Фтп + |
Чт.п) Ф а 1 Ч~ Ив./) |
---- Ь |
X |
|
x ( 4 r + ^ |
b |
+ e ~ - K |
- i + |
4 r h |
(|,-33) |
Выражения (11.32) и (11.33) записаны для анизотропно го тела. Рассмотрим упрощения, которые возникают для
50