This paper proposes a novel mean-field matrix-analytic method in the study of bike sharing systems, in which a Markovian environment is constructed to express time-inhomogeneity and asymmetry of processes that customers rent and return bikes. To achieve effective computability of this mean-field method, this study provides a unified framework through the following three basic steps.