книги / Механика трещин
..pdf(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) |
(х - |
|
|
|
|
|