一致性高阶无单元伽辽金法及裂纹扩展分析

来源 :大连理工大学 | 被引量 : 3次 | 上传用户:zhuyong006
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
无单元伽辽金法(Element-Free Galerkin method,EFG)是近二十余年逐步发展起来的一种无网格方法,由于其具有不依赖于网格单元建立近似函数、易于构造高阶近似以及形函数高度光滑等优点,因此在裂纹扩展、自适应分析、板壳计算以及大变形模拟等方面展现出显著的优势,极具发展潜力。然而,由于无单元伽辽金法的形函数为非多项式的有理函数,弱形式的数值积分难以准确计算,需要采用较多的数值积分点才能保证计算的稳定性,不仅导致计算效率低,而且积分精度也不够,不能精确通过对于保证收敛性具有重要意义的分片试验,这在很大程度上限制了其在工程上的广泛应用。如何理性地减少积分点数目,同时保证计算精度,从而显著提高无单元法的计算效率是具有重要研究意义的科学问题,这正是本文主要研究的问题。本文工作的前期基础是段庆林等于2012年提出的无网格法的“导数一致性框架”。采用由该框架确定的修正导数计算刚度阵,可达到减少积分点数目和提高计算效率的效果。然而,该方法将修正导数引入到刚度阵缺乏理论基础。本文基于Hu-Washizu三变量变分原理重新推导出该“导数一致性框架”,而且修正导数自然地出现在弱形式中,从而为该方法奠定了数学基础,并将其称之为一致性无单元伽辽金法(Consistent Element-Free Galerkinmethod,CEFG)。对于平面问题,本文分别建立了二阶和三阶CEFG方法,大幅度减少了 EFG所需的积分点数目,同时改善了计算精度和收敛性,因而显著提高了 EFG方法的计算效率。同时,本文还将该方法拓展到三维,建立了二阶一致四点积分方法(Quadratically Consistent 4-point integration method,QC4),显著改善了三维无单元法的计算精度、收敛性和计算效率。本文还发展了一致性无单元伽辽金法的一点积分方法。该方法基于所建立的导数修正一致性框架,通过引入泰勒展开技术进一步减少积分点数目,在每个积分子域上仅使用一个积分点,可精确再生线性应变场且无沙漏模式产生。所建立的一点积分方法能精确通过线性和二次分片试验,因而称之为二阶一致1点积分方法(QuadraticallyConsistent 1-point integration method,QC1)。与其它一点积分方法相比,所建立的QC1方法在计算精度、收敛性、稳定性以及计算效率等方面均展现出显著优势。裂纹扩展的数值模拟是无单元伽辽金法的重要应用领域之一。本文将所建立的一致性高阶无单元伽辽金法由连续体问题拓展到裂纹问题(即非连续体问题)。采用虚拟节点法描述裂纹处位移的强间断,提出了虚拟节点的引入算法和断裂“单元”的数值积分方法。相较于标准的高阶无网格法和低阶一致性无网格法,一致性高阶方法显著改善了应力强度因子的计算精度,同时能够准确预测裂纹扩展路径。进一步地,本文充分利用了无单元法易于局部加密节点的优点,随裂纹扩展在裂尖附近自适应地增加计算节点,实现了裂纹扩展过程的自适应模拟,显著减少了节点数目,缩小了计算规模。此外,本文还对无单元法精确施加本质边界条件开展了研究工作,分别考察了耦合有限元和耦合权函数的本质边界条件施加方法,并给出了改进措施。
其他文献
数控机床是国家装备制造业战略性产业,电主轴是数控机床的关键功能部件。电主轴控制系统性能影响数控机床的加工质量,改善其控制性能具有现实的意义。电主轴输出转矩脉动强弱和逆变器输出电流谐波大小对其控制性能有重要的影响。当定子电阻的辨识结果有偏差时,定子磁链轨迹就会发生畸变,使电主轴输出转矩脉动;电主轴采用逆变器供电,当电源中含有电流谐波时,会形成谐波效应,产生电磁谐波损耗。本文在深入研究生物地理优化(B
覆层附着性是评定覆层零件性能及其使役可靠性的核心问题,膜基结合力测量是覆层附着性评价的主要方法。针对划痕试验中软膜硬基体系的膜基结合力难以量化,划痕形貌和膜基结合力关系不完全清楚的难题,本文采用闭合场非平衡高功率调制脉冲磁控溅射(Modulated Pulsed Power Magnetron Sputtering,MPPMS)技术,在 Si(100)基体上沉积厚度0.2~1.7 μm的Cu薄膜,
近断层地震动具有显著异于远场地震动的运动特征,如大幅值、长周期速度脉冲和永久地面位移,这些独特的运动特征造成了严重的结构破坏、人员伤亡和经济损失。所以,近断层地震动的工程特性及其结构效应研究成为近二十年来地震工程和工程力学领域的重要课题。随着我国基础设施建设的不断推进以及城镇化的持续发展,设计建造了大量的高层/超高层建筑、大跨度桥梁和高速铁(公)路等工程结构。深入研究近断层地震动作用下建筑结构的地
三维流形拓扑理论主要研究三维流形的拓扑性质和结构。通常,利用三维流形中的一些曲面(如Heegaard曲面、不可压缩曲面、正则曲面等)来研究三维流形的拓扑性质和结构是行之有效的重要方法。柄体是3-流形的基本块。每个闭三维流形均可以表示成两个等亏格柄体沿其边界的并,这就是三维流形的Heegaard分解。Heegaard分解理论是研究三维流形的重要的方法。柄体中存在一组真嵌入的互不相交的圆片,使得沿这组
锥均衡约束的多目标优化是指均衡约束中含有闭凸锥定义的参数变分不等式或者广义方程的多目标优化问题.这类问题是均衡约束的数学规划(MPECs)和均衡约束的均衡问题(EPECs)的推广,它在经济、工程、能源等诸多领域有着广泛的应用.本文主要研究了锥均衡约束的多目标优化的最优性理论,包括参数变分不等式约束的多目标优化问题的最优性条件;二阶锥广义方程约束的多目标优化问题的最优性条件;参数变分不等式约束的随机
纳米器件具有灵敏度高、选择性强、功耗低等特点,在电子信息、生物医疗、环境能源等领域具有广阔应用前景。纳米尺度结构是纳米器件的基础,是实现纳米器件高性能的关键。目前,纳米尺度结构的制造过程普遍需要借助昂贵的聚焦离子束、电子束曝光、纳米压印等设备和多层套刻、结构转移等复杂工艺步骤,其工艺流程复杂、效率低、成本高,阻碍了纳米器件的广泛应用。本文基于电流体动力理论和多物理场复合作用,研究了电场-流体复合约
奥氏体不锈钢经等离子体基低能离子注渗工艺处理后可以形成氮超高过饱和(10-35 at.%)渗氮层。该渗氮层是非均匀的亚稳单相,耐磨抗蚀性能远优于奥氏体不锈钢,在核工业中具有广阔的应用前景。然而,与渗氮层优异性能相对的是,渗氮层晶体结构存在争议,缺乏原子尺度结构模型,制约了对渗氮层相稳定性和结构稳定性的认识,从而限制了这项技术在核工业关键零部件中的应用。本论文从奥氏体不锈钢过饱和渗氮层显微组织结构入
在过去的20多年里,非标准增长条件的偏微分方程和变分问题,以及相应的变指数空间理论是非常有吸引力的研究课题。这些研究都涉及了非弹性力学、电流变流学、图像处理等实际问题。本文的研究也密切此领域,具体地,本文研究了p(x)-Laplace方程Robin边界条件下的特征值问题,包括特征值的存在性问题和稳定性问题。另外,为了研究p(x)-Laplace方程Dirichlet边界条件下第一特征值的性质,证明
激波是可压缩流动中的一种重要间断现象,越过激波流动参数发生突跃,而且随之还有机械能的损失,是个不可逆的过程。在含有激波的流动模拟中,对于激波的处理是极为重要的。激波装配方法是一种古老的激波求解方法,其理论基础较为清晰合理。在计算流体力学(Computational Fluid Dynamics,CFD)发展的初期就已经被用于来处理含有激波的可压缩流动。在处理含有简单激波的流动时,激波装配有着较为突
典型李超代数可分为基本典型李超代数和奇异李超代数,其中基本典型李超代数包含:单李代数,A(m,n),A(n,n),B(m,n),C(n),D(m,n),G(3),F(4)和 DD(2,1,α);奇异李超代数包含两类:P(n)和Q(n).用К表示除A(n,n)和单李代数之外的基本典型李超代数.单李代数的每个局部导子或者2-局部导子都是导子,本文将李代数的局部导子和2-局部导子的定义推广到李超代数上,