Gd对ZK60镁合金组织与力学性能的影响
来源期刊:中国有色金属学报(英文版)2013年第9期
论文作者:黄正华 戚文军 徐 静
文章页码:2568 - 2576
关键词:ZK60镁合金;Gd变质;挤压;显微组织;力学性能;断口形貌
Key words:ZK60 magnesium alloy; Gd modification; extrusion; microstructure; mechanical property; fracture morphology
摘 要:分析了铸态和挤压态ZK60-xGd(x=0~4)合金的组织和相组成,测试了其拉伸力学性能。结果表明,随着Gd含量的增加,铸态组织逐渐细化,Mg-Zn-Gd新相逐渐增多,而MgZn2相逐渐减少直至消失,第二相趋于连续网状分布于晶界处;当Gd含量不超过2.98%时,铸态室温拉伸力学性能稍降低。经挤压比l=40和挤压温度 T=593 K的挤压后,组织显著细化,平均晶粒尺寸逐渐减至ZK60-2.98Gd合金的2 μm,破碎的第二相沿着挤压方向呈带状分布;挤压态的拉伸力学性能均显著提高:298和473 K时的抗拉强度分别从ZK60合金的355和120 MPa逐渐提高至ZK60-2.98Gd合金的380和164 MPa。挤压态拉伸断口呈现典型的韧性断裂特征。
Abstract: Microstructures and phase compositions of as-cast and extruded ZK60-xGd (x=0-4) alloys were investigated. Meanwhile, the tensile mechanical property was tested. With increasing the Gd content, as-cast microstructure is refined gradually. Mg-Zn-Gd new phase increases gradually, while MgZn2 phase decreases gradually to disappear. The second phase tends to distribute along grain boundary by continuous network. As-cast tensile mechanical property is reduced slightly at ambient temperature when the Gd content does not exceed 2.98%. After extrusion by extrusion ratio of 40 and extrusion temperature of 593 K, microstructure is refined further with decreasing the average grain size to 2 μm for ZK60-2.98Gd alloy. Broken second phase distributes along the extrusion direction by zonal shape. Extruded tensile mechanical property is enhanced significantly. Tensile strength values at 298 and 473 K increase gradually from 355 and 120 MPa for ZK60 alloy to 380 and 164 MPa for ZK60-2.98Gd alloy, respectively. Extruded tensile fractures exhibit a typical character of ductile fracture.
Trans. Nonferrous Met. Soc. China 23(2013) 2568-2576
Zheng-hua HUANG, Wen-jun QI, Jing XU
Guangzhou Research Institute of Non-ferrous Metals, Guangzhou 510650, China
Received 12 July 2012; accepted 15 November 2012
Abstract: Microstructures and phase compositions of as-cast and extruded ZK60-xGd (x=0-4) alloys were investigated. Meanwhile, the tensile mechanical property was tested. With increasing the Gd content, as-cast microstructure is refined gradually. Mg-Zn-Gd new phase increases gradually, while MgZn2 phase decreases gradually to disappear. The second phase tends to distribute along grain boundary by continuous network. As-cast tensile mechanical property is reduced slightly at ambient temperature when the Gd content does not exceed 2.98%. After extrusion by extrusion ratio of 40 and extrusion temperature of 593 K, microstructure is refined further with decreasing the average grain size to 2 μm for ZK60-2.98Gd alloy. Broken second phase distributes along the extrusion direction by zonal shape. Extruded tensile mechanical property is enhanced significantly. Tensile strength values at 298 and 473 K increase gradually from 355 and 120 MPa for ZK60 alloy to 380 and 164 MPa for ZK60-2.98Gd alloy, respectively. Extruded tensile fractures exhibit a typical character of ductile fracture.
Key words: ZK60 magnesium alloy; Gd modification; extrusion; microstructure; mechanical property; fracture morphology
1 Introduction
Magnesium alloys with high specific strength are used widely in automotive, communicated, electronic and aerial industries [1]. Wrought magnesium alloys exhibit higher strength, better ductility and comprehensive property; however, their plasticity is still relatively poor owing to the limit of crystal structure (hexagonal close packed, HCP) [2]. Therefore, it is necessary to develop new wrought magnesium alloys with high performance by investigating the deformation mechanism.
Rare earth (RE) with unique atomic electron and chemical property can purify alloy melt, ameliorate the microstructure and enhance the mechanical property effectively. Researchers had extensively studied the effects of Y [3-5], Ce [6], Nd [7], Yb [8] and Ho [9] on the microstructure and tensile mechanical property of ZK60 alloy, being one of the most widely used wrought magnesium alloys currently. Gd was commonly used to develop Mg-Gd-Y series heat-resistant magnesium alloys [10-13]; however, its effects on the microstructure and tensile mechanical property of ZK60 alloy had been rarely studied yet. As well as grain size, the morphology, distribution and volume fraction of second phase can also affect the tensile mechanical property of magnesium alloys significantly [7]. Therefore, the present work is focused on studying the microstructures and tensile mechanical properties of as-cast and extruded ZK60-xGd (x=0-4) alloys in detail, and then discussing the role of grain size and second phase on the strengthening and deformation mechanism.
2 Experimental
Chemical compositions of as-cast ZK60-xGd alloys measured by inductively coupled plasma analyzer (ICP, JY Ultima2) are listed in Table 1. Alloy ingots were prepared by melting pure Mg, pure Zn and Mg-30%Zr, Mg-20%Gd master alloys in electric resistance furnace under the mixed atmosphere of CO2 and SF6. When the melt temperature of pure Mg reached 1003 K, pure Zn and two master alloys were added into the melt in turn. Then, the melt was stirred twice within 1 h to ensure the compositional homogeneity. After adding refine agent (JDMJ), the melt was held at 1033 K for 30 min. When the temperature cooled to 988 K, the melt was poured into the wedge permanent mold with preheated temperature of 523 K, and then as-cast samples were obtained. Cast rods obtained under the same technology were homogenized at 673 K for 8 h in the heat-treating furnace, and then machined to a diameter of 100 mm. Extrusion rods were obtained by extruder with extrusion ratio l, extrusion temperature T and extrusion speed v of 40, 593 K and 1-2 m/min, respectively.
The as-cast and extruded specimens were etched with 4% HNO3 (volume fraction) in ethanol and a solution of 1.5 g picric acid, 25 mL ethyl alcohol, 5 mL acetic acid plus 10 mL distilled water, respectively. Microstructural observation was carried out on an optical microscope (OM, Leica DM IRM), scanning electron microscope (SEM, JEOL JXA-8100) with energy dispersive spectroscopy (EDS, OXFORD 7412) and transmission electron microscope (TEM, Tecnai G2 20). Phase analysis was carried out on an X-ray diffractometer (XRD, D/MAX-RC) with Cu Kα radiation. Phase change during the heating process was characterized by isochronally heating in a differential scanning calorimeter (DSC, Netzsch STA409) under flow of highly-purified Ar atmosphere at a rate of 20 K/min. Tensile test at the temperature ranging from 298 K to 473 K was performed on a test machine (GP-TS2000) at a rate of 2 mm/min. Tensile fractographs were observed using SEM.
Table 1 Chemical composition of ZK60-xGd alloys (mass fraction, %)
3 Results
3.1 Microstructure and mechanical property of as-cast alloy
Figures 1 and 2 show the optical graphs and SEM images of as-cast alloys, respectively. As-cast ZK60 alloy exhibits coarse microstructure and few dispersed second phases. With increasing the Gd content, as-cast microstructure is refined gradually. Meanwhile, the second phase increases gradually and tends to distribute along grain boundary by continuous network.
Fig. 1 Optical graphs of as-cast ZK60-xGd alloys
Fig. 2 SEM images of as-cast ZK60-xGd alloys
Figure 3 shows the EDS spectra of as-cast alloys and the results are listed in Table 2. For as-cast ZK60 alloy, there exist two different Mg-Zn phases containing relatively low and high Zr contents, respectively, while new Mg-Zn-Gd plus Mg-Zn and only Mg-Zn-Gd phases can be observed for as-cast ZK60-0.28Gd and ZK60-2.98Gd alloys, respectively. However, mole fraction in Mg-Zn phase is slightly different from MgZn2 phase determined by following XRD patterns owing to the affection by adjacent a-Mg matrix, and that in Mg-Zn-Gd phase changes slightly with most of the Gd content of about 5% in mole fraction. Majority of element Zn distributes along grain boundary to form MgZn2 phase for as-cast ZK60 alloy. Similar to element Zn, element Gd also mainly does along grain boundary to form Mg-Zn-Gd phase for as-cast ZK60-2.98Gd alloy (see Fig. 4).
Fig. 3 EDS spectra of as-cast ZK60-xGd alloys
Table 2 EDS results of as-cast ZK60-xGd alloys
In order to determine new phase owing to the addition of Gd, the comparison between the XRD spectra of as-cast ZK60 and ZK60-2.98Gd alloys was made (see Fig. 5(a)). XRD pattern of as-cast ZK60 alloy consists of the peaks of a-Mg and MgZn2 phases, while the peak of MgZn2 phase disappears completely and that of an unknown phase, not Mg-Gd binary phase, can be observed in the XRD pattern of as-cast ZK60-2.98Gd alloy. Combined with the EDS results, the phase is considered Mg-Zn-Gd phase. This indicates that added Gd atoms combine with Mg and Zn atoms to form Mg-Zn-Gd new phase on the priority.
Figure 5(b) shows the isochronal heating DSC curves of as-cast alloys and peak temperatures TPi of endothermic peaks are listed in Table 3. For as-cast ZK60 alloy, there exist two distinct endothermic peaks with TP1 and TP2 of 618 and 902 K, respectively, which are close to the melting point Tm (620 and 923 K) of MgZn2 phase and a-Mg matrix, respectively. Thus the two endothermic peaks should be corresponding to the melting of the two phases, respectively. For as-cast ZK60-0.28Gd alloy, the two endothermic peaks still exist with the significantly weak first one, and meanwhile two new ones occur at TP=717 and 859 K, respectively. For as-cast ZK60-2.98Gd alloy, the endothermic peak at TP≈620 K disappears completely and another new one occurs at TP=796 K on the basis of retaining others. These endothermic peaks existing at the peak temperature ranging from 717 K to 871 K for the two Gd-containing alloys should be corresponding to the melting of Mg-Zn-Gd phase. This indicates that Mg-Zn-Gd phase exhibits higher thermal stability than MgZn2 phase.
Fig. 4 Distribution of elements in as-cast microstructures of ZK60-xGd alloys
Fig. 5 XRD patterns (a) and isochronal heating DSC curves (b) of as-cast ZK60-xGd alloys
Table 3 Peak temperature TPi on isochronal heating DSC curves of as-cast ZK60-xGd alloys
For as-cast alloys at ambient temperature, tensile strength sb first decreases gradually from 225 MPa for ZK60 alloy to 194 MPa for ZK60-1.00Gd alloy, then slightly increases to 205 MPa for ZK60-2.98Gd alloy, and finally significantly decreases to 172 MPa for ZK60-3.81Gd alloy. Elongation d is not lower than 7.5% when the Gd content does not exceed 2.98% and significantly decreases to 3.5% for ZK60-3.81Gd alloy (see Fig. 6(a)). In summary, as-cast tensile mechanical property is reduced slightly at ambient temperature when the Gd content does not exceed 2.98%. All as-cast tensile fractures exhibit a complex mode of ductile and brittle fractures (see Figs. 7(a), (c) and (e)). Dimples occupy the majority, while cleavage steps are few for as-cast ZK60 alloy. It is contrary for as-cast ZK60-0.28Gd alloy, which is consistent with the slight decrease in elongation. No cleavage step exists for as-cast ZK60-2.98Gd alloy. Meanwhile, the increasing dimples are still less than those for as-cast ZK60 alloy, which is consistent with the slightly lower tensile mechanical property.
3.2 Microstructure and mechanical property of extruded alloy
ZK60, ZK60-0.28Gd and ZK60-2.98Gd alloys were chosen to study the extruded microstructure and tensile mechanical property in the present work. Optical, SEM and TEM images in the extrusion direction are shown in Figs. 8 and 9, respectively. Compared with as-cast state, extruded microstructure is refined significantly owing to the dynamic recrystallization during the hot deformation process. For extruded ZK60 alloy, dynamic recrystallization grains (DRGs) change at a size ranging from 2 μm to 10 μm with average grain size of 5 μm. Meanwhile, majority of MgZn2 phase dissolved through homogenization does not re-precipitate during the hot deformation process. With increasing the Gd content, DRGs become unique and gradually decreases to 2 μm for extruded ZK60-2.98Gd alloy. Meanwhile, broken second phase tends to distribute along the extrusion direction by zonal shape, especially for extruded ZK60-2.98Gd alloy. Diffraction spots to the particles in the two alloys are corresponding to Mg-Zn phase containing relatively high Zr content and Mg-Zn-Gd phase, respectively.
Fig. 6 Tensile mechanical properties of as-cast and extruded ZK60-xGd alloys
Fig. 7 SEM morphologies of tensile fractures at ambient temperature of as-cast and extruded ZK60-xGd alloys (Illustrations indicate the local enlarged regions)
Fig. 8 Optical and SEM images of extruded ZK60-xGd alloys
Compared with as-cast state, extruded tensile mechanical property is enhanced significantly at ambient temperature (see Fig. 6(b)). sb and d reach 355 MPa and 19.5% for extruded ZK60 alloy, respectively. As the Gd content increases, sb increases gradually to 380 MPa for extruded ZK60-2.98Gd alloy, while d gradually decreases to 7%, which is even slightly lower than as-cast value. With increasing the test temperature, sb decreases gradually, while d increases gradually. For extruded ZK60 alloy, sb decreases to 157 and 120 MPa at 448 and 473 K, respectively, indicating that its application temperature can be considered 448 K. Extruded ZK60-0.28Gd alloy exhibits slightly larger values of sb at each test temperature, while extruded ZK60-2.98Gd alloy has significantly larger values where sb reaches 164 MPa at 473 K, indicating that its application temperature is enhanced to 473 K. Meanwhile, the increase in the Gd content still reduces the value of d at each test temperature. In summary, the addition of Gd can enhance the tensile mechanical properties of extruded ZK60 alloy at ambient and elevated temperatures. All extruded tensile fractures at ambient temperature exhibit a typical character of ductile fracture (see Figs. 7(b), (d) and (f)). For extruded ZK60 alloy, uniform and deep dimples occupy the majority and meanwhile few fine second phase particles can be observed at the bottom of dimples. Dimples become less and shallow; meanwhile second phase particles increase gradually with increasing the Gd content, which is consistent with the gradual decrease in the elongation. Combined with as-cast results, the second phase particles at the bottom of dimples should be MgZn2 and Mg-Zn-Gd phases.
4 Discussion
4.1 Phase composition
Degree of difficulty in forming compounds between different elements can be judged by electronegativity difference △χ. The greater the value of Dc is, the larger the binding force is and then the easier the formation of compounds is [14]. Electronegativities c are 1.31, 1.65 and 1.20 for Mg, Zn and Gd, respectively [15], thus Zn-Gd exhibits the larger value of Dc (0.45) than Mg-Zn and Mg-Gd (0.34 and 0.11, respectively). Thus the binding force between Zn and Gd is greater and then Gd atoms combine with Zn and Mg atoms to form Mg-Zn-Gd new phase on the priority. During the process, majority of Zn and Gd atoms in the melt are consumed, which leads to no formation of MgZn2 and Mg-Gd binary phases for as-cast ZK60-2.98Gd alloy.
Fig. 9 TEM images of extruded ZK60-xGd alloys
4.2 Grain refinement and dynamic recrystallization
Grain refinement mechanism for magnesium alloys has not been entirely clear, and is different for different grain-refining methods. However, the basic starting point is to increase the nucleation rate and inhibit the growth of crystal nuclei. Solute with good segregation and effective nuclei are two essential factors. Solute with good segregation can generate a solutal undercooling at the liquid-solid interface, impede the growth of dendrite and provide the driving force to activate the nucleation. Nucleation capacity of particles determines the start of solidification and number of effective nuclei at the solutal undercooling region. Role of solute element can be expressed by growth restriction factor (GRF) [2]. Metal Gd exhibits a relatively large value of growth restriction factor coefficient m(k-1), i.e. 5.52, thus as-cast microstructure is refined effectively.
As a softening and grain-refining mechanism for magnesium alloys, dynamic recrystallization can ameliorate the deformation microstructure and enhance the mechanical property effectively. In the present work, large extrusion ratio (l=40) can generate large deformation degree and ensure the adequate refinement of microstructures and occurrence of dynamic recrystallization. Middle extrusion temperature (T=593 K) can guarantee not only the occurrence of dynamic recrystallization, but also no easy growth of dynamic recrystallization grains. Low extrusion speed (v=1-2 m/min) can guarantee the fine microstructure. As-cast microstructure is refined gradually with increasing the Gd content, thus the corresponding dynamic recrystallization grains become smaller gradually. Increasing Mg-Zn-Gd phase can hinder the growth of dynamic recrystallization grains. Taking the above mentioned aspects into account, extruded microstructure is refined obviously and decreases gradually to 2 μm for ZK60-2.98Gd alloy.
4.3 Strengthening mechanism
Grain-refinement strengthening is generally considered the most important strengthening mechanism for magnesium alloys [16,17]. In the present work, as-cast microstructure is refined gradually with increasing the Gd content. However, the second phase increases gradually and tends to distribute along grain boundary by continuous network. During the tensile deformation, dislocation piles at the interface between the second phase particles and a-Mg matrix, which can lead to stress concentration. When the stress reaches a certain degree, porosities and micro-cracks initiate. When it exceeds fracture strength, the matrix begins to fracture locally and plasticity decreases significantly. Taking the two aspects into account, as-cast tensile mechanical property is reduced slightly at ambient temperature when the Gd content does not exceed 2.98%.
Extruded microstructure is refined significantly and decreases gradually to 2 μm for ZK60-2.98Gd alloy, which leads to the prominent strengthening effect by grain refinement. Meanwhile, fine dispersed Mg-Zn-Gd phase can also play a certain role in dispersion strengthening. Therefore, sb of extruded alloys can be enhanced significantly at ambient temperature and meanwhile increases gradually with increasing the Gd content. In addition, grain refinement is also favorable for ameliorating the ductility. However, the increasing second phase tends to distribute along the extrusion direction parallel to the tensile direction by zonal shape. Thus, d decreases gradually and even is slightly lower than as-cast value for ZK60-2.98Gd alloy. In summary, it is necessary to avoid the distribution by continuous network in as-cast state or zonal shape in extruded state for the second phase.
For extruded ZK60 alloy, rare MgZn2 phase plays a weak role in pinning grain boundary and leads to relatively poor tensile mechanical property at elevated temperatures, where sb decreases to 120 MPa at 473 K. With increasing the Gd content, Mg-Zn-Gd phase with a higher thermal stability increases gradually. Grain boundary is pinned and its sliding is inhibited effectively. Therefore, the tensile mechanical property is enhanced effectively at elevated temperatures, where sb increases to 164 MPa at 473 K for extruded ZK60-2.98Gd alloy.
5 Conclusions
1) With increasing the Gd content, Mg-Zn-Gd phase increases gradually, and MgZn2 phase gradually decreases to disappear. The second phase tends to distribute along grain boundary by continuous network. As-cast microstructure is refined gradually, while its tensile mechanical property is reduced slightly at ambient temperature when the Gd content does not exceed 2.98%.
2) Extruded microstructure is refined obviously and decreases gradually to 2 μm for ZK60-2.98Gd alloy. Broken second phase tends to distribute along the extrusion direction by zonal shape. Extruded tensile mechanical property is enhanced significantly at ambient and elevated temperatures. sb values at 298 and 473 K gradually increase from 355 and 120 MPa for ZK60 alloy to 380 and 164 MPa for ZK60-2.98Gd alloy, respectively. Extruded tensile fractures exhibit a typical character of ductile fracture.
3) The grain size and the second phase play an important role in affecting the tensile mechanical property of magnesium alloys.
References
[1] POLMEAR I J. Magnesium alloys and applications [J]. Materials Science and Technology, 1994, 10(1): 1-16.
[2] CHEN Zhen-hua. Wrought magnesium alloys [M]. Beijing: Chemical Industry Press, 2005: 1-20. (in Chinese)
[3] LUO Z P, SONG D Y, ZHANG S Q. Strengthening effects of rare earths on wrought Mg-Zn-Zr-RE alloys [J]. Journal of Alloys and Compounds, 1995, 230(2): 109-114.
[4] XU D K, LIU L, XU Y B, HAN E H. Effect of microstructure and texture on the mechanical properties of the as-extruded Mg-Zn-Y- Zr alloys [J]. Materials Science and Engineering A, 2007, 443(1-2): 248-256.
[5] LUO Su-qin, TANG Ai-tao, PAN Fu-sheng, SONG Kai, WANG Wei-qing. Effect of mole ratio of Y to Zn on phase constituent of Mg-Zn-Zr-Y alloys [J]. Transactions of Nonferrous Metals Society of China, 2011, 21(4): 795-800.
[6] XIA C Q, WANG Y N, WU A R, GU Y. Effects of cerium on microstructure and mechanical properties of ZK60 alloy [J]. Journal of Central South University of Technolnolgy, 2005, 12(5): 515-520.
[7] LI Q, WANG Q D, WANG Y X, ZENG X Q, DING W J. Effect of Nd and Y addition on microstructure and mechanical properties of as-cast Mg-Zn-Zr alloy [J]. Journal of Alloys and Compounds, 2007, 427(1-2): 115-123.
[8] YU W B, LIU Z Y, HE H, CHENG N P, LI X L. Microstructure and mechanical properties of ZK60-Yb magnesium alloys [J]. Materials Science and Engineering A, 2008, 478(1-2): 101-107.
[9] LI D Q, WANG Q D, DING W J. Effects of Ho on the microstructure and mechanical properties of Mg-Zn-Ho-Zr magnesium alloys [J]. Rare Metals, 2011, 30(2): 131-136.
[10] LI L, ZHANG X M, TANG C P, DENG Y L, ZHOU N. Mechanical properties and deep drawability of Mg-Gd-Y-Zr alloy rolling sheet at elevated temperatures [J]. Materials Science and Engineering A, 2010, 527(4-5): 1266-1274.
[11] ZHANG X M, LI L, DENG Y L, ZHOU N. Superplasticity and microstructure in Mg-Gd-Y-Zr alloy prepared by extrusion [J]. Journal of Alloys and Compounds, 2009, 481(1-2): 296-300.
[12] LIU X B, CHEN R S, HAN E H. Effects of ageing treatment on microstructures and properties of Mg-Gd-Y-Zr alloys with and without Zn additions [J]. Journal of Alloys and Compounds, 2008, 465(1-2): 232-238.
[13] WANG J, MENG J, ZHANG D P, TANG D X. Effect of Y for enhanced age hardening response and mechanical properties of Mg-Gd-Y-Zr alloys [J]. Materials Science and Engineering A, 2007, 456(1-2): 78-84.
[14] HUANG Zheng-hua, GUO Xue-feng, ZHANG Zhong-ming. Effects of Ce on damping capacity of AZ91D magnesium alloy [J]. Transactions of Nonferrous Metals Society of China, 2004, 14(2): 311-315.
[15] DEAN J A. Lange’s handbook of chemistry [M]. 15th ed. New York: McGraw-Hill Press, 1999: 429-435.
[16] HALL E O. The deformation and ageing of mild steel: III. Discussion of results [J]. Proceedings of the Physical Society Section B, 1951, 64: 747-753.
[17] PETCH N J. The cleavage strength of polycrystals [J]. Journal of Iron Steel Institute, 1953, 174: 25-28.
黄正华,戚文军,徐 静
广州有色金属研究院,广州 510650
摘 要:分析了铸态和挤压态ZK60-xGd(x=0~4)合金的组织和相组成,测试了其拉伸力学性能。结果表明,随着Gd含量的增加,铸态组织逐渐细化,Mg-Zn-Gd新相逐渐增多,而MgZn2相逐渐减少直至消失,第二相趋于连续网状分布于晶界处;当Gd含量不超过2.98%时,铸态室温拉伸力学性能稍降低。经挤压比l=40和挤压温度 T=593 K的挤压后,组织显著细化,平均晶粒尺寸逐渐减至ZK60-2.98Gd合金的2 μm,破碎的第二相沿着挤压方向呈带状分布;挤压态的拉伸力学性能均显著提高:298和473 K时的抗拉强度分别从ZK60合金的355和120 MPa逐渐提高至ZK60-2.98Gd合金的380和164 MPa。挤压态拉伸断口呈现典型的韧性断裂特征。
关键词:ZK60镁合金;Gd变质;挤压;显微组织;力学性能;断口形貌
(Edited by Hua YANG)
Foundation item: Projects (2010A090200078, 2011A080403008) supported by the Major Science and Technology Project of Guangdong Province, China
Corresponding author: Zheng-hua HUANG; Tel/Fax: +86-20-37238039; E-mail: zhhuang@live.cn
DOI: 10.1016/S1003-6326(13)62769-7