Computational & Technology Resources
an online resource for computational,
engineering & technology publications
PROCEEDINGS OF THE EIGHTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY
Edited by: B.H.V. Topping, G. Montero and R. Montenegro
Discrete Green's Functions for Time-Harmonic 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 Time-Harmonic 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", Civil-Comp Press, Stirlingshire, UK, Paper 164, 2006. doi:10.4203/ccp.83.164
Keywords: discrete Green's functions, unbounded domain, periodic material.
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 .
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 . Somewhat similar approaches have been introduced by other scientists. The history can be found in reference . A breakthrough has recently been introduced by the authors in  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 . Here, the principles of Floquet theory  for the solution of partial differential equations with periodic coefficients are used. First, the fundamental exponential-like 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.
purchase the full-text of this paper (price £20)