Understanding the spatiotemporal distribution characteristics of cryoprotectant concentration in the cartilage tissue during loading is essential for the rational design of the loading procedure and the successful cryopreservation of articular cartilage (AC).To construct a diffusion mass transfer model to describe the permeation of cryoprotectant into AC.The model was based on the binary diffusion thermodynamic model to correlate the effective diffusion coefficient with the infinite dilution diffusion coefficient. The activity coefficient was estimated using the UNIFAC model. The concentration dependence of the diffusion coefficient was expressed using the Vignes formula, and the infinite dilution diffusion coefficient was calculated using the Siddiqi-Lucas correlation.For dimethyl sulfoxide (Me2SO), glycerol (GLY), ethylene glycol (EG) and propylene glycol (PG), the coefficient of determination and mean relative error of the average cryoprotectant concentration in the cartilage tissue calculated by the model and experimental values reported in the literature are 0.959-0.998 and 0.98-15.19%, respectively, for Me2SO; 0.990-0.995 and 13.56-19.19% for GLY; 0.969-0.988 and 8.89-22.09% for EG; and 0.971-0.992 and 5.35-23.76% for PG.The developed model could be used for these four widely used cryoprotectants to directly predict the spatiotemporal distribution of cryoprotectant concentration in the cartilage tissue during loading.