一维高导热C/C复合材料的制备研究
林剑锋1, 袁观明1, 李轩科1, 董志军1, 张江1, 张中伟2, 王俊山2
1. 武汉科技大学 化学工程与技术学院, 武汉430081
2. 航天材料及工艺研究所, 北京100076
李轩科, 教授. E-mail:xkli8524@sina.com

林剑锋(1986-), 男, 硕士研究生. E-mail:linjianfeng20066@126.com

摘要

以三种沥青作为基体前驱体, 实验室自制的AR中间相沥青基纤维为增强体, 通过500℃热压成型, 随后经炭化和石墨化处理制备出一维炭/炭(C/C)复合材料。研究了前驱体沥青种类和热处理温度对复合材料导热性能的影响, 并采用扫描电子显微镜和偏光显微镜对其石墨化样品的形貌和微观结构进行表征。结果表明; C/C复合材料在沿纤维轴向的室温热扩散系数和导热率均随热处理温度的升高而逐渐增大; 由AR沥青作为基体前驱体所制备的C/C复合材料具有更加明显的沿纤维轴向取向的石墨层状结构以及最好的导热性能, 其3000℃石墨化样品沿纤维轴向的室温热扩散系数和导热率分别达到594.5 mm2/s和734.4 W/(m·K)。

关键词: 一维高导热C/C复合材料; 热扩散系数; 导热率
中图分类号:TB332   文献标志码:A    文章编号:1000-324X(2013)12-1338-07
Preparation of 1D C/C Composites with High Thermal Conductivity
LIN Jian-Feng1, YUAN Guan-Ming1, LI Xuan-Ke1, DONG Zhi-Jun1, ZHANG Jiang1, ZHANG Zhong-Wei2, WANG Jun-Shan2
1. School of Chemical Engineering and Technology, Wuhan University of Science and Technology, Wuhan 430081, China
2. Aerospace Research Institute of Materials and Processing Technology, Beijing 100076, China
Abstract

Three kinds of pitches were used as matrix precursors, the self-fabricated mesophase (AR) pitch-based carbon fibers were used as the reinforcing agent, one dimension (1D) carbon/carbon (C/C) composites were prepared by a hot-pressing molding method at about 500℃ and subsequent carbonization as well as graphitization treatments. The influence of precursor pitches and heat-treatment temperatures (HTTs) on the thermal conductive properties of the C/C composites were investigated. The morphology and microstructure of graphitized C/C composites were characterized by scanning electron microscope and polarized-light microscope. The results show that both of the thermal diffusivities and the thermal conductivities of C/C composites along the longitudinal direction of fibers gradually increase with HTT. The C/C composite prepared by using AR pitch as a matrix precursor possesses obvious graphite-like layer structure, exhibiting the highest thermal conductivity. The room-temperature thermal conductivity and thermal diffusivity of the C/C composite treated at 3000℃ along the longitudinal direction of fibers are 734.4 W/(m·K) and 594.5 mm2/s, respectively.

Keyword: 1D C/C composite with high thermal conductivity; thermal diffusivity; thermal conductivity

随着大型计算机和笔记本电脑等电器装置性能的不断提升以及电子元件集成度的提高, 其单位容积电子器件的发热量迅速增大, 使系统产生的热量剧增。为了保证设备的正常和稳定运行, 需要及时散发产生的热量[ 1, 2, 3]。炭/炭(C/C)复合材料因具有高导热率、低密度、比金属低得多的热膨胀系数和良好的高温机械性能等特点, 成为近年来最具发展前景的一类导热材料。该类材料中, 以中间相沥青基炭纤维或气相生长炭纤维为功能纤维的高导热C/C复合材料的研究较多[ 4, 5, 6, 7, 8]。石墨化炭纤维或气相生长炭纤维的碳层面沿纤维长度方向晶体高度取向, 这是C/C复合材料沿纤维长度方向呈现出高导热率的主要原因[ 9, 10, 11, 12, 13, 14, 15, 16, 17]。虽然C/C复合材料的高导热率主要源于功能纤维特殊的结构和取向[ 18, 19, 20], 但是基体炭的引入无疑会对C/C复合材料的导热性能产生一定的影响。本研究拟以三种不同性能沥青为基体前驱体和实验室自制的中间相沥青基纤维为增强体制备了一维高导热C/C复合材料, 研究前驱体沥青种类以及热处理温度对复合材料导热性能的影响。

1 实验部分
1.1 原料

纤维是由本实验室自制的500℃热处理AR中间相沥青基圆形截面纤维, 直径约为 35 μm。三种前驱体沥青分别为日本三菱气体化学公司生产的AR光学各向异性中间相沥青、东岛碳素生产的SC光学各向同性沥青和济宁碳素生产的MP萘基光学各向异性中间相沥青, 其物理性能如表1所示。

表1 不同前驱体沥青的基本性能和组成 Table 1 Properties and compositions of different precursor pitches
1.2 一维高导热C/C复合材料的制备

以圆形截面纤维为C/C复合材料的增强体, 分别选用AR中间相沥青、SC各向同性沥青和MP萘基中间相沥青作为基体前驱体。将前驱体沥青细粉均匀涂覆在纤维表面, 然后将均匀涂覆有沥青的纤维单向平铺于耐热不锈钢模具内, 再将模具放入热压成型机中, 快速(10~30 ℃/min)升温至500℃进行热模压成型, 恒压保温5 h得到一维C/C复合材料坯体, 其尺寸约为80 mm×40 mm×25 mm,复合材料中炭纤维长度方向与块体长度方向一致(图1)。对所得坯体样品进行高温炭化(1000℃)和石墨化(2000~3000℃)处理。为简便起见以C/CAR、C/CSC和C/CMP分别表示AR、SC和MP沥青为基体前驱体制得的一维C/C复合材料。

图1 一维C/C复合材料的块体示意图Fig. 1 Schematic diagram of 1 dimension C/C composite prepared by a hot-pressing molding method

1.3 表征及性能测试

样品残炭率为不同高温下处理后样品质量与坯体质量(500℃热压成的C/C复合材料的质量)的百分比。样品体积收缩率为(1-高温处理后样品体积/坯体体积)×100%。

C/C复合材料的导热率采用间接法测定。将样品切成10 mm×10 mm×4 mm的块体, 用激光闪烁法测定材料的热扩散系数, 材料导热率根据导热率的计算公式 λ= α× ρ× Cp计算得到。式中 α为热扩散系数, ρ为体积密度, Cp为比热容。采用德国NETZSCH Instruments公司的LFA457Nanoflash TM激光热导仪测量材料的热扩散系数。

C/C复合材料的光学织构和微观结构及形貌采用Carl Zeiss AX10型偏光显微镜(PLM)和VEGE 3 SBH-EasyProbe型扫描电子显微镜(SEM)观察; 石墨化的纤维和前驱体沥青采用Philips χ′PERT MPD PRO 型转靶X射线衍射仪(XRD)测试, 以Cu Kα(发射波长为0.154056 nm)为辐射源, 管电压为40 kV, 管电流为30 mA。

石墨的晶面层间距( d002)由布拉格公式计算;

式中 λ=0.154056 nm, θ为(002)衍射角。

平均微晶尺寸 La Lc用谢乐公式计算:

其中, β为半峰宽, K为校正因子。(002)和(100)衍射峰分别用于计算 Lc La值, 计算 Lc La K值分别取1和1.84。

石墨化度( g)是根据Franklin模式, 用Mering和Maire公式[ 21]计算:

式中: 0.3440 nm表示完全未石墨化乱层结构炭的层间距; 0.3354 nm表示理想石墨晶体的层间距。

2 结果与讨论
2.1 中间相沥青基炭纤维的晶体结构与形貌

图2为石墨化炭纤维的SEM照片和纤维不同面(纤维轴向面与径向截面)的XRD图谱。由图2(a)可以观察到非常明显的石墨层片沿纤维径向的辐射状取向。且纤维沿纤维轴向出现了劈裂, 劈裂角度较大。纤维轴向面(如图2(b)所示)的XRD图谱出现了很强的炭(002)晶面衍射峰和较弱的(004)晶面衍射峰, 未观察到其它衍射峰, 这说明石墨微晶发育良好, 石墨尺寸较大且具有较好的晶体取向。从图2(c)径向截面的XRD图谱可发现较强的炭(100)和(110)晶面衍射峰, 而(002)晶面的衍射峰强度显著减弱,该结果进一步表明石墨化炭纤维微晶尺寸大且具有优异的晶体取向, 表现出高度的各向异性。石墨化炭纤维的形貌和晶体取向结构及各向异性会对其所制备的C/C复合材料导热性能有较大影响, 并使C/C复合材料的导热性能产生各向异性。

图2 3000℃石墨化炭纤维的SEM(a)以及沿纤维轴向面(b)和径向截面(c)的XRD图谱Fig. 2 SEM image of 3000℃ graphitized fibers (a) and XRD patterns of the plane along the fiber axis (b), or radial cross section (c)

2.2 高温热处理对一维C/C复合材料物理性能的影响

图3为三种C/C复合材料残炭率、体积收缩及体积密度随热处理温度的变化。由图3(a)可知复合材料具有较高的残炭率, 失重主要发生在1000℃热处理前, 这与炭化过程中前驱体沥青的脱氢缩聚脱除碳氢小分子以及低温处理炭纤维中氧、氢及碳氢小分子的脱除有关; 热处理温度进一步升高, 前驱体沥青焦和炭纤维中非碳原子会进一步脱除, 此阶段C/C复合材料内部石墨晶体会不断发育和取向生长, 因此材料略有失重但变化不大。由图3(b)可知, 三种沥青制备的C/C复合材料在炭化和石墨化过程中体积收缩率较为一致; 随着热处理温度升高, C/C复合材料体积均不断收缩, 这可能与坯体中炭纤维含量较高有关。1000℃以前体积收缩率为32%左右, 主要是沥青脱氢缩聚形成前驱体沥青焦和低温处理炭纤维在高温下的体积收缩。进一步热处理到3000℃,C/C复合材料体积继续收缩约8%,这与其石墨晶体发育和有序生长带来的体积收缩有关。复合材料较高的残炭率和热处理过程中较大的收缩率将导致复合材料密度的较快增加。图3(c)中不同前驱体沥青复合C/CAR、C/CSC和C/CMP的密度随热处理温度升高而增加也证实了此规律。三种复合材料经3000℃处理后密度在1.70~1.76 g/cm3之间, 相差不大, 这可能与三种前驱体沥青的残炭率较为接近(表1)以及复合材料中低温处理炭纤维体积比较高有关。

图3 不同基体前驱体所得C/C复合材料物理性能随热处理温度的变化Fig.3 Physical properties of C/C composites prepared from different matrix precursors varying with heat treatment temperatures(a) Carbon yield; (b) Volume shrinkage; (c) Bulk density

2.3 一维C/C复合材料的微观结构和形貌分析

图4a为所制备典型的3000℃石墨化C/C复合材料石墨化样品主压面和截面的SEM照片, 从图4(a1、a2)可以观察到C/CAR复合材料块体主压面上平铺的圆形纤维, 圆形纤维沿其C/C复合材料长度方向整齐和均匀排列, 而且纤维较少受到损伤。从其局部放大照片图4(a2)可以看到沿纤维表面似乎覆盖有一层炭且沿纤维表面呈现取向, 这可能是热压以及随后炭化及石墨化处理使得填充在纤维周边的AR沥青沿炭纤维表面炭化和石墨化而形成的。 图4(a3)为C/CAR复合材料横截面的局部SEM照片, 显示出单根纤维的径向截面形貌, 可以观察到石墨层片沿纤维径向的辐射状取向且出现劈裂, 石墨层片的劈裂可能与高温下炭纤维石墨晶体沿径向生长和发育以及收缩有关, 并与图2所示纤维径向的形貌一致。炭纤维的这种结构取向及其在复合材料中的排列方式可能会对C/C复合材料导热性能有较大影响并出现导热性能的各向异性。

图4 3000℃石墨化处理C/CAR复合材料主压面(a1和a2)和横截面(a3)以及不同基体前驱体所得C/C复合材料局部(b1、b2: C/CAR, c1、c2: C/CSC以及d1、d2: C/CMP)的SEM照片Fig. 4 SEM images of the main surface (a1 and a2) and transverse section (a3) of the C/CAR composite and the local sections of the C/C composites (b1-b2: C/CAR, c1-c2: C/CSC and d1-d2: C/CMP) prepared from different matrix precursors after graphitization at 3000℃

图4(b1、b2)、(c1、c2)和(d1、d2)分别对应AR沥青、SC沥青和MP沥青为基体前驱体所制得的3000℃石墨化C/C复合材料局部放大的SEM照片。从图中可以清楚地观察到三种前驱体沥青所得C/C复合材料中纤维表面都附有石墨颗粒形成的层状覆盖物。图4(b1、b2)和图4(d1、d2)所对应的C/CAR和C/CMP中圆形纤维与石墨颗粒形成的层状物结合良好, 而且其石墨颗粒所形成的覆盖物表面光滑, 呈现片层状结构。这可能与这两种前驱体沥青为中间相沥青有关, 中间相沥青高度各向异性平面状液晶分子在热压力的作用下可能沿高度取向的炭纤维外表面取向从而形成片层结构, 这种取向结构炭沿纤维长度方向较为连续, 有利于复合材料沿此方向导热率的提高。由图4(c1、c2)可观察到C/CSC中纤维表面的石墨呈现颗粒状且颗粒较多连续性较差, 所形成的覆盖物比较粗糙, 与圆形纤维结合较差, 这可能与C/CSC复合材料所用前驱体沥青为各向同性氧化沥青, 所形成的炭较难石墨化形成取向结构有关。这种因基体前驱体的不同而在纤维表面形成不同形态和结构的覆盖物可能会对C/C复合材料导热性能有较大影响。表2给出了这三种前驱体沥青分别3000℃石墨化后的微晶参数。 从表中可以看出, 三种沥青衍生石墨颗粒的 d002值较为接近, 其微晶大小和石墨化度各不相同, AR和MP沥青衍生石墨的微晶尺寸( La Lc)较大, 其石墨化度较高, SC沥青衍生石墨的微晶尺寸较小。这三种沥青在复合材料中衍生的石墨颗粒的发育结晶程度也会对复合材料的导热率产生影响。对于石墨材料而言, 其室温下的导热率主要由声子的平均自由程 L的大小来决定, 石墨微晶排列的规整程度和发育程度越好, 平均微晶尺寸越大, 声子的平均自由程就越大, 相应材料的导热性能越好[ 22]

表2 不同前驱体沥青3000℃石墨化样品的微晶参数 Table 2 Microcrystalline parameters of different precursor pitches after graphitization at 3000℃

图5为不同基体前驱体所得C/C复合材料3000℃石墨化样品横截面的偏光照片。从图5(a~c)可以观察到圆形截面炭纤维较为紧密的沿径向排列, 还可以明显看到较多孔洞的存在, 显示这种纤维堆积排列形成的空穴并未被前驱体沥青炭完全填充。这些小孔洞的出现与炭纤维以及沥青炭在石墨化过程中收缩开裂有关。从图5中还可看到复合材料的孔隙率C/CSC< C/CAR< C/CMP, 表明复合材料的体积密度的变化可能为C/CSC> C/CAR> C/CMP。此结果与图3(c)的结果一致, 这可能与SC沥青残炭率相对较高(表1)有关。

图5 不同基体前驱体所得C/C复合材料3000℃石墨化样品的偏光照片Fig. 5 Polarized-light microscope photographs of the transversal section of the C/C composites prepared from different matrix precursors graphitized at 3000℃

2.4 一维C/C复合材料的导热性能

图6为不同基体前驱体所得C/C复合材料沿纤维轴向的室温热扩散系数(a)和导热率(b)随热处理温度的变化。随热处理温度的升高, C/C复合材料的热扩散系数和导热率均逐渐升高。热处理温度较低时, 热扩散系数和导热率均较小, 这主要是因为C/C复合材料内部的晶体发育还不完整, 缺陷、孔隙和杂质等还较多, 声子的平均自由程较短, 因而C/C复合材料晶体热扩散系数和导热率较低。而随着热处理温度的升高, C/C复合材料的石墨化程度升高, 石墨微晶尺寸增大、结构渐趋完整, 晶体的缺陷减少, 晶体的不完整性降低, 这些变化都将导致声子的平均自由程增大, 自由电子数增多, 声子运动的平均速度增大, 声子导热与电子导热都增强, 从而导致热扩散系数和导热率逐渐升高[ 23, 24]

图6 不同基体前驱体所得C/C复合材料沿纤维轴向的室温热扩散系数(a)和导热率(b)与热处理温度的关系Fig.6 Thermal diffusivities (a) and thermal conductivities (b) of C/C composites prepared from different matrix precursors along longitudinal direction of carbon fiber varying with heat-treatment temperatures

图6可知三种3000℃石墨化处理C/C复合材料沿纤维轴向都具有高于670 W/(m·K)的室温导热率, 其中复合材料C/CAR样品的热扩散系数和导热率分别高达594.5 mm2/s和734.4 W/(m·K)。而C/CSC虽然具有最高的密度, 但其热扩散系数和导热率却最低, 这与各向同性氧化沥青炭较难石墨化形成取向结构及内部微晶尺寸较小有关。而C/CAR和C/CMP的导热率较高与中间相AR沥青和MP沥青炭所形成的石墨晶体尺寸较大、沿纤维取向更好并与纤维结合更好有关。与MP沥青相比, AR沥青所形成的石墨沿纤维轴向取向更好, 从而使C/CAR比C/CMP具有更好的导热性能, 这可能与低温处理炭纤维是由AR沥青纺丝制得, 因而与AR沥青具有较好的相容性密切相关。上述结果与图4分析的结果是一致的。这说明前驱体沥青的性质对所得C/C复合材料的传导性能具有一定的影响。

表3不同基体前驱体所得C/C复合材料3000℃石墨化样品不同方向的热扩散系数和导热率可以看出, C/C复合材料沿纤维轴向和纤维径向的热扩散系数和导热率明显不同, 沿纤维轴向测得的热扩散系数和导热率远远高于沿纤维径向的热扩散系数和导热率, 具有显著的各向异性导热性能, 这种导热性能的差异是由复合材料中取向排列炭纤维的各向异性结构决定的。由于炭纤维石墨晶体沿纤维轴取向程度以及石墨晶体尺寸是决定其导热的主要因素, 因此石墨化炭纤维的石墨晶体碳层面沿纤维长度方向高度取向, 使得其复合材料沿纤维长度方向也呈现出高的导热率, 这与图4(a)观察结果一致。虽然C/C复合材料因前驱体沥青的不同导致沿纤维轴向热扩散系数和导热率存在较大差异, 但在沿纤维径向上看不到类似轴向的这种区别, 从表3数据可知, 不同基体前驱体所得C/C复合材料在径向上的热扩散系数和导热率相差不大。此结果表明即使辐射状炭纤维在纤维径向晶体取向较好也具有较高的导热率, 但由于纤维与纤维之间界面和空隙较多, 从而使得复合材料在沿纤维径向的热扩散系数和导热率要远远低于沿纤维轴向的热扩散系数和导热率。

表3 不同基体前驱体所得C/C复合材料3000℃石墨化样品不同方向的热扩散系数和导热率 Table 3 Thermal diffusivities and thermal conductivities of the C/C composites made from different matrix precursors along different testing planes after graphitization at 3000℃
3 结论

1) 以三种沥青为基体前驱体, 高度晶体取向的中间相沥青炭纤维为增强体, 通过一步热压成型和随后的炭化和3000℃石墨化处理制得了导热性能高度各向异性的一维高导热C/C复合材料, 其沿纤维长度方向导热率高达671~734 W/(m·K), 体积密度为1.70~1.76 g/cm3

2) 前驱体沥青的性质对所得C/C复合材料的导热性能具有较大的影响。光学各向异性中间相沥青前驱体所形成的炭与炭纤维结合紧密以及沿纤维长度方向取向和较大的微晶尺寸, 提高了C/C复合材料整体在沿纤维轴向上的导热性能。

致谢: 炭/炭复合材料导热率的测试得到清华大学深圳研究生院杜鸿达博士的大力支持与帮助。

参考文献
[1] Luedtke A. Thermal management materials for high-performance application. Adv. Eng. Mater. , 2004, 6(3): 142-144. [本文引用:1] [JCR: 1.608]
[2] Zweben C. High performance thermal management materials. Adv. Packaging, 2006, 15(2): 1-5. [本文引用:1] [JCR: 1.276]
[3] Gao X Q, Guo Q G, Liu L, et al. The study progress on carbon materials with high thermal conductivity. Journal of Functional Materials, 2006, 37(2): 173-177. [本文引用:1] [CJCR: 0.588]
[4] 高晓晴. 高导热炭/炭复合材料的制备及其性能和结构的研究. 北京: 中国科学院研究生院博士论文, 2005. [本文引用:1]
[5] Ohno H. High PerformancePitch Based Carbon Fibers and Their Application. Verbundwerkstoffe: 17 Symposium Verbundwerkstoffe und Werkstoffverbunde, 2009: 265-269. [本文引用:1]
[6] Nysten B, Issi J P. Composites based on thermally hyperconductivecarbon fiber. Composites, 1990, 21(4): 339-343. [本文引用:1]
[7] Ting J M, Lake M L. Vapor-grown carbon-fiber reinforced carbon composites. Carbon, 1995, 33(5): 663-667. [本文引用:1] [JCR: 5.868]
[8] Ting J M, Lake M L, Dufrfy D R. Composites based on thermally hyper-conductive vapor grown carbon fiber. Journal of Materials Research, 1995, 10(6): 1478-1484. [本文引用:1] [JCR: 1.713]
[9] Yuan G M, Li X K, Dong Z J, et al. Preparation and characterization of highly oriented ribbon shape pitch-based carbon fibers. Journal of Inorganic Materials, 2011, 26(10): 1026-1030. [本文引用:1] [JCR: 0.531] [CJCR: 1.106]
[10] Adams P M, Katzman H A, Rellick G S, et al. Characterization of high thermal conductivity carbon fibers and a self-reinforced graphite panel. Carbon, 1998, 36(3): 233-245. [本文引用:1] [JCR: 5.868]
[11] Nysten B, Issi J P, Barton R D, et al. Determination of lattice defects in carbon fibers by means of thermal-conductivity measurements. Phys. Rev. B, 1991, 44(5): 2142-2148. [本文引用:1]
[12] Heremans J, Beetz C P. Thermal conductivity and thermopower of vapor-grown graphite fibers. Physi. Rev. B, 1985, 32(4): 1981-1988. [本文引用:1]
[13] Edie D D. Pitch and mesophase fibers. Carbon fibers Filaments and Composites. Kluwer Academic Publishers, Dordrecht, Netherland s, 1990: 43-72. [本文引用:1]
[14] Edie D D, Fox N K, Barnett B C, et al. Melt-spun non-circular carbon fibers. Carbon, 1986, 24(4): 447-482. [本文引用:1] [JCR: 5.868]
[15] Endo M. Structure of mesophase pitch-based carbon fiber. Journal of Materials Science, 1988, 23(2): 598-605. [本文引用:1] [JCR: 2.163]
[16] Edie D D, Stoner E G. Effect of Microstructure and Shape on Carbon Fiber Properties. Carbon Materials and Composites. Noyes Publishing. .1992: 41-69. [本文引用:1]
[17] Blanco C, Appleyard S P, Rand B. Study of carbon fibers and carbon carbon composites by scanning thermal microscopy. Journal of Microscopy, 2002, 205(1): 21-32. [本文引用:1] [JCR: 1.633]
[18] Chen Y M, Ting J M. Ultra high thermal conductivity polymer composites. Carbon, 2002, 40(3): 359-362. [本文引用:1]
[19] Edie D D. The effect of processing on the structure and properties of carbon fibers. Carbon, 1988, 36(4): 345-362. [本文引用:1] [JCR: 5.868]
[20] Paiva M C, Bernardo C A, Nardin M. Mechanical surface and interfacial characterization of pitch and PAN-based carbon fibers. Carbon, 2000, 38(9): 1323-1337. [本文引用:1] [JCR: 5.868]
[21] Inagaki M, Shiraishi M, Nakamizo Met al. The evaluation of graphitization degree. Carbon Techniques, 1991(5): 38-43. [本文引用:1] [CJCR: 0.4415]
[22] He F, Yang Y G. Super thermal conductivity mesophase pitch-based carbon fibers. Hi-Tech Fiber & Application, 2003, 28(5): 28-31. [本文引用:1] [CJCR: 0.5254]
[23] Zhao J G, Li K Z, Li H J, et al. Research on the thermal conductivity of C/C composites. Acta Aeronautica et Astronautica Sinica, 2005, 26(3): 501-504. [本文引用:1] [CJCR: 0.636]
[24] Manocha L M, Warrier A, Manocha S, et al. Thermophysical properties of densified pitch based carbon/carbon materials-I. Unidirectional composites. Carbon, 2006, 44(3): 480-487. [本文引用:1] [JCR: 5.868]