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

    能源樁熱-力學(xué)特性模型試驗(yàn)與數(shù)值模擬

    2022-06-14 09:13:36常虹朱萬(wàn)里王琰趙嵩穎
    科學(xué)技術(shù)與工程 2022年13期
    關(guān)鍵詞:樁體黏土阻力

    常虹, 朱萬(wàn)里, 王琰, 趙嵩穎

    (1.吉林建筑大學(xué)測(cè)繪與勘查工程學(xué)院, 長(zhǎng)春 130118; 2. 吉林建筑大學(xué)市政與環(huán)境工程學(xué)院, 長(zhǎng)春 130118)

    中國(guó)北方地區(qū)在冬季采暖所用的化石燃料燃燒是引起大氣質(zhì)量惡化的重要因素之一,利用可再生的清潔能源作為空調(diào)系統(tǒng)的主要能量來(lái)源是解決這一問(wèn)題的有效途徑。地源熱泵(ground source heat pumps,GSHP)技術(shù)的基本原理是利用埋設(shè)于地下的換熱管路實(shí)現(xiàn)建筑與地下相對(duì)恒溫層的熱交換,但由于其占用土地面積、鉆孔造價(jià)昂貴等缺點(diǎn)具有一定局限性。20世紀(jì)80年代,奧地利及瑞典等國(guó)工程師提出將換熱管與建筑物的混凝土樁基相結(jié)合,由此誕生能源樁的理念,又名能量樁或樁基埋管[1-2]。

    由于在傳統(tǒng)樁基礎(chǔ)中引入溫度場(chǎng),能源樁運(yùn)行過(guò)程中會(huì)引起樁周土溫度場(chǎng)的變化,樁身會(huì)產(chǎn)生附加應(yīng)力和變形,從而影響樁基承載性能。針對(duì)能源樁的研究需著重考慮如何克服循環(huán)變化的溫度場(chǎng)對(duì)樁基力學(xué)特性的影響。近年來(lái),針對(duì)能源樁的熱-力學(xué)特性,眾多學(xué)者通過(guò)試驗(yàn)與數(shù)值模擬等方法開(kāi)展系列研究并取得一定成果。Gunawan等[3]研究不同超固結(jié)比黏土地基中能源樁的承載力特性,得出隨著超固結(jié)比的增大樁頂位移逐漸減小的結(jié)論;王成龍等[4]開(kāi)展室內(nèi)模型試驗(yàn),對(duì)飽和砂土地基中能源樁熱-力耦合作用時(shí)的樁身承載力和傳熱機(jī)理進(jìn)行研究,得出隨溫度循環(huán)次數(shù)的增加樁頂產(chǎn)生累積沉降的結(jié)論;Yavari等[5]對(duì)飽和黏土地基中能源樁力學(xué)性能進(jìn)行研究,得出隨著工作荷載的增大黏土蠕變速率增加,多次溫度循環(huán)導(dǎo)致樁體產(chǎn)生不可逆的累積沉降的結(jié)論;郭浩然等[6]基于彈塑性理論建立能源樁數(shù)值模型,主要對(duì)能源樁運(yùn)行時(shí)的受力特性進(jìn)行分析,得出循環(huán)溫度對(duì)樁身的影響隨樁頂荷載的增大而減小的結(jié)論;黃旭等[7]開(kāi)展室內(nèi)模型試驗(yàn),針對(duì)不封底現(xiàn)澆鋼筋混凝土大直徑管(cast-in-place concrete large-diameter pipe,PCC)能源樁運(yùn)行過(guò)程中的熱傳導(dǎo)特性進(jìn)行研究,得出該樁型在夏季的換熱效率顯著高于普通樁型的結(jié)論;Wu[8]研究了飽和黏土地基中能源樁、相鄰的普通混凝土樁和承臺(tái)之間的相互作用,得出鄰樁和承臺(tái)可限制能源樁沉降的結(jié)論;費(fèi)康等[9]開(kāi)展模型試驗(yàn)與數(shù)值模擬,對(duì)正常固結(jié)黏土地基中能源樁長(zhǎng)期運(yùn)行的熱-力學(xué)特性進(jìn)行研究,得出樁頂沉降隨溫度循環(huán)次數(shù)的增多而不斷累積的結(jié)論;李小敏[10]研究了復(fù)合材料混凝土能源樁的承載力特性及傳熱效率,同時(shí)基于FLUENT軟件開(kāi)展數(shù)值模擬研究,得出銅礦渣細(xì)骨料大功率充電(high power charging,HPC)能源樁相較于傳統(tǒng)能源樁在換熱效率和承載力特性方面有顯著優(yōu)勢(shì);高磊等[11]采用多場(chǎng)耦合軟件建立能源樁三維數(shù)值模型,研究螺旋形埋管能源樁的熱-力學(xué)特性,得出樁身應(yīng)變隨樁長(zhǎng)縮短而線性減小等結(jié)論;楊濤等[12]采用數(shù)值模擬方法對(duì)熱-力耦合作用下的能源樁承載力特性進(jìn)行研究,設(shè)置溫度荷載按正弦函數(shù)變化,得出樁周土相較于樁身溫度變化較為滯后等結(jié)論;史宏財(cái)?shù)萚13]基于圓柱源換熱理論建立非飽和土地基中能源樁數(shù)值模型,提出樁周土飽和度與能源樁換熱比率的函數(shù)關(guān)系,并開(kāi)展室內(nèi)試驗(yàn)驗(yàn)證了模型適用性;趙蕾等[14]基于COMSOL軟件建立不同埋管形式能源樁的數(shù)值模型,對(duì)比研究其換熱性能、軸力、樁身位移等方面的不同,得出雙螺旋形埋管的能源樁換熱效率較高、W形埋管能源樁的樁身附加溫度荷載較大的結(jié)論。

    越來(lái)越多的學(xué)者將數(shù)值模擬方法運(yùn)用到能源樁研究領(lǐng)域,運(yùn)用有限元數(shù)值模擬軟件計(jì)算得出的結(jié)果可以一定程度上與試驗(yàn)結(jié)果相互印證。近年來(lái),針對(duì)能源樁熱-力學(xué)特性的研究主要集中于砂土和非飽和黏性土地基,而針對(duì)飽和黏土地基開(kāi)展的研究相對(duì)較少。鑒于飽和黏土具有高含水量、低強(qiáng)度、高壓縮性、低滲透性等特點(diǎn),且溫度改變引起的飽和黏土變形會(huì)改變樁土的位移模式,進(jìn)而改變樁側(cè)阻力和端阻力的發(fā)揮,現(xiàn)結(jié)合室內(nèi)模型試驗(yàn)與數(shù)值模擬,展開(kāi)飽和黏土地基能源樁在冷熱交替循環(huán)過(guò)程中的力學(xué)特性研究,為能源樁在飽和黏土地基的應(yīng)用給出建議。

    1 試驗(yàn)設(shè)計(jì)

    1.1 試驗(yàn)系統(tǒng)

    飽和黏土地基中能源樁模型試驗(yàn)系統(tǒng)包括模型槽、量測(cè)系統(tǒng)、換熱系統(tǒng)和加載設(shè)備。模型槽是由聚甲基丙烯酸甲酯板材圍成的直徑550 mm的圓桶,該材料具有延展性好、抗沖擊力強(qiáng)等特點(diǎn)。緊貼模型槽內(nèi)壁覆有透明塑料薄膜,以防止內(nèi)部土體水分流失?;鶞?zhǔn)梁橫跨模型槽兩側(cè),用以固定數(shù)顯百分表,測(cè)量樁頂位移。假定基準(zhǔn)梁與大地變形協(xié)調(diào),不考慮基準(zhǔn)梁形變對(duì)數(shù)據(jù)的影響。量測(cè)系統(tǒng)包括DH3818靜態(tài)應(yīng)變測(cè)試系統(tǒng)、120-50AA混凝土應(yīng)變片(量程0~20 000 με,1/4橋式)、電阻式微型土壓力盒、PT100鉑熱電阻溫度傳感器、溫度采集儀、數(shù)顯式百分表(精度0.01 mm)等。應(yīng)變片和土壓力盒與DH3818相連,溫度傳感器與溫度采集儀相連,可自定義采集頻率,數(shù)據(jù)自動(dòng)采集。換熱系統(tǒng)由數(shù)顯恒溫水浴鍋及全自動(dòng)自吸水泵組成,導(dǎo)熱液體為水,循環(huán)流速為22 L/min。加載設(shè)備為鐵制砝碼,工作荷載與樁頂之間墊有剛性加載板。

    1.2 樁周土體參數(shù)

    試驗(yàn)用土為飽和黏土,其物理力學(xué)參數(shù)如表1所示。

    表1 飽和黏土物理性質(zhì)Table 1 The physical properties of saturated clay

    1.3 模型樁制作

    試驗(yàn)所用模型樁采用C30強(qiáng)度的混凝土澆筑,模型樁直徑D=84 mm,樁長(zhǎng)L=500 mm,主筋由3根HRB400C8 mm鋼筋組成,箍筋為B6 mm@100 mm。樁內(nèi)換熱水管為PVC鋼絲軟管,水管外徑為14 mm,內(nèi)徑8 mm,試驗(yàn)采用單U形埋管方式,U形兩肢相距32 mm。

    1.4 試驗(yàn)方案

    1.4.1 準(zhǔn)備工作

    試驗(yàn)采用亞克力板圍成的直徑550 mm的模型槽,槽內(nèi)土體高度H=500 mm,分3層填筑,每層虛高20 cm,并以18 kPa荷載靜壓48 h再進(jìn)行下一層填筑,每填筑一層適當(dāng)對(duì)土體噴水養(yǎng)護(hù)以保證土體始終處于飽和狀態(tài)。樁端持力層填筑完成后,在樁底位置埋設(shè)土壓力盒。土體填筑結(jié)束后,將模型地基靜壓一周,以形成正常固結(jié)狀態(tài)的飽和黏土地基。

    量測(cè)系統(tǒng)包含的各傳感器布設(shè)如圖1所示,樁頂對(duì)稱布置兩個(gè)百分表B1、B2以測(cè)定樁頂豎向位移,讀數(shù)取二者平均值;在距離土體表面5、200、400 mm深度分別布置3層、每層布置4個(gè)共計(jì)12個(gè)溫度傳感器(T1~T12),同一深度處傳感器分別距樁軸線42 mm(1/2D)、84 mm(1D)、168 mm(2D)、252 mm(3D);土體表面距樁側(cè)30、130、230 mm處分別布設(shè)3個(gè)百分表B3、B4、B5,以測(cè)定土體沉降;土壓力盒布置在樁端下,距離模型槽底部10 cm,以測(cè)量樁端阻力;在樁身兩側(cè)對(duì)稱粘貼應(yīng)變片,每側(cè)4片等距布置,共計(jì)8片。

    圖1 測(cè)點(diǎn)布置圖Fig.1 Layout of measuring points

    1.4.2 能源樁靜載荷試驗(yàn)

    首先按照規(guī)范估算模型樁的極限承載力,隨后按照計(jì)算所得的極限承載力分10級(jí)加載。靜載試驗(yàn)采用慢速維持荷載法,加載過(guò)程中記錄樁頂沉降,當(dāng)每級(jí)荷載下樁頂沉降量小于0.1 mm/h時(shí),則認(rèn)為已趨于穩(wěn)定,可施加下一級(jí)荷載。當(dāng)某級(jí)荷載下樁頂沉降量為前一級(jí)荷載下沉降量的5倍,則立即終止加載。根據(jù)《建筑基樁檢測(cè)技術(shù)規(guī)范》(JGJ 106—2014)[15]的規(guī)定,荷載-沉降曲線上出現(xiàn)陡降段時(shí),單樁極限承載力取陡降段起點(diǎn)的荷載值為1.0 kN,取工作荷載為0.5 kN。

    1.4.3 試驗(yàn)工況

    根據(jù)溫度荷載大小、循環(huán)次數(shù)及樁頂荷載情況將試驗(yàn)分為5種工況,具體見(jiàn)表2。

    表2 試驗(yàn)工況Table 2 Test conditions

    設(shè)定工況一~工況三初始溫度為20 ℃以模擬室溫,3種工況下分別將換熱液體加熱至50、60、70 ℃,升溫過(guò)程約為15 min,達(dá)到目標(biāo)溫度后維持24 h,隨后將循環(huán)液體迅速冷卻至初始溫度并維持5 h。工況四和工況五將溫度循環(huán)設(shè)定為5 ℃→70 ℃→5 ℃,升溫過(guò)程約為30 min,達(dá)到70 ℃后維持24 h,隨后將循環(huán)液體迅速冷卻至5 ℃并維持5 h,如此循環(huán)3次。

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

    2.1 樁身和土體溫度變化規(guī)律

    在數(shù)據(jù)處理時(shí)發(fā)現(xiàn)T1測(cè)點(diǎn)的熱響應(yīng)較迅速,故選取工況一~工況三的T1測(cè)點(diǎn)數(shù)據(jù),分析循環(huán)水溫度的高低對(duì)樁身溫度的影響,從圖2可以看出,初始階段樁身溫度升高較明顯,超過(guò)5 h后逐漸趨于穩(wěn)定,3種工況下樁身最高溫度分別為30.0、34.8、38.5 ℃,可知樁身溫度隨入水溫度的升高而升高。

    圖3為工況四所得溫度隨時(shí)間變化曲線,可以看出,樁身各點(diǎn)溫度在初始階段升溫較快,超過(guò)5 h后逐漸趨于穩(wěn)定,測(cè)點(diǎn)T1、T5、T9的最高溫度分別為38.8、35.1、25.6 ℃,沿深度方向樁身溫度逐漸減小,隨循環(huán)次數(shù)的增加,樁身最高溫度先降低后升高,考慮是樁周土體溫度改變較為滯后,樁-土換熱量先增加后減小的緣故。

    圖4為距離填土表面200 mm深度處土體沿徑向的溫度變化,可以看出,距離樁越遠(yuǎn)土體溫度越低,遠(yuǎn)端T8點(diǎn)的溫度變化不明顯;隨循環(huán)次數(shù)的增加,土體溫度曲線逐漸變緩,溫度逐漸降低。

    圖2 不同循環(huán)溫度下樁身溫度變化曲線Fig.2 The change curves of pile temperature at different temperature cycles

    圖3 工況四樁身溫度隨時(shí)間變化曲線Fig.3 The change curves of pile temperature with time in the fourth condition

    圖4 工況四土體溫度隨時(shí)間變化曲線Fig.4 The change curves of soil temperature with time in the fourth condition

    2.2 樁頂和土體位移

    2.2.1 樁頂位移

    工況四多次溫度循環(huán)下樁頂沉降變化如圖5所示,規(guī)定向上的位移為正,向下的位移為負(fù),可以看出,升溫階段樁體發(fā)生膨脹,樁頂產(chǎn)生向上的位移,最大上升位移為0.089 mm;降溫時(shí)樁體收縮位移向下,最大沉降為0.052 mm,降溫所導(dǎo)致的沉降量大于升溫的膨脹量,考慮到樁端黏土由于樁體膨脹產(chǎn)生變形所致,隨著循環(huán)次數(shù)的增加,樁體的不可恢復(fù)沉降將逐漸累積,三次溫度循環(huán)后樁頂產(chǎn)生累積沉降量達(dá)-0.052 mm(0.6%D),此結(jié)論與王成龍等[4]、Yavari等[5]、Liao等[16]的結(jié)論一致,工程上須保證能源樁長(zhǎng)期運(yùn)行的穩(wěn)定性,而樁頂位移尤其是沉降位移是影響結(jié)構(gòu)安全的重要因素,尤其對(duì)于飽和黏土地基多次溫度循環(huán)后的樁頂累積沉降應(yīng)引起足夠重視。

    2.2.2 土體位移

    選取工況四土體表面豎向位移數(shù)據(jù)繪制位移-時(shí)間曲線,如圖6所示。百分表B3靠近樁身,由于升溫階段樁體膨脹產(chǎn)生擠土效應(yīng),土體表面隆起,故表現(xiàn)為較大的上升位移,不予考慮。升溫階段,由于土體發(fā)生熱固結(jié)現(xiàn)象產(chǎn)生沉降;降溫會(huì)加劇土體沉降,考慮是由于超靜孔隙水壓力的消散使土體固結(jié)程度增大的緣故[17]。受土體溫度場(chǎng)的影響,靠近樁身土體的沉降速率大于遠(yuǎn)端。隨循環(huán)次數(shù)的增加,土體沉降速率呈減小趨勢(shì),三次溫度循環(huán)后,B4、B5處土體的最終沉降量分別為-1.19 mm和-0.61 mm。因此,能源樁應(yīng)用于飽和黏土地基時(shí),設(shè)計(jì)階段需考慮到黏土發(fā)生熱固結(jié)現(xiàn)象,產(chǎn)生收縮變形對(duì)樁基承載力的影響。

    圖5 樁頂位移變化曲線Fig.5 The curve of displacement of pile top

    圖6 土體表面豎向位移變化曲線Fig.6 The curves of vertical displacement of soil surface

    3 試驗(yàn)與模擬力學(xué)特性對(duì)比分析

    3.1 三維數(shù)值模型建立

    基于COMSOL Multiphysics軟件建立與能源樁試驗(yàn)裝置等大的三維數(shù)值模型,利用軟件內(nèi)置的巖土工程模塊和偏微分方程(partial differential equation,PDE)引入物理力學(xué)場(chǎng)和溫度場(chǎng),為優(yōu)化求解過(guò)程,在建立模型之初做出如下假設(shè):①默認(rèn)樁周土壤中無(wú)地下水滲流運(yùn)動(dòng);②樁身及樁周土材料的物理力學(xué)參數(shù)保持恒定,不隨溫度場(chǎng)而變化;③將換熱管及換熱液體簡(jiǎn)化為線性熱源;④數(shù)值模型的初始溫度統(tǒng)一且均勻分布。

    考慮到計(jì)算難度隨網(wǎng)格劃分細(xì)密程度線性增長(zhǎng),對(duì)換熱管采用邊單元網(wǎng)格劃分法,最小單元尺寸0.12 mm,樁身及樁周土體采用自由四面體網(wǎng)格劃分法,最小單元尺寸分別為0.9 mm和2.4 mm,如圖7所示。

    樁身應(yīng)力測(cè)點(diǎn)分布越密集,越利于分析樁身力學(xué)性質(zhì),在數(shù)值模擬中沿樁身等距取9個(gè)測(cè)點(diǎn),相鄰測(cè)點(diǎn)間距為50 mm,運(yùn)用COMSOL軟件內(nèi)置的巖土工程模塊計(jì)算樁身各項(xiàng)力學(xué)指標(biāo)。

    圖7 模型網(wǎng)格劃分Fig.7 Model meshing

    3.2 樁端土壓力

    樁端阻力的變化表現(xiàn)為樁端土壓力的變化,5種工況下樁端土壓力隨時(shí)間變化的曲線如圖8所示。加熱和制冷階段開(kāi)始后,各工況下樁體的樁端土壓力在2 h內(nèi)分別迅速增加和減小,考慮是由于樁體溫度未達(dá)到穩(wěn)定狀態(tài),2 h之后樁端土壓力逐漸趨于穩(wěn)定。由于循環(huán)溫度的差別,樁端土壓力的大小也各不相同。

    工況一~工況三加熱階段結(jié)束后,樁端出現(xiàn)的最大壓應(yīng)力分別為0.064 1、0.086 1、0.109 5 MPa;制冷結(jié)束后,樁端最小壓應(yīng)力分別為0.017 6、0.020 6、0.021 8 MPa。加熱后樁體受熱膨脹,樁端土壓力升高,而制冷時(shí)樁體受冷收縮,樁端土壓力逐漸減少,這與 Yavari 等[5]的結(jié)論基本一致。循環(huán)水的溫度越高,樁端土壓力越大;加熱與制冷階段的溫差越大,樁端最大壓應(yīng)力與最小壓應(yīng)力的差值越大。由工況四和工況五的土壓力變化曲線可以看出,工作荷載作用下,加熱結(jié)束時(shí)的樁端壓應(yīng)力約為零荷載作用下的1.65倍,制冷結(jié)束時(shí)的樁端壓應(yīng)力約為零荷載作用下的14倍,變化規(guī)律與王成龍等[4]基本一致。加熱階段結(jié)束后,樁端出現(xiàn)的最大壓應(yīng)力分別為0.107 6 MPa和0.178 4 MPa;制冷結(jié)束后,樁端最小壓應(yīng)力分別為0.005 9 MPa和0.084 4 MPa。隨循環(huán)次數(shù)的增加,樁端最大壓應(yīng)力逐漸減小,考慮是由于樁周土體在多次冷熱交替后發(fā)生熱固結(jié)作用,土體強(qiáng)度增大的緣故。

    數(shù)值模擬所得樁端土壓力隨能源樁運(yùn)行時(shí)間的變化趨勢(shì)與試驗(yàn)結(jié)果基本一致,數(shù)值相對(duì)較大。數(shù)值模擬樁端土壓力平均值相較于試驗(yàn)結(jié)果偏大20%~25%,考慮是由于建立模型時(shí)簡(jiǎn)化荷載傳遞公式的緣故。

    圖8 試驗(yàn)與模擬樁端土壓力對(duì)比圖Fig.8 Comparison of pile end earth pressure between test and simulation

    3.3 樁身側(cè)摩阻力分布規(guī)律

    溫度循環(huán)引起的樁體膨脹和收縮受到約束時(shí),在樁體內(nèi)部會(huì)產(chǎn)生軸向附加應(yīng)力,此時(shí)樁身實(shí)際熱應(yīng)變小于自由應(yīng)變?chǔ)臫-free,產(chǎn)生熱應(yīng)力σT,進(jìn)而引起樁側(cè)摩阻力的改變,不同深度處的樁側(cè)摩阻力fs,mob,j計(jì)算公式推導(dǎo)過(guò)程如下。

    εT-free=αcΔT

    (1)

    σT=EεT-rstr=E(εT-free-εT-obs)

    (2)

    fs,mob,j=(σT,j-σT,j-1)D/(4Δl)

    (3)

    式中:σT為溫度荷載下樁身的附加熱應(yīng)力,MPa;E為彈性模量,GPa;εT-rstr為溫度荷載下樁身的附加熱應(yīng)變;εT-free為溫度荷載下樁身的自由應(yīng)變;εT-obs為溫度荷載下樁身的觀測(cè)應(yīng)變;αc為樁體熱膨脹系數(shù),樁體熱膨脹系數(shù)采用 10 με/℃;ΔT為測(cè)點(diǎn)溫差,℃;D為樁體直徑;Δl為相鄰應(yīng)變片的間距;j=1,2,3,4表示從土表面到樁端應(yīng)變片編號(hào)。定義側(cè)摩阻力向上為正,向下為負(fù)。

    工況一~工況三中24 h和29 h的樁身側(cè)摩阻力-深度分布曲線如圖9所示。對(duì)比分析模型試驗(yàn)數(shù)據(jù)與數(shù)值模擬數(shù)據(jù)可發(fā)現(xiàn)二者側(cè)摩阻力趨勢(shì)基本一致,樁體受熱時(shí)兩端分別向上和向下運(yùn)動(dòng),樁體上半部分產(chǎn)生負(fù)的側(cè)摩阻力,下半部分產(chǎn)生正的側(cè)摩阻力;樁身側(cè)摩阻力隨循環(huán)溫度的升高而逐漸增大;無(wú)論升溫階段還是降溫階段,位移零點(diǎn)均產(chǎn)生在樁身中部偏下位置;樁端附近產(chǎn)生的側(cè)摩阻力數(shù)值均小于樁體上半部分側(cè)摩阻力數(shù)值,考慮因?yàn)闃抖送恋募s束作用使樁-土相對(duì)位移較小的緣故。

    取工況四和工況五24 h和29 h的樁體側(cè)摩阻力進(jìn)行對(duì)比,以此分析工作荷載作用對(duì)側(cè)摩阻力的影響,試驗(yàn)與模擬結(jié)果對(duì)比如圖10所示。樁體升溫時(shí),樁體上部側(cè)摩阻力為負(fù),下部側(cè)摩阻力為正,降溫時(shí)恰好相反。模型試驗(yàn)升溫結(jié)束時(shí),零荷載和工作荷載作用下最大側(cè)摩阻力值分別為0.062 1 MPa和0.046 0 MPa;29 h時(shí)刻,二者最大側(cè)摩阻力值分別為0.029 9 MPa和0.047 3 MPa。工作荷載作用使樁身整體下沉,產(chǎn)生負(fù)摩阻力的區(qū)域逐漸變小,相較于零荷載作用時(shí),位移零點(diǎn)上移至樁中部偏上位置,此結(jié)論與數(shù)值模擬結(jié)果一致。

    選取工況五升溫和降溫結(jié)束時(shí)的樁側(cè)摩阻力數(shù)據(jù)繪制其沿深度分布圖,以分析溫度循環(huán)次數(shù)對(duì)樁體側(cè)摩阻力的影響,如圖11所示??梢钥闯?,無(wú)論升溫或降溫階段,靠近土體表面處的模擬結(jié)果均明顯小于試驗(yàn)結(jié)果,考慮是由于兩處位于溫度場(chǎng)的邊緣,溫度變化幅度較小且遲緩,不會(huì)產(chǎn)生明顯的熱量堆積效應(yīng),樁身附加熱應(yīng)力相對(duì)較小的緣故。因此,在設(shè)計(jì)模型試驗(yàn)時(shí)需注意在樁身兩端細(xì)分應(yīng)力測(cè)點(diǎn)。模型試驗(yàn)與數(shù)值模擬的樁身中下部側(cè)摩阻力值無(wú)較大差異,試驗(yàn)與模擬結(jié)果的最大側(cè)摩阻力值基本相同。模型試驗(yàn)中3次循環(huán)升溫階段的最大側(cè)摩阻力值依次為0.046 0、0.049 1、0.050 6 MPa;降溫時(shí)最大側(cè)摩阻力值依次為0.047 3、0.052 3、0.052 9 MPa。多次溫度循環(huán)后,降溫時(shí)位移零點(diǎn)下移至樁體中部偏下位置,升溫階段位移零點(diǎn)的位置始終保持在樁體上部。隨著循環(huán)次數(shù)的增加,側(cè)摩阻力逐漸增大,增量逐漸減小,考慮是多次冷、熱循環(huán)使樁側(cè)土體產(chǎn)生固結(jié)變形對(duì)樁體的約束作用增強(qiáng)的緣故。

    圖9 工況一~工況三側(cè)摩阻力沿深度分布Fig.9 The shaft friction distribution along depth in the first to third conditions

    圖10 工況四、工況五側(cè)摩阻力沿深度分布Fig.10 The shaft friction distribution along depth in the fourth and fifth conditions

    圖11 工況五試驗(yàn)與模擬側(cè)摩阻力對(duì)比圖Fig.11 The shaft friction distribution of test and simulation in the fifth condition

    4 結(jié)論

    通過(guò)開(kāi)展不同工況下飽和黏土地基中能源樁室內(nèi)模型試驗(yàn)和數(shù)值模擬研究,將模擬結(jié)果與試驗(yàn)結(jié)果相互印證,獲得樁土溫度、位移和樁基承載力及其變化規(guī)律,得到如下結(jié)論。

    (1)對(duì)換熱液體加熱時(shí),前5 h樁身溫度迅速升高,隨后逐漸趨于穩(wěn)定。降溫階段所引起的樁頂沉降量大于升溫產(chǎn)生的膨脹量,且多次溫度循環(huán)導(dǎo)致的累積變形可能會(huì)對(duì)上部建筑的安全造成影響,在設(shè)計(jì)階段需引起重視。多次溫度循環(huán)使土體發(fā)生熱固結(jié)現(xiàn)象,樁周土體表現(xiàn)為沉降,靠近樁身土體的沉降速率大于遠(yuǎn)端。隨循環(huán)次數(shù)的增加,樁周土體的沉降速率呈減小趨勢(shì)。

    (2)樁身側(cè)摩阻力隨溫度的升高而逐漸增大,任一工況樁端附近產(chǎn)生的側(cè)摩阻力數(shù)值均小于樁頂部分的側(cè)摩阻力數(shù)值,且側(cè)摩阻力隨循環(huán)次數(shù)的增加而逐漸增大。工作荷載的加入導(dǎo)致樁身產(chǎn)生負(fù)摩阻力的區(qū)域逐漸變小,位移零點(diǎn)逐漸上移。

    (3)采用數(shù)值模擬研究方法可以較好地模擬能源樁溫度循環(huán)下的承載力特性,數(shù)值模擬結(jié)果可與模型試驗(yàn)結(jié)果做到相互印證。數(shù)值模擬所得的側(cè)摩阻力在靠近土體表面處和樁端處明顯小于模型試驗(yàn)后果,在設(shè)計(jì)試驗(yàn)時(shí)需注意增加在樁身兩端的測(cè)點(diǎn)以便于更完整、細(xì)致地分析能源樁側(cè)摩阻力趨勢(shì)。

    猜你喜歡
    樁體黏土阻力
    盾構(gòu)隧道近接側(cè)穿引起橋梁樁體的變形規(guī)律*
    樁體模量和加筋體剛度對(duì)路堤穩(wěn)定性影響分析
    地鐵車站支撐拆除對(duì)圍護(hù)結(jié)構(gòu)的影響
    鼻阻力測(cè)定在兒童OSA診療中的臨床作用
    零阻力
    不一般的黏土插畫(huà)
    黏土多肉植物
    別讓摩擦成為學(xué)習(xí)的阻力
    報(bào)紙“黏土”等
    一起讀吧
    免费少妇av软件| 婷婷色综合www| 你懂的网址亚洲精品在线观看| 精品一品国产午夜福利视频| 精品亚洲乱码少妇综合久久| av片东京热男人的天堂| 亚洲美女黄色视频免费看| 男女免费视频国产| 一级片免费观看大全| 久久久久久久久久久久大奶| 亚洲熟女精品中文字幕| 成人手机av| 九草在线视频观看| 天堂中文最新版在线下载| 久久国产精品大桥未久av| 欧美日本中文国产一区发布| 国产精品av久久久久免费| 成年人免费黄色播放视频| 亚洲九九香蕉| 日韩av在线免费看完整版不卡| 18在线观看网站| 国产精品久久久久久人妻精品电影 | 久热爱精品视频在线9| 久久久久视频综合| 亚洲精品日韩在线中文字幕| 国产精品久久久人人做人人爽| 日韩制服丝袜自拍偷拍| 9191精品国产免费久久| 免费不卡黄色视频| 狂野欧美激情性xxxx| 女人精品久久久久毛片| 日韩av免费高清视频| 欧美大码av| 高清不卡的av网站| 久久精品国产a三级三级三级| 亚洲精品美女久久av网站| 欧美日韩视频精品一区| 看免费av毛片| 免费看av在线观看网站| 91精品三级在线观看| 嫁个100分男人电影在线观看 | 亚洲av美国av| 亚洲精品久久久久久婷婷小说| 两人在一起打扑克的视频| 欧美大码av| 国产成人av激情在线播放| 日本91视频免费播放| 一区福利在线观看| 蜜桃在线观看..| 男人舔女人的私密视频| 赤兔流量卡办理| av欧美777| 国产一区有黄有色的免费视频| 嫩草影视91久久| 国产成人系列免费观看| 久久国产精品男人的天堂亚洲| 成年人免费黄色播放视频| 免费高清在线观看视频在线观看| av片东京热男人的天堂| 黄色a级毛片大全视频| 手机成人av网站| 精品人妻1区二区| 国产精品久久久人人做人人爽| 中文字幕人妻丝袜制服| 精品国产乱码久久久久久男人| 最黄视频免费看| 亚洲欧洲日产国产| 国产免费一区二区三区四区乱码| 日日夜夜操网爽| 91九色精品人成在线观看| 黄色视频在线播放观看不卡| 亚洲精品日本国产第一区| 97人妻天天添夜夜摸| 热99国产精品久久久久久7| 99久久99久久久精品蜜桃| 最近手机中文字幕大全| 一级毛片我不卡| 欧美国产精品va在线观看不卡| 欧美国产精品va在线观看不卡| 亚洲欧美日韩另类电影网站| 成人国产一区最新在线观看 | 乱人伦中国视频| 国产97色在线日韩免费| 视频在线观看一区二区三区| 国产精品久久久久久人妻精品电影 | 国产亚洲一区二区精品| 老司机影院成人| 婷婷丁香在线五月| 国产成人欧美| 美女国产高潮福利片在线看| 曰老女人黄片| 丝袜在线中文字幕| 自线自在国产av| 国产亚洲av片在线观看秒播厂| 欧美精品高潮呻吟av久久| 黄色视频不卡| 精品人妻熟女毛片av久久网站| 欧美日本中文国产一区发布| 久久女婷五月综合色啪小说| 午夜福利,免费看| 免费人妻精品一区二区三区视频| 男的添女的下面高潮视频| 永久免费av网站大全| 成人国产av品久久久| 99久久综合免费| 国产片内射在线| 91麻豆av在线| 亚洲精品国产色婷婷电影| 99国产精品99久久久久| 97精品久久久久久久久久精品| 久久久久久久国产电影| 亚洲成av片中文字幕在线观看| 久9热在线精品视频| 久久 成人 亚洲| 黑人巨大精品欧美一区二区蜜桃| 狠狠婷婷综合久久久久久88av| 国产成人一区二区三区免费视频网站 | 亚洲国产精品成人久久小说| 国产一级毛片在线| 日韩 欧美 亚洲 中文字幕| 国产精品二区激情视频| 欧美日韩亚洲综合一区二区三区_| 亚洲国产欧美一区二区综合| 久久久久久免费高清国产稀缺| 国产精品99久久99久久久不卡| 国产亚洲精品第一综合不卡| 久久天躁狠狠躁夜夜2o2o | 人人妻,人人澡人人爽秒播 | 一区二区三区乱码不卡18| 在线观看一区二区三区激情| 丝袜人妻中文字幕| 亚洲成人手机| 亚洲成人免费av在线播放| 国产精品av久久久久免费| 日韩av免费高清视频| 亚洲三区欧美一区| 女性生殖器流出的白浆| 美女扒开内裤让男人捅视频| 悠悠久久av| 久久女婷五月综合色啪小说| 男女高潮啪啪啪动态图| 亚洲 国产 在线| 叶爱在线成人免费视频播放| 丰满迷人的少妇在线观看| 晚上一个人看的免费电影| 免费高清在线观看日韩| 激情五月婷婷亚洲| 亚洲av在线观看美女高潮| 狠狠婷婷综合久久久久久88av| 国产精品欧美亚洲77777| 国产精品一二三区在线看| 女性生殖器流出的白浆| 国产男人的电影天堂91| 日韩一区二区三区影片| 欧美日韩综合久久久久久| 亚洲欧洲日产国产| 久久久久久人人人人人| 欧美xxⅹ黑人| 50天的宝宝边吃奶边哭怎么回事| 无遮挡黄片免费观看| 久久久亚洲精品成人影院| 欧美久久黑人一区二区| 99精国产麻豆久久婷婷| 热99久久久久精品小说推荐| 人妻一区二区av| 狠狠婷婷综合久久久久久88av| 爱豆传媒免费全集在线观看| 亚洲欧美中文字幕日韩二区| 欧美+亚洲+日韩+国产| 国产成人系列免费观看| av视频免费观看在线观看| 人体艺术视频欧美日本| 各种免费的搞黄视频| 亚洲精品一卡2卡三卡4卡5卡 | 极品少妇高潮喷水抽搐| 美女午夜性视频免费| 高清视频免费观看一区二区| 久久毛片免费看一区二区三区| 中文字幕人妻丝袜制服| 中文精品一卡2卡3卡4更新| 国产精品二区激情视频| 青青草视频在线视频观看| 下体分泌物呈黄色| 女人高潮潮喷娇喘18禁视频| av在线app专区| 精品亚洲成a人片在线观看| 精品人妻1区二区| 亚洲欧美清纯卡通| 免费在线观看视频国产中文字幕亚洲 | 国产亚洲一区二区精品| 免费观看av网站的网址| 天堂8中文在线网| 国产亚洲av片在线观看秒播厂| 天天添夜夜摸| 成人18禁高潮啪啪吃奶动态图| 久久精品国产a三级三级三级| 电影成人av| 一区二区三区四区激情视频| av视频免费观看在线观看| av福利片在线| 国产精品免费视频内射| 叶爱在线成人免费视频播放| 亚洲av电影在线观看一区二区三区| 亚洲精品一区蜜桃| 男女边摸边吃奶| 一级毛片女人18水好多 | 精品一区在线观看国产| 国产高清视频在线播放一区 | 免费不卡黄色视频| 欧美成人精品欧美一级黄| 肉色欧美久久久久久久蜜桃| 精品一区二区三区av网在线观看 | 丝袜美足系列| 亚洲第一青青草原| 青草久久国产| 日本欧美国产在线视频| 久久天堂一区二区三区四区| 国产一区亚洲一区在线观看| 女人被躁到高潮嗷嗷叫费观| 中文字幕另类日韩欧美亚洲嫩草| 男女高潮啪啪啪动态图| 亚洲精品一区蜜桃| 亚洲欧美精品自产自拍| 国产精品免费大片| 妹子高潮喷水视频| 亚洲国产成人一精品久久久| 黄色毛片三级朝国网站| 日本av手机在线免费观看| 精品少妇一区二区三区视频日本电影| 国产xxxxx性猛交| a级片在线免费高清观看视频| 香蕉国产在线看| 极品人妻少妇av视频| 国产一区二区 视频在线| 亚洲成人手机| 啦啦啦啦在线视频资源| 岛国毛片在线播放| 久久免费观看电影| 久久ye,这里只有精品| 午夜激情久久久久久久| 亚洲国产欧美日韩在线播放| 久久人妻福利社区极品人妻图片 | 精品亚洲乱码少妇综合久久| 日本a在线网址| 搡老乐熟女国产| 国产视频一区二区在线看| 99国产精品免费福利视频| 可以免费在线观看a视频的电影网站| 亚洲成人免费电影在线观看 | 亚洲欧美日韩另类电影网站| 大型av网站在线播放| 高清不卡的av网站| 99国产精品免费福利视频| 熟女少妇亚洲综合色aaa.| 久久人人97超碰香蕉20202| 男女边摸边吃奶| 国产在视频线精品| 成人免费观看视频高清| 久久久久久久久免费视频了| 亚洲第一av免费看| 色婷婷av一区二区三区视频| 久热爱精品视频在线9| 欧美激情高清一区二区三区| 91精品伊人久久大香线蕉| 欧美老熟妇乱子伦牲交| www.精华液| 97在线人人人人妻| 精品国产一区二区久久| 国产成人欧美| 91国产中文字幕| 在线 av 中文字幕| 首页视频小说图片口味搜索 | 又大又黄又爽视频免费| 大香蕉久久成人网| 青春草亚洲视频在线观看| 女人被躁到高潮嗷嗷叫费观| 高清黄色对白视频在线免费看| 久久这里只有精品19| 脱女人内裤的视频| 悠悠久久av| 成人亚洲精品一区在线观看| 精品熟女少妇八av免费久了| 啦啦啦视频在线资源免费观看| 在线亚洲精品国产二区图片欧美| 欧美老熟妇乱子伦牲交| 新久久久久国产一级毛片| 国产一区二区激情短视频 | 在线观看免费午夜福利视频| 一级毛片我不卡| 日韩精品免费视频一区二区三区| 久久九九热精品免费| 久久国产精品人妻蜜桃| 亚洲图色成人| 视频在线观看一区二区三区| 国精品久久久久久国模美| av欧美777| 少妇人妻久久综合中文| 一本久久精品| 免费看不卡的av| 国产精品 国内视频| 日本wwww免费看| 国产精品99久久99久久久不卡| av有码第一页| 国产亚洲午夜精品一区二区久久| 国产视频一区二区在线看| 欧美激情 高清一区二区三区| 欧美日韩一级在线毛片| 91九色精品人成在线观看| 午夜福利一区二区在线看| 波多野结衣一区麻豆| 王馨瑶露胸无遮挡在线观看| 涩涩av久久男人的天堂| 欧美激情高清一区二区三区| 人人澡人人妻人| 精品视频人人做人人爽| 大香蕉久久成人网| 色播在线永久视频| 妹子高潮喷水视频| 人成视频在线观看免费观看| 国产一区二区 视频在线| 免费少妇av软件| 久久99一区二区三区| 伊人亚洲综合成人网| 国产成人av激情在线播放| 亚洲av成人不卡在线观看播放网 | 亚洲国产看品久久| 大香蕉久久成人网| 亚洲国产最新在线播放| 80岁老熟妇乱子伦牲交| 午夜老司机福利片| 2021少妇久久久久久久久久久| 高清视频免费观看一区二区| 大片电影免费在线观看免费| 欧美日韩av久久| 女人久久www免费人成看片| 精品久久蜜臀av无| 黄色 视频免费看| 精品免费久久久久久久清纯 | 黄色怎么调成土黄色| 久久人人97超碰香蕉20202| 在线看a的网站| 一级毛片女人18水好多 | 少妇的丰满在线观看| 免费一级毛片在线播放高清视频 | 极品人妻少妇av视频| 久久热在线av| 韩国精品一区二区三区| 波野结衣二区三区在线| 亚洲欧洲日产国产| 好男人电影高清在线观看| 午夜福利免费观看在线| 久久久国产一区二区| 欧美亚洲日本最大视频资源| 久久人妻福利社区极品人妻图片 | 久久国产精品大桥未久av| 在线观看免费高清a一片| 国产一区二区三区综合在线观看| 啦啦啦中文免费视频观看日本| 99国产精品免费福利视频| 精品一区二区三区av网在线观看 | 久久天堂一区二区三区四区| 亚洲国产日韩一区二区| 中文欧美无线码| 高清不卡的av网站| 国产精品国产av在线观看| 成人黄色视频免费在线看| 丁香六月欧美| 日韩视频在线欧美| 国产精品一国产av| 2021少妇久久久久久久久久久| 国产精品二区激情视频| xxx大片免费视频| 老司机靠b影院| 老鸭窝网址在线观看| 美女福利国产在线| 色94色欧美一区二区| 亚洲成人手机| 肉色欧美久久久久久久蜜桃| 在线 av 中文字幕| 久久久亚洲精品成人影院| 久久久久国产精品人妻一区二区| 母亲3免费完整高清在线观看| 黄色视频在线播放观看不卡| 麻豆av在线久日| 狂野欧美激情性bbbbbb| 国产精品国产三级专区第一集| 少妇 在线观看| 欧美亚洲 丝袜 人妻 在线| 一级黄片播放器| 日韩av在线免费看完整版不卡| 9热在线视频观看99| 久久久亚洲精品成人影院| 精品少妇一区二区三区视频日本电影| 丝袜人妻中文字幕| 亚洲精品在线美女| 国产欧美日韩综合在线一区二区| 亚洲少妇的诱惑av| 各种免费的搞黄视频| 精品一区二区三卡| 欧美精品av麻豆av| 日本黄色日本黄色录像| 精品国产国语对白av| 老司机影院毛片| 欧美日韩国产mv在线观看视频| 精品国产一区二区三区久久久樱花| 亚洲精品国产av成人精品| 日本a在线网址| 一边摸一边抽搐一进一出视频| 日韩人妻精品一区2区三区| 国产亚洲av高清不卡| 免费高清在线观看视频在线观看| 免费看av在线观看网站| 国产成人av教育| 人人妻,人人澡人人爽秒播 | 亚洲中文字幕日韩| 日本av免费视频播放| 晚上一个人看的免费电影| 又大又黄又爽视频免费| 91成人精品电影| 国产精品亚洲av一区麻豆| 国产一区有黄有色的免费视频| 日韩中文字幕视频在线看片| 成年动漫av网址| 亚洲欧洲精品一区二区精品久久久| 国产不卡av网站在线观看| 两个人看的免费小视频| 女人久久www免费人成看片| 99九九在线精品视频| 久久中文字幕一级| 精品福利观看| 在线观看免费日韩欧美大片| 婷婷成人精品国产| 免费在线观看黄色视频的| 久久人人97超碰香蕉20202| 桃花免费在线播放| 成人午夜精彩视频在线观看| 免费在线观看日本一区| 久久女婷五月综合色啪小说| 欧美精品啪啪一区二区三区 | 国产极品粉嫩免费观看在线| 新久久久久国产一级毛片| 又粗又硬又长又爽又黄的视频| 欧美黑人精品巨大| 男人舔女人的私密视频| 亚洲久久久国产精品| 欧美成人精品欧美一级黄| 免费在线观看视频国产中文字幕亚洲 | 久久精品国产综合久久久| 久久毛片免费看一区二区三区| 日本vs欧美在线观看视频| 亚洲精品在线美女| 国产女主播在线喷水免费视频网站| 国产成人a∨麻豆精品| 国产又色又爽无遮挡免| 成人三级做爰电影| 精品人妻一区二区三区麻豆| 国产精品九九99| 国产福利在线免费观看视频| 国产成人免费无遮挡视频| 免费看不卡的av| 欧美精品av麻豆av| avwww免费| 亚洲欧美成人综合另类久久久| 激情视频va一区二区三区| 亚洲一码二码三码区别大吗| 国产99久久九九免费精品| 两个人免费观看高清视频| 性少妇av在线| 美女视频免费永久观看网站| 最新的欧美精品一区二区| 老司机影院成人| 国产伦人伦偷精品视频| 侵犯人妻中文字幕一二三四区| 成年人黄色毛片网站| 国产在线一区二区三区精| 婷婷丁香在线五月| 深夜精品福利| 欧美乱码精品一区二区三区| 久久久亚洲精品成人影院| 亚洲国产欧美日韩在线播放| 桃花免费在线播放| 国产又色又爽无遮挡免| 中文字幕高清在线视频| 国产成人精品久久二区二区91| 黄网站色视频无遮挡免费观看| 亚洲欧美成人综合另类久久久| 国产xxxxx性猛交| 中国国产av一级| 色94色欧美一区二区| a级毛片黄视频| 国产色视频综合| 久久人人97超碰香蕉20202| 亚洲综合色网址| 亚洲av美国av| av电影中文网址| 色精品久久人妻99蜜桃| 老司机靠b影院| 国产成人91sexporn| 国产精品久久久av美女十八| 欧美久久黑人一区二区| 久久人妻福利社区极品人妻图片 | 人妻一区二区av| 国产av国产精品国产| 亚洲欧美精品自产自拍| 妹子高潮喷水视频| 一本色道久久久久久精品综合| 久久精品熟女亚洲av麻豆精品| 校园人妻丝袜中文字幕| 亚洲九九香蕉| 一二三四在线观看免费中文在| 免费在线观看黄色视频的| 美女午夜性视频免费| 欧美久久黑人一区二区| 欧美激情高清一区二区三区| 久久久精品区二区三区| 国产精品亚洲av一区麻豆| 欧美老熟妇乱子伦牲交| 人妻人人澡人人爽人人| 18在线观看网站| 精品亚洲乱码少妇综合久久| 国产视频首页在线观看| 大陆偷拍与自拍| 国产视频一区二区在线看| 亚洲伊人久久精品综合| 少妇精品久久久久久久| 伦理电影免费视频| 人人妻人人澡人人爽人人夜夜| 曰老女人黄片| 男女边摸边吃奶| av线在线观看网站| 中文字幕最新亚洲高清| 亚洲九九香蕉| 妹子高潮喷水视频| av又黄又爽大尺度在线免费看| 在线精品无人区一区二区三| 我要看黄色一级片免费的| 成人手机av| 亚洲国产毛片av蜜桃av| 久久久久久久久免费视频了| 男女之事视频高清在线观看 | 日韩一本色道免费dvd| www.自偷自拍.com| 久久ye,这里只有精品| 大陆偷拍与自拍| 亚洲成人免费电影在线观看 | 18禁国产床啪视频网站| 亚洲欧美一区二区三区久久| 久久精品亚洲av国产电影网| 国产精品久久久av美女十八| 啦啦啦中文免费视频观看日本| 男女下面插进去视频免费观看| 欧美日韩亚洲高清精品| 亚洲综合色网址| 欧美激情极品国产一区二区三区| 777久久人妻少妇嫩草av网站| 精品国产乱码久久久久久小说| 一区二区日韩欧美中文字幕| 王馨瑶露胸无遮挡在线观看| 婷婷色麻豆天堂久久| 亚洲精品一区蜜桃| 9色porny在线观看| 在线看a的网站| 999久久久国产精品视频| 两人在一起打扑克的视频| 国产在线观看jvid| 成人三级做爰电影| www.自偷自拍.com| 日韩制服骚丝袜av| 一个人免费看片子| 国产精品人妻久久久影院| 免费黄频网站在线观看国产| 亚洲免费av在线视频| 无遮挡黄片免费观看| 欧美精品人与动牲交sv欧美| 老司机影院成人| 又大又爽又粗| 悠悠久久av| 两人在一起打扑克的视频| 午夜免费观看性视频| av国产久精品久网站免费入址| 91国产中文字幕| 精品一区二区三卡| 男男h啪啪无遮挡| 超色免费av| 少妇精品久久久久久久| 一区二区三区四区激情视频| 日韩av在线免费看完整版不卡| 新久久久久国产一级毛片| 国产精品三级大全| 亚洲av电影在线进入| 精品国产超薄肉色丝袜足j| av又黄又爽大尺度在线免费看| 国产三级黄色录像| 国产片特级美女逼逼视频| 亚洲精品日韩在线中文字幕| 老司机靠b影院| 免费人妻精品一区二区三区视频| 无限看片的www在线观看| 欧美xxⅹ黑人| 在线观看www视频免费| 免费观看a级毛片全部| 国产深夜福利视频在线观看| 国产成人欧美在线观看 | 大片电影免费在线观看免费| 国产在线一区二区三区精| 亚洲人成电影免费在线| 亚洲第一av免费看| 黄色a级毛片大全视频| 中文字幕另类日韩欧美亚洲嫩草| 国产日韩欧美视频二区| 少妇人妻久久综合中文| 精品一品国产午夜福利视频| 久久人人爽av亚洲精品天堂| 亚洲三区欧美一区| 99精品久久久久人妻精品|