位势问题边界元法中几乎奇异积分的正则化

将一种通用算法应用于平面位势问题边界元法中近边界点几乎奇异积分的正则化。对线性单元,位势问题近边界点的几乎强和超奇异积分可归纳为两种形式。通过分部积分,将引起奇异的积分元素变换到积分号之外,从而对这两种积分分别给出了无奇异的正则化计算公式。除了线性元,二次元也应用于该算法。与近边界点临近的二次单元划分为两段线性单元,该算法仍然适用。算例证明了方法的有效性和精确性。对曲线边界问题,联合二次元和线性元可提高计算结果精确度。     

刚性基础上的弹性层在表面垂直集中简谐载荷作用下的动力响应

本文首次将文献[1]所提出的线载荷积分方程法应用于求解弹性动力学问题。导出了刚性基础上的弹性层在表面垂直集中简谐载荷作用下动力响应问题的一维非奇异积分方程组,并求得了数值解。     

圆内平面弹性问题的边界积分公式

根据双解析函数可以得到单位圆内平面弹性问题应力函数的边界积分公式,但式中包含强奇异积分,不能用于直接计算．将边界上的应力函数展开为Fourier级数,再利用广义函数论中的几个公式进行卷积计算,可以得到不含强奇异积分核的边界积分公式,通过边界的应力函数值和法向导数的积分,直接得到圆内应力函数值,并给出几个算例,表明该结果用于求解单位圆内平面弹性问题十分方便．     

Hermite三次样条多小波自然边界元法

针对应用自然边界元方法解上半平面的Laplace方程的Neumann边值问题时存在奇异积分的困难,本文提出了Hermite三次样条多小波自然边界元法.Hermite三次样条多小波具有较短的紧支集、很好的稳定性和显式表达式,而且它们在不同层上的导数还是相互正交的.因此,本文将它与自然边界元法相结合,利用小波伽辽金法离散自然边界积分方程,使自然边界积分方程中的强奇异积分化为弱奇异积分,从而降低了问题的复杂性.文中给出的算例表明了该方法的可行性和有效性.     

改进的无奇异局部边界积分方程方法

在局部边界积分方程方法中,当源节点位于分析域的整体边界上时,局部边界积分将出现奇异积分问题,这些奇异积分需要做特别的处理．为此,提出了对域内节点采用局部积分方程,而对边界节点直接采用移动最小二乘近似函数引入边界条件来解决奇异积分问题,这同时也解决了对积分边界进行插值引入近似误差的问题．作为应用和数值实验,对Laplace方程和Helmholtz方程问题进行了分析,取得了很好的数值结果．进而,在Helmholtz方程求解中,采用了含波解信息的修正基函数来代替单项式基函数进行近似．数值结果显示,这样处理是简单高效的,在高波数声传播问题的求解中非常具有前景．     

工程中一类拼接问题的复变方法

讨论工程中一类含边裂缝弹性材料补强的拼接问题 .根据平面弹性复变方法 ,问题归结为一类解析函数的边值问题 ,通过有效的分析方法和积分变换 ,进一步将问题简化为一类奇异积分方程 ,证明了方程解的存在唯一 ,并对方程解的简化进行了研究 ,得到了弹性材料体内应力分布的封闭形式解 ,并导出一直裂缝情况裂缝尖端应力强度因子的表达式     

双周期裂纹场平面弹性焊接的数学问题

本文讨论双周期胞腔中含任意形状裂纹的不同材料的弹性平面焊接(焊线为任意形状的封闭光滑曲线)的第二基本问题.运用Мусхелишвили复变函数方法,对这类弹性平面问题建立起了数学模型,将求解弹性平衡问题化归为寻求复应力函数满足一定边界条件的边值问题,然后构造其解的形式,再将其转化为正则型的奇异积分方程,数学上严格证明其解的存在与唯一.     

固体弹性三维问题统一解

依据三维弹性力学问题的Kelvin解,用三维虚边界元法来建立积分方程,从而使三维实体和各类板、壳等问题的求解思想得到统一.对各类三维问题采用统一的思想建立数学模型,更有利于程序模块化,增强了程序的通用性.另外,建立积分方程时直接引用Kelvin解,而未引入任何其它假设,使该方法的解更偏于实际,且使应用范围拓宽.再者,与边界元直接法相比,该方法的优点在于无需处理奇异积分,且系数阵是对称的.最后给出部分算例,以证明方法的有效性和计算精度.     

一类具有弱奇异核的偏积分微分方程的Chebyshev小波数值方法（英文）

本文提出一种基于第四类Chebyshev小波配置法,求解了一类具有弱奇异核的偏积分微分方程数值解.利用第四类移位Chebyshev多项式,在Riemann-Liouville分数阶积分意义下,导出Chebyshev的分数次积分公式.通过利用分数次积分公式和二维的第四类Chebyshev小波结合配置法,将具有弱奇异核的偏积分微分方程转化为代数方程组求解.给出了第四类Chebyshev小波的收敛性分析.数值例子证明了本文方法的有效性.     

嵌在弹性半空间的刚性变直径圆轴的扭转*

本文用线载荷积分方程法(LLIEM)研究嵌在弹性半空间的刚性变直径圆轴的轴对称扭转.利用将“半空间的点环力偶”(PRCHS)这一虚的基本载荷沿对称轴的圆轴区间中分布就能使本问题归结为一个一维的非奇异的Ferdholm第一种积分方程,且很易用数值求解.文中给出刚性圆锥,圆柱,圆锥柱嵌在弹性半空间的扭转的数值计算的例并和他人的已知结果相比较.文中也给出了刚性半球嵌在弹性半空间的扭转的精确解.     

奇异杂交边界点方法求解泊松方程

As a boundary-type meshless method,the singular hybrid boundary node method(SHBNM)is based on the modified variational principle and the moving least square(MLS)approximation,so it has the advantages of both boundary element method(BEM)and meshless method.In this paper,the dual reciprocity method(DRM)is combined with SHBNM to solve Poisson equation in which the solution is divided into particular solution and general solution.The general solution is achieved by means of SHBNM,and the particular solution is approximated by using the radial basis function(RBF).Only randomly distributed nodes on the bounding surface of the domain are required and it doesn't need extra equations to compute internal parameters in the domain.The postprocess is very simple.Numerical examples for the solution of Poisson equation show that high convergence rates and high accuracy with a small node number are achievable.     

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.     

A new model for the analysis of reinforced concrete members with a coupled HdBNM/FEM

As a truly boundary-type meshless method, the hybrid boundary node method (HdBNM) does not require ‘boundary element mesh’, either for the purpose of interpolation of the solution variables or for the integration of ‘energy’. In this paper, the HdBNM is coupled with the finite element method (FEM) for predicting the mechanical behaviors of reinforced concrete. The steel bars are considered as body forces in the concrete. A bond model is presented to simulate the bond-slip between the concrete and steels using fictitious spring elements. The computational scale and cost for meshing can be further reduced. Numerical examples, in 2D and 3D cases, demonstrate the efficiency of the proposed approach.     

含裂纹平面问题Erdogan基本解的显式表达

基本解是边界元法、基本解法和无网格法等数值方法的重要理论基础.在断裂问题中,采用含裂纹的基本解可以避免将裂纹表面作为边界条件,从而大大简化问题的求解.在复变函数表示的含裂纹平面问题Erdogan基本解的基础上,对Erdogan基本解的使用条件进行了注解,修正了Erdogan基本解的一些错误,并推导出Erdogan基本解中位移函数解答的显式表达形式.编写了基于Erdogan基本解显式表达的样条虚边界元法（spline fictitious boundary element method, SFBEM）计算程序,计算了具有复合边界条件平面问题的位移、应力和应力强度因子.数值算例结果表明了该文提出的Erdogan基本解显式表达形式的正确性.     

New method based on the HPM and RKHSM for solving forced Duffing equations with integral boundary conditions

This paper investigates the forced Duffing equation with integral boundary conditions. Its approximate solution is developed by combining the homotopy perturbation method (HPM) and the reproducing kernel Hilbert space method (RKHSM). HPM is based on the use of the traditional perturbation method and the homotopy technique. The HPM can reduce nonlinear problems to some linear problems and generate a rapid convergent series solution in most cases. RKHSM is also an analytical technique, which can solve powerfully linear boundary value problems. Therefore, the forced Duffing equation with integral boundary conditions can be solved using advantages of these two methods. Two numerical examples are presented to illustrate the strength of the method.     

Solving boundary value problems,integral, and integro-differential equations using Gegenbauer integration matrices

We introduce a hybrid Gegenbauer (ultraspherical) integration method (HGIM) for solving boundary value problems (BVPs), integral and integro-differential equations. The proposed approach recasts the original problems into their integral formulations, which are then discretized into linear systems of algebraic equations using Gegenbauer integration matrices (GIMs). The resulting linear systems are well-conditioned and can be easily solved using standard linear system solvers. A study on the error bounds of the proposed method is presented, and the spectral convergence is proven for two-point BVPs (TPBVPs). Comparisons with other competitive methods in the recent literature are included. The proposed method results in an efficient algorithm, and spectral accuracy is verified using eight test examples addressing the aforementioned classes of problems. The proposed method can be applied on a broad range of mathematical problems while producing highly accurate results. The developed numerical scheme provides a viable alternative to other solution methods when high-order approximations are required using only a relatively small number of solution nodes.     

A hybrid method using wavelets for the numerical solution of boundary value problems on the interval

In this work, various aspects of wavelet-based methods for second order boundary value problems under Galerkin framework are investigated. Based on the B-spline multiresolution analysis (MRA) on the line we propose a hybrid method on the interval which combines different treatments for interior and boundary splines. By using this procedure, the MRA structure was conserved and hierarchical representations of the solution at different scales were obtained without much computational effort. Numerical examples are given to verify the effectiveness of the proposed method and the comparison with other techniques is presented.     

薄板的局部Petrov-Galerkin方法

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

Application of reproducing kernel method to third order three-point boundary value problems

This paper investigates the analytical approximate solutions of third order three-point boundary value problems using reproducing kernel method. The solution obtained by using the method takes the form of a convergent series with easily computable components. However, the reproducing kernel method can not be used directly to solve third order three-point boundary value problems, since there is no method of obtaining reproducing kernel satisfying three-point boundary conditions. This paper presents a method for solving reproducing kernel satisfying three-point boundary conditions so that reproducing kernel method can be used to solve third order three-point boundary value problems. Results of numerical examples demonstrate that the method is quite accurate and efficient for singular second order three-point boundary value problems.     

The dual reciprocity boundary integral equation technique to solve a class of the linear and nonlinear fractional partial differential equations

In this paper, we apply the boundary integral equation technique and the dual reciprocity boundary elements method (DRBEM) for the numerical solution of linear and nonlinear time‐fractional partial differential equations (TFPDEs). The main aim of the present paper is to examine the applicability and efficiency of DRBEM for solving TFPDEs. We employ the time‐stepping scheme to approximate the time derivative, and the method of linear radial basis functions is also used in the DRBEM technique. This method is improved by using a predictor–corrector scheme to overcome the nonlinearity that appears in the nonlinear problems under consideration. To confirm the accuracy of the new approach, several examples are presented. The convergence of the DRBEM is studied numerically by comparing the exact solutions of the problems under investigation. Copyright © 2015 John Wiley & Sons, Ltd.