Gas recovery from hydrate-bearing sediments attracts worldwide interest due to the huge reserves of methane there, recognized as a promising future energy resource. However, the mechanical behavior of hydrate-bearing sediment is rarely studied. In this study, the drained shear test on hydrate-bearing sand with different hydrate saturation was simulated using the discrete element method. In particular, the pore habits of hydrate in sands were considered in the simulation, and the meso-failure behavior such as the evolution of contact-force chains and the occurrence of microcracks were thoroughly investigated. In general, the simulations yielded the mechanical response (e.g., deviatoric stress-strain curve) of hydrate-bearing sand similar to the laboratory experiments, and the existence of hydrates plays a crucial role in the failure behavior. It is found that “X-type” shear bands are more likely to form in hydrate-bearing sand with higher hydrate saturation because enough bonds are broken to facilitate microcrack connection. Furthermore, tensile failure, tensile-shear failure, and compression-shear failure exist simultaneously during shearing, of which tensile failure and tensile-shear failure are dominant.