A numerical method for elasticity in axisymmetric semi-infinite domains is proposed.
The method couples finite elements with Dirichlet-to-Neumann boundary conditions.
The lack of an explicit closed-form expression for the DtN map needs to be overcome.
This is done by using a procedure that combines analytical and numerical techniques.
The numerical experiments confirm the effectiveness and accuracy of the method.