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 171
A New Boundary Element Method Formulation in ThreeDimensional Exterior Elastodynamics C.G. Provatidis and N.K. Zafiropoulos
Department of Mechanical Engineering, National Technical University of Athens, Greece C.G. Provatidis, N.K. Zafiropoulos, "A New Boundary Element Method Formulation in ThreeDimensional Exterior Elastodynamics", 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 171, 2006. doi:10.4203/ccp.83.171
Keywords: boundary element method, elastodynamics, fictitious eigenvalues.
Summary
It is well known that when the boundary element method (BEM) is applied to
exterior acoustic problems, then unrealistic fictitious eigenfrequencies appear; thus
the overall numerical solution becomes inaccurate in certain large regions of
frequencies [1]. In brief, the fictitious eigenvalues of an exterior problem under
Neumann boundary conditions are those of its dual internal problem under Dirichlet
boundary conditions and viceversa. An effective remedy to this shortcoming is the
dominating technique of Burton and Miller [2] while a detailed literature survey of
most other known techniques has been previously reported by the authors [3].
Similar fictitious eigenvalues appear also in the BEM solution of exterior timeharmonic elastodynamics [4], which is the topic of this paper. However, apart from older theoretical comments, a detailed report on computational results related to the numerical treatment of this particular shortcoming in exterior elastodynamics is still lacking in the literature. Even for the simplest case of the combined Helmholtz integral equation formulation (CHIEF) method [1], which is well known in acoustics, numerical results in elastodynamics were not available in the literature. The disadvantage of CHIEF is that it uses additional equations that are written for internal nodes that should be carefully selected. As a result of to the larger number of equations than the boundary unknowns, a least squares scheme should be applied. Within this context, this paper proposes a new technique to circumvent the fictitious eigenfrequencies in the BEM solution of the exterior timeharmonic elastodynamics and thoroughly investigates its performance. The novel feature of this paper is the derivation of a new elastodynamical operator for the derivation of the initial integral equation, which was obtained through the following two steps. First, the NaviérCauchy equation was written in vector form and was compared with the equivalent acoustic equation. From this comparison a correspondence of partial differential operators between elastodynamics and acoustics was obtained. Furthermore, by expressing in the acoustical analogue the sound velocity in terms of the bulk modulus and mass density , the original internal equation and its derivative were written in a modified way so that material properties appear as well. As a result, the comparison yielded that the elasticity matrix corresponds to the bulk modulus . Second, in order to avoid the hypersingular integrals involved in the differentiated integral equations, in this work it was decided to apply a modified scheme previously proposed by Cunefare et al. [5], according to which the differentiation is performed perpendicularly to an internal surface, for example an offset of the real domain boundary. By intuitively extending the meaning of the normal derivative over the internal boundary to a more complex operator that includes both geometrical and material quantities, the aforementioned two steps were successfully combined, thus a new elastodynamical operator suitable for the differentiation of the initial integral equation was identified. The proposed technique was applied to a cubic enclosure of unit size excited by a singular source. The numerical study includes many alternative choices (real and imaginary) of the coupling coefficient between the original and the differentiated integral equation and finally proposes a specific formula for which the most accurate numerical solution is achieved. It was also found that the proposed method is more accurate than the usual CHIEF method [1]. This is an important conclusion that suggests that the new method might replace the CHIEF method; an additional reason for the replacement might be the instabilities coming from the least squares procedure used in the CHIEF method. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
