We review recent developments in the field of chemodynamical simulations ofelliptical galaxies, highlighting (in an admittedly biased fashion) the work conducted with our cosmological N-body/SPH code GCD+. We have demonstrated previously therecovery of several primary integratedearly-type system scalingrelations (e.g. colour-magnitude relation, LX − TX − [Fe/H]X) when employing a phenomenological AGN heating scheme in conjunctionwith a self-consistent treatment of star formation, supernovae feedback, radiative cooling, chemical enrichment, and stellar/X-raypopulation synthesis. Here we emphasise characteristics derivedfrom the full spatial information contained within the simulated dataset, including stellar and coronal morphologies, metallicity distributionfunctions, and abundance gradients.