二维结构化网格上的数值流形法

来源 :中国科学院大学 | 被引量 : 0次 | 上传用户:csmale
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
数值流形法利用两套覆盖实现了插值网格和积分背景网格的分离,计算结果对于网格的依赖性相对有限元减弱了很多。数值流形法的数学网格不需要适应求解域的边界和各种不连续面,规则的数学网格上可以取得高精度的插值,因此,本文中我们选择使用规则的结构化网格来建立数值流形法的数学覆盖,研究二维结构化网格上的数值流形法。本文的主要工作包括以下几个方面:  1)建立了Hermite形式的数值流形法空间,用于求解四阶偏微分方程控制的问题,例如结构中的Kirchhoff薄板问题。有限元中的Kirchhoff薄板问题的Zienkiewicz单元的形函数满足我们对于Hermite形式插值的要求,但这种单元对于网格构型有非常严格的要求,在有限元的背景下这种单元无法求解形状复杂的问题。我们建立满足要求的规则数学网格,在数值流形法的背景下求解了复杂形状的薄板受力问题。我们这里在数值流形法中实现Zienkiewicz单元,可以作为用数值流形法挽救那些对网格构型有严格要求的非协调元的范例。  2)建立了两种数值流形法中结构化网格加密的方法:多重覆盖法和加密物理片方法。多重覆盖数值流形法是基于在局部区域建立的多重局部覆盖实现的,可以对任意区域加密,具有非常强的通用性,唯一的不足是这种方法的刚度矩阵是奇异的。通过在物理片上布置结构化网格,并在布置的网格上建立局部插值,作为物理片上的局部近似函数,同样可以实现结构化网格加密,我们称这种方法为加密物理片的方法。文中大量算例显示,在二维情况下,如果没有奇异物理片,加密物理片方法的刚度矩阵在施加本质边界条件之前亏秩数都是3,这表明,本质边界条件正确施加后可以得到正定的刚度矩阵。  3)数值流形法利用两套覆盖包含了对于裂纹面等不连续面两侧的位移不连续的描述,在形成不连续面上位移的近似时,不需要做任何特殊处理,就可以表征位移的不连续。因此,数值流形法非常适合求解断裂力学问题。Williams级数是线弹性裂纹尖端渐进场的解析解,用Williams级数的各项来作为裂纹尖端附近区域场变量近似中的基函数,可以精确的表达裂纹尖端的奇异性。将线弹性断裂力学问题的求解域划分成裂纹尖端区域和外部区域,分别用Williams级数和数值流形法来近似不同区域,然后用拉格朗日乘子法、罚函数法、势能区和余能区分区变分三种方法建立了三种求解线弹性断裂力学问题的区域分解方法。由于后两种方法各有缺点,文中重点研究了基于拉格朗日乘子法构造的区域分解方法,研究了拉格朗日乘子的插值函数,并通过数值实验的方式分析了几种因素对于计算精度的影响。数值算例显示这种方法对于各种复杂裂纹问题有良好的收敛性。
其他文献
随着西部大开发的不断深入和相关地区经济社会的发展,在青藏高原修建高等级公路已十分必要,目前青藏高速公路非多年冻土区已进入了实施阶段。高等级公路的路面宽度一般大于20m,
该文从分析中国建设项目投资失控的原因入手,结合中国建设监理行业发展现状,针对项目实施阶段,即项目招投标阶段和施工阶段,从投资控制经济措施和组织管理措施两个不同方面展
针对工程设计中普遍反映的现行规范GBJ10-89对在正常室内环境条件下采用高强钢材作预应力筋的结构裂缝控制标准过严的问题和以裂缝为核心的简化设计过程、优化配筋设计的实际
该文对单桩的动、静刚度和承载力的机理作了一定的探讨,分析了各种因素对刚度与承载力的影响.从分析的结果来看,刚度与承载力之间并无直接的联系.
该文利用日本大林组株式会社改进的FLUSH程序,通过对40余栋具有不同刚度的上部结构分别在基础固定以及与地基、基础共同工作时的地震反应分析比较,研究了上部结构刚度变化对
该论文重点对在垂直荷载作用下的框架一筒体结构楼板的有限元计算结果进行了分析和研究,发现A区和C区板具有单向板的受力特性,B区板具有双向板的受力特性,提出支承板梁的概念
在实际工程中,连梁设计中往往由于抗剪截面不足而引起超筋的现象,并且有的设备管道需要穿过连梁,由此引出了对连梁进行开缝或者开洞处理的思路,即双连梁构件。从而我们需要对
该文概述了电凝聚法的特点,在水处理中的应用和发展状况,并针对该法的主要缺点--能耗和铁耗大,以电镀含铬废水为处理对象进行了一系列的实验研究.
该文在查阅了大量的国内外有关文献、对现有的各种水环境质量评价方法进行了分析、归纳和总结后发现:直接采用现有的水环境质量评价方法评价居住区的水环境质量都有这 样或那