收藏本站
《上海交通大学》 2010年
收藏 | 手机打开
二维码
手机客户端打开本文

四阶椭圆型方程若干有限元新方法和高效求解算法

黄学海  
【摘要】: 四阶椭圆型微分方程广泛应用于固体力学、材料科学和图像处理诸领域,因此对它的数值解研究不但具有重要的理论意义也具有直接应用价值.本文的主要工作是构造了重调和方程的一种基于Poisson求解子的快速求解器;对四阶椭圆偏微分方程的MWX元方法设计了两网格局部和并行求解算法并进行了误差分析;将C0间断有限元方法应用于四阶椭圆偏微分方程并给出求解由该方法离散得到的线性代数方程组的区域分解法;也提出了自适应有限元方法的一个抽象框架. 首先,构造了重调和方程的一种基于Poisson求解子的快速求解器.通过位能极小原理,建立了重调和方程和Stokes方程之间的等价性,即重调和方程等价于一个Stokes方程和两个Possion方程.类似地,在有限元离散情形,建立了重调和方程Morley元方法和Stokes方程非协调P1-P0元方法之间的等价性,即重调和方程Morley元方法等价于一个Stokes方程非协调P1-P0元方法和两个Possion方程Morley元方法.然后基于这个等价性,结合代数多重网格法,构造了重调和方程的一种基于Poisson求解子的快速求解器.最后,用数值试验验证了这个求解器的有效性. 其次,提出了三种用于求解任意空间维数的四阶椭圆偏微分方程的两重网格局部和并行的Morley-Wang-Xu (MWX)元方法.由于MWX元空间不是嵌套的,需要引进一些网格转移算子,以从全局的粗网格有限元解和局部的细网格修正来获得改进的全局有限元解.首先通过构造基于修正Argyris元的网格转移算子,提出了求解MWX元离散的四阶问题的三种局部和并行算法.然后,证明了相应的数值解的离散能量误差的阶为O(h + H2),其中H和h分别为有限元单形剖分的粗细网格的网格大小.同时,构造了基于算术平均网格转移算子的局部和并行算法,证明了相应的数值解的离散能量误差的阶为O(h + H2(H/h)(d?1)/2),其中d为求解区域的空间维数.进一步,我们提供了数值实验来说明这些算法的计算效果. 再次,提出了求解薄板弯曲问题的一类新的C0间断有限元方法,并给出了数值解在离散能量误差和H1范数下的误差估计.首先将四阶偏微分方程写成二阶系统的形式,获得了构造求解原四阶问题的C0间断有限元方法的框架.然后,建立一个离散的稳定性等式,并由此选择可行的数值迹,得到求解薄板弯曲问题的一类稳定的C0间断有限元方法.根据二阶椭圆问题间断有限元方法误差分析中的思想和推导技巧,我们得到了薄板弯曲问题LCDG方法数值解在离散能量范数和H1范数下的最优阶误差估计以及CDG方法数值解在离散能量范数和H1范数下的误差估计.最后,用数值实验验证了理论的收敛阶. 接着,构造了薄板弯曲问题的局部C0间断有限元方法的两水平加型Schwarz预条件子并估计了条件数.对于我们提出的求解薄板弯曲问题的LCDG方法,通过定义特殊的网格转移算子,给出了该方法的区域分解算法.进一步估计了经过预处理的线性代数方程组的系数矩阵的条件数,该上界估计在δ≈H时是最优的.在小部分重叠情形,即δ H时,利用更多的子区域的信息,改进了条件数的上界估计. 然后,提出了基于标准循环的自适应有限元方法的抽象框架并得到了拟最优收敛性和最优复杂度.首先提出一些假设,基于这些假设得到自适应有限元方法的拟最优收敛性.在作了进一步的假设后,通过引进总误差,证明了自适应有限元方法的最优复杂度.最后通过前面提出的假设,将这个抽象框架应用于各种问题,包括一般二阶椭圆偏微分方程的高阶有限元方法, 2m阶椭圆偏微分方程的Morley型有限元方法,不定时谐Maxwell方程和H(div)方程. 最后,对于薄板弯曲问题k = 0,1时的混合有限元方法,验证了自适应有限元方法抽象框架中提出的假设.利用Helmholtz分解,对所有的k≥0,构造了薄板弯曲问题的混合有限元方法(H-H-J方法)基于残差的后验误差估计子并证明了估计子的可靠性.利用泡函数技巧分析了误差估计子的有效性.而且根据力在内边上的跳跃可以被体积部分所控制的性质,改进了误差估计子.进一步,获得了离散的Helmoltz分解和离散的inf-sup条件,利用这两个工具和后验误差分析中得到的性质,建立了拟正交性和离散可靠性.对于其它的假设,在k = 0,1时容易得到.
【学位授予单位】:

知网文化
【相似文献】
中国期刊全文数据库 前20条
1 林鹏程;解多维热传导方程的网格法[J];福州大学学报(自然科学版);1962年02期
2 谢惠民,袁杰,戴福隆;高温全息光栅的制作工艺研究[J];光学技术;1993年03期
3 高毅;高尚;;均匀设计优化方法[J];科学技术与工程;2011年02期
4 张启富;网格法在多极值规划问题中的应用[J];运筹与管理;1998年02期
5 王连生;牛昭芳;董仕深;;座标网格法实验的计算机数据处理[J];太原科技大学学报;1993年01期
6 张立华,朱庆,陈跃;网格法追踪等深线的多义性问题研究[J];测绘科学;2005年04期
7 赵安;3种子像元分解分类的对比研究——用NOAA影像进行鄱阳湖洪水监测为例[J];遥感技术与应用;2005年02期
8 张国富;用网格法计算任意平面形状超、亚音速机翼阻力[J];南京航空航天大学学报;1979年01期
9 廖森,廖正福,陈超球;化学试验中多目标优化问题计算方法的初探[J];广西师院学报(自然科学版);1997年01期
10 肖峰;李雪燕;邵战林;;论城镇土地定级网格法合理运用[J];新疆农业大学学报;2006年04期
11 王宏;一类依赖时间的非线性偏微分方程的非协调元方法的稳定性和收敛性[J];山东大学学报(理学版);1985年04期
12 张立华;刘家伟;张晓明;王建文;;规则网格法自动勾绘等深线[J];东北测绘;1997年04期
13 阮青;;利用参量法求地形点上的坡度[J];科技信息;2010年26期
14 谢惠民,戴福隆,王欢,杨海强,刘宁,高聚宁,庞世瑾;纳米级变形的扫描隧道显微镜测量研究[J];力学学报;1997年03期
15 夏嘉生;数据百分比网格法在宁芜及溧水地区地质找矿中的应用[J];江苏地质;1998年01期
16 王海庆;;基于GIS和RS的矿山地质环境评价方法比选[J];国土资源遥感;2010年03期
17 刘寒松;焦予秦;;结合响应面法与网格法的低速翼型优化设计[J];航空工程进展;2011年02期
18 谢敦礼;盛淑云;姜仁豪;;网格法计算底缝拱坝的应力[J];浙江大学学报(理学版);1978年01期
19 姜子文,余德浩;Sobolev方程的矩形网格混合体积元方法[J];山东师范大学学报(自然科学版);2005年01期
20 党兰学;侯彦娥;;GRIB数据中气象要素等值线的绘制[J];河南大学学报(自然科学版);2010年06期
中国重要会议论文全文数据库 前10条
1 王言金;冯其京;张树道;郝鹏程;;基于动网格的欧拉守恒方程的间断有限元方法研究:1D[A];第十四届全国激波与激波管学术会议论文集(下册)[C];2010年
2 何洋洋;高静怀;;任意高阶间断有限元法在弹性波模拟中的应用研究[A];中国地球物理学会第二十四届年会论文集[C];2008年
3 包芸;姚冬静;;透水网格法模拟计算伶仃洋抛石促淤[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
4 黄秀声;余伟;陈银莉;;数字图像处理在网格法中的应用[A];全国炼钢连铸过程自动化技术交流会论文集[C];2006年
5 于颖;李永生;於孝春;胡毅;陈立苏;;网格法在波纹管结构优化设计中的应用[A];第九届全国膨胀节学术会议膨胀节实用技术精选集[C];2006年
6 宋少卿;陈亮;;摩擦阻力作用下货船经济航速的计算[A];2009航海技术理论研究论文集[C];2009年
7 王景焘;刘凯欣;;基于时空守恒元解元方法的流体弹塑性问题计算[A];中国力学学会学术大会'2005论文摘要集(上)[C];2005年
8 张蕾;李筠;;小型灯反射镜曲面定点法线方向求解的新方法[A];2008中国仪器仪表与测控技术进展大会论文集(Ⅲ)[C];2008年
9 买歌菲热提;菊池正纪;买买提明艾尼;庄茁;;基于有限元重合网格法的三维疲劳裂纹扩展模拟与验证[A];北京力学会第18届学术年会论文集[C];2012年
10 隋常玲;蔡德利;翟瑞常;汪春;;黑龙江精准农业示范区土壤养分空间变异研究[A];农业机械化与新农村建设——中国农业机械学会2006年学术年会论文集(下册)[C];2006年
中国重要报纸全文数据库 前4条
1 吴山青;知难而上为推广[N];光明日报;2001年
2 采写/刘观涛;挑选教辅的“实用指南”[N];中国图书商报;2004年
3 王碧珺;“被”平均的无奈?[N];21世纪经济报道;2010年
4 鲁文;山东青岛:胶州湾禁建制革等工业项目[N];消费日报;2010年
中国博士学位论文全文数据库 前10条
1 黄学海;四阶椭圆型方程若干有限元新方法和高效求解算法[D];上海交通大学;2010年
2 吴迪;自适应间断有限元方法求解三维欧拉方程的研究[D];中国工程物理研究院;2010年
3 陈韵骋;线性弹性力学问题的间断有限元方法研究[D];上海交通大学;2011年
4 李灿华;平均间断有限元的强超收敛性及其应用[D];湖南师范大学;2011年
5 廉西猛;海水入侵问题的间断有限元方法[D];山东大学;2010年
6 秦树杰;对流占优扩散问题的某些新型间断有限元方法[D];南开大学;2001年
7 宋伦继;应用于偏微分方程的增量未知元和类小波增量未知元方法[D];兰州大学;2009年
8 李东方;几类常及抛物型微分方程的数值算法研究[D];华中科技大学;2011年
9 赵国忠;变分迭代方法及Lagrange坐标系下的间断有限元方法研究[D];中国工程物理研究院;2011年
10 杨继明;渗流驱动问题间断有限元高效数值方法研究[D];湘潭大学;2007年
中国硕士学位论文全文数据库 前10条
1 张海湘;间断有限元求解奇异摄动反应扩散方程[D];湖南师范大学;2010年
2 房金伟;间断有限元方法中限制器的研究[D];哈尔滨工业大学;2011年
3 刘彦忠;流动与传热问题的间断有限元法数值模拟[D];上海交通大学;2012年
4 杨芳;局部间断有限元求解带有抛物边界层的奇异摄动方程[D];湖南师范大学;2011年
5 刘妍;对流扩散方程的特征间断有限元方法[D];南京大学;2012年
6 饶向阳;胡适方法:一种现代中国元方法的建构与实践[D];华中师范大学;2003年
7 李少芳;两类方程的有限体积元方法[D];山东师范大学;2012年
8 丛美芹;一类四阶抛物型积分—微分方程的混合有限体积元方法[D];山东师范大学;2012年
9 王星;半线性抛物问题基于应力佳点的一类二次和三次有限体积元方法[D];天津师范大学;2011年
10 贾保敏;两类非线性发展方程(组)的有限体积元方法[D];山东师范大学;2010年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978