收藏本站
《苏州大学》 2012年
收藏 | 手机打开
二维码
手机客户端打开本文

无单元Galerkin方法的研究及应用

曹阳  
【摘要】:无网格方法是一种新兴的偏微分方程数值求解方法,受到了国际科学与工程界学者的高度重视,具有重要的研究价值和应用前景。在已提出的众多无网格方法中,无单元Galerkin(Element-free Galerkin, EFG)方法是被研究最多、应用最广的方法。本文对EFG方法从形函数的构造、本质边界条件的施加、与其他数值方法相结合以及离散线性方程组的求解四个方面进行了系统的研究,提出了一些新方法并应用于二维弹性力学、Poisson方程等问题的计算。 无单元Galerkin方法是一种基于MLS近似的全局弱式无网格方法。MLS近似方法的主要缺点之一是在每个计算点处都需要计算力矩矩阵的逆,这使得计算量大大增加。针对这一缺点,一些学者提出了用加权正交基函数代替一般基函数的方法并被称为改进移动最小二乘近似(IMLS)方法。研究表明IMLS方法比MLS方法有更好的计算效率和计算精度,但不能克服可能产生的病态甚至奇异力矩矩阵。本文基于矩阵理论详细讨论了IMLS近似方法,指出了产生奇异性的原因,并且基于单项式基函数(可以推广到任意基函数),通过矩阵三角化过程,提出了一种新的构造形函数的方法。研究表明使用此算法不仅简单,更能有效地克服力矩矩阵的奇异性。 无网格径向基点插值(RPIM)方法也从本质上来讲也是一种无单元Galerkin方法,它与EFG方法的不同之处在于形函数的构造。RPIM方法中采用径向基函数点插值方法构造形函数,所得到的形函数满足Kronecker delta函数性质,但该方法比EFG方法的计算量大,这主要是由于每个计算点所得到的力矩矩阵的维数比MLS近似方法得到的力矩矩阵的维数大得多。目前在RPIM方法中降低计算量方面的研究很少,本文提出了一种加权节点径向基点插值无网格方法(WN-RPIM).这种方法通过求节点处所产生的力矩矩阵的逆来得到节点系数,再通过节点系数加权得到形函数,所得到的形函数满足Kronecker delta函数性质。与RPIM方法相比该方法大大减少了计算量。几个数值算例验证了该方法的有效性。 MLS近似方法的另一个主要缺点是构造出的形函数一般不满足Kronecker delta函数性质,因而在无单元Galerkin方法中本质边界条件不易施加。本文从求解区域、变分原理、形函数的构造、离散系统方程四个方面出发对已提出的十多种常用的施加本质边界条件的方法进行了总结和分析,并提出了一种EFG-RPIM耦合方法。EFG-RPIM耦合方法可以直接施加本质边界条件且无需如FEM-EFG耦合方法中的过渡区域,因而具有一定的优越性。 无单元Galerkin方法是一种基于全局弱式的无网格方法,由于在整个求解域上都需要积分使得该方法的计算量较大。基于Galerkin法和最小二乘法相结合的Galerkin最小二乘无网格法是传统Galerkin法和最小二乘法的改进,它具有计算量小、精度高等优点,但由于使用了MLS近似方法构造形函数使得本质边界条件需做特殊处理。而径向基函数点插值法具有形式简单、与空间维数无关、各向同性等优点,且由该插值方法构造的形函数具有Kronecker delta函数性质,因而本质边界条件易于施加。本文结合各自的优点,提出了基于径向基函数的Galerkin最小二乘无网格法。这种方法具有本质边界条件易于施加、精度高、稳定性好等优点。 研究表明EFG方法中用增广Lagrange乘子法施加本质边界条件具有较高的精度,但所得到的线性方程组通常是不定的,常被称为鞍点问题。预处理GMRES方法是有效求解这类线性方程组的方法之一,它需要预先构造一个预条件子,而一个好的预条件子既能减少总的迭代次数,又能减少每步迭代的工作量。本文构造了一类新的分裂预条件子和一种松弛HSS预条件子结合预处理GMRES方法对这类线性方程组进行求解,明显提高了收敛速度。此外,还对两类块三角预条件子所导出的预处理矩阵特征值的界做了更精确估计。数值结果验证了本文所得到的特征值界的更优估计以及两类新的预条件子的有效性。
【学位授予单位】:苏州大学
【学位级别】:博士
【学位授予年份】:2012
【分类号】:O241.82

手机知网App
【参考文献】
中国期刊全文数据库 前10条
1 陈文;;奇异边界法:一个新的、简单、无网格、边界配点数值方法[J];固体力学学报;2009年06期
2 刘欣,朱德懋,陆明万,张雄;基于流形覆盖思想的无网格方法的研究[J];计算力学学报;2001年01期
3 张建辉,邓安福;基于正交基函数的薄板弯曲无单元法MLS导函数及其应用[J];计算力学学报;2003年03期
4 肖毅华;胡德安;韩旭;;弹性静力问题的无网格弱-强形式结合法[J];计算力学学报;2010年05期
5 孟广伟;周立明;李锋;彭惠芬;沙丽荣;;局部正交无网格伽辽金法的研究及其在含多裂纹多孔结构中的应用[J];计算力学学报;2010年06期
6 赵敏;陈文;;基于径向基函数的加权最小二乘无网格法[J];计算力学学报;2011年01期
7 丁睿;魏学润;;基于Uzawa算法的弹塑性扭转问题的局部微分求积法[J];计算力学学报;2011年04期
8 蒋美群;曹阳;;广义鞍点问题的块三角预条件子[J];计算数学;2010年01期
9 龙述尧,许敬晓;弹性力学问题的局部边界积分方程方法[J];力学学报;2000年05期
10 程玉民,陈美娟;弹性力学的一种边界无单元法[J];力学学报;2003年02期
【共引文献】
中国期刊全文数据库 前10条
1 唐习龙;林可勃;贾妮;李忠;;结构计算软件中的模块化设计[J];四川建筑科学研究;2005年06期
2 钱坤;尹新生;朱珊;佴磊;;三维地基沉降计算公式的改进与比较[J];四川建筑科学研究;2009年03期
3 熊学玉;王绍辉;;现浇空心板的正交异性板理论分析方法的研究[J];四川建筑科学研究;2010年05期
4 郝婷玥;陈贵清;;流固耦合作用下埋地管道的振动分析[J];四川建筑科学研究;2011年03期
5 胡卫兵;韩广森;于海平;;古建筑榫卯节点刚度分析[J];四川建筑科学研究;2011年06期
6 ;Incompatible numerical manifold method for fracture problems[J];Acta Mechanica Sinica;2010年02期
7 宁作君;桩基侧向位移的探讨与分析[J];鞍山钢铁学院学报;2000年03期
8 楼芬;邓建;;无网格法及其在岩石力学与工程中的应用[J];地下空间与工程学报;2007年06期
9 周小义;邓安福;丁灿;;流形方法及其研究进展[J];地下空间与工程学报;2008年05期
10 贺广零;李倩妹;洪芳;彭双艳;;基于反演法的山区变电站高填方地基沉降分析[J];地下空间与工程学报;2008年05期
中国重要会议论文全文数据库 前10条
1 周晓敏;纪洪广;邓斌;;高水压基岩井壁厚度设计方法[A];2011年中国矿业科技大会论文集[C];2011年
2 潘小飞;张雄;陆明万;;配点型无网格方法的研究[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
3 高凌天;刘凯欣;刘颖;;无网格方法在冲击动力学中的应用[A];计算爆炸力学进展[C];2006年
4 朱训国;杨庆;栾茂田;;隧道围岩中全长注浆岩石锚杆的应力分布分析[A];第九届全国岩土力学数值分析与解析方法讨论会论文集[C];2007年
5 顾士坦;赵同彬;;SMW工法型钢-水泥土组合梁抗弯性能分析[A];第九届全国岩土力学数值分析与解析方法讨论会论文集[C];2007年
6 杨仲;张博明;孙新杨;;内压载荷下夹层复合材料圆筒开口分析[A];第十五届全国复合材料学术会议论文集(下册)[C];2008年
7 姚林泉;王璠;黄娟;;改进广义移动最小二乘近似无网格法及其在平面压电问题中的应用[A];第三届全国压电和声波理论及器件技术研讨会论文集[C];2008年
8 陈文涛;高存法;杨权权;;含圆孔功能梯度板在对称载荷作用下的应力分析[A];第三届全国压电和声波理论及器件技术研讨会论文集[C];2008年
9 郭瑞平;刘曙云;;热力耦合作用下含表面裂纹圆柱壳体毁伤分析[A];第七届全国工程结构安全防护学术会议论文集[C];2009年
10 任业军;李磊;;SPH方法在高速碰撞冲击波传播中的应用研究[A];第九届全国冲击动力学学术会议论文集(上册)[C];2009年
中国博士学位论文全文数据库 前10条
1 郭庆勇;非保守系统Fourier卷积型拟变分原理研究及应用[D];哈尔滨工程大学;2010年
2 张选利;柔性注压锚杆锚固特性研究[D];辽宁工程技术大学;2010年
3 张祖军;无网格法及其在注塑充填数值分析中的应用[D];华南理工大学;2010年
4 刘海锋;网格结构精细化有限元分析方法研究[D];浙江大学;2010年
5 杨溢;爆破荷载对蠕动边坡的累积效应及稳定性影响研究[D];昆明理工大学;2010年
6 卢广林;立方氮化硼仿生耐磨复合材料的制备及其研究[D];吉林大学;2011年
7 周立明;不确定性结构的无网格伽辽金法研究[D];吉林大学;2011年
8 张俊兵;基于波谱单元法的结构动力分析[D];华中科技大学;2011年
9 谭飞;杂交边界点法理论及其在薄板问题中的应用[D];华中科技大学;2011年
10 李坤;基于无网格方法的声学问题数值模拟研究[D];华中科技大学;2011年
中国硕士学位论文全文数据库 前10条
1 闫丽梅;高层建筑结构自由振动分析的ODE方法[D];河南理工大学;2010年
2 李朝;箱形桥梁结构的半解析分析[D];河南理工大学;2010年
3 张继华;急倾斜中厚煤层软底综采采场矿压规律及其控制研究[D];山东科技大学;2010年
4 银国伟;谐波摩擦传动柔轮径向力与径向变形关系的分析及实验[D];郑州大学;2010年
5 王士革;填方路基下长短桩复合地基沉降特性研究[D];郑州大学;2010年
6 董权;钢索试验装置改进研究[D];哈尔滨工程大学;2010年
7 陈娟;水下爆炸冲击载荷的SPH算法研究[D];哈尔滨工程大学;2010年
8 孙巧英;光纤陀螺捷联系统振动特性研究[D];哈尔滨工程大学;2010年
9 修孝廷;双层周期加筋板声学特性数值方法研究[D];大连理工大学;2010年
10 张宛君;齿轮传动接触问题的数值模拟研究[D];辽宁工程技术大学;2009年
【二级参考文献】
中国期刊全文数据库 前10条
1 ;AN ASSESSMENT OF THE MESHLESS WEIGHTED LEAST-SQUARE METHOD[J];Acta Mechanica Solida Sinica;2004年03期
2 刘欣,朱德懋,陆明万,张雄;基于流形覆盖思想的无网格方法的研究[J];计算力学学报;2001年01期
3 张建辉,邓安福;基于正交基函数的薄板弯曲无单元法MLS导函数及其应用[J];计算力学学报;2003年03期
4 张雄,宋康祖,陆明万;无网格法研究进展及其应用[J];计算力学学报;2003年06期
5 庞作会,葛修润,郑宏,王水林;一种新的数值方法——无网格伽辽金法(EFGM)[J];计算力学学报;1999年03期
6 贝新源,岳宗五;三维SPH程序及其在斜高速碰撞问题的应用[J];计算物理;1997年02期
7 宋康祖,陆明万,张雄;固体力学中的无网格方法[J];力学进展;2000年01期
8 龙述尧;弹性力学问题的局部Petrov-Galerkin方法[J];力学学报;2001年04期
9 程玉民,陈美娟;弹性力学的一种边界无单元法[J];力学学报;2003年02期
10 张雄,胡炜,潘小飞,陆明万;加权最小二乘无网格法[J];力学学报;2003年04期
【相似文献】
中国期刊全文数据库 前10条
1 张伟星,庞辉;无单元法计算钢筋混凝土筏板[J];计算力学学报;2000年03期
2 陈莘莘,李庆华;用无单元法求解稳态热传导问题[J];株洲工学院学报;2003年05期
3 秦雅菲,张伟星;无单元法分析薄板自由振动问题[J];力学与实践;2003年05期
4 周瑞忠,周小平,缪圆冰;小波基无单元法及其工程应用[J];工程力学;2003年06期
5 邹振祝,陈建;无单元法在孔洞应力集中问题中的应用[J];石家庄铁道学院学报;2000年02期
6 周小平;对进一步发展无单元法的几点设想[J];福州大学学报(自然科学版);2001年03期
7 唐少武,冯振兴;关于无单元法的若干注记[J];力学进展;2003年04期
8 白泽刚,杨元明;一种新型核函数下的无单元法及应用[J];西安建筑科技大学学报(自然科学版);2001年01期
9 张振南,卢爱红,茅献彪;无单元法在薄板弯曲问题中的应用[J];河北建筑科技学院学报(自然科学版);2001年03期
10 胡云进,周维垣,寇晓东;三维无单元法及其应用[J];岩石力学与工程学报;2004年07期
中国重要会议论文全文数据库 前10条
1 樊志斌;贾晓峰;;一种高效率无单元偏移成像算法[A];中国地球物理学会第二十七届年会论文集[C];2011年
2 刘素贞;商植桐;陈海燕;杨庆新;杨晓光;;无单元法在二维电磁场数值计算中的初步研究[A];全国电工理论与新技术学术年会(CTEE'2001)论文集[C];2001年
3 介玉新;李广信;高波;葛锦宏;;无单元法在应力应变分析和渗流分析中的应用[A];中国土木工程学会第九届土力学及岩土工程学术会议论文集(下册)[C];2003年
4 庞辉;张伟星;;钢筋混凝土双参数地基板分析的无单元法[A];“力学2000”学术大会论文集[C];2000年
5 王建国;M.R.Karim;;波浪作用下海床瞬态液化的无单元法分析[A];中国土木工程学会第九届土力学及岩土工程学术会议论文集(下册)[C];2003年
6 王燕昌;张虹;;无单元法在三维弹塑性问题中的应用[A];力学与西部开发会议论文集[C];2001年
7 陈虹;林建华;;用无单元法求解水流运动方程[A];第十四届全国水动力学研讨会文集[C];2000年
8 沈明;胡云进;毛根海;;应用无单元法模拟岩体水力劈裂[A];第15届全国结构工程学术会议论文集(第Ⅱ册)[C];2006年
9 刘冬欢;郑小平;刘应华;;间断Galerkin有限元方法在由接触热阻引起的热力耦合问题上的应用研究[A];北京力学会第十六届学术年会论文集[C];2010年
10 张建辉;马丽红;;弹性地基板无单元法权函数及其参数研究[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
中国博士学位论文全文数据库 前10条
1 曹阳;无单元Galerkin方法的研究及应用[D];苏州大学;2012年
2 张琰;高土石坝张拉裂缝开展机理研究与数值模拟[D];清华大学;2009年
3 曹中清;无单元法及其面向对象程序实施[D];西南交通大学;2006年
4 王立鹏;无单元法若干技术研究及其在电磁场计算中的应用[D];沈阳工业大学;2010年
5 王飞;椭圆变分不等式的间断Galerkin方法[D];浙江大学;2010年
6 周晓杰;堤防的渗透变形及其发展的研究[D];清华大学;2006年
7 王志亮;软基路堤沉降预测和计算[D];河海大学;2004年
8 谷照升;水库湖泊水质分析、模拟与预测的综合数学方法及其应用[D];吉林大学;2006年
9 刘敦平;双向增强体复合地基承载特性及其数值分析[D];湖南大学;2008年
10 王晓峰;多连通域上的解析函数空间及其算子[D];四川大学;2003年
中国硕士学位论文全文数据库 前10条
1 吴小伟;无单元法在弹性地基圆板弯曲中的应用研究[D];重庆大学;2012年
2 阎树文;无单元法用于高分辨率地震模拟与成像[D];中国科学技术大学;2010年
3 王国庆;岩体水力劈裂试验及裂纹扩展的无单元法计算[D];河海大学;2004年
4 刘婷;用无网格Galerkin法求解地下水流拟三维问题[D];辽宁师范大学;2010年
5 缪圆冰;小波无单元方法及其应用[D];福州大学;2002年
6 刘学文;基于插值型形函数构造理论的无单元法研究[D];宁夏大学;2003年
7 钟凤远;水库水体中BOD-DO迁移转化的数学模型及其应用[D];吉林大学;2005年
8 商敬铜;向前向后热传导方程的时间间断Galerkin方法[D];山东大学;2012年
9 李海龙;无网格伽辽金法精度分析与本质边界条件的处理[D];燕山大学;2010年
10 张闪闪;马氏体相变模型的不连续Galerkin方法[D];上海交通大学;2012年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026