Лекции по механике разрушения
Подождите немного. Документ загружается.
l
a
lK
I
arcsin
2
1
.
Рис. 50 Рис. 51 Рис. 52
Для примера на рис. 50 на правом и левом конце соответственно
lKlK
II
4
1
,
4
3
.
Для щели, располагающейся вдоль длины окружности (рис. 51), КИН в
наиболее опасной точке О имеют вид
2
3
2sin
2
sinsin
2
sin
2
cos
2
3
2cos
2
5
2cos
2
sin
2
cos
2
sin12
sin
2422
2
R
K
I
2
3
2cos
2
sinsin
2
sin
2
cos
2
3
2sin
2
5
2sin
2
sin
2
sin
2
sin12
sin
222
2
R
K
II
.
В случае перешейка между двумя полубесконечными щелями имеют место
следующие формулы
- у левого конца перешейка
;
2
,
2
2
2222
a
V
aca
acQ
K
aa
M
a
P
aca
acP
K
III
- у правого конца перешейка
.
2
,
2
2
2222
a
V
aca
acQ
K
aa
M
a
P
aca
acP
K
III
Разрез в полуплоскости. Задача о растяжении упругой полуплоскости с краевой
щелью (рис. 53) рассматривалась многими авторами. Ее решение получено многими
способами – как точными аналитическими методами, так и приближенными:
61
lK
I
12.1
.
Коэффициент 1.12 вычислен с погрешностью 1%. Как видно, влияние свободной
границы тела приводит к увеличению КИН на 12%.
Для задачи на рис. 54
lK
I
68.0
.
Коэффициент 0.68 вычислен с ошибкой 3%.
В случае полубесконечного разреза, приближающегося к свободному краю
полуплоскости (рис. 55) КИН равен
cc
M
c
P
K
I
8
84
8
124
22
.
Рис. 53 Рис. 54 Рис. 55
Разрезы, исходящие из круглого отверстия. Решена задача о всестороннем и
одностороннем растяжении плоскости с одной и двумя щелями, исходящими из
кругового отверстия (рис. 56).
Рис. 56 Рис. 57 Рис. 58
Имеем
rllFK
I
.
62
Решена задача о всестороннем растяжении плоскости со звездообразной щелью
(рис. 57). Окончательный результат имеет вид
lnK
I
,
где
(т) – коэффициент, зависящий от числа разрезов n.
Для круглого диска с внутренней щелью, раздавливаемого двумя
сосредоточенными силами (рис. 58),
862
64
3
4
3
2
3
1
R
l
R
l
R
l
R
l
R
P
K
I
.
Рис. 59.
Тела сложной формы. Изучены многие случаи различных криволинейных
отверстий со щелями. В случае всестороннего растяжения упругой плоскости с
гипоциклоидальным отверстием, контур которого описывается уравнениями
20,0
,sin
1
sin
1
,cos
1
cos
1
a
n
nn
na
y
n
nn
na
x
,
где n – целое положительное число, КИН равен
1
2
n
an
K
I
.
Гипоциклоида имеет n+1 точку возврата, каждая из которых с точки зрения
концентрации напряжений эквивалентна концу трещины (на рис. 59 изображена
астроида с n=3).
Пространственные задачи. В силу сложности решения пространственных
задач готовых решений накоплено значительно меньше, чем для плоского случая.
Приведем окончательные результаты вычисления КИН в пространственных задачах.
Дискообразная щель. В общем осесимметричном случае дискообразного разреза
вдоль z=0, x
2
+y
2
<a
2
в безграничном пространстве (рис. 60) имеем
63
Рис. 60.
a
z
III
a
z
r
II
a
z
I
ra
drrr
a
K
ra
drrr
a
K
ra
drrr
a
K
0
22
0
22
0
22
.
2
,
2
,
2
Здесь предполагается, что нагрузки симметрично приложены к верхнему и нижнему
берегам щели, а на бесконечности напряжения исчезают.
При чистом изгибе стержня с дискообразной щелью (рис. 61) КИН будет
следующим
cos3/21/2
1
baJabMK
yyI
.
Здесь J
y
– соответствующий момент инерции поперечного сечения стержня, b –
расстояние центра щели от нейтральной линии.
Рис. 61 Рис. 62
При скручивании круглого цилиндрического стержня с дискообразной щелью
(рис. 62) имеем (здесь предполагается, что центр щели лежит на оси стержня, а
плоскость щели перпендикулярна этой оси)
4
2
Ra
aM
K
z
III
.
64
Итак, имеется большое количество решенных аналитическим способом задач по
определению КИН для тел с трещинами различной конфигурации, но все же число
этих задач ограничено.
14.2. Метод конечных элементов
В настоящее время известны различные методы расчета КИН, использованные
на использовании МКЭ. Они имеют свои преимущества и недостатки. Рассмотрим
некоторые из них.
14.2.1.Прямой метод
В этом случае для определения КИН в асимптотические формулы для
компонент напряжений и перемещений в вершине трещины непосредственно
подставляют возникающие в окрестности вершины трещины напряжения или
перемещения, полученные с помощью МКЭ. При использовании напряжений этот
способ называется прямым методом напряжений, а при использовании
перемещений – прямым методом перемещений. Значения напряжений или
перемещений, которые необходимо подставить в формулы, следует выбирать
такими, чтобы они были определяющими для рассматриваемого типа
деформирования. Например, в случае трещины типа I необходимо подставить
компоненту напряжения
y
, которая действует на оси х в направлении оси у. По
сравнению с методом напряжений метод перемещений дает более надежные
результаты с точки зрения точности.
Расчет КИН прямым методом требует, чтобы решение для напряжений и
перемещений имело достаточно высокую точность в окрестности вершины
трещины. При использовании обычных элементов нельзя отразить особенность при
приближении непосредственно к вершине трещины. Поэтому нельзя ожидать, что
точность расчета напряжений и деформаций в окрестности вершины трещины будет
высокой. Следовательно, установленные по этим параметрам КИН также будут
иметь не слишком хорошую точность. Поэтому при определении КИН с помощью
обычных элементов необходимо предусмотреть меры, позволяющие улучшить
точность решения. К таким мерам можно отнести следующие.
(1) У вершины трещины, насколько это возможно, желательно использовать
разбиение на малые элементы (мелкую дискретизацию).
(2) Часто могут возникать ситуации, например при решении трехмерных задач,
когда разбиение на мелкие элементы оказывается нерациональным с точки зрения
вычислительных ресурсов. В таких случаях сначала находят решения для грубого
разбиения, а затем выделяют окрестность вершины трещины, выполняют мелкое
разбиение и решают задачу. При этом в качестве граничных условий используют
перемещения узлов и узловые силы, полученные в предыдущем решении. Такой
способ решения называют поэтапным.
(3) Используя асимптотические формулы для компонент напряжений и
перемещений в вершине трещины, в точках, расположенных на различных
расстояниях r от вершины трещины, расчетным путем определяют значения К и
65
строят график зависимости К от r. Искомое значение К, соответствующее
рассматриваемой трещине, принимают равным тому значению, которое получается
в результате экстраполяции при r 0 (рис. 63). Следует иметь в виду, что значения
К, соответствующие точкам, наиболее близко расположенным к трещине, обладают
не очень хорошей точностью. Поэтому при проведении экстраполяции эти значения
исключают.
Рис. 63. Оценка КИН путем экстраполяции
(4) Используют всевозможные разбиения от грубых до мелких и определяют
значения К, проводя так же, как в п. 3, экстраполяцию при r 0.
В пп. (1) и (3) дискретизация выполняется один раз, а в пп. (2) и (4) – несколько
раз, что затрудняет решение задачи. Пункты (1) и (3), а также их сочетание можно
отнести к таким методам оценки К, когда необходимо, используя обычные элементы
и не выделяя особенностей, довольно быстро определить КИН. В других случаях
необходимо использовать специальную дискретизацию.
14.2.2. Энергетический метод
Известно несколько энергетических методов, наиболее распространенным из
которых является метод полной энергии. Интенсивность G освобождения энергии,
которое имеет место при распространении трещины в упругом теле, связана с КИН
зависимостью
222
2
1
2
1
2
1
IIIIII
KKKG
.
При распространении трещины на величину А изменение энергии деформации
составляет U. Интенсивность освобождения энергии можно также представить как
G=U/А. Таким образом, если в результате использования МКЭ установить
интенсивность освобождения энергии G, можно определить КИН. Однако видно,
что найденный таким образом КИН должен относиться к трещине одного из типов,
т.е. для этих типов трещин расчет К возможен только порознь.
При определении интенсивности освобождения энергии можно воспользоваться
изложенными ниже двумя способами.
А. Методом конечных элементов проводят численное решение для трещины
площадью А. При этом определяют энергию упругой деформации U(A). Затем
рассматривают рост трещины на один элемент. При тех же граничных условиях
вновь МКЭ решают задачу и определяют U(A+А). Тогда G=( U(A+А) - U(A))/ А.
66
Энергию упругой деформации можно установить путем интегрирования энергий
деформации в отдельных точках по всей рассматриваемой области. Полученное
значение должно быть равно работе внешних сил. Однако расчет работы внешних
сил довольно прост и позволяет исключить сложности, связанные с
интегрированием.
Б. Можно не рассматривать рост трещины на один элемент, а за счет
соответствующего смещения координаты вершины трещины задавать ее
распространение (рис. 64). В результате такой операции будет происходить
изменение жесткости. При этом можно считать, что жесткость меняется лишь у тех
элементов, которые окружают вершину трещины. Если принять это во внимание и
путем исключения решать систему линейных уравнений, можно за одно решение
определить U= U(A+А) - U(A). При этом необходимо в процессе исключения
вначале соответствующим образом вводить указанную операцию распространения
трещины.
Рис. 64. Виртуальный рост трещины
В случае использования метода А необходимо по отдельности дважды
проводить решение. В рассматриваемом же методе в этом нет необходимости. Его
можно эффективно использовать при большом числе уравнений, в частности при
решении трехмерных задач. Метод Б носит название метода виртуального роста
трещины, а метод А – метода податливости.
14.2.3. Использование специальных элементов
Методы, использующие аналитические решения. На рис. 65 изображен
треугольный элемент с трещиной, предложенный Бисковым. В качестве решения
вводится решение, полученное из функции напряжений Мусхелишвили. При этом
аналитическое решение имеет особенность напряжений, характерную для вершины
трещины. В случае такого элемента восьми степеням свободы узлов соответствует
восемь независимых переменных: три переменные, характеризующие перемещения
твердого тела, три средних напряжения и КИН К
I
и К
II
.
Рис. 65. Специальные элементы Бискова и Уилсона
67
Руководствуясь аналогичными соображениями, Уилсон предложил круглый
элемент с трещиной (рис. 65) и нашел решение, используя до высоких порядков
разложение функции перемещения в окрестностях трещины.
Необходимо отметить, что при использовании этих элементов в реальных
условиях обычно их окружение составлено из треугольных элементов. При этом
обеспечиваются такие условия, при которых перемещения смежных элементов в
узлах совпадают, однако в общем случае не выполняется условие неразрывности
перемещений на границах элементов, за исключением узлов. Поэтому не
существует теоретического обоснования сходимости. Тем не менее, численные
эксперименты показали, что при этом не возникает существенных проблем.
Рис. 66. Специальный элемент Уилсона
Методы, использующие простые функции для описания особенностей. На рис.
66 показан предложенный Уилсоном элемент, представляющий собой
равнобедренный треугольник, вершина которого находится в вершине трещины.
Для этого элемента перемещение можно представить таким образом
e
j
ij
i
i
ij
j
r
r
uuuu
0
.
Аналогично можно определить перемещение v. Такой элемент для напряжений и
деформаций имеет особенность в вершине трещины.
Деформированные изопараметрические элементы. У обычного восьмиузлового
изопараметрического элемента узлы находятся на серединах сторон. Если сместить
узлы в точки, соответствующие ¼ длины стороны, то можно получить особенность
в вершине трещины (рис. 67). Однако в этом случае будет отсутствовать сходимость
решения и энергия деформации будет иметь особенность в вершине трещины. В
связи с этим предложено считать длину одной из сторон четырехугольника равной
нулю. В результате объединения смежных узлов удалось получить элемент, который
одновременно обладал всеми необходимыми свойствами (рис. 67).
Рис. 67. Деформированные элементы
68
Такой элемент может быть использован как для упругого тела, так и для идеально
пластического без деформационного упрочнения.
15. Управление поведением трещин на основе линейной механики
разрушения
Одно из наиболее важных для инженерной практики применение методов
механики разрушения – управление поведением трещин. Здесь речь идет о том, что
при обнаружении трещины есть возможность тормозить и (или) останавливать ее
развитие, а также принимать меры, чтобы в данном месте конструкции трещина
больше не возникала. Достижение этих целей требует, как правило, некоторых
конструктивных изменений в поврежденный узел или элемент конструкции. И тогда
необходимо вести речь о проектировании модернизации. Другими словами, имея
дело уже с возникшей при эксплуатации трещиной, или предполагая, исходя из
анализа напряженно-деформированного состояния, возможность ее возникновения,
необходимо принятие соответствующих проектно-конструкторских решений. В
большинстве случаев они сводятся к конструкционному торможению трещин с
помощью ребер жесткости, ремонтных заплат, разгружающих отверстий и
стопперов.
Конструкционное торможение трещины ребрами жесткости. Одной из
важнейших характеристик сопротивления материала трещинообразованию является
величина предельной нагрузки, связанная с началом развития трещины. Однако это
справедливо только в случае лавинообразного неустойчивого распространения. Во
многих случаях взаимодействия трещин с препятствиями и границами, как
показывают эксперименты и расчеты, на значительном участке изменения нагрузки
развитие трещины протекает устойчиво. Очевидно, что наличие устойчивых трещин
в конструкциях, работающих зачастую в определенных режимах изменения
внешних нагрузок, гораздо менее опасно, а искусственное усиление таких
конструкций (за счет постановки заклепок, пластин и стрингеров, высверливания
отверстий на пути распространения трещин и т.д.) может значительно продлить их
жизнь.
Для предотвращениия катастрофического развития трещины и разрушения
конструкции трещины подкрепляют ребрами жесткости, препятствующими их
распространению (рис. 68). Возможна следующая простейшая схематизация этой
задачи. Бесконечная пластина единичной толщины с трещиной растягивается на
бесконечности однородным напряжением (рис. 69); действие подкрепляющих ребер
заменяется четырьмя симметрично расположенными сосредоточенными силами,
приложенными в местах расположения ближайших к трещине заклепок; величина
этих сил считается заданной. Как показывают оценки, действием более удаленных
заклепок можно пренебречь.
69
Рис. 68. Панель с приклепанными ребрами жесткости
Рис. 69. Схематизация задачи о трещине в листе с ребрами жесткости
На рис. 70 представлены зависимости критического напряжения от длины
трещины и местоположения заклепок.
Рис. 70. Зависимость критического напряжения от длины трещины для различных
положений точек приложения сосредоточенных сил, имитирующих действие
заклепок. Линия 1 соответствует отсутствию заклепок. Величина равна: для линий
2 – 0.15; 3 – 0.25; 4 – 0.4; 5 – 0.5; 6 – 0.75
Здесь же построена кривая 1, соответствующая отсутствию подкрепляющих ребер.
Очевидно, что трещина устойчива, если напряжение, необходимое для ее
поддержания в критическом состоянии, возрастает с увеличением длины трещины.
Как видно, существует критическое значение
0
безразмерного параметра = y
0
/L
такое, что при >
0
кривая не имеет участков возрастания, так что трещина всегда
неустойчива, а при <
0
имеется участок возрастания, на котором трещина
устойчива.
В этом решении сила Р, характеризующая действие заклепок, считалась
постоянной. На самом деле величина этих сил вполне определяется упругими и
геометрическими характеристиками рассматриваемой задачи, а также величиной
70