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

    拱壩聯(lián)合壩肩巖體三維計(jì)算模型構(gòu)建

    2014-05-18 08:04:20羅先啟鄭安興畢金鋒翁永紅
    巖土力學(xué) 2014年5期
    關(guān)鍵詞:楔塊拱壩壩址

    沈 輝,羅先啟,鄭安興,畢金鋒,翁永紅

    (1.上海交通大學(xué) 船舶海洋與建筑工程學(xué)院,上海 200240;2.長江水利委員會 長江勘測規(guī)劃設(shè)計(jì)研究院,武漢 430010)

    1 引 言

    拱壩壩肩巖體穩(wěn)定性分析是一個(gè)受多種因素影響的地質(zhì)力學(xué)問題。通常需要在地質(zhì)勘探、現(xiàn)場調(diào)查研究的基礎(chǔ)上,通過地質(zhì)條件綜合分析建立地質(zhì)體的概念模型,通過地質(zhì)力學(xué)模型試驗(yàn)或數(shù)值模擬方法模擬地質(zhì)體力學(xué)環(huán)境演化的全過程,以達(dá)到預(yù)測分析的目的。

    隨著理論分析的成熟和計(jì)算機(jī)技術(shù)的發(fā)展,在大型水利水電工程及巖土工程等領(lǐng)域,數(shù)值模擬方法以其高效低成本的特點(diǎn),成為地質(zhì)力學(xué)問題分析的關(guān)鍵技術(shù)手段。從目前的研究情況來看,地勘資料的詳細(xì)程度及研究問題的不同,決定了數(shù)值模型的精細(xì)度及規(guī)模,也使得分析人員根據(jù)實(shí)際情況選擇合理的建模及分析軟件[1]。由于地質(zhì)模型建立、空間單元?jiǎng)澐值惹疤幚砉ぷ鳠┈崳坏貌徊捎煤喕刭|(zhì)模型的做法,對工程建筑物近域地基中復(fù)雜的地質(zhì)條件以及壩址河谷兩岸實(shí)際地形的影響考慮較少。文獻(xiàn)[2]采用ANSYS建立拱壩及概化密集節(jié)理模型,探討密集節(jié)理對壩肩穩(wěn)定的影響問題;文獻(xiàn)[3]以大型通用有限元程序ADINA為計(jì)算平臺,對比分析壩肩邊坡卸荷后,不同地震加速度及不同錨固長度和錨固噸位在卸荷條件下關(guān)鍵點(diǎn)的位移和應(yīng)力以及塑性區(qū)的變化規(guī)律;文獻(xiàn)[4]采用考慮開挖過程的有限元強(qiáng)度折減法分析錦屏一級水電站左岸壩肩開挖邊坡在施工期的穩(wěn)定性。也有許多學(xué)者[1,5]將三維地學(xué)建模與數(shù)值模擬進(jìn)行耦合,讓數(shù)值模擬直接繼承和利用三維地質(zhì)詳細(xì)的地勘數(shù)據(jù)。但需要掌握三維地學(xué)建模系統(tǒng)的格式和算法,編制計(jì)算程序復(fù)雜,不利于推廣。

    本文以烏東德超高拱壩壩肩穩(wěn)定性分析數(shù)值模型為例,對可研階段地質(zhì)勘探數(shù)據(jù)、二維地質(zhì)剖面和工程設(shè)計(jì)數(shù)據(jù)進(jìn)行預(yù)處理,提取包括地形地貌、區(qū)域地層、斷層構(gòu)造、巖溶系統(tǒng)等自然地質(zhì)對象界限,以及開挖面、拱壩等人工對象邊界?;贖ypermesh軟件對各類對象進(jìn)行插值、擬合和幾何建模,利用其強(qiáng)大的網(wǎng)格參數(shù)交互調(diào)整功能以及獨(dú)立單元與節(jié)點(diǎn)的處理功能,建立拱壩聯(lián)合壩址區(qū)巖體的三維數(shù)值精細(xì)模型,為壩址區(qū)地應(yīng)力反演分析,壩體載荷、工程開挖等因素影響下壩肩穩(wěn)定性分析提供了良好的分析平臺。

    2 工程地質(zhì)概況

    烏東德水電站大壩為混凝土雙曲拱壩,壩頂高程為988 m,壩底高程為723 m,最大壩高為265 m。壩址處河谷形態(tài)呈狹窄的“V”型,兩岸地形對稱,高陡險(xiǎn)峻,坡形完整。拱座巖體為落雪組第3段第1~4亞段(P t32-11~P t32-14)變質(zhì)灰?guī)r、大理巖,以厚層為主,夾中厚層、巨厚層及少量薄層;地層產(chǎn)狀走向?yàn)?60°~280°,傾向SW,傾角為70°~85°,巖層走向與河流流向交角為 60°~80°,為斜橫向谷。拱座巖體中沒有斜切邊坡、貫穿性強(qiáng)的斷層和裂隙性斷層發(fā)育。

    根據(jù)可研階段資料,大壩拱肩槽呈弧形開挖,左岸開挖邊坡走向從62°轉(zhuǎn)至100°,開挖走向與巖層走向接近平行,拱肩槽正面坡為橫向坡;右岸開挖邊坡走向從68°轉(zhuǎn)至 76°,拱肩槽正面坡為斜交順向坡。左岸開挖邊坡坡頂高程約為1163 m,最大坡高約為440 m;右岸開挖邊坡坡頂高程約為1120 m,最大坡高約為397 m。

    左壩肩巖體主要為 Pt2l3-1厚層灰?guī)r夾大理巖,局部有少量(高程為810 m以下)Pt2l2-3互層大理巖化白云巖。右壩肩中下部(高程為875~905 m以下)為Pt2l3-1厚層灰?guī)r夾大理巖;上部(高程875~905 m以上)為Pt2l3-2-1~Pt2l2-2-6巨厚層白云巖夾中厚層灰?guī)r、厚層石英巖。兩岸壩肩建基巖體以Ⅱ1級為主,巖體質(zhì)量優(yōu)良,左壩肩主要為Ⅱ1(95.2%),少量為Ⅲ1(4.8%),右壩肩主要為Ⅱ1(82.7%),少量為Ⅲ1(17.3%)。

    3 拱壩聯(lián)合壩址區(qū)巖體三維模型

    3.1 拱壩壩址區(qū)自然地質(zhì)條件三維模型

    拱壩壩址區(qū)自然地質(zhì)條件的三維建模主要包括計(jì)算范圍的選擇、地形地貌、關(guān)鍵地質(zhì)結(jié)構(gòu)(地層、斷層、溶洞系統(tǒng)、軟弱夾層等)。

    三維計(jì)算模型以壩體為中心,沿順河向、橫河向及壩底深部基巖以下一定的范圍內(nèi)生成。下游模擬范圍及左右兩岸模擬范圍均大于2.5倍壩高。整體模擬范圍為1225 m×1225 m×1440 m。

    壩址區(qū)會理群(Pt2H)地層主要出露因民組(Pt2y)和落雪組(Pt2l)兩組地層,巖性有8類。根據(jù)地質(zhì)勘探數(shù)據(jù)及二維地質(zhì)剖面信息,將設(shè)計(jì)單位提供的 CAD格式的工程地質(zhì)圖、地形圖、各高程平切圖等進(jìn)行預(yù)處理,剔除非幾何邊界的冗余信息,并直接輸出為 DXF矢量數(shù)據(jù)格式,包含地表及巖性分界線的DXF數(shù)據(jù)文件直接輸入Hypermesh軟件,對地表、地層及斷層進(jìn)行插值和擬合處理,壩址區(qū)巖層幾何模型如圖1所示。

    壩址區(qū)地表實(shí)測斷層有35條,規(guī)模較大。對壩肩穩(wěn)定性有影響的斷層有4條,從上游至下游分別為紅溝斷層F3、梁子斷層F7、白溝斷層F14、雷家灣溝斷層F15。F15走向NNE,與河谷斜交。其余斷層與河谷大角度相交,壩線選擇段地表斷層不發(fā)育,可見9條小斷層,模型考慮距離壩體較近的次一級斷層 f40和f42。右岸存在已停止發(fā)育的充填古巖溶系統(tǒng)(K25),巖溶系統(tǒng)發(fā)育于Pt2l3-3厚層灰?guī)r與Pt2l3-2厚層白云巖分界面附近 Pt2l3-3厚層灰?guī)r中,基本上順層發(fā)育,高程向下具向坡外發(fā)育和規(guī)模逐漸減小的趨勢,平切面上呈橢圓狀,長軸與巖層走向方面近于一致,短軸與巖層傾向近于一致。斷層和巖溶與壩體相對位置如圖3所示。

    圖1 巖層三維地質(zhì)模型Fig.1 3D geological model for rock mass

    圖2 斷層和巖溶與壩體的位置Fig.2 Positions of dam,faults and karst

    圖3 整體有限元網(wǎng)格模型Fig.3 Integral finite element model

    借助Hypermesh軟件在幾何體的布爾(Boolean)操作運(yùn)算方面的優(yōu)越性,模型三維實(shí)體模型采用圖2結(jié)構(gòu)面切割圖1巖層幾何體完成。異形幾何體的空間復(fù)雜關(guān)系使得有限元網(wǎng)格的劃分工作十分繁瑣。利用HyperMesh的網(wǎng)格劃分模塊交互調(diào)整特性,方便調(diào)整曲面或邊界的網(wǎng)格參數(shù),特別是對單個(gè)單元及節(jié)點(diǎn)的調(diào)整功能,提高不規(guī)則幾何體網(wǎng)格劃分效率。

    本文巖體全部采用8節(jié)點(diǎn)六面體和6節(jié)點(diǎn)五面體實(shí)體單元,F(xiàn)3、F14、F15及 K25采用實(shí)體單元模擬。小斷層f40和f42采用無厚度節(jié)理單元模擬,計(jì)算網(wǎng)格如圖 3所示。有限元網(wǎng)格劃分后單元總數(shù)為294992,節(jié)點(diǎn)總數(shù)為 300272(壩體單元總數(shù)為10040)。根據(jù)巖體力學(xué)參數(shù)建議值,地層巖性計(jì)算力學(xué)參數(shù)見表1。

    表1 地層巖性計(jì)算力學(xué)參數(shù)Table 1 Mechanical parameters of materials for calculation

    3.2 壩址區(qū)初始地應(yīng)力反演分析三維模型

    假定現(xiàn)今初始地應(yīng)力場主要是遠(yuǎn)古地應(yīng)力場在后期的地表剝蝕與河流侵蝕下切等外動(dòng)力地質(zhì)作用下逐漸形成。根據(jù)壩址區(qū)實(shí)測地應(yīng)力資料,假定應(yīng)力邊界條件,反演最優(yōu)地表剝蝕和河流侵蝕分層厚度[6-7]。計(jì)算需模擬從遠(yuǎn)古模擬夷平面地表剝蝕至現(xiàn)今壩址區(qū)河谷地貌,反演構(gòu)造應(yīng)力邊界條件。將經(jīng)過神經(jīng)網(wǎng)絡(luò)反演得到的邊界條件和重力加速度代入進(jìn)行正算,從而獲取壩址區(qū)初始地應(yīng)力場。

    本文壩址區(qū)夷平面的高程選取1800 m高程作為地表剝蝕和河流侵蝕下切的起點(diǎn),以壩址區(qū)現(xiàn)今河谷高程最低的點(diǎn)(高程值為723 m)為地表侵蝕下切的終點(diǎn),河谷的下切過程分為6次完成。現(xiàn)今地表與遠(yuǎn)古夷平面間平均插值5個(gè)侵蝕平面,利用Hypermesh軟件網(wǎng)格投影功能,地表面網(wǎng)格分別在5個(gè)侵蝕面及遠(yuǎn)古夷平面進(jìn)行二維投影。侵蝕面間實(shí)體網(wǎng)格形態(tài)由侵蝕面二維預(yù)置網(wǎng)格控制,由此方法可快速構(gòu)建地應(yīng)力反演分析數(shù)值模型,如圖4所示。

    圖4 初始地應(yīng)力反演分析網(wǎng)格模型Fig.4 Finite element model for initial geostress field inversion

    3.3 壩肩槽開挖分析三維模型

    模型高邊坡包括壩肩壩頂高程以上邊坡及拱肩槽邊坡。左岸高邊坡開口線最高高程為1155 m,平均開挖坡度為75°;右岸高邊坡開口線最高高程為1135 m,平均開挖坡度為65°。拱肩槽邊坡平面形態(tài)似弧形,略凸向上游側(cè)。根據(jù)上、下游邊坡不同形態(tài),可進(jìn)一步分為拱肩槽上游邊坡、正面邊坡及拱肩槽下游邊坡。左右岸開挖邊坡大致呈對稱分布,兩岸一般每隔20 m設(shè)置一馬道。開挖后地形地貌基本順原坡體走向,退坡深度不大。

    根據(jù)壩肩邊坡設(shè)計(jì)及壩體體型方案,在 CAD軟件中按照開挖線高程及壩體橫剖面高程,構(gòu)建人工對象邊界模型。導(dǎo)入到Hypermesh軟件,快速生成曲面,與壩址區(qū)三維地質(zhì)模型進(jìn)行體與面之間的布爾運(yùn)算,將開挖巖體部分切割分離。

    圖5 拱肩開挖邊坡網(wǎng)格模型Fig.5 Finite element model for abutment slope excavation

    3.4 壩肩巖體滑動(dòng)楔塊三維滑面模型

    地勘裂隙統(tǒng)計(jì)資料表明,拱座內(nèi)無定位塊體,控制拱座滑動(dòng)穩(wěn)定的主要是短小裂隙形成的非連通的滑動(dòng)邊界。裂隙可分為SN組和EW組,其中SN組與河流平行,兩岸均發(fā)育,占總數(shù)的54.1%,EW組占總數(shù)的 15.4%,與河流近垂直,兩岸均發(fā)育。裂隙以中傾角為主,次為陡傾角,緩傾角裂隙不發(fā)育,中傾角裂隙以SN組裂隙為主,優(yōu)勢產(chǎn)狀為傾向273°、傾角 51°[8]。根據(jù)裂隙產(chǎn)狀和規(guī)模及巖層分布情況,確定9個(gè)可能滑動(dòng)楔塊,典型楔塊滑面如圖6所示。

    圖6 典型楔塊滑面Fig.6 Sliding surface of typical wedge block

    本文研究的問題中存在多個(gè)可能滑動(dòng)楔塊的情況下,若三維建模時(shí)考慮所有滑面情況,則使得拱座內(nèi)巖體支離破碎,單元網(wǎng)格質(zhì)量差,操作繁瑣,不利于推廣。

    本文采用的是采用滑面多重網(wǎng)格的方法,即將有限元計(jì)算的應(yīng)力結(jié)果插值到重建的滑面網(wǎng)格上,從而適應(yīng)復(fù)雜滑面的組合[9]。在Hypermesh軟件中確定楔塊滑面后,利用卸荷帶界線切割楔塊滑面并進(jìn)行網(wǎng)格劃分。拱座內(nèi)滑動(dòng)楔塊由新鮮巖體、弱卸荷巖體及強(qiáng)卸荷巖體組成,如圖7所示。

    圖7 楔塊巖體卸荷帶界線Fig.7 Classification of unloading zones of wedge blocks

    4 計(jì)算結(jié)果

    將三維數(shù)值模型通過接口程序進(jìn)行單元數(shù)據(jù)格式轉(zhuǎn)換[10],模型導(dǎo)入到FLAC3D程序中進(jìn)行三維彈塑性數(shù)值計(jì)算。根據(jù)地應(yīng)力反演結(jié)果,壩址區(qū)地應(yīng)力分布可分為3個(gè)區(qū)域,沿壩軸線剖面初始地應(yīng)力s1分布如圖8所示。

    圖8 壩軸線剖面地應(yīng)力分布Fig.8 Stress distribution at dam axis section

    兩側(cè)河谷近地表地應(yīng)力帶:受地表剝蝕和金沙江長期沖蝕下切作用,最大水平主應(yīng)力量值一般為10 MPa左右,以小于10 MPa的低地應(yīng)力為主;總體上,左岸較右岸稍大,岸坡地形對地應(yīng)力有一定影響。河床深部地應(yīng)力區(qū):兩岸深部高程在 750~850 m范圍內(nèi),應(yīng)力量值屬中等地應(yīng)力水平;河床深部最大主應(yīng)力達(dá)20 MPa左右,局部達(dá)到高應(yīng)力水平。構(gòu)造應(yīng)力區(qū):在斷層區(qū)域,受斷層擠壓剪切影響,斷層兩側(cè)地應(yīng)力突變,以左岸雷家溝斷層F15最為明顯。因此,斷層穿越地段受斷層控制,為應(yīng)力集中的構(gòu)造應(yīng)力帶。

    計(jì)算模擬過程按照初始地應(yīng)力反演、高邊坡及壩肩槽開挖、壩體施工、壩體載荷的步驟進(jìn)行??紤]荷載不同組合,計(jì)算模擬分為4個(gè)方案。方案1荷載組合為正常蓄水+自重+泥沙壓力+浪壓力+冰壓力。方案2荷載組合為方案1荷載組合疊加溫升。方案3荷載組合為方案1荷載組合疊加溫降。方案4荷載組合為方案2荷載組合疊加地震荷載。壩肩槽開挖完成后,壩肩附近巖體順河向位移如圖9所示,右岸拱肩槽 800 m高程附近位移值最大達(dá)到12.9 mm。左右岸可能滑動(dòng)楔塊安全系數(shù)計(jì)算如表2所示。

    正常蓄水后(方案 1),各楔塊具有較高的安全裕度,整體上左岸穩(wěn)定性強(qiáng)于右岸,安全系數(shù)計(jì)算結(jié)果如表2所示。右岸滑面通過巖溶區(qū)域的楔塊相對穩(wěn)定性較差,楔塊 RM1安全系數(shù)最小值為5.399?;窘M合2條件下,溫降時(shí)拱圈縮短并偏向右岸下游方向變形,左岸楔塊相對方案 1,楔塊穩(wěn)定性有所增加。右岸部分楔塊穩(wěn)定性降低,最大降低幅度7.4%。地震荷載作用下所有楔塊安全系數(shù)均大幅下降,左岸楔塊LM2-1安全系數(shù)最小為5.511。

    圖9 拱肩開挖邊坡位移分布(沿河向)Fig.9 Displacement distribution of excavation slope

    表2 楔塊安全系數(shù)Table 2 Safety factors of blocks

    5 結(jié) 論

    采用通用軟件平臺,建立用于高拱壩壩肩穩(wěn)定分析所需的三維數(shù)值模型的方法,并以烏東德工程為例,構(gòu)建適用于壩肩穩(wěn)定性分析的三維精細(xì)數(shù)值模型。利用地勘資料和工程設(shè)計(jì)數(shù)據(jù),分析了影響壩肩穩(wěn)定的地質(zhì)特征和工程措施,采用面向?qū)ο蠹夹g(shù)、多種 CAD軟件協(xié)同作業(yè),建立自然地質(zhì)條件和工程對象的幾何模型。以壩肩部位為中心,基于CAE領(lǐng)域功能強(qiáng)大的Hypermesh軟件,較好地完成了復(fù)雜地質(zhì)結(jié)構(gòu)與工程對象的網(wǎng)格劃分。通過接口程序進(jìn)行三維數(shù)值模型單元格式轉(zhuǎn)換,在 FLAC3D完成三維彈塑性計(jì)算。

    本文將地應(yīng)力反演分析、壩體載荷、工程開挖等多因素影響下壩肩靜動(dòng)力穩(wěn)定性分析的問題統(tǒng)一到獨(dú)立模型中。計(jì)算結(jié)果表明,精細(xì)的三維模型最大程度地模擬了工程巖體復(fù)雜結(jié)構(gòu),為高拱壩壩肩巖體穩(wěn)定性分析提供可靠的基礎(chǔ)。

    [1]李明超,鐘登華,秦朝霞,等.基于三維地質(zhì)模型的工程巖體結(jié)構(gòu)精細(xì)數(shù)值建模[J].巖石力學(xué)與工程學(xué)報(bào),2007,26(9): 1893-1898.LI Ming-chao,ZHONG Deng-hua,QIN Zhao-xia,et al.Refined modeling for numerical simulation of engineering rock mass structures based on 3D geological model[J].Chinese Journal of Rock Mechanics and Engineering,2007,26(9): 1893-1898.

    [2]楊寶全,張林,陳建葉,等.復(fù)雜地質(zhì)條件下拱壩壩肩穩(wěn)定分析及密集節(jié)理影響研究[J].巖石力學(xué)與工程學(xué)報(bào),2010,29(增刊2): 3972-3978.YANG Bao-quan,ZHANG Lin,CHEN Jian-ye,et al.Analysis of arch dam’s abutment stability under complicated geological conditions and study of influence of dense joints[J].Chinese Journal of Rock Mechanics and Engineering,2010,29(Supp.2): 3972-3978.

    [3]劉杰,李建林,宛良朋,等.基于卸荷和抗震分析理論的大崗山壩肩邊坡錨固優(yōu)化研究[J].巖土力學(xué),2012,33(增刊 2): 275-282.LIU Jie,LI Jian-lin,WAN Liang-peng,et al.Research on Dagangshan dam abutment slope anchorage optimization based on theory of unloading and seismic analysis[J].Rock and Soil Mechanics,2012,33(Supp.2): 275-282.

    [4]漆祖芳,姜清輝,唐志丹,等.錦屏一級水電站左岸壩肩邊坡施工期穩(wěn)定分析[J].巖土力學(xué),2012,33(2): 531-538.QI Zu-fang,JIANG Qing-hui,TANG Zhi-dan,et al.Stability analysis of abutment slope at left bank of Jinping-I hydropower project during construction[J].Rock and Soil Mechanics,2012,33(2): 531-538.

    [5]侯恩科,吳立新,李建民,等.三維地學(xué)模擬與數(shù)值模擬的耦合方法研究[J].煤炭學(xué)報(bào),2002,27(4): 388-392.HOU En-ke,WU Li-xin,LI Jian-min,et al.Study on the coupling of 3D geoscience modeling with numerical simulation[J].Journal of China Coal Society,2002,27(4): 388-392.

    [6]張傳慶,周輝,馮夏庭,等.局域地應(yīng)力場獲取的插值平衡方法[J].巖土力學(xué),2008,29(8): 2016-2024.ZHANG Chuan-qing,ZHOU Hui,FENG Xia-ting,et al.Method for obtaining local stress field based on interpolation and numerical simulation method[J].Rock and Soil Mechanics,2008,29(8): 2016-2024.

    [7]張振華.深切河谷岸坡開挖過程動(dòng)態(tài)預(yù)警方法研究[博士學(xué)位論文 D].武漢: 中國科學(xué)院武漢巖土力學(xué)研究所,2008.

    [8]沈輝,羅先啟,李野,等.烏東德拱壩壩肩三維抗滑穩(wěn)定分析[J].巖石力學(xué)與工程學(xué)報(bào),2012,31(5): 1026-1033.SHEN Hui,LUO Xian-qi,LI Ye,et al.Three-dimensional stability analysis of dam abutment of Wudongde arch dam[J].Chinese Journal of Rock Mechanics and Engineering,2012,31(5): 1026-1033.

    [9]楊強(qiáng),朱玲,翟明杰.基于三維非線性有限元的壩肩穩(wěn)定剛體極限平衡法機(jī)理研究[J].巖石力學(xué)與工程學(xué)報(bào),2005,24(19): 3403-3409.YANG Qiang,ZHU Ling,ZHAI Ming-jie.On mechanism of rigid-body limit equilibrium method for abutment stability analysis of arch dam[J].Chinese Journal of Rock Mechanics and Engineering,2005,24(19): 3403-3409.

    [10]劉心庭,唐輝明.FLAC3D復(fù)雜網(wǎng)格模型的構(gòu)建及其工程應(yīng)用[J].金屬礦山,2010,(11): 108-111.LIU Xin-ting,TANG Hui-ming.FLAC3Dcomplex grid model construction and its engineering application[J].Metal Mine,2010,(11): 108-111.

    猜你喜歡
    楔塊拱壩壩址
    庫爾楚水庫壩址比選分析
    漸開線型面楔塊超越離合器接觸特性仿真分析
    Phytochemicals targeting NF-κB signaling:Potential anti-cancer interventions
    拉伸機(jī)夾頭結(jié)構(gòu)形式及楔塊角度的選取
    閘門楔塊的標(biāo)準(zhǔn)化與夾具設(shè)計(jì)
    楔塊重心位置對斜撐離合器性能的影響
    淺析步亭水庫壩址、壩線的方案比選
    淺議高拱壩壩踵實(shí)測與計(jì)算應(yīng)力差異原因
    酒泉市洪水河水庫壩址方案比選
    砌石雙曲拱壩拱冠梁設(shè)計(jì)的探討和實(shí)踐
    在线观看美女被高潮喷水网站| 免费观看无遮挡的男女| 观看免费一级毛片| 亚洲精品亚洲一区二区| av在线老鸭窝| 日韩电影二区| 男女无遮挡免费网站观看| 内地一区二区视频在线| 9色porny在线观看| 久久精品国产亚洲av涩爱| 乱人伦中国视频| 国产无遮挡羞羞视频在线观看| 国产黄色视频一区二区在线观看| 99热全是精品| 亚洲四区av| 一级毛片电影观看| 色哟哟·www| av一本久久久久| 国产精品女同一区二区软件| 精品久久国产蜜桃| 91精品一卡2卡3卡4卡| 在线天堂最新版资源| 国产免费一级a男人的天堂| 成年人免费黄色播放视频 | 国产精品一区二区三区四区免费观看| 国模一区二区三区四区视频| 国产成人免费观看mmmm| 国产成人精品久久久久久| 女人久久www免费人成看片| 亚洲精品一二三| 亚洲人成网站在线观看播放| 国产又色又爽无遮挡免| 99re6热这里在线精品视频| 夫妻性生交免费视频一级片| 欧美精品人与动牲交sv欧美| 哪个播放器可以免费观看大片| 最近最新中文字幕免费大全7| 少妇的逼好多水| 久久99蜜桃精品久久| 在线观看免费日韩欧美大片 | 国产黄频视频在线观看| 亚洲伊人久久精品综合| 国产一区二区在线观看日韩| 亚洲欧美成人综合另类久久久| 国产精品一区二区在线不卡| 一区二区三区乱码不卡18| 黄色怎么调成土黄色| 久久精品夜色国产| 国产成人精品一,二区| 夜夜爽夜夜爽视频| 国产日韩欧美视频二区| 国产亚洲5aaaaa淫片| 亚洲av国产av综合av卡| 久久久久久久久久久久大奶| 亚洲精品国产av蜜桃| 久久久a久久爽久久v久久| 又爽又黄a免费视频| 国产精品99久久99久久久不卡 | 久久毛片免费看一区二区三区| 成人综合一区亚洲| 亚洲三级黄色毛片| 免费看不卡的av| 又爽又黄a免费视频| 日韩免费高清中文字幕av| 大片电影免费在线观看免费| av福利片在线观看| 人人妻人人爽人人添夜夜欢视频 | 99热6这里只有精品| 春色校园在线视频观看| 天堂8中文在线网| 卡戴珊不雅视频在线播放| 国产成人免费无遮挡视频| 欧美日韩视频精品一区| 亚洲精品自拍成人| 熟女电影av网| 人体艺术视频欧美日本| 在线观看免费日韩欧美大片 | 国产精品不卡视频一区二区| 午夜影院在线不卡| 国产伦理片在线播放av一区| 国产淫片久久久久久久久| 激情五月婷婷亚洲| 国产精品人妻久久久影院| 欧美精品国产亚洲| a级毛片免费高清观看在线播放| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品一区三区| 午夜福利影视在线免费观看| 国产乱来视频区| 丁香六月天网| 午夜老司机福利剧场| 精品国产国语对白av| av线在线观看网站| 看免费成人av毛片| 我的女老师完整版在线观看| 99久国产av精品国产电影| 久久国产精品男人的天堂亚洲 | 亚洲,欧美,日韩| 中文欧美无线码| 精品午夜福利在线看| 国产老妇伦熟女老妇高清| 久久狼人影院| 久久免费观看电影| 亚洲情色 制服丝袜| 国产中年淑女户外野战色| 亚洲精品一二三| 欧美精品一区二区免费开放| 91精品国产国语对白视频| 丝袜在线中文字幕| 久久国产精品大桥未久av | 另类亚洲欧美激情| 哪个播放器可以免费观看大片| av在线播放精品| 纯流量卡能插随身wifi吗| 成人免费观看视频高清| 精品午夜福利在线看| 久久久久精品性色| 十八禁网站网址无遮挡 | 成年av动漫网址| .国产精品久久| 免费不卡的大黄色大毛片视频在线观看| 午夜91福利影院| 国产日韩欧美亚洲二区| 男女啪啪激烈高潮av片| 一区二区三区四区激情视频| tube8黄色片| 国产精品一区二区性色av| 一边亲一边摸免费视频| 久久国产亚洲av麻豆专区| 亚洲欧美成人综合另类久久久| 亚洲美女搞黄在线观看| 麻豆精品久久久久久蜜桃| 国产亚洲av片在线观看秒播厂| 中文字幕亚洲精品专区| 亚洲人成网站在线观看播放| 高清午夜精品一区二区三区| 亚洲精品日韩av片在线观看| 成人特级av手机在线观看| 精品少妇久久久久久888优播| 狂野欧美激情性bbbbbb| 中文在线观看免费www的网站| 中文字幕人妻丝袜制服| 国产日韩一区二区三区精品不卡 | 一级黄片播放器| 欧美激情国产日韩精品一区| 亚洲精品色激情综合| 国产精品久久久久久精品电影小说| av专区在线播放| 自拍偷自拍亚洲精品老妇| 啦啦啦中文免费视频观看日本| 99热这里只有精品一区| 日韩人妻高清精品专区| 午夜福利视频精品| 国产一区有黄有色的免费视频| 亚洲精品aⅴ在线观看| 青春草亚洲视频在线观看| 亚洲,一卡二卡三卡| 亚洲图色成人| 欧美精品人与动牲交sv欧美| 偷拍熟女少妇极品色| 18禁动态无遮挡网站| 国产熟女欧美一区二区| 天堂8中文在线网| 亚洲精品一区蜜桃| 欧美精品一区二区大全| 9色porny在线观看| 婷婷色av中文字幕| 日本黄色日本黄色录像| 国产一区二区在线观看av| 国产精品女同一区二区软件| 另类精品久久| 欧美区成人在线视频| 日本猛色少妇xxxxx猛交久久| 天堂俺去俺来也www色官网| 看十八女毛片水多多多| 97超视频在线观看视频| 最后的刺客免费高清国语| 亚洲精品456在线播放app| 最近中文字幕2019免费版| 大片电影免费在线观看免费| 成人二区视频| 欧美精品人与动牲交sv欧美| 久久精品久久精品一区二区三区| 一级毛片aaaaaa免费看小| 国产精品福利在线免费观看| av国产久精品久网站免费入址| 亚洲av成人精品一区久久| 亚洲精品日韩在线中文字幕| 性色avwww在线观看| 国产在线男女| 国产成人免费观看mmmm| 亚洲精品成人av观看孕妇| 久久鲁丝午夜福利片| 最近最新中文字幕免费大全7| 久久精品夜色国产| 夜夜看夜夜爽夜夜摸| 久久久久人妻精品一区果冻| 在线观看一区二区三区激情| 欧美bdsm另类| 女性被躁到高潮视频| 亚洲性久久影院| 亚洲伊人久久精品综合| 国产黄色视频一区二区在线观看| 黑丝袜美女国产一区| 一边亲一边摸免费视频| 亚洲不卡免费看| 两个人的视频大全免费| 99九九在线精品视频 | 黄色毛片三级朝国网站 | 色视频在线一区二区三区| 久久精品国产a三级三级三级| 边亲边吃奶的免费视频| 好男人视频免费观看在线| 国产伦精品一区二区三区四那| 五月伊人婷婷丁香| 一级黄片播放器| 国产高清不卡午夜福利| 观看美女的网站| 免费黄频网站在线观看国产| 精品久久久久久久久亚洲| 日韩欧美精品免费久久| 99国产精品免费福利视频| 久久久久视频综合| 久久精品久久久久久久性| 大陆偷拍与自拍| 日日撸夜夜添| 特大巨黑吊av在线直播| 大香蕉久久网| 国产亚洲最大av| 高清在线视频一区二区三区| 亚洲精品自拍成人| 在现免费观看毛片| 看免费成人av毛片| 亚洲国产精品一区三区| 老女人水多毛片| 国产一区二区三区av在线| 欧美 亚洲 国产 日韩一| 色5月婷婷丁香| 搡女人真爽免费视频火全软件| 精品人妻熟女av久视频| 亚洲第一区二区三区不卡| 国产精品一区二区在线不卡| 日韩精品免费视频一区二区三区 | 国产欧美亚洲国产| 一级毛片电影观看| 亚洲国产毛片av蜜桃av| 国产一区二区在线观看av| 亚洲国产精品成人久久小说| 成人亚洲欧美一区二区av| 国产精品欧美亚洲77777| 少妇 在线观看| 国产av码专区亚洲av| 免费av中文字幕在线| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕| 中文乱码字字幕精品一区二区三区| 精品酒店卫生间| 国产亚洲欧美精品永久| 伊人久久国产一区二区| 国产一区二区在线观看日韩| 欧美丝袜亚洲另类| 婷婷色麻豆天堂久久| 一区二区三区乱码不卡18| 国产精品99久久久久久久久| av一本久久久久| 久久久a久久爽久久v久久| 日本午夜av视频| 乱系列少妇在线播放| 国产男女超爽视频在线观看| 黄色毛片三级朝国网站 | 男男h啪啪无遮挡| 国产91av在线免费观看| 国产有黄有色有爽视频| 搡女人真爽免费视频火全软件| 黄色怎么调成土黄色| 51国产日韩欧美| 赤兔流量卡办理| 国产亚洲5aaaaa淫片| 久久青草综合色| 久久精品国产a三级三级三级| 免费观看av网站的网址| 菩萨蛮人人尽说江南好唐韦庄| 久久午夜福利片| 国产高清三级在线| 成人毛片60女人毛片免费| 精品少妇久久久久久888优播| 精品亚洲乱码少妇综合久久| 国产精品一区二区在线观看99| 亚洲国产精品一区三区| 精品视频人人做人人爽| 97超视频在线观看视频| 涩涩av久久男人的天堂| 九九久久精品国产亚洲av麻豆| 在线观看www视频免费| 男的添女的下面高潮视频| 国语对白做爰xxxⅹ性视频网站| 人妻一区二区av| av免费在线看不卡| 久久久久久久久久久丰满| 亚洲精品乱码久久久久久按摩| 美女cb高潮喷水在线观看| 伦理电影大哥的女人| 国产精品女同一区二区软件| a级片在线免费高清观看视频| 十分钟在线观看高清视频www | 欧美精品人与动牲交sv欧美| 老司机亚洲免费影院| 成人亚洲欧美一区二区av| 最近中文字幕高清免费大全6| 成人综合一区亚洲| 国产男人的电影天堂91| videos熟女内射| av福利片在线| 婷婷色av中文字幕| 如日韩欧美国产精品一区二区三区 | 人人妻人人澡人人看| 青春草国产在线视频| 国产女主播在线喷水免费视频网站| 日韩成人av中文字幕在线观看| 亚洲人成网站在线播| 国产精品国产三级专区第一集| videossex国产| 男人添女人高潮全过程视频| 男女免费视频国产| 2021少妇久久久久久久久久久| 免费av不卡在线播放| 国产黄片视频在线免费观看| 亚洲在久久综合| 亚洲av电影在线观看一区二区三区| 成人综合一区亚洲| 色婷婷av一区二区三区视频| 国产在线一区二区三区精| 3wmmmm亚洲av在线观看| 亚洲av不卡在线观看| 男女边摸边吃奶| 亚洲精品日本国产第一区| 国产有黄有色有爽视频| 亚洲欧美清纯卡通| 全区人妻精品视频| 夜夜骑夜夜射夜夜干| 黄色日韩在线| 丰满少妇做爰视频| 五月开心婷婷网| 久久青草综合色| 伦精品一区二区三区| 中文精品一卡2卡3卡4更新| 老女人水多毛片| 日本vs欧美在线观看视频 | 亚洲av福利一区| 亚洲国产精品国产精品| 久久久久久久久大av| 亚洲图色成人| 久久久久久久久大av| 久久精品国产a三级三级三级| 日韩伦理黄色片| 晚上一个人看的免费电影| 亚洲真实伦在线观看| 午夜激情福利司机影院| 在线观看av片永久免费下载| 久久精品夜色国产| 亚洲第一区二区三区不卡| 久久午夜综合久久蜜桃| 日本爱情动作片www.在线观看| 一边亲一边摸免费视频| 欧美xxxx性猛交bbbb| 黄色配什么色好看| 亚洲一级一片aⅴ在线观看| 国产精品99久久99久久久不卡 | 多毛熟女@视频| 久久人人爽人人片av| 欧美日本中文国产一区发布| 十分钟在线观看高清视频www | a级片在线免费高清观看视频| 亚洲怡红院男人天堂| 老司机影院成人| 成年人免费黄色播放视频 | 日本黄大片高清| 伦理电影免费视频| 亚洲不卡免费看| 日产精品乱码卡一卡2卡三| 精品国产国语对白av| 91精品一卡2卡3卡4卡| 一本色道久久久久久精品综合| 亚洲精品第二区| 亚洲国产日韩一区二区| 大码成人一级视频| 一级黄片播放器| 亚洲人成网站在线观看播放| 国产一区二区在线观看av| 香蕉精品网在线| 精品人妻偷拍中文字幕| 亚洲情色 制服丝袜| 亚洲av成人精品一区久久| 人人妻人人澡人人看| 久久99热6这里只有精品| 黄色日韩在线| 欧美日韩视频高清一区二区三区二| 丁香六月天网| 99热国产这里只有精品6| 在线 av 中文字幕| 熟女电影av网| 91精品伊人久久大香线蕉| 国产欧美日韩精品一区二区| 日日啪夜夜撸| 免费黄色在线免费观看| 精品国产一区二区三区久久久樱花| 一区二区三区四区激情视频| 亚洲欧美成人精品一区二区| 一级a做视频免费观看| 精品一区二区三卡| 亚洲精品乱久久久久久| 亚洲精品456在线播放app| 蜜臀久久99精品久久宅男| 丝袜脚勾引网站| 亚洲欧美一区二区三区黑人 | 97精品久久久久久久久久精品| 欧美日韩综合久久久久久| 日韩成人av中文字幕在线观看| 青春草视频在线免费观看| 精品国产一区二区久久| 性色av一级| 久久免费观看电影| 色婷婷久久久亚洲欧美| 国产av精品麻豆| 久久97久久精品| 国产成人一区二区在线| 国产国拍精品亚洲av在线观看| 少妇熟女欧美另类| 你懂的网址亚洲精品在线观看| 在线观看免费高清a一片| 中文资源天堂在线| 97在线视频观看| 午夜福利视频精品| 久久婷婷青草| 欧美亚洲 丝袜 人妻 在线| 美女大奶头黄色视频| 午夜福利视频精品| 久久亚洲国产成人精品v| 亚洲国产av新网站| 久热久热在线精品观看| 欧美三级亚洲精品| 国产高清有码在线观看视频| 久久狼人影院| 少妇人妻 视频| 欧美日韩av久久| 国产在线视频一区二区| 色哟哟·www| 又黄又爽又刺激的免费视频.| 日韩av不卡免费在线播放| 五月天丁香电影| 男女免费视频国产| 2021少妇久久久久久久久久久| 女性生殖器流出的白浆| 午夜免费鲁丝| av一本久久久久| 狠狠精品人妻久久久久久综合| 国产极品天堂在线| videossex国产| 女性生殖器流出的白浆| 久久久久视频综合| 99热6这里只有精品| 全区人妻精品视频| 亚洲欧美日韩卡通动漫| 国产毛片在线视频| √禁漫天堂资源中文www| 国产精品国产三级国产专区5o| 国产男女超爽视频在线观看| 亚洲,一卡二卡三卡| 又黄又爽又刺激的免费视频.| 一区二区三区精品91| 一级黄片播放器| 亚洲av中文av极速乱| 国产精品99久久99久久久不卡 | 97在线人人人人妻| 最新中文字幕久久久久| 国产一级毛片在线| 欧美精品一区二区免费开放| 日韩一区二区三区影片| 亚洲精品日韩在线中文字幕| 51国产日韩欧美| 伊人亚洲综合成人网| 简卡轻食公司| av网站免费在线观看视频| 亚洲国产精品成人久久小说| av天堂久久9| 欧美日韩国产mv在线观看视频| 国产成人一区二区在线| 97在线视频观看| 男女无遮挡免费网站观看| 啦啦啦在线观看免费高清www| 成人综合一区亚洲| 蜜臀久久99精品久久宅男| 亚洲国产精品成人久久小说| 国产女主播在线喷水免费视频网站| 免费观看的影片在线观看| 插逼视频在线观看| 青青草视频在线视频观看| 国产在线视频一区二区| 在线观看免费高清a一片| 久久久久国产网址| 久久午夜综合久久蜜桃| 午夜免费男女啪啪视频观看| 久久人人爽人人爽人人片va| a级毛色黄片| 少妇被粗大猛烈的视频| 99热6这里只有精品| 久久久亚洲精品成人影院| 一本一本综合久久| 夜夜看夜夜爽夜夜摸| 成年人午夜在线观看视频| 在线看a的网站| 2021少妇久久久久久久久久久| av.在线天堂| 久久国产亚洲av麻豆专区| .国产精品久久| 嫩草影院入口| 女人久久www免费人成看片| 午夜91福利影院| 亚洲av日韩在线播放| 在现免费观看毛片| 我的老师免费观看完整版| 我要看黄色一级片免费的| 国产精品一区二区性色av| 日韩欧美精品免费久久| 日本欧美视频一区| 免费av中文字幕在线| 男人添女人高潮全过程视频| 18禁裸乳无遮挡动漫免费视频| 伊人久久国产一区二区| 九九久久精品国产亚洲av麻豆| 乱人伦中国视频| 美女国产视频在线观看| 亚洲色图综合在线观看| 欧美丝袜亚洲另类| 纵有疾风起免费观看全集完整版| 日日爽夜夜爽网站| 久久久久国产网址| 毛片一级片免费看久久久久| 免费av不卡在线播放| 大片免费播放器 马上看| 美女脱内裤让男人舔精品视频| 国模一区二区三区四区视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲四区av| 午夜福利网站1000一区二区三区| 国产精品久久久久久精品电影小说| 婷婷色麻豆天堂久久| 嘟嘟电影网在线观看| 亚洲精品,欧美精品| 亚洲精品国产av蜜桃| 一区在线观看完整版| 如何舔出高潮| 日韩精品有码人妻一区| videos熟女内射| 国产亚洲欧美精品永久| 18+在线观看网站| 欧美性感艳星| 欧美人与善性xxx| av国产精品久久久久影院| 中文在线观看免费www的网站| 国产精品国产三级国产av玫瑰| videos熟女内射| 少妇的逼好多水| 亚洲无线观看免费| 在线观看免费高清a一片| 久久女婷五月综合色啪小说| 男人和女人高潮做爰伦理| 亚洲精品乱码久久久久久按摩| 香蕉精品网在线| 亚洲一级一片aⅴ在线观看| 赤兔流量卡办理| 国精品久久久久久国模美| 国产有黄有色有爽视频| 精品人妻熟女av久视频| 欧美日韩视频精品一区| 成年av动漫网址| av又黄又爽大尺度在线免费看| 成人综合一区亚洲| 黑人猛操日本美女一级片| 国产伦在线观看视频一区| 亚洲国产精品一区三区| 国产白丝娇喘喷水9色精品| 热99国产精品久久久久久7| 亚洲av欧美aⅴ国产| 国内精品宾馆在线| 亚洲国产欧美日韩在线播放 | 91在线精品国自产拍蜜月| 国产高清不卡午夜福利| 亚洲欧美清纯卡通| 天天操日日干夜夜撸| 黄色视频在线播放观看不卡| 日本与韩国留学比较| 老熟女久久久| 国产一区二区三区av在线| 日本猛色少妇xxxxx猛交久久| 精品一品国产午夜福利视频| 人人妻人人添人人爽欧美一区卜| 久久久久久久久久久免费av| 欧美日韩国产mv在线观看视频| 边亲边吃奶的免费视频| 久久精品久久久久久久性| 美女国产视频在线观看| 成人美女网站在线观看视频| 亚洲精品国产色婷婷电影| 久久久国产精品麻豆| 国产精品久久久久久精品电影小说| 欧美日韩在线观看h| 高清毛片免费看| 五月玫瑰六月丁香| 秋霞伦理黄片| 校园人妻丝袜中文字幕| 色94色欧美一区二区| 欧美xxxx性猛交bbbb| 亚洲欧美精品专区久久| av线在线观看网站|