Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 83
PROCEEDINGS OF THE EIGHTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY Edited by: B.H.V. Topping, G. Montero and R. Montenegro
Paper 164
Discrete Green's Functions for TimeHarmonic Wave Problems on Unbounded Domains with Periodic Variation of Material Properties B. Boroomand and F. Mossaiby
Civil Engineering Department, Isfahan University of Technology, Iran B. Boroomand, F. Mossaiby, "Discrete Green's Functions for TimeHarmonic Wave Problems on Unbounded Domains with Periodic Variation of Material Properties", in B.H.V. Topping, G. Montero, R. Montenegro, (Editors), "Proceedings of the Eighth International Conference on Computational Structures Technology", CivilComp Press, Stirlingshire, UK, Paper 164, 2006. doi:10.4203/ccp.83.164
Keywords: discrete Green's functions, unbounded domain, periodic material.
Summary
In this paper we present a numerical method suitable for solution of steady state
wave problems in which the material properties vary periodically throughout the
domain. Discrete Green's functions, in finite element sense, are selected as the
representatives of the problems. The Green's functions are evaluated on unbounded
domains and this involves satisfaction of radiation conditions in the solutions.
Having found the discrete Green's functions, one can use them in solution of more
general bounded or unbounded problems. To this end a simple superposition algorithm
is required, especially for general unbounded problems where the effects of the
radiation condition should be taken into account and the discrete Green's functions can
play a prominent role. For problems with homogenous materials, there are a number
of approaches for considering the effect of the radiation on the unbounded domain. The
boundary element method (BEM) is one of the approaches. Apart from the problem
of dealing with the inherent singularity of the functions in the BEM, the Green's functions
are not always available. An example is a model comprising a set of periodic
mechanical properties that we shall deal with in this paper.
Like many other engineering problems the finite element method (FEM) may also be used for general unbounded domains. Then the problem is considered to consist of a bounded part, discretized with elements, and an unbounded part. Nevertheless, for the unbounded part some special treatments are usually required. The use of infinite elements and artificial boundary conditions are recommended by many researchers during the past thirty years. For the such categories of treatments, the reader may consult review papers [1,2,4] or reference [3]. The aforementioned treatments for using the FEM are basically devised for models with homogenous mechanical properties. For models with periodic properties the available information, there is little to be found in the literature. Dealing with such problems needs behaviour prediction of the FEM when the entire unbounded domain is discretized with elements. For modelling unbounded domains with homogenous properties using the FEM, an early work by Thatcher is for Laplace's equations [5]. Somewhat similar approaches have been introduced by other scientists. The history can be found in reference [6]. A breakthrough has recently been introduced by the authors in [7] assuming that the entire unbounded domain is discretized into similar patterns of elements. The model is capable of solving problems with periodic material properties. The formulation given in this paper is the extension of the one given in [7]. Here, the principles of Floquet theory [8] for the solution of partial differential equations with periodic coefficients are used. First, the fundamental exponentiallike wave bases are obtained through the dispersion relations and then the radiation conditions are satisfied by selecting the wave bases. For selecting the wave bases, a quadrant of the main unbounded domain together with a set of appropriate boundary conditions is considered. For satisfaction of the boundary conditions a discrete transformation proposed by the authors in [7,9] is used. Application of the method is shown on a sample problem and the results are compared with those obtained from solution of a problem with homogenized material. Excellent agreement is seen between the amplitudes and the wavelengths of the solutions. However, the details of the deformation are just seen in the results of the numerical solution proposed. The method we present in this paper is suitable for more general complicated cases where homogenization cannot easily be used. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
