基于扩展有限元的应力强度因子的位移外推法

针对平面裂纹问题,阐述了扩展有限元法的单元位移模式、推导了扩展有限元法的控制方程、介绍了特殊单元的数值积分技术.基于最小二乘法,建立了应力强度因子位移外推法的计算公式.利用MATLAB编写计算程序,对平面裂纹问题用扩展有限元法进行了计算.基于扩展有限元法的计算结果,分别利用位移外推法和相互作用积分法,对平面裂纹的应力强度因子进行了计算.计算结果表明,位移外推法比相互作用积分法能更方便和准确地计算平面裂纹的应力强度因子.     

基于扩展有限元法的混凝土细观断裂破坏过程模拟

扩展有限元法（XFEM）是分析不连续力学问题（特别是断裂问题）的一种有效的数值方法。在常规的有限元位移模式中,基于单位分解的思想加入一个跳跃函数和渐进缝尖位移场来对不连续体附近的节点自由度进行局部加强,从而反映了位移的不连续性。介绍了扩展有限元的基本原理,给出了扩展有限元进行混凝土开裂及裂纹扩展的分析方法,最后采用扩展有限元法模拟了湿筛混凝土单轴拉伸作用下及WinklerL-型混凝土板的细观断裂破坏过程。分析了混凝土裂纹萌生、扩展的过程及破坏形态,数值结果与实验结果吻合良好。研究表明：扩展有限元法通过特定的位移模式,使裂纹两侧不连续位移场的表达独立于网格划分,能有效地模拟混凝土材料细观断裂破坏过程。     

基于扩展有限元的地质聚合物混凝土断裂过程研究

扩展有限元法是基于常规有限元框架分析裂纹等不连续力学问题的一种有效数值方法,在常规的有限元位移表达式中,增加了能够反映位移不连续性的跳跃函数和渐进缝尖位移场函数来对不连续结构附近的节点自由度进行局部加强。本文介绍了扩展有限元法及粘聚力模型的基本原理,给出了基于扩展有限元法的地质聚合物混凝土断裂过程分析方法。分别采用四种不同的软化曲线对I型缺口地质聚合物混凝土梁从裂纹萌生、扩展直至断裂破坏的全过程进行了模拟,并基于双K断裂准则分析了其断裂韧性。结果表明,Petersson模型与试验结果吻合较好,最后基于模拟结果进一步揭示了断裂过程区的演化过程。     

广义扩展有限元法及其在裂纹扩展分析中的应用

结合广义有限元法(GFEM)和扩展有限元法(XFEM)的特点,提出了一种新的数值方法——广义扩展有限元法(GXFEM)。阐述了广义扩展有限元法的基本原理,对相关公式进行推导,探讨数值实施中需注意的重要问题,给出利用广义扩展有限元法进行断裂分析时应力强度因子的计算方法,编写了广义扩展有限元法程序。通过算例进行了应力强度因子的计算,模拟了结构裂纹的扩展过程。算例结果表明,利用广义扩展有限元法计算裂纹扩展问题,不需要进行过密的网格划分,且网格在裂纹扩展后无需重新剖分,具有相当高的计算精度。     

模拟裂纹扩展的单位分解扩展无网格法

单位分解扩展无网格法(PUEM)是一种求解不连续问题的新型无网格方法.其基于单位分解思想,通过在传统无网格法的近似函数中加入扩展项来反映由裂纹所产生的不连续位移场.详细描述了水平集方法,PUEM不连续近似函数的构造及控制方程的离散.针对裂纹扩展问题,提出了一种十分简单的水平集更新算法;讨论了不同的节点数、高斯积分阶次以及围线积分区域对应力强度因子计算结果的影响,并给出了合理的参数;模拟了边裂纹和中心裂纹的扩展问题,并与XFEM的数值结果进行了比较.数值算例表明,本文方法具有较高的计算精度,是模拟裂纹扩展非常有效的方法,具有广阔的应用前景.     

含两类附加函数的扩展等参有限元法

基于扩展有限元的基本思想,提出一类指数型间断函数来模拟。由于裂纹或节理等非连续结构所导致的位移不连续现象,该附加函数是以到间断处的垂直距离为自变量,且随距离的增大而呈指数衰减,同时,在非连续结构末端引入能反映其奇异场特性的三角基函数。本文用弱解形式推导了扩展有限元格式,并论证了两类附加函数在单元公共边上能够保持位移连续性这一要求。最后,编制了二维4节点和三维8节点的扩展等参有限元程序,并分别给出了算例,结果表明在模拟裂纹追踪时,扩展有限元法可行且有效。     

弱不连续问题扩展有限元法的数值精度研究

主要研究了扩展有限元法(extended finite element method, XFEM)在处理弱不连续问题时不同改进函数形式对XFEM数值求解精度的影响,阐述了各种改进函数影响XFEM求解精度的关键因素,指出校正的扩展有限元法(corrected-XFEM)能够提高数值求解精度的实质在于它拓展了改进结点域,即将常规扩展有限元法(standard-XFEM)的改进结点域增加一层作为corrected-XFEM的改进结点域,文中建议延拓corrected-XFEM的改进结点域,即在corrected-XFEM的改进结点域基础上再增加一层改进结点. 利用水平集函数表征材料内部的不连续界面,推导了XFEM求解的支配方程,给出了一种改进单元的数值积分方案以及改进单元处高精度应力的求解方法. 含夹杂问题的数值计算结果表明:建议的延拓corrected-XFEM改进结点域的方法能够明显提高XFEM的数值求解精度.      

改进型扩展比例边界有限元法

结合了扩展有限元法(extended finite elementmethods,XFEM)和比例边界有限元法(scaled boundary finite elementmethods,SBFEM)的主要优点,提出了一种改进型扩展比例边界有限元法(improvedextended scaled boundary finite elementmethods,$i$XSBFEM),为断裂问题模拟提供了一条新的途径.类似XFEM,采用两个正交的水平集函数表征材料内部裂纹面,并基于水平集函数判断单元切割类型;将被裂纹切割的单元作为SBFE的子域处理,采用SBFEM求解单元刚度矩阵,从而避免了XFEM中求解不连续单元刚度矩阵需要进一步进行单元子划分的缺陷;同时,借助XFEM的主要思想,将裂纹与单元边界交点的真实位移作为单元结点的附加自由度考虑,赋予了单元结点附加自由度明确的物理意义,可以直接根据位移求解结果得出裂纹与单元边界交点的位移;对于含有裂尖的单元,选取围绕裂尖单元一圈的若干层单元作为超级单元,并将此超级单元作为SBFE的一个子域求解刚度矩阵,超级单元内部的结点位移可通过SBFE的位移模式求解得到,应力强度因子可基于裂尖处的奇异位移(应力)直接获得,无需借助其他的数值方法.最后,通过若干数值算例验证了建议的$i$XSBFEM的有效性,相比于常规XFEM,$i$XSBFEM的基于位移范数的相对误差收敛性较好;采用$i$XSBFEM通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合\较好.      

基于光滑有限元法的热-弹相场断裂研究

基于相场理论,论文构建了热力耦合作用下热-弹性材料的复杂断裂控制方程,采用光滑有限元法对该模型进行了数值离散;设计编写了相应的Matlab计算程序;开展了热-力载荷下的二维拉伸断裂模拟,将结果和有限元计算结果对比,吻合一致.同时采用该模型开展热冲击下陶瓷板的多裂纹扩展过程模拟,计算结果和实验结果吻合良好,验证了论文方法的有效性.     

近场动力学与有限元方法耦合求解热传导问题

求解含裂纹等不连续问题一直是计算力学的重点研究课题之一,以偏微分方程为基础的连续介质力学方法处理不连续问题时面临很大的困难. 近场动力学方法是一种基于积分方程的非局部理论,在处理不连续问题时有很大的优越性. 本文提出了求解含裂纹热传导问题的一种新的近场动力学与有限元法的耦合方法. 结合近场动力学方法处理不连续问题的优势以及有限元方法计算效率高的优势,将求解区域划分为两个区域,近场动力学区域和有限元区域. 包含裂纹的区域采用近场动力学方法建模,其他区域采用有限元方法建模. 本文提出的耦合方案实施简单方便,近场动力学区域与有限元区域之间不需要设置重叠区域. 耦合方法通过近场动力学粒子与其域内所有粒子(包括近场动力学粒子和有限元节点)以非局部方式连接,有限元节点与其周围的所有粒子以有限元方式相互作用. 将有限元热传导矩阵和近场动力学粒子相互作用矩阵写入同一整体热传导矩阵中,并采用Guyan缩聚法进一步减小计算量. 分别采用连续介质力学方法和近场动力学方法对一维以及二维温度场算例进行模拟,结果表明,本文的耦合方法具有较高的计算精度和计算效率. 该耦合方案可以进一步拓展到热力耦合条件下含裂纹材料和结构的裂纹扩展问题.