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

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

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

(5.17)

Видно, что при к = 0 (и только при к = 0), когда функция С2 опре­ делена выражением (5.16), однородные решения удовлетворяют нужным условиям. При этом функция f(t) остается неопределенной. Она соответствует сосредоточенному источнику энергии, движуще­ муся со сверхрэлеевой скоростью (сд < / < с2).

Как следует из формул (2.2.25), (4.31), (5.17), интенсивность излу­ чения энергии из края трещины связана с функцией f(t) следующим образом:

- г — (o/Cjî- i y 2m m n t ) =

(CR< и < Са).

§ 5.6. Автомодельная задача

Рассмотрим задачу о трещине, начальная длина которой равна нулю. В этом случае метод, основанный на последовательном учете взаимного влияния напряжений изложенный в §5.4, не может быть применен. Если задача автомодельна: и = tKu(t/xl9 t/x2), то ее можно решить другим способом, например на основе метода функционально-инвариантных решений [32, 31, 117, 122]. При этом используются решения уравнений теории упругости, определенные и вне плоскости трещины. Однако для приближенной модели (1.30) состояние вне указанной плоскости не определено, постулирована лишь связь (4.1) между перемещением и напряжением в плоскости трещины. С целью получить решение как для точной, так и для прибли­ женной моделей, воспользуемся другим методом, основанным на введении аналитических представлений, определяемых формулами (1.16), (1.17). Для реализации этого метода достаточно соотноше­ ния (4.1).

Полагаем, что функция SLF(s, q) (4.1) удовлетворяет условиям

SLF(s,sp) =-

sLF(l,p);

s

— SLF(1, -

it/z) = S0(z/t)----- const (z = x + iy, Ы -♦ °°);

z

z

 

(6.1)

50[(x + iO)/t] = - 5 0[(x - iO)/t] ( - u 0t < x < u 0t9 o 0 > 0),

причем функция S0(z, t) ограничена и отлична от нуля в некоторой

полосе + 0 ^ у <

а ( - а < у < -

0),

- (и 0 - e)t ^ х < (о 0 - e)f при любом

значении е > 0.

Введенные в

§

5.1 функции SLF(s, q)9 отвечающие

рассматриваемым задачам, удовлетворяют условиям (6.1). Для задач I, II и0 = ср , для задачи III о 0 = с2.

Пусть напряжение o(t, х) представимо в виде

о = Dn(t) * Q{x/i).

(6.2)

Обобщенная функция Dn{t)9п = 0, ± 1,. . ., определена так:

D0(f) = ô(t),

Dn(t) = 0

(t< 0),

Dn_m(t) = d™Dn(t)/dt™.

При этом свертка с Dn(t) означает - n-кратное дифференцирование по t (п < 0), n-кратное интегрирование (п > 0) или тождественное преобра­ зование (п = 0), Dn(t) * Dm(t) = Dn+m{t). Из представления (6.2) следует, что L F - преобразование функции o(t, х) имеет вид

0LF(s, q) = 5"n" 2Q0(s/q),

(6.3)

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

uLF{s, q) = s~n~3Q0(s/q)SLF(1, q/s),

(6.4)

и, следовательно, перемещение представимо в виде

 

u(t, х) = Dn+1(t) * Qi(x/t).

(6.5)

Значение п в (6.2) - (6.5) определяется зависимостью для напряже­ ния о_ или перемещения и+ (для них это значение предполагается

одинаковым).

Так, если о_ = ô(u0f - л:), то

п = - 1 . Действительно,

Q= H(o0-x/t),

о_ = d/dtQ{x/t) = D.^t) * Q(x/t).

Для o_ = const ( - и J <

< x < v +t) л = 0.

Используя формулы (1.16), (1.17), находим связь между аналитиче­ скими представлениями для функций и и о (обозначаем их теми же

буквами, но вместо переменной х для и и о пишем z для их аналити­ ческих представлений):

u{t, z) = ~ —-----D„+2(0 *

^ ' ■ - т И т

2niz

Отсюда

u{t, z) =Dn+2(t) * 5 0f — [[D.n.iW * o(f, z)][,

t[D-„-2(t) * u(t, z)]

0{t,z) = Dn+1{t) *

(6.6)

 

zS0(z/t)

Полагаем, что для функций u{t, х), o(t, х) существуют представле­ ния Коши, которые даются формулами (1.16), (1.17), и введенные здесь аналитические представления - представления Коши. В связи с этим

o(f, z) = 0(l/z), u{t, z) = 0(l/z) (Izl °°). (6.7)

Данное предположение, очевидно, выполняется в рассматриваемой задаче о трещине, поскольку при любом конечном значении t функ­ ции и, о - финитные.

Пусть заданы функции

 

 

u(t, х) = u+(f, х)

(х < -

иtt,

х > О2t);

o(f, х) = o.{t, х)

( - иxt < х < o 2t);

(6.8)

О^ U 1,2 ^ U 0>

 

 

 

и требуется найти функцию

 

 

u = u_(- i V < x < u 2f),

о = о+

( x < - r U tf, х > и 2t).

При сформулированных условиях (6.6) - (6.8) решение задачи не единственно, так как существуют нетривиальные однородные реше­ ния. Например, для п = - 1, положив и+ = о_ = 0, можно записать

u(t, z) = u0(z/t) = i[(z/f + о t)(z/t- u2] - 1/2 =

(

i[(xA + U j)(x A -u 2)]"1/2sgnx ( x < - t V , x > u 2t, y=0),

(

± [(*/f+ U i)(t> 2 -.x /il)]-1/a ( - u ^ < x < u 2f, y = ± 0).

При этом

 

 

 

 

uol^]=lim[uo{z/t)-uo(z/t)] = 0

( x < - i V ,

х > и 2t)

\ 11 у*+о

 

 

 

 

и, так как при -

< x < u2f

u0(z/t) = - u0(z/t) (у = + 0), в соответ­

ствии с равенствами

(6.1), (6.6)

на

том же

интервале a(t, z) = a(t, z)

(y = + 0), a следовательно, o_ = 0. Вместе с тем u_ Ф 0, o+ Ф 0, т. e. нетри­ виальное однородное решение действительно существует.

В связи с этим введем, как обычно, условие непрерывности пере­

мещения границы полуплоскости

 

u_(f, + u12t ± е)~ иЛ*> + ui 2f + е) "" 0 (е -►0),

(6.10)

которое для задачи о трещине ограничивает поток энергии. Уравнение (6.6) решается, по существу, тем же методом, что и аналогичное урав­ нение для статики (см. § 2.2). В соотношение (6.6) введем функцию fi

D_n. 2(t)*u(t,z)

S0(z/t)

i(0 * o(t, z)].

(6. 1 1 )

fi (t, z) =

 

[D.

zu0{z/t)

 

tu0(z/t)

 

 

Скачок функции fi

на

вещественной

оси fi(f, x) = fi(t, z) -

fi(f, z)

(y = + 0) при x < - u1t и

при x > u 2t известен. Функция [zu0(z/t)]~l

в указанных областях непрерывна, поэтому он определяется так:

 

fi(f, х) = [D_n_2(t) * u(t, x)]/[xu0(x/f)] =

 

 

= [D_n. 2(0 * u+(f, x)]/[xu0(x/t)\.

 

(6.12)

Второе равенство вытекает из первого, поскольку скорости концов трещины - и1 ^0, и2 >0, и если при некотором значении t = t x коорди­ ната данной точки x > u 2tx или х < - о xtl9 то эти неравенства сохраня­ ются и для t < t 19 вследствие чего свертку по времени с функцией и можно заменить (для указанных областей) на свертку с функцией и+.

На интервале - < х < u2f непрерывно отношение S0{z/t)/u0{z/i). Поэтому на том же интервале, как видно из равенства (6.11),

 

S0[{x + f0)/f]

(6.13)

O f t * ) -

W-n-i(t) * o(t, x)].

tuo[{x + i0)/t]

 

Здесь, вообще говоря, нельзя заменить о на 0_: если функция о интегрируется (п = - 2, - 3,. . . ), то для t = tx свертка будет зависеть от напряжений o(t9х)9 которые для t < t 19 точнее, для f < т а х ( - х / и ^ х/и2), обозначены как о+. Однако при п > - 1, когда свертка в правой части равенства (6.13) не меняет носителя функции o(t, х)9

D-n-i(fl * ®((> x) = D_n_j(l) * o_(f, x)

( - utt < x < o 2t),

поэтому скачок Q(t,z)

известен

и в указанном интервале. Именно

этими значениями п > -

1 и ограничимся.

Остается учесть возможные

скачки в точках х = - и £ x = u 2t9

существенные для дальнейшего лишь как сингулярные обобщенные функции. Поскольку последние сосредоточены в указанных точках,

а функция

ÏÎ(U), как следует из ее определения (6.11) и выражений

(6.5), (6.9),

представима в

виде О(t,z) = D_£t) *M(z/t) (М- некоторая

неопределенная функция), можно написать

 

 

fi(U )=

i

+

и2t)\

 

m = О

 

 

 

(X = - I V , х = и 2t),

 

 

(6.14)

где Am'Bm- неопределенные постоянные.

что

fi(f, z) = 0(l/z) при

Из выражений (6.5),

(6.9), (6.11) видно,

\z\ -+ 00 и,

следовательно,

fi(f, z) - представление

Коши для Q{t, х).

Таким образом, функция

Q{t, z) определяется интегралом типа Коши

по ее скачкам (6.12), (6.13), (6.14). Для перемещений и напряжений, используя равенства (6.11), можно написать следующие формулы:

u(t,x) = lim[ii(t,z)- u(t,z)];

у-*+О

a(f,x) = lim [o(f,z)- o(f,z)]; y-*+0

u(^ z) = Dn+2{t) *[t-'w(z/%

 

 

 

 

U2 t

 

 

 

 

 

 

 

K <z\

zujjz/t)

(

Г

S M l+ iW

 

 

 

 

W

4

2ni

l

J

u j(g + i0)/fl(6 -

z)

 

 

 

 

 

 

•Vit

 

 

 

 

 

 

 

x [D .„ _ i(f)* o .(U )]^

+ f j

 

 

d\~

 

 

 

 

b

 

 

 

v

 

 

 

 

 

 

 

X

m!(" У

 

 

 

 

 

B

 

 

 

-

(z + и ^ ) т+1

( z -

t>2i)m+1

 

 

 

 

m = 0

 

 

 

 

 

 

 

 

 

 

 

z) = Dn+1(0 * U[0 _n_2(f) * * u(t, z)]/[z50(z/i)]} =

 

 

 

= Dn+i(0 * M z/9/[z5o(z/0]} •

 

 

 

 

(6.15)

Здесь

y -

область

интегрирования:

£ < - u^, £ > u2f;

сумма

по

m - аналитическое представление обобщенной функции -

Q(t, х)

для

точекх = -

х = и2/(6.14). Из условия (6.10) следует к= п. Искомые

функции и, о определены формулами

(6.15) с точностью

до

2(n + 1)

произвольных постоянных (п = -

1,

0,

. . .)• Функция uQ(z/t)

и

все

слагаемые в выражении для w(z/Y), кроме интеграла по области у (6.15),

непрерывны на у.

Отсюда, как

непосредственно видно из формул

15 — 171

-

225-

(6.15), последние определяют перемещения на у

(и = и+)

независимо

от значений упомянутых

постоянных. Напряжения о на

интервале

- u 1t < x < u 2t при п> -

1 от этих постоянных

зависят,

поскольку

свертка с Dn+1(t) (6.15) означает (п + 1) кратное интегрирование по L По той же причине, по которой нельзя было заменить о на о . в форму­ ле (6.13) при п < - 1, можно утверждать, что вклад в интеграл на ука­ занном интервале будут давать значения о = о +, соответствующие однородным решениям.

Из сказанного следует, что произвольные постоянные должны определяться условием (6.8) относительно функции о. А именно, они должны быть выбраны так, чтобы на интервале - i)tt < х < u2t напря­ жение о, рассчитанное по формуле (6.15), совпадало с заданным: о = о_.

Из соотношений (6.15) следует, что при 0 < х < u2t x/U2

(6.16)

о

где Фх(хА) - скачок функции Ф^z/t) = w{z/t)/[(z/t)S0{z/t)] при x > u 2t.

Второе слагаемое в правой части выражения (6.16) должно быть равно

нулю. Следовательно, 1/и,

ш = 0 ,1,. . .,п .

(6.17)

о

Система линейных уравнений (6.17) определяет постоянные Вт. Точно так же из условия о = о_ при - и J < х < 0 получается система

уравнений относительно постоянных Ат: 1/U

оÎ

 

 

 

где

Ф2(- x/t) -

скачок функции

Ф2( - z/t) = w{z/t)/[(z/t)SQ(z/t)\ при

х < -

utt.

 

 

В симметричной задаче, т. е. в случае когда и2 = и2 = и, o_{t, -х ) =

= o_(t, х) имеем

Ат = Вт. При этом,

как следует из формул (6.15),

при х -►иt + 0 коэффициент интенсивности напряжений выражается следующим образом:

и

(n = - 1);

Vu2 - а 2

О

Пусть к берегам трещины приложены постоянные сосредоточенные силы: o_(f, х) = - b(x)H(t) (п = - 1). Тогда из формул (6.18) получаем

So(0)

1

1

s 0(u)

v W

/ 0(и)\/лГ’

где / = текущая полудлина трещины, функция / 0 определена формулами (4.35)- (4.37) (см. рис. 5.3). Сравнивая с формулой (2.2.17) для р(£) = ô(|), видим, что здесь отношение статического коэффициен­ та интенсивности напряжений к динамическому равно / 0(и).

В случае равномерно распределенной нагрузки o_(f, х )= — H (ut-

-Ixl)H(t) получаем зависимости

=- q ^ D ^ R ^ y fïw t;

КцI = q2]/im/E(q2);

 

 

R jo ) = (2 - blu2)2 - 4qtq2 ;

q1? = ^l I -

b ^ u 2 ;

Dt « [(4b2u2q2 + Ь*\)4Ж(Ч1) -

4blu2qîK(q2) +

+ 8q \E (q 2) ~ (4q \ + 4q§ + b4u4)E(qi)]-‘ ;

 

D2 = [(4^ u 2q| + Ь<и<)Щ2) -

A b ^ q ^ q ,) +

+ S d i^ i) - (8q§ + b^u4)E(q2) ] - i ,

 

Е(Я1?) - - Яу Я*(и) - 1

(u - 0).

(6.19)

Здесь К и £ - полные эллиптические интегралы первого и второго рода.

Для приближенного описания (1.30)

*1,и= Кц| = Ч(1,2)*Чэ

V ^ ;

 

 

= [Чз^(Ч(1^)*) + (и2/4 ) (1 -

cj?b212)*) X

 

 

 

.____ !____

(6. 20)

XK(4( V )*)]"1 ,

4(12)* = Vl - u2b(21>2)* ,

4з = 1 - и2/с| .

На рис. 5.11 представлены графики отношений К(и)/К(0)9построен­ ные по формулам (6.19) и пронумерованные соответственно рассматри­ ваемым трем задачам, а на рис. 5.12, 5.13отношения KffKi, Щ/Кц соответственно (v = 0,3). При расчетах принимались значения парамет­ ров Ь(1#2)*, определяемых формулами (1.31).

Как видно из рис. 5.4, 5.10, 5.13, в случае задачи И приближенное

описание приводит к довольно значительному отличию от точного результата даже для умеренных скоростей. Лучшее соответствие получается при выборе параметра Ь2* из условий совпадения функций 5 ц и 5i ю в среднем на интервале 0 ^ x/t < сх [86] (принимаются глав­ ные значения соответствующих интегралов). Графики зависимостей отношений btJb2 (кривая 1), b2Jb2 (кривая 2) от коэффициента Пуас­ сона показаны на рис. 5.14, а соответствующий этому график отноше­ ния коэффициентов интенсивности напряжений Кцц, определенных по формулам (6.19), (6.20), - на рис. 5.15.

Рис. 5.11.

§ 5.7. Некоторые задачи о динамике пояубесконечной трещины

Начнем, как обычно, с антиплоской задачи. Пусть, полубесконечная трещина < /(/^расширяется со скоростью I < с2 (/ > 0). Как было показано выше (3.9), функции S+, Р+ для антиплоской задачи сосре­ доточены на луче x = c 2t. При этом общее решение нестационарной задачи относительно напряжения на продолжении трещины о23 = о+, соответствующее импульсу о_ = ô ( / - f 0) ô (x - x 0), х 0 < /(0), дается формулой (4.9) при С = 0

С2

/(то) - *0

c2{t - t0)),

(7.1)

о+

ô (x - х 0 -

л(х - Х 0)

Х - /(Т 0)

 

 

где т0 подчиняется уравнению

 

 

/(т0) “ х0 "*

”” ^о)*

 

(7*2)

Итак, при действии сосредоточенного импульса напряжение на продолжении трещины представляет собой сосредоточенную силу, движущуюся со скоростью и и изменяющуюся во времени. Параметр т0, удовлетворяющий уравнению (7.2), соответствует моменту, когда фронт волны, излученной импульсом, достигает края трещины и по­ рождает указанную сосредоточенную силу.

Перемещение в плоскости х г = х, х 2 определяется сверткой фунда­ ментального решения (1.18) с напряжениями о_ + о+[о_ = ô (f - f0)ô (x -

- х 0), о+-

определяется зависимостью (7.1)]. Вычисления

приводят

к формуле

 

 

 

 

п = и3

-------- [f{t, х 19 х 2)] 1/2

х 19 х 2)] -

 

 

Л|1

 

 

 

fit, х 19 х 2) = it - tQ)2 - b2[ixt -

х 0)2 + X2];

b2 = 1/с2.

 

Таким образом, при действии на берега трещины указанного

импульса

влияние ее движения

сказывается

лишь через

параметр

т0 - время прихода фронта сдвиговой волны, излученной импульсом, к краю трещины. Последующее движение трещины (0 ^ / < с2) не вли­ яет на решение.

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

 

В круге,

ограниченном

ее

фронтом

 

( * 1

- , * 0

-

C 2 ( T 0 -

*о))2

+ А

< b l i t -

 

- т0)2, отраженная и прямая волны

 

взаимно

 

уничтожаются,

поэтому

 

перемещения

и

напряжения

там

 

равны нулю (см. рис. 5.16, где возму­

 

щенная область заштрихована).

про­

 

 

Решение,

соответствующее

 

извольной нагрузке берегов трещи­

 

ны,

получается

из приведенных

 

выше суперпозицией. В частности,

 

зависимость для напряжений на про­

 

должении трещины имеет вид

 

 

х - х 0

\ / / ( Т0) - Х 0

 

 

 

 

---------- , х0

------------------ dx0 ;

 

 

С2

 

I

 

%

х 0

 

 

 

 

 

 

 

 

 

 

 

 

(7.3)

При вычислении асимптотики напряжений

l(t) + 0) уравнение

относительно

легко разрешается. Действительно, при этом ig

t9

несли в данный момент скорость трещины l(t) непрерывна, то /(т0 - - l(t) l(t)(t - т0). Отсюда и из указанного уравнения следует

T0~f -------/(То) - 00 ~ (То -

0/(0

 

 

 

с2

- lit)

 

 

 

(7.4)

(х -* l{t) + 0).

 

 

 

 

 

 

 

 

 

Учитывая это, можно записать

 

 

 

 

________________ w

 

 

1

/

1 - / ( 0 / с 2 Г о

 

X

о,

У

х -

/ ( 0 J

0

*0

Л

 

 

>/f(0 —л:0

 

x~c2t

 

 

dx0

/(0 + 0).

 

 

X

 

 

 

X “* XQ

Заметим, что при постоянной скорости трещины асимптотические равенства (7.4) становятся точными.

Пусть берега трещины при t = 0 нагружаются постоянными равно­ мерно распределенными напряжениями о_ = - о = const. При этом из формулы (7.3) получаем

2о

fc2t -

х + /(т0)

fc2t - x

+ ljт0)

о+= -

 

 

— arctg

lit о)

1

 

X -

/ ( Т 0 )

X -

2о

tjc2-

/(f))

lit) + 0).

 

л

X -

Ut)

(х -