УДК 620.199.001.18 В. В. Носов
ОЦЕНКА РЕСУРСА ТЕХНИЧЕСКИХ ОБЪЕКТОВ НА ОСНОВЕ МОДЕЛИРОВАНИЯ ВРЕМЕННОЙ ЗАВИСИМОСТИ ПАРАМЕТРОВ АКУСТИЧЕСКОЙ ЭМИССИИ
UDC 620.199.001.18 V. V. Nosov
THE ASSESSMENT OF THE RESOURCE OF TECHNICAL OBJECTS BASED ON MODELING THE TIME DEPENDENCE OF PARAMETERS OF ACOUSTIC EMISSION
Аннотация
Рассмотрены проблемы определения ресурса технического объекта. Изложен метод нахождения времени до разрушения на основании анализа сигналов акустической эмиссии. В результате проведенных исследований увеличена точность методики определения ресурса на основе идентификации участка однородного разрушения.
Ключевые слова:
акустическая эмиссия, однородное разрушение, ресурс, диагностика, энергия, прочность, технический объект, моделирование.
Abstract
The paper deals with the problems of determining the resource of a technical object. The method of finding life time before failure based on the analysis of acoustic emission signals is presented. The research resulted in increasing the accuracy of methods for determining the resource based on the identification of the area of homogeneous destruction.
Key words:
acoustic emission, homogeneous destruction, resource, diagnostics, energy, strength, technical object, modeling.
Несмотря на то, что метод акустической эмиссии (АЭ) признан как наиболее информативный в отношении диагностики процесса разрушения, проблема оценки ресурса (времени до разрушения) ещё не решена [1]. В работе рассматривается метод нахождения времени до разрушения на основании развития изложенных в [2-6] идей.
Результаты разработок используются в экспертной деятельности при оценке промышленной безопасности объектов химической, нефтехимиче-
ской, нефтяной, газовой, пищевой отраслей [7].
Основной задачей являлась проверка возможности идентификации (выделения) участка однородного разрушения по статистическим признакам (вариациям значений параметров АЭ) и возможного повышения точности оценки ресурса. Программа, реализуемая в методике, использует для аппроксимации этап однородного разрушения, определение которого способствует увеличению точности метода.
©Носов в. В., 2013
Описание методики
Одной из основных причин нарушения работоспособности и уменьшения ресурса конструкций и деталей машин является развитие в их материале процесса накопления повреждений. Для описания этого процесса используется микромеханическая модель разрушения, согласно которой временная зависимость поврежденности материала в условиях мелкодисперсного роста концентрации микротрещин С имеет следующий вид:
®СР = Т0 ехР
С (/) = Со
1 - ехр[-{^-] , (1)
0 ®СР О )
Л
где I - текущее время; С0 - начальная концентрация структурных элементов в материале до разрушения; 0СР - среднестатистическое время ожидания разрушения одного структурного элемента, задаваемое формулой Журкова; при начальном условии С(0) = 0.
Формула Журкова в условиях растяжения при постоянстве напряжения имеет вид:
'и0-Г° ЯТ
(2)
где т0 - период атомных колебаний,
13
т0 ~ 10 с; и0 - энергия активации процесса разрушения, - структурно-чувствительный коэффициент; Я - универсальная газовая постоянная; Т - абсолютная температура.
Из-за неоднородности структурного и напряженного состояния материала среднестатистические значения времени разрушения структурных элементов тела в общем случае неодинаковы. Таким образом, для различных элементов параметр 0 принимает разные значения. Для характеристики неоднородности состояния материала вводится функция ¥(у) распределения параметра у по структурным элементам, в которой в зависимости от степени неоднородности материала используются различные распределения. В работе при описании разрушения образцов сварных соединений рассматриваются два из них: ограниченное вейбулловское распределение с параметрами к, X, д:
у ) =
0, 7 е[0, д),
(кЛ(7^к-1 (( ^^
Л
Л
ехр
7_
Л
, 7^[ q, ;
^ "У У
(3)
А =-
1,1 Л Л
чк-1 (г \к\
У\
ехр
Л
ё7
и логарифмически-нормальное с параметрами ог, ¡и:
¥(7) =
1
-ехр
■я<т2у
(Му) -М)2
2а\
. (4)
Зависимость концентрации микротрещин от времени для неоднородного материала можно представить в виде
1
/шах
С к) = Со I у(у)
/шт
Вестник Белорусско-Российского университета. 2013. № 2 (39) ( ( \\
1 - ехр
-1-
о
(
т0 ехр
и0-уа(г)
ЯТ
Л
)))
йу,
(5)
где уш;п и ушах - минимальное и максимальное значения коэффициента у для данного распределения у/(у).
Переход к макроскопическому разрушению материала происходит, ко-
гда число разрушенных структурных элементов составляет 1 % от числа начальных: С/С0 = 0,01. Время до разрушения находят из уравнения
\Л
шах / \
1 у(/)
1 - ехр
1
Уш\п
~т\\
т0 ехр
и-уо(х)
йу = 0,01.
(6)
ЯТ
)))
Для вычисления времени до разрушения технического объекта необходимо определить значения величин, входящих в (6), знать зависимость напряжения от времени, температуру, величину энергии активации и0 (постоянна для конкретного материала) и пара-
метры функции Ф(у). Для нахождения последних используется регистрация сигналов АЭ при нагружении объекта контроля в области упругих деформаций. Число сигналов АЭ N-¡0 пропорционально числу разрушенных структурных элементов:
N¡(0 = кЛЕС(!)
или
NAt ) = кЛЕС0 I у (г) 7шт
1 - ехр
-1-
\Л
т0 ехр
и0)
ЯТ
(7)
) ))
где кЛЕ - акустико-эмиссионный коэффициент (может быть принят постоянным при нагружении с постоянной скоростью в области упругих деформаций).
Параметры функции ¥(у) находят путем аппроксимации экспериментальной зависимости числа импульсов АЭ от времени теоретической кривой (7). Далее из (6) рассчитывается время до разрушения технического объекта.
В работе предполагается определе-
ние параметров функции Ф(у) для образцов стыковых сварных соединений. Материал пластин - сталь СтЗпс, сварка -ручная и автоматическая в среде углекислого газа. Образцы испытывались на растяжение с постоянной скоростью роста деформации, при этом регистрировались сигналы АЭ. Нагружение образцов производилось следующим образом (рис. 1).
Рис. 1. Схема нагружения образца: 1 - образец; 2 - сварной шов; 3 - верхний захват нагружающего устройства; 4 - нижний захват нагружающего устройства; 5 - палец; 6 - ПАЭ
Концы образца 1 закреплялись в захватах 3 и 4 испытательной установки посредством пальцев 5. Верхний захват 3 оставался неподвижным, а нижний захват 4 перемещался вертикально вниз с постоянной скоростью. Сила, приложенная к образцу, определялась по шкале установки (цена деления шкалы 3000 Н). Для регистрации сигналов АЭ была использована двухканальная акус-тико-эмиссионная система СДАЭ-16(2). Система состояла из двух преобразователей акустической эмиссии (ПАЭ) типа П113с, предварительных усилителей, блока регулирования усиления, блоков обработки сигналов АЭ, блока управления и компьютера. Принцип работы одного ее канала следующий: акустический сигнал, возникающий в материале образца при нагружении, превращается пьезоэлектрическим преобразователем в электрический импульс, который подвергается усилению и обработке, после чего информация
поступает в компьютер, где отображается в виде таблицы данных и графиков зависимости от времени различных параметров АЭ сигнала. Регистрировались число импульсов АЭ, их амплитуда, энергия, длительность, суммарный счет АЭ и время нарастания сигнала.
Рабочий диапазон частот аппаратуры составлял 20... 1000 кГц, эффективное значение напряжений собственных шумов усилительного тракта, приведенное ко входу предварительного усилителя, не превышало 5 мкВ. Нижний уровень порога дискриминации системы определялся необходимостью устранения электромагнитных помех и составлял 47 дБ на каждом канале. Для регистрации сигналов были выбраны преобразователи с диапазоном частот 20.500 кГц.
Для обработки первичных данных АЭ-испытаний, определения коэффициентов микромеханической модели раз-
рушения предназначена программа RG, состоящая из нескольких макросов Microsoft Excel.
Одним из этапов выполнения работы является нахождение интервала для аппроксимации. Требуется выбрать линейный участок зависимости логарифма числа импульсов от времени (этап однородного разрушения), а также разработать методологию нахождения данного участка. Были произведены испытания 30-ти образцов. Проанализированы временная зависимость числа импульсов, логарифм числа импульсов, длительность импульсов, коэффициент перекрытия импульсов (вычисляется как отношение суммарной длительности импульсов в течение интервала времени к величине этого интервала), амплитуда, коэффициент вариации амплитуды, параметр Knt (отношение активности АЭ к средней длительности импульсов). В результате анализа и выявления зависимостей выделены следующие этапы для определения начальной и конечной точек прямолинейного участка:
— найти на графике зависимости длительности импульсов АЭ от времени момент времени tb когда наблюдается резкий рост длительности импульсов
(свидетельствует о начале пластического деформирования);
— найти на графике зависимости коэффициента перекрытия от времени момент времени t2, соответствующий началу роста коэффициента перекрытия. Он свидетельствует о начале пластического деформирования, при котором наблюдаются движение дислокаций, мелкие перестройки в структуре материала, а значит, имеется большое количество «коротких» сигналов, которые, накладываясь друг на друга, сливаются, что отражается в виде скачка на графике зависимости коэффициента перекрытия от времени;
— выбрать меньшее из двух значений времени (ti и t2).
Пример определения интервала времени для аппроксимации.
1. Резкий рост длительности импульсов наблюдается в момент времени ti = 186 c (рис. 2).
2. Коэффициент перекрытия импульсов начинает расти в момент времени t2 = 175 c (рис. 3).
3. Для проведения дальнейших действий используется величина времени t2 < t1 (рис. 4).
о
го о .0
о о
X .0
CD
14000 12000 10000 SOOO 6000 4000 2000 О
--t,= li с,- 1 ►
U L — ♦
♦
_^ h , ць . л. ► ж
► Ч И t ♦
" -............
О 100 200 300 400 500 600 700 300
Время, с
Рис. 2. Определение момента времени ^ по графику зависимости длительности импульсов от времени
100 200 ЗОО +00 БОО
Время, с
600
700
аоо
Рис. 3. Определение момента времени ^ по графику зависимости коэффициента перекрытия импульсов от времени
.о
с
м
го
го
£ 0.060 со
аЗ 0,040
т о
)Р 200 300 400 500 600 700 800
Время, с
Рис. 4. Определение момента времени /1 по графику зависимости коэффициента вариации амплитуды от времени
Изменение характера сигналов АЭ указывает на появление пластической деформации материала образца в момент времени Ь = 175 с.
Далее требуется найти первоначальную точку линейного участка временной зависимости логарифма числа сигналов АЭ (начало этапа однородного разрушения). Для этого рассмотрим график зависимости коэффициента вариации амплитуды от времени. На участке от 0 до Ь (движение от Ь к 0) выделяем t1, в котором происходит аномаль-
ный скачок коэффициента вариации амплитуды.
Скачок вариации амплитуд объясняется разрушением связей с большим разбросом размеров, структур и напряжений. Разрушение неоднородно, следовательно, возникают аномально большие амплитуды.
Полученный интервал вносится в макрос программы, с помощью которого строится теоретическая кривая зависимости числа импульсов АЭ от времени. При этом используется в качестве функ-
ции ¥(у) ограниченное вейбулловское распределение. Изменяя величину параметра X, можно проследить за изменением среднего значения параметра у распределения ¥(у). Путем изменения параметров X, д и кАЕС0 (к = 3; одинаков для всех образцов) добиваемся совпадения теоретической и экспериментальной кривых N2,^). Величина средней ошибки аппроксимации должна быть минимальной. При совпадении кривых получаем искомые значения параметров микромеханической модели разрушения.
Вывод
В результате авторами был получен метод определения участка однородного разрушения. С помощью метода была увеличена точность определе-
ния времени до разрушения объекта на 2 % (по сравнению с ранее предложенной методикой). Сократилась погрешность с 15 до 13 %, что соответствует нормам. Таким образом, разработанная методика обладает повышенной точностью. На основе АЭ-показателей, напрямую связанных с определяющим прочность процессом и статистически устойчивых к помехам регистрации, получены файлы MS Excel, представляющие собой базу, позволяющую выполнять различные вычисления по желанию пользователя. Методика проста в применении, универсальна, т. к. диагностические параметры позволяют оценить ресурс и стандартные характеристики прочности различных по структуре и размерам объектов.
СПИСОК ЛИТЕРАТУРЫ
1. Акустико-эмиссионный контроль железнодорожных конструкций / Л. Н. Серьёзнов [и др.]. -Новосибирск : Наука, 2011. - 272 с.
2. Лахова, Е. Н. Моделирование процесса перестройки структуры при пластической деформации сварных соединений металлоконструкций / Е. Н. Лахова, В. В. Носов // Современное машиностроение. Наука и образование : материалы Междунар. науч.-практ. конф., Санкт-Петербург 14-15 июня 2011 г. -СПб. : СПбПУ, 2011. - С. 301-309.
3. Носов, В. В. Диагностика машин и оборудования : учеб. пособие / В. В. Носов. - 2-е изд., испр. и доп. - СПб. : Лань, 2012. - 384 с.
4. Nosov, V. V. The Effect of Strength State Nonuniformity on the Acoustic Emission of Structural Materials / V. V. Nosov, G. S. Elchaninov // Russian Journal of Nondestructive Testing. - 2011. - Vol. 47, № 12. - Р. 824-833.
5. Носов, В. В. Оценка качества технологического процесса на основе структурно-силового анализа прочностной неоднородности материала изделия / В. В. Носов, Е. Н. Лахова // Заводская лаборатория. Диагностика материалов. - 2012. - № 8. - С. 69-71.
6. Носов, В. В. Неразрушающий контроль качества заготовок для производства горячекатаной полосы методом акустической эмиссии / В. В. Носов, В. Г. Лаврин // Дефектоскопия. - 2012. - № 3. -С. 18-26.
7. Пат. 2445615 РФ. Способ определения прочности материала изделий / В. В. Носов, Е. Н. Лахова. -№ 2010145698/28 ; заявл. 09.11.10 ; опубл. 20.03.12, Бюл. № 8. - 4 с.
Статья сдана в редакцию 29 марта 2013 года
Виктор Владимирович Носов, д-р техн. наук, проф., Национальный минерально-сырьевой университет «Горный». Тел.: (904) 617-94-38. E-mail: nosovvv@list.ru.
Viktor Vladimirovich Nosov, DSc (Engineering), Prof., National Mineral Resources University (Mining University). Tel.: (904) 617-94-38. E-mail: nosovvv@list.ru.