УДК 620.199.001.18
Т.Д. Ахметшин, В.В. Носов
ОПРЕДЕЛЕНИЕ РЕСУРСА ТЕХНИЧЕСКИХ ОБЪЕКТОВ НА ОСНОВЕ ИДЕНТИФИКАЦИИ ЭТАПА ОДНОРОДНОГО РАЗРУШЕНИЯ
T.D. Akhmetshin, V.V. Nosov
ASSESSMENT OF THE STRENGTH OF CONSTRUCTION MATERIALS FOUNDED IN THE PHASE OF HOMOGENEOUS DESTRUCTION
Статья посвящена проблеме определения ресурса технических объектов, которая сводится к нахождению времени, оставшегося до наступления их предельного состояния, неразрушающим методом — на основании анализа сигналов акустической эмиссии (АЭ) с позиции микромеханической модели разрушения и временных параметров АЭ.
НЕРАЗРУШАЮЩИЙ КОНТРОЛЬ. АКУСТИЧЕСКАЯ ЭМИССИЯ. ОДНОРОДНОЕ РАЗРУШЕНИЕ. РЕСУРС. МИКРОМЕХАНИЧЕСКАЯ МОДЕЛЬ.
The article is devoted to the problem of determining a resource of technical objects, which boils down to finding the time remaining before their limiting state, non-destructive method for the analysis of acoustic emission (AE) from a position of micro-mechanical models of destruction and time parameters of AE.
NON-DESTRUCTIVE TESTING. ACOUSTIC EMISSION. HOMOGENEOUS DESTRUCTION. RESOURCE. MICROMECHANICAL MODEL.
Надежность технических объектов зависит от точности оценки их ресурса. Оценка ресурса технического объекта сводится к нахождению времени, оставшегося до наступления его предельного состояния. Существует большое количество методов разрушающего и неразрушаю-щего контроля для решения данной проблемы. К недостаткам существующих методик контроля разрушения можно отнести:
низкую точность оценки, поскольку методики основаны на неадекватных моделях (дают качественную оценку; занижают значение остаточного ресурса при неучете неоднородности состояния), статистически неустойчивы и имеют повышенную чувствительность к помехам (не учитывают изменения условий регистрации сигналов АЭ во время контроля, завышают значение ресурса в случае их отсутствия);
трудоемкость, так как требуют проведения обучающего эксперимента и определения регрессионных зависимостей вследствие недостаточной физической обоснованности и базирования на статистических моделях;
они не дают долгосрочного прогноза остаточного ресурса, поскольку основаны на регистрации только заключительной стадии разрушения (следствие пренебрежения микромеханическими аспектами разрушения и АЭ-излучения).
В предлагаемой статье рассматривается метод нахождения времени до разрушения на основании анализа сигналов акустической эмиссии, полученных в результате испытаний опытных образцов. Методика обладает повышенной точностью. Предложенные АЭ-показатели тесно связаны с определяющим прочность процессом, статистически устойчивы к помехам регистрации. Снижена неопределенность в оценке класса источника сигналов АЭ. Методика проста в применении, реализуется без обучающего эксперимента, универсальна и унифицирована, так как диагностические параметры позволяют оценить ресурс и стандартные характеристики прочности различных по структуре и размерам объектов. Главное достоинство данного метода — возможность оценки ресурса объекта на основании информации, полученной неразру-шающим методом на месте эксплуатации. Ис-
пользуемые физическая и математическая модели не зависят от множества внешних факторов, — это упрощает модель, делает ее наиболее точной.
Описание задачи
Основная задача — достичь максимально возможной точности определения времени до разрушения. Программа, разработанная применительно к данной задаче, использует для аппроксимации этап однородного разрушения, определение которого способствует увеличению точности метода.
Одна из основных причин нарушения работоспособности и уменьшения ресурса конструкций и деталей машин — развитие в их материале процесса накопления повреждений. Для описания этого процесса используется микромеханическая модель разрушения, согласно которой временная зависимость по-врежденности материала в условиях мелкодисперсного роста концентрации микротрещин С имеет следующий вид:
С (/) = Со
1 - ехр
-I
Ж
©»С)
0 ^ср
(1)
где ^ — текущее время; С0 — начальная концентрация структурных элементов в материале до разрушения; 0ср — среднестатистическое время ожидания разрушения одного структурного элемента, задаваемое формулой Журкова, при начальном условии С(0) = 0.
Формула Журкова в условиях растяжения при постоянстве напряжения имеет вид
ср
= т0 ехр
и
0
ЯТ
(2)
состояния материала вводится функция Т(у) распределения параметра у по структурным элементам. В качестве функции Т(у) в зависимости от степени неоднородности материала используются различные распределения. В нашей работе при описании разрушения образцов сварных соединений рассматриваются два из них (рис. 1):
ограниченное вейбулловское распределение с параметрами к, X, д
¥(у) =
0, X е [ 0, д),
к-1
4
ехр
к
У
(3)
А = -
Л XI X
к-1
ехр
й у
логарифмически-нормальное распределение с параметрами стг, ц
1 11111 и I - I ^2
у( У) = ~ГТ=— ехр
л/2П(
яст 2 у
(1п( 7) -I)2
2ст2
(4)
-■Фп
ф ф
где т0 « 10-13, с, — период атомных колебаний; и0 — энергия активации процесса разрушения; у — структурно-чувствительный коэффициент; Я — универсальная газовая постоянная; Т — абсолютная температура.
Из-за неоднородности структурного и напряженного состояний материала среднестатистические значения времени разрушения структурных элементов тела в общем случае неодинаковы. Таким образом, для различных элементов параметр © принимает различные значения. Для характеристики неоднородности
\ф ф\
-1
Рис. 1. Схема нагружения образца:
1 — образец; 2 — сварной шов; 3 — верхний захват нагружающего устройства; 4 — нижний захват нагружающего устройства; 5 — палец; 6 — ПАЭ
1
д
5
4
Зависимость концентрации микротрещин от времени для неоднородного материала можно представить в виде
Ymax
C(t) = C0 Jv(Y)
1 - exp
dt
' т0 exp
Uо -YP(t ) RT
d y,
(5)
где ут1п и утах — минимальное и максимальное значения у для данного распределения Т(у).
Переход к макроскопическому разрушению материала происходит, когда число разрушенных структурных элементов достигает 0,01 от
числа начальных: С/С0 = 0,01. Время до разру-
*
шения т находится из уравнения
max
J ¥(Y)
1 - exp
-J-
dt
0 , U0 — Y°'(t ) 0 т0 expl —0— 0 1 RT
dY = 0,01. (6)
Для вычисления времени до разрушения технического объекта необходимо определить значения величин, входящих в уравнение (6). Необходимо знать зависимость напряжения от времени, температуру, величину энергии активации и0 (постоянна для конкретного материала) и параметры функции Т(у). Для нахождения последних используется регистрация сигналов АЭ при нагружении объекта контроля в области упругих деформаций. Число сигналов АЭ N¡(1) пропорционально числу разрушенных структурных элементов:
N¡(0 = кАЕт,
или
^ (1) =
= kAEC0 J VM
1 - exp
ЛЛ
-J-
dt
1 т0 exp
U0 -Yo(t ) RT
d Y,
(7)
где кАЕ — акустико-эмиссионный коэффициент (может быть принят постоянным при нагруже-нии с постоянной скоростью в области упругих деформаций).
Параметры функции Т(у) находятся путем аппроксимации экспериментальной зависимости числа импульсов АЭ от времени теоретической кривой (7). Далее определяется время до разрушения технического объекта из уравнения (6).
Описание эксперимента
В работе предполагалось определение параметров функции Y(y) для образцов стыковых сварных соединений. Материал пластин — сталь СтЗпс; сварка — ручная и автоматическая в среде углекислого газа. Образцы испытывались на растяжение с постоянной скоростью роста деформации, при этом регистрировались сигналы АЭ. Нагружение образцов производилось следующим образом (рис. 1). Концы образца 1 закреплялись в захватах 3 и 4 испытательной установки посредством пальцев 5. Верхний захват
3 оставался неподвижным, а нижний захват
4 перемещался вертикально вниз с постоянной скоростью. Сила, приложенная к образцу, определялась по шкале установки (цена деления шкалы 3000 Н). Для регистрации сигналов АЭ была использована двухканальная акустико-эмис-сионная система СДАЭ-16(2), которая состояла из двух преобразователей акустической эмиссии (ПАЭ) типа П113с, предварительных усилителей, блока регулирования усиления, блоков обработки сигналов АЭ, блока управления и компьютера. Принцип работы одного ее канала следующий: акустический сигнал, возникающий в материале образца при его нагружении, превращается пьезоэлектрическим преобразователем в электрический импульс, который подвергается усилению и обработке, после чего информация поступает в компьютер, где отображается в виде таблицы данных и графиков зависимости от времени различных параметров АЭ сигнала. Регистрировались число импульсов АЭ, их амплитуда, энергия, длительность, суммарный счет АЭ и время нарастания сигнала.
Рабочий диапазон частот аппаратуры составлял 20—1000 кГц, эффективное значение напряжений собственных шумов усилительного тракта, приведенное ко входу предварительного усилителя, не превышало 5 мкВ. Нижний уровень порога дискриминации системы определялся необходимостью устранения электромагнитных помех и составлял 47 дБ на каждом канале. Для регистрации сигналов были выбраны преобразователи с диапазоном частот 20—500 кГц.
Для обработки первичных данных АЭ испытаний, определения коэффициентов микромеханической модели разрушения предназначена программа RG, состоящая из нескольких макросов Microsoft Excel.
Y
Y
min
Y
Одним из этапов выполнения работы является нахождение интервала для аппроксимации. Требуется выбрать линейный участок зависимости логарифма числа импульсов от времени (этап однородного разрушения), а также разработать методологию нахождения данного участка. Были произведены испытания 30-ти образцов и проанализированы: временная зависимость числа импульсов, логарифм числа импульсов, длительность импульсов, коэффициент перекрытия импульсов (вычисляется как отношение суммарной длительности импульсов в течение интервала времени к величине этого интервала), амплитуда, коэффициент вариации амплитуды, параметр Кт (отношение активности АЭ к средней длительности импульсов). В результате анализа и выявления зависимостей выделены следующие задачи этапов нахождения начальной и конечной точек прямолинейного участка:
а) найти на графике зависимости длительности импульсов АЭ от времени момент 1Ъ когда наблюдается резкий рост длительности импульсов (свидетельствует о начале пластического деформирования);
б) найти на графике зависимости коэффициента перекрытия от времени момент ¿2, соответствующий началу роста коэффициента перекрытия. Он свидетельствует о начале пластического деформирования, при котором происходит движение дислокаций, мелкие перестройки в структуре материала, а значит, большое количество «коротких» сигналов, которые накладываются друг на друга. Они сливаются, что
отражается в виде скачка на графике зависимости коэффициента перекрытия от времени;
в) из двух значений времени (^ и ¿2) выбрать меньшее.
Пример определения интервала времени для аппроксимации. Для одного из образцов стыковых сварных соединений:
а) резкий рост длительности импульсов наблюдается в момент времени ^ = 186 с (см. рис. 2).
б) коэффициент перекрытия импульсов начинает расти в момент времени = 175 с (рис. 3).
в) Ь < — для проведения дальнейших действий используем величину времени ¿2. Изменение характера сигналов АЭ указывает на появление пластической деформации материала образца в момент времени = 175 с.
Далее требуется найти первоначальную точку линейного участка (начало упругой деформации), для чего рассматриваем график зависимости коэффициента вариации амплитуды от времени. На участке от 0 до ¿2 (двигаемся от 12 к 0) выделяем момент в который происходит аномальный скачок коэффициента вариации амплитуды.
Скачок вариации амплитуд объясняется тем, что происходит разрушение связей с большим разбросом размеров, структур и напряжений. Разрушение неоднородно, а следовательно, возникают аномально большие амплитуды.
Полученный интервал вносится в макрос программы. Макрос строит теоретическую кривую зависимости числа импульсов АЭ от времени, используя в качестве функции ¥(у) ограниченное вейбулловское распределение.
Длительность импульса, мкс
12000
8000 6000 4000 2000 0
= 18 >6 с ►
>6 с Л ♦ .А ♦
♦
^ ♦ < ►
► V} № 1 ♦
1+
0
100
200
300
400
500
600
700 Время, с
Рис. 2. Определение момента времени ^ по графику зависимости длительности импульсов от времени
Коэффициент перекрытия
0,035 0,030 0,025 0,020 0,015 0,010 0,005 0
0
100
200
t
t2 = 17: м i A
№ 1
Л A, M & t
♦t • J In fr i
JL и ■ LU A П \
I^V mul
шшт Wr -w m—
300
400
500
600
700
Время, с
Рис. 3. Определение момента времени 12 по графику зависимости коэффициента перекрытия импульсов от времени
Изменяя величину параметра Я, мы можем проследить за изменением среднего значения параметра у распределения ¥(у).
Путем изменения параметров Я, # и кАЕС0 (к = 3 — одинаков для всех образцов) добиваемся совпадения теоретической и экспериментальной кривых N¡(1). При этом величина «Средней ошибки аппроксимации» должна быть минимальной. При совпадении кривых получаем искомые значения параметров микромеханической модели разрушения.
Таким образом, нами был получен метод определения участка однородного разрушения; с его помощью увеличена точность определения
времени до разрушения объекта на 2 %. Сократилась погрешность (с 15 до 13 %), данный результат в пределах нормы. Мы убедились, что разработанная методика обладает повышенной точностью. На основе АЭ-показателей, напрямую связанных с определяющим прочность процессом и статистически устойчивых к помехам регистрации, получены файлы MS Excel, представляющие собой базу, позволяющую выполнять различные вычисления по желанию пользователя. Мы убедились в простоте применения и универсальности методики, так как диагностические параметры позволяют оценить ресурс и стандартные характеристики прочности различных по структуре и размерам объектов.
Коэффициент вариации амплитуды
0,100 0,080 0,060 0,040 0,020 0,000
0 100 200 300 400 500 600 700 Время, с
Рис. 4. Определение момента времени по графику зависимости коэффициента вариации амплитуды от времени
СПИСОК ЛИТЕРАТУРЫ
1. Носов, В.В. Методология оценки прочности и ресурса технических объектов с помощью метода АЭ [Текст] / В.В. Носов // Труды СПбГТУ. № 504.— 2007.— С. 119-132.
2. Лахова, Е.Н. Моделирование процесса перестройки структуры при пластической деформации сварных соединений металлоконструкций [Текст] / Е.Н. Лахова, В.В. Носов // Современное машиностроение. Наука и образование: матер. Междунар. научно-практ. конф. 14-15 июня 2011 года, Санкт-Петербург.— СПб.: Изд-во Политехи. ун-та, 2011.— С. 301-309.
3. Регель, В.Р. Кинетическая природа прочности
твердых тел [Текст] / В.Р. Регель, А.И. Слуцкер, Э.Е. Томашевский.— М.: Наука, 1974.
4. Носов, В.В. Микромеханическая модель акустической эмиссии гетерогенных материалов [Текст] / В.В. Носов, И.Н. Бураков // Дефектоскопия.— 2004. № 2.— С. 53-61.
5. Бураков, И.Н. Оценка параметров процесса разрушения на основе компьютерного моделирования акустической эмиссии сварных соединений [Текст] / И.Н. Бураков // Моделирование акустической эмиссии гетерогенных материалов: Матер. меж-дун. науч. конф.— СПб.: Изд-во СПбГПУ, 2004.— С. 38-39.
АХМЕТШИН ТЬмур Дамирович — студент VI курса Санкт-Петербургского государственного политехнического университета.
195251, Политехническая ул., д. 29, Санкт-Петербург, Россия [email protected]
НОСОВ Виктор Владимирович — доктор технических наук, профессор кафедры машиноведения и основ конструирования Санкт-Петербургского государственного политехнического университета.
195251, ул. Политехническая, д. 29, Санкт-Петербург, Россия [email protected]
© Санкт-Петербургский государственный политехнический университет, 2013