Full Text

Введение Задача продольного расслоения полосы, составленной из двух полос меньшей толщины, обладающих различными упругими свойствами, решалась многими авторами (например, [1-5]), при этом в основном интерес был сосредоточен на поведении решения вблизи вершины трещины. Однако в ряде задач, например при исследовании потери устойчивости отслаивающимися покрытиями, значение имеет поведение поля смещения вдали от вершины трещины [6-18]. Одной из наиболее простых конфигурация является полоса, составленная из двух полос меньшей толщины. Если длина отслоения много больше толщины полосы, для ряда случаев отслоение может рассматриваться как полубесконечное. Если толщина одного из слоев много больше протяженности отслоения, то такой слой может рассматриваться как полуплоскость; эффективное и элегантное решение задачи для такого случая, полученное путем применения преобразования Лапласа и сведения ее к матричной задаче Римана, было дано в работах [19-22] для одинаковых упругих свойств слоя и полуплоскости и в работе [17] для различающихся свойств. Ключевым моментом решения являлась факторизация матричного коэффициента. В указанных работах получены асимптотические представления для трансформант поля напряжений на линии продолжения трещины, и в конечном счете вычислены коэффициенты интенсивности напряжений (КИН) в зависимости от главных вектора и момента напряжений, действующих на продолжении трещины. Однако интерес к задаче не ограничивается вычислением КИН. В частности, из указанного решения можно извлечь асимптотику смещений берегов трещины вдали от ее вершины [14]. При этом при наличии ненулевых главных вектора и момента действующих нагрузок ведущие члены асимптотики смещений берегов трещины будут, очевидно, соответствовать смещения балки (пластины), подверженной действию указанных главных вектора и момента при некоторых граничных условиях. Следует подчеркнуть, что эти граничные условия вовсе не являются условиями жесткой заделки. Они представляют собой [7, 10-15, 17] условия обобщенной упругой заделки, т.е. условия пропорциональности смещений и угла поворота в точке заделки действующим компонентам главных вектора и момента нагрузки. В работе [17] данное решение было обобщено на случай отслоения и основания, обладающих различными упругими свойствами, хотя и связанных некоторыми дополнительными ограничениями. В работе [16] рассматривалась возможность пренебрежения влиянием нормальных к линии трещины напряжений на сдвиговые компоненты смещения (и сдвиговых компонент напряжения на нормальные компоненты смещения), в результате было получено приближенное решение о сдвиговой трещине на границе раздела. Там же было показано, что подобное решение по крайней мере не хуже традиционного подхода - замены отслаивающейся части полосы одномерным объектом (стержнем). В настоящей статье рассматривается аналогичная задача для расслаиваемой составной полосы. 1. Общая постановка задачи о расслоении полосы, состоящей из двух упругих слоев по границе контакта В условиях плоской деформации рассматривается неоднородная упругая полоса , составленная из двух изотропных частей с различными упругими свойствами, контакт которых проходит вдоль линии . Вдоль части границы имеется полный контакт, вдоль оставшейся части границы контакт отсутствует (рис. 1). Все величины, относящиеся к нижней полосе , будут обозначаться индексом 1, все величины, относящиеся к верхней полосе , индексом 2. Модули Юнга и коэффициенты Пуассона (модифицированные для условий плоской деформации) материалов полуплоскости и полосы обозначим соответственно. С обычными модулями Юнга и коэффициентами Пуассона они связаны соотношениями (1) В постановке плоской деформации рассматривается однородная задача: предполагается, что все поверхности свободны от напряжений, при и при (2) а нагрузка с эквивалентными главным вектором и моментом приложена на бесконечности, так что (3) Условия сопряжения на границе имеют вид при (4) Здесь - компоненты вектора смещения; - компоненты тензора напряжений. Рис. 1. Расслоение составного слоя. Геометрия задачи В окрестности нуля поле напряжений должно быть интегрируемо. (5) Из механики разрушения известно, что это поле имеет корневую особенность: (6) Однако для построения решения можно довольствоваться условием (5), а условие (6) получится из решения. В случае различных значений упругих постоянных здесь могут присутствовать осциллирующие члены. Воспользуемся результатами работы [22], в которой предлагается связь между образами двустороннего преобразования Лапласа от производной компонент вектора смещения (7) и напряжениями (8) для полосы (при поставленных граничных условиях свободной внешней границы при ). Вывод данных выражений дан в Приложении 1. (9) (10) Здесь - мнимая ось. Рассмотрим преобразование Лапласа от следующих величин: (11) (12) Здесь использовано свойство равенства нулю скачка смещения для отрицательной полуоси и равенства нулю напряжений для положительной полуоси. На основании (4) подынтегральное выражение (11) тождественно равно нулю для и, следовательно, аналитична в правой полуплоскости Аналогично согласно (2) подынтегральное выражение (12) равно нулю для , и, следовательно, аналитична в левой полуплоскости Здесь аналогично [19-22] предполагается выполнение условия убывания напряжений на бесконечности при . В окрестности нуля выполняются условия (6). При смещения могут расти как полином третей степени (что соответствует ненулевой поперечной силе на бесконечности), следовательно, в Лаплас-образе возможно появление полюса в нуле до третьего порядка включительно. Комбинируя (7)-(12), получаем уравнение матричной задачи Римана, аналогичное [19-22]: (13) где (14) (15) (16) (17) (18) Основная сложность состоит в факторизации матричного коэффициента , т.е. в представлении его в виде (19) где функции аналитичны в правой и левой полуплоскости комплексного переменного р. в соответствующих плоскостях вплоть до границы. В настоящее время общее точное аналитическое решение указанной задачи неизвестно; можно указать несколько путей построения решения: - применение аппроксимаций для матричного коэффициента (14) задачи методами [23-26]; - построение приближенных решений в предположении возможностью пренебречь перекрестными членами в (14) (положить ) [16]; - построение решения для частных случаев, для которых факторизация может быть осуществлена [27-34]. Далее решение будет получено с использованием второго из указанных методов. 2. Задача о трещине нормального отрыва. Постановка С помощью предположения о возможности пренебрежения в (14) влиянием перекрестных членов, что соответствует , матричная задача распадается на две скалярные. Данное приближение заведомо не хуже традиционной замены слоя однородным объектом - стержнем либо балкой [35-37]. Рассмотрим задачу о нормальном отрыве. , (20) (21) (22) (23) Аналог (19) становится скалярным уравнением (24) После нахождения окончательное решение дается с помощью теоремы Лиувилля (25) Здесь - функция, подлежащая определению. Второе из уравнений (3) и условия (6) для интересующих нас величин после трансформации принимают соответственно вид (26) (27) Из элементарных соображений также следует, что производные от нормальных смещения для нижней кромки верхней полосы при будут расти как , а нижней - как , что в Лаплас-образах согласно (17) даст (28) Условие отсутствия скачка смещений в нуле в трансформантах записывается следующим образом (29) Данные условия будут использованы для определения 3. Решение задачи Римана Уравнения (23) могут быть записаны в виде (здесь - мнимая ось) , (30) , (31) С помощью стандартного приема факторизации арктангенса через гамма-функцию решение может быть выражено следующим образом: (32) (33) (34) Для получения решения осталось определить функцию число независимых параметров в которой определяется индексом задачи и поведением решения в опорных точках. Наличие опорных точек только в нуле и на бесконечности определяет вид как сумму слагаемых с положительными и отрицательными степенями р. Более детальное исследование показывает, что наиболее общим видом является (35) Здесь - постоянные, определяемые из граничных условий. Наличие членов со степенями большими минус двух привело бы к нарушению условия (29), а меньшими минус трех - к нарушению условия (28). Для определения постоянных потребуется разложение функции вблизи нуля (два члена). Однако для определения параметров упругой заделки необходимо четыре члена. Для внеинтегрального члена (32) имеем (36) Здесь - дзета-функция Римана, Для интегрального члена (32) из (34), (97), (101), (103), положив , имеем (37) (38) (39) (40) (41) (42) Выражение для интеграла (40) можно преобразовать следующим образом: (43) из (40), (38), (43) следует (44) (45) (46) Подстановка выражения для функции (35) и разложения вблизи нуля функции (результата подстановки (36)-(39) в (32)) в первое из равенств (25) дает разложение Лаплас-образа производной от скачка смещения вблизи нуля. С точностью до двух ведущих членов оно имеет вид (47) Сравнение (47) с (28) позволяет вычислить значение констант, входящих в (48) Таким образом, задача решена: подстановка значений констант (48) в (35), а затем вместе с (32)-(34) в (25) дает Лаплас-образы напряжений и производной от разности смещений вдоль границы раздела упругих свойств. 4. Определение асимптотики смещения берегов разреза вдали от вершины трещины Лаплас-образ указанной асимптотики может быть получен путем подстановки (48), (35), (32), (36)-(39) в (25). Результат данной подстановки может быть записан в виде (49) (50) (51) (52) (53) Первый член в (49) соответствует прогибу балки, нагружаемой силой и моментом. Оставшиеся члены могут интерпретироваться как поправки к элементарной балочной теории (типа модели балки Тимошенко, учитывающей вклад в смещения от поперечной силы) и дополнительные смещения, вызываемые граничными условиями, отличающимися от условий жесткой заделки. Величины, определяемые формулами (50)-(53), являются при этом коэффициентами пропорциональности между поперечной силой и моментом, с одной стороны, и углом поворота и нормальным смещением - с другой. Поскольку из приведенного решения непосредственно невозможно выделить указанные причины и при этом естественно ожидать, что матрица податливости должна быть симметрична, для модели эквивалентной балки естественно принять граничное условие в симметричном виде: (54) Зависимость от для различных отношений модулей представлена на рис. 2. 2 4 6 8 10 h 10 20 30 40 K u'M Рис. 2. Зависимость коэффициента упругой заделки от относительной толщины основания для различных отношений модулей: сплошная линия - ; пунктирная - ; точечная - ; штрих-пунктирная - Для толстых оснований справедлива асимптотика (55) (56) Здесь по-прежнему определяется (45). Зависимость от отношения модулей для больших представлена на рис. 3. Данная зависимость неплохо аппроксимируется формулой (57) η 50 100 150 200 10 15 20 25 30 35 40 45 K u'М Рис. 3. Зависимость коэффициента упругой заделки от отношения модулей для больших толщин подложек 5. Определение параметров разрушения При для (33) справедливо следующее асимптотическое разложение: (58) Подстановка (58), (35), (48) в (25) и разложение для больших отрицательных дает (59) Последнее выражение немедленно дает асимптотику напряжений при (60) Отсюда КИН (61) Скорость высвобождения энергии [38, 39] при этом имеет вид (62) Полученное выражение для скорости высвобождения энергии для совпадает с элементарным решением (например, [2]). При рассмотрении задачи о распространении интерфейсной трещины в классической постановке приведенный анализ для определения параметров разрушения может рассматриваться лишь как приближенный, поскольку в общем случае в точном решении будут присутствовать обе моды. Однако задача имеет и самостоятельное значение [13] для случая, когда прилегающие слои могут скользить друг относительно друга, но удерживаются силами адгезии (например, вандерваальсовыми). Заключение Получено и исследовано решение однородной задачи о полубесконечной трещине, проходящей вдоль границы, разделяющей два упругих слоя с отличающимися свойствами. Путем применения преобразования Лапласа задача сведена к однородной матричной задаче Римана. В предположении возможности пренебрежения влиянием нормальных напряжений на сдвиговые смещения и сдвиговых напряжений на нормальные смещения задача сведена к двум скалярным задачам Римана. Получена асимптотика поля напряжений вблизи вершины трещины, имеющего корневую особенность при отсутствии осциллирующих членов. Посчитан коэффициент интенсивности напряжений. Показано, что ведущие члены асимптотики смещений берегов трещины вдали от вершины трещины соответствуют смещению балки при граничных условиях типа упругой заделки, то есть условиях пропорциональности смещения в точке заделки вектору усилия и изгибающему моменту. Приложение 1. Вывод связи между производной от скачка смещения и напряжениями на границе слоя Рассмотрим плоскодеформированное состояние полосы в декартовой системе координат, верхняя граница которой свободна от напряжений. (63) а к нижней границе приложена некоторая система нагрузок (64) (65) Здесь - заданные функции. В основном изложение следует работам [19-22]. Найдем связь между действующей на нижней границе полосы системой нагрузок и производными от компонент смещений . Полная система уравнений теории упругости включает: 1) уравнения равновесия в отсутствие объемных сил (66) (67) где - компоненты тензора напряжений; 2) уравнения закона Гука, выраженные через компоненты смещений (68) (69) где - модуль Юнга и коэффициент Пуассона полосы; 3) уравнение совместности, выраженное через компоненты тензора напряжений (70) Для решения поставленной задачи подвергнем уравнение (69) следующим преобразованиям. Подставим значение из уравнения (68) в продифференцированное по уравнение (69). (71) Подстановка сюда уравнений равновесия (67) дает (72) Для решения системы уравнений (66), (67), (68), (72), (70) используем двустороннее преобразование Лапласа, определяющее образ функции (73) Обратное преобразование при этом определяется как (74) где направление обхода по контуру , соответствующему мнимой оси, осуществляется сверху вниз. Трансформанты от уравнений (66), (67), (68), (72), (70) и граничных условий (63), (64), (65) примут вид (75) (76) (77) (78) (79) (80) (81) (82) Здесь (83) а (84) есть трансформанты от напряжений на нижней границе. Решение системы обыкновенных дифференциальных уравнений (75)-(79) с граничными условиями (80)-(82) дает выражения для трансформант производных смещения через действующие на нижнем крае полосы напряжения. Для поставленной задачи интерес представляют значения трансформант производных смещения на этой же линии (85) (86) Здесь (87) Коэффициенты в полученных выражениях аналитичны на всей мнимой оси за исключением нуля, где они имеют полюса до третьего порядка включительно. Полюс третьего порядка в нуле трансформанты соответствует квадратичному росту оригинала на бесконечности. Следовательно, из полученных выражений, в частности, следует, что при отсутствии особенностей для напряжений производные от смещения могут расти как квадрат расстояния, а сами смещения как куб расстояния. Для полосы такое поведение является общеизвестным фактом, в частности, в рамках элементарной балочной теории, являющейся при определенных условиях строгой асимптотикой точной теории, смещения вдали от приложенной системы сил с ненулевым главным вектором растут как третья степень расстояния. Приложение 2. Разложение факторизующих функций вблизи нуля Пусть существует факторизация функции (88) (89) Пусть функция четная . Тогда функция тоже четная. Согласно формулам Сохоцкого-Племеля предельные значения выражаются следующим образом: (90) (91) Из последнего равенства с учетом четности для следует . Подставляя полученное равенство в , получаем (92) Для получения следующих членов разложения вычислим k-ю логарифмическую производную от обеих частей равенства (88): (93) На основании формул Сохоцкого-Племеля предельные значения функций в нуле следующие: (94) (95) причем если функция четная, то для всех нечетных в правой части равенства (94) будут обнуляться первые члены, а для всех четных - вторые. Выпишем несколько первых членов получаемого таким образом разложения. Для первой логарифмической производной (96) (97) Для второй логарифмической производной (98) С учетом (96) данное выражение преобразуется к виду (99) Для третьей производной получаем (100) (101) С учетом (96) и (99) выражение (100) преобразуется к виду (102) В результате для разложения факторизующих функций вблизи нуля может использоваться следующее представление: (103) Данный вывод не представляет сложности, однако в литературе (например, [40, 41]) обычно ограничиваются лишь двумя членами разложения.

About the authors

K B Ustinov

A.Yu. Ishlinsky Institute for Problems in Mechanics RAS

References

Statistics

Views

Abstract - 178

PDF (Russian) - 52