-
利用有限变形理论, 以无限长压电圆杆为研究对象, 考虑了在横向惯性、等效泊松比效应以及在热电弹耦合共同作用下, 基于Hamilton原理, 并引入Euler方程推导出压电圆杆的纵向波动方程. 采用Jacobi椭圆函数展开法, 求解压电圆杆的波动方程和对应的解. 最后, 通过Matlab软件得到不同波速比下的色散曲线, 以及温度场对压电圆杆的波形、波幅和波数的影响曲线. 数值分析结果表明: 随着温度的升高, 波速逐渐降低, 温度场的改变可影响和控制孤立波的传播特性.
-
关键词:
- 压电圆杆 /
- Hamilton原理 /
- 温度效应 /
- 孤波
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.-
Keywords:
- piezoelectric rod /
- Hamilton principle /
- temperature effect /
- solitary wave
1. 引 言
随着现代科学技术的飞速发展, 压电材料因其独特的性能而被广泛应用[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椭圆函数展开法, 推导出压电圆杆的波动方程和对应的解, 并讨论温度变化对压电介质的波形、波幅以及波速等的影响.
2. 建立压电圆杆的波动方程
图1为无限长压电圆杆示意图, 建立圆柱坐标系
(r,θ,z) , z是沿着杆的轴向, 也是波传播的方向. θ = [0, 2π], 0 ≤ r ≤ R, 其中R为压电圆杆半径. 为了研究方便, 假设: 1)变形后, 板中初始状态垂直于中心平面的点仍然垂直于中心平面; 2)杆的截面是轴对称的, 即Uθ=0 和∂/∂θ=0 , 其中Uθ 为θ方向位移; 3)考虑泊松比效应, 纵向位移U和径向位移Ur 之间满足Ur=veffr∂U/∂z , 其中veff 是有效泊松比.对于横观各向同性的压电材料圆杆, 在考虑温度效应时的本构方程如下[25]:
{σr=c11εr+c12εθ+c13εz−e31Ez−λ11Θ,σθ=c12εr+c11εθ+c13εz−e31Ez−λ11Θ,σz=c13εr+c13εθ+c33εz−e33Ez−λ33Θ,τrz=c44γrz−e15Er,τθz=c44γθz−e15Eθ,τ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=∂Ur∂r,εθ=∂Uθr∂θ+Urr,εz=∂U∂z+12(∂U∂z)2,γrθ=∂Urr∂θ+∂Uθ∂r−Uθr,γθz=∂Ur∂θ+∂Uθ∂z,γrz=∂U∂r+∂Ur∂z. (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变分原理可得
δ∫tt0∫zz0Ldzdt=δ∫tt0∫zz0(T−EP+We)dzdt=0, (5) 式中L为Lagrange密度函数, T为系统的动能, EP为系统的势能, We为系统的电能, 具体表达式分别为
T=12∫VρdV(∂U∂t)2+12∫VρdV(∂Ur∂t)2, (6) EP=12∫VSTQdV, (7) We=12∫VETDdV, (8) 这里
ρ 为压电材料的密度, V为压电材料的体积,S={εr,εθ,γrθ}T 和Q={σθ,σz,τθz}T 分别表示材料的应变向量和应力向量,E={0,0,Ez}T 表示电场向量,D={0,0,Dz}T 表示电位移向量.根据Euler方程, 若
L=L(U,Uz,Ut,Uzz,Utt,Uzt,⋅⋅⋅) , 则∂L∂U−∂∂z∂L∂Uz−∂∂t∂L∂Ut+∂2∂z2∂L∂Uzz+∂2∂t2∂L∂Utt+∂2∂z∂t∂L∂Uzt−⋅⋅⋅=0, (9) 此处,
Uz=∂U∂z,Ut=∂U∂t,Uzt=∂2U∂z∂t⋯ , 其余类同.根据以上关系可以得到如下表达式:
{∂L∂U=0,∂∂z∂L∂Uz=−12V[(2k1v2eff+2k2v2eff+2k3+k4+k5)(∂U∂z)+3(12k2veff+k5)(∂U∂z)2+k5(∂U∂z)3],∂∂t∂L∂Ut=ρV∂2U∂t2, ∂2∂z2∂L∂Uzz=0,∂2∂t2∂L∂Utt=ρV∂2∂t2(∂2U∂t2),∂2∂z∂t∂L∂Uzt=ρVv2eff∂2∂z∂t(∂2U∂z∂t), (10) 其中
k1=(c11−c212c11),k2=2(c13−c12c13c11),k3=(c13e31c11−e33)Ez,k4=(c13c11λ11−λ33)Θ,k5=(c33−c213c11). (11) 令位移梯度
u=∂U∂z , 由方程(9)可得出压电圆杆波导的纵波运动方程∂2u∂t2−c02∂2u∂z2=∂2∂z2(αu2+β∂2u∂t2), (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ε33−2c11e231+c12e231−c211ε33. (14) 3. 纵波方程的孤波解
假设方程(12)的行波解为
u=u(ξ),ξ=k(z−ct), (15) βc2k2d4udξ4−(c2−c20)d2udξ2+αd2u2dξ2=0, (16) 其中k和c分别是波数和波速. 用下列Jacobi椭圆函数表示方程(16)的解:
u(ξ)=n∑j=0ajcnjξ. (17) 可以看出
u(ξ) 的最高阶数是n, 即O(u(ξ))=n. (18) 为了讨论方便, 列出三种Jacobi椭圆函数之间的关系以及渐近值
sn2ξ+cn2ξ=1,dn2ξ+m2sn2ξ=1, (19) 其中
cnξ 为Jacobi椭圆余弦函数,dnξ 为第三类Jacobi椭圆函数, m为模数(0⩽m⩽1) .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ξ4−N1d2udξ2+N2d2u2dξ2=0, (23) 其中
N1=c2−c20β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(1−m2)−a1(1−2m2)cnξ+4a2(2m2−1)cn2ξ−2a1m2cn3ξ−6a2m2cn4ξ. (27) 将方程(26)和方程(27)代入方程(25), 并比较
cnξ 相同次幂的系数, 可得a0=N1−4k2(2m2−1)2N2,a1=0,a2=6k2m2N2. (28) 因此方程(23)的精确周期解为
u(ξ)=N1−4k2(2m2−1)2N2+6k2m2N2cn2(ξ,m), (29) 其中m为模数(0 < m < 1). 应该指出, 当
m→1 时,cnξ→sechξ , 于是非线性波动方程(29)的孤波解可以写成:u(ξ)=N1−4k22N2+6k2N2sech2ξ. (30) 由于此处只讨论波的特性, 令常数项为零. 即
k2=c2−c204βc2. (31) 将方程(31)代入方程(30), 得到方程(25)的一个标准孤波解:
u(ξ)=Asech2z−ctΛ, (32) 式中A为波幅, Λ是波长,
A=3(c2−c20)2α,Λ=2πk=4π√c2βc2−c20, (33) 其中c > c0是孤立波存在的条件.
利用Jacobi椭圆正弦函数展开法, 方程(17)的解有如下表达形式:
u(ξ)=a0+a1snξ+a2sn2ξ. (34) 将方程(34)对ξ微分两次求得
d2udξ2=2a2−a1(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) u(ξ)=4k2(1−2m2)+N12N2+6k2m2N2cn2(ξ,m), (38) 这就是方程(23)的另一个精确周期解, 可以观察到取m = 1时, 方程(38)退化为与方程(29)同样的形式.
4. 数值计算和讨论
本文通过Matlab软件, 对如图1所示的模型进行数值模拟, 取无限长压电圆杆为BaTiO3材料, 初始温度为T0 = 20 ℃, 杆的半径分别为R1 = 0.025 m, R2 = 0.05 m和R3 = 0.075 m. 将这些参数代入相关方程进行计算, 通过改变温度及波速比值, 得到相应的模拟结果.
根据表1所列参数, 通过改变温差Θ的大小, 可以计算出不同温度下波速c0的大小. 如表2所列, 可以看出, 当温差Θ = 10 ℃时, 波速c0最大, 随着温度的升高, 波速c0逐渐减小, 但也可以看出波速c0并没有存在大幅度的衰减.
参数 参数值 参数 参数值 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 × 109 d1/(C·km–2) –5.4831 × 10–6 c44/(N·m–2) 43 × 109 d3/(C·km–2) –5.4831 × 10–6 e31/(C·m–2) –4.4 α1 = 0.5α3 10 × 10–6 e33/(C·m–2) 18.6 β1/m2 0.5278 × 10–4 β2/m2 2.112 × 10–4 β3/m2 4.7502 × 10–4 α/(N·m·kg–1) 3.5561 × 107 veff 0.2906 ρ/(kg·m-3) 5.8 × 103 表 2 不同温度下波速比较Table 2. Comparison of the wave velocities at different temperature.Θ/℃ c0/(103 m·s–1) 10 5.4482 50 5.4367 90 5.4251 通过改变波速比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/c0 A/m Λ/m k 1.1 0.2567 0.5429 11.5675 1.2 0.5378 0.4092 15.3470 1.3 0.8434 0.3540 17.740 表 4 波速比c/c0 = 1.1时不同半径下参数比较Table 4. Comparison of parameters under different radii when c/c0 = 1.1.R/m Λ/m k 0.025 0.2705 23.2162 0.050 0.5429 11.5675 0.075 0.8116 7.7375 图2给出了在不考虑温度影响时不同波速比c/c0下, 孤波波速u和变量ξ的关系. 可以看出, 当ξ = 0时, u达到最大值, 并且幅值关于ξ = 0对称, 同时可以观察到, 随着c/c0的增大, 孤波幅值增大, 波长减小. 简而言之, 孤波振幅越大, 波长越小, 这体现了非线性孤波的弥散特性.
图3给出了在波速比c/c0 = 1.3, 时间固定在t = 0.01 s时, 当Θ分别取10, 50, 90 ℃三种不同温度下的波形[30], 观察到随着温度Θ的改变, 孤立波在传播过程中波形并没有发生改变, 体现了其稳定性. 同时可以看出, 当波速比c/c0一定时它们的波幅非常接近, 温度的改变对波幅的影响并不是很明显, 但是随着温度的逐渐升高, 波速却逐渐降低, 这一点也与表2的数据相符合.
图4给出了当温差Θ = 50 ℃时, 波速比c/c0分别取为1.1, 1.2和1.3时的波形, 可知, 孤立波的能量主要集中在中间有效区域, 并且没有因为波速比的改变而扩散, 体现了孤立波的稳定性. 而当温度一定时, 随着波速比c/c0的升高, 波幅增大, 波长减小. 这一点与图2表示的孤波特性相一致.
图5给出了孤立波的三维曲面, 图6给出了压电圆杆的孤波特性. 可以看出, 当给定某一时间t和z之后, 孤立子就会出现, 这说明了孤立波不是单独关于时间或空间的单一变量, 而是以时间和空间为组合的变量, 随着时间的变化, 波在传播. 并且相同的波形在给定t和z的组合后会重复出现, 且能量比较集中, 这与孤波的稳定性相符合. 因此它在当代通信技术、缺陷检测等许多方面有很大的研究意义和应用潜力.
图7给出了在三种不同温度下, 压电圆杆的波速c和波数k的关系图. 可以看出, 当温度一定时, 随着波数k的增加, 波速c也呈增加趋势, 同样地, 当波数k一定时, 随着温度的升高波速c反而呈减小趋势, 这与表2以及图3和图4的模拟结果相符合.
5. 结 论
孤立波在数学上是一类非线性偏微分方程的局部行波解, 这类方程可表示为运动项+色散项+非线性项+(耗散项) = 0. 从物理本质上讲, 非线性效应使得波形在传播过程中出现陡突(能量聚集). 孤立子就是由非线性场激发的、能量不弥散的、形态上稳定的粒子. 本文研究了温度效应下无限长压电圆杆的孤波问题, 由于结构的有限变形(如轴向压缩等)可引起非线性效应, 而二次运动和变形(如横向泊松效应等)可分散这些效应. 相互作用在非线性效应和色散效应之间, 在一定条件下就产生了孤波这一稳定传播的行波. 利用Hamilton变分及Euler方程推导出圆杆的波动方程, 并采用Jacobi椭圆余弦函数展开法和椭圆正弦函数展开法对推导出的波动方程进行求解. 通过Matlab软件进行数值模拟, 发现压电圆杆中不仅有孤波存在, 而且在不同温度下具有不同的性质. 可以得到如下结论:
1)由计算结果可分析得到, 当固定波速比时, 温度的改变对波速的影响比较明显, 随着温度的升高, 波速逐渐降低.
2)当固定温度时, 波速比的改变对孤波的幅值影响比较明显, 随着波速比的增大, 波幅逐渐升高, 这也是孤波的特点之一.
3)温度的改变虽然对孤立波有一定的影响, 但在传播的过程中, 孤立波仍然是关于ξ对称的钟型波, 这也体现了非线性和色散效应共同作用下孤立波的稳定特性. 另外, 利用Jacobi椭圆函数求解得到压电圆杆的精确周期解, 周期解可退化为孤波解, 从理论上也证明了压电圆杆中可能有稳定传播的孤波. 目前, 波动理论在结构的无损检测和提高信息传输质量等方面得到了较为广泛的应用. 因此, 将波动理论运用到压电材料中具有现实工程意义和理论研究价值.
[1] Janshoff A, Steinem C, Galla H J 2000 Angew. Chem. Int. Ed. 39 4004
Google Scholar
[2] 刘延柱, 薛纭, 陈立群 2004 物理学报 53 2424
Google Scholar
Liu Y Z, Xue Y, Chen L Q 2004 Acta Phys. Sin. 53 2424
Google Scholar
[3] He T H, Tian X G, Shen Y P 2002 Int. J. Eng. Sci. 40 1081
Google Scholar
[4] 范恩贵 2000 物理学报 49 1409
Google Scholar
Fan E G 2000 Acta Phys. Sin. 49 1409
Google Scholar
[5] 李向正, 张卫国, 原三领 2010 物理学报 59 744
Google Scholar
Li X Z, Zhang W G, Yuan S L 2010 Acta Phys. Sin. 59 744
Google Scholar
[6] 冯依虎 2019 应用力学与数学 40 1
Google Scholar
Feng Y H 2019 Appl. Math. Mech. 40 1
Google Scholar
[7] Guo J G, Zhou L J, Zhang S Y 2005 Appl. Math. Mech. 26 667
Google Scholar
[8] 李敏, 王博婷, 许韬, 水涓涓 2020 物理学报 69 010502
Google Scholar
Li M, Wang B T, Xu T, Shui J J 2020 Acta Phys. Sin. 69 010502
Google Scholar
[9] 李志斌, 潘素起 2004 物理学报 50 402
Google Scholar
Li Z B, Pan S Q 2004 Acta Phys. Sin. 50 402
Google Scholar
[10] Xia T C, Li B, Zhang H Q 2001 Appl. Math. E-Notes 1 139
[11] 钱存, 王亮亮, 张解放 2011 物理学报 60 064214
Google Scholar
Qian C, Wang L L, Zhang J F 2011 Acta Phys. Sin. 60 064214
Google Scholar
[12] 刘式达, 傅遵涛, 刘式适, 赵强 2002 物理学报 51 718
Google Scholar
Liu S D, Fu Z T, Liu S K, Zhao Q 2002 Acta Phys. Sin. 51 718
Google Scholar
[13] Zhang S Y, Zhang W 1987 Acta Mech. Sin. 3 64
Google Scholar
[14] 刘志芳, 张善元 2007 固体力学学报 28 55
Google Scholar
Liu Z F, Zhang S Y 2007 Acta Mech. Solid. Sin. 28 55
Google Scholar
[15] 刘志芳, 张善元 2006 物理学报 55 628
Google Scholar
Liu Z F, Zhang S Y 2006 Acta Phys. Sin. 55 628
Google Scholar
[16] Liu Z F, Zhang S Y 2006 Appl. Math. Mech. 27 1431
Google Scholar
[17] 邓庆田, 罗松南, 彭亮 2009 应用力学学报 26 519
Google Scholar
Deng Q T, Luo S N, Peng L 2009 Chin. J. Appl. Mech. 26 519
Google Scholar
[18] Seadawy A R, Manafian J 2018 Results Phys. 8 1158
Google Scholar
[19] Baskonus H M, Bulut H, Atangana A 2016 Smart Mater. Struct. 25 035022
Google Scholar
[20] Bulut H, Sulaiman T A, Baskonus H M 2018 Opt. Quantum. Electron. 50 2
Google Scholar
[21] Wang Q 2001 Int. J. Solids. Struct. 38 8207
Google Scholar
[22] Xue C X, Pan E 2013 Int. J. Eng. Sci. 62 48
Google 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 091701
Google Scholar
[25] Ansari R, Oskouie F M, Gholami R, Sadeghi F 2016 Compos. Part B-Eng. 89 316
Google Scholar
[26] Ootao Y, Tanigawa Y 2000 Int. J. Eng. Sci. 38 47
Google Scholar
[27] Shiv P J 1992 Smart Mater. Struct. 1 80
Google Scholar
[28] Xue C X, Pan E, Zhang S Y 2011 Smart Mater. Struct. 20 105010
Google 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)
-
参数 参数值 参数 参数值 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 × 109 d1/(C·km–2) –5.4831 × 10–6 c44/(N·m–2) 43 × 109 d3/(C·km–2) –5.4831 × 10–6 e31/(C·m–2) –4.4 α1 = 0.5α3 10 × 10–6 e33/(C·m–2) 18.6 β1/m2 0.5278 × 10–4 β2/m2 2.112 × 10–4 β3/m2 4.7502 × 10–4 α/(N·m·kg–1) 3.5561 × 107 veff 0.2906 ρ/(kg·m-3) 5.8 × 103 表 2 不同温度下波速比较
Table 2. Comparison of the wave velocities at different temperature.
Θ/℃ c0/(103 m·s–1) 10 5.4482 50 5.4367 90 5.4251 表 3 R = 0.05 m时不同波速比下参数比较
Table 3. Comparison of parameters under different wave velocity ratios when R = 0.05 m.
c/c0 A/m Λ/m k 1.1 0.2567 0.5429 11.5675 1.2 0.5378 0.4092 15.3470 1.3 0.8434 0.3540 17.740 表 4 波速比c/c0 = 1.1时不同半径下参数比较
Table 4. Comparison of parameters under different radii when c/c0 = 1.1.
R/m Λ/m k 0.025 0.2705 23.2162 0.050 0.5429 11.5675 0.075 0.8116 7.7375 -
[1] Janshoff A, Steinem C, Galla H J 2000 Angew. Chem. Int. Ed. 39 4004
Google Scholar
[2] 刘延柱, 薛纭, 陈立群 2004 物理学报 53 2424
Google Scholar
Liu Y Z, Xue Y, Chen L Q 2004 Acta Phys. Sin. 53 2424
Google Scholar
[3] He T H, Tian X G, Shen Y P 2002 Int. J. Eng. Sci. 40 1081
Google Scholar
[4] 范恩贵 2000 物理学报 49 1409
Google Scholar
Fan E G 2000 Acta Phys. Sin. 49 1409
Google Scholar
[5] 李向正, 张卫国, 原三领 2010 物理学报 59 744
Google Scholar
Li X Z, Zhang W G, Yuan S L 2010 Acta Phys. Sin. 59 744
Google Scholar
[6] 冯依虎 2019 应用力学与数学 40 1
Google Scholar
Feng Y H 2019 Appl. Math. Mech. 40 1
Google Scholar
[7] Guo J G, Zhou L J, Zhang S Y 2005 Appl. Math. Mech. 26 667
Google Scholar
[8] 李敏, 王博婷, 许韬, 水涓涓 2020 物理学报 69 010502
Google Scholar
Li M, Wang B T, Xu T, Shui J J 2020 Acta Phys. Sin. 69 010502
Google Scholar
[9] 李志斌, 潘素起 2004 物理学报 50 402
Google Scholar
Li Z B, Pan S Q 2004 Acta Phys. Sin. 50 402
Google Scholar
[10] Xia T C, Li B, Zhang H Q 2001 Appl. Math. E-Notes 1 139
[11] 钱存, 王亮亮, 张解放 2011 物理学报 60 064214
Google Scholar
Qian C, Wang L L, Zhang J F 2011 Acta Phys. Sin. 60 064214
Google Scholar
[12] 刘式达, 傅遵涛, 刘式适, 赵强 2002 物理学报 51 718
Google Scholar
Liu S D, Fu Z T, Liu S K, Zhao Q 2002 Acta Phys. Sin. 51 718
Google Scholar
[13] Zhang S Y, Zhang W 1987 Acta Mech. Sin. 3 64
Google Scholar
[14] 刘志芳, 张善元 2007 固体力学学报 28 55
Google Scholar
Liu Z F, Zhang S Y 2007 Acta Mech. Solid. Sin. 28 55
Google Scholar
[15] 刘志芳, 张善元 2006 物理学报 55 628
Google Scholar
Liu Z F, Zhang S Y 2006 Acta Phys. Sin. 55 628
Google Scholar
[16] Liu Z F, Zhang S Y 2006 Appl. Math. Mech. 27 1431
Google Scholar
[17] 邓庆田, 罗松南, 彭亮 2009 应用力学学报 26 519
Google Scholar
Deng Q T, Luo S N, Peng L 2009 Chin. J. Appl. Mech. 26 519
Google Scholar
[18] Seadawy A R, Manafian J 2018 Results Phys. 8 1158
Google Scholar
[19] Baskonus H M, Bulut H, Atangana A 2016 Smart Mater. Struct. 25 035022
Google Scholar
[20] Bulut H, Sulaiman T A, Baskonus H M 2018 Opt. Quantum. Electron. 50 2
Google Scholar
[21] Wang Q 2001 Int. J. Solids. Struct. 38 8207
Google Scholar
[22] Xue C X, Pan E 2013 Int. J. Eng. Sci. 62 48
Google 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 091701
Google Scholar
[25] Ansari R, Oskouie F M, Gholami R, Sadeghi F 2016 Compos. Part B-Eng. 89 316
Google Scholar
[26] Ootao Y, Tanigawa Y 2000 Int. J. Eng. Sci. 38 47
Google Scholar
[27] Shiv P J 1992 Smart Mater. Struct. 1 80
Google Scholar
[28] Xue C X, Pan E, Zhang S Y 2011 Smart Mater. Struct. 20 105010
Google 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)
计量
- 文章访问数: 5575
- PDF下载量: 85