Full Text

Введение Закритическая стадия деформирования материала характеризуется снижением уровня напряжений при растущих деформациях [1-8]. Процесс накопления повреждений в материале с развитием зон разупрочнения может быть устойчивым при достаточной жесткости нагружающей системы [9, 10]. Основные теоретические положения механики закритического деформирования представлены в [11, 12]. Экспериментально полная диаграмма деформирования материала может быть получена с использованием испытательных машин с высокой жесткостью [13-20]. В работе [21] получены опытные данные о закритическом поведении конструкционных сталей при различных видах напряженно-деформированного состояния, в частности в опытах на кручение и совместное растяжение с кручением образцов. Изучено влияние различных факторов на степень реализации стадии разупрочнения. Рассмотрено влияние повышенной температуры на проявление закритического поведения сталей [19, 21]. Исследовано влияние предварительной закритической деформации на процессы ползучести и релаксации при высоких температурах. Показано положительное влияние дополнительных вибрационных воздействий на устойчивость процессов закритического деформирования. Аналитические решения задач с учетом полной диаграммы деформирования материала получены в работах В.А. Ибрагимова, В.Д. Клюшникова, Л.В. Никитина, Е.И. Рыжака, С.Д. Волкова, Г.И. Дубровиной, Ю.П. Соковнина, В.В. Стружанова, В.Э. Вильдемана [11, 22-26]. Численные решения краевых задач механики закритического деформирования получены в работах В.П. Радченко, С.В. Горбунова, В.Э. Вильдемана и других авторов [27-30]. 1. Математическая постановка краевой задачи механики закритического деформирования Краевая задача механики закритического деформирования [11-12, 21] (без учета массовых сил) представляет собой систему из 15 дифференциальных уравнений: (1) где - радиус-вектор; - тензор напряжений; - тензор конечных деформаций; - вектор перемещений; - тензор упругих свойств; ; - дельта-символ Кронекера; - тензор поврежденности четвертого ранга. В случае изотропного материала физические соотношения можно записать в виде (2) где - скалярная функция поврежденности [31]; - интенсивность напряжений; - интенсивность деформаций; G - модуль сдвига. Данные уравнения дополняются граничными условиями контактного типа, которые учитывают жесткость (либо податливость) нагружающей системы: (3) где - коэффициенты жесткости нагружающей системы; - коэффициенты податливости нагружающей системы; - направляющие косинусы вектора нормали к поверхности тела в данной точке; - перемещения, номинально заданные на границе ; - силы, номинально заданные на границе . При бесконечно большой жесткости нагружающей системы данные граничные условия вырождаются в граничные условия в перемещениях, при бесконечно малой - в напряжениях. Все решения, приведенные ниже, получены в стандартном вычислительном комплексе Abaqus CAE. Рассмотрены модели материалов с изотропным упрочнением и упругой разгрузкой. Полные диаграммы деформирования задавались в виде таблицы зависимости эквивалентных напряжений по Мизесу от эквивалентных пластических деформаций; промежуточные значения получаются интерполированием точек. 2. Особенности экспериментального изучения закритической стадии деформирования Для экспериментального изучения закономерностей закритического деформирования материалов применяются испытательные машины, предназначенные для работы в режиме кинематического нагружения и характеризующиеся высокой жесткостью. Это требуется для реализации контролируемой скорости деформирования материала на стадии разупрочнения и обеспечения высокой жесткости нагружающей системы по отношению к испытываемому образцу. С той же целью используются образцы с малой длиной рабочей части. Были проведены испытания образцов различных конструкционных сталей на растяжение, в результате которых получены диаграммы деформирования с развитыми участками закритического поведения. На рис. 1 приведены диаграммы деформирования для сталей: 1 - 15Х2ГМФ, 2 - 40Х, 3 - 45, 4 - 25, 5 - Ст3, 6 - 20. Испытания проводились на универсальной двухосевой сервогидравлической испытательной системе Instron 8850, обладающей высокой жесткостью [21]. Данные на рис. 1 приведены в осях инженерных напряжений и инженерных деформаций, полученных путем отношения нагрузки и перемещений к начальной площади поперечного сечения и начальной длине рабочей части образца соответственно. Учет возниконовения локализации деформаций в виде «шейки» при интерпретации результатов испытаний в условиях реализации стадии закритического поведения является самостоятельным направлением расчетно-экспериментальных исследований и испытаний при совместном использовании машины и видеосистемы регистрации полей перемещений и деформаций, основанной на методе корреляции цифровых изображений. Рис. 1. Диаграммы деформирования сталей: 1 - 15Х2ГМФ; 2 - 40Х; 3 - 45; 4 - 25; 5 - Ст3, 6 - 20 Fig. 1. Strain curves of steels: 1 - 15CrMoV4, 2 - 40Cr4, 3 - C45E, 4 - C25E, 5 - A3, 6 - 20 Из рис. 1 видно, что закритическая стадия деформирования, начало которой соответствует достижению максимальных напряжений, является более протяженной по сравнению со стадией упругопластического деформирования. Очевидно, что учет этой стадии в расчетах должен обнаружить деформационный резерв элементов конструкций. Для дальнейших расчетов будет использована диаграмма 2 на рис. 1, полученная при испытании образца из стали 40Х. 3. Моделирование процесса деформирования пластины с концентратором напряжений для двухзвенной и трехзвенной диаграммы деформирования Рассмотрим моделирование процесса разрушения тонкой пластины с круглым концентратором при растяжении. Геометрия образца: ширина и длина пластины 100 мм, диаметр центрального отверстия 10 мм; к верхней грани приложено перемещение 0,1 мм за условную единицу времени, нижняя грань закреплена шарнирно. Принято, что пластина состоит из однородного материала, модуль Юнга Е = 200 ГПа, коэффициент Пуассона = 0,3, предел текучести = 200 МПа. Были рассмотрены следующие модели: упругопластическая с линейным упрочнением, модули упрочнения = 100 ГПа, 25 ГПа; идеальная упругопластическая модель; упругопластическая с линейным разупрочнением, модули спада D = 5 ГПа, 10 ГПа. Диаграммы деформирования приведены на рис. 2. Рис. 2. Диаграммы деформирования: 1 - = 100 ГПа; 2 - = 25 ГПа; 3 - = 0 ГПа; 4 - D = 5 ГПа; 5 - D = 10 ГПа Fig. 2. Strain curves: 1 - = 100 GPa, 2 - = 25 GPa, 3 - = 0 GPa, 4 - D = 5 GPa, 5 - D = 10 GPa В результате моделирования получим набор полей перемещений, деформаций и напряжений для разных моментов времени. На рис. 3 показано, каким образом меняется пластическая зона для разных моделей материала. Для материала с упрочнением наблюдается широкая пластическая зона, расположенная под углом к направлению нагружения. Рис. 3. Результаты расчета эволюции процесса деформирования пластины для диаграмм деформирования: а - с модулем упрочнения = 25 ГПа; б - с площадкой идеальной текучести = 0 ГПа; в - с модулем разупрочнения D = 10 ГПа Fig. 3. Calculation results of plate deformation evolution for strain curves: a - with hardening module = 25 GPa, b - with perfect flow area = 0 GPa, c - with softening module D = 10 GPa Для идеального упругопластического материала данная полоса становится более узкой, аналогичный результат получен для материала с модулем разупрочнения 5 ГПа. Для модуля спада D = 10 ГПа в некоторый момент пластическая зона начинает резко вытягиваться, приобретая вид трещины, происходит потеря устойчивости численного решения. Рассмотрим аналогичную задачу для пластины из упругопластического материала с участками линейного упрочнения ( = 50 ГПа) и разупрочнения (D = 60 ГПа и 120 ГПа). Диаграммы деформирования материала представлены на рис. 4. Рис. 4. Диаграммы деформирования: 1 - D = 60 ГПа; 2 - D = 120 ГПа Fig. 4. Strain curves: 1 - D = 60 GPa, 2 - D = 120 GPa На рис. 5 показано, как меняются пластическая зона и зона закритического деформирования для нагружения с модулем разупрочнения D = 60 ГПа (темно-серый цвет - зона упрочнения, черный цвет - зона разупрочнения), для разупрочнения с модулем D = 120 ГПа зоны аналогичны. При уменьшении модуля спада зона закритического деформирования становится более выраженной. При прохождении магистральной трещины в ее вершине находится небольшая зона ослабленного материала. Рис. 5. Результаты расчета эволюции процесса деформирования и разрушения пластины Fig. 5. Calculation results of plate deformation and destruction evolution Рассмотрим эпюру напряжений вдоль горизонтальной оси симметрии, изображенную на рис. 6. На ней можно выделить несколько характерных точек: 1 - точка, где полностью реализована закритическая стадия деформирования материала; 2 - точка, соответствующая частичной реализации закритической стадии; 3 - точка, соответствующая пределу прочности материала; 4 - точка в области линейного упрочнения; 5 - точка, соответствующая границе пластической зоны (предел текучести); 6 - точка в упругой области. Из данного рисунка видно, каким образом реализуется диаграмма деформирования с учетом закритической стадии при данном нагружении. Рис. 6. Эпюра напряжений σ22 с соответствующими точками на диаграмме деформирования материала Fig. 6. Stress diagram σ22 with corresponding points on the material deformation diagram Отмечено, что с уменьшением модуля разупрочнения происходит увеличение максимальной расчетной нагрузки, выдерживаемой конструкцией, система позже теряет устойчивость. Конструкция способна к увеличению прикладываемой внешней нагрузки даже при развитии зон закритического деформирования. Для диаграммы деформирования 1 увеличение нагрузки составило 2,1 % относительно расчета без учета закритической стадии деформирования. Следовательно, ее учет рационален для уточнения несущей способности конструкций. 4. Моделирование процесса деформирования пластины с концентратором напряжений при использовании реальных диаграмм деформирования Рассмотрим процесс деформирования образца (геометрия изображена на рис. 7, толщина пластины 3 мм) из стали 20. Диаграмма деформирования материала приведена на рис. 8, на ней обозначены точки, по которым строится кусочно-линейная аппроксимация диаграммы деформирования, которая будет использована в численном анализе. Было принято, что по достижении последней точки из указанных на диаграмме он продолжает нести маленькую нагрузку. К торцу (аналогично предыдущим задачам) прикладывалось перемещение 10 мм за условную единицу времени. Эволюция процесса деформирования пластины изображена на рис. 9. Рис. 7. Геометрия образца Fig. 7. The sample’s geometry Рис. 8. Диаграмма деформирования стали 20 Fig. 8. Steel 20 strain curve Рис. 9. Результаты расчета эволюции процесса деформирования пластины из стали 20 (светло-серый цвет - зона с начальным упрочнением, темно-серый цвет - зона с развитым упрочнением, черный цвет - зона разупрочнения) Fig. 9. Calculation results of a 20 steel plate deformation evolution (the light gray color shows the zone with initial hardening, the dark gray color shows the zone with developed hardening, the black color shows the zone of softening) Построена расчетная диаграмма нагружения (рис. 10), она качественно соответствует экспериментальной [21]. Видно, что учет закритической стадии деформирования в данном случае позволяет получить выигрыш 2 % по нагрузке и 20 % по перемещению. Рассмотрим процесс деформирования прямоугольной пластины 10×20 мм с центральным эллиптическим концентратором (большая и малая полуоси 2,5 и 1,25 мм соответственно) из стали 40Х. Диаграмма деформирования материала приведена на рис. 11, на ней обозначены точки, по которым строится кусочно-линейная аппроксимация диаграммы деформирования, которая будет использована в численном анализе. а б Рис. 10. а - Расчетная диаграмма нагружения: 1 - появление зоны с малым модулем упрочнения; 2 - появление зоны с большим модулем упрочнения; 3 - появление зоны закритического деформирования; б - экспериментальная диаграмма нагружения Fig. 10. Calculated loading diagram: 1) the appearance of a zone with a small hardening modulus; 2) the appearance of a zone with a large hardening modulus; 3) the appearance of a zone of postcritical deformation; b the experimental loading diagram На диаграмме деформирования можно выделить следующие участки: упругий, упрочнения, начального разупрочнения и развитого разупрочнения. Рис. 11. Диаграмма деформирования стали 40Х; отмеченные точки заложены в модель Fig. 11. Steel 40Cr4 deformation diagram; the marked points are embedded in the model Результаты расчета эволюции процесса деформирования пластины изображены на рис. 12. Видно, как по мере деформирования появляется зона разупрочнения вблизи вершины концентратора. На расчетной диаграмме нагружения (рис. 13) видно, что рост нагрузки не останавливается после появления зоны разупрочнения; выигрыш по нагрузке составил 4,8 %, деформационный ресурс - 50,0 %. После достижения максимального значения нагрузки на образце наблюдается локализация деформаций в узкую полосу. На рис. 14, 15 изображены поля напряжений σ11, σ22, σ12, σi вблизи вершины концентратора напряжений для момента максимальной расчетной нагрузки (точка в на рис. 13) и момента потери численного решения (точка г на рис. 13): видно, что по мере приближения к вершине концентратора происходит уменьшение эквивалентных напряжений в связи с реализацией закритической стадии деформирования. Аналогичные результаты были получены в работе [27]. Рис. 12. Результаты расчета эволюции процесса деформирования пластины с эллиптическим концентратором из стали 40Х (светло-серый цвет - упругая область, серый цвет - зона упрочнения, темно-серый цвет - зона разупрочнения) Fig. 12. Calculation results of the 40Cr4 steel plate with the elliptical concentrator deformation evolution (the light gray color shows the elastic area, the gray color shoes the hardening zone, the dark gray color shows the softening zone) Рис. 13. Расчетная диаграмма нагружения: а - появление зоны упрочнения; б - появление зоны разупрочнения; в - точка с максимальной нагрузкой; г - точка срыва Fig. 13. The calculated loading diagram: a) the appearance of a hardening zone; b) the appearance of a zone of softening; c) the point with a maximum load, d) the point of failure Аналогичные задачи были решены для концентраторов напряжений иной геометрии: эллиптического, круглого, узкого прямоугольного, широкого прямоугольного, ромбовидного. Зоны закритического деформирования для максимального значения расчетной нагрузки приведены на рис. 16. Значения выигрышей по нагрузке и перемещению приведены в таблице. Видно, что учет закритической стадии деформирования позволяет выявить значительные прочностные и деформационные ресурсы конструкций с концентраторами напряжений. Рис. 14. Поля напряжений у вершины концентратора для точки с максимальной нагрузкой: a - σ11; б - σ22; в - σ12; г - σi Fig. 14. Fields of stresses near the top of the concentrator for the point with the maximum load: a - σ11; b - σ22; c - σ12; d - σi Рис. 15. Поля напряжений у вершины концентратора для точки срыва: a - σ11; б - σ22; в - σ12; г - σi Fig. 15. Fields of stresses near the top of the concentrator for the point of failure: a - σ11; b - σ22; c - σ12; d - σi Рис. 16. Зона разупрочнения для разных концентраторов: 1, 2 - эллиптические; 3 - круглый; 4 - узкий прямоугольный; 5 - широкий прямоугольный; 6 - ромбовидный Fig. 16. The softening zone for different concentrators: 1, 2) elliptical; 3) round; 4) narrow rectangular; 5) wide rectangular; 6) diamond-shaped Выигрыши по нагрузке и перемещению для концентраторов разной геометрии Номер концентратора Выигрыш по нагрузке, % Выигрыш по перемещению, % 1 4,8 50,0 2 10,8 75,8 3 7,6 59,4 4 12,1 92,3 5 10,7 75,3 6 11,0 85,0 Заключение В работе получены новые численные решения задач о деформировании пластин с концентраторами напряжений различной геометрии. Рассмотрены двухзвенная и трехзвенная аппроксимации диаграммы деформирования, а также реальные диаграммы деформирования стали 20 и стали 40Х. Изучена эволюция процесса деформирования пластины с выявлением зон закритического деформирования материала. Построены расчетные диаграммы нагружения. Выявлены значительные прочностные и деформационные ресурсы при учете закритической стадии деформирования материала. Из всего вышеизложенного можно сделать вывод, что учет стадии разупрочнения материала целесообразен при проведении расчетов конструкций с концентраторами напряжений различной геометрии.

About the authors

V. E Wildemann

Perm National Research Polytechnic University

M. P Tretyakov

Perm National Research Polytechnic University

A. I Mugaratov

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 462

PDF (Russian) - 536