1. 引言
谐振式压力传感器因其高精度、高灵敏度得到越来越多的关注 [1] 。典型谐振式传感器多以谐振梁作为敏感元件,通过测量谐振梁振动频率的变化来测量待测压力。实际应用中谐振梁的高动态应用日益受到重视,在动态载荷的作用下,谐振梁需做出快速精确的响应,为保证测量结果的准确性,对谐振梁的动态响应进行分析显得尤为重要。目前对谐振梁振动特性的研究局限于移动力或只对梁的结构进行分析 [2] [3] [4] ,本文试图从理论上探讨动态载荷对欧拉梁振动特性的作用机理,建立更一般形式的动态轴向载荷下梁的振动的数学模型。该模型在特定情况下可以对应为Mathieu方程,然而Mathieu方程在求解过程中有一定的局限性,因此采用等效电路的方法对该模型进行分析 [5] 。此外,考虑到实际振动中谐振梁的阻尼损耗,分别通过对谐振梁边界条件的理论计算以及Comsol软件仿真分析了阻尼对谐振梁振动特性的影响。
2. 谐振梁振动数学模型建立
假设谐振梁在动态轴向载荷的作用下为伯努利-欧拉梁的模型 [6] ,受力模型如图1。以微元dx段为研究对象,受力分析如图2。
图中各项的物理意义:
F——梁承受的载荷;
Figure 1. Schematic diagram of bending vibration of double end fixed beam under dynamic axial force
图1. 动态轴向载荷作用下双端固支梁弯曲振动示意图
Figure 2. Force analysis of the micro sec- tion dx of the resonant beam
图2. 谐振梁微段dx受力分析
——梁的密度;
E——梁的弹性模量;
I——梁横截面对中性轴的惯性矩。
A、L、h、b分别是梁的横截面积、长度、高度、宽度。设
表示谐振梁上距原点为x处的截面在t时刻的横向位移,由力矩平衡公式与弯矩平衡公式共同联立得梁横向弯曲振动的方程为:
(1)
用分离变量法求该方程的解,设方程(1)的解为:
(2)
其中
为梁振动的固有频率,
为梁振动的模态函数。则式(1)可变形为:
(3)
解得梁自由振动振型公式:
(4)
谐振梁在实际振动时,总有能量通过谐振梁两端的固支端损失,因此将谐振梁的边界等效为弹簧支撑情况代入有阻尼的边界条件,其中k表示弹性系数
(5)
解得
(6)
谐振梁的各项参数取值如表1所示:
Table 1. Parameters of the resonant beam
表1. 谐振梁各项参数取值
代入各项已知数值后,计算出此时谐振梁频率为88 kHz,经Matlab计算可得,在仿真计算中发现,中所定k值的改变不影响两方程交点的改变,因此,得到一阶结果为
(7)
对应参数
(8)
得出梁自由振动–阶振型公式
(9)
对式(3)两边都乘以
,再对
进行积分,利用分部积分法以及齐次边界条件可得:
(10)
式(10)简化可写为:
(11)
其中A,B为常量,该式为一般载荷形式作用下梁的振动模型。
输入载荷为正弦时,通过如式(10)的变量代换:
(12)
将式(10)进行转换后,可得到如式(13)所示的标准形式的Mathieu方程。
(13)
将公式中的参数全部统一为国际单位制后,经过数值计算,式(11)为:
(14)
通常可以利用Mathieu方程可以推导出在动态输入载荷下谐振梁的振动模型,但是Mathieu方程对输入载荷的形式有要求,因此有局限性。故本文采取建立等效电路的方法模拟动态轴向载荷下谐振梁振动方程特性,来实现动态输入载荷的一般性 [7] 。
3. 谐振梁振动模型电路模拟
根据已经建立的谐振梁振动模型建立输出谐振梁振动信号的等效电路,如图3所示。
其中F(t)即对应系统中的输入载荷,通过改变输入来模拟动态力的加载模式。
根据数值要求及电路自激振荡的条件确定电路各项参数后,对电路输入不同的激励进行仿真,如图所示,分别为对电路输入载荷为0 mv的直流电压、−20 mv的直流电压、频率为10 kHz载荷为10 mv及20 mv的交流电压和频率为50 kHz载荷为10 mv的交流电压(如图4)。
根据电路仿真结果可以看出,当电压频率及幅值改变时,输出响应的频率和幅值也相应改变,且当输入信号增大时,输出响应出现频率被调制的现象。符合在机械系统中,
因外载荷增大而系统刚度增大,振荡频率升高而振幅下降的现象。
4. 阻尼对谐振梁振动特性影响
在谐振梁的分析中,阻尼的分类主要有以下几种:1) 锚点损耗,当谐振器振动时,谐振梁的动能通过边界传递到锚点,再通过锚点将动能传递到谐振器基底,最后造成能量耗散;2) 表面能量损耗,由表面晶格缺陷以及表面含有吸附物造成的表面应力损耗;3) 空气阻尼,当谐振器工作在空气中时,空气阻尼损耗是谐振器的主要能量损耗方式,为了提高谐振器Q值,必须对谐振器采用真空封装以减小阻尼损耗。然而绝对真空的封装是很难实现的,空气阻尼损耗也不可能完全消除;4) 热弹性阻尼,当谐振器振
图3. 谐振梁动态响应等效电路图
动时,谐振梁内部温度会发生变化,使得谐振梁内部的温度场为非均匀温度场,由于材料热胀冷缩的性质,会产生相应的应力、应变,阻碍谐振器的振动,造成能量损耗,可以将其称为热弹性阻尼 [7] [8] 。
目前很多传感器均采用真空封装,这样将会减小谐振器的空气阻尼损失,虽然绝对真空的封装是很难实现的,空气阻尼损耗也不可能完全消除,但是其影响在真空封装的作用下,已有很大程度的减少。热弹性阻尼及表面损耗阻尼均属于内部损耗,在此不多加讨论。因此,本部分主要讨论锚点损耗对谐振梁的影响。
① 等效电路方法分析
考虑在谐振梁的振动方程中加入阻尼项,则原振动方程的一般形式变为
(15)
其中c为阻尼系数,则原始的电路模型变为(如图5)。
加入由比例放大电路构成的阻尼模块,从第一级积分输出引出到加法器的输入端,设第一级积分输
图5. 有阻尼谐振梁动态响应等效电路图
出为
,阻尼模块放大倍数为
,根据输入输出关系列出如下方程:
(16)
整理后可得到,
(17)
带入相关参数后对电路进行仿真,设定阻尼系数,此时可得到如下仿真结果,其中红色实线表示加入阻尼后的电路输出情况,黑色实线表示未加入阻尼的电路输出情况(图6):
仿真结果显示,在电路中加入阻尼后几乎不会对电路的谐振频率产生影响,但幅值因阻尼的加入有所改变。
② COMSOL方法分析
实验中利用Comsol Multiphysics 5.0对谐振梁进行仿真分析,仿真模型中,设置了一个参数k,定义为阻尼系数,即在两侧的连接边界设置固定载荷,载荷T = kv,该式为矢量式,即受到来自外部的应力(Pa)正比于边界处的速度(m/s),那么如此定义的阻尼系数单位为[Pa*s/m]。仿真时,使用参数化扫描对k的值进行变化,k分别取值为0,50,100,150,200,250,300。其中,0则代表自由边界。
图7和图8所示坐标图为以上7个k值下,谐振梁中心点的振幅,已归一化到最大值为1。
仿真结果显示,在不同阻尼的影响下,所得谐振频率基本相同,随着两端阻尼系数的增加,谐振梁能量在锚点损耗,振幅逐渐减小,但谐振频率基本保持不变,由此可以得出结论,不同阻尼导致的边界条件下,阻尼的改变仅对谐振梁的幅值变化有影响,对谐振梁的谐振频率并不存在很大的影响。文献 [9]
Figure 7. Simulation of resonant beam’s amplitude by comsol
图7. Comsol对谐振梁仿真中心振幅变化
Figure 8. Simulation of resonant beam’s frequency by comsol
图8. Comsol对谐振梁仿真频率变化
利用小参数摄动法所得出的结论中,带阻尼类Mathieu方程与不带阻尼类Mathieu方程仅在幅值项上有区别,对频率项没有影响,与实验结果相符合 [9] 。
5. 结论
载荷的输入对谐振式压力传感器敏感元件谐振梁的振动频率会产生影响。在等效电路实验中,输出信号会根据输入电压信号的频率与幅值改变而产生相应的变化,电路中的调制现象与机械系统的振动特性相符合,即当增加外载荷时,机械系统振动频率变大,当外载荷达到最大值时,振动系统刚度最大,谐振频率最高,但会出现振幅下降的现象。此外,通过在等效电路中加入阻尼项及利用Comsol仿真后得到,阻尼仅影响谐振梁的振幅大小,对频率没有影响,因此对频率型谐振器来说,在工作过程中其准确性不会受到阻尼影响。
此外,本文中所设计的谐振梁模型及实验仿真都是在理想条件下进行,实际谐振器振动过程中将面临不同类型的输入叠加及梁体构造的不同,后期工作也将对此展开进一步的研究探讨。
致谢
感谢总装基金以及“长江学者和创新团队发展计划资助(IRT_16R02)”的资助。