The flow mechanism of shale gas reservoir is complex including viscous flow, adsorption-desorption, and Kundsen diffusion. The shale reservoir is usually developed by hydraulic fracturing which produces large-scale fractures. The porosity and permeability of the fractures are greatly affected by the geostress, so that the influence of the hydro-mechanical coupling needs to be taken into account in the numerical simulation. The reservoir geological model is usually described by corner-point grid and due to the existence of unmatched nodes in the corner-grid and complex geological structure which makes the numerical simulation of the hydro-mechanical coupling in the corner-point grid not be solved effectively. In this paper, based on corner-point grids, an effective hydro-mechanical coupling numerical simulation method for fractured shale gas reservoirs is proposed, in which the flow equations and elastic equations are discretized by the finite volume method and the virtual element method respectively, the fracture is characterized by embedded discrete fracture model and fixed-stress splitting algorithm is used to solve the flow and stress equations itetatively. The correctness of this method is verified and the effects of different parameters on gas production in fractured shale gas reservoirs are analyzed.