As shown in a recent study, estimation of the magnitude of the acoustic radiation force (ARF) in isotropic tissue-like solids is the key to measuring a third-order nonlinear modulus of elasticity (C)—responsible for coupling the shear and volumetric constitutive responses of a material. Due to experimental limitations, a direct measurement of the ARF in tissue is however difficult. To overcome this problem, we deploy 3D elastodynamic finite element (EFE) simulations of the ARF experiment in tissue-mimicking phantoms to estimate (i) the “background” shear modulus from the phase of the induced shear waves; (ii) the ARF magnitude from the shear wave amplitude, and (iii) the modulus C from the knowledge of the ARF. In this way, by interpreting the ARF-generated shear waves through the prism of 3D-EFE simulations, we facilitate local estimation of the nonlinear C modulus with spatial resolution equaling the size of the focal region. We call this method the C-Elastography (CE). In this study, the CE technique is applied to gelatin phantoms containing an agar-based inclusion. The profiles of the C-modulus in phantoms tested demonstrate good correlation with the phantom geometry, showing marked (and sharp) C-contrast at push points acting on the inclusion. [Acknowledgment: NIH-Grant EB23113.]