圧電解析（静的）
理论与物理
压电效应
老师，压电效应是把机械力转换成电压吗？
结构和电磁场是耦合的呢。
是的。在FEM中，需要同时求解位移$u$和电位$\phi$作为未知数，这就是压电耦合分析。
主要压电材料
材料 $d_{33}$ [pC/N] 用途 PZT (锆钛酸铅) 300〜600 执行器、传感器 BaTiO₃ 190 陶瓷电容器 PVDF -33 柔性传感器 AlN 5 MEMS谐振器 LiNbO₃ 6 SAW滤波器
总结
- 机械-电耦合 — 应力↔电场双向作用
- 在FEM中同时求解$u$和$\phi$ — 压电耦合分析
- PZT应用最广泛 — $d_{33} = 300$〜600 pC/N
Coffee Break 闲谈
压电效应的发现——皮埃尔·居里与保罗·居里在1880年的实验
压电效应（Piezoelectric Effect）于1880年由法国的皮埃尔·居里（居里夫人的丈夫）与其兄保罗·雅克·居里在石英晶体中发现。他们证明了在石英上施加力时表面会出现电荷的“正压电效应”，次年（1881年）加布里埃尔·李普曼从理论上预测了相反的“逆压电效应（施加电场会产生变形）”。压电效应的工业应用始于20世纪初的声呐（水下声波探测），第一次世界大战期间保罗·朗之万开发了使用压电石英的水下探测器。现代智能手机的摄像头防抖执行器、医用超声波诊断设备、喷墨打印机的喷头，都是居里兄弟这一发现的后代。
各项的物理意义
- 电场项 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：法拉第电磁感应定律。随时间变化的磁通密度会产生电动势。【日常例子】自行车的发电机（磨电灯）通过旋转磁铁使附近的线圈产生电压——这是磁场随时间变化会感应出电场这一定律的直接应用。IH电磁炉也基于相同原理，高频磁场的变化在锅底感应出涡流，通过焦耳热加热。
- 磁场项 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：安培-麦克斯韦定律。电流和位移电流会产生磁场。【日常例子】电线通电时周围会产生磁场——这就是安培定律。电磁铁根据此原理工作，通过线圈通电产生强磁场。智能手机的扬声器也是应用此原理：电流→磁场→振膜的力。在高频（GHz频段天线等）情况下，位移电流 $\partial D/\partial t$ 不可忽略，它描述了电磁波的辐射。
- 高斯定律 $\nabla \cdot \mathbf{D} = \rho_v$：表明电荷是电位移矢量的发散源。【日常例子】用垫板摩擦头发后，静电会使头发竖起——带电的垫板（电荷）向外辐射电力线，对轻的头发施加力。电容器（电容）的设计中，电极间的电场分布就是根据此定律计算的。ESD（静电放电）对策也基于高斯定律的电场分析。
- 磁通守恒 $\nabla \cdot \mathbf{B} = 0$：表示不存在磁单极子。【日常例子】将条形磁铁切成两半，也无法得到只有N极或只有S极的磁铁——N极和S极总是成对出现。这意味着磁力线描绘的是“没有起点和终点的闭合回路”。在数值分析中，为了满足此条件，会使用矢量势 $\mathbf{B} = \nabla \times \mathbf{A}$ 这种公式化方法，自动保证磁通守恒。
假设条件与适用范围
- 线性材料假设：磁导率·介电常数不依赖于磁场·电场强度（饱和区域需要非线性B-H曲线）
- 准静态近似（低频）：位移电流项可忽略（$\omega \varepsilon \ll \sigma$）。涡流分析中常用
- 2D假设（截面分析）：电流方向均匀、可忽略端部效应时有效
- 各向同性假设：各向异性材料（如硅钢板的轧制方向等）需要定义方向特性
- 不适用的情形：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
老师，压电效应是把机械力转换成电压吗？
结构和电磁场是耦合的呢。
是的。在FEM中，需要同时求解位移$u$和电位$\phi$作为未知数，这就是压电耦合分析。
| 材料 | $d_{33}$ [pC/N] | 用途 |
|---|---|---|
| PZT (锆钛酸铅) | 300〜600 | 执行器、传感器 |
| BaTiO₃ | 190 | 陶瓷电容器 |
| PVDF | -33 | 柔性传感器 |
| AlN | 5 | MEMS谐振器 |
| LiNbO₃ | 6 | SAW滤波器 |
- 机械-电耦合 — 应力↔电场双向作用
- 在FEM中同时求解$u$和$\phi$ — 压电耦合分析
- PZT应用最广泛 — $d_{33} = 300$〜600 pC/N
压电效应的发现——皮埃尔·居里与保罗·居里在1880年的实验
压电效应（Piezoelectric Effect）于1880年由法国的皮埃尔·居里（居里夫人的丈夫）与其兄保罗·雅克·居里在石英晶体中发现。他们证明了在石英上施加力时表面会出现电荷的“正压电效应”，次年（1881年）加布里埃尔·李普曼从理论上预测了相反的“逆压电效应（施加电场会产生变形）”。压电效应的工业应用始于20世纪初的声呐（水下声波探测），第一次世界大战期间保罗·朗之万开发了使用压电石英的水下探测器。现代智能手机的摄像头防抖执行器、医用超声波诊断设备、喷墨打印机的喷头，都是居里兄弟这一发现的后代。
各项的物理意义
- 电场项 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：法拉第电磁感应定律。随时间变化的磁通密度会产生电动势。【日常例子】自行车的发电机（磨电灯）通过旋转磁铁使附近的线圈产生电压——这是磁场随时间变化会感应出电场这一定律的直接应用。IH电磁炉也基于相同原理，高频磁场的变化在锅底感应出涡流，通过焦耳热加热。
- 磁场项 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：安培-麦克斯韦定律。电流和位移电流会产生磁场。【日常例子】电线通电时周围会产生磁场——这就是安培定律。电磁铁根据此原理工作，通过线圈通电产生强磁场。智能手机的扬声器也是应用此原理：电流→磁场→振膜的力。在高频（GHz频段天线等）情况下，位移电流 $\partial D/\partial t$ 不可忽略，它描述了电磁波的辐射。
- 高斯定律 $\nabla \cdot \mathbf{D} = \rho_v$：表明电荷是电位移矢量的发散源。【日常例子】用垫板摩擦头发后，静电会使头发竖起——带电的垫板（电荷）向外辐射电力线，对轻的头发施加力。电容器（电容）的设计中，电极间的电场分布就是根据此定律计算的。ESD（静电放电）对策也基于高斯定律的电场分析。
- 磁通守恒 $\nabla \cdot \mathbf{B} = 0$：表示不存在磁单极子。【日常例子】将条形磁铁切成两半，也无法得到只有N极或只有S极的磁铁——N极和S极总是成对出现。这意味着磁力线描绘的是“没有起点和终点的闭合回路”。在数值分析中，为了满足此条件，会使用矢量势 $\mathbf{B} = \nabla \times \mathbf{A}$ 这种公式化方法，自动保证磁通守恒。
假设条件与适用范围
- 线性材料假设：磁导率·介电常数不依赖于磁场·电场强度（饱和区域需要非线性B-H曲线）
- 准静态近似（低频）：位移电流项可忽略（$\omega \varepsilon \ll \sigma$）。涡流分析中常用
- 2D假设（截面分析）：电流方向均匀、可忽略端部效应时有效
- 各向同性假设：各向异性材料（如硅钢板的轧制方向等）需要定义方向特性
- 不适用的情形：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
数值解法与实现
压电FEM的公式化
离散化后:
$[K_{uu}]$: 机械刚度、$[K_{\phi\phi}]$: 介电刚度、$[K_{u\phi}]$: 压电耦合项。
结构和电学的自由度被整合到一个矩阵里了呢。
每个节点有位移DOF（$u_x, u_y, u_z$）和电位DOF（$\phi$），合计4个（3D）。Abaqus的压电单元C3D8E等支持此功能。
总结
- 耦合矩阵 — 同时求解机械与电学
- Abaqus C3D8E / COMSOL Piezoelectric — 商业实现
压电FEM的耦合设置——机械-电耦合方程的弱形式与材料张量的处理
压电材料的FEM分析是“机械-电耦合问题”，需要联立求解位移u（机械场）和电势phi（电场）。弱形式化将压电本构方程（应力张量＝弹性常数×应变－压电常数×电场）纳入虚功原理，用刚度矩阵K（机械）、介电矩阵（电）、耦合矩阵（压电）这三个块矩阵表示。实现上的陷阱是压电张量e的坐标变换——如果晶体轴方向与分析坐标系不匹配，耦合会变弱，导致电压输出被低估。COMSOL的Piezoelectric Devices Physics和ANSYS Mechanical的PIEZO模块会自动处理此变换，但用户自定义材料需要手动设置变换矩阵。
边单元（Nedelec单元）
专用于电磁场分析的单元。自动保证切向分量的连续性，排除伪模式。是3D高频分析的标准。
节点单元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM: 对应非线性材料·非均匀介质。BEM: 可自然处理无限域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿-拉弗森法处理B-H曲线的非线性。残差标准: $||R||/||R_0|| < 10^{-4}$ 是通用的。
频域分析
通过时间谐波假设归结为稳态问题。需要进行复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分可以使用更大的步长，但需注意精度。
频域与时域的使用区分
频域分析类似于“将收音机调谐到特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——可以再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实务
超声波换能器、压电执行器（喷墨头）、振动传感器、能量收集器的设计。
检查清单
超声波清洗机的压电换能器设计——空化分布与FEM优化
超声波清洗机利用压电换能器产生的超声波振动在液体中产生空化气泡，通过其溃灭冲击去除污垢。换能器的谐振频率设计（20〜40kHz）与清洗槽内的声压分布决定了清洗均匀性。通过FEM（压电耦合）分析换能器的振动模式与声场（亥姆霍兹方程），优化配置与形状以使声压均匀。在实际设计案例中，通过将4个圆盘型换能器非对称配置，使空化分布的均匀性改善了30%，将精密清洗产线的清洗不良率从0.8%降低到0.1%以下。FEM的设计优化直接关系到半导体晶圆清洗、医疗器械清洗的质量管理。
分析流程的比喻
电机的电磁场分析，其感觉接近于“给吉他调音”。调整琴弦的粗细（线圈匝数）和琴桥的位置（磁铁配置），以引出最美妙的音色（高效的扭矩特性）。改变一个参数，整体的平衡就会改变——所以参数化研究很重要。
初学者容易陷入的陷阱
“空气区域？为什么要用网格划分空气？”——这是几乎所有初次接触电磁场分析的人都会产生的疑问。答案是“因为磁力线也会扩散到铁芯之外”。如果将分析区域刚好设在铁芯边界，无处可去的磁通会“撞上”墙壁反射，产生实际中不可能出现的磁通集中。想象一下房间太小，球在墙上不断弹跳的状态。
边界条件的思考方式
远场边界条件虽然不起眼但超级重要。需要在数值上表达“从这里开始是无限广阔的空间”。如果设置错误，磁通就会像撞上“看不见的墙”一样被反射回来。
软件比较
工具
| 工具 | 特点 |
|---|---|
| COMSOL (Piezoelectric) | 最易用。GUI设置极化方向 |
| Abaqus | C3D8E等压电单元。研究标准 |
| ANSYS Mechanical | PIEZO模块。与ANSYS Electronics Suite集成 |
| OnScale (旧OnScale) | 云端压电FEM。擅长超声波器件 |
| OpenFOAM (piezoFoam) | 开源。需要编程知识 |
总结
- COMSOL — 快速原型设计
- Abaqus — 学术·研究标准
- OnScale — 超声波·MEMS专用
云端FEM的崛起——OnScale如何将压电超声波仿真时间从数小时缩短到数分钟
传统的压电FEM（如Abaqus）在本地工作站上运行，一个3D超声波换能器的瞬态分析可能需要数小时到一天。2016年成立的OnScale（现为Ansys OnScale）将FEM求解器移植到云端（AWS、GCP），利用近乎无限的CPU/GPU并行计算资源，将同样的分析缩短到数分钟。其秘诀是“降阶建模”与“参数化求解数据库”——用户改变几何参数时，系统并非从头计算，而是从预计算的响应曲面中插值出结果。一家助听器制造商利用OnScale优化其MEMS超声波距离传感器，在2周内完成了1,200种设计变体的评估，而传统方法需要6个月。云端FEM正在改变“仿真驱动设计”的速度与规模。
COMSOL的优点
- GUI直观，易于设置极化方向与电极边界条件
- 多物理场耦合（压电+声学+热）容易实现
- 丰富的材料库与后处理功能
Abaqus的优点
- 非线性材料模型（如铁电迟滞）支持强大
- 隐式/显式求解器选择，适合动态冲击分析
- 用户子程序（UMAT、UEL）允许高度自定义
OnScale的优点
- 云端求解，无需高性能本地硬件
- 专为超声波/MEMS器件优化，计算速度快
- 参数化扫描与DOE（实验设计）功能强大
开源现状
- OpenFOAM的piezoFoam：基础功能具备，但文档与GUI支持有限
- FEniCS、MFEM：研究用框架，需要较强的编程能力
- 商业软件在易用性与支持方面仍占优势
なった
詳しく
報告