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

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

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

АВТОРЕФЕРАТ

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

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

Работа выполнена в Федеральном государственном образовательном учреждении высшего профессионального образования «Южный федеральный университет» Научные консультанты: доктор физико-математических наук, профессор Белоконь А.В., доктор физико-математических наук, профессор Наседкин А.В.

Официальные оппоненты: доктор технических наук, профессор Панич А.Е., доктор физико-математических наук, профессор Соловьев А.Н., доктор технических наук, Заслуженный деятель науки РФ, профессор Шаповалов В.В.

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

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

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

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

Ученый секретарь диссертационного Совета канд. физ.-мат. наук, доцент Кренев Л.И.

ОБЩАЯ ХАРАКТЕРИСТ ИКА РАБОТЫ Актуальность темы. Возрастающая сложность объектов современной техники требует непрерывного совершенствования методов и средств диагностики их надежности, работоспособности, безопасности эксплуатации и прогнозирования критических состояний. В настоящее время проблема усугубляется в связи с выработкой ресурса работоспособности у значительной час ти действующих объектов машиностроения. Продление ресурса и экологическая безопасность эксплуатации становятся важнейшими экономико-социальными факторами.

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

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

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

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

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

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

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

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

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

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

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

- разработка необходимого программного инструментария для прямых численных методов расчета соответствующих задач теории упругости и электроупругости на основе метода конечных элементов (МКЭ);

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

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

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

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

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

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

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

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

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

- использовании строгого аппарата математической теории упругости и электроупругости;

проверке сходимости результатов аналитических методов;

исследовании сходимости численных методов в зависимости от степени дискретизации;

- сопоставлении результатов, полученных независимым использованием аналитических методов механики и МКЭ;

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Впервые:

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

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

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

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

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

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

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

7. Разработана модель и проведен расчет АЧХ приемных преобразователей с учетом их взаимодействия с ОК, исследованы изменения АЧХ в результате этого взаимодействия.

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

9. Разработаны варианты конс трукций селективных преобразователей.

Научная новизна результатов представленной диссертационной работы подтверждена 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 авт. свидетельс тва и патент. В автореферате приведен список основных публикаций из наименований. Большинство работ выполнено в соавторс тве, в них автору принадлежит: конкретная пос тановка граничных задач, анализ и обобщение результатов;

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

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

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

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

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

Структура и объем работы. Диссертационная работа состоит из введения, шести глав, заключения, списка литературы из наименований и приложений. В приложении приведены 6 актов внедрения и использования. Материал изложен на 345 страницах, содержит рисунков и 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, в качестве примера, представлена зависимость амплитудных нормированных значений модуля тока от час тоты для относительного радиуса R 0 = 5,516 (R0 = a/h) для ЦТС 19. Сплошной линией показаны теоретические значения, пунктирной – kh 2, – частота, к – плотность, с – экспериментальные, = c нормирующий коэффициент.

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

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

Следующей областью со слабой зависимостью резонансной час тоты от радиуса является область толщинного резонанса. Были исследованы формы колебаний в районе 3 на 10-ой и 11-ой гармониках на участках, на которых имеются плато. Во всех исследованных случаях поршневое движение отсутствует. Проведен анализ коэффициента электромеханической связи (КЭМС) по формуле Мэзона f a2 f r kd =, (1) f a где fr и fa – частоты резонанса и антирезонанса. Для области 2,96 3, определенный по указанной формуле КЭМС оказывается незначительным:

k d 0,03. Отсюда следует, что эта облас ть не соответствует толщинному резонансу, поскольку толщинный коэффициент связи для ЦТС-19 k t = 0,44.

uz ur а) R0 = 2,4;

б) R 0 = 3,114;

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

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

Полученное теоретическое распределение механических и электрических полей в объеме плиты позволило провести анализ составляющих плотности внутренней энергии в пьезокерамическом диске и получить значения энергетического КЭМС k эн wэу dV k эн = V 0, w у dV wэ dV V V Проведено сравнение результатов расчетов k эн, kd и экспериментальных значений k экс. В большинстве исследованных случаев результаты достаточно близки, особенно при прос тых формах колебаний, что подтверждает обоснованность использования результатов расчетов k d по приведенной методике для практических целей.

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

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

Рассмотрена задача моделирования в ANSYS динамического режима работы составного пьезопреобразователя, нагруженного, в общем случае, на твердую деформируемую среду. В такой модели имеется одно или несколько пьезоэлектрических тел pn, n = 1,2… и упругих сред ek, k =1,2,… Предполагается, что волновые процессы в средах p = U n en и e = Uk ek описываются классическими линейными теориями пьезоэлектричес тва и упругости соответственно. Волновые поля в различных средах согласуются по граничным контактным условиям. Для расчетов в осесимметричной постановке вводилась цилиндрическая система координат r,, z. Для реализации пьезоэлектрических сред был выбран плоский четырехугольный элемент с опцией осесимметричнос ти.

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

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

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

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

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

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

Расхождение между экспериментальными и теоретическими значениями, по-видимому, объясняется невозможностью дос тичь полного соответс твия между принятыми в расчете константами (средними для пьезоматериала) и их действительными значениями в экспериментальной Таблица Аср (%) Кср (%) ср (%) R 1,032 3,7 3,4 0, 2,563 2,7 2,5 0, 5,156 2.3 3,2 1, 10,450 1,7 2,0 0, партии. Это заключение подтверждается тем, что расхождение между двумя теоретическими методами меньше, чем между каждым из методов и экспериментальными значениями.

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

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

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

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

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

По этой же модели были исследованы высокотемпературные преобразователи (ВТ П). Кроме использования специальных пьезокерамических материалов с повышенной температурой Кюри, такие преобразователи имеют и некоторые конс труктивные особенности, в том числе и наличие керамических протекторов и технологических крышек, Ф 0, 0, 0, 184 0, 253 50 100 150 200 250 f, кГц Рис. 3. Рис. 4.

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

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

Погрешнос ть расчетов для резонансных частот не более 16%.

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

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

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

Ф 0, 0,4 0 200 f, кГц Рис. 7 Рис. протектор, инертную массу и демпфер. Рассчитывается реакция преобразователя на принимаемый сигнал с целью анализа характеристик источника и оптимизации преобразователя для данного типа источника.

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

Источник излучения представляет собой гармоническую сосредоточенную силу, заглубленную в изделие и действующую в направлении оси z (источник первого типа) или в перпендикулярном направлении – в направлении оси r (источник второго типа).

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

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

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

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

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

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

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

Исследовано влияние геометрии активных и пассивных элементов конструкции на АЧХ.

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

Упругая полуполоса (волновод) занимает облас ть |x| 1, 0 y.

Боковые грани полуполосы свободны от усилий. На торце y = рассматривается два типа граничных условий: силовые и кинематические.

Метод решения задач показан на примере изгибных колебаний. Решение для полуполосы строится в виде рядов по однородным решениям u ( x, y ) = Cnu n ( x) ei y, v( x, y) = Cn vn ( x )e i y, (2) n n n =1 n = где Cn – неизвестные коэффициенты разложения, определяемые из граничных условий, n – корни дисперсионного уравнения для изгибных или продольных волн Релея-Лэмба. При этом корни, по которым ведется суммирование, выбираются в соответствии с принципом энергетического излучения.

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

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

Поведение этой особенности исследовано О.К. Аксентян, Я.С. Уфляндом и другими и имеет вид {, r } ~ r, r 0, (3) находится из соотношения ( - коэффициент Пуассона) 2 (3 - 4)sin 0,5 = (1 - ) – (1 - 2), 0 0,5. (4) Первый подход основан на задании специального вида асимптотики постоянных в рядах (2). Используя асимптотический вид однородных решений и заданную асимптотику постоянных, получаем, что напряжения в угловой точке содержат бесконечные ряды следующего вида (1 x ) n s e in ny, j = 1, 2, 3,... (5) j n= N где =-(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 и зависимость погрешности в угловых точках от порядка системы на разных частотах. Погрешнос ть снижается с 0, уменьшением нормированной частоты = 2h/c2 (c2 = (/) ) и ростом порядка редуцированной системы. Показано, что при исследованных формах заданных на торце напряжений амплитуды первых мод, соответс твующих действительным корням, с ростом частоты убывают, и более эффективно возбуждаются высшие моды. Следовательно, для повышения эффективнос ти волноводной приемной сис темы при высоких частотах нужно ориентироваться на прием высших мод. Показано размывание формы заданной на торце нагрузки с удалением от торца на низких и высоких частотах. Однако тип волновых процессов (изгибный или продольный) дает информацию о волновом процессе в объекте контроля и, в конечном итоге, об источнике излучения. Изучено волновое поле и его модовый состав при кинематическом нагружении. Показано, что как для сдвигового, так и продольного напряжения в угловых точках имеется характерная особенность.

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

Внешние границы считались свободными.

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

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

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

Волновод занимает облас ть |х| 1, 0 y, полуплоскость – область y 0. На поверхности y = 0 задана набегающая на волновод волна Релея.

В угловых точках области сцепления полуполосы и полуплоскости (±1,0) имеется особеннос ть напряженного состояния. В полярных координатах с центром в угловой точке напряжения имеют вид { } (i ) = i ), r) ~ r 0 1, r ( (i (6) i = 1, 2 соответственно для полуполосы и полуплоскости. Показатель особенности 0 находится из уравнения, полученного О.К. Аксентян.

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

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

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

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

Исследована сходимость решений при увеличении порядка сис темы, оценены погрешности по выполнению условий сопряжения перемещений в области контакта. Предложенная модель позволяет оценить эффективность возбуждения колебаний в волноводе для разных типов волн и частот. При приеме волны Релея на низкой час тоте = 0,5, когда в волноводе существует одна распространяющаяся мода, отношение амплитуды продольного перемещения в волноводе к амплитуде (m) (m) соответс твующего перемещения в падающей волне Релея v /v R = 0,42.

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

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

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

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

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

Методика отличается от аналогов возможностью прямого определения АЧХ чувствительнос ти по свободному полю;

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

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

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

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

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

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

Исследование приемников продольных колебаний основано на применении упругих волноводов. Одна из конс трукций преобразователя для приема продольных колебаний на основе ХГС-2 защищена патентом.

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

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

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

Полученные результаты показывают, что разработанные нами преобразователи не уступают по основным параметрам приемникам известных зарубежных фирм.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

СПИСОК ОСНОВНЫХ ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ 1. Шихман В. М., Трипалин А. С., Бешенцев В. Д. и др. Исследование медленного роста трещин методом акустической эмиссии // Изв.

СКНЦ ВШ, Естеств. науки – 1982. № 2. - С. 3-6.

2. Дорошенко В. А., Трипалин А.С., Шихман В. М. Исследование сигналов акустической эмиссии в процессе переполяризации сегнетоэлектриков // Изв. СКНЦ ВШ. Ес теств. науки. - 1980. № 3. С. 8-11.

3. Пельц С.П., Трипалин А.С. Шихман В.М. Расчет акус тического поля в упругом волноводе конечной длины, нагруженном на торцах // Изв.

Сев.-Кавк. научн. центра Высш. школы. Естеств. науки. 1984. № 3. С. 7-10.

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

Сев.-Кавк. научн. центра Высш. школы. Естеств. науки. – 1985. № 4.

- С. 5-8.

5. Пельц С. П., Шихман В.М. О сходимости метода однородных решений в динамической смешанной задаче для полуполосы // Докл.

АН СССР. - 1987. Т. 295. № 4. - С. 821-824.

6. Пельц С. П., Шихман В.М. Распространение волн в крестообразном соединении бесконечных упругих полос // ПММ. - 1987. - Т. 51. Вып.

1. - С. 54-59.

7. Пельц С. П., Шихман В.М. Рассеяние волны Релея на упругой полуполосе, сцепленной на торце с упругой полуплоскостью // Д АН СССР. - 1987. Т. 292. № 2. - С. 299-303.

8. Жиров В.Е., Захарова С.В., Трипалин А.С., Шихман В. М. Расчет активных элементов преобразователя для акустико-эмиссионного контроля атомных реакторов // Изв. СКНЦ ВШ. Ес тественные науки.

- 1990. № 2. - С. 76-81.

9. Захарова С. В., Шихман В.М. Вынужденные изгибные колебания упругой полуполосы при смешанных граничных условиях. // ПММ. 1996. № 2.- С. 340-345.

10. Наседкин А. В., Шихман В.М., Захарова С.В., Иванилов И. В.

Применение конечно-элементных методов для расчета приемных систем акустико-эмиссионного контроля // Дефектоскопия. 2006. № 2. - С. 16 – 27.

11. Шихман В. М. Взаимодействие поверхностной волны с упругим препятствием // Экологический вестник научных центров Черноморского экономического сотрудничес тва. 2008. № 3. – С. 53 57.

12. Шихман В. М. Распространение волн Релея-Лэмба в Т-образном соединении упругих волноводов // Известия вузов Сев.-Кавк.

региона. Естес тв. науки. 2009. № 1. – С. 35 - 37.

13. Трипалин А.С., Шихман В.М. Ряд пьезоэлектрических преобразователей для приема сигналов акус тической эмиссии // Автоматическая сварка, 1984. № 5. - С. 33 - 37.

14. Трипалин А. С., Шихман В. М., Коваленко В. И. и др. Акустическая эмиссия при малоцикловых испытаниях сварных тавровых элементов // Диагностика и прогнозирование разрушения сварных конструкций. – 1985. Вып.1. - С. 89-93.

15. Трипалин А.С., Шихман В.М., Карпенко С.Б. Расчет характерис тик преобразователей из текстурных пьезоматериалов с волноводами // Техническая диагнос тика и неразрушающий контроль. Киев. 1989. № 3. - C. 52-56.

16. Шихман В.М., Гринева Л.Д. Преобразователи акустической эмиссии на основе современных пьезоматериалов // Техн. диагн. и неразр.

контроль. Киев. - 1994. № 1. - С. 34-38.

17. Трипалин А.С. Кузьмин Г. А., Елецкий С.А., Коломиец и. В., Шихман В.М. Измерение параметров акустической эмиссии при гидравлическом нагружении корпуса реактора ВВЭР- Нововоронежской АЭС // Энергомашиностроение. - 1982. № 3. - С.

33-36.

18. Трипалин А.С., Шихман В.М., Жиров В.Е., Захарова С.В. и др.

Расчет пьезокерамических преобразователей методами трехмерной теории электроупругости // Акус тическая эмиссия материалов и конструкций. Ч.1. Ростов-на-Дону.: Изд.-во РГУ, 1989. - С. 95 - 99.

19. Шихман В.М. Разработка преобразователей акустической эмиссии в НИИМ и ПМ РГУ (основные направления и результаты) // Современные проблемы мех. сплошной среды. Труды 6-ой Междун.

конф. Ростов-на-Дону. 12-14 июня 2000 г. Ростов-на-Дону: Изд-во СКНЦ ВШ. 2001. Т.1. - С. 243-247.

20. Захарова С.В., Шихман В.М. Колебания толстых пьезокерамических дисков на высокочастотных гармониках. Современные проблемы мех. сплошной среды // Труды 4-ой Международной конференции.

Ростов-на-Дону, 1998. Т. 1. - С. 160-164.

21. Иванилов И. В., Наседкин А. В., Шихман В.М. Расчет АЧХ преобразователей акустической эмиссии в осесимметричной постановке методом конечных элементов // Сборник докладов международной научно-практической конф. “Фундаментальные проблемы пьезоэлектрического приборостроения”. (Пьезотехника 2002). 17-21 сент. 2002 г. Тверь, 2002 г. - С. 164-170.

22. Богма К.К., Шихман В. М. Преобразователи акус тической эмиссии.

Инженерный анализ и исследования // Современные проблемы механики сплошной среды. Тр. VIII Межд. конф., г.Ростов-на-Дону, 14-18 октября 2002 г. Т. 1. / Ростов-на-Дону: изд-во “Новая книга”, 2003. С. 41 - 45. Труды 8-ой Международной конференции. г. Ростов на-Дону., г. Ростов-на-Дону. С. 41-45.

23. Шихман В.М., Иванилов И. В. Моделирование приема упругих колебаний акустико-эмиссионным преобразователем методом конечных элементов // Труды III Всероссийской конференции по теории упругости. г. Азов. 13-17 октября 2003. - С. 406-409.

24. Наседкин А. В., Шихман В. М., Иванилов И. В. Конечно-элементное моделирование распространения и приема акустических сигналов о т заглубленного источника в полубесконечной среде // Современные проблемы механики сплошной среды. Тр. VIII Межд. конф., г.Ростов-на-Дону. Т. 1. Ростов-на-Дону: изд-во “Новая книга”, 2003.

С. 116-120.

25. Наседкин А. В., Шихман В.М., Иванилов И. В. Совершенствование приемных устройс тв и методик акустико-эмиссионного контроля объектов железнодорожного и авиационного транспорта на основе теоретического и экспериментального моделирования волновых процессов // Техника, технология и перспективные материалы.

Межвуз. сб. научн. тр. – Москва, МГИУ. 2004. С. 251 - 256.

26. Захарова С. В., Наседкин А. В., Шихман В. М. Зависимость амплитудно-частотной характеристики преобразователей от граничных условий на торце // Современные проблемы механики сплошной среды. Труды Х Межд. конф., г.Ростов-на-Дону, 5- декабря 2006 г. Т. 1. Ростов-на-Дону: «Новая книга», 2007. - С. 128 132.

27. Захарова С.В., Наседкин А. В., Шихман В. М. Исследование многослойных высокочувствительных акусто-эмиссионных датчиков // Современные проблемы механики сплошной среды. Тр. XI Межд.

конф., г. Ростов-на-Дону, 26-29 ноября 2007 г. Т. 2. Ростов-на-Дону:

изд-во "ЦВВР", 2008. С. 88-92.

28. Шихман В. М., Наседкин А. В., Захарова С.В. Применение конечно элементных методов для расчета многослойных пьезопреобразователей в режиме приема акустических колебаний // Материалы 16-той международной конференции «Современные методы и средства неразрушающего контроля и технической диагностики». Ялта, 1 - 5 октября, 2008. - С. 66 - 69.

29. Шихман В.М. Перспективные направления в разработке преобразователей акустической эмиссии // ХVI Российская н.-т.

конференция “Неразрушающий контроль и диагнос тика”. – Москва.

28 июня - 2 июля 1999 г. Тез. докл. - Т.2. С. 111.

30. Шихман В.М. Расчет приемных ус тройс тв для акустических методов неразрушающего контроля // ХХ Уральская региональная конф.

«Контроль технологий, изделий и окружающей среды физическими методами». 15-16 мая 2001 г. Екатеринбург. Тез. докл. - С. 60-61.

31. Шихман В.М. Разработка нового поколения пьезоэлектрических преобразователей для акустико-эмиссионной диагностики // ХХ Уральская региональная конференция «Контроль технологий, изделий и окружающей среды физическими методами». 15-16 мая 2001 г. Екатеринбург. Тез. докл. - С. 58-59.

32. Патент № 1784095 (СССР). Козинкина А. И., Трипалин А.С. и Шихман В.М. Пьезоэлектрический преобразователь для приема сигналов акустической эмиссии // Опубл. в Б.И. № 47. 1992.

33. А.с. 918284 (СССР). Способ поляризации пьезокерамических элементов / Трипалин А.С., Шихман В.М.. - Опубл. в Б.И. № 13.1982.

34. А.с. № 1167497 (СССР). Пьезоэлектрический преобразователь (его варианты) / Трипалин А.С., Шихман В.М.. - Опубл. в Б.И. № 26.1985.

35. А.с. № 1260849 (СССР). Ультразвуковой преобразователь для калибровки систем акустико-эмиссионного контроля. / Кузьмин Г.А., Крейнер Б.Я., Шихман В.М. Опубл. в БИ. № 4, 1985.

36. А.с. № 1320734 (СССР). Способ акустикоэмиссионного контроля изделий / Мельцер Я.Е., Трипалин А.С., Шихман В. М., Колесников В. И. - Опубл. в Б.И. № 24, 1987.

37. 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.

38. 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.

39. 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.

40. 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. 41. 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.

42. Shikhman V.M., Grineva L.D. High-efficiency Transducers for Acoustic Diagnostics on the Bases of Novel Anisotropic Piezoceramic // Materials.

th Abstracts of 9 International Symposium on Application of Ferroelectrics.

- Penn. State Scantion conference Center, Pennsylvania.- USA. Aug.7-10.

1994. - P. 39.

43. 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.

44. 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. 45. 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.

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.