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

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

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

Выш е было показано, что решение задачи «А» является также обобщенным ее решением.

Покажем, что если обобщенное решение является достаточно гладким, то оно является решением задачи «А».

В самом деле, решение задачи «А * должно удовлетворять условиям (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) для всякой достаточно гладкой тензорфункции т € То. Другими словами, обобщенным решением задачи

«В » называется тензор-функция <х € Г , удовлетворяющая для всякой гладкой тензор-функции т€ Го интегральному тождеству

I а№ )тц ёУ = А ъ (3(°>).

(4.45)

V

 

Выше было показано, что решение задачи «В » является также обобщенным ее решением. Бели обобщенное решение задачи явля­ ется достаточно гладким, то оно является решением задачи «В ».

В самом деле, решение задачи « В » в односвязной области должно удовлетворять условиям (4.3), (4.40), (4.41). По определе­ нию обобщенного решения выполняются соотношения (4.3), (4.40)

ивторое из условий (4.4). Вводя систему гладких функций я$(х),

хе V (обобщенные множители Лагранжа), можно записать:

^

ёЕ - ^

ёУ - ^ г ч п к ёУ = 0.

(4.46)

2 ,

V

V

 

Применяя к (4.46) теорему Остроградского-Гаусса, в силу про­ извольности поля г € Го получим

ец = + *м ). (4.47)

Для того, чтобы существовало непрерывное поле х, необходи­ мо выполнение условий (4.41) причем из (4.47) следует выполнение первого из условий (4.4).

Предположим теперь, что тензор деформаций является потен­

циальным:

 

 

еу = &>(?) =

(4-48)

В этом случае можно ввести кастильяни&н К. по формуле

 

К (а) = -Ф (? ) +

(30), Ф(?) = 1 * ё У .

(4.49)

 

V

 

Рассмотрим теперь некоторый линейный тензор-оператор от

напряжений

 

ху = Пу(е),

(4.50)

такой, что в функциональном пространстве а € Т0 величина

 

= У *у(?)ту ёУ

(4.51)

V

удовлетворяет всем аксиомам скалярного произведения, так что рассматриваемое пространство 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)

существует обобщенное решение * задачи «В», единственное в этой окрестности, и при любом значении итерационного парат

метра 0 € (о,

к нему сходится начинающийся с <г(°) процесс

последовательных приближений (4.54)-(4.56), причем

 

 

11<Г(П) ~ 7*11* <

1 ^ Н 7 (1) - 7(0)11*.

(4.60)

где

 

 

 

 

 

 

9 = шах(|1 - 0по\, |1 - /?^о|) < 1.

(4 61)

 

Доказательство следует из рассмотрения тождества

I

*о'(?)г« ЛГ = ! *ц(<г)тч

- 0 У ёц (?)тч &

~ Аъ А ъ 8°)

V

 

V

 

(4.62)

 

 

 

 

и применения к нему процедуры, использованной при доказатель­ стве теоремы 4.1.

Итак, процесс последовательных приближений (4.54)-(4.56) сходится со скоростью геометрической прогрессии со знамена­ телем ?, который принимает наименьшее значение

_

N0

- По

0

=

2

9

N0

при

по + N0

+ 110

 

 

Заметим, что при выполнении условий теоремы 4.2 оператор Ф(<7(")) сходится к Ф(<г*), а следонательно, кастильяниан А^(<г(т )) — к кастильяниану &{<?*)■ В самом деле, полагая в (7.75) гл. 1 Т = <т(т\ получим

Ф(?(т)) - *(? *) ^ А ъ М (т) ~ ?*, 2°) + у ||?(т) - 7*||?. (4.63)

Пользуясь теоремами вложения Соболева [100], получим отсюда для перемещений и0, удовлетворяющих условиям (4.57),

ф(е(т)) - *(? *) ^ ( 5 + т ) ||- ( т ) "

(4'64)

где В — некоторая константа, зависящая только от области Е 1# Поэтому, используя (4.60), имеем

% (т)) - Ф(7*) ^ ( 5 + у ) ?2тЦ?(0) - 7*11? - 0

(4-65)

при т —* оо.

Рассмотрим теперь задачу « Б » М Д ТТ, которая заключается н решении шести уравнений

Щ ( ? ) = Й ц ( ? ) + О н { 5 } + ((ц - 6ц )<?*„ { 5 } = 0

(4 .6 6 )

(определение тензорам-операторам Я , 0 и тензорам 5 и (

дано в

§ 8 гл. 1) и удовлетворении на Б граничным условиям (4.43) и

=

(4.67)

Лля задачи «Б » можно организовать метод последовательных приближений, аналогичный процессу для задачи «В ». В этом случае в уравнениях (4.54) следует заменить тензор-оператор ?){П(?)} на М {П (?)}. (Обозначим уравнение (4.54) после такой замены через (4.54)1.) Соотношения (4.55) должны выполняться только на границе тела:

(<г‘; 7 1> + Х ,)е = 0.

(4.68)

Тогда итерационная схема для задачи «Б » состоит из уравнений (4.54)1, (4.56), (4.68).

§5. БЫ СТРО СХО ДЯШ И Й СЯ М ЕТОД

ПО С Л ЕД О ВА ТЕЛ ЬН Ы Х ПРИБЛИЖ ЕНИЙ

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

д 2<г%)

Ло*

< Ц Н цЬц?*3,

Ь > 0.

(5.1)

^дсыдбтп

 

 

 

 

Предположим далее, что пространство Пх с введенным скалярным произведением

( 5 , 3 ) 1

е« (3 )Л '’

(5.2)

 

V

 

является гильбертовым для функций « € 1*0, определенных в конечной области V. Тогда справедлива следующая теорема.

Теорема 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 ^

\ —СЧа

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