МЕТОДЫ РАСЧЕТА, МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ И ЭКСПЕРИМЕНТАЛЬНОГО ИССЛЕДОВАНИЯ ПРИЕМНЫХ УСТРОЙСТВ АКУСТИЧЕСКОЙ ДИАГНОСТИКИ 01.02.04 - механика деформируемого твердого тела 05.13.18 - математическое моделирование, численные методы и комплексы программ АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора технических наук Ростов-на-Дону 2010

Авторефераты по всем темам >> Авторефераты по техническим специальностям

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

ШИХМАН ВЛАДИМИР МАРКОВИЧ

МЕТОДЫ РАСЧЕТА,

МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ

И ЭКСПЕРИМЕНТАЛЬНОГО ИССЛЕДОВАНИЯ

ПРИЕМНЫХ УСТРОЙСТВ АКУСТИЧЕСКОЙ ДИАГНОСТИКИ

01.02.04 - механика деформируемого твердого тела

05.13.18 - математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

доктора технических наук

Ростов-на-Дону

2010

Работа выполнена в Федеральном государственном образовательном учреждении высшего профессионального образования

лЮжный федеральный университет

Научные консультанты: доктор физико-математических наук,

профессор

Белоконь А.В.,

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

профессор

Наседкин А.В.

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

профессор

Панич А.Е.,

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

профессор

Соловьев А.Н.,

доктор технических наук,

Заслуженный деятель науки РФ,

профессор

Шаповалов В.В.

Ведущая организация: Кубанский государственный университет.

Защита состоится 10 ноября 2010 г. в 15.00 часов на заседании диссертационного совета Д 212.058.03 в Донском государственном техническом университете (ДГТУ) по адресу: 344010, г. Ростов-на-Дону, пл. Гагарина, 1, ауд. 252.

С диссертацией можно ознакомиться в библиотеке ДГТУ

Автореферат разослан " " 2010 г.

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

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

канд. физ.-мат. наук, доцент Кренев Л.И.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

К наиболее перспективным и универсальным методам неразрушающего контроля прочности и диагностики состояния ответственных объектов относятся акустические методы: акустическая эмиссия (АЭ), акустическое течеискание, вибродиагностика и другие. Суть этих методов заключается в приеме и анализе акустических сигналов, генерируемых в процессе развития дефектов, истечения жидкостей или газов при испытаниях или функционировании объектов. Распространяясь в изделии, акустические сигналы достигают поверхности и могут быть зарегистрированы приемными устройствами в виде различных типов волн.

Преобразование механических колебаний в электрические сигналы осуществляется пьезоэлектрическими преобразователями. Часто, например, при контроле нагретых объектов: корпусов и толстостенных трубопроводов атомных реакторов, а также при наличии агрессивной или взрывоопасной среды, невозможен прямой контакт преобразователя с контролируемым объектом. В этих случаях применяют упругие волноводы, контактирующие с объектом контроля и выводящие приемные преобразователи из зоны нагрева. Широкое применение волноводы находят и в экспериментальных исследованиях приемных устройств. Пьезоэлектрический преобразователь или систему преобразователь с упругим волноводом будем называть приемным устройством. Волноводами являются также элементы изделий, в которых распространяются упругие колебания.

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

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

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

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

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

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

Объектом исследований является приемный тракт и его математические модели, которые включают: объект контроля, в котором распространяются акустические сигналы; источник акустического сигнала, заглубленный в объект контроля, приемное устройство, состоящее из акустического волновода и пьезоэлектрического преобразователя, устанавливаемое на поверхности изделия.

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

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

Достоверность результатов диссертационной работы основана:

на сопоставлении расчетных (аналитических и МКЭ) результатов с натурными экспериментами;

использовании строгого аппарата математической теории упругости и электроупругости; проверке сходимости результатов аналитических методов; исследовании сходимости численных методов в зависимости от степени дискретизации;
сопоставлении результатов, полученных независимым использованием аналитических методов механики и МКЭ;
использовании разработанных приемных устройств в практической диагностике объектов.

На защиту выносятся:

Разработанная в рамках метода однородных решений модель, включающая:

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

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

- исследование волновых полей в волноводах приемных устройств, контактирующих с объектом контроля, включая определение напряженного состояния во всей области торца и анализ эффективности приема волн Релея и Релея-Лэмба волноводом приемного устройства в зависимости от частоты и типа принимаемых волн.

2. Разработанные конечно-элементные модели и программный инструментарий для задач акустической диагностики, в том числе:

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

- расчет приемных преобразователей с учетом их взаимодействия с объектом контроля, анализ изменений АЧХ в результате этого взаимодействия;

- создание модели, включающей пьезопреобразователь, контактирующий с объектом контроля (ОК), и источник акустического излучения, заглубленный в ОК;

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

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

3. Создание экспериментальных методик исследования преобразователей для акустической диагностики и результаты разработки ряда преобразователей, в том числе:

- методики для исследования АЧХ преобразователей в условиях, приближенных к условиям реальной работы с учетом влияния параметров ОК, и контроль адекватности расчетных моделей реальным объектам на основе этой методики;

- преобразователей для раздельного приема составляющих волнового поля.

4. Комплекс методов и программ, дающих возможность более точного расчета приемных устройств акустической диагностики с заданными параметрами, в том числе при динамическом взаимодействии с ОК.

Научная новизна работы:

Впервые:

Создана численно-аналитическая модель и проведены исследования акустического и приемного трактов систем акустической диагностики в рамках строгих аналитических методов теории упругости и численных методов математического моделирования.
Разработан строгий подход к решению динамических задач теории упругости для полубесконечных волноводов при изгибных колебаниях с произвольными условиями на торце; для смешанной задачи предложено два подхода к учету особенностей в угловых точках;
В динамической задаче для упругих волноводов, сцепленных с объектом контроля, построено полное решение, включая определение напряжений в области контакта, и исследована эффективность приема в зависимости от частоты и типа набегающих волн.
При использовании метода однородных решений связанной теории электроупругости для свободных и вынужденных колебаний пьезокерамических плит исследованы новые особенности краевого и толщинного резонансов.
В рамках конечно-элементного подхода при математическом моделировании задач акустической диагностики использованы методы связанного динамического анализа задач пьезоэлектричества, методы учета затухания, внешних электрических цепей и моделирования не отражающих (лбесконечных) границ контролируемого объекта.
Разработаны конечно-элементные модели и программный инструментарий для расчетов приемного тракта акустических диагностических систем, выявлены закономерности, определяющие АЧХ приемных устройств.
Разработана модель и проведен расчет АЧХ приемных преобразователей с учетом их взаимодействия с ОК, исследованы изменения АЧХ в результате этого взаимодействия.
Сопоставлены для задач акустической диагностики результаты расчета пьезоэлектрических преобразователей методом однородных решений, методом конечных элементов и экспериментальные результаты.
Разработаны варианты конструкций селективных преобразователей.

Научная новизна результатов представленной диссертационной работы подтверждена 4 авторскими свидетельствами и 1 патентом.

Практическая значимость результатов работы состоит в повышении эффективности акустической диагностики за счет совершенствования модели приемного тракта; повышения технических характеристик приемных устройств путем оптимизации параметров (АЧХ, коэффициентов электромеханической связи, динамического согласования с ОК) для широкого ряда преобразователей из наиболее важных отечественных пьезоматериалов; создания экспериментальных методик исследования АЧХ преобразователей в условиях, приближенных к условиям контроля; создания на основе проведенных исследований ряда специализированных приемных устройств для акустической диагностики.

Основные результаты диссертации получены при выполнении НИР по следующим темам: Исследовать закономерности распространения и приема акустических сигналов в элементах изделий и создать приемные преобразователи с минимальным искажением информативных параметров акустической эмиссии. Регистрац. номер: 01813012227 М-59. Руководитель - Шихман В.М.; Разработать методы математического моделирования и оптимизации параметров пьезоэлектрических преобразователей для систем диагностики. Регистрац. номер: № 01880030093 М-89. Руководитель - Белоконь А.В.; Разработка теории и методов расчета вынужденных колебаний пьезоэлектрических преобразователей, используемых в качестве приемников акустических сигналов в устройствах неразрушающего контроля. Регистрац. номер: 01890075261 М-100. Руководители - Устинов Ю.А., Шихман В.М.; Создание методик, математическое и физическое моделирование излучения, распространения и приема упругих колебаний и разработка ряда преобразователей для систем акустико-эмиссионной диагностики с заданными амплитудно-частотными характеристиками, в том числе для повышенных температур. Регистрац. номер: 01980008200 ПМ-63. Руководитель - Шихман В.М.; Теоретическое и экспериментальное моделирование волновых процессов акустической диагностики прочности и создание специализированных высокоэффективных приемных устройств для различных частотных диапазонов. Регистрац. номер: 02930004033 М-11. Руководитель - Шихман В.М.; Разработка технических средств автоматизированной акустической системы диагностики и контроля течей оборудования и трубопроводов АЭС с РУ РБМК. Программа сотрудничества Мин. образования РФ и Мин. РФ по атомной энергии по направлению Научно-инновационное сотрудничество. Регистрац. номер: 3.2.03-20. Руководитель - Белоконь А.В.; Разработка программно-технического и методического комплекса акустической системы контроля состояния энергетического оборудования и трубопроводов АЭС. Программа сотрудничества Мин. образования РФ и Мин. РФ по атомной энергии по направлению Научно-инновационное сотрудничество. Регистрац. номер: 3.3.03-01. Руководитель - Белоконь А.В.; Разработка амплитудно-частотного метода и аппаратно-программных средств акустической системы диагностики и мониторинга течей оборудования и трубопроводов АЭС. Программа сотрудничества Мин. образования РФ и Мин. РФ по атомной энергии по направлению Научно-инновационное сотрудничество. Регистрац. номер: 3.4.03-09. Руководитель - Белоконь А.В.; Грант РФФИ Решение динамических задач и анализ волновых полей в контактирующих упругих телах канонической формы. Регистрац. номер 94-01-00220. Руководитель - Шихман В.М.; Грант Министерства образования РФ Разработка фундаментальных основ для совершенствования приемных устройств и методик акустико-эмиссионного контроля ответственных объектов железнодорожного и авиационного транспорта. Шифр гранта Т02-13.0-3841. Руководитель - Наседкин А.В.; Грант РФФИ Математическое моделирование и компьютерный дизайн новых видов активных композиционных материалов и устройств на их основе. Регистрац. номер 09-01-00875. Руководитель - Наседкин А.В.

Автор принимал непосредственное участие в указанных НИР в качестве исполнителя, ответственного исполнителя или научного руководителя работ.

Реализация результатов: Разработанные преобразователи использовались в НПО Молния для диагностики прочности крепления теплозащиты космического аппарата Буран, а также для акустико-эмиссионного контроля корпуса реактора ВВЭР-1000 на Нововоронежской АЭС. В ФГУП СибНИА им. С.А. Чаплыгина полосовые преобразователи акустической эмиссии внедрены при проведении акустико-эмиссионной диагностики состояния тонкостенных (авиационных) и литых (железнодорожных) конструкций. На объектах энергетики преобразователи использовались при испытаниях: паропроводов высокого давления и деаэраторов. ООО ИНТЕРЮНИС внедрил в системах диагностического мониторинга целостности горных пород низкочастотные преобразователи акустической эмиссии (ПАЭ) со встроенными предварительными усилителями. Использование разработанных преобразователей позволило увеличить расстояние между ПАЭ на объекте мониторинга и сократить необходимое количество измерительных каналов системы мониторинга. Высокотемпературные приемные устройства, работоспособные до 250 оС, использовались ЗАО НПФ Диатон для контроля элементов оборудования первого контура реактора Нововоронежской АЭС. В результате повышена достоверность диагностики трубопроводов первого контура и патрубков технологических каналов реактора. Указанные приемные устройства в составе многоканальной АЭ системы были переданы во ВНИИАЭС (г. Москва) для проведения работ по диагностике энергетического оборудования атомных электростанций. Разработанные пьезоэлектрические преобразователи нашли применение как первичные чувствительные элементы в системах регистрации выноса твердых фракций из газовых скважин, разработанных на ЗАО Объединение БИНАР г. Саров. С их использованием стало возможным построение передаточных функций систем. В ООО НПП "УЛЬТРАТЕСТ" разработанные преобразователи применяются для провендения АЭ контроля объектов экспертизы промышленной безопасности на предприятиях химической и пищевой промышленности, а также предприятиях других отраснлей, что обеспечило высокую достовернность получаемых результатов. Полосовые герметичные ПАЭ, а так же герметичные ПАЭ со встроенной электроникой применяются в Объединенном институте высоких температур РАН для измерения акустического поля при физическом моделировании электромагнитного воздействия на напряженные горные породы и моделирования сейсмоэлектрических и электросейсмических эффектов в горных породах.

Апробация работы: Результаты докладывались и обсуждались на Региональной конференции Динамические задачи механики сплошной среды (Краснодар, 1988); 2-м Всесоюзном семинаре Полимерные и композиционные сегнето,-пьезо, -пироматериалы и электреты в ускорении научно-технического прогресса (Москва, 1989); 12-th World Conference on NonDestruct Testing (Amsterdam, 1989); 1-ом Всесоюзном совещании Диэлектрические материалы в экстремальных условиях (Суздаль, 1990); Научно-технической конференции стран СНГ Производство и надежность сварных конструкций (Калининград, 1993); 9-th International Symposium on Application of Ferroelectrics.- (Pennsylvania, 1994); World Conference Review of Progress in Quantitative NDE (Washington, 1995; Brunswick, 1996); 15-й Российской научно-технической конференции "Неразрушающие физические методы и средства контроля" (Москва, 1999); 20-й Уральской региональной конференции Контроль технологий, изделий и окружающей среды физическими методами (Екатеринбург, 2001); Международной научно-практической конференции УПьезотехника - 94, 95, 2002Ф; на семинаре Акустико-эмиссионный метод диагностики на железнодорожном транспорте (Санкт-Петербург, 2003); 2-й, 3-й Всесоюзной и Всероссийской конференций по теории упругости (Тбилиси, 1984; Азов, 2003); 3-й Научно-технической конференции Научно-инновационное сотрудничество (Москва, 2004); Всероссийской конференции по волновой динамике машин и конструкций (Н. Новгород, 2004); на семинарах кафедры математического моделирования Ростовского госуниверситета (2001, 2004, 2009); 1-ой, 4-ой, 6-ой, 8-ой, 10-ой и 11-ой Международных конференциях Современные проблемы механики сплошной среды (Ростов-на-Дону, 1995, 1998, 2000, 2002, 2006, 2007); 5-ой Межд. научно-практич. конф. "Исследование, разработка и применение высоких технологий в промышленности" (Санкт-Петербург, 2008); 16-той международной конференции Современные методы и средства неразрушающего контроля и технической диагностики (Ялта, 2008).

В полном объеме диссертация заслушана и одобрена на заседаниях кафедр математического моделирования ЮФУ (Ростов-на-Дону, 2009) и кафедры информационных технологий ДГТУ (Ростов-на-Дону, 2010).

Публикации. По теме диссертации опубликовано работ 89, из них 16 в научных журналах по списку ВАК, получено 4 авт. свидетельства и 1 патент. В автореферате приведен список основных публикаций из 46 наименований. Большинство работ выполнено в соавторстве, в них автору принадлежит: конкретная постановка граничных задач, анализ и обобщение результатов; реализация метода однородных решений для динамических задач теории упругости с особенностями напряженного состояния в угловых точках; решение и анализ задач теории упругости для волноводов, контактирующих с объектом контроля, в режиме приема набегающих на область контакта волн Релея или волн Релея-Лэмба; разработка конечно-элементных моделей и программного инструментария для задач акустической диагностики, в том числе, расчет преобразователей с учетом их взаимодействия с объектом контроля и модель, включающая пьезоприемник, контролируемый объект и источник акустического излучения; создание экспериментальных методик исследования приемных устройств, планирование и участие в экспериментах, сопоставление экспериментальных результатов с теоретическими; практическое приложение результатов моделирования для разработки приемных устройств.

Структура и объем работы. Диссертационная работа состоит из введения, шести глав, заключения, списка литературы из 296 наименований и приложений. В приложении приведены 6 актов внедрения и использования. Материал изложен на 345 страницах, содержит 113 рисунков и 18 таблиц.

СОДЕРЖАНИЕ РАБОТЫ

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

Основными научными направлениями диссертационной работы являются строгие аналитические методы теории упругости и электроупругости и численные методы математического моделирования приемных устройств акустической диагностики с учетом их взаимодействия с объектами контроля.

Значительный вклад в указанную тематику внесли С.М.Айзикович, О.К.аАксентян, В.М.Александров, Н.П.Алешин, Н.Х.Арутюнян, В.А.Бабешко, В.М.Баранов, В.А.Белоконь, А.Н.Бескопыльный, Г.А.Бигус, Ю.И.аБобровицкий, М.В.Богуш, Н.В.аБоев, Бородин Ю.П., С.И.Буйло А.О.Ватульян, Л.П.аВовк, И.И.Ворович, Л.А.Галин, И.П.Гетман, Е.В.Глушков, Н.В.Глушкова, А.И.Гневко, А.Г.Горшков, И.В.Горячева, Д.В.Грилицкий, Г.А.Гринберг, В.Т.Гринченко, В.К.Доля, В.Г.Диденко, В.А.Еремеев, В.И.Иванов, Н.Ф.Ивина, Р.-Й.Ю.аКажис, И.Г.Кадомцев, В.В.Калинчук, Д.Н.Карпинский, В.В.Клюев, Е.В.Коваленко, Н.Н.Колоколова, В.В.Копасенко, А.С.Космодамианский, Ю.А.аКрамаров, Л.И.Кренев, Б.А.Кудрявцева, В.Д.Купрадзе, И.Я.Кучеров, А.И.Лурье, В.В. Мелешко, Н.Ф.аМорозов, В.И.Моссаковский, А.В.Наседкин, Б.М.Нуллер, М.Б.Оразов, А.Е.Панич, В.З.Партон, С.П.Пельц, Г.Я.Попов, О.Д.Пряхина, В.М.Сеймов, М.Г.Селезнев, И.Т.Селезов, А.С.аСкалиух, Б.И.Сметанин, Б.В.Соболь, А.Н.Соловьев, Л.Н Степанова, М.А.аСумбатян, С.А.Тарараксин, А.Н.Тихонов, А.С.Трипалин, А.Ф.Улитко, Ю.А.Устинов, Я.С.Уфлянд, В.Г.Харебов, М.И.аЧебаков, И.И.Шемякин, G.G.Adams, D.Berlincourt D.B.Bogy, E.P.Eer Nisse, A.G.Evans, T.D.Gerhardt, G.M. Gladvell, R.D.Gregory, R.аHolland, W.P.Mason, R.D.аMindlin, G.Muravin, M.Onoe, H.S.Paul, Shun Cheng, H.F.аTiersten, P.J.Torvik, Wu C.H. и др.

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

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

Во второй главе с помощью метода однородных решений теории связанной электроупругости, развитой Ю.А.Устиновым и В.В.Мадорским, исследовался важнейший элемент приемных устройств - круглая пьезокерамическая плита. В цилиндрической системе координат r, , z рассмотрены свободные и вынужденные колебания пьезокерамической плиты толщиной 2h и радиуса a, с электродированными торцами. Плита предварительно поляризована по толщине вдоль оси z.

Построенное решение использовалось для исследования собственных осесимметричных колебаний круглых плит из промышленной пьезокерамики ЦТС-19 и пьезокерамики ПКР-1 разработки РГУ (ЮФУ), широко применяемых в приемных преобразователях.

Исследованы важнейшие характеристики, широко применяемые при разработке преобразователей: частотные спектры для свободных колебаний пьезокерамической плиты из пьезокерамик ЦТС-19 и ПКР-1 в режимах короткого замыкания (к.з.) и холостого хода (х.х.) и вынужденные колебания в случае электрической гармонической нагрузки, задаваемой внешним источником. На рис. 1, в качестве примера, представлена зависимость амплитудных нормированных значений модуля тока от частоты для относительного радиуса R0 = 5,516 (R0 = a/h) для ЦТС-19. Сплошной линией показаны теоретические значения, пунктирной - экспериментальные, , - частота, к - плотность, с - нормирующий коэффициент.

Рис. 1

Рассчитанные спектры показали хорошее совпадение с экспериментальными исследованиями для пьезокерамической плиты с произвольным соотношением диаметра и толщины. Максимальная погрешность составляет 5,7 % при R0 = 1,03.

Изучение спектров собственных частот пьезоэлектрических плит показывает ряд особых частотных диапазонов, для которых характерно наличие пологих участков со слабой зависимостью резонансной частоты от радиуса (плато). Первый из таких участков соответствует частоте краевого резонанса. Выявлены некоторые дополнительные особенности краевого резонанса. Отмечено, что при изменении R0 и переходе к краевому резонансу изменяется фаза колебаний обеих составляющих вектора смещений, а в точке R0 = 3,114 вектор смещения уменьшается практически до нуля. При данном значении R0 и ток перестает носить резонансный характер. Аналогичные точки имеют место вблизи краевого резонанса и для других спектральных кривых. Таким образом, переход от формы колебаний, соответствующей какой-либо резонансной кривой, к краевому резонансу осуществляется через полное затухание колебаний диска и изменение их фазы.

Следующей областью со слабой зависимостью резонансной частоты от радиуса является область толщинного резонанса. Были исследованы формы колебаний в районе 3 на 10-ой и 11-ой гармониках на участках, на которых имеются плато. Во всех исследованных случаях поршневое движение отсутствует. Проведен анализ коэффициента электромеханической связи (КЭМС) по формуле Мэзона

(1)

где fr и fa - частоты резонанса и антирезонанса. Для области 2,96 < < 3,02 определенный по указанной формуле КЭМС оказывается незначительным: kd < 0,03. Отсюда следует, что эта область не соответствует толщинному резонансу, поскольку толщинный коэффициент связи для ЦТС-19 kt = 0,44.

а) R0 = 2,4; б) R0 = 3,114; н в) R0 = 4

Рис. 2

При 3,2 < < 3,7 также имеется тенденция к появлению горизонтальных участков на резонансных кривых. Форма перемещения приближается к поршневой, хотя и сохраняет некоторые неравномерности рельефа по радиусу. КЭМС здесь достигает 0,3. Именно этот резонанс, по нашему мнению, экспериментально идентифицируется как толщинный и имеет практическое значение.

Знание для каждой резонансной частоты эффективности преобразования механических колебаний в электрические имеет для разработки преобразователей столь же фундаментальное значение, как и знание резонансных частот. На примере пьезокерамики ЦТС-19 по формуле (1) проведен численный анализ зависимости kd от R0. В частности показано, что, начиная с третьей резонансной кривой, появляются точки, на которых КЭМС обращается в нуль. Эти точки соответствуют начальным участкам плато краевого резонанса, где происходит смена фазы колебаний.

Полученное теоретическое распределение механических и электрических полей в объеме плиты позволило провести анализ составляющих плотности внутренней энергии в пьезокерамическом диске и получить значения энергетического КЭМС kэн

Проведено сравнение результатов расчетов kэн, kd и экспериментальных значений kэкс. В большинстве исследованных случаев результаты достаточно близки, особенно при простых формах колебаний, что подтверждает обоснованность использования результатов расчетов kd по приведенной методике для практических целей.

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

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

Рассмотрена задача моделирования в ANSYS динамического режима работы составного пьезопреобразователя, нагруженного, в общем случае, на твердую деформируемую среду. В такой модели имеется одно или несколько пьезоэлектрических тел , n = 1,2Е и упругих сред , k =1,2,Е Предполагается, что волновые процессы в средах и описываются классическими линейными теориями пьезоэлектричества и упругости соответственно. Волновые поля в различных средах согласуются по граничным контактным условиям. Для расчетов в осесимметричной постановке вводилась цилиндрическая система координат r, , z. Для реализации пьезоэлектрических сред был выбран плоский четырехугольный элемент с опцией осесимметричности. При расчетах была использована каноническая сетка с необходимой плотностью конечных элементов, определяемой из апостеорного анализа погрешности. Демпфирование в пьезоэлектрике учитывалось только через упругие свойства, пьезоэлектрические и диэлектрические потери явно не учитываются, однако при этом с учетом связанности эффекты затухания проявляются и для электрических полей. Затухание в элементах системы задавалось с использованием усовершенствованной методики на основе учета демпфирования по Релею с учетом матриц масс и матриц жесткости сред через добротности Qj каждой среды. Рассматриваемые пьезоэлементы имели вид сплошных цилиндрических дисков, торцы полностью электродированы, и пьезоэлементы поляризованы по толщине вдоль оси z.

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

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

Для сопоставления результатов конечно-элементного моделирования, аналитических расчетов методом однородных решений и экспериментальных результатов исследованы собственные колебания пьезокерамической круглой плиты в режиме короткого замыкания. Максимальная погрешность (по отношению к эксперименту) - 5,7 % для аналитических значений и 5,5 % для конечно-элементных значений наблюдалась на первом резонансе при диаметре диска примерно равном толщине. Максимальное расхождение между двумя теоретическими методами равно 1,87 % при R0 = 5,16. В таблице для пьезокерамики ЦТС-19 приведены усредненные по разным резонансным частотам расхождения теоретических значений с экспериментальными для аналитического метода - ΔАср, такое же расхождение для МКЭ - ΔКср и расхождение теоретических частот, полученных разными методами между собой - Δср.

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

Таблица

R0	ΔАср (%)	ΔКср (%)	Δср (%)
1,032	3,7	3,4	0,25
2,563	2,7	2,5	0,21
5,156	2.3	3,2	1,1
10,450	1,7	2,0	0,35

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

Исследовано также распределение электрической индукции на торцах пьезокерамической плиты. Максимальное отклонение не превышает 2%.

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

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

Для исследования собственных частот преобразователя в режиме приема акустических сигналов разработана модель, в которой рассматривается преобразователь с установленной на его рабочей поверхности тонкой пьезопластиной, играющей роль излучателя. Данная модель использовалась и в расчетах и в экспериментах. Исследованы пьезоэлементы из ПКР-1 и ЦТС-83Г с излучателем из пьезокерамики на основе титаната висмута - ТВ. Пьезоматериалы ЦТС-83Г и ТВ разработаны и выпускаются НКТБ Пьезоприбор ЮФУ для работы при повышенных температурах. На рис. 3 и рис. 4 показаны результаты расчета АЧХ для пьезоэлемента из ПКР-1 диаметром R = 12 мм, толщиной h = 10 мм, нагруженного по плоской поверхности на излучатель в форме диска из ТВ толщиной hизл = 0,4 мм и радиусом, равным Rизл = 5 мм, а на рис. 4 - экспериментальные результаты для той же задачи. Ф - относительное значение потенциала. Отклонение в резонансных частотах не более 7 %.

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

Рис. 3. Рис. 4.

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

Реальный преобразователь всегда работает в контакте с объектом контроля (ОК), который оказывает значительное влияние на АЧХ. Метод конечных элементов позволяет моделировать такие сложные системы с учетом контакта приемника с ОК в пространственной осесимметричной постановке. Рассмотрены преобразователи, установленные на ОК в виде слоя толщиной h0 и полупространства, причем в обеих моделях присутствовал излучатель для возможности сопоставления с экспериментом. Для устранения отражений границы ОК в расчете покрывались демпфирующими слоями с теми же материальными свойствами, что и у объекта, но с высоким поглощением ультразвуковых колебаний (низкой добротностью). Это позволило добиться исчезновения резонансных явлений. При сравнении АЧХ преобразователя, установленного на ОК, с АЧХ такого же преобразователя при отсутствии ОК отмечено смещение первого резонанса в низкочастотную область, уменьшение добротности и общее расширение полосы пропускания. При увеличении толщины слоя и переходе к полупространству наблюдалось дальнейшее снижение частоты первого резонанса. Расчеты АЧХ при промежуточных значениях толщины слоя продемонстрировали плавность этих переходов. Экспериментальная проверка была проведена при h0 = 3,8 мм и для модели полупространства. На рис. 5 показана расчетная АЧХ приемника из пьезокерамики ПКР-1 (R = 6 мм, h = 10 мм) с ОК из стали, а рис. 6 демонстрирует экспериментальные результаты для той же задачи. Погрешность расчетов для резонансных частот не более 16%.

Рис. 5 Рис. 6

Описанная методика применена также к исследованию ВТП, установленных на полупространство. На рис. 7 представлены результаты расчета АЧХ и на рис. 8 - эксперимента для ВТП из пьезокерамики ЦТС- 83Г с R = 4 мм, h = 8.25 мм и Rпр = 8,2 мм. Максимальное расхождение не более 14%.

Дополнительные исследования показали, что провал на частоте 170 кГц связан с изгибными колебаниями протектора.

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

Рис. 7 Рис. 8

протектор, инертную массу и демпфер. Рассчитывается реакция преобразователя на принимаемый сигнал с целью анализа характеристик источника и оптимизации преобразователя для данного типа источника. Важной особенностью разработанной модели, как и предыдущей, является учет взаимодействия преобразователя с контролируемым объектом. Источник излучения представляет собой гармоническую сосредоточенную силу, заглубленную в изделие и действующую в направлении оси z (источник первого типа) или в перпендикулярном направлении - в направлении оси r (источник второго типа).

Были рассмотрены источники различного типа и расположенные на разной глубине. Исследования волновых полей на рабочей поверхности ОК из стали при отсутствии преобразователя показали, что источник первого типа преимущественно порождает продольные смещения на поверхности объекта (т.е. направленные вдоль оси z), а источник второго типа - смещения вдоль оси r, (т.е. поперечные смещения, лежащие в плоскости рабочей поверхности). Исследование пьезоэлементов разной толщины (от 1 до 15 мм) показало: колебания источника, действующие в направлении оси z, достаточно эффективно возбуждают резонансные радиальные колебания пьезоэлементов даже при малой относительной толщине (h/R) этих пьезоэлементов; увеличение толщины приводит к преимущественному приему продольных колебаний и расширяет полосу пропускания. Напряжение на выходе преобразователя из пьезокерамики ПКР-1 с диаметром элемента R = 5 мм, толщиной h = 15 мм при приеме сигналов из ОК при глубине источника Н = 40 мм показано на рис. 9 для источника с колебаниями, действующими в направлении оси z, и на рис. 10 - для источника с колебаниями в направлении оси r.

Можно отметить, что при переходе к колебаниям источника второго

Рис. 9 Рис. 10

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

Рассмотрены в качестве приемника преобразователи более сложного вида, содержащие, кроме пьезоэлемента, протектор, инертную массу и демпфер. Исследованы их АЧХ.

Проведено моделирование многослойных приемных устройств, содержащих дополнительные активные слои. Основное достоинство многослойных преобразователей заключается в повышении чувствительности при работе на длинный кабель. На примере преобразователя с 4-слойным активным элементом с чередующимися направлениями поляризации было подтверждено преимущество по чувствительности многослойного преобразователя перед соответствующим монолитным при работе на кабель с большой емкостью. Показано также отсутствие изгибных колебаний, ухудшающих АЧХ, и слабая зависимость резонансных частот от внешней емкости, что можно отнести к дополнительным преимуществам многослойного преобразователя перед монолитным. Исследована модель, в которой 4-слойный пьезоэлемент установлен на ОК с источниками, действующими в направлении оси z и в направлении оси r. В результате расчетов было выявлено, что в зависимости от частоты реакция на источник второго типа от 10 до 40 раз ниже, чем на источник первого типа. Подтверждено сделанное ранее заключение об относительном увеличении вклада высокочастотных составляющих при поперечных колебаниях источника.

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

Глава четвертая посвящена изучению в строгой постановке (в рамках метода однородных решений) возбуждения и распространения изгибных и продольных гармонических волн Релея-Лэмба в плоских волноводах приемных систем с силовой или кинематической нагрузкой на торце.

Упругая полуполоса (волновод) занимает область |x| 1, 0 y < ∞. Боковые грани полуполосы свободны от усилий. На торце y = 0 рассматривается два типа граничных условий: силовые и кинематические. Метод решения задач показан на примере изгибных колебаний. Решение для полуполосы строится в виде рядов по однородным решениям

(2)

где Cn - неизвестные коэффициенты разложения, определяемые из граничных условий, n - корни дисперсионного уравнения для изгибных или продольных волн Релея-Лэмба. При этом корни, по которым ведется суммирование, выбираются в соответствии с принципом энергетического излучения.

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

Проведен анализ динамического поведения полуполосы при заданном на торце кинематическом возбуждении. Исследовано два подхода к решению данной задачи, которая характеризуется наличием локальной особенности напряженного состояния в угловых точках. Поведение этой особенности исследовано О.К. Аксентян, Я.С. Уфляндом и другими и имеет вид

{, r} ~ r-γ, r → 0, (3)

γ находится из соотношения (ν - коэффициент Пуассона)

(3 - 4ν)sin20,5γ = (1 - γ) - (1 - 2ν)2, 0 < γ < 0,5. (4)

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

(5)

где =-(1-x), s1=0,5x - 0,5 + γ1, s2=0,5x Ц1,5 + γ1, s3=-0,5x - 0,5 + γ1.

Сопоставляя (5) и (3) для угловой точки, находим: γ1 = γ, где γ определяется из (4).

Анализируя напряженное состояние, можно отметить, что среди рядов (5), кроме сходящихся и дающих особенность лишь в угловой точке, есть также ряды, расходящиеся при y = 0 для области x ≥ 1 - 2γ. Для определения напряжений во всей области торца необходимо привлечь методы обобщенного суммирования расходящихся рядов, описанные Г. Харди. В результате, напряжения на торце представляются в виде конечных сумм, а напряжения вблизи угловой точки имеют асимптотическое поведение (3).

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

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

Для численных расчетов в качестве примера рассмотрен стальной волновод и исследована точность выполнения граничных условий при различных значениях x и зависимость погрешности в угловых точках от порядка системы на разных частотах. Погрешность снижается с уменьшением нормированной частоты = 2h/c2 (c2 = (μ/ρ)0,5) и ростом порядка редуцированной системы. Показано, что при исследованных формах заданных на торце напряжений амплитуды первых мод, соответствующих действительным корням, с ростом частоты убывают, и более эффективно возбуждаются высшие моды. Следовательно, для повышения эффективности волноводной приемной системы при высоких частотах нужно ориентироваться на прием высших мод. Показано размывание формы заданной на торце нагрузки с удалением от торца на низких и высоких частотах. Однако тип волновых процессов (изгибный или продольный) дает информацию о волновом процессе в объекте контроля и, в конечном итоге, об источнике излучения. Изучено волновое поле и его модовый состав при кинематическом нагружении. Показано, что как для сдвигового, так и продольного напряжения в угловых точках имеется характерная особенность.

Для исследования распространения акустических колебаний в волноводах применены также и численные конечно-элементные методы с ориентацией на программный инструментарий пакета ANSYS. Для моделирования упругих сред были выбраны четырехузловые прямоугольные конечные элементы. Рассмотрен упругий слой толщиной h конечной длины a, расположенный от x = 0 до x = a. Для устранения отражений от границы исследуемого объекта эта граница покрывалась демпфирующим слоем. Источник излучения моделировался гармонической силой, заданной на торце. Измерялась амплитуда вертикального перемещения uy на поверхности слоя (выходной сигнал). Внешние границы считались свободными.

Проведена экспериментальная проверка для стального слоя толщиной h = 3,8 мм. Сравнение теоретических и экспериментальных графиков подтверждает эффективность модели, особенно хорошо совпадает общая картина резонансных пиков. Показано, что метод конечных элементов при надлежащих конечно-элементных сетках обеспечивает достаточную для практических целей точность расчетов акустических полей в волноводах.

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

В первой задаче рассмотрено взаимодействие волны Релея, бегущей вдоль границы контролируемого объекта (полуплоскости), с упругим волноводом (полуполосой), сцепленным торцом с полуплоскостью. Волновод занимает область |х| 1, 0 y , полуплоскость - область y 0. На поверхности y = 0 задана набегающая на волновод волна Релея.

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

(6)

i = 1, 2 соответственно для полуполосы и полуплоскости. Показатель особенности γ0 находится из уравнения, полученного О.К. Аксентян.

Для решения используется метод, описанный в главе 4. Решение для волновода ищется в виде аналогичном разложению (2). Для определения напряжений во всей области контакта использовано обобщенное суммирование рядов. Напряжения вблизи угловой точки области контакта при y = 0 имеют вид

Отсюда видно, что (i)(ρ,) дают требуемую особенность (6). Решение задачи сводится к решению конечной алгебраической системы относительно коэффициентов разложения в рядах (2).

Во второй задаче исследовался прием волноводом волны Релея-Лэмба, распространяющейся в бесконечной полосе. Общий ход решения подобен решению задачи о приеме волны Релея. В итоге определено наряженное состояние в области контакта и получена алгебраическая система уравнений для нахождения постоянных в рядах (2).

Численные исследования поставленных задач проведены в условиях как одномодового, так и многомодового распространения в волноводе. Исследована сходимость решений при увеличении порядка системы, оценены погрешности по выполнению условий сопряжения перемещений в области контакта. Предложенная модель позволяет оценить эффективность возбуждения колебаний в волноводе для разных типов волн и частот. При приеме волны Релея на низкой частоте = 0,5, когда в волноводе существует одна распространяющаяся мода, отношение амплитуды продольного перемещения в волноводе к амплитуде соответствующего перемещения в падающей волне Релея v(m)/vR(m) = 0,42. Это свидетельствует о достаточно эффективном возбуждении первой моды на низких частотах. На высокой частоте при = 3,0 усредненная по сечению амплитуда vср2 = 0,46 при амплитуде вертикального перемещения в волне Релея vR2(m) = 2,8. Отношение средней по сечению амплитуды продольного перемещения в волноводе к амплитуде соответствующего перемещения в падающей волне Релея составляет vср2/ vR2(m) = 0,16. Из приведенных расчетов следует, что эффективность возбуждения на частоте = 3 в 2,6 раза меньше, чем на частоте = 0,5. Однако для расчета и оптимизации приемной системы в реальных условиях эксплуатации следует учесть частотную зависимость акустического шума в объекте контроля, если она имеется. В качестве примера рассмотрена задача оптимизации приемной волноводной системы для акустико-эмиссионной диагностики объекта, в котором с ростом частоты шумы уменьшаются обратно пропорционально частоте. Показано, что при диагностике указанного объекта эффективность приема волн Релея приведенная к шуму, для f = 600 кГц в 2,3 раза больше, чем для f = 100 кГц. В данном случае повышение эффективности приема с ростом частоты связано со спектральными особенностями акустических шумов объекта контроля.

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

В шестой главе описаны экспериментальные методы для проверки и дополнения результатов теоретического моделирования, изложены итоги практической реализации теоретических разработок, приводятся результаты совершенствования преобразователей в направлении оптимизации частотной характеристики, увеличения чувствительности и повышения селективности к рабочим типам колебаний. Экспериментальные исследования касаются, в основном, преобразователей акустической эмиссии (ПАЭ), однако большинство рассматриваемых вопросов имеют приложение и к другим акустическим методам.

Экспериментальная методика применяемая нами для исследования АЧХ преобразователей является дальнейшим развитием методики, рекомендуемой ГОСТ 23702-90 для преобразователей ультразвуковой дефектоскопии, и отличается от нее заменой диэлектрического излучателя тонким пьезокерамическим диском из анизотропной пьезокерамики. Методика отличается от аналогов возможностью прямого определения АЧХ чувствительности по свободному полю; проведением испытаний с акустической нагрузкой, имеющей геометрические и физические параметры, моделирующие реальные условия диагностики; возможностью конечно-элементного моделирования. Предложенный подход позволяет исследовать частотные характеристики ПАЭ раздельно на продольное и поперечное смещение.

Проведен сравнительный анализ перспективных пьезоматериалов, разработанных в Ростовском университете (ЮФУ) и описаны результаты макетирования широкого ряда преобразователей. Приведены примеры разработки преобразователей, основанные на использовании аналитических и численных методов расчета, и результаты сопоставления с экспериментом.

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

Изложены результаты разработки преобразователей для раздельного приема составляющих волнового поля (селективных преобразователей). Цель этой разработки - повышение акустической помехозащищенности и информативности контроля. Разработаны экспериментальные методики и исследованы характеристики селективных преобразователей. Исследование приемников продольных колебаний основано на применении упругих волноводов. Одна из конструкций преобразователя для приема продольных колебаний на основе ХГС-2 защищена патентом. Подтверждена высокая степень селективности преобразователя на основе ХГС-2. АЧХ селективных сдвиговых преобразователей исследованы при работе на поверхности массивной стальной плиты. Для устранения диаграммы направленности, отличающейся от круговой, предложены, защищенные авторским свидетельством, конструкции преобразователя, работающего на сдвиговых колебаниях, но имеющего круговую диаграмму направленности. Предложена конструкция преобразователя с круговой диаграммой направленности для раздельного приема двух компонент перемещения. Экспериментально исследован прием поверхностных волн Релея различными преобразователями.

Созданы высокотемпературные преобразователи различных температурных диапазонов, в том числе и до 500 0С. Описаны их конструктивные особенности и исследованы температурные характеристики.

Проведены сравнительные испытания разработанных преобразователей и лучших зарубежных аналогов, включая датчики ведущей в области АЭ фирмы Physical Acoustic Corporation (США). Полученные результаты показывают, что разработанные нами преобразователи не уступают по основным параметрам приемникам известных зарубежных фирм.

Заключение. В заключении кратко изложены основные результаты диссертации.

Основные результаты работы состоят в следующем:

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

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

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

1.3. На основе метода однородных решений связанной электроупругости проведен анализ свободных и вынужденных осесимметричных колебаний в пьезокерамической круглой плите в различных режимах электрического возбуждения, рассчитаны спектры, исследованы новые особенности краевого и толщинного резонансов, характеризующихся высокой стабильностью резонансных частот при изменении радиуса, изучены КЭМС для ряда низших и высокочастотных мод колебаний. Разработаны подходы к созданию приемных устройств резонансного и полосового типов на основе модели приемной системы, содержащей преобразователь из текстурного пьезоматериала и упругий цилиндрический волновод.

2. Созданы модели, методики расчета и разработан необходимый программный инструментарий для прямых численных методов расчета задач акустической диагностики на основе МКЭ для связанного динамического анализа задач пьезоэлектричества, теории упругости с учетом демпфирования, внешних электрических цепей и моделирования не отражающих (лбесконечных) границ ОК.

2.1. Разработаны подтвержденные экспериментами конечно-элементные модели и методики расчета пьезоэлектрического преобразователя, представленного многослойным пакетом, содержащим в общем случае протектор, многослойный пьезоэлемент, инертную массу и демпфер.

2.2. Созданы конечно-элементные модели расчета приемных преобразователей с учетом их взаимодействия с объектом контроля (ОК), которые продемонстрировали принципиальные изменения АЧХ преобразователей в результате этого взаимодействия: снижение частоты первого резонанса, существенное уменьшение добротности и значительное общее расширение полосы пропускания, зависящие, в частности, и от параметров ОК. Адекватность модели подтверждена экспериментальными исследованиями.

2.3. Разработана конечно-элементная модель приема колебаний от источников акустического излучения в объектах контроля, включающая пьезопреобразователь, контролируемый объект и источник, представляющий собой гармонические силы, заглубленные в ОК и действующие вдоль оси z или оси r. Обнаружено, что при переходе от источника с колебаниями вдоль оси z к источнику с колебаниями вдоль оси r наблюдается смещение спектра в сторону высоких частот, что дает дополнительную информацию о параметрах и динамике источника.

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

4. Создана методика экспериментального определения АЧХ чувствительности преобразователей с учетом взаимодействия преобразователя с объектом контроля заданной геометрии, позволяющая исследовать частотные характеристики ПАЭ раздельно на продольное и поперечное смещение.

5. На основе теоретического моделирования разработан ряд ПАЭ, в том числе: резонансные, полосовые, высокотемпературные, описаны их конструктивные особенности и АЧХ. Созданы преобразователи для раздельного приема составляющих волнового поля, конструкции которых защищены авторским свидетельством и патентом. Приведены примеры разработки преобразователей с заданными характеристиками с использованием аналитического и конечно-элементного моделирования. Анализ технической документации и сравнительные испытания показали, что преобразователи нашего производства не уступают по основным параметрам приемникам известных зарубежных фирм.

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

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

СПИСОК ОСНОВНЫХ ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ

Шихман В.М., Трипалин А. С., Бешенцев В. Д. и др. Исследование медленного роста трещин методом акустической эмиссии // Изв. СКН - ВШ, Естеств. науки - 1982. № 2. - С. 3-6.
Дорошенко В.А., Трипалин А.С., Шихман В.М. Исследование сигналов акустической эмиссии в процессе переполяризации сегнетоэлектриков // Изв. СКН - ВШ. Естеств. науки. - 1980. № 3. - С. 8-11.
Пельц С.П., Трипалин А.С. Шихман В.М. Расчет акустического поля в упругом волноводе конечной длины, нагруженном на торцах // Изв. Сев.-Кавк. научн. центра Высш. школы. Естеств. науки. 1984. № 3. - С. 7-10.
Пельц С. П., Трипалин А. С., Шихман В.М. Вынужденные колебания полубесконечного волновода, нагруженного на торце, при неразрушающем контроле изделий энергомашиностроения // Изв. Сев.-Кавк. научн. центра Высш. школы. Естеств. науки. - 1985. № 4. - С. 5-8.
Пельц С.П., Шихман В.М. О сходимости метода однородных решений в динамической смешанной задаче для полуполосы // Докл. АН СССР. - 1987. Т. 295. № 4. - С. 821-824.
Пельц С.П., Шихман В.М. Распространение волн в крестообразном соединении бесконечных упругих полос // ПММ. - 1987. - Т. 51. Вып. 1. - С. 54-59 .
Пельц С.П., Шихман В.М. Рассеяние волны Релея на упругой полуполосе, сцепленной на торце с упругой полуплоскостью // ДАН СССР. - 1987. Т. 292. № 2. - С. 299-303 .
Жиров В.Е., Захарова С.В., Трипалин А.С., Шихман В.М. Расчет активных элементов преобразователя для акустико-эмиссионного контроля атомных реакторов // Изв. СКН - ВШ. Естественные науки. - 1990. № 2. - С. 76-81.
Захарова С.В., Шихман В.М. Вынужденные изгибные колебания упругой полуполосы при смешанных граничных условиях. // ПММ. - 1996. № 2.- С. 340-345.
Наседкин А.В., Шихман В.М., Захарова С.В., Иванилов И.В. Применение конечно-элементных методов для расчета приемных систем акустико-эмиссионного контроля // Дефектоскопия. 2006. № 2. - С. 16 - 27.
Шихман В.М. Взаимодействие поверхностной волны с упругим препятствием // Экологический вестник научных центров Черноморского экономического сотрудничества. 2008. № 3. - С. 53-57.
Шихман В.М. Распространение волн Релея-Лэмба в Т-образном соединении упругих волноводов // Известия вузов Сев.-Кавк. региона. Естеств. науки. 2009. № 1. - С. 35 - 37.
Трипалин А.С., Шихман В.М. Ряд пьезоэлектрических преобразователей для приема сигналов акустической эмиссии // Автоматическая сварка, 1984. № 5. - С. 33 - 37.
Трипалин А. С., Шихман В.М., Коваленко В.И. и др. Акустическая эмиссия при малоцикловых испытаниях сварных тавровых элементов // Диагностика и прогнозирование разрушения сварных конструкций. - 1985. Вып.1. - С. 89-93.
Трипалин А.С., Шихман В.М., Карпенко С.Б. Расчет характеристик преобразователей из текстурных пьезоматериалов с волноводами // Техническая диагностика и неразрушающий контроль. Киев. 1989. № 3. - C. 52-56.
Шихман В.М., Гринева Л.Д. Преобразователи акустической эмиссии на основе современных пьезоматериалов // Техн. диагн. и неразр. контроль. Киев. - 1994. № 1. - С. 34-38.
Трипалин А.С. Кузьмин Г.А., Елецкий С.А., Коломиец и.В., Шихман В.М. Измерение параметров акустической эмиссии при гидравлическом нагружении корпуса реактора ВВЭР-1000 Нововоронежской АЭС // Энергомашиностроение. - 1982. № 3. - С. 33-36.
Трипалин А.С., Шихман В.М., Жиров В.Е., Захарова С.В. и др. Расчет пьезокерамических преобразователей методами трехмерной теории электроупругости // Акустическая эмиссия материалов и конструкций. Ч.1. Ростов-на-Дону.: Изд.-во РГУ, 1989. - С. 95 - 99.
Шихман В.М. Разработка преобразователей акустической эмиссии в НИИМ и ПМ РГУ (основные направления и результаты) // Современные проблемы мех. сплошной среды. Труды 6-ой Междун. конф. Ростов-на-Дону. 12-14 июня 2000 г. Ростов-на-Дону: Изд-во СКН - ВШ. 2001. Т.1. - С. 243-247.
Захарова С.В., Шихман В.М. Колебания толстых пьезокерамических дисков на высокочастотных гармониках. Современные проблемы мех. сплошной среды // Труды 4-ой Международной конференции. Ростов-на-Дону, 1998. Т. 1. - С. 160-164.
Иванилов И.В.,аНаседкин А.В.,аШихман В.М. Расчет АЧХ преобразователей акустической эмиссии в осесимметричной постановке методом конечных элементов // Сборник докладов международной научно-практической конф. УФундаментальные проблемы пьезоэлектрического приборостроенияФ. (Пьезотехника-2002). 17-21 сент. 2002 г. Тверь, 2002 г. - С. 164-170.
Богма К.К., Шихман В.М. Преобразователи акустической эмиссии. Инженерный анализ и исследования // Современные проблемы механики сплошной среды. Тр. VIII Межд. конф., г.Ростов-на-Дону, 14-18 октября 2002 г. Т. 1. / Ростов-на-Дону: изд-во УНовая книгаФ, 2003. С. 41 - 45. Труды 8-ой Международной конференции. г. Ростов-на-Дону., г. Ростов-на-Дону. С. 41-45.
Шихман В.М., Иванилов И.В. Моделирование приема упругих колебаний акустико-эмиссионным преобразователем методом конечных элементов // Труды III Всероссийской конференции по теории упругости. г. Азов. 13-17 октября 2003. - С. 406-409.
Наседкин А.В., Шихман В.М., Иванилов И.В. Конечно-элементное моделирование распространения и приема акустических сигналов от заглубленного источника в полубесконечной среде // Современные проблемы механики сплошной среды. Тр. VIII Межд. конф., г.Ростов-на-Дону. Т. 1. Ростов-на-Дону: изд-во УНовая книгаФ, 2003. С. 116-120.
Наседкин А.В., Шихман В.М., Иванилов И.В. Совершенствование приемных устройств и методик акустико-эмиссионного контроля объектов железнодорожного и авиационного транспорта на основе теоретического и экспериментального моделирования волновых процессов // Техника, технология и перспективные материалы. Межвуз. сб. научн. тр. - Москва, МГИУ. 2004. С. 251 - 256.
Захарова С.В., Наседкин А.В., Шихман В. М. Зависимость амплитудно-частотной характеристики преобразователей от граничных условий на торце // Современные проблемы механики сплошной среды. Труды Х Межд. конф., г.Ростов-на-Дону, 5-9 декабря 2006 г. Т. 1. Ростов-на-Дону: Новая книга, 2007. - С. 128-132.
Захарова С.В., Наседкин А.В., Шихман В.М. Исследование многослойных высокочувствительных акусто-эмиссионных датчиков // Современные проблемы механики сплошной среды. Тр. XI Межд. конф., г. Ростов-на-Дону, 26-29 ноября 2007 г. Т. 2. Ростов-на-Дону: изд-во "ЦВВР", 2008. С. 88-92.
Шихман В. М., Наседкин А.В., Захарова С.В. Применение конечно-элементных методов для расчета многослойных пьезопреобразователей в режиме приема акустических колебаний // Материалы 16-той международной конференции Современные методы и средства неразрушающего контроля и технической диагностики. Ялта, 1 - 5 октября, 2008. - С. 66 - 69.
Шихман В.М. Перспективные направления в разработке преобразователей акустической эмиссии // ХVI Российская н.-т. конференция УНеразрушающий контроль и диагностикаФ. - Москва. 28 июня - 2 июля 1999 г. Тез. докл. - Т.2. С. 111.
Шихман В.М. Расчет приемных устройств для акустических методов неразрушающего контроля // ХХ Уральская региональная конф. Контроль технологий, изделий и окружающей среды физическими методами. 15-16 мая 2001 г. Екатеринбург. Тез. докл. - С. 60-61.
Шихман В.М. Разработка нового поколения пьезоэлектрических преобразователей для акустико-эмиссионной диагностики // ХХ Уральская региональная конференция Контроль технологий, изделий и окружающей среды физическими методами. 15-16 мая 2001 г. Екатеринбург. Тез. докл. - С. 58-59.
Патент № 1784095 (СССР). Козинкина А.И., Трипалин А.С.и Шихман В.М. Пьезоэлектрический преобразователь для приема сигналов акустической эмиссии // Опубл. в Б.И. № 47. 1992.
А.с. 918284 (СССР). Способ поляризации пьезокерамических элементов / Трипалин А.С., Шихман В.М.. - Опубл. в Б.И. № 13.1982.
А.с. № 1167497 (СССР). Пьезоэлектрический преобразователь (его варианты) / Трипалин А.С., Шихман В.М.. - Опубл. в Б.И. № 26.1985.
А.с. № 1260849 (СССР). Ультразвуковой преобразователь для калибровки систем акустико-эмиссионного контроля. / Кузьмин Г.А., Крейнер Б.Я., Шихман В.М. Опубл. в БИ. № 4, 1985.
А.с. № 1320734 (СССР). Способ акустикоэмиссионного контроля изделий / Мельцер Я.Е., Трипалин А.С., Шихман В.М., Колесников В.И. - Опубл. в Б.И. № 24, 1987.
Shikchman V.M., Pelts S.P. Mathematical modeling of the Rayleigh wave reception by the system with elastic waveguide. - NY.: Plenum Press. QNDE. 1996. V. 15A. - P. 153-160.
Shikhman V.M., Grineva L.D. High-efficiency Transducers for Acoustic Diagnostics on the Bases of Novel Anisotropic Piezoceramic Materials // Proceedings of 9th IEEE International Symposium on Application of Ferroelectrics. Penn. State Scantion Сonference Center. Pensilvania. USA. Aug. 1994. - P. 172-174.
Tripalin A.S., Shikchman V.M. The modern mathematical models of piezoelectric transducers // Proceeding of 12th World Conference on NonDestruct. Testing. Amsterdam. April 23-28, 1989. V. 2. - P. 1501-1503.
Shikchman V.M. The Improvement of the Transducers for Acoustic Evaluation Methods // Review of Progress in Quantitative NDE. July 1995. Abstracts. University of Washington-Seattle. USA. - P. 151
Shikhman V. M. The Research of acoustic emission transducers, selective to the specific modes // 26 Annual Review of Progress in Quantitative Nondestructive Evaluation. July 25 - 30, 1999. Montreal Quebec, Canada. - P. 32.
Shikhman V.M., Grineva L.D. High-efficiency Transducers for Acoustic Diagnostics on the Bases of Novel Anisotropic Piezoceramic // Materials. Abstracts of 9th International Symposium on Application of Ferroelectrics. - Penn. State Scantion conference Center, Pennsylvania.- USA. Aug.7-10. 1994. - P. 39.
Shikchman V.M., Grineva L.D. Wide-band High-sensitive Transducers, Made of New Piezoelectric Materials // Review of Progress in Quantitative NDE. July 27. 1997. San Diego, California. USA. - P. 10.
Shikchman V.M., Pelts S.P. Mathematical Modeling of the Rayleigh Wave Reception by the System with Elastic Waveguide // Quantitative NDE Review of Progress in Quantitative NDE. University of Washington-Seattle. Abstracts. Washington. July 30. 1995. - P. 52
Shikchman V.M., Pelts S.P. Modeling of the Rayleigh-Lamb Wave Reception by the Elastic Waveguide // Review of Progress in Quantitative NDE. July28 - August 2. Brunswick, Maine. USA. 1996. - P. 46.
Shikhman V.M., Zakharova S.V. Analytical Approach to the Calculation of Three-Dimensional Transducers for the Acoustic Control Methods. Review of Progress in Quantitative NDE. July 28 - August 2, 1996. Brunswick, Maine, USA. P. 38

Авторефераты по всем темам >> Авторефераты по техническим специальностям

Blog

СПИСОК ОСНОВНЫХ ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ