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

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

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

ОСТ 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 кВ;

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

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

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

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

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

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

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

Вначале для определения первых семи модулей s₁₁ , s₁₂ , s₁₃ , s₃₃ , d₃₁ , d₃₃ , ε₃₃^T, кроме сдвиговых, в отличие от известных методов для повышения точности определения констант ПКМ измерения резонансных частот (трех радиальных с учетом толщины диска и двух толщинных мод) и первого низкочастотного антирезонанса проводятся только на одном образце в виде поляризованного по толщине кольца с электродами на торцах. Затем для определения сдвиговых модулей s₄₄ , d₁₅ , ε₁₁^T измеряют первые два толщинных резонансов на том же кольце, но с новыми электродами на боковых цилиндрических поверхностях. Старые электроды на торцах удаляются. Использование для определения полного набора совместимых констант только одного образца вместо четырех повышает точность измерения упругих модулей, т.к. отсутствуют погрешности, вызванные разбросом свойств образцов различных форм в процессах спекания и поляризации.

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

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

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

α² = tΩ² + γΩ⁴ + ηΩ⁶… U = U₀ + Ω²U₂ + Ω⁴ U₄ +… (1)

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

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

t = c₄₄ /c^P₁₁; γ = t²(c₁₃ /c₃₃)² /3; c₁₃ /c₃₃ = -s₁₃/(s₁₁+s₁₂);

c^P₁₁ = c₁₁ – c²₁₃ /c₃₃ = 1/(s₁₁(1– ν²)); ν = –s₁₂ /s₁₁;

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

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

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

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

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

a₁₁ a₂₂ – a₁₂ a₂₁ =0 (2)

где a₁₁ = F₁(x₁) [1+ α₀² c₁₃/c₃₃(1/3+(c₁₃/6c₃₃))–1/3 (α₀ c₁₃/c₃₃)²J₀(x₁)];

a₁₂ = F₂(x₁) [1+ α₀² c₁₃/c₃₃(1/3+(c₁₃/6c₃₃))–1/3 (α₀ c₁₃/c₃₃)²N₀(x₁];

a_{21 =} F₁(x₂) [1+ α₀² c₁₃/c₃₃(1/3+(6c₁₃/c₃₃))–1/3 (α₀ c₁₃/c₃₃)²J₀(x₂)];

a₂₂ = F₂(x₂) [1+ α₀² c₁₃/c₃₃(1/3+(c₁₃/6c₃₃))–1/3 (α₀ c₁₃/c₃₃)²N₀(x x₂)];

x₁= α/ε ; x₂= x₁ b/a;

F₁(x)= xJ₀(x)–(1–ν)J₁(x)]; F₂(x)= xN₀(x)–(1– ν)N₁(x); {x= x₁, x₂}.

В уравнение (2) ν – коэффициент Пуассона; J₀(x), J₁(x) – функции Бесселя первого рода; N₀(x), N₁(x) – функции Бесселя второго рода порядка нулевого и первого.

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

Ω² = α² {1 – (ε α c₁₃/ c₃₃)² /3} .

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

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

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

ν = – s₁₂ / s₁₁ ; с₁₃ / с₃₃ = – s₁₃ /( s₁₁ + s₁₂),

из данного трехчастотного метода можно численно определить эти три неизвестные величины s₁₁, s₁₂ , s₁₃.

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

ε₃₃^T = 2Ch/(π(a² - b²)).

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

1- k_p² + k_p²(1+ν) [Δ₁/Δ (J₁(x₁) – b/a J₁(x₂)) + Δ₂/Δ (N₁(x₁) – b/a N₁(x₂))]/(1 – b²/a²) = 0

где x₁ = α₀/ε ; x₂ = x₁ b/a; Δ= F₁(x₁) F₂(x₂) – F₂(x₁) F₁(x₂) ; (3)

Δ₁= F₂(x₂) – b/a F₂(x₁); Δ₂= b/a F₁(x₁) – F₁(x₂);

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

Пьезомодуль d₃₁ определяем по формуле d₃₁² = k_p²(s₁₁+s₁₂) ε₃₃^T/2.

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

1 – k_t² tg(β)/β =0; β² = (2πhf_t)²ρ/c₃₃^D . (4)

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

Имея два уравнения (2) с двумя неизвестными численными методами можно определить c₃₃^D и k_t². Далее используя соотношения

k_t²/(1– k_t²)= (e₃₃/c₃₃)² c₃₃^D/ (ε₃₃^T(1–k_p²)); d₃₃ = e₃₃/c₃₃– 2d₃₁ с₁₃/c₃₃

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

По определенным модулям исследуемой ПКМ s₁₁ , s₁₂ , s₁₃ , с₃₃^D можно рассчитать упругую податливость s₃₃.

с₃₃=(1– k_t²) c₃₃^D ; s₃₃=1/ с₃₃ + 2s₁₃/(s₁₁+s₁₂).

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

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

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

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

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

div D =0, (5)

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

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

T_rz⁰= Ae₁₅/(r ln(a/b)) (1– cos(Ως)/ cos(Ω));

T_rr⁰ =Ah(c₁₁ – c₁₂) e₁₅ sin(Ως)/( r² ln(a/b) c₄₄ Ω cos(Ω)); T_zz⁰ =0;

U⁰ = - h A e₁₅ sin(Ως) )/( r ln(a/b) c₄₄ Ω cos(Ω));

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

Y⁰= –iω C_k (1– k₁₅² + k₁₅² tg(Ω)/ Ω) (6)

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

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

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

1- k₁₅² + k₁₅² tg(Ω)/ Ω =0

где k₁₅² = с₄₄ d₁₅² / ε₁₁^T

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

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

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

d₃₃ = 2,8584⋅10^-10.

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

2. Для исследуемого кольца, но уже с электродами на боковых поверхностях, из рассчитанных с помощью ANSYS толщинных критических частот f_Y = 939 кГц и f_А =1190 кГц, по формулам (6) можно рассчитать с₄₄ и k₁₅² . Пьезомодуль d₁₅ легко рассчитать, используя формулу

k₁₅² = с₄₄ d₁₅² / ε₁₁^T.

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

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

s₁₁ = 0,123 (0,123); s₁₂ = –0,0405 (–0,0405); s₁₃ = –0,0528 (–0,0531); s₃₃ = 0,153 (0,155);

s₄₄ =0,38 (0,39); d₃₁ =–1,23 (–1,23); d₃₃ =2,8584 (2,89); d₁₅ =5,001 (4,96);

ε₃₃^T =114.95 (115,05); ε₁₁^T = 130 (130,54).

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