Mathematical models of cardiac function at the cellular level include three major components, such as electrical activity, Ca2+ dynamics, and cellular shortening. We developed a model for mouse ventricular myocyte contraction which is based on our previously published comprehensive models of action potential and Ca2+ handling mechanisms. The model was verified with extensive experimental data on mouse myocyte contraction at room temperature. In the model, we implemented variable sarcomere length and indirect modulation of the tropomyosin transition rates by Ca2+ and troponin. The resulting model described well steady-state force-calcium relationships, dependence of the contraction force on the sarcomere length, time course of the contraction force and myocyte shortening, frequency dependence of the contraction force and cellular contraction, and experimentally measured derivatives of the myocyte length variation. We emphasized the importance of the inclusion of variable sarcomere length into a model for ventricular myocyte contraction. Differences in contraction force and cell shortening for epicardial and endocardial ventricular myocytes were investigated. Model applicability for the experimental studies and model limitations were discussed.