Full Text

Введение Создание изделий для современной техники с повышенным уровнем механических и служебных свойств инициирует обращение к особому классу связанных задач механики деформируемого твердого тела, ориентированных на описание процессов формирования существенно неоднородной структуры при высокоэнергетических термомеханических воздействиях [1]. Актуальность данной проблемы подтверждается многочисленными публикациями с анализом ее различных аспектов в области сварки [2], литья [3], термической обработки [4], поверхностного упрочнения [5] и т. д. Однако решение таких задач с единых позиций осложнено необходимостью учета непрерывных и быстротекущих трансформаций структуры, свойств, высоко градиентных напряженно-деформированных состояний (НДС), эффектами динамики (с анализом необходимости учета инерционных членов при записи разрешающих уравнений) и связанности процессов, сложностью согласования временных сеток на различных этапах и масштабах решения тепловой, механической и структурной задач. Кроме того, требуется обращение к специфическим соотношениям вычислительного материаловедения при установлении связей между термонапряженным состоянием в точках тела и образующейся в них структурой. В работе рассматриваются особенности единой постановки и решения связанных задач на примере электромеханической обработки (ЭМО) стали. 1. Постановка задачи Решение указанной связанной технологической задачи определения НДС неоднородного тела базируется на совместном рассмотрении нескольких групп основных уравнений, соответствующих трем базовым моделям системы. Первую группу составляют уравнения движения и гиперболическое уравнение теплопроводности, соответствующие задаче термоупругопластичности при интенсивном термосиловом нагружении: (1) (2) где T - текущая температура; t - время; kT - коэффициент теплопроводности; c - коэффициент теплоемкости; r - плотность; trel - время релаксации теплового потока; ; f - мощность внутренних источников тепла; l, m - параметры Ламе; ; ekk - первый инвариант тензора деформаций (относительное изменение объема); ; ; L - энтальпия фазового перехода; Y - массовая доля новой фазы при фазовом переходе; ui - компоненты перемещения; - оператор Лапласа; a - коэффициент теплового расширения; DT = T-T0 - изменение температуры тела относительно начальной величины T0; D -относительное изменение размеров при структурном превращении; - дополнительные напряжения (по методу дополнительных деформаций). Для компактности записи уравнений (1), (2) запятая перед нижним индексом означает частную производную по соответствующей переменной. Особенность формирования системы связанных технологических задач заключается в необходимости анализа диаграмм состояния структур металлического сплава, подвергнутого ЭМО, для определения физико-механических характеристик фаз в данной точке в данный момент времени. Для этого система (1), (2) дополняется условиями формирования металлографических структур стали на основе существующих зависимостей [1, 6] критических температур Tкр и скоростей их изменения от текущего химического состава стали, определяющих форму строения металла в каждой точке тела: (3) где Tmin, tmin - температура и время минимальной устойчивости аустенита. Конкретизация зависимостей (3) для сталей различных классов дается в специальной материаловедческой литературе [1, 6-8], например: для температуры Ac3 конца аустенитного превращения (в °C) [7] для температуры Mн начала мартенситного превращения [8] и т.д. Кроме членов, учитывающих взаимное влияние термодеформационных процессов ( - в (1) и - в (2)), система (1), (2) включает компоненты, определяющие величину дополнительных деформаций ( ) и скрытую теплоту при фазовых переходах. При этом условия (3) дополняются специфическими факторами связанности, устанавливающими смещение критических температур стали в зависимости от скорости нагрева и компонент sij тензора напряжений и соответственно). (4) Как и в случае (3), соотношения (4) задаются для каждой из критических точек данной марки стали на стадии нагрева и охлаждения, например: для точки Ac1 углеродистой стали, согласно [9] где DTн - изменение температуры начала превращения, °С; a0 - структурный параметр, равный половине среднестатистического расстояния между зародышами аустенита, см; D - коэффициент диффузии углерода в аустените, см2/с; - константа материала; K - коэффициент, имеющий размерность температуры (для углеродистых сталей K = 230 °С); для мартенситной точки DMн [6, 10] где so - октаэдрические нормальные напряжения; J2 - второй инвариант тензора девиатора напряжений; A, B - эмпирические коэффициенты. На рис. 1 показана схема взаимодействия связанных задач механики, теплофизики и материаловедения в системе моделей неоднородного тела в условиях нестационарных термосиловых и структурных полей, применительно к технологии поверхностного упрочнения материалов ЭМО. Выделены три группы базовых моделей, соответствующих системе разрешающих соотношений (1)-(3), а также группа моделей текущего состояния тела (МТС). Последняя включает модели конфигурации, физико-механических свойств формирующегося неоднородного тела и оценки его НДС (например, модели изменения геометрии образующихся структурных зон и внешних контактных поверхностей с уточнением граничных условий связанных задач в ходе деформации; изменение тепловых и механических характеристик материала при фазовых превращениях и т.д.). В силу обширности, разнородности и частного характера моделей описание группы МТС вынесено за рамки данной статьи. Отмечены факторы связанности рассматриваемых задач, задаваемые соответствующими членами основных уравнений. Для базовых моделей показаны дополнительные модули: 1) согласования расчетных временных шагов (Dt', Dt'', Dt''') при анализе рассматриваемых разноскоростных процессов; 2) учета динамических эффектов, связанных с инерционными членами в уравнениях (1)-(3) (включая изменение tmin в локальных областях неоднородности аустенита, полученного с быстрого нагрева). Рис. 1. Система связанных задач формирования структуры и свойств материалов при высокоэнергетических технологических воздействиях Начальные и граничные условия для системы (1), (2) задаются применительно к импульсной (переменным током) ЭМО стальной детали твердосплавным инструментом в форме ролика [11], моделируемой расчетной областью в виде полупространства с подобластями (зонами) с различающимися физико-механическими характеристиками (рис. 2). Границы данных зон соответствуют областям металла с разной структурой, формирующейся в ходе термосилового нагружения (например, при ЭМО - фрагментам со структурой закалки, отпуска и исходного состояния; зоны 1-3 на рис. 2). Их текущее положение устанавливается на каждом расчетном шаге из решения рассматриваемой связанной задачи (1)-(4). Принимается гипотеза о недеформированном, ненагруженном состоянии тела в начальный момент времени: напряжения и деформации отсутствуют; температуры одинаковы по объему тела и равны начальной температуре технологической системы. Начальные скорости изменения температуры, с учетом синусоидального закона изменения модуля q вектора плотности теплового потока при нагреве металла переменным электрическим током, равны нулю. Скорости деформации в начальный момент времени также равны нулю в случае обработки материала инструментом в виде свободно вращающегося ролика, движущегося по поверхности вдоль координатной оси 0x с линейной скоростью v: (5) Рис. 2. Расчетная схема неоднородного тела и граничные условия при контактном температурно-силовом воздействии (в ходе ЭМО стали) В случае невращающегося или приторможенного (проскальзывающего) инструмента в точках его контакта с деталью должны быть заданы соответствующие начальные скорости вдоль направления движения (например, ). На внешней границе S устанавливаются поверхностные нагрузки, перемещения и тепловые потоки, моделирующие условия нагружения и контакта материала (с инструментом и оснасткой) при его обработке. При этом в подвижной области контакта Sк обеспечивается соответствие (см. рис. 2): 1) нормальных к поверхности напряжений - действующему давлению p инструмента на обрабатываемую деталь (остальные напряжения считаются равными нулю); 2) перемещений - профилю инструмента (считая его в первом приближении абсолютно жестким) в точках соприкосновения рабочих поверхностей; 3) плотности теплового потока - удельной мощности импульсного источника q, действующего в области Sк при пропускании электрического тока через зону контакта «инструмент-деталь» в ходе ЭМО (6). Необходимо отметить, что конфигурация зоны контакта Sк заранее неизвестна и устанавливается из итерационной процедуры на основе альтернирующего метода Шварца [12, 13] при поочередном задании на поверхностях контакта векторов перемещений и поверхностных сил с учетом движения зоны контакта и термоструктурных деформаций. В этом случае область Sк является фактором кинематической связанности задач, от которого зависят плотность теплового потока и давление в зоне контакта, дополняя схему физической связанности (см. рис. 1). За пределами зоны контакта Sк наружная поверхность тела охлаждается за счет теплоотдачи во внешнюю среду (с температурой Tср) и остается свободной от напряжений. На границе m и (m + 1)-й структурных областей задаются условия совместности перемещений и поперечных (к поверхности раздела заданной в локальной системе координат abg с нормалью g в данной точке (см. рис. 2) нормальных и касательных напряжений, на бесконечности - равенство перемещений и напряжений нулю. Перечисленные граничные условия записаны в виде системы (6) Физические соотношения записываются в рамках метода дополнительных деформаций с использованием соотношений теории пластического течения, в предположении изотропного упрочнения, отсутствия эффекта Баушингера и существования единой обобщенной многомерной поверхности неизотермического деформирования металла с трансформирующейся структурой [14, 15], представленной на рис. 3, а набором поверхностей для различных Y. Приращение полной деформации представляется в виде суммы приращений упругих , температурных , структурных и дополнительных деформаций , описывающих поведение при пластичности в условиях изменяющихся напряжений, температуры и структуры металла: (7) Приращения упругих деформаций определяются законом Гука: (8) для остальных параметров используются следующие соотношения: (9) где si - интенсивность напряжений; s - среднее напряжение; sij - компоненты девиатора напряжений; E = E (T, Y) - модуль упругости; Eк = Eк (si, T, Y) - касательный модуль; sт = sт (eip, T, Y) - мгновенный предел текучести, соответствующий интенсивности накопленной пластической деформации eip, температуре T и содержанию Y исследуемой фазы; , - изменение мгновенного предела текучести с изменением температуры и фазовых превращений материала; - символ Кронекера; ak - коэффициент теплового расширения; Dk - удельное изменение объема k-й фазы при превращении. а б в Рис. 3. Анализ условий приращения пластических деформаций на основе семейства поверхностей пластического деформирования (A' - исходное состояние материала на данном этапе нагружения, A'' - после изменения температуры при неизменном структурном состоянии, A''' - после изменения температуры и структуры материала): а - общая схема; б - при нагреве (B' - нагружение, B'', B''' - разгрузка); в - при охлаждении (B'', B''' - нагружение, B' - разгрузка) Зависимости диаграмм деформирования и механических свойств (E, Eк, sт) сталей в различных структурных состояниях от температуры принимаются на основе эмпирических данных и их аппроксимаций, приведенных в работах [4, 16, 17]. Например, согласно [16] для модуля Юнга и предела текучести аустенита (МПа), Для определения доли новых образующихся фаз Y также используются эмпирические зависимости, обобщающие на различные структурные состояния уравнение Колмогорова-Джонсона-Метла-Аврами [16]. Неизотермический характер исследуемых процессов учитывается на основе правила аддитивности Шайля [18]. 2. Процедура решения задачи Алгоритм анализа системы рассматриваемых задач в целом соответствует схеме их взаимодействия (см. рис. 1). Система уравнений (1)-(9) решается методом конечных разностей на основе равномерной по каждому направлению разностной сетки с аппроксимацией второго порядка. Применяются как явная (для быстротекущих процессов в микромасштабных областях), так и неявная (для макромасштабных процессов) двух- и трехслойные конечно-разностные схемы, обоснование и сравнительный анализ которых частично (для тепловой задачи) рассмотрены в [19]. При решении используются итерационные процедуры переменных направлений. Устойчивость указанных схем проверяется по известным условиям [20], а также по результатам серии вычислительных экспериментов, приведенным ниже. Основной особенностью рассматриваемого класса задач является необходимость описания на каждом временном шаге нового состояния неоднородного тела, вызванного перестройкой его структуры в ходе термо-силового нагружения. В этом случае из совместного анализа зависимостей типа (3), (4) и действующих полей напряжений и температур, найденных из (1), (2), устанавливаются текущие положения границ областей структурной неоднородности. Для каждой области из банка данных задаются параметры диаграммы деформирования, теплофизические характеристики, а также возникающие структурные деформации, соответствующие изменению удельного объема металла при фазовых переходах. Далее, согласно (9) вычисляются приращения деформаций , , для текущего приближения. При расчете процедура метода дополнительных деформаций обобщается на случай семейства поверхностей неизотермического пластического деформирования при различных структурных состояниях (см. рис. 3). В зависимости от соотношения между величинами приращений мгновенного предела текучести dW1, dW2 и интенсивности напряжений dsi в ходе последовательных приближений устанавливается положение точки, изображающей траекторию нагружения и разгрузки материала на множестве поверхностей деформирования, соответствующих текущему сочетанию НДС, температуры и структурных параметров (см. рис. 3). Для установленного таким образом механического состояния в данной точке материала вычисляется приращение пластической деформации (согласно (9)) относительно предыдущего расчетного шага. Нагружению материала может соответствовать как увеличение интенсивности напряжений - в случаях, если dsi больше, чем приращение мгновенного предела текучести за счет упрочнения материала при его охлаждении или фазовом переходе (точка B' на рис. 3, а, б), так и их снижение, когда изменение свойств металла приводит к уменьшению мгновенного предела текучести на величину, превышающую dsi (точка B'' на рис. 3, в). Разгрузка соответствует активному снижению напряжений на величину dsi большую, чем падение мгновенного предела текучести при разупрочнении (точка B' на рис. 3, в). Повышение механических характеристик материала также может сопровождаться разгрузкой, если приращение интенсивности напряжений меньше приращения предела текучести (точка B'' на рис. 3, б). Для повышения точности каждый из этапов нагружения разбивается на достаточно малые расчетные шаги. Это относится как к постепенному приложению действующих активных нагрузок, так и к дифференциации шага временной сетки. При решении связанной задачи возникает необходимость согласования вычислительных процедур с различными шагами по времени в рамках различных методов, используемых при анализе рассматриваемых разноскоростных процессов (например, полный цикл обработки ЭМО поверхности цилиндрической детали, соответствующий одному ее обороту в шпинделе станка, составляет около 4,5 с; длительность одного электрического импульса при обработке переменным током - 0,01 с; стадия нагрева (нагружения) до максимальной температуры - 0,0075 с; время выдержки на мартенсит закалки - 0,5 с, на сорбит - 3,5 с и т. д.). Базовый шаг вычислений, определяемый из условий работоспособности численной схемы (устойчивости, допустимой погрешности счета, накопленной ошибки и т.д.), может уточняться в зависимости от требуемой точности для различных приближений - термического, механического, структурного (рис. 4). Рис. 4. Схема коррекции шагов временной сетки при совместном анализе разноскоростных процессов (термоструктурная задача) в области критических значений по термокинетической диаграмме ( , - уточненные значения температуры на i+1 и i+2 итерациях; s0,2 - условный предел текучести материала в текущем структурном состоянии [21-23]) На рис. 4 для примера показана схема коррекции расчетной траектории охлаждения (от прямой oa при базовом шаге по времени Δt0 до ломаной obde с уточненными шагами Δti, Δti+1 и т.д.) в момент структурного (мартенситного) превращения в стали для уточнения расчетных величин температуры (Ti+1 и Ti+2 в i+1 и i+2 итерациях) относительно целевого значения Ts в критической точке (определяемого по термокинетической диаграмме распада аустенита) (см. рис. 4). В этом случае скорректированная величина шага временной сетки для ближайших итераций находится в зависимости от заданной погрешности размеров структурных зон и условия равенства текущей расчетной скорости охлаждения и критической скорости закалки для рассматриваемого структурного превращения в данный момент времени: где - критическая скорость изменения температуры на i итерации; (N - число дополнительных временных шагов, принятых для уточнения решения; k - номер текущего шага); grad Ti - модуль температурного градиента; Δhmin - минимальный размер конечно-разностной ячейки; ξ - допустимая относительная погрешность определения размера структурной зоны , в первом приближении ξ = 0,5). Таким образом, расчетный временной шаг, уточненный на каждой итерации решения связанной задачи, обеспечивает рациональное (по вычислительным затратам) описание высокоскоростных эффектов вблизи характерных точек при контролируемом уровне погрешности (траектория ob-bd-de-ea, рис. 4). При этом расчет с фиксированным (для одной из задач) шагом может приводить к заметным погрешностям при анализе сопряженных процессов (например, на рис. 4 траектории oa при шаге Δt0 или ob-bc-ca при Δti, приводящие к ошибкам в определении размеров и составе формирующихся структурных зон). 3. Сходимость и устойчивость метода. Адекватность модели Устойчивость решения по Адамару проверялась на основе известной процедуры многократного пересчета с вариацией в широком диапазоне параметров сетки, граничных условий, свойств материала с установлением степени отклонения полученных результатов (табл. 1 и рис. 5). Как видно из табл. 1, при изменении количества узлов сетки в 3-8 раз (с соответствующим измельчением или увеличением шагов) значения максимальной температуры изменяются не более чем на 5 %, модуль температурного градиента - не более чем на 2 %, величины максимальных линейных деформаций (по оси 0z на рис. 2) - на 30-40 %, интенсивности напряжений - 15-20 % и т. д. Увеличение удельной мощности теплового источника q от 5 до 50 % сопровождается пропорциональным повышением температуры от 3 до 34 %, модуля градиента - от 5 до 50 %, деформаций - от 8 до 80 %, интенсивности напряжений - от 4 до 45 %. Указанные величины отклонений являются максимальными из достигнутых в различных точках расчетной области, при этом качественный характер полученных закономерностей не изменяется во всем диапазоне изменения параметров задачи (см. рис. 5). Таблица 1 Изменение максимальных расчетных значений температуры (T), модуля температурного градиента (grad T), осевой деформации (ε33) и интенсивности напряжений (σi) в зоне обработки при ЭМО в зависимости от количества узлов сетки и удельной мощности (q) теплового источника № Параметры решения T grad T ε33 σi 1 100x100x50 = 500000 узлов q = 5,82·108 Вт/м2 1720 °C 4284 °C/м 0,285 % 1888 МПа 2 50x50x25 = 62500 узлов q = 5,82·108 Вт/м2 +3,8 % -1,1 % -36 % -19 % 3 125x125x100 = 1562500 узлов q = 5,82·108 Вт/м2 -2,3 % +0,2 % +29 % +17 % 4 100x100x50 = 500000 узлов q = 6,11·108 Вт/м2 (+5 %) +3,0 % +5,0 % +7,8 % +4,3 % 5 100x100x50 = 500000 узлов q = 7,28·108 Вт/м2 (+25 %) +15,8 % +25,3 % +39,6 % +21,9 % 6 100x100x50 = 500000 узлов q = 8,73·108 Вт/м2 (+50 %) +33,5 % +50,9 % +81,9 % +45,3 % а б Рис. 5. Изменение расчетных максимальных линейных деформаций (вдоль оси 0z образца) и интенсивности напряжений по глубине поверхностного слоя стального образца при ЭМО для различных параметров сетки и мощности источника (обозначения кривых 1-6 соответствуют данным табл. 1) Адекватность разработанных моделей проверялась на основе сопоставления результатов решения тестовых задач при помощи различных методов. В частности, ниже приводятся данные такого анализа для задачи о напряжениях в бесконечном однородном цилиндре при термоупругих деформациях (рис. 6, а) и в неоднородном цилиндре с упрочненным поверхностным слоем при термоупругопластических деформациях (рис. 6, б). В первом случае (см. рис. 6, а) результаты, полученные по разработанной модели, сопоставляются с известным [24] аналитическим решением задачи и численным (МКЭ) решением на базе Ansys, во втором (см. рис. 6, б) - сопоставление проводится с результатами расчетов по численно-аналитической модели, приведенной в [25]. Взаимное отклонение результатов решения указанных задач по различным методикам не превышает 5 % (рис. 6). а б Рис. 6. Сопоставление расчетных напряжений, полученных разными методами: а - цилиндр с линейным по радиусу распределением температуры (линия - аналитический расчет, кружки - МКЭ Ansys, квадраты - рассматриваемая модель); б - полупространство с равномерным распределением температуры в эллиптической области (линия - численно-аналитический расчет [25], точки - рассматриваемая модель) 4. Влияние связанности и инерции Для оценки необходимости учета эффектов связанности при записи и анализе основных уравнений термоупругости (1), (2) проведены тестовые расчеты по определению скоростей изменения температуры и деформаций (рис. 7, а) в исследуемой расчетной области (неоднородное полупространство) при характерном для ЭМО импульсном температурно-силовом воздействии (при действии на поверхность материала движущегося эллипсоидного штампа и импульсного теплового источника). Показаны приведенные к максимальным значениям величины ( ) температуры (i = 1; Tmax = 1640 °C), скорости нагрева/охлаждения (i = 2; °C/c), объемной деформации (i = 3; ) и скорости ее изменения (i = 4; c-1) в фиксированной точке на поверхности тела. При рассматриваемых воздействиях указанные величины изменяются по подобным закономерностям с некоторым опережением во времени процессов деформации над процессами изменения температуры за счет механического действия движущегося штампа (рис. 7, а). На основе полученных данных вычисляются параметры связанности, определяющие влияние тепловых источников при деформации на распределение температуры и соотношение скоростей деформаций и температур. Согласно [24, 26] связанность будет малой, если указанные безразмерные параметры (обозначенные как ω и υ) удовлетворяют следующим соотношениям , (10) (11) На рис. 7, б показан характер изменения ω и υ в различные моменты времени протекания характерного для ЭМО термосилового нагружения (при действии на поверхность материала движущегося эллипсоидного штампа и импульсного теплового источника). Для установления диапазона возможных значений параметров в (10), (11) принято равенство начальной температуры T0 текущему значению T в данный момент времени согласно рис. 7, а. Рассчитанные для различных T максимальные величины параметра ω (соответствующие физико-механическим характеристикам стали) не превышают 0,08 (кривая 1, рис. 7, б), что удовлетворяет условию (7). а б Рис. 7. Изменение относительных расчетных величин (а) температуры (i = 1), скорости нагрева / охлаждения (i = 2), объемной деформации (i = 3) и скорости ее изменения (i = 4), а также (б) безразмерных параметров ω (кривая 1) и υ (кривая 2) в центре пятна контакта за время одного теплового импульса (полупериод электрического тока) Параметр υ (кривая 2, рис. 7, б) принимает значения от 1,5 до 0,5 на интервалах наибольших скоростей деформации и температуры (0-6 и 8-10 мс; рис. 7, б), резко увеличиваясь при околонулевых скоростях теплового поля (в интервале 6-8 мс; рис. 7, б). При этом правая часть неравенства (11) за время одного температурного импульса меняется от 55 до 8, превышая текущие значения υ на большей части исследуемого интервала. Это означает, что скорости изменения температур и деформаций являются величинами одного порядка и дают возможность с учетом малости дополнительных слагаемых, оцениваемых параметром ω в условии (10), пренебречь членами связанности в уравнении теплопроводности. Как показали дополнительные расчеты, проведенные авторами [19], учет инерционных членов в уравнениях (1), (2) в исследуемом случае (при ЭМО металлов и покрытий) дает поправку в значениях температур и напряжений менее 1 %, при этом время расчета возрастает более чем на 20 % по каждому из этапов решения. 5. Результаты расчета Далее приводятся результаты расчета напряжений при ЭМО стального образца с учетом образующейся структурной неоднородности и упруго-пластических деформаций. Решение задачи искалось при следующих исходных данных: размеры расчетной области - параллелепипед размерами 5×5×1 мм, выделенный из полупространства (на сетке 100×100×50 узлов); толщина поверхностного слоя (покрытия) - 0,15 мм; материал - сталь 45 нормализованная и закаленная; скорость обработки (движения источника) - 5 м/мин; размеры зоны воздействия (большая и малая полуоси эллипса) - 1,4´0,9 мм; ток (амплитуда) - 400 А, плотность тока - 404 А/мм2, частота - 50 Гц; напряжение - 3 В; плотность теплового потока - 5,8∙108 Вт/м2; деформирующее усилие в зоне обработки - 1000 Н. Получены расчетные картины формирования упрочненных фрагментов в поверхностном слое закаленной стали 45 в различные моменты времени в ходе ЭМО (рис. 8). Иллюстрируется образование отдельных упрочненных фрагментов, соответствующих единичному тепловому импульсу электрического тока. Показано постепенное продвижение зоны разупрочнения (составленной из продуктов высокоскоростного отпуска исходной мартенситной структуры) в поверхностном слое на глубину около 1 мм (зона 3 на рис. 2) и формирование закаленной зоны на глубину до 0,2 мм с образованием регулярной дискретной структуры (зона 2 на рис. 2). а б в г Рис. 8. Расчетные картины формирования неоднородной структуры деформируемого тела на примере технологии упрочнения поверхностного слоя стали 45 в ходе ЭМО в различные моменты времени от начала процесса: а - 7,5 мс; б - 17,5 мс; в - 27,5 мс; г - 37,5 мс Для сравнения на рис. 9, а показана фотография, полученная на растровом электронном микроскопе, двух соседних упрочненных в ходе ЭМО фрагментов («белый слой») и зоны отпуска между ними и в исходном материале (закаленная сталь 45), а на рис. 9, б - их сопоставление с соответствующими расчетными картинами структуры упрочненного трека. а б Рис. 9. Фрагменты белого слоя (а) после ЭМО поверхности закаленной стали 45 (×400) и сопоставление (б) расчетных и экспериментальных данных о распределении структур в упрочненном треке (БС - белый слой; ЗО - зона отпуска) На рис. 10 показаны расчетные картины распределения интенсивности напряжений (в абсолютных и относительных (к пределу текучести σт) единицах) в объеме поверхностного слоя материала в ходе ЭМО для момента времени, соответствующего максимуму теплового потока (7,5 мс). Сопоставляются результаты решения термоупругопластической задачи для однородного тела (рис. 10, а, в) и для двухслойной композиции, моделирующей в первом приближении материал со сплошным (без учета отдельных фрагментов, треков и зон отпуска) (см. рис. 2) упрочненным поверхностным слоем (рис. 10, б, г). В объеме материала в ходе ЭМО формируется сложное НДС. В зоне обработки под действием контактной нагрузки возникают высокие сжимающие напряжения, максимальная интенсивность σi которых на поверхности тела достигает 300-400 МПа для однородного тела и около 700 МПа - для неоднородного материала с упрочненным слоем на поверхности. При этом для инструмента в форме эллипсоидного штампа наибольшая интенсивность достигается на периферии зоны контакта (тогда как для равномерно распределенной нагрузки - в центре данной зоны). По мере удаления от области контакта интенсивность напряжений снижается (до 1,5-2 раз на глубине 0,15 мм, соответствующей границе упрочненного слоя с исходным материалом), а с достижением предела текучести σт в данной точке упругопластические деформации распространяются на соседние области (прежде всего - по глубине материала), где уровень σi также приближается к σт. Формирование тонкого высокопрочного поверхностного слоя в ходе контактного температурно-силового воздействия приводит к перераспределению напряжений в неоднородном теле - упрочненная поверхность воспринимает наибольшие напряжения и деформации, разгружая подложку. На границе слоя и подложки уровень σi также несколько выше интенсивности напряжений в однородном теле на соответствующей глубине (0,15 мм) от поверхности (см. рис. 10, кривые 3 и 4). Рис. 10. Распределение интенсивностей напряжений в объеме расчетной области (а, б) и их относительных (к пределу текучести σт) величин (в, г) в поперечном сечении через центр зоны контакта инструмента и детали для однородного (а, в; кривые 1, 3, 5) и неоднородного тела (б, г; кривые 2, 4, 6): 1, 2 - на глубине z = 0 мм; 3, 4 - z = 0,15 мм; 5, 6 - z = 0,5 мм Расчеты, проведенные в рамках разработанной модели, позволили определить картины распределения температурных, фазовых и суммарных остаточных напряжений по глубине упрочненного ЭМО слоя (рис. 11). Анализ полученных зависимостей подтвердил, что формирование величины и знака остаточных напряжений в приповерхностных областях происходит в соответствии с двумя «конкурирующими» механизмами. При этом термические напряжения играют определяющую роль при высоких температурах на стадии нагрева, а также на начальном этапе охлаждения и приводят к формированию на поверхности материала растягивающих остаточных напряжений. Фазовые напряжения начинают сказываться при охлаждении в более низком интервале температур, соответствующих мартенситному превращению, и являются причиной возникновения в упрочненном слое сжимающих остаточных напряжений. Рис. 11. Сопоставление расчетных (линия 1) и экспериментальных (линия 2 - И.И. Колпакас; 3 - М.С. Нерубай, А.П. Овчинников; 4, 5, 6 - Б.М. Аскинази) данных о распределении окружных остаточных напряжений по сечению материала после ЭМО (глубина упрочненного слоя - 0,2 мм) Сопоставление (см. рис. 11) полученных теоретических результатов с экспериментальными данными (Б.М. Аскинази [27], А.Г. Григорьянца, И.И. Колпакаса [28], М.С. Нерубая [29] и др.) свидетельствует о достигнутом согласовании по характеру распределения напряжений по сечению материала. Выводы 1. Показано, что при решении рассматриваемого класса связанных технологических задач необходимо систему основных уравнений механики термоупругопластических деформаций дополнить специфическими соотношениями вычислительного материаловедения, устанавливающими связь между термо-напряженным состоянием в точке тела и образующейся в ней структурой. 2. Обоснованы важность и необходимость согласования шагов временных сеток для каждой из рассматриваемых связанных задач (тепловой, механической, структурной) для повышения точности описания формирующегося состояния тела. 3. Проведенный анализ решения задач (на основе известных условий оценки связанности и инерционных эффектов) в приложении к технологии поверхностного упрочнения импульсной ЭМО позволил ввести корректные упрощения в тепловой задаче (на базе параболического уравнения теплопроводности) и в описании тонких эффектов влияния деформаций на тепловыделение. 4. С учетом описанных особенностей разработан итерационный алгоритм численного описания прямых и обратных связей между тепловой, упругопластической задачами и процессами структурообразования. 5. На основе сопоставительного анализа собственных и известных расчетных и экспериментальных результатов выявлены особенности распределения полей напряжений и деформаций в телах с трансформирующейся структурой и показана существенная роль неоднородности строения в формировании конечного состояния изделия.

About the authors

V P Bagmutov

Volgograd State Technical University

Email: sopromat@vstu.ru

I N Zakharov

Volgograd State Technical University

Email: zaxap@mail.ru

D S Denisevich

Volgograd State Technical University

Email: adven148@ya.ru

References

Statistics

Views

Abstract - 230

PDF (Russian) - 69