
Citation: | ZHAO Kanglin, WANG Qiming, ZHANG Youjun, JIANG Gang, PENG Fang, LI Yanchun. Physico-Mechanical Behavior and Size Effect of Nano-Tungsten under High Pressure[J]. Chinese Journal of High Pressure Physics, 2024, 38(3): 030103. doi: 10.11858/gywlxb.20230756 |
旋耕刀是旋耕机作业最主要的受力部件, 也是最易受损的部件, 其性能直接影响旋耕机的作业效率和耕作质量。随着社会生产力的急剧发展, 大马力旋耕机的使用成为农业机械化发展的必然趋势, 而大耕深旋耕刀的加工制造成为关键技术难题。目前, 大耕深旋耕刀的制造工艺为同比例放大法, 即按比例将原有旋耕刀回转半径扩大至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沿x、z两个坐标轴分解成Fx、Fz两分力, 则有
{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, 旋耕刀工作时该处极易折断。
旋耕刀采用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为材料熔点; A、B、n、C、m为待定系数, A、B、n表征材料应变强化项系数, 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。
采用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。
利用超声波清洗器将试样在丙酮溶液中清洗15 min, 并用酒精反复冲洗, 最后吹干[18-19]。采用MSF-3M型X射线应力测定仪进行残余应力测量, 应力测试取点如图 3(b)所示。X射线分析测试条件为V靶材, 所用衍射晶面为α-Fe(211)衍射晶面, 交相关方法定峰[20]。
从ANSYS后处理输出各测试点的表面残余应力仿真值, 对比试验实测数据, 如图 4所示。各测点的仿真值与实测值曲线呈上下交替分布, 曲线拟合较好, 误差在±20 MPa内。仿真值呈下凹状态, 主要是因为ANSYS中施加的高斯热源为固定热源, 且集中施加在冲击区域的中心, 造成中心温度向边缘扩散, 结果表现为整个区域中心的残余压应力明显高于周边。同理, 实测值呈现下凹状态也是因为激光冲击时, 光斑对准区域中心, 中心受热过大, 从而残余应力的引入明显多于周边。
利用ANSYS模拟激光冲击强化, 其在材料表层引入的残余应力与试验实测的结果吻合得较好, 说明ANSYS可完成激光冲击强化试验, 且误差在允许范围内, 通过此类仿真分析研究, 可对激光冲击强化的工艺参数进行优化, 进一步提高材料性能。
图 5(a)、图 5(b)分别是激光冲击前后旋耕刀表面材料残余应力沿x、y轴的分布情况, 未经强化的平均残余压应力仅为146.90 MPa, 强化后的平均残余压应力达到390.70 MPa, 提高了166%。经激光冲击强化后, 最大残余压应力达到412.25 MPa, 在冲击区域的中心位置; 最小残余压应力为375.82 MPa, 在冲击区域的最边缘位置。同时, 冲击后沿x轴及y轴方向的残余应力误差大小在±20 MPa内, 说明整个冲击强化区域的残余应力分布均匀, 无应力集中现象, 说明激光冲击强化能有效提高材料表层残余应力。
由于旋耕刀工作环境恶劣, 且工作时间较长, 易存在应力集中, 引发疲劳裂纹, 激光冲击强化处理引入的残余压应力可消除应力集中的影响, 减慢或抑制裂纹的扩展, 能起到提高板料疲劳寿命的作用, 进而提高旋耕刀使用寿命。
(1) 利用ANSYS对旋耕刀进行应力分析, 得出旋耕刀柄外弯角处应力最为集中, 最大拉应力达到303 MPa, 与日常刀具易受损折断处的位置基本吻合。
(2) ANSYS可模拟激光冲击强化工艺, 其算得的残余应力引入值与试验实测的结果吻合很好, 误差在±20 MPa内, 通过仿真分析可优化相关工艺参数, 进一步提高材料性能。
(3) 经激光冲击强化处理的旋耕刀表层材料残余压应力明显增大, 最大残余压应力达412.25 MPa, 可消除旋耕刀在恶劣工况下应力集中的影响, 减慢或抑制裂纹的扩展, 提高材料疲劳寿命。
[1] |
CORDERO Z C, KNIGHT B E, SCHUH C A. Six decades of the Hall-Petch effect: a survey of grain-size strengthening studies on pure metals [J]. International Materials Reviews, 2016, 61(8): 495–512. doi: 10.1080/09506608.2016.1191808
|
[2] |
SCHIØTZ J, DI TOLLA F D, JACOBSEN K W. Softening of nanocrystalline metals at very small grain sizes [J]. Nature, 1998, 391(6667): 561–563. doi: 10.1038/35328
|
[3] |
YIP S. The strongest size [J]. Nature, 1998, 391(6667): 532–533. doi: 10.1038/35254
|
[4] |
SCHIØTZ J, JACOBSEN K W. A maximum in the strength of nanocrystalline copper [J]. Science, 2003, 301(5638): 1357–1359. doi: 10.1126/science.1086636
|
[5] |
ZHOU X L, FENG Z Q, ZHU L L, et al. High-pressure strengthening in ultrafine-grained metals [J]. Nature, 2020, 579(7797): 67–72. doi: 10.1038/s41586-020-2036-z
|
[6] |
GU Q F, KRAUSS G, STEURER W, et al. Unexpected high stiffness of Ag and Au nanoparticles [J]. Physical Review Letters, 2008, 100(4): 045502. doi: 10.1103/PhysRevLett.100.045502
|
[7] |
YANG J, DENG W, LI Q, et al. Strength enhancement of nanocrystalline tungsten under high pressure [J]. Matter and Radiation at Extremes, 2020, 5(5): 058401. doi: 10.1063/5.0005395
|
[8] |
HE D W, DUFFY T S. X-ray diffraction study of the static strength of tungsten to 69 GPa [J]. Physical Review B, 2006, 73(13): 134106. doi: 10.1103/PhysRevB.73.134106
|
[9] |
XIONG L, LI B, TANG Y M, et al. Radial X-ray diffraction study of the static strength and texture of tungsten to 96 GPa [J]. Solid State Communications, 2018, 269: 83–89. doi: 10.1016/j.ssc.2017.10.016
|
[10] |
邹勇, 陈立溁. 贵金属金的状态方程 [J]. 高压物理学报, 2006, 20(3): 308–312. doi: 10.3969/j.issn.1000-5773.2006.03.015
ZOU Y, CHEN L Y. Equation of state for gold [J]. Chinese Journal of High Pressure Physics, 2006, 20(3): 308–312. doi: 10.3969/j.issn.1000-5773.2006.03.015
|
[11] |
MUSTAPHA S, NDAMITSO M M, ABDULKAREEM A S, et al. Comparative study of crystallite size using Williamson-Hall and Debye-Scherrer plots for ZnO nanoparticles [J]. Advances in Natural Sciences: Nanoscience and Nanotechnology, 2019, 10(4): 045013. doi: 10.1088/2043-6254/ab52f7
|
[12] |
WU Z Q, LIN F. Evaluation of Pt and Au pressure scales based on MgO absolute pressure scale [J]. Science China Earth Sciences, 2017, 60(1): 114–123. doi: 10.1007/s11430-015-0232-4
|
[13] |
WANG Y, AHUJA R, JOHANSSON B. Reduction of shock-wave data with mean-field potential approach [J]. Journal of Applied Physics, 2002, 92(11): 6616–6620. doi: 10.1063/1.1518781
|
[14] |
HAMMERSLEY A P, SVENSSON S O, HANFLAND M, et al. Two-dimensional detector software: from real detector to idealised image or two-theta scan [J]. High Pressure Research, 1996, 14(4): 235–248. doi: 10.1080/08957959608201408
|
[15] |
SINGH T B, REY L, GARTIA R K. Applications of PeakFit software in thermoluminescence studies [J]. Indian Journal of Pure & Applied Physics, 2011, 49(5): 297–302.
|
[16] |
李春, 刘泽民, 陈恒杰, 等. PeakFit、1stOpt在弗兰克-赫兹实验数据处理中的应用 [J]. 大学物理实验, 2018, 31(5): 117–123. doi: 10.14139/j.cnki.cn22-1228.2018.05.028
LI C, LIU Z M, CHEN H J, et al. Application of PeakFit and 1stOpt in data processing of the Frank-Hertz experiment [J]. Physical Experiment of College, 2018, 31(5): 117–123. doi: 10.14139/j.cnki.cn22-1228.2018.05.028
|
[17] |
LANGFORD J. X-ray powder diffraction studies of vitromet samples [J]. Journal of Applied Crystallography, 1971, 4: 164–168. doi: 10.1107/S002188987100654X
|
[18] |
余建新, 王晓鹏, 崔喜平. 高温环境下材料泊松比测试方法研究[J]. 实验科学与技术, 2022, 20(1): 28–33.
YU J X, WANG X P, CUI X P, Material poisson’s ratio measurement method at elevated temperatures [J]. Experiment Science and Technology, 2022, 20(1): 28–33.
|
[19] |
EINARSDOTTER K, SADIGH B, GRIMVALL G, et al. Phonon instabilities in fcc and bcc tungsten [J]. Physical Review Letters, 1997, 79(11): 2073–2076. doi: 10.1103/PhysRevLett.79.2073
|
[20] |
KATAHARA K, MANGHNANI M, FISHER E. Pressure derivatives of the elastic-moduli of bcc Ti-V-Cr, Nb-Mo and Ta-W alloys [J]. Journal of Physics F: Metal Physics, 1979, 9(5): 773–790. doi: 10.1088/0305-4608/9/5/006
|
[21] |
ZHANG J Z, ZHAO Y S, PALOSZ B. Comparative studies of compressibility between nanocrystalline and bulk nickel [J]. Applied Physics Letters, 2007, 90(4): 043112. doi: 10.1063/1.2435325
|
[22] |
DEWAELE A, LOUBEYRE P, MEZOUAR M. Equations of state of six metals above 94 GPa [J]. Physical Review B, 2004, 70(9): 094112. doi: 10.1103/PhysRevB.70.094112
|
[23] |
HIXSON R S, FRITZ J N. Shock compression of tungsten and molybdenum [J]. Journal of Applied Physics, 1992, 71(4): 1721–1728. doi: 10.1063/1.351203
|
[24] |
CHIJIOKE A D, NELLIS W J, SILVERA I F. High-pressure equations of state of Al, Cu, Ta, and W [J]. Journal of Applied Physics, 2005, 98(7): 073526. doi: 10.1063/1.2071449
|
[25] |
MING L C, MANGHNANI M H. Isothermal compression of bcc transition metals to 100 kbar [J]. Journal of Applied Physics, 1978, 49(1): 208–212.
|
[26] |
MA Y Z, CUI Q L, SHEN L H, et al. X-ray diffraction study of nanocrystalline tungsten nitride and tungsten to 31 GPa [J]. Journal of Applied Physics, 2007, 102(1): 013525.
|
[27] |
CHEN B, PENWELL D, BENEDETTI L R, et al. Particle-size effect on the compressibility of nanocrystalline alumina [J]. Physical Review B, 2002, 66(14): 144101. doi: 10.1103/PhysRevB.66.144101
|
[28] |
GLEASON A E, MARQUARDT H, CHEN B, et al. Anomalous sound velocities in polycrystalline MgO under non-hydrostatic compression [J]. Geophysical Research Letters, 2011, 38: L03304.
|
[29] |
LE GODEC Y, KURAKEVYCH O O, MUNSCH P, et al. Effect of nanostructuration on compressibility of cubic BN [J]. Journal of Superhard Materials, 2012, 34(5): 336–338. doi: 10.3103/S1063457612050085
|
[30] |
AL-KHATATBEH Y, LEE K K M, KIEFER B. Compressibility of nanocrystalline TiO2 anatase [J]. Journal of Physical Chemistry C, 2012, 116(40): 21635–21639. doi: 10.1021/jp3075699
|
[31] |
JIANG J Z, OLSEN J S, GERWARD L, et al. Enhanced bulk modulus and reduced transition pressure in γ-Fe2O3 nanocrystals [J]. Europhysics Letters, 1999, 45(2): 275–275. doi: 10.1209/epl/i1999-00158-3
|
[32] |
TOLBERT S H. High-pressure structural transformations in semiconductor nanocrystals [J]. Annual Review of Physical Chemistry, 1995, 46(1): 595–625. doi: 10.1146/annurev.pc.46.100195.003115
|
[33] |
HONG X G, DUFFY T S, EHM L, et al. Pressure-induced stiffness of Au nanoparticles to 71 GPa under quasi-hydrostatic loading [J]. Journal of Physics: Condensed Matter, 2015, 27(48): 485303. doi: 10.1088/0953-8984/27/48/485303
|
[34] |
MIKHEYKIN A S, DMITRIEV V P, CHAGOVETS S V, et al. The compressibility of nanocrystalline Pt [J]. Applied Physics Letters, 2012, 101(17): 173111. doi: 10.1063/1.4758000
|
[35] |
CHEN B, PENWELL D, KRUGER M B, et al. Nanocrystalline iron at high pressure [J]. Journal of Applied Physics, 2001, 89(9): 4794–4796. doi: 10.1063/1.1357780
|
[36] |
WANG Z W, PISCHEDDA V, SAXENA S K, et al. X-ray diffraction and Raman spectroscopic study of nanocrystalline CuO under pressures [J]. Solid State Communications, 2002, 121(5): 275–279. doi: 10.1016/S0038-1098(01)00509-9
|
[37] |
GU Q F, KRAUSS G, GRAMM F, et al. On the compressibility of TiC in microcrystalline and nanoparticulate form [J]. Journal of Physics: Condensed Matter, 2008, 20(44): 445226. doi: 10.1088/0953-8984/20/44/445226
|
[38] |
WANG Q M, HE D W, PENG F, et al. Compression behavior of nanocrystalline TiN [J]. Solid State Communications, 2014, 182: 26–29. doi: 10.1016/j.ssc.2013.12.015
|