The instantaneous regional stress distribution within the myocardium, which cannot be directly measured, has been estimated using improved numerical methods and nonaxisymmetric biventricular geometry. To do this, we have employed computer-aided solid mathematical modeling to generate a three-dimensional representation for an ex vivo canine biventricular unit using magnetic resonance imaging. A two-dimensional transverse section was isolated from the solid mathematical model for regional stress analysis using p-version finite element analysis. Loading conditions and material property descriptions were taken from published reports. Analyses showed the maximum principal stresses to range from -1.76 × 105 to 8.52 × 105 dynes/cm2 during systolic loading, and from -3.85 × 104 to 1.13 × 105 dynes/cm2 during diastolic loading. This study demonstrates that magnetic resonance image-based solid mathematical biventricular models are suitable for regional stress analysis using p-version finite element analysis, p-Version finite element analysis using magnetic resonance image-based cardiac representations facilitates in vivo stress-strain analyses and may allow the clinical estimation of regional myocardial stress.