Math-Net.

Ru
Общероссийский математический портал

Ю. Н. Радаев, Л. В. Степанова, О влиянии удаленной локализованной

пластической зоны на раскрытие трещины нормального отрыва, Вестн.
Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 1999, выпуск 7, 70–85

DOI: http://dx.doi.org/10.14498/vsgtu210

Использование Общероссийского математического портала Math-Net.Ru подразумевает, что вы про-

читали и согласны с пользовательским соглашением
http://www.mathnet.ru/rus/agreement

Параметры загрузки:
IP: 93.75.62.85
4 февраля 2017 г., 20:58:43
 УДК 539.214

Ю. Н. Радаев, Л, В. Степанова

О ВЛИЯНИИ УДАЛЕННОЙ ЛОКАЛИЗОВАННОЙ ПЛАСТИЧЕСКОЙ ЗОНЫ

НА РАСКРЫТИЕ ТРЕЩИНЫ НОРМАЛЬНОГО ОТРЫВА

Нагружение трещины нормального отрыва в условиях плоского напряженного состояния в поле

остаточных напряжений (связанных с предыдущими циклами нагружения) может приводить к
образованию двух очагов пластического течения: непосредственно у кончика трещины и в зоне
максимального остаточного растяжения, которое в случае циклического нагружения достигает
одной трети предела текучести. Моделируя по схеме Дагдейла пластические зоны отрезками, для
определения трех безразмерных параметров, характеризующих положения пластических зон, по­
лучена система нелинейных уравнений, которая анализируется с помощью оригинального численно­
го алгоритма, разработанного специально для этой цели. Получена точная формула для вычисле­
ния раскрытия трещины при двухзонно локализованных пластических деформациях. Асимптоти­
ческий анализ величины раскрытия трещины для случая, когда линейный размер удаленной пласти­
ческой зоны мал по сравнению с длиной трещины, приводит к заключению, что влияние удаленной
пластической зоны на трещину проявляется в форме ее дополнительного закрытия.

1. Постановка задачи и основные уравнения.

Исследование плоского напряженного состояния трещины нормального отрыва бьшо

впервые вьшолнено в работе [1]. В соответствии с экспериментальными данными предполага­
лось, что пластическая зона представляет собой отрезок, находящийся на продолжении линии
трещины. С тех пор предложенная в [1] модель, известная как модель Дагдейла, широко ис­
пользуется для расчета напряженного состояния и вьшисления раскрытия трещины в стальньпс
пластинах [2].
Соответствующая подобной сильно локализовашюй у вершины трещины зоны текучести
модель трещины основывается на следующих основньпс предположениях.
1. Элементы пластической зоны подвержены воздействию только нормального растяги­
вающего напряжения, равного пределу текучести.
2. Толщина пластической зоны намного меньше, чем ее длина, так что внутренняя граница
упругой зоны может рассматриваться как сильно сплющенный эллипс с большой полуосью
равной 2(1 + с), где 21 - длина трещины, с - длина пластической зоны.
3. Длина пластической зоны в точности такова, что пластическая зона полностью погло­
щает сингулярность поля нормальньпс напряжений.
В представляемой работе рассматривается нагружение пластины с трещиной нормального
отрьша, при котором возникает дополнительная, удаленная от кончика трещины, зона пласти­
ческого растяжения. Подобная ситуация, как показывают расчеты, действительно реализуется
при циклическом нагружении с переменной амплитудой напряжений. Разгрузка трещины со­
провождается мгновеюшш образованием у кончика трещины сжатой пластической зоны, в то
время как на некотором удалении от кончика трещины возникает зона остаточного растяжения,
в которой напряжения могут достигать трети от предела текучести. В этой последней зоне в
течение последующих циклов нагружения могут локализоваться пластические деформации,
образуя новый очаг пластического течения.
Вьршсление параметров, определяющих локализацию зон пластического течения, сопря­
жено с решением сложной системы двух нелинейньпс уравнений и поэтому реализовано чис­
ленно. Величина раскрьггия трещины при повторном нагружении вьшисляется с помощью по­
лученной точной формулы, в которую входят только указанные параметры локализации. Эта
формула тем самым позволяет оценить влияние удаленной пластической зоны на величину
раскрытия трещины.
Рассмотрим трещину нормального отрьша < / в пластине в условиях плоского напря­
женного состояния. Обозначим через Y предел текучести при одноосном растяжении.
70
 Напряжения и перемещения представляются через единственный комплексный потенциал
по формулам Колосова-Мусхелишвили [3]:

- 2x,Im ^f(p*\, СГ22 =2Rel5'j«i?*J+2x2lm[^^^«?*| cr,*2 =-2x2Re[^^V* (1)

IGu* = (к-r)Re<p* -2x2lm[^,q>*], Ющ = ( / c + l)Im 9 Pz9 (2)

В которых
Cr*i = С Г ц ,a'*2 = <Ti2 , ^ 2 2 = ^22 - ^22 ' = W| , ^ 3 = ^ 2 "~ ^2^'^^S ' ЛГ = (3 - v)/(l + v) ,

9*(z) - комплексный потенциал, который является аналитической в верхней полуплоскости

Im ^ > О функцией комплексной переменной z = лг^ + ^ 2 ; G - модуль сдвига; v - коэффици­
ент Пуассона; Е - модуль Юнга; о*^ - нормальные растягивающие напряжения на бесконеч­
ности.
Обозначим через а^^., щ напряжения и перемещения, а через = / + с - координату гра­
ницы пластической зоны при первом нагружении напряжением ст^ . В соответствии с решени­
ем Дагдейла находим:

п
arcsm (l + ^)sec/? + arcsm~'-?5^^^^-4l
(l-^)sec/?J ^
>secy^);
^
(3)

TtEui
= -ln(^^ -l)-2lncosy3 + 2ln(sin>^ + ^ / l ~ ^ 2 cos^ ;^) + ^1п|~~ +
2YJ

(4)
sin^7l-i^'co77-^cos^ +1

y = sec;ff-l, (5)

где введены следующие безразмерные величины:

(6)
/ " 27
Для исследования напряженно-деформироваююго состояния после разгрузки пластины
введем следующие разности: AG^J = 5^^ -e^j и Ащ = Ц - Щ, где 5^^ я щ - остаточные напря­
жения и перемещения, которые наблюдаются в пластине с трещиной после разгрузки. Основой
последующего анализа является предположение о том, что после снятия нагрузки у вершины
трещины мгновенно происходит перераспределение напряжений: нормальные растягивающие
напряжения на пределе текучести мгновенно сменяются сжимающими, равными пределу теку­
чести при сжатии. Таким образом, поле остаточньпс напряжений сжимает вершину трещины,
что приводит к ее закрытию. Что касается перемещений, то поскольку снятие нагрузки за пре­
делами локализованной зоны сжатия приводит к обратному переходу из пластичности в упру­
гое состояние, то изменения нормальньпс перемещений вне указанной зоны остаточного сжатия
не происходит.
Исходя из приведенного выше анализа, для определения комплексного потенциала Аф*
граничные условия могут быть сформулированы в виде

Аа^з^а^ ( х , </, ^2 = 0); Аа^з = + (/^^i </ + с , Х 2 = 0 ) ;

71
 Ащ=0 (X, >l + c,X2=0), д,Ац,\2) = 0(2-') (z->c«), (7)

где AGIJ= AG22 + ^22' "^i = ^ + ^ - координата остаточной сжатой пластической зоны.
Предполагается (и это подтверждается полученным решением), что при повторном нагру­
жении пластины растягивающим напряжением образуются две локализованные пластиче­
ские зоны I <х^ <1'\-с и / + c ' < X i < / + с", первая из которьпс примыкает к вершине трещи­
ны, а вторая находится на некотором удалении от первой.
Относительно разностей Aa^j = - и Aw. = щ - щ снова получается упругая задача с
граничтши условиями вида

Аа;^ = 2 7 - 5 ^ ( / < | х , | < / + с,Х2 = 0 ) , Аа^^ = 2 7 - 5 ^ 2 ( / < | x i | < / + c, =0);

Аа22 =^'-522 - а ^ + <|х^| < / + Х2 = 0 ) ; (8)
А^2=0 (/ + c < X i | < / + c', Х2 = 0 ; Ixi > / 4 - с " , Х2 = 0 ) , а,Аф*(Ю = С>(2~^) (z->oo),

где Аа22 = Аа22 - 5 ^ .

Граничные задачи (7), (8) относительно потенциала Аф* сводятся к смешанной задаче
теории аналитических функций для верхней полуплоскости, которая может быть эффективно
разрешена с помощью интегралов типа Коши [4], анализ которьпс в случае дополнительной
пластической зоны представляет, однако, значительные трудности.
Решение граничной задачи (7) может быть получено в замкнутом виде. Для определения
комплексного потенциала Аф* может быть получено следующее уравнение:

dA(p\z) 1
-dt (9)
z-t -I z-t

где аналитическая ветвь корня определяется условием

4z^-{l + cf =z-- (Z->G0).
2z
Интегралы в формуле (9), по существу, аналогичные тем, которые входят в схему Дагдей­
ла, вычисляются явно как это приводится ниже:

'] ^^-dt = -жф^ - (/ + cf - z) ( z g [ - / - c , / + c]); (10)

l+c
-dt = -nxi (x, € [ - / - c , / + c ] ) ; (11)

'Jil + c)'-t' = _ 2 z a r c s i n - ^ + V(/ + c ) ^ - z M n ^ +

t-z I+c z+l

iiUcf-P^il + cf-z^+{l + cf-zl

x ^
\J^^lLjLdt = -2x, a r c s i n - ^ + J(/ + c)'-;c,Mn
J t-x, l+c x^ +1
72
 +4{l + cf-x^^'^l \ 7 \ ' L^^, </ + с); (13)
V(/ + cf -V(/ + cf -xl + (/ + cf -xj

•^(/ + с) ^ ai-csin—^ + Jxf - {l + cf X

t-x^ l+c

. (l + Zf + xJ . (l + cf-x^l] . ^, ^,
arcsm- (14)
m- - - L + arcsin-^^ --—} (\x, >l + c).
(l + xO(l + c) (/-x,)(/ + c)j

В приведенных вьппе формулах yj(l + c f - z ^ обозначает аналитическую ветвь, которая

принимает вещественные положительные значения, если z е (-/ - с, / + с ) , символ In обозна­
чает аналитическую ветвь, которая принимает вещественные значения, если / < z < / + с .
Длина остаточной пластической зоны вычисляется по формуле
— = sec — - ] . (15)
/ 2
Поле остаточных напряжений определяется в виде

5=22
—^ = -1, 1 < . &^с^
I arcsm Р+4
-—^+arcsm- • sec^y^-^
7t
(l + 4)sec^

. sec4/ff/2) + ^
2 arcsm . sec\j3/2)-^
——^^-^+arcsm • (16)
H—
(l + #)sec(/?/2) (l-^)sec(;6/2) j^|>sec>§);

a,2 , 2 . sec'(/?/2) + ^ . sec'(y9/2)-^ r P .\

(17)
= - 1 + —<^ arcsm ——V-^ + arcsm — — \ (sec—^ £\ <. sec B);
¥ 7r[ (l + ^)sec(/9/2) (\-<^)sec0/2)} 2
Остаточные перемещения определяются в виде
COS^(;5/2)
= 21п- + 21n ^ У ^ ^ - ^ I n ^ ^ + ln(^^ -1) +
~2Л COSР [sin(^/ 2) + V l - ^ ^ cos^(;5/2)]^ ^+*

(sinygVl - ^^os^ ;g + ^cos^ yg + l)(sin(;g/ 2)^1 - cos^ iP/2)-^cos^ (;g/ 2) +1)^

+ ^1п
(sinp^l\-^^cos^ p - ^ cos^ у5 + l)(sin(y6/2)^1-^^ cos^(у5/2) + ^cos^ (y^/2) +1)^

(l<^<sec(p/2)); (18)

^ = -ln(^^ -1)-2lncoSy^ + 2ln(sinp + •sjl-^^cos^ P) +

2У7

, a - l)(sin6-v/l - С c o s ' В + £ c o s ' В + 1 )

+ Un^^ PV S P S e_i (sec(p/2)<H<secP). (19)
(^ + l ) ( s i n p V l - cos' p - ^ c o s ' p +1)

Как следует из приведенного решения, вблизи кончика трещины поле остаточных напря­
жений является сжимающим и стремится закрьггь трещину. По мере удаления от трещины,
сжимающие напряжения сменяются растягивающими. Максимум остаточного растягивающего
напряжения достигается у правой грашщы пластической зоны, сформировавшейся непосредст­
венно перед разгрузкой и равен
73
 '22
max = -1 -h —^ arcsin : . + arcsin ——- (20)
2cos'(>5/2)

При малых значениях отношения / F этот максимум достигает значительной величи­

ны, равной одной трети предела текучести.
Безразмерная координата точки, в окрестности которой остаточное сжатие переходит в ос­
таточное растяжение, есть

Решению, таким образом, подлежит краевая задача с условиями (8), в которьпс остаточное
поле напряжений вьшисляется по формулам (16), (17).

2. Локализация пластических деформаций при повторном нагружении.

Решение граничной задачи (8) имеет вид

. / ^ ч t - z \ t - z J, t - z
-(/+с) -/ /+с'

В последней формуле использованы следуюш[ие обозначения:

G(z) = \f - (/ + с")'][z' - (/ + с')'][z' - (/ + с ) ' ] , 27Л(х,) = 7 - 5^2 - ^Гз • (22)

Помимо неизвестршпс параметров локализации пластических деформаций, определению

подлежит также и константа Cj.
Устраняя сингулярность в формуле (21), приходим к условиям на границах пластических
зон:
Ит(2лг7^(^^,А^*) = 0, lim(2;rVG(^^,A^*) = 0, И т ( 2 ; г 7 ^ ^ , Л ^ * ) = 0. (23)
2-^1+с

Введем три безразмерньпс параметра, характеризующих положение очагов пластических

деформаций

А~ / ' А' ~ / + с ' А" /+с

и еще одно безразмерное отношение

(25)
2Y
Условия устранения сингулярностей в граничньпс точках пластических зон примут тогда
форму
1 -2 А "-2
1-S (А'-'-9)(Д"-'-?) ' (А-'Л'-'-9)(А-А''--^)
dq +
ill

А'А'" <?(<?-А"') 2Y1'

-2 A » - 2
1-Е (l-qXA'-'-q) (A-'-qXA-^A'-'-q)
dq- dq +
qiA'-'-q) q(A-'A'-'-q)

74
 1-Z (l-g)(A'-'-^)
i i qiA'-^-q) оо «
g(A-2A-2-^)

(,-A-A-)(,-A--)
A'A'" 9(1-^) 217'

где h{q) = 1 - Б ( (AA" sec(P / 2 ) ) ' > ^ ) ;

sec'(^/2) + (AA'')-'^ sec'(y3/2)-(AA'')-'79

arcsm - +arcsm-
(l + (AA'')-'V?)sec фП) (l-(AA'')-'79^)sec (;ff/2)

( (AA"sec(p/2)f <q<(AA"seep)') ;

1 . 8ес'(р/2) + ( А А " ) " ' л Й " . sec'(p/2)-(AA")"'-J^

arcsin = ~ — + arcsm -
(1 + (AA")-' V^) see(P 12) (1 - (AA")"' 4^) sec(p / 2)

+ — <
. sec'p +
arcsm _
(AA")-'J^. sec'p-(AA")-'V^
+ arcsm -
2л (1 + (AA ")-' V^) seep (1 - (AA ")-' 4i) sec p

(^>(AA"secp)').

Полз^енная система уравнений после вычитания второго уравнения из третьего, а третьего

уравнения из первого приобретает вид (первое уравнение сохраняется неизменным, но стано­
вится последним по порядку следования)
1
(1-1) -dq +
JV 9(А'-' - 9)(А''-' - q) JV 9(А-'А'-' - ^)(А-'А''-' - q)

1 9-A«'
J. Kq).
Д'^Д' b(i-^)(q-A"'A")
dq = 0;

A'-'-, A-'A'-'-q
(1-2) -dq +
Л 9(l-9)(A''-' - 9 ) J\9(A-'-qXA-'A"-'-q)

^-A'-'A"'
KqX dq = 0; (26)
д'^д'-
\q{\-q)(q-A'')

1-Е (A'-'-^)(A'-'-9) (A-'A'-'-9)(A-'A''-'-9)

dq-j.
A^ i l l 9(1-9) 9(A-'-9)
1 (l-,)(,-A"^'A'-')^^^^^^
2 л "2
й(<7).
A'A Д'^Д-
q(q-A"') 2Yl'

75
 Полученная система уравнений (2.6) может быть также представлена в сжатой форме:
(1 - S) J , 3 - J„ + Z,2 = О; (1 - E)J23 - Л . + ^22 = 0;
(1-2)/зз - А ' / з , +A"-'S,2 +(27) Г А С, = О, (27)
или, предполагая, что длина удаленной пластической зоны значительно меньше длины трещи­
ны и заменяя функцию h(q) ее максимальным значением,
(1 - Z) J,3 - J„ + hj,2 = О; (1 - S)/^3 - Л . + ^сЛа = 0;
(1-Е)/зз - А ' / з , + А " - ' й , / з 2 +(2Г)-'Г'А'С, = 0 , (28)
где введены следующие обозначения:

,f2
A-'-q
-2 4 »
-dq;
qaA-'A'-'-qXA-'A'-'-q) ^^ll^d-^X^-A'-'A"')

q-A l-q
^2 = Й(9) -dq; J,3 = -dq ;
Д'^Д'-
qil-qXq-A'-'A"') q(A'-'-q)(A''-'-q)

r f A-'A-'-q . _ V ^-A'-'A"' ,
^^• = i u ( A - ' - , ) ( A - V - ' - / ^ ^ ^^^^.U,(l-,)(,-A-)^^^

q-A'-^A"^
"^22 "~
J
Д'^Д'
Kq) —dq;
^iq(l-q)(q-A"')
J23 = r;;^—
^qil-q)(A"-'-q)

-2 A «-2
(A-'A'-'-^)(A-'A"-'-

0 V 9(A-'-9) o^V 9(1-^)

(l-9)(^-A'"A'-') |(l-,)(,-A''A'-')^^.
^2 = i>2—J-dq;S,2= \ Kq)
q(q-A 9(9-A"')
д«д^
1.1
, ^ 1f . 3cos'(>0/2)-l . sec(;6/2)
=1 - Z < arcsin ^: ^ — + arcsin———-

Для определения константы Q (или безразмерного отношения A^Q / (27/^)) необходи­

мо дополнительное уравнение. Наиболее простой путь получить это дополнительное условие
описывается ниже.
Прежде всего заметим, что из непрерывности нормальньпс перемеш;ений на границах уда­
ленной пластической зоны следуют условия
Aulil + с',0) = О, Aulil + с",0) = 0. (29)
Интегрируя частную производную

ах,
вдоль удаленной пластической зоны и учитьшая формулы (29), получим уравнение

7^^.=0, (30)
1+с' 1
в котором в силу (2) подьштегральное выражение выражается в форме

76
 —т—^(^1,0) = —--\т\дМ (31)
ас, 2G
Подставляя (21) в это последнее уравнение и затем интегрируя (31) вдоль удаленной пла­
стической зоны, находим с помощью формулы (30), что (см. [5], формулы 235.00,235.17)

С, А' А'-'-1
2У7' К(5)

1 А-'Л"-'-?
Il(af,s)dq-
•1 q(A-'-q)iA-'A'-'-q)

A'-'-q (1-,)(,-А-А'-Ъ^^.
-(1-Z)
^qil-q){A'-'-q) n(a\s)dq qiq-A"')

(А'--,)(А"--,_)^^_^д. (A-^A'-^-qXA-^A'-'-q)
-(1-S) dq, (32)
q(i-q) q{A''-q)
где
1-g . „ 2 _ „ 2 A-'-qr q-A'" _ Д-^.д'-г
a =s —-Z ;a, =.y ,"2
• (33)
A"-'-!

Альтернативное представление может бьггь получено с помопц>ю формулы ([6], формула

1.2.27.5)

1 de П(й)',5),
(А"-'-9Х9-А'-^){в-\) в'-в _ 2 _Q>-^^« ^

А"-'-А'-'
где (О =

В результате преобразований вместо (32) получается уравнение (см. также [5], формула
235.00)
С,А'К(з)^^„_2 ' ( 9 - А ' " ) ( 9 - А ' " А ' - Ъ гTI(co^,s)dq
г.„2
Kq). +
2Y1'
qa-q) (2.14)

+ А' * - ' ^ ^ " ^ - ^ > ^ ^ ' ^ ^ " ^ - % ( . ^ . ) . . - ( 1 - Е ) ^'-fX-.'Wo^\s)dq, (34)

qiA-'A'-'-q) q[A -q)

^_ д''-2 _ д ' - 2 1-Д''2д'-2 ^ Д-'Д"-' -Д'-'Д-'

где Й; =
\-д д-'д"-'-^

С учетом уравнения (32) третье уравнение системы (27) примет вид

-(1-Z)J33+J3,+Хз2=0, (35)

или, в случае малого линейного размера удаленной пластической зоны,

-(1-1)/зз+/з,+/гУз2=0> (36)

77
где J^^ = jn(a^,s) -dq; 7зз= U{a2,s\\-:: -dq;
0 q(A-'-q)(A-'A"'-q) ^q(l-q)(A'-'-q)

l-q
A-^ \qiq- A'^Xq-A" 'A'
l-q
Z„ = J h(q)Uiai,s)
q{q-A-'){q-A- 'A'

Константа С, при этом вычисляется по формуле

С,Д' 2 А'-2 - д„л/А-2

i z l X i - A ^ - ' ) , ^ д2 V |(А-'А'-' а - ' АА
" -"-2
^-q)
217'
-А"-' f h{q) .2 '«^9+А'
i 9(А-'-9)

V/(A'-'-g)(A'-'--g)
dq. (37)
9(1-9)

Уравнения (35), (36) заменяют третье уравнение систем (27), (28) соответственно. Из по­
лученных таким образом систем уравнений может быть исключен параметр (1 - S) , § пара­
метр А может быть принят в качестве независимого параметра нагружения.
Итак, решению относительно параметров локализации пластических деформаций А', А"
подлежит система уравнений

Уз, - Л з 2 ^ +^22 Уц+Ап ^^1+A22

(38)
/33 — У32 + "^22 *^13 "^12 / 2 3 + "^22
1 2 A ' -2
=J
д<г2 1
^ - A ' ff ^A
где Л,2 Я(9Х ; -dq; Л„= f H{q\ ^(l-9)(9-A''^)
q{\-q){q-A"'A-') J . '^'"^
1

H{q)Il{a^,s\ 1-9 —dq; H(q) = l-i:-h{q),

Лз2 =
Д"Д"
^(9-A'")(^-A'"A'-^)
ИЛИ, в случае малого линейного размера удаленной пластической зоны, система

(39)
*^33 ~ *^32 •^23 "^2 Лз *^12 *^23 «^2

. 3cos'(p/2)-l . sec(p/2)
где /2 = — arcsin ^^-х-^^— + a r c s m — — -
71 2 cos'(р/2)

Укажем ряд канонических представлений интегралов, с помощью которьпс формулируют­

ся разрешающие системы уравнеьшй через эллиптические интегралы Лежандра:

2А'
A"-'n(Q, с т ' , А ' 5 ) - ( А " - ' - 1 ) F ( Q , A ' 5 )
V A ^ ^

2А'
Л, = {A"-'n(Q, a\A's)-(A"-'-A'-')F(Q,A's)};

78
 n{s\A's)-K{A's) ; (40)

2A'(A' -^-1) 1
2V(1-5^)(1-A'^)

_ 2A4A'-'-l) f , 71 1
J. U(s\A's) + KiA's) + -
VA"-'-I 2^(1-5')(1-A")

-1 A"-'-l
где a =
A" ^ - 1
, sin' a •• A-'A"-'-l

Следует отметить также следуюпще соотношения:

2А'(А'-'-1) 2АЧА' -^-1)

•^22 ^^\2 ~ •1Г(А'5); Л,-Л.=- F(Q.Ms)\
Л/А"-'-1 л/А"-^-1

7гАЧА'-'-1)
•^22 ^^23 - + Лз - i ^ ^ X ^(Г-/)(1-А")'

С помощью канонических представлений (40) и первых двух уравнений системы (28) на­
ходим

1
A"-'n(Q,o-',A'.)-(A''-' - A ' - ' ) F ( Q , A ' 5 ) - ( A ' - ' " l ^ f ^^^_^2)(1_д.2)

(А'-^-1) - П ( . ' , А ' 5 ) + ^

2^(1-.')(1-А")

(A""'-l)F(aA'5)

( А ' - ' - 1 ) - П ( . л' , АА'

' . )^ +^^ ' 1
2V(l-.^)(l-A")

A''-'n(Q,a',A'.)-(A"" -A'-')F(Q,A'.)-(A'-^ f i(i-.^)(l-А^')

(A'-'-l)[-n(.',A-.).K(A-.).fJ^^_^,;^_^,^^

Полученные уравнения затем после ряда преобразований приводятся к следующему виду:

K{A's)h = F{£l,A's); (41)

2д/(1-5')(1-А") {А'-^ - \)F{fl,A's)

К(А'5)-П(5',А'^) +
;г(А'-'-1) K{A's)

79
 +— + A^-'nCQ, cr\A's)-(A"-'-A'-')F(Q,A'.s) (42)
2 1 - A ,2
i \ J

Уравнение (36) на основании (42) преобразуется к виду

д-2д''-2_^ ii'(Q,A'5)'
/П(а^,5). -«9- П(«2,5), :71^—г^^=
19(А-'-9)(А-'А'-'-9) К(А'5) J ' •\(g(l-9)(A'-'-9)

2V(1-^')(1-A") f (A'-'-l)F(aA'5) -I

K(A'5) K(A'5) - П(5', A'i) +1 V(l - )"' (1 - A") +

я-(А'-'-1)

+ А"-'П(ао-',А'^)-(А''-' -A'-')F(fi,A'5)} [ЩагЛ dq-

q){A'-'-q)

l-q
(43)
q(q-A'')(q-A"'A'-')

Таким образом, уравнения (41), (43) образуют систему уравнений для определения пара­
метров локализации пластического течения, пригодную для использования в случае малого
линейного размера удаленной пластической зоны, т. е. А"^А'~^ « 1. Напомним также, что ве­
личина А должна выступать в качестве параметра нагружения.
С помощью уравнения (42) уравнение (34) в случае А"^А'"^ « 1 сводится к следующему:

С,А'Щз)_^2 (А-'-,)(А-'А'-'-,) 2 _mA'.)V {l-q){A'-^-q)^^,

2У7' q{A-'A'-'-q) ' " ' ^ К(А'5) qiA'-'-q) n(co\s)dq-

2V(1-^')(1-A")f (A'-'-l)F(Q,A'5)
K(A 5) -U(s\A's) + ^ V(l-^'r'(l-A")"' +
л-(А'-'-1) K(A's)

+ А''-'ЩП,с7\A's) - (A"-' - A'-' )F(Cl,A's)} М ^ П ( . ' , . ) ф .

qiA'-'-q)

it-2 b-A-'lb-A-'A-']^^ a Ь
-A (44)
A'^A' I'll-?)

Численный анализ локализации пластических деформаций проводился в случае малого

линейного размера удаленной пластической зоны. Система (39) была заменена эквивалентной
системой уравнений

(45)
«^22 «^23 h\ " ^22

23 (46)
•^32 ~ Лз /2,-^22
80
в которой

(47)
IqiA'-qXA "^Д' -q)

Численное решение системы уравнений (45) и (46) относительно параметров локализации

А' и А" сопряжено со значительными трудностями и требует значительньпс затрат времени.
Для оптимизации счета бьш специально разработан алгоритм ускореьшого вьршсления входя­
щих в уравнения квадратур, основанный на разложеьши эллиптических интегралов в ускоренно
сходящиеся ряды Фурье. Численные расчеты, которые проводились для сл>^ая, когда линей­
ный размер удаленной локализованной пластической зоны значительно меньше длины трещи­
ны, подтверждают разрешимость указанной системы и, таким образом, также обосновывают
предлагаемую схему локализации пластических деформаций.
Если характерный линейный размер удаленной пластической зоны нельзя считать малым
по сравнению с длиной трещины, то система уравнений (45) и (46) должна быть заменена на

•^22 ^^23 ^21 ~ ' ^22 ^^32 "^33 ^21 ~ ^22

где А 2 1 = j H{q)
А'^А-^ V

3. Вычисление раскрытия трещины при двухзонной

локализации пластических деформаций.

Раскрытие трещиш»! — одга из важнейших параметров механики разрушения. Поэтому

оценка влияния удаленной пластической зоны на трещину может быть вьшолнена на основе
вьршсления изменения раскрытия трещины, обусловленного влиянием второго очага локали­
зации пластических деформаций. В принципе, для точного определения раскрытия трещины
необходимо дополнительно учесть изменение геометрии берегов трещины, искаженньпс оста­
точными перемещениями. Однако перенос граничньпс условий на деформированные берега
трещины делает невозможным аналитический анализ проблемы. В дальнейшем предполагает­
ся, что граничные условия формулируются на неискаженньпс берегах трещины.
Рассмотрим уравнение

^(х„0) = '^1т[^А<Р*1 (48)

интегрируя которое вдоль пластической зоны, примыкающей к кончику трещины, получим

. 1+с
-Au2(h0) = ^^ {bn[3zA(p,]dxi . (49)
2G /
Раскрытие трещины поэтому может быть вьшислено по формуле

— f = 1п ^^--тр; m[dzA(p,]dxx. (50)

4// cos/? уI J

Второе слагаемое в правой части полученного уравнения может быть вьгчислено только
через три параметра локализации пластических деформаций А, А', А" , как показано ниже.
Используя формулу (21), получим

л- de

I JV(A"-'-^)(A'-^-^)(l-^)

81
 А' de

(А'-^-^')(А''^-^')(1-^') d9-

в'
(A'-^-^')(A'-^-g')(l-^'),„
•

Вводя обозначение
(51)
можно заключить, что
Im[^,A^*] ^, . _ L . F 0 ' , ^ ) -
ГА/ 2}7^

\-q
[(А'-2 - 1)П(г, , 5) + (1 - 9)F(>',5)]cf^-
q{q-A"^Xq-A'^A-^)

A-V-2-9
1|д(А-2-д)(Д-2д'-2-9)
[(А'-2 - 1)П(г, , S) + (1 -q)F(y,s)]dq -

-(1-2) [(А'-^ -l)П(^y5^.) + (l-q-)F(r,5)Wg, (52)

q(l-q)(A-^-q)

q-\ q-ix q-/\

Подставляя формулу (52) в (50) и выражение для ( 1 - 2 ) , согласно уравнеьшю (42), пред­
полагая тем самым малым размер удаленной пластической зоны, после многочисленньпс пре­
образований получим

4AYI cosp 2Yl^

A-'A'-'-q
[(A'-2 - , 5 ) + (1 - q)Fiy, sWq-
q{A-^-q\A-^A-^-q)

A"' 9 [i,A'-^-\)n{y,p\s) + {\-q)F{y,s)\dq~

K(A'.) и q(l-q)(A'-'-q)

2У(1-^Ъ(1-А'Ъ [ (A^-^-i)F(Q,A'^) к(A'.)-П(.^A'.)+|д/a-^'гЧl-A'Ъ"'

л-(А'-^ -1) I K(A's)

+ А'-^ЩД , A'^) - (A"-2 - A'-2 )F(Q, A'^). X

82
 1
Д -9 гМ'-2

l-q
[(A'-2 -m(r,pls) + (l-q)Fir,s)]dq\. (53)
^,}lq{q-A'')iq-A"'A'-')

Если исключить теперь константу С, из последнего уравнения с помощью соотношения

(44), то для вычисления раскрытия трещины при повторном нагружении можно получить фор­
мулу
У А ^ = А - ' > / А ^ 1 П ^ ^ ^ ' Ч ^ ^ ) Н-
4AI7 cosp

F(r,s) (А-^-^)(А-^А'-^-^)^, , F4r,s)

II{co\s)dq-
ВД ill ^(А-^А'-^-^) qiA-'-q)

2F{r,s) У(1-^Ъ(1-А'Ъ (A'-^-l)F(Q,A'5)

K(A'.) - П ( . ^ A'.)+|У(1-,2)-Ч1-А'Ъ-'
яА^Щз) (A'-^-l) K(A's)

( l - g ) ( A ' - ^ - g ) u{(o^,s)dq-
„/ 2
+A''-^n(Q, a\Ms)- (A"-2 - A'-^ ) F ( a A'5)
1^ 9(A"~'-9)

- A - J J
9(1-9)

1
A-^A"-^-q
-J \q(A-^-q)(A'^A'-^-q)
[(A'-2 - 1)П(г, J3^ ,s) + {l-q)F(r, s)]dq-

Fin, A's ) ' ^"^ ^ [{A'-^-\)n(r,P\s) + il-q)F(r,s)]dq-

K(A's) i^qd-qXA'-^-q)

2д/(1-^Ъ(1-А'Ъ (A'-2-i)F(0,A5)
K(A'5)
K(A'^) - П(5^ A'i)+1 Va - (1 - А'Ъ"'
;г(А'-^ -1)

+А''-2п(П, <т\ A'^) - (A'-^ - A'-^ )F(Q,i,A'.)L

A

[(A'-2 - 1)П(г, у9 ', 5) + (1 - q)F(r, s)]dq +

\qil-q)iA"^-q)

83
 [(A'-' - 1)П(г,Pi,s) + {\-q)F(r,s)]dq . (54)

Уравнение (54) позволяет вычислить раскрытие трещины, если три безразмерньпс пара­
метра локализации пластических деформаций А, А', А" определены или заданы. Тем самым
может бьп^ь охарактеризовано влияние удаленной пластической зоны на трещину.
В случае, если характерный линейный размер удаленной пластической зоны сравним с
длиной трещины, то на основании формул (50) и (52) и (37) находим

лг £ W2 (hO) ^ cos^(y?/2) ^ А
.F(r,s)-V,+V2-(l-i:W + V2)\, (55)
4Y1 2W
cos

J\gil-q)iA'-' -q)

q(A-' -q){A-'A'-' -q)

l-q
[(A'-2 - 1)Щг, Pi,5) + (1 - q)F{r, s)]dq ;
q{q-A'^){q-A"^A'-^)

i - = A''~^A-^S,2 + In - ( 1 - 2 ) А - ^ ( / з з + А"-^/з2);
lYl

(l-q)(q-A''A'-')^_ , v_-^11+^12
532 = H{q\ '-dq; 1-E =
qiq-A-") J,3 +7)2

Таким образом, если параметры локализации определены, то формула (55) позволяет вы­
числить раскрытие трещины при любой протяженности удаленной пластической зоны (для
этого необходимо вычислить только приведенные вьппе квадратуры).
Асимптотический анализ величины раскрытия трещины при малой (по сравнению с дли­
ной трещины) длине удаленной пластической зоны приводит к простой зависимости

Au2 = - г [ 2 - (4Г)-'яд/(А'-'-1)-'аг21] ( А'^А'-^ -> 1 - О ) (56)

417

для дополнительного закрытия трещины при повторном нагружении.

Таким образом, возникновение нового очага пластического течения при повторном нагру­
жении проявляется в форме дополнительного закрытия трещины, асимптотическое значение
которого может быть вьпшслено по формуле (56).

84
 БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Dugdale D. S. Yielding of steel sheets containing slits // J. Mech. Phys. SoHds. 1960. V. 8 . No. 2. P. 100-104.
2. Качанов Л. М. Основы механики разрушения. М.: Наука, 1974. 312 с.
3. Мусхелишвили Н. И. Некоторые основные задачи математической теории упругости. М.: Изд-во АН СССР, 1954.
648 с.
4. Штаерман И. Я. Контактная задача теории упругости. М.- Л.: Гостехтеориздат, 1949. 272 с.
5. Byrd Р. F., Friedman М. D. Handbook of elliptic integrals for engineers and physicist. Berlin-Gottingen-Heidelberg:
Springer-Verlag, 1954. P. 68-161.
6. Прудников A. n, Брычков Ю. A., Маричев О. И. Интегралы и ряды. М.: Наука, 1981. 800 с.

85