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

    基于VB的優(yōu)化雙圓弧擬合誤差分析及其算法改進(jìn)

    2014-03-29 10:04:04陳錦昌申藝杰
    圖學(xué)學(xué)報(bào) 2014年1期
    關(guān)鍵詞:法向圓弧斜率

    張 弛, 陳錦昌, 申藝杰

    (1.華南理工大學(xué)設(shè)計(jì)學(xué)院,廣東 廣州 510006;2.華南理工大學(xué)機(jī)械與汽車工程學(xué)院,廣東 廣州 510640)

    在數(shù)控加工中,刀具的路徑一般為直線和圓弧,所以大多數(shù)數(shù)控機(jī)床只能實(shí)現(xiàn)直線插補(bǔ)和圓弧插補(bǔ)。在加工帶有曲線輪廓的工件時(shí),按照刀具的路徑特點(diǎn),需要把曲線轉(zhuǎn)化為由直線或圓弧組成的樣條以方便刀具的加工。進(jìn)行直線插補(bǔ)時(shí)由于刀具在每個(gè)插補(bǔ)段都要經(jīng)歷加速、勻速、減速和停頓,使加工時(shí)間延長(zhǎng)且產(chǎn)生過(guò)多的振動(dòng)和磨損,而使用圓弧插補(bǔ)能實(shí)現(xiàn)勻速加工且加工曲線光滑還能減少加工時(shí)間提高效率,因此,使用圓弧逼近來(lái)代替曲線更為有利[1-2]。圓弧逼近一般分為單圓弧逼近和雙圓弧逼近。單圓弧逼近具有使用圓弧數(shù)少連續(xù)性好等優(yōu)點(diǎn),但同時(shí)穩(wěn)定性差、缺乏靈活性,不能處理有拐點(diǎn)的情況;雙圓弧逼近則具有較高的靈活性和穩(wěn)定性,更加適用于曲線輪廓工件的數(shù)控加工[3]。

    本文主要研究雙圓弧對(duì)曲線的逼近及其優(yōu)化設(shè)計(jì),運(yùn)用VB編程實(shí)現(xiàn)了曲線輪廓工件的雙圓弧擬合以及擬合誤差的分析計(jì)算,使用隨機(jī)數(shù)和判斷函數(shù)實(shí)現(xiàn)了簡(jiǎn)單的雙圓弧擬合優(yōu)化設(shè)計(jì),求出可行的擬合參數(shù)。對(duì)影響誤差的因數(shù)進(jìn)行了討論。并按照所得結(jié)論在優(yōu)化雙圓弧算法[4]基礎(chǔ)上對(duì)算法進(jìn)行了改進(jìn)。

    1 優(yōu)化雙圓弧擬合模型

    1.1 雙圓弧擬合方法

    雙圓弧的擬合方法常用的有3種:第1種是以三角形的內(nèi)心作為雙圓弧的聯(lián)接點(diǎn),稱為內(nèi)心法。第2種是將雙圓弧的聯(lián)接點(diǎn)選在樣條曲線上,稱為交點(diǎn)法。第3種方法是采用改變雙圓弧聯(lián)接點(diǎn)的切線斜率使雙圓弧曲線與樣條曲線之間的法向誤差達(dá)到最小,又稱優(yōu)化雙圓弧擬合。內(nèi)心法計(jì)算簡(jiǎn)便,但誤差不易控制在允許范圍內(nèi);交點(diǎn)法的誤差較小,但是由于要解四次非線性方程,計(jì)算復(fù)雜;優(yōu)化雙圓弧計(jì)算簡(jiǎn)便,線性,穩(wěn)定性好,十分適于運(yùn)用到數(shù)控加工。本文以優(yōu)化雙圓弧作為曲線輪廓工件加工時(shí)擬合算法,并通過(guò)調(diào)節(jié)分段數(shù)和每段聯(lián)接點(diǎn)斜率來(lái)減小誤差,達(dá)到擬合的優(yōu)化設(shè)計(jì)。為了找出較適合的斜率,將使用隨機(jī)數(shù)生成斜率值并計(jì)算相應(yīng)的誤差,通過(guò)誤差與最大允許誤差的比較選出適合的斜率。

    1.2 優(yōu)化雙圓弧擬合算法

    設(shè)由測(cè)量或?qū)嶒?yàn)所得的一組離散有序的點(diǎn)列Pti(xi,yi)(i=1,2,…,n)形成一條樣條曲線,稱目標(biāo)曲線。設(shè)目標(biāo)曲線相鄰兩個(gè)節(jié)點(diǎn)為Pti、Pti+1(i=1,2,…,n),現(xiàn)在不妨在曲線上取兩點(diǎn)來(lái)做研究。設(shè)已知兩節(jié)點(diǎn)坐標(biāo)為P1(x1,y1),P2(x2,y2),經(jīng)過(guò)P1、P2點(diǎn)的曲線切線斜率為T1、T2(圖1),并設(shè)T1>T2,兩切線的交點(diǎn)為P3,設(shè)ΔP1P3P2三邊長(zhǎng)a=,b=、c=,若已經(jīng)取定一數(shù)值TN(T1>TN>T2)作為雙圓弧擬合公切點(diǎn)的切線斜率,則有

    圖1 優(yōu)化雙圓弧擬合基本原理

    (1)過(guò)1P作斜率為TN的直線交P23P于P4,此時(shí)我們可以求出已知量d=、e=;

    (2)假設(shè)已經(jīng)找到雙圓弧公切點(diǎn)P5及公切線P6P7,注意P6∈P1P3、P7∈P2P3,則有關(guān)系P6P7//P1P4,下面我們確定P6,進(jìn)而確定P5的位置;

    (3)設(shè)m=,由關(guān)系P6P7//P1P4,由文獻(xiàn)[4]有設(shè)n=,易得

    當(dāng)d>e時(shí)有m>0,可在P6P7上求得符合條件的公切點(diǎn)P5點(diǎn),上式成立。當(dāng)d<e時(shí)則有m<0,這時(shí)P6在P1P3的延長(zhǎng)線上,無(wú)法求出公切點(diǎn)P5。此時(shí)應(yīng)該取P2為基點(diǎn),過(guò)P2點(diǎn)作斜率為TN的直線,同理相反方向作圖即可。若T1<TN<T2,出現(xiàn)d<e的情況同理。

    (4)設(shè)O1、O2分別為兩個(gè)圓弧的圓心,r1、r2為兩圓弧的半徑。有O1P1⊥P1P3、O2P2⊥P2P3、O1P5⊥P6P7。且有P1P6=P6P5、P5P7=P7P2,從而可以作出雙圓弧。

    按照以上結(jié)果,將優(yōu)化雙圓弧擬合算法作為曲線擬合時(shí)求解的數(shù)學(xué)模型。應(yīng)用VB編寫雙圓弧擬合計(jì)算子函數(shù),函數(shù)需要5個(gè)已知變量,在代碼中命名如下:兩個(gè)節(jié)點(diǎn)的坐標(biāo)值(x1,y1)、(x2,y2)及對(duì)應(yīng)的切矢ks和ke,TN值則由聯(lián)接點(diǎn)切線與節(jié)點(diǎn)連線的夾角rr=∠P3P1P2來(lái)間接表示。函數(shù)運(yùn)行后返回輸出對(duì)應(yīng)兩圓弧的圓心坐標(biāo)、聯(lián)接點(diǎn)坐標(biāo)和圓弧半徑,在代碼中表示為xo1,yo1,xo2,yo2,xn,yn,r1,r2。另外,函數(shù)中定義數(shù)組x(i)、y(i)用來(lái)存儲(chǔ)曲線上每段圓弧端點(diǎn)的坐標(biāo)值,用于函數(shù)中的計(jì)算及雙圓弧曲線的繪制。函數(shù)名:Function shuangyuanhu(x1,y1,x2,y2,ks,ke,rr,xo1,yo1,xo2,yo2,xn,yn,r1,r2)。再通過(guò)圓弧繪制命令:Call draw_Arc(x(i),y(i),xn,yn,xo1,yo1)Call draw_Arc(xn,yn,x(i+1),y(i+1),xo2,yo2)繪制兩段圓弧即可實(shí)現(xiàn)曲線輪廓工件的雙圓弧擬合。

    1.3 雙圓弧誤差分析及切線斜率TN選取

    曲線輪廓工件雙圓弧擬合分析誤差時(shí),計(jì)算原曲線和雙圓弧擬合曲線之間的法向誤差:

    以每相鄰兩型值點(diǎn)為一個(gè)單元進(jìn)行討論,可在型值點(diǎn)間取10個(gè)等分點(diǎn),若等分點(diǎn)在以O(shè)1為圓心的圓弧上時(shí),法向誤差為

    若等分點(diǎn)在以O(shè)2為圓心的圓弧上時(shí)

    設(shè)允許的最大誤差為Δ允,當(dāng)ΔX<Δ允時(shí),雙圓弧合格,反之則不合格[5]。若雙圓弧不合格要改變TN值重新擬合直到雙圓弧合格為止。

    2 雙圓弧擬合誤差分析

    2.1 雙圓弧擬合實(shí)例

    例:用數(shù)控車床加工某回轉(zhuǎn)體工件時(shí),設(shè)工件的母曲線方程為二次方程y=0.213157×x-0.50216×10-3×x2其中x∈(0,170)。

    這里不妨先按步長(zhǎng)為10,各段雙圓弧的TN值為隨機(jī)生成的方案進(jìn)行擬合。運(yùn)用編寫的程序畫出雙圓弧擬合曲線圖如圖2所示。

    圖2 雙圓弧擬合實(shí)例

    2.2 誤差分布和簡(jiǎn)單優(yōu)化

    在VB的窗口中繪制法向誤差圖表,以x值為橫坐標(biāo),每組雙圓弧對(duì)應(yīng)的最大誤差值為縱坐標(biāo),畫出誤差分布如圖3所示,并在圖中標(biāo)出允許誤差值Δ允,令Δ允=0.007μm。圖3可以直觀地觀察雙圓弧曲線各部分誤差情況。

    圖3 誤差分布圖

    誤差分布圖顯示出這次的擬合有4段曲線是不合格的,要重新選擇每段的TN值,并重新擬合。設(shè)置誤差判斷條件,當(dāng)誤差不滿足條件時(shí)重新生成TN值,如此重復(fù)循環(huán),直到法向誤差符合要求為止。擬合優(yōu)化程序框圖如圖4所示。

    圖4 擬合優(yōu)化程序框圖

    進(jìn)行TN值的重選和雙圓弧曲線再擬合就可以得出基本符合要求的雙圓弧擬合曲線如圖5所示。

    圖5 重選TN值后的誤差分布圖

    運(yùn)用重選TN值求解雙圓弧曲線計(jì)算量大且常會(huì)出現(xiàn)難以收斂導(dǎo)致程序無(wú)限循環(huán)的問(wèn)題。因此應(yīng)該從誤差的規(guī)律入手尋求更簡(jiǎn)單高效的擬合方案。

    3 影響誤差的因素

    上例是以步長(zhǎng)為10這一條件為前提進(jìn)行研究的,所以TN值成為影響曲線擬合精度的唯一變量。但由于雙圓弧擬合曲線的形狀是由各段基點(diǎn)位置及其TN值共同決定的,所以對(duì)雙圓弧擬合方案求最優(yōu)解將是一個(gè)多元非線性問(wèn)題。這里我們將問(wèn)題簡(jiǎn)化,固定其中一部分條件研究其他條件對(duì)雙圓弧擬合的影響。

    3.1 TN取值的對(duì)誤差的影響

    這里先固定各節(jié)點(diǎn)位置,以便直觀的研究TN值對(duì)擬合誤差的影響。由圖1可見端點(diǎn)及其對(duì)應(yīng)的切矢固定時(shí),取不同TN值將獲得不同的擬合曲線,在ΔP1P3P2包圍的區(qū)域內(nèi)有無(wú)數(shù)解。

    不妨分別按照P1P2斜率和黃金分割率來(lái)確定TN值來(lái)觀察擬合效果。

    (1)平行弦法確定TN值

    平行弦法是指取TN值等于P12P斜率,這時(shí)點(diǎn)4P與2P重合,直線P6P7平行于直線P12P。圖6(a)為上例曲線采用平行弦法雙圓弧擬合時(shí)得出的結(jié)果。該曲線使用平行弦法擬合的最大法向誤差為1.96E-6mm,情況比較理想。

    (2)黃金分割法確定TN值

    在T1,T2間按黃金分割率0.618取值作為TN值,圖6(b)為上例雙圓弧擬合時(shí)得出的結(jié)果。

    圖6 不同TN值對(duì)應(yīng)的誤差分布圖

    TN值確定了雙圓弧擬合曲線的形狀,進(jìn)而決定了誤差的分布。在確定基點(diǎn)的情況下,調(diào)節(jié)TN值可以減小誤差。

    3.2 分段數(shù)對(duì)誤差的影響

    (1)減少分段數(shù),增大步長(zhǎng)

    將上例中步長(zhǎng)10,改為步長(zhǎng)19進(jìn)行雙圓弧擬合,為了減小其他變量的影響,TN按照平行弦法取值,所得結(jié)果如圖7(a)與圖6(a)比較可知雙圓弧擬合的誤差增大,超出了允許值。

    (2)增加分段數(shù),減小步長(zhǎng)

    將上例步長(zhǎng)改為步長(zhǎng)15.4mm進(jìn)行雙圓弧擬合,TN按照平行弦法取值,所得結(jié)果如圖7(b),與圖6(a)比較,誤差減小,擬合精確度提高,法向誤差降低為6.89E-06mm,剛好達(dá)到要求。

    圖7 改變步長(zhǎng)時(shí)誤差分布圖

    經(jīng)過(guò)多次改變步長(zhǎng)獲得多組誤差分布圖后發(fā)現(xiàn)步長(zhǎng)15.4mm是臨界步長(zhǎng),相應(yīng)分段數(shù)為11段。當(dāng)所選步長(zhǎng)小于15.4mm時(shí)雙圓弧曲線都能滿足要求。

    4 改進(jìn)的雙圓弧擬合算法

    4.1 誤差值變化規(guī)律分析

    比較圖6(a)和圖7(a)時(shí)可以發(fā)現(xiàn),盡管所取的分段數(shù)不同,其誤差曲線波動(dòng)趨勢(shì)基本相似。這是由于采用平行弦法確定TN值時(shí),雙圓弧擬合曲線將只由兩節(jié)點(diǎn)及其切矢決定,也就是說(shuō)此時(shí)擬合誤差將只受節(jié)點(diǎn)選取的影響。

    改變步長(zhǎng)就相當(dāng)于改變了節(jié)點(diǎn)的選取。通過(guò)減小步長(zhǎng)來(lái)減小誤差雖然效果明顯,但是常常會(huì)引起大部分曲線段精度過(guò)高同時(shí)過(guò)小的步長(zhǎng)會(huì)使分段數(shù)增多。在具體的曲線擬合中,對(duì)于給定的誤差,往往要求用盡量少的圓弧段數(shù)[6],因此單純通過(guò)減小步長(zhǎng)來(lái)控制誤差不一定是可取的。

    4.2 基于步長(zhǎng)規(guī)劃的優(yōu)化算法

    結(jié)合平行弦法和改變步長(zhǎng)的方法來(lái)進(jìn)行優(yōu)化,在采用平行弦確定各段TN值的前提下,根據(jù)曲線的擬合誤差分布情況來(lái)重新規(guī)劃各分段所使用的步長(zhǎng),使擬合曲線的最大誤差值低于允差的同時(shí)避免出現(xiàn)精度過(guò)高的情況并降低分段數(shù)。步長(zhǎng)規(guī)劃算法:

    設(shè)第i段的兩個(gè)基點(diǎn)為Pti、Pti+1(i=1,2,…,n)。兩點(diǎn)間的步長(zhǎng)為L(zhǎng)i測(cè)得滿足擬合條件的臨界步長(zhǎng)Lc以及能使大部分曲線段不合格的較大步長(zhǎng)Lb,Li的取值范圍為L(zhǎng)c<Li<Lb。則兩基點(diǎn)的坐標(biāo)為:

    擬合曲線時(shí),比較該段雙圓弧曲線擬合的最大法向誤差,當(dāng)最大法向誤差Δmax>Δ允時(shí),選用更小的Li重新擬合。當(dāng)最大法向誤差Δmax過(guò)小時(shí)則增大Li,這樣將使Li自動(dòng)向理想值逼近,直到雙圓弧曲線符合要求為止并繼續(xù)擬合下一段曲線。步長(zhǎng)規(guī)劃的程序框圖如圖8所示。

    圖8 步長(zhǎng)規(guī)劃程序框圖

    該程序運(yùn)行后將自動(dòng)生成各分段的步長(zhǎng)并繪出相應(yīng)的雙圓弧曲線。

    對(duì)于上例中的曲線,按照第三部分的實(shí)驗(yàn)結(jié)果,Lc取13mm,Lb取20mm。

    運(yùn)行程序生成各段步長(zhǎng)14.11,14.62,15.21…19.39,共生成10段。擬合曲線誤差如圖9所示。

    圖9 上例改進(jìn)的擬合誤差分布圖

    與圖7(a)的情況相比,分段數(shù)同為10段,采用變化步長(zhǎng)的擬合曲線沒(méi)有超出允差。最大誤差為:6.77E-06mm。誤差較之前分布更加均勻且沒(méi)有出現(xiàn)精度過(guò)高問(wèn)題。

    實(shí)驗(yàn)證明規(guī)劃步長(zhǎng)擬合法可以用較少的分段數(shù)擬合出符合允差的雙圓弧擬合曲線,計(jì)算量少,更易于實(shí)現(xiàn)工件的輪廓曲線的加工。

    5 結(jié)束語(yǔ)

    本文運(yùn)用VB實(shí)現(xiàn)優(yōu)化雙圓弧算法的計(jì)算機(jī)模型,對(duì)擬合誤差進(jìn)行了計(jì)算和分析。討論了影響擬合誤差的因素,并在此基礎(chǔ)上提出了雙圓弧擬合的改進(jìn)算法,通過(guò)步長(zhǎng)規(guī)劃方法改進(jìn)雙圓弧曲線的擬合,保證了較低的計(jì)算量、較少的步長(zhǎng)以及允許的加工精度,使擬合曲線的獲得更簡(jiǎn)化,更有利于計(jì)算機(jī)的編程和硬件實(shí)現(xiàn)。

    [1]沈紀(jì)桂.數(shù)控加工圓弧擬合的優(yōu)化方法[J].機(jī)械工程學(xué)報(bào),1997,(6): 64-69.

    [2]蘇步青,劉鼎元.計(jì)算幾何[M].上海: 上海科學(xué)出版社,1981: 195-203.

    [3]丁克會(huì),席平原,周紅斌.參數(shù)方程曲線的最優(yōu)逼近算法及實(shí)現(xiàn)[J].機(jī)械傳動(dòng),2008,(6): 57-59.

    [4]沈紀(jì)桂,方志梅,李凌豐.優(yōu)化雙圓弧擬合[J].工程圖學(xué)學(xué)報(bào),1995,(1): 53-56.

    [5]郭文蘭,張 彤.曲線方程的雙圓弧擬合算法及實(shí)現(xiàn)[J].哈爾濱理工大學(xué)學(xué)報(bào),2004,9(6): 87-89.

    [6]虞銘財(cái),楊勛年,汪國(guó)昭.整體最優(yōu)雙圓弧擬合[J].高校應(yīng)用數(shù)學(xué)學(xué)報(bào)A輯,2004,19(2): 225-232.

    猜你喜歡
    法向圓弧斜率
    落石法向恢復(fù)系數(shù)的多因素聯(lián)合影響研究
    淺析圓弧段高大模板支撐體系設(shè)計(jì)與應(yīng)用
    物理圖像斜率的變化探討
    物理之友(2020年12期)2020-07-16 05:39:16
    外圓弧面銑削刀具
    求斜率型分式的取值范圍
    低溫狀態(tài)下的材料法向發(fā)射率測(cè)量
    基于子孔徑斜率離散采樣的波前重構(gòu)
    MMC-MTDC輸電系統(tǒng)新型直流電壓斜率控制策略
    落石碰撞法向恢復(fù)系數(shù)的模型試驗(yàn)研究
    數(shù)控車削圓弧螺紋程序編制與加工
    国产精品国产三级国产专区5o| 国产男人的电影天堂91| 一区二区三区高清视频在线| 只有这里有精品99| 免费播放大片免费观看视频在线观看| 亚洲国产欧美在线一区| 偷拍熟女少妇极品色| 久久国内精品自在自线图片| 女人被狂操c到高潮| 国产在线一区二区三区精| 国产视频内射| 亚洲成人一二三区av| 小蜜桃在线观看免费完整版高清| 婷婷色麻豆天堂久久| 久久鲁丝午夜福利片| 国产有黄有色有爽视频| 黑人高潮一二区| 91精品伊人久久大香线蕉| 干丝袜人妻中文字幕| 午夜视频国产福利| 日韩一本色道免费dvd| 2018国产大陆天天弄谢| 在线 av 中文字幕| 国产激情偷乱视频一区二区| 亚洲激情五月婷婷啪啪| 久久精品综合一区二区三区| 久久精品国产亚洲网站| 十八禁网站网址无遮挡 | 丝袜美腿在线中文| 九九爱精品视频在线观看| 成人高潮视频无遮挡免费网站| 日产精品乱码卡一卡2卡三| 国产亚洲午夜精品一区二区久久 | av在线老鸭窝| 男的添女的下面高潮视频| 亚洲av福利一区| 久久久久九九精品影院| 精品久久久噜噜| 欧美日韩一区二区视频在线观看视频在线 | 内射极品少妇av片p| 国产av在哪里看| 久久久精品欧美日韩精品| 国产麻豆成人av免费视频| 国产色婷婷99| 亚洲精品亚洲一区二区| 国产综合精华液| 韩国高清视频一区二区三区| 久久鲁丝午夜福利片| 五月伊人婷婷丁香| 好男人视频免费观看在线| 99热这里只有精品一区| 一个人看的www免费观看视频| 日韩欧美 国产精品| 亚洲精品一二三| 夫妻午夜视频| 热99在线观看视频| 18禁在线无遮挡免费观看视频| 欧美极品一区二区三区四区| 99久久九九国产精品国产免费| 五月伊人婷婷丁香| 一区二区三区乱码不卡18| 伊人久久国产一区二区| 寂寞人妻少妇视频99o| 成人国产麻豆网| 亚洲欧洲国产日韩| 青春草国产在线视频| 一个人观看的视频www高清免费观看| 一级毛片久久久久久久久女| 国产精品一区二区三区四区久久| 99久久中文字幕三级久久日本| 亚洲国产精品专区欧美| 久久久成人免费电影| 3wmmmm亚洲av在线观看| av免费观看日本| 日韩欧美精品v在线| 亚洲真实伦在线观看| 国产精品无大码| 欧美精品国产亚洲| 能在线免费观看的黄片| 免费av观看视频| 亚洲欧美精品专区久久| 久久精品久久精品一区二区三区| 卡戴珊不雅视频在线播放| 国产午夜精品一二区理论片| 成人一区二区视频在线观看| 久久久久免费精品人妻一区二区| 成人鲁丝片一二三区免费| 久久久久性生活片| 人妻系列 视频| 99热这里只有精品一区| 亚洲色图av天堂| 日本wwww免费看| 久久午夜福利片| 老女人水多毛片| 欧美区成人在线视频| 久久久久久久久久黄片| 少妇人妻精品综合一区二区| 内射极品少妇av片p| 久久精品久久久久久噜噜老黄| 亚洲aⅴ乱码一区二区在线播放| 女的被弄到高潮叫床怎么办| 国产熟女欧美一区二区| 激情五月婷婷亚洲| 啦啦啦啦在线视频资源| 国产精品1区2区在线观看.| 九九久久精品国产亚洲av麻豆| 精品久久久精品久久久| 久久久欧美国产精品| 蜜臀久久99精品久久宅男| 99九九线精品视频在线观看视频| 久久精品国产鲁丝片午夜精品| 熟女电影av网| 男女边吃奶边做爰视频| 国产精品久久久久久久久免| 精品国产露脸久久av麻豆 | 春色校园在线视频观看| 国产又色又爽无遮挡免| 又大又黄又爽视频免费| 国产精品嫩草影院av在线观看| 啦啦啦韩国在线观看视频| 在线观看av片永久免费下载| 丝袜美腿在线中文| 丝袜喷水一区| 汤姆久久久久久久影院中文字幕 | 精品久久久久久久人妻蜜臀av| 日产精品乱码卡一卡2卡三| 最近中文字幕高清免费大全6| 国产久久久一区二区三区| 最近中文字幕2019免费版| 国产男人的电影天堂91| 久久草成人影院| 99久久九九国产精品国产免费| 免费电影在线观看免费观看| 婷婷色麻豆天堂久久| 成人漫画全彩无遮挡| 欧美极品一区二区三区四区| 亚洲av日韩在线播放| 久久久久久久午夜电影| 久久精品国产亚洲av天美| 一级毛片黄色毛片免费观看视频| 一级二级三级毛片免费看| 中文字幕制服av| 最近中文字幕2019免费版| 一级片'在线观看视频| 久久久久性生活片| 高清午夜精品一区二区三区| 舔av片在线| 成年av动漫网址| 精品熟女少妇av免费看| 国产综合精华液| 色尼玛亚洲综合影院| 日本-黄色视频高清免费观看| 亚洲国产欧美人成| 内地一区二区视频在线| 亚洲av国产av综合av卡| 欧美日韩国产mv在线观看视频 | 国产 亚洲一区二区三区 | 七月丁香在线播放| 高清av免费在线| 最后的刺客免费高清国语| 午夜精品在线福利| 又爽又黄无遮挡网站| 欧美成人午夜免费资源| 中文字幕人妻熟人妻熟丝袜美| 精品一区二区三卡| 特级一级黄色大片| 亚洲国产高清在线一区二区三| 日韩欧美精品v在线| 久久人人爽人人爽人人片va| 91在线精品国自产拍蜜月| 国产亚洲精品久久久com| 18禁动态无遮挡网站| 欧美97在线视频| 大香蕉久久网| 欧美成人午夜免费资源| av在线天堂中文字幕| 2021少妇久久久久久久久久久| 欧美高清成人免费视频www| av免费观看日本| 大香蕉97超碰在线| 成年版毛片免费区| 乱人视频在线观看| 日韩视频在线欧美| 亚洲欧洲国产日韩| 免费黄色在线免费观看| 大香蕉久久网| 精品久久国产蜜桃| 精品一区二区三区视频在线| 男插女下体视频免费在线播放| 国产女主播在线喷水免费视频网站 | 日韩中字成人| 三级毛片av免费| 国产黄片美女视频| 免费大片18禁| 可以在线观看毛片的网站| 日韩欧美精品v在线| 日本一二三区视频观看| 免费看光身美女| 亚洲精品乱久久久久久| 成年av动漫网址| 久久久亚洲精品成人影院| 久久韩国三级中文字幕| 伦精品一区二区三区| 亚洲综合色惰| 中文字幕亚洲精品专区| 只有这里有精品99| 能在线免费观看的黄片| 亚洲av二区三区四区| 有码 亚洲区| 内地一区二区视频在线| 亚洲自拍偷在线| 日本一二三区视频观看| 国产黄片美女视频| 日韩成人av中文字幕在线观看| 国产 一区 欧美 日韩| 汤姆久久久久久久影院中文字幕 | 国产成人午夜福利电影在线观看| av国产免费在线观看| 日韩亚洲欧美综合| 国产白丝娇喘喷水9色精品| 精品久久久久久久久av| 中文在线观看免费www的网站| 2021天堂中文幕一二区在线观| 国模一区二区三区四区视频| av在线亚洲专区| 国产精品国产三级国产av玫瑰| 建设人人有责人人尽责人人享有的 | 日本欧美国产在线视频| 少妇熟女aⅴ在线视频| 久久久国产一区二区| 中文字幕制服av| 日韩一本色道免费dvd| 欧美成人一区二区免费高清观看| 中文字幕av在线有码专区| 菩萨蛮人人尽说江南好唐韦庄| 三级毛片av免费| 国产熟女欧美一区二区| 性插视频无遮挡在线免费观看| 18+在线观看网站| 亚洲av.av天堂| 一级av片app| 日本黄大片高清| 卡戴珊不雅视频在线播放| 少妇的逼水好多| 日日撸夜夜添| 国产久久久一区二区三区| 汤姆久久久久久久影院中文字幕 | 一个人看的www免费观看视频| 精品一区在线观看国产| 亚洲在线自拍视频| 亚洲伊人久久精品综合| 午夜福利成人在线免费观看| 男人舔奶头视频| 99久久精品国产国产毛片| 欧美激情久久久久久爽电影| 欧美性猛交╳xxx乱大交人| 亚洲av成人精品一区久久| 少妇丰满av| 麻豆乱淫一区二区| 亚洲激情五月婷婷啪啪| 国产人妻一区二区三区在| 黄色日韩在线| 国产高清有码在线观看视频| 亚洲国产日韩欧美精品在线观看| 亚洲欧洲日产国产| 国产探花在线观看一区二区| 日韩大片免费观看网站| 别揉我奶头 嗯啊视频| 亚洲成人精品中文字幕电影| 日韩一区二区三区影片| 男女啪啪激烈高潮av片| 日本一二三区视频观看| 精品国产露脸久久av麻豆 | 18+在线观看网站| 国产欧美另类精品又又久久亚洲欧美| 男女那种视频在线观看| 哪个播放器可以免费观看大片| 日本一二三区视频观看| 夫妻性生交免费视频一级片| 久久6这里有精品| 国产亚洲精品av在线| 国产成人精品婷婷| 最近最新中文字幕大全电影3| 亚洲熟女精品中文字幕| 麻豆av噜噜一区二区三区| 永久免费av网站大全| 大片免费播放器 马上看| 人人妻人人看人人澡| 国产欧美另类精品又又久久亚洲欧美| 内地一区二区视频在线| 又大又黄又爽视频免费| 亚洲精品国产av蜜桃| 免费看a级黄色片| 午夜精品在线福利| 欧美xxⅹ黑人| 亚洲av免费高清在线观看| 九草在线视频观看| 免费观看性生交大片5| 赤兔流量卡办理| 在线免费观看的www视频| 非洲黑人性xxxx精品又粗又长| 亚洲电影在线观看av| 国产免费视频播放在线视频 | 97热精品久久久久久| 日日摸夜夜添夜夜添av毛片| 2018国产大陆天天弄谢| 免费人成在线观看视频色| 日韩欧美三级三区| 观看免费一级毛片| 国产精品国产三级国产av玫瑰| 内地一区二区视频在线| 亚洲美女视频黄频| 直男gayav资源| 中文欧美无线码| 18禁裸乳无遮挡免费网站照片| 久久久久久伊人网av| 日本午夜av视频| 成人漫画全彩无遮挡| 亚洲成人一二三区av| 夫妻性生交免费视频一级片| 亚洲美女视频黄频| 欧美日韩国产mv在线观看视频 | 国产熟女欧美一区二区| 国产精品av视频在线免费观看| 少妇裸体淫交视频免费看高清| 少妇人妻精品综合一区二区| 伦理电影大哥的女人| 亚洲自偷自拍三级| 在线 av 中文字幕| 可以在线观看毛片的网站| 一个人看的www免费观看视频| 国产精品福利在线免费观看| 2021少妇久久久久久久久久久| 日日摸夜夜添夜夜添av毛片| 2021天堂中文幕一二区在线观| 一级毛片我不卡| 免费看不卡的av| av网站免费在线观看视频 | 美女高潮的动态| 欧美激情国产日韩精品一区| 国产爱豆传媒在线观看| 亚洲精品一区蜜桃| 全区人妻精品视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 黄色配什么色好看| 亚洲av免费在线观看| 亚洲va在线va天堂va国产| 久久精品国产鲁丝片午夜精品| 亚洲欧美一区二区三区黑人 | 综合色av麻豆| 日韩一区二区三区影片| 国产亚洲一区二区精品| 久久久久久久久中文| 亚洲精品aⅴ在线观看| 国精品久久久久久国模美| 久久久久久久久中文| 色综合站精品国产| 日韩成人伦理影院| 国产精品女同一区二区软件| 内射极品少妇av片p| 在线观看美女被高潮喷水网站| 中国美白少妇内射xxxbb| av一本久久久久| 国产成人午夜福利电影在线观看| 成人漫画全彩无遮挡| 中国美白少妇内射xxxbb| 青春草亚洲视频在线观看| 免费av观看视频| 舔av片在线| 日韩成人伦理影院| 最近中文字幕高清免费大全6| 伊人久久精品亚洲午夜| 一个人观看的视频www高清免费观看| 欧美不卡视频在线免费观看| 大片免费播放器 马上看| av福利片在线观看| 毛片一级片免费看久久久久| 色网站视频免费| 国语对白做爰xxxⅹ性视频网站| 三级男女做爰猛烈吃奶摸视频| 91精品伊人久久大香线蕉| 国产精品1区2区在线观看.| 国产精品久久视频播放| 国产熟女欧美一区二区| 男人爽女人下面视频在线观看| 国产免费福利视频在线观看| 亚洲欧美一区二区三区国产| 国产黄频视频在线观看| 亚洲成人一二三区av| 日本午夜av视频| 中文字幕制服av| 亚洲欧洲国产日韩| 久久久a久久爽久久v久久| 亚洲18禁久久av| 亚洲在线观看片| 亚洲av免费在线观看| 午夜精品一区二区三区免费看| av又黄又爽大尺度在线免费看| 中文字幕久久专区| 一级片'在线观看视频| 天堂av国产一区二区熟女人妻| 国产精品三级大全| 久久久久久久久中文| 亚洲av电影在线观看一区二区三区 | 国产精品99久久久久久久久| 国产免费视频播放在线视频 | 欧美xxxx性猛交bbbb| 国产极品天堂在线| 亚洲精品中文字幕在线视频 | 国产乱来视频区| 亚洲成人精品中文字幕电影| 99热这里只有精品一区| 丰满人妻一区二区三区视频av| 国产高清三级在线| 男女边吃奶边做爰视频| 欧美+日韩+精品| 日日摸夜夜添夜夜爱| av福利片在线观看| 成年人午夜在线观看视频 | 午夜福利视频1000在线观看| 成人毛片60女人毛片免费| 男人爽女人下面视频在线观看| 男女边摸边吃奶| 久久精品久久久久久久性| 精品久久久久久久久亚洲| 久久久久久久久久久丰满| 亚洲av免费在线观看| 秋霞伦理黄片| 我要看日韩黄色一级片| 亚洲最大成人中文| 日韩欧美精品免费久久| 永久免费av网站大全| 成年av动漫网址| 久久精品久久久久久噜噜老黄| 亚洲精品一区蜜桃| 2018国产大陆天天弄谢| 国产伦在线观看视频一区| 日韩av在线大香蕉| 午夜福利在线观看免费完整高清在| 国产片特级美女逼逼视频| 久久久午夜欧美精品| 女人被狂操c到高潮| 国产一区二区在线观看日韩| 免费黄频网站在线观看国产| 激情五月婷婷亚洲| 亚洲在久久综合| 日韩av在线大香蕉| 亚洲av免费在线观看| 亚洲三级黄色毛片| 国产综合精华液| 网址你懂的国产日韩在线| 一区二区三区乱码不卡18| 边亲边吃奶的免费视频| 狂野欧美白嫩少妇大欣赏| 只有这里有精品99| 久久精品国产亚洲av天美| 777米奇影视久久| 18+在线观看网站| 久久精品久久精品一区二区三区| 一级毛片黄色毛片免费观看视频| 在线免费观看不下载黄p国产| 国产精品福利在线免费观看| 亚洲国产精品国产精品| 国产一区二区亚洲精品在线观看| 国产一区二区三区综合在线观看 | 一区二区三区四区激情视频| 网址你懂的国产日韩在线| 春色校园在线视频观看| 亚洲av日韩在线播放| 免费黄色在线免费观看| 黄片wwwwww| 极品少妇高潮喷水抽搐| 色哟哟·www| 亚洲性久久影院| 毛片一级片免费看久久久久| 国内精品美女久久久久久| 欧美xxxx性猛交bbbb| 街头女战士在线观看网站| 久久久久久久国产电影| 爱豆传媒免费全集在线观看| 国产午夜精品一二区理论片| 国产av国产精品国产| 免费av不卡在线播放| 日韩亚洲欧美综合| 亚洲真实伦在线观看| 午夜福利成人在线免费观看| 美女黄网站色视频| 少妇的逼水好多| 日日干狠狠操夜夜爽| 大陆偷拍与自拍| 五月天丁香电影| 久久精品人妻少妇| 国产精品蜜桃在线观看| 日韩欧美三级三区| 久久久久久久久中文| 久久97久久精品| 夫妻性生交免费视频一级片| 在现免费观看毛片| 自拍偷自拍亚洲精品老妇| 亚洲人成网站在线播| xxx大片免费视频| kizo精华| 国产精品久久久久久av不卡| 边亲边吃奶的免费视频| 成人二区视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产伦精品一区二区三区四那| 最近中文字幕高清免费大全6| 欧美三级亚洲精品| 天堂俺去俺来也www色官网 | 亚洲久久久久久中文字幕| 天堂√8在线中文| 搞女人的毛片| 午夜福利视频1000在线观看| 黄片无遮挡物在线观看| 国产一区有黄有色的免费视频 | av专区在线播放| 两个人视频免费观看高清| 身体一侧抽搐| 久久久久久久大尺度免费视频| 777米奇影视久久| 久久精品国产亚洲av涩爱| 国产亚洲精品久久久com| 黄片无遮挡物在线观看| 国产成人a区在线观看| 在线观看美女被高潮喷水网站| 一级毛片黄色毛片免费观看视频| 成年免费大片在线观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲,欧美,日韩| 2021少妇久久久久久久久久久| 亚洲成人精品中文字幕电影| 国产一区有黄有色的免费视频 | 成人午夜精彩视频在线观看| 亚洲欧美日韩卡通动漫| 我的老师免费观看完整版| 国产单亲对白刺激| 亚洲一级一片aⅴ在线观看| 国产伦一二天堂av在线观看| 久久久久九九精品影院| 国产真实伦视频高清在线观看| 十八禁国产超污无遮挡网站| 国产亚洲91精品色在线| 久久这里有精品视频免费| 国产精品无大码| 校园人妻丝袜中文字幕| 国产成人一区二区在线| 综合色丁香网| 少妇高潮的动态图| 欧美成人a在线观看| 大陆偷拍与自拍| 在线 av 中文字幕| 日日啪夜夜撸| 校园人妻丝袜中文字幕| 国产成人freesex在线| 69av精品久久久久久| 国产av码专区亚洲av| 国产黄色视频一区二区在线观看| 国产精品国产三级专区第一集| 中文字幕制服av| 国产成人91sexporn| 国产乱人偷精品视频| 亚洲18禁久久av| 久久精品国产亚洲av涩爱| 男人爽女人下面视频在线观看| 日韩精品有码人妻一区| 18+在线观看网站| 午夜免费激情av| 哪个播放器可以免费观看大片| 免费人成在线观看视频色| 男女视频在线观看网站免费| 日产精品乱码卡一卡2卡三| 十八禁国产超污无遮挡网站| 啦啦啦中文免费视频观看日本| 国产精品综合久久久久久久免费| 欧美成人午夜免费资源| 亚洲av电影在线观看一区二区三区 | 久久精品国产亚洲av涩爱| 男人爽女人下面视频在线观看| 日韩不卡一区二区三区视频在线| 日韩在线高清观看一区二区三区| 精品久久久精品久久久| 久久国产乱子免费精品| 久久这里有精品视频免费| av国产久精品久网站免费入址| 免费高清在线观看视频在线观看| 国产亚洲av片在线观看秒播厂 | 久久久精品免费免费高清| 丰满少妇做爰视频| 中文欧美无线码| 成人一区二区视频在线观看| 午夜福利在线观看吧| 人人妻人人澡欧美一区二区| 精品久久久噜噜| 毛片一级片免费看久久久久| 国产三级在线视频| 久久精品熟女亚洲av麻豆精品 | 精品少妇黑人巨大在线播放| 国产老妇伦熟女老妇高清| 九九久久精品国产亚洲av麻豆| 国产伦在线观看视频一区| 亚洲国产精品成人久久小说| 亚洲国产精品专区欧美| 国产探花极品一区二区| 国产麻豆成人av免费视频| 美女黄网站色视频| 精品酒店卫生间| 建设人人有责人人尽责人人享有的 | 18+在线观看网站| 日本欧美国产在线视频| 草草在线视频免费看| 欧美激情国产日韩精品一区| 国产精品美女特级片免费视频播放器| 免费电影在线观看免费观看| 精品酒店卫生间|