A new method for calculating the direction of fracture propagation by stress numerical search based on the displacement discontinuity method