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

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

.pdf
Скачиваний:
0
Добавлен:
08.09.2026
Размер:
2 Мб
Скачать
☆
допускает, в смысле С. Ли, алгебру Ли размерности 10. Показано, что про-
x


A
A
странственное двумерное деформированное состояние это есть суперпози­ция плоского напряженного состояния и пластического кручения вокруг оси z. Построены два инвариантных решений уравнений, описывающих пространственное двумерное деформированное состояние.
Первое решения можно использовать для описания пластических те­чений между двумя жесткими плитами, которые сближаются с разными скоростями.
Второе решение служит для описания напряженно-деформи­рованного состояния материала внутри плоского канала, образованного сходящимися плитами.
В название пункта вынесено понятие «пространственное двумерное решение». Механикам известно плоское деформированное состояние – это случай, когда в двумерной декартовой системе координат две компоненты вектора скорости деформации и гидростатическое давление зависят от
., y
Известно также плоское напряженное состояние – это когда компо-
ненты тензора напряжений не зависят от
z
.
,,
равны нулю, а компоненты
yzxzz
,,
xyyx
В нашем случае все компоненты тензора напряжений не зависят от
z
, именно такой случай мы и назвали пространственным двумерным со-
стоянием.
Система пространственных уравнений пластичности в декартовой системе координат
,1xx
,2yx
zx
в стационарном случае имеет вид
3
, + (453)
здесь
,1uu
рицательная функция; ks – постоянная пластичности;
s,
– компоненты тензора и девиатора тензора напряжений;
ijij
wu
,2vu
3
– компоненты вектора скорости деформаций; – неот-
– гидростатическое
давление, по повторяющимся индексам проводится суммирование.
Исключая из системы уравнений (453) компоненты девиатора тензо­ра напряжений, получаем следующую нелинейную систему уравнений
2
k
2
p
u
2
2
k
3
22
.,0,
eeAuuee
(454)
mnmniinmjmnijijji
Известно, что система уравнений (454) имеет эллиптический тип. Опишем известные решения этой системы.
Решения этой системы построены Р. Хиллом в 1948 г. [1], В. Праге­ром в 1954 г. [2] и Д. Д. Ивлев в 1960 г. [3; 4], М. А. Задояном в 1964 г. [6–8], а также авторами данной книги [9–13]. Отметим также ряд точных решений, построенных Б. Д. Анниным [14], для уравнений пластичности в пространственном случае с условием текучести Треска.
111
Симметрии системы (454). Группа точечных симметрий системы
x
x
Z
уравнений (2.2) порождается следующими операторами
преобразования
X
,
ix
Tx x
12 3
Zx x
12 3xx
Операторам
',
xa ',
iii
uu
32
32
Zx x
31 2xx
XYNM
uub 'exp,
,
Y  ,
iu
i
uu
23
21
,,,
ii
iii
Nx
i
,
Tx x
23 1
uu
32

соответствуют следующие непрерывные
ii
,
Mu
ix
i
uu
13
,
Zx x
23 1xx
uu S

12
uu p
21
xa'exp,
iu
,
Tx x
31 2

13
,.
uu b
ii
1, 2, 3i
i

uu
21
uu

31
(455)
,
,
uu
13
1, 2, 3 .i
Это переносы по координатам и компонентам вектора скоростей де­формации, а также растяжения.
T
Преобразования няется при жестких перемещениях
',uuс x
2213
',uuс x
3321
Группы порождаемые операторами динатных осей
'cosxx
22 1
'cosxx
33 2
'cosxx
11 3
показывают, что система уравнений (454) не ме-
i
',uuс x
3312
',uuс x
2231
sin ,x
31
sin ,x
12
sin ,x
23
'sinxx
321
'sinxx
132
'sinxx
213
',uuс x
1123
'.uuс x
1132
– это вращения вокруг коор-
i



cos ,x
31
cos ,x
12
cos .x
23
Последнее преобразование описывает инвариантность гидростатиче­ского давления относительно сдвигов
Во всех этих формулах
'.ppd
,,, ,
abc d
iii i
групповые параметры, обычно
предполагается, что они непрерывно меняются в окрестности нуля.
Проведенные исследования показали, что все построенные решения Р.Хиллом, В. Прагером, Д. Д, Ивлевым и М. А. Задояном суть инвариант­ные, решения относительно некоторых одномерных подгрупп точечных преобразования порождаемых операторами (455). Инвариантность здесь означает, что решения не меняются при некоторых преобразованиях, по­рождаемых группой симметрий (455). Так решение Р. Хилла инвариантно относительно подалгебры порожденной оператором
2CS X Y T
0111
решение Д. Д. Ивлева инвариантно относительно подалгебры
,
112
2CS X Y
X


011
aS X T T
11 2
, решение Прагера инвариантно относительно подалгебры
, решение Задояна инвариантно относительно этой же
подалгебры. О чем говорит этот факт? Он говорит о том, что фактически все эти решения «двумерны», т. е. в подходящей системе координат их можно записать как функции только от двух независимых переменных. Это же можно сказать и о решениях построенных авторами этой работы. Тогда возникает вопрос: что надо понимать под пространственным реше­нием? Исходя из приведенных здесь решений, ответ получается такой: пространственные решения это такие решения, которые имеют три компо­ненты вектора скорости, давление, фактически зависящие от двух переменных в подходящей системе координат. Эти решения являются инвариантными решениями ранга 2. В этом случае задача поиска про­странственных решений может быть формально поставлена так: построе­ние новых инвариантных решений ранга 2 для пространственных уравне­ний идеальной пластичности. Вид таких решений можно без труда перечислить, если перечислить все различные, с точностью до внутренних автоморфизмов, одномерные подалгебры алгебры (455).
Подалгебры на которых можно построить инвариантные решения ранга 2 существует несколько. Перечислим их
,, ,, ,
SXZ SMNSNYSZ NYS
       
311 11 1
ZNMSXZMS
    
111
XY TSXZTSMNTS
  
11 1 111 1
,,,
,,
(456)
,,

111 1
.
Здесь
ZY T SX X T T S
  
11 1 1 2 2 3
XZY T S
,,
– произвольные постоянные, разным значениям этих постоян-
ных соответствуют не подобные подалгебры.
В этой работе рассмотрим только решения инвариантные относи­тельно подалгебры
3
.XS
Замечание. Другие решения, построенные на подалгебрах (456),
а также вид всех инвариантных решений, которые могут быть построены для системы уравнений (454) можно найти в [11].
Решения, инвариантные относительно этой подалгебры, следует ис­кать в виде
.),(),,(),,(),,( zyxppyxwwyxvvyxuu
(457)
Подставляя соотношения (457) в систему уравнений (454) получим
,0
xyyxx
2
)(
yx
,6)(6)(
2
)(
zx
,0
yyxyx
,
yzyxzx
,0 vu
yx
k
(458)
22222
syzxzxyzy
,0
,upxx
,vpyy
p
z
113
),(2 uv
yxxy
,2 w
xxz
.2 w
yyz
Сделаем замену переменных в (458) по следующим формулам




),2cossincos3(
),2cossincos3(
x
k
y
k
k
xy
(459)
,2sinsin
Здесь
kk
s
10,
– некоторая постоянная.
Подставляя эти соотношения в(2.6) получаем
)2cossin3cos(
x
)2cossin3(cos
y
,
yzyxzx
,upxx
),(2 uv
yxxy
,0 vu
yx
,vp
yy
,2 w
xxz
2sinsin3
y
2sinsin3
x
syzxz
p
z
x
Kkk
,0
.2 w
yyz
,0sin2
,0sin2
(460)
y
22222
,
Из (460) видим, что исходная система распалась на две подсистемы: два первых уравнения по сути совпадающими с уравнениями, описываю­щими плоское напряженное состояние
)2cossin3cos(
x
)2cossin3(cos
y
2sinsin3 ,0sin2
y
2sinsin3
x
y
x
(461)
,0sin2
и уравнениями, напоминающими уравнения, описывающими пластическое
0
кручение стержня, при
и другом пределе текучести
,
yzyxzx
syzxz
22222
Kkk
.
(462)
Решая уравнения (461) и (462), можно найти компоненты тензора на­пряжений, при этом
деление постоянной
p
.k
Единственная проблема – это опре-
).(2/1
yxz
Для определения компонент вектора скорости получаем следующие уравнения
2
)(
u
x
yx
w
x
x
ww
v
y
2
22
yx
xy
w
y
y
vv
yx
6
ww
,
xy
.)(
22
yx
Отметим, что последнее уравнение при правой части не равной нулю до их пор недостаточно исследовано, оно даже не вошло в справочник [15].
В итоге получилось, что решение системы уравнений (462) это фак­тически суперпозиция плоского напряженного состояния и пластического кручения вокруг оси
.z
Приведем некоторые другие решения уравнений (462). Для этого найдем группу точечных симметрий системы уравнений (462).
114
Эта группа порождается операторами





,



xxT
21 uu
,
Y
,
12
X
1 x
,
X
2 y
Ищем решение инвариантное относительно подалгебры, порождае­мой операторами
Решение следует искать в виде
Из уравнения несжимаемости получаем
Подставляем (464) в (462) и получаем систему обыкновенных диф­ференциальных уравнений
d
U
dy
d
,)'(
dy
Из (465) имеем
,
ui
i
2
xN
x
xxZ
213 xx
.
SYYTX
311
s
22
,
y
y
12
,
uM
ui
i
21 puu
Suu
12
.)(2/
Cyyxv
1
',' CWCyU
.
32
3,2,1i
(463)
.,
).(),(),(2/),(
yPxpyWxwyVxvyUxyxu
(464)
222222
).)'()'((6)1(26,0)'(
WUykW
(465)
,, CCC
Здесь
– произвольные постоянные.
321
Получаем
,/)('/'
CCyWU
32
s
2222
4
22
./,/),)')()(1((6)1(26
CCCCWyCyk
3243
Систему уравнений (465) удается свести к квадратурам, которые вы­ражаются через эллиптические интегралы.
Построенное решение можно использовать для описания пластиче­ского течения слоя, сжимаемого жесткими плитами ортогональными оси oz.
),//( rvrup
22222
,6)(6)(
srzzrz
zr ,,
(4.66)
Запишем систему (454) в цилиндрической системе координат
1
1
z
,2 w
rzr
rzzrzr
,0p
2
)(
r
rrrr
,0/1
rr
,0/)(
rr
)(
zr
rr
rr
2
,up
)),/(/(2 rvrru
z
rrr
1
,0/2
rr
,2
wr
k
0//
u
r
Ищем инвариантное решение на подалгебре
),(Uu
),(Vv
115
rvru
.
MZ
, оно имеет вид
3
.ln)(),( rPpWw
(467)
Подставляя (467) в (466) получаем систему обыкновенных диффе-

ренциальных уравнений
d
rr
r
2
)(
),/'/( rVrUp
z
2
)(
zr
p
z
1
Wr
d
,)(
Здесь штрих означает производную по
,02
r
,'2
VU
d
,0
.0'
.
zr
,/2 r
,0
rzr
r
rr
22222
,6)(6)(
k
srzzrz
,0p
)),/(/'(2 rVrrU
(468)
Отсюда получаем систему обыкновенных дифференциальных урав­нений
s
2222
.1)'()'(6,0)''(,)'(,))''((
WVUkWPVUVU
(469)
Система (469) решается полностью аналогично системе (466).
Найденное решение можно использовать для описания пластическо­го течения в сходящемся плоском канале с жесткими и шероховатыми стенками.
Другие решения системы уравнений (469) можно найти в [11].
Вывод. В этом параграфе изучен класс уравнений, который назван
уравнениями, описывающими пространственное двумерное деформиро­ванное состояние. Для этих уравнений найдена группа точечных симмет­рий, допускаемая ими в смысле Ли. Показано, что двумерное напряженное состояние это есть суперпозиция плоского напряженного состояния и кру­чения вокруг оси z. Построены несколько инвариантных решений этих уравнений.
4.3. НОВЫЕ КЛАССЫ РЕШЕНИЙ ДИНАМИЧЕСКИХ ЗАДАЧ
АНИЗОТРОПНОЙ ТЕОРИИ ПЛАСТИЧНОСТИ
Динамические задачи – это наименее изученная область теории пластичности. Динамические задачи возникают в самых разных областях техники и науки, но сложность исходных дифференциальных уравнений не позволяет строить точные решения и корректно численно решать крае­вые задачи. Это еще в большей степени касается динамических уравнений анизотропной пластичности. Анизотропия уменьшает группу симметрий, допускаемую уравнениями, а, следовательно, и сужает количество инвари­антных решений. Неплохо исследованы одномерные динамические задачи пластичности, но уже двумерные задачи вызывают непреодолимые мате­матические сложности, из-за нелинейности основных уравнений, даже в изотропном случае. Изучение симметрий уравнений пластичности позволило построить некоторые точные решения. Наиболее известное ре-
116
шение построил Б. Д. Аннин, описывающее сжатие жесткими плитами пластического слоя из изотропного материала. Решение Аннина линейно по двум пространственным переменным, и в него входят произвольные функции времени. В предлагаемой работе также используются симметрии. В статье впервые вычислены точечные симметрии для динамических урав­нений пластичности в анизотропном случае. Алгебра Ли, порождаемая найденными симметриями оказалась бесконечномерной. Это обстоятель­ство дало возможность применить методику построения новых классов нестационарных решений. Симметрии позволяют преобразовать точные решения стационарных динамических уравнений в нестационарные реше­ния. В построенные решения входят произвольные функции и произволь­ные постоянные. План работы следующий.
1. По методике Ли-Овсяннникова вычисляется группа точечных
симметрий, допускаемая уравнениями анизотропной пластичности.
2. Строятся два класса новых стационарных инвариантных решений. Эти стационарные решения, с помощью преобразований, порождаемых точечными симметриями, преобразуются в новые нестационарные ре­шения.
3. В заключение работы построено новое автомодельное решение не­стационарных уравнений анизотропной пластичности, а решение Аннина обобщено на анизотропный случай. Приведенные в работе решения можно использовать для описания сжатия пластического материала между жест­кими плитами, а также для тестирования программ, предназначенных для исследования анизотропных пластических задач.
В целом, теории пластичности посвящена обширная литература.
Это объясняется важностью и актуальностью рассматриваемых там задач.
нологических процессах использующих пластические деформации, при пробивании брони снарядом. В современных и классических статьях и мо­нографиях [1–8], рассматриваются в основном статические задачи и изо­тропные материалы. Это связано не с тем, что динамические задачи не важны для приложений, а с тем, что не разработаны методы для решения динамических задач. Впервые пространственное решение динамических уравнений построил Б. Д. Аннин [11] в 1978 г. Это решение было линейно по пространственным переменным и содержало несколько произвольных функций зависящих от времени. Это решение было построено Б. Д. Анни­ным на основе изучения группы точечных симметрий, допускаемых системой уравнений динамической теории пластичности. Позднее на осно­ве групповых свойств уравнений были построены точные решения некото­рых плоских динамических задач [12]. С тех пор, до работы [15], новых решений видимо не появилось. Здесь авторы ссылаются на уникальный обзор [14]. Оказалось, что анизотропным динамическим уравнениям пла­стичности можно применить новый способ – преобразовать стациона
рные
решения идеальной пластичности в нестационарные, так же, как в [15].
117
Групповые свойства можно использовать для разных целей. Чаще




всего они применяются для построения инвариантных решений. Это такие решения, которые не изменяются при непрерывных преобразованиях, до­пускаемых, в смысле Ли данной системой дифференциальных уравнений. Более подробно с инвариантными решениями уравнений пластичности и способами их построения можно ознакомиться, например, в [11] и цити­руемой там литературе. В работах авторов [16–18] показано, как можно с помощью точечных симметрий деформировать точные решения и реду­цировать одно точное решение в другое точное решение в случае плоских стационарных уравнений идеальной пластичности. В данной работе мы используем группу точечных симметрий для преобразования новых стационарных решений в новые нестационарные для трехмерных неста­ционарных уравнений пластичности. Впервые это было сделано для дру­гих решений в [15].
Постановка задачи. Пусть x = x
декартова система координат, u =
, y = x2, z = x3 – ортогональная
1
, v , w – компоненты вектора
скорости деформаций,
– компоненты тензора скоростей деформации,
– компоненты тензора напряжений. Компоненты тензора напряжений
и компоненты вектора скорости удовлетворяют уравнениям движения
ivvv
ijiijjit
3,2,1,
По повторяющимся индексам предполагается суммирование. Девиа-
тор тензора напряжений и тензор скоростей деформации сосны
2/)(
vvep
,. ijjiijijij
где
– символ Кронекера; – некоторая неотрицательная функция,
.
Среда предполагается несжимаемой, поэтому имеет место уравнение
несжимаемости
. (472)
Система уравнений (460)–(462) замыкается условием пластичности
Мизеса

1111
2
2222
2
3333
2
2
1212
2
1313
2
12
aaapapapa
2323
,
(470)
(471)
(473)
где
– текущие параметры анизотропии.
Групповые свойства уравнений динамической теории пластичности
Вычислим группу точечных симметрий допускаемых уравнениями (470)–(473). Сделаем это по известной методике Ли-Овсянникова, которая уже применяется к уравнениям пластичности в [11]. Она порождается сле­дующими операторами:
,)(,,
0 pxitt
118
i
tSxtMX
(474)
iii
,)()()(
tfxtftfT
piivixii
по i не суммировать,
Функции являются произвольными из класса . Поэтому операторы (474) порождают бесконечномерную алгебру Ли. Точка вверху
означает производную по переменной t.
Замечательным свойством точечных симметрий является то, что они решение системы (470)–(473) снова переводят в точные решения этой же системы. Пусть
исходные координаты, тогда они преобразуются
в новые координаты
с помощью преобразований, соответст-
вующих операторам (474)
, (475)
здесь
– групповые параметры, которые непрерывно меняются в некото-
рой окрестности нуля. Далее эти преобразования используем для построе­ния новых решений системы (470)–(473).
Стационарные решения системы (470)–(473).
Поскольку система (470)–(473)допускает оператор
то можно
искать инвариантные решения этой системы, не зависящие от перемен­ной t. Такие решения определяются из системы
2/)(

1111
2
2222
2
2
3333
vvep
2
1212
, ,
,. ijjiijijij
2
1313
2
12
aaapapapa
2323
(476)
Система (466) проще, чем исходная система уравнений, поскольку имеет меньше независимых переменных, но ее решений, кроме тривиаль­ных, насколько известно авторам, нет [1–3].
А. Будем искать инвариантное решение системы (476) относительно
двумерной подалгебры, порождаемой операторами Это решение следует искать в виде
(477)
Здесь A, B произвольные постоянные, а функции u, v, w, p определя­ются из системы (476). Подставляем соотношения (476) в (477). Имеем
(478)
Из (468) получаем систему обыкновенных дифференциальных урав­нений для определения функций v, w
.
(479)
119
Введем новые искомые функции по формулам
Тогда система (479) запишется так
Из (481) получаем
Интегрируя это уравнение, получаем
. (480)
. (481)
где
.
Отсюда получаем
,
h =
.
,
(482)
Исследуем поведение функций, входящих в формулы (482).
Пусть монотонно возрастает от нуля до тонно возрастает от нуля до
меняется от до в этом случае
, ) моно-
. При этом изменяется от до
Теперь, в силу формул (482) становится понятно, как меняются компонен­ты тензора напряжений.
Это небольшое исследование позволяет интерпретировать получен­ное решение следующим образом. Имеются две жесткие шероховатые плиты
. Пластический материал продавли-
вается между ними. На плитах заданы касательные напряжения
Б. Будем искать инвариантное решение системы (476) относительно
одномерной подалгебры, порождаемой оператором решение будем искать в виде
(483)
120
Это