CC BY
48
УДК 004.94:621.78
И.Е. РОГОВ, А.А. ЛАВРЕНТЬЕВ, С.Д. АНИСИМОВ
ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ В СИСТЕМАХ АВТОМАТИЗАЦИИ ДИАГНОСТИКИ ИЗДЕЛИЙ МАШИНОСТРОЕНИЯ
Описаны математическая модель образования сигнала токовихревого датчика и имитационная модель создания виртуальных образцов для настройки электромагнитного прибора неразрушающего контроля. Ключевые слова: моделирование, математическая модель, имитационная модель, неразрушающий контроль, токовихревой датчик.
Введение. Важные детали машин и механизмов, испытывающие значительные нагрузки, изготавливаются из сталей, и их поверхность упрочняется такими технологическими операциями, как поверхностная закалка токами высокой частоты, химико-термическое упрочнение, плазменнодиффузионное, лазерная обработка, упрочнение наклепом, нанообработка и др. Для обеспечения высококачественного производства необходимо иметь различные средства неразрушающего (в некоторых отраслях - стопроцентного) контроля, которые позволят не только выявлять и изолировать брак, но и проводить диагностику причин появления брака для корректировки технологического процесса.
Применяемые методы настройки приборов для неразрушающего контроля, использующие электромагнитные датчики, связаны с изготовлением и экспериментальным исследованием настроечных образцов с известными (а зачастую и заданными) свойствами. Причем количество этих образцов должно быть статистически значимым (иногда более 100 штук). Изготовление такого количества образцов и исследование их физико-механических свойств требует как больших финансовых затрат, так и затрат времени, что делает невозможной быструю переналадку системы контроля на новую задачу.
Проблему можно решить, если применить компьютерное моделирование [1]: математическая модель описывает связь между сигналом датчика и физико-механическими свойствами изделия, а имитационная модель позволяет создавать «виртуальные» образцы изделия с требуемыми свойствами. Таким образом, методом компьютерного эксперимента возможно определять оптимальные условия контроля и находить значения коэффициентов для настройки приборов неразрушающего контроля.
Физико-математическая модель образования сигнала датчика. В приборах неразрушающего контроля в качестве накладного электромагнитного преобразователя часто используется датчик, показанный на рис. 1.
Изделие с переменными по толщине (координате Х) свойствами представляется состоящим из n тонких слоев, имеющих неизменные электромагнитные свойства: толщину h, комплексную магнитную проницаемость х, электропроводность о, коэффициент Рэлея %. Датчик состоит
из ферромагнитного сердечника с намотанными на нем возбуждающей обмоткой W1 и измерительной обмоткой W2. Он расположен над контролируемым изделием с зазором 5 [2].
Электромагнитное поле в проводящей среде описывается уравнениями Максвелла:
rotН = oE, rotE = - ja хН .
Для схемы расположения датчика над контролируемой средой можно записать граничные условия:
= о, дЛ = о, hz = о.
dZ dZ Z
Рис.1. Электромагнитный датчик, расположенный над многослойной контролируемой средой
Следовательно,
Е_ = 0, Е,
Получаем уравнение:
дЙ дНл
ду дх
стЕ,
д_К
ду
= Нх
д2Ег д2 Ё2 . •
2 +----^ - ]®\мзЕг = 0 ,
ду
дх
которое является однородным уравнением Гельмгольца в декартовой системе координат:
д 2Е д 2 Е„
'*■ + т 2 Е,
0
где т
Дальнейшее решение уравнения позволяет определить ЭДС датчика [3]:
\
k л ~2
Е = -Ж1Ж2
7
у т1
1 1 А + С2(т) 4 ^ 1 А + С7
+-------------=^-
2 в+Аі_то) ■ л2 й k2 в+аіт)
Полученное уравнение можно использовать для определения сигнала электромагнитного датчика в имитационной модели. Для примера приведено сравнение моделированных и измеренных на образцах изделий сигналов датчика, из которого хорошо видно практически полное их совпадение (рис.2).
0
Рис.2. Сравнение амплитудно-частотных характеристик моделированного и измеренного на образцах изделий сигналов при различных относительных зазорах между датчиком и изделием
В ряде случаев, применение основных колебаний для электромагнитного контроля приводит к недопустимой погрешности. Например, ЭДС датчика в зависимости от толщины упрочненного слоя определяется, главным образом, отражением электромагнитной волны от поверхности, разделяющей упрочненный слой и основу. Упрощенно зависимость ЭДС от толщины упрочненного слоя в этом случае определяется коэффициентом отражения:
Ц 2С1
-1,
(*)
где ц1 и о1 соответствуют упрочненному слою, а ц2 и о2 - основе.
При всех известных видах упрочнения, вызывающих физические дефекты в идеальном кристалле железа, магнитная проницаемость Ц и удельная электропроводность с1 упрочненного слоя уменьшаются по сравнению с магнитной проницаемостью Ц 2 и удельной электропроводностью с2 основной стали. Если ц и с изменяются пропорционально, то выражение (*) стремится к нулю, и преобразователь теряет чувствительность к толщине упрочненного слоя. В действительности выражение (*) не равно нулю, но близко к нему. Сигнал отраженной волны на границе раздела слоев настолько мал, что его превосходят помехи, вызванные вариациями электромагнитных свойств изделия. Это и является причиной возникновения недопустимо большой ошибки при контроле некоторых видов упрочнения по основным гармоникам сигнала преобразователя.
Исследование зависимости информативности комбинационных колебаний по сравнению с основными показало, что в ряде случаев комбинационные колебания значительно информативнее. Несмотря на то, что на комбинационные колебания сильно влияет зазор между датчиком и контролируемым изделием, влияние толщины упрочненного слоя на них гораздо сильнее. Поэтому отношение полезного сигнала к помехе в них больше, чем в основных гармониках. Использование комбинационных колебаний повышает точность контроля толщины упрочненного слоя.
Например, при упрочнении борированием коэффициент отражения [см. формулу (*)] равен 0,058, тогда как аналогичный коэффициент для комбинационных колебаний равен 0,44 - это почти в десять раз выше.
Ц 1с2
ЭДС комбинационных колебаний в слабых полях рассчитывают, учитывая коэффициент Рэлея х и то, что электромагнитные свойства слоев изделия изменяются по глубине (см. рис.1).
Решение двумерного дифференциального уравнения Гельмгольца для каждого / -го слоя модели имеет вид:
Нуг =— [4 ехР(- А** х) - Сгк ехР( Да? х)]соэ ^У .
Ю Ц ,■ к
Граничные условия между слоями г и г +1 можно представить выражениями:
А,(х) г - - и
— [ Ак еХР(-АкX ) - Сгк ехр(X ) ] =
Ц 1 J
= [Д.+и ехр(-) - сг+и ехр() ],
ЦI+1
Агк еХР(- АХк X ) + Сгк еХР(Дк'1х. ) =
= А +1,к еХР(-Ай-1, кхг ) + Сг+1,к еХР( Ай1 кхг )-
В последней (п+1)-й области металла (т.е. на бесконечно большой глубине) электромагнитное поле полностью затухает, поэтому при х Нп+1у ^ 0 и Нп+1х ^ 0. Такое возможно
при Ап+1,к = 0 .
В модели выбирается большое количество настолько тонких слоев, что напряженность магнитного поля внутри каждого слоя изменяется менее чем на 1%. Большим количеством слоев достаточно точно дискретно отражается непрерывное распределение электромагнитных свойств упрочненного слоя и основного металла.
Напряженность поля в указанных точках определяется как функция времени:
Нуг = Ншг + Уг ) .
Приведенная последовательность повторяется дважды для двух частот ю1 и ю2 основных
колебаний тока возбуждения преобразователя, которые в нелинейной среде, взаимодействуя, образуют комбинационные колебания.
Чтобы комбинационные колебания отличались по частоте от основных колебаний и их высших гармоник, основные частоты должны относиться как достаточно большие целые взаимно простые числа N и М:
ю1 :ю2 = N:М .
Напряженность результирующего поля в точке хг:
Нг = Нш,1,г 8Ш(Ю^ + Уи ) + Нш,2,г + V2.)-
Получившееся колебание - периодическое; его период составляет
Т = — N = — М .
Ю Ю 2
Непрерывно изменяющаяся во времени напряженность вычисляется дискретно через интервалы времени Аг1. Если ограничить спектр сигнала суммарными комбинационными колебаниями третьего порядка, то, по теореме Котельникова, интервал времени
Т
А ? = -
2 N + 4М
Индукция магнитного поля в точке хг в момент времени tp определяется уравнением
Рэлея:
вгр = A iH ip +i,HfpsignH,p.
Рабочий поток магнитной индукции в контролируемом изделии рассчитывается заменой интеграла суммой:
* и+1
ф р = b\Bpdx * bZ Biphi >
0 i=1
где Фp есть дискретная функция времени.
Гармоническим анализом переменного магнитного потока можно определить вещественную и мнимую составляющие потока суммарных комбинационных колебаний третьего порядка:
2 1 2 m 9^
Фком = — p(0sin(2<B2 +<Bi)tdt * — ^Ф sin—(2M + N)p,
1 J m , m
о p=1
2 1 2 m 2^
Ф;Ком = — |O(0cos(2<i>2 +&i)tdt * — ^Фp cos—(2M + N)p,
1 J m , m
о p=1
где m - число интервалов времени A t, на которые разбит период.
Комплексная амплитуда рабочего потока:
Ф = Ф + /Ф1
ком ком J ком
частично сцеплена со вторичной обмоткой преобразователя, поэтому комплексная амплитуда комбинационных суммарных колебаний третьего порядка рассчитывается следующим образом [4]:
Э = - /(в, + 2ю 2) W2 Ф' .
ком J V 1 2 ' 2 ком
Подобным образом с помощью гармонического анализа магнитного потока в математической модели вычисляются ЭДС основных колебаний и разностных комбинационных колебаний
третьего порядка. Принципиально возможно и вычисление высших гармоник.
Имитационная модель создания виртуальных образцов. Предназначена для имитации виртуального прибора для неразрушающего контроля, виртуальных образцов, ЭДС сигналов измерительного преобразователя (датчика), вычисления линейной регрессионной модели по сигналам основных и комбинационных колебаний, а также определения ошибки контроля на виртуальных проверочных образцах.
Имитационная модель реализована в виде программы для ЭВМ. Структурная схема имитационной модели показана на рис.3. При моделировании вычисляются ЭДС датчика для набора значений либо зазора, либо толщины упрочненного слоя. Кроме того, при вычислении ЭДС происходит случайное варьирование электромагнитных параметров виртуального образца изделия, а также случайное варьирование того параметра, который не варьируется детерминирован-но. В конце моделирования к полученным значениям ЭДС добавляются случайные значения, имитирующие аддитивную и мультипликативную погрешности прибора для неразрушающего контроля. Такие вычисления ЭДС датчика производятся для заданного ряда частот как основных, так и комбинационных колебаний.
На рис.4 показана структурная схема блока вычисления ЭДС сигналов измерительного преобразователя при детерминированном варьировании зазора между датчиком и изделием. Аналогичный вычислительный блок используется при детерминированном варьировании толщины упрочненного слоя, тогда производится цикл не по значениям зазоров, а по значениям толщин упрочненного слоя.
Рис.3. Структурная схема имитационной модели
В системе создается два набора виртуальных образцов - настроечные и проверочные. Настроечные образцы используются для построения регрессионной модели. После получения коэффициентов уравнения регрессии, в это уравнение подставляют значения сигналов, полученные на проверочных образцах. То есть в этом случае имитируется процесс измерения контролируемого параметра виртуальным прибором для неразрушающего контроля. Результат этой операции -показания прибора, получающиеся при измерении образца, не участвовавшего в настройке. Оба набора виртуальных образцов могут различаться между собой.
Рис. 4. Структурная схема вычислительного блока ЭДС сигналов измерительного преобразователя
При моделировании возможно задавать различные начальные условия для массивов настроечных и проверочных образцов. Например, количество настроечных образцов можно задавать небольшим, чтобы оценить представительность выборки настроечных образцов в данных условиях контроля. Различной может быть и величина случайных вариаций электромагнитных свойств в каждом из наборов. Тогда появляется возможность определить, как повлияет ошибка в оценке параметров настроечных образцов, или что покажет прибор при контроле изделий из другой партии, отличающихся либо значениями технологических факторов, либо их разбросом.
Ошибка контроля зависит от количества настроечных образцов при различном количестве частот, составляющих спектр измерительных сигналов. С ростом количества настроечных образцов ошибка сначала падает, а потом перестает меняться. Кроме того, ошибка контроля зависит и от числа частот, составляющих спектр сигнала датчика. При малом количестве настроечных образцов с ростом числа частот, составляющих спектр, включенный в регрессионную модель, растет и ошибка контроля. Это происходит из-за снижения количества степеней свободы при построении модели и, следовательно, снижения представительности выборки сигналов.
Имитируемый образец разбивается на п слоев с неизменными электромагнитными свойствами внутри слоя. Магнитный поток в изделии вычисляется интегрированием по толщине изделия. Поэтому толщина одного слоя разбиения должна быть достаточно мала, чтобы изменение величины напряженности поля при переходе от слоя к слою также было малым. Однако при слишком большом числе слоев резко растет вычислительная ошибка и продолжительность вычислений, следовательно, требуется оптимизация числа слоев для минимизации погрешности моделирования.
Ошибка контроля зависит и от числа слоев, на которые разбивается виртуальный образец. При количестве разбиений более 1000, ошибка контроля начинает расти. Это происходит из-за накопления погрешностей округления при вычислениях с конечной точностью [5].
—■—/7=200 мкм —* /г =400 мкм —•—/7=600 мкм
—■—/г=800 мкм =1000 мкм —•—/7=2000 мкм
Рис.5. Частотные характеристики модуля передаточного коэффициента толщины упрочненного слоя для комбинационных колебаний
Для обеспечения высокой точности измерения контролируемого параметра должен быть установлен оптимальный спектр частот возбуждения преобразователя. Амплитудно-частотные характеристики передаточного коэффициента толщины упрочненного слоя (рис.5) позволили установить области частот возбуждения датчика, при которых он высокочувствителен к контролируемому параметру, и выбрать, таким образом, составляющие многомерного сигнала. Заключение. Разработана математическая модель образования сигнала электромагнитного датчика, расположенного над многослойной проводящей ферромагнитной средой для основных и комбинационных колебаний. Несмотря на то, что амплитуды комбинационных колебаний намного меньше, чем основных, они оказываются более информативными. Используя комбинационные колебания, можно существенно повысить точность неразрушающего контроля. Поскольку при упрочнении значительно изменяется коэффициент Рэлея, то комбинационные колебания оказываются намного чувствительнее к толщине упрочненного слоя, чем основные.
Разработана имитационная модель образования сигнала электромагнитного датчика для основных и комбинационных колебаний. Поскольку аналитическое выражение для комбинационных колебаний не очень удобно для прямого расчета, то предложен рекуррентный метод, используя который можно достаточно просто получить значения ЭДС датчика как для основных, так и для комбинационных колебаний (а также для высших гармоник, если такое понадобится) любого порядка. Генерация виртуальных образцов и их использование для построения и проверки точности регрессионной модели позволяет определить оптимальные условия контроля без изготовления и исследования реальных образцов изделий. Кроме того, значительно сокращается количест-
во реальных образцов изделий, требуемое для окончательной настройки прибора при его эксплуатации.
Усовершенствованы методы формирования и обработки сигнала датчиков в целях повышения точности контроля. Использование имитационной модели позволяет определить оптимальный спектр возбуждения датчика и параметры регрессионной модели для обработки его сигнала. Выбор оптимального спектра возбуждения значительно повышает точность контроля и уменьшает влияние на его результаты случайных вариаций свойств изделия, а также различных помех, воздействующих на устройство в целом.
Библиографический список
1. Рыжиков Ю.И. Имитационное моделирование: теория и технологии / Ю.И. Рыжиков. -М.: Альтекс, 2004. - 384 с.
2. Рогов И.Е. Модель электромагнитного датчика с немагнитной проводящей вставкой / И.Е. Рогов // Диагностика и управление в технических системах: межвуз. сб. науч. ст. /Донск. гос. техн. ун-т. - Ростов н/Д: Издательский центр ДГТУ, 1998. - С.92-96.
3. Электромагнитный неразрушающий контроль упрочненного слоя на поверхности стальных изделий / С.Д. Анисимов [и др.] // Дефектоскопия. - 2000. - №1. - С.18-30.
4. Создание математической модели и программы ЭВМ для расчета комбинационных колебаний сигналов электромагнитных датчиков: отчет о НИР по гранту 49Гр-98 / Донск. гос. техн. ун-т; рук. Анисимов С.Д.; исполн.: И.Е. Рогов [и др.]. - Ростов н/Д, 2001. - №ГР01980006915. -Инв. №02.20.0300133.
5. Оптимизация модели упрочненного слоя при поверхностном упрочнении / И.Е. Рогов [и др.] // Вестн. Донск. гос. техн. ун-та. - 2009. - Т.9, №2(41). - С.208-216.
Материал поступил в редакцию 03.11.10.
References
1. Ryjikov Yu.I. Imitacionnoe modelirovanie: teoriya i tehnologii / Yu.I. Ryjikov. - M.: Al'teks, 2004. - 384 s. - In Russian.
2. Rogov I.E. Model' elektromagnitnogo datchika s nemagnitnoi provodyaschei vstavkoi / I.E. Rogov // Diagnostika i upravlenie v tehnicheskih sistemah: mejvuz. sb. nauch. st. /Donsk. gos. tehn. un-t. - Rostov n/D: Izdatel'skii centr DGTU, 1998. - S.92-96. - In Russian.
3. Elektromagnitnyi nerazrushayuschii kontrol' uprochnennogo sloya na poverhnosti stal'nyh izdelii / S.D. Anisimov [i dr.] // Defektoskopiya. - 2000. - №1. - S.18-30. - In Russian.
4. Sozdanie matematicheskoi modeli i programmy EVM dlya rascheta kombinacionnyh kolebanii signalov elektromagnitnyh datchikov: otchet o NIR po grantu 49Gr-98 / Donsk. gos. tehn. un-t; ruk. Anisimov S.D.; ispoln.: I.E. Rogov [i dr.]. - Rostov n/D, 2001. - №GR01980006915. -Inv. №02.20.0300133. - In Russian.
5. Optimizaciya modeli uprochnennogo sloya pri poverhnostnom uprochnenii / I.E. Rogov [i dr.] // Vestn. Donsk. gos. tehn. un-ta. - 2009. - T.9, №2(41). - S.208-216. - In Russian.
I.E. ROGOV, A.A. LAVRENTYEV, S.D. ANISIMOV
SIMULATION MODELING IN AUTOMATION AND DIAGNOSTIC SYSTEMS OF ENGINEERING PRODUCTS
A mathematical model of generating an eddy current probe signal and a simulation model of creating virtual prototypes for setting the electromagnetic device of nondestructive testing are described.
Key words: modeling, mathematical mode, simulation model, nondestructive testing, eddy current probe.
РОГОВ Игорь Евгеньевич (р. 1964), старший преподаватель кафедры «Электротехника и электроника» Донского государственного технического университета. Окончил Ростовский-на-Дону институт сельскохозяйственного машиностроения (1986).
Область научных интересов: методы неразрушающего электромагнитного контроля, имитационное моделирование в структуроскопии.
Автор более 20 научных публикаций. Имеет 1 патент.
ЛАВРЕНТЬЕВ Анатолий Александрович (р. 1952), заведующий кафедрой «Электротехника и электроника» Донского государственного технического университета, доктор физикоматематических наук (2001), профессор (2003). Окончил Ростовский государственный университет (1974).
Область научных интересов: моделирование электронно-энергетической структуры сложных полупроводниковых соединений и некоторых металлических сплавов, имитационное моделирование в структуроскопии.
Автор более 150 научных публикаций. [email protected]
АНИСИМОВ Сергей Дмитриевич (р. 1935), доцент кафедры «Электротехника и электроника» Донского государственного технического университета, кандидат технических наук (1966), доцент (1971). Окончил Таганрогский радиотехнический институт (1959).
Область научных интересов: методы неразрушающего электромагнитного контроля.
Автор 120 научных публикаций. Имеет 1 патент и 7 авторских свидетельств.
Igor E. ROGOV (1964), Senior Lecturer of the Electrical and Electronics Engineering Department, Don State Technical University. He graduated from Rostov-on-Don Institute of Agricultural Engineering (1986).
Research interests: nondestructive electromagnetic testing methods, imitating modeling in structural analysis.
Author of over 20 publications and 1 invention.
Anatoly A. LAVRENTYEV (1952), Head of the Electrical and Electronics Engineering Department, Don State Technical University. PhD in Physics and Maths (2001), Professor (2003). He graduated from Rostov State University (1974).
Research interests: modeling of compound semiconductors and metal alloys electronic-power structures, simulation modeling in structural analysis.
Author of over 150 publications.
Sergey D. ANISIMOV (1935), Associate Professor of the Electrical and Electronics Engineering Department, Don State Technical University. Candidate of Science in Engineering (1966), Associate Professor (1971). He graduated from Taganrog Radio-Engineering Institute (1959).
Research interests: nondestructive electromagnetic testing methods.
Author of 120 publications, 1 patent and 7 certificates of authorship.
