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

книги / Численные методы в теории упругости и пластичности

..pdf
Скачиваний:
5
Добавлен:
20.11.2023
Размер:
10.57 Mб
Скачать

Тогда можно ввести две теплоемкости вещества: ср и с„ (кал/г ■ град), причем, как следует из (2.31) и (2.33),

 

дц_

_ т ^ ф _

а я

 

Су — дт

 

ат2 ~

дт’

(5.20)

дц 0

= -т а2Ф0 = т

ая0

(5.21)

 

дТ

 

я г 2

а г

 

Поэтому

т

(5.22)

То

По закону Дюлонга-Пти у твердых веществ для температуры выше некоторой, называемой температурой Лебая, ср = сопз1. В втом случае

ср 1п Т_

(5.23)

То

Бели а мало по сравнению с То, то

Г0

^

(524)

Сравнивая (5.16) и (5.17) и считая теплоемкость вещества за­ данной, получаем, что функция энтропии Я полностью определена законом связи между напряжениями и деформациями:

рН = рср 1п Т_

(5.25)

П + ?•?{§ —9^}-

 

Поэтому уравнение притока тепла (5.2) примет вид

 

рерТ = сЦу(Ат 8гас!Г) - Т [аГ {е - а#}] + Р9 + УГ.

(5.26)

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

<гц = /

- г)й[ец(т) - аы»?(т)],

(5 27)

о

 

 

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

(? =

1 ,...,^ )

выражаются

через ет по закону

 

 

 

 

<

 

 

 

 

р(?) = I М

$ , Ц - т)4[ек1(т) -

аыд(т)],

(5.28)

о

 

 

 

 

функция рассеивания выражается в виде [72, с. 229]

 

IV'

т•

яе1

 

(5.29)

 

 

 

 

У пражнение 5.6. Показать, что при ТУ = 1 из предыдущего упражнения следует, что

И" =

- ^ПОы(0)[<г0 «гы] ,

(5 30)

где Пцы(0 — тензор ядер, обратный к тензору Дуы(0> т е-

/ ^ ы {1 - г)Л1ытп(г) = Дут „ - % ы (0П м тп (0).

(5.31)

о

Упражнение 5.7. Пусть для изотропной квазилинейной среды с определяющими уравнениями в виде

^

*

 

50 = ^ К Ц -т ) <Те0-(т) + ^ Т(1 —

, ет)еу (т)Ат,

\

;

(5-32)

«г = У Ях(* -

т) А9(т) + ^ Г\(< - г, 0т,ет)9т(т) Ат,

о

о

 

где Г и Г[ являются функционалами от величин 0? = 9(т) — Зо^(т) и е = еу(т1)е0 (т2), причем Г(0,аг,у) = Гх(0,х,у) = 0. Пусть тензоры связаны с ет соотношениями, аналогичными (5.32), но с другими ядрами. Доказать, что в этом случае функция рассеивания имеет вид [72, с. 231]

 

N

 

Ж* = 5,-,-е0- + <т(9 - 3«*) -

+ 3 ^ » У ?) ].

(5.33)

1=1

И(я) = \ $ 6 ц ; (№ ? = ( $ } - и(я)6ц)(№ ~

<5 34)

Упражнение 5.8. Показать, что для предыдущего упражнения в случае N = 1

IV = в»,-в” + <г(9 - За»?)' - П(0)<г„<ги - Пх(0)<г<г',

(5.35)

где П(<) и Щ(1) — ядра, резольвентные к Щ1) и Дх(0 соответ­ ственно.

Упражнение 5.9. Показать, что если связь между напряжени­ ями и деформациями описывается соотношениями теории малых упругопластических деформаций (4.47), (4.35) и N = 1, то

V/* 8,у[»уи;(еи)] = <ти[шеи]

Ш

(5.36)

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

<т= К (в -Ъ а д ),

5 у = Л Д е),

(5.37)

и уравнение притока тепла (5.26)

приобретает вид

 

рсрТ = (АТГД,- -

ЗГа<г + р я + IV*.

(5.38)

Бели среда однородная, то

 

 

 

рсрТ = Ат Д Г -

ЗТа* +Д4 -ИУ*,

(5.39)

где А — оператор Лапласа.

Используя первое из соотношений (5.37), уравнение (5.39) мож­

но записать в виде

 

рсрТ = Ат АТ - 3аТ9 +ря + IV*,

(5.40)

где

 

рс„ = рср - 9а2КТ.

(5.41)

Бели среда обратимая, то IV* = 0. Бели, кроме того, величина не изменяется со временем, то получим из (5.39) линейное уравне­ ние теплопроводности, замкнутое относительно температуры:

рсрТ = ХтАТ + рд.

(5.42)

Заметим, что из (5.41) следует, что обе теплоемкости ср и с„ не могут быть одновременно постоянными (не зависящими от тем­ пературы). Поэтому очень часто принимается допущение о том, чго во втором слагаемом правой части уравнения (5.26) темпе­ ратуру Т можно заменить на температуру ТоТогда уравнение (5.26) принимает вид

рср = ут §гас1Т) — % [а Р {е — ад?}]1 + р^ + VV*,

(5.43)

а соотношение (5.41) —

рс„ = рср 9 а 2КТо.

(5.44)

§ 6. КЛАССИФИКАЦИЯ И ПОСТАН ОВКА ЗАДАЧ М Д ТТ

При изотермических процессах мы можем получить замкнутую систему уравнений движения, исходя из определяющих соотно­ шений (4.1). В самом деле, подставляя их в уравнение (2.9) и воспользовавшись уравнениями Коши (1.1) или (1.9), получим систему трех уравнений относительно трех неизвестных и,-:

/>«" = р Р + ЕЙу ^ {Б е Г и }.

(6 .1)

Будем выражение ^"{БеГ3 } сокращенно записывать в виде

<г{«} = ^ {Б еГ н }.

(6.2)

Таким образом, запись <т{й} означает, что проделаны все выклад­ ки с тем, чтобы напряжения выразить через перемещения. Точно также запись &{€) означает, что вместо напряжений по опреде­ ляющим соотношениям подставлены деформации. Аналогичным образом следует трактовать и запись | {?} (деформации по закону (4.60) выражены через напряжения). Тогда три уравнения (6.1) в компонентах можно записать в следующем виде:

РЩ = рРг + &Ц,] { « } .

(6.3)

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

[аО?Ч * { « } п* + ь,-1 Ч ь , = Ч ?)>

(6-4)

где а ^ , 6^ — некоторые положительно определенные матрицы,

— вектор контактных усилий. Мы будем в большинстве случаев рассматривать частный случай граничных условий (6.4), а именно: на части границы Ех будут заданы перемещения 3°, а на части Е 2 — нагрузки (усилия) 5 0:

«. |е ,= «®(х , 0 ; * ч {« } пУ |е ,= 5?(2,<).

(65)

Для классификации граничных условий принята такая терми­ нология. Если на границе заданы только перемещения (только кинематические условия), т.е. в (6.5) Е 2 = 0, то такие условия на­ зываются граничными условиями первого рода, а задача М Л ТТ , использующая эти условия, — первой краевой задачей. Если на границе заданы только усилия (только статические условия), т.е. в (6.5) Ех = 0, то граничные условия называются граничными условиями второго рода, а соответствующая задача М Л Т Т — второй краевой задачей. Условия (6.5) называются смешанными граничными условиями, а задача М Л Т Т , их использующая, — смешанной краевой задачей М Л ТТ . В некоторый момент I = 1а должны быть заданы и начальные условия. Например,

при ^ = ^о «,• = 1/*(х), «< = Ц(х).

(6.6)

Если тело не ограничено, то должны еще быть заданы условия на бесконечности. Итак, изотермическая задача М Л Т Т заключается в решении трех уравнений (6.3) с граничными условиями (6.4) или (6.5) и начальными данными (6.6). В этом и состоит постановка динамической задачи М Л Т Т в перемещениях. Сюда, разумеется, следует добавить требование гладкости разыскиваемого решения, границы и «входных данных»: рР{, 5 °, и®, Щ и У*.

Статическая задача М Л Т Т заключается в решении уравнений равновесия в перемещениях

«И.ЛЩ + р Ь = 0

(6.7)

при выполнении граничных условий (6.4) или (6.5).

Если определяющие соотношения или входные данные заданы так, что решение может зависеть от времени, то задача (6.7),(6.5) называется квазистатической задачей М Л Т Т .

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

Заметим, что решение динамической задачи ид (х, I) иногда может стремиться к решению статической задачи ис(х,<):

ие(х) = Иш Йд(г,<).

(6-8)

I—►ОО

 

Это явленйе связано с диссипацией, или рассеиванием энергии. Условие (6.8) мажет и не выполняться, как, например, в идеальных упругих средах. В таком случае статическая задача М Л Т Т является самостоятельной задачей о равновесии среды и никак не связана с динамической задачей, в которой рассматривается распространение волн и кинетической энергией среды пренебречь нельзя.

Л ля второй краевой задачи М Л Т Т , когда граничные условия имеют вид

*п |Ез= 5°(5,<). <гцп5 Ь ,= 5?0М).

(6.9)

можно дать постановку задачи в напряжениях (для одно.связной области V). Л ля этого в уравнения совместности (1.14) требуется подставить определяющие соотношения в форме (4.51):

1пк(/{<г} = 0, еш е,-тп**п,!т{е} =

(6Л0)

Тогда статическая (или квазистатическая) задача М Д ТТ заклю­ чается в решении системы уравнений (6.10) при удовлетворении уравнений равновесия в напряжениях

Югу <т + рР = 0, <гу,^ -I- рР% = 0

(6-1 1 )

и граничных условий (6.9).

Если массовые силы являются градиентом некоторой скаляр­

ной функции х> так что

 

р Р = %таАх, лР« = Х,»,

(6.12)

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

= 1пк Ф + х -[, &1] = *Нс1С]тп ф *„ ,1т + Х^>] •

(6.13)

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

напряжения заключается в следующем. Требуется найти реше­ ние системы из шести уравнений

1пк 0{1пк Ф + х 1 } = 0

(6.14)

относительно шести неизвестных Ф,; при удовлетворении гранич­ ным условиям

1пкф .п|Ег= 5 0 -х й | Е 2 .

(6.15)

Д ля динамической задачи М Л Т Т также можно дать постанов­ ку в напряжениях. Для этого применим к уравнениям движения

(2.9) операцию БеГ (1.11):

 

ре" = р БеГ(Г) + БеГБпг <т.

(6.16)

Подставляя сюда определяющие соотношения (4.60), получим сис­ тему шести уравнений

р[Ф{<т}]" = БеГ(р^) + БеГ Б1у <т .

(6.17)

Сформулируем теперь начальные условия. Применяя к ус­ ловиям (6.6) операцию БеГ и воспользовавшись определяющими соотношениями (4.60), получим

=

И ? }]

(6-18)

Тогда динамическая задача М Д ТТ состоит в решении системы уравнений (6.10), (6.17) при удовлетворении граничным условиям (6.9) и начальным данным (6.18).

В конце § 5 мы столкнулись с ситуацией, когда требуется решить «чистое» уравнение теплопроводности (5.43) без механи­ ческих членов

рСрТ = ХтА Т + рд.

(6.19)

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

РСрТ — (А ^ТД , + рд.

(6.20)

Для уравнений (6.19) или (6.20) также должны быть начальные данные:

при 1 — 0 Т = Г0.

(6.21)

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

Ат Д Г + р« = 0,

(6.22)

или

 

(\ Ъ Т ;) + рд = 0.

(6.23)

Если источники тепла отсутствуют, то в (6.22) температура стал новится гармонической функцией

А Т = 0.

(6.24)

На поверхности тела может быть задана температура (гранич­ ное условие первого рода, или условие типа Дирихле):

Г Ь = Г °(5 ,*).

(6.25)

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

\$Т;тц \ ъ = - Я°(х,1).

(6.26)

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

А $ 7 > , |2 = /)(Т - Т °) |Е,

(6.27)

где Т° — заданная температура окружающей среды, а /? — ко­ эффициент теплоотдачи (кал/см2 •с •град). М огут встретиться и более сложные (нелинейные) граничные условия, например для случал теплопередачи излучением. Однако мы будем рассматри­ вать только линейные граничные условия. Их можно объединить следующим образом. Пусть на части Е , поверхности Е

 

с^ Т .п ,- + Ь ^ Т = Т °(?),

(6.28)

где с<?>,

— некоторые величины, зависящие, вообще говоря,

от координат

и времени.

 

Итак, нестационарная задача теплопроводности заключается в интегрировании уравнения (6.20) при выполнении начального условия (6.2 1) и граничного условия (6.28) или условий на бес­ конечности.

Стационарная задача теплопроводности заключается в ин­ тегрировании уравнения (6.24) при удовлетворении граничного условия (6.28).

Разумеется, и здесь необходимо добавить требования, пре­ дъявляемые к гладкости разыскиваемого решения, границы и входных данных.

Как уже было отмечено здесь, решение нестационарной задачи Тн(х,1) связано с решением стационарной задачи Тс(х) условием

Тс(2 )= Н т Т .(2 ,0 .

(6.29)

1-^00

 

Если входные данные стационарной задачи зависят от времени как от параметра, то задача (6.23), (6.28) называется квазистационарной. Если Хт явно зависит от координат, то стационарная, квазистационарная или нестационарная задача называется неод­ нородной. В противном случае она будет однородной.

Если рассматривается изотермическое деформирование, то в уравнения движения (2.9) следует подставить вместо определяю­ щих соотношений (4.1) уравнения (5.7) или (5.8). Тогда получим уравнения движения

РЩ = рРг + агцл{и, Т )

(6.30)

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

<гц>;{ и ,Т } + рГй = 0

(6.31)

и граничные условия типа (6.4) или (6.5):

-*{й, т } » * + б (;Ч -Ь , = лг^О М ),

(6-32)

и»|е 1 = и»(*>^)!

{ и>-О пУ |ез = 5^*(*|4).

(6.33)

Если в уравнение притока тепла не входят механические чле­ ны, т.е. оно имеет вид (6.20) или (6.23), то его можно решить отдельно, найти температурное поле, а затем уже решать уравне­ ния (6.30) или (6.31), считая температуру известной. Такие задачи называются несвязанными задачами термомеханики. Например, несвязанная динамическая нестационарная задача термомеханики заключается сначала в решении уравнений (6.20) при выполнении

граничных условий (6.28) и начальных данных (6.2 1), а затем в ре­ шении уравнений (6.30) при выполнении граничных условий (6.32) и начальных данных (6.6).

Упражнение 6.1 . Д ать постановку несвязанных задач тер­ момеханики в перемещениях: статической (квазистатической) нестационарной, статической (квазистатической) стационарной (нестационарной).

Упражнение 6.2. Дать постановку указанных в упражнении 6.1 несвязанных задач термомеханики в напряжениях. ■

Если в уравнение притока тепла включены и механические члены, то его нужно решать совместно с уравнениями (6.30) или (6.31). Такая задача называется связанной задачей термомехани­ ки. У равнение притока тепла (5.26) можно переписать в виде

рсрТ = (АТ.ТД, - Г [ а

, {3, Т}\ + р 9 + Ж * .

(6.34)

При этом конкретное выражение

в виде зависимости от и и Т

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

ворении граничным условиям (6.32), (6.28) и начальным данным

(6.6), (6.21).

Упражнение 6.3. Описать постановку связанных задач термо­ механики в перемещениях: статической (квазистатической) неста­ ционарной и статической (квазистатической) стационарной (квазистационарной).

Упражнение 6.4. Описать постановку указанных в упражне­ нии 6.3 связанных задач термомеханики в напряжениях. ■

Заметим, что если оператор Т определяющих соотношений является линейным, то несвязанная задача термомехаиики после определения температурного поля может быть сведена к задаче М Д Т Т при изотермическом деформировании. В самом деле, в силу линейности оператора Т из (5.8) и (5.9) следует

<г = ? { е }

 

(6.35)

) в (б.ЗО)-(б.ЗЗ), имеем соответственно

 

/Ж| = рР- +

 

(6.36)

=

0,

(6.37)

 

 

(6.38)

. 1 = « ? (*,< ).

з?°(г »0 -

(6.39)

50

Соседние файлы в папке книги