CC BY
191
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, ЧИСЛЕННЫЕ МЕТОДЫ И КОМПЛЕКСЫ ПРОГРАММ
УДК 539.3:004.9
И. В. Александров, Р. Г. Чембарисова, М. И. Латыпов
АНАЛИТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НАПРЯЖЕНИЯ ТЕЧЕНИЯ СПЛАВА AL 1570, ПОДВЕРГНУТОГО ИНТЕНСИВНОЙ ПЛАСТИЧЕСКОЙ ДЕФОРМАЦИИ
В работе представлена развитая на основе уравнений дислокационной кинетики математическая модель, позволяющая описывать и анализировать прочностные характеристики металлов и сплавов. Составлены уравнения эволюции плотности дислокаций в сплаве А1 1570 с учетом его состояния и структурных особенностей, сформированных в результате интенсивной пластической деформации. Рассчитаны зависимости условного предела текучести от размера зерен в координатах Холла-Петча, а также интенсивности процессов динамического возврата от размера зерен. Ультрамелкозернистые материалы; закон Холла-Петча; моделирование деформационного поведения
ВВЕДЕНИЕ
Механические свойства являются важным аспектом привлекательности объемных ультра-мелкозернистых (УМЗ) материалов как нового поколения конструкционных материалов. Изначально надежды исследователей на высокую прочность УМЗ материалов связаны с широко известным законом Холла-Петча (ХП) [1], который предсказывает повышение прочности материала при уменьшении размера зерен:
т7 = тохп + Кй-1/2, (1)
где тТ - предел текучести или напряжение течения, т0 хп - сопротивление движению дислокаций, не зависящее от размера зерен, К - коэффициент ХП, й - средний размер зерен.
Несмотря на то, что первоначально это соотношение устанавливало связь исключительно между нижним пределом текучести и размером зерен в стали, в литературе указывают на справедливость выражения (1) для широкого круга металлических материалов [2].
Более того, согласно работе [3] выражение
(1) выполняется и для напряжения течения, значительно превышающего предел текучести.
1. ОТКЛОНЕНИЯ ОТ ЗАКОНА ХОЛЛА-ПЕТЧА
Начиная с 1989 г., в литературе активно обсуждаются отклонения от значений, предсказанных соотношением ХП [2, 4-10]. В работе
[2] в зависимости от того, насколько применим закон ХП для различных материалов, в результате анализа экспериментальных данных выделены три области размеров зерен:
а) область от монокристалла до зерен размером порядка 1 мкм, где соотношение (1) выполняется;
б) область зерен размером от 1 мкм до 30 нм, где закон ХП грубо выполняется, но классический показатель -0,5 отклоняется к значению, близкому к нулю;
в) область зерен с очень мелким размером менее некоторого критического (порядка 20.. .30 нм), для которой характерен нулевой показатель с отсутствием повышения или даже понижением прочности при уменьшении размера зерен.
Тип (в) уменьшения прочностных характеристик с уменьшением размера зерен называют «обратным» соотношением ХП [5, 7, 11].
2. МЕХАНИЗМЫ И МОДЕЛИ
На сегодняшний день основным механизмом деформации материалов с размером зерен, соответствующим группам (а) и (б), общепринято считать дислокационный механизм [5]. При этом пластическая деформация поликристалла отличается от таковой, характерной для монокристалла, по двум причинам [12]:
• барьерный эффект - непроницаемость границ зерен (ГЗ) для дислокаций;
• усложнение деформации (для сохранения сплошности материала в ходе деформации при разной ориентировке плоскостей скольжения в каждом зерне должно действовать больше систем скольжения, чем в монокристалле той же ориентировки).
В рамках дислокационного механизма для объяснения возникновения соотношения типа (1) были разработаны две классические модели
[3].
Контактная информация: 8(347)273-79-77
Модель скопления дислокаций. Так как ГЗ являются непроницаемыми барьерами для дислокаций, скольжение приводит к скоплению дислокаций у ГЗ, что вызывает появление концентрации напряжений и активизацию дислокационных источников в соседних зернах (эстафетный механизм деформации).
Модель деформационного упрочнения.
В этой модели предполагается, что при данной степени деформации плотность дислокаций в образце с мелкими зернами выше, чем с крупными, что приводит к большему внутреннему напряжению и, следовательно, к большему напряжению течения у такого материала. Следовательно, величина зерен влияет на напряжение течения материала косвенно благодаря воздействию масштабного фактора на плотность дислокаций.
Относительно механизмов деформации материалов с наиболее мелкими зернами (менее 30 нм) мнения исследователей расходятся. В литературе имеются упоминания о действии таких механизмов, как скольжение дислокаций, скольжение частичных дислокации, двойнико-вание, зернограничное проскальзывание, поворот зерен [13].
Несмотря на многообразие точек зрения, авторы приходят к единому мнению о возрастающей роли ГЗ при деформации УМЗ материалов [13, 5].
3. АЛЮМИНИЕВЫЕ СПЛАВЫ, ПОДВЕРГНУТЫЕ ИПД
Одним из привлекательных объектов исследований на сегодняшний день являются алюминиевые сплавы, в которых УМЗ состояния формируются в результате интенсивной пластической деформации (ИПД). Показатели прочности этих сплавов также обнаруживают отклонения от закона ХП, однако, в отличие от описанных ранее, - в большую сторону [9].
Авторы работ [9, 14] получили высокие значения предела текучести около 900 МПа в алюминиевых сплавах А1 1570 (Л1-5,7М§-0,328е-0,4Ми, вес. %) и А1 7475 (А1-5,72п-2,2Mg-1,6Cu-0,25Cr, вес. %), подвергнутых интенсивной пластической деформации кручением (ИПДК) при комнатной температуре (КТ).
Особенности структуры и механические свойства сплава А1 1570, подвергнутого ИПДК, подробно описаны в работе [14]. Образцы сплава, вырезанные из горячепрессованного прутка, нагревали до температуры 380 °С в течение 2 ч с целью получения предельной концентрации
магния в твердом растворе, достижения максимального эффекта твердорастворного упрочнения, а также обеспечения максимального измельчения зерен в процессе деформации. После ИПДК исходная крупнозернистая (КЗ) структура полностью трансформируется в УМЗ структуру со средним размером зерен 130 нм и преимущественно большеугловыми границами зерен. Были также получены свидетельства высокого уровня внутренних упругих напряжений в этом сплаве.
Выдающуюся прочность этих сплавов связывают, прежде всего, с выделением атомов магния из твердого раствора в области ГЗ в результате ИПДК. К выводу об образовании сегрегаций в области ГЗ авторы [9, 14] приходят, опираясь на полученные данные об изменении параметра решетки, и подтверждают этот факт результатами исследований методом пространственной атомной томографии.
Выделение растворенных атомов из твердого раствора будет приводить к уменьшению параметра решетки, в частности, для системы А1-Mg - на 0,0043 А или на 0,0046 А согласно [12] и [15] соответственно. Выделение атомов легирующих элементов из твердого раствора может быть связано с образованием сегрегаций или выделением частиц. В случае сплава А1 1570 полученное значение концентрации атомов магния в области шириной около 3 нм вдоль ГЗ, равное » 1,8 X 1028 м-3 (» 23,1 ат. %), намного ниже, чем характерное для частиц A13Mg2 значение (40 ат. %), что исключает второй вариант.
Похожие данные об образовании сегрегации атомов легирующих элементов были получены при исследовании сплава А1 6061 системы А1-Mg-Si (А1-1,0Mg-0,6Si-0,3Cu-0,25Cr-0,15Мп-
0,7Ре-0,252п-0,15Ті, вес. %) в работе [16]. Сплав нагревали до температуры 530 °С в течение 2 ч с последующим быстрым охлаждением в воде со льдом. В результате в сплаве была сформирована КЗ структура со средним размером зерен около 500 мкм. Механические испытания на растяжение показали существенное увеличение предела текучести с одновременным уменьшением пластичности до 5,5 %. Старение образцов сплава после ИПДК не способствовало дальнейшему увеличению предела текучести. Значение предела текучести (660 МПа) было намного выше, чем после традиционной обработки (~ 270 МПа) - искусственного старения в течение 12 ч при 160 °С после закалки с температуры 530 °С. Оно было также выше, чем предел текучести этого же сплава, обработанно-
го другими видами ИПД, такими как равнока-нально-угловое прессование (РКУП) (390 МПа
[17]), РКУП с последующим старением (415 МПа [17]), РКУП с последующей холодной прокаткой (475 МПа), многоцикловая прокатка (363 МПа) или РКУП, использованное для консолидирования порошка [16].
Изучение структуры сплава после ИПДК показало, что равноосные зерна со средним размером около 100 нм с высокоугловыми ГЗ формируются как в центре, так и на середине радиуса образцов. Средний размер зерен около 100 нм наблюдался также в образцах алюминиевых сплавов систем A1-Mg-Li и A1-Zn-Mg, в то время как в технически чистом алюминии средний размер зерен достигал лишь порядка 800 нм, что может говорить о том, что наличие легирующих элементов играет ключевую роль в измельчении зеренной структуры [16]. К другим особенностям структуры сплава А1 6061, подвергнутого ИПДК, можно отнести высокий уровень внутренних упругих напряжений и искажений кристаллической решетки, связанный с наличием высокой плотности дислокаций, а также наличие дисперсных частиц Mg2Si со средним размером 10...20 нм. Исследования методом 3Б атомной томографии показали уменьшение концентрации атомов легирующих элементов в твердом растворе в результате ИПДК аналогично случаю атомов магния в сплаве А1 1570, но в менее выраженной степени. Авторы связывают это уменьшение также с образованием сегрегаций, удлиненная форма которых может свидетельствовать об их преимущественном зарождении вдоль дислокаций.
Таким образом, при исследовании структуры алюминиевых сплавов, обнаруживающих весьма (сплав А1 6061) или даже рекордно высокую (сплав А1 1570) прочность, выявлены сегрегации примесных атомов вдоль ГЗ. Авторы работ [9, 18] закономерно отводят этому явлению ключевую роль в формировании таких сверхпрочных состояний.
4. ДИСЛОКАЦИОННЫЕ МОДЕЛИ
Одним из наиболее универсальных подходов к изучению деформационного поведения материалов является моделирование на основе уравнений дислокационной кинетики, которые позволяют получать зависимости средней плотности дислокаций и напряжения течения от степени деформации, размера зерен, температуры, и скорости деформации. Основное преимущество таких моделей состоит в том, что они спо-
собны описать большинство наблюдаемых эффектов, касающихся прочности и пластичности металлов в диапазоне от наноструктурных (НС) до КЗ состояний [11].
Как следует из литературных данных, УМЗ и НС материалы являются интересными и перспективными объектами исследований как с исследовательской (с целью изучения физических механизмов на субмикронном и нано уровнях), так и с практической точек зрения. УМЗ и НС материалы претендуют на роль материалов нового поколения с уникальным сочетанием свойств. В то же время, остается множество нерешенных проблем, в том числе касающихся выяснения причин формирования механических свойств этих материалов. Среди них — причины отклонения уровня прочностных свойств от предсказываемого законом ХП в диапазоне от ультрамелких до нанометрических зерен. Целью настоящей работы является разработка и применение дислокационной модели для анализа экспериментальных данных, касающихся прочности сплава Al 1570, подвергнутого ИПД, обнаруживающей отклонения от закона ХП.
ОПИСАНИЕ МОДЕЛИ
Учитывая представленные выше результаты экспериментальных исследований в развиваемой модели для анализа напряжения течения алюминиевых сплавов, подвергнутых ИПД, были учтены три возможных механизма деформационного упрочнения.
1. Закрепление дислокаций атмосферами примесных атомов
Поскольку, как показал эксперимент, ИПД может приводить к выделению сегрегаций в сплаве Al 1570, в развиваемой модели проведен учет закрепления дислокаций примесными атомами.
Скопление атомов примеси на дислокациях и дефектах упаковки (сегрегация) ограничивает возможности пластического течения. Отягощенная сегрегацией дислокация требует дополнительных усилий для своего перемещения oag, а если такое напряжение sag недостижимо, то это означает блокировку - исключение ее из числа подвижных. В обоих случаях сегрегации увеличивают сопротивление пластической деформации и упрочняют материал [12].
Напряжение закрепления tag определяется временным законом изменения добавочного напряжения, связанного с накоплением приме-
сей в зависимости от времени в окрестности краевой дислокации в алюминиево-магниевых сплавах [19]:
т „ = Ат 0(1 - 7''),
Г, = (6оо8Ь(АЖ/2кТ) Гс)-1,
Гс =Vпе-
(2)
Ат 0 = 2ас0 АЖ гапЬ(АЖ /2кТ )/7эЬ3
где Ат0 — асимптотическое значение напряжения закрепления, АЖ » 0,13 эВ — разность средней энергий связи между областями сжатия и растяжения в окрестности ядра дислокации, АНС ~ ~ 0,97 эВ - среднее значение энтальпии активации перехода из области растяжения в область сжатия, а' ~ 0,56 — параметр, зависящий от термодинамической движущей силы АЖ ~ 0,13 ±
0,02 эВ и величины АЖ, с0 — атомная концентрация раствора.
2. Торможение дислокаций при консервативном скольжении
На начальном этапе пластической деформации ИПД материалов с сегрегациями примесных атомов вдоль ГЗ дислокации могут оторваться от атмосфер примесных атомов и скользить в теле зерен. При этом дислокации при своем движении взаимодействуют друг с другом, а также с другими дефектами и препятствиями, которые затрудняют движение дислокаций. Для дальнейшего движения им требуется дополнительный подвод энергии, чему и соответствует упрочнение. В рамках развиваемой модели учтены как торможение дислокаций при их упругом взаимодействии, так и при их взаимодействии с твердым раствором.
Торможение дислокаций при их упругом взаимодействии. Между дислокациями существует упругое дальнодействие через поле напряжений [20]. По мере развития пластической деформации по разным причинам торможение дислокаций при их взаимодействии усиливается, что обусловливает наклеп. Экспериментальные данные и теоретический анализ с использованием разных моделей торможения дислокаций показал, что напряжение течения та растет пропорционально корню квадратному из плотности дислокаций [21] согласно соотношению
т а = а^р1'
(3)
где а - коэффициент, учитывающий термическую активацию дислокационного взаимодействия, G - модуль сдвига, Ь - вектор Бюргерса.
Упрочнение твердых растворов. В общем виде прирост предела текучести, благодаря твердорастворному упрочнению имеет вид
Ос = КС,
где КС - концентрационное упрочнение, с - весовая концентрация растворенного элемента, п -показатель степени.
В различных моделях п варьируется от 1/2 до 4/3 [12, 22, 23, 24]. Однако обычно прирост не настолько велик, чтобы отличить его от линейной зависимости с п = 1 [12].
Трение решетки. Для того чтобы продвинуть единственную прямолинейную дислокацию в совершенном кристалле, требуется приложить к нему некоторое конечное напряжение. Сила, препятствующая скольжению прямолинейной дислокации, называется силой Пайер-лса-Набарро [21, 22]. Вычислить напряжение Пайерлса-Набарро очень трудно, а расчет ее величины, проведенный в литературе для конкретных структур [22, 25], является недостаточно строгим.
3. Упрочнение за счет эволюции плотности дислокаций
В развиваемой модели использована известная связь между напряжением течения та, обусловленным взаимодействием дислокаций друг с другом, и плотностью дислокаций р (3) [3, 26].
Считалось, что увеличение плотности дислокаций в поликристалле (из-за наличия в нем ГЗ) происходит в результате ограничения свободного пробега дислокаций и несовместности деформации в различных зернах. Согласно работе [11], для учета влияния ГЗ на накопление дислокаций введем в уравнение эволюции плотности дислокаций слагаемое
ар+ _р_
ау ьа
с обобщенным коэффициентом в, который будет подгоночным.
Помимо влияния ГЗ дополнительно учитывалось возможное увеличение плотности дислокаций вследствие их размножения двойным поперечным скольжением. Этот процесс рассмотрен в работах [12, 25]. Различают размножение на препятствиях деформационного и недеформационного происхождения, в обоих случаях накопление дислокаций определяется длиной свободного пробега дислокаций между препятствиями - Ха и ХР соответственно [11]:
1 1
ау ЬХа ЬХр
Для расчета ХР воспользуемся эмпирической зависимостью, представленной в работе [11],
ОЬ
(4)
X р А-р
Тр :
где Ар - численный коэффициент, тР - упрочнение частицами выделений.
Известно, что сидячие дислокации являются эффективными барьерами для скользящих дислокаций (барьеры Ломера-Коттрелла) [21]. Поэтому их также учитываем как препятствия, на которых возможно размножение двойным поперечным скольжением. Расстояние между актами последнего Ха определяли согласно выражению, полученному на основе экспериментальных данных [11]:
Ха = Аа 4т, (5)
Р
где Аа - численный коэффициент.
В то же время предполагали, что плотность дислокаций может уменьшаться во время деформации в результате динамического возврата [12].
В рамках данной модели считалось, что указанные процессы происходят независимо. Это приводит к их линейному сложению в уравнении эволюции плотности дислокаций, как и в работе [26]:
ёр в 1 1 ,
- = — +---------------+-----------------кр .
(6)
ау ьа ьх
где к - коэффициент, учитывающий всевозможные процессы динамического возврата.
Аннигиляция краевых дислокаций в ГЗ определяется временем их растворения в границе ^ = 13 / \ЬДЬ [27, 11] и временем парной аннигиляции в них дислокаций противоположного знака t2 = I 2/4п2Д, [6, 11, 28, 29, 30],
П » » ОЬъ / кТ . Предполагая, что указанные
процессы протекают независимо друг от друга, а их характерные времена зависят от максимального расстояния в границе I = а, коэффициент к в уравнении (6) был записан в виде к = = к + к2, где
1
к = -
(7)
У
В выражениях (7) присутствует коэффициент самодиффузии в ГЗ Оь:
Вь = А0 ехР
кТ
(8)
где ДЬ0 - предэкспоненциальный множитель зернограничной диффузии, Qь - энергия активации самодиффузии в ГЗ.
Коэффициент аннигиляции винтовых участков дислокаций ка, согласно [31], равен
К =
5/2 (о л
п2 р
, Т , V а у
(9)
где - доля винтовых участков дислокационных петель, р - число систем скольжения; та -напряжение, определяющее критическое расстояние спонтанной аннигиляции винтовых дислокаций разного знака.
Напряжение та определяется суммой [31]
Та = ТС + ТIII ,
(10)
где тс - твердорастворное упрочнение.
Наличие слагаемого тс в выражении (10) обусловлено тем, что торможение дислокаций атомами растворенных элементов препятствует поперечному скольжению, так как примесь искривляет дислокацию, а для поперечного скольжения нужна абсолютно прямая дислокация [12].
Напряжение тш зависит от температуры и ЭДУ [12, 31]. Необходимое для расчетов значение было получено из соотношения [31]:
Г • Л
1п ^ = 1п-
о
,(0) кв
0(0) А(0)
1п
У
Т
(11)
где 0(0) - модуль сдвига при Т = 0, А(0) - константа, имеющая размерность энергии,
о (Т )Ь3(Т)
- безразмерная функция тем-
0(0)Ь3(0)
пературы для исключения влияния температурной зависимости модуля сдвига О и вектора Бюргерса Ь.
Пренебрегая небольшой температурной зависимостью этих величин (т. е. приняв F(T) ~ 1) и потенциируя выражение (11), было получено
Т і,,(0)
( • Л
о(0)
V У 0 У
А(0)
Используем зависимость тш(0) от ЭДУ у0 [31], представленную выражением
20(0) |
* ііі(0) = -
0,056 -
Ъ.
оь3
+
Т
квТ
У
п
где п - число дислокаций в скоплении. Формула для расчета критического напряжения тш, необходимого для поперечного скольжения, в зависимости от ЭДУ уД имеет вид:
20
0,056-
У Д оь3
кТ
< • 1 А(0)
1
V У 0 У
(12)
где отношение
кТ
определяется из графика
А(0)
[31], а у0 - из выражения [32]
У0 = Ь'ртVв ,
где рт - плотность подвижных дислокаций, VD -частота Дебая, определенная из выражения [22]
кВТД
к
(13)
где ТД - температура Дебая, к - постоянная Планка.
ЭДУ имеет «электронное» происхождение [20] и сильно зависит от состава твердых растворов, а в чистых непереходных металлах - от их валентности: высокая в поливалентных (А1, Zn, Mg), низкая - в одновалентных (Си, Аи, Ag). Растворение любой примеси замещения в любом чистом металле понижает ЭДУ [20]. В случаях ограниченной растворимости концентрационная зависимость ЭДУ чаще всего бывает
аУ Д
нелинейной, затухающей до
асл
0 у границы
растворимости. Ее зависимостью [20]
удачно аппроксимировали
У,
ХСа) = (0) =
ехр
Кдх (1 + х)2
где уД(с) - ЭДУ твердого раствора с концентрацией растворенного элемента сА, уД(0) - ЭДУ чистого металла, КД - некоторый коэффициент,
СА
х = —та^ - расстояние до границы растворимо-
С А
сти стах. Также в литературе встречается аппроксимация простой гиперболической зависимостью. В работе [33] подробно исследовано влияние концентрации основных легирующих элементов - Mg, Си, Ag конкретно на ЭДУ алюминия, поэтому при расчетах воспользуемся результатами этой работы.
С учетом аннигиляции винтовых участков дислокаций ка, коэффициент, учитывающий всевозможные процессы динамического возвра-
та к в уравнении (6), с учетом независимости всевозможных процессов динамического возврата принимает вид
к = к1 + к2 + ка .
СУПЕРПОЗИЦИЯ МЕХАНИЗМОВ УПРОЧНЕНИЯ
Наиболее общим и простым правилом сложения составляющих упрочнения является линейное сложение [12, 23, 34]:
*=Е , (14)
где т - напряжение течения или предел текучести, т - вклад /-го механизма упрочнения.
Линейная сумма дает лучшую аппроксимацию при наличии в материале небольшого количества прочных и большого числа мягких препятствий.
В рамках данной модели линейное сложение (14) можно считать удовлетворительным приближением, так как, например, при движении дислокаций и их взаимодействии на всю их длину действует сопротивление решетки и твердого раствора.
При определении эквивалентного одноосного напряжения течения поликристалла о [23]:
0= Их, (15)
где М - константа, называемая фактором Тейлора, использованы модель Тейлора [12, 23] и значение М = 3,07 [11]. Значения величин, использованных при расчетах, приведены в таблице.
РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ
Сплав А1 1570 является термически неуп-рочняемым, следовательно, размножение дислокаций в результате двойного поперечного скольжения у препятствий недеформационного происхождения осуществляться не будет.
Плотность дислокаций в сплаве А1 1570, подвергнутом ИПДК, достигает порядка 1014 м-2 [14]. Согласно выражению (5), длина свободного пробега дислокаций между сидячими дислокациями 1а ~10-5 м, что на два порядка превышает средний размер зерен а УМЗ материалов (порядка 10-7 м). Следовательно, этот механизм размножения дислокаций не будет вносить вклад в увеличение плотности дислокаций в УМЗ материале. Оценка величины 1а для КЗ состояния сплава А1 1570 дает тот же результат, но с меньшей разницей: на один порядок.
п
Т аблица
Значения величин, использованных при расчетах
Величина Значение Место Источник
использования
Константы материала
Тпл, К 933 [20]
То, К 0,46 Тпл (12) [20]
О, Па 26,8х109 (3), (12) [20]
Ь, м ал/2/2 (3), (6), (12) [21]
Значения, принятые в рамках данной работы
М 3,07 (15) [24]
°ьо, м2/с 8,7х10-5 (8) [36]
Яь, Дж 10,8 кВТпл. (8) [36]
п 10 (12) [12]
квТ/А (0) 0,006 (12) [37]
Ар 75 (4) [11]
АС 100 (5) [11]
Ю* 0,5 (9) [11]
Р 12 (9) [11]
Экспериментальные данные
КЗ сплав УМЗ сплав
С, м 60х10-6 130х10-9 (6) [14]
_2 р0, м 1011 0,51х1015 (17) [14]
00,2, Па 231х106 905х106 (19) [14]
То, Дж 90х10-3 (12) [33]
Условия испытаний
Т, К 298 [14]
• -1 1х10-4 [14]
8 , С
Подгоночные коэффициенты
а 0,54 (3)
в 1 (17)
В случае КЗ сплава считаем, что ГЗ находятся в равновесном состоянии, т. е. процессы уменьшения плотности дислокаций будут представлены только аннигиляцией винтовых дислокаций при поперечном скольжении: к = ка.
Характеристическое время «возврата» ГЗ быстро возрастает при увеличении размера зерен: ^ ~ А3 и для КЗ материалов принимает огромные значения, а соответствующий коэффициент к ~ 1 / ^ - ничтожно малые. Т. е. время ^ имеет эффективные значения и влияет на эволюцию плотности дислокаций только в материалах с очень мелкими зернами. С другой стороны, время ^ только для таких материалов и имеет физический смысл. Таким образом, несмотря на то, что в КЗ материалах ГЗ находятся в равновесном состоянии, в уравнение эволюции плотности дислокаций (6) можно внести коэффициент, учитывающий процессы «возврата» ГЗ для обеспечения единой формы этого
уравнения во всем диапазоне размеров зерен и построения в конечном счете зависимости вида ХП. С учетом вышесказанного уравнение эволюции плотности дислокаций для сплава А1 1570 во всем диапазоне размеров зерен принимает вид
- кр .
Ьс1
(16)
Интегрирование уравнения (16) позволяет получить выражение для эволюции плотности дислокаций в зависимости от степени деформации:
Р = Р 0ехр(- кУ )+-Л-(1 - ехр(- кг)), (17)
оак
где р0 - начальная плотность дислокаций в материале, т. е. плотность, имеющаяся в материале к началу испытаний на растяжение.
В общем случае условный предел текучести
о0,2 получается в результате действия
(3), упрочнения твердым раствором оС, сопротивления решетки о0, упрочнения дисперсными частицами т/. Под дисперсными частицами будем подразумевать сегрегации атомов магния вдоль ГЗ
а02 = а0 + ас + а/ + М ОЪ . (18)
Из-за недостаточной информации об упрочнении дисперсными частицами и твердым раствором, а также напряжении а0 в выражении (18) было принято
0 п = 0 0 + 0 С .
Значение напряжения оп оценивалось из экспериментальной зависимости предела текучести сплава в КЗ состоянии и состоянии после РКУП от значения размера зерна в степени -1/2 как предельное значение при аГ1/2®0 (рис. 1), оп » 200 МПа.
Величину о/ получали, вычитая другие составляющие упрочнения из экспериментально полученного предела текучести о0,2 эксп
а / а 0,2 эксп.
-ап - МаОЬр1
(19)
КЗ состояние. В случае КЗ состояния сплава А1 1570, полученного с помощью нагрева до 380 °С [14], когда атомы легирующих элементов
равномерно распределены в твердом растворе, а
11 —2
плотность дислокаций составляет 10 м , считаем, что
а 0,2 =а п+а + Ма°Ьр
(20)
УМЗ состояние. В этом случае помимо действия тех же упрочняющих механизмов, что и в КЗ состоянии сплава, учитывалось напряжение о/, а также дополнительное упрочнение атомами магния, концентрация которых повышена вдоль ГЗ (упрочнение в результате закрепления дислокаций примесными атомами т^). В результате
+ а + а г + МаОЬр1
(21)
При концентрации магния 23,1 (ат. %) (21,3 вес. %) в области ГЗ [14, 18] уравнение (2) предсказывает напряжение закрепления т^ -
- 121 МПа. В исходном отожженном состоянии (с0 = 5,7 вес. %, т^ - Дт0 = 32,4 МПа.)
Значения экспериментального и рассчитанного пределов текучести в зависимости от ат представлены на рис. 1.
Рис. 1. Зависимость условного предела текучести о0,2 от оТ1/2: ▲—
экспериментальные точки, 0 —
соответствующие модельные точки,
рассчитанные согласно уравнениям (17), (20). 1 — экспериментальные данные о пределах текучести алюминиевых сплавов А1 1560 и А1 5083 без дополнительного упрочнения атомами магния [35]. Нижняя экспериментальная точка при а = 130 нм соответствует состоянию после ИПДК при 100° С, ♦ — соответствующая модельная точка
Составляющие предела текучести для КЗ и УМЗ состояний сплава А1 1570 представлены на рис. 2, из которого следует, что высокая прочность УМЗ сплава А1 1570 во многом обусловлена упрочнением сегрегациями атомов магния вдоль ГЗ, а также деформационным упрочнением благодаря повышенной плотности дислокаций, накопленной в процессе ИПДК.
Рис. 2. Зависимость вклада в условный предел текучести сплава А1 1570 дислокационного упрочнения (♦) и
упрочнения за счет взаимодействия дислокаций с растворенными атомами магния 0^(0) от Ст
Варьирование размера зерен й в составленных уравнениях (17), (20), (21) позволяет проанализировать напряжение течения при степени истинной деформации гт = 0,2 % в координатах ХП. На рис. 3 видно, что оба состояния сплава А1 1570 попадают в область размеров зерен, где выполняется классический закон ХП, т.е. предел текучести возрастает при измельчении зерен. Ниже критического размера зерен йкр = 46 нм, получено «обратное» соотношение ХП: предел текучести падает при уменьшении размера зерен.
Рис. 3. График зависимости предела
текучести от размера зерен: 1 - данные о пределах текучести алюминиевых сплавов А1 1560 и А1 5083; 2 -теоретический предел текучести УМЗ сплава А1 1570 при комнатной температуре без учета дополнительного упрочнения атомами магния вдоль ГЗ
За такой вид кривой ответственно дислокационное упрочнение, пропорциональное квадратному корню из средней плотности дислокаций. Остальные составляющие предела текучести сдвигают вверх всю кривую, так как не зависят от размера зерен. Для контроля наклона кривой на график нанесены данные о пределах текучести сплавов А1 1560 и А1 5083 с размерами зерен около 400 нм, полученные в работе [35].
Ниспадающая часть кривой связана с ускорением процессов релаксации дислокаций в ГЗ
[18]. В целом, вышесказанное согласуется с распространенными в литературе представлениями о причине возникновения «обратного» закона ХП - возрастающей роли ГЗ в деформационном поведении материала.
Отклонение от закона ХП сплава А1 1570 в УМЗ состоянии, как и предложено в работах
[14, 9] и показано в работе [18], в значительной мере обусловлено влиянием сегрегаций атомов магния на процесс пластической деформации, концентрация которых в результате ИПДК повышена в области ГЗ, в свою очередь выполняющих основную роль «стоков» для дислокаций. Между точками, соответствующими условным пределам текучести сплава в КЗ и УМЗ состояниях, можно провести прямую, которая будет соответствовать классическому закону ХП, но с другими параметрами о0 хп и К (рис. 3). Несмотря на высокую прочность сплава А1 1570 с размером зерен ё = 130 нм, согласно полученной модели, предел текучести может быть еще выше при достижении «зерна максимальной прочности», размер которого ёкр = = 46 нм.
Суммируя вышесказанное, можно заключить: говоря о сплавах, подвергнутых ИПД, упоминать закон ХП не совсем уместно. В таких сплавах, «накачанных» упрочняющими факторами, действуют специфические механизмы упрочнения, не характерные для КЗ материалов, в то время как закон ХП установлен для материалов, прочность которых определяется в основном ограничением ГЗ свободного пробега дислокаций.
ВЫВОДЫ
1. Развита математическая модель, позволяющая анализировать деформационное поведение сплавов с учетом различных деформационных механизмов.
2. Составленные уравнения эволюции плотности дислокаций в рамках модели учитывают возможности:
а) действия механизмов увеличения плотности дислокаций в результате ограничения свободного пробега дислокаций ГЗ и несовместности деформации, двойного поперечного скольжения;
б) протекания процессов, ведущих к уменьшению плотности дислокаций в результате аннигиляции винтовых дислокаций при поперечном скольжении и релаксации неравновесных ГЗ.
3. Проанализированы возможные причины формирования высокопрочных состояний в УМЗ промышленном алюминиевом сплаве А1 1570, подвергнутом ИПДК при КТ. При этом учтены высокие значения ЭДУ алюминия и ее уменьшение из-за наличия магния в сплаве, что
оказывает существенное влияние на возможность протекания поперечного скольжения.
4. Рассчитаны составляющие условного предела текучести сплава Al 1570 в КЗ и УМЗ состояниях с учетом их микроструктурных особенностей. Установлено, что решающую роль в высокой прочности УМЗ сплава играет дополнительное упрочнение благодаря высокой концентрации магния вдоль ГЗ и увеличенное значение дислокационного упрочнения вследствие накопления высокой плотности дислокаций. Обе причины обусловлены применением ИПДК. Показано, что в УМЗ сплаве действует ряд упрочняющих факторов, не характерных для КЗ состояния, что приводит к невозможности выполнения закона ХП для этих материалов.
СПИСОК ЛИТЕРАТУРЫ
1. Hall E. O. The deformation and ageing of mild steel: iii discussion of results // Proceedings of the Physical Society. Section B. 1951. V. 64. № 9. P. 749-752.
2. Masumura R. A., Hazzledine P. M., Pande C. S. Yield stress of fine grained matterials // Acta Materia-lia. 1998. V. 46, № 13. P. 4527-4534.
3. Конрад Х. Модель деформационного упрочнения для объяснения величины зерна на напряжение течения металлов // Сверхмелкое зерно в металлах. М.: Металлургия, 1973. С. 207-219.
4. Meyers M. A., Mishra A., Benson D. J. Mechanical properties of nanocrystalline materials // Progress in Materials Science. 2006. V. 51. P. 427-556.
5. Wolf D., Yamakov V. MD Simulation of deformation mechanisms in nanocrystalline // Bulk Nanostructured Materials. Weinheim: WILEY-VCH. 2009. P. 165-199.
6. Малыгин Г. А. Нарушение закона Холла-Петча в микро- и нанокристаллических материалах // Физика металлов и металловедение. 1995. Т. 37, № 8. C. 2281-2292.
7. Kim H. S., Estrin Y., Bush M. B. Plastic deformation behaviour of fine-grained materials // Acta Materialia. 2000. V. 48. P. 493-504.
8. Zehetbauer M. J., Estrin Y. Modeling of strength and strain hardening of bulk nanostructured materials // Bulk Nanostructured Materials. Weinheim: WILEY-VCH. 2009. P. 109-136.
9. On the origin of the extremely high strength of ultrafine-grained al alloys produced by severe plastic deformation / R. Z. Valiev [et al.] // Scripta Materialia. 2010. V. 63. P. 949-952.
10. Pande C. S., Cooper K. P. nanomechanics of hall-petch relationship in nanocrystalline materials // Progress in Materials Science. 2009. V. 54. P. 689-706.
11. Малыгин Г. А. Пластичность и прочность микро- и нанокристаллических материалов // Физика твердого Тела. 2007. Т. 49, № 6. C. 961-982.
12. Штремель М. А. Прочность сплавов. Часть
II. Деформация: Учеб. для вузов. М.: МИСИС, 1997. 527 с.
13. Zhu Y. T., Han B. Q., Lavernia E. J. Defor-
mation mechanisms of nanostructured materials // Bulk Nanostructured. Weinheim: WILEY-VCH, 2009.
P. 89-108.
14. Мурашкин М. Ю., Кильмаметов А. Р., Валиев Р. З. Особенности структуры и механические свойства алюминиевого сплава 1570, подвергнутого интенсивной пластической деформации // Физика металлов и металловедение. 2008. Т. 106, № 1. С. 93-99.
15. Влияние легирования на распад алюминиево-магниевого твердого раствора / Л. Н. Гусева [и др.] // Известия АН СССР. Металлы. 1972. № 4. С. 208-213.
16. Nanostructure and related mechanical properties of an Al-Mg-Si Alloy Processed by severe plasric deformation / G. Nurislamova [et al.] // Philososphical Magazine Letters. 2008. V. 88, № 6. P. 459-466.
17. Enhancement of strength and superplasticity in a 6061 Al alloy processed by equal-channel-angular-pressing / W. J. Kim [et al.] // Metallurgical and Materials Transactions A. 2002. V. 33. P. 3155-3164.
18. Alexandrov I. V., Chembarisova R. G. The analysis of the effect of microstructure parameters on the ability to develop high - strength states and peculiarities of deformation states of Al 1570 alloy // Rev. Adv. Mater. Sci. 2010. V. 23, № 3. P. 209-218.
19. Curtin W. A., Olmsted D. L., Hector L. G. A predictive mechanism for dynamic strain ageing in alu-minum-magnesium alloys // Jr. Nature materials. 2006. V. 5, № 11. P. 875-880.
20. Штремель М. А. Прочность сплавов. Ч. I. Дефекты решетки: Учеб. пособие для вузов. М.: Металлургия, 1982. 280 с.
21. Новиков И. И. Дефекты кристаллического строения металлов: Учеб. пособие. М.: Металлургия, 1983. 232 с.
22. Фридель Ж. Дислокации. М.: Мир, 1967. 643 с.
23. Флейшер Р., Хиббард У. Упрочнение при образовании твердого раствора // Структура и механические свойства металлов: матер. конф., проведенной Науч. физ. лабораторией, Треддингтон, Мидлсекс, январь 1963 г. М.: Металлургия, 1967. 384 с.
24. Yan J. Strength modelling of Al-Mg-Cu type alloys: thesis for PhD. University of Southampton, 2006. 215 p.
25. Хирт Д., Лоте И. Теория дислокаций. М.: Атомиздат, 1972. 600 с.
26. Kocks U. F., Mecking H. Physics and phenomenology of strain hardening: the FCC case // Progress in Materials Science. 2003. V. 48. P. 171-273.
27. Valiev R. Z., Gertsman V. Yu., Kaiby-shev O. A. Non-equilibrium state and recovery of grain boundary structure. II. Energetics analysis // Phys. Stat. Sol. (a). 1983. V. 78. P. 177-186.
28. Alexandrov I. V., Chembarisova R. G. Development and application of the dislocation model for analysis of the microstructure evolution and deformation behavior of metals subjected to severe plastic deformation // Rev. Adv. Mater. Sci. 2007. V. 16, № 1/2. P. 51- 72.
29. Prinz F. B., Argon A. S. The evolution of plastic resistance in large strain plastic flow of single phase subgrain forming metals // Acta Metall. 1984. V. 32, № 7. P. 1021-1028.
30. Zehetbauer M. Cold work hardening in stages IV and V of fcc metals II: Model fits and physical results // Acta mater. 1993. V. 41 P. 589-599.
31. Малыгин Г. А. Анализ структурных факторов, определяющих образование шейки при растяжении металлов и сплавов с ГЦК-решеткой // Физика твердого тела. 2005. Т. 47, №2. C. 236-241.
32. Дубинко В. И., Клепиков В. Ф. Влияние неравновесных флуктуаций на пластичность металлов под облучением // Вюник Харьювського университету. 2005. Т. 710. C. 87-92.
33. Systematic study of stacking fault energies of random Al-based alloys / T. C. Schulthess [et al.] // Acta Materialia. 1998. V. 46, № 6. P. 2215-2221.
34. Мартин Д. У. Микромеханизмы дисперсионного твердения сплавов. М.: Металлургия, 1983. 167 с.
35. Маркушев М. В., Мурашкин М. Ю. Прочность и трещиностойкость промышленных алюминиевых сплавов 1560 и 5083 системы А1-М^-Мп после интенсивной пластической деформации угловым прессованием // Физика металлов и металловедение. 2004. Т. 98, № 2. С. 1-12.
36. Чувильдеев В. Н. Неравновесные границы зерен в металлах. М.: ФИЗМАТЛИТ, 2004. 304 с.
3 7 . Бернер Р., Кронмюллер Г. Пластическая деформация монокристаллов. М.: Мир, 1969. 272 с.
ОБ АВТОРАХ
Александров Игорь Васильевич, проф., зав. каф. физики, гл. науч. сотр. НИЧ УГАТУ. Дипл. физик (БГУ, 1976). Д-р физ.-мат. наук по физике твердого тела (ИФМ УрНЦ РАН, 1997). Иссл. в обл. физики прочности и пластичности материалов.
Чембарисова Роза Галиевна, доц., ст. науч. сотр. НИЧ УГАТУ. Дипл. физик (БГУ, 1980). Канд. техн. наук (Азерб. Индустриальн. ун-т им. М. Азизбекова, 1992). Иссл. в обл. нелинейной динамики, физики прочности и пластичности материалов.
Латыпов Марат Ильдарович, дипл. инженер-физик (УГАТУ, 2011). Иссл. в обл. физики прочно-
сти и пластичности материалов.
