The application relates to the technical field of
engineering design, in particular to a non-
Darcy seepage stress
coupling simulation method based on a fracture and matrix grid, information such as input fracture and
matrix unit grids, boundary conditions, material parameters and
control parameters is input; initial values and initial stresses of basic unknown quantities are given; at each
time step, matrix seepage and embedded fracture non-
Darcy flow iterative solving and matrix stress
balance equation solving are sequentially carried out; whether convergence is judged, if the
time step calculation does not converge, the next step iterative solving is continuously carried out; after the convergence condition is reached, the next
time step solving is entered until the time ends. The application provides a high-precision and high-efficiency non-
Darcy seepage stress
coupling simulation method, and improves the problems of great difficulty in modeling of
local structure surfaces, low efficiency and poor effect in traditional finite element solving, and can effectively simulate non-Darcy seepage characteristics and mechanical deformation behaviors of rock
mass under high seepage pressure.