Abaqus压电仿真:三个矩阵输错一个就白算,单位混用结果差出好几位

Abaqus 多场耦合

超声换能器、声呐、医用超声探头、加湿器、眼镜清洗机……这些产品的核心都是压电陶瓷:给它力能生电(压电效应),给它交变电压又能产生同频机械振动(逆压电效应)。做这类多物理场仿真,最难的不是画几何,而是各向异性材料的三个矩阵怎么输、单位怎么统一——很多人卡在这里,矩阵一错结果毫无意义。

这篇以「压电换能器粘接钢板」案例为背景,把 Abaqus 压电分析最关键的 5 点讲清楚:压电效应与换能原理、压电方程与三个矩阵、矩阵转置与下标对应、为什么必须用国际单位制、几何建模要点。


要点一:正逆压电效应,是换能器工作的物理基础


压电效应由材料微观晶体结构决定:对压电陶瓷施加机械应力(挤压/拉伸),材料变形的同时会产生电荷、形成电势差——力生电,即正压电效应;反过来施加电压(电势差),材料产生变形——电生力/变形,即逆压电效应。

工程上利用逆压电效应:在两片对贴的压电陶瓷(中极为正、上下为负)上施加交变电压,结构进入稳态后产生与激励同频率的正弦振动,形成超声波,用于声呐探测、医学成像、无损检测、超声清洗与加湿。材料有发射型(如 PZT-4)与接收型之分,需按场合选型。


要点二:压电本构由两个方程、三个矩阵描述


普通弹性材料只有一个本构关系:应力 = 刚度矩阵 × 应变(应力、应变各 6 个分量)。压电压多了电场耦合,由一对方程联立描述:应力方程中,应力一部分来自应变、一部分来自电场(减去压电系数矩阵 × 电场强度);电位移方程中,电位移 = 压电系数矩阵(转置形式)× 应变 + 介电系数矩阵 × 电场强度。

因此需要输入三个材料矩阵:①刚度(弹性)系数矩阵——因陶瓷各向异性,弹性用矩阵而非单一弹性模量;②压电系数矩阵——描述电场与应力的耦合倍数,大部分元素为 0,只在特定位置非零;③介电常数矩阵——对角阵,只需输入 11、22、33 三个对角项。

压电三个材料矩阵

Abaqus压电仿真:三个矩阵输错一个就白算,单位混用结果差出好几位的图2


图 1 压电本构与三个材料矩阵

踩坑点:刚度、压电、介电三个矩阵缺一个或输错一个,力电耦合就不成立,结果不可信。


要点三:文献数据要转置,下标对应关系必须搞清


查文献时常拿到 6×3 的压电系数矩阵,而 Abaqus 里输入的是其转置形式,不转置直接填会错位。可以借助软件中的下标理解:如「311」表示第 3 行、对应第 11 列的位置,「322」「333」依次对应;转置后行列对应关系就清晰了,与应力分量(σ11、σ22、σ33、τ12、τ13、τ23)一一对应。

踩坑点:文献矩阵与软件输入表排布方向不同,照抄不转置是压电仿真最常见的隐性错误,软件不一定报错但结果错。


要点四:带电学量必须用国际单位制,避免换算混乱


这是案例特别强调的一点:本案例全程采用国际单位制(m、kg、s、Pa)建模。原因是压电参数的单位是库仑每平方米(C/m²)、法拉每米(F/m)这类带电学单位,一旦改用毫米等工程单位,长度、力与电学量混在一起,小数点差多少位很容易算乱。

用国际单位制时,原始参数可直接录入、无需换算,例如:钢板 100mm×100mm×2mm 建成 0.1m×0.1m×0.002m;不锈钢弹性模量 200GPa 直接写 200×10⁹ Pa、泊松比 0.3、密度 7900 kg/m³。

国际单位制建模

Abaqus压电仿真:三个矩阵输错一个就白算,单位混用结果差出好几位的图4


图 2 以米为单位的几何建模

踩坑点:几何用毫米、材料用 GPa、电学量又用国际单位,三套单位混用必然出错;带电学耦合优先锁定 SI 单位制。


要点五:换能器几何按实际分层建模,回转件用旋转


模型组成:底部不锈钢板(100×100×2mm)、两片 PZT 陶瓷环片(外径 Φ37、内径 Φ12、每片厚 5mm,对贴共 10mm)、上部不锈钢实心头(Φ37、厚 20mm)、铝合金喇叭形放大段(旋转生成,端部 Φ45)、以及把上下夹紧的中心螺栓(Φ10、高 10mm),各部分最终粘接为一体。

建模要点:等截面圆环用拉伸(内径 0.006m、外径 0.0185m、厚 0.005m);喇叭形、螺栓等回转体用旋转 360°生成;建模时先设置工作目录(如 EX15)。金属件按各向同性赋参数,PZT 陶瓷单独走压电材料输入流程(下一阶段完成材料矩阵赋予与电边界)。

压电陶瓷环片建模

Abaqus压电仿真:三个矩阵输错一个就白算,单位混用结果差出好几位的图6


图 3 PZT陶瓷环片尺寸与拉伸建模

踩坑点:软件显示精度不足时尺寸数字可能「看着不对」,以实际输入的坐标值(如 0.006、0.0185)为准,不要被显示误导。


小结


①正压电效应力生电、逆压电效应电生振,换能器靠交变电压驱动同频振动;②压电压需输入刚度、压电系数、介电常数三个矩阵;③文献的 6×3 压电矩阵要按软件要求转置、对清下标;④带电学量务必用 SI 国际单位制,参数直录不换算;⑤分层建几何、圆环拉伸、回转体旋转,金属各向同性、陶瓷走压电输入。矩阵与单位这两关把住,压电仿真就成功了一半。


【系统学习 Abaqus 结构与多场耦合】


从静力、非线性到压电、热-结构等多物理场耦合,Abaqus 结构仿真体系课程配合工程师能力认证与真实项目实训,带你打通材料本构、多场耦合与结果评估。机械工程师转型仿真、想进阶多物理场方向可重点了解。


【配套资料领取】


扫码添加小助理,私信回复关键词【Abaqus压电】,领取本篇配套资料包:①压电效应与三个矩阵核心知识总结 ②材料矩阵与SI单位关键参数速查表 ③换能器分层建模实战案例(分步) ④矩阵转置与单位易错清单。

Abaqus压电仿真:三个矩阵输错一个就白算,单位混用结果差出好几位的图7
查看全文
我的头像
默认 最新
ansys结构交流群