книги / Численные методы в теории упругости и пластичности
..pdfоткуда, учитывая произвольность величин 6а, т.е. 6й, 6е и 6<т, получим (7.26)-(7-29).
Бели воспользуемся преобразованием Лежандра, то оператор
2 {а } (7.53) можно переписать в виде |
|
|
2 {« } = I |
+ «>,.) ~ «>{?} ~ |
ЛУ~ |
V |
|
|
~ У |
^ 5?ицИЗ, |
(7.56) |
Е1 |
Ез |
|
где под символом а понимается теперь совокупность величин « и <г. Упражнение 7.8. Исходя из оператора Рейсснера (7.56), сфор мулировать вариационный принцип и вывести из него соотноше
ния (7.26), (7.29) и соотношение
^ = 5 <«« + «/.«)• |
(Т.57) |
Упражнение 7.9. Показать, что лагранжиан является частным случаем оператора (7.53), если считать заранее справедливыми соотношения Коши (7.27) и определяющие уравнения (7.28).
Упражнение 7.10. Воспользовавшись преобразованием
\ ^ |
^ |
^ |
ацт цщ дУ , |
(7 .58) |
V |
V |
Е |
|
|
показать, что кастильяниан является частным случаем операто ра (7.56), если удовлетворяются уравнения равновесия (7.26) и статические граничные условия (7.29//).
Упражнение 7.11. Воспользовавшись принципом Д ’Аламбера, т.е. вводя инерционные силы заменой
рРх —» рЕ, + ри{ , |
(7.59) |
сформулировать вариационные принципы Лагранжа, Кастильяно, Рейсснера для случая, когда учитываются силы инерции. Ш
Предположим теперь, что определяющие соотношения (7.28) таковы, что для любого тензора Н
’д&ц |
М Нху ^ шЛуАу, тп > 0. |
(7.60) |
деы |
|
|
Тогда можно доказать, что стационарная точка лагранжиана (7.31) является точкой минимума. В самом деле, полагая в тождестве (7.15) й2 = 3 (любой кинематической системе), а Зх = 3 (решению задачи (7.26)-(7.29)), имеем, учитывал (7.60),
^ {3} = ^{3} - А(е\$) ^ ^ {3} - Л(е)(и)+
+ у I еу(3 - 3)еу(3 - 3) ё У > ф{Ь) - А ^ ( 3) = Г {3 }. (7.61)
V
При выполнении условия (7.60) справедлива также теорема един ственности обобщенного решения задачи (7.26)-(7.29).
Всамом деле, предположим противное: существуют решения
«1 и 32. Тогда из (7.3) следует, что они удовлетворяют тождеству
|
|
{32} - Ъц{«х})еу(3) ёУ = 0. |
(7.62) |
|
V |
|
|
Лалее, |
|
|
|
|
(^ ,{« 2 } - *<>{«1})еу(3) = |
|
|
= { |
+ |
“ а1)}[еы(Й2) - еы(их)]еу (3)1 <%. |
(7.63) |
о |
|
|
|
Поэтому, полагая в (7.63) еу(3) = еу (32)-е у (Йх), получим из (7.60)
®^ У[0у'{*Ь} —^»>{21}][еу(32) —еу(й\)]ёУ ^
у
> т ^ еу(«2 - «1)еу(32 - Зх) ё У |
(7.64) |
V
Отсюда следует, что
еу(«х) = еу(й2), |
(7.65) |
т.е. поля «х(5) и й2(х) могут различаться между собой только на смещение тела как жесткого целого. Однако в силу условий (7.297) такое смещение недопустимо. Отсюда следует единственность решения задачи (7.26)-(7.29).
Бели же граничные условия являются только статическими, то решение задачи М Д Т Т единственно с точностью до движения тела как жесткого целого:
тц(2) = «? + еук*,«2°> |
(7.66) |
где й° и й00 — два произвольных постоянных вектора. Исклю чить жесткое смещение тела можно, закрепив его какую-нибудь точку х°:
при х = х° й = О, 2 = 0, |
(7.67) |
или, например, потребовав, чтобы
/ й<*К = 0, / V х и (IV = 0. |
(7.68) |
Заметим также, что при статических граничных условиях о разрешимости задачи М Д Т Т можно говорить лишь в том случае, когда «система самоуравновешена», т.е. для всего тела выполня ются уравнения (2.2), (2.3), которые для случая равновесия можно записать в виде
! |
рРА У + ^ 5° |
= 0, |
(7.69) |
V |
Е |
|
|
I р [ х |
х Р ]й У + ! х |
х 5°<*Е = 0. |
(7.70) |
V |
Е |
|
|
Упражнение 7.12. Используя доказанную теорему единст венности, доказать, что точка минимума лагранжиана является единственной. ■
Потребуем теперь, чтобы для определяющих соотношений (7.44) выполнялось неравенство для произвольного симметрич ного тензора й:
Г Ы ц |
^ пйу Йу, п > |
0. |
(7.71) |
|
[д(гк1 Ьк1 |
||||
|
|
|
Тогда можно доказать,что стационарная точка кастильяниана (7.47) является точкой максимума.
§ 8. НОВАЯ ПОСТАНОВКА ЗАДАЧИ М ДТТ И НОВЫЙ ВАРИАЦИОННЫЙ ПРИНЦИП
Рассмотрим уравнения совместности в форме (1.40) гл. 1:
Нц = Де,;- + 0 у - е,-к,к] ~ ^]к,ы + (у(еы ,ы —Д0) = 0. |
(81) |
где ^ — произвольный симметричный тензор-константа второго ранга. Рассмотрим далее уравнения движения
5* = <7у>;- + рР{ —ри- = 0 |
(8.2) |
или в в случае квазистатических (статических) задач М ДТТ уравнения равновесия
5,- = 0 ^ 4 + рРг = 0. |
(8.3) |
Очевидно, вектор 5 получается из вектора 3* пренебрежением инерционными членами ри ". Применим оператор БеГ к вектору 5* (или вектору 3):
57; = 2 |
+ ^м) = 2 ^ к<к} + Ч к& ) + |
~ Реч = |
(8-4) |
Предположим, что ф является симметричным тензором-операто ром от тензора 5* (или тензора 3), таким, что выполняются условия (7.60), т.е. для произвольного симметричного тензора Л выполняется неравенство
рОй-(У) |
Ък1 Лц ^ 0*0ЛуЛу, то > 0. |
(8.5) |
I дЗк1 |
|
|
Тогда в силу доказанной в предыдущем параграфе теоремы един ственности следует, что уравнения
< Ы $ * } = 0 |
(8.6) |
имеют единственное тривиальное решение (8.2), если на границе области V заданы граничные условия
571= = о.
Пусть теперь нам заданы определяющие соотношения (4.60), граничные условия (6.5) и начальные условия (6.6). Будем счи тать область односвязной. Выразим уравнения (8.1) и (8.4) с
помощью определяющих соотношений (4.60) через напряжения. Тогда граничные условия (6.5) можно записать в виде
«<(?)1е , = « ?, |
= 5 ? , |
(8 .8 ) |
а начальные условия (6.6) — в виде
при < = *о <тч - <тцф\ е0-(о) = ец(У ). |
(8.9) |
О бразуй теперь шесть уравнений относительно компонент тен зора напряжений <х:
В ц (а) + 0 О{5 *} + (6,- - 5о )0 * * {5 * } = 0. |
(8.10) |
Новая постановка М ДТТ в напряжениях заключается в интегри ровании шести уравнений (8.10) при удовлетворении шести гра ничным условиям (8.8) и (8.7), а также начальным данным (8.9).
Покажем, что данная постановка эквивалентна классической постановке задачи в напряжениях, данной в § 6.
Всамом деле, если справедливы уравнения совместности (8.1)
иуравнения движения (8.2), то отсюда немедленно следует спра ведливость уравнений (8.10) и (8.7). Пусть теперь дана новая постановка задачи (8.10), (8.8), (8.7), (8.9).
Свернем (8.10) с единичным тензором 6у. Получим
(2 - 0 [Д *(г ) - Ч М (?)] + « - 2 ) 0 * * ( Г ) = 0, |
(8.11) |
где
о
о
Предположим, что 4 Ф 2. Применим оператор Ош к уравнени ям (8.10):
(6ц -&)[Д%)-ем.н(?)Ь + $ « Л Г ) + ( й = 0. (8.12)
Воспользовавшись соотношениями (8.11), получаем из (8.12) урав нения (8.6), а из них в силу (8.7) — уравнения движения (8.2). Тог да из уравнений (8.10) следуют уравнения (8.1), а из последних в силу произвольности тензора ^ — уравнения совместности в виде
Щ] — (»к1^тпЧп,1т —0. |
(8.13) |
Упражнение 8.1. Показать, что квазистатическая (статичес кая) задача М Д ТТ в напряжениях заключается в решении шести уравнений относительно компонент тензора напряжений
Нц (?) + Зо- {■?} + « у - «у )<?** {§ } = 0. |
(8.14) |
при выполнении граничных условий (8.8) и условий |
|
5<|2 = 0. |
(8.15) |
Упражнение 8.2. Показать, что динамическая задача М Д ТТ в деформациях заключается в решении шести уравнений относи тельно компонент тензора деформаций
Нц + Яц { 3 * } + (На - 8ц ) $ кк{ 3 * } = 0, |
(8.16) |
где компоненты тензора 3 * выражены с помощью определяющих соотношений (4.1) через тензор деформаций
Зц = | [Ъкм (§) + ё}км {0\ + Р(^,з + рз,*) ~ /*у |
(8-17) |
при удовлетворении граничным условиям, выраженным в дефор мациях,
“«(сЬг |
= |
°ц {е)пз\ъз = 5?, |
(8.18) |
|
= {в-у,Д§) + Р& ~ /»[«.■(€)]' }1е = 0. |
(8.19) |
|||
и начальным данным: |
|
|
|
|
при I = «о |
ец = |
(П), |
(У). |
(8.20) |
Упражнение 8.3. Показать, что квазистатическая (статичес кая) задача М Д ТТ в деформациях заключается в решении шести уравнений относительно компонент тензора деформаций
Нц + 0 ц {5 } + Цц - 8ц )0кк{3} = 0, |
(8.21) |
где компоненты тензора 3 выражены с помощью определяющих соотношений (4.1) через тензор деформации
ЗУ = |
(?) + ^*,м(§)] + />(Яд + *},<), |
(8-22) |
при удовлетворении граничным условиям (8.18) и
5<Ь = |
{ ^ ( е ) + / > ^ Ь = 0. ■ |
(8.23) |
Образуем теперь от векторов 5* или 5, |
соответствующие |
|
вектор-операторы Я(5*) |
и Я (5*), такие,что уравнение |
|
|
к (5 ) = 0 |
(8.24) |
имеет единственное решение |
|
|
|
5 = 0. |
(8.25) |
Тогда можно сформулировать основные уравнения в дивергентной форме. Для этого заметим, что
А ц(5) = Ни + к ^ (5 ) + Д;-,<(5) - (ц Л к,к(§) = 0. |
(8.26) |
Динамическая задача М ДТТ в напряжениях (8.26) заключается в решении шести уравнений
Ау(3*) = 0 |
(8.27) |
с использованием определяющих соотношений (4.60) при выпол нении граничных условий (8.7), (8.8) и начальных данных (8.9). Соответствующая квазистатическая (статическая) задача М ДТТ в напряжениях заключается в решении шести уравнений (8.26) с использованием определяющих соотношений (4.60) при удовлет ворении граничным условиям (8.8), (8.15). Покажем, что при веденная постановка задачи эквивалентна классической. Пусть, например, дана постановка квазистатической задачи (8.26), (8.8), (8.15). Свернем уравнения (8.26) с единичным тензором 6ц:
(2 - |
0(Д<? - еу,у) - (2 - 0 Д *,*(5) = 0. |
(8.28) |
Применим теперь оператор Бит к уравнениям (8.26): |
|
|
(&Н - 4ч)(Д0 — еы,ы)л + Д^«(5) + (Ьн — )Д*,ьу(3) = 0. |
(8.29) |
|
Из (8.29) и (8.28) |
о |
|
при 2 ф ^ следует, что |
|
|
|
ДЯ*(5) = 0. |
(8.30) |
Учитывал граничные условия (8.15) и свойство оператора (8.24), (8.25), получаем уравнения равновесия (8.3), что и требовалось доказать.
Введем теперь следующие обозначения:
еУ.*^У = |
еУ,**>* = с»> |
^ ^ |
= Ркг |
—Я» |
|
и рассмотрим тензор третьего ранга
Еик = Схз.к + бы^2^ “ е}^ + Ы ^2^» ~ е«^ +
+&Де* ~ ^*) + &(я)й)к + К](д)6ц. - &)&к(я)- |
(8.32) |
Тогда уравнения (8.26) можно записать в виде |
|
Е\]к,к + Уу —0, |
(8.33) |
где |
|
У» = & ,;(* ) + Ь А * ) - М * ,* ( * ) . |
(8.34) |
X = рР. |
(8-35) |
Пусть на границе заданы нагрузки |
|
<гг]п}Ы — |
(8.36) |
Условия (8.15) манено переписать в виде |
|
|
(8.37) |
Пусть, кроме того, заданы взаимнообратные определяющие со отношения (4.1) и (4.60):
<гц = ЯЖе), |
(8-38) |
®у — ?ч(с). |
(8.39) |
Очевидно, что квазистатическая (статическая) задача МДТТ в напряжениях заключается в решении уравнений (8.33) с учетом определяющих соотношений (8.39) при удовлетворении граничным условиям (8.36), (8.37). Та же задача в деформациях заключает ся в решении уравнений (8.33) при удовлетворении граничным
условиям (8.36) и (8.37), если воспользоваться определяющими соотношениями (8.38).
Ладим теперь вариационную постановку рассматриваемой за дачи. Для этого определим такой скалярный оператор 93, за висящий от градиентов напряжений, что выполняются условия
потенциальности тензора (8.32) |
|
Ец к — |
а<в |
(8.40) |
|
|
даЧ],к |
Назовем тензором потоков симметричный тензор второго ранга X, определенный на поверхности Е:
Ху — Б^кпк. |
(8.41) |
|
Запишем функционал |
|
|
I = / ( » - Уц°ц) ы |
- 1 |
<*Е+ |
V |
2 |
|
+ ^ ^(Л?,?,- + В<г„П]Сг{кпк) + |
|
- В5?ацп, «?Е, (8.42) |
где А я В — некоторые размерные постоянные.
Упражнение 8.4. Показать, что формулировка задачи (8.33), (8.36), (8.37), (8.39) вытекает из вариационного принципа
В1(сгу ,&го-) = 0, |
(8.43) |
причем при варьировании функционала (8.42) потоки х не варьи руются (считаются «замороженными»), а подставляются их вы ражения по формуле (8.41).
Упражнение 8.5. Показать, что для линейной изотропной среды при
|
б / - \ — 1 + а - |
^ _ ( 1 + ц)е —и . |
|
|
(8.44) |
тензор Е{]к (8.32) примет вид |
|
|
Еч к = 2^ |
+ 2и (р*6к] + |
~ еРкЬ] + Ь6цЧк + а(6к^ { + 6кщ ) |
|
|
(8.45) |