The application of the finite element method to the modeling of Cosserat solids is investigated in detail. In two- and three-dimensional elasticity problems, the nodal unknowns are the components of the displacement vector, which form a linear field. In contrast, when dealing with Cosserat solids, the nodal unknowns form the special Euclidean group SE(3), a nonlinear manifold. This observation has numerous implications on the implementation of the finite element method and raises numerous questions: (1) What is the most suitable representation of this nonlinear manifold? (2) How is it interpolated over one element? (3) How is the associated strain field interpolated? (4) What is the most efficient way to obtain the discrete equations of motion? All these questions are, of course intertwined. This paper shows that reliable schemes are available for the interpolation of the motion and curvature fields. The interpolated fields depend on relative nodal motions only, and hence, are both objective and tensorial. Because these schemes depend on relative nodal motions only, only local parameterization is required, thereby avoiding the occurrence of singularities. For Cosserat solids, it is preferable to perform the discretization operation first, followed by the variation operation. This approach leads to considerable computation efficiency and simplicity.