Full Text

1. Постановка задачи. Исходная конфигурация В первой части [1] была сформулирована задача теории упругости в постановке плоской деформации о полубесконечной трещине, проходящей вдоль границы, разделяющей два изотропных упругих слоя с произвольными соотношениями толщин и упругих постоянных. Самоуравновешенная нагрузка, характеризуемая главным вектором и изгибающим моментом, предполагалась приложенной на бесконечности со стороны трещины. Путем применения преобразования Лапласа задача была сведена к задаче Римана-Гильберта с матричным коэффициентом. Далее приводится изложение точного решения поставленной задачи для частного случая: слоев одинаковой толщины, на соотношение упругих констант которых наложено некоторое ограничение. Исходная конфигурация задачи представлена на рис. 1. Рис. 1. Расслоение составного слоя. Геометрия задачи Fig. 1. Delamination of contained layer. Geometry of the problem Здесь и далее - компоненты тензора напряжений, - компоненты вектора смещения; все величины, относящиеся к нижней полосе обозначаются индексом 1, все величины, относящиеся к верхней полосе - индексом 2. Для рассматриваемой однородной задачи предполагается, что все поверхности свободны от напряжений, при и при (1) а нагрузка с эквивалентными главным вектором и моментом приложена на бесконечности, так что (здесь по сравнению с [1] изменены знаки) (2) Условия сопряжения на границе имеют вид при (3) Модифицированные для условий плоской деформации модули Юнга и коэффициенты Пуассона материалов обозначим соответственно. С обычными модулями Юнга и коэффициентами Пуассона они связаны соотношениями (4) Путем применения преобразования Лапласа [1] задача сведена к задаче Римана-Гильберта с матричным коэффициентом. Запишем ее для условий одинаковых толщин полосы, равных единице и дополнительного ограничения на сочетание упругих слоев, необходимого для получения аналитического решения, а именно (5) Данное условие соответствует равенству нулю второго параметра Дундурса (параметра упругого несоответствия) [9, 10]; при выполнении данного условия при сингулярности отсутствует осцилляция. Итак, матричную задачу Римана-Гильберта на мнимой оси для указанных ограничений запишем в виде (без нарушения общности, будем полагать - данное условие позволит избежать необходимости рассмотрения различных случаев выделения ветвей функций, содержащих радикалы [5, 6]) (6) (7) (8) (9) (10) (11) с условиями в опорных точках (в нуле и на бесконечности) (12) (13) Значение соответствует расслоению однородной полосы [11-13]. Основная сложность состоит в факторизации матричного коэффициента , т.е. представлении его в виде (14) где функции аналитичны в правой и левой полуплоскости комплексного переменного , и в соответствующих плоскостях вплоть до границы. В настоящее время общее точное аналитическое решение задачи для произвольной матрицы неизвестно, однако для указанных ограничений задача решается применением метода [3-6] (см. также литературу, указанную в [1]). После нахождения окончательное решение задачи дается с помощью теоремы Лиувилля: (15) Здесь - векторная функция, подлежащая определению. 2. Решение задачи Римана-Гильберта Матрицы могут быть представлены в виде [3-6] (16) Здесь - единичная матрица, (17) Функции являются решениями скалярных задач Римана-Гильберта (18) (19) Здесь определитель и показатель матрицы определяются как (20) (21) Собственные числа матрицы могут быть представлены в виде (22) (23) Таким образом, исследуемая задача отличается от задач, решенных в [7, 8], видом функций , определяемых формулами (23). С учетом (20)-(22) скалярные задачи Римана-Гильберта могут быть записаны в виде (24) (25) Решение задачи (18) или (24) запишем в виде, аналогичном данному в работах [3-6] через интегралы типа Коши: (26) (27) внеинтегральные члены получены факторизацией тангенса через гамма-функцию. Решение задачи (19) также запишем в виде, аналогичном данному в работах [3-6]: (28) Вспомогательные функции (29) (30) введены таким образом, чтобы правая часть выражения , (31) удовлетворяла условиям Гёльдера. Тогда решение задачи [3-6]. (32) Информация о порядке ведения функций содержится в работах [3-6]. Посредством формул (16)-(32) матричная задача факторизации (14) является решенной. Остается определить векторную функцию , входящую в (15). Для этого, а также для получения асимптотических выражений для напряжений и смещений понадобятся асимптотические разложения интегралов (27), (32) и их комбинаций. 3. Разложение вблизи нуля Вблизи разложения имеют вид (33) (34) (35) (36) (37) (38) Для удобства проведения дальнейших преобразований введем обозначение (39) Подстановка результатов разложения (34), (37) в (16) дает (здесь выписаны только члены, дающие вклад в интересующие искомые величины) (40) 4. Разложение вблизи минус бесконечности По свойствам интегралов типа Коши (например, [6]) поведение из (32) при больших определяется как (41) С учетом (17) для больших отрицательных справедливы оценки (42) Разложение второй формулы (16) по большим отрицательным р с учетом (42) дает (43) 5. Определение векторной функции Для определения векторной функции (15) существенным является вырожденность главного члена разложения матрицы для больших . Действительно, для вектора (44) справедливо равенство (45) Следующий член разложения уже не вырожден. Из (13) следует, что функция не может иметь членов порядка , причем старший член должен иметь вид . Из требования ограниченности решения в нуле (12), следует наиболее общий вид для (46) Здесь - векторная и скалярная константы, которые можно определить из условия (12), используя (15), (40): (47) (48) 6. Определение параметров разрушения К параметрам разрушения относятся коэффициент интенсивности напряжений (КИН) и скорость высвобождения энергии при росте трещины. Подстановка (46)-(48) и (41), (42) во второе уравнение (15) для дает (49) С использованием теоремы абелевого типа [14, 15] немедленно получаем асимптотику напряжений при (50) Отсюда КИН (51) Здесь введены обозначения , (52) Полагая (53) формулы (52) представим в виде, аналогичном [16, 19]: (54) На основании обобщения формулы Ирвина [20, 21] следует, что скорость высвобождения энергии (55) При отсутствии нормальной силы выражение (55) совпадает с выражениями, даваемыми в [10, 19]. Зависимости входящих в данные выражения величин и от первого параметра Дундурса представлены на рис. 2 и 3. Рис. 2. Зависимость величины от первого параметра Дундурса Fig. 2. Dependence between value and Dundurs first parameter Рис. 3. Зависимость величины ω от первого параметра Дундурса Fig. 3. Dependence of value ω on Dundurs first parameter 7. Разложение решения матричной задачи Римана вблизи нуля и вычисление параметров эффективной упругой заделки Подстановка (46)-(48) и (16) в первое выражение (15) с последующим разложением по малым дает (56) Здесь , (57) (58) (59) (60) , (61) (62) (63) (64) (65) Здесь - дзета функция Римана ; функции находятся по формулам прил. 2 работы [1]. Для использующихся в дальнейшем величин получаем (66) (67) Величина определяется формулой (35). Выражение (56) представляет собой трансформанту Лапласа от асимптотики производной половины скачка смещений. В соответствии со свойствами преобразования Лапласа (напр. [22]) трансформанта разности между смещениями на бесконечности и смешениями в нуле определяется делением образа на (68) Асимптотики смещений вычисляются отсюда с помощью обратного преобразования Лапласа (69) Здесь - контур, проходящий параллельно мнимой оси плоскости комплексной переменной и расположенный между мнимой осью и ближайшей особой точкой. Этот контур может быть деформирован в соответствии с правилами интегрирования в комплексной плоскости. Находя оригинал по изображению, можно видеть, что скачок смещений согласно (57)-(69) представляет собой полином третей степени, соответствующий изгибу оси балки. Наличие членов нулевой и первой степени по (минус первой и минус второй степени по для трансформанты) свидетельствует о том, что граничные условия для эквивалентной балки не являются условиями жесткой заделки в месте сопряжения. Члены разложения вертикального смещения , содержащие , и горизонтального смещения , содержащие (57), соответствуют членам для смещения балки-стержня под действием силы и момента. Оставшиеся члены соответствуют постоянным смещениям и повороту. Как следует из полученного решения, эти члены пропорциональны действующему главному вектору усилия и изгибающему моменту. Это позволяет выписать матрицу коэффициентов податливости [17, 18]. Здесь следует помнить, что формулы (56)-(69) соответствуют смещениям нижней границы полосы (балки), горизонтальное же смещение центральной оси балки отличается на величину произведения поворота в данной точке на половину толщины полосы. Также следует учесть, что продольные силы приложены к краям полос, а не к центрам, как в балочных теориях. По-видимому, для матрицы коэффициентов податливости следует использовать симметризованную матрицу полученных коэффициентов [19]. Окончательно, выражение для матрицы коэффициентов упругой заделки записывается следующим образом: (70) (71) Зависимость коэффициентов упругой заделки от параметра , зависящего, в свою очередь, от отношения упругих модулей (9), представлены на рис. 4, 5 и 6 соответственно. На этих же рисунках представлены зависимости тех же коэффициентов, посчитанные в приближении упрощенных моделей [1, 2], построенных в предположении возможности пренебрежения перекрестными членами в (7), для различных величин толщины а для величины - для случая полосы, отслаивающейся от полуплоскости [7]. Рис. 4. Зависимость коэффициента упругой заделки от параметра . Сплошная линия - полученное решение; точечная линия - приближенное решение [1] для ; пунктирная линия (короткие штрихи) - то же для ; пунктирная линия (длинные штрихи) - то же для ; штрихпунктирная линия - решение для слоя на полуплоскости [7] Fig. 4. Dependence of coefficient of elastic constraint on parameter . Solid graph shows the resulted solution; dotted line shows an approximate solution [1] for ; dotted line (short dashes) is the same for ; the dotted line (long dashes) is the same for ; dash-dot line shows the solution for the layer on the half-plane [7] Рис. 5. Зависимость коэффициента упругой заделки от параметра . Сплошная линия - полученное решение; точечная линия - приближенное решение [2] для ; пунктирная линия (короткие штрихи) - то же для ; пунктирная линия (длинные штрихи) - то же для Fig. 5. Dependence of coefficient of elastic constraint on parameter . Solid graph shows the resulted solution; dotted line shows an approximate solution [2] for ; dotted line (short dashes) is the same for ; the dotted line (long dashes) is the same for Рис. 6. Зависимость коэффициента упругой заделки от параметра . Сплошная линия - полученное решение; точечная линия - приближенное решение [1] Fig.6. Dependence of coefficient of elastic constraint on parameter . Solid graph shows the resulted solution; dotted line shows an approximate solution [1] Из представленных графиков видно, в частности, что для случая расслоения полос одинаковой толщины модели [1, 2] дают точные решения для полос одинаковой жесткости и решения, близкие к точным для несильного отличия модулей (сплошная и точечные линии на рис. 4-6). С ростом различия в толщинах и приближении решения к случаю полосы на полуплоскости модели [1, 2] дают точные решения для весьма жестких полос (пунктирная линия с длинными штрихами и штрих-пунктирная линия на рис. 4). Отметим, что линия, соответствующая десятикратному различию в толщинах на рис. 4 неотличима от линии, соответствующей асимптотике модели [1] (формула (56), в которой необходимо исправить опечатку - убрать двойку перед ). На рис. 5, 6 не представлены линии, соответствующие решению слоя на полуплоскости, как точного [7], так и приближенного [23], полученного в тех же предположениях, что и решения [1, 2] ввиду наличия в них логарифмической зависимости от расстояния до вершины трещины [7, 23, 24]. Полученные аналитические выражения для компонент матрицы упругой заделки представляются полезными для решения задач, связанных с отслоением покрытий [10, 18, 19, 25-27], с работой усиливающих ребер жесткости [28], с определением адгезионных свойств [29-32]. Заключение Получено решение однородной задачи о полубесконечной трещине, проходящей вдоль интерфейса, разделяющего два упругих слоя одинаковой толщины и различных упругих свойств (однако на сочетание упругих констант наложено дополнительное ограничение, состоящее в равенстве нулю второго параметра Дундурса). Путем применения преобразования Лапласа задача сведена к матричной задаче Римана-Гильберта. Получены асимптотические выражения для смещений берегов трещины вдали от ее вершины. Показано, что ведущие члены асимптотики смещений берегов трещины соответствуют смещением балки (пластины) при граничных условиях типа обобщенной упругой заделки, т.е. условиях пропорциональности смещений и угла поворота в точке заделки действующим компонентам главных вектора и момента нагрузки. Получены выражения для компонент матрицы коэффициентов упругой заделки. Также получены выражения для коэффициентов интенсивности напряжений и скорости высвобождения энергии. Полученное решение представляется полезным для решения задач об отслоении покрытий. Приложение Решение матричной задачи Римана для рассматриваемого частного случая Изложение следует работам [3-6]. Пусть являются компонентами невырожденной матрицы , (72) определенной на мнимой оси комплексной плоскости . Задача состоит в нахождении такой кусочно-голоморфной матрицы-функции (2х2) в форме (73) такой, что и , а также их обратные матрицы и были бы голоморфными в и , и на контуре выполнялось бы условие (74) В общем случае решение в замкнутой форме поставленной задачи неизвестно. Однако для рассматриваемого далее частного случая эффективное решение было получено [3-6]. Под девиатором матрицы будем понимать матрицу, получаемую из исходной матрицы вычитанием единичной матрицы , умноженной на половину ее следа (75) При этом, поскольку определитель девиатора (76) выполняется следующее тождество: (77) Данное тождество может быть переписано в виде (78) Определитель и показатель матрицы определяются следующим образом: (79) (80) где - собственные числа матрицы (81) Отсюда получаются следующие соотношения: (82) (83) Подстановка последних двух выражений в (75) дает (84) Определив коммутанту матрицы как (85) перепишем выражение (84) как (86) Укажем на два важные свойства коммутанты: во-первых, она не меняет значения при умножении девиатора на произвольный ненулевой множитель; во-вторых, ее квадрат дает единичную матрицу (87) Пусть имеются две матрицы с одинаковыми коммутантами и, вообще говоря, различными определителями и показателями . Используя формулы (86) и (87), их произведение можно представить как (88) Из сравнения (86) и (88) следует, что показатель произведения матриц с одинаковыми коммутантами равен сумме показателей сомножителей. Факт равенства определителя произведения матриц произведению их определителей является общеизвестным. Вычисление обратной матрицы в представлении (86) также весьма удобно. Непосредственная проверка легко убеждает, что (89) Действительно, (90) Рассмотрим поставленную задачу. Ограничимся случаем, когда матричный коэффициент имеет коммутанту вида (91) где (92) а - полиномы. Матрицу можно рассматривать как граничное значение вдоль мнимой оси, определенной для любых точек комплексной плоскости: (93) Будем искать решение поставленной задачи в виде (94) Здесь - функции, голоморфные справа и слева от мнимой оси. Голоморфность матриц-функций соответственно в областях и доказывается, например, разложением гиперболических синуса и косинуса в степенные ряды, сходящиеся, как известно, в любых точках комплексной плоскости: (95) (96) Уравнения (16) удовлетворяются тождественно, если , (97) (98) После проведения факторизации матрицы в форме (94) решение задачи находится по теореме Лиувилля.

About the authors

K B Ustinov

A.Yu. Ishlinsky Institute for Problems in Mechanics RAS; Moscow Aviation Institute

Email: ustinov@ipmnet.ru

References

Statistics

Views

Abstract - 1098

PDF (Russian) - 49