CC BY
24
УДК 539.375
Равновесная плоская трещина с угловыми точками контура в упругом слое
Б. В. Соболь, Е. В. Рашидова, Е. В. Борисова, С. Б. Петренкова
(Донской государственный технический университет)
Рассмотрена трёхмерная статическая задача теории упругости о равновесии упругого слоя, ослабленного плоской прямоугольной трещиной. Трещина расположена в срединной плоскости слоя, поддерживается в раскрытом состоянии под действием нормальной нагрузки, приложенной к её берегам. Гоани слоя находятся в условиях гладкого контакта с двумя жёсткими основаниями. Применением двумерного интегрального преобразования Фурье к уравнениям равновесия задача сведена к решению известного сингулярного интег-ро-дифференциального уравнения относительно функции раскрытия трещины. Решение уравнения строится прямым вариационным методом. В окрестности угловых точек контура решение достраивается численно, с учётом заранее выделенной особенности. Получены значения коэффициента интенсивности нормальных напряжений в окрестности контура трещины. Установлены особенности поведения решения в окрестности прямолинейных участков и угловых точек контура.
Ключевые слова: трещина, контур, коэффициент интенсивности напряжений, слой, сингулярное интегральное уравнение, вариационный метод, силовой критерий разрушения, угловая точка, окрестность, относительная толщина слоя
Введение. Как было установлено в многочисленных работах отечественных и зарубежных исследователей (например, [1]), в точках гладкости контура трещины функция раскрытия трещины должна иметь асимптотику
где nix, у) — расстояние по нормали от внутренней точки области до контура трещины, N(x, у) — некоторая регулярная функция.
В случаях, когда контур трещины имеет особые точки: во-первых, нарушается состояние плоской деформации в окрестности этих точек, во-вторых, меняется показатель особенности поля напряжений и перемещений в окрестности контура вблизи таких точек. Определению показателей при особенности в зависимости от угла излома контура посвящён ряд работ отечественных и зарубежных исследователей [2].
Если в плоскости расположения трещины ввести полярную систему координат с полюсом в особой точке контура, то функция раскрытия трещины в окрестности этой точки асимптотически может быть, описана соотношением
где р, а — полярные координаты точки, ¿(а) — некоторая регулярная функция.
Рассмотрим наиболее распространённый частный случай, когда контур трещины в точке излома образует прямой угол. В этом случае установлен следующий показатель особенности в соотношении (2): е = 0,816.
Постановка задачи. Продемонстрируем методику решения задач теории упругости для трещин с особыми точками контура на примере задачи о прямоугольной трещине в слое.
Трещина расположена в срединной плоскости слоя, поддерживается в раскрытом состоянии под действием нормальной нагрузки, приложенной к её берегам. Грани слоя находятся в условиях гладкого контакта с двумя жёсткими основаниями. Применением двумерного интегрального преобразования Фурье к уравнениям равновесия в перемещениях задача о продольной равновесной
Y(x,y) = N(x,y)nm (х,у),
(1)
(2)
трещине в слое сводится к решению двумерного сингулярного интегро-дифференциального уравнения [3]:
А ifY (5/'л) ^+JJ Y (5, П) S (х, у, п) c£d] = ~ р(х, у),
2п
R
0
(3)
здесь (л> у) е О., О. — область расположения трещины в плане, (1 = х - ^)2 + (у - п)2/
з2 з2
0 = Е / [2(1 - у)], Л = —- +—-, р{х, у) — нагрузка, приложенная к берегам трещины, Е— модуль
у ,л дх ду
Юнга, V — коэффициент Пуассона, 5(х,у,^,г\) — регулярная часть ядра уравнения, характеризующая влияние геометрических параметров задачи: размеров трещин в плане, относительных расстояний между ними и (или) до границ рассматриваемых тел. Для задачи о трещине, лежащей в срединной плоскости слоя, имеем [3]:
1 °°
5(х,у,?„Ц) = -^Е№ /11), ^ (£) = \и[Н (и)- £/]_70 ,
" О
Н {и) = и (зЫи +2и)(сЬи -1) \ J0 {х) ~ функция Бесселя 1-го рода.
Решение уравнения (3) будем строить прямым вариационным методом, описанным в [4]. В рассматриваемом случае область интегрирования описывается неравенствами: 0:|>г|<а, \у\^Ь (рис. 1).
Метод решения. Выберем, систему координатных функций таким образом, чтобы установленная выше асимптотика (2) с соответствующим показателем явно учитывалась. Ограничимся рассмотрением, как наиболее характерного частного случая, равномерно распределённой нагрузки р (х,у) на берегах трещины, где р(х,у) = р = const. Решение задачи при этом строится в виде:
Ф (х,у) =
(
Y(*,У) = ^яФ1/2 (х,у) X X A™ cosT^L.cos
У /77=0 /7=0 & &
(4)
Коэффициенты Атп в разложении (4) определяются из условия минимума функционала:
'"У 2
J (у) = fjY {х,у) ||Y (5,п) К2 (5 - Х,п - y)dE,di]dxdy - JL JJ Y (х, у) pdxdy, (5)
Q Q Q
со со ______ _______
где К2 (а, 3) = [ j л/у2 +v2L (Ьл1и2 +v2 J cos а у cos 3 vdudv .
О О
Следует отметить, что в аналогичном виде вариационный метод был использован при исследовании задачи о равновесии упругого пространства ослабленного трещиной [5]. В качестве
координатных функций в этой работе были использованы степенные функции с множителем, явно содержащим требуемую асимптотику в окрестности контура трещины.
Минимизация функционала (5) в рассматриваемом случае приводит к следующей системе линейных алгебраических уравнений для определения коэффициентов Атп:
М N
=2гГ (п/л,гш), (6)
/=0 ;=0
где т = 0,1,2,п = 0,1,2,...,Л/.
В данной задаче /"(0,3) = ^—При вычислении коэффициентов системы (6) ис-
о 3
пользовано значение интеграла
С| ф1/2 (х, у)соз^-со5^-с1хс1у = ^-abf (а, 3),
полученного с помощью (3.752) из [6].
Вычисление элементов четырёхмерной матрицы <?..тп производится с использованием
представлений (5.4) из [4].
Сравнительный анализ и сопоставление численных результатов. Конечной целью исследования задачи теории упругости для тела с трещиной, как известно, является определение критических нагрузок, при которых трещина становится неустойчивой и начинает расти. Поскольку в качестве основного критерия разрушения в данной работе принят силовой критерий [7], перейдём к определению коэффициента интенсивности нормальных напряжений в окрестности контура трещины в каждом из рассмотренных случаев. Особое внимание в дальнейшем также будет уделяться сравнительному анализу применяемых здесь методов и сопоставлению численных результатов.
Контроль точности решения задачи вариационным методом приводился сопоставлением результатов при различном количестве используемых координатных функций (Ми IV). В частности, вычислялась амплитуда раскрытия берегов трещины в центре при различных соотношениях геометрических параметров задачи. Во всех рассмотренных ниже случаях расхождение результатов счёта при Л/ = /V = 4 и А/ = /V = 6 не превысило 3,5 %.
Кроме того, о достоверности получаемых результатов могут свидетельствовать и следующие качественные рассуждения. При вытягивании трещины вдоль одной из осей симметрии, амплитуда её раскрытия в центре, а также численные значения коэффициента интенсивности нормальных напряжений в средней части большей из сторон прямоугольного контура должны быть близки к известным соответствующим решениям задач о полосовой трещине и эллиптической трещине с такими же полуосями. Причём, поскольку упомянутая полосовая трещина ширины 2Ь содержит внутри себя прямоугольную область, а эллиптическая, в свою очередь, содержится в последней, решение задачи о прямоугольной трещине, по известному принципу сравнения [8], должно быть заключено между соответствующими решениями этих двух задач. В таблице 1 представлены приведённые значения амплитуды вертикальных перемещений у (0,0)0 /(ра) в центре
трещины для полосовой [3], эллиптической [9] и прямоугольной трещины. Вычисления проведены для предельного случая Л ^ оо. Для сопоставления, здесь же приведён соответствующий результат работы [5].
В случаях эллиптической и прямоугольной трещин при вычислениях принято следующее соотношение размеров Ь/а = 0,5. Вычисления показывают, что и при конечной толщине слоя установленные закономерности не нарушаются.
Таблица 1
Приведённые значения амплитуды вертикальных перемещений в центре полосовой,
прямоугольной и эллиптической трещины
(4) [5] [9] [3]
0,882 0,872 0,826 1,000
Перейдём к анализу численных значений коэффициента интенсивности нормальных напряжений К, при обходе по контуру прямоугольной трещины. В этой задаче более удобно рассматривать не относительные значения параметра К, (как в задачах об эллиптических трещинах, где точное аналитическое решение для предельного случая Л ^ оо существует), а абсолютные. Вдоль каждой из сторон прямоугольной области значения коэффициента интенсивности нормальных напряжений будем вычислять по формулам:
[у М N /Т7П V
Kl3(x) = paJ-/(x,0)££A„,„cos—-, (7)
V D /?7=0 /7=0 О
К,ь {У) = раЛ/(0,у)X X Amn cos^. (8)
1 а /77 = 0/7 = 0 о
Графики изменения коэффициента интенсивности нормальных напряжений вдоль большей и меньшей сторон прямоугольника представлены на рис. 2 (а, б). Вычисления проведены для ряда
значений параметра К = h / а при Ь / а = 0,5. Отметим, что максимальные амплитудные значения
К, достигает на большей стороне. Опираясь на силовой критерий разрушения, можно сделать вывод, что при достижении приложенной нагрузки своего критического значения, рост прямоугольной трещины начнётся с середин больших сторон. На рис. 3 представлены соответствующие графики, иллюстрирующие решение задачи для случая трещины квадратной формы (Ь = а).
С приближением к угловым точкам контура трещины коэффициент интенсивности нормальных напряжений стремится к нулю. Это обусловлено особым видом асимптотики (2) функции раскрытия трещины в окрестности этих точек, с соответствующим показателем. Перейдём к построению решения задачи у(х, у) в окрестности угловых точек.
Вдоль лучей, выходящих из угловой точки контура под различными углами а, вычисляются значения функции: ¿(а) = у(х,у) / (р /а).
Фрагменты результатов таких вычислений представлены в таблице 2 (при b/а = 0,5) и таблице 3 (при У а = 1).
В каждом из случаев, представленных в таблицах, нетрудно обнаружить интервал изменения отношения р/a e\_rur2~\, в котором функция с достаточно высокой точностью сохраняет
стабильность вдоль каждого луча и изменяется лишь при переходе с одного луча на другой. Упомянутая закономерность ранее была установлена для предельного случая рассматриваемой здесь задачи, соответствующего Л^оо [5], [10].
Расположение и размеры этой зоны зависят от а, отношения сторон области трещины b/а и относительной толщины слоя Л. Таким образом, определив для каждого значения величину ¿(а) на указанном отрезке [/i, г{[, представляется возможным достроить искомую функцию у(х, у) по асимптотической формуле (2) непосредственно до вершины угла трещины.
Графики функции ¿(о) при различных относительных толщинах слоя представлены на рисунке 4 (слева — Ь/а = 0,5, справа — Ь/а = 1).
Следует отметить симметричность функции у(х У) в окрестности угловой точки контура относительно биссектрисы. Между тем, сама прямоугольная область расположения трещины такой симметрией не обладает.
О 0,25 0,5 0,75 х/а
6) Щ
рл'я 2,5
2
1,5
1 0 0,33 0,66 Х/а
Рис. 2. Изменение коэффициента интенсивности нормальных напряжений вдоль большей (а) и меньшей (6) сторон прямоугольной трещины
Рис. 3. Изменение коэффициента интенсивности нормальных напряжений для случая трещины квадратной формы (Ь = а)
Таблица 2
Приведённые значения функции раскрытия вдоль лучей, выходящих из угловой точки
прямоугольной трещины (Ь/а = 0,5)
1 Л 2
а = п/4 а = п/8 а = Зп/8 р/а а = п/4 а = п/8 а = Зп/8
2,38 2,04 2,01 1/16 2,60 2,24 2,21
2,59 2, 21 2,18 1/8 2,84 2,48 2,39
2,67 2,29 2,22 3/16 2,92 2,52 2,43
2,66 2,30 2,20 1/4 2,92 2,53 2,41
2,60 2,29 2,14 5/16 2,36 2,53 2,36
2,55 2,24 2,07 3/8 2,80 2,48 2,28
2,42 2,18 1,99 7/16 2,67 2,41 2,19
2,31 2,07 1,90 1/2 2,57 2,33 2,10
Таблица 3
Приведённые значения функции раскрытия вдоль лучей, выходящих из угловой точки
квадратной трещины [Ь/а = 1)
1 Л 2
а = п/4 а = п/8 р/а а = п/4 а = п/8
1,91 2,27 1/16 2,28 2,70
2,00 2,51 1/8 2,51 3,00
2,20 2,61 3/16 2,63 3,13
2,24 2,64 1/4 2,69 3,18
2,25 2,64 5/16 2,71 3,19
2,23 2,63 3/8 2,73 3,19
2,19 2,58 7/16 2,67 3,15
2,14 2,51 1/2 2,63 3,09
b/a=0,5
_4Є вд ра
b/a=l
49 ,
Рис. 4. Изменение множителя —Lía) при особенности функции раскрытия трещины в окрестности
ра
угловой точки для прямоугольной (а) и квадратной (б) формы трещины
Выводы
1. В результате проведения численного эксперимента установлено, что решение задачи о прямоугольной трещине в упругом слое удовлетворяет принципу сравнения.
2. Коэффициент интенсивности напряжений достигает своего максимального значения при обходе по контуру в окрестности центра длинной стороны прямоугольной области.
3. Установлена симметричность функции раскрытия трещины в окрестности угловой точки контура, независимо от соотношения размеров прямоугольной области.
Библиографический список
1. Sneddon, I. N. The stress intensity factor for a flat elliptical crack in an elastic solid under uniform tension / I. N. Sneddon // Int. J. Eng. Sei. — 1979. —V. 17. — № 2. — p. 92—103.
2. Bazant, Z. Tree-dimentional harmonic functions near termination or interaction of gradient singularity lines: A general numerical method / Z. Bazant // Int. J. Eng. Sei. — 1974. — № 12. — p. 221—243.
3. Александров, В. М. Тонкие концентраторы напряжений в упругих телах / В. М. Александров, Б. И. Сметанин, Б. В. Соболь. — Москва: «Физматлит», 1993. — 224 с.
4. Сметанин, Б. И. Равновесие упругого слоя, ослабленного системой плоских трещин / Б. И. Сметанин, Б. В. Соболь // ПММ. — 1984. — Т. 48. — №. 6. — С. 1030-1038.
5. Гольдштейн, Р. В. Качественные методы в механике сплошных сред / Р. В. Гольдштейн, В. М. Енотов. — Москва: Наука, 1989. — С. 110—115.
6. Градштейн, И. С. Таблицы интегралов, сумм, рядов и произведений / И. С. Градштейн, И. М. Рыжик. — Москва: Наука, 1971. — 1108 с.
7. Irvin, G. R. Analysis of stress and strain near the end of a crack, traversing a plate/G. R. Irvin // J. Appl. Mech. — 1957. — № 3. — P. 361—364.
8. Гольдштейн, P. В. Вариационные оценки для коэффициента интенсивности напряжений на контуре плоской трещины нормального разрыва / Р. В. Гольдштейн, В. М. Енотов // Изв. АН СССР, МТТ. — № 3. — С. 59—64.
9. Пэрис, П. Анализ напряжённого состояния около трещин / П. Пэрис, Дж. Си // Прикладные вопросы вязкости разрушения. — Москва: Мир. — 1968. — С. 64—142.
10. Мураками, Ю. Справочник по коэффициентам интенсивности напряжений, (в 2-х томах) / Ю. Мураками. — Москва: Мир. — 1990. — Т. 1 — 448 с. — Т. 2 — 1014 с.
Материал поступил в редакцию 08.09.2012.
References
1. Sneddon, I. N. The stress intensity factor for a flat elliptical crack in an elastic solid under uniform tension / I. N. Sneddon // Int. J. Eng. Sci. — 1979. —V. 17. — № 2. — p. 92—103.
2. Bazant, Z. Tree-dimentional harmonic functions near termination or interaction of gradient singularity lines: A general numerical method / Z. Bazant // Int. J. Eng. Sci. — 1974. — № 12. — p. 221—243.
3. Aleksandrov, V. M. Tonkie koncentratory' napryazhenij v uprugix telax/ V. M. Aleksan-idrov, В. I. Smetanin, В. V. Sobol'. — Moskva: «Fizmatlit», 1993. — 224 s. — In Russian.
4. Smetanin, В. I. Ravnovesie uprugogo sloya, oslablennogo sistemoj ploskix treshhin / В. I. Smetanin, В. V. Sobol' // PMM. — 1984. — T. 48. — №. 6. — S. 1030—1038. — In Russian.
5. Gol'dshtejn, R. V. Kachestvenny'e metody' v mexanike sploshny'x sred / R. V. Gol'dshtejn, V. M. Enotov. — Moskva: Nauka, 1989. — S. 110—115. — In Russian.
6. Gradshtejn, I. S. Tablicy' integralov, summ, ryadov i proizvedenij / I. S. Gradshtejn, I. M. Ry'zhik. — Moskva: Nauka, 1971. — 1108 s. — In Russian.
7. Irvin, G. R. Analysis of stress and strain near the end of a crack, traversing a plate / G. R. Irvin // J. Appl. Mech. — 1957. — № 3. — P. 361—364.
8. Gol'dshtejn, R. V. Variacionny'eocenki dlya koe'fficienta intensivnosti napryazhenij na konture ploskoj treshhiny' normal'nogo razry'va / R. V. Gol'dshtejn, V. M. Enotov // Izv. AN SSSR, MTT. — № 3. — S. 59—64. — In Russian.
9. Pe'ris, P. Analiz napryazhyonnogo sostoyaniya okolo treshhin / P. Pe'ris, Dzh. Si // Prikladny'e voprosy' vyazkosti razrusheniya. — Moskva: Mir. — 1968. — S. 64—142. — In Russian.
10. Murakami, Yu. Spravochnik po koe'fficientam intensivnosti napryazhenij. (v 2-x tomax) / Yu. Murakami. — Moskva: Mir. — 1990. — Т. 1 — 448 s. — T. 2 — 1014 s. — In Russian.
EQUILIBRIUM PLANAR CRACK WITH CONTOUR ANGULAR POINTS IN ELASTIC LAYER
В. V. Sobol, E. V. Rashidova, E. V. Borisova, S. B. Petrenkova
(Don State Technical University)
The three-dimensional static problem of the elasticity theory on the elastic layer balance weakened by a flat rectangular crack is considered. The crack is located in the median layer plane; it is supported open under the normal loading applied to its edges. The layer edges are under the smooth contact with two rigid bases. The problem is converted to the solution of the known singular integro-differential equation concerning the crack opening function through the application of two-dimensional integral Fourier transformation to the equilibrium equation. The equation solution is constructed by the direct variational method. In neighborhood of the contour angular points, the solution completes the construction numerically, with regard for the early distinguished singularity. Values of the direct stress intensity factor in neighborhood of the crack periphery are obtained. The solution behavioral features in neighborhood of the contour rectilinear sites and angular points are established.
Keywords: crack, contour, stress intensity factor, layer, singular integral equation, variational method, force fracture criterion, salient point, neighborhood, layer relative thickness.
