книги / Численные методы в теории упругости и пластичности
..pdfВыш е было показано, что решение задачи «А» является также обобщенным ее решением.
Покажем, что если обобщенное решение является достаточно гладким, то оно является решением задачи «А».
В самом деле, решение задачи «А * должно удовлетворять условиям (4.1)-(4.4). По определению обобщенного решения вы полняются соотношения (4.1), (4.2) и первое из соотношений (4.4). Применяя к тождеству (4.12) теорему Остроградского-Гаусса, получим
^ + |
х >)ь* |
^ (<ТчЪ - 5°)«,- с/Х = |
0. |
(414) |
V |
|
Е2 |
|
|
В силу произвольности поля V 6 Уо отсюда следуют уравнения |
||||
равновесия (4.3) и статические граничные условия |
(4.4). |
|
Предположим теперь, что тензор напряжений является потен
циальным: |
д\У(е) |
|
|
—•^0'(§) — |
(4.15) |
дец |
Если соотношения (4.15) достаточно гладкие, то можно пост роить функциональные производные типа
(4.16)
и справедливо тождество (7.15) гл. 1.
Рассмотрим теперь некоторый линейный тензор-оператор от
деформаций |
|
Р.7 = */ (€ ). |
(417) |
такой, что в функциональном пространстве и € Уо величина
(4.18)
удовлетворяет всем аксиомам скалярного произведения [53], так что рассматриваемое функциональное пространство П является гильбертовым. Пусть, кроме того, оператор (4.17) таков, что для произвольного симметричного тензора Л выполняются нера венства
т 0ру(Л)Лу ^ |
Ьц 4: М>Ру(Л)Лу, 0 < т0 $ Мо. (4.19) |
|
ОСИ |
Заметим, что если
р1] (с) = 2 |
+ &и&]к)Еы —с у , |
(4.20) |
то первое из неравенств (4.19) эквивалентно при т0 = т неравен ству (7.60) гл. 1, причем гильбертовость пространства П в этом случае следует из неравенств Корна [58].
Если теперь существует единственное обобщенное решение задачи «А» для случая, когда оператор Р определяющих соотно шений (4.1) является оператором р,; (4.17) (задача «Ар»), можно организовать метод последовательных приближений
Ш * я+1)) =Р ч Л * п)) - /?(п)Кд(«(п)) +х ,] , |
(4.21) |
Ро(«(п+1))п/|Еа = Р о («(п))«/1е 2 - ^ п)[^ (« (л)К | е 3 -■??],(4.22)
начиная с некоторого нулевого приближения й(°) и полагая п = 0,1,2,....
Теорема 4.1. Пусть существует единственное обобщенное решение задачи «Ар», справедливы условия (4.19), объемные и поверхностные силы принадлежат пространствам Ь ^ У ), причем
х { е ь ,( У ) , Ч> 1 |
з ^ е ь ^ ) , |
? > ^. |
(4.23) |
Пусть, кроме того, для нулевого приближения |
для произволь |
||
ного симметричного тензора А выполняется условие |
|
||
[<М«(0)) -Р у (й (0))]Лу ^ тру(А)Ау. |
(4.24) |
||
Тогда в некоторой окрестности |
|
|
|
||« - |
й(0)||„ ^ Г |
|
(4.25) |
существует обобщенное решение «* задачи «А», единственное в этой окрестности, и при любом значении итерационного пара метра /? € (0, ^ ) к нему сходится начинающийся с й(°) процесс последовательных приближений (4.21), (4.22), причем
Цй(п)-г*11р ^ Л | 5(0>-5*||, |
(4.26) |
где |
|
Ч = тах(|1 - /?т|, |1 - 0М\) < 1. |
(4 27) |
Для доказательства рассмотрим тождество
I Рц («) 4У = ! ру (2)су (3) АV- Р ^ «ту (2)су (3) (IV-А* (г)
(4.28) Слева в (4.28) стоит скалярное произведение (2,2)р. Правая часть представляет собой согласно (4.19) линейный функционал от V . Используя теоремы вложения Соболева [100], легко устаг новить, что для этого необходимо, чтобы выполнялось условие (4.23). Тогда по теореме Рисса этот функционал мажет быть представлен в виде скалярного произведения (и',ь)р где 2' € П. Следовательно, некоторый оператор <3 ставит в соответствие каждой функции 2 Е П функцию и' 6 П. Поэтому вопрос нахож дения обобщенного решения задачи «А» заключается в решении операторного уравнения
2 = ф2. (4.29)
Для двух векторных полей 2х и 22 и их разности й> = щ —й2
имеем из |
(4.28) |
|
(<Э«1 - |
фй2, й)р = (й, й)Р - Р ! [*у 0*0 - ёу(«2)Ь»у № , |
(4.30) |
|
V |
|
где |
|
|
|
10у = Су(2 0 - еу (32). |
(4.31) |
Используя равенство (7.63) гл. 1 и условие (4.19), получим из (4.30)
|
|(<3«1 ~ 0«2,й)|>| < «1И1р, |
(4.32) |
где ^ определяется из (4.27). При этом при 0 < /? < |
т+Ш |
|
|
\ 1 -0 т \ *\ 1-0 М \, |
|
а ПРИ |
< Р < Ш |
|
|
\ 1 -0М \ *\ 1 -0т \ - |
|
Поэтому |
при |
|
|
0 </?< М2_ |
(4.33) |
выполняется условие для д < 1. Введем единичный оператор Е
и оператор Ё согласно (4.28): |
|
$ = Е - р Ё . |
(4.34) |
Тогда для функциональной производной оператора Ё, которую обозначим Ё 1, имеем
||<521-022|К8ир||^-^/(31 +^(22 - и 1))||||и1 - 2 2||^д||^ (4.35)
Заметим, что наименьшего значения |
величина д достигает |
||
при (3 = |
. Заметим также, что на каждом итерационном шаге |
||
можно изменять значения /3 так, чтобы Р ^ |
€ (0, ^ ) . |
|
|
Из неравенства (4.35) следует, что оператор $ осуществляет |
|||
в П сжатые отображения [53]. Далее имеем |
|
||
(<?3 - |
й(°), Й)р = (<Э« - <5#0), ь)р + (<$3(0) - й(0),3)р. |
(4.36) |
|
Но из тождества (4.28) имеем |
|
|
|
(<5«(°> - й(°\ ь)р = р ^ [о-у(3(°>) - Ру(3(0))]еу(3)<*К. |
(4.37) |
||
|
V |
|
|
Применяя к (4.37) условие (4.24) и полагая в (4.36) 3 = |
2 — «(°), |
||
получим |
||<33 - «(0)|^ (д + Ртп)г < |
г, |
(4.38) |
|
т.е. оператор С}, совершая сжатые отображения, не выводит ни одну точку из окрестности (4.25). Поэтому согласно принципу сжатых отображений [49] существует обобщенное решение задачи «А». Его единственность следует из соответствующего примене ния теоремы, доказанной в § 7. гл. 1.
Из формулы (4.35) следует, что процесс последовательных приближений сходится как геометрическая прогрессия со зна менателем д. Более практическое значение имеет следствие формулы (4.35):
|Й<В>-Й*| 1 - д 110й(0 )- й (0)||„ (4.39)
Теорема доказана.
Из сходимости й(") к й* следует также сходимость ^(й^")) к И й*) [76].
Предположим теперь, что операторные соотношения (4.1) од нозначно разрешимы относительно деформаций:
= <?,;(?)• |
(4.40) |
Условиями интегрируемости системы дифференциальных уравнений (4.2) относительно перемещений являются уравнения совместности Сен-Венана, обращающие в нуль симметричный тензор несовместности Г):
Vу = ик1^тпгкп,1т = 0 (т] = 1п1:е = 0). |
(4.41) |
Для односвязной области V условия (4.41) являются необходи мыми и достаточными для однозначной разрешимости уравнений (4.2) относительно перемещений, например, в виде, предложенном Чезаро. Таким образом, если выполняются условия (4.41), то найдется такой вектор и, для которого справедливы соотноше ния Коши.
Если подставим соотношения (4.40) в (4.41) и (1.19) гл. 1, то получим систему шести уравнений относительно компонент тензора напряжений
|
% (? ) = 0 |
(4.42) |
и граничные условия |
|
|
«.'(?) к = |
<туП] |г3 = 5,-0). |
(4.43) |
Таким образом, в (4.3), (4.40), (4.41),.(4.43) или (4.3), (4.42) и (4.43) дается постановка квазистатической (статической) задачи М Д ТТ
в напряжениях |
(задача «В »). |
Очевидно, что постановки задач |
||
«А » и «В » |
эквивалентны между собой. |
т € То и |
||
Помножим скалярно соотношения (4.2) на тензор |
||||
проинтегрируем по объему V. Тогда, используя теорему Остро- |
||||
градского-Гаусса и условия (7.21) гл. 1, получим |
|
|||
I 4 |
У |
= Ае ^ й0), |
АЕ1(и°) = у 7у и “п7<*>Е - |
(4-44) |
V |
|
|
2! |
|
Назовем обобщенным решением задачи « В » тензор а € Т, для которого справедливы соотношения (4.40) и который удовлетво ряет тождествам (4.44) для всякой достаточно гладкой тензорфункции т € То. Другими словами, обобщенным решением задачи
удовлетворяет всем аксиомам скалярного произведения, так что рассматриваемое пространство 5 является гильбертовым. Пусть, кроме того, оператор (4.50) таков, что для произвольного симмет ричного тензора Л выполняются неравенства
ЯоЯ’у (Л)Лу ^ ’дёц |
Ьы Лу ^ |
((*)Лу, |
0 < По ^ Nо. |
(4.52) |
д<ты |
|
|
|
|
Заметим, что если |
|
|
|
|
7Гу(?) = |
+ ЬцЬ}к)<*Ы = |
<гу, |
(4.53) |
то первое из неравенств (4.52) эквивалентно при п = По нера венству (7.71) гл. 1. При таком выборе оператора П обозначим гильбертово пространство 5 через 5<>.
Если теперь существует единственное обобщенное решение за дачи «В » для случая, когда оператор определяющих соотношений (4.40) является оператором П (4.50), то можно организовать метод последовательных приближений
*ЫП(?(т+1))} = % {П(?(т))} - М Лё(<Г(т>)}, |
(4.54) |
= |
(4-55) |
« , { П ( < г ( '" + 1) ) }| 21 = «,- {й (< г ( т ) )}| Б 1 - А » [ М $ ( < Г ( т ) ) } к |
- « ? ] . |
<гу ,+1>п1 = 5*? |
(4.56) |
начиная с некоторого нулевого приближения о^°) и полагая тп =
0,1.......
Теорема 4.2. Пусть существует единственное обобщенное решение задачи «В *», справедливы условия (4.52), заданные пе
ремещения удовлетворяют условиям |
|
«? € 1Р(Е), р > |
(4.57) |
Пусть, кроме того, для нулевого приближения |
выполняется |
условие |
|
[ёу(ог(0)) - *у (? (0))]лУ < «*у(Л)Лу, |
(4.58) |
где Л — произвольный симметричный тензор. |
|
Тогда в некоторой окрестности |
|
1 к - * (0)||х<г |
(4.59) |
Теорема 5.1. («быстросходящийся» метод). Пусть оператор (4.17) имеет вид
Р у (Л )= | ^ Л ы |
(5.3) |
и существует единственное обобщенное решение соответствующей задачи «Ая». Пусть выполнены неравенство (5.1) и неравенство
шхАуЛу ^ [ з ^ Лы] |
^ М*А0 Лч> 0 < т < Мх. |
(5.4) |
Пусть, кроме того, а — такое положительное число, что
У[*<Д2(0))-р{Ди(0))]е*Ди(0))</1' ^ тща^ су(3(0))ео(и(0))ЛУ. (5.5)
V |
V |
Тогда найдется такое число а, |
|
О < а ^ 1, |
(5.6) |
что задача «А» имеет единственное обобщенное решение «* в окрестности
|
||«(0)-«*||, < г 0) |
(5.7) |
если выполняется неравенство |
|
|
|
Я ^ а~а С, |
(5.8) |
где |
|
|
Я = |
' С н о ^ Ц - а ) -1^ ; |
(5 9) |
го — меньший, а п |
— больший корни уравнения |
|
|
Чг1+а - г + а = 0. |
(5.10) |
При /? = 1 к этому решению сходится начинающийся с «(°) процесс последовательных приближений (4.21), (4.22), причем
||й(п) - 3*111 ^ |
|
«о - и*\\?+а)п, |
(5.11) |
или |
|
|
|
и — |
- * ц . |
аС |
(5.12) |
1и(п) |
и ||1 ^ |
\ —СЧа |