О расслоении полосы по границе раздела упругих свойств часть 3. Точное аналитическое решение для частного случая Текст научной статьи по специальности «Физика»

Устинов К.Б. О расслоении полосы по границе раздела упругих свойств. Часть 3. Точное аналитическое решение для частного случая // Вестник Пермского национального исследовательского политехнического университета. Механика.-2016. - № 3. - С. 222-240. DOI: 10.15593/perm.mech/2016.3.15

Ustinov K.B. On delamination of a stripe along the boundary between two elastic layers. Part 3, exact analytical solution for a particular case. PNRPUMechanics Bulletin. 2016. No. 3. Рр. 222-240. DOI: 10.15593/perm.mech/2016.3.15

ВЕСТНИК ПНИПУ. МЕХАНИКА

№ 3,2016 PNRPU MECHANICS BULLETIN

http://vestnik.pstu.ru/mechanics/about/inf7

Э01 10.15593/регш.шесЬ/2016.3.15 УДК 593.3

0 РАССЛОЕНИИ ПОЛОСЫ ПО ГРАНИЦЕ РАЗДЕЛА УПРУГИХ СВОЙСТВ ЧАСТЬ 3. ТОЧНОЕ АНАЛИТИЧЕСКОЕ РЕШЕНИЕ ДЛЯ ЧАСТНОГО СЛУЧАЯ

К.Б. Устинов1, 2

1 Институт проблем механики им. А.Ю. Ишлинского РАН, Москва, Россия, Московский авиационныи институт, Москва, Россия

О СТАТЬЕ АННОТАЦИЯ

Получено аналитическое решение однородной задачи изотропной теории упругости в постановке плоской деформации о полубесконечной трещине, проходящей вдоль границы, разделяющей два упругих слоя одинаковой толщины с отличающимися свойствами, для равного нулю второго параметра Дундурса (безразмерной комбинации упругих констант материалов). Данная задача является частным случаем задачи о полубесконечной трещине, проходящей вдоль границы, разделяющей два упругих слоя с произвольными соотношениями толщин и упругих свойств. Решения для случаев трещин нормального отрыва и сдвига в предположении возможности пренебрежения влиянием нормальных напряжений на сдвиговые смещения и сдвиговых напряжений на нормальные смещения даны в первой и второй частях статьи (Вестник ПНИПУ. Механика. 2015. № 4; 2016, № 2). Рассматриваемое решение для частного случая получено путем применения преобразования Лапласа и сведения задачи к матричной задаче Римана-Гильберта. Ограничение, состоящее в необходимости равенства толщин, и ограничение на сочетание упругих постоянных материалов определяются методом решения, позволяющим осуществить факторизацию матричного коэффициента задачи в замкнутой форме.

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

© ПНИПУ

Получена: 22 марта 2016 г. Принята: 29 июля 2016 г. Опубликована: 30 сентября 2016 г.

Ключевые слова:

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

© Устинов Константин Борисович - доктор физико-математических наук, доцент, старший научный сотрудник, e-mail: ustinov@ipmnet.ru

Konstantin B. Ustinov - Doctor of Physical and Mathematical Sciences, Associate Professor, Senior Researcher, e-mail: ustinov@ipmnet.ru

ON DELAMINATION OF A STRIPE ALONG THE BOUNDARY BETWEEN TWO ELASTIC LAYERS PART 3. EXACT ANALYTICAL SOLUTION FOR A PARTICULAR CASE

K.B. Ustinov

1

A.Yu. Ishlinsky Institute for Problems in Mechanics RAS. Moscow, Russian Federation 2Moscow Aviation Institute, Moscow, Russian Federation

ARTICLE INFO

ABSTRACT

Received: 22 March 2016 Accepted: 29 Jule 2016 Published: 30 September 2016

Keywords:

delamination, interface crack, matrix factorization, elastic clamping, matrix of elastic compliance

The analytical solution has been obtained for a homogeneous problem of elasticity in plane strain setting of semi-infinite crack separating two isotropic layers of equal thicknesses and different but related with a particular restriction (corresponding to vanishing of the second Dundurs parameter) elastic constants. The problem under consideration is a particular case (allowing us to obtain the exact analytical solution) of the general problem on semi-infinite crack separating two isotropic layers of different thickness and arbitrary combination of elastic constants posed. The solutions for the cases of normal and shear crack were obtained under the assumption of the possibility to neglect cross-terms related to the influence of the normal stresses on the shear displacements and the shear stresses on the normal displacements. In this part of the work the solution for the case in question has been derived by means of Laplace transformation; and reducing it to a homogeneous Riemann-Hilbert problem with the matrix coefficient. The restriction imposed on elastic constants and the demand of equal thicknesses are determined by the used method allowing the factorization of the matrix coefficient.

The asymptotical expression has been derived for the relative displacements of the crack faces far from its tip, corresponding to a beam deflection with the boundary condition of the type of generalized elastic clamping. i.e. the proportionality of the displacement and angle of rotation of the clamping point to the total vector and bending moment of the applied load by means of the matrix of coefficients of compliance. The analytical expressions for these coefficients have been obtained. The asymptotical expression for the stress field near the crack tip (stress intensity factor and energy release rate) has also been derived.

© PNRPU

1. Постановка задачи. Исходная конфигурация

В первой части [1] была сформулирована задача теории упругости в постановке плоской деформации о полубесконечной трещине, проходящей вдоль границы, разделяющей два изотропных упругих слоя с произвольными соотношениями толщин и упругих постоянных. Самоуравновешенная нагрузка, характеризуемая главным вектором и изгибающим моментом, предполагалась приложенной на бесконечности со стороны трещины. Путем применения преобразования Лапласа задача была сведена к задаче Римана-Гильберта с матричным коэффициентом. Далее приводится изложение точного решения поставленной задачи для частного случая: слоев одинаковой толщины, на соотношение упругих констант которых наложено некоторое ограничение. Исходная конфигурация задачи представлена на рис. 1.

Рис. 1. Расслоение составного слоя. Геометрия задачи Fig. 1. Delamination of contained layer. Geometry of the problem

Здесь и далее а^., а , а - компоненты тензора напряжений, и, V - компоненты вектора смещения; все величины, относящиеся к нижней полосе -1 < у < 0, обозначаются индексом 1, все величины, относящиеся к верхней полосе 0 < у < 1 - индексом 2.

Для рассматриваемой однородной задачи предполагается, что все поверхности свободны от напряжений,

ауу =аху = °> при У = 1, У = —К и при у = 0, х > 0, (1)

а нагрузка с эквивалентными главным вектором (Т, N) и моментом М приложена на бесконечности, так что (здесь по сравнению с [1] изменены знаки)

-М = | хауу(х,0)ёх, N = |ауу(х,0)ск, Т = |аху(х,0)ёх (2)

—ад —ад —ад

Условия сопряжения на границе имеют вид

ауу = а®, = а%\ и(1) = и(2), V« = v(2) при у = 0, х < 0. (3)

Модифицированные для условий плоской деформации модули Юнга и коэффициенты Пуассона материалов обозначим Е(г),V(г), / = 1,2 соответственно. С обычными модулями Юнга Е0(г) и коэффициентами Пуассона V0(г) они связаны соотношениями

Е (') (')

Е (') =_Е_ Vг) = __(4)

Е (1—V» )(1)' " 1—V;« ()

Путем применения преобразования Лапласа [1] задача сведена к задаче Римана-Гильберта с матричным коэффициентом. Запишем ее для условий одинаковых толщин полосы, равных единице и дополнительного ограничения на сочетание упругих слоев, необходимого для получения аналитического решения, а именно

1 |(1) (1 — 2v(2)) —1(2) (1 — 2v(1) ) =

2 = 0 (5)

Данное условие соответствует равенству нулю второго параметра Дундурса (параметра упругого несоответствия) [9, 10]; при выполнении данного условия при сингулярности отсутствует осцилляция.

Итак, матричную задачу Римана-Гильберта на мнимой оси (р е Ь) для указанных ограничений запишем в виде (без нарушения общности, будем полагать а > 0 - данное условие позволит избежать необходимости рассмотрения различных случаев выделения ветвей функций, содержащих радикалы [5, 6])

( р ) = К ( р ) Ж ( р ), (6)

к (р ) =

^(п р соб р + р)/ ё ар2/ ё ^

ре Ь (7)

—ар2/ ё (п р соб р — р)/ ё

ё = бш2 р — р2, (8)

— —1 Е (2) — Е (!)

— = Е(2)/Е(1), а =——1 = Е(2) Е(1), (9)

л + 1 Е(2) + Е(1)

F ( ) Е(2) Г 5 [у(2) (х'0)"v(1) (х'0)1 - _ (10)

*(р) = 5 {„о (х,0)- и(1) (х,0)|е (10)

F-(P) = ? {°*(^-^ (11)

-Г [ах, (х,0)]

с условиями в опорных точках (в нуле и на бесконечности)

Р-(Р) = {г+-„«+ °(Р'}• КеР - 0-' (12)

Ж-(р)^-р=р | О (р-3/2 )• Яе Р--Г. (13)

Значение ^ = 1 (а = 0) соответствует расслоению однородной полосы [11-13]. Основная сложность состоит в факторизации матричного коэффициента К(Р), т.е. представлении его в виде

K ( Р ) = X+ ( Р ) X-1 ( Р ) • (14)

где функции X± (Р) аналитичны в правой и левой полуплоскости комплексного переменного Р , и ёе! X± (р) ^ 0 в соответствующих плоскостях вплоть до границы. В настоящее

время общее точное аналитическое решение задачи для произвольной матрицы неизвестно, однако для указанных ограничений задача решается применением метода [3-6] (см. также литературу, указанную в [1]).

После нахождения X± (Р) окончательное решение задачи дается с помощью теоремы Лиувилля:

¥+(р ) = X-1 ( Р ) П ( Р ), Яе ( Р )> 0 ¥-(р ) = X-1 (р )П (Р ), Яе (Р )< 0

Здесь П (р) - векторная функция, подлежащая определению. 2. Решение задачи Римана-Гильберта

Матрицы X±1 () могут быть представлены в виде [3-6]

X+1 (р) = Л±1 (р)[1 сЬ(ф(+)±Б(р) вЬ(ф(3+)], X±1 (р) = Л±1 (р)[I сЬ(фP-)±B(р) вЬ(ф (-)].

Здесь I - единичная матрица,

( 1 ар ^

(15)

(16)

Б (р ) = ф-1 (р) ар1 , ф(р ) = у11 -а2 р2 (17)

I-ар -1)

Функции Л±(р), Р±(р) являются решениями скалярных задач Римана-Гильберта

Л-1 (р)Л+(р) = А1/2 (р), р е I (18)

р+(р)-р-(р) = ф-1 (р)е(р), р е I (19)

Здесь определитель А( р) и показатель в( р) матрицы К (р ) определяются как

А(р) = ^1 (р)Х2 (р), (20)

е(р ) = 1]л Щ. (21)

2 X (Р)

Собственные числа Х1 (р), X 2 (р) матрицы К (р ) могут быть представлены в виде

X, (р) = -2с1в3 рА(р/,), Х2 (р) = 2с1в рВ(р/,), (22)

2А (э) = !Ь ^ (!Ь ^ + ^ 1 + а2я2 сЬ 2 я) / (1 - я2 бЬ 2 э),

. Х ___х 7 (23)

2В (э) = (1 - Ы1 + а2э2 БЬ 1 сЬ 1 э) / (1 - э2 БЬ 2 э).

Таким образом, исследуемая задача отличается от задач, решенных в [7, 8], видом функций A (s), B (s), определяемых формулами (23).

С учетом (20)-(22) скалярные задачи Римана-Гильберта могут быть записаны в виде

л:1 (p)Л+(р) = 2i ctg2 p^A(p /i)B (p /i), p e L, (24)

ß+(p ):ß:( p ) = Ф:1 (p )2ln

( лГ„/ЛЛ

: Ctg2 p

V

A ( p / i) B(p / i)

, p e L. (25)

Решение задачи (18) или (24) запишем в виде, аналогичном данному в работах [3-6] через интегралы типа Коши:

, ч лГ2 (1/2 + р/я) , ч , ч р2Г2 (1/2-р/ л) , ч

Л+(Р)^-Г2(-Р7рг'+(Р), Л-(Р)= Р2я,■ Г2(1-рР/я) ■/-(Р), (26) Л(р) = е*р{";Ц>" [А (»)В(»)]Д} ^ (27)

внеинтегральные члены получены факторизацией тангенса через гамма-функцию.

Решение задачи (19) также запишем в виде, аналогичном данному в работах [3-6]:

Р+ (р) = (0+ (Р)- Р* + (Р)- Рг + (Р), Р- (Р) = (0- (Р) - Р*- (р) (28)

Вспомогательные функции

1 dt 1П (ap + ^ 1 :а 2 p 2 )

ß»(p)=4iФт^7) - \jrap ' ß-(p)-:ß"(:p(29)

р ( ,_ 1г58п(/¡)Л = —1 У1 — а2р2 +1 Рт+(р)- 21 ф(0(, — р)~ ¿—Ор21П ар

введены таким образом, чтобы правая часть выражения

Р0+ (р) — р0— (р) = Ф—1 (р) е(р) + Р5 + (р) + Р5 — (р) + Рт + (р), р е Ь удовлетворяла условиям Гёльдера. Тогда решение задачи [3-6].

РЛр) = — ¿1(1 + а252)—1/21п ^

¡А (5) а5 сЬ 5

.2 „2

'В (5 )У1

+ а2 52 +1

15 — р

(31)

(32)

Информация о порядке ведения функций Р5± (р), Рт — (р) содержится в работах [3-6]. Посредством формул (16)-(32) матричная задача факторизации (14) является решенной. Остается определить векторную функцию П (р), входящую в (15). Для этого, а также для получения асимптотических выражений для напряжений и смещений понадобятся асимптотические разложения интегралов (27), (32) и их комбинаций.

3. Разложение вблизи нуля

Вблизи р = 0 разложения имеют вид

Л+(р ) = Л— (р ) = 51 (а) =

9а2

л/2

р

1 I 41п2

1 + 1-+ 51 I р

я

1 л/2-4/12 — 9а2 ё 1п J + ( р )

1 +

41п2

я

+ 81 I р

о (р 2) + о (р4),

ёр

11 Тс 'п [лА^РЙ

р=0

я •

?

5

ехр [±р0+ (р )] = [1 ±б2 р 1 + О (р 2)

ехр [±р0— (р)] = ^ЩГ[1 ± б2р] + О (р2) (1+а252)—1/21п

5, (а) = Т1пР+(р>

ёр

р=0

= 1 ад ё я 0 ёя

А(5) а сШ 5 В (5) V 1+а252 +1

сЪ

5

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

5 (а) = 41п2 +51 (а) + 52 (а) = я

41п2 1 ад ё

-+ -1—

я я 0 ёэ

1п^ А (5 )В (5) + (1 + а2 52)—1/21п

А (5) ая сЛ ^

1В (5 + а2 52 +1

(33)

(34)

(35)

(36)

(37)

(38)

(39)

Подстановка результатов разложения (34), (37) в (16) дает (здесь выписаны только члены, дающие вклад в интересующие искомые величины)

(1/Тэ 0 ^

-i л/4

0

0

ia

Í е3г'л/4

л/4 - 3a 3ia

2

Р2 +

е / е va г- 3i<

—^(a + 2i5) -V3 —¡=

2>/6a 7 6V2 V V4-

0(1) 0(1)

у

3a2 у

Р3 + 0 (р 4).

(40)

4. Разложение вблизи минус бесконечности

По свойствам интегралов типа Коши (например, [6]) поведение Р0-(р) из (32) при больших р определяется как

Ро-(р) = -— + 0(Р-2), Р p 4 '

, . 1 м / 2 2ч -1/2 A (s) ,(a) = — J(1 + a 2s2) ln l-L! V ' 2л -i ' ]¡ B (s)

¡A (s) as cth s

,22

(41)

IB(s) VT+a2s2 +1

С учетом (17) для больших отрицательных p справедливы оценки

ch фр0_ (p) = cos q + О (p-2), sh фР0_ (p) = -i sin q + О (p-2).

Разложение второй формулы (16) по большим отрицательнымр с учетом (42) дает

2>/2 Í cos q sin q V1 -i^ - sin q cos qу Vi 1 у

X-(p ) = -ip

Va

0(р 5/2).

(42)

(43)

5. Определение векторной функции П (р)

Для определения векторной функции П (р) (15) существенным является вырожденность главного члена разложения матрицы X-(р) для больших р . Действительно, для вектора

Co =

Vi у

(44)

справедливо равенство

X-(р)С0 = О (р 5/2). (45)

Следующий член разложения X- (р) уже не вырожден. Из (13) следует, что функция

П ( Р )

не может иметь членов порядка р , причем старший член должен иметь вид

Ьр 1С0. Из требования ограниченности решения Б—(р) в нуле (12), следует наиболее общий вид для П (р) :

П ( р ) = Ао р—2 + Ь Со р—1. (46)

Здесь А0, Ь - векторная и скалярная константы, которые можно определить из условия (12), используя (15), (40):

0 л/2а

31 я/4 Iа^^

V —1Т У

Ь =

„1 я/4 /17 е у/а

2^2

12л/3 (М + 5N) + а >/3 ( — 1Т) + л/ 4 — 3а2 Т

(47)

(48)

6. Определение параметров разрушения

К параметрам разрушения относятся коэффициент интенсивности напряжений (КИН) и скорость высвобождения энергии при росте трещины.

Подстановка (46)-(48) и (41), (42) во второе уравнение (15) для р ^ —ад дает

Б- (р ) =

1

4——2р

кмМ + 1 кт +02км IТ + кКN

(49)

С использованием теоремы абелевого типа [14, 15] немедленно получаем асимптотику напряжений при х 0 :

о

- (х) =

1

V—2ях

кмМ + | кт км |Т + к^

Отсюда КИН

(50)

^1 = кмМ + 1 кт + укм | Т + к^.

(51)

Здесь введены обозначения

км = V6

I соб ад — б1п ад ^

соб ад + бш ад

, кт

У

42

\

>/4 — 3а2 I—соб ад — б1п ад ч соб ад — б1п ад у

, =5к аУэ к

кк =5км I-г кт

3а

2

Полагая

я/4 — ад = ю,

формулы (52) представим в виде, аналогичном [16, 19]:

км

^ соб ю^

V Б1п Юу

, кт

3а2

(— б1П ю^

V соб ю у

(52)

(53)

(54)

На основании обобщения формулы Ирвина [20, 21] следует, что скорость высвобождения энергии

Е(1) + Е(2)

J = (( + )■

2 E(1) E(2)

E(1) + E(2) 2 E(1) E(2)

12СM + 5N--2Г I +(4-3a2)| T +

a

V3

2

(55)

л/ 4 - 3a2

N

При отсутствии нормальной силы выражение (55) совпадает с выражениями, даваемыми в [10, 19]. Зависимости входящих в данные выражения величин 5 и ю от первого параметра Дундурса а представлены на рис. 2 и 3.

Рис. 2. Зависимость величины 5 от первого параметра Дундурса a Fig. 2. Dependence between value 5 and Dundurs first parameter a

Рис. 3. Зависимость величины ю от первого параметра Дундурса a Fig. 3. Dependence of value ю on Dundurs first parameter a

7. Разложение решения матричной задачи Римана вблизи нуля и вычисление параметров эффективной упругой заделки

Подстановка (46)-(48) и (16) в первое выражение (15) с последующим разложением по малым р дает

0

n=-3

F+(p)= Z G(n)p-n + O(p).

(56)

Здесь

—6 N Л

V 0 у

ч(-2) _

6М — 3аТЛ 3аN

С(-1) =

65М

л/12 — 9а2

Л г 3 г

— 35 аТ + 2 + 352 — 3 (р "0 + J "0) + - а2 1 + 1п-

3а

2 Л Л

4 — 3а2

N

уу

С(0) =

2Т — 3аМ

г £ШМ + ^ + gттT Л

^2 ММ + g 2 NN + g 2ТТ

glM =—2 + 352 +3 (ри0 + J -0) —4 а2

(

3а

2 Л — 1

1п -

V 4 — 3а у

(

g 2М а

л/12 — 9а2

— 35

g1T 3а(522 +Р "0 + J '' ) — 3 Г

—а"

2 2

Г ^ + 511 —(2 — 3а2 )5—^ а'УПа7, V я у 2

1п-

3а

2

V 4 — 3а2 у

g2т =(4 — 3а 2)

glN = 253 —

24 С(3)

+3а2

я

41п2

-5 — 2 я

— 51 I — (J'''0 + Р '"0),

g 2 N

а

л/3

л/4 — 3а2 5 +—| — 2 — 52 +Р ''0 + J ''0 | — 3

(

—а

1п

3а

2Л - +1

^^ 4 — 3а у

(58)

(59)

(60) (61) (62)

(63)

(64)

(65)

2 2 V 3

Здесь С,(р) - дзета функция Римана С(3)~ 1,2021; функции J''0, Р''0, J'''0, Р'''0 находятся по формулам прил. 2 работы [1]. Для использующихся в дальнейшем величин получаем

1 ад ё3 (, гтт^ттттт 1 ёэ

5

1''' = —

0

Р -0=—-ч

тт •>

1 ад ё3

я 0 ёэ

11 СС~3 }

(1 + а2 52)—1/21п

(66)

А (я) аэ йЬ я В (я+ а2я2 +1

5

(67)

Величина 51 определяется формулой (35).

Выражение (56) представляет собой трансформанту Лапласа от асимптотики производной половины скачка смещений. В соответствии со свойствами преобразования Лапласа (напр. [22]) трансформанта разности между смещениями на бесконечности и смешениями в нуле определяется делением образа на р :

0 [V(2) (х,0) — V(1) (х,0)1 2 (Е(1) + Е(2) К () 1 /ч

Н \е—рхёх = У (1) (2) ' ^ С(п)р—п—1 + О(р),

—1ад|и(2) (х,0) — и(1) (х,0)^ Е(1) Е(2) У

Е(1) Е

(68)

п=—3

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

\у(2) (х,0)- V(1) (х,0) [и(2) (х,0)- и(1) (х,0)

2 ( Е

(1)

Е

(2)

1

(1) (2)

Е(1) Е

| ^ С(п)р-п-1ерх

2л ьп=-з

'(1)-

(69)

Здесь Ь - контур, проходящий параллельно мнимой оси плоскости комплексной переменной р и расположенный между мнимой осью и ближайшей особой точкой. Этот контур может быть деформирован в соответствии с правилами интегрирования в комплексной плоскости. Находя оригинал по изображению, можно видеть, что скачок смещений согласно (57)-(69) представляет собой полином третей степени, соответствующий изгибу оси балки. Наличие членов нулевой и первой степени по х (минус первой и минус второй степени по р для трансформанты) свидетельствует о том, что граничные условия для эквивалентной балки не являются условиями жесткой заделки в месте сопряжения.

Члены разложения вертикального смещения V, содержащие С(_3), С(_2), и горизонтального смещения и, содержащие С(_1) (57), соответствуют членам для смещения балки-стержня под действием силы и момента. Оставшиеся члены соответствуют постоянным смещениям и повороту. Как следует из полученного решения, эти члены пропорциональны действующему главному вектору усилия и изгибающему моменту. Это позволяет выписать матрицу коэффициентов податливости [17, 18].

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

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

г и ^

V V у

2 (Е(1) + Е(2))

(1)77(2)

К

£2

а

а

>/12 - 9а2

- 35

а V

2

4 - 3а2 -V125)

Е(1) Е

>/12 - 9а2

65

-а2 + 652 2

К

Г т Л м

VР у

(70)

- 35

а

(V4 - 3а2 -Тш)2

-а2 + 652

ш

(71)

Зависимость коэффициентов упругой заделки К22, К11, К33 от параметра а, зависящего, в свою очередь, от отношения упругих модулей ^ (9), представлены на рис. 4, 5 и 6 соответственно. На этих же рисунках представлены зависимости тех же коэффициентов, по-

считанные в приближении упрощенных моделей [1, 2], построенных в предположении возможности пренебрежения перекрестными членами в (7), для различных величин толщины И, а для величины К22 - для случая полосы, отслаивающейся от полуплоскости [7].

■— --

:

'•v 4

0,2 0,4 0,6 0,8 1,0 а

Рис. 4. Зависимость коэффициента упругой заделки K22 от параметра а . Сплошная линия - полученное решение; точечная линия - приближенное решение [1] для h = 1; пунктирная линия (короткие штрихи) - то же для h = 2 ; пунктирная линия (длинные штрихи) - то же для h = 10 ;

штрихпунктирная линия - решение для слоя на полуплоскости [7] Fig. 4. Dependence of coefficient of elastic constraint K22 on parameter а . Solid graph shows the resulted solution; dotted line shows an approximate solution [1] for h = 1; dotted line (short dashes) is the same for h = 10 ; the dotted line (long dashes) is the same for ; dash-dot line shows the solution for the layer on the half-plane [7]

Рис. 5. Зависимость коэффициента упругой заделки K11 от параметра а . Сплошная линия - полученное решение; точечная линия - приближенное решение [2] для h = 1; пунктирная линия (короткие штрихи) - то же

для h = 2 ; пунктирная линия (длинные штрихи) - то же для h = 10 Fig. 5. Dependence of coefficient of elastic constraint K11 on parameter а . Solid graph shows the resulted solution; dotted line shows an approximate solution [2] for h = 1; dotted line (short dashes) is the same for h = 2 ; the dotted line (long dashes) is the same for h = 10

5 4 3 2 1

0,2 0,4 0,6 0,8 1,0 a

Рис. 6. Зависимость коэффициента упругой заделки K33 от параметра а .

Сплошная линия - полученное решение; точечная линия - приближенное решение [1] Fig.6. Dependence of coefficient of elastic constraint K33 on parameter а .

Solid graph shows the resulted solution; dotted line shows an approximate solution [1]

Из представленных графиков видно, в частности, что для случая расслоения полос одинаковой толщины (h = 1) модели [1, 2] дают точные решения для полос одинаковой жесткости и решения, близкие к точным для несильного отличия модулей (сплошная и точечные линии на рис. 4-6). С ростом различия в толщинах (h ^да) и приближении решения к случаю полосы на полуплоскости модели [1, 2] дают точные решения для весьма жестких полос (пунктирная линия с длинными штрихами и штрих-пунктирная линия на рис. 4). Отметим, что линия, соответствующая десятикратному различию в толщинах (h = 10) на рис. 4 неотличима от линии, соответствующей асимптотике (h ^да) модели [1] (формула (56), в которой необходимо исправить опечатку - убрать двойку перед ^). На рис. 5, 6 не представлены линии, соответствующие решению слоя на полуплоскости, как точного [7], так и приближенного [23], полученного в тех же предположениях, что и решения [1, 2] ввиду наличия в них логарифмической зависимости от расстояния до вершины трещины [7, 23, 24]. Полученные аналитические выражения для компонент матрицы упругой заделки представляются полезными для решения задач, связанных с отслоением покрытий [10, 18, 19, 25-27], с работой усиливающих ребер жесткости [28], с определением адгезионных свойств [29-32].

Заключение

Получено решение однородной задачи о полубесконечной трещине, проходящей вдоль интерфейса, разделяющего два упругих слоя одинаковой толщины и различных упругих свойств (однако на сочетание упругих констант наложено дополнительное ограничение, состоящее в равенстве нулю второго параметра Дундурса). Путем применения преобразования Лапласа задача сведена к матричной задаче Римана-Гильберта. Получены асимптотические выражения для смещений берегов трещины вдали от ее вершины. Показано, что ведущие члены асимптотики смещений берегов трещины соответствуют смещением балки (пластины) при граничных условиях типа обобщенной упругой заделки, т.е. условиях пропорциональности смещений и угла поворота в точке заделки действующим компонентам главных вектора и момента нагрузки. Получены выражения для компонент

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

Приложение

Решение матричной задачи Римана для рассматриваемого частного случая

Изложение следует работам [3-6]. Пусть а^ (/, ] = 1,2) являются компонентами невырожденной матрицы Л(/),

А( 2) =

а11 а12

а21 а22

(72)

определенной на мнимой оси Ь комплексной плоскости 2 . Задача состоит в нахождении такой кусочно-голоморфной матрицы-функции (2х2) X (г) в форме

Х( 2) =

(73)

Х+ (2), 2 е Я+, Х_ (2), 2 е Я -,

такой, что X + (2) и X_ (2), а также их обратные матрицы Х+-1(2) и Х -1(2) были бы голоморфными в Я+ и Я_, и на контуре Ь выполнялось бы условие

X + (Г)X_ (Г) = А (Г), I е Ь. (74)

В общем случае решение в замкнутой форме поставленной задачи неизвестно. Однако для рассматриваемого далее частного случая эффективное решение было получено [3-6].

Под девиатором матрицы ёеу А будем понимать матрицу, получаемую из исходной матрицы А вычитанием единичной матрицы I, умноженной на половину ее следа я :

ёеу А = А _ я I =

1

а11 а22

2

а,

а.

12

а22 а11

21

(75)

Я = ^ К + а22 ).

При этом, поскольку определитель девиатора

ёе1ёеу А = _ 4(а11 _ а22 )2 _

выполняется следующее тождество:

ёе! А = я2 + ёе! ёеу А. Данное тождество может быть переписано в виде

а12а21,

(сСе! А )1/2

(-с1е1ёеу А ) (сСе! А )1/2

= 1.

(76)

(77)

(78)

2

Определитель и показатель матрицы определяются следующим образом:

Д = ёй А = ^ 2, (79)

8 = ±1п (V ^), (80)

где , X2 - собственные числа матрицы А,

Хх,2 = 5 + (- ё^ёеу А )1/2. (81)

Отсюда получаются следующие соотношения:

А-1/2 5 = оЬ 8, (82)

А-1/2 (-ёе1ёеу А )1/2 = бЬ 8, (83)

Подстановка последних двух выражений в (75) дает

А = А

1/2

I оЬ 8 + (-ёе1ёеу А) ёеу А бЬ 8 , (84)

Определив коммутанту матрицы А как

Я = сот А = (- (Мёеу А )-1/2<ёеу А, (85)

перепишем выражение (84) как

А-1/2 А = [1 оЬ 8 + Я бЬ 8]. (86)

Укажем на два важные свойства коммутанты: во-первых, она не меняет значения при умножении девиатора на произвольный ненулевой множитель; во-вторых, ее квадрат дает единичную матрицу

Я2 = I. (87)

Пусть имеются две матрицы А1, А 2 с одинаковыми коммутантами Я1 = Я 2 = Я и, вообще говоря, различными определителями А1, А 2 и показателями 81, 8 2. Используя формулы (86) и (87), их произведение можно представить как

(А1А2 )-12 А1А2 = (А1А2 )-12 А2 А1 = [I оЬ 81 + Я бЬ 81 ][1 оЬ 82 + Я БЬ 82 ] =

= I оЬ 81 оЬ 82 + Я2 бЬ 81 бЬ 82 + Я[Ь 81 оЬ 82 + оЬ 81 бЬ 82 ] = (88)

= I оЬ ( + 82 ) + Я бЬ (81 + 82 ).

Из сравнения (86) и (88) следует, что показатель произведения матриц с одинаковыми коммутантами равен сумме показателей сомножителей. Факт равенства определителя произведения матриц произведению их определителей является общеизвестным.

Вычисление обратной матрицы в представлении (86) также весьма удобно. Непосредственная проверка легко убеждает, что

А 1 = А-1/2 [IоЬ8-ЯБЬ8]. (89)

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

АЛ-1 = А-1А = АА-1 =[1 еЬ 8 + Я бЬ е][1 еЬ в-Я бЬ е] = I еЬ2 в-Я2бЬ2 е= I. (90)

Рассмотрим поставленную задачу. Ограничимся случаем, когда матричный коэффициент А(/) имеет коммутанту К/) вида

R(t) = y:1/2(t)

l (t) m(t)

n(t) :l (t )

(91)

где

V (t) = 12(t) + m (t) n(t), (92)

а l(t), m(t), n(t) - полиномы.

Матрицу R(t) можно рассматривать как граничное значение R (z) вдоль мнимой оси, определенной для любых точек комплексной плоскости:

l(z) m(z)

R( z) = Г1/2( z)

(93)

(94)

п(г) -¡(I)

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

X+1 (г) = Л±1 (г) [I еЬ [(р+)(!)] ± К (г) бЬ [(Р+ ) (г)

Х-1 (г) = Л- (г)[1 еЬ[(^1/2Р-)г)]±К(г) бЬ[(у1/2р-))г)]

Здесь Л±( г), Р±( г) - функции, голоморфные справа и слева от мнимой оси. Голоморфность матриц-функций Х±1 (г) соответственно в областях £ + и доказывается, например, разложением гиперболических синуса и косинуса в степенные ряды, сходящиеся, как известно, в любых точках комплексной плоскости:

ад

еЬ [М)(г )] = !

Vj р2j

r1/2sh [(vP)(z )] = р£

j=о (2j)!'

Vj Р2 j

j=0 (2j +1)!'

Уравнения (16) удовлетворяются тождественно, если

Л:1 (z )Л+( z) = А1/2 (z), z е L, Р+(z):Р:(z) = V:1/2 (z)в(z), z е L.

(95)

(96)

(97)

(98)

После проведения факторизации матрицы в форме (94) решение задачи находится по теореме Лиувилля.

Работа выполнена при финансовой поддержке Программы Президиума РАН № 1.33.

Библиографический список

1. Устинов К.Б. О расслоении полосы по границе раздела упругих свойств. Часть 1. Постановка задачи, случай нормального отрыва // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2015. - № 4. - С. 226-245.

2. Устинов К.Б. О расслоении полосы по границе раздела упругих свойств. Часть 2. Случай сдвиговой трещины // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2016. - № 2. - С. 131-142. DOI: 10.15593/perm.mech/2016.2.09

3. Златин А.Н., Храпков A.A. Полубесконечная трещина, параллельная границе упругой полуплоскости // Докл. АН СССР. - 1986. - Т. 31. - С. 810-813.

4. Златин А.Н., Храпков A.A. Упругая полуплоскость, ослабленная трещиной, параллельной ее границе // ЛГУ. Исследования по упругости и пластичности. - 1990. - T. 16. Проблемы современной механики разрушения. - С. 68-75.

5. Златин А.Н., Храпков A.A. Векторная задача Римана с ненулевым индексом показателя матрицы-коэффициента // Изв. ВНИИГ им. Б.Е. Веденеева. - 1985. - Т. 181. - С. 12-16.

6. Khrapkov A.A. Wiener-Hopf method in mixed elasticity theory problems. - S.-P., 2001.

7. Устинов К.Б. Еще раз к задаче о полуплоскости, ослабленной полубесконечной трещиной, параллельной границе // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2013. - № 4. - C. 138-168.

8. Устинов К.Б.. Об отслоении слоя от полуплоскости; условия упругой заделки для пластины, эквивалентной слою // Изв. РАН. МТТ. - 2015. - № 1. - С. 75-95.

9. Dundurs J., Discussion: "Edge-bonded dissimilar orthogonal elastic wedges under normal and shear loading" (Bogy, D. B., 1968, ASME J. Appl. Mech., 35, pp. 460-466) // Trans ASME J. Appl. Mech. - Vol. 36. - Ко. 3. - P. 650-652.

10. Suo Z., Hutchinson J.W. Interface crack between two elastic layers // Int. J. Fract. - 1990. -Vol. 43. - P. 1-18.

11. Ентов В.М., Салганик Р. Л. О балочном приближении в теории трещин // Изв. АН СССР. Механика. - 1965. - № 5. - С. 95-102.

12. Fitcher W.B. The stress intensity factor for the double cantilever beam // Int. J. Fract. 1983. -Vol. 22. - P. 133-143.

13. Foote R.M.L., Buchwald V.T. An exact solution for the stress intensity factor for a double cantilever beam // Int. J. Fract. 1985. - Vol. 29. - P. 125-134.

14. Нобл Б. Применение метода Винера-Хопфа для решения дифференциальных уравнений в частных производных. - М.: Изд-во иностр. лит. 1962. - 279 с.

15. Гахов Ф.Д. Краевые задачи. - М.: Наука, 1977. - 640 с.

16. The edge cracking and spalling of brittle plates / M.D. Thouless, A.G. Evans, M.F. Ashby, J.W. Hutchinson // Acta Metallurgica. - 1997. - Vol. 35. - P. 1333-1341.

17. Ustinov K.B., Dyskin A.V., Germanovich L.N. Asymptotic analysis of extensive crack growth parallel to free boundary // 3rd Int. Conf. Localized Damage 94. - 1994. - P. 623-630.

18. Cotterell B., Chen Z. Buckling and cracking of thin film on compliant substrates under compression // Int. J. Frac. - 2000. - Vol. 104. - No 2. - P. 169-179.

19. Hutchinson, J.W., Suo, Z., Mixed Mode Cracking in Layered Materials // Advances in Applied Mechanics / eds J.W. Hutchinson and T.Y. Wu. - 1992. - Vol. 29. - P. 63-191.

20. Салганик Р.Л. О хрупком разрушении склеенных тел // ПММ. - 1963. - Т. 27, вып. 5. -C. 957-962.

21. Malyshev B.M., R.L. Salganik. The strength of adhesive joints using the theory of crack // Int. J. Fracture Mechanics. - 1965. - Vol. 1. - No. 2. - P. 114-128.

22. Дёч Г., Руководство к практическому применению преобразования Лапласа. - М.: Гос. изд-во физ.-мат. лит., 1958. - 207 с.

23. Устинов К.Б.. О сдвиговом отслоении тонкой полосы от полуплоскости // Изв. РАН. МТТ. - 2014. - № 6. - С. 141-152.

24. Салганик Р.Л., Устинов К.Б. Задача об упругозаделанной пластине, моделирующей частично отслоившееся от подложки покрытие (плоская деформация) // Изв. РАН МТТ. - 2012. -№ 4. - C. 50-62.

25. Yu H.-H., Hutchinson J.W. Influence of substrate compliance on buckling delamination of thin films // Int. J. Fract. - 2002. - Vol. 113. - P. 39-55.

26. Гольдштейн Р.В., Устинов К.Б., Ченцов А.В. Оценка влияния податливости подложки на напряжения потери устойчивости отслоившегося покрытия // Вычисл. мех. спл. сред. - 2011. -Т. 4, № 3. - С. 48-57.

27. Кургузов В.Д. Моделирование отслоения тонких пленок при сжатии // Выч. мех. спл. сред. - 2014. - Т. 7, № 1. - С. 91-99.

28. Попов Г.Я. Изгиб полубесконечной плиты, лежащей на линейно-деформируемом основании // ПММ. - 1961. - № 25, вып. 2. - С. 342-355.

29. Dannenberg H. Measurement of Adhesion by a Blister Method // J. Appl. Pol. Sci. - 1961. -Vol. 5. - No. 14. - P. 125-134.

30. Williams M.L. The continuum interpretation for fracture and adhesion // J. Appl. Pol. Sci. -1969. - Vol. 13. - P. 29-40.

31. Karihaloo B.L., Stang H. // Composites: Part B. - 2008. - Vol. 39. - P. 386-395.

32. Storakers B., Anderson B., Nonlinear plate theory applied to delamination in composites // J. Mech. Phys. Solids. - 1988. - Vol. 36. - P. 689-718.

References

1. Ustinov K.B. O rassloenii polosy po granice razdela uprugih svojstv. Chast' 1. Postanovka zadachi, sluchaj normal'nogo otryva [On delamination of a stipe along the boundary between two elastic layers. Part 1, problem formulation, the case of normal crack] Vestnik Permskogo nacional'nogo issledovatel'skogopolitekhnicheskogo universiteta. Mekhanika 2015, no. 4, рр. 226-245 [in Russian].

2. Ustinov K.B. On delamination of a stipe along the boundary between two elastic layers. Part 1, the case of shear crack. PNRPU Mechanics Bulletin, 2016, no. 2, pp. 131-142. DOI: 10.15593/perm.mech/2016.2.09

3. Zlatin A. N., Khrapkov A. A., A Semi-Infinite Crack Parallel to the Boundary of the Elastic Half-Plane, Sov. Phys. Dokl, 1986, Vol. 31, pp. 1009-1010.

4. Zlatin A. N., Khrapkov A. A., Uprugaya poluploskost, oslablennaya treshinoi, parallelnoi eyo granitse [Elastic Half-Plane Weakened by a Crack Parallel to Its Boundary], Issledovaniyapo uprugosti i plastichnosti [Studies in Elasticity and Plasticity], vol. 16: problem sovremennoi mechaniki razrusheniya [Problems of Contemporary Fracture Mechanics], (LGU, 1990), pp. 68-75 [in Russian].

5. Zlatin A. N., Khrapkov A. A., Vectornaya zadacha Rimana s nenulevym indeksom pokazatelya matritsy-koeffitsienta [Vector Riemann Problem with Zero Index of Matrix-Coefficient Exponent] Izv. VNIIG im. B. E. Vedeneeva, 1985, vol. 181, pp. 12-16.

6. Khrapkov A.A. Wiener-Hopf method in mixed elasticity theory problems. S.-P. 2001.

7. Ustinov K.B. Eshe raz k zadache o poluploskosti, oslablennoi polubeskonechnoi treshinoi, parallelnoi granitse [Once more to the problem on a half-plane, weakened by a semi-infinite crack, parallel to its boundary] VestnikPNIPUMechanika, 2013, no. 4. pp. 138-168 [In Russian]

8. Ustinov K.B., On separation of a layer from the half-plane: elastic fixation conditions for a plate equivalent to the layer. Mech. Solids, 2015, vol. 50, no. 1, pp. 62-80.

9. Dundurs J., Discussion: "Edge-bonded dissimilar orthogonal elastic wedges under normal and shear loading" (Bogy, D. B., 1968, ASME J. Appl. Mech., 35, pp. 460-466) Trans ASME J. Appl. Mech. vol. 36, по. 3, pp. 650-652.

10. Suo Z., Hutchinson J.W. Interface crack between two elastic layers. Int. J. Fract. 1990, vol. 43, pp. 1-18.

11. Entov V.M., Salganik R.L. O balochnom priblizhenii v teorii treshin [On beam approximation in crack theory/. Izv. ANSSSR. Mehanika, 1965, no. 5. pp. 95-102 [In Russian].

12. Fitcher W.B. The stress intensity factor for the double cantilever beam. Int. J. Fract, 1983, vol. 22, pp. 133-143.

13. Foote R.M.L., Buchwald V.T. An exact solution for the stress intensity factor for a double cantilever beam. Int. J. Fract, 1985, vol. 29, pp. 125-134.

14. Noble B., Methods Based on the Wiener-Hopf Technique for the Solution of Partial Differential Equations (Pergamon Press, London etc., 1958; Mir, Moscow, 1962) 279 p.

15. Gakhov F.D. Boundary Value Problems. (International Series of Monographs in Pure and Applied Mathematics, Volume 85) Oxford/London/Edinburgh/New York/Paris/Frankfurt 1966. Pergamon Press 564 p.

16. Thouless M.D., Evans A.G., Ashby M.F., Hutchinson, J.W., The edge cracking and spalling of brittle plates. Acta Metallurgica, 1997, vol. 35, pp. 1333-1341.

17. Ustinov K.B., Dyskin A.V., Germanovich L.N. Asymptotic analysis of extensive crack growth parallel to free boundary. 3rd Int. Conf. Localized Damage 94, 1994, pp. 623-630.

18. Cotterell B., Chen Z. Buckling and cracking of thin film on compliant substrates under compression. Int. J. Frac, 2000, vol. 104, no. 2, pp. 169-179.

19. Hutchinson, J.W., Suo, Z., Mixed Mode Cracking in Layered Materials. Advances in Applied Mechanics edited by J.W. Hutchinson and T.Y. Wu, 1992, vol. 29, pp. 63-191.

20. Salganik R. L., The Brittle Fracture of Cemented Bodies. J. Appl. Math. Mech, 1963, vol. 27, no.5, pp.1468-1478.

21. Malyshev B.M., R.L. Salganik. The strength of adhesive joints using the theory of crack. Int. J. Fracture Mechanics, 1965, vol. 1, no. 2, pp. 114-128.

22. Doetsch G., Rukovodstvo k prakticheskomu primeneniiu preobrazovaniia Laplasa [Handbook der Laplace-Transformation]. Birkhauser, Basel, 1946; Fizmatlit, Moscow, 1958.

23. Ustinov K.B., On Shear Separation of a Thin Strip from the Half-Plane. Mech. Solids, 2014, vol. 49, no. 6, pp. 713-724.

24. Salganik R.L., Ustinov K.B. Deformation problem for an elastically fixed plate modeling a coating partially delaminated from the substrate (plane strain). Mech. of Solids, 2012, vol. 47, no. 4, pp.415-425

25. Yu H.-H., Hutchinson J.W. Influence of substrate compliance on buckling delamination of thin films. Int. J. Fract, 2002, vol. 113, pp. 39-55.

26. Goldstein R.V., Ustinov K.B., Chentsov A.V. Otsenka vliianiia podatlivosti podlozhki na napriazheniia poteri ustoichivosti otsloivshegosia pokrytiia [Influence of substrate compliance on stresses generating stability loss of a delaminated coating]. VychislitelnayaMechanica Sploshnykh Sred, 2011, vol. 4, no. 3, pp. 48-57.

27. Kurguzov V.D. Modelirovanie otsloeniyz tonkih plenok pri szhatii [Imulation of thin films de-lamination under compression]. Vychislitelnaya Mechanica Sploshnykh Sred, 2014, vol. 7, no. 1, pp. 91-99.

28. Popov G.Ya. Bending of a semi-infinite plate resting on a linearly deformable foundation J Appl Math Mech, 1961, vol. 25, is. 2, pp. 502-520.

29. Dannenberg H. Measurement of Adhesion by a Blister Method, J. Appl. Pol. Sci., 1961, vol. 5, no. 14, pp.125-134.

30. Williams M.L. The continuum interpretation for fracture and adhesion, J. Appl. Pol. Sci., 1969, vol. 13, pp. 29-40.

31. Karihaloo B.L., Stang H., Composites: Part B, 2008, vol. 39, pp. 386-395.

32. Storakers B., Anderson B., Nonlinear plate theory applied to delamination in composites. J. Mech. Phys. Solids, 1988, vol. 36, pp. 689-718.