The invention provides a numerical simulation method for generating a three-dimensional damaged pattern. The method comprises the steps of establishing a three-dimensional geometric model for an object described in a three-dimensional damaged pattern, dissecting mesoscopic grids, setting a boundary condition, an initial condition and material parameters for the numerical simulation, selecting a yield surface with angular points, considering a strength-softening constitutive relationship and a non-associated flow rule, conducting the numerical simulation on the geometric model based on the finite-element dynamic relaxation method to get a solution, and outputting a three-dimensional damaged pattern. In this way, the three-dimensional damaged pattern can be generated directly based on a mechanical equation without the need of pre-assuming any physical mechanism. The above method provided in the invention is used for directly generating three-dimensional damaged patterns. Three-dimensional damaged patterns, generated based on the above method, can better reflect the actual damage situations of materials. Therefore, the mechanical performances of materials can be reflected more intuitively and more truly.