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

Упруго-пластичность и законы сохранения. Монография

.pdf
Скачиваний:
0
Добавлен:
08.09.2026
Размер:
2 Мб
Скачать
☆
На основе законов сохранения построим решений краевой задачи

r
r
r
r
r



для уравнений плоской теории упругости.
Известно, что
1
)1(4
0211
GyGxGGu )(
x
0212
,)(
GyGxGGv
y
0211
GyGxGG
есть решение системы уравнений (231).
i
G
Здесь
– произвольные гармонические функции, - коэффициент
Пуассона, который можно выразить через параметры Ламе.
Из (239) получаем
21
.
GG
xy
Пусть Положим
21
yx
2
)ln(( xxG .0,ln))(
0
GG
(
vu
yx
, yx
– произвольная точка области с гладкой границей Г.
00
1
),2/1)(
0
vu
xy
222
Gryy
Тогда получаем
xx
0
,2
vu
yx
2
vu
xy
yy
0
.2
2
Из (233) по формуле Грина имеем
,)(
x
(239)
(240)
(241)

1
yx
))(
yx
Пусть теперь контур Г радиуса
:
2
0
Полагая в (242) вместо Г контур Г получим
1
Г
1
2
)2((
Запишем уравнение окружности в параметрическом виде:
yyxx
00
(243) запишется так

окружность
2
0
1
1
Г
1
dyuv
состоит из двух контуров – Г и окружности
1
22
)()(
yyxx
0
xx
0
2
1

.
1
yx
cos2)2((
dsin)sin2
=
BdxAdydxdyBA )(
dxuv
20,sin,cos

vu
yy
0
dy
)2
2

Г окружность
, тогда второй из интегралов в формуле
1
1
yx
1
Г
))(2((
yx
и с учетом соотношения (241),
1
)(,(2
00
1
1
.0))(
dcos)sin2
yx
)2((
1
))(2((
vu
yx
+
.0))(
dxuv
xx
0
2
.
))2(
+
2
cos2)2((
(242)
(243)
(244)
51
В формуле (244) была использована теорема о среднем.
r
r
dxr
r
r
r
r
r
r
r
r
r
Окончательно из (244) получаем
1
)2((2
1
+
)(2)2((
1
Г
xx
0
2
1
xx
0
)(2)2((
dy
2
r
yy
0
))(2
2
1
yy
0
)(2
2
1
).,(
yx
00
Теперь в уравнениях (240) положим
122
.0,ln))(
Gryy
+
0
2
В этом случае аналогичные рассуждения и вычисления дают форму­лу для
1
).,(
yx
00
Имеем
Формулы (245), (246) позволяют вычислить
),(00yx
из области, если мы знаем значение этих функций на границе Г.
А так как
1
)2((2
1
+
)(2)2((
1
Г
yy
0
2
dx
yy
0
2
r
1
2
1
)(2)2((
dy
xx
0
))(2
11
,
является решением уравнений (231), то мы получаем метод
xx
0
)(2
2
1
,
).,(
yx
00
11
в любой точке
для построения решений системы (231) по их граничным условиям.
А именно: пусть поставлена следующая краевая задача
Имеется область
, с границей Г. На Г заданы компоненты вектора
смещения
),,(|
Г
yxuu
0
Г
).,(|
yxvv
0
(247)
(245)
2
)ln(( xxG
0
(246)
Тогда, на основании формул (245),(246) решение задачи (247) для системы уравнений (231) можно записать в виде
1
)2((2
u
0
Г
v ))(2
dx
0
xx
0
)(2)2((
v )(2
dy
2
r
yy
0
2
0
yy
0
2
).,(00yxu
v
)(2)2((
0
. (248)
xx
0
2
1
u
0
)2((2
Г
u ))(2
dx
0
yy
0
)(2)2((
v )(2
dy
2
r
xx
0
2
0
xx
0
2
v
)(2)2((
0
yy
0
2
).,(
yxv
. (249)
000
Заключение. Найдены законы сохранения, зависящие от производ-
ных первого порядка и линейные по ним. Показано, что два коэффициента
52
при производных являются решениями уравнений упругости в плоском







случае. Это позволило свести решение краевых задач в перемещениях к вычислению контурных интегралов по границе области. Предложенная методика без труда может быть обобщена на многосвязные области, а также на трехмерный случай.
3.2. ГРУППОВОЕ РАССЛОЕНИЕ УРАВНЕНИЙ
ДВУМЕРНОЙ ТЕОРИИ УПРУГОСТИ И ЗАКОНЫ СОХРАНЕНИЯ
На основе группового расслоения: разложения системы на авто­морфную и разрешающую системы, которые являются системами диффе­ренциальных уравнений первого порядка. Построены бесконечные серии законы сохранения для разрешающей системы уравнений и автоморфной системы. Эти законы позволили решить краевые задачи для асимметрич­ной теории упругости в двумерном случае. Решения эти построены в виде квадратур, которые вычисляются по контуру исследуемой области.
Предлагается следующая связь тензоров напряжений и тензора де­формаций
,)2(
221111
2
1212
,)2(
(250)
112222
– компоненты тензора де-
ij
где
– компоненты тензора напряжений;
ij
формаций;
0,0
– постоянные Ламе, то (250) есть классический закон
Гука для изотропного, однородного случая.
Подставляя (250) в уравнения равновесия, то в случае отсутствия массовых сил получаем
где
)2(,0)(
)(
vu
xyxy
vu,
– компоненты вектора деформаций, индексы внизу, если не указа-
vu
xyxx
vu
xyyy
uv
xyyy
(251)
,0)2(
но иное, означают производные по соответствующим переменным.
Групповые свойства уравнений асимметричных уравнений упруго­сти в динамическом случае изучены в [7], там, в частности, показано, что система (251) допускает бесконечную группу точечных преобразований, порождаемую операторами
где
21
,
hhX
(252)
vu
21
,hh
– произвольное решение уравнений Коши- Римана
hh
,021
yx
53
.021xyhh
(253)
Сделаем групповое расслоение системы уравнений (251) по методу




[1] на подалгебре, порождаемой (252). Для этого продолжим операторы (252) на первые производные. Имеем
21
hhXX
vyux
1
yx
21
,
hh
(254)
vxuy
xy
Дифференциальные инварианты для (251), с учетом (254), имеют вид
,
,1xI
,2yI
3 yx
vuI
4 xy
.
vuI
(255)
Тогда автоморфная система уравнений имеет вид
yx
),,( yxvu
xy
).,( yxvu
(256)
Напомним некоторые свойства автоморфных систем. Любое реше­ние автоморфной системы может быть получено из одного решения этой системы с помощью преобразований, порождаемых оператором (252).
Подставляя (257) в (251) получаем разрешающую систему
F
1

yx
2
F
,0)2(
,0)2(
(257)
xy
Повторяя почти дословно рассуждения из [7], можно утверждать, что система (257) равносильна системе уравнений (251).
Поэтому построив решение системы (257) мы получим решение сис­темы (251).
Пусть для системы (257) поставлена следующая краевая задача:
),,(|
L
yx
0
L
),,(|
yx
0
(258)
где L – некоторая гладкая замкнутая кривая.
Для решения этой задачи построим законы сохранения для системы уравнений (257).
В силу линейности системы (257) она будет иметь бесконечное чис­ло законов сохранения. В работе будут найдены только те законы сохране­ния, которые позволят решить краевую задачу (258).
Определение. Законом сохранения для системы уравнений (2.8) назо­вем выражение вида
,
где
– некоторые функции, которые не равны тождественно нулю од-
новременно;
yxA
x
BA,
– называются компонентами сохраняющегося тока.
),,,(
y
FFyxB
,0),,,(
(259)
21
Более подробная информация по построению законов сохранения для произвольных систем дифференциальных уравнений может быть най­дена в главе 2 или в [8–10].
54
Предположим, что компоненты сохраняющегося тока имеют вид

где
21
aa
21
,,
,bb
– некоторые функции от x, y.
Тогда из (259) получаем
1
Отсюда получаем
Из (259) следует
Пусть
21
,
aaA
a
00
Syx ),(
, тогда из (259) следует
2
),2(

a
11
,0
ba
yx
,0
yx
yx
1
,
0

LS

L
21
,
bbB
(260)
b
.022yxba
.0
(262)
xy
BdxAdydxdyBA
BdxAdyBdxAdy
.)(
(263)
.
(264)
2
),2(
,

a
(261)
где
2
0
Syx ),(
00
, если функции A, B имеют особенности в этой точке. Вычислим
)()(:
0
22
yyxx
– окружность радиуса вокруг точки
интеграл в правой части (264) для разных решений уравнений Коши­Римана. Пусть
yy
xx
BdxAdy
0
0
2
yyxx
0
,
2
)()(
0
2
0
,
(265)
2
)()(
yyxx
0

dxdy
.))2(())2((
(266)
Подставляем (265) в (264) и сделаем замену переменных
xx
0
,cos
yy
0
BdxAdy
,sin
получаем
2
0

d)]cossin2))2(([
(267)
)],(),()2[(2
yxyx
.
0000
В формуле (267) устремили
0
и использовали теорему о среднем.
Теперь сделаем аналогичные вычисления, положив
xx
yy
0
0
2
,
2
)()(
yyxx
0
55
0
2
0
,
2
yyxx
)()(
0
В результате получим



BdxAdy
).,(2
yx
(268)
00
Формулы (267) и (268) позволяют, с учетом граничных условий (259) определить значения функций θ и ω в произвольной точке
Syx ),(
.
00
Теперь после восстановления решений разрешающей системы, най­дем решения автоморфной системы, т. е. решения исходной системы урав­нений (251).
Имеем
,0),(
3
yx
yxvuF
4
xy
.0),(
yxvuF
(269)
Здесь в правой части стоят известные функции. Найдем законы со­хранения уравнений (269) в следующем виде
где
,,43aa
Имеем
Расщепляя систему уравнений (271) получаем
Отсюда получаем
Пусть для системы (269) поставлена следующая краевая задача:
143
,
caaA
2143
,,, ccbb
– некоторые функции от x, y.
x
3
),,,( vuyxA
4
,
a
,033yxba
0
L
y
3
,
a
yx
),,(|
yxuu
b
,0
L
243
,
cbbB
4
b
,
2144
ccba
yxyx
.0
(273)
xy
),,(|
yxvv
0
(270)
,0),,,(
FFvuyxB
(271)
43
,
(272)
.,0

(274)
Рассмотрим закон сохранения в виде
L
BdxAdy
.
BdxAdy
(275)
Пусть решение уравнений (272) имеет вид
xx
0
0
2
,
2
yyxx
)()(
0
56
yy
0
2
0
,
(276)
2
yyxx
)()(
0
Подставляем (276) в правую часть (275), получаем
Пусть решение уравнений (273) имеет вид
Подставляем (278) в правую часть (275), получаем
деформации
В результате получаем формулы для вычисления компонент вектора
BdxAdy
BdxAdy
yy
0
2
0
,
2
yyxx
)()(
0
2
0
21
dxcvudycvu )()(
21
xx
0
. (278)
2
)()(
yyxx
0
21
dxcvudycvu )()(
),(2)cossin()sincos(
yxudxcdyc
. (277)
00
21
),(2)sincos()cossin(
yxvdxcvudycvu
. (279)
00
,),(2
BdxAdyyxu
00
L
,),(2
BdxAdyyxv
(280)
где
1
dxc

2
dxc
,

00
L
.
Таким образом, получены новые бесконечные серии законов сохра­нения, которые позволили решить краевую задачу для разрешающей системы уравнений, а также для автоморфной системы. Все это позволило построить аналитическое решение краевой задачи для уравнений асиммет­ричной теории упругости.
3.3. ЗАКОНЫ СОХРАНЕНИЯ
ДЛЯ ПРОСТРАНСТВЕННЫХ УРАВНЕНИЙ ЭЛАСТОСТАТИКИ
В параграфе построены новые законы сохранения для уравнений Ламе в трехмерном стационарном случае. При этом предполагается, что сохраняющийся ток закона сохранения линейно зависит от первых производных компонент вектора перемещений по пространственным ко­ординатам.
Законы сохранений для дифференциальных уравнений были введены Э. Нетер более ста лет назад. Она показала, что для систем дифференци­альных уравнений, которые выводятся из вариационного принципа, точеч­ные симметрии позволяют генерировать законы сохранения. Долгое время это служило основным методом построения законов сохранения. Позднее в работах А. М. Виноградова [9; 11] было показано, что законы сохранения
57
имеют двойственную природу по отношениям к симметриям и генериру-


ются симметриями только для уравнений, операторы линеаризации кото­рых являются самосопряженными. Это позволило строить законы сохра­нения и для уравнений, которые не выводятся из вариационного принципа. Новая методика позволила существенно расширить возможности для по­строения законов сохранения, но законы сохранения практически не ис­пользовались для исследования и решения дифференциальных уравнений. А. М. Виноградовым и С. И. Сенашовым в 1987 г на конференции по Групповому анализу дифференциальных уравнений в г. Красноярске была выдвинута гипотеза, что законы сохранения можно использовать для ре­шения краевых задач для систем дифференциальных уравнений. Эта гипо­теза основывалась на том, что законы сохранения по своей природе глобальны и лучше подходят для решения краевых задач, чем точечные симметрии. Первый такой результат был получен в работе [12] для уравне­ний двумерной теории пластичности. Позднее было показано, что законы сохранения можно использовать и для решения упругопластических задач. Ссылки и результаты по решению упруго пластических задач можно найти в работе [9]. Для уравнений линейной теории упругости законы сохране­ния, на основании теорем Нетер и их модификации были впервые построе­ны в [3; 4].
Постановка задачи. Рассмотрим уравнения линейной теории упру-
гости в пространственном случае, в предположении, что массовые силы отсутствуют
11321
,0)())(()2(
Здесь
1
2
3
,
– упругие постоянные Ламе,
i
uuuuuF
zzyyxzxyxx
22312
uuuuuF
zzxxyzxyyy
33213
uuuuuF
yyxxyzxzzz
)3,2,1( iu
компоненты вектора
(281)
,0)())(()2(
.0)())(()2(
деформации, индекс внизу означает производную по соответствующим пе­ременным.
Определение. Законом сохранения для системы уравнений (281) на-
зовем выражение вида
где
CBA ,,
– называются компонентами сохраняющегося тока;
zyx
FFFCBA
(282)
332211
i
)3,2,1( i
–
некоторые линейные дифференциальные операторы, одновременно не равные тождественно нулю.
Замечание. Более общие и точные формулировки определения зако-
нов сохранения и методов их вычисления, подходящие и для других сис­тем дифференциальных уравнений можно найти в главе 2 и работах [8–11] и цитируемой там литературе.
58
В этой работе мы будем предполагать, что компоненты сохраняюще-
zyx
zyx
гося тока имеют следующий вид
где от
от
)2(
,,
.
В этом случае
,,
.
111111
uuuwA
zyx
111111
ucubuaB
zyx
111111
uDuBuAC
zyx
iiiiiiiiii
DBAcbaw ,,,,,,,,,
)3,2,1( i
i
тоже будут некоторыми функциями
222222
uuuw
,
zyx
222222
ucubua ,
zyx
222222
uDuBuA
zyx
)3,2,1(i
– функции зависящие только
uuuw
ucubua
1333333
2333333
zyx
(283)
3333333
,
zyx
zyx
uDuBuA
Законы сохранения для систем (281). Будем искать сохраняющийся ток для системы уравнений (281) в форме (282). Для этого подставим (281) и (283) в (282) получим, после некоторых преобразований и исключе­ния
, следующие выражения
i
31322331
,)(,,)(
wAAcwA
(284)
32311232
,)(,,)(
cwBcBcwB
(285)
221133
,,,)2(
wDwDwD
111
cw
zyx
ww
,0
222
cw
ww
xyz
www
111
cw
yzx
2212
,0)(
zyyx
2232
cww
zxzy
,0
(293)
3313
,0)(
zyzx
(286)
33212121
,,)(,)(
awawa
(287)
332211
,,)2(,
wbwbwb
(288)
13121
,0)()()2(
zzyyx
(289)
,0
(290)
(291)
,0)()2(
(292)
(294)
zxzy
(295)
333
.0)2(
cw
yzx
(296)
3323
,0)(
cww
59
Из уравнений (284)–(296) следует, что система уравнений (281) до-
–
пускает бесконечную серию законов сохранения вида (282). Выпишем в явном виде некоторые.
Дифференцируем первое уравнение (290) по y, а второе уравнение (290) по z. Получаем
111
cw
zyyyxy
Продифференцируем уравнение (289) по x и сложим его с уравне­ниями (287), в результате получим
Аналогично из уравнений (291)–(293) получаем
А из (294)–(296) имеем
Из уравнений (298)–(300) следует, что решение системы уравнений (281).
В этом случае сохраняющийся ток имеет вид
)2(
111111
uuuwA
zyx
,0
uuuw
111
.0
cw
yzzzxz
222222
zyx
(297)
11321
,0)())(()2(
wwwww
zzyyxzxyxx
22312
wwwww
zzxxyzxyyy
33213
wwwww
yyxxyzxzzz
321
,, www
(298)
.0)())(()2(
(299)
.0)())(()2(
(300)
– есть произвольное
uuuw
zyx
1333333
,
(301)
где ременные
Уравнения (302) можно кратко записать в виде
))((
))2((
321
,, www
– есть произвольное решение системы уравнений (281), а пе-
iii
)3,2,1(,, ic
удовлетворяют следующим уравнениям.
1111112
ucuwuwB
zyx
uwucuwC
rot
2333333
,
ucuwu
zyx
1111111
zyx
)2()(
))((
uwucwuw
zyx
uwuсwu
333332331
,)2())(())((
wcwcwww
(303)
,
ii
2222322
zyx
11111132111
zyxyzxzyxzy
2222222
ucuwuw
zyx
+
(302)
,,,)()()2(
22232222122
,,)()2(,)(
wcwwcww
zyxzyzxxyzy
33332333133
.)2(,)(,)(
wcwwcww
zyxyzzxxzzy
60