Упрощенное определение параметров пьезоматериалов на образцах элементов в форме диска
Аннотация
Получены приближенные формулы для определения планарного коэффициента электромеханической связи и пьезомодуля на образцах пьезоэлементов в форме диска, в которых отсутствуют коэффициент Пуассона и наименьший корень частотного уравнения. Это упрощает процедуру измерений, поскольку не требуется измерение частоты первого обертона.
Ключевые слова: пьезоматериал, диск, планарный коэффициент электромеханической связи, пьезомодуль, коэффициент Пуассона
В качестве активного элемента пьезоэлектрических датчиков часто используют пьезокерамические элементы (ПКЭ) в форме диска. Такие датчики находят широкое применение для контроля состояния технически сложных объектов.
Поскольку качество ПКЭ и датчика в целом напрямую зависит от параметров пьезоматериала, в частности, от его планарного коэффициента электромеханической связи и пьезомодуля , определение указанных параметров является часто выполняемой процедурой.
Однако, эта процедура в стандартом варианте достаточно сложная, поскольку помимо измерения частот резонанса fp и антирезонанса fа на основной моде колебаний необходимо дополнительно проводить измерения частоты первого обертона fp1 (отношение частоты этого обертона к частоте резонанса fp1/fp обычно лежит в интервале 2,54 – 2,57). Затем рассчитывать коэффициент Пуассона σ, а также находить наименьший корень частотного уравнения η, содержащий функции Бесселя нулевого и первого порядков по таблицам [1]. Формула для расчетов получается громоздкой и не удобной для инженерных расчетов.
Рассмотрим возможность определения указанных выше параметров по приближенным формулам, в которых отсутствует зависимость от η и σ.
Будем использовать известные исходные соотношения [2,3]:
, (1)
, (2)
, (3)
, (4)
где – компонента упругой податливости, – плотность материала, d – диаметр диска, – диэлектрическая проницаемость, рассчитываемая через емкость ПКЭ на низкой частоте .
С учетом выражений (2) и (3), преобразуем (4) к виду:
. (5)
Вначале определим простую аналитическую зависимость между η и σ. Для этого методом наименьших квадратов выполним аппроксимацию данных для η, заданных в табличной форме.
График зависимости η от σ, построенный по табличным значениям показан на рис. 1 линией 1 и по внешнему виду близок к прямой. Поэтому для аппроксимации использовалась функция вида ηа=aσ+b и была получена следующая аппроксимирующая функция:
. (6)
Результаты расчетов по формуле (6) приведены на рис. 1 линией 2 и практически совпадают с табличными данными для η.
Рис. 1. Зависимость η от σ
Подставляя (6) в формулу (3) для определения планарного коэффициента электромеханической связи, и выполняя несложные преобразования, получим
. (7)
Если представить формулы (3) и (7) в виде , то можно построить, приведенные на рис. 2 зависимости .
Обращают на себя внимание следующие свойства формулы (7). Во-первых, отсутствует зависимость kp от η, во-вторых, видна очень слабая зависимость планарного коэффициента электромеханической связи от коэффициента Пуассона: при изменении σ от 0,24 до 0,5 изменение kp составляет всего несколько процентов.
1 - 2 - , 3 – 1,17
Рис. 2. Зависимость
Последнее обстоятельство позволяет сделать вывод о том, что в практике измерений планарного коэффициента электромеханической связи иногда можно вообще отказаться от измерений частоты первого обертона fp1 и определения σ, а расчет kp вести по приближенной формуле
, (8)
в которой коэффициент 1,17 соответствует среднему значению при изменении σ в формуле (7) от 0,24 до 0,5 (линия 3 на рис. 2). При этом дополнительная методическая погрешность определения kp на краях диапазона изменения σ не будет превышать ±1,5%. В середине же этого диапазона, что чаще всего встречается на практике, дополнительная погрешность составляет доли процента.
Из (8) также следует простая зависимость между планарным коэффициентом электромеханической связи пьезоматериала и эффективным коэффициентом электромеханической связи ke для ПКЭ в форме диска: .
Перейдем теперь к определению пьезомодуля.
Обозначим в (5) выражение в скобках буквой и решим задачу путем графического построения зависимости от σ. Затем, методом наименьших квадратов, выполним аппроксимацию данных.
График зависимости от σ показан на рис. 3 линией 1 и по внешнему виду близок к прямой. Поэтому для аппроксимации использовалась функция вида Fd=aσ+b и была получена следующая аппроксимирующая функция:
. (9)
Результаты расчетов по формуле (9) приведены на рис. 3 линией 2 и практически совпадают с линией 1.
1 – формула (5), 2 – формула (9), 3 – Fd = 1,88.
Рис.3. Зависимость Fd от σ
Поскольку функция меняется незначительно, то ее можно приближенно заменить константой, равной 1,88 (линия 3 на рис. 3). Методическая погрешность определения пьезомодуля при таком приближении не превысит ±2%, а формула (5) будет иметь простой вид:
. (10)
В таблице 1 представлены результаты экспериментальных исследований 3-х образцов ПКЭ из разного пьезоматериала, имеющих одинаковый диаметр d=2 см и разную толщину t.
Таблица 1
Результаты измерений на образцах ПКЭ в форме диска
|
t = 2 мм |
t = 3 мм |
t = 1 мм |
, нФ |
2,954 |
1,61 |
5,783 |
, кГц |
105,105 |
104,542 |
127,142 |
, кГц |
122,690 |
117,951 |
137,672 |
, кГц |
267,763 |
259,600 |
330,429 |
|
0,41 |
0,49 |
0,34 |
|
2,115 |
2,14 |
2,074 |
х10-9, Ф/м |
18,8 |
15,37 |
18,4 |
x10-12, м2 /н |
16,5 |
18,7 |
10,2 |
|
0,55 |
0,48 |
0,41 |
|
0,56 |
0,5 |
0,42 |
x10-12, Кл/Н |
167,4 |
132,4 |
103,6 |
x10-12, Кл/Н |
169,1 |
136,5 |
103,4 |
Таким образом, при определении планарного коэффициента электромеханической связи и пьезомодуля пьезокерамических материалов на образцах элементов в форме диска зачастую можно отказаться от значений коэффициента Пуассона и наименьшего положительного корня частотного уравнения и вести расчет по приближенным формулам. При этом нет необходимости дополнительно измерять частоту первого обертона, лежащую существенно выше частоты основного резонанса ПКЭ.
Кроме этого, поскольку измерения проводятся только в одной резонансной области, приближенные формулы позволяют реализовать оперативный контроль пьезомодуля пьезоматериала в партии образцов ПКЭ в форме диска.
Такая возможность обусловлена линейной зависимостью между квадратом пьезомодуля и интегралом от активной составляющей проводимости (площадью под кривой активной составляющей проводимости) ПКЭ [4]:
. ( )
Константа зависит от упругих свойств материала и геометрии ПКЭ.
На рис. 4 представлены результаты измерений активной составляющей проводимости 3-х ПКЭ, которые позволяют проверить их на идентичность по величине пьезомодуля в относительных единицах через площадь под кривой активной составляющей проводимости:
Рис. 4. Зависимость активной составляющей проводимости от частоты
В таблице 2 приведены результаты сравнения. Величина квадрата пьезомодуля первого ПКЭ взята за 100% и относительно этого образца определяют в процентном отношении значение квадратов пьезомодуля других ПКЭ
Таблица 2
Сравнение относительных значений величины пьезомодуля
для трех образцов ПКЭ
№ |
|
Относ. |
I |
Относ. |
1 |
22052 |
100 % |
10,46 |
100 % |
2 |
15836 |
72 % |
7,37 |
70 % |
3 |
19902 |
90 % |
9,54 |
91 % |
Простая аппаратная реализация оперативного контроля пьезоматериала по величине пьезомодуля в виде приставки к анализатору частотных характеристик может состоять в следующем.
К выходу генератора качающейся частоты (ГКЧ) подключают последовательно соединенные ПКЭ и вспомогательный резистор Rв, величиной значительно меньшей сопротивления ПКЭ на частоте резонанса, как это обычно делается при исследовании частотных характеристик проводимости. Параллельно вспомогательному резистору подключают синхронный детектор, причем, один вход синхронного детектора подключен к вспомогательному резистору, а другой – к выходу ГКЧ. За синхронным детектором следует интегратор.
В зависимости от величины добротности ПКЭ устанавливают время качания частоты, и включают ГКЧ, который начинает формировать на выходе синусоидальный сигнал с линейно меняющейся частотой. Одновременно линейно нарастающее (пилообразное) напряжение, используемое для работы ГКЧ, включает интегратор. На выходе синхронного детектора формируется частотная характеристика активной составляющей проводимости. К моменту завершения сканирования частоты на выходе интегратора появляется сигнал постоянного уровня, соответствующий интегралу от активной составляющей проводимости в заданном диапазоне частот (соответствующий площади под кривой активной составляющей проводимости). Этот сигнал, пропорциональный квадрату пьезомодуля, можно использовать непосредственно для управления процессом разбраковки ПКЭ при их производстве.
Работа выполнена при финансовой поддержке Минобрнауки России в рамках ФЦП «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России на 2007—2013 годы», госконтракт № 16.552.11.7024.
Литература
1.ГОСТ 12370 - 80. Материалы пьезокерамические. Методы испытаний. М.: Издательство стандартов, 1980.
2.ОСТ 11 0444-87. Материалы пьезокерамические. Технические условия. М.: Электростандарт, 1987.
3.Смажевская Е.Г., Фельдман И.Б. Пьезоэлектрическая керамика. М.: Советское радио, 1971, 182с.
4.Земляков В.Л., Земляков В.В Новый подход к измерению пьезо-модуля пьезокерамических материалов в динамическом режиме // Измерительная техника. 2002. № 4. С. 52−55.