Search

In this paper, we study the relation between the error estimate of the bilinear interpolation on a general quadrilateral and the geometric characters of the quadrilateral. Some explicit bounds of the interpolation error are obtained based on some sharp estimates of the integral over 1/∣J∣p-1 for 1 ≤ p≤∞ on the reference element, where J is the Jacobian of the nonaffine mapping. This allows us to introduce weak geometric conditions (depending on p) leading to interpolation error estimates in the W1,p norm, for any p ϵ [1,∞), which can be regarded as a generalization of the regular decomposition property (RDP) condition introduced in [G. Acosta and R. G. Durán, SIAM J. Numer. Anal., 38 (2000), pp. 1073-1088] for p = 2 and new RDP conditions (NRDP) for p ≠ 2. We avoid the use of the reference family elements, which allows us to extend the results to a larger class of elements and to introduce the NRDP condition in a more unified way. As far as we know, the mesh condition presented in this paper is weaker than any other mesh conditions proposed in the literature for any p with 1 ≤ p≤∞. [ABSTRACT FROM AUTHOR]

Abstract: The hybrid power flow analysis (PFA) is an analytic method proposed for the effective prediction of vibrational and acoustic responses of low-damping system in the medium-to-high frequency ranges by using the PFA algorithm and statistical energy analysis (SEA) coupling concepts. This paper presents the hybrid boundary condition on 3-D case for hybrid PFA in addition to 1-D and 2-D cases which are derived in the other companion paper, and formulates the hybrid power flow finite-element method (PFFEM) including coupling loss factor (CLF) of SEA to extend the application area of hybrid PFA to built-up structures. To verify the derived boundary condition and hybrid PFFEM, numerical analyses were successfully performed for various analytic models and reverberance factors. [Copyright &y& Elsevier]

Abstract: In this paper, wave interaction with a bottom-mounted uniform porous cylinder of arbitrary shape is studied using SBFEM (Scaled Boundary Finite-element Method) and domain decomposition technique. The non-homogeneous equation due to the complex configuration of the structure is processed by introducing a porous effect parameter G. The resulting final scaled boundary finite-element equation is homogeneous and can be solved in conventional manners. Only the boundary is spatially discretized with curved surface finite elements, leading a reduction of the spatial dimension by one. The solution in the radial direction is analytical; while in the circumferential direction, the solution is numerical. High accuracy is achieved with a small number of surface finite elements. [Copyright &y& Elsevier]

A recent study [D. N. Arnold, D. Boffi, and R. S. Falk, SIAM J. Numer. Anal., 42 (2005), pp. 2429-2451] reveals that convergence of finite element methods using H(div , O)-compatible finite element spaces deteriorates on nonaffine quadrilateral grids. This phenomena is particularly troublesome for the lowest-order Raviart-Thomas elements, because it implies loss of convergence in some norms for finite element solutions of mixed and least-squares methods. In this paper we propose reformulation of finite element methods, based on the natural mimetic divergence operator [M. Shashkov, Conservative Finite Difference Methods on General Grids, CRC Press, Boca Raton, FL, 1996], which restores the order of convergence. Reformulations of mixed Galerkin and leastsquares methods for the Darcy equation illustrate our approach. We prove that reformulated methods converge optimally with respect to a norm involving the mimetic divergence operator. Furthermore, we prove that standard and reformulated versions of the mixed Galerkin method lead to identical linear systems, but the two versions of the least-squares method are veritably different. The surprising conclusion is that the degradation of convergence in the mixed method on nonaffine quadrilateral grids is superficial, and that the lowest-order Raviart-Thomas elements are safe to use in this method. However, the breakdown in the least-squares method is real, and there one should use our proposed reformulation. [ABSTRACT FROM AUTHOR]

Abstract: In this paper, we study the superconvergence of the frictionless Signorini problem. When approximated by bilinear finite elements, by virtue of the information on the contact zone, we can derive a superconvergence rate of under a proper regularity assumption. Finally, a numerical test is given to verify our result. [Copyright &y& Elsevier]

Abstract  The purpose of the paper is to apply monotone multigrid methods to static and dynamic biomechanical contact problems. In space, a finite element method involving a mortar discretization of the contact conditions is used. In time, a new contact-stabilized Newmark scheme is presented. Numerical experiments for a two body Hertzian contact problem and a biomechanical application are reported. [ABSTRACT FROM AUTHOR]

A cohesive interface element is presented for the finite element analysis of crack growth in thin specimens. In this work, the traditional cohesive interface model is extended to handle cracks in the context of three-dimensional shell elements. In addition to the traction-displacement law, a bending moment-rotation relation is included to transmit the moment and describe the initiation and propagation of cracks growing through the thickness of the shell elements. Since crack initiation and evolution are a natural outcome of the cohesive zone model without the need of any ad hoc fracture criterion, this model results in automatic prediction of fracture. In particular, this paper will focus on cases involving mode I/III fracture and bending, typical of complex cases existing in industrial applications in which thin-walled structures are subjected to extreme loading conditions (e.g., crashworthiness analysis). Finally, we will discuss how the three-dimensional effects near the crack front may affect the determination of the cohesive parameters to be used with this model. [ABSTRACT FROM AUTHOR]

A finite element discretization is proposed and analyzed for a linear parabolic problems with discontinuous coefficients. Due to low global regularity of the solution, it seems difficult to achieve optimal order of convergence with classical finite element methods [Numer. Math., 79 (1998), pp. 175--202]. In this paper, we have used a finite element discretization, where interface triangles are assumed to be curved triangles instead of straight triangles as in classical finite element methods. Optimal order error estimates in L2 and H1 norms are shown to hold even if the regularity of the solution is low on the whole domain. While the continuous time Galerkin method is discussed for the spatially discrete scheme, the discontinuous Galerkin method is analyzed for the fully discrete scheme. The interfaces and boundaries of the domains are assumed to be smooth for our purpose. [ABSTRACT FROM AUTHOR]

The combination of geometrically nonlinear analysis with structural optimization against instability opens many possibilities for new formulations of the optimization problem. One of them is optimization against instability in the large. When a critical state does not exist and instability occurs at finite displacements either the lower critical loading is maximized or the generalized displacement for that load is minimized. This paper undertakes these new optimization problems providing formulations and numerical solutions for selected elements. The simple finite dimensional model is analysed first showing characteristic features of the problem, and then a real elastic element is optimized. [ABSTRACT FROM AUTHOR]

We consider strategies to construct optimal order two- and multilevel hierarchical preconditioners for linear systems as arising from the finite element discretization of self-adjoint second order elliptic problems using non-conforming Crouzeix–Raviart linear elements. In this paper we utilize the hierarchical decompositions, derived in a previous work by the same authors (Numerical Linear Algebra with Applications 2004; 11:309–326) and provide a further analysis of these decompositions in order to assure robustness with respect to anisotropy. Finally, we show how to construct both multiplicative and additive versions of the algebraic multilevel iteration preconditioners and show robustness and optimal order convergence estimates. Copyright © 2005 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

A partial word of length n over a finite alphabet A is a partial map from {0,…, n − 1} into A. Elements of {0,…, n − 1} without image are called holes (a word is just a partial word without holes). A fundamental periodicity result on words due to Fine and Wilf [1] intuitively determines how far two periodic events have to match in order to guarantee a common period. This result was extended to partial words with one hole by Berstel and Boasson [2] and to partial words with two or three holes by Blanchet-Sadri and Hegstrom [3]. In this paper, we give an extension to partial words with an arbitrary number of holes. [Copyright &y& Elsevier]

In this paper, a mortar finite element method for parabolic problem is presented. Multigrid method is used for solving the resulting discrete system. It is shown that the multigrid method is optional, i.e., the convergence rate is independent of the mesh size L and the time step parameter Τ. [ABSTRACT FROM AUTHOR]

A finite-element implementation for the analysis of propagation constants and fields in a guiding structure is presented in this paper. The solution explicitly takes into account the field divergence on the edges of the internal wedges by exploiting substitutive divergent bases. © 2005 Wiley Periodicals, Inc. Microwave Opt Technol Lett 44: 327–328, 2005; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/mop.20624 [ABSTRACT FROM AUTHOR]