摘要:为提高求解含裂纹压电体的能量释放率的计算效率和精度以及减少程序编写与调试的工作量，针对四节点平面压电单元提出了机电虚拟裂纹闭合法计算公式.基于通用商业有限元软件ABAQUS，开发了哑节点断裂压电单元，编写用户自定义子程序UEL，该单元可独立求解总能量释放率分量，对不同材料和裂纹长度的含裂纹压电体的能量释放率分量进行了求解，与理论解作了对比，讨论了不同形式网格对求解精度的影响.数值算例结果表明，哑节点断裂压电单元具有精度高、简单方便和对网格尺寸不敏感等优点.