磁驱动中心对称挠曲电夹层板力电耦合性能分析.pdf

sandwich plate:Chinese Journal of Theoretical and Applied Mechanics,2023,55(7):1517-1525Guo Ziwen,Zhang Gongye,Mi Changwen.On the magnetically induced electromechanical coupling of centrosymmetricflexoelectric引用格式：郭子糜长稳.磁驱动中心对称烧由电夹层板力电耦合性能分析2355(7):1517-1525固体力学Chinese Journal of Theoretical and Applied MechanicsJul.,20232023年月Vol.55,No.7力第5 5 卷第7 期报学学磁驱动中心对称曲电夹层板力电耦合性能分析郭子文章公也2)糜长稳(东南大学土木工程学院，南京2 10 0 96)摘要王现代工业的发展对材料性能和结构尺寸提出更高的要求，机电器件的设计越来越偏向于小型化、高频化和智能化.最新研究成果表明，磁电耦合复合材料不仅能够以较强的磁电转换效率实现磁能、机械能和电能之间的相互转换，还可以避免结构与机械驱动源的直接接触，实现非接触调控，这对于制备多功能微纳米器件具有重要意义.文章基于Mindlin所发展的多物理场结构理论分析方法，结合宏观压磁理论和偶应力挠曲电理论，研究由单个挠曲电电介质层和两个对称压磁层构成的三明治型夹层板在外部横向磁场驱动下的动态力电耦合响应，其中通过引入曲率将经典力电耦合理论拓宽到中心对称材料.夹层板在正弦型全局磁场和均布局部磁场驱动下的动态数值算例表明：位移和电势具有一定的频率依赖性，当激振频率达到固有频率时，振幅达到最大值；此外，对称式驱动压磁层分布方式趋于提高多层复合板的力电耦合性能.文章理论模型和研究结果可为磁控机电器件的优化设计提供新的改进思路.关键词挠曲电效应，压磁驱动，偶应力理论，夹层板，中心对称中图分类号：TB381文献标识码：Adoi:10.6052/0459-1879-23-103ON THE MAGNETICALLY INDUCED ELECTROMECHANICAL COUPLING OFCENTROSYMMETRIC FLEXOELECTRIC SANDWICHPLATE!)GuoZiwenZhang Gongye 2)Mi Changwen(School of Civil Engineering,Southeast University,Nanjing 210096,China)AbstractThe development of modern industry inspires higher requirements for material properties and structuraldimensions.The design of electromechanical devices is increasingly biased towards miniaturization,high frequency andintelligence.The most recent studies demonstrate that composite materials with magnetoelectric coupling can not onlyachieve mutual conversion of magnetic,mechanical,and electrical energy with high magnetoelectric conversionefficiencies,but can also avoid direct contact between the structure and the mechanical driving source to achieve non-contact control,which is crucial for the creation of multifunctional micro and nanoscale devices.Based on the multi-physics structural analysis framework developed by Mindlin,this paper studies the dynamic electromechanical couplingresponse of a sandwich plate composed of a flexoelectric dielectric layer and two symmetric piezomagnetic layersinduced by external transverse magnetic fields.The macroscopic piezomagnetic and curvature-induced flexoelectrictheories are employed and the classical electromechanical coupling theory is extended to centrosymmetric materials.Thedynamic numerical examples of the sandwich plate driven by a sinusoidal global magnetic field and a uniformlydistributed local magnetic field show that the magnitudes of displacement and potential are frequency dependent.When2023-03-21收稿，2 0 2 3-0 5-0 9录用，2 0 2 3-0 5-10 网络版发表，1)国家自然科学基金（12 0 0 2 0 8 6)和中央高校基本科研基金(2 2 42 0 2 2 R40040)资助项目.2)通讯作者：章公也，副研究员，主要研究方向为高阶弹性理论和多场耦合问题.E-mail:力15182023年第5 5 卷报学学the excitation frequency reaches the natural frequency,the amplitude reaches the maximum.In addition,the distributionof symmetrical piezomagnetic layer tends to improve the electromechanical coupling performance of multilayercomposite plates.Both the theoretical model and numerical results provide new ideas for the optimization design ofmagnetic-controlled electromechanical devices.Keywordsflexoelectric effect,piezomagnetically driven,couple stress theory,sandwich plate,centrosymmetric引言力电耦合效应是广泛存在于各类电介质材料中的多场耦合效应之一，常见的形式有压电效应、电致伸缩效应、铁电效应等1-5 .一些线性力电耦合效应，如压电效应和挠曲电效应，由于可以实现电能和机械能之间的相互转化而被广泛应用于多种机电器件和设备中，例如谐振器、传感器、发电机和滤波器等6-10 挠曲电效应是电介质中应变梯度（非均匀应变)与电极化之间的耦合(正挠曲电效应),或者电极化梯度与应力的耦合关系(逆挠曲电效应).正挠曲电效应将机械能转化为电能，可用于制造传感元件或发电元件.逆挠曲电效应能够将电能转化为机械能，可应用于制造驱动元件.在传统压电材料中，均匀应变(如单轴拉压)能够通过压电效应激发电极化,其中，非中心对称晶体是材料表现压电性的必要条件；而挠曲电效应即使在中心对称晶体中，也可通过非均匀应变(如弯曲变形)局部破坏材料内部正负电荷的对称性并诱发挠曲电极化.近年来，挠曲电效应逐渐引起了研究者们的广泛关注.相比于压电效应，挠曲电效应具有明显的优势.一方面，压电效应仅存在于2 0 种非中心对称的介电材料中，这极大限制了压电材料的进一步发展与应用.而非均匀应变能够在所有介电材料中产生电极化，即挠曲电材料的来源比压电材料更广泛.另一方面，由于压电材料在高于居里温度时会发生相变，导致材料对称性增高，压电性消失，因而压电效应很难应用于高温场合，而挠曲电效应没有这个限制，故挠曲电材料的工作温度范围更广.除此之外，挠曲电效应会随着尺度的降低而急剧增强，当材料特征尺寸达到微纳米量级时，挠曲电效应会变得十分显著，因此在微纳米材料中挠曲电效应变得更加不可忽视11-12 ,在考虑挠曲电效应的连续介质力学模型中，机械载荷引起的弯曲变形和曲率引起的挠曲电极化是大多数研究人员关心的问题13-2 0 .除此之外，利用热弹性、热电性或磁致伸缩效应引起的热或磁致变形和电极化在多场耦合领域也已得到广泛应用2 1-2 4.特别地，与直接施加机械载荷相比，利用磁电耦合材料进行磁调控是一种非接触式的智能驱动方法，通常出现在以压磁性材料作为驱动器的复合结构中，能够实现磁-力-电的多场耦合.磁电耦合材料具有很多重要的应用，包括磁传感器、数据存储、医学中的远程药物传递和无线能量采集器等2 5-2 7 .虽然自然界中的磁电耦合单相材料确实存在，如BiFeO3,Cr2O3和YMnO3，但是，它们具有较低的居里温度并且在室温下的磁电耦合性能很弱，这无疑限制了磁电耦合材料的应用2 5 .为此，人们转而采用压电和磁致伸缩复合结构，用以规避单一构型磁电材料的缺点,如BaTiO3和CoFe2O4磁电复合材料2 8 、PZT压电陶瓷与稀土合金Terfenol-D磁电复合材料等2 9.近期,Li等3 0 针对一维非压电悬臂梁的研究表明,通过磁场也可以激发挠曲电极化，这进而拓宽了挠曲电器件的应用范围.现有文献对于挠曲电效应的结构理论研究大多只针对弯曲和屈曲等静态问题，缺乏对于较为复杂的受迫振动等动态问题的研究，尤其是将磁场作为驱动源时的动态振动机制尚不明晰.有鉴于中心对称挠曲电介电材料在微纳米尺度下显著的力电耦合效应和磁场的非接触驱动特性，本文以具有单个挠曲电层和两个对称压磁层组成的夹层板为例，通过多物理场结构理论建模，研究夹层板在全局和局部磁场驱动下的动态力电耦合响应，以期为磁控机电器件的优化设计提供新的改进思路.1结构理论框架从连续介质力学的角度看，对于特定的模型结构，可以将三维框架简化为一维的梁、杆或二维的板模型，再通过解析法对这些模型进行求解.因此，本节将基于Mindlin发展的多物理场结构理论分析1519第7 期郭子文等：磁驱动中心对称挠曲电夹层板力电耦合性能分析方法建立挠曲电夹层板的二维理论分析框架1.1单层曲电板如图1所示的单层挠曲电薄板，长度和宽度分别为和b,厚度为2 h.采用直角坐标系(x1,x2,x3)，在外部横向力q(x1,x2)的作用下发生变形.根据Mindlin结构理论假设，板内任意质点的位移可沿板的厚度方向按级数形式展开，其一阶截断的表达式为3 1-3 2 ui(x,t)=uo(x1,x2,)+x31(x1,X2,t)(0)(1)u2(x,t)(x1,X2,t)+X3u2)(x1,X2,t)(1)(0)u3(x,t)(X1,X2,t)其中，u1,u2和u3分别表示点(x1,x2,x3)在t时刻沿x1,2和方向的位移，和u表示板中面上一点(x1,x2,0)沿面内两个方向的位移，u)表示板中面上一点的挠度，）和u）代表薄板直法线在两个竖直坐标面内的转角.与此相应，板内的电势也可以表达为9(x,t)=g(0)(x1,X2,t)+x3g(1)(x1,X2,t)(2)其中，s(0)和(1)表示沿板面的零阶电势和沿板横向的一阶电势.根据经典偶应力理论3 3-3 6 1,无穷小应变张量定义为8=(Vu+uV)/2；曲率张量表达式为=V,其中(=(Vu)/2）为旋转矢量.结合式（1)可得非0应变和曲率分量分别为(0)(1)811（O）8128212813831(3)2(0)8222,2823三832一(0)(0)X1123,12X1223,22(1)(0)1(1)(0)X212X2221,23,12(0)(0)(1）X31(4)22,11.11(0)(0)1X32（u.X3212121X3312hX(O)3，Ya(O)u2h图1单层板及采用的坐标系Fig.1A single-layer plate and coordinate system此外，静电场强度可由电势的负梯度计算得到，即E=-Vg.结合式(2)可得其分量形式为E1=-0-0.1 x3,E2=-02)-0-9 x3,E3=-s(1)(5)(0)()1)根据偶应力挠曲电理论3 6-3 7 ,单层挠曲电板的Gibbs 自由能(G)在0,T时间内的一阶变分为SGdt=JJ,(oijoeij+mijoxij-D;oEi)dVdt(6)其中，和mi分别表示与应变共轭的对称应力和与曲率共轭的偶应力，D;为电位移，V代表板的体积，dV为板内任意微分单元，dt表示时间微分增量.对应于上式的本构方程为3 8 gi,=Cijkiekl,mij=-fkijEk，D;=S i j Ej+f i k i Xk l (7)其中，Cijkl为挠曲电板的4阶弹性刚度张量、Sij为介电系数张量、fikl为3 阶挠曲电系数张量，反映电极化和曲率之间的力电耦合效应1.2夹层板的三维本构关系考虑图2 所示的三明治式压磁驱动挠曲电夹层板结构，该模型建立在直角坐标系(x1,x2,x3)之中，由中间的挠曲电层和沿着挠曲电层的中面呈上下对称分布的两个压磁层组成.其中，挠曲电层选用具有中心对称的立方晶体结构(m3m点群)材料，具有挠曲transversemagneticfieldH,X3X2c/22h-c_c/20(l)piezomagneticlayer(2)flexoelectriclayer(3)piezomagneticlayer图2 磁场作用下的夹层板模型Fig.2A sandwich plate model under a magnetic field力15202023年第5 5 卷报学学电性.压磁层选用具有横观各向同性的六方晶体结构（6 mm点群）材料，不含压电性和磁电耦合效应，仅考虑压磁性能.两侧六方晶体压磁层的c轴虽均沿x3方向，但是二者极化方向相反.在横向磁场H3作用下，当其中一个压磁层在面内由磁致伸缩效应产生拉伸变形时，另一个压磁层则会产生压缩变形,夹层板结构整体的拉压变形得以相互抵消，反之亦然.对于立方晶体结构材料，挠曲电层的三维本构关系为3 8 i=C11e11+C12822,02=C12811+C118220i;=2C44e13,02=2C4e23 0li2=2C412m13=-m31=-fi6E2,m23=-m32=fi6E(8)m12=-m21=fi6E3,Dl=S1Ei+fi6X32D2=S11E2-fi6X31,D3=S11E3+fi6(21-X12)其中,C1l=C-Ci2/C1,C12=C12-Ci2/Cl1 为薄板的等效弹性刚度系数，二者通过令%0 求得,基于该等效方法(应力松弛处理)能够有效解释沿板厚度方向的泊松效应.对于横观各向同性结构材料，上部压磁层的三维本构关系为3 90il=C111+C122-h31H3(9)Di=SinE1,D2=SifE2,D3=$33E3其中,Ci1=Cii-Ci3/C3,C12=Ci2-Ci3/C33 为压磁层的等效弹性刚度系数，h31（=h 3 1-C i 3 h 3 3/C 3 3）对压磁系数hs,进行修正.式(9)中的上标“表示该材料参数仅适用于压磁层下部压磁层的c轴方向为负时，压磁系数的符号发生改变，因而其面内正应力本构方程变为0il=C11811+C12822+h31H3)(10)0i2=C(1211+C(122+h31 H J1.3夹层板的控制方程和边界条件根据Mindlin的电介质变分原理可同时推导出夹层板的控制方程和边界条件3 8 ,而已有研究表明3 8,40 ,一阶电势由电极化与曲率(由弯曲变形引起)的耦合直接产生，零阶电势再与一阶电势耦合.由于结构沿厚度方向具有对称性，横向磁场作用下三明治夹层板的零阶电势s(0)与一阶电势s(1)得以解耦，面内位移和相互抵消并与挑度和转角和发生解耦.因此，下文将着重关注横向磁场驱动引起的挠度、转角和横向电势分布，夹层板的控制方程为3 8 S(0)S(O)1(0)1(0)1(0)Ommm13,123,21221,11十11,121212,2221(0)f(O)(0);(0)m十m222,123S(O)S(1)S(1)一(0)1(0)1(0)mm22.2m+1311,112,212121233,2221,11(1)1(1)(1)(2);(1)12十mm121232,221S(0)S(1)S(1)1(0)1(0)1(0)mm十2312.111.112.2m22,222233,11(1)一(1)1）(2);(1)mmm231,111232,122DQ)(1)D(0)2,23(11)对应的边界条件为S(O)1(0)1(0)1(0)(0)(0)十mssmsn.nmnn,s0orn322S23S(1)1(1)十=msn1(0)0orun(1)认n(1)nn2m3n,s2S(1)1(0)1(1)1(1)1(0)mnn0or12mm(12)33123n,nmns23s,s12(1)(1)(0)(0)(0)Ws0orUSsn3n3,n(1)orus,n(1)(1)m0us,n3nDl)=0 or g(1)=g(1)其中，顶标“.表示物理量对时间的二阶导数，“_ 表示给定值.s为n阶应力，mf为n阶偶应力，Dm为n阶电位移，Jin)为n阶外力，mn)表示n阶转动惯量，它们被定义为(0)S(1),m(0),m(2)=,mij,x3mij,Di,xDi,p,xgpdx3(13)其中，表示质量密度.式（13）中的物理量建立在（n,S,x 3）下，由坐标系（x1,x2,x3）转换得到41),将式(3)式(5)、式(8)式(10)代入式(13),再将所得二维本构代入式(11)，可得到以位移u，w，电势s以及磁场强度H表示的夹层板控制方程，即1521郭子文等：磁驱动中心对称挠曲电夹层板力电耦合性能分析第7 期0（0)+(O)1)+C2+66(0)之(1)2-s1(0)22(14)其中(0)332(2)2h+C32+311一11113232(15)(2)2C32C-h+31231221+312(2)2C32C3-h+一C4h44663662人32232(2)1(-h+%C+h3+一S111132322受迫振动分析为了论证上节建立的挠曲电夹层板二维理论模型，本节将分别对其在全局和局部横向磁场作用下的动态力电耦合响应展开研究.以简支板闭合电路状态为例（边界配置如图3 所示），边界条件可根据式(12)得到,即(0)S(1)1(1)1(0)0,m十=msn3nn23n,s2(16)(1)(0)(1)us0,msn0,m3n0,二ba图3 横向磁场驱动下夹层板的边界配置Fig.3Boundary configuration of sandwich plate under a transversemagnetic field对于矩形板,沿x1=0和x1=a,边界条件为(1)0,(0)(1)=0(17)uYS(1)1(1)1m2(0)(0)(1)0,m210,m31=0m11231,22(18)将二维本构代入式（18)可进一步得到(19)沿x2=0和x2=b,边界条件为(1)=0,(0)=0,6(1)=0(20)S(1)1(1)1(0)()=0,(1)m/20,m/2m32：022m232,12(21)将二维本构代入式(2 1)可进一步得到+(2),(1)(2),(1)(22)本节主要分析外部磁场对模型的驱动作用，不考虑体力项，即令=()=)=0.夹层板受到时谐磁场的作用，磁场强度表达式为88H3=Hmn sin(5mx1)sin(Snx2)eior(23)m=1 n=1其中，5 m=m元/a，S n=n 元/b，是板的激振频率，i是虚数单位并满足i2=-1.Hmn代表磁场强度正弦级数展开式的系数，可通过积分计算得到，即Am=J1H3 sin(smx1)sin(nx2)dxidx2(24)取决于外加驱动磁场的形式，位移场和电势场也必然具有相应的时空分布特征，即88ww(0sin(tmx1)sin(nx2)eiorl(0)mnm=1 n=188V(1)cos(5mx1)sin(5nx2)eiotUmnm=1n=1(25)88(1）V(1)sin(smx1)cos(Cnx2)eiotmnm=1 n=188V(1)6(1)mnsin(smx1)sin(Cn.x2)eiorm=ln=1其中，ww(),U(),v和ow(是与m,n有关的傅mnmnmnmn力1522学报2023年第5 5 卷学里叶级数待定系数.对于任意Hmm，W v(O)，U v(1)mnv(1)和 v(),由式(2 3)和式(2 5)给出的展开式均能mn满足边界条件(17)和式(19)式(2 1).将式(2 3)和式(2 5)代入控制方程（14)得线性方程组,V(O)mn0UV(1)Smmn(Kij+w?Mi)(26)VV(1)mn3.V(1)0mn其中，Kij和Mij分别为44的刚度矩阵和质量矩阵，二者的非零分量为KiI=-C(0(s%+mK12=K21=-C(O)44SmKi3=K31=-C0snK14=K41=-fl0(Cm+$%)/2K2=-C-c(2)2-C2211Sm66SnK23=K32=-(C（2)(27)66K24=K42=fl5m/2K33=-C(0)-C(2C(2)24466 Sm11SnK34=K43=fflon/216(2)2+S33(O)M11=m(0),M22=M33=m(2)通过对方程组（2 6)中的激振频率进行扫频分析，可以求出满足条件的wv(O)，v()，v v)和 v(1),mn再将这些求出的傅里叶系数代回式(2 5)，进而可计算出u，u ，u 和(),并完成求解为了求解数值结果，本节选取的板尺寸参数为：h=10nm,c=0.4h,a=b=40h.压磁层和挠曲电层材料分别选取CoFe2O4和Si,二者的物理参数如表1所示40,42 .值得注意的是，由于CoFe2O4的挠曲电系数无法在现有文献中得出，本文忽略了压磁层的挠曲电效应43 ,此外，为了避免数值计算的奇异性，本节在数值计算中采用复刚度法，通过引入复数弹性刚度，以考虑机械阻尼的影响.即令C)=(1-i/Q)C）（=0,2)，并以c)代替刚度矩阵Kij中的c进行以上数值计算.这里Q为材料品质因数，本文取为Q=10044,表1 CoFe2O4和Si 的材料参数Table1Material parameters of CoFe2O4 and SiProperty/UnitCoFe204SiCi=286Ci2=173C11=165.7elasticconstants/GPaCis=170C12=63.9C%;=269.5C44=79.56C44=45.3S11=0.08dielectric constants/(C2(Nm)-l)x10-9511=0.10353=0.093h5i=580.3piezomagnetic constants/(N:(A-m)-l)h33=699.7flexoelectriccoefficient/(nC-ml)fi6=0.4mass density/(kg:m-3)p=5300p=23322.1全局磁场驱动本小节研究全局磁场驱动下复合板的受迫振动问题，考虑幅值为ho=1.010A/m的正弦型磁场强度分布形式HY=ho,Hn=0(m#1 or n#1)(28)图4展示了单侧压磁层驱动挠曲电双层板（蓝色线条)(参见文献40 )与当前夹层板模型(红色线条)在正弦型全局磁场驱动下的一阶电势幅频关系对比，两种模型的材料组分占比和几何尺寸均保持一致.从图中可以看出，一阶电势表现出明显的频率依赖性，当激振频率远离固有频率时，一阶电势的振动幅值几乎为零，而当激振频率靠近固有频率时，一10615A:(7.53 109,1.35 107)12B:(7.51 109,1.19 107)(i_.A)(z/qz/D)n)0l9current sandwich plate6doublelayercompositeplate300246810o/(rads-l)109图4全局磁场驱动下的一阶电势幅频关系对比Fig.4 Comparison of amplitude-frequency curve of first-order potentialunder the global magnetic field1523期郭子文等：磁驱动中心对称曲电头电耦合性能分析第阶电势的振动幅度急剧增加，在固有频率处达到最大值.此外，相较于“双层复合板”、“当前夹层板”的一阶电势明显增大，峰值提高了10%以上，这表明对称式驱动压磁层分布方式趋于提高多层复合板的力电耦合性能.2.2局部磁场驱动本小节采用H3=hoexp(iwt)的均布局部磁场驱动复合板，磁场的驱动范围是一个长度为2 d（d a/2),宽度为2 e(eb/2)的矩形区域，如图5 所示的灰色区域,关于矩形板的中心(xo,yo)对称.对于该均布局部磁场,式(2 4)中的傅里叶系数可计算为4hoHmn=(cos5m(xo-d)-cos 5m(xo+d).mn元2(cosSn(yo-e)-cosn(yo+e)sin(Emx1)sin(nx2)(29)局部磁场的具体驱动范围可通过改变d和e进行调整，为了便于进行数值计算，本小节取d=e=25nm.图6 展示了局部磁场驱动下复合板的挠度振动幅值与激振频率关系,其中,ho=1.010A/m.在式(23）中共保留了3 0 项，此时的计算结果与保留60项时在小数点后第3 位依然一致，表明在驱动磁场的傅里叶级数中取3 0 项可以满足收敛性要求.从图6 中可以看出，在0 5.0 10 1rad/s的激振频率范围内，一共出现了3 阶固有频率，相应产生了3 个共振峰，并且，随着固有频率的增加，共振峰值越来越低.受到复合板振动阻尼的影响，固有频率不再表现为单一的频率值,而是表现为一小段固有频率区间.与图6 的前两阶固有频率相对应，图7 给出了夹层板的挠度、转角和横向电势振型.虽然这些振型实际上是复值，但由于它们的虚部远小于实部，因X2b2dyo2eXXoa图5 局部驱动磁场作用区域(灰色区域)Fig.5 Local magnetic field area(grey area)2.42.00,=5.4610rad/swu/l(z/qz/0)(ojnl1.61.2W2=2.58 101 rad/s0.80.4W;=4.65101rad/s0012345o/(rads-1)1010图6 局部磁场驱动下夹层板的挠度幅频关系Fig.6Deflection amplitude-frequency curve of the sandwich plateunder local magnetic field此可略去虚部.由于矩形板几何尺寸和磁场作用范围的对称性,夹层板的转角)与心 的分布类似,出于简化的目的，必 的分布未做展示.从图7 中可以看出，这些振型具有很强的频率依赖性，低阶振型具有更大的振动幅值，而高阶振型沿复合板长度和宽度方向具有更短的周期.特别地，从图7(b)和图7(f)中还可以明显看出，由于磁场的局部驱动特征，在作Reugo/nmReus/nm0.8-0.2-0.60.4uru/55400-1.044003004003004003000200200300 x./nm200-1.4X./nm20010010010010000 x/nm00 x/nm-1.8-0.4(a)挠度(=o.)(b)挠度(a=2)(a)Deflection(=o,)(b)Deflection(o=2)图7 局部磁场驱动下夹层板的前两阶挠度、转角和一阶电势振型Fig.7The first two order deflection,rotation angle and first-order potential mode of the sandwich plate under the local magnetic field力1524学报学2023年第5 5 卷Reu0.015Reul0.0100.0060.005Wu/5wu/0.002-54004003004003004000.002+./nm200300-0.005X./nm2003002002001001000.00610010000 x/nm00 x/nm-0.010-0.015(c)转角(=)(d)转角(=2)(c)Rotation angle(=)(d)Rotation angle(=w2)105105Reo()/(V.m-)Reo(/(V-m-l)3.5-12.5-3Wu/5u/1.5400400300400-53004000.5300200300X2/nm200X./nm200200100100100100-0.500 x/nm一700 x/nm-1.5(e)一阶电势(=w)(f)一阶电势(a=2)(e)First-order potential(o=w,)(f)First-order potential(o=2)图7 局部磁场驱动下夹层板的前两阶挠度、转角和一阶电势振型（续）Fig.7The first two order deflection,rotation angle and first-order potential mode of the sandwich plate under the local magnetic field(continued)用区域中心附近产生挠度、转角和电势峰值.图7所展示的结果对板型磁驱动器（磁能转换为机械能）和板型磁传感器（磁能转换为电能）的优化设计具有重要的参考价值，3结论本文以探索新型磁控纳米机电器件为研究背景，提出了一种新的三明治式压磁驱动挠曲电夹层板模型，构建了夹层板的二维控制方程和边界条件，并分析了矩形简支板在全局和局部磁场激励下的动态力电耦合响应.基于本文的数值研究结果可以获得以下两点重要结论.(1)在外部时谐磁场驱动下，夹层板的位移和电势振动幅值依赖于激振频率.随着固有频率的增加，振型幅值迅速衰减，振型表现出周期变化特征.(2)驱动压磁层的分布方式对复合板结构的力电耦合响应具有重要影响，本文设计的对称式驱动压磁层分布方式趋于提高复合板的力电耦合性能.本文的理论模型和研究结果可为磁控机电器件的优化设计提供新的改进思路.参考文献1李东，袁惠群.磁致伸缩换能器耦合磁弹性模型与振动特性分析.固体力学学报，2 0 11,3 2(4)：3 6 5-3 7 1(LiDong,YuanHuiqun.Ana-lysis on coupling magneto-elastic characteristic of a giant magneto-strictive transducer.Chinese Journal of Solid Mechanics,2011,32(4):365-371(in Chinese)2王歆钰，储瑞江，魏胜男等.应力作用下EuTiO，铁电薄膜电热效应的唯象理论研究.物理学报,2 0 15,6 4(11):117 7 0 1（WangXinyu,Chu Ruijiang,Wei Shengnan,et al.Phenomenological theory for in-vestigation on stress tunable electrocaloric effect in ferroelectric Eu-TiO films.Acta Physica Sinica,2015,64(11):117701(in Chinese)3 Qu Y,Jin F,Yang J.Temperature-induced potential barriers inpiezoelectric semiconductor films through pyroelectric and ther-moelastic couplings and their effects on currents.Journal of AppliedPhysics,2022,131:0945024张伟，刘爽，毛佳佳等.磁耦合式双稳态宽频压电俘能器的设计和俘能特性.力学学报，2 0 2 2,5 4(4):110 2-1112(ZhangWei,LiuShuang,Mao Jiajia,et al.Design and energy capture characteristicsof magnetically coupled bistable wide band piezoelectric energy har-vester.Chinese Journal of Theoretical and Applied Mechanics,2022,54(4):1102-1112(in Chinese)5 Qu Y,Pan E,Zhu F,et al.Modeling thermoelectric effects in piezo-electric semiconductors:New fully coupled mechanisms for mech-anically manipulated heat flux and refrigeration.InternationalJournal of Engineering Science,2023,182:1037756 Kagawa Y,Yamabuchi T.Finite element simulation of two-dimen-sional electromechanical resonators.IEEE Trans.Sonics Ultrason.,1525期郭子文等：磁驱动中心对称挠曲电夹层板耦合性能分析第1974,21(4):275-2827 Devoe DL.Piezoelectric thin film micromechanical beam resonat-ors.Sensors and Actuators A-Physical,2001,88(3):263-27288 Jiang X,Huang W,Zhang S.Flexoelectric nano-generator:Materi-als,structures and devices.Nano Energy,2013,2(6):1079-10929Ray MC.Analysis of smart nanobeams integrated with a flexoelec-tric nano actuator layer.Smart Materials and Structures,2016,25(5):05501110 Chu Z,Pourhosseiniasl M,Dong S.Review of multi-layered magne-toelectric composite materials and devices applications.Journal ofPhys