книги / Численные методы в теории упругости и пластичности
..pdfинтегралу в (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, причем каждый из них пересекает только одну поверхность 5а [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), то выполняются условия обращения в нуль всех компонент сим метричного тензора второго ранга Я :
Яу = Аеу + |
- с и м ~ *}>=№+ &](ек1,к1 —Д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ц = 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А[е) + 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) можно переписать в виде
рС/ 6У = |
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) |