Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
книги / Циклические нагружения элементов конструкций..pdf
Скачиваний:
0
Добавлен:
20.11.2023
Размер:
33.39 Mб
Скачать

нились области упругих деформаций, будем иметь

Oij =

fij (я, 0S, Fi

F i , Ri

R i ,

stj =

(Pii (x,

F[ -

Fl, R\ -

R l at),

*

*

 

 

t t

т. e. величины 0ij,

можно получить, используя 0ij, 8ij, путем

замены в последних а* на оа, Ft HaFi — Fix Ri на Ri — Ri и кон­

стант Ofe на константы ah.

Таким образом, возникающие при переменном нагружении

силами F iy Ri напряжения

0 ij

и деформации е^ равны разно­

стям

"

/

*

п

/

*

 

 

&ij =

O ij —

0 { j ,

8 ij =

8 i j

8 i j

соответственно напряжений а^и деформаций e!j, существовавших перед началом разгрузки, и некоторых напряжений Oij и дефор­ маций 8{j, которые возникли бы в рассматриваемом теле при нагружении его силами F% — /4 , Ri Ri при условии, что предел текучести материала равен ов, а его упругопластические свойства

характеризуются константами Ограничения на внешние нагрузки при условии существова­

ния перед началом разгрузки областей упругих деформаций (§ 1 2 ) здесь сохраняются.

Приведем, наконец, систему экспериментов с целью определе­ ния универсальных функций и констант, введенных в этом па­ раграфе.

Случай чистого сдвига при наличии гидростатического напря­

жения Он/ 3 =

'-р имеет

место при кручении тонкостенных трубок

в камере с Давлением

р. При

этом

5,7 = ^ кроме

Si2 = a1 2 ^ 0 ;

за==0, кроме

э12 = 81 2 ^ 0 ; а =

У2а12,

э = У2 е12; о» =

—Зр. Здесь

0*2 ^ касательное напряжение, 2е42 — относительный сдвиг. В опы­ те задаются, Например, а12 и р и находятся соответствующие ве­

личины

е12 и

В этом случае, согласно (15.2),

могут быть по­

строены

графики

функции У2 е12 = ф(У2 ст12,__—Зр)

при_ различных

значениях давления р ы график функции У2а12 =

/(V2ei2, е**) при

различных значениях относительного изменения объема е**, най­

денных в

результате

экспериментов.

Кроме

того, соотношения

 

~ ел* могут быть Построены согласно (15.3):

 

 

— Р = Kzkhy Ык, / 2 е 12),

ehh= —

X (— 3р, /2<х12).

 

В случае 1)астя®С‘пня цилиндрических образцов в камере с

давлением

р Цмаем

au = а0 — р, o.i2

—• о33

=

— р; еи = е 0 — 0/3,

е*2

== е83 «

— ve0 0 /3 ; ahh = — 3 р +

а0, ehh =

в + ( 1 2 v) е0; su =

~

2 V 3>s22

= «33 =

ao/3 ; a =

{W )'U o0, эп

=

2 ( 1 -f- v) e0/3, э22 =

7 H. В. Москвитин

 

 

 

 

 

= э83 = — ( 1 v) 80/3, э = (2/3)

%( 1 - f v) е0. Здесь

с0> 0 — при­

лагаемое осевое

напряжение,

ео ^ 0 — деформация, вызывае­

мая напряжением а0, 0 — относительное

изменение

объема, выз­

ванное давлением

р, v — коэффициент

поперечной

деформации.

В опыте задаются, например, а0 и р и в результате эксперимента находятся е0, 0, v, а затем э и е^. При этом могут быть построе-

где величины a, е** и а** приведены выше.

Определенную информацию об упругопластических свойствах рассматриваемых здесь материалов могут дать и результаты опы­ тов по растяжению (аАА> 0 ) и сжатию (а * * < 0 ) цилиндрических образцов при естественном атмосферном давлении. Соответствую­ щие соотношения в случае растяжения следуют из приведенных выше при р = 0 , 0 = 0 , а в случае сжатия необходимо еще у а0 и в0 поменять знаки на обратные.

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

Укажем в заключение некоторые работы, в которых исследу­ ются вопросы, рассмотренные в этом параграфе [1, 26, 90, 184, 199, 238].

§ 16. Термопластичность при переменных нагружениях

Рассмотрим в рамках теории малых упругопластпческих де­ формаций процесс неизотермического нагружения из естествен­ ного состояния внешними силами F*, при граничном перемеще­ нии и0{. Температура Т (ж, t) отсчитывается от некоторого началь­ ного состояния Т0.

Примем, что зависимость между средним напряжением, отно­ сительным изменением объема и температурой имеет вид

 

т+т0

 

 

Он = 3/L (Т) (вц За0Т),

« 0 = 1 J

a (T )d T ,

(16.1)

где а (Г ) — коэффициент линейного

теплового

расширения.

Мо­

дуль объемной деформации К может, вообще говоря, зависеть от температуры.

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

ствующее соотношение имеет вид

Sij - slj = 2 G (Т) ( » « - 4 ) ,

(16.2)

/#

где sij, эц — девиаторы напряжений и деформаций перед началом разгрузки. Здесь отмечено, что модуль сдвига G может зависеть от температуры.

В области

где не

появлялись пластические

деформации,

справедливо соотношение

 

 

 

 

=

(16.3)

В области Qp активного нагружения из естественного состоя­ ния справедливо соотношение

Sij = -J Sij, а — Ф (э, Т), (16.4)

причем здесь учтено, что универсальная зависимость а ~ э может содержать температуру Т.

Наконец, в области

й р

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

примем соотношения ( 1 0 .1 1

), записав их в следующей форме:

S{j

Sij

(pij

9ij),

 

a* — Ф* (э*, T ),

 

с * =

((<„. -

( » „

-

 

8. „

_ ,■ „ )

_ * ;;) ) ‘'. . (,6 '5)

При

э * ^ э 3

о* = 2бэ*

и зависимость (16.5) переходит в (16.2).

Запишем соотношения на границах указанных выше областей.

На

границе

1\ зон

Q*

и Qp должно

выполняться условие э =

= эа(Т) = o*(T)/(2G(T)). На границе Г2 областей £2Р и Qe начи­ нается разгрузка; это условие запишем в виде ds/dt< 0. На гра­

нице областей Q* и Qp должно выполняться условие э* = эа(Т) = = G8(T)/2G. Отметим еще, что в областях активного нагружения соответствующие модули деформации э и э* не должны убывать: дэ/ d tX ) в области Qp,da*/d£^0 в области £2р. Приведенные соотношения позволяют различать процессы нагружения и раз­ грузки.

К выписанным уравнениям состояния присоединим очевид­ ные соотношения и граничные условия:

Gij j+

= 0,

Gijlj = Ri на Sa,

Ui = uoi на Su,

(16.6)

 

 

2e,j = uitj + ujt i.

 

 

 

 

 

Температура T(x, t)

в случае несвязанной задачи

термопла-

стичпости должна удовлетворять уравнению теплопроводности

pc ^ - =

div (х (Г ) grad Т) + W*,

Т(х,0 ) = Т°,

(16.7)

7 *

где р — плотность, с — удельная теплоемкость, к — коэффициент теплопроводности, ТУ* — интенсивность источников тепла. К (16.7) следует присоединить граничные условия, которые мы здесь не конкретизируем.

Так ставится задача термопластичности с учетом разгрузок и вторичных пластических деформаций в рамках теории малых упругопластических деформаций. При этом на изменения во вре­ мени внешних нагрузок, граничных перемещений и температуры накладываются ограничения, требующие, чтобы соответствующие траектории нагружения не относились к классу существенно сложных нагружений (п. 3 § 1 ).

Заметим, что здесь рассмотрен общий случай, когда осуще­ ствляется взаимодействие внешних нагрузок с нестационарными температурными полями. Представляет также значительный ин­ терес случай свободных тел, в которых переменные нагружения осуществляются исключительно за счет циклического изменения во времени градиентов температурных полей. При этом возникает явление термической усталости (§ 57), когда элементы конструк­ ций разрушаются после небольшого числа циклов изменения тем­ пературы. Пример циклического изменения температуры в упру­ гопластической пластинке приведен в § 18.

Рассмотрим теперь один случай переменного нагружения, ис­ следование которого сводится к задаче изотермического нагруже­ ния из естественного состояния. Пусть, начиная со времени t = внешние силы изменяются таким образом, что во всех точках

пластически

деформируемых областей тела

й р

 

происходит раз­

грузка и последующее переменное нагружение

объемными сила-

—/

 

 

„//

_.

 

при граничном

ми f i и поверхностными силами

хц (на

о0)

перемещении

Щ\ (на Su)

При этом

будем

предполагать, что за

все время разгрузки и последующего переменного нагружения температура во всех точках тела остается неизменной, совпадаю­ щей с полем температуры Т'(х) = Т(х1 tY) к моменту U начала разгрузки. (Такое условие реализуется, например, в случае, когда разгрузка и переменное нагружение осуществляются за ограни­ ченный интервал времени, в течение которого температура в каж­ дой точке тела практически не изменяется.) Обозначим соответ­ ствующие напряжения, деформации и перемещения, как и рань­

ше, через

Для

величин

перед

началом

разгрузки

сохраним обозначения а ^ , е м * .

 

 

 

 

Запишем уравнения состояния. В зонах

й* и йв

разгрузки

и упругого деформирования

справедливы

соотношения (16.2),

которые в нашем случае запишутся в виде

 

 

* « =

2G {J')a ih

э * < 5 , ( Г ) ,

(До.М

*

'

п

*

'

/

S{j = Sij

S ij f

9{j = 9\j

9\j.

В области Qpt где в процессе переменного нагружения про­ исходит изменение пластических деформаций, должны быть спра­ ведливыми зависимости (16.5)

4* = £ 4 * ,

о* = Ф * (э * ,Г ).

(16.9)

Во всех точках тела зависимость (16.1) сохраняется. Перепи­ шем ее для состояния перед началом разгрузки и для текущего состояния

<Г« = зК ( Г ) (ей - 3« 0Г

), Си = 3К ( Г ) (е'« -

За0Г '),

откуда следует

 

 

Ои =

ЗК(Т')е*{.

(16.10)

Уравнения равновесия, граничные условия и соотношение Ко­

ши для величин

 

 

даются зависимостями (10.15), (10.16)

o*jJ + F *i =

0,

F*

= Fi (х, tx) Fi (x, t), t >

t±;

o*jlj = R* на SOJ

U {=U oi на Su\

(16.11)

Ft{ = R\ (fj)

R\ (£), UQ\= UQI J) UQX(/),

(16.12)

2sij i

-f- Ujj.

 

 

 

Оценивая соотношения (16.8)— (16.12), заключаем, что ве­ личины tfij, еij есть напряжения и деформации, которые возни­ кают в некотором неоднородном упругопластическом фиктивном теле при его изотермическом нагружении из естественного со­ стояния внешними усилиями Ft , R{ при граничном перемеще­

нии и0|. Фиктивное тело геометрически совпадает с рассматрива­ емым, его упругопластические свойства характеризуются

переменными

по координатам модулем сдвига G (T'(x )),

объем­

ным модулем

K iT'ix)), универсальной зависимостью

о* =

= ф*(э*, Т') и пределом текучести эв(Т'). Если указанная зада­ ча о деформации неоднородного упругопластического тела ре­ шена, то искомые величины найдутся из соотношений

C{j Oij

C\jj 8 jj —

Sij

8ij*

(16.13)

Применим теперь

этот результат к

определению остаточных

напряжений и деформаций,

сохранившихся

в теле

после

его

полной разгрузки. В

этом

случае внешние

усилия

F ixRi

и

граничное перемещение w0i изменяются в пределах от соответ­

ствующих значений F M ^),

ДД^), zz0i(*i) до нуля, а следователь­

но, величины Fi, Ri, uoi

изменяются от нуля до значений

Fiiti), jRtUi),

WoiUt).

При этом следует различать два случая:

в

первом

случае

при F\ =

F { {t^), R\ = i? i(*i)»woi = uoi(ti)

в

указанном выше

фиктивном

теле_не возникают области пла­

стических деформаций, т. е. э* ^ эа(Т ') во всех точках тела; во втором случае при указанных выше значениях внешних па­ раметров возникают области вторичных пластических деформа­

ций,

где э* >

эа(Т').

 

 

 

 

 

В

первом

случае

остаточные

напряжения ojj и остаточные

деформации

ejj определяются формулами (16.13)

 

 

 

o?j =

o 'ij — afj,

e?j =

г'ц — efh

(16.14)

в которых

0 ij, е[j — напряжения

и деформации перед началом

разгрузки, а afj,

— напряжения и

деформации,

возникаю­

щие в указанном выше неоднородном фиктивном теле при его

упругом деформировании

внешними

силами F\ (^),

R{ (tx) при

граничном перемещении u0i(£i).

 

 

Во

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

и

e?j

также

определяются

формулами

(16.14),

однако теперь

afj

и

ef есть напряжение

и деформа­

ции, возникающие в неоднородном фиктивном теле

при его на­

гружении силами Fiiti), Riiti) с граничным перемещением Uoi(ti) при условии возникновения областей пластических де­ формаций. Условие возникновения вторичных пластических де­ формаций: а* ^ эа(Т ') хотя бы для одной точки х тела, при­ чем э* подсчитывается в указанном выше случае.

Как видим, для определения остаточных напряжений и де­ формаций в рассматриваемом случае термопластичности имеет место аналог соответствующих теорем, приведенных в § 13.

§ 17. Моментные напряжения в циклически деформируемых упругопластических телах

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

сти — предельного состояния при циклическом

изменении пла­

стических

деформаций.

Учет

влияния моментных напряжений

при оценке прочности приведен в работе

[87].

 

 

Обозначим через Оц и

тпц компоненты тензоров соответствен­

но силовых

и моментных

напряжений.

Симметричная

часть Сц

и кососимметричная часть ки тензора <5ц имеют вид

 

 

1

CTji)j

1

 

®ji)r

(17.1)

 

Cij — ~2

kij = ~2 (<7jj

 

оa — сц + ku-

 

 

(17.2)

Для тензора изгиба — кручения Ху, тензоров Сц и Ау запишем соотношения [108]

 

 

 

 

1

 

 

 

 

X {j

2

 

 

(17.3)

 

 

 

 

 

 

 

 

C\j

= c8jj ~f- S{jy

3 C =

C { { 7 Aij =

A8ij ~f“ Aij,

^A == A«.

Здесь

(Dj — компоненты

вектора вращения,

Ецк— альтернативный

тензор

[108].

 

 

 

 

 

Уравнения равновесия

с учетом

(12.2)

и

(17.3), граничные

условия и соотношение для тензора Ау запишутся в следующем виде [108]:

 

 

Cmntm

£

(jftij'jm“Ь Pitm) “I” ^n =

0j

2

&jhi l^hn ,n “b

(^aP^a^p).ft ”b

jPA]J

=

 

 

=

— *2 &jki ((7a^a),fc^j

на

TH\jlj =

Qi HB SQJ (17.4)

 

U\ == UQ%^

^iii

tt>i

 

— ®oi Ha S%n

 

 

 

Amn =

*2" ®imn (^ij.j ~f" Pi).

 

Здесь Л-,

ДЛ| как

и раньше,— объемные и поверхностные силы;

Pi, qi — массовые

и поверхностные моменты; ц0ь

©0i — заданные

на Su перемещения и вращения. (По а и р здесь тоже произво­ дится суммирование.)

Запишем теперь соотношения типа теории малых упругопла­ стических деформаций с учетом моментных напряжений при на­ гружении из естественного состояния:

a

 

 

 

 

О

т1 о

(17.5)

Sij э

 

Оц — 3КЕЦш mb = Tr^ih

где

 

 

 

 

 

xi

 

 

 

 

 

 

 

 

mfj =

у

{mi} + тя),

Mij

~2 (Mi]

^ji)»

 

 

2

0

0

^ 2 =

mij mij;

 

 

 

 

= mlniij7

 

x>ij =

~2 (Xij + Xji),

X?j = -J- (Xij

Xji)*,

2

 

0

0

 

X2 =

 

Xi =

XijXij,

 

 

К соотношениям (17.5) следует присоединить универсальные

зависимости

 

 

 

 

 

 

 

о

=

ФДэ, X i,

Х 2),

=

Ф 2(х „ х2, э),

 

 

 

т 2 =

Ф3(х 2, Xi, э).

(17.6)

При этом должны быть выполнены следующие условия [1101:

до _ д т 1

до __ д п г2

^т 1

д т 2

 

(17.7)

дх1

дэ 1

д х 2

дэ

дх2

 

*

 

 

Будем считать, что приведенная выше

задача решена, т. е.

известны величины

utl Оц,

mi}.

 

 

 

 

 

Рассмотрим процесс разгрузки и последующего переменного нагружения. Введем разности

где, как и выше, величины перед началом разгрузки отмечены штрихом, а двумя штрихами отмечены текущие величины.

Аналогично соотношениям (17.5) составим зависимости для разностей (17.8)

где

Предполагаются справедливыми следующие универсальные зависимости:

причем эти функции должны удовлетворять соотношениям типа (17.7) .

Как видим, (17.9) и (17.10) есть обобщение зависимостей (10.11) — (10.13) на случай наличия при переменных нагруже­ ниях моментных напряжений. Аналогично могут быть записаны уравнения состояния для любого л-го нагружения [1681.

Уравнения равновесия и граничные условия для разностей (17.8) записываются в форме (1J.4)

Стп,т

2

4” P i,m ) 4“ ^ п — 0»

 

 

 

 

Sa;

(17.11)

Здесь

 

 

 

 

 

*

*

*

1

*

(17.12)

Tiij

ji

 

— *2*

 

c*j = c*6ij +

Sih

Sc* =

c*u s*i = o.

 

Соотношение для тензора fc?j имеет вид

 

ктп =

*гр 6imn

 

Ч~ Pi )•

(17.13)

Опуская вопросы экспериментального определения универ­ сальных функций и условия активного нагружения (см., напри­ мер, [ 1 1 0 , 1681), сформулируем теорему, являющуюся аналогом теоремы о переменном нагружении, приведенной в § 1 2 .

Представим зависимости (17.10) в следующем виде:

(17.14)

где Ф1? Ф2, Ф3— универсальные функции (17.6), характеризую­ щие свойства материала при нагружении из естественного состояния, а}к) — константы материала. В случае отсутствия моментных напряжений соотношения (17.14) переходят в обобщен­ ный принцип Мазинга.

Сравнивая теперь соотношения (17.9) — (17.14) с соответст­ вующими зависимостями (17.4) — (12.6), справедливыми при на­ гружении из естественного состояния, заключаем, что с учетом (17.8) текущие величины оцх &ijxтцх*ij при переменном нагруже-

пп

пп

 

*' п

нии внешними усилиями t {,

а |И внешними моментами р\, q\ рав­

ны разностям соответствующих величин ст*;* в*;, fliij,

перед на­

чалом разгрузки и некоторых фиктивных сг^*,

т^х

причем

эти последние есть решения задачи о нагружении рассматривае­ мого тела из естественного состояния усилиями^— F{%R { —

и моментамир' рпх q\— дчпри условий, что масштабы осей

а, э,

772-!, т 2,

х 2 при

определении

поверхностей (17.6) изме­

нены

соответственно в а (1), а (2) и а (3)

раз.

Приведенная

теорема

позволяет избежать решения дополни­

тельной краевой задачи при переменном нагружении, если соот­ ветствующая задача при первом нагружении решена.

§18. Приложения

1 . Циклическое кручение призматического стержня овального поперечного сечения. Пусть под действием крутящего момента М' > 0 в призматическом стерже возникают напряжения а13, сг23 и деформации е13, е23.Согласно известной гипотезе, перемещения

выражаются формулами

ui =

0 х2тз,

U '2=

Q'X 1X 3,

(18.1)

из =

6 / .(*ц

^г)»

где 0 ' — крутка, /'Ui, хг) — функция, характеризующая депланацию поперечного сечения.

Уравнение равновесия и ус­ ловие совместности деформаций имеют вид

 

 

 

 

 

 

013,1 + 023,2 = 0 ,

^ g

2^

в стержне

овального

поперечного

 

— 8 13,2 ~Ь е 23,1 =

 

 

 

сечения.

 

 

На

 

контуре

поперечного

се­

 

 

 

 

чения

 

должно

выполняться ус­

ловие tg ф =

сг23/сг13, где ф — угол наклона касательной к контуру

сечения (рис. II.3).

 

 

 

 

 

 

 

 

Выпишем уравнения состояния

 

 

 

 

 

 

/

а'

а' t

.

. у

, v

(18.3)

 

а1 з — 8i3 > 023 =

®23> а

= Ф (а ).

Момент М' связан с напряжениями соотношением

 

 

 

М* = j j

(а^сггз — я2а^3) dxxdx2,

 

(18.4)

где F — поперечное сечение стержня.

Если после удаления крутящего момента стержень вновь на­ гружать моментом М " < 0 , то для определения соответствующих напряжений и деформаций можно воспользоваться теоремой в переменных нагружениях. В частности, для определения зависи­ мости, связывающей М " и соответствующую крутку 0 ", следу­

ет воспользоваться соотношением М ' = Af'(0', ав). При этом, со­ гласно теореме о переменных нагружениях,

Л/"(е") = М'(е', a J - M 'O '- e ", a20.).

(18.5)

Рассмотрим более подробно задачу о переменном кручении стержня почти эллиптического поперечного сечения, контур ко­ торого задан следующим уравнением в параметрической форме

>36):

Xi [а + Ь(cos 2 ф + 2 )] sin q>, х2= —la + Mcos 2 ф — 2 ) 1 cos ф. (18.6)

При Ь = 0 это — круг радиуса а.

Если рассматриваемый стержень деформируется только упру­

го, то для определения напряжений можно принять

 

о[ 3 = Ax2t а23 = Вхг.

(18.7)

При этом удовлетворяется уравнение равновесия (18.2). Из усло­ вия совместности деформаций (18.2) следует, что В А = 2G0'.

Напряжения (18.7) удовлетворяют точно граничным условиям лишь в четырех точках контура поперечного сечения. Удовлет­ ворим граничным условиям еще в четырех точках контура, для которых

 

 

__

я

Зя

 

 

 

 

Ф “

7 * "41

"4"’

7 е

 

 

При этом 4

’ оА =

““

Таким образом,

 

 

-A a — 26

 

 

 

 

 

 

 

 

'

670' (д

26)

*

6г0л(д — 26)

(18.8)

 

<*13 = ------------«2

 

 

a 23 =

------ ^

--------- «1 .

Так как

 

г

 

 

 

г

 

 

 

 

 

 

 

 

 

 

'

_ а 13

,

Л/

'

а 23

Л/

 

 

ИЗД — “д +

У * 2»

^8,2 =

“Я* — “ ^1 *

 

ТО

щ = 2Q’x1x2bla.

Наконец, согласно (18.4) и (18.8),

Пусть 0, — значения

крутки, при

которой на границе

конту­

ра впервые

появятся

пластические

деформации. Это

будет в

точках Xi =

0, хг = ± (я — Ъ). В этих точках |a|3|= хв1 с** = 0, по­

этому согласно (18.8)

Соответствующее значение крутящего момента будет

М а=

И - 5 ^

+ 51 $)■2 (а + 2 Ъ) (а — Ь)

Пусть теперь

М' > М„

причем будем считать, что материал

не обладает упрочнением и пластическая область полностью охва­ тывает контур поперечного сечения. В такой постановке задача весьма оригинальным методом решена В. В. Соколовским [245]. Согласно этому решению, граница, определяющая области упру­

гих и пластических деформаций, есть эллипс с полуосями сг и (рис. Н.З):

‘ ! } - 2 6 [ т а т

+ (1 + (т и ’бГ),/’ *

1 ]-

<18-9>

Внутри эллипса

 

 

 

 

 

 

+

(‘ +

Ш

)

- w

r

l

 

 

 

 

 

 

(18.10)

Вне эллипса

 

 

 

 

 

 

Г

 

/

TSsm ф,

 

(18.11)

а13 = cos ф,

а2з =

 

причем по всему сечению

х2+ xt ctg ф = 4 Ь cos ф.

Осевое смещение в упругой зоне будет

 

u3 = —

| (l +

( 4С0'Ь)

)

— 4С(ПГ]«

(18.12)

в пластической области

 

 

 

 

 

 

 

 

uJ = 80'&2c o s ( p { j - s m ( p [ l +

( l

+

( - ^ ^ - ]

)

] } . (18.13)

Крутящий момент равен

 

 

 

 

 

 

М ' =

яг, (а3 - 1-аЬ3 +

463) -

 

 

 

 

 

 

 

^ 6*{* + 2(тшУ-

 

 

 

 

 

 

*/.)

- т

-

< w

r ) l [ 1 + ( w

r )

J 1 - 118-14»

Будем теперь стержень после его разгрузки вновь нагружать моментом М " < 0. Для определения напряжений и перемещений

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

стических деформаций, определяется полуосями clf с2, выраже­ ния для которых следуют из (18.9):

и

Cl

2Ъ . 2<2(0' — %")Ъ

+ I1+

(г е Н в '- В ')))

±1\

 

=

 

При этом в соотношениях (18.9) величина

0'

заменена на

0 '—

0 " и т* на 2 т„ поскольку принято а2 = 2 .

 

 

 

 

 

 

Область

изменения пластических

деформаций

появится

при

условии с2^

а — 6 пли

 

 

 

 

 

 

 

 

 

 

 

Тг____

, r4

, (

т*

v i x/* <

i

( ±

,

л

 

2bG(0' — 0")

L

\2bG(0' - 0")/

.

^

2 l b ^

/•

 

Эта область полностью охватывает контур поперечного сече­

ния, если сг

 

+ Ь или

 

 

 

 

 

 

 

 

 

 

2bxa +

[т* +

4G2&2 (0' -

0*)2]1/s <

-

6) G (0' -

0").

 

Если 0" =

0 ',

то

с1=с'и с”2 =

с2.

Если

же

10"1<0 ', то

ci > cii С2 > С2-В этом случае в поперечном сечении будут три области.

а. Область вне эллипса с полуосями съ с2. В этой области со­ гласно (18.11)

а*3 =

2 т, cos ф,

а23 = 2 т, sin ф.

 

Поэтому

 

 

 

 

13 =

— т, cos фА

а23 =

т, sin ф.

 

Аналогичным образом, используя (18.13); получим

 

Us = — 40'Ь2 sin 2<р ( l

+ [ l + ( - щ ъ ) ]

} + Wbxy cos <p +

+ 4(0

0 ) Ь2 sin 2ср |l +

J^2C6(0/_ Q»))

+ l] j-

б. Область между эллипсами с п о л у о с я м и с 2 и си

с2. В этой

области согласно (18.10)

 

 

 

° 1 3 = T; COS ф + G (0' — 0") *2 [ 1 +

( i +

 

 

 

 

 

 

If____ 1 .

 

 

 

 

) ' - 0")J’

023 —

S1H

 

 

 

1 — l

1 + [ 2Gb (0' — 0 "j) )

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

*.

1

н согласно (18.12) и (18.13)

 

 

 

 

 

2Gb (6' — 0")J

 

 

 

 

 

 

 

 

 

Из =

80'b2 cos *

f e ~ |1+ (1+ (w )2)/*]sin<p}+

 

 

 

 

+

(0

— 0 ) Xx* 2 [ ( l +

( 2Cb(0' _

0")) )

2Gb(0' — 0 ")]'

в. Область внутри эллипса с полуосями

с'2.

В

этой области

согласно (18.10)

 

 

 

 

 

 

 

 

 

 

 

0^3 = - GVX2[1 + (1 + (-45^) )

-4ёьг] +

 

 

 

 

 

+ G (0>- 0 >) ^ [ i +

(i +

U

(0 > - r ))2) '

' 26b (0' — 0 ").

аа =

GQ'Xl [ l

-

( l

+ (-jgggr) )

+

4 c ^ r ] -

 

 

 

 

 

 

 

 

 

 

 

4GbQ

 

 

 

 

 

 

 

— G (0 ' — 0 ' ) * x [ l — ( l

 

h — Y )

 

2Gb (O' — 0").

 

 

 

 

 

 

 

2Gb (0' -

Q")j )

 

и согласно (18.12)

 

 

 

 

 

 

 

 

 

 

 

 

..................... vV,

 

 

 

 

 

 

 

 

us =

_ 0'Xlx2 ,(1 + ( ш Й )

 

4GbQ'

+

 

 

 

 

 

 

+

(0'

0 ) x xx 2

( l +

 

lGb

 

)

 

2Gb (0' — 0")..

Если 0" =

—0', то вторая область исчезает, и по всему сечению

 

 

п

 

/

//

 

t

 

tt

/

 

 

 

013 =

— 013,

023 =

— СГ23,

 

^3 = — U3-

 

 

Принимая

во

внимание

(18.5)

и учитывая

формулу

 

(18.14),

получим

 

 

 

 

 

 

 

 

 

 

 

 

м' = - тят*(а3- таЬг)~ Iят*ь#{к[~wrУ-

 

+

 

- т ( б б ( в ' - р ) ) _ [ 1 _ T ( w F ) ] [ 1 + ( w r ) ]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 11 /2

 

 

 

 

[2 — ( вь(0' — е")) j f 1

”*■ т {~Ш Ж ^ ¥ )

 

Приведенные выше соотношения, в соответствии с ограниче­ ниями применимости теоремы о переменных нагружениях, спра­ ведливы при \М" |< М'

На рис. II.4 приведена графическая зависимость между без­ размерным крутящим моментом и безразмерной круткой за один

м

5

Рпс. П.4. Зависимость между крутящим моментом и круткой при цикличе­ ском кручении.

цикл кручения стержня при Ъ/а = 0,05, Ъ/а = 0,15. Точками отме­ чены границы участков, на которых пластическая зона не охва­ тывает полностью внешний контур сечения. Эти участки на гра­

фике построены на

основании условия непрерывности кривых

М = Ж 0).

 

2 .

Вторичные

Пластические деформации в толстостенной сфе­

рической оболочке. Толстостенная оболочка (полый шар) с радиу­ сом внутренней поверхности а и внешней поверхности Ъ нагружа­ ется сначала внутренним давлением р. Пусть материал шара не­ сжимаем и обладает линейным упрочнением. Если через и = = и'(г) обозначить радиальное перемещение, то вследствие усло­ вия несжимаемости

откуда

Соотношение для модуля вектора напряжений

(18.16)

вместе с уравнением равновесия

 

 

 

 

а4 + т

^ - ° > ) = °

 

 

 

 

 

дают возможность определить

и о2 с точностью до произволь­

ных постоянных:

 

 

 

 

 

 

 

 

 

 

 

« ; -

« + / «

j ® '

( з

 

 

 

 

(18.17)

Для случая линейного упрочнения

 

 

 

 

 

 

ф,('^‘)= 2<?

г><т<ъ’

 

 

 

(18.18)

 

 

 

 

 

 

 

 

 

 

 

 

Ф ' ( №

\ = ко* + 2 & ( 1 - к ) ^ ,

а ^ г ^

г,.

 

При этом в области упругих деформаций (г, *£ г <

Ь)

согласно

(18.17)

и

(18.18)

 

 

 

 

 

 

 

 

 

 

 

 

a’i — с — AG

 

 

 

 

 

а после удовлетворения условия при г = b

 

 

 

 

 

 

 

 

 

 

,3

3

 

 

 

 

 

 

 

 

 

= -

AGc' Ц ^ - .

 

 

 

 

(18.19)

 

 

 

 

 

Ь3г3

 

 

 

 

 

В

области

пластических

деформаций

<

г «S г,),

согласно

(18.17), (18.18),

после

удовлетворения условия при г =

а,

можем

записать

 

 

 

 

 

 

г 3 аЯ

 

 

 

 

о'х — р +

Y <6ка‘ In — -\- AGc' (1 — к)

 

(18.20)

 

 

л 8

 

 

Постоянные

с' и г, определим

из условия,

что

при

г = г,

э — э, и что напряжения (18.19) и

(18.20)

должны совпадать. В

итоге

при

a «S г ^ г,

получим

 

 

 

 

 

 

 

 

о \’ = — Р + 2к о, In + т ( I — к) a, L z iL ± t

 

 

 

 

 

+ kot + ( 1 к) аг - f,

 

 

 

 

(18.21)

 

 

 

 

 

Г

 

 

 

 

 

При г, ^ г ^

Ъ

 

 

 

 

 

 

« ; —

( l -

£ ) й щ

= 1 с, (1 + £ )

£

(18.22).

При этом радиус сферы г„ отделяющий

области

упругих и

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

 

 

U , -

1 - ? ■

+ « г ■+( 1 -

*> ( ? ■- ■>)•

 

Я 8-23»'

В приведенных выше формулах

аа=

'yf

о8, вв =

j / * э*.

Если пластические деформации распространились по всей тол­

щине шара, то

 

 

 

 

 

 

 

<*i = - p

+ 2кояIn £ + (р -

2кояIn |

 

 

 

о'2 = о [ +

ко$+

2kosIn

 

g

 

(18.24)

8* = m h k ) ( P ~

2ha>ln T )

 

ei = “

2s*-

 

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

о

'

♦ о

'

* о

' *

о

/ *

o j =

о г , do =

(J‘2 —

сг2 ? ei =

e i — 8i 1

8*2 =

е2 — е2,

где ои а2, еь е2 — напряжения и деформации перед началом раз­ грузки, определяемые формулами (18.21) и (18.22) при условии сохранения области упругих деформаций и формулами (18.24), если эта область исчезла; вторые слагаемые определяются форму­ лами (18.22) при г0 ^ г ^ Ъ и формулами (18.21) при а ^ г ^ г0> в которых следует о, заменить на 2аа и г, на г0. Отсюда при а ^ г ^ г9

 

 

 

 

 

гз

ol

=

— Р + 4колIn £

+

го

(1 — к) о5

о*

=

CTi + 2 /сОз +

2 ( 1

к) oirl/r31

еГ =

— 2etro/r3,

е* =

e5rj/r*;

3 в. В. Москвитян

при Г0 < Г ^ Ъ

(Ti = — j O s

1 +

* !

1

 

+

Ь*

J г3 •

Радиус сферической поверхности г0, отделяющей область вто­ ричных пластических деформаций, согласно (18.23), удовлетворя­ ет уравнению

3 р Ъ3- г 3 _ . г0 . ..

_ . г3 - а 3

4ас

Отсюда, в частности, можно определить условие появления вторичных пластических деформаций, положив г0 = а,

р ^ -j os (fe3 а3)lb3.

Ha рис. II.5 представлен график остаточных напряжении

с® при следующих значениях

параметров: Ь/а = 2,5, р = 2,227о«,

о

к = 0,95.

Штрихами

отмечены ос­

 

таточные

напряжения, вычислен­

 

ные по теореме об упругой раз­

 

грузке (без учета вторичных пла­

 

стических деформаций). Из при­

 

веденного графика следует, что в

 

случае

наличия областей вторич­

 

ных

пластических

деформаций

 

распределение остаточных напря­

 

жений существенно отличается от

 

распределения напряжений, под­

 

считанных по теореме об упругой

 

разгрузке. Такое различие харак­

 

терно для всех случаев, когда в

напряжения в толстостенной сфе­

процессе разгрузки элементов кон­

рической оболочке.

струкций появляются области вто­

 

ричных пластических деформаций.

3. Циклический изгиб круглой пластины. Круглая шарнирно опертая пластина толщины h нагружена равномерным давлением

<7 > 0 . Пусть М г — радиальный, а М 2— тангенциальный изгибаю­ щие моменты. Уравнение равновесия имеет вид

 

dMл

Мл— М9

1

f

/ v

,

 

 

 

 

 

 

 

(18.25)

Радиальная

и тангенциальная

%2

кривизны пластины вы-

ражаются через прогиб и/(г):

 

 

 

 

> ,

 

d2u>'

/

1

dw'

 

/

 

dr

* 2 =

7 1 7 *

*! =

** + г 2 ?.

 

 

 

 

При этом

м'г = D (1 - Q') +4- *;), м [ = D (1 - Q') ( « ;+ i- *1),

(18.27)

где D = GAV3; в случае линейного упрочнения £У = О при э ^

= Э |г|*Л/2 ^ *«*

fi' = fcfl — -| = + ^ s ] ПРИ ~э> э>-

Подставив эти выражения в (18.25) и учитывая (18.26), получим уравнение для определения х2[81]:

L

+ 4 < i - a ' ) ^ = s - (i8 -28>

dr

После опеределения х2 неизвестную кривизну хх найдем из (18.26), а изгибающие моменты по формулам (18.27). Прогиб w'(r) определится в результате интегрирования согласно (18.26).

Если при нагружении равномерным давлением q в пластине не появились пластические деформации, то соответствующий про­ гиб будет [81]

wе Ида*

( ,

14 Г 2

, З г4 ^

(18.29)

192Z)

^ 1

11а2 11а4 /*

 

где а — радиус граничной окружности пластины.

 

Пластические деформации (в

точке

г = 0, z = ±h/2)

появятся

в случае, если давление достигнет величины

 

7.

 

16*4

 

(18.30)

 

21а2

 

 

 

 

 

Пусть теперь q > qa. Для определения прогиба в этом случае можно воспользоваться, например, приближенным решением А. А. Ильюшина [81], согласно которому

w* (г) =s= Wo ( 1 -

 

3 г4

'

1 1 1 / 2

(18.31)

11

На*

=

п

 

М У З ™ 1

 

где т = э(0)/з,. Построен график величин т ~ п, где обозначена

гс = 21 q a 2/ ( 1 6 o t h z ).

Предположим, например, что пластические области распрост­ раняются вблизи граничных плоскостей по всей пластине, в этом случае щ = 3,5, а из указанного графика при & = 0,95 получаем тг = 2,175. Соответствующие зоны пластических деформаций на

8 *

рис. I I .6 заштрихованы. При этом

2

,2

=

д = 1 ,6 6 0 -^ , z' (0) = 0,286А/2, (18.32)

Л

CL

где 2 в(г) — уравнение границы, отделяющей области упругих и пластических деформаций.

После полной разгрузки пластина сохранит некоторый оста­ точный прогиб = w°(r). Для его определения воспользуемся сначала теоремой об упругой разгрузке (§ 13), согласно которой

=

 

w°(0) =

0 ,4 1 8 ^ .

(18.33)

При этом использованы формулы (18.29) и

(18.31). Заметим, что

и;°(0 ) составляет около 38%

от максимального прогиба

ш0, кото-

I

 

рый имела пластина перед нача­

Г *

 

лом разгрузки.

 

 

 

 

 

 

Если

теперь подсчитать

мо­

 

 

дуль а*,

то

окажется,

что

э* =

'^пкЛ Ш Й Ш Ш т т т ^

I

= 2,18эв> 2э3.

Это

означает,

что

в процессе разгрузки появляются

 

 

1а--

 

области

вторичных

пластических

Рис. II.6. Области упругих, плас­

деформаций.

Однако

вследствие

тических и вторичных пластиче­

того, что в нашем случае область

ских деформаций при изгибе и

вторичных

пластических

дефор­

разгрузке пластины.

 

маций

мала

(эта

область

на

 

 

рис. I I .6 заштрихована вертикальными линиями, при этом zj (0) = 0,917Л/2, rj = 0,313а) , фактический максимальный оста­ точный прогиб мало отличается от соответствующего прогиба, подсчитанного по формуле (18.33).

Рассмотрим теперь переменное нагружение пластины давле­ нием q" < 0 , под действием которого прогиб сначала будет умень­ шаться до нуля, а затем увеличиваться в направлении действия нагрузки q " . Для определения в этом случае прогиба w" (г) вос­ пользуемся теоремой о переменном нагружении, согласно кото­ рой

wn = w‘ — w* = w

14

г2

З гМ

11

аа +

(18.34)

где

Н а * Г

 

 

 

wo = н?о— wo.

 

(18.35)

аш 0 — максимальный прогиб в пластине, нагруженной из естественного состояния давлением q — q" при условии, что предел текучести равен а2а,.

Воспользуемся представлением прогиба в форме (18.34) п определим величину давления q l, которое необходимо приложить к пластине с тем, чтобы после его удаления остаточный прогиб пластины был равен нулю (правка пластины). Пусть перед на­

гружением

давлением q0 пластина имела остаточный

прогиб

U7° = w°(r),

сохранившийся после нагружения давлением

q' При

разгрузке из состояния, в котором находилась пластина под дей-

ствием искомого давления

/

имеет место теорема о переменном

q0l

нагружении

(при /г =

3 ),

согласно

которой

w = w" + w, где w"

определяется

формулой (18.34) при g" = д<ь

a w есть

прогиб в

упругой пластине при изгибе

ее давлением q — g0 .Поскольку тре­

буется, чтобы при q =

0 w = 0 , то

 

 

 

 

IV = —

W =

11?оа

14

_3г_

 

 

192D

 

11а4

 

Учитывая (18.34) и (18.35), получим

 

 

 

 

 

/

*

ll?oV

 

 

 

 

“ 'О

w* -

1920 •

 

 

Очевидно

 

 

 

 

 

 

 

 

 

W0 - w° (0) + J 2D ,

 

 

поэтому

 

 

 

 

 

 

 

 

(0 ) + ^

(g' — go) -

 

(18.36)

Величина

определяется решением задачи об упругопласти­

ческом изгибе пластины давлением д' — g0 при условии,

что пре­

дел текучести материала равен а 2ав. При этом уравнение (18.36) определяет единственное неизвестное — искомую величину давле­

ния д0. Используя числовые данные из нашего примера

 

а2а,

q’ =

Л2ав

 

(0) = 0,418—

1,660 —2^,

получим ql =

— l,240A2crs/a2, что составляет 75%

от qf При этом

принято а 2 =

2 .

 

 

 

Вьшишем, наконец, выражение

для

прогиба,

который будет

иметь пластина из циклически упрочняющегося материала после значительного числа циклических нагружений давлением ± д ', воспользовавшись теоремой о предельном состоянии (§ 14),

причем величина в^ = {21Ъ)1^ э 1

при г = О и

Ы = h /2

связана с

давлением

q' определенным

соотношением,

которое следует из ре­

шения задачи об упругопластическом изгибе пластины при пер­

вом нагружении.

 

 

 

 

 

 

 

 

В рассматриваемом нами примере

 

 

 

 

 

ш =

(0.518Р +

0,857) ^

( l -

+

^ 7

) .

(18-37)

Поскольку Р <

1 , ТО W <

w '

 

 

 

 

 

 

4.

Циклические температурные воздействия в упругопластичс-

ской пластине.

Пластина

толщиной

h произвольной

формы в

плане свободна от внешних нагрузок. На обеих граничных

плоскостях пластины осуществляется подвод тепла </(£), причем

будем считать, что д(0) =

0. Оси х2 и х3 прямоугольной декарто­

вой системы

координат

 

расположим

в срединной

плоскости.

Распределение

температуры

Т(х^ t) симметрично

относительно

срединной

плоскости Xi =

0.

Рассматривается область

пластины

на4 достаточном удалении

от

ее

краев. Все

константы

материала

будем считать независимыми от температуры.

Предположим, что все компоненты напряжений равны нулю, за исключением о22 = Озз — сг0. При этом из соотношений Оц ~ Вц следует, что отличными от нуля будут лишь компоненты дефор­ маций Вц и е22 = е33 s е0. Поскольку все компоненты напряже­ ний и деформаций можно рассматривать как функции только xt

и £, из условия совместности и условия симметрии

следует, что

8 #= BQ{t) [34]. В нашем случае

 

 

 

 

 

 

 

 

5ц —; 2а0/3,

S22 == S33 === Оо/З,

 

Эц

8 ц

8 ,

 

^22 ==s Э33 80

8 ,

Зе == 8 ц

2е0.

 

 

При этом из уравнения (16.1) находим

 

 

 

 

 

 

2о0(х, t) =

9К(г(х,

t) — аТ(х,

£)),

х =

х {.

(18.38)

В области упругих деформаций, согласно

(16.3)

и

(18.38),

о0,

t) = 6 G(e0(£) — e(x,

t)).

 

 

 

(18.39)

В области пластических деформаций из

(16.4)

и

(18.38)

следует

« . (* ,< > - ■ £ < • . « > - « < * ,о > .

 

 

 

 

 

 

(1840)

а — ko‘ + 2 G ( 1 к) $,

а = У 2 1 а0 \/У3.

 

Здесь принято линейное упрочнение материала. Из (18.40) на­ ходим

2|о0|< ? (!- * ) I °о |~ ка<

Исключим г(х, t) из (18.38) и (18.39)

 

 

 

а0=

(во № — а-Т (*. 0).

э < э*

(18.42)

— модуль Юнга,

v — коэффициент Пуассона) и из

(18.38) и

(18.40)

 

 

 

 

6C|an| (l - f e )(„

- « Г ) ,

э > э 8.

(18.43)

<*о

 

|% I ~ ка* \ ° 9К

 

 

 

Уравнения (18.42) и (18.43) определяют связь а0 ~ е0 соответ­ ственно в упругой и пластической областях пластины.

Предположим теперь, что тепловой поток q(t) достаточно мед­ ленно меняется со временем; в этом случае можно принять [34, 3451

 

Т (х ,г) = Щ

£ ,

 

 

(18.44)

 

х/Г

 

 

 

 

где и — коэффициент теплопроводности.

 

 

 

Осталось

использовать условие

отсутствия

результирующих

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

 

 

 

 

 

 

/ 1/2

 

 

 

 

 

 

I* о0 (х, t)dx =

0.

 

 

(18.45)

 

—/1/2

 

 

 

 

 

Введем безразмерные величины

 

 

 

 

 

 

оо

 

 

 

 

(18.46)

 

а5

 

 

 

 

 

 

 

 

 

 

Из анализа упругого решения можно заключить, что впервые

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

при х =

± h l2 , причем

здесь

о0 < 0 . При

этом соотношения (18.42)

и (18.43) преобразуются

в безразмерных величинах (18.46) с учетом (18.44):

 

а = - | - ( е ( « ) - е И * ) £ 2) ^ ; .

 

 

 

 

3? =* Z> [2(1 — * ) ( в (*) — e j ( 0 6*) — Л].

s > 5 „

(18.47)

г»

1 + v

 

0 _

a q

 

 

U

— 1 + v + 2 (1 к) (1 2v)’

ът ~

2xes’

 

 

гДе 1 — граница, отделяющая области упругих и пластиче­ ских деформаций.

Дальнейшие выкладки проделаем дли несжимаемого материа­ ла (v «=* 1 /2 ). При этом (18.47) преобразуется к виду

S = , 2 ( l - * ) ( ? ( f ) - e j ( f ) & * ) - * , 6 > ь .

Из условия сопряжения при %= |, находим

е (t) =

8 т ( 0

( * )

- 1/2

(к ф 0).

(18.49)

Удовлетворим условию

(18.45), используя (18.48),

 

1

— Б.) - i - d

- i . ) -

 

(18.50)

1 * ( 1

 

 

 

Исключив из (18.49) и (18.50) е, находим уравнение для опреде ления

___________1 __________

(18.51)

2 k g — (! — *) (1 -36*)’

Формулы (18.48)— (18.51) дают решение задачи о распреде­ лении напряжений и деформаций при температурном нагруже­ нии в случае, когда пластические области примыкают к гранич­ ным плоскостям пластины (до возникновения пластической зо­ ны в центре пластины). Из (18.51) находим функцию |aU), поел©

чего согласно (18.50) становится^ известной деформация в(О,

а из (18.48) находим напряжение а.

Предположим теперь, что тепловой поток монотонно возра­

стает от нуля (д(0 ) =

0 ) до значения qm = q ( t m)

и затем монотон­

но убывает до нуля

и что при t = tm еще не возникла централь­

ная пластическая зона. При этом,

начиная с

t = tmi в каждой

точке пластины

будет происходить

разгрузка.

Соответствующие

величины перед

началом разгрузки

обозначим

через е(т = £т(*т)%

о', &'Л$- Они находятся, как отмечалось, по формулам (18.48)— (18.51). Закон деформирования при разгрузке (16.2) приводит в нашем случае к соотношению

а(х, t) - о’ (х)

=

2 (e (f) -

в' - g* (eS-(f) - е°г')),

(18.52)

где сг, е, ег — текущие

величины. Удовлетворяя условию

(18.45),

получим

 

 

 

 

8 (f)

=

8 ' + 4 -

(е°г (<) - е°г).

(18.53)

Отсюда мы можем найти остаточные деформации ев, соответ­ ствующие нулевому значению теплового потока q. Посколь­

ку при q = 0 Вт = 0, из (18.53) находим остаточные дефор­ мации

б° = в1— вт/Зд

а из (18.52)— остаточные напряжения

а°(|) = а' ® + A eo,'(3 |2 _ 1)i

или с учетом обозначений (18.46) и (18.47)

о°0(х,) = о; ( *х) +

-

1 ) •

<18-54)

Мы привели аналитическое решение задачи о циклических температурных нагружениях пластины с линейным упрочнением (О < it < 1), включая определение остаточных напряжений. Как

О

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

Г Л А В A III

ТЕОРИЯ ТЕЧЕНИЯ И ЦИКЛИЧЕСКИЕ НАГРУЖ ЕН И Я

§19. Уравнения теории течения. Соотношения o*h~9k

Как уже отмечалось (§ 1), в пятимерном пространстве напря­ жений S5 предполагается существование поверхности нагруже­ ний Fioh) = 0, обладающей тем свойством, что при продолжении нагружения вне этой поверхности будет осуществляться активная деформация, а если траектория нагружения продолжается внутрь поверхности, то имеет место разгрузка. При постановке задач, приведенной в первых двух главах, понятие поверхности нагру­ жения было использовано лишь для того, чтобы установить ус­ ловие возникновения нового пластического состояния при после­ дующем нагружении после разгрузки.

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

В основе всех теорий течения лежит так называемый принцип градиентальности, согласно которому [85]

йэр + od(E-i) = dX grad F ,

где d(E-i) — матрица, характеризующая деформационную анизо­ тропию — изменение коэффициентов упругости вследствие пла­ стического деформирования (§ 1). Если деформационной ани­ зотропией пренебречь, то [299]

dap = dX grad F

(19.1)

или

(19.2)

Скалярный множитель dX может быть выражен через прираще­ ния пластических деформаций

dX

dabd3b j 1/z

Т_ °f °f

(19.3)

 

 

* i i dsi j

 

Как известно, принцип градиентальности можно получить как следствие постулата Друккера [299], который утверждает, что работа дополнительных напряжений Ао,^ = As,-, + б^До^/З на лю-

бом замкнутом цикле по напряжениям неотрицательна, т. е.

J hsijddij^ 0,

(19.4)

(а )

причем она равна нулю только в том случае, когда траектория нагружения целиком располагается внутри поверхности F = 0.

Учитывая (19.2), запишем систему уравнений состояния

,

dsij

, dF

dsa .j

(19.5)

ddij *2Q

h яс

 

 

dsi3

SK

 

 

 

 

 

Таким образом, согласно теории течения, зависимость между приращением деформаций и напряженным состоянием будет установлена, если в (19.5) выбрана функция F и указаны ее ар­ гументы. Приведем варианты теории течения и отметим возмож­ ность их применения для исследования процессов циклического нагружения.

В случае изотропного упрочнения принимается

F = F(o),

о = (sijSij)in.

(19.6)

При этом из (19.2) следует

ddva = h (a) sudX,

h (а) = 4 " 5 ? ’

(19-7)

причем функция упрочнения Мо) определяется, например, диаг­ раммой одноосного деформирования образца. Если, в частности, 2F — о2 — р2 = 0 (р = о*), мы приходим к известной теории Прандтля — Рейсса (# = <*>)

ds>

(19.7a)

d$ij = -ggrSijdK d&ii = 0.

Заметим, что при F = F(o) начальная поверхность

нагруже­

ния на активном участке траектории изменяет изотропно свои размеры так, что если при начальном нагружении условие пла­ стичности было о = о", то при разгрузке из состояния $ij, и при последующем нагружении по любому п ути пласти чески е де­

формации будут изменяться в том случае, если о " >

о ', где а" —

модуль вектора текущих напряжений nar =

Это

будет, в частности, и при изменений направления

нагружения.

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

В случае трансляционного упрочнения принимается

 

2 JF (si}) = <р (sij — sij) — р2 (к) = О,

(19.8)

где Sij — смещение поверхности нагружения,

р М

характеризует

изотропное изменение поверхности

текучести

в

зависимости

от

X на активных участках траектории

нагружения.

Очевидно,

что

тензор s?j определяется компонентами пластических деформаций

efj = sfj, причем,

вообще говоря, их

предысторией; при

этом

простейшей зависимостью будет

 

 

 

 

sfj = H0dfj,

Н 0=

const.

(19.9)

Для условия текучести Губера — Мизеса

 

2F = (si} -

sfj) (sn -

sfj) -

p* (X), p2 (0) = -§-"?«

(19Л0>

откуда, учитывая (19.2), получим

 

 

 

 

ddfj =

(s{j — s^j) dX.

(19.11)

Из (19.10) и (19.11) следует, что

 

 

 

 

р (X) dX = (dsR&Rj)4'.

(19.12)

Скалярный множитель dX в (19.2) при условии (19.8) может быть исключен, и, таким образом, приращение пластических деформаций будет связано с напряженным состоянием через функцию F. В самом деле, так как точка s{j + dsi} остается на поверхности нагружения, должно выполняться следующее со­ отношение:

{ds^ -

dsfj) j£- = 0 .

 

 

 

ij

 

Исключая отсюда и из (19.11) ds\j =

H Qddijt

получим искомое

соотношение

 

 

 

 

 

4

ь

V7

 

л ___ !_

 

 

aN ~

Н л

dF

0F

 

или при условии (19.10)

 

dskl dshl

 

 

 

 

 

# 0р2 (X) dX = ( si5-

s°ij) dSij.

(19.13)

В зависимости от вида соотношений Sij dij

и р = р (X) фор­

мулируются различные теории пластичности при трансляционном упрочнении.

Втеории пластичности, предложенной А. Ю. Ишлинским

[97]и независимо В. Прагером [208], предполагается, что s[j

связаны с э\) соотношением (19.9) и что р = р (0) = (2/3) ^ а, = а*. При этом из (19.11) следует

ddvij = (sij — H 03?j) dX.

(19.14)

Таким образом, в рассматриваемом случае не учитывается возможность изотропного изменения поверхности нагружения, т. е. имеет место так называемое кинематическое упрочнение. При этом, например, в случае чистого сдвига из (19.10) и (19.9) получим

а12 = ATS + 2G (1 — к) е12,

^ = ^ =

(19.15)

Это есть линейное соотношение между напряжением и деформа­ цией (линейное упрочнение).

Втеории пластичности, предложенной В. В. Новожиловым и

Ю.И. Кадашевичем [99, 100], предполагается, что

Р = Р М , sb = H (J *)9?h J° = s^ij,

(19.16)

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

реальных материалов. X. Циглер [349] предложил для

следую­

щее дифференциальное соотношение:

 

 

 

 

 

дР

 

Л

/

АЧ

ds..dsij

(19.17)

dsij =

( si} s°j) dp, dp =

------ - -----7 ТГ > 0.

ifhi ~ shl)

А. А. Вакуленко и P. А. Арутюнян [7, 8 ] с целью описания некоторых качественных особенностей последующих поверхностей нагружения (§ 8 ) предложили следующие соотношения:

р « р ( Х ) ,

ds^j = Я (J2) dafj,

J2 =

sijSi},

daP = (da^da?,)1' 1.

 

 

 

 

 

(19.18)

При этом для

установления

связи

~ э,7

следует

определить

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

функции Я(/а)

и рШ . При

Я= Н0= const из (19.18) следует представление (19.9).

Г.Бакхауз [288], используя принцип градиентальности (19.2)

и начальное условие пластичности Губера — Мизеса

 

= - § - ° г»»

(19.19)

предложил следующее интегральное соотношение:

 

S{j ~ s{j + si} | А (эр) da? + osJВ (эр) daft +

 

+ " « j ,C (a P )^ d a ? j . (19.20)

Сюда входят три экспериментально определяемые функции А(эр), В(эр), С(эр), для нахождения которых используются и ре­ зультаты опытов с циклическим изменением пластических дефор­ маций.

Заметим, что одним из возможных

соотношений s?j ~ э\j, оп­

ределяющих

закон смещения центра

поверхности текучести в

результате

пластического

деформирования с учетом

влияния

предыстории, может быть следующий функционал:

 

 

 

s

 

 

 

* « =

j

 

(19.21)

 

 

s T

 

 

где sT— длина дуги, начиная с которой появляются

пластиче­

ские деформации. Если ядро R(s, s{) представимо в виде R(s, st) =

=т о и з (19.21) следует

Sh - sb = Ri (s) R* (s) efi, (19.22)

где точкой обозначена производная по длине дуги s траектории

нагружения. В случае Ri(s) =

R 0= const из

(19.22)

следует

ds\3= R0R2 (S) <1э?5.

 

 

Если R(s, St) = Ho = const, то

зависимость

(19.21)

преобразуется

к соотношению (19.9).

Если же Rzisi) = = const, то из (19.21) находим s°ij = R°R1(s )9 l.

Как видим, представление (19.21) обладает достаточной общ­ ностью.

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

Пусть уравнение поверхности текучести в пространстве sfj при наличии температуры Т есть Fisa, Т) = 0. Перепишем урав­ нение (19.2) в виде

> Ь (М-23)

где точкой обозначена производная по времени. При этом [ 3 4 ]

• р

9F

тг

г\

&F

* .

dF rrt

л

— л

 

ъ

если/г = °»

d T sv +

d f T =

° ’

(19.24)

 

 

 

хз

 

 

 

=

0, если

F < 0

или

если

F =

0 и

s»j +

j f Т < 0.

В случае условия текучести Губера — Мизеса

 

 

Р (5iji Т) = ~2 ~svSij

3" о* (Г).

(19.25)

При этом (19.24) преобразуется к виду

 

3?j =

если

о2 =

Os (Г ),

/•

р > 0;

SijSij = 4 oto'tT,

 

 

 

 

 

 

 

(19.26)

= 0,

если

а2 <

af (Т)

или

если

 

 

о2 =

 

(Г)

и

 

/*

 

 

 

SijSij < - | - 0«а5!Г.

 

Здесь

 

 

 

 

 

 

 

 

ff

da‘ (T)

,

а - З $- э ? ±

 

 

G&

jr„

Р — oStJ9tj

 

 

 

 

dr

 

 

 

 

В приведенных соотношениях учтено, что изменение пласти­ ческих деформаций происходит в том случае, если точка, харак­ теризующая напряженное состояние, находится на поверхности текучести и если dF/dt = 0 при условии р > 0 (случай равенства соответствует нейтральному нагружению). Из последнего условия

с л е д у е т = 2Сэ^),что вместе с dF/dt = 0 приводит

к критерию нагружения SijS*j>-^-05(Tsr . При этом закон течения

(19.26) принимает вид

 

 

 

 

 

*

еСЛИ

2

$ (Г),

2

/•

Эij == P^ii»

0 “ == ~77~ 0

SijSij

о

*

= 0,

если

2

,

или

если

(«.2 7 >

о2С - ^ - а ; ( Т )

 

а* = 4

с; (Г ) и

s ^ j <

ог^Г

 

(случай равенства соответствует р = 0).

С целью записи соотношений (19.27) в единой форме иногда вводят функцию g(x, t), которая равна единице, если в данной точке (х, t) не происходит изменения пластических деформаций, и равна нулю в противном случае.

При этом (19.27) можно преобразовать с учетом соотношений

между коэффициентами упругости и связи (16.1) [34]:

 

£{j = ”4 ^

^ij 4 ahh + (1 ~ g) P (CTij —

4 )

(19.28)

Это соотношение удобно при численных реализациях.

Основные законы теории течения проверялись эксперименталь­ но. Наряду с удовлетворительными результатами здесь отмеча­ ются и отклонения, с чем можно познакомиться, например, в об­ стоятельном обзоре В. С. Ленского [136].

Приведенными выше уравнениями теории течения мы огра­ ничимся. В заключение рассмотрим вопрос об использовании соотношений теории течения в случае, когда траектория перемен­ ного нагружения в пространстве напряжений есть прямая (§ 3).

Запишем закон градиентальности (19.2) в пространстве 2 5

аэ1 = ^д -дХ ,

0.

(19.29)

dak

За функцию F(ck) выберем (19.10), что соответствует транс­ ляционному упрочнению в случае условия текучести Губера — Мизеса,

2 F =

(ак - о*) (ah-

o j) -

р2 (X),

 

р2 (0 ) =

- f

(19-30)

Пусть для определенности координаты а®

центра

 

поверхности

текучести связаны с э*

соотношением (19.9)

 

 

 

d

=

Нэ1,

 

Я

= # 0 =

const.

 

(19.31)

С учетом

(19.30)

и

(19.31)

соотношение

(19.29)

преобразуется

к виду

 

 

dal =

 

 

Hat) dk.

 

 

 

 

 

 

(о* -

 

 

(19.32)

Уравнение прямой в пространстве 2*

запишется в виде

 

O k = lc h,

6 = а*,

chck = 1,

ch = const.

(19.33)

Поскольку текущие напряжения а* связаны с о* соотношениями

ок = —- tf/iукомбинируя (19.32)

и (19.33), получим

 

ddl =

(ак %chНэ%) Adi,

Л = const,

(19.34)

где принято dk =

Ad%.

 

 

 

Соотношение (19.34) есть дифференциальное уравнение от­

носительно неизвестной э£(£);

это уравнение может быть проин­

тегрировано, поскольку величины ок и ch не зависят от

 

Произвольную постоянную с* определим из условия, что при о * = о, (6 = 6«) пластические слагаемые э* равны эк — a j (2 G)

пластическим слагаемым перед началом разгрузки. При этом

4 = (*• ~ 5 - ТГ - ZP + Ьтг) «ч> 1-НА« - EJ1 +

 

 

 

+ тг+ -г(ж-г)-

(19Ж)

Воспользуемся теперь зависимостями

 

 

Эh -— Эк

a k - ° k

Эк = эк— э*,

<**•

(19.36)

2G

 

 

 

 

Комбинируя (19.35) и (19.36) и учитывая (19.33), после неслож­ ных преобразований получим

(19.37)

где обозначено

а*

+1 = (-Jr — э') I 1 — СХР ( — НА (°* — а*))]-

(19.38)

D 2G ' Н '

Как отмечалось, предположение (19.31) соответствует линей­ ному упрочнению, поэтому при первом нагружении

о' = ко• + 2G (1 - к) э\

Н = 2G Ц р .

(19.39)

Отсюда

 

 

 

 

кэ*

 

D

1 -

к'

 

При этом формула для г|?2 преобразуется к виду

 

= Т ~1: t1 — ехр

НА(а* ~ О ) ] •

(19.40)

Как видим, (19.37) совпадает с соотношениями (3.24), кото­ рые мы получили при анализе переменного нагружения по

прямолипейньш

траекториям в пространстве 2*. Более

того,

для функций

(<**, и*) и Ф2 (а*>а*)

здесь

получены явные вы­

ражения (19.3$) и (19.40) с точностью до

константы А.

(Если

же, в отличие от приведенного выше,

считать

величину

А не

константой, а функцией Л(£), мы также придем

к (19.37), только

ij?i и % будут выражаться через о* и <JS с помощью квадратур.) 9 В. П- москвитиц

Вместе с тем здесь следует заметить, что соотношения (3.24), введенные для рассмотренной программы нагружения, являются более общими, чем (19.37) при условии (19.38) и (19.40) (или с учетом А = Л (| )), поскольку, во-первых, универсальные функции

( a*j °s) и ^ 2 (а*» ) в (3.24) определяются опытным путем, и они не обязательно описываются формулами (19.38) и (19.40) (или в более общем случае при Ж £ )) и, во-вторых, прп выводе уравнений (3.24) не использовались гипотезы теории течения —

закон градиеитальиости, соотношения о® ~

и поверхность те­

кучести Мизеса.

 

§ 20. Постановка задачи о циклическом нагружении

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

При нагружении из исходного состояния должны удовлетво­ ряться соотношения

+ F { = 0; Oulj = Ri на £ а, щ = uoi на Su; 2в„ = J + UJA

(20.1)

Принимая за поверхность нагружения (19.10), запишем

- Н'0э%) (4- - Я ; $ ) = Р2 (*'). р2 (0) = X а‘ - <20-2)

При этом имеет место (19.11), поэтому

ds'tj = ^ + (* « - Н М dk',

(20.3)

где скалярный множитель dk' может быть представлен в виде (19.12). В дальнейшем для определенности воспользуемся обоб­ щенным принципом Мазинга, вследствие чего р2 = а\оЦ§. Поэтому

 

Щ ё к ' =

(dB?;da?j )1/*.

(20 А)

При

< р2 (0) имеют место соотношения обобщенного зако­

на Гука

 

 

 

 

зц = 2G/v,

Gu = 3Ke'w

(20.5)

//

причем соотношение Оц ~ справедливо в процессе всего нагру­ жения.

Будем

считать,

что

записанная выше

задача

решена,

т. е.

 

 

t t

&ij нам известны. Пусть после разгрузки в процессе

величиныcTjj,

последующего

нагружения

силами

F\,

R{

в

теле возникают

напряжения Оц, деформации

ц перемещения U{ , удовлетворяю­

щие соотношениям

 

 

 

и

"

 

 

 

 

"

п

н

_

п

 

/

/

П" г» "

"

° i j j

+ Fi

= 0;

Oijlj

Ri на Sa,

щ =

uQi

на Su;

2еь =

uitj + u jti.

Поверхность нагружения по-прежнему имеет вид (20.2)

 

 

 

(*'и - н"оэ?;) (4 -

н ’’э?;) = р* (хя) =%.<£,

 

(20.6)

остаются справедливыми и соотношения (20.3),

(20.4) и

Оц~&ц-

 

 

d d i j

— -ттг +

(^ij —

 

) d \ " ,

о а — 3К е,ц ,

 

 

 

 

 

2G

 

 

 

 

 

 

 

 

 

(20.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Щ ч v

= ( « ы Ж

У '*

 

 

 

 

 

 

 

 

Новое пластическое состояние появится теперь при условии,

которое следует из (20.2):

 

 

 

 

 

 

 

 

 

 

(* « - я * # ) (* « - я ^ г ;) > р2 ( х ') =

 

 

 

(20.8)

а до его появления справедливы соотношения

 

 

 

 

 

s'a ~ s'lj - 2G {эц - эи),

 

ан ~

Ои = Ж

(е « - е «).

(20.9)

 

Как видим, упругопластические характеристики при втором

нагружении а3 (масштабный коэффициент) и

Н 0 (коэффициент

линейного

упрочнения)

отличны

от

соответствующих величин

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

Обобщая приведенное выше на любое число п нагружений, запишем соответствующие уравнения равновесия, граничные ус­ ловия и соотношения Коши:

оЪл + Fi - 0; а?;/; = Hi на 5“ и? = и?* па

2е% = u?j + uli. (20.10)

Поверхность нагружений описывается соотношением (19.10) при условии (19.9) с использованием обобщенного принципа Мазинга

(*« - ВДГ) (*Ъ - н 0э !? ) = р2 (хп)

п-И

«

(20 . 1 1 )

■б~а-

 

 

При этом

+ ( * Ъ - Щ ф ) а к я,

(2 0 . 1 2 )

а?, = зКгЪ, ^

dKn =

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

( s i - Я Г М Г 1) (*« - Я Г М Г 1) = ^ <*!

(20.13)

В зоне упругой разгрузки

 

sir1 - s i = 2G (air1 - э1), a ir1 - опн = ЪК (е?Г1- в?,).

(20.14)

Приведенные соотношения являются замкнутой системой для определения су” , e?j. Заметим, однако, что, вообще говоря, при ре­ шении конкретной задачи необходимо последовательное интегри­ рование системы уравнений, начиная с первого нагружения, по­ скольку соответствующие соотношения при данном нагружении содержат параметры предшествующего нагружения. Исключе­ ние составляет случай полной пластичности, когда при каждом

нагружении для определения су*7}, 8Гу достаточно соотношений (20.10) — (20.12), не содержащих величин с индексами п — 1.

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

Как отмечалось, параметры а„ и #о определяются экспери­ ментально. В частности, если из опытов найдены коэффициенты линейного упрочнения fcn, соответствующие величины могут быть определены по формуле (19.15)

H n0 = 2 G \ - k n

k n

§ 21. Теорема о переменном нагружении в рамках теории течения

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

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

Однако, как будет показано ниже, такая возможность пред­ ставляется, если: 1) несколько видоизменить уравнения теории течения для последующих нагружений; 2) в качестве поверхно­ сти пагружения использовать условие текучести Треска — СенВенана; 3) наложить некоторые ограничения на траектории на­ гружения. t

Пусть при исходном (первом) нагружении по произвольной программе силами Fi, Я* в рассматриваемом теле возникают напряжения Оц, деформацинв^- и перемещения и\. В области уп­ ругих деформаций справедлив обобщенный закон Гука

 

 

 

s\j =

2Gdih

<r« =

3ffe«.

(21.1)

За

условие

текучести

примем

условие Треска — Сен-Венана

тт ах =

Та, поэтому на границе области упругих и пластических

деформаций должно выполняться одно из трех соотношений:

 

 

I

<?21 = 2 Ts,

2 0Г31

= 2 Ts, |сгз — cFi |= 2 rs,

(2 1 .2 )

где о{,

о.!,

Оз—

главные напряжения, тв— предел текучести при

чистом сдвиге.

Как известно, условие

текучести Треска — Сеш-

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

напряжении Oi, а2, а3, что и делается ниже.

 

 

 

В случае условия текучести

Треска — Сен-Венана

за

уравне­

ние

поверхности

пагружения

при

трансляционпом

упрочнепии

следует принять

 

 

 

 

 

 

 

AF[ =

[do — аз — Я 0 (е£' — е£')]2 — a h ; = О,

 

 

4^2 =

Ws — o'i н о(еЗГ — е ?')]2 — а Ы = 0,

(21.3)

 

АУэ =

[oi- от; - Я 0 (в?' -

8?')]2 - « М

= О,

 

где

(ос = i , 2, 3) — пластические

слагаемые

главных

дефор­

маций.

 

 

 

 

 

 

 

При этом использованы гипотезы (19.9) и обобщенный прин­

цип Мазинга.

 

 

 

 

 

 

 

Как видим, в случае (21.3)

имеем три различных

аналитиче­

ских выражения для поверхности нагружения, поэтому следует составить три соотношения типа (19.1) для приращения пласти­

ческих деформаций

dF*

 

 

del' =

-? < & ;,

 

 

0,

а =

1,2,3.

 

(21.4)

 

 

даа

 

 

 

 

 

 

 

 

 

 

 

К соотношениям (21.4) следует присоединить

соотношения

для упругих слагаемых деформаций и 0ц ~

е«:

 

 

 

 

 

 

— 2G

ОЦ

vak A j

(Tii --

 

 

 

(21.5)

 

 

1 +

V

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача определения o>ij, e*j замыкается,

если

записать еще

не­

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

 

 

 

 

 

G i j J Ч" P i = 0 ,

2 S i j = U i j + U j fi ,

 

 

 

(21.0)

 

 

f

 

f

 

 

,

,

 

 

 

 

 

 

Ojjlj

на

 

Hi -— Mflj на

 

 

 

 

Пусть теперь после разгрузки

осуществляется

переменное

нагружение

п

силами

/'ь

в процессе которого возникают

на-

 

 

 

 

н

 

 

 

 

 

 

 

 

 

пряжения 0Г|; и деформации е^-. Введем разности

 

 

 

 

 

 

*

=

'

//

 

*

f

 

"

 

 

(21.7)

 

 

O i j

O i j -

 

 

8 ij =

8 i j

 

8 i j

 

 

 

 

*

 

 

 

 

 

главные

напряжения

<Ti,

и определим по o*j соответствующие

о*, Од и по б*; — главные деформации

е*^ е^, е3. Для пластиче­

ских слагаемых деформаций

составим уравнения теории

те­

чения в форме (21.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

del* =

-9^ d \ *,

 

 

0,

а =

1,2,3,

 

(2 1 . 8)

 

 

 

да.

 

 

 

 

 

 

 

 

 

 

 

и присоединим к ним формулы дляе^-

и а*

 

*

 

 

 

 

 

 

 

 

 

 

 

1

( ~ *

vakh^ij\

« * _ Q

|

 

 

(21.9)

 

 

eij 2G

 

1 +

v

) ’

**** ~ oKZjj.

 

 

 

С целью построения функций Fa

заметим,

что

для

простей­

ших напряженных состояний имеют место соотношения

 

 

 

 

OVr — O’er — On

 

*

'

"

 

 

 

(21. 10)

 

 

 

8а =

еа

 

 

где аа, еа — компоненты,

определенные

по

формулам

главных

напряжений и главпых деформаций через Oij и

соответственно.

Справедливы соотношения, аналогичные (21.3), например,

 

[в{ — о"г — Н 0(е?" — е?")]2 — а|т? = О,

которые вместе с соотношениями (21.3) при условии (21.10) для

некоторых путей нагружения сводятся к следующим:

[о^ — о* Н 0(е Г — е£*)]2 — (« 2

+

0С3)2Т* =

О,

[а* _ а3* - Н 0(е Г - е?*)]2 - К

+

а*)2 т* =

0,

(21.11)

[а3 — at — Н 0(е Г — е Г )]2 — (а2 +

а3)2 ^2 =

0.

 

На границе области, где возникает новое состояние пластичности,

8i* = = ез* = 0, поэтому согласно (21.11) при замене а3 на а2 имеем

Iа* — а * |=

2OC2T s, I я* — I =

2а2т5, I аз — о* |= 2а2т5.

 

 

(2 1 .1 2 )

К условию

пластичности (21.12)

мы могли бы прийти и дру­

гим путем. В самом деле, пусть, например, уравнения двух про­

тивоположных

граней

шестигранника Треска — Сен-Венана пе­

ред началом разгрузки имеют вид

/

/

л о1

' — ег )

= а2т5, а'г — о'« — Н 0(ef' — е ? ) -= — а2т5

Ci

— а2

 

 

 

 

(21.13)

Предположим, что траектория разгрузки и последующего переменного нагружения такова, что ее начало находится на первой грани, а точка с коордипатами а1? а2 находится на вто­ рой грани. В этом случае

с[ — Оо — Н 0(e j" — е2') = — a2xs.

Учитывая при этом уравнение первой грани и соотношение (21.7), получим

*

*

о

Ci

Со =

2а2т5,

в чем и следовало убедиться.

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

выполнялось условие

at ~ol—H0( e f — е Г ) = — <х3т8.

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

странстве главных напряжений на параллельную (по отношению к исходному нагружению) плоскость.

Так как функции Fa по своей природе являются универсаль­ ными (не зависящими от вида напряженного состояния), учиты­ вая (21.11), принимаем в общем случае

4Ft = [а* - вз - Я „ (е Г - e f* )]2 - (а, + а 3)2 xj = 0, (21.14)

В области разгрузки до появления нового пластического со­ стояния справедливы соотношения обобщенного закона Гука для

разностей (21.7)

 

S - J = 2Gat}, о « = зКеи,

(21.15)

причем последнее соотношение имеет место по всему объему де­ формируемого тела.

Выпишем, наконец, уравнения равновесия, граничные усло­ вия и соотношения Коши

+

F 'i - F 'i= 0;

2e*j =

u*iti + и*л;

o*jlj =

R{ на Sa,

иi =

(21.16)

uoi UQI на Su

Сравнивая теперь соотношения (21.8), (21.9), (21.14), (21.15) и (21.16) с соответствующими соотношениями предшествующего нагружения (21.4), (21.5), (21.3), (21.1) и (21.6), заключаем, что задача определения Cij, £* сводится к задаче нахождения напря­ жений и деформаций при нагружении рассматриваемого тела из исходного состояния силами F{ F { , К* — R\ при граничном перемещении и0\— и условии, что предел текучести материала т, заменен на а2т* и, кроме того, константа а2, если она входит в решение при первом нагружении, заменена на as. После опре­ деления сг*j и e*j искомые величины a*j, определятся с по­ мощью соотношений (21.7).

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

Сделаем несколько замечаний.

1.Из приведенного выше доказательства следует, что теорема

опеременном нагружении не имеет места, если области пласти­ ческих деформаций при переменном нагружении распространя­

ются на области упругих деформаций предшествующего нагру­ жения (см. также § 12).

2. Приведенное выше доказательство может быть распрост­ ранено на случай любого числа переменных нагружений совер­ шенно так же, как это сделано в § 12. Напряжения a?j и. дефор-

мации e?j определяются по формулам

 

оЬ = o'i} - S (- 1)к *h

е « = в Ь - 2 ( -

1)Ае: 1

 

fc=2

 

 

 

2

 

где

, Bij,. согласно

теореме,

совпадают с напряжениями и де­

формациями, возникающими при

нагружении

из

естественного

состояния силами

 

 

 

 

 

 

 

 

(

- i)h ( R t 1 -

я*)

 

с учетом граничного

перемещения (— l)ft («о Г 1 — «о») при уело-

вии, что предел текучести равен aftTe.

3. Выше использовалось условие пластичности Треска — СенВенана. Однако приведенные выводы будут справедливы и в слу­ чае условия пластичности максимального приведенного напряже­ ния, которое записывается в виде [78]

шах [||аг — сги/31|, |сг2 — сг«/31|, |ст3 — сг«/3||] = const. (21.16')

4. Мы рассматривали процесс нагружения в трехмерном пространстве главных напряжений. Пусть теперь компоненты 0{j образуют пространство напряжений. В случае кинематическо­ го упрочнения с изотропным изменением размеров поверхности нагружения уравнение поверхности будет

2F = Ф (ои — Oij) ц2(Я) = О,

где Л — параметр пагружения. При этом, согласно гипотезе градиелтальпости,

(21.17)

ij

В случае условия пластичности Треска — Сен-Венана соотно­ шения (21.2) определяют в пространстве Оц гиперповерхности, уравнения которых следуют из (21.2), если компоненты главных

напряжений оь <т2, о3 выразить через Оц. Это будут

уравнения

начальной поверхности нагружения.

 

 

 

Уравнение последующей поверхности

нагружения

мы полу­

чим, если в уравнении начальной поверхности

нагружения на­

пряжения Oij заменим разностями Да^ = а^ — а?,. Например,

 

[ф (Acrii) — <у2(До^-)]2 =

aU*.

 

(21.18.)

Функции нагружения будут

 

 

 

41\ =

[ст2 (До„) - а3 (Да;;)]2 -

а\х] =

О,

 

4 К

= [<т3'(Астъ) - ах (Да;;)12 - «1т? = 0,

(21.19)

=

[ ° i (A ffij) — о 2 (AcTij)]2 — a 2t j =

0.

 

Поясним это на примере плоского

напряженного

состояния.

В этом случае

 

 

 

 

 

 

 

Oof

±

4

((^и -

м

2+

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

Oi — а, =

((<хи — сг22)2 -| 4а\2) /г,

 

и из (21.18) получаем

 

 

 

 

 

 

 

(Д аа — Да22) 2 + 4Да22 =

а 2т2.

 

Соответствующее уравнение поверхности нагружения (21.19)

будет

 

 

 

 

 

 

 

4F3= (ои — Оп — о22 +

ст22) 2 +

4 (а12— о?2) 2 — а\х1 = 0.

Используя (21.17) и (21.19)

при условии <т?7- =

мы мо-

жем повторить доказательство приведенной выше теоремы.

§ 22. Приложения

 

 

 

 

 

 

 

1. Задача Нордгрена

и

Нахдп

[328].

Кольцевая

пластинка

с внутренним радиусом а и внешним

радиусом Ъ нагружается

по внутреннему контуру

независимо

изменяющимися внутрен­

ним давлением рШ и моментом M(t). В цитированной выше ра­ боте эта задача обстоятельно исследована с учетом разгрузок на базе уравнений теории течения.

Запишем уравнения равновесия

д о г

ст0 -

-|г(агег2) =-= 0,

(22.1)

дг

г

 

 

геометрические соотношения

диг

ит

e a - _ L f t §

“0|

Br=z~dF>

е0 = — .

Етв - 2 V дг —

г

и граничные условия

or {a,t) = — p(f),<yr (6,f) = 0, ue(b,t) = 0.

(22.2)

(22.3)

Уравнения теории течения запишем в Ф^рме (19.2), но не для девиаторных величин, а для текущих напряжений Оц и пласти­ ческих деформаций e fj:

i3

<22-4>

 

причем в пластических областях для матОР^алов без упрочнения (ограничимся этим случаем)

F= 0, i >О, F= 0

(22.5)

К (22.4) следует присоединить соотношения для упругих со­ ставляющих деформаций

еЪ' “ 2G

1 -f v )‘

(22.6)

Выберем функцию нагружения F, используя условие текуче­

сти Треска — Сен-Венана,

 

 

F ~ ~ Y max [ | — <тг |,

|<т2 — сга|, |ог3 — ^ | ] — та,

(22.7)

где rs— предел текучести при чистом сдвиге.

 

В пашем случае

 

 

1

1/2

(22.8)

= — (<Уг + <*о) dh

- г (о г — СТо)2 + OrQ

При этом для случая сгга0 < Оиз

можно принять

 

F3 т (°г — а«)2 + а'ге — т!.

Заметим, что из второго уравнения (22.1) следует

ОгвГ2= c(f)x5,

или учитывая, что Л/ = J сг,.or2d0 = 2лста, получим для

о

(Тг0 = А .

1пг

(22.9)

(22. 10)

а< г < Ь

(22.11)

Первые пластические деформации появятся в зоне, прилегаю­ щей к г = а; в пластической области < г < гв) функция = О, поэтому из (22.9) и (22.10) находим

t f e - f f r - - = ^ ( r 4- c 2)1/2,

(22.12)

Г

 

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

условию оо >

ог.

 

 

 

 

 

После интегрирования

первого уравнения (22.1)

с учетом

(22.12) и граничного условия (22.3) находим

 

ог + р

г 2 + (г 4 -

с 2)1

'2

,

(«4 - с У ’2

 

=

In я2 -J- (а4 —

с2) 1

/2

I

а2

(22.13)

В упругой области га< г *

От

= В

 

 

 

ОвГе}~(л+^Ь’ АЬ2

- - ( l + v ) ( i - i \ с

(22.14)

^ = ( l - v M + l± -VB. ^

 

 

 

\ г

Ъ)

 

При этом мы удовлетворим граничным условиям (22.3). Используем теперь условие непрерывности ог и условие теку­

чести (22.12) при г =

г,:

 

 

 

 

 

 

 

А1 - ^ = - £

+

14

r'i

+ ( ri

 

С2) 1/2

(а4 _ с2)1/2

 

a2 +

(о4

 

с2)Х/2 +

в*

8

В= (г,4 -

 

 

 

 

 

с-у:\

 

 

 

 

 

(22.15)

 

 

 

 

 

 

 

 

 

откуда с учетом соотношения АЬг = В находим

 

 

 

 

 

 

 

 

 

 

2\ 1/2

in +

+

? +

+

 

I

 

 

(22.16)

>2 +

(я4- с

2)U /2

 

 

 

 

Поскольку величина с определяется значениями момента М , уравнение (22.16) содержит единственную неизвестную вели­ чину г*.

Перемещения иг и ие в пластической области a < г ^ гв на­ ходятся путем интегрирования соотношений (22.4) при условии (22.9). Запишем их в виде

 

ег ------— "9*(сгг — сто),

eJ?o =

к&го, 6? — 0.

(22.17)

откуда находим

 

 

 

 

 

 

 

 

 

е? =

-в ?,

с? =

0.

 

 

f(22.18)

 

Учитывая (22.18) и используя (22.2) и (22.6), определим пере­

мещение Ur*

 

 

 

 

 

 

 

 

 

 

Еиг

 

 

Ее

 

 

(22.19)

 

 

— = (1 - ^ ) а г + ^ ,

 

 

 

 

 

 

 

где

величина

сМ ) находится из

условия

сопряжения

иг при

г =

гш:

 

 

 

 

 

 

 

 

 

Теперь нз

(22.17) находится

неопределенный

множитель X:

 

X = + [_ £ £ _

+

 

2гагз -

сс

С2)] 1/2

1

(22.20)

 

[ ( ^ - с 2) ( . 4-

 

 

е [ г 4- ,

 

 

Остается определить окружное перемещение uQ. С этой целью воспользуемся соотношением (22.2) для еге

_д_ / и Л

_ 2 е гв

 

дг \

г /

г 1

 

откуда находим

 

 

 

t

ь

Зег6 (Г- О drdt.

 

щ

 

(2 2 .2 1)

 

 

dt

 

При этом удовлетворяется условие (22.3).

Выражение для еге определяется с помощью формул (22.6),

(22.10),

(22.17):

 

 

 

- С Т s

, 1 +

V СТ5

 

8 г0 — А —

I------77---- Г »

 

г

^

г

где X дается формулой (22.20). Теперь правая часть (22.21) мо­ жет быть вычислена, и и0 становится известным.

Приведенное решение справедливо при условиях

1/2

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

Если теперь внешние усилия pit) и Mit) с некоторого момента времени начинают убывать, то в пластине будет происходить разгрузка, причем до появления вторичных пластических дефор­ маций решение Может быть получено с помощью теоремы об упругой разгрузке (§ 13). Если же один из внешних параметров pit), Mit) продолжает возрастать, а другой убывает, то следует определить условие начала разгрузки с последующим исследова­ нием решения, в том числе и при возникновении нового пласти­ ческого состояния [328].

2. Циклическое нагружение стержня с вырезом. При анализе напряженного и деформированного состояния с использованием уравнений теории течения получил распространение так назы­ ваемый метод дискретных элементов. Принимается линейное со­ отношение между приращениями пластических деформаций и напряжений, используется линейное матричное уравнение, опи­ сывающее упругое поведение исследуемых тел, в которых пласти­ ческие деформации интерпретируются как начальные. Не имея здесь возможности изложить существо метода и привести алго­ ритм для числейПой реализации, ограничимся приведением одного результата конкретного расчета [287].

17у80

1

^

—----------------- 7 Р

Рис. III.1. Геометрические

характеристики стержня с вырезом; р=

 

— 1,672 дюйма.

(5,нг с/см 1'7,03

Рис. 111.2. Зависимость местпых напряжений от деформаций у дна выреза для первого цикла нагружения.

Рис. III.3. Кривые нагрузка — деформация у дпа выреза для трех циклов.

Стержень

с вырезом, геометрические

параметры

которого

представлены

на рис. II 1.1,

испытывает

циклическое

растяже­

ние — сжатие

осевой силой

Р. Материал

стержня — алюминие­

вый сплав 2024-ТЗ. Использовался вариант теории

течения,

предложенный Циглером (см. (19.17)). Начальная

поверхность

текучести Губера — Мизеса

принималась смещенной

относитель­

ной нулевой точки пространства напряжений, и, таким образом, учитывалось различие в упрочнении сплава при растяжении и сжатии.

Результаты вычисления для одного цикла растяжение — сжа­ тие представлены на рис. III.2. Здесь показана зависимость меж­ ду местными напряжениями и деформациями у дна выреза. Раз­ личные кривые соответствуют указанным на рисунке максималь­ ным значениям номинального напряжения отах= Р т&т/Ро, где F0— площадь наименьшего сечения. Штриховая линия — кривая остаточных напряжений. Отмечается достаточно хорошее совпа­ дение с указанными на рисунке экспериментальными данными.

На рис. Ш .З приведены

результаты

для

трех циклов при

постоянной

амплитуде

поминального

напряжения

Стал=

= 351,5 кгс/см2.

 

 

 

 

Как видим,

соответствующие кривые

при

втором и

третьем

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

Г Л А В А IV

ВЯЗКОУПРУГИЕ И ВЯЗКОПЛАСТИЧЕСКИЕ СРЕДЫ ПРИ ЦИКЛИЧЕСКИХ Н АГРУЖ ЕН И ЯХ

§ 23. Вязкоупругие среды наследственного типа. Простые нагружения. Комплексные модули

Девиатор э{5 тензора деформации представляется в виде суммы

Эц = 9ij + э\$, причем слагаемые э% связаны с девиатором тензора напряжений склерономными соотношениями нелинейной теории упругости типа теории малых упругопластических де­ формаций

2G03ij = fi (о) o’ = (SijSij)

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

11

2GQ9*j = J Г (t, т) /2 (a) sijdx.

О

В итоге

t

2G03ij = Sij/i (о) - f | Г (t, т) /2 (a) s^dx.

(23.1)

О

Здесь G0— мгновенный модуль сдвига, Г(£, т) — ядро ползучести. Аналогичное соотношение предполагается справедливым и для ~ Оц:

ЪКфц = %1 ( |оц |) -f- j Гг (£, т) %2 ( I I)

(23.2)

о

 

где К 0— мгновенный модуль объемной деформации. Альтернативными к (23.1) и (23.2) будут следующие соотно­

шения, выражающие напряжения через деформации: t

=

»«Фх (э) — J*R (*. т) Фг (э)

(23.3)

0

О

 

 

 

t

 

 

 

( I eii I ) ^ ^1

*0 Фа ( |eii |)

(23.4)

о

 

 

Э = {ЭцЭц)и-.

Функции /л, фл, Хл» % предполагаются универсальными, не за­ висящими от вида напряженного состояния.

Соотношения типа (23.1) рассматривали Н. X. Арутюнян 15J, М. А. Колтунов [111, 112], X. Лидерман [314], Ю. Н. Работпов [215, 216], М. И. Розовский [231, 232] и др. Эти соотношения являются следствием квазилинейной теории вязкоупругости А. А. Ильюшина и П. М. Огибалова [91J.

Можно показать [166], что в случае подобия кривых ползу­

чести универсальные функции /До)

и /2(о)

совпадают /До) =

= /Д о)=/(о); при этом соотношения

(23.1) разрешаются относи­

тельно напряжений

 

 

=

x )3ijdx,

(23.5)

причем правая часть становится линейной функцией компонент Эц. Здесь R(t, т) — резольвента ядра Г(£, т). Аналогичным обра­ зом можно поступить и в том случае, когда определяемые из опы­ та универсальные функции связаны соотношением /До) — А/До), где А = const.

В уравнениях (23.1) может быть выделена линейная часть, ес­

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

функции

соДо) и соДо) соот­

ношениями

 

 

 

/До) = 1 +

со До),

/2(о) = 1 +

0)До).

При этом

 

 

 

2G Q9{J =

Sij -f- j* F

( t , T ) Sijdx

-j- f i j ,

 

° t

 

(23.6)

fi5== sijCOi (о) +

Г (t, x) CO2(O) Sijdx.

Аналогичным образом можно записать и соотношения (23.2):

ЗК^Е\\ — ои “Ь J Г*! (£, т) Gadx )- хп?

 

= аийх ( |Он I) + j (t, т) й 2 ( |ои I ) aHdx,

(23.7)

^ 1 ( I

I ) = Xl ( I

I ) — 1 » ^ 2 ( I

I ) = %2 ( I ° ii I ) —

 

Если

функции (Di =

G)2= Xi = %2 =

0 при всех интересующих

значениях аргументов, из (23.6) и (23.7) получим соотношения

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

t

t

26r0a^j = Sij + j Г (t, т) Sijdx, ЗК0гц = oi{ +

j Гх (t, т) 6ndx. (23.8)

10 В. В. Москвитин

В приложениях часто принимается ядро Г 4= 0, что отражает отсутствие реологических свойств в объемных характеристиках.

Соотношения (23.8) при Г = Г(£ — т), r t = Г 4(^ — т) могут быть представлены в виде интегралов Стилтьеса, например,

 

 

t

 

 

 

 

 

 

9ii =

f П (t -

х) dSij (х),

Г (t) =

2G0

 

 

 

0

 

 

 

 

 

 

e «

=

n

 

 

don (т),

Гx (t) =

dU

I t )

) U 1( t - x )

3K0

(23.9)

 

 

0

 

 

 

 

 

 

или в общем

случае

(23.1)

теории

нелинейной

вязкоупругости

при условии h =

/2 =

/:

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9l} =

$ H ( t - x ) d ( j ( o ) 8 ii).

(23.10)

 

 

 

 

О

 

 

 

 

Совершенно аналогичным образом могут быть преобразованы и соотношения (23.3) и (23.4); например, в случае теории ли­ нейной вязкоупругости

t

 

 

t

 

 

^ r = 9ij — ^ R (t — x)9ijdx,

=

еи — | Rx (t — т) z{idx

(23.11)

или в виде интегралов Стилтьеса

 

 

 

 

 

S

 

 

 

 

4j =

J* К ( t - x ) dai},

R (t) =

d£ %

 

 

0

 

 

0

 

 

t

 

 

dK

 

oii =

!iK 1( t - x ) d B ii,

=

 

(23.12)

 

0

 

 

о

 

 

К (0) = 2G0l

К г (0) = ЗК0.

 

 

В приведенные выше соотношения между напряжениями и деформациями как в случае линейных, так и нелинейных физи­ ческих состояний входит операция однократного интегрирования. При переменных нагружениях в случае существенно немонотон­ ных изменений внешних нагрузок однократного интегрального оператора может оказаться недостаточно для описания всех осо­ бенностей реологических свойств среды в этих условиях. При этом могут быть использованы уравнения квазилинейной теории вязкоупругости [92] с интегральными операторами более высоко­

го порядка, например, t

2G09tJ =

 

Sij +

j г (*, т) Sijd x

+

 

 

 

 

 

 

 

 

О

 

 

 

 

 

 

+

*

1 \

т1 2.тз) / (Т,)) s (Ti, т2 Sij

(т3

 

йххйх^х3, (23.13)

I I j Гз

 

 

0

0

0

 

 

 

 

 

 

 

ГДе 5

(Ti, Т2 =

5 tj(Ti)5 tj(T2 .

 

 

 

 

 

Здесь дополнительно введена функция /(а) с тем, чтобы в

частном случае эти соотношения с в о д и л и с ь

к

выписанным выше.

Задача определения о«, е«, и,- замыкается, если к физическим

соотношениям

 

теории

вязкоупругости присоединить

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

 

 

 

 

Оц + F, =

0;

2e0 =

u,iJ+M j>i;

 

 

 

 

 

 

7

D

С

 

о

 

(23.14)

 

 

 

 

0 ,-j/j =

/i,

на о0, Ui = uoi на

о и.

 

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

 

Оц(х,

t) = %(t)Oij(x),

Uiix, t) =

y\(t)iii(x).

(23.15)

Воспользуемся соотношениями

(23.1)

при условии

/2(а) =

= A f i ( o ) :

 

 

i

 

 

 

 

 

 

 

 

 

 

 

2G0Sij =

Sjj/j (о) + A j

Г (t, x) /x (a) si;dx,

(23.16)

 

 

 

 

0

 

 

 

где A >

0 — известная константа.

 

 

 

Что

же касается

соотношения

а„ — е,„

будем считать, что

имеет место одно из следующих:

 

 

 

 

ЕЦ =

0,

 

 

 

 

 

 

 

 

 

t

 

 

 

 

ЗК0г н =

<х»Хх ( |о » |) +

A j

Г (t, т) xi ( |Оц |) o^dx,

(23.17)

 

 

 

 

о

 

 

 

иначе говоря, предполагается, что имеет место или условие не­ сжимаемости, или нелинейное соотношение наследственного типа

(23.2) при условии х2(10«П =

£ x id a«l) и

x ) = A T ( t 1 т)/В,

В > 0.

 

 

 

 

Предположим

теперь, что

универсальные

функции

/До) и

5Ci(la«l) удовлетворяют условиям

 

 

/i(l^la)

= l^la/1(o),

xidfcl l o « l ) « l£la0C i(Q ),

(23.18)

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

При условии (23.15) и (23.18) из (23.16) п (23.17) следует

t

П (t) = I («) Ш а +

A J Г (*, х) 5 (г) 15 Г dx,

(23.19)

 

О

 

2G07ij = /х(о )

Sij, о = (sijSij)1/2,

(23.20)

=

0,

(23.21)

3/LOBU

( |(Jit |)•

 

Здесь учтено, что o(x, |) = |£|a(x).

Пусть внешние силы F,-, R( изменяются в процессе нагруже­

ния пропорционально одному параметру |(t):

 

Ftix, t) = & t)F((x),

R ({x, t) = %{t)Rt{x).

(23.22)

При условии (23.15) и (23.22) уравпения и граничные усло­

вия (23.14) перепишутся в виде

 

 

^ + У , - 0 ;

2 i„ - ~ и ,,+ щ.й

(2 з м )

Gijlj = Ri на So, Ux= 0 на Stt.

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

Таким образом, при сформулированных выше условиях точ­ ное решение задачи теории нелинейной вязкоупругости пред­ ставляется в виде (23.15), где функция £(£) характеризует закон (23.22) изменения внешних сил со временем; функция цШ опре­ деляется интегральным соотношением (23.19) и становится из­ вестной, если заданы %{t) и ядро Г(£, т).

Величины Оц, гц есть решение задачи для склерономной не­ линейной упругой среды, уравнения состояния которой пред­ ставлены соотношениями (23.20) и одним из (23.21); при этом

внешние силы равны Ft и Ru граничное перемещение и, = 0.

Рассмотрим теперь отдельно

случаи

линейной вязкоупруго­

сти, для которой /i = %i = 1,

А = В =

1. В этом случае из

(23.19)— (23.21) следует

 

 

 

 

 

Л (0

== I

(t) +

\Г (t,

т) I (т) dx;

(23.24)

 

 

 

О

 

 

 

 

 

 

2G0dij

Sijy

]

(23.25)

 

Ец

0,

3 K QEU z= O JJ,

^

 

 

а с учетом (23.15)

 

 

 

 

 

 

Sij =

2 G * ( t ) 3 ih о ,ч = 3К t(t)e и,

(23.26)

где

 

 

 

 

 

 

 

 

 

 

 

 

(23.27)

причем £U) и riU) связапы соотношением (23.24). Отсюда следу­ ет, что напряжения a,j и деформации е0, возникающие в линей­ ном вязкоупругом теле, совпадают с соответствующими напряже­ ниями и деформациями в линейном упругом теле с модулем сдви­ га G* и модулем объемной деформации /£*, определяемыми фор­ мулами (23.27). В этом смысле в случае (23.15) не существует самостоятельной проблемы теории линейной вязкоупругости. Заметим, однако, что при доказательстве (23.15) предполагалось, как одна из возможностей, выполнение условия совпадения ядер r tU, т) = Г(£, т), что если и реализуется, то для ограниченного числа вязкоупругих материалов. Вместе с тем остается еще це­ лый класс несжимаемых материалов (К0= °°), для которых ус­ ловия (23.15) имеют место.

Пусть теперь внешние силы изменяются со временем по за­

кону симметричных циклов

 

£U) = £0ехр (Ш ).

(23.28)

Воспользуемся интегральным уравнением (23.24) и найдем соот­

ветствующую функцию т](£), поменяв в (23.24) пижний

предел

нуль на —«>. Функцию г}(t)

представим в виде

 

г\(t) =

т|0exp (mt — йр).

(23.29)

Подставив (23.28) и (23.29) в (23.24), получим

 

По ехр (Ш — icp) = £0 ехр (Ш ) + Е0 j Г (t — т) ехр (йот) dx.

(23.30)

 

— оо

 

Здесь принято, что Г(£, т) =

Г(£ — т). Поскольку

 

ехр (mt) = cos a>t + i sin cof,

 

выражение (23.30) преобразуется к следующему виду:

 

Т)о ехр (m t г'ф) = |0(1 + Гс — гГ.) ехр (mt),

(23.31)

где

 

оо

 

 

 

оо

 

 

 

Г,

= J Г (z) sin cozdz,

Гс =

 

I* Г (z) coscozdz.

(23.32)

 

о

 

 

 

о

 

 

 

Учитывая (23.28) и (23.29), перепишем выражение (23.31) в

следующем виде:

_

1______1 +

г с + 1Г8

 

 

 

К О

 

 

 

Т (0 — 1 + гс -

(1 +

гс а +

Г2 ’

 

что, в свою очередь, позволяет

представить

соотношения

(23.26)

в виде

 

 

 

 

 

 

 

 

 

si} = 2 (G' +

iG") эц,

аи =

3 (К' +

iK") еи.

(23.33)

Здесь введены следующие обозначения:

 

 

 

 

G' (со =

G0

1 +

Гс

G" (to

=

Go

 

Г,

 

Гс)2+ Г 2 ’

 

Гс)2 + Г2 ’

 

 

(1 +

 

 

(1 +

(23.34)

 

 

1+ г1с

 

 

 

 

 

К ' (ю) =

К 0

К!' (со

=

Ко

 

Tis

 

rlc)2+ r 2;

 

T i c f + r 2;

 

 

(i +

 

 

(i +

 

где Г|, и Г|Свыражаются через ядро

 

(? — т) формулами (23.32).

(В действительности, согласно выражению (23.27) для К *, в фор­ мулы для К ' и К " должны входить не Г 1з и Г^, а Г5 и Гс; это вызвано предположением о совпадении ГД?) и Г(?), формулы же (23.34) записаны в общем случае.) Если ГД?) = 0, то из (23.34)

следует, что К ' = К е, К "

= 0 .

 

Комбинации G' + iG"

и K' + iK "

называют комплексными

модулями; их введение

позволяет в

ряде случаев эффективно

решать задачи о периодических возмущениях в теории лтгаейиой вязкоупругости [41].

Формулы (23.34) упрощаются, если вместо ядра Г(?) ввести

ядро R(t),

выразив

|(?)

с

помощью интегрального уравнения

(23.24),

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 (?) = ц (?) —

f R (? — т) г| (т) dx.

(23.35)

Повторив все выкладки, приведенные выше, получим

 

С'(со) =

С0(1 -Я с ,

G"(co) =

G0fis,

(23.36)

 

К ’ (со

=

Ко (1 -

Д1С,

К" (со

=

K 0RU.

 

 

Здесь R u

и R l3 выражаются

формулами,

аналогичными (23.32),

через ядро ЯД?) интегрального соотношения

(23.4)

при ф, = ф2=

= 1. Приравнивая выражения (23.34) и (23.36), получим

 

1 - R c =

 

1 + гс

Rc =

 

 

(23.37)

 

(1 +

Гс)2+ Г 2 ’

(1 + гс)2 +

 

 

 

г

и, в свою очередь,

 

 

 

 

 

 

 

 

 

 

 

1 +

Г

=

1 ~

Д<!

 

г

=

 

R*

 

 

 

 

 

с

 

( 1 - R cf

+

R * '

1

( 1 - Л

е)* +

Л У

 

Аналогичными

соотношениями

связаны R 1с, R is и

Г1с, Г„. Заме­

тим еще, что нз (23.31) следует

 

 

 

 

 

 

•По sin Ф =

 

 

 

 

 

 

 

 

 

 

Г

Q "

£оГ 8)

п « COS ф =

g0 (1 +

1’е),

tg ф =

=

7 ? ’

 

 

 

 

 

 

 

 

 

 

 

 

 

(23.38)

Поэтому,

если

через s?;,

 

и о?*, е?{ обозначить величины s,j,

a,j и а«, е«, соответствующие

значениям

£0 и

г)о. то из (23.26) и

(23.27)

получим

 

 

 

 

 

 

 

 

 

 

 

s\i =

2 (G '2 + G"0- ) ' 2B\h

ои = 3(К'2+

К"’) 1' 2гйц.

(23.39)

§ 24. Обратная ползучесть. Переменные нагружения

 

Под

обратной ползучестью

 

будем

понимать деформирование

во времени твердого тела после полного снятия напряжений. Рас­

смотрим сначала частный пример — деформацию

чистого

сдвига.

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

(23.8)

получаем

 

t

 

 

 

 

 

 

 

2G0e12 = а12+ j

Г (t — т) <т12 (т) dx.

 

 

 

(24.1)

 

о

 

 

 

 

 

 

 

Если напряжение о12 задано условием o12(t) =

<JJ2A(£),

где h(t)

единичная функция

Хевисайда (h = 0

при

£ <0,

h = 1

при

t > 0), то

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

2G0ela(t) = a°V2J (t),

/ ( f ) » 1 + J Г (£) d\.

 

 

(24.2)

 

 

 

0

 

 

 

 

 

Пусть теперь в некоторый момент времени tt мгновенно пол­

ностью снято напряжение (Ji2, т. е. о12 =

а?2 — <*12 ^ (t

*i)-

При

этом из (24.1) для t >

ti следует

 

 

 

 

 

 

 

2G0e12 (t) = а?2 (J (t) - J {

t - tx)),

 

 

 

(24.3)

в частности,

 

 

 

 

 

 

 

 

2Сови (h ) =

ai2 {J Ui)

!)•

 

 

 

 

 

Таким образом, при мгновенном снятии напряжения ai2 деформа­ ция падает на величину ej2 = aj2/(2G) и затем при t > ti умень­ шается, стремясь асимптотически к нулю (рис. IV.1).

Выражение (24.3) может быть переписано в следующем виде:

е 12 (0 = е 12 (£) — е 12 (^*)» 2G0812 = O12J (£*),

(24.4)

=0,

где е12 (t) — деформация при обратной

ползучести,

( 0

деформация ползучести при напряжении

ст?2 , е*2 (**) деформа-

ция ползучести при напряжении

Рис. 1V.1. Обратная ползучесть при чистом сдвиге.

Отсюда следует, что при t > tt

0x2 с началом нагружения в точке t —*1 U *=0).

Если напряжение снимается неполностью, а падает до значе­

ния £xai2

( O ^ l i ^ l ) , то

анало­

гично соотношению

(24.3)

полу­

чим

 

 

 

2G08X2 (0 =

0?2 [ J (0

 

- ( 1 - E l ) / ( * - * ! ) ] .

( 2 4 . 5 )

 

2 G 08X2 ( 0 = Ех^12^ (О ИЛИ 8 х2 (* ) =

(* ),

(24.6)

т. е. если кривая ползучести 8x2 = вх2 ( 0 имеет асимптоту

8 х2 (оо),

то деформация ползучести после частичной разгрузки асимптоти­

чески стремится к значению ?1 8i2 ( 00)*

Заметим, что формулы (24.5) и (24.6) остаются справедливы­

ми не только при 0 < | 4< 1

, но и при любых значениях |4, в том

числе и отрицательных. Например, при

= —1

 

 

2

(t) =

о°12 [J (l) - 2

J ( t - t,)],

 

и асимптотическое значение деформации равно — £12 (°°)-

 

Пусть, наконец,

при t = t2> t{ напряжение

мгновенно

изме­

няется до значения £2а?2, причем параметр

имеет любое

зна­

чение. При этом

 

 

 

 

 

2G0e12 = <т° 2 [J (0 + (%г -

1) J ( t - /,) -h & -

lx) J ( t - L)}.

 

 

 

 

 

(24.7)

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

Ei<*i2

при t = U:

п

 

 

2<?0Ц

^ = / (0 +

- £ i - i ) / ( f - * 0 , So = 1, * > « » . (24.8)

а 12

i=l

Если напряжение изменяется по ступенчатому закону в пре­

делах от £aj2 до а °2 (0 < S < 1) через равные отрезки времени #0,

то из (24.8) следует

2G„ ^

= J (t) + (1 - I) 2 ( - l)1j (t - it0), t > nt0. (24.9)

012

i= l

Эта формула может быть использована при исследовании виб­ роползучести (§ 25) линейных вязкоупругих сред.

Воспользуемся теперь соотношением (23.1) и определим де­ формацию в случае обратной ползучести физически нелинейной вязкоупругой среды. При чистом сдвиге из (23.1) получим

t

2G09(t) =

a/i (a) + J Г (f — т) /2 (a) a (T)

dx,

 

 

0

 

 

 

 

 

* = V 2|e13|,

a = / 2 | a 13|. (24.10)

Если напряжение

o(t) прикладывается мгновенно, т. е. oU) =

= o°h(t), то из (24.10) следует

 

t

 

 

 

 

2G0 V

= /1 (°°) +

U Ю (J (0 - 1),

/

(t) = 1 + j г ( I) d\.

 

 

 

 

о

 

 

 

 

(24.11)

Пусть

теперь при

t = tL напряжение

а полностью снимается.

Соответствующая деформация обратпой ползучести определится формулой

2G0 ^ r

= и

 

 

(24.12)

а в случае изменения напряжения при t =

t{ до значения £о° —

формулой

 

 

 

 

2G0o " /а0=

-

/2(£о0)) + /2(о°)Д*) +

 

 

 

+

Ш Ь о ° ) - h(o0m t - t j.

(24.13)

Согласно

(24.12)

деформация обратной

ползучести,

так же

как и в случае линейной вязкоупругости, асимптотически стре­

мится к нулю, а в случае частичного снятия

напряжения при

t > ti из (24.13) имеем

 

 

2G09"/O° = 6/i(go°) + I h i W m -

1).

(24.14)

Эта формула совпадает с (24.11), если в последней величину о0 заменить на £о0, т. е. при £ > ^ деформация ползучести совпадает с той, которую имел бы образец, находясь под напряжением |о°, начиная с £ = 0; иначе говоря, материал как бы «забывает», что он некоторое время 0 ^ ^ ^ деформировался при напряжении о0.

Мы рассмотрели деформацию обратной ползучести для случая чистого сдвига. Пусть теперь имеется твердое тело произвольной геометрии, на которое действуют внешние силы Ft и являю­ щиеся произвольными функциями координат точек тела, но из­ меняющиеся со временем пропорционально одному общему пара­ метру £(0. Будем считать тело несжимаемым, а функцию /До) удовлетворяющей условию (23.18). В этом случае напряжения и перемещения представляются в форме (23.15), а функции £U) и т](£) связаны интегральным соотношением (23.19) для случая фи­ зической нелинейности и соотношением (23.24) для линейного вязкоупругого материала.

Если все внешние силы, достигнув при t = 0 величин F\ = = l 0Fi (х), Ri = £>oRi(x), в дальнейшем остаются неизменными, т. е. если |U) = |0MO, то компоненты тензора деформации в каждой точке тела будут со временем изменяться в соответствии с формулой

М *, 0 = М * ) Ц +а/ 0(0, £о>0,

(24.15)

/0(<) = 1 + л ] ‘ г (Е )^ .

Здесь E i j ( x ) и соответствующие напряжения О ц ( х ) связаны скле­ рономными соотношениями (23.20), причем должны быть удов­ летворены уравнения равновесия и граничные условия (23.23).

Если при t = tt удалить все внешние силы, то напряжения станут равными нулю, в точках тела будет происходить обратная ползучесть, при этом

e"j (х, t) = е~ •(х) Ц+а (/„ (<) - /0 (f - h)).

(24.16)

Очевидно, что деформации Eij (#, t) со временем

асимптотически

стремятся к нулю.

 

 

 

В случае частичного снятия внешних сил (или увеличения до

значений, соответствующих параметру £ =

^

0)

е « (*, 0 =

(х) (Jo (0 + (£1+a - 1) J0(f -

 

f,)) io1+a. (24.17)

Если внешние силы испытывают скачки в моменты времени ih

до значений параметра £ft£0, соответствующие деформации будут

определяться формулами

 

г и (х, t) = ei; (х) l i +a Jo (t) + £ (й +а - й й ) J 0(t - th) I,

ft=l

J

S i = l , t > t n .

(24.18)

В случае линейной среды здесь следует положить

а = 0 и

JQ(t) = J (t ), поскольку при этом А = 1.

 

§ 25. Виброползучесть

Было замечено (см., например, [ 12, 182, 2441), что если к напряжению о0= const, под действием которого осуществляется ползучесть е0Ш, присоединить напряжение Дo(f), изменяющееся по гармоническому закону с высокой частотой со,

До = До0sin cof,

(25.1)

то при этом деформация ползучести e(f) станет отличной от e0(f) при сколь угодно малых амплитудах До0по сравнению с с0.

Прежде всего заметим, что если ввести разность напряжений

Oi — а2 = 2Да0

(о! > а2)

(25.2)

и определить деформации ползучести et(f) u e2(f) при напряже­ ниях соответственно щ = const и о 2= const, то в случае линейной вязкоупругости окажется

ei(f) =

 

e2(t) =

o2^(f),

(25.3)

где = тр(^)— функция (24.2),

характеризующая

вязкоупругие

свойства материала

 

 

 

 

\\) (t) = D yl +

j Г (£) d lj,

D = const.

(25.4)

Пусть теперь напряжение о(£) изменяется со временем по за­

кону

 

 

 

(25.5)

a(t) =

о0+

До0sin о)£,

причем напряжение а0= (щ + а2)/2, а амплитуда Да0 удовлетво­ ряет условию (25.2). При этом, используя соотношения линейной вязкоупругости (23.8), получим

е (г) = е0(г) + D —

(25.6)

е0(0 = - -t — ^(0 -

Пусть ядро ползучести имеет вид

Г(|) = Г0е х р (- ^ ).

При этом в случае со > у и t -*■ °о имеют место следующие оценки:

3 , 1

.

1 , 3

шахе = ~г 1 + "4’ е2 >

mine =

е2,

т. е.

 

(25.7)

е2^

в ^ 6i,

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

Иными словами, если напряжение оШ изменяется со време­ нем по закону (25.5) и о2< о < щ, то соответствующая деформа­ ция, согласно (25.7), не выходит из интервала ei — е2. Таким об­ разом, явление виброползучести не может быть объяснено в рам­ ках линейной теории вязкоупругости.

Эффект увеличения скорости ползучести при наложении на статическую нагрузку малой циклической составляющей для ря­ да полимерных материалов может быть убедительно объяснен разогревом материала вследствие рассеяния энергии (гл. V I).

Соответствующие исследования выполнены в работах [12, 182].

Зависимость скорости деформации ползучести р от темпера­ туры Т принималась в соответствии с известным кинетическим уравнением

р = р ( р ) ехр — ^ кТ Y— •

 

(25-8>

Здесь р — деформация ползучести, U — энергия

активации,

о — напряжение, к — постоянная Больцмана,

— константа ма­

териала. Из (25.8), считая функцию F(p) известной, находим

откуда

^5.9)

Формула (25.9) дает возможность определить деформацию ползу­ чести p(t), если известны функции а(£) и T(t).

Ha рис. IV.2 приведены экспериментальные и расчетные (по формуле (25.9)) кривые виброползучести капролона. Температу­ ра Г в (25.9) определялась как решение соответствующего урав­

нения теплопроводности при

условии,

что

диссипация энергии

дается формулами (35.2) и

(35.3), а коэффициент г|э определяет­

ся первой формулой (35.4)

при а = 0.

Как

видим, расчетные и

экспериментальные кривые различаются незначительно.

Оценка циклической ползучести жаропрочных никелевых

сплавов ЭИ 867 и Э П 109 приведена в работе [98]. Напряжение

а со временем изменялось в пределах упругости по закону (25.5),

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

Л = Ло0/о0 варьировалось в преде­

лах от 0 до °о. На рис. IV.3

представлепа серия кривых ползуче­

сти при различных значениях параметра А\ температура образ­ цов составляла 950°С. Как видим, имеется участок ускоренной

Рис. IV.2. Кривые виброползучести капролана: 1 — эксперимент, 2 — расчет, 3 — изменение температуры [182].

Рис. 1V.3. Кривые ползучести сплава ЭИ 867 при различных значениях па­ раметра А : о — 900°С, • — 950°С [98].

ползучести с интенсивным накоплением деформации. Условием интенсивного развития процесса циклической ползучести, по мне­ нию авторов, является

Т > (0,55 -

0,60) Т плавл»

^ -0 < Л р ,

 

где ТПлавл — температура

плавления, A llv — некоторое

(зависящее

от температуры испытания) критическое значение.

 

§ 26. Гипотезы теории ползучести и переменные нагружения

Рассмотрим вопрос об использовании некоторых

вариантов

известных теорий ползучести — теорий

старения, течения и уп­

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

Теория старения. Согласно этой теории предполагается, что, например, в случае чистого сдвига деформация г 12 есть функция

напряжения Oi2U) и времени t

 

e12 = 0(ai2, t),

(26.1)

причем вид этой функции определяется кривыми ползучести при постоянном напряжении aj2

81 2(£) = Ф (о 12 , £)•

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

Однако можно попытаться несколько видоизменить уравнения

(26.1) с целью использования их

и в случае изменения направ­

ления нагружения. В самом деле,

пусть

при

t = ti

мгновенно

снимается напряжение ai2(£i) > 0

и

прикладывается напряжение

п

 

 

 

 

 

 

 

,,

и пусть из опыта определено соответствующее семейство

кривых ползучести 8 1 2 ( 0

(£ >

0 )

для

различных

значений

a12 = const. Введем разности (£* = t t{ >

0)

 

 

 

е 12 (О)

e 12 ( 0 = е 12 (£*)»

^ 12

(0 )

&12( 0

= ^ 12

(£*)’

(26.2)

где, в отличие

от (24.4),

деформация е,2

берется в точке

t = tt.

Будем считать, что

е?2 = Ф * (<т?2) е12(^ ), t*), t* = t — £x > 0 , (2 6 .3 )

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

«Гг (**) = Ф* (о*2°, е12 (fx), i*). а*2 = а12(^ ) — о22 = const. (26.4)

В приведенных соотношениях о12— напряжения при обратном нагружении, е12 — соответствующая деформация. В качестве ар­ гумента функции введена величина деформации e i2U i), поскольку возможно, вообще говоря, что соотношение е12 ~ cr12 ~ £* будет зависеть от того, при каком значении е,2(^) построено соответ­ ствующее семейство кривых ползучести.

Зависимость (26.3) описывает и процесс обратной ползучести,

при этом or*® =<JI2(£I). Она может быть использована для исследо-

п

вания процессов переменного нагружения, если напряжение ° 12после изменения направления нагружения слабо меняется во времени.

В теориях ползучести часто полная деформация е12 представ­ ляется в виде суммы мгновенной деформации ej2, связанной линейным или нелинейным соотношением с напряжением 0|2, и деформации ползучести /?12(а12, t)

е12 =

е12 (^12 ) "Ь Pl2 (^12» Oi

(26.5)

причем в случае подобия кривых ползучести

 

е12 =

е12 (^12 ) ф (^12 ) ^ (0 е

(26.6)

Аналогичным образом может быть записано и соотношение (26.3)

е12 = е12 (^12 ) ~\r Pl2 (^12»

(26.7)

Приведем теперь уравнения теории старения для случая сложного напряженного состояния. Компоненты тензора дефор­ мации представим в виде суммы мгновенных деформаций и де­ формаций ползучести

= eij + Piji

(26.8)

причем Sjj и Oij связаны соотношениями, вообще говоря, нели­ нейной теории упругости, например,

о

_э

Э° = э° (а),

аа = 3Кг%

эа 5 " sih

о

О

 

(26.9)

6ifikk/3, эо =

(в ? ;*? / 7*

 

) —=Ci ;

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

Рц = i (оц ~ 8ijakh/3), р = (pijpij)11* = Ф(a, t), (26.10)

причем эти соотношения выписаны для случая, когда деформа­ ции ползучести рц не вызывают изменения объема (р,г = 0). Здесь, в соответствии с теорией старения для одноосного случая,

предполагается существование универсального

соотношения р ~

~ a ~ t, определяемого результатами

опытов

па

ползучесть для

простейших напряженных состояний;

например,

в случае .(26.6)

можно записать р = У2'ф(а/У2)£‘2Ш.

Если же отказаться от представления (26.8), то аналогично (26.9) и (26.10) полные деформации и напряжения будут связа­ ны конечными соотношениями деформационной теории пластич­ ности

Эц =

Sij, э = Ф (a, t), аи = ЗКен,

(26.11)

причем время t в эти соотношения входит как параметр.

 

В случае, если

найдена зависимость (26.1), можно

записать:

а = У2Ф(о/У2^ t).

Выпишем теперь соотношения теории старения для сложного напряженного состояния при переменных нагружениях. Введем

разности

 

 

 

 

 

 

Sij (^1 )

Sij (£) =

Sij (£*),

tfij (£]_) Gij (t) =■ Gij (f*)

(26.12)

и воспользуемся соотношениями, аналогичными (26.11),

 

9ij =

Э* =

/*

Э (^l)j

°ii = 3 *в „,

^0 ^

= Sij

6ijS^/3,

 

Sij =

G\j

^ijGhk!3.

 

Здесь универсальная

зависимость э * ~ а * ~ £ * определяется

соотношениями типа (26.3) при одноосном напряженном состоя­ нии, построенными по результатам опытов на ползучесть при мгновенной разгрузке и последующем ступенчатом переменном нагружении.

Если найдена экспериментальная зависимость (26.3), то

э* = У2Ф*(о*/У2, а'(^)/У2, £*).

Заметим, что уравнения (26.13), если в них исключить пара­ метр t*, совпадают с уравнениями переменного нагружения скле­ рономных упругопластических тел, с которыми мы встречались в гл. II.

Присоединим к (26.13) недостающие уравнения и граничные условия

a*jJ Fi ih) F\ ( у = 0,

2е*;- = ut.j +

(26.14)

o*jlj = R[ (tx) R\ (£*) на S.

Отметим, что приведенная выше постановка задачи приемле­ ма лишь в том случае, когда изменение направления нагружения происходит во всех точках деформируемого тела одновременно.

Теория течения. В теории течения предполагается, что ско­

рость деформации ползучести р12 является функцией от напряже­ ний Oi2U) и времени t

7З12 = /(а12| t).

(26.15)

Будем считать, что мгновенная деформация ej2 является, вообще говоря, нелинейной функцией о12(£). Поэтому для полной деформации справедливо соотношение

е » = ^ У ^ а 1 1 + /(а18>

t).

(26.16)

12

 

 

Функция / определяется экспериментально, например, путем обработки семейства кривых ползучести pi2 ~ t. Если кривые пол­ зучести pi2~ t подобны, то /(а12, t) = / i(a i2)/2(t).

В случае изменения при t = tl направления нагружения мы поступим здесь так же, как и выше. Для величин е12(£*)

ai2(£*)v введенных соотношениями (26.2), сформулируем гипо­ тезу (26.16):

Р*а = /* (<Ji2. О . еГг = е ?2 (о£>) + /* (аГ2)

=

(a^, f*);

 

 

(26.17)

здесь точкой отмечена производная по £*. Универсальная функ­

ция г|э* (ст12, t j)

определяется, например, путем обработки

экспе­

риментальных кривых ползучести е*2 = е*2 (a*2 , £*)>

=

const.

В случае сложного напряженного состояния соответствующие

уравнения теории течения запишем в следующем виде:

 

 

Pij =

«и. 3 = Ф (a,, О»- Я = (РиРа)Ч\

 

(26.18)

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

деформирования рн = 0 .

 

 

Универсальная зависимость q ~ о ~ t определяется

результа­

тами соответствующих экспериментов. Если,

например,

найдено

соотношение (26.15), то искомая зависимость

будет следующей:

q = V2 /(o/V2 , t).

 

 

И В. В. Москвитии

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

эи = Ш + 17 sij’ ч = ф (Ст- *). °И — ЗКги-

Если же не выделять компоненты тензора деформации ползу­ чести, а составить уравнения теории течения для полных компо­ нент тензора деформации e,j, то соответствующие уравнения при­ мут вид

Эц = е = <PJ (ст, t), Он = 3Кеи, е = (эцэц)1г. (26.19)

В этой форме уравнения теории течения для девиаторных ве­ личин совпадают с соответствующими уравнениями (7.6) теории

пластичности для траектории малой кривизны,

с которыми мы

уже встречались.

 

 

Универсальная зависимость е ~ о ~ t

может

быть найдена,

если известно, например, соотношение

(26.16),

в этом случае

е = У2 ф(а/У2 , t).

Воспользуемся соотношениями типа (26.18) и (26.19) и выпи­ шем уравнения состояния для процессов переменного нагруже­

ния. Уравнения, аналогичные (26.18), имеют вид

 

Pij = £ s*i, q* = 9* (о*, t*),

q* = {ptjptj) 1*,

(26.20)

причем, если использовать (26.17), то q* = У2 Ф*(о*/У2 , t tJ.

 

Здесь p*j = e*j e*f; величины o*j

и ег* введены соотноше­

ниями (26.12). Мгновенные деформации е*° связаны с ог?,-, вооб­ ще говоря, нелинейными соотношениями, например, уравнения­ ми деформационной теории пластичности при активном нагру­ жении.

Уравнения теории течения при переменном нагружении типа (26.19) запишутся в следующем виде:

 

е* — Ф* (ст*> t*)i

°и ~

ЗКе*{.

(26.21)

Здесь

 

 

 

 

V.

 

*

*

bij4k/3,

е*

(

 

 

= Eij

4j)

 

Если использовать

(26.17),

то соотношение

 

е* ~ а* ~ t*

запи­

шется в форме е* = У2 ф*(о*/У2 , ).

Теория упрочнения. Согласно теории упрочнения предпола­ гается, что существует зависимость

Как видим, в отличие от первых двух теорий, в это соотноше­ ние явно не входит время t, и только в случае установившейся

ползучести, когда pl2= Put, pl2 = const, из (26.22) формально сле­ дует уравнение теории течения (26.15).

Широкое распространение получило соотношение (26.22), за­ писанное в конкретной форме [215J,

Р12Р12 = / (ai2

а — const,

(26.23)

причем функция / часто выбирается степенной / = Ао\2. Из (26.23) для описания кривых ползучести получаем

Р12 =

т

_ i __

а?2 = const-

1 + а’

Уравнения теории упрочнения прн сложном напряженном состоянии запишем в следующей форме [215J:

.

S-J

Ри = 0

Pij = 9 1 ? =

 

 

(26.24)

( р =

( Р а Р ц ) 1*, д =

( р ц Р ц ) 1/\ a ^ { s i j S i j ) 12) .

Если известно, например, соотношение (26.22) при чистом сдвиге, то универсальная зависимость q ~ о ~ р запишется в виде

Ф(р/У2, д/У2, о/У2) = 0.

При изменении направления нагружения постулируется су­ ществование универсальной зависимости (26.22)

 

 

 

Ф* (р*2, Ри, о*2) =

0,

(26.25)

где,

как и раньше, р*2 =

е*2 — е*2, е*2

=

е*2° (о*2).

Величины сг*2 и

е12

введены выше соотношениями (26.2).

 

 

В случае сложного напряженного состояния выпишем следую­

щие уравнения состояния:

 

 

 

 

 

 

* ♦

а* *

q* ^ ф* (о*, р*),

 

ри -■= О

 

 

Р и

= 5*si.b

 

(26.26)

 

[р* =

 

 

q* =

{риРи)Ч\

 

 

( P U P */)

,

о* =

 

 

Для зависимости q* ~ о* ~

р *

может быть использовано соот­

ношение типа (26.23)

 

 

 

 

 

 

 

 

 

q*p*& =

Во*1;

В,

6 ,

7 = const.

 

Приведенное выше распространение известных гипотез тео­ рий ползучести на случай переменных нагружений может ока-

11*

заться полезным при исследовании циклических нагружений эле­ ментов конструкций (металлы при повышенных температурах). Более того, при некоторых условиях исследование переменного нагружения в этих случаях может быть сведено к анализу нагру­ жения из естественного состояния, как это будет ниже (§ 28) продемонстрировано на примере вязкопластических сред.

§ 27. Циклические нагружения вязкопластических наследственных тел

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

Уравнения (23.1), (23.2) при наличии температурного поля Tix, t) можно записать в следующем виде:

 

 

*

 

 

2(?о (^) Эг) ~

s i j f i (&i Т)

J

Г (t

т, Т) / 2 (сг, Т) S i j d x ,

 

 

0

 

(27.1)

зк 0 (Т) (г и -

ЗаА Г) =

а«;

а =

(si}si})1/*.

Альтернативными к первым шести уравнениям будут следую­ щие, выражающие напряжения через деформации:

2Q (т) =

(^> ^)

^) Ф2 (^> Т)

 

0

(27.2)

он = 3К0 (Т) (еа — ЗосЛТ);

э = (эиэ^)1^.

Здесь отмечено, что от температуры могут зависеть мгновен­ ные модули упругости Gо и К0, ядра ползучести Tit, Т) и релак­ сации Rit, Т ), а также функции fkio, Т ), срл(э, Т) (/с = 1, 2 ), ха­ рактеризующие нелинейные физические свойства среды.

Если же воспользоваться соотношениями (23.10), то при нали­ чии температурного поля можно записать

t

 

= .[ П (t — т, Т) d [f (а, Т) SlJ\.

(27.3)

О

 

Влияние температуры на физико-механические свойства мож­ но в ряде случаев учесть, как известно [92J, с помощью темпера­ турно-временной аналогии.

Если соответствующее модифицированное время обозначить через то

v

2G ^ j (t' )

=

sa (t' )

Д (a) + \ Г (*'—т') / 2 (a) s{j (x') d x ',

 

 

 

 

 

 

 

(27.4)

(0

=

2G0

(O <Pi(э) — j (*' — т') ф2(3

CO dx'j .

Здесь G0 пе зависит от T; модифицированное время t'

связано

с физическим временем t соотношением

 

 

 

 

 

t

 

 

 

 

 

 

*1

х'

 

(27.5)

 

 

 

 

 

М Д Е ))’

где ат(Т) — функция температурно-временного смещения.

При этом, как видим из соотношений (27.4), в случае мгно­

венного нагружения

(f

0 , t' -*■ 0 ) влияние температуры не учи­

тывается.

 

 

 

 

 

 

 

 

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

состояния присоединим к уравнениям состояния недостающие

соотношения и граничные условия:

 

 

 

 

 

 

ay, 1

+ ^

= 0;

2еч= atii + и

, -

,

(27.6)

 

 

a(jl} = Ri

на S„,

и( = uoi на

Su.

 

(27.7)

Рассмотрим вопрос о возможности представления решения

сформулированной выше задачи в следующем виде:

 

 

 

 

Оц(х, t) = t(t)on{x), и((х, t) = r\(t)Ui(x).

 

(27.8)

Предположим, что

 

 

 

 

 

 

1 )

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

 

 

 

 

 

 

Пх, t) =

r\(t)T{x);

 

 

(27.9)

2 )

материал механически несжимаем

 

 

 

 

 

 

 

е« = ЗаТ;

АТ »

Т;

 

 

(27.10)

3)

универсальные

функции фДэ, Г), ф2(э, Т)

и

мгновенный

модуль сдвига G0(T) удовлетворяют условиям:

 

 

 

ф ,(э, т) =

ф ш ./ т )\

ф2(э, т) = ф ш ,/ т )\

 

(27Л1)

ф(т]э) =

т^фЫ,

GO(TJT’) = n\pG0(T);

4 , б, р, р, Ts=

const.

4)ядро релаксации R(t, Т) представимо в виде

R(t, Т) = Д,(М 27Г.)*;

Т„ q = const.

(27.12)

При условии

(27.9)— (27.12) уравнения

(27.2)

преобразуются

к виду

 

 

 

ЪаТ,

 

Sij =2G0 (f)ф ( » ) ~эц( TjT)\

е «

=

(27.13)

1 ( * ) = Л Р (0 ^

г (0 - | я ( « - т ) т 1г( т ) ^ ;

 

Y = S — q,

(27.14)

 

г = 1 — у + р.

 

 

 

 

Предположим еще, что объемные силы F,

равны нулю, а по­

верхностные силы Ri и граничное перемещение uoi представимы

в виде № , t) = 'E>(t)Ri(x),

uoi(x, t) = i\(t)uoi(x). При этом, учиты­

вая (27.8), преобразуем (27.7), (27.6)

 

o,it j

0;

26ij

iiitj +

Ujt i,

Oijlj =

Ri

на Sa,

щ =

(27.15)

uoi на Su

Таким образом, при сформулированных выше условиях точ­ ное решение задачи (27.2), (27.7) представляется в виде (27.8), где функции времени £(£) и г\(t) связаны интегральным соотно­

шением (27.14). Величины су и Ец согласно (27.13) и (27.15) представляют собой решение задачи термопластичностн в случае

активного нагружения поверхностными силами при наличии

граничного перемещения uoi.

Заметим, что поскольку функции £U) и т\(t) между собой свя­ заны, процессы изменения внешних сил -Я* = |(Ш?<(#), граничных перемещений uoi = цШи0{(х) и температуры Т = цШТ{х) на гра­ нице тела оказываются в рассматриваемом случае связанными. Если же исследуемое тело свободно, то единственным (и незави­ симым) внешним параметром будет т](£).

Предположим теперь, что сформулированная в начале этого параграфа задача (27.7) и (27.1) (или (27.2)) решена и соответ­ ствующие величины, которые, как и раньше, мы обозначим

штрихом, нам известны. Пусть далее,

начиная со времени t = tly

во

всех

точках

деформируемого тела

осуществляется

разгрузка

и

последующее

переменное нагружение

объемными

силами

F{ (х, t) и поверхностными силами

Ri

{х, t)

при граничном

перемещении

u0i(x,t).

Возникающие

при

этом напряжения, де­

формации п

перемещения обозначим

через

{х, t), Bij {х, t) и

U\

(х, t).

Поскольку нас

интересуют реономные

упругопластиче­

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

введем разности

o*j = d'j (s, t) Gij (*, /), ei* = elj (ar, 0 — e", (*, 0

(f > *')

(27.16)

и заметим, что если нагружение, последующие разгрузка и пере­ менное нагружение в случае линейной вязкоупругости осущест­ вляются в соответствии с уравнениями (23.8), то, как легко убе­ диться, уравнения (23.8) при£* = £ — t ^ O эквивалентны сле­

дующим соотношениям:

t*

2G0s*j = s*j + J Г (f* — т*) s*j (т*) dx*.

(27.17)

О

 

В самом деле, запишем сначала соотношения (23.8) для некото­ рого монотонного процесса изменения напряжений s\j (t)

t

 

 

2 <V ij ( 0 =

4

 

( 0 +

j

г (£ — x) Sjj (x) dx.

 

 

 

 

 

 

 

 

 

0

 

 

 

 

Рассмотрим теперь

другой

процесс Sij (t)

такой,

что при t ^

^ i t1

Sij(t)^Sij(t),

а

при

t > t x

s'ij (t) Ф s[j (t).

При

этом также

справедлива зависимость (23.8)

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2G03ij (t) =

s"j (t) +

J Г (£ t) s"j (x) dx.

 

 

 

 

 

 

 

 

 

0

 

 

 

 

Вычитая одно из другого, получим при t > t t

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

2G0( 4 -

4 )

=

4

-

4

+

f Г (t -

т) {sij - s"j) dx.

 

 

 

 

 

 

 

 

 

*1

 

 

 

Если

теперь в

этом

соотношении сделать

замену аргументов

t =

+

£*, т =5 +

т*

и учесть обозначения

(27.16), мы придем

к (27.17), в чем и следовало убедиться.

Обобщением соотношений (27.17) на случай физически нели­

нейных сред

будут

следующие

уравнения типа

(27.1) при

t^ = t — t1^ 0 :

 

1*

 

 

 

 

 

 

 

 

 

 

2G0( Т) 4

- stjft (о*, Т) +

С Г* (£* -

т, Т) 1* (о*,

Т) s*5dx,

 

 

 

 

О

 

 

 

ЗЛГ0 {Т) [е*

— За (АТ (£г) — АТ (£*))] =

стГ*

(27.18)

или

 

 

 

 

 

 

 

*

 

 

'*

 

 

 

 

§5 ^ 7

) =

э*№ (э*, Т) — [

Я* (*»—■т, Т) ср£ (а*, Т) afjdx,

<4 =

ЗК0(Г ) [е*; -

За(ДГ (fx) -

А Т (£*))].

(27.19)

 

Здесь функции (а*, Т) и щ (э*, Т) (к = 1, 2) предполагают­ ся универсальными, не зависящими от вида напряженного со­

стояния. Соотношения

Оц ~

следуют из соответствующих со-

'

/

/

П

у

отношений Сц ~ £ц и ои ~

£ц

в соответствии с формулами

(27.16).

 

 

 

 

Уравнения (27.18) или (27.19) в случае отсутствия интеграль­ ных слагаемых совпадают с соответствующими уравнениями тео­ рии малых упругопластических деформаций при переменных нагружениях (§ 1 0 ).

С целью общей постановки задачи присоединим к указанным выше соотношениям недостающие уравнения и граничные ус­ ловия:

j "Ь (^i)

P i (^*) = 0; 2& ij

= U i j -f- U j t\]

^7 20)

Cijlj == P-i (^i) Pi (^sjc)

^a Siy,

U{ = UQI(^I)

U0i J|J) Ha

 

После определения

<Уу,еу

искомые

величины Су,

най­

дутся из соотношений (27.16).

Как видим, приведенная здесь постановка задачи отличается от той, что приведена в § 26: в зависимостях (27.16) &ij и сгц берутся в точке £, а не в точке как было принято в (26.2). Указанное отличие объясняется спецификой уравнений состоя­ ния наследственного типа, что было выше продемонстрировано на примере линейного уравнения (27.17). Это обстоятельство ис­ пользовалось Шойи и для характеристики склерономной упруго­ пластической среды [163, 338].

Обратим еще внимание на задачу определения си(х, t), £ц(х, t) при t > tt: эти деформации и напряжения соответствуют внешним параметрам ^ (£), и0* (£), монотонно изменяющимся на отрезке ( X f ^ ^ и остающимся неизменными и равными P i{tj)D P i(t i),u 0i(£i) при t < tu что и отмечено в (27.20). Если при t > > ti не происходит со временем перераспределения напряжений

Gij, то CTij = cTij (fx);

деформации 8|j( 0

находятся из соотношений

(27.1), в которых напряжение следует заменить па Oij(t±):

(fх) U {o', Тг) ^1 + f Г (t - т) dxj

+ 1

1

 

оQ

 

H - f r ( t - T )

si} (т) / 2 {o’, T)dx +

y -

6 цо'кк {tj + 2G0bijaATi,

0

 

 

 

7*1* Г (*i).

Если осуществляется произвольное п-е нагружение, при ко­

тором возникают напряжения cry и деформации бу, то введя раз­ ности

 

=

( ~ 1)” W T 1 - <*«),

=

( - 1)п

- e"j),

(27.21)

по аналогии с (27.18) запишем следующие соотношения:

 

 

 

 

*

 

 

 

 

2

=

s tn ln «

Т) + J

Г5 (t* -

т) й п (an*, Т) s*? (т) dx,

\

 

 

О

 

 

 

 

3 ^ [е*Г - ЗаДГ (*n-i) + ЗаАГ (£)] = at?,

 

(27.22)

* п =

* —

* n - i , On =

( s y 1s,*n) 1/2.

 

 

 

Здесь через £n-i обозначено время начала разгрузки и после­ дующего n-го нагружения.

При отсутствии интегральных слагаемых соотношения (27.22) совпадают с уравнениями теории малых упругопластических де­ формаций в случае гс-го переменного нагружения при наличии температурного поля (§ 1 0 ).

Запишем еще недостающие уравнения и граничные условия:

+

( - I f

(Р ’Г 1(tn-i) -

П ( О ) =

0 ;

2 еу* = и% +

и*Ъ

a t j h

=

( - I f

( И Г 1 (tn-i) - Щ ( О )

па

S a,

(27.23)

Щ п =

( - 1 ) " ("оГ1 (tn-i) -

uoi (0«) на Su.

 

Приведем теперь несколько примеров определения универ­ сальных функций, введенных в этом параграфе.

В работе [247] были обработаны результаты экспериментов, проведенных Л. Б. Гецовым [45]. Материал — сплав на никеле­ вой основе ЭИ 607 А; опыты — на растяжение образцов. При мгновенном нагружении (£ -*• 0) из (27.1) следует

h (tf,

Т) = 2G0 (Т)

э = (зуЗу)1/2.

(27.24)

Зависимость

G0 ~ Т

аппроксимировались соотношением GB—

= G°(TJT)P. Оказалось,

что

= 5,21 • 103 кгс/мм2, р ~ 1;

за Т,

было принято 923 *К. Результаты представлены на рис. IV.4. На рис. IV.5 приведен график функции /До, Т), построенный по фор­ муле (27.24). Зная Д(о, Т), функцию /2(о, Т) находим из (27.1) с использованием кривых ползучести

Э

(27.25)

/ .(*, Л “ » 4 п ( 2й. ('о О

где

#

ф = Г (£ — т, 71) dx -{- 1 .

о

1

Ядро Г(£, Т) определялось па участке физической линейности (/t = /2 = 1). График функции /2(о, Т ), построенный по формуле (27.25), приведен на рис. IV .6 .

ЗРр

нес

 

 

 

703

мм1

 

 

 

76-

 

 

 

 

19 -

 

 

 

 

72 -

 

I

i_____

 

70_____i_____i

 

0,62

0,86 0,90 0,99

0,98 TQ

 

 

 

 

Т

 

Рис. IV.4. Зависимость

мгновенного

Рис. IV.5. График функции

модуля

сдвига от

температуры

П(о,Т) [247]. ‘

 

[247].

 

 

 

fz(aj)

Рис. IV.6. График функции

Рис. IV.7.

Графики функций

/2(0, Т) [247].

/1 (0 , T),f2(o,T) [167].

Как видно из рис. IV.5 и

IV .6 , кривые

/ 4 и /2, соответствую­

щие различным температурам

, могут быть приближенно сов­

мещены с кривыми /i и /2 прп температуре Тв путем их смещения по оси а па отрезке ф(7\)(ф(Гв) = 0 ):

П (а, Т„) (а + Ф (7*)), А (а) m /, (<r, Tt) (i = 1, 2).

Этот результат при Т = Tit) дал возможность Л. X. Талыблы [247] свести задачу нелинейной термовязкоупругости (27.1), (27.7) с существенно зависящими от температуры свойствами ма­ териала к соответствующей задаче нелинейной термовязкоупру­ гости для среды, свойства которой не зависят от температуры.

Результаты опытов Е. Робинсона [215] с ползучестью образ­ цов из хромоникелемолибденовой стали были использованы в [167] для определения функций Д и /2 по формулам (27.24) и (27.25). Соответствующие кривые представлены на рис. IV.7.

Рис. 1V.8. Кривые ползучести при прямом а п обратном в кручешш [167].

Кривые ползучести при кручении тонкостенных трубок из алюминиевого сплава при температуре 150°С (В. С. Наместников [176]) приведены на рис. IV .8 (без мгновенных составляющих де­ формаций). Здесь же отмечены расчетные точки, причем прини­

малось /х = /2, fi = /а [167].

§ 28. Теорема о переменных нагружениях вязкопластических тел

Воспользуемся уравнениями состояния, введенными в преды­ дущем параграфе, и докажем теорему, которая позволит без ре­ шении дополнительной задачи определить напряженное и дефор­ мируемое состояние, возникающее в вязкопластическом теле при его переменном нагружении, если соответствующее решение при нагружении из естественного состояния известно. Иначе говоря, будет доказан аналог теоремы, приведенной в § 1 2 .

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

аргумент Т универсальных функций Д, /2, Д , /2 ниже опущен. Пусть при t < ti на рассматриваемое тело действуют внешние

силы F i(x 1t)l

R i(x , t),

которые вызывают напряжения

и деформации

гц (х, t).

При этом должны удовлетворяться

уравнения равновесия, граничные условия н соотношения Коши

Gijt3 Fi = 0 ; Oijlj = i?i на S ; 2 &ij = u%%j 4- Wj.i* (28.1)

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

о' ^ о3

и напряжения и деформации связаны соотношениями

 

t

 

 

 

2G0ajj =

s'ij + J Г (£ — r) sjj (т) dt, 3K0(е'и — ЗаДГ') =

а'ц. (28.2)

 

о

 

 

 

В области физической нелинейности (о' > о*) имеют место за­

висимости (27.1)

i

 

 

 

 

 

 

 

2G03i} =

Sij/x (а') + J Г (* — т) /» (а') sij (т) dr,

 

 

 

6

 

(28.3)

 

3/f0 (ей — ЗаДГ') =

сг«,

 

причем

на границе

указанных

выше областей /До*) =

/2(о*) = 1 ,

поскольку уравнения состояния (28.3) и (28.2) должны переходить одни в другие непрерывным образом.

Экспериментально определяемая величина о* есть напряже­ ние, при превышении которого кривые ползучести становятся фи­ зически нелинейными или, что то же самое, все изохронные кри­ вые о ~ э при с > ов становятся нелинейными.

Предположим теперь, что при t > tt в каждой точке рассмат­ риваемого тела происходит разгрузка и последующее знакопере­

менное нагружение силами F i(x ,t), Ri(x,

t), которые

вызывают

напряжения

(я, t) и деформации гц (х, t). При этом

 

Ч~ F\ ^

G\jlj = R{

на S\ 2e$j =

и и -J-

(28.4)

В области

физической

линейпости (а* ^

а*)

имеют

место со­

отношения (27.17)

и (27.19) для

~ е**:

 

 

 

 

 

 

 

t*

 

 

 

 

 

 

2G0d*j = s*j +

j Г* (f* — т) s*j (т) dr,

 

 

 

 

 

.

 

°

Г ) )

= а«,

t*

 

 

(28.5)

3K0 (e*i -

За ( Г -

=

t - f , >

О,

где величины, отмеченные звездочкой, введены зависимостями

(27.16).

 

 

 

 

В области физической нелинейности

( о * > а * ) справедливы

соотношения (27.18)

 

 

 

 

t*

 

 

 

2G04 =

*?,/? (о*) + J Г* (t* -

т) fi

(0 *) 8*и dr,

 

 

О

 

 

 

ЗЛ'0 ( е * -

З а ( Г - Г ) ) = а*,

о* =

(SJV U) 1/s,

(28.6)

причем на границе зон физической линейности и физическом не-

г

'* /

*\

г* / *\

л

линейности

/ 1

J =

/ 2 \о8) =

1 .

Экспериментально определяемая величина о3, как и о\ есть напряжение, при превышении которого соответствующие кривые ползучести на плоскости э* ~ £* (в случае изменения направле­ ния нагружения) становятся физически нелинейными, или, что то же самое, соответствующие изохронные кривые о* ~ э* при 0 * > ов становятся нелинейными.

Если теперь учесть обозначения (27.16), соотношения и гра­

ничные условия (28.1) и (28.4), то придем к (27.20):

 

<*ijj +

Fi (*i) — Fi (**) = 0; o*jls = R\ (tx) R” (t*) на

S;

O *

^ I ♦

(28.7)

 

“f"

 

Пусть для определенности универсальные функции представ­

лены в виде следующих конечных сумм:

 

 

 

 

 

 

м

о =

2

 

«2 ( - У * .

(28-8)

 

 

 

 

fc=о

/

 

 

 

 

 

7712

 

 

 

 

/**(<T*) =

2 « » k( S H

(28-9)

 

 

 

 

fc=0

V W

 

где а£га£, ai k%a2hi oc/t, Ра, аЛ,

— постоянные величины.

 

Сравнивая теперь соотношения (28.5) — (28.7) и (28.9)

с соот­

ветствующими соотношениями

(28.1)— (28.3)

и

(28.8), заключаем,

что величины Оц (х, £*) и 8 i;- (я, £*) совпадают

с

соответствующими

величинами 0 *j(:r, t) и е^(я, £),

если

только

в

последних

осуще­

ствить следующие замены:

 

 

 

 

 

 

о*-* а?; а?, а*,

aft,

а?\

 

а£,

р£;

 

F; (ж,

- К

(ar. f*);

Г (g)

 

Г* (|);

 

Ri-*- («, О — Rl {x, t*).

Таким образом, если

Oil = fii («, t, a*, a?, . .. ) ,

e-j = cpij lx, t, a*, a}, . .. ) ,

то, согласно приведенному выше утверждению, при Г(^) = Г <:(?) имеем

f ij (*£>

># 1 »•••)»

— Tij (#* ^ а3* , а*к

(28.10)

После определения величин a*j и

по формулам (28.10) ис-

комые напряжения ai;* и деформации

Sij находятся по формулам

(27.16), в чем и следовало убедиться.

 

 

Совершенно аналогичным образом доказывается теорема для любого и-го нагружения. Уравнения состояния запишем в форме

(27.22):

*

*п

 

 

 

 

2С0о*’1= st? J г ; {£ -

т) s*?dx,

а* <

а,*";

О

 

 

 

 

2G03*;1

— т) /

( О

S i? (1т,

On Os ,

ЗК0(г*? -ЗаТ*п)

tn t

tn—1

 

Оц ,

 

Кэтим уравнениям присоединяются соотношения (27.23). Универсальные функции представим в виде (28.9)

где а£А, а£, — постоянные величины. Теперь следует по­ вторить рассуждения, приведенные выше для случая переменного нагружения после первой разгрузки.

При этом, например,

 

 

, (

.*

ПН

\

 

 

ОЦ

j

j

1Os , d\ ,

)

Используя (27.21)

 

 

 

 

 

о*? =

( - I f

" ^ 1 -

оЪ),

г*? =

( -

1 )п ( е Г - еЪ), (28.11)

искомые напряжения Оц

и деформации

гц найдем по формулам

оЪ =

а'и

2 (— l) ho*j,

=

elj — S ( — I)'* <4*•

 

 

h=2

 

 

 

 

k=2

§29. Приложения

1.Обратная ползучесть толстостенного цилиндра. Полый ци­ линдр кругового сечения внутреннего радиуса а, внешнего ра­ диуса Ъ нагружается внутренним давлением pit). При этом в слу­

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

» =

С(0

,

3

гр

Т0 = const,

g (О I 3 гр

+

-2

a 7V%

ег = ---- 0>

 

 

 

 

 

 

(29.1)

ео =

+

4

а7,о, » s =

2 ( ^ ) * +

- fa*tf

Здесь cU) — неизвестная пока функция времени, a — коэффи­ циент линейного температурного расширения, Т0— однородная

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

Из уравнений (27.2) при

= ф2 =

Ф находим

 

 

 

 

 

t

 

 

a\ G—

=

(ее — ег) Ф(з) — J R (t — т) (ее — ег) ф (э) dx.

(29.2)

О

 

 

 

о

 

 

Подставляя в уравнение равновесия

 

 

 

 

d<5,.

\

 

 

 

 

 

= _ ( а е - а г)

 

выражение (29.2), с учетом (29.1) при условии

 

 

 

фЫ = ф0з~в,

ф0> 0,

б > 0 — const,

(29.3)

получим после интегрирования

 

 

2^стг= с0 + 2ф0 |^c(f)x(r. с)~ J Д(* — т)с(т)х(г, с) drj,

 

 

Г

 

 

 

 

(29.4)

Х(г, с) = j

-

 

 

 

 

 

а

- 3

2

 

r o J

 

 

 

г 3

г4 + 2 а

 

 

 

 

 

 

 

Удовлетворяя

граничным

условиям

ог = —р, г = а? определим

произвольную постоянную Со, после чего

 

а (г , t) -f*

р (*)

 

 

(29.5)

 

------ = с (*)х(г. с) ~ J Д (*—*)*(*) Х(г»c)dT-

 

О^О

 

п

 

 

Из граничного условия аг(Ь, £) = 0 с учетом (29.5) находим интегральное уравнение для %i(£) = c(t)%{b, cU)):

*

® ^ ; = X i W - J f l ( * - T ) X i ( T ) d T .

Его решение запишем в виде

t

4G0<p0c(t)x(b, с (0) = P (0 + .f Г (* — т)р (т)*г,

(29.6)

6

 

где Г(|) — резольвента ядра Д(|).

После определения из (29.6) функции c(t) согласно формулам (29.1), (29.2) и (29.5) находятся перемещение, деформации и на­ пряжение. Как видим, задача свелась к вычислению интеграла (29.4) и решению конечного уравнения (29.6) относительно c(t).

Пусть теперь при t * * t t давление в цилиндре мгновенно сни­ мается, при этом в точках цилиндра будет происходить разгрузка;

требуется определить

п

п

п п

определенности

и , ег, ее, сгг, OQ. Пусть для

функция <р* (э*) = срх =ф *,

входящая в (27.19), имеет вид (29.3)

 

 

Ф* =

фо* (а*)-6*.

(29.7)

Для нахождения величин

 

и* (г, t*), е* (г, £*), ее (г, £*), с* (г, £*),

ое (г, £*), связанных

с

искомыми величинами

соотношениями

(27.16), воспользуемся теоремой § 28. Согласно этой теореме из

(29.1)

находим

 

 

 

 

 

 

 

 

 

 

 

Р

1

Б0 ----- ег —

^2

1

а* = /

2

(29.8)

Здесь учтено, что

Т% =

То — То = 0 .

Функция <:*(£*)

удов­

летворяет уравнению (29.6)

 

 

 

 

 

 

 

 

 

 

0Фос*Х* (Ь, с*) =

(*i) V

(**)>

(29.9)

 

 

 

<*

 

 

 

 

 

 

 

 

 

i|>* =

1 +

J Г* (|) d£,

р*

=

р (h) =

const,

 

где согласно (29.4) 0^>-Т£ =

0)

 

 

 

 

 

 

%*(Ъ, с*) =

 

1,2—26* _

2- 26.

 

 

п*

6 * < 1 .

л

——--------- ^------— п - =

^ ,

(29.10)

A

v > +/

2в»/г2 (1 —

6*) с®* (аЬ)2-26*

 

с**

 

 

Из (29.9) и (29.10) определим с*(£*):

 

 

 

 

 

 

 

=

/р (Ч) Г ( М

\тгз;

 

б

ф 1#

(29.11)

 

 

*

^ AG^

D *

)

 

 

 

 

 

Запишем еще формулу для а*, используя (29.5):

g? + P(fi)

1

6 .)

4СоФо*

2б*/22 (1 -

 

X

Л . « , - , ) с - * - ( г ) * }

Напряжение ае

находится из соотношения (29.2), в кото­

ром следует заменить ее — ег ->-ее — ег, <р-»-ф* и £-»-£*. Искомое величины найдутся из зависимостей (27.16), например, формула для перемещения будет

и" (г , t) = и1( г , t) у\

: iGo < D*

I

где, согласно (29.1) и (29.6),

4<?оФ0с'(*)х (Ь,с') =

 

= P(*i) +

<1

 

 

t

 

j T(i —

х)р(т)^т + р(^) J r ( i - x ) d T .

При

и " {г,

t ) стремится

к

некоторому остаточному пе­

ремещению u°(r):

 

 

 

 

 

гц° (г) = с0 —

 

 

б!1!< 1 ,

 

 

 

 

 

 

оо

4СоФосоХ (Ь. со) = РФо.

Фо =

1 + J г (£)dli

 

р =

const,

i|)J =

1

+

j Г* (I) dl,

 

 

 

 

 

 

О

Если принять

JT1о=

0, то для и°(г)

получим конечную формулу

D =

1 / 1

1 \ 1

1 _______м ________м

 

b2-26 J 1 — 6 Э D* =

2б/22 U 2“ 26

2 б*/г2 (1

 

 

(1 — 6*) U 2-26* б2“ 26*у

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

дует

положить 6 * = 0 , фо = 1 .

 

Не

вызывает трудности рассмот­

Рис. IV.9. Изгиб части кругового

рение и общего случая, когда при

цилиндра.

нагружении и разгрузке сущест­ вуют и области физически нелинейного, и области физически ли­ нейного состояний [2 ].

2. Циклический изгиб части кругового цилиндра. 10. Н. Работпов обратил внимание на задачу о чистом изгибе части кругового цилиндра [215], являющуюся в некотором смысле аналогом за­ дачи об изгибе кривого бруса. Воспользуемся этой постановкой и приведем решение задачи о переменном изгибе части цилиндра

12 В. В. Москвнтнн

'(рис. IV.9) с использованием уравнений состояния для вязкопла­ стической среды.

Отличные от нуля компоненты напряжений аг и а0 не зависят от угла 0 и удовлетворяют уравнению равновесия

± (гаг) = ае,

(29.12)

из которого следует, что по любому сечению 0 = const суммарное усилие равно нулю.

Компоненты деформаций связаны с перемещениями следую­ щими формулами:

__

dii

 

__ и

,

1

dv

V =

1

ди

,

dv

v

 

®г =

Иг'

 

 

 

 

 

90 :

г

т

+

дг

г

 

Из условия совместности деформаций и условия несжимаемо­

сти при ez = 0

 

следует

 

 

 

 

 

 

 

 

 

 

 

 

 

ее — ег =

 

г

Л (г),

ее +

£г =

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим через с(£) значение г, при котором е0 — ег =

0. При

- этом

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ee _ e r =

4 ^

-

l ) ,

А > 0,

 

=

 

 

1 )“

(29.13)

З&аюльзуя

(29.2), (29.3) и (29.13), получим

 

 

 

"г) =

2 6/2ф0

A1- 6

-

l )

-

 

 

 

 

 

 

 

— т )Л 1 _в(т)(- ^- — i )

 

d t j

 

(8 < 1 )

(29.14)

в растянутой области

(г <

с)

и

 

 

 

 

 

 

 

щ (<*е - Or) =

-

2 lV2(p0 1 Л 1 " 6 ( 1 -

j

 

 

 

 

 

 

 

 

-

 

 

— т )Л 1_б(т) ^ - - ^ - +

l j

J t j

(29.15)

JB сжатох! области > с).

Подставляя (29.14) и (29.15) в уравнение равновесия (29.12), получим в растянутой области

х

су.

= А' 6 (0 Xi(r - c)

-

— т) Лг

6 (т) Xi (Г, с) dx,

26/22СФ

 

 

 

 

(29.16)

с> “Я ? -1) Т-

 

 

 

и в сжатой ооластн

 

 

 

а

 

 

(* - т) ^

6 СО X* с) dx,

Р72

= -41 6 (О X* (г>с) -

f

^ ZCrOfpQ

 

О

 

(29.17)

 

 

 

 

При этом мы удовлетворяем условиям о, = 0 при г = а и при

г ~Ъ .

Уравнение для определения c(t) мы получим, приравнивая ог при г = с:

из которого следует, что с не зависит от t.

Для определения о из (29.18) можно воспользоваться методи­ кой, предложенной в [215].

Единственная оставшаяся неизвестной величина A U) найдет­ ся из условия

ъь

М (t) = J <r0 (г, t)r d r = § вт(г, о Г dr,

а

а

которое мы преобразовали с учетом

(29.12) и условии аг = 0 при

г = а и при /• = 6 . Используя при этом (29.16) и (29.17), получим

 

= л '~‘ (,) - . ( « < « - " ) л '~‘ (т)

откуда находим

^

^1_б(0 *

3/5“ ----( Af W + f Г(* - T )A f (Т)(*Л (29.19)

2б'-2СоФох\

^

/

где обозначено

с

Ь

X j x 1(r,c)rdr + j

%t(rxc)rdr.

ас

После определения величин с и A it) становятся известными все искомые величины.

Пусть теперь при t = tx начинается разгрузка и последующий

переменный

изгиб моментом

M "(t).

Согласно

(29.13)

можем

записать

 

 

 

 

 

ее — е*

— 1],

е 2 + е *

= 0, t* =

t — tv

(29.20)

Постоянная

с* при условии

(29.7)

удовлетворяет уравнению

<29.18)

 

 

 

 

 

 

аш - r t - kс *

- s r t

Функция

А * (£*)

определяется формулой (29.19)

Л 1 " в* { h )

= 2a* /32G0tp*X* ( М * ( * * } +

! Г * { t * ~ Т ) М * (Т ) Л ) «

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

 

С* , ^

v . д

ь

 

- j ( 4 - ‘ Г ’ ’ * ■+К 1 - - ^ Г * ' * . м * « . ) -

а

 

Сф

 

= м (tj - М" (t«).

Аналогично определяются и напряжения, отмеченные звездоч­

кой. Например, используя (29.16), определим сгг (гг £*) в области г < с * :

*

Ч

('*) XI (г) -

Xi (г) J П* (h - Т) А\~ъ* (т) dxt

" "

"

"

*

#

Искомые величины аг, а©,

ег и е©

после определенияогг,

а©»,

er , 8Q находятся по формулам (27.16).

3. Переменные нагружения трехслойных пластин*). Рассмот­ рим сначала задачу о нагружении из естественного состояния трехслойной пластины. Для несущих слоев (с толщинами ht и hz)9 которые считаются вязкопластическими, принимаются гипотезы Кирхгофа — Лява. Для вязкоупругого заполпителя (h3= 2с) при­ нимается гипотеза о прямолинейности нормали и недеформируемости в направлении оси х3 [52] (оси xt и хг расположены в срединной плоскости заполнителя).

Используем следующие уравнения состояния для слоев плас­ тины U\ j = 1, 2):

4 ; = 2Gk0 (1 - со*) _ J д * {t _ т) ф* {эи) э\. (т)

окц = 3/4 (4; - 3ahT0),

 

(29.21)

т4= G[3) ^ )i (1 — co(3)) — j

i?(3) (t — т) ф<3)'ф1(т) dx

 

Здесь к = 1, 2, 3 — номер слоя,

—разность между

углами по-

ворота нормали в заполнителе и несущих слоях, Х{ =

(з)

si3 — каса­

тельные напряжения в заполнителе.

Согласно принятым выше гипотезам, перемещения и дефор­ мации в слоях можно выразить через перемещения точек средин­

ной плоскости заполнителя i\, w и поперечные сдвиги

[521:

г?з =

+ u\i =

(«!, хг, t).

(29.22)

Отсюда после интегрирования по х3 получим

w(i3) = Щ+ Z3l|>i — x3wti.

В плоскостях соприкосновения слоев

u (i1} = Vi + c\\>i cw%i, u(i2) = Vi c\pi + cwti.

Учитывая гипотезы Кирхгофа — Лява, находим перемещения точек внешинх слоев, расположенных на расстоянии х3 от опор­ ной поверхности:

l4 X) — Vi + С\l>i — x~wti, и ? = Vi — C\|)i — x3wti.

Для деформаций в слоях имеем следующие выражения:

=

Ч~

Ч~

4“

(29.23)

 

8ij

= &ij Ч" C^ij Ч~ ^ij)j

^i3 =

 

*) Этот пункт параграфа иаппсан Э. И. Старовойтовым.

Здесь

= X Ohi +

y-ij ■-= — w.ij, eij = -^ -(i\j + Vj.i), (29.24)

e,j — малые деформации срединной поверхности заполнителя. Для шаровой части тензоров деформаций находим

е(1)

 

1

 

 

 

е(1) — -Y

= ~Т

+

*з*гО,

 

( 2)

 

_1_ (е .. _ с^н +

 

 

е(2) =

=

;r,Xj,),

(29.23а)

е'.я)

 

,

 

Х3У.ц).

 

е(3)= -у- =

— (eii 4" * аЦ>« +

 

Введем усилия и моменты

 

 

 

Ми = f o'ihdx,,

 

Tij =

f а?;Аг3,

QI - f o-fcZxg.

(29.25)

Л/i

 

 

,lft

^3

 

Здесь интегралы берутся по толщине соответствующего слоя /гА, причем h = x3- с, ga = х3+ с, £3= х3.

Уравнения равновесия и граничные условия можно получить из вариационного принципа [52]

Т ц j + Pi = 0f

Я у, j — Q i = 0, Л/у. у + g =

0,

( 29.26)

где введены обобщенные усилия и моменты

 

 

 

Tij = 2 Т\и

Hi} = 3/W +

с ( 71У -

rtf),

м и =

 

 

к

 

 

 

 

 

 

 

 

 

 

 

= 2

М ки +

С ( 71У -

Т\У ). (29.27)

 

 

 

к

 

 

 

 

Граничные условия, например, вдоль линии

 

имеют вид

Та =

Я и = Я „,

Л/и =

Л/и, Л/n.i + 2Л/12>2=

^ 1# (29.28)

Здесь р,-, g — поверхностная

нагрузка,

Tih Я п, Л/ц,

Qi — контур­

ная нагрузка.

 

 

 

 

выделим в Тц, Д/у, IIih Qt

Учитывая

(29.27), (29.25)

и (29.21),

линейные слагаемые

Тij = Тij Ту, Л/у = Л/у+ Л/у,

Н ц = Н ц + Н ъ <?с = <?«+ <?7. (29.29)

Здесь

= 2 f [ г с ^ + б ^ Ц е ^ - З а ^ о ) ] ^ , h L

- Z*j =

2 j

 

^2Gji |э?У + j‘ Я* (< -

T ) q > 4 * ) ] dx3l

М ц =

2

f

[2Gj,''i +

ЬцКк0 (e\\ - 3aftr 0)] \hdx3,

 

h4

 

 

 

 

 

=

2 2

 

j

[ G '‘

+ j' Д* ( * -

т) q>4dT hdx„ (29.30)

=

Arg> +

c ( r t f -

H f ), На =

Щ Р +

с (T $ - Y $ ) ,

0

 

c-

tyidx3.

 

 

 

Qi =

Go3)

j

 

 

 

 

 

 

—C

 

 

 

Qi = GQ3) j

 

£%со(з) +

j R(3)(t — x) t|)j<p(3)dTj

dx3.

Преобразуем теперь уравнения равновесия (29.26), учитывая (29.29) и (29.30),

T ij,j ~\г P i

~Г Т- i j j =

0, H ij,j

Qi “Ь

^

О,

 

 

 

 

 

(29.31)

"Ь ? +

= 0.

 

 

 

Аналогично преобразуются граничные условия (29.28)

Тц = Тц Т\1,

=

Н Мч

=

Л/ц

Л/11?

 

 

 

 

 

(29.32)

А/ИД “Ь 2М12,2 =

— -Л/цд — 2л/12,2 «

 

 

 

Используя соотношения

(29.22)— (29.26),

можно выразить все

усилия и моменты, входящие в уравнения равновесия (29.31) и граничные условия (29.32), через перемещения vu w и сдвиги ф,-. В результате получим замкнутую систему уравнений и граничные условия относительно пяти неизвестных vu v2, го, \|)i, \|)2. Эту си­ стему ввиду ее громоздкости мы выписывать пока не будем, а рас­ смотрим один метод линейных приближении для ее реализации.

Предположим, что интегральные слагаемые в (29.21) имеют малый параметр, например, <р*<1. Поскольку все о )*< 1 , возмо­ жен при этом следующий метод итерации: для любого ттг-го при­ ближения имеем задачу

Т?},} + p i + T i p x = О, Я ш — @ t + H i j j 1 = о ,

(29.33)

M ? j,ij + Q + M i } * = = О

с граничными условиями (29.32)

 

Т?Х= ТЛ - Т Т \

=

М ? х = М п - Мп~\

м?1л + 2Д/Гг,2 = Qt -

MuTi - 2МТы

(29.34)

причем величины Г у 1

и т. д. при тп — 1 принимаются равными

нулю.

 

 

Усилия и моменты выражаются через искомые перемещения

по формулам

 

 

 

 

 

 

 

 

 

 

Tii =

(vb + v?,i) 2h Tfc+ с (Vi — 7г) «

i

+

C i ) -

 

 

 

 

 

— 2 (®i — ai)

 

+ 6у

2 b ft + c(bi — b2)

 

 

 

 

- ( с 1- с 2)ш™11- З Т 0'2 а Ь к Щ ,

(29.35)

=

(% — a2) (v™j +

vfti)

+ c ^

+

o2 +

-y-j (фу,- +

 

 

— 2

+ d2+ -J J ]

+

Sy [^(ci — c2) yu + c ( ci + c2 +

"^•)

 

— (?1 + £ 2 +

-§-) WJI — ЗТ0[а (1)К(о \

+

 

 

 

 

 

 

-

Л

<

% ( с

+

^ ) ) ] ,

(29.36)

H™ =

(Vi — Y 2) (v™j +

v™i) + c2

+

Y2 +

“j f ) ('Фи +

Ф]?*) ~

 

2c ^ax + a2+ - y j w?,%+

6y [(^i — b2) v™i+ c2 ^bx +

b2+

-y-j 'Фы

c (cx + c2+ -y-j Wji — ЗсГ0 (a^K^h,! — a(2)^o2)^2)]»

(29-37)

 

 

 

^ Г =

О з Ф Г .

 

 

 

 

 

Здесь введены следующие обозначения:

 

 

 

 

 

 

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

-

=

2 2

Gh0 J Гэ\ТC0hm +

f Rk(t -

т) cpfc (экт) э\Г (т) d r] dx,,

 

 

h

ч L

 

о

 

J

- Щ

=

2 S

Gh0 J Г

+

f Д* (* -

T ) cp'1(з*ж)

(т) dx] l hdx„

 

 

h

i'h L

 

о

 

J

 

 

 

Щ = M li)m + с (ГЙ )т - r (i])m) ,

(29.38)

- Q ?

=

G[s) J jVV

3)m + j

R(3) (t -

г) Ф(3) (э(3)т) Ч»Г (T ) drjdx?.

Таким образом, в каждом приближении имеем задачу линей­ ной термоупругости с известными внешними нагрузками.

Пусть теперь, начиная с t = tu осуществляется разгрузка и последующее нагружение усилиями обратного знака. Следуя (27.16), введем для всех слоев разности

 

(^*) ^ °ij (0 Gij (О» ®ij (^*) = ®ij (0

(£)» (29.39)

где

^ ^ 0. Аналогично для перемещений и углов сдвига:

U *

и т. д.

 

Ограничимся случаем, когда материал заполнителя обладает свойствами линейной вязкоупругости. При этом для величин, от­ меченных звездочкой, запишем уравнения состояния типа (29.21)

s lf

=

2G h0 [ *

(1 -

CD*ft) -

j R

l ( f , - T) ф** (э*й)

 

=

3K0 (s*l -

3ah (T 0 -

To)),

(29.40)

Ti

-

G <?> [ <

(1 -

B-«>)-J

« » -

г» Ф- « ч ; (т) jл .

#

 

 

 

 

 

 

 

 

Деформации & i j ,

перемещения

и г|)| связаны между собой

соотношениями (29.22),

(29.23),

(29.23а)

и (29.24), которые

мы здес*> не повторяем.

Уравнения равновесия для величин со звездочкамп записыва­

ются аналогично (29.26)

 

 

+ P i =

0,

Н i j j Qi = 0 ,

М i j tij + g* = 0,

а гранитные условия — аналогично (29.28)

7ц =

?и,

# £ = / / ,*,

Л/и = .1/Гг?

 

VI/11,1 + 2Л/Г2,о

$1.

Усилия и моменты, отмеченные звездочками, связаны с пере­ мещениями Vi, w* и сдвигами соотношениями, которые сов­ падают по написанию с соответствующими соотношениями между Ту, Hij, Мц и Vi, w, ^ при нагружении из естественного состояния, если в уравнениях (29.21) принять со(3) = 0, ф(3) = 1. Предположим, наконец, что универсальные функции в (29.21) и (29.40) различа­ ются только константами щ, bt и а*, Ъх :

со*1=

со71(aft, а{),

cpft =

cp/l(a7i, bt),

co*ft =

co*ft (a*ft, я*),

cp*ft =

ф*л (a*A, b*).

Как видим, при указанных выше предположениях имеет мес­ то теорема о переменных нагружениях, приведенная в § 28. Поль­ зуясь этой теоремой, можно найти величины ь\ , и>*,%,если опре­ делены vh w, ^ — перемещения и углы сдвига при нагружении из естественного состояния. Для этого достаточно в последних за­

менить

давление q{t{) на q* = q ( t j q" (f*) и константы сц, bt

7 ^

на я;,

bt.

Г Л А В А У ДИНАМИЧЕСКИЕ НАГРУЖ ЕН И Я (КОЛЕБАНИЯ)

УПРУГОПЛАСТИЧЕСКИХ И ВЯЗКОУПРУГИХ СИСТЕМ

§30. Колебания упругопластических систем

сконечным числом степеней свободы

Рассмотрение циклических нагружении с учетом сил инерции мы начнем с постановки задачи о колебаниях упругопластических систем. Будем рассматривать только случаи одноосных напряжен­ ных и деформированных состояний, поэтому для линейно упроч­ няющегося материала упругопластические свойства можно описать с помощью диаграммы зависимости напряжений о и де­ формации е, представленной на рис. V .I. Согласно этой диаграм­ ме, зависимость между напряжением и соответствующей деформа­ цией в любом случае нагружения пли разгрузки можно выразить одним уравнением

где

 

 

 

 

 

 

о = о °

+

г)Ее,

(30.1)

 

 

 

Е ф — 1) е£, i|> =

1 — /ф, е£ = к (е0 + ае8).

 

а° =

оф/ссг* +

(30.2)

Здесь Е — модуль упругости,

е0

пластическое слагаемое

е0

деформации перед началом раз­

 

 

 

грузки,

к — коэффициент

уп­

 

 

 

рочнения материала (к =

0 при

 

 

 

упругом

нагружении

0 <

к < 1

 

 

 

при наличии упрочнения, к = 1

 

 

 

для идеально

пластических

ма­

 

 

 

териалов),

05

и Es — величины,

 

 

 

характеризующие

предел

теку­

 

 

 

чести

материала

при

растяже­

 

 

 

нии

из

естественного

состоя­

 

 

 

ния. В приведенные соотноше­

 

 

 

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

 

 

 

Р, которые

дают

возможность

Рпс. V.I. Значения параметров ос и

выделить

 

в

уравнении

(30.1)

 

Р при циклическом нагружении.

определенную

прямую;

зпаче-

 

 

 

нпя а и р

для различных участков нагружения и разгрузки при­

ведены на рис. V.I.

 

 

 

 

 

 

Представление зависимости о - е в единой форме (30.1) имеет то преимущество, что позволяет избежать написания в общем слу­ чае четырех уравнении, характеризующих последовательно про-

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

Перед тем как перейти к рассмотрению конкретных систем, за­ пишем еще зависимости для переменного предела текучести (рис. V.1):

е* = еo + 2aes>

Os = a0 + 2aas.

(30.3)

Продольные колебания упругопластического стержня [158]. Рассмотрим стержень постоянного поперечного сечения F с мас­ сой единицы длины т и длиной Z, жестко закрепленный, напри­ мер, левым концом. Разделим I на п равных частей и поместим массы mil71 в середину каждой часты. В результате получим си­ стему, положение которой можно описать п параметрами. Примем за эти параметры uhit) — продольные перемещения точек, содер­ жащих массу, относительно положения, в котором все элементы стержня не имеют ни напряжений, ни деформаций.

Будем считать, что на к-ю точку системы действует внешняя сила Pk(t) и, кроме того, на правый конец стержня действует си­ ла pit). Выделив к-й элемент двумя близкими сечениями, получим

mh-rjr = Ph + F (pk+x — <rft).

(30.4)

at

где ак— напряжение на участке между точками А: — 1 п к. Де­ формация на этом участке равна

eft = -у- К — Wfc-i).

(30-5)

Для составления разности а®+1 — сг£ воспользуемся представ­ лением (30.1), при этом, согласно (30.2),

a®+ i —

<j{[ = has(cXfc+iPft+i — « аРа) + Е [(P A+ I — 1) &oh+i — (Р* — * ) e®ft] •

 

 

(30.6)

С учетом (30.1), (30.5) и (30.6) соотношение (30.4) преобра­

зуется к виду

 

d?u

а2п* [и* (ф* + Фа+I) Uh-i'h ~ Uh+ityh+i] =

- X +

dt

 

 

 

— 4h(t) -j- ка?n2

а— «ft+i + 2a/t+1p/l+1 — 2CC*PA) +

+ (Pft+l — 1) ( Mft+1 — uft) — (Pft — 1) i uh — W ft-»)]' $0.7)

Здесь иь— перемещения, связанные с предшествующей макси­ мальной деформацией боя соотношениями (30.5). Через аг и qk

обозначены величины а2= EF/(ml2), qh= npk/(ml).

Практически, по-видимому, удобно будет использовать диф­ ференциальные уравнения, записанные в несколько иной формег если сохранить величины е?ь в соотношении (30.6):

+

а2« г [« а + $h+i) — “ fc-i^ft — Щ+Лк+i] = qh(*) +

где

(30.8)

 

7ft =

— [A'BS(aft+ifWi — ocfcPft) + (Pft+i — 1) £ok~i — (Pft — 1) ejft].

 

(30.9)

Заметим, что возможна и такая постановка задачи, когда из­ вестен закон движения u(t) правого конца стержня; в этом случае

ип+1 = ип.

Системы (30.7) н (30.8) являются по виду линейными диффе­ ренциальными уравнениями, однако наличие параметров ак и принимающих различные значения в процессе упругопластиче­ ского нагружения и разгрузки (различные участки диаграммы (рис. V.1), не позволяет написать решения uhU), имеющие одно и то же аналитическое выражение за все время движения. Поэто­ му придется интегрировать уравнения по интервалам времени, используя условия сопряжения перемещений и скоростей на гра­ нице интервалов.

Пусть на некотором интервале по времени параметры ак и для всех значений к остаются неизменными. Представил! частное

решение однородны* уравнений (30.8) в виде

 

u'h{t) = CfcSin(|i£ +

Ф).

(30.10)

При этом

 

 

— И-2** + a2n2 [с#(t|>ft + >ft+1) —

— cft+1i|)ft+1] = 0.

(30.11)

Это — система ддгебраических однородных линейных урав­ нений относительно с*. Поскольку нас интересует случай, когда хотя бы одно значение ск было отлично от нуля, определитель системы (30.11) должен равняться нулю:

2 ^ + % - «'

= 0,

(30.12)

% —

где о 2== fiVaru2.

Отсюда следует, что рассматриваемая механическая система относится к типу хорошо изученных штурмовых систем [42].

Для функции Дл(со2) справедливо следующее рекуррентное соот­ ношение

Да (ю2) = (Фа + Фа+ i — (°2) Да — i (°^) — Фа^ а- з (<D )•

Пусть 0 8— пронумерованные в порядке возрастания корни уравнения (30.12) (положительные вещественные числа). Под­

ставив со? в уравнения (30.11), определим из последних п зна­ чений csh/c[ = Ak,s после чего общее решение однородных уравне­ ний с учетом <(30.10) запишется в виде

и'к (t) = 2 AhCssin

+ <ps)

(cs = c0-

(30.13)

S=1

 

 

 

Частное решение системы неоднородных уравнений (30.8) бу­ дем искать в виде

u'h(t) = B h+ D hsin vt,

(30.14)

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

quit) = qlsinvt.

При этом

Вк(грй + th+i) — Bk-ityk — Bh+ityk = Yft,

(30.15)

Dkv2+ fl2/i2 [О й('I’/t+ ^/t+i) — Dk-ityk— ^ft+i^fc+il = g®,

к = 1,2,

n. (30.16)

Определитель системы (30.15) отличен от нуля, поскольку он получается из определителя (30.12) приравнением со2 нулю, а п

значений cof являются отличными от нуля корнями уравнения (30.12). Определитель системы (30.16) совпадает с Дп в случае, если частота изменения внешних сил v совпадает с одним из зна­ чений р,в, и, следовательно, всегда отличен от нуля, за исключе­ нием последнего случая. В этом последнем случае решение может быть построено по методу Б. В. Булгакова [25].

Учитывая (30.13) и (30.14), запишем общее решение системы (30.8)

п

 

ик ( t ) = 2 Aicssin (yL,t + Ф.) + B h + Dhsin vt.

(30.17)

3=1

 

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

Запишем теперь дифференциальные уравнения упругопласти­ ческих колебаний п масс тк, расположенных на вертикальном не­ весомом стержне [1581:

£ 4 + all

%

 

 

 

 

dt2

xh

 

 

 

 

где xh— расстояние между массами к — l

и к; g — ускорение си­

лы

тяжести; и0 = 0;

qk =

al = EF/(tnhiy, %k = \|)fe+1eJ(+1 —

+ (Pft+i — 1)ejft+x — (Рь — 1)Boh+ kes (ah+ifib+i

=

 

 

 

 

П

 

 

=

8* t+

(Uk — uk- 1 ),

e*t =

2 ™h <

es.

 

 

 

k

 

h

 

 

 

Крутильные колебания при

наличии

пластических

деформа­

ций [160]. Рассмотрим колебания п дисков, насаженных на не­ весомый вал кругового поперечного сечения, под действием произвольных внешних сил. Эти силы приводятся к моментам Mh(t), приложенным к соответствующим дискам. Пусть

..., /п — моменты инерции дисков относительно оси вращения. Уравнение движения к-то диска будет

d \ b

ь

= 1, 2, ... , п). (30.18)

h - ~ г - = м ш -

М к + .\Ih

Зависимость между крутящим моментом Mh и соответствую­ щей круткой 0h при наличии пластических деформаций аппрок­ симируем линейными соотношениями (при упругих деформациях и разгрузке это — точные линейные зависимости):

М h =

Ml. -f tykBffik,

 

(30.19)

где

 

 

 

 

 

M l - афkkhM sk+

(pfc -

1) ВЖь,

(30.20)

Msk = Bh®sh, 9o/t = kh (02 +

afc0Sfc),

ф/i =

1 — АгйРл,

(30.21)

Bh— жесткость на кручение

участка стержня

между {к — 1)-м и

&-м дисками, к — коэффициент, характеризующий наклон прямой ^*(0J при пластическом деформировании. Значения параметров а* и Рд в различных случаях нагружения и разгрузки те же, что и на рис. V.I.

Линеаризация (30.19) позволяет, как и в случае продольных колебаний, с помощью введенных параметров ак и написать единые системы дифференциальных уравнений, характеризующих движения как при нагружении, так и в стадии разгрузки. При этом решение строится по участкам, в каждом из которых пара­ метры <хк и р* не изменяют своих значений, и используются

условия сопряжения соответствующих углов поворота и скоро­ стей на границе участков.

Пластические деформации после разгрузки появятся при условии, если

 

Мк = M oh + 2ahM 3h,

ek = Q°h + 2ahQthi

M ok = M k(Q°h).

Поскольку

 

 

 

 

 

 

 

 

 

 

efe --- Ф* T -t* - 1

 

(30.22)

 

 

 

 

 

lh

 

 

 

(lh— расстояние между дисками

к — l

и к),

используя (30.19) и

(30.20), преобразуем уравнения

(30.18)

 

 

^2_

 

 

 

 

 

 

 

 

— г

+ аI [—

+

Hft^ft+i) Фа — Mft'l’ft+ifl’ft+il =

 

 

=

Vft +

^

 

(* =

1.2........n).

(30.23)

где

 

 

 

 

 

 

 

 

Ук ~

ак [otft+lPfc-|-lA*A+lH'A^+l®afc+l

 

°&ftPк^к^к^зк

 

 

 

 

 

— (Ра— 1) ®okh + (Pft+i— 1) H'/A+I^OA+I] 9

 

Ук

Вк+1

 

ak

 

 

(30.24)

 

 

B h

 

I klk

 

 

 

 

lh+i

 

 

Заметим, что если концы вала свободны, то <р0= ф1, фп+1 = фп;

если же концы жестко закреплены, то ф0=

фп+ 1 = 0. Возможны

и другие условия для ф0и фп+1.

 

 

 

 

Предположим теперь, что в течение некоторого отрезка вре­

мени параметры ак и [J* не изменяют своих

значений. Если при

этом частные решения

однородных

уравнений (30.23)

искать

в виде

 

 

 

 

 

 

 

то

Фh(t) = /а sin Ш

+ х),

 

(30.25)

 

 

 

 

 

 

 

 

— ©2/к + а1[— 'Фа/а- т +

а +

И-аФа+I ) /ft — M >ft+i/ft+il = о.

(30.26)

Поскольку fh не тождественпо равны нулю, соответствующий определитель будет равен нулю:

СО

А » (со2) =

= 0.

(30.27)

 

 

■ T

+ ’I’n

 

При этом предполагается, что концы вала свободны и ц* = 1.

Для развертывания определителя Ап(со2) можно воспользо­ ваться рекуррентными формулами

А* (ю2) = |фй +

Фл+х — j Aft-i (ю2) —ФлДк-г (®2);

 

 

 

А0 = 1,

Фх = Фп+i = 0.

Формы главных колебаний

f\ определяются по формулам

Гн =

е, Aft-x 2)

(к, s = 1 , 2 , . . . ,

п),

 

Фг^з-'-Фк

 

 

где сов— корни уравнения (30.27), сл—- произвольные постоянные. Частные решения неоднородных уравнений (30.23) будем ис­

кать в виде

 

 

фА(£) =

Ak+ Dhsin v£,

(30.28)

предположив при этом, что

 

 

M h ( t ) = M o S i n v t .

 

Постоянные Ah и Dk удовлетворяют уравнениям

 

а\I Фл-^л-х + (фл +

фл+х) ^ft —Фл+х-^л+х! — Vft>

(30.29)

M h

v2D k + а\ [— tykDk-i + (tyh+ tyk+т) Dk

(30.30)

Определитель системы (30.29), составленный из коэффициен­ тов при A h,, отличен от нуля, поскольку он равен Ап(0). Опреде­ литель же системы (30.30) равен нулю лишь при v = ©,. В этом последнем случае для построения общего решения удобно ввести нормальные координаты zh соотношениями

 

п

 

 

 

фh 2

/ftz«

(30.31)

причем, так как

 

 

 

 

 

s^ P i

 

 

 

s = р,

 

то

п

 

 

 

 

 

 

z s 2

Ikfk4>k'

 

 

ft=X

 

 

Уравнения (30.23) в переменных z, имеют простой вид:

 

+

+

(« « 1, 2, .. .л п),

(30.32)

a t

 

 

 

1 3 В, В. Москвитин

ГДе

 

 

 

А* =

i t

Ih4kfl

D* =

2

Miti.

 

 

 

 

 

 

 

h=l

 

h=1

 

 

Если

теперь

частота

v совпадает

c s-м корнем уравнения

(30.27), то

г

 

 

 

/

As

Dst

 

 

 

 

Z$=

 

(Dst +

 

 

 

 

Cssin

CSCOS CDs* + — — -— COS 0)s*.

 

 

 

 

 

 

 

 

 

 

(D^

 

 

 

Для тех значений со3, которые отличны от v,

 

 

 

 

 

/ .

 

 

п

A S

Ds

 

 

 

 

zs — CgSmojgf -f- Cs costosf-J— -

+

—-----5 sinvi,

 

 

 

/

П t

 

/

 

 

 

a ;

 

<BJ- V

 

 

где

 

 

 

 

 

 

 

опреде­

c5,

c5, Cs, Cs — произвольные постоянные. После

ления z3 искомые величины cpAU)

находятся по формулам (30.31).

В случае, если v отлично от всех сов, общее решение уравне­

ния

(30.23) согласно

(30.25) и (30.28) будет

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

фл(*)

=

2

fk sin (cos* + %s) +

Ak+ Dk sin v*.

 

 

 

 

 

 

8=1

 

 

 

 

 

 

 

Если по истечении некоторого отрезка времени один

или

несколько

параметров

ak и

изменят

свои значения,

то

для

новых значений параметров по написанным выше формулам строится решение. Новые произвольные постоянные находятся,

как указывалось, из условия

сопряжения срА

и dcpft/d*, причем

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

дисками

(к — 1 )-м

и /с-м должно выполняться

условие dcpA/d* =

dcpfe_i/d*,

которое

следует из (30.22) при dQJdt = 0. Таким методом решение может

быть

построено для любого значения аргумента *, начиная

с * =

0. Разумеется, при значительном п следует воспользоваться

ЭВМ с применением приведенного алгоритма.

§ 31. Динамические нагружения линейных вязкоупругих систем. Метод усреднения. Метод замораживания

Воспользуемся соотношениями теории линейной вязкоупруго­

сти, записанными в форме (23.11), для

случая

равного

нулю

ядра объемной релаксации Н{:

 

 

 

t

 

 

 

щ = Э{’ ~ . f R ~ т) ^ dx>

=

3^ °еи

(3 1

Oij = 2GQejj + G0j — 2G0J R (t т) эц (x) dx. (31.2)

0

Уравнения движения

ри.<, tt = F i + a ijtj

с учетом (31.2) и соотношений Коши запишутся в следующем виде:

Ри%м — Fi + GQuitjj + ( к 0+

Uj,ij

t

— G0^ R {t — r) [uUj + ju j'ij') dr; (31.3)

0

O j JRJ на

Wj W’Oi на j

Ui(x, 0 )= / tU ), Uit J i=0 = cp,U).

(31.4)

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

Как известно, применение к задаче (31.3), (31.4) вариацион­ ных принципов, процедуры Бубнова — Галеркина приводят ее к интегрированию системы интегро-дифференциальных уравнений относительно конечного числа п некоторых функций zk(t)

2

t

 

 

d

+ aUi = /ft ( 0 + J R (t T )

biZi (T ) dr.

(31.5)

 

0

 

 

Здесь а\г b\

известные константы, их

2 я2, fk— известные

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

В случае одного уравнения из (31.5) получим

t

d? + ®2Z = ^ j* R (t — r)z (т) dr + f (t),

(31.6)

—oo

 

причем здесь нижний предел взят равным — <». Пусть

 

/(*) = A sin pt + В cos pt.

(31.7)

Будем при этом искать решение уравнения (31.6) в виде

z(t) = С sin pt + D cos pt.

(31.8)

Внесем это выражение в уравнение (31.6) и приравняем коэф­ фициенты при sin pt и cos pt соответственно:

А + KCRc+ W R C= Ciсо2 - р2), В - XCRS+ XDRC= Жсо2 - р2),

13*

откуда

 

 

 

 

^ _ в ш , - а ( м с+ р2~<»2)

п

,ш ?4 - м ( * д с+ р 2- © 2)

Я2Д,2 + ( Х

Я с + Р 2 - а > 2) 2 »

 

Х2Д 2 + (Х Л С + Р 2 _

о 2) 2 ’

где использованы обозначения (23.32)

 

(31.9)

 

 

 

00

 

ОС

 

л , = j R{l)sinpld%,

Л с =

ji? (E )c o s p id i

(31.10)

 

о |

 

6

 

Как видим,

точное решение

(31.8), (31.9) уравнения (31.6)

при условии (31.7) представляет собой установившееся вынуж­ денное колебание линейной вязкоупругой системы. Если выраже­

ние (31.8) представить в виде z(t) = a sin (pt +

<р), то

 

А2 +

В2'.

 

В (XRC+

р2 —гчд)2) +

AXRS

k2R] + ( Щ

+ р2 - <о2) 2 ’

(ХЯС+

р2 ~ ш2) _

В Ы г *

Отсюда в случае резонанса =

<в) при В — 0 имеем

 

а =

А

.

В,

 

 

 

tg T =

s .

 

 

Ц я 2 + Д2) 1 / а 1

Рассмотрим теперь свободные колебания, используя уравнение

(31.6) при /(f) = О,

X J R (t — т) z (т) dx

 

-f o)2z =

(31.11)

 

—оо

 

при условии z(0 ) = z0, dz/dt =

v0 при t = 0 .

 

Точное решение уравнения (31.11) записывается в виде

z(t) = A exp (pt), А, р — const.

(31.12)

Внесем это выражение в уравнение (31.11) и введем новую переменную интегрирования % = t т; после сокращения на A exp (pt) получим

 

p2 + o 2 = x J i ? ( i ) d i .

 

(31.13)

Это уравнение служит для

определения ри р2,

Теперь

общее

решение согласно (32.12) будет

 

 

z (0 =

2 Ah ехР (Put) % Ah — const.

 

 

 

h

 

 

 

Если уравнение

(31.13)

имеет кратный корень, то,

кроме

(3 1 .1 2 ), решением будет функция

 

 

 

z(t)

= At exр (pt).

 

(31.14)