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.     

三维热传导型半导体问题的特征混合元方法和分析

本文研究三维热传导型半导体态问题的特征混合元方法及其理论分析,其数学模型是一类非线性偏微分方程的初边值问题,对电子位势方程提出混合元逼近,对电子,空穴浓度方程笔挺表限元逼近;对热传导方程采用对时间向后差分的Ｇａｌｅｒｋｉｎ逼近,应用微分方程先验估计理论和技巧得到了最优阶Ｌ＾２误差估计。     

Rigid‐plastic hybrid element analyses of the plane strain upsetting

A rigid‐plastic hybrid element method (HEM) for simulation of metal forming is developed. This method is a mixed approach of the rigid‐plastic domain‐BEM and the rigid‐plastic FEM based on the theory of compressible plasticity. Because the compatibilities of not only velocity but also velocity's derivative between the adjoining boundary elements and finite elements can be met, the velocities and the derivatives of the velocity can be calculated with the same precision for the rigid‐plastic HEM. Then, it is considered that the rigid‐plastic HEM is a more precise method in formulation than the conventional rigid‐plastic FEMs for which the compatibilities of velocity's derivative cannot be met. The plane strain upsetting processes with two friction factors are analyzed by the rigid‐plastic HEM in this article. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 726–737, 2002; Published online in Wiley InterScience (www.interscience.wiley.com); DOI 10.1002/num.10031.     

三维热传导型半导体问题的交替方向特征有限元方法及理论分析

１．引言 三维热传导型半导体器件瞬态问题的数学模型由四个非线性偏微分方程描述 ［１,２］．工程研究中一般考虑绝流边条件,由于绝流条件可以看作一反射条件来处理、为了数值分析方便,我们在此考虑三维周期问题： 其中,　＝［０,１］３,未知函数是电子位势　;电子,空穴浓度ｅ,ｐ;温度函数Ｔ．方程（１,１）－（１.４）中出现的系数均有正的上下界,且是　周期的． ａ＝Ｑ／ε,Ｑ,ε分别表示电子负荷和介电系数,均为正常数．Ｎ（ｘ）是给定的函数．Ｄｓ（ｘ）为扩散系数,μｓ（ｘ）为迁移率,ｓ＝ｅ,Ｐ．Ｒ（ｅ,ｐ,Ｔ）…     

Numerical solution for the model of vibrating elastic membrane with moving boundary

In this work, we are interested in obtaining an approximated numerical solution for the model of vibrating elastic membranes with moving boundary. The model is an extension of Kirchhoff’s model, which takes into account the change of size during the vibration. We apply the finite element method with a finite difference method in time to obtain an approximated numerical solution. Some numerical experiments are presented to show the effect of moving boundary effects in vibrating elastic membranes.     

Three-dimensional simulation of flat rolling through a combined finite element–boundary element approach

This paper presents an innovative approach for analysing three-dimensional flat rolling. The proposed approach is based on a solution resulting from the combination of the finite element method with the boundary element method. The finite element method is used to perform the rigid–plastic numerical modelling of the workpiece allowing the estimation of the roll separating force, rolling torque and contact pressure along the surface of the rolls. The boundary element method is applied for computing the elastic deformation of the rolls. The combination of the two numerical methods is made using the finite element solution of the contact pressure along the surface of the rolls to define the boundary conditions to be applied on the elastic analysis of the rolls. The validity of the proposed approach is discussed by comparing the theoretical predictions with experimental data found in the literature.     

三维热传导型半导体问题的特征有限元方法和分析

本文研究三维热传导型半导体瞬态问题的特征有限元方法及其理论分析,其数学模型是一类非线性偏微分方程的初边值问题,对电子位势方程提出Ｇａｌｅｒｋｉｎ逼近;对电子,空穴浓度方程采用特征有限元逼近;对热传导方程采用对时间向后差分的Ｇａｌｅｒｋｉｎ逼近．应用微分方程先验估计理论和技巧得到了最优阶Ｌ＾２误差估计。     

A finite element method for granular flow through a frictional boundary

A finite element method for the flow of dry granular solids through a domain involving a frictional contact boundary is formulated. The granular material is assumed as a compressible viscous-elastic–plastic continuum. Based on the principles of continuum mechanics, a complete set of equations is developed. The resulting boundary value problem is solved by the finite element method in space and by the finite difference method in time. The derivation of the finite element equations and the mathematical framework of the numerical technique are presented, together with two illustrative examples to demonstrate the validity of the technique.     

Well-Posedness and Finite Element Approximations for Elliptic SPDEs with Gaussian Noises

The paper studies the well-posedness and optimal error estimates of spectral finite element approximations for the boundary value problems of semi-linear elliptic SPDEs driven by white or colored Gaussian noises. The noise term is approximated through the spectral projection of the covariance operator, which is not required to be commutative with the Laplacian operator.Through the convergence analysis of SPDEs with the noise terms replaced by the projected noises, the well-posedness of the SPDE is established under certain covariance operator-dependent conditions. These SPDEs with projected noises are then numerically approximated with the finite element method. A general error estimate framework is established for the finite element approximations. Based on this framework, optimal error estimates of finite element approximations for elliptic SPDEs driven by power-law noises are obtained. It is shown that with the proposed approach, convergence order of white noise driven SPDEs is improved by half for one-dimensional problems, and by an infinitesimal factor for higher-dimensional problems.     

FINITE ELEMENT NONLINEAR GALERKIN COUPLING METHOD FOR THE EXTERIOR STEADY NAVIER-STOKES PROBLEM

1.IntroductionNonlinearGalerkinmethodsaremultilevelschemesforthedissipativeevolutionpartialdifferentialequations.Theycorrespondtothesplittingsoftheunknownu:u=y z)wherethecomponentsareofdifferentorderofmagnitudewithrespecttoaparameterrelatedtothespati...     

ALTERNATING-DIRECTION MULTISTEP PRECONDITIONED ITERATIVE METHODS FOR SEMICONDUCTOR PROBLEM WITH HEAT-CONDUCTION

The model of transient behavior of semiconductor with heat-conduction is an initial and boundary problem. Alternating-direction multistep preconditioned iterative methods and theory analyses are given in this paper. Electric potential equation is approximated by mixed finite element method, concentration and heat-conduction equations are approximated by Galerkin alternating-direction multistep methods. Error estimates of optimal order in L2 are demonstrated.     

Hybrid boundary element formulation for acoustic wave propagation

The so‐called hybrid stress boundary element method (HSBEM) is introduced in a frequency domain formulation for the computation of acoustic radiation and scattering in closed and in finite domains. Different from other boundary element formulations, the HSBEM is based on an extended Hellinger‐Reissner variational principle and leads to a Hermitian, frequency‐dependent stiffness equation. Due to this, the method is very well suited for treating fluid structure interaction problems since the effort for the coupling the structure, discretized by a finite elements, and the fluid, discretized by the HSBEM is strongly reduced. To arrive at a boundary integral formulation, the field variables are separated into boundary variables, which are approximated by piecewise polynomial functions, and domain variables, which are approximated by a superposition of singular fundamental solutions weighed by source strength. This approximation cancels the domain integral over the equation of motion in the hybrid principle and leads to a boundary integral formulation, incorporating singular integrals. Comparing to previous results published by the authors, new considerations concerning the interpretation of singular contributions in the stiffness matrix for exterior domain problems are communicated here.     

A characteristics-mixed finite element method for Burgers’ equation

In this paper, we propose a new mixed finite element method, called the characteristics-mixed method, for approximating the solution to Burgers’ equation. This method is based upon a space-time variational form of Burgers’ equation. The hyperbolic part of the equation is approximated along the characteristics in time and the diffusion part is approximated by a mixed finite element method of lowest order. The scheme is locally conservative since fluid is transported along the approximate characteristics on the discrete level and the test function can be piecewise constant. Our analysis show the new method approximate the scalar unknown and the vector flux optimally and simultaneously. We also show this scheme has much smaller time-truncation errors than those of standard methods. Numerical example is presented to show that the new scheme is easily implemented, shocks and boundary layers are handled with almost no oscillations. One of the contributions of the paper is to show how the optimal error estimates inL ²(Ω) are obtained which are much more difficult than in the standard finite element methods. These results seem to be new in the literature of finite element methods.     

三维热传导型半导体问题的交替方向特征有限元方法

提出交替方向特征有限元方法,对电场位势方程采用混合元格式,对电子,空穴浓度方程采用交替方向特征有限元格式,对温度方程提出交替方向格式.应用向量积计算及先验估计理论和技巧,得到最佳的L2误差估计.     

n-sided polygonal hybrid finite elements with unified fundamental solution kernels for topology optimization

In topology optimization, the optimized design can be obtained based on spatial discretization of design domain using natural polygonal finite elements to reduce the influence of mesh geometry on topology optimization solutions. However, the natural polygonal finite elements require separate interpolants for each type of elements and involve troublesome domain integrals. In this study, an alternative n-sided polygonal hybrid finite element possessing multiple-node connection is formulated in a unified form to compress the checkerboard patterns caused by numerical instability in topology optimization. Different from the natural polygonal finite elements, the present polygonal hybrid finite elements involve two sets of independent displacement fields. The intra-element displacement field defined inside the element is approximated by the linear combination of the fundamental solution of the problem to achieve the purpose of the local satisfaction of the governing equations of the problem, but not the specific boundary conditions and the inter-element continuity conditions. To overcome such drawback, the inter-element displacement field defined over the entire element boundary is independently approximated by means of the conventional shape function interpolation. As a result, only line integrals along the element boundary are involved in the computation, whose dimension is reduced by one compared to the domain integrals in the natural polygonal finite elements, and more importantly, allowing us to flexibly construct any polygons from Voronoi tessellations in discretizing complex design domains using same fundamental solution kernels. Numerical results obtained indicate that the present n-sided polygonal hybrid finite elements can produce more accurate displacement solutions and smaller mean compliance, compared to the standard finite elements and the natural polygonal finite elements.     

三维热传导型半导体问题的交替方向有限元方法

1　引　　言三维热传导型半导体器件瞬态问题的数学模型由四个非线性偏微分方程描述[1 ,2 ] ,记 Ω为 Ω=[0 ,1 ] 3的边界 ,三维问题-Δψ =α( p -e+ N( x) ) ,　　 ( x,t)∈Ω× [0 ,T] ,( 1 .1 ) e t= . ( De( x) e-μe( x) e ψ) -R( e,p,T) ,　 ( x,t)∈Ω× ( 0 ,T] ,( 1 .2 ) p t= . ( Dp( x) p +μp( x) p ψ) -R( e,p,T) ,　 ( x,t)∈Ω× ( 0 ,T] ,( 1 .3 )ρ( x) T t-ΔT =[( Dp( x) p +μp( x) p ψ) -( De( x) e-μe( x) e ψ) ] . ψ,　　　　　　 ( x,t)∈Ω× ( 0 ,T] . ( 1 .4 )ψ( x,t) =e( x,t) =p( …     

On the symmetric boundary element method and the symmetric coupling of boundary elements and finite elements

The finite element method and the boundary element method areamong the most frequently applied tools in the numerical treatmentof partial differential equations. However, their propertiesappear to be complementary: while the boundary element methodis appropriate for the most important linear partial differentialequations with constant coefficients in bounded or unboundeddomains, the finite element method seems to be more appropriatefor inhomogeneous or even nonlinear problems. but is somehowrestricted to bounded domains. The symmetric coupling of thetwo methods inherits the advantages of both methods. This paper treats the symmetric coupling of finite elementsand boundary elements for a model transmission problem in twoand three dimensions where we have two domains: a bounded domainwith nonlinear, even plastic material behaviour, is surroundedby an unbounded, exterior, domain with isotropic homogeneouslinear elastic material. Practically. the coupling is performedsuch that the boundary element method contributes a macro-element,like a large finite element, within a standard finite elementanalysis program. Emphasis is on two-dimensional problems wherethe approach using the Poincaré-Steklov operator seemsto be impossible at first glance. E-mail: cc{at}numerik.uni-kiel.de     

Recovery Type a Posteriori Error Estimates of Fully Discrete Finite Element Methods for General Convex Parabolic Optimal Control Problems

This paper is concerned with recovery type a posteriori error estimates of fully discrete finite element approximation for general convex parabolic optimal control problems with pointwise control constraints. The time discretization is based on the backward Euler method. The state and the adjoint state are approximated by piecewise linear functions and the control is approximated by piecewise constant functions. We derive the superconvergence properties of finite element solutions. By using the superconvergence results, we obtain recovery type a posteriori error estimates. Some numerical examples are presented to verify the theoretical results.     

Particle Finite Element Simulations of Cutting Processes in Manufacturing

The cutting of metals is an important process in manufacturing and challenges established methods in the field of computational mechanics. The particle finite element method (PFEM) combines the benefits of particle based methods and the standard finite element method (FEM) to account for large deformations and separation of material. In cutting simulations the workpiece is realised as a set of particles, whose boundary is detected by the α-shape method. After the boundary detection, the particles are meshed with finite elements. Since metals show a plastic behavior under large deformations, a suitable material model needs to be considered. Numerical examples show the effect of the choice of the parameter α on the cutting force. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)