It is of great significance to study the damage and failure law of hydrate-bearing sediments for for exploration and development, as well as for warning secondary disasters such as tsunami and earthquake. The discrete element modeling and simulation method has the advantages of low cost, strong repeatability and accurate response to the microstructure of samples, therefore , the discrete element method is used to simulate and analyze the damage and failure law of hydrate-bearing sediments in this paper. First, a triaxial undrained shear teat model of hydrate-bearing sediments is established. by discrete element simulation software; Then, the effects of different influencing factors on the fracture characteristics of hydrate hydrate-bearing sediments is studied; Finally, the effects of different factors on the initiation stress and damage stress of hydrate sediments are analyzed, and the damage law of hydrate-bearing sediments is obtained. The results show that: (1) The cementation of hydrate particles is greater than its bearing effect, thus the lithology of sediments changes from loose sand body to brittle rock with the increase of hydrate saturation, resulting in the corresponding change of fracture characteristics from loose sand body failure to brittle rock failure. (2) With the increase of hydrate saturation, the initiation stress of sediment shows a step-by-step increase law, and the damage stress is positively correlated with the shear modulus. (3) The heterogeneity of hydrate distribution is closely related to the failure mode of sediments. With the increase of hydrate distribution heterogeneity, the initiation stress displays an increase-decrease-increase pattern while the damage stress shows a law of increasing and then decreasing, which has the high sensitivity to heterogeneity of hydrate distribution. The research results have certain theoretical reference significance for the exploitation practice of natural gas hydrate.