发明名称 微机电系统的角度参数化宏建模方法
摘要 本发明公开了一种微机电系统的角度参数化宏建模方法,属于微机电系统设计与模型降阶领域。该方法首先从相同几何形状的微结构中任意选取一个作为参考微结构,采用数值降阶方法提取其宏模型;然后根据微结构方位之间的关系,采用矩阵坐标变换方法建立其它方位微结构相对参考微结构的力与位移之间的转换方程;最后,将力与位移之间的转换方程代入宏模型的行为方程中,实现MEMS的角度参数化宏建模。本发明采用矩阵坐标变换方法实现角度参数化的宏建模,通过设置相应的角度参数,即可避免结构完全相同,初始方位不同的MEMS微结构的重复宏建模过程,从而使基于宏模型的MEMS系统级建模与仿真速度加快,进而提高整个MEMS的设计效率。
申请公布号 CN101551835A 申请公布日期 2009.10.07
申请号 CN200910134186.9 申请日期 2009.04.15
申请人 西北工业大学 发明人 苑伟政;徐景辉;常洪龙;吕湘连;姜澄宇
分类号 G06F17/50(2006.01)I;B81C5/00(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西北工业大学专利中心 代理人 夏维力
主权项 1.一种微机电系统的角度参数化宏建模方法,其特征在于,包括如下步骤:步骤一:进行参考微结构的宏建模,对于几何形状完全相同,初始方位不同的微结构,任意选择其中一个部件作为参考微结构,建立其宏模型:<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>C</mi><mi>r</mi></msub><msub><mover><mi>x</mi><mo>&CenterDot;</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是时间变量,x<sub>r</sub>(t)是宏模型的状态变量,C<sub>r</sub>、G<sub>r</sub>为宏模型的系统矩阵,B<sub>r</sub>是宏模型的输入矩阵,F<sub>c</sub>为输入载荷向量,u<sub>c</sub>为输出位移向量,L<sub>r</sub>是宏模型的输出矩阵;步骤二:根据微结构方位之间的关系,采用矩阵坐标变换方法建立其它方位微结构相对参考微结构的力与位移转换方程,首先定义两个坐标系以描述微结构在三维空间内的运动,把与参考微结构初始方位重合的坐标系称为总体坐标系,用c-x<sub>c</sub>y<sub>c</sub>z<sub>c</sub>表示,简称c坐标系,把固连于参考微结构某点p的坐标系称为局部坐标系,用p-x<sub>p</sub>y<sub>p</sub>z<sub>p</sub>表示,简称p坐标系,采用以下转动顺序引入确定p坐标系在c坐标系的方位角,即p坐标系绕c坐标系x、y、z坐标轴的相对转角,也即其它微结构相对参考微结构的转角,用(α,β,γ)表示:最初两坐标系完全重合,而后顺序通过三次转动达到p坐标系所在方位,这三次转动为:绕z<sub>c</sub>转一个γ角,使p-x<sub>p</sub>y<sub>p</sub>z<sub>p</sub>由最初与c-x<sub>c</sub>y<sub>c</sub>z<sub>c</sub>重合的位置转到p-x<sub>p1</sub>y<sub>p1</sub>z<sub>c</sub>,第二次绕y<sub>p1</sub>轴转β角,使p-x<sub>p1</sub>y<sub>p1</sub>z<sub>c</sub>到达新的p-x<sub>p2</sub>y<sub>p1</sub>z<sub>p2</sub>位置,第三次绕x<sub>p2</sub>转a角,使p-x<sub>p2</sub>y<sub>p1</sub>z<sub>p2</sub>到达p-x<sub>p</sub>y<sub>p</sub>z<sub>p</sub>的最终位置,把这种转动方式记为<maths num="0002"><![CDATA[<math><mrow><mi>c</mi><munderover><mo>&RightArrow;</mo><msub><mi>z</mi><mi>c</mi></msub><mi>&gamma;</mi></munderover><mn>1</mn><munderover><mo>&RightArrow;</mo><msub><mi>y</mi><mn>1</mn></msub><mi>&beta;</mi></munderover><mn>2</mn><munderover><mo>&RightArrow;</mo><msub><mi>x</mi><mn>2</mn></msub><mi>&alpha;</mi></munderover><mi>p</mi><mo>;</mo></mrow></math>]]></maths>变换过程中c坐标系到p坐标系的转换矩阵T<sub>c</sub><sup>p</sup>为:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>T</mi><mi>c</mi><mi>p</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mi></mi><mi>&gamma;</mi><mi>cos</mi><mi>&beta;</mi></mtd><mtd><mi>sin</mi><mi></mi><mi>&gamma;</mi><mi>cos</mi><mi>&beta;</mi></mtd><mtd><mo>-</mo><mi>sin</mi><mi>&beta;</mi></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><mi></mi><mi>&gamma;</mi><mi>cos</mi><mi>&alpha;</mi><mo>+</mo><mi>cos</mi><mi></mi><mi>&gamma;</mi><mi>sin</mi><mi></mi><mi>&beta;</mi><mi>sin</mi><mi>&alpha;</mi></mtd><mtd><mi>cos</mi><mi></mi><mi>&gamma;</mi><mi>cos</mi><mi>&alpha;</mi><mo>+</mo><mi>sin</mi><mi></mi><mi>&gamma;</mi><mi>sin</mi><mi></mi><mi>&beta;</mi><mi>sin</mi><mi>&alpha;</mi></mtd><mtd><mi>cos</mi><mi></mi><mi>&beta;</mi><mi>sin</mi><mi>&alpha;</mi></mtd></mtr><mtr><mtd><mi>sin</mi><mi></mi><mi>&gamma;</mi><mi>sin</mi><mi>&alpha;</mi><mo>+</mo><mi>cos</mi><mi></mi><mi>&gamma;</mi><mi>sin</mi><mi></mi><mi>&beta;</mi><mi>cos</mi><mi>&alpha;</mi></mtd><mtd><mo>-</mo><mi>cos</mi><mi></mi><mi>&gamma;</mi><mi>sin</mi><mi>&alpha;</mi><mo>+</mo><mi>sin</mi><mi></mi><mi>&gamma;</mi><mi>sin</mi><mi></mi><mi>&beta;</mi><mi>cos</mi><mi>&alpha;</mi></mtd><mtd><mi>cos</mi><mi></mi><mi>&beta;</mi><mi>cos</mi><mi>&alpha;</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其它方位微结构相对参考微结构的力与位移转换方程为:<maths num="0004"><![CDATA[<math><mrow><msub><mi>u</mi><mi>cp</mi></msub><mo>=</mo><mfenced open='[' 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>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>F</mi><mi>cp</mi></msub><mo>=</mo><mfenced open='[' 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>]]></maths>式中u<sub>cp</sub>=[u<sub>x</sub>,u<sub>y</sub>,u<sub>z</sub>,a<sub>x</sub>,a<sub>y</sub>,a<sub>z</sub>]<sup>T</sup>,表示与p坐标系固连的微结构在c坐标系的位移与转角,F<sub>cp</sub>=[F<sub>x</sub>,F<sub>y</sub>,F<sub>z</sub>,T<sub>x</sub>,T<sub>y</sub>,T<sub>z</sub>]<sup>T</sup>,表示微结构所受力及转矩,u<sub>x</sub>,u<sub>y</sub>,u<sub>z</sub>,a<sub>x</sub>,a<sub>y</sub>,a<sub>z</sub>分别为沿x、y、z方向的平动位移和绕x、y、z轴的转角,F<sub>x</sub>,F<sub>y</sub>,F<sub>z</sub>,T<sub>x</sub>,T<sub>y</sub>,T<sub>z</sub>分别为沿x、y、z方向的力和绕x、y、z轴的转矩,u<sub>c</sub>为参考微结构在c坐标系中的位移与转角,F<sub>c</sub>为参考微结构在c坐标系中的力与转矩;步骤三:将步骤二所得力与位移之间的转换方程(1)代入到步骤一所得的参考微结构的宏模型中,同时根据其它微结构相对参考微结构的方位,设置相对转角参数值(α,β,γ),代入方程(1)中,提取其它方位微结构的宏模型,实现角度参数化的宏建模。
地址 710072陕西省西安市友谊西路127号