CN113962093A - 一种基于谱元法的声表面波滤波器设计方法 - Google Patents
一种基于谱元法的声表面波滤波器设计方法 Download PDFInfo
- Publication number
- CN113962093A CN113962093A CN202111243054.7A CN202111243054A CN113962093A CN 113962093 A CN113962093 A CN 113962093A CN 202111243054 A CN202111243054 A CN 202111243054A CN 113962093 A CN113962093 A CN 113962093A
- Authority
- CN
- China
- Prior art keywords
- acoustic wave
- surface acoustic
- wave filter
- filter model
- model
- 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.)
- Pending
Links
- 238000010897 surface acoustic wave method Methods 0.000 title claims abstract description 129
- 238000000034 method Methods 0.000 title claims abstract description 88
- 230000003595 spectral effect Effects 0.000 title claims abstract description 54
- 238000013461 design Methods 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 66
- 238000003780 insertion Methods 0.000 claims abstract description 29
- 230000037431 insertion Effects 0.000 claims abstract description 29
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 22
- 230000008878 coupling Effects 0.000 claims abstract description 20
- 238000010168 coupling process Methods 0.000 claims abstract description 20
- 238000005859 coupling reaction Methods 0.000 claims abstract description 20
- 230000005284 excitation Effects 0.000 claims abstract description 11
- 239000000758 substrate Substances 0.000 claims description 25
- 239000000463 material Substances 0.000 claims description 20
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 230000005686 electrostatic field Effects 0.000 claims description 7
- 239000007787 solid Substances 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 6
- 238000010521 absorption reaction Methods 0.000 claims description 4
- 239000004020 conductor Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 11
- 238000004891 communication Methods 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- XLOMVQKBTHCTTD-UHFFFAOYSA-N Zinc monoxide Chemical compound [Zn]=O XLOMVQKBTHCTTD-UHFFFAOYSA-N 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000005684 electric field Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000000306 component Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000011787 zinc oxide Substances 0.000 description 2
- 239000004411 aluminium Substances 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011438 discrete method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005669 field effect Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Surface Acoustic Wave Elements And Circuit Networks Thereof (AREA)
Abstract
本发明提出了一种基于谱元法的声表面波滤波器设计方法，包括：对声表面波滤波器进行几何建模，得到声表面波滤波器模型；设置所述模型的参数以及激励信息，以确定所述模型的计算区域；对所述模型进行离散化处理；将离散后的所述模型导入谱元算法中，并对所述模型的所述计算区域进行多物理场的耦合；在所述模型的终端添加输入电压，利用谱元算法计算公式在不同频段下对所述模型进行求解，得到所述声表面波滤波器在每个频段下的输出电压，并计算所述声表面波滤波器的插入损耗。本发明基于谱元法的声表面波滤波器设计方法，在计算内存小的前提下能够极大的提高计算效率，并且计算结果的稳定性好。
Description
技术领域
本发明涉及5G射频滤波器的技术领域，具体涉及一种基于谱元法的声表面波滤波器设计方法。
背景技术
随着5G通信、物联网和智能汽车的发展为射频器件行业带来新机，未来几年渗透率将快速提升，且射频前端行业高速增长趋势明确。SAW(声表面波)射频滤波器作为5G通信系统的核心组件，它主要功能是解决不同频段和通信系统之间产生的信号干扰问题。SAW滤波器的工作原理是：当输入电信号激励时，压电材料由于压电效应会引发弹性形变，引起固体质点的振动，形成沿固体表面的传播的声表面波。当表面波传到输出电极时，由于逆压电效应，会在SAW滤波器的输出端产生电信号，该过程将实现射频信号的过滤。然而，随着5G通信行业的快速发展，移动设备对SAW滤波器的要求尺寸、滤波的频段也越来越高。因此准确模拟出尺寸小、高频段的SAW滤波器对通信领域具有重要意义。
通过数值仿真来设计SAW滤波器是最为直接的方法，常用的方法有：δ函数模型、等效电路模型、耦合模(COM)理论、P矩阵方法和有限元方法。其中，δ函数模型和等效电路模型对滤波器模型进行了简化，忽略了反射、电极阻抗等特性；COM理论考虑了电极厚度，波速色散和衰减损耗。该理论具有明确的物理意义，同时也能考虑寄生效应和非线性效应；P矩阵是COM理论的扩展形式，将入射波和反射波振幅、激励电压和电流用矩阵的形式联立求解；随着三维数值计算的发展，利用有限元方法直接求解压电方程，不仅能准确得到器件中场分布，同时也考虑各种复杂的几何模型，仿真精度和可靠性高，也是目前主流的计算方法，但计算复杂度远大于上述其它方法。随着滤波器频率的不断提高,传统的近似理论已无法满足计算精度需要，且三维有限元方法计算效率低，需要开发更高效的三维全波仿真算法，实现对SAW滤波器的快速设计。
发明内容
为了解决现有技术中的有限元算法对于高频声表面波滤波器的仿真中计算量过大、计算时间过长、限制模型参数设置等技术问题，本发明提出了一种基于谱元法的声表面波滤波器设计方法。
本申请提出了一种基于谱元法的声表面波滤波器设计方法，包括：
S1、对声表面波滤波器进行几何建模，得到声表面波滤波器模型；
S2、设置所述声表面波滤波器模型的参数以及激励信息，以确定所述声表面波滤波器模型的计算区域；
S3、对所述声表面波滤波器模型进行离散化处理；
S4、将离散后的所述声表面波滤波器模型导入谱元算法中，并对所述声表面波滤波器模型的所述计算区域进行多物理场的耦合；
S5、在所述声表面波滤波器模型的终端添加输入电压，利用谱元算法计算公式在不同频段下对所述声表面波滤波器模型进行求解，得到所述声表面波滤波器在每个频段下的输出电压，并计算所述声表面波滤波器的插入损耗。
通过建立声表面波滤波器模型，对声表面波滤波器模型设置参数和激励信息，以及添加多物理场进行耦合，再基于谱元法对声表面波滤波器模型的计算区域进行数值仿真计算，计算出声表面波滤波器的插入损耗。这种基于谱元法的声表面波滤波器的设计方法，在计算内存小的前提下能够极大的提高计算效率，并且计算结果的稳定性好。
优选的，所述声表面波滤波器模型的所述计算区域包括由压电材料构成的衬底和金属导电材料构成的电极。
优选的，所述步骤S2中设置所述声表面波滤波器模型的参数以及激励信息具体包括：
设定所述衬底的长度、宽度、厚度以及材料属性；
对所述衬底的四周及底面施加吸收边界条件，对所述衬底的上表面施加自由边界条件；
确定所述电极的长度、宽度、厚度、材料属性以及数量。
进一步优选的，所述电极的数量为15个。声表面波滤波器的电极数量越多，插入损耗在一定程度上也会有所减小降低，但计算量也会更大，因此将电极的数量设置成15个，既使得插入损耗较低，同时也避免其计算量过大，从而方便通过其它方法来对声表面波滤波器模型进行对比仿真计算，以验证本申请设计方法的可靠性。
优选的，所述步骤S3具体为：通过COMSOL软件将所述声表面波滤波器模型离散成多个六面体网格单元。通过对声表面波滤波器模型进行离散，从而方便后续利用谱元算法对声表面波滤波器模型进行求解。
优选的，所述步骤S4中对所述声表面波滤波器模型的所述计算区域进行多物理场的耦合具体包括：
在所述电极上添加静电场；
在整个所述声表面波滤波器模型上添加固体力学场。
通过在声表面波滤波器模型的计算区域上添加多物理场，从而对声表面波滤波器模型进行多物理场的耦合。
优选的，所述步骤S4中对所述声表面波滤波器模型的所述计算区域进行多物理场的耦合具体还包括：将所述衬底设置为压电区域，并在所述压电区域上添加模拟损耗。通过添加模拟损耗，进一步提高仿真结果的准确性。
进一步优选的，所述模拟损耗包括机械损耗、耦合损耗和介电损耗。
优选的，所述步骤S5中的谱元算法计算公式具体为：
其中，Ω为整个计算区域，Γ为边界面，n为边界面法向，w为测试函数，D为电位移矢量，ρe为外加电荷密度，n·D表示输出电压。
优选的，所述步骤S5中的所述插入损耗的计算公式具体为：
其中，IL表示插入损耗，Vin是输入电压，Vout是输出电压。
本申请提出了一种基于谱元法的声表面波滤波器设计方法，通过建立声表面波滤波器模型，对声表面波滤波器模型设置参数和激励信息以确定计算区域，然后对声表面波滤波器模型进行离散并导入谱元算法中，在声表面波滤波器模型的计算区域上添加多物理场和模拟损耗，利用谱元算法对声表面波滤波器模型的计算区域进行耦合计算，得到声表面波滤波器在不同频率下的输出电压，最后通过插入损耗计算公式计算出声表面波滤波器的插入损耗。该设计方法采用基于谱元算法的声表面波滤波器设计方法，在计算内存小的前提下能够极大的提高计算效率，且计算结果的稳定性较好。
附图说明
包括附图以提供对实施例的进一步理解并且附图被并入本说明书中并且构成本说明书的一部分。附图图示了实施例并且与描述一起用于解释本发明的原理。将容易认识到其它实施例和实施例的很多预期优点，因为通过引用以下详细描述，它们变得被更好地理解。附图的元件不一定是相互按照比例的。同样的附图标记指代对应的类似部件。
图1是根据本发明实施例的基于谱元法的声表面波滤波器设计方法的流程图；
图2是根据本发明一个具体实施例的声表面波滤波器模型的结构示意图；
图3是根据本发明一个具体实施例的声表面波滤波器模型的网格剖分图；
图4是根据本发明一个具体实施例的声表面波滤波器在300MHz到400MHz之间使用谱元法与COMSOL得到插入损耗对比图；
图5是根据本发明另一个具体实施例的声表面波滤波器在650MHz到800MHz之间使用谱元法得到的插入损耗结果图。
附图标记说明：1、衬底；2、电极。
具体实施方式
下面将详细描述本发明的各个方面的特征和示例性实施例，为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细描述。应理解，此处所描述的具体实施例仅被配置为解释本发明，并不被配置为限定本发明。对于本领域技术人员来说，本发明可以在不需要这些具体细节中的一些细节的情况下实施。下面对实施例的描述仅仅是为了通过示出本发明的示例来提供对本发明更好的理解。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括......”限定的要素，并不排除在包括要素的过程、方法、物品或者设备中还存在另外的相同要素。
图1示出了根据本发明实施例的基于谱元法的声表面波滤波器设计方法的流程图，如图1所示，基于谱元法的声表面波滤波器设计方法包括以下步骤：
S1、对声表面波滤波器进行几何建模，得到声表面波滤波器模型。
通过建立声表面波滤波器模型，实现后续对声表面波滤波器的数值仿真。
S2、设置声表面波滤波器模型的参数以及激励信息，以确定声表面波滤波器模型的计算区域。
图2示出了根据本发明一个具体实施例的声表面波滤波器模型的结构示意图，如图2所示，声表面波滤波器模型的计算区域包括由压电材料构成的衬底1和金属导电材料构成的电极2的组成。因此，所设置的声表面波滤波器模型的参数包括：设定衬底1的长度、宽度、厚度以及材料属性；设定电极2的长度、宽度、厚度、材料属性以及数量。所设置的激励信息包括：对衬底1施加边界条件。
在优选的实施例中，衬底1的长度、宽度、厚度分别为272um、60um、12um，衬底1的材料具体为ZnO(氧化锌)。电极2的长度、宽度、厚度分别为24um、2um、0.24um，电极2的材料为Al(铝)，电极2的数量为15个。
需要说明的是，除了衬底1的尺寸、材料属性以及电极2的尺寸、材料属性和数量会影响到声表面波滤波器的插入损耗值外，电极2的叉指与叉指间的间距以及声表面波滤波器的孔径同样会影响到声表面波滤波器的插入损耗值。本实施例中只对衬底1的尺寸、材料属性以及电极2的尺寸、材料属性和数量进行限定，但本领域技术人员应该认识到，在其它实施方式中，对电极2的叉指与叉指间的间距以及声表面波滤波器的孔径作简单限定的，仍在本申请的保护范围之内。
在优选的实施例中，所设置的声表面波滤波器模型的激励信息为固定力学边界条件。根据声表面波滤波器的应用，固定力学边界条件可分为三种类型：1、自由边界：n·τ|Γ＝0，即声表面波滤波器表面的牵引力为零；2、固定边界条件：u|Γ＝0，即声表面波滤波器表面的位移场为零；3、吸收边界条件：τ|Γ＝ρcpvn+ρcsvt，即模拟无限大的声表面波滤波器衬底模型。本实施例中，对衬底1的四周及底面施加吸收边界条件，对衬底1的上表面施加自由边界条件。
需要说明的是，本实施例中只针对衬底1添加固定力学边界条件，在其它实施方式中，同样可以根据声表面波滤波器电极2的连接方式的不同，在电极2上添加静电场方程边界条件。根据声表面波滤波器电极的连接方式的不同，静电场方程边界条件可分为五种类型：1、零电荷边界，即电极上无电流通过；2、理想电压源，即电压源直接加载到电极上；3、理想电流源，即电流源直接加载到电极上；4、带电阻的电压源，即电压源存在电阻；5、带电阻的电流源，即电流源存在电阻。
S3、对声表面波滤波器模型进行离散化处理。
图3示出了根据本发明一个具体实施例的声表面波滤波器模型的网格剖分图，如图3所示，在具体的实施例中，通过COMSOL软件对声表面波滤波器模型进行离散，从而将声表面波滤波器模型剖分成多个六面体网格单元。由于谱元算法本身就是一种离散方法，因此将声表面波滤波器模型离散以便后续通过谱元算法对声表面波滤波器模型进行计算。在其它实施方式中。也可通过Trelis等其它软件对声表面波滤波器模型进行离散，这里不作进一步限定。
S4、将离散后的声表面波滤波器模型导入谱元算法中，并对声表面波滤波器模型的计算区域进行多物理场的耦合。
在具体的实施例中，在谱元算法程序中对电极上添加静电场，对整个声表面波滤波器模型上添加固体力学场，并将衬底设置成压电区域，并在压电区域上添加机械损耗、耦合损耗、介电损耗等三种模拟损耗，以进一步的提高仿真的准确性。
S5、在声表面波滤波器模型的终端添加输入电压，利用谱元算法计算公式在不同频段下对声表面波滤波器模型进行求解，得到声表面波滤波器在每个频段下的输出电压，并计算声表面波滤波器的插入损耗。
在具体的实施例中，采用频域谱元法对声表面波滤波器模型进行求解。时域谱元法相较于频域谱元法来说，虽然时域谱元法能通过傅里叶变换得到不同频率的信号，而频域谱元法则需要单独计算每个频点。但频域谱元法在处理复杂边界条件(如PML、ABC)、外接RCL复杂电路模型、以及为各向异性材料添加衰减具有明显优势。因为频域方法能直接处理复数媒质和边界，而时域谱元法则需要通过复杂的卷积运算才能够实现。相比较而言，频域谱元法在实现这些功能时更为简便，也能够极大的降低模型的计算时间和空间复杂度。
具体的，声表面波滤波器是基于电声多物理场效应实现电与声信号之间的相互转化。声信号和电信号的传播分别满足固体力学方程和静电场方程：
其中，式(1.1)、式(1.2)中的ρ为压电材料密度，ω为滤波工作频率，u为发生形变的位移矢量，τ为材料的应力张量，D为电位移矢量，ρe为外加电荷密度。在τ和D的本构关系下，将实现电声信号的相互耦合：
D＝e∶ε+∈·E (1.3)
τ＝c∶ε-e·E (1.4)
其中，式(1.3)、式(1.4)中的c为四阶弹性系数矩阵，e为三阶压电耦合系数，∈为二阶介电常数，ε为应变张量，E为电场强度。
式(1.3)为压电效应，即对物体施加应力场将产生电信号；式(1.4)为逆压电效应，即对材料施加电场将产生声信号。电场和应变张量与电势和位移场的关系如下：
其中，V为声表面波滤波器的电压场。因此，方程(1.1)–(1.6)实现了电声信号的转化关系。
而频域谱元法是基于加权余量法来计算声表面波滤波器方程，需要将式(1.1)和(1.2)转化为弱形式，并在式的两端同时乘上测试函数(w)，并对其进行积分，再利用分部积分，得：
其中，式(2.1)和(2.2)中的Ω为整个计算区域，Γ为边界面，n为边界面法向，n·D表示输出电压。将声表面波滤波器的计算区域Ω剖分为K个不重叠的六面体网格单元，再将物理空间上的六面体网格单元通过坐标变换映射到参考单元中，就可直接在参考单元中进行数值积分，求导和插值。本实施例中的谱元法使用Legrange函数作为基函数，并且选取Gauss-Lobatto-Legendre(GLL)点作为插值点。在参考单元中，N阶基函数，具有N+1个GLL点，定义为：
其中，ξ∈[-1,1]，i＝1,2,…,N+1,LN(ξi)和L′N(ξi)分别是N阶Legrange多项式和一阶导数。因此在三维问题中，物理场将利用基函数就行展开，例如位移场x分量为：
这样一来，在已知K个不重叠的六面体网格单元的坐标以及输入电压的情况下，采用谱元法对K个六面体网格单元的计算区域进行离散计算，便可获得对应的系数矩阵A，以及对应输入电压的矩阵B，再利用直接法或迭代法求解线性方程组(Ax＝b)，即可得到整个模型位移场和电压场，从而计算出输出电压。
在得到输出电压后，利用插入损耗计算公式计算出声表面波滤波器的插入损耗。
在具体的实施例中，插入损耗计算公式具体为：
其中，IL表示插入损耗，Vin是输入电压，Vout是输出电压。
在具体的实施例中，为了验证本设计方法的准确性，本实施例中还将声表面波滤波器模型导入到其它软件中进行求解，从而形成参照对比。由于在步骤S3中是通过COMSOL软件对声表面波滤波器模型进行离散的，因此在这里同样采用COMSOL软件进行对比仿真，保证了两种计算方法具有相同的未知量。具体的，设定声表面波滤波器模型的电极数量为15个，在声表面波滤波器模型的终端添加1V的输入电压，通过COMSOL软件和谱元算法分别对300MHz到400MHz下的声表面波滤波器模型进行数值仿真。
需要说明的是，由于COMSOL软件的计算量有限，因此经过多次实验之后，优选电极的数量为15个，保证COMSOL软件能够对声表面波滤波器模型进行计算的同时，声表面波滤波器模型的插入损耗值能够较低。
图4示出了根据本发明一个具体实施例的声表面波滤波器在300MHz到400MHz之间使用谱元法与COMSOL得到插入损耗对比图，如图4所示，该声表面波滤波器可以对300MHz到400MHz这个频段的信号进行过滤，插入损耗值能够达到-22dB左右，并且谱元法(SEM)与COMSOL软件仿真出来的插入损耗的相对误差仅为3.77％，从而证明了基于谱元算法的声表面波滤波器设计方法的可靠性。
进一步的，对于相同的模型和计算结点，测试后发现谱元法代码在计算时间上要比COMSOL有限元法快2倍以上，计算内存为COMSOL的一半不到，证明了谱元法在计算内存小的前提下能够极大的提高计算效率。具体测试参数如表1所示：
| 算法 | 自由度 | 计算时间 | 计算内存 |
| 谱元法 | 500万 | 10小时 | 148G |
| 有限元法 | 500万 | 21.2小时 | 355G |
表1
此外，本实施例中还基于谱元法对650MHz到800MHz这个频段的声表面波滤波器模型进行仿真，该声表面波滤波器模型的衬底的长度、宽度、厚度分别为183um、29um、5.7um，电极的长度、宽度、厚度分别为11.4um、0.95um、0.114um，电极的个数为22个。图5示出了根据本发明另一个具体实施例的声表面波滤波器在650MHz到800MHz之间使用谱元法得到的插入损耗结果图，如图5所示，该声表面波滤波器能够实现对690MHz到720MHz这个频段的信号进行过滤，通带宽度能够达到30MHz，插入损耗的值能够达到-21dB。
本发明提出了一种基于谱元法的声表面波滤波器设计方法，通过建立声表面波滤波器模型，设置声表面波滤波器模型的参数信息和激励信息，对声表面波滤波器模型进行离散并导入谱元算法程序中，对声表面波滤波器模型设置静电场和固体力学场，将声表面波滤波器模型的衬底设置为压电区域，为压电区域添加模拟损耗，然后利用谱元算法计算公式对整个声表面波滤波器模型的计算区域进行多物理场的耦合计算，得出该声表面波滤波器在不同频率下的输出电压。最后通过插入损耗计算公式算出该声表面波滤波器的插入损耗。本发明提供了基于谱元法的声表面波滤波器设计方法，相比于传统的COMSOL，计算时间快2倍以上，且计算内存不到COMSOL的一半，从而在计算内存小的前提下能够极大的提高计算效率，并且计算结果的稳定性较好。
在本申请实施例中，应该理解到，所揭露的技术内容，可通过其它的方式实现。其中，以上所描述的装置/系统/方法实施例仅仅是示意性的，例如所述单元的划分，可以为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，单元或模块的间接耦合或通信连接，可以是电性或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
显然，本领域技术人员在不偏离本发明的精神和范围的情况下可以作出对本发明的实施例的各种修改和改变。以该方式，如果这些修改和改变处于本发明的权利要求及其等同形式的范围内，则本发明还旨在涵盖这些修改和改变。词语“包括”不排除未在权利要求中列出的其它元件或步骤的存在。某些措施记载在相互不同的从属权利要求中的简单事实不表明这些措施的组合不能被用于获利。权利要求中的任何附图标记不应当被认为限制范围。
Claims (10)
1.一种基于谱元法的声表面波滤波器设计方法，其特征在于，包括以下步骤：
S1、对声表面波滤波器进行几何建模，得到声表面波滤波器模型；
S2、设置所述声表面波滤波器模型的参数以及激励信息，以确定所述声表面波滤波器模型的计算区域；
S3、对所述声表面波滤波器模型进行离散化处理；
S4、将离散后的所述声表面波滤波器模型导入谱元算法中，并对所述声表面波滤波器模型的所述计算区域进行多物理场的耦合；
S5、在所述声表面波滤波器模型的终端添加输入电压，利用谱元算法计算公式在不同频段下对所述声表面波滤波器模型进行求解，得到所述声表面波滤波器在每个频段下的输出电压，并计算所述声表面波滤波器的插入损耗。
2.根据权利要求1所述的方法，其特征在于，所述声表面波滤波器模型的所述计算区域包括由压电材料构成的衬底和金属导电材料构成的电极。
3.根据权利要求2所述的方法，其特征在于，所述步骤S2中设置所述声表面波滤波器模型的参数以及激励信息具体包括：
设定所述衬底的长度、宽度、厚度以及材料属性；
对所述衬底的四周及底面施加吸收边界条件，对所述衬底的上表面施加自由边界条件；
确定所述电极的长度、宽度、厚度、材料属性以及数量。
4.根据权利要求3所述的方法，其特征在于，所述电极的数量为15个。
5.根据权利要求1所述的方法，其特征在于，所述步骤S3具体为：通过COMSOL软件将所述声表面波滤波器模型离散成多个六面体网格单元。
6.根据权利要求2所述的方法，其特征在于，所述步骤S4中对所述声表面波滤波器模型的所述计算区域进行多物理场的耦合具体包括：
在所述电极上添加静电场；
在整个所述声表面波滤波器模型上添加固体力学场。
7.根据权利要求2所述的方法，其特征在于，所述步骤S4中对所述声表面波滤波器模型的所述计算区域进行多物理场的耦合具体还包括：将所述衬底设置为压电区域，并在所述压电区域上添加模拟损耗。
8.根据权利要求7所述的方法，其特征在于，所述模拟损耗包括机械损耗、耦合损耗和介电损耗。
9.根据权利要求1所述的方法，其特征在于，所述步骤S5中的谱元算法计算公式具体为：
其中，Ω为整个计算区域，Γ为边界面，n为边界面法向，w为测试函数，D为电位移矢量，ρe为外加电荷密度，n·D表示输出电压。
10.根据权利要求9所述的方法，其特征在于，所述步骤S5中的所述插入损耗的计算公式具体为：
其中，IL表示插入损耗，Vin是输入电压，Vout是输出电压。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111243054.7A CN113962093A (zh) | 2021-10-25 | 2021-10-25 | 一种基于谱元法的声表面波滤波器设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111243054.7A CN113962093A (zh) | 2021-10-25 | 2021-10-25 | 一种基于谱元法的声表面波滤波器设计方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113962093A true CN113962093A (zh) | 2022-01-21 |
Family
ID=79466881
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111243054.7A Pending CN113962093A (zh) | 2021-10-25 | 2021-10-25 | 一种基于谱元法的声表面波滤波器设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113962093A (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050273996A1 (en) * | 2004-06-15 | 2005-12-15 | Yuji Mitsui | Method for manufacturing surface acoustic wave element, as well as surface acoustic wave element manufactured by the same method |
| CN105022850A (zh) * | 2014-04-25 | 2015-11-04 | 南京理工大学 | 一种声表面波器件的建模方法 |
| CN112532202A (zh) * | 2020-11-25 | 2021-03-19 | 重庆邮电大学 | 一种基于bvd模型的saw谐振器及滤波器的设计方法 |
| CN113228628A (zh) * | 2015-08-20 | 2021-08-06 | 高通股份有限公司 | 用于转换非拜尔图案颜色滤波器阵列图像数据的系统和方法 |
| CN113361218A (zh) * | 2021-06-04 | 2021-09-07 | 安徽大学 | 一种基于多物理场耦合建模的体声波滤波器及其设计方法 |
-
2021
- 2021-10-25 CN CN202111243054.7A patent/CN113962093A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050273996A1 (en) * | 2004-06-15 | 2005-12-15 | Yuji Mitsui | Method for manufacturing surface acoustic wave element, as well as surface acoustic wave element manufactured by the same method |
| CN105022850A (zh) * | 2014-04-25 | 2015-11-04 | 南京理工大学 | 一种声表面波器件的建模方法 |
| CN113228628A (zh) * | 2015-08-20 | 2021-08-06 | 高通股份有限公司 | 用于转换非拜尔图案颜色滤波器阵列图像数据的系统和方法 |
| CN112532202A (zh) * | 2020-11-25 | 2021-03-19 | 重庆邮电大学 | 一种基于bvd模型的saw谐振器及滤波器的设计方法 |
| CN113361218A (zh) * | 2021-06-04 | 2021-09-07 | 安徽大学 | 一种基于多物理场耦合建模的体声波滤波器及其设计方法 |
Non-Patent Citations (1)
| Title |
|---|
| YIWEN YAO: "Design of 3D SAW Filters Based on the Spectral Element Method", IEEE ACCESS, 16 March 2023 (2023-03-16) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11626852B2 (en) | Hierarchical cascading in two-dimensional finite element method simulation of acoustic wave filter devices | |
| Yilmaz et al. | Time domain adaptive integral method for surface integral equations | |
| CN106712738B (zh) | 根据频率响应要求设计声微波滤波器的方法 | |
| Kowalczyk et al. | Formulation of locally reacting surfaces in FDTD/K-DWM modelling of acoustic spaces | |
| US7844421B2 (en) | Hybrid finite element method for traveling surface acoustic waves with thickness effect | |
| US8001515B2 (en) | Simultaneous optimization of analog design parameters using a cost function of responses | |
| CN113962093A (zh) | 一种基于谱元法的声表面波滤波器设计方法 | |
| CN112232001A (zh) | 一种集成电路超宽频谐振响应的自适应确定方法及系统 | |
| Belforte et al. | Frequency domain behavior of S‐parameters piecewise‐linear fitting in a digital‐wave framework | |
| CN112380737B (zh) | 一种基于表面阻抗边界的薄层电磁结构的时域分析方法 | |
| Elsherbini et al. | Towards a simple model for saw delayline using cad | |
| Knockaert et al. | Laguerre-based bandlimited reduced-order modeling | |
| Grigoroscuta-Strugaru et al. | A modified discontinuous Galerkin method for solving efficiently Helmholtz problems | |
| CN117723894B (zh) | 基于微弱行波信号特征提取的故障检测方法和装置 | |
| WO2023032764A1 (ja) | 多端子キャパシタの等価回路モデル作成方法、等価回路モデル作成プログラム、等価回路モデル作成プログラムを記憶した記憶媒体、シミュレーション方法およびシミュレーション装置 | |
| Li | Efficient numerical modeling of miniature RF devices based on acoustic wave technologies | |
| Lalléchère et al. | Electrical stochastic modeling of cell for bio-electromagnetic compatibility applications | |
| CN116106637B (zh) | 一种时域电磁散射分析的有限冲激偶序列表征方法 | |
| Janeliauskas et al. | Initial design and quick analysis of SAW ultra–wideband HFM transducers | |
| Tasinkevych et al. | On electrostatics of finite planar system of strips applied in surface acoustic wave interdigital transducers | |
| Mei et al. | Prediction of Circuit Response to an Electromagnetic Environment (ASC IC FY2020 Milestone 7179) | |
| Ventura et al. | Analysis of SAW transducer having aperiodic multi-electrode cells using a coupled FEM/BIE numerical model | |
| Pólik et al. | RF inductor development by using the FEM | |
| Halpern | Schwarz waveform relaxation algorithms | |
| Tasinkevych | Methods of IDT charge spatial spectrum evaluation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |