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

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

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

интегралу в (1.16) формулу (1.8) и произведем интегрирование по частям:

М

М

 

 

! Ь>Ц <%] =

& \м0 - I $•»«,»

(11?)

м0

Л/п

 

 

Мо

 

 

Тогда получим

мм

«. = + “ о *; - «2»*® + ^ е^<цт - ^ Ььщ,т <%т-

(1-18)

Мо Мо

Добавляя и вычитая в правой части (1.18) выражение и^х,- и используя тождество (1.14), получим

м

Щ= «< + <*$(*7 - * “) + ^ [с,ж + (х; -

- ет ;,.)] <^т. (119)

Мо

Э т и формулы известны под названием формул Чезаро [62]. Они

позволяют для односвязной области однозначно выразить переме­ щения по заданным деформациям в данной точке М .

Выберем теперь замкнутый контур, т.е. предположим, что точ­ ка М совпадает в (1.19) с точкой Мо- Тогда интеграл в (1.19) должен для односвязной области обратиться в ноль. Применим к нему теорему Стокса о роторе [84]:

^

^ ^рдт^гт^^р 6^2.

(1.20)

 

Е

 

Тогда получим, используя формулу (1.15),

 

Срдт ^ \&тд,1 + (*/

^)(пк1^Ы]^тк,д!]^р —- 0.

(1.21)

2

 

 

достаточное условие обращения в нуль интеграла (1.21) имеет

мид

 

Т] = 1пке = 0, Т}рп = (рдт^пк1^тк,д1 — 0>

(1-22)

I (

должен обращаться в нуль компонент некоторого симметрич-

.....

тензора второго ранга Т), который носит название тензора

несовместности. Уравнения (1.22) называются условиями совмес­ тности. В конце настоящего параграфа будет доказано, что для односвязной области они являются необходимыми и достаточными условиями интегрируемости уравнений (1.9). Для многосвязной области они являются только необходимыми.

Упражнение 1.2. Доказать, что для многосвлзной области

V, в которой можно провести поверхности

8 а , а =

1,2,...,<*,

так, чтобы сделать область V односвязной, для интегрируемости

уравнений (1.9) требуются дополнительные условия

 

■+■(х) ~ ^})(пк!еп>}^тк,1]

= 0,

(1-23)

У СЦк^т.к <%т = 0,

 

(1.24)

Г„

 

 

где Га (а = 1 ,2 ,...,9) — замкнутые контуры в V, причем каждый из них пересекает только одну поверхность [62, с. 491].

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

«,- = ц? + (ы?- + с0 )(х,- - х°).

(1.25)

Упражнение 1.4. Показать, что если в условиях упражнения 1.3 содержится предположение о «жестком» защемлении точки Мо, т.е. при X] — х°

щ = 0, ш4 = т^цк^к,] = 0,

(1.26)

то

 

и»- — с,-,(х,- —х^ . И

(1.27)

Мы будем часто пользоваться разложением тензора деформа­ ции на девиаторную и шаровую части:

€ = 5 +

а * = сч +

1

(1-28)

где 0 — первый инвариант тензора деформации, описывающий изменение объема среды:

0 = (е) = еи = <Иу «.

(1.29)

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

' * /(®5) —

(1.30)

Тензор

а*

е ** =

(1.31)

назовем тензором скорости! деформации. Из соотношений (1.9) и (1.4) следует, что

«'« =

+ *!.<)•

(1.32)

Упражнение 1.5. Используя формулу (1У.15) приложения ГУ, показать, что тензор несовместности т) макет быть записан в виде

ЧУ — Сгк1(»тп€кп,1т

(1.33)

6 { } ( е Пт ,пт А -0) +

+ Дс,4^ &{к,к] ~- ^ к .к г -

Упражнение 1.6. Доказать, что тензор несовместности Т) мож­ но разбить на шаровую составляющую и девиатор г} следующим

образом:

1

 

 

Чу =

1 с

 

(1.34)

з*7*»у + Чу>

где

 

 

 

1) = щ б ц

= €к1,к1 -

Д0,

(1.35)

9у = ^.У + Д«У ~ € *к,к] + е ]к,к>

+ §(ек!,к1 - Д в ) 6 ц .

(1.36)

Упражнение 1.7. Доказать, что если справедливы уравнения

совместности (1.22)

 

 

 

ЧУ = 0,

(1.37)

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

 

 

 

= Д0

(1.38)

и

 

 

0,х] +

^гк,к] + €]к,к>-

(1.39)

Упражяея ие 1.8. Доказать, что уравнения (1.37) и (1.39) экви­ валентны, т.е. если справедливы одни, то справедливы и другие.

Упражнение 1.9. Доказать, что если выполняется хотя бы одно из условий — (1-37) или (1.39), то справедливы и условия (1.38). ■

Очевидно, что если выполняются условия совместности (1.37), то выполняются условия обращения в нуль всех компонент сим­ метричного тензора второго ранга Я :

Яу = Аеу +

- с и м ~ *}>=№+ &](ек11 Д0) = 0,

(1-40)

где ( — произвольный симметричный тензор-константа второго ранга. Очевидно, что

Я,-у =

щ

при (,у — у,

(1.41)

Я<у н

щ

при (у = |«<у.

(1.42)

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

В самом деле, каждое тензорное поле г (необязательно сим­ метричное), обращающееся на бесконечности в нуль, можно од­ нозначно представить в виде [113]

т = V ® а + V x А , Ту = а,-,,- + иыАц,к-

(143)

Но в силу той же теоремы такое представление справедливо и для тензора А:

А — 6® V + В х V,

А Ь]^ + сутп В{т)П

(1.44)

Подставляя (1.44) в (1.43), получим

 

 

Ту — Оу,« ■+■С<у -(- €,'ц€утп В|т,пк!

(1.45)

где

 

 

 

(1.46)

С = V X Ь,

 

Суу = Ик1^1^к-

Разобьем теперь тензоры т,

В

на симметричные и антисиммет-

ричные составляющие

 

 

 

 

Ту —с,у + ?у,

 

—су,-,

— ?у»'

(1-47)

В ц — •З’у + ^»у,

= 5 ц >

=

(1-48)

и, кроме того,

положим

 

 

 

а* •+■с,-.= и,-, а, —с,- =

(1.49)

 

 

 

Тогда из

(1.45)

имеем

 

<г,у+9*у =

^(и>1<+»у|,)+ ^ (« ,-,;-« о ')+с,к' с>тп^ ,т >,‘* + ^ ,я’'пЬ^ ^15°^

Отсюда

 

 

 

 

 

€>] = —(и,-у + Муу) + ^гк1()тп^1т,пк!

(1.51)

 

 

 

 

 

Щ — -(ц у —«у,!1) + С,ЫСутп<2/т,п*.

(1.52)

 

 

 

Следовательно, согласно (1.51) всякий симметричный тензор

второго ранга е представим в виде

 

 

 

е = Бе1« + 1пк$.

(153)

Упражнение

1.10. Доказать тождества

 

 

 

1пкБеГ = 0,

(1-54)

 

 

БЬЛпк = 0.

(1.55)

Упражнение 1.11. Доказать, что если для тензора е выпол­

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

 

 

 

 

Б1уе = О,

(1.56)

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

 

 

 

е = 1пк5.

(1.57)

Заметим, что если справедливы соотношения Коши (1.11), то

 

 

1пк 5 = 0

(1.58)

и из (1.54) и (1.53) следует (1.22).

 

Пусть теперь выполнены условия (1.22). Тогда в

силу (1.53)

имеем

 

 

 

 

 

1пк1пк5 = 0,

(1.59)

что эквивалентно

Д (1пк5) = 0,

(1.60)

а в силу поведения тензора поля на бесконечности отсю да вы­ текает (1.58), а потому из (1.53) следует (1.11), что и требовалось доказать.

§ 2. О С Н О ВН Ы Е ЗА К О Н Ы М С С

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

 

 

 

(2.1)

а

= /

рР А У + 1 5 (п) сЕ ,

(2.2)

I

х щ а у

= /

р[х х Р ]< 1 У + / * X

(2.3)

 

V

Г

 

Здесь р — плотность вещ ества (г/см3), Р — вектор массовых сил (см/с2), — вектор поверхностных сил, действующих перпен­ дикулярно к площадке, характеризующейся единичным вектором нормали п. Известным образом можно ввести тензор напряжения о, причем

§(п) = ^ * п , = <ту пу. (2.4)

Э тот тензор может быть разложен на дивиаторную и шаровую составляющие:

 

О- = 5 -Ь 0-7,

<гу = 5у + <г$у,

(2.5)

где

первый инвариант

тензора напряжения,

называемый

средним

напряжением

 

 

 

0 = 3<г=(?) = <гу.

(2.6)

* Мы здесь не будем рассматривать полярные среды [69], т.е. среды, в которых учитываются массовые и поверхностные распределения моментов.

В качестве второго инварианта тензора напряжения чаше всего будет приниматься интенсивность тензора напряжения

= = у/ • (2.7)

Следствием закона сохранения масс является уравнение нераз­ рывности, которое при лагранжевом способе описания имеет вид

р<№ = сопз!.

(2.8)

Таким образом, в М Л Т Т изменение плотности может быть найде­ но после определения изменения кинематических характеристик среды.

Упражнение 2.1. Пользуясь соотношениями (2.4), (2.8) и (1.4), вывести из закона об изменении количества движения (2.2) урав­

нения движения сплошной среды:

 

рй" = рР + БЬго’, ри- = рР{ + гг; ^-.

(2.9)

Упражнение 2.2. Пользуясь соотношениями (2.4), (2.8) и (2.9), а также законом об изменении момента количества движения (2.3), доказать симметричность тензора напряжения:

<х = ёг,

(2-10)

Умножим левую и правую часть уравнения (2.9) скалярно на вектор V и проинтегрируем по объему, занимаемому средой:

I рьь<1У = I р Р ь < 1 У + у р п г ? ).5</К

(2.11)

Последний интеграл левой части (2.11) преобразуем следующим образом:

^(IV = ^ <ГцЬ(п; <1%- ^ <т,7и,); ЛУ. (2.12)

Подынтегральное выражение последнего члена в (2.12) на основа­ нии (2.10) и (1.32) можно представить в виде

= \(<П]Ъ,з + «•*«*,*) =

(*'«•,1

(213)

Обозначим ЛЕ изменение кинетической энергии, 6 А ^ изменение работы-внешвих .хил, состоящей из изменения работы внешних

массовых сил

и внешних поверхностных сил 6А^\ а

изменение работы внутренних сил:

 

 

АЕ

рь

 

 

(2.14)

6А\е) = А1

1 рР - у АУ =

1 рЕ-АйАУ,

(2.15)

 

 

У

 

У

 

 

 

V

 

V

 

8А^^ = А1

[

АЕ =

[ 5 ( п ) АиАЕ,

(2.16)

УУ

ЕЕ

6А& = -А1 [ <гУ5у АУ = - 1 <г^Ае^ АУ.

(2.17)

УУ

V V

Тогда из (2.11) получаем теорему живых сил

АЕ = <е>+ 6А&, 6А<е>= 6А[е) + 6А^е).

(2.18)

Заметим, что АЕ является полным дифференциалом кинетической энергии Е :

Е = Ц р ь 2 АУ,

(2.19)

V

чего нельзя сказать об изменении работы внешних и внутренних сил. Этим и объясняется различнее обозначениях, т.е в символах

А и 6. В самом деле, например, выражение будет полным дифференциалом работы внешних массовых сил

А[е) = у рР Ы У

(2.20)

V

только при условии потенциальности массовых сил, т.е. сущест­ вования такой скалярной функции хг(*> й), что

Е =

дх р

(2

.2 1 )

дщ

 

 

То же замечание относится и к работе внешних поверхностных сил

А(2е) = 1 5<п>-пАЕ,

=

(2.22)

Е

Заметим, что потенциалы и Х з существуют, в частности, если

Ги 5 (п} не зависят от перемещений 2.

Если рассматриваются неизотермические процессы, то требу­

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

<Е. (2.23)

Здесь ^ — массовый приток тепла, ^п — проекция вектора внеш­ него потока тепла д на нормаль к поверхности Е:

9(п) = 9 -п = 9,-п,-.

(2.24)

Внутренняя энергия, как и любая функция состояния, зависит от внешнего термодинамического параметра состояния, которым является температура Т, и некоторых внутренних параметров состояния, которые характеризуют рассматриваемую среду.

Второй закон термодинамики постулирует существование не­ которой функции состояния Я , называемой энтропией, и записы­ вается в интегральном виде следующим образом:

=

 

+

(2.25)

V

V

Е

V

 

где IV* — так называемая функция рассеивания, причем IV*

0.

Если функция рассеивания равна нулю, то среда называется обратимой; если она строго больше нуля, то среда необратимая. Последнее слагаемое в правой части уравнения (2.25) называют

производством энтропии Я *,

причем

 

 

Я *

У2

0.

(2.26)

 

 

 

Заметим, что все пять законов (2.1), (2.2), (2.3), (2.23) и (2.25) имеют одинаковую структуру «законов сохранения». Законом изменения величины а называется выражение вида

^ рАОУ + I В<п>оК + / С (IV,

(2.27)

где все величины а, А, в ( п> и С имеют одинаковую тензорную структуру; А называется источником величины а; В*”) — потоком величины а:

В<"> = В п,

(2.28)

причем В — тензор на единицу большей валентности, чем а; С на­ зывается производством величины а, если задано дополнительное требование С ^ 0. Если А = В (п) = С = 0, то (2.27) называется законом сохранения величины ра. Пользуясь выражением (2.28) и применяя теорему О строградского-Гаусса к выражению (2.27), получим закон изменения величины а в дифференциальной форме:

ра = рА +

В + С.

(2.29)

Воспользуемся теперь теоремой живых сил (2.18) и обозначениями (2.14)—(2.17). Тогда первый закон термодинамики (2.23) можно переписать в виде

рС/ =

I

+ рд] Ы - { ?<п><Е.

(2.30)

V

V

Е

 

Упражнение 2 .3 . Проанализировав уравнения (2.1)—(2.3), (2.30) и (2.23), установить, об изменении какой величины в каждом из них идет речь, что является источником, что — потоком и что — производством этой величины.

Упражнение 2 .4 . Пользуясь выражениями (2.29) и (2.25), пре­ дставить первый закон термодинамики (2.30) в дифференциальной форме

рН = <г$ еи +

Р

Я

- (2.31)

Упражнение 2.5 . Пользуясь выражениями (2.29) и (2.25), пре­ дставить второй закон термодинамики (2.25) в дифференциальной форме

рТ Н = р д ~ 9М + 1У'.

(2.32)

Упражнение 2 .6 . Вводя функцию состояния Ф, называемую, свободной энергией, по формуле

Ф = V - Т Н

(2.33)

получить из выражений (2.31) и (2.32) закон сохранения в форме

рФ + РН Т + № * =

.

(2.34)

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