亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        軟土三軸剪切蠕變試驗研究及模型分析

        2012-11-05 07:25:22孟慶山
        巖土力學(xué) 2012年1期
        關(guān)鍵詞:數(shù)字圖像測量方法徑向

        楊 超,汪 稔,孟慶山

        (中國科學(xué)院武漢巖土力學(xué)研究所 巖土力學(xué)與工程國家重點實驗室,武漢 430071)

        1 引 言

        在自重或坡頂荷載作用下可能會發(fā)生基坑邊坡滑移,導(dǎo)致新開挖基坑工程附近的建筑開裂或傾斜,嚴重危害基坑周圍建筑物的安全。上海是典型的濱海相沉積軟土地區(qū),21世紀以來上海地區(qū)的基坑工程呈現(xiàn)出一些新的特點:基坑規(guī)模越來越大;開挖深度越來越深;場地臨近建筑紅線范圍,周圍含有大量管線等其他構(gòu)筑物。工程實踐表明,上海地區(qū)的基坑和建筑物地基等在長期荷載作用下變形具有明顯的時間效應(yīng),即流變現(xiàn)象,這種現(xiàn)象嚴重危害基坑及基坑周圍建筑物和構(gòu)筑物的安全使用,因此,研究上海軟土的剪切流變特性對上海地區(qū)工程建設(shè)有重要的現(xiàn)實意義。

        軟土剪切流變試驗常采用三軸流變儀進行,由于測量技術(shù)上的限制,人們難以得到試樣徑向變形的準確值,直接影響對土的剪切流變特性的認識。通常徑向變形通過試樣的體積變化換算獲得,試樣徑向變形的不均勻使測量結(jié)果存在一定誤差,影響剪切應(yīng)變的計算結(jié)果,且在不排水或非飽和三軸試驗中更加難以獲得試樣的徑向變形,無法獲得試樣的剪切應(yīng)變。研究人員已經(jīng)對傳統(tǒng)三軸儀的徑向量測方法進行了很多改進[1],但這些測量方法受到操作難度或者測量精度的影響,多為接觸式測量,會對試樣產(chǎn)生約束,影響試驗結(jié)果。

        數(shù)字圖像測量技術(shù)能夠有效地彌補這些測量方法中的不足,在不接觸試樣的基礎(chǔ)上實現(xiàn)試樣變形的精確測量。這種方法是通過對圖像的采集、處理分析和識別等一系列過程,量測圖像中試樣大小,實現(xiàn)對試樣變形的測量。數(shù)字圖像測量技術(shù)在土工三軸試驗試樣徑向變形量測中已經(jīng)取得了很好的效果[1-5],實現(xiàn)了對三軸試樣徑向變形的直接測量,具有較高的測量精度。

        2 數(shù)字圖像測測量方法及誤差分析

        三軸蠕變試驗的測量精度直接關(guān)系到試驗結(jié)果的可靠性和準確性,因此,試驗過程中應(yīng)采取一定措施消除或減少測量系統(tǒng)帶來的系統(tǒng)誤差。

        2.1 數(shù)字圖像測量原理

        數(shù)字圖像通常表示為M×N個像素,在 Matlab中存儲為一個多維數(shù)據(jù)矩陣。三軸試驗數(shù)字圖像測量時采用 unit8型數(shù)組的灰度圖像來描述,像素灰度值用 0~255表示黑~白之間灰度變化[6]。在灰度圖像下,物體邊界表現(xiàn)為灰度值的不連續(xù)性,通過灰度函數(shù)的一階或二階導(dǎo)數(shù)變化即可判斷出物體的邊界位置。在實際圖像處理過程中,常采用直方圖或空域濾波器等手段增強圖像的對比效果,以便在圖像處理時準確判別試樣邊界。

        2.2 三軸蠕變試驗數(shù)字圖像測量方法

        利用數(shù)字圖像測量方法進行三軸蠕變試驗的試驗裝置主要由三軸蠕變試驗系統(tǒng)、圖像采集系統(tǒng)和數(shù)字圖像處理系統(tǒng)三部分組成。三軸蠕變系統(tǒng)為應(yīng)力控制式三軸儀,采用分級加載的方式進行;用數(shù)碼相機進行三軸試樣變形的圖像采集,選擇合適的位置固定相機三角支架,保證在整個試驗過程中相機的位置不變,在試驗開始之前進行相機對焦,整個試驗過程中保持焦距不變,且鏡頭不發(fā)生偏移;按照試驗的設(shè)計時長進行圖像采集,并傳入計算機,再應(yīng)用 Matlab強大的數(shù)字圖像處理功能對圖像進行數(shù)字化處理分析,得到三軸蠕變試驗中試樣的變形信息。

        2.3 系統(tǒng)標定

        數(shù)字圖像中每個像素點都表示一定的真實長度,通過對像素點的統(tǒng)計分析就可以準確獲得物體的長度值,從而得到物體變形的真實值。

        對測量系統(tǒng)標定可以獲得每個像素代表長度的真實值。在距離物體 D處安放好圖像采集儀后,計算長度為L0的物體所占的像素值N0,即可以得出每個像素在距離采集儀 D處代表的實際長度l:

        測量時,可根據(jù)物體在圖像中所占像素值N換算得到該物體實際的長度:L=lN。在標定和測量過程中保持采集儀和物體之間的距離D不變。

        數(shù)字圖像測量中,單位像素代表的真實長度越小,測量中產(chǎn)生的誤差就越小,即l越小,測量精度越高。當(dāng)測量精度要求很高時應(yīng)優(yōu)先選用像素較高的圖像采集系統(tǒng)以保證試驗結(jié)果的可靠性。

        2.4 成像原理及誤差分析

        進行數(shù)字圖像測量試樣徑向變形,實際測得的是弦AB的長度而不是直徑2R的長度,因此,需要通過換算才能得到試樣直徑的長度2R。

        由成像原理可知,成像時試驗不同階段中試樣的弦AB、A′B′,鏡頭距試樣中心的距離為D,其中 A、B、A′、B′分別為變形前和變形后圖像物體的邊緣(見圖1),由三角函數(shù)關(guān)系有

        由式(2)可以看出,當(dāng)∠AOM 越接近π/2時,即試樣半徑R與鏡頭距試樣中心的距離D比值R/D越小,弦AB與直徑2R越接近。三軸試驗試樣直徑常采用39.1 mm,當(dāng)鏡頭距離試樣中心1.0、1.5、2 m時,誤差分別為0.0064、0.0032、0.0018 mm,相對誤差僅為0.19%、0.08%、0.05%,此時測量誤差與真值相比已達10-3量級,誤差相對試樣真實直徑而言是一個非常小的值。因此,當(dāng)鏡頭距離試樣足夠遠時,可以用弦長AB近似代替試樣直徑2R。同樣,當(dāng)距離D一定時,半徑R越小,弦長AB越接近直徑2R長度,誤差也越小。

        圖1 物體成像示意圖Fig.1 Sketch of imaging

        系統(tǒng)標定時,若采用較小直徑圓柱進行標定,誤差會隨著試驗中試樣直徑的增大而增大。因此,采用直徑較大的圓柱徑向進行標定,能夠有效地降低在試驗中試樣鼓脹而產(chǎn)生的誤差,保證試驗質(zhì)量。

        3 三軸蠕變試驗的數(shù)字圖像測量成果

        本次試驗試樣采用上海五號溝地區(qū)15 m深度處灰褐色淤泥質(zhì)黏土進行三軸剪切蠕變試驗,土的物理力學(xué)性質(zhì)見表 1。上海五號溝地區(qū)臨近長江口,與上海長江隧道工程相鄰。長江口位于長江三角洲沖積平原的東南前緣,成陸較晚。60 m深度以上主要為飽和的黏性土和粉性土,具有含水率高、抗剪強度低、滲透性小、固結(jié)時間長、流變特性顯著等特點。

        試驗采用分級加載方法,考慮所取土層深度的實際情況,選用圍壓σ3= 100 kPa進行三軸不排水和排水剪切蠕變試驗和圍壓σ3= 200 kPa進行三軸排水剪切蠕變試驗。試驗共分5級進行加載,加載等級按照試樣在相同條件下的三軸試驗強度στ的1/5進行加載。

        表1 土的物理力學(xué)性質(zhì)參數(shù)Table1 Physico-mechanical properties of soil

        3.1 三軸蠕變試驗成果

        圖2~4為三軸剪切蠕變試驗成果曲線(圖中CU為三軸排水剪切試驗,CD為三軸不排水剪切試驗)??傮w上看,較低偏應(yīng)力水平D(偏應(yīng)力水平指試樣上施加的偏應(yīng)力σ與破壞時最大偏應(yīng)力στ之比,即D=σ/στ,D=0~1)下,試樣變形穩(wěn)定較快,表現(xiàn)為衰減蠕變。隨著應(yīng)力水平D的提高,逐漸表現(xiàn)為先穩(wěn)定蠕變一段時間后發(fā)生衰減。當(dāng)試樣受到的D較高時,試樣蠕變并未表現(xiàn)出明顯的衰減現(xiàn)象,在很長時間內(nèi)變形仍處于增長狀態(tài),此次加載歷時14 d,試樣變形還在緩慢增長。當(dāng)試樣偏應(yīng)力σ達到三軸剪切試驗強度στ時,試樣在較短時間內(nèi)產(chǎn)生大量變形,隨后形成剪切面,剪切變形迅速增加發(fā)生剪切破壞。

        不排水時,試驗在前3級偏應(yīng)力作用下,變形較小,土的結(jié)構(gòu)較完整,偏應(yīng)力主要由土體結(jié)構(gòu)承受。當(dāng)偏應(yīng)力到達第四級時,土體結(jié)構(gòu)已經(jīng)被破壞,試樣在加載后較短時間內(nèi)產(chǎn)生大量變形且后期蠕變量也較大(見圖2)。在排水條件下,剪切蠕變試驗有明顯的固結(jié)硬化現(xiàn)象,由于排水作用的影響,試樣發(fā)生固結(jié)硬化強度提高,各級偏應(yīng)力水平下產(chǎn)生的變形差異不如不排水條件顯著(見圖3)。

        加載時,孔隙水壓力在短時間內(nèi)完成上升和下降過程(如圖 2(b)~4(b)孔隙水壓力與時間關(guān)系曲線)。不排水條件下,試樣在受到豎向壓力作用時,試樣孔隙被壓縮,由于無法排水,試樣無法固結(jié),孔隙水壓力只能在試樣內(nèi)部孔隙之間實現(xiàn)平衡。加載過程中,孔隙水壓力并無下降的趨勢,經(jīng)過短時間的平衡后,隨著試樣變形和孔隙的壓縮,孔隙水壓力呈緩慢上升趨勢。由Terzaghi和Biot固結(jié)原理可知,排水時,加載瞬間孔隙水壓力急劇上升,試樣中的水在高壓力下迅速排出產(chǎn)生固結(jié)。由于試樣僅從試樣頂部的透水石處排水,試樣頂部土體很快發(fā)生固結(jié),于是產(chǎn)生 Manadei-Cryer效應(yīng),孔隙水壓力在一段時間內(nèi)依然處于上升階段,持續(xù)時間約為10 min。同時,在進行三軸排水剪切蠕變試驗時排水路徑受到限制,試樣上部率先發(fā)生固結(jié)硬化,下部固結(jié)較慢,導(dǎo)致試樣上部鼓脹通常較下端小,剪切面常常在試樣下部形成。

        土粒周圍水膜的厚度決定著土體黏滯系數(shù)的大小,對軟土流變起著十分重要的作用[7]。不排水試驗中,土粒周圍的水膜未發(fā)生明顯的變化,試樣的黏滯系數(shù)基本為一個定值,蠕變主要是由于土顆粒和土骨架的移動或變形產(chǎn)生的。排水時,受到豎向壓力作用時,孔隙被壓縮,試樣中水被排出,土顆粒相互靠近,使得土粒之間水膜變薄,主要表現(xiàn)為土的密度增大,黏滯性增強,變形速率變小。隨著試樣中水進一步排出,土顆粒之間水膜越來越薄,土的黏滯性也越來越強。在試樣加載初期,孔隙水壓力迅速增大使得試樣中水排擠較快,試樣中土粒水膜厚度發(fā)生較大變化,土的固結(jié)和黏滯性變化明顯。隨著時間的推移,試樣中排水速度逐漸減小,此時,土顆粒間的水膜厚度變化較小,相對穩(wěn)定,試樣中的黏滯系數(shù)變化緩慢。對照蠕變試驗成果圖可以看出,試驗加載初期試樣中厚度孔隙水壓力和排水量較大,水膜變化也較大,與之對應(yīng)的試樣變形也較大為穩(wěn)定蠕變;隨后試樣的孔隙水壓力和排水量逐漸穩(wěn)定,此時土顆粒間水膜厚度也變化不大,使得試樣的黏滯系數(shù)基本為一個穩(wěn)定的數(shù)值,試樣變形逐漸穩(wěn)定,即衰減蠕變,可見不排水條件下,隨著加載等級的提高,試樣的蠕變曲線的斜率較排水條件下大。

        排水時,試樣蠕變和固結(jié)相互影響。一方面,固結(jié)作用會使得土體發(fā)生硬化,從而減小試樣的變形。另一方面,試樣蠕變也會阻塞排水通道,影響土體的滲透性。試驗中,固結(jié)和蠕變二者同時進行,并在一段時間內(nèi)達到穩(wěn)定平衡,整個過程中試樣孔壓消散與試樣固結(jié)變形時間基本一致。

        對比不同圍壓下排水時的蠕變曲線(見圖3、4)可以發(fā)現(xiàn),偏應(yīng)力σ近乎相等時,圍壓為 100 kPa時試樣變形量較圍壓為200 kPa時大,表明周圍壓力對土樣的蠕變有較大的影響,但在相同偏應(yīng)力水平 D下圍壓為 200 kPa時試樣的變形量較圍壓為100 kPa時大。

        3.2 圖像測量法與排水體積法的比較

        圖5為采用圖像測量方法和排水體積法測量得到的剪切蠕變曲線。剪切應(yīng)變-時間對數(shù)曲線與垂直應(yīng)變-時間曲線變化相似,是由于試樣的剪切應(yīng)變是通過軸向應(yīng)變和徑向應(yīng)變計算得到,而徑向的變形相對軸向小很多,但通過兩種不同測量方法獲得的剪切應(yīng)變有一定的差別。在偏應(yīng)力水平較低時,圖像測量方法與排水體積方法得到的剪切應(yīng)變值相差不大,曲線基本重合;隨著偏應(yīng)力水平D增大,二者差值逐步增大,圖像測量得到的剪切應(yīng)變值大于排水體積方法得到的剪切應(yīng)變值。這一現(xiàn)象符合實際情況,在三軸試驗過程中,試樣變形是不均勻的,此次圖像測量方法得到的是試樣1/2高度截面的變形值,而排水體積方法得到的是整個試樣的平均變形。在偏應(yīng)力水平較低時,試樣徑向變形較小,試樣1/2高度截面處的變形值與整個試樣的平均變形相差不大;隨著偏應(yīng)力水平的增大,試樣發(fā)生鼓脹變形,試樣1/2高度截面的變形值明顯比試樣整體變形的平均值大,隨著變形發(fā)展這個差值逐漸增大(見圖 6)。采用試樣排水體積的方法進行測量獲得的剪切應(yīng)變值也隨試樣變形的發(fā)展誤差在不斷擴大。

        圖2 σ3=100 kPa時CU三軸蠕變試驗成果Fig.2 Results of CU triaxial creep test (σ3=100 kPa)

        圖3 σ3=100 kPa時CD三軸蠕變試驗成果Fig.3 Results of CD triaxial creep test (σ3=100 kPa)

        圖4 σ3=200 kPa時CD三軸蠕變試驗成果Fig.4 Results of CD triaxial creep test (σ3=200 kPa)

        圖5 不同測量方法的剪切應(yīng)變與時間關(guān)系曲線Fig.5 Relationships of shear strain and time in different measuring methods

        圖6 不同偏應(yīng)力下試樣的二值圖Fig.6 Binary images in different deviatoric stresses

        3.3 圖像測量與排水體積法得到的Singh-Mitchell模型參數(shù)比較

        Singh-Mitchell提出能夠恰當(dāng)描述多種土在20%~80%范圍內(nèi)偏應(yīng)力水平的應(yīng)變速率-時間關(guān)系特征方程[8-10]:

        圖7為 3種不同情況下三軸剪切蠕變試驗的lnε-lnt關(guān)系曲線和lnε-D關(guān)系曲線。由圖可以看出,不排水條件下各級偏應(yīng)力作用下剪切應(yīng)變和時間在雙對數(shù)坐標下線性關(guān)系良好,且斜率基本一致。排水條件下,除第一級偏應(yīng)力下直線斜率差別較大以為外,其他各級偏應(yīng)力水平下lnε-lnt曲線基本平行,且圖像測量方法和排水體積方法獲得的lnε-lnt曲線基本重合。圖中選用t1=1 d時lnε-D曲線進行 Singh-Mitchell模型參數(shù)的確定,由圖可以看出,lnε和D之間線性關(guān)系明顯,線性相關(guān)系數(shù)分別為0.9867、0.9740和0.9947。因此,采用具有廣泛應(yīng)用的 Singh-Mitchell模型是適合的,其中λ為lnε-lnt的斜率。各參數(shù)見表3。

        通過兩種不同的測量方法獲得的 Singh-Mitchell模型參數(shù)的差別不大(見表2),應(yīng)力水平較低且難以準確測量徑向變形時可近似用排水體積法進行參數(shù)確定。同時,Singh-Mitchell模型很好的描述上海軟土在偏應(yīng)力水平在20%~80%范圍內(nèi)的流變特性,對于工程應(yīng)用基本可以滿足要求。排水時,對于偏應(yīng)力水平大于80%時上海軟土剪切流變特性也可以較好的描述,而對于偏應(yīng)力水平低于20%時軟土的剪切流變特性的描述差別較大。

        圖7 不同條件下lnε-lnt曲線和lnε-D曲線Fig.7 Curves of lnε-lnt and lnε-D in different test conditions

        表2 Singh-Mitchell模型參數(shù)Table2 Parameters of Singh-Mitchell model

        4 結(jié) 論

        (1)數(shù)字圖像測量技術(shù)簡單實用,且測量精度較高,能夠滿足高精度土工試驗測量的需要。在三軸蠕變試驗中能夠有效的測得試樣的徑向變形,實現(xiàn)試樣的非接觸式測量,減少接觸式測量對試驗結(jié)果的影響,為難以測量試樣體積變化的不排水試驗和非飽和試驗提供了有效的徑向變形測量途徑。

        (2)進行三軸蠕變試驗圖像測量時,必須嚴格保證圖像采集系統(tǒng)的鏡頭在標定和試驗過程中不發(fā)生偏移或不產(chǎn)生相對位置變化,避免由于鏡頭偏移而產(chǎn)生測量誤差;試樣鼓脹會使測量誤差增大,采用較大直徑的圓柱進行圖像測量標定能夠有效的抑制這一現(xiàn)象的產(chǎn)生,提高測量精度。

        (3)在三軸排水剪切蠕變試驗中,在短時間內(nèi)孔隙水壓力會持續(xù)上升產(chǎn)生 Manadei-Cryer效應(yīng),整個過程歷時約10 min。試驗過程中,排水路徑的影響使得試樣上部率先發(fā)生固結(jié)硬化,試樣鼓脹和剪切面通常偏向試樣下部產(chǎn)生。

        (4)軟土的流變性質(zhì)與土中水的含量密切相關(guān)。土中含水率高時,土顆粒周圍水膜較厚,黏滯系數(shù)較低,土顆粒移動相對較容易;反之,水膜較薄,黏滯系數(shù)較大。排水剪切流變試驗時,隨著隨著土體中水的排出土樣中水含量減少,顆粒周圍的水膜變薄,土的密度增大、黏滯性增強,土的流變性減弱。

        (5)由于試樣變形的不均勻,獲得土的剪應(yīng)變時圖像測量方法和排水體積測量方法存在一定的差異。圖像測量方法在獲得土的剪應(yīng)變時排水體積法要大,隨著荷載等級的增加、試樣鼓脹變形增大,這一差別更加明顯。在進行 Singh-Mitchell模型參數(shù)確定時,兩種方法獲得的參數(shù)差別不大,在徑向變形不大且難以準確測量時可用排水體積法代替。

        本次試驗僅從一個方向進行試樣變形的圖像采集,當(dāng)試樣產(chǎn)生上不均勻變化或形成剪切面時不能準確獲取試樣滑移的試驗信息,建議在試驗時至少從兩個方向獲得試驗變形圖像以便準確合理的獲取試樣變形信息。

        [1]邵龍?zhí)?王助貧,韓國城,等. 三軸試驗土樣徑向變形的計算機圖像測量[J]. 巖土工程學(xué)報,2001,23(3): 337-341.SHAO Long-tan,WANG Zhu-pin,HAN Guo-cheng,et al.Digital image processing technique for measurement of the radial deformation of specimen in triaxial test[J].Chinese Journal of Geotechnical Engineering,2001,23(3): 337-341.

        [2]Ali Hakan ?ren,Okan ?nal,Gürkan ?zden,Abidin Kaya.Nondestructive evaluation of volumetric shrinkage of compacted mixtures using digital image analysis[J].Engineering Geology,2006,85: 239-250.

        [3]劉永祿. 三軸實驗土樣變形數(shù)字圖像測量的實現(xiàn)[D].大連: 大連理工大學(xué),2002.

        [4]邵龍?zhí)?孫益振,王助貧,等. 數(shù)字圖像測量技術(shù)在土工三軸試驗中的應(yīng)用研究[J]. 巖土力學(xué),2006,27(1):29-34.SHAO Long-tan,SUN Yi-zhen,WANG Zhu-pin,et al.Application of digital image processing technique to triaxial test in soil mechanics[J]. Roch and Soil Mechanics,2006,27(1): 29-34.

        [5]董建軍,邵龍?zhí)?劉永祿,等. 基于圖像測量方法的非飽和壓實土三軸試樣變形測量[J]. 巖土力學(xué),2008,29(6): 1618-1622.DONG Jian-jun,SHAO Long-tan,LIU Yong-lu,et al.Measurement of deformation of unsaturated compacted soil triaxial specimen based on digital image measurement method[J]. Rock and Soil Mechanics,2008,29(6): 1618-1622.

        [6](美)岡薩雷斯(Gonzalez),等. 數(shù)字圖像處理(MATLAB版)[M]. 阮秋琦,等譯. 北京: 電子工業(yè)出版社,2005.

        [7]何俊,肖樹芳. 結(jié)合水對海積軟土流變性質(zhì)的影響[J].吉林大學(xué)學(xué)報,2003,33(2): 204-207.HE Jun,XIAO Shu-fang. Some influence of bond water on rheological properties of marine soft soils[J]. Journal of Changchun University of Science and Technology,2003,33(2): 204-207.

        [8]SINGH A,MITCHELL J K. General stress-strain-time function for soils[J]. Journal of Soil Mechanics and Foundation Division,ASCE,1968,94(1): 21-46.

        [9]李軍世,林詠梅. 上海淤泥質(zhì)粉質(zhì)黏土的 Singh-Mitchell蠕變模型[J]. 巖土力學(xué),2000,21(4): 363-366.LI Jun-shi,LIN Yong-mei. Singh-Mitchell creep model of Shanghai very soft silty clay[J]. Rock and Soil Mechanics,2000,21(4): 363-366.

        [10]王常明,王清,張淑華. 濱海軟土蠕變特性及蠕變模型[J]. 巖石力學(xué)與工程學(xué)報,2004,23(2): 227-230.WANG Chang-ming,WANG Qing,ZHANG Shu-hua.Creep characteristics and creep model of marine soft soils[J]. Chinese Journal of Rock Mechanics and Engineering,2004,23(2): 227-230.

        猜你喜歡
        數(shù)字圖像測量方法徑向
        淺探徑向連接體的圓周運動
        RN上一類Kirchhoff型方程徑向?qū)ΨQ正解的存在性
        基于PID+前饋的3MN徑向鍛造機控制系統(tǒng)的研究
        重型機械(2020年3期)2020-08-24 08:31:40
        一類無窮下級整函數(shù)的Julia集的徑向分布
        ARGUS-100 藝術(shù)品鑒證數(shù)字圖像比對系統(tǒng)
        基于迭代稀疏分解的介損角測量方法
        基于應(yīng)變原理隨鉆鉆壓測量方法研究
        基于塊效應(yīng)測度的JPEG數(shù)字圖像盲取證
        一種高效的PCB翹曲度測量方法
        基于壓電激振的彈性模量測量方法
        亚洲av资源网站手机在线 | 八戒网站免费观看视频| 国产高清精品自在线看| 一级午夜理论片日本中文在线| 东京热日本av在线观看| 亚洲国产av无码专区亚洲av| 99这里只有精品| 中文字幕乱码亚洲无线| 婚外情长久的相处之道| 亚洲国产成人片在线观看| 91日本精品国产免| av在线免费播放网站| 福利视频一区二区三区| 熟女性饥渴一区二区三区| 国产av成人精品播放| 91精品国产综合久久久蜜臀九色| 青青草 视频在线观看| 欧美怡红院免费全部视频| 99热成人精品国产免国语的| 在线观看免费视频发布白白色| 4455永久免费视频| 吃奶摸下的激烈视频| 亚洲成AV人片在一线观看| 蕾丝女同一区二区三区| 国产激情综合在线观看| 免费国产交换配乱淫| 国产精品视频免费一区二区三区 | 台湾佬中文网站| 午夜无码片在线观看影院| 人人爽亚洲aⅴ人人爽av人人片| 日韩美女av一区二区| 三年中文在线观看免费大全| 日韩高清毛片| 激情在线视频一区二区三区| 人人妻人人澡人人爽欧美一区| 无遮挡边摸边吃奶边做视频免费| 九一成人AV无码一区二区三区| 国产精品久久av色婷婷网站| 99精品国产一区二区三区不卡 | 日韩有码在线一区二区三区合集| 国产永久免费高清在线|