999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

振動碾作用下筑壩堆石料動力響應特性

2025-07-19 00:00:00張威余懿凡
人民長江 2025年6期
關鍵詞:振動質量

DOI:10.16232/j.cnki.1001-4179.2025.06.022

0 引言

堆石壩的壓實質量對其穩定性和耐久性至關重要[1-2],壓實不當會導致堆石壩強度與承載能力降低以及更大的沉降、滲透性和體積變化[3]。堆石壩碾壓中最常用的壓實機械是振動碾[4-5],這是一種振動和碾壓相結合的壓實機械,可有效壓縮堆石料,增加被壓實材料的剛度[6]

堆石壩施工質量通常采用規范要求的試坑取樣方法來控制[7-8],通過試坑試驗獲取采樣點的料性結果(如干密度、級配、孔隙率、含水率等)來代表一定方量填筑體的壓實質量。這種方式簡單易行,且較為直觀,但同時也具有片面性[9]。因此,國內學者針對傳統的原位試驗評估方法采取了多項改進措施。如王敬[10]提出通過建立沉降量與固體體積率回歸模型的方法,對壓實質量進行評估;徐平等[11]發現貫入儀可間接反映施工質量,可在測試的基礎上建立模型,得到貫入度與壓實質量之間的關系,實現壓實質量的快速評估;王東偉[12在討論不同壓實評估試驗的有效性及相關性的基礎上,提出了通過分析比較來選擇合理試驗的必要性;張智[13]通過試驗建立了載荷沉降的關系模型,從而獲得土壤的變形模量,以達到評估壓實度的目的。以上研究均針對傳統方法的弊端做出了改進,但是仍無法做到實時、連續的壓實質量評估。國外施工碾壓的質量評價研究主要集中在道路施工的路基壓實質量評價方面。Anderegg等[14]結合壓實機械GPS定位建立了智能壓實的反饋控制系統,提出以壓實機械的振動狀態評估路基的碾壓質量;Mooney等[15研究了碾壓控制參數的時空分布規律,嘗試采用連續壓實控制和智能碾壓機(CCC/IC)2種方法,通過參數分布客觀評價道路施工碾壓質量的控制效果。

從國內外開展的一系列研究和應用來看,傳統的面板堆石壩填筑碾壓質量控制方法受人為因素影響較大,存在管理不到位的環節,且試坑取樣檢測具有較大的離散性,以此來反映整個工作區的壓實狀況存在較大誤差,同時試坑試驗占用時間長,對壩體施工造成干擾[16]。為此本文以老撾南公1水電站面板堆石壩工程為例,通過構建振動碾和筑壩堆石料四自由度動力耦合分析模型,結合現場碾壓作用下堆石料動力響應原位監測試驗,開展振動碾-堆石體的動力響應分析,通過振動監測反饋堆石料碾壓質量,提出了一種基于堆石料動力響應的碾壓質量評價新方法,以避免在傳統的壩體填筑碾壓質量評價中將隨機取樣結果作為評價標準的缺陷,保證壩體碾壓質量并實現快速施工。

1 工程概況

1.1 南公1水電站概況

南公1水電站位于老撾南部阿速坡省內的南公河上,是中老兩國清潔能源戰略合作開發的重點項目。壩址距省城阿速坡公路里程約 53km ,距萬象約958km ,距越南峴港約 464km ,距泰國曼谷約1167km[17] 。南公1水電站裝機容量 160MW ,保證出力55.17MW,工程等別為Ⅱ等大(2)型,主要任務是發電,兼顧防洪、灌溉、交通、航運。工程由首部樞紐及引水發電系統組成,首部樞紐包括面板堆石壩、左岸溢洪道、左岸導流隧洞和右岸電站進水口[18];引水發電系統由電站進水口、引水隧洞、地下廠房洞室群、尾水隧洞等組成,引水發電建構筑物布置于右岸,為單洞、兩機布置形式。

南公1水電站主擋水建構筑物為混凝土面板堆石壩,布置于主河床,壩頂長 400m ,壩頂寬 8.8m ,壩頂高程 325.0m ,最大壩高 90m 。壩體上游壩坡1:1.4,下游壩坡1:1.35。電站的溢洪道布置于大壩左岸,主要由引渠段、閘室段和泄槽段(共設三級消力池)等組成。

1.2 面板堆石壩的堆石區施工工藝

南公1水電站采用混凝土面板堆石壩方案,溢洪道開挖料為大壩填筑供料。混凝土面板堆石壩以堆石體作為支承結構,堆石材料的質量以及施工質量是衡量混凝土面板堆石壩安全運行的主要指標,面板作為防滲結構,既能滿足抗滲性和耐久性的需求,又具有一定的柔性,可以很好地適應堆石體的變形,更好地保證壩體結構的穩定性和運行的安全性[19]。

溢洪道開挖料巖性單一,為三疊系中一下統流紋巖,干抗壓強度平均為 126.13MPa ,濕抗壓強度平均為 76.8MPa ,干彈性模量約為 84GPa 。堆石料填料的設計要求如下:最大粒徑不超過 800mm ,小于 5mm 的顆粒含量不超過 20% ,小于 0.075mm 的顆粒含量不超過 5% ,含泥量不大于 5% ,孔隙率不大于 22% ,壓實干密度不小于 2.12g/cm3 ,不均勻度系數不小于15,曲率系數為1\~3,主要填充參數如表1所列。

表1堆石料主要填充參數Tab.1Main filling parameters of rockfill

面板堆石壩是由堆石區、過渡區以及防滲體組成的,其中堆石區施工控制最為關鍵。該工程采用薄層碾壓法進行堆石體碾壓,這種方式能夠有效提升堆石體密實度,同時對于堆石體穩定性以及防滲結構變形的影響比較小。筑壩堆石料的水平碾壓選用 32t 振動碾,工作參數如表2所列。堆石料填筑層厚 110cm ,碾壓10遍,每臺振動碾的臺班生產效率為 2400~ 3 000m3. /臺班,需要碾壓設備不少于2臺。

表2振動碾主要工作參數Tab.2Main working parameters of vibrating roller

2 集總參數模型

利用集總參數模型分析碾壓層剛度與振動峰值之間的關系,論證基于堆石料動力響應進行碾壓質量評價的可行性。

2.1 模型建立

假設滾碾與堆石料在振動過程中始終處于接觸狀態,引入集總參數模型,建立四自由度DOF系統,如圖1所示。圖中 mf 為上部車輛由車架在其前軸上傳播的重量, md 為滾碾和內車架傳播的車輛重量; ms 為被壓縮的堆石質量; kdf 是橡膠阻尼器剛度; Cdf 為橡膠阻尼器的阻尼, Ks 和 Ks0 為堆石料壓實剛度, Cs 和 Cs0 為堆石料壓實阻尼; F0 為激振力; Zf 為車架位移; Zd 為滾碾位移; Zs 為堆石料位移。

該模型包含兩部分:上部為振動碾的車架和滾碾;下部分為正在壓縮的堆石料。振動激勵是由繞滾碾軸循環的偏心配重結構產生的,離心力 F 可由公式 F0= mecrecωec2 確定,其中 mecrec 為偏心塊的固定偏心力矩,ωec2 為振動壓實機的圓周激勵頻率。YZ32d型振動碾總質量 32t ,滾碾寬度接近 2 200mm ,直徑接近1700mm 。振動碾的滾碾質量( md )為 9120kg ,車架質量Ξ(mf) 為 12 000kg 。圓形激勵頻率為 28Hz ,垂直滾筒幅值為 1.83mm 。滾碾與車架之間的剛度系數 kdf= 3.4MN/m ,將 kdf 乘以損耗因子( n=0.16 ),再除以偏心轉動頻率,得到滾碾與車架之間的阻尼系數 Cdf 。堆石料模型參數 ms=0.2md,Ks=72MN/m , Cs=230 kN?s/m,Ks0=100MN/m,( Cs0 通過將 Ks0 乘以損耗因子( n=0.16; ,再除以偏心轉動頻率確定。

圖1振動碾與堆石料的四自由度集總參數模型Fig.1 Four-degree-of-freedomlumpedparametermodelofvibratingrollerand rockfill

通過力平衡法計算了接觸行為的運動方程,滾碾和堆石料的自由體圖見圖1。式(1)~(2)反映滾碾/堆石料接觸模式振動時滾碾和框架的動作。式(3)表示滾碾堆石料相互作用時的堆石料作用反力。

2.2 模型結果

振動峰值與碾壓層剛度的關系如圖2(a)所示,可以看出振動峰值與剛度呈正相關關系,且兩者同時增加。圖2(b)表示隨著振動碾碾壓遍數的增加,堆石料的密實度和振動速度也隨之增加,達到6次時趨于穩定和恒定。圖2(c)展示了振動壓實機與堆石料連續接觸的正弦波,振動速度隨著碾壓遍數的增加而逐漸增大。

圖2集總參數模型結果分析Fig.2Analysis of lumped parameter model results

通過對四自由度系統集總參數模型的分析,碾壓層的振動峰值隨著剛度的增加而逐漸增大,兩者之間存在對應關系。堆石料碾壓層的剛度在實際施工中測量比較困難,但振動測量比較方便。因此,可以用碾壓層的振動峰值來表示堆石料的剛度。隨著碾壓過程的進行,堆石料的密實度增加,最終趨于穩定。將振動碾和碾壓層簡化為一個四自由度的運動系統,碾壓層的剛度隨碾壓過程逐漸增大,趨于穩定。因此,在評價壓實質量時,堆石料壓實到一定程度后,壓實度難以繼續增加,剛度趨于穩定值,相應的振動峰值也趨于穩定,即可以認為,當振動峰值相同時,堆石料的碾壓效果是相同的。

3堆石料動力響應的測試方案

3.1 測試背景

面板堆石壩的碾壓質量對壩體的安全運行至關重要,如果碾壓密實度不夠,可能造成嚴重的不均勻沉降并形成滲漏通道[20]。且當壩體上下游填筑速度不同時,上下游會產生臨時內坡,內坡附近需預留一定寬度不碾壓,待下游填筑到相應標高后再碾壓,臨時邊坡的碾壓厚度厚于常規區域,更難滿足工程碾壓質量要求。

3.2 試驗程序

(1)傳感器預處理。由于傳感器的尺寸太小,為保證傳感器在鋪裝和找平過程中不會有較大的位移或傾斜,將傳感器包裹在尺寸為 15cm×15cm×15cm 的混凝土試件中,以確保傳感器位置和監測方向固定,同時還可以保護傳感器不受損壞,并在振動監測前進行傳感器功能試驗,確保傳感器可有效監測振動。

(2)傳感器埋設。傳感器埋設在軋制層之間的界面上,用鋼管保護數據線并引出斜坡。

(3)上下游兩側分層鋪裝。堆石壩進行填筑時需要對同一層進行階段性填筑,試驗方案及傳感器布設如圖3所示。堆石區分為I1區,I2區,Ⅱ1區,Ⅱ2區,Ⅲ區, ΔN 區, ΔV 區,碾壓順序為:I1區 區 I 2區 區 $$ Ⅲ區 區 V 區;測點1\~5號共5個,每個測點監測3個方向,其中1,2,3號布置在I區,4號,5號布置在Ⅱ區。不同層鋪裝完成后,根據振動監測結果開展后續分析。

圖3碾壓試驗方案示意(尺寸單位:m) Fig.3Schematic diagram of rolling test scheme

3.3 振動碾主要參數

原位監測試驗選取的振動碾型號為YZ32d,主要工作參數如表2所列,與現場施工選用機械一致。試驗區域碾壓遍數為10次,各區域均從左岸開始進行振動碾壓,振動碾運行速度如表3所列。

表3各區域振動碾運行速度 Tab.3 Vibration roller running speed in each area

3.4 振動曲線

不同層鋪裝完成后,根據測點的位置和上面的覆蓋層數,可以得到各區域的振動監測數據,為后續分析提供數據支撐。典型的振動時程曲線如圖4所示。

4堆石料的動力響應特性

4.1 I1區振動試驗

對I1區進行堆石填筑和碾壓,在末端邊界處,3號,2號,1號測點從邊界線處依次以 0.3,1.3,2.7m 的間隔排列,試驗碾壓范圍為 2. 1m×10. 0m (寬 × 長),1號測點位于碾壓區下方深度 1.0m 、距試驗區短邊 4.0m 長邊 0.4m 處。分別統計各碾壓次數下各測點的振動速度峰值,見表4和圖5。由監測結果可知,1號測點的振動峰值明顯大于2號和3號測點,1號測點的振動速度峰值主要在 10~20cm/s 之間,2號和3號測點的振動速度峰值主要在 2~4cm/s 之間。

圖4典型的振動時程曲線(豎直向) Fig.4Typical vibration history curve (vertical)
表4碾壓振動速度峰值與I1區域碾壓次數的關系 Tab.4Relationship between the peak value of rolling vibration velocity and thenumber ofrollingin I1 region cm/

隨著碾壓次數的增加,各測點的振動峰值波動較大。分析其主要原因如下: ① 堆石料最大粒徑為80cm,不是均勻物料,即使密實度增加,對監測點振動峰值的影響也較小; ② 32t振動壓實機在碾壓過程中振動幅值不穩定,對碾壓層產生的激振載荷波動較大。但在豎直向上,隨著碾壓次數的增加,振動峰值呈現波動上升后趨于穩定的規律。如圖5所示,在正下方1號測點的豎直向上,當碾壓次數達到6次時,振動峰值從 10cm/s 增大到 17cm/s ,然后在 17cm/s 上下波動;同樣斜向下的2號測點豎直向振動峰值從2.8cm/s 增加到 3.0cm/s ,并趨于穩定,3號測點豎直向振動峰值從 2.5cm/s 增加到 3.0cm/s ,并趨于穩定。

4.2 Ⅱ1區域振動試驗

4號測點位于1號測點正上方,5號測點位于2號測點正上方,試驗碾壓范圍為 2.1m×10.0m (寬 × 長),4號測點位于碾壓區下方深度 1.0m 、距試驗區短邊 4.0m 長邊 1.1m 處。對 I1 區進行堆石填筑和碾壓,統計各碾壓次數下各測點的振動峰值,見表5和圖6。由監測結果可知,在水平徑向上,測點1的峰值為 2.36~4.12cm/s ,測點4的峰值為 4.88~7.65cm/ s,測點5的峰值為 3.22~10.56cm/s ;在水平切向上,測點1的峰值為 4.32~5.95cm/s ,測點4的峰值為7.99~11.20cm/s ,測點5的峰值為 3.29~6.23 cm/s ;在豎直向上,測點1的峰值為 7.46~9.89cm/s 測點4的峰值為 14.11~25.51cm/s ,測點5的峰值為3.87~5.40cm/s, 0

圖5 I1區域碾壓振動峰值與碾壓次數的關系 Fig.5Relationship between the peak value of rolling vibrationand thenumber ofrollingvibrationinI1 region

對于正下方的測點1和測點4,在豎直向上,隨著碾壓次數的增加,振動峰值有先增大后趨于穩定的趨勢

4.3 Ⅱ2區振動試驗

對 I2 層堆石料進行碾壓,試驗碾壓面積為 2.1m ×10.0m (寬 × 長),5#測點位于碾壓區下方深度1.0m 、距試驗區短邊 4.0m 、長邊 0.8m 處。統計各碾壓次數下各測點的振動峰值,見表6和圖7。在豎直向,2號測點的峰值為 9.05~10.48cm/s,3 號測點的峰值為 10.65~14.26cm/s,5 號測點的峰值為 22.64~ 31.40cm/s 。其中,2號和3號測點的峰值隨碾壓次數的增加有先上升后趨于穩定的趨勢,而5號測點的這一規律并不明顯。

表5碾壓振動速度峰值與II1區域碾壓次數的關系 Tab.5Relationship between the peak value of rolling vibration velocity and the number of rollinginII1 region cm/:
圖6 Ⅱ1區域碾壓振動速度峰值與碾壓次數的關系 Fig.6Relationship between the peak value of rolling vibration velocity and the number of rolling vibrationinⅡ1 region
表6碾壓振動峰值與ⅡI2區碾壓次數統計
圖7 Ⅱ2區域碾壓振動速度峰值與碾壓次數的關系 Fig.7 Relationship between the peak value of rolling vibrationvelocityandthenumberofrollingvibration in Ⅱ 2 region

4.4 碾壓質量評價

分別總結各區域的碾壓振動情況,如表7和圖8所示。碾壓I1區時的1號測點數據、碾壓Ⅱ1區時的4號測點數據、碾壓 區時的5號測點數據,均為碾壓厚度為 1.0m 時的監測數據,且測點均處于正常碾壓狀態,測點振動速度波動劇烈。與 1.0m 的鋪裝厚度相比,堆石料最大粒徑為 0.8m ,因此堆石料填筑的均勻性及振動壓實機的振動幅值穩定性對碾壓振動峰值的影響更為復雜和隨機。

表7各區域振動速度峰值Tab.7 Vibration velocitypeak in each region cm/s
圖8各區域振動速度峰值比較 Fig.8Comparison of vibration velocity peak value in each region

碾壓Ⅱ1區時的1號測點數據和碾壓 I2 區時的2號測點數據,均為碾壓厚度為 2.0m 時的監測數據,且測點均處于正常碾壓狀態。由監測結果可以看出當碾壓次數達到6次時,振動峰值趨于穩定,Ⅱ1區進行碾壓時,1號測點 7~10 次的振動峰值平均值為9.37cm/s , I2 區進行碾壓時,2號測點 7~10 次的振動峰值平均值為 9.85cm/s ,兩者之間的差異不明顯。從振動峰值來看,1號測點與2號測點處的振動峰值一致,即堆石料正常碾壓處與臨時邊坡處的碾壓質量相同,滿足碾壓質量要求。

5結論

本文提出了在堆石壩碾壓過程中,振動壓實機對堆石料的動態響應。通過確定不同堆石區域及不同深度的振動峰值,評價堆石料壓實質量,主要結論如下:

(1)通過模型分析可知,隨著振動碾碾壓次數的增加,堆石料的壓實質量增加,最終趨于穩定。碾壓層振動峰值與剛度呈正相關關系,兩者均隨著碾壓次數的增加而增大,反之亦然

(2)對比堆石料正常碾壓處與臨時邊坡處的碾壓振動峰值,發現兩者之間的差值隨著碾壓次數的增加而減小。對Ⅱ1區進行碾壓時,1號測點 7~10 次的振動峰值平均值為 9.37cm/s ,對 I2 區進行碾壓時,2號測點 7~10 次的振動峰值平均值為 9.85cm/s 。從振動峰值來看,臨時邊坡處的碾壓質量與其余堆石區一致,符合質量要求。

(3)本文提出了一種基于堆石料動力響應的碾壓質量評價新方法,有助于保證壩體碾壓質量并實現快速施工。

參考文獻:

[1] 崔博,鐘航,王佳俊,等.基于CA-RL的堆石壩填筑倉面碾壓作業動態路徑規劃[J].水利學報,2024,55(3):253-265.

[2] 唐德勝,謝鵬飛,馬子茹,等.考慮施工質量與填筑進度的面板壩精細有限元分析[J].人民長江,2024,55(10):141-148,165.

[3] 劉東海,王愛國,柳育剛,等.基于碾輪振動性態分析的土石壩壓實質量實時監測與評估[J].水利學報,2014,45(2):163-170.

[4] 王櫻峻,趙琳,雷顯陽,等.抽水蓄能電站高瀝青混凝土面板堆石壩設計[J].人民長江,2022,53(1):148-153.

[5] 楊衛甲,岳朝俊,陳道春,等.萬花溪水庫瀝青混凝土心墻現場鋪筑試驗[J].人民長江,2022,53(增1):96-100.

[6] 肖浩波.面板堆石壩變形特性及分析模型研究[D].南京:河海大學,2006.

[7] 蔡加興,方德揚.堆石料相對密度控制法試驗及檢測方法應用研究[J].人民長江,2019,50(12):136-141.

[8] 劉東海,吳優.實時監控下堆石壩壓實質量模糊綜合評估[J].水力發電學報,2019,38(3):142-153.

[9]王曉玲,周龍,任炳昱,等.基于實時監控的堆石壩碾壓質量二元耦合評價[J].水力發電學報,2015,34(2):164-170.

[10]王敬.土石混填路基壓實質量快速評定方法研究[D].上海:上海交通大學,2007.

[11]徐平,王復明,陳忠平.DCP在高速公路路基壓實質量評價中的應用[J].路基工程,2008(4):10-12.

[12]王東偉.瀝青路面現場壓實質量評估標準探討[J].重慶交通大學學報(自然科學版),2009,28(2):268-272.

[13]張智.客運專線鐵路路基變形模量檢測與路基壓實質量評價分析[J].鐵道建筑,2010(7):84-86.

[14]ANDEREGGR,KAUFMANN K.Intelligent compaction with vibratoryrollers :feedback control systems in automatic compaction and compac-tion control[J].TransportationResearch Record,2004,1868(1) :124-134.

[15]MOONEYAM,FACASWN,RINEHARTVR.Developmentand eval-uation of relative compaction specifications using roller-based meas-urements[J].Geotechnical TestingJournal,2011,34(6):634-642.

[16]劉彪,趙宇飛,陳祖煜,等.基于碾壓波速的堆石壩壓實質量實時監測指標的研究[J].中國水利水電科學研究院學報(中英文),2022,20(1) :20 -29.

[17]陳育權,劉廣鑫,柯昌佳,等.南公1水電站溢洪道保護層開挖孔底消能-聚能結構優選研究[J].爆破,2020,37(4):10-17.

[18]吳世然,劉廣鑫,陳育權,等.老撾南公1水電站非襯砌溢洪道巖體爆破開挖技術[J].爆破,2021,38(3):31-37.

[19]田聰,汪華東.老撾南公1水電站大壩填筑施工質量控制[J].四川水力發電,2020,39(增2):73-76.

[20]陳祖煜,趙宇飛,鄒斌,等.大壩填筑碾壓施工無人駕駛技術的研究與應用[J].水利水電技術,2019,50(8):1-7.

(編輯:郭甜甜)

Abstract:Poorcompactionqualitycanleadtotheformationofseepagechannelsandunevensetlementofdams,weakening its stabilityandlong-termdurability.Therefore,itiscrucialtoevaluatethecompactionqualityofrockfilldamsBasedontheNam Kong-1HydropowerStationproject inLaos,afour-degree-of-freedomdynamiccouplinganalysis modelwasestablishedto studytheinteractionbetweenthevbratingolerandterockfiling materials.This model,combined within-situmontoringtests of thedynamicresponseofockfllmaterialsunderrolercompaction,enabledadynamicresponseanalysisandcompactioquality evaluationthrough vibrationcharacteristicsofthecompactedlayer.ThemodelanalysisresultsindicatethatthevibrationpeakvaluesofrockfillmaterialswrepositivelycorelaedwithstiffssAsthenumbrofollngpassesireases,tedensityofecompactedlayerrises,leadingtoagradualincreaseintheiffessofterockfillmaterialsandacorespondingincreaseinvration peak values.Fieldtestresultsfurtherdemonstratedthatasthenumberofroling pases increase,thevibration peaksinthecompactionarea progresivelyrise.Thecompactionqualitymets therequired standardswhen the numberofrolingpassesexceds seven.Thisstudyproposesanovelmethodforevaluatingcompactionqualitybasedonthedynamicresponseofrockfillmaterials, which helps ensure dam compaction quality and facilitates rapid construction.

Key Words:faced rockfill dam;rockfill material;compacting;vibration monitoring;dynamic model

猜你喜歡
振動質量
振動的思考
科學大眾(2023年17期)2023-10-26 07:39:14
噴水推進高速艇尾部振動響應分析
“質量”知識鞏固
質量守恒定律考什么
This “Singing Highway”plays music
做夢導致睡眠質量差嗎
振動攪拌 震動創新
中國公路(2017年18期)2018-01-23 03:00:38
中立型Emden-Fowler微分方程的振動性
關于質量的快速Q&A
質量投訴超六成
汽車觀察(2016年3期)2016-02-28 13:16:26
主站蜘蛛池模板: 亚洲最新地址| 国产精品无码制服丝袜| 精品国产免费观看| 欧美一区二区三区国产精品| 高清久久精品亚洲日韩Av| 最新亚洲人成无码网站欣赏网 | 国产免费精彩视频| 无码'专区第一页| 伊人五月丁香综合AⅤ| 亚洲成年人片| 无码中文字幕精品推荐| 91香蕉国产亚洲一二三区 | 呦女精品网站| 中文字幕第1页在线播| 青青青视频91在线 | аv天堂最新中文在线| 国产精品无码AⅤ在线观看播放| 免费一级无码在线网站| 日韩精品一区二区三区大桥未久| 国产欧美视频在线| 亚洲人在线| 伊人婷婷色香五月综合缴缴情| jizz在线观看| 美女啪啪无遮挡| 亚洲无限乱码一二三四区| 日韩欧美国产中文| 999精品视频在线| 亚洲无码视频图片| 国产欧美日韩18| 国产无码精品在线播放| 视频二区亚洲精品| 欧美国产日韩在线播放| 亚洲国产成人综合精品2020| 欧美精品色视频| 五月综合色婷婷| 亚洲中文字幕23页在线| 55夜色66夜色国产精品视频| 在线欧美a| 精品午夜国产福利观看| 狼友视频一区二区三区| 伊人久久久久久久久久| 国产91线观看| 91青青草视频| 欧美日韩免费观看| 国模视频一区二区| 91青青在线视频| 97视频免费看| 9啪在线视频| 欧美日韩高清在线| 毛片免费视频| 国产免费福利网站| 天天躁夜夜躁狠狠躁躁88| 国产综合另类小说色区色噜噜| 亚洲精品第1页| 无码精品国产dvd在线观看9久| 色偷偷男人的天堂亚洲av| av一区二区三区在线观看| 日本高清免费一本在线观看| 免费人成在线观看视频色| 蝴蝶伊人久久中文娱乐网| 亚洲第一精品福利| av尤物免费在线观看| 伊人成人在线视频| 国产精品七七在线播放| 国产麻豆福利av在线播放| 亚洲国产精品美女| 亚洲一区二区精品无码久久久| 国内精品91| 国产精品久久久久久久伊一| 国产区免费| 国产乱人伦AV在线A| 又粗又硬又大又爽免费视频播放| 亚洲综合二区| 伊人久热这里只有精品视频99| 亚洲毛片网站| 亚洲成肉网| 国产在线拍偷自揄观看视频网站| 精品国产电影久久九九| 浮力影院国产第一页| 国产精品毛片一区视频播| 人妻21p大胆| 久久婷婷国产综合尤物精品|