发明名称 一种声表面波器件的建模方法
摘要 本发明公开了一种声表面波器件的建模方法。对于声表面波器件,该方法采用曲六面体进行剖分,因而可以实现任意网格建模。通过求解克里斯托夫方程,可以满足晶体介质的压电性,描述压电晶体内部的耦合场,并准确地求得介质中的电场分布和声场分布。在求解过程中,对声表面波器件输入端强加冲激脉冲,将器件输出端的电势信息通过做差转化为电压输出,获得系统冲击响应,并进行傅立叶变换,最终得到描述器件特性的频率响应。本发明利用时域谱元法能精确求解偏微分方程的特性,同时采用了高效的计算机算法,能实现声表面波器件的高精度建模与仿真。
申请公布号 CN105022850A 申请公布日期 2015.11.04
申请号 CN201410173063.7 申请日期 2014.04.25
申请人 南京理工大学 发明人 陈如山;丁大志;樊振宏;盛亦军;顾闻阳
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 朱显国
主权项 一种声表面波器件的建模方法,其特征在于步骤如下:第一步,建立声表面波器件的三维模型,并采用曲六面体对声表面波器件整体剖分,得到器件的网格信息,包括六面体的个数和每个六面体上的点的坐标;第二步,将电势值和位移值定义在六面体的每个点上,并用时域谱元法中的GLL多项式作为基函数对电势和位移进行展开,对电势的展开式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&Phi;</mi><mo>=</mo><msub><mi>N</mi><mi>&Phi;</mi></msub><mover><mi>&Phi;</mi><mo>^</mo></mover><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>N</mi><mn>1</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>N</mi><mn>27</mn></msub></mtd></mtr></mtable></mfenced><mrow><mo>(</mo><msub><mi>&Phi;</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000496646790000011.GIF" wi="659" he="87" /></maths>其中,N<sub>i</sub>即为展开基函数,(Φ<sub>i</sub>)为六面体元胞内各个GLL点上的电势值;将声表面波器件内部的位移在XYZ三个方向进行展开,对位移的展开式如下:<img file="FDA0000496646790000015.GIF" wi="1284" he="236" />其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>u</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mrow><mn>3</mn><mi>i</mi></mrow></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000496646790000013.GIF" wi="123" he="233" /></maths>为六面体元胞内各个GLL点上的XYZ三个方向的位移值;第三步,将展开后的电势和位移代入克里斯托夫方程,采用哈密顿变分原理,填充系数矩阵,得到矩阵方程;第四步,对声表面波器件输入端强加冲激脉冲,脉冲的表达形式如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>V</mi><mi>i</mi></msub><mfenced open='{' close=''><mtable><mtr><mtd><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mn>0</mn></msub><mi>sin</mi><mrow><mo>(</mo><mn>4</mn><mi>&pi;</mi><msub><mi>f</mi><mn>0</mn></msub><mi>t</mi><mo>)</mo></mrow></mtd><mtd><mi>if</mi><mn>0</mn><mo>&lt;</mo><mi>t</mi><mo>&le;</mo><mn>1</mn><mo>/</mo><mn>4</mn><msub><mi>f</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>ift</mi><mo>></mo><mn>1</mn><mo>/</mo><mn>4</mn><msub><mi>f</mi><mn>0</mn></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000496646790000014.GIF" wi="836" he="159" /></maths>其中,f<sub>0</sub>为声表面波器件的特征频率;根据第三步中得到的矩阵方程,求得输出端时域电势信息,进而提取声表面波器件的冲击响应,对冲击响应进行傅立叶变换,即可获得对应的频率响应。
地址 210094 江苏省南京市孝陵卫200号