Efficient computation of order and mode of corner singularities in 3D-elasticity

A general numerical procedure is presented for the e#cient computation of corner singularities, which appear in the case of non-smooth domains in three-dimensional linear elasticity. For obtaining the order and mode of singularity, a neighbourhood of the singular point is considered with only local boundary conditions. The weak formulation of the problem is approximated by a Galerkin–Petrov =nite element method. A quadratic eigenvalue problem (P + Q+ 2R) u=0 is obtained, with explicitly analytically de=ned matrices P;Q;R. Moreover, the three matrices are found to have optimal structure, so that P;R aresymme tric and Q is skew symmetric, which can serve as an advantage in the following solution process. On this foundation a powerful iterative solution technique based on the Arnoldi method is submitted. For not too large systems this technique needs only one direct factorization of the banded matrix P for =nding all eigenvalues in the interval Re( ) ∈ (−0:5; 1:0) (no eigenpairs can be ‘lost’) as well as the corresponding eigenvectors, which is a great improvement in comparison with the normally used determinant method. For large systems a variant of the algorithm with an incomplete factorization of P is implemented to avoid the appearance of too much =ll-in. To illustrate the eDectiveness of the present method several new numerical results are presented. In general, they show the dependence of the singular exponent on diDerent geometrical parameters and the material properties