Таким образом, разрушение металла может вызвать его наводороживание и в конце концов новое водородное охрупчивание. Естественно, что происходят и обратные процессы разводороживания металла, которые могут быть прерваны на любой стадии наводороживания.
В наших работах [10, 11] на примере системы колесо - рельсы железнодорожного транспорта исследовались процессы наводороживания поверхности катания колеса, распределение диффузионно-подвижного водорода по глубине и возможности предотвращения неизбежного, но более отдаленного по времени эксплуатации, наводороживания колесной стали. Предлагаемые в этих работах технологические приемы позволили увеличить пробег колеса в 3 раза.
ЛИТЕРАТУРА
1. Варенцов Е.А., Хрусталев Ю.А. Механоэмиссия и механохимия молекулярных кристаллов органических веществ // Успехи химии.
1995. Т. 64. № 8. С. 834-849.
2. Гай А.П., Глазунов М.П., Журавлев Л.Т., Хрусталев Ю.А., Шенге-лия К.Я. О причине газовыделения в процессе механического раз-
рушения твердых тел // Доклады АН СССР. 1984. Т. 277. № 2.
С. 388-390.
3. Хрусталев Ю.А., Симаков Ю.С., Глазунов М.П., Губин В.В. Образование водорода при сухом трении металлов // Журнал физической химии. 1989. Т. 63. № 5. С. 1355-1357.
4. Симаков Ю.С., Коган Б.М., Ляхов Б.Ф., Тененбаум М.М., Хрусталев Ю.А., Кауфман С.М. Наводороживание металла при фрикционном взаимодействии с биологической массой // Трение и износ. 1989. Т. 10. № 4 .С. 749-750.
5. Хрусталев Ю.А. Физико-химическая концепция наводороживания металлов. Роль эмиссии электронов с поверхности разрушения конструкционных материалов в процессе образования водорода // Эффект безызносности и триботехнологии. 1997. №2. С. 19-35.
6. Розенфельд И.Л., Персианцева В.П. Ингибиторы атмосферной коррозии. М.: Наука, 1985. 278 с.
7. Белоглазов С.М. Наводороживание стали при электрохимических процессах. Л.: ЛГУ, 1973. 411 с.
8. Лебедев В.В. Физико-химические основы процессов получения водорода из воды. М.: Наука, 1969. 134 с.
9. Пикаев А.К. Сольватированный электрон в радиационной химии. М.: Наука, 1969. 457 с.; его же. Современная радиационная химия. Радиолиз газов и жидкостей. М.: Наука, 1986. 440 с.
10. Хрусталев Ю.А., Ляхов Б.Ф., Балабанов В.И., Мамыкин С.М. Наводороживание поверхности качения колесной пары // Вестн. машиностроения. 1997. № 11. С. 23-26.
11. Хрусталев Ю.А., Ляхов Б.Ф., Гребенюк М.Н., Терегеря В.В. Влияние металлоорганической присадки Гретерин-З на наводороживание стали при трении // Вестн. машиностроения. 2000. № 5. (В печати.)
УДК 539.4:537.5
МИГРАЦИЯ ТОЧЕЧНЫХ ДЕФЕКТОВ У ВЕРШИНЫ ТРЕЩИНЫ ПОД ДЕЙСТВИЕМ МЕХАНИЧЕСКОГО НАПРЯЖЕНИЯ И ЭЛЕКТРИЧЕСКОГО ТОКА
© Д.Н. Карпинский, С.В. Санников
Россия, Ростов-на-Дону, НИИ Механики и прикладной математики при Ростовском Госуниверситете
Karpinsky D.N., Sannikov S.V. The electrical current influence on the point defects migration at the crack tip. The estimation is realised of the constant electrical current influence on the point defect migration at the crack tip in an alpha-Fe crystal. The computation takes into account the plastic deformation at the crack tip in the loaded sample, Joule heat and the gas-exchange in crack banks. The plastic deformation is caused by dislocation motion in the active sliding crystal planes. Stress relaxation and asymmetry are revealed in point defect distribution induced by an electrical current at the crack tip.
1. В последние годы достоверно установлено наличие резкого снижения сопротивления металла деформированию и повышение его пластичности под влиянием электрического тока большой плотности [1]. В связи с этим представляет значительный интерес исследование эволюции точечных дефектов у вершины трещины при совместном действии сил механического и электрического происхождения. Временная зависимость распределения точечных дефектов в отсутствие электрического тока оценена в [2]. Однако для разработки методов электрохимического наводороживания [3] и других задач важно понимать закономерность влияния электрического тока на прочность образцов, которая существенно зависит от миграции точечных дефектов в них. В этом направлении исследований авторам известна только работа [4], в которой выполнен расчет эволюции распределения атомов водорода, растворенных вблизи надреза в токонесущем стержне, но в этом расчете предполагалось, что внешняя нагрузка не вызывает пластической деформации, и отсутствует накопление газообразной примеси в полости надреза.
Цель настоящей работы заключается в расчете эволюции распределения атомов водорода у вершины тре-
щины в кристалле при решеточной диффузии под действием постоянного электрического тока и механической нагрузки с учетом газообмена на берегах трещины.
2. Исследование эволюции концентрации с(г, Г) ме-жузельных атомов примеси и соотношения вкладов дислокационных механизмов в процесс переноса примесных межузельных атомов выполнено в [2]. В этой работе рассмотрена трещина длиной I, расположенная в плоскости скола (010) бесконечного кристалла с объемно-центрированной кубической (ОЦК) решеткой. К плоскостям кристалла у^+<х приложено однородное напряжение растяжения ст„(Г) (мода I), монотонно
возрастающее до некоторого значения стшах пластического деформирования кристалла, но недостаточного для роста трещины. Дополним [2] предположением о протекании через кристалл перпендикулярно плоскости трещины постоянного электрического тока плотностью /(г) в удаленной от трещины части кристалла (/х(г) = 0, ]'у(г) = /о при г^<х). Расчет эволюции пластической деформации у вершины трещины аналогичен [2] за исключением вычисления тепловыделения в пластической зоне у вершины трещины. Будем предполагать, что трещина расположена вдоль оси абсцисс, а ее
вершина - в начале координат. Температура Т в верхней полуплоскости в момент времени Г определяется по формуле [5]:
Т (х,у,г) = | йу1 | Т0 (х1,у1 )О2 (x,y,x1,y1,г)dx1 -
0 -да
2 Р І* (1)
- а I ^ \ Ф(х1,т)О2(х,у,х1,у1 = 0, г - г)dX1 +
0 -да да да да
+| йу11Е(х1,у1,г)О2(х,у,х1,у1, г - т)dx1 .
0 0 -да
В (1) О2 - функция Грина.
02(х,у,х1,у1,г) = 0(х,у,х1,у1,г) + 0(х,у,хь- уі,г),
^ 1 \ (х - х1)2 + (у - у1)2
0(х,у,х1,у1,г) = д л ехр<! ---------------------------------
4%кі
4а 2г
а = X / рС - температуропроводность кристалла, X -теплопроводность, р - плотность и С - теплоемкость,
^(х,у,Г) = р/(г) + //'(г), Р = Р0(1 + а(Т- 273 К)) - удельное электросопротивление кристалла, а - температурный коэффициент электросопротивления, / - коэффициент Томсона [6]. Используя симметрию задачи относительно у = 0, получим граничное условие Ф(х,Г) = дТ
= --- = 0. Выражение для плотности электрического
дУ
тока /(г) определяется аналогично [7, 8], но в нашем случае ограничимся выражением для ее асимптотики
вблизи вершины трещины, /х = /0 VI / 2г 8Шф, ]у =
= /0 VI / 2г С08ф, где ф - полярный угол.
3. Рассмотрим теперь основные уравнения электромеханодиффузии [4, 9]:
дс -_.т т Ос ^
= -У/, У =——Уд,
дг
квТ
(2)
где для слабого раствора примесей в упругом поле напряжений химический потенциал ц = квТ1п(с/с0) -
- и(г,г), О - коэффициент диффузии, и(г,г) = Ауо(г,г) -
- е2*ф, где Ау - изменение объема кристалла за счет искажения кристаллической решетки межузельным атомом водорода, ст(г,г) - шаровая компонента тензора напряжений, е - элементарный заряд, 2* - эффективное зарядовое число, а ф - электрический потенциал. Тогда для диффузионного потока:
У = О Ус-—Уо
квТ
(3)
Из (2, 3) получаем (учитывая, что Аи = 0 [4]) дс О
— = оу2 с( г, г) +--Ус( г, г )Уи (г, г). (4)
дг квТ
Начальное условие выбрано в виде с = с0 при Г = 0, дс
а граничные условия О— = кт (с - (Г')-Р(Г))
дУ
при у = 0, х < 0 (кт - коэффициент массообмена на
границе газовой и твердой фазы, Г' - модифициро-
дс е2* /0рс II
ванная постоянная Зивертса [2]) и — =-------------л —
ду квТ V 2г
при у = 0, х >0.
Обсудим теперь подробнее выбор граничных условий при у = 0, х < 0. Предположим, что газ в полости трещины идеальный, тогда его давление Р(Г) внутри трещины равно:
Р(г) = (^а'(Г)/2)[(1 + 40ТМ(Г) / [п(1 - v)(2/)2а'а2(Г)])1/2 - 1], (5)
где О - модуль сдвига, V - коэффициент Пуассона, N -количество молекул газа в полости трещины. Предположим, что основной поток газа, растворенного в объеме кристалла, попадает в полость трещины через ее
вершину. В этом случае величина Щ(Г) = 13(Г)ёГ , где
3(Г) - поток атомов газа, растворенных в кристалле, которые переносятся к вершине трещины.
Метод решения (4) - (5) аналогичен [2], а для коэффициента интенсивности напряжения (КИН) трещины имеет место представление [2]:
К(г) = Кс(г) + кр(г),
(6)
где Кс(г) - слагаемое КИН трещины, не учитывающее влияния на нее пластической деформации, а Кр(г) -определяющееся исключительно этим влиянием (выражение для Кр(г) см. в [2]).
Расчеты выполнены для кристалла а - Ее при следующих значениях постоянных: 21 = 10-3 м, Ау = = 2,06-10-30 м3, О = 60 ГПа, V = 0,3, р0 = 8,6-10-8 Ом-м, Ь = -22,8-Ю-6 В/К, а2 = 2-10-3 м2/с, а = 3,3-10-3 К-1, О = 4,88-10~12 м2/с, кт = 4,88-10-15 м/с. [10]. Скорость нагружения кристалла выбиралась так, чтобы максимальная скорость деформации єшах в пластической зоне была равна 0,1 с-1. После достижения внешней нагрузкой своего верхнего предела стшах = 5 МПа в дальнейшем она оставалась постоянной.
Оценкой процесса релаксации напряжения в пластической зоне является эволюция КИН. Эта зависимость вычислена согласно (6) и показана на рис. 1.
Распределение концентрации атомов водорода у вершины трещины в последовательные моменты времени показаны при двух величинах плотности тока на рис. 2, 3.
В заключение приведем основные выводы по результатам расчетов:
1. Постоянный электрический ток является дополнительным механизмом релаксации напряжений у вершины трещины.
2. Увеличение плотности тока приводит к перераспределению концентрации межузельных атомов водорода у вершины трещины, обуславливая асимметрию ее распределения относительно линии продолжения трещины.
3. Величина плотности электрического тока слабо влияет на сток атомов водорода в полость трещины и ее нагружение за счет «распирания» берегов газом.
Рис. 1. Временная зависимость коэффициентов интенсивности напряжения КИН: Кс - хрупкой трещины, К - пластичной трещины с учетом протекающего через образец постоянного электрического тока 5-108 А/м2 (сплошная линия) и без него (пунктир)
У'/""'
2- Ц-
и Г
\ /
J Ц Эсуит
8-
/ /-
^2
$
- ц -
Рис. 2. Распределение безразмерной концентрации с(г,/)/с0 межузельных атомов водорода при решеточной диффузии при j0 = 108 А/м2 в моменты времени а) I = 0,27 с; б) t = 5,7 с. На а) 1 - 0,997, 2 - 1,01, 3 - 1,03, 4 - 1,04, 5 - 1,05; б) 1 - 1,02, 2 - 1,04, 3 - 1,07, 4 - 1,10
Рис. 3. Распределение безразмерной концентрации с(г^)/с0 межузельных атомов водорода при решеточной диффузии при jo = 5-10 А/м в моменты времени а) t = 0,14 с; б) t = 4,1 с. На а) 1 - 0,932, 2 - 0,957, 3 - 0,983, 4 - 1,01, 5 - 1,03, 6 - 1,06; б) 1 - 0,936, 2 - 0,967, 3 - 0,999, 4 - 1,03, 5 - 1,06
ЛИТЕРАТУРА
1. Спицин В.И., Троицкий О.А. Электропластическая деформация металлов. М.: Наука, 1985. 159 с.
2. Карпинский Д.Н., Санников С.В. Миграция точечных дефектов у вершины газонаполненной трещины // ФММ. 1998. Т. 85. № 1. С. 10-18.
3. Белоглазое С.М. Наводораживание стали при электрохимических процессах. Л.: Изд-во Ленингр. ун-та, 1975.
4. Угодчиков Н.А., Берендеева А.Г. Численное моделирование низкотемпературной нестационарной направленной диффузии водорода // Прикладные проблемы прочности и пластичности. Методы решения: Межвуз. сб. М.: Тов-во научн. изд. КМК, 1997. С. 69-74.
5. Будак Б.М., Самарский А.А., Тихонов А.Н. Сборник задач по математической физике. М.: Наука, 1972. 687 с.
6. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошной среды. М.: ГИТТЛ, 1957. 532 с.
7. Максимов И.Л., Свирина Ю.В. Диссипативные неустойчивости разрушения в проводящих материалах с транспортным током.
1. Критерий неустойчивости, качественный анализ // ЖТФ. 1996. Т. 66. № 9. С. 64-74.
8. Максимов И.Л., Свирина Ю.В. Диссипативные неустойчивости разрушения в проводящих материалах с транспортным током.
2. Эволюционные уравнения, диаграммы неустойчивости // ЖТФ.
1996. Т. 66. № 9. С. 75-85.
9. Герцрикен С.Д., Дехтяр И.Я. Диффузия в металлах и сплавах в твердой фазе. М.: Физматгиз, 1960.
10. Физические величины. Справочник / Под ред. И.С. Григорьева и Е.З. Мейлихова. М.: Энергоатомиздат, 1991. 1232 с.