A computational model which integrates a complex circulatory model and a finite element to determine the dynamics of a left ventricle continuously over consecutive cardiac cycles. The model includes determining a LV pressure (plv) using a circulatory model, using a finite element model plv as input and determining a LV volume (vlv), computing a LV elastance according to: elv=plv/vlv, driving the circulatory model with the elv; and returning to determining a LV pressure and starting the next iteration, wherein the steps continue at a sufficient time resolution for a desired number of entire cardiac cycles. The dynamic Young's modulus functions are assigned to individual finite elements, resulting in a time-varying left ventricular elastance that drives the circulatory model.