We have developed an approximate method for simulating the three-dimensional contact of soft biphasic tissues in diarthrodial joints under physiological loading. Input to the method includes: (i) kinematic information describing an in vitro joint articulation, measured while the cartilage is deformed under physiological loads, (ii) geometric properties for the relaxed (undeformed) cartilage lavers, obtained for the analyses in this study via stereophotogrammetry, and (iii) material parameters for the biphasic constitutive relations used to represent cartilage. Solid models of the relaxed tissue layers are assembled in physiological positions, resulting in a mathematical overlap of the cartilage layers. The overlap distribution is quantified and converted via the biphasic governing equations into applied traction boundary conditions for both the solid and fluid phases for each of the contacting layers. Linear biphasic, three-dimensional, finite element analysis is performed using the contact boundary conditions derived for each of the contacting layers. The method is found to produce results consistent with the continuity, requirements of biphasic contact. Comparison with results from independent, biphasic contact analyses of axisymmetric problems shows that the method slightly underestimates the contact area, leading to an overestimation of the total traction, but yields a good approximation to elastic stress and solid phase displacement.