V. Mallardo¹, M.H. Aliabadi²

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.2, pp. 77-96, 2012, DOI:10.3970/cmes.2012.087.077

Abstract The present paper intends to couple the Fast Multipole Method (FMM) with the Boundary Element Method (BEM) in the 2D scalar wave propagation. The procedure is aimed at speeding the computation of the integrals involved in the governing Boundary Integral Equations (BIEs) on the basis of the distance between source point and integration element. There are three main contributions. First, the approach is of adaptive type in order to reduce the number of floating-point operations. Second, most integrals are evaluated analytically: the diagonal and off-diagonal terms of the H and G matrices by consolidated techniques, More >

M.R. Hematiyan^1,2, A. Khosravifard¹, Y.C. Shiah³, C.L. Tan⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.1, pp. 55-76, 2012, DOI:10.3970/cmes.2012.087.055

Abstract An inverse technique, based on the boundary element method (BEM) and elastostatic experiments for identification of elastic constants of orthotropic and general anisotropic 2D bodies is presented. Displacement measurements at several points on the boundary of the body, obtained by a few known load cases are used in the inverse analysis to find the unknown elastic constants of the body. Using data from more than one elastostatic experiment results in a more accurate and stable solution for the identification problem. In the inverse analysis, sensitivities of displacements of only boundary points with respect to the More >

Yoshihiro Ochiai¹, Vladimir Sladek², Jan Sladek²

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.1, pp. 41-54, 2012, DOI:10.3970/cmes.2012.087.041

Abstract The conventional boundary element method (BEM) requires a domain integral in unsteady thermal stress analysis with heat generation or an initial temperature distribution. In this paper it is shown that the three-dimensional unsteady thermal stress problem can be solved effectively using the triple-reciprocity boundary element method without internal cells. In this method, the distributions of heat generation and initial temperature are interpolated using integral equations and time-dependent fundamental solutions are used. A new computer program was developed and applied to solving several problems. More >

Y.C. Shiah¹, C. L. Tan², C.Y. Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.1, pp. 1-22, 2012, DOI:10.3970/cmes.2012.087.001

Abstract The fundamental solution, or Green's function, for 3D anisotropic elastostatics as derived by Ting and Lee (1997) [Q.J. Mech. Appl. Math.; 50: 407-426] is one that is fully explicit and algebraic in form. It has, however, only been utilized in boundary element method (BEM) formulations quite recently even though it is relatively straightforward and direct to implement. This Green's function and its derivatives are necessary items in this numerical analysis technique. By virtue of the periodic nature of the angles when it is expressed in the spherical coordinate system, the present authors have very recently… More >

Miaoxia Xie¹, Yueming Li¹, Hualing Chen^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.1, pp. 65-78, 2012, DOI:10.3970/cmes.2012.085.065

Abstract Energy finite element method(EFEM) is a promising method to solve high-frequency vibro-acoustic problem. Energy boundary element method (EBEM) is an effective way to compute high-frequency sound radiation in the unbounded medium. Vibro-acoustic coupling of cavity structure in anechoic chamber includes both the interior acoustic field and unbounded exterior acoustic field. In order to predict this kind of high-frequency vibro-acoustic coupling problem in anechoic chamber, an approach combined EFEM and EBEM is developed in this paper. As a numerical example, the approach is applied to solve the high-frequency vibro-acoustic coupling response of a cubic cavity structure More >

Shuncai Li^1,2,3, Shichuang Zhuo⁴, Qiang Zhang⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.84, No.3, pp. 283-296, 2012, DOI:10.3970/cmes.2012.084.283

Abstract Based on the complex functions theory in elastic mechanics, the bending deflection formula expressed by the complex Fourier series is derived for the infinite plate with a circular opening at first, then the boundary conditions of the circular opening are expanded in Fourier Series, and the unknown coefficients of the Fourier series are determined by comparing coefficients method. By means of the convolution of the complex Fourier series and some basic formulas in the generalized functions theory, the natural boundary integral formula or the analytical deflection formulas expressed by the boundary displacement or loads are… More >

M. Koyuncu¹, F. Y. Ikikat¹, G. C. Icoz², B. Baranoglu³, A. Yazici²

CMES-Computer Modeling in Engineering & Sciences, Vol.84, No.1, pp. 13-26, 2012, DOI:10.3970/cmes.2012.084.013

Abstract A parallel dual reciprocity boundary element method solution to thermoelasticity and thermoviscoelasticity problems is proposed. The DR-BEM formulation is given in Fourier Transform Space where the Time Space solutions are obtained through inverse Fourier Transform. The parallellization of the code is achieved through solving each frequency at a distinct computational node. The implemented parallel code is tested on 64-core IBM blade servers and it is seen that a linear speed-up is achieved. More >

C. J. ZHENG, H. B. CHEN, T. MATSUMOTO, T. TAKAHASHI

The International Conference on Computational & Experimental Engineering and Sciences, Vol.19, No.4, pp. 121-122, 2011, DOI:10.3970/icces.2011.019.121

Abstract A fast multipole boundary element method is presented for three dimensional acoustic shape sensitivity analysis in this study. The Burton-Miller formula which is a linear combination of the conventional boundary integral equation and the normal derivative boundary integral equation is adopted to conquer the fictitious eigenfrequency problem associated with the conventional boundary integral equation method in solving exterior acoustic problems. The continuous adjoint variable method is implemented in the sensitivity analysis and the concept of material derivative is used in the derivation. Constant elements are employed to discretize the boundary so that the hypersingular boundary More >

Zai You YAN

The International Conference on Computational & Experimental Engineering and Sciences, Vol.17, No.3, pp. 77-78, 2011, DOI:10.3970/icces.2011.017.077

Abstract In the boundary element method, treatment of all the possible singular integrals is very important for the correctness and accuracy of the solutions. Generally, the directional derivative of a weakly singular integral is computed by an integral in the sense of Cauchy principle value if the directional derivative of the weakly singular integral kernel is strongly singular or in the sense of Hadamard finite part integral if the the directional derivative of the weakly singular integral kernel is hypersingular. We will try to discover how the strongly singular and hypersingular integrals are generated and propose… More >

H.B. Chen

The International Conference on Computational & Experimental Engineering and Sciences, Vol.20, No.4, pp. 105-106, 2011, DOI:10.3970/icces.2011.020.105

Abstract The calculation of field values and their derivatives near the domain boundary through the boundary element method (BEM) will meet the nearly singularity problem, i.e. the boundary layer effect problem. The tangential derivatives of field values on the boundary often meet an obvious deduction of calculation accuracy. An effective algorithm was proposed by Chen et al. [1,2] to treat these two problems in the same time in elastic BEM and it was recently extended to calculate the second derivative values in potential problem [3]. This algorithm is based on the regularized formulations and is now More >