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

    軟基水閘底板脫空反演的響應(yīng)面模型優(yōu)化

    2022-01-07 05:59:50劉雙平柯賢勇李火坤李輝唐義員涂源方靜
    關(guān)鍵詞:水閘反演底板

    劉雙平,柯賢勇,李火坤,李輝,唐義員,涂源,方靜

    (南昌大學(xué)建筑工程學(xué)院,江西 南昌 330031)

    長期處于激振交替、動靜懸殊的工作環(huán)境中的軟基水閘容易出現(xiàn)底板脫空[1],目前,基于動力學(xué)特性對水閘結(jié)構(gòu)進(jìn)行反分析被廣泛研究[2-4]。尋找準(zhǔn)確可靠的響應(yīng)面模型代替有限元模型,是進(jìn)行軟基水閘底板脫空反演的關(guān)鍵。

    響應(yīng)面方程是決定脫空參數(shù)與模態(tài)參數(shù)之間非線性關(guān)系的重要因素,當(dāng)響應(yīng)面方程過于簡單時,響應(yīng)面模型難以保證其精度;當(dāng)響應(yīng)面方程過于復(fù)雜時,響應(yīng)面模型的效率會降低[5]。于景飛等[6]采用線性多項式函數(shù)構(gòu)建響應(yīng)面方程,對大跨徑混凝土斜拉橋進(jìn)行了優(yōu)化預(yù)測研究;Fang等[7]指出參數(shù)間相互效應(yīng)對響應(yīng)面模型總方差的貢獻(xiàn)非常小;張運濤等[8]基于逐步回歸的方法對變量的選取進(jìn)行顯著性分析,提高了所選取不含交叉項二次多項式響應(yīng)面模型的精度和計算效率,并將其運用于大跨連續(xù)剛構(gòu)橋長期變形預(yù)測;李延強(qiáng)等[9]采用F檢驗法進(jìn)行參數(shù)的顯著性分析,以二次多項式擬合響應(yīng)面模型,對比了考慮交叉項和不考慮交叉項對斜拉橋模型的修正結(jié)果;宗周紅、劉杰等[10-11]采用高階響應(yīng)面模型應(yīng)用于橋梁結(jié)構(gòu)的損傷識別。張建偉等[12]采用基于方差分析的F檢驗法,對響應(yīng)面模型參數(shù)進(jìn)行顯著性檢驗,并將優(yōu)化后的響應(yīng)面模型應(yīng)用于某大型渡槽結(jié)構(gòu)的有限元模型修正;蒙偉等[13]分析了F檢驗和t檢驗的特點,采用了一種能直觀方便地判斷回歸模型和回歸系數(shù)顯著性的p值法檢驗,并將其應(yīng)用于初始地應(yīng)力場的反演分析。

    為獲得軟基水閘底板脫空反演中最佳的響應(yīng)面模型,本文提出了一種基于顯著性檢驗的響應(yīng)面模型優(yōu)選與優(yōu)化方法。首先,采用拉丁超立方抽樣(LHS)獲得脫空參數(shù)樣本集,并將其輸入水閘有限元模型計算相應(yīng)的模態(tài)參數(shù)集,以此構(gòu)建表征脫空參數(shù)與模態(tài)參數(shù)之間非線性關(guān)系的多種不同形式的響應(yīng)面模型;其次,對不同形式的響應(yīng)面模型進(jìn)行F檢驗,綜合考慮精度和效率,優(yōu)選形式最佳的響應(yīng)面模型;然后,對形式最佳的響應(yīng)面模型進(jìn)行t檢驗,優(yōu)化響應(yīng)面模型的參數(shù),剔除對模態(tài)參數(shù)影響不足的脫空參數(shù)及其組合(贅余項)以獲得最佳的響應(yīng)面模型;最后,將最佳的響應(yīng)面模型應(yīng)用于軟基水閘物理模型的底板脫空反演,驗證方法的有效性。

    1 響應(yīng)面模型優(yōu)選與優(yōu)化

    1.1 拉丁超立方抽樣試驗設(shè)計

    試驗樣本的選取方法與數(shù)量關(guān)系到響應(yīng)面模型的精度和效率,以及分析結(jié)果的準(zhǔn)確性[14]。采用拉丁超立方抽樣試驗設(shè)計,能夠以少量的脫空參數(shù)樣本,全面而準(zhǔn)確地反映軟基水閘底板脫空的分布情況。采用LHS生成脫空參數(shù)樣本集思路如下[15-16]:

    然后,根據(jù)脫空參數(shù)的概率分布情況,對每一個子區(qū)間進(jìn)行獨立的等概率抽樣,第j個子區(qū)間抽樣獲取的脫空參數(shù)值dij滿足以下關(guān)系:

    (1)

    式中:dij表示第i個脫空參數(shù)在第j個子區(qū)間的抽樣值;r表示服從[0,1]均勻分布的隨機(jī)數(shù);m表示脫空參數(shù)的個數(shù),即抽樣維度;n表示子區(qū)間數(shù),即抽樣組數(shù)。

    最后,將各子區(qū)間生成的樣本點匯總成樣本集,對各脫空參數(shù)的抽樣值隨機(jī)組合形成多維脫空參數(shù)樣本集。

    1.2 響應(yīng)面模型

    有限元模型結(jié)構(gòu)復(fù)雜,單元龐大,難以進(jìn)行反分析計算。合理構(gòu)建脫空參數(shù)與模態(tài)參數(shù)(固有頻率和節(jié)點振型)之間的響應(yīng)面模型代替有限元模型,是進(jìn)行軟基水閘底板脫空反演的基礎(chǔ)。本文采用的響應(yīng)面方程為基于Taylor展開的多項式函數(shù),響應(yīng)面方程如下:

    (2)

    (3)

    (4)

    (5)

    (6)

    (7)

    (8)

    (9)

    式中:α表示頻率響應(yīng)面方程的待定系數(shù)。

    響應(yīng)面模型的構(gòu)建思路:首先,通過LHS獲取脫空參數(shù)樣本集;其次,將脫空參數(shù)樣本集輸入有限元模型中,計算相應(yīng)的模態(tài)參數(shù)集;然后,根據(jù)脫空參數(shù)集和模態(tài)參數(shù)集,以上述7種不同形式的多項式函數(shù)構(gòu)建響應(yīng)面方程組;最后,根據(jù)多元回歸分析確定響應(yīng)面方程組的待定系數(shù),從而建立7種不同形式的響應(yīng)面模型。

    1.3 顯著性檢驗

    為獲得軟基水閘底板脫空反演中最佳的響應(yīng)面模型,采用顯著性檢驗對響應(yīng)面模型進(jìn)行優(yōu)選與優(yōu)化。首先,基于回歸模型的顯著性檢驗(F檢驗),對上述7種不同形式的響應(yīng)面模型進(jìn)行顯著性分析,優(yōu)選形式最佳的響應(yīng)面模型;隨后,基于回歸系數(shù)的顯著性檢驗(t檢驗),對形式最佳的響應(yīng)面模型的回歸系數(shù)進(jìn)行顯著性分析,優(yōu)化響應(yīng)面模型的參數(shù),得到最佳的響應(yīng)面模型。

    1.3.1F檢驗

    H0:αi=αii=αiii=αiiii=αij=αijk=αiij=0

    (10)

    式中:H0表示頻率響應(yīng)面模型的檢驗狀態(tài)。

    利用方差分析法,計算上述7種響應(yīng)面方程的F統(tǒng)計量,公式如下:

    (11)

    (12)

    (13)

    (14)

    1.3.2t檢驗

    H0:αI=0,I∈{i,ii,iii,iiii,ij,ijk,iij}

    (15)

    式中:H0表示頻率響應(yīng)面模型回歸系數(shù)的檢驗狀態(tài)。

    可由以下公式計算響應(yīng)面方程中每一項回歸系數(shù)的TI值:

    (16)

    (17)

    (18)

    (19)

    1.4 響應(yīng)面模型評價指標(biāo)

    在響應(yīng)面方程形式持續(xù)復(fù)雜化的過程中,響應(yīng)面模型的精度必然會一步步提高,與此同時,擬合響應(yīng)面以及進(jìn)行軟基水閘底板脫空反演的效率必然會有所降低。因此,在優(yōu)選形式最佳的響應(yīng)面方程時,需綜合考慮精度和效率兩方面的影響。本文采用3個精度評價指標(biāo)和1個效率評價指標(biāo)來衡量響應(yīng)面模型的優(yōu)劣。

    1.4.1 精度評價指標(biāo)

    (20)

    (21)

    相對均方根誤差RMSE′反映響應(yīng)面模型模態(tài)參數(shù)的輸出值與有限元模型模態(tài)參數(shù)的計算值之間的誤差,其表達(dá)式如下:

    (22)

    相對均方根誤差RMSE′絕對值(經(jīng)歸一化之后的振型存在負(fù)值)越小,表明響應(yīng)面模型模態(tài)參數(shù)的輸出值與有限元模型模態(tài)參數(shù)的計算值之間的誤差越小,響應(yīng)面模型精度越高。

    以上兩種評價指標(biāo)均表示響應(yīng)面模型的整體精度,而每一組脫空參數(shù)樣本點對應(yīng)的響應(yīng)面模型模態(tài)參數(shù)的輸出值與有限元模型模態(tài)參數(shù)的計算值之間的相對誤差Rsp,能夠有效地評價響應(yīng)面模型點對點的擬合精度,其表達(dá)式如下:

    (23)

    對所有的脫空參數(shù)樣本點對應(yīng)的模態(tài)參數(shù)由式(23)遍歷,設(shè)定閾值為5‰[17]。對LHS的所有樣本點,超出閾值的個數(shù)越少,Rsp max越小,響應(yīng)面模型精度越高。

    1.4.2 效率評價指標(biāo)

    效率是檢驗響應(yīng)面模型的擬合過程以及在代替有限元模型進(jìn)行軟基水閘底板脫空反演過程中耗時是否經(jīng)濟(jì)合理的重要指標(biāo)。考慮到不同計算機(jī)之間的差異性,計算耗時并不適合直接作為效率評價指標(biāo)。在利用多元回歸分析擬合響應(yīng)面方程中,待定系數(shù)越多,擬合效率越低,此外,在進(jìn)行軟基水閘底板脫空反演中,響應(yīng)面模型的參數(shù)越多,反演效率越低,因此,以響應(yīng)面方程展開后的項數(shù)K作為響應(yīng)面模型的效率評價指標(biāo)是合理可行的。響應(yīng)面方程展開后的項數(shù)K越小,響應(yīng)面模型效率越高。

    1.5 響應(yīng)面模型優(yōu)化流程

    基于顯著性檢驗的軟基水閘底板脫空反演中響應(yīng)面模型優(yōu)選與優(yōu)化流程圖如圖1所示。

    圖1 軟基水閘底板脫空反演中響應(yīng)面模型優(yōu)選與優(yōu)化流程圖Fig.1 Flowchart of response surface model selection and optimization in the inversion for the floor void of sluice on the soft foundation

    2 軟基水閘物理模型驗證

    2.1 軟基水閘物理模型基本概況

    以某水閘工程為背景實例,制作一單孔軟基水閘室內(nèi)物理模型[4]。基于水閘物理模型,采用人工掏空的方式對水閘模型底板順?biāo)鞣较蜻M(jìn)行了單側(cè)型和雙側(cè)型隨機(jī)脫空模擬,并測量脫空縱深,如圖2所示。在描述底板脫空范圍的數(shù)學(xué)模型中[3],采用5個參數(shù)描述單側(cè)脫空(即m=5);10個參數(shù)描述雙側(cè)脫空(即m=10)。脫空模擬共分為3種工況:其中,工況一和工況二為單側(cè)脫空,工況三為雙側(cè)脫空。不同工況下的脫空參數(shù)di如表1所示。

    (a) 物理模型脫空

    表1 水閘物理模型底板脫空工況及脫空參數(shù)Tab.1 Void conditions and void parameters of the sluice physical model bottom plate

    2.2 構(gòu)建響應(yīng)面模型

    根據(jù)水閘物理模型尺寸,結(jié)合工程經(jīng)驗,將底板脫空的抽樣范圍確定為上游側(cè)[0,0.65]m,下游側(cè)[0,0.4]m。單側(cè)脫空和雙側(cè)脫空的抽樣維度分別為5和10,采用LHS生成2 000組脫空參數(shù)樣本。將脫空參數(shù)樣本集輸入水閘有限元基準(zhǔn)模型中,計算與之相應(yīng)的模態(tài)參數(shù),提取前四階頻率和振型,分別構(gòu)建脫空參數(shù)與頻率和振型之間的響應(yīng)面模型。

    2.3 優(yōu)選與優(yōu)化響應(yīng)面模型

    2.3.1 優(yōu)選響應(yīng)面模型的形式

    以雙側(cè)脫空,頻率響應(yīng)面模型為例,對響應(yīng)面模型進(jìn)行F檢驗,F(xiàn)統(tǒng)計量及其臨界值如表2所示,可知,7種響應(yīng)面方程的F統(tǒng)計量遠(yuǎn)大于臨界值,表明7種響應(yīng)面模型均能準(zhǔn)確有效地代替有限元模型。此外,第三階響應(yīng)面模型的F統(tǒng)計量比其余三階都小,表明第三階模態(tài)參數(shù)樣本變異性更大,對脫空參數(shù)的敏感性更高。

    表2 7種頻率響應(yīng)面模型的F統(tǒng)計量及其臨界值(雙側(cè)脫空)Tab.2 F statistics of seven frequency response surface models and their critical values (two-sided void)

    表3 7種頻率響應(yīng)面模型的整體精度(雙側(cè)脫空)Tab.3 Overall accuracy of seven frequency response surface models (two-sided void)

    表4 7種頻率響應(yīng)面模型的樣本點精度(雙側(cè)脫空)Tab.4 Sample point accuracy of seven frequency response surface models (two-sided void)

    表5 7種響應(yīng)面模型的效率Tab.5 Efficiency of seven response surface models

    2.3.2 優(yōu)化響應(yīng)面模型的參數(shù)

    (a) 第一階

    根據(jù)圖3可以直觀地看出對頻率響應(yīng)影響不足的脫空參數(shù)及其組合,將各階頻率響應(yīng)面方程中的贅余項剔除,常數(shù)項及保留項所構(gòu)成的響應(yīng)面方程為最佳的響應(yīng)面方程,經(jīng)參數(shù)優(yōu)化的響應(yīng)面模型即為最佳的響應(yīng)面模型。

    2.4 結(jié)果驗證

    為驗證經(jīng)顯著性檢驗的響應(yīng)面模型的優(yōu)越性,將最佳的響應(yīng)面模型應(yīng)用于軟基水閘底板脫空反演,并將反演結(jié)果與參考文獻(xiàn)[4]對比。不同工況下水閘底板脫空控制參數(shù)反演值與實際值如表6所示、識別的脫空面積及趨勢與實際脫空如圖4所示、識別的脫空面積與實際的脫空面積之間的相對誤差如表7所示,同時,為更準(zhǔn)確衡量反演結(jié)果的準(zhǔn)確性,本文提出了一種更為科學(xué)合理的評價指標(biāo)——面積不重合度,即識別脫空與實際脫空不重合部分的面積在實際脫空的面積占比。從反演識別結(jié)果來看,識別脫空與實際脫空能較好地吻合,3種工況下水閘底板脫空反演結(jié)果與實際脫空的面積不重合度分別為14.79%,5.46%,7.36%,相比參考文獻(xiàn)[4]的反演結(jié)果有了較大的提升,驗證了本文方法的有效性。

    表6 不同工況下脫空參數(shù)反演結(jié)果Tab.6 Inversion results of void parameters under different conditions

    (a) 工況一

    表7 識別脫空面積與實際脫空面積的相對誤差Tab.7 Relative error between the identify void area and the actual void area

    3 結(jié)論

    (1) 基于顯著性分析,采用F檢驗,從7種不同形式的響應(yīng)面模型中優(yōu)選出形式最佳的響應(yīng)面模型;隨后,對形式最佳的響應(yīng)面模型的脫空參數(shù)及其組合進(jìn)行t檢驗以刪除贅余項,建立了能夠準(zhǔn)確高效地代替有限元模型的響應(yīng)面模型。

    (2) 提出了采用面積不重合度來衡量反演結(jié)果的準(zhǔn)確性,使得軟基水閘底板脫空識別的精度評價更加科學(xué)合理。

    (3) 將經(jīng)顯著性檢驗獲得的最佳的響應(yīng)面模型應(yīng)用于軟基水閘底板脫空反演,結(jié)果表明:識別脫空與實際脫空能夠較好地吻合,3種工況下的水閘底板脫空反演結(jié)果與實際脫空的面積不重合度分別為14.79%,5.46%,7.36%,提高了反演精度,驗證了本文方法的有效性。

    猜你喜歡
    水閘反演底板
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    2016年河南省已建成水閘數(shù)量
    全省已建成水閘數(shù)量
    軟土地基上超長底板節(jié)制閘沉降及應(yīng)力分析
    底板巷一巷兩用的回采工作面防火實踐
    底板隔水層破壞的力學(xué)模型及破壞判據(jù)
    河南省2014年已建成水閘數(shù)量
    疊前同步反演在港中油田的應(yīng)用
    亚洲国产精品国产精品| 欧美色欧美亚洲另类二区| 国产成年人精品一区二区| 国产免费一级a男人的天堂| av在线亚洲专区| 午夜精品一区二区三区免费看| 在线免费观看不下载黄p国产| 国产一区二区三区在线臀色熟女| 欧美日本亚洲视频在线播放| 一个人看的www免费观看视频| 欧美高清性xxxxhd video| 亚洲av免费高清在线观看| 久久久久久久午夜电影| 日本成人三级电影网站| 欧美成人免费av一区二区三区| 亚洲高清免费不卡视频| 免费电影在线观看免费观看| 成年版毛片免费区| a级一级毛片免费在线观看| 国产黄色视频一区二区在线观看 | av天堂中文字幕网| 国产精品乱码一区二三区的特点| 成人美女网站在线观看视频| 亚洲国产精品国产精品| 丝袜喷水一区| 国产成人aa在线观看| 亚洲国产精品成人综合色| 亚洲18禁久久av| eeuss影院久久| 精品熟女少妇av免费看| 91在线观看av| 免费人成在线观看视频色| 国产在视频线在精品| 日本欧美国产在线视频| 一个人免费在线观看电影| 免费电影在线观看免费观看| 搡老岳熟女国产| 国产精品久久久久久亚洲av鲁大| av.在线天堂| 亚洲精品乱码久久久v下载方式| 色播亚洲综合网| 亚洲18禁久久av| 成人国产麻豆网| 插逼视频在线观看| 一级毛片aaaaaa免费看小| 一本一本综合久久| 国产国拍精品亚洲av在线观看| av国产免费在线观看| 国产午夜福利久久久久久| 久久久久免费精品人妻一区二区| 精品久久久久久久人妻蜜臀av| 国产乱人偷精品视频| 女的被弄到高潮叫床怎么办| 成年av动漫网址| 色视频www国产| 禁无遮挡网站| 国产高清有码在线观看视频| 免费高清视频大片| 一级a爱片免费观看的视频| 91av网一区二区| 亚洲国产高清在线一区二区三| 亚洲,欧美,日韩| 青春草视频在线免费观看| 欧美另类亚洲清纯唯美| 能在线免费观看的黄片| 毛片一级片免费看久久久久| 九色成人免费人妻av| 菩萨蛮人人尽说江南好唐韦庄 | av在线观看视频网站免费| 身体一侧抽搐| 亚洲中文字幕一区二区三区有码在线看| 欧美日韩国产亚洲二区| 午夜视频国产福利| 老女人水多毛片| 性欧美人与动物交配| 一本一本综合久久| 嫩草影院入口| 一夜夜www| 天堂√8在线中文| av在线观看视频网站免费| 国产精品综合久久久久久久免费| 在线观看免费视频日本深夜| 国语自产精品视频在线第100页| 联通29元200g的流量卡| 一卡2卡三卡四卡精品乱码亚洲| 桃色一区二区三区在线观看| 欧美成人一区二区免费高清观看| 狠狠狠狠99中文字幕| 亚洲av成人av| 精品乱码久久久久久99久播| 在线免费观看的www视频| 日韩欧美三级三区| 久久99热这里只有精品18| 久久精品国产亚洲av香蕉五月| 日韩中字成人| 91av网一区二区| 免费在线观看影片大全网站| 一个人看的www免费观看视频| 国产伦一二天堂av在线观看| 欧美日韩国产亚洲二区| 亚洲中文日韩欧美视频| 精品久久久久久久久久免费视频| 久久久精品94久久精品| 色哟哟·www| 欧美精品国产亚洲| 天堂√8在线中文| 免费av毛片视频| 男女之事视频高清在线观看| 成人特级av手机在线观看| 舔av片在线| 亚洲欧美日韩卡通动漫| 变态另类丝袜制服| 久久99热这里只有精品18| 久久这里只有精品中国| 村上凉子中文字幕在线| 人人妻,人人澡人人爽秒播| 色哟哟哟哟哟哟| 国产视频一区二区在线看| 人人妻人人澡人人爽人人夜夜 | 精品人妻熟女av久视频| 日本免费a在线| av在线天堂中文字幕| 九九热线精品视视频播放| 男女视频在线观看网站免费| 在线观看午夜福利视频| 精品欧美国产一区二区三| 国产亚洲av嫩草精品影院| 国产精品一区www在线观看| 国产视频内射| 日韩欧美国产在线观看| 别揉我奶头 嗯啊视频| 亚洲在线自拍视频| 欧美高清成人免费视频www| 亚洲最大成人中文| .国产精品久久| 亚洲电影在线观看av| 97碰自拍视频| 黄片wwwwww| 深夜a级毛片| 一a级毛片在线观看| 久久精品久久久久久噜噜老黄 | 97超级碰碰碰精品色视频在线观看| 久久99热6这里只有精品| 成年av动漫网址| 97超级碰碰碰精品色视频在线观看| 亚洲av不卡在线观看| 天堂网av新在线| 欧美极品一区二区三区四区| 亚洲成人久久爱视频| 夜夜夜夜夜久久久久| 最新在线观看一区二区三区| 丝袜喷水一区| 亚洲性久久影院| 午夜福利在线观看免费完整高清在 | 午夜福利视频1000在线观看| 夜夜夜夜夜久久久久| 成人一区二区视频在线观看| 国产老妇女一区| 91av网一区二区| 国产精品综合久久久久久久免费| 欧美最新免费一区二区三区| 成人高潮视频无遮挡免费网站| 赤兔流量卡办理| 直男gayav资源| 日本一本二区三区精品| 中文字幕久久专区| 又黄又爽又刺激的免费视频.| 国产在线精品亚洲第一网站| 亚洲国产精品合色在线| 一个人看视频在线观看www免费| 看免费成人av毛片| 国产精品一区二区性色av| 一进一出抽搐gif免费好疼| 亚洲成人精品中文字幕电影| 国产蜜桃级精品一区二区三区| 国产精品一区二区三区四区久久| 中文字幕熟女人妻在线| 国产成人a∨麻豆精品| av卡一久久| 亚洲av免费高清在线观看| 亚洲成人久久爱视频| 亚洲av不卡在线观看| 十八禁国产超污无遮挡网站| 亚洲人成网站在线观看播放| 亚洲精华国产精华液的使用体验 | 国产精品一区二区免费欧美| 人妻丰满熟妇av一区二区三区| 国产伦在线观看视频一区| 99久久九九国产精品国产免费| av.在线天堂| av在线观看视频网站免费| 色吧在线观看| 一级毛片久久久久久久久女| 日本一二三区视频观看| 久久精品夜色国产| 久久草成人影院| 亚洲av成人精品一区久久| 真人做人爱边吃奶动态| 婷婷精品国产亚洲av| 免费在线观看影片大全网站| 少妇人妻一区二区三区视频| 国产精品av视频在线免费观看| 最新中文字幕久久久久| 波多野结衣高清作品| 一级毛片aaaaaa免费看小| 久久久久久久久中文| 精品福利观看| 神马国产精品三级电影在线观看| eeuss影院久久| 亚洲国产欧洲综合997久久,| 我要搜黄色片| 午夜免费激情av| 啦啦啦啦在线视频资源| 又黄又爽又刺激的免费视频.| 久久久久免费精品人妻一区二区| 亚洲国产高清在线一区二区三| 村上凉子中文字幕在线| 我的女老师完整版在线观看| 国国产精品蜜臀av免费| 别揉我奶头~嗯~啊~动态视频| 国产大屁股一区二区在线视频| 国产亚洲精品综合一区在线观看| 亚洲性久久影院| 精品不卡国产一区二区三区| 波多野结衣高清无吗| 我的女老师完整版在线观看| av在线天堂中文字幕| 一本一本综合久久| 国产v大片淫在线免费观看| 亚洲欧美日韩高清在线视频| 日韩国内少妇激情av| 12—13女人毛片做爰片一| 久久精品国产亚洲av香蕉五月| 一个人免费在线观看电影| 美女cb高潮喷水在线观看| 麻豆av噜噜一区二区三区| 久久久久精品国产欧美久久久| 久久亚洲国产成人精品v| 永久网站在线| 亚洲av中文av极速乱| 亚洲av.av天堂| 日本爱情动作片www.在线观看 | 天堂av国产一区二区熟女人妻| 99久久精品一区二区三区| 高清午夜精品一区二区三区 | 成人特级黄色片久久久久久久| 色哟哟哟哟哟哟| 黄片wwwwww| 亚洲精品国产成人久久av| 国模一区二区三区四区视频| 99热6这里只有精品| 国产精品久久电影中文字幕| 亚洲图色成人| 国产av麻豆久久久久久久| 婷婷亚洲欧美| 日韩av不卡免费在线播放| 蜜臀久久99精品久久宅男| 午夜a级毛片| 国产精品精品国产色婷婷| 日韩亚洲欧美综合| 国产精品久久久久久久久免| 淫妇啪啪啪对白视频| 婷婷亚洲欧美| 中文亚洲av片在线观看爽| 久久国产乱子免费精品| 国产精品精品国产色婷婷| 国产欧美日韩一区二区精品| 国产 一区 欧美 日韩| 久久亚洲国产成人精品v| 亚洲av二区三区四区| 亚洲精品在线观看二区| 天堂动漫精品| 岛国在线免费视频观看| 午夜精品在线福利| 欧美高清性xxxxhd video| 国产精品一区二区免费欧美| 看免费成人av毛片| 午夜福利在线在线| 天堂√8在线中文| av视频在线观看入口| 99热6这里只有精品| 成人美女网站在线观看视频| 白带黄色成豆腐渣| 久久精品夜夜夜夜夜久久蜜豆| 乱人视频在线观看| 插阴视频在线观看视频| 蜜桃亚洲精品一区二区三区| 天堂网av新在线| 欧美人与善性xxx| 不卡一级毛片| 欧美日韩精品成人综合77777| 久久精品影院6| 日本欧美国产在线视频| av在线老鸭窝| 国产老妇女一区| 国产伦精品一区二区三区视频9| 精品久久久久久久久av| 久久韩国三级中文字幕| 日韩 亚洲 欧美在线| 精品久久久久久久久久免费视频| 国内精品一区二区在线观看| 久久久午夜欧美精品| 国产成人91sexporn| 一级毛片电影观看 | 亚洲精品国产av成人精品 | 熟妇人妻久久中文字幕3abv| 日本成人三级电影网站| 国产男靠女视频免费网站| 亚洲精品色激情综合| 日日干狠狠操夜夜爽| 麻豆国产av国片精品| 久久精品国产99精品国产亚洲性色| 十八禁国产超污无遮挡网站| 精品少妇黑人巨大在线播放 | 看黄色毛片网站| 欧美三级亚洲精品| 91麻豆精品激情在线观看国产| 亚洲婷婷狠狠爱综合网| 精品人妻熟女av久视频| 精品久久久久久久久久免费视频| 在线观看免费视频日本深夜| 热99在线观看视频| 免费看美女性在线毛片视频| 国产中年淑女户外野战色| 成人综合一区亚洲| 欧美三级亚洲精品| 丰满人妻一区二区三区视频av| 午夜福利在线观看吧| 国产aⅴ精品一区二区三区波| av在线蜜桃| 国产一区二区亚洲精品在线观看| 18禁裸乳无遮挡免费网站照片| 免费在线观看成人毛片| 国产探花极品一区二区| 日韩中字成人| 亚洲av熟女| av在线老鸭窝| 22中文网久久字幕| 国产精品精品国产色婷婷| 极品教师在线视频| 看十八女毛片水多多多| 菩萨蛮人人尽说江南好唐韦庄 | 一区二区三区免费毛片| 中文字幕精品亚洲无线码一区| 国内精品久久久久精免费| 成人特级黄色片久久久久久久| 精品久久久久久久久亚洲| 男女做爰动态图高潮gif福利片| 国产视频内射| 色尼玛亚洲综合影院| av在线亚洲专区| 午夜福利18| 精品免费久久久久久久清纯| 久久草成人影院| 久久久久国产精品人妻aⅴ院| 美女大奶头视频| 久久综合国产亚洲精品| 中文字幕免费在线视频6| 黄色配什么色好看| 最好的美女福利视频网| 女人十人毛片免费观看3o分钟| 真实男女啪啪啪动态图| 国产精品综合久久久久久久免费| 丰满乱子伦码专区| 国产日本99.免费观看| 国产亚洲精品综合一区在线观看| 中文字幕av成人在线电影| 国产私拍福利视频在线观看| 美女xxoo啪啪120秒动态图| 秋霞在线观看毛片| 成人精品一区二区免费| 神马国产精品三级电影在线观看| 中文字幕熟女人妻在线| 99热全是精品| 亚洲,欧美,日韩| 无遮挡黄片免费观看| 乱人视频在线观看| 精品无人区乱码1区二区| 亚洲中文日韩欧美视频| 日韩欧美 国产精品| 免费高清视频大片| 欧美色视频一区免费| 亚洲美女搞黄在线观看 | 精品久久久久久久久亚洲| 天美传媒精品一区二区| 欧美色欧美亚洲另类二区| 国产成人a区在线观看| 日日摸夜夜添夜夜添小说| 九九在线视频观看精品| 99视频精品全部免费 在线| 搡女人真爽免费视频火全软件 | 成人鲁丝片一二三区免费| 日日摸夜夜添夜夜添小说| av中文乱码字幕在线| 中国美女看黄片| 午夜精品一区二区三区免费看| 久久欧美精品欧美久久欧美| 久久综合国产亚洲精品| 亚洲精品乱码久久久v下载方式| 国产成人精品久久久久久| 免费大片18禁| 黄色欧美视频在线观看| 亚洲真实伦在线观看| 日本a在线网址| 亚洲天堂国产精品一区在线| 男女啪啪激烈高潮av片| 国产黄色小视频在线观看| 国内久久婷婷六月综合欲色啪| 国产免费一级a男人的天堂| 深爱激情五月婷婷| 免费人成在线观看视频色| 精品熟女少妇av免费看| 久久久久久大精品| 午夜精品国产一区二区电影 | 床上黄色一级片| 国产单亲对白刺激| 全区人妻精品视频| 噜噜噜噜噜久久久久久91| 亚洲人成网站在线播| 精品一区二区三区人妻视频| 美女xxoo啪啪120秒动态图| 内地一区二区视频在线| 久久精品91蜜桃| 综合色av麻豆| 午夜精品国产一区二区电影 | 久久精品影院6| 免费在线观看成人毛片| 色综合色国产| 九九久久精品国产亚洲av麻豆| 亚洲欧美清纯卡通| 美女被艹到高潮喷水动态| 亚洲精品在线观看二区| 日韩欧美免费精品| 亚洲不卡免费看| 国产一级毛片七仙女欲春2| 麻豆国产97在线/欧美| 非洲黑人性xxxx精品又粗又长| 国产成人一区二区在线| 久久精品国产自在天天线| 99riav亚洲国产免费| 中文字幕久久专区| 亚洲精品一区av在线观看| 狠狠狠狠99中文字幕| 国产乱人偷精品视频| av在线蜜桃| 俄罗斯特黄特色一大片| 内地一区二区视频在线| 此物有八面人人有两片| 嫩草影院入口| 人人妻人人看人人澡| www日本黄色视频网| 美女内射精品一级片tv| 人妻丰满熟妇av一区二区三区| 伦精品一区二区三区| 日韩一本色道免费dvd| 18禁在线无遮挡免费观看视频 | 久久亚洲精品不卡| 51国产日韩欧美| 一级黄色大片毛片| 最近视频中文字幕2019在线8| 亚洲精品乱码久久久v下载方式| 久久精品人妻少妇| 久久天躁狠狠躁夜夜2o2o| a级一级毛片免费在线观看| avwww免费| 国产中年淑女户外野战色| 国产精品一区二区性色av| 观看美女的网站| 真实男女啪啪啪动态图| 日本与韩国留学比较| 久久久色成人| h日本视频在线播放| 国国产精品蜜臀av免费| 久久久久久久久久久丰满| 亚洲精品一区av在线观看| 久久久久国内视频| 亚洲欧美成人精品一区二区| 久久久久久久久大av| 日本撒尿小便嘘嘘汇集6| av.在线天堂| 免费搜索国产男女视频| 久久久欧美国产精品| 香蕉av资源在线| 免费黄网站久久成人精品| 亚洲精品在线观看二区| 给我免费播放毛片高清在线观看| 亚洲国产日韩欧美精品在线观看| 国内揄拍国产精品人妻在线| 国产真实乱freesex| 久久久久免费精品人妻一区二区| 日本精品一区二区三区蜜桃| 亚洲一区高清亚洲精品| 国产极品精品免费视频能看的| 免费大片18禁| 大又大粗又爽又黄少妇毛片口| 最新中文字幕久久久久| 午夜激情欧美在线| 亚洲中文字幕日韩| 日韩欧美免费精品| 国产成人a区在线观看| 久久精品久久久久久噜噜老黄 | 国语自产精品视频在线第100页| 两个人视频免费观看高清| 人人妻人人看人人澡| 国产一区亚洲一区在线观看| 国产av一区在线观看免费| 亚洲乱码一区二区免费版| 亚洲国产精品成人久久小说 | 成人高潮视频无遮挡免费网站| 嫩草影院精品99| 亚洲中文日韩欧美视频| 少妇的逼好多水| 国产综合懂色| 在线观看一区二区三区| 干丝袜人妻中文字幕| 22中文网久久字幕| 最近手机中文字幕大全| 国产精品国产高清国产av| 激情 狠狠 欧美| 人妻丰满熟妇av一区二区三区| 国产蜜桃级精品一区二区三区| 又爽又黄无遮挡网站| 又粗又爽又猛毛片免费看| 给我免费播放毛片高清在线观看| 婷婷亚洲欧美| 成人漫画全彩无遮挡| 国产精品一区二区免费欧美| 少妇丰满av| 最新在线观看一区二区三区| 一进一出抽搐gif免费好疼| 18禁裸乳无遮挡免费网站照片| 久久鲁丝午夜福利片| 久久久精品欧美日韩精品| 亚洲三级黄色毛片| 国产伦精品一区二区三区视频9| 成年免费大片在线观看| 久久久久久九九精品二区国产| 日韩高清综合在线| 国产精品乱码一区二三区的特点| 亚洲内射少妇av| 精品人妻一区二区三区麻豆 | 国产精品爽爽va在线观看网站| 淫妇啪啪啪对白视频| 日韩三级伦理在线观看| 男女那种视频在线观看| 久久久久久伊人网av| 别揉我奶头 嗯啊视频| 日本一二三区视频观看| 国产成人freesex在线 | 欧美日本视频| 最近视频中文字幕2019在线8| 天堂av国产一区二区熟女人妻| 国产精品一区二区三区四区久久| 国产又黄又爽又无遮挡在线| 菩萨蛮人人尽说江南好唐韦庄 | 精品人妻偷拍中文字幕| 别揉我奶头~嗯~啊~动态视频| 精品一区二区三区视频在线观看免费| 久久天躁狠狠躁夜夜2o2o| 国产白丝娇喘喷水9色精品| 91在线精品国自产拍蜜月| 熟妇人妻久久中文字幕3abv| 中文字幕熟女人妻在线| 人人妻人人看人人澡| 久久久精品94久久精品| 成人精品一区二区免费| 国产精品三级大全| 一级毛片aaaaaa免费看小| 天堂√8在线中文| 大又大粗又爽又黄少妇毛片口| 天堂动漫精品| 人妻丰满熟妇av一区二区三区| АⅤ资源中文在线天堂| 国产探花在线观看一区二区| 国产欧美日韩精品亚洲av| 精品久久久久久久人妻蜜臀av| 嫩草影院精品99| 伦精品一区二区三区| 亚洲最大成人av| 成人精品一区二区免费| 嫩草影院精品99| 婷婷亚洲欧美| 嫩草影院新地址| 国产白丝娇喘喷水9色精品| 亚洲精品久久国产高清桃花| 免费av不卡在线播放| 久久6这里有精品| 一个人看视频在线观看www免费| 午夜亚洲福利在线播放| 麻豆一二三区av精品| 九色成人免费人妻av| 欧美最黄视频在线播放免费| 国产男人的电影天堂91| 成人一区二区视频在线观看| 亚洲熟妇熟女久久| 联通29元200g的流量卡| 激情 狠狠 欧美| 午夜日韩欧美国产| 欧美性猛交黑人性爽| 午夜免费男女啪啪视频观看 | 欧美丝袜亚洲另类| 成人毛片a级毛片在线播放| 99在线视频只有这里精品首页| 91久久精品电影网| 亚洲人与动物交配视频| 日韩一本色道免费dvd| 51国产日韩欧美| 又粗又爽又猛毛片免费看| 最近的中文字幕免费完整| 日韩av不卡免费在线播放| 国产成人aa在线观看| 大型黄色视频在线免费观看| 国产成人aa在线观看| 男女啪啪激烈高潮av片| 久久久久性生活片|