方 剛,徐向濤,全恩臣,戴特力,范嗣強,張 鵬
摻Y(jié)b3+雙包層光纖激光器的研究進展
方 剛,徐向濤,全恩臣,戴特力,范嗣強,張 鵬*
(重慶師范大學物理與電子工程學院重慶市高校光學工程重點實驗室,重慶400047)
光纖激光器以其獨特的優(yōu)勢得到快速發(fā)展,其應用范圍已經(jīng)擴展到工業(yè)加工、國防軍事、醫(yī)療等領域。綜述了連續(xù)、脈沖摻Y(jié)b3+雙包層光纖激光器的國內(nèi)外研究進展,介紹了利用能承受高功率的合束器、光纖光柵的全光纖激光器,利用種子光主振蕩光纖放大技術產(chǎn)生高光束質(zhì)量、高平均功率、高峰值功率的脈沖光纖激光器。分析了影響光纖激光器功率提高的因素,如光纖的損傷、非線性性效應、熱效應。最后,對摻Y(jié)b3+雙包層光纖激光器的發(fā)展前景進行了展望。
激光器;摻Y(jié)b3+;光纖激光器;主振蕩功率放大器;雙包層光纖;非線性效應;熱效應
近年來,隨著大功率半導體激光器抽運技術和雙包層光纖制作工藝的不斷提高,光纖激光器的輸出功率已經(jīng)從最初的幾百毫瓦上升到數(shù)千瓦水平[1]。光纖激光器以摻稀土元素(鉺、鐿、銩等)光纖作為增益介質(zhì),鐿離子(Yb3+)具有很高的吸收截面、較寬的吸收光譜(典型值為910nm~980nm)以及較寬的發(fā)射光譜(典型值為1040nm~1100nm),使其具有90%以上的量子效率,并且可以高濃度摻雜[2];雙包層光纖的出現(xiàn)使得抽運光的耦合效率相當高,因此,摻Y(jié)b3+雙包層光纖激光器的綜合電光效率高達20%以上,光光轉(zhuǎn)換效率超過60%;雙包層光纖激光器的輸出光束質(zhì)量由纖芯直徑d和數(shù)值孔徑dNA決定,只要光纖的歸一化頻率ν<2.404,就可以實現(xiàn)基橫模運轉(zhuǎn),獲得近衍射極限的激光輸出;雙包層光纖的表面積/體積比很大,使得散熱性能非常好,對于百瓦級的光纖激光器,采用風扇制冷即可,對于千瓦級的光纖激光器,需采用水冷[3];能承受高功率的合束器、光纖光柵的出現(xiàn),使得雙包層光纖激光器實現(xiàn)全光纖化成為可能,極大地簡化了光纖激光器的設計和制作,而且光纖具有極好的柔繞性,因此光纖激光器的結(jié)構(gòu)簡單、緊湊靈活。本文中對近年來國內(nèi)外連續(xù)、脈沖摻Y(jié)b3+雙包層光纖激光器的研究進展進行了綜述,對影響光纖激光器功率提高的影響因素進行了分析。
1.1分離器件連續(xù)摻Y(jié)b3+雙包層光纖激光器進展
普通單模光纖纖芯直徑很小,一般在10μm以下,使得高功率抽運光很難耦合到纖芯。由于要將半導體激光抽運入單模光纖的纖芯,要求半導體激光也必須是單模,而單模半導體激光能達到的功率小于1W,這使得普通單模摻雜光纖激光器難以實現(xiàn)高功率輸出。為了提高抽運轉(zhuǎn)換效率和光纖激光的輸出功率,1988年,美國Polariod公司的SNITZER等人提出雙包層光纖[4],雙包層光纖比普通光纖多了一個內(nèi)包層,內(nèi)包層可以傳輸多模抽運光,其提高了抽運光的入纖功率和耦合效率,同時也降低了對抽運源的要求。隨后,光纖激光器的輸出功率不斷攀升。2004年,英國南安普頓大學的JEONG等人報道了高效包層抽運的摻Y(jié)b3+光纖激光器,其在波長1100nm處獲得了1.36kW接近衍射極限連續(xù)激光輸出,斜率效率達83%[5]。實驗裝置如圖1所示。該摻鐿光纖激光器使用雙端抽運結(jié)構(gòu),抽運光通過準直聚焦后耦合到有源光纖,耦合效率大于90%,采用光纖端面自由輸出并用V型槽對光纖兩端進行散熱。2009年,該研究小組又報道了利用同樣的抽運結(jié)構(gòu)在波長1.1μm處獲得了大于2.1kW的連續(xù)激光輸出,斜率效率達74%,光束質(zhì)量M2<1.2[6]。2010年,HE等人報道了采用優(yōu)化結(jié)構(gòu)參量的大芯徑雙包層光纖,雙端抽運在波長1090nm處獲得了l.75kW的連續(xù)波激光輸出,斜率效率為76%[7]。
圖1 兩個二極管堆抽運的摻鐿光纖激光器實驗裝置圖
1.2全光纖連續(xù)摻Y(jié)b3+雙包層光纖激光器進展
由于全光纖激光技術具有效率高、光束質(zhì)量好、穩(wěn)定性高、免維護等優(yōu)點,近幾年來,利用能夠承受高功率的光纖合束器、光纖光柵等無源器件使摻Y(jié)b3+雙包層光纖激光器全光纖化成為最受關注的研究熱點之一。IPG公司利用其發(fā)明的光纖側(cè)向抽運技術,不斷刷新高功率光纖激光器的記錄,其產(chǎn)品涵蓋50kW以內(nèi)連續(xù)光纖激光器和200W以內(nèi)脈沖光纖激光器。英國SPI公司的產(chǎn)品涵蓋1000W以內(nèi)連續(xù)光纖激光器和平均功率70W、峰值功率20kW以內(nèi)脈沖光纖激光器。IPG高功率全光纖激光器的結(jié)構(gòu)如圖2所示。二極管激光器(laser diode,LD)的抽運光通過光纖合束器進入有源光纖,諧振腔鏡分別由高反射率的光纖光柵和低反射率的光纖光柵組成。來看看國內(nèi)的情況:2011年,國防科學技術大學的MA等人報道了一臺全光纖結(jié)構(gòu)的光纖激光器,采用雙端抽運結(jié)構(gòu),獲得了1008W的高功率輸出,輸出激光的中心波長為1082nm,光轉(zhuǎn)換效率為68%[8];武漢銳科公司目前已經(jīng)實現(xiàn)了4000W以內(nèi)連續(xù)和100W以內(nèi)脈沖光纖激光器產(chǎn)品。
圖2 IPG高功率全光纖激光器結(jié)構(gòu)圖
單根光纖激光的輸出功率達數(shù)千瓦時,其熱效應及非線性效應成為提高功率的限制因素,采用光纖放大器的方式可以實現(xiàn)更高功率的激光輸出[9]。IPG公司己經(jīng)實現(xiàn)了10kW的單模、50kW的多模連續(xù)激光輸出[8]。IPG公司10kW單模光纖激光器結(jié)構(gòu)如圖3所示。主振蕩器輸出功率為1kW、波長為1070nm的種子光,注入纖芯為30μm、長15m的摻鐿光纖,放大級由45個波長為1018nm、功率為300W的光纖激光器進行抽運,采用1018nm的光纖激光器作為抽運源,雖然摻鐿光纖在1018nm的吸收峰比975nm的吸收峰小約1個量級,但是1018nm光纖激光器的亮度要高約100倍(相比30W多模半導體激光器),而且量子缺陷低于5%,從而進一步減少了熱量的產(chǎn)生。再看看國內(nèi)的情況:2009年,中國科學院西安光學精密機械研究所在全光纖激光技術研究方面取得重要階段性進展,采用兩級放大輸出功率超過1000W,光光轉(zhuǎn)換效率為62%[10];2011年,清華大學的YAN等人報道了采用主振蕩器-二級放大的全光纖激光器,在波長1085nm處輸出功率1.1kW,光光效率達66.5%[11]。
圖3 IPG 10kW單模光纖激光器結(jié)構(gòu)圖
1.3脈沖摻Y(jié)b3+雙包層光纖激光器進展
與連續(xù)光纖激光器相比,脈沖光纖激光器具有高的峰值功率,能夠與物質(zhì)高效地進行相互作用,在激光加工、激光雷達、光譜學等領域有著廣泛的應用。目前摻Y(jié)b3+雙包層光纖激光器產(chǎn)生脈沖的方式主要通過調(diào)Q、鎖模、種子光主振蕩功率放大(master oscillator power amplifier,MOPA)3種技術。調(diào)Q能產(chǎn)生高能量脈沖,脈沖寬度在納秒級的激光;鎖模輸出激光脈沖寬度可達飛秒級。種子光主振蕩光纖功率放大是用高光束質(zhì)量較低功率的脈沖激光作為種子光源,注入到二極管激光器抽運的有源光纖內(nèi)進行功率放大,實現(xiàn)高能量的脈沖激光輸出,產(chǎn)生的激光波長、重復頻率與種子光相同[12]。調(diào)Q、鎖模獲得的平均功率往往并不一定很大,因此,采用種子光主振蕩光纖功率放大技術是實現(xiàn)高光束質(zhì)量、高平均功率、高峰值功率的超短脈沖激光的理想方法。典型主振蕩光纖功率放大器結(jié)構(gòu)如圖4所示,一般采用三級放大,既可以抑制光纖放大器產(chǎn)生的自發(fā)輻射,還可以實現(xiàn)高功率超短脈沖激光輸出,光纖放大器之間的隔離器防止反向傳輸光影響種子光和前一級光纖放大器的工作特性,避免損害種子源和各級放大器的抽運源。
圖4 典型主振蕩光纖功率放大器結(jié)構(gòu)圖
2002年,德國耶拿大學的LIMPERT等人報道了在重復頻率50kHz時,獲得平均功率達100W;在重復頻率3kHz時,獲得脈沖能量4mJ,接近衍射極限的脈沖激光輸出[13]。2005年,CHENG等人報道了在脈沖寬度為500ns和50ns時,脈沖能量分別達82mJ和27mJ;在脈沖寬度4ns時,峰值功率達2.4MW的摻Y(jié)b3+光纖放大器系統(tǒng)[14];LIMPERT報道了在重復頻率47MHz時,獲得了平均功率97W、峰值功率200kW、脈沖寬度10ps的超短脈沖輸出[15];ROSER報道了在重復頻率73MHz時,獲得了平均功率為131W、脈沖寬度為220fs、峰值功率8.2MW的超短脈沖輸出[16]。2006年,英國南安普頓大學DUPRIEZ等人報道了利用主振蕩功率放大獲得了平均功率超過300W、重復頻率為1GHz、脈沖寬度為20ps的超短脈沖激光輸出[17]。2007年,ROSER報道了采用聲光調(diào)制器獲得了平均功率90W,脈沖能量100mJ,重復頻率900kHz,脈沖寬度500fs超短激光輸出[18]。2010年,CHEN等人報道了在重復頻率908MHz時,獲得了平均功率100W、脈沖寬度21ps的線偏振、單模摻鐿皮秒光纖激光放大系統(tǒng)[19]。
2006年,中國科學院上海光學精密機械研究所KONG等人報道了在重復頻率100kHz時,獲得了平均功率133.8W、脈沖寬度400ns的輸出;在重復頻率20kHz時,獲得了脈沖寬度15ns,峰值功率達300kW[20]。清華大學YE等人報道了利用固體激光器種子源在摻Y(jié)b3+雙包層保偏光纖進行放大,在重復頻率40kHz時,獲得了脈沖能量大于1mJ、脈沖寬度30ns、平均功率53.1W的線偏振脈沖激光[21]。2009年,LI等人報道了利用激光二極管抽運聲光開關調(diào)制的固體激光器作為種子源,獲得了平均功率為103W的脈沖激光輸出,重復頻率50kHz,脈寬12.7ns,峰值功率達162kW[22]。ZHOU等人報道了采用MOPA方式,獲得了平均功率150W的脈沖放大激光輸出[23]。2011年,YANG等人報道了在重復頻率40MHz時,平均功率150W、脈沖寬度273fs的超短脈沖輸出[24]。2012年,SU等人報道了主振蕩功率放大MOPA結(jié)構(gòu)的單頻脈沖全光纖激光器,對重復頻率10MHz、脈寬約8ns、平均功率約0.5mW的單頻脈沖種子光進行級聯(lián)放大,獲得了平均功率300.8W的高功率激光輸出[25]。
隨著摻Y(jié)b3+雙包層光纖激光器功率的不斷提高,光纖的損傷、非線性效應、熱效應已經(jīng)成為光纖激光器功率提高的制約因素。
2.1光纖的損傷
雙包層光纖激光器要實現(xiàn)單模激光輸出,其纖芯相關參量必須滿足歸一化頻率:
式中,d為纖芯的直徑,λ為信號光波長。由(1)式可知,雙包層光纖激光器能否實現(xiàn)單模輸出,是由纖芯的直徑和數(shù)值孔徑?jīng)Q定的。純石英光纖的損傷閾值非常高,在脈沖激光下的損傷閾值約為100W/μm2(峰值功率)。對于摻雜石英光纖,由于摻雜降低了石英的均勻性和純度,大大降低了光纖的損傷閾值,為了保證雙包層光纖激光器性能穩(wěn)定可靠避免燒毀光纖,在實際光纖激光器設計時一般取1.5W/μm2作為光纖的損傷閾值[1]。對于滿足歸一化頻率實現(xiàn)單模輸出的雙包層光纖,6μm~10μm的纖芯直徑可實現(xiàn)的最高功率在百瓦量級,這就限制了雙包層光纖激光器在滿足單模激光輸出的條件下實現(xiàn)高功率輸出。大模場光纖的出現(xiàn),有效地解決了這一問題:通過增大纖芯直徑和減小數(shù)值孔徑來實現(xiàn)低階模、高光束質(zhì)量的激光輸出。但是由于光纖材料的限制,纖芯數(shù)值孔徑也不可能很小,目前技術上可實現(xiàn)其典型值為0.06,對應的最大光纖單模直徑在17μm左右。
2.2非線性效應
在雙包層光纖激光器中,由于纖芯所能承受的激光功率密度很高,且光纖長度較長,容易引起非線性效應,其主要包括受激布里淵散射(stimulated Brillouin scattering,SBS)、受激喇曼散射(stimulated Raman scattering,SRS)[26]。對于這兩種非線性效應的產(chǎn)生,都存在一個閾值,當激光功率大于閾值時,才會出現(xiàn)非線性散射。SBS和SRS作為損耗都會影響高功率的激光輸出[27]。
SBS散射的閾值表達式為:
SRS散射的閾值表達式為:
式中,gB和gR分別是布里淵增益系數(shù)和喇曼增益系數(shù),其典型值分別為gR=1×10-13m/W,gB=5× 10-11m/W,Aeff是光纖纖芯的有效模場面積,Leff是光纖的有效長度,ΔνB是SBS增益帶寬,一般為30MHz~50MHz,Δνs是信號光的帶寬??梢姽饫w中受激布里淵散射的閾值比受激喇曼散射的閾值要低約兩個數(shù)量級,根據(jù)(2)式可知,通過增大纖芯面積和減小光纖長度,增大受激布里淵散射的閾值,有利于高功率激光的輸出。
2.3熱效應
由于雙包層光纖激光器使用細長的摻雜光纖作為增益介質(zhì),有大的表面積/體積比,因此散熱性能很好。但是光纖激光器的輸出功率達到幾百瓦以上時,熱效應問題就成為一個不可忽視的重要問題[28]。
雙包層光纖激光器的徑向溫度為:
雙包層光纖激光器的軸向溫度為:
式中,Tc為環(huán)境溫度;a,b,c分別為雙包層光纖的纖芯、內(nèi)包層和外包層的半徑;相應的熱導率分別為κ1,κ2,κ3;λp,λs分別為抽運光和信號激光的波長;α0為光纖對抽運光的吸收系數(shù);P0為抽運光功率;L為光纖的長度。由(4)式與(5)式比較可知,雙包層光纖激光器的徑向溫度梯度遠大于軸向溫度梯度,這表明在熱傳導過程中,熱量主要是沿徑向傳輸?shù)摹S桑?)式可知,雙包層光纖激光器的效應與環(huán)境溫度、光纖對光纖對抽運光的吸收系數(shù)、抽運功率、包層直徑等參量有關。光纖對抽運光的吸收系數(shù)越小,光纖內(nèi)的溫度相對較低,高功率的雙包層光纖激光器,可以選擇對抽運光吸收系數(shù)小的雙包層光纖;由于雙包層光纖的外包層為聚合物材料,其導熱系數(shù)κ3比內(nèi)包層和纖芯的導熱系數(shù)要低很多,因此選擇大半徑尺寸的纖芯和內(nèi)包層、小半徑的外包層有利于光纖激光器的散熱。由(5)式可知,雙包層光纖的溫度在抽運光注入端最高,隨著光纖長度L的增大,雙包層光纖的溫度越來越小。因此,選擇合適的雙包層光纖對激光器的性能至關重要。
與傳統(tǒng)的固體激光器相比,光纖激光器具有轉(zhuǎn)換效率高、熱管理方便、輸出光束質(zhì)量好、結(jié)構(gòu)簡單緊湊等優(yōu)點,在通信、醫(yī)療、工業(yè)、軍事等領域有著廣闊的應用前景。隨著半導體抽運激光、特種光纖、能承受高功率低損耗的合束器、光纖光柵、功率提升的相干合成、單頻激光的功率提升、摻鐿特殊波長光纖激光器975nm,1018nm,1030nm,1100nm等以及摻鐿光纖激光器的頻率變換技術等關鍵器件技術的不斷提高,光纖激光器將克服光纖損傷、非線性效應、熱效應等困難,實現(xiàn)全光纖化,輸出功率從百瓦級向千瓦級甚至萬瓦發(fā)展,高功率連續(xù)激光向高平均功率、高峰值功率的脈沖光纖激光發(fā)展。可以預見,光纖激光器將成為其它激光器的有力競爭對手,必將在未來各個領域中發(fā)揮重要作用。
[1] LOU Q H.High power fiber lasers and its applications[M].Beijing:Press of University of Science and Technology of China,2010:8-12(in Chinese).
[2] NILLSSON J,PAYNE D N.High-power fiber lasers[J].Science,2011,332(6032):921-922.
[3] CHEN Z L,HOU J,JIANG Z F.Estimation of maximum output power of double clad fiber laser[J].High Power Laser and Particle Beams,2007,19(4):577-580(in Chinese).
[4] SNITZER E,PO H,HAKIMIF,et al.Double clad,offset core Nd fiber laser[C]//Optical Fiber Sensors.Washington DC,USA:Optical Society of America,1988:PD5.
[5] JEONG Y,SAHU J K,PAYNE D N,et al.Ytterbium-doped large-core fiber laser with 1.36kW continuous-wave output power[J].Optics Express,2004,12(25):6088-6092.
[6] JEONG Y,BOYLAND A J,SAHU JK,et al.Multi-kilowatt single-mode ytterbium-doped large-core fiber laser[J].Journal of the Optical Society of Korea,2009,13(4):416-422.
[7] HE B,ZHOU J,LOUQH,et al.1.75-kilowatt continuous-wave output fiber laser using homemade ytterbium-doped large-core fiber[J].Microwave and Optical Technology Letters,2010,52(7):1668-1671.
[8] MA Y X,XIAO H,ZHOU P,et al.kW level fiber laser in all fiber format[J].High Power Laser and Particle Beams,2011,23(5):1137-1138(in Chinese).
[9] WUW M,XIAOH,XU JM,et al.Research progress of tandempumped fiber lasers[J].Laser&Optoelectronics Progress,2011(9):11-18(in Chinese).
[10] DUAN K L,ZHAOB Y,ZHAOW,et al.1000W all-fiber laser[J].Chinese Journal of Lasers,2009,36(12):3219(in Chinese).
[11] YAN P,YIN Sh P,HE JW,et al.1.1kW ytterbium monolithic fiber laser with assembled end-pump scheme to couple high brightness single emitters[J].IEEE Photonics Technology Letters,2011,23(11):697-699.
[12] KOECHNERW.Solid-state laser engineering[M].Beijing:Beijing Publishing House,2002:410-505(in Chinese).
[13] LIMPERT J,HOFER S,LIEM A,et al.100W average-power,high-energy nanosecond fiber amplifier[J].Applied Physics,2002,B75(4/5):477-479.
[14] CHENGM,CHANG Y,GALVANAUSKASA,etal.High-energy and high-peak-power nanosecond pulse generation with beam quality control in 200μm core highly multimode Yb-doped fiber amplifiers[J].Optics Letters,2005,30(4):358-360.
[15] LIMPERT J,DEGUIL-ROBIN N,MANEK-H?NNINGER I,et al.High-power picosecond fiber amplifier based on nonlinear spectral compression[J].Optics Letters,2005,30(7):714-716.
[16] ROSER F,ROTHHARD J,ORTAC B,et al.131W 220 fs fiber laser system[J].Optics Letters,2005,30(20):2754-2756.
[17] DUPRIEZ P,PIPER A,MALINOWSKIA,et al.High average power,high repetition rate,picosecond pulsed fiber master oscillator power amplifier source seeded by a gain-switched laser diode at 1060nm[J].IEEE Photonics Technology Letters,2006,18(9):1013-1015.
[18] ROSER F,SCHIMPFD,SCHMIDTO,et al.90W average power 100μJ energy femtosecond fiber chirped-pulse amplification system[J].Optics Letters,2007,32(15):2230-2232.
[19] CHEN K K,PRICE JV,ALAM V,et al.Polarisation maintaining 100W Yb-fiber MOPA producingμJpulses tunable in duration from 1 to 21ps[J].Optics Express,2010,18(14):14385-14394.
[20] KONG L F,LOU Q H,ZHOU J,et al.133W pulsed fiber amplifier with large-mode-area fiber[J].Optical Engineering,2006,45(1):010502.
[21] YE Ch G,GONGM L,YAN P,et al.Linearly-polarized singletransverse-mode high-energy multi-ten nanosecond fiber amplifier with 50W average power[J].Optics Express,2006,14(17):7604-7609.
[22] LIY,ZHU Ch,WANG X F,et al.Experimental study on hundred-watt output power high repetition rate narrow pulse duration fiber lasers[J].Chinese Journal of Lasers,2009,36(2):281-284(in Chinese).
[23] ZHOU J,DU ST,LIU X,etal.150W high repetition frequency and narrow pulse fiber lasers with the MOPA way[J].Chinese Journal of Lasers,2009,36(4):861(in Chinese).
[24] YANG KW,HAO Q,LIW X,et al.High-power fem tosecond pulses fiber laser system[J].Infrared and Laser Engineering,2011,40(7):1254-1256(in Chinese).
[25] SU R T,ZHOU P,WANG X L,etal.Single-frequency nanosecond pulsed laser with output power of 300W in all fiber format[J].High Power Laser and Particle Beams,2012,24(5):1009-1010(in Chinese).
[26] YANG Ch B,LENG JY,LU Q Sh.The theoretical analysis of SBSof Yb-doped double-clad in single-frequency fiber amplifier[J].Laser Technology,2011,35(1):117-121(in Chinese).
[27] AGRAWAL G.Applications of nonlinear fiber optics[M].2nd ed.New York,USA:Academic Press,2001:179-244.
[28] BROWN D C,HOFFMAN H J.Thermal,stress,and thermo-optic effects in high average power double-clad silica fiber lasers[J].IEEE Journal of Quantum Electronics,2001,37(2):207-217.
Research progress of Yb-doped double-clad fiber lasers
FANG Gang,XUXiangtao,QUANEnchen,DAI Teli,F(xiàn)ANSiqiang,ZHANGPeng
(Optical Engineering Key Laboratory of Chongqing City,College of Physics and Electronic Engineering,Chongqing Normal University,Chongqing 400047,China)
Fiber lasers have grown rapidly for their unique advantages,which become very competitive in many applications such as industrial processing,defense and military,medical and other fields.The progress of continuous wave and pulse Yb-doped double-clad fiber laser at home and abroad was reported.All-fiber lasers with combiners and fiber gratings were introduced.Simultaneously,pulse fiber lasers with high beam quality,high average power and high peak power based on master oscillator fiber power amplifier technology were presented.The influencing factors of increasing power of fiber laser,including damage of the optical fiber,nonlinear effects,thermal effects,were analyzed.Finally,the development trends of the Yb-doped double-clad fiber lasers were predicated.
lasers;Yb-doped;fiber laser;master oscillator fiber power amplifier;double-clad fiber;nonlinear effect;thermal effect
TN248.1
A
10.7510/jgjs.issn.1001-3806.2014.02.028
1001-3806(2014)02-0278-05
重慶市自然科學基金資助項目(cstcjjA40029);重慶師范大學博士啟動基金資助項目(11XLB014)
方 剛(1989-),男,碩士研究生,主要從事微納光子器件及高功率光纖激光器方面的研究。
*通訊聯(lián)系人。E-mail:zhangpeng2010@cqnu.edu.cn
2013-04-24;
2013-05-15