We simulate the chemodynamical evolution of a hundred elliptical galaxies using our GRAPE-SPH code, and succeed in reproducing the radial metallicity gradients and the global scaling relations such as the fundamental plane. These observations cannot be explained by either monolithic collapse or by major merger alone. Rather it requires a model in which both formation processes arise, such as the present CDM scheme.