• <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ù)控車削圓弧螺紋程序編制與加工
    亚洲少妇的诱惑av| 久久精品国产亚洲av高清一级| 中文字幕av电影在线播放| svipshipincom国产片| aaaaa片日本免费| 麻豆国产av国片精品| 久久精品91无色码中文字幕| 午夜福利免费观看在线| 久久精品国产清高在天天线| 涩涩av久久男人的天堂| av天堂在线播放| 欧美一区二区精品小视频在线| 亚洲五月色婷婷综合| 99久久99久久久精品蜜桃| 级片在线观看| 可以免费在线观看a视频的电影网站| 亚洲国产欧美一区二区综合| 电影成人av| 亚洲国产欧美一区二区综合| 大陆偷拍与自拍| 免费女性裸体啪啪无遮挡网站| 亚洲五月婷婷丁香| 午夜久久久久精精品| 亚洲av五月六月丁香网| 午夜福利欧美成人| 一个人观看的视频www高清免费观看 | 久久久精品国产亚洲av高清涩受| 国产日韩一区二区三区精品不卡| 亚洲欧美日韩另类电影网站| 宅男免费午夜| 一边摸一边抽搐一进一小说| www.精华液| 精品久久久精品久久久| 波多野结衣av一区二区av| 久久久久久久午夜电影| 99久久精品国产亚洲精品| 怎么达到女性高潮| 色哟哟哟哟哟哟| 国产精品免费一区二区三区在线| 久久午夜亚洲精品久久| 我的亚洲天堂| 夜夜看夜夜爽夜夜摸| 国产免费av片在线观看野外av| 国产成人精品久久二区二区91| 国产成人欧美在线观看| 在线十欧美十亚洲十日本专区| 嫩草影视91久久| 变态另类丝袜制服| 午夜老司机福利片| 好看av亚洲va欧美ⅴa在| 黑人巨大精品欧美一区二区蜜桃| 老汉色av国产亚洲站长工具| 99在线视频只有这里精品首页| 精品第一国产精品| 亚洲电影在线观看av| 久久天堂一区二区三区四区| 国产色视频综合| av天堂久久9| 久热爱精品视频在线9| 亚洲片人在线观看| a在线观看视频网站| 色尼玛亚洲综合影院| 免费无遮挡裸体视频| 国产一区在线观看成人免费| 日本a在线网址| 中文字幕色久视频| 777久久人妻少妇嫩草av网站| 日韩精品免费视频一区二区三区| 亚洲黑人精品在线| 日本在线视频免费播放| 精品一区二区三区视频在线观看免费| 黄色视频不卡| 97超级碰碰碰精品色视频在线观看| 国产成人精品久久二区二区免费| 日韩成人在线观看一区二区三区| 免费女性裸体啪啪无遮挡网站| 久久久久国产一级毛片高清牌| 成人永久免费在线观看视频| 日韩精品免费视频一区二区三区| 国产色视频综合| 一级毛片女人18水好多| 亚洲伊人色综图| 免费高清在线观看日韩| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩瑟瑟在线播放| 欧美午夜高清在线| 男人的好看免费观看在线视频 | 在线av久久热| 99国产精品一区二区三区| 国产在线观看jvid| 免费看十八禁软件| 国产成人av激情在线播放| 久久久久久免费高清国产稀缺| 国产成年人精品一区二区| 一区在线观看完整版| 免费在线观看日本一区| av中文乱码字幕在线| 叶爱在线成人免费视频播放| 欧美中文综合在线视频| 久久久久久免费高清国产稀缺| 欧美日本中文国产一区发布| 18禁黄网站禁片午夜丰满| 男男h啪啪无遮挡| 亚洲精品在线观看二区| 久久久久久久精品吃奶| 丰满人妻熟妇乱又伦精品不卡| 亚洲三区欧美一区| 十八禁网站免费在线| 久久久国产成人免费| 一级毛片精品| 欧美日本亚洲视频在线播放| 一边摸一边抽搐一进一小说| 精品高清国产在线一区| 美女午夜性视频免费| cao死你这个sao货| 淫妇啪啪啪对白视频| 国产成人欧美| 亚洲七黄色美女视频| 欧美另类亚洲清纯唯美| 亚洲五月色婷婷综合| 国产精品一区二区三区四区久久 | 精品久久蜜臀av无| 中文亚洲av片在线观看爽| 99精品欧美一区二区三区四区| 麻豆国产av国片精品| 中文字幕高清在线视频| 男人舔女人下体高潮全视频| 视频在线观看一区二区三区| 黄色a级毛片大全视频| 国产高清视频在线播放一区| 中文字幕另类日韩欧美亚洲嫩草| 九色国产91popny在线| 麻豆久久精品国产亚洲av| 亚洲中文av在线| 极品人妻少妇av视频| 老司机深夜福利视频在线观看| 一a级毛片在线观看| 午夜福利在线观看吧| 久久久久久久久免费视频了| 最好的美女福利视频网| 91国产中文字幕| 日日爽夜夜爽网站| 欧美在线黄色| 国产av又大| 一区二区三区国产精品乱码| 日本 av在线| 18禁裸乳无遮挡免费网站照片 | 老司机深夜福利视频在线观看| 两个人免费观看高清视频| 禁无遮挡网站| 香蕉久久夜色| 亚洲成av片中文字幕在线观看| 国产欧美日韩一区二区精品| 欧美性长视频在线观看| 久久精品国产亚洲av香蕉五月| 久久性视频一级片| 国产成+人综合+亚洲专区| 久久精品aⅴ一区二区三区四区| 两性夫妻黄色片| 一区二区三区国产精品乱码| 国产精品久久久久久精品电影 | 窝窝影院91人妻| 精品久久久久久,| 视频在线观看一区二区三区| 91精品国产国语对白视频| 久久久国产欧美日韩av| 亚洲成av人片免费观看| 又黄又粗又硬又大视频| 天堂√8在线中文| e午夜精品久久久久久久| 曰老女人黄片| 少妇被粗大的猛进出69影院| 桃色一区二区三区在线观看| 色哟哟哟哟哟哟| 日韩欧美免费精品| 亚洲av电影不卡..在线观看| 日日爽夜夜爽网站| 久久天躁狠狠躁夜夜2o2o| 90打野战视频偷拍视频| 51午夜福利影视在线观看| 99riav亚洲国产免费| 国产极品粉嫩免费观看在线| 两人在一起打扑克的视频| 一区二区三区国产精品乱码| 国产av一区二区精品久久| 久9热在线精品视频| av电影中文网址| 黄色女人牲交| av福利片在线| 一二三四在线观看免费中文在| 国产一区二区三区综合在线观看| 熟女少妇亚洲综合色aaa.| 亚洲精品国产区一区二| 久久亚洲真实| 成人永久免费在线观看视频| 国产一区二区在线av高清观看| 色综合欧美亚洲国产小说| 亚洲 欧美 日韩 在线 免费| 久久精品国产亚洲av高清一级| 久久人人爽av亚洲精品天堂| 午夜福利欧美成人| 精品无人区乱码1区二区| 国产精品国产高清国产av| 国产成年人精品一区二区| 首页视频小说图片口味搜索| 每晚都被弄得嗷嗷叫到高潮| 国产精品 国内视频| 成人欧美大片| 精品人妻在线不人妻| 一级作爱视频免费观看| 亚洲专区国产一区二区| 香蕉丝袜av| 久久午夜综合久久蜜桃| 一级a爱片免费观看的视频| 男人操女人黄网站| 又紧又爽又黄一区二区| 90打野战视频偷拍视频| 日韩欧美国产在线观看| 波多野结衣av一区二区av| 国产亚洲av嫩草精品影院| 精品国产国语对白av| 一区二区日韩欧美中文字幕| 国产区一区二久久| 12—13女人毛片做爰片一| 日韩欧美一区视频在线观看| 一区二区三区高清视频在线| 色综合欧美亚洲国产小说| 九色国产91popny在线| 激情在线观看视频在线高清| 成年版毛片免费区| 色在线成人网| 久久中文看片网| 搡老妇女老女人老熟妇| 亚洲成av人片免费观看| 法律面前人人平等表现在哪些方面| 欧美国产精品va在线观看不卡| 无遮挡黄片免费观看| 亚洲av电影在线进入| 很黄的视频免费| 国产亚洲精品综合一区在线观看 | 香蕉国产在线看| 夜夜躁狠狠躁天天躁| 国产一区二区三区视频了| 美女高潮到喷水免费观看| 18禁黄网站禁片午夜丰满| 精品卡一卡二卡四卡免费| 成人国产综合亚洲| 国产精品九九99| 精品午夜福利视频在线观看一区| 黑人巨大精品欧美一区二区蜜桃| 亚洲五月婷婷丁香| 又紧又爽又黄一区二区| 制服丝袜大香蕉在线| 黄片播放在线免费| 国产高清激情床上av| 人人妻人人澡人人看| 欧美成人午夜精品| 亚洲黑人精品在线| www日本在线高清视频| 国产黄a三级三级三级人| 美女扒开内裤让男人捅视频| 男人舔女人的私密视频| 欧美乱妇无乱码| 不卡av一区二区三区| 亚洲精品美女久久av网站| 精品久久久久久成人av| 村上凉子中文字幕在线| 如日韩欧美国产精品一区二区三区| 在线十欧美十亚洲十日本专区| 琪琪午夜伦伦电影理论片6080| 欧美大码av| 天天躁狠狠躁夜夜躁狠狠躁| 非洲黑人性xxxx精品又粗又长| 在线观看免费视频日本深夜| 69av精品久久久久久| 精品国产国语对白av| 天堂动漫精品| 欧美大码av| 无遮挡黄片免费观看| 亚洲七黄色美女视频| av有码第一页| 精品国产一区二区三区四区第35| 淫秽高清视频在线观看| 黄色片一级片一级黄色片| 日日夜夜操网爽| 一级黄色大片毛片| 亚洲第一青青草原| 亚洲精品中文字幕一二三四区| 99精品久久久久人妻精品| 中出人妻视频一区二区| 久久青草综合色| 精品电影一区二区在线| 国产精品自产拍在线观看55亚洲| 亚洲熟妇熟女久久| 欧美色视频一区免费| 亚洲人成电影观看| 精品久久蜜臀av无| 欧美最黄视频在线播放免费| 国产视频一区二区在线看| 亚洲欧美日韩无卡精品| 亚洲人成77777在线视频| 久9热在线精品视频| 久久影院123| 国产99久久九九免费精品| 成人特级黄色片久久久久久久| www.999成人在线观看| 成年人黄色毛片网站| 制服人妻中文乱码| 国产一区二区三区视频了| 国产欧美日韩一区二区三| 日韩高清综合在线| 欧美最黄视频在线播放免费| 又紧又爽又黄一区二区| 国产xxxxx性猛交| 69av精品久久久久久| 久久久久精品国产欧美久久久| 一本大道久久a久久精品| 午夜影院日韩av| 亚洲欧美一区二区三区黑人| 亚洲成a人片在线一区二区| 精品高清国产在线一区| 男女午夜视频在线观看| 两个人视频免费观看高清| 国产精品久久电影中文字幕| 一级a爱视频在线免费观看| 久久精品影院6| 中文字幕色久视频| 亚洲中文字幕日韩| 国产成人精品在线电影| 亚洲国产中文字幕在线视频| 亚洲精品中文字幕在线视频| 黑丝袜美女国产一区| 成人三级做爰电影| 日韩中文字幕欧美一区二区| 视频区欧美日本亚洲| 国产精品久久电影中文字幕| 好男人电影高清在线观看| 最近最新免费中文字幕在线| 啦啦啦韩国在线观看视频| 亚洲精品国产区一区二| 丝袜人妻中文字幕| 黑人巨大精品欧美一区二区mp4| 精品免费久久久久久久清纯| 欧美精品啪啪一区二区三区| 亚洲七黄色美女视频| 熟妇人妻久久中文字幕3abv| 99国产综合亚洲精品| 亚洲国产欧美网| 色综合站精品国产| 波多野结衣巨乳人妻| 变态另类丝袜制服| 一区二区三区国产精品乱码| 免费高清在线观看日韩| 一本综合久久免费| 少妇裸体淫交视频免费看高清 | 99国产精品一区二区蜜桃av| 久久精品亚洲熟妇少妇任你| 69av精品久久久久久| 欧美日韩瑟瑟在线播放| 夜夜爽天天搞| 777久久人妻少妇嫩草av网站| cao死你这个sao货| 又紧又爽又黄一区二区| 久久性视频一级片| 悠悠久久av| 香蕉国产在线看| 国产午夜福利久久久久久| 欧美成人一区二区免费高清观看 | 波多野结衣一区麻豆| 日韩精品中文字幕看吧| 国产91精品成人一区二区三区| 日韩视频一区二区在线观看| 欧美色视频一区免费| 男男h啪啪无遮挡| 欧美日本亚洲视频在线播放| 人成视频在线观看免费观看| 久久影院123| av天堂在线播放| 亚洲精品久久国产高清桃花| 久久精品国产综合久久久| 精品国产美女av久久久久小说| 性少妇av在线| 变态另类成人亚洲欧美熟女 | 精品久久久久久久人妻蜜臀av | 黄片大片在线免费观看| 亚洲自偷自拍图片 自拍| 国产午夜精品久久久久久| 国产片内射在线| 波多野结衣av一区二区av| 麻豆国产av国片精品| 亚洲专区国产一区二区| 性色av乱码一区二区三区2| 国产成+人综合+亚洲专区| 18禁美女被吸乳视频| 超碰成人久久| 国产精品久久久久久人妻精品电影| 女性被躁到高潮视频| 免费搜索国产男女视频| 亚洲中文字幕日韩| 欧美日韩福利视频一区二区| 韩国精品一区二区三区| 精品一区二区三区四区五区乱码| 欧美av亚洲av综合av国产av| 多毛熟女@视频| 国产91精品成人一区二区三区| 9色porny在线观看| 变态另类丝袜制服| 91在线观看av| 国产成人精品无人区| 777久久人妻少妇嫩草av网站| 中国美女看黄片| 动漫黄色视频在线观看| 中文字幕人妻熟女乱码| 女人高潮潮喷娇喘18禁视频| 中出人妻视频一区二区| 在线永久观看黄色视频| 日韩大尺度精品在线看网址 | 岛国在线观看网站| 99riav亚洲国产免费| 妹子高潮喷水视频| 国产麻豆69| 天天躁狠狠躁夜夜躁狠狠躁| 日本在线视频免费播放| 十八禁网站免费在线| 麻豆一二三区av精品| 大陆偷拍与自拍| 日韩国内少妇激情av| 母亲3免费完整高清在线观看| 国产视频一区二区在线看| 制服诱惑二区| 午夜免费成人在线视频| 婷婷丁香在线五月| 国产精品免费一区二区三区在线| 久久 成人 亚洲| 国产av精品麻豆| 岛国在线观看网站| 午夜福利视频1000在线观看 | 视频在线观看一区二区三区| 一进一出好大好爽视频| 精品一区二区三区视频在线观看免费| 色综合亚洲欧美另类图片| 免费在线观看亚洲国产| 麻豆久久精品国产亚洲av| 丝袜美腿诱惑在线| 国产在线精品亚洲第一网站| avwww免费| 精品国产一区二区久久| 在线观看www视频免费| 久久亚洲真实| 亚洲最大成人中文| 极品教师在线免费播放| 国产激情欧美一区二区| 国产91精品成人一区二区三区| 欧美另类亚洲清纯唯美| 国产私拍福利视频在线观看| 成人精品一区二区免费| 亚洲av第一区精品v没综合| 丝袜美足系列| 99国产精品免费福利视频| 国产av一区在线观看免费| 国产精品二区激情视频| 亚洲avbb在线观看| 亚洲精品在线观看二区| 久久久久九九精品影院| 国产免费男女视频| av有码第一页| 久久久久久久久免费视频了| 亚洲国产毛片av蜜桃av| 人人妻人人澡人人看| 成人av一区二区三区在线看| 这个男人来自地球电影免费观看| 在线十欧美十亚洲十日本专区| 免费看美女性在线毛片视频| 亚洲av成人一区二区三| 亚洲无线在线观看| 日韩大尺度精品在线看网址 | 97超级碰碰碰精品色视频在线观看| 91字幕亚洲| 欧美日韩福利视频一区二区| 欧美亚洲日本最大视频资源| 丝袜人妻中文字幕| 侵犯人妻中文字幕一二三四区| 亚洲国产欧美一区二区综合| 久久久久国产精品人妻aⅴ院| 久久精品国产亚洲av高清一级| 国产亚洲精品av在线| 女生性感内裤真人,穿戴方法视频| 正在播放国产对白刺激| 国产av一区二区精品久久| 成熟少妇高潮喷水视频| 嫩草影视91久久| 如日韩欧美国产精品一区二区三区| 黄色成人免费大全| 午夜福利免费观看在线| 我的亚洲天堂| 在线观看一区二区三区| 三级毛片av免费| 成人三级黄色视频| 人人澡人人妻人| av在线天堂中文字幕| 黄色成人免费大全| 午夜老司机福利片| 男女做爰动态图高潮gif福利片 | 国产精品久久久久久精品电影 | 91字幕亚洲| 午夜福利成人在线免费观看| 长腿黑丝高跟| 少妇熟女aⅴ在线视频| 久久九九热精品免费| 国产精品久久久久久亚洲av鲁大| 国产精品久久电影中文字幕| 久久香蕉国产精品| 激情在线观看视频在线高清| 久9热在线精品视频| 一级a爱片免费观看的视频| 亚洲情色 制服丝袜| 亚洲狠狠婷婷综合久久图片| 精品不卡国产一区二区三区| aaaaa片日本免费| 国产免费av片在线观看野外av| 看片在线看免费视频| 两个人视频免费观看高清| 亚洲精品国产区一区二| 中出人妻视频一区二区| 亚洲精品美女久久久久99蜜臀| 香蕉久久夜色| 国产色视频综合| 日韩成人在线观看一区二区三区| 亚洲在线自拍视频| 欧美日韩亚洲综合一区二区三区_| 亚洲精品一卡2卡三卡4卡5卡| 中文字幕人妻熟女乱码| 在线观看免费视频网站a站| 中文字幕av电影在线播放| 69av精品久久久久久| 纯流量卡能插随身wifi吗| 亚洲精品国产一区二区精华液| 亚洲男人天堂网一区| 无限看片的www在线观看| 国产亚洲av嫩草精品影院| 久久国产乱子伦精品免费另类| 亚洲情色 制服丝袜| 午夜亚洲福利在线播放| 精品国产亚洲在线| 大码成人一级视频| www.999成人在线观看| 在线十欧美十亚洲十日本专区| 国产欧美日韩一区二区精品| 在线观看一区二区三区| 麻豆一二三区av精品| 久久天躁狠狠躁夜夜2o2o| 久久久久亚洲av毛片大全| svipshipincom国产片| 一个人免费在线观看的高清视频| 亚洲欧美日韩另类电影网站| 桃色一区二区三区在线观看| 人人妻人人爽人人添夜夜欢视频| 欧美乱色亚洲激情| 久久九九热精品免费| 国产又爽黄色视频| 免费一级毛片在线播放高清视频 | 18禁观看日本| 成人国语在线视频| 欧美黑人精品巨大| 精品一区二区三区av网在线观看| 国产黄a三级三级三级人| 在线国产一区二区在线| 国产精品爽爽va在线观看网站 | 亚洲精品国产一区二区精华液| 国产av一区二区精品久久| 欧美大码av| 变态另类丝袜制服| 国产亚洲精品av在线| 亚洲一码二码三码区别大吗| 一级黄色大片毛片| 制服诱惑二区| 此物有八面人人有两片| 成年女人毛片免费观看观看9| 欧美在线黄色| 日韩欧美国产在线观看| 久久久水蜜桃国产精品网| 在线观看舔阴道视频| 亚洲aⅴ乱码一区二区在线播放 | 欧美中文日本在线观看视频| avwww免费| 精品一区二区三区av网在线观看| 亚洲av日韩精品久久久久久密| netflix在线观看网站| 亚洲av第一区精品v没综合| 侵犯人妻中文字幕一二三四区| 真人一进一出gif抽搐免费| 亚洲精品在线观看二区| 自拍欧美九色日韩亚洲蝌蚪91| 在线永久观看黄色视频| 欧美最黄视频在线播放免费| 久久久水蜜桃国产精品网| 欧美日本视频| 最近最新中文字幕大全免费视频| 亚洲专区国产一区二区| 久久久国产成人免费| 日韩精品青青久久久久久| 亚洲一区高清亚洲精品| 精品一品国产午夜福利视频| 久久国产精品人妻蜜桃| 精品免费久久久久久久清纯| 亚洲 国产 在线| 欧美久久黑人一区二区| 国产精品久久久久久人妻精品电影| 99国产极品粉嫩在线观看| 9色porny在线观看| 亚洲精品中文字幕一二三四区| 欧美+亚洲+日韩+国产| 久久久久久久久免费视频了| 一个人观看的视频www高清免费观看 | 亚洲成av片中文字幕在线观看|