Inferring the mechanical properties of soft tissues from measured deformations is a fundamental challenge in elastography. A rarely examined assumption underlying existing approaches is that the assumed constitutive law correctly describes the imaged material. When it fails, inversion still yields plausible-looking estimates - an illusion of fit with no indication of local model invalidity, which can mislead clinical interpretation. We propose a probabilistic framework that transforms constitutive model validity from an implicit assumption into an explicit, spatially resolved inference target. The key is to treat the stress field as an independent latent variable rather than deriving it from the constitutive law. This enables a pointwise comparison between the stress required by mechanical equilibrium and the stress predicted by the assumed constitutive model. Both governing equations enter the probabilistic learning objective as virtual observables with separate precision hyperparameters: the conservation law precision is set a priori to a small value reflecting its undisputed validity, while the constitutive precision is inferred under a sparsity-promoting prior. The resulting constitutive precision field provides an uncertainty-aware map of where the assumed model is supported by the data and where it is not. Inference is carried out via stochastic variational inference and is forward-model-free. We validate the framework on synthetic harmonic elastography experiments on a brain-slice geometry with an anisotropic inclusion. The inferred precision field identifies the inclusion with a five-order-of-magnitude precision contrast against the valid domain, robustly across 25-35 dB noise and four-fold sparser observations. A phantom experiment with ultrasound measurements on a linear elastic material yields no false-positive violations and recovers the true stiffness contrast.