接触颤振（状态振动） — 故障排除
故障排除步骤
颤振出现时，应该按什么顺序具体检查？
首先确认收敛历史。对于Abaqus，可以在.sta文件中追踪接触状态变化次数，对于Ansys，可以在solution.out中查看。
步骤1：接触对的确认
特定发生颤振的接触对。对于Abaqus，接触状态变化记录在.msg文件中SEVERE DISCONTINUITY发生位置。对于Ansys，可以使用CNCHECK命令可视化接触状态。
如果只在特定对中发生，那就可以集中修复那里。
步骤2：网格质量的确认
确认接触面网格。检查以下内容：
- slave面的单元大小是否比master面更细密
- 接触面是否有扭曲单元（宽高比 < 5:1）
- 节点是否集中在相对面的单元边缘上
步骤3：罚函数刚度的调整
对于Abaqus，使用*SURFACE BEHAVIOR, PRESSURE-OVERCLOSURE=LINEAR显式指定罚函数刚度。对于Ansys，逐步提高FKN从0.01开始。Marc 则调整CONTACT TABLE的BIAS因子。
降低太多会不会产生穿透？
好问题。降低罚函数刚度确实会增加接触面贯穿量（penetration）。监测COPEN（接触开口）和CPRESS（接触压），确保贯穿量保持在单元大小的10%以下。
步骤4：稳定化阻尼的引入
如果仍未解决，引入接触稳定化。对于Abaqus：
```
*CONTACT CONTROLS, STABILIZE=0.001, TANGENTIAL FACTOR=0.001
```
对于Ansys：
```
STABILIZE,ON,ENERGY,0.001
```
0.001这个值是从哪里来的？
这是经验初值。值过大会影响结果，过小则无稳定化效果。在调整散逸能量比例时检查。如果散逸能量比小于1%，结果影响可以忽略。
预防措施
*CONTACT PAIR, ADJUST预先调整初始接触状态接触问题看似平凡，但实务中最容易卡壳。今后我会先从统一网格开始。
花在求解器错误原因查找上的时间应该更少。 — Project NovaSolver以改进错误诊断体验为研究主题之一。
请告诉我们您在接触颤振（状态振动）实务中遇到的问题
Project NovaSolver致力于解决CAE工程师日常面临的问题——设置的复杂性、计算成本、结果解释——以实现更好的工具。您的实务经验将成为开发更好工具的动力。
联系我们（准备中）相关主题
价值
详细
错误