A new coupling of mixed finite element and boundary element methods for an exterior Helmholtz problem in the plane

This paper deals with the scattering of time harmonic electromagnetic waves by an infinitely long cylinder containing a non-homogeneous conducting medium. More precisely, we study the transverse magnetic field that solves an interface problem holding between the cross section of the cylinder and the exterior two-dimensional free space. We apply a dual-mixed variational formulation in the obstacle coupled with a boundary integral equation method in the unbounded homogeneous space. A Fredholm alternative is utilized to prove that this continuous formulation is well posed. We define the corresponding discrete scheme by using the lowest order rotated Raviart-Thomas finite elements for the magnetic field and spectral elements for the boundary unknown. Then, we show that the resulting Galerkin scheme is uniquely solvable and convergent, and prove optimal error estimates. Finally, we illustrate our analysis with some results from computational experiments. This research was partially supported by CONICYT-Chile through the FONDAP Program in Applied Mathematics, by the Dirección de Investigación of the Universidad de Concepción through the Advanced Research Groups Program, and by the Ministerio de Educación y Ciencia of Spain, through the project No. MTM2004-05417.     

A mathematical model for beams on geosynthetic reinforced earth beds under strip loading

In this study, an attempt has been made to analyze a beam on geosynthetic reinforced earth beds subjected to strip loading. Geosynthetic layer has been assumed to have finite bending stiffness and therefore idealized as a beam. The foundation beam has been placed on compacted granular soil layer overlying the geosynthetic layer below which lies on the original weak/loose soil deposit. The upper dense and lower loose soil layers have been idealized as Winkler springs of different stiffnesses. Governing differential equations for the flexural response of the system have been derived and presented in non-dimensional form. These equations have been solved using appropriate boundary and continuity conditions. It was possible to obtain a closed form analytical solution for such a foundation system.     

A posteriori error analysis for a boundary element method with nonconforming domain decomposition

We present and analyze an a posteriori error estimator based on mesh refinement for the solution of the hypersingular boundary integral equation governing the Laplacian in three dimensions. The discretization under consideration is a nonconforming domain decomposition method based on the Nitsche technique. Assuming a saturation property, we establish quasireliability and efficiency of the error estimator in comparison with the error in a natural (nonconforming) norm. Numerical experiments with uniform and adaptively refined meshes confirm our theoretical results. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 947–963, 2014     

A new boundary element technique for solving plane problems of linear elasticity: 1. Theory

A new boundary elements technique for solving plane problems of linear elasticity theory is described. The method is based upon the Muskhelishvili complex variable representation of the displacement and stress fields involving two independent complex functions. These functions are represented by complex Cauchy integrals where the path of integration is taken around the external boundary of the solid. Two complex density functions appearing in the integrands of the Cauchy integrals are represented by spline functions and these are determined by the application of appropriate boundary conditions. The theory presented is suitable only for bounded simply-connected regions.     

Mathematical analysis of a viscoelastic problem with temperature‐dependent coefficients—Part II: Regularity

In this paper, we study local regularity properties of the stress solution of a quasistatic thermoviscoelastic problem whose behaviour law is of the Maxwell–Norton type with temperature‐dependent coefficients. Copyright © 2007 John Wiley & Sons, Ltd.     

Nonlinear heat equation for nonhomogeneous anisotropic materials: A dual‐reciprocity boundary element solution

A dual‐reciprocity boundary element method is presented for the numerical solution of initial‐boundary value problems governed by a nonlinear partial differential equation for heat conduction in nonhomogeneous anisotropic materials. To assess the validity and accuracy of the method, some specific problems are solved. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

A new element for analyzing large deformation of thin Naghdi shell model. Part II: Plastic

In this paper a new element is developed that is based on Cosserat theory. In the finite element implementation of Cosserat theory shear locking can occur, especially for very thin shells. In the present investigation the director vector is constrained to remain perpendicular to the mid surface during deformation. It will be shown that this constraint yields accurate results in very large deformation of thin shells also the rate of convergency is very good. For plastic formulation, the model introduced by Simo is used and it has been reduced for constrained director vector and the consistent elasto-plastic tangent moduli is extracted for finite element solution. This model includes both kinematic and isotropic hardening. For numerical investigations an isoparametric nine node element is employed then by linearization of the principle of virtual work, material and geometric stiffness matrices are extracted. The validity and the accuracy of the proposed element is illustrated by the numerical examples and the results are compared with those available in the literature.     

A FE‐BE coupling for a fluid‐structure interaction problem: Hierarchical a posteriori error estimates

In this article, we establish a hierarchical a posteriori error estimate for a coupling of finite elements and boundary elements for a fluid‐structure interaction problem posed in two and three dimensions. These methods combine boundary elements for the exterior fluid and finite elements for the elastic structure. We consider two weak formulations, a nonsymmetric one and a symmetric one, which are both uniquely solvable. We present the reliability and efficiency of the error estimates. For the two dimensional case, we compute local error indicators which allow us to develop an adaptive mesh refinement strategy on triangles. For the three dimensional case, we use hexahedrons as elements. Numerical experiments underline our theoretical results. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012