Study on size-dependence of tensile mechanical properties and plastic deformation mechanism in CoCrFeNiMn high entropy alloy
-
摘要: 采用分子动力学模拟方法研究了不同横截面尺寸的CoCrFeNiMn高熵合金在拉伸载荷下的力学性能。模拟采用修正的嵌入原子势函数描述Co、Cr、Fe、Ni和Mn原子之间的相互作用。研究表明，CoCrFeNiMn高熵合金的拉伸性能和变形机制依赖于横截面的尺寸。峰值应力在横截面尺寸为7.00 nm时，出现了临界值，这与其中位错成核点的数目相关。随着横截面尺寸增加，塑性变形机制从非晶化主导转变为位错滑移和FCC结构向HCP结构的相变主导。该研究结果对于设计和制备高性能的高熵合金具有一定的科学价值和指导意义。Abstract: The mechanical properties of CoCrFeNiMn high entropy alloys (HEAs) with different cross section sizes were investigated by molecular dynamics simulation under tensile loading. Modified embedding atomic potential was employed to describe the interactions between Co, Cr, Fe, Ni and Mn atoms. The simulation results indicated that the tensile properties and deformation mechanism of CoCrFeNiMn HEAs depended on the cross section size. A critical value of peak stress occurred when the cross section size was 7.00 nm, which depended on the number of dislocation nucleation sites. With the increase of cross section size, the plastic deformation mechanism changed from the formation of amorphous structures to dislocation slip and transformation from FCC structure to HCP structure. The results of this study can provide a guidance for design and preparation of high entropy alloy with high performance.
-
Key words:
- high entropy alloy /
- mechanical behavior /
- size effect /
- deformation mechanism
-
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标准进行扩孔试验。
图 1 高速拉伸试样示意（单位：mm）Figure 1. Schematic diagram of the sample for high-speed stretching tests2. 试验结果与分析
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效应）提供了足够动力，使变形区域的塑性进一步提升，从而延缓裂纹的形成与扩展，提高材料的变形能力。
图 7 强塑性指标与应变速率的关系Figure 7. Relationship between strength ＆ plasticity index and strain rate残余奥氏体的引入使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 CoCrFeNiMn高熵合金模型
（a）初始构型 ；（b）经公共近邻分析得到的原子结构
Figure 1. Schematic diagram of CoCrFeNiMn high-entropy alloy model
图 2 不同横截面尺寸高熵合金的应力-应变曲线
Figure 2. Stress-strain curves of CoCrFeNiMn high-entropy alloys with different cross-section sizes
图 3 高熵合金的峰值应力和平均流动应力随横截面尺寸大小的变化关系
Figure 3. Variation of peak stress and average flow stress of CoCrFeNiMn high-entropy alloys with different cross-section sizes
图 4 横截面尺寸为3.00 nm的高熵合金在不同拉伸应变下的原子结构快照
Figure 4. Atomic configuration evolutions of the CoCrFeNiMn high-entropy alloys with cross-section of 3.00 nm at different tensile strains
图 5 横截面尺寸为7.00 nm的高熵合金在不同拉伸应变下的结构快照
Figure 5. Atomic configuration evolutions of the CoCrFeNiMn high-entropy alloys with cross-section of 7.00 nm at different tensile strains
图 6 横截面尺寸为13.00 nm的高熵合金在不同拉伸应变下的原子结构快照
Figure 6. Atomic configuration evolutions of the CoCrFeNiMn high-entropy alloys with cross-section of 13.00 nm at different strains
图 7 不同横截面尺寸高熵合金在拉伸过程中的不同结构原子比例变化
Figure 7. Atomic fraction evolutions of different structures of the CoCrFeNiMn high-entropy alloys with different cross section size
-
[1] Yeh J W, Chen S K, Lin S J, et al. Nanostructured high-entropy alloys with multiple principal elements: Novel alloy design concepts and outcomes[J]. Advanced Engineering Materials, 2004,6:299−303. doi: 10.1002/adem.200300567 [2] George E P, Raabe D, Ritchie R O. High-entropy alloys[J]. Nature Reviews Materials, 2019,4:515−534. [3] 乔珺威. 面心立方结构高熵合金[M]. 北京: 冶金工业出版社, 2021.Qiao Junwei. High entropy alloy with face-centered cubic structure [M].Beijing: Metallurgical Industry Press, 2021. [4] 王龙. 高熵合金涂层的微结构与耐腐蚀性能研究[D]. 自贡: 四川理工学院, 2015.Wang Long. Study on microstructure and corrosion resistance of high entropy alloy coatings[D]. Zigong: Sichuan University of Science and Technology, 2015. [5] 鲁一荻, 张骁勇, 侯硕. 高熵合金的发展及工业应用展望[J]. 稀有金属材料与工程, 2021, 50(1): 333-341.Lu Yidi, Zhang Xiaoyong, Hou Shuo, et al. Development and industrial application prospect of high entropy alloys[J]. Rare Metal Materials and Engineering, 2021, 50(1): 333-341. [6] 张荣, 祁文军, 张爽. AlxCoCrFeNi拉伸力学性能的分子动力学模拟[J]. 钢铁钒钛, 2022, 43(6): 173-179.Zhang Rong, Qi Wenjun, Zhang Shuang. Molecular dynamics simulation of tensile mechanical properties of AlxCoCrFeNi[J], Iron Steel Vanadium Titanium, 2022,43(6):173-179. [7] Zhao Haichao, Qiao Yulin, Liang Xiubing, et al. Research progress and prospect of lightweight high entropy alloys[J]. Rare Metal Materials and Engineering, 2020,49(4):1457−1468. (赵海朝, 乔玉林, 梁秀兵, 等. 轻质高熵合金的研究进展与展望[J]. 稀有金属材料与工程, 2020,49(4):1457−1468.Zhao Haichao, Qiao Yulin, Liang Xiubing et al. Research progress and prospect of Lightweight High Entropy Alloys [J]. Rare Metal Materials and Engineering, 2020, 49(04): 1457-1468. [8] Csikor F F, Motz C, Weygand D, et al. Dislocation avalanches, strain bursts, and the problem of plastic forming at the micrometer scale[J]. Science, 2017,318:251−254. [9] Xie K, Sirestha S, Cao Y, et al. The effect of pre-existing defects on the strength and deformation behavior of a-Fe nanopillars[J]. Acta Materialia, 2013,61:439−452. doi: 10.1016/j.actamat.2012.09.022 [10] Zhang J Y, Lei S Y, Liu Y, et al. Length scale-dependent deformation behavior of nanolayered Cu/Zr micropillars[J]. Acta Materialia, 2012,(60):1610−1622. [11] Zhang J Y, Liu G, Lei S Y, et al. Transition from homogeneous-like to shear-band deformation in nanolayered crystalline Cu/amorphous Cu–Zr micropillars: Intrinsic vs. extrinsic size effect[J]. Acta Materialia, 2012,(60):7183−7196. [12] Yu Zou, Soumyadipta Maiti, Walter Steurer, et al. Size-dependent plasticity in an Nb25Mo25Ta25W25 refractory high-entropyalloy[J]. Acta Materialia, 2014,65:85−97. doi: 10.1016/j.actamat.2013.11.049 [13] Indranil Basu, Vaclav Ocelík, Jeff Th M, et al. Size dependent plasticity and damage response in multiphase body centered cubic high entropy alloys[J]. Acta Materialia, 2018,150:104−116. doi: 10.1016/j.actamat.2018.03.015 [14] An Minrong, Li Silan, Su Mengjia. Molecular dynamics simulation of size-dependent plastic deformation mechanism of CoCrFeNiMn crystal/amorphous biphase high entropy alloy[J]. Acta Physica Sinica, 2022,(24):146−157. (安敏荣, 李思澜, 宿梦嘉. 尺寸依赖的CoCrFeNiMn晶体/非晶双相高熵合金塑性变形机制的分子动力学模拟研究[J]. 物理学报, 2022,(24):146−157. doi: 10.7498/aps.71.20221368An Minrong, Li Silan, Su Mengjia. Molecular dynamics simulation of size-dependent plastic deformation mechanism of CoCrFeNiMn crystal/amorphous biphase high entropy alloy. Acta Physica Sinica. doi: 10.7498/aps.71.20221368. [15] Xiao L L, Zheng Z Q, Guo S W, et al. Ultra-strong nanostructured CrMnFeCoNi high entropy alloys[J]. Materials and Design, 2020,194:108895. doi: 10.1016/j.matdes.2020.108895 [16] Ji W, Wu M S. Nanoscale origin of the crystalline-to-amorphous phase transformation and damage tolerance of Cantor alloys at cryogenic temperatures[J]. Acta Materialia, 2022,226:117639. doi: 10.1016/j.actamat.2022.117639 [17] Zhao S, Li Z, Zhu C, et al. Amorphization in extreme deformation of the CrMnFeCoNi high-entropy alloy[J]. Science Advance, 2021,7:eabb3108. [18] Li J, Bai Z, Max Powers. et al. Deformation mechanisms in crystalline-amorphous high-entropy composite multilayers[J]. Materials Science & Engineering A, 2022, 848: 143144. [19] Ren J, Sun J, Xiao L, et al. Size-dependent of compression yield strength and deformation mechanism in titanium single-crystal nanopillars orientated [0001] and [11-20][J]. Materials Science & Engineering A, 2014,615:22−28. [20] Shuichi Nosé. A unified formulation of the constant temperature molecular dynamics methods[J]. The Journal of Chemical Physics, 1984, 81: 511. [21] Won-Mi Choi, Yong Hee Jo, Seok Su Sohn, et al. Understanding the physical metallurgy of the CoCrFeMnNi high-entropy alloy: an atomistic simulation study[J]. NPJ Computational Materials, 2018,4:1. doi: 10.1038/s41524-017-0060-9 [22] Utt D, Lee S, Xing Y, et al. The origin of jerky dislocation motion in high-entropy alloys[J]. Nature Communications, 2022,13:4777. doi: 10.1038/s41467-022-32134-1 [23] Hayakawa S, Xu H. Temperature-dependent mechanisms of dislocation-twin boundary interactions in Ni-based equiatomic alloys[J]. Acta Materialia, 2021,211:116886. doi: 10.1016/j.actamat.2021.116886 [24] Plimpton S, Fast parallel algorithms for short-range molecular dynamics[J].Journal of Computational Physics, 1995, 117: 1-19. [25] Stukowski A. Visualization and analysis of atomistic simulation data with OVITO-the open visualization tool [J].Modelling and Simulation in Materials Science and Engineering, 2009, 18: 015012. [26] Faken D, Jónsson H. Systematic analysis of local atomic structure combined with 3D computer graphics [J]. Computational Materials Science, 1994, 2: 279-286. [27] Li J, Chen H, Feng H, et al. Microstructure evolution and deformation mechanism of amorphous/crystalline high-entropy-alloy composites[J]. Journal of Materials Science & Technology, 2020,54:14−19. [28] Shen Tianzhan, Song Haiyang, An Minrong. Molecular dynamics simulation of the effect of twin boundary on mechanical behavior of Cr26Mn20Fe20Co20Ni14 high entropy alloy[J]. Acta Physica Sinica, 2021,70(18):186201. (申天展, 宋海洋, 安敏荣. 孪晶界对Cr26Mn20Fe20Co20Ni14高熵合金力学行为影响的分子动力学模拟[J]. 物理学报, 2021,70(18):186201. doi: 10.7498/aps.70.20210324Shen Tianzhan, Song Haiyang, An MinRong. Molecular dynamics simulation of the effect of twin boundary on mechanical behavior of Cr26Mn20Fe20Co20Ni14 high entropy alloy[J]. Acta Physica Sinica, 2021, 70(18): 186201. doi: 10.7498/aps.70.20210324 [29] Qi Y M, Zhao M, Feng M L. Molecular simulation of microstructure evolution and plastic deformation of nanocrystalline CoCrFeMnNi high-entropy alloy under tension and compression[J]. Journal of Alloys and Compounds, 2021,851:156923. doi: 10.1016/j.jallcom.2020.156923 [30] Fang Q H, Chen Y, Li J, et al. Probing the phase transformation and dislocation evolution in dual-phase high-entropy alloys[J]. International Journal of Plasticity, 2019,114:161−173. doi: 10.1016/j.ijplas.2018.10.014 [31] Chowdhury P, Canadinc D, Sehitoglu H. On deformation behavior of Fe-Mn based structural alloys[J]. Materials Science and Engineering R, 2017, 122: 1-28. 期刊类型引用(0)
其他类型引用(1)
-
下载:
下载:
点击查看大图
计量
- 文章访问数: 256
- HTML全文浏览量: 59
- PDF下载量: 17
- 被引次数: 1
