局部最優小線段路徑解析光順的實時前瞻插補方法及系統的制作方法

            文檔序號:10686340閱讀:350來源:國知局
            局部最優小線段路徑解析光順的實時前瞻插補方法及系統的制作方法
            【專利摘要】本發明提供了一種局部最優小線段路徑解析光順的實時前瞻插補方法及系統,包括:路徑光順模塊、速度規劃模塊和實時插補模塊。其中,所述路徑光順模塊,通過讀取數控G代碼G01段信息,提取代碼中的小線段的坐標信息,每相鄰四個坐標信息作為一個光順單元,依次對每個光順單元實現最優的轉接光順;所述速度規劃模塊又稱為躍度有限的S型速度規劃;所述實時插補模塊,用于根據速度規劃模塊得到速度生成離散插補點。本發明中的方法和系統計算效率高,計算量小,編程實現簡單,可廣泛適用于高速高精的數控機床。
            【專利說明】
            局部最優小線段路徑解析光順的實時前瞻插補方法及系統
            技術領域
            [0001] 本發明涉及數控加工技術領域的在線系統,具體地,涉及一種局部最優小線段路 徑解析光順的實時前瞻插補方法及系統。
            【背景技術】
            [0002] 在CNC加工中,采用一系列線段近似表示自由曲線(Bezier曲線、Bspline曲線等)。 為了獲得自由曲線的精確近似,線段都非常短。通常情況下,將這種采用短線段近似表示的 刀具路徑稱為小線段刀具路徑。一方面,由于小線段路徑僅位置連續,需要頻繁加減速,該 速度的大小及變化對零件加工精度、表面質量和加工效率都有著重要影響;另一方面,加工 精度越高,生成的坐標數據越多,需要存儲空間越大。因此,針對大型且高精加工零件的局 部最優小線段路徑解析光順的實時前瞻插補系統在工業應用中十分重要。
            [0003] 目前常用的光順方法雖然考慮了光順精度要求和曲率連續性且利于速度規劃,但 是光順對象通常為所有小線段路徑,對計算機硬件系統提出了很高的要求,為此,通常采用 滑動的前瞻窗口實現光順,窗口寬度同樣依賴于計算機的硬件系統。具體的實現過程如下:
            [0004] 根據窗口寬度依次讀入數據,窗口數據讀入完畢,對窗口內的小線段路徑實施光 順算法,光順對象為窗口內的全部小線段路徑,光順后,從窗口的出口處輸出一個速度規劃 單元,在窗口的入口處讀入一個新的坐標信息,繼續重復上述過程,直至數據輸出完畢。但 是上述方法僅限于數控加工領域的離線系統。
            [0005] 對現有文獻的技術檢索發現,雖然速度規劃中的雙向掃描算法被證明是一種非常 有效的方法,但窗口內的每個速度規劃單元都需經過反向躍度有限加速和正向躍度有限加 速過程,數據的計算量很大且不適于在線系統。

            【發明內容】

            [0006] 針對現有技術中的缺陷,本發明的目的是提供一種局部最優小線段路徑解析光順 的實時前瞻插補方法及系統。
            [0007] 根據本發明提供的局部最優小線段路徑解析光順的實時前瞻插補方法,包括如下 步驟:
            [0008] 路徑光順步驟:在前瞻窗口入口處讀入小線段刀具路徑的坐標信息,對坐標信息 依次進行路徑光順,獲得速度規劃單元并保存在前瞻窗口中;
            [0009] 速度規劃步驟:獲得速度規劃單元的弧長,并進行反向躍度有限加速和正向躍度 有限加速,獲得速度規劃單元銜接點處的最終速度,執行躍度有限的S型規劃獲得相應的運 動學曲線;
            [0010] 實時插補步驟:由于光順后的刀具路徑是由線段和參數曲線組成的混合路徑,因 此采用混合插補方式生成離散插補點,并將該插補點從前瞻窗口的出口處輸出。
            [0011]優選地,所述路徑光順步驟包括:
            [0012]步驟A1:從前瞻窗口的入口處依次讀入小線段刀具路徑的坐標信息,坐標信息的 個數取決于前瞻窗口的寬度,所述寬度是指前瞻窗口被坐標信息填滿時所能夠容納的字符 數;依次選取在前瞻窗口的入口處讀入的毗鄰的四個坐標信息,作為一個光順單元;
            [0013] 步驟A2:每個光順單元能夠生成對應的兩條光順曲線,并獲得對應光順曲線的曲 率極大值;
            [0014] 步驟A3:以光順單元生成的兩條光順曲線的曲率極大值和達到最小為優化目標, 以小線段刀具路徑的長度和設定的光順誤差為約束條件,建立曲率優化模型;
            [0015] 步驟A4:將步驟A3中的優化目標轉化為關于控制多邊形邊長的一元二次方程,在 光順誤差約束條件下,所述一元二次方程的解唯一,最終,解析地獲得當前光順單元的兩條 光順曲線;
            [0016]步驟A5:在步驟4中,得到兩條光順曲線,即具有五個控制點的三次Bspline曲線, 對第一條光順曲線執行節點u = 0.5的重復插入操作,其中,u表示三次B樣條曲線的曲線參 數,亦可稱為節點,其取值范圍為[0,1 ],將Bspline曲線分割為兩條對稱的貝齊爾曲線,將 當前光順單元中前兩個線性刀路坐標信息及相應的Bezier保存在前瞻窗口;
            [0017]步驟A6:將在步驟A4中獲得的兩條光順曲線的交點作為兩條光順曲線的最優銜接 占.
            [0018] 步驟A7:判斷小線段刀具路徑的坐標信息在前瞻窗口的入口處是否全部光順完 畢;
            [0019] 若否,則將步驟A6中的最優銜接點將作為下一個光順單元的起始坐標信息,并取 最優銜接點后面依次毗鄰的未被光順的三個線性刀路的坐標信息,組成新的光順單元,并 則執行步驟A2至步驟A6;若是,則執行速度規劃步驟。
            [0020] 優選地,所述速度規劃步驟包括:
            [0021] 步驟B1:根據辛普森公式依次計算前瞻窗口內速度規劃單元的弧長;并根據最大 進給速度、加速度、躍度以及弓高誤差進行反向躍度有限加速和正向躍度有限加速,獲得速 度規劃單元銜接點處的最終速度;
            [0022]步驟B2:由步驟B1中獲得的弧長和銜接點處的速度,針對每個速度規劃單元執行 躍度有限的S型速度規劃,得到運動學曲線,并對速度規劃單元的速度曲線執行實時插補步 驟;其中,所述運動學曲線包括:速度、加速度、躍度曲線;
            [0023]步驟B3:速度規劃單元生成插補點后,前瞻窗口輸出第一個速度規劃單元對應的 插補點信息,即入口處的坐標信息向前瞻窗口的出口方向移動一個坐標信息的寬度;
            [0024]步驟B4:判斷整個線性刀路是否光順完畢,若否,在前瞻窗口的入口端讀入一個新 的坐標信息,并與前一個光順單元得到的最優銜接點以及該銜接點后面依次毗鄰兩個坐標 信息組成新的光順單元,執行步驟A4至A7;若是,則結束。
            [0025] 優選地,光順曲線采用傳統的具有五個控制點的三次B樣條,光順曲線的曲率極大 值能夠根據控制多邊形長度和轉角解析獲得。
            [0026] 優選地,步驟A1中的線段刀具坐標信息取自數控G代碼G01段,前瞻窗口的寬度依 賴于計算機硬件。
            [0027] 根據本發明提供的局部最優小線段路徑解析光順的實時前瞻插補系統,包括:
            [0028] 路徑光順模塊:用于在前瞻窗口入口處讀入小線段刀具路徑的坐標信息,對坐標 信息依次進行路徑光順,獲得速度規劃單元并保存在前瞻窗口中;
            [0029] 速度規劃模塊:用于獲得速度規劃單元的弧長,并進行反向躍度有限加速和正向 躍度有限加速,獲得速度規劃單元銜接點處的最終速度,執行躍度有限的S型規劃獲得相應 的運動學曲線;
            [0030] 實時插補模塊:通過混合插補方式生成離散插補點,并將該插補點從前瞻窗口的 出口處輸出。
            [0031] 優選地,所述路徑光順模塊包括:
            [0032] 坐標信息讀入子模塊,用于從前瞻窗口的入口處依次讀入小線段刀具路徑的坐標 信息,坐標信息的個數取決于前瞻窗口的寬度,所述寬度是指前瞻窗口被坐標信息填滿時 所能夠容納的字符數;依次選取在前瞻窗口的入口處讀入的毗鄰的四個坐標信息,作為一 個光順單元;
            [0033] 光順曲線生成子模塊,用于根據每個光順單元生成對應的兩條光順曲線,并獲得 對應光順曲線的曲率極大值;
            [0034] 曲率優化模型構建子模塊:用于以光順單元生成的兩條光順曲線的曲率極大值和 達到最小為優化目標,以小線段刀具路徑的長度和設定的光順誤差為約束條件,建立曲率 優化模型;
            [0035] 多邊形頂點坐標獲取子模塊,用于將優化目標轉化為關于控制多邊形邊長的一元 二次方程,在光順誤差約束條件下,所述一元二次方程的解唯一,最終,解析地獲得當前光 順單元的兩條光順曲線;
            [0036] 貝齊爾曲線生成子模塊,用于當得到兩條光順曲線,即具有五個控制點的三次 Bspl ine曲線,對第一條光順曲線執行節點u = 0.5的重復插入操作,其中,u表示三次B樣條 曲線的曲線參數,亦可稱為節點,其取值范圍為[0,1],將Bspline曲線分割為兩條對稱的貝 齊爾曲線,此時,將當前光順單元中前兩個線性刀路坐標信息及相應的Bezier保存在前瞻 窗口;
            [0037] 最優銜接點獲取子模塊:用于將獲得的兩條光順曲線的交點作為兩條光順曲線的 最優銜接點;
            [0038] 判斷子模塊,用于判斷小線段刀具路徑的坐標信息在前瞻窗口的入口處是否全部 光順完畢;若否,則將最優銜接點將作為下一個光順單元的起始坐標信息,并取最優銜接點 后面依次毗鄰的未被光順的三個線性刀路的坐標信息,組成新的轉接單元。
            [0039]優選地,所述速度規劃模塊包括:
            [0040] 銜接點處的最終速度獲取子模塊:用于根據辛普森公式依次計算前瞻窗口內速度 規劃單元的弧長;并根據最大進給速度、加速度、躍度以及弓高誤差進行反向躍度有限加速 和正向躍度有限加速,獲得速度規劃單元銜接點處的最終速度;
            [0041] 運動學曲線生成子模塊:用于根據的弧長和銜接點處的速度,針對每個速度規劃 單元執行躍度有限的S型速度規劃,得到運動學曲線,并對速度規劃單元的速度曲線執行實 時插補步驟;其中,所述運動學曲線包括:速度、加速度、躍度曲線;
            [0042] 插補點信息輸出子模塊:用于在速度規劃單元生成插補點后,前瞻窗口輸出第一 個速度規劃單元對應的插補點信息,即入口處的坐標信息向前瞻窗口的出口方向移動一個 坐標信息的寬度;
            [0043]新坐標錄入子模塊:用于判斷整個線性刀路是否光順完畢,若否,在前瞻窗口的入 口端讀入一個新的坐標信息,并與前一個光順單元得到的最優銜接點以及該銜接點后面依 次毗鄰兩個坐標信息組成新的光順單元;若是,則結束。
            [0044] 優選地,小線段刀具坐標信息取自數控G代碼G01段,前瞻窗口的寬度依賴于計算 機硬件。
            [0045] 與現有技術相比,本發明具有如下的有益效果:
            [0046] 1、本發明提供局部最優小線段路徑解析光順的實時前瞻插補方法計算效率高,計 算量小,編程實現簡單,尤其適用于高速高精的數控機床。
            [0047] 2、本發明提供局部最優小線段路徑解析光順的實時前瞻插補方法能夠加工大型 復雜零件,且加工精度高。
            【附圖說明】
            [0048]通過閱讀參照以下附圖對非限制性實施例所作的詳細描述,本發明的其它特征、 目的和優點將會變得更明顯:
            [0049] 圖1為本發明提供的局部最優小線段路徑解析光順實時前瞻插補方法的流程示意 圖;
            [0050] 圖2為局部最優的解析光順算法示意圖;
            [0051] 圖3(a)至圖3(f)為速度規劃單元中刀具路徑的六種類型的示意圖;
            [0052]圖4(a)、圖4(b)為躍度有限加速模式的兩種類型的示意圖;
            [0053]圖5(a)至圖5(c)為混合插補中線性插補和參數插補的轉換示意圖;
            [0054]圖6為二維金魚小線段刀具路徑示意圖;
            [0055] 圖7為光順前和光順后速度和加速度曲線對比圖;
            [0056] 圖8為光順前和光順后X軸和Y軸速度和加速度曲線對比圖;
            [0057] 圖9為采用混合插補方式得到的光滑插補軌跡示意圖。
            【具體實施方式】
            [0058] 下面結合具體實施例對本發明進行詳細說明。以下實施例將有助于本領域的技術 人員進一步理解本發明,但不以任何形式限制本發明。應當指出的是,對本領域的普通技術 人員來說,在不脫離本發明構思的前提下,還可以做出若干變化和改進。這些都屬于本發明 的保護范圍。
            [0059] 根據本發明提供的最優小線段路徑解析光順實時前瞻插補方法,包括如下步驟:
            [0060] 路徑光順步驟:在前瞻窗口入口處讀入小線段刀具路徑的坐標信息,對坐標信息 依次進行路徑光順,獲得速度規劃單元并保存在前瞻窗口中;
            [0061] 速度規劃步驟:獲得速度規劃單元的弧長,并進行反向躍度有限加速和正向躍度 有限加速,獲得速度規劃單元銜接點處的最終速度,執行躍度有限的S型規劃獲得相應的運 動學曲線;
            [0062] 實時插補步驟:由于光順后的刀具路徑是由線段和參數曲線組成的混合路徑,因 此采用混合插補方式生成離散插補點,并將該插補點從前瞻窗口的出口處輸出。
            [0063]所述路徑光順步驟包括:
            [0064]步驟A1:從前瞻窗口的入口處依次讀入小線段刀具路徑的坐標信息,坐標信息的 個數取決于前瞻窗口的寬度,所述寬度是指前瞻窗口被坐標信息填滿時所能夠容納的字符 數;依次選取在前瞻窗口的入口處讀入的毗鄰的四個坐標信息,作為一個光順單元;
            [0065] 步驟A2:每個光順單元能夠生成對應的兩條光順曲線,并獲得對應光順曲線的曲 率極大值;
            [0066] 具體地,光順曲線采用傳統的具有五個控制點的三次B樣條,光順曲線的曲率極大 值能夠根據控制多邊形長度和轉角解析獲得;五個控制點的三次B樣條作為光順曲線,因其 曲率極大值可解析計算且計算量小;執行節點插補操作,方便確定速度規劃單元。
            [0067] 步驟A3:以光順單元生成的兩條光順曲線的曲率極大值和達到最小為優化目標, 以小線段刀具路徑的長度和設定的光順誤差為約束條件,建立曲率優化模型;
            [0068] 步驟A4:將步驟A3中的優化目標轉化為關于控制多邊形邊長的一元二次方程,在 光順誤差約束條件下,所述一元二次方程的解唯一,即能夠通過解析獲得控制多邊形的頂 點坐標;
            [0069]步驟A5:在步驟4中,得到兩條光順曲線,即具有五個控制點的三次Bspline曲線。 對第一條光順曲線執行節點u = 0.5的重復插入操作,其中,u表示三次B樣條曲線的曲線參 數,亦可稱為節點,其取值范圍為[0,1 ],將其分割為兩條對稱的貝齊爾曲線,此時,將當前 光順單元中前兩個線性刀路坐標信息及相應的Bezier保存在前瞻窗口。
            [0070] 步驟A6:將在步驟A4中獲得的兩條光順曲線的交點作為兩條光順曲線的最優銜接 占.
            [0071] 步驟A7:判斷小線段刀具路徑的坐標信息在前瞻窗口的入口處是否全部光順完 畢;
            [0072] 若否,則將步驟A6中的最優銜接點將作為下一個光順單元的起始坐標信息,并取 最優銜接點后面依次毗鄰的未被光順的三個線性刀路的坐標信息,組成新的光順單元,并 則執行步驟A2至步驟A6;若是,則執行速度規劃步驟。
            [0073]所述速度規劃步驟包括:
            [0074]步驟B1:根據辛普森公式依次計算前瞻窗口內速度規劃單元的弧長;并根據最大 進給速度、加速度、躍度以及弓高誤差進行反向躍度有限加速和正向躍度有限加速,獲得速 度規劃單元銜接點處的最終速度;
            [0075]步驟B2:由步驟B2中獲得的弧長和銜接點處的速度,針對每個速度規劃單元執行 躍度有限的S型速度規劃,得到運動學曲線,并對速度規劃單元的速度曲線執行實時插補步 驟;其中,所述運動學曲線包括:速度、加速度、躍度曲線;
            [0076] 步驟B3:速度規劃單元生成插補點后,前瞻窗口輸出第一個速度規劃單元對應的 插補點信息,即入口處的坐標信息向前瞻窗口的出口方向移動一個坐標信息的寬度;
            [0077] 步驟B4:判斷前瞻窗口內的速度規劃單元是否全部規劃完畢,若否,則選取下一個 速度規劃單元執行步驟B1至步驟B3;若是,則結束。
            [0078] 步驟A1中的小線段刀具坐標信息取自數控G代碼G01段,前瞻窗口的寬度依賴于計 算機硬件。
            [0079]所述最優銜接點是局部最優的,該最優銜接點能夠保證每個小線段刀具路徑轉角 僅被光順兩次,能夠適用于在線系統。
            [0080]所述躍度有限的S型速度規劃,即指生成的速度曲線躍度有限。
            [0081] 在實時插補步驟中,根據規劃速度,進行插補,將插補所得信息直接用于位置的實 時控制。
            [0082] 具體地,如圖1所示:首先設定實時前瞻滑動窗口,窗口的寬度依賴于計算機的硬 件;其次,從窗口入口端依次讀入小線段刀具路徑的坐標信息,毗鄰的四個坐標信息點作為 一個光順單元,每個光順單元生成一個局部最優的銜接點,作為下一個光順單元的起始坐 標信息,從窗口的出口端至入口端依次進行光順;然后,對窗口內的光順曲線進行節點插入 操作,將一條B樣條曲線分割成兩條對稱的貝奇爾曲線,即生成可用于速度規劃的速度規劃 單元,從窗口入口端至出口端進行反向躍度有限加速和正向躍度有限加速,獲得速度規劃 單元銜接點處的最終速度,該最終速度可以保證機床的整成運轉,在此基礎上,進行最終的 躍度有限的S型速度規劃;最后,采用線性插補和參數插補的混合插補模式,生成離散插補 點。
            [0083] 具體描述如下:
            [0084]步驟S1:從前瞻窗口入口處依次讀入小線段刀具路徑的坐標信息,坐標信息的個 數依賴于窗口的寬度,直至坐標數據將窗口填滿。
            [0085] 步驟1)中小線段刀具的定義亦可參考文獻[Zhang L B,You Y P,He J,et al .The transition algorithm based on parametric spline curve for high-speed machining of continuous short line segments.The International Journal of Advanced Manufacturing Technology 2011;52(1-4):245-254],[Wei Fan,Chen-Han Lee,Ji-Hong Chen.A real-time curvature-smooth interpolation scheme and motion planning for CNC machining of short line segments.International Journal of Machine Tools and Manufacture 2015;96:27-46].
            [0086] 小線段刀具路徑的坐標信息可表示為集合形式 0, . . .,m},其中Pi表示第i個坐標信息,三個分量xi,yi,zi分別表示刀具參考點在工件直角 坐標系中的對應X軸、Y軸、Z軸的坐標;設前瞻窗口的寬度為m,那么窗口可讀入的坐標信息 也為m〇
            [0087] 步驟S2:在前瞻窗口的出口端選取毗鄰的四個坐標點作為光順單元,相應地,每個 光順單元應生成兩條光順曲線。
            [0088]如圖2所示,首先選取四個點?〇,?1,內和?3構成第一個速度規劃單元,其中 £為設定 的光順誤差,Q/為第i個光順單元光順后得到的局部最優銜接點。
            [0089]步驟S3:光順曲線采用傳統的具有五個控制點的三次B樣條,曲率極大值可根據控 制多邊形長度和轉角解析獲得。
            [0090]曲率極大值Kmax為
            [0092] 其中,山為光順曲線控制多邊形邊長,a表示矢量PoPjPP^夾角的一半。
            [0093] 步驟3)中公式的推導過程可參考文獻[Huan Zhao,YaoAn Lu,LiMin Zhu,Han Ding.Look-ahead interpolation of short line segments using B-spline curve fitting of dominant points.Proc IMechE Part B:Journal of Engineering Manufacture 2014;1-13].
            [0094]步驟S4:在步驟S3中得到的曲率極大值的基礎上,以光順單元內光順曲線的曲率 極大值和最小為目標,小線段刀具路徑的長度和設定的光順誤差為約束,建立曲率優化模 型。
            [0095]局部最優光順的優化模型為
            [0097] Li^| iQi-i'Pil
            [0098] Li彡 | |PiPi+11
            [0099] Li^;2e(1+c)/sinai [0100] s.t.Li+i^| |PiPi+i|
            [0101] Li+i^| |Pi+iPi+2|
            [0102] Li+i^;2e (l+c)/sinai+i
            [0103] Li+Li+i^ | |PiPi+i|
            [0104]式中:|卜| |表示歐式范數;c表示任意的一個正常數;n表示線性刀路的坐標信息 額個數; 表示第i個光順單元第一條光順曲線的曲率極大值;表示第i個光順單元第 二條光順曲線的曲率極大值;U表示第i個光順單元第一條光順曲線的轉接長度的最大值; L1+1表示第i個光順單元第二條光順曲線的轉接長度的最大值;e表示指定的最大光順誤差; ai表不矢量Qi-i'Pi和PiPi+i夾角的一半。
            [0105]步驟S5:將步驟S4得到的優化目標轉化為關于控制多邊形邊長的一元二次方程, 在步驟S4的約束條件下,該方程的解唯一,即可解析地獲得控制多邊形的頂點坐標。
            [0106]直接求解優化問題是很困難的,若h(l+c)辛A1+1,可將目標函數表示為關于控制 多邊形長度的一元二次方程:
            [0108] 其中,人i和人i+i表示引入的輔助變量,人i = 4sinai/(3cos2ai),Ai+1 = 4(l+c)sinai+1/ (3cos2ai+i) 〇
            [0109] 光順單元內的兩條光順曲線的銜接點一定在中間兩個坐標點的連線上,即 0<4<丨|伙+1丨卜
            [0110]其解為
            [0112] 上式僅保證目標函數取最小值,并未考慮小線段路徑的長度和設定的光順誤差, 優化問題的約束模型確定的可行區間可表示為兩點所在的閉區間,如圖2(a)的[Q^Qs]和圖 2(c)的[Qi,Qi+i]所示。
            [0113] 最終,得到局部最優的銜接點
            [0115] 根據該銜接點,可計算控制多邊形的邊長及控制頂點坐標。
            [0116] 步驟S6:通過節點重復插入操作,將步驟S5得到的第一條光順B樣條分割為兩條對 稱的貝奇爾曲線,并將第一、二個坐標信息及相應的貝奇爾曲線作為速度規劃單元保存在 前瞻窗口。
            [0117] 如圖3所示,通過節點重復插入操作,生成用于速度規劃的速度規劃單元,每個速 度規劃單元由線段和貝奇爾曲線組成,共由6種類型。
            [0118] 步驟S7:在步驟S5確定控制多邊形頂點坐標的同時,也確定了兩條光順曲線的最 優銜接點,該銜接點將作為下一個光順單元的起始坐標信息,再選后面毗鄰的三個坐標信 息,組成新的光順單元并進行光順,依次類推,直至將前瞻窗口內的小線段刀具路徑光順完 畢。
            [0119]步驟S8:根據辛普森公式依次計算前瞻窗口內每個速度規劃單元的弧長;在此基 礎上,根據最大進給速度、加速度、躍度和弓高誤差進行反向躍度有限加速和正向躍度有限 加速,獲得速度規劃單元銜接點處的最終速度;
            [0120]具體的辛普森方法計算速度規劃單元弧長方法請參考文獻[Lei WT,Sung MP,Lin LY,Huang JJ.Fast real-time NURBS path interpolation for CNC machine tools.International Journal of Machine Tools and Manufacture 2007:47(10)1530-41].
            [0121] 如圖4所示,躍度有限加速分為兩種情況。首先按如圖4(a)所示的加速模式加速, 獲得加速后速度的同時,還獲得本速度規劃單元能夠達到的最大切向加速度;其次,判斷這 個最大加速度是否超過切向加速度約束,若沒超過,此速度最為躍度有限加速獲得的最終 速度;反之,采用如圖4(b)所示的加速模式重新加速,此時獲得的速度方可作為躍度有限加 速獲得的最終速度。
            [0122] 對于反向躍度有限加速,前瞻窗口入口處的速度規劃單元末速度始終設置為零; 首先從末速度為零的速度規劃單元開始反向加速,加速后獲得的速度與最大進給速度、以 及法向加速度、法向躍度、弓高誤差確定的速度進行比較,取最小值作為該速度規劃單元的 起始速度;然后再以這個速度作為前一個速度規劃單元的末速度,為了將該算法應用于在 線系統,并不是每個速度規劃單元都進行反向加速,以新分配末速度的速度規劃單元為例, 最大進給速度以及根據法向加速度、法向躍度、弓高誤差計算得到的速度的最小值若小于 等于末速度,則這個最小值作為該速度規劃單元根據反向掃描確定的起始速度;反之,需要 執行反向加速,確定最終的起始速度。
            [0123] 由于前瞻窗口執行先進先出,故對于正向躍度有限加速,掃描方向為前瞻窗口出 口至入口方向。在反向掃描的基礎上進行正向掃描,故在反向掃描中考慮的最大進給速度、 法向加速度、法向躍度以及弓高誤差在這里不予考慮,只要考慮切向加速度和切向躍度即 可。注意,在正向掃描中,僅于整個刀具路徑的第一個速度規劃單元的起始速度設置為零。
            [0124] 步驟S9:由步驟S8中獲得的弧長和銜接點處的速度,針對每個速度規劃單元執行 躍度有限的S型速度規劃,得到運動學曲線,并將第一個速度規劃單元的速度曲線傳遞給插 補模塊。
            [0125] 具體的速度規劃算法請參考文獻[Xu Du,Jie Huang,Li_Min Zhu.A complete S-shape feed rate scheduling approach for NURBS interpolator.Journal of Computational Design and Engineering 2015:206-217].
            [0126] 步驟S10:光順后的刀具路徑是由線段和參數曲線組成的混合路徑,采用混合插補 方式生成離散插補點,并將該插補數據從前瞻出口輸出。
            [0127] 如圖3所示,共有6種類型的刀具路徑,每種都是由線段和貝奇爾曲線組成,線性插 補對于線段軌跡有效,參數插補對于參數段有效,故采用混合插補模式。
            [0128] 對于線段路徑,采用線性插補
            [0130]其中,Ts表不插補周期;C((k+1)TS)表不第k+1個插補點;A Sj表不第j個插補周期 對應的弧長增量;Pi表示線段路徑的起始坐標信息;I為方向矢量。
            [0131]對于Bezier曲線,采用參數插補
            [0133] 其中,C(u)為參數曲線,Uk表示當前插補點對應的曲線參數,Uk+i表示后一個插補 點對應的曲線參數,(:'(111〇、(:"(1 11〇分別表示為曲線在一階、二階導數,81{表示第1^個插補 周期對應的弧長增量,〈?>表示兩個矢量的內積,| | ? | |表示歐式范數。
            [0134] 混合插補的難點在于插補模式的有效切換,如圖5所示。
            [0135] 切換方式共分為三種:
            [0136] 1.參數插補轉換為線性插補
            [0137] 如圖5(a)所示,轉換條件為
            [0138] 卜' (3) lut+l >1
            [0139] 第i+1個插補點坐標信息為
            [0140] C((k+l)Ts)=b3+AskI (4)
            [0141] 其中Aa = VkTs-| |b3_C(kTs)| |,b3表示光順曲線的最后一個控制點。
            [0142] 2.線性插補轉換為線性插補
            [0143] 如圖5(b)所示,轉換條件為
            [0145] 式中:bo表示光順曲線的第一個控制點,vk表示當前的進給速度,C(kTs)表示當前 的插補點;
            [0146] 第i+1個插補點對應的曲線參數為
            [0148] 其中,弧長增量 ASk = VkTs-| |C(kTs)_bo|
            [0149] 3.參數插補轉換為參數插補
            [0150]如圖5(c)所示,轉換條件可根據公式(3)確定,第i + 1個插補點對應的曲線參數根 據公式(6)計算可得。
            [0151]以上對本發明的具體實施例進行了描述。需要理解的是,本發明并不局限于上述 特定實施方式,本領域技術人員可以在權利要求的范圍內做出各種變化或修改,這并不影 響本發明的實質內容。在不沖突的情況下,本申請的實施例和實施例中的特征可以任意相 互組合。
            【主權項】
            1. 一種局部最優小線段路徑解析光順的實時前瞻插補方法,其特征在于,包括如下步 驟: 路徑光順步驟:在前瞻窗口入口處讀入小線段刀具路徑的坐標信息,對坐標信息依次 進行路徑光順,獲得速度規劃單元并保存在前瞻窗口中; 速度規劃步驟:獲得速度規劃單元的弧長,并進行反向躍度有限加速和正向躍度有限 加速,獲得速度規劃單元銜接點處的最終速度,執行躍度有限的S型規劃獲得相應的運動學 曲線; 實時插補步驟:由于光順后的刀具路徑是由線段和參數曲線組成的混合路徑,因此采 用混合插補方式生成離散插補點,并將該插補點從前瞻窗口的出口處輸出。2. 根據權利要求1所述的局部最優小線段路徑解析光順的實時前瞻插補方法,其特征 在于,所述路徑光順步驟包括: 步驟A1:從前瞻窗口的入口處依次讀入小線段刀具路徑的坐標信息,坐標信息的個數 取決于前瞻窗口的寬度,所述寬度是指前瞻窗口被坐標信息填滿時所能夠容納的字符數; 依次選取在前瞻窗口的入口處讀入的毗鄰的四個坐標信息,作為一個光順單元; 步驟A2 :每個光順單元能夠生成對應的兩條光順曲線,并獲得對應光順曲線的曲率極 大值; 步驟A3:以光順單元生成的兩條光順曲線的曲率極大值和達到最小為優化目標,以小 線段刀具路徑的長度和設定的光順誤差為約束條件,建立曲率優化模型; 步驟A4:將步驟A3中的優化目標轉化為關于控制多邊形邊長的一元二次方程,在光順 誤差約束條件下,所述一元二次方程的解唯一,最終,解析地獲得當前光順單元的兩條光順 曲線; 步驟A5:在步驟4中,得到兩條光順曲線,即具有五個控制點的三次Bspline曲線,對第 一條光順曲線執行節點u = 0.5的重復插入操作,其中,u表示三次B樣條曲線的曲線參數,亦 可稱為節點,其取值范圍為[〇,1],將Bspline曲線分割為兩條對稱的貝齊爾曲線,將當前光 順單元中前兩個線性刀路坐標信息及相應的Bezier保存在前瞻窗口; 步驟A6:將在步驟A4中獲得的兩條光順曲線的交點作為兩條光順曲線的最優銜接點; 步驟A7:判斷小線段刀具路徑的坐標信息在前瞻窗口的入口處是否全部光順完畢; 若否,則將步驟A6中的最優銜接點將作為下一個光順單元的起始坐標信息,并取最優 銜接點后面依次毗鄰的未被光順的三個線性刀路的坐標信息,組成新的光順單元,并則執 行步驟A2至步驟A6;若是,則執行速度規劃步驟。3. 根據權利要求2所述的局部最優小線段路徑解析光順的實時前瞻插補方法,其特征 在于,所述速度規劃步驟包括: 步驟B1:根據辛普森公式依次計算前瞻窗口內速度規劃單元的弧長;并根據最大進給 速度、加速度、躍度以及弓高誤差進行反向躍度有限加速和正向躍度有限加速,獲得速度規 劃單元銜接點處的最終速度; 步驟B2:由步驟B1中獲得的弧長和銜接點處的速度,針對每個速度規劃單元執行躍度 有限的S型速度規劃,得到運動學曲線,并對速度規劃單元的速度曲線執行實時插補步驟; 其中,所述運動學曲線包括:速度、加速度、躍度曲線; 步驟B3:速度規劃單元生成插補點后,前瞻窗口輸出第一個速度規劃單元對應的插補 點信息,即入口處的坐標信息向前瞻窗口的出口方向移動一個坐標信息的寬度; 步驟M:判斷整個線性刀路是否光順完畢,若否,在前瞻窗口的入口端讀入一個新的坐 標信息,并與前一個光順單元得到的最優銜接點以及該銜接點后面依次毗鄰兩個坐標信息 組成新的光順單元,執行步驟A4至A7;若是,則結束。4. 根據權利要求2所述的局部最優小線段路徑解析光順的實時前瞻插補方法,其特征 在于,光順曲線采用傳統的具有五個控制點的三次B樣條,光順曲線的曲率極大值能夠根據 控制多邊形長度和轉角解析獲得。5. 根據權利要求2所述的局部最優小線段路徑解析光順的實時前瞻插補方法,其特征 在于,步驟A1中的線段刀具坐標信息取自數控G代碼G01段,前瞻窗口的寬度依賴于計算機 硬件。6. -種局部最優小線段路徑解析光順的實時前瞻插補系統,其特征在于,包括: 路徑光順模塊:用于在前瞻窗口入口處讀入小線段刀具路徑的坐標信息,對坐標信息 依次進行路徑光順,獲得速度規劃單元并保存在前瞻窗口中; 速度規劃模塊:用于獲得速度規劃單元的弧長,并進行反向躍度有限加速和正向躍度 有限加速,獲得速度規劃單元銜接點處的最終速度,執行躍度有限的S型規劃獲得相應的運 動學曲線; 實時插補模塊:通過混合插補方式生成離散插補點,并將該插補點從前瞻窗口的出口 處輸出。7. 根據權利要求6所述的最優小線段路徑解析光順的實時前瞻插補系統,其特征在于, 所述路徑光順模塊包括: 坐標信息讀入子模塊,用于從前瞻窗口的入口處依次讀入小線段刀具路徑的坐標信 息,坐標信息的個數取決于前瞻窗口的寬度,所述寬度是指前瞻窗口被坐標信息填滿時所 能夠容納的字符數;依次選取在前瞻窗口的入口處讀入的毗鄰的四個坐標信息,作為一個 光順單元; 光順曲線生成子模塊,用于根據每個光順單元生成對應的兩條光順曲線,并獲得對應 光順曲線的曲率極大值; 曲率優化模型構建子模塊:用于以光順單元生成的兩條光順曲線的曲率極大值和達到 最小為優化目標,以小線段刀具路徑的長度和設定的光順誤差為約束條件,建立曲率優化 豐旲型; 多邊形頂點坐標獲取子模塊,用于將優化目標轉化為關于控制多邊形邊長的一元二次 方程,在光順誤差約束條件下,所述一元二次方程的解唯一,最終,解析地獲得當前光順單 元的兩條光順曲線; 貝齊爾曲線生成子模塊,用于當得到兩條光順曲線,即具有五個控制點的三次Bspline 曲線,對第一條光順曲線執行節點u = 0.5的重復插入操作,其中,u表示三次B樣條曲線的曲 線參數,亦可稱為節點,其取值范圍為[〇,1],將Bspline曲線分割為兩條對稱的貝齊爾曲 線,此時,將當前光順單元中前兩個線性刀路坐標信息及相應的Bezier保存在前瞻窗口; 最優銜接點獲取子模塊:用于將獲得的兩條光順曲線的交點作為兩條光順曲線的最優 銜接點; 判斷子模塊,用于判斷小線段刀具路徑的坐標信息在前瞻窗口的入口處是否全部光順 完畢;若否,則將最優銜接點將作為下一個光順單元的起始坐標信息,并取最優銜接點后面 依次毗鄰的未被光順的三個線性刀路的坐標信息,組成新的轉接單元。8. 根據權利要求6所述的最優小線段路徑解析光順的實時前瞻插補系統,其特征在于, 所述速度規劃模塊包括: 銜接點處的最終速度獲取子模塊:用于根據辛普森公式依次計算前瞻窗口內速度規劃 單元的弧長;并根據最大進給速度、加速度、躍度以及弓高誤差進行反向躍度有限加速和正 向躍度有限加速,獲得速度規劃單元銜接點處的最終速度; 運動學曲線生成子模塊:用于根據的弧長和銜接點處的速度,針對每個速度規劃單元 執行躍度有限的S型速度規劃,得到運動學曲線,并對速度規劃單元的速度曲線執行實時插 補步驟;其中,所述運動學曲線包括:速度、加速度、躍度曲線; 插補點信息輸出子模塊:用于在速度規劃單元生成插補點后,前瞻窗口輸出第一個速 度規劃單元對應的插補點信息,即入口處的坐標信息向前瞻窗口的出口方向移動一個坐標 信息的寬度; 新坐標錄入子模塊:用于判斷整個線性刀路是否光順完畢,若否,在前瞻窗口的入口端 讀入一個新的坐標信息,并與前一個光順單元得到的最優銜接點以及該銜接點后面依次毗 鄰兩個坐標信息組成新的光順單元;若是,則結束。9. 根據權利要求7所述的最優小線段路徑解析光順的實時前瞻插補系統,其特征在于, 小線段刀具坐標信息取自數控G代碼G01段,前瞻窗口的寬度依賴于計算機硬件。
            【文檔編號】G05B19/414GK106054817SQ201610404868
            【公開日】2016年10月26日
            【申請日】2016年6月8日
            【發明人】杜煦, 黃杰, 朱利民
            【申請人】上海交通大學
            網友詢問留言 已有0條留言
            • 還沒有人留言評論。精彩留言會獲得點贊!
            1
            婷婷六月激情在线综合激情,亚洲国产大片,久久中文字幕综合婷婷,精品久久久久久中文字幕,亚洲一区二区三区高清不卡,99国产精品热久久久久久夜夜嗨 ,欧美日韩亚洲综合在线一区二区,99国产精品电影,伊人精品线视天天综合,精品伊人久久久大香线蕉欧美
            亚洲精品1区 国产成人一级 91精品国产欧美一区二区 亚洲精品乱码久久久久久下载 国产精品久久久久久久伊一 九色国产 国产精品九九视频 伊人久久成人爱综合网 欧美日韩亚洲区久久综合 欧美日本一道免费一区三区 夜夜爽一区二区三区精品 欧美日韩高清一区二区三区 国产成人av在线 国产精品对白交换绿帽视频 国产视频亚洲 国产在线欧美精品 国产精品综合网 国产日韩精品欧美一区色 国产日韩精品欧美一区喷 欧美日韩在线观看区一二 国产区精品 欧美视频日韩视频 中文字幕天天躁日日躁狠狠躁97 视频一二三区 欧美高清在线精品一区二区不卡 国产精品揄拍一区二区久久 99久久综合狠狠综合久久aⅴ 亚洲乱码视频在线观看 日韩在线第二页 亚洲精品无码专区在线播放 成人亚洲网站www在线观看 欧美三级一区二区 99久久精品免费看国产高清 91麻豆国产在线观看 最新日韩欧美不卡一二三区 成人在线观看不卡 日韩国产在线 在线亚洲精品 亚洲午夜久久久久中文字幕 国产精品成人久久久久久久 精品国产一区二区在线观看 欧美精品国产一区二区三区 中文在线播放 亚洲第一页在线视频 国产午夜精品福利久久 九色国产 精品国产九九 国产永久视频 久久精品人人做人人综合试看 国产一区二区三区免费观看 亚洲精品国产电影 9999热视频 国产精品资源在线 麻豆久久婷婷国产综合五月 国产精品免费一级在线观看 亚洲国产一区二区三区青草影视 中文在线播放 国产成人综合在线 国产在线观看色 国产亚洲三级 国产片一区二区三区 久久99精品久久久久久牛牛影视 亚洲欧美日韩国产 四虎永久免费网站 国产一毛片 国产精品视频在 九九热在线精品 99精品福利视频 色婷婷色99国产综合精品 97成人精品视频在线播放 精品久久久久久中文字幕 亚洲欧美一区二区三区孕妇 亚洲欧美成人网 日韩高清在线二区 国产尤物在线观看 在线不卡一区二区 91网站在线看 韩国精品福利一区二区 欧美日韩国产成人精品 99热精品久久 国产精品免费视频一区 高清视频一区 精品九九久久 欧美日韩在线观看免费 91欧美激情一区二区三区成人 99福利视频 亚洲国产精品91 久热国产在线 精品久久久久久中文字幕女 国产精品久久久久久久久99热 成人自拍视频网 国产精品视频久久久久久 久久影院国产 国产玖玖在线观看 99精品在线免费 亚洲欧美一区二区三区导航 久久久久久久综合 国产欧美日韩精品高清二区综合区 国产精品视频自拍 亚洲一级片免费 久久久久久九九 国产欧美自拍视频 视频一区二区在线观看 欧美日韩一区二区三区久久 中文在线亚洲 伊人热人久久中文字幕 日韩欧美亚洲国产一区二区三区 欧美亚洲国产成人高清在线 欧美日韩国产码高清综合人成 国产性大片免费播放网站 亚洲午夜综合网 91精品久久一区二区三区 国产无套在线播放 国产精品视频网站 国产成人亚洲精品老王 91在线网站 国产视频97 欧美黑人欧美精品刺激 国产一区二区三区免费在线视频 久久久国产精品免费看 99re6久精品国产首页 久久精品91 国产成人一级 国产成人精品曰本亚洲 日本福利在线观看 伊人成综合网 久久综合一本 国产综合久久久久久 久久精品成人免费看 久久福利 91精品国产91久久久久久麻豆 亚洲精品成人在线 亚洲伊人久久精品 欧美日本二区 国产永久视频 国产一区二 一区二区福利 国产一毛片 亚洲精品1区 毛片一区二区三区 伊人久久大香线蕉综合影 国产欧美在线观看一区 亚洲国产欧洲综合997久久 国产一区二区免费视频 国产91精品对白露脸全集观看 久久亚洲国产伦理 欧美成人伊人久久综合网 亚洲性久久久影院 久久99国产精一区二区三区! 91精品国产欧美一区二区 欧美日韩亚洲区久久综合 日韩精品一二三区 久久久夜色精品国产噜噜 国产在线精品福利91香蕉 久久久久久久亚洲精品 97se色综合一区二区二区 91国语精品自产拍在线观看性色 91久久国产综合精品女同我 日韩中文字幕a 国产成人亚洲日本精品 久久国产精品-国产精品 久久国产经典视频 久久国产精品伦理 亚洲第一页在线视频 国产精品久久久久三级 日韩毛片网 久久免费高清视频 麻豆国产在线观看一区二区 91麻豆国产福利在线观看 国产成人精品男人的天堂538 一区二区三区中文字幕 免费在线视频一区 欧美日韩国产成人精品 国产综合网站 国产资源免费观看 亚洲精品亚洲人成在线播放 精品久久久久久中文字幕专区 亚洲人成人毛片无遮挡 国产一起色一起爱 国产香蕉精品视频在 九九热免费观看 日韩亚洲欧美一区 九九热精品在线观看 精品久久久久久中文字幕专区 亚洲欧美自拍偷拍 国产精品每日更新 久久久久国产一级毛片高清板 久久天天躁狠狠躁夜夜中文字幕 久久精品片 日韩在线毛片 国产成人精品本亚洲 国产成人精品一区二区三区 九九热在线观看 国产r级在线观看 国产欧美日韩精品高清二区综合区 韩国电影一区二区 国产精品毛片va一区二区三区 五月婷婷伊人网 久久一区二区三区免费 一本色道久久综合狠狠躁篇 亚洲综合色站 国产尤物在线观看 亚洲一区亚洲二区 免费在线视频一区 欧洲精品视频在线观看 日韩中文字幕a 中文字幕日本在线mv视频精品 91精品在线免费视频 精品国产免费人成在线观看 精品a级片 中文字幕日本在线mv视频精品 日韩在线精品视频 婷婷丁香色 91精品国产高清久久久久 国产成人精品日本亚洲直接 五月综合视频 欧美日韩在线亚洲国产人 精液呈暗黄色 亚洲乱码一区 久久精品中文字幕不卡一二区 亚洲天堂精品在线 激情婷婷综合 国产免费久久精品久久久 国产精品亚洲二区在线 久久免费播放视频 五月婷婷丁香综合 在线亚洲欧美日韩 久久免费精品高清麻豆 精品久久久久久中文字幕 亚洲一区网站 国产精品福利社 日韩中文字幕免费 亚洲综合丝袜 91精品在线播放 国产精品18 亚洲日日夜夜 伊人久久大香线蕉综合影 亚洲精品中文字幕乱码影院 亚洲一区二区黄色 亚洲第一页在线视频 一区二区在线观看视频 国产成人福利精品视频 亚洲高清二区 国内成人免费视频 精品亚洲性xxx久久久 国产精品合集一区二区三区 97av免费视频 国产一起色一起爱 国产区久久 国产资源免费观看 99精品视频免费 国产成人一级 国产精品九九免费视频 欧美91精品久久久久网免费 99热国产免费 久久精品色 98精品国产综合久久 久久精品播放 中文字幕视频免费 国产欧美日韩一区二区三区在线 精品久久蜜桃 国产小视频精品 一本色道久久综合狠狠躁篇 91在线免费观看 亚洲精品区 伊人成综合网 伊人热人久久中文字幕 伊人黄色片 99国产精品热久久久久久夜夜嗨 久久免费精品视频 亚洲一区二区三区高清不卡 久久久久国产一级毛片高清板 国产片一区二区三区 久久狠狠干 99久久婷婷国产综合精品电影 国产99区 国产精品成人久久久久 久久狠狠干 青青国产在线观看 亚洲高清国产拍精品影院 国产精品一区二区av 九九热在线免费视频 伊人久久国产 国产精品久久久久久久久久一区 在线观看免费视频一区 国产精品自在在线午夜区app 国产精品综合色区在线观看 国产毛片久久久久久国产毛片 97国产免费全部免费观看 国产精品每日更新 国产尤物视频在线 九九视频这里只有精品99 一本一道久久a久久精品综合 久久综合给会久久狠狠狠 国产成人精品男人的天堂538 欧美一区二区高清 毛片一区二区三区 国产欧美日韩在线观看一区二区三区 在线国产二区 欧美不卡网 91在线精品中文字幕 在线国产福利 国内精品91久久久久 91亚洲福利 日韩欧美国产中文字幕 91久久精品国产性色也91久久 亚洲性久久久影院 欧美精品1区 国产热re99久久6国产精品 九九热免费观看 国产精品欧美日韩 久久久久国产一级毛片高清板 久久国产经典视频 日韩欧美亚洲国产一区二区三区 欧美亚洲综合另类在线观看 国产精品自在在线午夜区app 97中文字幕在线观看 视频一二三区 精品国产一区在线观看 国产欧美日韩在线一区二区不卡 欧美一区二三区 伊人成人在线观看 国内精品91久久久久 97在线亚洲 国产在线不卡一区 久久久全免费全集一级全黄片 国产精品v欧美精品∨日韩 亚洲毛片网站 在线不卡一区二区 99re热在线视频 久久激情网 国产毛片一区二区三区精品 久久亚洲综合色 中文字幕视频免费 国产视频亚洲 婷婷伊人久久 国产一区二区免费播放 久久99国产精品成人欧美 99国产在线视频 国产成人免费视频精品一区二区 国产不卡一区二区三区免费视 国产码欧美日韩高清综合一区 久久精品国产主播一区二区 国产一区电影 久久精品国产夜色 国产精品国产三级国产 日韩一区二区三区在线 久久97久久97精品免视看 久久国产免费一区二区三区 伊人久久大香线蕉综合电影网 99re6久精品国产首页 久久激情网 亚洲成人高清在线 国产精品网址 国产成人精品男人的天堂538 香蕉国产综合久久猫咪 国产专区中文字幕 91麻豆精品国产高清在线 久久国产经典视频 国产精品成人va在线观看 国产精品爱啪在线线免费观看 日本精品久久久久久久久免费 亚洲综合一区二区三区 久久五月网 精品国产网红福利在线观看 久久综合亚洲伊人色 亚洲国产精品久久久久久网站 在线日韩国产 99国产精品热久久久久久夜夜嗨 国产综合精品在线 国产区福利 精品亚洲综合久久中文字幕 国产制服丝袜在线 毛片在线播放网站 在线观看免费视频一区 国产精品久久久精品三级 亚洲国产电影在线观看 最新日韩欧美不卡一二三区 狠狠综合久久综合鬼色 日本精品1在线区 国产日韩一区二区三区在线播放 欧美日韩精品在线播放 亚洲欧美日韩国产一区二区三区精品 久久综合久久网 婷婷六月激情在线综合激情 亚洲乱码一区 国产专区91 97av视频在线观看 精品久久久久久中文字幕 久久五月视频 国产成人福利精品视频 国产精品网址 中文字幕视频在线 精品一区二区三区免费视频 伊人手机在线视频 亚洲精品中文字幕乱码 国产在线视频www色 色噜噜国产精品视频一区二区 精品亚洲成a人在线观看 国产香蕉尹人综合在线 成人免费一区二区三区在线观看 国产不卡一区二区三区免费视 欧美精品久久天天躁 国产专区中文字幕 久久精品国产免费中文 久久精品国产免费一区 久久无码精品一区二区三区 国产欧美另类久久久精品免费 欧美精品久久天天躁 亚洲精品在线视频 国产视频91在线 91精品福利一区二区三区野战 日韩中文字幕免费 国产精品99一区二区三区 欧美成人高清性色生活 国产精品系列在线观看 亚洲国产福利精品一区二区 国产成人在线小视频 国产精品久久久久免费 99re热在线视频 久久久久久久综合 一区二区国产在线播放 成人国产在线视频 亚洲精品乱码久久久久 欧美日韩一区二区综合 精品久久久久免费极品大片 中文字幕视频二区 激情粉嫩精品国产尤物 国产成人精品一区二区视频 久久精品中文字幕首页 亚洲高清在线 国产精品亚洲一区二区三区 伊人久久艹 中文在线亚洲 国产精品一区二区在线播放 国产精品九九免费视频 亚洲二区在线播放 亚洲狠狠婷婷综合久久久久网站 亚洲欧美日韩网站 日韩成人精品 亚洲国产一区二区三区青草影视 91精品国产福利在线观看 国产精品久久久久久久久99热 国产一区二区精品尤物 久碰香蕉精品视频在线观看 亚洲日日夜夜 在线不卡一区二区 国产午夜亚洲精品 九九热在线视频观看这里只有精品 伊人手机在线视频 91免费国产精品 日韩欧美中字 91精品国产91久久久久 国产全黄三级播放 视频一区二区三区免费观看 国产开裆丝袜高跟在线观看 国产成人欧美 激情综合丝袜美女一区二区 国产成人亚洲综合无 欧美精品一区二区三区免费观看 欧美亚洲国产日韩 日韩亚州 国产欧美日韩精品高清二区综合区 亚洲午夜国产片在线观看 精品久久久久久中文字幕 欧美精品1区 久久伊人久久亚洲综合 亚洲欧美日韩精品 国产成人精品久久亚洲高清不卡 久久福利影视 国产精品99精品久久免费 久久久久免费精品视频 国产日产亚洲精品 亚洲国产午夜电影在线入口 精品无码一区在线观看 午夜国产精品视频 亚洲一级片免费 伊人久久大香线蕉综合影 国产精品久久影院 久碰香蕉精品视频在线观看 www.欧美精品 在线小视频国产 亚洲国产天堂久久综合图区 欧美一区二区三区不卡 日韩美女福利视频 九九精品免视频国产成人 不卡国产00高中生在线视频 亚洲第一页在线视频 欧美日韩在线播放成人 99re视频这里只有精品 国产精品91在线 精品乱码一区二区三区在线 国产区久久 91麻豆精品国产自产在线观看一区 日韩精品成人在线 九九热在线观看 国产精品久久不卡日韩美女 欧美一区二区三区综合色视频 欧美精品免费一区欧美久久优播 国产精品网址 国产专区中文字幕 国产精品欧美亚洲韩国日本久久 日韩美香港a一级毛片 久久精品123 欧美一区二区三区免费看 99r在线视频 亚洲精品国产字幕久久vr 国产综合激情在线亚洲第一页 91免费国产精品 日韩免费小视频 亚洲国产精品综合一区在线 国产亚洲第一伦理第一区 在线亚洲精品 国产精品一区二区制服丝袜 国产在线成人精品 九九精品免视频国产成人 亚洲国产网 欧美日韩亚洲一区二区三区在线观看 在线亚洲精品 欧美一区二区三区高清视频 国产成人精品男人的天堂538 欧美日韩在线观看区一二 亚洲欧美一区二区久久 久久精品中文字幕首页 日本高清www午夜视频 久久精品国产免费 久久999精品 亚洲国产精品欧美综合 88国产精品视频一区二区三区 91久久偷偷做嫩草影院免费看 国产精品夜色视频一区二区 欧美日韩导航 国产成人啪精品午夜在线播放 一区二区视频在线免费观看 99久久精品国产自免费 精液呈暗黄色 久久99国产精品 日本精品久久久久久久久免费 精品国产97在线观看 99re视频这里只有精品 国产视频91在线 999av视频 亚洲美女视频一区二区三区 久久97久久97精品免视看 亚洲国产成人久久三区 99久久亚洲国产高清观看 日韩毛片在线视频 综合激情在线 91福利一区二区在线观看 一区二区视频在线免费观看 激情粉嫩精品国产尤物 国产成人精品曰本亚洲78 国产成人精品本亚洲 国产精品成人免费视频 国产成人啪精品视频免费软件 久久精品国产亚洲妲己影院 国产精品成人久久久久久久 久久大香线蕉综合爱 欧美一区二区三区高清视频 99热国产免费 在线观看欧美国产 91精品视频在线播放 国产精品福利社 欧美精品一区二区三区免费观看 国产一区二区免费视频 国产午夜精品一区二区 精品视频在线观看97 91精品福利久久久 国产一区福利 国产综合激情在线亚洲第一页 国产精品久久久久久久久久久不卡 九色国产 在线日韩国产 黄网在线观看 亚洲一区小说区中文字幕 中文字幕丝袜 日本二区在线观看 日本国产一区在线观看 欧美日韩一区二区三区久久 欧美精品亚洲精品日韩专 国产日产亚洲精品 久久综合九色综合欧美播 亚洲国产欧美无圣光一区 欧美视频区 亚洲乱码视频在线观看 久久无码精品一区二区三区 九九热精品免费视频 久久99精品久久久久久牛牛影视 国产精品成久久久久三级 国产一区福利 午夜国产精品视频 日本二区在线观看 99久久网站 国产亚洲天堂 精品国产一区二区三区不卡 亚洲国产日韩在线一区 国产成人综合在线观看网站 久久免费高清视频 欧美在线导航 午夜精品久久久久久99热7777 欧美久久综合网 国产小视频精品 国产尤物在线观看 亚洲国产精品综合一区在线 欧美一区二区三区不卡视频 欧美黑人欧美精品刺激 日本福利在线观看 久久国产偷 国产手机精品一区二区 国产热re99久久6国产精品 国产高清啪啪 欧美亚洲国产成人高清在线 国产在线第三页 亚洲综合一区二区三区 99r在线视频 99精品久久久久久久婷婷 国产精品乱码免费一区二区 国产在线精品福利91香蕉 国产尤物视频在线 五月婷婷亚洲 中文字幕久久综合伊人 亚洲精品一级毛片 99国产精品电影 在线视频第一页 久久99国产精品成人欧美 国产白白视频在线观看2 成人精品一区二区www 亚洲成人网在线观看 麻豆91在线视频 色综合合久久天天综合绕视看 久久精品国产免费高清 国产不卡一区二区三区免费视 欧美国产中文 99精品欧美 九九在线精品 国产中文字幕在线免费观看 国产一区中文字幕在线观看 国产成人一级 国产精品一区二区制服丝袜 国产一起色一起爱 亚洲精品成人在线 亚洲欧美精品在线 国产欧美自拍视频 99精品久久久久久久婷婷 久99视频 国产热re99久久6国产精品 视频一区亚洲 国产精品视频分类 国产精品成在线观看 99re6久精品国产首页 亚洲在成人网在线看 亚洲国产日韩在线一区 久久国产三级 日韩国产欧美 欧美在线一区二区三区 国产精品美女一级在线观看 成人午夜免费福利视频 亚洲天堂精品在线 91精品国产手机 欧美日韩视频在线播放 狠狠综合久久综合鬼色 九一色视频 青青视频国产 亚洲欧美自拍一区 中文字幕天天躁日日躁狠狠躁97 日韩免费大片 996热视频 伊人成综合网 亚洲天堂欧美 日韩精品亚洲人成在线观看 久久综合给会久久狠狠狠 日韩精品亚洲人成在线观看 日韩国产欧美 亚洲成aⅴ人片在线影院八 亚洲精品1区 99久久精品免费 国产精品高清在线观看 国产精品久久久免费视频 在线亚洲欧美日韩 91在线看视频 国产精品96久久久久久久 欧美日韩国产成人精品 91在线亚洲 热久久亚洲 国产精品美女免费视频观看 日韩在线毛片 亚洲永久免费视频 九九免费在线视频 亚洲一区网站 日本高清二区视频久二区 精品国产美女福利在线 伊人久久艹 国产精品久久久久三级 欧美成人精品第一区二区三区 99久久精品国产自免费 在线观看日韩一区 国产中文字幕一区 成人免费午夜视频 欧美日韩另类在线 久久99国产精品成人欧美 色婷婷中文网 久久天天躁夜夜躁狠狠躁2020 欧美成人伊人久久综合网 国产精品福利资源在线 国产伦精品一区二区三区高清 国产精品亚洲综合色区韩国 亚洲一区欧美日韩 色综合视频 国语自产精品视频在线区 国产高清a 成人国内精品久久久久影 国产在线精品香蕉综合网一区 国产不卡在线看 国产成人精品精品欧美 国产欧美日韩综合精品一区二区三区 韩国电影一区二区 国产在线视频www色 91中文字幕在线一区 国产人成午夜免视频网站 亚洲综合一区二区三区 色综合视频一区二区观看 久久五月网 九九热精品在线观看 国产一区二区三区国产精品 99久热re在线精品996热视频 亚洲国产网 在线视频亚洲一区 日韩字幕一中文在线综合 国产高清一级毛片在线不卡 精品国产色在线 国产高清视频一区二区 精品日本久久久久久久久久 亚洲国产午夜精品乱码 成人免费国产gav视频在线 日韩欧美一区二区在线观看 欧美曰批人成在线观看 韩国电影一区二区 99re这里只有精品6 日韩精品一区二区三区视频 99re6久精品国产首页 亚洲欧美一区二区三区导航 欧美色图一区二区三区 午夜精品视频在线观看 欧美激情在线观看一区二区三区 亚洲热在线 成人国产精品一区二区网站 亚洲一级毛片在线播放 亚洲一区小说区中文字幕 亚洲午夜久久久久影院 国产自产v一区二区三区c 国产精品视频免费 久久调教视频 国产成人91激情在线播放 国产精品欧美亚洲韩国日本久久 久久亚洲日本不卡一区二区 91中文字幕网 成人国产在线视频 国产视频91在线 欧美成人精品第一区二区三区 国产精品福利在线 久久综合九色综合精品 欧美一区二区三区精品 久久国产综合尤物免费观看 久久99青青久久99久久 日韩精品免费 久久国产精品999 91亚洲视频在线观看 国产精品igao视频 色综合区 在线亚洲欧国产精品专区 国产一区二区三区在线观看视频 亚洲精品成人在线 一区二区国产在线播放 中文在线亚洲 亚洲精品第一国产综合野 国产一区二区精品久久 一区二区三区四区精品视频 99热精品久久 中文字幕视频二区 国产成人精品男人的天堂538 99精品影视 美女福利视频一区二区 久久午夜夜伦伦鲁鲁片 综合久久久久久久综合网 国产精品国产欧美综合一区 国产99视频在线观看 国产亚洲女在线精品 婷婷影院在线综合免费视频 国产亚洲3p一区二区三区 91成人爽a毛片一区二区 亚洲一区二区高清 国产欧美亚洲精品第二区首页 欧美日韩导航 亚洲高清二区 欧美激情观看一区二区久久 日韩毛片在线播放 亚洲欧美日韩高清中文在线 亚洲日本在线播放 国产精品一区二区制服丝袜 精品国产一区二区三区不卡 国产不卡在线看 国产欧美网站 四虎永久在线观看视频精品 国产黄色片在线观看 夜夜综合 一本色道久久综合狠狠躁篇 欧美亚洲综合另类在线观看 国产91在线看 伊人久久国产 欧美一区二区在线观看免费网站 国产精品久久久久三级 久久福利 日韩中文字幕a 亚洲午夜久久久久影院 91在线高清视频 国产亚洲一区二区三区啪 久久人精品 国产精品亚洲午夜一区二区三区 综合久久久久久 久久伊人一区二区三区四区 国产综合久久久久久 日韩一区精品视频在线看 国产精品日韩欧美制服 日本精品1在线区 99re视频 无码av免费一区二区三区试看 国产视频1区 日韩欧美中文字幕一区 日本高清中文字幕一区二区三区a 亚洲国产欧美无圣光一区 国产在线视频一区二区三区 欧美国产第一页 在线亚洲欧美日韩 日韩中文字幕第一页 在线不卡一区二区 伊人久久青青 国产精品一区二区在线播放 www.五月婷婷 麻豆久久婷婷国产综合五月 亚洲精品区 久久国产欧美另类久久久 99在线视频免费 伊人久久中文字幕久久cm 久久精品成人免费看 久久这里只有精品首页 88国产精品视频一区二区三区 中文字幕日本在线mv视频精品 国产在线精品成人一区二区三区 伊人精品线视天天综合 亚洲一区二区黄色 国产尤物视频在线 亚洲精品99久久久久中文字幕 国产一区二区三区免费观看 伊人久久大香线蕉综合电影网 国产成人精品区在线观看 日本精品一区二区三区视频 日韩高清在线二区 久久免费播放视频 一区二区成人国产精品 国产精品免费精品自在线观看 亚洲精品视频二区 麻豆国产精品有码在线观看 精品日本一区二区 亚洲欧洲久久 久久中文字幕综合婷婷 中文字幕视频在线 国产成人精品综合在线观看 91精品国产91久久久久福利 精液呈暗黄色 香蕉国产综合久久猫咪 国产专区精品 亚洲精品无码不卡 国产永久视频 亚洲成a人片在线播放观看国产 一区二区国产在线播放 亚洲一区二区黄色 欧美日韩在线观看视频 亚洲精品另类 久久国产综合尤物免费观看 国产一区二区三区国产精品 高清视频一区 国产精品igao视频 国产精品资源在线 久久综合精品国产一区二区三区 www.五月婷婷 精品色综合 99热国产免费 麻豆福利影院 亚洲伊人久久大香线蕉苏妲己 久久电影院久久国产 久久精品伊人 在线日韩理论午夜中文电影 亚洲国产欧洲综合997久久 伊人国产精品 久草国产精品 欧美一区精品二区三区 亚洲成人高清在线 91免费国产精品 日韩精品福利在线 国产一线在线观看 国产不卡在线看 久久99青青久久99久久 亚洲精品亚洲人成在线播放 99久久免费看国产精品 国产日本在线观看 青草国产在线视频 麻豆久久婷婷国产综合五月 国产中文字幕一区 91久久精品国产性色也91久久 国产一区a 国产欧美日韩成人 国产亚洲女在线精品 一区二区美女 中文字幕在线2021一区 在线小视频国产 久久这里只有精品首页 国产在线第三页 欧美日韩中文字幕 在线亚洲+欧美+日本专区 精品国产一区二区三区不卡 久久这里精品 欧美在线va在线播放 精液呈暗黄色 91精品国产手机 91在线免费播放 欧美视频亚洲色图 欧美国产日韩精品 日韩高清不卡在线 精品视频免费观看 欧美日韩一区二区三区四区 国产欧美亚洲精品第二区首页 亚洲韩精品欧美一区二区三区 国产精品视频免费 在线精品小视频 久久午夜夜伦伦鲁鲁片 国产无套在线播放 久热这里只精品99re8久 欧美久久久久 久久香蕉国产线看观看精品蕉 国产成人精品男人的天堂538 亚洲人成网站色7799在线观看 日韩在线第二页 一本色道久久综合狠狠躁篇 国产一区二区三区不卡在线观看 亚洲乱码在线 在线观看欧美国产 久久福利青草精品资源站免费 国产玖玖在线观看 在线亚洲精品 亚洲成aⅴ人在线观看 精品91在线 欧美一区二三区 日韩中文字幕视频在线 日本成人一区二区 日韩免费专区 国内精品在线观看视频 久久国产综合尤物免费观看 国产精品系列在线观看 一本一道久久a久久精品综合 亚洲免费播放 久久精品国产免费 久久人精品 亚洲毛片网站 亚洲成a人一区二区三区 韩国福利一区二区三区高清视频 亚洲精品天堂在线 一区二区三区中文字幕 亚洲国产色婷婷精品综合在线观看 亚洲国产成人久久笫一页 999国产视频 国产精品香港三级在线电影 欧美日韩一区二区三区四区 日韩国产欧美 国产精品99一区二区三区 午夜国产精品理论片久久影院 亚洲精品中文字幕麻豆 亚洲国产高清视频 久久免费手机视频 日韩a在线观看 五月婷婷亚洲 亚洲精品中文字幕麻豆 中文字幕丝袜 www国产精品 亚洲天堂精品在线 亚洲乱码一区 国产日韩欧美三级 久久999精品 伊人热人久久中文字幕 久热国产在线视频 国产欧美日韩在线观看一区二区三区 国产一二三区在线 日韩国产欧美 91精品国产91久久久久 亚洲一区小说区中文字幕 精品一区二区免费视频 国产精品视频免费 国产精品亚洲综合色区韩国 亚洲国产精品成人午夜在线观看 欧美国产日韩精品 中文字幕精品一区二区精品