CC BY
41
УДК 539.3
Диаграммы квазихрупкого разрушения при усталости (двухчастотное нагружение)
В.М. Корнев
Институт гидродинамики им. М.А. Лаврентьева СО РАН, Новосибирск, 630090, Россия
Рассматривается распространение трещины скачками в казихрупких материалах при цикловом нагружении. Изучается стационарное и нестационарное нагружение при пульсирующем приложении нагрузки, для последнего типа нагрузки особое внимание уделено двухчастотному нагружению.
Предлагается использовать для анализа указанного процесса диаграммы квазихрупкого разрушения тел при циклическом нагружении. Одна из кривых предлагаемой диаграммы напоминает диаграмму Китагава-Такахаши.
В явном виде получены оценки безразмерной средней скорости продвижения вершины трещины за один цикл нагружения при скачкообразном подрастании трещины. Полученные соотношения для средней скорости можно рассматривать как структурные формулы для построения кривых Пэриса. Для кривых Пэриса проведен подробный анализ процесса продвижения вершины трещины. Рассмотрены два предельных случая и один промежуточный случай: предельные случаи соответствуют началу и завершению процесса продвижения вершины трещины, а промежуточный случай — установившейся стадии процесса.
Ключевые слова: хрупкое разрушение, квазихрупкое разрушение, усталость, кривые Пэриса
Diagrams of quasibrittle fatigue fracture (two-frequency loading)
V.M. Komev
Lavrentiev Institute of Hydrodynamics SB RAS, Novosibirsk, 630090, Russia
The paper considers jump-like propagation of a crack in quasibrittle materials under stationary and nonstationary cyclic loading. In the nonstationary mode, particular attention is given to two-frequency load. For analyzing the crack propagation, it is proposed to use a diagram of quasibrittle fracture of solids under cyclic loading. One of the curves of the proposed diagram resembles the Kitagawa-Takahashi diagram. Explicit relations for the dimensionless average velocity of an abruptly extending crack per loading cycle are derived. The derived relations for average velocity can be considered as structural formulae for plotting Paris curves. The advance of the crack tip for the Paris curves is analyzed in detail covering two limiting cases and one intermediate case: the onset and the end of the process and its steady-state development, respectively.
Keywords: brittle, quasibrittle fracture, fatigue, Paris curves
1. Введение
Теоретические предсказания и практические рекомендации по малоцикловой усталости дисков компрессоров авиадвигателей могут отличаться в несколько раз [1]. В данной работе предпринята попытка сблизить указанные оценки, когда рассматриваются нестационарные пульсирующие режимы нагружения, причем особое внимание уделяется двухчастотному нагруже-нию. Пусть одна из частот двухчастотного режима нагружения соответствует режиму «взлет-посадка», а вторая частота этого режима описывает либо вынужденные стационарные колебания лопаток компрессора,
либо колебания лопаток компрессора при попадании птицы. Колебания лопаток передаются на диск компрессора через специальные крепления. В окрестности этих креплений реализуется двухчастотный пульсирующий режим нагружения материала. Простейший вариант двухчастотного нагружения приведен на рис. 1, где Р — суммарное нагружение при двухчастотном пульсирующем режиме; р — нагружение, соответствующее режиму «взлет-посадка»; Р2 — нагружение от колебания лопаток. Справа на рис. 1 приведено преобразование исходного режима нагружения к двум пульсирующим режимам.
О Корнев В.М., 2012
Допустим, что в окрестности крепления лопаток к диску имеется внутренний трещиноподобный дефект. Получим оценки числа циклов при двухчастотном нагружении, когда режиму «взлет-посадка» соответствует такой уровень нагружения диска компрессора, при котором реализуется малоцикловая усталость. Для получения оценок числа циклов при двухчастотном нагружении используется диаграмма квазихрупкого разрушения при малоцикловой усталости [2], построенная для одночастотного пульсирующего нагружения на макроуровне. Использование указанной диаграммы оправдано, когда наложение второго режима нагруже-ния на первый не выводит комбинированный режим нагружения за рамки малоцикловой усталости.
2. Диаграммы квазихрупкого разрушения при однократном нагружении образца конечной ширины
Примем простейшую аппроксимацию а-е-диа-граммы материала, когда эта диаграмма аппроксимируется двухзвенной ломаной, описывающей поведение упругопластических материалов. Параметры этой аппроксимации Е, а^ ад, е0, е1: Е — модуль упругости; а у — предел текучести или «теоретическая» прочность гранулированного материала; ад — постоянные напряжения, действующие согласно модифицированной модели Леонова-Панасюка-Дагдейла (ау ^ ад) [3, 4]; е0 — максимальное упругое удлинение материала (аY = Е е0); е1 — максимальное удлинение материала. Пусть для гранулированного материала с регулярной структурой г — диаметр зерна.
Подход Нейбера-Новожилова [5, 6] позволяет для сред со структурой использовать решения, имеющие сингулярную составляющую, поскольку бесконечные напряжения в вершине модельной трещины, не допускаемые классическими континуальными критериями прочности, не противоречат дискретным критериям, если сингулярная составляющая решения имеет интегрируемую особенность.
Ради простоты рассматривается внутренняя трещина нормального отрыва. Диаграммы квазихрупкого
разрушения тел с краевыми трещинами получены в [7]. Пусть внутренняя плоская трещина нормального отрыва распространяется прямолинейно. Кроме реальной внутренней прямолинейной трещины-разреза длиной 210 введем в рассмотрение модельные трещины-разрезы длиной 21 = 210 + 2Д. Каждая из зон предраз-рушения Д расположена на продолжении реальной трещины (2/, Д — длины модельных трещин и зон предразрушения).
Таким образом, наша задача об усталостном разрушении имеет два линейных масштаба: если диаметр зерна г определяется структурой материала, то второй линейный масштаб вырабатывается самой системой. Этим вторым линейным масштабом при малоцикловой усталости являются длины зон предразрушения Д, которые изменяются в соответствии с тем, как меняются длина реальной трещины, которая подрастает скачками, и интенсивность нагружения при двухчастотном режиме. За изменением второго линейного масштаба легко проследить, анализируя изменение расстояний между гребнями усталостных бороздок на изломе детали. Подчеркнем, что при однократном нагружении критическая длина зоны предразрушения Д* — вполне определенный параметр [8].
При построении диаграмм квазихрупкого разрушения в условиях малоцикловой усталости в [2] используются достаточные критерии разрушения, когда рассматриваются трещины нормального отрыва:
1 пг
— \ау (х,0^ <а¥, (1)
кг 0
2у(х,0) <5, -Д < х < 0. (2)
Здесь ау (х, 0) — нормальные напряжения на продолжении трещин; Оху — прямоугольные системы координат, ориентированные относительно правых частей трещин (начало координат совпадает с вершиной модельной трещины в модифицированной модели Лео-нова-Панасюка-Дагдейла [2, 8]); п, к — целые числа (п > к, 1 < п < 4); пг—интервал осреднения; (п - к)/п — коэффициент поврежденности исходного материала на интервале осреднения (1/2 < (п - к)/п < 1); 2у = = 2у(х, 0) — раскрытие трещин (х < 0).
Рис. 1. Двухчастотный режим нагружения диска, соответствующий вращающемуся диску при наличии колебаний лопаток
В модифицированной модели Леонова-Панасюка-Дагдейла [3, 4] поле нормальных напряжений сту (х, 0) на продолжении модельных трещин х > 0 можно представить в виде суммы двух слагаемых:
ст у (х, 0) = К:/(2 пх)1/2 + 0(1), (3)
К1 = К1ш + К1Д , К1ш > 0, К!Д < 0,
где Кг = К (I, Д) — суммарные коэффициенты интенсивности напряжений в вершинах модельных трещин; К1ш — коэффициенты интенсивности напряжений, порождаемые напряжениями стш, заданными на бесконечности; К1Д — коэффициенты интенсивности напряжений, порождаемые постоянными напряжениями стД. Первое и второе слагаемые в соотношении (3) — сингулярная и гладкая части решения соответственно.
При описании зон предразрушения используется [2, 8] первый класс решений, соответствующий хрупкому и квазихрупкому типам разрушения, когда в соотношении (3) суммарные коэффициенты интенсивности напряжений положительны, т.е.
КI > 0. (4)
Для критических параметров неравенства (1) и (2) достаточного критерия превращаются в равенства, причем первое равенство контролирует достижение напряжениями на продолжении модельной трещины после осреднения предела текучести или «теоретической» прочности сту, а второе описывает затупление в вершине реальной трещины.
Первое равенство достаточного критерия (1), (2) представляет типичный силовой критерий разрушения [9], а второе равенство этого критерия — деформационный критерий [9]. Таким образом, достаточный критерий разрушения (1), (2) одновременно учитывает силовой и деформационный критерии разрушения в характерных точках зоны предразрушения. При таком подходе [2, 8] к описанию процесса разрушения отпадает необходимость обсуждать преимущества и недостатки силовых или деформационных критериев [10-12]. По Новожилову [6] предлагаемый критерий (1), (2) является достаточным, он сочетает преимущества силового и деформационного критериев и частично нивелирует их недостатки. Выбранные характерные точки зоны предразрушения хорошо приспособлены к описанию скачкообразного продвижения вершин трещин и накопления повреждений в материале зоны предразрушения при усталости [2].
Простейшее аналитическое представление поля напряжений сту (х,0) на продолжении модельной трещины имеет вид [8], когда выполнено неравенство (4):
К ™ + К Д
s y ( x,0) = -
\l2nx
ID
+ S =
ст^ IT + 2CTWÄ+CT^.
•Jx v 2 njx
Введем в рассмотрение критические напряжения,
рассчитанные по необходимому ст° и достаточному
si критериям разрушения. Эти критические напряжено *
ния стш и стш представляют нижние и верхние оценки напряжений sœ. Наиболее интересная область нагру-жения с точки зрения изучения малоцикловой усталости есть область ст° < sœ < si. Критическая кривая разрушения ст° = si (l0, sY ) описывает разрушение хрупких материалов, когда зона предразрушения отсутствует, т.е. Д = 0, а sY — «теоретическая» прочность хрупкого материала. Критическая кривая разрушения sL = sL[l*, Sy, (Si - So)/So, sa/sy] описывает разрушение квазихрупких материалов, когда Д */10 << 1, при-
Д*
— критическая длина зоны предразрушения; 21 * = 2l0 + 2D* — критическая длина макротрещины; sY — предел текучести упругопластического материала.
После решения системы (1), (2) с учетом неравенства (4) и соотношения (5) получим критические параметры разрушения для квазихрупких материалов:
SL = SL [l *, SY,(S1 — S0 VS0 , Sд/SY], * * * (6)
Д* = Д*[Г, Sy,(Si —SoVSo]. Равенство (1) — необходимый критерий разрушения при Д = 0, используется для отыскания критического параметра
s°=s° (lo, Sy). (7)
Диаграммами разрушения для квазихрупких материалов являются пары кривых
o o п \ Si = Si ('0, SyX
sL = sL [l *, s y , (Si —SoV So, Sd/Sy]. Аналитические представления критических параметров разрушения для квазихрупких материалов с учетом представления решения (5) имеют вид:
s
(
1 — -
5 sD s — S
Л
8n s
2 С
д* =
S — S
y o / \2 ( * Л2
S„
2l * \fn + n r k k
s
l , 2l = 2lo + 2Д*
Выражения (8) и (9) имеют смысл, если
1 — A SD s1 — So 8п ss
> o.
(8)
(9)
(1o)
Y bo
Неравенство (1o) — ограничение, которое выполняется только для хрупких и квазихрупких материалов. Аналитическое представление кривой si = si (lo, sY ) получается из соотношения (8) при Д* ^ o, точнее при s1 ^ so (для однократного нагружения приняты обозначения lo = lo): _o ( / 2lo yfn
sL
n ■ + —
k k
Y
(11)
Очевидно, что si < sL при s1 > so, см. (8), (11).
s
Y
Введем в рассмотрение критические коэффициенты интенсивности напряжений, полученные по необходимому К° и достаточному К*м критериям разрушения с учетом представления К1м = УКХ = Уа^^/ЛТ" в полосе конечной ширины:
К0 = У ° К0 = У0 а^Тч,
Г- (12)
К*, = У * К* = У *а^л/л/*.
Здесь У 0 = У °(/0/ м), У * = У * (/*/м) — коэффициенты; м — ширина полосы.
Критические коэффициенты интенсивности напряжений К° и из соотношений (12) используем для получения оценки влияния ширины образца на вид диаграмм разрушения. Величины коэффициента У0 = = У0 (/0 / м) для разных соотношений /0 / м и условий нагружения можно найти в стандартных справочниках [13, 14]. Целесообразно обратить внимание на приближенные аналитические соотношения для коэффициента У0 = У0 (/0/м) из [13], которые допускают предельный переход для трещин конечной длины /0:
У0 (/0/м) ^ 1 при /0/м ^ 0,
У0 (/0/м )^ад при 2 /0/ м ^ 1.
Предельные значения коэффициента У0 = У0 (/0/м) из соотношений (13) таковы, что критические напряжения а^=а^[/*, аY,(e1 -80)/80, ад/аY, п, к], а° = = а^ (/0, а^ п, к) из (8) и (11) можно использовать в практических расчетах, когда выполняется сильное неравенство /0/м ^ 1. Для полосы ко-нечой ширины критические напряжения =
= [[* а, (81 - 80 V80 , аД /аY , м n, к] и =
= (/0, аY, м, п, к) отыскиваются с использованием представления (5), см. [6], и имеют вид:
(13)
s
0
UY
w Ф ад,
s*
Y 0 ( у w )
2l0 4n
r к
n
1
к 1 - 2 l
'w
1-
5 sD
8n Sv
Y * (l */ w)>
21 * Vn к
1
к 1 - 21 *
'w
-1
w Ф ад.
(15)
Последние сомножители 1/(1 - 2 /0 /м), 1/(1 - 2 /*/м) в скобках соотношений (14), (15) отражают влияние изменения работающего сечения образца на критические параметры , . Поскольку для величин коэффициента У * = У * (/*/ м) отсутствуют общепринятые данные в справочниках [13, 14], используем приближенное соотношение У* « У0 для квазихрупких материалов, когда Д*//0 << 1.
Результаты расчетов по влиянию конечности ширины образца на критические кривые разрушения
= S
= s,
[l* Sy, (81 - So V80, sa/sy , w, n, к], (lo, sY, w, n, к)
0,
на плоскости «длина трещины - внешняя нагрузка» представлены на рис. 2 в двойных логарифмических координатах для w/r = 500. На этом рисунке сплошная кривая 1 построена по соотношению (11) для s^ пунктирная кривая 4 — по соотношению (14) для s| сплошные кривые 2 и 3 — по соотношению (8) для s^ соответственно при (Sj -e0)/e0 = 1.5, 2.5 (sa/sy = 1), пунктирные кривые 5 и 6 — по (15) для s^w соответственно при (Sj -e0)/e0 = 1.5, 2.5 (sa/sy = 1). При построении кривых принято, что отсутствуют повреждения материала, т.е. n = к = 1.
Пары кривых 1-2, 1-3 и 4-5, 4-6 на рис. 2 представляют диаграммы квазихрупкого разрушения полос бесконечной и конечной ширины соответственно, причем пары кривых 1-2, 4-5 и 1-3, 4-6 построены при (Sj -e0)/e0 = 1.5 и 2.5 (sa/sy = 1). Влияние параметра sa/sy на указанные диаграммы, когда sa/sy ф 1, легко анализируется. Это влияние сводится к исследованию влияния комбинированного параметра (s a /s Y) х х (sj - s0)/e0 в соотношениях (8), (15).
Отклонение кривых 4, 5, 6 для пластин конечной ширины от соответствующих кривых 1, 2, 3 наблюдается в интервале 21/г > 100 при w/r = 500. Соотношения (8), (11), (14), (15) для образцов с внутренней трещиной позволяют указать интервал относительных длин трещин 2l0/w < 0.2, в котором можно пренебречь влиянием ширины образца на критические параметры s^, s^.
Получим из соотношений (8) и (11) оценки для предельных случаев:
lim s^ = lim s^ = ksY/n, (16)
(14)
10 ^0
s_2
s Y
*
s^
l
к 4n
w к
s Y yjn
7*
при l ^
^ при l0
1-
5 sa
Л-1
8n sY
(17)
ад.
1.00
0.50
>-ь
b 0.20
>-ь
:0.10
0.05
........1 ........1 1 ■ ■ ч \ \ \ \\\ \ \\\ \ .......
0.1
1
10
100
21/r
Рис. 2. Диаграммы квазихрупкого разрушения, учитывающие ширину образца м/г = 500, при однократном нагружении (сплошные и пунктирные кривые для образцов бесконечной и конечной ширины)
Когда трещина отсутствует, точнее /0 ^ 0 или /* ^ 0 в соотношениях (16), прочность тела прямо пропорциональна напряжениям сту с учетом повреж-денности материала на интервале осреднения как для достаточного (8), так и для необходимого (11) критериев разрушения. Для длинных трещин, точнее 10 ^ да или / ^-да в соотношениях (17), прочность тела прямо пропорциональна напряжениям сту с учетом повреж-денности материала на интервале осреднения и обратно пропорциональна корню квадратному из относительных длин трещин 2 /0/г или 2/*/г для необходимого (11) или достаточного (8) критериев разрушения. Для достаточного критерия надо принять во внимание поправочный коэффициент (стд/ст-у)(е1 -е0)/е0, зависящий от характеристик упругопластического материала зоны предразрушения. Оценки (16), (17) согласуются с видом кривых 1, 2, 3 на рис. 2.
Критические коэффициенты интенсивности напряжений К{° или К* для пластины бесконечной ширины не являются постоянными материала, так как вычисляются по известным параметрам материала со структурой. В рамках построенной модели имеются всего два таких основных параметра для коэффициента интенсивности напряжений К:0 = К:0(ст^ 2 /0/г, п, к), а коэффициент интенсивности напряжений
К * = К* [ст Y , 2/*/Г , (б! -80 )/80 , Стд /ст Y , п, к ] включает уже четыре основных параметра, п и к — второстепенные параметры.
Принимая во внимание стандартные представления коэффициентов интенсивности напряжений К° = = стда^/Ц,, К* = ст^л/п/ *, легко получить предельные соотношения для критических коэффициентов интенсивности напряжений К°, К*, см. оценки (17):
nr
\— пРи l0 ^ ^
K0 к s KI ~~Т SY sjn
K'
(18)
при l ^ да.
Только для достаточно длинных трещин критические коэффициенты интенсивности напряжений можно считать постоянными, так как KI и const, K' и const. Причем эти постоянные зависят от поврежденности материала на интервале осреднения, а последняя постоянная зависит еще и от свойств упругопластического материала зоны предразрушения.
Первые оценки из соотношений (17) и (18) для ct^/cty и K° при l0 ^да соответствуют линейной механике разрушения. Вторые оценки из соотношений (17) и (18) для <^/cty и K' при l* ^да очень напоминают соотношения линейной механики разрушения, если принять во внимание влияние поправочного коэффициента (ctd/cty) (Sj -s0)/80 для упругопла-
стического материала, имеющего квазихрупкий тип разрушения.
Предлагаемая диаграмма квазихрупкого разрушения рис. 2 при однократном нагружении состоит из двух кривых, полученных по необходимому и достаточному критериям разрушения, причем всегда критические напряжения для трещин конечной длины /* > /0 > 0 достаточного критерия стда превосходят критические напряжения необходимого критерия ст^-,, т.е. стда >ст|^. Часть энергии теряется при пластическом деформировании материала зоны предразрушения. Построение диаграмм разрушения обсуждается в монографиях [9, рис. 3.12 на с. 76; 15, рис. 27.1-27.4 на с. 233-239]. Особенно много внимания уделено в [15] построению кривых разрушения исходной трещины /0 при последовательном догружении, однако в предлагаемых исследованиях используется другая параметризация, отличная от параметризации в аналогичной задаче, см. [15, с. 232].
В предлагаемой модификации модели Леонова-Панасюка-Дагдейла кривая разрушения для исходной трещины /0 на плоскости (/, стда) есть составная кривая с предельной (критической) точкой, когда стда < стда. Кривая разрушения состоит из вертикальной прямой / = /0, когда стда < стда, кривой, соответствующей соотношению (9), когда стда<стда <стда (е0 < е <е1), и заканчивается критической точкой, когда стда = стда (Д = Д *). Максимальная длина зоны предразрушения определяется величиной параметра неупругого деформирования материала е0 < е < е1 (0 < Д < Д *); этот параметр в задачах усталости есть переменная величина из-за охрупчивания материала зоны предразрушения. Используем естественное параметрическое представление кривой разрушения / = /(е), стда = стда (е) на плоскости (/, стда), где е—текущее значение параметра неупругого удлинения, поскольку е0 < е < е1. Это параметрическое представление второго участка кривой разрушения при квазихрупком приближении Д//0 << 1 имеет простой вид:
2 (
2l = 2l0 + 2D,
д=4
8-8
2
l
Vu Y
0 < д < д '
(
1 -
(19)
5 СТЛ 8-8о
8п ct-v
4n n
— + —
к к
ст < ст < ст
да — да — да •
На рис. 3 изображены диаграммы квазихрупкого разрушения, представляющие из себя две пары кривых 1-2 или 1-3. При построении кривой 1 для стда/сту использовалась формула (11), по формулам (8) построены кривые 2 и 3 для ст^сту при (е1 -е0)/е0 = 2.5 и 3.0 соответственно. Диаграммы квазихрупкого
2
Рис. 3. Диаграммы квазихрупкого разрушения и кривые разрушения для разных длин трещин, штриховые и штрихпунктирных кривые
разрушения на рис. 3 дополнены вторыми участками кривых разрушения для конкретных длин трещин 2l0/ r, см. соотношения (19) (вертикальные участки кривых разрушения для < а^ не изображены). Длины исходных трещин 2l0/r подобраны так, чтобы получилась «спина динозавра» по Р. Томсону [8, 16, 17]. На рис. 3 эти кривые разрушения изображены в виде сплошных кривых, упирающихся в штриховую или штрихпунктирную кривую, последние кривые соответствуют параметрам (e1 -e0 )/e0 = 2.5 и 3.0. Чем больше параметр неупругого удлинения материала, тем больше значения критических длин зон предразруше-ния. Непрерывные кривые разрушения из [16] имеют множественные локальные минимумы и максимумы и описывают решеточный захват. В рассматриваемом случае точки экстремумов непрерывных кривых разрушения [16] превращаются в точки разрыва первого рода из-за вида аппроксимации стандартной а - e-диаграммы упругопластического материала. Точки разрыва первого рода кривых разрушения являются причиной скачкообразного продвижения вершин трещин при усталости. За счет материала зоны предразрушения в окрестности вершины трещины формируется композиционная структура с переменными свойствами при усталости. Материал зоны предразрушения охрупчивается. Основным параметром, описывающим этот процесс, является параметр неупругого удлинения e (e0 <e<e1); при охрупчивании этот параметр уменьшается.
3. Диаграмма квазихрупкого разрушения при усталости (двухчастотный режим нагружения)
Исходная диаграмма квазихрупкого разрушения при усталости совпадает с диаграммой квазихрупкого разрушения при однократном нагружении, если рассматривается одночастотный пульсирующий режим нагружения с постоянной амплитудой [2]. Ниже, как и в работе [2], рассматривается циклически стабильный материал, т.е. предел текучести не зависит от номера цикла а Y = const. Одна из кривых диаграммы квазихрупкого разрушения преобразуется с учетом накопле-
ния повреждений в зоне предразрушения для последующих циклов нагружения до первого скачка вершины трещины, так как материал этой зоны охрупчивается. После скачка вершины трещины все повторяется, но для новой длины трещины.
Сравним предлагаемую диаграмму квазихрупкого разрушения при малоцикловой усталости для одно-частотного стационарного пульсирующего режима нагружения [2] с диаграммой Китагава-Такахаши [18, 19]. Приведенная на рис. 4 диаграмма Китагава-Такахаши заимствована из [19, рис. 2 на с. 749]. Диаграмма Китагава-Такахаши [18] состоит из одной кривой, аналитическое представление которой можно найти на рис. 4; эта кривая разделяет область усталостного разрушения и область, где такое разрушение отсутствует. Для очень коротких трещин эта кривая непосредственно связана с пределом усталости по напряжениям, а для длинных трещин включает пороговое значение коэффициента интенсивности напряжений при усталости [20]. В диаграмме Китагава-Такахаши [18] не прослеживается какая-либо связь с упругопластическими свойствами материала (как исходного материала, так и модифицированного материала зоны предразрушения). Все рассуждения по построению диаграммы Китагава-Та-кахаши [18] содержат некоторое противоречие, связанное с одновременным использованием критического коэффициента интенсивности напряжений и предела усталости по напряжениям.
Предлагаемая диаграмма квазихрупкого разрушения при малоцикловой усталости [2] не связана с использованием коэффициентов интенсивности напряжений. Ее построение непосредственно зависит как от упругопластических свойств материала, так и от длины трещины. Эта диаграмма, состоящая из пары кривых, разбивает плоскость «длина трещины - амплитуда внешней нагрузки» (21/r, /аY ) на три подобласти: 1) левее и ниже кривой а^/аY расположена область, где усталостное разрушение отсутствует; 2) правее и выше кривой а^/аY расположена область
Рис. 4. Диаграмма Китагава-Такахаши [19]
разрушения при однократном нагружении; 3) между кривыми ct°/cty и ст^/сту расположена область скачкообразного продвижения вершины трещины из-за охрупчивания материала зоны предразрушения при циклическом нагружении. При повторных нагружениях кривая ст*/сту преобразуется, при преобразованиях учитывается охрупчивание материалы зоны предразрушения. Предлагаемая диаграмма усталости в двойных логарифмических координатах ограничена горизонтальной и наклонной прямыми, см. кривые 1, 2, 3 на рис. 2 для пластин бесконечной ширины. Расположение горизонтальной прямой связано с пределом текучести сту упругопластического материала, угол наклона и расположение наклонной прямой определяются комплексным параметром (стд/сту )(ех -e0)/e0 неупругого деформирования материала. Предлагаемая диаграмма имеет два пороговых значения ст°/сту и af/aY , см. соотношения (8) и (11).
Предлагаемая модифицированная модель Леонова-Панасюка-Дагдейла, примененная к усталостному разрушению, очень напоминает модель с гипотетическими элементами усталости [21, рис. 1], расположенными в вершине реальной трещины. Какая-либо связь длины этих элементов с размером пластической зоны в вершине реальной трещины отсутствует, но подчеркивается, что пластическая зона связана с пиковыми величинами нагрузок и коэффициентов интенсивности напряжений [21, рис. 2].
Остановимся более подробно на двухчастотном режиме нагружения, который иллюстрирует правая часть рис. 1. Пусть двухчастотный режим нагружения диска соответствует раскрученному диску при наличии колебаний лопаток. Предположим, что частоты первого («взлет-посадка») режима f отличаются от второго на несколько порядков, т.е. f << f2. Пусть для частот справедливо соотношение
Л/ f =w. (20)
Рассмотрим такое двухчастотное нагружение, при котором второй режим накладывается на первый основной режим нагружения (второй режим нагружения не может существовать без первого режима). Таким образом, как и в [21], выбирается конкретная программа циклического нагружения. При осуществлении программы циклического нагружения, приведенной на рис. 1, имеет место полная разгрузка после посадки. Эта программа отличается от программы нагружения из [21, рис. 2(5)], где изучается единичная перегрузка.
Вместо ограничения
ст°/сту <ст+/сту <ст^/сту , (21)
которое справедливо для одночастотного стационарного режима нагружения [2], когда амплитуда одночастотно-го пульсирующего нагружения постоянна (af = const), рассмотрим другие ограничения:
af/aY <af/aY <af/aY ,
ct°/CTy <CT++/CTY <CT^/CTY , (22)
soo = af + 2af2 ' af = af 1 — af2 .
Здесь af — суммарная амплитуда колебаний [21]; af j = const и af 2 = const — амплитуды, соответствующие первому и второму стационарным режимам; af — амплитуда с учетом разгрузки от второго колебательного режима, причем при отсутствии колебаний лопаток имеем
lim af = lim af = afP
af2 ^0 af 2 ^0
Подчеркнем, что второе ограничение и последнее равенство из (22) соответствуют пульсирующим режимам нагружения с учетом преобразований правой части рис. 1. При выполнении ограничений (22) имеем малоцикловую усталость как при одночастотном основном режиме нагружения с амплитудой af, так и при двухчастотном режиме нагружения с амплитудой af. Двухчастотный режим нагружения рассматривается как своеобразная перегрузка на фоне постоянной нагрузки af, какие-либо дополнительные ограничения на амплитуду af 2 не налагаются. Второе ограничение из (22) очень напоминает соотношение (21). Очевидно, что при выполнении соотношения af/aY >af/aY имеет место катастрофическое развитие событий: происходит разрушение диска, например при попадании в двигатель крупной птицы.
Рассмотрим двухчастотное нагружение, при котором N j — номер цикла в режиме «взлет-посадка»; N2 j — номер цикла в режиме стационарных колебаний лопаток. При скачкообразном продвижении вершины трещины меняется ее длина, поэтому надо в номерах циклов N j и N2 j ввести дополнительный индекс j, такой что j = 1, 2, ..., j* — число скачков вершины трещины; j* — критическое число скачков, когда объект разделяется на части, причем образец имеет исходную трещину длиной l0 до первого цикла нагружения.
Исследуем изменения подобласти диаграммы квазихрупкого разрушения для циклов нагружения N2 j между скачками вершины трещины, когда принимается во внимание охрупчивание материала зоны предразру-шения. Сначала рассматривается простейший случай, ср. (20), когда
N2 < ю. (23)
Здесь N2 — критическое число циклов колебаний лопаток, т.е. до первой посадки диск компрессора теряет целостность при наличии трещины. Изменения в подобласти диаграммы квазихрупкого разрушения для циклического нагружения, соответствующего режиму стационарных колебаний лопаток, с учетом накопления повреждений описываются неравенствами
*( s)
0 S
s Y s Y s Y
1 < s, < N2j -1, j = 1,2,...,/-1,
<s:(s) <...<s
•(2)
<s
•(1) - ,
(24)
^ЭД > ^ЭД
Sj = N2 j>
у = 1,2,..., у ,
когда критическое число циклов (колебаний лопаток) нагружения N2 подсчитывается так:
N2= 1+1
j -1
N
N > 2,
Л 1 2, j 2, j
j = 1,2,..., j *-1, N2, j.= 1.
(25)
Здесь аЭДs^ — критическая нагрузка, полученная с использованием достаточного критерия квазихрупкого разрушения на s-м цикле нагружения до j = 1, 2, ..., j * -1 скачка; s, — номер цикла между j - 1 и j скачками, Sj = 1 соответствует состоянию поставки материала после каждого скачка вершины трещины; N2 , — число (группа) циклов колебаний лопаток между j - 1 и j скачками вершины трещины; N2 — критическое число циклов нагружения, после j скачков длина исходной трещины l0 изменяется так, что система разрушается за один цикл N2 * = 1. Подчеркнем, что, как правило, N2 >> j*, так как N21 > N22 > ... > N > 2 и N * = 1. Таким образом, критическое (общее) число циклов нагружения N2 разбивается на j * групп, в каждой из которых N2j циклов (j = 1, 2,..., j *).
Первые три строки соотношений (24) описывают накопление повреждений в зоне предразрушения. В двух последних строках соотношений (24) записано условие, при котором происходит скачок вершины трещины. Напомним, что аЭД — амплитуда напряжений в каждом цикле, а материал образца рассматривается как своеобразный композит. В соотношениях (24) используется только одно значение аЭД для критической нагрузки по необходимому критерию (11), это соответствует рассмотрению циклически стабильных материалов, когда s Y = const. Соотношения (24) описывают перестройку диаграмм квазихрупкого разрушения композиционной структуры с переменными свойствами при циклическом нагружении. Имеет место кумулятивный эффект накопления повреждений в зоне предразрушения.
Остановимся на процессе продвижения вершины трещины скачками. В работах [2, 22, 23] была высказана гипотеза об останове макротрещины раскалывания при стационарном пульсирующем нагружении, эта гипотеза была подтверждена экспериментально [24-26]. Исходная острая трещина длиной 2lj послеj-го скачка распространяется только по охрупченному материалу зоны предразрушения на A j (A j — длина зоны пред-
разрушения перед вершиной реальной трещины длиной 2lj-1). После останова перед новой вершиной реальной трещины расположен исходный материал, имеющий квазихрупкий тип разрушения.
Длины модельных 2lj трещин отличаются от длин реальных 2lj-1 трещин на две длины зон предразрушения 2Aj после каждого скачка:
2lj = 2lj-1 + 2Aj, j = 1, 2, ..., j* (26)
причем 2l0 — исходная длина реальной трещины;
2l .* = 2lS* — критическая длина модельной трещины
j S / при циклическом нагружении, когда Y = const.
На рис. 5 представлена схема, поясняющая скачкообразное продвижение вершин трещин при двухчастот-ном нагружении, см. рис. 1 и соотношения (24)-(26). Исходная диаграмма квазихрупкого разрушения при усталости изображена на рис. 5 в двойных логарифмических координатах: кривая 1 соответствует напряжениям аЭД/sY из (11), кривая 2 — напряжениям аЭД/sY из (12), горизонтальная прямая аЭД/sY = const со стрелками отражает продвижение вершин реальных трещин. На кривой 1 отмечены точки (2lS0/r, аЭД /sY), (2l+°/r, s++/sY), соответствующие длинам трещин
2lS0 и 2l+0 для заданных уровней нагружения аЭД /<
s Y
и s+^/sY. На кривой 2 отмечена критическая точка (2lS*/r, аЭД/sY), в которой образец разделяется на части 2Г* = 2lS*. Чтобы выполнить оба ограничения из (22) исходная длина 2l0 реальной трещины должна
+0 / / S* /
подчиняться ограничениям 2l /г < 2l0/ r < 2l /г .На горизонтальной прямой аЭД /sY множественные стрелки от точки (2l0/r, s^;/sY) к точке (2lS*/r, s;^/sY) отражают скачкообразное продвижение вершины реальной трещины.
Из соотношений (8), (9) по заданному уровню нагружения аЭД/а Y = const оцениваются критическая длина
7 / S*
2l /r трещины и критическая длина A зоны предразрушения для последнего цикла нагружения:
-1
sY
1 -
5 s A ei e0
8n sY
2l
S*
sfn П -+ —
k k
(27)
Рис. 5. Скачкообразное продвижение вершин трещин на диаграмме квазихрупкого разрушения при усталости (образец бесконечной ширины)
DS* = L-29
2 Л
\2 Л
S*
(28)
Подчеркнем, что полученные оценки (27), (28) не зависят от критического числа циклов нагружения N.
Второе ограничение из (24) выполняется для последнего цикла в каждой группе циклов. Изучим накопление повреждений в материале зон предразрушения на каждом цикле нагружения, принимая во внимание к какой из групп циклов] = 1, 2, ..., ]* рассматриваемый цикл принадлежит, см. соотношение (26) для длин модельных трещин 21 ^. Имеет место кумулятивный эффект накопления повреждений в зоне предразрушения. Это накопление повреждений можно описать линейным или нелинейным законами суммирования повреждений.
На каждом цикле (1 < < Ы2-1, ] = 1, 2, ..., ]* -1) для реальной трещины длиной 21 ^г отыскиваются относительные (е^ -е0 )/е0 и -е0 )/е0 неупругие удлинения материала в вершинах реальных трещин при заданных уровнях нагружения ст^ /ст
3/cty в рамках модифицированной модели Ле-онова-Панасюка-Дагдейла. Для уровня нагружения
(29)
j . (30)
(31)
V /
Выражения (29), (30) имеют смысл, если
1 -
5 ст
Д j e0
8п sY sn
> 0, j = 1, 2,
j.
Для уровня нагружения ст^ /ст ■ получаются такие же соотношения, см. (29)-(31), если в последних верхний значок Е заменить на значок +. Очевидно, что
де >д+. (32)
Поэтому в соотношении (26) длина скачка вершины трещины равна Д^ = Д;-: разрушение имеет место в охрупченном материале. Соотношение (32) влияет на выбор длины скачка вершины трещины, далее величина меньшей из длин зон предразрушения Д+ нигде не используется.
Переходим к оценке накопления повреждений в материале зоны предразрушения. По известным значени-
ям ст^/ст
/с
2L
r и ctd/cty или стт/ст
/с
21
j-1
Ir
стД /ст■ из соотношения (29) отыскиваются неупругие удлинения (е^^ -е0) или (е+^-1 -е0), а из соотношения (30) подсчитываются длины Д^ зон предразрушения. Разность между указанными неупругими удлинениями (е^Е-1 - е+-1) провоцирует накопление повреждений ма-
териала зоны предразрушения. Л.Ф. Коффин [27, 28] и его последователи экспериментально установили зависимость между неупругой деформацией материала в цикле и числом циклов. При другом подходе для учета накопления повреждений в материале зоны предраз-рушения [29] вводится понятие «критическая величина диссипируемой работы», при которой материал разрушается. Соотношения первого [27, 28] и второго [29] подходов могут быть записаны единообразно. Число циклов колебаний лопаток Ы2 ^ между] - 1 и] скачками вершины трещины подсчитывается как
(N2, j )C =-
j = 1, 2, ..., j*-1.
(33)
j Sj-1
Здесь 0.2 < C < 1 — некоторые постоянные (постоянные Коффина), численные значения которых зависят от свойств материала [27, с. 76-77]. Число циклов колебаний лопаток N2 j в соотношении (33) определяется отношением максимального неупругого удлинения материала к разности (размаху) неупругих удлинений при двух режимах нагружения с учетом постоянной Коффина. Таким образом, в уравнение Коффина входят параметры исходного материала (Sj - e0 ) и неупругие удлинения (sj^ -s++-1 ) при двух режимах нагружения, причем последние получены при решении задачи о деформированном состоянии зоны предразрушения, см. (29). Соотношение (33) при линейном C = 1 или нелинейном C < 1 суммировании повреждений описывает процесс охрупчивания материала в вершинах реальных трещин. Эти постоянные C подбираются по натурному эксперименту [27, 28]. В отличие от Коффина [28] в числителе соотношения (33) использовано максимальное пластическое удлинение материала (s1 - e0 ) при однократном нагружении, что согласуется с предлагаемой моделью. Соотношение (33) описывает процесс охрупчивания материала, когда реализуется частичная разгрузка. Напомним, что при полной разгрузке материала в знаменателе вместо s++-1 появляется e0, см. [2].
4. Продвижение вершины трещины при малоцикловом нагружении (кривые Пэриса)
В предыдущем разделе были получены оценки для числа циклов N2 . и длин скачков вершин трещины
S *
Д j при j = 1, 2, ..., j -1. При усталостном разрушении естественно получить оценки безразмерной средней
S
скорости V = (Д^/r )/n2
"j г ч " " 2 j продвижения вершины трещины за один цикл нагружения
V=
Dj/r.
5
N
2, j
( S -Sj-1 S0
-0
Л2 As
Sj-1 Sj-1
1
\V с
Y /
^ , j = 1, 2, ..., j*-1.
VCTY J Г
x
9
2
0
Эта средняя скорость V = (Ду /г)/^2 у достаточно легко измеряется в натурных экспериментах [30].
Переходим к анализу полученной оценки (34). Так как соотношение (34) содержит параметры ст«, /у-1 в нужной комбинации, то перепишем его в виде соотношения
Ду/г
( Е -еу-1 е0
N
2, у
2
(„Е
0
у/ с
(К^ -1 )2
псту 2 г
(35)
у = 1, 2, ..., У*-1,
в котором используется обозначение К^, у для амплитуды коэффициента интенсивности напряжений при пульсирующем нагружении квазихрупкого материала:
-Е
К12, у-1 = ст;
у-!, у = 1, 2, ..., у-1, (36)
когда длина 2/у-1 реальной трещины меняется скачками, как в (26).
Соотношениям (34), (35) можно придать другой вид, если воспользоваться следствием, полученным из (29):
ДЕ/г
N
2,
(-Е V /
у
8п 5
у-1
г
( \2+1/с (
ст
у/ с
2/, 1 л/П
V с
1 -
2/у-1 -Щ
Е
(37)
У = 1, 2, ..., у-1, ДЕ/г П (8пУ/С (ст. Г/С ( е 0 ^
N
2,
(К!,,-: )2
V Д У
V Ь1 0 У
ЛСТу г
Ст^л/г
1 с
V К12,у-1 К12,у-1 у
1 /П к а ■л/г
2 КЕ,-!
У = 1, 2, ..., ] -1.
Е
к ст^
(38)
Возможно и другое представление члена в последней квадратной скобки в соотношении (38) в виде:
1 -
п ст-
. л/г
2 К12,]-1
И 2/
У-1
Любое из представлений этой скобки не влияет на окончательный вывод: полученный закон продвижения вершины трещины (38) не сводится к классическому уравнению Пэриса.
Компактные (34), (35) или громоздкие (37), (38) соот-
12,]-1
= СТ„
ношения с учетом обозначений (36) и К (у = 1, 2, ..., у] -1) описывают осредненный процесс продвижения вершины реальной трещины при заданном двухчастотном режиме нагружения. Режим на-гружения оказывает существенное влияние на процесс продвижения вершины трещины скачками, этот режим описывается параметрами ст«, ст«, /у-1 в соотношении (37) и параметрами К12 у-1, К+ у-1, ст« или К12 у-1, К12у-1, /у-1 в соотношении (38). Двухчастотный режим нагружения отражает первая квадратная скобка. Второстепенные члены в (37), (38) собраны в последней квадратной скобке.
Упомянутые четыре соотношения имеют смысл для следующих реальных длин трещин
(39)
2/+0 < 2/,. < 2/... = 2/Е* у = 1, 2, ., у*-1
у у
при заданных уровнях нагружения ст«, ст« или для следующих амплитуд коэффициентов интенсивности напряжений (рис. 5)
'+0 <КЕ,у-1 =а.
Е 'п/у-1 <СТ«
тпр ,
(40)
у = 1, 2, ..., у*-1.
Полученные соотношения (3 7), (3 8) с учетом ограничений (39), (40) можно рассматривать как структурные формулы для построения кривых Пэриса [31] для заданного двухчастотного режима нагружения. По сравнению со структурной формулой (8) к двум основным параметрам добавились три новых параметра: два параметра нагружения ст«/^ , ст«/^ и постоянная Коффина С. Параметры нагружения связаны с двухчастотным режимом нагружения при малоцикловой усталости, а третий параметр описывает накопление повреждений в материале зоны предразрушения. Основными функциями, определяющими структуру кривых Пэриса, целесообразно признать сомножители (ст« /сту )2 /у^/г в соотношении (37) и (К^у^ )2/сту2г в соотношении (38), дополнительные сомножители в квадратных скобках каждого из выражений вносят существенные коррективы. При построении кривых Пэриса очень часто используются двойные логарифмические координаты.
Обратим внимание на то, что соотношение (38) кроме коэффициентов интенсивности напряжений К12 у-1, К12у-1 содержит параметр длины трещины /у-1. Только для достаточно длинных трещин, когда можно пренебречь членом (п/к)^г/2/у-1 во второй квадратной скобке, это соотношение более или менее соответствует классической кривой Пэриса [31].
Соотношения (37), (38) не описывают резкое уве-
личение средней скорости (Ду /г)/N2 у продвижения вершины трещины на последних циклах нагружения [27, 32]. Это увеличение связано с тем, что при натурных экспериментах испытываются образцы конечной ширины. Напомним, что соотношения (37), (38) полу-
9
2
X
чены для трещин в пластинах бесконечной ширины. Чтобы учесть влияние конечности ширины образца на процесс продвижения вершины трещины при усталости, надо сначала перестроить диаграмму разрушения [7], как на рис. 2, см. разд. 7. Приложение.
Так как для реальных длин 21 у (у = 1, 2, ..., у*-1) макротрещин выполняются строгие ограничения (39), то рассмотрим два предельных случая и один промежуточный случай. Каждый из изучаемых случаев соответствует характерным точкам на кривых Пэриса [31] для пластин бесконечной ширины. Предельные случаи соответствуют началу и завершению процесса продвижения вершины трещины, а промежуточный случай — установившейся стадии процесса: 1) начало процесса 21, — 21+0, но 21, > 21+0; 2) завершение про-
цесса 2/, - 2/= 21
•J J
S*
j
но 2/
, < 2/.. = 2/
i ]
■S*.
3) промежуточная стадия процесса 21,-. .
Проведем анализ для каждого случая [1]. 1. Для начала процесса продвижения вершины трещины имеем оценку: £ /,.
lim
А-
2/, ^2/+°+0 N
- = a Ф 0,
(41)
2, j
где a — малая величина (т.е. для начала процесса характерны малые скорости продвижения вершины трещины). Порог усталости отсутствует из-за двухчастотного режима нагружения, соответствующего ограничениям (22).
2. При завершении процесса продвижения вершины трещины имеем оценку (надо принять во внимание соотношения (27), (28) и то, что N * = 1):
lim
^2/S*
AS / r
2j
-о N
2, j
2
2
e -e
= lim
2/j ^2/S*-
Л2 ( S A2
A'
S*
о r
,S*
(42)
то есть при завершении процесса, когда образец бесконечной ширины разделяется на части, имеет место наибольшая скорость V = (А^r)jN2, продвижения вершины трещины за последний цикл нагружения. Для этой скорости получено аналитическое представление (42), которое зависит от трех параметров /S* /г, ст^ /стY , (Sj -s0)/s0, т.е. важна суммарная амплитуда ct^/cty нагружения, а не ее составляющие. Наибольшая скорость (AS /r)/n2 j в соотношении (38) не зависит от постоянной Коффина и параметра /cty , связанного с амплитудой колебаний лопаток.
3. При описании промежуточной стадии процесса, когда сформировалась достаточно длинная усталостная трещина 2/jjr » (/+0 + /S*)/r >> 1, скорость (AS /r)jN2 j продвижения вершины трещины можно представить в виде:
AS/r
N
2, j
(/ +0 +/S* )/2
1 с
sA
2+1/ С /
\V с
( s a
1 -
/ +0 + /S*
>(/+0 + / S*)4n
V с
'2 (/+0 + /S*
(43)
Соотношение (43) не соответствует простейшей аппроксимации кривой Пэриса [31, 33]. Во время всего процесса и на промежуточной стадии процесса первый член в квадратной скобке соотношений (38) и (43) отражает эффект «остаточного раскрытия вершины трещины» после снятия части нагрузки, см. «Residual crack opening stretch» в [33, с. S37-S40].
Итак, двухчастотный режим нагружения отражается на скорости V продвижения вершин трещин, эта скорость при двухчастотном режиме существенно отличается от аналогичных зависимостей [2] для одно-частотного режима нагружения.
Результаты расчета по соотношениям (37) и (46) приведены на рис. 6, а, б соответственно при линейном С = 1 и нелинейном С = 0.5 суммировании повреждений, когда k = n = 1, sY/sA = 1, (Sj -e0)/s0 = = 2.5, ct^/cty = 0.12, ct^/cty = 0.1. Кривые 1 и 2 каждого рисунка соответствуют пластинам бесконечной (ю ^ œ) и конечной (ю = 4 000) ширины. Отличие кривых рис. 6 связано с разным типом суммирования повреждений в материале зоны предразрушения: нарастание скорости V при нелинейном суммировании намного меньше, чем при линейном суммировании. Для установившейся стадии процесса скорости V при линейном и нелинейном суммировании повреждений различаются почти на порядок. Кривые 2 каждого из рисунков имеют характерную перевернутую форму S-образной кривой. Лавинообразное нарастание скорости V продвижения вершин трещин перед разделением образца на части связано с конечностью ширины этого образца. Имеется существенная зависимость начала третьего участка от ширины образца. Это нарастание скорости V не имеет никакого отношения к материалу образца, так как определяется геометрией образца. При двухчастотном нагружении отсутствует пороговое значение напряжений или пороговый коэффициент интенсивности напряжений, а при пульсирующем одночастотном нагружении такие пороговые значения существуют [2]. При реализации расчетов по формулам (38), (45) возникают определенные трудности, так как кроме коэффициентов интенсивности напряжений в них входят параметры, описывающие конкретный режим нагружения. Частично этим можно объяснить
r
Рис. 6. Кривые скорости роста длин трещин при линейном (а) и нелинейном (б) суммировании повреждений материала зоны предразрушения
возникновение полосы разброса экспериментальных данных, поскольку для каждого режима нагружения возникает своя кривая Пэриса.
5. Обсуждение
Скачкообразное продвижение вершин трещин при циклическом нагружении подробно обсуждается в монографии [34, с. 122-132], где кинетика пластической зоны в вершине трещины непосредственно связывается с особенностями роста усталостных трещин. В той же монографии много внимания уделяется вопросам автомодельности процесса продвижения вершин трещин. Естественным продолжением и обобщением подхода Баренблатта-Ботвиной является работа [35]. В заключении этой работы было обращено внимание на другие подходы, в которых предлагаются законы распространения трещин, отличные от классических кривых Пэриса, см., например, [36-38].
В предлагаемом подходе, см. (37), и в работе [2] процесс скачкообразного роста трещин описывается законами вида:
=
I ^
V
г а Л
-о
- С
а У
(44)
У
в которые входят текущая длина трещины 1/г
раметры нагрузки аа
У '
па-
параметры а-е-
диаграммы квазихрупкого материала (е1 -е0)/е0, аА/аУ и параметр Коффина С. Рассматриваемый закон (44) обобщает линейное (43) и нелинейное (44) представления из седьмого раздела работы [35]. Влияние конечности ширины образца на процесс разрушения удалось приближенно оценить, учитывая поправку к коэффициентам интенсивности напряжений, см. соотношения (45), (46) в приложении. Только после преобразований закона вида (44), см. соотношения (37), (46) данной работы, удалось получить кривую (38), (45), напоминающую классическую кривую Пэриса. Аналогичное утверждение справедливо и по результатам работы [2]. Заранее предугадать, в каком виде необходимо выбирать закон, описывающий рост трещин, не представляется возможным [35], так как сценарии нагружения при усталости могут быть весьма разноо-
бразными, а условия автомодельности при скачкообразном росте трещин, вообще говоря, не выполняются. Кроме того, надо рассматривать переходные процессы на стадии начала продвижения и на завершающей стадии перед разделением образца конечной ширины на части, а не только вторую стадию процесса, соответствующую установившемуся процессу продвижения вершины трещины. Предпочтение в представлении скорости изменения длины трещины при скачкообразном ее росте в материалах со структурой желательно отдать закону в виде (44). Было бы желательно отыскивать эти законы, опираясь на использование более или менее правдоподобных моделей деформирования материалов и накопления повреждений в них около вершин трещин при циклическом нагружении, причем, по мнению Ю. Мураками и К. Миллера, «усталостные повреждения должны быть выражены в терминах трещин»: «... Модели, которые игнорируют реальность усталостных повреждений в терминах трещин, не должны быть использованы для предсказания времени жизни при усталости» [36, с. 991].
При усталостном разрушении металлов «принципиальным следует считать вопрос о нелинейности накопления повреждений ..., которые могут быть реализованы на разных масштабных уровнях» [39, с. 14]. При нелинейном деформировании материалов в зонах предразрушения имеют место процессы самоорганизации системы на мезоуровне [40].
Подходы механики сплошной среды, соответствующие макромасштабному уровню описания, позволяют установить эмпирические закономерности развития трещин, опираясь на богатый экспериментальный материал [27, 41]. Например, основная часть перевернутой S-образной кривой диаграммы усталостного разрушения хорошо может быть представлена эмпирическими соотношениями Пэриса [41]. «Однако все попытки ввести единообразное описание кинетического процесса (роста трещины) до настоящего времени не дали положительного результата» ([42, с. 21]). «... Особенно в случае малоцикловой усталости, линейное суммирование накопленных повреждений не отражает реального, нелинейного процесса накопления повреждений» ([42, с. 37]).
6. Заключение
Предлагаемая диаграмма квазихрупкого разрушения при усталости описывает скачкообразный докритичес-кий рост трещин и предсказывает время жизни образца при усталости, используя аппроксимации стандартных CT-s-диаграмм квазихрупкого материала и характеристики накопления повреждений в материале зоны предразрушения. Модифицированная модель Леоно-ва-Панасюка-Дагдейла позволяет выразить в терминах трещин усталостные повреждения материала, расположенного в зоне предразрушения, как при линейном, так и при нелинейном суммировании повреждений. Время жизни образца существенно зависит от программы на-гружения. При рассмотренных режимах нагружения полученные соотношения (законы), описывающие рост трещин, не сводятся к простым законам роста трещин при усталости. Эти законы, в общем случае, не являются автомодельными процессами.
Работа выполнена при финансовой поддержке РФФИ (грант № 10-08-00220).
7. Приложение
Воспользуемся представлением соотношения (38) в таком виде, который учитывает ширину образца так, как это было выполнено в (12), (14), (15). Соотношения (12), учитывающие поправку на ширину образца через коэффициенты интенсивности напряжений, — достаточно грубая аппроксимация реальных условий нагружения. С учетом (12) для квазихрупкого приближения, когда
S* S*
A << / , имеем:
AS /r П (8nV/C (sv Г/С
N
2, j
(Y S KI2, j-1 )2
у/с
TCCTV r
yfn
r4r
r4r
V+ Vs
V Y KI2, j-1 Y KI2, j -1 y
1 с
1 jn k sY4r
2 л/n YS к I2, j-1
1-
1
(45)
к сту 1 - 2V
у = 1, 2, ..., /-1,
или в более приемлемом виде в соответствии с представлением закона в форме (44), ср. исходное соотношение (37):
AS/r П (8пУ/с (av Г/с N2, j ""î
у/ с
(YSsS A
v S Y y
/j-1
2/j-1 -Jn
Y+s+ Y SsS
V œ y
1 с
1 -
2/j-1 4n
Y scts k 1 - 2/,,,
œ j-1/
(46)
j = 1, 2, ..., j*-1
В соотношениях (45), (46) коэффициенты Y0 = Y0 (l0/w), Y+ = Y+ (l+/w), YS = YS (lE/w) зависят как от ширины образца, так и от длины трещин. Так как рассматривается квазихрупкое приближение, когда Д+/1+ << 1 и ДЕ/1S << 1, то при практических расчетах принято Y0 » Y+ » YS Вторые квадратные скобки не являются второстепенными сомножителями в соотношениях (45), (46) при завершении процесса, когда 2l ^ w.
Литература
1. ШанявскийА.А. Масштабные уровни бимодального распределения
усталостной долговечности титановых сплавов // Тезисы докладов Межд. конф. по физической мезомеханике, компьютерному конструированию и разработке новых материалов, 5-9 сентября 2011 г., Томск, Россия. - Томск: ИФПМ СО РАН, 2011. - С. 88-89.
2. КорневВ.М. Диаграммы квазихрупкого разрушения тел с иерархией
структур при малоцикловом нагружении // Физ. мезомех. - 2011. -Т. 14. - № 5. - С. 31-45.
3. ЛеоновМ.Я., ПанасюкВ.В. Развитие мельчайших трещин в твердом
теле // Прикл. механика. - 1959. - Т. 5. - № 4. - С. 391-401.
4. Dugdale D.S. Yielding of steel sheets containing slits // J. Mech. Phys.
Solids. - 1960. - V. 8. - No. 2. - P. 100-108.
5. Neuber G. Kerbspannunglehre: Grunglagen fur Genaue Spannungsrechnung. - Berlin: Springer-Verlag, 1937. - 160 s.
6. Новожилов В.В. О необходимом и достаточном критерии хрупкой
прочности // ПММ. - 1969. - Т. 33. - № 2. - С. 212-222.
7. Корнев В.М., ДемешкинА.Г. Диаграмма квазихрупкого разрушения
тел со структурой при наличии краевых трещин // ПМТФ. - 2011. -Т. 52. - № 6. - С. 152-164.
8. Корнев В.М. Оценочная диаграмма квазихрупкого разрушения тел с иерархией структур. Многомасштабные необходимые и достаточные критерии разрушения // Физ. мезомех. - 2010. -Т. 13. - № 1. - С. 47-59.
9. Керштейн И.М., КлюшниковВ.Д., Ломакин Е.В., Шестериков С.А.
Основы экспериментальной механики разрушения. - М.: Изд-во МГУ, 1989. - 140 с.
10. Leguillon D. Strength or toughness? A criterion for crack onset at a notch // Eur. J. Mech. A. Solids. - 2002. - V. 21. - P. 61-72.
11. Newman J.C., James M.A., Zerbst U. A review of the CTOA/CTOD fracture criterion // Eng. Fract. Mech. - 2003. - V. 70. - P. 371-385.
12. Castrodeza E.M., Perez Ipina J.E., Bastian F.L. Fracture toughness evaluation of unidirectional fiber metal laminates using traditional CTOD (5) and Schwalbe (S55) methodologies // Eng. Fract. Mech. -2004. - V. 71. - P. 1127-1138.
13. СаврукМ.П. Механика разрушения и прочность материалов. Т. 2. Коэффициенты интенсивности напряжений в телах с трещинами. -Киев: Наукова думка, 1988. - 618 с.
14. Справочник по коэффициентам интенсивности напряжений. В 2-х т. / Под ред. Ю. Мураками. Т. 1. - М.: Мир, 1990. - 448 с.
15. Партон В.З., Морозов Е.М. Механика упругопластического разрушения: Основы механики разрушения. - М.: Изд-во ЛКИ, 2008. - 352 с.
16. Thomson R., Hsieh C., Rana V. Lattice trapping of fracture cracks // J. Appl. Phys. - 1971. - V. 42. - No. 8. - P. 3154-3160.
17. Корнев В.М. Зона предразрушения для трещин продольного сдвига в материалах со структурой // Изв. РАН. МТТ. - 2011. - № 3. -С. 102-111.
18. Kitagawa H., Takahashi S. Applicability of Fracture Mechanics to Very Small Cracks or the Cracks in the Early Stages // Proc. II Int. Conf. Mechanical Behavior of Materials. - Metals Park, OH: American Society for Metals, 1976. - P. 627-631.
19. Kruzic J.J., Ritchie R.O. Kitagawa-Takahashi diagrams define the limiting conditions for cyclic failure in human dentin // J. Biomed. Mater. Res. A. - 2006. - P. 747-751.
20. ElHaddadM.H., Toper T.H., Smith K.N. Prediction of non-propagating cracks // Eng. Fract. Mech. - 1979. - V. 11. - P. 573-584.
21. Stouffer D.C., Williams J.F. A model for fatigue crack growth with a variable stress intensity factor // Eng. Fract. Mech. - 1979. - V. 11. -P. 525-536.
22. Корнев В.М. Двухмасштабная модель малоцикловой усталости. Переход от квазивязкого разрушения к хрупкому // Деформация и разрушение материалов. - 2008. - № 2. - С. 2-11.
23. Kornev V.M. Two-scale model of low-cycle fatigue. Embrittlement of pre-fracture zone material // Procedia Eng. - 2010. - V. 2. - No. 1. -P. 453-463.
24. ДемешкинА.Г., КарповЕ.В., КорневВ.М. Малоцикловая усталость образцов с краевой трещиной из сталей с разными степенями предварительного деформирования // Физ. мезомех. - 2009. -Т. 2. - № 3. - С. 91-99.
25. Kornev V., Karpov E., Demeshkin A. Damage accumulation in the pre-fracture zone under low-cyclic loading of specimens with the edge crack // Procedia Eng. - 2010. - V. 2. - No. 1. - P. 465-474.
26. Демешкин А.Г., Карпов Е.В., Корнев В.М. Накопление повреждений в образцах с краевой трещиной в зоне предразрушения при нестационарном малоцикловом нагружении // Изв. РАН. МТТ. -2011. - № 4. - С. 141-154.
27. Coffin L.F., Schenectady N.Y. A Study of the effects of cyclic thermal stresses on a ductile metal // Transactions ASME. - 1954. - V. 76. -No. 6. - P. 931-950.
28. Романив О.Н., Ярема С.Я., Никифорчин Г.Н. и др. Механика разрушения и прочность материалов. Т. 4. Усталость и циклическая трещиностойкость конструкционных материалов. - Киев: Наукова думка, 1990. - 680 с.
29. Никитенко А.Ф. Ползучесть и длительная прочность металлических материалов. - Новосибирск: Изд-во ИГИЛ СО РАН и НГАСУ, 1997. - 278 с.
30. Kornev V., Karpov E., Demeshkin A. Interrelation between Failure and Damage Accumulation in the Pre-fracture Zone under Low-cycling Loading // Recent Trends in Casting, Welding and Degradation of
Aluminium Alloys / Ed. by Zaki Ahmad. - Croatia: Tech. Rijeeka, 2011. - P. 407-422.
31. Paris P.C., GomezM.P., Anderson W.E. A rational analytic theory of fatigue // The Trend in Engineering. - 1961. - V. 13. - P. 9-14.
32. Броек Д. Основы механики разрушения. - М.: Высш. школа, 1980. - 368 с.
33. Paris P.C., Tada H., Donald J.K. Service load fatigue damage — a historical perspective // Int. J. Fatigue. - 1999. - V. 21. - P. S35-S46.
34. Ботвина Л.Р. Разрушение: кинетика, механизмы, общие закономерности. - М.: Наука, 2008. - 334 с.
35. CiavarellaM., PaggiM., CarpinteriA. One, no one, and one hundred thousand crack propagation laws: A generalized Barenblatt and Botvina dimensional analysis approach to fatigue crack growth // J. Mech. Phys. Solids. - 2008. - V. 56. - P. 3416-3432.
36. Murakami Y., Miller K.J. What is fatigue damage? A view point from the observation of low fatigue process // Int. J. Fatigue. - 2005. -V. 27. - P. 991-1005.
37. Nisitani H., GotoM., Kawagoshi N. A small-crack growth law and its related phenomena // Eng. Fract. Mech. - 1992. - V. 41. - P. 499-513.
38. Todinov M.T. Necessary and sufficient condition for additivity in the sense of the Palmgren-Miner rule // Comput. Mater. Sci. - 2001. -V. 21. - P. 101-110.
39. Шанявский А.А. Моделирование усталостных разрушений металлов. Синергетика в авиации. - Уфа: Монография, 2007. - 435 с.
40. ПанинВ.Е., Гриняев Ю.В., Данилов В.И. и др. Структурные уровни пластической деформации и разрушения / Под ред. В.Е. Панина. -Новосибирск: Наука, 1990. - 255 с.
41. Махутов Н.А., Гаденин М.М., Бурак М.И. и др. Механика малоциклового разрушения. - М.: Наука, 1986. - 264 с.
42. Шанявский А.А. Безопасное усталостное разрушение элементов авиаконструкций: Синергетика в инженерных приложениях. -Уфа: Монография, 2003. - 800 с.
Поступила в редакцию 26.03.2012 г.
Сведения об авторе
Корнев Владимир Михайлович, д.ф.-м.н., проф., гнс ИГИЛ СО РАН, [email protected]
