CN104809279B - 矩形薄膜压膜阻尼等效电路 - Google Patents
矩形薄膜压膜阻尼等效电路 Download PDFInfo
- Publication number
- CN104809279B CN104809279B CN201510182481.7A CN201510182481A CN104809279B CN 104809279 B CN104809279 B CN 104809279B CN 201510182481 A CN201510182481 A CN 201510182481A CN 104809279 B CN104809279 B CN 104809279B
- Authority
- CN
- China
- Prior art keywords
- inductance
- secondary inductance
- film
- primary
- circuit units
- 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.)
- Expired - Fee Related
Links
Landscapes
- Coils Or Transformers For Communication (AREA)
Abstract
本发明公开了一种矩形薄膜压膜阻尼等效电路,包括初级电感回路和次级电感回路单元矩阵。将矩形薄膜均匀划分为(n+1)×(m+1)个薄膜单元,每个薄膜单元对应矩阵中的一个次级电感回路单元,各次级电感回路单元的结构和元件参数完全相同。初级电感回路的电感和各次级电感回路单元的电感形成耦合,各次级电感回路单元的电感间没有耦合。由离散近似的雷诺方程计算得到各薄膜单元的压膜阻尼压强值,将计算得到的阻尼压强值作为初级电感和各次级电感的耦合系数。不同薄膜运动速度下的各次级电感回路单元的次级电感上的电压值分布表征了薄膜下压膜阻尼的阻尼压强分布,通过本发明,可将流体域的问题转换为电域问题,且精度高、速度快。
Description
技术领域
本发明属于微机电系统(MEMS)系统级建模领域,尤其是一种矩形薄膜压膜阻尼等效电路。
背景技术
在MEMS器件中,有很多利用微机械结构振动性能工作,尤其是沿垂直方向的振动,矩形薄膜结构是采用该方式工作的MEMS器件的典型单元。薄膜振动时,其下部空间吸入和压出空气而产生的压膜阻尼效应是影响器件性能的重要物理效应之一。随着微结构特征尺寸的缩小,其表面积与体积比逐渐增大,空气阻尼效应变得异常显著。
目前的微机电器件系统级建模中,常将阻尼近似为一个阻尼系数,以阻尼系数和运动速度的乘积表征阻尼力的大小,而实际的空气阻尼行为是比较复杂的,近似的方法存在较大的误差。
发明内容
发明目的:提供一种矩形薄膜压膜阻尼等效电路,以解决现有技术存在的上述问题。
技术方案:一种矩形薄膜压膜阻尼等效电路,包括初级电感回路和次级电感回路单元矩阵,
所述初级电感回路包括依次连接的信号源V、电阻R0和初级电感L0,其中信号源V一端接地,另一端和电阻R0的一端连接,电阻R0的另一端与初级电感L0的一端连接,初级电感L0的另一端接地;
所述次级电感回路单元矩阵由若干个行列排布的次级电感回路单元组成,任意两个次级电感回路单元的结构和参数相同,任一次级电感回路单元包括依次连接的次级电感L1和电阻R1,电阻R1的一端接地,另一端与次级电感L1连接,次级电感L1的另一端接地;
所述初级电感回路中的初级电感L0和各次级电感回路单元的次级电感L1形成耦合,且初级电感L0与次级电感L1的匝比为1:1,各次级电感回路单元的次级电感L1间彼此独立。
优选的,矩形薄膜的x方向均匀划分为(n+1)份,y方向均匀划分为(m+1)份,构成(n+1)×(m+1)个薄膜单元,其中n、m为不小于0的正整数;所述次级电感回路单元矩阵包括(n+1)×(m+1)个初级电感回路单元,且薄膜单元与初级电感回路单元一一对应。每个薄膜单元的尺寸为:hx=Lx/(n+1)、hy=Ly/(m+1),其中,Lx和Ly分别为矩形薄膜在x、y方向上的尺寸。各个次级电感回路单元中的电阻阻值相同,各个次级电感回路单元中的次级电感的电感值相同。
根据位置不同将次级电感L1标记为次级电感Li,j,其中,i=0,1,…n,j=0,1,…m;所述初级电感回路中的初级电感L0和各次级电感回路单元的次级电感Li,j的耦合系数由下式获得:
式中,x、y为坐标方向,μeff为空气粘滞系数,d为薄膜下空气层的厚度,v Z为薄膜运动速度,Ω={(x,y)|0≤x≤Lx,0≤y≤Ly},pi,j为初级电感和位置(i,j)处的次级电感Li,j的耦合系数,其中,i=0,1,…n;j=0,1,…m。
有益效果:本发明由离散近似的雷诺方程计算得到各薄膜单元的压膜阻尼压强值,将计算得到的阻尼压强值作为初级电感和各次级电感的耦合系数。不同薄膜运动速度下的各次级电感回路单元的次级电感上的电压值分布表征了薄膜下压膜阻尼的阻尼压强分布,通过本发明的等效电路模型将流体域的问题转换为电域问题,具有精度高和速度快的优点。
附图说明
图1a和图1b分别是本发明的两种实施例。
图2是本发明的等效电路的初级电感回路结构。
图3是本发明的等效电路的次级电感回路单元结构。
具体实施方式
微机电器件系统级模型采用等效电路或硬件描述语言进行描述,为了能将实际阻尼以系统级模型的常用方法表示,需要对阻尼进行系统级建模,符合系统级模型的描述要求。阻尼是一种阻碍薄膜运动的力,并且是一种分布力,即在膜的中间和四周,其阻尼力的大小是不一样的。因此,经过上述创造性的分析思考,申请人认为可以将具有分布特性的阻尼力类比于呈一定分布的电压。
为便于理解,采用图1a和图1b两种方式描述本发明的等效电路。如图1a和图1b所示,矩形薄膜压膜阻尼等效电路由两部分组成:初级电感回路101和次级电感回路单元矩阵102(图1中虚线框住的部分)。
所述初级电感回路101由信号源V 101-1、电阻R0101-2和初级电感L0101-3依次连接而成,信号源V 101-1一端接地,另一端和电阻R0101-2的一端连接,电阻R0101-2的另一端与初级电感L0101-3的一端连接,初级电感L0101-3的另一端接地。
所述次级电感回路单元矩阵102由(n+1)×(m+1)个结构和元件参数完全相同的次级电感回路单元组成。各次级电感回路单元由次级电感Li,j201,和电阻Ri,j202连接而成,其中,i=0,1,…n;j=0,1,…m,下同。电阻Ri,j202的一端接地,电阻Ri,j202的另一端与次级电感Li,j201的一端连接,次级电感Li,j201的另一端接地。各个次级电感回路单元中的电阻(Ri,j,i=0,1,…n;j=0,1,…m)阻值相同,各个次级电感回路单元中的次级电感(Li,j,i=0,1,…n;j=0,1,…m)的电感值相同。
所述初级电感回路(101)中的初级电感L0101-3和各次级电感回路单元的次级电感Li,j,i=0,1,…n;j=0,1,…m形成耦合,且初级电感L0101-3与次级电感Li,j的匝比为1:1,i=0,1,…n;j=0,1,…m,。各次级电感回路单元的次级电感(Li,j,i=0,1,…n;j=0,1,…m)各单元彼此独立即没有耦合。
描述压膜阻尼运动的基本方程是雷诺方程,对于以一定速度运动的平板,简化后的雷诺方程为椭圆方程:
其中,μeff为空气粘滞系数,p为气体压膜压强,d为薄膜下空气层的厚度,x、y为坐标方向,vZ为薄膜运动速度。
假设薄膜x方向的尺寸是Lx,薄膜y方向的尺寸是Ly,将矩形薄膜沿x方向均匀划分为(n+1)份,沿y方向均匀划分为(m+1)份,得到(n+1)×(m+1)个薄膜单元,薄膜单元x和y方向的尺寸为hx=Lx/(n+1),hy=Ly/(m+1)。
将雷诺方程式(1)用二阶差分进行离散近似,得到差分方程:
式中,pi,j为位置(i,j)处的压膜压强,Ω={(x,y)|0≤x≤Lx,0≤y≤Ly}。
由式(2)可以计算得到各薄膜单元的阻尼压强。采用类比的方法将计算得到的阻尼压强值作为初级电感101-3和各次级电感(Li,j,i=0,1,…n;j=0,1,…m)的耦合系数。对于给定的初级电感回路信号源信号电压,即可获得各次级电感的输出电压。因初级电感和次级电感的匝比为1:1,同时因为耦合系数直接得自于阻尼压强,因此,各次级电感输出电压的分布与阻尼压强分布相同。通过本发明的等效电路模型,将流体域的问题转换为电域问题。
总之,目前对压膜阻尼的分析方法主要采用有限元方法,对于简单的情况也可以采用解析求解的方法。不论是那种方法,都无法将分析得到的阻尼压强分布直接用于包含了微机电器件和信号处理系统的系统级仿真。目前在系统级仿真中采用的方法是利用一个阻尼系数来近似,存在着较大的误差,而且,对于不同空气隙其阻尼系数也存在较大的差异。
而本发明提供的矩形薄膜压膜阻尼等效电路,显然,其与信号处理系统具有描述上的一致性。同时,其计算的基本依据又是雷诺方程,因此,本发明的模型既具有有限元分析的精度又具有电路分析的速度。
以上详细描述了本发明的优选实施方式,但是,本发明并不限于上述实施方式中的具体细节,在本发明的技术构思范围内,可以对本发明的技术方案进行多种等同变换,这些等同变换均属于本发明的保护范围。另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合。为了避免不必要的重复,本发明对各种可能的组合方式不再另行说明。此外,本发明的各种不同的实施方式之间也可以进行任意组合,只要其不违背本发明的思想,其同样应当视为本发明所公开的内容。
Claims (1)
1.一种矩形薄膜压膜阻尼等效电路,其特征在于,包括初级电感回路(101)和次级电感回路单元矩阵(102),
所述初级电感回路(101)包括依次连接的信号源V、电阻R0和初级电感L0,其中信号源V一端接地,另一端和电阻R0的一端连接,电阻R0的另一端与初级电感L0的一端连接,初级电感L0的另一端接地;
所述次级电感回路单元矩阵(102)由若干个行列排布的次级电感回路单元组成,任意两个次级电感回路单元的结构和参数相同,任一次级电感回路单元包括依次连接的次级电感L1和电阻R1,电阻R1的一端接地,另一端与次级电感L1连接,次级电感L1的另一端接地;
所述初级电感回路中的初级电感L0和各次级电感回路单元的次级电感L1形成耦合,且初级电感L0与次级电感L1的匝比为1:1,各次级电感回路单元的次级电感L1间彼此独立;
矩形薄膜的x方向均匀划分为(n+1)份,y方向均匀划分为(m+1)份,构成(n+1)×(m+1)个薄膜单元,其中n、m为不小于0的正整数;
所述次级电感回路单元矩阵包括(n+1)×(m+1)个次级电感回路单元,且薄膜单元与初级电感回路单元一一对应;
每个薄膜单元的尺寸为:hx=Lx/(n+1)、hy=Ly/(m+1),其中,Lx和Ly分别为矩形薄膜在x、y方向上的尺寸;
各个次级电感回路单元中的电阻阻值相同,各个次级电感回路单元中的次级电感的电感值相同;
根据位置不同将次级电感L1标记为次级电感Li,j,其中,i=0,1,…n,j=0,1,…m;所述初级电感回路中的初级电感L0和各次级电感回路单元的次级电感Li,j的耦合系数由下式获得:
式中,x、y为坐标方向,μeff为空气粘滞系数,d为薄膜下空气层的厚度,vZ为薄膜运动速度,Ω={(x,y)|0≤x≤Lx,0≤y≤Ly},pi,j为初级电感和位置(i,j)处的次级电感Li,j的耦合系数,其中,i=0,1,…n;j=0,1,…m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510182481.7A CN104809279B (zh) | 2015-04-16 | 2015-04-16 | 矩形薄膜压膜阻尼等效电路 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510182481.7A CN104809279B (zh) | 2015-04-16 | 2015-04-16 | 矩形薄膜压膜阻尼等效电路 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104809279A CN104809279A (zh) | 2015-07-29 |
CN104809279B true CN104809279B (zh) | 2017-11-03 |
Family
ID=53694096
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510182481.7A Expired - Fee Related CN104809279B (zh) | 2015-04-16 | 2015-04-16 | 矩形薄膜压膜阻尼等效电路 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104809279B (zh) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3262048B2 (ja) * | 1997-10-01 | 2002-03-04 | 株式会社村田製作所 | 圧電共振子およびそれを用いた電子部品 |
EP1517443B1 (en) * | 2003-09-12 | 2011-06-29 | Panasonic Corporation | Thin film bulk acoustic resonator, method for producing the same, filter, composite electronic component device, and communication device |
CN104015406B (zh) * | 2014-06-11 | 2015-12-30 | 青岛理工大学 | 嵌入式共固化网格阻尼复合材料结构及其制作工艺 |
-
2015
- 2015-04-16 CN CN201510182481.7A patent/CN104809279B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104809279A (zh) | 2015-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106094525A (zh) | 一种基于分数阶微积分的终端滑模控制器及控制方法 | |
CN101873210B (zh) | 网状形多涡卷混沌电路及产生多涡卷的方法 | |
CN102623062A (zh) | 一种忆阻器仿真模型 | |
CN108920811A (zh) | 一种用于直升机飞行仿真的模拟方法及系统 | |
CN104809279B (zh) | 矩形薄膜压膜阻尼等效电路 | |
CN103729518A (zh) | 一种简单的忆阻器仿真器 | |
CN105892297A (zh) | 一种自适应分数阶动态滑模控制算法 | |
Yaghmaie et al. | Computational modeling of finite deformation piezoelectric material behavior coupling transient electrical and mechanical fields | |
Wu et al. | Novel microfluidic device integrated with a fluidic‐capacitor to mimic heart beating for generation of functional liver organoids | |
CN206042010U (zh) | 网格多翼混沌电路 | |
CN101051328A (zh) | 空间连续变形mems的可重用参数化组件建模方法 | |
Lin et al. | Two-cavity MEMS variable capacitor for power harvesting | |
CN108732941A (zh) | 一种多维磁悬浮隔振平台的控制系统及仿真方法 | |
CN110457791A (zh) | 基于特征线的龙格-库塔时间离散的流体力学有限元算法 | |
CN100570616C (zh) | 微机电变截面固支梁系统级宏模型的建立方法 | |
CN101554990B (zh) | 微机电系统的跨时间尺度仿真方法 | |
Toshiyoshi et al. | A mixed-design technique for integrated MEMS using a circuit simulator with HDL | |
CN108875182B (zh) | 模拟压电陶瓷作动器迟滞非线性的饱和电容模型建模方法 | |
CN109214058B (zh) | 一种机电惯容器结构选型与参数设计方法 | |
CN106980740B (zh) | 一种电力电容器芯子振动模型建模方法 | |
CN103838927B (zh) | 电阻模块的spice电路仿真模型、spice仿真方法和装置 | |
Acciani et al. | Multi-physics simulation of a wind piezoelectric energy harvester validated by experimental results | |
CN113094928B (zh) | 一种人行单步荷载模型构建方法 | |
CN104392038B (zh) | 不同介质柔性充液结构的有限元建模方法 | |
CN103336440B (zh) | 一种含噪变频振荡衰减信号极值群的往复搜索定迹法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
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: 20171103 |