Modal analysis of plates using the dual reciprocity boundary element method

This paper is concerned with a dual reciprocity boundary element method (DRBEM) applied to the transient heat conduction problem of anisotropic materials. The integral equation formulation employs the fundamental solution of Laplace equation for isotropic materials, and hence from the anisotropic part of the governing differential equation a domain integral arises in the boundary integral equation. This domain integral is transformed into boundary integrals by using new radial basis functions. Furthermore, time derivative is approximated by the finite difference method, and the domain integral also appears from this approximation. Such a domain integral at each time step is also transformed into boundary integrals via the same dual reciprocity method. The details of the proposed DRBEM are presented, and a computer code is developed for two-dimensional problems. The usefulness of the present method is demonstrated through comparison of the results obtained by the computer code with the analytical ones.

Comparison between continuous and discontinuous boundary elements in the multidomain dual reciprocity method for the solution of the two-dimensional Navier–Stokes equations

The dual reciprocity method (DRM) and the radial integration method (RIM) are techniques used in BEM to transform domain integrals into boundary integrals. They are suitable for boundary element formulations where a complete fundamental solution is either unavailable or very complex, because in these cases one or more terms can remain as domain integrals in order to use a simpler fundamental solution. Thus, a large number of problems can be solved with the knowledge of a few number of fundamental solutions and additional terms as inertia or non-linear effects, can be treated as body forces and taken to the boundary. In both methods, the remaining terms are approximated through a finite series expansion involving proposed approximating functions and coefficients to be determined. This expansion is substituted in the generated domain integrals that are, subsequently, tr ansformed into boundary integrals. In this paper, a boundary element formulation for thin shallow shells with no domain discretization is presented. The domain integrals due to the curvature of the shells are tran sformed into boundary integrals using the radial integration method. Two approximation functions are used in this work. The first is the radial ba sis function that has been used extensively in the DRM given by:

The Kriging integration method applied to the boundary element analysis of Poisson problems

PurposeThe purpose of this paper is to describe further developments on a novel formulation of the boundary element method (BEM) for inelastic problems using the dual reciprocity method (DRM) but using object-oriented programming (OOP). As the BEM formulation generates a domain integral due to the inelastic stresses, the DRM is employed in a modified form using polyharmonic spline approximating functions with polynomial augmentation. These approximating functions produced accurate results in BEM applications for a range of problems tested, and have been shown to converge linearly as the order of the function increases.Design/methodology/approachA programming class named DRMOOP, written in C++ language and based on OOP, was developed in this research. With such programming, general matrix equations can be easily established and applied to different inelastic problems. A vector that accounts for the influence of the inelastic strains on the displacements and boundary forces is obtained.FindingsThe C++ DRMOOP class has been implemented and tested with the BEM formulation applied to classical elastoplastic problem and the results are reported at the end of the paper.Originality/valueAn object-oriented technology and the C++ DRMOOP class applied to elastoplastic problems.

A dual reciprocity boundary element method based on global shapes functions belonging to a complete set to transform domain integrals into boundary ones is presented. Governing equations of Poisson and Helmholtz type are tested. It is demonstrated that the complete set shape functions are more efficient and accurate in approximating domain integrals than those based on incomplete set. INTRODUCTION The application of Boundary Element Method to linear partial differential operators, or to nonlinear and time-dependent ones, often leads to domain integrals that are cumbersome to handle. The dual reciprocity method pioneered by Nardini and Brebbia* is a powerful numerical technique which converts domain integrals into boundary ones using a set of global shape functions to approximate the integrand. It has been successfully applied in many problems involving dynamic, diffusion and nonlinear type operators. See Partridge, et al. for a comprehensive bibliographical review. The success of the dual reciprocity hinges on the ability to express the global shape functions into the Laplacian of a set of known functions such that the divergence theorem can be applied to resolve the domain integral into a boundary one. Such a constraint has so far limited the selection of shape functions to a set of radial functions of the form 1 + ar -f br^ -f • • •• Transactions on Modelling and Simulation vol 1, © 1993 WIT Press, www.witpress.com, ISSN 1743-355X 344 Boundary Elements These radial functions do not form a complete set hence may lack convergence properties. Although practical evidences have indicated that those shape functions work quite satisfactorily, it is still desirable to use shape functions from a set that is complete in the mathematical sense. This paper present a dual reciprocity BEM based on global shape functions derived from a polynomial series. Perhaps due to its mathematical completeness, we are able to demonstrate that the numerical results are more efficient and accurate than the traditional radial set. DUAL RECIPROCITY For many non-standard, transient, or nonlinear partial differential equations it is possible to extract a standard, time-independent, linear partial differential operator C and lump the remainder to the right-hand-side: jC% = 6(z, ?/, Z, f , %, . . . , %,., %„ , %,, %,, %= . • •) (1) For the purpose of BEM implementation, these non-standard portions can be formally regarded as body forces. By the application of generalized Green's theorem = / Jr (2) where v is a conjugate variable, £* is the adjoint operator of £, B» and B* are operators resulting from integration by parts, which may be viewed as generalized normal derivatives, F is the solution boundary, and fi the solution domain, a boundary integral equation can be formed as: cu = / (u*q uq*) j bu* dSl (3) where quantities denoted by asterisk are free-space Green's functions satisfying £*u* = —6, with 6 a Dirac delta function, c is a constant dependent on the Cauchy principal value integration of the singularity at the base point, q and q* are generalized normal derivatives of u and u*. To circumvent the domain integration, the function b can be approximated by an expansion based on a set of global shape functions /, such that &*fX£ (4) j=l in which a, are coefficients determined by collocation, and n^ is the number of collocation nodes. The shape functions are chosen such that particular solutions of the set of inhomogeneous equations Ciij = /,• (5) Transactions on Modelling and Simulation vol 1, © 1993 WIT Press, www.witpress.com, ISSN 1743-355X Boundary Elements 345 are explicitly known. The domain integral in (3) can be transformed as [ / j^bu da = i ŷ /ju (6) j=l L ^1 J The last expression was a consequence of the application of (2). Substitution of the above into (3) yields / (u*q uq} dY = Y a, \cu, I (%*& titf) dT (7) •/r fr(. L 7r J cu We hence in principal have achieved the elimination of the domain integral. However, the detail of the implementation is somewhat dependent on the type of body force terms. We shall illustrate below the implementation by examples. EXAMPLES Operator of V^u = b(x, y) type For the governing equation of Poisson type, q in (7) is du/dn, and u* is — lnr/2%for 2-D or l/4%r for 3-D geometries. For simplicity, only two-dimensional problems will be studied herein. If the right-hand-side 6(x, y) is a mathematical expression whose particular solution is known as Up, the solution of (8) can be expressed as % = %& + %? (9) where u^ is the homogeneous solution and V2%& = 0 (10) It is apparent that Uh = u — Up satisfies the boundary integral equation c(u «„) = / [u(q q,} -(uu,)q'] (11) If we approximate the right-hand-side b(x,y) as (4), and find the particular solution term by term as V = /,(12) Transactions on Modelling and Simulation vol 1, © 1993 WIT Press, www.witpress.com, ISSN 1743-355X 346 Boundary Elements it is clear that j=l By the above interpretation, we have provided an alternative way for deriving (7), which is rewritten below with explicit reference to coordinates as c(x)u(x)f(u*(x',x)q(x')-u(x')q-(x',x)] dT(x') j r

Isogeometric dual reciprocity BEM for solving non-Fourier transient heat transfer problems in FGMs with uncertainty analysis

On the convergence of the dual reciprocity boundary element method

Boundary Element Analysis of Cracks in Shear Deformable Plates and Shells. Topics in Engineering, Vol 43

A subregion DRBEM formulation for the dynamic analysis of two-dimensional cracks

IG-DRBEM of three-dimensional transient heat conduction problems

Due to the application of cyclic loads, the existing crack length in a structure increases with time, resulting in greater stress concentration around the crack tip and therefore causing an increase in propagation velocity and a decrease in residual strength of the structure which, at a given time, becomes so low that the structure can no longer support the service loads. This paper presents a new approach to evaluate two-dimensional elastoplastic models in a crack propagation scenario, based on the boundary element method and its dual formulations. The methodology adopted consists of two stages: firstly, the simulation of the elastoplastic behavior, which considers the process of initial stresses and allows the treatment of several flow criteria, without taking into account the incompressibility of inelastic deformations. In this step, the domain integral due to non-homogeneous terms in the plastification region is transformed into a boundary integral by the Dual Reciprocity Method (DRM) using the polyharmonic splines functions due to its local behavior and, from the plastic calculus, the J-Integral is used to compute the SIF's. The second stage aims to simulate, in an incremental way, the crack propagation path using homemade software for crack modeling and analysis based on ElastoPlastic Fracture Mechanics (EPFM) and the Dual Boundary Element Method (DBEM). To validate the adopted methodology, as well as correctly simulate the mechanical behavior of cracks in the plastic regime of the material, three 2D models with straight and inclined cracks are used. Initially, the von Mises criterion for perfectly plastic materials is used and the numerical results compared with classic models from the literature, showing the effectiveness of the program in crack treatment and in the prediction of the propagation path, as well as the non-linear algorithm in the region plastification using the DRM.

The Dual Reciprocity Method is a well established technique for taking domain integrals to the boundary in Boundary Element analysis. It has been successfully applied to Diffusion Convection analysis on rectangular channels in a series of recent papers. Here the method is applied to effluent dispersion in the eastern part of the Solent, between Southampton and Portsmouth.

In this paper a 2D boundary element model is used to characterize the transient response of a piezoelectric based structural health monitoring system for cracked beam. The BE model is written for piezoelectric non-homogeneous problem employing generalized displacements. The dual reciprocity method is used to write the mass matrix in terms of boundary parameters only. The multidomain boundary element technique is implemented to model non-homogeneous and cracked configuration, unilateral interface conditions are also considered to prevent the physical inconsistence of the overlapping between interface nodes belonging to the crack surfaces. To assess the reliability and the effectiveness of the model numerical analyses are carried out on the modal and dynamic response of undamaged beam and results are compared with finite element calculations. Electrical response of piezoelectric sensors are then reported for different crack configurations in comparison with the undamaged case.

This paper describes the solution of a steady state natural convection problem in porous media by the dual reciprocity boundary element method (DRBEM). The boundary element method (BEM) for the coupled set of mass, momentum, and energy equations in two dimensions is structured by the fundamental solution of the Laplace equation. The dual reciprocity method is based on augmented scaled thin plate splines. Numerical examples include convergence studies with different mesh size, uniform and non-uniform mesh arrangement, and constant and linear boundary field discretizations for differentially heated rectangular cavity problems at filtration with Rayleigh numbers of Ra*=25, 50, and 100 and aspect ratios of A=1/2, 1, and 2. The solution is assessed by comparison with reference results of the fine mesh finite volume method (FVM). Copyright © 2000 John Wiley & Sons, Ltd.