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

    葉片軸截面曲線擬合方法

    2019-04-10 07:01:04郭子健
    科技資訊 2019年35期
    關鍵詞:曲線擬合渦輪發(fā)動機

    郭子健

    摘 ?要:渦輪葉片作為飛機發(fā)動機的關鍵零部件,渦輪葉片和發(fā)動機高溫區(qū)的零部件壽命往往是制約現(xiàn)代高性能航空器壽命周轉的決定性因素。航空發(fā)動機零部件修復可謂是高投入、高附加值、高層次政策下的高技術產(chǎn)業(yè),是一項軍民通用技術。其巨大的經(jīng)濟利益吸引了許多國家和地區(qū)紛紛涉足此領域。逆向工程是一種高效率、低成本的新型設計方式。該文基于逆向工程技術,通過對點云數(shù)據(jù)進行測量,并對點云數(shù)據(jù)進行預處理,進而將曲面重構等技術對渦輪葉片的軸截面特性進行分析以及葉型進行擬合,為飛機發(fā)動機葉片的逆向重構提供了前提。這為航空發(fā)動機渦輪葉片損耗的有效降低和高效維修奠定了基礎,具有一定的理論價值和實際意義。

    關鍵詞:葉片軸截面特性 ?葉片修復 ?逆

    中圖分類號:TP391 ? ? ? ? ? ? ? ?文獻標識碼:A 文章編號:1672-3791(2019)12(b)-0001-05

    1 ?緒論

    1.1 葉片修復的發(fā)展現(xiàn)狀和經(jīng)濟意義以及該課題的意義和目的

    1.1.1 航空發(fā)動機葉片修復的發(fā)展現(xiàn)狀和經(jīng)濟意義

    航空發(fā)動機葉片,特別是渦輪葉片的工作環(huán)境非常惡劣,要承受高溫、高壓燃氣的嚴重腐蝕,還要承受發(fā)動機高速轉動所帶來的葉片自身的離心力、氣動力、熱應力以及振動負荷;當發(fā)動機工況不斷變化時,葉片還要經(jīng)受冷熱疲勞。因此,渦輪葉片是發(fā)動機各部件中受力和受熱最為嚴重的零件之一,也經(jīng)常出現(xiàn)一些失效的現(xiàn)象。失效原因一般包括斷裂、疲勞裂紋、伸長或縮短量超標、外物打傷程度超過了修理允許量、葉冠變形、使用過熱或超溫、腐蝕、涂層破壞等。由于渦輪葉片是航空發(fā)動機內(nèi)主要的損傷葉片,也是價格最為昂貴的部件。而且通過國內(nèi)各大航空公司(中國國際航空公司、中國南方航空公司等)以及發(fā)動機維修廠(AMECO飛機維修公司、廈門太古飛機維修公司等)的調(diào)研而知,國內(nèi)對渦輪葉片的維修技術都只停留在非常低級的水平,甚至有些連渦輪葉片的清洗工作都沒有授權或者做不了。

    最先于20世紀70年代中期,美國惠普公司就開始研究發(fā)動機渦輪葉片的修理技術,并取得了顯著的成效,推向了市場。之后,英國、法國等歐洲國家也共同制定了COST50計劃,進行了大量的基礎性研究,并對羅羅公司生產(chǎn)的RB211發(fā)動機和其他渦輪部件的修復技術進行研究。這些年來,隨著材料科學的研究、質(zhì)控技術以及零件生產(chǎn)制造工藝的迅速發(fā)展,發(fā)動機渦輪部件維修也得到了很大的發(fā)展。雖然渦輪葉片的修復價格高昂,但是她可以顯著提高葉片的使用壽命,從而使葉片的每小時使用費顯著降低。比如,預防性維修可以在渦輪葉片服役3500h之后進行,則維修費用為葉片成本的5%,葉片的維修后還可以服役3500h,則分攤到葉片每個工作小時的成本僅為新葉片價格的0.015%;如果對服役3500h后的葉片進行大修,并對其進行熱等靜壓處理,此修理成本為新葉片價格的20%。一般認為若葉片修復的成本小于新葉片價格的70%則修復工作是經(jīng)濟可行的[1]。

    1.1.2 該課題的意義和目的

    該課題基于逆向工程技術,通過對點云數(shù)據(jù)進行測量,并對點云數(shù)據(jù)進行預處理,進而將曲面重構等技術對渦輪葉片的軸截面特性進行分析以及葉型進行擬合,以實現(xiàn)飛機發(fā)動機葉片的逆向重構,為飛機發(fā)動機葉片的高效修復奠定了基礎,具有一定的經(jīng)濟意義和實際價值。

    1.2 論文主要研究的內(nèi)容

    該文所需討論研究的內(nèi)容主要是對葉片軸截面的曲線進行擬合和分析。

    在葉片設計時,雖然有很多可以解決強度、剛度和振動問題的先進設計方法,但是葉片的燒蝕和損傷在發(fā)動機中仍然很常見,為了解決葉片損壞的故障就必須在計算它的動態(tài)靜態(tài)應力和應變場分布振動特性的同時深入了解葉片的動態(tài)靜態(tài)特性。為了計算葉片的力學性能,首先必須了解葉片幾何特性,而葉片的幾何特性指的是與葉片形狀和尺寸有關的數(shù)據(jù),如葉片的截面面積、質(zhì)心坐標、葉片的最大厚度、葉片的弦長等,而曲線的擬合是計算以上參數(shù)的必備前提。

    2 ?葉片點云的輸入與處理

    點云模型的獲取是一切研究的開端,也是逆向工程的開端。該章主要會詳細介紹點云從獲取、初步處理、數(shù)據(jù)擬合等方面的具體過程,以及相關的概念,為接下來的軸截面幾何特性奠定基礎。

    2.1 點云的輸入

    先提取出葉片坐標點數(shù)據(jù),對于復雜曲面數(shù)據(jù)的獲取方法初步計劃為才用數(shù)字化測量,可分為兩大類:接觸式測量和非接觸式測量[5]。

    其主要原理和特點如下。

    2.1.1 觸發(fā)式數(shù)據(jù)測量的原理

    當采樣探頭首先接觸到樣品表面時,探頭尖端因為產(chǎn)生力發(fā)生細微形變,從而使采樣開關被觸發(fā),這也就可以讓探頭記錄下當前的坐標,并通過每一個點逐次采樣移動來采集樣品表面的坐標數(shù)據(jù)。這個方式具有非常高的精確度,但是效率低下,用來測算軟性材料通常具有巨大誤差,并且不宜測量凹型曲面,該方法常用于測量回轉體的表面數(shù)據(jù)。

    2.1.2 連續(xù)式數(shù)據(jù)測量的原理

    使用模擬開關探頭來進行連續(xù)數(shù)據(jù)的測量,其實是利用懸掛在三維彈簧系統(tǒng)中的探針位置偏移所產(chǎn)生的電感或者電容的變化,與此同時對機電模擬量進行轉換。當探針以切向速度沿樣品表面移動時,它就會發(fā)出與坐標位置相關偏移量相對應的電流信號或者電壓信號。由于采集過程是連續(xù)的,所以這個方法的速度比點出探針要快很多,所以會有比較高的采樣精度,可以用來進行大量的數(shù)據(jù)采集工作。另外,因為接觸壓力很小,小直徑探頭可以掃描模型的細小零件或者軟材料模型。

    該文所研究的點云大概提取到了137萬個點,之后便可以對收集到的葉片點云利用Solidwork或者Imageware進行初步的還原并且可初步觀察一下葉片的形狀和與坐標軸的位置關系。

    2.2 點云的降噪(誤差點的處理)

    點云取樣時產(chǎn)生的誤差點(即指噪點,以下稱作誤差點)關系到葉片截面幾何特性計算的精度,也直接影響在截面幾何特性修復過程之后的維修工作,決定了修復之后的飛機發(fā)動機葉片的空氣動力學性能。該章將對獲取點云時產(chǎn)生的誤差以及誤差對截面幾何特性計算時的影響進行分析。

    2.2.1 誤差點產(chǎn)生的主要形式

    (1)采集數(shù)據(jù)發(fā)生缺失或者數(shù)據(jù)的密度不符合要求。利用這種不完整的數(shù)據(jù)來進行擬合會出現(xiàn)較大的誤差,很難達到所需測量精度的要求。

    (2)當被測物體是大工件或透明件時,通常會發(fā)生同一個表面采集的數(shù)據(jù)結果得到的是多層復合點云。

    (3)單次的數(shù)據(jù)采集不準確而對整體的測量精度產(chǎn)生影響。

    (4)由于誤差的過量積累而導致錯誤的測量結果。

    (5)點云在幾何分析上的錯誤拼接而導致測量誤差出現(xiàn)。

    (6)測量數(shù)據(jù)中的噪點過多。

    2.2.2 誤差點對葉片截面幾何特性分析的影響

    誤差點會導致截面曲線擬合上的出錯,以至于無法完成平滑曲線的擬合。在面積的插值算法上也會導致輸出面積結果錯誤。此外,在定量分析中,其分析結果一般不是直接通過一步測算就可以得到的,而是在經(jīng)過多次的測量后通過計算得到的。每次的測量都存在一定的誤差,這些誤差均會帶入分析結果中并產(chǎn)生影響,這些誤差是由傳遞性的,這將導致點云失真而導致所有計算的幾何參數(shù)失去可信度。所以,在點云處理之前需要進行消除誤差點的預處理。

    2.2.3 誤差點的識別和去除

    在所獲取的點云中存在的“跳點“以及”壞點“對點云數(shù)據(jù)使用MATLAB進行葉片軸截面幾何特性運算的結果影響相當大。如果不把這些誤差點去除將會使計算結果還有誤差分析的結果的可信度大打折扣。一般的,誤差點的判斷方法有以下幾個。

    (1)直接觀察法。通過圖形終端將偏離點云圖像較大的點或者明顯的孤立點直接去除。此方法比較適合作點云數(shù)據(jù)的初步檢查,可以從點云點集中直接篩選出一些有著明顯誤差的異常點。

    (2)利用曲線擬合檢查的方法。通過提取出分析的截面點云的起點和終點而得到的曲線采用最小二乘法來進行擬合,曲線擬合所需的階次可根據(jù)飛機發(fā)動機葉片的具體形狀來選?。ㄒ话闱闆r選取3~4階)然后分別計算中間數(shù)據(jù)點到擬合曲線的歐氏距離,若點P的∣e∣≥[ε],其中[ε]為給定的允差,則認定點P是壞點,應給去除。

    (3)弦高差法。連接要檢查的前后兩個點,計算P到弦的距離,若點P的∣e∣≥[ε],其中[ε]為給定的允差,則認為P是壞點,應該去除。這種方法適合用于點云均衡分布但是分布比較密集的情況,特別是在點云曲率明顯較大的情況。

    2.3 點云的初步處理

    在點云數(shù)據(jù)獲取之后,將點云數(shù)據(jù)導入Imageware可初步觀察到點云姿態(tài)。之后在Imageware中利用斜線的梯度化現(xiàn)象以地面為基準擺正工件的榫頭(即擺正葉片點云)。擺正時要從正視圖、左視圖、右視圖分別調(diào)整,最終得出一個擺放工整的工件點云。

    由于后續(xù)工作都要在正位點云的狀態(tài)下進行的,所以正位之后的點云的坐標尤其重要。在微調(diào)擺正之后利用Imageware導出正位點云坐標,并另存為Blade.txt格式文本文檔等待后續(xù)處理。

    在Imageware環(huán)境下觀察葉片點云的具體形狀,并確定整個葉片點云中的葉片部分還有底座榫頭部分。通過Imageware對點云邊界點的分析功能可知,該葉片樣品在Z軸坐標15.00以上的部分為葉片部分,其以下的部分為榫頭部分。我們主要的研究任務是以葉片的軸截面的幾何特性作為研究對象,所以我們以Z軸坐標15.00及其以上部分作為主要的研究對象來進行分組以及處理。

    之后我們利用SPSS數(shù)學統(tǒng)計軟件的篩選數(shù)據(jù)的功能對點云數(shù)據(jù)進行篩選。SPSS有著強大的數(shù)據(jù)管理功能,其中包括數(shù)據(jù)的排序、篩選和轉置等,其中篩選功能非常強大,可以根據(jù)用戶的需要進行自定義篩選,也可以根據(jù)各類函數(shù)、條件進行數(shù)據(jù)篩選并且進行刪除,轉存等工序,而且還可以自動輸出數(shù)據(jù)處理日志以便進行數(shù)據(jù)操作分析。

    2.4 葉片軸截面曲線的擬合

    2.4.1 數(shù)據(jù)點的分組

    將之前用SPSS提取出來的一個點云薄片除去Z軸坐標之后近似看作一個平面上的點,這樣就可以對該提取點云進行擬合處理。實際上,葉片的截面曲線是由兩條型線組成的,分別是上型線和下型線,這個曲線結果通常是由流場分析得出的。

    因此,在進行曲線擬合之前必須要把之前得到的數(shù)據(jù)對分成兩組,也就是一組上型線點和一組下型線點。具體的實施方式如下所述。

    (1)利用SPSS篩選出所截取的數(shù)據(jù)點中,X坐標的兩個極點,即最大和最小的兩個點,此兩點分別記為Dmax和Dmin。

    (2)再利用點到距離的算法篩選出到Dmax和Dmin的連線的距離最遠的點Dfar。

    (3)根據(jù)坐標分別推算出Dmax和Dfar的連線的直線方程以及Dmin和Dfar的連線的直線方程。

    (4)對(1)所篩選出來的點云的其余點與(3)中所聯(lián)立出的兩條直線方程分別進行比較,這兩條線上方的點屬于葉片剖面輪廓線的上輪廓即上型線,而此兩線下方的點屬于葉片剖面輪廓線的下輪廓即下型線。Dfar是位于截面輪廓上型線的點,Dmax和Dmin既是上型線的點又是下型線的點,它們是葉片截面上型線和下型線的分界點。

    如此便可將提取出來的數(shù)據(jù)點分為兩組,以便接下來對整個葉片截面曲線進行擬合研究。

    2.4.2 葉片截面曲線的擬合

    對于飛機發(fā)動機葉片截面曲線擬合一般才用的方法有兩種,分別如下。

    (1)最小二乘多項式法。

    當在研究兩個自變量(x,y)之間的關系時,我們通常能夠得到一系列數(shù)據(jù)對(x1,y1,x2,y2...xn,yn);在把上述數(shù)據(jù)對在X-Y直角坐標系里繪制時,如果這些數(shù)據(jù)點在一條直線附近時,就可以得出線性方程式,如式(1)。

    (1)

    其中:a0、a1是任意實數(shù)。

    為了建立這個線性方程,我們必須確定a0和a1,那么我們就可以應用“最小二乘法原理”,把實測值Yi和使用方程式(1)的計算值Yj=a0+a1X的離差(Yi-Yj)的平方和小作為“優(yōu)化準則”。

    令:φ= ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?(2)

    把方程式(1)代入方程式(2)中得:

    φ= ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?(3)

    在的值為最小值時,可以用上述函數(shù)式(1-003)計算其對a0、a1的偏導數(shù),使這兩個偏導數(shù)的數(shù)值為0。

    (4)

    (5)

    亦即: ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?(6)

    (7)

    如此我們就得到了兩個關于a0、a1兩個方程組,通過對這兩個方程組求解,我們可以得到以下結果:

    (8)

    (9)

    此時,將a0、a1代進式(1)式子里去,那么式(1)便為回歸的元線性方程,也就是數(shù)學模型。

    在回歸過程中,其回歸本身的相關性也不會精準地通過每個回歸數(shù)據(jù)點(x1,y1,x2,y2...,xm,ym),為了確定其關聯(lián)性的質(zhì)量,可以用相關系數(shù)“R”、統(tǒng)計量“F”以及剩余的標準差“S”進行判斷;“R”的值越接近1越佳;“F”的絕對值越大越適合;“S”越逼近0越相關。

    (10)

    在方程式(1)中,m是樣本本身的體量即容量的大小,也就是實驗進行的次數(shù),而Xi、Yi是任意一組實驗數(shù)據(jù)對應的X、Y的值[6]。

    (2)正交多項式配合回歸曲線法。

    設存在一組觀察值(xt,yt)t=1,2,…,n,那么如果該觀察值存在非線性關系,那么其多項式回歸方程為:

    (11)

    為了讓離回歸平方和,也就是SSQ=∑(y-)2最小,可以從最小二乘法原理推導出下面的正規(guī)方程組:

    (12)

    解上述方程組可得:b0,b1,b2,…,bp。

    若令x1=x,x2=x2,…,xp=xp,或φ1(x)=x,φ2(x)=x2,…,φp(x)=xp,則(11)可改寫成 :

    (13)

    (14)

    如此,將xi或Φi(x)視為一個新變量,(13)或(14)方程式就為一個元數(shù)為p的多元線性回歸方程,每一個偏回歸系數(shù)di仍然可以根據(jù)下列正規(guī)方程組求得相應的解。

    (15)

    其中,(i,j=1,2,…,p)。

    同樣,對于多元多項式回歸,也可以化為多元線性回歸來分析,例如,對于多變量的任意多項式回歸方程:

    (16)

    只要令x1=z1,x2=z2,x3=z12,x4=z1z2,x5=z22…可化為多元線性回歸方程:

    (17)

    其偏回歸系數(shù)的計算,回歸方程的顯著性檢驗,各偏回歸平方和的計算及顯著性檢驗,都與多元線性回歸分析相似。

    一般來講對于飛機發(fā)動機葉片截面的擬合,一般才用的是最小二乘多項式法擬合,也有才用正交多項式配合回歸發(fā)動機性能曲線擬合以減少計算工作量。因為計算和測量所存在的誤差是不可避免的,我們獲得的數(shù)據(jù)也不可能準確地反映真實的結果。因此,用上述方法擬合實驗結果曲線或者是流場計算結果是合理的。所以最小二乘法是相對較好的處理選擇,最小二乘法準則可以使測量值與其擬合值之差的加權平方和最小,如此可以盡量消除誤差的影響。擬合的曲線是一條盡可能平滑地連接每個測量點的曲線。假定擬合曲線必需準確經(jīng)過型值點,則最小二乘多項式是有點艱難的。因為只有當擬合的階數(shù)很高時,我們所擬合出的曲線才能穿過全部的類型值點,但該擬合曲線在相鄰的類型值點之間會劇烈波動,所以由此得到的擬合曲線也不合理。因此,在需要精確傳遞類型值點的曲線擬合中,最常用的是三次樣條函數(shù)。三次樣條函數(shù)是一種分段代數(shù)多項式,它簡單易擬合,擬合能力強。它可以通過平面上任意數(shù)量的類型點,保證曲線的第一導數(shù)和第二導數(shù)的連續(xù)性。所以對于某些曲線,擬合三次插值的效果可能優(yōu)于三次樣條插值 。

    在曲線擬合中,MATLAB給出了多種擬合措施,以便于積分計算,同時并考慮了曲線擬合的精度,所以我們使用傅里葉級數(shù)極限逼近曲線擬合。葉片截面點云分布非常復雜,所以我們利用傅立葉級數(shù)逼近式:

    f(t)=a0+(b1sinωt+a1cosωt)+(b2sin2ωt+a2cos2ωt)+…+(bksinkωt+akcoskωt) ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? (18)

    作為葉片截面曲線擬合的依據(jù)。多次試驗并且取誤差最小的擬合方程,最后選擇了有限項傅里葉級數(shù)作為擬合的方法,并且根據(jù)這種擬合方法編寫了MATLAB程序,利用cftool工具箱對截面輪廓的上型線進行了擬合(如圖2-11)。

    由于下型線近似為一條直線,所以我們用兩條直線來進行擬合,上下型線共同構成封閉的截面曲線。

    采用上述的方法,我們就可以對之前分組的葉片截面曲線進行擬合(如圖2-12)。

    參考文獻

    [1] 白瑞金,張利國.渦輪葉片修復及其市場分析.北京:北京航空制造工程研究所,2002.

    [2] 陳裕芹.反求工程在發(fā)動機葉片檢測中的應用研究:(碩士學位論文).廣東:廣東工業(yè)大學,2011.

    [3] 石辛民,郝整清. 基于 MATLAB 的實用數(shù)值計算[M]. 北京:清華大學出版社,2006(2)

    [4] (美)Delores M. Etter David C. Kuncicky Holly Moore 著,邱李華譯. MATLAB 7 級工程問題解方案[M]. 北京:機械工業(yè)出版社,2006(3).

    [5] 呂秋娟,侯根良.反求工程中的幾種數(shù)字化測量方法.西安:第二炮兵工程學院,2000.

    [6] 王穗輝.誤差理論與測量平差:同濟大學出版社,2010.

    [7] 蔡山,張浩,陳洪輝,沙基昌. 基于最小二乘法的分段三次曲線擬合方法研究. 科學技術與工程[J],2007(2)

    [8] R.A.WalliS.Axial Flow Fans and Ducts.John wiley and Sons,Inc,1983.

    猜你喜歡
    曲線擬合渦輪發(fā)動機
    2014款寶馬525Li渦輪增壓壓力過低
    發(fā)動機空中起動包線擴展試飛組織與實施
    曲線擬合的方法
    價值工程(2017年31期)2018-01-17 00:34:27
    基于曲線擬合的投棄式剖面儀電感量算法
    電子測試(2017年12期)2017-12-18 06:35:46
    Matlab曲線擬合工具箱在地基沉降預測模型中的應用
    Matlab曲線擬合法在地基沉降預測中的應用
    渦輪增壓發(fā)動機與雙離合變速器的使用
    新一代MTU2000發(fā)動機系列
    Opel公司新型1.0L渦輪增壓直接噴射汽油機
    新型1.5L-Eco-Boost發(fā)動機
    aaaaa片日本免费| 国产精品亚洲一级av第二区| 老师上课跳d突然被开到最大视频| 俺也久久电影网| 亚洲一级一片aⅴ在线观看| 国产成人精品久久久久久| 国产精品亚洲美女久久久| 淫妇啪啪啪对白视频| 免费看av在线观看网站| 精品欧美国产一区二区三| 亚洲av熟女| 嫩草影院精品99| 久久精品国产亚洲av香蕉五月| av中文乱码字幕在线| 亚洲真实伦在线观看| 亚洲成人久久性| 国产免费男女视频| 青春草视频在线免费观看| 久久久国产成人精品二区| av女优亚洲男人天堂| 激情 狠狠 欧美| 精品国产三级普通话版| 我的老师免费观看完整版| 中文资源天堂在线| 久久精品国产鲁丝片午夜精品| 简卡轻食公司| 日本a在线网址| 国产精品国产三级国产av玫瑰| 人人妻人人澡欧美一区二区| 亚洲专区国产一区二区| 波多野结衣巨乳人妻| 色尼玛亚洲综合影院| 国产美女午夜福利| 亚洲国产精品成人综合色| 十八禁网站免费在线| 少妇人妻一区二区三区视频| 女同久久另类99精品国产91| 国产亚洲精品av在线| 欧美成人免费av一区二区三区| 国产片特级美女逼逼视频| 69人妻影院| 亚洲av成人av| 精品久久久久久久久亚洲| 国产精品一区二区性色av| 看免费成人av毛片| 联通29元200g的流量卡| 国产精品久久久久久久久免| 国产真实乱freesex| 六月丁香七月| 亚洲最大成人av| 国产成人91sexporn| 亚洲成人精品中文字幕电影| 亚洲中文字幕日韩| 一个人看视频在线观看www免费| 丝袜喷水一区| 男人狂女人下面高潮的视频| 国产片特级美女逼逼视频| 一区二区三区四区激情视频 | 国产精品国产高清国产av| 一级毛片aaaaaa免费看小| 精品日产1卡2卡| 91精品国产九色| 久久精品国产亚洲av天美| 熟女人妻精品中文字幕| 俺也久久电影网| 国产片特级美女逼逼视频| 亚洲精品国产av成人精品 | 99久久精品一区二区三区| 天堂网av新在线| 日韩欧美国产在线观看| 欧美极品一区二区三区四区| 国产视频内射| 色av中文字幕| 午夜福利成人在线免费观看| 老司机福利观看| 桃色一区二区三区在线观看| 久久精品夜色国产| 最后的刺客免费高清国语| 欧美激情在线99| 欧美性感艳星| 免费观看人在逋| 女人十人毛片免费观看3o分钟| 91久久精品国产一区二区成人| 欧美一区二区亚洲| 久久久久久久午夜电影| 亚州av有码| 国产精品国产三级国产av玫瑰| 亚州av有码| 久久久久久大精品| 亚州av有码| 亚洲成av人片在线播放无| 中文字幕久久专区| 级片在线观看| 亚洲不卡免费看| 蜜桃亚洲精品一区二区三区| 久久久久性生活片| 韩国av在线不卡| 又爽又黄无遮挡网站| 国产亚洲欧美98| 精品无人区乱码1区二区| 久久久国产成人免费| 国产亚洲精品久久久com| 久久九九热精品免费| 国产亚洲av嫩草精品影院| 99久久精品国产国产毛片| 最近手机中文字幕大全| 给我免费播放毛片高清在线观看| 性色avwww在线观看| 在线观看午夜福利视频| 69av精品久久久久久| 最近中文字幕高清免费大全6| 偷拍熟女少妇极品色| 有码 亚洲区| 一级毛片电影观看 | 国产一区亚洲一区在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品色激情综合| 在线观看一区二区三区| 99久久中文字幕三级久久日本| 国产精品国产三级国产av玫瑰| 久久中文看片网| 精品一区二区免费观看| 美女免费视频网站| 亚洲美女黄片视频| 搡老熟女国产l中国老女人| 老女人水多毛片| 国产淫片久久久久久久久| 日韩欧美精品免费久久| 六月丁香七月| 久久久色成人| 国产精品美女特级片免费视频播放器| 国产黄色视频一区二区在线观看 | 国内精品一区二区在线观看| 内射极品少妇av片p| 人妻制服诱惑在线中文字幕| 亚洲最大成人中文| 日韩国内少妇激情av| 两性午夜刺激爽爽歪歪视频在线观看| 看黄色毛片网站| 日韩在线高清观看一区二区三区| 麻豆一二三区av精品| 国产伦一二天堂av在线观看| 国产毛片a区久久久久| 黄色配什么色好看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲最大成人中文| 人妻制服诱惑在线中文字幕| 国产精品乱码一区二三区的特点| 日韩国内少妇激情av| 又黄又爽又免费观看的视频| 看免费成人av毛片| 大香蕉久久网| 日韩,欧美,国产一区二区三区 | 亚洲18禁久久av| 国产精品一区二区三区四区久久| 亚洲18禁久久av| 成人漫画全彩无遮挡| 搡老妇女老女人老熟妇| 亚洲aⅴ乱码一区二区在线播放| 六月丁香七月| 久久九九热精品免费| 久久精品夜色国产| 久久精品国产99精品国产亚洲性色| 欧美国产日韩亚洲一区| av专区在线播放| 精品少妇黑人巨大在线播放 | 国产一区二区激情短视频| 观看美女的网站| 久久韩国三级中文字幕| 一边摸一边抽搐一进一小说| 男人狂女人下面高潮的视频| 欧洲精品卡2卡3卡4卡5卡区| 看黄色毛片网站| 久久久久久九九精品二区国产| 免费不卡的大黄色大毛片视频在线观看 | 亚洲成a人片在线一区二区| 日韩亚洲欧美综合| 欧美xxxx性猛交bbbb| 国产视频内射| 中出人妻视频一区二区| 老司机午夜福利在线观看视频| 国产视频内射| 国产v大片淫在线免费观看| 国产精品久久电影中文字幕| 中文在线观看免费www的网站| 精品久久久噜噜| 美女被艹到高潮喷水动态| 中国美白少妇内射xxxbb| 日韩,欧美,国产一区二区三区 | 最后的刺客免费高清国语| 狂野欧美白嫩少妇大欣赏| 色视频www国产| 级片在线观看| 国产av麻豆久久久久久久| 日日摸夜夜添夜夜添av毛片| 国语自产精品视频在线第100页| 熟女电影av网| 久久久午夜欧美精品| 亚洲美女视频黄频| 亚洲国产欧美人成| 久久综合国产亚洲精品| 男女视频在线观看网站免费| 秋霞在线观看毛片| 欧美最黄视频在线播放免费| 一区二区三区高清视频在线| 国产成年人精品一区二区| 好男人在线观看高清免费视频| 国产精品野战在线观看| 神马国产精品三级电影在线观看| 国产成人福利小说| 国产高清有码在线观看视频| 在线免费十八禁| 99热这里只有是精品在线观看| 丝袜喷水一区| 成人特级黄色片久久久久久久| 久久久午夜欧美精品| 久久亚洲精品不卡| 欧美性猛交╳xxx乱大交人| 欧美成人一区二区免费高清观看| 男人狂女人下面高潮的视频| 两性午夜刺激爽爽歪歪视频在线观看| av中文乱码字幕在线| a级毛片免费高清观看在线播放| 久久久色成人| 欧美极品一区二区三区四区| 啦啦啦观看免费观看视频高清| 99热精品在线国产| 免费无遮挡裸体视频| 热99在线观看视频| 禁无遮挡网站| 国产av麻豆久久久久久久| 欧美又色又爽又黄视频| 我要搜黄色片| 美女 人体艺术 gogo| 国产精品亚洲美女久久久| 国内少妇人妻偷人精品xxx网站| 日本 av在线| av在线播放精品| 别揉我奶头~嗯~啊~动态视频| 天堂网av新在线| 国产精品1区2区在线观看.| 老熟妇仑乱视频hdxx| 亚洲精品一区av在线观看| 香蕉av资源在线| 男插女下体视频免费在线播放| 观看免费一级毛片| 国产精品一区二区三区四区免费观看 | 精品少妇黑人巨大在线播放 | 日韩三级伦理在线观看| 女生性感内裤真人,穿戴方法视频| 在线观看午夜福利视频| 干丝袜人妻中文字幕| 国内精品久久久久精免费| 伦理电影大哥的女人| 精华霜和精华液先用哪个| 亚洲熟妇中文字幕五十中出| 欧美xxxx黑人xx丫x性爽| a级毛色黄片| 丝袜喷水一区| 亚洲最大成人中文| 日韩欧美 国产精品| 在线观看一区二区三区| 国产黄a三级三级三级人| 两个人视频免费观看高清| 国产亚洲精品久久久com| 蜜桃亚洲精品一区二区三区| 精品一区二区三区av网在线观看| 亚洲无线观看免费| 夜夜看夜夜爽夜夜摸| 国产 一区精品| 少妇裸体淫交视频免费看高清| 成年女人毛片免费观看观看9| 亚洲中文日韩欧美视频| av福利片在线观看| 欧美日韩在线观看h| 国产精品国产三级国产av玫瑰| 我要搜黄色片| 午夜激情福利司机影院| 国产伦一二天堂av在线观看| 村上凉子中文字幕在线| 黄色视频,在线免费观看| 少妇人妻精品综合一区二区 | 亚洲av第一区精品v没综合| 亚洲人与动物交配视频| 中文字幕人妻熟人妻熟丝袜美| 日韩,欧美,国产一区二区三区 | 久久久久九九精品影院| 干丝袜人妻中文字幕| 国产亚洲精品av在线| 最好的美女福利视频网| 国产伦在线观看视频一区| 久久热精品热| 久久久久久久久久黄片| 国产精品久久视频播放| 国产精品伦人一区二区| 欧美日本视频| 久久久午夜欧美精品| 国产一区亚洲一区在线观看| 精品久久久久久久末码| 女人被狂操c到高潮| 99久久九九国产精品国产免费| 在线观看一区二区三区| 日本一二三区视频观看| 久久久久免费精品人妻一区二区| 日本 av在线| 免费人成视频x8x8入口观看| 欧洲精品卡2卡3卡4卡5卡区| 两个人的视频大全免费| 麻豆国产97在线/欧美| 日韩 亚洲 欧美在线| 搡老熟女国产l中国老女人| 亚洲成人精品中文字幕电影| 欧美一区二区亚洲| 欧美3d第一页| 国产视频一区二区在线看| 国产一区二区亚洲精品在线观看| 国产中年淑女户外野战色| 亚洲国产日韩欧美精品在线观看| 久久精品久久久久久噜噜老黄 | 少妇人妻一区二区三区视频| 日韩欧美免费精品| 久久精品国产亚洲网站| 国产成人一区二区在线| 草草在线视频免费看| 亚洲欧美中文字幕日韩二区| 久久久久免费精品人妻一区二区| 久久人人精品亚洲av| 欧美zozozo另类| 国产aⅴ精品一区二区三区波| 精品人妻偷拍中文字幕| 久久久久久久午夜电影| 99久国产av精品国产电影| av.在线天堂| 美女免费视频网站| 麻豆国产97在线/欧美| videossex国产| 国产单亲对白刺激| 天堂av国产一区二区熟女人妻| 小说图片视频综合网站| 国产男人的电影天堂91| 亚洲成人久久爱视频| 欧美不卡视频在线免费观看| 草草在线视频免费看| 99国产精品一区二区蜜桃av| 校园春色视频在线观看| 成年版毛片免费区| 美女大奶头视频| 男女之事视频高清在线观看| 最近2019中文字幕mv第一页| 人妻丰满熟妇av一区二区三区| 免费av不卡在线播放| eeuss影院久久| 国产精品精品国产色婷婷| 国产黄色视频一区二区在线观看 | 国产乱人视频| 天堂√8在线中文| 久久午夜亚洲精品久久| 午夜久久久久精精品| 国产精品国产高清国产av| 免费看a级黄色片| 91久久精品国产一区二区三区| 国产亚洲精品久久久com| 99久国产av精品| 欧美高清成人免费视频www| 日韩精品青青久久久久久| 长腿黑丝高跟| 99在线视频只有这里精品首页| 老女人水多毛片| 亚洲中文字幕一区二区三区有码在线看| 国产一区二区三区av在线 | 亚洲av.av天堂| 日日干狠狠操夜夜爽| 国产精品美女特级片免费视频播放器| 亚洲丝袜综合中文字幕| 国产极品精品免费视频能看的| 国产高清激情床上av| 成人午夜高清在线视频| 亚洲精品乱码久久久v下载方式| 国产精品不卡视频一区二区| 91午夜精品亚洲一区二区三区| 欧美精品国产亚洲| 美女cb高潮喷水在线观看| 国产精品伦人一区二区| 欧美+日韩+精品| 国产精品无大码| 最近中文字幕高清免费大全6| 国产黄色视频一区二区在线观看 | 麻豆国产97在线/欧美| 久久国内精品自在自线图片| 中文字幕免费在线视频6| 内射极品少妇av片p| 2021天堂中文幕一二区在线观| 亚洲精品国产成人久久av| 免费在线观看成人毛片| 尾随美女入室| 3wmmmm亚洲av在线观看| 午夜亚洲福利在线播放| 国产精品久久久久久久久免| 国产精品99久久久久久久久| avwww免费| 久99久视频精品免费| 国产精品不卡视频一区二区| 久久这里只有精品中国| 可以在线观看的亚洲视频| 内射极品少妇av片p| 国产精品一区二区三区四区久久| 亚洲成人久久爱视频| 少妇熟女aⅴ在线视频| 日韩欧美一区二区三区在线观看| 好男人在线观看高清免费视频| 成人二区视频| 波野结衣二区三区在线| 男女边吃奶边做爰视频| 99久国产av精品| 又黄又爽又免费观看的视频| 久久婷婷人人爽人人干人人爱| 美女内射精品一级片tv| 精品久久久久久久久久久久久| 精品不卡国产一区二区三区| 香蕉av资源在线| 国产亚洲av嫩草精品影院| 国产探花在线观看一区二区| 国内少妇人妻偷人精品xxx网站| 精品久久久久久久久av| 一级黄片播放器| 国产伦一二天堂av在线观看| 美女免费视频网站| 一区二区三区免费毛片| 中文在线观看免费www的网站| 午夜福利在线观看免费完整高清在 | av女优亚洲男人天堂| 小蜜桃在线观看免费完整版高清| 18禁黄网站禁片免费观看直播| 午夜免费激情av| av专区在线播放| 久久精品国产亚洲av香蕉五月| 国产精品久久久久久久久免| 国产在线精品亚洲第一网站| 嫩草影视91久久| 麻豆av噜噜一区二区三区| 亚洲丝袜综合中文字幕| 亚洲精品亚洲一区二区| 亚洲熟妇中文字幕五十中出| 国产日本99.免费观看| 成人亚洲欧美一区二区av| 男女边吃奶边做爰视频| 国产日本99.免费观看| 亚洲精华国产精华液的使用体验 | 国产免费男女视频| 男人舔女人下体高潮全视频| 日韩av在线大香蕉| 亚洲在线观看片| 婷婷精品国产亚洲av在线| 乱人视频在线观看| 麻豆精品久久久久久蜜桃| 国产精品99久久久久久久久| 18禁黄网站禁片免费观看直播| 日韩高清综合在线| 小蜜桃在线观看免费完整版高清| 亚洲乱码一区二区免费版| 两个人视频免费观看高清| 国产精品一区二区免费欧美| 蜜臀久久99精品久久宅男| 免费看美女性在线毛片视频| 精品午夜福利视频在线观看一区| 在线a可以看的网站| 舔av片在线| 99久久成人亚洲精品观看| 国内少妇人妻偷人精品xxx网站| 亚洲精品粉嫩美女一区| 国产中年淑女户外野战色| 久久久久久久久久成人| 自拍偷自拍亚洲精品老妇| 99热只有精品国产| 一本久久中文字幕| 啦啦啦观看免费观看视频高清| 久久精品人妻少妇| 国产男人的电影天堂91| 亚洲av免费高清在线观看| 一进一出抽搐动态| 午夜久久久久精精品| av天堂中文字幕网| 日本精品一区二区三区蜜桃| 非洲黑人性xxxx精品又粗又长| 又黄又爽又免费观看的视频| 丰满的人妻完整版| 国产成人aa在线观看| 亚洲最大成人中文| 在线观看一区二区三区| 亚洲四区av| 91狼人影院| 亚洲欧美成人综合另类久久久 | 看黄色毛片网站| 日本 av在线| 日本一二三区视频观看| 精品久久久久久成人av| 免费av不卡在线播放| av国产免费在线观看| 免费观看的影片在线观看| 亚洲精华国产精华液的使用体验 | 成人三级黄色视频| 最近在线观看免费完整版| 色5月婷婷丁香| 国产高清激情床上av| 国产精华一区二区三区| 内射极品少妇av片p| 亚洲美女视频黄频| 夜夜夜夜夜久久久久| 日本黄大片高清| 欧美在线一区亚洲| 成人欧美大片| 日韩精品有码人妻一区| 国产美女午夜福利| 成人亚洲欧美一区二区av| 露出奶头的视频| 中文字幕人妻熟人妻熟丝袜美| 免费看a级黄色片| 一级毛片我不卡| 婷婷亚洲欧美| 99久久无色码亚洲精品果冻| 国产一区亚洲一区在线观看| 欧美色欧美亚洲另类二区| 看十八女毛片水多多多| 国国产精品蜜臀av免费| 亚洲精品国产av成人精品 | 午夜福利在线在线| 毛片女人毛片| 一个人看的www免费观看视频| 91久久精品国产一区二区三区| 天堂影院成人在线观看| 18禁黄网站禁片免费观看直播| 麻豆av噜噜一区二区三区| 国产女主播在线喷水免费视频网站 | 亚洲国产精品久久男人天堂| 国产爱豆传媒在线观看| 性插视频无遮挡在线免费观看| 亚洲av中文字字幕乱码综合| 国产精品av视频在线免费观看| 桃色一区二区三区在线观看| 悠悠久久av| 1024手机看黄色片| 亚洲中文字幕日韩| 99国产精品一区二区蜜桃av| 搡老熟女国产l中国老女人| 国产亚洲91精品色在线| 九九爱精品视频在线观看| 日日啪夜夜撸| 成人性生交大片免费视频hd| 国产精品久久久久久久电影| 两个人视频免费观看高清| 小蜜桃在线观看免费完整版高清| 国产精品免费一区二区三区在线| 蜜桃亚洲精品一区二区三区| 国产精品福利在线免费观看| 国产精品一区www在线观看| 国产视频一区二区在线看| 亚洲国产精品成人综合色| 91在线观看av| 男女边吃奶边做爰视频| 国产高清三级在线| 久久久成人免费电影| 国产91av在线免费观看| 日韩中字成人| 女人被狂操c到高潮| 大型黄色视频在线免费观看| 一本久久中文字幕| 亚洲av五月六月丁香网| 俄罗斯特黄特色一大片| 国产成人a区在线观看| 国产高清激情床上av| 亚洲一区高清亚洲精品| 亚洲国产高清在线一区二区三| 男女视频在线观看网站免费| 天天一区二区日本电影三级| 在线观看美女被高潮喷水网站| 在线国产一区二区在线| 少妇的逼好多水| 国产高清三级在线| av免费在线看不卡| 直男gayav资源| 国内揄拍国产精品人妻在线| 国产单亲对白刺激| 91午夜精品亚洲一区二区三区| 免费观看的影片在线观看| 免费av不卡在线播放| 成人午夜高清在线视频| 又黄又爽又免费观看的视频| 一区二区三区四区激情视频 | 又爽又黄a免费视频| 一区福利在线观看| 免费黄网站久久成人精品| 国产不卡一卡二| 亚洲av电影不卡..在线观看| 国产精品福利在线免费观看| 亚洲av免费高清在线观看| 在线观看免费视频日本深夜| 免费大片18禁| 久久久久久久亚洲中文字幕| 色噜噜av男人的天堂激情| 成人亚洲精品av一区二区| 久久国产乱子免费精品| 人妻少妇偷人精品九色| 老熟妇仑乱视频hdxx| 国产精品美女特级片免费视频播放器| 老熟妇仑乱视频hdxx| 日韩 亚洲 欧美在线| 亚洲精品日韩在线中文字幕 | 久久精品国产亚洲网站| 日韩欧美精品v在线| 精品日产1卡2卡| 色噜噜av男人的天堂激情| 久久精品国产亚洲av天美| 婷婷六月久久综合丁香| 国产久久久一区二区三区| 国产蜜桃级精品一区二区三区|