Traditional Euler angle-based methods for the kinematic and dynamic modeling of spherical joints involve highly complicated formulas that are numerically sensitive, with complex bookkeeping near local coordinate singularities. In this regard, exponential coordinates are known to possess several advantages over Euler angle representations. This paper presents several new exponential coordinate-based formulas and computational procedures that are particularly useful in the modeling of mechanisms containing spherical joints. Computationally robust procedures are derived for evaluating the forward and inverse formulas for the angular velocity and angular acceleration in terms of exponential coordinates. We show that these formulas simplify the parametrization of joint range limits for spherical joints, and lead to more compact equations in the forward and inverse dynamic analysis of mechanisms containing spherical joints.