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

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

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

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

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

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

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

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

瞬态热传导问题的复变量重构核粒子法

在重构核粒子法的基础上,引入复变量,讨论了复变量重构核粒子法.复变量重构核粒子法的优点是在构造形函数时采用一维基函数建立二维问题的修正函数.然后,将复变量重构核粒子法应用于瞬态热传导问题的求解,结合瞬态热传导问题的Galerkin积分弱形式,采用罚函数法引入本质边界条件,建立了瞬态热传导问题的复变量重构核粒子法,推导了相应的计算公式.与传统的重构核粒子法相比,复变量重构核粒子法具有计算量小、精度高的优点.最后通过数值算例证明了该方法的有效性. 关键词： 重构核粒子法 复变量重构核粒子法 修正函数 瞬态热传导问题     

轴对称变形核中配对粒子态的构成

量子多粒子体系的配对粒子态的构成与体系对称性密切相关.轴对称变形核中的配对粒子态可有三种构成方式,即单极对,±Ω和Signature对,分别是(J2,J3),J3,和(J32,R1(π))的共同本征态.在只计及对力时,三种形式的对力得出的能谱完全相同,常被等价地混用.但在推转壳模型中,由于三种配对态在空间转动下的性质不尽相同,能谱及有关性质随ω的变化将有不同程度的差异.     

高速碰撞数值计算中的光滑粒子法

扼要讨论了改进的光滑粒子法的离散思想,给出了二维轴对称问题中连续介质力学守恒方程的离散过程及离散格式,提出了轴对称坐标下确定影响域内粒子数的方法.最后通过高速碰撞的系列算例说明,光滑粒子法不但适宜于计算大变形冲击力学问题,而且有着其它网格法所无法替代的优势.     

基于平面问题的位移压力混合配点法

引入压力变量,将弹性力学控制方程表达为位移和压力的耦合偏微分方程组,采用重心插值近似未知量,利用重心插值微分矩阵得到平面问题控制方程的矩阵形式离散表达式.采用重心插值离散位移和应力边界条件,采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法求解过约束方程组,得到平面问题位移数值解.数值算例验证了所提方法的有效性和计算精度.     

再生核粒子法在辐射传输无网格计算中的应用

基于光滑粒子动力学(SPH)方法发展起来的再生核粒子法(RKPM)是一种拉格朗日型无网格方法。该方法在处理复杂形体内流动、导热、相变、应力应变等多物理场耦合问题时获得了较广泛的应用,但在辐射传输中的适用性尚无研究。而高温半透明材料的相变过程,除以上多物理过程外还需要考虑辐射传输过程的影响。本文基于建立高温半透明颗粒相变过程光热力耦合过程统一的无网格模拟方法的目的,重点研究了RKPM模拟辐射传输的可行性问题。文中首先建立了辐射强度的RKPM拟合公式,然后建立了辐射传输方程的直接配点离散格式。计算了一维参与性灰介质的无量纲热流和无因次温度分布。通过和文献值对比,初步验证了RKPM能够模拟辐射传输方程,且具有较好的计算精度。     

物理质团法:一个普适的数值方法体系

无网格方法根据分布于近邻空间各个方向的微元体物理信息构造离散方程,显著降低了空间导数计算对微元体本身及微元体之间拓扑结构的条件限制,极大提高了拉氏方法的大变形计算能力.由于不能利用微元体的完备几何信息,不容易构造符合物理的无网格算法,对那些物理参数存在间断的模型对象,难以获得稳定和准确的计算结果.本文基于对物理规律及数值模拟发展趋势的分析,提出符合物理且具有强普适性的无网格方法体系.基于该方法的一维算例表明,即使物理参数存在强烈间断,数值结果也能很好地逼近问题的真解.     

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.     

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.     

An improved interpolating element-free Galerkin method with nonsingular weight function for two-dimensional potential problems

In this paper, an improved interpolating moving least-square (IIMLS) method is presented. The shape function of the IIMLS method satisfies the property of Kronecker δ function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. And the number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square (MLS) approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin (IIEFG) method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin (EFG) method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has a higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.     

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.     

Two-dimensional fracture analysis of piezoelectric material based on the scaled boundary node method

A scaled boundary node method(SBNM) is developed for two-dimensional fracture analysis of piezoelectric material,which allows the stress and electric displacement intensity factors to be calculated directly and accurately. As a boundarytype meshless method, the SBNM employs the moving Kriging(MK) interpolation technique to an approximate unknown field in the circumferential direction and therefore only a set of scattered nodes are required to discretize the boundary. As the shape functions satisfy Kronecker delta property, no special techniques are required to impose the essential boundary conditions. In the radial direction, the SBNM seeks analytical solutions by making use of analytical techniques available to solve ordinary differential equations. Numerical examples are investigated and satisfactory solutions are obtained, which validates the accuracy and simplicity of the proposed approach.     

An improved interpolating element-free Galerkin method with a nonsingular weight function for two-dimensional potential problems

In this paper, an improved interpolating moving least-square (IIMLS) method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker δ function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square (MLS) approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin (IIEFG) method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin (EFG) method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.     

An improved boundary element-free method （IBEFM） for two-dimensional potential problems

The interpolating moving least-squares (IMLS) method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method (BEFM), combining the boundary integral equation (BIE) method with the IMLS method, the improved boundary element-free method (IBEFM) for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker δ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker δ function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.     

A complex variable meshless local Petrov-Galerkin method for transient heat conduction problems

On the basis of the complex variable moving least-square （CVMLS） approximation, a complex variable meshless local Petrov-Galerkin （CVMLPG） method is presented for transient heat conduction problems. The method is developed based on the CVMLS approximation for constructing shape functions at scattered points, and the Heaviside step function is used as a test function in each sub-domain to avoid the need for a domain integral in symmetric weak form. In the construction of the well-performed shape function, the trial function of a two-dimensional （2D） problem is formed with a one-dimensional （1D） basis function, thus improving computational efficiency. The numerical results are compared with the exact solutions of the problems and the finite element method （FEM）. This comparison illustrates the accuracy as well as the capability of the CVMLPG method.