Present in both plutonic and volcanic rocks, enclaves are inclusions of magma into a compositionally distinct magmatic host. Classical links between their shapes and the dynamical conditions that prevailed during their formation have been drawn from fluid-fluid analogies and from solid rock mechanics. Magmas, however, are hydrogranular suspensions with a rheology distinct from these two-end-members. This work presents results from computational fluid dynamics with discrete element modeling (CFD-DEM) aimed at deforming crystal-rich enclaves in pure shear. The CFD-DEM approach explicitly resolves solid-solid interactions such as contact and friction while taking into account fluid coupling. The first series of deformation involved only pure fluids to validate the setup. The second series comprised seven runs aimed at reproducing magmatic conditions. Enclaves were made of a cylindrical suspension of particles embedded into a host with different characteristics. In both media, particles and fluids had densities, viscosities, elastic characteristics, and sizes tailored to the geological constraints of the Adamello batholith, Italy. Each run corresponds to a temperature along the two respective crystallization paths and span crystal contents from 10 to 62 vol.%. Results show that, to first order, deformation does not depend on differences in melt viscosities, crystal contents, or bulk viscosity contrast. This is due to the formation of force chains parallel to the main compression direction, which transmits stress across the enclave. A simple, first-order relationship could be fitted to our data to relate shear and enclave deformation, which we applied to the case of the Adamello pluton. There is a second-order dependence of deformation on the onset of particles contacts and force chains, which are both related to particle concentration. The main control of these second-order effects lies in the host crystal content. Enclave particles pack early, quickly erasing differences in initial content and building force chains parallel to the compression axis that transmit stresses to the host. Whether the host is able to transmit those stresses across its own volume is controlled by host crystal content.