弹塑性大变形分析的一致性高阶无单元伽辽金法

采用无单元伽辽金法求解弹塑性大变形问题。充分利用无单元法易于建立高阶近似函数的优点,位移采用二阶移动最小二乘近似。在更新拉格朗日方法的框架下,通过对控制方程弱形式的线性化建立了内力率的表达式,并区分为材料和几何两部分。采用Hughes-Winget算法更新应力,建立了Newton-Raphson迭代求解所需的一致切线刚度阵。刚度阵的数值积分采用近来针对小变形分析建立的二阶一致三点积分格式QC3（Quadratically Consistent 3-point integration scheme）。数值结果证明了本文方法分析弹塑性大变形问题的有效性和优越性。     

自适应一致性高阶无单元伽辽金法

近来提出的一致性高阶无单元伽辽金法通过导数修正技术大幅度减少了所需积分点数目,并能够精确地通过线性和二次分片试验,显著改善标准无单元伽辽金法的计算效率、精度和收敛性.本文在此基础之上,充分利用无单元法易于在局部区域添加节点的优势,发展了一致性高阶无单元伽辽金法的h型自适应分析方法.根据应变能密度梯度该方法自适应地确定需节点加密的区域,基于背景积分网格的局部多层细化要求生成新的计算节点,同时考虑了节点分布由密到疏渐进过渡的情形.采用相邻两次计算的应变能的相对误差作为自适应过程的停止准则,将所发展自适应无网格法应用于由几何外形、边界外载和体力等因素造成的应力集中问题的计算分析.数值结果表明,所发展方法能够自适应地对高应力梯度区域进行节点加密,自动给出合理的计算节点分布.与已有的标准无网格法的自适应分析相比,所发展方法在计算效率、精度和应力场光滑性等方面均展现出显著优势.与采用节点均匀分布的一致性高阶无单元伽辽金法相比,它大幅度地减少了计算节点数目,有效提高了一致性高阶无单元伽辽金法在分析应力集中等存在局部高梯度问题时的计算效率和求解精度.     

弹塑性体脆断相场模型的一致性无单元伽辽金法

采用无单元伽辽金法（EFG）对弹塑性体脆性断裂的相场模型进行了数值实现。利用无单元法便于构建高阶近似函数的优势,位移和相场均采用二阶移动最小二乘（MLS）近似。刚度阵的数值积分采用更为高效的二阶一致三点积分格式QC3（Quadratically Consistent 3-point integration scheme）。本构算法采用Newton-Raphson迭代和弹塑性一致性切线模量。数值结果表明了本文方法模拟弹塑性体脆性断裂的有效性。     

无单元伽辽金法的并行计算

对无单元伽辽金法的并行计算进行了详细研究,并将其应用于弹性动力学问题。使用并行桶搜索算法进行节点搜索,使用并行几何搜索算法进行样点搜索,讨论了移动最小二乘MLS(Moving Least Squares)形函数及其导数的并行计算和方程组的并行求解,并利用多层图形划分实现负载平衡。最后给出了并行无单元伽辽金法应用于弹性动力学的计算流程和实例。计算结果表明无单元伽辽金法具有很高的并行性和很好的并行效率,对其进行并行计算具有非常重要的意义。     

反平面断裂问题的无单元伽辽金比例边界法

将比例边界法与无单元伽辽金法相结合,建立了反平面断裂分析的无单元伽辽金比例边界法。这是一种边界型无网格法,在环向方向上采用无单元伽辽金法进行离散,因此计算时仅需要边界上的节点信息,不需要边界元所要求的基本解。为了便于施加本质边界条件,通过建立节点值和虚拟节点值之间的关系给出了修正的移动最小二乘形函数。在径向方向上,该方法利用解析的方法求解,因此是一种半解析的数值方法。最后,给出了数值算例,并验证了所提方法后处理简单和计算精度高的特点,适合于求解反平面断裂问题。     

粘弹性问题的插值型无单元伽辽金比例边界法

作为一种最近发展起来的半解析数值方法,插值型无单元伽辽金比例边界法不仅无需基本解,且在处理应力奇异性问题和无限域问题时十分有效．为了更有效地求解粘弹性问题,对插值型无单元伽辽金比例边界法应用于此类问题进行了研究,并发展了相应的算法. 通过时域分段展开,将时空耦合的初边值问题转化为一系列递推形式的边值问题,然后采用插值型无单元伽辽金比例边界法进行自适应计算．在径向保持解析特性的基础上,环向采用无单元伽辽金法离散可简化前处理和后处理工作量．此外,改进的插值型移动最小二乘法形函数具有插值性,有效地解决了本质边界条件不能直接施加的困难．最后给出了数值算例,并验证了所提方法的有效性和正确性.     

质点积分无单元伽辽金法及其在金属挤压过程中的应用

无单元伽辽金法需要在背景网格上积分,计算量大.节点积分无单元伽辽金法把对求解域的积分转化为对节点的求和,效率高,但因零能模态不受控制而会产生不稳定现象,需要采取一定的稳定化方案.本文采用应力点思想,通过Newtor-Cotes法计算积分,建立了质点积分无单元伽辽金法,并通过小变形弹性静力学问题说明了该方法具有良好的稳定性,且计算效率远高于无单元伽辽金法.最后本文将质点积分无单元伽辽金法成功地应用于三维金属挤压成型过程的数值模拟,显示了该方法在分析此类问题时的优势和潜力.     

无单元伽辽金法在船体开孔板格弹性屈曲分析中的应用

针对船体开孔板格弹性屈曲问题,评估了无单元伽辽金法(EFG)的适用性。编程计算分别基于MLS近似和移动Kriging插值两种形函数构造方式进行,以ABAQUS有限元解为基准确定了两种算法最佳内置参数,对比分析了计算耗时和相对误差。结果表明:MLS-EFG算法在不同的板格尺寸、开孔位置、载荷条件下的结果都比较准确,具有较好的适用性;而MK-EFG算法的计算速度更优,但在纯剪切屈曲分析时计算结果存在一定程度的正偏差。本研究可为无网格法在船体结构设计中的推广应用提供参考。     

基于三变量伽辽金无单元法的杆系模型动力计算与参数分析

将移动最小二乘法与广义移动最小二乘法相结合,提出同时考虑节点轴向变形、切向变形和转角的三变量无单元插值形函数,应用于框架结构杆系模型的动力计算。运用伽辽金弱形式推导离散动力方程,采用罚函数法引入边界条件,建立三变量无单元法动力计算公式,采用Newmark-β法进行动力时程分析。讨论了节点数量、影响半径以及基函数阶次对计算结果的影响。研究结果表明,当节点的间距≤10%杆件长度,影响半径scale乘子界于2.5~4.0,采用三阶多项式基函数时,无单元的动力计算可获得较好的精度。     

几何非线性分析的无网格伽辽金算法

利用几何非线性的应变-位移关系,在小应变假设的条件下,推导出二维几何非线性问题中的无网格伽辽金法的计算格式。由于无网格方法中的形函数不具备Kronecker delta性质,文中采用罚方法来实现本质边界条件。数值实例表明,无网格伽辽金法在处理几何非线性问题时,具有较高的计算精度,是一种有效的数值计算方法。     

波前法在无网格伽辽金法中的应用

有别于有限元法,无网格法采用基于点的近似,可彻底或部分地去除网格(只保留积分所需的背景网格),在保证计算精度同时降低计算难度.无网格伽辽金法(Element Free Galerkin Method,EFG)是一种基于移动最小二乘近似(Moving Least-Squares,MLS)的全局弱式无网格法,广泛应用于计算力学等领域,该方法的一个缺点是:计算过程中产生的系数矩阵含有的非零元数量比有限元法多,即使处理中等规模模型时,也要求计算机有很大的存储空间,并且计算时间长.波前法在有限元法中已有很成熟的应用,但至今没有应用于无网格方法.论文介绍了波前法在无网格伽辽金法中的应用方法,编写了相应的计算程序,并以弹性力学为例做了验算.     

无网格伽辽金法求解平面偶应力问题

提出采用无网格伽辽金法(EFGM)求解偶应力问题,以避免有限元求解中因C¹连续要求可能引起的不便。推导了基于二次基和移动最小二乘技术的EFGM计算公式,通过计算受轴向均匀拉伸的带中心小孔无限大板和细长杆的偶应力问题,对所提方法进行了数值验证,与解析解相比结果令人满意,此外还讨论了节点密度和权函数对计算结果的影响。数值结果表明,所提算法可有效地求解平面偶应力问题。     

复合材料层合板的谱随机无网格伽辽金法

针对基于摄动理论的计算方法不适用于分析随机结构参数大变异问题,提出谱随机无网格伽辽金法,该方法基于随机场正交分解理论,将随机场采用Karhunen-Loève级数展开为一系列不相关随机变量,再引入结构位移随机响应的混沌多项式分解,结合无网格伽辽金法,从而导出含随机变量的复合材料层合板的谱随机无网格伽辽金法,给出结构响应的统计特征值的计算公式,该方法既适用于随机参数大变异情况,又具有无网格法的优势;数值算例结果表明该方法是正确有效的。     

动力弹塑性分析的无网格自然单元法

基于无网格自然单元法,提出了结构动力弹塑性响应分析的一条新途径.自然单元法是一种新兴的无网格数值计算方法,其实质是基于自然邻近插值的伽辽金法.自然单元法在本质边界条件的施加上较采用移动最小二乘法的无网格法具有明显的优势.在空间域上采用自然单元法离散,并运用加权余量法推导了动力弹塑性分析的离散控制方程.然后,采用预校正形式的Newmark法在时间域上进行求解.最后给出了数值算例,并验证了所提方法的有效性和正确性.     

伽辽金型无网格法的数值积分方法

无网格法直接通过节点信息构造形函数,不依赖于节点之间的有序单元连接,能够建立任意高阶连续的整体协调形函数.与传统的有限元法相比,无网格法对大变形问题、移动边界问题和高阶问题的求解有比较明显的优势.伽辽金型无网格法是目前应用最为广泛的一类无网格法.虽然无网格形函数本身不依赖于单元,但伽辽金型无网格法需要采取合适的方法进行弱形式的数值积分.由于无网格形函数一般不是多项式,具有非插值性且影响域与背景积分网格通常不重合,伽辽金型无网格法通常需要采用高阶的高斯积分进行数值积分,导致了计算效率低下,难于求解大型实际问题.因此,如何通过建立高效积分方法提高无网格法的计算效率成为无网格法研究领域的一个核心问题.论文总结了伽辽金型无网格法中若干常用的数值积分方法,并对伽辽金型无网格法的数值积分方法领域存在的一些问题进行了探讨.     

混合变换法在无网格伽辽金方法中的应用

移动最小二乘近似的非插值特征给无网格伽辽金方法的应用带来了一定的的困难，本文将再生核质点法中的混合变换法与无网格伽辽金方法相结合，通过对移动最小二乘近似进行局部修正，实现了无网格伽辽金方法中本质边界条件的精确施加。对权函数、影响半径、积分阶等对计算精度的影响进行了有益的探讨。数值计算结果表明了方法的可行性。     

伽辽金法在悬臂梁大挠度问题中的应用

伽辽金法在悬臂梁大挠度问题中的应用胡辉（邵阳高等专科学校,邵阳４２２００４）扁平弹簧一类的柔性杆,即使在载荷不大的情况下,其最大挠度与杆长也是同一数量级的,这时如仍用小挠度理论计算,将会产生很大的误差,因此研究梁的大挠度问题是有重要意义的．本文讨论如...     

关于伽辽金法的一点注记

伽辽金法作为一种数值方法,广泛用于各种数学物理工程问题.教科书和一些文献关于伽辽金法试函数的选取存在一些争议,本文通过应用该法求解轴向力作用下悬臂梁的静动力学问题为例,证实仅选取满足位移而不满足力边界条件的试函数时,即使增加试函数的个数,可能也无法获得正确的计算结果.     

完全变换法在无网格伽辽金方法中的应用

由于移动最小二乘形函数一般不具有常规有限元或边界元形函数所具有的插值特征，本质边界条件的处理成为无网格伽辽金法实施中的一个难点。本文通过建立节点位移和广义位移之间的关系对移动最小二乘形函数进行修正，给出了修正的移动最小二乘形函数；以二维问题为例，对完全交换法在无网格伽辽金方法中的应用进行了研究，实现了本质边界条件在节点处的精确施加。数值计算结果表明该方法不仅简单合理，而且具有较高的精度、收敛性和稳定性。     

二维弹性力学问题的光滑无网格伽辽金法

计算效率低的问题长期阻碍着无网格伽辽金法(element-free Galerkin method, EFGM)的深入发展.为了提高EFGM的计算速度,本文提出一种求解二维弹性力学问题的光滑无网格伽辽金法.该方法在问题域内采用滑动最小二乘法(moving least square, MLS)近似、在域边界上采用线性插值建立位移场函数;基于广义梯度光滑算子得到两层嵌套光滑三角形背景网格上的光滑应变,根据广义光滑伽辽金弱形式建立系统离散方程.两层嵌套光滑三角形网格是由三角形背景网格本身以及四个等面积三角形子网格组成.为了提高方法的精度,由Richardson外推法确定两层光滑网格上的最优光滑应变.几个数值算例验证了该方法的精度和计算效率.数值结果表明,随着光滑积分网格数目的增加,光滑无网格伽辽金法的计算精度逐步接近EFGM的,但计算效率要远远高于EFGM的.另外,光滑无网格伽辽金法的边界条件可以像有限元那样直接施加.从计算精度和效率综合考虑,光滑无网格伽辽金法比EFGM具有更好的数值表现,具有十分广阔的发展空间.