Prediction of primary carbide size in high carbon chromium bearing steel
-
摘要: 采用Thermo-Calc热力学计算软件对轴承钢凝固过程相转变规律、初生碳化物物相类型、初生碳化物析出温度进行热力学分析。基于Scheil偏析模型和Goto模型对轴承钢凝固过程中初生碳化物粒子尺寸进行预测,并通过热模拟试验和碳化物显微形貌定量分析对预测结果进行试验验证。计算结果表明,轴承钢凝固过程初生碳化物为(Fe,Cr)C型碳化物,析出温度为
1158 ℃。其析出行为主要是凝固末端残留液相中C和Cr元素偏析所引起的,当固相率为0.92时,残留液相中C元素的质量分数达到4.12%;残留液相中的Cr元素质量分数达到2.59%,初生碳化物尺寸为8.2 μm。此外,钢液中C元素的增加可以使初生碳化物析出温度提高,初生碳化物尺寸增大。但是,凝固过程中初生碳化物的尺寸主要是由冷却速率决定的,当冷却速率由0.2 ℃/s增加到3 ℃/s时,初生碳化物最大尺寸由19.25 μm减小至5.02 μm,其影响机制是冷却速率增大会使初生碳化物晶体各界面附近原子扩散速率降低,使得碳化物晶体各界面各向异性生长受阻,最终导致初生碳化物尺寸减小。Abstract: A Thermo-Calc thermodynamic calculation software was used to conduct thermodynamic analysis on the phase transformation, primary carbide phase type, and primary carbide precipitation temperature during solidification process of bearing steel. Based on both Scheil segregation model and Goto model, the size of primary carbide particles during solidification process of bearing steel was predicted, and the predicted results were experimentally verified through thermal simulation tests and quantitative analysis of carbide microscopic morphology. The results show that primary carbides during solidification process of bearing steel are (Fe, Cr) C carbides and its precipitation temperature is1158 ℃. The precipitation behavior is mainly caused by the segregation of C and Cr in the residual liquid phase at the end of solidification. When the solid phase ratio is 0.92, the mass fraction of C in the residual liquid phase reaches 4.12%; The mass fraction of Cr element in the residual liquid phase reached 2.59%, and the size of primary carbides reached 8.2 μm. In addition, the precipitation temperature of primary carbides and the size of primary carbides can be increased with the increase of solute C in molten steel. However, the size of primary carbides during solidification is mainly determined by the cooling rate. When the cooling rate is increased from 0.2 ℃/s to 3 ℃/s, the maximum size of primary carbides increases from 19.25 μm reduced to 5.02 μm. The influence mechanism is that increasing cooling rate will reduce the diffusion rate of atoms near the interfaces of primary carbide crystal, which will hinder the anisotropic growth of carbide crystal interfaces, and ultimately lead to decreasing primary carbide size.-
Key words:
- bearing steel /
- primary carbides /
- goto model /
- growth behavior /
- cooling rate /
- solute element /
- dynamics
-
0. 引言
冷轧双相钢(DP钢)作为第1代先进高强钢(Advanced High Strength Steel),具有屈服强度低、初始加工硬化率高以及良好强塑性匹配的特点,广泛应用于白车身零部件[1]。传统DP钢显微组织由铁素体(F)与马氏体(M)组成,其中F/M的比例决定了强度/塑性的匹配[2]。然而,随着强度的提升、塑性衰减较快,限制了复杂零件的成形与应用。随着汽车制造业的发展与升级,安全与节能成为汽车评价的主要指标,对汽车用钢的生产提出了更高要求。面对汽车行业新发展需求及高强双相钢的应用限制,首钢联合北京科技大学成功开发了增强成形性双相钢(DH钢),在传统双相钢两相组织基础上引入少量亚稳态残余奥氏体,在受力变形状态下,残余奥氏体发生相变诱导塑性(即TRIP效应),使材料强度与塑性得到提升[3],更适合用于复杂结构件与安全件的加工成形,作为新型先进高强钢极具市场应用前景。
目前,关于DH钢的研究主要集中在退火热处理工艺对组织性能的影响方面[4-6],关于DH钢与DP钢显微组织、力学性能及形变机制的差异报道较少。笔者选取典型牌号DH780与DP780作为研究对象,对不同应变速率下钢的力学性能及扩孔性能进行了对比分析,探究了增强成形性双相钢强韧性机制,对汽车选材与DH钢的应用推广具有重要参考意义。
1. 试验材料与方法
试验材料为某钢厂采用260 t转炉→LF+RH双精炼→板坯连铸→2250热连轧→2180酸连轧→2030连续退火工艺路径生产的厚度为1.2 mm 的DH780与DP780冷轧高强钢,其主要化学成分如表1所示。传统DP780采用低C-Mn-Si-Cr-Nb-Ti成分体系,低碳设计保证良好焊接性能;Mn为奥氏体稳定元素及发挥固溶强化作用[1];Si为铁素体强化元素,抑制碳化物的生成;Cr发挥固溶强化作用,同时提高淬透性;固溶Nb在高温奥氏体晶界偏聚,对高温奥氏体晶界具有拖曳作用,可细化形变奥氏体晶粒,进而细化相变组织;Nb、Ti具有第二相析出强化作用[7]。增强成形性双相钢DH780较DP780具有更高C含量,目的是实现两相区均热过程具有足够C原子向奥氏体富集,提升奥氏体(残余奥氏体)的稳定性;Al元素与Si元素作用相似,促进铁素体相变、C原子扩散与抑制碳化物析出,同时Al具有推迟珠光体相变的作用[8]。
表 1 试验钢主要化学成分Table 1. Main chemical compositions of experimental steels% 牌号 C Mn Si P S Cr Al Nb Ti DH780 0.17~0.19 2.1~2.3 0.4~0.5 ≤0.010 ≤0.005 0.18~0.21 0.7~0.9 0.02~0.04 - DP780 0.09~0.12 2.0~2.4 0.5~0.6 ≤0.010 ≤0.005 0.28~0.30 - 0.02~0.03 0.02~0.04 利用Empyrean型X射线衍射仪(XRD)对试验钢中残余奥氏体含量进行检测。试验采用Co靶、步宽为0.02°、扫描速率1°/min、扫描角度30°~130°、管电压40 kV、管电流200 mA。选取γ相中(200)、(220)、(311)衍射线与α相中(200)、(211)衍射线,利用五峰法对各晶面累计衍射强度进行计算得到残余奥氏体含量[9]。试验钢经抛光、4%硝酸酒精溶液侵蚀后分别在Zeiss金相显微镜与ZeissUItra55 型场发射扫描电子显微镜(SEM)下进行显微组织观察。经20%高氯酸酒精溶液电解抛光后进行电子背散射衍射(EBSD)表征。
为了研究板料冲压成形过程中的受力状态,同时考虑拉伸数据准确性,采用CMT5305型拉伸试验机分别以0.001、0.01、0.1 s−1的应变速率进行准静态拉伸试验,采用HTM5020型高速拉伸试验机分别以1、10、50、100、200、500、
1000 s−1的应变速率进行动态拉伸试验。按照ISO 26203-2-2011《金属材料.高应变率拉伸试验.第2部分:伺服液压和其他试验系统》要求,采用钼丝切割机将钢板加工成高速拉伸试样,线切割加工后用砂纸打磨去除线切割切割痕,以保证试样切割面平整光滑,试样尺寸如图1所示。扩孔试验能够反映板料凸缘翻边性能,为了对比成形过程凸缘翻边能力,按照ISO16630标准进行扩孔试验。2. 试验结果与分析
2.1 XRD对比分析
图2为DH780与DP780试样XRD衍射图谱,可知,DH780存在(200)γ、(220)γ、(311)γ特征峰,根据公式(1)[10]计算,可知DH780中残余奥氏体含量为5.1%,而DP780中几乎不存在残余奥氏体。
$$ {V}_{{\rm{A}}}=\frac{1-{V}_{{\rm{C}}}}{1+G\dfrac{{I}_{{\mathrm{M}}_{{\left(\mathrm{h}\mathrm{k}\mathrm{l}\right)}_{\mathrm{i}}}}}{{{{I}_{\mathrm{A}}}_{\left(\mathrm{h}\mathrm{k}\mathrm{l}\right)}}_{\mathrm{j}}}} $$ (1) 式中,VA为奥氏体相的体积分数,%;VC为碳化物相总量的体积分数,%;IM(hkl)i为马氏体(hkl)i晶面衍射线的累计强度;IA(hkl)i为奥氏体(hkl)i晶面衍射线的累计强度;G为奥氏体(hkl)晶面与马氏体(hkl)晶面所对应的强度因子之比。
2.2 显微组织对比分析
图3为DH780与DP780光学显微组织形貌,二者均为铁素体与马氏体,其中DH780晶粒均匀细小,马氏体组织分布更为弥散,而DP780组织中存在大块状多边形铁素体与明显的碳化物析出,带状组织较明显。由于光学显微镜无法分辨是否存在残余奥氏体与马氏体的分布形态,需要借助扫描电镜与EBSD进一步表征。在扫描电镜下发现,DH780与DP780显微组织中马氏体均为板条马氏体,其中DP780中马氏体以淬火态为主,而DH780中存在少量回火马氏体,分析与过时效过程中马氏体中碳化物析出有关,如图4、5所示。
综合SEM与EBSD结果显示,DH780中残余奥氏体呈块状、薄膜状、链状与细小粒状在铁素体基体中分布[11],其主要以3种形式分别位于铁素体界面(F/F,γⅠ)、位于铁素体或马氏体晶粒中(M&F,γⅡ)、位于铁素体和马氏体交界(F/M,γⅢ),且多数位于相界面与铁素体晶界处,如图4与图5所示。点链状或者薄膜长条状分布的残余奥氏体较稳定,一般在应力加载过程的中期才发生应力诱导马氏体相变;块状形式存在的残余奥氏体在施加应力初期优先发生TRIP效应。残余奥氏体γⅢ在F/M边界呈亮白边圈,主要因为两相区退火过程中Mn元素短程扩散致使在F/M边界形成富Mn区,富Mn区淬透性较高,局部区域Ms点下降,奥氏体稳定性提高,以残余奥氏体形式保留至室温,与岛内马氏体组织具有不同的色差效应。
2.3 力学性能及扩孔性能对比分析
2.3.1 不同应变速率下力学性能对比分析
钢铁材料变形本质是位错滑移与增殖的过程。冷轧双相钢以马氏体与铁素体组织为主,两相中位错主要以无钉扎自由态存在,在低应变速率下(含准静态),铁素体内可动位错优先开启进行滑移与增殖,故DP780与DH780应力-应变曲线均无明显屈服现象,如图6所示。应变速率增加初期,DP780和DH780的屈服强度、抗拉强度随应变速率的增加均呈现小幅度增加的趋势。在准静态条件下,当应变速率达到10 s−1时,屈服强度、抗拉强度出现较为明显的提升,之后随着应变速率增加,强度不断升高。
图7为应变速率对强塑性指标的影响。应变速率的提高改变了位错滑移与增殖机制,高应变速率下材料强塑性均明显增加。在应变速率由0.001 s−1增加至
1000 s−1过程中,DP780屈服强度由531 MPa增加到724 MPa,增加了36%,抗拉强度由876 MPa增加到1021 MPa,增加了16%;而DH780屈服强度由500 MPa增加到690 MPa,增加了38%,抗拉强度由796 MPa增加到997 MPa,增加了26%,说明DH780较DP780表现出较强的应变速率敏感性特征。由于DH780基体存在一定量的残余奥氏体,奥氏体面心立方结构(FCC)中可动滑移系多,有利于位错滑移,更为重要的是变形过程应变速率的提高为残余奥氏体转变为马氏体(即TRIP效应)提供了足够动力,使变形区域的塑性进一步提升,从而延缓裂纹的形成与扩展,提高材料的变形能力。残余奥氏体的引入使DH780具有TRIP效应的增塑机制,可有效降低位错运动阻力,较传统DP780表现出更高强塑特性。随着应变速率提升,DH780强塑积由27.06 GPa·%增加至38.83 GPa·%,其材料吸能性能显著增强。表2为试验钢应变速率0.01 s−1条件下的准静态力学性能,DP780和DH780垂直轧制方向(横向)的强度均高于轧制方向,DH780的断后延伸率均达到30%以上,整体断后延伸率较DP780高6.9~12.5个百分点,塑性指标显著优于DP780。
表 2 准静态力学性能(ε=0.01 s−1)Table 2. Quasi-static mechanical properties(ε=0.01 s−1)牌号 方向 屈服强度/MPa 抗拉强度/MPa 断后伸长率/% 强塑积/(GPa·%) DH780 0° 501 823 31.3 25.76 45° 487 804 31.6 25.41 90° 522 831 30.8 25.59 DP780 0° 511 818 19.5 18.41 45° 529 821 24.7 17.82 90° 541 843 18.3 16.27 2.3.2 扩孔性能对比分析
根据GB/T 24524-2021《金属材料 薄板和薄带 扩孔试验方法》分别对试验钢DH780和DP780进行扩孔率检测,测量3次取平均值。DP780试样扩孔率为50.32%,而DH780试样扩孔率达到74.61%,较DP780提升48.27%。扩孔试验过程中与锥头接触的板料下表面(凸缘翻边后为内侧)受压应力作用,而板料上表面(凸缘翻边后为外侧)受张应力作用。扩孔过程中板料上表面变形程度大于下表面,随着扩孔直径的增大,凸缘翻边外侧受到的张应力越大,当超过材料强度极限时出现裂纹萌生及扩展,裂纹由凸缘翻边外侧向内侧延伸,扩展路径与板料厚度方向呈45°分布。由于板带横向强度高,塑性指标低于轧制方向,裂纹多在垂直轧向的两侧出现,具体如图8所示。
宏观而言由于铁素体与马氏体硬度、弹性模量不同,受力过程两相变形能力存在显著差异,导致在铁素体/马氏体两相界面处产生应力集中,当相界面处应力超过两相结合力时开始萌生裂纹,并且随着变形过程裂纹沿铁素体/马氏体相界面处扩展。微观组织观察可知,DH780与DP780显微组织均以铁素体与弥散分布马氏体为主,组织细小均匀,扩孔变形初期,位错在铁素体中滑移,并在晶界与相界处产生塞积,随着应力提高,马氏体位错开启与增殖,在相界面处不断交织与缠结,形成较大应力集中,进而发展成为裂纹源。试验钢DH780中存在约5%残余奥氏体,奥氏体具有更多滑移系可有效减缓位错塞积,同时应力作用下残余奥氏体TRIP效应更为显著,从而延缓了应力集中与裂纹源的形成[12]。
图9为试验钢扩孔试验断口微观形貌,两者均为韧性断裂,DH780韧窝较大且深、撕裂脊明显,而DP780韧窝相对较浅、无明显撕裂脊。亦说明扩孔过程中DH780抗局部变形能力更强,具有更高的塑性与扩孔性能。
3. 结论
1) 试验钢DH780由铁素体、马氏体、残余奥氏体组成,其中残余奥氏体含量约5.1%,呈块状、薄膜状、链状与细小粒状,位于相界面与铁素体晶界处。
2) 随着应变速率的提升,材料屈服强度、抗拉强度均呈现增强趋势,DH780较DP780具有更高的应变速率敏感性特征,DH780强塑积增加至38.83 GPa·%,吸能性能显著增强。
3) DH780残余奥氏体在塑性变形过程中转变为马氏体,TRIP效应显著实现材料塑性提升;同时奥氏体具有更多滑移系可有效减缓位错塞积、延缓应力集中与裂纹源的形成,使得DH780较DP780具有更高扩孔率,扩孔率达到74.61%。
-
表 1 GCr15轴承主要化学成分
Table 1. Main chemical composition of GCr15 bearing steel
% C Si Mn P S Cr Al 0.97 0.21 0.32 0.017 0.001 1.49 0.0019 -
[1] 顾超. 高品质轴承钢疲劳寿命预测模型及夹杂物影响规律研究[D]. 北京: 北京科技大学, 2019.Gu Chao. Study on the fatigue life prediction model and the influence of inclusions on high quality bearing steel[D].Beijing: Beijing University of Science and Technology, 2019. [2] Wang Kun, Hu Feng, Zhou Wen, et al. Research status and development trend of bearing steel[J]. China Metallurgy, 2020,30(9):119−128. (王坤, 胡锋, 周雯, 等. 轴承钢研究现状及发展趋势[J]. 中国冶金, 2020,30(9):119−128.Wang Kun, Hu Feng, Zhou Wen, Wu Kaiming. Research Status and Development Trend of Bearing Steel[J]. China Metallurgy, 2020, 30 (09): 119-128 [3] Lü Haotian, Yang Liang, Chen Hao, et al. Long life design of bearing steel[J]. China Metallurgy, 2020,30(11):16−22. (吕皓天, 杨亮, 陈浩, 等. 轴承钢的长寿命化设计[J]. 中国冶金, 2020,30(11):16−22.Lv Haotian, Yang Liang, Chen Hao, et al. Long Life Design of Bearing Steel[J]. China Metallurgy, 2020, 30 (11): 16-22 [4] 关健. 航空滚动轴承用M50钢的接触疲劳损伤行为研究[D]. 哈尔滨: 哈尔滨工业大学, 2019.Guan Jian. Study on contact fatigue damage behavior of M50 steel for aerospace rolling bearings[D]. Harbin: Harbin Institute of Technology, 2019. [5] Zheng Kai, Cao Wenquan, Yu Feng, et al. Research and development status and progress of high temperature stainless carburizing bearing steel[J]. Iron and Steel, 2022,57(7):125−136. (郑凯, 曹文全, 俞峰, 等. 高温不锈渗碳轴承钢的研发现状与进展[J]. 钢铁, 2022,57(7):125−136.Zheng Kai, Cao Wenquan, Yu Feng, et al. Research and Development Status and Progress of High Temperature Stainless Carburizing Bearing Steel[J]. Iron and Steel, 2022, 57 (07): 125-136 [6] 呙如兵. 高速列车轴箱轴承材料微观结构及力学性能演化规律研究[D]. 北京: 北京交通大学, 2019.Guo Rubing. Study on microstructure and mechanical property evolution of high-speed train axle box bearing materials[D].Beijing: Beijing Jiaotong University, 2019. [7] 马泽铭. 高速列车轴承微观组织与力学性能演化规律及断裂机理研究[D]. 石家庄: 石家庄铁道大学, 2021.Ma Zeming. Study on the evolution law and fracture mechanism of microstructure and mechanical properties of high speed train bearings[D]. Shijiazhuang: Shijiazhuang Railway University, 2021. [8] 韩伟. 高速精密机床球轴承力学特性分析[D]. 武汉: 武汉科技大学, 2016.Han Wei. Mechanical characteristics analysis of ball bearings in high speed and precision machine tools[D]. Wuhan: Wuhan University of Science and Technology, 2016. [9] 安立愿, 李朋, 刘德义, 等. 352226X2-2RZ轴承渗碳内圈剥离失效分析[J]. 大连交通大学学报, 2020, 40(1): 80-84.An Liyuan, Li Peng, Liu Deyi, et al. Failure analysis on peeling of carburized inner ring of 352226X2-2RZ bearing[J]. Journal of Dalian Jiaotong University, 2020, 40 (1): 80-84. [10] Fu Hanwei, Cui Yinan, Zhang Chi, et al. Progress in research on rolling contact fatigue of bearing steel[J]. China Metallurgy, 2020,30(9):11−23. (付悍巍, 崔一南, 张弛, 等. 轴承钢滚动接触疲劳研究进展[J]. 中国冶金, 2020,30(9):11−23.Fu Humiwei, Cui Yinan, Zhang Chi, et al. Progress in research on rolling contact fatigue of bearing steel[J]. China Metallurgy, 2020, 30 (09): 11-23 [11] Liu Ye, Yin Qing, Li Feng, et al. Quality evaluation of ultra long fatigue life bearing steel[J]. China Metallurgy, 2020,30(9):37−40. (刘烨, 尹青, 李锋, 等. 超长疲劳寿命轴承钢的质量评价[J]. 中国冶金, 2020,30(9):37−40.Liu Ye, Yin Qing, Li Feng, et al. Quality Evaluation of Ultra Long Fatigue Life Bearing Steel[J]. China Metallurgical Journal, 2020, 30 (09): 37-40 [12] Khan F A. The effect of soaking on segregation and primary-carbide dissolution in ingot-cast bearing steel[J]. Metal, 2018,8(10):1−18. [13] Li G, Ke L, Peng W, et al. Effects of natural aging and variable loading on very high cycle fatigue behavior of a bearing steel GCr15[J]. Theoretical and Applied Fracture Mechanics, 2022,(119):103360. [14] Zhang J W, Shiozawa K, Lu L T, et al. Fatigue fracture behavior of bearing steel GCr15 in very high cycle regime[J]. Advanced Materials Research, 2008,44-46:119−126. [15] Yu Feng, Chen Xingpin, Xu Haifeng, et al. Current status of metallurgical quality and fatigue performance of rolling bearing steels and development direction of high-end bearing steels[J]. Journal of Metals, 2020,56(4):513−522. (俞峰, 陈兴品, 徐海峰, 等. 滚动轴承钢冶金质量与疲劳性能现状及高端轴承钢发展方向[J]. 金属学报, 2020,56(4):513−522.Yu Feng, Chen Xingpin, Xu Haifeng, et al. Current status of metallurgical quality and fatigue performance of rolling bearing steels and development direction of high-end bearing steels[J]. Journal of Metals, 2020, 56 (04): 513-522 [16] Zhu Zuchang, Yang Yitao. Research progress in the first, second, and third generation bearing steels and their heat treatment technologies (III)[J]. Heat Treatment Technology and Equipment, 2019,40(2):67−74. (朱祖昌, 杨弋涛. 第一、二、三代轴承钢及其热处理技术的研究进展(三)[J]. 热处理技术与装备, 2019,40(2):67−74.Zhu Zuchang, Yang Yitao. Research Progress in the First, Second, and Third Generation Bearing Steels and Their Heat Treatment Technologies (III)[J]. Heat Treatment Technology and Equipment, 2019, 40 (2): 67-74 [17] Xiao Maoguo, Li Donghui, Lü Xinyang, et al. Effect of heat treatment on the microstructure and properties of high CrCoMo bearing steel[J]. Journal of Material Heat Treatment, 2018,39(8):75−81. (肖茂果, 李东辉, 吕新杨, 等. 热处理对高Cr-Co-Mo轴承钢组织与性能的影响[J]. 材料热处理学报, 2018,39(8):75−81.Xiao Maoguo, Li Donghui, Lv Xinyang, et al. Effect of heat treatment on the microstructure and properties of high Cr Co Mo bearing steel[J]. Journal of Material Heat Treatment, 2018, 39 (08): 75-81 [18] 王升千. GCr15轴承钢低倍检验孔洞的形成机理及控制研究[D]. 北京: 北京科技大学, 2016.Wang Shengqian. Study on the formation mechanism and control of macroscopic inspection holes in GCr15 bearing steel[D]. Beijing: Beijing University of Science and Technology, 2016. [19] 库尔兹, 费希尔, 李建国, 等. 凝固原理[M]. 北京: 高等教育出版社, 2010.Kurtz, Fisher, Li Jianguo, et al. Fundamentals of solidification[M]. Beijing: Higher Education Press, 2010. [20] Goto H, Miyazawa K, Honma H. Effect of the primary oxide on the behavior of the oxide precipitating during solidification of steel[J]. ISIJ International, 1996,36(5):537−542. doi: 10.2355/isijinternational.36.537 [21] 熊辉辉. 钢中碳化物析出及其界面行为的第一性原理研究[D]. 上海: 上海大学, 2019.Xiong Huihui. First-principles study on carbide precipitation and interfacial behavior in steel[D]. Shanghai: Shanghai University, 2019. [22] Du N, Liu H, Cao Y, et al. Formation mechanism of MC and M2C primary carbides in as-cast M50 bearing steel[J]. Materials Characterization, 2021,174:111011-111024. doi: 10.1016/j.matchar.2021.111011 [23] Li Shanshan, Chen Yun, Gong Tongzhao, et al. Effect of cooling rate on solidification and precipitation mechanism of liquid precipitated carbide in high carbon chromium bearing steel[J]. Journal of Metals, 2022,58(8):1024−1034. (李闪闪, 陈云, 巩桐兆, 等. 冷速对高碳铬轴承钢液析碳化物凝固析出机制的影响[J]. 金属学报, 2022,58(8):1024−1034.Li Shanshan, Chen Yun, Gong Tongzhao, et al. Effect of cooling rate on solidification and precipitation mechanism of liquid precipitated carbide in high carbon chromium bearing steel[J]. Journal of Metals, 2022, 58 (08): 1024-1034 [24] Song W, Choi P P, Inden G, et al. On the spheroidized carbide dissolution and elemental partitioning in high carbon bearing steel 100Cr6[J]. Metallurgical and Materials Transactions A, 2014,45A(2):595−606. doi: 10.1007/s11661-013-2048-5 期刊类型引用(0)
其他类型引用(1)
-