矿用压电俘能器建模与俘能特性研究1)

张晓宇 张旭辉,†,

* (西安科技大学机械工程学院,西安 710054)

† (陕西省矿山机电装备智能监测重点实验室,西安 710054)

引言

由于煤矿环境恶劣、危险系数大,为了避免发生人员伤亡事故,将工作人员从井下撤离出来,让采煤机自主运行,此时需要在采煤机上安装大量无线监测节点对运行状态进行实时监测,并将监测信息反馈给工作人员,保证采煤机正常运行.由于化学电池不能长时间为无线监测节点供电,供电成为瓶颈问题.压电俘能器可以将振动能转化为电能,研究俘能特性非常重要.

张伟等[1]发现增大压电俘能器的正弦激励强度可以增大功率输出.张颖等[2]研究提高计算俘能器磁力精度的方法.钱有华等[3]研究压电俘能器的阻尼比、势阱深浅程度和负载电阻对俘能效果的影响.曹东兴等[4]研究俘能器的俘能特性,加入磁力可以降低共振频率,提高电压输出.吴娟娟等[5]研究压电俘能器在不同中心频率的高斯白噪声激励下的输出电压.刘琦等[6-7]设计一种能使驱动力放大的双固支梁压电俘能结构,发现在弱高斯白噪声激励下可以发生突跳运动.马天兵等[8]研究不同形状压电悬臂梁在正弦激励下的输出功率.张忠华等[9]分析俘能器在不同强度正弦激励下的输出电压峰值.满大伟等[10-12]研究发现俘能器在合适的起振位置可以产生大幅值运动,并研究磁距、磁铁质量、负载电阻对输出响应的影响.李颍等[13]发现通过改变激励幅值可以使三稳态压电能量采集器实现阱间运动,提高功率输出.张强等[14-16]采用压电俘能装置为刨煤机的刨刀受力检测系统、采煤机扭矩检测系统和刮板输送机张力检测系统提供电量,通过实验发现在采煤机、刨煤机和刮板输送机工作时能够产生电量.Xie 等[17]发现具有非对称势阱的压电能量采集器在高斯白噪声激励下的输出电压比传统的压电能量采集器的输出电压大.Chen 等[18]发现双耦合梁的压电能量收集器加入弹簧可以降低势垒,在较低的激励下能够实现大幅值运动,选择合适的弹簧刚度可以提高功率输出.赵丽娟等[19-21]采用EDEM建立煤层模型,采用EDEM 与RecurDyn 耦合提取采煤机的振动加速度.Xia 等[22]发现增加压电能量采集器的端部质量可以降低共振频率,有利于提高能量收集效果.张宇等[23]发现悬臂梁在不同频率的简谐激励下存在跳跃现象.张广义等[24]发现梯形压电梁比矩形压电梁的输出功率大.Lu 等[25]发现加入多层压电材料可以加大能量采集器的输出电压.高扬等[26]采用振动能量采集装置为矿用无线监测系统供电,发现在谐振频率处的发电量最大,发电量与激励幅度呈正相关.刘建政等[27]采用压电-电磁-摩擦能量采集器实现矿用无线监测系统自供电.解胜东[28]采用压电叠堆俘能器为煤机检测系统提供电量,发现压电片的厚度与发电量呈正相关.

综上所述,目前研究压电俘能器受到采煤机采掘激励下的俘能特性非常少.因此,本文采用实验与数据拟合方法建立恢复力模型,磁化电流法建立磁力模型,拉格朗日函数建立动力学模型,RecurDyn提取滚筒、摇臂的截割方向加速度,龙格库塔法求解动力学模型,分析在不同磁距的俘能特性,并开展实验研究.

1 压电俘能器结构

压电俘能器结构,如图1 所示.在悬臂梁上固定压电薄膜,悬臂梁固定在基座上,磁铁A,B 分别固定在悬臂梁末端与基座上.目前压电梁大多采用直线结构,文献[29-30]发现曲线结构压电梁与直线结构压电梁受到同等作用力,采用曲线结构可以提高电压与功率输出,因此,本文压电梁采用曲线结构.压电俘能器分别受到滚筒、摇臂的截割方向加速度影响,由于有煤落在滚筒上,导致俘能器无法安装在滚筒上,现分别安装在摇臂上端、摇臂质心位置,如图2 的1,2,3,4 位置.

图1 压电俘能器结构Fig.1 Structure of piezoelectric energy harvester

图2 压电俘能器安装图Fig.2 Installation figure of piezoelectric energy harvester

悬臂梁与压电薄膜的线形段长度、拱形段半径分别为2 cm 和1 cm,悬臂梁与压电薄膜的宽度为0.8 cm,悬臂梁与压电薄膜的厚度分别为0.02 cm 与0.0 1 1 c m,悬臂梁与压电薄膜的密度分别为8300 kg/m3和1780 kg/m3,悬臂梁与压电薄膜的弹性模量分别为1.28×1011Pa 和3.0×109Pa,磁铁A,B 的长、宽、厚分别为0.5 cm,1 cm 和1 cm,磁化强度、真空磁导率和密度分别为5.5×105A/m,4π×10-7H/m 和7500 kg/m3.

2 压电俘能器数学模型与采煤机截割煤层仿真模型

2.1 压电俘能器数学模型

2.1.1 恢复力模型

采用文献[18]实验方法测量恢复力,如图3 所示,将线性端固定在基座上,采用YLK-10 测力计多次推动曲线端,对不同位移的恢复力取平均值,恢复力随位移变化,如图4 所示,对恢复力进行拟合,得到的模型为

图3 恢复力测量实验Fig.3 Restoring force measuring experiment

图4 恢复力随位移变化Fig.4 The change of restoring force with displacement

2.1.2 磁力模型

采用磁化电流法建立磁力模型,悬臂梁转动φ,磁铁A 同样转动φ.磁铁A,B 的位置图,如图5所示.

图5 磁铁A,B 的位置图Fig.5 Location figure of magnets A,B

在磁铁B 的中心建立坐标系,磁铁A 上下表面的中心点分别为O1和O2.∠AOC=φ,BC=d,AC=e.sinφ=e/(L+lA/2),lA和L分别为磁铁A 的长度和悬臂梁的水平长度.

磁铁A 受到磁铁B 的磁力Fm为

式中,µ0为真空磁导率,MA,SA,e,hA为磁铁A 的磁化强度、面积、位移、高度,Hx2(O2),Hx1(O1)为磁铁B 在O2,O1点的磁场强度,d为磁铁A,B 的水平距离.

磁铁B 在O处的磁场强度Hx(O)为

式中,MB,hB 和wB为磁铁B 的磁化强度、高度和宽度,x和y为O处的坐标值.

2.1.3 动力学模型

动力学模型为

式中,L,T,U和W为总能量、动能、弹性势能与磁场势能之和、电能.

动能为

式中,ρS和ρP为悬臂梁和压电薄膜的密度,AS和AP为悬臂梁和压电薄膜的侧面积,(x,t)为悬臂梁的x点在t时刻相对于基座的速度,和z(t)为基座的速度和位移,(L,t),w(L,t),I和m为磁铁A 相对基座的速度、位移、转动惯量和质量,t为时间,x为距悬臂梁基点的水平距离.

磁铁A 绕通过O点的Z轴做旋转运动,磁铁A 转动惯量[31]为

弹性势能与磁场势能之和为

电能为

式中,bP,Lp和为压电薄膜的宽度、水平长度和介电常数.

压电俘能器主要以一阶振动模态为主,故仅考虑一阶振动模态.对w(x,t)进行离散,采用一端夹紧、一端自由的悬臂梁容许函数近似振型函数ϕ(x)

拉格朗日方程为

式中,R为电阻,C为等效阻尼[32],w为谐振频率,ζ为阻尼比,M为等效质量.

将式(10)分别代入式(6)、式(8)和式(9),将分别算出的结果代入式(5),将代入式(5)算出的结果分别代入式(11)和式(12),得出动力学模型为

式中,K为等效刚度,HS为基础激励系数,ϑ 为机电耦合系数,CP为电容,(t)为加速度.

让振动台以加速度为20 m/s2,频率为13 Hz 的激励振动,振动一段时间后,停止振动,采用激光测振仪测量压电俘能器速度,采用coco80 采集速度信号.速度原始信号和滤波信号,如图6 和图7,速度信号衰减线与包络线,如图8 所示.

图6 原始信号Fig.6 Original signal

图7 滤波信号Fig.7 Filtering signal

图8 衰减线与包络线Fig.8 Attenuation and envelope line

对速度衰减信号的振幅进行指数函数拟合[33],为下式,其中ζw=0.934 7.

式(13)的Kr(t)项为恢复力,表示出恢复力随位移变化为线性关系,而压电俘能器结构为曲线型,恢复力随位移变化为非线性关系,Kr(t)项采用式(1)代替

在Matlab 中,采用龙格库塔法对式(13)和式(14)求解,得出电压

其中,r1(t)为位移,r2(t)为速度,r3(t)为电压.

2.2 采煤机截割煤层仿真模型

2.2.1 煤层模型

在EDEM 建立硬度为f4 和f6 的煤粒,半径均为50 mm,f4 煤的剪切模量和泊松比为9.56×108Pa和0.15,f6 煤的剪切模量和泊松比为1.72×109Pa和0.22,f6 煤的硬度相对于f4 煤的硬度大,将煤层几何模型导入EDEM,煤粒在煤层几何模型中生成,并采用黏结键黏结煤粒.黏结键参数计算公式为

式中,Sn,St和 σy分别为法向刚度、切向刚度和极限法向应力,G,μ,r和f为煤粒的剪切模量、泊松比、半径和普氏系数.

根据式(25)和式(26),f4 煤层的Sn,St和 σy为1.998 91×1010N/m3,8.690 9×109N/m3和4.0×107Pa,极限切向应力1.6×107Pa;f6 煤层的Sn,St和 σy为3.815 27×1010N/m3,1.563 6×1010N/m3和6.0×107Pa,极限切向应力2.4×107Pa.煤粒模型、煤层几何模型和煤层模型,分别如图9～图11 所示.

图9 煤粒模型Fig.9 Coal particle model

图10 煤层几何模型Fig.10 Coal seam geometric model

图11 煤层模型Fig.11 Coal seam model

2.2.2 采煤机动力学模型

(1)仿真模型的建立.将采煤机和刮板输送机中部槽几何模型导入到RecurDyn,仿真模型,如图12.

图12 仿真模型Fig.12 Simulation model

(2)约束的建立.为提取滚筒载荷,在滚筒与摇臂、摇臂与机身、摇臂与油缸、活塞与机身的连接处建立旋转副,在油缸与活塞、平滑靴与刮板输送机中部槽的连接处建立移动副,在导向滑靴与机身、平滑靴与机身、大地与刮板输送机中部槽的连接处建立固定副;为提取截割方向加速度,在滚筒与摇臂、摇臂与机身、摇臂与油缸、活塞与机身、导向滑靴与机身、平滑靴与机身的连接处建立弹性与阻尼元件,在油缸与活塞的连接处建立弹簧元件,在平滑靴与刮板输送机中部槽连接处建立移动副,在大地与刮板输送机中部槽的连接处建立固定副.

(3)驱动的建立.为提取滚筒载荷和截割方向加速度,在平滑靴与刮板输送机中部槽的移动副处建立移动速度驱动,大小为4 m/min,方向沿X轴负方向;为提取滚筒载荷,在滚筒与摇臂的旋转副处建立转速驱动,大小为3π rad/s,方向为逆时针.

采用提取滚筒载荷的动力学模型与煤层模型做采煤机截割煤层的仿真,可以提取滚筒载荷,截割煤层图,如图13 所示.将滚筒载荷导进提取截割方向加速度的动力学模型,可以提取截割方向加速度.

3 采煤机截割方向加速度的提取

截割f4 和f6 煤层的滚筒载荷分别如图14 和图15.

图14 截割f4 煤层的滚筒载荷Fig.14 Load of roller when exploiting f4 coal seam

图15 截割f6 煤层的滚筒载荷Fig.15 Load of roller when exploiting f6 coal seam

由图14 和图15 可知,截割f4 煤层的前滚筒的X向、Y向和Z向载荷有效值分别为3.694×105N,3.309×105N 和3.521×105N,后滚筒的X向、Y向、Z向载荷有效值分别为1.826×105N,2.905×105N 和3.404×105N.截割f6 煤层的前滚筒的X向、Y向和Z向载荷有效值分别5.344×105N,4.832×105N 和6.501×105N,后滚筒的X向、Y向、Z向载荷有效值分别2.881×105N,3.960×105N 和5.324×105N.

截割f4 和f6 煤层的采煤机各部位截割方向加速度,分别如图16 和图17 所示.

图16 截割f4 煤层的各部位截割方向加速度Fig.16 the exploiting direction acceleration of each part when exploiting f4 coal seam

图17 截割f6 煤层的各部位截割方向加速度Fig.17 the exploiting direction acceleration of each part when exploiting f6 coal seam

由图16 和图17 可知,振动从滚筒传到摇臂会减小,前滚筒3 向载荷相对于后滚筒大,且前、后滚筒的3 向载荷随煤层硬度的增大而增大,使前滚筒、后滚筒、前摇臂、后摇臂的截割方向加速度在依次减小且均随煤层硬度的增大而增大.截割f4 煤层的前滚筒、后滚筒、前摇臂和后摇臂的截割方向加速度的有效值分别为21.257 m/s2,19.580 m/s2,9.004 m/s2和8.455 m/s2,截割f6 煤层的前滚筒、后滚筒、前摇臂、后摇臂的截割方向加速度的有效值分别为30.724 m/s2,28.538 m/s2,13.095 m/s2和12.370 m/s2.

4 压电俘能器的俘能特性分析

压电俘能器受到采煤机各部位截割方向加速度发生振动,压电薄膜产生变形并输出电压.将截割方向加速度导入到动力学模型中,采用龙格库塔法求解,得出电压随时间变化趋势.

4.1 在截割f4 煤层的俘能特性分析

当受到前滚筒的截割方向加速度,在不同磁距的电压如图18 所示.

图18 受到前滚筒截割方向加速度的电压Fig.18 The voltage with exploiting direction acceleration of front roller

由图18 可知,在无磁力时,电压有效值为2.662 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为5.107 V,4.049 V 和3.845 V,加入磁力可以加大电压,磁距为12 mm,俘能特性较好,电压较大.

当受到前摇臂的截割方向加速度,在不同磁距的电压如图19 所示.

图19 受到前摇臂截割方向加速度的电压Fig.19 The voltage with exploiting direction acceleration of front swing arm

由图19 可知,在无磁力时,电压有效值为0.513 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为0.998 V,0.802 V 和0.851 V,加入磁力可以加大电压,磁距为12 mm,俘能特性较好,电压较大.

当受到后滚筒的截割方向加速度,在不同磁距的电压,如图20 所示.

图20 受到后滚筒截割方向加速度的电压Fig.20 The voltage with exploiting direction acceleration of after roller

由图20 可知,在无磁力时,电压有效值为2.655 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为3.078 V,4.224 V 和3.038 V,加入磁力可以加大电压,磁距为16 mm,俘能特性较好,电压较大.

当受到后摇臂的截割方向加速度,在不同磁距的电压,如图21 所示.

图21 受到后摇臂截割方向加速度的电压Fig.21 The voltage with exploiting direction acceleration of after swing arm

由图21 可知,在无磁力时,电压有效值为0.477 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为0.882 V,0.736 V 和0.756 V,加入磁力可以加大电压,磁距为12 mm,俘能特性较好,电压较大.

截割方向加速度呈随机波动,使电压呈随机波动.截割f4 煤层时,前滚筒、后滚筒、前摇臂和后摇臂的截割方向加速度在依次减小,使受到前滚筒、后滚筒、前摇臂和后摇臂的截割方向加速度,在俘能特性较好的磁距时,电压在依次减小.

4.2 在截割f6 煤层的俘能特性分析

当受到前滚筒的截割方向加速度,在不同磁距的电压,如图22 所示.

图22 受到前滚筒截割方向加速度的电压Fig.22 The voltage with exploiting direction acceleration of front roller

由图22 可知,在无磁力时,电压有效值为4.910 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为5.825 V,7.298 V 和5.571 V,加入磁力可以加大电压,磁距为16 mm,俘能特性较好,电压较大.

当受到前摇臂的截割方向加速度,在不同磁距的电压,如图23 所示.

图23 受到前摇臂截割方向加速度的电压Fig.23 The voltage with exploiting direction acceleration of front swing arm

由图23 可知,在无磁力时,电压有效值为1.060 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为1.307 V,1.592 V 和1.048 V,加入磁力可以加大电压,磁距为16 mm,俘能特性较好,电压较大.

当受到后滚筒的截割方向加速度,在不同磁距的电压,如图24 所示.

图24 受到后滚筒截割方向加速度的电压Fig.24 The voltage with exploiting direction acceleration of after roller

由图24 可知,在无磁力时,电压有效值为2.872 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为3.571 V,6.747 V 和3.798 V,加入磁力可以加大电压,磁距为16 mm,俘能特性较好,电压较大.

当受到后摇臂的截割方向加速度,在不同磁距的电压,如图25 所示.

图25 受到后摇臂截割方向加速度的电压Fig.25 The voltage with exploiting direction acceleration of after swing arm

由图25 可知,在无磁力时,电压有效值为0.918 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为1.125 V,1.397 V 和0.928 V,加入磁力可以加大电压,磁距为16 mm,俘能特性较好,电压较大.

截割f6 煤层时,前滚筒、后滚筒、前摇臂和后摇臂的截割方向加速度在依次减小且均随煤层硬度的增大而增大,使受到前滚筒、后滚筒、前摇臂和后摇臂的截割方向加速度,在俘能特性较好的磁距时,电压依次减小且均随煤层硬度的增大而增大.

5 压电俘能器俘能特性实验验证

采用RecurDyn 提取截割f4,f6 煤层的前滚筒截割方向加速度的功率谱密度,分别如图26 和图27所示.

图26 截割f4 煤层的功率谱密度Fig.26 Power spectral density when cutting f4 coal seam

图27 截割f6 煤层的功率谱密度Fig.27 Power spectral density when cutting f6 coal seam

对俘能特性进行实验验证.实验过程为: 将功率谱密度数据导进计算机中,计算机控制振动控制器产生激励信号,激励信号经过振动控制器、功率放大器传递到振动台,振动台产生振动.示波器与压电俘能器相连,采集电压信号.实验测试平台、压电俘能器,分别如图28 和图29 所示.

图28 实验测试平台Fig.28 Experimental testing platform

图29 压电俘能器Fig.29 Piezoelectric energy harvester

当受到截割f4 煤层的前滚筒截割方向加速度,在不同磁距的电压,如图30 所示.

图30 受到截割f4 煤层的前滚筒截割方向加速度的电压Fig.30 The voltage with exploiting direction acceleration of front roller when exploiting f4 coal seam

由图30 可知,在无磁力时,电压有效值为2.219 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为3.340 V,3.092 V 和2.456 V,相应的理论结果为: 在无磁力时,电压有效值为2.662 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为5.107 V,4.049 V 和3.845 V,将理论结果与实验结果对比发现,磁距均为12 mm 的俘能特性较好,电压较大,加入磁力可以加大电压,理论结果与实验结果定性吻合,验证了理论分析结果具有正确性.

当受到截割f6 煤层的前滚筒截割方向加速度,在不同磁距的电压,如图31 所示.

图31 受到截割f6 煤层的前滚筒截割方向加速度的电压Fig.31 The voltage with exploiting direction acceleration of front roller when exploiting f6 coal seam

由图31 可知,在无磁力时,电压有效值为2.304 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为4.379 V,4.959 V 和2.856 V,相应的理论结果为: 在无磁力时,电压有效值为4.910 V,磁距分别为12 mm,16 mm 和20 mm,电压有效值分别为5.825 V,7.298 V 和5.571 V,将理论结果与实验结果对比发现,磁距均为16 mm 的俘能特性较好,电压较大,加入磁力可以加大电压,理论结果与实验结果定性吻合,并通过实验结果发现,当受到截割f4 和f6 煤层的前滚筒截割方向加速度,在俘能特性较好的磁距时,电压随煤层硬度的增大而增大,通过理论结果分析发现同样有此现象,验证了理论分析结果具有正确性.

6 结论

(1)设计一种矿用压电俘能器,采用实验与数据拟合的方法建立恢复力模型,采用磁化电流法建立磁力模型,采用拉格朗日函数建立动力学模型.采用RecurDyn 提取滚筒、摇臂的截割方向加速度.

(2)压电俘能器受到前滚筒、后滚筒、前摇臂和后摇臂的截割方向加速度,在俘能特性较好的磁距时,电压依次减小且均随煤层硬度的增大而增大.截割f4 煤层时,磁距分别为12 mm,16 mm,12 mm 和12 mm 的俘能特性较好,电压有效值分别为5.107 V,4.224 V,0.998 V 和0.882 V,截割f6 煤层时,磁距均为16 mm 的俘能特性较好,电压有效值分别为7.298 V,6.747 V,1.592 V 和1.397 V.

(3)通过实验研究发现,压电俘能器受到截割f4 和f6 煤层的前滚筒截割方向加速度,在俘能特性较好的磁距时,电压随煤层硬度的增大而增大,磁距分别为12 mm 和16 mm 的俘能特性较好,电压有效值分别为3.340 V 和4.959 V,加入磁力可以加大电压,验证了理论分析结果具有正确性.