In this paper, we develop two least-squares approaches for the solution of
the Stokes equations perturbed by a Laplacian term. (Such perturbed Stokes
equations arise from finite element approximations of the Reissner Mindlin
plate.) Both are two-stage algorithms that solve rst for the curls of the r
otation of the fibers and the solenoidal part of the shear strain, then for
the rotation itself (if desired). One approach uses L-2 norms and the othe
r approach uses H-1 norms to de ne the least-squares functionals. It is sho
wn that the H-1 norm approach, under general assumptions, and the L-2 norm
approach, under certain H-2 regularity assumptions, admit optimal performan
ce for standard finite element discretization and either standard multigrid
solution methods or preconditioners. These methods do not degrade when the
perturbed parameter (the plate thickness) approaches zero. We also develop
a three-stage least-squares method for the Reissner Mindlin plate, which r
st solves for the curls of the rotation and the shear strain, next for the
rotation itself, and then for the transverse displacement.