带裂纹方形截面杆扭转问题的自然边界元与有限元耦合法

根据基于区域分解的自然边界元与有限元的耦合法,研究了带裂纹的方形截面杆的扭转问题,编制了耦合法计算程序,计算了几种尺寸截面的抗扭刚度、截面各点的应力及裂纹的应力强度因子,并绘出了裂纹尖端的应力分布图.计算中,还探索了松弛因子对迭代收敛速度的影响.从实践上证实了自然边界元与有限元的耦合法所具有的优点.     

无界区域上Stokes问题的自然边界元与有限元耦合法

§1.引言 对于用有限元方法求解平面有界区域上的Stokes问题,国内外已有大量工作,例如可见[2]、[9]及其所引文献.但对无界区域上的这一问题,由于区域的无界性给有限元方法带来了困难,边界元方法及边界元与有限元的耦合法便显示其优越性.本文提出用自然边界元与有限元的耦合法求解无界区域上的Stokes问题.这一耦合法早在作者以前的工作中被应用于求解调和问题、重调和问题和平面弹性问题,但将它用于求解     

三维Helmholtz方程外问题的自然边界元与有限元耦合法

１．引言 设Γ０是空间闭曲面,Ω是Γ０外部的无界区域,考虑三维Ｈｅｌｍｈｏｌｔｚ 方程外Ｎｅｕｍａｎｎ问题其中 是波数,ｗ是频率,ｃ０为波在均匀介质中的传播速度,ｖ是区域Ω的边界Γ０的外法线方向,即指向由 Γ０包围的内部区域, 　　　　　　　为 Γ０上的已知函数．为了保证问题（１．１）和（１．２）的解的存在唯一性,必须附加上无穷远边界条件,即所谓的 Ｓｏｍｍｅｒｆｅｌｄ辐射条件其中ｉ是虚数单位,　　　　　． 许多数学物理问题,例如时间调和声波对不可穿透的障碍物的散射,海洋水下声波的传播,电磁波的绕射与辐…     

一种有限元-边界元耦合分域算法

提出了一种有限元-边界元耦合分域算法．该算法将所分析问题的区域分解成有限元和边界元子域,在满足两子域界面上位移和面力协调连续的条件下,通过迭代求解得到问题的解．在迭代求解过程中,引入动态松弛系数,使收敛得以加速．该方法在两子域界面上有限单元结点和边界单元结点的位置相互独立,无需协调一致,对诸如裂纹扩展过程的模拟具有独特的优势．用所提出的耦合算法分析算例,得到的结果与有限元法、边界元法和另一种耦合算法的数值计算结果一致,验证了这种算法的正确性和可行性．     

Signorini传输问题的有限元-边界元耦合

对一类带Signorini接触条件的非线性传输问题, 提出了新的有限元-边界元耦合框架.为求解所得到的耦合变分不等式, 设计了一种区域分解型迭代方法, 并对其做了完全的收敛性分析.     

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

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

解一类半线性外边值问题的自然边界元与有限元耦合法

In this paper, we combine a finite element approach with the natural boundary element method to stduy the weak solvability and Galerkin approximations of a class of semilinear exterior boundary value problems. Our analysis is mainly based on the variational formulation with constraints. We discuss the error estimate of the finite element solution and obtain the asymptotic rate of convergence O（h^n) Finally, we also give two numerical examples.     

热传导方程的有限元与边界积分方法

0 引言 设Ω是R~2中具有光滑边界гΩ的有界区域,Ω~cR~2\Ω.对任意实数T>0,记I[0,T].我们考虑如下初边值问题: u-Δu=f(x,t),(x,t)∈Ω~c×I, u(x,t)=0,(x,t)∈г×I, (0.1) u(x,0)=φ(x),x∈Ω~c. 现在我们引入一条嵌入Ω~c中的人工边界г_0(г_0可以是光滑的,也可以是不光滑的),г_0将Ω~c分为两部分:无界区域Ω_2;有界区域Ω_1,且假定suppfΩ_1,suppφΩ_1.用n=(n_1,n_2)表示г_0的单位外法向量.     

无界区域瞬时涡流问题有限元-边界元耦合的A-φ法的误差分析

针对三维无界区域带有凸多边形导体的瞬时涡流问题,本文提出了一种基于势场的有限元-边界元耦合的方法,从理论上讨论了其能量模误差估计.虽然电场被分解为电矢势A与磁标势φ的梯度之和后增加了方程与未知量的个数,但这种分解可以很好地处理不同介质间的间断.与传统的A-φ法不同,本文讨论了一种全离散的A-φ解耦形式,这样不仅可以避免传统格式所产生的鞍点问题的求解,又可以减少计算量.     

平面弹性力学问题边界元直接法中边界积分的解析处理

确立平面位势和弹性力学问题的边界元直接法中边界积分的解析计算框架系统,从而避免了传统的高斯近似求积分,数值算例表明它具有较高的精度和效率,特别是在边界量和边界附近区域内点物理量的计算可获得较高的精度。     

FE2 Modelling of Hybrid Sandwich Composites

This paper presents a procedure to simulate the behaviour of a sandwich composite material. The classical solution to model this kind of structure is the plate theory. Nevertheless, because of the incapacity to take non-linear material behaviour into account, this method cannot be used to model our problem. Indeed, the sandwich structure contains metallic layers, and their behaviour is elasto-plastic. For these reasons, a numerical homogenisation will be presented. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the h-adaptive coupling of FE and BE for viscoplastic and elasto-plastic interface problems

This paper presents a posteriori error estimates for the symmetric finite element and boundary element coupling for a nonlinear interface problem: A bounded body with a viscoplastic or plastic material behaviour is surrounded by an elastic body. The nonlinearity is treated by the finite element method while large parts of the linear elastic body are approximated using the boundary element method. Based on the a posteriori error estimates we derive an algorithm for the adaptive mesh refinement of the boundary elements and the finite elements. Its implementation is documented and numerical examples are included.     

A Hybrid Method for Nonlinear Least Squares Problems

A negative curvature method is applied to nonlinear least squares problems with indefinite Hessian approximation matrices. With the special structure of the method, a new switch is proposed to form a hybrid method. Numerical experiments show that this method is feasible and effective for zero-residual, small-residual and large-residual problems.     

An Improved Hybrid Method for Inverse Obstacle Scattering Problems

An improved hybrid method is introduced in this paper as a numerical method to reconstruct the scatterer by far-field pattern for just one incident direction with unknown physical properties of the scatterer. The improved hybrid method inherits the idea of the hybrid method by Kress and Serranho which is a combination of Newton and decomposition method, and it improves the hybrid method by introducing a general boundary condition. The numerical experiments show the feasibility of this method.     

A Hybrid Smoothing Method for Mixed Nonlinear Complementarity Problems

In this paper, we describe a new, integral-based smoothing method for solving the mixed nonlinear complementarity problem (MNCP). This approach is based on recasting MNCP as finding the zero of a nonsmooth system and then generating iterates via two types of smooth approximations to this system. Under weak regularity conditions, we establish that the sequence of iterates converges to a solution if the limit point of this sequence is regular. In addition, we show that the rate is Q-linear, Q-superlinear, or Q-quadratic depending on the level of inexactness in the subproblem calculations and we make use of the inexact Newton theory of Dembo, Eisenstat, and Steihaug. Lastly, we demonstrate the viability of the proposed method by presenting the results of numerical tests on a variety of complementarity problems.     

An Adaptive Hybrid Spectral Method for Stochastic Helmholtz Problems

The implementation of an adaptive hybrid spectral method for Helmholtz equations with random parameters is addressed in this work. New error indicators for generalized polynomial chaos for stochastic approximations and spectral element methods for physical approximations are developed, and systematic adaptive strategies are proposed associated with these error indicators. Numerical results show that these error indicators provide effective estimates for the approximation errors, and the overall adaptive procedure results in efficient approximation method for the stochastic Helmholtz equations.     

An Improved Hybrid Method for Inverse Obstacle Scattering Problems

An improved hybrid method is introduced in this paper as a numerical method to reconstruct the scatterer by far-field pattern for just one incident direction with unknown physical properties of the scatterer. The improved hybrid method inherits the idea of the hybrid method by Kress and Serranho which is a combination of Newton and decomposition method, and it improves the hybrid method by introducing a general boundary condition. The numerical experiments show the feasibility of this method.     

微极弹性动力学中非保守力场问题的变分方法

本文利用卷和卷的交换性质给出并证明了微极弹性动力学中非保守力场问题的几种拟变分原理。本文结果还可以推广到非局部弹性介质和非局部微极弹性介质力学中去。     

Parallel Strategies Applied to FE Simulations of Multiphasic Problems

While the theoretical background of various porous media models is well understood, it is still a demanding task to deal with these models numerically. In this contribution, a triphasic model is presented, which is capable of describing partially saturated soils. In quasi‐static conditions, this model results in the primary variables solid displacement, pore‐liquid pressure and pore‐gas pressure. For a stable numerical implementation, Taylor‐Hood elements are required, which need quadratic ansatz functions for the displacement and linear ansatz functions for the pressure terms. Looking at numerical simulations in 2‐d, challenging finite element calculations have already been realized in combination with adaptivity in time and space [1]. Nevertheless, new strategies have to be considered for a realization of applications of the model in 3‐d in order to handle the huge amount of unknowns arising from the discretization with Taylor‐Hood elements. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On a Hybrid Extragradient-Viscosity Method for Monotone Operators and Fixed Point Problems

The purpose of this article is to give a more general scheme for approximating a common element of the fixed-point set of a certain mapping and the set of solutions of a variational inequality problem. This scheme is inspired by the recent work of Maingé [A hybrid extragradient-viscosity method for monotone operators and fixed point problems, SIAM J. Control Optim. 47, 1499–1515 (2008)]. We also show that some assumption imposed in his result can be relaxed. Moreover, our scheme is a genuine generalization of Maingé's result because there is a class of mappings to which our scheme is applicable, but which is beyond the scope of his result.