青藏高速公路路基地震動力響應分析

趙濤，吳志堅，梁慶國，趙文琛，陳拓，徐世民

(1. 蘭州交通大學 甘肅省道路橋梁與地下工程重點實驗室，甘肅 蘭州 730070；2. 中國地震局蘭州地震研究所 黃土地震工程重點實驗室，甘肅 蘭州 730000；3. 中國市政工程西北設(shè)計研究院有限公司，甘肅 蘭州 730000；4. 中國礦業(yè)大學 深部巖土力學與地下工程國家重點實驗室，江蘇 徐州 221116)

青藏高速公路路基地震動力響應分析

趙濤1,2，吳志堅2，梁慶國1，趙文琛3，陳拓4，徐世民2

(1. 蘭州交通大學 甘肅省道路橋梁與地下工程重點實驗室，甘肅 蘭州 730070；2. 中國地震局蘭州地震研究所 黃土地震工程重點實驗室，甘肅 蘭州 730000；3. 中國市政工程西北設(shè)計研究院有限公司，甘肅 蘭州 730000；4. 中國礦業(yè)大學 深部巖土力學與地下工程國家重點實驗室，江蘇 徐州 221116)

為了研究即將修筑的青藏高速公路在地震荷載作用下的動力響應特征，擬合出該區(qū)50 a超越概率為63.5%，10.0%和2.0%的人工地震波。在數(shù)值模擬的基礎(chǔ)上，運用二維非線性動力有限元分析方法，建立整體式和分離式2種典型的高速公路路基橫斷面結(jié)構(gòu)，對比分析不同超越概率地震荷載作用下2種形式路基的動力響應特征。研究結(jié)果表明：分離式路基在路基表面中心和天然地表中心2處的加速度放大系數(shù)均小于整體式路基，而其路基表面的水平、垂直峰值位移均大于整體式路基；2種路基左右兩側(cè)坡腳下方土體內(nèi)剪應力集中分布；隨著輸入地震波峰值加速度的增大，2種路基表面的水平峰值位移均呈線性增長趨勢，垂直峰值位移呈指數(shù)函數(shù)增長趨勢。分析結(jié)果可為多年凍土區(qū)高等級公路抗震設(shè)計提供相應參考。

青藏高速公路；地震波；典型路基；動力響應

青藏高原地區(qū)是世界性高海拔多年凍土區(qū)的典型代表[1]。在寬度不足10 km范圍內(nèi)的青藏凍土工程走廊內(nèi)，擬建與已建的國家重大基礎(chǔ)工程密集分布，高原多年凍土的脆弱性決定了這些重大工程的穩(wěn)定性。加之青藏高原地區(qū)又是地震多發(fā)區(qū)，其構(gòu)造活動強烈，地震活動頻繁，強震較多[2-3]。而該區(qū)域內(nèi)的場地和建(構(gòu))筑物抗震設(shè)防的基礎(chǔ)研究基本處于空白，青藏高速公路也已開展了可行性研究，立項在即。規(guī)劃的青藏高速公路從西寧市至拉薩市全長約1 900 km，屬于京藏高速公路(國家高速公路網(wǎng)編號G6)的一部分。目前，因地質(zhì)條件、施工設(shè)備，尤其核心技術(shù)瓶頸是面臨跨越500多km的多年凍土區(qū)等因素制約項目建設(shè)，格爾木-拉薩長約1 100 km路段尚待修建。因為青藏高原地理位置的特殊性、氣候環(huán)境的復雜性以及多年凍土的溫度敏感性，使青藏高速公路的鋪筑一貫被認為是難以攻克的問題。目前，國內(nèi)外鮮有多年凍土地區(qū)修建高等級公路的經(jīng)驗可供借鑒，同時，相關(guān)學術(shù)文獻也少有提及。相關(guān)學者們對于機車振動荷載[4-8]和地震動荷載[9-10]作用下青藏鐵路路基結(jié)構(gòu)的動力響應特征已經(jīng)做了初步研究，其成果對于青藏高原多年凍土地域修筑窄幅路基等構(gòu)筑物提供重要的參考依據(jù)。然而，路基寬度會直接影響凍土路基的傳熱過程、強度以及穩(wěn)定性。青藏鐵路、青藏公路的路基寬度較窄，而高速公路路基相對較寬，路基吸熱強度將成倍提高，導致凍土退化速率和范圍的增加，對天然凍土的擾動加大[11-12]。共和-玉樹高速公路穿越多年凍土區(qū)，現(xiàn)僅有集中于溫度效應方面[13-15]的研究分析。另外，符進等[16]遴選出分離式和整體式路基適當路段進行了相關(guān)試驗示范研究。但是對于多年凍土區(qū)修筑的高等級公路受地震動力荷載作用的影響還未引起關(guān)注。因此，對于青藏高原多年凍土區(qū)即將修筑的高速公路，在地震荷載作用下的變形特征及動力響應特性，已成為必須考慮和亟待解決的問題。鑒于此，筆者根據(jù)青藏高原內(nèi)的地震地質(zhì)條件、地球物理特征以及地震活動性條件，擬合出北麓河多年凍土區(qū)50年超越概率為63.5%，10.0%和2.0%的人工地震波，在數(shù)值模擬的基礎(chǔ)上，運用二維非線性動力有限元分析方法，建立整體式和分離式2種典型的高速公路路基橫斷面結(jié)構(gòu)，對比分析不同超越概率地震荷載作用下多年凍土區(qū)高速公路路基的動力響應特征。研究結(jié)果對于即將修建的青藏高速公路格拉段典型路基結(jié)構(gòu)形式在地震荷載作用下的動力破壞失穩(wěn)機理具有迫切的現(xiàn)實意義和參考作用。

1 北麓河試驗區(qū)概況

青藏高原北麓河多年凍土地區(qū)海拔高程4 620～4 640 m，年平均氣溫-5.2 ℃，該試驗段地下冰分布廣泛，全區(qū)域內(nèi)有含土冰層分布，其體積含冰量超過50%，多為懸浮狀結(jié)構(gòu)，主要分布在凍土上限以下1～3 m范圍內(nèi)，厚度一般為1～3 m。凍土上限一般為2～3 m，最深可達3.4 m，最淺為1.7 m。除少數(shù)區(qū)域外，該試驗段多年凍土區(qū)的年均稱地溫為-1.8～-0.4 ℃。[9]

2 合成模擬地震波

盡管青藏高原地區(qū)地震頻發(fā)，但由于前些年高原區(qū)域人口分布稀疏，重要建筑較少，同時監(jiān)測技術(shù)相對滯后等原因，青藏地震波實測記錄實屬罕見。因此，根據(jù)青藏高原區(qū)域內(nèi)的地震地質(zhì)條件、地球物理特征以及地震活動性條件，加速度時程曲線按照人工合成地震動的方法，擬合出50 a超越概率分別為63.5%，10.0%和2.0%的人工地震動加速度時程，如圖1所示。人工擬合地震波加速度曲線中，需引入擬合地震動的強度包絡(luò)函數(shù)，用以控制時程的隨機性和收斂約束性，采用逐步逼近目標譜的方法，使合成的加速度時程精度滿足目標峰值加速度，并近似滿足目標加速度反應譜，時程間隔為0.02 s，擬合60個目標譜控制點，分布于0.04～6 s之間，擬合地震波加速度曲線的反應譜與目標譜之間的相對誤差限值為5%。3種不同超越概率的加速度時程特征參數(shù)見表1。

(a)50 a超越概率63.5%;(b)50 a超越概率10%;(c)50 a超越概率2%圖1 不同超越概率的人工地震動加速度時程Fig.1 Artificial acceleration time histories of ground motion for different exceeding probability

表1 加速度時程特征參數(shù)Table 1 Characteristic parameters of acceleration time histories

3 路基模型的建立

3.1 2種典型高速公路路基橫斷面設(shè)計

依據(jù)已開通的北京至格爾木段高速公路路基斷面形式，青藏高速公路格拉段的設(shè)計行車速度為100 km/h，采取雙向4車道。據(jù)文獻[17]，設(shè)計出2種標準路基斷面形式，分別為寬度24 m的整體式和寬度12 m的分幅式路基。其中，行車道3.75 m，中央分隔帶寬2 m(僅整體式)，路緣帶為0.75 m，行車道兩側(cè)的硬路肩寬2.5 m，硬路肩外側(cè)為0.75 m的土路肩，分幅路基兩內(nèi)側(cè)坡腳間距10 m。2種路基形式頂面均設(shè)置2%的橫坡坡度，邊坡坡率均為1∶1.5，路基高度4.8 m。根據(jù)現(xiàn)場波速鉆孔資料，試驗區(qū)土層有成層分布特征，分別建立2種形式路基的有限元分析模型如圖2所示。

單位：m(a)整體式路基橫斷面圖;(b)分離式路基橫斷面圖圖2 路基橫斷面Fig.2 Cross-section of embankment

3.2 土體參數(shù)及模型分析方法

各土層和路基填料的厚度和物理力學參數(shù)依據(jù)試驗區(qū)現(xiàn)場鉆孔資料[5]及文獻[11,17-18]確定，見表2。模型土體采用平面應變單元，同時為了減小邊界效應的影響，地基兩側(cè)采用無限元邊界，底部為固定邊界。計算過程中，首先進行結(jié)構(gòu)自重計算，將所得到的應力作為初始應力條件，再建立動力分析步，通過在模型結(jié)構(gòu)底部輸入3種不同超越概率的地震波，對2種不同結(jié)構(gòu)形式路基進行動力響應特征分析。

表2 數(shù)值計算土層物理力學參數(shù)Table 2 Physical and mechanical parameters of soil in numerical calculation

4 路基的動力響應分析

基于ABAQUS有限元軟件，通過二維非線性動力有限元分析方法，在3種不同超越概率地震荷載作用下，對整體式和分離式2種填土路基結(jié)構(gòu)進行動力分析，得到了2種形式路基的加速度，應力及位移動力響應計算結(jié)果。

4.1 加速度響應特征

圖3為整體式和分離式2種結(jié)構(gòu)路基在3種超越概率地震動作用下路基表面中心位置處(圖2中A，A'點)的加速度時程曲線。由圖3可知，隨著地震波峰值加速度的增大，2種形式路基表面加速度均有明顯增加趨向。圖4為不同超越概率地震動作用下，2種結(jié)構(gòu)形式路基分別在路基表面中心和路基下伏的天然地表中心(圖2中B，B'點)2處的峰值加速度相對于基巖輸入峰值加速度的放大系數(shù)。由圖4可以得出，在峰值加速度分別為51，150和252 cm/s2的地震作用下，整體式路基A，B2點的加速度放大系數(shù)分別為5.3，6.7，5.7，3.9，6.5和5.1。分離式路基A'，B' 2點的加速度放大系數(shù)分別為5.3，6.4，4.7，3.7，5.4和4.8。由特征點的加速度放大特征曲線規(guī)律可知，在3種超越概率的地震動作用下，分離式路基在路基表面和天然地表中心2處的加速度放大系數(shù)均小于整體式路基。其中，在10%設(shè)防地震作用下路基加速度放大系數(shù)值大于罕遇地震和多遇地震，可認為在多年凍土區(qū)修筑高等級公路路基，通過參考整體式路基的設(shè)防地震加速度放大系數(shù)值以滿足設(shè)防要求。

為便于觀察在3種地震荷載作用下加速度放大效應沿土層高度方向的變化規(guī)律，分別選取圖2斷面的4個特征點A，B，C和D(其中C點深度2 m，D點深度5.9 m)來對比分析其放大系數(shù)曲線。圖5(a)和5(b)分別為2種路基加速度放大系數(shù)沿高度方向的傳播情況。從圖中可以看出，沿著土層高度方向，2種路基的加速度放大系數(shù)均呈增大趨勢；地震波在天然地表以下的C和D點的加速度放大系數(shù)比較接近，受不同峰值地震荷載的影響變化不大，當?shù)卣鸩▊鞑サ教烊坏乇碇敝恋竭_路基表面的過程中，其放大系數(shù)受荷載峰值變化比較明顯，超越概率10%的地震波變化尤為顯著。

(a)整體式路基表面中心(A點);(b)分離式路基表面中心(A'點)圖3 不同路基表面中心加速度時程曲線Fig.3 Curves of acceleration time histories at the center of surface for different embankment

(a)路基表面中心;(b)天然地表中心圖4 不同路基特征點的加速度放大系數(shù)Fig.4 Amplification coefficient of acceleration of characteristic points for different embankment

(a)整體式路基；(b)分離式路基圖5 不同路基沿高度方向的加速度放大系數(shù)Fig.5 Amplification coefficient of acceleration along height for different embankment

4.2 應力響應特征

圖6為超越概率10%的青藏人工地震波(150 cm·s-2)作用下，2種不同結(jié)構(gòu)形式路基的最大應力響應云圖，包括Mises應力和剪應力的響應云圖。從應力響應云圖可以看出，地震荷載作用下2種形式路基的應力值成層分布，并且有相同數(shù)量級的數(shù)值，其中分離式路基表面(3.62 kPa)的Mises應力值略大于整體式(3.22 kPa)。在兩側(cè)坡腳下方土體內(nèi)均有明顯的剪應力分布。從3種工況下路基的應力響應分布可以得出，隨著地震波峰值加速度的增大，2種結(jié)構(gòu)形式路基的最大應力響應均有增加趨勢。

單位：Pa(a)整體式路基Mises應力云圖；(b)整體式路基剪應力云圖；(c)分離式路基Mises應力云圖；(d)分離式路基Mises應力云圖圖6 不同路基的應力響應云圖Fig.6 Nephogram of stress response for different embankment

4.3 位移響應特征

為了進一步探究3種超越概率地震荷載作用下2種路基表面土體的變形特征，分別提取路基表面特征點(圖2中A，A'點)和模型底部(圖2中E，E'點)的水平和垂直位移時程曲線，并做差值得到路基表面的水平絕對位移和垂直絕對位移的時程曲線。在絕對位移時程曲線上分別提取A，A'點的水平峰值位移U1，垂直峰值位移U2，其與輸入地震波峰值加速度的變化關(guān)系如圖7(a)和7(b)所示。由圖7可知，在相同地震荷載作用下，整體式路基表面水平、垂直峰值位移均小于分離式路基。隨著輸入地震波峰值加速度的增大，2種路基表面的水平峰值位移均呈線性增長趨勢，約相差3mm；垂直峰值位移呈較好的指數(shù)函數(shù)增長趨勢。因而，在類似強度特征的地震荷載作用下，如若地震動持時較長，分離式路基表面位移的變化特征較整體式更為顯著，在宏觀上將體現(xiàn)于路基結(jié)構(gòu)更容易發(fā)生水平開裂。

(a)水平位移；(b)垂直位移圖7 水平、豎直位移與地震波峰值加速度關(guān)系曲線Fig.7 Relationship curves betweenhorizontal displacement, vertical displacement and seismic peak acceleration

5 結(jié)論

1)分離式路基在路基表面中心和天然地表中心2處的加速度放大系數(shù)均小于整體式路基。沿著土層高度方向，2種路基的加速度放大系數(shù)均呈增大趨勢；地震波在天然地表以下加速度放大系數(shù)受不同地震荷載峰值的影響變化不大，當?shù)卣鸩▊鞑サ教烊坏乇碇敝恋竭_路基表面的過程中，其放大系數(shù)受荷載峰值變化比較明顯。

2)整體式和分離式路基的應力均為成層分布，在兩側(cè)坡腳下方土體內(nèi)剪應力集中分布，隨著地震波峰值加速度的增大，2種結(jié)構(gòu)形式路基的最大應力均有增加趨勢。

3)整體式路基表面水平、垂直峰值位移均小于分離式路基。隨著輸入地震波峰值加速度的增大，2種路基表面的水平峰值位移均呈線性增長趨勢；垂直峰值位移呈較好的指數(shù)函數(shù)增長趨勢。

4)本文僅考慮地震荷載作用下整體式和分離式2種典型結(jié)構(gòu)形式路基的動力響應特征。由于青藏高原工程走廊內(nèi)受多因素耦合疊加影響，下一步應該結(jié)合溫度場、滲流場等多場耦合情況進行分析。此外，由于青藏工程走廊帶地形、地質(zhì)條件的復雜性，對于地震荷載作用下，有必要進一步將多年凍土上限變化，路基設(shè)計高度，邊坡坡率和分幅間距等影響因素聯(lián)合進行探討分析。

[1] 馬巍，牛富俊，穆彥虎．青藏高原重大凍土工程的基礎(chǔ)研究[J]．地球科學進展，2012，27(11)：1185-1191． MA Wei, NIU Fujun, MU Yanhu. Basic research on the major permafrost projects in the Qinghai-Tibet plateau[J]. Advances in Earth Science, 2012, 27(11): 1185-1191.

[2] 吳志堅，孫軍杰，王蘭民，等.青藏鐵路沿線凍土場地地震動特征研究[J].巖石力學與工程學報，2007,26(12):2466-2472. WU Zhijian, SUN Junjie, WANG Lanmin, et al. Study on characteristics of ground motion at permafrost sites along Qinghai-Tibet railway[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(12): 2466-2472.

[3] 吳志堅，王蘭民，孫軍杰，等.青藏高原典型多年凍土場地地面脈動特征[J].巖石力學與工程學報，2008,27(11)：2316-2323. WU Zhijian, WANG Lanmin, SUN Junjie, et al. Characteristics of microseism at typical permafrost sites in Qinghai-Tibet plateau[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(11): 2316-2323.

[4] ZHU Zhanyuan, LING Xianzhang, CHEN Shijun, et al. Analysis of dynamic compressive stress induced by passing trains in permafrost subgrade along Qinghai-Tibet Railway [J]. Cold Regions Science and Technology, 2011, 65: 465-473.

[5] 吳志堅，陳拓，馬巍.機車荷載作用下青藏鐵路多年凍土區(qū)普通路基的蠕變分析[J].巖土力學，2011，32(增2)：83-87. WU Zhijian, CHEN Tuo, MA Wei. Creep analysis of plain fill embankment at permafrost regions along Qinghai-Tibet Railway under train dynamic load[J]. Rock and Soil Mechanics, 2011, 32(Supp.2): 83-87.

[6] CHEN Tuo, MA Wei, WU Zhijian, et al. Characteristics of dynamic response of the active layer beneath embankment in permafrost regions along the Qinghai-Tibet railroad[J]. Cold Regions Science and Technology, 2014, 98: 1-7.

[7] 田立慧，凌賢長，王立娜，等.青藏鐵路高溫多年凍土區(qū)列車行駛路基長期永久變形數(shù)值模擬研究[J]．地震工程學報，2014，36(4)：850-856. TIAN Lihui, LING Xianzhang, WANG Lina, et al. Numerical simulation of permanent deformation of embankment induced by train load in permafrost regions along the Qinghai-Tibet Railway[J]. China Earthquake Engineering Journal, 2014, 36(4): 850-856.

[8] 陳拓，吳志堅，林碧蒼，等.青藏高原多年凍土地區(qū)鐵路路基工程動力穩(wěn)定性分析[J].地震工程學報，2014，36(4)：939-945. CHEN Tuo, WU Zhijian, LIN Bicang, et al. Dynamic stability analysis of railway embankment engineering in permafrost regions of the Qinghai-Tibet Plateau[J]. China Earthquake Engineering Journal, 2014, 36(4): 939-945.

[9] 陳拓.動荷載作用下青藏鐵路凍土路基工程動力響應及變形特征分析[D].北京:中國科學院大學,2014. CHEN Tuo. The dynamic response and deformation characteristics of railway embankments under train dynamic load in permafrost regions[D]. Beijing: University of Chinese Academy of Sciences, 2014.

[10] LI Shuangyang, LAI Yuanming, ZHANG Mingyi, et al. Seismic analysis of embankment of Qinghai-Tibet Railway[J]. Cold Regions Science and Technology, 2009, 55: 151-159.

[11] 湯濤，馬濤，黃曉明，等．青藏高速公路寬幅路基溫度場模擬分析[J]．東南大學學報(自然科學版)，2015，45(4)：799-804. TANG Tao, MA Tao, HUANG Xiaoming, et al. Simulation analysis of temperature field of wide subgrade of Qinghai-Tibet Highway[J]. Journal of Southeast University (Natural Science Edition), 2015, 45(4): 799-804.

[12] 俞祁浩，樊凱，錢進，等．我國多年凍土區(qū)高速公路修筑關(guān)鍵問題研究[J]．中國科學：技術(shù)科學，2014，44：425-432. YU Qihao, FAN Kai, QIAN Jin, et al. Key issues of highway construction in permafrost regions in China[J]. Science in China(Technological Sciences), 2014, 44: 425-432.

[13] 馮子亮，盛煜，陳繼，等．青海省共和-玉樹高速公路新建塊石路基下的溫度狀況分析[J]．冰川凍土，2014，36(4)：969-975. FENG Ziliang, SHENG Yu, CHEN Ji, et al. Analyzing the temperature regime within the block stone embankment of the newly constructed gonghe-yushu expressway in qinghai province[J]. Journal of Glaciology and Geocryology, 2014, 36(4): 969-975.

[14] 陳佳宇，田亞護，池秀靜.多年凍土區(qū)寬幅公路路基的溫度特性模擬分析[J]．路基工程，2015(3)：22-29. CHEN Jiayu, TIAN Yahu, CHI Xiujing. Simulation analysis on thermal characteristics of large-width subgrade in permafrost regions[J]. Subgrade Engineering, 2015(3): 22-29.

[15] 劉志強，辛建，喻文兵.寒區(qū)高等級公路寬幅路基的隨機溫度場[J]．冰川凍土，2013，35(6)：1499-1504. LIU Zhiqiang, XIN Jian, YU Wenbing. Random temperature fields of high-level highway wide embankment in cold regions[J]. Journal of Glaciology and Geocryology, 2013, 35(6): 1499-1504.

[16] 符進，李俊，唐曉星，等. 高海拔高寒地區(qū)高速公路建設(shè)技術(shù)研究試驗示范工程的選擇[J]．中外公路，2016，36(1)：6-11. FU Jin, LI Jun, TANG Xiaoxing, et al.The choose of experimental example projects of study on construction technology for expressway in high altitude and cold regions[J]. Journal of China & Foreign Highway, 2016, 36(1): 6-11.

[17] 交通部第二公路勘察設(shè)計院.路基[M].北京：人民交通出版社，1996. The Second Highway Survey and Design Institute of the Ministry of Communications. Roadbed [M]. Beijing: China Communications Press, 1996.

[18] JTGD50—2006，公路瀝青路面設(shè)計規(guī)范[S]. JTGD50—2006，Specifications for design of highway asphalt pavement[S].

Seismic dynamic response analysis of embankment along the Qinghai-Tibet Expressway

ZHAO Tao1,2, WU Zhijian2, LIANG Qingguo1, ZHAO Wenchen3, CHEN Tuo4, XU Shimin2

(1. Key Laboratory of Road & Bridge and Underground Engineering of Gansu Province, Lanzhou Jiaotong University, Lanzhou 730070, China;2. Key Laboratory of Loess Earthquake Engineering, Lanzhou Institute of Seismology, CEA, Lanzhou 730000, China;3. CSCEC AECOM Consultants Co., Ltd, Lanzhou 730000, China;4. The State Key Laboratory for GeoMechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou 221116, China)

The artificial earthquake waves with the exceeding probability for 63.5%, 10% and 2% in 50 years related to this areas were established to study the dynamic response characteristics of proposed Qinghai-Tibet Expressway under seismic loading. Based on numerical simulation, and the two-dimensional nonlinear dynamic finite element method, cross sections of integral and separated embankment were respectively established for the analysis of the dynamic response characteristics under different exceeding probability of seismic wave. The results of this study are summarized as follows. Acceleration amplification coefficients of separated type embankment at the center of both embankment surface and natural ground surface are less than that of integral type embankment. However, horizontal displacement and vertical displacement of separated type embankment are greater than that of integral type embankment. Meanwhile, the shear stress of both integral and separated embankment around two sides are concentrated under the slope feet. With the increase of peak acceleration of seismic waves, horizontal peak displacement on the surface of the embankment increases as the law of linear trend, and the vertical peak displacement increases as the law of exponentially function. These results could provide references and bases for seismic design of expressway in permafrost regions.

the Qinghai-Tibet Expressway; seismic wave; typical embankment; dynamic response

2016-08-25

國家自然科學基金資助項目(41562013,41472297)；甘肅省青年科技基金計劃資助項目(1606RJYA239);甘肅省基礎(chǔ)研究創(chuàng)新群體資助項目(145RJIA332)

吳志堅(1974-)，男，安徽安慶人，研究員，博士，從事巖土地震工程與凍土工程領(lǐng)域的研究；E-mail：zhijianlz@163.com

U416.1

A

1672-7029(2016)12-2381-07