**Abstract**

Based on the concept of embedded discrete elements (EDEs), the discrete elementembedded finite element model (DEFEM) is extended in this work. The new method can be used to calculate the motion and stress variation of particles. This work discusses its application in granular flow simulation for particle motions with small deformations. The updated Lagrangian finite element method is used to obtain the coupling solution of the internal stress and the overall motion of particles in the DEFEM. The computation of deformation displacement is based on the concepts of displacement decomposition (translational and rotational motions and deformation displacement). The deformation displacement is the difference between particles and template particles [rigid body, using the discrete element method (DEM) to calculate translational and rotational displacements]. It is used to calculate the dynamic stress distribution of particles and the internal force of the node. Therefore, it has a wide scope of application (for example, it can be extended to non-spherical particles). The software validation proves the accuracy of this method. The application of the DEFEM in the accumulation process of particles is given. The motion characteristics and deformation of particles are discussed, and the stress distribution and force chain structure in particle accumulation are obtained.

**Keywords:** discrete element method, finite element method, dynamic stress distribution, pebble bed, pebble flow, coupled method
