2. 压电材料参数体系
在 COMSOL 中做压电耦合仿真，第一步不是画几何，而是把材料参数“喂对”。压电材料本构关系同时涉及力学量与电学量，参数之间并非彼此独立，而是通过晶体对称性和热力学关系相互约束。理解这套参数体系，才能判断材料库数据是否可用、缺失参数如何补齐、仿真结果是否可信。
一、四类核心参数及其物理含义
压电本构方程常用应力-电荷型（d 型）或应变-电荷型（e 型）表达。以 d 型为例：
S = sE · T + dᵗ · E
D = d · T + εT · E其中 S 为应变，T 为应力，E 为电场，D 为电位移。弹性柔度矩阵 sE 描述恒定电场下应力产生应变的能力，是 6×6 对称矩阵，因晶体对称性可大幅约简：PZT 陶瓷（6mm 点群）只有 5 个独立分量（s11、s12、s13、s33、s44），且 s66=2(s11−s12)。压电应变常数 d 是 3×6 矩阵，表示单位电场产生的应变（逆压电效应）或单位应力产生的电位移（正压电效应），常用分量 d33、d31、d15。介电常数矩阵 εT 为 3×3 对角阵，恒应力下只有 ε11 与 ε33 两个独立分量，工程中常以相对介电常数 εr=ε/ε0 给出。
压电应力常数 e 与 d 描述同一物理效应，只是自变量不同，二者通过弹性刚度矩阵 cE 换算：e = d · cE。COMSOL 的“压电材料”节点允许选择应变-电荷型或应力-电荷型，若材料数据以 e 型给出，可直接输入；若只有 d 型，软件会依据 cE 自动转换，但前提是弹性矩阵完整且单位一致。
二、耦合系数与机电转换效率
耦合系数 k 是衡量机电能量互换能力的无量纲指标，定义为转换能量与输入能量之比的平方根。对 33 方向：
k33² = d33² / (s33E · ε33T)k 越大，同一电场下产生的应变或同一应力下产生的电荷越多。典型软 PZT（如 PZT-5H）k33≈0.75、k31≈0.39；硬 PZT（如 PZT-8）k33≈0.64、k31≈0.30。需注意 k 并非效率：谐振状态下机电转换效率 η 还受机械品质因数 Qm 和介电损耗 tanδ 影响，η 可近似写为
η ≈ k² / (1 + (π/2)·(tanδ/k²)·(1/Qm))因此高 k 材料若损耗大，实际效率未必高。仿真中若只关心静态位移或电荷，k 足够；若做谐振换能器，必须同时给出 Qm 与 tanδ，否则导纳曲线峰值和带宽会失真。
三、参数获取途径与 COMSOL 落地
参数来源主要有四类。第一，厂商数据表：如 PZT-5H 给出 d33=593 pC/N、d31=−274 pC/N、ε33T/ε0=3400、s33E=2.07×10⁻¹¹ m²/N，但往往缺少完整 sE 矩阵和 e 矩阵，需按对称性补全。第二，COMSOL 内置材料库：Piezoelectric 目录下含 PZT-4、PZT-5H、BaTiO₃、石英等，已按 6mm 或 32 点群填好 d、sE、εT，可直接调用，但需核对极化方向与坐标轴是否一致。第三，文献与标准：IEEE 176 标准给出测量与换算规范，谐振-反谐振法可测 k 与 sE，准静态法可测 d。第四，反演与拟合：用阻抗分析仪测导纳曲线，再通过优化反推完整参数集，适合自研材料。
在 COMSOL 中操作时，建议按以下顺序自检：先确认材料节点选择的是“应力-电荷”还是“应变-电荷”形式；再检查弹性矩阵是否为正定（否则求解器报错或出现非物理模态）；然后核对 d 与 e 是否同时被误填（二者只能选其一，否则重复耦合）；最后用一个小立方体做静态验证——施加 1 V 电压，读取自由应变，与 d33·E 手算值对比，偏差应在 1% 以内。若做谐振分析，还需在材料节点补充 Qm 与介电损耗，否则导纳峰过于尖锐，与实际换能器不符。
掌握这套参数体系后，后续的压电耦合建模、极化方向设置、谐振与瞬态分析才有可靠基础。参数不是孤立的数字，而是一张相互约束的网络，理解换算链与物理边界，比记住单个数值更重要。