New discrete element models for elastoplastic problems

The discrete element method （DEM） has attractive features for problems with severe damages, but lack of theoretical basis for continua behavior especially for nonlinear behavior has seriously restricted its application. The present study proposes a new approach to developing the DEM as a general and robust technique for modeling the elastoplastic behavior of solid materials. New types of connective links between elements are proposed, the interelement parameters are theoretically determined based on the principle of energy equivalence and a yield criterion and a flow rule for DEM are given for describing nonlinear behavior of materials. Moreover, a numerical scheme, which can be applied to modeling the behavior of a continuum as well as the transformation from a continuum to a discontinuum, is obtained by introducing a fracture criterion and a contact model into the DEM. The elastoplastic stress wave propagations and the tensile failure process of a steel plate are simulated, and the numerical results agree well with those obtained from the finite element method （FEM） and corresponding experiment, and thus the accuracy and efficiency of the DEM scheme are demonstrated.