CN101046371A - 基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法 - Google Patents
基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法 Download PDFInfo
- Publication number
- CN101046371A CN101046371A CN 200710057124 CN200710057124A CN101046371A CN 101046371 A CN101046371 A CN 101046371A CN 200710057124 CN200710057124 CN 200710057124 CN 200710057124 A CN200710057124 A CN 200710057124A CN 101046371 A CN101046371 A CN 101046371A
- Authority
- CN
- China
- Prior art keywords
- amount
- bow
- curvature
- finite difference
- micro
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
本发明公开一种基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法,是由如下过程完成:第一步,应用相移显微干涉法测量被测微/纳悬臂梁的表面弯曲数据;第二步,应用有限差分法解析出被测微/纳悬臂梁的弯曲量;第三步,将第一步得到的表面弯曲数据与第二步得到的弯曲量进行比较,从而求出测量弯曲量与解析弯曲量之间差值最小时的曲率值。本发明采用相移显微干涉法测量悬臂梁表面的弯曲量,具有高精度和高分辨率,对待测器件是非破坏的无损测量,应用有限差分法解析出悬臂梁弯曲量的实现简单,建模更加接近于真实情况,结果准确可靠,可用于微悬臂梁或纳悬臂梁(具有微米宽度和纳米厚度的纳悬臂梁)曲率的测量。
Description
技术领域
本发明涉及一种微/纳机电系统的微/纳悬臂梁曲率测量方法。特别是涉及一种实现简单,结果准确可靠,具有高精度和高分辨率的基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法。
背景技术
微/纳机电系统(Micro/Nano ElectroMechanical System,MEMS/NEMS)是在微电子技术的基础上发展起来的关键尺寸在微米/纳米量级的集传感、信息处理和执行于一体的集成系统。MEMS/NEMS的研究涉及多门学科,其制造涉及到设计、材料、工艺、测量、控制、能源以及系统集成和封装等技术,而其应用更是囊括了汽车、通讯、医疗、生化、航空航天、科学仪器等领域。在MEMS/NEMS器件从设计加工到投入应用的过程中,对其结构的力学特性进行测量是至关重要的一步。
由于MEMS/NEMS器件往往是由硅和金属等材料通过化学气相沉积、溅射和电镀等方法形成薄膜后,再经过光刻、刻蚀和牺牲层腐蚀等方法制作而成的,因此这种通过表面微机械加工制备的薄膜往往存在着不可忽视的残余内应力和应力梯度,它极大地影响着器件的成品率和使用性能,如何实现对微/纳结构残余应力和应力梯度的精确、在线测量是MEMS/NEMS测量的重要任务之一,而在对薄膜的应力特性进行测量时,精确确定薄膜结构的曲率值往往是首要的一步。因此,如何简便而精确地实现曲率测量对MEMS/NEMS器件应力特性的分析来说具有非常重要的意义。
微/纳悬臂梁是MEMS/NEMS器件的基本构件之一,它是谐振器、RF开关以及生化传感器等许多MEMS/NEMS器件中的关键结构。微/纳悬臂梁的结构简单,加工制造方法容易,并且建模方法相对比较成熟,是微/纳测试中最常用的测试结构之一,通过微/纳悬臂梁测试结构可以很好地评价薄膜结构的应力特性。显微干涉法利用光学干涉来实现微小器件表面的轮廓测量,特别是相移显微干涉的应用,极大地提高了测量的垂直分辨率和精度,已成为微/纳测量中最常用的测量手段之一。而有限差分法是一种重要的数值算法,在工程分析中已经得到了广泛的应用。
发明内容
本发明所要解决的技术问题是,提供一种实现简单,结果准确可靠,具有高精度和高分辨率的基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法。
本发明所采用的技术方案是:一种基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法,测量过程包括:应用相移显微干涉法测量被测微/纳悬臂梁的表面弯曲数据;应用有限差分法解析出被测微/纳悬臂梁的弯曲量;将前两个过程得到的表面弯曲数据和弯曲量进行比较,从而求出测量弯曲量与解析弯曲量之间差值最小时的曲率值。
所述的应用相移显微干涉法测量被测微/纳悬臂梁的表面弯曲数据,其过程包括:应用相移显微干涉测量系统采集被测微/纳悬臂梁的干涉图,通过相位提取和相位展开后得到悬臂梁表面的原始高度信息;在得到的原始高度信息上以支撑部分的表面作为零基准面去除倾斜误差并得到悬臂梁表面的垂直弯曲数据。
所述的应用有限差分法解析出被测微/纳悬臂梁的弯曲量,是在悬臂梁的二维模型上应用有限差分法解析梁的弯曲量与曲率之间的关系,并考虑了初始端转角这个非理想边界条件的影响。
所述的求出测量弯曲量与解析弯曲量之间差值最小时的曲率值,是采用拟牛顿算法和最小二乘法两种算法中的一种,通过比较悬臂梁的测量弯曲量与解析弯曲量,从而得到使得两者差值最小时的最佳曲率匹配值。
本发明的基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法,采用相移显微干涉法测量悬臂梁表面的弯曲量,具有高精度和高分辨率,对待测器件来讲是非破坏的无损测量;模型解析时考虑了不为零的初始端转角这个非理想边界条件的影响,使得建模更加接近于真实情况;使用有限差分法进行近似计算,当结点取得越多时精度越高,这在以像素为单位的微/纳结构表面测量结果分析中可以达到很高的精度;提供两种数值分析方法实现测量和解析弯曲量的比较,可以将两种方法得到的结果进行对比,验证结果的准确性。此方法实现简单,结果准确可靠,可用于微悬臂梁或纳悬臂梁(具有微米宽度和纳米厚度的纳悬臂梁)曲率的精确测量。
附图说明
图1是本发明的相移显微干涉示意图;
图2是本发明的微/纳悬臂梁模型;
图3是本发明的微/纳悬臂梁模型结点示意图;
其中:
1:被测表面 2:参考镜
3:分光镜 4:干涉图
5:光源 6:悬臂梁
7:初始端转角 8:梁的支撑部分
9:结点 -1,0,1……N:结点的序号
L:梁长 λ:相邻结点之间的间隔
具体实施方式
下面结合附图和实施例对本发明的基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法做出详细说明。
本发明的基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法,其测量过程包括:
1、应用相移显微干涉法测量被测微/纳悬臂梁的表面弯曲数据:
应用相移显微干涉测量系统采集被测微/纳悬臂梁的干涉图,通过相位提取和相位展开后得到悬臂梁表面的原始高度信息;在得到的原始高度信息上以梁支撑部分的表面为基准面进行坐标旋转去除倾斜误差,得到悬臂梁表面相对于支撑平面的垂直高度信息。
具体作法是:首先将由原理显示于图1中的相移显微干涉系统得到的被测微/纳悬臂梁的表面干涉图像4采集下来并保存;然后,对干涉图像进行相位提取和相位展开得到被测表面的初始轮廓图像;最后,以图2中梁的支撑部分8的表面为基准面进行坐标旋转消除倾斜误差,并将支撑部分8的表面高度作为零高度基准,从而得到悬臂梁6表面的垂直弯曲量,取梁的中轴线上的一列数据备用。
2、应用有限差分法解析出被测微/纳悬臂梁的弯曲量:
该步骤通过在悬臂梁的二维模型上应用有限差分法解析梁的弯曲量与曲率之间的关系,并考虑了初始端转角这个非理想边界条件的影响。
具体作法是:按照图3中的结点9的分布将悬臂梁分成N个单元,用有限差分法计算出各结点9的垂直高度与曲率及图2中初始端转角7之间的关系。
微/纳悬臂梁的变形符合小变形原理,满足Bernoulli-Euler方程:
式中M为弯矩,E为杨氏模量,I为惯性矩,k为曲率,θ为转角。由于悬臂梁在z向的变形相对于梁的长度实际上是非常小的,所以我们将梁上各点的曲率k视为相对于x的常数。将整个悬臂梁用N+1个结点分割成N个小单元,并引入虚结点-1,代入边界条件:梁初始端的位移为零,而初始端的转角往往不为零,将方程(1)和(2)应用一阶和二阶有限差分法公式(3)和(4)进行近似计算,求得梁的弯曲量(挠度z(x))与曲率k及梁初始端转角θ0(结点0处的转角)之间的解析关系如式(5)所示。式(3)和(4)中的λ为相邻结点之间的间隔。
3、计算最佳曲率匹配值:
将前面得到的测量弯曲量与解析弯曲量,通过用拟牛顿算法或者最小二乘法来进行比较,求出测量弯曲量与解析弯曲量之间差值最小时的曲率值,作为最佳曲率匹配值,该步骤有两种具体作法:
(1)采用拟牛顿算法:
拟牛顿算法是一种数值最优化方法,可用于求解函数的最小值问题。将梁的中轴线的弯曲量与公式模拟的弯曲量逐点对比,取使得均方根差值最小的k和θ0值为最终解。均方根差值的表达式如式(6)所示,其中f为弯曲量的测量值与解析模拟值之间的均方根差值,z(xi)是悬臂梁xi处的弯曲测量值。
只要从某个初值出发,拟牛顿算法就可以得到使得目标函数f的值最小时k和θ0的最优解。
(2)采用最小二乘法:
最小二乘法是如果将式(6)构造成一个超定方程组,如式(7)所示,等式左边为与k和θ0有关的解析弯曲量,等式右边为弯曲测量值,由此可以使用最小二乘法求出使得两个弯曲量近似相等时的k和θ0解。
本发明公开和揭示的方法可借鉴本文公开内容,尽管本发明的方法已通过较佳实施例进行了描述,但是本领域技术人员明显能在不脱离本发明内容、精神和范围内对本文所述的方法进行拼接或改动,或增减某些环节的,更具体地说,所有相类似的替换和改动对本领域技术人员来说是显而易见的,它们都被视为包括在本发明精神、范围和内容中。
Claims (4)
1.一种基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法,其特征在于,测量过程包括:应用相移显微干涉法测量被测微/纳悬臂梁的表面弯曲数据;应用有限差分法解析出被测微/纳悬臂梁的弯曲量;将前两个过程得到的表面弯曲数据和弯曲量进行比较,从而求出测量弯曲量与解析弯曲量之间差值最小时的曲率值。
2.根据权利要求1所述的基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法,其特征在于,所述的应用相移显微干涉法测量被测微/纳悬臂梁的表面弯曲数据,其过程包括:应用相移显微干涉测量系统采集被测微/纳悬臂梁的干涉图,通过相位提取和相位展开后得到悬臂梁表面的原始高度信息;在得到的原始高度信息上以支撑部分的表面作为零基准面去除倾斜误差并得到悬臂梁表面的垂直弯曲数据。
3.根据权利要求1所述的基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法,其特征在于,所述的应用有限差分法解析出被测微/纳悬臂梁的弯曲量,是在悬臂梁的二维模型上应用有限差分法解析梁的弯曲量与曲率之间的关系,并考虑了初始端转角这个非理想边界条件的影响。
4.根据权利要求1所述的基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法,其特征在于,所述的求出测量弯曲量与解析弯曲量之间差值最小时的曲率值,是采用拟牛顿算法和最小二乘法两种算法中的一种,通过比较悬臂梁的测量弯曲量与解析弯曲量,从而得到使得两者差值最小时的最佳曲率匹配值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100571243A CN100547347C (zh) | 2007-04-11 | 2007-04-11 | 基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100571243A CN100547347C (zh) | 2007-04-11 | 2007-04-11 | 基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101046371A true CN101046371A (zh) | 2007-10-03 |
CN100547347C CN100547347C (zh) | 2009-10-07 |
Family
ID=38771167
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007100571243A Expired - Fee Related CN100547347C (zh) | 2007-04-11 | 2007-04-11 | 基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100547347C (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107894200A (zh) * | 2017-09-28 | 2018-04-10 | 东南大学 | 一种多双端固支梁结构的曲率传感器 |
CN110132473A (zh) * | 2019-05-21 | 2019-08-16 | 中国石油大学(华东) | 弯矩测量机构及测量方法、抽油杆弯矩测量装置及测量方法 |
CN111442753A (zh) * | 2020-03-24 | 2020-07-24 | 石家庄铁源工程咨询有限公司 | 一种高速铁路连续梁线性的测量方法、处理装置及存储介质 |
CN112556906A (zh) * | 2020-10-29 | 2021-03-26 | 瑞声新能源发展(常州)有限公司科教城分公司 | 一种测量薄膜不同方向应力梯度的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3130177B2 (ja) * | 1992-07-07 | 2001-01-31 | 松下電子工業株式会社 | 結晶評価装置および結晶評価方法 |
CN1170164C (zh) * | 2002-11-14 | 2004-10-06 | 上海交通大学 | 测量纳米级晶须材料杨氏模量的方法 |
-
2007
- 2007-04-11 CN CNB2007100571243A patent/CN100547347C/zh not_active Expired - Fee Related
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107894200A (zh) * | 2017-09-28 | 2018-04-10 | 东南大学 | 一种多双端固支梁结构的曲率传感器 |
CN110132473A (zh) * | 2019-05-21 | 2019-08-16 | 中国石油大学(华东) | 弯矩测量机构及测量方法、抽油杆弯矩测量装置及测量方法 |
CN111442753A (zh) * | 2020-03-24 | 2020-07-24 | 石家庄铁源工程咨询有限公司 | 一种高速铁路连续梁线性的测量方法、处理装置及存储介质 |
CN111442753B (zh) * | 2020-03-24 | 2021-08-31 | 石家庄铁源工程咨询有限公司 | 一种高速铁路连续梁线性的测量方法、处理装置及存储介质 |
CN112556906A (zh) * | 2020-10-29 | 2021-03-26 | 瑞声新能源发展(常州)有限公司科教城分公司 | 一种测量薄膜不同方向应力梯度的方法 |
CN112556906B (zh) * | 2020-10-29 | 2021-12-24 | 瑞声新能源发展(常州)有限公司科教城分公司 | 一种测量薄膜不同方向应力梯度的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100547347C (zh) | 2009-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101520385A (zh) | 微机电系统薄膜材料力学性能与可靠性测试方法和装置 | |
WO2013188131A1 (en) | Microelectromechanical system and methods of use | |
CN111103073A (zh) | 一种多参量协同敏感的谐振式压力传感器及其制备方法 | |
CN100402978C (zh) | 微纳结构3d轮廓测量中基于坐标变换的倾斜误差补偿方法 | |
CN101046371A (zh) | 基于显微干涉和有限差分的微/纳悬臂梁的曲率测量方法 | |
Muntwyler et al. | Three-axis micro-force sensor with sub-micro-Newton measurement uncertainty and tunable force range | |
WO2010017409A1 (en) | Force, pressure, or stiffness measurement or calibration using graphene or other sheet membrane | |
Foland et al. | One-dimensional nanograting-based guided-mode resonance pressure sensor | |
CN1289893C (zh) | 弯曲曲率和弯曲方向的同时测量方法 | |
McShane et al. | Young's modulus measurement of thin-film materials using micro-cantilevers | |
Mariani et al. | On-Chip testing: A miniaturized lab to assess sub-micron uncertainties in polysilicon MEMS | |
CN101046374A (zh) | 基于原子力显微镜力曲线的纳米梁厚度的非破坏测量方法 | |
CN104568236A (zh) | 一种薄膜应力测试装置及方法 | |
Zhao et al. | Temperature-insensitive silicon resonant pressure sensor by thermal stress control | |
Clark et al. | Practical techniques for measuring MEMS properties | |
US20150012230A1 (en) | Method of measuring micro- and nano-scale properties | |
Tsuchiya | Evaluation of mechanical properties of MEMS materials and their standardization | |
CN100442007C (zh) | 量测薄膜应变的方法及其装置 | |
Chen et al. | Micromachined SiO 2 microcantilever for high sensitive moisture sensor | |
Zhang et al. | Mechanical characterization of post-buckled micro-bridge beams by micro-tensile testing | |
Schalko et al. | Improved load-deflection method for the extraction of elastomechanical properties of circularly shaped thin-film diaphragms | |
Davies et al. | Engineering in-and out-of-plane stress in PECVD silicon nitride for CMOS-compatible surface micromachining | |
Shi-Hao | A novel method for determining Young's modulus of thin films by micro-strain gauges | |
Zhang et al. | Microgripping force measuring device based on SU-8 microcantilever sensor | |
Chang-Soo | Development and verification of PZT actuating micro tensile tester for optically functional materials |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091007 Termination date: 20210411 |