弹性力学的复变量重构核粒子法

在重构核粒子法的基础上,提出了复变量重构核粒子法.复变量重构核粒子法的优点是采用一维基函数建立二维问题的修正函数.然后,将复变量重构核粒子法应用于弹性力学,提出了弹性力学的复变量重构核粒子法,并推导了相关公式.与传统的重构核粒子法相比,复变量重构核粒子法具有计算量小、效率高的优点.最后给出了数值算例证明了该方法的有效性. 关键词： 重构核粒子法 复变量重构核粒子法 弹性力学 无网格方法     

弹性力学的重构核粒子边界无单元法

将重构核粒子法(RKPM)和边界积分方程方法结合，提出了一种新的边界积分方程无网格方法——重构核粒子边界无单元法(RKP-BEFM).对弹性力学问题，推导了其重构核粒子边界无单元法的公式，研究其数值积分方案，建立了重构核粒子边界无单元法离散化边界积分方程，并推导了重构核粒子边界无单元法的内点位移和应力积分公式.重构核粒子法形成的形函数具有重构核函数的光滑性，且能再现多项式在插值点的精确值，所以本方法具有更高的精度.最后给出了数值算例，验证了本方法的有效性和正确性. 关键词： 重构核粒子法 弹性力学 边界无单元法     

弹性力学的插值型重构核粒子法

采用具有离散点插值特性的重构核粒子法形函数, 较精确地重构弹性体 变形的位移试函数, 再与弹性力学的最小势能原理相结合, 形成新的分析弹性力 学平面问题的插值型重构核粒子法. 由于插值型重构核粒子法形函数具有点插值特性和不低于核函数 的高阶光滑性, 因而既克服了多数无网格方法处理本质边界条件的困难, 也保证了较高的数值精度. 与早期的无网格方法相比, 本方法具有精度高、解题规模较小、可直接施加边界条件等优点. 通过对典型弹性力学问题数值模拟, 验证了所提方法的有效性和正确性.     

重构核粒子法对光滑粒子法的改进效果

光滑粒子法通过核函数进行近似估计,在计算域边界附近,核估计的精度明显下降.重构核粒子法通过校正函数对核函数进行重新构造,提高核估计方法在边界点和内点上对函数的估计精度以及计算的稳定性.研究发现,虽然校正函数的构造立足于对函数的精确估计,但这个优势同样能在对导数的估计中继续保持.通过理论研究及数学、物理模型的模拟,展示重构核粒子法的改进效果,揭示其能够提高精度的原因.     

改进的物理粘性SPH方法及其在溃坝问题中的应用

在低雷诺数物理粘性SPH方法基础上引入再生核粒子法进行密度重构,既避免了用人工粘性所导致的数值耗散问题,又提高了低雷诺数物理粘性SPH方法的数值稳定性;以溃坝问题为例,对比分析低雷诺数物理粘性SPH方法和本文方法的仿真结果表明,本文方法可有效消除数值不稳定,压强和速度分布更加光滑,粒子秩序更好,可应用于雷诺数较高或粘性不可忽略的流动问题.     

轴对称构件受力分析的插值粒子法

面对土木工程与机械工程中广泛存在的轴对称力学问题, 采用具有离散点插值特性的无网格方法形函数, 结合弹性力学空间轴对称问题的最小势能原理, 建立了轴对称构件力学分析的插值粒子法. 本文无网格法方法构造形函数不依赖网格, 也具有像有限元法一样可直接施加边界条件的优点. 本方法能直接获得全域连续应力场, 避免了有限元法应力后处理二次拟合带来的计算误差. 最后通过实例分析, 验证了所建立的无网格方法的有效性.     

瞬态热传导问题的一阶对称SPH方法模拟

为提高传统光滑粒子动力学(smoothed particle hydrodynamics, SPH)方法模拟瞬态热传导问题的精度和稳定性,本文提出了一种一阶对称光滑粒子动力学(first order symmetric SPH, FO-SSPH)方法.该方法将具有二阶热传导方程分解成两个一阶偏微分方程,然后基于梯度离散和Taylor级数展开思想,对一阶核梯度形式进行修正,并将得到的局部矩阵对称化.数值结果表明:与传统SPH方法相比,FO-SSPH方法精度高、数值稳定性好; 该方法能较准确地直接施加混合边值 关键词： 瞬态热传导 光滑粒子动力学 非线性     

弹性力学的复变量无网格局部 Petrov-Galerkin 法

复变量移动最小二乘法构造形函数, 其优点是采用一维基函数建立二维问题的试函数, 使得试函数中所含的待定系数减少, 从而有效提高计算效率. 文章基于复变量移动最小二乘法和局部Petrov-Galerkin弱形式, 采用罚函数法施加边界条件, 推导相应的离散方程, 提出弹性力学的复变量无网格局部Petrov-Galerkin法. 数值算例验证了该方法的有效性.     

基于修正并行光滑粒子动力学方法三维变系数瞬态热传导问题的模拟

为了解决传统光滑粒子动力学(SPH)方法求解三维变系数瞬态热传导方程时出现的精度低、稳定性差和计算效率低的问题,本文首先基于Taylor展开思想拓展一阶对称SPH方法到三维热传导问题的模拟,其次引入稳定化处理的迎风思想,最后基于相邻粒子标记和MPI并行技术,结合边界处理方法得到一种能够准确、高效地求解三维变系数瞬态热传导问题的修正并行SPH方法.通过对带有Direclet和Newmann边界条件的常/变系数三维热传导方程进行模拟,并与解析解进行对比,对提出的方法的精度、收敛性及计算效率进行了分析;随后,运用提出的修正并行SPH方法对三维功能梯度材料中温度变化进行了模拟预测,并与其他数值结果做对比,准确地展现了功能梯度材料中温度变化过程.     

弹性力学的复变量无网格方法

在移动最小二乘法的基础上，提出了复变量移动最小二乘法.复变量移动最小二乘法的优点是采用一维基函数建立二维问题的逼近函数，所形成的无网格方法计算量小.然后，将复变量移动最小二乘法应用于弹性力学的无网格方法，提出了复变量无网格方法，推导了复变量无网格方法的公式.与传统的无网格方法相比，复变量无网格方法具有计算量小、精度高的优点.最后给出了数值算例. 关键词： 移动最小二乘法 复变量移动最小二乘法 无网格方法 弹性力学 复变量无网格方法     

The complex variable reproducing kernel particle method for two-dimensional elastodynamics

On the basis of the reproducing kernel particle method (RKPM), a new meshless method, which is called the complex variable reproducing kernel particle method (CVRKPM), for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.     

Analysis of variable coefficient advection-diffusion problems via complex variable reproducing kernel particle method

The complex variable reproducing kernel particle method （CVRKPM） of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape function of a two-dimensional problem is formed with a one-dimensional basis function. The Galerkin weak form is employed to obtain the discretized system equation, and the penalty method is used to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional variable coefficient advection-diffusion problems are obtained. Two numerical examples are given to show that the method in this paper has greater accuracy and computational efficiency than the conventional meshless method such as reproducing the kernel particle method （RKPM） and the element- free Galerkin （EFG） method.     

The complex variable reproducing kernel particle method for elasto-plasticity problems

On the basis of reproducing kernel particle method(RKPM),using complex variable theory,the complex variable reproducing kernel particle method(CVRKPM) is discussed in this paper.The advantage of the CVRKPM is that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is formed.Then the CVRKPM is applied to solve two-dimensional elasto-plasticity problems.The Galerkin weak form is employed to obtain the discretized system equation,the penal...     

Meshless analysis of three-dimensional steady-state heat conduction problems

Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson’s equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.     

Combining the complex variable reproducing kernel particle method and the finite element method for solving transient heat conduction problems

In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.     

An interpolating reproducing kernel particle method for two-dimensional scatter points

An interpolating reproducing kernel particle method for two-dimensional(2D) scatter points is introduced. It eliminates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating reproducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.