Full Text

Введение Метод инфракрасной термографии (ИКТ) является бесконтактным методом измерения температуры поверхности объектов и применяется как для контроля температуры образцов во время механических испытаний [1], так и для определения механических и термодинамических характеристик материала, таких как предел усталости [2, 3], критическая величина деформации для быстрых трещин [4], температуропроводность [5] и т.д. В России основной областью применения метода ИКТ является неразрушающий контроль [6]. В последние два десятилетия метод ИКТ активно используется для детектирования возникновения и развития повреждений (трещин, пор и других неоднородностей) материала и исследования процесса усталостного разрушения. Исследование поведения материала в условиях многоцикловой усталости и прогнозирование его выносливости на основе метода ИКТ представлено в [7, 8]. В работе [9] показана возможность использования ИКТ для прогнозирования усталостной долговечности сварных соединений. На основе данных метода ИКТ решается и ряд задач, посвященных определению энергетических характеристик материала, таких как скорость диссипации и накопления энергии в процессе деформирования [10, 11, 12, 13, 14]. Оценка значений различных энергий (пластической, накопленной, диссипированной), вовлеченных в процесс деформирования материала, играет важную роль по нескольким причинам. Во-первых, их измерение дает возможность глубже проникнуть в суть механизмов деформации и разрушения. Во-вторых, энергетическая концепция описывает изменение внутренней энергии материала и, таким образом, вносит вклад в построение и валидацию моделей, базирующихся на принципах термодинамики. В-третьих, информация об эволюции различных энергетических параметров может быть использована для идентификации коэффициентов (например, коэффициента Тейлора-Куинни), часто использумых при численном моделировании адиабатических процессов для оценки количества работы пластической деформации, перешедшей в тепловую энергию [15]. Явными преимуществами метода ИКТ являются возможность бесконтактного измерения температуры поверхности материала в различных условиях и получение поля мощности источников тепла на поверхности образца. Несмотря на это, данный метод обладает рядом ограничений, связанных с излучательной способностью исследуемого материала, зашумлением сигнала, вызванным влиянием внешних факторов, условиями теплообмена образца и окружающей среды, а также ограничением по точности определения мощности источников тепла. Указанные факторы не позволяют использовать метод инфракрасной термографии в эксплуатационных условиях материалов и конструкций для определения их энергетического состояния. В данной работе сделана попытка верифицировать величину рассчитываемой на основе метода ИКТ мощности источников тепла на поверхности материала, возникающих при его деформировании. Для этого использовался разработанный авторами контактный датчик потока тепла [16], устройство которого основано на эффекте Зеебека. 1. Материалы и условия эксперимента Экспериментальное исследование проводилось на серии образцов из титанового сплава ОТ4-0. Геометрия образцов соответствовала ГОСТ 1497-84, размеры образцов (мм) представлены на рис. 1. Химический состав приведен в таблице 1. Рис. 1. Геометрия образцов Таблица 1 Химический состав сплава титана ОТ4-0, % Элемент AL Mn Zr* Si* Fe C* AL O2* Требования ОСТ1-90013-81 £1,4 £1,3 £0,30 £0,15 £1,5 £0,10 £1,4 £0,18 Факт. 1,1 0,9 0,3 0,004 1,12 0,008 1,1 0,082 Элемент N2* H2* S прочих Cr+Mn* Cu+Nr* Требования ОСТ1-90013-81 £0,05 £0,015 £0,30 £0,15 £0,10 Факт. 0,012 0,003 0,022 - 0,020 Примечание. Ti - основа, * - примеси. Пробоподготовка осуществлялась в два этапа. Сначала поверхность образцов, предназначенная для инфракрасной (ИК) съемки, полировалась абразивной бумагой (на завершающей стадии полировки размер абразивных частиц не превышал 20 мкм), затем поверхность покрывалась тонким слоем аморфного углерода. Поверхность образцов, на которой устанавливался контактный датчик потока тепла, не требовала дополнительной подготовки. Перед испытанием на нее крепились метки для регистрации перемещения в ходе механического тестирования. На область примыкания датчика потока тепла и поверхности образца наносилась термопаста для улучшения характеристик теплового контакта. На рис. 2 схематично изображена экспериментальная установка, включающая в себя испытательную электромеханическую машину Shimadzu AG-X Plus (300 кН), контактный датчик потока тепла, ИК-камеру FLIR SC5000 и видеоэкстензометр Shimadzu TRViewX240S. ИК-камера имеет следующие характеристики: спектральный диапазон 3-5 мкм, максимальный размер кадра 320×256 точек, пространственная чувствительность составляет 10-4 м, минимальная погрешность измерения температуры 25 мК при температуре образца 300 К. После установки образца в захваты испытательной машины вся система экранировалась ширмой с целью исключения влияния факторов внешней среды и паразитных переотражений. ИК-съемка и запись данных датчика потока тепла, установленного на образце, велись непрерывно во время всего механического теста. Частоты записи ИК-съемки составляла 25 Гц. Рис. 2. Схема экспериментальной установки по исследованию баланса энергии в материале при квазистатическом растяжении: 1 - образец; 2 - контактный датчик потока тепла; 3 - видеокамера для регистрации смещения меток; 4 - ИК-камера; 5 - захваты испытательной машины; 6 - метки для регистрации перемещения Контактный датчик потока тепла предварительно калибровался. Калибровочная функция показывает связь между разностью потенциалов U (В), возникающей из-за разности температур двух поверхностей элементов Пельтье, из которых состоит контактный датчик потока тепла, и мощностью потока тепла в Вт. (1) После настройки всех систем и экранирования установки для каждого образца проводились дополнительные эксперименты по определению условий теплообмена образца с окружающей средой. Эксперименты заключались в регистрации поля температуры поверхности образца после его импульсного точечного нагрева и построения кривой остывания, по которой рассчитывалась константа теплообмена. После этого проводились механические испытания на квазистатическое растяжение плоских образцов с одновременной регистрацией поля температуры поверхности образцка (ИК-съемка), потока тепла (контактным датчиком) и перемещения рабочей зоны образца (по видеоэкстензометру). Скорость движения захватов варьировалась и в различных испытаниях составляла 5 и 10 мм/мин, что соответствовало скорости деформации 7,7∙10-4 с-1 и 15∙10-4 с-1. На рис. 3 показана характерная диаграмма деформирования образцов, построенная на основе данных видеоэкстензометра. Рис. 3. Характерная диаграмма деформирования образцов из титанового сплава ОТ4-0 2. Оценка мощности источников тепла на основе уравнения теплопроводности Расчет поля мощности источников тепла по данным ИКТ осуществлялся на основе уравнения теплопроводности (2) где T(x,y,z,t) - поле температуры; ρ - плотность материала; c - удельная теплоемкость; k - коэффициент теплопроводности; Q(x,y,z,t) - поле мощности источников тепла; x,y,z - декартовы координаты; t - время. Метод ИКТ позволяет регистрировать поле температуры поверхности образца. Контроль распределения температуры по толщине образца данным методом невозможен, поэтому в экспериментальных исследованиях используются достаточно тонкие образцы, и предполагается, что в этом случае распределение температуры по толщине образца однородно [17, 18]. Для оценки интегральной мощности источника тепла используется уравнение (2), осредненное по объему. Была проведена стандартная процедура осреднения. Разницу θ´(t) между средней по объему температурой образца и начальной температурой образца T0, находившегося в тепловом равновесии с окружающей средой, определим следующим образом: (3) где T0 - начальная температура образца, находящегося в тепловом равновесии с окружающей средой; a, b, h - длина, ширина и толщина образца соответственно; V - объем. Запишем граничные условия в виде (4) где gx - коэффициент теплообмена образца с окружающей средой на соответствующей грани образца. Граничные условия по остальным двум направлениям имеют аналогичный вид за исключением коэффициентов gi, i = {x,y,z}. Интегрируя уравнение (2) по объему с учетом выражения (3) и граничных условий (4), получим усредненное по объему уравнение теплопроводности для расчета поля мощности источников тепла: (5) где θ - средняя температура исследуемой поверхности; T0 - начальная температура исследуемой поверхности, находившейся в тепловом равновесии с окружающей средой; m - масса области, по которой берется средняя температура (плотность титанового сплава ОТ4-0 ρ = 4457 кг/м3 [19, 20]); c - удельная теплоемкость (для ОТ4-0 456 Дж/(кг∙К) [19, 20]); S(t) - мощность источников тепла (Вт); β - параметр материала, определяющий потери тепла, связанные с теплообменом с окружающей средой. Параметр β определяется экспериментально из тестов по остыванию образцов после их импульсного точечного нагрева. 3. Экспериментальное определение констант теплообмена образца с окружающей средой Для расчета мощности источников тепла на основе уравнения (5) по экспериментально полученному полю температур необходимо определить параметр β, отвечающий за потери тепла при взаимодействии образца с окружающей средой. Для этого были проведены дополнительные эксперименты по остыванию образца после его импульсного точечного нагрева. На рис. 4 представлена характерная ИК-картина, получаемая при импульсном точечном нагреве образца и экспериментальные данные средней температуры области нагрева. Остывание образца после импульсного точечного нагрева можно описать усредненным уравнением теплопроводности (5) с нулевой мощностью источника: (6) Решением этого уравнения будет являться функция вида (7) где C1 и C2 - константы; ρ - плотность. Для нахождения константы β необходимо аппроксимировать экспериментальные данные изменения средней температуры образца после импульсного точечного нагрева. В качестве аппроксимирующей функции выберем функцию вида (8) Время, с Рис. 4. Характерное ИК-изображение образца при его импульсном точечном нагреве и график зависимости средней температуры области нагрева от времени Для оценки параметра теплообмена образца с окружающей средой аппроксимировалась только линейная в полулогарифмических координатах часть кривой остывания образца после импульсного точечного нагрева (выделена зеленой линией на рис. 5), то есть часть, которая соответствует экспоненциальной зависимости температуры от времени (7). Параметр β зависит от условий окружающей среды и не является постоянной величиной для данного материала, поэтому его необходимо определять непосредственно перед началом механических испытаний. В табл. 2 приведены характерные данные, полученные для одного из испытанных образцов: величина произведения плотности и теплоемкости материала, взятых из справочных данных [19, 20], экспериментально полученная величина степени экспоненты B из уравнения (8) и рассчитанный по этим значениям параметр β. In(t) Рис. 5. Аппроксимация участка кривой остывания образца после импульсного точечного нагрева: Test1: Test2: Test3: Таблица 2 Результат аппроксимации кривой остывания образца после импульсного точечного нагрева Номер теста Произведение плотности и коэффициента теплопроводности ρc, Дж/(м3К) Степень экспоненты в аппроксимации B β 1 2,04∙106 0,0080 1,632∙104 2 2,04∙106 0,0061 1,239∙104 3 2,04∙106 0,0065 1,331∙104 Среднее: 1,4∙104 В результате согласно описанной выше методике был проведен эксперимент по остыванию образца после его импульсного точечного нагрева, и определены величины параметра β для каждого образца перед его деформированием. 4. Экспериментальное определение величины интегрального источника тепла по данным метода ИКТ На рис. 6 представлено характерное ИК-изображение поверхности образца на завершающей стадии деформирования непосредственно перед моментом разрушения и экспериментальная зависимость средней по выделенной области температуры от времени всего механического теста. На начальном этапе в результате термоупругого эффекта средняя температура поверхности образца падает, затем начинает превалировать термопластический эффект, и средняя температура образца начинает расти вплоть до момента образования шейки и полного разрушения. На основе данных об изменении средней температуры поверхности образца во время механического испытания и полученной экспериментально величине параметра β по формуле (5) была определена зависимость мощности источника тепла от времени. Время, с а б Рис. 6. Характерное ИК-изображение поверхности образца при квазистатическом растяжении в момент перед разрушением (а) и зависимость средней по выделенной области температуры от времени всего испытания (б) Согласно данным, приведенным в табл. 2, колебания параметра β составляют порядка 30 %. На рис. 7 представлено характерное изменение величины мощности источника тепла от времени с указанием доверительного интервала. На начальном этапе нагружения (0-70 с) колебание параметра β практически не сказывается на величине мощности источника тепла. Начиная с 70-й секунды, влияние параметра β становится значительным. В этот период времени максимальная величина стандартного отклонения мощности источника тепла, рассчитанного при средней величине β, составляет 0,27. Рис. 7. Величина мощности источника тепла, рассчитываемая по данным ИКТ при среднем значении β с указанием доверительного интервала На рис. 8, 9 представлены зависимости средней по области на поверхности образца температуры, рассчитанной мощности источников тепла и силы прикладываемой нагрузки от времени для двух скоростей деформирования 7,7∙10-4 и 15∙10-4 с-1. а б в Рис. 8. Характерная зависимость средней по области температуры от времени при квазистатическом растяжении (а); соответствующая зависимость интегральной величины мощности источника тепла от времени по данным ИКТ (б); диаграмма нагрузки (в). Скорость деформации 7,7∙10-4 с-1 а б в Рис. 9. Характерная зависимость средней по области температуры от времени при квазистатическом растяжении (а); соответствующая зависимость интегральной величины мощности источника тепла от времени по данным ИКТ (б); диаграмма нагрузки (в). Скорость деформации 15∙10-4 с-1 Согласно полученным зависимостям практически постоянный тепловой поток устанавливается в промежуток времени 40-250 с при скорости деформирования 7,7∙10-4 и 25-120 с при скорости деформирования 15∙10-4. На диаграмме нагружения в этот момент наблюдается площадка текучести и слабовыраженное разупрочнение. К моменту образования шейки тепловой поток резко возрастает, что свидетельствует о переходе материала к разрушению. 5. Сравнение результатов и обсуждение Контактный датчик потока тепла непрерывно фиксировал величину мощности источника тепла в той области образца, где он был закреплен. Такая возможность регистрации интегрального потока является одним из преимуществ разработанного датчика, которое позволяет отслеживать энергетическое состояние материала in situ. На рис. 10 представлены характерные графики изменения мощности источника тепла во время механических испытаний при различных скоростях деформирования. Графики получены с помощью контактного датчика потока тепла и на основе обработки данных бесконтактного метода регистрации температуры (ИКТ). а б Рис. 10. Сравнение данных, получаемых на основе метода ИКТ, и датчика потока тепла: а - скорость деформации 7,7∙10-4 с-1; б - скорость деформации 15∙10-4 с-1 Метод ИКТ демонстрирует бόльшую чувствительность к изменению теплового потока, чем контактный способ его измерения, но требует значительного времени для обработки данных, учета условий теплообмена и нахождения параметра β, расчета мощности источника тепла по разностной схеме уравнения (5). С другой стороны, контактный датчик, регистрируя непрерывно интегральную величину потока тепла, не позволяет оценить пространственное распределение источников и отслеживать место локализации пластической деформации материала. Не абсолютное соответствие данных теплового потока, получаемых двумя методами, можно объяснить различной чувствительностью приборов, погрешностями численной обработки данных метода ИКТ и определения параметра теплообмена образца с окружающей средой. Открытым остается вопрос о теплоинерционности контактного датчика. В статье [21] приведены сведения о точности показаний датчика при изменяющемся потоке тепла, вызванном источником слабой мощности. При больших значениях потоков тепла инерция датчика должна быть более заметна и может достигать нескольких секунд. Эта задача требует дальнейшего исследования и создания дополнительных алгоритмов обработки данных. Несмотря на это, полученные результаты позволяют говорить о возможности использования контактных и бесконтактных измерений как в комплексе - для верификации величины источников тепла, их распределения на поверхности материала и определения условий теплообмена образца с окружающей средой, так и по отдельности - в качестве экспресс методики оценки состояния материала на различных этапах нагружения.

About the authors

A Yu Iziumova

Institute of continuous media mechanics UB RAS

Email: fedorova@icmm.ru

A N Vshivkov

Institute of continuous media mechanics UB RAS

Email: vshivkov.a@icmm.ru

A E Prokhorov

Institute of continuous media mechanics UB RAS

Email: prokhorov.a@icmm.ru

O A Plekhov

Institute of continuous media mechanics UB RAS

Email: poa@icmm.ru

B Venkatraman

Indira Gandhi Centre for Atomic Research

Email: bvenkat@igcar.gov.in

References

Statistics

Views

Abstract - 272

PDF (Russian) - 97