CC BY
42
УДК 620.193
ОЦЕНКА ОПАСНОСТИ КАПЛЕУДАРНОЙ ЭРОЗИИ МЕТОДАМИ ТЕОРИИ КОНТАКТНЫХ УПРУГИХ ДЕФОРМАЦИЙ
© 2016 М.А. Скотникова, В.В. Елисеев, А.А. Москалец Санкт-Петербургский политехнический университет Петра Великого
Статья поступила в редакцию 02.03.2016
Удары капель жидкости по поверхности тела, вызывающие эрозию, рассматриваются как упругое контактное взаимодействие. Средствами теории упругости и компьютерной математики определено напряженное состояние тела и найдены параметры нагрузки, при которой начинаются пластические деформации. Выведена формула для скорости капель. при которой начинается эрозия.
Ключевые слова: каплеударная эрозия, контактное взаимодействие, напряжения в полупространстве, условие текучести, компьютерная математика
Каплеударная эрозия (КЭ) - это поверхностное разрушение тел под действием струи газа с каплями жидкости. Большое практическое значение изучения этого сложного и вредного явления обусловлено, в частности, возникновением КЭ в лопаточных аппаратах [1-8]. Выполнено много серьезных исследований как экспериментального, так и теоретического характера [9-17]. Литература по теме представлена не только множеством статей, но и монографий [16, 18-20]. Для полного понимания сущности КЭ желательно разобраться в двух процессах: гидродинамике капли при ударе о поверхность тела и деформировании тела локальной ударной нагрузкой. Многие авторы считают причиной КЭ накопление пластических контактных деформаций. Такая точка зрения принимается и в данной работе.
Цель работы: математическое моделирование возникновения КЭ посредством решения контактной задачи теории упругости.
Тело считается упругим полупространством, а капля - шариком радиусом R, падающим с некоторой скоростью v (рис. 1). Предпринята попытка построения на плоскости параметров v, R границы области, в которой КЭ не возникает. Работа опирается на литературу по контактным задачам теории упругости [21, 22], а также компьютерную математику [23]. Именно последняя (Mathcad) позволила довести классическое решение контактной задачи по Герцу до полной картины контактных напряжений в теле.
Формулы Герца. В результате решения задачи линейной теории упругости о контактном взаимодействии полупространства и шара (рис. 1) выведены следующие известные формулы Герца [20, 21]:
I_Ш 5_ГI23К-Уг.
Я 4а3 ' 4 ) '
Q = Ро2 П a
ß, ^+k
(i)
Рис. 1. Расчетная схема контакта Для контактного давления и силы имеем фор-
мулы
Р(г)= Ро
1 -
, Q = 2лJ p(r) rdr
где ро - значение в центре.
Нелинейное соотношение между силой и сближением можно выразить через потенциальную энергию взаимодействия П :
4 1/ з/ О _ — Я125/2 _ П'(5),
^ зр
2- - >3 (3ß У3 Q
П = — R>22572=2R-731 3ß 3ß 5 5 I 4
(2)
Эти формулы можно использовать для расчета максимального сближения и силы при ударе. По закону сохранения энергии (m - масса шарика, ро - плотность)
1 2 4 з
П = — mv , m = — %Rpn
2 3 0
(3)
где Q - прижимающая сила, a - радиус контактного круга, § - сближение тел; E, Eo - модули Юнга материалов (нолик относится к шару), v, vo - коэффициенты Пуассона.
Скотникова Маргарита Александровна, доктор технических наук, профессор. E-mail: [email protected] Елисеев Владимир Васильевич, доктор физико-математических наук, профессор. E-mail: [email protected] Москалец Артем Анатольевич, аспирант. E-mail: [email protected]
Однако формул Герца недостаточно для определения напряженно-деформированного состояния (НДС) тела.
Контактные напряжения в полупространстве. Для бесконечного тела нельзя применять распространенные программы конечноэлементного анализа (А^УБ). Зато весьма эффективны классические аналитические подходы. Например, решение Папковича-Нейбера [21, 24], представляющее вектор перемещения и через гармонические функции:
о
и = ^ ув - + [(3 - 4у) в - 2д;В - аВ ]к (4)
Здесь к - орт оси г (рис. 1), У - оператор Гамильтона в перпендикулярной плоскости, В, Во - гармонические скаляры ( ДВ = ДВо=0 - с оператором Лапласа);
д -д|дz .
Определив соответствующие перемещениям (4) деформации и далее напряжения (по закону Гука), обращаются к граничным условиям. На плоскости г=0 отсутствуют касательные напряжения, а нормальное напряжение равно нулю вне контактного круга и давлению с минусом - в круге. В результате, во-первых, находится соотношение
(1 - 2у) В = дгВ ^ В = -С1 - 2у)| В dz
(5)
Г -
- (р/а)2
^Jr2—22rpcosQT~pFTZ
(7)
ОСz, %) - 1п(z2 + ^z2 + %) % - г2 - 2грео8 0 + р2
, М ^го;
(8)
8г =-д 2 (zB + В] ),еф =-1 дг (zB + В0).
г
8z = 2(1 - 2у)дВ - zд2в -д2В],
8^ = (1 - 2у)дВ - zдгдВ -дгдzBo
(9)
Введены компактные выражения частных производных. Подставив (9) в закон Гука, определим напряжения , , Gz , Тrz . Далее найдем компоненты девиатора напряжений и норму Мизеса т [24]:
^ = (2°г -% -огуз,...,= Тгг;
Т = <[-(^2+^ф2+)+т2
(10)
Во-вторых, для гармонической функции В получается задача Неймана с заданной нормальной производной на границе, решаемая с потенциалом простого слоя (контактного круга) [21]:
В(г) = -!- Г РЫ^
4пцГ |г- ^ (6)
Здесь ^ = Е/ 2(1 + V) - модуль сдвига, г - радиус -
вектор в полупространстве; г1 - радиус-вектор на площадке контакта - по ней производится интегрирование.
Рассматриваемая задача осесимметричная, поэтому используем цилиндрические координаты г, ф, г (рис. 1). При этом компоненты векторов и тензоров будут функциями лишь двух переменных г, г. Формулу (6) можно записать в виде
2п а
В(г, z) = | de| /(г, z, р, 0^р,
4пЦ ] ]
Выражения компонент Бф, Бг отличаются заменой индексов. Начало пластических деформаций определяется условием текучести
т(г, z) = к
(11)
где к - предел текучести материала. Для титана, например, к=500 МПа. При достаточно слабом контактном воздействии (малая скорость капли) везде в материале будет т<к. Эрозия начнется, предположительно, при выполнении условия (11) в некоторой точке. Данный расчет должен показать, где и когда это произойдет.
В силу симметрии достаточно рассмотреть норму Мизеса на оси г: т(0, г)=Т(г). Построенный средствами Mathcad график этой функции показан на рис. 2.
При интегрировании в круге по переменной г1 используются полярные координаты р, 0. Для радиус-вектора г остаются г, ф, г, причем в осесимметричном случае можно считать ф=0. Двойной интеграл (7), зависящий от параметров г, г, надежно вычисляется в МаЬксай [23]. Во избежание расходимости на оси г интегрируем по р не от нуля, а от некоторого малого значения (например, Х=10-2а). Дифференцировать по г, г оказывается можно под интегралом, используя символьные вычисления - понадобится далее при вычислении перемещений и деформаций.
Определяя вторую функцию Во согласно (5), получим для нее интеграл как в (7), но вместо f будет
/] (г, z, Р,0) = -(1 - 2v)pN/l - (р/а)2 О( ^ %),
Рис. 2. Зависимость нормы Мизеса от глубины
На графике видно, что максимум нормы Мизеса достигается не на поверхности, а очень близко от нее в глубине тела. График типичный, проведены десятки расчетов с различными параметрами. Радиус контактного круга для графика а=10-4 м. Другие параметры: £=10п Па; У=0,3; к=5 108 Па (титан). Считалось ро=1, поскольку решение пропорционально этому множителю. Максимум на рис. 2 равен Т1=5,06 104 МПа и достигается на глубине а 3,96 10-3 мм. При начале пластической деформации будет
рТ = к ^ р0 = 9,88 -10-
(12)
Для сходимости интеграла в расчетах принималось М=10 а. Интегрирование и дифференцирование в МаЬксай надежно работали и для Во. По найденным функциям В, Во найдем деформации:
Оказалось, что значение ро не зависит от радиуса контактного круга а. Этот впечатляющий результат, впрочем, можно было предвидеть по соображениям подобия.
Расчет для удара капли. Все формулы для этого расчета уже представлены выше. Сначала из (1) находим силу О при начале пластической деформации:
2
3
3Q
2пРо
= (QR ]13 ^ Q = ( 3ß Y ( 2npo
I 4 I
(13)
Далее из (2, 3) определим соответствующую скорость капли
V =
ß2 п2
-Ро/2 =
1.56 • ß2
k
T V T
(14)
турбинных лопаток (участок минимальной скорости износа). Для определения эрозионной стойкости металла большое значение имеет первый (начальный) период, продолжительность которого в значительной степени характеризует сопротивляемость материала микроударному разрушению.
Константа Т определена выше (подчеркнуто). Величина скорости не зависит от радиуса капли К (удивительный результат) и по степенному закону зависит от предела текучести к и параметра р. Но возникает серьезная проблема - свойства воды как упругого тела. Для любого изотропного тела имеем пару упругих модулей; например, модуль Юнга и коэффициент Пуассона. Или же модули изменения объема К=£/3(1-2^ и сдвига р. Модуль К для воды определяется известной скоростью распространения звука. А вот модуль сдвига обычно принимают равным нулю. Но тогда вся логика наших расчетов рушится - Однако найдены экспериментальные данные по упругим свойствам воды [9, 25], где сообщается: v=0,5; р=1,3 10-5. Используя эти значения вместе с результатом (12) в формуле (14), получим критическое значение скорости капли v=177 м/с.
Обсуждение. При экспериментальном определении эрозионной стойкости различных материалов установлено, что процессы эрозии (глубина эрозионного износа 1э) протекают во времени в соответствии с рис. 3. Различают три типичных этапа протекания процессов эрозии (глубины эрозионного износа 1э) во времени. На первом этапе Т1, так называемом инкубационном периоде, видимых повреждений поверхности нет, потерь массы материала зафиксировать не удаётся. В ряде работ показано, что собственно эрозии материала предшествует рост микронапряжений, увеличение плотности дислокаций в металле предельного значения 1012-1013 см-1, зарождение под поверхностных микротрещин. Второй этап Т2 характеризуется тем, что имеет место максимальная скорость эрозии и в течение этого отрезка времени она остаётся практически постоянной. Во время третьего этапа тз по скорость эрозии снова ослабевает. Поэтому для оценки эрозионной стойкости материалов обычно используют значение скорости эрозии на первом или втором этапе, которое зависит от скорости встречи лопатки с частицами воды, формы и размеров этих частиц, частоты соударений, эрозионной стойкости материала лопаток и рельефа эродирующей поверхности. Приближённо глубину эрозионного износа 1э (или износ металла Аш) в зависимости от времени т (или массы выпадающей влаги) можно на участках т2 и т3 представить прямой линией.
В литературе известны лишь качественные объяснения вида кривой износа, основанные на представлениях об изменении характера взаимодействия между каплей и поверхностью в результате изменений рельефа эродированной поверхности. Согласно этим представлениям, по истечении инкубационного периода на поверхности появляются и накапливаются небольшие разрушения (начальный этап износа), возрастают глубина впадин (этап высокой скорости износа), углубления впадин до нескольких диаметров капель и заострение перемычек между кавернами (замедление износа), формирование «сотовой» и иглообразной формы поверхности, типичной для эрозии
Рис. 3. Характерная кривая каплеударного разрушения металла и отвечающие её участкам формы эродированной поверхности
Выводы: основным результатом работы является формула (14) для скорости капли при начале эрозии. Ее величина определяется упругими свойствами и пределом текучести материала, но не зависит от размера капли. Содержащаяся в формуле константа Т определяется из решения задачи теории упругости средствами компьютерной математики. Хотя работа опирается на точные решения, ее представления являются модельными и нуждаются в экспериментальной проверке.
Работа выполнена в рамках государственного задания при финансовой поддержке Минобрнауки России. Коды проектов: № 933-2014, № 1972-2014.
СПИСОК ЛИТЕРАТУРЫ:
1. Аверкина, Н.В. Влажно-паровая эрозия дисков и валов паровых турбин / Н.В. Аверкина, И.В. Железняк, Ю.Я. Ка-чуринер и др. // Электрические станции. 2010. № 8. С. 2735.
2. Дергачев, К.В. Особенности разработки и программной реализации имитационной модели эрозионного изнашивания рабочих лопаток мощных влажно-паровых турбин / К.В. Дергачев, Д.А. Коростелев // Вестник Брянского государственного технического университета. 2008. № 4. С. 49-57.
3. Ланина, А.А. Исследование высокоскоростного каплеударного воздействия на поверхность лопаток паровых турбин // Инструмент и технологии. № 28-29. 2008. С. 8487.
4. Селянский, А.С. Ударная эрозия рабочих лопаток цилиндра низкого давления // Вкник НТУУ «КП1». Машинобу-дування. 2009. № 55. C. 37-42.
5. Шубенко, А.Л. Влияние эрозии на основные эксплуатационные характеристики рабочей лопатки последней ступени цилиндра низкого давления мощной паровой турбины / А.Л. Шубенко, А.Э. Ковальский, Ю.С. Воробьев и др. // Проблемы машиностроения. 2010. Т. 13. № 1. С. 3-11.
6. Byeong-Eun, Lee. Development of a water droplet erosion model for Large Steam Turbine Blades / Byeong-Eun Lee, Kap-Jong Riu, Se-Hyun Shin, Soon-Bum Kwon // KSME International Journal, 2003. Т. 17. № 1, С. 114-121.
7. Sandeep, Soni. Erosion behaviour of steam turbine blades of glass-epoxy / Sandeep Soni, J.P. Pandey // International Journal of Advanced Engineering Technology, 2011. С. 110-117.
2.5
5
8. Qulan Zhou. Water Drop Erosion on Turbine Blades: Numerical Framework and Applications / Qulan Zhou, Na Li, Xi Chen et al. // Materials Transactions, 2008. Т. 49. № 7. С. 16061615.
9. Апакашев, Р.А. Определение предела прочности и модуля сдвига воды при малых скоростях течения / Р.А. Апакашев, В.В. Павлов // Изестия РАН. Механика жидкости и газа. 1997. № 1. С. 3-7.
10. Варавка, В.Н. Прочность и механизмы разрушения высокопластичных материалов при воздействии дискретного водно-капельного потока / В.Н. Варавка, О.В. Кудряков // Вестник Донского государственного технического университета. 2011. Т. 11. № 8-2 (59). С. 1376-1384.
11. Варавка, В.Н. Прогнозирование эрозионного износа титановых сплавов в условиях длительного каплеудар-ного воздействия / В.Н. Варавка, О.В. Кудряков, А.Ф. Медников // Вестник Донского государственного технического университета, 2011. № 9. С. 1574-1585.
12. Варавка, В.Н. Закономерности и параметры каплеудар-ной эрозии титановых сплавов / В.Н. Варавка, О.В. Куд-ряков // Известия высших учебных заведений. СевероКавказский регион. Серия: Технические науки. 2011. № 6. С. 92-98.
13. Кудряков, О.В. Механизмы формирования эрозионного износа металлических материалов при высокоскоростных капельных соударениях / О.В. Кудряков, В.Н. Варавка // Материаловедение. 2012. № 5. С. 36-43.
14. Рыженков, В.А. О механизме зарождения и развитии разрушения при каплеударном воздействии двухфазного потока / В.А. Рыженков, Г.В. Качалин, А.Ф. Медников //
Естественные и технические науки. 2012. № 5 (61). С. 289-295.
15. Тарелин, A.A. Электрофизические аспекты каплеударно-го разрушения элементов проточной части паровых турбин / A.A. Тарелин, Н.В. Сурду, А.В. Нечаев // Вестник НТУ «ХПИ». 2012. № 7. С. 88-96.
16. Фаддеев, И.П. Эрозия влажнопаровых турбин. - Л.: Машиностроение, 1974. 208 с.
17. Skotnikova, M.A. Transformation in Two Phase Titanium Alloys under High speed Mechanical Loading / M.A. Skotnikova, N.A. Krylov, G.D. Motovilina et al. // Вопросы материаловедения. 2007. № 4 (52). С. 359-365.
18. Спринжер, Дж. Эрозия при воздействии капель жидкости. - М.: Машиностроение, 1981. 200 с.
19. Ташлицкий, Н.И. Эрозия рабочих лопаток последней ступени влажно-паровой турбины. Библиографические ссылки // Вестник машиностроения. 2011. № 1. С. 95-96.
20. Перельман, Р.Г. Эрозия элементов паровых турбин / Р.Г. Перельман, В.В. Пряхин. - М.: Энергоатомиздат, 1986. 184 с.
21. Лурье, А.И. Теория упругости. - М.: Наука, 1970. 940 с.
22. Джонсон, К. Механика контактного взаимодействия. -М.: Мир, 1989. 510 с.
23. Кирьянов, Д.В. Mathcad 14. - СПб: БХВ-Петербург, 2007. 704 с.
24. Елисеев, В.В. Механика деформируемого твердого тела. -СПб: Изд-во Политехн. ун-та, 2006. 231 с.
25. Коренченко, А.Е. Определение модуля сдвига воды в экспериментах с плавающим диском / А.Е. Коренченко, В.П. Бескачко // ПМТФ. 2008. Т. 49, № 1. С. 100-103.
RISK ASSESSMENT FOR WATER DROP EROSION WITH APPLICATION OF ELASTIC CONTACT DEFORMATION THEORY
© 2016 M.A. Skotnikova, V.V. Eliseyev, A.A. Moskalets
Peter the Great St. Petersburg Polytechnic University
The impact of a water drop is considered to be an elastic contact interaction. The stressed condition and characteristics of load that cause plastic deformation are defined by means of the theory of elasticity and computer mathematics. The formula for droplet's velocity that gives rise to the erosion is derived.
Key words: water droplet erosion, contact interaction, stress in the semi-infinite space, yield condition, computer mathematics
Margarita Skotnikova, Doctor of technical Sciences, Professor. E-mail: [email protected]
Vladimir Eliseev, Doctor of Physics and Mathematics, Professor.
E-mail: [email protected]
Artem Moskalets, Post-graduate Student. E-mail:
