This paper evaluates the degree of inelasticity in hydraulic fracture propagation in weak sandstone by using models obtained by the Discrete Element Method as a basis for quantifying Failure Assessment Diagram (FAD) and J-integral. Although current engineering models use Linear Elastic Fracture Mechanics (LEFM), including improvements like fracture tip damage zone where fracture tip plasticity is small enough to fit into the near-field stress zone, this paper hypothesizes that hydraulic fracture propagates in an inelastic regime in certain conditions typically for underground rock reservoirs. Discrepancies within field or laboratory data, including microcrack clouds identified as acoustic emissions that induce plastic deformations, motivate the investigation of inelasticity and implementation of elastoplastic fracture propagation theories. Discrete Element Method (DEM) has been widely used to approach rock failure problems from the micromechanical level, where explicit modeling of local particle-bond breakage enables insights into propagating fracture branching, damage, microcrack coalescence, stress–strain re-distribution, irreversible deformation, and fracture arrest. Results show that inelastic J-integral increases dramatically with rock confinement, especially its non-elastic portion, and the elastic portion of the total J-integral remains a relatively small part. Higher confinement stresses and higher confinement stress contrast enhance the inelastic fracture propagation. Additionally, shorter fracture lengths and smaller dry tip zones characterize inelasticity. Rock stiffness increase leads to total J-integral increase and decrease of J-elastic, leading to pronounced inelasticity. Therefore, results indicate that LEFM is rarely applicable for describing fracture propagation through rock at higher confinement stresses and stiffness.