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

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

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

ТО П р и и -*■ Ср

1

4

\Л -

сЦс\ 4 + y/l - сЦс\ K\i

 

(2.16)

Г(и)--------—

c i

 

 

—У оо

2 ц

R*(CR)(CR - и)

 

 

 

 

 

 

 

Для антиплоской задачи критической является скорость волн

сдвига

 

 

 

 

 

 

 

1

Кш ;

1

 

 

lim Г =

 

г(и)

 

С2).

и-»о

 

 

¥

л/2(1 - и/с.)

(2.17)

Отношения -

 

Т(и)/Г0

для о < с2

совпадают с

функциями

a(v)S22

(задача I), a(v)51A

(задача II), а(0)533 (задача III)-

см. формулы (1.18),

(1.29), (2.13) и рис. 5.3.

Таким образом, в динамике особенно ярко проявляется некор­ ректность рассматриваемой континуальной модели сплошной упру­ гой среды. Оставаясь в рамках этой модели, мы получаем вполне определенную связь между силовыми и энергетическими параметра­ ми механизма разрушения (2.13) - (2.17), причем зависимость этой связи от скорости оказывается одинаковой для любого материала. Но в теории упругости не содержится сведений о механизме разруше­ ния. Этот парадокс, как уже отмечалось, нельзя разрешить, не выходя за рамки континуальной модели без внутренней структуры.

Естественно принять силовой критерий разрушения, т. е. полагать, что разрушение упругого материала происходит при некотором значе­ нии силы, действующей на элемент определенной (малой) площади. В данном случае это эквивалентно введению критического значения коэффициента интенсивности напряжений. Предположим, что поверх­ ностная энергия также фиксирована. Тогда она не может превышать того значения, которое получено выше (2.15), (2.17) для о = 0 (иначе критическая сила не приводила бы к разрушению, т. е. не была бы кри­ тической). При этом избыток энергии Т(о) - 2у, неограниченно увели­ чивающийся с ростом скорости трещины, должен уноситься упругими волнами высокой частоты. Эти волны, после того как их интенсивность станет достаточной, могут приводить к дополнительным поврежде­ ниям материала у берегов трещины. Возможно, именно этим объясня­ ется явление, наблюдаемое в опытах на некоторых материалах: вна­ чале, пока скорость трещины мала, ее берега оказываются гладкими, а после того, как скорость становится достаточно большой, - шерохо­ ватыми (см. рис. 5.1). При большей скорости (около половины скорости волн сдвига) может наблюдаться ветвление трещины. Этот факт обыч­ но связывают с тем, что при такой скорости направление, на котором растягивающие напряжения максимальны, не совпадает с продолже­ нием трещины, а составляет с ней некоторый угол [79].

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

Считая, что решение рассматриваемой „однородной” стационарной задачи существует как предел решения некоторой нестационарной задачи, можно перейти от соотношений (1.26) к соответствующим зави­ симостям для стационарной задачи, положив 5 = 5' + iqu, домножив обе части равенства (1.26) на s' и устремив 5' к нулю [93,115]. В результате получим

UFI *(Q) = P M

0Fil{q),

u^*(q) = P2{q)o[* (q)

ит*(я) = S *♦0

 

+ »«и, q);

 

=

Gmn

(2.18)

 

s'-* 0

 

 

P

(q) = liinSff (5r + iqu, q),

9

S'* 0

9

 

где символом F* обозначено нат, движущейся вдоль оси = их{х\ - и0, то

преобразование Фурье в системе коорди­ x t со скоростью и. Например, если и1 =

if*

= 5

-

о

 

 

-СО

 

 

 

Для приближенной модели (1.30)

 

1

_

1

ц(1 - i>7cj?) д2

iO)- 1/ 2

 

 

 

r(q + iO)_1/2( q -

^1,2

^10,20

(1 ' v ) V l - bf2)1>u:

(2.19)

Положим

 

 

 

 

 

u ^ ^ A ^ y/ ü T ^ xlH io t-x t),

 

F

=

yfn

i-3l2( q - m - 3'2.

(2.20)

 

 

Обращаясь к формулам (1.26) и учитывая найденные выше зависимо­ сти (2.18), (2.19), получаем

F*

1

F*

 

 

 

 

<>12,22=-;------- « *

=

 

 

 

 

 

г 10,20

 

 

 

 

 

1-1

 

1 -

о 2/CR

t

Viïi-1/2(q + iO)-1/2.

(2.21)

2(1 -

V)

Vl -

bf2,i> о 2

1,2

 

 

 

Отсюда, вводя коэффициенты интенсивности напряжений, находим

°12,22

Ifu,l/ )/2лх

(x = x t -

uf >0,

*2 = °)>

 

 

 

1 - V

Vl ~ ь(2,р* ц2 f

2x

 

«1.2 = ± Кид---------

1 - u2/ 4

v

---- (x < 0,

 

 

 

 

Л

 

^12,22 ^ (*^ ^

^)>

~ ^ (x ^ О, X2 "

0),

 

1 -

V

 

bl* u2).

(2.22)

Г=7\

(Kf Vl - b 2, u 2 + Kft ^1 -

2|i(l -

ua/cfc)

 

 

 

Отношение потока энергии Г* , вычисленного по формулам (1.31), (2.22), к его точному значению Т (2.13) при одинаковых коэффициен­ тах интенсивности напряжений совпадает с отношением S220/S2g (задача О или 5 .^ /5 ^ (задача II). Графики указанных отношении показаны на рис. $.4.

Рассмотрим еще задачу I для распространения разреза с межзву­ ковой скоростью [87,114]: с2 < u < Cj. При этом р = - у < 0 и уравнение относительно ф (2.2) становится волновым

0

0

(2.23)

дх\

= 0;

^ дх2

 

так что дозвуковое решение (2.3) -

(2.6) не подходит.

С учетом симметрии можем ограничиться решением для верхней полуплоскости. Положим, удовлетворяя уравнениям (2.2), (2.23),

(p = ARe(x+ îV ây)VK,

ф = В ( - х -

y/yyŸKH (- х - у ),

(2.24)

где А, В-вещественные

постоянные;

vK= const Определяя

отсюда

напряжения и подчиняя их граничным условиям

 

°12=О22=0> * < 0 ,

У“ 0,

 

 

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

у -

1

 

 

; к + со,

1

(У - 1 )а

;

 

В= А

cos лVfl-,

 

0) = —

arctg — 1=

4 у а у

2у V

 

 

 

 

л

 

 

 

г—У + 1

V - 1 .

п „ .

х < 0;

 

(2.25)

UI - A , , V 7 —

( - , ) «

»

 

^22 A\JLVK{ V k ~“

1)(у -

1) х

к

,

х >0.

 

 

 

Из условия ограниченности потока энергии у края разреза следует к > 2, что отвечает отсутствию потока энергии в точку х = 0. При мини­ мальном значении к 2 < v2 <5/2. Заметим, что к частному решению неоднородной задачи для устранения слишком сильной особенности (если такая возникает в точке х = 0) прибавляется однородное реше­ ние, для которого к < 2.

Положив для четных к А(у - 1) > 0 (у Ф1), имеем и2 ^ 0, о22 ^ 0. Для к нечетных о22 < 0 при и2 ^ 0.

Характерным отличием от до- и сверхзвукового режимов является зависимость значений от скорости и.

Таким образом, в отличие от докритического режима энергия в край разреза не стекает. Вместе с тем, волне расширения (потенциалу ф) соответствует поток энергии к берегам разреза:

A2\iu3ba\2(\K- I) 2 х

2

х (- Х)2У>*~А sin2nvK, *<0, 7>0,

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

точке

^ = 0^ здесь она вся уносится

сдвиговой волной

с носителем

lyl < -

х/у/у, * ^ 0, отраженной от свободных берегов разреза.

Если край разреза движется в

противоположном

направлении

(трещина укорачивается: - с, < и < - с2) в отличие от (2.24) следует положить (у > 0) ф = (- х + Vv уУкЯ (- х + у yfÿ). При этом = к - 0),

к >2.

Для задачи II 0) в (2.25) заменяется на 0) - 1/2.

§ 5.3. Расщепление фундаментальных решений

Решение обобщенной плоской задачи о динамике плоской трещины в линейно-упругом теле, так же как и в статике, можно представить суперпозицией решений трех основных задач (см. задачи I, И, III). Каждая из таких задач является смешанной: на части границы полу­ пространства задана одна из компонент перемещения, на остальной части - соответствующая компонента напряжения (еще одна компо­ нента напряжения задана на всей границе) - см. условия (2.1.10) — (2.1.12), (2.2.1) - (2.2.6). Однако в отличие от статики здесь точки раздела граничных условий (края трещины) движутся, причем в общем случае с переменной скоростью.

Кроме указанных условий на границе полупространства задана линейная связь типа (1.19), (1.26), например в задаче I

uL/(s, 4 = s ^ s ,

4 o £ g s , 4 ;

ui t, x j = S2£

(3.1)

o2£ t, x£.

Для того чтобы решить такую задачу, достаточно найти компонен­ ты, связанные соотношением типа (3.1), на границе полупространства, поскольку еще одна компонента напряжения задана (для задачи I о12= 0). Полное решение выражается свертками. Для той же задачи

*19 -*2) Х\9 * 2) * * ^22^ -*l)>

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

Однако „уравнение в свертках” (3.1), где S22(t, x t) = S22(t, x v + 0),

содержит две неизвестные функции:

и2(хг ^ о) и

е(о). Здесь

0) - область, где задана функция о22, -

берег трещины. Метод решения

таких уравнений для случая, когда область w не является изолирован­ ной точкой (0 < t), основан на факторизации - разбиении изображения

S на множители

SLF= S^FS^ F9 что эквивалентно представлению

5= 5+ * * 5_. Цель

факторизации - расщепление фундаментального

решения на функции, отвечающие направленным волновым воз­ мущениям.

Рассмотрим изображение SLF(s, q) - аналитическую функцию

переменных 5и Q с особыми точками

 

 

 

sm= iqum 9

Qm= “

is/um, um = const Ф 0,

 

 

гп=1, 2, . . . , п ;

Ü1 < Ü2 < . . .

< и п.

 

(3.2)

Предположим, что ее можно представить в виде произведения

 

 

SLF(s, q) = П S%F(s, q),

 

 

 

 

где каждый сомножитель

S%Fимеет свои особые точки, отвечающие

скоростям о т в диапазоне

0

Полагаем, что функции

SFF

на плоскости q- медленного роста

(существует такое число

V,

что

IS^I < IQIVпри I9 I -

~).

 

 

 

 

Покажем, что 5 ^ t, х) -

функция с носителем в секторе и~Kt ^ x ^

на плоскости

х, L Как следует из условий (3.2), при

R e s >0

на вещественной оси (/особых точек нет. В формуле обращения

 

 

1+о° a+i<x>

=Г Î $ Sifts, q) et-iVdsdq

4 л12?4* 0 . IOO

положим s= i(/c+ д Тогда

J

+ 0 0 0L*i°°

 

SH{t,$ =

j S SFI{ p+ iqc, q) <bdpdq.

4 Ti

t -00 a-ioo

 

<ï>=exp[pf-K j(x-c4].

(3.3)

При этом особые точки на плоскости q(3.2) переместятся:

Чт=~ »(Р+

ЩпА^т,

ip

Ят=-------------- (Re р >0).

 

 

с

Если с <и~ то, деформируя путь интегрирования в формуле

обращения для преобразования Фурье (3.3) в верхнюю полуплоскость q, где особых точек нет, обнаруживаем, что SK= 0 при х < et. Если же с> и £ , то, деформируя путь в нижнюю полуплоскость q (при с > особые точки расположены в верхней полуплоскости), получаем: SK= О при х > e t. Поскольку значение с произвольно, приходим к выводу

SK= 0 (x < v z t, х > и *t).

Если Ок = ик = ик, носитель SKсосредоточен на луче х = о Kt:

= 2 £кп(0^П<Ч* ““

п=0

SKF= Z ( - *4 )"^ n (s "

W K)-

 

 

 

n=0

 

 

 

 

 

В частности, если

 

 

 

 

S%F= (s -

iquK)\

 

 

 

 

ТО

 

 

 

 

 

с * - 1

 

 

 

(а. =#од ,. . . ) ;

Г(-

X) б ( М " х)

 

 

V

 

 

 

 

(3.4)

Щ+х/ик

 

 

(X = л = 0,1,. . .),

 

N

 

 

 

 

Щ+х/ок=

I

ô ( x / n" m^

r

m

2

 

т=0

 

 

 

 

 

/ а

д

и

 

х).

щ +х/и к* */(f,x )= I—

+ и к—

1 / а

Носитель

обобщенной

функции

D^+jc/Ukточка £ = х = 0, которая

лежит на луче х = о Kt .

Рассмотрим теперь свертку нескольких функций SK- произведе­ ние их изображений. Поскольку произведению принадлежат все особые точки сомножителей, справедливо утверждение: носитель свертки (по t и х ) представляет собой объединение носителей сверты­ ваемых функций SKи областей (секторов) между ними.

Далее, пусть Q{t, х) = 0 при х < l(t), f(t, х) = 0 при х < uf, / = dl/dt < о

(и = consj). Тогда

 

 

Q(t, х) * * f{t, х) = 0

[х < l{t)].

(3.5)

Аналогично, если Qit, х) = 0 при х > lit), fit, х) = 0 при х > of,

/ 2* и, то

Справедливость утверждения (3.5) можно доказать так. Представим свертку (3.5) в виде

о * * м

s

т, х - ю><

-о о

О

 

Х Щ - / ( т ) № - É- u(f- т)]</тс£,

где множители - функции Хевисайда - введены с целью подчеркнуть расположение носителей свертываемых функций. Вклад в интеграл будут давать те области интегрирования, где подынтегральное выраже­ ние отлично от нуля: £ > /(т), х - - ит. Отсюда для x < /(f) получаем

0 <É- /(т)</(0- /(т)- u(f- т).

(3.7)

Но / ^ о, поэтому правая часть соотношения (3.7) не больше нуля. Противоречивость неравенств (3.7) означает, что при х < l(t) подынте­ гральное выражение везде равно нулю и, следовательно, равна нулю свертка (3.5). Точно так же доказывается и равенство (3.6).

Для антиплоской задачи факторизация проводится без труда.

Действительно,

функция S%%(s, g),

как видно

из выражения (1.13)

при х 2= 0, определяется так:

 

 

 

5 LF= — — (g2+ bjS2) -172.

 

 

 

33

д

 

 

 

 

Ее можно представить в виде

 

 

 

S ^ S ^ S ^ ;

 

 

 

 

 

 

J T

+iqcj'1'2.

(3.8)

S^F= (s -

iqcj-1'2 V c T , S F =- - - --( s

Функция SLFимеет две особые точки g = ± isb2, первая из которых

сохранилась для функции SFF, вторая - для функции SFF

 

Для дальнейшего потребуются

функции P+F= 1/SiF Оригиналы

указанных функций определяются формулой (3.4) при к = -

1/2, к = 1/2:

у/с7

п 1'2

 

 

 

S--------- ---- — т=" б(с2/+ х);

 

 

 

ц

Vu

 

 

 

V eT c1

à(c2t~ х);

 

 

(3.9)

5+ =

Г

 

 

7

 

 

 

цt; 3/2

Р_ =—=г

г— è(c2t+ х);

 

£—3/2

 

 

 

Р+= ------7 = -

ô(c2f - х)

(ô(—С2Г -

х) = 0(с21 + х)).

 

2у лс2

 

 

 

В результате факторизации функции SFF фундаментальное реше­

ние (1.18) длях2 =0

 

 

S33

1

H (t-b 2\x\)

 

-----

—= = = = - = S. * * S

(х = х 1)

 

лр

у *2 Ь2х2

 

оказалось расщепленным на две функции S+, S_, каждая из которых

отвечает

направленному

волновому

возмущению. Если носитель

533 - расширяющийся отрезок - c2t < х < c2t (сектор на плоскости х,

t), то носитель 5+ - точка x = c2t, носитель S_ -

точка х = -

c2t (лучи

на плоскости х, t).

I функция

S22

(1.26)

Перейдем к плоской задаче. Для задачи

(в дальнейшем нижние индексы опускаются)

имеет особые

точки

q = ± isb1', q = ± isb2 (точки ветвления) и q = ± islcR (полюса). Предста­ вим ее в виде

SLF{s, q) = SL/ D l f ,

 

y[b\s2+ q2

S%F(S, q) = -----------

 

 

s2/c2R+q2

D L F {S, q ) = JL

s2(s2/c2 + q2)

(3.10)

\DLF\eid.

1 - v

R

 

При таком представлении указанные выше полюса не являются особыми точками функции DLF (устранимые особые точки). Для нее сохраняются лишь точки ветвления. Функция SFF, которая после заме­ ны параметра Ьх на blic совпадает с приближенным описанием (1.30), легко факторизуется разложением на множители. Основная задача состоит в факторизации функции DLF.

Полагая s > 0, проведем разрезы на мнимой оси q: [ - ib2s, - ibts]9 [ib15, ib2s]. Тогда R ~ b 2s2q2!(1 - v) и, как видно из (3.10), \DLF\ 1, 0 “»0 (\q\ -* °°). Кроме того, функция DLF вне особых точек не обра­ щается в нуль. Поэтому факторизацию DLF= DFFDFF [где имеет лишь две особые точки из четырех указанных - в нижней (верхней) полуплоскости q] можно провести с помощью интеграла типа Коши [61], положив

оо

d +F= exp [Ф+ (s, g)], Ф (s, q) = ± ~

f

l n [ D s^£)]s -, ^

2Л1 J

t - q/s

 

-CO

 

(при вычислении интеграла Imq > 0 для Ф+ и Imq < 0 для Ф _).

Учитывая, что In DLF(s, g) -* 0 при q -> ± 00 и что подынтегральное выражение представляет со­ бой аналитическую функцию комплексной пере­ менной £ с особыми точками £ = ± ib19 £ = ± ib2 (точки ветвления) и £ = q/s (полюс), контур инте­ грирования можно деформировать, проведя его по берегам отрезка [- ib2, - ib1] для функции Ф+, когда Img > 0, и по берегам отрезка [ib19 ib2] для функции Ф_ (Img < 0). Контуры и направление их обхода показаны на рис. 5.6.

Заметим, что на берегах указанных отрезков, как следует из выражения для функции R (1.24),

R = {bis2 + 2g2)2 + 4ig2 V b 2s+2 g2 7 - bfs2 - g2,

где радикалы положительны, знак минус относит­ ся к правому берегу в верхней полуплоскости и к левому берегу в нижней полуплоскости, а знак плюс - к остальным двум берегам. Учитывая это и полагая £ = ia, можно преобразовать равенство (3.11) к виду

ibo

§

К/ ibf

S N -ibf

~ib0

Рис. 5.6.

 

 

1 Г

<p(oc)doc

J

cp (a )=

4a2 yjb\- a 2 Va2- b\

D+f = exp

л

J + iq/sa

arctg-

(2a 2- b2) 2

 

 

 

 

 

 

 

 

 

 

(3.12)

В соответствии с

формулами (1.16),

(1.17),

где

gL= 1, g = 6'( f),

 

 

 

 

b2

 

 

 

 

 

Л

 

I

1

Г

<p(a)da

 

 

 

 

= -------exp

J

a +t/z

 

 

 

 

 

2nz

л

 

 

 

 

1

Г

<p(a)da

ь, .

у

 

Г <p(«)da

[

л

1

J

n J

a

+ f/z

 

b.

+ t/x

a

 

 

bl

 

 

 

 

 

 

 

 

 

 

+ i <p{t/x)H(± t/x-

b/)H(b2 + t/x)

(y = + 0),

+

± iif{t/x)H{± t/x- b1)H(b2 + t/x)

(y = -

0)

( < p

( - a ) =

<p(a));

 

 

 

 

 

 

D°±= à(x)H(t) + Q±(t/x);

 

 

 

 

 

Q+-

 

 

1

. .

<P(a)daГ

 

 

 

 

 

exp —

V.p.

a + t/x

ф +;

 

 

 

nlxl

л

 

J

 

 

 

ь.

 

4

a2 \/bl -

а

2 / а 2 -

Ь2.

- ;а 2) ( а

2

-

b 2)

sin ср(а) =

-V, ( 2 a 2

-

b l ) 4

+

1 6 a 4(b 2

D± = -J -

D ® - f i ( x ) 6 ( f ) + 4 -

Q± .

 

 

 

 

 

(3.13)

ot

 

 

ot

 

 

 

 

 

 

 

 

Видно, что носитель функции D+ -

сектор

c2t^ x^ \ ^f,

функции

D_ - сектор - c Л ^ х < -

c2t,

как

это

и определяется

положением

особых точек

на плоскости q.

 

 

 

 

 

 

 

Функцию 5 ^ (ЗЛО) можно представить в виде

 

 

 

S F =

 

 

 

 

 

 

 

 

 

 

 

т = 1

 

 

 

 

 

 

 

 

 

 

 

501 = ------- 7— Vs + 'СхЧ ,

$ 02

----- Г— ;

 

 

 

 

й VCt

 

 

 

 

s+ !дсй

 

 

 

 

çLF_.

5 ^ = —T=r Vs-lc^q".

ù03 ~ s - iqCjj

VCi

Оригиналы, отвечающие этим функциям, будут

 

1 - V

Г 3' 2

ô(c^ + x);

■*01

2р V q

V ÎT

 

 

S02 = Cflôfof + x),

503 = Cflôfof - x),

 

t-3,2

0(ctf - x ) .

$04 Г---

 

2 vnc1

 

 

Функциям P0^ f = 1/S0^ f соответствуют

Vc7

c 1/2

Л ,1= " i - v

ôOV + x);

VÎT

1

P02 = — [ô '(f)ô (x)-Cjî6'(x)ô(0];

(3.14)

(3.15)

P 03 = 1— [ ô '( f ) ô W + C tfô 'U )ô (f)];

V c7 c 1/2