亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于關(guān)聯(lián)成像的氣缸蓋內(nèi)壁影像檢測方法研究

        2021-12-10 06:46:28張皓宇楊風暴王肖霞
        內(nèi)燃機工程 2021年6期
        關(guān)鍵詞:熱場單缸氣缸蓋

        張皓宇,楊風暴,王肖霞

        (中北大學 信息與通信工程學院,太原 030051)

        0 概述

        氣缸蓋作為內(nèi)燃機氣缸的核心部件之一,其主要功能是將氣缸密閉形成封閉空間,承受高溫高壓的工作環(huán)境。氣缸蓋上分布有多個氣門、噴油孔、螺栓和其他組件,在氣缸非穩(wěn)定工況下氣缸蓋內(nèi)壁極易產(chǎn)生磨損、凹陷、破裂等損壞,這將導致氣缸各部件相互間的協(xié)調(diào)紊亂,影響內(nèi)燃機氣缸的工作性能[1-3],因此需對氣缸蓋內(nèi)壁的狀況進行檢測。目前常規(guī)的氣缸蓋內(nèi)壁檢測方法需將氣缸拆卸進行水壓、氣壓及目視等專項檢測再重新組裝調(diào)試,存在檢測時間長、耗費人力多、成本消耗較大等缺陷[4-5];而借助計算機模擬軟件對氣缸蓋進行有限元分析,通過溫度的變化規(guī)律及應力應變分析預測氣缸蓋內(nèi)壁容易出現(xiàn)危險點的位置,難以直觀體現(xiàn)影響氣缸正常運行的主要原因,且建立氣缸蓋模型時缺乏切實可靠的邊界條件使得模擬結(jié)果與實際存在差距[6-12]。綜上所述,需要一種新型的更為直觀的方法來實現(xiàn)對氣缸蓋內(nèi)壁的檢測。

        影像技術(shù)在醫(yī)學、夜間探測等領(lǐng)域逐步獲得廣泛的應用[13-14]:將影像技術(shù)應用于醫(yī)學X光診斷中,可實現(xiàn)對人體內(nèi)部器官、組織、腫瘤等的成像,輔助醫(yī)生對患者進行診斷,提高治療方案的準確性;將紅外檢測技術(shù)用于夜間行車探測,通過紅外熱像儀獲取路面目標圖像,識別行人、車輛等目標,可為夜間行車安全提供保障[15]。如果在氣缸蓋內(nèi)壁的檢測中引入影像技術(shù),直觀呈現(xiàn)氣缸蓋內(nèi)壁影像,將利于氣缸蓋內(nèi)壁的有效探測。然而,氣缸蓋內(nèi)壁位于氣缸內(nèi)部,處在高溫高壓的惡劣環(huán)境中,若要獲得相關(guān)圖像信息需將探測器放入氣缸內(nèi)部,而傳統(tǒng)的影像技術(shù)無法對氣缸蓋內(nèi)壁完美呈現(xiàn),因而選擇關(guān)聯(lián)成像。關(guān)聯(lián)成像[16-17]又稱鬼成像(ghost imaging, GI),是一種新型的成像方式。其原理與傳統(tǒng)光學成像不同,目標圖像無需直接獲取,而是由兩部分探測數(shù)據(jù)進行關(guān)聯(lián)計算得到,且探測過程與重構(gòu)影像過程完全獨立。光場關(guān)聯(lián)成像的探測過程收集了不同的光場分布信息及相應作用于目標后的總光強,而重構(gòu)過程則是通過光場分布與總光強值關(guān)聯(lián)計算來還原出目標影像。利用關(guān)聯(lián)成像原理分離熱探測過程與影像重構(gòu)過程,在氣缸蓋上布置傳感器測點,收集內(nèi)壁的熱能量場分布和外壁的總熱量,再通過計算機對數(shù)據(jù)進行關(guān)聯(lián)重構(gòu)影像,這種方式恰好可以解決傳統(tǒng)影像技術(shù)無法利用圖像探測器直接獲得氣缸蓋內(nèi)壁影像的問題。

        本文中針對目前氣缸蓋內(nèi)壁檢測效率低且難以直觀體現(xiàn)的問題,利用關(guān)聯(lián)成像原理來實現(xiàn)氣缸蓋內(nèi)壁影像的有效探測。采用壓縮感知關(guān)聯(lián)成像方法(compressed sensing ghost imaging, CSGI)對氣缸蓋內(nèi)外壁兩部分數(shù)據(jù)進行關(guān)聯(lián)運算,實現(xiàn)氣缸蓋內(nèi)壁影像的間接獲取。通過數(shù)值仿真,利用峰值信噪比及結(jié)構(gòu)相似性等定量評價指標,結(jié)合重構(gòu)影像的直觀性分析驗證了關(guān)聯(lián)成像應用于氣缸蓋內(nèi)壁影像檢測的有效性。該方法可針對發(fā)動機缸內(nèi)高溫高壓環(huán)境中無法對發(fā)動機缸蓋等熱負荷嚴重的關(guān)鍵部件進行直觀分析的問題提供一種新的分析手段,利于對缸蓋內(nèi)壁磨損區(qū)域進行有效判定。

        1 結(jié)合氣缸蓋工作規(guī)律的關(guān)聯(lián)成像分析

        內(nèi)燃機在高溫高壓的工作環(huán)境下,氣缸蓋易在熱應力的作用下產(chǎn)生磨損、開裂等問題。圖1為氣缸蓋底板的熱能傳遞示意圖。在氣缸蓋獲得來自燃燒室的熱能量后,熱能由底板內(nèi)壁火力面①區(qū)域向外壁②區(qū)域傳播,在這一過程中,熱能在氣缸蓋內(nèi)主要以熱傳導的方式傳遞[18-19]。

        圖1 氣缸蓋底板的熱能傳遞

        熱傳導過程遵循導熱基本定律——傅里葉定律,如式(1)所示。

        (1)

        式中,Φ為熱流量;A為傳熱面積;λ為熱導率;T為溫度;x為壁面厚度。

        氣缸蓋內(nèi)壁的熱能量場關(guān)聯(lián)成像[20]與傳統(tǒng)光場的關(guān)聯(lián)成像原理相同,需要通過氣缸蓋內(nèi)壁熱場分布和氣缸蓋底板外壁總能量值兩部分數(shù)據(jù)來實現(xiàn)成像,如圖2所示。其中,氣缸蓋內(nèi)壁熱場分布記錄的是氣缸運行的N個時刻缸蓋內(nèi)壁每一點的熱能值,記作Ia(x,y;N);氣缸蓋底板外壁總能量值則是內(nèi)壁熱場與氣缸蓋相互作用后底板外壁對應每個點透過的熱能值的疊加值,記作Ib(t)。由于內(nèi)壁影像是由氣缸蓋內(nèi)壁每一點的計算結(jié)果疊加而成,那么就需要不斷改變熱場Ia(x,y;N),同時收集Ib(t),以獲得足夠的數(shù)據(jù)量。最后將收集完成的數(shù)據(jù)信息導入計算機,通過關(guān)聯(lián)運算重構(gòu)出氣缸蓋內(nèi)壁影像。

        圖2 氣缸蓋熱場關(guān)聯(lián)成像原理圖

        結(jié)合以上分析,將氣缸蓋底板熱透過系數(shù)函數(shù)用T(x,y)表示,那么熱場與氣缸蓋底板相互作用后的總能量值Ib可表示為式(2)形式。

        (2)

        再通過外壁總能量值Ib與探測的內(nèi)壁熱場分布Ia(x,y) 相互關(guān)聯(lián),重建目標對象TGI(x,y)如式(3)所示。

        (3)

        為獲得高質(zhì)量的氣缸蓋內(nèi)壁影像,整個探測過程需獲得多組不同的氣缸蓋內(nèi)壁熱場分布Ia(x,y)及對應的外壁總能量值Ib,所需的數(shù)據(jù)可使用熱電偶在內(nèi)燃機不同工況下按時刻連續(xù)收集。結(jié)合文獻[21-22],在氣缸加速工作再到減速工作這一過程中,通過氣缸蓋內(nèi)壁不同區(qū)域的溫度場分布計算氣缸蓋內(nèi)壁熱場分布。缸蓋內(nèi)壁測點位置如圖3所示。由于氣缸蓋內(nèi)壁溫度分布具有從中心沿徑向遞減的特點,以四氣門氣缸蓋為例,在各氣門圍繞的中間位置布置測點,其中點1~3為縱向從排氣門側(cè)向進氣門側(cè)分布的測點,點4~8為進排氣門間橫向分布的測點。圖4為各測點的溫度隨時間的變化曲線。

        圖3 缸蓋內(nèi)壁測點位置分布

        圖4 缸蓋不同位置測點的溫度變化曲線

        分析圖4可知,在氣缸由加速工作到減速工作這一過程中,缸蓋內(nèi)壁溫度隨著時間的推進呈現(xiàn)先上升后下降的連續(xù)變化,最高溫度約為400 ℃。但缸蓋內(nèi)壁各測點的變化量是不同的,這樣的變化也體現(xiàn)為相應氣缸蓋內(nèi)壁各點熱能不同,從而使得每次采集的熱場Ia(x,y)之間具有較大的差異性。本研究正是利用這樣的差異性特征為后續(xù)熱場矩陣的數(shù)據(jù)模擬及關(guān)聯(lián)成像技術(shù)的運用提供了基礎條件。

        2 原理和方法

        獲得高質(zhì)量氣缸蓋內(nèi)壁影像的關(guān)鍵在于數(shù)據(jù)的采集及關(guān)聯(lián)方法的選擇。數(shù)據(jù)采集的目的是獲取氣缸蓋內(nèi)壁熱場分布及氣缸蓋外壁總能量值,為關(guān)聯(lián)成像提供足夠的計算數(shù)據(jù)。關(guān)聯(lián)方法會對關(guān)聯(lián)計算采用的數(shù)據(jù)量、關(guān)聯(lián)計算的時間及最終成像質(zhì)量產(chǎn)生影響。

        2.1 關(guān)聯(lián)數(shù)據(jù)的采集

        由于氣缸蓋內(nèi)壁分布的熱場可由氣缸蓋內(nèi)壁的溫度分布情況計算,本研究中利用傳統(tǒng)的熱電偶測量方式對氣缸蓋熱場進行收集。內(nèi)燃機運行由 700 r/min 加速至4 000 r/min,在這一過程中按不同時刻使用陣列式熱電偶對氣缸蓋內(nèi)壁進行一系列的熱場收集;同時在氣缸蓋外壁區(qū)域布置貼片式熱敏電阻,將各個區(qū)域的能量值疊加作為氣缸蓋外壁的總能量值。在使用關(guān)聯(lián)成像對氣缸蓋影像進行計算獲取時,遵循“內(nèi)壁熱場分布與氣缸蓋底板熱透過率函數(shù)的乘積等于氣缸蓋外壁總能量值”的原則,將系列的氣缸蓋內(nèi)壁熱場與外壁總能量值一一對應。

        選擇在不同工況下采集數(shù)據(jù)的主要原因為內(nèi)燃機在不同轉(zhuǎn)速、不同功率及不同負荷的工況下工作時,氣缸蓋內(nèi)壁受到的影響形式不同,會引起內(nèi)壁熱場的分布變化較大,收集的數(shù)據(jù)組實現(xiàn)線性關(guān)聯(lián)運算時容易得到更多像素點位置的解,這將明顯提高重構(gòu)影像的質(zhì)量。

        2.2 壓縮感知關(guān)聯(lián)成像法

        采用壓縮感知關(guān)聯(lián)成像法[23]對探測的數(shù)據(jù)進行關(guān)聯(lián)計算能夠提高氣缸蓋內(nèi)壁重構(gòu)影像的質(zhì)量。作為一種新型的信號處理方法,壓縮感知理論可在獲取可壓縮信號的同時重建稀疏信號,相較于奈奎斯特采樣理論,它通過開發(fā)信號的稀疏特性用隨機采樣獲取信號的離散數(shù)據(jù)樣本,提高數(shù)據(jù)間的不相干性,用較少的采樣次數(shù)通過非線性重建算法還原出氣缸蓋內(nèi)壁重構(gòu)影像。其關(guān)鍵在于信號的稀疏結(jié)構(gòu)和不相關(guān)性。該過程涵蓋稀疏表示、壓縮觀測和信號優(yōu)化重構(gòu)3個部分。

        使用壓縮感知關(guān)聯(lián)成像方法重構(gòu)氣缸蓋內(nèi)壁影像的過程如圖5所示,分為以下幾個步驟:(1) 將氣缸蓋底板熱透過系數(shù)T(x,y)(后文簡寫作T)展開為M×1的列向量,在稀疏變換域下表示為式(4)形式。

        圖5 壓縮感知關(guān)聯(lián)成像方法重構(gòu)氣缸蓋內(nèi)壁影像步驟

        T=ΨS

        (4)

        式中,稀疏矩陣Ψ為M×M維的正交完備基矩陣;S為T在稀疏矩陣Ψ上的系數(shù)矩陣。(2)選取N個適合的氣缸蓋內(nèi)壁熱場強度分布矩陣Ia(x,y)按行展開并組合作為測量矩陣Φ,確保能最大限度保留信號T中的信息,大小為N×M(N?M)。那么,氣缸蓋外壁探測結(jié)果Ib可構(gòu)成一個N×1的數(shù)據(jù)向量Ib,表示為式(5)形式。

        Ib=ΦT=ΦΨS

        (5)

        (3)通過壓縮感知測量矩陣Φ重構(gòu)信號,由Ib來關(guān)聯(lián)計算目標對象TCS,即可轉(zhuǎn)為求解最小l1范數(shù)下的最優(yōu)化,滿足式(6)所示條件。

        TCS=argmin‖ΨS‖l1s.t.Ib=ΦT=ΦΨS

        (6)

        根據(jù)正交匹配追蹤(orthogonal matching pursuit, OMP)算法,重構(gòu)出稀疏系數(shù)S的近似解S′后,再由相應的正交稀疏基Ψ,通過公式(6)即可獲得氣缸蓋內(nèi)壁信息T的近似解TCS,完成目標的重構(gòu)。

        3 仿真試驗分析

        為了驗證氣缸蓋內(nèi)壁關(guān)聯(lián)成像技術(shù)的有效性,本研究中選用了3幅128像素×128像素的氣缸蓋模型二值圖,即單缸雙氣門、單缸四氣門和雙缸雙氣門模型圖作為目標進行試驗。仿真中使用的N個熱場分布矩陣Ia(x,y)是由MATLAB軟件根據(jù)氣缸蓋內(nèi)壁熱場分布規(guī)律區(qū)域性增強后隨機生成,氣缸蓋外壁探測結(jié)果Ib則通過式(5)獲得。對于氣缸蓋內(nèi)壁區(qū)域的影像圖還原,只需使得像素點的尺寸遠小于探測區(qū)域,以保證氣缸蓋內(nèi)壁影像的還原質(zhì)量。這里為了客觀地評估重構(gòu)的氣缸蓋內(nèi)壁影像的質(zhì)量,引入峰值信噪比(peak signal to noise ratio, PSNR)和結(jié)構(gòu)相似性(structural similarity, SSIM)兩個指標進行定量分析,計算方法如式(7)和式(8)所示。

        (7)

        (8)

        (9)

        式中,IMSE為原始圖像與重構(gòu)圖像的均方誤差(mean square error, MSE);f(x,y)與f′(x,y)為相同位置對應的第(i,j)個像素點的原圖和重構(gòu)影像;M=n×m,為總像素點個數(shù);T表示原圖像信息;T′表示重構(gòu)影像信息;μT、μT′為T、T′的平均值;σT、σT′為T、T′的方差;σTT′為T、T′的協(xié)方差;c1、c2為常數(shù),是為了避免分母為0而維持穩(wěn)定。

        圖6~圖8為單缸雙氣門、單缸四氣門和雙缸雙氣門氣缸蓋的二值模型圖及不同采樣次數(shù)下的重構(gòu)影像圖。從圖中可以直觀看出,重構(gòu)出的氣缸蓋內(nèi)壁影像可將氣缸蓋內(nèi)壁上進排氣門口、噴油孔、圓環(huán)凹槽及各受熱區(qū)域完整呈現(xiàn),且隨著關(guān)聯(lián)計算采樣次數(shù)的逐漸增加,氣缸蓋內(nèi)壁影像的清晰度也迅速提高。這是由于在達到全采樣前,每一次采樣數(shù)據(jù)都可能會為氣缸蓋內(nèi)壁影像的重構(gòu)補充有效信息,而隨關(guān)聯(lián)計算結(jié)果的不斷完善即圖片像素點數(shù)據(jù)的不斷補充,重構(gòu)圖像也將越來越接近實際圖像。但關(guān)聯(lián)計算采樣次數(shù)增多的同時也增加了對大量的冗余信息的處理,會使計算時間大大延長。試驗結(jié)果表明,4組采樣次數(shù)下3種氣缸蓋模型圖的重構(gòu)時間分別為143.4 s、263.8 s、379.7 s、509.9 s。根據(jù)以上分析,關(guān)聯(lián)成像的重構(gòu)時間受到原模型圖像素尺寸大小及關(guān)聯(lián)計算采樣次數(shù)的影響而與氣缸蓋類型無關(guān),因此在實際試驗時可選擇合適的采樣率以減少關(guān)聯(lián)計算消耗的時間。

        圖6 單缸雙氣門模型圖及不同采樣次數(shù)下的重構(gòu)影像圖

        圖7 單缸四氣門模型圖及不同采樣次數(shù)下的重構(gòu)影像圖

        圖8 雙缸雙氣門模型圖及不同采樣次數(shù)下的重構(gòu)影像圖

        為了更好地體現(xiàn)方案的有效性,圖9給出了在不同采樣次數(shù)下氣缸蓋內(nèi)壁影像的PSNR和SSIM兩個指標的曲線圖。峰值信噪比PSNR的大小可判斷重構(gòu)圖像的失真程度,結(jié)構(gòu)相似性SSIM指標則定量地從數(shù)值上評判重構(gòu)圖像恢復的質(zhì)量。由圖 9(a) 可知,在采樣次數(shù)達到4 000后,單缸雙氣門與單缸四氣門重構(gòu)影像的PSNR在62 dB左右,雙缸雙氣門的PSNR則達到65 dB以上,且隨著采樣次數(shù)增加逐漸升高,PSNR數(shù)值升高,這代表著重構(gòu)圖像的失真越來越小。再對圖9(b)中SSIM變化曲線進行分析,采樣次數(shù)低于4 000次時,關(guān)聯(lián)計算獲得的有效數(shù)據(jù)較多,結(jié)構(gòu)相似性指標迅速上升,在采樣次數(shù)到達10 000次后,各像素點數(shù)據(jù)逐步完善,冗余信息較多,重構(gòu)影像質(zhì)量提升緩慢并趨于平衡。綜合兩項性能指標曲線圖分析,隨著采樣次數(shù)的增加,重構(gòu)圖像的PSNR和SSIM都隨之提升,且采樣次數(shù)達到8 000次后,PSNR最高可達66 dB且SSIM的值也都在0.4以上,利于氣缸蓋影像的高質(zhì)量重構(gòu)。結(jié)合圖6~圖9進行分析,單缸雙氣門和單缸四氣門氣缸蓋的SSIM性能指標高于雙缸雙氣門,重構(gòu)影像中雙缸雙氣門少部分區(qū)域?qū)毠?jié)的顯現(xiàn)不夠完美。這是由于重構(gòu)影像圖質(zhì)量受到原圖尺寸大小、像素點個數(shù)及有效識別區(qū)域復雜程度的綜合限制,壓縮感知關(guān)聯(lián)成像法未能對每個像素點精確恢復,但總體效果仍比較理想。若使用更高像素的氣缸蓋模型圖并采集足夠多的試驗數(shù)據(jù),那么還原的氣缸蓋內(nèi)壁影像將更清晰,同時選擇更優(yōu)化的關(guān)聯(lián)成像方法也會提升重構(gòu)影像質(zhì)量。

        圖9 不同采樣次數(shù)下氣缸蓋內(nèi)壁重建影像的性能指標曲線

        氣缸蓋常見的損傷主要有變形、開裂及腐蝕,而對氣缸蓋內(nèi)壁影像磨損區(qū)域進行分析時,主要考慮火力面的凹陷、裂紋和凸起損傷。在出現(xiàn)這些問題時,氣缸蓋底板上的熱能傳遞距離、溫度梯度、導熱系數(shù)等參數(shù)變化使得對應位置點的熱透過系數(shù)T(x,y) 發(fā)生改變,這也將直接顯現(xiàn)為磨損區(qū)域的重構(gòu)影像變化。

        為了通過仿真試驗驗證缸蓋磨損區(qū)域的評估的有效性,將各類氣缸蓋模型圖作缺陷處理,再對其進行模擬分析以驗證實際成像效果。圖10、圖11分別為單缸雙氣門和單缸四氣門氣缸蓋處理后的原模型圖及采樣次數(shù)為4 000次、8 000次、12 000次、16 000次的重構(gòu)影像圖,對應的相同采樣次數(shù)下各組的重構(gòu)時間平均為141.2 s、264.5 s、381.3 s、511.8 s。由于重構(gòu)影像圖是通過關(guān)聯(lián)成像法求解計算出各像素點位置的數(shù)據(jù)并疊加還原得到的,那么為使磨損區(qū)域影像能夠清晰顯現(xiàn),則要求構(gòu)成氣缸蓋內(nèi)壁區(qū)域的像素點個數(shù)足夠多,氣缸蓋內(nèi)壁模型圖的像素尺寸就不能太小。根據(jù)圖10、圖11各組仿真試驗結(jié)果進行分析,單缸雙氣門氣缸蓋進排氣門間出現(xiàn)細小裂紋時的情況,4種采樣次數(shù)下的SSIM指標分別為0.366、0.396、0.413、0.424;單缸四氣門氣缸蓋表面出現(xiàn)氣泡等凸起或凹陷情況,4種采樣次數(shù)下的SSIM指標分別為0.375、0.408、0.417、0.436。從重構(gòu)影像圖來看,可以清晰分辨氣缸蓋內(nèi)壁的缺陷區(qū)域的位置和尺寸信息,為氣缸蓋內(nèi)壁具體缺陷類型的判定提供了重要的依據(jù)。

        圖10 單缸雙氣門氣缸蓋磨損情況下的原模型圖及重構(gòu)影像圖

        圖11 單缸四氣門氣缸蓋磨損情況下的原模型圖及重構(gòu)影像圖

        4 結(jié)論

        (1) 基于關(guān)聯(lián)成像的氣缸蓋內(nèi)壁檢測方法能夠在氣缸運行過程中對氣缸蓋進行檢測,能直接呈現(xiàn)氣缸蓋內(nèi)壁影像,直觀性強,利于氣缸蓋火力面的檢測分析。

        (2) 對3種類型氣缸蓋重構(gòu)影像的仿真試驗結(jié)果表明,相應的PSNR和SSIM指標可達66.79 dB、0.427。結(jié)合存在缺陷區(qū)域的氣缸蓋重構(gòu)影像分析,關(guān)聯(lián)成像法仍能清晰地顯現(xiàn)內(nèi)壁影像并能準確辨識缺陷區(qū)域信息。

        猜你喜歡
        熱場單缸氣缸蓋
        動力系統(tǒng)永磁同步電動機流體和熱場的計算分析
        單缸風冷汽油機高壓縮比燃燒室的優(yōu)化設計仿真研究
        熱場中Cu-Nb-C球磨納米粉的尺寸生長與合金化研究
        云南化工(2021年6期)2021-12-21 07:30:54
        柴油發(fā)動機氣缸蓋噴油器孔開裂故障分析及應對
        離心式空壓機國內(nèi)外發(fā)展趨勢簡述
        專利名稱:一種藍寶石晶體生長爐的鎢鉬熱場防短路裝置
        用于高品質(zhì)發(fā)動機的熱場集成模擬方法
        汽車文摘(2017年2期)2017-03-08 07:25:46
        單缸柴油機廢氣渦輪增壓試驗研究
        單缸柴油機齒輪搭配
        人間(2015年23期)2016-01-04 12:47:50
        重型車用氣缸蓋濕型砂鑄造氣孔缺陷的降低
        一区二区三区在线乱码 | 亚洲学生妹高清av| 国产无码十八禁| 牛仔裤人妻痴汉电车中文字幕 | 国产大片黄在线观看| 国产精品成人一区二区三区| 中文字幕第一页亚洲观看| 国产精品中文字幕日韩精品| 国产av无码专区亚洲精品| 男男车车的车车网站w98免费| 久久福利青草精品资源| 亚洲视频在线免费观看一区二区| 日本精品一区二区三区二人码 | 中文字幕爆乳julia女教师| 日本一区二区三区激情视频| 亚洲精品一区二区三区麻豆| 人妻丝袜中文无码av影音先锋专区| 日日碰狠狠躁久久躁9| 国产精品日日摸夜夜添夜夜添| 91精品国产自拍视频| 成人午夜福利视频| 国产剧情国产精品一区| av亚洲在线一区二区| 成人国产精品一区二区八戒网| 午夜无码伦费影视在线观看| 亚洲人成精品久久久久| 中文片内射在线视频播放| 亚洲欧美v国产一区二区| 亚洲精品国精品久久99热一| 亚洲AV无码日韩一区二区乱| 日本黑人乱偷人妻在线播放| 性大毛片视频| 亚洲成av人在线观看无堂无码| 日本人妻高清免费v片| 久久伊人少妇熟女大香线蕉| 欧美人与物videos另类xxxxx| 手机在线免费看av网站| 国模gogo无码人体啪啪| 亚洲熟妇丰满大屁股熟妇| 亚洲国产字幕| 午夜一区二区视频在线观看|