插值型无单元Galerkin比例边界法与有限元法的耦合在压电材料断裂分析中的应用

插值型无单元Galerkin比例边界法是一种只需在边界上采用插值型无单元Galerkin法离散且无需基本解的半解析方法,能有效求解压电材料的断裂问题.为进一歩提高这种方法的适用性,该文提出了一种用于压电材料断裂分析的插值型无单元Galerkin比例边界法耦合有限元法(finite element method,FEM)的分析方法.裂纹周边一定范围的计算域采用插值型无单元Galerkin比例边界法离散,其余区域采用FEM离散.插值型无单元Galerkin比例边界法方程和FEM方程的耦合可利用界面两侧广义位移的连续条件方便地实现.最后,给出了两个数值算例验证了该文所提方法的有效性.     

三维稳态对流扩散问题的无网格求解算法研究

无网格法是一种不需要生成网格就可模拟复杂形状流场计算的流体力学问题求解算法.为了提高基于Galerkin弱积分形式的无网格方法求解三维稳态对流扩散问题的计算效率,提出了在空间离散上采用基于凸多面体节点影响域的无网格形函数,并通过选取适当节点影响半径因子避免节点搜索问题,同时减少系统刚度矩阵带宽.计算中当节点影响因子为1.01时,无网格方法的形函数近似具有插值特性且本质边界条件的施加与有限元一样简单.三维立方体区域的稳态对流扩散数值算例表明:在保证计算精度的同时,采用凸多面体节点影响域的无网格方法比传统无网格方法最高可节省计算时间42%.因此从计算效率和精度考虑,在运用无网格方法求解三维问题时建议采用凸多面体节点影响域的无网格方法.     

无网格Chebyshev配点法求解二维位势问题

基于Chebyshev正交多项式插值理论和无网格配点技术,提出一种新型的无网格数值离散方法,称之为Chebyshev配点法.所提方法采用Chebyshev多项式的零点(Gauss-Lobatto节点)为插值节点,可最大限度地降低龙格现象,并且提供插值多项式的最佳一致逼近.数值算例表明,本文算法稳定,效率高,并可达到很高的计算精度.     

耦合Navier-Stokes/Darcy模型多重网格方法的有限元实现(英文)

本文研究耦合Navier-Stokes/Darcy模型问题.构造一种从粗网格到细网格的有限元空间插值方法,不但简化了数值积分的单元匹配,也保证了数值积分的精度.利用基于有限元空间的多重网格方法,获得与直接法求解耦合问题误差相同的收敛阶,推广两重网格方法的结果.     

基于虚节点的多边形有限元法

虚节点法是一种新的基于单位分解理论的多边形有限元法.将虚节点法应用于求解弹性力学问题,并且通过大量数值实验测试虚节点法的计算效果.因为虚节点法具有多项式形式,所以有效地降低了传统多边形有限元法的积分误差.数值实验证明,在分片实验中虚节点法能得到比包括Wachspress法和mean value法在内的传统多边形有限元法更精确的数值结果.在收敛性试验中,虚节点法在相同节点数的条件下能取得比三角形一次单元更精确的数值结果.因为虚节点法能适应任意边数的多边形单元,所以对网格具有很强的适应性,在几何条件复杂、网格生成困难的问题中具有良好的应用价值.为了展示虚节点法潜在的应用价值,用虚节点法求解断裂力学应力强度因子和模拟裂纹扩展.同时,基于多边形单元的网格重划分技术和网格加密技术也应用于求解断裂力学应力强度因子和模拟裂纹扩展.     

计及时变演化特征的硅泡沫垫层非线性黏弹性模型研究

基于黏弹性基本理论,引入材料非线性特征,考虑了材料加载、保载应力松弛历史、老化效应以及黏弹性模型各运动单元退化的差异性,并从两种老化机制出发,获得了老化硅泡沫垫层力学模型以及长时应力松弛硅泡沫垫层接续加载力学模型.模型机理清晰,能够反映材料服役历史信息及其对力学效应的影响.     

直角坐标系下层状黏弹性地基Biot固结解析刚度矩阵解

基于直角坐标系下Biot固结的基本控制方程,并考虑软土土骨架的黏弹性特性,通过Fourier-Laplace积分变换、解耦变换、微分方程组理论和矩阵理论,推导了黏弹性地基Biot固结三维空间问题和平面应变问题在积分变换域的解析解,进而得到对应问题的单元刚度矩阵.然后根据对号入座原则组装得到层状黏弹性地基Biot固结对应问题的总体刚度矩阵.通过求解总体刚度矩阵形成的线性代数方程,得到层状黏弹性地基Biot固结对应问题在积分变换域内的解答.最后应用Fourier-Laplace逆变换得到其物理域内的解.对比求解黏弹性Biot固结问题退化的弹性Biot固结问题与已有解答,验证了刚度矩阵计算方法的正确性,为层状黏弹性地基Biot固结问题提供了理论基础.     

二维Helmholtz方程的插值型边界无单元法

针对二维Helmholtz方程的内外边值问题,提出了插值型边界无单元法(interpolating boundary element-free method).在间接位势理论的基础上,利用Laplace方程基本解的特性,建立了求解Helmholtz方程Neumann边值内外问题的正则化形式,有效消除了强奇异积分的计算.其次通过引入全局距离展开成局部距离的幂级数,详细推导了距离函数的导数和法向导数差值的极限表达式.最后给出了4个插值型边界无单元法的数值算例,表明了该方法可取得较高的可行性和有效性.     

轴对称弹性动力学问题的重构核插值法

重构核插值法是近年来提出的一种新型无网格方法.该方法的形函数具有点插值性和高阶光滑性,不仅能够直接施加本质边界条件,而且能保证较高的计算精度.为了更有效地求解三维轴对称弹性动力学问题,对重构核插值法(reproducing kernel interpolation method, RKIM)应用于此类问题进行了研究,并发展了相应的数值模拟方法.由于几何形状和边界条件的轴对称性,计算时只需要横截面上离散节点的信息,因而前处理变得简单.采用Newmark-β法进行了时域积分.数值算例表明,轴对称弹性动力学分析的重构核插值法既有无网格方法的优势,又有较高的计算精度.     

非定常扩散方程的单元中心型有限体积格式

本文基于调和平均点建立了一种新的单元中心型有限体积格式,用以求解非定常扩散方程.在网格边上离散法向流时,选择该网格边两端点和该边上的一个调和平均点作为辅助插值点,并将这些辅助插值点上的未知量用网格单元中心点的未知量进行替换,最终得到一个只含网格单元中心未知量的有限体积格式.该格式满足线性精确性质和局部守恒性,且适用于任意多边形网格.在六种不同的多边形网格上进行四个数值实验,分别考虑扩散系数是连续的和间断的以及非线性的情况,数值结果表明:本文所构造的格式在六种网格上的L2误差均可达到二阶收敛精度,对于不同类型的扩散系数,该格式保持良好的鲁棒性,并且从编程实现的角度来说,该格式更易于向三维情况推广.     

A meshless method based on boundary integral equations and radial basis functions for biharmonic-type problems

This paper presents a meshless method, which replaces the inhomogeneous biharmonic equation by two Poisson equations in terms of an intermediate function. The solution of the Poisson equation with the intermediate function as the right-hand term may be written as a sum of a particular solution and a homogeneous solution of a Laplace equation. The intermediate function is approximated by a series of radial basis functions. Then the particular solution is obtained via employing Kansa’s method, while the homogeneous solution is approximated by using the boundary radial point interpolation method by means of boundary integral equations. Besides, the proposed meshless method, in conjunction with the analog equation method, is further developed for solving generalized biharmonic-type problems. Some numerical tests illustrate the efficiency of the method proposed.     

Dual reciprocity hybrid radial boundary node method for the analysis of Kirchhoff plates

A meshless method of dual reciprocity hybrid radial boundary node method (DHRBNM) for the analysis of arbitrary Kirchhoff plates is presented, which combines the advantageous properties of meshless method, radial point interpolation method (RPIM) and BEM. The solution in present method comprises two parts, i.e., the complementary solution and the particular solution. The complementary solution is solved by hybrid radial boundary node method (HRBNM), in which a three-field interpolation scheme is employed, and the boundary variables are approximated by RPIM, which is applied instead of moving least square (MLS) and obtains the Kronecker’s delta property where the traditional HBNM does not satisfy. The internal variables are interpolated by two groups of symmetric fundamental solutions. Based on those, a hybrid displacement variational principle for Kirchhoff plates is developed, and a meshless method of HRBNM for solving biharmonic problems is obtained, by which the complementary solution can be solved.     

关于薄板的无网格局部边界积分方程方法中的友解

无网格局部边界积分方程方法是最近发展起来的一种新的数值方法,这种方法综合了伽辽金有限元、边界元和无单元伽辽金法的优点,是一种具有广阔应用前景的、真正的无网格方法．把无网格局部边界积分方程方法应用于求解薄板问题,给出了薄板无网格局部边界积分方程方法所需要的友解及其全部公式．     

薄板的局部Petrov-Galerkin方法

利用薄板控制微分方程的等效积分对称弱形式和对变量(挠度)采用移动最小二乘近似函数进行插值,研究了薄板弯曲问题的无网格局部Petrov-Galerkin方法．这是一种真正的无网格方法,它不需要任何有限元或边界元网格,不管这种网格是用于能量积分还是进行插值的目的．所有的积分都在规则形状的子域及其边界上进行,并用罚因子法施加本质边界条件．数值例子表明,无网格局部Petrov-Galerkin法不但能够求解二阶微分方程的边值问题,而且求解四阶微分方程的边值问题也很有效,也具有收敛快、稳定性好、对挠度和内力都具有精度高的特点．     

Meshless method with ridge basis functions

Meshless collocation methods for the numerical solutions of PDEs are increasingly adopted due to their advantages including efficiency and flexibility, and radial basis functions are popularly employed to represent the solutions of PDEs. Motivated by the advantages of ridge basis function representation of a given function, such as the connection to neural network, fast convergence as the number of terms is increased, better approximation effects and various applications in engineering problems, a meshless method is developed based on the collocation method and ridge basis function interpolation. This method is a truly meshless technique without mesh discretization: it neither needs the computation of integrals, nor requires a partition of the region and its boundary. Moreover, the method is applied to elliptic equations to examine its appropriateness, numerical results are compared to that obtained from other (meshless) methods, and influence factors of accuracy for numerical solutions are analyzed.     

The closed-form particular solutions for Laplace and biharmonic operators using a Gaussian function

Particular solutions play a critical role in solving inhomogeneous problems using boundary methods such as boundary element methods or boundary meshless methods. In this short article, we derive the closed-form particular solutions for the Laplace and biharmonic operators using the Gaussian radial basis function. The derived particular solutions are implemented numerically to solve boundary value problems using the method of particular solutions and the localized method of approximate particular solutions. Two examples in 2D and 3D are given to show the effectiveness of the derived particular solutions.     

Spectral meshless radial point interpolation (SMRPI) method to two‐dimensional fractional telegraph equation

H. Ammari In this article, an innovative technique so‐called spectral meshless radial point interpolation (SMRPI) method is proposed and, as a test problem, is applied to a classical type of two‐dimensional time‐fractional telegraph equation defined by Caputo sense for (1 < α≤2). This new methods is based on meshless methods and benefits from spectral collocation ideas, but it does not belong to traditional meshless collocation methods. The point interpolation method with the help of radial basis functions is used to construct shape functions, which play as basis functions in the frame of SMRPI method. These basis functions have Kronecker delta function property. Evaluation of high‐order derivatives is not difficult by constructing operational matrices. In SMRPI method, it does not require any kind of integration locally or globally over small quadrature domains, which is essential of the finite element method (FEM) and those meshless methods based on Galerkin weak form. Also, it is not needed to determine strict value for the shape parameter, which plays an important role in collocation method based on the radial basis functions (Kansa's method). Therefore, computational costs of SMRPI method are less expensive. Two numerical examples are presented to show that SMRPI method has reliable rates of convergence. Copyright © 2015 John Wiley & Sons, Ltd.     

Scattered data fitting of Hermite type by a weighted meshless method

In many practical problems, it is often desirable to interpolate not only the function values but also the values of derivatives up to certain order, as in the Hermite interpolation. The Hermite interpolation method by radial basis functions is used widely for solving scattered Hermite data approximation problems. However, sometimes it makes more sense to approximate the solution by a least squares fit. This is particularly true when the data are contaminated with noise. In this paper, a weighted meshless method is presented to solve least squares problems with noise. The weighted meshless method by Gaussian radial basis functions is proposed to fit scattered Hermite data with noise in certain local regions of the problem’s domain. Existence and uniqueness of the solution is proved. This approach has one parameter which can adjust the accuracy according to the size of the noise. Another advantage of the weighted meshless method is that it can be used for problems in high dimensions with nonregular domains. The numerical experiments show that our weighted meshless method has better performance than the traditional least squares method in the case of noisy Hermite data.     

On the analysis of a kind of nonlinear Sobolev equation through locally applied pseudo‐spectral meshfree radial point interpolation

In this study, we develop an approximate formulation for two‐dimensional nonlinear Sobolev problems by focusing on pseudospectral meshless radial point interpolation (PSMRPI) which is a kind of locally applied radial basis function interpolation and truthfully a meshless approach. In the PSMRPI method, the nodal points do not need to be regularly distributed and can even be quite arbitrary. It is easy to have high order derivatives of unknowns in terms of the values at nodal points by constructing operational matrices. The convergence and stability of the technique in some sense are studied via some examples to show the validity and trustworthiness of the PSMRPI technique.     

Deformation of a heavy viscoelastic cylinder supported on a rigid plate

Plane strain of a viscoelastic cylinder supported arcwise on a rigid plate is considered. By means of Laplace transformation the viscoelastic problem is reduced to the perfectly elastic problem. The elastic solution is found by methods of complex variables by means of the stress function. Inverse transformation is carried out by the method of Il'yushin approximations. the first approximation of the stated problems is found in the method of approximations.