专利名称：：微机电系统的角度参数化宏建模方法
技术领域：
：本发明涉及一种微机电系统(MEMS)的角度参数化宏建模方法，属于微机电系统设计与模型降阶领域。
背景技术：
：随着MEMS设计技术及工艺水平的迅速发展，MEMS器件结构也变得越来越复杂，大量新型器件相继产生。然而，结构对称性依然是发展中MEMS器件的一大特点，g卩在MEMS的设计中，大量存在几何形状完全相同，而初始方位不同的功能部件或微结构。对这些微结构的有限元模型进行自由度縮减，建立能够表征微结构输入输出行为的低阶模型的过程称为宏建模过程，所提取的低阶模型也称为宏模型。现有的宏建模方法由于无法实现角度参数化，必须对这些结构分别处理，各自独立地进行宏建模，而宏建模本身需要大量耗时的数值计算，从而导致这种方法计算代价高，效率低。因此，实现这些结构的角度参数化宏模型对于避免重复的宏建模过程，提高整个MEMS的设计效率具有重要意义。CarnegieMellon大学的QiJing在其博士论文"ModelingandSimulationforDesignofSuspendedMEMS"中，针对MEMS集总参数模型，采用坐标变换方法实现了角度参数化，但其集总参数模型采用矩阵结构分析和静电场等经典理论方法，经过解析推导而提取，与本发明将角度参数化应用于数值降阶的宏建模方法不同。此外，近年来，参数化宏建模方法取得较大进展，德国IMTEK的E.B.Rudnyi教授，美国MIT的LucaDaniel博士以及UCDavis的BaiZhaojun教授等人对参数化宏建模方法进行了广泛而深入的研究，然而其主要针对的是MEMS结构的几何尺寸及材料属性参数化，对于角度参数化宏建模方法方面的研究尚未涉及。
发明内容本发明的目的在于提出一种微机电系统的角度参数化宏建模方法，对于MEMS设计中几何形状完全相同，初始方位不同的微结构，只需完成其中一个微结构的宏建模，通过设置相应的角度参数，即可立即实现其它任意方位微结构的宏建模。本发明提出的微机电系统的角度参数化宏建模方法首先从相同几何形状的微结构中任意选取一个作为参考微结构，采用数值降阶方法提取其宏模型；然后根据微结构方位之间的关系，采用矩阵坐标变换方法建立其它方位微结构相对参考微结构的力与位移之间的转换方程；最后，将力与位移之间的转换方程代入宏模型的行为方程中，实现MEMS的角度参数化宏建模。该方法主要步骤如下步骤一进行参考微结构的宏建模，对于几何形状完全相同，初始方位不同的微结构，任意选择其中一个部件作为参考微结构，采用三维造型工具建立微结构的三维实体模型，并使用有限元方法进行模态分析，从分析结果中提取参考微结构的质量、刚度矩阵，建立微结构的二阶动力学行为方程wc=五(0式中f为时间，g(O为节点位移向量，尸为载荷矩阵，《为输入载荷向量，^为输出位移向量，￡为输出矩阵，M,D,Keiwxw分别为参考微结构的质量、阻尼和刚度矩阵，MEMS结构一般采用Rayleigh结构阻尼模型，即Z)-aM+ZX，其中"、々是和结构模态阻尼比相关的参数。上述二阶动力学行为方程经等价变换后，可转换为一阶状态方程卜c=￡rx(0(2)其中X(0:<formula>formulaseeoriginaldocumentpage5</formula>『是一个JVxW的非奇异矩阵，在结构力学中，通常取『-M，这样上述状态方程的系统矩阵C、G保持了原二阶状态方程的系数矩阵的对称特性。从上述变换过程可知，二阶动力学行为方程变换为一阶状态方程后，其方程规模也由W增加到2iV。以^二-G一C、i-G—力作为输入，应用Arnoldi算法对状态方程(2)进行降阶，得到参考微结构的宏模型斗r(3)其中,是时间变量，XfW是宏模型的状态变量，C,《为宏模型的系统矩阵，A是宏模型的输入矩阵，；是宏模型的输出矩阵，r为宏模型的规模，由于，<<^，所以宏模型的求解效率相对原始状态方程大幅提高。步骤二根据微结构方位之间的关系，采用矩阵坐标变换方法建立其它方位微结构相对参考微结构的力与位移转换方程。参阅图1，本发明首先定义两个坐标系以描述微结构在三维空间内的运动。把与参考微结构初始方位重合的坐标系称为总体坐标系，用C-Xc3^c表示，简称c坐标系。为了便于描述微结构的转动行为、把固连于参考微结构某点/7的坐标系称为局部坐标系，用，Xpj;^p表示，简称户坐标系。根据定义可知，c坐标系是静止坐标系，而户坐标系随部件运动。这样，通过p坐标系的转动可描述微结构的转动行为。本发明采用以下转动顺序引入确定p坐>示系在c坐标系的方位角，即p坐标系绕c坐标系x、_y、z坐标轴的相对转角，也即其它微结构相对参考微结构的转角，用(《,/，/)表示最初两坐标系完全重合，而后顺序通过三次转动达到P坐标系所在方位，这三次转动为绕Zc转一个^角，使^-X^pZp由最初与C-JCchZc重合的位置转到p-X^^/2c，第二次绕轴转^角，使，xw>^zc到达新的户-x^2^W位置，第三次绕xp2转a角，使p-Xp教;^到达j9-xp>ypzp的最终位置，把这种转动方式记为c^"^l""^24—p，其转动过程参阅图2。设r/为变换过程中c坐标系到;坐标系的转换矩阵，则根据矩阵坐标变换原理，可知J1/7=cos"0-sin々010sin"0coscos;ksin;k0—sinycosy00011000cosasina0-sinacos"」^..」__cos;kcos/sinycos〃一sin々-sinycosa+c0s7sin/sin"cosycos+sin;ksin>9sin"cos/sin"sin;rsin"+cos/sin/cos"—cos7sin"+sin;rsin"cos"cos々cos"设与户坐标系固连的微结构在c坐标系的位移与转角为^=[^,"y,"z,^，s，"::T，微结构所受力及转矩为Kp=[^,Fv,F:,7;，rv,7;]\则可得其它方位微结构相对参考微结构的力与位移转换方程c_p_c_(4)其中^,，W—-A,S,a:分别为沿x、_y、Z方向的平动位移和绕X、》Z轴的转角，K,Fy,F:,7;，7;,7;分别为沿x、_y、z方向的力和绕x、；;、z轴的转矩。"￡为参考微结构在c坐标系中的位轻与转角，《为参考微结构在c坐标系中的力与转矩，式(4)表明，微结构在c坐标系与p坐标系中的动力学关系可通过矩阵坐标变换来表征。步骤三将力与位移之间的转换方程代入宏模型的行为方程中，实现宏模型的角度参数化。将步骤二所得力与位移之间的转换方程(4)代入到步骤一所得的参考微结构的宏模型(3)中，同时根据其它微结构相对参考微结构的方位，设置相对转角参数值(《,/，/)，迅速提取其它方位微结构的宏模型，实现角度参数化的宏建模。6本发明的有益效果是采用矩阵坐标变换方法实现角度参数化的宏建模，通过设置相应的角度参数，即可避免结构完全相同，初始方位不同的MEMS微结构的重复宏建模过程，从而使基于宏模型的MEMS系统级建模与仿真速度加快，进而提高整个MEMS的设计效率。图1c坐标系与户坐标系的定义示意图图2c坐标系到坐标系的转换过程示意图图3本发明具体实施例中平面加速度计的结构示意图图中，l一左上角折叠梁，2—右上角折叠梁，3—左下角折叠梁，4一右下角折叠梁，5—质量块具体实施例下面结合一种平面加速度计的折叠梁角度参数化的宏建模过程对本发明进行进一步说明。参照附图3，加速度计的框架结构由左上角折叠梁l、右上角折叠梁2、左下角折叠梁3、右下角折叠梁4和质量块5构成，它关于坐标系o;^z的原点o中心对称以及x、y坐标轴的轴对称性，在加速度计的左上、右上、左下、右下方均有一个结构复杂的折叠梁。折叠梁的角度参数化宏建模过程具体步骤如下步骤一选择左上角折叠梁l作为参考结构，采用有限元软件ANSYS，选择beaml88单元进行模态分析，其有限元模型含有356个节点，其中包含两个输入输出节点，由于beaml88的每个节点包含6个自由度，即沿x、y、z坐标轴的平动和绕x、y、z坐标轴的转动，因此左上角折叠梁1共包含12个输入输出端口，模型的质量刚度矩阵的规模为#=356^6=2136阶。提取有限元模型的质量、刚度矩阵，并令《=0,^=6><1(^，可得左上角折叠梁l的二阶动力学行为方程为-M柳+6x10_6■+聰=尸i^(0(5)其中，？为时间，^)为节点位移向量，12是输入矩阵，M,《ef2136分别是质量、刚度阵，i^,[^,^，i^7;,，^，7^i^，F卢i^，7;"7;2，7;2:ref是输入载荷向量，表示作用于左上角折叠梁1的全部12个端口上的力与转矩，￡ei2136x12为输出矩阵，Wcpl^["d,"卢":l'"ri,"卢A""rf，"W":2,^2,^2,^-e是输出位移向量'表示左上角折叠梁l7的全部12个端口上的输出位移及转角。某些情况下还可能包括部件的速度向量，取决于所关注的物理量。将方程(6)转换为一阶状态方程如下fC雄)+GxO^(f)以J^-G^C、i二G"S作为输入，应用Arnoldi算法对状态方程(6)进行降阶，并取宏模型的阶数为4，得到左上角折叠梁l的宏模型jCA(/)+W)=5(7)其中Z是时间变量，x,(/)eW是宏模型的状态变量，C;,《ei^s是宏模型的系统矩阵，&ei4x12是输入矩阵，^ei4"是输出矩阵。步骤二定义ax^为c坐标系，参阅图3，它是一个静止坐标系。由于参考微结构所在p坐标系和c坐标系重合，也即左上角折叠梁l相对其本身的方位角为(0，0，0)。而右上角折叠梁2、左下角折叠梁3和右下角折叠梁4所在方位相对左上角折叠梁1的转角分别下表所示<table>tableseeoriginaldocumentpage8</column></row><table>因此它们所在的p坐标系相对O^Z坐标系的转换矩阵分别如下:<formula>formulaseeoriginaldocumentpage8</formula>其中772、T/3、774分别表示右上角折叠梁2、左下角折叠梁3和右下角折叠梁4所在的户坐标系相对a^z坐标系的转换矩阵，将它们代入方程(4)中，可得右上角折叠梁2、左下角折叠梁3和右下角折叠梁4相对左上角折叠梁1的力与位移之间的转换方程分别为<formula>formulaseeoriginaldocumentpage8</formula><formula>formulaseeoriginaldocumentpage9</formula>步骤三将步骤二所得力与位移的转换方程(8)、(9)、(10)依次代入方程(7)输入载荷向量和输出位移向量中，即可获取右上角折叠梁2、左下角折叠梁3和右下角折叠梁4的宏模型行为方程，进而实现角度参数化的宏建模。这样，本实例首先选取左上角折叠梁为参考部件并进行宏建模，然后通过设置相应的角度参数，以提取其它方位折叠梁的宏模型，从而实现了角度参数化的宏建模。权利要求1.一种微机电系统的角度参数化宏建模方法，其特征在于，包括如下步骤步骤一进行参考微结构的宏建模，对于几何形状完全相同，初始方位不同的微结构，任意选择其中一个部件作为参考微结构，建立其宏模型<mathsid="math0001"num="0001"><math><![CDATA[<mfencedopen='{'close=''><mtable><mtr><mtd><msub><mi>C</mi><mi>r</mi></msub><msub><mover><mi>x</mi><mo>·</mo></mover><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>G</mi><mi>r</mi></msub><msub><mi>x</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>B</mi><mi>r</mi></msub><msub><mi>F</mi><mi>c</mi></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>c</mi></msub><mo>=</mo><msubsup><mi>L</mi><mi>r</mi><mi>T</mi></msubsup><msub><mi>x</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math></maths>其中t是时间变量，xr(t)是宏模型的状态变量，Cr、Gr为宏模型的系统矩阵，Br是宏模型的输入矩阵，Fc为输入载荷向量，uc为输出位移向量，Lr是宏模型的输出矩阵；步骤二根据微结构方位之间的关系，采用矩阵坐标变换方法建立其它方位微结构相对参考微结构的力与位移转换方程，首先定义两个坐标系以描述微结构在三维空间内的运动，把与参考微结构初始方位重合的坐标系称为总体坐标系，用c-xcyczc表示，简称c坐标系，把固连于参考微结构某点p的坐标系称为局部坐标系，用p-xpypzp表示，简称p坐标系，采用以下转动顺序引入确定p坐标系在c坐标系的方位角，即p坐标系绕c坐标系x、y、z坐标轴的相对转角，也即其它微结构相对参考微结构的转角，用(α，β，γ)表示最初两坐标系完全重合，而后顺序通过三次转动达到p坐标系所在方位，这三次转动为绕zc转一个γ角，使p-xpypzp由最初与c-xcyczc重合的位置转到p-xp1yp1zc，第二次绕yp1轴转β角，使p-xp1yp1zc到达新的p-xp2yp1zp2位置，第三次绕xp2转α角，使p-xp2yp1zp2到达p-xpypzp的最终位置，把这种转动方式记为<mathsid="math0002"num="0002"><math><![CDATA[<mrow><mi>c</mi><munderover><mo>→</mo><msub><mi>z</mi><mi>c</mi></msub><mi>γ</mi></munderover><mn>1</mn><munderover><mo>→</mo><msub><mi>y</mi><mn>1</mn></msub><mi>β</mi></munderover><mn>2</mn><munderover><mo>→</mo><msub><mi>x</mi><mn>2</mn></msub><mi>α</mi></munderover><mi>p</mi><mo>;</mo></mrow>]]></math></maths>变换过程中c坐标系到p坐标系的转换矩阵Tcp为<mathsid="math0003"num="0003"><math><![CDATA[<mrow><msubsup><mi>T</mi><mi>c</mi><mi>p</mi></msubsup><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><mi>cos</mi><mi></mi><mi>γ</mi><mi>cos</mi><mi>β</mi></mtd><mtd><mi>sin</mi><mi></mi><mi>γ</mi><mi>cos</mi><mi>β</mi></mtd><mtd><mo>-</mo><mi>sin</mi><mi>β</mi></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><mi></mi><mi>γ</mi><mi>cos</mi><mi>α</mi><mo>+</mo><mi>cos</mi><mi></mi><mi>γ</mi><mi>sin</mi><mi></mi><mi>β</mi><mi>sin</mi><mi>α</mi></mtd><mtd><mi>cos</mi><mi></mi><mi>γ</mi><mi>cos</mi><mi>α</mi><mo>+</mo><mi>sin</mi><mi></mi><mi>γ</mi><mi>sin</mi><mi></mi><mi>β</mi><mi>sin</mi><mi>α</mi></mtd><mtd><mi>cos</mi><mi></mi><mi>β</mi><mi>sin</mi><mi>α</mi></mtd></mtr><mtr><mtd><mi>sin</mi><mi></mi><mi>γ</mi><mi>sin</mi><mi>α</mi><mo>+</mo><mi>cos</mi><mi></mi><mi>γ</mi><mi>sin</mi><mi></mi><mi>β</mi><mi>cos</mi><mi>α</mi></mtd><mtd><mo>-</mo><mi>cos</mi><mi></mi><mi>γ</mi><mi>sin</mi><mi>α</mi><mo>+</mo><mi>sin</mi><mi></mi><mi>γ</mi><mi>sin</mi><mi></mi><mi>β</mi><mi>cos</mi><mi>α</mi></mtd><mtd><mi>cos</mi><mi></mi><mi>β</mi><mi>cos</mi><mi>α</mi></mtd></mtr></mtable></mfenced></mrow>]]></math></maths>其它方位微结构相对参考微结构的力与位移转换方程为<mathsid="math0004"num="0004"><math><![CDATA[<mrow><msub><mi>u</mi><mi>cp</mi></msub><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><msubsup><mi>T</mi><mi>c</mi><mi>p</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>T</mi><mi>c</mi><mi>p</mi></msubsup></mtd></mtr></mtable></mfenced><msub><mi>u</mi><mi>c</mi></msub><mo>,</mo></mrow>]]></math>id="icf0004"file="A2008101901110002C4.tif"wi="25"he="13"top="238"left="42"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0005"num="0005"><math><![CDATA[<mrow><msub><mi>F</mi><mi>cp</mi></msub><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><msubsup><mi>T</mi><mi>c</mi><mi>p</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>T</mi><mi>c</mi><mi>p</mi></msubsup></mtd></mtr></mtable></mfenced><msub><mi>F</mi><mi>c</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math>id="icf0005"file="A2008101901110002C5.tif"wi="97"he="13"top="238"left="71"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>式中ucp＝[ux，uy，uz，ax，ay，az]T，表示与p坐标系固连的微结构在c坐标系的位移与转角，Fcp＝[Fx，Fy，Fz，Tx，Ty，Tz]T，表示微结构所受力及转矩，ux，uy，uz，ax，ay，az分别为沿x、y、z方向的平动位移和绕x、y、z轴的转角，Fx，Fy，Fz，Tx，Ty，Tz分别为沿x、y、z方向的力和绕x、y、z轴的转矩，uc为参考微结构在c坐标系中的位移与转角，Fc为参考微结构在c坐标系中的力与转矩；步骤三将步骤二所得力与位移之间的转换方程(1)代入到步骤一所得的参考微结构的宏模型中，同时根据其它微结构相对参考微结构的方位，设置相对转角参数值(α，β，γ)，代入方程(1)中，提取其它方位微结构的宏模型，实现角度参数化的宏建模。全文摘要本发明公开了一种微机电系统的角度参数化宏建模方法，属于微机电系统设计与模型降阶领域。该方法首先从相同几何形状的微结构中任意选取一个作为参考微结构，采用数值降阶方法提取其宏模型；然后根据微结构方位之间的关系，采用矩阵坐标变换方法建立其它方位微结构相对参考微结构的力与位移之间的转换方程；最后，将力与位移之间的转换方程代入宏模型的行为方程中，实现MEMS的角度参数化宏建模。本发明采用矩阵坐标变换方法实现角度参数化的宏建模，通过设置相应的角度参数，即可避免结构完全相同，初始方位不同的MEMS微结构的重复宏建模过程，从而使基于宏模型的MEMS系统级建模与仿真速度加快，进而提高整个MEMS的设计效率。文档编号B81C99/00GK101554989SQ20081019011公开日2009年10月14日申请日期2008年12月30日优先权日2008年12月30日发明者吕湘连,姜澄宇,常洪龙,徐景辉,苑伟政申请人:西北工业大学