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

    軌道隨機不平順的譜表示-隨機函數(shù)模擬及應用

    2021-07-09 01:12:20劉付山曾志平王衛(wèi)東
    鐵道學報 2021年5期
    關(guān)鍵詞:概率密度平順級數(shù)

    劉付山,曾志平,2,王衛(wèi)東,2

    (1.中南大學 土木工程學院,湖南 長沙 410075;2.中南大學 重載鐵路工程結(jié)構(gòu)教育部重點實驗室,湖南 長沙 410075)

    軌道隨機不平順是由不同波長、不同幅值、不同相位的不同不平順波疊加而成的隨線路里程變化的復雜隨機過程[1]。軌道隨機不平順常通過線路實地測量獲得,并采用功率譜密度函數(shù)表示其統(tǒng)計特征,我國及英國、日本、德國等都測定了各自的軌道不平順譜密度和相關(guān)函數(shù)。然而,在非線性車輛-軌道耦合系統(tǒng)中,功率譜密度函數(shù)無法直接輸入動力系統(tǒng),最有效的辦法是以軌道隨機不平順的時域樣本作為輸入。翟婉明等[2]提出一種根據(jù)軌道隨機不平順功率譜求出頻譜的幅值和隨機相位,再通過傅里葉逆變換得到軌道隨機不平順時域樣本的數(shù)值模擬方法。除此以外,三角級數(shù)法也是采用較多的一種軌道隨機不平順時域樣本數(shù)值模擬方法[3]。眾所周知,隨機分析計算效率隨著所需樣本數(shù)增加而增加,為了更加高效的研究軌道隨機不平順激勵下的車輛-軌道耦合系統(tǒng)隨機振動特性,需要用更少的隨機樣本來反映軌道隨機不平順的統(tǒng)計特性。無論是常用的三角級數(shù)法還是傅里葉逆變換法,在模擬軌道隨機不平順時都需要大量的隨機變量,而分析所需的隨機樣本數(shù)與隨機變量數(shù)正相關(guān)。因而用更少的隨機變量來描述軌道隨機不平順,可極大提高的車輛-軌道耦合系統(tǒng)隨機振動分析的效率及難度。

    文獻[4]中,通過隨機諧和函數(shù)法生成軌道隨機不平順時域樣本,并通過頻率與相位分別隨機以減少所需的隨機變量數(shù),其中隨機頻率數(shù)為50個,考慮軌道高低、方向、水平、軌距不平順情況下,各種不平順相互獨立,隨機頻率和隨機相位共計需400個隨機變量,所需隨機變量依然很多,且單個隨機樣本中頻率分量減少必然會導致模擬的軌道不平順與實際偏差增大,并不能很好地反映不平順的時域特征,該方法可用于線性車輛-軌道耦合系統(tǒng)中軌道隨機不平順的模擬,但不適用于非線性車輛-軌道耦合動力學系統(tǒng)。

    本文基于譜表示方法,將譜表示方法中的標準正交隨機變量表達為基本隨機變量的正交函數(shù)形式,從而實現(xiàn)用較少隨機變量來描述軌道隨機不平順的目的。減少隨機變量數(shù)的同時,保留更多的頻率及相位分量,可更準確反映軌道隨機不平順統(tǒng)計特性以及時域特征,提高非線性車輛-軌道耦合系統(tǒng)隨機振動計算效率。

    1 隨機過程的譜表示-隨機函數(shù)模擬

    一維、單變量、功率譜密度函數(shù)為S(ω)的隨機過程η(t)可以表示為[5]

    (1)

    式中:Ut(ω)和Vt(ω)為隨機過程η(t)的譜過程,其增量dUt(ω)和dVt(ω)滿足隨機過程譜表示的基本條件,即

    其中,E[·]表示數(shù)學期望;ω、ω′為隨機過程中的頻率;S(·)為雙邊功率譜密度函數(shù)。

    對式(1)進行離散,并引入一組標準正交隨機變量{Xk,Yk},當頻率間隔Δω足夠小時,式(1)可表述為

    (2)

    式中:dUt(ω)≈ΔUt(ω)=σkXk; dVt(ω)≈ΔVt(ω)=σkYk;σk=[2S(ωk)Δω]1/2,S為軌道不平順功率譜密度函數(shù);Nω為頻率離散點數(shù);頻率分量均勻離散情況下,ωk=ωl+k·Δω,Δω=(ωu-ωl)/Nω,ωu、ωl為頻率離散點上、下限值。

    (3)

    式中:cas(x)=cos(x)+sin(x)為Hartley正交基函數(shù);隨機變量Θ1、Θ2相互獨立且服從區(qū)間[0,2π)上均勻分布。直接生成的標準正交變量并不能直接代入式(2),還需要構(gòu)造隨機映射將標準正交變量映射到標準正交隨機變量{Xk,Yk},由此將描述隨機過程的隨機變量數(shù)降低為2。

    2 軌道隨機不平順的譜表示-隨機函數(shù)模擬方法

    以TB/T 3352—2014《高速鐵路無砟軌道不平順譜》[6]為例,軌道不平順譜包括高低、方向、軌距和水平四種軌道不平順譜,計算不平順譜時,將實測得到的軌道不平順按1 024 m劃分單元,剔除異常值和趨勢項,加窗處理后進行傅里葉變換(FFT)計算,對所有單元的軌道不平順譜進行擬合,得到軌道不平順功率譜密度函數(shù)。我國的軌道不平順譜擬合公式為

    (4)

    式中:f為空間頻率;A、k為擬合系數(shù)。

    統(tǒng)計誤差可表示為

    (5)

    需要指出的是,軌道不平順功率譜密度函數(shù)是在大量實測數(shù)據(jù)基礎(chǔ)上,按一定單元長度組成軌道不平順樣本集,由樣本集計算得到的均值譜。單個軌道不平順分段的隨機誤差達到了100%,隨著分段次數(shù)的增加隨機誤差逐漸減小,如果分段100次則隨機誤差可以降到10%[6]。因而采用單個軌道隨機不平順樣本或單元計算得到的不平順功率譜與現(xiàn)有軌道不平順功率譜對比,并不能說明數(shù)值模擬方法的準確性,需采用多個不平順樣本或分段所組成的樣本集不平順功率譜與現(xiàn)有軌道不平順功率譜對比驗證。

    對于多樣本模擬,常采用蒙特卡羅法選擇隨機變量,但其計算效率偏低。擬蒙特卡羅法是通過數(shù)論方法產(chǎn)生高度均勻化的擬隨機數(shù)序列代替蒙特卡羅隨機數(shù)模擬,以減少隨機樣本數(shù)從而提高計算效率,但當隨機變量較多,即高維度情況下,數(shù)論選點法選取的隨機數(shù)列不可避免的會產(chǎn)生叢聚現(xiàn)象[7]。采用譜表示-隨機函數(shù)法模擬生成軌道不平順時域樣本,可減少所需的隨機變量數(shù),避免高維數(shù)論選點中叢聚現(xiàn)象對計算結(jié)果的影響。

    綜上所述,采用譜表示-隨機函數(shù)法生成軌道隨機不平順樣本集過程可表示為

    Step1根據(jù)生成的軌道隨機不平順類型確定隨機變量數(shù)N,當需同時生成軌道高低、方向、軌距、水平不平順時,由于各類型不平順間相互獨立,因而需要8個隨機變量(各類型不平順均為2個)。

    Step2根據(jù)所需隨機變量數(shù),采用擬蒙特卡羅法生成隨機變量的低偏差序列Θ,常用的低偏差序列有Corput序列、Halton序列、Sobol序列等[7]。本文采用Sobol序列生成隨機變量[5],采用Matlab直接計算Sobol的命令。

    Θ={θ1j,θ2j,…,θNj}j=1,2,…,Npt

    式中:N為隨機變量數(shù);Npt為樣本集合的樣本數(shù)。此時生成的隨機變量序列服從[0,1)上均勻分布, 需將其變換為[0,2π)上均勻分布隨機變量Θ′

    Θ′=2πΘ

    (6)

    高低不平順ηV(t)為

    方向不平順ηA(t)為

    水平不平順ηC(t)為

    軌距不平順ηG(t)為

    式中:SV、SA、SC、SG分別為高低、方向、水平、軌距不平順功率譜密度:ωk為軌道隨機不平順中不同的頻率分量,其表達式為

    (7)

    采用譜表示-隨機函數(shù)模擬方法生成樣本數(shù)為300的軌道隨機不平順樣本集,其中軌道不平順功率譜計算方法中單元長1 024 m,采樣間隔0.25 m,單個單元共計4 096個點,因而模擬軌道隨機不平順時頻率離散點數(shù)Nω取4 096。

    對生成的軌道不平順樣本集中各樣本,參考文獻[8]中高速鐵路無砟軌道不平順譜的生成方法及軌道不平順譜頻率分布規(guī)律檢驗方法,采用周期圖法計算得到各樣本的不平順譜值,并求得樣本集不平順譜均值,對各頻點不平順譜值Sk,j按下式進行變換

    (8)

    采用Kolmogorov-Smirnov檢驗方法,對λk,j進行顯著性水平α=0.05、自由度為2的χ2分布檢驗。軌道不平順單樣本情況、單樣本功率譜、樣本均值、樣本標準差、部分頻點功率譜概率分布CDF等評價指標以及樣本集總體功率譜PSD與目標譜對比情況如圖1~圖4所示。

    圖1 譜表示-隨機函數(shù)法高低不平順樣本評價指標對比

    圖2 譜表示-隨機函數(shù)法方向不平順樣本評價指標對比

    圖3 譜表示-隨機函數(shù)法水平不平順樣本評價指標對比

    圖4 譜表示-隨機函數(shù)法軌距不平順樣本評價指標對比

    從圖1~圖4中可以看出,計算得到的軌道高低、方向、水平、軌距不平順單個樣本功率譜密度值與軌道不平順功率譜存在較大的誤差,而軌道不平順樣本集計算得到的功率譜均值與目標譜吻合較好,其中各種不平順樣本總體均值與目標譜最大相差分別為0.18、0.10、0.07、0.08 mm,標準差最大偏差分別為4.68%、5.77%、6.97%、6.99%,樣本集功率譜最大偏差分別為14.54%、10.23%、10.07%、9.87%,各不平順98.00%、96.33%、98.00%、98.33%的頻點功率譜概率分布滿足自由度為2的χ2分布??梢钥闯鰳颖炯髦笜伺c目標譜吻合良好,且僅需8個隨機變量便可同時模擬生成軌道高低、方向、水平、軌道不平順時域樣本,同時還具有大量頻率分量。

    3 軌道隨機不平順的數(shù)值模擬方法對比分析

    目前在工程中應用較多的軌道隨機不平順數(shù)值模擬方法有三角級數(shù)和傅里葉逆變換法。三角級數(shù)法模擬軌道隨機不平順可表示為[2]

    (9)

    式中:ak為相互獨立的均值為0、標準差為σk的高斯隨機變數(shù);φk為相互獨立的隨機相位。

    傅里葉逆變換法模擬軌道隨機不平順可表示為[1]

    η=ifft(X(k))

    (10)

    式中:X(k)為時間序列,其為復序列,實部關(guān)于Nr/2偶對稱,虛部關(guān)于Nr/2奇對稱。設(shè)模擬總時間為Ts,則時域采樣點數(shù)為Ts/ΔT,一般需要在末尾添0以保證采樣點數(shù)為2的整數(shù)次冪,即為Nr。

    X(k)表達式為

    k=0,1,…,Nr-1

    (11)

    式中:ξn為獨立隨機相位序列,ξn=cosφn+isinφn=exp(iφn),φn服從[0,2π]上的均勻分布;Δf=1/(NrΔT)。

    三角級數(shù)法中,ak為服從高斯分布的隨機變量,φk為服從均勻分布的隨機變量,因而三角級數(shù)模擬軌道隨機不平順時所需隨機變量數(shù)為2Nω,且存在高斯分布與均勻分布的兩種類型隨機變量。傅里葉逆變化法中,僅有φn服從正態(tài)分布且相互獨立,而所需隨機變量數(shù)與時域采樣點數(shù)相關(guān),一般為2的整數(shù)次冪,當采樣間隔0.25 m,模擬1 km軌道隨機不平順時所需隨機變量數(shù)為4 096個。三種方法模擬軌道隨機不平順時域樣本時,不平順頻率點計算方法基本一致,采用相同頻率數(shù)時,具有相同的不平順波長成分,因而文中僅對不平順時域樣本集的幅值進行對比分析。

    采用擬蒙特卡羅法生成低偏差隨機向量序列,通過三角級數(shù)法、傅里葉逆變換法模擬得到300個軌道隨機不平順時域樣本。計算樣本集均值、標準差、功率譜,并檢驗樣本集各頻點軌道譜分布,各項對比情況如圖5和圖6所示。各種方法生成的軌道不平順樣本集不平順幅值范圍情況見表1。

    圖5 三角級數(shù)法高低不平順樣本評價指標對比

    圖6 傅里葉逆變換法高低不平順樣本評價指標對比

    表1 樣本集不平順幅值范圍 mm

    可以看出,除了傅里葉逆變換法計算得到高低不平順幅值范圍偏小外,三種方法計算得到的軌道隨機不平順幅值范圍大體相當,最大相差5.51%。三角級數(shù)法計算得到的軌道不平順單個樣本功率譜與目標功率譜存在較大誤差,各方向軌道不平順樣本集與目標譜均值最大相差分別為0.22、0.09、0.11、0.08 mm,標準差最大偏差分別為5.14%、6.78%、4.95%、7.25%,功率譜密度最大偏差分別為15.82%、13.93%、10.78%、9.65%,各不平順97.33%、96.67%、97.33%、98.67%的頻點功率譜概率分布均滿足自由度為2的χ2分布。傅里葉逆變換法計算得到的軌道隨機不平順,無論是單樣本還是樣本集,其功率譜值均能很好滿足目標功率譜,最大偏差僅分別為0.30%、0.24%、0.35%、0.13%,但各頻點功率譜概率分布均不滿足文獻[8]中關(guān)于軌道不平順譜分布的論述。除此以外,傅里葉逆變換法各方向軌道不平順樣本集與目標譜均值最大相差分別為0.28、0.11、0.13、0.07 mm,標準差最大偏差分別為11.38%、8.53%、8.73%、6.57%,與其他兩種方法模擬結(jié)果基本一致。綜上可以看出,譜表示-隨機函數(shù)法在模擬軌道不平順時域樣本時其幅值范圍基本一致。

    4 非線性車輛-軌道-橋梁耦合系統(tǒng)隨機振動分析

    采用譜表示-隨機函數(shù)方法僅需8個隨機變量就可以描述四種軌道不平順,這為非線性車輛-軌道-橋梁耦合系統(tǒng)隨機振動提供了極大的方便。文獻[9]中指出,由于涉及輪軌空間動態(tài)接觸幾何關(guān)系、輪軌切向蠕滑力和輪軌法向力等多種復雜非線性求解過程,車輛-軌道耦合系統(tǒng)最有效的隨機振動研究方法是將軌道不平順譜轉(zhuǎn)換為時域樣本,采用一段足夠長的時域樣本來代替整個隨機過程,即采用蒙特卡洛法選取足夠長的時域樣本。

    分別采用三角級數(shù)法、傅里葉逆變換法、譜表示-隨機函數(shù)法模擬得到的軌道隨機不平順樣本,利用三維非線性車輛-軌道-橋梁耦合動力學模型,計算得到系統(tǒng)動力響應,對不同距離上的動力響應結(jié)果進行統(tǒng)計分析,以輪重減載率、脫軌系數(shù)為例,樣本計算結(jié)果截口概率分布PDF見圖7。

    圖7 不同計算法的響應截口概率分布

    從圖7中可以看出,在300個樣本集的情況下,三種方法計算得到的動力響應其截口概率分布基本相同,都具有較大的離散性,速度300 km/h、不平順波長2~120 m時,輪重減載率變化范圍為±0.25,脫軌系數(shù)變化范圍為±0.12。由于文中采用中國高速鐵路無砟軌道譜,平順性較其他軌道譜好,因而輪重減載率與脫軌系數(shù)的計算結(jié)果偏小。

    對于車輛-軌道耦合系統(tǒng)動力響應,采用三種軌道不平順譜數(shù)值模擬方法計算得到結(jié)果基本一致,但所需的隨機變量數(shù)卻存在巨大差異。隨機變量數(shù)的減少方便了樣本空間概率賦予,同時也可在一定程度上減少所需的樣本數(shù),為采用廣義概率密度演化理論對車輛運行安全性進行評價提供了方便。

    對于非線性車輛-軌道-橋梁耦合系統(tǒng)的任一動力響應ξ(t),車輛運行過程中不考慮其他系統(tǒng)隨機的情況下,其隨機性來源于軌道隨機不平順中隨機參數(shù)Θ的隨機性,該系統(tǒng)為一個概率保守隨機系統(tǒng)。根據(jù)隨機系統(tǒng)的概率守恒原理,其應滿足廣義概率密度演化方程[10]

    (12)

    研究中常通過與大量蒙特卡洛法(MCM)計算結(jié)果對比,對廣義概率密度演化方法(PDEM)計算結(jié)果進行驗證[4,12]。采用三角級數(shù)法計算模擬軌道隨機不平順樣本,結(jié)合MCM法計算得到3 000個樣本動力響應統(tǒng)計情況,并采用譜表述隨機函數(shù)法,結(jié)合廣義概率密度演化理論,計算得到300個樣本下的動力響應概率密度演化結(jié)果,兩種方法車輛運行距離中間點的截口概率密度對比情況見圖8,由于三角級數(shù)法所需的隨機變量數(shù)較多,因而在3 000個樣本情況下,動力響應結(jié)果依然存在一定離散性,但兩種方法計算得到的截口概率密度基本一致,因而采用譜表示-隨機函數(shù)法模擬軌道隨機不平順,結(jié)合廣義概率密度演化理論可極大提高非線性車輛-軌道-橋梁耦合系統(tǒng)隨機振動分析效率,相較于采用三角級數(shù)等方法模擬軌道隨機不平順并結(jié)合MCM法進行車輛-軌道-橋梁耦合系統(tǒng)隨機振動分析,其計算效率至少可提高一個量級。參考“3σ”準則,取E±3σ對應值為動力響應指標最值。車輛運行過程中輪重減載率及脫軌系數(shù)最大值、等概率密度演化曲線見圖9。由此,譜表示-隨機函數(shù)法軌道不平順數(shù)值模擬方法結(jié)合廣義概率密度演化理論可對非線性車輛-軌道-橋梁耦合系統(tǒng)車輛運行安全性進行評價。

    圖8 MCM與PDEM截口概率密度對比

    圖9 等概率密度演化曲線

    5 結(jié)論

    軌道不平順譜是描述軌道隨機不平順最常用的統(tǒng)計函數(shù),本文以軌道隨機不平順時域樣本數(shù)值模擬方法為研究對象,采用譜表示-隨機函數(shù)法分別由2個隨機變量模擬生成軌道高低、水平、方向、軌距不平順集。

    對比樣本集與目標譜均值、標準差、軌道譜以及各頻點軌道譜概率分布,驗證了譜表示-隨機函數(shù)法模擬軌道隨機不平順時域樣本準確性;與三角級數(shù)法、傅里葉逆變化法相比,模擬生成軌道不平順幅值范圍基本一致,誤差在5.51%以內(nèi),其中傅里葉逆變化法生成的不平順樣本功率譜能很好地吻合目標譜,但各頻點功率譜分布并不滿足與文獻[8]所述頻點功率譜分布;采用譜表示-隨機函數(shù)法模擬生成軌道隨機不平順時域樣本集,結(jié)合廣義概率密度演化理論,可用相對較少樣本計算得到車輛運行過程中概率演化情況,相較于采用三角級數(shù)法模擬軌道隨機不平順,并結(jié)合MCM法進行隨機振動分析,計算效率至少提高了1個量級,為非線性車輛-軌道-基礎(chǔ)結(jié)構(gòu)耦合系統(tǒng)隨機振動分析中軌道隨機不平順的模擬提供了新思路。

    猜你喜歡
    概率密度平順級數(shù)
    基于Simulink的汽車行駛平順性研究
    連續(xù)型隨機變量函數(shù)的概率密度公式
    平順植保站:開展粟灰螟防治
    Dirichlet級數(shù)及其Dirichlet-Hadamard乘積的增長性
    幾個常數(shù)項級數(shù)的和
    p級數(shù)求和的兩種方法
    Hunt過程在Girsanov變換下的轉(zhuǎn)移概率密度的表示公式
    隨機變量線性組合的分布的一個算法
    隨機結(jié)構(gòu)-TMD優(yōu)化設(shè)計與概率密度演化研究
    Dirichlet級數(shù)的Dirichlet-Hadamard乘積
    少妇高潮的动态图| 亚洲一级一片aⅴ在线观看| 最新中文字幕久久久久| 国产av国产精品国产| 少妇人妻精品综合一区二区| 成人欧美大片| 久久久久久九九精品二区国产| av免费在线看不卡| 在现免费观看毛片| 99热这里只有是精品50| 国产免费福利视频在线观看| 日韩精品有码人妻一区| 午夜福利在线观看吧| 国产精品久久久久久久久免| 欧美不卡视频在线免费观看| 精品国产三级普通话版| 欧美日韩国产mv在线观看视频 | 听说在线观看完整版免费高清| 一级毛片我不卡| 九九爱精品视频在线观看| 内射极品少妇av片p| 免费高清在线观看视频在线观看| 国产精品爽爽va在线观看网站| 国产片特级美女逼逼视频| 国产综合懂色| 日本黄大片高清| 国产精品国产三级国产av玫瑰| 美女内射精品一级片tv| 最近最新中文字幕大全电影3| 一个人看的www免费观看视频| 人人妻人人澡欧美一区二区| 精品人妻偷拍中文字幕| 国产精品精品国产色婷婷| 高清av免费在线| 一边亲一边摸免费视频| 好男人视频免费观看在线| 国产亚洲精品久久久com| 啦啦啦中文免费视频观看日本| 两个人视频免费观看高清| 乱人视频在线观看| av在线天堂中文字幕| 日本av手机在线免费观看| 91在线精品国自产拍蜜月| 国产精品久久久久久av不卡| 高清欧美精品videossex| 老司机影院成人| 国产中年淑女户外野战色| 成人二区视频| 午夜精品一区二区三区免费看| 人妻一区二区av| 日本免费在线观看一区| 99久久精品热视频| 男女视频在线观看网站免费| 尤物成人国产欧美一区二区三区| av福利片在线观看| 又粗又硬又长又爽又黄的视频| 2021少妇久久久久久久久久久| 美女大奶头视频| 国产 亚洲一区二区三区 | 少妇的逼好多水| 国产精品嫩草影院av在线观看| 国产一区有黄有色的免费视频 | 少妇裸体淫交视频免费看高清| 建设人人有责人人尽责人人享有的 | 中文欧美无线码| 国产视频首页在线观看| 秋霞在线观看毛片| www.av在线官网国产| 精华霜和精华液先用哪个| 一夜夜www| 色视频www国产| 日韩在线高清观看一区二区三区| 国产 一区精品| 免费大片18禁| 亚洲精品乱码久久久v下载方式| 日韩,欧美,国产一区二区三区| 日本色播在线视频| 一二三四中文在线观看免费高清| 十八禁国产超污无遮挡网站| 午夜福利高清视频| 久久久精品94久久精品| 久热久热在线精品观看| 熟女电影av网| 国产亚洲午夜精品一区二区久久 | 美女被艹到高潮喷水动态| 男女啪啪激烈高潮av片| 欧美日韩视频高清一区二区三区二| 韩国高清视频一区二区三区| 伦理电影大哥的女人| 五月伊人婷婷丁香| 高清欧美精品videossex| 精品久久久噜噜| 日本黄大片高清| 一级a做视频免费观看| 人妻少妇偷人精品九色| 国产成人一区二区在线| 26uuu在线亚洲综合色| 色播亚洲综合网| 免费黄网站久久成人精品| 亚洲av二区三区四区| 99久久中文字幕三级久久日本| 最近手机中文字幕大全| 午夜免费男女啪啪视频观看| 男女下面进入的视频免费午夜| 久久午夜福利片| 午夜久久久久精精品| 国产精品久久久久久久电影| 亚洲精品国产成人久久av| 国产精品一二三区在线看| 人妻制服诱惑在线中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费在线观看成人毛片| 少妇高潮的动态图| 国产黄a三级三级三级人| 午夜激情福利司机影院| 亚洲欧美日韩东京热| 亚洲精品日韩av片在线观看| 成人无遮挡网站| 国产一区有黄有色的免费视频 | 日本一二三区视频观看| a级一级毛片免费在线观看| 少妇熟女aⅴ在线视频| 黄片wwwwww| 男女边吃奶边做爰视频| 亚洲成人久久爱视频| 国产午夜精品论理片| 日韩中字成人| 亚洲四区av| 亚洲av中文字字幕乱码综合| 日韩不卡一区二区三区视频在线| 在现免费观看毛片| 97超碰精品成人国产| 国产真实伦视频高清在线观看| 成人午夜高清在线视频| 精品一区二区三卡| 精品久久久久久久久亚洲| 性插视频无遮挡在线免费观看| 一级爰片在线观看| 亚洲精品自拍成人| 毛片一级片免费看久久久久| 日韩一区二区视频免费看| 自拍偷自拍亚洲精品老妇| 日日摸夜夜添夜夜爱| 国产综合懂色| 亚洲成人精品中文字幕电影| 毛片一级片免费看久久久久| 国产亚洲av片在线观看秒播厂 | 99热6这里只有精品| 国产高清有码在线观看视频| 午夜老司机福利剧场| 国产精品国产三级国产专区5o| 久久久久性生活片| 最近最新中文字幕大全电影3| 91精品一卡2卡3卡4卡| 成年女人在线观看亚洲视频 | 老司机影院成人| 91精品伊人久久大香线蕉| 国产精品久久久久久av不卡| 欧美zozozo另类| 简卡轻食公司| 婷婷色av中文字幕| 亚洲av成人精品一区久久| 欧美不卡视频在线免费观看| 亚洲精品久久久久久婷婷小说| 最近最新中文字幕免费大全7| 欧美性感艳星| 在现免费观看毛片| av线在线观看网站| 伦精品一区二区三区| 中文天堂在线官网| 日韩中字成人| 久久99热这里只频精品6学生| 国产精品美女特级片免费视频播放器| 国产亚洲av片在线观看秒播厂 | 青青草视频在线视频观看| av线在线观看网站| av播播在线观看一区| 97热精品久久久久久| 成人综合一区亚洲| 国内揄拍国产精品人妻在线| 午夜激情福利司机影院| 亚洲国产高清在线一区二区三| 午夜精品国产一区二区电影 | 午夜福利高清视频| 亚洲国产最新在线播放| 97超碰精品成人国产| 日韩制服骚丝袜av| 久久久久性生活片| 最近中文字幕2019免费版| 一级毛片电影观看| 国产 一区 欧美 日韩| 久久精品人妻少妇| 日日摸夜夜添夜夜爱| 一级爰片在线观看| av播播在线观看一区| 亚洲va在线va天堂va国产| 亚洲欧美成人精品一区二区| 中文精品一卡2卡3卡4更新| 超碰97精品在线观看| 美女脱内裤让男人舔精品视频| 国产视频首页在线观看| 黄色欧美视频在线观看| 成人国产麻豆网| 国产精品一区www在线观看| 九草在线视频观看| 激情 狠狠 欧美| 内地一区二区视频在线| 日日啪夜夜撸| 国模一区二区三区四区视频| 丝袜喷水一区| 你懂的网址亚洲精品在线观看| 欧美三级亚洲精品| h日本视频在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 爱豆传媒免费全集在线观看| 免费观看无遮挡的男女| 高清毛片免费看| 精品国产一区二区三区久久久樱花 | 久久97久久精品| 乱码一卡2卡4卡精品| av网站免费在线观看视频 | 国产一级毛片在线| 99久久九九国产精品国产免费| 丰满人妻一区二区三区视频av| 国产精品久久久久久av不卡| 十八禁国产超污无遮挡网站| 亚洲精品视频女| 观看美女的网站| 熟妇人妻不卡中文字幕| 国产午夜福利久久久久久| 亚洲欧洲日产国产| 国产69精品久久久久777片| 特级一级黄色大片| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品国产av成人精品| 精品国内亚洲2022精品成人| 亚洲真实伦在线观看| 国产精品久久久久久av不卡| 亚洲在线观看片| av网站免费在线观看视频 | a级毛色黄片| 免费看光身美女| 99久久精品热视频| 亚洲av免费高清在线观看| 久久久精品免费免费高清| 国产三级在线视频| 一级毛片久久久久久久久女| 午夜福利成人在线免费观看| 精品国产露脸久久av麻豆 | 最近手机中文字幕大全| 日韩视频在线欧美| 亚洲精品日韩在线中文字幕| 亚洲av一区综合| 少妇熟女欧美另类| 久久亚洲国产成人精品v| 神马国产精品三级电影在线观看| 激情 狠狠 欧美| 日韩一区二区三区影片| 亚洲精品一二三| 久久精品久久久久久噜噜老黄| videossex国产| 免费观看精品视频网站| 亚洲欧美精品专区久久| 亚洲欧美日韩卡通动漫| 男女那种视频在线观看| 国产免费福利视频在线观看| 色综合站精品国产| 国产精品一二三区在线看| 看免费成人av毛片| 国产成人免费观看mmmm| 亚洲最大成人手机在线| 成人综合一区亚洲| 日韩精品青青久久久久久| 校园人妻丝袜中文字幕| 蜜桃久久精品国产亚洲av| 大话2 男鬼变身卡| 国产精品爽爽va在线观看网站| 国产精品国产三级国产专区5o| 高清视频免费观看一区二区 | 国产精品爽爽va在线观看网站| 国产亚洲精品久久久com| 久热久热在线精品观看| 午夜精品一区二区三区免费看| 久久久久国产网址| 午夜福利网站1000一区二区三区| 亚洲怡红院男人天堂| 少妇裸体淫交视频免费看高清| 日韩欧美精品免费久久| 成人二区视频| 男女那种视频在线观看| 久久久久久久久久久免费av| 亚洲精品亚洲一区二区| 久久久久久久久久黄片| 日日摸夜夜添夜夜爱| 丝瓜视频免费看黄片| av天堂中文字幕网| 亚洲av成人精品一二三区| 插阴视频在线观看视频| 中文字幕久久专区| 日韩一本色道免费dvd| 欧美性感艳星| 中文字幕亚洲精品专区| 白带黄色成豆腐渣| 国产成人免费观看mmmm| 亚洲av免费高清在线观看| av专区在线播放| 国产一区二区三区av在线| 大片免费播放器 马上看| 一级爰片在线观看| 国产三级在线视频| 一区二区三区乱码不卡18| 午夜免费观看性视频| 亚洲精品日韩av片在线观看| 最后的刺客免费高清国语| 日韩欧美一区视频在线观看 | 精品国产三级普通话版| 国产亚洲精品久久久com| 特级一级黄色大片| 五月玫瑰六月丁香| 成人鲁丝片一二三区免费| 伦理电影大哥的女人| av免费在线看不卡| 精品久久久久久成人av| 又大又黄又爽视频免费| 亚洲久久久久久中文字幕| av免费观看日本| 天堂俺去俺来也www色官网 | 人妻少妇偷人精品九色| ponron亚洲| 亚洲av免费高清在线观看| 99视频精品全部免费 在线| 国产成人a∨麻豆精品| 一个人看视频在线观看www免费| 91久久精品国产一区二区三区| 成人欧美大片| 午夜激情久久久久久久| 啦啦啦韩国在线观看视频| 国产精品一区二区三区四区久久| 青春草视频在线免费观看| 爱豆传媒免费全集在线观看| 免费少妇av软件| 一级毛片 在线播放| 久久久久久久久久久丰满| 久久精品国产自在天天线| 极品教师在线视频| 欧美bdsm另类| 久久精品国产亚洲av天美| 精品久久久久久久久亚洲| 亚洲精品456在线播放app| av播播在线观看一区| a级毛色黄片| 午夜福利在线在线| 亚洲一区高清亚洲精品| 欧美日韩精品成人综合77777| 免费看日本二区| 亚洲欧美精品专区久久| 99久国产av精品| 日本色播在线视频| 亚洲久久久久久中文字幕| 男女下面进入的视频免费午夜| 18禁动态无遮挡网站| 亚洲成人精品中文字幕电影| 国产精品蜜桃在线观看| 少妇丰满av| 欧美极品一区二区三区四区| 九九久久精品国产亚洲av麻豆| 狠狠精品人妻久久久久久综合| 美女大奶头视频| 蜜桃亚洲精品一区二区三区| 26uuu在线亚洲综合色| 成人亚洲精品av一区二区| 26uuu在线亚洲综合色| 国产91av在线免费观看| 国产伦精品一区二区三区视频9| 亚洲精品乱码久久久v下载方式| 亚洲av成人精品一二三区| 别揉我奶头 嗯啊视频| 99久久精品热视频| 久久久久久久亚洲中文字幕| 2018国产大陆天天弄谢| 18+在线观看网站| 青春草亚洲视频在线观看| 亚洲欧美成人精品一区二区| 日本av手机在线免费观看| 亚洲欧美日韩东京热| 久久97久久精品| 少妇高潮的动态图| 一区二区三区乱码不卡18| 夜夜看夜夜爽夜夜摸| 亚洲最大成人手机在线| 在线观看av片永久免费下载| 国产视频内射| 国产真实伦视频高清在线观看| 99久国产av精品| 久久久久精品久久久久真实原创| 亚洲天堂国产精品一区在线| 国产成人精品婷婷| 精品久久久久久久人妻蜜臀av| 狂野欧美激情性xxxx在线观看| 91精品一卡2卡3卡4卡| 深爱激情五月婷婷| 精品酒店卫生间| 一区二区三区免费毛片| 免费观看av网站的网址| 精品一区二区三区视频在线| 日韩av在线免费看完整版不卡| 亚洲人与动物交配视频| 黄片无遮挡物在线观看| 蜜桃久久精品国产亚洲av| 久久久久久伊人网av| 99热这里只有是精品50| 久久久精品94久久精品| 久久精品国产亚洲av涩爱| 日韩av在线大香蕉| 美女脱内裤让男人舔精品视频| 亚洲最大成人手机在线| 日韩精品青青久久久久久| 久久99热6这里只有精品| 亚洲欧美日韩无卡精品| 直男gayav资源| 国产成人免费观看mmmm| 亚洲丝袜综合中文字幕| 18禁在线播放成人免费| 网址你懂的国产日韩在线| 国语对白做爰xxxⅹ性视频网站| 午夜福利在线观看免费完整高清在| 男人和女人高潮做爰伦理| 边亲边吃奶的免费视频| 国产亚洲一区二区精品| 男人舔女人下体高潮全视频| 国产亚洲av片在线观看秒播厂 | 极品少妇高潮喷水抽搐| 国产av不卡久久| 少妇的逼水好多| 成人漫画全彩无遮挡| 伊人久久精品亚洲午夜| 插阴视频在线观看视频| 日韩人妻高清精品专区| av免费在线看不卡| 久久久久久久亚洲中文字幕| 精品一区在线观看国产| 免费大片黄手机在线观看| 亚洲国产最新在线播放| 天天躁夜夜躁狠狠久久av| 天堂俺去俺来也www色官网 | 99热这里只有是精品在线观看| 五月伊人婷婷丁香| 丝袜喷水一区| or卡值多少钱| 伦理电影大哥的女人| 最近手机中文字幕大全| 日韩精品青青久久久久久| 嘟嘟电影网在线观看| 九九久久精品国产亚洲av麻豆| 精品久久久噜噜| 噜噜噜噜噜久久久久久91| 亚洲最大成人中文| 春色校园在线视频观看| 国产黄色免费在线视频| 精品99又大又爽又粗少妇毛片| av免费在线看不卡| 久久久久久久久中文| 久久99热这里只频精品6学生| 国产男女超爽视频在线观看| 成人av在线播放网站| 蜜臀久久99精品久久宅男| 18禁动态无遮挡网站| 免费电影在线观看免费观看| 国产又色又爽无遮挡免| 在线观看免费高清a一片| 亚洲欧美精品自产自拍| 天堂俺去俺来也www色官网 | 如何舔出高潮| 国产精品国产三级国产av玫瑰| 亚洲av中文字字幕乱码综合| 成人午夜高清在线视频| 身体一侧抽搐| 欧美3d第一页| 菩萨蛮人人尽说江南好唐韦庄| 成人综合一区亚洲| 日本免费a在线| 国产精品1区2区在线观看.| 久久久久网色| 国产精品无大码| 精品久久久久久久久久久久久| 日韩av免费高清视频| av在线亚洲专区| 好男人视频免费观看在线| 午夜福利网站1000一区二区三区| 亚洲精品456在线播放app| av免费在线看不卡| 日韩在线高清观看一区二区三区| 亚洲高清免费不卡视频| 久久精品人妻少妇| 高清毛片免费看| 少妇的逼好多水| 十八禁网站网址无遮挡 | 国产乱人视频| 午夜精品国产一区二区电影 | 偷拍熟女少妇极品色| 一个人看的www免费观看视频| 麻豆成人午夜福利视频| 美女主播在线视频| 国产成人精品一,二区| 噜噜噜噜噜久久久久久91| 亚洲经典国产精华液单| 狂野欧美白嫩少妇大欣赏| 听说在线观看完整版免费高清| 免费无遮挡裸体视频| 久久久色成人| 91aial.com中文字幕在线观看| 国产乱人偷精品视频| 久久久精品94久久精品| 国产黄片美女视频| 亚洲成人久久爱视频| 欧美bdsm另类| av专区在线播放| 日日摸夜夜添夜夜爱| 永久免费av网站大全| 欧美成人精品欧美一级黄| 精品午夜福利在线看| 久久久久久久久久成人| 亚洲四区av| 亚洲熟女精品中文字幕| 国产日韩欧美在线精品| 免费观看性生交大片5| 大话2 男鬼变身卡| 如何舔出高潮| 一边亲一边摸免费视频| 岛国毛片在线播放| 精品人妻熟女av久视频| 久久这里只有精品中国| 免费看不卡的av| 国模一区二区三区四区视频| 汤姆久久久久久久影院中文字幕 | 91av网一区二区| 国产一区有黄有色的免费视频 | 亚洲av中文字字幕乱码综合| 欧美成人a在线观看| 久久精品人妻少妇| 亚洲精品一区蜜桃| 搡老乐熟女国产| 九九在线视频观看精品| 韩国高清视频一区二区三区| 亚州av有码| 久久精品国产亚洲av涩爱| 老司机影院成人| 黄片wwwwww| 秋霞伦理黄片| 18禁在线无遮挡免费观看视频| 亚洲av国产av综合av卡| 国产亚洲av片在线观看秒播厂 | 亚洲av中文字字幕乱码综合| 亚洲高清免费不卡视频| 久久99精品国语久久久| 亚洲美女搞黄在线观看| 高清日韩中文字幕在线| 联通29元200g的流量卡| 成人特级av手机在线观看| 中国美白少妇内射xxxbb| 亚洲av一区综合| 久久久久国产网址| 毛片一级片免费看久久久久| .国产精品久久| 婷婷色综合www| av在线蜜桃| 亚洲av电影不卡..在线观看| 国产亚洲午夜精品一区二区久久 | 精品一区二区免费观看| 国产午夜精品论理片| 久久久色成人| 国内揄拍国产精品人妻在线| 久久久久久久久中文| 亚洲综合精品二区| 亚洲国产精品sss在线观看| 欧美一级a爱片免费观看看| 亚洲精品中文字幕在线视频 | 国产91av在线免费观看| 综合色av麻豆| 波野结衣二区三区在线| 国产片特级美女逼逼视频| 亚洲精品久久久久久婷婷小说| 超碰97精品在线观看| 国产不卡一卡二| 日日啪夜夜爽| 爱豆传媒免费全集在线观看| 日韩电影二区| 在线观看美女被高潮喷水网站| 三级国产精品欧美在线观看| 欧美高清性xxxxhd video| 亚洲av电影在线观看一区二区三区 | 精品酒店卫生间| 国产精品嫩草影院av在线观看| freevideosex欧美| 晚上一个人看的免费电影| 美女国产视频在线观看| 精品一区二区三卡| 国产一区二区亚洲精品在线观看| 午夜视频国产福利| 视频中文字幕在线观看| 久久久久久久久久久免费av| 国产精品一二三区在线看| 一级片'在线观看视频| 亚洲经典国产精华液单| 欧美3d第一页| 夫妻午夜视频| 精品不卡国产一区二区三区| 成人二区视频| 国产精品日韩av在线免费观看| 欧美激情久久久久久爽电影| 日韩一区二区三区影片| 久久人人爽人人爽人人片va| 夫妻性生交免费视频一级片| 免费看日本二区| 淫秽高清视频在线观看|