CAE如何进行解析？

压电材料需要联合求解力学方程和静电场方程。ANSYS、COMSOL、ABAQUS可以输入压电材料的本构关系（电气和机械耦合矩阵），进行联合FEA解析。用于超声波探头的振动模式分析和MEMS传感器灵敏度分析等。

超声波探伤试验中使用的超声波探针设计也可以用CAE进行吗？

可以。通过FEM计算压电振荡子的共振频率、位移灵敏度、阻抗特性，确定最优形状是当代的设计标准。对于排列成阵列的相控阵探伤器，需要考虑各素子间的干涉进行联合解析，没有仿真无法进行设计。

压电素子的非线性行为也能进行仿真吗？

可以。PZT在强电场和应力下会显示极化的滞后（非线性行为）。对此进行建模需要非线性压电材料本构关系，可以通过CAMSYS或MATEC专用工具，或者通过ABAQUS的UMAT来实现。在执行器精密定位控制系统设计中，非线性模型变得重要。

压电效应 — CAE术语解释

Q: 压电效应是指施加压力就会产生电力吗？

正确。对晶体施加力时会产生电荷（正压电效应），反之施加电压时会产生变形（逆压电效应），这是一种双向转换现象。石英和PZT（钛酸铅锆酸铅）是代表性的压电材料，用于超声波传感器、喷墨喷嘴、振动发电等领域。

分类: 术语表 | 2026-01-15

CAE visualization for piezoelectric effect - technical simulation diagram

压电效应

机械-电气转换的物理现象

🧑‍🎓

压电效应是指施加压力就会产生电力吗？

压电效应的理论基础

基本概念和支配方程

🧑‍🎓

压电效应，教科书上写的是"施加力产生电压的现象"，但具体是用什么材料，产生多少电压呢？

🎓

好问题。代表性的压电陶瓷PZT（钛酸铅锆酸铅）例如，施加10 MPa的压缩应力时，沿材料分极方向会产生约200 V/mm的电场。这是因为材料常数压电常数

$$ d_{33} $$

约为300 pC/N。具体电压与材料厚度成正比。厚度为1 mm的PZT素子施加10 MPa时，计算得约200 V的电压。

🧑‍🎓

反过来，施加电压产生变形的"逆压电效应"，公式就是简单的倒数关系吗？

🎓

本质上是对称的，但定式化略有不同。压电材料的基本线性本构方程按IEEE规范表示如下。

正压电效应（传感器）:

$$ D_i = d_{ijk} T_{jk} + \varepsilon_{ik}^T E_k $$

逆压电效应（执行器）:

$$ S_{ij} = s_{ijkl}^E T_{kl} + d_{kij} E_k $$

其中，

$$ D $$

是电感应强度，

$$ T $$

是应力，

$$ S $$

是应变，

$$ E $$

是电场。关键点是相同的压电常数

$$ d $$

联系两个方程。例如PZT的

$$ d_{33} $$

为300 pC/N，意味着"1 N的力产生300 pC的电荷"，同时也意味着"1 V/m的电场产生300 pm的应变"。

🧑‍🎓

需要区分"正压电效应"和"逆压电效应"进行仿真吗？还是用一个模型就能计算两个效应？

🎓

CAE通常将结合的本构方程作为一个"多物理场联合解析"同时求解。也就是说，一个模型可以计算由机械负荷产生的发电（传感）和由电压驱动的变形/振动（驱动）两种效应。只需改变边界条件和加载条件，就能评估两种效应。

压电效应的数值计算方法

FEM离散化和求解器设置

🧑‍🎓

用有限元法求解压电效应时，位移和电位都作为未知数。位移和电位的形状函数可以用相同的吗？

🎓

这是一个陷阱。很多实现中，位移使用二次要素（例如8节点六面体要素），而电位使用一次要素的"混合要素"被推荐。原因是电场

$$ E $$

是作为电位的梯度计算的。如果位移和电位都用相同的高阶形状函数，电场的精度会不必要地提高，引起数值"锁定"，导致刚度被高估。Abaqus手册中提及了C3D8E（一次电位）和C3D20E（二次电位）的区分使用。

🧑‍🎓

总体刚度矩阵是机械部分和电气部分耦合的，求解器需要特别的设置吗？

🎓

需要。耦合矩阵可能变成非对称的。但是，如果压电本构方程在能量上是保守系，矩阵就是对称的。许多商业软件（Ansys、COMSOL）默认使用对称求解器，但当考虑非线性或损失时，需要切换到非对称求解器。另外，机械自由度（位移）和电气自由度（电位）的量纲差异很大（例如位移为μm，电位为V），容易导致数值条件数恶化。为了防止这种情况，应启用缩放选项或谨慎选择单位制。

🧑‍🎓

求共振频率等频率响应时，特征值解析怎么做？用普通结构分析的模态分析手法可以吗？

🎓

需要特殊的"压电耦合模态分析"手法。支配方程使用质量矩阵

$$ \mathbf{M} $$

、刚度矩阵

$$ \mathbf{K}_{uu} $$

、介电矩阵

$$ \mathbf{K}_{\phi\phi} $$

和耦合矩阵

$$ \mathbf{K}_{u\phi} $$

，可写成：

$$ \begin{bmatrix} \mathbf{M} & \mathbf{0} \\ \mathbf{0} & \mathbf{0} \end{bmatrix} \begin{Bmatrix} \ddot{\mathbf{u}} \\ \ddot{\boldsymbol{\phi}} \end{Bmatrix} + \begin{bmatrix} \mathbf{K}_{uu} & \mathbf{K}_{u\phi} \\ \mathbf{K}_{u\phi}^T & -\mathbf{K}_{\phi\phi} \end{bmatrix} \begin{Bmatrix} \mathbf{u} \\ \boldsymbol{\phi} \end{Bmatrix} = \begin{Bmatrix} \mathbf{F} \\ \mathbf{Q} \end{Bmatrix} $$

求解这个特征值问题，会得到两组不同的共振频率，分别对应"短路条件"（电极间电位为零）和"开路条件"（电荷为零）。在Ansys中，执行`ANTYPE,MODAL`并指定`PIEZO`选项。

压电效应的实务应用

工作流程和检查清单

🧑‍🎓

请告诉我进行压电执行器位移仿真的初始建模步骤。该从哪里开始？

🎓

首先是输入材料常数。最低限度需要弹性矩阵

$$ \mathbf{s}^E $$

或

$$ \mathbf{c}^E $$

、压电常数矩阵

$$ \mathbf{d} $$

或

$$ \mathbf{e} $$

、介电率矩阵

$$ \boldsymbol{\varepsilon}^T $$

或

$$ \boldsymbol{\varepsilon}^S $$

这三组数据。从制造商目录（例如富士钛业的PZT-5H）或标准规范（IEEE Std 176-1987）获取数据。常数注意"定应变"还是"定电场"，以及是"d常数"还是"e常数"的表示，需要转换成软件要求的格式。

🧑‍🎓

边界条件要注意什么？像结构分析一样指定固定面和施加电压的面就行？

🎓

除此之外，电气的"接地"条件很重要。作为电极的面应用"电位"边界条件。通常，一个电极设为0 V（接地），另一个电极设为施加电压（例如100 V）。电极面被视为导体，面上电位为常数的"等电位条件"在大多数软件中自动应用。结构固定条件应反映实际的支持条件。例如，双层梁执行器的悬臂梁模型，固定端所有位移被完全约束。

🧑‍🎓

怎样验证结果才能令人信服？当位移量与目录值不符时，首先怀疑什么？

🎓

首先怀疑"材料常数输入错误"和"分极方向定义错误"。压电常数由三个轴向（1、2、3）和剪切（4、5、6）的组合定义。检查模型的材料坐标系（通常是单元坐标系）的3个方向是否与实际分极方向一致。其次检查边界条件，特别是电极面的指定是否正确。对于简单的平板执行器，存在理论解。施加电场为

$$ E_3 = V/t $$

，产生的应变为

$$ S_3 = d_{33} E_3 $$

，自由端位移为

$$ \delta = S_3 \cdot L $$

。先检查这个简单计算是否与FEM结果一致，从这里开始最好。

压电效应的软件比较

Ansys/Abaqus/COMSOL的特点

🧑‍🎓

用Ansys、Abaqus、COMSOL进行压电分析时，基本方法有本质区别吗？

🎓

最大的区别是耦合的处理方式。COMSOL Multiphysics以"多物理场"为核心理念，通过GUI将物理接口（结构力学和静电场）联接来定义压电效应。灵活性强，但用户有责任正确设置耦合。相反，Ansys Mechanical（APDL或Workbench）和Abaqus/Standard预定义了"压电要素"，只需输入材料常数就能使耦合方程内置。更像黑箱，但设置简单，使用优化的求解器。

🧑‍🎓

材料库的充实度如何？PZT等标准材料的数据从一开始就有吗？

🎓

COMSOL的材料库比较充实，包含几个通用PZT材料（例如PZT-5A、PZT-5H）的常数。Ansys和Abaqus的标准材料库中虽然通用金属和塑料很多，但压电材料几乎没有。用户通常需要从Morgan Electro Ceramics或PI Ceramics等制造商目录、或NIST数据库中查找常数自己输入。这是因为材料常数很大程度上取决于成分和制造工艺。

🧑‍🎓

频率响应和共振分析的功能有什么区别吗？

🎓

这里各软件的特点差异最大。Abaqus可以用`*FREQUENCY`步直接计算直接法和模态法频率响应，用`*SELECT EIGENFREQUENCIES`提取特定模式。Ansys APDL用`ANTYPE,HARMIC`和`HROPT,MSUP`可进行模态叠加法频率响应分析，能在共振周频附近追踪详细行为。COMSOL是选择"频率域"研究步，通过参数扫描改变频率，更直观。任何一个软件的关键都是能否正确设置"短路"和"开路"边界条件。

压电效应的故障排除

常见错误和对策

🧑‍🎓

运行分析时出现"刚度矩阵为奇异"或"对角项为负"的错误。可能是什么原因？

🎓

最多见的原因有两个。一是"电气浮动节点"。如果压电体中所有节点的电位自由度都没有被某些边界条件（电极指定或接地）约束，电位会有无穷大的解，刚度矩阵不会成为正定。检查模型中所有压电要素是否至少有一个节点的电位被固定。二是"材料常数输入错误"，特别是介电率矩阵

$$ \boldsymbol{\varepsilon} $$

的对角项是否以正值输入。负值或零在物理上不现实，会数值不稳定。

🧑‍🎓

共振频率分析中，计算频率远大于（或小于）目录值。即使细化网格也没有变化。

🎓

如果网格无关，问题在材料常数或边界条件。首先确认使用的弹性常数是"定电场下"

$$ \mathbf{c}^E $$

，不要误用"定电通密度下"

$$ \mathbf{c}^D $$

。错用后刚度被高估，共振频率会偏高。其次检查特征值分析的边界条件。如果要求"短路共振频率"，两个电极都接地为0V；如果要求"开路共振频率"，一个电极的电荷（或电流）设为0。设置错误会计算出完全不同的频率。在Abaqus中，用`*ELECTRODE`关键字定义电极，用`CURRENT=0`表示开路条件。

🧑‍🎓

计算执行器位移时，结果只有理论值的一半。可能是什么原因？

🎓

首先怀疑机械约束的影响。理论公式

$$ \delta = d_{33} V (L/t) $$

假定"完全自由变形"。检查FEM模型中是否有意外的位移约束边界条件或周边部件的接触条件。其次，实际执行器不仅产生

$$ d_{33} $$

模式，也产生横向应变（泊松比效应）。如果这被周边结构约束，见观位移会减少。最后，材料的非线性性。高电场下压电常数本身随电场变化，会比低电场值更小（饱和现象）。施加电压较高时，应考虑使用非线性压电材料模型。

压电效应 — CAE术语解释

压电效应

机械-电气转换的物理现象

压电效应的理论基础

基本概念和支配方程

压电效应的数值计算方法

FEM离散化和求解器设置

压电效应的实务应用

工作流程和检查清单

压电效应的软件比较

Ansys/Abaqus/COMSOL的特点

压电效应的故障排除

常见错误和对策

相关主题