呂松濤，陳杰東，張　暉

(長沙理工大學(xué)　公路養(yǎng)護(hù)技術(shù)國家工程實(shí)驗(yàn)室,湖南　長沙　410114)

?

水泥穩(wěn)定碎石拉壓彎靜態(tài)模量與動(dòng)態(tài)模量比較分析

呂松濤，陳杰東，張暉

(長沙理工大學(xué)公路養(yǎng)護(hù)技術(shù)國家工程實(shí)驗(yàn)室,湖南長沙410114)

為了揭示水泥穩(wěn)定碎石半剛性基層材料的動(dòng)、靜態(tài)抗拉模量、抗壓模量和彎拉模量之間的相關(guān)關(guān)系，開發(fā)了同時(shí)測試3種模量的新方法，推導(dǎo)出了在四點(diǎn)彎曲受力狀態(tài)下水泥穩(wěn)定碎石中的抗拉模量、抗壓模量和考慮剪切效應(yīng)影響的彎拉模量計(jì)算公式，開展了不同加載頻率下的動(dòng)態(tài)模量和加載速率為1 mm/min條件下靜態(tài)模量的測試試驗(yàn)；揭示了水泥穩(wěn)定碎石動(dòng)態(tài)模量隨加載頻率的變化規(guī)律，建立了動(dòng)、靜態(tài)模量的轉(zhuǎn)化關(guān)系，揭示了拉、壓、彎3種動(dòng)、靜態(tài)模量的差異性。研究結(jié)果表明：水泥穩(wěn)定碎石材料不是各向同性材料，呈現(xiàn)出明顯的拉、壓模量不同的特性；3種動(dòng)態(tài)模量均隨加載頻率的增長呈現(xiàn)冪函數(shù)增長規(guī)律，動(dòng)態(tài)模量試驗(yàn)條件更加符合路面實(shí)際受力狀態(tài)；相同荷載級(jí)位下，3種動(dòng)態(tài)模量都大于對應(yīng)的靜態(tài)模量，比值區(qū)間約為：1.11～1.91。

道路工程；水泥穩(wěn)定碎石；室內(nèi)試驗(yàn)；動(dòng)態(tài)模量；靜態(tài)模量；加載頻率

0　引言

我國半剛性基層材料回彈模量的測定多以靜態(tài)試驗(yàn)為主，靜態(tài)試驗(yàn)方法中試件受力狀態(tài)與路面結(jié)構(gòu)真實(shí)應(yīng)力狀態(tài)差異較大，不能真實(shí)反映重復(fù)荷載作用下基層材料的動(dòng)態(tài)響應(yīng)性質(zhì)。采用動(dòng)態(tài)試驗(yàn)方法進(jìn)行材料模量測試勢必將成為我國路面結(jié)構(gòu)設(shè)計(jì)的發(fā)展趨勢[1-3]。另一方面，經(jīng)典彈性理論描述的是拉伸、壓縮彈性模量相同的材料，而在工程實(shí)際中，許多材料如混凝土等表現(xiàn)出拉、壓不同的材料性質(zhì)已成為大家的共識(shí)[4-6]。同樣，半剛性基層是拉、壓模量差異較大的結(jié)構(gòu)層，我國在進(jìn)行瀝青路面結(jié)構(gòu)設(shè)計(jì)時(shí)簡單采用無側(cè)限抗壓回彈模量來進(jìn)行結(jié)構(gòu)的荷載響應(yīng)分析[7-8]，造成計(jì)算結(jié)果與實(shí)際情況存在較大的偏差。在路面結(jié)構(gòu)設(shè)計(jì)時(shí)如不考慮以上兩方面因素，勢必導(dǎo)致設(shè)計(jì)精度的降低，亦無法準(zhǔn)確把握交通荷載作用下路用性能的衰變規(guī)律，導(dǎo)致養(yǎng)護(hù)決策的盲目性。

本文利用材料試驗(yàn)系統(tǒng)MTS(Material Test System)，提出了一種拉、壓、彎拉模量同步測試的新方法，開展了不同加載頻率下水泥穩(wěn)定碎石動(dòng)態(tài)彎拉模量、抗壓模量、抗拉模量和1 mm/min加載速率下3種靜態(tài)模量的試驗(yàn)，揭示了3種動(dòng)態(tài)模量隨加載頻率的變化規(guī)律，對比分析了3種動(dòng)、靜態(tài)模量之間的差異性，為優(yōu)化我國瀝青路面結(jié)構(gòu)剛度設(shè)計(jì)參數(shù)提供了試驗(yàn)依據(jù)和理論支撐。

1　配合比設(shè)計(jì)與試件制作養(yǎng)生

1.1原材料

水泥：采用P·C32.5級(jí)水泥，對其進(jìn)行細(xì)度檢測、凝結(jié)時(shí)間測試、安定性檢測以及抗折強(qiáng)度試驗(yàn)，得到其各項(xiàng)技術(shù)性質(zhì)如表1所示。

表1　水泥技術(shù)指標(biāo)Tab.1　Cement technical indexes

集料：采用某碎石廠生產(chǎn)的石灰?guī)r集料，公稱最大粒徑為19 mm，原材料各項(xiàng)技術(shù)指標(biāo)均滿足規(guī)范規(guī)程要求。

表2　集料性質(zhì)試驗(yàn)結(jié)果Tab.2　Test result of aggregate property

1.2配合比試驗(yàn)

采用的礦料級(jí)配如表3所示。

表3　水泥穩(wěn)定碎石礦料合成級(jí)配Tab.3　Synthetic aggregate gradation of cement stabilized macadam

根據(jù)表3確定的礦料合成級(jí)配，進(jìn)行水泥穩(wěn)定碎石的重型擊實(shí)試驗(yàn)，首先確定不同水泥劑量下的最佳含水量與最大干密度，而后通過各水泥劑量下的無側(cè)限抗壓強(qiáng)度試驗(yàn)確定滿足要求的水泥劑量，最終試驗(yàn)結(jié)果為：水泥穩(wěn)定碎石材料的外摻水泥劑量為4.5%，最大干密度為2.35 g/cm3，最佳含水量為4.5%。

1.3試件制備養(yǎng)生

據(jù)此，按照試驗(yàn)規(guī)程的要求[9]，制作彎拉強(qiáng)度及彎拉模量梁式試件，其尺寸為100 mm×100 mm×400 mm，試件成型均采用靜壓成型的試驗(yàn)方法，成型后，試件在溫度為(20±2)℃、濕度≥95%的標(biāo)準(zhǔn)養(yǎng)生室內(nèi)進(jìn)行90 d標(biāo)準(zhǔn)齡期的養(yǎng)生。

2　模量試驗(yàn)方法與計(jì)算公式

2.1動(dòng)靜模量試驗(yàn)方法

路面材料靜態(tài)模量是在逐級(jí)加、卸載條件下，根據(jù)材料受力和相應(yīng)的回彈變形計(jì)算得出的模量?，F(xiàn)行試驗(yàn)規(guī)程對于水泥穩(wěn)定碎石材料彎拉模量的測試采用的是三分點(diǎn)加載的方式。將最大施加荷載(本論文采用極限破壞荷載的30%)均分為5～6級(jí)，持續(xù)穩(wěn)定地施加第1級(jí)荷載的同時(shí)，利用LVDT位移計(jì)持續(xù)采集試件中部撓度變化，荷載施加時(shí)間為1 min；然后全部卸去荷載并保持0.5 min，讓試件充分回彈，LVDT全程記錄試件跨中撓度變化；之后逐級(jí)加載、卸載，得到多級(jí)荷載作用下的荷載-回彈變形曲線，加載、卸載速率均為1 mm/min。

在進(jìn)行彎拉模量測試的基礎(chǔ)上，在試驗(yàn)梁跨中上下表面分別對稱粘貼4片應(yīng)變片，共8片，并分別連接動(dòng)態(tài)應(yīng)變采集儀，以記錄上表面的壓應(yīng)變和下表面的拉應(yīng)變，用于抗壓模量和抗拉模量的計(jì)算。由于水穩(wěn)碎石試件表面一般比較粗糙，試驗(yàn)前需在粘貼應(yīng)變片處用水泥漿填平表面縫隙，待水泥漿凝固后用砂紙打磨處理，避免水泥漿在試件表面成層。測試原理圖和應(yīng)變片粘貼方式分別如圖1、圖2所示。

1—LVDT位移計(jì)；2—應(yīng)變片；3—圓柱形支座；4—加載頭圖1　動(dòng)靜態(tài)拉壓彎模量測試原理圖Fig.1　Testing principle of tensile, compression, and flexural static and dynamic moduli

圖2　應(yīng)變片粘貼示意圖Fig.2　Schematic diagram of strain gauges pasting

動(dòng)態(tài)模量為材料在特定加載頻率下受到荷載循環(huán)往復(fù)作用，根據(jù)對應(yīng)的應(yīng)力-應(yīng)變關(guān)系計(jì)算得出的模量。其與靜態(tài)模量測試方法的區(qū)別在于加載方式的不同，試驗(yàn)采用半正矢波加載，加載頻率分為1，4，7，10 Hz 4個(gè)等級(jí)，荷載等級(jí)采用30%試件破壞荷載。荷載循環(huán)作用200次，取最后10次加載循環(huán)的平均荷載和變形振幅用于計(jì)算3種動(dòng)態(tài)模量。

2.2動(dòng)靜模量計(jì)算方法

當(dāng)梁的跨高比大于5時(shí)，純彎曲的正應(yīng)力公式可以推廣應(yīng)用于橫力彎曲的正應(yīng)力計(jì)算。然而對于跨高比小于5的試驗(yàn)梁，用材料力學(xué)細(xì)長梁的純彎曲理論及假設(shè)計(jì)算時(shí)，誤差會(huì)隨高跨比減小而迅速增大[10]?，F(xiàn)行無機(jī)結(jié)合料試驗(yàn)規(guī)程針對彎拉模量采用的是跨高比l/h為3的比例，并未考慮剪切作用對梁撓度的影響，故與實(shí)際情況不符[9]。本文在規(guī)程原公式的基礎(chǔ)上，推導(dǎo)出了考慮剪切作用影響的彎拉模量計(jì)算修正公式(后文提及的彎拉模量均為考慮剪切影響后的彎拉模量)。

根據(jù)文獻(xiàn)[9]，不考慮剪切影響的彎拉模量與跨中撓度存在以下關(guān)系：

(1)

式中，ω為跨中撓度；Ef為不考慮剪切影響的彎拉模量；P為施加荷載；L為試件跨徑；b為跨中截面寬度；h為跨中截面高度。

當(dāng)考慮剪切作用影響時(shí)，剪切作用引起的附加撓度在圖1所示的四點(diǎn)彎曲試驗(yàn)條件下有

(2)

(3)

則彎、剪共同作用引起的撓度ω′可表示為：ω′=ω+Δω。將以上各式代入可得：

(4)

(5)

對于試件拉、壓模量，文獻(xiàn)[10]由平面假設(shè)和平衡條件，利用頂部壓應(yīng)變和底部拉應(yīng)變得到了四點(diǎn)彎曲作用下，梁式試件的抗壓模量、抗拉模量計(jì)算公式：

(6)

(7)

3　試驗(yàn)結(jié)果與分析

3.1動(dòng)靜態(tài)模量試驗(yàn)結(jié)果

針對3種靜態(tài)模量進(jìn)行5組試驗(yàn)分析，分析時(shí)泊松比μ取值0.25。計(jì)算分析結(jié)果如表4所示。

表4　靜態(tài)拉壓彎模量試驗(yàn)結(jié)果Tab.4　Test result of tensile, compression, and flexural static moduli

從表4可以看出：(1)拉模量均值為7 964 MPa，壓模量均值為14 620 MPa，彎拉模量均值為777 MPa，壓模量遠(yuǎn)遠(yuǎn)大于彎拉模量和拉模量，彎拉模量遠(yuǎn)遠(yuǎn)小于拉、壓模量。(2)壓模量與彎拉模量比值約為18.81，壓模量與拉模量比值約為1.84，拉模量與彎拉模量比值約為10.25。這表明水泥穩(wěn)定碎石材料是一種拉、壓模量不同的各向異性材料。

實(shí)驗(yàn)室條件下的加載頻率取決于外部荷載對路面結(jié)構(gòu)的振動(dòng)頻率?，F(xiàn)行半剛性材料動(dòng)態(tài)模量測試借鑒國外經(jīng)驗(yàn)，采用10 Hz的標(biāo)準(zhǔn)來表征車輛60～80 km/h速度下材料的受荷頻率，然而對于其他車速對應(yīng)的加載頻率并未加以考慮，因此有必要開展不同加載頻率下的動(dòng)態(tài)模量測試，并揭示動(dòng)靜模量兩者之間的對應(yīng)關(guān)系。

不同加載頻率下3種動(dòng)態(tài)模量的測試結(jié)果如表5所示。

表5　不同加載頻率下的動(dòng)態(tài)拉壓彎模量試驗(yàn)結(jié)果Tab.5　Test result of tensile, compression, and flexural moduli under different loading frequencies

由表5可看出，不同加載頻率下壓模量與彎拉模量比值分別為：13.40，12.17，11.63，11.13，比值均值為12.08；壓模量與拉模量比值分別為：2.01，2.00，2.02，2.02，比值均值為2.01;拉模量與彎拉模量比值分別為：6.66，6.07，5.77，5.52，比值均值6.01。在動(dòng)態(tài)荷載下，壓模量亦遠(yuǎn)遠(yuǎn)大于拉模量和彎拉模量，水泥穩(wěn)定碎石仍然表現(xiàn)出明顯的拉、壓模量不同的力學(xué)性質(zhì)。

3種動(dòng)態(tài)模量隨加載頻率的變化規(guī)律可用冪函數(shù)表征。同時(shí)，因?yàn)殪o態(tài)模量測試時(shí)，對應(yīng)的加載頻率應(yīng)接近于0。因此，將靜態(tài)模量視為加載頻率為0時(shí)的取值考慮其中，得到了更為普適的模量-加載頻率經(jīng)驗(yàn)關(guān)系：

(8)

式中，E為動(dòng)態(tài)拉、壓、彎模量；E0為靜態(tài)拉、壓、彎模量；f為加載頻率；a，b為擬合參數(shù)。

擬合過程如圖3所示。

擬合結(jié)果匯總見表6。

研究表明造成動(dòng)態(tài)模量增長的原因是多樣的：有材料自身因素的影響，如試件本身壓密成型，存在進(jìn)一步的壓緊密實(shí)的過程，導(dǎo)致在循環(huán)荷載作用下回彈變形的減小[12]；也存在試驗(yàn)條件的影響，在相同荷載條件下，荷載作用頻率越高，對材料造成的激振頻率越快，試件來不及產(chǎn)生相應(yīng)的變形，導(dǎo)致試件的實(shí)際變形與響應(yīng)回彈變形存在的滯后現(xiàn)象[13]；另外試件與設(shè)備之間的摩擦也會(huì)導(dǎo)致試件的結(jié)構(gòu)阻尼增大，即材料的黏性變形越來越大[14]，以上因素綜合導(dǎo)致了動(dòng)態(tài)模量的逐級(jí)增加。從以上擬合結(jié)果可見，隨著加載頻率的增長，水泥穩(wěn)定碎石動(dòng)態(tài)模量呈現(xiàn)出較好的冪函數(shù)增長規(guī)律。

表6　不同加載頻率下的3種動(dòng)態(tài)模量擬合結(jié)果Tab.6　Fitting result of 3 dynamic moduli under different loading frequencies

3.2動(dòng)靜態(tài)模量比較分析

(1)靜態(tài)荷載作用下的3種模量均小于加載頻率為1 Hz下的動(dòng)態(tài)模量，考慮到車速60～80 km/h的汽車荷載對應(yīng)于10 Hz下的動(dòng)態(tài)模量值，可見靜態(tài)模量值不適合表征高速行駛車輛荷載作用下的路面真實(shí)受力狀態(tài)，動(dòng)態(tài)模量更能反映路面的實(shí)際情況。

(2)不同加載頻率下動(dòng)、靜態(tài)壓模量比值分別為：1.19，1.21，1.24，1.24，比值均值為1.22；拉模量比值分別為：1.08，1.11，1.13，1.14，比值均值為1.11；彎拉模量比值分別為：1.67，1.87，2.01，2.11，比值均值為1.91。這表明不同加載頻率下，3種動(dòng)態(tài)模量都大于靜態(tài)模量，用靜態(tài)模量分析路面結(jié)構(gòu)力學(xué)響應(yīng)是不合理的。

(3)從大量試驗(yàn)中可知路面常用材料在受力狀態(tài)不同時(shí)，會(huì)表現(xiàn)出不同的應(yīng)力-應(yīng)變關(guān)系，即表現(xiàn)出不同的拉、壓、彎模量特征。對于路面材料而言，早有學(xué)者提出路面設(shè)計(jì)中簡單采用無側(cè)限抗壓回彈模量作為結(jié)構(gòu)層模量的弊端，推導(dǎo)出了雙模量理論的基本公式及試驗(yàn)方法[11]。從表4、表5的結(jié)果中可以看出，無論動(dòng)態(tài)荷載還是靜態(tài)荷載作用下，水泥穩(wěn)定碎石材料的壓模量都遠(yuǎn)遠(yuǎn)大于拉模量和彎拉模量，進(jìn)一步佐證了水泥穩(wěn)定碎石材料是一種拉、壓模量不同的各向異性材料。這表明現(xiàn)規(guī)范用單一的無側(cè)限抗壓回彈模量來進(jìn)行路面結(jié)構(gòu)設(shè)計(jì)計(jì)算是不科學(xué)的，而應(yīng)根據(jù)結(jié)構(gòu)層不同位置所處的實(shí)際應(yīng)力狀態(tài)選用不同的模量類型。

4　結(jié)論

(1)水泥穩(wěn)定碎石材料在動(dòng)、靜態(tài)荷載作用下，都表現(xiàn)出拉壓模量不同的各向異性特性。用單一的靜態(tài)模量來表征路面受力狀態(tài)是不合適的，而應(yīng)根據(jù)路面材料的實(shí)際受力狀態(tài)選取相應(yīng)的模量參數(shù)來進(jìn)行結(jié)構(gòu)應(yīng)力分析。

(2)3種動(dòng)態(tài)模量隨加載頻率的增長呈現(xiàn)冪函數(shù)增長關(guān)系。

(3)3種動(dòng)態(tài)模量都大于靜態(tài)模量，動(dòng)、靜態(tài)壓模量比值約為1.22，拉模量比值約為1.11，彎拉模量比值約為1.91。動(dòng)態(tài)模量更適合于分析路面結(jié)構(gòu)在車輛荷載作用下的力學(xué)響應(yīng)。

[1]沙愛民,賈侃,陸劍卿，等.半剛性基層材料動(dòng)態(tài)模量的衰變規(guī)律[J].中國公路學(xué)報(bào),2009,22(3):1-6.

SHA Ai-min,JIA Kan,LU Jian-qing，et al. Deterioration Laws of Dynamic Modulus of Semi-rigid Base Course Materials[J]. China Journal of Highway and Transport,2009,22(3):1-6.

[2]劉偉,嚴(yán)金海,李鋒，等.乳化瀝青冷再生混合料動(dòng)靜模量相關(guān)性[J].公路交通科技,2015,32(5):1-6.

LIU Wei,YAN Jin-hai,LI Feng, et al. Correlation between Dynamic Modulus and Static Modulus of Cold Recycled Mixture with Emulsified Asphalt [J]. Journal of Highway and Transportation Research and Development, 2015,32(5):1-6.

[3]孫建.瀝青混合料動(dòng)態(tài)模量研究[D].西安:長安大學(xué)，2007.

SUN Jian. Research on Dynamic Modulus of Asphalt Mixture[D]. Xi’an: Chang’an University,2007.

[4]張洪武,張亮，高強(qiáng).拉壓不同模量材料的參變量變分原理和有限元方法[J]. 工程力學(xué),2012,29(8):22-17.

ZHANG Hong-wu, ZHANG Liang, GAO Qiang. The Parametric Variational Principle and Finite Element Method for Material with Different Modulus in Tension and Compression[J]. Engineering Mechanics, 2012,29(8):22-17.

[5]孔娟,袁聚云,潘曉明，等.不同拉壓模量厚壁球殼彈塑性分析[J].力學(xué)與實(shí)踐,2010,32(1):41-45.KONG Juan,YUAN Ju-yun,PAN Xiao-ming,et al. Elasto-plastic Analysis of Thick Spheric Shell with Different Elastic Moduli for Tensile and Compressive Deformations [J].Mechanics in Engineering, 2010,32(1):41-45.

[6]PATEL B P, KHAN K, NATH Y. A New Constitutive Model for Bimodular Laminated Structures: Application to Free Vibrations of Conical/Cylindrical Panels[J]. Composite Structures,2014,110(3):183-191.

[7]JTG D50—2006,公路瀝青路面設(shè)計(jì)規(guī)范[S].

JTG D50—2006, Specifications for Design of Highway Asphalt Pavement[S].

[8]姚祖康.瀝青路面結(jié)構(gòu)設(shè)計(jì)[M].北京：人民交通出版社，2011.

YAO Zu-kang. Design of Asphalt Pavement Structure [M].Beijing: China Communications Press,2011.

[9]JTG E51—2009,公路工程無機(jī)結(jié)合料穩(wěn)定材料試驗(yàn)規(guī)程[S].

JTG E51—2009, Test Methods of Materials Stabilized with Inorganic Binders for Highway Engineering [S].

[10]王中正,朱軍祚,諶磊，等.集中力作用下深梁彎剪耦合變形應(yīng)力計(jì)算方法[J].工程力學(xué),2008,25(4):115-119.

WANG Zhong-zheng, ZHU Jun-zuo, CHEN Lei,et al. The Stress Caculation Method for Deep Beams with Shear-bending Coupling Distortion under Concentrated Load [J]. Engineering Mechanics, 2008,25(4):115-119.

[11]劉金蓮,應(yīng)榮華.雙模量理論在柔性路面設(shè)計(jì)中的應(yīng)用研究[J].湖南交通科技,2001,27(1):16-23.

LIU Jin-lian,YING Rong-hua. Application Research of Double Modulus Theory in Flexible Pavement Design[J].Hunan Communications Science and Technology,2001,27(1):16-23.

[12]程箭,許志鴻,李淑明，等. 水泥穩(wěn)定碎石靜態(tài)模量與動(dòng)態(tài)模量比較[J].建筑材料學(xué)報(bào),2009,12(1):63-75.

CHEN Jian,XU Zhi-hong,LI Shu-ming,et al. Comparison of Compressive Resilient Modulus and Dynamic Modulus of Cement-Stabilized Macadam[J].Journal of Building Materials, 2009,12(1):63-75.

[13]郭大進(jìn),黑鵬,韓立志.半剛性基層材料動(dòng)態(tài)特性研究[J].公路交通科技,2005,22(12):42-61.

GUO Da-jin, HEI Peng, HAN Li-zhi. Study on Dynamic Characteristics of Semi-rigid Base Layer [J].Journal of Highway and Transportation Research and Development,2005,22(12):42-61.

[14]ORUC S, CELIK F, AKPINAR M V. Effect of Cement on Emulsified Asphalt Mixtures[J]. Journal of Materials Engineering & Performance, 2007, 16(5):578-583.

Comparative Analysis of Tensile, Compression, Flexural Static Modulus and Dynamic Modulus of Cement-stabilized Macadam

Lü Song-tao，CHEN Jie-dong，ZHANG Hui

(State Engineering Laboratory of Highway Maintenance Technology, Changsha University of Science & Technology,Changsha Hunan 410114, China)

In order to reveal the relations of tensile, compression and flexural resilient moduli with dynamic moduli of cement-stabilized macadam semi-rigid base, a new test method for measuring 3 moduli simultaneously is developed. The formulas of tensile, compression, and flexural modulus which considering the influence of shearing effect for cement-stabilized macadam under 4-point bending are derived. The measurment experiments of dynamic modulus under different loading frequencies and static modulus at 1 mm/min loading speed are conducted. The change rules of tensile, compression and flexural moduli with loading frequency are revealed, and the transformation relationships of static and dynamic moduli are established. The differences of tensile, compression and flexural static modulus and dynamic modulus are revealed. The result indicates that (1) cement-stabilized macadam is not isotropic material, the tensile and compression moduli show a significant difference obviously; (2) the variations of 3 dynamic moduli show a power function growth relationship with the increase of loading frequency, the testing condition of dynamic modulus is closer to the actual stress condition of the pavement; (3) three dynamic moduli are larger than the corresponding static modulus under the same load level, the ratio range is from 1.11 to 1.91 approximately.

road engineering; cement stabilized macadam; laboratory test; dynamic modulus; static modulus; loading frequency

2015-11-12

國家自然科學(xué)基金項(xiàng)目(51208066，51578081)

呂松濤(1979-)，男，河南鄭州人，博士,副教授.(58003611@qq.com)

10.3969/j.issn.1002-0268.2016.10.007

U416.217

A

1002-0268(2016)10-0039-05