Dual reciprocity hybrid radial boundary node method for Winkler and Pasternak foundation thin plate

An efficient dual reciprocity hybrid radial boundary node method is developed for the analysis of Winkler and Pasternak foundation thin plate, in which a hybrid displacement variational principle, radial point interpolation method (RPIM) and dual reciprocity method (DRM) are combined. Firstly, the hybrid displacement variational principle is developed, in which the domain variables are interpolated by two groups of symmetric fundamental solutions, while the boundary variables are interpolated by RPIM instead of the traditional moving least square, and the shape function obtained by RPIM satisfies the delta function property, so boundary conditions can be applied directly. Besides, DRM is exploited to evaluate the particular solutions of inhomogeneous terms, which can be used to transform the domain integrals arising from the inhomogeneous term into equivalent boundary integrals. Finally, some additional equations based on the DRM theory are proposed to overcome the problem that the boundary integral equations are not enough to solve all variables. This method has the advantages of both no element mesh of meshless method and dimensionality reduction of boundary element method. Numerical examples of Winkler and Pasternak foundation plates are given to illustrate that the present method is effective, accurate and it can be further expanded into practical engineering.     

弹性动力学的双互易杂交边界点法

将双互易法同杂交边界点法相结合,提出了求解弹性动力问题的新型数值方法------双互易杂交边界点方法. 该算法在求解弹性动力问题时,将控制方程非齐次项的域内积分转化为边界积分. 该方法将问题的解分为通解和特解两部分,通解使用杂交边界点法求得,特解则使用局部径向基函数插值得到,从而实现了使用静力问题的基本解来求解动力问题. 计算时仅仅需要边界上离散点的信息,无论积分还是插值都不需要网格,域内节点仅用来插值非齐次项,因此该算法仍是一种边界类型的无网格方法. 数值算例表明,该方法后处理简单,计算精度高,适合于求解弹性动力问题.      

????????????????????????е????

提出了将杂交边界点法和双重互易法结合求解势问题的一种新的算法. 将势问题的解分为通解和特解两部分，通解使用 杂交边界点方法求解，特解则利用局部径向基函数近似. 该方法输入数据只是求解域上离散 的点，不需要额外的方程来计算域内物理量，后处理十分简便. 数值算例表明了该方法的稳 定性和有效性.     

NEW TREATMENT OF ESSENTIAL BOUNDARY CONDITIONS IN EFG METHOD BY COUPLING WITH RPIM

One of major difficulties in the implementation of meshfree methods using the moving least square (MLS) approximation, such as element-free Galerkin method (EFG), is the imposition of essential boundary conditions as the approximations do not pass through the nodal parameter values. Another class of meshfree methods based on the radial basis point interpolation can satisfy the essential boundary conditions exactly since its approximation function passes through each node in an influence domain and thus its shape functions possess the properties of delta function. In this paper, a coupled element-free Galerkin(EFG)-radial point interpolation method (RPIM) is proposed to enhance their advantages and avoid their disadvantages. Discretized equations of equilibrium are obtained in the RPIM region and the EFG region, respectively. Then a collocation approach is introduced to couple the RPIM and the EFG method. This method satisfies the linear consistency exactly and can maintain the stiffness matrix symmetric. Numerical tests show that this method gives reasonably accurate results consistent with the theory.     

弹性薄板弯曲问题的双互易杂交边界点法

基于一种板的修正变分泛函,将杂交边界点法与双互易法结合,用于薄板弯曲问题的分析。该方法将问题的解分为齐次方程的通解和非齐次的特解两部分,特解采用径向基函数插值得到,而通解则使用杂交边界点法求解。在杂交边界点法用于求解通解的列式过程中,边界变量采用移动最小二乘近似,域内变量则采用基本解插值。与有限元法相比,该方法仅需要边界上离散点的信息,无论插值还是积分都不需要网格,域内点仅用来插值非齐次项,因而仍是一种纯边界类型的无网格方法。数值算例表明,本文方法能以很少的计算自由度获得与其它方法同样的计算精度,且具有前后处理简单、收敛速度快等优点,适合于求解工程中各种薄板的弯曲问题。     

轴对称薄壁结构自由振动的边界元分析

采用双互易法分析薄壁轴对称结构自由振动的特征频率以及特征模态.首先,采用径向基函数插值域积分里的位移,利用双互易法将域积分转化为子午面边界的积分.然后,将边界物理量、基本解和特解展开为傅里叶级数,沿环向积分后得到的边界积分方程可用于轴对称结构带体积力问题和受非对称载荷的动力学分析,其积分域为轴对称结构子午面边界上的线积分,进一步降低了问题的维度和离散的难度.文章详细探讨了源点处于对称轴的特殊情况,根据基本解和特解的退化形式,针对无体积力和有体积力分别给出了处理奇异矩阵的方案.对于薄壁结构,采用双曲正弦变换处理近奇异积分有效提高积分精度.最后将双互易法和双曲正弦变化应用于薄壁轴对称结构带体积力的静力学和自由振动分析.数值结果表明,文章提出的处理奇异矩阵的方法能够有效处理源点处于对称轴的情况;当圆筒厚高比为$10^{-3}$,边界元计算的特征频率的相对误差为$10^{-3}$,且优于有限元的结果.      

基于径向基点插值法的旋转Mindlin板高次刚柔耦合动力学模型

将无网格径向基点插值法(radial point interpolation method,RPIM)用于中心刚体?旋转柔性板的动力学分析.基于浮动坐标系方法和一阶剪切变形理论即Mindlin板理论,考虑剪切变形的影响,并计入板面内变形的非线性耦合变形项,采用径向基点插值法描述板的变形场,保留动能中有关非线性耦合变形项...     

断裂力学问题的杂交边界点方法

提出了一种求解断裂力学的新的边界类型无网格方法-杂交边界点法.以修正变分原理和移动最小二乘近似为基础,同时具有边界元法和无网格法的优良特性,求解时仅仅需要边界上离散点的信息.该文将杂交边界点方法应用到弹性断裂问题中,将移动最小二乘方法中的基函数扩充,能更好的模拟裂纹尖端应力场的奇异性,推导了求解断裂力学的杂交边界点法方程,与传统的元网格方法相比,文中方法具有后处理简单,计算精度高的优点.数值算例表明了该方法的稳定性和有效性.     

二维水下声辐射问题的改进径向基点插值法研究

径向基点插值法是一种典型的无网格数值计算方法,在分析声学问题时,相比于传统有限元法能更好地抑制频散误差,且在相同的节点分布下通常可以得到更精确的数值解。本文提出一种改进的节点选取方案用于构造插值形函数,即改进径向基点插值法。该方案采取一个简单而直接的格式,可确保在进行数值积分时同一背景积分单元中的被积函数是连续可微的,从而减小数值积分误差,得到比原始径向基点插值法更精确的数值解。同时,为了处理外声场问题,本文采用DtN映射技术将无限域截断为有界计算域,满足索默菲尔德辐射条件。数值试验表明,相比于传统有限元法和原始径向基点插值法,本文改进方法具有更高的计算精度和计算效率,在研究水下声辐射问题时具有良好的应用前景。     

中厚板的耦合多项式基径向点插值无网格法

采用Mindlin平板理论,通过最小位能原理建立了各向同性中厚板的伽辽金整体弱式方 程,形函数采用耦合多项式基的径向点插值法构造,可以直接施加本质边界条件. 算例表明, 用耦合多项式基的径向点插值无网格法分析中厚板问题,具有效率高、精度高和易于实现等 优点,可以避免薄板弯曲时的剪切自锁现象.     

配点型广义节点无网格法的基本原理及其应用

借鉴流形方法思想,引入广义节点的概念,对传统的无网格法进行了改进,建立了可具有任意高阶多项式插值函数的广义节点无网格方法。同时采用径向插值函数构造具有插值特性的逼近函数;采用配点法建立系统的离散方程。在阐述了这种方法基本原理的同时,针对线弹性力学问题给出了这种方法的数值计算列式。与传统无网格方法相比,这种方法更具有一般性;同时由于采用了配点法而不需要背景积分网格,所以可以认为这种方法是某种真正意义上的无网格法。当选取0阶广义节点位移插值函数时便可得到传统的无网格法;在不增加支持域内节点数目的条件下,通过选取高阶广义节点位移插值函数可以提高计算精度。最后通过算例分析,对0阶、1阶及2阶广义节点无网格法与现有的有关解答进行了对比,论证了其合理性。     

APPLICATION OF THE DUAL RECIPROCITY BOUNDARY ELEMENT METHOD TO SOLUTION OF NONLINEAR DIFFERENTIAL EQUATION

In this paper the dual reciprocity boundary element method is employed to solve nonlinear differential equation ∇² u+u+ɛu ³ =b. Results obtained in an example have a good agreement with those by FEM and show the applicability and simplicity of dual reciprocity method(DRM)in solving nonlinear differential equations.     

径向基点插值法在旋转柔性梁动力学中的应用

将无网格径向基点插值法用于旋转柔性梁的动力学分析. 利用无网格方法对柔性梁的变形场进行离散,考虑梁的纵向拉伸变形和横向弯曲变形,并计入横向弯曲变形引起的纵向缩短,即非线性耦合项,运用第二类拉格朗日方程推导得到系统刚柔耦合动力学方程. 将无网格径向基点插值法的仿真结果有限元法和假设模态法进行比较分析,说明假设模态法的局限性,并表明其作为一种柔性体离散方法在刚柔耦合多体系统动力学的研究中具有可推广性,并讨论了径向基形状参数的影响. 同时运用3 种求解系统动力学方程的方法:纽马克方法、4阶龙格库塔法、亚当姆斯预报校正法,并比较各方法的计算效率, 结果表明纽马克方法最快.      

BEM for transient 2D elastodynamics using multiquadric functions

A new boundary element model for transient dynamic analysis of 2D structures is presented. The dual reciprocity method (DRM) is reformulated for the 2D elastodynamics by using the multiquadric radial basis functions (MQ). The required kernels for displacement and traction particular solutions are derived. Some terms of these kernels are found to be singular; therefore, a new smoothing technique is proposed to solve this problem. Hence, the limiting values of relevant kernels are computed. The validity and strength of the proposed formulation are demonstrated throughout several numerical applications. It is proven from the results that the present formulation is more stable than the traditional DRM, which uses the conical (1 + R) function, especially in predicting results in the far time zone.     

An efficient RPIM for simulating wave motions in saturated porous media

An efficient meshless technique for simulating wave motions in saturated porous media is introduced in this paper. Using radial point interpolation method (RPIM) and considering the precision requirement in wave motion simulation, a lumped-mass RPIM for saturated porous media with clear physical concepts is derived. Combining explicit time integration with the lumped-mass RPIM leads to a decoupled radial point interpolation method (decoupled RPIM) for dynamic analysis without solving algebra equation set. A two-dimensional problem is calculated by the proposed decoupled RPIM, the ordinary RPIM, and the finite element method (FEM). The accuracy and efficiency of these methods are compared, which demonstrate that the proposed decoupled RPIM has comparable accuracy with the ordinary RPIM and the finite element method, but has high efficiency in simulating wave motions in saturated porous media.     

AN IMPROVED HYBRID BOUNDARY NODE METHOD IN TWO-DIMENSIONAL SOLIDS

The hybrid boundary node method （HBNM） is a promising method for solving boundary value problems with the hybrid displacement variational formulation and shape functions from the moving least squares（MLS） approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the latter. Following its application in solving potential problems, it is further developed and numerically implemented for 2D solids in this paper. The rigid movement method is employed to solve the hyper-singular integrations. Numerical examples for some 2D solids have been given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method are studied through numerical examples.     

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

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

复合材料层合板波动特性分析的一种快速计算方法

将径向点插值方法应用在结构波动特性的快速分析中.在波数与方位角构成的二维参数空间中,通过初始规则节点布点和梯度自适应过程,对表征结构波动特性之一的群速度进行了插值计算,并把插值计算的结果与直接计算的结果进行比较,通过误差分析验证了径向插值方法在结构波动特性的快速计算中的有效性.     

不可压黏性流问题的无网格法分析

不可压缩黏性流问题一般采用Navier-Stokes方程来描述，基于加权残值法，推导了问题的无网格伽辽金法（EFGM）离散Navier-Stokes方程，在时间域上采用分步方法计算，速度和压力由相互独立的方程以解耦的形式求解，并采用同阶移动最小二乘近似，在每一时间步中，对压力解和速度解采用了Newton-Raphson迭代法进行修正，最后将所得到的方法应用到剪切驱动空腔流问题中，验证了方法的有效性，且解的精度高、稳定性好。     

一种高效的局部径向基点插值无网格方法

提出了一种弹性动力分析的高效局部径向基点插值无网格方法(MLRPI).该方法采用径向基点插值形函数近似解变量,运用局部Petrov-Galerkin法推导出了相应的离散方程,并根据波动模拟的精度要求,得到某一结点的动力方程.然后采用Newmark常平均加速度法和中心差分法相结合的显式积分格式进行时域积分,得到每个自由度的一种解耦递推格式.最后,对一平面应变问题进行了求解,比较了该文提出的解耦MI.RPI方法、常规MLRPI方法和ANSYS有限元方法的精度和计算时间,结果表明解耦MLRPI方法与常规MLRPI方法的精度相当,但计算效率大大提高.