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

книги / Основы математической теории термовязкоупругости

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

Величины Р (а0, I) нам известны из опыта 1 , величины Р2 (а0, I) определя­

ются

из опыта 1

и соотношений (32.17). Отсюда мы

определяем ядро

Рг «

О

 

 

 

 

Р г (0°, 0 -

Р г (0°, 0) = 27 к ,. (*) -

_

 

 

 

1

Сц

 

 

-

-§--[/>(6®. 0 + р 2(о0, 0 -

Р (о°, 0) -

Р2 (0°, 0)]}. (40.17)

Таким образом, мы получили все материальные функции, используемые в теории как функции параметра о0в виде набора экспериментальных кри­ вых. Так как все эти функции монотонно убывают по I вместе со своими производными, то их можно хорошо аппроксимировать набором экспонент

(40.18)

2=1

хотя бы методом наименьших квадратов. Число п в (40.18) принимают обычно равным от 1 до 3.

Сделаем теперь замечание относительно метода экспериментального определения ядер К , К 2, К г, К±, входящих в уравнения (32.23), (32.24) главной кубичной теории ползучести [3].

Легко убедиться, что в одномерных опытах, когда только одна компо­ нента девиатора является произвольной (например, в случае растяже­ ния — сжатия или чистого сдвига), все пять ядер не могут быть определе­ ны, а могут быть найдены лишь четыре независимых:

* (*). («) = К г (*) + (О, (40 19)

(*) —2^а(0 + К2{1), Ка(1).

При этом опыты на ползуч,есть позволяют найти только ядро К (I) и сумму ядер: К 1Ъ(I) + К12 (г) + Къ {I). Одномерные опыты при постоян­ ной скорости нагружения достаточны для определения всей группы ком­ бинаций ядер (40.19). Следовательно, если в экспериментальном обоснова­ нии теории ограничиться только результатами одномерных опытов, то ядра К г, К2 или К 3 можно отбросить. Например, при К3 = 0 из (32.23)

имеем

I

ем(0 =

+ А (*; *)]

(0 + § 1%(* т) +

2%г У — X) 8 (г, т) +

 

 

О

 

+ К* (* *)»(т. т)]

(т) 6.x,

(40.20)

причем А имеет прежнее выражение (32.24).

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

Рассмотрим связь между тензорами напряжений и деформаций самого об­ щего вида. Пусть, например, имеют место операторные соотношения (33.7) и (33.21)

$ = ^ (Е ),

Е = <?($),

(41.1)

вой производной по Фреше следует, что оператор Р может быть представ­ лен в виде

Р ( Е ) = Р ' ( Е )

+ В г (Е),

(41.2)

где || Вг (Е) || /|| Е || -*■ 0

при || Е [| -> О, а Р0 (Е) — оператор, линейный по е

и поэтому представимый в интегральном виде. Если

он инвариантен

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

Будем рассматривать квазистатические

задачи. Уравнения равновесия

имеют вид (15.2)

 

 

+

- 0,

(41.3)

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

ниями (15.3) и видом оператора Р , получим выражение тензора напряже­

ния как некоторого оператора от вектора перемещений

и

<*« = *«{“}•

(41-4)

В соотношениях (41.1) и (41.4) можно учесть влияние температуры. А имен­

но, если обозначим через е*у тензор, состоящий из суммы тензоров соб­ ственно деформаций е^ и деформаций, обусловленных температурным воз­

действием а^Т, где

— тензор теплового расширения,

 

&ъ]

^г

=

2 (^>3 ^7,0»

(41,5)

где

 

ЫцТ1

 

1,7*=

 

 

то соотношения (41.1) и (41.4) изменяются заменой и и е на

и и е соответ­

ственно:

 

 

 

 

8 = Р ( Е),

Е

= 0 ( 8 ) ,

(41.1)'

Ог1 =

Ри {и}.

 

(41.4у

Пусть, кроме того, заданы граничные условия (15.7), вообще говоря, кон­ тактного типа

«ЙV , г № \ Ь (х)/1 = Мх).

(41.6)

В силу разложения (41.2) оператора Р соотношения (41.4) и (41.4)' мож­ но записать в следующем виде:

и } + в „ { и ) .

 

(41.7)

Тогда уравнения равновесия (41.3) примут вид

 

-5Г *«{»>=

3

(41.8)

3

 

а граничные условия (41.6)

 

 

и}1, + № ик- ^

= Мх>-а1?Вк}{и}1,.

(41.9)

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

(41.8) и (41.9) Вг= 0, т. е. В^{и} = 0. Тогда получается линейная задача,

заключающаяся в определении

трех функций

перемещений гг(°) по трем

уравнениям

равновесия

 

1

^ П

- { « (0)} = - Р ^

 

(41.8)'

3

 

 

 

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

 

« М * {и(0>} Ц 4- РЙЧ

М х).

(41.9)'

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

Пусть и(0) — решение соответствующей линейной задачи, так называе­ мое нулевое приближение поставленной нелинейной задачи. Определяем теперь «фиктивные» массовые и поверхностные силы как выражения, стоящие в правых частях соотношений (41.8) и (41.9) соответственно. Тогда

для получения

первого приближения

нужно решить линейную задачу

для определения вектора перемещения из уравнений

 

щ - П

{и(1)} = -

рР ? =

Вц {и(0>},

(41.8)"

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

 

 

 

{и(1)} Ц+

=

М Г (1) = Мх) -

{и(0)} (41.9)"

Решение последующих приближений осуществляется с помощью введения «фиктивных» массовых и поверхностных сил по схеме

■57 П' (“<“>) -

- рЯ"> = - РР. - 4 - Ви

(41.8)"

3

^

 

а(^ П 3{и(п)}гз-

ЛТ)(п} = М х) - о $ В к} {и(п-1)} Ц. (41.9)”

Указанный метод применим для соотношений довольно

общего характера

и справедлив для нелинейной среды с анизотропией произвольного вида. Если, например, уравнения связи между напряжениями и деформациями имеют вид суммы интегралов возрастающей кратности (34.30), то для по­ лучения операторов Ри (А1Л) в этом случае нужно в соотношениях (34.30)

всюду под Ец понимать е17*, которые необходимо заменить на их выраже­ ние через вектор перемещения и и температуру Т по формуле (41.5).

Если рассматривается изотропная среда, то в этом случае тензором теплового расширения будет дельта Кронеккера, умноженная на скаляр­ ную величину а (см. § 4). Поэтому соотношения (41.5) для изотропной среды можно записать в следующем виде:

1

^

еЦ =

ег] = ~~2 ~ ( иг, з Щ , г)

и Ь> А л

(41.10)

0= 0—Заг = Щ'Ь — ъ&т.

Итак, операторы Рц (41.4) для изотропной термовязко-упругой среды получаются подстановкой в соотношения (35.17) — (35.19) выражений

еи и 0 (на которые следует заменить еи и 0 соответственно) по формулам (41.10). В качестве примера выпишем операторы Р ^ для квазилинейной

теории термовязко-упругости. Имеем из (37.4)

 

 

Рц (и) — “

Гп 06^ +

Г12

(иг,] +

щл ) ---- д-

+

 

 

ОО

 

 

СО

V

 

 

 

+

2 ^ ' гл (-з- ®) ^

“Ь 2

2

Гп р-г1 |( ге'— 2р +

1) X (41.11)

 

п=2

'

'

п=2 р=1

 

 

 

х

0"-2рер6ъ. +

 

 

[4 - К ; + «,-.*) - ~

^ - дА Л

,

где скалярные величины 0 и е определяются из выражений

 

 

0 = 1Д,* — ЗаТ,

(0~г =

0),

 

 

(41.12)

е =

 

+ ииЩл) $-(и*,к)2-

 

 

Два первых слагаемых в правой части (41.11) образуют оператор Р#,

ос­

тальные слагаемые — оператор В и (41.7).

 

 

 

В частности,

операторы

для главной квазилинейной теории термо­

вязко-упругости, квадратичной по девиаторам с мгновенной линейной упругостью, будут на основании (32.18), (32.19) иметь вид

Ра (и) =

Р%{и} +

{и},

г

 

г

 

 

 

РЬ {«} =

6гЗ 5

(* — *)

(т) +

^ (< — т) ^ [4 “ (Щ,з (т) + Щ,1 (т))—

 

е«

щ

 

 

(41.13)

 

 

I

 

 

I

в и {и} =

{$ <?[ (0Х, I -

т) 02(т) йх + $ Й (0., I - г) е (х, х) с*т} +

 

I

о

 

 

о

 

 

 

 

 

+

§<?'(&. 1 — т) 9 (*) [4 -

(“{,>(*)+ «;,»СЮ)----

 

о

 

 

 

 

где величины 0 (т) и е (т, т) определяются по формулам (41.12). Операторы 7^- и7?^{и} можно записать в сокращенной записи в следующем виде:

Р%з —

+ Н \~2Г (и1'1

иэ*1) ----§“

»

 

 

 

(41.13)'

Вц — Ьц [(?102 + &2е] +

^~2Г

и'и^ ----3~ ик>к^и\ *

Для случая квазилинейной теории термовязко-упругости несжимае­ мой среды операторы Рц имеют следующий вид (38.13):

Рц {и} =

Р% {и} + В15{и},

 

 

Р%0{и} =

1\ ^“2“ (щ,з “Ь

»

(41.14)

 

ОО

 

 

Вц {и} =

4" 2 Г2т+Гет (ии

+

и1л),

где

&— ~2~

(41.15)

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

ми выбрать уравнения (37.17)

кубичной теории вязко-упругости, то опе­

раторы

примут вид

 

 

 

 

Рц {и} = Р°ц{и} +

{и},

 

 

/'у {и} =

4" ^ в ({ ~

х) а Iм*.* (т) + и>л ООЬ

(41.16)

 

 

О

 

 

 

 

в а (и} =

**I

 

т2) т3) е (ть т2) [ии

(т3) + и}л (т3)] й х ^ х ^ .

 

- |4 ^ Г 3(г,

00 0

Если же воспользоваться уравнениями (39.17) главной кубичной теории

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

 

Рц {и} =

Р%{«} +

ВИ{и}.

 

 

 

г

 

 

 

Р%{ъ} =

4"5 Я $ ~ т) Й

(Т)Ь

(41.17)

 

о

 

 

 

в а («} — 4 " ^Гз

т) е ( 1

т) 1“*.з (т) + ии ( х)]

 

 

0

 

 

 

Соотношения

(41.16)

и (41.17)

можно объединить,

воспользовавшись

операторной записью для Р^{Щ и Вц {и},

 

Р%{и} =

4" & (ии + и;,г)>

 

(4118)

 

.

 

 

Ва {и} = 4" Гзв(«{,; + Щ,г),

где оператор Вгбудет одним и тем же для обоих случаев, а под оператором

Г3 следует понимать тройной интеграл в (41.16), если рассматривается кубичная теория вязко-упругости, и одинарный интеграл в последнем вы­ ражении (41.17), если рассматривается главная кубичная теория вязко­ упругости с мгновенной линейной упругостью.

Если добавить к соотношениям (37.17) или (39.17) еще линейный упру­

гий закон для связи шаровых ч а^ ш

тензоров напряжений и деформа­

ций

 

 

 

 

 

о = Кв,

е = 4 - а ,.

 

 

(41.19)

то уравнения (41.8) примут вид

 

 

 

— [-^ +

Г3е] (А^|) +

4

-{- Г3е)^ 9+

= — рР

где вектор

 

 

 

 

(41.20)

 

 

 

 

 

с 1 = Т3е,} [4 - (ии +

ии1) — 4 - 0^4 .

(41.21)

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

§42. Задача о расширении сферической полости

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

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

~ В щ) = ( в , +

А -7*) 0-* =

-

р Ъ -

в 1и5{и),

(42.1)

а краевые условия (41.9) —

 

 

 

 

 

 

 

 

 

 

« Я ?

я

 

( “ * . ; +

 

- I -

 

( Л х -

Л ) е } Ц +

Р Й Ч

- ^ -

=

 

 

=

м х) -

а №

в к} й

-|- 3 4 ? Л х » ^ .

 

 

 

 

(42.2)

где оператор В {й} (41.13) запишем несколько по-иному

 

 

 

В X ) { “ }

4

( “ М

” 1 ' Ц

м )

А

 

-----------1

"

------------0 3

+

 

 

+

4

- 6 ^

2

 

+

И к . г М к . О 4 - 6 . 3. { ^

— 6 &

} а г е

-

 

<?аГ (щг, +

щл) +

Оа'Ь&Т*.

 

 

 

 

(42.3)

В качестве

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

задачу о давлении

на

вязко-упругое

пространство со стороны сферы 170]. Пусть на сфере радиуса г = а дейст­ вует давление р (г), изменяющееся со временем. Предположим, на беско­ нечности давление отсутствует. Пусть температурное поле нам известно из решенной задачи о распределении температур по заданным гранич­ ным и начальным условиям. (Решение последней задачи для сферической полости в неограниченной среде имеется, например, в [87].)

Ввиду сферической

симметрии

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

компонентой вектора перемещения будет радиальная его

составляющая

и, для которой из уравнений (42.1) получается одно уравнение

| ^

( &и , 2 Зи

2и\

д, ^

 

 

+ ^

Л)\-*Г + — ф Г - -*) =

 

 

— Т<в '

йе9)+

3 / ^ 4 4 •

 

(42.4)

Граничные условия

(42.2) имеют вид

 

 

' Г Г

|у==а =

^

+ (Пг- Л) 0 4- в „ -

3Л1а71}г=а =

— р(1), (42.5)

где отличные от нуля компоненты тензора В ^ находятся из (42.3)

В„ =

 

+ {<?, - 4 & - 4- $} е* + <?,[$)’+ 2 (4-)’] +

4- [<? -

б ^ а Т в - З ^ а Т ^ . -)- Эа2^ - 2,

(42.6)

Лв9 =

=

№ ~ + {& - 4 - & - 4" С} 02 + С*[($-)’ +

+ 2 (-г)"] + И? ~ 6 <?а] а Г б - 3 $ а Т ^ - +

2<?1772.

Обозначим обращение линейного оператора

 

(г 1 + 4

г) - = П 8.

(42.7)

Пусть функция П3 (/) найдена каким-либо образом. Тогда метод упругих решений позволяет получить каждое приближение поставленной задачи в квадратурах.

Заметим, что функцию П3 (2) можно найти приближенно методом ап­ проксимаций или точно, если известно ядро ^2/з (г) для данного материала.

Для этого обозначим через 6 итерированный оператор П1Г (6 = ПХГ). Тогда оператор П3 можно представить следующим образом (см. § 20):

Пз = Й * .А = ---- Ц —

(42.8)

1 + 1Гй)

 

Если ядро ^ /з (г) неизвестно, то можно воспользоваться методом аппрокси­ маций

П3 ж и г (а + Ьсо + ск),

(42.9)

где а, Ъ, с находятся методом наименьших квадратов [28]. Таким образом, ядро П3 (I) имеет следующий вид:

I

П3(0

=

$ Пх (I — х) йвч, (т) +

§ч, (°) Пх (<) Ж аПх (*) -|-

(42.10)

 

 

0

 

 

 

 

 

 

1

 

 

 

 

 

+

^ Пх (I — т) д. [Ьсо (т) +

сп (т)] +

[Ьсо (0) +

сл (0)] Пх (*).

 

 

о

 

 

 

 

Выберем за

нулевое приближение

решение

задачи

(42.4),

(42.5) при

Вгг = Воо =

0.

Получаем неоднородное уравнение Эйлера, общее реше­

ние которого

имеет вид

 

 

 

 

 

 

г

 

 

 

и = с 1г + с %± - + зп .Л х -^ а Д Г (X, 0 ах,

 

(42.И)

а

где константы Сх и С2 определяются из граничных условий. Удовлетворяя линейным краевым условиям (42.5), получим в качестве

нулевого приближения

М(0) =

-%г[П р +

ЗП 3Пх*Та-

З П ^ П а Т * .]

+■

 

Г

 

 

+

ЗП3^ 1 - ^ - ^ 2Г (ж, 1)йх,

 

(42.12)

 

а

 

 

 

0(О) =

ЗПзЛ^Т1,

 

 

 

где Та — значения температуры на

сфере радиуса а. Если температура

в теле не меняется =

0), то соотношения (42.12)

упрощаются

«(0) = П р - |^ ,

0<о>= О,

 

(42.13)

или в полной записи

 

 

 

■*№>= т И п <‘ -

х)йр(%),

0«» = 0.

(42.13)'

При выполнении второго равенства (42.13) нелинейные операторы релак­

сации (?, (?!, ()2не будут зависеть от координат и поэтому они перестано­ вочны с операторами дифференцирования.

Подсчитываем нулевое приближение тензора В ^

в%} = в = - § - & № ) 25 -

(42.14)

Подставляя эти значения в правую часть уравнений (42.4) и в граничные условия (42.5), получим неоднородное уравнение Эйлера с правой частью

/ (г) = - Пз [ 4 В» + -Г (Яг, - Дее)] ,

(42.15)

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

Сх=

---- ^/(г)< 1г,

С2 = — -^-^г3/(г)Лг.

(42.16)

Получим решение первого приближения в виде

 

“(1) =

-^ -{ (4 - Г'1 - Г)

(ПрГ -I- р \ ~ 4 г М * (IV)2.

(42.17)

По заданному значению 0М подсчитываем нелинейные ядра (), ()1, (?2, затем по формулам (42.6) находим величину Сгг и Сее, подсчитываем / (г) (42.15) и сразу находим решение, пользуясь (42.16) и удовлетворяя изме­ ненным граничным условиям. И так далее. Отсюда видно, что любое при­ ближение поставленной задачи находится в квадратурах.

§ 43. Вязко-упругая деформация тонкостенных конструкций

Решения задач, относящихся к равновесию и деформированному состоянию вязко-упругих пластин и оболочек, так же, как и упругих, основываются на гипотезе Кирхгофа — Лява.

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

оболочек формулам

деформации

слоя, расположенного на расстоянии

%от срединной поверхности

 

 

 

 

Н"

 

 

 

(43.1)

Последнюю формулу можно записать в виде двух равенств,

одно для шаро-

вых составляющих тензоров, другое для девиаторов

 

где

 

 

е* =

е

гх,

(43.2)

1

*

 

 

 

 

х

кг] ~

Иг2

б^Х.

(43.3)

з

 

*) Всюду в этом параграфе (г, / = 1, 2).

Если оболочка достаточно тонка, так что отношением толщины ее к радиусу кривизны можно пренебречь по сравнению с единицей, то мы получим следующие формулы для определения усилий Т1}-, действующих в плоско­ сти оболочки:

ь_

2

н

(43.4)

 

и перерезывающих сил

()*:

А

 

2

 

<?{ = \ ои йъ.

(43.5)

2

Перерезывающие силы ( , несмотря на равенство нулю соответствующих им по физическому закону компонент тензора деформаций, согласно гипотезе Кирхгофа — Лява, отличны от нуля и определяются только из уравнений равновесия.

Аналогично выпишем формулы для моментов

А

АГ„=

^

 

 

 

(43-6)

 

к

 

 

 

 

 

2~

 

 

 

 

Девиатор тензора напряжений имеет в этом случае выражение

8и =

йгз<5,

,б =

"Г (а11 + бгг)1-

 

(43.7)

Для каждой компоненты 8^ из (43.7) имеем

 

 

 

2

1

522 “ “3“^22

2

1

511 — "3“ ^11

з~ ^22»

3” ^11»

$12 = $12* (43.8)

Поэтому, если мы введем дополнительно к формулам (43.4), (43.6) еще выражения для усилий и моментов в виде

 

А

 

 

2

 

8 ц =

$

(43.9)

 

к

 

 

2~~

 

 

А

 

 

2

 

Нц =

^ гз^ д.2,,

(43.10)

 

к

 

 

2

 

то получим связь усилий 8 ^ и Тн из (43.8), а моментов Н^ и Мц по следу­ ющим формулам:

== “з~(^11 2~ ^ 22) * ^11 “ “з” и ----2~ ^22^ |

-$22 = -д- ^22

2" ^ П) 1

-^22 = “3“

22------

?Г -^ 11у *

(43.11)

и обратно

тлл

2 ( ^11 + “У ^22)1

Л/ц =

2

( н 11

+ ^

22

 

 

 

# ) *

 

^*22 —2 \ 4---2” ^11^ Ж

М22 = 2 (# 22 + 4

- Яц) ,

(43.12)

Т'12 =

 

МVI

Н12'

 

 

 

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

«у (0 = $ г Ц — Г) ец (т) йх +

О

5^ /

~ 4 5 § Гз^ ’ т1’ т2’ тз)е*(т1. тг) ву (т3)

<1Хо йхъ.

(43.13)

 

%0

 

 

 

Выражение е\ (тх, т2) благодаря тому, что материал является несжи­

маемым, можно представить, учитывая (43.2), в виде [2]

 

е (тх, т2) =

е*3- (т^ е\-3(т2) =

еу (тх) еу (т2) = еу (тх) еу (т2) +

(43.14)

 

+

2 [б у (Тх) Х у (Т2)

+ 6 у (Т2) Х у (Тх)] +

22Х у (т х) Х у (Т2).

Если введем три билинейные формы е (т2, т2), х (тх, т2)] и к (тх, т2) по

формулам

 

 

 

 

 

е С*ь *2) =

в*; (тх) е*,. (та) = еп (тх) 8ц (т2) + е22 (Тх) е22 (та) +

 

 

2е12 (тх) е12 (т2),

 

 

 

X (Тх,

т2) =

х*,- (Тх) х4,- (т2) -

Хц х) Хц (т 2) 4~ х22 (Тх) х22 (т2) +

 

“Ь 2х12 (т х) х 12 (т 2),

 

 

(43.15)

* (ТЬ

Т2) =

е4,- (Тх) х 4 ,• (т2) +

8 Ъ- (т2) х 45 (Тх) = 8 ц

(Тх) Х ц (т2) +

 

 

+

в22 (тх) х22 (т2) 4 " 2 е12 (тх) х12 (т2) 4 ~ хХ (тх) Ъц

4-

 

 

Х22 (тх) 822 (т2) 4" 2х12 ( Т х) 812 (т2),

 

 

то инвариант

е* (т1? т2) можно представить в виде

 

 

е* (Тх, т2) =

е (Тх, т2) 4- т,к(т1х т2) +- 22х (тх, т2).

 

(43.16)

Подставим теперь выражение (43.13) с учетом (43.16)

в формулы (43.3) и

(43.10). После интегрирования по ъ получим для усилий

^

I

 

(43.17)

 

 

8^ = к ^ Г (2 — т) 8^ (т) йт == ^ В (/ — т)

(т),

(43.18)

о

о