论文部分内容阅读
采用有限元方法解决结构静力学问题的关键是线性方程组的求解。随着结构的复杂化和单元划分的细致化,完全依靠内存求解线性方程组的方法已不能满足要求。因此,求解巨大方程组的算法研究成为当前的一个研究热点。本文在综述钻柱静力学的发展、钻柱有限元研究现状、有限元分析中C语言应用研究现状的基础上,研究了国内外使用分块思想求解大型线性有限元方程组的算法,重点研究了分块算法和大型有限元软件利用内存调用外存的波前法,分析了两种算法的优缺点。以减少元素存储量和提高计算效率为目的,提出了分布算法。分布算法采用一维变带宽存储方式存储刚度矩阵。在存储外存文件时,不是像分块算法那样,而是将下三角矩阵组成的一维数组的指定个数调入外存文件,而后运用调用函数对外存文件进行调进和调出,减少内存的需求量和节点编号不合理的影响;外存文件存储下三角矩阵时,无法用一般高斯消去法求解方程组。针对该问题,参考改进的平方根求解方法,实现了边调用外存文件边进行改进的乔里斯基分解,分解完成后再放入外存文件,最后反向调入外存文件进行回代求解,达到了内外存交换次数较少、程序实现简单的目的。以钻柱为研究对象,用C语言编制了分布算法程序,研究了钻柱在多种载荷作用下的弹性变形情况,并采用有限元分析软件ANSYS的参数化设计语言APDL对钻铤与钻杆结合处以上小段进行了仿真,对计算程序进行了验证。数值实验证明,分布算法具备分块算法和波前法的优点,编程方便、计算结果正确,并且随着单元划分的扩大,存储效率越来越大,计算用时随结点总数按抛物线规律分布,分析得出分布算法的计算速度是一般算法的4倍左右。本文实现的算法编程对于打破国外软件的垄断具有一定的理论价值。