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

    基于時變AR模型的非平穩(wěn)非高斯隨機過程的數(shù)值模擬

    2016-01-15 03:24:46李錦華,陳水生,吳春鵬
    振動與沖擊 2015年17期
    關(guān)鍵詞:數(shù)值模擬

    第一作者李錦華男,博士, 講師, 碩士生導師,1981年12月生

    通信作者陳水生男,博士, 教授, 碩士生導師,1968年6月生

    基于時變AR模型的非平穩(wěn)非高斯隨機過程的數(shù)值模擬

    李錦華1, 2,陳水生2,吳春鵬2,李建豐2

    (1.華東交通大學鐵路環(huán)境振動與噪聲教育部研究工程中心,南昌330013; 2.華東交通大學土木建筑學院,南昌330013)

    摘要:為了有效地模擬具有目標非平穩(wěn)、非高斯特征的隨機過程,提出了基于時變AR模型的非平穩(wěn)非高斯隨機過程的模擬方法。該方法首先需要建立實現(xiàn)非高斯與高斯隨機過程之間相互轉(zhuǎn)換的非線性平移關(guān)系,然而該非線性平移也會導致平移前后高斯與非高斯隨機過程的功率譜發(fā)生變化。因此該方法還需要進一步建立平移前后高斯與非高斯隨機過程的功率譜或相關(guān)函數(shù)的轉(zhuǎn)換關(guān)系。然后,通過已建立的非線性平移,以及功率譜或相關(guān)函數(shù)的轉(zhuǎn)換關(guān)系,可將非平穩(wěn)非高斯隨機過程的模擬轉(zhuǎn)化成對非平穩(wěn)高斯隨機過程的模擬。而非平穩(wěn)高斯隨機過程可通過建立的時變AR模型進行有效的模擬。最后將具有目標非平穩(wěn)、非高斯特征的脈動風速模擬作為數(shù)值算例,驗證了該方法模擬非平穩(wěn)非高斯隨機過程的有效性。

    關(guān)鍵詞:時變AR模型;數(shù)值模擬;非平穩(wěn)隨機過程;非高斯;脈動風速

    基金項目:國家自然科學基金項目(11162005);江西省自然科學基金項目(20132BAB216003)

    收稿日期:2015-01-07修改稿收到日期:2015-03-11

    中圖分類號:TU311

    文獻標志碼:A

    DOI:10.13465/j.cnki.jvs.2015.17.024

    Abstract:In order to simulate effectively a stochastic process possessing given non-stationary non-Gaussian features, a method to simulate a non-stationary non-Gaussian stochastic process based on a time-varying AR model was proposed here. Firstly, it was necessary to establish a nonlinear translation relationship to realize mutual conversion between a non-Gaussian and a Gaussian one. Meanwhile, the power spectra both the non-Gaussian and the Gaussian random processes were different due to the nonlinear translation. Furthermore, the transformation relationship between their power spectra or correlation functions was established. Then, the simulation of a non-stationary non-Gaussian stochastic process was converted into a simulation of a non-stationary Gaussian random process by utilizing the nonlinear translation relationship and the transformation relationship of their power spectra or correlation functions constructed above. The non-stationary Gaussian random process was effectively simulated with a presented time-varying AR model. Finally, taking the simulation of a fluctuating wind velocity possessing target non-stationary non-Gaussian characteristics as a numerical example, the effectiveness of the method to simulate a non-stationary non-Gaussian random process was verified.

    Simulation of non-stationary non-Gaussian stochastic process based on a time-varying AR model

    LIJin-hua1, 2,CHENShui-sheng2,WUChun-peng2,LIJian-feng2(1. MOE Engineering Research Center of Railway Environmental Vibration and Noise, East China Jiaotong University, Nanchang 330013, China;2. College of Civil Engineering,East China Jiaotong University, Nanchang 330013, China)

    Key words:time-varying AR model; numerical simulation; non-stationary stochastic process; non-Gaussian; fluctuating wind velocity

    計算機技術(shù)的快速發(fā)展,使得運用計算隨機力學解決不確定性系統(tǒng)的分析成為可能[1]。在進行不確定性系統(tǒng)分析過程中,需要首先獲取輸入的不確定性隨機參量?;诿商乜錗onte-Carlo的隨機模擬技術(shù)[2],使得人工模擬系統(tǒng)的隨機輸入得到了廣泛的應用。在模擬產(chǎn)生具有指定目標特征的隨機過程中,不確定性隨機參量常被假定為平穩(wěn)高斯、非平穩(wěn)高斯隨機過程來模擬,這個假定極大地簡化模擬計算過程。在自然界中,不確定性系統(tǒng)的隨機參量如工程結(jié)構(gòu)上的風壓、結(jié)構(gòu)與機械工程中的幾何特性和材料特性、巖土工程中的土壤特性、海洋波浪的隨機激勵等均具有非高斯的特征[3]。為了有效地研究分析隨機激勵下的系統(tǒng)響應,具有目標特征的平穩(wěn)非高斯乃至非平穩(wěn)非高斯隨機過程的有效模擬逐漸受到廣泛關(guān)注。非平穩(wěn)非高斯隨機過程不僅具有相關(guān)函數(shù)或功率譜密度隨時間變化的特征,而且每一時刻的概率密度函數(shù)分布呈現(xiàn)出非高斯分布的特征。高斯隨機過程的概率密度函數(shù)可以完全由前兩階統(tǒng)計參數(shù)( 即均值和方差) 來描述。而對于非高斯隨機過程,特別是非平穩(wěn)非高斯隨機過程,可采用高階統(tǒng)計參數(shù)(特別是三階和四階統(tǒng)計參數(shù),即偏度和峰度)對非高斯概率密度函數(shù)的特征進行描述。模擬非平穩(wěn)非高斯隨機過程,需要模擬具有目標時變特征的相關(guān)函數(shù)或功率譜,以及目標非高斯特征。因此,非平穩(wěn)非高斯隨機過程的模擬,可先通過模擬生成某一非平穩(wěn)高斯隨機過程,然后再將該非平穩(wěn)高斯隨機過程經(jīng)過非線性平移生成具有目標特征的非平穩(wěn)非高斯隨機過程。

    運用Monte-Carlo思想模擬隨機過程主要有兩大類方法,即譜表示法和線性濾波法。譜表示方法的模擬精度較高,但模擬效率較低[4]。而線性濾波法盡管模擬精度較低,但計算量小、速度快,被廣泛用于工程領(lǐng)域隨機過程的模擬[6]。Decolatis等[6]對桅桿結(jié)構(gòu)由不同方法得到的風速時程作用下的結(jié)構(gòu)響應結(jié)果作出比較后,建議采用AR模型模擬互相相關(guān)多重脈動風速時程。Owen等[7]采用AR模型模擬了平穩(wěn)隨機風荷載,并用于斜拉橋的風振響應分析。Saramas等[8]采用ARMA模型模擬了多變量平穩(wěn)隨機過程。以上模擬主要針對的是平穩(wěn)高斯隨機過程的模擬。張文福等[9]通過時變函數(shù)均勻調(diào)制基于AR模型的平穩(wěn)脈動風速來生成了下?lián)舯┝鞣瞧椒€(wěn)脈動風速。Li等[10]通過建立AR和ARMA模型輸入與輸出統(tǒng)計特征之間的關(guān)系,建立了基于AR和ARMA模型的非高斯平穩(wěn)隨機過程模擬,并通過非高斯平穩(wěn)脈動風荷載的模擬驗證了模擬的有效性。本文將進一步探討基于時變AR模型的非平穩(wěn)非高斯隨機過程的有效模擬。

    1非平穩(wěn)高斯隨機過程的時變AR模型

    設(shè)時變AR模型階數(shù)為p,非平穩(wěn)高斯隨機過程模擬樣本點數(shù)為N,采樣時間間隔為Δt,則基于時變AR(p)模型的非平穩(wěn)高斯隨機過程模擬公式可表示為:

    (1)

    1.1時變模型系數(shù)Ai(t)的確定

    在某時刻t=t′時,有

    (2)

    對于Ai(t′)的確定,將式(2)兩邊同時右乘f(t′-jΔt),并取數(shù)學期望,有

    E[f(t′)f(t′-jΔt)]=

    E[L(t′)w(t′)f(t′-jΔt)]

    (3)

    Rff(t′,jΔt)=

    對于右二項,Rwf(t′,jΔt)可理解為w(t′)與f(t′-jΔt)的互相關(guān)函數(shù)。w(t)為該系統(tǒng)的輸入,而f(t)為系統(tǒng)的輸出。當前的輸出只依賴于當前和過去的輸入,而與將來的輸入無關(guān),因此w(t′)與f(t′-jΔt)互相獨立,故

    Rwf(t′,jΔt)=0

    (5)

    將式(5)代入式(4),可確定Ai(t′)

    (6)

    其展開式為

    1.2時變模型系數(shù)L(t)的確定

    在某時刻t=t′時,將式(2)兩邊同時右乘w(t′),并取數(shù)學期望,有

    Rfw(t′,0)=

    (7)

    當前輸入w(t′)與過去輸出f(t′-iΔt)互相獨立,則

    Rfw(t′-iΔt,iΔt)=0

    (8)

    將式(8)代入式(7),可得

    Rfw(t′,0)=L(t′)Rww(t′,0)

    (9)

    根據(jù)白噪聲的特性,有Rww(t′,0)=1,則

    Rfw(t′,0)=L(t′)

    (10)

    再一次將式(2)兩邊同時右乘f(t′),并取數(shù)學期望,有

    L(t′)Rwf(t′,0)

    (11)

    又因為

    Rwf(t′,0)=E[w(t′)f(t′)]=

    E[f(t′)w(t′)]=Rfw(t′,0)

    (12)

    將式(10)和(12)代入式(11),有

    Rff(t′,0)=

    (13)

    令A0(t′)=-1,則

    (14)

    因此,時變模型系數(shù)Ai(t′)和L(t′),t′=1,2,3,…,可分別根據(jù)式(6)和(14)來確定。

    2非線性平移過程

    非高斯隨機過程y(t)可通過潛在高斯隨機過程f(t)進行非線性平移轉(zhuǎn)換。

    h4(f3(t)-3f(t))]

    (15)

    (16)

    (17)

    為了使得式(15)的非線性平移可逆,形狀系數(shù)需滿足

    (18)

    則非高斯隨機過程的偏度和峰度需滿足下列關(guān)系

    (19)

    當具有目標功率譜的標準化高斯隨機過程,經(jīng)過式(15)非線性平移前后,生成的標準化非高斯隨機過程的功率譜必然發(fā)生非線性變化。因此,標準化非高斯隨機過程的目標功率譜不能直接作為高斯隨機過程目標功率譜進行模擬。為此,需要建立標準化非高斯隨機過程的目標功率譜與高斯隨機過程目標功率譜的轉(zhuǎn)化關(guān)系。根據(jù)Pristley的“進化譜”理論[11],一零均值非平穩(wěn)隨機過程可表示為:

    (20)

    式中:A(ω,t)為非均勻調(diào)制函數(shù);Z(ω)為一正交過程,且滿足

    E[dZ(ω)]=0

    (21)

    E[dZ*(ω)dZ(ω′)]=J(ω)δ(ω-ω′)dωdω′

    (22)

    其中:*表示共軛,δ為狄拉克函數(shù)。則非平穩(wěn)隨機過程的均值為

    (23)

    相關(guān)函數(shù)為

    Rff(t,τ)=E[f*(t)f(t+τ)]=

    τ)e-iωteiω′(t+τ)E[dZ*(ω)dZ(ω′)]=

    (24)

    當τ=0時,

    (25)

    因此,進化譜S(ω,t)可通過一時頻函數(shù)A(ω,t)對功率譜的非均勻調(diào)制來表示[12],即

    (26)

    根據(jù)式(24)和式(26),可建立非平穩(wěn)隨機過程的相關(guān)函數(shù)與時變功率譜之間的關(guān)系

    Rff(t,τ)=E[f*(t)f(t+τ)]=

    (27)

    根據(jù)相關(guān)函數(shù)與功率譜之間的相互轉(zhuǎn)換關(guān)系,非高斯隨機過程與高斯隨機過程的目標功率譜之間的非線性關(guān)系,可通過目標相關(guān)函數(shù)之間的非線性關(guān)系進行表達。

    根據(jù)式(15)的非線性平移,標準化的非高斯隨機過程的相關(guān)函數(shù)可表示為

    ρx(t,τ))·dx(t)·dx(t+τ)=

    α[f(t+τ)+h3(f2(t+τ)-1)+

    h4(f3(t+τ)-3f(t+τ))]×

    φ(f(t),f(t+τ),ρf(t,τ))·df(t)·df(t+τ)

    (28)

    其中:ψ為具有相關(guān)系數(shù)ρx的兩非高斯隨機過程的聯(lián)合概率密度函數(shù);φ為具有相關(guān)系數(shù)ρf的兩高斯隨機過程的聯(lián)合概率密度函數(shù)。

    (29)

    而相關(guān)函數(shù)與相關(guān)系數(shù)的轉(zhuǎn)化關(guān)系為

    R(t,τ)=ρ(t,τ)·σ2

    (30)

    根據(jù)式(28)~(30),可建立非高斯隨機過程與高斯隨機過程的相關(guān)系數(shù)或相關(guān)函數(shù)之間的轉(zhuǎn)換關(guān)系。

    3非平穩(wěn)非高斯風荷載的數(shù)值模擬

    為了驗證上述非平穩(wěn)非高斯隨機過程模擬方法的有效性,本節(jié)將以具有非平穩(wěn)非高斯特征脈動風速的模擬作為數(shù)值算例。在數(shù)值模擬過程中,脈動風速的目標非平穩(wěn)特征主要表現(xiàn)在功率譜隨時間變化的時變功率譜S(t,ω),該時變功率譜可采用Davenport非均勻調(diào)制函數(shù)[13]

    (31)

    調(diào)制Davenport譜

    (32)

    目標非高斯特征主要通過偏度γ3、峰度γ4來考慮,分別假定為偏度γ3=-1.2、峰度γ4=8。根據(jù)式(16)、(17)可進一步確定非線性平移變換的形狀系數(shù)h3=-0.140 6、h4=0.075 3,以及參數(shù)α=0.965 1。由該非線性平移變換可知,在生成具有目標非高斯特征的非平穩(wěn)隨機過程之前,首先需要模擬潛在的非平穩(wěn)高斯隨機過程。根據(jù)式(27)~(30),可建立非平穩(wěn)非高斯隨機過程與潛在的非平穩(wěn)高斯隨機過程之間的相關(guān)系數(shù)轉(zhuǎn)換關(guān)系如1圖所示。

    圖1 高斯與非高斯隨機過程相關(guān)系數(shù)之間的轉(zhuǎn)換關(guān)系 Fig.1 The transformation relationship of the correlation coefficients between Gaussian and non-Gaussian random processes

    基于時變AR模型生成潛在的非平穩(wěn)高斯隨機過程見圖2,并將其經(jīng)過非線性平移變換生成的非平穩(wěn)非高斯脈動風速見圖3。對多組脈動風速樣本進行功率譜估計,其估計譜(見圖4)明顯具有如圖5所示的目標時變功率譜的時變特征。任意時刻估計譜與目標譜的對比見圖6,在該圖中四個任意時刻的估計譜均與目標譜基本吻合,且相應的相關(guān)函數(shù)也與目標相吻合見圖7。在時間上統(tǒng)計平均的功率譜、相關(guān)函數(shù)也均與目標基本吻合見圖8、9。因此,模擬的脈動風速具有目標特征的非平穩(wěn)特性。

    圖2 模擬的非平穩(wěn)高斯隨機過程 Fig.2 Thesimulated non-stationary Gaussian stochastic process

    圖3 非線性平移生成的非平穩(wěn)非高斯脈動風速 Fig.3 Thegenerated non-stationary non-Gaussian fluctuating wind velocity through nonlinear translation

    圖4 脈動風速樣本時變譜 Fig.4 The time-varying power spectrum of the fluctuating wind velocity

    圖5 脈動風速的目標時變譜 Fig.5 The target power spectrum of the fluctuating wind velocity

    圖6 脈動風速瞬時估計譜與目標譜的對比 Fig.6 Instantaneous power spectrums of the fluctuating wind velocity with regard to the corresponding targets

    圖7 脈動風速的瞬時相關(guān)函數(shù)與目標相關(guān)函數(shù)的對比 Fig.7 Instantaneous correlation functions of the fluctuating wind velocity with regard to the corresponding targets

    為了進一步說明模擬的有效性,圖10、11展示了脈動風速高階統(tǒng)計特征值(即偏度值、峰度值)的估計。圖中所示的脈動風速的任意時刻偏度和峰度盡管都在其各自均值附近有一定的離散性,但其方差均較小分別為0.2和2.4;而其偏度均值為-1.18接近于目標偏度值-1.2,峰度均值7.7也接近于目標峰度值8,說明了模擬的脈動風速具有目標非高斯特征。因此,模擬的脈動風速具有目標非平穩(wěn)、非高斯特征,說明了非平穩(wěn)非高斯隨機過程模擬方法的有效性。該方法建立在AR模型基礎(chǔ)上,受到AR模型本身模擬精度的影響,盡管在模擬中存在一定誤差,但在工程中是可以接受的。因此,該方法可應用于精度要求較低的工程中。

    圖8 脈動風速統(tǒng)計平均估計譜與目標譜的對比 Fig.8 Temporal average power spectrum of the fluctuating wind velocity with regard to the corresponding target

    圖9 脈動風速統(tǒng)計平均相關(guān)函數(shù)與目標相關(guān)函數(shù)的對比 Fig.9 Temporal average correlation function of the fluctuating wind velocity with regard to the corresponding target

    圖10 脈動風速的偏度 Fig.10 The skewness of the fluctuating wind velocity

    圖11 脈動風速的峰度 Fig.11 The kurtosis of the fluctuating wind velocity

    4結(jié)論

    為了有效地模擬具有目標非平穩(wěn)、非高斯特征的隨機過程,本文提出了基于時變AR模型的非平穩(wěn)非高斯隨機過程的模擬方法。該方法通過建立非高斯與高斯隨機過程之間相互轉(zhuǎn)換的非線性平移關(guān)系,以及非線性平移前后高斯與非高斯隨機過程的功率譜或相關(guān)函數(shù)的轉(zhuǎn)換關(guān)系,可將非平穩(wěn)非高斯隨機過程轉(zhuǎn)化為非平穩(wěn)高斯隨機過程的模擬。而非平穩(wěn)高斯隨機過程可通過建立的時變AR模型進行有效的模擬。為了驗證該方法的有效性,文中進行了具有目標非平穩(wěn)非高斯特征的脈動風速模擬。脈動風速的時變功率譜估計表明:模擬生成的脈動風速樣本的時變功率譜和相關(guān)函數(shù)均與目標基本吻合,模擬的脈動風速具有目標特征的非平穩(wěn)特性。脈動風速高階統(tǒng)計特征值(即偏度值、峰度值)的估計表明:非平穩(wěn)脈動風速的任意時刻偏度和峰度在其各自均值附近有一定的離散性,但其方差均較??;而其偏度均值、峰度均值都非常接近于目標值。因此,模擬的脈動風速具有目標非平穩(wěn)與非高斯特征,說明了非平穩(wěn)非高斯隨機過程模擬方法的有效性。

    參考文獻

    [1]Acharjee S, Zabaras N. Uncertainty propagation in finite deformations-a spectral stochastic Lagrangian approach [J]. Computer Methods in Applied Mechanics and Engineering, 2006, 195: 2289-2312.

    [2]Spanos P D, Zeldin B A. Monte Carlo treatment of random fields: a broad perspective [J]. Applied Mechanics Reviews, 1998, 51:219-237.

    [3]李錦華, 陳水生. 非高斯隨機過程模擬與預測的研究進展[J].華東交通大學學報, 2011, 28(6):1-6.

    LI Jin-hua, CHEN Shui-sheng. Advances in simulation and prediction of non-gaussian stochastic processes[J]. Journal of East China Jiaotong University, 2011, 28(6):1-6.

    [4]李錦華, 李春祥, 申建紅. 非平穩(wěn)脈動風速的數(shù)值模擬[J]. 振動與沖擊, 2009, 28(1):18-23.

    LI Jin-hua, LI Chun-xiang, SHEN Jian-hong. Numerical simulation of non-stationary fluctuating wind velocity[J]. Journal of Vibration and Shock, 2009, 28(1):18-23.

    [5]舒新玲, 周岱. 風速時程AR模型及其快速實現(xiàn)[J].空間結(jié)構(gòu),2003,9(4):27-32.

    SHU Xin-ling, ZHOU Dai. AR model of wind speed time series and its rapid implementation[J]. Spatial Structures, 2003, 9(4):27-32.

    [6]Deodatis G, Micaletti R C. Simulation of highly skewed non-Gaussian stochastic processes [J]. Journal of Engineering Mechanics ASCE, 2001, 127(12): 1284-1295.

    [7]Owen J S, Eccles B J, Choo B S, et al. The application of auto-regressive time series modeling for the time-frequency analysis of civil engineering structures [J]. Engineering Structures, 2001,23:521-536.

    [8]Samaras E, Shinozuka M, Tsurui A. ARMA representation of random process [J].Journal of Engineering Mechanics ASCE, 1985, 111(3): 449-461.

    [9]張文福, 謝丹, 劉迎春, 等. 下?lián)舯┝骺臻g相關(guān)性風場模擬[J]. 振動與沖擊, 2013, 32(10): 12-16.

    ZHANG Wen-fu, XIE Dan, LIU Ying-chun, et al. Simulation of downburst wind field with spatial correlation[J]. Journal of Vibration and Shock, 2013, 32(10): 12-16.

    [10]Li J H, Li C X. Simulation of non-gaussian stochastic process with target power spectral density and lower-order moments [J]. Journal of Engineering Mechanics ASCE, 2012, 138(5):391-404.

    [11]Priestley M B. Power spectral analysis of non-stationary random processes [J]. Journal of Sound and Vibration, 1967, 6:86-97.

    [12]Liang J, Chaudhuri S R, Shinozuka M. Simulation of non-stationary stochastic processes by spectral representation [J]. Journal of Engineering Mechanics ASCE, 2007, 133 (6): 616-627.

    [13]李錦華, 吳春鵬, 陳水生. 下?lián)舯┝鞣瞧椒€(wěn)脈動風速的數(shù)值模擬[J].振動與沖擊, 2014, 33(14):54-60.

    LI Jin-hua, WU Chun-peng, CHEN Shui-sheng. Simulation of non-stationary fluctuating wind velocity in downburst[J]. Journal of Vibration and Shock, 2014, 33(14):54-60.

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學中的應用
    科教導刊(2016年28期)2016-12-12 06:22:00
    基于氣象信息及風場信息的風機輪轂處風速預測
    鉆孔灌注樁樁底沉渣對樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測試、預處理和加工研究報告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    97超视频在线观看视频| 熟妇人妻久久中文字幕3abv| 国内精品久久久久久久电影| 日韩av在线大香蕉| 亚洲无线在线观看| 国产精品三级大全| 国产一区二区激情短视频| 国产一区在线观看成人免费| 三级国产精品欧美在线观看| 久久久久免费精品人妻一区二区| 亚洲在线观看片| 久久久精品大字幕| 免费人成在线观看视频色| 男女那种视频在线观看| 久久久国产成人精品二区| 18禁黄网站禁片午夜丰满| 怎么达到女性高潮| 亚洲av不卡在线观看| 国产成人aa在线观看| 好男人在线观看高清免费视频| 欧美在线一区亚洲| а√天堂www在线а√下载| 天天添夜夜摸| 午夜福利免费观看在线| 亚洲内射少妇av| 99精品欧美一区二区三区四区| 男女做爰动态图高潮gif福利片| 国产成人系列免费观看| 国产成人系列免费观看| 中文字幕人成人乱码亚洲影| 国内揄拍国产精品人妻在线| 国产免费av片在线观看野外av| 婷婷六月久久综合丁香| 亚洲av成人av| 在线免费观看的www视频| 日韩欧美一区二区三区在线观看| 美女免费视频网站| 午夜福利在线在线| 国产一区二区三区视频了| 婷婷精品国产亚洲av| 99国产精品一区二区三区| 看片在线看免费视频| 亚洲在线自拍视频| 国内久久婷婷六月综合欲色啪| av片东京热男人的天堂| 亚洲欧美日韩高清专用| 99久久无色码亚洲精品果冻| 51午夜福利影视在线观看| 免费大片18禁| 精品久久久久久久毛片微露脸| 舔av片在线| 久久久久久久久久黄片| 亚洲激情在线av| 精品国产美女av久久久久小说| 午夜免费观看网址| 日韩 欧美 亚洲 中文字幕| 亚洲精品在线美女| 亚洲av二区三区四区| 国产精品一区二区三区四区免费观看 | 亚洲色图av天堂| 国产欧美日韩精品一区二区| 99国产精品一区二区三区| 激情在线观看视频在线高清| 小说图片视频综合网站| 成人特级黄色片久久久久久久| 亚洲 欧美 日韩 在线 免费| 午夜福利在线观看免费完整高清在 | 一区福利在线观看| 在线视频色国产色| 波多野结衣高清无吗| 精品电影一区二区在线| 日韩欧美国产一区二区入口| 国产午夜精品久久久久久一区二区三区 | 欧美3d第一页| 十八禁人妻一区二区| 精品久久久久久久末码| 最后的刺客免费高清国语| 国产淫片久久久久久久久 | 久久久久久久精品吃奶| 高清毛片免费观看视频网站| 老司机午夜十八禁免费视频| 国产精品永久免费网站| 91在线观看av| 国产亚洲精品久久久久久毛片| 久久久久亚洲av毛片大全| 男女做爰动态图高潮gif福利片| www.999成人在线观看| 91久久精品国产一区二区成人 | 99热6这里只有精品| 内地一区二区视频在线| 99热精品在线国产| 国产亚洲精品av在线| 一个人看视频在线观看www免费 | 噜噜噜噜噜久久久久久91| or卡值多少钱| 丰满乱子伦码专区| 俄罗斯特黄特色一大片| 天堂动漫精品| 欧美bdsm另类| 久久精品国产自在天天线| 大型黄色视频在线免费观看| 三级男女做爰猛烈吃奶摸视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产伦精品一区二区三区四那| 国产在视频线在精品| 麻豆国产97在线/欧美| 日本与韩国留学比较| 男女午夜视频在线观看| 国产精品久久久久久精品电影| 99国产综合亚洲精品| 亚洲国产精品成人综合色| 国产精品 国内视频| 嫩草影院入口| 欧美成人a在线观看| 久久久国产精品麻豆| 久久中文看片网| 午夜福利欧美成人| 精品人妻1区二区| 悠悠久久av| 在线a可以看的网站| 国产亚洲精品综合一区在线观看| 国产激情偷乱视频一区二区| 熟妇人妻久久中文字幕3abv| 老司机在亚洲福利影院| 女同久久另类99精品国产91| 午夜福利在线观看免费完整高清在 | 国产主播在线观看一区二区| 色综合亚洲欧美另类图片| 久久精品国产清高在天天线| 美女高潮喷水抽搐中文字幕| 97超视频在线观看视频| 十八禁网站免费在线| 老汉色av国产亚洲站长工具| 精品乱码久久久久久99久播| 97超视频在线观看视频| 久久香蕉国产精品| 亚洲狠狠婷婷综合久久图片| 90打野战视频偷拍视频| 亚洲熟妇熟女久久| 国产精品日韩av在线免费观看| 18禁美女被吸乳视频| www.www免费av| 午夜精品久久久久久毛片777| 无人区码免费观看不卡| 欧美黄色淫秽网站| 国产精品永久免费网站| 观看免费一级毛片| 国产在线精品亚洲第一网站| 午夜视频国产福利| 欧美日韩国产亚洲二区| 亚洲国产精品合色在线| 动漫黄色视频在线观看| 日韩欧美精品v在线| 久久这里只有精品中国| 日本 av在线| 深夜精品福利| 超碰av人人做人人爽久久 | 精品熟女少妇八av免费久了| 欧美色欧美亚洲另类二区| 国产成人aa在线观看| 无人区码免费观看不卡| 少妇人妻精品综合一区二区 | 天堂av国产一区二区熟女人妻| 丁香六月欧美| 真人做人爱边吃奶动态| 91在线观看av| 波野结衣二区三区在线 | 国产精品一及| 欧美性猛交╳xxx乱大交人| 99riav亚洲国产免费| 国产精品美女特级片免费视频播放器| 亚洲精品456在线播放app | 国产免费男女视频| 久久久久免费精品人妻一区二区| 免费av观看视频| 免费av观看视频| 国产激情偷乱视频一区二区| 国产精品野战在线观看| 中文字幕高清在线视频| 精品福利观看| 无遮挡黄片免费观看| 国产91精品成人一区二区三区| 亚洲内射少妇av| 我要搜黄色片| 一级作爱视频免费观看| 黄色片一级片一级黄色片| 亚洲男人的天堂狠狠| 国产欧美日韩精品一区二区| 亚洲国产精品合色在线| 成人av一区二区三区在线看| 国产一区在线观看成人免费| 亚洲国产欧洲综合997久久,| 免费人成视频x8x8入口观看| 精品欧美国产一区二区三| 在线看三级毛片| 色av中文字幕| 国产精品电影一区二区三区| 亚洲无线观看免费| 国产黄片美女视频| 国产高潮美女av| 久久久久久大精品| 一边摸一边抽搐一进一小说| 性欧美人与动物交配| 看片在线看免费视频| 久久精品91蜜桃| 国产精品香港三级国产av潘金莲| 国产三级黄色录像| 最近最新免费中文字幕在线| 熟女人妻精品中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 黄色丝袜av网址大全| 好看av亚洲va欧美ⅴa在| av在线蜜桃| www日本黄色视频网| 一二三四社区在线视频社区8| 制服丝袜大香蕉在线| 亚洲人与动物交配视频| 看片在线看免费视频| 国产精品1区2区在线观看.| 日韩免费av在线播放| 深夜精品福利| 香蕉丝袜av| 欧美成人a在线观看| 午夜激情福利司机影院| 午夜福利视频1000在线观看| 国产亚洲精品久久久com| 日韩av在线大香蕉| 久久精品91无色码中文字幕| 午夜精品久久久久久毛片777| 99久久综合精品五月天人人| 精品国内亚洲2022精品成人| 亚洲最大成人中文| 午夜福利高清视频| 免费av毛片视频| 舔av片在线| 在线观看免费午夜福利视频| 国产久久久一区二区三区| 小蜜桃在线观看免费完整版高清| 国语自产精品视频在线第100页| 日本免费一区二区三区高清不卡| 草草在线视频免费看| 听说在线观看完整版免费高清| 嫁个100分男人电影在线观看| 99久国产av精品| 天堂av国产一区二区熟女人妻| h日本视频在线播放| 波多野结衣高清无吗| 亚洲av免费在线观看| 欧美在线一区亚洲| 在线看三级毛片| 国产色爽女视频免费观看| 内射极品少妇av片p| 人妻夜夜爽99麻豆av| 国产免费av片在线观看野外av| 欧美一区二区亚洲| 禁无遮挡网站| 亚洲男人的天堂狠狠| 欧美+日韩+精品| 国产精品综合久久久久久久免费| www.色视频.com| 国产伦人伦偷精品视频| 18禁黄网站禁片午夜丰满| 亚洲天堂国产精品一区在线| 又黄又爽又免费观看的视频| 国产精品久久久久久亚洲av鲁大| 五月玫瑰六月丁香| 午夜免费激情av| 国产精品久久电影中文字幕| 好看av亚洲va欧美ⅴa在| 国产野战对白在线观看| 最近最新中文字幕大全免费视频| 亚洲不卡免费看| 日本成人三级电影网站| 久久久成人免费电影| 超碰av人人做人人爽久久 | 日本撒尿小便嘘嘘汇集6| 99热精品在线国产| 不卡一级毛片| 欧美黑人欧美精品刺激| 欧美成人a在线观看| 内地一区二区视频在线| 免费看光身美女| www日本黄色视频网| 国产精品女同一区二区软件 | 精品久久久久久,| 超碰av人人做人人爽久久 | 小蜜桃在线观看免费完整版高清| 国产极品精品免费视频能看的| 少妇丰满av| 久久久久久久精品吃奶| 首页视频小说图片口味搜索| 搞女人的毛片| 亚洲自拍偷在线| 亚洲av成人av| 亚洲欧美日韩高清专用| 日韩大尺度精品在线看网址| 三级国产精品欧美在线观看| 国产精品免费一区二区三区在线| 一个人观看的视频www高清免费观看| 最新在线观看一区二区三区| 国产激情偷乱视频一区二区| 每晚都被弄得嗷嗷叫到高潮| 五月玫瑰六月丁香| 黄片大片在线免费观看| 国产色爽女视频免费观看| 欧美最新免费一区二区三区 | 狂野欧美白嫩少妇大欣赏| 国产成人a区在线观看| 悠悠久久av| 90打野战视频偷拍视频| 中文字幕人妻熟人妻熟丝袜美 | 国产日本99.免费观看| a级一级毛片免费在线观看| 精品人妻偷拍中文字幕| 国产黄色小视频在线观看| 午夜精品久久久久久毛片777| 9191精品国产免费久久| 久久精品国产亚洲av香蕉五月| 国产一区二区三区在线臀色熟女| 亚洲成a人片在线一区二区| 老司机深夜福利视频在线观看| 精品国产美女av久久久久小说| 精品熟女少妇八av免费久了| 91九色精品人成在线观看| 国产伦精品一区二区三区视频9 | 国产亚洲精品综合一区在线观看| 色噜噜av男人的天堂激情| 最近在线观看免费完整版| 国产精品国产高清国产av| 岛国视频午夜一区免费看| x7x7x7水蜜桃| 伊人久久精品亚洲午夜| 搡女人真爽免费视频火全软件 | 十八禁人妻一区二区| 两个人看的免费小视频| 此物有八面人人有两片| 亚洲七黄色美女视频| 精品久久久久久久久久久久久| 在线观看免费视频日本深夜| 啦啦啦免费观看视频1| 亚洲av日韩精品久久久久久密| 一进一出抽搐动态| 亚洲最大成人中文| 午夜免费观看网址| 给我免费播放毛片高清在线观看| 亚洲成a人片在线一区二区| 黄色丝袜av网址大全| 亚洲无线观看免费| 亚洲精品一卡2卡三卡4卡5卡| 国产爱豆传媒在线观看| 国产私拍福利视频在线观看| 深爱激情五月婷婷| 最近最新免费中文字幕在线| 亚洲国产色片| 精品日产1卡2卡| 欧美日韩国产亚洲二区| 久久久久久人人人人人| 亚洲精品美女久久久久99蜜臀| 在线观看日韩欧美| netflix在线观看网站| 97碰自拍视频| a级毛片a级免费在线| 午夜激情福利司机影院| 99国产精品一区二区三区| 久久婷婷人人爽人人干人人爱| 成年人黄色毛片网站| 国内精品一区二区在线观看| 香蕉av资源在线| 中亚洲国语对白在线视频| 国产成人av教育| 成人精品一区二区免费| 日韩精品青青久久久久久| 久久精品影院6| 日本 欧美在线| 99久国产av精品| 午夜日韩欧美国产| 一进一出抽搐gif免费好疼| 精品乱码久久久久久99久播| 欧美黑人巨大hd| 97超级碰碰碰精品色视频在线观看| 黑人欧美特级aaaaaa片| 国产熟女xx| av在线蜜桃| 国产伦精品一区二区三区视频9 | 综合色av麻豆| 亚洲中文字幕日韩| 亚洲成人精品中文字幕电影| 天美传媒精品一区二区| 欧美极品一区二区三区四区| 亚洲最大成人中文| 久久这里只有精品中国| 18美女黄网站色大片免费观看| 老司机午夜福利在线观看视频| 国产成+人综合+亚洲专区| 亚洲自拍偷在线| 亚洲性夜色夜夜综合| 国产成人aa在线观看| 国产精品一区二区免费欧美| 岛国在线观看网站| 日本五十路高清| 国产精品爽爽va在线观看网站| 给我免费播放毛片高清在线观看| 99久久综合精品五月天人人| 国产亚洲精品一区二区www| 不卡一级毛片| 色吧在线观看| 久久午夜亚洲精品久久| 欧美日韩亚洲国产一区二区在线观看| 在线十欧美十亚洲十日本专区| 国产精品影院久久| 一个人免费在线观看的高清视频| 1024手机看黄色片| 男女那种视频在线观看| 亚洲电影在线观看av| 男女做爰动态图高潮gif福利片| 日日干狠狠操夜夜爽| 精品久久久久久成人av| 99国产精品一区二区三区| 91av网一区二区| 国产黄色小视频在线观看| 国产成人av教育| 亚洲不卡免费看| 在线观看免费视频日本深夜| 一a级毛片在线观看| 免费大片18禁| 欧美日韩福利视频一区二区| 成人18禁在线播放| 国产免费男女视频| 日韩欧美一区二区三区在线观看| 亚洲五月天丁香| 亚洲天堂国产精品一区在线| 亚洲国产高清在线一区二区三| 国产精品久久久人人做人人爽| 18美女黄网站色大片免费观看| 国产久久久一区二区三区| 全区人妻精品视频| 色综合亚洲欧美另类图片| 噜噜噜噜噜久久久久久91| 日韩亚洲欧美综合| 国产精品久久电影中文字幕| 中文字幕久久专区| 免费在线观看影片大全网站| 热99在线观看视频| 亚洲av熟女| 欧美极品一区二区三区四区| 亚洲av成人精品一区久久| 国产毛片a区久久久久| 色吧在线观看| 麻豆国产av国片精品| 啪啪无遮挡十八禁网站| 国产美女午夜福利| 久久久久久久久中文| 精品久久久久久,| 亚洲精华国产精华精| 欧美最黄视频在线播放免费| 亚洲精品一区av在线观看| 99国产精品一区二区蜜桃av| 日本一本二区三区精品| av天堂在线播放| 亚洲精品影视一区二区三区av| 欧美大码av| 一进一出抽搐gif免费好疼| 深爱激情五月婷婷| 国产成人啪精品午夜网站| 人妻夜夜爽99麻豆av| 18禁美女被吸乳视频| 久久国产精品人妻蜜桃| 我的老师免费观看完整版| 亚洲性夜色夜夜综合| e午夜精品久久久久久久| 少妇熟女aⅴ在线视频| 18禁国产床啪视频网站| 热99re8久久精品国产| 成年版毛片免费区| 天堂√8在线中文| 国产精品亚洲av一区麻豆| av福利片在线观看| 国产一区二区激情短视频| 校园春色视频在线观看| 精品人妻偷拍中文字幕| 香蕉丝袜av| a级一级毛片免费在线观看| 美女高潮喷水抽搐中文字幕| 高清在线国产一区| 一级a爱片免费观看的视频| 老司机午夜十八禁免费视频| 全区人妻精品视频| 乱人视频在线观看| 日本一本二区三区精品| 色在线成人网| 国产亚洲精品久久久com| 免费在线观看日本一区| 91av网一区二区| 男人和女人高潮做爰伦理| 内射极品少妇av片p| 国产精品一及| av女优亚洲男人天堂| 淫妇啪啪啪对白视频| 美女大奶头视频| 免费观看的影片在线观看| 欧美乱色亚洲激情| 国产精品av视频在线免费观看| 亚洲欧美日韩高清专用| 亚洲无线观看免费| 国产精品综合久久久久久久免费| 日本a在线网址| 国产欧美日韩精品一区二区| 国产主播在线观看一区二区| 国产精品美女特级片免费视频播放器| h日本视频在线播放| 免费高清视频大片| 亚洲人与动物交配视频| 99久久九九国产精品国产免费| 天堂√8在线中文| 亚洲精品乱码久久久v下载方式 | 亚洲狠狠婷婷综合久久图片| 99久久综合精品五月天人人| 日韩免费av在线播放| 欧美成人a在线观看| 国产色婷婷99| 国产乱人视频| 亚洲精品久久国产高清桃花| 在线国产一区二区在线| 亚洲片人在线观看| 免费观看的影片在线观看| 欧美高清成人免费视频www| 国产私拍福利视频在线观看| 亚洲成a人片在线一区二区| 日韩中文字幕欧美一区二区| 免费大片18禁| 久久久精品大字幕| 国产主播在线观看一区二区| 无人区码免费观看不卡| 亚洲欧美日韩东京热| 五月玫瑰六月丁香| 成熟少妇高潮喷水视频| av天堂中文字幕网| 国产黄色小视频在线观看| a级一级毛片免费在线观看| av中文乱码字幕在线| 一本久久中文字幕| 日韩国内少妇激情av| 亚洲熟妇熟女久久| 老司机福利观看| 亚洲真实伦在线观看| 给我免费播放毛片高清在线观看| 日本熟妇午夜| 国产精品久久久久久精品电影| 欧美色视频一区免费| 最近最新中文字幕大全电影3| 亚洲真实伦在线观看| 免费看十八禁软件| 国产欧美日韩精品一区二区| 色综合亚洲欧美另类图片| 女人被狂操c到高潮| 中文字幕av成人在线电影| 亚洲五月天丁香| 国产精品国产高清国产av| 香蕉丝袜av| 久久精品国产综合久久久| 国产国拍精品亚洲av在线观看 | 亚洲成人久久爱视频| 变态另类成人亚洲欧美熟女| 18禁在线播放成人免费| 日日夜夜操网爽| 91麻豆精品激情在线观看国产| 久久久色成人| 中文亚洲av片在线观看爽| 亚洲av电影在线进入| 国产伦精品一区二区三区四那| 天堂动漫精品| 可以在线观看的亚洲视频| 国产99白浆流出| 久久久久国内视频| 黄色日韩在线| 成人午夜高清在线视频| 国产精品99久久久久久久久| 欧美日韩中文字幕国产精品一区二区三区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲av电影不卡..在线观看| 亚洲国产精品sss在线观看| 色视频www国产| 亚洲精品在线观看二区| 99热6这里只有精品| 999久久久精品免费观看国产| 国产成人影院久久av| 欧美区成人在线视频| 免费人成视频x8x8入口观看| 国内精品一区二区在线观看| 午夜亚洲福利在线播放| 搡老熟女国产l中国老女人| 欧美午夜高清在线| 亚洲欧美激情综合另类| 精品人妻偷拍中文字幕| 成人性生交大片免费视频hd| 99久久久亚洲精品蜜臀av| 久久久国产成人免费| 亚洲一区二区三区色噜噜| 久久性视频一级片| 精品国产超薄肉色丝袜足j| 给我免费播放毛片高清在线观看| 夜夜看夜夜爽夜夜摸| 俺也久久电影网| 久久香蕉精品热| 十八禁人妻一区二区| 国产老妇女一区| 性色avwww在线观看| 国产视频一区二区在线看| 成人无遮挡网站| 久久国产精品影院| 免费看十八禁软件| 国产精品久久久久久久电影 | 久久久久久久久中文| 久久久精品欧美日韩精品| 亚洲国产日韩欧美精品在线观看 |