边界元方法的一些研究进展

本文旨在综述我们小组近二十年来在边界元方法这一领域的一些研究成果,在简要介绍边界元方法的基本思想后,主要介绍了一类非线性界面问题的有限元-边界元耦合方法、求解电磁散射问题的有限元-边界元耦合方法和超奇异积分的一类计算方法.     

“混合元”的有限元法

大家知道,用有限元求解微分方程时,如果微分方程的解足够光滑,边界比较接近拆线,采用高阶元比采用低阶元网格可以稀得多,因此这时采用高阶元比低阶元有利得多,但在实际计算中,边界常常比较复杂,解在边界附近变化较大,这时只采用高阶元在边界     

“混合元”的有限元法

大家知道,用有限元求解微分方程时,如果微分方程的解足够光滑,边界比较接近拆线,采用高阶元比采用低阶元网格可以稀得多,因此这时采用高阶元比低阶元有利得多,但在实际计算中,边界常常比较复杂,解在边界附近变化较大,这时只采用高阶元在边界     

边界元方法的数学分析

本文介绍边界元方法与其它数值计算方法的关系.用拟微分算子给边界元法中所遇到的各种类型的边界积分方程以统一的数学描述.由拟微分算子的强椭圆性,得出边界积分方程的可解性.本文还介绍边界元空间的建立和几种求解边界积分方程的离散化方法.对于边界积分方程的变分形式,给出边界元近似解的收敛性和渐近误差估计.     

非齐次弹性力学问题双互易边界元方法研究北大核心CSCD

基于弹性力学边界元方法理论,将边界元法与双互易法结合,采用指数型基函数对非齐次项进行插值得到双互易边界积分方程.将边界积分方程离散为代数方程组,利用已知边界条件和方程特解求解方程组,得出域内位移和边界面力.指数型基函数的形状参数是由插值点最近距离的最小值决定,采用这种形状参数变化方案,分析径向基函数(RBF)插值精度以及插值稳定性.再次将指数型基函数应用到双互易边界元法中,分析双互易边界元方法下计算精度及稳定性,验证了指数型插值函数作为双互易边界元方法的径向基函数解决弹性力学域内体力项问题的有效性.     

不同边界条件对路基结构动力学响应的影响

分别采用静力人工边界、无限元边界、粘弹性动力人工边界和远置静力人工边界对承受交通载荷的典型路基结构进行有限元分析。以远置静力人工边界计算获得的路面竖向位移和面层底部拉应力时程曲线为参考解,揭示边界条件对路基结构动力学响应的影响。分析结果表明:边界条件对拉应力的影响很小,它们与参考解的最大误差为0.82%;静力人工边界下的位移时程曲线一直在参考解附近震荡,无限元边界下的位移解高于参考解,粘弹性动力人工边界下的路面位移与参考解最为接近。综合考虑位移和应力的精确度与CPU的计算时间,粘弹性动力人工边界为最佳的路基结构动力学分析边界条件。     

求解变系数非齐次亥姆霍茨方程的边界单元法

本文利用拉普拉斯方程的基本解作为权函数，给出求解交系数非齐次亥姆霍茨方程的迭代格式，进而得到求解这类方程的边界元迭代法．文中给出了算例．最后，把本文给出的边界元迭代法与作者早些时候提出的边界元耦合法进行了比较．     

关于弹塑性边界元分析的迭代收敛性

迭代过程是弹塑性边界元法的核心部分,文中对边界元弹塑性分析的迭代收敛性问题作了详细地理论分析和讨论.     

曲边有限元和自然边界元求解外边值问题的耦合法

1　引　　言边界元与有限元耦合法在科学和工程计算中有着独特的作用 .由于区域的无限性给人们常用的有限元方法带来困难 ,边界元方法又难以独立处理非线性和非均质的问题以及具有不规则边界的区域上的问题 ,而两者相结合却可以克服各自的缺点 ,故边界元与有限元耦合法在处理一般区域问题特别是无界区域问题时便得到科学与工程界的青睐 ,获得了比较广泛的应用 .自然边界元方法并不引入新的变量 ,属于直接边界元方法[2 ] [8] .它保持能量不变和原边值问题的许多有用性质 ,例如双线性型的对称性和强制性 ,从而自然积分方程的解的存在唯一性及…     

抛物型初边值问题的有限元与边界积分耦合的离散化及其误差分析

1．引言边界元方法是近二十几年来迅速发展起来的一类新的偏微分方程的数值方法．它的独特之处是将空间的维数降低一维,从而倍受工程技术人员的青睐,并在工程技术与计算数学领域得到越来越广泛的重视和研究．对椭圆型问题,边界元方法的理论与应用研究已取得丰硕成果;对发展型问题,近年来在理论方面的研究也已取得重要进展［6－11］．但边界元方法难以处理非均质问题,而有限元对各类问题及各种区域具有较好的适应性,将两者结合起来可充分发挥各自的优点．文山提出了一种抛物方程初边值问题的有限元与边界积分的耦合方法,其主要思想是…     

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.     

A numerical algorithm for solving a class of linear nonlocal boundary value problems

In order to solve a class of linear nonlocal boundary value problems, a new reproducing kernel space satisfying nonlocal conditions is constructed carefully. This makes it easy to solve the problems. Furthermore, the exact solutions of the problems can be expressed in series form. The numerical results demonstrate that the new method is quite accurate and efficient for solving fourth-order nonlocal boundary value problems.     

A new error estimation approach based on fuzzy logic system for H-R adaptive boundary element method

Conventional adaptive boundary element method cannot be universally applied to solve many more problems than the subject it discussed, and different error estimation formulas need to be designed for varied problems. This paper put forward a new error analysis method based on the fuzzy logic system, which is able to make error estimation effectively using human expert experience, and solve the two classical elasticity problems in conjunction with the H-R adaptive boundary element method. Numerical examples have illustrated the effectiveness, superiority and potential of a fuzzy logic approach in the adaptive boundary element method.     

An asymptotic initial value method for second order singular perturbation problems of convection-diffusion type with a discontinuous source term

In this paper a numerical method is presented to solve singularly perturbed two points boundary value problems for second order ordinary differential equations consisting a discontinuous source term. First, in this method, an asymptotic expansion approximation of the solution of the boundary value problem is constructed using the basic ideas of a well known perturbation method WKB. Then some initial value problems and terminal value problems are constructed such that their solutions are the terms of this asymptotic expansion. These initial value problems are happened to be singularly perturbed problems and therefore fitted mesh method (Shishkin mesh) are used to solve these problems. Necessary error estimates are derived and examples provided to illustrate the method.     

Singular boundary integral equations of boundary value problems of elastic dynamics in the case of subsonic running loads

We consider boundary value problems for an elastic medium bounded by a cylindrical surface on which a load with a constant-in-time form is moving at a constant subsonic velocity. This class of problems is a model class for the dynamics of underground structures like transport tunnels and ground transport. The method of generalized functions is developed for the solution of boundary value problems. We construct dynamic analogs of Green formulas and use them to derive singular boundary integral equations, which solve the considered boundary value problems.     

二阶线性齐次微分方程边值问题相似构造解的应用及Matlab图版分析

在二阶线性齐次微分方程边值问题相似构造解式的基础上,首先利用相似构造法求解Bessel方程和变型的Bessel方程边值问题的解,然后建立了均质油藏的渗流规律的数学模型,再将均质油藏的渗流数学模型转换成变型的Bessel方程的边值问题,利用二阶线性齐次微分方程边值问题的相似构造法求解均质储层渗流的数学模型.最后通过Matlab编程进行图版分析,展示实例的函数解.这将极大地方便试进分析软件的编制,也提高了石油工作者的效率.     

带边界层的线性两点边值问题的打靶——小波配点法

本文将打靶法和小波配点法相结合，提出了打靶－小波配点数值算法，用于求解带边界层的常微分方程边值问题。文中给出了数值算例，并进行了分析，验证了这种方法对处理边界层问题的有效性。     

A boundary element approach for solving plane elastostatic equations of anisotropic functionally graded materials

A boundary element method is proposed for the numerical solution of an important class of boundary value problems governed by plane elastostatic equations of anisotropic functionally graded materials. The grading function of the material properties may be any general function that varies smoothly from point to point in the material. The proposed boundary element method is applied to solve some specific problems to check its validity and accuracy.     

基于再生核三次样条基底的构造及其应用

将三次样条理论与再生核理论相结合,利用再生核函数巧妙地构造了三次样条函数空间的一组基底.基于三次样条插值的高收敛特点,得到了微分方程边值问题近似解的一种新的求解方法.数值算例展现出算法简单、有效.     

Complete second order differential equations in Banach spaces with dynamic boundary conditions

In this paper, we deal with the mixed initial boundary value problem for complete second order (in time) linear differential equations in Banach spaces, in which time-derivatives occur in the boundary conditions. General wellposedness theorems are obtained (for the first time), which are used to solve the corresponding inhomogeneous problems. Examples of applications to initial boundary value problems for partial differential equations are also presented.