利用边界元法求解一类重调和方程

边界元法(BEM)和多重互易法(MRM)相结合求解一类重调和方程.通过重调和基本解序列给出的MRM-方法和BEM, 推导出该类问题的MRM-边界变分方程, 用边界元法求解该变分方程, 从而得到重调和方程的近似解, 并给出了解的存在唯一性证明.通过数值算例说明了MRM-方法具有收敛速度快、计算精度高, 易编程等优点, 为使用边界元法数值求解重调和方程提供了方法和理论依据.适合于工程中的实际运算.     

带有摩擦的单边接触大变形问题的研究(Ⅰ)——增量变分方程

本文以非线性连续体几何场论为基本理论和方法,建立了拖带坐标下弹塑性大变形增量变分方程的更一般表示式.给出了二维、三维连续体接触边界变化率公式,得到了变边界接触大变形增量变分公式和速率型变分不等式,为有限元计算求解带有摩擦弹塑性大变形接触问题提供了理论基础.     

第二类混合变分不等式的MRM-边界元分析

本文以弹性力学中的摩擦问题为背景,采用多重互易方法(MRM方法),边界元方法,将摩擦问题中的第二类混合变分不等式化解为MRM-边界混合变分不等式,给出了MRM-边界混合变分不等式解的存在唯—性,通过引入变换将原MRM-边界混合变分不等式化解为标准的凸极值问题,采用正则化方法处理后,给出了MRM-边界混合变分不等式的迭代分解方法。文末给出了数值算例。     

有限元广义伽略金方程,边界变分方程,边界积分方程

在[1]的基础上,我们进一步应用可动边界的变分原理于固体体系的离散分析,得到有限元广义伽略金方程,边界变分方程,边界积分方程.这些方程描述了待解函数在元素内部与元素的边界上应满足的方程.当对固体体系进行离散分析时,可以应用这些方程去建立不同情况下的求解待解函数的离散方程.亦可作为相应情况下的简化计算的依据.由本文得到的边界积分方程可知,在[2]中提出的J积分形式,应用于内部元素边界的围道积分计算是不适宜的.     

偏微分方程Δ2u-sΔu+k2u=0边值问题的MRM边界变分方程

导出边值问题Δ2u-sΔu+k2u=o;x∈Ω∪Ω'(R2;u|г=uo;аu/аn|г=go的定解问题,MRM边界变分方程,全平面解的表达式.从中可以看出,MRM边界变分方程中只包含弱奇异积分核,并且自动消除了原第一、二MRM边界积分方程中出现的强奇异积分核.问题解的表达式后并不加任何多项式,因而也不需要引入Lagrange乘子求解该项,这给边界元数值求解过程带来极大的方便.数值分析结果表明该方法具有明显优势.     

带惩罚的有限元方法的收敛性

在用有限元方法解二阶椭圆型方程的边值问题时,首先将边值问题化为一个等价的变分问题,即泛函的极值问题。对于第二,第三边值问题,在相应的变分问题中边界条件被吸收到泛函的表示式中。因而在求泛函极值时不再对允许函数类附加边界条件的约束。在这种情况下边界条件就称为自然边界条件,相应的变分问题称为无约束变分问题。而对于第一边值问题即所谓狄氏问题则不然,相应的等价变分问题是带约束的变分问题。求泛函极值时的允许函数类必须满足强加的边界条件。例如我们考虑二维有界区域Ω上的方程     

含开边界二维Stokes问题的Galerkin边界元解法

本文推导了含有开边界的二维有限域上Stokes问题的边界积分方程, 得出基于单层位势的第一类间接边界积分方程.对与之等价的边界变分方程用Galerkin边界元求解以得出单层位势的向量密度. 对于含有开边界端点的边界单元,采用特别的插值函数, 以模拟其固有的奇异性.论文用若干数值算例模拟了含有开边界的有限区域上不可压缩粘性流体的绕流.       

变系数导热方程的Robin系数反演问题

本文研究基于非局部边界附加条件下,一类变系数抛物型方程的Robin系数确定问题,这里的Robin系数仅与时间相关.首先给出了变分公式,并利用变分公式证明了解的唯一性,其次给出了时间离散模型,基于线性离散化的变分形式,导出了一系列先验估计,证明了弱解的存在性,并对其进行了误差分析.     

论微分与积分方程以及有限与无限元

1 椭圆微分方程的边界值问题可以有种种不同的数学成型,在理论上等价,但在实践上不等效。有限元方法成功的一个关键就是合理选取了变分的数学型式。举例来说,取调和方程的第二类边界问题,定义于区域Ω,具有光滑边界Г:     

椭圆外区域上Helmholtz问题的耦合法

以椭圆外区域上Helmholtz方程为例,研究一种带有椭圆人工边界的自然边界元与有限元耦合法,给出了耦合变分问题的适定性及误差分析并给出数值例子.理论分析及数值结果表明,用方法求解椭圆外问题是十分有效的.为求解具有长条型内边界外Helmholtz问题提供了一种很好的数值方法.     

The equilibrium boundary element method in boundary-value problems of elasticity theory

An equilibrium boundary element method is proposed for solving boundary-value problems in the theory of elasticity, thermo-elasticity, the dynamical theory of elasticity, bar torsion calculations, and the bending of a plate. The idea is to use simultaneously the method of constructing bundles of functions which exactly satisfy the equilibrium equations, the boundary variational equations of mechanics, and the methods of discrete finite-element approximation. The variational method of constructing the resolving boundary equations ensures that the linear system is symmetric and easily coupled to the finite-element method. Since volume integrals are eliminated the dimensions of the problem are reduced by one, but, unlike the boundary element method, there is no need to know the fundamental solutions. The solution of some bar torsion and plate bending problems confirms the high numerical efficiency of the method.     

The h-p version of the coupling of finite element and boundary element methods for transmission problems in polyhedral domains

Summary. This paper analyzes the rate of convergence of the h-p version of the coupling of the finite element and boundary element method for transmission problems with a linear differential operator with variable coefficients in a bounded polyhedral domain and with constant coefficients in the exterior domain . This procedure uses the variational formulation of the differential equation in and involves integral operators on the interface between and . The finite elements are used to obtain approximate solutions of the differential equation in and the boundary elements are used to obtain approximate solutions of the integral equations. For given piecewise analytic data we show that the Galerkin solution of this coupling procedure converges exponentially fast in the energy norm if the h-p version is used both for finite elements and boundary elements. Received February 10, 1996 / Revised version received April 4, 1997     

On the mixed finite element method with Lagrange multipliers

In this note we analyze a modified mixed finite element method for second‐order elliptic equations in divergence form. As a model we consider the Poisson problem with mixed boundary conditions in a polygonal domain of R ². The Neumann (essential) condition is imposed here in a weak sense, which yields the introduction of a Lagrange multiplier given by the trace of the solution on the corresponding boundary. This approach allows to handle nonhomogeneous Neumann boundary conditions, theoretically and computationally, in an alternative and usually easier way. Then we utilize the classical Babu?ka‐Brezzi theory to show that the resulting mixed variational formulation is well posed. In addition, we use Raviart‐Thomas spaces to define the associated finite element method and, applying some elliptic regularity results, we prove the stability, unique solvability, and convergence of this discrete scheme, under appropriate assumptions on the mesh sizes. Finally, we provide numerical results illustrating the performance of the algorithm for smooth and singular problems. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 192–210, 2003     

Construction of locally conservative fluxes for the SUPG method

We consider the construction of locally conservative fluxes by means of a simple postprocessing technique obtained from the finite element solutions of advection diffusion equations. It is known that a naive calculation of fluxes from these solutions yields nonconservative fluxes. We consider two finite element methods: the usual continuous Galerkin finite element method for solving nondominating advection diffusion equations and the streamline upwind/Petrov‐Galerkin method for solving advection dominated problems. We then describe the postprocessing technique for constructing conservative fluxes from the numerical solutions of the general variational formulation. The postprocessing technique requires solving an auxiliary Neumann boundary value problem on each element independently and it produces a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. We provide a convergence analysis for the postprocessing technique. Performance of the technique and the convergence behavior are demonstrated through numerical examples including a set of test problems for advection diffusion equations, advection dominated equations, and drift‐diffusion equations. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1971–1994, 2015     

Multiplicity of solutions for nonlinear second order impulsive differential equations with linear derivative dependence via variational methods

This paper uses critical point theory and variational methods to investigate the multiple solutions of boundary value problems for second order impulsive differential equations. The conditions for the existence of multiple solutions are established. An example is constructed to illustrate the proposed result.     

BOUNDARY ELEMENT ANALYSIS FOR A CLASS OF MULTI-DIMENSIONAL LINEAR PARABOLIC EQUATIONS

In this paper,by me as of beundary element method,we try to deal with the initial -boundary value problem for a class of linear parunolic equations,which is a linear heat conduction equation. We tresent a boundary integral equation for the solution to the problem and its variational formalation The well-posedness of the variational formulation is proved. And the error estimates for the approsutate solutions are provided. The results of this paper are more general than those of[1]     

Energetic BEM-FEM coupling for the numerical solution of the damped wave equation

Time-dependent problems modeled by hyperbolic partial differential equations can be reformulated in terms of boundary integral equations and solved via the boundary element method. In this context, the analysis of damping phenomena that occur in many physics and engineering problems is a novelty. Starting from a recently developed energetic space-time weak formulation for the coupling of boundary integral equations and hyperbolic partial differential equations related to wave propagation problems, we consider here an extension for the damped wave equation in layered media. A coupling algorithm is presented, which allows a flexible use of finite element method and boundary element method as local discretization techniques. Stability and convergence, proved by energy arguments, are crucial in guaranteeing accurate solutions for simulations on large time intervals. Several numerical benchmarks, whose numerical results confirm theoretical ones, are illustrated and discussed.     

On the p-version approximation in the boundary element method for a variational inequality of the second kind modelling unilateral contact and given friction

In this paper we complement recent work of Maischak and Stephan on adaptive hp-versions of the BEM for unilateral Signorini problems, respectively on FEM-BEM coupling in its h-version for a nonlinear transmission problem modelling Coulomb friction contact. Here we focus on the boundary element method in its p-version to treat a scalar variational inequality of the second kind that models unilateral contact and Coulomb friction in elasticity together. This leads to a nonconforming discretization scheme. In contrast to the work cited above and to a related paper of Guediri on a boundary variational inequality of the second kind modelling friction we take the quadrature error of the friction functional into account of the error analysis. At first without any regularity assumptions, we prove convergence of the BEM Galerkin approximation in the energy norm. Then under mild regularity assumptions, we establish an a priori error estimate that is based on a novel Céa–Falk lemma for abstract variational inequalities of the second kind.     

三维Helmholtz方程边值问题的新型边界积分-微分方程及其边界元解法

用双层位势表示的二维Neumann边值问题的边界归化方法,将原始问题归化为新型边界积分-微分方程,由此导出一种新的既能保持原始问题的自伴性,又具有可积弱奇性积分核的边界变分方程.本文将此法推广到三维Helmholtz方程Neumann边值问题,并给出最优能量模误差估计和内部最大模超收敛估计.     

Variational approach for Sturm-Liouville boundary value problems on time scales

In this paper, we consider a class of second order Sturm-Liouville boundary value problems with positive parameter λ on time scales. By using variational method and critical point theory, we obtain that the boundary value problem has solutions for λ being in some different intervals. Recent results in the literature are generalized and improved.