Loading [MathJax]/jax/output/SVG/jax.js

大块高纯Fe3C的高温高压合成

杨俊 杨刚 陈星 赵彬 冷春蔚 刘勋 黄海军

奚小波, 缪宏, 赵荔, 朱河霖, 张瑞宏, 金亦富. 大耕深旋耕刀激光冲击强化残余应力研究[J]. 高压物理学报, 2015, 29(1): 52-58. doi: 10.11858/gywlxb.2015.01.009
引用本文: 杨俊, 杨刚, 陈星, 赵彬, 冷春蔚, 刘勋, 黄海军. 大块高纯Fe3C的高温高压合成[J]. 高压物理学报, 2019, 33(4): 043302. doi: 10.11858/gywlxb.20190715
XI Xiao-Bo, MIAO Hong, ZHAO Li, ZHU He-Lin, ZHANG Rui-Hong, JIN Yi-Fu. Study on Residual Stress of Big Tilling Depth Rotary Blade for Laser Shock Peening[J]. Chinese Journal of High Pressure Physics, 2015, 29(1): 52-58. doi: 10.11858/gywlxb.2015.01.009
Citation: YANG Jun, YANG Gang, CHEN Xing, ZHAO Bin, LENG Chunwei, LIU Xun, HUANG Haijun. Synthesis of Pure Bulk Fe3C under High Temperature and High Pressure[J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 043302. doi: 10.11858/gywlxb.20190715

大块高纯Fe3C的高温高压合成

doi: 10.11858/gywlxb.20190715
基金项目: 国家自然科学基金(41874103)
详细信息
    作者简介:

    杨 俊(1994-),男,硕士研究生,主要从事高压物理研究. E-mail:18207241437@163.com

    通讯作者:

    黄海军(1976-),男,博士,教授,主要从事高压物理与地球物理研究. E-mail:hjhuang@whut.edu.cn

  • 中图分类号: O521.3; P311

Synthesis of Pure Bulk Fe3C under High Temperature and High Pressure

  • 摘要: Fe3C(渗碳体)是一种用途广泛、力学和磁学性能优良的材料,因而其制备方法一直受到人们的关注。通过高温高压下铁和碳的固相反应,成功制备出大块、致密、高纯的Fe3C样品,探索了原料种类、颗粒大小、合成温度、合成压力、保温时间等因素对烧结样品的影响。结果表明:当Fe粉粒径为9 μm,石墨粒径为1.3 μm,在4 GPa、1000 ℃条件下烧结的Fe3C最为致密。

     

  • 旋耕刀是旋耕机作业最主要的受力部件, 也是最易受损的部件, 其性能直接影响旋耕机的作业效率和耕作质量。随着社会生产力的急剧发展, 大马力旋耕机的使用成为农业机械化发展的必然趋势, 而大耕深旋耕刀的加工制造成为关键技术难题。目前, 大耕深旋耕刀的制造工艺为同比例放大法, 即按比例将原有旋耕刀回转半径扩大至350 mm, 这就带来刀面是否加厚的问题。事实证明[1], 当刀面厚度不增加时, 机具耕作100 h内的刀具折断率高达30%以上; 当刀面厚度同比例增加时, 机具耕作100 h内的刀具折断率仍在20%左右, 且两种情况下后者的能耗增加了15%。为节省能耗、减少刀具材料使用, 不宜增加刀面厚度, 通过合适的表面处理提高刀具性能是一种有效方法[2-5]

    激光冲击强化(Laser Shock Peening, LSP)技术是一种新型的表面技术, 目前已广泛应用于材料表面改性的研究。激光冲击强化采用几十纳秒的短脉冲高峰值功率密度(> 109 W/cm2)激光辐射向金属表层, 使金属表面涂覆的保护层吸收激光能量并发生爆炸性气化蒸发, 产生大于1 GPa高压等离子体冲击波, 利用冲击波的力效应使表层材料微观组织发生变化引入残余压应力, 从而提高金属材料抗疲劳、耐磨损和防应力腐蚀等性能[6-8]

    旋耕刀切土时受到的土壤阻力大多集中在刀具的正切削刃上, 会导致旋耕刀刀柄处产生较大的集中应力, 刀柄在循环交变应力的作用下极易形成疲劳裂纹并失效[9-12]。为解决这一问题, 可在刀柄应力集中部位引入残余压应力以抵消旋耕刀作业时的外界应力, 抑制或减缓疲劳裂纹的产生, 提高刀具使用寿命。本研究借助ANSYS分析出旋耕刀应力集中区域, 并模拟激光冲击试验得出残余应力引入情况, 同时采用激光冲击强化技术处理应力集中部位, 利用X射线衍射法对残余应力进行测试, 讨论激光冲击强化处理后刀面的残余应力的影响。

    旋耕机作业时, 旋转刀对土壤进行切削、破碎及抛掷, 土壤对旋耕刀的反作用力构成了土壤阻力[13]。刀具在切削土壤时, 其耕作深度及切削土壤面均先由小到大, 后由大到小, 其在土壤中的位置也不断变化, 所以刀具受到土壤阻力的大小、方向和作用点在其切削土壤时是不断变化的。在切土过程中, 刀具基本不受轴向力的作用, 可将阻力F沿xz两个坐标轴分解成FxFz两分力, 则有

    {F=Fx+FzF=F2x+F2z (1)

    旋耕刀的回转半径为350 mm, 耕作深度为21 cm, 刀面切土宽度为6 cm, 土壤阻力计算公式为

    Rx=0.7abK (2)

    式中:Rx为土壤阻力, N;a为刀面切土宽度, cm;b为耕作深度, cm;K为土壤比阻, N·cm-2。一般土壤的K为4~5 N·cm-2, 粘土的K为6~8 N·cm-2, 本研究K取7 N·cm-2, 根据公式(2)可求出旋耕刀所受土壤阻力Rx约为620 N。

    在PRO/E 5.0中建立回转半径为350 mm、刀面厚度为7 mm的旋耕刀三维模型, 并将模型导入ANSYS中进行有限元分析。65Mn旋耕刀的弹性模量E为208 GPa, 泊松比n为0.288。采用Solid 168单元将模型自动划分网格, 共计86 021个节点、53 902个单元。在旋耕刀柄两个侧面及顶面与连接孔处施加约束, 并在受力处施加620 N的集中载荷, 最后运算求解。

    图 1(a)为旋耕刀的位移分布等值线图, 工作时旋耕刀的最大位移可以达到0.508 mm, 发生在旋耕刀切削土壤时的正切削刃和侧切削刃处, 而且位移变形是向两边等距递减的, 由于旋耕刀柄固定, 切土时正切削刃和侧切削刃受力, 因而此处变形较大。图 1(b)为旋耕刀的应力分布等值线图, 旋耕刀柄外弯角处应力最为集中, 最大拉应力达到303 MPa, 旋耕刀工作时该处极易折断。

    图  1  (a) 位移分布等值线图(b)应力分布等值线图
    Figure  1.  (a) Isoline of displacement distribution (b) Isoline of stress distribution

    旋耕刀采用65Mn为材料, 其主要化学成分为(质量分数, %):C 0.65, Si 0.26, Mn 1.10, S 0.002, P 0.019, Cr 0.02, Ni 0.01, Cu 0.01, Fe余量。材料在激光冲击过程中产生动态变形, 因此在选择材料应力-应变本构方程时, 静态方程已经不能反映材料的真实响应, 必须用动态本构方程替代。一般材料的动态本构方程, 目前较多采用Johnson-Cook方程[14], 其本构方程为

    σ=(A+Bεn)[1+Cln(˙ε˙ε0)][1(θθrθmθr)m] (3)

    式中:εθ分别为应变和温度为应变率分别为参考应变率和参考温度; θm为材料熔点; ABnCm为待定系数, ABn表征材料应变强化项系数, C表征材料应变率强化项系数, m表征材料热软化系数。

    段志勇等[15]在实验基础上, 提出了一个与实际比较相符的半经验型激光冲击波估计模型公式

    pmax=0.8ρI0 (4)

    式中:pmax为激光峰值压力, GPa; ρ为折合密度, g/cm3; I0为激光功率密度, GW/cm2。折合密度公式为

    2ρ=1ρ1+1ρ2 (5)

    式中:ρ1为约束层密度, g/cm3; ρ2为靶材密度, g/cm3。激光功率密度公式为

    I0=Eπr2τ (6)

    式中:E为激光能量, J;τ为激光脉宽, ns;r为光斑半径, mm。

    已知水密度ρ1=1 g/cm3, 靶材密度ρ2=7.85 g/cm3, 算得折合密度ρ=1.774 g/cm3; 激光能量E=10 J, 脉冲宽度τ=30 ns, 光斑半径r=3 mm, 算得激光功率密度I0=1.179 GW/cm2, 激光峰值压力pmax=1.16 GPa。激光冲击波作用时间为激光脉宽的3倍左右甚至更高[16], 实验采用的激光脉宽为30 ns, 取激光冲击波作用时间为100 ns, 其冲击波作用压力与时间的关系如图 2所示, 在50 ns时的压力最大, 最大压力值为1.16 GPa。

    图  2  激光冲击波加载曲线
    Figure  2.  Shock wave of laser loading curve

    采用Solid 70单元划分网格并定义材料属性, 设定环境温度为22 ℃, 设定对流及辐射边界条件, 并在应力集中区域的旋耕刀另一侧面设定绝热边界条件。在应力集中部位施加高斯表面热源, 由于区域面积较小, 高斯热源固定施加在待冲击区域中心, 设定载荷步选项和分析选项运算求解, 得到温度场模型。

    进行应力分析时, 通过前处理模块定义温度场, 按图 2中激光冲击波加载曲线定义载荷函数, 采用Solid 45单元重新划分网格并定义材料属性, 在应力集中区域定义加载单元组元, 将定义的冲击压力载荷施加在定义的单元组元。同时在单元组元对应的模型各个表面施加透射条件将应力波透射, 防止过大的应力波使模型发生弯扭变形影响仿真结果。最后施加约束并设置求解选项进行运算, 输出表面残余应力结果。

    激光冲击强化前, 将激光冲击强化区域进行抛光, 利用超声波清洗器将试样在乙醇中清洗15 min, 除去表面油污等杂质, 最后烘干试样[17]。采用铝箔作为激光冲击的金属涂敷层, 铝箔涂层长为25 mm, 宽为10 mm, 厚为0.1 mm。将铝箔粘贴在旋耕刀应力最集中的刀柄外弯角处, 如图 3(a)所示。使用专用夹具将贴有铝箔的旋耕刀夹于激光冲击强化处理工作台。用水作约束层, 水流厚度为0.5 mm。试验采用江苏大学Gaia-R系列高能量脉冲灯抽运YAG激光器, 激光器调整冲击参数为:激光能量10 J, 波长1.06 μm, 脉冲宽度30 ns, 频率0.1 Hz, 光斑直径6 mm。

    图  3  (a) 旋耕刀冲击强化前准备(b)残余应力测点位置
    Figure  3.  (a) The preparation before LSP (b) Residual stress measuring point location

    利用超声波清洗器将试样在丙酮溶液中清洗15 min, 并用酒精反复冲洗, 最后吹干[18-19]。采用MSF-3M型X射线应力测定仪进行残余应力测量, 应力测试取点如图 3(b)所示。X射线分析测试条件为V靶材, 所用衍射晶面为α-Fe(211)衍射晶面, 交相关方法定峰[20]

    从ANSYS后处理输出各测试点的表面残余应力仿真值, 对比试验实测数据, 如图 4所示。各测点的仿真值与实测值曲线呈上下交替分布, 曲线拟合较好, 误差在±20 MPa内。仿真值呈下凹状态, 主要是因为ANSYS中施加的高斯热源为固定热源, 且集中施加在冲击区域的中心, 造成中心温度向边缘扩散, 结果表现为整个区域中心的残余压应力明显高于周边。同理, 实测值呈现下凹状态也是因为激光冲击时, 光斑对准区域中心, 中心受热过大, 从而残余应力的引入明显多于周边。

    图  4  激光冲击仿真与实测残余应力
    Figure  4.  Residual stress datum of simulation and experiment after LSP

    利用ANSYS模拟激光冲击强化, 其在材料表层引入的残余应力与试验实测的结果吻合得较好, 说明ANSYS可完成激光冲击强化试验, 且误差在允许范围内, 通过此类仿真分析研究, 可对激光冲击强化的工艺参数进行优化, 进一步提高材料性能。

    图 5(a)、图 5(b)分别是激光冲击前后旋耕刀表面材料残余应力沿xy轴的分布情况, 未经强化的平均残余压应力仅为146.90 MPa, 强化后的平均残余压应力达到390.70 MPa, 提高了166%。经激光冲击强化后, 最大残余压应力达到412.25 MPa, 在冲击区域的中心位置; 最小残余压应力为375.82 MPa, 在冲击区域的最边缘位置。同时, 冲击后沿x轴及y轴方向的残余应力误差大小在±20 MPa内, 说明整个冲击强化区域的残余应力分布均匀, 无应力集中现象, 说明激光冲击强化能有效提高材料表层残余应力。

    图  5  冲击前、后的残余应力
    Figure  5.  Residual stress before and after LSP

    由于旋耕刀工作环境恶劣, 且工作时间较长, 易存在应力集中, 引发疲劳裂纹, 激光冲击强化处理引入的残余压应力可消除应力集中的影响, 减慢或抑制裂纹的扩展, 能起到提高板料疲劳寿命的作用, 进而提高旋耕刀使用寿命。

    (1) 利用ANSYS对旋耕刀进行应力分析, 得出旋耕刀柄外弯角处应力最为集中, 最大拉应力达到303 MPa, 与日常刀具易受损折断处的位置基本吻合。

    (2) ANSYS可模拟激光冲击强化工艺, 其算得的残余应力引入值与试验实测的结果吻合很好, 误差在±20 MPa内, 通过仿真分析可优化相关工艺参数, 进一步提高材料性能。

    (3) 经激光冲击强化处理的旋耕刀表层材料残余压应力明显增大, 最大残余压应力达412.25 MPa, 可消除旋耕刀在恶劣工况下应力集中的影响, 减慢或抑制裂纹的扩展, 提高材料疲劳寿命。

  • 图  Fe3C的晶体结构

    Figure  1.  Crystal structure of cementite

    图  样品组装示意图

    Figure  2.  Schematic of sample assembly

    图  烧结工艺曲线

    Figure  3.  Sintering process curve

    图  烧结样品的XRD谱

    Figure  4.  XRD patterns of sintered samples

    图  样品的背散射电子图和能谱分析图

    Figure  5.  Backscattered electron map and energy spectrum analysis of the sample

    图  碳的Kα峰值计数率随碳质量分数的变化(标准样品Fe3C和样品4中C的质量分数分别为(6.67±0.01)%和(6.58±0.01)%)

    Figure  6.  The Kα peak count rate of carbon varies with carbon mass fraction (The mass fraction of C in the standard sample Fe3C and Sample 4 are (6.67±0.01) % and (6.58±0.01)%, respectively.)

    图  Fe3C的横波声速(a)和纵波声速(b)的超声测量结果

    Figure  7.  (a) Shear wave and (b) compressional wave of Fe3C measured using ultrasonic equipment

    图  烧结样品的SEM图像

    Figure  8.  SEM image of sintered sample

    图  Fe3C的磁滞回线

    Figure  9.  Magnetic hysteresis loop diagram of Fe3C

    表  1  Fe3C原料及烧结参数

    Table  1.   Fe3C raw material and sintering parameters

    No.dFe/μmdC/μmT/℃t/minp/GPaProductDensity/(g·cm–3)
    114974 (Acticarbon) 800104Fe3C6.869
    21491.3 (Graphite) 800104Fe-C
    391.3 (Graphite) 800104Fe-C
    491.3 (Graphite)1000104Fe3C7.533
    50.031.3 (Graphite)1000104Fe3C6.931
    下载: 导出CSV
  • [1] 蔡红, 卢军, 王志明. 渗碳体的结构和一次渗碳体形态的研究 [J]. 热处理, 2010, 25(2): 50–53. doi: 10.3969/j.issn.1008-1690.2010.02.011

    CAI H, LU J, WANG Z M. Study on the structure of cementite and the morphology of primary cementite [J]. Heat Treatment, 2010, 25(2): 50–53. doi: 10.3969/j.issn.1008-1690.2010.02.011
    [2] 何强. 氧化物添加剂对两步法制备碳化铁的影响[D]. 武汉: 武汉科技大学, 2008: 3–5.

    HE Q. Effect of oxide additives on preparation of iron carbide by two-step process [D]. Wuhan: Wuhan University of Science and Technology, 2008: 3–5.
    [3] DAVID B, SCHNEEWEISS O, PIZÚROVÁ N, et al. Nanopowders with superparamagnetic Fe3C particles and their annealing behaviour [J]. Surface and Interface Analysis, 2010, 42(6/7): 699–702.
    [4] TERADA M, ITOH M, LIU J R, et al. Electromagnetic wave absorption properties of Fe3C/carbon nanocomposites prepared by a CVD method [J]. Journal of Magnetism and Magnetic Materials, 2009, 321(9): 1209–1213. doi: 10.1016/j.jmmm.2008.11.001
    [5] LIPERT K, KAZMIERCZAK J, PEŁECH I, et al. Magnetic properties of cementite (Fe3C) nanoparticle agglomerates in a carbon matrix [J]. Materials Science-Poland, 2007, 25(2): 399–404.
    [6] PRESCHER C, DUBROVINSKY L, BYKOVA E, et al. High Poisson’s ratio of Earth’s inner core explained by carbon alloying [J]. Nature Geoscience, 2015, 8(3): 220–223. doi: 10.1038/ngeo2370
    [7] WEERASINGHE G L, NEEDS R J, PICKARD C J. Computational searches for iron carbide in the Earth’s inner core [J]. Physical Review B, 2011, 84(17): 174110. doi: 10.1103/PhysRevB.84.174110
    [8] DODD S P, SAUNDERS G A, CANKURTARAN M, et al. Ultrasonic study of the temperature and hydrostatic-pressure dependences of the elastic properties of polycrystalline cementite (Fe3C) [J]. Physica Status Solidi (A), 2003, 198(2): 272–281. doi: 10.1002/(ISSN)1521-396X
    [9] DAVID B, ZBORIL R, MASHLAN M, et al. Single ferromagnetic behaviour of nanopowders with Fe3C [J]. Journal of Magnetism and Magnetic Materials, 2006, 304(2): e787–e789. doi: 10.1016/j.jmmm.2006.02.224
    [10] KOMOGORTSEV S V, ISKHAKOV R S, BALAEV A D, et al. Magnetic properties of Fe3C ferromagnetic nanoparticles encapsulated in carbon nanotubes [J]. Physics of the Solid State, 2007, 49(4): 734–738. doi: 10.1134/S1063783407040233
    [11] YANG K Y, XU W, ZHANG Y, et al. Synthesis and characteristics of Fe3C nanoparticles embedded in amorphous carbon matrix [J]. Chemical Research in Chinese Universities, 2010, 26(3): 348–351.
    [12] GAO L. Density, magnetic properties and sound velocities of iron-rich materials at high temperature and high pressure [D]. University of Illinois at Urbana-Champaign, 2011.
    [13] STEWART A J. Planetary cores: effect and behaviour of minor elements on the Fe-S system to 40 GPa [D]. Switzerland: ETH Zurich, 2006: 1–4.
    [14] DASGUPTA R, WALKER D. Carbon solubility in core melts in a shallow magma ocean environment and distribution of carbon between the Earth’s core and the mantle [J]. Geochimica et Cosmochimica Acta, 2008, 72(18): 4627–4641. doi: 10.1016/j.gca.2008.06.023
  • 加载中
图(9) / 表(1)
计量
  • 文章访问数:  5484
  • HTML全文浏览量:  3391
  • PDF下载量:  50
出版历程
  • 收稿日期:  2019-01-18
  • 修回日期:  2019-03-26

目录

/

返回文章
返回