1. 压电耦合仿真入门
压电结构耦合仿真的核心，是把“电”与“力”之间可逆的双向转换关系放进同一个数值模型中求解。理解这一点，是后续所有建模操作的前提。本章从物理本质出发，建立正逆压电效应的数学描述，梳理 COMSOL 多物理场耦合界面的组织逻辑，并给出典型应用场景的对照，帮助你在动手建模前先建立正确的物理图像。
一、压电效应的物理本质
压电效应来源于某些晶体材料内部结构的非中心对称性。以石英、锆钛酸铅（PZT）、铌酸锂等为代表的压电材料，其晶胞在无外场时正负电荷中心重合，宏观上不显电性；一旦受到机械应力，晶胞发生畸变，正负电荷中心分离，形成电偶极矩，宏观上表现为材料两端出现束缚电荷，这就是正压电效应——机械能转换为电能。反之，施加外电场会使带电离子发生位移，晶格产生应变，宏观上表现为材料伸缩，这就是逆压电效应——电能转换为机械能。
需要强调三点：其一，压电效应是线性可逆的，正逆效应由同一组压电常数描述，只是作用方向相反；其二，压电材料同时具有介电性和弹性，因此其本构关系是电学量与力学量的耦合；其三，压电效应只在非中心对称晶体中出现，且存在居里温度，超过该温度材料退极化，压电性消失。这三条决定了仿真中必须同时考虑电场、结构场以及二者的耦合项。
二、正逆压电效应的数学描述
在 COMSOL 中，压电材料采用应力-电荷型本构方程，其张量形式为：
S = sE : T + d^T · E
D = d : T + εT · E其中 S 为应变，T 为应力，E 为电场强度，D 为电位移；sE 为恒电场下的弹性柔顺系数，εT 为恒应力下的介电常数，d 为压电应变常数矩阵。第一式描述逆压电效应（电场产生应变），第二式描述正压电效应（应力产生电位移）。d 矩阵的转置出现在第一式中，正是“可逆”在数学上的体现。
以常用的 PZT-5H 为例，其压电常数 d33 ≈ 593 pC/N，d31 ≈ −274 pC/N，相对介电常数 εr33 ≈ 3400。d33 为正表示沿极化方向施加电场时材料沿同方向伸长；d31 为负表示垂直极化方向收缩。这些符号约定必须在建模时与坐标系定义保持一致，否则会得到方向相反的位移结果。
三、COMSOL 多物理场耦合界面总览
COMSOL 处理压电耦合的标准路径是“固体力学 + 静电 + 压电效应多物理场节点”。在模型开发器中，先添加固体力学接口（提供位移场 u），再添加静电接口（提供电势场 V），然后在“多物理场”下添加压电效应耦合节点，将两个物理场绑定。该节点会自动把本构方程中的耦合项写入方程组，无需手动输入 d 矩阵。
材料属性方面，压电材料需在“压电材料”节点下填写弹性矩阵、耦合矩阵和介电矩阵，可直接调用内置的 PZT-5H、PZT-4、石英等库材料。边界条件上，结构场常用固定约束、自由、指定位移；电场常用接地、终端（施加电压或电荷）、零电荷。求解器通常选择稳态或频域研究，频域下可进一步做阻抗、导纳和共振分析。
四、典型应用场景
压电结构耦合仿真的应用大致可分为四类，其关注量与建模侧重点差异明显：
| 应用类型 | 典型对象 | 关注输出 | 建模要点 |
|---|---|---|---|
| 传感器 | 加速度计、水听器 | 输出电压 / 电荷 | 正压电为主，关注灵敏度与频响 |
| 执行器 | 压电微位移台、喷墨头 | 位移 / 出力 | 逆压电为主，关注行程与迟滞 |
| 换能器 | 超声探头、清洗振子 | 谐振频率、导纳 | 频域求解，关注机电耦合系数 |
| 能量收集 | 振动俘能梁 | 输出功率 | 正压电 + 电路耦合，关注阻抗匹配 |
以压电微位移台为例，一个典型教学示例可设为：PZT-5H 长方体 20 mm × 5 mm × 1 mm，极化方向沿厚度（z 向），底面接地，顶面施加 100 V 电压，四周自由。此时逆压电效应使厚度方向伸长，理论自由应变约为 d33·E = 593e-12 × (100/1e-3) ≈ 5.93e-5，对应厚度变化约 59 nm。仿真中可在“终端”节点施加电压，稳态求解后读取顶面平均 z 向位移，与理论值对比即可验证模型是否正确。若结果偏差过大，应优先检查极化方向、单位制和材料参数是否与坐标系一致。