en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
参考文献 1
倪萌,朱惠人,裘云,等. 航空发动机涡轮叶片冷却技术综述[J]. 燃气轮机技术, 2005, 18(4): 25-33.
参考文献 2
BOGARDD G, THOLEK A. Gas Turbine Film Cooling[J]. Journal of Propulsion and Power, 2006, 22(2): 249-270.
参考文献 3
王富强,张志强. 涡轮冷却叶片气膜孔加工方法研究[C]. 第16届全国特种加工学术会议,苏州:中国机械工程学会特种加工分会, 2015: 405-408.
参考文献 4
陈贵林,柳万珠. 航空发动机先进制造技术[M]. 西安:西北工业大学出版社, 2015: 97-100.
参考文献 5
何英,朱红钢,韩野. 航空发动机小孔特种加工技术[J]. 航空制造技术, 2011(4): 56-60.
参考文献 6
齐歆霞,朱海南. 航空发动机叶片气膜孔加工工艺分析[C]. 第14届全国特种加工学术会议,苏州:中国机械工程学会特种加工分会, 2011: 160-163.
参考文献 7
崔智勇,刘遵友,王轶,等. 叶片气膜孔加工再铸层控制[J]. 机械工程师, 2013(11): 208-209.
参考文献 8
KONONENKOT V, FREITAGC, SOVYKD N, et al. Influence of pulse repetition rate on percussion drilling of Ti-based alloy by picosecond laser pulses[J]. Optics and Lasers in Engineering, 2018(103): 65-70.
参考文献 9
彭国平,彭秀云. 发动机涡轮导向叶片气膜孔的加工工艺[J]. 航天制造技术, 2008(6): 26-29.
参考文献 10
ZHANGY, XUZ Y, ZHUY, et al. Machining of a film-cooling hole in a single-crystal superalloy by high-speed electrochemical discharge drilling[J]. Chinese Journal of Aeronautics, 2016, 29(2): 560-570.
参考文献 11
LIANGJ W, AIX, WENZ X, et al. Experimental investigation on low cycle fatigue of DZ125 with film cooling holes in different processes of laser drilling[J]. Engineering Failure Analysis, 2016(59): 326-333.
参考文献 12
IZQUIERDOB, PLAZAS, SANCHEZJ A, et al. Numerical prediction of heat affected layer in the EDM of aeronautical alloys[J]. Applied Surface Science, 2012(259): 780-790.
参考文献 13
MUHAMMADA, MIZRAJ, JUNAIDA A, et al. Parametric analysis of recast layer formation in wire-cut EDM of HSLA steel[J]. International Journal of Advanced Manufacturing Technology, 2016(87): 713–722.
参考文献 14
WANGC C, ZHOUH M, YANGL D, et al. Recast layer removal after electrical discharge machining via Taguchi analysis: A feasibility study[J]. Journal of Materials Processing Technology, 2009(209): 4134-4140.
参考文献 15
WENZ, PEIH, ZHANGC, et al. Analysis of surface quality of multi-film cooling holes in nickel-based single crystal superalloy[J]. Metal Science Journal, 2016, 32(18): 1845-1854.
目录 contents

    摘要

    针对目前航空发动机涡轮叶片气膜孔加工精度低和重熔层难去除的问题,提出了“电火花打孔、磨削扩孔”的新型气膜孔加工工艺,研制出小孔磨削专用微细CBN砂轮并对电火花气膜孔进行了磨削工艺实验。实验结果表明:经磨削加工后气膜孔圆度降低50.9%,孔径尺寸标准差降低90.7%,表面粗糙度降低65.9%,重熔层被全部去除,证明了航发涡轮叶片气膜孔磨削加工的可行性。

    Abstract

    For the problem that the machining precision of film cooling hole of aero-engine turbine blade is low and the recast layer is difficult to be removed. A new film cooling hole process which is mainly based on the thought of using EDM to drill and using grinding to expand is proposed. Grinding experiment was carried out on the film cooling hole machined by EDM using the special fine CBN grinding wheel. The result shows that after the grinding process, the roundness of the film hole is reduced by 50.9%, the standard deviation of the hole size is reduced by 90.7%, the surface roughness is reduced by 65.9%, and the recast layer is completely removed, proving the feasibility of the grinding approach of film cooling hole of aero-engine turbine blade

  • 0 引言

    随着航空发动机推重比的不断上升,燃烧室出口的温度也随之升高,涡轮前温度达到1 300~1 500 K,为了适应不断提高的工作温度,一方面需要不断改善材料性能,另一方面则需要采用冷却技[1]。在航空发动机的涡轮叶片、燃烧室、火焰筒等热端部件加工冷却气膜孔是降低零件表面温度的有效手[2]。涡轮叶片材料难加工,气膜孔孔径尺寸小,且多分布于曲面区域,材料和结构的变革对涡轮叶片等零件制孔技术提出无重熔层、无裂纹、无热影响区的质量要[3,4]

    重熔层厚度影响金属表层的均匀性,微裂纹等表面缺陷增加疲劳源产生的概率,加工精度影响气膜孔的散热效率,上述三项是衡量气膜孔加工质量的重要指[5]。目前航空发动机气膜孔的加工工艺主要有电液束打孔、激光打孔、电火花打孔以及上述方法的复合打孔方[6],电液束打孔具有无重熔层、无微裂纹、无加工应力等特点,但其加工效率低,尺寸精度重复性差,距离实际工程化应用还有一定距[7]。传统纳秒激光打孔效率高,加工能力较强,但其加工圆度差,精度低,气膜孔表面粗糙度差,且会产生难以去除的重熔[8]。电火花打孔效率较高,通过调整工艺参数可以获得较小的重熔层厚度,但是电火花电极磨损较快,且加工精度不[9]。国内外学者对上述特种加工方式进行了大量研究,Y.ZHANG等通过正交实验获得电液束加工气膜孔的优化参数,并对比了电液束和电火花加工气膜孔的表面质量和重熔层厚[10],J.W.LIANG等研究了气膜孔分布和激光打孔参数对DZ125低周疲劳寿命的影响规[11]。B.IZQUIERDO等通过分析电火花加工过程中的热源传递和材料去除建立了电火花工艺参数与重熔层厚度的预测模型并对预测模型进行了检[12]。A.MUHAMMAD等进行了电火花加工HSLA的正交试验,确定了脉冲开启时间和线速度是对重熔层影响较显著的参[13]

    本文针对气膜孔加工精度低和重熔层难去除的问题,在目前行业内最广泛采用的电火花打孔工艺基础上,通过自主研制的微型电镀CBN砂轮对电火花加工的气膜孔进行磨削工艺试验,并对电磨复合孔试样和电火花孔试样进行了圆度、表面形貌、粗糙度、重熔层等方面的对比分析,验证了气膜孔磨削加工的可行性,以期为后续航发真实叶片的气膜孔磨削加工提供数据支持和理论基础。

  • 1 实验

  • 1.1 气膜孔磨削机理分析

    在电火花加工过程中,极间介质被击穿后形成放电通道,形成的电场将使通道中的正负离子向正负极高速运动从而使电能转换为动能,正负离子高速撞击到电极表面又使动能转换为热能,从而使正负极达到高温。高温不仅使冷却液汽化和分解,更可将金属材料熔化甚至沸腾汽化,未被抛出的金属冷却凝固后附着在基体表面形成重熔[14],如图1所示。

    图1
                            重熔层

    图1 重熔层

    Fig.1 Recast layer

    为消除电火花打孔产生的重熔层,本文提出了“电火花打孔,磨削扩孔”的新型气膜孔加工工艺,该方法将气膜孔加工分为两步,如图2所示,首先通过电火花加工出比最终尺寸小的中心孔,然后用微细CBN砂轮进行气膜孔的精加工并得到最终尺寸。

    图2
                            气膜孔电火花磨削打孔复合工艺

    图2 气膜孔电火花磨削打孔复合工艺

    Fig.2 Compound EDM-grinding process for film cooling hole

    磨削工艺去除量稳定,在上一步工序重熔层厚度不大的前提下微细砂轮磨削可以实现重熔层的定量去除,并通过控制磨削参数和磨削层数来保证气膜孔的尺寸精度合格,其作用机理如图3所示。同时磨削过程磨削力小,在合理的磨削参数下不会产生表面缺陷,从而避免疲劳源在表面缺陷处产生,且电镀CBN砂轮磨削会产生一定的表面残余压应力,残余压应力可以和叶片服役时所受的拉伸载荷相抵消从而降低表面所受载荷,使得高温环境交变应力下疲劳源产生的可能性降低,有利于零件疲劳寿命的提高。

    图3
                            微细砂轮作用机理

    图3 微细砂轮作用机理

    Fig.3 Mechanism of fine grinding wheel

  • 1.2 实验设计

    试验所用材料为航空发动机所用的固溶强化型镍基高温合金GH3128,试件的尺寸为30 mm×30 mm×4 mm,预先在试件上通过电火花打出一定数量的微孔,电火花打孔工艺与某型号航空发动机叶片气膜孔加工工艺一致,之后选择一部分孔进行磨削加工。

    由于气膜孔的尺寸较小,普通砂轮无法满足要求,本文设计并研制出气膜孔磨削专用的微细CBN砂轮,砂轮基体采用直径为0.2 mm的硬质合金刀杆,将被CBN颗粒完全覆盖的基体浸入含有沉积金属的金属盐溶液电解槽中,通过电沉积的方式将金属结合剂沉积在金属基体上固定CBN磨粒,从而在刀杆的头部均匀镀上一层CBN砂粒,见图4所示。

    图4
                            微细CBN砂轮

    图4 微细CBN砂轮

    Fig.4 Fine CBN grinding wheel

    磨削采用三轴立式磨床;磨削参数为主轴转速n=60 000 r/min,进给速度f=40 mm/min,行宽Ap=0.2 mm,磨深Ae=0.005 mm。

  • 1.3 测试过程

    对试件表面使用王水进行腐蚀,腐蚀后在Zeiss Aviovert 200MAT 显微镜下观察孔边缘的金相组织和轮廓特征。

    将试件沿孔深方向进行线切割并采用日本电子株式会社JSM-6510扫描电子显微镜获取孔内剖面表面形貌,采用Taylor Hobson轮廓仪测试孔内剖面粗糙度。

    采用扫描电子显微镜配套的EDS电子元素分析仪对磨削后的内孔表面进行元素谱分析。

  • 2 结果与讨论

  • 2.1 圆度检测

    5为电磨复合孔和电火花孔在100倍显微镜头下的图像,电火花孔边缘凹凸不平,有附着或半剥离的材料,而电磨复合孔边缘非常光滑,无附着无杂质。

    图5
                            图像对比

    图5 图像对比

    Fig.5 Contrast of image

    孔的轮廓可视为孔的表面边缘和图像阴影的边界曲线,通过提取原始轮廓线并按照轮廓线计算各个孔的圆度,如图6所示,圆度计算公式为:

    r=(D-d)/2
    (1)
    图6
                            气膜孔圆度计算

    图6 气膜孔圆度计算

    Fig.6 Calculation of the roundness of hole

    式中,D为孔轮廓的外接圆直径,d为孔轮廓的内切圆直径,r为孔的圆度误差值。

    电火花孔和电磨复合孔的圆度检测结果见表1。电火花孔和电磨复合孔的圆度平均值分别为16.93和8.31 μm,电磨复合孔的圆度相比于电火花孔降低50.9%,电火花孔和电磨复合孔的孔径尺寸标准差分别为8.86和0.82 μm,电磨复合孔的孔径尺寸标准差相比于电火花孔降低90.7%。

    表1 电火花孔圆度检测结果

    Tab.1 Test result of EDM hole roundness

    工艺d/μmD/μmΦ/μmr/μm
    电火花打孔333.33362.86348.1014.77
    341.02368.16354.5913.57
    352.25383.13367.6915.44
    343.88395.85369.8725.99
    339.05371.43355.2416.19
    355.71386.94371.3315.62
    电磨复合打孔544.35562.86553.619.26
    544.08561.90553.008.91
    546.26558.78552.526.26
    546.63563.21554.428.79

    由对比结果可知电磨复合工艺下的气膜孔尺寸精度和圆度相比于电火花工艺都有较大的提升,通过控制磨削参数和磨削层数可以保证尺寸精度达到要求。

  • 2.2 表面形貌检测

    电火花孔和电磨复合孔的内孔剖面SEM图如图7所示。电火花孔剖面呈倒锥形,从入口处到出口处直径均匀减小,孔内壁粗糙且全部覆盖有重熔层,并且可以观测到图8所示的微裂纹,微孔等表面缺陷。而经过磨削加工后的孔身尺寸均匀,入口处尺寸稍大,有倒角特征,孔内壁光滑,纹路均匀,纹理方向与磨削进给方向相同,表面未发现有表面缺陷。

    图7
                            内孔剖面放大图像

    图7 内孔剖面放大图像

    Fig.7 Enlarged image of hole profile

    图8
                            电火花孔表面缺陷

    图8 电火花孔表面缺陷

    Fig.8 Defect of EDM drilling

    分别沿孔深方向对电火花孔和电磨复合孔内表面进行粗糙度测试,粗糙度取样长度为0.08 mm,评定长度为0.4 mm,粗糙度测试结果分别为0.85和0.29 μm,电磨复合孔的表面粗糙度相比于电火花孔降低65.9%。轮廓结果如图9所示,电磨复合孔相比于电火花孔表面波纹平滑。

    图9
                            表面轮廓对比

    图9 表面轮廓对比

    Fig.9 Contrast of surface profile

  • 2.3 重熔层检测

    电火花孔和电磨复合孔边缘腐蚀后的图像如图10所示,可以看出电火花孔的孔边缘经腐蚀后出现贝壳状的重熔层,可见厚度为20 μm,考虑到金相制样中小孔边缘凹陷效应,实际重熔层厚度会大于20 μm,而经过磨削加工后已观察不到孔边缘重熔层的痕迹。

    图10
                            孔边缘金相图

    图10 孔边缘金相图

    Fig.10 Metallographic picture of hole edge

    电火花加工过程中由于介质水的分解会产生氧气,熔化的材料与氧反应从在基体表面生成氧化层,因此电火花加工形成的重熔层中可以检测到一定的氧含[15]。对磨削后的内孔表面进行EDS分析结果见图11,磨削表面的氧含量几乎为0,说明由电火花打孔产生的重熔层经过磨削加工后已被全部去除。

    图11
                            电磨复合孔EDS分析结果

    图11 电磨复合孔EDS分析结果

    Fig.11 EDS analysis result of compound EDM-grinding hole

  • 3 结论

    通过电火花磨削复合工艺得到了高精度无重熔层的气膜孔,验证了航发涡轮叶片气膜孔磨削加工的可行性。

    (1)电火花孔边缘凹凸不平,经磨削加工后孔边缘光滑,圆度误差降低50.9%,孔径尺寸标准差降低90.7%。

    (2)电火花孔内壁粗糙,表面存在微裂纹、微孔洞等缺陷,经磨削加工后孔身尺寸均匀,内壁光滑,无表面缺陷,表面粗糙度降低65.9%。

    (3)电火花孔的重熔层厚度大于20 μm,磨削去量稳定可控,磨削加工可以将重熔层全部去除。

    后续将通过优化磨削参数和研制尺寸更小的CBN磨具来提高磨削效率和孔径可加工范围。

  • 参考文献

    • 1

      倪萌,朱惠人,裘云,等. 航空发动机涡轮叶片冷却技术综述[J]. 燃气轮机技术, 2005, 18(4): 25-33.

    • 2

      BOGARD D G, THOLE K A. Gas Turbine Film Cooling[J]. Journal of Propulsion and Power, 2006, 22(2): 249-270.

    • 3

      王富强,张志强. 涡轮冷却叶片气膜孔加工方法研究[C]. 第16届全国特种加工学术会议,苏州:中国机械工程学会特种加工分会, 2015: 405-408.

    • 4

      陈贵林,柳万珠. 航空发动机先进制造技术[M]. 西安:西北工业大学出版社, 2015: 97-100.

    • 5

      何英,朱红钢,韩野. 航空发动机小孔特种加工技术[J]. 航空制造技术, 2011(4): 56-60.

    • 6

      齐歆霞,朱海南. 航空发动机叶片气膜孔加工工艺分析[C]. 第14届全国特种加工学术会议,苏州:中国机械工程学会特种加工分会, 2011: 160-163.

    • 7

      崔智勇,刘遵友,王轶,等. 叶片气膜孔加工再铸层控制[J]. 机械工程师, 2013(11): 208-209.

    • 8

      KONONENKO T V, FREITAG C, SOVYK D N, et al. Influence of pulse repetition rate on percussion drilling of Ti-based alloy by picosecond laser pulses[J]. Optics and Lasers in Engineering, 2018(103): 65-70.

    • 9

      彭国平,彭秀云. 发动机涡轮导向叶片气膜孔的加工工艺[J]. 航天制造技术, 2008(6): 26-29.

    • 10

      ZHANG Y, XU Z Y, ZHU Y, et al. Machining of a film-cooling hole in a single-crystal superalloy by high-speed electrochemical discharge drilling[J]. Chinese Journal of Aeronautics, 2016, 29(2): 560-570.

    • 11

      LIANG J W, AI X, WEN Z X, et al. Experimental investigation on low cycle fatigue of DZ125 with film cooling holes in different processes of laser drilling[J]. Engineering Failure Analysis, 2016(59): 326-333.

    • 12

      IZQUIERDO B, PLAZA S, SANCHEZ J A, et al. Numerical prediction of heat affected layer in the EDM of aeronautical alloys[J]. Applied Surface Science, 2012(259): 780-790.

    • 13

      MUHAMMAD A, MIZRA J, JUNAID A A, et al. Parametric analysis of recast layer formation in wire-cut EDM of HSLA steel[J]. International Journal of Advanced Manufacturing Technology, 2016(87): 713–722.

    • 14

      WANG C C, ZHOU H M, YANG L D, et al. Recast layer removal after electrical discharge machining via Taguchi analysis: A feasibility study[J]. Journal of Materials Processing Technology, 2009(209): 4134-4140.

    • 15

      WEN Z, PEI H, ZHANG C, et al. Analysis of surface quality of multi-film cooling holes in nickel-based single crystal superalloy[J]. Metal Science Journal, 2016, 32(18): 1845-1854.

李潜彤

机 构:北京航空航天大学机械工程及自动化学院,北京 100191

Affiliation:School of Mechanical Engineering and Automation, Beihang University,Beijing 100191

角 色:第一作者

Role:First author

邮 箱:lqt_1996@163.com

第一作者简介:李潜彤,1996年出生,硕士,研究方向为难加工材料加工表面完整性。E-mail:lqt_1996@163.com

陈志同

机 构:北京航空航天大学机械工程及自动化学院,北京 100191

Affiliation:School of Mechanical Engineering and Automation, Beihang University,Beijing 100191

全芳

机 构:北京航空航天大学机械工程及自动化学院,北京 100191

Affiliation:School of Mechanical Engineering and Automation, Beihang University,Beijing 100191

马磊

机 构:北京航空航天大学机械工程及自动化学院,北京 100191

Affiliation:School of Mechanical Engineering and Automation, Beihang University,Beijing 100191

html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image001.jpeg
html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image003.jpeg
html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image005.jpeg
html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image007.jpeg
html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image009.jpeg
html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image011.jpeg
工艺d/μmD/μmΦ/μmr/μm
电火花打孔333.33362.86348.1014.77
341.02368.16354.5913.57
352.25383.13367.6915.44
343.88395.85369.8725.99
339.05371.43355.2416.19
355.71386.94371.3315.62
电磨复合打孔544.35562.86553.619.26
544.08561.90553.008.91
546.26558.78552.526.26
546.63563.21554.428.79
html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image013.jpeg
html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image015.jpeg
html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image017.jpeg
html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image019.jpeg
html/yhclgy/201902017/media/628517c7-2dca-4aa8-92cc-c5f96c93093c-image021.jpeg

图1 重熔层

Fig.1 Recast layer

图2 气膜孔电火花磨削打孔复合工艺

Fig.2 Compound EDM-grinding process for film cooling hole

图3 微细砂轮作用机理

Fig.3 Mechanism of fine grinding wheel

图4 微细CBN砂轮

Fig.4 Fine CBN grinding wheel

图5 图像对比

Fig.5 Contrast of image

图6 气膜孔圆度计算

Fig.6 Calculation of the roundness of hole

表1 电火花孔圆度检测结果

Tab.1 Test result of EDM hole roundness

图7 内孔剖面放大图像

Fig.7 Enlarged image of hole profile

图8 电火花孔表面缺陷

Fig.8 Defect of EDM drilling

图9 表面轮廓对比

Fig.9 Contrast of surface profile

图10 孔边缘金相图

Fig.10 Metallographic picture of hole edge

图11 电磨复合孔EDS分析结果

Fig.11 EDS analysis result of compound EDM-grinding hole

image /

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

  • 参考文献

    • 1

      倪萌,朱惠人,裘云,等. 航空发动机涡轮叶片冷却技术综述[J]. 燃气轮机技术, 2005, 18(4): 25-33.

    • 2

      BOGARD D G, THOLE K A. Gas Turbine Film Cooling[J]. Journal of Propulsion and Power, 2006, 22(2): 249-270.

    • 3

      王富强,张志强. 涡轮冷却叶片气膜孔加工方法研究[C]. 第16届全国特种加工学术会议,苏州:中国机械工程学会特种加工分会, 2015: 405-408.

    • 4

      陈贵林,柳万珠. 航空发动机先进制造技术[M]. 西安:西北工业大学出版社, 2015: 97-100.

    • 5

      何英,朱红钢,韩野. 航空发动机小孔特种加工技术[J]. 航空制造技术, 2011(4): 56-60.

    • 6

      齐歆霞,朱海南. 航空发动机叶片气膜孔加工工艺分析[C]. 第14届全国特种加工学术会议,苏州:中国机械工程学会特种加工分会, 2011: 160-163.

    • 7

      崔智勇,刘遵友,王轶,等. 叶片气膜孔加工再铸层控制[J]. 机械工程师, 2013(11): 208-209.

    • 8

      KONONENKO T V, FREITAG C, SOVYK D N, et al. Influence of pulse repetition rate on percussion drilling of Ti-based alloy by picosecond laser pulses[J]. Optics and Lasers in Engineering, 2018(103): 65-70.

    • 9

      彭国平,彭秀云. 发动机涡轮导向叶片气膜孔的加工工艺[J]. 航天制造技术, 2008(6): 26-29.

    • 10

      ZHANG Y, XU Z Y, ZHU Y, et al. Machining of a film-cooling hole in a single-crystal superalloy by high-speed electrochemical discharge drilling[J]. Chinese Journal of Aeronautics, 2016, 29(2): 560-570.

    • 11

      LIANG J W, AI X, WEN Z X, et al. Experimental investigation on low cycle fatigue of DZ125 with film cooling holes in different processes of laser drilling[J]. Engineering Failure Analysis, 2016(59): 326-333.

    • 12

      IZQUIERDO B, PLAZA S, SANCHEZ J A, et al. Numerical prediction of heat affected layer in the EDM of aeronautical alloys[J]. Applied Surface Science, 2012(259): 780-790.

    • 13

      MUHAMMAD A, MIZRA J, JUNAID A A, et al. Parametric analysis of recast layer formation in wire-cut EDM of HSLA steel[J]. International Journal of Advanced Manufacturing Technology, 2016(87): 713–722.

    • 14

      WANG C C, ZHOU H M, YANG L D, et al. Recast layer removal after electrical discharge machining via Taguchi analysis: A feasibility study[J]. Journal of Materials Processing Technology, 2009(209): 4134-4140.

    • 15

      WEN Z, PEI H, ZHANG C, et al. Analysis of surface quality of multi-film cooling holes in nickel-based single crystal superalloy[J]. Metal Science Journal, 2016, 32(18): 1845-1854.