• Publisher-DOI

In this work, a very efficient numerical solution of three‐dimensional boundary value problems of linear elasticity including stress singularities is discussed, focussing on its convergence behaviour. For the employed scaled boundary finite element method, a discretization is only needed at the boundary, while the solution is considered analytically in a scaling coordinate. This presents a major advantage for two‐dimensional problems, when the scaling center is placed at a stress singularity. Unfortunately, three‐dimensional problems usually do not only include point singularities but also line singularities, which results in singular gradients in the boundary coordinates and thereby diminishes the method's original advantages. To alleviate this drawback, this work discusses an enrichment of the separation of variables representation with analytical asymptotic near fields of the line singularities. In contrast to previous works, besides the near‐field functions with λ=0.5, also those with λ=1.5 were determined and used for enrichment. This leads to a high accuracy and it is shown that this approach is required to recover the convergence properties of smooth boundary value problems without singularities when using quadratic Lagrange shape functions. In order to recover the convergence rates for higher order shape functions, near‐field functions with higher singularity exponent have to be included for enrichment.