A contacting constitutive model was established by using the particle flow theory and a numerical simulation of biaxial test on GBS (graded broken stone) was proposed based on PFC2D. The reliability of the numerical method was verified by dynamic triaxial laboratory test. The plastic deformation characteristics of GBS were studied by simulation. The destruction of critical stress and that of critical strain of GBS were studied by the equation of accumulative plastic deformation. The results show that numerical test results of plastic deformation law of GBS meet the actual shakedown law. The rational damage critical strain of GBS is 2.5%. Under this critical strain, the relationship between critical stress and confining pressure is proportional with a proportional coefficient of 4.95. The increasing of confining pressure has a great effect on the resistance performance of plastic deformation of graded broken stone.