无网格方法在动力学中的应用

再生核质点法是近几年研究出的一种新型无网格方法,该方法具有只需要质点信息无需划分单元的无网格特性,具有计算优势.介绍了再生核质点法,并将其应用到非线性动力学研究.动力学过程涉及多重非线性,假设动力学中变形属于小变形情况,同时考虑材料非线性的前提下,通过引入增量型的材料本构模型,采用完全Lagrangian计算格式,推导了动力学的再生核质点法计算控制方程.通过算例验证了该方法在动力学问题中的有效性.     

无网格伽辽金法在施工导流二维数值模拟中的应用

针对有限元法等传统数值计算方法存在受单元网格限制、前后处理工作复杂的问题,提出应用一种数值计算方法--无网格伽辽金法,对具有简单边界条件的水利水电工程施工导流的恒定二维浅水流动问题进行了分析、计算.同时利用有限元法进行了对比计算,从流速、水位的计算结果来看,两种计算方法结果相近、误差较小,表明采用无网格伽辽金法解决此类问题是可行的.     

具有动态接触边界自动识别功能的刚(粘)塑性无网格RKPM法

将刚(粘)塑性流动理论与再生核质点法(RKPM)相结合,提出了基于刚(粘)塑性不可压缩材料的无网格RKPM法.分别采用边界奇异权法和修正的罚函数法处理本质边界条件和体积不可压缩条件,根据刚(粘)塑性材料的不完全广义变分原理,推导了金属塑性成形过程无网格RKPM法数值模拟的刚度方程.将动态接触边界自动识别技术引入无网格法.最后,对金属塑性成形过程进行了数值模拟,并将模拟结果与刚塑性有限元商品软件Deform2D及实验结果作了比较.无网格RKPM法计算结果与有限元结果和实验结果均吻合良好,表明了该方法的正确性.     

再生核质点法在非线性冲击动力学中的应用

利用再生核方法来构造计算域函数的插值函数，通过区域内的质点务数来实现数值计算，满足一致性要求；在大变形、大转动、大应变的情况下，引入Jaumann应力和应变速率张量来表征实际应力和应变速率，进而可以表征出材料冲击条件下的本构关系；将本构关系引入冲击过程虚功原理，采取Lagrangian描述的方法，并通过再生核质点方法将能量积分方程离散，推导了非线性冲击过程的再生核质点法计算控制方程，给出了方程求解的Newmark递推公式和Newton-Raphson迭代方法．通过具体实例来说明了再生核质点法在非线性动力学中的应用过程，并验证了计算方法的正确性．     

轴对称几何非线性问题中的无网格算法

应用无网格伽辽金法对轴对称几何非线性问题进行了分析。在小变形假设的条件下，利用几何非线性的应变-位移关系，基于线性弹性本构关系，推导了无网格法的计算控制方程，并采用Newton—Raphson迭代法来求解非线性方程，初步计算了压力管道的几何非线性问题。由于无网格方法中的形函数不具备Kroneckerdelta性质，采用罚方法来实现本质边界条件。数值实例表明．无网格伽辽金法在处理轴对称几何非线性问题时，具有较高的计算精度，是一种有效的数值计算方法。     

无网格伽辽金-有限元耦合方法在温度场中的应用

为了充分利用无网格法和有限元法的优点,将无网格伽辽金-有限元耦合方法用于分析温度场问题.根据无网格伽辽金-有限元耦合计算原理得出了耦合区域的形函数,从能量泛函弱变分形式中得到控制方程,从而求出数值解.EFGM-FE耦合法克服了单纯使用无网格法带来的边界条件难处理及计算效率较低的缺点.数值算例表明了这种方法是可行的,有效的.     

侧向挤压过程的耦合无网格-有限元法数值模拟

应用刚塑性耦合无网格-有限元法,对平面侧向挤压过程进行数值模拟.采用修正的罚函数法约束体积不可压缩条件,避免体积闭锁现象.采用边界奇异核方法,直接、准确地施加本质边界条件.针对侧向挤压工艺的特点,将无网格节点布置在变形剧烈的区域,充分发挥其处理大变形问题的能力;在变形较小的区域采用有限元单元,利用其较高的计算效率.数值模拟结果表明,该耦合方法对于侧向挤压这类局部成形问题的求解具有良好的效果.     

Euler梁自由振动的Hermite再生核无网格分析

基于Hermite再生核无网格近似,建立了Euler梁自由振动分析的伽辽金无网格离散方程.针对常见的几种典型边界条件的Euler梁自由振动问题,详细分析了前两阶频率的误差和收敛性.结果表明,与传统仅采用挠度近似的伽辽金无网格法和Hermite有限元法相比,考虑节点转角对挠度近似影响的Hermite无网格方法具有更高的精度,为Euler梁振动分析提供了一种高精度的数值方法.     

无网格RKPM法及其在体积成型中的应用

介绍了无网格再生核质点法(RKPM)的基本原理及其在刚塑性可压缩材料体积成型中的应用·采用罚函数满足本质边界条件,模拟了理想塑性材料的二维平面无摩擦镦粗过程,计算了不同坯料高宽比在不同可压缩因子下的变形情况,将计算结果与解析解进行比较,结果表明变形工件形状和应变参数与解析解吻合较好,应力计算的误差可能与点积分有关·模拟结果说明了再生核质点法在体积成型中应用的可行性和有效性·     

次弹性材料的无网格分析方法研究

次弹性材料在实际工程中是常见的,传统计算中大多数采用有限元方法。利用无网格伽辽金法对次弹性材料进行数值计算,并通过罚参数来实现本质边界条件,推导了求解此类问题的无网格伽辽金法离散格式。算例结果表明,无网格伽辽金法处理次弹性材料时,具有较高的计算精度,是一种有效的数值计算方法。     

基于RKPM梁固有频率及薄板柔度的拓扑优化

以节点相对密度为设计变量,以固有频率最大为目标函数,通过修正低密度区质量矩阵建立了基于重构核粒子法(RKPM)的结构动力拓扑优化数学模型.采用罚函数法施加本质边界条件,利用直接微分法推导了结构固有频率灵敏度方程,同时研究了受横向载荷弯曲的基尔霍夫薄板柔度最小的拓扑优化问题.最后对比分析了节点依赖性以及设计变量对最优拓扑结构的影响,并结合以上算法和优化准则法编写程序完成了2个拓扑优化算例.优化结果表明:所建立的模型不仅能有效抑制局部模态和重特征频率的出现,而且因通过重构核近似提高了计算点密度场的连续性,棋盘格现象得以消除,可以得到清晰光滑的拓扑边界.     

无网格伽辽金法(EFGM)求解接触问题

本文考虑到无网格伽辽金法只需节点信息,不需将节点连成单元,并且有精度高,后处理方便等优点,从而用它求解接触问题。这里采用的方法是将Ｋａｔｏｎａ界面单元引入ＥＦＧＭ,迭代求得两物体间的接触状态。算例表明,本文方法基本可行。     

具有混合边界条件的地下水污染模型的混合元——特征有限元方法

本文对具有混合边界条件的地下水污染模型问题提出了混合元—特征有限元全离散格式，即对水量方程采用混合元格式，而对浓度方程沿特征方向有限元离散。利用椭圆投影及其误差估计，建立了计算格式在Ｈ模意义下的最优误差估计；数值算例及结果分析验证了该方法的正确有效性     

Signorini问题的插值型边界无单元法

本文将改进的移动最小二乘插值法和边界积分方程结合,提出了求解Signorini问题的一种新的边界类型无网格方法——插值型边界无单元法.该方法用投影算子处理Signorini问题中的非线性边界不等式条件,然后将Signorini问题归化为边界积分方程,并用改进的移动最小二乘插值法近似未知的边界变量,然后本文分析了该方法的收敛性.数值算例表明该方法在求解Signorini问题时的可行性和有效性,相对于边界元方法也具有更好的精度和收敛速度.     

采用无网格伽辽金法求解电容层析成像正问题

由于有限元法求解电容层析成像正问题的计算准备及后处理非常费时，对正问题的三维求解造成了瓶颈，为此，提出采用无网格伽辽金法求解电容层析成像正问题，获得正问题的弱变分形式，并用拉格朗日乘子法施加边界条件，从而得到数值解．在同样的仿真条件下，2种方法的计算时间分别为14．046S和5．078S．对5种典型流型进行仿真，结果表明，2种方法计算结果的最大相对误差为2．25％．因此，无网格伽辽金法与有限元法具有相当的精度，且计算速度有较大提高．     

非线性热传导问题的基于滑动Kriging插值的MLPG法

利用基于滑动Kriging插值的无网格局部Petrov-Galerkin(MLPG)法来求解二维非线性稳态和瞬态热传导问题,Heaviside分段函数作为局部弱形式的权函数,并通过加权余量法推导相应的离散方程.该问题考虑了材料热传导系数随温度的线性变化,并通过拟线性法来求解非线性问题的解,时间域的离散通过向后差分法来实现.基于滑动Kriging插值构造MLPG中的形函数由于满足克罗内克δ性质,因此可以直接准确地施加本质边界条件.在构造刚度矩阵过程中,只涉及边界积分,不涉及区域积分和奇异积分.将数值计算结果与有限元法得到的结果加以对比可以看出,基于滑动Kriging插值的MLPG法能够很好地解决此类热传导问题.