由于其微小的质量、极高的共振频率和品质因子, 微纳米机械系统(NMES/MEMS)广泛应用于许多科学技术领域[1-3].在20世纪90年代中期开始出现了基于MEMS的传感器, 它的出现使得量取用其他测量方法很难测得的微观物理量(如极小的质量、力、变形、温度等)成为可能, 因而成为探索微观领域的有力工具[4].基于微悬臂梁结构的质量传感器就是其广泛应用的一个典型代表.有研究表明, 微悬臂梁质量传感器的质量分辨率已经达到了10-18~10-21 g的量级[5].这样的质量分辨率使得其能够用于生物金属原子[6]和生物化学分子[7]的检测, 甚至是同位素的鉴定[8].在其众多应用中, 即使是单层的吸附层或者谐振器表面的不同材料喷涂层都会产生很大的表面应力, 影响悬臂梁的振动频率[9-10], 进而影响质量的精确检测; 然而在悬臂梁质量传感器的研究中却很少考虑表面应力的影响.
本文考虑悬臂梁表面应力对其振动频率的影响.将悬臂梁的表面应力等效分解, 通过引入应力系数对悬臂梁的振型函数加以修正, 再应用瑞利法建立悬臂梁共振频率和吸附微粒的位置和质量之间的关系.根据这个关系便可以利用悬臂梁的共振频率同时求得吸附微粒的位置和质量.
1 传感器模型的建立 1.1 悬臂梁的振型函数悬臂梁的表面应力可以按集中载荷模型分解为作用于悬臂梁自由端的等效轴向力和等效弯矩[11-12], 如图 1所示:σs为因为吸附作用而产生的悬臂梁表面应力; Te, Me分别为表面应力的等效轴向力和等效弯矩.等效弯矩对悬臂梁的振动频率没有影响[13], 这样受表面应力悬臂梁的振动问题就转化成为受单一轴向力的悬臂梁振动的一般问题.
现有一微悬臂梁, 在其任意位置lΔm处吸附有一质量为Δm微纳颗粒.如考虑吸附层表面应力引起的轴向力, 忽略其引起的其他影响(如对弯曲刚度和弯曲截面转动惯量的影响), 并且假设吸附层的厚度非常薄而可以将其忽略.根据欧拉-伯努利梁理论, 可得承受轴向力的悬臂梁的控制方程, 见式(1).
(1) |
式中:t为时间; E, I(对于方形截面梁I=bh3/12), T(T> 0时为轴向拉力, T < 0为轴向压力, 限于篇幅文中只讨论T> 0的情况)和m(m=ρS, ρ为密度, S=bh为弯曲截面面积)分别为杨氏模量, 弯曲截面转动惯量, 轴向力和梁的单位长度质量.为方便计算引入无量纲系数ξ=x/l(0≤ξ≤1), 设悬臂梁的第n(n=1, 2, 3, …)阶标准化挠度函数的形式为
(2) |
式中:ωn为悬臂梁的第n阶共振频率; Un(ξ)为悬臂梁的第n阶标准化的振型函数.将式(2)代入式(1)可得
(3) |
其中
(4) |
式中
(5) |
则可得悬臂梁横向振动的动力学特征方程:
(6) |
式中:
Ωn即为满足式(6)的正根, 根据式(6)可以得到任意应力系数情况下悬臂梁的共振频率系数Ωn, 进而由式(7)得到相应的共振频率:
(7) |
根据式(6)可以得到应力系数λ2和各阶的频率系数Ωn的关系, 如图 2所示.
由图 2可知, 随着应力系数λ2的增大, 各阶的频率系数Ωn也随之单调地增大.图中的虚线表示应力系数为零, 表征悬臂梁不受轴向力的一般情况; 虚线左边区域表示悬臂梁处于压曲状态; 右边区域表示悬臂梁处于拉伸状态.值得指出的是, 第一阶曲线与横轴的交点表示此时处于压曲状态的悬臂梁的一阶频率为零(图中箭头所指处), 对于更高阶的模态也有类似的结论.表 1列出了应力系数λ2分别为0, 10, 20和40时的悬臂梁的前3阶频率系数.
将由式(6)求得的频率系数Ωn代入式(4)中, 可得到悬臂梁第n阶标准化振型函数:
(8) |
式中:an为满足条件
图 3为悬臂梁在应力系数λ2为0, 10, 20和40的情况下的前三阶振型.
图 4所示为一基于微悬臂梁的质量传感器, Δm为吸附微粒质量.如果Δm远小于悬臂梁的质量M(M=lρS), 在微粒吸附前后可以假设悬臂梁的振型函数保持不变[14].根据瑞利-里兹法, 振动系统的最大动能Ekin等于其最大应变能Estrain.对于一个吸附微小质量的悬臂梁系统, 认为其最大势能不变, 因此可得Estrain=Ekin+Ekin, Δm, Ekin, Δm为质点质量的动能.其中
(9) |
将式(9)中各项代入瑞利方程并化简, 可得
(10) |
其中
因为α«1, 因此对式(10)进行泰勒展开得到
(11) |
式中
改写式(11)可以得到
(12) |
其中Δωn=ωn-ωn, Δm, ωn和ωn, Δm为吸附微粒前后悬臂梁的共振角频率.
在传感器的实际应用中, 式(12)中质量比α和吸附位置ξΔm均为待确定的未知量且与悬臂梁的振动模态和振型函数Un无关; 频率比Δωn/ωn为试验中的测量值, 其值和振动模态相关.对于第i阶振动模态, 频率比Δωi/ωi为已知定值, 由式(12)可知, 有无数的点(ξΔm, α)可以满足关系式Δωi/ωi=αUi2ξΔm/2.因此要确定两个未知量ξΔm和α, 则至多需要三阶振动模态的频率移动Δωn/ωn来唯一确定微粒的位置和质量.
2 仿真与验证用有限元仿真结果作为虚拟试验数据对本文方法进行验证.借助有限元软件ANSYS14.5建立一个三维的悬臂梁有限元模型.悬臂梁模型的材料参数和尺寸参数都为任意选定, 杨氏模量、泊松比和密度分别为100 GPa, 0.28和2 850 kg/m3; 梁的长、宽和厚度分别设定为100, 10和1 μm.吸附微粒为六面体, 其尺寸设为长、宽为2 μm高为1 μm.三维模型采用10节点的soild 86四面体单元, 该单元支持塑性、超弹性、蠕变、应力刚化、大挠度和大应变能力等分析.在模态分析之前进行静态力分析, 以便在有限模型中施加初始应力.不失一般性, 将悬臂梁所受轴向力的应力系数取为λ2=10, 对应于受拉伸轴向力的悬臂梁.在有限元仿真中, 将质量比Δm/M设为0和1×10-3(质量比为0时表示悬臂梁无吸附质量), 标准化位置ξΔm设为0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 共进行10次仿真(一次为吸附质量为0的情况, 其余9次为Δm/M=1×10-3, ξΔm=0.1~0.9的情况).每次仿真提取前三阶共振频率ω1, ω2和ω3(或ω1, Δm, ω2, Δm和ω3, Δm), 最后可以得到9组频率比数据Δω1/ω1, Δω2/ω2和Δω3/ω3.将这些数据代入到式(12)中便可以求得吸附微粒的位置ξΔm和质量比α.验证结果如图 5所示, 各阶曲线的交点(ξΔm, α)(箭头所指处)即为根据频率移动得出的计算结果, 图 5a~图 5i依次对应于(0.1, 10-3), (0.2, 10-3), …, (0.9, 10-3).这里需要指出的是, 如图 5a~图 5e中所示, n=1和n=2的曲线仅有唯一的交点, 所以针对这种情况, 只需Δω1/ω1和Δω2/ω2便可唯一确定ξΔm和α.而图 5f~图 5i, n=1和n=2的曲线有两个交点, 因此为唯一确定ξΔm和α需要引入n=3来确定其真实值.
本文提出一个同时确定单一吸附微粒位置和质量的方法, 该方法为单个细胞、药物分子等微小粒子的位置和质量的测定提供了更加简单便捷的方法.不同于以往的很多其他研究, 本文将表面应力对悬臂梁质量传感器的影响考虑在内, 对悬臂梁的振型函数进行了修正, 因此可以提高传感器的准确性和可靠性.利用该方法, 通过测量悬臂梁的最多三阶共振频率就可以确定出微粒的位置和质量.
[1] |
Stachiv I, Fedorchenko A I, Chen Y L.
Mass detection by means of the vibrating nanomechanical resonators[J]. Applied Physics Letters, 2012, 100(9): 093110.
DOI:10.1063/1.3691195 |
[2] |
Gil-Santos E, Ramos D, Martinez J, et al.
Nanomechanical mass sensing and stiffness spectrometry based on two-dimensional vibrations of resonant nanowires[J]. Nature Nanotechnology, 2010, 5(9): 641–645.
DOI:10.1038/nnano.2010.151 |
[3] |
Park J, Qin H, Scalf M, et al.
A mechanical nanomembrane detector for time-of-flight mass spectrometry[J]. Nano Letters, 2011, 11(9): 3681–3684.
DOI:10.1021/nl201645u |
[4] |
Boisen A, Dohn S, Keller S S, et al.
Cantilever-like micromechanical sensors[J]. Reports on Progress in Physics, 2011, 74(3): 036101.
DOI:10.1088/0034-4885/74/3/036101 |
[5] |
Dohn S, Svendsen W, Boisen A, et al.
Mass and position determination of attached particles on cantilever based mass sensors[J]. The Review of Scientific Instruments, 2007, 78(10): 103303.
DOI:10.1063/1.2804074 |
[6] |
Jensen K, Kim K, Zettl A.
An atomic-resolution nanomechanical mass sensor[J]. Nature Nanotechnology, 2008, 3(9): 533–537.
DOI:10.1038/nnano.2008.200 |
[7] |
Hanay M S, Kelber S, Naik A K, et al.
Single-protein nanomechanical mass spectrometry in real time[J]. Nature Nanotechnology, 2012, 7(9): 602–608.
DOI:10.1038/nnano.2012.119 |
[8] |
Chiu H Y, Hung P, Postma H W C, et al.
Atomic-scale mass sensing using carbon nanotube resonators[J]. Nano Letters, 2008, 8(12): 4342–4346.
DOI:10.1021/nl802181c |
[9] |
Berger R, Delamarche E, Lang H P, et al.
Surface stress in the self-assembly of alkanethiols on gold[J]. Science, 1997, 276(5321): 2021–2024.
DOI:10.1126/science.276.5321.2021 |
[10] |
Koev S T, Powers M A, Yi H, et al.
Mechano-transduction of DNA hybridization and dopamine oxidation through electrodeposited chitosan network[J]. Lab on a Chip, 2007, 7(1): 103–111.
DOI:10.1039/B609149K |
[11] |
Ren Q, Zhao Y P.
Influence of surface stress on frequency of microcantilever-based biosensors[J]. Microsystem Technologies, 2004, 10(4): 307–314.
DOI:10.1007/s00542-003-0329-4 |
[12] |
Cherian S, Thundat T.
Determination of adsorption-induced variation in the spring constant of a microcantilever[J]. Applied Physics Letters, 2002, 80(12): 2219–2221.
DOI:10.1063/1.1463720 |
[13] |
McFarland A W, Poggi M A, Doyle M J, et al.
Influence of surface stress on the resonance behavior of microcantilevers[J]. Applied Physics Letters, 2005, 87(5): 053505.
DOI:10.1063/1.2006212 |
[14] |
Nieradka K, Stegmann H, Gotszalk T.
Focused ion beam milling and deposition techniques in validation of mass change value and position determination method for micro and nanomechanical sensors[J]. Journal of Applied Physics, 2012, 112(11): 114509.
DOI:10.1063/1.4768715 |