3. COMSOL压电接口详解

在COMSOL Multiphysics中建立压电仿真，核心是正确理解并配置“压电效应”多物理场接口。该接口本质上把固体力学与静电学两个物理场耦合起来，通过本构方程描述应力、应变、电场与电位移之间的相互转换。本节围绕方程形式、材料参数输入、坐标系设定、几何非线性开关、参考温度与热膨胀耦合五个关键配置项展开，帮助你把接口设置与真实器件的物理行为对应起来。

一、接口的方程形式：应力-电荷与应变-电荷

压电接口提供两种本构形式。默认的“应力-电荷”形式以应力T和电场E为自变量，本构关系写作：

S = sE : T + dT · E
D = d : T + εT · E

其中sE是恒定电场下的柔度矩阵，d是压电耦合矩阵，εT是恒定应力下的介电矩阵。这种形式适合已知外加载荷、需要求解位移和电势的常规驱动与传感分析。另一种“应变-电荷”形式以应变S和电场E为自变量：

T = cE : S − eT · E
D = e : S + εS · E

其中cE为恒定电场下的刚度矩阵，e为压电应力常数，εS为恒定应变下的介电矩阵。当材料数据手册给出的是刚度与压电应力常数，或需要做模态与频域分析时，选用应变-电荷形式更自然。两种形式在数学上等价，但输入参数不同，切换后必须重新核对材料数据，否则会因单位或常数含义不一致导致结果偏差。

二、材料参数输入方式

压电材料通常具有6mm或4mm对称性，其材料参数以矩阵形式给出。在COMSOL中，你可以选择“各向异性”并直接填写完整矩阵，也可以使用“压电材料”节点下的预定义材料。以PZT-5H为例，典型输入包括：柔度矩阵sE的12个独立分量、压电常数矩阵d的3个独立分量（d31、d33、d15）、以及相对介电常数εr的2个独立分量（ε11、ε33）。输入时注意单位：sE常用10⁻¹² m²/N，d常用10⁻¹² C/N或pC/N，εr为无量纲。若从文献复制数据，务必确认其参考坐标系与你的极化方向一致。

三、坐标系设定：极化方向与材料主方向

压电材料的各向异性由极化方向决定。COMSOL默认以全局坐标的Z轴作为极化方向，即材料主方向3对应Z轴。若你的器件极化方向沿X或Y，必须新建一个旋转坐标系或使用“材料坐标系”节点，把材料主方向映射到实际极化方向。常见做法是：在“压电材料”节点下选择“坐标系”，指定一个绕某轴旋转的坐标系，旋转角由极化方向与全局Z轴的夹角确定。例如极化沿X轴时，绕Y轴旋转−90°。坐标系设置错误是压电仿真中最隐蔽的错误之一，表现为耦合方向反号或耦合强度异常，务必在建模初期确认。

四、几何非线性开关

压电接口的“几何非线性”选项控制是否在变形后的构型上计算应变与电场。对于小变形、线性压电分析（如常规谐振器、小信号传感器），保持关闭即可，计算量小且收敛稳定。当器件发生大位移、大转动或屈曲，例如压电弯曲驱动器、超声换能器的大振幅振动、压电能量采集器的大挠度梁，必须打开几何非线性。打开后，应变采用Green-Lagrange应变，电场与电位移在变形构型中定义，本构关系中的应力与电位移也随之更新。注意：几何非线性会显著增加计算成本，且需要更细的载荷步和更严格的收敛控制。建议先用线性分析获得初步结果，再打开非线性做验证。

五、参考温度与热膨胀耦合

压电材料在温度变化时会产生热应变，同时压电常数和介电常数也随温度漂移。COMSOL压电接口允许设置“参考温度”Tref，并可选启用热膨胀耦合。当启用时，总应变分解为机械应变、压电应变和热应变：

S_total = S_mech + dT·E + α(T − Tref)

其中α为热膨胀系数。参考温度应设为器件无热应力时的温度，通常为制造温度或室温。若忽略热膨胀，则只考虑压电耦合，适用于温度变化小的场景。若做温度循环或高温下的压电响应分析，必须同时启用热膨胀并输入α(T)的温度依赖关系。此外，压电常数d和介电常数ε也常随温度变化，可在材料节点中定义为温度的函数，例如d33(T) = d33_0·(1 + k·(T − Tref))。这样在频域或瞬态分析中，温度场变化会自动更新压电耦合强度。

六、配置检查清单

完成上述设置后，建议按以下顺序核对：本构形式是否与材料数据匹配；sE、d、εr的单位与数值是否与极化方向一致；材料坐标系是否把主方向3对准实际极化方向；几何非线性是否与预期变形量级匹配；参考温度是否等于无应力温度；热膨胀系数是否已输入且与温度相关。若做频域分析，还需确认介电损耗和机械损耗是否已加入，否则谐振峰品质因数会偏高。把这些配置项逐一确认后，压电接口的物理行为才能与真实器件对应，后续的网格、求解器和后处理才有可靠基础。