DOI: 10.18721/JPM.10208 УДК 539.3
симметризация тензора эффективных напряжений для сред с анизотропной поврежденностью
А.С. Семенов
Санкт-Петербургский политехнический университет Петра Великого Санкт-Петербург, Российская Федерация
Предложена единообразная форма представления процедуры симметризации тензора эффективных напряжений, содержащая классические аддитивную и мультипликативную схемы как частные случаи. Проведено сравнение различных вариантов симметризации тензора эффективных напряжений для однонаправленного поврежденного материала с параллельными микротрещинами и для двунаправленного поврежденного материала с системой ортогональных микротрещин. Предложена процедура идентификации анизотропной поврежденности на основе методов акустической эмиссии.
Ключевые слова: континуальная механика повреждений; тензор поврежденности; тензор эффективных напряжений; анизотропия; моделирование; идентификация
Ссылка при цитировании: Семенов А.С. Симметризация тензора эффективных напряжений для сред с анизотропной поврежденностью // Научно-технические ведомости СПБГПУ. Физико-математические науки. Т. 10. № 2. С. 82-98. DOI: 10.18721/JPM.10208
SYMMETRIZATION OF THE EFFECTIVE STRESS TENSOR FOR ANISOTROPIC DAMAGED CONTINUA
A.S. Semenov
Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russian Federation
A uniform formulation of the effective stress tensor symmetrization procedure has been proposed. This form contains the classical additive and multiplicative symmetrization schemes as particular cases. Different options of symmetrization of the effective stress tensor for the unidirectionally damaged material with parallel microcracks and for the bidirectionally damaged material with a system of orthogonal microcracks were compared. The differences in any forms of the effective stress tensor were second-order infinitesimals for low damage levels. The differences in predictions from considered symmetrization schemes increased with the growth of the damage level and with that of the differences between the eigenvalues of damage. An identification procedure for anisotropic damage was proposed on the basis of acoustic emission methods and then it was discussed.
Key words: continuum damage mechanics; damage tensor; effective stress tensor; anisotropy; simulation; identification
Citation: A.S. Semenov, Symmetrization of the effective stress tensor for anisotropic damaged continua, St. Petersburg Polytechnical State University Journal. Physics and Mathematics. 10(2) (2017) 82-98. DOI: 10.18721/JPM.10208
Введение
Тензор эффективных напряжений (ТЭН) является одним из центральных понятий механики континуальных повреждений
(МКП). Он был введен для одномерного случая в первых работах по механике поврежденного континуума Л.М. Качановым [1] и Ю.Н. Работновым [2]. Возникновение
и эволюция микродефектов (микротрещины и микропоры) приводит к сокращению эффективной площади несущих элементов квазиконтинуума, что в свою очередь приводит к эффекту возрастания внутренних напряжений (называемых эффективными), которые на макроуровне могут рассматриваться в качестве осреднения реально действующих на микроуровне напряжений в материале. Поэтому формулировки феноменологических определяющих уравнений, используемые в рамках МКП, как правило, связаны с заменой традиционного тензора напряжений Коши тензором эффективных напряжений.
Среди многочисленных обзоров отечественных и зарубежных исследований, выполненных с использованием МКП, следует упомянуть работы [3 — 7], отражающие современное состояние и обладающие полнотой охвата ключевых аспектов. Возможности МКП были продемонстрированы при анализе процессов ползучести [5 — 10], пластичности [6, 11, 12], усталости [13 — 17, 6], комбинации усталости и ползучести [6, 18], термоусталости [18] и хрупкого разрушения [19, 20, 6, 7].
В зависимости от характера пространственного и ориентационного распределения микродефектов могут использоваться скалярные [1, 2, 8, 21 и др.], векторные [22 — 24] и тензорные (второго [8, 25 — 28 и др.], четвертого [29 —31] и восьмого [32] рангов) параметры поврежденности. Возможны также различные подходы к определению тензора поврежденности:
на основе спектрального представления в собственном базисе, с заданием главных значений как относительных площадей дефектов в соответствующих направлениях [26];
на основе анализа морфологии микротрещин [25];
как к коэффициенту разложения в обобщенном ряду Фурье для ориентационного распределения плотности микродефектов
[5, 28];
путем введения эквивалентной нагруженной псевдонеповрежденной конфигурации и анализа соответствующего градиента деформаций [28] (геометрический подход).
Следует отметить, что использование тензора поврежденности второго ранга позволяет лишь описать поведение ортотроп-ных материалов, но не более сложных групп симметрии. Однако ортотропного приближения оказывается вполне достаточно для широко распространенного на практике случая, когда повреждения вызваны микрорастрескиванием, ориентационные свойства которого определяются направлением действия максимальных растягивающих напряжений.
В общем случае ТЭН, соответствующий тензору поврежденности второго ранга, оказывается несимметричным [26]. Корректное использование несимметричного тензора напряжений приводит к необходимости рассмотрения моделей сред Коссера и/или микрополярных сред. Для того чтобы избежать значительного усложнения модели, в ряде работ предлагается симметризировать полученный несимметричный тензор. В литературе предложены различные варианты симметризации ТЭН [26, 33, 34].
Целью данного исследования является сравнительный анализ различных вариантов симметризации ТЭН и последующий поиск единой однопараметрической формулировки, содержащей известные определения в качестве частных случаев.
Тензор эффективных напряжений
Изотропная поврежденность. В простейшем случае при одноосном растяжении силой Г прямолинейного стержня поперечной площади Б, часть которой БВ занимают микродефекты, эффективное напряжение ст (действует в неповрежденном материале вне пустот) определяется из условия
Г = ст Я = ст Б = ст (Б - БВ).
Если ввести поврежденность как долю площади поперечного сечения стержня, занятую дефектами В = БВ / Б, то эффективное напряжение ст = ¥ / Б выразится через традиционное понятие напряжения ст = Г /Б равенством
ст = ст / (1 - В).
В общем случае при многоосном напряженном состоянии изотропного по-
врежденного материала также используется скалярный параметр поврежденности Б (0 < Б < 1). Тогда ТЭН определяется следующим выражением:
1
о =
1 - Б
о.
(1)
Поврежденность при этом представляет собой поверхностную плотность микродефектов и определяется равенством
Б = 85Б / 55,
где дSD — площадь дефектов на произвольно ориентированной элементарной площадке 55*. При отсутствии повреждений (Б = 0) ТЭН о идентичен тензору напряжений Коши о.
При использовании гипотезы эквивалентности деформаций (г = г) тензор упругих модулей поврежденной среды 4 Се (для которого в упругом состоянии справедливо выражение о = 4 Се • • г) связан с тензором начальных модулей упругости (неповрежденного материала) 4 Се (для которого о = 4 Се • • г) равенством
4 Се = (1 - Б)4 Се. (2)
Альтернативное представление можно получить в случае использования гипотезы эквивалентности (упругих) энергий; оно имеет вид
4 Се = (1 - Б)24 Се. (3)
Анизотропная поврежденность. Элементарная ориентированная площадка п55 в актуальной конфигурации трансформируется в эквивалентной нагруженной псевдонеповрежденной конфигурации в п85. Аналогично рассмотренному выше одномерному случаю, симметричный тензор по-врежденности Б можно ввести как элемент оператора линейного преобразования п55 в п85 :
(1 - Б) • п85 = п85,
где 1 — единичный тензор.
Использование поврежденности в виде тензора второго ранга позволяет учесть анизотропию процесса накопления повреждений (возникновение и рост микротрещин перпендикулярно направлению
действия максимальных главных напряжений). Вследствие условия совпадения результирующих векторов «традиционных» и эффективных напряжений, приходим к соотношению
о • п8 5 = о пб 5, откуда можем определить ТЭН [9]:
о = о • (1 - Б)-1. (4)
В общем случае ТЭН, введенный определением (4), оказывается несимметричным. Симметричность наблюдается при изотропной поврежденности Б = Б1, а также при соосности тензора напряжений о и тензора поврежденности Б.
В общем случае, не ограниченном рассмотрением только тензоров поврежден-ности второго ранга, связь между эффективными и традиционными напряжениями осуществляется с помощью тензора влияния поврежденности четвертого ранга 4 М и представляет собой линейное преобразование общего вида [32]:
о = 4М •• о. (5)
Для рассмотренных случаев изотропной (1) и анизотропной (4) поврежденностей тензоры влияния определяются следующими выражениями:
1
4 М =
2(1 - Б)
(1<8>1 +1<1),
(6)
1
4 М = ^ [1®(1 - Б)-1 + 1<(1 - Б)"1], (7)
где использованы символы непрямого диа-дного умножения
(А®Б),И =4В,„ = 4^.
Появление двух слагаемых в правых частях выражений (4) и (5) обусловлено учетом симметрии тензора
о = о5 = (1/2)(о + о7),
где (...)5 и (...)7 — операции симметризации и транспонирования (перестановки индексов) тензора второго ранга.
Например, при получении выражения (4), как следствие формулы (1) имеем: (1 - Б) о = о = 1о1 = 1-(1/2)(о + о7)*! =
= е, 0 е.(1 / 2)(5Ла,и^ + 5Ла^5т-) = = е, ® е,(1/2)(5Л5иу +5т5к,-)а^ = = (1 / 2)( 101 + 101) • • о,
где 5 к — символ Кронекера, 0 — символ прямого диадного умножения
(А 0 в)J^fcí =4А •
При использовании гипотезы эквивалентности деформаций тензор упругих модулей поврежденной среды 4 Се связан с тензором начальных модулей упругости 4 Се равенством
4 Се = 4 М-1 ••4 Се • (8)
В случае использования гипотезы эквивалентности энергий получается равенство
4 Ce = 4 М^1 ••4 С ••4 M-т • (9)
Симметризация тензора эффективных напряжений
Явные схемы. Простейший вариант симметризации ТЭН, предложенный С. Мура-ками и Н. Оно [26], связан с вычислением симметричной части от тензора второго ранга, фигурирующего в правой части равенства (4):
о = 2[(1 - Б)-1 • о + о • (1 - D)-1]. (10)
Тензор влияния поврежденности может быть получен как прямое следствие формулы тензора (10) в виде
4М = 1 [10(1 - D)-1+ 10(1 - D)-1 +
+ (1 - Б)-101 + (1 - Б)-101].
Матрица, соответствующая тензору влияния поврежденности (11), в главных осях п1 , п2, п3 тензора поврежденности
D = Дп
имеет вид:
[M] =
п1 + а2п2
1 п2 + а3п3
М11 0 0 0 0 0
0 М22 0 0 0 0
0 0 М33 0 0 0
0 0 0 М44 0 0
0 0 0 0 М55 0
0 0 0 0 0 Мб6
,(12)
где диагональные элементы, связывающие нормальные и касательные компоненты напряжений, определяются выражениями различной структуры:
1
Мц =
1 Г
2 I1 -
1 - а
1
Л
(13)
1 - а
Элементы матрицы М22 ляются аналогично М,
2
и М
33
вычис-
Элементы М55 и на осно-индексов
М66 находятся аналогично М44, ■ ве циклической перестановки 1 ^ 2 ^ 3 ^ 1 ... .
Матрица (12) является диагональной, благодаря выбору системы координат, связанной с главными осями Б. В общем случае для произвольной системы координат она будет отлична от диагональной.
При получении матрицы [М] (12) была принята очередность индексов, соответствующая следующему порядку перечисления компонент тензора напряжений:
к
•"1Р ^22' ^33' Ч2> "23> "31)
Альтернативный вариант симметризации (4) был предложен в работе [33]:
о = (1 - Б)-1/2 -о (1 - Б)-1/2. (14)
В случае соосности тензора напряжений и тензора поврежденности оба определения (10) и (14) совпадают.
Тензор влияния, соответствующий процедуре мультипликативной симметризации (14), имеет вид:
1
М = - [(1 - Б)-1/20(1 - Б)
-1/2
(15)
+ (1 - Б)-1/2 0(1 - Б)-1/2]
1/2
Матрица, соответствующая тензору влияния (15), в главных осях тензора повреж-денности допускает представление также в виде диагональной матрицы (12), компоненты которой вычисляются на основе выражений:
1
М11 = ,
11 1 - а
(16)
М44 =
7(1 - АХ! - б2) ■
4
Так же, как и в предыдущем случае, элементы матрицы Ы22 и М33 вычисляются аналогично М11, а элементы М55 и М66 находятся аналогично М44, на основе циклической перестановки индексов 1 ^ 2 ^ 3 ^ 1 ... . Следует отметить, что выражения для М11 в (13) и (16) совпадают.
В дальнейшем, для удобства изложения материала, преобразование, связанное с равенством (10), будем называть аддитивной симметризацией, а с равенством (14) — мультипликативной симметризацией.
В данной работе предлагается использовать единое универсальное представление для симметризованного ТЭН, содержащее варианты (10) и (14) как частные случаи. Предлагаемое представление имеет следующий вид:
о = ^[(1 - Б)-а • о• (1 - Б)
(1 -а)
(17)
+ (1 - Б)-(1 -а) • о • (1 - Б)-а ],
где а — безразмерный скалярный параметр, изменяющийся в диапазоне от 0 до 1.
При а = 0 или 1 получаем процедуру аддитивной симметризации (10), при а = 1/2 — процедуру мультипликативной симметризации (14).
Тензор влияния, соответствующий обобщенной процедуре симметризации (17), имеет вид
1
4 М= [(1 - Б)-а ® (1 - Б)
,-(1 -а)
-(1-а)
+ (1-Б)-а®(1-Б; + (1 - Б)-(1-)) ®(1-Б)-а + + (1 - Б)--1-'0®(1-Б)-а].
(18)
Матрица для тензора влияния (18), в главных осях тензора поврежденности определяется матрицей диагонального вида (12), где
1
М44 = ^
М11 = , 11 1 - Б1
1
(1 - Б1)а(1 - Б2)1-а 1
(19)
(1 - Б1 )1-а(1 - Б2)а
М22 и М33 вычисляются аналогично Мп, так же как и М55 и М66 находятся аналогично М44, на основе циклической перестановки индексов 1 ^ 2 ^ 3 ^ 1 ... .
В качестве дальнейшего обобщения представления (17) можно рассмотреть такое двухпараметрическое представление:
о = ^[(1 - Б)-а •о• (1 - Б)- +
2 (20)
+ (1 - Б)- • о • (1 - Б)-а].
Оно содержит выражение (17) при Ь = 1 — а как частный случай. Процедура симметризации (20) позволяет учесть вариант ТЭН, предложенный Беттеном в статье [27], имеющий вид
о = (1 - Б)-1 • о • (1 - Б)-1,
а также независимость в выборе показателей в степенях напряжений о и тензора сплошности Т = 1 - Б в обобщенных формулировках определяющих уравнений ползучести [8].
Тензор влияния, соответствующий двух-параметрической процедуре симметризации (20), определяется следующим равенством:
4 М = 1[(1- Б)"в®(1 -Б)-Ь +
+ (1 - Б)-®(1 - Б)-Ь + (1 - Б)-Ь®(1 - Б)- + (21) + (1 - Б)-Ь®(1 - Б)-].
Неявные схемы. В работах [33, 34] рассматривается возможность неявного задания ТЭН в виде решения следующего уравнения:
о = |[(1 -Б) -о + о • (1 -Б)]. (22)
Условия для определения симметричного ТЭН, получаемого путем неявной аддитивной симметризации (22), обсуждаются в работе [34].
При введении предположения о симметричности ТЭН, которое выражается как
о = о5 = (1/2)(о + о7),
тензор влияния определяется в следующем виде:
4 М =
1(1®1 + 1®1) - (23)
-1 (10Б + 10Б + Б01 + Б01)
(23)
Матрица, соответствующая тензору влияния (23), в главных осях тензора поврежден-ности имеет вид (12), где
Ми =
1
М44 =
1 - а1 1
(24)
1 - а + А
2
Как и в предыдущих случаях, элементы матрицы М22 и М33 вычисляются аналогично М11, а элементы М55 и М66 находятся аналогично М44, на основе циклической перестановки индексов 1 ^ 2 ^ 3 ^ 1 ... .
Неявная мультипликативная симметризация определяется выражением
о = (1 - Б)-1/2 • о • (1 - Б)-1/2. (25)
Очевидно, что данное определение не отличается от варианта явной мультипликативной симметризации (14).
В качестве обобщения неявных схем симметризации (22) и (25) может рассматриваться однопараметрическое представление следующего вида:
о = ^[(1 -Б)" • о • (1 - Б)
1-я
+ (1 - Б)1-" ■ о • (1 - Б)"],
(26)
где а — безразмерный скалярный параметр, изменяющийся в диапазоне от 0 до 1 (как и для представления (17)).
При а = 0 (или а = 1) получаем процедуру аддитивной симметризации (22), при а = 1/2 — процедуру мультипликативной симметризации (25).
Тензор влияния, соответствующий обобщенной процедуре симметризации (26), имеет вид:
4М = 4[(1-Б)а0(1-Б)1-0 + + (1 - Б)а 0(1 - Б)1-а + (1 - Б)1-а 0(1 - Б)а + (27) + (1 - Б)1-а 0(1-Б)а ]-1.
Матрица влияния, записанная в главных осях тензора поврежденности, имеет вид (12), где
М11 =
1 - а1
М44 = 2[(1 - АГ(1 - А)1 + (1 - Б,)1-"(1 - Б2)а]-1
(28)
Как и для вышеприведенных случаев, М22 и М33 вычисляются аналогично М11, так же как и М55 и М66 находятся аналогично М 44, на основе циклической перестановки индексов 1 ^ 2 ^ 3 ^ 1 ... .
Возможность явного выполнения инверсии, присутствующей в тензоре (27), реализована при вычислении элементов матрицы (28), благодаря диагональному виду матрицы влияния (12) в главных осях тензора поврежденности.
Сравнительный анализ различных вариантов симметризации
Сравнение результатов, получаемых с использованием различных процедур явной (см. формулы (10), (14), (17)) и неявной (формулы (22), (25), (26)) симметризации, осуществлялось для двух идеализированных вариантов поврежденного материала с различной микроструктурой:
система одинаково ориентированных (параллельных) микротрещин (вертикальных или наклонных, см. рис. 1, а, Ь), периодически повторяющихся в двух направлениях;
система ортогональных микротрещин различной длины в каждом из направлений (см. рис. 1, с), периодическая в двух направлениях.
На макроуровне материал рассматривается как однородный. Информация о микроструктуре на макроуровне учитывается только путем введения тензора поврежден-ности. Тензорный характер поврежденно-сти обеспечивает возможность учета анизотропии свойств начально изотропного материала.
Система параллельных микротрещин. Тензор поврежденности в рассматриваемом случае (см. рис. 1, Ь) имеет вид
б = а1п1
п
(29)
где Б — максимальное собственное значение тензора поврежденности (0 < Б1 < 1),
1
Рис. 1. Идеализированные варианты материала с повреждениями:
представлены системы вертикальных (а), односторонне направленных, наклонных (Ь) и ортогональных (с)
микротрещин
определяемое соотношением между длинои микротрещины и размером неповрежденной перемычки (расстояние между ближайшими вершинами соседних коллинеарных трещин); п1 — нормаль к берегам микротрещин.
Особый интерес представляет простейший случай системы вертикальных микротрещин (рис. 1, а, где ф = 0, Б = Х)11®1), допускающий наглядную интерпретацию результатов. При одноосном растяжении макрообразца в направлении, перпендикулярном ориентации микротрещин, в нем на макроуровне будет реализовано одноосное однородное напряженное состояние; тензор напряжений выражается как
о = аи1 ® 1.
Этот тензор является соосным с тензором поврежденности. В данном случае прогнозы всех рассмотренных выше аддитивных и мультипликативных явных и неявных процедур симметризации ТЭН совпадают:
" (30)
1 - д
Этот очевидный результат можно было получить и из сопоставления матриц влияния [М] для рассматриваемых вариантов определения ТЭН. Все матрицы (13), (16), (19), (24) и (28) в главных осях тензора по-врежденности имеют одинаковые первые диагональные элементы М11.
Аналогичное утверждение справедливо также для второго и третьего диагональных
элементов. Отличия наблюдаются только в четвертом, пятом и шестом диагональных элементах матриц влияния. Это указывает на возможность появления различий при сдвиговом воздействии.
При плоском чистом сдвиге, тензор напряжений которого имеет вид
О = Т12<1 ® ] + ] ® 1),
наблюдается различие в прогнозах ТЭН.
Формулы для вычисления отношений т12 / т12, полученных как величины М44 при Д1 ^ 0 и Д2 = 0, для всех рассмотренных выше аддитивных и мультипликативных, явных и неявных процедур симметризации представлены в таблице. На основе полученных соотношений проведено сравнение различных схем симметризации, результаты которого показаны на рис. 2.
Наибольшие абсолютные значения отношения т12 / т12 наблюдаются для явной аддитивной схемы симметризации (10) (см. рис. 2, а) во всем диапазоне изменения поврежденности. Наименьшие значения Т2 / т12 соответствуют неявной аддитивной схеме симметризации (22). Кривые, соответствующие мультипликативной симметризации (явная (14) и неявная (25) схемы дают в этом случае совпадающий результат), занимает промежуточное положение
(рис. 2, а).
В отличие от явной, неявная схема аддитивной симметризации не демонстрирует сингулярности т12 / т12 при стремлении Д1 к единице.
Таблица
Отношения т 12 / т12 при плоском чистом сдвиге в системе вертикальных микротрещин
Процедура симметризации ТЭН Базовая формула Т2 / Т12 Номер кривой на рис. 2
Явная аддитивная (10) 1 (1 + - ) 2 ( ^ 41 1
мультипликативная (14) 1 >А - 4 2
обобщенная (17) 1Г 1 1 " 21_(1 - Д)а ' (1 - А)1-а 4
Неявная аддитивная (22) 1 1 - 4/2 3
мультипликативная (25) 1 л/1 - 4 -
обобщенная (26) 2 (1 - А)а + (1 - А)1-а 5
а)
# / # / 1 1
1 , / / / 1 / / / # / / /2
/ / / * ¿г ** ^^^^ »•
о.оо
0.25
0.50
0.75
1.00 А
Ь)
4
••-""Г
0.0
0.2
0.4
0.6
0.8
1.0 а
Рис. 2. Зависимость величины Т12 / т12 от поврежденности В1 (а) и параметра а (при значении 4 = 0,8) (Ь) в системе вертикальных микротрещин для различных вариантов симметризации тензора эффективных напряжений (см. таблицу)
Следует отметить, что сравнение вторых инвариантов ТЭН для рассмотренных схем симметризации со вторым инвариантом несимметричного ТЭН (4) указало на то, что наиболее близкий к (4) прогноз дается явной аддитивной схемой (10).
Результаты анализа влияния параметра а на отношение Т12 / т12 при фиксированном значении Л1 = 0,8 показаны на рис. 2, Ь.
Наибольшие различия значений х12 / т12 в явной (17) и неявной (26) схемах симметризации ТЭН наблюдаются при а = 0 и а = 1 (соответствуют случаю аддитивной симметризации) и составляют 45 %. Различия между аддитивной и мультипликативной процедурами явной симметризации составляет 25 %. С ростом уровня поврежденности Л1 указанные различия будут расти
(см. рис. 2, а). Важно отметить, что отношение т12 / т12 при явной симметризации всегда не меньше, чем при неявной.
В отличие от случая вертикальных микротрещин (см. рис. 1, а, ф = 0), при анализе материала с наклонными микротрещинами (рис. 1, Ь, ф ф 0), при одноосном растяжении представительного объема, ТЭН будет содержать не одну, а три ненулевых компоненты, и для каждой из них будут наблюдаться различия в прогнозах процедур симметризации. Результаты вычисления отношения СТП / стп при одноосном нагружении, когда
о = сти1 ® 1, Б = Д1п1 ® и1,
для рассмотренных явных (17) и неявных (26) методов симметризации определяются соответственно следующими уравнениями:
СТ,, 4 1
11 4
11 = соз ф
<СТ
+ 8Ш2 ф С082 ф
СТ
1 - Д 1
+ 81П ф +
1
(1 - А)а (1 - А)
1-а
11 4 —- = С08 ф-
СТ 1 - Д
1 • 4
+ 81П ф +
+ 81п2 ф С082 ф
(1 - Д)а + (1 - А)1
Предельный анализ соотношений (31) и (32) при Д ^ 0 показал, что различия для любых вариантов ТЭН имеют второй порядок малости О (Д2).
На рис. 3 представлено сравнение значений СТ11 / ст11, полученных при помощи процедур явной аддитивной, явной мультипликативной и неявной аддитивной симметризации при различных значениях угла ориентации микротрещин ф и фиксированном значении Д1 = 0,8. Максимальное значение СТ11 / ст11 (как и для рассмотренного выше случая т12 / т12 (см. рис. 2, а)) наблюдается для явной аддитивной схемы симметризации (10); наименьшее значение СТ11 / ст11 соответствует неявной аддитивной схеме симметризации (22); мультипликативная симметризация (явная и неявные схемы дают совпадающий результат) занимает промежуточное положение (см. рис. 3). Максимальные различия значений СТ11 / ст11 для явной и неявной схем симметризации наблюдаются при ф = 45°.
Характер изменения отношения стп / ст11 с ростом поврежденности Д1 для материала с наклонными микротрещинами в целом соответствует наблюдаемому для отношения (32) Т2 / т12 (сравните рис. 4, а и 2, а). Однако различия между прогнозами для различных схем симметризации ТЭН в рассматривае-
(31)
Рис. 3. Зависимости отношения СТ11 / ст11 от угла ф ориентации наклонных микротрещин для трех вариантов симметризации тензора эффективных напряжений:
явная аддитивная, представлена формулой (10) (кривая 1); явная мультипликативная, ф-ла (14) (2); неявная аддитивная, ф-ла (22) (3)
Рис. 4. Зависимости величины СТ11 / ст11 от поврежденности Д1 (а) и параметра а (при значении Д1 = 0,8) (Ь) в системе наклонных микротрещин (ф = 45°) для различных вариантов симметризации тензора эффективных напряжений:
явная аддитивная, представлена формулой (10) (кривая 1); явная мультипликативная, ф-ла (14) (2); неявная аддитивная, ф-ла (22) (3); явная, ф-ла (31) (4); неявная, ф-ла (32) (5)
мом случае (ф = 45°) менее значительные. Неявная аддитивная схема симметризации для СТ11 / ст11 в отличие от т12 / т12 демонстрируют сингулярное поведение СТ11 / ст11 при стремлении Д1 к единице.
Результаты анализа влияния параметра а на отношение СТ11 / ст11 при фиксированных значениях параметров Д1 = 0,8 и ф = 45° для системы наклонных трещин показаны на рис. 4, Ь. Так же, как и в случае вертикальных микротрещин для т12 / т12 (см. рис. 2, Ь), наибольшие различия в значениях СТ11 / ст11 для явной (31) и неявной (32) схем симметризации наблюдаются при а = 0 и а = 1 (соответствуют случаю аддитивной симметризации ТЭН). Отношение СТ11 / ст11 при явной симметризации ТЭН всегда не меньше, чем при неявной.
Система ортогональных микротрещин. Тензор поврежденности в рассматриваемом случае системы ортогональных микротрещин (см. рис. 1, с) имеет вид
Б = Дп ® п1 + Д2п2 ® п2. (33)
Величины Д и Д2 определяются соотношением между длинами трещин и промежутков между ними (расстояния между вершинами соседних коллинеарных трещин) в направлениях п1 и п2 соответственно. п1 и п2 —
это главные оси тензора поврежденности. Компоненты тензора поврежденности Д1 и Д2 удовлетворяют ограничениям
0 <Д <1 0 <Д2 <1.
Результаты вычисления отношения СТ11 / ст11 для материала, поврежденного системой ортогональных микротрещин (см. рис. 1, с) с тензором Б, определяемым по формуле (33), при одноосном нагружении, когда тензор напряжений выражается как
О = СТц1 ® 1 ,
для рассмотренных ранее явных (17) и неявных (26) процедур симметризации определяются соответственно следующими уравнениями:
СТ
СТ
11 С0й4 ф й1п4 ф . 2 2 11 = -- +-- + 81и ф С082 ф X
11 1 - Д 1 - Д у у
1 - Д 1
1
(1 - Д1)а(1 - Д2)1-а (1 - Д2)а(1 - Д1)
СТ11 С084 ф 81п4 ф . 2 2
11 - + --г11 + йт2 ф С082 ф :
1-а
(34)
СТ
11
1 - Д 1 - Д
2
(35)
^^ - Д)^^^ ^ Д^1- + (1 - Д2)а(1 - Д)1-а Различие уравнений (34) и (35) состоит
X
Рис. 5. Зависимость отношения СТ11 / ст11 от параметров поврежденности Д1 и Д2 в системе
ортогональных микротрещин при значении Ф = 45° для явной аддитивной симметризации (см. ф-лу (10))
только в выражениях, стоящих в квадратных скобках.
Используя неравенство
1 1 4
- + ->-,
X у X + у
справедливое при х > 0, у > 0, несложно показать, что значение отношения СТ11 / ст11,
вычисленное на основе формулы (34), всегда больше или равно соответствующего значения, найденного по формуле (35). Ранее введенные равенства (31) и (32) входят как частные случаи в формулы (34) и (35) при Д2 = 0.
Сравнение значений аи / сти, найденных с использованием формул (34) и (35), показало, что максимум отношения СТ11 / ст11 наблюдается для явной аддитивной схемы симметризации (10) (частный случай (34) при а = 0) при всех значениях угла ф ; минимальное значение а11 / ст11 соответствует неявной аддитивной схеме симметризации (22) (частный случай (35) при а = 0); мультипликативная симметризация (явная и неявная схемы дают совпадающий результат) занимает промежуточное положение. Максимальные различия значений СТ11 / ст11 для явной и неявной схем симметризации ТЭН наблюдаются при ф = 45°.
Характер изменения а11 / сти с ростом поврежденностей Д1 и Д2 соответствует наблюдаемому для односторонне ориентированных микротрещин (см. рис. 4, а). Зависимость СТ11 / ст11(Д1, Д2) для явной аддитивной схемы симметризации ТЭН (10) показана на рис. 5.
Влияние Д2 на угловое распределение
11 / ст11 от угла ориентации микротрещин для различных
Рис. 6. Зависимости величины ст
соотношений Д1 и Д2 при фиксированном значении Д1 = 0,8. Значение Д2= 0,7 (кривая 1), 0,5 (2), 0,0 (3)
стп / сти при фиксированном значении Б1 показано на рис. 6 для явной аддитивной схемы симметризации ТЭН (10). При стремлении значения Б2 к Б1 кривая стп / стп стремится к окружности, соответствующей случаю изотропной поврежденности.
Изолинии стп / стп на плоскости Б1 — Б2 (рис. 7), построенные с использованием различных схем симметризации ТЭН, показывают, что различие в прогнозах увеличивается с ростом разницы между Б1 и Б2. При Б1 = Б2 прогнозы совпадают. Наибольшее расхождение наблюдается при односторонне ориентированных микротрещинах с Б1 = 0 или Б2 = 0 (см. рис. 7). Разброс по углам ориентации микротрещин сглаживает расхождение в прогнозах схем симметризации ТЭН .
Влияние скалярного параметра а на отношение стп / сти для материала, поврежденного системой ортогональных микротрещин, имеет характер, аналогичный ранее рассмотренному варианту односторонне ориентированных микротрещин (см. рис. 4, Ь). Значения отношения стп / а11 при
0.0 0.2 0.4 0.6 0.8 7>
Рис. 7. Изолинии значений а11 / а11 на плоскости Б1 — Б2 для различных схем симметризации ТЭН: явная аддитивная, ф-ла (10) (кривая 1); явная мультипликативная, ф-ла (14) (2); неявная аддитивная, ф-ла (22) (3)
явной симметризации ТЭН всегда больше или равны соответствующим значениям при его неявной симметризации при любых допустимых значениях Б1, Б2 и ф, что непосредственно следует из анализа вида соотношений (34) и (35).
Идентификация компонент тензора поврежденности
Известные методы оценки поврежден-ности основаны на измерении физических величин [35]: модулей упругости и плотности материала, его микротвердости, скорости распространения в материале ультразвуковых волн, деформаций ускоренной стадии ползучести, электрического сопротивления материала и др.
При анализе анизотропной поврежденности представляются перспективными методы акустической эмиссии, основанные на измерении анизотропии акустоупругих свойств твердого тела [36, 37].
Фазовые скорости V распространения плоской монохроматической волны
u = A ехр Iк ^ • г - VI)
в направлении N в упругой среде, характеризуемой тензором модулей упругости 4 С1, определяются на основе решения уравнения
с1е1;(N • 4С • N - pv22 = 0,
представляющего собой условие определения собственных чисел акустического тензора
N • 4 е • N.
Выше нами введено обозначение р для плотности материала. При использовании гипотезы эквивалентности деформаций (8) уравнение для определения скоростей распространения волн V в поврежденной среде с тензором модулей упругости
4 Се = 4 M-1 ••4 C
имеет вид
ёе^ • 4M-1 ••4С1 • N -pV2) = 0. (36)
Решение уравнения (36) в общем случае имеет громоздкий вид. Упрощения можно достичь, если рассматривать распростране-
ние волны в направлении одной из главных осей тензора поврежденности, например вдоль третьей оси N = п3. В этом случае для начально изотропного или для начально ортотропного материала (при условии совпадения осей упругой анизотропии с главными осями тензора поврежденности) акустический тензор имеет диагональный вид и решение уравнения (36) допускает простое компактное представление:
^3 = ТМСРЛ < V, = л/м^р-1, (37)
V = М-С-р-1,
где у3 — скорость распространения продольной волны в направлении п3; у1, у2 — скорости поперечных волн в направлениях п1 и п2 соответственно.
Для начально изотропного упругого материала, характеризуемого константами Ляме 1, ц и определяемого тензором упругих модулей вида
4 Се = И 0 1 + ц(101 + 101), (38)
при использовании явной схемы симметризации (17), с учетом формулы (19), как следствие соотношений (37) получаем следующие выражения для скоростей:
>3 = {(1 - Д3)(Х + 2ц)/ р}1/2, V, = {[(1 - Д3)-"(1 - Д2)"-1 + <+ (1 - Д3)а-1(1 - Д2)-"]-12ц / Р}1/2, (39) у; = {[(1 - д3) -а (1 - д1) а 1 +
+ (1 - Д3)а-1(1 - ДГ"]- 12ц / Р}1/2.
При использовании неявной схемы симметризации (26), с учетом формулы (28), как следствие соотношений (37) получаем альтернативные выражения для скоростей:
V = {(1 - Д)(Х + 2ц)/ Р}1/2, У2 = {[(1 - Д3)а(1 - Д)1-а + <+ (1 - Д)1-а (1 - Д2)а ]/2 -ц / р}1/2, (40) у = {[(1 - Д3)а(1 - Д)1-" + + (1 - Д3Га (1 - Д)а ]/2 -ц / Р}1/2.
Выражения для скоростей у1,У2,У3 (39)
или (40) могут рассматриваться как система трех нелинейных алгебраических уравнений относительно трех неизвестных главных значений тензора поврежденности Д1, Д2, Д3. В ряде случаев решение можно получить аналитически в замкнутом виде.
В случае явной аддитивной симметризации главные поврежденности будут определяться как решения системы (39) при а = 0:
Д = 1 - (V У3)2,
Д = 1 - [2(у / У2)-2 - (V У3)-2]-1, (41)
Д = 1 - [2(У1/ у)-2 - (уъ/ уг)-2]-1,
где у3 = ^(X + 2ц) / р — скорость распространения продольной волны ^ в неповрежденном материале, у = = / р — скорости распространения поперечных волн в этом материале.
Выражения (41) можно непосредственно использовать для идентификации анизотропной поврежденности на основе данных измерения скоростей ультразвуковых сигналов [36 — 38].
В случае явной мультипликативной симметризации, в качестве следствия системы уравнений (39), при а = 1/2 имеем:
Д = 1 - (V У3)2, <Д2 = 1 - (У2/ У2)4 / (V, / (42)
Д = 1 - (*;/ ^/ф/ У3)2.
В случае неявной аддитивной симметризации ТЭН, поврежденности будут определяться на основе поиска решения системы уравнений (40) при а = 0 :
Д = 1 - Ф/ У3)2,
• Д2 = 1 - [2(У2/ У2)2 - (У3/ У3^], (43) Д = 1 - [2(\/ у)2 - (уъ/ У3)2].
Аналогичные выражения для повреж-денностей можно получить в случае использования гипотезы эквивалентности энергий (9). Скорости распространения волн в предположении постоянства поврежденности определяются в этом случае как решения характеристического уравнения:
ёе^ - 4М-1 - -4Се - -4М-т - N - ру2) = 0.(44)
Идентификация главных значений тензора поврежденности при использовании уравнения (44) производится аналогично рассмотренному выше случаю анализа (36). Для явной аддитивной симметризации (10) при этом имеем:
Б3 = 1 - V, / Vз,
Б2 = 1 - [2(V2/ V!) - (Vз/ У3)]-1, (45)
Б = 1 - [2(Vl/ V!) - (Vз/ Vз)]-1;
в случае явной мультипликативной симметризации (14) получаем выражения:
Б3 = 1 - V; / Vз,
Б2 = 1 - (V2 / 72)2/(7з/ Vз), (46)
Б = 1 - (Vl/ Vl)2/(Vз/ Vз);
в случае неявной аддитивной симметризации (22) приходим к уравнениям:
Б3 = 1 - 7з / Vз,
<Б2 = 1 - [2(72 / V2) - (Гз / Vз)], (47)
Б = 1 - [2(71/ 71) - (7з/ 7з)].
Произвол в выборе схемы симметризации ТЭН и гипотезы эквивалентности может быть устранен на основе анализа экспериментальных данных для конкретного вида материалов и условий нагружения. Необходимо обеспечить выполнение неравенств
0 <Б1 <1, 0 < Б2 < 1, 0 < Б3 <1,
монотонность увеличения поврежденности с ростом нагрузки, возможность получения термодинамически согласованной и наиболее простой формы уравнений эволюции для тензора поврежденности.
Заключение
Предложенная единая обобщенная од-нопараметрическая форма представления симметризованного ТЭН позволяет ввести в рассмотрение новые варианты ТЭН, а также анализировать свойства поврежденного континуума в общем виде для различных ТЭН.
На основе решения ряда тестовых задач для параллельных и ортогональных микротрещин в настоящей работе показано, что явная аддитивная схема симметризации дает верхнюю границу тензора эффективных напряжений (т. е. является консервативной оценкой), в то время как неявная аддитивная — его нижнюю границу. Различия в прогнозах рассмотренных схем симметризации увеличиваются с ростом уровня поврежденности, а также с ростом различий между главными значениями тензора поврежденности. При низких уровнях по-врежденности различия в любых вариантах ТЭН являются малыми второго порядка.
Вариативность в определении симме-тризованного ТЭН приводит к различным эквивалентным определениям поврежден-ности. В ряде случаев вопрос о выборе наиболее подходящего варианта ТЭН может быть решен на основе экспериментальных исследований в процессе идентификации параметров уравнения эволюции повреждений. В работе также предложена процедура идентификации анизотропной поврежден-ности на основе широко известных акустических методов.
Работа выполнена при финансовой поддержке гранта Российского научного фонда (проект №15-19-00091).
список литературы
1. Качанов л.М. О времени разрушения в условиях ползучести // Известия АН СССР. Отд-ние техн. наук. 1958. № 8. С. 26-31.
2. Работнов Ю.Н. О механизме длительного разрушения // Вопросы прочности материалов и конструкций. М.: Изд-во АН СССР, 1959. С. 5-7.
3. Волегов П.С., Грибов Д.С., Трусов П.В.
Поврежденность и разрушение: классические континуальные теории // Физическая мезоме-ханика. 2015. Т. 18. № 4. С. 68-86.
4. Докощенко А.М. Применение кинетической теории при анализе длительного высокотемпературного разрушения металлов в условиях сложного напряженного состояния. Обзор // Прикладная механика и техническая физика. 2012. T. 53. № 4. С. 149-164.
5. Murakami S. Continuum damage mechanics. Heidelberg, London, New York: Springer, 2012. 423 p.
6. Lemaitre J., Desmorat R. Engineering damage mechanics: ductile, creep, fatigue and brittle failures. Berlin: Springer, 2005. 380 p.
7. Krajcinovic D. Damage mechanics. Amsterdam: Elsevier, 1996. 761 p.
8. Работнов Ю.Н. Ползучесть элементов конструкций. М.: Наука, 1966. 752 с.
9. Huit J. Creep in continua and structures. Topics in applied continuum mechanics. Ed. by J.L. Zeman and F. Ziegler. Vienna: Springer, 1974. Pp. 137-155.
10. Leckie F., Hayhurst D. Creep rupture of structures // Proc. R. Soc. Lond. А. 1974. Vol. A340. No. 1622. Pp. 323-347.
11. Dragon A. Plasticity and ductile fracture damage: study of void growth in metals // Eng. Fract. Mech. 1985. Vol. 21. No. 4. Pp. 875-885.
12. Lemaitre J. A continuous damage mechanics model for ductile fracture // J. Eng. Mater. Technol. 1985. Vol. 107. No. 1. Pp. 83-89.
13. Lemaitre J. Evaluation of dissipation and damage in metals under dynamic loading // Proc. Int. Conf. on the Mech. Behavior of Materials I. Japan. 1971. Pp. 540-549.
14. Chaboche J.L. Une loi différentielle d'en-dommagement de fatigue avec cumulation non linéaire // Revue Française de Mécanique. 1974. No. 50-51. Pp. 71-82.
15. Семенов А.С. Исследование процесса образования зигзагов при распространении усталостной трещины методами континуальной механики повреждений // Нелинейные проблемы механики и физики деформируемого твердого тела. Под ред. К.Ф. Черных. Сборник трудов научной школы академика В.В. Новожилова. СПб.: Изд-во Политехнического университета, 2000. С. 186-211.
16. Семенов А.С. Механика нелокальных континуальных повреждений // Научно-технические ведомости СПБГПУ. 2006. № 1 (48). С. 148-159.
17. Мельников Б.Е., Семенов А.С. Энергетическая модель накопления повреждений при малоцикловой усталости // XXI Петербургские чтения по проблемам прочности. К 100-летию со дня рождения Л.М. Качанова и Ю.Н. Работнова. Сб. материалов. СПб.: Соло, 2014. С. 59-61.
18. Семенов А.С., Гецов Л.Б. Критерии термоусталостного разрушения монокристаллических жаропрочных сплавов и методы определения их параметров // Проблемы прочности. 2014. № 1. С. 50-62.
19. Krajcinovic D., Fonseka G.U. The continuous damage theory of brittle materials. Parts I, II // ASME J. Appl. Mech. 1981. Vol. 48. No. 4. Pp. 809-824.
20. Бенин А.В., Семенов А.С., Семенов С.Г., Мельников Б.Е. Математическое моделирова-
ние процесса разрушения сцепления арматуры с бетоном. Ч. 2. Модели без учета несплошности соединения // Инженерно-строительный журнал. 2014. № 1 (45). С. 23-40.
21. Павлов П.А. Основы инженерных расчетов элементов машин на усталость и длительную прочность. Ленинград: Машиностроение, 1988. 251 с.
22. Качанов Л.М. Основы механики разрушения. М.: Наука, 1974. 312 с.
23. Петреня Ю.К. Физико-механические основы континуальной механики повреждаемости. СПб.: Политехника, 1997. 253 с.
24. Davison L., Stevens A. Thermomechanical constitution of spalling elastic bodies // J. Appl. Phys. 1973. Vol. 44. No. 2. Pp. 667-674.
25. Вакуленко A.A., Качанов М.Л. Континуальная модель среды с трещинами // Изв. АН СССР. Механика твердого тела. 1971. № 4. С. 159-166.
26. Murakami S., Ohno N. A continuum theory of creep and creep damage // Creep in structures. A.R.S. Ponter, D.R. Hayhurst (eds.). Berlin: Springer, 1981. Pp. 422-443.
27. Betten J. Damage tensors in continuum mechanics // Journal de Mécanique théorique et appliqué. 1983. Vol. 2. No. 1. Pp. 13-32.
28. Радаев Ю.Н. Тензорные меры повреж-денности и гармонический анализ тонкой структуры поврежденности // Вестник Самарского гос. ун-та, 1998. № 2 (8). С. 79-105.
29. Chaboche J. Continuum damage mechanics: Part I. General concepts // J. Appl. Mech. 1988. Vol. 55. No. 1. Pp. 59-64.
30. Krajcinovic D. Damage mechanics // Mech. Mater. 1989. Vol. 8. No. 2-3. Pp. 117-197.
31. Leckie F.A., Onat E.T. Tensorial nature of damage measuring internal variables // Physical Non-Linearities in Structural Analysis. Berlin: Springer, 1981. Pp. 140-155.
32. Chaboche J.-L. Le concept de contrainte effective applique a l'elasticite et ala viscoplasti-cite en presence d'un endommagement anisotrope // Mechanical Behavior of Anisotropic Solids. J.P. Boehler (ed.). Hague: Martinus Nijhoff, 1982. Pp. 737-760.
33. Cordebois J.P., Sidoroff F. Damage induced elastic anisotropy // Mechanical Behavior of Anisotropic Solids, J.P. Boehler (ed.). Hague: Martinus Nijhoff, 1982. Pp. 761-774.
34. Betten J. Net-stress analysis in creep mechanics // Ingenieur-Archiv. 1982. Vol. 52. No. 6. Pp. 405-419.
35. Lemaitre J., Dufailly J. Damage measurements // Eng. Fract. Mech. 1987. Vol. 28. No. 5-6. Pp. 643-661.
36. Никитина Н.Е. Акустоупругость. Опыт практического применения. Н. Новгород: Та-лам, 2005. 208 с.
37. Беляев А.К., Лобачев А.М., Модестов В.С., Пивков А.В., Полянский В.А., Семенов А.С., Третьяков Д.А., Штукин Л.В. Оценка величины пластических деформаций с использо-
Статья поступила в редакцию 18.03.2017, принята
ванием акустической анизотропии // Известия Российской АН. Механика твердого тела. 2016. № 5. C. 124-131.
38. Belyaev A.K., Polyanskiy V.A., Lobachev A.M., et al. Propagation of sound waves in stressed elasto-plastic material // Proc. Int. Conf. Days on Diffraction. 2016. № 7756813. Pp. 56-61.
к публикации 09.04.2017.
сведения об авторе
СЕМЕНОВ Артем Семенович — кандидат физико-математических наук, доцент Санкт-Петербургского политехнического университета Петра Великого, Санкт-Петербург, Российская Федерация.
195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29 [email protected]
references
[1] L.M. Kachanov, O vremeni razrusheniya v usloviyakh polzuchesti [On the creep-rapture life], Izv. AN SSSR. Ser. Techn. Sci. (8) (1958) 26-31.
[2] Yu.N. Rabotnov, O mekhanizme dlitelnogo razrusheniya [On the delayed fracture mechanics], Voprosy prochnosti materialov i konstruktsiy [On the strength of materials and constructions], Moscow, AN SSSR (1959) 5-7.
[3] P.S. Volegov, D.S. Gribov, P.V. Trusov, Povrezhdennost i razrusheniye: klassicheskiye kontinualnyye teorii [The damage and fracture: Classical continuum theories], Fizicheskaya mezomekhanika. 18 (4) (2015) 68-86.
[4] A.M. Lokoshchenko, Primeneniye kineticheskoy teorii pri analize dlitelnogo vysokotemperaturnogo razrusheniya metallov v usloviyakh slozhnogo napryazhennogo sostoyaniya. Obzor [An application of the kinetic theory in the analysis of the delayed high-temperature fracture of metal under complex stress condition (Review)], Priklad. Mekhanika i Tekhnich. Fizika. 53 (4) (2012) 149-164.
[5] S. Murakami, Continuum damage mechanics, Springer, Heidelberg, London, New York (2012).
[6] J. Lemaitre, R. Desmorat, Engineering damage mechanics: ductile, creep, fatigue and brittle failures, Springer, Berlin (2005).
[7] D. Krajcinovic, Damage mechanics, Elsevier, Amsterdam (1996).
[8] Yu.N. Rabotnov, Polzuchest elementov konstruktsiy [Creep of structural elements], Moscow, Nauka (1966).
[9] J. Hult, Creep in continua and structures, Topics in Applied Continuum Mechanics, Ed. by J.L. Zeman and F.Ziegler, Vienna, Springer (1974) 137-155.
[10] F. Leckie, D. Hayhurst, Creep rupture of structures, Proc. R. Soc. Lond. A340 (1622) (1974) 323-347.
[11] A. Dragon, Plasticity and ductile fracture damage: study of void growth in metals, Eng. Fract. Mech. 21 (4) (1985) 875-885.
[12] J. Lemaitre, A continuous damage mechanics model for ductile fracture, J. Eng. Mater. Technol. 107 (1) (1985) 83-89.
[13] J. Lemaitre, Evaluation of dissipation and damage in metals under dynamic loading, Proc. Int. Conf. on the Mech. Behavior of Materials I. Japan, (1971) 540-549.
[14] J.L. Chaboche, Une loi différentielle d'en-dommagement de fatigue avec cumulation non linéaire, Revue Française de Mécanique. (50-51) (1974) 71-82.
[15] A.S. Semenov, Issledovaniye protsessa obrazovaniya zigzagov pri rasprostranenii ustalost-noy treshchiny metodami kontinualnoy mekhaniki povrezhdeniy [A study of zigzag formation during the fatigue crack propagation by methods of the continuum damage mechanics], Nelineynyye pro-blemy mekhaniki i fiziki deformiruyemogo tverdogo tela, Ed. Chernykh K.F., Sbornik Trudov Nauchnoy Shkoly akademika V.V. Novozhilova, St. Petersburg (2000) 186-211.
[16] A.S. Semenov, Mekhanika nelokalnykh kontinualnykh povrezhdeniy [Nonlocal continuum damage mechanics], St. Petersburg Polytechnical State University Journal. No. 1 (48) (2006) 148-159.
[17] B.E. Melnikov, A.S. Semenov,
Energeticheskaya model nakopleniya povrezhdeniy pri malotsiklovoy ustalosti [An energetic model of damage accumulation during low-cycle fatigue], XXI Peterburgskiye chteniya po problemam prochnosti. K 100-letiyu so dnya rozhdeniya L.M. Kachanova i Yu.N. Rabotnova. Sb. Materialov, St. Petersburg, Solo (2014) 59-61.
[18] A.S. Semenov, L.B. Getsov, Kriterii termoustalostnogo razrusheniya monokristal-
licheskikh zharoprochnykh splavov i metody opredeleniya ikh parametrov [Thermal fatigue fracture criteria of single crystal heat-resistant alloys and methods for identification of their parameters], Strength of Materials. 46(1) (2014) 38-48.
[19] D. Krajcinovic, G.U. Fonseka, The continuous damage theory of brittle materials, Parts I and II, ASME, J. Appl. Mech. 48 (4) (1981) 809-824.
[20] A.V. Benin, A.S. Semenov, S.G. Semenov, B.E. Melnikov, Matematicheskoye modelirovaniye protsessa razrusheniya stsepleniya armatury s betonom. Chast 2. Modeli bez ucheta nesploshnosti soyedineniya [The simulation of bond fracture between reinforcing bars and concrete. Part 2. Models without taking the bond discontinuity into account], Magazine of Civil Engineering. 45(1) (2014) 23-40.
[21] P.A. Pavlov, Osnovy inzhenernykh raschetov elementov mashin na ustalost i dlitelnuyu prochnost [A basis of design studies of machine parts for fatigue and long-term strength], Leningrad, Mashinostroyeniye (1988).
[22] L.M. Kachanov, Osnovy mekhaniki razrusheniya [A basis of the failure mechanics], Moscow, Nauka (1974).
[23] Yu.K. Petrenya, Fiziko-mekhanicheskiye osnovy kontinualnoy mekhaniki povrezhdayemosti [Physical and mechanical foundations of continuum damage mechanics], St. Petersburg, Politekhnika (1997).
[24] L. Davison, A. Stevens, Thermomechanical constitution of spalling elastic bodies, J. Appl. Phys. 44 (2) (1973) 667-674.
[25] A.A. Vakulenko, M.L. Kachanov, Kontinualnaya model sredy s treshchinami [Continuous model of the medium with cracks], Izv. AN SSSR. Mekh. Tverd. Tela. (4) (1971) 159-166.
[26] S. Murakami, N. Ohno, A continuum theory of creep and creep damage, Creep in Structures, A.R.S. Ponter, D.R. Hayhurst (eds.), Berlin, Springer (1981) 422-443.
[27] J. Betten, Damage tensors in continuum mechanics, Journal de Mécanique Théorique et Ap-
Received 05.03.2017, accepted 22.03. 2017.
plique. 2 (1) (1983) 13-32.
[28] Yu.N. Radayev, Tenzornyye mery povrezhdennosti i garmonicheskiy analiz tonkoy struktury povrezhdennosti [Tensor damage measures and a harmonic analysis of the damage fine structure], Vest. Samarskogo Gos. Un-ta. No. 2 (8) (1998) 79-105.
[29] J. Chaboche, Continuum damage mechanics: Part I, General concepts, J. Appl. Mech. 55(1) (1988) 59-64.
[30] D. Krajcinovic, Damage mechanics, Mech. Mater. 8 (2-3) (1989) 117-197.
[31] F.A. Leckie, E.T. Onat, Tensorial nature of damage measuring internal variables, Physical Non-Linearities in Structural Analysis, Springer, Berlin (1981) 140-155.
[32] J.-L. Chaboche, Le concept de contrainte effective applique a l'elasticite et a la viscoplasticite en presence d'un endommagement anisotrope, Mechanical Behavior of Anisotropic Solids, Boehler J.P. (ed.). Martinus Nijhoff, Hague (1982) 737-760.
[33] J.P. Cordebois, F. Sidoroff, Damage induced elastic anisotropy, Mechanical Behavior of Anisotropic Solids, Boehler J.P. (ed.). Martinus Nijhoff, Hague (1982) 761-774.
[34] J. Betten, Net-stress analysis in creep mechanics, Ing.-Archiv. 52 (6) (1982) 405-419.
[35] J. Lemaitre, J. Dufailly, Damage measurements, Eng. Fract. Mech. 28 (5-6) (1987) 643-661.
[36] N.E. Nikitina, Akustouprugost. Opyt prakticheskogo primeneniya [Acoustoelastisity. Practical experience], Nizhniy Novgorod, TALAM (2005)
[37] A.K. Belyayev, A.M. Lobachev, V.S. Modestov, et al., Otsenka velichiny plasticheskikh deformatsiy s ispolzovaniyem akusticheskoy anizotropii [Estimating the plastic strain with the use of acoustic anisotropy], Mechanics of Solids. 51(5) (2016) 606-611.
[38] A.K. Belyaev, V.A. Polyanskiy, A.M. Lobachev, et al., Propagation of sound waves in stressed elasto-plastic material, Proc. Int. Conf. Days on Diffraction. No. 7756813 (2016) 56-61.
the author
SEMENOV Artem S.
Peter the Great St. Petersburg Polytechnic University
29 Politechnicheskaya St., St. Petersburg, 195251, Russian Federation
© Санкт-Петербургский политехнический университет Петра Великого, 2017