Способ измерения полного набора модулей пьезоэлектрической керамики на одном образце Российский патент 2018 года по МПК G01N29/12 

Описание патента на изобретение RU2663271C1

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

Известны способы определения полного набора модулей (констант) для пьезкерамической керамики:

ОСТ 11.0444-87. Материалы пьезокерамические. Технические условия. - Москва, 1987.- 141 с.,

IRE standards on piezoelectric crystals, measurements of piezoelectric ceramic. – Proc. IRE, 1961. V.49 - P. 1162.,

European standard EN 50324-2:2002 – CENELEC European Committee for Electrotechnical Standardization, Piezoelectric components, Part 2: Metods of measurent – Low poper, Accessed on:2002.

Недостатками вышеуказанных способов являются:

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

- сложно определить (измерить) на измерительном стенде частоты «антирезонанса» для высокочастотных толщинных колебаний из-за наличия вблизи антирезонансных частот множества других мод колебаний [1] ;

- для поляризации стержней с межэлектродным расстоянием порядка 10-20 мм часто требуются электрические напряжения до 100 - 200 кВ;

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

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

Наиболее близким является способ «Способ определения упругих податливостей s11E,s12E,s13E,s33E и пьезоэлектрических модулей d31,d33 на одном образце в виде диска» [патент RU № 2629927, МПК G01N 29/12, опубл. 04.09.2017], где на одном образце в виде диска, по которому возбуждают колебания и измеряют частоты резонансов и первого антирезонанса образца, возбуждают радиальные и толщинные колебания, дополнительно возбуждают и измеряют частоту первого и второго обертонов радиальных колебаний и вторую частоту резонансных колебаний образца по толщине, и по этим частотам определяют модули исследуемого пьезокерамического материала.

Недостатками данного способа являются:

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

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

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

Вначале для определения первых семи модулей s11 , s12 , s13 , s33 , d31 , d33 , ε33T, кроме сдвиговых, в отличие от известных методов для повышения точности определения констант ПКМ измерения резонансных частот (трех радиальных с учетом толщины диска и двух толщинных мод) и первого низкочастотного антирезонанса проводятся только на одном образце в виде поляризованного по толщине кольца с электродами на торцах. Затем для определения сдвиговых модулей s44 , d15 , ε11T измеряют первые два толщинных резонансов на том же кольце, но с новыми электродами на боковых цилиндрических поверхностях. Старые электроды на торцах удаляются. Использование для определения полного набора совместимых констант только одного образца вместо четырех повышает точность измерения упругих модулей, т.к. отсутствуют погрешности, вызванные разбросом свойств образцов различных форм в процессах спекания и поляризации.

Способ измерения полного набора констант на одном образце осуществляется следующим образом.

Из исследуемого материала изготовляют кольцо (высотой 2h, внутреннего радиуса b и внешнего a, плотностью ρ) с соотношением сторон: диаметр/толщина не менее 40:1, а отношение a/b не менее 5. Сначала наносят электроды на торцевые поверхности кольца и поляризуют по толщине. Возбуждают низкочастотные радиальные колебания в исследуемом образце при подаче электрического сигнала на электроды. Измеряют три первые резонансные частоты fr (основной резонанс и два его обертона), а упругие податливости s11,s12,s13 (здесь и далее верхний индекс E у податливостей опущен) определяют из решения частотного уравнения радиальных колебаний кольца с учетом его относительной толщины.

Для рассматриваемого кольца решение электроупругой задачи [1] ищем в виде разложения по малому параметру - относительной частоте Ω2 = ρ(ωh)2/c44

α2 = tΩ2 + γΩ4 + ηΩ6… U = U0 + Ω2U2 + Ω4 U4 +… (1)

где α – безразмерное волновое число ; t, γ, η - неизвестные постоянные, зависящие от модулей пьезокерамики; Un – вектор, составляющие которого U, W- механические перемещения, φ – электрический потенциал; сij – упругие модули; ω – круговая частота.

Опуская громоздкие выкладки и ограничиваясь членами с Ω4 для определения волнового числа α2 и Ω2 для смещений U в (1), приведем конечный результат рассматриваемой краевой задачи с нулевыми (свободными) механическими граничными условиями на торцах и электрическими условиями короткого замыкания :

t = c44 /cP11; γ = t2(c13 /c33)2 /3; c13 /c33 = -s13/(s11+s12);

cP11 = c11 – c213 /c33 = 1/(s11(1– ν2)); ν = –s12 /s11;

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

kp2 – планарный коэффициент связи;

α20 = tΩ2 – приближенное волновое число нулевого порядка;

α22 = tΩ2 + γΩ4 - приближенное волновое число второго порядка.

Частотное уравнение, соответствующее низкочастотным колебаниям, для колебаний кольца с учетом его относительной толщины ε=h/a имеет вид:

a11 a22 – a12 a21 =0 (2)

где a11 = F1(x1) [1+ α02 c13/c33(1/3+(c13/6c33))–1/3 (α0 c13/c33)2J0(x1)];

a12 = F2(x1) [1+ α02 c13/c33(1/3+(c13/6c33))–1/3 (α0 c13/c33)2N0(x1];

a21 = F1(x2) [1+ α02 c13/c33(1/3+(6c13/c33))–1/3 (α0 c13/c33)2J0(x2)];

a22 = F2(x2) [1+ α02 c13/c33(1/3+(c13/6c33))–1/3 (α0 c13/c33)2N0(x x2)];

x1= α/ε ; x2= x1 b/a;

F1(x)= xJ0(x)–(1–ν)J1(x)]; F2(x)= xN0(x)–(1– ν)N1(x); {x= x1, x2}.

В уравнение (2) ν – коэффициент Пуассона; J0(x), J1(x) – функции Бесселя первого рода; N0(x), N1(x) – функции Бесселя второго рода порядка нулевого и первого.

Частотное уравнение (2) для тонких ПКМ в форме кольца зависит только от упругих податливостей s11 , s12 , s13 и относительной толщины ε через α2. а также в неявном виде через соотношение Ω(α) входит толщинная поправка ε=h/b к определению резонансных частот

Ω2 = α2 {1 – (ε α c13/ c33)2 /3} .

И именно этим соотношением приведенное выше частотное уравнение (2) для радиальных колебаний кольца конечной толщины отличается от известного частотного уравнения для радиальных колебаний с «нулевой толщиной» [3].

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

Если измерить три резонансные частоты, получим три уравнения (2) относительно трёх неизвестных. Учитывая соотношения

ν = – s12 / s11 ; с13 / с33 = – s13 /( s11 + s12),

из данного трехчастотного метода можно численно определить эти три неизвестные величины s11, s12 , s13.

Диэлектрическую константу ε33T свободного образца рассчитаем из измеренной емкости C для диска на частоте 1 кГц [2]

ε33T = 2Ch/(π(a2 - b2)).

Для определения величин пьзомодуля d31 или коэффициента связи kp2 необходимо воспользоваться другим типом однородных электрических условий – равенством нулю тока на электродах. Этот тип граничных условий соответствует условиям антирезонанса. Ограничиваясь приближением нулевого порядка в (1), получим из [3] условие для первого низкочастотного антирезонанса:

1- kp2 + kp2(1+ν) [Δ1/Δ (J1(x1) – b/a J1(x2)) + Δ2/Δ (N1(x1) – b/a N1(x2))]/(1 – b2/a2) = 0

где x1 = α0/ε ; x2 = x1 b/a; Δ= F1(x1) F2(x2) – F2(x1) F1(x2) ; (3)

Δ1= F2(x2) – b/a F2(x1); Δ2= b/a F1(x1) – F1(x2);

При b → 0 уравнение (3) вырождается в известное уравнение антирезонансных частот диска нулевой толщины [1,2].

Пьезомодуль d31 определяем по формуле d312 = kp2(s11+s12) ε33T/2.

Упругую константу с33D и коэффициент связи kt определяют из первых двух измеренных резонансных частот ft толщинных колебаний согласно соотношениям [2]

1 – kt2 tg(β)/β =0; β2 = (2πhft)2ρ/c33D . (4)

Для этого повышаем частоту возбуждения колебаний исследуемого ПКМ. Толщинный резонанс легко определяется, т.к. его амплитуда значительно превышает амплитуды сигналов других мод колебаний. Измерение второго толщинного резонанса повышает точность измерения параметров с33D и kt , т.к. в отличии от известных ранее методов определения модулей керамики в рассматриваемом способе нет измерений высокочастотных антирезонансов с электродами на торцах: применение толщинных мод колебаний пластин в области антирезонансных частот для определения констант пьезокерамики сопряжено с массой сложностей: технически трудно решить какой минимум соответствует какому максимуму для высокочастотных колебаний - и поэтому их использование нежелательно [1].

Имея два уравнения (2) с двумя неизвестными численными методами можно определить c33D и kt2. Далее используя соотношения

kt2/(1– kt2)= (e33/c33)2 c33D/ (ε33T(1–kp2)); d33 = e33/c33– 2d31 с13/c33

легко рассчитать пьезомодуль d33 .

По определенным модулям исследуемой ПКМ s11 , s12 , s13 , с33D можно рассчитать упругую податливость s33.

с33=(1– kt2) c33D ; s33=1/ с33 + 2s13/(s11+s12).

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

Решение данной краевой антисимметричной задачи (1-3) состоит из суммы двух решений при нулевых граничных условий на неэлектродированных торцах при z = ±h :

- однородного решения [4];

- частного решения.

Однородное решение в явном виде не участвует в расчете проводимости Y, что следует из теоремы Гаусса для пьезокерамики

div D =0, (5)

где D – вектор электрической индукции.

Частное решение не составляет труда построить. Пусть компоненты вектора электрической индукции частного решения равны Dz0 = 0 и Dr0 11SA/r , где величина A не зависит от координаты r и определяется из электрических граничных условий при r=a,b. Данный выбор компонент вектора электрической индукции позволяет автоматически удовлетворить уравнению (5) . Далее потребуем, чтобы механические смещения соответственно были равны W = 0 и U = Bh sin(Ως)/r. Тогда механические составляющие частного решения для антисимметричных колебаний, удовлетворяющие нулевым граничным условиям на неэлектродированных торцах кольца имеют вид:

Trz0= Ae15/(r ln(a/b)) (1– cos(Ως)/ cos(Ω));

Trr0 =Ah(c11 – c12) e15 sin(Ως)/( r2 ln(a/b) c44 Ω cos(Ω)); Tzz0 =0;

U0 = - h A e15 sin(Ως) )/( r ln(a/b) c44 Ω cos(Ω));

Используя закон Ома I=2V*Y, формулу Dr = (e15 U,z – ε11 φ,r), где запятая обозначает дифференцирование по координате следующей за запятой, и проинтегрировав по dS ( где S - площадь электродов на боковой поверхности кольца), получим проводимость частного решения Y0 рассматриваемого кольца

Y0= –iω Ck (1– k152 + k152 tg(Ω)/ Ω) (6)

где величина Ck для низших частот равна 2h ε11T/ ln(a/b).

Из выражения (6) при низких частотах Ω→0 получим формулу Y= –iωCk , из которой легко рассчитать ε11T. Проводимость Y0 → max при Ω = π/2, что соответствует резонансной частоте

fY = 923760 Гц для PZT4. Антирезонанс fA из (18) получим при условии Y0 =0, или

1- k152 + k152 tg(Ω)/ Ω =0

где k152 = с44 d152 / ε11T

Для проверки корректности методики используется численное решение краевых задач методом конечных элементов. В качестве программы, осуществляющей реализацию метода конечных элементов, используется программа ANSYS [5].

Алгоритм проверки следующий:

1. Для кольца из PZT4 с электродами на торцах (a=15 мм, b=3 мм, h= 1 мм) из рассчитанных с помощью ANSYS первых трех низкочастотных резонансных частот (f0=66,719, f1 = 174,460, f2 = 303,480 КГц) были численно найдены из (2) упругие модули гибкости s11 = 0,1229⋅10-10; s12 = –0,04049⋅10-10; s13 = –0,052844⋅10-10. Ошибка для модулей гибкости не превышает 1%. Используя рассчитанные на ANSYS первые два толщинных резонанса ft1=2,02 МГц и ft2 = 6,728 МГц , из (4) получим упругую константу c33D = 15,427⋅1010 и пьзомодуль

d33 = 2,8584⋅10-10.

Пьезомодуль d31 определим по известной методике [2,3] из первого антирезонанса. Модуль равен –1,23⋅10-10 для частоты 73959 Гц.

2. Для исследуемого кольца, но уже с электродами на боковых поверхностях, из рассчитанных с помощью ANSYS толщинных критических частот fY = 939 кГц и fА =1190 кГц, по формулам (6) можно рассчитать с44 и k152 . Пьезомодуль d15 легко рассчитать, используя формулу

k152 = с44 d152 / ε11T.

3. Диэлектрические проницаемости ε33T и ε11T рассчитываем из измерений емкостей на частоте 1 КГц согласно [2] для колец с электродами на торцах и на боковых поверхностях соответственно.

Приведем все модули для PZT4, определенные новым методом на одном образце с помощью программы ANSYS:

s11 = 0,123 (0,123); s12 = –0,0405 (–0,0405); s13 = –0,0528 (–0,0531); s33 = 0,153 (0,155);

s44 =0,38 (0,39); d31 =–1,23 (–1,23); d33 =2,8584 (2,89); d15 =5,001 (4,96);

ε33T =114.95 (115,05); ε11T = 130 (130,54).

Здесь в скобках указаны табличные справочные данные для PZT4 [6]. Порядковый коэффициент для модулей соответственно 1010 или 10-10 опущен. Ошибка не превышает 1 %.

Похожие патенты RU2663271C1

название год авторы номер документа
Способ определения упругих податливостей s11Е, s12Е, s13Е, s33Е и пьезоэлектрических модулей d31,d33 на одном образце в виде диска 2016
  • Мадорский Виктор Вениаминович
RU2629927C1
ПЬЕЗОЭЛЕКТРИЧЕСКИЙ КЕРАМИЧЕСКИЙ МАТЕРИАЛ 2014
  • Резниченко Лариса Андреевна
  • Таланов Михаил Валерьевич
  • Вербенко Илья Александрович
  • Шилкина Лидия Александровна
RU2561439C2
Способ определения пьезомодулей 1991
  • Земляков Виктор Леонидович
SU1800406A1
Низкочастотный пьезоэлектрический керамический материал на основе ниобата натрия 2020
  • Резниченко Лариса Андреевна
  • Андрюшин Константин Петрович
  • Глазунова Екатерина Викторовна
  • Дудкина Светлана Ивановна
  • Андрюшина Инна Николаевна
  • Вербенко Илья Александрович
RU2751323C1
Пьезоэлектрический керамический материал на основе метаниобата лития 2019
  • Резниченко Лариса Андреевна
  • Дудкина Светлана Ивановна
  • Разумовская Ольга Николаевна
  • Андрюшин Константин Петрович
  • Андрюшина Инна Николаевна
  • Вербенко Илья Александрович
RU2712083C1
ПЬЕЗОЭЛЕКТРИЧЕСКИЙ КЕРАМИЧЕСКИЙ МАТЕРИАЛ 2013
  • Резниченко Лариса Андреевна
  • Разумовская Ольга Николаевна
  • Павленко Анатолий Владимирович
  • Вербенко Илья Александрович
  • Дудкина Светлана Ивановна
  • Самойлова Виктория Игоревна
RU2542008C1
Высокочастотный пьезоэлектрический керамический материал на основе титаната-цирконата свинца 2021
  • Андрюшин Константин Петрович
  • Андрюшина Инна Николаевна
  • Глазунова Екатерина Викторовна
  • Дудкина Светлана Ивановна
  • Мойса Максим Олегович
  • Вербенко Илья Александрович
  • Резниченко Лариса Андреевна
RU2764404C1
ПЬЕЗОЭЛЕКТРИЧЕСКИЙ КЕРАМИЧЕСКИЙ МАТЕРИАЛ 2010
  • Резниченко Лариса Андреевна
  • Разумовская Ольга Николаевна
  • Андрюшин Константин Петрович
  • Вербенко Илья Александрович
  • Павелко Алексей Александрович
  • Таланов Михаил Валерьевич
  • Павленко Анатолий Владимирович
RU2440955C2
ПЬЕЗОЭЛЕКТРИЧЕСКИЙ КЕРАМИЧЕСКИЙ МАТЕРИАЛ 2013
  • Резниченко Лариса Андреевна
  • Таланов Михаил Валерьевич
  • Шилкина Лидия Александровна
  • Вербенко Илья Александрович
RU2551156C1
СПОСОБ ИЗМЕРЕНИЯ МАГНИТОСТРИКЦИИ ТОНКИХ ПЛЕНОК 2020
  • Филиппов Дмитрий Александрович
  • Лалетин Владимир Михайлович
  • Поддубная Наталья Никитична
RU2739161C1

Реферат патента 2018 года Способ измерения полного набора модулей пьезоэлектрической керамики на одном образце

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

Формула изобретения RU 2 663 271 C1

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

Документы, цитированные в отчете о поиске Патент 2018 года RU2663271C1

Способ определения упругих податливостей s11Е, s12Е, s13Е, s33Е и пьезоэлектрических модулей d31,d33 на одном образце в виде диска 2016
  • Мадорский Виктор Вениаминович
RU2629927C1
СВЕРЛИЛЬНЫЙ ПАТРОН 1927
  • Омельяненко Т.И.
SU6210A1
Способ измерения констант упругости в пьезоэлектрической керамике 1983
  • Акопьян Владимир Акопович
SU1193573A1
US 3954015 A, 04.05.1976
JP 2005201816 A, 28.07.2005
CN 106353404 A, 25.01.2017.

RU 2 663 271 C1

Авторы

Мадорский Виктор Вениаминович

Даты

2018-08-03Публикация

2017-10-12Подача