Серия 2. Технология машиностроения и материалы. зов при функционировании по причине потери устойчивости и разрушения.
Литература
1. Лясников А.В., Агеев Н.П., Кузнецов Д.П., Данилин Г.А., Дриго А.В. и др. Сопротивление материалов пластическому деформированию в приложениях к процессам обработки металлов давлением. СПб.: «Внешторгиздат-Петербург», 1995. 527 с.
2. Агеев Н.П., Данилин Г.А., Огородников В.П. Технология производства патронов стрелкового оружия. Ч.2. Процессы штамповки. СПб.: Балт. гос. техн. ун-т., 2006. 533 с.
3. Попов Е.А., Основы теории листовой штамповки. М.:Машиностроение, 1968. 283 с.
Феноменологический подход к определению показателя напряженного состояния для диаграммы пластичности
к.т.н. доц. Грушко А.В. Винницкий национальный технический университет, Украина +38-0432-59-84-65, [email protected]
Аннотация. Предложен параметр напряженного состояния в виде суммы относительных главных напряжений с коэффициентами их влияния. Коэффициенты влияния являются механическими свойствами материала. Диаграмму пластичности при помощи данного показателя предлагается аппроксимировать экспоненциальной функцией. Проверка для различных материалов показала достаточную близость опытных точек к диаграмме пластичности. Даны рекомендации по определению коэффициентов при произвольных видах испытаний. Минимальное количество испытаний - 4.
Ключевые слова: параметр напряженного состояния, диаграмма пластичности, механические свойства материала, разрушение
Введение
Пластичность металлов зависит от ряда факторов, среди которых, кроме природы самого материала, основными являются термомеханические параметры процесса: вид напряженного состояния, температура, история деформирования и пр. Наибольшее распространение получил показатель вида напряженного состояния, характеризующий пластичность, определяется показателем жесткости напряженного состояния (относительным гидростатическим
(71 + а 2 + 7 3
давлением) по В.А. Бабичкову в виде Т] = е-, где 7^ - интенсивность нормаль-
7i
ных напряжений (напряжение Мизеса), с - нормировочный коэффициент; с],72> а3 - главные нормальные напряжения [1-3].
Считается, что пластичность металла достаточно хорошо описывается показателем ^ при его отрицательных значениях. Однако в положительной области для пластичных материалов наблюдается немонотонность зависимости пластичности от ^, что вносит неопределенность при аппроксимации и экстраполяции экспериментальных данных. Поиск удобного показателя и функции диаграммы пластичности, отвечающие условиям монотонности, нулевой пластичности при равномерном трехосном растяжении и удовлетворительной корреляции с экспериментальными данными до сих пор является дискуссионным вопросом. Не исключено, что универсальный показатель напряженного состояния, инвариантный по отношению к механическим свойствам материала, может вообще не существовать, и для каждого материала, вообще говоря, может быть свой параметр, от которого монотонно зависит его пластичность [3, 4].
Целью работы является обоснование показателя вида напряженного состояния с учетом свойств материала при аппроксимации диаграммы пластичности монотонной функцией, поиск его констант и проверка для различных материалов и критериев деформируемости.
Исследование. Результаты
Основываясь на феноменологическом подходе, предположим, что пластичность зависит от величины
к1а1 + к?&2 + кзаз т г
-^^-а1 + к2а2 + кзаз , (1)
а1 а2 а3 7 7 7
где а1 - —,а2 --,аз -- - относительные главные напряжения, к1,к2,кз - коэффициенты влияния относительных главных напряжений на пластичность металла - являются механическими характеристиками материала, могут быть как положительными так и отрицательными, а также равными нулю. В дальнейшем показатель (1) нормирован к едини_ _ к о
це по величине к1 -1, тогда к2 - —, кз - — , т. е. при одноосном растяжении С1 - к1 - 1
к1 к1
Если к1 - к2 - кз -1 - влияние всех относительных главных напряжений равносильно, то получим С - ^ . Отметим также, что практически все известные показатели можно выразить через величины а1, а2, аз с учетом к},к2,кз что дает возможность ставить в соответствие различные показатели через соответствующие формулы перерасчета [1-5, 8]. Однако однозначное соответствие корректно лишь при плоском напряженном состоянии или частных случаях осесимметричного или плоского деформированного состояний.
Удачно подобранная аппроксимирующая функция дает возможность строить последнюю по результатам испытаний с минимальным числом опытов и возможностью экстраполяции.
Отметим удобную для практических расчетов аппроксимацию диаграммы пластичности по В.А. Огородникову [1, 2]
ер - ер0ехР(, (2)
где ер0 - пластичность металла при сдвиге,
X - чувствительность пластичности металла к изменению схемы напряженного состояния.
В областях положительных и отрицательных гидростатических давлений предполагается различная чувствительность пластичности металла X. Функция (2) дает необходимую точность при 0 и сильно зависящую от вида металла при г/> 0 [8].
В результате анализа ряда литературных данных [5-7] (алюминий, сталь), а также проведенных экспериментов для сталей У8, У12, титанового сплава ВТ6 установлено, что диаграмма пластичности и для показателя С достаточно точно описывается экспоненциальной функцией вр - й • ехр(— q • С);
вр - й • ехр (— q • ( + к2а2 + кзаз )),
где - константы диаграммы пластичности, определяемой показателем С .
Выбор такой функции связан также с тем, что при трехосном растяжении пластичность равна нулю, при трехосном сжатии - бесконечности. Очевидно, что должно выполняться условие: при q > 0, 1 + к2 + кз > 0; при q < 0, 1 + к2 + кз < 0, иначе при трехосном растяжении пластичность металла будет равна бесконечности.
Если С-1, то выражение й• ехр(— q) определяет пластичность при растяжении вр1;
если положить к2 - кз - 1 и обозначить q - X, й - вр0 то получим зависимость (2). Тогда, по аналогии с X, q есть чувствительность пластичности металла к изменению схемы напря-
женного состояния, определяемой показателем С. Когда к 2 = кз Ф 1, коэффициенты q , d теряют ясный физический смысл. В частности, параметр d есть условная пластичность металла, когда комплекс <1 + к2<2 + кз<з равен нулю (или С = 0).
Аппроксимация (3) содержит 4 неизвестных величины, следовательно, необходимо по крайней мере 4 экспериментальные точки для определения диаграммы пластичности. Наибольший интерес представляет область положительного гидростатического давления ( а> 0), поскольку именно здесь наблюдается немонотонность зависимости пластичности от ^. Таким образом, ограничиваясь областью а> 0 для калибровки модели выделим 4 характерные точки - испытанием образцов на растяжение, сдвиг, плоское деформированное состояние и равномерное двухосное растяжение. Если количество экспериментальных точек для построения диаграммы пластичности превышает 4, то следует искать параметры модели (3) методом наименьших квадратов.
Проверку предложенной модели выполним на основе ряда литературных данных [5-8]. Так, на рисунке 1 показаны диаграммы пластичности алюминиевого сплава 6061-Т6 (лист 2 мм) [6], построенные в зависимости от показателей напряженного состояния ^ и С • Диаграмма пластичности в координатах ^ - ер имеет сильную нелинейность и характерные экстремальные зоны, а диаграмма пластичности в координатах С - ер имеет свойство гладкости и тесную корреляционную связь с экспериментом. Ее функция имеет вид:
ер = 2578 ■ ехр(- 8,124 ■ С) = ехр(7,85 - 8,124 ■ С), (4)
а1 -0,063 ■ а2 -0,761 ■а3
где С =
<1
Чистый сдвиг Одноосное растяжение Двухосное растяжение
+ а2 + а3
деформирования
Рисунок 1 - Диаграммы пластичности алюминиевого сплава 6061-Т6 (лист 2 мм)
Аналогичные результаты дала обработка экспериментов, изложенных в [5-8]. Для дальнейшей проверки исследовалась пластичность высокоуглеродистых сталей У 8 и У12 (состояние поставки) и титанового сплава ВТ6 (отжиг). Для получения расчетных точек ис-пытывались плоские образцы в условиях сдвига, с выкружками, в условиях плоского изгиба и на двухосное растяжение по пробе Эриксена. Параметры напряженного состояния и величины предельных деформаций рассчитывались по методу сеток. В таблице 1 приведены коэффициенты аппроксимации (3).
Неудобство использования параметра напряженного состояния С состоит в его уникальности для каждого материала. Можно предположить, что для определенной группы ме-
таллов некоторые константы, входящие в показатель, будут иметь определенные корреляционные связи. Однако данное предположение требует глубокой экспериментальной проверки.
Следует отметить высокую степень корреляции рассчитанных диаграмм с экспериментом, а также слабую зависимость коэффициента кз от свойств исследуемых материалов (-0,75...-0,82).
Таблица 1
Коэффициенты аппроксимации, согласно (3)_
Материал а Я к2 к3
алюминиевый сплав 2025-Т351 [5] 58130 11,81 -0,108 -0,786
алюминиевый сплав 6061-Т6 (лист 2 мм) [6] 2578 8,124 -0,063 -0,761
сталь 1№ВР [7] 202 -0,134 0,028 -0,822
титановый сплав ВТ6 38,4 5,09 0,2 -0,75
сталь У12 23,33 3,63 0,02 -0,79
сталь У8 11,24 3,153 0,08 -0,78
Использование предложенного показателя напряженного состояния С и соответствующих аппроксимаций (3) в моделях разрушения состоит в том, что вместо традиционной функции диаграммы пластичности вр (77) следует применять функцию диаграммы разрушения вр (с) , поскольку физический смысл накопления повреждений введенный показатель С
не изменяет. Область применения того или иного критерия остается в соответствие с рекомендациями его использования в зависимости от вида процесса и истории нагружения.
Таким образом, аналитические зависимости расчета степени использования запаси пластичности и деформационные критерии разрушения запишем с учетом вр (С). Например, критерий В. Л. Колмогорова [9]
*
¥- } —¡Ц < 1. (5)
0
Тензорный критерий Г. Д. Деля [4]
в
4)
ае? а — 1 + дЯ1 — а)2 + уууу 2аф)—^, у--*-У < 1. (6)
а I 2аф , , ч , у -
вр (тт) 2а
5р'
Направляющий тензор скоростей деформаций Ру также следует поставить в зависимость от параметра С
2
У
^ - ГУ С (7)
где Гу (С) - функция, связывающая историю нагружения в виде направления с параметром С . Г у (С), может быть рассчитана с использованием ЭВМ по напряженно-
деформированному состоянию.
Рассмотрим процесс, характеризующийся линейностью пути деформирования в координатах 7 - ер. В координатах С - ер для алюминиевого сплава 6061-Т6 (лист 2 мм) тот же
путь деформирования будет иметь достаточно сложную форму (см. рисунок 1). Для примера
расчета степени использования запаса пластичности (СИЗП) применим (5):
* *
¥-1^ - 0,87, - 1,1. (8)
0 вр т) 0 вр (7)
Погрешность расчета СИЗП составила А = 2б%, что существенно при оценке деформируемости, поскольку с использованием менее точной модели разрушение должно наступить, а по предложенной - имеется достаточный запас.
Предложенную аппроксимацию можно расширить до области одноосного сжатия, т.е. для области — i <q< 2, в которой реализуются большинство процессов штамповки. В этом случае диаграмму разрушения описывает одна аналитическая функция, что удобно при моделировании процессов разрушения. Следует отметить, что для проверки и калибровки модели использовались преимущественно данные для плоского и частично для осесимметрич-ного напряженного состояния. При объемном напряженном состоянии проверка требует постановки экспериментов в камере высокого давления [1, 3].
Выводы
Использование параметров напряженного состояния, с учетом свойств материала, позволяет аппроксимировать диаграммы пластичности монотонными функциями с существенно меньшим отклонением от опытных данных, чем с использованием параметра q. Одним из таких параметров может быть С - сумма относительных главных напряжений с весовыми коэффициентами их влияния на пластичность в виде (3). Больше всего влияет на пластичность ai = ai/ai , второе по силе влияния - aß = aß/ai (ki > |kß| > ). Установлено, что
коэффициент kß для исследованных алюминиевых сплавов равен -G,76...-0,79, углеродистых
сталей, титанового сплава ВТ6 -G,S2...-0,85. Для получения диаграммы разрушения необходимо провести 4 испытания для разных схем напряженного состояния. С учетом постоянства коэффициента kß - три. Точность диаграммы разрушения в области положительных растягивающих напряжений существенно выше, чем диаграммы пластичности ep 7—) при проведении такого же количества опытов. Использование диаграммы разрушения ep 7С) в известных деформационных критериях для процессов, протекающих в области положительных гидростатических давлений, может существенно уточнить значение СИЗП до 2-х раз.
Литература
1. Огородников В. А. Оценка деформируемости металлов при обработке давлением. - К. : Выща шк., 1983. - 175 с.
2. Огородников В. А. Энергия. Деформации. Разрушение (задачи автотехнической экспертизы). I В. А. Огородников , В. Б. Киселев , И. О. Сивак - Винница : Универсум-Винница, 2005. - 204 с.
3. Дель Г.Д. Технологическая механика. - М.: Машиностроение, 1978. - 174 с.
4. Dell, H. CrachFEM A Comprehensive Approach For The Prediction Of Sheet Metal Failure I Dell, H.; Gese, H. and Oberhofer, G. II NUMIFORM 'G7, Materials Processing and Design: Modeling, Simulation and Applications, Porto, 1S - 21 June 2GG7, American Institute of Physics, CP9GS, pp. 165-17G.
5. Wierzbicki, T. Calibration and evaluation of seven fracture models I Wierzbicki, T.; Bao, Y.; Lee, Y.-W.; Bai, Y. II International Journal of mechanical Sciences, 47 (2GG5) pp. 719-743.
6. Luo, M. Ductile Fracture Calibration and Validation of Anisotropic Aluminum Sheets I M. Luo, T. Wierzbicki, II By Massachusetts Institute of Technology Conference: 2GG9 SEM Annual Conference & Exposition on Experimental & Applied Mechanics Proceedings, June 1-4, 2GG9.
7. Kolleck, R. Failure prediction method for hydro forming simulation of thick walled tubes I Kolleck, Ralf; Auer, Peter; Auer, Gerfried II The 14th International Esaform Conference On Material Forming: ESAFORM 2G11. AIP Conference Proceedings, Volume 1353, 2G11. - рр. 295-3GG.
8. Грушко А. В. Параметр напряженного состояния, учитывающий свойства материала, и его влияние на пластичность I А. В. Грушко II Вюник Нацюнального техтчного утверситету Украши "КШ". Серк: "Машиностроение". - 2G12. - №64. - С. 220-22б. - ISSN 23G5-9GG1.
9. Колмогоров В.Л. Напряжения, деформации, разрушение. - М.: Металлургия, 1970. - 229 с.