样条有限元方法和比例边界有限元方法的若干研究

来源 :大连理工大学 | 被引量 : 0次 | 上传用户:xunmengya
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
有限元方法是数值计算的有力工具,也是处理复杂工程问题的重要手段.样条有限元方法和比例边界有限元方法都是基于有限元方法发展起来的新的数值方法.样条有限元方法发挥了样条函数满足一定分片连续性,逼近精度高的优点.比例边界有限元方法是由Wolf和Song提出来的一种新的半离散半解析方法.该方法结合了有限元方法和边界元方法,不仅继承了两者的优点,而且拥有自己的特点,更加灵活和有效.目前,样条有限元方法和比例边界有限元方法在数值计算和工程领域发展迅速,应用广泛.本文在数学理论上对两种方法以及两者之间的结合进行了研究,具体研究工作主要包括以下三部分内容.1.二阶椭圆问题的比例边界有限元的高阶完备性分析迄今为止,比例边界有限元方法在工程应用方面取得了许多丰富的研究成果,但在数学理论方面的研究工作较少.对有限元理论,单元的完备性分析是比例边界有限元理论基础中非常重要的一部分.本文在第3章严格地给出了比例边界有限元方法对于二阶问题的高阶完备性的理论证明.我们从比例边界有限元的齐次解(无体力项)和非齐次解(有体力项)的表达形式入手,引入了比例边界坐标下的环向插值算子,给出了多项式在比例边界坐标下的表达式.通过分析得到,比例边界坐标下的插值算子能否精确重构多项式的关键在于多项式表达形式中的特定整数幂次r(多项式的次数)能否存在.一方面,齐次解中的r是由特征值分解法的求解过程中的[Z]矩阵的特征值决定;另一方面,非齐次解的特解中的r由体力项的逼近精度决定.本章严格证明了,对于封闭的S单元,特定的整数幂次r总是可以在比例边界有限元的计算中得到并且和S单元的形状无关,即比例边界有限元方法具有高阶完备性.另外,在完备性分析中,我们也发现了一些比例边界有限元在求解过程中的相关理论问题并给出了必要的证明.2.三维二阶椭圆问题的样条比例边界有限元比例边界有限元方法对网格剖分适应性好,尤其是它在多边形或者多面体单元下的方程推导过程和三角形、四边形或者六面体单元没有差别.实际上,对于三维问题,多面体单元的构造是相当困难的.考虑到比例边界有限元方法中多面体单元构造简单,通用性好,本文第4章将样条单元和其相结合,分别构造了二次和三次样条比例边界有限单元SBFEM-L8和SBFEM-L12.该单元是将四边形样条单元L8或L12作为面单元,用于三维比例边界有限元方法中.由第3章的理论分析和第4章的数值实验可得,SBFEM-L8具有二阶完备性,SBFEM-L12具有三阶完备性.除此之外,数值实验也表明,SBFEM-L8和SBFEM-L12既保留了比例边界有限元对网格适应性好的特点,又发挥了样条单元节点少、精度高且对网格畸变不敏感的优势.3.四阶椭圆问题的超收敛非协调四边形样条单元对于样条有限元方法,由于样条函数在单元内连续性灵活,对网格适应性强,本文第5章针对四阶椭圆问题,构造了一个12自由度的非协调的四边形样条单元NCQS12.该单元基于Ⅱ型三角剖分上的样条空间S31(QT),用B网方法选择了一个包含完整三次多项式的S31(QT)的子空间作为样条有限元空间.自由度选为顶点的函数值、每条边上函数的积分值以及法向导数的积分值.理论分析得到,样条单元NCQS12的插值误差为O(h2),相容误差为O(h1).特别地,若网格为平行四边形网格,相容误差可以达到O(h2),即该单元具有超收敛性.数值实验验证了我们的理论结果.此外,对于两种退化网格:二分网格和渐近规则的平行四边形网格,数值实验表明,单元NCQS12仍具有超收敛性.
其他文献
光学显微镜被广泛应用于细胞学、肿瘤学、遗传工程学及植物学等各领域的生物医学检测。传统高性能光学显微镜存在体积较大、价格昂贵且需要专业人员在实验室操作完成等不足,限制了显微检测技术面向个人用户及基层的现场检测应用。因此开发一种微米级分辨率、小畸变的便携式显微镜,满足现场实时检测的需求,具有重要的现实意义和广阔的应用前景。针对以上检测需求,光学透镜由于性能稳定、造价低廉且原理简单等优点被作为市面上商用
随着人类对世界探索的扩展与深入,所涉猎环境愈发复杂化和多样化,这对光电设备提出愈来愈多的挑战,传统理论所构建的仪器已趋于技术极限,人们进而研究开发新的技术,偏振技术便在近些年来应运而生。然而,相关技术研究不足却逐渐限制着偏振技术的发展,因此,亟需在关键技术研究方面有所进展。光在散射环境中传输时的偏振特性及其相关技术作为能够促进偏振技术发展的重要问题,近年来备受各国研究机构重视。在散射环境中,Mie
本文的研究内容分为以下两个方面.首先,探讨了有限群不变式中的超曲面问题.在研究有限群与其子群不变式的关系时,T-函子被证明是一有力工具.本文研究了T-函子对超曲面的作用.我们证明了在Steenrod代数上的不稳定代数范畴中,T-函子不会增加作用对象的嵌入维数.进而有T-函子在不稳定代数范畴内保持超曲面这一结论.从而可知如果有限群的不变式是超曲面,那么它稳定子群的不变式也是超曲面.反之则不然,本文给
风机基础结构在复杂荷载作用下的疲劳分析是海上风电场建设中的关键技术问题之一。对于海上风机三桩基础结构而言,其疲劳分析关键管节点为三平面Y型管节点,该节点为空间管节点,且在服役期承受多种疲劳荷载(如风、浪、流、水位、冲刷等),因此该节点的热点应力往往受反应叠加效应和多平面相互作用影响。各国规范和已发表文献中缺乏对复杂荷载同时作用下热点应力计算方法的系统探讨,也没有对多平面相互作用的定量结论,更没有适
随着我国经济不断加速发展,基础建设规模不断扩大,尤其在地质条件及其复杂的西南地区,其开发力度和建设规模都不断增加,许多大规模的地下洞室群开挖不可避免。岩体的开挖就会破坏其原有的力学平衡状态,受扰动岩体应力将会重新分布,从而会形成局部应力集中区域,导致围岩较大变形或者诱发一系列的地质灾害,如坍塌、岩爆等。其中,岩爆是脆性岩石地区高应力环境下开挖过程中产生的特殊破坏形式,其发生时伴随的强烈能量释放及围
地震的发生和地震的强度具有高度不确定性,现行建筑抗震设计规范采用了多地震水平和多设防目标的抗震设计思想,但这些方法仍属于基于力的弹性设计方法,难以全面考虑结构进入非线性阶段后的损伤破坏及性能变化,因此近年又提出了基于位移的抗震设计方法。基于位移的抗震设计方法通过限制结构在不同地震水平下的侧向位移,实现对地震下结构性能的控制,是下一代结构抗震设计规范的发展方向。然而,对于多层、多跨的钢筋混凝土(RC
导截流是一项非常复杂的系统工程,其成败直接影响着水利枢纽工程的施工安全、施工工期和工程造价。因此,对导截流工程的研究一直是水利枢纽建设中十分重要的课题。分期导截流工程中,后期导流需要通过前期已建成坝体泄流建筑物进行泄流。泄流建筑物内水流流态可能为明渠流、有压流或明满流过渡,使得二维模型无法对分期导流工程进行计算模拟。三维模型虽然能够重现分期导截流工程中导流河道及泄流建筑物内流场,但泄流建筑物尺度通
中国不仅是粮食生产大国,还是粮食消费大国,农业在我国国民经济中占有重要地位。然而近年来我国城镇化与工业化的快速发展,使得农业可用水量不断压缩,这对粮食生产安全与水资源优化配置带来了严峻的考验。我国农业用水量占总用水量的60%以上,农田灌溉又占农业用水量的90%左右,而我国目前基于地面站点的农业用水监测方法并不能完全实现农业灌溉用水的全面监测,也不能完全支撑农业水资源精细化管理工作的开展。遥感技术与
地震液化是一种常见的灾害现象。为研究其致灾程度和致灾机理,动力数值模拟是一种省时省力的研究手段,可以有效地重现、预测地震液化过程。然而,在动力数值模拟计算过程中,为了提高计算精度,时间步长设置过小,导致计算成本过高;为了节省计算时间,时间步长设置过大,导致计算结果发散甚至计算突然中断;为了提高计算效率,寻找一个合适的时间步长,在节省计算成本的同时提高计算精度,却难以实现。为了解决上述问题,本文开发
我国西南地区横跨中国地形三大台阶,地势落差大且河流密布,是我国水能资源最丰富的地区。近20年来,西南地区相继建成以金沙江、澜沧江、雅砻江、乌江、大渡河、红水河等干流梯级水电站为代表的巨型梯级水电站群。这些梯级水电站群呈现出大规模、高水头、大容量、巨型机组、面临复杂调度需求,面临灵敏水力耦合等复杂特征。这些特征进一步加剧了西南地区水电优化调度问题约束建模与求解难度。为此,本文以西南地区水电优化调度突