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

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

.pdf
Скачиваний:
0
Добавлен:
08.09.2026
Размер:
2 Мб
Скачать
☆
Перепишем уравнение (624) для функций A и B:
Г

Разобьем границу
Г
– окружность
3

l
 
 
 
l
2
vdy
l
1
2
1
 
 
dy
yl
l
1
l
2
Г
dxuvdyvuBdxAdy
l
1
2
ГГГ
dxvudy 0
 
 
 
 
dxdx
на части (рис. 5.5), т. е.
2

0
2
2
.
Ryyxx
0
l
l
dy
2
1
dxdx
ГГ
x
 
 
xl
y
.
vdy
dx
 
 
 
 
l
1
l
2
dxudx
ГГГГ
(632)
ГГГГГ
, где
4321
Рис. 5.5. Граница
Тогда

  
Очевидно, что
условия уравнение (633) примет вид

Вычислим интеграл
Г
1
, где
Г
dxvudydxvudydxv
ГГ
43
ГГ
ГГ
13
– окружность радиуса R.
1
dxvudydxvudy
dxvudydxvudy
udydxvudydxvudyBdxAdy
ГГГГ
21
.0
.042
С учетом этого
(633)
(634)
161
тогда
Пусть
xx
1

0
Введем полярную систему координат
 
В результате вычислений при
2
sin
cos
0
1
dRdx
dRdy
Г
Аналогично при
1
,
2
yyxx
0
arctga
  
,
,
0
0
cos
,
1

yy
0
xx
0
Rxx
cos
Ryy
sin
,
sin
R
0R
, получим

,
2
,
2
2
.
yy
0
2
0
,
2
yyxx
0
(635)
(636)
.
a
R
.,
yxudxvudy
00
(637)
(638)
udy v dx v x y
Г
В результате из () имеем
1
Г
Зададим кривую

где

t'
,' tf
– соответственно производные функций
Тогда функции по следующим формулам

yxu
,
00
T
1

yxv
,
00
0
1
 
Для получения этих соотношений использованы решения (631) и (635) соответственно.
3
udyvdxvxy
222 00
Г
3
Г
3

T
 
0


0
2
0
,.

00

,,
yxudxvudy
00111
,.
в параметрическом виде:

0
,tfx
,,00yxu
2
0
,ty
,,00yxv
из (640), (641) будут вычисляться
22
ttfxtf
''
2
''
ytxtf
0
ytxtf
0
22
ttfyt
a
2
2
,0 Tt

arctgtaf
'

 
ln'
(639)
(640)
(641)
(642)


tf
yt
xtf
и
0
0

.t
 
dt
,
 
(643)
2
0
2
dtytxtftf
0
.
 
162
Теперь вычисляем значение выражения


в точке

.,
yx
Те точки, в которых (610) больше или равно единице при-
00
22
vu
(644)
надлежат пластической области, а те в которых выражение (610) меньше единицы – упругой.
На основе формул (634), (635) созданы программы, которые позво­ляют с любой точность строить пластические и упругие области в скручи­ваемом стержне [8; 9; 10].
Решение тестовых задач показало хорошее совпадение с известными решениями.
5.10. УПРУГО-ПЛАСТИЧЕСКАЯ ГРАНИЦА В ИЗГИБАЕМОЙ КОНСОЛИ
Рассмотрим консоль постоянного сечения, ограниченного конту­ром Г. Консоль находится под действием сосредоточенной силы P на сво­бодном конце параллельно главным осям (рис. 5.6).
Рис. 5.6. Упруго пластическая граница
в сечении консоли
Пусть компонента тензора напряжений консоли как и в случае чистого изгиба
z
Пусть компоненты тензора напряжений таковы
Тогда оставшиеся компоненты тензора напряжений удовлетворяют уравнениям
xz
yz
zz
xz
,0
распределена вдоль оси
z
xzip
.
l
.0
xyyx
yz
xx
px
.
l
(645)
163
Обычно уравнения (645) заменяют двумя уравнениями совместности
p
xz

1 vl
,
,0
yz
где
– оператор Лапласа;
v
но решается полуобратным методом Сен-Венана.
Запишем систему (645) в терминах вектора деформаций (u, v; w). Для полученной системы будет поставлена краевая задача, и она будет решена с помощью законов сохранения.
Используя формулы, связывающие компоненты тензора напряжений и тензора деформаций, имеем
где
u
e
x
xy
,
u
y
– коэффициенты Ламе;
,02
x
y
v
,0
x
v
e
,02
y
z
e
w
z
zlp
,2
l
(646)
xz
w
u
z
,
x
e
,
1
v
u
x
y
v
yz
z
w
,,
z
w
y
– искомые функции
21
2
от x, y.
Из первых трех уравнений (646) мы получаем
u
x
где постоянные
BA ,
ii
y
могут быть легко выражены через
v
,11xBxzA
w
,
xBxzA
22
z
,33xBxzA
.,,, lp
(647)
Из уравнений (647) мы получаем
3
zBzA
3

2
21
1
,,/
yxVv
VxUzy
x
.0//
xxyyxy
где
2
xzA
3
w
2
VU ,,
– искомые функции.
3
Из соотношения
 
v
u
x
y
3

u
,,
yxxzB
3
y
0
xy
2
находим

Соотношения
u
 
подставляем в предыдущее уравнение, в результате получаем

x
v

y

xxxx
yyyy
xAxA
21
yyxx
Px
l
xBxzAUz
111
xBxzAUz
222
.
,/
./
,,/
yxUz
164
Предположим, что боковая поверхность бруса свободна от напряже-




ний. Это означает, что
,012 nn
(648)

где
,nn
yzxz
– внешняя нормаль к контуру Г. Мы предполагаем, что пласти-
21
ческое течение начинается с внешней стороны боковой поверхности бруса. Условие пластичности Мизеса в этом случае имеет вид
где k – постоянная пластичности. Решая систему (648), (649) получаем
xz
222
,
k
yzxz
1
,
kn
yz
.2kn
(649)
Выбирая в этих соотношениях верхний знак мы получаем следую­щую задачу. Необходимо решить уравнение
При следующих условиях на Г:
 
 
x
1
kn
1
2
yyxx
2
xA
 
,
y
,ax
xyAkn
22
.
(650)
(651)
Замечание. В сечении бруса есть две зоны пластическая и упругая. Если в некоторой точке выполнено условие
,122
yzxz
Она попадает в упругую зону. Остальные точки, включая границу контура Г, принадлежат пластической зоне.
Соотношение в форме
BA
,0
yx
(652)
Назовем законом сохранения (652) для уравнения (10.6) если уравне­ние (652) верно для всех решений уравнения (650). Пусть сохраняющийся ток имеет вид
 
Из (10.6) и (10.8) мы получаем
yx
a
,,,, yxyxyxA

111
yx
11111
.0

yyyyyxyxyxyxxyxxyy
.,,,
yxyxyxB
(653)
Соотношение (653) верно для всех решений уравнения (650), поэто­му из (653) мы имеем
,0
yx
,0
.0
yx
a
1
yx
(654)
Закон сохранения (652) может быть записан с помощью формулы Грина следующим образом
1

165
yxyx
.0
dxdy
Рассмотрим два решения уравнений (654).
.
Г
Первое решение возьмем в виде
1

xx
0
2
0
yyxx
0
1
,
2

0
yy
0
2
,
2
yyxx
0
1
,0
xa
x
11
,
y
 
arctgax
 
yy
0
.
xx
0
Второе решение пусть будет таким
yy
1

0
2
0
1
,
2
yyxx
0

xx
0
2
0
yyxx
0
,
,0
2
1
y
2
xa
.
x
Используя закон сохранения и применяя его к контуру изображен­ному на рис.6, аналогично предыдущему параграфу, мы получаем

1
dxdy

yxyx

2

0
2
2
Ryyxx
0
1
dxdy

yxyx
Вычислим второй интеграл для первого и второго решений уравне­ния (10.10), в результате получим формулы для нахождения

y
состояние в произвольной точке
.,
yx
00
Формулы для

x
,,
yx
00
y
, yx
позволяют найти напряженное
00
, yx
. Это означает, что мы можем
00

x
,,00yx
для каждой точки области определить в пластической или упругой зоне она находится. Изложенный здесь метод позволяет построить границу между упругой и пластической зонами с любой степенью точности на ком­пьютере. Предварительные вычисления подтверждают этот вывод.
5.11. УПРУГО-ПЛАСТИЧЕСКАЯ ГРАНИЦА
ДЛЯ ОБЛАСТЕЙ БОЛЬШИХ РАЗМЕРОВ
В этом параграфе мы ограничимся областями, граница которых гладкая выпуклая кривая.
Рассмотрим область, ограниченную контуром
(рис. 5.7).
Рис. 5.7. Область больших размеров
166
Пусть на Г заданы граничные условия, в этом случае они запишутся


Г
L
L
так
xyx
yxy
,Yml
(655)
,Xml
Кроме того предположим, что нагрузки на границе таковы, что вся граница и близкие к ним точки области находятся в пластическом состоя­нии. Тогда можно ввести переменные
x
k
,2cos
y
k
по следующим формулам
,
,2cos
k
xy
.2sin
(656)
В этом случае условия (655) запишутся так

Эти условия можно записать в следующем виде
,2sin2cos Xmklk
.
,'X
'
Y
.2cos2sin Ymklk
(657)
Тем самым на границе Г мы имеем задачу Коши в пластической об­ласти (рис. 5.8). Решая эту задачу по формулам §9 вдоль
получаем два
семейства характеристик.
Рис. 5.8. Упруго-пластическая область
Для этих характеристик строим огибающую
(линию возврата). Эта огибающая (линия возврата) и будет искомой упруго-пластической границей. Теперь достаточно решить упругую задачу внутри области огра-
ниченной кривой
и упруго- пластическая задача будет полностью решена.
Теперь можно поставить задачу Коши для системы пластичности
(594)–(595).
Пусть на части граничной линии Г, которую обозначим SP заданы
начальные данные:
'
,
0
SP
SP
Выпустим характеристику
ристику
:
constRS
0
из некоторой точки S. Тогда решение задачи Коши
определено в криволинейном треугольнике
'
,
0
constPR
:
0
из точки P и характе-
SPR
(рис. 5.9).
Таким образом, необходимо определить координаты точки
пересечения характеристик. Зная координаты точки
,
определим в ней значение функций
00
167
.
,
yxR ,
и значения
RR
(658)
R
Рис. 5.9. Задача Коши
SPR
Запишем интеграл по замкнутому контуру
. Используя соотношение (596) в теореме Стокса для плоскости [4], заключаем, что этот интеграл равен нулю:
1. Для координаты
Интегралы
 
 
Полагая
CD
PR
и
:
0
dy
 
dx
dy
CD
dx
,1
.0
CdyDdx
имеем:
x
R
dy
CDCdyDdx
dx
SPRSPR
берутся по частям:
:
RS
0
 

PRPR

RSRS
получим граничные условия для системы (604)
,0

::
RSPRSPR SP
 
00
.dx
xx
R
xx
P
xx
S
xx
R
dxxdxdxCctgDdx
,

PRPR

RSRS
.
dxxdxdxCctgDdx
вида:
При таких граничных условиях, окончательное выражение для
координаты
принимает вид:
x
R
SPR
SP
R
PR
'
xx
R
xxCdyDdxCdyDdx
xx
P
'
.
xCdyDdxx
S
,1'
RS
.0

SP
168
(659)
xx
S
,0
xx
R
(660)
2. Аналогично, для координаты
Интегралы
dy
 
D
dx
dyC
 
dx
D
 
dy
и
:
PR
0
dy
CD
dx
dy
dy
:
RS
dx
 
dyC
 
tg
берутся по частям:
0
CctgD
ctg
CtgD
dy
имеем:
y
R
dx
DCdyDdx
dy
SPRSPR
dy
RSRS
tg
PRPRPR
dy
ctg
tg
 
.dyC
yy
R

PRPR
yy
P
y
tg
dy
y
y
ctg
xx
S

RSRS
xx
R
y
 
.
d
d
ctg
Полагая
tg
,
получим граничные условия для системы
,0
(8.10) вида:
'
0
tg
RS
2
,
.0
PR
'
(661)
При таких граничных условиях, окончательное выражение для
координаты
принимает вид:
y
R
SPR
yy
R
yy
P
'
.
yCdyDdxy
S

SP
R
CdyDdxCdyDdx
SP
x
ctg
tg
yy
S
y
yy
,0
R
(662)
Остается только решить задачи (604), (659) и (604), (661). Эти задачи для системы можно свести к следующим задачам для одного уравнения (8.11). А именно, учитывая, что
cos
,
связаны с функцией  так:
'
,
,
sin
2
(663)
'
получим задачи:
2
0
cos
2
'
,0
4
'
0
,
0
(664)
0
и
2
0
sin
2
'
,0
4
'
0
,
0
(665)
0
,
169
Общим решением первой задачи (664) (для координаты

функция:

1
2
причем
1
1 4
0
1
cos
2


0
Решением второй задачи (664) (для координаты функция:

1
2
причем
2
1 4
Во всех решениях фунция типа мнимого аргумента,
Функции
0
1
sin
2


0
,
восстанавливаются по формулам (604). Из
) является
x
R

I


00
00


I
2
01

I


00
00


I
2
01

I
0
,10
sin
cos,
0001
0
2
001
sin
2
sin,
0002
0
0
0
cos
2
001
0
cos
0
I
является функцией Бесселя первого
0
'
.00
I
0
2
00
2
'
,
dI
0
0
0
.
dI
) является
y
R
00
(666)
(667)
2
'
,
dI
0
0
0
.
dI
(668)
(669)
соотношений (604) определим компоненты законов сохранения:
D
C
ctgtg
ctgtg
ctgtg

22
,cossin
(670)
.cossin
Подставляя найденные C и D в (660) и в (662), определим координаты точки
R
. Таким образом, задача Коши для системы пластичности (594), (595) с начальными данными (658) полностью и конструктивно решена.
170