The finite element method using the principle of virtual work was applied to the biphasic theory to establish a numerical routine for analyses of articular cartilage behavior. The matrix equations that resulted contained displacements of the solid matrix (u) and true fluid pressure (p) as the unknown variables at the element nodes. Both small and large strain conditions were considered. The algorithms and computer code for the analysis of two-dimensional plane strain, plane stress, and axially symmetric cases were developed. The u-p finite element numerical procedure demonstrated excellent agreement with available closed-form and numerical solutions for the configurations of confined compression and unconfined compression under small strains, and for confined compression under large strains. The model was also used to examine the behavior of a repaired articular surface. The differences in material properties between the repair tissue and normal cartilage resulted in significant deformation gradients across the repair interface as well as increased fluid efflux from the tissue.