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

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

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

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

Пользуясь соотношениями (13.11) и (13.15), можно найти скалярные ядра тензора ядер ползучести. Например, выделяя сингулярную часть тензора

ядер

ползучести

по

формуле

(12.5),

получим следующий

о

тензор К ^м

в случае ортотропии:

 

 

 

 

 

 

 

1

 

^21 _

Г31

 

0

0

0

 

 

Е \

 

#2

Е 3

 

 

 

 

 

 

 

 

 

 

*\Ч2

 

1

\^32

 

0

0

0

 

 

Ех

 

 

Е з

 

 

 

 

 

 

 

 

 

 

 

^13

 

^23

1

 

0

0

0

(13.16)

 

Е \

 

Е 2

Е з

 

 

 

 

 

 

 

 

 

0

 

0

0

1/С и

0

0

 

 

0

 

0

0

 

0

1/^13

0

 

 

0

 

0

0

 

0

0

1/С,

 

Здесь

Е г — так

называемые

модули

упругости;

6 и — модули сдвига;

^ — коэффициенты

Пуассона ортотропной анизотропии, между кото­

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

 

 

 

 

^12-^2 — ^21^1>

^13-^3 — ^ЗГ®Ъ

^23^3 — ^32-^2-

(13.17)

Коэффициенты, входящие в выражения (13.16), можно определить из динамических экспериментов, замеряя скорости распространения волн различных типов.

§ 14. Тензор теплового расширения

Для того чтобы учесть влияние температуры на напряженное состояние вязко-упругого тела, рассмотрим свободное деформирование этого тела под действием однородного поля температур. Если '0' — разность между температурой в рассматриваемый момент и температурой исходного со­ стояния, то деформация должна быть, очевидно, однородной и все компо­

ненты тензора деформаций обязаны быть пропорциональны

т. е.

е*; =

Ъ,

(14.1)

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

Элементы симметрии любого физического свойства кристалла должны включать элементы симметрии точечной группы кристалла.

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

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

Если в основных линейных соотношениях, данных в § 1, выберем вместо шестерки функций е ( I ) шестерку функций (14.2), то с учетом темпера­ турной изотропии получи^ обобщение соотношений (12.1) и (12.8), а именно

(

 

 

«у (0 = $ГуМ(< — *) [««К*) — <**№)] * =

 

0

 

 

1

 

 

= ^ д у« (* ~ т)

[е« (*) — *«д (*)1. .

(14-3)

О

г

 

г

 

еу(0 — <*уд (0 = $ К ш

{I — х) оы(т)йх = $ П цы (^— т)

(*)

о

о

(14.4)

В силу сформулированного выше постулата Неймана мы можем для среды с известной симметрией построить тензор теплового расширения

в явном виде.

Пусть среда является изотропной. В этом случае, мы знаем, тензор­ ный базис состоит из одного тензора Кронеккера 6*/. Поэтому для изо­ тропной среды тензор ац имеет вид

«у = «йу.

(14-5)

где а — константа,

называемая коэффициентом теплового расширения.

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

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

имеет вид

 

«у = «гГу + «2йз1йзз»

(14.6)

т. е. уже имеется два коэффициента теплового расширения аг и а 2. Поэтому связь между напряжениями и деформациями для трансвер­

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

 

 

 

I

 

*) + 2А2 (<-*)] [8ц(г) - А&(г)] +

 

би (0 = $ {[Ах(« -

 

+

0

 

 

 

 

(т)]} йх,

(I —х) [е22 (т) — ахд (г)] •+. А3 (I — тг) [е33 (т) — а2д

 

1

 

 

 

 

 

<*22(0 = $

Ц — X) [еп (тг) — ахд (т)] [Ах{I — т) + 2А2 (Ь — тг)] X

 

О

 

 

А3 {I — т) [е33 (г) — а2д (тг)]}йтг,

(14.7)

X [е22 (г) — ахд (т)] +

 

*

 

 

 

 

 

<ЗзЗ (0 ^

о

(^ ’— ^) [е11 (^) “Ь &22(тО —

(т)] -Ц [Ах (I — Т)

+

 

 

 

 

 

2Л2 {I — т)] [е33 (х) — а2й (тг)]} йх,

 

 

 

г

 

 

 

 

 

(523 (0 =

$Лб (I — т) е23 (х) йх,

 

 

' о

 

 

 

г

 

 

*

 

 

 

 

<313 (0 =

5 Лб (^

Т) 813 (Т)

°12 (*) =

5 А2 (* — *) 212 (Т) йх.

О

О

Для ортотропной среды тензорный базис состоит из трех тензоров (13.12).

Поэтому тензор теплового расширения а ^ имеет вид

а« = а ^ у + а2Р»; + а3Г«,

,

(14.8)

а связь между напряжениями и деформациями будет выглядеть так:

х

«и (0 =

5 {Аа Ц — т) [8ц (т) — «!<> (т)] +

Л4 X) [е22 (т) —

 

О

 

(*)] + Ав (* — т) [е33 (т) — ос3# (т)]}Ах,

 

х

 

0*2 (*) =

5 (Л1 (* т) 1®и (Т) — «1* (*)] +

Л2(* — X) [822(X) —

 

о

 

«гй (*)] + л8(* — х) [833 (т) — а3# (*)]} Ах,

X

«33 (<) =

$ (Лв (* -г) [8ц (т) —

(-г)] + Л8 {г — -г) [е22(х) —

 

О

[е33 О») — Яз®(^)1}

 

(^)] “Ь ^3

 

X

 

 

«23(0 =

$Л9 (г — т) е23(т) Ах,

 

 

 

О

 

 

 

X

 

 

«13 (0 = $ л 7 {I х) 813 (г) Ах,

 

 

 

О

 

 

«12 (*) =

X

 

 

^ А6 (« — г) е12 (т) йт.

 

(14.9)

 

О

 

 

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

кций релаксации Где в (12.10) можно представить в виде

Т„м (<) = Гун Г (0 + Гун г х (*),

(14.10)

где Где , Где — тензоры-константы структуры, определяемые, как и Г в (12.10), из опытов и обладающие свойствами симметрии (12.5), (12.6);

г (0 = — Д' ({), г, (г) = — К' (0 .

Упрощения (14.10) для частных видов анизотропии в § 13 уже указа­ ны. Например, для ортотропных стеклопластиков будут справедливы со­ отношения (13.15), в которых надо положить для любого к = 1,2, ..., 9:

А* (*) =

К

б (*) —

Г

(0 -

Гг (<),

(14.11)

причем А* .

%1

%1 (к — 1,2,

..., 9)— упругие постоянные

структуры.

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

§ 15. Основные типы задач линейной теории вязко-упругоети

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

Ои,] + Р^г =

(15.1)

где щ — вектор перемещения, — массовые силы; индексы, стоящие по­ сле запятой, означают дифференцирование по соответствующим коор­ динатам

д

8 . -

(*\ /, к, 1= 1, 2, 3).

дх,

Чъ ы — дхидхг

 

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

оги +

= 0.

(15.2)

Связь между тензором деформации и вектором перемещений, как и в линейной теории упругости, задается соотношениями Коши

=

(^1,7 ~Ь ^7.г)/2,

 

(15.3)

причем для

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

этой системы

шести дифференциальных

уравнений (г, / = 1,2,3) для щ при заданных

ен необходимо выполнение

условий совместности Сен-Венана

 

1^(г7) __

€^* п ^>1т^кп ~

(15.4)

где в1к1 — альтернирующий тензор.

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

Пусть поверхность 2 ,

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

двух частей 2 Хи Е2 и пусть на части поверхности 2 Хзаданы нагрузки

°ц1}121 = *$чо.

(15.5)

где 1$ — косинусы углов между нормалью к поверхности 2 и /-й коор­ динатной осью. На части поверхности 2 2 заданы перемещения щ0

Более общими граничными условиями будут так называемые «контактные» граничные условия, имеющие следующий вид:

с4т ОтзЬ +

Р $ Ите Е Л =

 

 

(15.7)

где Еус — модули

Юнга

контактирующих

упругих*

тел; оь*т,

постоянные величины;

I — характерный

линейный

размер;

N4?

так называемые «контактирующие» силы; х — индекс площадки, причем поверхность 2 состоит из суммы всех таких 2 Х.

Нетрудно видеть, что условия (15.5) и (15.6) являются частным слу­ чаем условия (15.7), когда х принимает значения 1 и 2, причем

«»т =

в«т,

Р " =

О,

М ” =

^гО,

<х$ =

о,

РШ =

б4т,

Я ?Ц Е %=

и,ю.

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

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

1.Задачи, в которых можно пренебречь изменением механических свойств материала (эти свойства считаются определенными заданием ядер ползучести и релаксации) за счет изменения температур.

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

3.Задачи, в которых также учитывается существенное влияние тем­

пературы на механические свойства

материала, но температурное поле

в теле является нестационарным и

неоднородным.

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

 

*

 

 

 

 

* « (« )= $ я ( * - * ) * « (* ) ,

(15.8)

 

О

 

 

 

 

 

И

 

 

 

 

 

<3(0 = ^ 1 ( « — Т)с/0т(<),

 

 

О

 

 

 

 

 

0 Г = 0 — Зад,

0 =

у и,

(15.9)

где д =

Т — Г0 — температура, представленная по шкале,

в которой

имеет место условие: при

I =

0, д = 0, т. е. в начальный момент темпе­

ратура

тела постоянна

=

Г0).

 

Соотношения (15.8) и (15.9) можно обратить, т. е. выразить деформации

через напряжения

 

*«(0 =

5 п

(15.10)

 

О

 

 

 

I

 

 

8т (0 =

5

^ ( 4

т)Ла (т).

 

о

 

 

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

§ 16. Преобразование Лапласа —Карсона

Для решения задач первого и второго типов линейной теории используем преобразование Лапласа — Карсона [26] и предпосылки метода аппрок­ симаций [28], [35]. Это преобразование заключается в том, что всякой функции / (I), определенной на интервале 0 ^ I < оо и называемой в дальнейшем оригиналом, ставится в соответствие функция /* (р) комп­ лексной переменной р по следующему закону:

оо

 

/*(р) = />$ е-Р*/у)(И.

(16.1)

О

 

Функция /* (р) носит название

изображения функции / (0* Как видно

из (16.1), преобразование Лапласа — Карсона является линейным.

Если известна функция /* (р),

то / (2) находится с помощью обратного

преобразования Мелина

 

с-Ноо ^

 

/ М - - Я Г - $

(46.2)

О—200

 

причем интегрирование ведется в комплексной плоскости вдоль вертикаль­ ной прямой с асбциссой о.

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

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

между функциями времени

/ (г), § {I) и ф (I) в

следующем виде:

г

 

 

Ф( *) = $ /(« - О * *

(О-

(16-3>

О

 

 

От левой и правой частей соотношения (16.3) возьмем интеграл вида (16.1). Получим соотношение в изображениях

ооI

Ф* (Р) = Р$ ерШ ^/(* — О (О*

Оо

Поменяв местами порядок интегрирования в соотношениях (в предполо­ жении, что функции / и д допускают такую замену) и вводя новую пере­ менную под знаком интеграла х — I т, получим

ОО оо оо оо

•ф*(р) — р \ ^8 (О ^ / (* —О е~р1сИ=

^ в.§ (т) е-Р*р^ / (х) е~рхйх —

о

о

о

 

 

ОО

 

оо

 

= 1 * (р ) 5 е~рх^8 (О = /*(р ) (°) + Р $ е~рх8

 

О

**

О

 

Принимая, что при I — 0 функция § равна нулю, получаем изображение (16.3)

У(Р)=Г(Р)8*(Р)-

(16.4)

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

Применим формулы (16.3) и (16.4) к соотношениям (15.8) — (15.11). Получим в изображениях связь напряжений с деформациями

*Ь (р ) = •к * (р) еЬ (р )> а* (р) = Я* (р) [0* (р) — Зад* (р)],

еЬ (р) = П » зЬ (р),

0* (Р) ~ 3 ад* (р) = Пх* (р) а* (р).

(16.5)

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

Я> Ср) П* (р) = 1,

(р) Пх (р) = 1.

(16.6)

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

ной теории упругости удвоенный модуль сдвига 26?, а величине Яг (р) — модуль всестороннего сжатия К. Чтобы сделать аналогию более полной, будем наряду с уже введенными обозначениями применять также следую­ щие:

К 9 = 26?*,

Я[ = К \

(16.7)

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

. . З К ' - 2 С * _ З Д ;~ Д *

(16.8)

г

6К* + 20*

6Д* + Д* ’

 

Я* =

К* -

2С*/3 = Дх* - Д*/3,

 

[X* =

С* =

Л 72;

 

(16.9)

... _

9Д*С* _

9д;д*

(16.10)

 

зя* + е*

бд*+ д* ’

 

 

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

Отметим важную особенность оператора Пуассона V*. В случае, если оператор Пуассона есть число V* = V , то, как видно из (16.8), операторы

Я* и Яг подобны и их отношение представляет собой число

1 + у

ю*

3 (1 — 2V)

*

И обратно, если между операторами Кг и Л* существует зависимость = аК*, то оператор V* есть число. В самом деле, из (16.8) следует

»За — 6а + 1

Соотношения (15.8) — (15.11) можно представить в другом виде, согласно предыдущим главам. А именно

 

г

*

 

Ъ) (0 = 5 Г {I т) (т) йх = 20ец (I) — $ Г (< — т) Дег}(т),

 

0

о

 

 

1

 

 

о («) = 51\ {I — т) [0 (т) — 3x0 (г)] ах =

К [0 (0 — 3x0

 

0

 

 

 

I

 

 

— т) [0 (т) — 3x0 (т)] а%

 

(16.12)

 

0

 

 

и обратно

г

г

 

 

 

ец (0 = {\К(1 — т)81)(т)ах = -1^ 8 ц (<) +

х)81}(х)ах,

 

о

о

г

 

г

 

0 (*) — Зад (I) = ^ К г(Ьт) а (т) йх =

 

а (г) + ^ Кг— т) в (т)йх.

 

о

 

о

 

 

 

(16.13)

Здесь, как и ранее, ядра с волной вверху означают регулярную часть со­ ответствующих ядер без волны.

Для применения преобразования Лапласа — Карсона к соотношениям (16.12), (16.13) проделаем с соотношением

 

г

 

■Ф(0 =

§/(* — ^)в (х) Лх

(16.14)

 

О

 

(где / и # — регулярные функции, причем § (0) =

0) те же выкладки, ко­

торые были проведены для (16.3). Тогда получим в изображениях

'ГСр) =

- ^ - П р )?*(р ),

(16.15)

Применив формулы (16.14) и (16.15) к соотношениям (16.12) — (16.13), получим в изображениях

ец {р) = \ 2 С ~ ^ М Л ^ { р ) ,

(16.16)

о* (р) = - ^ - ] [0*(р) - 3x0* (р)]

и обратно

1

К* (Р)

Разрешая уравнения (16.16) относительно еу (р), получим

4 ( Р ) =

------ ------------- = - щ - * ч { р )+

20 * р-------4 (Р ). (16.18)

 

 

20-

— Г (р)

2С - — г (р)

Сравнивая (16.17) и (16.18) и обозначив

 

 

1

Г Ф( Р)

-

<2*(р),

 

(16.19)

г20

р

 

имеем

 

 

 

 

 

2 0 Т ( р )

^

(?*(/>)

 

(16.20)

 

Р

 

1 - <2* (р)

 

 

 

 

 

Разложим выражение (16.20) в ряд по степеням <2* (р); получим в ориги­ налах следующее выражение:

ооI

еа (0 = -дг *« ( 0 + 2

Г(п)(< - т) «у (т) <?т,

(16.21)

п = 0

0

 

где

г (0)(() = г (о ,

I

Г(1) (<) = 5 г(0)(< —т)Г(т)«?т,

(16.22)

Г(п) (0 = ^ Г(п_1) (* — т) Г (т) <7т.

о

Как видно из соотношений (16.19) и (16.21), в разложении правой части (16.20) в ряд величину 1/2 С можно использовать как малый параметр.

§17.

Постановка квазистатических задач первого и второго типов

 

и общее решение в изображениях

Будем считать, что температурное поле в теле задано однородным, это существенно только при учете влияния его на ядра ползучести П и ре­ лаксации 7?, и нестационарным, т. е. температура Т ~ Г0 + д задана функцией истинного времени I. Приведенное по — ^-аналогии время

I

I' = ^ 6.1)а т

о

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

точках тела. Следовательно, задаваемые на поверхности силы *$^0, Л^х) и перемещения к*0, а также объемные силы рР г и температуру'6 можно рас­ сматривать, как заданные функции координат и приведенного времени Понятно, что любая функция истинного времени ф (I) есть некоторая дру­ гая функция ф(0> но, учитывая обычное правило сохранять обозначе­ ние функции при замене независимой переменной, мы пишем ф (I') вместо

* (О-

Обращаем на это внимание: во всех задачах второго типа внешние силы и перемещения считаются заданными функциями приведенного времени I'; кроме того, штрих сверху писать не будем (вместо I' всегда пишем 2).

Таким образом, в математической записи всех соотношений, в которых

не используется конкретный вид заданных функций времени 5 ^, мг0, Ф, задачи первого и второго типов неразличимы, т. е. решение какой-

нибудь задачи первого типа всегда есть решение некоторой соответствую­

щей задачи второго типа. Постановки и общие решения

таких задач оди­

наковы.

 

 

е^

и скалярами а,

Соотношения Больцмана между девиаторами

0 т = 0 — ЗаФ в изображениях (16.5)

 

 

 

=

4 = П*4-,

Д *П* = П*Д* =

1,

(17.1)

 

 

 

 

о* =

е; = п 1 4

я[ п* = п Ж = 1

 

вместе с соотношениями Коши для изображений компонент деформаций

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

и*

 

 

 

2бг; ===

г >

 

 

_| 0 ^г^/3,

 

0* = и*, 1,

 

4* =

4

— 6**и/3

(17•2)

замыкают систему дифференциальных уравнений равновесия

 

$ “Ь Р ^г =

0,

бу =

+

3*6^-.

(17 *3)

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

щихся исключением в (17.3) оц и Ец на основе (17.1) и (17.2):

 

— ^ е! * + «I, П+ 2р1ГТ?; = 0.

(17.4)

Здесь множитель — оператор, содержащий «коэффициент Пуассона» V%

выражается через основные операторы, входящие в (17.1)

 

5 - 1 ^ = 7 3 + 2 0 ^ ;.

(17.5)

Если задачу мы хотим решать в напряжениях, то уравнения (17.3) замыкаются шестью условиями совместности деформаций Сен-Венана>

вытекающими из

(17.2)

 

0,1‘Э= —

кк Н" &гк9дк &]к, гкч

(17.Ь)

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

еу = П* [о4, -

.

 

(17.7)

Следовательно, в этом случае уравнения (17.4)

замыкаются шестью урав­

нениями совместности Бельтрами

[27]

 

 

бгл кк + | 1 ^ * б, г] +

 

Р и Г$к] + Р и

= 0.

(17.8)