автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Математические модели и алгоритмы компьютерной обработки цифровых изображений реакции датчиков деформации интегрального типа

кандидата технических наук
Змызгова, Татьяна Рудольфовна
город
Курган
год
2007
специальность ВАК РФ
05.13.18
Диссертация по информатике, вычислительной технике и управлению на тему «Математические модели и алгоритмы компьютерной обработки цифровых изображений реакции датчиков деформации интегрального типа»

Автореферат диссертации по теме "Математические модели и алгоритмы компьютерной обработки цифровых изображений реакции датчиков деформации интегрального типа"

На правах рукописи

Змызгова Татьяна Рудольфовна

МАТЕМАТИЧЕСКИЕ МОДЕЛИ И АЛГОРИТМЫ КОМПЬЮТЕРНОЙ ОБРАБОТКИ ЦИФРОВЫХ ИЗОБРАЖЕНИЙ РЕАКЦИИ ДАТЧИКОВ ДЕФОРМАЦИИ ИНТЕГРАЛЬНОГО ТИПА

05 13 18 — математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

1 003162281

Тюмень-2007

003162281

Работа вьтполнена на кафедре технологии и автоматизации сварочного производства ГОУ ВПО Курганский государственный университет

Научный руководитель заслуженный деятель науки РФ,

доктор технических наук, профессор Сызранцев Владимир Николаевич

Официальные оппоненты доктор технических наук, профессор

Шапцев Валерий Алексеевич

доктор физико-математических наук, профессор

Борзых Владимир Эрнестович

Ведущая организация Институт механики ИжГТУ (г Ижевск)

Защита диссертации состоится 13 ноября 2007 года в 16 часов на заседании диссертационного совета К21227401 при Тюменском государственном университете по адресу 625003, г Тюмень, ул Перекопская, 15А, ауд 217

С диссертацией можно ознакомиться в библиотеке Тюменского государственного университета

Автореферат разослан «У » октября 2007 года

диссертационного совета

Ученый секретарь

Бутакова Н Н

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы. Один из новых подходов к диагностике усталости различных деталей, узлов и конструкций машин базируется на обработке информации, полученной при помощи датчиков деформации интегрального типа (ДДИТ), которые представляют собой новые средства экспериментального исследования Существенный вклад в развитие данного научного направления внесли В Н Сызранцев, А С Терехов, А Ю Розенберг, Д А Троценко, А Ю Удовикин, С Л Голофаст, О А Штин и др

Интенсификация научных исследований требует широкого применения и развития существующих методов хранения, обработки и анализа экспериментальной информации, полученной при помощи ДДИТ, на базе современных вычислительных средств Осуществление сложных методов обработки изображений реакции ДДИТ в режиме реального времени требует применения эффективных базовых моделей (фильтрации, идентификации), экономично использующих соответствующие информационные ресурсы

Таким образом, в настоящее время существует достаточно важная научно-техническая проблема разработки новых моделей и устройств формализации и анализа цифровых изображений реакции ДДИТ в рамках единой автоматизированной системы обработки изображений их реакции

Целью работы является синтез математических моделей и алгоритмов эффективной формализации изображений реакции ДДИТ, которые позволят решать задачи оценки нагруженности и прогнозирования усталостной долговечности деталей и конструкций машин при стендовых и эксплуатационных испытаниях

Для достижения поставленной цели в диссертации автором были поставлены следующие задачи

1 Рассмотреть различные модели шумовых воздействий и провести сравнительный анализ известных моделей фильтрации для предварительной обработки изображений реакции ДДИТ

з

2 Разработать новые модели обработки показаний ДДИТ, адаптированные к особенностям исследуемых изображений, которые обеспечат построение эффективного формального описания и калибровку ДДИТ по относительной площади «темных пятен»

3 Разработать новую модель автоматической сегментации информативных объектов изображений реакции ДДИТ, исследовать ее эффективность и точность

4 Разработать новые способы линеаризации бинарных изображений реакции ДДИТ (скелетизация, выделение контура)

5 Установить зависимость между изменением вида напряженного состояния образцов и ориентацией скелетных линий изображений реакции ДДИТ Разработать методику идентификации условий циклического нагружения на основе анализа скелетизированных изображений

6 Разработать новые математические модели восстановления одно- и двумерных тарировочных зависимостей с применением строгих математических критериев оптимальности

7 Разработать специализированное программное обеспечение, объединяющее все предложенные модели и обеспечивающее алгоритмизацию и автоматизацию решения задач структурного и ситуационного анализа показаний ДДИТ

Объект и предмет исследования. Исследования выполнены в области цифровой обработки и анализа изображений Предметом исследования являются объекты на полутоновых изображениях реакции ДДИТ в виде "темных пятен", которые являются следствием необратимых смещений атомов кристаллической решетки в результате накопления пластических деформаций в материале ДДИТ при циклическом нагружении

Методы исследований. В качестве теоретических методов исследования использовались методы цифровой обработки изображений, теории вероятности, математической статистики, математического анализа и математического моделирования Для программной реализации

разработанных математических моделей и алгоритмов использовались методы создания программных систем, методы оптимизации программных комплексов и программирование на языках высокого уровня (Visual С++)

Моделирование и основной вычислительный эксперимент проводились с использованием математических пакетов (MathCAD, MatLab для предварительных исследований и проверки гипотез) и пакета программ, реализованных в среде Visual С++ (основной вычислительный эксперимент)

Научная новизна диссертационной работы заключается в следующем

1 Модель комбинированной фильтрации представляет собой новый перспективный подход к дифференцированному сглаживанию фрагментов изображений реакции ДДИТ Предложенные способы модификации данной модели базируются на свойствах локальной адаптации комплекса применяемых преобразований к особенностям изображений Задачи построения моделей фильтрации, связанные с получением качественных формальных описаний изображений реакции ДДИТ как объектов анализа, решаются при помощи изучения структуры и содержания изображений

2 Модель сегментации, разработанная на основе маркировки связных объектов («темных пятен»), представляет собой иерархическую кластер-процедуру и обеспечивает автоматическое разделение малоконтрастных объектов, расположенных на сложном фоне изображений реакции ДДИТ Эта модель является основой для реализации последующих процедур структурного и морфологического анализа изображений

3 Новые модели скелетизации бинарных изображений реакции ДДИТ позволяют решить проблему сохранения геометрических и топологических характеристик изображений На основе компьютерного анализа структуры скелетизированного изображения реализована методика идентификации некоторых условий нагружения образцов деталей и конструкций машин в месте наклейки ДДИТ (изгиб, кручение, изгиб с кручением, растяжение)

4 Разработаны новые математические модели восстановления одно- и двумерных тарировочных зависимостей, которые позволяют реализовать

степень оптимального соотношения между точностью аппроксимации эмпирического материала и надежностью полученного результата в условиях ограниченного объема данных (и < 20)

5 Разработан и внедрен программный комплекс 1та§еАпа1узег, который, реализуя построенные модели и алгоритмы, объединяет в себе результаты проведенных теоретических и экспериментальных исследований и позволяет осуществлять научные исследования усталости на основе компьютерной обработки изображений реакции ДДИТ

На защиту выносятся следующие научные положения:

1 Разработанная модель комбинированной фильтрации изображений представляет собой способ эффективной формализации изображений реакции ДДИТ, который дает возможность работать с представлениями, отражающими семантику изображения, т е информацию, заключенную в его внутренней структуре и структуре внешних связей с изменениями, происходящими в материале исследуемой детали Использование данной модели позволяет получить количественные оценки измененного состояния поверхности ДДИТ и осуществить калибровку датчиков по величине относительной площади "темных пятен"

2 Новая модель маркировки бинарных изображений реакции ДДИТ решает проблему автоматического обнаружения и сегментации малоконтрастных объектов на сложном фоне изображений При помощи данной модели можно выделить различные объекты изображений, выявить пространственные связи между ними, оценить их количество и размеры

3 Разработанные модели линеаризации изображений (выделение контура, скелетизирование) не нарушают условия связности (непрерывности) объектов и сохраняют их топологию Компьютерный анализ топологической структуры скелетизированных изображений реакции ДДИТ определяет ориентацию "темных пятен" на поверхности датчика, идентифицирует условия циклического нагружения, что дает возможность получать

качественную информацию о характере усталостных повреждений в процессе испытаний объектов на выносливость

4 Математические модели восстановления одно- и двумерных тарировочных зависимостей позволяют решить целый комплекс нетривиальных проблем интерпретации экспериментальных данных, одна из важнейших задач которого состоит в определении функциональных связей между характеристиками процесса нагружения и фиксируемыми результатами наблюдений по заданному критерию оценки реакции ДДИТ Использование этих моделей позволяет с заданной надежностью (порядка 0,95) добиться оптимальной степени близости аппроксимирующей функции к экспериментальным значениям

Практическая ценность диссертации.

1 Разработанные математические модели и алгоритмы обработки изображений реакции ДДИТ обеспечивают однозначность идентификации выявленных параметров деформирования исследуемой детали или конструкции (число и амплитуда нагружения, поврежденность материала по заданному критерию оценки реакции датчика в виде относительной площади «темных пятен»), повышают информационность и достоверность методов диагностики усталости по показаниям ДДИТ Простота и вычислительная эффективность разработанных моделей обеспечивают обработку экспериментальных данных в режиме реального времени

2 Модули пакета программ, разработанного автором на языке Visual С++ и оптимизированного по критерию минимума показателя вычислительной сложности, легко интегрируются в другие цифровые системы анализа и обработки изображений для решения практических задач оценки остаточной работоспособности деталей и конструкций машин по заданному критерию усталостного повреждения на основе показаний ДДИТ

Реализация результатов исследования Программный комплекс ImageAnalyser, включающий в себя разработанные математические модели, используется для обработки показаний ДДИТ и оценки усталостного

повреждения элементов ходовой части легковых автомобилей марки ВАЗ-2110, ВАЗ-2111, ВАЗ-2112, проходящих регулярное техническое обслуживание и ремонт в «Автоцентре на Омской» (ЗАО «Апрель», г Курган), что подтверждается соответствующим актом о внедрении

Личный вклад соискателя. Все предлагаемые математические модели и алгоритмы были разработаны и программно реализованы лично автором Научный руководитель принимал участие в постановке задач, определении возможных путей решения и их предварительном анализе

Апробация работы. Основные результаты и положения работы обсуждались на следующих конференциях международная научно-техническая конференция «Проблемы развития ТЭК Западной Сибири на современном этапе», посвященная 40-летию Тюменского государственного нефтегазового университета (г Тюмень, ТюмГНГУ, 2003 г), 2-я международная научно-техническая конференция «Новые технологии в нефтяной отрасли и образовании» (г Тюмень, ТюмГНГУ, 2006 г )

Публикации. По теме диссертации опубликовано 12 печатных работ, из них 2 - в центральных журналах, 7 статей - в сборниках научных трудов, 3 доклада на международных конференциях

Структура и объем диссертации. Диссертационная работа состоит из введения, четырех глав, заключения, библиографии и приложений Материал работы изложен на 149 страницах основного текста, содержит 24 рисунка, 8 таблиц, 1 приложение Библиография включает 84 наименования

Автор выражает благодарность доценту кафедры «Технология и автоматизация сварочного производства» Курганского государственного университета, ктн Троценко ДА за консультации и помощь в процессе работы над диссертацией

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ Во введении обоснована актуальность выполненного исследования Дана краткая характеристика состояния проблемы, поставлены цель и задачи

исследования, сформулированы научная новизна и практическая значимость работы, аннотируются основные положения, выносимые на защиту

В первой главе диссертационной работы выполнен анализ состояния предметной области Проведен обзор существующих методов исследования распределения напряжений и деформаций на поверхности деталей и конструкций машин в условиях циклического нагружения (дифференциальный и интегральный) Установлено, что новые методики, основанные на применении ДДИТ, обладают потенциально высокими возможностями На основе анализа существующих методов оценки регистрации и оценки структурных изменений ДДИТ (по моменту появления первых «темных пятен», по первым зернам измененной структуры, по измерению микротвердости датчика, с помощью оригинальных оптических средств измерения отражательной способности поверхности датчика, по относительной площади «темных пятен») сделан вывод о том, что использование возможностей компьютерных технологий для обработки и идентификации цифровых изображений реакции ДДИТ является перспективным направлением исследований в этой области

Рассмотрены основные задачи и особенности обработки изображений реакции ДДИТ, среди которых можно выделить массовый характер решаемой задачи, наличие количественных критериев точности, сложную модель наблюдения объекта, ограничение на вычислительные ресурсы и на время ответа, необходимость объединения различных математических моделей и алгоритмов в единую систему математического обеспечения

Во второй главе «Разработка новых математических моделей и алгоритмов компьютерной обработки и анализа цифровых изображений реакции датчиков деформаций интегрального типа» рассмотрена проблематика и особенности задач обработки и идентификации показаний датчиков Показано, что изображения реакции ДДИТ можно отнести к идентифицируемым объектам с плохо определенной структурой (нечетким образам) Удаление шумов на таких изображениях без внесения искажений в

информативные детали и структуры представляет собой одну из главных проблем автоматизации процессов обработки и анализа изображений

Сравнительный анализ известных моделей фильтрации (фильтр Винера, медианный фильтр и его модификации, статистический фильтр, низкочастотный и высокочастотный пространственные фильтры, фильтрация с помощью преобразования Фурье, вейвлет-преобразование и др ) показал, что использование этих моделей обеспечивает формирование слабоконтрастных изображений реакции ДДИТ, приводит к размыванию границ объектов и к стиранию значимых элементов Кроме того, при тарировании ДДИТ по относительной площади «темных пятен» наблюдается неконтролируемый разброс экспериментальных данных в зависимости от выбранного режима фильтрации

Разработана модель маркировочного фильтра Использование в рамках этой модели дифференцированного подхода к сглаживанию различных областей изображений реакции ДДИТ позволяет достичь компромиссного решения известного противоречия между сглаживанием шума и сохранением контрастных структур различной ориентации

Модель сегментации объектов изображений реакции ДДИТ основана на элементах кластерного анализа с применением иерархических методов классификации Особенность данной модели заключается в сегментации объектов бинарного изображения по принципу их связности и маркировке выделенных максимальных связных компонент, что дает возможность реализовать методики автоматической классификации и распознавания объектов изображений реакции ДДИТ

Рассмотрены способы выделения контуров изображений (методы высокочастотной фильтрации, пространственного дифференцирования, функциональной аппроксимации, прослеживания границ) Показано, что на основе этих моделей контуры объектов изображений реакции ДДИТ получаются незамкнутыми либо утолщенными, поэтому необходимость сохранения связности автоматически порождает проблемы, относящиеся к

ю

топологии и связанные с порождением ложных линий Разработанная модель выделения контура осуществляется на основе эффективного подбора весовых множителей специального оператора К Выходной эффект данного преобразования определяется линейной комбинацией входных данных

11 ):

Здесь ^ = \/1;}, < = 1,т, у = 1,п, - дискретная матрица изображения, е.?, 5- апертура фильтра, /<(!,,у,) - весовые коэффициенты, совокупность которых представляет собой двумерную импульсную характеристику оператора, О - сумма весовых коэффициентов

Модель скелетизации изображения реакции ДДИТ основана на метрическом определении скелета точка Р является скелетной точкой множества Л, если Л(Р,В) = с1{Р,А/,) = ¿(ДМ.,), где М,,М2 - точки, лежащие на границе В данного множества Я, причем с1(Р, В) = т^сЦР, У), У е В}

Компьютерная реализация разработанной модели заключается в постепенном удалении одного слоя граничных пикселов каждой промаркированной связной компоненты изображения

Ф, = /®Ф. , Ф, = /®Ф, , Ф„ = /®Ф„, Ф„ = /®Ф„ Здесь Ф - операция сложения по модулю 2, /- 1-я строка изображения,

фв =а, л/}, =Зл7л/;(2A5v4лlvlл5), фг =а, л/?, = 5 лlл/,(4л7v6лЗvЗл7),

Ф„ =«, * А, =7лЗл /,(6л1у8л5у5л1), ф„ = а„ л Д, =Та5л/,1(8л2у2л6у6лЗ)

- булевы функции, определяющие крайние верхние, левые, нижние и правые

пикселы каждого промаркированного объекта

изображения соответственно Для удобства

математического описания процедуры элементы

8-связной окрестности рабочего пиксела

обозначим в соответствии с рис 1

Разработаны модели восстановления одно-

и двумерных тарировочных зависимостей, которые определяют

функциональную связь между характеристиками процесса нагружения

п

г4 3 2")

5 /„ 1

,6 7 8, Рис 1 Нумерация 8-связных пикселов

(амплитуда, число циклов деформирования) и фиксируемыми результатами наблюдений (относительная площадь «темных пятен») Модели основаны на теории минимизации функционала эмпирического риска

/,(«) = -¿(у, -Р(х„а))\

который соответствует минимуму среднеквадратической ошибки построения приближающей функции Для зависимостей 8 = ¿>(М) функционал /-,(<?) имеет вид

п -ЕГ а;

где (N,,6,), 1 = 1, ,п- независимые пары значений аппроксимируемой функции 8 (относительная площадь "темных пятен") и ее аргумента N (число циклов нагружения образца), <х(2 - дисперсии замеров 3,

Построение одномерных тарировочных зависимостей реализовано в

классе алгебраических полиномов (полиномиальная регрессия)

к

и в классе полиномов Чебышева

<5(ЛО = 2>,е,(Л>),Где е/(Ю = со8(у агссо8(Л0)

Для восстановления многомерной тарировочной зависимости используется модифицированный метод всех регрессий В его основе лежит использование многофакторной модели следующего вида

у = ао /,(*,) /г(х2) /„(*„) Искомая функция представляет собой произведение одномерных базовых функций /,(*,)> параметры которых рассчитываются на основе разработанной модели построения полиномиальной регрессии

В третьей главе «Компьютерная реализация и результаты тестирования математических моделей и алгоритмов обработки изображений реакции

датчиков деформаций интегрального типа» рассмотрены вопросы программной реализации разработанных моделей

Приведено описание возможностей и интерфейса программного комплекса ImageAnalyser, разработанного для проведения вычислительного эксперимента, для создания которого была использована интегрированная среда разработки программ - Microsoft Visual С++ Программный комплекс содержит меню пользователя, архивы данных, графическое представление получаемых результатов. Используется многодокументный интерфейс, что позволяет вести обработку большого количества изображений в рабочем поле программы Для ряда процедур осуществляется вызов специальных диалоговых окон — мастеров Работа мастеров демонстрируется при реализации конкретных процедур обработки изображений

Для улучшения результатов фильтрации при помощи маркировочного фильтра осуществлена классификация изображений реакции ДДИТ по особенностям представления информативных объектов (размер объектов и неравномерность основного фона) и рассмотрены модификации маркировочного фильтра, разработанные для каждого класса изображений Модели фильтрации обладают важным свойством локальной адаптации к характеристикам обрабатываемого изображения, что обеспечивает достоверность последующих процедур их идентификации

Установлено, что эффективность и точность разработанной модели маркировки находится в непосредственной зависимости от промежуточного значения счетчика неполных компонент, неконтролируемое увеличение которого влечет за собой большую погрешность, поэтому при обработке пикселов изображения нельзя забывать о коррекции результата по яркостному диапазону (0 255)

Компьютерная реализация моделей линеаризации изображений реакции ДДИТ подтвердила их высокую эффективность Для реализации модели скелетизации требуется не более /,ъм / 2 циклов, где /,„„ — максимальная толщина линии скелетизированного изображения

Проведен анализ изображений реакции ДДИТ при различных видах циклического нагружения (изгиб, кручение, изгиб с кручением, растяжение) Разработана модель слежения скелетной линии, методика для идентификации данных условий нагружения

Дня компьютерного тестирования математических моделей восстановления одно- и двумерных тарировочных зависимостей использованы задачи построения полиномиального приближения к регрессии, значения которой заданы со случайной помехой Качество аппроксимации на основе экспериментальной выборки ограниченного объема п с вероятностью 1-77 = 095 определялось по равномерной оценке величины среднего риска для заданного класса функций емкости И

В четвертой главе разработана общая стратегия получения прогнозируемых оценок накопления усталостных повреждений в материале конструкций магистральных трубопроводов при помощи программного комплекса ImageAnalyser Рассмотренная расчетно-экспериментальная методика обработки изображений реакции ДДИТ была использована для исследования прочностных свойств образцов конструкции магистрального трубопровода (диаметр 720 мм, толщина стенки 9 мм, материал трубы -сталь 14 ХГС) Испытания проведены на машине усталостных испытаний МУИ-6000 при циклическом растяжении до наработки соответствующего числа циклов (N-86000) или до поломки образцов Используемые в эксперименте ДДИТ были изготовлены из алюминиевой фольги в виде прямоугольной полоски размером 4x60 мм

Выполнена калибровка ДДИТ, деформируемых при изменении длительности нагружения образцов, по относительной площади «темных

пятен» (рис. 2). Эти данные фиксируют не только момент появления реакции в виде первых «темных пятен», но и динамику их распространения по поверхности ДДИТ.

а) ЛГ = 30-10', <5 = 1.4045 б) N = 50-10', <5 = 3.6250 в) N = 81 ■ 10', 8 = 12.9071

Рис. 2. Обработка изображений реакции ДДИТ

Структурный анализ скелетных линий изображений реакции ДДИГ, полученных в условиях циклического растяжения образца, подтверждает гипотезу о том, что эти линии располагаются вдоль оси образца. Для построения основных расчетных зависимостей метода прогнозирования восстановлена тарировочная зависимость (рис. 3), которая обеспечивает необходимое качество аппроксимации с вероятностью 0,95 при а = 265.06 МПа и 0 <ы <81000:

8 = 0.3701 + 0.0662- х + 0.9116 л2 + 6,1916 + 5.0771- х1, где х = 1%Ы.

Коэффициент корреляции экспериментальных данных равен 0.9907. Построены зависимости при усталостных испытаниях идентичных образцов трубы для значений напряжения ст = 235.2 МПа и сг = 280 МПа.

Проведен дисперсионный анализ результатов тарирования ДЦИТ. Определена величина остаточного рассеивания: дт„ = -д'(х,))г = 4.3654 и

остаточная дисперсия линии регрессии:

4,3654

» 0,2425, которые

подтвердили, что построенные математические модели хорошо согласуются с экспериментальными данными.

. .IVay.ii.iHpyn.,,,,,». «•«•.mi»r|,<„»,»-

■и ■■ - ~ " '•' ' -Л»,■>•-." ^ • vf ■»•''■•л'i

х = 3.3010 у- 0.0040

0 у

« = ЗЛ771 у = 0.0056 х = 4.0000 у = 0.058В

. * 4 2.104 у 0 00«

. --«Ж?

: •Ж1'.'"".''." ii'" • • •'«. '

/ с^пень построенного полинома в

' si

. х Л-'-Hllh у = 0.186?

X 4.11,17 у ЮПИ!

х 4.4149 у 0.5337

■• X = 4.4771 у = 0.8541 *

х - 4.5314 у = 1.2Э07 :-

X 4.568? у 1.5454

■ т. ■ life,

х = 4.6627 у = 2.6138 ix = 4.6989 у- 3.1371: х- 4.7243; у = 3.5468

.X 4.7922 у 4.11311

: " ' х = 4.8325 у = 5.7362 .

. х = 4.8579, у = 6.3664

х = 4.8921 у = 7.2926

X 4.9031 у =7.6103

х = 4.9084 у = 7.7670

КГ-дт'ТГГТп - i м*<

• •Л «С

• Я*;;

t- в виде обычного полиноиа -: ...... *•■■ ; " '

Рис. 3. Коэффициенты тарировочной зависимости

Для каждого значения доверительной вероятности р = 0,95 и р = 0,99 определены границы доверительных интервалов:

д* - = £(*,) ±0.2313, 6* = = <?(*,)+0.2313

sin ~Jn

Здесь tnp - статистика, имеющая /-распределение Стьюдента с к - п-2 степенями свободы; S , 5' - нижняя и верхняя граница доверительного интервала соответственно; п- объем экспериментальной выборки.

Определена доля объясненной вариации

= 4 -100% = 3816193 я 98 48У 51, 393.6195

где = £> = ¿(<5,- сумма квадратов отклонений, обусловленная общим

разбросом параметра fi относительно его среднего значения;

,)-з)г - сумма квадратов отклонений, обусловленная

направленным трендом

Установлена значимость влияния двух одновременно действующих факторов (.V,, <т1) на величину усталостного деформирования образцов 8

Восстановлена двумерная тарировочная зависимость 6 = 8(Ы,а) при О <N<80000, 235 2 < о-<282 0 (МПа) Коэффициент множественной корреляции составил 0 9973

Комплекс разработанных математических моделей адекватно интерпретирует экспериментальные данные, полученные при помощи ДДИТ С его помощью реализована методика прогнозирования величины напряжений, накопленного усталостного повреждения и числа циклов нагружения образцов конструкции магистрального трубопровода (при наработке до 86000 тыс циклов) на заданном диапазоне тарирования ДДИТ

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИОННОЙ РАБОТЫ

1 Осуществлена классификация изображений реакции ДДИТ и разработаны дифференцированные модели фильтрации, адаптированные к особенностям каждого класса изображений Эти модели обеспечивают генерализацию изображений реакции ДДИТ и калибровку ДДИТ по величине накопленного усталостного повреждения в локальных зонах металлоконструкций в виде относительной площади «темных пятен»

2 Разработана модель маркировки «темных пятен» бинарных изображений реакции ДДИТ Промаркированное изображение является основой для подсчета ряда морфологических параметров (площадь, периметр, длина, ширина, координаты центра масс, моменты инерции) для формирования признакового пространства исследуемых объектов

3 Модели линеаризации изображений реакции ДДИТ выделяют непрерывные контуры и скелеты информативных объектов, не нарушая их

топологии Разработанная модель анализа топологической структуры скелетизированных изображений реакции ДДИТ определяет ориентацию «темных пятен» на поверхности ДДИТ и идентифицирует некоторые условия циклического деформирования образцов (изгиб, кручение, изгиб с кручением, растяжение)

4 Разработаны модели восстановления оптимальных тарировочных зависимостей в условиях малого объема экспериментальных данных (« < 20) с заданной надежностью 0,95 Аналитическое выражение зависимостей, связывающих микроструктуру материала с его физическими и механическими свойствами, представляет собой простой метод прогнозирования работоспособности деталей и конструкций машин Этот метод базируется на традиционном подходе - сопоставлении результатов тарирования ДДИТ и данных усталостных испытаний деталей с ДДИТ при одной и той же реакции датчика Число циклов нагружения образцов до реакции ДДИТ составляет до 10 тыс циклов, а прогнозируемое число циклов до поломки - сотни тысяч - десятки миллионов циклов, поэтому глубина прогноза превышает на два-три порядка интервал наблюдения

5 Разработан и внедрен программный комплекс 1та§еАпа1у8ег, который, реализуя разработанные математические модели, объединяет в себе результаты проведенных теоретических и экспериментальных исследований и осуществляет научные исследования усталости на основе обработки цифровых изображений реакции ДДИТ Реализация построенных математических моделей представляет собой новую методику расшифровки реакции ДДИТ, обеспечивает алгоритмизацию и автоматизацию решения задач структурного и ситуационного анализа показаний ДДИТ

4 Эффективность созданных моделей подтверждена серией стендовых экспериментов, проведенных для проверки адекватности разработанных математических моделей на примере образцов нефтегазового трубопроводного оборудования Результаты испытаний подтверждают, что изложенная методика обработки показаний ДДИТ работоспособна и дает

возможность решать задачи оценки состояния различных металлоконструкций в условиях фактического нагружения с минимальным количеством экспериментов ( и 5 20)

5 Разработанные математические модели и алгоритмы могут служить основой для разработки новых процедур использования ДДИТ в практике экспериментальных работ, учитывающих особенности функционирования исследуемых металлоконструкций и возможности получения информации, требуемой для расшифровки и анализа показаний датчиков Идеология построения математических моделей легко может быть распространена на целый ряд подобных вопросов В частности, разработанное программное обеспечение может быть использовано при решении вопросов обеспечения надежности сложных диагностических систем и комплексов

СПИСОК РАБОТ ПО ТЕМЕ ДИССЕРТАЦИИ

1 Сызранцев В Н, Змызгова Т Р, Змызгов С Е Одномерная аппроксимация результатов тарирования датчиков деформаций интегрального типа // Проблемы развития ТЭК Западной Сибири на современном этапе Труды международной научно-технической конференции, посвященной 40-летию Тюменского государственного нефтегазового университета Т 2 - Тюмень ТюмГНГУ, 2003 —С 105-108

2 Сызранцев В Н, Змызгова Т Р, Змызгов С Е Алгоритмы аппроксимации результатов тарирования датчиков деформаций интегрального типа // Проблемы развития ТЭК Западной Сибири на современном этапе Труды Международной научно-технической конференции, посвященной 40-летию Тюменского государственного нефтегазового университета Т 2 - Тюмень ТюмГНГУ, 2003 -С 126-129

3 Сызранцев В Н , Змызгова Т Р Алгоритмы бинарного квантования цифровых изображений реакции датчиков деформаций интегрального типа // Современные методы изучения пластов и скважин при решении задач разработки нефтяных и газовых месторождений Сборник трудов кафедры

«Разработка и эксплуатация газовых и газоконденсатных месторождений» Вып 1 -Тюмень ТюмГНГУ, 2004 -С 214-220

4 Сызранцев В H , Змызгова Т Р Структурный анализ цифровых изображений реакции датчиков деформаций интегрального типа // Современные методы изучения пластов и скважин при решении задач разработки нефтяных и газовых месторождений Сборник трудов кафедры «Разработка и эксплуатация газовых и газоконденсатных месторождений» Вып 1 -Тюмень ТюмГНГУ, 2004 -С 221-227

5 Сызранцев В H , Троценко Д А , Змызгова Т Р Предварительная обработка цифровых изображений реакции датчиков деформаций интегрального типа для оценки поврежденности сварных соединений // Современные методы изучения пластов и скважин при решении задач разработки нефтяных и газовых месторождений. Сборник трудов кафедры «Разработка и эксплуатация газовых и газоконденсатных месторождений» Вып 1 -Тюмень ТюмГНГУ, 2004 - С 228-233

6 Змызгова Т Р Алгоритм сегментации цифровых бинарных изображений // Аспирант и соискатель Вып 2 — M Спутник-!-, 2005 - С 159-162

7 Змызгова Т Р Методы фильтрации цифровых изображений реакции ДДИТ с целью количественной оценки накопления усталостных повреждений в сварных соединениях // Сборник научных трудов аспирантов и соискателей Курганского государственного университета Вып 7 — Курган Изд-во Курганского гос ун-та, 2005 -С 18-19

8 Змызгова Т Р Проблематика и выбор средств экспериментальной оценки накопленных усталостных повреждений в материале магистральных трубопроводов // Сборник научных трудов аспирантов и соискателей Курганского государственного университета Вып 8 - Курган Изд-во Курганского гос ун-та, 2006 - С 21-22

9 Змызгова Т Р Исследование процесса накопления пластических деформаций в датчиках деформаций интегрального типа // Естественные и технические науки Вып 2-М Спутник+, 2006 - С 110-112

10 Сызранцев В Н, Змызгова Т Р, Невелев Я П Компьютерная обработка показаний датчиков деформаций интегрального типа // Новые технологии в нефтяной отрасли и образовании Сборник материалов 2-й международной научно-технической конференции - Тюмень ТюмГНГУ, 2006 - С 73-76

11 Сызранцев В Н , Голофаст С Л , Змызгова Т Р Обработка данных тарирования датчиков деформаций интегрального типа по критерию изменения их микротвердости //Известия вузов Машиностроение - 2006 -№12 - С 12-15

12 Сызранцев В Н , Голофаст С Л, Обакшин П А , Змызгова Т Р // Восстановление многомерных функциональных зависимостей по выборкам ограниченного объема // Вестник машиностроения - 2007 - № 8 - С 30-32

Змызгова Татьяна Рудольфовна

МАТЕМАТИЧЕСКИЕ МОДЕЛИ И АЛГОРИТМЫ КОМПЬЮТЕРНОЙ ОБРАБОТКИ ЦИФРОВЫХ ИЗОБРАЖЕНИЙ РЕАКЦИИ ДАТЧИКОВ ДЕФОРМАЦИИ ИНТЕГРАЛЬНОГО ТИПА

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Подписано к печати 5 10 2007 Формат 60x84 1/16 Бумага тип № 1

Печать трафаретная Уел печ л 1,5 Уч-издл 1,5

Заказ № 933 Тираж 100 экз Бесплатно

РИЦ Курганского государственного университета, 640669, г Курган, ул Гоголя, 25 Курганский государственный университет

Оглавление автор диссертации — кандидата технических наук Змызгова, Татьяна Рудольфовна

ВВЕДЕНИЕ.

1. ОБЗОР СОСТОЯНИЯ ПРОБЛЕМЫ И ПОСТАНОВКА ЗАДАЧ ИССЛЕДОВАНИЙ.

1.1. Краткая характеристика датчиков деформаций интегрального типа.

1.2. Методы регистрации и оценки структурных изменений датчиков деформаций интегрального типа.

1.3. Основные средства и задачи построения математических моделей и алгоритмов обработки изображений.

1.4. Постановка задач исследования.

2. РАЗРАБОТКА НОВЫХ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ И АЛГОРИТМОВ КОМПЬЮТЕРНОЙ ОБРАБОТКИ И АНАЛИЗА ЦИФРОВЫХ ИЗОБРАЖЕНИЙ РЕАКЦИИ ДАТЧИКОВ ДЕФОРМАЦИЙ ИНТЕГРАЛЬНОГО ТИПА.

2.1. Проблематика создания математических моделей для компьютерной обработки и идентификации изображений реакции датчиков.

2.2. Модели предварительной обработки изображений реакции ДДИТ.

2.2.1. Сравнительный анализ известных моделей фильтрации.

2.2.2.Модель маркировочной фильтрации.

2.3. Модель автоматической сегментации объектов бинарных изображений реакции ДДИТ на основе маркировки связных компонент.

2.4. Модели линеаризации изображений реакции ДДИТ.

2.4.1. Методы выделения контуров. Линейная модель построения непрерывного контура объектов изображений.

2.4.2. Модель скелетизации бинарных изображений.

2.4.3. Способы улучшения скелетизированных изображений.

2.5. Модели восстановления тарировочных зависимостей.

2.5.1. Задача восстановления функциональных зависимостей по выборкам ограниченного объема.

2.5.2. Метод структурной минимизации эмпирического риска.

2.5.3. Восстановление одномерных тарировочных зависимостей для ДДИТ в классе алгебраических полиномов.

2.5.4. Модифицированная многофакторная модель построения полиномиального приближения многомерной тарировочной зависимости.

2.6. Выводы.'.

3. КОМПЬЮТЕРНАЯ РЕАЛИЗАЦИЯ И РЕЗУЛЬТАТЫ ТЕСТИРОВАНИЯ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ И АЛГОРИТМОВ ОБРАБОТКИ ИЗОБРАЖЕНИЙ РЕКЦИИ ДАТЧИКОВ ДЕФОРМАЦИЙ ИНТЕГРАЛЬНОГО ТИПА.

3.1. Структурная схема системы автоматизированной обработки цифровых изображений реакции ДДИТ.

3.2. Прикладное программное обеспечение и исследование эффективности разработанных математических моделей и алгоритмов.

3.2.1. Улучшение визуального качества изображений реакции ДДИТ при помощи модели маркировочной фильтрации. Классификация изображений реакции ДДИТ. Модифицированная модель маркировочной фильтрации.

3.2.2. Сегментация объектов бинарных изображений реакции ДДИТ по принципу их связности.

3.2.3. Тестирование моделей выделения контуров и скелетизации бинарных изображений реакции ДДИТ.

3.2.4. Модель слежения скелетной линии. Идентификация вида напряженного состояния образцов нз;'основе структурного анализа скелетизированных изображений реакции ДДИТ.

3.2.5. Модель построения полиномиального приближения одномерной тарировочной зависимости и ее компьютерное тестирование.

3.2.6. Компьютерная реализация модели восстановления двумерной тарировочной зависимости.

3.3. Выводы по разделу.

4. ОБЩАЯ СТРАТЕГИЯ ПОЛУЧЕНИЯ ПРОГНОЗИРУЕМЫХ ОЦЕНОК НАКОПЛЕНИЯ УСТАЛОСТНЫХ ПОВРЕЖДЕНИЙ В МАТЕРИАЛЕ КОНСТРУКЦИЙ МАГИСТРАЛЬНЫХ ТРУБОПРОВОДОВ НА ОСНОВЕ КОМПЬЮТЕРНОЙ ОБРАБОТКИ ЦИФРОВЫХ ИЗОБРАЖЕНИЙ РЕАКЦИИ ДЦИТ.Г.

4.1. Проблематика и выбор средств оценки усталостных повреждений в конструкциях нефтегазового трубопроводного оборудования.

4.2. Расчетно-экспериментальная методика регистрирования и оценки накопленного усталостного повреждения с помощью ДЦИТ.

4.2.1. Краткое изложение методологии решения задачи диагностики остаточного ресурса металлоконструкций по показаниям ДЦИТ. Построение кривой выносливости.

4.2.2. Постановка эксперимента для исследования прочностных свойств конструкций магистральных трубопроводов при циклическом растяжении по показаниям ДЦИТ.

4.2.3. Результаты тарировочных испытаний. Взаимосвязь реакции ДЦИТ с процессом накопления усталостных повреждений в материале образцов при испытаниях на выносливость.

4.2.4. Оценка величины накопленного усталостного повреждения в образце трубопровода на основе цифровых изображений реакции ДЦИТ. Калибровка ДЦИТ по относительной площади "темных пятен".

4.2.5. Определение режима нагружения тарировочного образца на основе компьютерного анализа топологической структуры скелетизированного изображения реакции ДЦИТ.

4.3. Математическая обработка результатов тарировочных испытаний. Основные расчетные зависимости метода прогнозирования.

4.3.1. Определение одномерной функциональной зависимости между одним из параметров циклического нагружения и количественным выражением реакции ДДИТ. Прогнозирование числа циклов нагружения образцов на основе корреляционной связи с реакцией ДЦИТ.

4.3.2. Регрессионный анализ результатов тарирования ДЦИТ. Проверка значимости одномерной тарировочной зависимости. Оценка метрологической точности измерений.

4.3.3. Дисперсионный анализ оценки влияния двух одновременно действующих факторов на величину усталостного деформирования тарировочного образца.

4.3.4. Построение полиномиального приближения двумерной тарировочной зависимости. Определение величины напряжений и накопленного усталостного повреждения по показаниям ДЦИТ.

4.4. Выводы по разделу.

Введение 2007 год, диссертация по информатике, вычислительной технике и управлению, Змызгова, Татьяна Рудольфовна

Анализ ситуации, сложившейся на современном этапе развития отечественного машиностроения, показывает, что создание конкурентоспособных образцов продукции неразрывно связано с необходимостью технической диагностики изделий, конечная цель которой заключается в получении данных для оценки их нагруженности и работоспособности в заданных условиях эксплуатации. Данная проблема особенно актуальна для строительных, авиационных конструкций, для нефтепроводов, газопроводов, где развитие усталостных трещин может повлечь за собой крупные аварии и катастрофы. Чтобы их исключить, необходимо иметь надежные средства и методы регистрации фактического усталостного повреждения.

Многочисленные теоретические и экспериментальные методы, разработанные к настоящему времени, несмотря на свой достаточно высокий уровень, не позволяют с высокой степенью достоверности предсказывать индивидуальный ресурс деталей машин и металлоконструкций [7-10, 20, 27-28, 33, 37, 42, 58, 59, 63-69, 72]. Основными причинами сложившегося положения являются неопределенность характеристик важнейших величин, используемых в расчетах на длительную прочность: предела выносливости материала деталей и предела выносливости деталей, напряжений, возникающих в местах вероятного разрушения деталей, а также ошибки, вносимые использованием гипотез суммирования усталостных повреждений. Поэтому разработка новых средств и методов, позволяющих экспериментально оценивать нагруженность и накопление усталостных повреждений в деталях различных конструкций, является актуальным направлением в науке и технике.

Существенный вклад в развитие данного направления внесли ученые В.Н. Сызранцев [64-69], В.В. Болотин [9], С.А. Тимашев [71], В.И. Бойко [78], Э.Г. Керимов [27,28], В.А. Копнов [33], Д.А. Троценко [73] и другие.

Одним из перспективных направлений решения указанной проблемы является использование методик, основанных на применении датчиков деформации интегрального типа (ДЦИТ). ДДИТ представляют собой новые средства экспериментального исследования. Благодаря универсальности и большим потенциальным возможностям датчиков, появилась уникальная возможность решения широкого спектра задач по исследованию характера изменения напряжений и деформаций на поверхностях деталей, восстановлению закона их распределения и определению эквивалентных напряжений, оценке мест и сроков вероятного разрушения деталей при стендовых и эксплуатационных испытаниях машин.

Интенсивное развитие и интеграция информационных и компьютерных технологий в научные исследования, актуальность развития существующих способов определения усталостных повреждений и прогнозирования ресурса деталей и конструкций машин на основе показаний ДДИТ привели к необходимости создания новых автоматизированных методов обработки экспериментальных данных, полученных при помощи датчиков. Проведенные исследования подтверждают, что это новое направление позволит не только существенно расширить область применения ДДИТ, но и существенно повысить информационность и достоверность методов диагностики усталости деталей машин в заданных условиях эксплуатации, обеспечить однозначность идентификации выявленных параметров.

Решение задач компьютерной обработки показаний реакции ДЦИТ во многом зависит от интенсификации научных исследований и практического применения их результатов. В свою очередь интенсификация научных исследований требует широкого применения и развития методов хранения, обработки и анализа экспериментальной информации на базе современных вычислительных средств. Среди комплекса проблем, которые должны решаться в этом направлении, особенное значение приобретает создание теоретических и эмпирических предпосылок для построения автоматизированных систем обработки экспериментальной информации, представленной изображениями реакции ДДИТ. Несмотря на огромный интерес к цифровой обработке экспериментальных данных, в литературе вопросам анализа и идентификации изображений уделяется недостаточное внимание. Анализ существующих методов обработки изображений свидетельствует о том, что все они применительно непосредственно к изображениям реакции ДДИТ, помимо ограниченных возможностей, неэффективны в реализации и поэтому не могут непосредственно использоваться и требуют либо конструктивного изменения, либо создания принципиально новых подходов. Поэтому задача компьютерной обработки и распознавания изображений реакции ДДИТ, которой посвящена данная работа, представляется актуальной и имеет научную новизну.

Для решения поставленной задачи необходимо разработка новых математических моделей и алгоритмов, которые предназначены для реализации в автоматических системах анализа изображений и позволят осуществлять научные исследования усталости на основе компьютерной обработки изображений реакции ДДИТ. При этом для того, чтобы предотвратить отрыв между теоретически разработанными методами и их практической реализацией, необходима тщательная разработка соответствующего математического и программного обеспечения. Принимая во внимание эти факторы, в диссертации предлагается развитие известных методов и разработка новых математических моделей, предназначенных для компьютерной обработки и анализа изображений реакции ДДИТ. Практическая реализация этих моделей обеспечивает получение необходимых оценок уровня и характера распределения деформаций на исследуемых участках деталей машин в процессе ходовых йЬпытаний, что в свою очередь дает возможность целенаправленно осуществлять конструкторско-технологическую доработку.

Для достижения данной цели в диссертации автором были поставлены следующие задачи:

1. Рассмотреть различные модели шумовых воздействий и провести сравнительный анализ известных моделей фильтрации для предварительной обработки изображений реакции ДДИТ.

И'

2. Разработать новые модели обработки показаний ДДИТ, адаптированные к особенностям исследуемых изображений, которые обеспечат построение эффективного формального описания и позволят осуществить калибровку ДДИТ по относительной площади «темных пятен».

3. Разработать новую модель автоматической сегментации информативных объектов изображений реакции ДДИТ, исследовать ее эффективность и точность.

4. Разработать новые способы линеаризации бинарных изображений реакции ДДИТ (скелетизация, выделение контура).

5. Установить зависимость между изменением вида напряженного состояния образцов и ориентацией скелетных линий изображений реакции ДДИТ. Разработать методику идентификации условий циклического нагружения на основе анализа скелетизированных изображений.

6. Разработать новые математические модели восстановления одно- и двумерных тарировочных зависимостей с применением строгих математических критериев оптимальности.

7. Разработать специализированное и ^программное обеспечение, объединяющее все предложенные модели и обеспечивающее алгоритмизацию и автоматизацию решения задач структурного и ситуационного анализа показаний ДД ИТ.

Решение поставленных задач ориентировано на разработку соответствующих математических моделей и алгоритмов и подразумевает их адекватность исследуемым свойствам физической модели, под которой имеется в виду реакция ДДИТ в виде измененного состояния его поверхности, зафиксированная на соответствующем цифровом изображении. Эффективное использование перспективных новинок в технологии программирования позволяет обеспечить преемственность и правильность получаемых результатов и осуществить рациональное проектирование прикладного математического обеспечения.

Объект и предмет исследования. Исследования выполнены в области предварительной обработки и анализа цифровых изображений и экспериментальных данных, полученных при помощи ДДИТ. Предметом исследования являются объекты на полутоновых изображениях реакции датчиков деформации интегрального типа в виде "темных пятен", которые являются следствием необратимых смещений атомов кристаллической решетки в результате накопления пластических деформаций в материале ДДИТ при циклическом нагружении.

Методы исследований. В качестве теоретических методов исследования использовались методы цифровой обработки сигналов и изображений, теории вероятностей, математической статистики, теории случайных функций, математического анализа и математического моделирования. Для программной реализации разработанных математических моделей и алгоритмов использовались методы создания программных систем, методы оптимизации программных комплексов и программирование на языках высокого уровня (Visual С++).

Моделирование и основной вычислительный эксперимент проводились с использованием математических пакетов (Matchcad, MatLab для предварительных исследований и проверки гипотез) и пакета программ, реализованных в среде Visual С++ (основной вычислительный эксперимент).

Научную новизну диссертационной работы определяют впервые полученные и перечисленные ниже результаты:

1. Модель комбинированной фильтрации представляет собой новый перспективный подход к дифференцированному сглаживанию фрагментов изображений реакции ДДИТ. Предложенные способы модификации данной модели базируются на свойствах локальной адаптации комплекса применяемых преобразований к особенностям изображений. Задачи построения моделей фильтрации, связанные с получением качественных формальных описаний изображений реакции ДДИТ, как объектов анализа, решаются при помощи изучения структуры и содержания изображений.

2. Модель сегментации, разработанная на основе маркировки связных объектов («темных пятен»), представляет собой иерархическую кластер-процедуру и обеспечивает автоматическое разделение малоконтрастных объектов, расположенных на сложном фоне изображений реакции ДДИТ. Эта модель является основой для реализации последующих процедур структурного и морфологического анализа изображений.

3. Новые модели скелетизации бинарных изображений реакции ДДИТ, позволяют решить проблему сохранения геометрических и топологических характеристик изображений. На основе компьютерного анализа структуры скелетизированного изображения реализована методика идентификации некоторых условий нагружения образцов деталей и конструкций машин в месте наклейки ДДИТ (изгиб, кручение, изгиб с кручением, растяжение).

4. Разработаны новые математические модели восстановления одно- и двумерных тарировочных зависимостей, которые позволяют реализовать степень оптимального соотношения между точностью аппроксимации эмпирического материала и надежностью полученного результата, равной 0.95, в условиях ограниченного объема данных (п < 20).

5. Разработан и внедрен программный комплекс Image Analyser, который, реализуя построенные модели и алгоритмы, объединяет в себе результаты проведенных теоретических и экспериментальных исследований и позволяет осуществлять научные исследования усталости на основе компьютерной обработки изображений реакции ДДИТ.

Научные положения диссертации, выносимые на защиту:

1. Разработанная модель комбинированной фильтрации изображений представляет собой способ эффективной формализации изображений реакции ДДИТ, который дает возможность работать с представлениями, отражающими семантику изображения, т.е. информацию, заключенную в его внутренней структуре и структуре внешних связей с изменениями, происходящими в материале исследуемой детали. Использование данной модели позволяет получить количественные оценки измененного состояния поверхности ДДИТ и осуществить калибровку датчиков по величине относительной площади "темных пятен".

2. Новая модель маркировки бинарных изображений реакции ДДИТ позволяет решить проблему автоматического, обнаружения и сегментации малоконтрастных объектов на сложном фоне изображений. При помощи данной модели можно выделить различные объекты изображений, выявить пространственные связи между ними, оценить их количество и размеры.

3. Разработанные модели линеаризации изображений (выделение контура, скелетизирования) не нарушают условия связности (непрерывности) объектов и сохраняют их топологию. Компьютерный анализ топологической структуры скелетизированных изображений реакции ДДИТ позволяет определить ориентацию "темных пятен" на поверхности датчика, идентифицировать условия нагружения, что дает возможность получать качественную информацию о характере усталостных повреждений в процессе испытаний объектов на выносливость.

4. Математические модели восстановления одно- и двумерных тарировочных зависимостей позволяют решить целый комплекс нетривиальных проблем интерпретации экспериментальных данных, одна из важнейших задач которого состоит в определении функциональных связей между характеристиками процесса нагружения и фиксируемыми результатами наблюдений по заданному критерию оценки реакции ДДИТ. Использование этих моделей позволяет с заданной надежностью (порядка 0,95) добиться оптимальной степени близости аппроксимирующей функции к экспериментальным значениям.

Личный вклад соискателя. Все предлагаемые математические модели и алгоритмы были разработаны и программно реализованы лично автором.

Научный руководитель принимал участие в постановке задач, определении возможных путей решения и их предварительном анализе.

Достоверность результатов диссертационной работы обеспечивается строгим характером использованных методов, адекватностью и корректностью примененного математического аппарата и подтверждается результатами практической реализации и внедрения разработанных математических моделей.

Практическая значимость работы.

1. Разработанные математические модели и алгоритмы обеспечивают однозначность идентификации выявленных параметров деформирования исследуемой детали или конструкции (число и амплитуда нагружения, поврежденность материала по заданному критерию оценки реакции датчика в виде относительной площади «темных пятен»), повышают информативность и достоверность методов диагностики усталости по показаниям ДДИТ. Простота и вычислительная эффективность разработанных моделей обеспечивают обработку экспериментальных данных в режиме реального времени.

2. Модули пакета программ, разработанного автором на языке Visual С++ и оптимизированного по критерию минимума показателя вычислительной сложности, легко интегрируются в другие цифровые системы анализа и обработки изображений для решения практических задач оценки остаточной работоспособности деталей и конструкций машин по заданному критерию усталостного повреждения на основе показаний ДДИТ.

Реализация результатов исследования. Программный комплекс ImageAnalyser, включающий в себя разработанные математические модели, используется для обработки показаний ДДИТ и оценки усталостного повреждения элементов ходовой части легковых автомобилей марки ВАЗ-2110, ВАЗ-2111, ВАЗ-2112, проходящих регулярное техническое обслуживание и ремонт в «Автоцентре на Омской» (ЗАО «Апрель», г. Курган), что подтверждается соответствующим актом о внедрении. if

Заключение диссертация на тему "Математические модели и алгоритмы компьютерной обработки цифровых изображений реакции датчиков деформации интегрального типа"

4.4. Выводы по разделу

1. Подтверждена практическая эффективность построенной модели оценки величины накопленного усталостного повреждения в локальных зонах металлоконструкций в виде относительной площади "темных пятен",

14. определенной на основе компьютерной обработки цифровых изображений реакции ДДИТ при помощи программного комплекса 1та§еАпа1узег. Учитывая, что калибровка ДДИТ по относительной площади первых "темных пятен" на стальных образцах осуществляется при числах циклов нагружения на два-три порядка меньших, чем разрушение образцов от усталости, разработанные способы построения различных математических моделей позволяют реализовывать методы оценки усталостного повреждения и работоспособности конструкций на начальной стадии эксплуатации

2. Показано, что использование моделей топологического анализа скелетизированных изображений реакции ДДИТ, позволяет успешно решать задач идентификации условий циклического деформирования деталей и конструкций машин (изгиб, кручение, изгиб с кручением, растяжение).

3. Аналитическое выражение тарировочных зависимостей, связывающих микроструктуру материала с его физическими и механическими свойствами, дает возможность адекватно (с надежностью 0,95) интерпретировать экспериментальные данные и позволяет получить простой и легко воспроизводимый метод контроля работоспособности деталей и конструкций машин. Решение этих задач базируется на традиционном подходе -сопоставлении результатов тарирования ДДИТ и данных усталостных испытаний деталей с ДДИТ при одной и той же реакции датчика. Все построенные математические модели в соответствии с алгоритмом их реализации имеют место на заданном диапазоне тарирования ДДИТ.

4. Примеры практической реализации созданных методик подтверждают их эффективность и позволяют существенно расширить область функциональных возможностей использования ДДИТ при проведении экспериментальных исследований на основе разработанных моделей.

139

ЗАКЛЮЧЕНИЕ

В рамках проблемы по созданию программного комплекса, обеспечивающего цифровую обработку изображений и экспериментальных данных, полученных при помощи ДДИТ, в диссертационном исследовании были решены следующие задачи и сделаны выводы:

1. Осуществлена классификация изображений реакции ДДИТ и разработаны дифференцированные модели фильтрации, адаптированные к особенностям каждого класса изображений. Эти модели обеспечивают генерализацию изображений реакции /ЩИТ и калибровку ДДИТ по величине накопленного усталостного повреждения в локальных зонах металлоконструкций в виде относительной площади "темных пятен".

2. Разработана модель маркировки «темных пятен» бинарных изображений реакции ДДИТ. Промаркированное изображение является основой для подсчета ряда морфологических параметров (площадь, периметр, длина, ширина, координаты центра масс, моменты инерции) для формирования признакового пространства исследуемых объектов.

3. Модели линеаризации изображений реакции ДДИТ выделяют непрерывные контуры и скелеты информативных объектов, не нарушая их топологии. Разработанная модель анализа топологической структуры скелетизированных изображений определяет ориентацию "темных пятен" на поверхности ДДИТ и идентифицирует некоторые условия циклического деформирования образцов (изгиб, кручение, изгиб с кручением, растяжение).

4. Разработаны модели восстановлениям оптимальных тарировочных зависимостей в условиях малого объема экспериментальных данных (п < 20) с заданной надежностью 0,95. Аналитическое выражение зависимостей, связывающих микроструктуру материала с его физическими и механическими свойствами, представляет собой простой метод прогнозирования работоспособности деталей и конструкций машин. Этот метод базируется на традиционном подходе - сопоставлении результатов тарирования ДДИТ и данных усталостных испытаний деталей с ДДИТ при одной и той же реакции датчика. Число циклов нагружения образцов до реакции ДЦИТ составляет до 10 тыс. циклов, а прогнозируемое число циклов до поломки - сотни тысяч -десятки миллионов циклов, поэтому глубина прогноза превышает на два-три порядка интервал наблюдения.

5. Разработан и внедрен программный комплекс 1та§еАпа1узег, который, реализуя разработанные математические модели, объединяет в себе результаты проведенных теоретических и экспериментальных исследований и осуществляет научные исследования усталости на основе обработки цифровых изображений реакции ДДИТ. Реализация построенных математических моделей, представляет собой новую методику расшифровки реакции ДДИТ, обеспечивает алгоритмизацию и автоматизацию решения задач структурного и ситуационного анализа показаний ДДИТ.

4. Эффективность созданных моделей подтверждена серией стендовых экспериментов, проведенных для проверки адекватности разработанных математических моделей на примере образцов нефтегазового трубопроводного оборудования. Результаты испытаний подтверждают, что изложенная методика обработки показаний ДДИТ работоспособна и дает возможность решать задачи оценки состояния различных металлоконструкций в условиях фактического нагружения с минимальным количеством экспериментов (и < 20).

5. Разработанные математические модели и алгоритмы могут служить основой для разработки новых процедур использования ДДИТ в практике экспериментальных работ, учитывающих особенности функционирования исследуемых металлоконструкций и возможности получения информации, требуемой для расшифровки и анализа показаний датчиков. Идеология построения математических моделей легко может быть распространена на целый ряд подобных вопросов. В частности, разработанное программное обеспечение может быть использовано при решении вопросов обеспечения надежности сложных диагностических систем и комплексов. пг

Библиография Змызгова, Татьяна Рудольфовна, диссертация по теме Математическое моделирование, численные методы и комплексы программ

1. Алгоритмы и программы восстановления зависимостей / Под ред. Вапника В. Н. М.: Наука, Главная редакция физико-математической литературы, 1984. - 816 с.

2. Александров В.В., Булкин Г.А., Поляков А.О. Автоматизированная обработка информации на языке предикатов М.: Наука, 1982. - 240 с.

3. Александров В.В., Горский Н.Д. Алгоритмы и программы структурного метода обработки данных. JL: Наука, 1983. - 208 с.

4. Александров В.В., Шнейдеров B.C. Рисунок, чертеж, картина на ЭВМ. JL: Машиностроение, 1987. - 128 е.: ил.

5. Ахо А., Хонкрофт Дж., Ульман Дж. Построение и анализ вычислительных алгоритмов. М.: Мир, 1979. - 519 с.

6. Бойко В.И. О выборе материала, конструкции и технологии изготовления датчика усталостного повреждения // Завод, лаб. 1981. - №1, с.79-82.

7. Бойко В.И., Коваль Ю.Н. Анализ неразрушающих методов оценки усталостного повреждения металлов: Обзор.- Киев: Препринт АН УССР, 1982.-35 с.

8. Болотин В.В. Прогнозирование ресурса машин и конструкций. М.: «Машиностроение», 1984. - 312 с.

9. Борисенков Е.П., Романов М.А. Алгоритмы и программы статистической обработки информации. Ленинград: Гидрометеорологическое издательство (Гидрометеоиздат), 1969.-365 с.

10. П.Бутаков Е.А., Островский В.И., Фадеев И. Л. Обработка изображений на ЭВМ. М.: «Радио и связь», 1987. 240 с.

11. Быстрые алгоритмы в цифровой обработке изображений / Т. С. Хуанг, Дж-О. Эклунд, Г.Дж. Нуссбаумер и др.; Под ред. Т. С. Хуанга: Пер. с англ. М.: Наука, Радио и связь, 1984. - 224 е.: ил.

12. Виттих В.А., Сергеев В.В., Сойфер В.А. Обработка изображений в автоматизированных системах научных исследований. М.: Наука, 1982. 216 с.

13. Воскобойников Ю.Е., Преображенский Н.Г., Седельников А.И. Математическая обработка эксперимента в молекулярной газодинамике. -Новосибирск: Наука, 1984.-240 с.

14. Гренандер У. Лекции по теории образов: Синтез образов: Пер. с англ. М.: Мир, 1979. - 384 е.: ил.

15. Гренандер У. Лекции по теории образов: Анализ образов: Пер. с англ. -М.: Мир, 1981.-448 е.: ил. Ь

16. Гренандер У. Лекции по теории образов: Регулярные структуры: Пер. с англ. М.: Мир, 1983. - 448 с.

17. Григорович В.К. Твердость и микротвердость металлов. М.: Наука, 1976.-230 с.

18. Дворецкий В.И., Троценко Д.А., Мясоедов М.И. Исследование накопления усталостного повреждения в сварных соединениях с помощью медных гальванических пленок // Автоматическая сварка. 1982. - № 6. - с. 5-8.

19. Дель Г.Д. Определение напряжений в пластической области по распределению твердости. М.: Машиностроение, 1971. - 199 с.

20. Дрейпер Н., Смит Г. Прикладной регрессионный анализ. М.: Статистика, 1973. - 816 с.

21. Дуда Р., Харт П. Распознавание образов и анализ сцен: Пер. с англ. -М.: Мир, 1976.-512 е.: ил.

22. Елисеева И.И., Юзбашев М.М. Общая теория статистики: Учебник / Под ред. чл.-корр. РАН И.И. Елисеевой. М.: Финансы и статистика, 1996. -368 е.: ил.

23. Иосилевич Г.Б. Концентрация напряжений и деформаций в деталях машин. М.: Машиностроение, 1981. - 224 с.

24. Калинина В.Н., Панкин В.Ф. Математическая статистика: учеб. для студ. сред. спец. заведений. 3-е изд., испр. - М.: Высш. шк., 2001. - 336 е.: ил.

25. Керимов З.Г., Мамед-Заде O.A. Применение метода гальванического меднения к оценке степени повреждаемости деталей машин при нестационарном нагружении // Машиноведение, 1976. № 4. - с. 94-96.

26. Керимов З.Г., Мамед-Заде O.A. Способы определения напряжений и деформаций гальваническим меднением // Проблемы прочности, 1977. № 5.- с. 116-119.

27. Компьютеры на СБИС:В2-х кн. Кн.2: Пер. с япон./ Макино Т, Охаси М., Докэ X., Макино К. М.: Машиностроение, 1984. - 312 с.

28. Коллакот Р. Диагностика повреждений: Пер. с англ. М.: Мир, 1989.-512 с.

29. Коллинз Дж. Повреждение материалов в конструкциях: анализ, предсказание, предотвращение: Пер. с англ.- М.: Мир, 1984. 624 с.

30. Контроль качества с помощью компьютеров. Пер. с япон. / Мотоока Т., Хорикоси X., Сакаути М и др.- М.: Мир, 1991. 224 е.: ил.

31. Копнов В.А., Тимашев С.А. Датчики усталостного повреждения (обзор). Екатеринбург: «Наука». Уральское отделение, 1992. - 86 с.

32. Корн Г., Корн Т. Справочник по математике для научных работников. Пер. с англ. М., 1968. 720 с.

33. Линдли К. Практическая обработка изображений на языке Си/ Пер. с англ. М.: Мир, 1996. - 512 с.

34. Лукомский Я.И. Теория корреляции и ее применение к анализу производства М.: Госстатиздат, 1961, 231 с.

35. Маленков А.И. Разработка методов диагностики и прогнозирования работоспособности деталей передач при циклическом нагружении: Дис. канд. техн. наук. Курган, 1996. - 187 с.

36. Марчук Г.И. Методы вычислительной математики. М.: Наука, 1989.-608 с.

37. Математика и САПР: В 2-х кн. Кн.1. Пер. с франц. / Шенен П., Коснар М., Гардан И. И др. М.: Мир, 1988. - 204 с.

38. Математика и САПР: В 2-х кн. Кн.2. Пер с франц. / Жермен-Лакур П., Жорж П.Л., Пистр Ф., Безье П. М.: Мир, 1989. - 264 с.

39. Математическая статистика: Учебник/Иванова В.М., Калинина В.Н., Нешумова Л.А. и др. 2-е изд., перераб. и доп. - М.: Высш.школа, 1981. -371 е.: ил.

40. Махутов Н. А., Пашков Ю.И. Применение механики разрушения для оценки трещиностойкости трубопроводов // Проблемы машиностроения и автоматизации. 1991. №1. - с. 43-52.

41. Методы вычислений на ЭВМ: Справочное пособие / Иванов В.В. -Киев: Наук, думка, 1986.- 584 с.

42. Методы и алгоритмы аппроксимации решений на неравномерной сетке: Методические рекомендации. Минск: Институт технической кибернетики АН БССР, 1991. - 142 с.

43. Методы и средства натурной тензометрии: Справочник/ Дайчик М.Л., Григоровский Н.И., Хуршудов Г.Х.-М.: Машиностроение, 1989 240 с.

44. Методы компьютерной обработки изображений / Под. ред. В.А. Сойфера. 2-е изд., испр. - М.: ФИЗМАТГИЗ, 2003. - 784 с.

45. Носач В.В. Решение задач аппроксимации с помощью персональных компьютеров. М.: МИКАП, 1994. - 382 с.

46. Окубо X. Определение напряжений гальваническим меднением / Пер. с японск. М.: Машиностроение, 1968. - 152 с.

47. Павлидис Т. Алгоритмы машинной графики и обработки изображений: Пер. с англ.- М.: Радио и связь, 1986. 400 е.: ил.

48. Пинскер И.Ш., Ракчеева Т.А. Алгоритм аппроксимации гладкой функции двух переменных и анализ его погрешности.- В кн.: Поиск зависимости и оценка погрешности. М.: Наука, 1985. с. 32-50.

49. Поллард Дж. Справочник по вычислительным методам статистики / Пер. с англ. М.: Финансы и статистика, 1982. - 344 с.

50. Поляков А.Ю., Брусенцев В.А. Методы и алгоритмы компьютерной графики в примерах на Visual С++, 2-е изд., перераб. и доп. СПб.: БХВ -Петербург, 2003. - 576 е.: ил.

51. Прочность и надежность нефтегазового оборудования: Сборник докладов 20-го межотраслевого семинара (19-?2 сентября 2001 г.). М.: ГУП НИКИЭТ, 2001. с. 161-178.

52. Разработка приложений на Microsoft Visual С++6.0. Учебный курс: Официальное пособие Microsoft для самостоятельной подготовки / Пер. с англ. М.: Издательско-торговый дом «Русская редакция», 2000. - 576 е.: ил.

53. Резисторы: справочник / Дубровский В.В. Иванов Д.М., Пратусевич Н.Я. и др. 2-е изд., перераб. и доп. - М.: Радио и связь, 1991. - 527 с.

54. Роджерс Д. Алгоритмические основы машинной графики: Пер. с англ. М.: Мир, 1989. - 512 е.: ил.

55. Роджерс Д., Адаме А. Математические основы машинной графики: Пер. с англ. М.: Машиностроение, 1980. - 240 с.

56. Розенберг А.Ю. Методы экспериментальной оценки нагруженности и долговечности зубчатых передач с помощью гальванических медных датчиков циклических деформаций: Дис. канд. техн. наук. Курган, 1985, 223 с. И

57. Ряхин В.А., Мошкарев Г.Н. Долговечность и устойчивость сварных конструкций строительных и дорожных машин. М.: Машиностроение, 1984.-232 с.

58. Садыков С.С., Кан В.Н., Самандаров И.Р. Методы выделения структурных признаков изображений. Ташкент: Фан, 1990. -104 с.

59. Семенков О.И., Абламейко C.B. Методы и алгоритмы обработки растровой графической информации. Минск: ИТК АН БССР, 1984. 116 е.,

60. Семенков О.И. и др. Методы обработки и формирования растровых изображений. Минск: ИТК АН БССР, 1986. 98 с.

61. Семенков О.И., Митяев И.С. Диагностический анализ систем проектирования. Минск: Наука и техника, 1975. 136 е.: ил.

62. Сызранцев В.Н. Методы экспериментальной оценки концентрации циклических деформаций и напряжений на поверхностях деталей машин: Учебное пособие. Курган: РИО КМИ, 1993. - 83 с.

63. Сызранцев В.Н., Голофаст C.JL^ Измерение циклических деформаций и прогнозирование долговечности деталей по показаниям датчиков деформаций интегрального типа. Новосибирск: Наука, 2004- 206 с.

64. Сызранцев В.Н., Голофаст C.JL, Маленков А.И., Колпакова К.В. Диагностика усталости несущих систем и элементов транспортных машин с помощью датчиков деформаций интегрального типа: Учебное пособие. -Курган: КГУ, 1996. 87 с.

65. Сызранцев В.Н., Голофаст C.JL, Сызранцева К.В. Диагностика нагруженности и ресурса деталей трансмиссий и несущих систем машин по показаниям датчиков деформаций интегрального типа. Новосибирск: Наука, 2004.-188с.

66. Сызранцев В.Н., Добрынько А.И. Методы прогнозирования долговечности деталей по показаниям датчиков деформаций интегрального типа: Учебное пособие. Курган: КМИ, 1993. - 103 с.

67. Сызранцева К.В. Методическое и программное обеспечение измерения напряжений в деталях машин датчиками деформаций интегрального типа: Дис. канд. тех. наук.- Курган, 1998. 154 с.

68. Тимашев С.А. Надежность больших механических систем. М.: Наука, 1982. 184с.

69. Томпсон Н. Секреты программирования трехмерной графики для Windows 95/ Пер. с англ. СПб.: Питер, 1997. - 325 с.

70. Троценко Д. А. Разработка метода количественной оценки накопления усталостных повреждений в сварных соединениях с помощью гальванодатчиков: Автореф. дис. канд. техн. наук. Челябинск, 1986. - 16 с.

71. Фень Юань. Программирование графики для Windows. СПб.: Питер, 2002. - 1072 е.: ил.

72. Фу К. Структурные методы в распознавании образов: Пер. с англ. -М.: Мир, 1977.-320 с.

73. Хардле В. Прикладная непараметрическая регрессия: Пер. с англ. М., Мир, 1993. 349 с.

74. Холзнер С. Visual С++ 6: учебный курс./ Пер. с англ. СПб.: Питер, 2001.-576 е.: ил. н>

75. Чавкин A.M. Методы и модели рационального управления в рыночной экономике. М.: Финансы и статистика, 2001, 320 с.

76. Чукин Ю.В. Структуры данных для представления изображений // Зарубежная радиоэлектроника. 1983. №8, с. 87-103.

77. Черносвитов A. Visual С++ 7: учебный курс. СПб.: Питер, 2002. -528 е.: ил.

78. Шлихт Г.Ю. Цифровая обработка цветных изображений. М.: ЭКОМ, 1997.-336 с.

79. Энджел Э. Интерактивная компьютерная графика. Вводный курс на базе OpenGL К.: Вильяме, 2001. - 592 с.

80. Экспериментальные методы исследования деформаций и напряжений: Справочное пособие / Б.С.Касаткин, А.Б.Кудрин, Л. М. Лобанов и др. Киев: Наукова думка, 1981. - 584 с.

81. Ярославский Л.П. Введение в цифровую обработку изображений. -М.: Сов. радио, 1979.-312 с.