Multivariable dynamic ankle mechanical impedance in two coupled degrees-of-freedom (DOFs) was quantified when muscles were active. Measurements were performed at five different target activation levels of tibialis anterior and soleus, from 10% to 30% of maximum voluntary contraction (MVC) with increments of 5% MVC. Interestingly, several ankle behaviors characterized in our previous study of the relaxed ankle were observed with muscles active: ankle mechanical impedance in joint coordinates showed responses largely consistent with a second-order system consisting of inertia, viscosity, and stiffness; stiffness was greater in the sagittal plane than in the frontal plane at all activation conditions for all subjects; and the coupling between dorsiflexion–plantarflexion and inversion–eversion was small—the two DOF measurements were well explained by a strictly diagonal impedance matrix. In general, ankle stiffness increased linearly with muscle activation in all directions in the 2-D space formed by the sagittal and frontal planes, but more in the sagittal than in the frontal plane, resulting in an accentuated “peanut shape.” This characterization of young healthy subjects' ankle mechanical impedance with active muscles will serve as a baseline to investigate pathophysiological ankle behaviors of biomechanically and/or neurologically impaired patients.