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

    應(yīng)用實(shí)測資料的橋址區(qū)設(shè)計風(fēng)速推算方法改進(jìn)

    2022-03-28 07:47:54孫銘悅
    關(guān)鍵詞:模擬退火準(zhǔn)則風(fēng)速

    張 田,孫銘悅,何 理

    (大連海事大學(xué) 交通運(yùn)輸工程學(xué)院,遼寧 大連 116026)

    0 引言

    抗風(fēng)設(shè)計是橋梁設(shè)計的重要組成部分,橋址區(qū)設(shè)計風(fēng)速的推算是橋梁抗風(fēng)設(shè)計的基礎(chǔ),是保障橋梁抗風(fēng)安全首要解決的問題.目前,推算橋址區(qū)風(fēng)速的一般方法是按實(shí)測資料序列計算統(tǒng)計參數(shù),包括均值、標(biāo)準(zhǔn)差σ、離差系數(shù)Cv、偏差系數(shù)Cs,由以上參數(shù)確定一條理論頻率分布曲線(極值-I型曲線或Pearson-III型曲線)[1-4],再計算不同重現(xiàn)期下的設(shè)計風(fēng)速值.然而,由于實(shí)測資料觀測年限與總體相比相差很遠(yuǎn),而且經(jīng)驗(yàn)頻率點(diǎn)分布不是一條光滑連續(xù)曲線,直接由統(tǒng)計參數(shù)得到的理論頻率分布曲線與經(jīng)驗(yàn)頻率點(diǎn)偏離較大.因此,借鑒水文學(xué)中推算設(shè)計流量的方法[5],即適線法,通過調(diào)整統(tǒng)計參數(shù)值,以一定的準(zhǔn)則選配出一條與經(jīng)驗(yàn)頻率點(diǎn)擬合最好的理論頻率分布曲線,進(jìn)而推算設(shè)計風(fēng)速.

    1 橋址區(qū)風(fēng)速理論頻率分布曲線

    由于氣象記錄中的極大風(fēng)速值是觀測期內(nèi)的極大值,工程中可能遇到的極端風(fēng)速情況不能簡單地取氣象資料中觀測到的極大值,而需要利用最大風(fēng)速的頻率分布來推算橋位處的風(fēng)速極值.常用的風(fēng)速頻率分布曲線形式有極值-I型和Pearson-III型[6-9].

    1.1 極值-I型分布曲線

    極值-I型分布函數(shù)為

    其超過保證率的函數(shù)為

    式中,x為隨機(jī)變量;xp為頻率p時的隨機(jī)變量取值;α和μ分別為分布的尺度參數(shù)和位置參數(shù),只要利用已知的樣本序列x1,x2,…,xn-1,xn,合理估計出參數(shù)α和μ的值,則F(x)就可確定.

    重現(xiàn)期為T年(累積頻率為p=1/T)的風(fēng)速值可表示為

    對于參數(shù)α和μ的估計,可以采用矩法.由數(shù)學(xué)期望和方差的計算公式,進(jìn)一步推導(dǎo)得

    一般采用有限樣本容量的均值和標(biāo)準(zhǔn)差作為理論值Ex和σ的近似估計.

    1.2 Pearson-III型分布曲線

    設(shè)風(fēng)速值為隨機(jī)變量x,則Pearson-III型曲線的密度函數(shù)為

    式中,Γ(α)為α的伽馬函數(shù),α、β和x0為曲線的參數(shù),可以用樣本序列的3個統(tǒng)計參數(shù)即均值、離差系數(shù)Cv和偏差系數(shù)Cs來表示

    式中,

    式中,n為樣本數(shù),σ為樣本推算總體的標(biāo)準(zhǔn)差.

    計算時不計算α、β和x0,只計算、Cv和Cs,以及給定工程設(shè)計的標(biāo)準(zhǔn),即重現(xiàn)期T的設(shè)計要求,則累積頻率p=1/T,于是有

    對式(8)進(jìn)行積分計算,就可以得到頻率分布曲線縱坐標(biāo)值,即累積頻率p對應(yīng)的隨機(jī)變量值為

    式中,Φ為離均系數(shù),是頻率和偏差系數(shù)的函數(shù),

    2 風(fēng)速頻率分布曲線的適線調(diào)整

    2.1 適線法及其準(zhǔn)則

    適線法為選定最優(yōu)理論頻率曲線的方法,即通過不斷優(yōu)化調(diào)整統(tǒng)計參數(shù),繪制多條理論頻率曲線,選擇其中與實(shí)測資料吻合最好的曲線作為最終的理論頻率曲線.適線之前需要繪制經(jīng)驗(yàn)頻率點(diǎn)圖,其步驟為:將已知的風(fēng)速序列不論年序按數(shù)值大小遞減次序排列x1≥x2≥…≥xi≥…≥xn-1≥xn;按Weibull公式,即式(10),計算各風(fēng)速值的經(jīng)驗(yàn)頻率,并在海森機(jī)率格紙上繪出經(jīng)驗(yàn)頻率點(diǎn)群.

    式中,P為風(fēng)速系列中第i項(xiàng)隨機(jī)變量的經(jīng)驗(yàn)頻率;n為風(fēng)速系列的容量;i為計算隨機(jī)變量的序號.

    根據(jù)概率論與數(shù)理統(tǒng)計相關(guān)計算方法,適線準(zhǔn)則可以采用離(殘)差平方和最小準(zhǔn)則(OLS準(zhǔn)則)和離(殘)差絕對值和最小準(zhǔn)則(ABS準(zhǔn)則).

    ①由OLS準(zhǔn)則得目標(biāo)函數(shù)為

    ②由ABS準(zhǔn)則得目標(biāo)函數(shù)為

    2.2 對極值-I型曲線的適線調(diào)整

    對極值-I型分布曲線,由上述兩種不同的適線準(zhǔn)則給出的目標(biāo)函數(shù),將式(3)代入式(11)、式(12)得

    以O(shè)LS準(zhǔn)則為例給出計算過程,若Q1取最小值,則令和?Q1?α=0,則有

    式(15)是關(guān)于α和μ的非線性方程組,可采用迭代法求解.令

    設(shè)迭代的初始值為 1α和1μ,將式(4)計算結(jié)果作為初始值.將 1φ和 2φ在(1α,1μ)的鄰域作Taylor展開,略去高次項(xiàng)后的迭代方程組為

    直到滿足下式時迭代終止,即

    式中,1ε和 2ε為預(yù)先給定的允許誤差.

    2.3 對Pearson-III型曲線的適線調(diào)整

    對Pearson-III型分布曲線,將式(9)代入式(11)、式(12)得

    以O(shè)LS準(zhǔn)則為例,若Q1取最小值,則令和,有

    由于經(jīng)驗(yàn)頻率pi為已知,所以給定一個Cs,可以查表求出與pi相對應(yīng)的Φi,由式(22)按上述迭代求解方法可解出此時使Q1最小的和Cv值.為此,將Cs從0開始,按步長Δ增長(取Δ為0.001),從而得到不同組值,每一組值都對應(yīng)一條Pearson-III型理論頻率曲線,也對應(yīng)一個按式(20)計算的Q1值,最后選用使Q1最小的一組作為計算結(jié)果.

    2.4 模擬退火算法計算統(tǒng)計參數(shù)

    對極值-I型或Pearson-III型分布曲線的參數(shù),雖然可以由上述的迭代方法求解,但從推導(dǎo)過程可以看出,該方法需要推求偏導(dǎo)數(shù)矩陣,且涉及到矩陣的求逆過程,求解復(fù)雜,收斂速度較慢.可以采用模擬退火算法[10-13],該方法對優(yōu)化目標(biāo)要求較低,一般只需要知道其數(shù)值關(guān)系,不受函數(shù)連續(xù)性、光滑性的限制,避免了大量求導(dǎo)數(shù)計算.

    采用模擬退火算法對于一般無約束優(yōu)化求最小值問題的求解過程[14-17]如下.

    (1)給定起始溫度T0、終止溫度Tfinal和初始解x0.

    (2)在x0的鄰域內(nèi)模擬產(chǎn)生隨機(jī)擾動Δx.

    (3)計算擾動引起的目標(biāo)函數(shù)值的變化ΔE.

    (4)若ΔE≤0,接受新值x→x0+Δx,轉(zhuǎn)向步驟(5);否則,若exp(ΔE/T)>rand(0,1),rand(0,1)表示0~1內(nèi)的隨機(jī)數(shù),也接受新值x→x0+Δx,轉(zhuǎn)向步驟(5);否則就拒絕x→x0+Δx,轉(zhuǎn)向步驟(6).

    (5)按退火策略進(jìn)行降溫T→update(T).

    采用指數(shù)退溫方式,設(shè)tk、tk-1分別為第k和k-1次迭代時的溫度,t0為初始溫度,η為參數(shù),tk=η·tk-1.

    (6)是否滿足T>Tfinal?若滿足,則輸出求解結(jié)果,結(jié)束迭代計算;否則,轉(zhuǎn)向步驟(2).

    3 計算實(shí)例

    作為實(shí)例,以貴州省興仁縣氣象站1966-1998年的實(shí)測風(fēng)速資料按不同的推斷方法推算氣象站10 m高度處的設(shè)計風(fēng)速值,經(jīng)高度換算后的實(shí)測資料[7]列于表1,其中1975年、1992年的資料缺失.

    表1 興仁縣10 m高度,年最大10 min平均風(fēng)速 Tab.1 annual maximum 10 min average wind speed data at 10 m altitude in Xingren County

    根據(jù)表1的觀測資料可以計算經(jīng)驗(yàn)頻率點(diǎn)群圖,同時計算各統(tǒng)計參數(shù)的初始值,即α0、μ0,、Cv0、Cs0.為了使理論頻率曲線與經(jīng)驗(yàn)頻率點(diǎn)群吻合得更好,采用文中提出的模擬退火算法,按2種適線準(zhǔn)則逐步調(diào)整統(tǒng)計參數(shù)以獲得最優(yōu)解,計算的統(tǒng)計參數(shù)和目標(biāo)函數(shù)結(jié)果列于表2.

    表2 不同方法獲得的參數(shù)及目標(biāo)函數(shù)值 Tab.2 parameters and objective function values obtained by different methods

    從表2可以發(fā)現(xiàn),無論采用那種準(zhǔn)則,應(yīng)用模擬退火算法計算出的統(tǒng)計參數(shù)均優(yōu)于直接計算方法.利用上述方法求出的統(tǒng)計參數(shù),來擬合理論頻率曲線,見圖1,進(jìn)而求出不同重現(xiàn)期的設(shè)計風(fēng)速列于表3.分布曲線的吻合程度可以通過直接觀察分布曲線與經(jīng)驗(yàn)頻率點(diǎn)群圖像或判定系數(shù)檢驗(yàn)、誤差平方和檢驗(yàn)等方法確定.在圖1中將理論頻率分布曲線與經(jīng)驗(yàn)頻率點(diǎn)群分布比較,可以直觀地看出Pearson-III型曲線作為理論頻率分布曲線與經(jīng)驗(yàn)頻率點(diǎn)群吻合得更好.

    圖1 經(jīng)驗(yàn)頻率點(diǎn)及擬合頻率 Fig.1 empirical frequency point and fitting frequency

    表3 不同重現(xiàn)期時采用模擬退火算法推算的設(shè)計風(fēng)速 Tab.3 design wind speed calculated by simulated annealing algorithm in different return periods

    由表3和圖1可以發(fā)現(xiàn),若采用Pearson-III型曲線作為理論頻率分布曲線,運(yùn)用模擬退火算法計算出的統(tǒng)計參數(shù),所得的不同重現(xiàn)期時的設(shè)計風(fēng)速差別不大,可見模擬退火算法計算統(tǒng)計參數(shù)是有效可行的.

    4 結(jié)論

    (1)提出適線法,通過不斷調(diào)整統(tǒng)計參數(shù),使理論頻率曲線與經(jīng)驗(yàn)頻率點(diǎn)群盡可能吻合良好,同時給出相應(yīng)的適線準(zhǔn)則及采用模擬退火算法優(yōu)化統(tǒng)計參數(shù)的步驟.

    (2)通過實(shí)例計算,得到:①基于經(jīng)驗(yàn)頻率點(diǎn)群,采用適線法優(yōu)化統(tǒng)計參數(shù)比由觀測數(shù)據(jù)直接計算統(tǒng)計參數(shù)能使理論頻率曲線與經(jīng)驗(yàn)頻率點(diǎn)群擬合得更好;②給出適線準(zhǔn)則:離(殘)差平方和最小準(zhǔn)則(OLS準(zhǔn)則)和離(殘)差絕對值和最小準(zhǔn)則(ABS準(zhǔn)則),兩種目標(biāo)函數(shù)都能很好用于優(yōu)化統(tǒng)計參數(shù);③采用模擬退火算法,可避免求解統(tǒng)計參數(shù)時大量的求導(dǎo)運(yùn)算,迭代過程簡單,對兩種適線準(zhǔn)則,都能使目標(biāo)函數(shù)快速收斂.

    猜你喜歡
    模擬退火準(zhǔn)則風(fēng)速
    基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測
    基于最優(yōu)TS評分和頻率匹配的江蘇近海風(fēng)速訂正
    海洋通報(2020年5期)2021-01-14 09:26:54
    具非線性中立項(xiàng)的二階延遲微分方程的Philos型準(zhǔn)則
    模擬退火遺傳算法在機(jī)械臂路徑規(guī)劃中的應(yīng)用
    基于Canny振蕩抑制準(zhǔn)則的改進(jìn)匹配濾波器
    基于GARCH的短時風(fēng)速預(yù)測方法
    基于模糊自適應(yīng)模擬退火遺傳算法的配電網(wǎng)故障定位
    SOA結(jié)合模擬退火算法優(yōu)化電容器配置研究
    一圖讀懂《中國共產(chǎn)黨廉潔自律準(zhǔn)則》
    考慮風(fēng)速分布與日非平穩(wěn)性的風(fēng)速數(shù)據(jù)預(yù)處理方法研究
    亚洲成人国产一区在线观看| 欧美日韩视频精品一区| 国产91精品成人一区二区三区| а√天堂www在线а√下载| 成人手机av| 欧美成人免费av一区二区三区| 国产成人精品久久二区二区免费| 亚洲精品国产区一区二| 一二三四社区在线视频社区8| 最近最新免费中文字幕在线| 亚洲成国产人片在线观看| 欧美大码av| 免费久久久久久久精品成人欧美视频| 中文字幕人妻熟女乱码| 亚洲精华国产精华精| 69av精品久久久久久| 日本黄色日本黄色录像| 亚洲,欧美精品.| 精品乱码久久久久久99久播| 免费在线观看完整版高清| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看免费午夜福利视频| 国产成人一区二区三区免费视频网站| 免费在线观看亚洲国产| 日本撒尿小便嘘嘘汇集6| www国产在线视频色| 精品欧美一区二区三区在线| 国产高清videossex| 我的亚洲天堂| 一级毛片高清免费大全| 免费在线观看影片大全网站| 十八禁网站免费在线| 青草久久国产| 亚洲 欧美一区二区三区| 久久精品91蜜桃| 国产精品久久久人人做人人爽| 69av精品久久久久久| 99精品久久久久人妻精品| 日韩高清综合在线| 国产精品av久久久久免费| 久久中文字幕一级| 日本精品一区二区三区蜜桃| 日本 av在线| 美女高潮喷水抽搐中文字幕| 女生性感内裤真人,穿戴方法视频| 91成年电影在线观看| 久久中文字幕人妻熟女| 一夜夜www| 国产成人影院久久av| 久久香蕉激情| 日本撒尿小便嘘嘘汇集6| 老司机福利观看| 制服人妻中文乱码| 看片在线看免费视频| 精品国产国语对白av| 亚洲欧美激情在线| 国产精品免费一区二区三区在线| 制服诱惑二区| 久久久久久久久中文| 日韩免费av在线播放| 国产一区二区三区综合在线观看| 国产熟女午夜一区二区三区| 午夜日韩欧美国产| 亚洲国产毛片av蜜桃av| 亚洲精华国产精华精| 国产人伦9x9x在线观看| 天堂俺去俺来也www色官网| 91麻豆精品激情在线观看国产 | 老司机靠b影院| 天堂动漫精品| 久久久久久亚洲精品国产蜜桃av| 欧美乱码精品一区二区三区| 精品一品国产午夜福利视频| 老鸭窝网址在线观看| 校园春色视频在线观看| 欧美性长视频在线观看| 精品人妻在线不人妻| 亚洲aⅴ乱码一区二区在线播放 | 亚洲av日韩精品久久久久久密| 国产成人一区二区三区免费视频网站| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产欧美一区二区综合| 日韩欧美在线二视频| 老司机亚洲免费影院| 亚洲五月色婷婷综合| 男人舔女人下体高潮全视频| 国产精品爽爽va在线观看网站 | 久久天堂一区二区三区四区| 成人手机av| 脱女人内裤的视频| 国产精品久久久av美女十八| 中国美女看黄片| 亚洲欧洲精品一区二区精品久久久| 亚洲欧美日韩另类电影网站| 国产黄a三级三级三级人| av福利片在线| 欧美乱色亚洲激情| 久久精品aⅴ一区二区三区四区| av网站免费在线观看视频| 一级片'在线观看视频| 欧美人与性动交α欧美软件| 精品国产美女av久久久久小说| 久久精品91无色码中文字幕| 可以在线观看毛片的网站| 天天躁夜夜躁狠狠躁躁| 俄罗斯特黄特色一大片| 多毛熟女@视频| 又紧又爽又黄一区二区| 色综合站精品国产| 国产精品久久电影中文字幕| 亚洲一区二区三区色噜噜 | 黄色丝袜av网址大全| 亚洲欧美日韩另类电影网站| 成人亚洲精品一区在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 激情在线观看视频在线高清| 国产精品免费一区二区三区在线| 日韩一卡2卡3卡4卡2021年| 国产熟女xx| av在线播放免费不卡| 亚洲七黄色美女视频| av网站在线播放免费| 久久精品亚洲熟妇少妇任你| avwww免费| 午夜福利在线观看吧| 国产高清国产精品国产三级| 超碰97精品在线观看| 亚洲精品美女久久久久99蜜臀| 99香蕉大伊视频| 久久国产亚洲av麻豆专区| 久久久久亚洲av毛片大全| 日本a在线网址| 日本三级黄在线观看| 亚洲欧美一区二区三区久久| 他把我摸到了高潮在线观看| 一级作爱视频免费观看| 又黄又粗又硬又大视频| 精品人妻在线不人妻| 午夜免费激情av| 亚洲色图综合在线观看| 日韩中文字幕欧美一区二区| 欧美一区二区精品小视频在线| 每晚都被弄得嗷嗷叫到高潮| 国产成人系列免费观看| 久久亚洲精品不卡| 国产精品久久久人人做人人爽| 精品久久久久久久毛片微露脸| 久久久久久久精品吃奶| 午夜精品在线福利| 可以免费在线观看a视频的电影网站| 精品久久久久久久久久免费视频 | 亚洲av成人不卡在线观看播放网| 国产欧美日韩一区二区三| 国产有黄有色有爽视频| 中文欧美无线码| 精品日产1卡2卡| 久久青草综合色| 变态另类成人亚洲欧美熟女 | 乱人伦中国视频| 国产一区二区激情短视频| 久久久国产成人精品二区 | 精品熟女少妇八av免费久了| 色老头精品视频在线观看| 国产免费av片在线观看野外av| 成年人黄色毛片网站| 色综合婷婷激情| av在线播放免费不卡| 国产成人影院久久av| 日日夜夜操网爽| 国产av在哪里看| 亚洲国产精品一区二区三区在线| 亚洲 欧美一区二区三区| 波多野结衣av一区二区av| 日韩av在线大香蕉| 日日爽夜夜爽网站| avwww免费| 在线观看免费日韩欧美大片| 色婷婷久久久亚洲欧美| 日韩欧美一区二区三区在线观看| 久久久精品国产亚洲av高清涩受| 黑人巨大精品欧美一区二区mp4| 久久久水蜜桃国产精品网| 99久久精品国产亚洲精品| 9热在线视频观看99| 18禁国产床啪视频网站| 无人区码免费观看不卡| 国产精品99久久99久久久不卡| 亚洲av成人av| 国内久久婷婷六月综合欲色啪| 午夜精品国产一区二区电影| 国产精品久久视频播放| 在线观看舔阴道视频| 国产主播在线观看一区二区| 久9热在线精品视频| 好看av亚洲va欧美ⅴa在| 精品国产一区二区三区四区第35| 后天国语完整版免费观看| 99久久99久久久精品蜜桃| 在线观看一区二区三区| 天天影视国产精品| a在线观看视频网站| 黄色毛片三级朝国网站| 在线观看一区二区三区激情| 在线观看午夜福利视频| 欧美人与性动交α欧美软件| 天堂动漫精品| 免费久久久久久久精品成人欧美视频| 亚洲精品一区av在线观看| 亚洲黑人精品在线| 欧美最黄视频在线播放免费 | 国产成人一区二区三区免费视频网站| 亚洲午夜精品一区,二区,三区| 国产欧美日韩一区二区三| 精品久久久久久电影网| 日本黄色视频三级网站网址| 香蕉丝袜av| 美女扒开内裤让男人捅视频| 伊人久久大香线蕉亚洲五| 人成视频在线观看免费观看| 黑人操中国人逼视频| 一进一出好大好爽视频| 久久午夜亚洲精品久久| 一边摸一边抽搐一进一出视频| 免费在线观看完整版高清| 免费不卡黄色视频| 动漫黄色视频在线观看| 69av精品久久久久久| 亚洲 欧美一区二区三区| 欧美在线黄色| 免费在线观看亚洲国产| 一级a爱视频在线免费观看| www.999成人在线观看| 女人精品久久久久毛片| 成人av一区二区三区在线看| 人成视频在线观看免费观看| 日日摸夜夜添夜夜添小说| 午夜a级毛片| 在线观看免费午夜福利视频| 亚洲欧美激情在线| 日韩av在线大香蕉| 久久久久久久久中文| 国产av一区二区精品久久| 黄色片一级片一级黄色片| 美女国产高潮福利片在线看| 国产精品爽爽va在线观看网站 | a在线观看视频网站| 欧美乱码精品一区二区三区| 国产精品一区二区三区四区久久 | 精品国产乱码久久久久久男人| 久久久久精品国产欧美久久久| 久久人人97超碰香蕉20202| cao死你这个sao货| 一本综合久久免费| 亚洲精品美女久久av网站| 久久久国产一区二区| 亚洲专区中文字幕在线| 一区二区三区精品91| 高清欧美精品videossex| 天堂动漫精品| 日本三级黄在线观看| 乱人伦中国视频| 99riav亚洲国产免费| 超碰成人久久| 亚洲av日韩精品久久久久久密| 黄色丝袜av网址大全| 波多野结衣av一区二区av| 波多野结衣一区麻豆| 美女国产高潮福利片在线看| 欧美日韩亚洲综合一区二区三区_| 国产精品一区二区免费欧美| 久久精品国产亚洲av香蕉五月| 国产国语露脸激情在线看| 又紧又爽又黄一区二区| 久久久久久久精品吃奶| 一本大道久久a久久精品| 十八禁网站免费在线| 免费久久久久久久精品成人欧美视频| 99国产极品粉嫩在线观看| 久久久久国内视频| 欧美日韩瑟瑟在线播放| 久久人妻福利社区极品人妻图片| 精品高清国产在线一区| 在线天堂中文资源库| 法律面前人人平等表现在哪些方面| 午夜精品国产一区二区电影| 久久性视频一级片| 国产精品免费视频内射| 欧美在线一区亚洲| 国产精品久久久久成人av| 国产一区二区三区综合在线观看| 久久人妻av系列| 波多野结衣av一区二区av| bbb黄色大片| 亚洲专区国产一区二区| 亚洲av第一区精品v没综合| 免费一级毛片在线播放高清视频 | 又大又爽又粗| 国产精品99久久99久久久不卡| 久久人妻福利社区极品人妻图片| 久久久久久久久中文| 亚洲男人天堂网一区| 高清av免费在线| 日本五十路高清| 国产欧美日韩一区二区精品| 欧美亚洲日本最大视频资源| 亚洲一区二区三区欧美精品| av电影中文网址| 国产欧美日韩一区二区三| 久9热在线精品视频| 欧美久久黑人一区二区| 国产精品一区二区精品视频观看| av在线播放免费不卡| 极品教师在线免费播放| 很黄的视频免费| 9色porny在线观看| 亚洲熟女毛片儿| 美女高潮到喷水免费观看| 99久久国产精品久久久| 一本综合久久免费| 国产亚洲精品第一综合不卡| 99国产精品一区二区蜜桃av| 亚洲国产中文字幕在线视频| 亚洲人成网站在线播放欧美日韩| 亚洲成av片中文字幕在线观看| 欧美成狂野欧美在线观看| 黄色怎么调成土黄色| 嫩草影视91久久| 久久中文字幕人妻熟女| 亚洲欧美精品综合久久99| 美女国产高潮福利片在线看| 岛国视频午夜一区免费看| 日韩有码中文字幕| 男男h啪啪无遮挡| 日韩欧美国产一区二区入口| 免费搜索国产男女视频| 级片在线观看| 午夜精品国产一区二区电影| 天堂中文最新版在线下载| 国产激情久久老熟女| 美女国产高潮福利片在线看| 国产精品一区二区三区四区久久 | 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品一二三| 99香蕉大伊视频| 欧美乱色亚洲激情| 两个人免费观看高清视频| 成人黄色视频免费在线看| www日本在线高清视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av成人一区二区三| www日本在线高清视频| 国产有黄有色有爽视频| 日韩国内少妇激情av| 国产成人精品在线电影| 欧美日韩国产mv在线观看视频| 日韩三级视频一区二区三区| 麻豆一二三区av精品| 亚洲国产毛片av蜜桃av| 99在线视频只有这里精品首页| 国产精品爽爽va在线观看网站 | 午夜免费激情av| 久久精品国产综合久久久| 一进一出好大好爽视频| 在线av久久热| 每晚都被弄得嗷嗷叫到高潮| 黑丝袜美女国产一区| 日韩一卡2卡3卡4卡2021年| 亚洲色图av天堂| 国产亚洲欧美在线一区二区| 国产单亲对白刺激| 亚洲熟女毛片儿| 大陆偷拍与自拍| 国产97色在线日韩免费| 亚洲欧美日韩另类电影网站| 国产一区二区在线av高清观看| 亚洲欧美激情综合另类| 国产精品98久久久久久宅男小说| 亚洲成av片中文字幕在线观看| 黄色片一级片一级黄色片| 丁香六月欧美| 一区二区三区国产精品乱码| netflix在线观看网站| 亚洲国产精品999在线| 欧美色视频一区免费| 欧美精品啪啪一区二区三区| 午夜福利,免费看| 美女 人体艺术 gogo| 午夜两性在线视频| 69av精品久久久久久| 9色porny在线观看| 波多野结衣高清无吗| 日韩av在线大香蕉| 欧美中文日本在线观看视频| 国产亚洲精品久久久久久毛片| 国产一区二区在线av高清观看| 久久久久国产一级毛片高清牌| 老司机靠b影院| 免费女性裸体啪啪无遮挡网站| 校园春色视频在线观看| 国产成人欧美| 亚洲人成77777在线视频| 国产精品久久电影中文字幕| 999精品在线视频| 国产成人精品在线电影| 国产激情久久老熟女| 少妇粗大呻吟视频| 午夜亚洲福利在线播放| 久久人人爽av亚洲精品天堂| 欧美在线黄色| 熟女少妇亚洲综合色aaa.| 久久精品国产亚洲av香蕉五月| 欧美中文综合在线视频| 精品一品国产午夜福利视频| 黑人猛操日本美女一级片| av欧美777| 国产精品 欧美亚洲| 两个人免费观看高清视频| 午夜久久久在线观看| 免费在线观看日本一区| 日韩有码中文字幕| 一二三四社区在线视频社区8| 国产黄a三级三级三级人| 久久香蕉激情| 丁香六月欧美| 老汉色∧v一级毛片| 午夜免费观看网址| 久久天躁狠狠躁夜夜2o2o| 亚洲avbb在线观看| 国产亚洲av高清不卡| 1024香蕉在线观看| 亚洲色图 男人天堂 中文字幕| av天堂久久9| 久久精品91蜜桃| 99国产精品免费福利视频| 亚洲色图 男人天堂 中文字幕| 淫妇啪啪啪对白视频| 91字幕亚洲| 免费高清在线观看日韩| 免费在线观看黄色视频的| 亚洲精品av麻豆狂野| 满18在线观看网站| 国产精华一区二区三区| 国产日韩一区二区三区精品不卡| 亚洲精品在线美女| 国产欧美日韩精品亚洲av| 日韩一卡2卡3卡4卡2021年| 亚洲专区字幕在线| 国产成人系列免费观看| 香蕉丝袜av| 国产精品 欧美亚洲| 又紧又爽又黄一区二区| 黄片小视频在线播放| 99riav亚洲国产免费| 国产片内射在线| 亚洲一区二区三区色噜噜 | 午夜成年电影在线免费观看| 人妻丰满熟妇av一区二区三区| 欧美性长视频在线观看| 国产精品免费视频内射| 一区二区三区激情视频| 成人国产一区最新在线观看| 精品电影一区二区在线| www.自偷自拍.com| 新久久久久国产一级毛片| 精品乱码久久久久久99久播| 国产成人影院久久av| 成人亚洲精品av一区二区 | 99国产精品99久久久久| 好看av亚洲va欧美ⅴa在| 岛国在线观看网站| 美女国产高潮福利片在线看| 好看av亚洲va欧美ⅴa在| 免费一级毛片在线播放高清视频 | 18禁裸乳无遮挡免费网站照片 | 三级毛片av免费| 丰满迷人的少妇在线观看| 精品国产超薄肉色丝袜足j| 18禁国产床啪视频网站| 午夜免费激情av| 日日摸夜夜添夜夜添小说| 亚洲精品在线美女| 99在线人妻在线中文字幕| 操出白浆在线播放| 一本大道久久a久久精品| 久久精品aⅴ一区二区三区四区| 一个人观看的视频www高清免费观看 | 国产精品99久久99久久久不卡| 精品一品国产午夜福利视频| 亚洲一区二区三区欧美精品| 在线观看免费高清a一片| 婷婷精品国产亚洲av在线| 亚洲精品中文字幕一二三四区| 中出人妻视频一区二区| 久久香蕉精品热| 国产精品爽爽va在线观看网站 | 老司机福利观看| 亚洲熟女毛片儿| 国产成人精品在线电影| 成人亚洲精品av一区二区 | 欧美日韩乱码在线| 中出人妻视频一区二区| 97人妻天天添夜夜摸| 午夜免费观看网址| 99精品久久久久人妻精品| 久久精品国产综合久久久| 国产真人三级小视频在线观看| 国产黄色免费在线视频| 人人澡人人妻人| 午夜免费成人在线视频| 国产成人精品久久二区二区91| 法律面前人人平等表现在哪些方面| 国产精品一区二区精品视频观看| 日本黄色日本黄色录像| 亚洲一卡2卡3卡4卡5卡精品中文| 精品国产乱子伦一区二区三区| 一边摸一边抽搐一进一出视频| 久久欧美精品欧美久久欧美| 岛国视频午夜一区免费看| 少妇裸体淫交视频免费看高清 | 国产一区在线观看成人免费| 国产精品二区激情视频| 啪啪无遮挡十八禁网站| 99久久精品国产亚洲精品| 亚洲欧美日韩高清在线视频| 国产免费男女视频| 国产高清视频在线播放一区| 少妇被粗大的猛进出69影院| 国产成人精品在线电影| 电影成人av| 老司机午夜十八禁免费视频| 久久久久国产一级毛片高清牌| 亚洲情色 制服丝袜| 不卡av一区二区三区| 12—13女人毛片做爰片一| 国产一区在线观看成人免费| 久久久水蜜桃国产精品网| 亚洲一区二区三区欧美精品| 十分钟在线观看高清视频www| 久久影院123| 9热在线视频观看99| 波多野结衣高清无吗| 88av欧美| 久99久视频精品免费| 热re99久久国产66热| 一a级毛片在线观看| 大香蕉久久成人网| 丁香六月欧美| 国产精品国产av在线观看| 国产成人av教育| 一级a爱视频在线免费观看| 90打野战视频偷拍视频| 人妻丰满熟妇av一区二区三区| 美女大奶头视频| x7x7x7水蜜桃| 在线观看www视频免费| 18禁美女被吸乳视频| 午夜日韩欧美国产| 久久精品影院6| 日本三级黄在线观看| 亚洲第一av免费看| 夜夜夜夜夜久久久久| 精品人妻1区二区| 欧美成狂野欧美在线观看| 亚洲av五月六月丁香网| 少妇裸体淫交视频免费看高清 | 高清毛片免费观看视频网站 | 欧美日本亚洲视频在线播放| 制服诱惑二区| 久久国产精品影院| 亚洲av熟女| 91在线观看av| 国产免费现黄频在线看| 久久青草综合色| 精品乱码久久久久久99久播| 在线播放国产精品三级| 午夜a级毛片| www国产在线视频色| 亚洲精品一区av在线观看| 欧美日韩一级在线毛片| 动漫黄色视频在线观看| 男人的好看免费观看在线视频 | 国产精品国产高清国产av| 淫秽高清视频在线观看| 国产99久久九九免费精品| 国产精品二区激情视频| 精品第一国产精品| 日韩欧美一区二区三区在线观看| www.精华液| 国产一区在线观看成人免费| 久久精品国产99精品国产亚洲性色 | 久久久久九九精品影院| 在线观看66精品国产| 天堂中文最新版在线下载| 国产免费男女视频| 久久人妻熟女aⅴ| 超碰成人久久| 精品免费久久久久久久清纯| 在线观看舔阴道视频| 亚洲性夜色夜夜综合| 91在线观看av| 午夜福利在线免费观看网站| 色综合婷婷激情| 午夜福利免费观看在线| 人人妻人人添人人爽欧美一区卜| 亚洲视频免费观看视频| 伦理电影免费视频| 亚洲精品一二三| 国产日韩一区二区三区精品不卡| a级片在线免费高清观看视频| 久久久久精品国产欧美久久久| av视频免费观看在线观看| 黄色丝袜av网址大全| 91麻豆精品激情在线观看国产 | 巨乳人妻的诱惑在线观看| 久久久久国产一级毛片高清牌|