CN101554990B - 微机电系统的跨时间尺度仿真方法 - Google Patents

微机电系统的跨时间尺度仿真方法 Download PDF

Info

Publication number
CN101554990B
CN101554990B CN2008101901122A CN200810190112A CN101554990B CN 101554990 B CN101554990 B CN 101554990B CN 2008101901122 A CN2008101901122 A CN 2008101901122A CN 200810190112 A CN200810190112 A CN 200810190112A CN 101554990 B CN101554990 B CN 101554990B
Authority
CN
China
Prior art keywords
simulation
frequency group
time
mems
low
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
Application number
CN2008101901122A
Other languages
English (en)
Other versions
CN101554990A (zh
Inventor
苑伟政
秦子明
常洪龙
谢建兵
吕湘连
姜澄宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nantong Tongyang Mechanical and Electrical Manufacture Co., Ltd.
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN2008101901122A priority Critical patent/CN101554990B/zh
Publication of CN101554990A publication Critical patent/CN101554990A/zh
Application granted granted Critical
Publication of CN101554990B publication Critical patent/CN101554990B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种微机电系统的跨时间尺度仿真方法，属于微机电系统建模与仿真领域。该方法将微机电系统系统级模型中的组件分为高频组和低频组，在第一次仿真中将低频组作为仿真对象，用较大步长进行时域仿真，在第二次仿真中将低频组仿真结果与高频组组件作为仿真对象，用较小步长进行时域仿真，当需要用到低频组中组件的信号值时，对第一次仿真得到的仿真结果进行插值。使用该方法仿真时，低频组仿真时不考虑高频组的影响进行单独仿真，仿真规模减小；此外，对高频组进行仿真时只利用低频组仿真的结果，不耗费资源求解低频组中组件的模型，仿真规模减小，从而解决了现有技术中由时间尺度的差异带来的低效问题。

Description

微机电系统的跨时间尺度仿真方法
所属领域
本发明涉及一种微机电系统的跨时间尺度仿真方法，属于微机电系统(MEMS)建模与仿真领域。
现有技术
在微机电系统的系统级仿真中，微机电系统被分割成多个简单的组成部分，这些组成部分叫做组件；每个组件有若干端口，端口用来实现组件间信号的传输，相互连接的组件在对应连接端口处的信号值相等；用系统级求解器对由组件相互连接构成的网络进行数值求解即进行仿真，仿真对象可以是整个微机电系统，也可以是整个微机电系统的一部分；在系统级进行时域仿真时需要指定初始步长，步长指对系统模型进行数值求解时相邻两个求解点之间的时间间隔，初始步长指求解最初的几个点所用的步长，为了提高仿真效率，求解最初几个点之后的点时，步长通常由仿真器在仿真过程中进行动态调整。
在系统级进行微机电系统的仿真时，若组件在时间尺度上有很大差异，会给时域仿真带来困难。以微机械陀螺为例，其包含用于敏感角速度的机械部分和用于检测机械部分位移的电路部分，其机械部分频率在kHz量级，电路部分正弦载波频率在MHz量级，仿真时间由被测宏观物体的角速度决定，在1s量级。若将整个微机电系统作为仿真对象，为了准确地仿真其系统级行为，步长需设置成电路部分正弦载波频率的1/100，即0.01μs量级，产生1s的仿真时间需要仿真108步，非常耗费内存和硬盘空间等计算机资源。
现有的系统级仿真方法对这种跨时间尺度问题未提出有效解决方法，对时间尺度差异较大的组件没有分别对待，而是将它们都作为仿真对象的一部分，在一次仿真中按照同一步长进行仿真，这样就无法避免由时间尺度的差异带来的低效问题。Darrell Teegarden等人在文献How to model and simulate microgyroscope systems，如何进行微陀螺系统的建模与仿真，IEEE spectrum，July，1998，Qi Jing在文献Modeling and simulation for design of suspendedMEMS，悬浮MEMS的建模与仿真，Carnegie Mellon University，Ph.D.thesis，May 21，2003中对这种方法都有描述。
发明内容
本发明的目的在于公开一种微机电系统的跨时间尺度仿真方法，以解决现有技术中由时间尺度的差异带来的低效问题。该方法将微机电系统的系统级模型中的组件分为高频组和低频组，在第一次仿真中将低频组作为仿真对象，用较大步长进行时域仿真，在第二次仿真中将低频组仿真结果与高频组组件作为仿真对象，用较小步长进行时域仿真，当需要用到低频组中组件的信号值时，对第一次仿真得到的仿真结果进行插值。
本发明所采用的技术方案包括以下步骤：
步骤一：将微机电系统的系统级模型中的组件分到低频组和高频组，使得：低频组不包含高频组件，高频组对低频组的影响可以忽略；通常，微机电系统的机械组件、不受载波影响的电学组件和不受载波影响的机械-电学耦合组件等属于低频组件，光学组件、受载波影响的电学组件等属于高频组件，机械-光学耦合组件和受载波影响的机械-电学耦合组件等也属于高频组件。
步骤二：在系统级行为求解器中，将低频组作为仿真对象，用较大步长进行时域仿真，初始步长根据精度要求设置，取如下值中的最小值：最小时间常数的1/10、方波或脉冲信号源的最短上升或下降时间、正弦信号源周期的1/100，产生的仿真结果是与时间点对应的信号值。
步骤三：在系统级行为求解器中，将步骤二中时域仿真的结果与高频组作为仿真对象，用较小步长进行时域仿真，初始步长根据精度要求设置，取如下值中的最小值：最小时间常数的1/10、方波或脉冲信号源的最短上升或下降时间、正弦信号源周期的1/100；当需要用到低频组中组件的信号值时，对步骤二中时域仿真的结果进行插值，得到所需时间点处的信号值；插值可以采用分段线性插值、分段三次样条插值等方法。
本发明的有益效果是：1)由于高频组对低频组的影响很小，低频组仿真时可以不考虑高频组的影响进行单独仿真，这样，仿真规模减小，步长可以按照低频组中组件进行设置，防止了低频组和高频组在一次仿真中按同一步长仿真时，低频组步长被迫减小；2)对高频组进行仿真时只利用低频组仿真的结果，不耗费资源求解低频组中组件的模型，仿真规模减小。
下面结合附图和实施实例对本发明进一步说明。
附图说明
图1是本发明中仿真方法流程图
图2是实施例中加速度计机械部分的俯视图
图3是该加速度计机械部分的正视图
图4是分组前该加速度计的系统级模型
图5是对加速度计低频组进行时域仿真的示意图
图6是将低频组仿真结果和高频组作为仿真对象，进行时域仿真的示意图
图中，1-锚点，2-梁，3-质量块，4-电极，5-锚点芯核模型，6-梁芯核模型，7-质量块芯核模型，8-芯片运动芯核模型，9-端口联结线，10-仿真结果插值模型，11-电路模型，12-电容器芯核模型。
具体实施方式
下面结合一种微机械加速度计对本发明进行进一步说明。
参照图2，其结构包括锚点1、梁2、质量块3、电极4，图中实体部分表示硅，空心部分表示金属，当在x方向有加速度时，由于质量块3的质心不在上下两根梁2的连线上，质量块3会绕梁2发生转动，质量块3与两个电极4之间的电容就会发生变化，用检测电路检测这一变化，即可反映出加速度的大小。采用本发明中技术方案对其进行时域仿真的过程包括以下步骤：
步骤一：将该加速度计的系统级模型中的组件分到低频组和高频组。参阅图4，该加速度计的系统级模型由六种组件联结而成，其中，锚点芯核模型5、梁芯核模型6、质量块芯核模型7分别对应锚点1、梁2、质量块3，电路模型11对应检测电路，电容器芯核模型12对应质量块3与电极4之间的电容；芯片运动芯核模型8用于输入被测加速度，与其他组件之间采用虚拟联结。将该加速度计的系统级模型中的组件分到低频组和高频组：锚点芯核模型5、梁芯核模型6、质量块芯核模型7、芯片运动芯核模型8用于仿真机械结构对加速度的响应，属于低频组件中的机械组件，其中，外界输入正弦加速度100Hz，机械结构工作模态-绕z轴转动模态的谐振频率为679Hz，将这些组件分到低频组；电路模型11用于仿真检测电路，由于包含频率50kHz的载波，属于高频组件，电容器芯核模型12是机械-电学耦合组件，受电路模型11的影响，也属于高频组件，将这些组件分到高频组；这样，低频组不包含高频组件，由于检测电路频率离机械结构谐振频率很远，高频组对低频组的影响可以忽略。
步骤二：在系统级行为求解器Saber中，将低频组作为仿真对象，参阅图5，由于高频组对低频组的影响可以忽略，本应联结高频组件电容器芯核模型12的端口联结线9悬空；用较大步长进行时域仿真，由于只有正弦信号源，初始步长取低频组中频率最高组件的周期的1/100，即按照加速度计机械结构工作模态谐振频率的1/100设置为14.7μs，仿真时间0.05s；产生的仿真结果是与时间点对应的信号值，即与时间点对应的端口联结线9处的信号值，其中包括质量块芯核模型7的信号值。
步骤三：在系统级行为求解器Saber中，将步骤二中时域仿真的结果与高频组作为仿真对象，参阅图6，本应与低频组件质量块芯核模型7联结的电容器芯核模型12与仿真结果插值模型10联结；用较小步长进行时域仿真，由于只有正弦信号源，初始步长取高频组中频率最高组件的周期的1/100，即按电路模型11的周期的1/100取为0.2μs，仿真时间0.05s；当需要用到低频组中质量块芯核模型7的信号值时，对步骤二中时域仿真的结果进行插值，得到所需时间点处的信号值；插值采用分段线性插值。
结果表明，本发明所用方法的仿真效率大大高于现有技术中所用方法，本发明所用方法需要的时间和硬盘空间分别为460.8s和24.1M，现有技术需要的时间和硬盘空间分别为9466.7s和122.7M。前者需要的时间和硬盘空间约是现有技术中的1/21和1/5。

Claims (2)

1.一种微机电系统的跨时间尺度仿真方法，其特征在于，包括如下步骤：
步骤一：将微机电系统的系统级模型中的组件分到低频组和高频组，微机电系统的机械组件、不受载波影响的电学组件和不受载波影响的机械-电学耦合组件属于低频组件，光学组件、受载波影响的电学组件、机械-光学耦合组件和受载波影响的机械-电学耦合组件属于高频组件；
步骤二：在系统级行为求解器中，将低频组作为仿真对象进行时域仿真，初始步长根据精度要求设置，取低频组中如下值中的最小值：最小时间常数的1/10、方波或脉冲信号源的最短上升或下降时间、正弦信号源周期的1/100，产生的仿真结果是与时间点对应的信号值；
步骤三：在系统级行为求解器中，将步骤二中时域仿真的结果与高频组作为仿真对象进行时域仿真，初始步长根据精度要求设置，取高频组中如下值中的最小值：最小时间常数的1/10、方波或脉冲信号源的最短上升或下降时间、正弦信号源周期的1/100；当需要用到低频组中组件的信号值时，对步骤二中时域仿真的结果进行插值，得到所需时间点处的信号值。
2.一种如权利要求1所述的微机电系统的跨时间尺度仿真方法，其特征在于，所述的步骤三中所用的插值方法为分段线性插值方法或分段三次样条插值方法。
CN2008101901122A 2008-12-30 2008-12-30 微机电系统的跨时间尺度仿真方法 Expired - Fee Related CN101554990B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008101901122A CN101554990B (zh) 2008-12-30 2008-12-30 微机电系统的跨时间尺度仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008101901122A CN101554990B (zh) 2008-12-30 2008-12-30 微机电系统的跨时间尺度仿真方法

Publications (2)

Publication Number Publication Date
CN101554990A CN101554990A (zh) 2009-10-14
CN101554990B true CN101554990B (zh) 2011-05-25

Family

ID=41173304

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008101901122A Expired - Fee Related CN101554990B (zh) 2008-12-30 2008-12-30 微机电系统的跨时间尺度仿真方法

Country Status (1)

Country Link
CN (1) CN101554990B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104133966A (zh) * 2014-08-01 2014-11-05 电子科技大学 一种基于信号特征提取的故障电路行为建模方法

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112289368B (zh) * 2020-10-27 2023-07-25 英业达科技有限公司 通用型硬盘仿真配重治具
CN112949240B (zh) * 2021-03-22 2023-08-01 梁文毅 一种集中参数模型多物理场耦合仿真方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1673751A (zh) * 2004-03-24 2005-09-28 西北工业大学 微型惯性传感器件的芯核建模方法及芯核库
US7272801B1 (en) * 2003-03-13 2007-09-18 Coventor, Inc. System and method for process-flexible MEMS design and simulation
CN101136037A (zh) * 2007-08-15 2008-03-05 西北工业大学 微机电系统的宏模型建立方法
CN101315649A (zh) * 2008-07-17 2008-12-03 西北工业大学 含大量输入端口的微机电系统降阶建模方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7272801B1 (en) * 2003-03-13 2007-09-18 Coventor, Inc. System and method for process-flexible MEMS design and simulation
CN1673751A (zh) * 2004-03-24 2005-09-28 西北工业大学 微型惯性传感器件的芯核建模方法及芯核库
CN101136037A (zh) * 2007-08-15 2008-03-05 西北工业大学 微机电系统的宏模型建立方法
CN101315649A (zh) * 2008-07-17 2008-12-03 西北工业大学 含大量输入端口的微机电系统降阶建模方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
D.Teegarden et al.How To Model And Simulate Microgyroscope Systems.《IEEE Spectrum》.1998,第35卷(第7期),第66-75页. *
常洪龙等.支持任意流程的微机电系统设计工具.《机械工程学报》.2008,第44卷(第11期),第39-46页. *
霍鹏飞等.一种基于多端口组件网络的MEMS系统级建模方法.《机械科学与技术》.2005,第24卷(第2期),第192-194,198页. *
霍鹏飞等.基于组件网络方法的微加速度计建模与仿真.《航空学报》.2003,第24卷(第5期),第466-470页. *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104133966A (zh) * 2014-08-01 2014-11-05 电子科技大学 一种基于信号特征提取的故障电路行为建模方法
CN104133966B (zh) * 2014-08-01 2017-03-29 电子科技大学 一种基于信号特征提取的故障电路行为建模方法

Also Published As

Publication number Publication date
CN101554990A (zh) 2009-10-14

Similar Documents

Publication Publication Date Title
CN203275441U (zh) 传感器电路
CN102262699B (zh) 基于无网格伽辽金与质点弹簧耦合的软组织形变仿真方法
CN105260581B (zh) 舰船机电控制设备虚拟振动和冲击试验方法
CN103364590B (zh) 传感器电路和对微机电系统传感器进行测试的方法
CN101554990B (zh) 微机电系统的跨时间尺度仿真方法
CN102854998A (zh) 惯性感测装置
CN104091002B (zh) 一种车身接头模型的建立方法及装置
CN102840910B (zh) 具有振动测试功能的电子装置及建立振动测试算法的方法
Wei et al. The design of quadcopter frame based on finite element analysis
Gonuleri et al. CFD-Informed Rain Drop Impact Damage Predictions at Hypersonic Conditions
CN106844991B (zh) 气浮式振动控制系统空气弹簧刚度自平衡迭代识别方法
CN204679079U (zh) 一种mems三轴陀螺仪
CN105205250A (zh) 一种建模方法及建模系统
CN104020664A (zh) 一种基于偏微分方程的柔性机械臂干扰观测器设计方法
CN105103155B (zh) 分析装置以及分析方法
CN103592851A (zh) 多维波形模拟修正方法及控制装置
CN113507539B (zh) 一种获取终端设备跌落信息的方法和终端设备
CN102160057B (zh) 用于粒子法的界面粒子的判定方法及装置
CN105404702A (zh) 周视观察镜设备结构减振优化设计方法
CN106599395A (zh) 一种油浸变压器的噪声数值仿真计算方法
CN106055746A (zh) 自动识别最大模态组件的方法
CN111536995A (zh) 陀螺仪测试系统及方法
CN106595584A (zh) 天线姿态数据获取装置、获取方法及天线装置
CN112016234A (zh) 差容式地震计机械摆、测试系统及测试方法
CN108920733A (zh) 一种gil三支柱绝缘子及金属焊点松动故障仿真检测方法

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
ASS Succession or assignment of patent right

Owner name: NANTONG NANYANG ELECTROMECHANICAL MANUFACTURING CO

Free format text: FORMER OWNER: NORTHWESTERN POLYTECHNICAL UNIVERSITY

Effective date: 20140826

Owner name: NORTHWESTERN POLYTECHNICAL UNIVERSITY

Effective date: 20140826

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 710072 XI'AN, SHAANXI PROVINCE TO: 226600 NANTONG, JIANGSU PROVINCE

TR01 Transfer of patent right

Effective date of registration: 20140826

Address after: Haian Haian County town of 226600 in Jiangsu province Nantong City Road No. 9

Patentee after: Nantong Nanyang electromechanical Manufacturing Co., Ltd.

Patentee after: Northwestern Polytechnical University

Address before: 710072 Xi'an friendship West Road, Shaanxi, No. 127

Patentee before: Northwestern Polytechnical University

C56 Change in the name or address of the patentee

Owner name: NANTONG TONGYANG MECHANICAL AND ELECTRICAL MANUFAC

Free format text: FORMER NAME: NANTONG NANYANG ELECTROMECHANICAL MANUFACTURING CO., LTD.

CP01 Change in the name or title of a patent holder

Address after: Haian Haian County town of 226600 in Jiangsu province Nantong City Road No. 9

Patentee after: Nantong Tongyang Mechanical and Electrical Manufacture Co., Ltd.

Patentee after: Northwestern Polytechnical University

Address before: Haian Haian County town of 226600 in Jiangsu province Nantong City Road No. 9

Patentee before: Nantong Nanyang electromechanical Manufacturing Co., Ltd.

Patentee before: Northwestern Polytechnical University

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110525

Termination date: 20181230

CF01 Termination of patent right due to non-payment of annual fee