Сапунов Теория пластичности 2011 (1)
.pdfВ |
рассматриваемой задаче напряженияsr и |
||||
sj |
являются главными |
( trj = 0 ), причем |
|||
sj > sr . Соответственно, будем иметь: |
|||||
|
q = a - |
p |
= j + |
p |
. |
|
4 |
|
4 |
Рис. 2.15 |
|
|
|
|
|
|
Подставим полученное значение q в уравнение для линии скольжения a . После некоторых преобразований получим:
|
|
|
|
d y |
æ |
|
p ö |
Þ |
|
d y |
|
1 |
+ tg j |
. |
||||||||||
|
|
|
|
|
|
|
= tg ç j + |
|
÷ |
|
|
|
= |
|
|
|
|
|
|
|
||||
|
|
|
|
d x |
4 |
|
|
|
|
|
- tg j |
|||||||||||||
|
|
|
|
è |
|
ø |
|
|
d x 1 |
|
||||||||||||||
Поскольку y = r sin j , |
x = r cos j , то |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
d y = dr sin j + r cos jdj , |
d x = dr cos j - r sin jdj , |
|||||||||||||||||||||
и уравнение для линии скольжения a принимает вид |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
dr sin j + r cos jdj |
= |
1+ tg j |
. |
|
|
|
|||||||||||
|
|
|
|
|
|
|
dr cos j - r sin jdj |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
1- tg j |
|
|
|
||||||||||||
Выделяя слагаемые при dr |
и dj , будем иметь: |
|
|
|
|
|
|
|||||||||||||||||
|
1 |
[(1- tgj sin) j - 1+ tg(j cos j]dr) = - [(1+ tgj sin) j + 1 - tg(j cos j]dj) |
||||||||||||||||||||||
|
|
|||||||||||||||||||||||
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
dr |
= - |
[(1 + tgj)sin j + (1 - tgj)cos j] |
|
dj . |
|||||||||||||||||
|
|
|
|
[(1 - tgj sin) j - 1 + tg(j cos j] |
||||||||||||||||||||
|
|
|
r |
|
) |
|
После некоторых преобразований в квадратных скобках и интегрирования уравнение, определяющее линию скольжения a , принимает вид
ln r = j Þ r = Aej . A
где A - постоянная интегрирования.
41
Аналогичным образом получим уравнение для линии скольженияb . Будем иметь:
r = Ae -j .
Таким образом, в задаче о нагружении толстостенной трубы внутренним давлением линии скольжения представляются семейством взаимно ортогональных логарифмических спиралей.
Перейдем к определению напряжений и нагрузки в предельном состоянии. Решение поставленной задачи здесь можно существенно упростить, принимая во внимание, что касательное напряжение trj равно нулю, а напряжения sr и sj
являются функциями только переменной r , причем sj > sr . В этом случае ус-
ловие пластичности (2.13) принимает вид
(sr - sj )2 + 4t2rj = 4t2т Þ sj - sr = 2 tт .
Добавив к полученному уравнению условие, выполняемое вдоль линии скольжения a , а именно:
s |
- q = const = x |
q =j + (p / 4) |
sr + sj |
= ln |
r |
+ x + |
p |
, |
|
|
¾¾¾¾¾¾® |
|
|
|
|||||
2tт |
4tт |
A |
4 |
||||||
|
|
|
|
|
где x - постоянная вдоль линий скольжения a , будем иметь два уравнения, ко-
торые рассматриваем как уравнения с двумя неизвестнымиsr и sj . Отметим,
что вместо условия, выполняемого вдоль линии скольжения a , можно взять аналогичное условие вдоль линии b .
Решение представленной системы уравнений позволяет получить:
æ |
|
r |
|
|
p |
|
|
1 |
ö |
|
sr = 2tт ç ln |
|
|
+ x + |
|
- |
|
÷ , |
|||
A |
4 |
2 |
||||||||
è |
|
|
|
|
|
ø |
||||
æ |
|
r |
|
|
p |
|
|
1 |
ö |
|
sj = 2 tт ç ln |
|
|
+ x + |
|
|
+ |
|
÷ . |
||
A |
|
4 |
|
2 |
||||||
è |
|
|
|
|
|
ø |
||||
Поскольку напряжения sr и sj |
являются функциями только переменной r , |
|||||||||
величина x постоянна для всех линий скольжения, |
т.е. во всей области попереч- |
ного сечения цилиндра. Выполняя граничное условие на внутреннем контуре при r = a , будем иметь:
42
sr r =a = - p пр
Þ
Þ |
|
æ |
a |
|
|
|
p |
|
|
1 |
ö |
= - p пр |
Þ |
||||
2tт ç ln |
|
|
+ x + |
|
|
- |
|
|
÷ |
||||||||
A |
4 |
2 |
|||||||||||||||
|
|
è |
|
|
|
|
|
ø |
|
|
|||||||
x = - |
p пр |
- ln |
|
a |
- |
p |
+ |
1 |
|
. |
|
|
|||||
2tт |
|
A |
4 |
2 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
С учетом найденного значения постоянной x напряжения принимают вид
|
r |
|
æ |
r |
|
ö |
|
sr = - p пр + 2tт ln |
|
, |
sj = - p пр + 2 tт çln |
|
+1 |
÷ . |
|
a |
a |
||||||
|
|
è |
|
ø |
Выполняя теперь граничное условие на внешнем контуре приr = b , получим соотношение, определяющее предельное давление p пр :
sr |
|
r =b |
= 0 |
Þ p пр = 2 tт ln |
b |
. |
|
||||||
|
|
|||||
|
|
|
|
a |
Окончательно, для напряжений будем иметь:
|
r |
|
æ |
r |
|
ö |
|
æ |
r |
|
1 ö |
||
sr = 2 tт ln |
|
, |
sj = 2 tт ç ln |
|
+1 |
÷ |
, |
sz = 2tт ç ln |
|
+ |
|
÷ . |
|
b |
b |
b |
2 |
||||||||||
|
|
è |
|
ø |
|
è |
|
ø |
Отметим, что полученные соотношения для напряжений полностью совпадают с теми, которые следуют из решения упругопластической задачи о нагружении толстостенной трубы внутренним давлением для случая, когда материал не имеет упрочнения.
2.4. Тело с цилиндрической полостью, растягиваемое на бесконечности (решение Л.А. Галина).
Рассмотрим |
задачу |
о |
растяжении |
|
||
направлении осей x и y |
бесконечного |
|
||||
массива |
с |
цилиндрической |
поло |
|
||
(рис. 2.16). В направлении |
оси x |
массив |
|
|||
растягивается |
усилиями |
p |
(при |
r ® ¥ |
|
|
sx ® p ), в направлении оси y - усилия- |
|
|||||
ми q ³ p (при r ® ¥ sy ® q ). |
|
|
||||
Отверстие свободно от нагрузок. |
|
|||||
При небольших значениях усилий p и |
|
|||||
q массив целиком находится в упругом |
Рис. 2.16 |
43
состоянии, и решение находится достаточно просто методом Колосова- Мусхелишвили1 с использованием функций комплексного переменного.
При достаточной величине нагрузок p и q вокруг цилиндрического отверстия возникают пластические области (зоны). Не вникая в историю постепенного развития этих областей, рассмотрим случай, когда пластическая зона I полностью охватывает отверстие, тогда как остальная часть телаII остается в упругом состоянии. Поверхность, разграничивающую области I и II , обозначим через C .
Решение поставленной задачи, как задачи о плоской деформации, сводится к изучению напряженного состояния в упругой и пластической областях в плоскости x y и определению формы кривой C , которая заранее неизвестна.
Поскольку пластическая зона I расположена вокруг кругового отверстия, при рассмотрении напряженного состояния в пластической зоне целесообразнее применять полярную систему координат.
Для плоского деформированного состояния при отсутствии объемных сил ре-
шение задачи в полярной системе координатr , |
j |
сводится к интегрированию |
||||||||||||||
двух дифференциальных уравнений равновесия |
|
|
||||||||||||||
|
¶s |
r |
+ |
1 |
|
¶trj |
+ |
sr - sj |
= 0 |
, |
||||||
|
¶r |
r |
|
|
¶j |
|
r |
|||||||||
|
|
|
|
|
|
|
|
|||||||||
|
|
¶t rj |
|
+ |
|
1 ¶sj |
+ |
2t rj |
= 0 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
¶r |
|
|
|
r ¶j |
r |
|
||||||||
|
|
|
|
|
|
|
|
|
|
при выполнении условия пластичности
(sr - sj )2 + 4t2rj = 4t2т .
Совместное решение представленных уравнений достаточно сложно, однако в рассматриваемой задаче оно существенно упрощается за счет того, что границей тела является круговой контур. Рассматривая прямолинейную свободную границу в прямоугольной системе координат (разд. 2.3), показали, что вблизи такой границы всегда имеет место равномерное одноосное растяжение или сжатие, параллельное линии границы.Точно так же легко доказать, что у круговой границы (свободной или равномерно нагруженной) всегда реализуется осесимметричное поле напряжений.
Для осесимметричной задачи касательное напряжение trj равно нулю, а на-
пряжения sr и sj являются функциями только переменной r , причем sj > sr ,
и из трех уравнений остается только два:
________________________________________
1 Сапунов В.Т. Прикладная теория упругости. Ч.2. М.: МИФИ, 2008. – 140 с.
44
d s |
r |
+ |
sr - sj |
= 0 |
, |
|
|
r |
|||
d r |
|
|
sj - sr = 2 tт ,
совместное решение |
которых |
при выполнении |
граничного |
условияsr |
|
r =a = 0 |
||||||
позволяет получить: |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
r |
æ |
r |
|
|
ö |
(2.14) |
||||
sr |
= 2 tт ln |
|
|
, sj = 2tт ç ln |
|
+1 |
÷ . |
|||||
a |
a |
|||||||||||
|
|
è |
|
|
ø |
|
|
|
Отметим, что осесимметричное поле напряжений уже строилось при решении задачи о нагружении толстостенной трубы внутренним давлением, но с использованием свойств линий скольжения. Естественно, полученные там соотношения для напряжений совпадают с представленными при p пр = 0 .
Введем в рассмотрение«пластическую» функцию напряжений F p . Напом-
ним, что для осесимметричной задачи в полярных координатах компоненты -на пряжений через функцию напряжений представляются формулами:
sr = |
1 dF |
, |
sj = |
d |
2F |
. |
||
r |
|
dr |
dr 2 |
|||||
|
|
|
|
|
С учетом найденных значений sr и sj для пластической функции напряже-
ний F p будем иметь:
|
|
|
æ |
2 ln |
r |
|
r |
2 |
ö |
|
F |
p |
= t |
ç r |
- |
|
÷ . |
(2.15) |
|||
a |
|
|
||||||||
|
|
т ç |
|
2 |
÷ |
|
||||
|
|
|
è |
|
|
|
|
|
ø |
|
В упругой области II решение построим с использованием метода Колосова - Мусхелишвили. Бигармоническая функция упругих напряжений Fe может быть
представлена в виде
Fe = Re [ zj(z)+ c(z)] ,
где j(z) и c(z) - аналитические функции комплексной переменнойz = x + iy
( z = x - iy - сопряженная комплексная переменная).
45
Компоненты упругих напряжений при этом вычисляются по формулам:
sx + sy = 4 Re F1(z) ,
s y - sx + 2i txy = 2[z F1¢ (z)+ Y1(z)] ,
где F1(z )= j¢(z) , Y1 (z)= y¢(z)= c¢¢(z).
Для определения аналитических функций F1(z) и Y1 (z) нужно использовать граничные условия на бесконечности и условия непрерывности напряжений на
линии C , разграничивающей пластическую I |
и упругую II |
области. Можем за- |
||||||||||
писать: |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
ì |
2 t |
|
æ |
1+ 2ln |
r |
ö |
- на С , |
|
|
|
|
|
ï |
т |
ç |
|
|
÷ |
|
|||
|
|
|
|
|
|
|||||||
|
|
4 Re F1(z )= í |
|
è |
|
a ø |
|
|
||||
|
|
|
ï |
|
|
|
p + q |
|
|
|
- на ¥ ; |
(2.16) |
|
|
|
î |
|
|
|
|
|
|
|||
|
|
|
|
ì |
|
|
|
-2ij |
|
|
||
2[ |
z |
F¢ |
(z )+ Y (z )]= íï |
|
2tт e |
|
|
- на C , |
|
|||
1 |
1 |
ï-( p - q =)q - p - на ¥ . |
|
|||||||||
|
|
|
|
î |
|
|
|
|
|
|
|
При записи представленных соотношений использовались известные формулы теории упругости:
- формулы преобразования напряжений при замене системы координат
sx + sy = sr + sj ,
s y - sx + 2i txy =(sj - sr + 2i trj )e - 2ij ; - граничные условия на бесконечности
sx + s y = N1 + N2 ,
s y - sx + 2i txy = - (N1 - N 2 )e -2i a ,
где N1 и N2 - значения главных напряжений на бесконечности; a - угол между направлением N1 и осью x .
Отметим, что поскольку линия раздела упругой и пластической зон(кривая C ) неизвестна и подлежит определению наряду с комплексными функциями, соотношений (2.16) недостаточно для решения упругой задачи.
Решение Л.А. Галина строится несколько иначе, используя то обстоятельство, что введенная пластическая функция напряжений F p является бигармонической.
Действительно, можем записать:
46
|
|
|
|
|
é |
|
æ |
z |
|
1 |
öù |
, |
|
F p = tт Re |
ê z z ç ln |
|
- |
|
÷ú |
||||||||
a |
2 |
||||||||||||
|
|
|
|
|
ë |
|
è |
|
øû |
|
|||
æ |
z |
|
1 ö |
|
|
c(z)= 0 . |
|
|
|||||
где принято, что j(z )= tт z ç ln |
|
- |
|
÷ |
и |
|
|
||||||
a |
|
|
|
||||||||||
è |
|
2 ø |
|
|
|
|
|
|
|
|
Введем в рассмотрение новую бигармоническую функцию F = Fe - F p , кото-
рая будет определять новое поле напряжений:
sx = |
¶2 F |
, |
s y = |
¶ |
2 F |
|
, |
txy = - |
¶2 F |
. |
||||||||||
¶y |
2 |
|
¶x2 |
¶x¶y |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Соответственно, можем записать: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
sx + s y |
= |
¶2 F |
|
+ |
¶2 F |
= L (F ) , |
|
|
|||||||||||
|
¶x 2 |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
¶y 2 |
|
|
|
|
|
|||||
s y - sx + 2i txy = |
¶2 F |
|
- |
¶2 F |
|
- 2i |
¶2 F |
= M (F ) , |
||||||||||||
¶x2 |
|
¶y2 |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
¶x¶y |
|
|
где операторы L (F ) и M (F ) введены для упрощения дальнейших записей.
Очевидно, что вместо бигармонической функции упругих напряженийFe
можем разыскивать новую бигармоническую функциюF . Учитывая соотношения (2.16), граничные условия для функции F будут иметь вид:
- на бесконечности
|
æ |
|
r |
ö |
||
L (F )= p + q - 2 tт ç |
1+ 2ln |
|
|
÷ , |
||
a |
||||||
|
è |
|
ø |
|||
M (F )= q - p - 2 tт e -2ij ; |
||||||
- на контуре C |
|
|
|
|
|
|
L (F )= 0 , |
M (F )= 0 . |
|
||||
Отметим, что без операторов L (F ) |
и M (F ) граничные условия для функции |
|||||
напряжений F (для определяющих |
ее |
аналитических комплексных функций |
F2 (z) и Y2 (z )) имеют форму соотношений (2.16):
47
|
|
|
|
|
ì |
|
0 |
|
|
|
|
- на С , |
|
4 Re F2 (z )= íï |
|
æ |
|
|
r ö |
|
|||||
|
|
|
|
|
ï p + q - 2tт ç |
1+ 2ln |
|
÷ |
- на ¥ ; |
|||
|
|
|
|
|
|
|||||||
|
|
|
|
|
î |
|
è |
|
|
a ø |
(2.17) |
|
|
|
|
|
|
|
ì |
|
0 |
|
|
|
|
2 |
[ |
z |
F¢ (z |
)+ Y |
2 |
z (]= )íï |
- p - 2 t |
|
|
|
- на С , |
|
|
2 |
|
ï q |
т |
e -2ij - на ¥ . |
|||||||
|
|
|
|
|
|
î |
|
|
|
|
|
Так же, как и раньше, представленных условий (2.17) недостаточно для решения задачи, поскольку контур C неизвестен.
Проведем конформное отображение упругой области в плоскостиz на плоскость z , так чтобы контур C перешел в окружность g единичного радиуса, а
бесконечно удаленная точка плоскости z - в бесконечно удаленную точку плоскости z . В этом случае отображающая функция w(z) будет иметь структуру:
¥ |
bn |
|
|
z = w (z )= b z + å |
, |
||
zn |
|||
n =1 |
|
где b - вещественная положительная постоянная; bn - коэффициенты, вещест-
венные вследствие симметрии контура C относительно оси x .
Преобразования граничных условий (2.17) при конформном отображении проведем, учитывая следующие соотношения теории упругости:
F2 [w (z)]= F (z) , Y2 [w (z)]= Y (z) , F¢2 [w (z)]= F¢(z)/ w¢(z) .
Соответственно, для определения неизвестных функцийF(z) , Y(z) и w (z) будем иметь:
ì |
|
0 |
|
|
|
|
|
|
4 Re F(z )= íï |
æ 1 |
|
|
b |
|
|
ö |
|
|
|
|
|
|||||
ï |
p + q - 4tт ç |
|
|
+ ln |
|
+ ln |
z |
÷ |
2 |
|
a |
||||||
î |
è |
|
|
|
|
ø |
- на g ,
- на ¥
|
( где ln |
r |
® ln |
|
z |
® ln |
|
bz |
® ln |
b |
+ ln |
|
z |
|
) |
; |
||||||||||
|
|
|
|
|
||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
a |
|
|
a |
|
|
ì |
|
a |
0 |
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
é w(z ) |
|
|
|
|
|
ù |
|
|
|
|
|
|
|
- на |
g |
||||||||||
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
||||||||||||
2 |
ê |
|
|
|
|
F¢(z )+ Y z(ú |
=) |
í |
|
|
- 2 t |
|
e -2ij1 |
|
|
- на |
||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
ë w¢(z ) |
|
|
|
|
|
û |
|
ï q - p |
т |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
î |
|
|
|
|
|
|
|
|
|
|
|
|
( j1 = arg z ) .
(2.18)
,
¥
48
Соотношения (2.18) позволяют определить функции F (z ) , Y (z ) и w (z ).
Не останавливаясь на промежуточных выкладках, приведем окончательный вид функции w (z ):
æ |
k ö |
|
|
q - p |
|
æ q + p |
|
1 |
ö |
||
w (z )= bç z + |
|
÷ |
, |
k = |
|
, |
b = a expç |
|
- |
|
÷ . |
|
|
|
|
||||||||
ç |
÷ |
|
|
2tт |
|
ç |
4 tт |
|
2 |
÷ |
|
è |
z ø |
|
|
|
è |
|
ø |
Как известно, представленная функция отображает внешность эллипса. Следовательно, контур C , разграничивающий пластическую I и упругую II области, является эллипсом. Уравнение эллипса имеет вид
x 2 y 2
b2 (1+ k2 + b2 (1)- k2 =1 . )
Поскольку решение задачи построено в предположении, что пластическая зона полностью окружает круговое отверстие радиусомa , необходимо выполнение
условия: b (1 - k)> a .
Приведем |
некоторые |
результаты вычисления поля напряжений для случая |
||||
p = 2,4 tт и q = 3,0 tт (рис. 2.17). |
||||||
В |
этом случае |
b = 2,34 a и |
||||
полуоси эллипса будут, соответ- |
||||||
ственно, равны 3,04 a |
и 1,64 a . |
|||||
Для сравнения пунктиром пока- |
||||||
зана |
окружность |
|
радиусом |
|||
R = b = 2,72a , |
которая |
является |
||||
линией раздела упругой и пла- |
||||||
стической зон при |
p = q = 3,0 tт . |
|||||
Сплошными |
линиями |
пока- |
||||
заны |
распределения |
интенсив- |
||||
ности |
касательных |
напряжений |
||||
вдоль |
осей x |
и |
y |
для случая |
p = 2,4tт и q = 3,0tт .
Пунктиром - |
распределение |
|
интенсивности |
касательных |
|
напряжений |
вдоль |
радиуса- |
вектора при p = q = 3,0tт . |
Рис. 2.17 |
При заданном значении, например, нагрузки q = 3,0tт можно найти наименьшее значение нагрузки p ( p £ q ), которое следует из условия, что меньшая полуось эллипса будет равна a . В этом случае получим, что p =1,95tт .
49
3.Плоское напряженное состояние
3.1.Общие положения и определяющие уравнения
Взадаче о плоском напряженном состоянии будем рассматривать упругое тело в форме тонкой пластины постоянной толщины
h , нагруженное по боковой поверхности силами, параллельными плоскости пластины и распределенными симметрично относительно ее срединной плоскости, которую совместим с координатной плоскостью x y . Будем считать, что объемные силы отвечают аналогичным условиям.
Торцы (основания) пластины z = ± h / 2 свободны от нагрузки.
Граничные условия на торцах в поставленной задаче принимают вид
sz = tzx = tzy = 0 при z = ± h / 2 .
Если пластину считать тонкой, то при рассматриваемом нагружении с достаточной степенью точности можно принять напряже-
ния sz , |
tzx и tzy равными нулю во всех точках пластины. Ос- |
|||
тальные |
компоненты |
тензора |
напряженийs x , s y и txy |
можно |
считать |
функциями |
только |
переменныхx и y , усреднив |
их по |
толщине пластины. Такое напряженное состояние пластины будем называть плоским напряженным состоянием. Очевидно, что при нагрузке на боковой поверхности, не меняющейся по толщине, процедуру усреднения напряжений можно опустить.
|
Отметим, что при плоском напряженном состоянии имеем |
sz |
= 0 , а ez ¹ 0 (в отличие от плоской деформации, где ez = 0 , а |
sz |
¹ 0 ). Наличие поперечной деформации ez влечет за собой ис- |
кривление плоских оснований пластины, однако, поскольку задача |
симметрична, точки срединной плоскости после деформирования пластины остаются на месте. Данное обстоятельство позволяет утверждать, что при малой толщине пластины перемещение w будет весьма мало и что изменения перемещений u и v по толщине будут незначительны. Соответственно, можно считать, что переме-
50