Prediction and optimization of tensile strength of 7A52 aluminum alloy friction stir welding joints based on response surface methodology
-
摘要: 为了研究7A52铝合金搅拌摩擦焊的焊接速度、搅拌头转速及轴肩压深对接头抗拉强度的影响,采用响应面法的中心复合试验设计法设计20组试验,并建立抗拉强度响应函数关系式. 为了验证响应函数关系式的精确性,通过方差分析和回归分析确定该回归模型为显性,相关性系数R2的偏差为3.17%. 通过单一焊接参数因素和双因素焊接参数对抗拉强度的影响分析,进一步验证了模型的准确性,最后通过拉伸试验验证. 结果表明,基于响应面法拟合的搅拌摩擦焊焊接速度、搅拌头转速及轴肩压深与接头抗拉强度响应函数关系式能精确的预算不同焊接参数组合所对应的接头抗拉强度,并获得接头最佳参数组合为焊接速度110 mm/min、搅拌头转速1 436 r/min和轴肩压深0.55 mm,得到最大预测抗拉强度为380 MPa.Abstract: In order to study the effects of welding speed, stirring head rotation speed and pressure deep of shaft shoulder on tensile strength of 7A52 aluminum alloy friction stir welding. 20 groups of tests were designed by response surface methodology based on central composite test design, and response function relationship were established. In order to verify the accuracy of the response function relationship, variance analysis and regression analysis were used to determine the dominance of the regression model, and the deviation of correlation coefficient R2 was only 3.17%. The accuracy of the model was verified by analyzing the influence of single welding parameter and double welding parameter on tensile strength. Finally, the model was verified by tensile test. The results show that the joint tensile strength can be predicated based on response function relationship of response surface methodology fitting, and the best combination of welding parameter (welding speed 110 mm/min, stirring head rotation speed 1 436 r/min, pressure deep of shaft shoulder 0.55 mm) was gained. The maximum predication tensile strength was 380 MPa.
-
0. 序言
激光焊接由于其焊接速度快、热输入量低、焊接变形小等优势,已经广泛应用于汽车制造、航空航天、能源装备及船舶制造等工业生产焊接制造过程中. 激光焊接过程涉及到材料对激光能量的吸收过程、材料的熔化与蒸发、等离子体(羽烟)的形成、熔池与匙孔的形成、凝固与相变等复杂的多物理场问题,建立合理的反映激光焊接热过程特点的数值模型,定量描述和澄清激光焊接过程中的熔化、蒸发、凝固等物理现象,对于更好地理解和掌握激光热过程对焊缝成形、组织及应力演化的影响规律及内在机理十分重要,故数值模拟方法是焊接关键问题科学研究的强有力分析手段[1].
激光焊接自身快速熔凝的特点导致其焊接过程中更易产生驼峰、气孔、飞溅、侧壁未熔合等焊接缺陷[2-5],直接影响激光焊缝成形质量,另外激光焊缝组织梯度大及残余应力等问题也成为影响激光焊接构件服役性能的重要因素.
文中从影响激光焊缝成形质量及服役性能的焊接缺陷、组织演变、残余应力几个方面入手,阐述数值模拟技术在激光焊接基础问题方面应用及发展.
1. 激光焊缝成形问题的数值模拟
激光焊缝成形过程及相关缺陷的产生与熔池的流动行为、匙孔的动态波动及熔池的凝固过程密切相关. 此部分主要是针对熔池流动及匙孔行为进行阐述.
1.1 焊缝表面成形
在激光高速焊接时,焊缝表面成形主要涉及驼峰、咬边、未熔合、焊穿等问题,这与熔池流动及匙孔行为密切相关. 在表面驼峰成形方面,Otto等人[6]模拟研究了激光焊接表面驼峰的形成过程,发现金属蒸气的剪切力导致熔化金属向后加速,形成的熔体支流沿侧壁向后铺展,两侧支流汇聚形成驼峰.
在激光熔透焊接时,根部焊缝成形控制非常重要. Zhang等人[7]利用自行研发的三维激光焊接模型,利用锐利界面方法处理了影响匙孔动态行为的关键物理因素,发现小孔根部的流体受到强烈的蒸气反作用力导致根部驼峰的形成,如图1所示. Bachmann等人[8-9]研究表明激光熔透焊接过程中电磁力可以改变熔池的流动方向,进而阻止熔池下坠,利于熔透焊缝根部成形.
在激光焊接咬边、焊穿等成形方面,Cho等人[10]发现增加激光束震动频率改变了能量的传输及熔池的流动行为,从而明显改善焊缝成形质量,咬边和焊穿情况也得到明显改善,如图2所示.
1.2 焊接气孔
在焊缝气孔缺陷形成方面,许多学者分别在气泡形成、影响因素及气孔的演化等方面做了大量的数值模拟工作来研究气孔形成过程及其抑制机理. Zhao等人[11]认为当匙孔深度突然变大时容易造成匙孔坍塌,气体进入熔池形成气泡. Pang等人[12-13]自行开发了考虑匙孔、金属蒸气羽烟及熔池动态行为的三维多相自洽模型,发现蒸气羽烟摆动频率与匙孔波动频率近乎一致,阐明了蒸气羽烟形成的负压导致了匙孔后壁形成涡流,从而诱发环境气体卷入匙孔,促进了气泡的形成,如图3所示.
气泡能否演化为气孔涉及到熔池流动行为及凝固过程,Lu等人[14-15]模拟发现气泡迁移过程中存在消失、合并等现象,且被凝固界面捕获后才会演变成气孔,如图4所示. 另外低沸点元素蒸发对匙孔动态行为有着重要的影响,Huang等人[16-17]模拟考虑了Mg元素含量对激光焊接铝合金气孔生成的影响,认为Mg元素含量增加使得蒸气反作用力增加,从而使得匙孔更加稳定,降低了匙孔型气孔的生成. 激光震荡在消除气孔方面的作用机理也得到了模拟澄清[18-19],其根本原因除了激光震荡过程可以提高气泡逃逸速度外,模拟还发现气泡在震荡过程中被匙孔合并导致气泡消失,进而消除气孔的机理.
1.3 焊接飞溅
飞溅的产生会带来金属的元素损失,并影响焊缝成形质量. Hugger等人[20]模拟了钢、铝激光焊接过程,发现熔体首先在匙孔的边缘产生驼峰,在表面张力作用下颈缩,在冲击力的作用下离开表面形成飞溅,如图5所示. Chang等人[21]研究认为铝合金激光焊接时飞溅的产生与熔池的流动速度和涡流密切相关. Wu等人[22]模拟发现匙孔后壁的熔池沿熔合线向上流动,导致上表面飞溅产生.由于飞溅产生的影响因素较为复杂,考虑金属蒸气行为对于理解飞溅产生的过程十分必要.
激光焊接过程中低沸点元素蒸发诱发飞溅产生的情况,也引起了学者的兴趣. Hao等人[23-25]模拟研究了锌蒸气对匙孔与熔池动态行为的影响,分析了飞溅产生的过程,如图6所示,模拟还发现激光折线扫描路径下,由于锌蒸气分阶段进入匙孔,冲击力减小,使得匙孔和熔池更加稳定,明显降低了飞溅情况. Qi等人[26]模拟环形激光热源作用下镀锌钢板的飞溅情况,认为产生飞溅的主要原因是熔池中的回流和涡流导致,同轴环形和高斯热源的匹配增加了回流的阻力,进而控制了飞溅的产生.
2. 激光焊缝组织演化数值模拟
宏、微观组织的模拟也是激光焊接数值模拟的重要组成部分,近年来,各国学者已对激光焊缝组织的仿真进行了相关尝试,本节主要针对这方面的研究进展进行介绍.
2.1 宏观偏析
目前,针对焊缝宏观组织的数值模拟主要集中于异种材料因元素混合不均匀引起的宏观偏析方面. Gu等人[27]利用元素分离法实现了不同材料在激光选区熔化过程中的界面相迁移行为的模拟. Yao等人[28]利用该方法,并结合微尺度流体动力学和纳秒级别的热扩散分析,实现了316L不锈钢/Inconel 718异种材料激光熔池中鱼鳞状宏观偏析组织的模拟,如图7所示,这种宏观偏析组织的仿真重点在于界面参数的优化.
2.2 微观组织
焊缝微观组织的数值模拟的研究主要有以下3种方法:相场法、元胞自动机法以及蒙特卡罗法,且每种方法的侧重点不同.
Fallah等人[29]采用相场法,针对不同条件下Ti-Nb合金激光熔池中的一次枝晶臂间距进行了模拟计算. Mi等人[30]同样利用这一方法对Al-Cu合金激光焊熔池中的竞争生长行为进行了研究,发现在非最优取向的晶粒会影响到相邻晶粒的生长,晶界处的枝晶更易于在横向生长,通过阻碍其它方向枝晶的生长路径,进而改变不同取向晶粒的尺寸. Geng等人[31-32]发现Al-Mg较Al-Cu合金会在凝固过程中更早地发生广泛的枝晶连接现象,从而使其裂纹敏感性更低,这可以归因于其一次枝晶较低的溶质偏析,如图8a,8b所示.
图 8 溶质分布及晶粒形貌模拟结果Figure 8. Simulation results of solute distribution and microscopic morphologies. (a) solute distribution in columnar growth of Al-4.0%Cu alloy; (b) solute distribution in columnar growth of Al-4.0%Mg alloy; (c) microscopic morphologies with different pre-heating temperatures in a single track相较于相场法,元胞自动机法以其计算效率高、所需的模型参数少、试验匹配度较好的特点,得到了越来越多的关注. Ao等人[33]利用该方法对激光熔池的微观组织进行了模拟,结果如图8c所示,可见等轴晶比例会随着预热温度的增加和激光扫描速度的降低而显著增加,而等轴晶形貌会随着相邻两道激光距离的增加变得更加长而窄. Shi等人[34]利用该模型结合有限元的方法,实现了含气孔的不同柱状晶和等轴晶比例的熔池的仿真,并发现熔池晶粒生长形态取决于熔池的宽度、深度以及过冷度. 近期,Liu等人[35]提出了结合三维元胞自动机法和一维相场法的新模型,通过结合元胞自动机法的高效性和相场法的准确性,成功实现了大尺度铝合金凝固前沿枝晶形貌的预测.
对于蒙特卡罗法,该模型并不考虑凝固过程中的枝晶形貌、元素偏聚和过冷等因素,目前被广泛地用于晶粒长大、再结晶等晶粒演化行为的数值模拟. Zhang等人[36]利用该模型成功实现了12%Cr铁素体不锈钢在激光-电弧复合焊条件下热影响区晶粒的模拟,并发现顶部区域热影响区的晶粒尺寸要大于根部. Gleason等人[37-38]利用该方法针对激光冲击焊接过程中1100铝合金和304不锈钢的不均匀的界面微观组织进行了仿真,如图9所示,实现了对不同屈服面、层错能以及晶界滑动条件下的晶粒伸长行为的预测.
3. 焊接残余应力数值模拟
与传统的弧焊相比,激光焊对焊接残余应力和变形的影响主要体现在制造工艺方面[39],以下将对近几年来在激光焊接残余应力方面取得的进展进行介绍. Sun等人[40]以Q235钢平板对接接头为研究对象,采用试验手段和热-弹-塑性有限元方法研究了激光焊和电弧焊的温度场、残余应力与焊接变形. 图10是激光焊(Case A)与电弧焊(Case C)得到的焊缝中央位置的温度循环曲线,可见激光焊的熔化面积小,而且熔深贯穿整个板厚,而电弧焊尽管熔化面积更大,激光焊的加热速度更快,峰值温度更高,高温停留时间更短而且冷却速度也更快.
图11是薄板接头中央断面上的纵向残余应力与横向残余应力的数值模拟结果,可见在激光焊条件下接头的上、下表面的纵向残余应力的大小与分布几乎完全一致,这是因为激光焊产生的熔化区域在板厚方向分布较均匀所致,拉伸残余应力峰值与材料常温屈服强度基本一致. 在电弧焊条件下的横向残余应力的峰值要远高于激光焊,单就焊接残余应力的峰值大小而言,激光焊并不能有效减缓残余应力,该研究结果也得到了其它学者的支持[41-42].
在激光厚板焊接残余应力方面,Xu等人[43]模拟比较了板厚为12 mm的Q460钢在激光焊和多层多道气体保护焊条件下接头的残余应力,发现激光焊接头的高拉伸纵向残余应力范围要明显小于气体保护焊接头,但峰值应力没有明显的差异. Yan等人[44]研究了板厚为10 mm 的316L不锈钢在激光条件下的残余应力,如图12所示,由于板厚较厚,横向残余应力的大小与薄板接头相比有显著增加. Deng等人[45]采用考虑加工硬化及退火软化的热-弹-塑性有限元方法研究了SUS304管-管激光焊对接接头的残余应力,并讨论了由固溶处理引起的初期残余应力对焊接残余应力的影响.
Elmesalamy等人[46]以板厚为10和20 mm 316L钢对接接头为研究对象,采用轮廓法测量了多道窄间隙激光焊对接接头和多层多道TIG焊对接接头的残余应力. 图13是距焊接接头上表面1.5 mm位置的纵向残余应力分布,可见窄间隙激光焊和多层多道TIG焊产生的峰值应力前者为310 MPa,而后者为520 MPa,且激光焊纵向高拉伸残余应力区域明显较窄,故窄间隙激光焊方法可能是一个控制残余应力的有效途径,且多道窄间隙激光焊接方法将是用于控制加工硬化显著材料厚板接头焊接残余应力一种有效的方法.
近年来,激光焊残余应力的报道在铝合金、镁合金、钛合金[47]以及异种金属激光熔钎焊[48]方面也越来越多.
4. 结束语与展望
(1) 在激光焊缝成形模拟方面,匙孔的动态行为决定了焊接过程的稳定性,现有数值模型已经实现了匙孔壁与激光能量的实时耦合,但是激光束与等离子体(羽烟)的耦合行为很少涉及,考虑金属蒸气流动行为与激光热源的相互作用,是激光焊接过程建模发展的方向.
(2) 在激光焊缝组织模拟方面,熔池的流动对元素宏观偏析行为方面的模拟结果取决于界面参数的选取与优化;而微观组织方面,进一步完善3种模型需要准确获得新晶粒形核的相关参数十分重要,实现宏观与微观组织多尺度模拟还具有很高的挑战性.
(3) 在激光焊接残余应力方面,对于低合金高强钢及超高强钢而言,如何建立完备的“热-组织-力学”多场耦合模型,高精度的模拟激光焊条件下残余应力将是一个挑战性的课题. 此外,对加工硬化较显著材料优化工艺降低残余应力、降低累积塑性应变和减缓敏化程度是值得深入探索的方向.
-
表 1 7A52铝合金的化学成分(质量分数,%)
Table 1 Chemical compositions of 7A52 aluminum alloy
Zr Zn Mg Fe Mn Si Cr Cu Al 0.12 4.2 2.3 0.33 0.31 0.24 0.23 0.13 余量 表 2 焊接工艺参数
Table 2 Welding parameter
焊接速度
v/(mm·min−1)转速
n/(r·min−1)轴肩压深
h/mm60 ~ 200 600 ~ 2 500 0.3 ~ 0.8 表 3 基于CCD试验设计的FSW接头抗拉强度
Table 3 Results of tensile strength of FSW joints based on CCD
因素 实际焊接参数 标准化焊接参数 抗拉强度
Rm/MPa焊接速度v/(mm·min−1) 转速n/(r·min−1) 轴肩压深h/mm 焊接速度X1 转速X2 轴肩压深X3 1 70 700 0.3 −1 −1 −1 275.10 2 70 700 0.8 −1 −1 1 328.92 3 70 1500 0.3 −1 1 −1 251.61 4 70 1500 0.8 −1 1 1 305.43 5 150 700 0.3 1 −1 −1 167.58 6 150 700 0.8 1 −1 1 221.40 7 150 1500 0.3 1 1 −1 247.96 8 150 1 500 0.8 1 1 1 301.77 9 42.73 1 100 0.55 −1.68 0 0 267.45 10 177.27 1 100 0.55 1.68 0 0 266.72 11 110 427.28 0.55 0 −1.68 0 260.62 12 110 1 772.72 0.55 0 1.68 0 371.91 13 110 1 100 0.13 0 0 −1.68 141.16 14 110 1 100 0.97 0 0 1.68 231.66 15 110 1 100 0.55 0 0 0 368.89 16 110 1 100 0.55 0 0 0 378.68 17 110 1 100 0.55 0 0 0 348.72 18 110 1 100 0.55 0 0 0 388.89 19 110 1 100 0.55 0 0 0 375.59 20 110 1 100 0.55 0 0 0 386.58 表 4 标准条件下回归系数及P值
Table 4 Coefficients and their P values in coded condition
项目 预测系数 误差e F值 Prob > F值 X1 −16.371 1 6.236 805 6.890 166 0.025 4 X2 22.035 61 6.236 805 12.483 19 0.005 4 X3 26.907 62 6.236 805 18.613 42 0.001 5 X12 −35.301 9 6.071 366 33.808 2 0.000 2 X1 X2 25.966 25 8.148 775 10.153 92 0.009 7 X1 X3 −0.001 25 8.148 775 2.356 53 0.999 9 X22 −17.914 1 6.071 366 8.705 962 0.014 5 X2X3 −0.001 25 8.148 775 2.354 33 0.999 9 X32 −63.824 9 6.071 366 110.511 3 < 0.000 1 表 5 模型的方差分析结果
Table 5 Results of ANOVA test
项目 自由度 f 平方和 SS 均方值 MS F值 Prob > F值 X1 1 3 660.20 3660.20 8.268 0.014 X2 1 6 631.33 6631.33 14.980 0.002 X3 1 9 887.83 9887.83 22.336 0.001 X12 1 17 959.60 17959.60 40.570 <0 X1X2 1 5 393.97 5393.97 12.185 0.005 X22 1 4 624.78 4624.78 10.447 0.007 X32 1 58 705.87 58705.87 132.614 <0 模型 7 97 667.90 13952.56 31.518 <0 表 6 模型预测值与试验值比较
Table 6 Comparison of model predicted value and experimental value
因素 转速
n/(r·min−1)焊接速度
v/(mm·min−1)轴肩压深
h/mm抗拉强度Rm/MPa 误差e(%) 预测值 试验值 预测值 1100 110 0.55 374 369 1.33 优化值 1436 110 0.55 380 377 0.79 -
[1] 陈芙蓉, 贾翠玲. 7A52铝合金焊接及其接头表面纳米化研究现状[J]. 华东交通大学学报, 2019, 36(1): 1 − 11. Chen Furong, Jia Cuiling. Review on 7A52 aluminum alloy welding and its welded joint surface nanocrystallization[J]. Journal of East China Jiaotong University, 2019, 36(1): 1 − 11.
[2] Jia Y, Qin Y, Ou Y, Wang K, et al. The influence of microstructural heterogeneity on mechanical properties of friction stir welded joints of T6-treated Al-Zn-Mg alloy 7A52[J]. Metals, 2018, 8: 527 − 538. doi: 10.3390/met8070527
[3] Eivani A R, Vafaeenezhad H, Jafarian H R, et al. A novel approach to determine residual stress field during FSW of AZ91 Mg alloy using combined smoothed particle hydrodynamics/ neuro-fuzzy computations and ultrasonic testing[J]. Journal of Magnesium and Alloys, 2021, 9(4): 1304 − 1328. doi: 10.1016/j.jma.2020.11.018
[4] 赵军军, 张平, 王卫欣, 等. 7A52铝合金搅拌摩擦焊的焊缝成形[J]. 焊接学报, 2005, 26(5): 61 − 64. doi: 10.3321/j.issn:0253-360X.2005.05.016 Zhao Junjun, Zhang Ping, Wang Weixin, et al. Weldbead shaping of friction stir welded 7A52 aluminum Alloy[J]. Transactions of the China Welding Institution, 2005, 26(5): 61 − 64. doi: 10.3321/j.issn:0253-360X.2005.05.016
[5] Su H, Wu C S. Numerical simulation for the optimization of polygonal pin profiles in friction stir welding of aluminum[J]. Acta Metallurgica Sinica(English Letters), 2021, 34(8): 1065 − 1078. doi: 10.1007/s40195-021-01198-1
[6] Chen D G, Liu J H, Ma Z H, et al. Microstructure and properties of welding joints of 7A52 aluminum alloy by the friction stir welding[J]. Applied Mechanics and Materials, 2014, 2948: 292 − 296.
[7] 刘红伟, 周琦, 朱军, 等. 7A52铝合金厚板搅拌摩擦焊接头性能研究[J]. 兵器材料科学与工程, 2006, 29(3): 57 − 60. doi: 10.3969/j.issn.1004-244X.2006.03.016 Liu Hongwei, Zhou Qi, Zhu Jun, et al. Research on joint properties of 7A52 aluminum alloy thick plate by friction stir welding[J]. Ordnance Material Science and Engineering, 2006, 29(3): 57 − 60. doi: 10.3969/j.issn.1004-244X.2006.03.016
[8] 周鹏展, 钟掘, 贺地求. 7A52铝合金厚板搅拌摩擦焊[J]. 中国有色金属学报, 2006, 16(6): 964 − 969. doi: 10.3321/j.issn:1004-0609.2006.06.006 Zhou Pengzhan, Zhong Jue, He Diqiu. Friction-stir welding on thick plate of 7A52 aluminum alloy[J]. The Chinese Journal of Nonferrous Metals, 2006, 16(6): 964 − 969. doi: 10.3321/j.issn:1004-0609.2006.06.006
[9] Sivabalan S, Sridhar R, Parthiban A, et al. Experimental investigations of mechanical behavior of friction stir welding on aluminium alloy 6063[J]. Materials Today:Proceedings, 2021, 37(2): 1678 − 1684.
[10] 郝利新,贾瑞灵,张慧霞,等. 微弧氧化膜对7A52铝合金搅拌摩擦焊接头腐蚀不均匀性的影响[J]. 焊接学报, 2019, 40(3): 145 − 150. Hao Lixin, Jia Ruiling, Zhang Huixia, et al. Influence of micro-arc oxidation film on corrosion of inhomogeneity of 7A52 aluminum alloy friction stir welding joint[J]. Transactions of the China Welding Institution, 2019, 40(3): 145 − 150.
[11] Yuvaraj K P, Ashoka V P, Boopathiraja K P. Optimization of process parameters on friction stir welding of AA7075-T651 and AA6061 joint using response surface methodology[J]. Materials Research Express, 2019, 6(9): 6558 − 6578.
[12] Thanatkij S, Rapeepan P, Kanchana S, et al. Combined response surface method and modified differential evolution for parameter optimization of friction stir welding[J]. Processes, 2020, 8: 1080 − 1102. doi: 10.3390/pr8091080
[13] Wasif S, Salman H, Ahmad W, et al. Predicting the tensile strength, impact toughness, and hardness of friction stir-welded AA6061-T6using response surface methodology[J]. International Journal of Advanced Manufacturing Technology, 2016, 87: 1765 − 1781. doi: 10.1007/s00170-016-8565-9
[14] Ramachandran K K, Murugan N, Kumar S S. Performance analysis of dissimilar friction stir welded aluminium alloy AA5052 and HSLA steel butt joints using response surface method[J]. International Journal of Advanced Manufacturing Technology, 2016, 86: 2373 − 2392. doi: 10.1007/s00170-016-8337-6
-
期刊类型引用(8)
1. 范文学,宋飞. 互补模拟法对FSW焊接工艺参数优化的可行性分析. 当代化工研究. 2024(17): 173-175 . 百度学术
2. 陈晓,闫成旗,代俊辉. 基于IRSM的7075铝合金焊接接头抗拉强度可靠性优化设计. 焊接. 2024(10): 32-41 . 百度学术
3. 慕铠丞,王进,王勇强,李宝阁,张会. 基于响应面法的5052铝合金搅拌摩擦焊接工艺参数研究. 精密成形工程. 2023(02): 37-43 . 百度学术
4. 张开铭,王克鲁,鲁世强,柳木桐,钟平,田野. 基于物理参数和响应面法的S280超高强度不锈钢本构模型. 塑性工程学报. 2023(05): 105-115 . 百度学术
5. 王群,余洋,钱志强. 基于响应面法的HR-2抗氢钢电子束插接焊工艺参数优化. 焊接学报. 2023(04): 50-57+132 . 本站查看
6. 陈兴惠,张洪申. 基于主成分及灰色关联度分析的5083铝合金FSW接头工艺参数优化. 焊接学报. 2023(05): 62-69+132-133 . 本站查看
7. 王龙权,周海涛. 7xxx高强铝合金搅拌摩擦焊研究进展. 焊接. 2023(10): 47-54 . 百度学术
8. 李朔晗,李充,郝文波,姚肖洁,韩美玲,王依敬,邢松龄. 6082-T6铝合金三铣平面螺纹锥形针FSW工艺研究及优化. 电焊机. 2022(07): 76-83 . 百度学术
其他类型引用(3)