CC BY
79
2013 Механика № 2
УДК 539.4
М.В. Банников1, А.Ю. Федорова1, А.И. Терехина2, O.A. Плехов1
Институт механики сплошных сред УрО РАН, Пермь, Россия 2Пермский национальный исследовательский политехнический университет, Пермь, Россия
ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ ФРАКТАЛЬНЫХ ЗАКОНОМЕРНОСТЕЙ РОСТА УСТАЛОСТНОЙ ТРЕЩИНЫ И ДИССИПАЦИИ ЭНЕРГИИ В ЕЕ ВЕРШИНЕ
Статья посвящена исследованию зарождения и распространения усталостных трещин в титановых сплавах методом инфракрасной термографии при циклической нагрузке и с помощью анализа поверхности разрушения. Проведены две серии экспериментов на плоских гладких и на образцах с предварительно выращенной усталостной трещиной. Исследована пространственная и временная эволюция температуры в вершине трещины, определена форма пластической деформации и интенсивность тепловыделения. На основании результатов сравнительного анализа экспериментальных данных и уравнений линейной механики разрушения показано, что пространственная форма и характер диссипации тепла в зоне вершины трещины не соответствуют простой теоретической модели. В работе экспериментально получен эффект охлаждения за счет упругой деформации материала, а также исследованы особенности распределения напряжений в вершине трещины. Высокоскоростная съемка (на частоте 1 кГц) позволила определить интенсивность и форму зоны рассеиваемой энергии из-за пластической деформации в вершине трещины, а также сравнить скорость диссипации энергии для различных уровней напряжения. Разрушенные образцы исследовались с помощью электронного микроскопа и интерферометра New View для определения скейлинговых закономерностей образования и роста трещины. Полученные результаты позволяют нам проверить существующие модели неупругой деформации в вершине трещины, а также улучшить методы контроля накопления микроповреждений в процессе усталостного деформирования.
Ключевые слова: усталость, диссипация энергии, термоупругий эффект, инфракрасная термография, фрактография.
M.V. Bannikov1, A.Y. Fedorova1, A.I. Terekhina2, O.A. Plekhov1
institute of Continuous Media Mechanics UrB RAS, Perm, Russian Federation 2Perm National Research Polytechnic University, Perm, Russian Federation
EXPERIMENTAL STUDY OF FRACTAL PROPERTIES OF FATIGUE CRACK GROWTH AND ENERGY DISSIPATION IN CRACK TIP
Paper is devoted to investigation of initiation and propagation of cracks in titanium alloys by infrared thermography study of heat generation under cyclic loading and fractography analyses of fracture surface. Two series of experiments on flat smooth specimens and flat specimens with preliminary grown fatigue crack are carried out. Spatial and time temperature evolution into crack tip is investigated, the shape and dissipative intensity at crack process zone are estimated. Based on a result of comparative analysis of the experimental data and the linear fracture mechanics equations it is shown that the spatial distribution and character of heat dissipation zone into the crack tip doesn't correspond to the conventional models. High-speed shooting (at a frequency of 1 kHz) allowed us to determine the intensity and shape of zone of energy dissipation caused by plastic deformation at the crack tip, as well as to compare the rate of energy dissipation for different stress levels. Fractured specimens were analyzed by interferometer microscope and SEM to verify the existing models of inelastic deformation at the crack tip and to improve methods of monitoring of damage accumulation during fatigue test.
Keywords: fatigue, energy dissipation, thermo elasticity effect, infrared thermography, fractography.
Введение
Диссипация тепла, вызванная эволюцией структуры материала при циклическом деформировании, является объектом интенсивных исследований на протяжении последних десятилетий. В настоящее время достоверно известно, что в условиях циклического деформирования процессы локализации деформации сопровождаются интенсивным выделением тепла, что делает возможным их раннее обнаружение методами инфракрасной термографии [1].
Благодаря своей универсальности метод инфракрасной термографии в последнее время активно применяется при проведении механических испытаний как с целью получения детальной информации о процессе зарождения и распространения усталостных трещин [2-3], так и для исследования закономерностей преобразования и накопления энергии в процессе деформирования [4-9]. Возможности метода инфракрасной термографии позволяют в режиме реального времени исследовать процессы изменения температуры, вызванные термоупругостью и локализацией деформации в вершине трещины, а также эффектами трения на берегах трещины в процессе ее распространения.
Данная работа посвящена исследованию термоупругого и термопластического эффектов в вершине трещины, распространяющейся под действием растягивающего циклического напряжения, приложенного по нормали к плоскости трещины. В работе экспериментально получены эффекты охлаждения, вызванные упругой деформацией материала в вершине трещины, и исследованы особенности распределения напряжений в вершине трещины. Высокоскоростная съемка позволила определить интенсивность и форму зоны диссипации энергии, вызванную локализацией пластической деформации в вершине трещины, а также сопоставить интенсивность диссипации энергии для различных уровней напряжения.
1. Материалы и условия эксперимента
В работе исследованы особенности процесса тепловыделения при циклическом нагружении образцов титана ОТ-4 с частотами от 1 до 10 Гц [7]. Квазистатические испытания исследуемого материала позволяют определить его механические характеристики: модуль Юнга 64 ГПа, а0,2=683 МПа, ав=790 МПа.
Химический состав сплава ОТ-4
Элемент Весовой, %
Al 3,68
95,1
Mn 1,22
Испытания проводились на растяжной 100 кН сервогидравличе-ской машине Bi-00-100. Для определения величины локальной деформации во время эксперимента использовался осевой экстензометр Bi-06-304 с погрешностью ±1,5 мкм, установленный на берегах трещины. Исследуемые образцы изготавливались из листа толщиной 3 мм, размеры рабочей зоны 200x55 мм. Для получения начальной трещины образец был предварительно ослаблен тремя перекрывающимися отверстиями в центре рабочей зоны диаметром 4 и 3 мм, геометрия представлена на рис. 1.
На начальном этапе эксперимента с помощью повышенной нагрузки ((215±238) МПа, частота нагружения 20 Гц) создавалась усталостная трещина длиной 5 мм. В процессе регистрации эволюции поля
температуры нагрузка уменьшалась с целью замедления скорости распространения трещины (менее 1 мм/мин) и детального исследования процессов генерации тепла в ее вершине.
Рис. 1. Геометрия исследуемых образцов для изучения: эффекта термоупругости (а); тепловых эффектов в вершине трещины (б)
Рис. 2. Внешний вид экспериментальной установки: 1 - инфракрасная камера БЫК БС 5000; 2 - захваты сервогидравлической нагружающей машины В1-00-100; 3 - исследуемый образец
Исследование поля температур осуществлялось инфракрасной камерой БЫЯ БС 5000 (частота съемки до 20 кГц, пространственное разрешение до 2• 10-4 м). Поверхность образцов полировалась в не-
сколько этапов абразивной бумагой (на завершающей стадии полировки размер абразивных частиц не превышал 3 мкм), перед экспериментом полированная поверхность покрывалась тонким слоем аморфного углерода.
2. Теоретическое описание изменения температуры металла
В вершине усталостной трещины происходит интенсивная диссипация энергии, вызванная локализацией пластической деформации. Характерный размер зоны диссипации энергии в рамках линейной механики разрушения определяется величиной коэффициента интенсивности напряжений. Величину коэффициента интенсивности напряжений с учетом геометрии образца можно оценить, используя выражение
где Ж - ширина образца; а - полудлина трещины.
Эволюция температуры Т при циклическом деформировании
в отсутствие структурных переходов и пластической деформации, может быть описана уравнением Кельвина
где Р - коэффициент температурного расширения; р - плотность; с -удельная теплоемкость; V - коэффициент Пуассона; ю - угловая частота нагружения; (.) - производная по времени.
Решение уравнения (2) имеет вид
разложение в ряд соотношения (3) с учетом малости соотношения, стоящего в показателе экспоненты, дает следующие соотношения для амплитуды первой и второй гармоники:
в процессе циклического деформирования
(1)
аА +Аа8Іп(юґ) в предположении однородности ее распределения,
Да С08ЮҐ,
(2)
(3)
Р(1-2У) *„ Т Р(1-2у)л2
=ТР-Да, А2 =-0
р с 4
Р с
Да2. (4)
у
Анализ соотношений (3), (4) позволяет сделать вывод о том, что в классическом случае (соотношения (3)) амплитуда второй гармоники имеет порядок 10- и практически не может наблюдаться экспериментально, амплитуды первой и второй гармоник не зависят от частоты и являются линейными функциями амплитуды напряжения и квадрата амплитуды напряжения соответственно.
В настоящее время высказываются предположения о том, что эффект термоупругости является существенно нелинейным [1]. Заметный вклад в зависимость температуры от времени вносит процесс изменения упругих свойств материала от температуры. В предположении зависимости упругих модулей материала от температуры (Хт, цт ) изменение температуры описывается уравнением г
и>%т = -
К1"2у) +х ^0 + - СТо(1,5^2 + 2Хц + ц2)
У рс Т (3Х + 2ц)2 рс т ц2(3Х + 2ц)2 рс J
хДаюсоБ ю? +
г Хт (1,5Х2 + 2Хц + ц2)'
(5)
Да2ю8іп2юґ.
ч2(3Х + 2ц)2 рс ц2(3Х + 2ц)2 рс )
Анализ поведения амплитуд первой и второй гармоник проводился численно в предположении линейной зависимости упругих свойств материала от температуры. Величины дифференциалов упругих свойств Хт , цт подбирались из условия совпадения спектров зависимости температуры от времени для одного из экспериментов (при амплитуде напряжений 160 МПа и частоте 1 Гц) и использовались для остальных условий экспериментов. При расчетах принимались следующие значения физических констант: Р = 1х10~5 - коэффициент температурного расширения; р = 4300 кг/м3; с = 500 Дж/(кг-К); у = 0,34, коэффициенты зависимости упругих свойств от температуры получились соответственно: Хт =-9,4х107 Па/К, цг =-4,6х107 Па/К.
Радиус зоны пластической деформации на поверхности пластины
гг=*4, (6)
Р о2.
где к - коэффициент, зависящий от вида напряженного состояния и принимаемой модели пластической деформации; ау - напряжение
течения.
Форму зоны пластической деформации в вершине трещины при квазистатическом растяжении можно описать соотношениями: с учетом критерия Мизеса
На рис. 3 показано изменение максимальной температуры в вершине трещины, напряжения и раскрытия трещины в процессе нагружения при амплитуде напряжения 212 МПа, среднем напряжении 212 МПа и частоте 5 Гц. Данные датчика перемещения, закрепленного на берегах трещины, позволяют утверждать, что раскрытие берегов трещины изменяется синфазно с приложенным напряжением.
Анализ температурных данных позволяет утверждать, что максимум приложенных напряжений и максимум интенсивности тепловыделения в вершине усталостной трещины не совпадают во времени. В начале эксперимента образец нагружен средним напряжением и находится в состоянии термодинамического равновесия. На каждом цикле нагружения наблюдается участок падения температуры, вызванный термоупругим эффектом, переходящий в участок роста температуры, вызванный локальным переходом через предел пропорциональности и образованием зоны пластического деформирования. При падении напряжения тепловыделение в вершине трещины продолжается. Геометрия зоны пластического деформирования представлена на рис. 4. При уменьшении напряжения в вершине трещины тепловыделение возрастает, и температура достигает максимума практически при нулевом напряжении. Затем в начале следующего цикла температура падает за счет термоупругого эффекта, и процесс продолжается.
(7)
с учетом критерия Треска-Сен-Венана
3. Результаты экспериментального исследования изменения температуры от напряжения
Время, с
Рис. 3. Изменение максимальной температуры в вершине трещины (1), напряжения (2) и раскрытия трещины (3) в процессе нагружения при амплитуде напряжения 212 МПа, среднем напряжении 212 МПа и частоте 5 Гц
90
0.0006
270
Рис. 4. Форма зоны интенсивной диссипации энергии в вершине усталостной трещины на первом цикле деформирования и зоны пластической деформации, построенных по критериям Мизеса и Треска-Сен-Венана
Метод инфракрасной термографии позволяет с высокой точностью визуализировать зону интенсивной диссипации энергии в вершине усталостной трещины (см. рис. 4). Распределение температуры в вершине трещины в процессе деформирования может отличаться от формы зоны пластических деформаций из-за процессов теплопроводности, поэтому для анализа геометрии области интенсивного тепловыделения, вызванного пластической деформацией, логично использовать только первый цикл деформирования. На рис. 4 показано сравне-
ние наблюдаемой формы зоны интенсивной диссипации энергии на первом цикле деформирования и классических решений (3), (4) для формы зоны пластической деформации в вершине трещины (в расчетах использовались следующие данные: половина длины трещины 4 мм, приложенное напряжение 300 МПа). Анализ результатов позволяет сделать вывод только о качественном соответствии формы зоны пластической деформации в вершине усталостной трещины моделям Ми-зеса и Треска-Сен-Венана.
4. Фрактальные особенности роста усталостной трещины
Рост усталостной трещины описывается законом Пэриса:
= С(ЛК )т
Ш
(9)
где &а/Ш - скорость роста трещины; АК - размах коэффициента интенсивности напряжений; С и т константы, зависящие от свойств материала.
В рамках этого соотношения, построены кинетические кривые роста трещины, которые представлены в логарифмических координатах на рис. 5.
Рис. 5. Кинетические кривые роста усталостных трещин в титановом сплаве ОТ-4
В то время, как коэффициент С зависит от механических свойств материала и условий нагружения, о коэффициенте т ряд ученых имеет свое мнение. Так, в работах [10-13] его связывают с фрактальным рос-
том трещины. В работе [13] коэффициент т представлен следующим образом:
где В - фрактальная размерность. В настоящей работе проверено данное соотношение, и для определения фрактальной размерности трещины В использовалась технология ЯБ-анализа. Эта методика позволяет определить значение коэффициента Херста, который связан с фрактальной размерностью по следующей формуле:
Технология ЯБ-анализа заключается в следующем [14]:
- исследуемый ряд длиной N (Х_к:к = 1,...,«) делится на К непе-ресекающихся блоков,
- далее вычисляются отношения Я(_',й)/£(_',й) для каждого блока, где Я - размах накопленных сумм для каждого блока; £ - корень из среднеквадратичного отклонения для каждого блока.
Здесь ti - первые точки блоков, удовлетворяющие неравенству (ti - 1) + + й < N где d - количество элементов в блоке.
Я , ¿) = шах{0,Ж ,1),...,Ж , ¿)}~ шт{0,Ж ,1),...,Ж ,й)}, (13)
Затем находится значение дисперсии Б2 (ti ,й) для каждого блока, т.е.
для X.,...,X. . Отсюда для каждого блока вычисляются значения Я/Б;
ti ti+d-l
- после этого изменяем количество элементов d в блоке, т.е. ¿.¡+1=т-^, где т>1; для получения достоверного результата начальное количество элементов в блоке d0 берется равным 10, увеличение й ведется до тех пор, пока количество блоков будет не менее 6. Для каждого й повторяется процедура нахождения отношения Я(.,й)/Б(.,й);
т = -^, 1< В < 2, 2 - В
(10)
В = 2 - Н.
(11)
(12)
где
- используя полученные результаты, можно строить график зависимости log(R(ti,d)/S(ti,d)) относительно log d (рис. 6);
- применяется метод наименьших квадратов для определения углового коэффициента прямой, проходящей максимально близко к полученной зависимости, угловой коэффициент для данной прямой будет определять показатель Херста.
Исследуемые профили были получены с помощью цифрового интерферометра-профилометра New View, который позволяет получать 3-мерные картины профилей поверхности с точностью по горизонтали 0,5 мкм и по вертикали 1 нм. Из данных 3-мерных профилей вырезались 1-мерные, которые впоследствии анализировались R/S-методом. В результате вычислений получено следующее значение фрактальной размерности: D = 1,01 ± 0,01, которое при подстановке в формулу (10) дает близкий к экспериментальному показатель т.
R/S-анализ
| Linear model Poly 1:Дх)=0,9963 х-0,405 |
1 1,2 1,4 1,6 1,8 2 2,2
1оё10(<1)
Рис. 6. Зависимость Я/Б от количества элементов в блоке Л в логарифмических координатах
Заключение
Получены соотношения, описывающие изменение температуры на поверхности образца в процессе одноосного циклического деформирования с учетом линейного и нелинейного термоупругого эффекта. Исследование процессов изменения температуры при циклическом деформировании в титановом сплаве ОТ-4 на гладких образцах показало, что обратимое изменение температуры при упругом деформировании может достигать одного градуса
Показано, что процесс тепловыделения носит существенно нелинейный характер. Использование метода инфракрасной термографии позволяет эффективно исследовать процессы, связанные как с локализацией пластической деформации в вершине трещины, так и с трением на ее берегах.
Экспериментально показано, что зона пластической деформации не совпадает с предсказаниями линейной механики разрушения, а максимум тепловыделения достигается на нисходящей ветви нагружения.
Рассчитана фрактальная размерность профиля трещины, которая может быть представлена в кинетическом уравнении для роста усталостной трещины и согласуется с экспериментальными данными для сплава титана ОТ-4
Работа выполнена при финансовой поддержке РФФИ (грант № 12-01-31145 мол_а).
Библиографический список
1. Вавилов В.П. Динамическая тепловая томография (обзор) // Заводская лаборатория. Диагностика материалов. - 2006. - Т. 72, № 3. -C.26-36.
2. Fatigue crack initiation and growth in a 35CrMo4 steel investigated by infrared thermography / O. Plekhov, T. Palin-Luc, O. Naimark, S. Uva-rov, N. Saintier // Fatigue and fracture of engineering materials and structures. - 2005. - Vol. 28, iss. 1. - P. 169-178
3. Theoretical analysis, infrared and structural investigation of energy dissipation in metals under quasi-static and cyclic loading / O. Plekhov, N. Saintier, T. Palin-Luc, S. Uvarov, O. Naimark // Material Science and Engineering A. - 2007. - Vol. 462, no. 1. - P. 367-370.
4. Плехов О.А. Экспериментальное исследование термодинамики пластического деформирования методом инфракрасной термографии // ЖТФ. - 2011. - Т. 81, выи. 2. - С. 144-146.
5. Oliferuk W., Maj M., Raniecki B. Experimental analysis of energy storage rate components during tensile deformation of polycrystals // Materials Science and Engineering A. - 2004. - Vol. 374. - P. 77-81.
6. Плехов О.А., Santier N., Наймарк О.Б. Экспериментальное исследование накопления и диссипации энергии при упруго пластическом переходе // ЖТФ. - 2007. - Т. 77. - Вып. 9 - C. 1236-1238.
7. Экспериментальное исследование генерации тепла в вершине усталостной трещины / А.И. Терехина, М.В. Банников, О.А. Плехов,
Э.В. Плехова // Письма в ЖТФ. - 2012. - Т. 38. - Вып. 16.
8. Экспериментальное исследование аномалий диссипации энергии в нанокрокристаллическом титане при циклическом нагружении / О. Плехов, О. Наймарк, И. Семенова, Р. Валиев [и др.] // ПЖТФ. -2008. - Т. 34, выи. 13. - C. 33-40.
9. Исследование особенностей диссипации и накопления энергии в нанокристаллическом титане при квазистатическом и динамическом нагружении / О.А. Плехов, В.В. Чудинов, В.А.Леонтьев, О.Б. Наймарк // Вычислительная механика сплошных сред. - 2008. - Т. 1, № 4. -С. 69-78.
10. Walker K. The effect of stress ratio during crack propagation and fatigue for 2024-T3 and 7075-T6 Aluminium // Effects of environments and complex load history on fatigue life. ASTM STP 462. - Philadelphia, PA, 1970. - P. 1-14.
11. Carpinteri A., Paggi M. Self-similarity and crack growth instability in the correlation between the Paris’ constants // Eng. fracture mech. -2007. - No. 74. - P. 1041-1053.
12. Wnuk M.P., Yavari A. On estimation stress intensity factors and modulus of cohesion for fractal cracks // Eng. fracture mech. - 2003. -No. 70. - P. 1659-1674.
13. Paggi M., Carpinteri A. Fractal and multifractal approaches for the analysis of crack-size dependent scaling laws in fatigue // Chaos, Solitons and Fractals. - 2009. - No. 40. - P. 1136-1145.
14. Rose O. Estimation of the Hurst Parameter of Long-Range Dependent Time Series. Report. No. 137. February 1996. - Institute of Computer Science, University of Wurzburg, Germany, 1996.
References
1. Vavilov V.P. Dinamicheskaya teplovaya tomographiya (obzor) [Dynamic thermal tomography (Review)]. Zavodskaya laboratoriya. Diag-nostika materialov. 2006, vol. 72, no. 3, pp. 26-36.
2. Plekhov O., Palin-Luc T., Naimark O., Uvarov S., Saintier N. Fatigue crack initiation and growth in a 35CrMo4 steel investigated by infrared thermography. Fatigue and fracture of engineering materials and structures, 2005, vol. 28, iss. 1, pp. 169-178
3. Plekhov O., Saintier N., Palin-Luc T., Uvarov S., Naimark O. Theoretical analysis, infrared and structural investigation of energy dissipation in metals under quasi-static and cyclic loading. Material Science and Engineering A, 2007, vol. 462, no. 1, pp. 367-370.
4. Plekhov O.A. Experimental study of thermodynamics of plastic deformation by infrared thermography. Technical Physics, February 2011, vol. 56, iss. 2, pp. 301-304.
5. Oliferuk W., Maj M., Raniecki B. Experimental analysis of energy storage rate components during tensile deformation of polycrystals. Materials Science and Engineering A, 2004, vol. 374, pp. 77-81.
6. Plekhov O.A., Santier N., Naimark O.B. Experimentalnoe issledo-vanie nakopleniya i dissipacii energii pri uprugo plasticheskom perehode [Experimental study of the accumulation and dissipation of energy in elastic-plastic transition]. Zhurnal Tekhnicheskoi Fiziki, 2007, vol. 77, iss. 9, pp. 1236-1238.
7. Terekhina A.,I., Bannikov M.V., Plekhov O.A., Plekhova E.V. Experimental study of heat generation at the vertex of a fatigue crack. Technical Physics Letters, August 2012, vol. 38, iss. 8, pp. 736-738
8. Plekhov O.A., Naimark O.B., Valiev R.Z., Semenova I.P., Saintier N., Palin-Luc T. Experimental investigations of anomalous energy absorption in nanocrystalline titanium under cyclic loading conditions. Technical Physics Letters, July 2008, vol. 34, iss. 7, pp. 557-560.
9. Plekhov O., Chudinov V., Leont’ev V., Naimark O., Experimental investigations of the laws of energy dissipation during dynamic deformation of nanocrystalline titanium. Technical Physics Letters, January 2009, vol. 35, iss. 1, pp 92-95
10. Walker K. The effect of stress ratio during crack propagation and fatigue for 2024-T3 and 7075-T6 Aluminium. In: Effects of environments and complex and complex load history on fatigue life. ASTM STP 462. Philadelphia, PA; (1970), pp. 1-14.
11. Carpinteri A., Paggi M. Self-similarity and crack growth instability in the correlation between the Paris’ constants. Eng Fracture Mech., 2007, no. 74, pp.1041-1053.
12. Wnuk M.P., Yavari A. On estimation stress intensity factors and modulus of cohesion for fractal cracks. Eng. fracture mech., 2003, 70, pp. 1659-1674.
13. Paggi M., Carpinteri A., Fractal and multifractal approaches for the analysis of crack-size dependent scaling laws in fatigue. Chaos, Solitons and Fractals, 2009, no. 40, pp. 1136-1145.
14. Rose O. Estimation of the Hurst Parameter of Long-Range Dependent Time Series. Report No. 137. February 1996. Institute of Computer Science, University of Wurzburg, Germany.
Об авторах
Банников Михаил Владимирович (Пермь, Россия) - аспирант лаборатории физических основ прочности Института механики сплошных сред УрО РАН (614013, г. Пермь, ул. Академика Королева, 1, e-mail: mbannikov@icmm.ru).
Федорова Анастасия Юрьевна (Пермь, Россия) - аспирант лаборатории физических основ прочности Института механики сплошных сред УрО РАН (614013, г. Пермь, ул. Академика Королева, 1, e-mail: fedorova@icmm.ru).
Терехина Алена Игоревна (Пермь, Россия) - студентка Пермского национального исследовательского политехнического университета (614990, г. Пермь, Комсомольский пр., 29, e-mail:
tai_mm_9 0 @mail. ru).
Плехов Олег Анатольевич (Пермь, Россия) - доктор физикоматематических наук Института механики сплошных сред УрО РАН (614013, г. Пермь, ул. Академика Королева, 1, e-mail: poa@icmm.ru).
About the authors
Bannikov Mikhail Vladimirovich (Perm, Russian Federation) -postgraduate student in laboratory of physical foundations of strength, Institute of Continuous Media Mechanics, Ural Branch of Russian Academy of Sciences (1, Akademik Korolev st., 614013, Perm, Russian Federation, e-mail: mbannikov@icmm.ru).
Fedorova Anastasia Yurievna (Perm, Russian Federation) - postgraduate student in laboratory of physical foundations of strength, Institute of Continuous Media Mechanics, Ural Branch of the Russian Academy of Sciences (1, Akademik Korolev st., 614013, Perm, Russian Federation, e-mail: fedorova@icmm.ru).
Terekhina Alyona Igorevna (Perm, Russian Federation) - student, Perm National Research Polytechnic University (29, Komsomolsky av., 614990, Perm, Russian Federation, e-mail: tai_mm_90@mail.ru).
Plekhov Oleg Anatol’evich (Perm, Russian Federation) - Doctor of Physical and Mathematical Sciences, Institute of Continuous Media Mechanics, Ural Branch of Russian Academy of Sciences (1, Akademik Korolev st., 614013, Perm, Russian Federation, e-mail: poa@icmm.ru).
Получено 28.02.13
