二阶一致无网格与有限元耦合离散研究

为准确方便地施加本质边界条件,在连续掺混法(Continuous Blending Method,CBM)的框架下,通过增加一个边中节点,发展了采用二阶基底的无网格与二阶有限元的耦合离散方法.Galerkin弱形式的数值积分采用具二阶一致性的3点积分方法(Quadratically Consistent 3-point integration method,QC3).与原本在QC3中采用的Nitsche法相比,所发展的耦合离散方法可像有限元法一样简单高效地施加本质边界条件,不向弱形式中引入额外项,也不依赖于任何人工参数.而且,数值结果还表明,QC3的计算精度也得到进一步提高.     

断裂过程的有限元模拟

讨论了材料断裂过程的有限元模拟技术。基于自适应有限元的一般原理，并针对多相材料的裂纹扩展的特点，提出了一种简化的高精度和高效率有限元网格的动态重新划分策略。裂纹被假设沿着单元之间的路径连续扩展，利用节点力释放技术生成新的裂纹自由表面，发展了一种可随裂尖连续移动的网格动态加密和释放方法。这种方法已在各种裂纹问题中得以实现与应用。     

无网格方法的研究进展与展望

目前正在发展的无网格方法采用基于点的近似,可以彻底或部分地消除网格,因此在处理不连续和大变形问题时可以完全抛开网格重构。无网格方法是目前科学和工程计算方法研究的热点,也是科学和工程计算发展的趋势。本文首先简单地阐述了无网格方法,然后详细叙述了目前提出的各种无网格方法的研究进展,最后对目前无网格方法存在的问题进行了探讨,提出了今后的研究方向。     

板壳畸变单元的无网格和有限元自动耦合方法

板壳大变形时单元的严重畸变会使计算精度降低。无网格局部Petrov-Galerkin法是一种真正的无网格方法,能够消除网格畸变,但比有限元法计算效率低。根据板壳网格畸变的局部性特点,利用过渡单元法,基于板壳网格质量,建立了板壳的网格严重畸变区域由有限元分析切换为无网格分析的自动耦合算法,实现了有限元法和无网格局部彼得罗夫．迦辽金法的耦合。应用实例表明：通过自适应耦合,既能发挥有限元法计算效率高的特点,又能发挥无网格法适合大变形分析、没有网格畸变造成计算困难的特点。     

无单元Galerkin法和边界元耦合法

无网格方法与有限元法或边界元法耦合是无网格方法处理边界条件的方法之一,在无网格方法中研究无网格方法与有限元法或边界元法耦合的研究显得非常重要.本文在无单元Galerkin法和边界元法的基础上,基于无单元Galerkin法子域和边界元法子域的界面上位移连续和面力平衡条件,提出了一种新的无单元Galerkin法和边界元法的直接耦合方法,对弹性力学问题详细推导了在整个求解域上的耦合公式.与以往的耦合法相比,这种方法简单直观,不需要增加新的耦合区域,也不需要建立新的逼近函数来保证界面位移的连续性.算例结果表明,该方法具有较好的计算精度.     

无网格Galerkin法与非协调元方法的耦合

无网格Galerkin法(EFGM)处理不可压缩问题时不存在自锁现象,有限元方法(FEM)也常被用来与其耦合以方便地施加边界条件和提高计算效率。在有限元方法中使用等参元,EFGM与FEM的耦合方法在处理不可压缩问题时仍然存在自锁现象。本文在有限元方法中,采用非协调元,将无网格kGalerkin法与非协调元耦合,保留了耦合方法的优点,且避免了求解不可压缩问题时的自锁现象。算例显示本文方法在分析平面应变不可压缩问题时能得到合理的结果。     

基于不连续场修正权的无网格断裂分析

论文采用不连续场修正权函数法来处理二维无网格断裂力学问题.相较于目前常用的无网格裂纹不连续性处理方案,采用修正权函数处理裂纹附近不连续场时只需要对原权函数进行修正,算法简便易实现.本文基于不连续场修正权函数的无单元Galerkin方法(EFGM),对在边界上施加Ⅰ-Ⅱ混合型裂纹位移场的斜裂纹板和带直裂纹受剪矩形板进行了数值分析,并与可视性准则、衍射法和透射法等不连续准则对比了裂尖位移场、应力场和应力强度因子解的数值精度,并分析了不同的M积分域尺寸、节点及背景网格等计算参数对裂尖奇异场数值解精度和形函数计算时间的影响.     

无网格Galerkin法GPU加速并行计算及其应用

针对无网格Galerkin法计算耗时的问题,采用逐节点对法来组装刚度矩阵、共轭梯度法求解基于CSR格式存储的稀疏线性方程组,提出了一种利用罚函数法施加本质边界条件的EFG法GPU加速并行算法,给出了刚度矩阵和惩罚刚度矩阵的统一格式,以及GPU加速并行算法的流程图。编写了基于CUDA构架平台的GPU程序,且在NVIDIA GeForce GTX 660显卡上通过数值算例对所提算法进行了性能测试与分析比较,探讨了影响加速比的因素。算例结果验证了所提算法的可行性,并在满足计算精度的前提下,其加速比最大可达17倍;同时线性方程组的求解对加速比起决定性影响。     

三维无网格Galerkin结构分析方法及其应用研究

针对无网格Galerkin法在三维复杂几何形状的结构分析中存在的刚度矩阵的稀疏存储实现难、六面体背景网格适应性差等问题,本文采用逐节点对法组装刚度矩阵,利用CSR格式存储刚度矩阵,提出了一种基于四面体背景积分的改进的三维无网格Galerkin法。通过采用罚函数法施加对称约束和给定位移值约束,并推导出了施加这两种位移约束的统一格式。利用所提算法完成了三维悬臂梁的计算,所得结果与其理论解相吻合;完成了轴流式风机轮毂的结构分析,得到的位移与应力分布结果与其有限元解相吻合。这表明本文所提方法能满足工程应用中的计算要求,并适用于具有复杂形状的几何模型分析。     

基于无网格Galerkin法求解矩形绝缘管道内的磁流体流动

无量纲磁流体流动控制方程中的哈特曼数较大将导致数值计算发散或误差过大。将无网格Galerkin法引入绝缘管道内的稳定磁流体流动计算中,针对磁流体控制方程中大哈特曼数导致计算误差大的情况,对无网格Galerkin法添加了稳定项。计算结果表明,同等条件下,添加了稳定项的无网格Galerkin法总体相对误差远小于标准无网格Galerkin法的结果,且可以计算哈特曼数最大达50000绝缘管道内的磁流体流动。     

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

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

伽辽金有限元素法对旋翼气弹稳定性的应用

采用伽辽金加权余数有限元素法发展了一种悬停状态下无铰旋翼桨叶气弹稳定性的分析方法。分析模型包括预锥角、下垂角、预掠角、总距角、桨根预安装角、桨叶预扭角、变距轴偏置、根部外伸量和操纵线系刚度等结构参数，对无铰旋翼桨叶气弹稳定性研究有普遍适用意义。试验证明该理论可行并能用于研究无铰旋翼结构参数对桨叶气弹稳定性的影响，也能用于直升机旋翼的型号设计。     

NEW ALGORITHM OF COUPLING ELEMENT-FREE GALERKIN WITH FINITE ELEMENT METHOD

IntroductionMeshless methods, as a special numerical method, originated from1970s. Since thediffuse element method was proposed by Nayroleset al.[1]in1992, the meshless methodshave received wide attentions in the mechanics area, and have shown obvious adv…     

冲击载荷下含表面裂纹圆柱壳体的动态断裂

动态载荷下含表面裂纹的有限尺寸构件的断裂问题在工程实践中有着重要意义,但由于此类问题非常复杂,目前还不能求得解析解。本文针对含轴向半椭圆盘状表面裂纹的圆柱壳体,应用有限元法研究了动态载荷下其断裂问题,计算了动态应力强度因子与静态应力强度因子的比值KdIyn(t)/KsIta。从计算结果可以得出,比值KdIyn(t)/KsIta与结构和裂纹的尺寸有关,而与冲击载荷的大小无关。本文所得结果在一定程度上揭示了圆柱壳体表面、裂纹面、物质惯性和弹性波的相互作用及其对动态断裂的影响。     

无单元法（EFGM）——在岩土工程上有限元法的有力补充

作为一种新的岩土工程数值计算方法，无网格伽辽金法具有局部化技术的无单元特性。该法与传统FEM相比，在精度和后处理等方面有着明显的计算优势。本文全面介绍了EFGM，并结合岩土工程领域中的土体固结变形问题，推导了EFGM的一维、二维固结问题的离散方程；同时对算例进行计算，将其结果与精确解和传统FEM解进行比较后，得出结论：精度是较高的，处理边界是准确的。尤其当与FEM联合使用时，既能发挥FEM的经验优质，又能体现EFGM的无网格特点，计算效果十分显著，因此，EFGM是FEM的有力补充。     

瞬态耦合热弹塑性接触有限元分析方法

在弹性接触问题有限元混合法的基础上,把材料非线性和表面非线性两种迭代过程耦合,在瞬态温度场分析中将伽辽金法和向后差分法结合,并用混合法进行热接触迭代,把瞬态温度场分析和弹塑性接触分析耦合。提出了一种瞬态耦合热弹塑性接触有限元分析方法,并已成功地用于核容器的密封分析。     

一种裂纹扩展的有限元仿真分析方法及其实现

为分析含多裂纹体在任意载荷下的断裂过程问题,给出了一种基于有限元的仿真方法及程序.基于商业有限元软件的应力分析结果和合适的断裂力学准则,判断裂纹的扩展以及扩展角;通过对每一个裂纹尖端扩展影响区域进行单元重划分,实现了多裂纹任意扩展的几何描述有限元模型更新.由于只对裂纹扩展区局部进行有限单元更新,网格划分和前期网格结果到新网格的映射的工作量相对较小,典型的裂纹扩展分析实例表明本文方法和程序的有效性.     

移动荷载作用下地基动力分析的有限元方法

通过对地基动力问题的基本方程进行变换,把基本方程变换到随荷载移动的运动坐标系中,通过加权残数法推导了相应的单元刚度矩阵,从而建立了移动问题的有限元格式,并发现移动荷载问题的单元刚度矩阵是对相应静力问题单元刚度矩阵的修正,在静力单元刚度矩阵的主对角元素上增加与移动速度有关的项,即可得到移动问题有限元的单元刚度矩阵,这样就将动力学问题转化为“拟静力”问题处理。文中用移动问题有限元方法计算了地基的动力响应,并与解析解进行了对比,以说明本方法具有较好的精度。     

应力波载荷作用下线弹性断裂过程的动态分析方法研究

利用Hopkinson单压杆实验装置 ,对材料的线弹性动态断裂特性进行了研究 ,建立了应力波载荷作用下动态裂纹起裂及扩展过程的动态分析方法 ,采用该方法可同时测得材料的动态裂纹起裂时间、断裂韧性和裂纹扩展速度。40Cr钢三点弯曲试样的实验结果表明 :该钢的动态裂纹扩展过程主要是减速过程 ,在2 2 5TPam /s的加载速率下 ,起裂时间为 2 8 0 0 s,最大裂纹扩展速度为 478 91m/s ,动态断裂韧性为6 3 12MPam。     

弹性力学轴对称问题的有限元线法

给出了解弹性力学空间轴对称问题的有限元线法的基本理论。该法包括了２－４条结线的等参数单元，沿结线方向的两点边值问题采用插值矩阵法解之。算例表明，本法具有良好的收敛性和较高的计算精度。