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

книги / Механика трещин

..pdf
Скачиваний:
9
Добавлен:
12.11.2023
Размер:
11.85 Mб
Скачать

£ = —=- ((n2 + Q2)I LF~ liqn^1*),

(1.24)

ЦК

 

 

 

R - (n2 + я2)2 ~ Ц 2п1п2 .

 

 

Если в выражении

для

R положить s = iqc, то

уравнение R = О

при q ФО

 

 

 

Я*(с)= (2 - Ь2с2)2 -

4 л/l -

Ь\с2 УГ - Ь2с2 = 0

(1.25)

определит скорость поверхностных волн - волн

Рэлея: с = сд > 0.

Скорость волн Рэлея близка к с к о рости во л н сдвига (сд < с2). В частно­ сти, при V = 1/4 сд = с2 у2 - 2/>Д « 0,9194с2.

Исследование рассматриваемой задачи - задачи Лэмба - приведе­ но в [93]. Здесь мы ограничимся определением перемещений на гра­ нице полупространства х 2 = 0. На основании соотношений (1.23), (1.24) для х 2= 0 получаем соотношения аналогичные (1.20), но для плоской задачи

(1.26)

Для функции Sf'*' в обозначениях, принятых в формуле обраще­ ния (1.16),

(1.27)

гг»! = л/bf + р2,

т 2= л/ь| + р2;

§F(s) = 1/s, g' (f) = ô(0,

tv= tv0 = 0.

Радикалы, входящие в равенства (1.27), и функция R0 определя­ ются так:

(Ix l^ c^ ,

y= 0),

(W < c,< ,

y = ± 0 );

f o i - t2/x2

(\x\>c2t,

y=0),

 

m2

 

 

 

 

 

± i y/t2/x2-

(Ixl < c2t,

y= ±0);

R o= (b% - 2t2/x2)2 +

 

 

 

 

 

4(t2/x2) y/b\- t2/x2 y/b2t2/x2

(\x\ > Clt,

y = 0),

+< ± 4i(t2/x2) y/t2lx2-

b2 y/bl~

t2/x2 sgnx (c2f < Ixl < Cjf, y= ± 0),

- 4(t2/x2) y/t2/x2-

b2 yjt2!x 2-

(Ixl < c2t,

y 0) (x= x 1 Ф0).

Отсюда и из соотношений (1.16), (1.17), (1.27) находим = x v c=x/t)

5 i 2= - / / й= - / 0/^ = 0

(Ixl >с^),

 

 

 

_

2bfc2

( 2 -

bfc2)

V 1 - b jc2

yj b\c2- 1 sgnx

 

12

щГГ”

( 2 -

b2c2)4+ 16(1-

b2c2)(b2c2-

1)

 

(c2t <

Ixl ^Cjl),

 

 

 

 

 

 

s i2=

0 ^ 2 -

b2c2)

*(<#“

lxl)sgnx

(Ixl < c2f),

 

 

* 2p, / У

 

 

 

2\ic2R*(CR)

 

 

 

 

 

 

*'*(<*) “ Vo

4CP

 

M b? + b2 - 2d tfb ÿ - b2(2 -

b2c2) 3].

(1.28)

242

 

( 2 -

O2Cr )

 

 

 

 

 

 

A

При выводе последнего равенства учтено, что функция Д2) = //(2n£j является аналитическим представлением ô-функции Дирака 6(х).

Аналогичным образом определяются и остальные функции:

5 ц -

S22S21

0

 

(Ixl ^ Cj/),

 

 

 

 

4b2c2(l -

 

bfc2) Vl -

bjc2

 

 

5“ = л ц ф -

b|c2) 4+ 16(1 - Ь2с*)(Ь2с2-

1)]

 

 

b2c\2 -

 

b^c2) 2 Vl -

b?c2

 

 

Sz2=

Л(Х4(2- bfc^4+ 16(1 - bicubic2-

1)] ’

S2 l~ - S12

 

(c2t < Ixl < q # ;

 

 

 

b\c2 y/l

-

b\c2

 

_ bjc2 Vl -

bfc2

 

Л|11й*(с)

5з2

лц!й *(с)

$21

^12

 

(0 < Ixl < c2f);

 

 

(1.29)

 

 

 

- 182 -

 

 

 

 

 

 

 

 

 

*$ ц - $ 2 2

Л|it

'

*$12

*$21 ®

( *

 

 

 

 

 

(с = x/f,

R^c) ~ -

 

2Ь|с2(1 -

cllcl) = -

Ь|с2/(1 - v) (с - 0)).

Функции Sll9. . . ,S 22>как видно из первых двух равенств (1.26), определяют перемещения при сосредоточенном импульсном воздей­ ствии, например u1=Sll (т = ô(f)ô(x), T L F =1, о = 0). В общем случае решение выражается двойной сверткой:

u i = * $ и * * 1 + *$12 * * 0 >

U2 = 5 2 1 * * T + 5 2 2 * * ° -

Характерная особенность приведенных решений - изменение зна­ ка перемещений, когда с становится больше скорости волн Рэлея. В области Ixl < cRt соответствующая компонента перемещения грани­ цы полупространства направлена в ту же сторону, что и действующий на границу импульс. Действительно, если oLF= l, т = 0, то импульс действует вдоль отрицательного направления оси х 2 и перемещение u2= S 22 если же 0 = Û, TLF = 1, то аналогичное заключение можно сделать относительно импульса напряжений т и компоненты переме­ щения их = Slx <0. При тех же условиях, но в области 1x1 > сд*, указан­ ные компоненты положительны (кроме тех точек на оси с= x/f,. где они обращаются в ноль). Графики функций Я*(с) (кривая 1) и а(\>) SxJ[a(v) = = - v)] при с Ф CR (кривая 2) показаны на рис. 5.2, а графики а(v)5n , Û( V)5 22, а(0), S33на рис. 5.3 (S33перемещение в задаче III (1.18) для х 2=0). В расчетах принято v=0,3 ( с ^ ! ^0,5345, сд/с2 ^ % 0,9274).

Для определения напряжений на продолжении трещины и переме­ щений ее берегов достаточно знать функции Smn на оси х= х г Однако выражения для них (1.28), (1.29) слишком сложны, что затрудняет решение конкретных задач. В связи с этим предложено приближенное

описание функций S119 S22, сохраняющее основные черты точного, но существенно более простое [86, 96,148].

В приближенном описании функции

 

заменяются на

cLF

1 -

V

 

Vb^s2 + q2

 

 

 

 

ю ~

И

 

 

s2/c2R + q2

1

 

 

 

 

 

 

 

 

 

 

 

1 -

V

 

Jbli'S2 + g2

 

 

 

(1.30)

cLF = -

Ц

 

s2/c| + q 2

5

 

 

 

°2 2 0

 

 

 

 

 

При

этом

сохраняются

полюса s = ± iqcR,

отвечающие волнам

Рэлея, и асимптотика при q

°°:

 

 

 

 

сLF ъ cLF

ъ cLF^ çf.F

 

1 - V

y/q+ r'O

-

 

ц

q2

 

 

°11

°110

 

°22

°220

 

 

 

R ~ 2q2s2(b| -

Ь?) =

b|q2s2

( q - 00)

,

 

 

 

 

 

1 -

V

 

 

 

 

а значения постоянных

можно взять либо равными Ь^2, либо

определить их из условия' лучшей апроксимации в том или ином

смысле [86]. В частности, можно потребовать,

чтобы функции

и Sjjfffc (т = п) асимптотически совпадали в

окрестности полюсов,

отвечающих волнам Рэлея, т. е. при s

± iqc^ . При этом Smn ~ SmnQпри

с = x/t -►± Сд .

 

 

 

Формула обращения (1.16) или (1.17) приводит к функциям

*^110= ® ^ ^

*^220= 0 ( ^

^

 

1 -

V

\/l - b|*c2

S ll0 ‘ W

 

c V c t-

(Ixl < c2J=t/b2J ;

 

1

1 -

V

л/l - bf*C2

5220 ■ W

 

с2I c i -

(Ixl ^ cl1tt = f/b1#).

 

1

Отношения S22Q/S22,

Sll0/Slx (кривые 1, 2 соответственно) для

v = 0,3 показаны на рис. 5.4, где принято указанное выше требование асимптотического совпадения, которое определяет следующие значе­ ния постоянных:

1

(

( 2 - ^ с » - ( 1 - ь ; ас Щ | - б У

) и ,

CR

I

" Н-

2(ь ;ьь |;с 0

-

Ь}(2 - Ь’ с ® » ]’

)

Зависимость

отношений

2;*/Ь2

от

коэффициента Пуассона

показана на рис. 5.5 (кривые 1, 2 соответственно).

§5.2. Состояние у края трещины

ипоток энергии

Прежде чем переходить к решению конкретных задач о динамике трещины, найдем распределение напряжений и перемещений у ее края, а тем самым - соотношение между интенсивностью напряжений и по­ током энергии в край трещины в зависимости от ее скорости. Как пока­ зано в статье [121] и как это будет видно ниже, распределение напря­ жений и перемещений у края трещины зависит лишь от ее текущей скорости. Исключением являются те моменты времени, в которые к краю трещины приходят разрывные фронты волн напряжений от внешних сосредоточенных источников. Учитывая сказанное, состоя­ ние у края трещины можно определять, исследуя стационарную задачу

l(t) = ut,

u = u(x1- o t ,

х 2), о = const

> 0,

(2.1)

где х г = /(f) -

координата

правого края

трещины,

расположенной

•на осих 1(х1 < /, х 2 = 0).

 

 

 

При этом для того чтобы распространить результаты на нестацио­ нарную задачу, достаточно значение постоянной скорости трещины и принять равным значению скорости l{t) = dl/dt в нестационарной задаче в данный момент времени, т.е. положить и = и {t) = l(t). Подчеркнем, однако, что здесь речь идет лишь о распределении напряжений и пере­ мещений, а не об их амплитудах. Для определения последних, напри­ мер коэффициента интенсивности напряжений, необходимо рассматри­ вать конкретную задачу, и если эта задача нестационарна, то решение будет зависеть не только от скорости роста трещины в данный момент времени, но и от всей предыстории ее движения.

Обозначив х г - ut= х и полагая, что того же типа зависимость (2.1) имеет место и для потенциалов ф = ф(х, х 2), ф = ф(х, х 2), из уравнений (1.4) получим

а 2Ф

д2ф

а= 1 -------;

а -------+ --------= О,

дх2

дХо

о ’

 

Р

а2Ф

а2ф

т

т

Р = 1“

- •

(2.2)

дх2т

тдх\+

=

° .

 

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

неравенствам 0 < и < с2. Тогда а, р > 0 и можно положить

 

 

\А*~= P,

х 2 \/р = g,

 

(2.3)

в результате чего уравнения (2.2) преобразуются к виду

 

а2(р

а2«р

 

 

аф

а2Ф

(2.4)

+ — = о,

 

 

-----+ ---------- 0.

дх2

др2

 

 

дх2

dq2

 

Решения этих уравнений - гармонические функции переменных х, р(для ф) и х, g (для ф) - представим в виде

f=Re(Azp,

ф= Im(Bz3,

(2.5)

zp= х + ip, zq= х+

/g; А= а + ib, В= с+ id= const

 

Для антиплоской деформации, как следует из уравнения (1.7),

можно положить

 

 

u3= Im(Cz£-‘);

с= е+ '/= const

(2.6)

Дальнейшая конкретизация выписанных выражений проводится, как и в случае статики, из условий на берегах трещины (отсутствие внешних напряжений) и на ее продолжении: симметрично относитель­ но оси х х перемещение U j- задача I, перемещение и2задача II, антисимметрично перемещение и3задача IIL Кроме того, учитываем условие ограниченности и отличия от нуля потока энергии в край трещины. Из последнего условия для трещины, берега которой свобод­ ны от внешних напряжений, следует s= 3/2.

Из условий

симметрии

находим,

что для

задачи I b= d= 0, для

задачи II а= с= 0 и для задачи III /= 0.

 

 

Подчиняя решения условию отсутствия напряжений на берегах

трещины, приходим к уравнениям (задача 1):

а2Ф\

1дих

ди2\

I а2ф

а2Ф

°i2= I1 -------+ --------1= Д

2 ----------- + ------------------ =

\дх2 дх J

\ дх2дх дх\ дх2)

.

д2ф

д2ф

д2ф

 

 

 

 

= Д 2 ^ — — +Р

дх2

 

 

 

 

дхдр

dq2

 

 

 

3

 

6

 

 

 

6

= 0

2 т/аГр1'2 a sin -^

+ г"1/2 с(1 + P)sin

V

 

 

 

 

 

 

 

(0р =е„=л);

 

 

 

 

 

 

=2ц

ди 2

, ди1

ди2\

II -

р

ди2

■+ X

ÔXj

= Д

а

дх,

 

дх,

ÔXj

1 -

2 а - р - 1

d u , \

/1 - р

д 2ф

+

2 а - P - 1

+ --------------------

а

—— = д [ -----------

а ----------

 

х

1 -

дх

1 - а

др2

 

1

- а

д2ф

<— д2ф \

 

1 +,— Р■а cos —

+

х ------------

2 V P -------------

1= ---------

д

дх'2

 

dxdq J

4

ФЪ

 

2

 

Р

с cos

0,

 

 

 

 

 

+ 2 / —

0

 

( 0 р - 0q

л);

 

Гр=у[-X2+ р2 ,

rq = y/x2 + q2;

 

 

 

0р = arctg (р/х),

0q = arctg(q/x).

 

 

 

Из первого уравнения находим с = - 2 \fâal{l + Р). Второе уравнение удовлетворяется автоматически (вследствие принятого значения s). Определяя тем же путем напряжение о1Х и компоненты перемещения, получаем следующие зависимости.

Задача I:

 

 

 

cos (0р/2)

4 /а ]Г

cos (0д/2)

О ц = — да (1 + 2 а - Р)

1 + P

JrZ

4

 

 

 

 

- (1 + p)

COS (0p/2) 4 Va P

cos (0q/2)

да

V "ô T

 

1 + P

<>i2 = y N

i / «

‘ sin (0p/2) sin (0q/2)

(2.7)

 

 

 

 

 

,—

 

2 V5T

/—

 

 

V Гр cos (0p/2)----------—

v -rq cos (0д/2)

 

 

 

1 + Р

 

 

u2 = — а

-

у[7рsin (0р/2) + ----- — / г Г sin (0д/2)

 

 

 

1 + Р

 

°и = у

ц</VF

 

Sin(0p/2)

SinfGgÆ)

 

 

 

 

1 + p

\/"rp

\/rq

°22 = y

H \ /F

sin(0p/2^ sin(0y2)‘

 

3

 

 

 

T ^ T

 

 

4 V«P

cos(0p/2)

 

î W

o ., = — pd

1 + p

--------- ( I . t l ü

 

 

v fp

 

Vr*

 

 

 

 

 

 

 

2

 

 

 

«1 = —

л/F ^ р -\ Л > 5 т (0 р /2 )-

лЛ~зт(0,/2)

 

2 л/аЭ~

г—

 

 

u, = — d 7 7 J -

V гр cos(0^ 2) -

cos(0g/2)

(b = - 2d д/Р/(1 + P));

задача Ш:

1

sin(0g/2)

COS(0q/2)

° 13=' T

e|1_7 ^ ~

°23 = ~r en A/P

V Гр

u3= e

sin(0q/2).

 

(2.8)

(2.9)

Для трещины конечной длины эти формулы определяют асимпто­ тику перемещений и напряжений у ее правого края.

Если в полученных зависимостях положить a, d= const о " 2и устре­ мить скорость и к нулю, получим в пределе аналогичные соотношения

для неподвижной трещины (2.2.19) -

(2.2.21).

 

 

Заметим,

что

на

оси

х

гр= rq= л;

0р= 6<j= 0 (х >0,

х2 = 0),

rp= rq= - х>

0Р= Йд= ± л (х <0,

х2 = ±0).

Учитывая это и

вводя

новые постоянные -

коэффициенты интенсивности напряжений (2.2.16),

из формул (2.7) - (2.9) находим следующие зависимости:

 

задача I

 

 

 

 

 

 

 

 

 

о22= К\/ ^2лх= ity у/х

( х >0,

х 2 = 0),

 

 

 

2

х

л/0<1

"

 

Р)

 

( . )

 

 

-----

----------

=

± Mi у - х

 

± к 1. ~

 

 

1»К*(и)

 

 

 

 

2 10

 

 

 

 

 

 

 

(х < 0,

х2 = ± 0),

 

 

 

 

 

 

 

 

 

 

 

где функция i?*(u) определена равенством (1.25);

012 ^il/ л/2л~x~-N\\!yfx

(х> О, х2 = 0);

(2. 11)

< О, х2 = ± О);

задача III:

°22 к\\\1/ 2лх = ЛГШ/ /х~ (х > 0, х2 = 0);

(2.12)

(х < О, х2 = ± 0).

Здесь, как и выше, в пределе при и 0 получаются соотношения для статики (2.2.15), (2.2.16).

Покажем, что формула (1.3.4), выведенная для квазистатики, справедлива и для динамики [104]. Пусть вектор напряжения, дей­ ствующий на продолжение нижнего берега трещины, 0+ ( х - of), пере­ мещение верхнего берега и. (х - иf), перемещение нижнего отличается

лишь

знаком. Приложим к нижнему берегу напряжение ао+ (х -

uf + т)

и тот

же

вектор, но другого знака, - к верхнему берегу

(т > 0,

0 ^ а

^ 1,

0+ (х) = 0 при х < 0, и_ (х) = 0 при х > 0). Ввиду линейности

задачи дополнительное перемещение (на интервале - i ^ x - o f ^ O ) иа — аи7 (т 0). Сингулярность в точке х = of, порождающая поток энергии, приобретает множитель 1 - а, а в точке х - of = - т появляет­ ся та же сингулярность, но с множителем а. Потоки энергии через указанные точки пропорциональны квадратам соответствующих множителей, а поток энергии (на единицу приращения длины), созда­ ваемый нагрузкой, приложенной на интервале - т < х - uf ^ 0 к бере­ гам трещины, равен

(1 - а) lim [и_ ( - т) *о+ (т)]

Заметим теперь, что напряжения ±о+, приложенные к бере­ гам трещины на малом отрезке у ее края, не влияют на общий поток энергии (при т -*• 0 их влияние на фиксированном расстоянии от края трещины исчезает), они лишь перераспределяют его из

точки х - 0 на отрезок - т < х ^ 0. Поэтому справедливо уравнение

Г = (1 - а)2Г+ а2Г+ 2а(1 - a) lim[u_(- т) * о+(т)]',

из которого и следует формула (1.3.4).

Влияние скорости проявляется лишь в том, что отношения Afj/Wj, Mn/Nu, Miu/Niii, как видно из формул (2.10) - (2.12), зависят от скорости.

Итак, формулы (1.3.4), (2.10) - (2.12) приводят к следующему соотношению для потока энергии в край трещины, распространяющей­ ся в условиях плоской и антиплоской деформаций упругого тела

[36,121]:

 

 

(1 - Р)(7«

VFКц) + Kill

(2.13)

а

д

 

где а = 1 - и2/с2, р = 1 - о 2/с2,

a функция Я*(и) определена равен­

ством (1.25).

 

 

При переходе к плоскому напряженному состоянию, как уже отме­ чалось в §5.1, следует изменить значение коэффициента Пуассона - заменить параметр с2 на с\ (1.6). Учитывая неравенства Д*(и) < 0 (0 < о < сд), Я*(и) > 0 (сд < о ^ с2), видим, что перемещение верхнего берега трещины в задачах I, II того же знака, что и напряжение на ее продолжении

U2°22 > 0, Ul°12 >0, (2.14)

если скорость трещины меньше скорости волн Рэлея. При этом поток энергии в край трещины Г > 0 независимо от знаков указанных компо­ нент (2.13). В противном случае (сд < и < с2) неравенства (2.14) меняют­ ся на обратные, т. е. если в задаче I трещина раскрывается, то на ее продолжении напряжения о22 должны быть сжимающими. Растягиваю­ щим напряжениям соответствует „перехлест” берегов трещины (допу­ стимый лишь при такой постановке задачи, когда на перемещения берегов не накладывается ограничений, что, конечно, неприемлемо). При этом независимо от знака компонент о22, о 12 > 0, х 2= 0) поток

энергии Г < 0

(Т-+ 0 при о~+с2). Последнее

неравенство

приводит

к следующему

выводу: если эффективная

поверхностная

энергия

положительна, то трещина не может распространяться со скоростью, превышающей скорость волн Рэлея (сд < и < с2).

Из формулы (2.13) вытекает еще одно важное заключение: в пло­ ской задаче при фиксированном коэффициенте интенсивности напря­ жений поток энергии в край трещины неограниченно растет с увели­ чением ее скорости (0 < о < сд). Так, если в статике

1 -

V

ton Г = Г0 =

(2.15)