• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    傾斜煤系地層大斷面客專隧道大變形原因分析及處置

    2024-05-23 09:22:20康寶祥宋戰(zhàn)平張慶賀許曉靜周冠南潘紅偉
    土木建筑與環(huán)境工程 2024年2期
    關(guān)鍵詞:隧道工程數(shù)值模擬

    康寶祥 宋戰(zhàn)平 張慶賀 許曉靜 周冠南 潘紅偉

    DOI:?10.11835/j.issn.2096-6717.2022.001

    收稿日期:2021?08?16

    基金項目:國家自然科學(xué)基金(52178393);陜西省創(chuàng)新能力支撐計劃-創(chuàng)新團隊(2020TD-005).

    作者簡介:康寶祥(1989-?),男,主要從事隧道與地下工程研究,E-mail: kangbx@163.com。

    通信作者:宋戰(zhàn)平(通信作者),男,博士,教授,E-mail: songzhpyt@xauat.edu.cn。

    Received: 2021?08?16

    Foundation items: National Natural Science Foundation of China (No. 52178393); Innovation Capability Support Plan of Shaanxi Province Innovation Team (No. 2020TD-005)

    Author brief: KANG Baoxiang (1989-?), main research interests: tunneling and underground engineering, E-mail: kangbx@163.com.

    corresponding author:SONG Zhanping (corresponding author), PhD, professor, E-mail: songzhpyt@xauat.edu.cn.

    摘要:針對滬昆鐵路劉家莊隧道穿越煤系地層施工中發(fā)生的大變形現(xiàn)象,基于理論分析、數(shù)值模擬和現(xiàn)場監(jiān)測方法,分析隧道發(fā)生大變形的原因,推導(dǎo)并驗證了產(chǎn)生大變形的起始位置,得到了掌子面上方煤層單元體沿傾斜面方向的應(yīng)力狀態(tài)與隧道進入含傾斜煤系地層水平距離之間的變化規(guī)律,對傾斜煤層單元體與初期支護結(jié)構(gòu)的應(yīng)力與變形進行了分析。結(jié)果表明:掌子面上方煤體單元應(yīng)力狀態(tài)變化隨掌子面進入傾斜煤層下方水平距離的增加,圍巖經(jīng)歷擠壓、壓剪和剪切滑移3個變形階段。圍巖大變形與初期支護結(jié)構(gòu)破壞均發(fā)生在剪切滑移階段的初期,應(yīng)在壓剪變形階段結(jié)合現(xiàn)場監(jiān)測數(shù)據(jù)對掌子面進行噴射混凝土封閉、注漿加固掌子面上方松散煤層和加強初期支護結(jié)構(gòu)剛度等措施以預(yù)防隧道產(chǎn)生大變形。

    關(guān)鍵詞:隧道工程;煤系地層;圍巖大變形;破壞機理;數(shù)值模擬

    中圖分類號:U455.7 ????文獻標志碼:A ????文章編號:2096-6717(2024)02-0089-11

    Causal analysis and disposal on large deformation of large section passenger dedicated line tunnel in inclined coal strata

    KANG Baoxiang1a,1b?SONG Zhanping1a,1b,1c?ZHANG Qinghe2,?XU Xiaojing1a,1b?ZHOU Guannan1b,3?PAN Hongwei4

    (1a. School of Civil Engineering;?1b. Key Laboratory of Rock and Underground Space Engineering;?1c. Institute of Tunnel and Underground Structure Engineering, Xi,an University of Architecture and Technology, Xi,an 710055, P. R. China;?2. School of Architectural Engineering, Xuzhou College of Industrial Technology, Xuzhou 221000, Jiangsu, P. R. China;?3. China Railway Construction Bridge Engineering Group Co., Ltd., Tianjin 300300, P. R. China;?4. China Railway Beijing Engineering Bureau Group First Engineering Co., Ltd., Xi,an 710100, P. R. China)

    Abstract: For Liujiazhuang railway tunnel of Shanghai-Kunming through the coal strata in construction of large deformation phenomenon, based on theoretical analysis, numerical simulation and field monitoring method, this study investigates the causes of significant deformation in tunnels, focusing on the initial locations of this deformation. It examines the impact of stress states on unit cells aligned with the inclined direction of a coal seam and explores the relationship between tunnel depth and horizontal distance in tilted strata. Additionally, the research analyses the stress and deformation experienced by the element body and the primary support structures in an inclined coal seam. The result showed that the stress state of the coal element above the palm face increases with the horizontal distance between the palm face and the bottom of the inclined coal seam, and the rock mass is subjected to three deformation stages including extrusion, compression shear and shear slip. The large deformation of surrounding rock and the failure of primary support structure both occur in the early stage of shear slip stage, so shotcreting to palm face, grouting reinforcement to coal seam up the palm face and strengthening the stiffness of primary support should be taken to prevent the large deformation of the tunnel in the stage of compression and shear deformation combined with the field monitoring data.

    Keywords: tunnel engineering;?coal strata;?large deformation of surrounding rock;?failure mechanism;?numerical simulation

    隧道圍巖大變形破壞機理及處置技術(shù)的研究一直是隧道與地下工程建設(shè)的難題之一。近年來,隨著隧道和地下工程建設(shè)的發(fā)展,在復(fù)雜地質(zhì)條件下的隧道施工工法比選[1-2]、支護結(jié)構(gòu)的性能分析與優(yōu)化設(shè)計[3-5]以及隧道大變形的原因分析[6-7]與處治措施[8]等方面的研究均取得了一系列的研究成果。在軟弱圍巖施工方法與處治措施方面,李生杰等[9]探討了隧道穿越煤系地層時圍巖的應(yīng)力與變形規(guī)律,提出了圍巖大變形的處治方法。彭鵬等[10]提出了基于綜合評價指標的加固參數(shù)設(shè)計方法并應(yīng)用于實踐,可有效減小含軟弱夾層圍巖的松動區(qū),改善支護結(jié)構(gòu)應(yīng)力集中現(xiàn)象?;魸櫩频?sup>[11]研究了定西黃土隧道圍巖松動區(qū)內(nèi)應(yīng)力與位移的分布規(guī)律,為同類型黃土隧道的安全施工提供了借鑒和指導(dǎo)。宋戰(zhàn)平等[12]基于數(shù)值模擬試驗確定了寶蘭客專小墁坪淺埋偏壓大斷面隧道的施工方案。Huang等[13]通過模型試驗與數(shù)值模擬相結(jié)合的方法,得到了軟弱夾層傾角、厚度、距離、位置等因素與圍巖損傷區(qū)域之間的關(guān)系。馬時強[14]通過數(shù)值模擬分析得到了軟巖大變形條件下圍巖壓力與變形之間的關(guān)系曲線,為大變形讓壓設(shè)計提供依據(jù)。張德華等[15]以阜川隧道支護系統(tǒng)現(xiàn)場試驗為基礎(chǔ),運用數(shù)值模擬方法揭示了軟巖大變形隧道雙層支護的作用機理。陳建勛等[16]對寶漢高速公路連城山隧道不同初期支護結(jié)構(gòu)的現(xiàn)場試驗研究,揭示了圍巖與支護結(jié)構(gòu)的動態(tài)作用機制,提出了有效控制軟巖隧道大變形的支護體系與施工方法。在軟弱圍巖變形機制方面,He等[17]研究了45°傾斜巖層中反梯形隧道的斷裂機理。Qiu等[18]分析了淺埋隧道圍巖的變形機理,提出了控制圍巖變形的有效措施。Li等[19]基于突變理論得到了隧道軟弱圍巖穩(wěn)定性的判別方法。楊忠民等[20]通過研究不同深度隧道開挖的模型試驗,給出了預(yù)防隧道大變形的重點支護部位,明確了大變形處治時的支護措施。以上研究為治理隧道大變形積累了豐富的經(jīng)驗,一些研究通過對特征斷面進行開挖模擬,確定了開挖方法與支護措施;另一些研究通過模型試驗分析了隧道大變形的破壞機理,給出了支護的重點部位,但目前針對隧道開挖過程中軟弱夾煤層受力變化引起圍巖由小變形向大變形發(fā)展的研究較少。

    筆者對滬昆鐵路劉家莊隧道穿越煤系地層圍巖產(chǎn)生大變形災(zāi)害的現(xiàn)象進行分析,通過數(shù)值模擬試驗分析隧道初期支護結(jié)構(gòu)的應(yīng)力與位移變化,建立簡化力學(xué)模型推導(dǎo)大變形起始位置的判據(jù),研究掌子面上方軟弱夾層單元體隨隧道掌子面移動的受力變化過程,分析不同變形階段圍巖的受力與變形特征,對大變形段處置方案與處置效果進行了介紹。

    1 工程概況

    1.1 工程概況及地質(zhì)條件

    劉家莊隧道是位于貴州省境內(nèi)關(guān)嶺—普安區(qū)間的單洞雙線鐵路隧道,線路設(shè)計全長7 583 m,設(shè)計時速350 km/h。隧道出口段穿越東西走向的西隴斷裂帶與傾斜煤層,斷裂帶巖體破碎,地質(zhì)構(gòu)造復(fù)雜,煤層主要位于洞身D1K923+230~D1K924+140段,厚度約2 m,瓦斯壓力為0.85 MPa,地質(zhì)剖面如圖1所示。隧道出口段從大樁號向小樁號掘進,隧道穿越煤層段初期支護完成后產(chǎn)生大變形導(dǎo)致初期支護侵限,掌子面開挖顯示出該段巖體由砂巖、泥巖、頁巖夾煤層構(gòu)成,巖層走向與隧道軸向近平行,節(jié)理裂隙發(fā)育,開挖后部分巖體快速風化為粉末,自穩(wěn)能力極差,掌子面巖性狀況如圖2所示。

    隧道穿越V級圍巖,按V級加強進行設(shè)計。原設(shè)計的初期支護為I25a型鋼拱架,D1K923+905~D1K924+090段采用V[b]型復(fù)合(全封閉)式襯砌,拱部42徑向小導(dǎo)管L=4.0 m,壁厚3.5 mm,環(huán)向間距0.4 m的28根小導(dǎo)管置于I25a型鋼拱架上,拱架縱向間距2.4 m。初期支護采用厚度28 cm的C30氣密性噴射混凝土,預(yù)留變形量15 cm,二次襯砌采用厚度50 cm的C35氣密性鋼筋混凝土。

    1.2 隧道圍巖變形情況

    劉家莊隧道出口段開挖方式為上下臺階法,開挖進尺為2 m,鋼拱架距離掌子面2 m,上下臺階開挖間距2倍的開挖進尺,下臺階與仰拱間距3倍的上臺階開挖進尺。采用原設(shè)計的初期支護方案施工至D1K923+905斷面時,掌子面施工至傾斜煤層附近,強烈的擾動導(dǎo)致巖層中松散破碎的巖塊擠壓錯動,加之地層中富存含瓦斯的大量氣體,在禁止焊接的施工條件下,初期支護結(jié)構(gòu)中的鎖腳錨管與鋼拱架采用預(yù)制環(huán)的連接方式,難以形成牢固的整體受力體系,最終導(dǎo)致初期支護結(jié)構(gòu)在拱肩部位發(fā)生多處較大面積的噴射混凝土剝落破壞,鋼拱架扭曲變形嚴重,節(jié)點部位剛度較大,荷載傳遞過程中在鋼拱架節(jié)點上部產(chǎn)生應(yīng)力集中,導(dǎo)致拱架屈曲破壞,出現(xiàn)傾斜錯動及撕裂破壞現(xiàn)象;部分噴射混凝土沿鋼拱架翼緣脫落,產(chǎn)生環(huán)向開裂現(xiàn)象,如圖3所示。

    從現(xiàn)場監(jiān)測數(shù)據(jù)分析可知,在監(jiān)測時間段內(nèi),斷面D1K923+925進入傾斜煤層水平距離約20 m處,由于該范圍內(nèi)圍巖巖性相對較好,測點變形量最大值為69 mm,且拱頂沉降與凈空收斂逐步趨于穩(wěn)定;斷面D1K923+915進入傾斜煤層水平距離約30 m處,掌子面前方與傾斜煤層水平距離約20 m,因圍巖巖性變化復(fù)雜且?guī)r體相對破碎,測點變形量急劇增加,拱頂沉降與凈空收斂均大于設(shè)計預(yù)留變形量150 mm,穩(wěn)定跡象不明顯。由此可見,該隧道施工中圍巖發(fā)生小變形與大變形的臨界點位于D1K923+925~ D1K923+915內(nèi),該范圍內(nèi)圍巖內(nèi)部的應(yīng)力與支護結(jié)構(gòu)的受力均發(fā)生較為明顯的變化。

    為對比研究,選取D1K923+915為特征斷面,對現(xiàn)場監(jiān)測數(shù)據(jù)進一步分析,得到變形速率-時間關(guān)系曲線如圖4、圖5所示,該斷面前方與傾斜煤層水平距離約20 m,上臺階開挖共計16 d,施工前10 d內(nèi),變形速率較為平緩,測點變形量增長緩慢,10~16 d施工過程中,變形速率急劇增加,最大收斂速率達15 mm/d,上臺階初期支護施工完成后變形速率逐漸減小,測點變形量緩慢增加;下臺階開挖后,變形速率陡增,拱頂沉降與凈空收斂速率最大值分別為10、13 mm/d,測點變形量急劇增加,即使在下臺階初期支護施工完成后,變形速率雖有減小,但測點變形量未呈現(xiàn)出收斂態(tài)勢,在下臺階施工12 d后,凈空收斂值為148 mm,已接近設(shè)計預(yù)留變形量;在經(jīng)歷70 d左右的時間后穩(wěn)定跡象仍不明顯,到隧道侵限處置前拱頂沉降累計值為165 mm,凈空收斂累計值為278 mm,二者均超過設(shè)計預(yù)留變形量150 mm。

    從該斷面變形監(jiān)測中可以發(fā)現(xiàn),上下臺階施工的過程中變形速率變化頻繁、不穩(wěn)定,尤其是變形量急劇增加階段的變形速率變化幅值相對較大,證明了開挖斷面在接近傾斜煤層的施工過程中,由于含煤系地層圍巖內(nèi)部構(gòu)造比一般巖體復(fù)雜,軟弱夾層礦物組成中含有蒙脫石和伊利石,強度相差大,巖體內(nèi)部應(yīng)力變化劇烈,圍巖自穩(wěn)能力差,容易導(dǎo)致初期支護結(jié)構(gòu)產(chǎn)生應(yīng)力集中現(xiàn)象;凈空收斂值大于拱頂沉降值,表明水平應(yīng)力σH大于垂直應(yīng)力σZ,反映出地應(yīng)力以構(gòu)造應(yīng)力為主的特征。

    2 數(shù)值模擬分析

    2.1 模型幾何尺寸

    針對含傾斜煤層的地質(zhì)特征,考慮隧道開挖對圍巖的有效影響范圍,借助FLAC 3D有限差分軟件,模擬劉家莊隧道的施工力學(xué)行為,研究隧道開挖過程中掌子面上方煤層單元體應(yīng)力與位移變化過程及初期支護的力學(xué)性狀,為預(yù)防隧道大變形提供切實可行的理論基礎(chǔ)。數(shù)值模擬的里程樁號范圍為D1K923+945~D1K923+905,隧道凈高12.3 m,凈寬14.9 m,綜合考慮地層因素與邊界約束效應(yīng),模型橫向方向(X軸)由隧道中心線位置向兩側(cè)各取50 m,上邊界(Z軸)至地表,下邊界(Z軸)距隧道底部37 m,沿隧道縱向(Y軸)取110 m,其中煤層厚度為2 m,傾角為55°,煤層起始位置對應(yīng)里程樁號為D1K923+945。隧道與煤層的空間關(guān)系:隧道與傾斜煤層走向交角90°,即正交穿越傾斜煤層,隧道與煤層的空間關(guān)系如圖6所示。模型側(cè)面及底部邊界條件采用固定約束,上表面自由,只考慮巖體的自重應(yīng)力,忽略構(gòu)造應(yīng)力,巖體與初期支護均為實體單元,采用摩爾–庫倫(Mohr-Coulomb)強度準則模擬巖體材料,彈性本構(gòu)(Elastic)模擬初期支護材料。

    2.2 計算參數(shù)

    圍巖的力學(xué)參數(shù)參考所依托隧道工程地質(zhì)勘察資料,不考慮巖層分界面及煤層內(nèi)部的結(jié)構(gòu)面作用,結(jié)合規(guī)范要求,給出圍巖、煤層及初期支護結(jié)構(gòu)的物理力學(xué)參數(shù)見表2。

    2.3 計算結(jié)果及分析

    2.3.1 掌子面正上方煤層單元應(yīng)力分析

    圖7為掌子面開挖至D1K923+905時煤層的應(yīng)力云圖,煤層應(yīng)力整體呈層狀近似均勻分布,兩側(cè)拱肩附近煤層單元體應(yīng)力值達到3.15 MPa,該部位應(yīng)力集中明顯,說明開挖斷面離煤層越近,拱肩附近越容易出現(xiàn)破壞,對該部位巖體應(yīng)進行注漿加固并增加初期支護的剛度;隧道開挖過程中掌子面與煤層水平距離如圖8所示,不同里程段掌子面上方煤層單元體的應(yīng)力與位移變化曲線如圖9所示。

    從圖9給出的不同掌子面上方煤層單元體的應(yīng)力與位移變化曲線可以看出,掌子面正上方煤體單元的最大主應(yīng)力與最大剪應(yīng)力隨著開挖距離的增加呈線性增長,最大主應(yīng)力的增長速率明顯高于最大剪應(yīng)力的增長速率;當掌子面與傾斜煤層起始位置水平距離S1=20 m時,掌子面上方煤體單元的豎向位移增長平緩,累計位移為56 mm,與斷面D1K923+925處實測沉降值較為接近,由此可知,開挖距離小于S1時,煤層的最大正應(yīng)力與最大剪應(yīng)力呈線性增大,其應(yīng)力狀態(tài)均小于極限應(yīng)力,體現(xiàn)了巖體與煤層中以擠壓變形為主的特征。

    當掌子面與傾斜煤層起始位置水平距離S2=30 m時,煤體單元的累計位移變化量74 mm,掌子面介于S1S2之間時,其上方煤體單元的豎向位移增長速率逐漸變大,煤體單元所受的最大正應(yīng)力與最大剪應(yīng)力繼續(xù)增大,壓力作用導(dǎo)致內(nèi)部裂隙閉合、能量積蓄持續(xù)增加。當掌子面與傾斜煤層起始位置水平距離超過S2=30 m時,掌子面上方煤體單元的豎向位移增長速率激增,當S3=40 m時,煤體單元的累計位移增量為244 mm。

    圖9還表明,煤層最大主應(yīng)力、最大剪應(yīng)力在S=38 m處出現(xiàn)突變,前者逐漸趨于平緩,后者逐漸趨于零。因此,當掌子面逐漸靠近煤層時,趨于極限受壓狀態(tài)下的煤體單元在開挖卸荷過程中最大剪應(yīng)力逐漸接近其極限抗剪強度,由于圍巖等級較差、開挖過程中掌子面爆破對軟弱圍巖的擾動、初期支護施作不及時等因素,導(dǎo)致掌子面附近約束不足,圍巖擠出變形大,處于極限受力狀態(tài)下的煤體單元抗剪強度達到極限值,該狀態(tài)下的煤體單元因壓潰導(dǎo)致體積膨脹,破碎的塊體沿巖體破裂面產(chǎn)生剪切滑移現(xiàn)象,該階段傾斜煤層以剪切滑移變形為主。

    通過分析掌子面正上方煤體單元應(yīng)力與位移曲線圖,可以得到含傾斜煤層的隧道在開挖過程中煤體單元的應(yīng)力狀態(tài)主要經(jīng)歷了受壓、壓剪和剪切3個變化階段;單元體的豎向位移變化主要經(jīng)歷了變形速率較小的平緩增長階段(小變形階段)、變形速率逐漸增大的過渡變形階段和變形速率激增的大變形階段;煤層單元體的破壞經(jīng)歷了從壓應(yīng)力為主的裂隙閉合過程到剪應(yīng)力為主的剪切滑移過程,并伴隨著煤體單元壓剪應(yīng)力狀態(tài)下的體積膨脹現(xiàn)象。

    2.3.2 隧道初期支護受力分析

    掌子面隨著開挖里程的變化,隧道初期支護不同部位的最大主應(yīng)力、最小主應(yīng)力與最大剪應(yīng)力沿隧道縱向分布如圖10所示。

    隨著開挖距離的增加,初期支護結(jié)構(gòu)拱頂、拱肩和拱腳處的最大主應(yīng)力呈現(xiàn)出先緩后陡的增長趨勢,最大主應(yīng)力均為壓應(yīng)力,且拱頂受到的壓應(yīng)力較大,拱肩與拱腳受到的壓應(yīng)力值較??;初期支護各部位最大主應(yīng)力在掌子面開挖至S1=20 m時出現(xiàn)急劇增長趨勢,拱頂與拱肩處最大壓應(yīng)力值均達到1.87 MPa,拱腳處的最大主應(yīng)力值為1.67 MPa;在S2附近拱肩與拱腳的最大主應(yīng)力值波動較大,說明該斷面附近巖層應(yīng)力變化較為復(fù)雜;初期支護受到的最大壓應(yīng)力值為5.27 MPa,未超過C30混凝土的抗壓強度設(shè)計值,支護結(jié)構(gòu)未產(chǎn)生壓裂破壞。此時拱肩與拱腳最大壓應(yīng)力值分別為4.01、3.50 MPa。

    隨著開挖距離的增加,初期支護拱頂、拱肩和拱腳處的最小主應(yīng)力呈現(xiàn)出先錯峰遞增后快速減小,再逐漸增大的趨勢,拱頂與拱肩的最小主應(yīng)力由壓應(yīng)力轉(zhuǎn)變?yōu)槔瓚?yīng)力,且拱頂處拉應(yīng)力最大,拱腳最小主應(yīng)力為壓應(yīng)力;初期支護各部位最小主應(yīng)力在掌子面開挖至S1=20 m后變化明顯,拱頂與拱肩最小主應(yīng)力方向改變,應(yīng)力變化幅值分別達到1.80、1.50 MPa;最小主應(yīng)力突變說明該斷面開挖后巖層內(nèi)部應(yīng)力變化劇烈,對初期支護結(jié)構(gòu)產(chǎn)生較大的沖擊作用,致使初期支護拱頂與拱肩部位由受壓狀態(tài)轉(zhuǎn)變?yōu)槭芾瓲顟B(tài);此時支護結(jié)構(gòu)上的最大拉應(yīng)力位于拱頂,其值為1.30 MPa,未超過C30混凝土的抗拉強度設(shè)計值,結(jié)構(gòu)未產(chǎn)生壓裂破壞。

    隨著開挖距離的增加,初期支護結(jié)構(gòu)拱頂、拱肩和拱腳處的最大剪應(yīng)力呈現(xiàn)出先緩后快的增長趨勢,且拱頂與拱肩處的最大剪應(yīng)力方向向下,拱腳處所受的剪應(yīng)力方向向上;掌子面開挖至S1=20 m時,初期支護拱頂、拱肩與拱腳處的剪應(yīng)力值分別為-0.69、-0.59、0.57 MPa,增長速率平緩;掌子面開挖至S1=20 m后,拱頂、拱肩與拱腳處的最大剪應(yīng)力快速增長,其中拱頂增長速率最大、拱肩次之、拱腳最??;當掌子面開挖至S2=30 m時,拱頂最大剪應(yīng)力值為2.86 MPa,超過C30混凝土抗剪強度設(shè)計值2.85 MPa,支護結(jié)構(gòu)在拱頂處產(chǎn)生剪切破壞現(xiàn)象。

    2.3.3 隧道初期支護變形分析

    圖11給出了隧道初期支護不同部位的位移監(jiān)測曲線,從圖11中可以看出,隨著開挖距離的增加,初期支護拱頂沉降、拱肩變形和拱腳的收斂變形呈現(xiàn)出先平緩后激增,再緩慢的增長趨勢,位移激增階段位于S1=20 m與S2=30 m之間。

    表3給出了該段隧道數(shù)值計算結(jié)果與現(xiàn)場實測結(jié)果的對比分析結(jié)果。由表3可見,拱頂沉降最大值為278 mm,凈空收斂最大值為213 mm;當掌子面開挖至S=26 m時,拱頂沉降值已經(jīng)超過隧道設(shè)計預(yù)留變形量150 mm;當掌子面開挖至S2=30 m時,拱肩與拱腳收斂值均超過設(shè)計預(yù)留變形量,初期支護產(chǎn)生侵限變形,隧道發(fā)生大變形破壞,因此應(yīng)在掌子面開挖至S1=20 m處開始采取相應(yīng)的施工措施,防止隧道繼續(xù)開挖引起初期支護破壞。

    由于數(shù)值模擬的圍巖屬于均勻介質(zhì)且地層關(guān)于隧道中軸線對稱,忽略了實際巖層中的瓦斯壓力、節(jié)理裂隙、地下水以及施工過程中爆破振動等因素的影響,導(dǎo)致數(shù)值計算結(jié)果與現(xiàn)場實測值存在一定差異,但數(shù)值計算得到的位移分布規(guī)律與監(jiān)測結(jié)果趨勢一致,表明數(shù)值模擬基本可以反映實際隧道圍巖和初期支護的應(yīng)力、變形與破壞規(guī)律,可為預(yù)防隧道大變形災(zāi)害提供理論指導(dǎo)和參考。

    綜上分析,隧道掌子面開挖距離小于S1=20 m時,初期支護各部位的最大主應(yīng)力值、最小主應(yīng)力值及最大剪應(yīng)力值均小于C30混凝土強度設(shè)計值,支護結(jié)構(gòu)各部位的位移變化量均小于設(shè)計預(yù)留變形量,初期支護未發(fā)生破壞;當掌子面開挖距離大于S2=30 m時,拱頂處的最大剪應(yīng)力值達到了C30混凝土的抗剪強度設(shè)計值,拱頂沉降值與拱肩收斂值超過隧道設(shè)計預(yù)留變形量,初期支護發(fā)生剪切破壞,導(dǎo)致隧道出現(xiàn)變形侵限現(xiàn)象。

    3 圍巖大變形機理分析及處置

    3.1 變形產(chǎn)生原因分析

    劉家莊隧道軟弱圍巖大變形破壞的主要原因有如下幾種:

    1)地質(zhì)原因,該段洞身穿越二疊系龍?zhí)督M(P2l),龍?zhí)督M(P2l)地層含多層煤,屬高瓦斯地層,瓦斯壓力為0.85 MPa;該段由泥巖、砂巖夾煤層構(gòu)成,圍巖屬于整體較破碎的強風化軟巖。D1K923+905掌子面附近自穩(wěn)能力較差的破碎巖塊對拱腳基礎(chǔ)的支撐強度不足,導(dǎo)致拱架不穩(wěn)定而引起拱頂下沉。煤系地層巖體內(nèi)存在光滑的節(jié)理面,巖層之間有夾泥現(xiàn)象,相對較弱的膠結(jié)作用導(dǎo)致開挖后的圍巖極易沿薄弱部位產(chǎn)生滑移,加快了變形破壞。

    2)施工原因,該里程段穿越煤層,包含瓦斯在內(nèi)的大量氣體儲存在圍巖中,掌子面爆破時產(chǎn)生的沖擊作用,降低了圍巖的整體性,尤其是拱頂及線路方向右側(cè)拱肩處為極度破碎的塊狀煤矸石,夾雜有大量粉末,現(xiàn)場渣樣確定為煤。由于該段隧道瓦斯?jié)舛雀?,施工時采用預(yù)制環(huán)的方式連接鋼拱架與鎖腳錨管,與剛結(jié)點相比,鉸接點的作用較弱,降低了初期支護的承載能力,受爆破過程中產(chǎn)生的沖擊作用,初期支護產(chǎn)生嚴重變形。

    3.2 圍巖大變形的位置預(yù)測與力學(xué)分析

    3.2.1 大變形起始位置的判據(jù)推導(dǎo)

    為了分析傾斜煤系地層中軟弱夾層對隧道變形的影響,僅考慮左右對稱的地層狀況,建立簡化模型如圖12(a)所示,將隧道穿越的巖體分為3層,中間層②為煤層,煤層兩側(cè)的①層為相對完整的泥巖和砂巖層。從應(yīng)力傳遞的角度對掌子面上部煤層單位厚度的薄片單元進行受力分析,如圖12(b)所示。當薄片單元在豎直方向所受的摩擦力與支護壓力的合力小于重力時,煤層單元產(chǎn)生滑動現(xiàn)象,由此可定義隧道大變形的起始位置函數(shù)

    式中:S為隧道掌子面進入傾斜煤層下方的水平距離;Z為單元體的埋深;γ、c、分別為煤層單元體的重度、黏聚力和內(nèi)摩擦角;α為煤層傾角;T為初期支護的支護壓力。

    式中:σH為水平主應(yīng)力;K0為側(cè)壓力系數(shù)。

    基于劉家莊隧道地勘資料中煤層的力學(xué)參數(shù):煤層傾角α=55°,煤層重度γ=15 kN/m3,煤層黏聚力c=52 kPa,煤層內(nèi)摩擦角?=25.8°,側(cè)壓力系數(shù)K0=0.48,根據(jù)式(7)和式(8)得到隧道產(chǎn)生大變形的起始位置為S1=20.9 m,即隧道掌子面進入傾斜煤層的水平距離超過20.9 m后會產(chǎn)生大變形現(xiàn)象,與實際產(chǎn)生大變形的起始區(qū)間(D1K923+925~ D1K923+915)吻合。

    3.2.2 圍巖變形過程的受力分析

    以材料力學(xué)[21]為基礎(chǔ),分析小變形條件下的掌子面上方煤塊單元體的應(yīng)力狀態(tài)變化過程,揭示軟弱夾層地質(zhì)條件下隧道開挖過程中產(chǎn)生大變形現(xiàn)象的作用機理。取掌子面正上方②層煤體單元為研究對象進行受力分析,巖層分界面與煤層內(nèi)部的單元體受力狀態(tài)如圖13所示。巖體中初始豎向應(yīng)力及水平應(yīng)力分別為σZσH。

    根據(jù)簡化力學(xué)模型,以傾斜面的法線ζ和切線η為參考軸建立直角坐標系分析單元體的應(yīng)力平衡狀態(tài),可得

    式中:σS為傾斜面上煤體單元所受的正應(yīng)力;τS為傾斜面上體單元所受的剪應(yīng)力;β為結(jié)構(gòu)面傾角。

    式(10)為平面應(yīng)力狀態(tài)下,隧道開挖過程中掌子面上方煤層單元體沿傾斜面方向的正應(yīng)力與剪應(yīng)力(應(yīng)力狀態(tài)),反映了隧道進入含傾斜煤系地層水平距離與掌子面上方煤體單元應(yīng)力狀態(tài)之間的變化規(guī)律。結(jié)合摩爾庫倫強度準則得到煤體單元應(yīng)力狀態(tài)與水平距離S之間的關(guān)系曲線如圖14所示。

    3.2.3 擠壓滑移的形成

    隧道開挖前巖層中初始地應(yīng)力處于平衡狀態(tài),隧道開挖過程中周圍巖層處于卸載狀態(tài),隨著掌子面與傾斜煤層水平距離的變化,煤層內(nèi)部及邊界處應(yīng)力傳遞比較復(fù)雜,結(jié)合煤體單元應(yīng)力狀態(tài)與強度包線之間的關(guān)系曲線(圖14),將圍巖變形分為3個階段,即擠壓變形階段(0<S<S1)、壓剪變形階段(S1S<S2)和剪切滑移階段(S2S);變形階段分區(qū)如圖15所示。掌子面施工至S1前,隧道開挖對掌子面上方與前方傾斜煤層的擾動較小,開挖段上方楔形巖體在自重作用下存在向掌子面方向移動的趨勢,煤層受兩側(cè)巖體的擠壓作用,煤體單元沿傾斜面的正應(yīng)力與剪應(yīng)力增大,但剪應(yīng)力值均小于結(jié)構(gòu)面與煤體單元抗剪強度極限值,即τ<τf1τ<τf2,單元體未產(chǎn)生剪切破壞,巖層處于受壓狀態(tài),巖體以小變形為主,該階段為擠壓變形階段。

    掌子面施工至S1S2之間時,煤層單元體沿傾斜面的正應(yīng)力與剪應(yīng)力同時增大,煤體單元剪應(yīng)力值介于結(jié)構(gòu)面(巖層分界面與煤層內(nèi)部)與煤體抗剪強度極限值之間,即τf1<τ<τf2,受壓狀態(tài)下的單元體出現(xiàn)剪切破壞,節(jié)理裂隙在擠壓閉合的過程中伴隨有剪切變形,巖層處于極限抗剪條件下的受壓狀態(tài),塑性區(qū)率先在煤層中產(chǎn)生并隨著S的增加逐漸擴大,煤體單元所受壓應(yīng)力逐漸趨于抗壓強度極限值;由于初期支護的支撐作用及臺階土體對掌子面擠出變形的反壓作用,煤層整體沿傾斜面保持相對穩(wěn)定的狀態(tài),但圍巖中未釋放的形變能在巖層內(nèi)部逐漸聚集,巖體以延性變形為主,該階段為壓剪變形階段。

    隧道開挖至S2后,隨著掌子面與傾斜煤層距離的減小,其上方單元體均達到抗剪強度極限值,即τ=τf2;當煤體單元所受壓應(yīng)力達到抗壓強度極限值時,極限狀態(tài)下的壓剪作用導(dǎo)致單元體產(chǎn)生碎脹力與碎脹變形,由于臺階土體對掌子面的反壓作用減弱,施工對圍巖的強烈擾動與煤層內(nèi)形變能的釋放引起掌子面附近巖體的擠出變形過大,極限狀態(tài)下的煤體單元因壓潰而碎散,巖層分界面上的摩阻力因接觸面減小而減弱,碎散的巖塊沿傾斜面向下產(chǎn)生剪切滑移現(xiàn)象,最終導(dǎo)致掌子面附近圍巖失穩(wěn),破壞瞬間釋放的能量(形變能與高壓力瓦斯共同作用)對初期支護產(chǎn)生強烈的沖擊效應(yīng),致使作用在初期支護上的荷載激增,拱架之間因縱向抗剪強度不足產(chǎn)生沖切破壞,初期支護錯動為兩段,出現(xiàn)整體下沉;拱架內(nèi)部連接處由剛度不同導(dǎo)致應(yīng)力集中,鋼拱架扭曲變形位于連接板上部;巖體由延性變形轉(zhuǎn)變?yōu)榇嘈云茐?,圍巖產(chǎn)生大變形導(dǎo)致初期支護變形侵限,該階段為剪切滑移階段。

    綜上所述,通過分析隧道掌子面上方煤體單元應(yīng)力狀態(tài)變化的過程,根據(jù)掌子面進入傾斜煤層下方的水平距離將隧道施工變形分為3個階段,利用地質(zhì)勘察資料提供的巖土參數(shù)可初步估算位于傾斜煤層下方的隧道開挖時圍巖不同變形階段的起始距離;結(jié)合現(xiàn)場監(jiān)測數(shù)據(jù),在應(yīng)力變化較為復(fù)雜的壓剪變形階段通過調(diào)整初期支護的設(shè)計參數(shù)與開挖工法,避免支護結(jié)構(gòu)破壞或隧道產(chǎn)生大變形現(xiàn)象。式(10)可作為含傾斜煤系地層隧道施工中劃分變形階段的參考依據(jù),由于巖體結(jié)構(gòu)面上的強度參數(shù)精確值較難獲取,可根據(jù)煤巖的強度參數(shù)利用式(7)計算壓剪變形的起始距離,為預(yù)防隧道產(chǎn)生大變形的部位提供理論參考。

    3.3 變形處置方案及效果

    針對煤系地層剪切滑移導(dǎo)致的大變形特征,采取的相應(yīng)措施如下:

    1)對上臺階掌子面進行噴射混凝土封閉處理,阻止圍巖的進一步風化;同時采用渣土逐層夯填至開挖節(jié)點位置,抑制變形發(fā)展速率,并為侵限處理提供施工作業(yè)面;在施工前在大變形段上部拱架間采用濕式鉆孔方式設(shè)置3排引排初期支護結(jié)構(gòu)背部氣體的瓦斯溢出孔,確保施工安全。

    2)對掌子面上方的傾斜煤層采取深部主動支護原則,利用注漿加固方式對煤層內(nèi)部的松動巖塊進行加固,加強煤層松散塊體及巖層分界面上的黏聚力,增強圍巖的自穩(wěn)能力。上臺階增設(shè)42徑向小導(dǎo)管注漿加固圍巖,管長6 m,壁厚3.5 mm,縱、橫向間距1 m,端部1 m以下按梅花形鉆設(shè)孔徑6~8 mm、間距10~20 mm的注漿孔,壓注水泥漿。

    3)侵限段按Vd型全封閉復(fù)合襯砌換拱,根據(jù)支護特性原理,及時施加大剛度支護封閉的隧道結(jié)構(gòu),加密型鋼拱架間距,及時施作仰拱及二次襯砌,將圍巖變形控制在擠壓變形階段。

    4)換拱遵循“弱爆破、逐圈換拱、勤測量”原則,減少爆破裂隙的產(chǎn)生;換拱順序為先上后下,拆除噴射混凝土及變形鋼架,用I25a工字鋼替換原有變形拱架,并用錨桿固定,使新設(shè)拱架與未變形段拱架連接可靠;換拱過程中設(shè)置臨時支撐,并確保臺架與臨時支撐底部連接牢固,其后噴射混凝土完成初期支護。侵限段換拱如圖16所示。

    經(jīng)過噴射混凝土封閉掌子面,注漿加固松散巖體,回填渣土形成反壓平臺,更換變形的鋼拱架以及及時施作仰拱與二次襯砌后,拱頂沉降與凈空收斂監(jiān)測數(shù)據(jù)如圖17所示,在拆換時,由于圍巖失去反作用力,導(dǎo)致鋼拱架變形速率較大,在鋼拱架封閉與噴射混凝土完成后,隧道變形逐漸趨于穩(wěn)定,圍巖大變形經(jīng)處置后初期支護變形總量較小,趨于穩(wěn)定的時間較短,變形速率趨近于零,說明該處置措施能夠有效加固含傾斜煤層的圍巖,阻止圍巖持續(xù)變形。

    4 結(jié)論

    以滬昆鐵路劉家莊隧道為背景,針對隧道穿越煤系地層發(fā)生的大變形侵限現(xiàn)象,通過數(shù)值模擬與簡化力學(xué)模型,分析煤層單元體與初期支護結(jié)構(gòu)的受力和變形的變化,得到以下結(jié)論:

    1)煤系地層屬典型的軟弱圍巖,其力學(xué)性質(zhì)差,且含較大的氣體壓力,高瓦斯軟弱圍巖地層隧道初期支護過程中不允許焊接致使鋼拱架與鎖腳錨管不能形成整體傳力構(gòu)件,降低了初期支護結(jié)構(gòu)的整體承載能力,加之掌子面爆破過程中產(chǎn)生的沖擊作用,導(dǎo)致初期支護產(chǎn)生了嚴重變形。

    2)在假設(shè)地層均勻連續(xù)的條件下,基于應(yīng)力傳遞思想,給出了傾斜煤系地層中隧道開挖引起圍巖大變形的起始位置的判據(jù),并通過計算驗證了該判據(jù)的有效性;通過煤體單元的豎向位移變化與巖體力學(xué)參數(shù)的變化將軟弱圍巖的破壞過程分為擠壓變形階段(0<S<S1),壓剪變形階段(S1S<S2)和剪切滑移階段(S2S)。

    3)模擬隧道開挖過程中掌子面上方煤體單元與初期支護的應(yīng)力與位移變化,當掌子面開挖距離小于S1時,初期支護各部位的最大主應(yīng)力值、最小主應(yīng)力值及最大剪應(yīng)力值均小于C30混凝土強度設(shè)計值,各部位的位移變化量均小于設(shè)計預(yù)留變形量,初期支護不發(fā)生破壞;當掌子面開挖距離大于S2時,拱頂處的最大剪應(yīng)力值達到了C30混凝土抗剪強度設(shè)計值,拱頂沉降與拱肩收斂值超過隧道設(shè)計預(yù)留變形量,初期支護發(fā)生剪切破壞。

    4)通過噴射混凝土封閉掌子面、注漿加固松散巖體和堆載反壓的措施加強掌子面附近圍巖的自穩(wěn)能力,適時更換侵限段拱架,及時施作仰拱和二次襯砌等能有效控制穿越含傾斜煤系地層隧道的圍巖大變形。

    參考文獻

    [1] ?郗舉科, 吳軍, 宋戰(zhàn)平. 軟弱破碎帶大斷面隧道的開挖與支護技術(shù)[J]. 西北水力發(fā)電, 2003, 19(2): 36-39.

    XI J K, WU J, SONG Z P. The technology of excavation and set supporting for tunnel with large cross section in katamorphic zone [J]. Journal of Northwest Hydroelectric Power, 2003, 19(2): 36-39. (in Chinese)

    [2] ?牛澤林, 劉朝鵬, 鄭蕾, 等. 大斷面公路板巖隧道開挖與支護數(shù)值模擬研究[J]. 西安建筑科技大學(xué)學(xué)報(自然科學(xué)版), 2020, 52(3): 376-383, 403.

    NIU Z L, LIU C P, ZHENG L, et al. Numerical simulation study on excavation and support of slab tunnel in large section highway [J]. Journal of Xi,an University of Architecture & Technology (Natural Science Edition), 2020, 52(3): 376-383, 403. (in Chinese)

    [3] ?霍潤科, 王艷波, 宋戰(zhàn)平, 等. 黃土隧道初期支護性能分析[J]. 巖土力學(xué), 2009, 30(Sup2): 287-290.

    HUO R K, WANG Y B, SONG Z P, et al. Analysis of initial support properties for loess tunnel [J]. Rock and Soil Mechanics, 2009, 30(Sup2): 287-290. (in Chinese)

    [4] ?ZUO J P, WANG R K, WU A M, et al. Optimization support controlling large deformation of tunnel in deep mine based on discontinuous deformation analysis [J]. Procedia Environmental Sciences, 2012, 12(1): 1045-1054.

    [5] ?宋戰(zhàn)平. 基于位移反分析的將軍溝隧道結(jié)構(gòu)穩(wěn)定性評價技術(shù)[J]. 公路交通科技(應(yīng)用技術(shù)版), 2011, 7(9): 144-147.

    SONG Z P. Analysis structural stability evaluation technology based on displacement reverse of Jiangjungou tunnel [J]. Journal of Highway and Transportation Research and Development (Applied Technology), 2011, 7(9): 144-147. (in Chinese)

    [6] ?CAO C Y, SHI C H, LEI M F, et al. Squeezing failure of tunnels: a case study [J]. Tunnelling and Underground Space Technology, 2018, 77(2): 188-203.

    [7] ?ZHAO Y D, ZHANG Y X, YANG J S, et al. Geotechnical monitoring and analyses on the stability of highway tunnel constructed in carbonaceous slate stratum [J]. Geotechnical and Geological Engineering, 2019, 37(5): 4613-4625.

    [8] ?YANG Z M, WU S C, GAO Y T, et al. Time and technique of rehabilitation for large deformation of tunnels in jointed rock masses based on FDM and DEM numerical modeling [J]. Tunnelling and Underground Space Technology, 2018, 81: 669-681.

    [9] ?李生杰, 謝永利, 吳丹澤, 等. 穿越煤系地層隧道圍巖大變形機制及處治研究[J]. 巖石力學(xué)與工程學(xué)報, 2013, 32(Sup2): 3501-3508.

    LI S J, XIE Y L, WU D Z, et al. Research on mechanism of large deformation and countermeasures of surrounding rocks of tunnel through coal strata [J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(Sup2): 3501-3508. (in Chinese)

    [10] ?彭鵬, 張頂立, 孫振宇. 含軟弱夾層隧道圍巖變形特性與加固參數(shù)設(shè)計方法[J]. 巖石力學(xué)與工程學(xué)報, 2021, 40(11): 2260-2272.

    PENG P, ZHANG D L, SUN Z Y. Deformation characteristics of surrounding rock and reinforcement parameter design of weak interlayer tunnels [J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(11): 2260-2272. (in Chinese)

    [11] ?霍潤科, 黃沛, 宋戰(zhàn)平, 等. 隧道施工對黃土松動區(qū)內(nèi)力和位移影響研究[J]. 巖土力學(xué), 2009, 30(Sup2): 273-276.

    HUO R K, HUANG P, SONG Z P, et al. Research on effect of tunnel construction on internal force and deformation at loosed loess area [J]. Rock and Soil Mechanics, 2009, 30(Sup2): 273-276. (in Chinese)

    [12] ?宋戰(zhàn)平, 王童, 周建軍, 等. 淺埋偏壓大斷面隧道施工優(yōu)化及受力特征分析[J]. 地下空間與工程學(xué)報, 2017, 13(2): 459-468.

    SONG Z P, WANG T, ZHOU J J, et al. Analysis of construction optimization and mechanical characteristic of shallow large section tunnel [J]. Chinese Journal of Underground Space and Engineering, 2017, 13(2): 459-468. (in Chinese)

    [13] ?HUANG F, ZHU H H, XU Q W, et al. The effect of weak interlayer on the failure pattern of rock mass around tunnel -?Scaled model tests and numerical analysis [J]. Tunnelling and Underground Space Technology, 2013, 35(3): 207-218.

    [14] ?馬時強. 軟巖大變形條件下圍巖變形與圍巖壓力的關(guān)系研究[J]. 現(xiàn)代隧道技術(shù), 2020, 57(1): 44-50.

    MA S Q. On relationship between surrounding rock deformation and surrounding rock pressure under large deformation of soft rock [J]. Modern Tunnelling Technology, 2020, 57(1): 44-50. (in Chinese)

    [15] ?張德華, 雷可, 譚忠盛, 等. 軟巖大變形隧道雙層初期支護承載性能對比試驗研究[J]. 土木工程學(xué)報, 2017, 50(Sup2): 86-92.

    ZHANG D H, LEI K, TAN Z S, et al. Field test on bearing and deformation properties of double primary layers applied to large deformation tunnel in weak rock [J]. China Civil Engineering Journal, 2017, 50(Sup2): 86-92. (in Chinese)

    [16] ?陳建勛, 劉偉偉, 陳麗俊, 等. 綠泥石片巖地層大跨度公路隧道大變形控制及合理支護形式現(xiàn)場試驗[J]. 中國公路學(xué)報, 2020, 33(12): 212-223.

    CHEN J X, LIU W W, CHEN L J, et al. In-situ experimental study on large-deformation control and reasonable support forms for a large-span highway tunnel in chlorite schist [J]. China Journal of Highway and Transport, 2020, 33(12): 212-223. (in Chinese)

    [17] ?HE M C, PENG Y Y, ZHAO S Y, et al. Fracture mechanism of inversed trapezoidal shaped tunnel excavated in 45° inclined rock strata [J]. International Journal of Mining Science and Technology, 2015, 25(4): 531-535.

    [18] ?QIU H Z, CHEN X Q, WU Q H, et al. Deformation mechanism and collapse treatment of the rock surrounding a shallow tunnel based on on-site monitoring [J]. Journal of Mountain Science, 2020, 17(12): 2897-2914.

    [19] ?LI W J, ZHANG S M, HAN X M. Judgement method for surrounding rock stability of Guanjiao tunnel based on catastrophe theory [J]. Applied Mechanics and Materials, 2011, 90-93: 2307-2312.

    [20] ?楊忠民, 高永濤, 吳順川, 等. 隧道大變形機制及處治關(guān)鍵技術(shù)模型試驗研究[J]. 巖土力學(xué), 2018, 39(12): 4482-4492.

    YANG Z M, GAO Y T, WU S C, et al. Physical model test on large deformation mechanism and key treatment techniques of tunnel [J]. Rock and Soil Mechanics, 2018, 39(12): 4482-4492. (in Chinese)

    [21] ?孫訓(xùn)方, 方孝淑, 關(guān)來泰. 材料力學(xué)[M]. 5版. 北京: 高等教育出版社, 2009.

    SUN X F, FANG S X, GUAN L T. Mechanics of materials [M]. 5th Edition. Beijing: Higher Education Press, 2009. (in Chinese)

    (編輯??王秀玲)

    猜你喜歡
    隧道工程數(shù)值模擬
    隧道工程與水環(huán)境的相互作用
    卷宗(2016年10期)2017-01-21 16:32:59
    某公路隧道二次襯砌施工控制技術(shù)
    隧道工程施工成本管理策略探析
    隧道工程中的軟巖支護技術(shù)
    居業(yè)(2016年9期)2016-12-26 08:39:32
    隧道工程教學(xué)體系改革方法探索
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    亚洲久久久久久中文字幕| 国产真实伦视频高清在线观看| 人人妻人人看人人澡| 尾随美女入室| 日本一本二区三区精品| 建设人人有责人人尽责人人享有的 | 久久久成人免费电影| 国产真实伦视频高清在线观看| 亚洲自拍偷在线| 亚洲av中文av极速乱| 精品酒店卫生间| 欧美区成人在线视频| av在线蜜桃| 肉色欧美久久久久久久蜜桃 | 久久热精品热| 国产真实伦视频高清在线观看| 黄片wwwwww| 九九久久精品国产亚洲av麻豆| 丰满人妻一区二区三区视频av| 天堂√8在线中文| 久久久精品94久久精品| 综合色丁香网| 一个人免费在线观看电影| 亚洲va在线va天堂va国产| 久久久久精品性色| 亚洲精品日本国产第一区| 日韩亚洲欧美综合| 日韩 亚洲 欧美在线| 国产v大片淫在线免费观看| 国产高清三级在线| 人妻一区二区av| 人体艺术视频欧美日本| 亚洲精品自拍成人| 国产精品久久久久久精品电影| 噜噜噜噜噜久久久久久91| 毛片一级片免费看久久久久| 亚洲欧美一区二区三区国产| 亚洲国产成人一精品久久久| 免费少妇av软件| 日本熟妇午夜| 三级国产精品欧美在线观看| 亚洲欧美日韩无卡精品| 国产欧美另类精品又又久久亚洲欧美| 亚洲av电影在线观看一区二区三区 | 免费看不卡的av| 波野结衣二区三区在线| 两个人的视频大全免费| 日韩一本色道免费dvd| 街头女战士在线观看网站| 免费看a级黄色片| 国产伦精品一区二区三区视频9| 亚洲欧洲日产国产| 久久久久久国产a免费观看| 一个人免费在线观看电影| 亚洲精品久久午夜乱码| 性色avwww在线观看| 51国产日韩欧美| 亚洲精品国产av蜜桃| 免费人成在线观看视频色| 99久久九九国产精品国产免费| 中文天堂在线官网| 亚洲av成人精品一二三区| 亚洲精品国产成人久久av| 成人无遮挡网站| 欧美丝袜亚洲另类| videos熟女内射| 欧美+日韩+精品| 国产免费一级a男人的天堂| 1000部很黄的大片| 亚洲精品,欧美精品| 在线观看一区二区三区| 麻豆国产97在线/欧美| 久久综合国产亚洲精品| 最后的刺客免费高清国语| 国产精品美女特级片免费视频播放器| 国产精品久久久久久久久免| 国产亚洲av片在线观看秒播厂 | freevideosex欧美| 亚洲av成人精品一二三区| 男女国产视频网站| 国产成人a∨麻豆精品| 国产黄色视频一区二区在线观看| 免费看日本二区| 国产成人a区在线观看| 欧美区成人在线视频| 高清日韩中文字幕在线| 国产精品美女特级片免费视频播放器| 日韩成人伦理影院| 亚洲av成人精品一二三区| 最后的刺客免费高清国语| 又黄又爽又刺激的免费视频.| 欧美日韩一区二区视频在线观看视频在线 | 亚洲乱码一区二区免费版| 国产精品精品国产色婷婷| 最近的中文字幕免费完整| 免费看不卡的av| 欧美日韩在线观看h| 91在线精品国自产拍蜜月| 精品久久久久久成人av| 国产高清不卡午夜福利| 国产精品日韩av在线免费观看| 国产淫片久久久久久久久| 日韩一区二区视频免费看| 亚洲av国产av综合av卡| 亚洲aⅴ乱码一区二区在线播放| 别揉我奶头 嗯啊视频| 亚洲国产高清在线一区二区三| 午夜激情久久久久久久| 欧美3d第一页| 国产有黄有色有爽视频| 综合色av麻豆| 午夜福利成人在线免费观看| 精品熟女少妇av免费看| av在线天堂中文字幕| 国产一级毛片在线| 亚洲欧美一区二区三区黑人 | 大香蕉久久网| 九色成人免费人妻av| 国产麻豆成人av免费视频| 国产伦一二天堂av在线观看| av免费在线看不卡| 国产乱来视频区| 国产成人freesex在线| 又大又黄又爽视频免费| 精品欧美国产一区二区三| 亚洲av中文字字幕乱码综合| 精品国产三级普通话版| 免费黄色在线免费观看| 国产精品99久久久久久久久| 亚洲av免费在线观看| 能在线免费观看的黄片| 欧美不卡视频在线免费观看| 亚洲国产欧美人成| 亚洲图色成人| 亚洲精品自拍成人| 国产高清有码在线观看视频| 久久精品熟女亚洲av麻豆精品 | 亚洲国产精品成人久久小说| 亚洲欧美成人精品一区二区| 最近最新中文字幕免费大全7| 精品欧美国产一区二区三| 看黄色毛片网站| 91av网一区二区| 国产精品久久视频播放| 一级爰片在线观看| 亚洲精品,欧美精品| 少妇猛男粗大的猛烈进出视频 | 国产午夜福利久久久久久| 2018国产大陆天天弄谢| 狠狠精品人妻久久久久久综合| 熟女电影av网| 激情 狠狠 欧美| 午夜精品国产一区二区电影 | 免费看a级黄色片| 在线观看免费高清a一片| 最近最新中文字幕大全电影3| 在线观看av片永久免费下载| 美女内射精品一级片tv| 国产伦理片在线播放av一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜精品国产一区二区电影 | 亚洲成人精品中文字幕电影| 老司机影院成人| 欧美高清成人免费视频www| av国产免费在线观看| 精品国内亚洲2022精品成人| 美女大奶头视频| 午夜老司机福利剧场| 成人一区二区视频在线观看| 搡女人真爽免费视频火全软件| 亚洲av免费高清在线观看| 亚洲国产日韩欧美精品在线观看| 五月玫瑰六月丁香| 蜜桃久久精品国产亚洲av| 三级国产精品欧美在线观看| 免费观看无遮挡的男女| 午夜久久久久精精品| 天天躁夜夜躁狠狠久久av| 国产av国产精品国产| 国产成人精品一,二区| 婷婷色麻豆天堂久久| 欧美zozozo另类| 黄色欧美视频在线观看| 最新中文字幕久久久久| 午夜激情福利司机影院| 美女脱内裤让男人舔精品视频| 亚洲av成人av| 亚洲美女搞黄在线观看| 亚洲精品视频女| 嫩草影院入口| 日韩成人伦理影院| 在线免费观看的www视频| 99热这里只有是精品在线观看| 精品酒店卫生间| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品国产三级国产专区5o| 男人狂女人下面高潮的视频| 伊人久久精品亚洲午夜| 人人妻人人澡欧美一区二区| 国产一区亚洲一区在线观看| 一区二区三区四区激情视频| 又爽又黄无遮挡网站| 91精品一卡2卡3卡4卡| 久久热精品热| 国产精品一区二区三区四区久久| 在现免费观看毛片| 日本一本二区三区精品| 亚洲国产高清在线一区二区三| 一级毛片 在线播放| 两个人视频免费观看高清| 91精品国产九色| 日日摸夜夜添夜夜添av毛片| freevideosex欧美| 一级片'在线观看视频| 最近的中文字幕免费完整| 精品久久久久久久人妻蜜臀av| 久久久久精品久久久久真实原创| 色视频www国产| 一级毛片黄色毛片免费观看视频| 91久久精品国产一区二区成人| 国产日韩欧美在线精品| 亚洲欧洲国产日韩| 男人狂女人下面高潮的视频| 日韩一本色道免费dvd| 好男人视频免费观看在线| 80岁老熟妇乱子伦牲交| 在线免费观看的www视频| 男人狂女人下面高潮的视频| 中文字幕av成人在线电影| 视频中文字幕在线观看| 免费黄网站久久成人精品| 九九爱精品视频在线观看| 日日干狠狠操夜夜爽| or卡值多少钱| 国产精品久久久久久av不卡| 国产成人aa在线观看| 久久久久久久亚洲中文字幕| 国内精品美女久久久久久| 亚洲最大成人av| 蜜桃亚洲精品一区二区三区| 国产毛片a区久久久久| 在线观看一区二区三区| 欧美日韩视频高清一区二区三区二| 日本一二三区视频观看| 免费观看精品视频网站| freevideosex欧美| 18禁在线播放成人免费| 婷婷色av中文字幕| 日韩,欧美,国产一区二区三区| 国产一区二区三区av在线| 国产亚洲午夜精品一区二区久久 | 舔av片在线| 97精品久久久久久久久久精品| 精品久久久噜噜| 亚洲精品自拍成人| 久久97久久精品| 中文在线观看免费www的网站| 高清午夜精品一区二区三区| 日韩三级伦理在线观看| 亚洲精品日本国产第一区| 赤兔流量卡办理| 成人毛片60女人毛片免费| 噜噜噜噜噜久久久久久91| 日本黄色片子视频| 国产成人freesex在线| 人人妻人人澡欧美一区二区| 一级二级三级毛片免费看| a级毛色黄片| 欧美人与善性xxx| 一级毛片 在线播放| 在线免费十八禁| 夫妻性生交免费视频一级片| 久久99热这里只频精品6学生| 亚洲三级黄色毛片| 色综合亚洲欧美另类图片| 婷婷色av中文字幕| 免费看a级黄色片| 男人爽女人下面视频在线观看| 欧美不卡视频在线免费观看| 一区二区三区四区激情视频| 在线播放无遮挡| 国产av不卡久久| 亚洲av在线观看美女高潮| av播播在线观看一区| 亚洲国产色片| 国产女主播在线喷水免费视频网站 | av专区在线播放| av免费在线看不卡| 天堂影院成人在线观看| 国内揄拍国产精品人妻在线| 欧美成人精品欧美一级黄| 人妻夜夜爽99麻豆av| av在线天堂中文字幕| 国产黄色视频一区二区在线观看| 精品亚洲乱码少妇综合久久| 国产黄片美女视频| 免费不卡的大黄色大毛片视频在线观看 | 人人妻人人看人人澡| 3wmmmm亚洲av在线观看| 性色avwww在线观看| 夜夜爽夜夜爽视频| 日韩欧美精品v在线| 欧美成人一区二区免费高清观看| 成人综合一区亚洲| 日本黄大片高清| 国产乱来视频区| 亚洲经典国产精华液单| 亚洲av成人av| 亚洲国产色片| 99re6热这里在线精品视频| 亚洲精品成人久久久久久| 性色avwww在线观看| 成人综合一区亚洲| 成人毛片60女人毛片免费| 一个人免费在线观看电影| 亚洲精品国产av蜜桃| 国产成人精品久久久久久| 国产一区亚洲一区在线观看| 超碰97精品在线观看| 一个人看视频在线观看www免费| 久久这里只有精品中国| 最近手机中文字幕大全| 免费大片黄手机在线观看| 有码 亚洲区| 内地一区二区视频在线| 国产永久视频网站| 精品人妻熟女av久视频| 1000部很黄的大片| 能在线免费观看的黄片| 国产女主播在线喷水免费视频网站 | 搞女人的毛片| 最近中文字幕高清免费大全6| 色5月婷婷丁香| 日韩强制内射视频| 国产精品人妻久久久久久| 国产亚洲午夜精品一区二区久久 | 91久久精品电影网| 简卡轻食公司| 国产永久视频网站| 少妇人妻精品综合一区二区| freevideosex欧美| 国产69精品久久久久777片| 亚洲精品aⅴ在线观看| 一二三四中文在线观看免费高清| 精品国内亚洲2022精品成人| 少妇人妻一区二区三区视频| 男的添女的下面高潮视频| 久久热精品热| 免费观看无遮挡的男女| 亚洲人成网站在线观看播放| 国产精品av视频在线免费观看| 水蜜桃什么品种好| 久久久久国产网址| 夫妻午夜视频| 久久草成人影院| 永久免费av网站大全| 亚洲av电影在线观看一区二区三区 | 成人漫画全彩无遮挡| 亚洲av成人精品一区久久| eeuss影院久久| 人妻少妇偷人精品九色| 全区人妻精品视频| 欧美日韩视频高清一区二区三区二| 国产欧美日韩精品一区二区| 草草在线视频免费看| 成人av在线播放网站| 欧美人与善性xxx| 精品午夜福利在线看| 亚洲精品456在线播放app| av.在线天堂| 久久久久久久午夜电影| 女的被弄到高潮叫床怎么办| 国产午夜福利久久久久久| 啦啦啦中文免费视频观看日本| 久久久久久久亚洲中文字幕| 中文欧美无线码| 能在线免费看毛片的网站| 午夜免费观看性视频| 久久99热这里只有精品18| 男人和女人高潮做爰伦理| 国产成人一区二区在线| 国产欧美另类精品又又久久亚洲欧美| 一级av片app| 日本与韩国留学比较| 国产亚洲一区二区精品| 3wmmmm亚洲av在线观看| 国产在线一区二区三区精| www.色视频.com| 中文字幕人妻熟人妻熟丝袜美| 美女黄网站色视频| 欧美潮喷喷水| 99九九线精品视频在线观看视频| 看免费成人av毛片| 性插视频无遮挡在线免费观看| 国产精品一区二区在线观看99 | 少妇的逼水好多| 淫秽高清视频在线观看| 伦精品一区二区三区| 搡老妇女老女人老熟妇| 国产精品1区2区在线观看.| 午夜老司机福利剧场| 精品熟女少妇av免费看| 国产黄色视频一区二区在线观看| av播播在线观看一区| 久热久热在线精品观看| www.色视频.com| 又黄又爽又刺激的免费视频.| 亚洲av福利一区| 成人一区二区视频在线观看| 青春草亚洲视频在线观看| 日韩成人伦理影院| 女人十人毛片免费观看3o分钟| 一个人看的www免费观看视频| 日韩欧美精品免费久久| 日本黄大片高清| 国产成人aa在线观看| 免费高清在线观看视频在线观看| 中文字幕制服av| 热99在线观看视频| 久久久久久久亚洲中文字幕| 久久久精品94久久精品| 伊人久久精品亚洲午夜| 丰满少妇做爰视频| av网站免费在线观看视频 | 国产一级毛片七仙女欲春2| 在线免费十八禁| 2022亚洲国产成人精品| 18禁动态无遮挡网站| 日韩精品有码人妻一区| 亚洲在线自拍视频| 激情五月婷婷亚洲| 日韩,欧美,国产一区二区三区| 亚洲精品影视一区二区三区av| 大陆偷拍与自拍| 亚洲欧美清纯卡通| av网站免费在线观看视频 | ponron亚洲| a级毛色黄片| eeuss影院久久| 国产在视频线在精品| 又粗又硬又长又爽又黄的视频| 久久这里有精品视频免费| 99九九线精品视频在线观看视频| 三级经典国产精品| 午夜精品在线福利| 婷婷色麻豆天堂久久| 国产v大片淫在线免费观看| 日本-黄色视频高清免费观看| 国产在视频线精品| 亚洲精品视频女| 久久草成人影院| www.色视频.com| 免费观看精品视频网站| 一级毛片久久久久久久久女| av天堂中文字幕网| 在线免费十八禁| 亚洲乱码一区二区免费版| 99热这里只有精品一区| 爱豆传媒免费全集在线观看| 97超碰精品成人国产| 女人十人毛片免费观看3o分钟| 最近中文字幕2019免费版| 日本欧美国产在线视频| 婷婷色麻豆天堂久久| 欧美 日韩 精品 国产| 国产黄频视频在线观看| 麻豆乱淫一区二区| 国产午夜精品久久久久久一区二区三区| 18禁动态无遮挡网站| 精品人妻视频免费看| 日本一本二区三区精品| 国产一区二区亚洲精品在线观看| 精品少妇黑人巨大在线播放| 一级毛片我不卡| 天天躁夜夜躁狠狠久久av| 国产精品伦人一区二区| 丰满人妻一区二区三区视频av| 波多野结衣巨乳人妻| 中国国产av一级| 99久久九九国产精品国产免费| 大香蕉97超碰在线| 久久精品人妻少妇| 精品99又大又爽又粗少妇毛片| 国产一区亚洲一区在线观看| 国产午夜精品论理片| 好男人视频免费观看在线| 边亲边吃奶的免费视频| 如何舔出高潮| 国产精品一区二区三区四区久久| 麻豆成人av视频| 欧美成人a在线观看| 日本-黄色视频高清免费观看| 国产精品久久久久久精品电影小说 | 中文字幕人妻熟人妻熟丝袜美| 人妻夜夜爽99麻豆av| 春色校园在线视频观看| 国产男女超爽视频在线观看| 69av精品久久久久久| 久久久久国产网址| 中国国产av一级| 色5月婷婷丁香| 天美传媒精品一区二区| 成人高潮视频无遮挡免费网站| 欧美激情在线99| 黄色日韩在线| 能在线免费看毛片的网站| 成人亚洲欧美一区二区av| 亚洲精品影视一区二区三区av| 成年版毛片免费区| 亚洲aⅴ乱码一区二区在线播放| 国产精品99久久久久久久久| 狂野欧美激情性xxxx在线观看| 美女被艹到高潮喷水动态| 国产成人福利小说| 成人欧美大片| 久久久国产一区二区| 黄色欧美视频在线观看| 国产精品一区二区三区四区久久| 一级爰片在线观看| 日韩三级伦理在线观看| 久久久久久久亚洲中文字幕| 爱豆传媒免费全集在线观看| 99热这里只有是精品50| xxx大片免费视频| 亚洲在线自拍视频| 亚洲18禁久久av| 欧美97在线视频| 亚洲欧洲日产国产| 免费av观看视频| 综合色av麻豆| 久久精品国产亚洲av涩爱| 一个人免费在线观看电影| 久久久久国产网址| 欧美zozozo另类| 免费无遮挡裸体视频| 日本av手机在线免费观看| 国产伦在线观看视频一区| 久久综合国产亚洲精品| 不卡视频在线观看欧美| 69av精品久久久久久| 国产老妇女一区| 亚洲四区av| 在线观看av片永久免费下载| 搡老乐熟女国产| 精品一区二区三区视频在线| 午夜日本视频在线| 秋霞在线观看毛片| 免费高清在线观看视频在线观看| 免费看日本二区| 精品久久久久久电影网| 亚洲三级黄色毛片| 日韩强制内射视频| 色5月婷婷丁香| av在线播放精品| 国产男女超爽视频在线观看| 午夜福利视频1000在线观看| 国产精品一区二区性色av| 最近2019中文字幕mv第一页| 免费看a级黄色片| 如何舔出高潮| 精品一区二区免费观看| 日本一本二区三区精品| 精品久久久久久久久av| 久久久久免费精品人妻一区二区| 一级二级三级毛片免费看| 国产探花极品一区二区| 欧美成人a在线观看| 你懂的网址亚洲精品在线观看| 久久久久网色| 久久久午夜欧美精品| 特大巨黑吊av在线直播| 草草在线视频免费看| 十八禁国产超污无遮挡网站| 色哟哟·www| 欧美激情国产日韩精品一区| 国产黄片视频在线免费观看| 成人高潮视频无遮挡免费网站| 色播亚洲综合网| 看十八女毛片水多多多| 2021少妇久久久久久久久久久| 日本免费在线观看一区| 韩国av在线不卡| www.色视频.com| a级一级毛片免费在线观看| 国国产精品蜜臀av免费| 91精品一卡2卡3卡4卡| 久久精品久久精品一区二区三区| 亚洲成人av在线免费| 51国产日韩欧美| 日韩不卡一区二区三区视频在线| 国产探花在线观看一区二区| 一个人免费在线观看电影| 一级毛片黄色毛片免费观看视频| 2021少妇久久久久久久久久久| 大陆偷拍与自拍| 亚洲欧美中文字幕日韩二区| 日产精品乱码卡一卡2卡三| 日日啪夜夜爽| 国产精品爽爽va在线观看网站| 不卡视频在线观看欧美| 亚洲va在线va天堂va国产| 成人毛片a级毛片在线播放| 麻豆精品久久久久久蜜桃| 亚洲精品久久久久久婷婷小说| 麻豆av噜噜一区二区三区| eeuss影院久久| 日本一本二区三区精品| 淫秽高清视频在线观看| 国产成人精品福利久久| 99久国产av精品国产电影| 男女边摸边吃奶| 女的被弄到高潮叫床怎么办| 日韩三级伦理在线观看| 亚洲精品一二三| 六月丁香七月|