Finite element (FE) modeling is becoming increasingly prevalent in the world of cardiac mechanics; however, many existing FE models are phenomenological and thus do not capture cellular-level mechanics. This work implements a cellular-level contraction scheme into an existing nonlinear FE code to model ventricular contraction. Specifically, this contraction model incorporates three myosin states: OFF-, ON-, and an attached force-generating state. It has been speculated that force-dependent transitions from the OFF- to ON-state may contribute to length-dependent activation at the cellular level. The current work investigates the contribution of force-dependent recruitment out of the OFF-state to ventricular-level function, specifically the Frank-Starling relationship, as seen through the end-systolic pressure-volume relationship (ESPVR). Five FE models were constructed using geometries of rat left ventricles obtained via cardiac magnetic resonance imaging. FE simulations were conducted to optimize parameters for the cellular contraction model such that the differences between FE predicted ventricular pressures for the models and experimentally measured pressures were minimized. The models were further validated by comparing FE predicted end-systolic strain to experimentally measured strain. Simulations mimicking vena cava occlusion generated descending pressure volume loops from which ESPVRs were calculated. In simulations with the inclusion of the OFF-state, using a force-dependent transition to the ON-state, the ESPVR calculated was steeper than in simulations excluding the OFF-state. Furthermore, the ESPVR was also steeper when compared to models that included the OFF-state without a force-dependent transition. This suggests that the force-dependent recruitment of thick filament heads from the OFF-state at the cellular level contributes to the Frank-Starling relationship observed at the organ level.