The paper discusses an iterative scheme for solving large-scale three-dimen
sional linear elasticity problems, discretized on a tensor product of two-d
imensional and one-dimensional meshes. A framework is chosen of the additiv
e AMLI method to develop a preconditioner of a 'black-box' type which is ro
bust with respect to discontinuities of the problem coefficients and impose
s only weak (and acceptable in practice) restrictions on the choice of the
meshing procedure. The preconditioner works on a hierarchical sequence of n
ested finite element spaces to solve the problem with arithmetic cost, near
ly proportional to the number of degrees of freedom on the finest mesh. It
is particularly well suited for the case when the solution is known to be s
trongly varying in certain subregions of the domain and the mesh is locally
prerefined there to reduce the discretization error. Copyright (C) 1999 Jo
hn Wiley & Sons, Ltd.