AbstractIt is well known that there are many benefits in solving axisymmetric solids subjected to non‐axisymmetric loading by using Fourier series expansion of field variables. Although many researchers have utilized this technique to solve problems of cooling towers and rocket nozzles, etc. by the Finite Element Method, the BEM formulation and implementation of such a method have been available only recently, mainly due to the efforts of Rizzo and Shippy3 and Mayr et al.4 In this paper these early developments have been generalized further by including more expansion terms, an efficient decomposition of kernels into symmetric and antisymmetric forms, improved evaluation of matrix coefficients and developing more general computer implementation to include surface stress calculations, effects of body forces and extensive substructuring facilities to analyse multi‐zone problems.