Abstract: The compression molding process of salix psammophila granules was simulated based on finite element and discrete element numerical methods. Furthermore,the forming mechanism of salix psammophila granules was revealed by analyzing the macro and meso mechanical properties at the initial,middle,middle,late and the end stages. The results show that macroscopically the stress of the material column increases gradually from zero,and in the beginning the stress concentrates on the upper surface,then gradually spreads to the upper part,and finally the stress of the whole material column is evenly distributed;Microscopically,the particle distribution changes from loose to dense,the distribution of normal and tangential contact forces changes form the initial disorder to the formation of rhomboid,and then to the large ring,and finally to the small triangular stable structure. The particle velocity field goes straight down and forms a separation zone,which moves down until it disappears.