УДК 537.9
ТЕОРИЯ МАГНИТОЭЛЕКТРИЧЕСКОГО ЭФФЕКТА В ДВУХСЛОЙНЫХ МАГНИТОСТРИКЦИОННО-ПЬЕЗОЭЛЕКТРИЧЕСКИХ СТРУКТУРАХ С УЧЕТОМ НЕОДНОРОДНОСТИ ДЕФОРМАЦИЙ ПО ТОЛЩИНЕ ОБРАЗЦА
Д.А.Филиппов, Т.А.Галичян
THEORY OF MAGNETOELECTRIC EFFECT IN MAGNETOSTRICTIVE-PIEZOELECTRIC BILAYER STRUCTURES TAKING INTO ACCOUNT INHOMOGENEITY OF DEFORMATIONS OVER THE SPECIMEN THICKNESS
D.A.Filippov, T.A.Galichian
Политехнический институт НовГУ, [email protected]
Представлена теория линейного магнитоэлектрического эффекта в двухслойной магнитострикционно-пьезоэлектрической структуре в форме прямоугольной пластинки с учетом неоднородностей амплитуды деформаций механических колебаний по толщине образца. На основе совместного решения уравнений эластодинамики и электростатики для магнитострикционной и пьезоэлектрической фаз получено выражение для частотной зависимости магнитоэлектрического эффекта в области электромеханического резонанса. Представлены теоретические зависимости смещений и напряжений в магнитострикционной и пьезоэлектрической фазе по толщине образца. Зависимости имеют нелинейный характер, и их учет приводит к заметному вкладу в величину эффекта.
Ключевые слова: слоистая структура, магнитоэлектрический эффект, магнитострикция, пьезоэлектричество
The theory of linear magnetoelectric effect in magnetostrictive-piezoelectric bilayer structure is presented taking into account the dependence of mechanical oscillations' strain amplitude on the specimen thickness for structures in the shape of a rectangular plate. The expression for frequency dependence of the ME effect in the region of electromechanical resonance was derived, using simultaneous solution of elastodynamic and electrostatic equations for magnetostrictive and piezoelectric phases. The theoretical dependencies of displacement and stress distributions over the specimen thickness in magnetostrictive and piezoelectric phases are presented. These dependencies are nonlinear, and their consideration leads to a significant contribution to the magnitude of the effect. Keywords: layered structure, magnetoelectric effect, magnetostriction, piezoelectricity
Введение
Магнитострикционно-пьезоэлектрические структуры привлекают к себе большое внимание благодаря тому, что в них возникают эффекты, которые отсутствуют по отдельности и в магнитострикционной, и пьезоэлектрической фазах. Одним из таких эффектов является магнитоэлектрический (МЭ) эффект. Он заключается в возникновении напряжения на обкладках конденсатора, диэлектриком которого является магнитострик-ционно-пьезоэлектрический композит, при помещении его в магнитное поле.
Величина МЭ эффекта в двухслойных магни-тострикционно-пьезоэлектрических структурах значительно больше, чем в объемных композитах [1], что позволяет их рассматривать как перспективные материалы для практического применения. При теоретическом описании МЭ эффекта в магнитоэлектрических композитах в настоящее время наибольшее распространение получили два метода: метод эффективных параметров и метод, основанный на совместном решении уравнений отдельно для магнитострик-ционной и пьезоэлектрической фаз с учетом условий на границе между фазами. Теория МЭ эффекта в объемных и многослойных композитах на основе метода эффективных параметров представлена в работах [27]. В этих работах получено выражение для МЭ коэффициента по напряжению и изучена частотная за-
висимость. Одним из недостатков метода эффективных параметров является ограниченность его применения. Он может использоваться тогда, когда размеры структурных единиц композита много меньше длины акустической волны. В этом случае композит можно рассматривать как однородную среду с некоторыми эффективными параметрами. С другой стороны, возникает проблема определения эффективных параметров. Для их расчетов используются различные модели [3-5], которые в той или иной степени отражают реальные факты.
Теория МЭ эффекта в слоистых структурах, основанная на совместном решении уравнений эла-стодинамики и электростатики отдельно для магни-тострикционной и пьезоэлектрической фаз, развита в работах [8-13]. При этом межслоевое соединение на границе раздела учитывался формально либо введением коэффициента связи между слоями [8-11], либо предполагалось, что связь идеальная и смещения пьезоэлектрической и магнитострикционной фаз одинаковы и не изменяются по толщине слоя [12-14].
Недавно в работах [15,16] была представлена теория МЭ эффекта в двухслойных магнитострикци-онно-пьезоэлектрической структурах, в которой учитывалась неоднородность смещений по толщине образца. Однако в этих работах не было представлено пространственное распределение деформаций и напряжений, а также недостаточно полно проанализи-
ровано влияние неоднородного распределения на величину эффекта.
В данной работе на основе совместного решения уравнений эластодинамики и электростатики получены выражения для зависимости амплитуды смещений и деформаций по толщине образца. На основе полученных зависимостей, с учетом условия разомкнутой цепи, получено выражение для МЭ коэффициента по напряжению при поперечной и продольной ориентациях электрического и магнитного полей. Показано, что учет неоднородности амплитуды по толщине образца приводит к заметному вкладу в величину эффекта.
Модель
В качестве модели рассмотрим структуру, состоящую из механически взаимодействующих между собой пьезоэлектрической и магнитострикционной фаз в форме прямоугольной пластинки. Будем считать, что длина пластинки L много больше ее ширины
W, в то время как толщину феррита mt и пьезоэлек-
трика pt будем считать конечной величиной (рис.1).
\ 0
/
X
Рис.1. Схематичное изображение структуры: 1 — магнитострик-ционный слой, 2 — пьезоэлектрический слой, 3 — электроды
Выберем систему координат таким образом, чтобы начало совпадало с границей раздела феррит-пьезоэлектрик, а ось Z направим вертикально вверх, перпендикулярно границе раздела. Пьезоэлектрический слой предварительно поляризован перпендикулярно плоскости контактов (ось 2). Переменное магнитное поле с частотой ю возбуждает в магнитост-рикционной компоненте упругие колебания, которые посредством сдвиговых напряжений передаются через границу раздела в пьезоэлектрическую компоненту, что приводит к возникновению взаимосвязанных колебаний магнитострикционной и пьезоэлектрической подсистем. Поскольку имеется резкая граница, через которую осуществляется взаимодействие между двумя слоями, то амплитуда колебаний будет неоднородной перпендикулярно границе раздела. Полагая пластинку узкой, в первом приближении можно считать, что вдоль оси У смещения будут однородными и отличными от нуля компонентами будут только компоненты напряжений Т и Т . Уравне-
* хх хг *
ние движения для х — проекции вектора смещения среды ам запишем в виде:
„ дм даТ даТ
а р_х = хх +__ хг
где индекс а равен соответственно т для магнитост-рикционного и р для пьезоэлектрического слоя, а р —
ат
плотность, 1 .. — напряжение магнитострикционного
V
и пьезоэлектрического слоя соответственно.
Основные уравнения для магнитострикцион-ной и пьезоэлектрической фазы имеют следующий вид:
pS.=ps.pT.+pdpEt, (2)
г V . кг к 4 /
pDk = pеk 1'Е +pdгpT,
к кп п кг г'
=т&.тТ.+mqkгmИk, г г. . 1кг к
(3)
(4)
тВк=тЧгтТ1+т^кптИп, (5)
где ^ — тензор деформаций, pEk и тИк — компоненты векторов напряженности электрического и магнитного полей, pDk и тВк — компоненты векторов электрической и магнитной индукции, 1. — компоненты тензора напряжений, ая . — коэффициенты податливости, pdki и тды — пьезоэлектрический и
пьезомагнитный коэффициенты, p екп и т цкп — компоненты тензоров диэлектрической и магнитной проницаемости.
Решение уравнения для вектора смещения среды представим в следующем виде:
а м(х, ¿^^[МсобО^-kx)+aBsm(юt -Ах)], (6) где аА и аВ — постоянные интегрирования, ag (г) — некоторая функция. В зависимости от ориентации электрического и магнитного полей далее рассмотрим два случая, которые отличаются ориентацией полей.
1. Продольная ориентация электрического и магнитного полей
В случае продольной ориентации электрического и магнитного полей (продольный МЭ эффект) магнитные поля (поле подмагничивания Иь и переменное И) сонаправлены с вектором поляризации Р. Ограничимся рассмотрением планарных колебаний, распространяющихся вдоль оси X. Переменное магнитное поле с частотой ю возбуждает в магнитост-рикционном слое упругие колебания, которые посредством сдвиговых напряжений передаются через границу раздела в пьезоэлектрический слой, что приводит к возникновению электрического поля посредством пьезоэффекта.
При данной ориентации полей уравнения тензора деформаций для пьезоэлектрической и магнито-стрикционной фаз аS ,., и г — проекция вектора электрического смещения pD имеют следующий вид:
pS =—Т + Ч pE .
хх pY ..........
хх хх,г г7
дt2
дх
дг
(1)
pS —"Т
хг pQ хг
(7)
(8)
2
У
1
3
2
pD =P е РЕ +У РТ ,
^ 22 2 ХХ,2 XX'
^ =—тТ +mq тН ,
XX тV XX *ХХ,2 2
ns = тт
Х2 тX
(9) (10)
(11)
где аУ=—— — модули Юнга, аО=—— — модули
а 5 ' ' а 5
XX X
сдвига магнитострикционной и пьезоэлектрической фазы соответственно.
Условия механического равновесия на свободных боковых поверхностях образца в точках X=^/2 дают следующие граничные условия: о т
| ТД+L/2,z)dz+|тТза(±Ь/2)йг=0. (12)
о
Используя эти граничные условия и уравнение (6), для смещения магнитострикционной и пьезоэлектрической среды получим следующие выражения:
тих = [ехр(-2т к)ехр(тхг)+exp(-mxz)]Ssin(^x), (13) \^^^х;?) - tg(p к^ш^ %2))(1+ехр(-2тк)^т(Ь:), (14)
где В=-
{тн >+pYptpd <рЕ >
1XX,2S 2' XX,2 4 2'
kcos(к)(l+exp(-2mк)) тУП 1ЬПтК)+pYpt tgppк) у к к
к = ^ /2 и ак=ах°^ — безразмерные параметры.
тХ2 =-2(1+у)
ю
"У 2
2
V =2(1+у)
г-Г
аY'
аVL — скорости продольных волн в магнитострикци-
онной и пьезоэлектрической фазе соответственно, V — коэффициент Пуассона.
Выражая компоненты тензора напряжений
тТ^ и 1'Т через компоненты тензора деформаций
из (7) и (10), используя уравнение (6), получим для них следующие выражения:
X=тУ ^В^^)(ехр(-2т к)ехр(т ^)+ехр(-т
-mqxx,2mH2], (15)
pTxx=П ^В^^)^ехр(-2т к))х
х(<га( " Х2)-tg(p к)sin(p Х2))-pdxX/E2 ], (16)
Как видно из уравнений (13) и (14), решения представляют собой плоские волны, амплитуда которых изменяется по толщине образца. Это зависимость имеет нелинейный характер и в общем случае зависит от частоты колебаний. Легко показать, что в случае низких частот, когда безразмерные параметры 1?к и тк меньше единицы, амплитуда перестает зависеть по толщине образца. Таким образом, полученные ранее в работе [12] результаты в предположении о том, что амплитуда смещений одинакова по толщине образца, имеют место только для низких частот и тонких слоев. На рис. 2 представлена рассчитанная по уравнениям (13) и (14) теоретическая зависимость смещений магнитострикционной и пьезоэлектрической фазы для двухслойной структуры никель — цирконат-
титанат свинца (№-Р2Т). При расчетах использовались следующие параметры:
для никеля — ^ = 0,3 мм, тУ=204 ГПа, тр=8900 кг/м3, mqx:z =1156-10-12м/А;
для ЦТС — ^ = 0,73 мм, Т=65 ГПа, ^=7600 кг/м3, pd =-175-10-12м/В, p е / е=1750.
г I ' XX, 2 > 22 0
Величина напряженности переменного магнитного поля при расчетах принималась равной Н = 100 Ое, частота переменного магнитного поля f = 300 кГц.
На рис.3 и 4 представлены зависимости механических напряжений в магнитострикционной и пьезоэлектрической фазах, рассчитанные с использованием уравнений (15) и (16).
Магнитоэлектрический коэффициент по напряжению определяется как отношение среднего значения напряженности электрического поля в структуре к среднему значению напряженности внешнего магнитного поля, ее вызвавшей, т.е.:
(аЕ>=(Е>/Н, (17)
где (Е>=и/ (mt+ pt) — среднее значение напряженности электрического поля в структуре, и — возникающая разность потенциалов между электродами.
Подставляя в уравнение (9) полученные выражения для тензора напряжений (16) с использованием условия разомкнутой цепи для МЭ коэффициента по напряжению, при продольном эффекте получим следующее выражение:
XX,2 1XX,2
qx
пуп
7е Дг
тут11Ь( т к) + pYpt 13 к)
,tg(к)tg(p к) ^
* к p к т+pt'
где ДL =1-Кp'
(18)
л
1—
tg(к)tg(p к)
mYmtth(mк)+pYpttg(Pк) к
к
к
/
2 Т ( Ч^ )2
К =-— — квадрат коэффициента электро-
p p е
22
механической связи.
2.Поперечная ориентация электрического и магнитного полей
В случае поперечной ориентации магнитные поля (постоянное Нь и переменное Н) направлены перпендикулярно направлению поляризации вдоль оси X. В случае поперечного эффекта в уравнениях для тензора деформаций магнетика вместо переменной тНк более удобно использовать тВк. Это объясняется тем, что при подстановке в уравнения движения среды можем использовать тот факт, что дВ
=0, поскольку divB=0. В этом случае уравнения
дт
для тензора деформаций запишутся в виде: 1
XX туВ XX
Т +т?
XX о -
В ,
XX, X X'
тВ =тд "Т +> тН ,
X IXX,X XX г XX X'
(19)
2
а E,L
X
1
ю
К
к
p
к
где
т в =msxx — податливость пьезомагнитной фазы
д^
при постоянной магнитной индукции, g =-—
хх,х дтВ
х
— пьезомагнитный коэффициент.
Проводя аналогичные расчеты, для поперечного МЭ коэффициента по напряжению получим выражение в следующем виде:
а =-
Е,Т
хх,г ° хх,х
^х
цтувт
p е ДТ
гг 1
ту В mt ^(т к) + pYpt tg( p к)
т к p к
где ДТ =1-К
,tg(к)tg(p к) ^
* к p к т+pt'
(21)
Л
tg(к)tg(p к)
туВтг 1Ь(т к) + pYpt tg(P к) к P к
кк ц — проницаемость магнитной фазы.
Результаты и обсуждение
На рис.2 представлена зависимость величины амплитуды смещения по толщине образца. Как следует из рисунка, амплитуда колебаний магнитост-рикционной фазы практически не изменяется по толщине магнитострикционного слоя. Это объясняется тем, что переменное магнитное поле возбуждает колебания магнитострикционной среды одновременно по всей толщине магнитострикционного слоя. В пьезоэлектрическом слое колебания возбуждаются посредством сдвиговых деформаций через границу раздела феррит — пьезоэлектрик. Это приводит к тому, что амплитуда колебаний пьезоэлектрического слоя значительно уменьшается с глубиной.
24 22 201816-
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 Z (тт)
Рис.2. Зависимость смещений среды магнитострикционной и пьезоэлектрической фаз по толщине образца. 1 — магни-тострикционная фаза, линия 2 — пьезоэлектрическая фаза. Частота переменного магнитного поля f = 300 кГц
На рис.3 и 4 представлена зависимость величины амплитуды механических напряжений по тол-
щине образца. В полном соответствии с теорией, они имеют максимальные значения на границе раздела и равны нулю на свободных поверхностях образца.
0,15 2 (тт)
Рис.3. Зависимость напряжения сдвига в магнитострикционной фазе по толщине феррита. Частота переменного магнитного поля f = 300 кГц
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7
2 (тт)
Рис.4. Зависимость напряжения сдвига в пьезоэлектрической фазе по толщине пьезоэлектрика. Частота переменного магнитного поля f = 300 кГц
Расчет МЭ коэффициента по напряжению с использованием полученного в данной работе уравнения (18) дает для него значение на частоте 300кГц равное аЕ^ = 513 мВ/см-Э, в то время как расчеты, проведенные по формулам, представленным в работе [12] дают завышенное значение, равное аЕ^ = 860 мВ/см-Э.
Таким образом, на высоких частотах учет неоднородности приводит к значению МЭ коэффициента по напряжению значительно отличающегося от рассчитанного в предположении, что амплитуда смещений одинакова по толщине образца.
В области низких частот и тонких слоев обе модели приводят к одинаковым значениям МЭ коэффициента по напряжению.
0,4-
0,35-
0,30-
0,25
X
0,20
0,15
0,10
0,05
1
0
0
0,05
0,10
0,20
0,25
0,3
0,4-
0,35-
0,30-
0,25-
0,20-
0,15
0,10
0,05
0
4
Заключение
Неоднородность структуры, связанная с наличием границы раздела феррит-пьезоэлектрик, приводит к неоднородному распределению амплитуды смещений и механических напряжений по толщине образца. В области низких частот эту неоднородность можно не учитывать и считать, что амплитуда не изменяется по толщине образца. Однако в области частот порядка сотен килогерц это изменение становится существенным и вносит значительный вклад в величину МЭ коэффициента по напряжению.
1. Филиппов Д.А., Лалетин В.М., Srinivasan G. Низкочастотный и резонансный магнитоэлектрические эффекты в объемных композиционных структурах феррит никеля - цир-конат-титанат свинца // ЖТФ. 2012. Т.82. №1. С.47-51.
2. Harshe G., Dougherty J.P., Newnham R.E. Theoretical modelling of multilayer magnetoelectric composites // Int. J. Appl. Electromagn. Mater. 1993. V.4. P.145-159.
3. Harshe G., Dougherty J.P., Newnham R.E. Theoretical modelling of 3-0/0-3 magnetoelectric composites // Int. J. Appl. Electromagn. Mater. 1993. V.4. P.161-171.
4. Bichurin M.I., Petrov. V.M., Srinivasan G. Theory of low-frequency magnetoelectric effects in ferromagnetic-ferroelectric layered composites // J. Appl. Phys. 2002. V.92. P.7681-7683.
5. Osaretin I.A., Rojas R.G. Theoretical model for the magnetoelectric effect in magnetostrictive/piezoelectric composites // Phys. Rev. B. 2010. V.82. P.174415.
6. Bichurin M.I., Filippov D.A., Petrov V.M., et al. Resonance magnetoelectric effects in layered magnetostrictive-piezoelectric composites // Phys. Rev. B. 2003. V.68. P. 132408.
7. Филиппов Д.А., Бичурин М.И., Петров В.М., и др. Резонансное усиление магнитоэлектрического эффекта в композиционных феррит-пьезоэлектрических материалах // ФТТ. 2004. Т.46. №9. С.1621-1627.
8. Филиппов Д.А. Теория магнитоэлектрического эффекта в двухслойных ферромагнет-пьезоэлектрических структурах // Письма в ЖТФ. 2004. Т.30. №23. С. 24-31.
9. Филиппов Д.А. Теория магнитоэлектрического эффекта в гетерогенных структурах на основе ферромагнетик-пьезоэлектрик // ФТТ. 2005. Т.47. №6. С.1082-1084.
10. Chang C.M., Carman G.P. Analytically evaluating the properties and performance of layered magnetoelectric composites // J. Intell. Mater. Syst. Struct. 2008. V.19. P.1271-1280.
11. Filippov D.A., Srinivasan G., Gupta A. Magnetoelectric effects in ferromagnetic films on ferroelectric substrates // J. Phys.: Condens. Matter. 2008. V.20. P.425206.
12. Бичурин М.И., Петров В.М., Аверкин С.В., Филиппов А.В. Электромеханический резонанс в магнитоэлектрических слоистых структурах // ФТТ. 2010. Т.52. №10. С.1975-1980.
13. Wang Y., Hasanyan D., Li M., et al. Theoretical model for geometry-dependent magnetoelectric effect in magnetostric-tive/piezoelectric composites // J. Appl. Phys. 2012. V.111. P.124513.
14. Hasanyan D., Wang Y., Gao J., et al. Modeling of resonant magneto-electric effect in a magnetostrictive and piezoelectric laminate composite structure coupled by a bonding material // J. Appl. Phys. 2012. V.112. P.064109.
15. Филиппов Д.А., Лалетин В.М., Galichyan T.A. Магнитоэлектрический эффект в двухслойной магнитострикци-онно-пьезоэлектрической структуре // ФТТ. 2013. Т.55. №9. С.1728-1733.
16. Filippov D.A., Galichyan T.A., Laletin V.M. Magnetoelectric effect in bilayer magnetostrictive-piezoelectric structure. Theory and experiment // Appl. Phys. A. 2013. doi:10.1007/s00339-013-7957-z
References
1. Filippov D.A., Laletin V.M., Srinivasan G. Nizkochastotnyi i rezonansnyi magnitoelektricheskie effekty v ob"emnykh kompozitsionnykh strukturakh ferrit nikelia - tsirkonat-titanat svintsa [Low-frequency and resonance magnetoelectric effects in nickel ferrite-PZT bulk composites]. ZhTF -Technical Physics, 2012, vol. 82, no. 1, pp. 47-51.
2. Harshe G., Dougherty J.P., Newnham R.E. Theoretical modelling of multilayer magnetoelectric composites. Int. J. Appl. Electromagn. Mater., 1993, vol. 4, pp. 145-159.
3. Harshe G., Dougherty J.P., Newnham R.E. Theoretical modelling of 3-0/0-3 magnetoelectric composites. Int. J. Appl. Electromagn. Mater., 1993, vol. 4, pp. 161-171.
4. Bichurin M.I., Petrov V.M., Srinivasan G. Theory of low-frequency magnetoelectric effects in ferromagnetic-ferroelectric layered composites. J. Appl. Phys., 2002, vol. 92, pp. 7681-7683.
5. Osaretin I.A., Rojas R.G. Theoretical model for the magneto-electric effect in magnetostrictive/piezoelectric composites. Phys. Rev. B, 2010, vol. 82, p. 174415.
6. Bichurin M.I., Filippov D.A., Petrov V.M., Laletsin V.M., Paddubnaya N.N., Srinivasan G. Resonance magnetoelectric effects in layered magnetostrictive-piezoelectric composites. Phys. Rev. B, 2003, vol. 68, p. 132408.
7. Filippov D.A., Bichurin M.I., Petrov V.M., Laletin V.M., Srinivasan G. Rezonansnoe usilenie magnitoelektricheskogo effekta v kompozitsionnykh ferrit-p'ezoelektricheskikh mate-rialakh [Resonant amplification of the magnetoelectric effect in ferrite-piezoelectric composites]. FTT - Phys. Solid State,
2004, vol. 46, no. 9, pp. 1621-1627.
8. Filippov D.A. Teoriia magnitoelektricheskogo effekta v dvukhsloinykh ferromagnet-p'ezoelektricheskikh strukturakh [Theory of magnetoelectric effect in ferromagnetic-piezoelectric bilayer structures]. Pis'ma v ZhTF - Tech. Phys. Lett., 2004, vol. 30, no. 23, pp. 983-986.
9. Filippov D.A. Teoriia magnitoelektricheskogo effekta v get-erogennykh strukturakh na osnove ferromagnetik-p'ezoelektrik [Theory of magnetoelectric effect in ferromagnetic-piezoelectric heterostructures]. FTT - Phys. Solid State,
2005, vol. 47, no. 6, pp. 1118-1121.
10. Chang C.M., Carman G.P. Analytically evaluating the properties and performance of layered magnetoelectric composites. J. Intell. Mater. Syst. Struct., 2008, vol. 19, pp. 12711280.
11. Filippov D.A., Srinivasan G., Gupta A. Magnetoelectric effects in ferromagnetic films on ferroelectric substrates. J. Phys.: Condens. Matter, 2008, vol. 20, p. 425206.
12. Bichurin M.I., Petrov V.M., Averkin S.V., Filippov A.V.: Elektromekhanicheskii rezonans v magnitoelektricheskikh sloistykh strukturakh [Electromechanical resonance in mag-netoelectric layered structures]. FTT - Phys. Solid State, 2010, vol. 52, no. 10, pp. 2116-2122.
13. Wang Y., Hasanyan D., Li M., Gao J., Li J., Viehland D., Luo H. Theoretical model for geometry-dependent magneto-electric effect in magnetostrictive/piezoelectric composites. J. Appl. Phys., 2012, vol. 111, p. 124513.
14. Hasanyan D., Wang Y., Gao J., Li M., Shen Y., Li J., Viehland D. Modeling of resonant magneto-electric effect in a magnetostrictive and piezoelectric laminate composite structure coupled by a bonding material. J. Appl. Phys., 2012, vol. 112, p. 064109.
15. Filippov D.A., Laletin V.M., Galichyan T.A. Magnito-elektricheskii effekt v dvukhsloinoi magnitostriktsionno-p'ezoelektricheskoi strukture [Magnetoelectric effect in a magnetostrictive-piezoelectric bilayer structure]. FTT -Phys. Solid State, 2013, vol. 55, no. 9, pp. 1840-1845.
16. Filippov D.A., Galichyan T.A., Laletin V.M. Magnetoelectric effect in bilayer magnetostrictive-piezoelectric structure. Theory and experiment. Appl. Phys. A, 2013. doi:10.1007/s00339-013-7957-z