We develop a magneto-elastic (ME) coupling model for the interaction between the vortex lattice and crystal elasticity. The anisotropies in superconductivity and elasticity are simultaneously included in the GL theory frame. Under this consideration, the expression of the free energy unifies the different forms of the classical results. Concerning the ME effect on the magnetization, the theory can give a satisfying description for the field dependence of magnetization near the upper critical field. The contribution of the ME interaction to the magnetization is comparable to the vortex-lattice energy, in materials with relatively strong pressure dependence of the critical temperature. While the magnetization components along different vortex frame axes are strain dependent, the magnetization ratio is independent of the ME interaction. It is stressed that the GL description of the magnetization ratio is applicable only if the applied field moderately close to the upper critical field.