搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于温度效应的无限长压电圆杆纵波分析

陈琼 薛春霞 王勋

Li Yong-Dong, Yang Wen-Jin, Zhang Na, Cui Wan-Zhao, Liu Chun-Liang. A combined phenomenological model for secondary electron emission. Acta Phys. Sin., 2013, 62(7): 077901. doi: 10.7498/aps.62.077901
Citation: Li Yong-Dong, Yang Wen-Jin, Zhang Na, Cui Wan-Zhao, Liu Chun-Liang. A combined phenomenological model for secondary electron emission. Acta Phys. Sin., 2013, 62(7): 077901. doi: 10.7498/aps.62.077901

基于温度效应的无限长压电圆杆纵波分析

陈琼, 薛春霞, 王勋

Longitudinal wave analysis of infinite length piezoelectric circular rod based on temperature effect

Chen Qiong, Xue Chun-Xia, Wang Xun
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 利用有限变形理论, 以无限长压电圆杆为研究对象, 考虑了在横向惯性、等效泊松比效应以及在热电弹耦合共同作用下, 基于Hamilton原理, 并引入Euler方程推导出压电圆杆的纵向波动方程. 采用Jacobi椭圆函数展开法, 求解压电圆杆的波动方程和对应的解. 最后, 通过Matlab软件得到不同波速比下的色散曲线, 以及温度场对压电圆杆的波形、波幅和波数的影响曲线. 数值分析结果表明: 随着温度的升高, 波速逐渐降低, 温度场的改变可影响和控制孤立波的传播特性.
    Piezoelectric elements have been commonly used because of their wide applications in sensors, transducers, and some micro intelligent structures. However, in the fields of aviation, aerospace, and automation, some relevant equipment works in a harsh environment and is susceptible to the temperature change, thereby leading its performances to be greatly affected. Therefore, the problem of nonlinear wave relating to piezoelectric circular rods in different temperature fields is studied by modeling and numerical analysis. Firstly, based on the theory of finite deformation, we take infinite piezoelectric circular rod as a research object and consider the effects of transverse inertia and equivalent Poisson's ratio under the thermoelectric coupling action. Using the Hamilton principle and introducing the Euler equation, the longitudinal wave equation of piezoelectric circular rod is obtained. Secondly, Jacobi elliptic cosine function and Jacobi elliptic sine function expansion method are used to solve the wave equation of the piezoelectric circular rod, and the solitary wave solution and the exact periodic solution of the wave equation are obtained. It is found that the periodic solution can be reduced into a solitary wave solution under certain conditions, and it is proved theoretically that there may be solitary wave stably propagating in a piezoelectric circular rod. Finally, the dispersion curves of different wave velocity ratios and the curves about influences of temperature field on the waveform, amplitude and wave number of the piezoelectric rod are obtained by Matlab. The numerical results show that the wave velocity decreases with the increase of temperature when the wave velocity ratio is constant. Given the temperature is constant, it can be found that with the increase of the ratio, the amplitude of solitary wave gradually increases while the wavelength gradually decreases. In addition, the images obtained show that although temperature change can cause the characteristics of solitary waves to change, the solitary waves are always symmetrical bell shaped waves in the propagation process, reflecting the stability characteristics under the combined action of nonlinear and dispersion effects. Therefore, the variation of temperature field can influence and control some propagation characteristics of solitary waves. Moreover, the wave theory has been widely used in the nondestructive testing of structures and the improving of information transmission quality due to its special stability.
      通信作者: 薛春霞, xuechunxia@nuc.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11202190)资助的课题
      Corresponding author: Xue Chun-Xia, xuechunxia@nuc.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11202190)

    随着现代科学技术的飞速发展, 压电材料因其独特的性能而被广泛应用[1]. 在航空航天、智能结构等领域, 装备工作环境恶劣, 差异大, 特别是温度变化比较大, 对实施控制有着非常大的影响. 因此在实际工况条件的精确建模过程中, 需要考虑电场和温度场的耦合作用. 而杆作为常用的构件, 也吸引了不少学者的关注. 刘延柱等[2]和He等[3]利用广义热-弹耦合理论求解并研究了半无限压电杆的边值问题. 目前对于波动的问题, 主要采用多尺度法、齐次平衡法等展开了对压电圆杆波动的研究[4,5]. 冯依虎[6]利用泛函分析变分迭代的方法求出各次孤子波近似解, 进而研究了强非线性波动方程的行波解. Guo等[7]利用Hamilton变分原理, 根据有限变形理论的拉格朗日描述, 推导出弹性细杆的非线性波动方程, 利用多尺度法得到了稳定的行波解. 李敏等[8]通过对薛定谔方程的相平面分析, 约化得到其同异宿轨道, 并在相应条件下得到方程的明、暗孤立波解.

    然而上述求解方法具有一定的局限性, 只能求出波动方程的冲击波解、孤波解和初等函数的周期解[9-11]. 但采用Jacobi椭圆函数法便可求出波动方程的广义周期解和对应的孤立波解[12,13], 刘志芳和张善元[14-16]利用Jacobi椭圆函数展开法求得了无限长圆杆的非线性扭转波解、孤波解以及非圆截面杆的行波解、周期解.

    由于压电结构在工程领域具有广泛应用, 压电介质中波的传播吸引了很多学者的关注. 邓庆田等[17]用逐步近似法对位移函数进行了假设并通过变动参数法求解, 对压电层和圆杆中的几何非线性波进行了研究. Seadawy和Manafian[18]通过扩展尝试方程法和积分方法, 推导出了磁电弹圆杆纵波方程的暗孤子、亮孤子、孤波、周期孤波、有理函数解和椭圆函数解等不同形式的新的显式精确解. Baskonus等[19,20]利用Sin-gordon展开法对磁电弹性圆杆的纵波方程的解析解进行了研究, 得到了更多新的解析解, 给出了所有解的数值模拟, 很好地解释了一些实际物理问题. Wang[21]研究了压电耦合圆柱壳结构中波的传播, 从理论上得到了双模壳模型的频散曲线, 推导出波数极限情况下的截止频率和相速度. Xue和Pan[22]考虑几何非线性以及横向泊松比引起的弥散效应对无限长磁电弹圆杆进行了研究, 建立了纵波方程, 并通过Jacobi椭圆函数法对其进行了求解. Samsonov[23]首先报道了杆中存在孤波的实验研究, 利用聚苯乙烯的弹性介质, 设计了一套通过光学原理构造和记录孤波的实验方法, 用全息照相法记录下了孤波轨迹, 从事实上证实了弹性固体中孤波的存在. 2013年, Toffoli等[24]在一个大的定向波池实验中探究了平面波对斜摄动的调制和有限水深下异常波动的产生, 并采集了实验数据, 对流体中的波动进行了实验探究.

    综上所述, 由于非线性波的激发和观测是非常困难的, 导致实验上工作比较少. 所以本文采用建立模型和数值分析的方法, 研究不同温度场下的非线性波动问题. 通过Hamilton变分原理, 引入Euler方程, 采用Jacobi椭圆函数展开法, 推导出压电圆杆的波动方程和对应的解, 并讨论温度变化对压电介质的波形、波幅以及波速等的影响.

    图1为无限长压电圆杆示意图, 建立圆柱坐标系(r,θ,z), z是沿着杆的轴向, 也是波传播的方向. θ = [0, 2π], 0 ≤ rR, 其中R为压电圆杆半径. 为了研究方便, 假设: 1)变形后, 板中初始状态垂直于中心平面的点仍然垂直于中心平面; 2)杆的截面是轴对称的, 即Uθ=0/θ=0, 其中Uθθ方向位移; 3)考虑泊松比效应, 纵向位移U和径向位移Ur之间满足Ur=veffrU/z, 其中veff是有效泊松比.

    图 1 压电圆杆示意图\r\nFig. 1. Schematic diagram of piezoelectric rod.
    图 1  压电圆杆示意图
    Fig. 1.  Schematic diagram of piezoelectric rod.

    对于横观各向同性的压电材料圆杆, 在考虑温度效应时的本构方程如下[25]:

    {σr=c11εr+c12εθ+c13εze31Ezλ11Θ,σθ=c12εr+c11εθ+c13εze31Ezλ11Θ,σz=c13εr+c13εθ+c33εze33Ezλ33Θ,τrz=c44γrze15Er,τθz=c44γθze15Eθ,τrθ=c66γrθ,Dr=e15γrz+ε11Er+d1Θ,Dθ=e15γθz+ε11Eθ+d1Θ,Dz=e31εr+e31εθ+e33εz+ε33Ez+d3Θ,
    (1)

    其中σi是法向应力; τij是切向应力; εi是法向应变; γij是切向应变; Ei是电场; Di是电位移; cij是弹性常数; εij是介电常数; eij是压电耦合系数; di是热电耦合系数; λii是热机耦合系数[26,27], λ11=(c11+c12+c13)α1, λ33=(2c13+c33)α3; αi为热膨胀系数; Θ为相对于初始温度T0的温度增量.

    有限(非线性)弹性应变位移关系为

    εr=Urr,εθ=Uθrθ+Urr,εz=Uz+12(Uz)2,γrθ=Urrθ+UθrUθr,γθz=Urθ+Uθz,γrz=Ur+Urz.
    (2)

    由于这是一维问题, 杆的横向边界的牵引力应该为零. 因此可以得出σr=0, τrz=0, τrθ=0, Dr=0, 从中可以得到如下关系:

    Dθ=0,γrz=γθz=0,Er=Eθ=0,
    (3)
    εr=e31Ez+λ11Θc12εθc13εzc11.
    (4)

    根据广义Hamilton变分原理可得

    δtt0zz0Ldzdt=δtt0zz0(TEP+We)dzdt=0,
    (5)

    式中L为Lagrange密度函数, T为系统的动能, EP为系统的势能, We为系统的电能, 具体表达式分别为

    T=12VρdV(Ut)2+12VρdV(Urt)2,
    (6)
    EP=12VSTQdV,
    (7)
    We=12VETDdV,
    (8)

    这里ρ为压电材料的密度, V为压电材料的体积, S={εr,εθ,γrθ}TQ={σθ,σz,τθz}T分别表示材料的应变向量和应力向量, E={0,0,Ez}T表示电场向量, D={0,0,Dz}T表示电位移向量.

    根据Euler方程, 若L=L(U,Uz,Ut,Uzz,Utt,Uzt,), 则

    LUzLUztLUt+2z2LUzz+2t2LUtt+2ztLUzt=0,
    (9)

    此处, Uz=Uz,Ut=Ut,Uzt=2Uzt, 其余类同.

    根据以上关系可以得到如下表达式:

    {LU=0,zLUz=12V[(2k1v2eff+2k2v2eff+2k3+k4+k5)(Uz)+3(12k2veff+k5)(Uz)2+k5(Uz)3],tLUt=ρV2Ut2,  2z2LUzz=0,2t2LUtt=ρV2t2(2Ut2),2ztLUzt=ρVv2eff2zt(2Uzt),
    (10)

    其中

    k1=(c11c212c11),k2=2(c13c12c13c11),k3=(c13e31c11e33)Ez,k4=(c13c11λ11λ33)Θ,k5=(c33c213c11).
    (11)

    令位移梯度u=Uz, 由方程(9)可得出压电圆杆波导的纵波运动方程

    2ut2c022uz2=2z2(αu2+β2ut2),
    (12)
    c20=k1v2eff+k3+12k4+k5ρ,α=3k2veff+6k54ρ,β=v2effr2,
    (13)

    其中, c0是压电圆杆的线性纵向波速; α为耗散系数, β为弥散系数, 两者都由材料的性质和几何参数决定. 值得说明的是, 如果忽略电场、热-电耦合, 考虑纯弹性杆, 结果与文献[16]一致.

    在推导方程(12)时, 使用了有效泊松比veff[28]. 可以得到

    veff=εr/εz=ε33(c12c13+c11c13)e31e33(c11+c12)2c212ε332c11e231+c12e231c211ε33.
    (14)

    假设方程(12)的行波解为

    u=u(ξ),ξ=k(zct),
    (15)

    将方程(15)代入方程(12)可以转化为一个常微分方程:

    βc2k2d4udξ4(c2c20)d2udξ2+αd2u2dξ2=0,
    (16)

    其中kc分别是波数和波速. 用下列Jacobi椭圆函数表示方程(16)的解:

    u(ξ)=nj=0ajcnjξ.
    (17)

    可以看出u(ξ)的最高阶数是n, 即

    O(u(ξ))=n.
    (18)

    为了讨论方便, 列出三种Jacobi椭圆函数之间的关系以及渐近值

    sn2ξ+cn2ξ=1,dn2ξ+m2sn2ξ=1,
    (19)

    其中cnξ为Jacobi椭圆余弦函数, dnξ为第三类Jacobi椭圆函数, m为模数(0m1).

    ddξsnξ=cnξdnξ,ddξcnξ=snξdnξ,ddξdnξ=m2snξdnξ.
    (20)

    根据上述微分关系, 很容易有

    O(dudξ)=n+1.
    (21)

    类似地, 可以推出:

    O(u2)=2n,O(d2udξ2)=n+2.
    (22)

    因此, 将方程(16)进一步化为

    k2d4udξ4N1d2udξ2+N2d2u2dξ2=0,
    (23)

    其中

    N1=c2c20βc2,N2=αβc2.
    (24)

    ξ积分两次, 为计算方便, 令积分常数为零. 可以得到

    k2uξξN1u+N2(u2)=0.
    (25)

    通过谐波平衡法, 使方程(25)中的非线性项次数和微分项最高阶数相等, 结合方程(22)可以确定方程(18)中的最高阶数n = 2.

    根据Jacobi椭圆余弦函数展开法, 方程(17)的解有如下表达形式:

    u(ξ)=a0+a1cnξ+a2cn2ξ.
    (26)

    方程(26)对ξ微分两次可求得

    d2udξ2=2a2(1m2)a1(12m2)cnξ+4a2(2m21)cn2ξ2a1m2cn3ξ6a2m2cn4ξ.
    (27)

    将方程(26)和方程(27)代入方程(25), 并比较cnξ相同次幂的系数, 可得

    a0=N14k2(2m21)2N2,a1=0,a2=6k2m2N2.
    (28)

    因此方程(23)的精确周期解为

    u(ξ)=N14k2(2m21)2N2+6k2m2N2cn2(ξ,m),
    (29)

    其中m为模数(0 < m < 1). 应该指出, 当m1时, cnξsechξ, 于是非线性波动方程(29)的孤波解可以写成:

    u(ξ)=N14k22N2+6k2N2sech2ξ.
    (30)

    由于此处只讨论波的特性, 令常数项为零. 即

    k2=c2c204βc2.
    (31)

    将方程(31)代入方程(30), 得到方程(25)的一个标准孤波解:

    u(ξ)=Asech2zctΛ,
    (32)

    式中A为波幅, Λ是波长,

    A=3(c2c20)2α,Λ=2πk=4πc2βc2c20,
    (33)

    其中c > c0是孤立波存在的条件.

    利用Jacobi椭圆正弦函数展开法, 方程(17)的解有如下表达形式:

    u(ξ)=a0+a1snξ+a2sn2ξ.
    (34)

    将方程(34)对ξ微分两次求得

    d2udξ2=2a2a1(1+2m2)snξ4a2(1+m2)sn2ξ+2a2m2sn3ξ+6a2m2sn4ξ.
    (35)

    将方程(34)和方程(35)代入方程(25), 并比较snξ相同次幂的系数, 可得

    a0=4k2(1+m2)+N12N2,a1=0,a2=6k2m2N2.
    (36)

    因此方程(23)的精确周期解为

    u(ξ)=4k2(1+m2)+N12N2+6k2m2N2sn2(ξ,m).
    (37)

    将方程(19)代入到方程(37)求得

    u(ξ)=4k2(12m2)+N12N2+6k2m2N2cn2(ξ,m),
    (38)

    这就是方程(23)的另一个精确周期解, 可以观察到取m = 1时, 方程(38)退化为与方程(29)同样的形式.

    本文通过Matlab软件, 对如图1所示的模型进行数值模拟, 取无限长压电圆杆为BaTiO3材料, 初始温度为T0 = 20 ℃, 杆的半径分别为R1 = 0.025 m, R2 = 0.05 m和R3 = 0.075 m. 将这些参数代入相关方程进行计算, 通过改变温度及波速比值, 得到相应的模拟结果.

    根据表1所列参数, 通过改变温差Θ的大小, 可以计算出不同温度下波速c0的大小. 如表2所列, 可以看出, 当温差Θ = 10 ℃时, 波速c0最大, 随着温度的升高, 波速c0逐渐减小, 但也可以看出波速c0并没有存在大幅度的衰减.

    表 1  钛酸钡材料参数[29]
    Table 1.  Material parameters of BaTiO3[29].
    参数参数值参数参数值
    c11/(N·m–2)166 × 109 e15/(C·m–2)11.6
    c12/(N·m–2)77 × 109 ε11/(C2·N–2·m–1)11.2 × 10–9
    c13/(N·m–2)78 × 109ε33/(C2·N-2·m–1)12.9 × 10–9
    c33/(N·m–2)162 × 109d1/(C·km–2)–5.4831 × 10–6
    c44/(N·m–2)43 × 109d3/(C·km–2)–5.4831 × 10–6
    e31/(C·m–2)–4.4α1 = 0.5α310 × 10–6
    e33/(C·m–2)18.6β1/m20.5278 × 10–4
    β2/m22.112 × 10–4β3/m24.7502 × 10–4
    α/(N·m·kg–1)3.5561 × 107veff0.2906
    ρ/(kg·m-3)5.8 × 103
    下载: 导出CSV 
    | 显示表格
    表 2  不同温度下波速比较
    Table 2.  Comparison of the wave velocities at different temperature.
    Θ/℃c0/(103 m·s–1)
    105.4482
    505.4367
    905.4251
    下载: 导出CSV 
    | 显示表格

    通过改变波速比c/c0的大小, 可求得在不同比值下的波幅A、波长Λ和波数k. 如表3所列, 可以看出, 随着波速比的增大, 波幅A和波数k增大, 而波长Λ随之减小. 同样地, 如表4所列, 固定波速比c/c0 = 1.1, 改变压电圆杆的半径, 分别取R1 = 0.025 m, R2 = 0.05 m和R3 = 0.075 m. 能够发现随着压电圆杆半径R的增大, 波长Λ随之增大, 而波数k随之减小.

    表 3  R = 0.05 m时不同波速比下参数比较
    Table 3.  Comparison of parameters under different wave velocity ratios when R = 0.05 m.
    c/c0A/mΛ/mk
    1.1 0.2567 0.5429 11.5675
    1.2 0.5378 0.4092 15.3470
    1.3 0.8434 0.3540 17.740
    下载: 导出CSV 
    | 显示表格
    表 4  波速比c/c0 = 1.1时不同半径下参数比较
    Table 4.  Comparison of parameters under different radii when c/c0 = 1.1.
    R/mΛ/mk
    0.0250.270523.2162
    0.0500.542911.5675
    0.0750.81167.7375
    下载: 导出CSV 
    | 显示表格

    图2给出了在不考虑温度影响时不同波速比c/c0下, 孤波波速u和变量ξ的关系. 可以看出, 当ξ = 0时, u达到最大值, 并且幅值关于ξ = 0对称, 同时可以观察到, 随着c/c0的增大, 孤波幅值增大, 波长减小. 简而言之, 孤波振幅越大, 波长越小, 这体现了非线性孤波的弥散特性.

    图 2 不同波速比c/c0下孤波波速u与变量ξ的关系\r\nFig. 2. Relationship between solitary wave u and variable ξ under different wave velocity ratio c/c0 values.
    图 2  不同波速比c/c0下孤波波速u与变量ξ的关系
    Fig. 2.  Relationship between solitary wave u and variable ξ under different wave velocity ratio c/c0 values.

    图3给出了在波速比c/c0 = 1.3, 时间固定在t = 0.01 s时, 当Θ分别取10, 50, 90 ℃三种不同温度下的波形[30], 观察到随着温度Θ的改变, 孤立波在传播过程中波形并没有发生改变, 体现了其稳定性. 同时可以看出, 当波速比c/c0一定时它们的波幅非常接近, 温度的改变对波幅的影响并不是很明显, 但是随着温度的逐渐升高, 波速却逐渐降低, 这一点也与表2的数据相符合.

    图 3 波速比c/c0 = 1.3时三种不同温度下的波形\r\nFig. 3. Three waveforms at different temperatures when the velocity ratio of c/c0 = 1.3.
    图 3  波速比c/c0 = 1.3时三种不同温度下的波形
    Fig. 3.  Three waveforms at different temperatures when the velocity ratio of c/c0 = 1.3.

    图4给出了当温差Θ = 50 ℃时, 波速比c/c0分别取为1.1, 1.2和1.3时的波形, 可知, 孤立波的能量主要集中在中间有效区域, 并且没有因为波速比的改变而扩散, 体现了孤立波的稳定性. 而当温度一定时, 随着波速比c/c0的升高, 波幅增大, 波长减小. 这一点与图2表示的孤波特性相一致.

    图 4 当Θ = 50 ℃时不同波速比c/c0下的波形\r\nFig. 4. Waveforms under different wave velocity ratio c/c0 values when Θ = 50 ℃.
    图 4  Θ = 50 ℃时不同波速比c/c0下的波形
    Fig. 4.  Waveforms under different wave velocity ratio c/c0 values when Θ = 50 ℃.

    图5给出了孤立波的三维曲面, 图6给出了压电圆杆的孤波特性. 可以看出, 当给定某一时间tz之后, 孤立子就会出现, 这说明了孤立波不是单独关于时间或空间的单一变量, 而是以时间和空间为组合的变量, 随着时间的变化, 波在传播. 并且相同的波形在给定tz的组合后会重复出现, 且能量比较集中, 这与孤波的稳定性相符合. 因此它在当代通信技术、缺陷检测等许多方面有很大的研究意义和应用潜力.

    图 5 当波速比c/c0 = 1.3, Θ = 50 ℃时三维曲面图\r\nFig. 5. Three-dimensional surface figure when the wave ratio c/c0 = 1.3, Θ = 50 ℃.
    图 5  当波速比c/c0 = 1.3, Θ = 50 ℃时三维曲面图
    Fig. 5.  Three-dimensional surface figure when the wave ratio c/c0 = 1.3, Θ = 50 ℃.
    图 6 压电圆杆的孤波特性\r\nFig. 6. Solitariness of piezoelectric rod.
    图 6  压电圆杆的孤波特性
    Fig. 6.  Solitariness of piezoelectric rod.

    图7给出了在三种不同温度下, 压电圆杆的波速c和波数k的关系图. 可以看出, 当温度一定时, 随着波数k的增加, 波速c也呈增加趋势, 同样地, 当波数k一定时, 随着温度的升高波速c反而呈减小趋势, 这与表2以及图3图4的模拟结果相符合.

    图 7 三种不同温度下波速c和波数k的关系图\r\nFig. 7. Graph of wave velocity c and wave number k at three different temperatures.
    图 7  三种不同温度下波速c和波数k的关系图
    Fig. 7.  Graph of wave velocity c and wave number k at three different temperatures.

    孤立波在数学上是一类非线性偏微分方程的局部行波解, 这类方程可表示为运动项+色散项+非线性项+(耗散项) = 0. 从物理本质上讲, 非线性效应使得波形在传播过程中出现陡突(能量聚集). 孤立子就是由非线性场激发的、能量不弥散的、形态上稳定的粒子. 本文研究了温度效应下无限长压电圆杆的孤波问题, 由于结构的有限变形(如轴向压缩等)可引起非线性效应, 而二次运动和变形(如横向泊松效应等)可分散这些效应. 相互作用在非线性效应和色散效应之间, 在一定条件下就产生了孤波这一稳定传播的行波. 利用Hamilton变分及Euler方程推导出圆杆的波动方程, 并采用Jacobi椭圆余弦函数展开法和椭圆正弦函数展开法对推导出的波动方程进行求解. 通过Matlab软件进行数值模拟, 发现压电圆杆中不仅有孤波存在, 而且在不同温度下具有不同的性质. 可以得到如下结论:

    1)由计算结果可分析得到, 当固定波速比时, 温度的改变对波速的影响比较明显, 随着温度的升高, 波速逐渐降低.

    2)当固定温度时, 波速比的改变对孤波的幅值影响比较明显, 随着波速比的增大, 波幅逐渐升高, 这也是孤波的特点之一.

    3)温度的改变虽然对孤立波有一定的影响, 但在传播的过程中, 孤立波仍然是关于ξ对称的钟型波, 这也体现了非线性和色散效应共同作用下孤立波的稳定特性. 另外, 利用Jacobi椭圆函数求解得到压电圆杆的精确周期解, 周期解可退化为孤波解, 从理论上也证明了压电圆杆中可能有稳定传播的孤波. 目前, 波动理论在结构的无损检测和提高信息传输质量等方面得到了较为广泛的应用. 因此, 将波动理论运用到压电材料中具有现实工程意义和理论研究价值.

    [1]

    Janshoff A, Steinem C, Galla H J 2000 Angew. Chem. Int. Ed. 39 4004Google Scholar

    [2]

    刘延柱, 薛纭, 陈立群 2004 物理学报 53 2424Google Scholar

    Liu Y Z, Xue Y, Chen L Q 2004 Acta Phys. Sin. 53 2424Google Scholar

    [3]

    He T H, Tian X G, Shen Y P 2002 Int. J. Eng. Sci. 40 1081Google Scholar

    [4]

    范恩贵 2000 物理学报 49 1409Google Scholar

    Fan E G 2000 Acta Phys. Sin. 49 1409Google Scholar

    [5]

    李向正, 张卫国, 原三领 2010 物理学报 59 744Google Scholar

    Li X Z, Zhang W G, Yuan S L 2010 Acta Phys. Sin. 59 744Google Scholar

    [6]

    冯依虎 2019 应用力学与数学 40 1Google Scholar

    Feng Y H 2019 Appl. Math. Mech. 40 1Google Scholar

    [7]

    Guo J G, Zhou L J, Zhang S Y 2005 Appl. Math. Mech. 26 667Google Scholar

    [8]

    李敏, 王博婷, 许韬, 水涓涓 2020 物理学报 69 010502Google Scholar

    Li M, Wang B T, Xu T, Shui J J 2020 Acta Phys. Sin. 69 010502Google Scholar

    [9]

    李志斌, 潘素起 2004 物理学报 50 402Google Scholar

    Li Z B, Pan S Q 2004 Acta Phys. Sin. 50 402Google Scholar

    [10]

    Xia T C, Li B, Zhang H Q 2001 Appl. Math. E-Notes 1 139

    [11]

    钱存, 王亮亮, 张解放 2011 物理学报 60 064214Google Scholar

    Qian C, Wang L L, Zhang J F 2011 Acta Phys. Sin. 60 064214Google Scholar

    [12]

    刘式达, 傅遵涛, 刘式适, 赵强 2002 物理学报 51 718Google Scholar

    Liu S D, Fu Z T, Liu S K, Zhao Q 2002 Acta Phys. Sin. 51 718Google Scholar

    [13]

    Zhang S Y, Zhang W 1987 Acta Mech. Sin. 3 64Google Scholar

    [14]

    刘志芳, 张善元 2007 固体力学学报 28 55Google Scholar

    Liu Z F, Zhang S Y 2007 Acta Mech. Solid. Sin. 28 55Google Scholar

    [15]

    刘志芳, 张善元 2006 物理学报 55 628Google Scholar

    Liu Z F, Zhang S Y 2006 Acta Phys. Sin. 55 628Google Scholar

    [16]

    Liu Z F, Zhang S Y 2006 Appl. Math. Mech. 27 1431Google Scholar

    [17]

    邓庆田, 罗松南, 彭亮 2009 应用力学学报 26 519Google Scholar

    Deng Q T, Luo S N, Peng L 2009 Chin. J. Appl. Mech. 26 519Google Scholar

    [18]

    Seadawy A R, Manafian J 2018 Results Phys. 8 1158Google Scholar

    [19]

    Baskonus H M, Bulut H, Atangana A 2016 Smart Mater. Struct. 25 035022Google Scholar

    [20]

    Bulut H, Sulaiman T A, Baskonus H M 2018 Opt. Quantum. Electron. 50 2Google Scholar

    [21]

    Wang Q 2001 Int. J. Solids. Struct. 38 8207Google Scholar

    [22]

    Xue C X, Pan E 2013 Int. J. Eng. Sci. 62 48Google Scholar

    [23]

    Samsonov A M 2001 Strain Solitons in Solids and How to Construct Them (New York: Chapman and Hall/CRC) p111

    [24]

    Toffoli T, Fernandez L, Monbaliu J, Benoit M, Gagnaire-Renou E, Lefèvre J M, Cavaleri L, Proment D, Pakozdi C, Stansberg C T, Waseda T, Onorato M 2013 Phys. Fluids 25 091701Google Scholar

    [25]

    Ansari R, Oskouie F M, Gholami R, Sadeghi F 2016 Compos. Part B-Eng. 89 316Google Scholar

    [26]

    Ootao Y, Tanigawa Y 2000 Int. J. Eng. Sci. 38 47Google Scholar

    [27]

    Shiv P J 1992 Smart Mater. Struct. 1 80Google Scholar

    [28]

    Xue C X, Pan E, Zhang S Y 2011 Smart Mater. Struct. 20 105010Google Scholar

    [29]

    常瑞鼎 2014 硕士学位论文 (湘潭: 湘潭大学)

    Chang R D 2014 M. S. Thesis (Hunan: Xiangtan University) (in Chinese)

    [30]

    贾菲B 著 (林声和 译) 1976 压电陶瓷 (北京: 科学出版社)第125−145页

    Jaffe B (translated by Lin S H) 1976 Piezoelectric Ceramics (Beijing: Science Press) pp125−145 (in Chinese)

  • 图 1  压电圆杆示意图

    Fig. 1.  Schematic diagram of piezoelectric rod.

    图 2  不同波速比c/c0下孤波波速u与变量ξ的关系

    Fig. 2.  Relationship between solitary wave u and variable ξ under different wave velocity ratio c/c0 values.

    图 3  波速比c/c0 = 1.3时三种不同温度下的波形

    Fig. 3.  Three waveforms at different temperatures when the velocity ratio of c/c0 = 1.3.

    图 4  Θ = 50 ℃时不同波速比c/c0下的波形

    Fig. 4.  Waveforms under different wave velocity ratio c/c0 values when Θ = 50 ℃.

    图 5  当波速比c/c0 = 1.3, Θ = 50 ℃时三维曲面图

    Fig. 5.  Three-dimensional surface figure when the wave ratio c/c0 = 1.3, Θ = 50 ℃.

    图 6  压电圆杆的孤波特性

    Fig. 6.  Solitariness of piezoelectric rod.

    图 7  三种不同温度下波速c和波数k的关系图

    Fig. 7.  Graph of wave velocity c and wave number k at three different temperatures.

    表 1  钛酸钡材料参数[29]

    Table 1.  Material parameters of BaTiO3[29].

    参数参数值参数参数值
    c11/(N·m–2)166 × 109 e15/(C·m–2)11.6
    c12/(N·m–2)77 × 109 ε11/(C2·N–2·m–1)11.2 × 10–9
    c13/(N·m–2)78 × 109ε33/(C2·N-2·m–1)12.9 × 10–9
    c33/(N·m–2)162 × 109d1/(C·km–2)–5.4831 × 10–6
    c44/(N·m–2)43 × 109d3/(C·km–2)–5.4831 × 10–6
    e31/(C·m–2)–4.4α1 = 0.5α310 × 10–6
    e33/(C·m–2)18.6β1/m20.5278 × 10–4
    β2/m22.112 × 10–4β3/m24.7502 × 10–4
    α/(N·m·kg–1)3.5561 × 107veff0.2906
    ρ/(kg·m-3)5.8 × 103
    下载: 导出CSV

    表 2  不同温度下波速比较

    Table 2.  Comparison of the wave velocities at different temperature.

    Θ/℃c0/(103 m·s–1)
    105.4482
    505.4367
    905.4251
    下载: 导出CSV

    表 3  R = 0.05 m时不同波速比下参数比较

    Table 3.  Comparison of parameters under different wave velocity ratios when R = 0.05 m.

    c/c0A/mΛ/mk
    1.1 0.2567 0.5429 11.5675
    1.2 0.5378 0.4092 15.3470
    1.3 0.8434 0.3540 17.740
    下载: 导出CSV

    表 4  波速比c/c0 = 1.1时不同半径下参数比较

    Table 4.  Comparison of parameters under different radii when c/c0 = 1.1.

    R/mΛ/mk
    0.0250.270523.2162
    0.0500.542911.5675
    0.0750.81167.7375
    下载: 导出CSV
  • [1]

    Janshoff A, Steinem C, Galla H J 2000 Angew. Chem. Int. Ed. 39 4004Google Scholar

    [2]

    刘延柱, 薛纭, 陈立群 2004 物理学报 53 2424Google Scholar

    Liu Y Z, Xue Y, Chen L Q 2004 Acta Phys. Sin. 53 2424Google Scholar

    [3]

    He T H, Tian X G, Shen Y P 2002 Int. J. Eng. Sci. 40 1081Google Scholar

    [4]

    范恩贵 2000 物理学报 49 1409Google Scholar

    Fan E G 2000 Acta Phys. Sin. 49 1409Google Scholar

    [5]

    李向正, 张卫国, 原三领 2010 物理学报 59 744Google Scholar

    Li X Z, Zhang W G, Yuan S L 2010 Acta Phys. Sin. 59 744Google Scholar

    [6]

    冯依虎 2019 应用力学与数学 40 1Google Scholar

    Feng Y H 2019 Appl. Math. Mech. 40 1Google Scholar

    [7]

    Guo J G, Zhou L J, Zhang S Y 2005 Appl. Math. Mech. 26 667Google Scholar

    [8]

    李敏, 王博婷, 许韬, 水涓涓 2020 物理学报 69 010502Google Scholar

    Li M, Wang B T, Xu T, Shui J J 2020 Acta Phys. Sin. 69 010502Google Scholar

    [9]

    李志斌, 潘素起 2004 物理学报 50 402Google Scholar

    Li Z B, Pan S Q 2004 Acta Phys. Sin. 50 402Google Scholar

    [10]

    Xia T C, Li B, Zhang H Q 2001 Appl. Math. E-Notes 1 139

    [11]

    钱存, 王亮亮, 张解放 2011 物理学报 60 064214Google Scholar

    Qian C, Wang L L, Zhang J F 2011 Acta Phys. Sin. 60 064214Google Scholar

    [12]

    刘式达, 傅遵涛, 刘式适, 赵强 2002 物理学报 51 718Google Scholar

    Liu S D, Fu Z T, Liu S K, Zhao Q 2002 Acta Phys. Sin. 51 718Google Scholar

    [13]

    Zhang S Y, Zhang W 1987 Acta Mech. Sin. 3 64Google Scholar

    [14]

    刘志芳, 张善元 2007 固体力学学报 28 55Google Scholar

    Liu Z F, Zhang S Y 2007 Acta Mech. Solid. Sin. 28 55Google Scholar

    [15]

    刘志芳, 张善元 2006 物理学报 55 628Google Scholar

    Liu Z F, Zhang S Y 2006 Acta Phys. Sin. 55 628Google Scholar

    [16]

    Liu Z F, Zhang S Y 2006 Appl. Math. Mech. 27 1431Google Scholar

    [17]

    邓庆田, 罗松南, 彭亮 2009 应用力学学报 26 519Google Scholar

    Deng Q T, Luo S N, Peng L 2009 Chin. J. Appl. Mech. 26 519Google Scholar

    [18]

    Seadawy A R, Manafian J 2018 Results Phys. 8 1158Google Scholar

    [19]

    Baskonus H M, Bulut H, Atangana A 2016 Smart Mater. Struct. 25 035022Google Scholar

    [20]

    Bulut H, Sulaiman T A, Baskonus H M 2018 Opt. Quantum. Electron. 50 2Google Scholar

    [21]

    Wang Q 2001 Int. J. Solids. Struct. 38 8207Google Scholar

    [22]

    Xue C X, Pan E 2013 Int. J. Eng. Sci. 62 48Google Scholar

    [23]

    Samsonov A M 2001 Strain Solitons in Solids and How to Construct Them (New York: Chapman and Hall/CRC) p111

    [24]

    Toffoli T, Fernandez L, Monbaliu J, Benoit M, Gagnaire-Renou E, Lefèvre J M, Cavaleri L, Proment D, Pakozdi C, Stansberg C T, Waseda T, Onorato M 2013 Phys. Fluids 25 091701Google Scholar

    [25]

    Ansari R, Oskouie F M, Gholami R, Sadeghi F 2016 Compos. Part B-Eng. 89 316Google Scholar

    [26]

    Ootao Y, Tanigawa Y 2000 Int. J. Eng. Sci. 38 47Google Scholar

    [27]

    Shiv P J 1992 Smart Mater. Struct. 1 80Google Scholar

    [28]

    Xue C X, Pan E, Zhang S Y 2011 Smart Mater. Struct. 20 105010Google Scholar

    [29]

    常瑞鼎 2014 硕士学位论文 (湘潭: 湘潭大学)

    Chang R D 2014 M. S. Thesis (Hunan: Xiangtan University) (in Chinese)

    [30]

    贾菲B 著 (林声和 译) 1976 压电陶瓷 (北京: 科学出版社)第125−145页

    Jaffe B (translated by Lin S H) 1976 Piezoelectric Ceramics (Beijing: Science Press) pp125−145 (in Chinese)

  • [1] 吕杰, 方贺男, 吕涛涛, 孙星宇. MgO基磁性隧道结温度-偏压相图的理论研究. 物理学报, 2021, 70(10): 107302. doi: 10.7498/aps.70.20201905
    [2] 王云天, 曾祥国, 杨鑫. 高应变率下温度对单晶铁中孔洞成核与生长影响的分子动力学研究. 物理学报, 2019, 68(24): 246102. doi: 10.7498/aps.68.20190920
    [3] 金鑫, 杨春明, 滑文强, 李怡雯, 王劼. PS3000-b-PAA5000球形胶束温度效应的原位小角X射线散射技术研究. 物理学报, 2018, 67(4): 048301. doi: 10.7498/aps.67.20172167
    [4] 李明林, 万亚玲, 胡建玥, 王卫东. 单层二硫化钼力学性能温度和手性效应的分子动力学模拟. 物理学报, 2016, 65(17): 176201. doi: 10.7498/aps.65.176201
    [5] 杨永锋, 冯海波, 陈虎, 仵敏娟. 柔性杆与凸轮斜碰撞特性分析. 物理学报, 2016, 65(24): 240502. doi: 10.7498/aps.65.240502
    [6] 胡雪兰, 罗阳, 赵若汐, 胡艳敏, 张艳峰, 宋庆功. NiAl中Ni空位对杂质C原子的多重俘获及温度效应的第一性原理研究. 物理学报, 2016, 65(20): 206101. doi: 10.7498/aps.65.206101
    [7] 汪维刚, 林万涛, 石兰芳, 莫嘉琪. 非线性扰动时滞长波系统孤波近似解. 物理学报, 2014, 63(11): 110204. doi: 10.7498/aps.63.110204
    [8] 许永红, 韩祥临, 石兰芳, 莫嘉琪. 薛定谔扰动耦合系统孤波的行波近似解法. 物理学报, 2014, 63(9): 090204. doi: 10.7498/aps.63.090204
    [9] 郭巧能, 曹义刚, 孙强, 刘忠侠, 贾瑜, 霍裕平. 温度对超薄铜膜疲劳性能影响的分子动力学模拟. 物理学报, 2013, 62(10): 107103. doi: 10.7498/aps.62.107103
    [10] 欧阳成, 石兰芳, 林万涛, 莫嘉琪. (2+1)维扰动时滞破裂孤波方程行波解的摄动方法. 物理学报, 2013, 62(17): 170201. doi: 10.7498/aps.62.170201
    [11] 强蕾, 姚若河. 非晶硅薄膜晶体管沟道中阈值电压及温度的分布. 物理学报, 2012, 61(8): 087303. doi: 10.7498/aps.61.087303
    [12] 周先春, 林万涛, 林一骅, 莫嘉琪. 大气非均匀量子等离子体孤波解. 物理学报, 2012, 61(24): 240202. doi: 10.7498/aps.61.240202
    [13] 陈琼, 杨先清, 赵新印, 王振辉, 赵跃民. 周期型二元颗粒链中孤波传播的二体碰撞近似分析. 物理学报, 2012, 61(4): 044501. doi: 10.7498/aps.61.044501
    [14] 许永红, 姚静荪, 莫嘉琪. (3+1)维Burgers扰动系统孤波的解法. 物理学报, 2012, 61(2): 020202. doi: 10.7498/aps.61.020202
    [15] 宋柏, 吴晶, 过增元. 基于热质理论的Hamilton原理. 物理学报, 2010, 59(10): 7129-7134. doi: 10.7498/aps.59.7129
    [16] 吴亚敏, 陈国庆. 带壳颗粒复合介质光学双稳的温度效应. 物理学报, 2009, 58(3): 2056-2060. doi: 10.7498/aps.58.2056
    [17] 莫嘉琪, 张伟江, 陈贤峰. 一类强非线性发展方程孤波变分迭代解法. 物理学报, 2009, 58(11): 7397-7401. doi: 10.7498/aps.58.7397
    [18] 段文山, 洪学仁. 弱相对论等离子体横向扰动下的离子声孤波. 物理学报, 2003, 52(6): 1337-1339. doi: 10.7498/aps.52.1337
    [19] 徐桂琼, 李志斌. 构造非线性发展方程孤波解的混合指数方法. 物理学报, 2002, 51(5): 946-950. doi: 10.7498/aps.51.946
    [20] 李志斌, 潘素起. 广义五阶KdV方程的孤波解与孤子解. 物理学报, 2001, 50(3): 402-405. doi: 10.7498/aps.50.402
计量
  • 文章访问数:  5575
  • PDF下载量:  85
出版历程
  • 收稿日期:  2020-05-21
  • 修回日期:  2020-09-14
  • 上网日期:  2021-01-18
  • 刊出日期:  2021-02-05

/

返回文章
返回