No abstract

No abstract

a b s t r a c tIt is well known that the boundary element method may induce spurious eigenvalues while solving eigenvalue problems. The finding that spurious eigenvalues depend on the geometry of inner boundary and the approach utilized has been revealed analytically and numerically in the literature. However, all the related efforts were focused on eigenproblems involving circular boundaries. On the other hand, the extension to elliptical boundaries seems not straightforward and lacks of attention. Accordingly, this paper performs an analytical investigation of spurious eigenvalues for a confocal elliptical membrane using boundary integral equation methods (BIEM) in conjunction with separable kernels and eigenfunction expansion. To analytically study this eigenproblem, the elliptic coordinates and Mathieu functions are adopted. The fundamental solution is expanded into the separable kernel by using the elliptic coordinates and the boundary densities are expanded by using the eigenfunction expansion. The Jacobian terms may exist in the separable kernel, boundary density and boundary contour integration and they can cancel each other out. Therefore, the orthogonal relations are reserved in the boundary contour integration. In this way, a similar finding about the mechanism of spurious eigenvalues is found and agrees with those corresponding to the annular case. To verify this finding, the boundary element method and the commercial finite-element code ABAQUS are also utilized to provide eigensolutions, respectively, for comparisons. Good agreement is observed from comparisons. Based on the adaptive observer system, the present approach can deal with eigenproblems containing circular and elliptical boundaries at the same time in a semi-analytical manner. By using the BIEM, it is found that spurious eigenvalues are the zeros of the modified Mathieu functions which depend on the inner elliptical boundary and the integral formulation. Finally, several methods including the CHIEF method, the SVD updating technique and the Burton & Miller method are employed to filter out the spurious eigenvalues, respectively. In addition, the efficiency of the CHIEF method is better than those of the SVD updating technique and the Burton & Miller approach, since not only hypersingularity is avoided but also computation effort is saved.

Following the success of using the null-field integral approach to determine the torsional rigidity of a circular bar with circular inhomogeneities (Chen and Lee in Comput Mech 44(2):221-232, 2009), an extension work to an elliptic bar containing elliptic inhomogeneities is done in this paper. For fully utilizing the elliptic geometry, the fundamental solutions are expanded into the degenerate form by using the elliptic coordinates. The boundary densities are also expanded by using the Fourier series. It is found that a Jacobian term may exist in the degenerate kernel, boundary density or boundary contour integral and cancel out to each other. Null-field points can be exactly collocated on the real boundary free of facing the principal values using the bump contour approach. After matching the boundary condition, a linear algebraic system is constructed to determine the unknown coefficients. An example of an elliptic bar with two inhomogeneities under the torsion is given to demonstrate the validity of the present approach after comparing with available results.

No abstract