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

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

.pdf
Скачиваний:
0
Добавлен:
08.09.2026
Размер:
2 Мб
Скачать
☆
где
x
x

x

x
x
x

x
x
x
f
L
L
kk k k

El
,
– положительные корни трансцендентных уравнений
mn
[]
nm
12 1 2
1
2222
()()
kk
12
n
 
nn
()()
kk k k
34 3 4
[],
l
2
2222
()( )
kk
24
2

()()
2

m
 
nm
2
12 1 2 1
()(),kk k k ctg l
2
 
34 3 4 2
()().kk k k ctgl
Таким образом, приведен класс новых решений уравнений классиче­ских уравнений упругости. Эти решения можно использовать не только для описания распространения сдвиговых волн в упругих телах, но и для описания напряженно – деформированного состояния упругого параллеле­пипеда скручиваемого вокруг трех ортогональных осей, как это сделано в [16] в случае пластической деформации.
3.9. ТОЧЕЧНЫЕ СИММЕТРИИ УРАВНЕНИЙ ИЗОТРОПНОЙ УПРУГОСТИ
Рассмотрим уравнения упругости в изотропном случае, без массовых сил см. гл. 1, формула (281).
W W grad divW
tt
()
(350)
Ищем точечные симметрии этой системы уравнений согласно мето­дике, описанной в гл. 2.
Получаем, что точечные симметрии порождаются следующими опе­раторами, обозначения взяты из [6], имеем

 
01 2 3
www

123
,Pw w w
 
Pf f f
12 3
f

www
123

,
(351)


Zx x
12 3

32


Zx x
31 2
где

21
,,
ff
123
– произвольное решение системы уравнений (350).
1

x
x
23


12

,X
1
2
,X
2

,ww
Zx x

ww
32
23 1

,ww
ww
21
0
,X
3
(352)
3
x
13
,X
Rx
t
1
1


31
ww

13
x
2
2
3
3
,ww
,x
Алгебра Ли, порождаемая операторами (351), (352) бесконечномерная.
Операторы (351) порождают нормальный делитель этой алгебры, фактор алгебра по нормальному делителю порождается восемью операто­рами (9.3) и порождают алгебру Ли
подалгебр
…,
,
,
1
2
.L
7
L
. Укажем оптимальную систему
8
81
Оптимальная система одномерных подалгебр имеет вид
R
R
Z
Z
Z
Z
Z
Z
Z
,
Z
3
03
,XZ

033 0
,, ,XXZ RX
03
Оптимальная система двумерных подалгебр имеет вид
,,XRZ
03
03 3
,,XX Z
03 3
,,XXRZ 
,,XX
03
,,XXX
013
,,XR3,,
0
Z
,.XXXZ

0303
Оптимальная система трехмерных подалгебр имеет вид
,, ,XX RZ
03 3
01 3
,, ,XX X X Z 
12 0 3 3
,, ,XRZ
03
,, ,XXRX
,,,XXR
03
,, ,XXXX
0123
,, ,XX X Z
120 3
,, ,XXX
013
03 3
,, ,XX RZ
12 3
,, .
ZZ
123
Оптимальная система четырехмерных подалгебр имеет вид
,,, ,XXRZ
03 3
,,, ,XXXX
0123
,,, ,XXRZ
12 3
,,, ,XXRX
03 3
0312 3
,,, ,XXXXRZ 
012 3
,,, ,XXXRX
,,, ,XXX RZ
013 3
,,,.
ZZR
123
,,, ,XXX RZ
012 3
031203
,,, ,XXXXXZ
,,, .
ZZX
123 0
Оптимальная система пятимерных подалгебр имеет вид
,,,,,XXXXR
0123
,,,, ,XXXRZ
012 3
,,,, ,XXXX RZ
0123 3
03
,,XXR
,, ,XXRZ
,,,, ,XXXXRZ
0312 3
,,,,.
ZZXR
123 0
Оптимальная система шестимерных подалгебр имеет вид
,,,,, ,XXXXRZ
0123 3
,,,, , .
ZZXX X
123 1 2 3
Оптимальная система семимерных подалгебр имеет вид
,,,,, , ,
ZZX XXX
123 0 1 2 3
,,,, , ,,
ZZXX XR
123 1 2 3
Построенная система оптимальных подалгебр позволяет строить ин­вариантные и частично инвариантные решения любого ранга, которые до­пускаются необходимыми условиями соответствующих теорем.
Замечание. Точечные симметрии для систем линейных уравнений можно вычислять и по-другому. Опишем это на примере системы уравне­ний (350)
Систему уравнений (350) запишем так

      

Lw w

  
(2) ( ) ( )
tt xy xz
() (2) () 0. () () (2)
     
xy tt yy yz
xz yz tt zz
w

1

2
 
w
3

82
где
L
A
i
нений
Ищем допускаемый оператор в виде


1
x

X

2
y
 

3
z

1, 2, 3
i
есть функции только от
Тогда необходимо решить на
[] ,
XLXLX L
,,,txyz
i
i
 (353)
,
1, 2, 3
.
следующую систему урав-
где
искомая функция от
,,,txyz
.
Ясно, что из системы уравнений (353) мы не сможем найти операто­ры вида
,fPP
, которые допускаются любой системой линейных диффе-
0
ренциальных уравнений, но зато остальные операторы вычисляются гораз­до проще. Поэтому если пользоваться способом, описанным здесь, кото­рый подходит для большинства линейных систем дифференциальных уравнений, то необходимо помнить об этих двух операторах, образующих нормальный делитель в полученной алгебре Ли.
3.10. ОБЩИЕ РЕШЕНИЯ И ВЫСШИЕ СИММЕТРИИ УРАВНЕНИЙ
ЛИНЕЙНОЙ ТЕОРИИ УПРУГОСТИ
В этом параграфе предлагается метод, позволяющий свести уравне­ния теории упругости к более простой, диагональной системе. При этом обратное преобразование не требует интегрирования, а задается транспо­нированной матрицей. Приводятся формулы производства новых решений из известных решений. Подробно рассмотрен случай изотропного мате­риала. Для него указан способ построения высших симметрий и законов сохранения. Этот параграф написан на основе работ [17; 18].
Запишем уравнения упругости при произвольной анизотропии и от­сутствии массовых сил в декартовой системе координат
j
где
123
(, , )wwww
Lw L L A
ij ij ji i kl j kl tt
– вектор смещения,
0, ,
стоянный тензор модулей упругости,
– символ Кронекера.
ij
L
Пусть матрица невырожденная матрица
подобна некоторой матрице D, т. е. существует
T
такая, что
LT TD
()
()/2,
()
i kl j iklj ilkj iklj
AA A
(354)
– постоянная плотность материала,
(355)
83
– по-
что
L
L
L
D
L
L
Так как
'DD
, то
'LL
(штрих означает транспонирование) и предполагаем,
Пусть
0,Du
где
а w – решение системы (354). Тогда с учетом (355), (356) име-
ем
Таким образом, по формулам (359) решения уравнений (358), (359) переходят одно в другое.
Так как
'LL
, то D – диагональная матрица. Тогда из (358) имеем
ttD
1111
ij j i
Т. е. получили задачу на собственные операторы
,,
торы
tt t
123
jj j
для матрицы операторов
''.TL DT
''0.Du DT w T Lw
,
ttD
2222
ij j i
(356)
,wTu ,uTw


(357)
0,Lw LTu TDu
(358)
(359)
,
ttD
3333
ij j i
.
В силу симметричности
ij
,
,,
DD
11 22 33
и век-
.
собственные векторы
ij
Если
t
– числа, то формулы (355) соответствуют ортогональному
jp
,,
tt t
123
jj j
ортогональны.
преобразованию координат, и система (354) становится диагональной для специального ортотропного материала, когда
где
2211 1221 3311 1331 3322 2332
Для изотропного материала
L    
ij ij kk tt ij
,, .AAAAAA  
t
– операторы. В этом случае имеем
jp
()( ),
, – коэффициенты Ламе. В [11] показано, что
() ,
DDD   
123
kk tt kk tt
(360)
(361)
cуть собственные операторы для матрицы (360), а собственные векторы с точностью до произвольных множителей будут
tt ctc c     
12 3
jjjjpspsjmmjjkk
,, ,
(362)
где
– символы Леви–Чивиты;
ijk
c
– произвольный ненулевой числовой
j
или операторный вектор. Таким образом, система (354), (360) для изотроп-
84
ного материала сводится к диагональной системе с операторами (361).
D
B
L
L
B
f
L
s
f
f
Векторы (362) ортогональны
tt cc cc t t            
ip ig p g jj m m sn s n p g sn k k p g ll
Причем
tt cc c c c c            
ji ii jl jps lps p g sr m m j j kk n n l i hh
11 2 2 22 3 3
Ttt cccc     и
k k ll s n m m sn sn ll
Если выполняются (363), (364) и
() ,
22
()
()().
0,Lv
то
– также реше-
'wTTv
(363)
(364)
ние этого уравнения. Аналогичная формула связывает решения уравнения
0, ' .
vuTTv
Соотношение так как из любого решения
является формулой производства решений,
'wTTv
получаем новое решение .w Эту формулу
v
можно применять многократно. Для изотропного материала матрицы
'TT
и
имеют вид (363), (364).
Другой способ производства решений состоит в отыскании матриц
, коммутирующих с оператором ,L т. е. таких матриц, что .LB BL
[]',
За такую матрицу можно взять нение элемента
в матрице .L Для изотропного материала одна из мат-
jk
BL
где
jk
– алгебраическое допол-
jk
'TT
риц
здесь
имеет вид:
 
 
[( 2 ) ) [( 2 ) ] ,
Bb
kk kk
 
  
[( 2 ) ) [( 2 ) )
kk kk

,,b 
11
b
11 33 3 22 2
33 3 11 11 1
22 2 11 1 11
– произвольные операторы с постоянными коэффициентами.
[( 2 ) ] [( 2 ) ]
    
kk kk
b
Результаты первого пункта позволяют вычислить высшие симметрии для уравнений упругости. С помощью симметрий опишем еще один спо­соб производства решений для уравнения (354).
Напомним, что симметрия порядка s для уравнения (354) есть функ­ция
(10.7), где
ff f
ii
такая, что
(, ,..., ), 1,2,3
fxw wi w
f
на всех решения системы уравнений
0
– совокупность всех производ-
s
123
(, , )
ных порядка s.
Формула производства решений
для симметрий означает
'wTTv
наличие оператора рекурсии, который позволяет строить симметрии F сис­темы (10.7) из известных симметрий f, повышая их порядок
(365)
BF TT F
'.
Из этой формулы следует утверждение 1.
85
Утверждение 1. Система уравнений (354) допускает бесконечную
x
x
x


f
алгебру высших симметрий.
Стартуя с точечных симметрий, которые приведены в предыдущем параграфе, получаем требуемое.
Покажем, как можно использовать некоторые из построенных сим­метрий для производства решений. С помощью (365) получаем, что урав­нения (354) допускают симметрии второго порядка вида
111
Aw Bw Cw

xyyzz

fAwBwCw
222
xyyzz
333

Aw Bw Cw
xyyzz
,
 
где A, B, C – произвольные постоянные.
Из определения симметрии следует, что функция
, которая опреде-
w
ляется из решения задачи Коши для системы эволюционных уравнений
dw
,| ,

d
ww
00

есть решение системы уравнений (354), при условии, что вольное решение этой же системы. Здесь
(366)
w
есть произ-
0
– новая независимая перемен-
ная.
В частности, при A = B = C = 1 производство новых решений систе­мы (10.7) сводится к решению задачи Коши для трех уравнений теплопро­водности.
В силу того, что оператор L систем уравнений (354) самосопряжен, следует, что система допускает бесконечную серию законов сохранения. Некоторые из них уже были указаны в предыдущих параграфах.
3.11. ТОЧНЫЕ РЕШЕНИЯ УРАВНЕНИЙ
АСИММЕТРИЧНОЙ ТЕОРИИ УПРУГОСТИ
При решении групповой задачи классификации уравнений, описы­вающих движение чисто механического континуума, появились некоторые новые системы дифференциальных уравнений, которые можно использо­вать для описания реальных физических процессов. Одна из таких новых систем: асимметричная теория упругости. Эта система может быть исполь­зована для материалов имеющих малый модуль Юнга, а также для мате­риалов, которые работают при нагрузках близких к критическим. В этом параграфе на основе группового расслоения: разложения системы на авто­морфную и разрешающую системы, которые являются системами диффе­ренциальных уравнений первого порядка. Построены бесконечные серии законы сохранения для разрешающей системы уравнений и автоморфной системы. Эти законы позволили решить краевые задачи для асимметрич-
86
ной теории упругости в двумерном случае. Решения эти построены в виде





x
x
v
x
x
x
I
I
квадратур, которые вычисляются по контуру исследуемой области.
В работе [6] предлагается следующая связь тензоров напряжений и тензора деформаций
,)(2)(
221201111
– компоненты тензора
ij
(367)
221201122
где
,)(2)(
– компоненты тензора напряжений;
ij
,2
2201211012
деформаций;
0, 0  – постоянные Ламе;
венное число, если
0
, то (367) есть классический закон Гука для изо-
0
– произвольное вещест-
0
тропного, однородного случая.
Подставляя (367) в уравнения равновесия, то в случае отсутствия массовых сил получаем
(2) ( )
uuvv   
0
xxyxxxy
uuv v
00
xy yy xy yy
() 0,
 
(368)
uuv v
00
u,
где
()
yxyxy xy
(2) ( ) 0,
vuvu   
yy yy xy xy
0
- компоненты вектора деформаций, индексы внизу, если не указано
иное, означают производные по соответствующим переменным.
Групповые свойства уравнений асимметричных уравнений упруго­сти в динамическом случае изучены в [7], там, в частности, показано, что система (368) допускает бесконечную группу точечных преобразований, порождаемую операторами
где
Xh h
12
,hh
– произвольное решение уравнений Коши–Римана
12
 
uv
,
(369)
12
hh
xy
12
0,
hh
yx
0.
(370)
Сделаем групповое расслоение системы уравнений (368) по методу [21] на подалгебре, порождаемой (369). Для этого продолжим операторы (369) на первые производные. Имеем
1
XXh h
12
uyv
12
hh

yu x v
y
yx
(371)
,
Дифференциальные инварианты для (371), с учетом (370), имеют вид
,Ix
1
,Iy
2
uv
3
,
y
uv
4
.
yx
(3.72)
Тогда автоморфная система уравнений имеет вид
uv xy (, ).
xy
(, ),
87
uv xy
yx

(373)
Напомним некоторые свойства автоморфных систем. Любое реше-
x
x
A
B
ние автоморфной системы может быть получено из одного решения этой системы с помощью преобразований, порождаемых оператором (369).
Подставляя (373) в (368) получаем разрешающую систему
(2) 0,
F       
100
xyy x
(374)
(2) 0,
F 
200
yx yx
Повторяя почти дословно рассуждения из [7], можно утверждать, что система (368) равносильна системе уравнений (374).
Поэтому, построив решение системы (374), мы получим решение системы (368).
Пусть для системы (374) поставлена следующая краевая задача:
|(,),
L
y
|(,),
0

L
y
0
(375)
где L – некоторая гладкая замкнутая кривая.
Для решения этой задачи построим законы сохранения для системы уравнений (374).
В силу линейности системы (374) она будет иметь бесконечное число законов сохранения. В работе будут найдены только те законы сохранения, которые позволят решить краевую задачу (375).
Определение. Законом сохранения для системы уравнений (374) назовем выражение вида
(, ,, )
Axy
x
(, ,, ) 0,
Bxy F F
y
  
12
(376)
где , – некоторые функции, которые не равны тождественно нулю од- новременно;
,AB – компоненты сохраняющегося тока.
Более подобная информация по построению законов сохранения для произвольных систем дифференциальных уравнений может быть найдена в гл. 2 или в работах [8–12].
Предположим, что компоненты сохраняющегося тока имеют вид
где
12
,,aa
12
,bb
– некоторые функции от x, y.
Тогда из (376) получаем
1
(2) ,a    
2
Отсюда получаем
 0.
12
aa
0
xy
2
  
11
,a    
ab
xy
0
0,
,
88
12
bb

1
,a
0
22
0,
ab
xy
 
yx
,
(377)
(2) ,b  
0.
(378)
0
(379)
Из (376) следует
A
x
A
x
A
A
x
A
x
x
Пусть
где
(, )
yS
00
() .

SL
(, )
yS
00
:( ) ( )xx yy  
222
00
, тогда из (380) следует
 

L
, если функции A,B имеют особенности в этой точке. Вычислим
Bdxdy Ady Bdx
xy
dy Bdx Ady Bdx
– окружность радиуса вокруг точки
(380)
,
(381)
интеграл в правой части (381) для разных решений уравнений Коши­Римана. Пусть
xx

()( )
xx yy
dy Bdx


00

0
22
((( 2) ) ( ))dy
;

()( )
xx yy
00
yy
0
22

00
;
(382)
(383)
(( ) (( 2) )dx    .
00
Подставляем (382) в (383) и сделаем замену переменных
cos ,xx 
0
2
[(( 2 ) ( )) 2sin cos( ))]d
           
0
2[( 2 )( , ) ( ) ( , )]
В формуле (384) устремили
sin ,yy
0
00 0
получаем
dy Bdx
 
000 0 00
yxy       
0
и использовали теорему о среднем.
.
(384)
Теперь сделаем аналогичные вычисления, положив

0
()()
xx yy
22

00
;

()()
xx yy
xx
0
22

00
.
yy
В результате получим
dy Bdx
 
2( )( , ).
 
000
Формулы (384) и (385) позволяют, с учетом граничных условий (375) и равенства (381) определить значения функций
(, )
ке
yS
00
.
89
y
(385)
в произвольной точ-
и
Теперь после восстановления решений разрешающей системы, най-
A
B
A
A
A
A
дем решения автоморфной системы, т. е. решения исходной системы урав­нений (368). Имеем
Fuv xy
3
xy
(, ) 0,
Fuv xy
4
Здесь в правой части стоят известные функции. Найдем законы со­хранения уравнений (386) в следующем виде
где
34
,,aa
3412
,,,bbcc
– некоторые функции от x, y.
341
aa c
,
34 2
bb c
Имеем
(, ,,)
xyuv
x
(, ,,) 0.
Bxyuv F F
y
Расщепляя систему уравнений (388) получаем
3
33
ab
xy
4
,a 
44 12
0,
ab cc
xy xy
3
,a
0, .

Отсюда получаем
 0.
xy
0,
 
yx
Пусть для системы (368) поставлена следующая краевая задача:

yx
 
 
34
4
,b
(, ) 0.
,
(387)
(386)
(388)
,b
(389)
(390)
|(,),
uuxy
0
L
|(,),
vvxy
0
L
(391)
Рассмотрим закон сохранения в виде
dy Bdx 
L

dy Bdx
.
(392)
Пусть решение уравнений (390) имеет вид
xx

()( )
xx yy

00
0
22
,

()( )
xx yy
yy
0
22

00
Подставляем (393) в правую часть (392), получаем
dy Bdx
 
(cos sin ) (sin cos ) 2 ( , ).cdy c dx ux y
      
()( )uvcdy u vcdx
     
12
12
Пусть решение уравнений (390) имеет вид
,
(393)
00
(394)
yy


()( )
xx yy
00
0
22
,
90

()()
xx yy
xx
0
22

00
.
(395)