收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

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

肖映雄  
【摘要】:随着计算机日新月异的发展,科学计算特别是大规模、高性能科学计算越来越成为推动技术革命的强劲动力,计算方法是科学计算的核心,因此,计算机与计算方法的发展程度已成为科学计算能力提高的决定性因素。在物理、力学等应用领域里,有很多问题常常归结于偏微分方程(组)的求解,但由于实际问题的复杂性,根本无法得到它们的解析解,有限元数值求解已成为求解这类问题最为有效的途径和方法。有限元方法通常包括三个过程,即网格生成及优化、有限元离散代数系统的形成以及离散系统的代数求解,其中第一个和第三个过程是影响有限元分析整体求解效率的主要因素。多重网格法是求解偏微分方程(组)大规模离散化方程最为有效的方法,它一般可分为几何多重网格法和代数多重网格法。由于实际应用问题的错综复杂性,以及数值商业软件对“即插即用”型求解器的要求,使得几何多重网格法的应用变得越来越困难,而代数多重网格(AMG)法的“高效性”和“鲁棒性(robustness)”,使之成为了当今多重网格法的研究热点。 AMG方法的关键技术是建立生成各个粗网格层和相应的插值算子的代数方法。对大量的偏微分方程离散代数系统,其背景问题的许多重要特征仅仅通过总刚度矩阵是很难重构出来的。因此,借助于部分几何或分析信息,再通过代数途径来构造相应的高效AMG法是一种十分自然的想法,我们称这种AMG法为基于部分几何和分析信息的AMG法。目前,对标量椭圆型偏微分方程,AMG方法发展比较成熟,但对方程组情形,通常的AMG方法在求解效率上将变差,有时甚至失效。其主要原因是,对方程组情形,定义在每一个节点的自由度往往大于1,而且这些自由度又是相互耦合在一起的,总刚度矩阵所对应的代数网格图与相应的几何网格图根本不一致,利用通常的网格粗化技术,很难控制粗网格的规模及合理地设计插值(提升)算子,从而大大降低了其求解效率。特别是对三维情形,由于网格图的复杂度较二维情形有本质性的增加,相应的网格粗化技术在计算效率和鲁棒性等方面还不尽如人意。因此,我们需要发展新的网格粗化技术及插值算子的构造方法,以提高AMG方法的收敛速度,使其能应用于更多实际问题的快速求解。 本文主要针对固体力学计算中的几类应用问题,对其相应的AMG算法进行了深入的研究和探讨,借助于部分几何或分析信息,提出了几种有效的AMG算法,并进行了大量的数值实验与结果比较,得到了一些有意义的数值结果。这些研究进一步丰富和充实了AMG算法,拓宽了AMG方法在一些应用领域中的研究,具有重要的理论和工程应用价值。主要内容和结果包括: (一) 针对一类(多尺度)离散应变原子模型的数值求解,将基于标量椭圆型偏微分方程的AMG方法推广应用于方程组情形,给出了相应的AMG算法,并详细介绍了其中的网格粗化算法及插值(提升)算子的构造。据我们所知,这是首次尝试设计快速方法求解离散应变原子模型。对二维问题作了大量的数值实验,并与工程计算中常用的数值方法进行了比较,结果表明,本文设计的AMG算法特别是AMG-CG方法对求解应变模型(包括多尺度耦合模型)是有效的,具有很好的“鲁棒性”。 (二) 针对含间断系数弹性结构力学问题的数值求解,建立了一类界面保持粗化多重网格方法,这样,只需构造简单的插值算子及选取点块Gauss-Seidel作磨光迭代,就可得到一类相当有效的多重网格方法。数值结果表明,这种界面保持粗化多重网格方法


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 黄佩奇;陈金如;;MORTAR型旋转Q_1元的多重网格方法[J];计算数学;2008年01期
2 张士杰,袁新,叶大均;全三维粘性流场计算的快速收敛方法[J];清华大学学报(自然科学版);2000年10期
3 蔡宗熹;张培琴;张庆合;蔡越虹;;二重网格方法用于曲面位场转换[J];物探化探计算技术;1992年04期
4 牟斌;肖中云;周铸;陈作斌;刘刚;;多重网格技术在复杂粘性流场计算中的应用及研究[J];空气动力学学报;2006年01期
5 刘伊克,常旭;盆地模拟水动力油气二次运移隐式多重网格法[J];地球物理学报;1998年03期
6 王怀柱;;三维非稳态导热问题的高效稳定数值解法[J];河北北方学院学报(自然科学版);2009年05期
7 姜亚琴;;一类Mortar型旋转Q_1元的多重网格方法[J];南京师大学报(自然科学版);2011年02期
8 张培昌,李建通,顾松山;用多重网格法解赫姆霍兹型欧拉方程[J];南京气象学院学报;1995年02期
9 张正科,李凤蔚,罗时钧;欧拉方程多重网格方法及三维跨声速大迎角流动数值计算[J];空气动力学学报;1996年04期
10 马廷福;金涛;葛永斌;;三维对流扩散方程的高精度多重网格方法[J];赤峰学院学报(自然科学版);2009年12期
11 韩研;陈金如;;平面线弹性纯位移边值问题低阶非协调矩形元的Robust多重网格方法[J];高等学校计算数学学报;2011年02期
12 葛永斌;曹富军;马廷福;郝庆一;;三维对流扩散方程的多重网格算法研究[J];计算机工程与应用;2009年13期
13 高彦伟,郭华,周伟;利用多重网格方法研究波动方程正反演问题进展[J];世界地质;2002年01期
14 岳锋;李凤蔚;;三维Euler/N-S方程的多重网格求解方法[J];民用飞机设计与研究;2005年02期
15 葛永斌,田振夫,吴文权;三维波动方程的隐式多重网格方法[J];华中师范大学学报(自然科学版);2005年01期
16 梅立泉;张明;;各向异性椭圆方程的多重网格局部Fourier分析[J];工程数学学报;2006年05期
17 尹凤娥;;多重网格方法及其在高超声速进气道中的应用[J];航空计算技术;2010年02期
18 蒋长锦;;常微分方程初值问题的多重网格解法[J];中国科学技术大学学报;1989年04期
19 张淼;陈庆文;;基于套迭代技术的弹性梁动力响应分析的多重网格方法[J];长春工程学院学报(自然科学版);2008年02期
20 杨中;黄典贵;徐建中;;低速流动的预处理多重网格求解[J];工程热物理学报;2008年06期
中国重要会议论文全文数据库 前10条
1 陈小宏;牟永光;;地震波动方程反演的多重网格方法[A];1995年中国地球物理学会第十一届学术年会论文集[C];1995年
2 徐建东;陈国定;李东紫;;等温线接触弹流问题的两重网格法分析[A];第二届全国青年摩擦学学术会议论文专辑[C];1993年
3 刘伊克;;油气二次运移隐式多重网格法[A];1997年中国地球物理学会第十三届学术年会论文集[C];1997年
4 杨爱明;翁培奋;乔志德;;悬停旋翼跨声速无粘流场计算中的多重网格方法[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
5 葛永斌;田振夫;冯秀芳;;二维对流扩散方程的高精度紧致隐格式及其多重网格算法[A];第七届全国水动力学学术会议暨第十九届全国水动力学研讨会文集(下册)[C];2005年
6 王方飞;孟惠荣;;点接触弹性变形的多重网格算法[A];第二届全国青年摩擦学学术会议论文专辑[C];1993年
7 谢剑斌;王晖;刘通;李沛秦;;基于多重网格的光流场快速计算方法[A];第六届全国信息获取与处理学术会议论文集(2)[C];2008年
8 张×;周立伟;金伟其;张智诠;;多重网格法求解三维静电场分布[A];面向21世纪的科技进步与社会经济发展(上册)[C];1999年
9 何光洪;乔志德;;多重网格法在三维Euler方程/NS方程求解中的应用[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
10 阳莺;;一类静电势问题的多重网格算法[A];第八届中国青年运筹信息管理学者大会论文集[C];2006年
中国博士学位论文全文数据库 前10条
1 肖映雄;代数多重网格算法研究及其在固体力学计算中的应用[D];湘潭大学;2006年
2 鲁晶津;地球电磁三维数值模拟的多重网格方法及其应用研究[D];中国科学技术大学;2010年
3 刘群锋;最优化问题的几种网格型算法[D];湖南大学;2011年
4 张淼;结构动力响应分析的多重网格方法研究[D];吉林大学;2008年
5 张东良;起伏地表条件下的直流电场有限差分数值模拟[D];吉林大学;2011年
6 黄筱云;空间自适应的水波数值模型研究[D];天津大学;2012年
7 王立鹏;无单元法若干技术研究及其在电磁场计算中的应用[D];沈阳工业大学;2010年
8 郭荣文;贝叶斯MT反演的非线性和不确定度分析[D];中南大学;2011年
9 王德华;小波分析及其应用的几个问题[D];湖南大学;2006年
10 李莉;多重网格法在SF_6高压断路器电场与气流场数值计算中应用的研究[D];沈阳工业大学;2008年
中国硕士学位论文全文数据库 前10条
1 王燕龙;多重网格方法的大地电磁正演模拟与分析[D];中南大学;2011年
2 吴东阳;多重网格方法在工程电磁场数值计算中的应用研究[D];沈阳工业大学;2009年
3 王丫;基于小波的多重网格方法[D];中南大学;2008年
4 包常利;多重网格方法在风轮机流场数值模拟中的应用及研究[D];沈阳工业大学;2011年
5 乔中林;直流电法三维正演多重网格算法研究[D];中国地质大学(北京);2007年
6 张磊;Mortar元的多重网格方法[D];南京师范大学;2003年
7 赵钍焱;求解二维弹性力学问题高次有限元方程的代数多重网格法[D];湘潭大学;2006年
8 王军会;空气净化机高压静电场数值分析[D];华南理工大学;2010年
9 姜丽丽;基于多重网格方法的焊接熔池中耦合场的数值模拟研究[D];沈阳工业大学;2012年
10 王骁;基于多重网格方法的变分光流计算[D];西安电子科技大学;2008年
中国重要报纸全文数据库 前10条
1 本报记者 胡宇芬 通讯员 赵猛 蒋海文;寝室也是好课堂[N];湖南日报;2005年
2 瞿鹤鸣 赵猛;激情跨越[N];中国教育报;2003年
3 徐梦和;为了伟人的嘱托[N];湖南日报;2005年
4 黄海;彩虹总在风雨后[N];中国档案报;2004年
5 本报记者 杨期仁;湘潭大学学术成果丰硕[N];湘潭日报;2007年
6 记者 陈晓东 李伦娥;宁夏大学湘潭大学学生住得不错[N];中国教育报;2000年
7 赵湘明;新的跨越[N];光明日报;2003年
8 记者 唐湘岳 通讯员 蒋海文;湘潭大学:辅导员、班主任为学生成长成才导航[N];光明日报;2005年
9 刘文韬;高校里的“拥军书记”[N];湖南日报;2007年
10 祁金华;湘潭大学院系共建HPC系统[N];网络世界;2007年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978