Recently, a novel two-dimensional spin structure with non-trivial topological properties, called magnetic skyrmion, has been found in many chiral magnets. In most cases, magnetic skyrmions assemble spontaneously and form a lattice structure, called magnetic skyrmion crystal (SkX). SkX, as a novel macroscopic magnetic phase, may interact with different types of external fields through the intrinsic multi-field coupling of the material, resulting in many peculiar physical phenomena. It is found that due to the intrinsic magnetoelastic coupling of chiral magnets, SkX not only influences the mechanical properties of the materials, but also has emergent elastic properties when subjected to external forces. In this review, we first introduce and categorize various types of SkX-related magnetoelastic phenomena, and then introduce a unified theoretical framework to analyze these magnetoelastic phenomena. Specifically, we establish the Landau-Ginzburg free energy functional with a comprehensive description of the magnetoelastic effect for B20 chiral magnets obtained through symmetry analysis, and prove that SkX should be described by a Fourier series due to its wave nature. We show quantitative agreement between theoretical results and experimental results for three types of phenomena:1) the temperature-magnetic field phase diagrams of MnSi suffering uniaxial compression, it is found that uniaxial compression in the direction0, 0, 1T constricts the stable region of the skyrmion phase in the phase diagram, while uniaxial compression in the direction1, 1, 0T extends the stable region of the skyrmion phase in the phase diagram; 2) the emergent elastic behavior of SkX, it is found that this property derives from the magnetoelastic effect of the underlying material, and the linear constitutive equation (with coefficient matrix ) which determines the emergent deformation of SkX, is briefly introduced; 3) the variations of elastic coefficients C11, C33, C44, and C66 with the external magnetic field for MnSi, and the predictions of the variation of C12 and C13 are provided by the theory. Based on the theoretical framework, the analytical solutions of the eigenstrain problems for chiral magnets hosting SkX and the surface configuration of SkX in a half-space magnet are introduced. In this process, we show how to use the theoretical framework to deal with different problems. Finally, we make a summary and suggest several directions for the future development of this field.