一种基于同时消元回代法的大型线性方程组迭代求解法

来源 :2013中国力学大会 | 被引量 : 0次 | 上传用户:vgbin2
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
基于SEBSM,提出了一种求解超大型线性方程组的迭代求解法.该方法将系数矩阵分解为下矩阵和上矩阵两部分:下矩阵是迭代矩阵,沿主对角线具有一定的带宽,数据储存空间和迭代收敛速度完全由带宽来决定;上矩阵要比下矩阵有相当程度的小,在迭代中用于形成方程的右端项.由于在本方法中,所需数据储存空间只由所选带宽中的非零元素来决定,而不是由所有下矩阵中的非零元素决定,因此本文方法可以求解超大型工程问题形成的线性代数方程组.
其他文献
基于POD方法与区间数学理论将随机场分解为依赖于时间的区间不确定量和依赖于空间的协方差矩阵模态的线性组合.首先应用POD方法结合区间截断理论得到随机场协方差矩阵,然后采用Lancos方法提取协方差矩阵前几阶占优本征模态作为随机场的空间分布形式.根据计算流体力学方法得到不确定载荷,计算了三维机翼的位移响应,并对此方法的精度及效率进行了讨论.
会议
运用非惯性坐标系中的力学问题的理论来研究动力刚化典型实例,建议了在非惯性系中研究动力刚化问题的一种途径,讨论了正确进行零次建模的方法.本文合理地将典型实例发展为梁微元动力刚化模型.研究梁微元动力刚化问题的力学模型,得到具有明确物理意义的结果.以梁微元为例,研究动力刚化的零次建模.在此基础上,研究航天动力学中梁式太阳能帆板的动力刚化问题,得到问题的解析解.最后,研究航天动力学中板式太阳能帆板的动力刚
会议
工程优化最终结果有两种不能容忍的问题:(1)某些约束违背了约束限;(2)"最大约束"的数值远小于约束限.为了解决它们,以往的研究是:(1)寻求近似范围大且精度高的显函数;(2)提出运动极限在小范围内避免可行域被错估.本文突破思维定势的束缚,提出在寻优过程中自适应调整约束限的方法,控制迭代过程使"最大约束"成为等式约束,并使最优解满足所有约束.在位移约束下结构拓扑优化问题的求解中予以实现,计算量小,
会议
发展了一种形状自由的多节点裂纹奇异HSF单元,然后与已有的平面HSF单元组合,共同模拟断裂和裂纹扩展问题.通过数值实验,给出了建议的建模策略.数值算例和模拟结果表明,奇异HSF单元形状自由、精度高,整体计算网格单元数量很少且不受网格畸变影响,能够高效模拟准静态裂纹扩展问题.
会议
研究了附加阻尼材料层的壳体结构在受到简谐激励作用时阻尼材料的最优分布问题.采用带有惩罚的人工阻尼材料模型构造了拓扑优化模型,优化目标是在给定阻尼材料用量的情况下最小化结构的动响应.由于结构呈现非比例阻尼特性,采用了降阶的状态空间内的复模态叠加法计算结构的稳态响应,并给出了伴随变量法动响应灵敏度分析的计算格式.拓扑优化算例验证了优化模型和方法的正确性和有效性.
会议
基于三维质点法法程序MPM3D,通过减聚力模型描述材料失效界面的力学行为,发展了模拟延性断裂的有效工具.物质点法具有不需要固定网格连接的特性,相比有限单元法对模拟不连续问题具有较明显的优势.减聚力模型将断裂行为的描述引入到本构方程中,当材料应力达到阈值后破坏界面力通过减聚力模型卸载,可以考虑Ⅰ型、Ⅱ型以及Ⅰ型Ⅱ型混合失效.
会议
报告首先介绍百亿亿次计算机硬件设计面临的挑战,接着探讨新型计算机体系结构下高效数值算法应具有的特征,然后讨论领域应用如何应对来自极大规模的算法方面的挑战,最后介绍一些我们在国产千万亿次超级计算机以及Titan上的大规模计算实践与经验.
会议
介绍了这方面的最新研究进展,包括曲线网格上适合连续流域的高精度多维动理学格式的构造以及格式在可压缩湍流模拟中的应用.研究着重讨论了初场重构以及多维效应对格式精度的影响.研究结果表明了高阶动理学格式在可压缩流动精细模拟中的良好性能,为后续改进和优化工作奠定了基础.
会议
提出了一种结构振动分析的新型超收敛等几何方法.其基本思想是通过构造高阶质量矩阵来提高结构振动频率的计算精度.文中证明高阶质量可以通过一致质量矩阵和缩减带宽质量矩阵的最优组合进行构造,其频率计算精度比传统的一致质量矩阵提高两阶,具有几何精确和超收敛的特点.系列典型算例验证了理论结果.
会议
数值模拟在未知领域探索、物理本质认识等方面展现出巨大作用.随着科学问题的日益复杂,研发CAE软件有效地反应各种结果,重现多尺度下不同物理现象,已成为当前科学研究中亟待解决的问题.我国在计算力学程序开发方面做了诸多工作,但大多不具备良好的图形用户界面和便捷的功能拓展能力,难以形成通用性强、易被用户接受的高质量软件,计算力学软件市场仍被国外软件把持.我国为何研发不出优秀的商业有限元软件昵?
会议