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

    基于稀疏正則化的穩(wěn)態(tài)熱源識(shí)別*

    2020-02-26 08:31:56潘天成呂中榮汪利
    關(guān)鍵詞:熱傳導(dǎo)正則熱源

    潘天成,呂中榮,汪利

    (中山大學(xué)航空航天學(xué)院,廣東 廣州 510006)

    在實(shí)際工程中,存在著許多熱現(xiàn)象,及時(shí)監(jiān)測(cè)結(jié)構(gòu)系統(tǒng)的溫度,識(shí)別熱源的位置和強(qiáng)度,對(duì)整個(gè)熱傳導(dǎo)系統(tǒng)的安全監(jiān)測(cè)和溫度控制具有重要意義。例如,通過測(cè)量核反應(yīng)堆外部的溫度,從而得到其內(nèi)壁溫度,實(shí)時(shí)監(jiān)測(cè)和控制整個(gè)反應(yīng)的溫度;在電路中容易產(chǎn)生局部的熱源,識(shí)別熱源的位置可以用于對(duì)電路進(jìn)行檢修;航天飛機(jī)在進(jìn)入大氣層時(shí),經(jīng)常會(huì)產(chǎn)生氣動(dòng)加熱,航天器的局部結(jié)構(gòu)產(chǎn)生熱源,因此根據(jù)局部溫度場(chǎng)來識(shí)別熱源的位置和強(qiáng)度就顯得尤為重要。熱源識(shí)別屬于熱傳導(dǎo)反問題的研究領(lǐng)域[1],目的在于識(shí)別熱源的位置和強(qiáng)度,及時(shí)掌握系統(tǒng)中熱源的屬性。研究表明,熱源識(shí)別問題一般是非適定的[2],導(dǎo)致識(shí)別結(jié)果對(duì)溫度測(cè)量誤差十分敏感。

    實(shí)際上,熱源識(shí)別的研究主要分為兩類:一是給定熱源的位置或者空間分布,識(shí)別熱源的強(qiáng)度;二是同時(shí)識(shí)別熱源的空間位置和強(qiáng)度。過去幾十年,很多學(xué)者已經(jīng)對(duì)熱源識(shí)別問題做了大量的研究。為求解此類非適定問題,BECK提出在最小二乘方程中添加正則化項(xiàng)求解,發(fā)展并使用了Tikhonov正則化方法。在一維的熱傳導(dǎo)問題中,BLACK WELL提出了序列估計(jì)方法[3]。HUANG等通過共軛梯度法來研究了一維和二維模型中的一個(gè)或兩個(gè)點(diǎn)熱源的強(qiáng)度識(shí)別問題[4]。但是由于算法迭代的特性,計(jì)算時(shí)間比較長(zhǎng)。GENG等采用變分迭代算法求得拋物型熱方程無離散化問題的數(shù)值解并確定熱源的強(qiáng)度[5]。NETO等完成了矩形域內(nèi)特定位置的線熱源隨時(shí)間變化的強(qiáng)度識(shí)別[6]。LI等提出了一種并行蟻群的優(yōu)化算法來識(shí)別二維域中的一個(gè)點(diǎn)熱源[7]。第一類熱源識(shí)別的研究工作比較多,但對(duì)于第二類熱源的識(shí)別研究卻十分有限。主要難點(diǎn)在于測(cè)量得到的溫度與熱源強(qiáng)度和位置之間呈復(fù)雜、難以顯式表達(dá)的非線性關(guān)系,而且在有限元的框架下,測(cè)量數(shù)據(jù)量將小于未知參數(shù)數(shù)目,導(dǎo)致問題的嚴(yán)重非適定性。幸運(yùn)的是,稀疏正則化可以很好地解決這一問題[8-9]。

    本文將提出一種新的方法來研究熱源識(shí)別第二類問題??紤]到點(diǎn)熱源在空間上的分布是稀疏的,引入稀疏正則化[10-12]來對(duì)目標(biāo)函數(shù)進(jìn)行約束, 并且采用交替優(yōu)化方法[13]分別對(duì)溫度和熱源兩個(gè)分離的變量進(jìn)行迭代求解。在尋求正則化參數(shù)的過程中,提出了閾值法。本文所提算法能快速準(zhǔn)確地識(shí)別熱源位置和強(qiáng)度,這給實(shí)際工程的應(yīng)用提供了一定的參考。

    1 稀疏正則化目標(biāo)函數(shù)建立

    1.1 熱傳導(dǎo)問題有限元格式

    考慮一個(gè)包含邊界?Ω的二維域Ω的熱傳導(dǎo)問題,一般的控制方程為

    +Q(x,y,t),

    x,y∈Ω,t≥0

    (1)

    其中,x,y是笛卡爾坐標(biāo),T為溫度,α為不隨時(shí)間變化的熱傳導(dǎo)率,Q是熱源項(xiàng),本文主要考慮為點(diǎn)熱源的情形。

    邊界條件為

    =h(Te(x,y,t)-T(x,y,t))

    在?qΩ邊界上T(x,y,t)=g(x,y,t),在?TΩ邊界上

    ?qΩ∪?TΩ=?Ω,

    ?qΩ∩?TΩ=φ,?TΩ≠φ

    (2)

    其中,nx、ny分別是二維域邊界的法向方向余弦,h是對(duì)流換熱系數(shù),Te(x,y,t)是周圍的環(huán)境溫度,?TΩ和?qΩ分別為本質(zhì)邊界和對(duì)流邊界,g(x,y,t)是在本質(zhì)邊界?TΩ給定的溫度分布。

    初始條件為

    T(x,y,0)=T0(x,y)

    (3)

    其中,T0(x,y)為二維域內(nèi)的初始溫度分布。

    +Q(x,y,t)=0,

    x,y∈Ω,t≥0

    (4)

    有限元方法是求解熱傳導(dǎo)問題(1)-(4)的有效數(shù)值方法,將二維域Ω劃分為有限個(gè)微小的單元dΩ,通過分段線性插值得到的有限元形函數(shù),可以將域內(nèi)的溫度用各節(jié)點(diǎn)的溫度表示。那么,穩(wěn)態(tài)情況下的有限元格式如下:

    (5)

    1.2 建立目標(biāo)函數(shù)

    (6)

    (7)

    1.3 稀疏正則化

    (8)

    其中,λ≥0是一個(gè)正則化參數(shù)。至此,基于稀疏正則化的有限元格式穩(wěn)態(tài)熱源識(shí)別目標(biāo)函數(shù)已經(jīng)建立。

    2 交替優(yōu)化方法

    2.1 交替優(yōu)化方法求解目標(biāo)函數(shù)

    2.1.1 溫度恢復(fù)步 完整的溫度數(shù)據(jù)從測(cè)量得到的溫度中得到恢復(fù)。

    (9)

    +(F(k-1))TF(k-1)

    +(F(k-1))TF(k-1)

    其中,公式(10)中的上標(biāo)S和U分別對(duì)應(yīng)測(cè)量的節(jié)點(diǎn)和未測(cè)量的節(jié)點(diǎn),上標(biāo)US對(duì)應(yīng)耦合的部分,I是一個(gè)單位矩陣。因此可以得到公式(9)的解,如下

    2.1.2 熱源識(shí)別步

    從恢復(fù)得到的完整溫度數(shù)據(jù)中識(shí)別熱源位置和強(qiáng)度,如下

    (12)

    (13)

    其中,(ν)j表示向量ν的第j個(gè)分量。顯然,問題(13)是具有顯式解。因此,熱源識(shí)別步,即公式(12)的解為

    (14)

    可以發(fā)現(xiàn),公式(14)容易得到稀疏正則化目標(biāo)函數(shù)的最優(yōu)解,而且目標(biāo)函數(shù)是解耦的,無需迭代。進(jìn)一步的研究表明,稀疏正則化參數(shù)λ可能會(huì)對(duì)熱源的強(qiáng)度識(shí)別造成一定的偏差,因此應(yīng)避免,本文對(duì)公式(14)進(jìn)行了修正,消除λ對(duì)識(shí)別的熱源強(qiáng)度的影響。

    =fm(λ;aj)

    (15)

    關(guān)于λ的選取將在接下來進(jìn)行討論。

    2.2 閾值法尋求正則化參數(shù)

    在執(zhí)行每一熱源識(shí)別步之前,都必須選取一個(gè)合適的正則化參數(shù)。根據(jù)公式(15),可選取一個(gè)閾值集τ。

    τ={λjcr:=2|aj|,j=1,2,…n}

    (16)

    當(dāng)λ≥λjcr時(shí),公式(15)的解為0,即表示在第j個(gè)節(jié)點(diǎn)內(nèi)部沒有點(diǎn)熱源存在;λ<λjcr時(shí),結(jié)果表明點(diǎn)熱源在相應(yīng)的節(jié)點(diǎn)內(nèi)部。將閾值集進(jìn)行降序排序

    (17)

    那么就可以通過閾值法來選取一個(gè)正則化參數(shù),基本算法如圖1所示。其中,lmax和γ是兩個(gè)閾值參數(shù),lmax表示點(diǎn)熱源存在的最大個(gè)數(shù);γ是一個(gè)判別比,通常在lmax比較小的情況下,取值較大。

    圖1 閾值法確定正則化參數(shù)流程圖Fig.1 The flow chart for determining regularization parameters by threshold setting method

    3 數(shù)值算例

    3.1 模型建立

    本文考慮一個(gè)二維穩(wěn)態(tài)的數(shù)值算例。通常,測(cè)量的溫度是從有限元模擬并添加一定的噪聲得到。在穩(wěn)態(tài)情況下,噪聲施加下

    (18)

    其中,Tj是從公式(5)求得的在第j個(gè)節(jié)點(diǎn)的溫度數(shù)據(jù)。ε是噪聲水平,一般為0%、2%、5%等等。r是一個(gè)標(biāo)準(zhǔn)正態(tài)分布的隨機(jī)數(shù)。

    如圖1所示,薄板為鋁合金材料,長(zhǎng)和寬均為0.5 m,密度ρ=2 787 kg/m3,比熱容c=883 J/(kg·K),熱傳導(dǎo)率α=194 W/(m·℃), 對(duì)流換熱系數(shù)h=10 W/(m2·K),環(huán)境溫度為0 ℃,一個(gè)或多個(gè)熱源在薄板的內(nèi)部。每個(gè)熱源的面積為0.004 m×0.004 m,相對(duì)于板的面積是非常小的,因此可以看作是點(diǎn)熱源。

    薄板的有限元單元和節(jié)點(diǎn)信息,如圖2所示。

    圖2 二維穩(wěn)態(tài)模型圖Fig.2 Two-dimensional steady-state model

    考慮以下幾種工況,來驗(yàn)證所提方法的有效性和準(zhǔn)確性,熱源的強(qiáng)度Q15=1.2×106W/m2,Q20=1.0×106W/m2。

    3.2 熱源識(shí)別結(jié)果

    在識(shí)別表1四種工況的過程中,考慮噪聲的隨機(jī)性,進(jìn)行蒙特卡洛(monte carlo)試驗(yàn),每種工況取100組數(shù)據(jù)進(jìn)行研究,最后取其均值和標(biāo)準(zhǔn)差作為識(shí)別結(jié)果。在算法實(shí)現(xiàn)過程中,閾值法選取的參數(shù)分別為lmax=4,γ=5,β=1×106,fb的值則在每一次迭代的過程中選取矩陣ASS的2范數(shù),即fb=norm(ASS)。

    四種工況稀疏正則化與無正則化的識(shí)別結(jié)果如圖3-6所示,每種工況第20節(jié)點(diǎn)內(nèi)部的熱源強(qiáng)度識(shí)別結(jié)果及相對(duì)誤差如表2所示。在無噪聲的情況下,目標(biāo)函數(shù)中有稀疏正則化與無正則化均能準(zhǔn)確識(shí)別熱源的位置和強(qiáng)度,無正則化識(shí)別的強(qiáng)度誤差也為0。當(dāng)有噪聲水平的時(shí)候,稀疏正則化的識(shí)別結(jié)果明顯比無正則化的結(jié)果要好,而且能保證解的稀疏性,而無正則化的時(shí)候不能保證。同時(shí)發(fā)現(xiàn),噪聲水平的提高,蒙特卡洛實(shí)驗(yàn)得到的熱源識(shí)別標(biāo)準(zhǔn)差也會(huì)增大,但不影響熱源的定位。工況4表明,僅通過邊界點(diǎn)的測(cè)量,本文所提方法也能準(zhǔn)確定位結(jié)構(gòu)內(nèi)部的熱源,且識(shí)別的強(qiáng)度誤差不大,誤差為3.62%,而無正則化的算法甚至無法準(zhǔn)確定位。在測(cè)量數(shù)據(jù)量比較少的情況下,稀疏正則化的優(yōu)勢(shì)明顯地體現(xiàn)出來了。

    圖3 工況1有無稀疏正則化識(shí)別結(jié)果對(duì)比Fig.3 Identified results without and with regularization for scenario 1

    表2 第20節(jié)點(diǎn)熱源強(qiáng)度識(shí)別結(jié)果Table 2 The identified results of heat source strength at 20th node

    圖4 工況2有無稀疏正則化識(shí)別結(jié)果對(duì)比Fig.4 Identified results without and with regularization for scenario 2

    圖5 工況3有無稀疏正則化識(shí)別結(jié)果對(duì)比Fig.5 Identified results without and with regularization for scenario 3

    圖6 工況4有無稀疏正則化識(shí)別結(jié)果對(duì)比Fig.6 Identified results without and with regularization for scenario 4

    4 結(jié) 論

    本文考慮到點(diǎn)熱源在空間分布上的稀疏性,將稀疏正則化與交替優(yōu)化方法結(jié)合,提出了一種新的熱源識(shí)別方法,很好地克服了熱源識(shí)別非適定性的問題。數(shù)值算例的結(jié)果表明:該方法能同時(shí)識(shí)別熱源的位置和強(qiáng)度,稀疏正則化能加強(qiáng)解的稀疏性,具有較好的抗噪性,并且交替優(yōu)化方法的求解和閾值法尋求正則化參數(shù)使得識(shí)別快速準(zhǔn)確,甚至僅通過邊界的測(cè)量也可識(shí)別內(nèi)部的熱源,這給實(shí)際工程提供了一定的參考意義。

    猜你喜歡
    熱傳導(dǎo)正則熱源
    一類三維逆時(shí)熱傳導(dǎo)問題的數(shù)值求解
    橫流熱源塔換熱性能研究
    煤氣與熱力(2021年3期)2021-06-09 06:16:20
    剩余有限Minimax可解群的4階正則自同構(gòu)
    類似于VNL環(huán)的環(huán)
    熱傳導(dǎo)方程解的部分Schauder估計(jì)
    一類非線性反向熱傳導(dǎo)問題的Fourier正則化方法
    基于啟發(fā)式動(dòng)態(tài)規(guī)劃的冷熱源優(yōu)化控制
    電子制作(2017年19期)2017-02-02 07:08:31
    中部槽激光-MAG復(fù)合熱源打底焊焊接工藝研究
    焊接(2015年8期)2015-07-18 10:59:13
    多類型熱源點(diǎn)共存下的區(qū)域熱力統(tǒng)籌供應(yīng)探討
    有限秩的可解群的正則自同構(gòu)
    在线免费观看不下载黄p国产| 久久精品久久久久久久性| 亚洲精品中文字幕在线视频| 久久精品国产a三级三级三级| 日韩电影二区| 麻豆成人av视频| 久久女婷五月综合色啪小说| 九九久久精品国产亚洲av麻豆| 国产精品免费大片| 99久久综合免费| 欧美成人精品欧美一级黄| 免费黄网站久久成人精品| 国产精品麻豆人妻色哟哟久久| a 毛片基地| 久久久久久久亚洲中文字幕| 黄色配什么色好看| 国产黄色免费在线视频| 国产高清三级在线| 亚洲熟女精品中文字幕| 久久影院123| 亚洲婷婷狠狠爱综合网| 午夜激情福利司机影院| a级片在线免费高清观看视频| 欧美日韩综合久久久久久| 成人综合一区亚洲| 午夜免费观看性视频| 亚洲国产av新网站| 永久网站在线| 国产欧美亚洲国产| 熟女av电影| 亚洲精品一二三| 久久这里有精品视频免费| 亚洲三级黄色毛片| av在线老鸭窝| 97在线人人人人妻| 日韩电影二区| 国产色婷婷99| 日韩电影二区| av专区在线播放| 国产精品麻豆人妻色哟哟久久| 久久久久久久久久成人| 国产一区二区三区av在线| 国产欧美另类精品又又久久亚洲欧美| 色视频在线一区二区三区| 老司机亚洲免费影院| 亚洲不卡免费看| 婷婷色av中文字幕| 在线观看一区二区三区激情| 看十八女毛片水多多多| 全区人妻精品视频| 精品久久久久久久久av| 2021少妇久久久久久久久久久| 丝瓜视频免费看黄片| 亚洲精华国产精华液的使用体验| 久久热精品热| 最黄视频免费看| 校园人妻丝袜中文字幕| 丰满饥渴人妻一区二区三| 国产乱人偷精品视频| 人人妻人人澡人人爽人人夜夜| 久久免费观看电影| 精品午夜福利在线看| 久久青草综合色| 久久久久精品性色| 建设人人有责人人尽责人人享有的| 欧美变态另类bdsm刘玥| 色网站视频免费| 只有这里有精品99| 中国国产av一级| 大香蕉97超碰在线| 中文字幕免费在线视频6| 国产日韩欧美在线精品| 久久女婷五月综合色啪小说| 一区二区三区精品91| 中文字幕最新亚洲高清| 黄色视频在线播放观看不卡| 亚洲av成人精品一区久久| 欧美xxⅹ黑人| 日韩精品有码人妻一区| 久久99热这里只频精品6学生| 热re99久久精品国产66热6| 中文字幕久久专区| 夜夜看夜夜爽夜夜摸| 一区二区三区精品91| 国产精品人妻久久久影院| 天堂8中文在线网| 最近2019中文字幕mv第一页| 欧美97在线视频| 十八禁网站网址无遮挡| 国产精品国产三级国产av玫瑰| 国产成人av激情在线播放 | 天堂8中文在线网| 亚洲精品456在线播放app| av专区在线播放| 国产有黄有色有爽视频| 国产片内射在线| 2021少妇久久久久久久久久久| 欧美3d第一页| 夫妻性生交免费视频一级片| 日韩不卡一区二区三区视频在线| av不卡在线播放| 18禁观看日本| 日本黄色日本黄色录像| 日韩大片免费观看网站| 国产在视频线精品| 纯流量卡能插随身wifi吗| 看十八女毛片水多多多| 搡女人真爽免费视频火全软件| 免费观看无遮挡的男女| 三上悠亚av全集在线观看| 18+在线观看网站| 国产精品 国内视频| 亚洲精品自拍成人| 自线自在国产av| 免费日韩欧美在线观看| 熟女人妻精品中文字幕| 久久人人爽人人片av| 日本av手机在线免费观看| 亚洲成人手机| 青春草国产在线视频| 在线亚洲精品国产二区图片欧美 | 免费大片黄手机在线观看| 国产精品.久久久| 精品人妻偷拍中文字幕| 欧美亚洲 丝袜 人妻 在线| 免费播放大片免费观看视频在线观看| 少妇人妻 视频| 亚洲av男天堂| 国产一区亚洲一区在线观看| 久久久精品94久久精品| 一区二区三区乱码不卡18| 国产精品熟女久久久久浪| 亚洲欧美中文字幕日韩二区| 国产精品嫩草影院av在线观看| av不卡在线播放| 亚洲欧美中文字幕日韩二区| 久久精品久久久久久久性| 老熟女久久久| tube8黄色片| 18+在线观看网站| 80岁老熟妇乱子伦牲交| 久久精品人人爽人人爽视色| 夫妻性生交免费视频一级片| 91午夜精品亚洲一区二区三区| 欧美精品一区二区免费开放| 欧美日韩视频精品一区| 免费观看a级毛片全部| 少妇人妻 视频| 精品一区二区免费观看| 综合色丁香网| 嫩草影院入口| 熟女人妻精品中文字幕| 美女国产视频在线观看| 日本猛色少妇xxxxx猛交久久| 中文字幕最新亚洲高清| 最后的刺客免费高清国语| 久久久久久久国产电影| 免费黄频网站在线观看国产| 久久精品国产亚洲av涩爱| 日本黄色日本黄色录像| 免费人妻精品一区二区三区视频| 国产精品女同一区二区软件| 国产免费又黄又爽又色| 伊人亚洲综合成人网| 日韩中文字幕视频在线看片| 免费看av在线观看网站| 五月开心婷婷网| 在线天堂最新版资源| 亚洲国产欧美在线一区| 成年人午夜在线观看视频| 内地一区二区视频在线| 日韩一本色道免费dvd| 国产精品免费大片| 麻豆精品久久久久久蜜桃| 伦理电影大哥的女人| 国产伦理片在线播放av一区| 久久精品久久久久久噜噜老黄| 国产高清有码在线观看视频| 在线亚洲精品国产二区图片欧美 | 精品久久久久久电影网| 一本大道久久a久久精品| 蜜桃在线观看..| 亚洲精品亚洲一区二区| 国产精品一国产av| 高清欧美精品videossex| tube8黄色片| 99国产综合亚洲精品| 男女边吃奶边做爰视频| 少妇 在线观看| 青春草视频在线免费观看| 在线观看免费日韩欧美大片 | 91精品国产九色| 久久精品人人爽人人爽视色| 中国三级夫妇交换| 欧美日韩国产mv在线观看视频| 国产av一区二区精品久久| 天天躁夜夜躁狠狠久久av| 久久久精品94久久精品| 久久久a久久爽久久v久久| 久久国产精品大桥未久av| 国产一区二区三区av在线| 午夜老司机福利剧场| av线在线观看网站| 精品国产一区二区三区久久久樱花| 激情五月婷婷亚洲| 超碰97精品在线观看| 国产成人一区二区在线| 我的女老师完整版在线观看| 一二三四中文在线观看免费高清| av卡一久久| 久久精品久久精品一区二区三区| 男人操女人黄网站| 又粗又硬又长又爽又黄的视频| 啦啦啦视频在线资源免费观看| 中文字幕av电影在线播放| av女优亚洲男人天堂| 在线观看www视频免费| 在线精品无人区一区二区三| 你懂的网址亚洲精品在线观看| 两个人的视频大全免费| 哪个播放器可以免费观看大片| 人人妻人人澡人人看| 国产成人精品久久久久久| 国产精品一区二区三区四区免费观看| 中文字幕人妻丝袜制服| 免费高清在线观看视频在线观看| 少妇被粗大猛烈的视频| 亚洲精品色激情综合| 欧美日韩视频高清一区二区三区二| 99久国产av精品国产电影| 中文字幕制服av| 午夜免费观看性视频| 欧美最新免费一区二区三区| 人成视频在线观看免费观看| 亚洲国产欧美在线一区| 久久97久久精品| 狂野欧美白嫩少妇大欣赏| 中文字幕人妻熟人妻熟丝袜美| 久久亚洲国产成人精品v| 欧美xxⅹ黑人| 欧美另类一区| av在线老鸭窝| 国产 一区精品| 老司机亚洲免费影院| 两个人免费观看高清视频| 国产黄色免费在线视频| 国产极品粉嫩免费观看在线 | 99热全是精品| 欧美最新免费一区二区三区| 如日韩欧美国产精品一区二区三区 | av.在线天堂| 久久人人爽人人爽人人片va| 亚洲精品成人av观看孕妇| 插逼视频在线观看| 美女xxoo啪啪120秒动态图| 建设人人有责人人尽责人人享有的| 九九在线视频观看精品| 久久午夜综合久久蜜桃| 久久久久久久国产电影| 大香蕉久久成人网| 我要看黄色一级片免费的| 国产精品99久久久久久久久| 中国三级夫妇交换| 熟女人妻精品中文字幕| 高清午夜精品一区二区三区| 中国国产av一级| 婷婷成人精品国产| 22中文网久久字幕| 亚洲情色 制服丝袜| 一边摸一边做爽爽视频免费| 国产女主播在线喷水免费视频网站| xxx大片免费视频| 成人国语在线视频| 九九爱精品视频在线观看| 婷婷成人精品国产| 国产色爽女视频免费观看| 国产一区二区三区综合在线观看 | 久久精品久久久久久久性| 久久国产精品男人的天堂亚洲 | 黄色欧美视频在线观看| 两个人的视频大全免费| 日韩制服骚丝袜av| 一区二区三区四区激情视频| 一边摸一边做爽爽视频免费| 久久这里有精品视频免费| 全区人妻精品视频| 久久久久久伊人网av| 亚洲综合精品二区| a级片在线免费高清观看视频| 少妇被粗大猛烈的视频| 性色av一级| 一级毛片黄色毛片免费观看视频| 9色porny在线观看| 国产成人午夜福利电影在线观看| 男人操女人黄网站| 久久久精品免费免费高清| 母亲3免费完整高清在线观看 | 免费日韩欧美在线观看| 熟妇人妻不卡中文字幕| av免费观看日本| 能在线免费看毛片的网站| 不卡视频在线观看欧美| 亚洲av福利一区| 亚洲av综合色区一区| 22中文网久久字幕| 少妇猛男粗大的猛烈进出视频| 男女无遮挡免费网站观看| 亚洲熟女精品中文字幕| av女优亚洲男人天堂| 欧美成人午夜免费资源| 22中文网久久字幕| 一区二区三区免费毛片| 国产成人freesex在线| 交换朋友夫妻互换小说| 午夜免费鲁丝| 久久精品熟女亚洲av麻豆精品| 亚洲av中文av极速乱| 欧美日韩国产mv在线观看视频| 赤兔流量卡办理| 日韩中文字幕视频在线看片| 在线亚洲精品国产二区图片欧美 | 精品人妻熟女av久视频| 亚洲av二区三区四区| 久久精品久久精品一区二区三区| 精品人妻在线不人妻| 亚洲熟女精品中文字幕| 国产视频内射| 欧美三级亚洲精品| 七月丁香在线播放| 丝瓜视频免费看黄片| 久久久亚洲精品成人影院| 国产成人免费观看mmmm| 国产无遮挡羞羞视频在线观看| 少妇高潮的动态图| 国产成人av激情在线播放 | 午夜免费男女啪啪视频观看| 自拍欧美九色日韩亚洲蝌蚪91| 99热这里只有是精品在线观看| 国产成人精品福利久久| 啦啦啦中文免费视频观看日本| 精品一品国产午夜福利视频| 亚洲美女视频黄频| 中文字幕免费在线视频6| 我的老师免费观看完整版| 丝袜美足系列| 国产毛片在线视频| 久久人人爽人人爽人人片va| 黄色欧美视频在线观看| 日韩伦理黄色片| 黄片播放在线免费| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲精品日本国产第一区| 欧美+日韩+精品| 黑人欧美特级aaaaaa片| 亚洲av成人精品一区久久| 亚洲精品久久久久久婷婷小说| 免费黄色在线免费观看| 高清午夜精品一区二区三区| 久久久久人妻精品一区果冻| 91成人精品电影| 大香蕉久久成人网| 欧美亚洲日本最大视频资源| 一级,二级,三级黄色视频| 久久久久久久精品精品| 一个人看视频在线观看www免费| 免费不卡的大黄色大毛片视频在线观看| 另类精品久久| 久久免费观看电影| 国内精品宾馆在线| 高清黄色对白视频在线免费看| 亚洲色图综合在线观看| 国产精品人妻久久久久久| 一级片'在线观看视频| 亚洲图色成人| 男女边摸边吃奶| 丝袜美足系列| 亚洲精品成人av观看孕妇| 精品少妇久久久久久888优播| 一本大道久久a久久精品| 特大巨黑吊av在线直播| 少妇被粗大的猛进出69影院 | 精品一区二区三区视频在线| 亚洲丝袜综合中文字幕| 美女大奶头黄色视频| 色哟哟·www| 丝瓜视频免费看黄片| 欧美日韩综合久久久久久| 一本—道久久a久久精品蜜桃钙片| 中国国产av一级| 青春草国产在线视频| 国产精品不卡视频一区二区| 人妻一区二区av| 最新中文字幕久久久久| 性色av一级| 中文天堂在线官网| 最近手机中文字幕大全| av在线app专区| 久久这里有精品视频免费| 久久青草综合色| 嫩草影院入口| 国产爽快片一区二区三区| 不卡视频在线观看欧美| 中文字幕人妻丝袜制服| 国产精品麻豆人妻色哟哟久久| 看免费成人av毛片| 国产午夜精品久久久久久一区二区三区| 国产男人的电影天堂91| 少妇人妻精品综合一区二区| 精品人妻熟女毛片av久久网站| 国产精品一国产av| 亚洲成人手机| 日韩成人伦理影院| 人人妻人人爽人人添夜夜欢视频| av在线观看视频网站免费| 久久午夜综合久久蜜桃| 韩国高清视频一区二区三区| 久久毛片免费看一区二区三区| 欧美精品一区二区免费开放| 99九九线精品视频在线观看视频| 我要看黄色一级片免费的| 久久久久久久久久久丰满| 成人亚洲欧美一区二区av| 另类亚洲欧美激情| 日本-黄色视频高清免费观看| 午夜久久久在线观看| 亚洲五月色婷婷综合| 亚洲人成网站在线观看播放| 99热国产这里只有精品6| 插阴视频在线观看视频| 国产色婷婷99| 一区二区av电影网| 亚洲四区av| 国产精品国产av在线观看| 亚洲精品视频女| 黄色怎么调成土黄色| 久久鲁丝午夜福利片| 成人18禁高潮啪啪吃奶动态图 | 51国产日韩欧美| 亚洲精品,欧美精品| 国产成人aa在线观看| 制服丝袜香蕉在线| 最近的中文字幕免费完整| 亚洲国产精品一区二区三区在线| 最近的中文字幕免费完整| 久久免费观看电影| 天堂8中文在线网| 岛国毛片在线播放| 视频中文字幕在线观看| 人妻 亚洲 视频| 成人影院久久| 国产毛片在线视频| 欧美性感艳星| 一个人看视频在线观看www免费| 丰满迷人的少妇在线观看| 久久这里有精品视频免费| 国产精品女同一区二区软件| 亚洲av二区三区四区| 如日韩欧美国产精品一区二区三区 | av又黄又爽大尺度在线免费看| 亚洲五月色婷婷综合| 三级国产精品欧美在线观看| 亚洲高清免费不卡视频| 国产成人精品无人区| 亚洲欧美清纯卡通| 两个人免费观看高清视频| 高清av免费在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 婷婷色综合www| 永久免费av网站大全| 插阴视频在线观看视频| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久精品电影小说| 黑丝袜美女国产一区| 国产精品嫩草影院av在线观看| 免费观看av网站的网址| 韩国av在线不卡| freevideosex欧美| 如何舔出高潮| 欧美bdsm另类| 99久久人妻综合| 这个男人来自地球电影免费观看 | av在线老鸭窝| 欧美成人午夜免费资源| 国产淫语在线视频| 满18在线观看网站| 男男h啪啪无遮挡| 熟女av电影| 亚洲av福利一区| 特大巨黑吊av在线直播| 日产精品乱码卡一卡2卡三| av国产久精品久网站免费入址| 熟女av电影| 成人国语在线视频| 婷婷成人精品国产| 伊人亚洲综合成人网| 精品午夜福利在线看| 国产欧美日韩一区二区三区在线 | 亚州av有码| 亚洲精品第二区| 熟女av电影| 在线观看免费日韩欧美大片 | 99久久精品国产国产毛片| 一个人免费看片子| 亚洲av中文av极速乱| 久久午夜福利片| 蜜桃在线观看..| 少妇被粗大猛烈的视频| 最近手机中文字幕大全| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 老司机影院毛片| 久久久精品94久久精品| 亚洲欧美日韩另类电影网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 男女边摸边吃奶| 啦啦啦在线观看免费高清www| 久久精品国产亚洲网站| 80岁老熟妇乱子伦牲交| 国产精品不卡视频一区二区| 亚洲国产精品专区欧美| 狂野欧美激情性bbbbbb| 欧美+日韩+精品| av免费在线看不卡| 婷婷色av中文字幕| 国产精品久久久久久久电影| 久久影院123| 在线观看人妻少妇| 久久99蜜桃精品久久| 我要看黄色一级片免费的| 亚洲,欧美,日韩| 大片电影免费在线观看免费| 欧美 日韩 精品 国产| 久久99一区二区三区| 一区二区三区四区激情视频| 亚洲丝袜综合中文字幕| 哪个播放器可以免费观看大片| 成人无遮挡网站| 丁香六月天网| 精品熟女少妇av免费看| 国产免费视频播放在线视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产伦理片在线播放av一区| 中文字幕最新亚洲高清| 国产av精品麻豆| 亚洲欧美清纯卡通| 免费观看在线日韩| 香蕉精品网在线| 国产一区二区在线观看日韩| 一边摸一边做爽爽视频免费| 日本与韩国留学比较| 亚洲四区av| 美女主播在线视频| 亚洲美女搞黄在线观看| 一区在线观看完整版| 国产日韩欧美在线精品| 欧美国产精品一级二级三级| 久久久国产精品麻豆| 2021少妇久久久久久久久久久| 久久国产精品大桥未久av| 久久人人爽人人片av| 成人黄色视频免费在线看| 纯流量卡能插随身wifi吗| 欧美精品国产亚洲| 免费看不卡的av| 精品少妇久久久久久888优播| 欧美激情 高清一区二区三区| 午夜福利影视在线免费观看| 91久久精品电影网| 午夜激情福利司机影院| 国产女主播在线喷水免费视频网站| 国产黄色视频一区二区在线观看| freevideosex欧美| 人妻系列 视频| 久久久久视频综合| av视频免费观看在线观看| 精品国产国语对白av| 欧美xxⅹ黑人| 久久久久国产精品人妻一区二区| www.色视频.com| 亚洲第一av免费看| 亚洲精品国产色婷婷电影| 久久久久久久久久人人人人人人| 在线 av 中文字幕| 午夜日本视频在线| 欧美丝袜亚洲另类| 日韩中字成人| 日韩一区二区视频免费看| 制服人妻中文乱码| 岛国毛片在线播放| 久久韩国三级中文字幕| 亚洲av中文av极速乱| 亚洲国产欧美日韩在线播放| 久久久久精品久久久久真实原创| 国产午夜精品一二区理论片| 久久久久久久精品精品| 岛国毛片在线播放| 少妇熟女欧美另类| 精品国产乱码久久久久久小说| 妹子高潮喷水视频| 2018国产大陆天天弄谢| 亚洲性久久影院| 91在线精品国自产拍蜜月| 男的添女的下面高潮视频| 在线观看www视频免费| 亚洲天堂av无毛| 99久久精品一区二区三区| 丰满迷人的少妇在线观看| 成人免费观看视频高清| 91精品伊人久久大香线蕉| 91国产中文字幕| 亚洲精品国产av蜜桃| 久久久久久久久大av| 日韩,欧美,国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 日本免费在线观看一区| 亚洲av不卡在线观看|