代数多重网格算法研究及其在固体力学计算中的应用

来源 :湘潭大学 | 被引量 : 10次 | 上传用户:frog_t
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
随着计算机日新月异的发展,科学计算特别是大规模、高性能科学计算越来越成为推动技术革命的强劲动力,计算方法是科学计算的核心,因此,计算机与计算方法的发展程度已成为科学计算能力提高的决定性因素。在物理、力学等应用领域里,有很多问题常常归结于偏微分方程(组)的求解,但由于实际问题的复杂性,根本无法得到它们的解析解,有限元数值求解已成为求解这类问题最为有效的途径和方法。有限元方法通常包括三个过程,即网格生成及优化、有限元离散代数系统的形成以及离散系统的代数求解,其中第一个和第三个过程是影响有限元分析整体求解效率的主要因素。多重网格法是求解偏微分方程(组)大规模离散化方程最为有效的方法,它一般可分为几何多重网格法和代数多重网格法。由于实际应用问题的错综复杂性,以及数值商业软件对“即插即用”型求解器的要求,使得几何多重网格法的应用变得越来越困难,而代数多重网格(AMG)法的“高效性”和“鲁棒性(robustness)”,使之成为了当今多重网格法的研究热点。   AMG方法的关键技术是建立生成各个粗网格层和相应的插值算子的代数方法。对大量的偏微分方程离散代数系统,其背景问题的许多重要特征仅仅通过总刚度矩阵是很难重构出来的。因此,借助于部分几何或分析信息,再通过代数途径来构造相应的高效AMG法是一种十分自然的想法,我们称这种AMG法为基于部分几何和分析信息的AMG法。目前,对标量椭圆型偏微分方程,AMG方法发展比较成熟,但对方程组情形,通常的AMG方法在求解效率上将变差,有时甚至失效。其主要原因是,对方程组情形,定义在每一个节点的自由度往往大于1,而且这些自由度又是相互耦合在一起的,总刚度矩阵所对应的代数网格图与相应的几何网格图根本不一致,利用通常的网格粗化技术,很难控制粗网格的规模及合理地设计插值(提升)算子,从而大大降低了其求解效率。特别是对三维情形,由于网格图的复杂度较二维情形有本质性的增加,相应的网格粗化技术在计算效率和鲁棒性等方面还不尽如人意。因此,我们需要发展新的网格粗化技术及插值算子的构造方法,以提高AMG 方法的收敛速度,使其能应用于更多实际问题的快速求解。   本文主要针对固体力学计算中的几类应用问题,对其相应的AMG 算法进行了深入的研究和探讨,借助于部分几何或分析信息,提出了几种有效的AMG 算法,并进行了大量的数值实验与结果比较,得到了一些有意义的数值结果。这些研究进一步丰富和充实了AMG 算法,拓宽了AMG 方法在一些应用领域中的研究,具有重要的理论和工程应用价值。主要内容和结果包括:   (一)针对一类(多尺度)离散应变原子模型的数值求解,将基于标量椭圆型偏微分方程的AMG 方法推广应用于方程组情形,给出了相应的AMG 算法,并详细介绍了其中的网格粗化算法及插值(提升)算子的构造。据我们所知,这是首次尝试设计快速方法求解离散应变原子模型。对二维问题作了大量的数值实验,并与工程计算中常用的数值方法进行了比较,结果表明,本文设计的AMG 算法特别是AMG-CG 方法对求解应变模型(包括多尺度耦合模型)是有效的,具有很好的“鲁棒性”。   (二)针对含间断系数弹性结构力学问题的数值求解,建立了一类界面保持粗化多重网格方法,这样,只需构造简单的插值算子及选取点块Gauss-Seidel 作磨光迭代,就可得到一类相当有效的多重网格方法。数值结果表明,这种界面保持粗化多重网格方法的收敛性不依赖于网格规模及间断系数的大小,具有很好的数值稳定性,非常适合于含规则界面区域弹性力学问题的数值求解。然后,对不规则界面情形,通过从纯代数角度构造相应的界面保持粗化算法,或者通过选取合适的光滑迭代子的方法,大大提高了AMG 方法的计算效率,并使之适合于求解任意形状的夹杂问题。本研究成功地解决了将AMG 方法应用于非均匀介质(如:夹杂,局部软化)这一挑战问题。   (三)针对三维弹性问题的数值求解,建立了一类三维等代数结构面网格剖分下的AMG 算法和相应的预处理共轭梯度法。这种方法的基本思想是将三维插值算子的构造本质地转化为二维插值算子的构造,开辟了新的算法设计途径。数值结果表明,本文建立的AMG 方法对求解三维弹性问题有限元方程是十分有效的,具有很好的“鲁棒性”,较直接解法SuperLU和其它常用迭代方法具有明显的优越性。   (四)将AMG 方法引入到岩体力学有限元计算领域,论述了基于单元聚集和能量极小意义下适于岩体力学有限元求解的AMG 粗化策略与插值算子,并详细描述了相应的AMG 算法。本研究工作进一步开拓了AMG 方法应用研究的新领域,通过将其应用到采矿与岩土工程有限元分析中,提高了有限元模拟岩土工程大规模问题的数值效率,有助于解决采矿与岩土工程实践中不断出现的大规模科学计算问题,具有重要的科学价值和重大的工程应用价值。数值试验表明:在岩体力学与工程问题的有限元数值计算中,AMG 求解法是高效的、适用的,较直接法和其它常用迭代方法具有明显的优越性。   (五)针对线弹性中半正定问题的数值求解,首次提出了一类能重构刚体模式的AMG算法。将所建立的AMG 方法应用于一类带牵引自由边界条件的应变问题和一类多尺度应变原子模型的数值求解,结果表明,本文提出的AMG 方法较一般的AMG 方法具有更好的收敛特性。同时,利用正交分解技术,我们还得到一类能保持商空间的AMG方法,当原问题的真解落在所考虑算子的商空间时,由该方法求得的数值解也完全落在相应的商空间里,且收敛于原问题的真解。   (六)针对线弹性中的几类典型问题,研究了几种有代表性的复杂各向异性网格(包括贯穿型各向异性网格和具有局部各向异性的混合网格)对AMG 方法的影响;进而,通过利用问题和网格的部分几何和分析信息,引入“距离矩阵”的概念,设计了一种新的网格粗化算法,数值实验表明:基于新的网格粗化算法的AMG 法对求解线弹性中的各向异性问题是非常有效的,具有很好的数值稳定性。
其他文献
结构优化设计是近四十年发展起来的一门新技术,与传统设计相比,由于其具有设计合理、节省原材料等特点日益受到国内外工程界的广泛重视。结构优化设计可以降低结构重量和材料成
目的 探讨高迁移率组蛋白B1(HMGB1)与Ⅲ期直肠癌患者临床病理特征及术后3年无疾病生存(DFS)的关系.方法 选取我院2011年2月至2015年12月诊治的Ⅲ期直肠癌术后患者组织标本46
超音速客机巡航时,舵面偏转会产生铰链附近的激波诱导分离现象,由于分离流动会降低舵面效率,因此该现象在超音速客机设计中有重要的研究价值。本文通过CFD数值方法,研究了超音速
复合材料梁在飞机结构中的应用越来越广,而常常由于梁腹板失稳而发生破坏,因此对梁腹板的稳定性进行分析以及其对其破坏载荷进行预报是非常必要的。由于复合材料是一种多相材
为适应煤矿安全生产的要求,结合煤矿生产实际,本研究开发了一套井下人员射频识别定位管理系统、该系统具有性能稳定,人员识别率高,误码率小,使用方便,成本低等优点,对于煤矿安全生产
美国军用数据总线标准MIL-STD-1553B是一种集中控制式、飞机内部时分指令/响应型多路串行数据总线标准,具有高可靠性和灵活性,已经成为现代航空机载系统设备互联的关键技术,广泛
微梁谐振器由于其尺寸小,比表面积大,在振动的过程中,压膜阻尼会显著影响谐振器的动态响应特性,因此研究微梁谐振器中的气体压膜阻尼对提高器件的工作性能有重要意义。在研究微谐
基于重叠网格技术,本文应用有限体积法数值求解三维定常Navier-Stokes方程,没有附加任何尾迹模型,模拟了直升机旋翼在悬停状态下的粘性流场。同时,结合:Kirchhoff方程,研究了旋翼在
2009年8月21日,由中国工商银行票据营业部和《金融论坛》编辑部共同举办的票据市场发展的挑战与机遇研讨会在上海成功举办。人民银行、银监会、商业银行、高校研究机构的专家
风力发电是一种新型的绿色能源,正成为世界各国争相开发的新技术领域。随着工业的发展,随即出现供电环境恶化以及电力紧缺的问题,并日渐严重。针对此问题,各国都制定了相关政策,对风能进行深入分析研究,其可再生性可以有效解决上述问题,在未来发展过程中将发挥出重要作用。当前,变速双馈风力发电技术应用较为广泛,由于其成本低、效率高等特点成为当今风力发电技术中的主流。风力机含有较多的旋转机械,而且风力发电受自然风