劉國(guó)輝,肖 峰,李英杰
(1.中建交通建設(shè)集團(tuán)有限公司,北京 100142;2.北京交通大學(xué) 土木建筑工程學(xué)院,北京 100044)
?
地鐵車站蓋挖逆作法施工結(jié)構(gòu)內(nèi)力和變形數(shù)值模擬
劉國(guó)輝1,肖峰2,李英杰1
(1.中建交通建設(shè)集團(tuán)有限公司,北京 100142;2.北京交通大學(xué) 土木建筑工程學(xué)院,北京 100044)
摘要:依托北京地鐵某車站蓋挖逆作法施工過(guò)程,通過(guò)有限元數(shù)值模擬研究頂板、中板、底板和側(cè)墻的受力和變形狀態(tài)隨施工過(guò)程的變化規(guī)律,分析中間樁(柱)受力隨施工過(guò)程的變化過(guò)程。計(jì)算表明,隨施工過(guò)程進(jìn)行,中柱的應(yīng)力變化情況為先處于受拉狀態(tài),在上覆土回填之后又轉(zhuǎn)為受壓狀態(tài),而頂板上側(cè)在靠近側(cè)墻和中間柱附近位置處于受拉狀態(tài),下側(cè)則相反。此外,還研究了施工過(guò)程中圍護(hù)結(jié)構(gòu)、中間柱差異沉降規(guī)律,表明中柱間最大差異沉降以及中間柱與側(cè)墻最大差異沉降均小于規(guī)范允許值。最后對(duì)車站施工期間的安全性以及對(duì)周圍建筑物的影響進(jìn)行了評(píng)價(jià)。
關(guān)鍵詞:蓋挖逆作法;數(shù)值計(jì)算;結(jié)構(gòu)內(nèi)力;地層沉降
地鐵車站蓋挖逆作法施工過(guò)程中,車站結(jié)構(gòu)會(huì)涉及到復(fù)雜的受力和變形狀態(tài),以及力學(xué)轉(zhuǎn)換過(guò)程。已有學(xué)者采用有限元等數(shù)值計(jì)算方法,對(duì)逆作法施工進(jìn)行了模擬[1],并對(duì)不同設(shè)計(jì)方案及其支護(hù)結(jié)構(gòu)進(jìn)行穩(wěn)定性分析和參數(shù)優(yōu)化[2]。此外,一些研究[3]還對(duì)地鐵施工現(xiàn)場(chǎng)的車站主體土方開(kāi)挖、中間樁柱施工、差異沉降控制、防水處理等問(wèn)題進(jìn)行了分析,并分析結(jié)構(gòu)和地層的最終沉降量[4]。蘇潔等[5]以北京地鐵4號(hào)線動(dòng)物園車站施工過(guò)程為例,考慮蓋挖逆作施工過(guò)程中地層及結(jié)構(gòu)動(dòng)態(tài)力學(xué)響應(yīng)過(guò)程,采用有限元計(jì)算軟件對(duì)施工過(guò)程中車站結(jié)構(gòu)變形進(jìn)行模擬分析,提出車站施工過(guò)程的分階段控制方法。魏玉超等[6]就城市地鐵車站工程的蓋挖逆作法施工技術(shù)進(jìn)行了介紹。而鄭騏[7]和趙育紅[8]重點(diǎn)分析了在地鐵車站采用蓋挖逆作法進(jìn)行中間柱施工的工藝原理、工藝流程以及質(zhì)量和安全措施??傮w而言,由于地鐵車站蓋挖逆作法施工過(guò)程十分復(fù)雜,施工力學(xué)過(guò)程仍有待于進(jìn)行深入討論。
本文以北京地鐵14號(hào)線東湖渠站為工程背景,通過(guò)數(shù)值模擬手段,對(duì)車站結(jié)構(gòu)在施工過(guò)程中的受力和變形狀態(tài)進(jìn)行分析,進(jìn)而對(duì)蓋挖逆作法施工中的一些關(guān)鍵受力狀態(tài)進(jìn)行研究。
1.1物理模型
東湖渠站位于廣順北大街下方,北至望京北路,南至利澤中街,為14號(hào)線標(biāo)準(zhǔn)站。東湖渠站為地下兩層,采用12 m島式站臺(tái),雙柱三跨混凝土結(jié)構(gòu)。車站共設(shè)4個(gè)出入口、1個(gè)消防通道和2組風(fēng)亭。車站主體采用蓋挖逆作法施工,頂板分東西半幅施工,附屬結(jié)構(gòu)采用暗挖配合明挖施工。
用Midas GTS有限元軟件模擬車站施工過(guò)程。根據(jù)現(xiàn)場(chǎng)地質(zhì)條件和水文條件確定材料參數(shù),建立車站和周圍土體的三維模型。三維模型140 m長(zhǎng),75 m寬,50 m深。模擬區(qū)域?yàn)榭v向:車站靠北側(cè)中間柱軸20號(hào)至軸36號(hào)122 m外加北端基坑外側(cè)18 m的土體;橫向:基坑內(nèi)側(cè)25 m加上基坑?xùn)|西兩側(cè)25 m的土體;豎向從地面往下50 m深土體。
1.2本構(gòu)模型和參數(shù)選取
概化后的土層參數(shù)見(jiàn)表1,土體單元類型為實(shí)體單元,采用摩爾-庫(kù)倫本構(gòu)模型。中間柱、中間樁單元類型為梁?jiǎn)卧?,頂板、中板、底板、?cè)墻的單元類型為板單元,懸臂樁單元類型為板單元,圍護(hù)樁在建模過(guò)程中通過(guò)修改單元屬性的手段來(lái)添加,單元類型為實(shí)體單元。懸臂樁單元類型為板單元,圍護(hù)樁在建模過(guò)程中通過(guò)修改單元屬性的手段來(lái)添加,單元類型為實(shí)體單元。上述結(jié)構(gòu)材料假定為彈性材料,見(jiàn)表2。
表1 地層土性參數(shù)
表2 車站結(jié)構(gòu)的幾何尺寸和材料參數(shù)
1.3懸臂樁和圍護(hù)樁的等效方法
為簡(jiǎn)化模型,可以按照剛度等效的原則,將密集布置的懸臂樁和圍護(hù)樁轉(zhuǎn)化成一定厚度的地下連續(xù)墻(圖1)。具體方法為:設(shè)鉆孔樁的樁徑為D,樁凈距為t,則單根樁應(yīng)等價(jià)為長(zhǎng)D+t的地下連續(xù)墻,設(shè)等價(jià)后的地下連續(xù)墻厚度為h,則由等剛度轉(zhuǎn)換的原則得:
(1)
于是
(2)
圖1等效剛度轉(zhuǎn)換
對(duì)于攪噴樁,可以視作對(duì)土體灌漿,灌漿后土體剛度E=2E0(E0為灌漿前土體剛度)。由于車站采用蓋挖逆作法施工,需先施作車站頂板,頂板基坑四周懸臂樁采用Φ600@1000鉆孔灌注樁,即D=0.6,t=0.4,得到懸臂樁等效地下連續(xù)墻厚度h=0.42 m。
車站主體圍護(hù)結(jié)構(gòu)為Φ1000@1500鉆孔灌注樁,采用長(zhǎng)螺旋高壓攪噴樁與圍護(hù)樁咬合成墻的方案。攪噴樁按灌漿后土體來(lái)計(jì)算剛度,即
(3)
式中,攪噴樁和圍護(hù)樁直徑、間距均相同D=1,t=0.5,混凝土彈性模量Ec=2.17×107kPa,灌漿前土體彈性模量E0=2×104kPa,則計(jì)算得到圍護(hù)樁與攪噴樁等效地下連續(xù)墻厚度h=0.75 m。
1.4荷載確定
車站施工期間,基坑兩側(cè)存在雙向三車道車輛荷載,城市道路主要車輛為小客車。根據(jù)城市道路工程設(shè)計(jì)規(guī)范(CJJ-37-2012)第3.3.1條:小客車軸距為3.8 m;第3.6.1條:道路路面結(jié)構(gòu)設(shè)計(jì)應(yīng)以雙輪組單軸載100 kN為標(biāo)準(zhǔn)軸載??紤]前后兩輛車之間軸距為10 m,車道寬度為3.3 m,則等效的路面荷載:P=2×100/(3.8+10)×3.3=4.4 kPa。
圖2為中間柱的計(jì)算結(jié)果。圖中位置1位于中柱內(nèi)距頂板1 m處,位置2位于中柱內(nèi)距頂板2.5 m處,位置3位于中柱內(nèi)距頂板4 m處,位置4位于中柱內(nèi)距中板1 m處,位置5位于中柱內(nèi)距中板3.5 m處,位置6位于中柱內(nèi)距中板6 m處??梢钥闯?,在西半幅中柱施做后一段時(shí)間內(nèi)由于西側(cè)上覆土開(kāi)挖,導(dǎo)致西中柱周圍土體回彈,對(duì)中柱產(chǎn)生向上的負(fù)摩阻力,使柱體處于受拉狀態(tài)。頂板施工后,西側(cè)上覆土回填,中柱在承受上覆土荷載的情況下逐漸變?yōu)槭軌籂顟B(tài),并且壓應(yīng)力隨著施工過(guò)程逐漸增大。數(shù)值模擬最大拉應(yīng)力為1.9 MPa,最大壓應(yīng)力為8.3 MPa。
圖2中間柱應(yīng)力隨施工步驟變化過(guò)程
針對(duì)西半幅頂板,由西向東沿頂板橫向設(shè)置5個(gè)計(jì)算點(diǎn),分別距離西側(cè)墻1 m,3 m,5 m,7 m,9 m處。由圖3可以看出,頂板上側(cè)在施工過(guò)程中位置2、3、5處于受壓狀態(tài),即側(cè)墻與中間柱之間中間位置和中間柱與中間柱之間中間位置處于受壓狀態(tài),最大壓應(yīng)力值4 MPa。位置1處于頂板與側(cè)墻連接處附近,在施工過(guò)程中一直處于受拉狀態(tài),最大拉應(yīng)力4 MPa。位置4由于處在中間柱附近,豎向位移受到約束,施工過(guò)程中應(yīng)力狀態(tài)反轉(zhuǎn),由受壓狀態(tài)轉(zhuǎn)為受拉狀態(tài),拉應(yīng)力值較小。頂板下側(cè)在施工過(guò)程中位置2、3、5則處于受拉狀態(tài),最大值3 MPa;此時(shí)位置1則處于受壓狀態(tài),最大壓應(yīng)力值5 MPa;同樣的位于中間柱附近的位置4由受拉狀態(tài)轉(zhuǎn)為受壓狀態(tài),壓應(yīng)力值也較小。
圖3表明,頂板上側(cè)在靠近側(cè)墻和中間柱附近位置處于受拉狀態(tài),其余位置處于受壓狀態(tài);而頂板下側(cè)則相反,在靠近側(cè)墻和中間柱附近位置處于受壓狀態(tài)。應(yīng)該看到,數(shù)值分析中頂板上、下的最大拉應(yīng)力值都較大,說(shuō)明頂板的這些位置是受力較大的地方,需要重點(diǎn)注意。
圖3頂板應(yīng)力隨施工步驟變化過(guò)程
圖4表示的是數(shù)值模擬中板上側(cè)和中板下側(cè)應(yīng)力計(jì)算結(jié)果。位置4~7處于中間柱附近,在施工過(guò)程前期上側(cè)處于受拉狀態(tài),拉應(yīng)力隨著施工過(guò)程逐漸增大,最大值達(dá)到1.6 MPa,在地下二層土開(kāi)挖后,相應(yīng)位置處中板上側(cè)應(yīng)力發(fā)生轉(zhuǎn)變,與中板兩側(cè)其他位置處一樣,處于受壓狀態(tài),最大壓應(yīng)力2 MPa。中板兩側(cè)位置在施工過(guò)程前期下側(cè)處于受拉狀態(tài),拉應(yīng)力隨著施工過(guò)程逐漸增大,最大值達(dá)0.6 MPa,在地下二層土開(kāi)挖后,相應(yīng)位置處中板下側(cè)應(yīng)力發(fā)生轉(zhuǎn)變,與中板中間位置處一樣,處于受壓狀態(tài),最大壓應(yīng)力2.3 MPa。
圖4中板應(yīng)力隨施工步驟變化過(guò)程
由圖4可以看出,中板上側(cè)靠近中間柱位置及中板下側(cè)靠近側(cè)墻位置處應(yīng)力均會(huì)隨施工過(guò)程而發(fā)生應(yīng)力轉(zhuǎn)變,由最初的受拉狀態(tài)轉(zhuǎn)為受壓狀態(tài),最終中板各個(gè)位置上、下側(cè)均處于受壓狀態(tài)。導(dǎo)致中板應(yīng)力反轉(zhuǎn)的原因是由于地下二層土開(kāi)挖后,基坑外側(cè)土的側(cè)向土壓力釋放,原來(lái)由地下二層土承擔(dān)的水平應(yīng)力在地下二層土挖走后,通過(guò)側(cè)墻轉(zhuǎn)移到中板承擔(dān)。豎直方向中板主要承擔(dān)自重荷載,而自重荷載產(chǎn)生的彎矩相對(duì)較小,截面上的拉應(yīng)力和壓應(yīng)力大大小于基坑外側(cè)土側(cè)向土壓力轉(zhuǎn)移到中板上產(chǎn)生的壓應(yīng)力,二者疊加后使得整個(gè)中板處于受壓狀態(tài)。
對(duì)側(cè)墻表面沿豎向方向依次設(shè)置3個(gè)計(jì)算點(diǎn),位置由上至下分別為距離頂板1 m處、2.5 m、4 m處。圖5表明,側(cè)墻內(nèi)側(cè)施工后處于受拉狀態(tài),且隨施工過(guò)程拉應(yīng)力逐漸增大,在側(cè)墻下部的地下二層土開(kāi)挖后,應(yīng)力隨即變?yōu)槭軌籂顟B(tài),最大壓應(yīng)力0.45 MPa。側(cè)墻外側(cè)施工后處于受壓狀態(tài),在地下二層土開(kāi)挖后,壓應(yīng)力迅速增大,最大值達(dá)到0.23 MPa。另外,側(cè)墻內(nèi)側(cè)位置3處以及側(cè)墻外側(cè)位置2處在地下二層土開(kāi)挖后均處于受拉狀態(tài),最大拉應(yīng)力在0.05 MPa左右。
圖5側(cè)墻應(yīng)力隨施工步驟變化過(guò)程
此外,由圖5(a)可看出,在施工過(guò)程中靠上側(cè)的側(cè)墻內(nèi)表面施做后處于受拉狀態(tài),在地下二層土開(kāi)挖后轉(zhuǎn)為受壓狀態(tài)。圖6表示數(shù)值模擬的中板上表面應(yīng)力變化情況,顯示在施工過(guò)程中中板上表面大部分處于受壓狀態(tài),只有在中柱和側(cè)墻附近可能出現(xiàn)局部受拉的情況。
圖6中板上表面應(yīng)力隨施工步驟變化過(guò)程
蓋挖逆作法中樁(柱)、連續(xù)墻在施工時(shí)間上存在先后順序,故而中柱與連續(xù)墻之間存在差異沉降,當(dāng)這種差異沉降過(guò)大,超過(guò)規(guī)定值時(shí),往往會(huì)導(dǎo)致結(jié)構(gòu)開(kāi)裂,甚至產(chǎn)生破壞。
現(xiàn)分析車站結(jié)構(gòu)東西側(cè)墻和東西中柱的沉降,其隨施工步驟的沉降值見(jiàn)表3。東、西中柱的沉降值(或回彈)在相應(yīng)上覆土層開(kāi)挖時(shí)較大,隨著上覆土回填及后續(xù)施工,回彈值逐漸減小,在地下二層土開(kāi)挖時(shí)達(dá)到最大。兩側(cè)側(cè)墻在地下二層土開(kāi)挖時(shí)回彈最大,隨后減小。由于施工的先后差異,東、西中柱,東、西側(cè)墻之間存在差異沉降。例如,東、西中柱間最大差異沉降發(fā)生在西上覆土開(kāi)挖后,達(dá)到8 mm;西側(cè)墻與西中柱、東側(cè)墻與東中柱最大差異沉降均發(fā)生在底板側(cè)墻施工階段,分別為3.1 mm和0.5 mm。
表3 側(cè)墻、中柱沉降隨施工步驟變化過(guò)程 單位:mm
根據(jù)以往的施工經(jīng)驗(yàn),地下結(jié)構(gòu)中相鄰中柱之間以及相鄰的樁柱和連續(xù)墻之間的沉降差值不能超過(guò)20 mm,或者不能超過(guò)柱間距長(zhǎng)度的1/400。在本研究中,側(cè)墻與相鄰中柱之間的距離為7 m,東、西中柱間距離為6 m,柱間距長(zhǎng)度的1/400即15 mm,側(cè)墻與中柱間最大差異沉降3.1 mm,小于20 mm;東、西中柱間最大差異沉降8 mm,小于柱間距長(zhǎng)度1/400(圖7)。因此,施工過(guò)程中側(cè)墻、中柱的差異沉降能滿足要求。
圖7側(cè)墻與中柱差異沉降隨施工過(guò)程變化
在基坑開(kāi)挖過(guò)程中,圍護(hù)結(jié)構(gòu)的水平位移一直是工程關(guān)注的一個(gè)重點(diǎn),當(dāng)其水平位移過(guò)大時(shí),往往預(yù)示著基坑有破壞的可能。在本車站中采用蓋挖逆作法施工,開(kāi)挖車站上覆土?xí)r,基坑外側(cè)土靠基坑外圍懸臂樁支持。在上覆土開(kāi)挖和頂板施工過(guò)程中,基坑?xùn)|側(cè)懸臂樁樁身的水平位移如圖8,在步驟3東側(cè)上覆土開(kāi)挖時(shí),樁體水平位移變化較大,樁頂位移達(dá)到4 mm,在上覆土回填前,樁頂位移最大6.6 mm。規(guī)范中規(guī)定基坑側(cè)壁水平位移小于30 mm或者2‰H(H為基坑開(kāi)挖深度),本車站中上覆土層厚5 m,即H=5 000 mm,2‰H=10 mm??梢?jiàn),懸臂樁最大水平位移小于規(guī)范允許值。
圖9為東圍護(hù)樁隨施工過(guò)程水平位移變化情況??梢钥闯觯诘叵乱粚油灵_(kāi)挖過(guò)程中,樁身水平位移顯著增大,由于樁頂受已施工的頂板約束,水平位移幾乎不發(fā)展。地下二層土開(kāi)挖過(guò)程中,相應(yīng)位置樁體水平位移發(fā)展最快,并超過(guò)了樁頂?shù)乃轿灰?,而圍護(hù)樁上部由于受到頂板和中板的約束,水平位移發(fā)展很緩慢。整個(gè)樁體的變形趨勢(shì)也由原來(lái)的懸臂式變形曲線逐漸變?yōu)镾形變形曲線,在圖9中位置8處,位移曲線存在反彎點(diǎn),實(shí)際工程中此處對(duì)應(yīng)中板位置。由此可見(jiàn),雖然圍護(hù)樁埋深較懸臂樁深,相應(yīng)的側(cè)向土壓力更大,但由于頂板和中板的側(cè)向支撐約束作用,圍護(hù)樁在施工過(guò)程中的水平位移也較小。東側(cè)圍護(hù)樁在整個(gè)施工過(guò)程中最大水平位移6.3 mm,樁頂最大水平位移4.7 mm,同樣均小于規(guī)范允許值。這也反映出地鐵車站蓋挖逆作法施工的一個(gè)優(yōu)勢(shì),那就是施工過(guò)程中圍護(hù)結(jié)構(gòu)變形小,能夠有效控制周圍土體的變形和地表沉降,有利于保護(hù)臨近建筑物和構(gòu)筑物。
圖8東懸臂樁水平位移隨施工過(guò)程變化
圖9東圍護(hù)樁水平位移隨施工過(guò)程變化
(1)施工過(guò)程中,中柱的應(yīng)力變化情況為先處于受拉狀態(tài),在上覆土回填之后,轉(zhuǎn)為受壓狀態(tài),并且壓應(yīng)力逐漸增大。
(2)頂板上側(cè)在靠近側(cè)墻和中間柱附近位置處于受拉狀態(tài),其余位置處于受壓狀態(tài);而頂板下側(cè)則相反,在靠近側(cè)墻和中間柱附近位置處于受壓狀態(tài),其余位置處于受拉狀態(tài)。
(3)中板上側(cè)靠近中間柱位置以及中板下側(cè)靠近側(cè)墻位置處應(yīng)力均會(huì)隨著施工過(guò)程而發(fā)生應(yīng)力轉(zhuǎn)變,由最初的受拉狀態(tài)轉(zhuǎn)為受壓狀態(tài),最終中板各個(gè)位置上、下側(cè)均處于受壓狀態(tài)。
(4)側(cè)墻受力較為復(fù)雜,在側(cè)墻對(duì)應(yīng)位置的地下二層土未開(kāi)挖前,側(cè)墻內(nèi)側(cè)均受拉,外側(cè)均受壓;在地下二層土開(kāi)挖后,側(cè)墻內(nèi)外側(cè)大部分處于受壓狀態(tài),只在側(cè)墻下部存在局部受拉區(qū)域,拉應(yīng)力值相對(duì)較小。
(5)中柱間最大差異沉降8 mm,中間柱、側(cè)墻最大差異沉降3.1 mm;基坑最大水平位移分別為6.6 mm和6.3 mm,均小于規(guī)范允許值,施工過(guò)程有效的控制了周圍地層的變形。
參考文獻(xiàn):
[1]白冰,聶慶科,吳剛,等.考慮空間效應(yīng)的深基坑雙排樁支護(hù)結(jié)構(gòu)計(jì)算模型[J].建筑結(jié)構(gòu)學(xué)報(bào),2010,31(8):118-124.
[2]鄭剛,劉慶晨,鄧旭,等.天津站地下?lián)Q乘中心基坑工程2標(biāo)段蓋挖逆作法的實(shí)測(cè)分析[J].天津大學(xué)學(xué)報(bào),2012,45(10):930-937.
[3]趙玉璽.蓋挖逆作法地鐵車站方形中間樁柱定位技術(shù)[J].隧道建設(shè),2013,33(7):607-614.
[4]丁克勝,李旭,易順建.仿真蓋挖逆作法土方開(kāi)挖過(guò)程中地下連續(xù)墻的變形分析[J].建筑技術(shù),2014,45(3):204-206.
[5]蘇潔,張頂立,高自友,等.蓋挖逆作法施工地鐵車站結(jié)構(gòu)變形及其控制[J].中國(guó)鐵道科學(xué),2010,31(1):59-65.
[6]魏玉超,張志軍.城市地鐵車站工程的蓋挖逆作法施工技術(shù)[J].城市建設(shè)理論研究,2013,(24):27-34.
[7]鄭騏.蓋挖逆作法地鐵車站的中間樁施工技術(shù)[J].市政技術(shù),2009,27(5):498-500.
[8]趙育紅.地鐵車站蓋挖逆作法中間柱施工技術(shù)[J].四川建筑科學(xué)研究,2014,(3):327-330.
TheNumericalSimulationoftheStressandDeformationontheStructuresofSubwayStationsInducedbyTop-downExcavationMethod
LIU Guo-hui1,XIAO Feng2,LI Ying-jie1
(1.ChinaConstructionCommunicationsEngineeringGroupCorp.,Ltd.,Beijing100142,China;2.SchoolofCivilEngineering,BeijingJiaotongUniversity,Beijing100044,China)
Abstract:In the construction process,the stress and deformation characteristics of the main structure induced by top-down excavation method of a subway station in Beijing were studied by adopting numerical caculation.By using finite element numerical simulation,the stress state and deformation status of the top,middle and bottom plate and side walls in response to the construction process of the station were studied to analyze the changing process of the middle pile(or column)under stress during the construction.The results showed that the stress value of the middle pile changed continuously with the construction of the subway station,which was tensile stress at first and then turned into compressive stress when the pit was backfilled;whereas the the upper part of the top plate which was adjacent to the side walls and the middle pile was under tensile stress,but the bottom part was the opposite.Besides,the differential settlement of the middle pile and enclosing supporting structure were studied,which proved that the maximum differential settlement of the middle pile and that between the middle pile and the sidewalls were all less than the allowable value according to the specification.In the end,the safety of the construction of the subway station and its impact on the adjacent buildings were evaluated.
Keywords:top-down excavation method;numerical calculation;internal stress of the structure;settlement of soil layer
DOI:10.3969/j.issn.1672-1144.2014.06.039
中圖分類號(hào):TU473
文獻(xiàn)標(biāo)識(shí)碼:A
文章編號(hào):1672—1144(2014)06—0193—06
作者簡(jiǎn)介:劉國(guó)輝(1977—),女,河北承德人,高級(jí)工程師,主要從事巖土工程設(shè)計(jì)與施工方面的研究工作。
收稿日期:2014-08-11修稿日期:2014-09-07