在为现代通信系统设计压电谐振器时，高保真的多物理场仿真对于捕捉复杂的器件行为是不可或缺的。但在实际的射频前端中，通常是多个组件级联，或构建外部匹配网络。在这种情况下，进行完整的有限元（FEM）仿真将带来极大的计算开销。修正型 Butterworth–Van Dyke (mBVD) 等效电路模型正是为解决这一问题而生。通过将有限元数据凝练为紧凑的集总模型，您可以在不牺牲物理准确度的前提下，实现近乎实时的系统级仿真。
谐振器基础：物理机制、导纳、阻尼与机电耦合
在深入讨论等效电路形式之前，我们先从物理层面了解谐振器的基本行为。标准的体声波（BAW）或声表面波（SAW）谐振器利用压电效应，将声能束缚在微尺度腔体内，并使其在固有频率处谐振。为了从电学角度分析这一行为，我们考察端口导纳Y=I/V，其中：
- 实部，即电导G=real(Y)，代表能量转移
- 虚部，即电纳B=imag(Y)，追踪能量存储
这种机电相互作用会产生两个主要的谐振频率。在串联谐振频率f_s处，声学驻波发生相长干涉，器件内部机械阻抗降至最低，导致电导率B出现尖锐峰值。反谐振或并联谐振频率f_p略高于此。在此频率下，机械应变产生的压电电荷与电极上的静电板电荷相互抵消。器件汲取的净电流极小，使得总导纳幅值降至最低。
f_s与f_p之间的间隔决定了有效机电耦合系数k^2_{eff}的大小，它度量了器件的能量转换效率：k^2_{eff} \propto (f_p^2 – f_s^2)/f_s^2。频率间隔越大，k^2_{eff}越高，这通常意味着滤波器可以实现更高的带宽。有效机电耦合系数的准确定义取决于所遵循的标准或近似方法。对于 BAW 和 SAW 器件，其表达式可能略有不同；当系统满足弱耦合近似条件时，所使用的公式也可能有所简化。
在理想无损的谐振器中，f_s处的电导率峰值将趋于无穷大，而f_p处的导纳谷值将降至绝对零。然而，现实中压电器件会因为材料自身阻尼（机械的、介电的，以及其他机制，如锚点向基板的能量泄露）而产生能量损耗。这种总耗散决定了器件的品质因子，即Qf。在导纳曲线上，更高的材料阻尼会直接抑制并展平响应，使尖锐的谐振峰和谐振谷变得圆钝，降低f_s处的最大电导并抬高f_p处的最小导纳。对于射频滤波器设计人员而言，最小化这种阻尼以实现高 Q 至关重要，这正是确保低插损和高陡降的关键因素。
导纳幅值 |Y| 曲线示意图，其中标出了串联谐振和并联反谐振频率点。
用有限元方法表示压电谐振器
您可以在 COMSOL Multiphysics®中轻松构建一个经典的压电谐振器有限元模型，并充分利用该软件强大的多物理场能力。具体步骤如下：首先选用压电接口，耦合静电与固体力学接口。对于固体力学接口，选择压电材料材料模型，并通过其子节点（如机械阻尼和介电损耗）指定相关的损耗行为。
进行频域研究时，您可以在静电接口中添加终端条件来给器件施加激励。除了其他优点外，该特征可以自动推出终端导纳，并用于绘图和进一步的结果分析。结果分析则包括上述谐振频率和品质因子 Q 的提取。或者，您也可以选择进行特征频率分析以提取谐振频率和 Q。
从有限元到集总参数
集总参数电路模型的优势在于其简洁性：求解由少量电阻、电感和电容组成的网络，仅需不到一秒便可以得到所需的导纳曲线。本文所关注的机电等效框架，其渊源可追溯至 Stephen Butterworth 的滤波器设计以及 Walter Guyton Van Dyke 在 20 世纪初对石英晶体行为的等效工作。“修正型”现代版本增加了寄生元件，以满足高频 MEMS 谐振器的更为精细的建模要求。
现代射频 MEMS 谐振器需要快速电路模型与现代数值求解器之间的工程协同。电路模型的好坏与输入参数强相关。您无法仅凭猜测就得到复杂多层 SAW 叉指拓扑的等效电容，或者一个新型压电薄膜的动态支路电阻。通过使用 COMSOL Multiphysics®进行仿真，您可以还原器件的真实几何与材料复杂性，还可以提取填充 mBVD 电路所需的精确集总参数。总而言之，有限元法提供了完全符合物理真实的基准数据，而电路形式等效则实现了系统级的仿真速度。
注：本文讨论的经典 mBVD 形式主要针对单端口谐振器，即具有一个信号端和一个接地端的器件。如果您的设计归结为多端口设计或类似情况，则需要将该框架扩展为更复杂的形式。
同样，为清晰起见，我们将重点讨论单谐振模态。对于需要额外考虑的每一个谐振模态，您都需要添加一条并联的动态支路。
拆解 mBVD 电路模型
传统 mBVD 模型通过将初始串联路径之后的电学行为拆分为两条并联支路，来描述典型的压电谐振器行为：一条静态支路描述物理介电层的堆叠关系，另一条动态支路描述机电耦合的声学谐振。电路原理图如下所示。
单端口 mBVD 电路拓扑示意图，其中串联电阻R_s与由动态支路和静态支路构成的并联网络相连。动态支路由R_m、L_m和C_m串联组成；静态支路则由C_0和R_0二者并联组成。
下表逐一梳理了每个元件的物理含义、这些元件如何影响导纳响应，以及如何从有限元仿真中精准地提取对应参数值。
| 电路元件 | 物理含义 | 对导纳的影响 | 提取方法 |
|---|---|---|---|
| 串联电阻，R_s | 该元件代表金属电极、母线板以及通向谐振器有源区的互连导线的纯欧姆电阻。 | 它限制了串联谐振f_s处导纳幅值的最大峰值。 | R_s可直接根据模型中金属布线的电导率和几何形状计算，例如通过辅助的电流设置。通常无需在主压电仿真中显式包含它。 |
静态电容，C_0 |
C_0代表由重叠电极和中间压电材料形成的电学板电容。C_0与机械运动无关。 | 它建立了导纳虚部曲线的基线“底部”，并决定了串联谐振f_s与并联反谐振f_p之间的频率间隔。 | 在远低于f_l的任意频率点f_s处运行频域仿真，此时器件在机械上“静止”，动态支路相当于开路。计算导纳的虚部以分离出静态电容：C_0=B(f_l)/(2 \pi f_l) |
| 静态电阻，R_0 | R_0代表体介电材料损耗以及穿过压电层的电学泄露路径。 | 它在低频处为导纳实部曲线引入一条干净、平坦的基线电导。 | 假设您的模型具有非零介电损耗，在用于求C_0的同一低频区间（即远低于声学谐振频率），测量导纳的实部，直接分离出平坦的基线介电泄露R_0=1/G(f_l)。对于小介电损耗的情况，R_0应得到非常大的值。 |
| 动态电容，C_m | 该元件代表谐振器结构的机械弹性或柔度，与压电晶体的有效机电耦合系数k^2_{eff}密切相关。 | 它决定了谐振的整体强度。更高的耦合系数会增大C_m。进而，它拓宽了f_s与f_p之间的窗口，这是宽带射频滤波器的关键特性。 | C_m可利用提取的C_0值和谐振频率比，通过C_m=C_0 ( f_p/f_s)^2-1)直接求得。 |
| 动态电感，L_m | 该元件代表振动晶格或薄膜膜片的有效声学质量或惯性。 | 与C_m一起，它决定了串联谐振峰的确切频率位置。该计算遵循谐振频率的经典公式f_s=1/ (2 \pi \sqrt{L_m C_m}). | 一旦已知C_m并从仿真数据中识别出f_s，即可定义动态电感为L_m=1/((2\pi f_s)^2 C_m)。 |
| 动态电阻，R_m | 该电阻考虑了器件中的总力学和声学能耗。 | 它决定了谐振处导纳峰的高度和锐度。较低的R_m值意味着较低的损耗，从而获得更高的品质因子 Q。 | 在频域扫描中，f_s处的峰值导纳由动态支路主导。峰值电导率简化为R_m \approx 1/G(f_s)或者，您可以从特征频率研究中获取器件的机械品质因子Q_{fm}，并使用以下关系：R_m = \sqrt{L_m/C_m}/Q_{fm}. |
本文中，我们选择了静态支路中C_0和R_0并联的 mBVD 电路形式。行业内也存在将这些元件串联的惯例，如果采用这种形式，其参数提取方法和物理解释也会随之不同。这两种方法是互易的，因为只要能够正确、谨慎地提取集总参数，输出结果将完全一致。
参数提取
既然我们已经掌握了理论知识和提取逻辑，接下来我们就来提取整套集总参数，并在 COMSOL Multiphysics®中构建 mBVD 电路模型。
请注意，本篇博客文章是该模型的主要参考资料。如果您有兴趣查找利用类似模型和绘制导纳曲线的分步操作说明，请查看我们为MEMS 模块提供的配套示例模型。
在这里，我们选择了一个三维层状单端口薄膜体声波谐振器（FBAR），如下所示。我们假设R_s是已知的，为了便于后续比较，仅将其手动添加到器件导纳中作为贡献。
为了从频域研究中准确提取谐振频率和峰值电导率，您可能需要细化频率步长，并使用图形标记来绘制全局导纳图。如果考虑R_0，请不要忘记指定介电损耗。最后，为了提取R_0和C_0，您有必要新建一个研究，以求解一个比f_s小几个数量级的频点。
COMSOL Multiphysics®用户界面中显示了三维 FBAR 模型，其导纳通过频域研究计算得出。
一旦所有参数计算完毕，您可以将这些值直接代入同一模型文件中内置的电路接口，以验证提取结果。随后，将新创建的集总 mBVD 电路接入电压源，并将其阻抗频响曲线与原始有限元数据进行叠加对比。顺应当前便捷的 AI 趋势，您不妨利用 COMSOL Multiphysics®中内嵌的智能助手窗口，让它生成添加电路及完成各项设置的 COMSOL API 代码，此处的演示就采用了这种方式。
提取的 mBVD 参数用于电路接口，绘图描述了完整有限元模型与等效电路之间的良好相关性。
两条曲线结果在工作带宽附近非常吻合。从实际角度来看，两者存下的微小差异可以解释为有限元模型捕捉了系统中存在的高阶模式的贡献。一旦您的 mBVD 电路模型与基于物理场的模型匹配，就不再受限于器件的物理几何形状，并可以在 COMSOL Multiphysics®中直接扩展电路网络，添加外部调谐电容、电感或传输线，或级联多个相同模块以形成复杂的梯形滤波器。综上，您最终获得了一个经过验证的、运行速度极快的集总模型，它实现了微观多物理场与宏观系统设计的无缝衔接，使您能够在几秒钟内完成大规模参数扫描。
下一步
最后，分享本篇博客文章中展示的模型以及相关示例：
- 查看本篇博客文章中讨论的模型：从三维 BAW 谐振器提取 mBVD 电路参数
- 了解如何从特征频率角度处理同一任务，请参阅：SAW 基本单元的频域分析
- 如果您需要完美匹配和/或补偿缺失数据（例如阻尼比）的方法，请参阅此示例，其中利用了参数估计特征：具有等效电路的薄膜体声波谐振器
- 查看使用MEMS 模块的教学案例：案例下载
评论 (0)