Compacting relative density and force chain analysis of powders with different particle size ratios based on discrete element
-
摘要: 利用PFC三维数值模拟软件，通过改变粉末颗粒粒径分布建立各组冷压模型，得到压制过程中相对密度变化规律与力链分布情况。在特定粉末粒径配比下，能够得到相对密度最高的压坯。结果表明：在大、中、小粒径颗粒质量比为60:15:25的粒径配比下，压坯相对密度最高，压坯相对密度并不会随着细粉比例不断增加而一直提高；在压制过程中，随着附加细粉占比上升，压制方向上能产生更大应变。侧压系数与泊松比受粉末粒径分布影响较小，且在压制后期，在压坯已获得较高相对密度的情况下，会因缺乏足够的驱动力与位移空间发生下降；混合粒径粉末试样的力链数量远大于单一粒径粉末试样，在强力链数目充足的前提下，结合大量的弱力链能获得更高的压坯相对密度。Abstract: Using the PFC three-dimensional numerical simulation software, the law of relative density change and force chain distribution during the pressing process were obtained by changing the particle size distribution of powders to establish the cold pressing models. Under the specific particle size ratio of powders, the compacts with the highest relative density could be obtained. In the results, the large, medium, and fine particles in the mass ratio of 60:15:25 show the highest compact relative density, and the compact relative density does not increase with the increase of the fine powder ratio. As the proportion of the additional fine powders rises during the pressing process, the greater strain can be generated in the pressing direction. The lateral pressure coefficient and Poisson's ratio are less affected by the particle size distribution of powders, which can decline due to the lack of sufficient driving force and displacement space in the case that the compacts have obtained the higher relative density in the later stage of the compaction. The number of force chains for the powders with the mixed particle size is much greater than that of the powders in the single particle size. In the case that the number of strong force chains is sufficient, the higher compact relative density can be obtained, combining a large number of weak force chains.
-
金属粉末的粒径以及粒径组成不同时，就会表现出不同的压制行为。通常粒径小的颗粒粉末成形性能比较好，但是压缩性能比较差，这就会对压坯的相对密度造成影响。相对密度是粉末冶金制品性能的关键指标，主要受金属粉末材料性能和压制技术的影响，高的相对密度一直是行业追求的目标。生产实践表明，非单一粒径组成的粉末压制性能较好[1‒3]，因此针对粉末粒径分布对粉末压制成品性能的影响开展研究是非常有必要的，尤其是开展粉末粒径分布对粉末压制成品相对密度的影响研究，不仅能探讨粉末性能与粉末压制成品性能之间的机理关系，同时也能促使行业标准的改进，为获得性能更好的高相对密度粉末冶金产品提供一定的行业基础，有利于我国粉末冶金行业的发展。
对非连续体力学行为的研究主要使用离散单元法（discrete element method，DEM）[4‒7]，特别是针对求解大数量且不连续结构的问题。Skrinjar和Larsson[8]采用离散单元法对球形粉末的冷压过程进行了模拟研究，分析了粉末粒径和压坯尺寸比对压坯性能的影响，结果表明，在小颗粒占比高的情况下，尺寸效应是要考虑的关键因素。叶先勇等[9]通过压制多孔过滤材料研究分析了粉末粒径和压制力对孔隙率和渗透率的影响规律，提出相关本构模型。朱鹏程[10]制备了不同粒径的初始氢爆粉，利用他们在完全相同的工艺情况下制备了烧结铷铁硼试样，测取样品的抗弯强度及断裂韧性，分析了粉末粒径差异对磁体力学性能的影响。闫志巧等[11‒12]分别用不同平均粒径的Ti粉和Ti6Al4V合金粉末为原料，研究了粉末粒径对压坯密度、最大压制力和脱模力的影响，发现压坯密度与粉末粒径及松装密度有关，粒径小的粉末更难压制，相对密度也更小，但是烧结密度更高。杨志超[13]通过压制不同粒径分布的粉末研究得到粉末粒径分布对最终产品微观结构的形成具有非常大的影响，细而均匀的粉末所得压坯晶粒也是细而均匀，但是粒径过小的粉末如果工艺不合适，会产生晶粒长大的现象。潘诗琰等[14]利用群体平衡模型定量预测了相邻颗粒聚并效应作用下的粒径分布演化规律，研究了从瞬态到稳态的转变过程中，粒径上限及其变化率，发现颗粒瞬态粗化之后将得到稳态粒径分布。
我国粉末冶金工业已经在产量上领先日本等国家，但高质量产品的供应和推广还有待进步。已经有学者开展了粉末粒径对压制产品性能影响的研究，也有实践表明非单一粒径组成的粉末压制产品性能较好，但还没有研究人员针对这种情况展开细化研究。实际生产中粉末性能与产品性能直接相关，考虑到目前我国粉末原料标准不严格的现状，严重影响了粉末冶金行业对更高相对密度产品的追求。因此，探讨粉末粒径分布对压坯相对密度的影响具有较大意义，而且能为粉末原料行业标准的改进提供参考。为了得到不同粒径分布对压坯相对密度的影响规律，本文按照标定的仿真参数，通过改变粉末粒径分布参数来建立模型并进行冷压模拟。
1. 不同粒径分布的粉末模型建立
1.1 离散元建模分析
1.1.1 理论依据
PFC三维数值模拟软件运行原理是基于两组控制方程，一组是运动方程（牛顿运动方程），另一组是接触力方程（力和位移或位移增量的关系）。在PFC软件中，由接触模型来定义对象互相作用时的力和位移关系，是影响数值宏观力学性质的重要因素之一，球形颗粒接触理论严密，求解过程相当繁琐。因此在不产生显著误差的情况下，有必要对接触模型进行简化处理。采用线性接触模型，即把颗粒间的法向力视为弹簧和阻尼器，切向力视为弹簧、阻尼器和滑动器[13]。等效惯性质量（mc）由式（1）求得。
$$ {m}_{\text{c}}=\left\{\begin{array}{l}\dfrac{{m}^{(1)}{m}^{(2)}}{{m}^{(1)}+{m}^{(2)}},\text{ball-ball}\\ {m}^{(1)},\text{ball-facet}\end{array} \right. $$ (1) 式中：m(1)和m(2)分别表示接触的两个单元的质量，根据判断接触是位于球和球单元之间（球-球，ball-ball）还是球和墙（球-墙，ball-facet）之间来选择计算。
利用便捷的Method方法赋予模型参数，通过设定模型的等效弹性模量（E*）和法向切向刚度比（k*），由软件根据式（2）和式（3）求解。
$$ {k_{\text{n}}} = \frac{{AE^*}}{L} $$ (2) $$ {k_s} = \frac{{{k^n}}}{{{k^*}}} $$ (3) 式中：kn代表接触的法向刚度；A代表接触截面积，计算时取接触者中小的单元半径；L代表接触长度，接触为球-球类型取两个单元的半径和，接触为球-墙类型取球的半径值；ks代表接触的切向刚度。图1所示为接触模型中接触单元关系，把球与球单元的接触视为存在一个圆柱状（虚线框内）的联系，把球与墙的接触视为一个圆柱状的联系。
![]()
图 1 接触模型中接触单元关系：（a）球-球类型；（b）球-墙类型Figure 1. Relationship of the contact unit in the contact model: (a) ball-ball type; (b) ball-facet type1.1.2 动态测量球的设置
为测量颗粒模型内的孔隙率变化，通常会采取在模型内部建立测量球的方式。但随着压制的进行，试样的高度势必会发生变化，如果采用传统固定大小的测量球，只能在初期设置测量球监测直径小于最终压制后的试样高度，但是这样会导致测量范围仅仅位于试样的局部。为解决这一问题，通过编写控制命令，使得测量球的大小随着试样高度降低而逐渐变小，测量球测量的范围始终是在试样高度方向上的最大测量范围。动态测量球的分布方式如下图2所示，为尽可能的覆盖试样范围，采用十字交错的方式排布了9个测量球，测量球的初始直径为圆柱形压制模型的高度，测量球的位置和半径随着压制进行会动态变化，使得其每一个测量球的球心在竖直方向上始终在试样的中心，测量球的直径始终刚好等于压坯的高度。
![]()
图 2 动态测量球分布图：（a）初始测量球俯视图；（b）初始测量球正视图；（c）压制后测量球俯视图；（d）压制后测量球正视图Figure 2. Distribution of the dynamic measuring ball: (a) top view of the initial measuring ball; (b) front view of the initial measuring ball; (c) top view of the measuring ball after pressing; (d) front view of the measuring ball after pressing1.2 模型参数选取与建立
1.2.1 粒径分布标准对比
目前国内外还没有统一的粉末粒径标准，每个企业都有自己的粒径分布定义和表示方式。一直以来，粉末冶金企业都是通过生产试验积累经验，获得发展，缺少理论研究，这也说明了开展粉末粒径分布对压坯相对密度影响的研究是必要的。
查阅国标及国内外具有代表性的几家粉末冶金原料企业对粉末粒径分布的标准，可以发现这些标准主要关注粉末的压缩性，压缩性代表粉末在压制过程中被压紧的能力，在规定的模具和加载条件下，用一定单位压制压力（500 MPa）压制粉末后的压坯密度表示，是粉末冶金成品性能的一项重要指标[1]。通常使用筛网目数（简称目）来表示筛网的孔径和粉末的粒径，所谓目数就是指1英寸（25.4 mm）筛网长度上的网孔数，目数越大，网孔越细，筛孔尺寸与网丝的直径有关，以m代表目数，a代表筛孔尺寸，d代表网丝直径，则其间关系式如式（4）所示[1]。
$$ m = \frac{{25.4}}{{a + d}} $$ (4) 我国工业上对粉末进行筛分通常选用80目、100目、150目、200目、250目、325目等组成一套标准筛，其目数与筛孔尺寸、过筛粒径大小对照如表1[1,14‒15]所示。在目数前添加一个正负号表示粉末颗粒是否能通过该目数的网孔，负数表示粉末颗粒可以通过该目数的网孔，即粒径小于筛网孔径，正数表示粉末颗粒不能通过该目数的网孔，即粒径大于筛网孔径。
表 1 工业标准筛Table 1. Industrial standard screen目数 筛网孔径 / μm 粒径 / μm +80 +180 ≥180 −80～+100 −180～+150 150～180 −100～+150 −150～+106 106～150 −150～+200 −106～+75 75～106 −200～+250 −75～+63 63～75 −250～+325 −63～+45 45～63 −325 −45 ≤45 通过对比相关国内外标准和国内外几家代表性粉末冶金原料企业[16‒18]后发现，压缩性好的粉末，其粒径分布更为具体，尤其是针对目数牌号粉末的主体基础粉和附加细粉这两个范围的粉末，有着具体的标准。将查阅的一些国内外标准[19‒20]中压缩性最高的牌号取出，合并在一起，结果见表2，通过对比发现，压缩性好的粉末，其粒径分布更为具体。相同标准下，同一目数级别的基础粉末比例一般为60%，而对粒径小于45 μm的附加细粉末和过渡区间内的附加中径粉末有不同的标准。
表 2 各标准粒径分布对比Table 2. Comparison of the standard particle size distributions牌号 压缩性 / (g·cm‒3) 粒径分布（质量分数） / % ＞250 μm ＞180 μm ＞150 μm ＞75 μm ＞63 μm ＞45 μm ＜45 μm FHY100·270 ≥6.70 0 0 3 — — — 10～30 300NH ≥7.05 1 15 15 60 5～20 8～23 10～30 HAP100.30H ≥7.10 0 1 10 60 — 20～40 30～45 ABC100.30 ≥7.28# 0 * 10 * * * 20 FSW 100·30H 7.15# 0 1 10 — — — 15～30 注：“—”表示不作具体要求；“*”表示参数保密；压缩性的单位压制压力除特别标注“#”为600 MPa外，其余均为500 MPa。 1.2.2 压制模型参数设置
在离散元数值模拟中，需要通过定义颗粒的微观材料参数来建立与所模拟材料具有相同宏观力学性能的离散元模型。当生成颗粒模型后，需要赋予接触模型，其中接触模型中的参数可以通过离散元理论公式计算获取，这样能够将模型的微观参数和材料宏观参数联系起来，使建立的模型更准确。
首先验证设置的微观材料参数是否能正确反应材料的宏观力学性能。本次模拟的粉末材料为粒径分布不同的颗粒群，颗粒模型如图3所示，按照表3中分布要求，在墙体内部用distribute方式生成了圆形颗粒293319个，颗粒不能透过墙体，数量远超10万个，具有较高的置信度。
表 3 200目纯铁粉试样粉末粒径分布（质量分数）Table 3. Particle size distribution of the pure iron powder samples in 200 mesh% ＞150 μm ＞106 μm ＞75 μm ＞45 μm ＜45 μm 0 5 15 40 40 为了更好的探究不同粒径粉体压坯的压缩性，设置模型中模具尺寸参数为半径2.5 mm，高度为3 mm，由周围一圈墙体作为阴模模壁，顶部墙体和底部墙体分别作为上下模冲，采用铁粉颗粒，其物理性能如表4所示，通过式（2）和式（3）计算得到颗粒间法向刚度与切向刚度比为0.785，颗粒和墙体间法向刚度与切向刚度比为0.196，颗粒间接触模型为线性平行模型，初始孔隙率设置为0.5，在模拟压制前对颗粒进行重排，消除随机生成的颗粒间不平衡力，直到颗粒达到相对稳定的状态。
表 4 铁粉相关参数Table 4. Related parameters of Fe powder材料 密度 /
(kg·m−3)弹性模量 /
GPa法向切向
刚度比摩擦系数 铁粉 7800 190 0.785 0.3 在压制过程中，关注轴向即压制方向上的应力‒应变关系，如图4所示。由图可知，轴向应力随应变增大而增大，增速从平缓到加速上升。在压制初期，由于粉末颗粒间存在较大孔隙，受压制力作用，极易发生大量位移，引起孔隙率快速降低，故轴向应力随应变增大的增速较为平缓；在应力曲线与孔隙率曲线相交后，即应变达到0.16后，应力进入快速增加阶段，在该时间段的颗粒位移已经开始减少，粉末密度达到了一定值，出现了压缩阻力，孔隙率降速度变缓。最终当轴向应变为0.22时，轴向应力达到了370 MPa，实验测得的数据[21]为轴向应变0.22时，轴向应力达到369 MPa，两者一致，说明该接触参数下的压制模型符合实际情况。
1.2.3 模型建立
尝试以100目纯铁粉为研究对象，设定了6组试样的粒径分布，如表5所示，其中N组为单一粒径组成的粉末，全部由106～150 μm范围基础粉末颗粒组成，其余五组混合粒径粉末的粒径构成参照现有标准，保持基础粉末（粒径范围：106～150 μm）质量占比60%不变的前提下，分别改变附加中径粉末（粒径范围：45～106 μm）和附加细粉末（粒径小于45）的占比。
表 5 各组试样粉末粒径分布Table 5. Particle size distribution of the samples试样编号 粒径分布（质量分数） / % ＞150 μm 106～150 μm 45～106 μm ＜45 μm N 0 100 0 0 A 0 60 30 10 B 0 60 25 15 C 0 60 20 20 D 0 60 15 25 E 0 60 10 30 参考胡仙平[21]的研究，选取径高比为1:1.6的模型来探究，便于从泊松比对比中更清楚的体现粒径分布差异对粉体压缩性能的影响，综合考虑后将试样设置为半径3.75 mm、高度为12 mm的圆柱形，接触参数设置与之前一致，通过压制过程模拟分析压制过程中孔隙率与压制力的变化规律。为提高计算速度，每一分级的颗粒粒径都取其最大值，即分别取150 μm、106 μm、45 μm，采用distribute方式生成各组模型，各粒径颗粒数量分布见表6所示。图5为所建各组模型压制前后对比，可以明显看到，由于不同粒径混合导致各组试样的粉末颗粒数量不同。
表 6 各粒径颗粒数目Table 6. Number of the particles in each size试样编号 粒径分布（个数） 150 μm 106 μm 45 μm 总数 N 18076 0 0 18076 A 10822 15626 68923 95371 B 10822 13073 103432 127327 C 10822 10497 137742 159061 D 10822 7897 172301 191020 E 10822 5275 206605 222702 ![]()
图 5 各组粉末压制模型：（a）N组压制前；（b）N组压制后；（c）A组压制前；（d）A组压制后；（e）B组压制前；（f）B组压制后；（g）C组压制前；（h）C组压制后；（i）D组压制前；（j）D组压制后；（k）E组压制前；（l）E组压制后Figure 5. Each group of powder pressing model: (a) before N group suppression; (b) after N group suppression; (c) before A group suppression; (d) after A group suppression; (e) before B group suppression; (f) after B group suppression; (g) before C group suppression; (h) after C group suppression; (i) before D group suppression; (j) after D group suppression; (k) before E group suppression; (l) after E group suppression2. 压制模型结果对比与分析
2.1 粉末孔隙率
压制过程中粉末体内部孔隙率的变化如图6所示，通过孔隙率可得压坯相对密度。由图可知，各组孔隙率都随着压制的进行逐渐下降，但单一粒径组（N组）的孔隙率在整个压制过程中异常高于其余5组混合粒径粉末，各组试样在压制过程中会经历一段孔隙率下降速度减缓的过程，当最终压制力达到500 MPa时，最终孔隙率与缩小比关系如表7所示，其中缩小比为减少量与初始量的比值。
表 7 各组粉末孔隙率Table 7. Powder porosity of each group试样编号 孔隙率 / % 初始 最终 减少 缩小比 N组 55.00 30.00 25.00 0.45 A组 55.00 24.25 30.75 0.56 B组 55.00 21.61 33.39 0.61 C组 55.00 19.98 35.02 0.64 D组 55.00 18.08 36.92 0.67 E组 55.00 18.63 36.37 0.66 从表7可以看出，最终孔隙率最低的是D组混合粉末，达到了18.08%，并且其孔隙率减少最多，缩小比最大达到了0.67，为最终相对密度最高的一组粉末试样，其压缩性最佳。单一粒径粉末N组的最终孔隙率为30.00%，远高于其他混合粉末，最终缩小比最低，为0.45，说明按照单一粒径进行模拟的结果跟实际情况存在较大区别。混合粉末的压缩性整体趋势是随着附加细粉末（粒径小于45 μm）的质量占比增加，其压缩性增加，体现在压缩性从差到好排序为A组、B组、C组、E组、D组。D组反而比E组的压缩性能更好，是由于过多的细微粉末导致粉末体流动性变差，不利于压缩过程中的粉末颗粒位移。
2.2 压坯应力
2.2.1 侧压力与轴向应力
金属粉末在受到轴向应力时，压坯会延着径向膨胀，模壁会给压坯一个大小相等、方向相反的反作用力，叫侧压力，属于径向应力。在压制过程中，压坯与模壁发生相对运动，必然受到摩擦力作用，径向应力的大小影响着摩擦力大小，从而影响压制过程，因而有必要研究径向应力的大小，这对研究压制过程和压坯致密性都有意义。
径向应力与轴向应力关系如图7所示，把压坯简化为一个立方体，当压坯受到轴向应力时，会使其在x方向上产生膨胀，如式（5）所示，膨胀值（
$ {\Delta _{{L_{x1}}}} $ ）与材料泊松比（v）和轴向应力（P轴向）成正比，与弹性模量（E）成反比；如式（6）所示，y方向的径向应力（P径向）也使压坯沿x方向膨胀（$ {\Delta _{{L_{{\text{x2}}}}}} $ ）；如式（7）所示，沿x方向的径向应力（P径向）使压坯沿x方向压缩（$ {\Delta _{{L_{x3}}}} $ ）。$$ {\Delta }_{{L}_{x1}}=v\frac{{P}_{轴向}}{E} $$ (5) $$ {\Delta }_{{L}_{x2}}=v\frac{{P}_{径向}}{E} $$ (6) $$ {\Delta }_{{L}_{x3}}=\frac{{P}_{径向}}{E} $$ (7) 实际上压坯在压膜内不能径向膨胀，因此在x方向膨胀值总和等于其压缩值，如式（8）所示。由式（5）～式（8）可得到径向应力和轴向应力的关系，如式（9）和式（10）所示，两者之比即粉末体侧压系数（ξ）。
$$ {\Delta _{{L_{x1}}}}{\text{ + }}{\Delta _{{L_{x2}}}} = {\Delta _{{L_{x3}}}} $$ (8) $$ \xi =\frac{{P}_{径向}}{{P}_{轴向}}=\frac{\nu }{1-\nu } $$ (9) $$ \nu = \frac{\xi }{{1 + \xi }} $$ (10) 2.2.2 轴向应力对比
测得压制方向应力‒应变曲线如图8所示。由图可知，单一粒径组（N组）粉末试样在压制过程中的应力远高于其他混合粉末试样，其余各组混合粉末试样随着附加细粉末（粒径小于45 μm）质量占比的上升，粉末在发生相同应变的情况下受到的轴向应力是下降的，这说明混合粒径有利于减少压制过程中的能耗。轴向应力按从大到小排序为A组、B组、C组、D组、E组，其中D组和E组非常接近，但整体还是D组应变最大。
2.2.3 侧压系数与泊松比对比
根据式（8）利用6组粉末的径向应力数据与轴向应力数据可求得侧压系数，进一步通过式（9）求得泊松比，结果如图9和图10所示。由图9可知，侧压系数随着应变增加而增加，在压制初期，测压系数快速提升，到稳定阶段后，各组粉末侧压系数都能达到0.6左右，随着压制继续进行又逐渐减低。粉末粒径分布越单一，侧压系数越早下降，随着附加细粉末（粒径小于45 μm）在粉末中质量占比的增加，侧压系数缓慢下降，当进行到压制后期，各组孔隙率变化都不大，各组侧压系数值均接近与0.5，各组侧压系数变化趋势基本相同，说明侧压系数的大小受粉末粒径分布的影响较小。
压制过程中的粉末泊松比变化曲线如图10所示。由图可知，各组曲线都比较接近，平稳期泊松比保持在0.36，说明粉末体泊松比变化主要是受粉末材料参数的影响，粒径分布对于泊松比影响较小。值得注意的是，压制后期的粒径分布单一粉末在径向位移上缺乏足够的驱动力与位移空间，径向压制量较轴向压制先弱化，导致泊松比下降。
2.3 力链分布
各组粉末试样在外载荷作用下，颗粒间发生接触，形成了诸多强度不同的力链，分别用强弱力链表示，定义高于平均接触力所构成的力链是强力链，低于则称为弱力链。各组粉末试样压制后力链分布如图11所示，其中深色表示弱力链，亮色表示强力链。由图可知，单一粒径组（N组）粉末的力链数量远小于混合粉末组，随着粉末颗粒数量增加，颗粒之间的接触数量增加，力链数量也随之增加。E组粉末试样力链数量最多，达到了815648个。各组力链数量以及力链强度如表8所示。由表可知，单一粒径粉末的力链数量虽然少，但强力链的数目较多，混合粉末中有数目巨大的弱力链均匀分布在颗粒体系内，参与强力链的衔接和断后重构等辅助作用。起主要支撑作用的强力链数目充足，同时搭配一定数目的弱力链，才能保证压坯相对密度足够高。
![]()
图 11 各组粉末试样力链分布：（a）N组；（b）A组；（c）B组；（d）C组；（e）D组；（f）E组Figure 11. Force chain distribution of the powders for each group: (a) group N; (b) group A; (c) group B; (d) group C; (e) group D; (f) group E表 8 各组粉末试样力链数目与强度Table 8. Number and strength of force chain distribution of the powders for each group试样组 力链数目 强度范围 / MPa N 67453 2.0561×10−2～1.5216×103 A 288769 1.8707×10−4～1.3296×103 B 416107 3.0393×10−4～1.3683×103 C 549172 5.7130×10−4～1.4514×103 D 704764 5.1223×10−4～1.4251×103 E 815648 5.0515×10−4～1.3222×103 对比压坯上表面与侧壁的力链数目，可以明显发现上表面有大量树干状强力链，侧壁上的力链基本处于强度较低的强力链末端或者是成网状的弱力链，这会导致压坯表面粗糙度的不同。
3. 结论
（1）随着压制力的增大，粉末压坯相对密度逐渐提高，但在400 MPa时出现拐点。当粉末粒径分布质量比为（基础粒径：中粒径：细粒径）为60:15:25时，粉末压坯孔隙率达到了18.08%，且孔隙率减少最多，缩小比最大达到了0.67，最终相对密度最高，压缩性最大。证明了并不是细粒径粉末质量占比越高，粉末压坯相对密度越高，当细粉末质量占比达到一定比例时，会导致粉末流动性变差，影响相对密度提高。
（2）在基础粉末粒径为100目，各组混合粒径粉末试样轴向应力达到500 MPa的情况下，通过调整附加细粉末（粒径小于45 μm）的质量占比，压制方向上能产生更大应变。当基础粉末质量分数为60%，细粒径粉末质量分数为25%时，应变最大，有助于减少压制能量的消耗。侧压系数与泊松比受粉末粒径分布影响较小。在压制后期，当压坯已经获得较高相对密度时，侧压系数与泊松比会发生下降。
（3）在基础粒径粉末质量占比不变的情况下，随着中径粉末与细经粉末质量分数的增加，粉末颗粒数量会增加，力链数量也随之增加，细粉质量分数最高的试样力链数量最多，达到了814927个。在混合粒径粉末中有大量弱力链均匀分布在颗粒体系内，参与强力链的衔接和断后重构等辅助作用。在强力链数目充足的前提下，配合大量弱力链能获得更高的压坯相对密度。
-
![]()
图 1 接触模型中接触单元关系：（a）球-球类型；（b）球-墙类型
Figure 1. Relationship of the contact unit in the contact model: (a) ball-ball type; (b) ball-facet type
![]()
图 2 动态测量球分布图：（a）初始测量球俯视图；（b）初始测量球正视图；（c）压制后测量球俯视图；（d）压制后测量球正视图
Figure 2. Distribution of the dynamic measuring ball: (a) top view of the initial measuring ball; (b) front view of the initial measuring ball; (c) top view of the measuring ball after pressing; (d) front view of the measuring ball after pressing
![]()
图 5 各组粉末压制模型：（a）N组压制前；（b）N组压制后；（c）A组压制前；（d）A组压制后；（e）B组压制前；（f）B组压制后；（g）C组压制前；（h）C组压制后；（i）D组压制前；（j）D组压制后；（k）E组压制前；（l）E组压制后
Figure 5. Each group of powder pressing model: (a) before N group suppression; (b) after N group suppression; (c) before A group suppression; (d) after A group suppression; (e) before B group suppression; (f) after B group suppression; (g) before C group suppression; (h) after C group suppression; (i) before D group suppression; (j) after D group suppression; (k) before E group suppression; (l) after E group suppression
![]()
图 11 各组粉末试样力链分布：（a）N组；（b）A组；（c）B组；（d）C组；（e）D组；（f）E组
Figure 11. Force chain distribution of the powders for each group: (a) group N; (b) group A; (c) group B; (d) group C; (e) group D; (f) group E
表 1 工业标准筛
Table 1 Industrial standard screen
目数 筛网孔径 / μm 粒径 / μm +80 +180 ≥180 −80～+100 −180～+150 150～180 −100～+150 −150～+106 106～150 −150～+200 −106～+75 75～106 −200～+250 −75～+63 63～75 −250～+325 −63～+45 45～63 −325 −45 ≤45 
下载: 导出CSV
表 2 各标准粒径分布对比
Table 2 Comparison of the standard particle size distributions
牌号 压缩性 / (g·cm‒3) 粒径分布（质量分数） / % ＞250 μm ＞180 μm ＞150 μm ＞75 μm ＞63 μm ＞45 μm ＜45 μm FHY100·270 ≥6.70 0 0 3 — — — 10～30 300NH ≥7.05 1 15 15 60 5～20 8～23 10～30 HAP100.30H ≥7.10 0 1 10 60 — 20～40 30～45 ABC100.30 ≥7.28# 0 * 10 * * * 20 FSW 100·30H 7.15# 0 1 10 — — — 15～30 注：“—”表示不作具体要求；“*”表示参数保密；压缩性的单位压制压力除特别标注“#”为600 MPa外，其余均为500 MPa。
下载: 导出CSV
表 3 200目纯铁粉试样粉末粒径分布（质量分数）
Table 3 Particle size distribution of the pure iron powder samples in 200 mesh
% ＞150 μm ＞106 μm ＞75 μm ＞45 μm ＜45 μm 0 5 15 40 40
下载: 导出CSV
表 4 铁粉相关参数
Table 4 Related parameters of Fe powder
材料 密度 /
(kg·m−3)弹性模量 /
GPa法向切向
刚度比摩擦系数 铁粉 7800 190 0.785 0.3
下载: 导出CSV
表 5 各组试样粉末粒径分布
Table 5 Particle size distribution of the samples
试样编号 粒径分布（质量分数） / % ＞150 μm 106～150 μm 45～106 μm ＜45 μm N 0 100 0 0 A 0 60 30 10 B 0 60 25 15 C 0 60 20 20 D 0 60 15 25 E 0 60 10 30
下载: 导出CSV
表 6 各粒径颗粒数目
Table 6 Number of the particles in each size
试样编号 粒径分布（个数） 150 μm 106 μm 45 μm 总数 N 18076 0 0 18076 A 10822 15626 68923 95371 B 10822 13073 103432 127327 C 10822 10497 137742 159061 D 10822 7897 172301 191020 E 10822 5275 206605 222702
下载: 导出CSV
表 7 各组粉末孔隙率
Table 7 Powder porosity of each group
试样编号 孔隙率 / % 初始 最终 减少 缩小比 N组 55.00 30.00 25.00 0.45 A组 55.00 24.25 30.75 0.56 B组 55.00 21.61 33.39 0.61 C组 55.00 19.98 35.02 0.64 D组 55.00 18.08 36.92 0.67 E组 55.00 18.63 36.37 0.66
下载: 导出CSV
表 8 各组粉末试样力链数目与强度
Table 8 Number and strength of force chain distribution of the powders for each group
试样组 力链数目 强度范围 / MPa N 67453 2.0561×10−2～1.5216×103 A 288769 1.8707×10−4～1.3296×103 B 416107 3.0393×10−4～1.3683×103 C 549172 5.7130×10−4～1.4514×103 D 704764 5.1223×10−4～1.4251×103 E 815648 5.0515×10−4～1.3222×103
下载: 导出CSV
-
[1] 黄培云. 粉末冶金原理. 1版. 北京: 冶金工业出版社, 1982 Huang P Y. Theory of Power Metallurgy. 1st Ed. Beijing: Metallurgical Industry Press, 1982
[2] Amherd Hidalgo A, Frykholm R, Ebel T, et al. Powder metallurgy strategies to improve properties and processing of titanium alloys: A review. Adv Eng Mater, 2017, 19(6): 1600743 DOI: 10.1002/adem.201600743
[3] 郎利辉, 王刚, 黄西娜, 等. 粉末粒度对热等静压法制备2A12铝合金组织与性能的影响. 粉末冶金材料科学与工程, 2016, 21(1): 85 DOI: 10.3969/j.issn.1673-0224.2016.01.012 Lang L H, Wang G, Huang X N, et al. Effect of powder size on microstructure and properties of 2A12 aluminium alloy prepared by hot isostatic pressing. Mater Sci Eng Powder Metall, 2016, 21(1): 85 DOI: 10.3969/j.issn.1673-0224.2016.01.012
[4] 林立, 刘军, 周纯, 等. 金属粉末冲击压制过程中冲模质量及颗粒模型的优化分析. 粉末冶金技术, 2018, 36(3): 182 Lin L, Liu J, Zhou C, et al. Optimization analysis of die mass and particle model in metal powder impact compaction. Powder Metall Technol, 2018, 36(3): 182
[5] 丁义超, 尹红, 姜自莲. 铁基复合材料的制备技术与研究进展. 热加工工艺, 2013, 42(24): 22 Ding Y C, Yin H, Jiang Z L. Preparation technology and research progress of iron-based composite. Hot Working Technol, 2013, 42(24): 22
[6] 耿学文, 赵洪波, 樊振军. 铁基复合材料的研究进展综述. 中国科技信息, 2009(6): 34 Geng X W, Zhao H B, Fan Z J. Study of ferrous matrix composites. China Sci Technol Inf, 2009(6): 34
[7] 张超, 刘军, 罗晓龙, 等. 基于离散元法的金属粉末压制加载速度对压力分布影响. 粉末冶金技术, 2019, 37(2): 98 Zhang C, Liu J, Luo X L, et al. Effect of loading speed on pressure distribution in metal powder pressing based on discrete element method. Powder Metall Technol, 2019, 37(2): 98
[8] Skrinjar O, Larsson P L. On discrete element modelling of compaction of powders with size ratio. Comput Mater Sci, 2004, 31(1-2): 131 DOI: 10.1016/j.commatsci.2004.02.005
[9] 叶先勇, 刘京雷, 徐宏, 等. 粉末粒径和压制压力对316L不锈钢多孔材料结构特性的影响. 粉末冶金材料科学与工程, 2013, 18(3): 409 DOI: 10.3969/j.issn.1673-0224.2013.03.017 Ye X Y, Liu J L, Xu H, et al. Effects of powder size and molding pressure on structural characterization of 316L stainless steel porous material. Mater Sci Eng Powder Metall, 2013, 18(3): 409 DOI: 10.3969/j.issn.1673-0224.2013.03.017
[10] 朱鹏程. 粉末粒度与制备工艺对烧结钕铁硼性能的影响[学位论文]. 南京: 南京理工大学, 2012 Zhu P C. Effects of the Powder Size and Preparation Process on the Performance of Sintered NdFeB [Dissertation]. Nanjing: Nanjing University of Science and Technology, 2012
[11] 闫志巧, 陈峰, 蔡一湘. 不同粒径Ti粉的高速压制行为和烧结性能. 金属学报, 2012, 48(3): 379 DOI: 10.3724/SP.J.1037.2011.00612 Yan Z Q, Chen F, Cai Y X. High velocity compaction behavior and sintered properties of Ti powders with different particle sizes. Acta Metall Sin, 2012, 48(3): 379 DOI: 10.3724/SP.J.1037.2011.00612
[12] Yan Z Q, Chen F, Cai Y, et al. Influence of particle size on property of Ti‒6Al‒4V alloy prepared by high-velocity compaction. Trans Nonferrous Met Soc China, 2013, 23(2): 361 DOI: 10.1016/S1003-6326(13)62470-X
[13] 杨志超. 粉末粒度分布对烧结钕铁硼微观结构形成的影响. 山西冶金, 2014, 37(5): 11 DOI: 10.3969/j.issn.1672-1152.2014.05.004 Yang Z C. Effect of particle size distribution on the formation of the microstructure of sintered NdFeB. Shanxi Metall, 2014, 37(5): 11 DOI: 10.3969/j.issn.1672-1152.2014.05.004
[14] 潘诗琰, 代文杰, 周子豪, 等. 液相烧结过程中粉末粒径分布演化模拟研究. 粉末冶金技术, 2018, 36(6): 409 Pan S Y, Dai W J, Zhou Z H, et al. Simulation study on size distribution evolution of powder particles in liquid phase sintering. Powder Metall Technol, 2018, 36(6): 409
[15] Cundall P A, Strack O D L. Discussion: A discrete numerical model for granular assemblies. Géotechnique, 1980, 30(3): 331
[16] 中华人民共和国国家质量监督检验检疫总局. GB/T 1480-2012金属粉末干筛分法测定粒度. 北京: 中国标准出版社, 2013 General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China. GB/T 1480-2012 Metal Powder Dry Screening Method. Beijing: China Standard Press, 2013
[17] Japanese Industrial Standards Committee. JIS Z 2510-2004 Metallic Powders—Determination of Particle Size by Dry Sieving. Tokyo: Japan Marine Standard Association, 2004
[18] 中华人民共和国工业和信息化部. YB/T 5308-2011粉末冶金用还原铁粉. 北京: 冶金工业出版社, 2012 Ministry of Industry and Information Technology of the People's Republic of China. YB/T 5308-2011 Reduction Iron Powder for Powder Metallurgy. Beijing: Metallurgical Industry Press, 2012
[19] 全国生铁及铁合金标准化技术委员会. 20153628-T-602粉末冶金用水雾化纯铁粉、合金铁粉国家标准编制说明. 北京: 中国标准出版, 2015 National Pig Iron and Ferroalloy Standardization Technical Committee. 20153628-T-602 Water Atomized Pure Iron Powder for Powder Metallurgy, Specification for the Preparation of National Standards for Alloy Iron Powder. Beijing: China Standard Press, 2015
[20] 国家市场监督管理总局. GB/T 19734-2018粉末冶金用水雾化纯铁粉、合金钢粉. 北京: 中国标准出版社, 2018 State Administration of Market Regulation. GB/T 19734-2018 Water Atomized Pure Iron Powder, Alloy Steel Powder for Powder Metallurgy. Beijing: China Standard Press, 2018
[21] 胡仙平. 高速加载下金属颗粒接触过程的影响因素分析[学位论文]. 宁波: 宁波大学, 2016 Hu X P. Analysis on the Influential Factors of Metal Particle Contact Process under High Speed Loading [Dissertation]. Ningbo: Ningbo University, 2016
-
期刊类型引用(7)
1. 伍威，邓汲坤，任伯勇，刘君武，仲洪海，蒋阳. 注射成形低钨-钨连接度的90W-6Ni-4Mn的高比重合金及其组织与性能研究. 粉末冶金工业. 2024(03): 21-27+33 . 
百度学术
2. 陈新宇，李奋强，蒋继帅. 数值模拟在粉末压制细观分析中的应用与发展. 粉末冶金技术. 2024(04): 418-426 . 本站查看
3. 张炜，萧伟健，袁传牛，张宁，陈荣昕. 基于三维离散元模型粉末压制中力链对阻塞行为的影响机制. 粉末冶金技术. 2024(04): 403-410+417 . 本站查看
4. 王梦奇，刘越，肖春林，刘春明. SiCp粒径级配对55%SiCp/6061Al复合材料组织和性能的影响. 东北大学学报(自然科学版). 2024(06): 802-807 . 百度学术
5. 张炜，袁传牛，萧伟健，陈荣昕，张宁，刘焜. 粉末压制中粉末形状对粉体力学行为的影响. 锻压技术. 2023(07): 138-148 . 百度学术
6. 黄丽容，王云，谢俊杰，黄必华. 超声振动辅助压制粗晶WC-10Co硬质合金. 粉末冶金技术. 2023(06): 586-592 . 本站查看
7. 张炜，萧伟健，袁传牛，张宁，刘焜. 离散元法铁粉末压制中粒径分布对力链演化机制的影响. 力学学报. 2022(09): 2489-2500 . 百度学术
其他类型引用(5)
计量
- 文章访问数: 793
- HTML全文浏览量: 190
- PDF下载量: 117
- 被引次数: 12
