CC BY
90
____________УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА
Том 156, кн. 2 Физико-математические науки
2014
УДК 532.5.296:534.12
УДАРНОЕ ВОЗДЕЙСТВИЕ СТРУИ ЖИДКОСТИ НА УПРУГОПЛАСТИЧЕСКОЕ ТЕЛО
А.А. Аганин, Н.А. Хисматуллина
Аннотация
Исследована динамика приповерхностного слоя тела из никелевого сплава в случае, когда к его поверхности приложена нагрузка, характерная для импульсного воздействия микроструи жидкости, образующейся при схлопывании касающегося тела кавитационного пузырька в воде. Рассмотрен диапазон скоростей струи, при которых интенсивность напряжений в теле кратковременно находится на уровне предела текучести. Представлены используемая математическая модель динамики упругопластического тела и основные положения методики расчета. Показано, что в результате увеличения области нагружения на поверхности тела, неоднородности нагружения и скачкообразного уменьшения нагрузки на периферии области нагружения предельные значения интенсивности напряжений в теле возникают при значительно (примерно в два раза) меньших скоростях струи, чем в соответствующем одномерном приближении. Неоднородность нагрузки не вызывает качественных изменений в динамике тела, ее влияние с ростом скорости струи уменьшается.
Ключевые слова: кавитационное разрушение, кавитационная эрозия, ударное воздействие струи, упругопластическое тело.
Введение
С воздействием струй (и капель) жидкости на поверхность деформируемых тел связана в первую очередь такая проблема, как эрозионное разрушение деталей машин и механизмов, контактирующих с жидкостью. Основным фактором кавитационного разрушения многие авторы (см., например, [1] и цитируемую там литературу) считают удар струи, образующейся на поверхности пузырька при его асимметричном коллапсе вблизи тела. Удар высокоскоростной струи жидкости по поверхности тела может привести к ее деформации и разрушению. Характер и масштабы повреждений зависят от плотности жидкости, радиуса и скорости струи, прочностных характеристик тела и т. д. Например [2], капля воды радиуса 1 мм при скорости 750 м/с может вызвать разрушение и эрозию таких твердых материалов, как алмаз и карбид вольфрама, и пластическую деформацию высокопрочных сплавов, таких как мартенситные стали. При более умеренных скоростях, порядка 200 м/с, одиночный удар может не вызывать видимых изменений поверхности, но многократные удары приводят к эрозии материала. В то же время применение эрозионных технологий в промышленности обусловливает необходимость контролируемого использования разрушительного потенциала ударного воздействия на тело струи жидкости. Для понимания природы разрушения поверхностей тел струями важно изучить особенности однократного воздействия струи. Этому вопросу посвящено много работ (см., например, [2-12]), большинство из которых являются экспериментальными. Между тем возможности экспериментального исследования происходящих внутри тела изменений весьма ограничены. Наблюдению в общем случае, и это отмечают многие авторы, доступны результаты воздействия
72
УДАРНОЕ ВОЗДЕЙСТВИЕ СТРУИ ЖИДКОСТИ...
73
Рис. 1. Цилиндрическая струя жидкости с полусферическим концом в начале (при t = = 0) ударного воздействия на тело и используемая цилиндрическая система отсчета r, z. Прямоугольником выделена область, где влияние струи наиболее существенно
на тело, а не сам его процесс. Поэтому для изучения особенностей процесса воздействия большое значение приобретает численное моделирование. В большей части работ, посвященных численному исследованию ударного воздействия капель (см., например, [3, 4, 7, 8]), рассматривается в основном динамика жидкости в капле, определяется характер нагружения тела. Работ, моделирующих воздействие струи на динамику тела, довольно мало. Так, динамика тела численно исследуется при нагружениях, соответствующих стационарному растеканию струи [13], и пульсирующем нагружении [6] поверхности тела.
В настоящей работе изучается динамика приповерхностной части тела при нагрузке, характерной для ударного воздействии струи жидкости, возникающей на поверхности касающихся тела кавитационных пузырьков при их сжатии [14].
1. Постановка задачи
Рассматривается динамика приповерхностного слоя тела при столкновении с его поверхностью микроструи жидкости. Струя цилиндрическая, с полусферическим концом, ударяется о неподвижное тело ортогонально его плоской поверхности (рис. 1). Давление на поверхности тела вне области воздействия струи ро = = 0.1 МПа. Материал тела - никелевый сплав с плотностью р = 8000 кг/м3, модулем упругости E = 196000 МПа, коэффициентом Пуассона v = 0.3, пределом текучести при одноосном растяжении Уо = 125 МПа. Такие сплавы используются для покрытия поверхностей деталей машин и механизмов, контактирующих с жидкостью, для их защиты от разрушительного воздействия кавитации. Жидкостью является вода с плотностью р = 1000 кг/м3, скоростью звука со = = 1500 м/c. Радиус струи Rj = 20 мкм, ее скорость Vj варьируется в диапазоне 65 м/с < Vj < 650 м/c. Подобные струи возникают на поверхности касающихся тела кавитационных пузырьков в процессе их схлопывания [15].
В результате столкновения струи с поверхностью тела в струе возникает распространяющаяся от тела куполообразная ударная волна. Сначала край этой волны совпадает с границей круговой быстро расширяющейся области контакта струи и тела, но вскоре отрывается от нее, после чего все точки фронта ударной волны удаляются от поверхности тела. До отрыва края ударной волны от поверхности тела давление жидкости в струе между телом и ударной волной является весьма высоким (~ pocoVj), так что поверхность тела испытывает большие локальные нагрузки, рассматриваемые в качестве одной из наиболее вероятных причин повреждений. При отрыве края ударной волны от поверхности тела на боковой
74
А.А. АГАНИН, Н.А. ХИСМАТУЛЛИНА
поверхности струи между краем ударной волны и поверхностью тела возникает цилиндрическая волна разрежения, сходящаяся в жидкости к оси струи. В области этой волны образуется пристеночная цилиндрическая струйка жидкости, разлетающейся со скоростью, в несколько раз превышающей скорость натекающей струи. По мере радиального схождения волны разрежения в струе давление жидкости на поверхности тела уменьшается. Со временем на поверхности тела устанавливается профиль давления, характерный для стационарного натекания струи с максимумом в центральной точке (точке торможения жидкости на оси натекающей струи). В момент столкновения струи и тела давление в центре области контакта скачком возрастает от ро до p* « ро cqVj . На периферии области контакта оно принимает более высокие значения. В теоретических работах [16] отмечается, что эти значения имеют порядок ~ 2р* , хотя при численном моделировании в очень узкой области получаются и намного большие величины [17]. Максимум давления при стационарном натекании струи на тело, которое устанавливается после завершения переходных нестационарных процессов в центре области контакта, равен p*t « 0.5роv2. В рассматриваемом диапазоне скоростей струи давление 2р* значительно выше, чем p*t. Если момент столкновения струи с телом принять за начало отсчета времени t = 0, что и делается в настоящей работе, то максимум давления 2р* достигается при t « т = RjVj/2с0, а указанная выше волна разрежения в струе сходится в центре области контакта при t « 3т, после чего давление на поверхности тела имеет порядок p*t. В настоящей работе рассматривается начальный промежуток времени, в котором максимум давления на поверхности тела сначала скачком возрастает от ро до ~ р* , затем плавно увеличивается до ~ 2р* , после чего падает до ~ p*t. Пластические деформации в теле из рассматриваемого никелевого сплава появляются при повышении давления на поверхности тела до величин более рр « 2.25 кбар . В рассматриваемом диапазоне скоростей струи 65 м/c < Vj < 650 м/c такие нагрузки возникают при 2р* > рр, то есть при скорости струи Vj > 75 м/c. Рассматривается наиболее часто встречающийся случай, когда диаметр струи намного меньше линейных размеров тела. Поэтому тело моделируется как упругопластическое полупространство.
2. Математическая модель и основные положения методики расчета
Аппроксимация закона нагружения поверхности тела принимается в виде рь(трЬ) = ро, 0 < r < ж при t< 0,
рь{грЪ)
\ре(t) - ра(Ь)] (r/Rb)2 + рс(), 0 < Г < Rb(t) ро, Г > Rb(t)
при t > 0,
(1)
где r - радиальная координата, отсчитываемая от расположенного на оси симметрии струи центра круговой области нагружения радиуса Rb(t). Входящие в (1) функции рс(t), ре(t) представляют собой изменение нагрузки в центре r = 0 и на краю r = Rb(t) области нагружения соответственно. Они определяются выражениями
рс(t)
р* (bt2 + 1) , если р* (bt2 + 1) < р^уЬ), Ve(t), если р* (bt 2 + 1) > ре^),
ре (t)
^aiR^1 + lj р*, если 0 < t < т, aR^р*, если т < t < 3т,
р*
. Vjросорс1
ро +------:----,
росо + рс1
t
j rl
Rj
RL(t) = Rj vj
Rj , = 2c2 ,
(2)
УДАРНОЕ ВОЗДЕЙСТВИЕ СТРУИ ЖИДКОСТИ...
75
Рис. 2. Радиальные профили нагрузки (1) для шести моментов времени при Vj = 100 м/с (а) и Vj = 300 м/с (б)
Рис. 3. Зависимости безразмерного радиуса Rb/Rj (сплошная кривая) круговой области нагружения поверхности тела и безразмерной скорости ее расширения Vb/vj (штриховая кривая) от безразмерного времени. Точки соответствуют приведенным на рис. 2 профилям нагрузок
где ci - скорость распространения продольных возмущений в теле, a\, a,2, ai, а2, b - безразмерные константы. При указанных входных данных скорости Vj = = 100 м/c соответствует давление р* = 145 МПа, скорости Vj = 300 м/с-давление р* = 436 МПа, а скорости Vj = 500 м/с - давление р* = 726 МПа. На рис. 2 представлены рассчитанные по формуле (1) радиальные профили нагрузки для шести характерных моментов времени для скоростей струи Vj = 100 м/с (ai = = 50.07, ai = 1.38, a2 = 1.5 • 10-9, a2 = —7.4, b = -1.1 ■ 10-7) и Vj = 300 м/с
(a1 = 14.82, a1 = 1.38, a2 = 4.8 ■ 10-5, a2 = —5.5, b = —2.7 ■ 10-8).
Известно [3], что на начальном отрезке времени 0 < t < т радиус Rl круговой области нагружения поверхности тела определяется выражением
Rl = \J 2RVj t — V?t2. (3)
Это выражение используется в настоящей работе как при 0 < t < т, так и при t > т. При t > т оно, строго говоря, неверно. Однако наиболее важным аспектом при t > т, на наш взгляд, является то, что на этой стадии нагрузка на стенку тела с течением времени уменьшается (до p*t). В аппроксимации нагружения (1) этот аспект учитывается.
Изменение радиуса области нагружения Rl и скорости ее расширения vl = = dR^/dt с течением времени иллюстрирует рис. 3. Видно, что на начальном интервале времени радиус Rl быстро растет от нуля, а скорость vl еще быстрее уменьшается от бесконечности. В последующем скорость vl уменьшается медленнее, соответственно, замедляется и рост радиуса Rl .
76
А.А. АГАНИН, Н.А. ХИСМАТУЛЛИНА
Тело считается изотропным, деформации и перемещения - малыми. Для описания динамики тела используются следующие уравнения:
ди
Рт =
dSrr
~дТ
dSvv
dt
дР _ ~дЬ _
d(Srr — Р) дS,,
дг
+
дг
+
Srr — S,
ду _ clSrZ , д(Szz — Р) S,
2 7 2 ди ду и\
3Р \ дг дг г ) ,
2g ( ди ду 2и А
3 1 дг + дг г )
ди ду и
-К[дТ + дг + г
г
дS.
P8t _
дг
+
дг
zz 2 ди ду и
W _ - з Ч дг — 2 дг + 71’
дSrz (ди ду \
~зт _ р{дг + дт)’
+
г
(4)
Здесь и, у - компоненты скорости по осям г и г соответственно; Srr, Szz, ^‘эфф, Srz - компоненты девиатора S тензора напряжений a, записанные в цилиндрической системе координат; Р - всестороннее (гидростатическое) давление; S _ а + Pg; g - метрический тензор; А _ p(cf — 2е2,) и g _ pc2 - параметры Ламе; К _ А + 2/3g - коэффициент объемного расширения; ci и С2 - скорости распространения продольных и сдвиговых возмущений.
В пластических зонах выполняется условие текучести Мизеса ai _ Yo, ai -интенсивность напряжений,
ai _ (az — ar)2 + (az — аф)2 + (ar — аф)2 + бт^
где az _ Szz — P’ ar _ Srr — P, aф _ Sфф — P, Trz _ Srz - компоненты тензора
напряжений a в цилиндрических координатах.
Если условие текучести нарушается, то есть если ai в некоторой точке расчетной области превышает предел текучести, то девиатор тензора напряжений корректируется [18]
Scorr _ — S. (5)
ai
На поверхности тела г _ 0 полагается
az _ —рь(гр)’ Trz _ О, (6)
где pb(r,t) определяется выражением (1).
В начальный момент времени t _ 0
az _ —po, ar _ aф _ —А(А + 2g)-1po, Trz _ 0, и _ у _ 0. (7)
На бесконечном удалении от места приложения нагрузки параметры тела остаются невозмущенными.
В расчетах бесконечная область, занимаемая телом, 0 < г < ж, —ж < г < 0 заменяется конечной цилиндрической областью 0 < г < г8, hs < г < 0 (рис. 4). Внешние границы г _ г8 и г _ hs являются «искусственными». Они выбираются так, чтобы отраженные от них волны в интересующем промежутке времени не сильно изменяли напряжения в области влияния струи (на рис. 1 эта область выделена прямоугольником). Удаленность искусственных границ зависит от принимаемых условий. В настоящей работе на искусственных границах ставятся неотражающие условия [19], что позволяет размеры расчетной области принять равными г3 _ —hs _ 2.2R [20].
УДАРНОЕ ВОЗДЕЙСТВИЕ СТРУИ ЖИДКОСТИ...
77
Область нагружения / R,. гя
1 т 1
тело 1 1 1 1 1 1 1 1
1 1 1 1
Рис. 4. Расчетная область
Для численного решения задачи (4), (6), (7) применяется классический метод
С.К. Годунова [21]. При этом используется равномерная расчетная сетка с квадратными ячейками с Дг = Az = R/640 . Такая сетка выбрана с учетом сходимости численного решения по мере измельчения ячеек. Шаг по времени определяется из условия устойчивости [21]
S Дг • Az
т =----------
с\ Дг + Az
при S = 0.95.
3. Динамика приповерхностного слоя тела при нагрузках, характерных для воздействия струи
Перед тем, как перейти к рассмотрению динамики приповерхностного слоя тела при нагрузке (1), характерной для воздействия струи жидкости, полезно проанализировать соответствующие одномерные приближения, в которых изменение давления на поверхности тела обладает рядом особенностей нагружения (1). На рис. 5 приведено изменение напряжения az в теле при ударном повышении и последующем резком сбросе давления рь на всей поверхности тела (z = 0, 0 < г < ж):
Рь = Ро при t < 0, рь = р* при 0 < t < т, рь = ро при t > т. (8)
Показаны результаты для Vj = 145 м/с (т = 0.45 нс, р* = 210 МПа) при ti = = 0.51 нс, t2 = 1.29 нс; для Vj = 300 м/с (т = 0.67 нс, р* = 435 МПа) при t1 = 0.78 нс, t2 = 1.68 нс и для Vj = 434 м/с (т = 0.76 нс, р* = 630 МПа) при t1 = 1.29 нс, t2 = 2.07 нс.
При нагружении поверхности тела с амплитудой р* = 210 МПа (наименьшей из трех рассматриваемых) предел текучести материала тела не достигается. Поэтому конфигурация импульса по мере его распространения вглубь тела не изменяется. Пластические деформации в теле возникают при р* > рр = 225 МПа, что соответствует удару струи со скоростью Vj > 155 м/с. Поэтому при нагружении с р* = 300 МПа в теле возникают пластические деформации. В этом случае передний фронт распространяющегося в тело импульса расщепляется на упругий предвестник, перемещающийся со скоростью ci = [(Л + 2р)/р]1/2, и пластическую волну, распространяющуюся со скоростью с* = (K/p)1/2, меньшей в 1.27 раз скорости упругого предвестника. При разгрузке в этом случае пластические деформации не возникают, поэтому задний фронт импульса представляет собой упругую волну, распространяющуюся со скоростью ci . Поскольку ci > c* , область пластических деформаций со временем уменьшается. При нагружении поверхности тела с амплитудой р* = 630 МПа, максимальной среди трех рассматриваемых, пределы текучести превышаются как при нагрузке, так и при разгрузке. Поэтому расщеплению подвергаются как передний, так и задний фронты уходящего в тело
78
А.А. АГАНИН, Н.А. ХИСМАТУЛЛИНА
Рис. 5. Распространяющийся в тело импульс напряжений az при нагружении (8) с р* = = 210 (а), 300 (б) и 630 МПа (в) в указанные в тексте моменты времени ti, t2 (кривые с символами 1,2). Тонкие кривые без символов - точное решение
импульса. Причем как в области переднего, так и в области заднего фронта первым со скоростью ci следует упругий предвестник, а за ним со скоростью c* - фронт пластической волны. В области переднего фронта импульса амплитуда упругого предвестника значительно меньше, чем пластической волны, а в области заднего фронта, наоборот, значительно больше.
Следует отметить, что в момент t2 вся часть тела между его поверхностью и фронтом пластической волны в области переднего фронта уходящего в тело импульса является зоной текучести. Эта зона состоит из двух частей. Первая часть, между передним фронтом пластической волны и задним упругим предвестником, возникает в результате сильного сжатия, а вторая часть, между задним упругим предвестником и поверхностью тела, формируется из-за сильной разгрузки. С течением времени первая часть уменьшается.
Пусть теперь давление на всей поверхности тела (z = 0, 0 < r < ж) изменяется так же, как и в центре конечной области нагружения при законе (1):
pL = po при t < 0, pL = pc(t) при t > 0, (9)
где pc(t) определяется согласно (2). На рис. 6 представлено изменение давления на поверхности тела согласно (9) и приведены пространственные распределения напряжений az и интенсивности напряжений в момент времени t = 0.8 нс для скорости струи Vj = 100 м/с и в момент времени t =1.8 нс для Vj = 300 м/с. В изменении согласно (9) давления pL можно выделить два участка: начальный скачок до p* и последующий монотонный спад до уровня p*t. Имеем p* = 145 МПа при Vj = 100 м/с и p* = 436 МПа при Vj = 300 м/с. В первом случае в теле появляются лишь упругие деформации (интенсивность напряжений не достигает предела текучести Y0 = 125 МПа), так что вглубь тела уходит импульс, форма которого со временем не изменяется. Во втором случае в теле появляются пластические деформации. В этом случае передний фронт уходящего в тело импульса возникает в результате начального ударного повышения давления от po до значения pc(0) = = 436 МПа, превышающего уровень пластических деформаций p*pl = 225 МПа. В результате передний фронт импульса расщепляется на упругий предвестник и следующую за ним волну пластических деформаций. Далее за фронтом этой волны следует участок импульса с упругими деформациями. На этом участке из-за сильной разгрузки интенсивность напряжений сначала уменьшается до нуля, а затем вновь возрастает до значений, близких к пределу текучести.
Рассмотрим теперь динамику приповерхностного слоя тела при нагрузке (1), характерной для воздействия струи жидкости. Рис. 7 характеризует изменение поля
УДАРНОЕ ВОЗДЕЙСТВИЕ СТРУИ ЖИДКОСТИ...
79
Рис. 6. Изменение давления pb(t) на поверхности тела (а, г) по закону (9) и соответствующие импульсы напряжений az (б, д) и интенсивности напряжений ai (в, е) в теле для скорости струи Vj = 100 и 300 м/с (верхний и нижний ряды соответственно) в моменты времени, отмеченные на кривых pb(i) точками
интенсивности напряжений в теле при скорости струи Vj = 100 м/с. В этом случае предел текучести материала тела достигается вскоре после начала ударного воздействия (при t « ti) на краю r = Rb(t{) круговой области нагружения. Затем область, где интенсивность напряжений находится на пределе текучести, постепенно увеличивается, приобретая вид купола, сначала с отверстием в окрестности оси симметрии. Радиус основания этого купола примерно равен радиусу Rl нагруженной области поверхности тела. Постепенно радиус основания этого купола и его толщина увеличиваются, отверстие в окрестности оси симметрии исчезает (рис. 7, а). В момент ts (рис. 7, б) величина нагрузки на периферии области ее приложения (при r = Rl ) достигает своего максимума (рис. 2). После этого к моменту t4 величина нагрузки резко уменьшается (примерно в 2 раза). В результате такой разгрузки область с интенсивностью напряжений на пределе текучести отрывается от поверхности тела и постепенно удаляется от нее (рис. 7, в). По мере удаления размеры этой области уменьшаются (рис. 7, г), а в последующем она полностью исчезает.
Предел текучести в теле достигается, когда максимум нагрузки на его поверхности составляет около 200 МПа, что меньше, чем pp = 225 МПа - значение, необходимое для достижения предела текучести в одномерном приближении (см. предыдущий раздел). Это означает, что возникновение зон с интенсивностью напряжений на пределе текучести обусловлено в данном случае неодномерными эффектами, такими как неоднородность нагрузки, расширение области ее приложения и ее скачкообразное изменение на периферии. Представляется, что последнее обстоятельство оказывает на достижение предела текучести наиболее существенное влияние. В частности, в двумерном случае предел текучести достигается уже при Vj = 65 м/с, тогда как в одномерном приближении - лишь при Vj = 155 м/с.
а) б) в) г)
Рис. 7. Поля Oi в осевом сечении тела при нагрузке (1), соответствующей ударному воздействию струи со скоростью Vj = 100 м/с в четыре указанных на рис. 2 момента времени t2 (а), t3 (б), t4 (в) и t5 (г). В закрашенных областях, границы которых выделены жирной линией, oi = Y0
а) б) в) (г)
Рис. 8. То же, что на рис. 7, но при Vj = 300 м/с и в моменты ti, Ьз, t5 и te
При этом в двумерном случае он достигается уже, когда максимум нагрузки на поверхности тела составляет около 200 МПа, тогда как в одномерном приближении -при 225 МПа.
Рис. 8 характеризует изменение полей интенсивности напряжений ^ при нагрузке (1), соответствующей ударному воздействию струи со скоростью Vj = = 300 м/с. В этом случае предел текучести достигается в самом начале приложения нагрузки. В дальнейшем до момента времени t6 зона с интенсивностью напряжений на пределе текучести расширяется как вдоль поверхности тела в результате расширения области нагружения, так и вдоль оси симметрии из-за распространения импульса вглубь тела. На отрезке 0 < t < t5 эта зона примыкает к поверхности тела, а затем вследствие уменьшения нагрузки отрывается от нее и постепенно удаляется (рис. 8, г). При этом зона с интенсивностью напряжений на пределе текучести не является односвязной.
На рис. 9 для Vj = 100 и 300 м/с представлены области, в которых в промежутке 0 < t < 3т хотя бы в один из моментов времени был достигнут предел текучести. Видно, что при скорости струи Vj = 100 м/c в окрестности центра приложения нагрузки имеется небольшая конусообразная зона, в которой предел текучести не достигается (рис. 9, а). По мере увеличения скорости струи размеры этой зоны уменьшаются, а при скорости Vj = 300 м/с и выше она не возникает (рис. 9, б).
4. Влияние радиальной неоднородности нагрузки
При нагружении (1) радиальное распределение давления в области нагружения при t ~ t3 существенно неоднородно. В частности, в момент t3 давление на периферии области воздействия струи (r = RL) более чем в 2 раза выше, чем в ее
УДАРНОЕ ВОЗДЕЙСТВИЕ СТРУИ ЖИДКОСТИ...
81
б)
Рис. 9. Области, в которых при Vj = 100 м/с (а) и Vj = 300 м/с (б) хотя бы в один из моментов времени промежутка 0 < t < 3т был достигнут предел текучести
Рис. 10. Изменение давления на поверхности тела при радиально-однородном нагружении (11) при Vj = 100 м/с. Представлено шесть последовательных моментов времени, соответствующих рис. 2
центре (r = 0). Для оценки влияния неоднородности рассмотрим эволюцию волн напряжений в теле при соответствующей (1) однородной нагрузке вида
Pl (r, t) = pQ при t < 0, 0 < r < ж, pL = pc(t) при t < 0, (10)
)°.5[pc (t)+ Pe (t)] , Pl (r,t) = <
Ipq 5
0 < r < Rl (t) r > Rl (t)
при t > 0,
(11)
где pc(t), pe(t) вычисляются по (2). Согласно (11) давление pL(r,t) в круговой области нагружения 0 < r < Rl (t) в фиксированный момент времени представляет собой среднее по площади давление, соответствующее закону (1). Изменение давления по формуле (11) при Vj = 100 м/c показано на рис. 10.
Как и в предыдущем разделе, прежде чем перейти к рассмотрению двумерной динамики приповерхностного слоя тела при нагрузке (11), рассмотрим соответствующее ей одномерное приближение. В этом приближении давление на всей поверхности тела (г = 0, 0 < r < ж) изменяется по закону
Pl (t) =
Pq
0.5 [pc (t) + Pe (t)]
при t < 0, при t > 0.
(12)
На рис. 11 представлено изменение давления на поверхности тела по формуле (12), приведены пространственные распределения напряжений az и интенсивности напряжений gi в момент времени t = 0.5 нс для скоростей струи Vj = 100 и 300 м/с. При Vj = 100 м/с максимальное значение pL равно pmax = 210 МПа, а при Vj = 300 м/с оно равно pmax = 632 МПа. В первом случае в теле интенсивность
82
А.А. АГАНИН, Н.А. ХИСМАТУЛЛИНА
Рис. 11. То же, что и на рис. 6, но для нагрузки (12)
напряжений не достигает предела текучести Yo = 125 МПа, так что вглубь тела распространяется импульс, форма которого со временем не изменяется. Во втором случае в теле появляются пластические деформации. В этом случае в изменении давления на поверхности тела можно выделить три участка: начальный скачок от po до рь = 450 МПа, далее монотонный рост до ртах = 632 МПа и последующий плавный спад на уровень рь = 47 МПа. Передний фронт уходящего в тело импульса возникает в результате начального ударного повышения давления до значения рь = 450 МПа, превышающего уровень пластических деформаций рр = 225 МПа. В результате передний фронт импульса расщепляется на упругий предвестник и следующую за ним ударную волну пластических деформаций. Далее за фронтом этой волны находится участок импульса с пластическими деформациями. В нагружении (12) этот участок соответствует промежутку с увеличением давления рь до р-max = 632 МПа. За участком импульса с пластическими деформациями следует участок с упругими деформациями. В нагружении (12) этот участок соответствует промежутку уменьшения давления рь (то есть при разгрузке) от pmax = 632 МПа до рь = 166 МПа. За этим участком импульса с упругими деформациями имеется участок с пластическими деформациями. Он возникает в результате сильной разгрузки (давление рь становится ниже значения 166 МПа, при котором при разгрузке возникают пластические деформации). Таким образом, при нагружении (12) с максимальным давлением ртах = 632 МПа, значительно превышающим предельное для текучести значение р*р1 = 225 МПа, область пластических деформаций в уходящем в тело импульсе не является односвязной (в ней возникает упругое «включение»).
Рис. 12 характеризует изменение полей напряжений в теле при радиально-однородной нагрузке (11) для скорости струи Vj = 100 м/с. Видно, что результаты для однородного (рис. 12) и неоднородного (рис. 7) нагружений во многом подобны. В частности, при однородном нагружении (11), как и при неоднородном (1) (рис. 7), вскоре после начала ударного воздействия возникает зона с предельными
УДАРНОЕ ВОЗДЕЙСТВИЕ СТРУИ ЖИДКОСТИ...
83
а) б) в) г)
Рис. 12. То же, что и на рис. 7, но при радиально-однородной нагрузке (11)
г/к и.. ^ 1 11.4 > 0 r/к и.. д 11.4 5 0 Г/К V.. 3 0 г/к и.: ^mwi
а) б) в) г)
Рис. 13. То же, что и на рис. 8, но при радиально-однородной нагрузке (11)
значениями интенсивности напряжений. Сначала она также имеет вид окружности на поверхности тела радиуса r и Rl , а затем постепенно трансформируется в куполообразную форму (рис. 7, а). По мере расширения области нагружения радиус основания этой куполообразной области с интенсивностью напряжений на пределе текучести и ее толщина (здесь она несколько больше, чем при неоднородном нагружении (1)) постепенно увеличиваются (рис. 12, б, в). С уменьшением нагрузки эта область отрывается от поверхности тела и постепенно удаляется от нее (рис. 12, г). По мере удаления размеры области с предельными значениями интенсивности напряжений уменьшаются, так что в конце концов она полностью исчезает.
Рис. 13 характеризует изменение полей напряжений в теле при радиально-однородной нагрузке (11) для скорости струи Vj = 300 м/с. Видно, что эти поля изменяются примерно так же, как и в при неоднородном нагружении (1). Имеются лишь небольшие отличия в конфигурации области с предельными значениями интенсивности напряжений.
Области приповерхностного слоя тела в окрестности места приложения нагрузки, в которых хотя бы в один из моментов времени промежутка 0 < t < 3т был достигнут предел текучести, при однородном нагружении (11) для Vj = 100 и 300 м/с показаны на рис. 14. Видно, что конфигурации этих областей качественно (а для Vj = 300 м/с и количественно) близки к наблюдаемой при неоднородном нагружении (1) (рис. 9). При Vj = 100 м/с в обоих случаях внутри указанной области (рис. 9, а и рис. 14, а) в малой окрестности центра нагружения имеется небольшая зона, где предел текучести не достигается. Вместе с тем при однородной нагрузке эта зона больше вытянута вдоль поверхности тела, а при неоднородной нагрузке - вдоль оси симметрии.
84
А.А. АГАНИН, Н.А. ХИСМАТУЛЛИНА
б)
Рис. 14. То же, что и на рис. 9, но при радиально-однородной нагрузке (11)
Заключение
Представлены математическая модель и основные положения методики расчета динамики приповерхностного слоя упругопластического тела в случае, когда к его поверхности приложена нагрузка, характерная для импульсного воздействия микроструи жидкости, образующейся при схлопывании касающегося тела кавитационного пузырька. Струя цилиндрическая с полусферическим концом бьет по телу ортогонально его поверхности. Тело моделируется изотропным линейно-упругим полупространством. Его пластическое состояние описывается условием Мизеса. Исследована динамика приповерхностного слоя тела из никелевого сплава (р = = 8000 кг/м3, ci = 5742.9 м/c, с2 = 3069.7 м/c, предел текучести У/ = 125 МПа) при ударе по его поверхности водяной струи (р0 = 1000 кг/м3, со = 1450 м/c, Rj = 20 мкм) со скоростью 65 м/c < Vj < 600 м/c. Показано, что в результате увеличения области нагружения поверхности тела, неоднородности нагружения и скачкообразного уменьшения нагрузки на периферии области нагружения предельные значения интенсивности напряжений в теле достигаются при значительно (примерно в два раза) меньших скоростях струи, чем в соответствующем одномерном приближении. При относительно небольших скоростях в небольшой конусообразной приповерхностной части тела в окрестности центра области приложения нагрузки, соответствующего центру ударного воздействия струи, предельные значения интенсивности напряжений, а значит, и пластические деформации не возникают. При больших скоростях струи зоны с интенсивностью напряжений на пределе текучести не являются односвязными. Неоднородность нагрузки не оказывает качественных изменений в динамике тела, ее влияние с ростом скорости струи уменьшается.
Работа выполнена в рамках Программы РАН № 25П.
Summary
A.A. Aganin, N.A. Khismatullina. Liquid Jet Impact on an Elastic-Plastic Body.
The dynamics of the subsurface layer of a body made of nickel alloy has been studied for the case when the body surface is exposed to a load represented by a pulse action of a liquid microjet arising from collapse of an attached-to-the-body cavitation bubble in water. The velocity range in which the stress intensity in the body briefly attains the level of the yield strength has been considered. The mathematical model of an elastic-plastic body and the fundamentals of the numerical technique used in the study are given. It is shown that due to the extension of the loading domain on the body surface, the non-uniformity of the load and the abrupt decrease of the load at the edge of the loading domain, the stress intensity limit in the body is reached at the considerably (nearly two times) lesser jet velocities than in the
УДАРНОЕ ВОЗДЕЙСТВИЕ СТРУИ ЖИДКОСТИ...
85
corresponding one-dimensional approximation. The non-uniformity of the load does not lead
to qualitative changes in the body dynamics; its influence decreases with an increase in the jet
velocity.
Keywords: cavitation damage, cavitation erosion, jet impact, elastic-plastic body.
Литература
1. Preece C.V., Brunton J.H. A comparison of liquid impact erosion and cavitation erosion // Wear. - 1980. - V. 60, No 2. - P. 269-284.
2. Эрозия / Под ред. К. Прис. - М.: Мир, 1982. - 464 с.
3. Haller K.K., Ventikos Y., Poulikakos D., Monkewitz P. Computational study of highspeed liquid droplet impact // J. Appl. Phys. - 2002. - V. 92, No 5. - P. 2821-2828.
4. Чижов А.В., Шмидт А.А. Высокоскоростной удар капли о преграду // Журн. техн. физики. - 2000. - Т. 70, Вып. 12. - С. 18-27.
5. Obreschkow D, Dorsaz N., Kobel P., de Bosset A., Tinguely M., Field J., Farha M. Confined shocks inside isolated liquid volumes: A new path of erosion? // Phys. Fluids. -2011. - V. 23, No 10. - P. 101702-1-101702-4.
6. Reboud J.L., Guelin P. Impact response of an elastoplastic medium // Mech. Res. Commun. - 1988. - V. 15, No 4. - P. 253-260.
7. Kennedy C.F., Field J.E. Damage threshold velocities for liquid impact // J. Mater. Sci. - 2000. - V. 35, No 21. - P. 5331-5339.
8. Field J.E., Camus J.-J., Tinguely M., Obreschkow D, Farhat M. Cavitation in impacted drops and jets and the effect on erosion damage thresholds // Wear. - 2012. - V. 290-291. -P. 154-160.
9. Sanada T., Ando K., Colonius T. A computational study of high-speed droplet impact // Fluid Dynam. Mater. Process. - 2011. - V. 7, No 4. - P. 329-340.
10. Петушков В.А. Численные исследования нелинейных волновых процессов в жидкости и деформируемом теле при высокоскоростном ударном воздействии // Прикл. механика и техн. физика. - 1991. - № 2. - С. 134-143.
11. Li N., Zhou Q., Chen X., Xu T., Hui Sh., Zhang D. Liquid drop impact on solid surface with application to water drop erosion on turbine blades, Part I: Nonlinear wave model and solution of one-dimensional impact // Int. J. Mech. Sci. - 2008. - V. 50, No 10-11. -P. 1526-1542.
12. Adler W.F. Waterdrop impact modeling // Wear. - 1995. - V. 186-187, Pt. 2. - P. 341351.
13. Hill D.J., Pullin D., Ortiz M., Meiron D. An Eulerian hybrid WENO centered-difference solver for elastic-plastic solids // J. Comp. Phys. - 2010. - V. 229, No 24. - P. 9053-9072.
14. Аганин А.А., Ильгамов М.А., Халитова Т.Ф. Ударное воздействие струи на жесткую стенку // Актуальные проблемы механики сплошной среды. К 20-летию ИММ КазНЦ РАН. - Казань: Фолиант, 2011. - Т. 1. - С. 134-145.
15. Аганин А.А., Ильгамов М.А., Косолапова Л.А., Малахов В.Г. Схлопывание кавитационного пузырька в жидкости вблизи твердой стенки // Вестн. Башк. ун-та. -2013. - Т. 18, № 1. - C. 15-21.
16. Lesser M.B. Analytic solutions of liquid-drop impact problems // Proc. R. Soc. Lond. A. -1981. - V. 377, No 1770. - P. 289-308.
17. Nicolici S., Prisecaru I., Ghitescu P. Study of fluid-structure interaction in liquid droplet impingement phenomena // U.P.B. Sci. Bull. Ser. D. - 2012. - V. 74, No 1. - P. 147-154.
86
А.А. АГАНИН, Н.А. ХИСМАТУЛЛИНА
18. Уилкинс М.Л. Расчет упругопластических течений // Вычислительные методы в гидродинамике. - М.: Мир, 1967. - С. 212-263.
19. Ильгамов М.А., Гильманов А.Н. Неотражающие условия на границах расчетной области. - М.: ФИЗМАТЛИТ, 2003. - 240 с.
20. Аганин А.А., Малахов В.Г., Халитова Т.Ф., Хисматуллина Н.А. Расчет силового воздействия кавитационного пузырька на упругое тело // Вестн. Тат. гос. гуманит.-пед. ун-та. - 2010. - № 4 (22). - С. 6-13.
21. Годунов С.К., Забродин А.В., Иванов М.Я., Крайко А.Н., Прокопов Г.П. Численное решение многомерных задач газовой динамики. - М.: Наука, 1976. - 400 с.
Поступила в редакцию
15.04.14
Аганин Александр Алексеевич - доктор физико-математических наук, профессор, заведующий лабораторией, Институт механики и машиностроения КазНЦ РАН, г. Казань, Россия.
E-mail: [email protected]
Хисматуллина Наиля Абдулхаевна - кандидат физико-математических наук, старший научный сотрудник, Институт механики и машиностроения КазНЦ РАН, г. Казань, Россия.
E-mail: [email protected]
