CC BY
24
Скрипняк В.А., Скрипняк В.В., Козулин А.А., Иохим К.В. Влияние концентраторов напряжений на механическое поведение магниевого сплава при высокоскоростной деформации в температурном диапазоне от 295 К до 673 К // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2019. - № 1. - С. 151-160. DOI: 10.15593/perm.mech/2019.1.13
Skripnyak V.A., Skripnyak V.V., Kozulin A.A., Iohim K.V. The influence of stress concentrators on the magnesium alloy mechanical behavior under deformation at high strain rates in the temperature range from 295 to 673 K. PNRPUMechanics Bulletin, 2019, no. 1, pp. 151-160. DOI: 10.15593/perm.mech/2019.1.13
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 1, 2019 PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
Б01: 10.15593/регш.шееЬ/2019.1.13 УДК 539.3: 539.375: 539.4.019.1: 539.42
ВЛИЯНИЕ КОНЦЕНТРАТОРОВ НАПРЯЖЕНИЙ НА МЕХАНИЧЕСКОЕ ПОВЕДЕНИЕ МАГНИЕВОГО СПЛАВА ПРИ ВЫСОКОСКОРОСТНОЙ ДЕФОРМАЦИИ В ТЕМПЕРАТУРНОМ ДИАПАЗОНЕ ОТ 295 ДО 673 К
В.А. Скрипняк, В.В. Скрипняк, А.А. Козулин, К.В. Иохим
Национальный исследовательский Томский государственный университет, Томск, Россия
АННОТАЦИЯ
В работе представлены результаты экспериментальных исследований механического поведения сплава Mg-3% Al-1% Zn при высокоскоростном растяжении при комнатной и повышенных температурах. Для испытаний применялись плоские образцы с гладкой рабочей частью и с надрезами радиусом 10,5 и 2,5 мм. Экспериментальные исследования проведены на высокоскоростном сервогидравлическом стенде «Инстрон» VHS 40/50-20. Нагрев образцов плоскими керамическими инфракрасными излучателями до заданных температур занимал в среднем от 60 до 160 с. Контроль за температурой в рабочей части образцов осуществлялся в реальном времени с помощью хромель-алюмелевой термопары. Получены данные о влиянии параметра трехосности напряженного состояния в диапазоне от 0,33 до 0,5 на напряжение течения и величину предельной деформации до разрушения магниевого сплава при растяжении со скоростями деформации 100 и 1000 1/с при температурах 295, 473 и 673 К. Обнаружено, что с ростом параметра трехосности напряженного состояния в диапазоне от 0,33 до 0,5 величина предельной деформации до разрушения при растяжении сплавов магния уменьшается в два раза. Указанный эффект реализуется в широком диапазоне скоростей деформации и температурах T/Tm от 0,32 до 0,73 (Tm = 923 K - температура плавления магниевого сплава Mg-3% Al-1% Zn). Полученные данные были использованы для калибровки моделей пластического деформирования и разрушения магниевого сплава МА2-1.
©ПНИПУ
О СТАТЬЕ
Получена: 11 февраля 2019 г. Принята: 22 марта 2019 г. Опубликована: 30 марта 2019 г.
Ключевые слова:
высокоскоростная деформация, повреждение, разрушение при растяжении, параметр трехосности напряженного состояния, повышенные температуры, магниевые сплавы.
© Скрипняк Владимир Альбертович - д. ф.-м. н., проф., зав. каф., email: skrp2006@yandex. ru, : 0000-0001 -7162-3983 Скрипняк Владимир Владимирович - к. ф.-м. н., с.н.с., e-mail: skrp2012@yandex.ru, ID: 0000-0003-2238-5856 Козулин Александр Анатольевич - к. ф.-м. н., доц., e-mail: kozulyn@ftf.tsu.ru Иохим Кристина Владимировна - магистрант, e-mail: super-girl_1996.20@mail.ru
Vladimir A. Skripnyak - Doctor of Physical and Mathematical Sciences, Professor, Head of the Department, e-mail: skrp2006@yandex.ru, : 0000-0001 -7162-3983
Vladimir V. Skripnyak - CSc of Physical and Mathematical Sciences, Senior Researcher, e-mail: skrp2006@yandex.ru, ¡D: 0000-0003-2238-5856
Alexander A. Kozulyn - CSc of Physical and Mathematical Sciences, Associate Professor, e-mail: kozulyn@ftf.tsu.ru Kristina V. Iohim - Postgraduate Student, e-mail: super-girl_1996.20@mail.ru
Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
THE INFLUENCE OF STRESS CONCENTRATORS ON THE MAGNESIUM
ALLOY MECHANICAL BEHAVIOR UNDER DEFORMATION
AT HIGH STRAIN RATES IN THE TEMPERATURE RANGE FROM 295 TO 673 K
V.A. Skripnyak, V.V. Skripnyak, A.A. Kozulin, K.V. Iohim
National Research Tomsk State University, Tomsk, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 11 February 2019 Accepted: 22 March 2019 Published: 30 March 2019
Keywords:
high-strain rates, damage, fracture under tension, the stress triaxiality factor, elevated temperatures, and magnesium alloys.
The paper presents the experimental results of the mechanical behavior of Mg-3% Al-1% Zn alloy at high strain rates at room and elevated temperatures. The flat samples with smooth working parts and notches with a radius of 10 mm, 5 mm and 2.5 mm were used. The experimental studies were carried out using the high-velocity servo hydraulic test machine Instron vHs 40 / 5020. The samples were heated with flat ceramic infrared emitters on average from 60 seconds to 160 seconds. The temperature control in the working part of the samples was carried out in real time using a chromel-alumel thermocouple.
Data on the stress triaxiality effect on the plastic flow stress and the strain to fracture of the magnesium alloy under tension with strain rates of 100 and 1000 1 / s at temperatures of 295 K, 473 K and 673 K were obtained. The stress triaxiality factor was varied within the range from 0.33 to 0.5.
It was found that the value of tensile strain to fracture of the magnesium alloy decreased twice when the stress triaxiality factor increased from 0.33 to 0.5. This effect is realized in a wide range of strain rates and homologous temperatures T/Tm from 0.32 to 0.73 (Tm= 923 K is the melting point of Mg-3% Al - 1% Zn magnesium alloy). The obtained data were used for calibration of plastic deformation and fracture models which are used in computer-aided design of products from MA2-1magnesium alloy.
©PNRPU
Введение
Легкие сплавы на основе магния, относящиеся к изомеханической группе сплавов с гексагональной плотноупакованной решеткой, обладают повышенными удельными прочностными характеристиками и применяются для создания металлических легких и надежных конструкций [1-4]. Повреждаемость сплавов, обусловливающая развитие вязкого разрушения, накладывает ограничения на технологические режимы формования изделий, равноканального углового прессования, рифления при прессовании для получения сплавов с ультрамелкозернистой структурой [5-7]. Развитие вязкого разрушения магниевых сплавов в широком интервале температур и скоростей деформации является основой технологических процессов обработки изделий, включая резание и раскрой листового проката [8-9]. Для прогноза вязкого разрушения и влияния повреждений на пластическую деформацию магниевых сплавов предложены многопараметрические модели и критерии разрушения в широких условиях нагружения [10-21]. Получены экспериментальные данные о закономерностях деформации магниевых сплавов в широком диапазоне температур и скоростей нагружения [2, 3, 8-22]. Установлено, что магниевые сплавы обладают различным сопротивлением пластическим деформациям при растяжении и сжатии. Этот эффект обусловлен механизмами пластического течения, связанными с двойни-кованием [19-21]. Обнаружено, что в магниевых сплавах при пластической деформации возникает деформа-
ционная анизотропия, обусловленная формированием текстуры [18-20]. При пластической деформации магниевых сплавов в температурном диапазоне от 423 до 673 К происходит интенсивная динамическая рекристаллизация, сопровождающаяся изменением механических свойств [16, 22]. Экспериментальные данные свидетельствуют о том, что механическое поведение магниевых сплавов в широком интервале температур и скоростей деформации отличается от поведения сплавов, относящихся к другим изомеханическим группам [1]. В связи с этим актуальной задачей является разработка адекватных моделей механического поведения магниевых сплавов, учитывающих закономерности пластического течения и разрушения.
В настоящее время большое внимание уделяется разработке моделей со связанными критериями повреждения и пластичности [24-33].
В этих моделях описание пластической деформации связано с расчетными значениями параметра повре-жденности, зависящего от относительного объема пор в материале. Данный класс моделей развивается в рамках подхода Гурсона - Твергарда-Нидельмана [25-30]. Для калибровки этих моделей используются экспериментальные данные, полученные при испытаниях образцов на растяжение, сжатие, сдвиг, продавливание [28-29]. Методики испытаний образцов с различной геометрией были разработаны для получения экспериментальных данных о влиянии трехосности напряженного состояния на пластическое течение и развитие повреждений в области положительных напряжений [20, 28, 32-35].
Цель данной работы заключалась в получении новых экспериментальных данных о влиянии параметра трехосности напряженного состояния на механическое поведение магниевых сплавов Mg-Al-Zn в условиях высокоскоростной деформации и в широком диапазоне температур для калибровки определяющих соотношений, а также в получении новых экспериментальных данных о скоростной чувствительности напряжения течения сплава при повышенных температурах.
1. Материал и методы исследования
1.1. Материал
Исследования проведены на коммерческом тонколистовом прокате магниевого сплава МА2-1 (аналог сплава AZ31 или Mg - 3%Al - 1%Zn). Сплав находился в поликристаллическом состоянии. Зеренная структура сплава исследована на сканирующем микроскопе TESCAN методом дифракции вторичных электронов (EBSD). Анализ полученных сканов проведен с помощью лицензионного программного обеспечения HKL Channel 5 Oxford Instruments. Структура являлась мелкокристаллической, однородной и равноосной. Средний размер зерна сплава составлял 40 мкм.
Образцы для испытаний вырезались электроискровым методом на станке с числовым программным управлением DK7732, с ориентацией оси вдоль направления проката из листа толщиной 2,4 мм, что обеспечивало высокую точность воспроизводства геометрических параметров рабочей части образцов. Геометрические параметры образцов для высокоскоростных испытаний были выбраны в соответствии с рекомендациями INSTRON и стандартов ASTM E8 и DIN 50125-H. Геометрия образцов показана на рис. 1.
б
Рис. 1. Плоские образцы магниевого сплава с надрезами: а - конфигурация надрезов рабочей части образцов; б - фотографии образцов с надрезами радиусом 2,5, 5, 10 мм
Fig. 1. Flat samples of the magnesium alloy with notches: geometry of notches (a); photos of the samples with notches having a radius of 2.5 mm, 5 mm, and 10 mm, respectively (b)
1.2. Механические испытания и анализ данных
Для исследования механического поведения тонколистового проката магниевого сплава были проведены испытания по осевому растяжению образцов с постоянной скоростью деформации. Испытания проведены на высокоскоростном испытательном стенде INSTRON 8800 (конфигурация VHS 40/50-20). Для высокоскоростных испытаний применялись образцы, показанные на рис.1. Для проверки воспроизводимости результатов испытывали не менее пяти образцов с одинаковой конфигурацией рабочей части при комнатной температуре (295 К) и при повышенных начальных температурах (473±2) К и (673±2) К. Повышенные температуры во время испытаний были достигнуты путем нагрева рабочей части образцов в теплоизолированном боксе с помощью плоских инфракрасных керамических нагревателей. Нагрев образцов до заданных начальных температур занимал в среднем от 60 с до 160 с, что препятствовало протеканию рекристаллизации сплава. Однородность поля температуры в рабочей части образца обеспечивалась применением потоков инфракрасного излучения по нормали к плоскости образца. Температура на границе рабочей части образца измерялась непрерывно в реальном времени с помощью контактной хромель-алюмелевой термопары с разрешением 0,1 K и точностью 1,5 K.
Испытания выполнены при постоянных скоростях растяжения 20 и 2 м/с, которые обеспечивали скорости деформации 1000 с-1 и 100 с-1 соответственно. Скорость растяжения поддерживалась постоянной с точностью до 0,01 м/с и контролировалась в процессе нагружения образцов в каждом эксперименте. Механические отклики, полученные в результате испытаний на растяжение, получены в виде зависимостей усилий деформирования от перемещений. Эти данные были использованы для определения зависимостей истинных напряжений от истинных деформаций в соответствии со стандартом ASTM E8.
Истинное напряжение Ctu определялось по формуле
АР
(1+ А /„/10),
(1)
где AP - приращение усилия; Fo - площадь начального поперечного сечения рабочей части образца; Д^ -удлинение рабочей части образца; ^ - начальная длина рабочей части образца.
Истинная деформация е^ определялась по формуле
е, = ln(1 + А1 /L ).
true V cp 0 '
(2)
Относительное равномерное удлинение (5p) [%] определялась по формуле
Sp = ln(A cp/ lo)100%
(3)
где Дф - удлинение рабочей части образца до образования шейки; lo - начальная длина рабочей части образца.
а
=
true
При проведении испытаний на растяжение датчики позиционирования и ускорения испытательного стенда Instron VHS 40/50-20 обеспечивали определение пластических деформаций в диапазоне от 0 до 0,2 с точностью 0,5 %. В интервале пластических деформаций от 0,2 до 1 точность определения пластической деформации составляла ~2 % от измеренной величины.
Для определения деформации до разрушения е/ плоских образцов применялась формула
=7273^8? + 8? +[-(ej + 82 )]2
(4)
где
8 = ln
w,
s2 = ln
vSo
w1 - наименьшая изме-
ренная ширина сечения в области разрушения; -начальная ширина рабочей части образца в зоне надреза; 51 - среднее арифметическое измеренных значений толщины в отдельных сечениях области разрушения; 5о -начальная толщина рабочей части образца в зоне надреза.
Величины w1 и 51, 50 определялись с использованием измерительного цифрового микроскопа БМЦ-1М при анализе фотографий зоны разрушения образцов (рис. 2, а, б). Ширина зоны и толщина образца были измерены до и после разрушения образца с точностью 1 мкм (см. рис. 2, в). Толщина 51 определялась как среднее арифметическое измеренных значений толщины в отдельных сечениях области разрушения.
При определении деформации до разрушения е/ плоских образцов по формуле (4) предполагалось выполнение постулата пластической несжимаемости и малости упругих деформаций по сравнению с пластическими:
8+^+8 = 0, 8 =-(8 +е2),
(5)
где е, - главные компоненты тензора деформаций, г =
= 1, 2, 3.
Истинные напряжения рассчитывались для каждого испытания вплоть до разрушения образца. После начала формирования шейки принималось во внимание, что напряженное состояние становится неоднородным в рабочей части образцов из-за изменения формы рабочей части.
Параметр трехосности напряженного состояния п и параметр Лоде Ь определялись по формулам [14]
п = -р/Оед, Ь = (2 Оц - О1 - Сш)/(С1 - Ош), (6)
где а: = а» , ап = (1/2)(ст„.а. -а.а.) , стш = ё^а) =
= СТ11СТ22СТ33 + 2а12а23а31 -а22а33 — а2з а11 — ама22 - первый, второй и третий инварианты тензора напряжений Коши соответственно; р - давление, р = - о/3; а -
эквивалентное напряжение, aeq Gj -1 ст»att =
=J3стп ; Qij - компоненты тензора напряжений Коши.
д е ж з и
Рис. 2. Конфигурация образцов магниевого сплава в зоне разрушения: а - схема оцифровки области разрушения; б - контур области разрушения; в - схема измерения геометрических параметров зоны разрушения; г - начальная конфигурация сечения рабочей части образцов; фотографии зоны разрушения образцов с надрезами, имеющими радиус; д - 2,5 мм; е - 5 мм; ж - 10 мм;
з - с гладкой частью при 473 К; и - с гладкой частью при 673 K
Fig. 2. The geometry of the magnesium alloy samples in the zone of fracture: the scheme of digitizing the fracture area (a), a contour the fracture area (b), geometric parameters of the fractured zone and initial geometry of cross section samples (c); initial configuration of samples; photographs of the specimens fracture zones with notches with radius of 2.5 mm (d), f 5 mm (e), 10 mm f), with smoothed
specimen at 473K (g), and with a smoothed specimen at 673 K (h)
б
а
в
г
Начальное значение параметра трехосности напряженного состояния п оценивалось по аналитической формуле для образцов с надрезами в плоском напряженном состоянии [29]:
■Л=(1 + 2a)/(b4ar+a + i), a = ln[1 + w/(4Д)], (7)
где w - ширина образца; R - значение радиуса кривизны надреза.
Значения п при радиусах надрезов 10, 5, 2,5 мм и для образцов с рабочей частью без надрезов равны 0,3997; 0,4455; 0,5029; 0,33 соответственно.
В рамках подхода Гурсона - Твергаарда-Нидель-мана критерий пластичности имеет вид [25, 26]
(a2g / о2) + 2qif cosh (-q p /2а,) -1 - q3 (f )2 = 0, (8)
где Qeq - эквивалентное напряжение; os - напряжение течения; p - давление; q1 = 1,5, q2 = 1 и q3 = 2,25 - параметры материала; f - параметр поврежденности материала.
При отсутствии повреждений критерий (8) эквивалентен критерию пластичности Мизеса. В случае решения динамических задач девиатор тензора напряжений определяется из решения уравнения
пластическая
щ 2%¿
Dt Зе„„ ip
(9)
где D/Dt - производная по времени Яуманна; Sy - де-виатор тензора напряжений Коши, Sv = - (1 / 3)ст№81;/ ; iy ~ компоненты тензора скорости деформации.
Определение параметров материала в уравнениях (7), (8) и (9) осуществляется с использованием результатов экспериментальных исследований механического поведения в соответствующих диапазонах изменения температуры, скорости деформации и параметра трех-осности напряженного состояния.
Напряжения течения для сплавов с гексагональной плотноупакованной решеткой (ГПУ) описывались в виде суммы атермического напряжения и теплового напряжения [27]. Модификация определяющего уравнения Зерилли-Армстронга для магниевых сплавов использована в форме
- С2 ехр {-Сзг + САТ ln (ё^ / ёед0)},
(10)
где ctj0 , C, П , Кр, C, C, C - параметры материала;
dg - средний размер зерна; ё = [(2 / 3)ё,..ё,.. ] ~; 8^ = 1,0^; < = £[(2/3) ЦГЛ -
«7 Lv ' V V -
эквивалентная
деформация, T = T0 + (0.9/pCp)J a¿Jdsp ;
To - начальная температура нагружения; р - массовая плотность; Ср - удельная теплоемкость.
Атермическое напряжение, задаваемое вторым слагаемым в правой части уравнения (10), описывается соотношениями без учета зависимости от скорости деформации и температуры, но с учетом зависимости от степени пластической деформации, отражающей структурную эволюцию сплавов.
Связь между эквивалентной пластической деформацией и параметром повреждения в случае вязкого разрушения материала задается функцией пластического потенциала, которая зависит от параметра повре-жденности материала /:
f '/nucí ^/growth '
Lci =(f,KK exp{-0,5[sf, -s,)/^]2}, (11)
/growth — У)
где f - параметр поврежденности; е^ и sn представляют собой средние деформации нуклеации и стандартное отклонение соответственно; fN - параметр, контролирующий количество зарождающихся пор.
При использовании моделей поврежденной среды локальное разрушение происходит при достижении эквивалентной деформацией критического значения, величина которого зависит от температуры, скорости деформации и параметра трехосности напряженного состояния. Для прогноза предельной деформации до разрушения ef магниевых сплавов, относящихся к изо-механической группе сплавов с ГПУ-решеткой, использовались феноменологические соотношения
Е/= Ъ (éeJé0)F2[(T-TREF)/(Tm-TREF)l (12)
где F\, /•':. /':, - феноменологические функции; s0 = 0,01 сГ1; Т и '/„, - температура деформации и температура плавления соответственно; Tref = 295 К.
В данной работе были приняты следующие варианты феноменологических функций:
F, = А +Д ехр(/уо. F2 =(l + (É4 /ёХ), F3=(l + D5(T-TllEF)/(Tm-TREF)),
где Di, D2, D3, D4, D5 - коэффициенты соотношений, являющиеся параметрами материала.
Для определения коэффициентов Di, D2, D3, D4, D5 в соотношениях (13) необходимы экспериментальные данные о механическом поведении материала при варьировании скорости деформации, температуры и параметра трехосности напряженного состояния.
2. Результаты и обсуждение
На рис. 3 показаны полученные при обработке экспериментальных данных зависимости истинных напряжений от истинных деформаций при растяжении образ-
(13)
цов магниевого сплава Mg-3%Al-1%Zn с гладкой рабочей частью (линии 1), с рабочими частями, имевшими надрезы с радиусом 10 мм (линии 2), 5 мм (линии 3), 2,5 мм (линии 4).
Для каждого условия нагружения проводилось по 20 испытаний (4 типа геометрии рабочей части по 5 образцов).
0.05 0,10 0,15 0,20 0,25 0,30 Истинная деформация а
300 250 200 150 100 50 0
100 с-', 475 К
\ \ \
1 1 1 1 1 I
г I'1 4 1 \ ■ 2 I 1
1' / ll 1 • 1 1
Осциллирующий характер диаграмм деформирования для магниевого сплава МА2-1, показанных линиями 1 на рис 3, а, б, в, аналогичен волнообразному характеру диаграмм деформирования, зарегистрированных для сплава AZ31B в [30]. Формирование осциллирующего характера диаграмм деформирования обусловлено задержкой пластического течения.
250 200 150 100 50
0,0 0,1 0,2 0.3 Истинная деформация
0 0,0
100 с1, 673 К
\ \
// ч 7 \ У \ Л
/ 4 2\\!
0,1 0,2 0,3 Истинная деформация
0.4
Рис. 3. Зависимости истинных напряжений от истинных деформаций при нагружении образцов магниевого сплава Mg-3%Al-1% Zn со скоростью деформации 1000 с-1 - (а), (б) и 100 с-1 - (в), (г); линии 1 - для образцов с гладкой частью, линии 2 - для образцов с надрезами, имеющими радиус 10 мм; линии 3 - для образцов с надрезами, имеющими радиус 5 мм; линии 4 - для образцов
с надрезами, имеющими радиус 2,5 мм
Fig. 3. True stress - true strain curves of Mg - 3% Al - 1% Zn magnesium alloy at strain rates of 100 s-1 and 1000 s-1, the samples with a smooth part (a), the samples with notches having a radius of 10 mm (b), the samples with notches radius of 5 mm (c), (g), the samples with notches
radius of 2.5 mm (d)
При скоростях деформации, меньших 100 с-1, осциллирующего характера диаграмм деформирования образцов с гладкой рабочей частью не наблюдается. Осциллирующий характер диаграмм деформирования для образцов с надрезами не наблюдался при всех исследованных режимах нагружения. Анализ полученных результатов, приведенных на рис. 3, свидетельствует о том, что при увеличении скорости деформации от 100 с-1 до 1000 с-1 предел текучести сплава МА2-1 при растяжении увеличивается почти в два раза. Высокая скоростная чувствительность напряжения течения магниевого сплава получена для исследованного диапазона температур от 473 до 673 К. Рост скоростной чувствительности напряжения течения магниевых сплавов с увеличением температуры в диапазоне от 473 до 673 К обусловлен сменой доминирующих механизмов пластического течения с базисного скольжения дислокаций при температуре 473 К на пирамидальное поперечное скольжение дислокаций при температуре 673 К [21].
Диаграммы истинных напряжений от истинных деформаций, показанные на рис. 3 для МА2-1, качественно подобны диаграммам, полученным в [13, 30]. Усредненные значения напряжения течения для сплавов AZ31B и МА2-1 при скорости деформации 1000 с-1 при температуре 473 К различаются на 12 %, а деформации до разрушения - менее чем на 16 %.
Полученные экспериментальные диаграммы истинных напряжений от истинных деформаций, показанные на рис. 3, были использованы для определения численных значений параметров материала в уравнении (10). При определении компонент тензора пластической деформации использовались значения модуля Юнга, Е = 43 [1 - 5,3 10-4(Г 295 К)] ГПа [1]. При анализе результатов экспериментов по высокоскоростному растяжению образцов магниевого сплава Mg-3%Al-1%Zn, температура в уравнении (10) оценивалась с использованием линейной зависимости удельной теплоемкости Ср от температуры в интервале от 1005,9 Дж/кг К при
б
в
г
295 К до 1130 Дж/кг K при 523 К [13]. Начальная массовая плотность ро принималась равной 1780 кг/м3. При определении параметров уравнения (10) было использовано значение khp, полученное ранее в [37]. Полученные значения параметров модифицированной модели Зерилли-Армстронга приведены в таблице.
Параметры модели Зерилли-Армстронга для магниевого сплава Mg-3%Al-1%Zn
Zerilli-Armstrong model parameters for magnesium alloy Mg-3%Al-1%Zn
Параметр модели ГПа khp, Па нм112 C2, ГПа C3, K-1 C4, K-1 C5, ГПа П1
0,141 6,2 0,315 0,0029 0,000389 0,505 0,2514
Параметры модельных уравнений кинетики повреждения (11) подбирались для получения согласия прогнозируемых и экспериментальных значений напряжения течения на стадии предразрушения. Были получены следующие значения параметров ед = 0,1, sn = 0,2, /д = 0,4.
Значения деформации до разрушения материала е/, полученные при анализе конфигурации зон разрушения образцов, испытанных в указанных на рис. 1 условиях нагружения, были использованы для определения численных значений коэффициентов соотношений (13), являющихся параметрами материала.
Для определения коэффициентов A, D2, D3, D4, D5 в соотношениях (13) был использован метод регрессионного анализа, реализованный в стандартной процедуре User-Defined Regression Wizard в пакете SigmaPlot 12.5. Показано, что полученные экспериментальные данные могут быть описаны соотношениями (12)—(13).
Зависимость деформации до разрушения для магниевых сплавов Mg-3%Al-1%Zn принимает вид
ef = (0,0199 + 0,5 e^4n)(l + (ée(; /0,01)"0'0355)) x(l +1,6 (T - 295) / (923 - 295)),
(14)
где Т - абсолютная температура; ё - эквивалентная
скорость деформации, имеющая размерность с-1.
На рис. 4 показаны зависимости деформации до разрушения от параметра трехосности напряженного состояния при заданных температуре и скорости деформации, описываемые соотношением (14).
Полученные в данной работе результаты не имеют точных аналогов, опубликованных ранее в реферируемых журналах.
Отметим, что полученные в данной работе зависимости истинных напряжений от истинных деформаций для МА2-1 качественно согласуются с аналогичными диаграммами для Л231, полученными авторами [16] при испытаниях цилиндрических образцов магниевого сплава А231 ^-2,9%Л1-0,85%2п-0,3%Мп) с гладкой рабочей частью в диапазоне температур от 423 до 673 К со скоростями деформации 10-3, 10-2, 10-1, 1, 10, 100 с-1 при растяжении. Для скорости деформации 100 с-1
напряжения течения сплавов Л231Б и МЛ2-1 в диапазоне деформаций до 20 % при температурах 473 и 673 К отличаются менее чем на 5,5 %. Предельные деформации до разрушения в указанных условиях отличаются на 27 %. Это может быть обусловлено различиями в зеренных структурах сплавов Л231Б и МЛ2-1.
Рис. 4. Зависимость предельной деформации до разрушения при растяжении 8/ магниевого сплава Mg-3%Al-1%Zn от параметра трехосности напряженного состояния. Линии 1, 2 -температуры 673 и 473 К соответственно; скорость деформации 100 с-1; 3 и 4 - температуры 673 и 473 К соответственно, скорость деформации 1000 с-1; 5 - скорость деформации 0,001 с-1, температура 295 К. Символами показаны экспериментальные данные, полученные в [16, 22, 30] и в данной работе
Fig. 4. Strain to fracture depending on the triaxiality stress factor of Mg - 3% Al - 1% Zn magnesium alloy under tension predicted at different temperatures and strain rates (curve 1) - 100 s-1, 673 K, (2) - 100 s-1, 473 K, (3) - 1000 s-1, 673 K, (4) - 1000 s-1, 473 K, and (5) 0.001 s-1 and 295 K. The symbols show the experimental data obtained in [16,22,30] and this research
Результаты, представленные на рис. 4, показывают, что в диапазоне исследованных температур при увеличении скорости деформации от 0,001 до 1000 с-1 значения деформация до разрушения при растяжении сплавов Mg-3%Al-1%Zn уменьшаются. С увеличением параметра трехосности напряженного состояния от 0,33 до 0,5 величина деформации до разрушения при растяжении сплавов магния Mg-3%Al-1%Zn в исследованных условиях нагружения уменьшается в два раза. Указанный эффект реализуется в широком диапазоне скоростей деформации от 0,001 до 1000 с-1 при гомологических температурах TITm от 0,32 до 0,73.
Заключение
Полученные экспериментальные данные свидетельствуют о том, что с ростом параметра трехосности напряженного состояния от 0,33 до 0,5 в диапазоне гомологических температур T/Tm от 0,32 до 0,73 деформация до разрушения сплавов магния Mg-3%Al-1%Zn уменьшается почти в два раза при скоростях деформации 100 до 1000 с-1.
При увеличении скорости деформации от 100 до 1000 с-1 предел текучести сплава МА2-1 при растяжении увеличивается почти в два раза при температурах 473 и 673 К.
Полученные экспериментальные данные о механическом поведении магниевого сплава МА2-1 при скоростях деформации 100 и 1000 с-1 и при температурах 473 и 673 К позволили определить значения параметров материала модифицированной модели Зерилли-Армстронга для ГПУ-сплавов.
Определены численные значения параметров материала для модели, позволяющей описывать зависимость деформации до разрушения магниевого сплава МА2-1 от скорости деформации, температуры и параметра трехосности напряженного состояния.
Полученные результаты могут быть использованы для разработки вычислительных моделей механического поведения конструкций из магниевых сплавов, под-
Библиографический список
1. Frost H.J., Ashby M.F. Deformation Mechanism Maps. -Oxford: Pergamon Press,1983. - 166 p.
2. Mechanical behavior and microstructural evolution of a Mg AZ31 sheet at dynamic strain rates / I. Ulacia, N.V. Duda-mell, F. Galvez, S. Yi, M.T. Perez-Prado, I. Hurtado // Acta Mate-rialia. - 2010. - Vol. 58. - P. 2988-2998.
3. Tensile characterization and constitutive modeling of AZ31B magnesium alloy sheet over wide range of strain rates and temperatures / I. Ulacia, C.P. Salisbury, I. Hurtado, M.J. Worswick // Journal of Materials Processing Technology. - 2011. - Vol. 211. -Iss. 5. - P. 830-839.
4. Johnson-Cook based criterion incorporating stress triaxiali-ty and deviatoric effect for predicting elevated temperature ductility of titanium alloy sheets / B. Valoppi, S. Bruschi, A. Ghiotti, R. Shivpur // International Journal of Mechanical Sciences. -2017. -Vol. 123. - P. 94-105.
5. Получение мелкодисперсных материалов на основе магния. Результаты численного моделирования и эксперимент / В.Н. Аптуков, П.В. Романов, Н.Е. Скрябина, Д. Фрушар // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2017. - № 3. - С. 5-16.
6. Влияние структуры на сопротивление пластической деформации алюминиевого сплава 1560 после обработки методом прессования рифлением / Е.Н. Москвичев, В.А. Скрип-няк, Д.В. Лычагин, А.А. Козулин, В.В. Скрипняк // Письма о материалах. - 2016. - Т. 6, № 2. - С. 141-145.
7. The flow behavior in as-extruded AZ31 magnesium alloy under impact loading / В. Zhu, X. Liu, C. Xie, W. Liu, C. Tang, L. Lu // Journal of Magnesium and Alloys. -2018. - Vol. 6. -P. 180-188.
8. Palumbo G., Sorgente D., Tricarico L. A numerical and experimental investigation of AZ31 formability at elevated temperatures using a constant strain rate test // Materials and Design. -2010. - Vol. 31. - P. 1308-1316.
9. Mechanical response and texture evolution of AZ31 alloy at large strains for different strain rates and temperatures / A.S. Khan, A. Pandey, T. Gnaupel-Herold, R.K. Mishra // International Journal of Plasticity. - 2011. - Vol. 27. - P. 688-706.
10. Effect of structural factors on mechanical properties of the magnesium alloy MA2-1 under quasi-static and high strain rate deformation conditions / G.V. Garkushin, S.V. Razorenov, V.A. Krasnoveikin, A.A. Kozulin, V.A. Skripnyak // Physics of the Solid State. - 2015. - Vol. 57. - No. 2. - P. 337-343.
11. Shock response of magnesium single crystals at normal and elevated temperatures / G.I. Kanel, G.V. Garkushin, A.S. Savi-
вергающихся динамическим воздействиям, пластическим деформациям при повышенных температурах.
Благодарности
Работа выполнена при частичной поддержке гранта РФФИ (проект №16-08-00037-а), гранта Президента РФ (МК-2690.2017.8).
Acknowledgements
This work was supported partly within a grant from the Russian Foundation for Basic Research (Project No. 16-08-00037-a), the Russian Presidential Grant (MK-2690.2017.8).
nykh, S.V. Razorenov, T. De Resseguier, W.G. Proud, M.R. Tyu-tin // Journal of Applied Physics. - 2014. - Vol. 116. - No. 14. -P. 143504.
12. Effect of temperature and heating rate on mechanical properties of magnesium alloy AZ31 / Y. Quan, C. Bin, L. Ji, Z. Ding-fei, Q. Guo-Zheng // Transactions Nonferrous Metal Society of China. - 2010. - Vol. 20. - P. 426-429.
13. Experimental study on tensile property of AZ31B magnesium alloy at different high strain rates and temperatures / F. Feng, S. Huang, Z. Meng, J. Hu, Y. Lei, M. Zhou, D. Wu, Z. Yang // Materials and Design. - 2014. - Vol. 57. - P. 10-20.
14. High temperature deformation and microstructural instability in AZ31 magnesium alloy / S. Spigarelli, O.A. Ruano, M. El Mehtedi, J.A. Del Valle // Materials Science and Engineering A. - 2013. - Vol. 570. - P. 135-148.
15. Hot compressive flow stress modeling of homogenized AZ61 Mg alloy using strain-dependent constitutive equations / H. Wu, J. Yang, F. Zhu, C. Wu // Materials Science and Engineering A. - 2013. - Vol. 574. - P. 17-24.
16. Effect of temperature and strain rate on the deformation behavior and microstructure of a homogenized AZ31 magnesium alloy / G. Bajargan, G. Singh, G. Sivakumar, U. Ramamurty // Materials Science and Engineering A. - 2013. - Vol. 579. - P. 26-34.
17. Chapuis A., Liu Q. Investigating the temperature dependency of plastic deformation in a Mg-3Al-1Zn alloy // Materials Science and Engineering A. - 2018. -Vol. 725. - P. 108-118.
18. Maksoud I.A., Ahmed H., Rödel J. Investigation of the effect of strain rate and temperature on the deformability and microstructure evolution of AZ31 magnesium alloy // Materials Science and Engineering A. - 2009. - Vol. 504. - P. 40-48.
19. Plastic deformation and ductility of magnesium AZ31B-H24 alloy sheet from 22 to 450 °C / A.R. Antoniswamy, E.M. Taleff, L.G. Hector Jr., J.T. Carter // Materials Science and Engineering A. - 2015. - Vol. 631. -P. 1 -9.
20. Selvarajou B., Joshi S.P., Benzerga A.A. Three dimensional simulations of texture and triaxiality effects on the plasticity of magnesium alloys // Acta Materialia. - 2017. - Vol. 127. -P. 54-72.
21. Strain rate sensitivities of deformation mechanisms in magnesium alloys / H. Wang, P. Wu, S. Kurukuri, M.J. Worswick, Y. Peng, D. Tang, D. Li // International Journal of Plasticity. -2018. - Vol. 107. - P. 207-222.
22. Twinning, dynamic recrystallization, and crack in AZ31 magnesium alloy during high strain rate plane strain compression
across a wide temperature / X. Liu, B.W. Zhu, C. Xie, J. Zhang, C.P. Tang, Y.Q. Chen // Materials Science and Engineering A. -2018. - Vol. 733. -P. 98-107.
23. Численное исследование эволюции напряженно- деформированного состояния структурно-неоднородного материала при одноосном нагружении / С.В. Смирнов, А.В. Коновалов, М.В. Мясникова, Ю.В. Халевицкий, А.С. Смирнов // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2016. - № 3. -
C. 175-187.
24. Verleysen P., Peirs J. Quasi-static and high strain rate fracture behaviour of Ti6Al4V// International // Journal of Impact Engineering. - 2017. - Vol. 108. - P. 370-388.
25. Neilsen K.L., Tvergaard V. Ductile shear failure or plug failure of spot welds modelled by modified Gurson model // Engineering Fracture Mechanics. - 2010. - Vol. 77. - P. 1031-1047.
26. Tvergaard V. Study of localization in a void-sheet under stress states near pure shear // International Journal of Solids and Structures. - 2015. - Vol. 60-61. -P.28-34.
27. Gao C.Y., Zhang L.C., Yan H.X. A new constitutive model for HCP metals // Materials Science and Engineering A. -2011. - Vol. 528. - P. 4445-4452.
28. Bai Y., Wierzbicki T. A new model of metal plasticity and fracture with pressure and Lode dependence // International Journal of Plasticity. - 2008. - Vol. 24. -P. 1071- 1096.
29. Selini N., Elmeguenni M., Benguediab M. Effect of the Triaxiality in Plane Stress Conditions // Engineering, Technology and Applied Science Research. - 2013. - Vol. 3. - No. 1. -P. 373-380.
30. A constitutive and fracture model for AZ31B magnesium alloy in the tensile state / F. Feng, S. Huang, Z. Meng, J. Hu, Y. Lei,
References
1. Frost H.J., Ashby M.F. Deformation Mechanism Maps. Oxford, Pergamon Press, 1983, 166 p.
2. Ulacia I., Dudamell N.V., Galvez F., Yi S., Perez-Prado M.T., Hurtado I. Mechanical behavior and microstructural evolution of a Mg AZ31 sheet at dynamic strain rates, Acta Materialia, 2010, vol.58, pp.2988-2998.
3. Ulacia I., Salisbury C.P., Hurtado I., Worswick M.J. Tensile characterization and constitutive modeling of AZ31B magnesium alloy sheet over wide range of strain rates and temperatures, Journal of Materials Processing Technology, 2011, vol.211. Issue 5, pp. 830-839.
4. Valoppi B., Bruschi S., Ghiotti A, Shivpur R. JohnsonCook based criterion incorporating stress triaxiality and deviatoric effect for predicting elevated temperature ductility of titanium alloy sheets, International Journal of Mechanical Sciences, 2017, vol.123, pp.94-105.
5. Aptukov V.N., Romanov P.V., Skryabina N.E., Fruchart
D. Polucheniye melkodispersnykh materialov na osnove magniya [Processing fine magnesium materials. numerical simulation and experimental analyses], PNRPU Mechanics Bulletin, 2017, n 3, pp. 5-16.
6. Moskvichev E.N., Skripnyak V.A., Skripnyak V.V., Kozu-lin A.A., Lychagin D.V. Influence of structure to plastic deformation resistance of aluminum alloy 1560 after groove pressing treatment, Letters on materials, 2016, vol. 6, n. 2, pp. 141-145.
7. Zhu В., Liu X., Xie C., Liu W., Tang C., Lu L. The flow behavior in as-extruded AZ31 magnesium alloy under impact loading, Journal of Magnesium and Alloys, 2018, vol.6, pp. 180-188.
M. Zhou, Z. Yang // Materials Science and Engineering A. -2014. - Vol. 594. - P. 334-343.
31. Skripnyak V.A., Skripnyak V.V., Skripnyak E.G. Ductility of titanium alloys in a wide range of strain rates // Proceedings IRF2018: 6th International Conference on Integrity-Reliability-Failure Lisbon/Portugal, 22-26 July 2018. - Porto, Portugal: FEUP-INEGI. - 2018. - P. 813-822.
32. Li Q. Mechanical properties and microscopic deformation mechanism of polycrystalline magnesium under highstrain-rate compressive loadings // Materials Science and Engineering A. - 2012. - Vol. 540. - P. 130-134.
33. Prasad N.S., Narasimhan R., Suwas S. Effect of notch acuity on the fracture behavior of AZ31 Mg alloy // Engineering Fracture Mechanics. - 2018. - Vol. 187. - P. 241-261.
34. Напряженно-деформированное состояние и разрушение элементов конструкций с острыми концентраторами напряжений при изгибе / С.Б. Сапожников, М.А. Иванов, С.И. Ярославцев, И.А. Щербаков // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2017. - № 4. - С. 40-55.
35. Mechanical Behavior and Microstructure Evolution of a Magnesium Alloy AZ31B Under Low Stress Triaxiality / H. Pan, F. Wang, L. Jin, M. Feng, J. Dong // Journal of Materials Science and Technology. - 2016. - Vol. 32. - P. 1282-1288.
36. Kubik P., Sebek F., Petruska J. Notched specimen under compression for ductile failure criteria // Mechanics of Materials 2018. - Vol. 125. - P. 94-109.
37. Failure mechanisms of light alloys with a bimodal grain size distribution / V.A. Skripnyak, E.G. Skripnyak, N.V. Skrip-nyak, I.K. Vaganova, V.V. Skripnyak // 11th. World Congress on Computational Mechanics (WCCM XI). 25-29, July, 2014. - Barcelona, Spain, 2014. - IV. - P. 3915-3925.
8. Palumbo G., Sorgente D., Tricarico L. A numerical and experimental investigation of AZ31 formability at elevated temperatures using a constant strain rate test, Materials and Design, 2010, vol.31, pp. 1308-1316.
9. Khan A.S., Pandey A., Gnaupel-Herold T., Mishra R.K. Mechanical response and texture evolution of AZ31 alloy at large strains for different strain rates and temperatures, International Journal of Plasticity, 2011, vol.27, pp.688-706.
10. Garkushin G.V., Razorenov S.V., Krasnoveikin V.A., Kozulin A.A., Skripnyak V.A Effect of structural factors on mechanical properties of the magnesium alloy MA2-1 under quasi-static and high strain rate deformation conditions, Physics of the Solid State, 2015, vol.57, № 2, pp.337-343.
11. Kanel G.I., Garkushin G.V., Savinykh A.S., Razorenov S.V., De Resseguier T., Proud W.G., Tyutin M.R. Shock response of magnesium single crystals at normal and elevated temperatures, Journal of Applied Physics, 2014, vol.116, pp.143504.
12. Quan Y., Bin C., Ji L. , Ding-fei Z., Guo- Zheng Q. Effect of temperature and heating rate on mechanical properties of magnesium alloy AZ31, Transactions Nonferrous Metal Society of China, 2010, vol.20, pp. s426-s429.
13. Feng F., Huang S., Meng Z., Hu J., Lei Y., Zhou M., Wu D., Yang Z. Experimental study on tensile property of AZ31B magnesium alloy at different high strain rates and temperatures . Materials and Design, 2014, vol.57, pp. 10-20.
14. Spigarelli S., Ruano O.A., El Mehtedi M., del Valle J.A. High temperature deformation and microstructural instability in
AZ31 magnesium alloy, Materials Science and Engineering A, 2013, vol. 570, pp. 135-148.
15. Wu H., Yang J., Zhu F., Wu C. Hot compressive flow stress modeling of homogenized AZ61 Mg alloy using strain-dependent constitutive equations, Materials Science and Engineering A, 2013, vol.574, pp. 17-24.
16. Bajargan G., Singh G., Sivakumar G., Ramamurty U. Effect of temperature and strain rate on the deformation behavior and microstructure of a homogenized AZ31 magnesium alloy, Materials Science and Engineering A, 2013, vol. 579, pp. 26-34.
17. Chapuis A., Liu Q. Investigating the temperature dependency of plastic deformation in a Mg-3Al-1Zn alloy, Materials Science and Engineering A, 2018, vol.725, pp. 108-118.
18. Maksoud I.A., Ahmed H., Rödel J. Investigation of the effect of strain rate and temperature on the deformability and microstructure evolution of AZ31 magnesium alloy, Materials Science and Engineering A, 2009, vol.504, pp. 40-48.
19. Antoniswamy A.R., Taleff E.M., Hector Jr. L.G., Carter J.T. Plastic deformation and ductility of magnesium AZ31B-H24 alloy sheet from 22 to 450 oC, Materials Science and Engineering A, 2015, vol. 631, pp. 1-9.
20. Selvarajou B., Joshi S.P., Benzerga A.A. Three dimensional simulations of texture and triaxiality effects on the plasticity of magnesium alloys, Acta Materialia, 2017, vol. 127, pp. 54-72.
21. Wang H., Wu P., Kurukuri S., Worswick M.J., Peng Y., Tang D., Li D. Strain rate sensitivities of deformation mechanisms in magnesium alloys, International Journal of Plasticity, 2018, vol.107. , pp. 207-222.
22. Liu X., Zhu B.W., Xie C., Zhang J., Tang C.P., Chen Y.Q. Twinning, dynamic recrystallization, and crack in AZ31 magnesium alloy during high strain rate plane strain compression across a wide temperature, Materials Science and Engineering A, 2018, vol.733, pp.98-107.
23. Smirnov S.V., Konovalov A.V., Myasnikova M.V., Khalevitsky Yu.V., Smirnov A.S. Chislennoye issledovaniye evo-lyutsii napryazhenno- deformirovannogo sostoyaniya strukturno-neodnorodnogo materiala pri odnoosnom nagruzhenii [A numerical study of stress-strain state evolution in structurally inhomoge-neous materials subjected to uniaxial loading], PNRPU Mechanics Bulletin, 2016. No. 3. pp. 175-187.
24. Verleysen P., Peirs J. Quasi-static and high strain rate fracture behaviour of Ti6Al4V, International Journal of Impact Engineering, 2017, vol.108, pp.370-388.
25. Neilsen K.L., and Tvergaard V. Ductile shear failure or plug failure of spot welds modelled by modified Gurson model, Engineering Fracture Mechanics, 2010, vol.77, pp.1031-1047.
26. Tvergaard V. Study of localization in a void-sheet under stress states near pure shear, International Journal of Solids and Structures, 2015, vol. 60-61, pp. 28-34.
27. Gao C.Y., Zhang L.C., and Yan H.X. A new constitutive model for HCP metals, Materials Science and Engineering A, 2011, vol.528, pp. 4445-4452.
28. Bai Y., Wierzbicki T. A new model of metal plasticity and fracture with pressure and Lode dependence, International Journal of Plasticity, 2008 , vol.24. , pp.1071- 1096.
29. Selini N., Elmeguenni M., Benguediab M. Effect of the Triaxiality in Plane Stress Conditions, Engineering, Technology and Applied Science Research, 2013, vol. 3, n. 1, pp. 373-380.
30. Feng F., Huang S., Meng Z., Hu J., Lei Y., Zhou M., Yang Z. A constitutive and fracture model for AZ31B magnesium alloy in the tensile state, Materials Science and Engineering A, 2014, vol. 594, pp. 334-343.
31. Skripnyak V.A., Skripnyak V.V., Skripnyak E.G. Ductility of titanium alloys in a wide range of strain rates, Proceedings IRF2018: 6th International Conference on Integrity-Reliability-Failure Lisbon/Portugal, 22-26 July 2018. FEUP-INEGI: Porto, Portugal. 2018, pp.813-822.
32. Li Q. Mechanical properties and microscopic deformation mechanism of polycrystalline magnesium under highstrain-rate compressive loadings, Materials Science and Engineering A. - 2012, vol.540. - P. 130-134.
33. Prasad N. S., Narasimhan R., Suwas S. Effect of notch acuity on the fracture behavior of AZ31 Mg alloy, Engineering Fracture Mechanics, 2018, vol.187, pp. 241-261.
34. Sapozhnikov S.B., Ivanov M.A., Yaroslavtsev S.I., Shcherbakov I.A. Napryazhenno-deformirovannoye sostoyaniye i razrusheniye elementov konstruktsiy s ostrymi kontsentratorami napryazheniy pri izgibe [The stress-strain state and failure of structural elements with sharp stress concentrators under bending], PNRPU Mechanics Bulletin, 2017, № 4, pp. 40-55.
35. Pan H., Wang F., Jin L., Feng M., Dong J. Mechanical behavior and microstructure evolution of a magnesium alloy AZ31B under low stress triaxiality, Journal of Materials Science and Technology, 2016, vol. 32, pp.1282-1288.
36. Kubik P., Sebek F., Petruska J. Notched specimen under compression for ductile failure criteria. Mechanics of Materials 2018, vol. 125 , pp.94-109.
37. Skripnyak V.A., Skripnyak E.G., Skripnyak N.V., Vaganova I.K., and Skripnyak V.V. Failure mechanisms of light alloys with a bimodal grain size distribution, in: 11th. World Congress on Computational Mechanics (WCCM XI), Barcelona, Spain, 25-29, July. 2014, 2014, IV, pp. 3915-3925.
