发明名称 MEMS系统级设计中一种等腰梯形截面梁的建模方法
摘要 本发明公开了一种微机电系统(MEMS)系统级设计中的等腰梯形截面梁建模方法,属于微机电系统设计领域。该方法的主要过程为:首先建立等腰梯形截面梁的局部坐标系,然后根据结构力学和结构矩阵分析理论推导出局部坐标系下等腰梯形截面梁刚度矩阵、阻尼矩阵和质量矩阵,并利用这些矩阵建立等腰梯形截面梁的二阶动力学方程,最后利用硬件描述语言对二阶动力学方程进行编码,实现等腰梯形截面梁的参数化系统级组件建模。本发明提出的方法解决了当考虑加工误差时系统级无法对微梁结构精确仿真的问题,同时也实现了对具有等腰梯形截面梁结构的MEMS器件进行系统级仿真设计的功能。
申请公布号 CN101694673A 申请公布日期 2010.04.14
申请号 CN200910024394.3 申请日期 2009.10.20
申请人 西北工业大学 发明人 苑伟政;郝星;常洪龙
分类号 G06F17/50(2006.01)I;B81C5/00(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 夏维力
主权项 1.一种MEMS系统级设计中等腰梯形截面梁的建模方法,其特征在于,包括如下步骤:步骤一:建立等腰梯形截面梁的三维局部坐标系xyz;已知等腰梯形截面梁的弹性模量E,泊松比为υ,密度为ρ,抗扭量为J<sub>i</sub>,绕x轴扭转和沿x轴平动的梁单位长度上阻尼系数为B<sub>x</sub>,绕y轴转动和沿y轴平动的梁单位长度上阻尼系数为B<sub>y</sub>,绕z轴转动和沿z轴平动的梁单位长度上阻尼系数为B<sub>z</sub>,并以等腰梯形截面梁的长度方向为x轴,宽度方向为y轴,厚度方向为z轴建立局部坐标系xyz,l为等腰梯形截面梁长度,t为厚度,w<sub>1</sub>为截面上底的宽度,w<sub>2</sub>为截面下底的宽度,x<sub>C</sub>y<sub>C</sub>z<sub>C</sub>为全局坐标系;步骤二:计算等腰梯形截面梁的主惯性矩I<sub>y</sub>、I<sub>z</sub>:<maths num="0001"><![CDATA[<math><mrow><msub><mi>I</mi><mi>y</mi></msub><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><msubsup><mi>w</mi><mn>1</mn><mn>3</mn></msubsup><mo>+</mo><msub><mrow><mn>4</mn><mi>w</mi></mrow><mn>1</mn></msub><msub><mi>w</mi><mn>2</mn></msub><mo>+</mo><msubsup><mi>w</mi><mn>2</mn><mn>2</mn></msubsup><mo>)</mo></mrow><msup><mi>t</mi><mn>3</mn></msup></mrow><mrow><mn>36</mn><mrow><mo>(</mo><msub><mi>w</mi><mn>1</mn></msub><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>I</mi><mi>z</mi></msub><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><msubsup><mi>w</mi><mn>1</mn><mn>3</mn></msubsup><mo>+</mo><msubsup><mi>w</mi><mn>1</mn><mn>2</mn></msubsup><msub><mi>w</mi><mn>2</mn></msub><mo>+</mo><msub><mi>w</mi><mn>1</mn></msub><msubsup><mi>w</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>w</mi><mn>2</mn><mn>3</mn></msubsup><mo>)</mo></mrow><mi>t</mi></mrow><mn>48</mn></mfrac></mrow></math>]]></maths>步骤三:计算局部坐标系下等腰梯形截面梁的刚度矩阵K、质量矩阵M:<img file="F2009100243943C0000013.GIF" wi="800" he="564" /><img file="F2009100243943C0000021.GIF" wi="800" he="515" />其中G=E/(2(1+υ))为剪切模量,I<sub>p</sub>=I<sub>y</sub>+I<sub>z</sub>,A为等腰梯形截面梁的截面积;步骤四:计算局部坐标系下等腰梯形截面梁的阻尼矩阵B:<maths num="0003"><![CDATA[<math><mrow><mi>B</mi><mo>=</mo><mfrac><mi>M</mi><mi>&rho;Al</mi></mfrac><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>B</mi><mi>x</mi></msub></mtd><mtd><msub><mi>B</mi><mi>y</mi></msub></mtd><mtd><msub><mi>B</mi><mi>z</mi></msub></mtd><mtd><msub><mi>B</mi><mi>x</mi></msub></mtd><mtd><msub><mi>B</mi><mi>y</mi></msub></mtd><mtd><msub><mi>B</mi><mi>z</mi></msub></mtd><mtd><msub><mi>B</mi><mi>x</mi></msub></mtd><mtd><msub><mi>B</mi><mi>y</mi></msub></mtd><mtd><msub><mi>B</mi><mi>z</mi></msub></mtd><mtd><msub><mi>B</mi><mi>x</mi></msub></mtd><mtd><msub><mi>B</mi><mi>y</mi></msub></mtd><mtd><msub><mi>B</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow></math>]]></maths>步骤五:根据结构矩阵分析理论将局部坐标系下等腰梯形截面梁的刚度矩阵K、阻尼矩阵B和质量矩阵M,转换成全局坐标系下的刚度矩阵K<sub>C</sub>、阻尼矩阵B<sub>C</sub>和质量矩阵M<sub>C</sub>:<maths num="0004"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>K</mi><mi>C</mi></msub><mo>=</mo><msup><mi>&lambda;</mi><mi>T</mi></msup><mi>K&lambda;</mi></mtd></mtr><mtr><mtd><msub><mi>B</mi><mi>C</mi></msub><mo>=</mo><msup><mi>&lambda;</mi><mi>T</mi></msup><mi>B&lambda;</mi></mtd></mtr><mtr><mtd><msub><mi>M</mi><mi>C</mi></msub><mo>=</mo><msup><mi>&lambda;</mi><mi>T</mi></msup><mi>M&lambda;</mi></mtd></mtr></mtable></mfenced></math>]]></maths>式中λ为局部坐标系向总体坐标系投影的坐标变换矩阵,其表达式为<maths num="0005"><![CDATA[<math><mrow><mi>&lambda;</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mi></mi><mi>&lambda;</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>其中α、β、γ分别为局部坐标系x、y、z轴与全局坐标系x<sub>C、</sub>y<sub>C</sub>、z<sub>C</sub>轴之间的夹角;步骤六:利用全局坐标系下的刚度矩阵K<sub>C</sub>、阻尼矩阵B<sub>C</sub>和质量矩阵M<sub>C</sub>建立结构的二阶动力方程<maths num="0006"><![CDATA[<math><mrow><msub><mi>K</mi><mi>C</mi></msub><msub><mi>X</mi><mi>C</mi></msub><mo>+</mo><msub><mi>B</mi><mi>C</mi></msub><msub><mover><mi>X</mi><mo>&CenterDot;</mo></mover><mi>C</mi></msub><mo>+</mo><msub><mi>M</mi><mi>C</mi></msub><msub><mover><mi>X</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mi>C</mi></msub><mo>=</mo><msub><mi>F</mi><mi>C</mi></msub></mrow></math>]]></maths>其中X<sub>C</sub>、F<sub>C</sub>为全局坐标系下梁节点产生的位移和所受的合外力;步骤七:应用硬件描述语言对步骤六建立的二阶动力方程编码,即得等腰梯形截面梁的系统级组件模型。
地址 710072 陕西省西安市友谊西路127号