劉小麗,馬 悅,郭冠群,陶 韜,周 賀
(中國海洋大學(xué)1.海洋環(huán)境與生態(tài)教育部重點實驗室;2.工程學(xué)院,山東青島266100)
PLAXIS2D模擬計算基坑開挖工程的適用性分析*
劉小麗1,馬 悅2,郭冠群1,陶 韜2,周 賀1
(中國海洋大學(xué)1.海洋環(huán)境與生態(tài)教育部重點實驗室;2.工程學(xué)院,山東青島266100)
巖土工程有限元分析軟件PLAXIS2D已經(jīng)被廣泛應(yīng)用于基坑開挖的數(shù)值模擬計算中,但關(guān)于其模擬計算基坑開挖工程的適用性問題,目前還未見有相關(guān)文獻進行深入分析。鑒于此,利用PLAXIS2D對3個實際深基坑開挖工程進行了數(shù)值模擬計算,通過計算結(jié)果與實測結(jié)果的對比,分析了PLAXIS2D在基坑開挖工程中的適用性;在此基礎(chǔ)上,對應(yīng)用PLAXIS2D進行基坑開挖數(shù)值模擬計算的適用條件進行了探討,研究結(jié)論對正確應(yīng)用PLAXIS2D進行基坑開挖工程的模擬計算具有較好的參考價值。
PLAXIS2D;有限元;基坑開挖;數(shù)值模擬;適用性
巖土工程有限元分析軟件PLAXIS的研制始于1987年,初始目的是為了進行建立在軟土上的河堤分析。此后,PLAXIS一直不斷發(fā)展,直到今天,已經(jīng)成為功能強大的專門針對巖土工程變形與穩(wěn)定計算的有限元分析軟件,二維的PLAXIS有限元程序(PLAXIS2D)能夠計算平面應(yīng)變和軸對稱兩類工程問題[1]。
許多文獻都用巖土工程有限元分析軟件PLAXIS2D對基坑開挖工程進行了計算分析[2-7],但一般都是直接應(yīng)用該軟件對所研究的問題進行計算,沒有針對PLAXIS2D本身在基坑開挖工程中的適用性問題進行較為詳細的分析。鑒于此,利用PLAXIS2D對幾個具體的深基坑開挖工程實例進行計算,并將計算結(jié)果與工程實測結(jié)果進行對比,分析探討PLAXIS2D在基坑開挖工程中的適用性及適用條件,對正確應(yīng)用PLAXIS2D進行基坑開挖的數(shù)值模擬計算提供參考和建議。
上海軌道交通10號線某車站基坑工程,基坑寬度約17.8 m,標準段開挖深度約14.6 m,基坑支護結(jié)構(gòu)采用700 mm厚的地下連續(xù)墻,墻深28.0 m。地下連續(xù)墻采用C30混凝土,容重25 k N/m3,彈性模量E=30 GPa,泊松比v=0.2。設(shè)置4道內(nèi)支撐,第1道為鋼筋混凝土內(nèi)支撐,位于樁頂,混凝土標號C30,文中未提供該支撐截面尺寸,參考同類基坑設(shè)計資料,混凝土內(nèi)支撐構(gòu)件截面寬1 000 mm,高800 mm,支撐間距8.0 m(在常用混凝土支撐截面尺寸范圍內(nèi),對計算結(jié)果無實質(zhì)性影響);其余均為鋼管內(nèi)支撐,采用609 mm,厚度t=16 mm的鋼管,分別位于距樁頂5.0、8.7和11.7 m的位置,支撐水平間距6.0 m。基坑內(nèi)開挖面以下4 m深、3 m寬范圍內(nèi)土體采用高壓旋噴樁進行加固。
由上至下基坑地層分布及其基本計算參數(shù)如表1所示,其中的內(nèi)摩擦角和粘聚力為固結(jié)不排水指標(以下各例同)??觾?nèi)土體加固,根據(jù)經(jīng)驗,相應(yīng)加固區(qū)土體的力學(xué)性質(zhì)提高約2.5倍(以下各例同)。
表1 實例1的土層計算參數(shù)Table 1 Soil parameters of case 1
上海麗晶苑大廈基坑工程,基坑寬度42.0 m,基坑開挖深度12.65 m。該基坑工程支護結(jié)構(gòu)采用直徑1 000 mm,樁凈間距150 mm,長27 m的鉆孔灌注樁擋土,水泥攪拌樁止水。設(shè)置3道混凝土水平內(nèi)支撐,支撐構(gòu)件截面寬1 000 mm,高800 mm。3道水平內(nèi)支撐分別位于距樁頂1.45、5.25和9.20 m處?;觾?nèi)側(cè)開挖面以下4 m深度、6 m寬度范圍內(nèi)的土體進行了注漿加固。
由上至下基坑地層分布及其基本計算參數(shù)如表2中所示。
表2 實例2的土層計算參數(shù)Table 2 Soil parameters of case 2
上海徐家匯基坑工程,基坑面積約9 000 m2,開挖深度13.4 m?;訃o結(jié)構(gòu)采用地下連續(xù)墻及鋼筋混凝土平面框架內(nèi)支撐體系。地下連續(xù)墻厚0.8 m,墻深約26.0 m,3道鋼筋混凝土內(nèi)支撐構(gòu)件截面尺寸分別為800 mm×800 mm,1 200 mm×1 000 mm,1 200 mm×1 000 mm,分別位于地面下1.6、6.0和10.0 m處。連續(xù)墻內(nèi)側(cè)開挖面以下3 m深度、5 m寬度范圍內(nèi)土體采用水泥攪拌樁進行了加固。
由上至下基坑地層分布及其基本計算參數(shù)如表3所示。
表3 實例3的土層計算參數(shù)Table 3 Soil parameters of case 3
基坑沿寬度方向幾何對稱,因此在數(shù)值計算時只選取基坑寬度的1/2進行分析。建立平面應(yīng)變有限元模型,巖土體采用15節(jié)點三角形平面單元;土體材料模型采用硬化土模型(Hardening-soil模型);支護樁和內(nèi)支撐采用線彈性模型,根據(jù)剛度等效原則,支護樁和地下連續(xù)墻用等效厚度的板來模擬,用錨錠桿來模擬內(nèi)支撐;通過界面強度折減因子Rinter對相應(yīng)土體的粘聚力和內(nèi)摩擦角進行折減,以此來考慮結(jié)構(gòu)與相鄰?fù)馏w之間的粘聚力和內(nèi)摩擦角[1]。模型的底部邊界施加完全固定約束,左側(cè)邊界施加水平向約束,右側(cè)邊界施加對稱邊界條件;計算過程中不考慮基坑地下水位。具體數(shù)值計算模型見圖1。
圖1 實例1基坑的PLAXIS2D計算模型(1/2基坑寬度)Fig.1 PLAXIS2D model for case 1(half of the width)
由于計算模型為平面應(yīng)變模型,因此選擇接近基坑長邊中部位置剖面的相應(yīng)實測值進行比較。
實例1開挖到基坑底面時,支護結(jié)構(gòu)水平位移和基坑周邊地表沉降計算值和實測值的對比如圖2所示。
圖2 實例1數(shù)值計算結(jié)果與實測結(jié)果對比Fig.2 Comparison of numerical computing results and the monitored ones for case 1
由圖2(a)可知,計算的支護結(jié)構(gòu)水平位移最大值為43 mm,發(fā)生在距墻頂16.0 m(基坑開挖面下約1.4 m)的位置;實測支護結(jié)構(gòu)水平位移最大值約為39~41 mm,最大值位置在開挖面附近。此例中計算得到的地下連續(xù)墻水平位移最大值的位置較實測位置偏下,支護樁水平位移分布模式及其數(shù)值與實測結(jié)果比較符合。
根據(jù)圖2(b)可知,計算得到的基坑周邊地表沉降最大值為18 mm,開挖至基底的實測地表最大沉降值約16~25 mm,計算結(jié)果在實測值范圍之內(nèi)。上述數(shù)據(jù)對比分析表明,Plaxis2D計算結(jié)果具有較好的預(yù)測性。
開挖到基坑底面時,實例2支護結(jié)構(gòu)水平位移和基坑周邊地表沉降計算值和實測值的對比如圖3所示。
由圖3(a)可知,支護結(jié)構(gòu)的計算水平位移最大值43 mm,最大值位置距樁頂約14.7 m(基坑開挖面下約2 m);實測支護結(jié)構(gòu)水平位移最大值為37 mm,距樁頂約13 m。PLAXIS2D的計算結(jié)果稍偏大,支護樁的水平位移分布模式與實測值一致。
圖3 實例2數(shù)值計算結(jié)果與實測結(jié)果對比Fig.3 Comparison of numerical computing results and the measured ones for case 2
圖3(b)是PLAXIS2D計算的基坑周邊地表沉降分布曲線,沉降最大值為24mm,由于無基坑周邊地表沉降實測數(shù)據(jù),故未將計算結(jié)果和實測結(jié)果進行比較,但可以看出,計算得到的地表沉降分布模式符合一般規(guī)律。對于本例,PLAXIS2D的計算結(jié)果稍偏于安全,能夠較好的預(yù)測實際基坑開挖工程的變形。
對于實例3,分別對基坑底部土體加固和不加固2種情況進行了模擬計算,開挖到基坑底面時,實例3支護結(jié)構(gòu)水平位移和基坑周邊地表沉降計算值和實測值的對比如圖4所示。圖4中圖例顯示的“加固”表示基坑底部土體加固條件下的計算結(jié)果,“不加固”表示基坑底部土體不加固條件下的計算結(jié)果。圖4中鑒于距基坑邊0~20 m范圍內(nèi)的現(xiàn)場沉降監(jiān)測點數(shù)據(jù)較少,因此未把實測值數(shù)據(jù)點連線,以避免形成不正確的地表沉降模式。
由圖4(a)可知,基底土體加固條件下支護結(jié)構(gòu)水平位移計算最大值為89 mm,距墻頂約14.9 m(基坑開挖面下約1.5 m),不考慮基坑土體加固條件下(即不加固)支護結(jié)構(gòu)水平位移計算最大值為98 mm,距墻頂約14.9 m;實測支護結(jié)構(gòu)水平位移最大值為103 mm,發(fā)生在距墻頂約17.8 m的位置(基坑開挖面下約4 m)。
如圖4(b),加固條件下地表沉降計算最大值為52 mm,不加固條件下地表沉降計算最大值為61 mm,實測數(shù)據(jù)中的地表沉降值最大值約為60 mm。
由圖4可以看出,支護結(jié)構(gòu)水平位移和基坑周邊地表沉降的計算結(jié)果與實測結(jié)果相比,其分布模式一致性較好,加固條件下計算結(jié)果偏小,不加固條件下計算結(jié)果與實測結(jié)果符合較好,其原因可能是該基坑加固施工效果不理想,加固區(qū)土體沒有起到相應(yīng)的作用,或者是模型計算參數(shù)取值問題導(dǎo)致加固條件下PLAXIS2D計算結(jié)果偏小,亦或是其他不確定性因素如基坑邊超載、施工超挖等??傮w而言,PLAXIS2D的計算結(jié)果基本能反映基坑的實際變形情況,其計算結(jié)果能夠?qū)嶋H基坑工程的設(shè)計和施工起到較好的指導(dǎo)作用。
圖4 實例3數(shù)值計算結(jié)果與實測結(jié)果對比Fig.4 Comparison of numerical computing results and the monitoved ones for case 3
由于無相應(yīng)支護結(jié)構(gòu)內(nèi)力監(jiān)測資料,因此無法對PLAXIS2D計算的支護結(jié)構(gòu)內(nèi)力進行分析比較,但通過以上3個工程實例中支護結(jié)構(gòu)水平位移與基坑周邊地表沉降的計算結(jié)果與相應(yīng)實測結(jié)果的對比分析,可以看出,PLAXIS2D計算得到的支護結(jié)構(gòu)水平位移和基坑周邊地表沉降的數(shù)值和分布模式與實測結(jié)果具有較好的一致性,PLAXIS2D計算的支護結(jié)構(gòu)水平位移最大值位置一般位于開挖面以下1~2 m,在實測最大值位置上下約2 m的范圍內(nèi)。上述結(jié)果表明,PLAXIS2D在基坑開挖工程的模擬計算中具有較好的適用性,其模擬計算結(jié)果能夠?qū)嶋H基坑工程變形進行較好的預(yù)測,進而可以有效的指導(dǎo)實際基坑工程的設(shè)計與施工。
在實際的基坑開挖工程中,一般在基坑周邊設(shè)置止水帷幕后,進行基坑內(nèi)降水,而基坑外的水位則基本保持不變或有少量下降,因此,在基坑開挖過程中,基坑內(nèi)側(cè)開挖范圍內(nèi)無地下水存在,而基坑外側(cè)的地下水是存在的。
當基坑模型設(shè)置地下水位時,PLAXIS2D的計算過程中即單獨考慮了水壓力作用,將水土分算得到的水土壓力作為側(cè)向壓力計算支護結(jié)構(gòu)的內(nèi)力變形;當不考慮地下水位,即不設(shè)水位,并將地下水位以下的土體容重設(shè)置為飽和容重時,PLAXIS2D計算過程中相當于將水土合算得到的水土壓力作為側(cè)向壓力計算支護結(jié)構(gòu)的內(nèi)力變形。
對上述3個工程實例進行計算可知,當不設(shè)置地下水位時,計算結(jié)果與實測結(jié)果一致性較好;反之,當基坑設(shè)置相應(yīng)的地下水位時,PLAXIS2D計算得到的支護結(jié)構(gòu)水平位移和基坑周邊地表沉降均遠大于實測值,與實際情況嚴重不符(此處限于篇幅,未給出相應(yīng)的計算結(jié)果)。因此,對于軟土地區(qū)的基坑開挖工程(基坑土層絕大部分為滲透性差的黏土層或粘質(zhì)粉土層),采用PLAXIS2D進行基坑開挖模擬計算時,宜不設(shè)置基坑地下水位,而將地下水位以下的土層容重設(shè)置為飽和容重,即相當于按照水土合算計算支護結(jié)構(gòu)上的側(cè)壓力。
在用PLAXIS2D計算以滲透性較好的砂層為主的基坑開挖工程時,原則上應(yīng)設(shè)置相應(yīng)的基坑地下水位,即基于水土分算原理計算支護結(jié)構(gòu)上的側(cè)壓力,但應(yīng)多注重計算經(jīng)驗的積累,當計算結(jié)果較實測結(jié)果大較多時,在排除其他因素后宜考慮將水壓力進行折減,具體的折減幅度應(yīng)根據(jù)類似工程經(jīng)驗確定。
在利用PLAXIS2D計算兼具以黏土層和砂層為主的基坑開挖工程時,根據(jù)一般計算經(jīng)驗,計算時可不設(shè)置基坑外地下水位,將地下水位以下的土層容重設(shè)置為飽和容重進行計算即可。具體怎樣計算更加符合基坑工程實際,還有待積累更多實際工程計算經(jīng)驗。
關(guān)于水土壓力的分算與合算問題,一直以來在工程界都存在著異議[12]。雖然水土分算的概念清晰,理論嚴密,然而由于基坑開挖問題的復(fù)雜性,其應(yīng)用于基坑支護結(jié)構(gòu)上側(cè)壓力的計算時,其適用性并不好[13-14],其中的深層原因還有待于探討。本文不討論水土分算或水土合算在理論上的合理性,而是從工程實際出發(fā),根據(jù)計算結(jié)果與實測結(jié)果的對比分析,討論PLAXIS2D模擬基坑開挖工程的適用性條件,因此計算經(jīng)驗及實際工程經(jīng)驗的積累也就顯得尤為重要。
在PLAXIS2D中包含多種土體材料模型,其中的硬化土模型(HS模型)和摩爾-庫倫理想彈塑性模型(MC模型)應(yīng)用較多。摩爾-庫倫理想彈塑性模型主要包含彈性模量E,泊松比υ,粘聚力c,內(nèi)摩擦角和剪脹角ψ,這些參數(shù)對于大多數(shù)巖土工程師來說都是熟悉的,并且容易獲得[1]。硬化土模型是1種改進了的模擬巖土行為的模型,涉及到的土模量參數(shù)主要包含標準固結(jié)排水三軸試驗中的割線模量Eref50,主固結(jié)儀加載中的切線模量Erefoed,卸載-重新加載的剛度Erefur。模型中的參數(shù)m表示剛度應(yīng)力水平相關(guān)冪指數(shù),取值在0.5~1.0之間,土質(zhì)越軟,m越接近1.0。
以實例1的計算結(jié)果為例,分析比較2種土體模型在基坑開挖工程中的適用性。有限元模型如前所述,土體模型分別采用摩爾-庫倫理想彈塑性模型和硬化土模型,開挖至基坑底面時2種模型條件下的計算結(jié)果如圖5所示,內(nèi)支撐軸力計算結(jié)果如表4所示。
由圖5(a)可知,由硬化土模型(HS模型)計算出支護結(jié)構(gòu)水平位移最大值為43 mm,發(fā)生在距墻頂16.0 m的位置(基坑開挖面下約1.4 m);摩爾-庫侖模型(MC模型)計算出支護結(jié)構(gòu)水平位移最大值為41 mm,發(fā)生在距墻頂17.3 m的位置(在基坑開挖面下約2.7 m);實測支護結(jié)構(gòu)水平位移最大值約為39~41 mm,最大值位置在開挖面附近。總體而言,由摩爾-庫侖模型和硬化土模型計算得到的支護結(jié)構(gòu)水平位移數(shù)值與分布模式均與實測值符合較好,只是摩爾-庫侖模型計算得到的支護結(jié)構(gòu)水平位移在連續(xù)墻底部較實測值偏大。
圖5 實例1兩種數(shù)值計算結(jié)果與實測結(jié)果的對比Fig.5 Comparison of two kinds of numerical model results and the monitored ones for case 1
根據(jù)圖5(b),開挖至基底時的實測地表沉降最大值約16~25 mm,硬化土模型計算得到的基坑周邊地表沉降最大值為18 mm,其分布模式和數(shù)值與實測結(jié)果具有一致性;摩爾-庫侖模型計算得到沉降量最大值為6 mm,且距坑邊20 m范圍內(nèi)土體有過大隆起,該結(jié)果與實測值相差較大,表明摩爾-庫倫模型計算得到的基坑周邊地表沉降值偏小,不能對實際結(jié)果起到較好的預(yù)測作用。
圖5的(c)和(d)分別為使用2種土體模型得到的連續(xù)墻的剪力和彎矩對比圖。從圖中可以看出,利用2種土體模型計算得到的連續(xù)墻的內(nèi)力分布模式一致,剪力大小也較一致,利用摩爾-庫倫模型得到的連續(xù)墻最大彎矩小于硬化土模型的相應(yīng)結(jié)果。
由表4可知,2種土體模型計算得到的內(nèi)支撐軸力數(shù)值基本一致,硬化土模型的計算結(jié)果略大于摩爾-庫倫模型的相應(yīng)計算結(jié)果。
綜上可知,對于基坑開挖工程,硬化土模型的計算結(jié)果優(yōu)于摩爾-庫倫模型的相應(yīng)計算結(jié)果,尤其是基坑周邊的地表沉降計算值,摩爾-庫倫模型的計算結(jié)果明顯偏小。以上計算結(jié)果與該2種土體模型的理論基礎(chǔ)是相對應(yīng)的,硬化土模型中考慮了基坑土體的卸荷回彈特性,能更好的反映土體卸荷性質(zhì)的影響,因此,相對于摩爾-庫倫模型而言,硬化土模型在基坑開挖工程的模擬計算中具有更好的適用性。
利用PLAXIS2D模擬軟土地區(qū)基坑開挖的計算結(jié)果,對硬化土模型的參數(shù)進行敏感性分析可知,m和Rinter對計算結(jié)果的敏感性較?。粌?nèi)摩擦角和對支護結(jié)構(gòu)變形較為敏感,內(nèi)摩擦角和粘聚力對支護結(jié)構(gòu)的內(nèi)力較為敏感,同時這3個參數(shù)的變化范圍也較其他參數(shù)更大。因此,為了保證PLAXIS2D在基坑開挖工程中具有更好的適用性,應(yīng)對硬化土模型的參數(shù)尤其是內(nèi)摩擦角和粘聚力進行詳細的試驗研究,并注重積累相關(guān)地區(qū)經(jīng)驗。
利用巖土工程有限元分析軟件PLAXIS2D對3個基坑開挖工程實例進行了計算,將計算值與實測值進行了對比分析,并探討了PLAXIS2D應(yīng)用于基坑開挖工程的適用性條件,得到了以下主要結(jié)論。
(1)3個基坑開挖工程實例的計算結(jié)果表明,PLAXIS2D在模擬計算基坑開挖工程方面具有較好的適用性,其計算結(jié)果可用于對實際工程結(jié)果的預(yù)測。
(2)對于軟土地區(qū),利用PLAXIS2D模擬基坑開挖工程時,宜不設(shè)置基坑地下水位;對于其他地區(qū)的基坑開挖工程,應(yīng)根據(jù)具體條件進行考慮。
(3)利用PLAXIS2D模擬基坑開挖工程時,對于土體材料模型,硬化土模型較摩爾-庫倫理想彈塑性模型具有更好的適用性。
(4)硬化土模型中,內(nèi)摩擦角、Eref50和粘聚力對計算結(jié)果較為敏感,應(yīng)重視對模型參數(shù)的取值研究和相關(guān)地區(qū)經(jīng)驗的積累。
[1] Brinkgreve R B J.PLAXIS Version 8 Reference Manual[M].Netherlands:A.A.Balkema Publishers,2002.
[2] 張如林,徐奴文.基于PLAXIS的深基坑支護設(shè)計的數(shù)值模擬[J].結(jié)構(gòu)工程師,2010,26(2):131-136.
[3] 付先進,林作忠.基于PLAXIS的超深基坑開挖彈塑性有限元數(shù)值計算與分析[J].中國水運,2010,10(4):190-192.
[4] 林剛,徐長節(jié),蔡袁強.不平衡堆載作用下深基坑開挖支護結(jié)構(gòu)性狀研究[J].巖土力學(xué),2010,31(8):2592-2597.
[5] 李斯海,張玉軍.深基坑開挖與支護過程的平面有限元模擬[J].巖石力學(xué)與工程學(xué)報,1999,18(3):342-345.
[6] 蔡袁強,李碧青,徐長節(jié).挖深不同情況下基坑支護結(jié)構(gòu)性狀研究[J].巖土工程學(xué)報,2010,32(Supp.1):28-31.
[7] 王艷峰.軟土地區(qū)基坑工程變形性能分析[D].汕頭:汕頭大學(xué),2008.
[8] 李青.軟土深基坑變形性狀的現(xiàn)場試驗研究[D].上海:同濟大學(xué),2008.
[9] 陳忠漢,程麗萍.深基坑工程[M].北京:機械工業(yè)出版社,1999.
[10] 覃志毅.軟土基坑墻后地表沉降估計研究[D].武漢:武漢大學(xué),2004.
[11] 簡艷春.軟土基坑變形估算及其影響因素研究[D].南京:河海大學(xué),2001.
[12] 李廣信.基坑支護結(jié)構(gòu)上水土壓力的分算與合算[J].巖土工程學(xué)報,2000,22(3):348-352.
[13] 郭競宇,趙其華,張建剛.圍護結(jié)構(gòu)上土壓力實例分析[J].巖土工程學(xué)報,2003,25(2):246-248.
[14] 劉國彬,黃院雄,侯學(xué)淵.水及土壓力的實測研究[J].巖石力學(xué)與工程學(xué)報,2000,19(2):205-210.
[15] 劉暢.考慮土體不同強度與變形參數(shù)及基坑支護空間影響的基坑支護變形與內(nèi)力研究[D].天津:天津大學(xué),2008.
Applicability of PLAXIS2D Used for Numerical Simulation in Foundation Pit Excavations
LIU Xiao-Li1,MA Yue2,GUO Guan-Qun1,TAO Tao2,ZHOU He1
(Ocean University of China 1.The Key Lab of Marine Environment and Ecology,Ministry of Education;2.College of Engineering,Qingdao 266100,China)
The geotechnical engineering finite element analysis software,PLAXIS2D,has been widely used for numerical simulation in foundation pit excavations.However,few literatures had made a detail analysis on applicability of the software PLAXIS2D used for foundation pit excavations at present.In view of this,three actual foundation pit excavation engineerings were simulated by PLAXIS2D.By comparison of computing results and measured ones,applicability of PLAXIS2D has been investigated.Further,the application conditions of PLAXIS2D used for foundation pit excavation engineerings have been discussed.The conclusion can provide some references for the proper use of PLAXIS2D in foundation pit excavations.
PLAXIS2D;finite element;foundation pit excavation;numerical simulation;applicability
P642
A
1672-5174(2012)04-019-07
山東省自然科學(xué)基金項目(ZR2009FQ028);中國海洋大學(xué)SRDP項目(1012011004)資助
2011-05-17;
2011-08-25
劉小麗(1974-),女,副教授。E-mail:lxl4791@ouc.edu.cn
責任編輯 龐 旻