Управление микроструктурой силуминов термической обработкой в твердо-жидком состоянии Текст научной статьи по специальности «Технологии материалов»

Journal of Siberian Federal University. Engineering & Technologies 5 (2015 8) 582-593

УДК 669.017

Manipulation of Silumin Microstructure in Liquid-Solid Condition, Using the Heat Treatment

Aleksey I. Anikin*, Sergey V. Belyaev, Vladimir P. Zhereb and Valentina I. Anikina

Siberian Federal University 79 Svobodny, Krasnoyarsk, 660041, Russia

Received 30.03.2015, received in revised form 21.05.2015, accepted 12.07.2015

In this paper we examine the effect of heat treatment on the structure of the eutectic aluminum alloys, a theoretical study process on the basis of the microstructures. The phenomenon of clustering of liquid metal solutions and develop in following this process, heterogenization solid phases in chemical composition.

Keywords: silumin, Al-Si, eutectic, heterogenization, dendrites, microstructure. DOI: 10.17516/1999-494X-2015-8-5-582-593.

Управление микроструктурой силуминов термической обработкой в твердо-жидком состоянии

A.И. Аникин, С.В. Беляев,

B.П. Жереб, В.И. Аникина

Сибирский федеральный университет Россия, 660041, Красноярск, пр. Свободный, 79

В статье рассмотрено влияние режимов термической обработки на структуру эвтектических алюминиевых сплавов, проведено теоретическое обоснование процесса на основе полученных микроструктур. Изучено явление кластеризации жидких металлических растворов и развивающаяся вслед за этим процессом гетерогенизация твердых фаз по химическому составу.

Ключевые слова: силумин, Al-Si, эвтектика, гетерогенизация, дендриты, микроструктура.

© Siberian Federal University. All rights reserved

* Corresponding author E-mail address: anikin24081989@mail.ru

Введение

В статье рассмотрено явление кластеризации жидких металлических растворов и развивающаяся вслед за этим процессом гетерогенизация твердых фаз по химическому составу. Использование этого явления позволяет разделить процесс кристаллизации жидких растворов на два этапа: диффузионное разделение и непосредственно кристаллизация гетерогенных кластеров.

Эвтектические алюминиево-кремниевые сплавы широко используются в автомобильной и авиационной промышленности благодаря своей низкой плотности и низкому коэффициенту теплового расширения. В зависимости от химического состава при традиционном способе кристаллизации алюминиево-кремниевого сплава могут образоваться как иглообразные кристаллы кремния, так и крупные первичные кристаллы кремния и эвтектические аккумулятор -ные плиты, что проводит к резкому увеличению хрупкости сплавов. Следовательно, возникают сложности, если требуется улучшить микроструктуру после литья.

Двойные металлические расплавы, или двойные жидкие растворы металов, - это макроскопические однородные конденсированные фазы, характеризующиеся атомным диспергированием компонентов и ближним порядком в расположении атомов. В подавляющем большинстве двойных металлических систем компоненты проявляют неограниченную взаимную растворимость в жидком состоянии. Если такая же растворимость наблюдается и в кристаллическом состоянии, т.е. существует непрерывный ряд твердых растворов от чистого компонента А до чистого компонента Б, то затвердевание расплава представляет собой однофазную кристаллизацию: жидкий АБ-раствор любой концентрации переходит в а-твердый раствор.

Часто растворимость компонентов в кристаллическом состоянии ограничена, и если концентрация расплава превышает предел растворимости, то его затвердевание осуществляется путем многофазной кристаллизации. Одним из распространенных типов такой кристаллизации является эвтектическое превращение: диффузионное разделение расплава на две, кристаллические фазы образующиеся при одной температуре.

Необходимой предпосылкой эвтектического превращения служит, таким образом, отсутствие полной растворимости компонентов в кристаллическом состоянии.

Эвтектика и дендритные структуры сплавов

В эвтектических сплавах уже в жидком состоянии имеются области, обладающие структурой чистых компонентов, т.е. вблизи эвтектической точки можно проводить аналогию между структурой жидких и твердых эвтектических сплавов [1].

Смешение жидких компонентов - энергетически невыгодный процесс, который поэтому и не идет до конца. Расплавившиеся кристаллики компонентов эвтектики при температуре, немного превышающей температуру плавления, сохраняют до некоторой степени свою индивидуальность [2].

Таким образом, в подавляющем большинстве случаев дендритная кристаллизация сопровождается образованием эвтектических включений между ветвями дендритов. Можно даже сказать, что при полном затвердевании сплава дендритная кристаллизация почти всегда заканчивается кристаллизацией эвтектики [3].

Кристаллизация эвтектики

Рассмотрим вопрос о кинетике кристаллизации бинарных сплавов эвтектического типа.

При анализе эвтектического структурообразования прежде всего используют диаграммы состояния. Это не всегда делается корректно в том отношении, что линии термодинамического равновесия фаз принимают за линии начала или окончания кристаллизации. Описание кристаллизации с помощью равновесной фазовой диаграммы должно предусматривать возникновение температурно-концентрационных условий, являющихся отклонением от равновесных и нарушающих равенство химических потенциалов компонентов в фазах, что стимулирует фазовый переход. Эти отклонения связаны:

• с торможением диффузионного перераспределения компонентов между фазами и в фазах, заключающихся в образовании неоднородных твердых растворов;

• с появлением эвтектической структурной составляющей в сплавах, в которых еще не достигнут предел растворимости;

• с увеличением доли эвтектики против рассчитываемой для данного сплава по диаграмме состояния.

В морфологическом плане важно то, что в процессе диффузионного распада жидкости кристаллы эвтектических фаз могут зарождаться и расти порознь или совместно [4].

Изменение характера кристаллизации - переход от равноосной, хорошо ограненной формы роста к игольчатой - можно объяснить влиянием присутствия хотя бы и очень малочисленных, но всегда существующих примесей. Всякий кристаллик начинает расти в хорошо ограненной форме, но по мере роста вокруг него скапливаются оттесняемые на периферию инородные атомы и молекулы, которые в конце концов должны были бы создать пленку другого вещества, мешающего нормальному росту кристалла. Диффузия и конвекция не дают образоваться сплошному слою примесей и позволяют кристаллу расти. Однако при низких температурах, особенно при повышенном содержании примесей, влияние диффузии и конвекции недостаточно, чтобы обеспечить правильное питание кристалла по всем направлениям. При этом рост кристалла по нормалям к его граням затрудняется сильнее, чем тангенциальный рост граней путем перемещения ребер. Это объясняется большей омываемостью кристалла на ребрах, чем в центрах граней. Наибольшую скорость роста можно наблюдать там, где сильнее всего омывание кристалла жидкостью, т.е. по ребрам, и притом по ребрам меньших линейных размеров [3].

Следовательно, можно сделать вывод, что обычно наблюдаемые эвтектические структуры представляют собой частный случай дендритных структур.

Диаграмма состояния Al-Si (рис. 1) разработана давно, исследования последних лет касаются главным образом уточнения растворимости кремния в твердом алюминии.

Диаграмма состояния Al-Si относится к эвтектическому типу. Эвтектика образуется между алюминием и кремнием при температуре 577 °С и концентрации кремния 12,2 %. Максимальная растворимость кремния при эвтектической температуре составляет 1,65 %. По мнению Хансена и Андерко, наиболее точная растворимость кремния в алюминии получена в работе Дикса и Хиса, так как при исследовании был использован алюминий высокой чистоты и время отжига было достаточным для достижения равновесия. Растворимость кремния в алюминии, по их данным, следующая:

Температура, °С 577 550 500 450 400 350 300 200.

Растворимость, % 1,65 1,30 0,80 0,48 0,29 0,17 0,10 0,05.

Близкие данные по растворимости кремния в алюминии получила М.И. Захарова.

Некоторые уточнения температуры и состава эвтектики алюминия с кремнием приведены в работе Эллиота:

- эвтектическая температура составляет 577,2 °С, эвтектическая концентрация кремния 12,7 %;

- растворение кремния в алюминии уменьшает параметр решетки алюминия, каждый процент кремния уменьшает ее на 0,0017 %;

- сплавы Al-Si не упрочняются после термической обработки;

- сплавы Al-Si характеризуются хорошими литейными свойствами: высокой жидкотеку-честью, хорошей герметичностью, а также высокой коррозионной стойкостью.

Значительным преимуществом алюминиевых сплавов с кремнием является сравнительно низкий коэффициент линейного расширения, благодаря чему эти сплавы нашли широкое применение для деталей точных приборов.

Влияние способа литья на уровень свойств определяется главным образом скоростью кристаллизации и функционально связанной с ней дисперсностью микроструктуры. При литье в песчаные формы и ио выплавляемым моделям отливки охлаждаются и крссталлизуются медленнее, чем ири литье в кокиль и под давлением. В результате дендритные ячейки (Al) и кристаллы избыточных фаз в отливках, закристаллизованные в металлических формах, оказываются заметно больше разветвленными, чем кристаллы кремния и других фаз в составе эвте ктических колоний.

При старении после литья количество и морфология фаз кристаллизационного происхождения изме няются слабо. Происходит снятие литейных напряжений и распад алюминиевого

t.r

Si, T. / ne нйсс() S W te M W 50 (0 78 Sû sù tov

t3ùs

tm

sqd

7Щ

cf

5BU

jao

t't

-,set

SU!

4.___iT}\

{M}

/ wffii

H

в I

ЯХ {ami

Sel* il '

&

577-

_

ioad

sus

*

(Л)

X'tSi) pHWtti JTA

njt л

Si t V* fOJn) -1—-

wtt*

0 lù t§ $Ê 40 Se 60 ?0 SÛ SB 1BO Ai St, % (вяг) Si

Рис. 1. Диаграмма состояния Al-S i

твердого раствора в тех микрообъемах, которые пересыщены в результате неравновесной кристаллизации относительно предельной растворимости при температуре старения [5].

Типичная микроструктура силумина состоит из мягкой пластичной фазы - а-твердого раствора кремния в алюминии и эвтектики а+Si, состоящей из а-твердого раствора и иголок кремния.

Двойные сплавы не упрочняются термической обработкой. Это объясняется высокой скоростью распада твердого раствора, который частично происходит уже при закалке, а также большой склонностью ic коагуляции стабильных выделений кремния [6].

В качестве исходной заготовки для всех исследованных образцов использованы специально отлитые слитки.

Все образцы в исходном литом состоянии имели типичную структуру немодифицирован-ного алюминий-кремниевого сплава (рис. 2).

В структуре видна типично игольчатая эвтектика (Al-Si) и избыточные дендритные кристаллы алюминиевого твердого раствора, что характерно ник для доэвтектического сплава, так и для лвтектического, закристаллизованного в условиях сущнственного охлаждения (литье в металлическую форму).

Закалочный экспоримент продусматривал изу чение хлрактера изменения гетерогенного строения слабоперегретого расплава над точкой эвтектики в течение длительного времени - 10 ч - при температуре 580 °С, т.е. на 3 °С выше температуры эвтектического равновесия (577 °С). Немедленное очень быстрое охлаждение в холодной воде, предотвращая диффузионное перераспределение компонентов за тысячные доли секунды (скорость охлаждения ~104 °С/с), сохраняет в закаленном образце степень гетерогенности раствора, возникшего при нагреве и длительном выдерживании. Структура закаленного образца (рис. 3) отражает неоднородное распределение компонентов в таком растворе.

Следовательно, в жидком растворе в результате гетерогенизации возникают кластеры, обогащенные алюминием (светлыр округлые дендриты), а такжр рластерные участки, обогащенные кремнием (области на фотонрафии, оанятые высокодисперсной смесью темных и с ветлых кристаллов, располагающихся между кластерами, обогащвнными алюминием). Такой характер структуры закаленного образца является важным свиде тельством того, что области

Рис. 2. Микроструктура сплава АК12 (исходное Рис. 3. Микроструктура сплава АК12 в закаленном литое состояние) состоянии (Тз = 580 °С, т =10 ч, охлаждение в

воде)

- 58(5 -

в жидком растворе, обогащенном кремнием, не становятся сплошными кремниевыми кристаллами при сколь угодно длительной выдержке в субкристаллическом интервале температур выше точки эвтектического равновесия.

Поэтому все структуры, которые будут представлены, для различных условий длительности выдержки и охлаждения являются лишь отражением преобразования разнородных и однородных кластеров и взаимодействия при преобразовании их для формирования зародышей кристаллизации и роста кристаллов соответствующих фаз.

Конечная форма образующихся кристаллов зависит как от степени гомогенизации и ге-терогенизации образующегося эвтектического расплава над точкой эвтектики в начальные периоды выдержки при температурах слабого перегрева, так и от условий взаимодействия сформированных кластеров между собой при последующем охлаждении с разной скоростью, определяя рост кристаллов различной формы и размеров (от мелких ориентированных, разо-риентированных компактных до крупноигольчатых и крупных компактных кристаллов).

Использование минимальной степени перегрева над точкой эвтектического равновесия (температура под закалку (Тз) составляет 0,580 °С, т.е. на 3 °С выше 577 °С), а также кратковременной выдержки (5 мин), обеспечивая переход в новое агрегатное состояние исходного образца и формирование его кластерного строения, не позволяет еще реализоваться каким-либо изменениям в форме и распределении областей в расплаве, обогащенных тем или иным компонентом, возникших при плавлении в микрообъемах, обогащенных алюминием и кремнием.

Поэтому характер распределения компонентов в расплаве повторяет картину их расположения в фазах исходного образца, из которого этот расплав образовался при перегреве.

Относительно быстрое охлаждение (на воздухе) такого расплава приводит к формированию структуры, напоминающей исходную структуру образца до расплавления.

Из кластеров, обогащенных алюминием, вновь формируются овальные зерна твердого раствора на его основе. Из кластеров, обогащенных кремнием, вновь формируются кремниевые кристаллы, располагающиеся в структуре в виде ориентированных игольчатых образований, наследующих свою форму также от исходной структуры первичного образца (рис. 4).

Увеличение степени перегрева на 8 °С (585 °С) с той же длительностью выдержки в состоянии кластеризации (5 мин) приводит к более коренным преобразованиям в кластерной структуре только что полученного расплава.

Более высокая температура выдержки позволяет реализоваться частичному выравниванию кластерного строения за счет выравнивающей диффузии в жидкости и перераспределения кластеров между объемами, обогащенными алюминием и кремнием.

При охлаждении на воздухе более однородное кластерное строение расплава обеспечивает получение структур, слабо напоминающих исходное строение эвтектического сплава с более однородным распределением мелких кремниевых кристаллов, окруженных твердым раствором на базе алюминия (рис. 5).

Увеличение длительности выдержки кластеризованного сплава до 10 мин независимо от степени перегрева над точкой эвтектического равновесия (3 и 8 °С, т.е. при температуре 580 и 585 °С) приводит к усилению взаимодействия между кластерами, обогащенными кремнием, что, в свою очередь, приводит к образованию в структуре при охлаждении на воздухе компактных кремниевых кристаллов микронного размера, почти полностью свободных от грубых

Рис. 4. Микроструктура сплава АК12 в термообработаном состоянии (Тз = 580 °С, т = 5 мин, охлаждение на воздухе)

Рис. 5. Микроструктура сплава АК12 в термообработанном состоянии (Тз = 585 °С, т = 5 мин, охлаждение на воздухе)

Рис. 6. Микроструктура сплава АК12 в термообработанном состоянии (Тз = 580 °С, т =50 мин, охлаждение на воздухе)

■».J /V тЧЛ ; • /

Рис. 7. Микроструктура сплава АК1 2 в термообработанном состоянии (Тз = 585 °С, т= 10 мин, охлаждение на воздухе)

игольчатых кремниевых выделений, окруженных матрицей из пластичного твердого раствора на основе алюминия (рис. 6 - 7), все микроструктуры получены при увеличении х500.

Дальнейшее увеличение времени выдержки в условиях кластеризации до 15 мин не позволяет сохранить компактную форму при последующей кристаллизации с охлаждением на воздухе (рис. 8 - 9).

Наряду с сохранением некоторого числа компактных кристаллов в структуре наблюдается преимущественно игольчатая пластинчатая форма кристаллов кремния. Этому способствуют процессы самоорганизации в кластерной структуре расплава, которые происходят за счет стремления системы к минимуму свободной энергии, когда на уровне наноразмерных кластеров они ориентируются в жидкости относительно друг друга одноименными кристаллографическими плоскостями. Большие скорости охлаждения (охлаждение на воздухе) способствуют реализации только одного из семейства кристаллографических направлений, что приводит к образованию пластинчатых (иггльчатых на выявленных структурах) кристаллов.

Аналогичная форма кристаллов появляется при выдержке 15 мин (рис. 8) и при всех остальных исследованных, вплоть до 10 ч (рис. 10 - 13).

Причем эта закономерновть сохраняется при темпаратуре клвстеризации как 580°С, так и 585 °С (рис. 14 - 16П, увеличение х500.

Рис. 8. Микроструктура сплава АК12 в Рис;. 9. Микроструктура сплава АК12 в термообработанном состоянии (Тз = 58=0 °С, термооброботанном состоянии (Тз = 585 °С,

т =15 мин, охлаждение на воздухе)

т = 15 мин, охлаждение на воздухе)

1JF

Рис. 10. Микроттруктура сплава иАБС12 в Рис;. 11. М^ЮикЕсраострруукст^цра сплав- АК12 в термообработанном состоьнии (IT3 = 5 А0 °С, термообработанном состоянии (Тз = К580 °С,

т = 20 мин. охлажднние на воздухеб

т = 30 мин, охлнждениа на втндухе)

Рис. 12. Микроструктура нплава АК12 в Рис. 13. Микроструктура сплава АК12 в термо обрабо тайном состоянии (Тз = 5Н0 °СН, термообработанномсосооянии(Тз=580 °С, т= 1 0 ч,

т = В0 мин, охлаждение на воздухе)

охлаждение на воздухе)

Формироионие окончательной с;трэу1сн51р}ы споавов, прошедших стадию разделительной кластеризации в расплаве выыше точки эвтектического равновесии, опредебается не длительно -стью выедержки, а скоростью охлаждения кластеризованного расплава в интервале изменения агрегатного оостояния сплава. Поэтому быьло получено алияние споаобв охмаждения [в воде, на воздухе и 5 печью - в интерврле кристаллизации) на втруктрру сплава АК12.

-5- 589 -а

Рис. 14. Микроструктура сплава АК12 в термообработанном со стоянии (Тз = 55555 °С, т = 25 тмин, охлаждение на воздухе)

Рис;. 15. Микроструктура сплава АК12 в термообработанном состоянии (Тз = 585 °С, т = 350 мин, охлаждение на воздухе)

щт*

£ * [ 'VI

1 V- у Д'^Х'

Рис. 16, Микроструктура сплава АК12 в термообработанном состоянии (Тз = 585 °С, т = 90 мин, охл аждение на воздухе)

Роль этога параметра заключается в то м, что при использовании разных скоростей охлаждения процессы кристаллизации (взаимодлйствие с формированным кластеров между собой при образовании злролыша а рост кристаллов) вынуждтны проааоать за рааллчные промежутки времени.

Рас смттривая клдлтерныа взаимодействия кол некий особый механизм диффуаии, когда в перемещении масс вещества участвуют не отдельные атомы, а целые комплексы кристаллографически связанных между собой групп атомов, эти д ффузионные процессы могут проходить с большими скоростями и приводито к резко различным фирмам обязующихся кртсталлов в зависимости от внешних условий (температу ра, время, разница химических потенциалов, концентрацией! и т.д.)).

Испольаование слорости охлаждоная ~а04 аС/с обеспечивает выдкление кристаллографически нз ориентирссанноех керае;мниевых кристаллов (видимо, нонокристаллав) , размер которых можат быть саизмерим а размером кластеров в гетерогенном расплаве.

Увеличение длительности кристаллизационных процессов за счет применения охлаждения образцов сплава на воздухе (скорость охлаждения 1-3 град/с) приводит к формированию кристаллов микронных (несколько микронов) размеров. При этом форма кристаллов может быть различной - от компактных кристаллов до пластинчатых (игольчатых), в зависимости от

степени гетерогенизации расплава с применением времени выдержки в интервале температур кластеризации и гетерогенизации жидкости.

Применение технологической операции охлаждения кластеризованного расплава с печью (скорость охлаесдения 0,01-0,03 °С/с), т.е. на два порядка меньше, чем охлаждение на возду ие и на 55-6 порядков меньше, чем охлаждение в воде, должно было привести к тому, что формирование зародышей и их кристаллографически ориентированный рост будут происходить путем присоединения диффундирующих кластеров кремния не только в одном преимущественном кристаллографическом направлении, но и в других, менее благоприятных для роста. Поэтому в процессе роста будут формироваться компактные кристаллы кремния, достаточно крупные, имеющие криаталлографическую ограноу /рис. 17 - 1_9).

Нагрев эвтектических сплавов до температуры, на несколько градусов превышающей точку эвтектического раоновосия, вызывает расплавление, резулетатом охлиждения которого становится формировение кластерной смеси. Такая смесь состоит из микрообъемов, характеризующихся наличием ближнего порядпа, отвечающего криеталлическовпу строению фзз, при плавлении воторых эти кластеры образоволиаь.

При нсличии разницы в строеяии соответетвующих кластерных оброеований в жидкости происходит перераспределение компонентов, участвующих в образовании таких кластеров,

Рис. 17. Микроструктура сплава АК12 в термообработанном состоянии (Тз = 578 °С, т = 30 мин, ехлаждение с печью до 400 °С, в воду)

Рис. 18. Мбикроструктура сплава АК12в термообртботанном состоянии (Тз = 580 °С, т= 15 опин, охлаждение с печью до 450°С, в воду)

Рис. 19. МСикроструктура пплава АК12 в термообработанном состоянии (Тз = 585 °С, т = 120 мин, охлаждение с печью до 430°С, в воду)

которые мы фиксируем при резком охлаждении в воде или при более малых скоростях охлаждения (на воздухе, с печью).

Выводы

Таким образом, в настоящей работе обнаружено явление влияния предкристаллизацион-ных процессов, развивающихся в температурной области на несколько градусов выше начала температур фазового превращения, на формирование конечной структуры сплавов. Сущность обнаруженного явления - предкристаллизационная гетерогенизация исходного металлического раствора в связи с перераспределением компонентов сплава между кластерными образованиями, возникающими в таком растворе до развития фазового превращения. Это явление может найти разнообразные области применения при получении слитков и отливок из сплавов эвтектического типа.

Показано, что степени кластеризации и гетерогенизации расплава зависят от температуры и длительности выдержки в условиях перегрева выше точки эвтектического равновесия, что сказывается на структуре формирующегося при кристаллизации сплава. Установлены условия формирования структуры с компактными выделениями кремниевых кристаллов, которые могут быть рекомендованы при изготовлении изделий различного назначения.

Возможны три различные технологические схемы кластеризационной обработки эвтектических сплавов с различной реализацией.

Наиболее трудно реализуется схема - нагрев до температуры кластеризации расплава с различной длительностью выдержки, с очень быстрым закалочным охлаждением. Недостаток - трудность сохранения формы изделия.

Вторая схема - с использованием ограниченной длительности кластеризации 5-10 мин, с охлаждением на воздухе. Недостаток - возможная нестабильность результатов при изготовлении изделии, опасность сохранения пластинчатых (игольчатых) кристаллов кремния.

Третья схема - с использованием различных длительностей кластеризации, с охлаждением с печью. Недостаток - относительно крупные кремниевые кристаллы. Положительным фактором является возможность применения для сильно заэвтектических сплавов, где новые кристаллы кремния будут соизмеримыми с первичными, но сплавы будут характеризоваться повышенной пластичностью, обусловленной устранением непрерывности эвтектического кремниевого кристалла.

Список литературы

[1] Пригунова А.Г., Белов Н.А., Таран Ю.Н. Силумины. Атлас микроструктур и фракто-грамм промышленных сплавов. М.: МИСиС, 1996. 175 с.

[2] Белов Н.А., Савченко С.В., Хван А В. Фазовый состав и структура силуминов: справочное издание. М.: МИСиС, 2008. 283 с.

[3] Беляев А.И., Романова О.А., Бочвар О.С. и др. Алюминиевые сплавы. Металловедение алюминия и его сплавов: Справочное руководство. М.: Металлургия, 1971. 352 с.

[4] Таран Ю.Н., Мазур В.И. Структура эвтектических сплавов. М.: Металлургия, 1978. 312 с.

[5] Золоторевский В.С., Белов Н.А. Металловедение литейных алюминиевых сплавов. М.: МИСиС, 2005. 376 с.

[6] Грачев С.В., Бараз В.Р., Богатов А.А., Швейкин В.П. Физическое металловедение: учеб. для вузов. Екатеринбург: УГТУ-УПИ, 2001. 534 с.

[7] Аникина В.И., Жереб В.П., Аникин А.И. и др. Пат. 2525872 РФ (2013) / Б. И. 2014. № 23

[8] Биронт В.С., Блохин И.В., Аникин А.И. // Современные технологии ресурсов. Вып. 8.