首先，利用标定好的接触模型参数在DEM程序中生成四种不同尺寸（2μm, 10μm, 20μm, 100μm）的W粉的初始堆积结构。由于受重力作用，粉体内部的相对密度分布呈现分层趋势，并将密度分布情况作为初始条件导入FEM程序中。研究粉体整体在冷等静压压制过程的宏观性质相对密度的变化。

First, use the calibrated contact model parameters to generate the initial accumulation structure of W powder of four different sizes (2μm, 10μm, 20μm, 100μm) in the DEM program. Due to gravity, the relative density distribution inside the powder presents a stratification trend, and the density distribution is imported into the FEM program as an initial condition. The changes in the relative density of the macroscopic properties of the powder as a whole during the cold isostatic pressing process are studied.