Effect of loading speed on pressure distribution in metal powder pressing based on discrete element method
-
摘要: 基于离散元原理,采用三维离散元软件PFC3D模拟冲头加载过程,研究以不同冲头速度撞击金属粉末颗粒所产生的透射波以及对侧壁的压力分布。结果表明,对于离散性介质,冲击加载下的应力波传播不单受到颗粒材质波阻抗的影响,还受到颗粒运动和力链形成的很大影响。通过对比不同冲击加载速度下阴模底部的受力,发现随着冲击速度的提高,透射波的峰值成线性增长;通过均分阴模侧壁的方式得到侧壁不同位置所受的侧压力峰值,发现侧壁所受压力呈波谷式分布,从而分析得到颗粒间力链的传播和分布规律。Abstract: The punch loading process was simulated by three-dimensional discrete element software PFC3D based on discrete element method (DEM), the effects of punch loading speeds on the transmission wave of metal powders generated by the impact and the pressure distribution on the sidewall were investigated. In the results, for the discrete media, the stress wave propagation under impact loading is determined not only by the wave impedance of particle material, but also by the particle motion and force chain formation. By comparing the force on the bottom of female mold in the different loading speeds, it is found that the peak of transmission wave increases linearly with the increase of loading speed. The lateral pressure peak on different positions is obtained by means of the side wall of female mold. It is found that the pressure on the sidewalls is present as the valley-like distribution, and the propagation and distribution of the inter-particle force chains are analyzed.
-
Keywords:
- discrete element /
- impact loading /
- transmission wave /
- lateral pressure /
- stress wave
-
钨是一种难熔过渡族稀有金属,熔点高达3410 ℃,密度高达19.25 g·cm−3,具有高硬度、高强度、高弹性模量、良好的导热性和低蒸汽压等优异的物理性能,在一些特殊装备领域展现出不可替代的作用。钨板、钨坩埚是重要的深加工、高附值产品,在热屏、靶材、蓝宝石晶体生长、稀土冶炼、真空蒸镀等行业获得广泛应用[1–5]。随着产品尺寸不断增大,提高钨板、钨坩埚的整体组织均匀性和密度均匀性成为重要的研究内容[6–7]。
等静压–烧结法是当今钨产品最主要的制备技术[8]。原料粉体的颗粒特性对烧结过程中的致密化和晶粒生长具有显著影响。研究表明,粉末细颗粒具有易团聚、形状不规则的特点,在压制过程中会造成大孔隙和不均匀堆积等现象,导致烧结过程中的不均匀致密化,甚至引起整体相对密度降低和晶粒异常长大[9–12]。因此,对原料粉体进行改性是解决上述问题的关键方法之一。
气流磨分级技术是一种高效的物理改性粉体的方法,它利用高速涡流气体带动粉体颗粒运动,使粉末颗粒之间相互冲击、碰撞、摩擦,达到研磨并使团聚颗粒分散的目的,并通过涡轮分级机将粉末颗粒分级[13–14]。本文以商业微米级纯钨粉为原料,采用气流磨分级技术对粉末进行改性处理,获得了粒度分布窄、分散性好和球形度好的钨粉,并对改性钨粉进行粒度搭配,研究了不同粒度分布钨粉烧结体的显微组织和力学性能。
1. 实验材料及方法
实验所用原料粉末是厦门金鹭公司的商业钨粉(费氏粒度为5 μm,纯度>99.9%),采用对喷式气流磨装置,其工作原理如图1所示[13]。以氮气为研磨介质,设定研磨腔压力为0.7 MPa,在分选轮频率10~50 Hz范围内对原料钨粉进行改性处理,并对不同分选轮频率处理得到的钨粉进行粒度搭配,以获得不同粒度分布的粉末,探究气流磨前后的粉体状态以及粉末粒径搭配对粉体烧结性能的影响。实验共设置四种批次粉末进行对照实验:原料钨粉(原粉W0)、20 Hz处理的气流磨粉末(气流磨W1)、W1与粗颗粒混合搭配粉末(气流磨W2)以及W1与细颗粒混合搭配粉末(气流磨W3),实验所述粗颗粒由气流磨10 Hz条件下处理而得,细颗粒由气流磨50 Hz条件下处理而得。
分别对四种批次钨粉进行冷等静压+中频感应烧结处理。冷等静压设置压制力为200 MPa,保压时间为5 min。将压坯置于中频感应炉中进行烧结,全程通入氢气气氛,最高烧结温度为2000 ℃,保温2 h。四种批次钨粉的压制烧结实验全程在相同条件下完成。
采用BT-9300ST激光粒度仪测试粉末的粒度分布及比表面积;利用Rigaku-D/MAX-RB旋转阳极X射线衍射仪(X-ray diffraction,XRD)测定粉末的物相;使用SUPRA 55场发射扫描电子显微镜(scanning electron microscope,SEM)观察粉末形貌以及烧结态样品的显微组织和孔隙分布;通过阿基米德排水法测定烧结态样品的密度,为避免残余孔隙影响,测量前将样品表面涂抹凡士林;利用WDW-50万能试验机对烧结态样品进行室温三点抗弯测试,测试样品尺寸为3 mm×3 mm×21 mm,跨距设置为10 mm,测试前样品表面被机械抛光至光亮。
2. 结果与讨论
2.1 气流磨处理前后粉末性能变化
图2所示为气流磨处理前后粉末的扫描电子显微形貌。由图可知,气流磨处理前的原始粉体颗粒呈现出不规则的多面体形状,且存在显著的团聚现象,能观察到较大的团聚体,如图2(a)所示;经气流磨处理后,粉末原有的团聚体打开,呈现良好的分散状态,单个颗粒呈现近球形,如图2(b)~图2(d)所示。这是因为原始粉末在高速氮气的带动下,粉末颗粒之间发生剧烈碰撞,碰撞时产生的冲击力使团聚体打开并使单个颗粒分散,分散后的颗粒进一步碰撞、摩擦,使得原始颗粒的棱角被磨平,从而处理后的颗粒呈现出近球形,然后在涡轮分级过程中把大颗粉末或碰撞碎屑筛分出去,使粉末颗粒粒径趋于同一尺寸[14]。气流磨粉末W1为未进行粒度搭配的单一粒径分布粉体,气流磨粉末W2和W3经过粉末粒度搭配,从图2(c)和图2(d)中可以观察到人为掺入的较大和较小颗粒。
图3直观地显示了四种批次粉末的粒度分布图,具体粒度分布参数如表1所示。常采用D10、D50、D90三个参数表示粉末粒径分布特征,其物理意义为累计粒径分布分别达到10%、50%、90%时的粒径值。原粉W0的粒径分布最宽,D10=5.95 μm,D50=12.57 μm,D90=23.03 μm;气流磨W1的粒径分布最窄,D10=5.30 μm,D50=8.02 μm,D90=11.79 μm。比表面积在一定程度上反映了粉末的分散程度,原粉的比表面积为187.2 m2·kg−1,小于气流磨处理的钨粉(231.6~348.7 m2·kg−1),这也印证了气流磨处理后的粉末具有更好的分散性。
表 1 不同批次粉末的粒度分布参数Table 1. Particle size distribution parameters of the different powders批次 D10 / μm D50 / μm D90 / μm 比表面积 / (m2·kg−1) 原粉W0 5.95 12.57 23.03 187.2 气流磨W1 5.30 8.02 11.79 259.5 气流磨W2 5.49 8.91 15.72 231.6 气流磨W3 2.99 7.14 11.53 348.7 图4(a)为气流磨处理前后粉末的X射线衍射分析。由图可知,气流磨处理前后粉末的衍射谱均只观察到钨的衍射峰,没有其他杂质物相产生,说明气流磨处理能达到改性粉末效果的同时也能保证粉末的纯净度,这是因为整个气流磨过程在高纯氮气的保护下进行,未发生氧化反应或引入杂质相。Williamson-Hall方法被广泛用于计算金属材料的晶粒尺寸、微应变和位错密度[15–16],Williamson-Hall公式如式(1)所示。
$${\delta _{hkl}}\cos {\theta _{hkl}} = \frac{{K\lambda }}{D} + 4\varepsilon \sin {\theta _{hkl}}$$ (1) 式中:δhkl为各个衍射峰的半高宽(full width at half maximum,FWHM),θhkl为布拉格衍射角,K为常数,λ为衍射波长,D表示晶粒尺寸,ε表示微观应变值。绘制δhklcosθhkl和4sinθhkl函数并进行线性拟合,结果如图4(b)所示,其中斜率大小可表示样品的晶格畸变程度。由图可知,气流磨后粉末的衍射峰半高宽要略大于原粉的半高宽。经定量处理,气流磨处理后粉末的晶格畸变要显著大于原粉,这反映了气流磨过程中粉末颗粒间的机械碰撞可能会使粉末颗粒的晶格畸变增大,使粉末的缺陷密度增大,从而影响其烧结性能。
2.2 烧结态样品的显微组织
对四种不同批次的粉末在完全相同的压制和烧结条件进行固结,从而对比观察烧结态样品的显微组织。根据实验结果讨论两个问题:一是气流磨工艺改性粉末对烧结组织的影响;二是颗粒尺寸分布对烧结组织的影响。
由于所用钨粉粒径较大,即使在2000 ℃温度下烧结仍难以实现完全致密化,存在大量残余孔洞。为了清晰观察烧结孔隙的形貌和分布状态,利用扫描电子显微镜观察烧结态样品的抛光面。图5显示了不同批次粉末烧结态的孔隙结构。首先定义两种孔隙结构,即一次粗孔和二次细孔。一次粗孔尺寸较大,孔径大于10 μm,形状呈现不规则放射状;二次细孔呈现规则的球形或椭球形,尺寸远小于一次粗孔,约为5 μm以下。由于一次粗孔对材料性能的恶化影响要远大于二次细孔,因此烧结过程需要避免一次粗洞的形成。由图5可知,W0、W2和W3样品存在一些一次粗孔,而具有更窄粒度分布的W1样品(未进行粒度搭配的气流磨处理粉末)基本消除了不规则形状的一次粗孔,孔隙结构以均匀分布的二次细孔为主。对孔隙的尺寸和形状进行定量分析,图6显示了不同批次粉末烧结后的孔隙孔径和长径比散点分布图,散点分布越集中说明孔隙分布越均匀规则。由图6可见,气流磨处理的三个样品的散点分布图比原粉W0样品均更加集中,孔径和长径比分布要更窄,统计计算得到平均孔径和平均长径比也略小,这说明了经过气流磨处理的粉末烧结后的孔隙结构更加均匀且规则。
不同的孔隙结构是由生坯中不同的颗粒堆积状态造成的[17]。Kingery和Francois[18]首次提出孔洞配位数的概念,即孔隙周围晶粒的数量。配位数较大的孔洞趋于长大,而配位数较小的孔洞趋于收缩。颗粒堆积紧密的区域在烧结致密化过程中产生配位数较小的孔洞,形成规则的闭孔;颗粒堆积松散的区域产生配位数较大的孔洞,即使通过后续的高温烧结也很难消除。气流磨粉末W1粒度分布更窄,分散性更好,从而实现良好的堆积均匀性,消除了团聚体在烧结过程中形成的大孔隙,从而提升了结构均匀性,提高了烧结后期整体相对密度。以上结果表明,消除一次粗孔、提高烧结后期组织结构均匀性的主要因素是颗粒尺寸分布而不是颗粒形状。
2.3 烧结态样品的性能
对上述四种批次粉末烧结态样品进行了室温抗弯实验。烧结态样品的各项性能参数如表2所示,相比于原粉W0样品,气流磨处理后的样品具有更高的抗弯强度。图7为四种批次粉末烧结态样品抗弯实验的载荷–位移曲线,四种样品均展现出脆性断裂的特征,在同等样品尺寸下,气流磨样品比原粉样品展现出更高的断裂载荷。力学性能的提升大致可归因于两点:一是气流磨处理的样品具有更高的相对密度;二是气流磨处理的样品具有更加均匀的孔隙结构。值得注意的是,粉末粒度分布最窄的样品W1具有最高的相对密度(92.8%)和最高的抗弯强度(292.4 MPa),这可以归因于其优良的粉末特性。如前面所述,烧结样品中的一次粗孔对性能可能产生更大的恶化影响,在这里得到印证。气流磨W1样品由于其较好的粉末分散性和窄粒径分布,避免了团聚体在烧结过程中形成一次粗孔,优化了力学性能。但是四种烧结样品的总体抗弯强度并不高,这是因为其整体相对密度不高所导致的。
表 2 不同批次粉末的烧结后性能参数Table 2. Properties parameters of the different sintered samples批次 相对密度 / % 平均晶粒尺寸 / μm 抗弯强度 / MPa 原粉W0 90.7 10.5 238.5±21.3 气流磨W1 92.8 10.1 292.4±10.5 气流磨W2 92.6 10.3 263.1±14.6 气流磨W3 92.3 10.4 248.3±15.1 3. 结论
(1)通过气流磨处理得到了近球形、分散性好、粒度分布窄的微米级钨粉,改性后的钨粉纯净度高,无杂质相生成。经气流磨处理后,商业钨粉的D10,D50,D90分别由处理前的5.95 μm、12.57 μm、23.03 μm降低为5.30 μm、8.02 μm、11.79 μm,比表面积由187.2 m2·kg−1增大到259.5 m2·kg−1。
(2)对比研究了气流磨处理和不同粒度搭配钨粉的烧结孔隙特征。经烧结得到两种孔隙结构,一种是孔径大于10 μm的不规则粗孔,另一种是孔径小于5 μm的近球形细孔。气流磨处理的钨粉烧结后的孔隙结构整体更加均匀且规则,且粒度分布较窄的处理粉末可完全消除第一类孔隙,说明消除烧结体不规则粗孔的关键因素是窄化粉末粒度分布。
(3)对比研究了气流磨处理和不同粒度搭配钨粉的力学性能。气流磨处理后的钨粉比处理前的钨粉具有更好的烧结性能,而不同粒度搭配的气流磨粉末中颗粒粒度分布最窄的粉末烧结性能最佳,烧结相对密度由90.7%提高至92.8%,相应地,抗弯强度由238.5 MPa提升至292.4 MPa。烧结形成的不规则粗孔是导致力学性能恶化的主要原因,气流磨分级处理并窄化颗粒粒度分布可消除烧结钨中的粗孔并提升组织结构均匀性,从而优化其性能。
-
表 1 金属粉末冲击加载模拟主要参数
Table 1 Main simulation parameters of metal powder impact loading
密度/(kg·m-3) 颗粒直径/mm 颗粒间摩擦系数 模壁摩擦系数 刚度系数 冲击速度/(m·s-1) 2700 0.15 0.3 0.1 108 10 表 2 冲击速度与对应的运行时间
Table 2 Impact speed and corresponding running time
冲模速度/(m·s-1) 5 7 10 14 20 运行总时间/μs 70 50 35 25 17.5 表 3 不同冲击加载速度和压制时段的透射波振动次数
Table 3 Vibration frequency of transmission wave in different impact loading speeds at different suppression periods
冲击加载速度/(m·s-1) 透射波振动次数 前期 前中期 中期 中后期 后期 5 55 50 34 30 23 7 43 35 17 21 11 10 30 20 14 15 7 14 19 12 12 6 2 20 12 9 2 4 1 表 4 不同颗粒填充高度下侧压力峰值与冲击加载速度的关系
Table 4 Relationship of impact loading speeds and peak values of lateral pressure on die wall in different particle-filled heights
冲击速度/(m·s-1) 最大压力/(102N) 一层 二层 三层 四层 五层 5 0.52 1.03 1.20 1.52 0.92 7 0.58 1.31 1.20 1.59 1.40 10 0.84 2.12 1.40 2.10 0.70 14 0.86 2.44 2.04 2.44 1.50 20 1.50 3.55 3.60 3.65 2.75 -
[1] 杨霞. 粉末冶金高速压制致密化机制的研究进展. 粉末冶金工业, 2016, 26(5): 57 https://www.cnki.com.cn/Article/CJFDTOTAL-FMYG201605019.htm Yang X. Research progress on densification mechanism of powder metallurgy highvelocity compaction. Powder Metall Ind, 2016, 26(5): 57 https://www.cnki.com.cn/Article/CJFDTOTAL-FMYG201605019.htm
[2] 王建忠, 曲选辉, 尹海清, 等. 铁粉的高速压制成形. 材料研究学报, 2008, 22(6): 589 DOI: 10.3321/j.issn:1005-3093.2008.06.006 Wang J Z, Qu X H, Yin H Q, et al. High velocity compaction of ferrous powder. Chin J Mater Res, 2008, 22(6): 589 DOI: 10.3321/j.issn:1005-3093.2008.06.006
[3] 何杰, 肖志瑜, 关航健, 等. 纯钛粉高速压制行为及其烧结性能研究. 粉末冶金技术, 2016, 34(3): 178 DOI: 10.3969/j.issn.1001-3784.2016.03.004 He J, Xiao Z Y, Guan H J, et al. High velocity compaction behavior and sintered properties of pure Ti powder. Powder Metall Technol, 2016, 34(3): 178 DOI: 10.3969/j.issn.1001-3784.2016.03.004
[4] 尤萌萌, 潘诗琰, 申小平, 等. 粉末压制过程数值模拟的研究现状及展望. 粉末冶金工业, 2017, 27(4): 49 https://www.cnki.com.cn/Article/CJFDTOTAL-FMYG201704015.htm Yu M M, Pan S Y, Shen X P, et al. Current progress and prospect of numerical simulation in powder compaction. Powder Metall Ind, 2017, 27(4): 49 https://www.cnki.com.cn/Article/CJFDTOTAL-FMYG201704015.htm
[5] 胡建华, 尚会森, 程呈, 等. 金属粉末压制成形理论与工艺进展. 热加工工艺, 2012, 41(20): 45 https://www.cnki.com.cn/Article/CJFDTOTAL-SJGY201220012.htm Hu J H, Shang H S, Cheng C, et al. Research progress of metal powder compression theory and technology. Hot Working Technol, 2012, 41(20): 45 https://www.cnki.com.cn/Article/CJFDTOTAL-SJGY201220012.htm
[6] 孙婧, 何佩珊, 齐梦菊. 关于颗粒流软件PFC的离散元数值模拟参数标定. 山东工业技术, 2016(10): 42 https://www.cnki.com.cn/Article/CJFDTOTAL-SDGJ201610040.htm Sun J, He P S, Qi M J. Parameter calibration of discrete element numerical simulation by particle flow PFC software. Shandong Ind Technol, 2016(10): 42 https://www.cnki.com.cn/Article/CJFDTOTAL-SDGJ201610040.htm
[7] 徐泳, 孙其诚, 张凌, 等. 颗粒离散元法研究进展. 力学进展, 2003, 33(2): 251 DOI: 10.3321/j.issn:1000-0992.2003.02.010 Xu Y, Sun Q C, Zhang L, et al. Advances in discrete element methods for particulate materials. Adv Mech, 2003, 33(2): 251 DOI: 10.3321/j.issn:1000-0992.2003.02.010
[8] 程远方, 果世驹, 赖和怡. 球形颗粒随机排列过程的计算机模拟. 北京科技大学学报, 1999, 21(4): 387 DOI: 10.3321/j.issn:1001-053X.1999.04.020 Cheng Y F, Guo S J, Lai H Y. Computer simulation of random packing of spherical particles. J Univ Sci Technol Beijing, 1999, 21(4): 387 DOI: 10.3321/j.issn:1001-053X.1999.04.020
[9] 郑洲顺, 朱远鹏, 裴朝旭, 等. 高速压制成形中应力波传播的特征. 系统仿真学报, 2009, 21(增刊2): 226 https://www.cnki.com.cn/Article/CJFDTOTAL-XTFZ2009S2060.htm Zheng Z S, Zhu Y P, Pei Z X, et al. Characteristics of stress wave propagation in high velocity compaction. J Syst Simul, 2009, 21(Supple 2): 226 https://www.cnki.com.cn/Article/CJFDTOTAL-XTFZ2009S2060.htm
[10] 杜欣, 曾亚武, 高睿, 等. 用离散元方法研究颗粒外形对摩擦机理的影响. 西南交通大学学报, 2012, 47(2): 252 https://www.cnki.com.cn/Article/CJFDTOTAL-XNJT201202015.htm Du X, Ceng Y W, Gao R, et al. Effects of particle shape on friction mechanism based on discrete element method. J Southwest Jiao Tong Univ, 2012, 47(2): 252 https://www.cnki.com.cn/Article/CJFDTOTAL-XNJT201202015.htm
[11] 曹秒艳, 董国疆, 赵长财. 基于离散元法的固体颗粒介质传力特性研究. 机械工程学报, 2011, 47(14): 62 https://www.cnki.com.cn/Article/CJFDTOTAL-JXXB201114012.htm Cao M Y, Dong G J, Zhao C C. Research on pressure-transfer characteristics in the solid granule medium forming based on the discrete element method. J Mech Eng, 2011, 47(14): 62 https://www.cnki.com.cn/Article/CJFDTOTAL-JXXB201114012.htm
[12] 易明军, 尹海清, 曲选辉, 等. 力与应力波对高速压制压坯质量的影响. 粉末冶金技术, 2009, 27(3): 207 http://pmt.ustb.edu.cn/article/id/fmyjjs200903012 Yi M J, Yin H Q, Qu X H, et al. Influence of force and stress wave on the quality of green compacts in high velocity compaction. Powder Metall Technol, 2009, 27(3): 207 http://pmt.ustb.edu.cn/article/id/fmyjjs200903012
[13] Huang Z M, Ma Z G, Zhang L, et al. A numerical study of macro-mesoscopic mechanical properties of gangue backfill under biaxial compression. Int J Min Sci Technol, 2016, 26(2): 309 DOI: 10.1016/j.ijmst.2015.12.018
[14] 张炜, 周剑, 于世伟, 等. 离散元法金属粉末高速压制过程中力链特性量化研究. 机械工程学报, 2018, 54(10): 85 https://www.cnki.com.cn/Article/CJFDTOTAL-JXXB201810012.htm Zhang W, Zhou J, Yu S W, et al. Quantitative investigation on force chains of metal powder in high velocity compaction by using discrete element method. J Mech Eng, 2018, 54(10): 85 https://www.cnki.com.cn/Article/CJFDTOTAL-JXXB201810012.htm
[15] 孙其诚, 王光谦. 静态堆积颗粒中的力链分布. 物理学报, 2008, 57(8): 4667 DOI: 10.3321/j.issn:1000-3290.2008.08.007 Sun Q C, Wang G Q. Force distribution in static granular matter in two dimensions. Acta Phys Sin, 2008, 57(8): 4667 DOI: 10.3321/j.issn:1000-3290.2008.08.007
-
期刊类型引用(9)
1. 陈新宇,李奋强,蒋继帅. 数值模拟在粉末压制细观分析中的应用与发展. 粉末冶金技术. 2024(04): 418-426 . 本站查看
2. 张炜,萧伟健,袁传牛,张宁,陈荣昕. 基于三维离散元模型粉末压制中力链对阻塞行为的影响机制. 粉末冶金技术. 2024(04): 403-410+417 . 本站查看
3. 李化蓥,刘军,张超,张璐栋,王海陆,柯建忠. 基于离散元锥形零件冲击加载的相对密度. 粉末冶金技术. 2023(04): 322-329 . 本站查看
4. 王树强,陈钊,陈昊雷,周游. 机械镀锌过程中锌粉运动的模拟分析. 电镀与涂饰. 2022(01): 23-27 . 百度学术
5. 张炜,谈健君,张帅,陈荣昕,陈丙三,张宁,刘焜. 基于颗粒物质力学的铁粉末压制中摩擦特性对力链演化影响. 摩擦学学报. 2022(02): 386-395 . 百度学术
6. 薛亚军,贺福强,陈发江,徐浩然. 基于离散元法的木质板材压制成形过程仿真与分析. 林业机械与木工设备. 2021(05): 52-57 . 百度学术
7. 王海陆,刘军,林立,张超,张璐栋,柯建忠,李化蓥. 基于离散元的不同粒径配比粉末压制相对密度与力链分析. 粉末冶金技术. 2021(06): 490-498 . 本站查看
8. 刘义伦,曾洋. 不同成型工艺下钕铁硼模压成型过程的力学行为分析. 粉末冶金技术. 2020(04): 262-268 . 本站查看
9. 张璐栋,刘军,罗晓龙,张超,林立,周纯,王海陆. 基于离散元法的颗粒高速压制模拟及动态力学分析. 粉末冶金技术. 2020(05): 350-356+362 . 本站查看
其他类型引用(4)