对于结构稳定性分析中超大规模矩阵正定性判定, 必须采用并行计算方法, 传统方法如计算特征值、主子式行列式及LDLT等直接方法难以实现. 本文提出了一些可适用于并行的迭代判定算法. 借鉴力学系统中能量下降的思想, 发展了一种判定矩阵正定性的新思路, 即将矩阵的正定性判定问题转化为一个优化问题, 并基于优化算法来判定矩阵的正定性. 提出了基于最速下降法和共轭梯度法来进行矩阵正定性判定的算法. 然后考虑到力学系统刚度矩阵的稀疏性和结构刚失稳状态的弱非正定性, 提出可以先截超平面后解方程求驻值点的方法来判定弱非正定矩阵的正定性. 为了保证对强非正定矩阵判定的准确性, 本文提出可以高效混杂使用截平面法和共轭梯度法. 数值实验结果表明, 本文提出的算法具有准确性和高效性. 对于强非正定矩阵而言, 共轭梯度算法更加高效; 而对于弱非正定矩阵, 则是截平面法和混杂算法更加高效. 这些算法都容易在机群上实现并行计算, 能够快速判定大规模矩阵的正定性.
2017, 49(6): 1223-1230.
doi: 10.6052/0459-1879-17-292
徐辉, 刘彬. 固体结构大规模矩阵正定性判定的快速算法. 力学学报, 2017, 49(6): 1223-1230. doi: 10.6052/0459-1879-17-292.