книги / Моделирование на ЭВМ дефектов в металлах
..pdfгде Cef — термодинамически равновесная концентрация дефектов в материале F.
4.1.Средипе концентрации дефектов
Чаще всего вместо уравнений (24)—(27) пользуются алгеб раическими уравнениями для средних стационарных концентраций Cj. В беспримесном случае (С*.0) при этом имеем (см., например,[20])
С% = С0Н,. ( У Т + Ц - ч0) + С}, |
(34) |
где |
|
со = ]/ i W v |
Hi = Y DJ‘V D^ = Hv - |
Если в материале имеются примеси (пр), может происходить диссоциация (D ) комплексов fa и могут существовать бпваканспп (2v), тогда вклад от каждого из перечисленных факторов (АС) в концентрацию С*, можно представить следующим образом [49]:
с + » |
- |
дс;.|р +A c f + |
07A c f + с;, |
где |
|
|
|
0 . = 1, |
0„ = |
- 1 , АС]* = |
аЯвСвЯ//2(гО<1 |
ACf = аЯ дС аЯ//2
аС„
АС1? = а,2v 1 — [С +С +-ЬС ^(Д С ?-Д С „ПР) +
DA + *>M
+ дс+(дс{, -дсур)],
а величина а2р = V £v0 exp (—з\?jkT), V[ — объем, содержащий по
ложения вакансий около одной выделенной вакансии для образо
вания |
бпваканспп; |
v0 — дебаевская частота; е"‘ — энергия акти |
вации |
миграции |
вакансии; Т — температура; к — постоянная |
Больцмана. |
|
|
В случае выпадения из твердого раствора, например, когерент |
ных предвыделений с ег < 0 (сжатых) при облучении материала генерируемые межузлия внутри предвыделений (I) будут стре миться выйти в растянутую область (II), т. е. можно записать для усредненных эффективных скоростей появления дефектов
el “ в1. |
eJ = eI + 6„. |
— ^ + |
г" « я 11, |
= (l ~~vp) 1J Ip (^р) fp (^р) ^р» |
(30) |
||
& v ~ vp |
[ Ip * (Я р ) f p (^р) d R p . |
(37) |
Тогда для стационарных концентраций CJ11 получим
131 |
9* |
С +И = |
C $H t ([/1 + |
Л 2 - i i ) + |
Cet, |
(38) |
|
c + " = |
(Ctn _ Щ |
D t f / D j l - |
b j D vk l + C®, |
||
|
c j 1 = Y |
$ h iDi’ |
^ = |
ч 0 0 |
— pbilDJ%*t) Y |
g II^ ' ‘I • |
|
|
При малой эффективности рекомбинации дефектов |
из (38) имеем |
||||||
С у 1 Cej 4~ gy/Djk}, |
а |
из |
(34) |
получаем |
Гл II |
C'i + е 'Ч о м . |
|
°уо |
|||||||
Следовательно, в нашей модели можно записать |
|
||||||
С+“ : |
c j 0 + V s" ) + с;: |
с+" |
45 |
*v' |
(39) |
Рассматривая средние нестационарные концентрации дефектов в областях I и II п учитывая влияние когерентных предвыделений, также удается показать перераспределение вещества между мат рицей и формирующимся выделением [50].
4.2. Распределения точечных дефектов вблизи формирующихся выделений второй фазы
Стационарные распределения точечных дефектов как внутри предвыделений, так и в матрице получим из уравнений (24)—(27) при следующих предположениях: 1) внутри формирующихся вы делений точечные дефекты не образуют кластеров и у них нет дру гих стоков, кроме поверхности выделения; 2) стоки для точечных дефектов распределены в областях II равномерно и к у 1=Щ]
в) рекомбинация одиночных вакансий п межузлий мала. Тогда имеем [51]
СИ'> = С/ |
Si |
•*). |
|
|
+-6Т57(* |
|
|||
|
|
|
R P sh [{Rp + Lp — r)kj) |
|
с " {г)=~с," + - щ |
г |
sh (Lpkj) |
||
, ■gen |
RP |
sh (rkj) — sh [(r — RP) kj] |
|
|
_rL ' |
r |
sh (Rpkj) |
|
(40a)
(406)
Из (40) видна зависимость CF не только от г, но и от энергии
взаимодействия Е 3р |
через Се/ |
■ Так, например, |
езли выделение |
||
сжато, т. е. ех< 0 , |
то Sp3p > 0 , |
S p cJ^ O . При этом |
# р < 0 , |
||
Ее > 0 , Е р ' > 0 , |
Е ‘р" < 0 , |
а |
значит х; > 1 , |
Х; < 1 , |
х Г < |
xj1 > 1. Следовательно, в предвыделение будут поступать преиму щественно вакансии, а из него выходить межузельные атомы, тем самым повышая в областях II пересыщение по межуалиям и спо собствуя усилению рекомбинации разноименных дефектов в этих областях.
132
В твердых растворах, находящихся в состоянии распада под облучением, когда выделения обособились, для описания их эво люции необходимо знать распределения точечных дефектов около них. Этим вопросам посвящен ряд работ (например, [16—18, 52—56 J). Обычно различают два предельных случая: 1) диффузия одиночных атомов примеси к выделениям; 2) диффузия атомов при меси в составе комплексов (/а).
Рпс. 5. Схема решений кинетических урапненпн баланса.
Если энергия связи комплекса е*0 < е™, т. е. время жизни комп
лекса ja порядка времени его перескока на одно межатомное рас
стояние, то концентрации атомов примеси С0, собственных точечных дефектов С . и комплексов CJa находятся в локальном равновесии,
т. е.
aCj(r, t)C a (r, t) = xCj a (w, t). |
(41) |
В этом случае система уравнений (24)—(27) расщепляется на две не связанные друг с другом системы (рпс. 5), одна из которых опи сывает эволюцию собственных точечных дефектов, а вторая — диффузию одиночных примесей и комплексов ja к своим стокам.
Поскольку в реальных твердых растворах концентрация атомов
133
Л Т Р р Т Г fpX* + Y V(pX,lY ] = bL p { <*)•
В работе [44] подробно изложена процедура получения реше
ния задачи (44)—(45), зная которое, запишем |
выражения для кон |
центраций дефектов: |
|
для медленного истощения |
|
С,а № ) « ( > - - J L)lC )A * + DiJ?ia + A ’ V |
P M ' - '.) /'* ] : <«) |
для быстрого истощения при больших временах t > te и в пре
небрежении диссоциацией комплексов и рекомбинацией разно
именных |
дефектов |
|
|
|
|
Cj (г, |
t) & c e. + g t |
— Сд (1 — ~f~) *~t/tRp , |
|
(47а) |
|
С,-а {г, 0 « Ceja + |
С°а (1 - |
е~Ч'*Р . |
|
(476) |
|
Здесь Cja— термодинамически |
равновесная |
концентрация |
ком |
||
плексов; время истощения tc = |
(аС} + Djakja+ |
y)laCjDjaa]a; |
Cl — |
начальная концентрация примеси; tpp — LP(LP— R pfjSD ^p; 2Lp—
расстояние между преципитатами; |
1 = |
DJa[р.®- г1(Ьр — i?/>)]2; {AJ- 1 |
удовлетворяет уравнению tg р.®*1 = |
|а°> 1LPl(Lp — R P); RP— средний |
|
радиус преципитата. |
|
|
4.3. Распределения точечных дефектов |
около дислокаций |
|
и скоплений |
|
|
В распадающемся твердом растворе рассмотрим распределе ния точечных дефектов около краевых дислокаций п скоплений точечных дефектов (например, мелкпх дислокационных петель). При этом удается учесть налпчпе когерентных предвыделенпй, создающих упругие поля напряжений впда (28).
Аналогично работе [43] в [58] решена задача по отысканию ста
ционарных |
распределений |
дефектов около |
краевой дпслокацпп |
|
в присутствии когерентных предвыделенпй в впде |
|
|||
|
|
|
|
(48) |
Здесь Up (г, |
G) = E D ( г , Ь)/кТ = Щ sin 6/г; Ер (г, 8) — энергия |
изаи |
||
модействня дислокации с точечным дефектом |
/ х — R^ 2r, |
Хр — |
||
— ^ /2 г0; |
= gRflAD ji |
|
|
|
Фщ (*) = |
OliD + Лги) V0 (■*) — a0£)pfl (.v); |
|
|
|
у / ч |
7о М «о М ~~ {0 (хр ) К» (А‘г>) , |
|
|
|
lo (XR) « о(XD) ~ |
(*D) К0(лл) ’ |
|
|
135
ные, зависящие от параметров материала п задачи, Li (z) — ин тегральный логарифм.
Таким образом, в разделе 4 получены аналитические выражения для стационарных концентраций дефектов при учете различных факторов распределения точечных дефектов как около формирую щихся выделений второй фазы, так и около краевых дислокаций и скоплений в материале с когерентными предвыделениями.
5. Скорости и эффективности поглощения точечных дефектов стоками
Подставляя полученные выражения для распределений де фектов около когерентных предвыделенпй (40), краевых дислока
ций (48) и скоплений (50) в (17), |
а затем в (15) и |
(16), |
запишем |
|||||||||||||
|
|
|
|
|
|
|
+ Ч ) < |
|
|
|
|
|
|
|
<51> |
|
«Ь (О = <4> Di № |
|
+ |
*i>). |
|
|
|
|
|
|
|
<52> |
|||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ACT = c f - |
С% |
|
|
|
|
|
|
|
|
|
|
|
|
|||
«/„ = |
4n, 4 |
= |
( Ц + |
|
|
(1 |
- |
^ |
|
|
- |
C |
f ( ' - |
|||
kjLP |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
sh k'-Lp /) ’ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
2л |
|
- |
, |
g^ ' [ p (x |
x |
:v* |
|
1 |
|
l _ c ? |
||||
__________ к/ |
|
|
|
|||||||||||||
OLL — |
ln(2 * / * { ) ’ |
D |
|
4 ^ 1 |
o(-J?) |
kjR |
dx | ^ |
R| |
1- |
|||||||
D “ |
|
|
||||||||||||||
v)L = |
4лzlL, |
zvL « 1, |
zlL « |
3 p ' | E{r (Rl ) IIkTI Г |
* |
|
|
|
||||||||
/ _ |
|
grl |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H^~ 6 D T ’ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
R — радиус |
области |
|
I |
около |
|
дислокации, |
где |
|
(R )/kT < 1, |
|||||||
но еще пе выполняется условие E iJkT |
1. Добавки х> определя |
ются неоднородностью распределений дефектов / вблизи стока q.
В случае сжатого предвыделенпя выражение (51) имеет смысл скорости поглощения вакансий, в то время как для испускания межузлий подстановка (40а) в (15) дает
jv |
l p t ( R p ) ~ S i ^ p (^р)> |
где |
Vp (Rp) — объем предвыделенпя радиуса R P. |
При сильном взаимодействии точечных дефектов с дислокацией, т. е. при больших значепиях XD=R{J2r0 ^> 1, и с учетом того, что
137
захватов межузельных атомов в областях между предвыделениямп, так что на дислокации из области II уходит меньше межузлий, чем вакансий. Такие зависимости в материале с когерентными пред-
выделениями могут приводить с ростом vp к изменению знака раз
ности КЪ (г*)—КЪ (г*), которая определяет вклад дислокаций в из менения свойств материалов при облучении.
Приведем также при |
е, = |
—0.02, |
g-=1016 см-3 |
• с"1, |
||||
pD= 1 0 10 см-2 вклады vJh |
в сравнении |
с |
ACJ: |
х^/Д£+ = |
72.2 %, |
|||
x£(bC t= 7& %. Отсюда видно, что нельзя |
пренебрегать |
слагае |
||||||
мыми y.DJ в (52). |
|
|
|
|
|
|
|
|
Для скоплений дефектов при тех же |
параметрах материала |
|||||||
и облучения можно получить z£« 1.22, |
|
а прп |
Щ=Щ можно |
|||||
оцепить |
роль вкладов v.j : xjjkCt = 36.3 %, |
|
= |
36.4 %. |
Это |
|||
означает, |
что, как и для |
дислокаций, роль |
слагаемого xj |
в |
(51) |
существенна.
Таким образом, удалось получить в случае стационарных рас пределений около стоков различной природы общие выражения (51) и (52) для таких важных характеристик материала, как ско рости и эффективности поглощения стоками точечных дефектов. Оценки показывают также существенный вклад неоднородных до бавок х£.
6. Эволюция предвыделенпй второй фазы под облучением
Полученные для концентраций дефектов вблизи от формирую щего выделения выражения (46) п (47) позволяют, пользуясь уравнением (13), получить соответствующие аналитические зави симости для радиусов выделений R? как функций времени.
В случае медленного |
истощения |
твердого раствора (te > t%) |
из (46) и (13) получим |
|
|
Rp (0 да |
1 —ехр |
(53) |
где BJtl= a0 (1-|-|)2)^в, Д р* — радиус |
выделения в момент t= t |
а £ описывает стехиометрический состав формирующихся выделе ний (т. е. на один атом примеси приходится | атомов матрицы), при чем параметр £ учитывает также разность удельных объемов мат рицы и выделения и возможные присоединения одиночных соб ственных точечных дефектов. Результаты численного расчета по (53) для твердого раствора Си—Be п экспериментальные данные из [61 ] приведены на рис. 7.
В случае быстрого истощения (tc < f*), подставляя (476) в (13),
имеем
139
Яо |
(f) л* " |/”Я ре + 2В^а Сд (/ |
*с)> |
(54) |
где R PC= R P (tc)- Сравнивая результаты для сплава А1—Zn, |
полу |
||
ченные |
в работах [17, 18], |
с (54), можно заметить, что при ком |
натных температурах выражение (54) хорошо описывает экспери менты на начальной стадии старения сплавов (см. [17]). Поэтому некорректность приближения локального равновесия (см. [62]) осталась незамеченной в [17]. Между тем отказ от условия локаль
ного |
равновесия |
(41) |
приводит при |
повышенных температурах |
||||||||||
|
|
|
|
|
|
|
|
к более быстрому росту преципи |
||||||
|
|
|
|
|
|
|
|
татов, |
поскольку |
комплексы |
по |
|||
10 |
|
|
|
|
|
|
пути |
следования |
к |
выделению |
||||
|
|
|
|
* |
* |
практически |
не диссоциируют. |
|
||||||
|
|
|
|
|
|
|
||||||||
|
|
/ |
|
I |
|
|
|
|
|
|
|
|
||
£ |
5 |
|
|
|
|
|
|
|
|
|
|
|
||
' |
|
1 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
1 |
|
|
|
Рис. 7. Зависимость радиуса выделе |
||||||
tn |
|
|
|
1 |
|
|
|
|||||||
|
|
|
1 |
|
|
|
ния R p от времени t для неоднородного |
|||||||
|
2 |
|
|
1 |
|
|
|
распределения |
по объему |
комплексов |
||||
|
|
|
1 |
|
|
|
атом |
прпмесп—межузельный |
атом |
|||||
1П1 __I |
|
1__ 1 |
К |
' |
1 |
|
в твердом растворе |
Си—0.18 вес.% |
Be |
|||||
|
■' |
|
при Г=583 К |
и |
g=4*1019 см~3-с-1 |
|||||||||
|
10г 2 |
|
5 10s 2 |
5 |
10* |
|
|
(Ф=2.9*1018 эл-см"2-с"1). |
|
|||||
|
|
|
Ig |
t |
(c) |
|
|
Сплошная линия — расчет по |
(53); |
точки — |
||||
|
|
|
|
|
|
|
|
экспериментальные данные LG1]. |
|
7. Распад [твердых растворов и радиационное распухание
Изучение радиационного распухания материалов как экс периментально, так и теоретически ведется уже достаточно давно (см. литературу к [6, 7]). Получен обширный фактический мате риал, разработаны отдельные модели описания этого явления, по строены основы теории. Однако, как справедливо отмечено в ра ботах [7, 63], ключевые вопросы теории радиационного распу хания требуют дальнейших исследований. Теоретические иссле дования в этой области в основном касались однородных материа лов. Экспериментально же известно [2—5, 45, 46, 64, 65], что на личие примесей и выделений второй фазы на различных стадиях их формирования и эволюции может приводить к снижению распуха ния материалов. В связи с этим в продолжение работы [5] были обобщены постановка задачи для теоретической оценки скорости
радиационного распухания S в гетерогенном материале и про
цедура ее расчета [12, 13] (см. рис. 1).
Радиационное распухание материалов обусловлено поступле нием на «далекие» (дислокации, границы, внешние поверхности) и «близкие» (кластеры собственных дефектов, примеси, выделения второй фазы) стоки выбитых облучением межузельных атомов, которые увеличивают объем материала на AV. Радиационное поро
образование развивается в материале, когда имеется избыточная концентрация нескомпенсированных вакансий, поскольку меж-
140