在ANSYS中,接触刚度过大通常会导致收敛困难,甚至导致求解失败。这是因为过大的接触刚度会使系统的整体刚度矩阵出现病态,导致迭代求解过程不稳定或无法达到平衡。
接触刚度与收敛及计算精度的矛盾关系
接触刚度是罚函数法等接触算法中的关键参数,它直接影响了仿真结果的精度与求解过程的收敛性,两者之间存在需要权衡的矛盾关系。
精度要求:接触刚度越大,两个接触表面之间的允许穿透量就越小,计算结果的精度理论上越高。
收敛要求:然而,过大的接触刚度会使刚度矩阵病态,显著增加迭代次数,并极易引起求解震荡,从而导致收敛困难或失败。
平衡原则:因此,设置接触刚度的通用原则是:选择一个足够大的值以保证穿透量在可接受范围内,同时又足够小以避免引起收敛问题。
接触刚度过大导致收敛困难的具体机制与表现
刚度设置不当是导致非线性接触分析不收敛的常见原因之一,其影响体现在求解过程的多个层面。
引发矩阵病态与迭代激增:过大的接触刚度会直接导致总体刚度矩阵的条件数变差(病态),这使得求解器(特别是迭代求解器如PCG)需要更多的迭代次数,且常常无法满足收敛准则。
在收敛图中表现为力残差震荡:如果求解不收敛是由于接触刚度过大引起的,在力收敛图中,计算的力残差(紫色线)可能会在收敛容差(蓝色线)附近持续震荡,无法稳定地低于容差。
作为一系列收敛问题的触发因素:接触刚度过大可能与其他问题耦合,共同导致不收敛,例如加剧系统的力不平衡、或与快速加载、材料不稳定等问题叠加。
如何判断与调整接触刚度以促进收敛
当遇到收敛困难时,应系统性地检查并调整接触刚度,这是一个需要结合监控指标与经验尝试的过程。
诊断收敛问题的根源:检查穿透量与迭代次数:首先应后处理查看穿透量(如使用PLESOL,CONT,PENE命令),并检查求解信息中的平衡迭代次数。
区分原因:如果收敛困难主要由过大的穿透量引起,可能是接触刚度(FKN)被低估或允许穿透(FTOLN)设置过小;如果困难是由于力或位移残差需要过多迭代,则可能是FKN值被高估。
调整刚度的方法与经验值:尝试法(Trial-and-Error):对于使用罚函数法的接触,通常建议从一个较小的FKN值(例如0.1)开始计算,因为较小的值有助于初始收敛,然后逐步增加FKN进行一系列计算,直至获得满意的穿透结果。
参考经验因子:对于面-面接触单元,接触刚度(FKN)常定义为基体材料刚度的一个比例因子。初始估算时可参考:大面积实体接触使用1.0;对于弯曲主导的较柔软部件,使用0.01至0.1。
考虑更换接触算法或求解器:算法选择:不同接触算法的收敛性不同。例如,罚函数法通常收敛性最好,但需人工设置刚度;拉格朗日乘子法虽无需定义刚度,但可能因接触状态突变而导致收敛困难;增强拉格朗日法则是自动更新刚度的罚函数法。
求解器选择:对于因刚度问题导致的病态矩阵,可以考虑使用直接法求解器(如稀疏求解器)替代迭代求解器,因为直接法更能有效处理病态问题。
武汉格发信息技术有限公司 | 许可分析,许可优化,许可管理,许可授权,软件授权