mohsen kamalian, Mohammad Saffar


Since the topography and surface irregularities of the terrain have a significant effect on the seismic response of the earth's surface, the numerical analysis using boundary element method in elastodynamics, which leads to a significant increase in the number of DOF and stiffness matrix dimension, as well as the sparse and unsymmetrical structure of stiffness matrix, indicates the inefficiency of the conventional method. This paper presents a fast boundary element algorithm in seismic analysis of two-dimensional elastic media for the first time. In the proposed method, instead of the usual node-to-node, a method of cell-to-cell relation method is implemented by hierarchy tree structure. In addition, the Plane Wave Time Domain algorithm and the Iterative method has been used to solve a system of equations which increases the speed of network convergence, especially in high degrees of freedom that has not been used in the study of the seismic waves’ response yet. Nowadays, this new algorithm adopts for investigating the response wave fields of electrical, magnet [1], acoustic and thermal conductivity which are discussed in details in various articles and books.


Plane Wave Time Domain (PWTD), Fast Boundary Element, Iterative Solver, Hierarchy tree, Degree of Freedom (DOF)

Full Text:



Lu M. Michielssen E. Shanker B. (2004)Fast time domain integral equation solvers for analysing two-dimentional scattering phenomena;part I,J. Electromagnetics, 24:425-449

Lu M. Michielssen E. Shanker B. (2004)Fast time domain integral equation solvers for analysing two-dimentional scattering phenomena;part II,J. Electromagnetics, 24:451-470

Nishimura n.,(2002) Fast multipole accelerated boundary integral equation methods. J. American society of mechanical Engineers, 55(4):299-324

Yarlin n. and Rokhlin V.(1999), an improved fast multipole algorithm for potential fields on the line,J. Numer. Anal.36,629-666

Greengard L. and Helsing J.,(1998), on the numerical evaluation of elastostatic fields in locally isotropic two-dimensional composits, J. Mech. Phys. Solids, 46, 1441-1462

Greengard L. and Kropinski MC.,(1996), Integral equation methods for stokes flow and isotropic elasticity in the plane, J. Comput. Phys.,125,403-414

Ergin A. Michielssen E. Shanker B. (1999), Fast transient analysis of acoustic wave scattering from rigid bodies using a two-level plane wave time domain algorithm, J. Acoust. Soc. Am., 106,2405-2416

Kamalian M. Gatmiri B. Sohrabi bidar A. (2003) On Time-Domain Two-Dimensional Site Response Analysis of Topographic Structures by BEM ,J JSEE, 5(2):35-45

Abramowitz M.,(1972) Handbook of mathematical functions: with formulas, graphs, and mathematical tables,. New York: Dover Publications

Ergin A. Michielssen E. Shanker B. (2000)Fast analysis of transient acoustic wave scattering from rigid bodies using the multilevel plane wave time domain algorithm. Acoustical Society of America, 107(3):1168-1178

H. Cheng, L. Greengard, and V. Rokhlin, (1999). “A fast adaptive multipole algorithm in three dimensions,” J. Comput. Phys. 155, 468–498

Ergin A. Michielssen E. Shanker B. (1998)Fast Evaluation of three-Dimensional Transient Wave Fields Using Diagonal Translation Operators. J. Journal of computational physics, 146,157-180

Takahashi T. Nishimura N. Kobayashi S.,(2001), Fast Boundary Integral Equation method for Elastodynamic Problems in 2D in Time Domain,661(67):1409-1416

Knab J. (1979), Interpolation of band-limited functions using the approximate prolate Series, J. IEEE Trans. Inform. Theory, 25,717-726


  • There are currently no refbacks.

Copyright (c) 2017 mohsen kamalian, Mohammad Saffar

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.