4. 压电结构静力学仿真
悬臂梁式压电驱动器是最典型的压电结构耦合问题之一：一端固定、另一端自由，压电层在电压激励下产生面内应变，通过层间约束把应变转化为弯曲曲率，从而在自由端形成可观的横向位移。本节以“PZT-5H 压电层 + 铜弹性层”的双层悬臂梁为教学对象，完整走通几何建模、材料与物理场设置、电压激励、位移场求解、应力提取、网格收敛性验证与后处理这条链路。所有数值均为教学示例参数，用于说明量级与判断方法，不代表实测结果。
一、几何与材料：先定参照系，再定叠层
建模前必须先固定坐标系约定，否则位移方向、应力正负号会在后处理中互相矛盾。本节统一约定：梁沿 x 方向延伸，长度 L = 20 mm；厚度沿 y 方向，宽度沿 z 方向，w = 5 mm；固定端位于 x = 0，自由端位于 x = L。压电层厚度 t_p = 0.2 mm，铜层厚度 t_s = 0.1 mm，压电层位于铜层上方，二者共节点贴合。压电层极化方向沿 y 轴（厚度方向），即 d31 模式工作。
材料参数示例：PZT-5H 的相对介电常数 εr33 ≈ 3400，压电常数 d31 ≈ -274 pC/N，弹性柔度 s11 ≈ 16.4×10⁻¹² m²/N；铜的杨氏模量 E ≈ 110 GPa，泊松比 0.34。COMSOL 中可直接调用内置压电材料库，但需核对极化方向与坐标轴是否一致——这是新手最常见的错误来源：极化方向设错，位移会整体反向或量级异常。
二、物理场与边界：压电耦合的“三件套”
在 COMSOL 中新建“压电效应”多物理场接口，它会自动耦合“固体力学”与“静电”两个物理场。关键设置如下：
- 固体力学：固定约束施加在
x = 0端面；其余表面自由。压电层与铜层之间使用“连续性”或共享边界，保证位移与电势连续。 - 静电：压电层上表面设为“终端”，施加电压
V = 100 V；下表面（与铜层界面）设为“接地”。铜层本身作为等势体，可整体设为接地或终端。 - 压电材料：在压电层域上指定极化方向为
+y，并确认d31的符号与坐标系一致。
这里有一个容易忽略的细节：压电层下表面与铜层贴合，若铜层设为接地，则压电层下表面电势被强制为 0，上表面为 100 V，厚度方向电场 E_y ≈ V/t_p = 100/0.2 mm = 500 V/mm。电场方向与极化方向一致时，d31 为负，压电层在 x 方向收缩；由于铜层不收缩，双层结构产生弯曲，自由端向上或向下位移取决于叠层顺序与极化方向。建议在第一次求解后先检查位移方向是否符合预期，再继续做参数化扫描。
三、网格与收敛性：不要只看“能不能算”
压电层厚度只有 0.2 mm，而梁长 20 mm，长厚比达到 100:1。若用自由四面体网格直接剖分，压电层厚度方向至少需要 2～3 层单元才能捕捉弯曲应变梯度，否则应力结果会严重失真。推荐做法：在厚度方向使用“扫掠”网格，压电层和铜层各划分 3 层；在长度方向用“分布”控制，初始单元数设为 40，然后做网格收敛性验证。
收敛性验证的判据不是“网格越密越好”，而是关注两个量：自由端位移 δ_tip 和压电层上表面最大 σxx。以单元数 N 为变量，依次取 20、40、80、160，记录 δ_tip 与 σxx_max。若从 80 到 160 时 δ_tip 变化小于 1%、σxx_max 变化小于 3%，即可认为网格收敛。教学示例中，N = 40 时 δ_tip ≈ 0.42 mm，N = 80 时 δ_tip ≈ 0.45 mm，N = 160 时 δ_tip ≈ 0.455 mm，说明 80 单元已接近收敛。注意：这些数值是示例量级，实际结果受材料参数、层间连接方式、边界条件细节影响。
δ_tip，再用“一维绘图组”画 δ_tip 随单元数的变化曲线。若曲线趋于水平，说明收敛；若仍单调上升或下降，需继续加密或检查是否存在应力奇异。
四、应力分布提取：看什么、在哪看
压电悬臂梁的应力分布有两个关键区域：固定端附近的压电层上表面（拉应力最大）和压电层与铜层界面（剪应力集中）。后处理时，建议依次查看：
- 三维表面图：选择压电层上表面，绘制
σxx分布，观察从固定端到自由端的衰减趋势。固定端σxx最大，自由端趋近于零。 - 厚度方向线图：在
x = 2 mm处沿y方向取一条线，绘制σxx随y的变化，可清晰看到压电层受压/受拉与铜层的应力梯度。 - 界面剪应力：在压电层与铜层界面绘制
σxy，检查是否存在异常集中。若界面网格过粗，剪应力会出现锯齿状振荡，这是网格不足的典型信号。
位移场方面，除了自由端 δ_tip，还应检查梁的挠曲线形状是否光滑。若挠曲线出现局部折点，通常意味着层间连接设置不当或网格在界面处不连续。
五、结果后处理与工程判断
完成求解后，后处理的目标不是“看图好看”，而是回答三个工程问题：位移是否达到设计行程？应力是否低于材料许用值？网格是否足够可信？
| 检查项 | 示例结果 | 判断依据 |
|---|---|---|
| 自由端位移 δ_tip | ≈ 0.45 mm @ 100 V | 与设计行程对比，偏差 < 10% 可接受 |
| 压电层最大 σxx | ≈ 28 MPa | PZT-5H 抗拉强度约 50 MPa，安全系数 > 1.5 |
| 网格收敛性 | 80→160 单元 δ_tip 变化 < 1% | 位移收敛，应力变化 < 3% |
| 界面剪应力 | 无锯齿振荡 | 界面网格连续，无应力奇异 |
y 轴一致，以及电压施加面是否与极化方向匹配。若应力结果在固定端出现异常高值，可能是固定约束施加在了单个节点而非端面，导致应力奇异；应改为端面固定约束。
最后，建议把上述检查项整理成一个可复用的后处理模板：在 COMSOL 中建立“结果 → 派生值 → 全局计算”节点，分别定义 δ_tip、σxx_max、σxy_max，并保存为方法文件。这样每次修改电压或材料参数后，只需重新求解即可一键输出关键指标，大幅提升迭代效率。