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

    高速鐵路CFG樁-土復(fù)合地基減振性能研究

    2023-08-08 01:04:24高廣運張繼嚴畢俊偉
    中南大學學報(自然科學版) 2023年6期
    關(guān)鍵詞:振動效果模型

    高廣運 ,張繼嚴 ,畢俊偉

    (1.同濟大學 土木工程學院,上海,200092;2.同濟大學 巖土及地下工程教育部重點實驗室,上海,200092;3.華南理工大學 土木與交通學院,廣東 廣州,510641;4.廣州市設(shè)計院集團有限公司 巖土與地下空間院,廣東 廣州,510620)

    在我國四通八達的高速鐵路網(wǎng)中,樁-土復(fù)合地基憑借其穩(wěn)定性好、工期短、造價低等優(yōu)勢而成為軟基處理的主要形式,其中,CFG樁-土復(fù)合地基具有較大占比[1]。源于高鐵線路對軌道平順性的嚴苛要求,樁-土復(fù)合地基的研究多聚焦于土拱效應(yīng)、動應(yīng)力分布及沉降[2-3]問題?,F(xiàn)階段,軌道交通引起的環(huán)境振動已成為軌道交通發(fā)展中亟需解決的關(guān)鍵問題,高速鐵路樁-土復(fù)合地基的地面振動問題開始受到關(guān)注。

    現(xiàn)場實測可準確反映高鐵荷載下樁-土復(fù)合地基應(yīng)力分布和振動的傳播衰減規(guī)律。肖宏等[4-5]在遂渝線無砟軌道綜合試驗段進行動車組和貨車荷載下的地基振動測試,探討了樁-網(wǎng)復(fù)合地基動應(yīng)力和振動加速度沿深度的衰減規(guī)律;ZHAI等[6-7]分別在京滬高鐵蘇州東站和蚌埠南站附近開展環(huán)境振動實測,分析不同車速下樁-土復(fù)合地基的地面振動特性及傳播衰減規(guī)律;WANG等[8]在京滬高鐵線昆山段進行地面振動實測,發(fā)現(xiàn)無砟軌道線列車運行引起的振動主頻為相鄰車廂中間距離對應(yīng)頻率的整數(shù)倍?,F(xiàn)場實測可為相關(guān)研究提供可靠的數(shù)據(jù)支持,但現(xiàn)場實測活動會受到限制而不能大規(guī)模開展。模型試驗在荷載施加以及傳感器布置方面更具靈活性,可以開展針對性的研究。劉漢龍等[9-11]利用大比尺模型試驗,研究高鐵荷載下X 形樁-網(wǎng)復(fù)合地基的動力特性。陳仁朋等[12]基于全比尺試驗?zāi)P?,研究長期列車移動荷載下樁網(wǎng)加筋復(fù)合地基土拱效應(yīng)和樁土動應(yīng)力演化規(guī)律。模型試驗的模型尺寸受到試驗場地限制,無法開展距軌道較深或較遠處的振動研究。

    在數(shù)值研究方面,翟婉明等[13]將攪拌樁復(fù)合地基簡化為等效復(fù)合土層,發(fā)現(xiàn)樁-土復(fù)合地基能有效減小地面振動。FENG等[14]采用相同的簡化計算方法,建立ABAQUS 三維模型,分析了復(fù)合地基加固深度和加固范圍對地面振動的影響;THACH 等[15-17]利用三維有限元數(shù)值模型研究高速鐵路現(xiàn)澆大直徑混凝土管樁(PCC)復(fù)合地基的振動特性;高廣運等[18-19]利用ABAQUS 建立CFG 樁-土復(fù)合地基三維有限元計算模型,分析了土體阻尼和樁體模量等對振動的影響;TANG等[20]建立三維有限元模型,研究準靜態(tài)列車荷載下車速和布樁形式(樁間距)對樁-土復(fù)合地基動應(yīng)力和豎向變形的影響。李丹陽等[21]將地基土視為彈塑性材料,研究了樁徑、樁長和樁間距對樁承式地基基床表層振動的影響;LI等[22]采用LS-DYNA建立三維耦合模型,分析了不同車速條件下地基加固前后的地面振動響應(yīng),指出樁-土復(fù)合地基可減小地面振動。FERNANDEZ-RUIZ 等[23]建立三維仿真模型,研究設(shè)置石柱樁對軟土地基臨界速度的提升效果,著重分析了樁長對臨界速度的影響;CASTANHEIRA-PINTO等[24]基于樁-土復(fù)合地基頻散圖和軌道頻散曲線提出軌道-樁土復(fù)合地基耦合系統(tǒng)臨界速度簡化計算方法;WANG 等[25]建立了三維數(shù)值模型,對比了樁加固地基與未加固地基振動速度差異,并探討了地面局部放大現(xiàn)象的成因。

    基于商業(yè)數(shù)值軟件的三維有限元分析對計算機性能要求較高,計算時間長,且上述研究中均未全面討論樁土復(fù)合地基設(shè)計參數(shù)對地面振動和減振性能的影響。YANG 等[26]將更為高效的2.5 維有限元法引入軌道交通環(huán)境振動研究;高廣運等[27]采用等代樁墻對樁-土復(fù)合地基進行等效計算,建立2.5維有限元計算模型,對比分析了高鐵荷載下CFG樁-土復(fù)合地基與無樁地基的地面振動特性,但文中未討論設(shè)計參數(shù)對地面振動的影響。

    為分析CFG樁設(shè)計參數(shù)對樁-土復(fù)合地基減振性能的影響,本文基于2.5維有限元方法,在沿列車運行方向采用等代樁墻對樁-土復(fù)合地基進行等效,建立軌道-路堤-CFG 樁-土復(fù)合地基系統(tǒng)的2.5 維有限元計算模型,分析CFG 樁-土復(fù)合地基和無樁地基的振動特性,研究樁徑、相鄰樁中心距(樁間距)和樁體剛度對CFG 樁-土復(fù)合地基減振性能的影響。

    1 模型的建立及驗證

    1.1 CFG樁-土復(fù)合地基2.5維有限元控制方程

    單次列車運行引起的地基振動為小應(yīng)變問題[28]。本文視路堤、CFG 樁和地基土體為均質(zhì)黏彈性介質(zhì),其本構(gòu)關(guān)系如下:

    式中:σij為應(yīng)力張量;εij=(ui,j+uj,i)/2,為應(yīng)變張量;ui為位移;δij為克羅內(nèi)克函數(shù);e=ui,i為體應(yīng)變;i和j表示空間坐標x、y和z軸;λd和μd為考慮材料阻尼的Lamé常數(shù),定義如下:

    式中:E和v分別為彈性模量和泊松比;β為阻尼比;i為虛數(shù)單位。

    采用位移表示的運動方程,

    式中:ρ為密度,上標“··”表示對時間的二次導(dǎo)數(shù)。

    2.5 維有限元法是將三維瞬態(tài)動力問題通過雙重Fourier變換,轉(zhuǎn)變?yōu)轭l域-波數(shù)域內(nèi)的問題進行解答。圖1所示為軌道-路堤-CFG樁-土復(fù)合地基系統(tǒng)2.5維有限元模型,高速列車以車速c沿x方向行駛。本文定義對時間t和坐標x的雙重Fourier 變換分別如下:

    圖1 軌道-路堤-CFG樁-土復(fù)合地基系統(tǒng)2.5維有限元模型Fig.1 Schematic of 2.5D FE model of track-embankment-CFG pile-soil composite subgrade system

    式中:上標“-”和“~”分別表示頻域和波數(shù)域內(nèi)的變量;ω為圓頻率;ξx為對應(yīng)x方向的波數(shù)。

    按式(7)對式(6)進行雙重Fourier 變換,采用四節(jié)點等參單元對計算模型離散,基于Galerkin法得頻域-波數(shù)域內(nèi)的有限元控制方程,其矩陣形式如下:

    實際工程中,CFG 樁沿列車運行方向呈周期性分布,無法滿足2.5維有限元方法關(guān)于幾何形狀和材料性質(zhì)沿x向保持不變的基本假設(shè)。為此,本文采用如圖2所示的等代樁墻對CFG樁-土復(fù)合地基進行等效計算。等代樁墻的等效彈性模量EPS計算公式如式(10)所示。單次列車運行引起的樁土相對滑移可忽略不計[22]。因此,模型中等代樁墻單元與地基土體單元采用共節(jié)點方式耦合。為消除有限元模型截斷邊界處反射波的影響,本文采用無限元邊界[29]模擬振動波在模型邊界外的傳播。

    圖2 等代樁墻示意圖Fig.2 Schematic of equivalent pile wall

    式中:EP和ES分別為樁體和樁間土的彈性模量。

    1.2 軌道模型與列車荷載

    研究表明,是否考慮軌道結(jié)構(gòu)各部分耦合對列車荷載下振動響應(yīng)影響較小[30]。因此,本文假定板式軌道結(jié)構(gòu)在列車荷載下發(fā)生整體豎向變形,將其視為鋪設(shè)在路堤上的歐拉梁,其頻域-波數(shù)域內(nèi)的控制方程如下:

    式中:ErIr為軌道抗彎剛度;mr為軌道單位長度質(zhì)量;ur為軌道的豎向位移;為軌道與路堤間的相互作用力;為列車移動荷載。

    列車荷載是由一系列的輪對荷載組成,對于由N節(jié)車廂組成的列車,考慮自振的連續(xù)軸重移動荷載χ(ξx)在頻域波數(shù)域內(nèi)的表達式為[31]

    式中:pn為第n節(jié)車廂的軸重;an和bn分別為同一轉(zhuǎn)向架下的軸距和同一車廂下的轉(zhuǎn)向架間距;Lm為第m節(jié)車廂長度;ω0為列車荷載自振頻率。

    1.3 模型驗證

    ZHAI等[6]對京滬高鐵蘇州東站附近的樁-土復(fù)合地基開展地面振動實測,列車為CRH380AL 動車組,車速為350 km/h。為驗證模型的正確性,依據(jù)實測路段的軌道結(jié)構(gòu)、路基形式和地質(zhì)條件[6]建立2.5維有限元計算模型。本文模型計算得到的地面振動加速度峰值與實測結(jié)果如圖3 所示。由圖3可知,計算得到的地面振動加速度峰值和衰減規(guī)律與實測結(jié)果較吻合,驗證了模型的正確性和可靠性。

    圖3 地面豎向加速峰值沿地表的衰減曲線Fig.3 Attenuation curves of peak ground vertical acceleration with distance from track centerline

    2 CFG樁-土復(fù)合地基減振特性

    選取單線路堤進行建模分析。地基由2層土體組成,總深度為20.0 m,寬度為80.0 m。路堤高度為3.0 m,頂面寬度為8.6 m,坡度為1.0∶1.5。CFG樁-土復(fù)合地基中,墊層為厚度為0.6 m 的碎石墊層夾1層土工格柵,土工格柵彈性模量設(shè)為碎石墊層彈性模量的1.4 倍[27]。CFG 樁長為15.0 m,樁徑為0.5 m,樁間距為2.0 m,采用正方形布樁。如圖1 所示,基于問題的對稱性,選取半結(jié)構(gòu)進行建模,計算參數(shù)如表1所示。模型左側(cè)采用軸對稱邊界,右側(cè)采用無限元邊界,底部采用固定邊界。列車荷載選取CRH380AL動車組列車。

    表1 動力計算參數(shù)Table 1 Dynamic calculation parameters

    當車速為300、350、400、450、470和500 km/h時,無樁地基與CFG樁-土復(fù)合地基地面豎向振動加速度峰值沿地表的衰減曲線如圖4 所示。由圖4可知:當車速低于地基表層土的Rayleigh波速(cR=398.25 km/h)時,CFG 樁-土復(fù)合地基與無樁地基的地面振動加速度峰值衰減曲線較為平滑。當車速接近Rayleigh波速時,2 類地基中地面振動加速度峰值沿地表的衰減表現(xiàn)出波動性,且隨車速提高,波動性愈加強烈,其中CFG樁-土復(fù)合地基的波動性明顯弱于無樁地基的波動性。對于無樁地基(圖4(a)),當車速為400 km/h 時(接近表層土Rayleigh 波速),距軌道較近處的振動加速度驟然增大,出現(xiàn)類共振現(xiàn)象,說明無樁地基的臨界速度為表層土的Rayleigh波速。然而,對于CFG樁-土復(fù)合地基(圖4(b)),其類共振現(xiàn)象發(fā)生在車速為470 km/h時,表明CFG樁-土復(fù)合地基提高地基的臨界速度。

    圖4 各車速條件下豎向加速度峰值沿地表的衰減曲線Fig.4 Attenuation curves of peak vertical acceleration along ground surface at various train speeds

    距軌道中心3.0 m 和30.0 m 處CFG 樁-土復(fù)合地基和無樁地基地面振動豎向加速度峰值隨車速的變化曲線如圖5所示。從圖5(a)可見:距軌道中心3.0 m 處,無樁地基和CFG 樁-土復(fù)合地基分別在車速400 km/h 和470 km/h 時軌道近處出現(xiàn)類共振現(xiàn)象。當車速遠未達到兩類地基的臨界車速時,CFG樁-土復(fù)合地基的豎向振動加速度比無樁地基的略小,減振效果較弱。當車速接近無樁地基臨界速度時,CFG樁-土復(fù)合地基的豎向振動加速度明顯小于無樁地基的豎向振動加速度,此時減振效果最顯著。在CFG樁-土復(fù)合地基臨界車速一定范圍內(nèi),CFG樁-土復(fù)合地基的豎向振動加速度峰值大于無樁地基的豎向振動加速度峰值,此時無減振效果。因此,在距軌道較近處,CFG樁-土復(fù)合地基的減振效果由車速和共振條件共同主導(dǎo)。

    圖5 兩類地基距軌道中心3.0 m和30.0 m處地面豎向加速度峰值隨車速的變化Fig.5 Variation of peak vertical accelerations with vehicle speed at 3.0 m and 30.0 m from track centerline for two types of subgrade

    從圖5(b)可見:距軌道中心30.0 m處,兩類地基的地面豎向振動加速度峰值均隨車速提高不斷增大,振動響應(yīng)不再受共振條件影響。各車速條件下CFG樁-土復(fù)合地基的豎向振動加速度峰值均明顯比無樁地基的小。當車速為300.0 km/h時,無樁地基加速度峰值為0.072 m/s2,CFG 樁加固后加速度峰值減小62.5%;當車速達到500.0 km/h 時,無樁地基加速度峰值為0.913 m/s2,CFG 樁加固后加速度峰值減小53.1%。這表明CFG 樁-土復(fù)合地基能有效減小距軌道較遠處的地面振動。

    圖6所示為各車速條件下無樁地基與CFG樁-土復(fù)合地基地面豎向加速度振級隨與軌道中心間距的變化曲線,其中灰色虛線代表《城市區(qū)域環(huán)境振動標準》[32]規(guī)定的鐵路干線兩側(cè)區(qū)域環(huán)境振動限值80.0 dB。由圖6 可見:兩類地基的地面振動振級隨與軌道中心間距的變化規(guī)律與加速度峰值曲線相似。由圖6(a)可知:在距軌道中心30.0 m范圍內(nèi),各車速條件下無樁地基的地面振動振級均高于振動限值80.0 dB,不滿足振動控制要求,需要采取減振措施。從圖6(b)可見,CFG 樁-土復(fù)合地基的減振作用導(dǎo)致各車速條件下在距軌道中心一定距離處的振級降至規(guī)范限值以下?,F(xiàn)階段,我國高速鐵路的運營速度一般不高于400 km/h,當車速≤400 km/h 時,距CFG 樁-土復(fù)合地基軌道中心25.0 m 以外的地面振動均滿足振動控制要求。因此,采用CFG樁-土復(fù)合地基能顯著減小高速列車運行引起的地面振動,可有效降低高鐵沿線各類建筑物的環(huán)境振動控制成本。

    圖6 各車速下地面振動豎向加速度振級隨與軌道中心距離的變化曲線Fig.6 Variations of vertical acceleration levels with distance from track centerline at different train speeds

    3 CFG樁設(shè)計參數(shù)對減振性能影響

    選取不同樁徑、樁間距和樁體剛度,分析CFG 樁設(shè)計參數(shù)對樁-土復(fù)合地基地面振動的影響,各計算工況如表2所示。參考我國高速列車現(xiàn)階段運營速度及下階段的發(fā)展規(guī)劃,選取車速為300、350、400 和450 km/h。計算模型及其他計算參數(shù)同第2節(jié)所述。

    表2 各工況CFG樁設(shè)計參數(shù)Table 2 CFG pile design parameters for each cases

    3.1 樁徑

    選取樁間距2.0 m,對比分析樁徑0.3、0.5 和0.8 m這3種工況(即表2工況1~3,樁土剛度比RPS=EP/ES,EP和ES分別為樁體和土體的彈型模量)下的地面豎向振動加速度峰值,研究樁徑對CFG樁-土復(fù)合地基減振性能的影響。

    圖7 所示為各車速條件下,3 種樁徑工況地面豎向振動加速度峰值隨與軌道中心間距的變化曲線,此外,添加無樁地基工況進行對比分析。由圖7(a)和7(b)可知:當車速低于地基表層土Rayleigh波速時,各樁徑條件下豎向振動加速度峰值沿地表的衰減曲線較為光滑。樁徑越大,地面各處的豎向加速度峰值越小,減振效果越強。各樁徑條件下,加速度峰值減小幅度隨與軌道中心距離增加而呈先增大后減小的趨勢,這表明CFG樁-土復(fù)合地基的減振效果隨與軌道中心距離增加而先增強后略有減弱。整體上,距軌道較遠處的減振效果優(yōu)于軌道近處的減振效果。另外,隨樁徑增大,振動加速度峰值的減小幅度逐漸降低。當樁徑≥0.5 m(1/4倍樁間距)時,進一步增大樁徑難以提升減振效果。

    圖7 不同樁徑條件下地面振動豎向加速度峰值Fig.7 Variations of peak ground vertical accelerations with distance from track centerline for different pile diameters

    如圖7(c)和7(d)所示,當車速接近或高于地基表層土Rayleigh波速時,各樁徑條件下振動加速度峰值的衰減曲線呈現(xiàn)出不同于無樁地基的波動性,使得減振效果隨與軌道中心距離的變化同樣具有波動性,且各樁徑條件下減振效果的波動性各有差異。整體上,減振效果隨樁徑增大逐漸增強,且樁徑限值仍為0.5 m(1/4 倍樁間距)。值得注意的是,軌道近處的減振效果由車速和共振條件共同主導(dǎo),導(dǎo)致車速450 km/h時CFG樁-土復(fù)合地基在距軌道近處的振動加速度峰值比無樁地基的略大。此外,由圖7分析得出,車速450 km/h時各樁徑條件下的減振效果均小于其他車速的減振效果。

    3.2 樁間距

    考慮樁徑為0.5 m,對比分析樁間距為3 倍樁徑(1.5 m)、4倍樁徑(2.0 m)和5倍樁徑(2.5 m)3種工況(表2 工況4~6)的地面振動豎向加速度峰值變化曲線(圖8),研究樁間距對CFG 樁-土復(fù)合地基減振性能的影響。

    圖8 不同樁間距條件下地面振動豎向加速度峰值Fig.8 Variations of peak ground vertical accelerations with distance from track centerline for different pile spacing

    由圖8 可知,各車速條件下,隨樁間距減小,地面豎向振動加速度峰值逐漸減小,減振效果逐漸增強。同時,隨樁間距不斷減小,對減振效果的提升逐漸減小,當樁間距≤4 倍樁徑時,進一步縮減樁間距難以提高CFG 樁-土復(fù)合地基的減振效果。

    3.3 樁體剛度

    設(shè)置樁徑0.5 m,樁間距2.0 m,考慮樁體彈性模量分別為1.0、5.0、10.0 和15.0 GPa,對應(yīng)樁土剛度比RPS分別為13.9、69.4、138.9和208.3(表2工況7~10),通過樁土剛度比分析樁體剛度對CFG樁-土復(fù)合地基減振性能的影響。

    圖9 所示為不同樁土剛度比RPS條件下,地面豎向振動加速度峰值沿地表的衰減曲線。各車速條件下,隨樁土剛度比增大,地面加速度峰值逐漸減小,減振效果逐漸增強。當車速低于地基表層土Rayleigh 波速時(見圖9(a)和(b)),各樁土剛度比RPS條件下的減振效果隨與軌道中心距離的變化規(guī)律與樁徑工況(圖7)相似。當車速接近或高于Rayleigh波速時(如圖9(c)和(d)),各樁土剛度比RPS條件下振動加速度峰值沿地表衰減曲線的波動性相似,導(dǎo)致不同樁土剛度比下的減振效果沿地表的波動性趨于一致。當樁土剛度比RPS≥138.9 時,進一步提高樁土剛度比RPS難以顯著增強減振效果。值得注意的是,當車速為450 km/h 時,CFG樁-土復(fù)合地基在距軌道近處的振動加速度比無樁地基的大,且當RPS=69.4(EP=5.0 GPa)時,振動加速度最大,說明此條件下復(fù)合地基的臨界車速最接近450 km/h,表明CFG樁-土復(fù)合地基的臨界車速受樁體剛度影響。

    圖9 不同樁土剛度比時地面振動豎向加速度峰值Fig.9 Variations of peak ground vertical accelerations with distance from track centerline for different pile-soil stiffness ratio

    4 結(jié)論

    1) CFG樁-土復(fù)合地基能有效減小高鐵運行引起的地面振動,具有良好的減振效果,且在距軌道較遠處的減振效果優(yōu)于軌道近處的減振效果。在現(xiàn)階段高鐵運營車速下(≤400 km/h),CFG樁-土復(fù)合地基可將距軌道中心≥25.0 m的地面振動降至滿足相關(guān)規(guī)范的振動限值要求(80.0 dB)。

    2) 在軌道近處,CFG樁-土復(fù)合地基的減振效果由車速和共振條件共同主導(dǎo)。當車速小于地基表層土Rayleigh 波速時,CFG 樁-土復(fù)合地基的減振效果隨與軌道中心距離增大而先增強后略有減弱。當車速接近或大于Rayleigh波速時,減振效果沿地表呈波動性變化。

    3) CFG 樁-土復(fù)合地基的減振效果隨樁徑d、樁土剛度比RPS增大而增強,隨著樁間距s減小而增強。當樁徑、樁間距或樁土剛度比RPS超過某一限值時(本文樁徑≥1/4 倍樁間距,樁間距≤4 倍樁徑,樁土剛度比RPS≥138.9),繼續(xù)增大樁徑、增大樁土剛度比RPS或減小樁間距均難以顯著提高CFG樁-土復(fù)合地基的減振效果。

    4) CFG樁-土復(fù)合地基提高了引發(fā)類共振的臨界速度,且該臨界速度受樁體剛度影響。當車速接近臨界速度時,CFG樁-土復(fù)合地基難以發(fā)揮較好的減振效果。因此,實際工程中,高速列車運行速度應(yīng)盡量避開CFG 樁-土復(fù)合地基的臨界速度。

    猜你喜歡
    振動效果模型
    一半模型
    振動的思考
    科學大眾(2023年17期)2023-10-26 07:39:14
    按摩效果確有理論依據(jù)
    重要模型『一線三等角』
    振動與頻率
    天天愛科學(2020年6期)2020-09-10 07:22:44
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    迅速制造慢門虛化效果
    抓住“瞬間性”效果
    中華詩詞(2018年11期)2018-03-26 06:41:34
    中立型Emden-Fowler微分方程的振動性
    3D打印中的模型分割與打包
    欧美成狂野欧美在线观看| 亚洲精品在线美女| 亚洲国产日韩欧美精品在线观看 | 熟女少妇亚洲综合色aaa.| 性色avwww在线观看| 国产一区二区激情短视频| 久久人人精品亚洲av| 亚洲色图av天堂| 国产亚洲精品一区二区www| 欧美日韩亚洲国产一区二区在线观看| 国产成人影院久久av| 欧美激情在线99| 国产一区二区在线av高清观看| 人人妻人人澡欧美一区二区| 国产精品久久久久久久电影 | 国产蜜桃级精品一区二区三区| 国产爱豆传媒在线观看| 精品福利观看| 亚洲国产欧美网| 99热精品在线国产| 高清在线国产一区| 亚洲人成伊人成综合网2020| 免费观看的影片在线观看| 日本免费a在线| 噜噜噜噜噜久久久久久91| 男女做爰动态图高潮gif福利片| 好男人电影高清在线观看| 黄色成人免费大全| 中文字幕av在线有码专区| 国内久久婷婷六月综合欲色啪| 夜夜夜夜夜久久久久| 51国产日韩欧美| 日韩欧美精品免费久久 | 天堂动漫精品| 久久99热这里只有精品18| 日本免费一区二区三区高清不卡| 色老头精品视频在线观看| 成人欧美大片| 窝窝影院91人妻| 制服人妻中文乱码| 欧美成人免费av一区二区三区| 亚洲成人久久爱视频| 午夜老司机福利剧场| 夜夜爽天天搞| 俄罗斯特黄特色一大片| 天堂√8在线中文| 一区福利在线观看| 精品久久久久久,| 亚洲专区中文字幕在线| av在线蜜桃| 免费av不卡在线播放| 久久国产精品影院| 在线国产一区二区在线| 一级毛片女人18水好多| xxx96com| 国产在线精品亚洲第一网站| 国产成人a区在线观看| 九色国产91popny在线| 国产色爽女视频免费观看| 欧美日韩亚洲国产一区二区在线观看| 婷婷精品国产亚洲av在线| 国产一区二区三区视频了| 国产精品亚洲美女久久久| 午夜日韩欧美国产| 亚洲在线自拍视频| 亚洲av成人精品一区久久| 色尼玛亚洲综合影院| 99久国产av精品| 日本a在线网址| 亚洲七黄色美女视频| 亚洲七黄色美女视频| av视频在线观看入口| 欧美激情在线99| 午夜福利在线观看吧| 十八禁人妻一区二区| 美女cb高潮喷水在线观看| 国产亚洲精品久久久com| 日本一二三区视频观看| 欧美成狂野欧美在线观看| 男人的好看免费观看在线视频| 亚洲精品影视一区二区三区av| 欧美成人一区二区免费高清观看| 午夜老司机福利剧场| 午夜免费激情av| 99在线视频只有这里精品首页| 日本一二三区视频观看| 99在线视频只有这里精品首页| 日韩有码中文字幕| 一级黄片播放器| 性欧美人与动物交配| 国产v大片淫在线免费观看| 色综合欧美亚洲国产小说| 午夜福利成人在线免费观看| 国产成人啪精品午夜网站| 国产精品电影一区二区三区| 亚洲国产日韩欧美精品在线观看 | 婷婷六月久久综合丁香| 久久精品国产亚洲av香蕉五月| 国产一区在线观看成人免费| 国产欧美日韩一区二区三| 久久精品国产清高在天天线| АⅤ资源中文在线天堂| 国产精品爽爽va在线观看网站| 国内毛片毛片毛片毛片毛片| 精品久久久久久久末码| 1024手机看黄色片| 国产中年淑女户外野战色| 欧美中文日本在线观看视频| 国产精品久久久久久久久免 | a级毛片a级免费在线| 舔av片在线| 日日干狠狠操夜夜爽| 久久天躁狠狠躁夜夜2o2o| 久久天躁狠狠躁夜夜2o2o| 亚洲国产色片| 日本撒尿小便嘘嘘汇集6| 欧美另类亚洲清纯唯美| 淫妇啪啪啪对白视频| 中文资源天堂在线| 在线a可以看的网站| 午夜福利高清视频| 国产真实伦视频高清在线观看 | 亚洲第一欧美日韩一区二区三区| 蜜桃久久精品国产亚洲av| 国产 一区 欧美 日韩| 啦啦啦韩国在线观看视频| 又黄又粗又硬又大视频| 亚洲精品在线美女| 亚洲精品一区av在线观看| 久久婷婷人人爽人人干人人爱| 亚洲精品美女久久久久99蜜臀| 精品久久久久久久毛片微露脸| 日本 欧美在线| 欧美日韩福利视频一区二区| 成熟少妇高潮喷水视频| 精品久久久久久久久久免费视频| 一级毛片高清免费大全| 又黄又爽又免费观看的视频| 久久精品91蜜桃| 亚洲精品一区av在线观看| 午夜a级毛片| 国产一区二区激情短视频| 免费人成视频x8x8入口观看| 日韩高清综合在线| 一进一出抽搐动态| 亚洲欧美精品综合久久99| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本在线视频免费播放| 国产成人a区在线观看| 好看av亚洲va欧美ⅴa在| 中文亚洲av片在线观看爽| 日本五十路高清| 色综合亚洲欧美另类图片| 综合色av麻豆| 免费看十八禁软件| 亚洲人成电影免费在线| 精品99又大又爽又粗少妇毛片 | 嫩草影院精品99| 在线免费观看的www视频| 日本精品一区二区三区蜜桃| 俄罗斯特黄特色一大片| 亚洲精品影视一区二区三区av| 免费观看的影片在线观看| 亚洲第一欧美日韩一区二区三区| 免费人成在线观看视频色| 99久久九九国产精品国产免费| 精品一区二区三区视频在线观看免费| 久久久久国产精品人妻aⅴ院| 成人高潮视频无遮挡免费网站| 欧美一区二区精品小视频在线| 日韩 欧美 亚洲 中文字幕| 久久久久久国产a免费观看| 女人十人毛片免费观看3o分钟| aaaaa片日本免费| 中文在线观看免费www的网站| 亚洲va日本ⅴa欧美va伊人久久| 精品国产超薄肉色丝袜足j| 国产成人系列免费观看| 九色成人免费人妻av| 母亲3免费完整高清在线观看| 精品无人区乱码1区二区| 亚洲天堂国产精品一区在线| 欧美高清成人免费视频www| 麻豆成人av在线观看| 午夜福利18| 国产高清视频在线观看网站| 午夜亚洲福利在线播放| 精品欧美国产一区二区三| 亚洲av不卡在线观看| 国产亚洲精品综合一区在线观看| 久久国产精品影院| 午夜影院日韩av| 美女高潮喷水抽搐中文字幕| 天天添夜夜摸| av女优亚洲男人天堂| 级片在线观看| 99久久精品一区二区三区| 日韩精品中文字幕看吧| 亚洲av成人不卡在线观看播放网| 欧美在线黄色| 一个人免费在线观看的高清视频| 欧美成狂野欧美在线观看| 日韩欧美在线二视频| 久久精品国产综合久久久| 1000部很黄的大片| 3wmmmm亚洲av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 一边摸一边抽搐一进一小说| 免费看日本二区| 日日摸夜夜添夜夜添小说| 亚洲精品一卡2卡三卡4卡5卡| 婷婷精品国产亚洲av在线| 欧美黄色片欧美黄色片| 一夜夜www| 无人区码免费观看不卡| 两个人看的免费小视频| 亚洲最大成人中文| 国产三级中文精品| avwww免费| 天美传媒精品一区二区| 黄色片一级片一级黄色片| 国产成人av激情在线播放| 亚洲色图av天堂| 国产野战对白在线观看| 国产一区二区激情短视频| 国产精华一区二区三区| 色尼玛亚洲综合影院| 国产精品亚洲美女久久久| 亚洲精品在线观看二区| 亚洲内射少妇av| 亚洲av熟女| 最新在线观看一区二区三区| 美女被艹到高潮喷水动态| 亚洲人成网站在线播放欧美日韩| 色在线成人网| 国产成人福利小说| 俺也久久电影网| 欧美国产日韩亚洲一区| 99国产精品一区二区三区| 久久久久亚洲av毛片大全| 岛国视频午夜一区免费看| 日本 欧美在线| 51国产日韩欧美| 欧美一级a爱片免费观看看| 欧美最新免费一区二区三区 | 午夜精品一区二区三区免费看| 免费观看精品视频网站| 国产视频一区二区在线看| 亚洲精品成人久久久久久| 极品教师在线免费播放| 国产视频内射| 窝窝影院91人妻| 国产高潮美女av| 一进一出好大好爽视频| 久久6这里有精品| 99久久综合精品五月天人人| 久久久久免费精品人妻一区二区| 免费在线观看影片大全网站| 日本精品一区二区三区蜜桃| a在线观看视频网站| 成人18禁在线播放| 精品福利观看| 国产精品av视频在线免费观看| 狠狠狠狠99中文字幕| 99久久精品国产亚洲精品| 国产色爽女视频免费观看| 一个人看的www免费观看视频| 国产av在哪里看| 丰满乱子伦码专区| 日韩有码中文字幕| 日日摸夜夜添夜夜添小说| 又紧又爽又黄一区二区| 久久99热这里只有精品18| 日韩欧美免费精品| www国产在线视频色| 中文字幕av在线有码专区| 免费在线观看影片大全网站| 亚洲人成电影免费在线| 国产在视频线在精品| 黄色成人免费大全| 国产免费av片在线观看野外av| 国产精品久久久久久人妻精品电影| 成年免费大片在线观看| 欧美日本视频| 好男人电影高清在线观看| 久久久久久国产a免费观看| 热99re8久久精品国产| 一区二区三区国产精品乱码| 精品人妻偷拍中文字幕| 国产一区二区激情短视频| 女警被强在线播放| 国产精品野战在线观看| 日韩免费av在线播放| 国产av不卡久久| 真人一进一出gif抽搐免费| 观看免费一级毛片| 亚洲精品在线美女| 国产成年人精品一区二区| 十八禁人妻一区二区| 九色国产91popny在线| 欧美三级亚洲精品| 亚洲精品色激情综合| av女优亚洲男人天堂| 精品久久久久久久人妻蜜臀av| 日韩欧美国产一区二区入口| 国产一区二区亚洲精品在线观看| 成人性生交大片免费视频hd| 国产激情偷乱视频一区二区| 又粗又爽又猛毛片免费看| 国产精品精品国产色婷婷| 免费在线观看亚洲国产| 欧美日韩国产亚洲二区| 亚洲成人中文字幕在线播放| 欧美激情在线99| 99久久成人亚洲精品观看| 99视频精品全部免费 在线| 岛国视频午夜一区免费看| 真人一进一出gif抽搐免费| 成人一区二区视频在线观看| 日本五十路高清| 丝袜美腿在线中文| xxx96com| 亚洲专区中文字幕在线| 乱人视频在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品99久久久久久久久| 欧美av亚洲av综合av国产av| 亚洲精华国产精华精| 99热精品在线国产| 国产三级在线视频| 性色av乱码一区二区三区2| 午夜福利成人在线免费观看| 波野结衣二区三区在线 | 亚洲精华国产精华精| 久久国产精品人妻蜜桃| 国产高潮美女av| 噜噜噜噜噜久久久久久91| 亚洲avbb在线观看| 亚洲中文字幕日韩| 国产精品日韩av在线免费观看| e午夜精品久久久久久久| 白带黄色成豆腐渣| 91麻豆精品激情在线观看国产| 热99re8久久精品国产| 日本熟妇午夜| 他把我摸到了高潮在线观看| 亚洲精华国产精华精| 一级毛片女人18水好多| 大型黄色视频在线免费观看| 超碰av人人做人人爽久久 | 欧美绝顶高潮抽搐喷水| 国内揄拍国产精品人妻在线| 国产亚洲精品久久久com| 高潮久久久久久久久久久不卡| 日韩免费av在线播放| 免费人成视频x8x8入口观看| 黑人欧美特级aaaaaa片| 97碰自拍视频| 性色av乱码一区二区三区2| 亚洲国产精品久久男人天堂| 久久精品91蜜桃| 99国产综合亚洲精品| 欧美日韩亚洲国产一区二区在线观看| 黄色日韩在线| 国产成人福利小说| 成人鲁丝片一二三区免费| av天堂中文字幕网| 亚洲精品粉嫩美女一区| 一级黄片播放器| 少妇熟女aⅴ在线视频| 97超视频在线观看视频| 亚洲人成电影免费在线| 亚洲欧美日韩东京热| 午夜免费观看网址| 很黄的视频免费| 一级黄色大片毛片| 久久香蕉国产精品| av片东京热男人的天堂| 伊人久久精品亚洲午夜| 丰满的人妻完整版| 成人av在线播放网站| 欧美日韩国产亚洲二区| 97超视频在线观看视频| 最新在线观看一区二区三区| 国产精品1区2区在线观看.| 国产精品一区二区三区四区久久| 亚洲人与动物交配视频| 免费在线观看日本一区| 国产高清videossex| 一个人观看的视频www高清免费观看| 婷婷精品国产亚洲av在线| 免费av观看视频| 国产精品永久免费网站| 99热这里只有精品一区| 亚洲国产欧美网| 亚洲中文字幕日韩| 长腿黑丝高跟| 成人鲁丝片一二三区免费| 日韩中文字幕欧美一区二区| 欧美性猛交黑人性爽| www.熟女人妻精品国产| 此物有八面人人有两片| 69av精品久久久久久| 免费人成在线观看视频色| 久久精品夜夜夜夜夜久久蜜豆| 久9热在线精品视频| 国内精品一区二区在线观看| 波多野结衣高清作品| www国产在线视频色| 国产午夜福利久久久久久| 日韩欧美免费精品| 两个人看的免费小视频| 国产男靠女视频免费网站| 成人高潮视频无遮挡免费网站| 国产精品一及| 欧美黄色淫秽网站| 色哟哟哟哟哟哟| 99国产精品一区二区蜜桃av| 脱女人内裤的视频| 午夜福利欧美成人| 在线观看日韩欧美| 成年版毛片免费区| 亚洲性夜色夜夜综合| 熟女少妇亚洲综合色aaa.| 不卡一级毛片| 亚洲aⅴ乱码一区二区在线播放| 午夜视频国产福利| 午夜免费观看网址| 国产中年淑女户外野战色| 99久久九九国产精品国产免费| 国产成年人精品一区二区| xxxwww97欧美| 国产99白浆流出| 成人午夜高清在线视频| 丁香六月欧美| 中文字幕人成人乱码亚洲影| 黄色片一级片一级黄色片| 小说图片视频综合网站| 日韩欧美免费精品| 人妻丰满熟妇av一区二区三区| 美女大奶头视频| 成人av在线播放网站| 精品不卡国产一区二区三区| 少妇裸体淫交视频免费看高清| 国产色婷婷99| 国产真人三级小视频在线观看| 午夜精品久久久久久毛片777| 亚洲五月天丁香| 天天躁日日操中文字幕| 真人做人爱边吃奶动态| 久久99热这里只有精品18| 亚洲欧美日韩高清专用| h日本视频在线播放| 亚洲熟妇熟女久久| 99热6这里只有精品| 国产99白浆流出| 久久久久久久精品吃奶| 少妇的逼水好多| 天堂√8在线中文| 午夜福利成人在线免费观看| 久久欧美精品欧美久久欧美| 最近最新免费中文字幕在线| 中文字幕av在线有码专区| 国产成年人精品一区二区| 精品99又大又爽又粗少妇毛片 | 久久九九热精品免费| 国产精品一区二区三区四区免费观看 | 青草久久国产| 国产伦一二天堂av在线观看| 长腿黑丝高跟| 成人18禁在线播放| 日韩成人在线观看一区二区三区| 欧美日韩福利视频一区二区| 欧美最新免费一区二区三区 | 他把我摸到了高潮在线观看| 精品日产1卡2卡| 国产精品久久久久久亚洲av鲁大| 男插女下体视频免费在线播放| 欧美激情久久久久久爽电影| 国产精品精品国产色婷婷| 啦啦啦观看免费观看视频高清| 香蕉av资源在线| 成人18禁在线播放| 午夜福利欧美成人| 国产97色在线日韩免费| 亚洲成人久久性| 禁无遮挡网站| 国产精品国产高清国产av| 国产精品精品国产色婷婷| 麻豆一二三区av精品| 久久精品影院6| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 怎么达到女性高潮| 美女大奶头视频| 亚洲成a人片在线一区二区| 国产爱豆传媒在线观看| 亚洲av免费高清在线观看| 黄色成人免费大全| 国产精品久久久久久人妻精品电影| 熟妇人妻久久中文字幕3abv| 有码 亚洲区| 国内精品久久久久精免费| 亚洲片人在线观看| 欧美日韩福利视频一区二区| 天天一区二区日本电影三级| 免费av毛片视频| 亚洲成av人片免费观看| 日韩欧美在线乱码| 午夜福利欧美成人| 黄片大片在线免费观看| 脱女人内裤的视频| 亚洲真实伦在线观看| 在线观看美女被高潮喷水网站 | 精华霜和精华液先用哪个| 午夜福利在线观看吧| 国产在线精品亚洲第一网站| 91久久精品电影网| 18禁在线播放成人免费| 人妻久久中文字幕网| 国产极品精品免费视频能看的| 99在线人妻在线中文字幕| x7x7x7水蜜桃| 国产乱人伦免费视频| 国产av不卡久久| 久久久国产成人精品二区| av天堂中文字幕网| 亚洲国产欧美人成| 婷婷丁香在线五月| 亚洲中文字幕一区二区三区有码在线看| 有码 亚洲区| 成年女人毛片免费观看观看9| 小说图片视频综合网站| 午夜福利免费观看在线| 黑人欧美特级aaaaaa片| 97超视频在线观看视频| 日韩欧美国产在线观看| 久久99热这里只有精品18| 国产精品永久免费网站| 欧美bdsm另类| 校园春色视频在线观看| 亚洲自拍偷在线| 色av中文字幕| 国产真人三级小视频在线观看| www国产在线视频色| 五月伊人婷婷丁香| 欧美成人免费av一区二区三区| www日本黄色视频网| 亚洲精品456在线播放app | 国产高清videossex| 色综合亚洲欧美另类图片| 尤物成人国产欧美一区二区三区| 日本 av在线| 天堂影院成人在线观看| 国产私拍福利视频在线观看| 天堂影院成人在线观看| 久久精品国产清高在天天线| www.熟女人妻精品国产| 婷婷六月久久综合丁香| 法律面前人人平等表现在哪些方面| 啦啦啦免费观看视频1| 成人亚洲精品av一区二区| 一a级毛片在线观看| 日本三级黄在线观看| 九九久久精品国产亚洲av麻豆| 国产高清videossex| 国产精品,欧美在线| 久久人妻av系列| 国产成年人精品一区二区| 欧美一级a爱片免费观看看| 美女高潮喷水抽搐中文字幕| 制服人妻中文乱码| 香蕉av资源在线| 日本一本二区三区精品| 哪里可以看免费的av片| 窝窝影院91人妻| 久久久久久大精品| 韩国av一区二区三区四区| 日本免费a在线| 婷婷亚洲欧美| 亚洲精品影视一区二区三区av| 黄色成人免费大全| 亚洲五月婷婷丁香| 日韩欧美 国产精品| 在线观看日韩欧美| 又粗又爽又猛毛片免费看| 在线看三级毛片| 日本与韩国留学比较| 99久国产av精品| 久久这里只有精品中国| svipshipincom国产片| 丰满的人妻完整版| 国产成人a区在线观看| 亚洲熟妇熟女久久| 午夜日韩欧美国产| 一夜夜www| 变态另类成人亚洲欧美熟女| www.www免费av| 色哟哟哟哟哟哟| 非洲黑人性xxxx精品又粗又长| 亚洲欧美日韩无卡精品| 99久久综合精品五月天人人| 久久婷婷人人爽人人干人人爱| 黑人欧美特级aaaaaa片| 女生性感内裤真人,穿戴方法视频| 亚洲午夜理论影院| 69人妻影院| 亚洲中文日韩欧美视频| 一a级毛片在线观看| 最新在线观看一区二区三区| 午夜福利在线观看免费完整高清在 | 亚洲专区国产一区二区| 久久香蕉国产精品| www.熟女人妻精品国产| 97人妻精品一区二区三区麻豆| 99热6这里只有精品|