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

        ?

        基于產(chǎn)品DSM 自動(dòng)構(gòu)建的復(fù)雜產(chǎn)品模塊劃分方法

        2023-12-22 06:47:46諶炎輝陳劍濤班華電梁帆陳炳森
        裝備制造技術(shù) 2023年10期
        關(guān)鍵詞:產(chǎn)品方法

        諶炎輝,陳劍濤,班華電,梁帆,陳炳森

        (1.廣西科技大學(xué) 機(jī)械與汽車工程學(xué)院,廣西 柳州 545006;2.廣西水利電力職業(yè)技術(shù)學(xué)院 機(jī)電工程學(xué)院,廣西 南寧 530105;3.廣西徐沃工程機(jī)械設(shè)備有限公司,廣西 南寧 530105)

        產(chǎn)品模塊劃分方法一直是國(guó)內(nèi)外學(xué)者研究的熱點(diǎn)問題之一,Stone 等[1]基于能量流和信息流提出過一種啟發(fā)式的模塊劃分方法,Tseng 等[2]提出從零件的接口類型和特征出發(fā)對(duì)產(chǎn)品進(jìn)行模塊劃分的思想,Lee等[3]人從綠色環(huán)保的角度出發(fā),對(duì)產(chǎn)品零部件做了獨(dú)立性分析建立綜合交互矩陣。文獻(xiàn)[4-5]在對(duì)產(chǎn)品模塊化研究中分別提出了一種基于產(chǎn)品物料清單(BOM)的適用于復(fù)雜產(chǎn)品的模塊劃分方法和最大最小模塊劃分的方法,主要分析產(chǎn)品內(nèi)部零部件間的幾何、物理、功能相關(guān)性對(duì)產(chǎn)品進(jìn)行模塊劃分。程賢福等[6]針對(duì)產(chǎn)品的模塊劃分提出了適用于對(duì)稱和非對(duì)稱的數(shù)值型或布爾型設(shè)計(jì)結(jié)構(gòu)矩陣的模塊劃分準(zhǔn)則;姜芊芊等[7]在應(yīng)對(duì)復(fù)雜機(jī)電產(chǎn)品模塊劃分難度大的特點(diǎn)提出了通過功能-行為-結(jié)構(gòu)映射模型分析產(chǎn)品功能和結(jié)構(gòu)的關(guān)聯(lián)關(guān)系;文獻(xiàn)[8-9]提出從產(chǎn)品生命周期和綠色制造的角度出發(fā)對(duì)產(chǎn)品進(jìn)行模塊劃分。蘇銘等[10]結(jié)合計(jì)算機(jī)輔助設(shè)計(jì)技術(shù)提出了模塊自動(dòng)裝配策略,并以此為基礎(chǔ)開發(fā)了數(shù)控機(jī)床模塊化設(shè)計(jì)系統(tǒng),顯著提高了產(chǎn)品模塊化設(shè)計(jì)周期;李中凱等[11]通過分析產(chǎn)品的結(jié)構(gòu)功能構(gòu)建了產(chǎn)品的量化信息模型,通過分配量化指標(biāo)不同的權(quán)重來獲得產(chǎn)品各組件之間的綜合關(guān)聯(lián)矩陣。顧新建等[12]利用群體智能、大數(shù)據(jù)智能、人機(jī)協(xié)同智能等技術(shù),對(duì)大量產(chǎn)品模塊數(shù)據(jù)進(jìn)行智能分析,能夠識(shí)別零部件間的一些關(guān)聯(lián)關(guān)系來幫助快速完成產(chǎn)品的模塊劃分工作。此前已有相關(guān)學(xué)者提出通過對(duì)設(shè)計(jì)結(jié)構(gòu)樹進(jìn)行訪問從中得到裝配體組件間的約束配合類型進(jìn)而判斷組件間是否存在連接關(guān)系構(gòu)建產(chǎn)品設(shè)計(jì)結(jié)構(gòu)矩陣(Design Structure Matrix:DSM)用于產(chǎn)品的模塊劃分工作[13]。但是這種方法存在一定的缺陷:其一,設(shè)計(jì)結(jié)構(gòu)樹中記錄的裝配體組件間的配合關(guān)系不全;其二,僅憑組件間的配合關(guān)系難以判斷其連接關(guān)系;其三,這種方法的計(jì)算量龐大,且需要耗費(fèi)一定的時(shí)間進(jìn)行人工分析,不利于產(chǎn)品的快速設(shè)計(jì)研發(fā)。

        綜上,在研究中參考前人的研究成果的基礎(chǔ)上,提出了一種基于產(chǎn)品DSM 自動(dòng)構(gòu)建的復(fù)雜產(chǎn)品模塊劃分方法。該方法通過對(duì)三維制圖軟件進(jìn)行二次開發(fā)直接獲取產(chǎn)品零部件間的空間距離以及產(chǎn)品的屬性名稱并自動(dòng)判斷零部件間的連接關(guān)系從而自動(dòng)獲取產(chǎn)品的DSM 矩陣,解決了傳統(tǒng)模塊劃分方法全依賴人工賦值導(dǎo)致模塊劃分方案主觀性過強(qiáng)的問題。對(duì)比現(xiàn)有的模塊劃分方法操作性強(qiáng)、工作量小,模塊劃分結(jié)果相對(duì)客觀,一定程度上提高了模塊劃分的準(zhǔn)確性與合理性。且該方法中更多地使用了計(jì)算機(jī)代替人工計(jì)算,使產(chǎn)品模塊劃分智能化水平得到了提高,在一定程度上縮短了產(chǎn)品模塊劃分的設(shè)計(jì)周期,減少了在對(duì)復(fù)雜機(jī)電產(chǎn)品進(jìn)行模塊劃分時(shí)容易出現(xiàn)錯(cuò)誤的情況。

        1 產(chǎn)品DSM 的自動(dòng)構(gòu)建

        1.1 基于零部件連接關(guān)系的設(shè)計(jì)結(jié)構(gòu)矩陣

        常見DSM 有兩種類型:數(shù)值型與布爾型,在產(chǎn)品模塊化的過程中數(shù)值型DSM 在判斷零部件間的連接關(guān)系時(shí)基本依靠人工賦值,導(dǎo)致賦值結(jié)果帶有一定的主觀性,因此在實(shí)驗(yàn)中采用布爾型DSM 對(duì)零部件的關(guān)聯(lián)關(guān)系進(jìn)行判斷,使模塊劃分結(jié)果更具有客觀說服力。

        1.2 基于NX 二次開發(fā)的產(chǎn)品DSM 自動(dòng)構(gòu)建

        研究中所使用的二次開發(fā)平臺(tái)版本為NX10.0,二次開發(fā)命令中測(cè)量?jī)蓪?shí)體間的距離的函數(shù)是UF_MODL_ask_minimum_dist。通過三維裝配體模型自動(dòng)構(gòu)建基于產(chǎn)品零部件連接關(guān)系的DSM 需要用到的輔助軟件是Visual Studio 和Excel,其中前者主要用來搭建編程環(huán)境、生成程序插件等功能,后者的作用是記錄產(chǎn)品屬性名稱、組件連接關(guān)系等數(shù)據(jù)。自動(dòng)構(gòu)建產(chǎn)品DSM 流程,如圖1 所示。

        圖1 自動(dòng)構(gòu)建產(chǎn)品DSM 流程圖

        流程步驟可以概括為:

        1)用NX 打開所需測(cè)量產(chǎn)品裝配體模型點(diǎn)擊功能菜單并對(duì)模型進(jìn)行框選;

        2)從裝配體模型中提取各組件的名稱以及屬性編號(hào)等數(shù)據(jù),并對(duì)組件從1 到N排序,在Excel 表格中創(chuàng)建一個(gè)行數(shù)與列數(shù)均為產(chǎn)品零部件數(shù)N的方形矩陣A;

        3)設(shè)置參數(shù)i=1,設(shè)置判斷閾值x;

        4)設(shè)置參數(shù)j=i;

        5)判斷條件i>N,若i>N則轉(zhuǎn)至步驟11,若i<=N則進(jìn)入下一步驟;

        6)判斷條件i=j,若i=j則替換矩陣A[i,j]為1 在轉(zhuǎn)至步驟9,若i≠j則進(jìn)入下一步驟;

        7)計(jì)算組件i與組件j之間的距離d;

        8)判斷條件d>x,若d>x則替換矩陣A[i,j]為0在轉(zhuǎn)至下一步驟,若d<=x則替換矩陣A[i,j]為1 在轉(zhuǎn)至下一步驟;

        9)參數(shù)j=j+1;

        10)判斷條件j>N,若j>N則參數(shù)i=i+1 后轉(zhuǎn)至步驟4,若j<=N則轉(zhuǎn)至步驟6;

        11)將各組件名稱以及所得數(shù)據(jù)保存至Excel 中。

        這種方法不僅僅只適用于NX 的二次開發(fā),在AutoCAD、SolidWorks 等三維設(shè)計(jì)軟件中也同樣適用,且開發(fā)思路基本一致,只需替換不同的軟件所對(duì)應(yīng)的開發(fā)函數(shù)即可。

        2 模塊劃分方法

        2.1 其他模塊劃分影響因素

        把產(chǎn)品的DSM 作為產(chǎn)品模塊劃分唯一的參考數(shù)據(jù)是不嚴(yán)謹(jǐn)?shù)?,在?duì)產(chǎn)品進(jìn)行模塊劃分時(shí)還應(yīng)該考慮到其他因素的影響,在查詢了一些相關(guān)文獻(xiàn)后發(fā)現(xiàn)學(xué)者在研究機(jī)電產(chǎn)品的模塊劃分確定產(chǎn)品零部件的關(guān)聯(lián)性時(shí)通常會(huì)從功能相關(guān)性、裝配關(guān)系、能量接口、幾何相關(guān)性以及物理相關(guān)性等方面分析[4,14-15]。由于上述通過對(duì)NX 的二次開發(fā)自動(dòng)獲取的產(chǎn)品DSM 矩陣已經(jīng)能夠體現(xiàn)產(chǎn)品零部件的連接關(guān)系,在一定程度上能夠表達(dá)產(chǎn)品的裝配關(guān)系以及幾何相關(guān)性等信息,考慮到模塊劃分方法應(yīng)該盡可能的簡(jiǎn)潔便于操作,通過咨詢相關(guān)專家的意見以及小組討論研究最終選取產(chǎn)品零部件的功能相關(guān)性結(jié)合產(chǎn)品的DSM 作為產(chǎn)品模塊劃分的依據(jù),并總結(jié)出如表1 所示的產(chǎn)品零部件模糊關(guān)聯(lián)度值。

        表1 產(chǎn)品零部件模糊關(guān)聯(lián)度值

        產(chǎn)品的零部件數(shù)量越多對(duì)其進(jìn)行模塊劃分的過程就越困難,過程中產(chǎn)生的計(jì)算量也會(huì)大幅度提高,通過前往工廠實(shí)地考察并咨詢了相關(guān)工程師后總結(jié)出以下幾條針對(duì)模塊劃分時(shí)產(chǎn)品零部件的簡(jiǎn)化原則:

        整體焊接件簡(jiǎn)化原則:產(chǎn)品裝配時(shí)焊接在一起無法分開的零部件在模塊劃分時(shí)看作一個(gè)部件處理。

        特殊整體模塊簡(jiǎn)化原則:特殊外購(gòu)件如馬達(dá)、電機(jī)等在裝配時(shí)不會(huì)再進(jìn)行拆分的零部件在模塊劃分時(shí)看作一個(gè)部件處理。

        國(guó)標(biāo)零件弱化原則:在產(chǎn)品生產(chǎn)時(shí)國(guó)標(biāo)零件通常會(huì)進(jìn)行大批量的生產(chǎn)且可替換性高,因此在模塊劃分時(shí)作弱化處理,不考慮國(guó)標(biāo)零件的模塊劃分。

        2.2 關(guān)聯(lián)度矩陣求解

        通過完成上述步驟分別可以得到產(chǎn)品的DSM 以及產(chǎn)品零部件的功能相關(guān)性矩陣將其分別記為N1、N2,接下來將會(huì)對(duì)這兩個(gè)矩陣進(jìn)行分析處理進(jìn)而得出產(chǎn)品零部件間的綜合關(guān)聯(lián)度。

        綜合關(guān)聯(lián)度得到矩陣計(jì)算式為:

        式中:N為計(jì)算所得的產(chǎn)品綜合關(guān)聯(lián)度矩陣;N1為產(chǎn)品的DSM;N2為產(chǎn)品功能相關(guān)性矩陣;ω1、ω2分別為上述兩個(gè)矩陣的關(guān)聯(lián)權(quán)重。

        通過式(1)求解所得的綜合關(guān)聯(lián)度矩陣:

        式中:α12表示第一個(gè)零件與第二個(gè)零件的綜合關(guān)聯(lián)度值;αij依次表示第i個(gè)零件與第j個(gè)零件的綜合關(guān)聯(lián)度值;對(duì)角線元素為零件與自身的關(guān)聯(lián)度值記為1。

        2.3 模塊劃分

        對(duì)矩陣進(jìn)行求解的目的是要對(duì)產(chǎn)品的零部件進(jìn)行模塊劃分,針對(duì)這一目的研究中采取層次聚類的方法對(duì)產(chǎn)品零部件綜合關(guān)聯(lián)度矩陣進(jìn)行求解,層次聚類在分類前不需預(yù)先設(shè)置聚類數(shù)目可生成層次聚類樹狀圖,能夠清晰直觀地展示不同類之間的層次關(guān)系,因此層次聚類適合用于產(chǎn)品的模塊劃分。完成層次聚類得到層次聚類樹狀圖后只需選擇不同的聚類閾值即可得到不同的聚類結(jié)果,依據(jù)模塊劃分的原則:模塊內(nèi)部相似度盡可能高,模塊間相似度盡可能低,提出了兩個(gè)模塊劃分評(píng)價(jià)指標(biāo):模塊內(nèi)平均聚合度D1以及模塊間平均聚合度D2。

        模塊內(nèi)平均聚合度為:

        式中:n為模塊劃分個(gè)數(shù);nk表示第k個(gè)模塊包含的零部件個(gè)數(shù)表示第k個(gè)模塊內(nèi)元素組成的模塊內(nèi)部關(guān)聯(lián)度矩陣第i行第j列的元素值。

        模塊間平均聚合度D2

        式中:n為模塊劃分個(gè)數(shù);α(Mi,Mj)表示子模塊Mi和Mj零部件之間總的關(guān)聯(lián)度值;ni和nj分別表示子模塊Mi和Mj包含的零部件個(gè)數(shù)。

        模塊劃分?jǐn)?shù)量會(huì)影響到產(chǎn)品模塊劃分的質(zhì)量,有研究表明模塊劃分的最佳模塊數(shù)量應(yīng)該在附近[16],m表示參與模塊劃分的總元素個(gè)數(shù),在模塊劃分時(shí)圍繞上下分別選取兩組以上模塊劃分方案,通過分析其模塊內(nèi)平均聚合度D1和模塊間平均聚合度D2綜合判斷出最佳模塊劃分方案。

        3 應(yīng)用實(shí)例

        在實(shí)驗(yàn)中以某企業(yè)生產(chǎn)的伐木機(jī)為試驗(yàn)對(duì)象,對(duì)其進(jìn)行模塊劃分工作以驗(yàn)證上述模塊劃分方法的可行性,按照上述2.1 中的產(chǎn)品模塊劃分簡(jiǎn)化原則對(duì)其零部件進(jìn)行簡(jiǎn)化處理后剩余零部件數(shù)量為76 個(gè)。

        為上述1.1 中所設(shè)計(jì)開發(fā)的基于產(chǎn)品三維裝配模型的產(chǎn)品DSM 自動(dòng)構(gòu)建程序的開始界面,如圖2所示。圖2 第一個(gè)選項(xiàng)為選取需要構(gòu)建的三維裝配體模型,可以逐個(gè)選取也可利用鼠標(biāo)進(jìn)行框選,第二個(gè)選項(xiàng)為選擇導(dǎo)出產(chǎn)品DSM 表格的文件保存路徑。

        圖2 DSM 自動(dòng)導(dǎo)出程序

        考慮到產(chǎn)品零部件存在裝配公差以及一些間隙連接的情況,因此設(shè)置0.1 mm 為伐木機(jī)零部件連接關(guān)系判斷閾值,首先通過產(chǎn)品的簡(jiǎn)化三維模型自動(dòng)獲取產(chǎn)品的DSM 表格以及自動(dòng)讀取產(chǎn)品的零部件名稱并對(duì)其進(jìn)行排序,處理結(jié)果如圖3、圖4 所示。圖3 為產(chǎn)品零部件序號(hào)對(duì)照?qǐng)D,圖4 中第一行與第一列為產(chǎn)品各零部件的序號(hào),圖中對(duì)應(yīng)表格為相關(guān)零部件之間的關(guān)聯(lián)關(guān)系判斷,其中著色單元格的數(shù)值為1 表示對(duì)應(yīng)零部件存在連接關(guān)系,不著色單元格數(shù)值為0 表示對(duì)應(yīng)零部件不存在連接關(guān)系,參考表中數(shù)據(jù)可建立產(chǎn)品的DSM 矩陣。

        圖3 伐木機(jī)零部件對(duì)照

        圖4 伐木機(jī)DSM

        根據(jù)表1 產(chǎn)品零部件模糊關(guān)聯(lián)度值和產(chǎn)品零部件排序建立產(chǎn)品功能相關(guān)性矩陣將其記為N2,分別選取權(quán)重ω1=0.5、ω2=0.5 通過式(1)對(duì)矩陣N1、N2進(jìn)行處理得到產(chǎn)品綜合關(guān)聯(lián)度矩陣。

        將綜合關(guān)聯(lián)度矩陣N導(dǎo)入MATLAB 中利用層次聚類算法求解伐木機(jī)零部件關(guān)聯(lián)度矩陣,得到如圖5所示的聚類樹狀圖。

        圖5 伐木機(jī)零部件聚類樹狀圖

        本次模塊劃分試驗(yàn)中伐木機(jī)包含的零部件數(shù)量為76 個(gè),根據(jù)上述2.3 中分析可知最佳模塊數(shù)取值應(yīng)該在左右,故選取模塊劃分?jǐn)?shù)量為7、8、9、10、11 五組方案,計(jì)算各方案模塊內(nèi)平均聚合度和模塊間平均聚合度,由于模塊間聚合度計(jì)算所得數(shù)值較小,為了便于分析觀察繪圖數(shù)據(jù)為原計(jì)算數(shù)據(jù)的十倍如圖6 所示。

        圖6 模塊聚合度對(duì)比圖

        為了更加直觀分析對(duì)上述數(shù)據(jù)可用式(5)處理。

        式中:φ為模塊聚合度增長(zhǎng)率,xi表示模塊劃分?jǐn)?shù)為時(shí)模塊內(nèi)聚合度或模塊間聚合度參考值,xmax、xmin則表示上述各組模塊數(shù)據(jù)中模塊內(nèi)聚合度或模塊間聚合度的最大值和最小值。

        計(jì)算各模塊數(shù)組模塊內(nèi)與模塊間平均聚合度增長(zhǎng)率φ,得到表2 的不同模塊數(shù)模塊劃分方案之間內(nèi)部模塊的平均聚合度變化情況。

        表2 聚合度增長(zhǎng)率

        從表2 可以看出隨著模塊劃分?jǐn)?shù)目的增加,模塊內(nèi)平均聚合度呈現(xiàn)持續(xù)增長(zhǎng)的趨勢(shì),這也是一個(gè)必然的結(jié)果,可以想象當(dāng)模塊劃分?jǐn)?shù)與產(chǎn)品零部件數(shù)相同時(shí)模塊內(nèi)平均聚合度即為零部件與自身的聚合度取得最大值1,然而這樣也就失去了模塊劃分的意義,因此要結(jié)合模塊間平均聚合度進(jìn)行綜合分析,由表中數(shù)據(jù)可知選取模塊劃分?jǐn)?shù)為10 時(shí)模塊間平均聚合度呈現(xiàn)下降趨勢(shì),而當(dāng)模塊劃分?jǐn)?shù)為11 時(shí)模塊間平均聚合度呈現(xiàn)顯著的上升趨勢(shì),結(jié)合圖8 可知,模塊劃分?jǐn)?shù)為10 時(shí),模塊間聚合度取值為0.1328,明顯低于9 組模塊和11 組模塊時(shí)的模塊間聚合度,此時(shí)模塊內(nèi)聚合度為0.4810,處于一個(gè)較高的水平,且10 較為接近最佳模糊聚類數(shù)。通過綜合考慮決定取10個(gè)做為最佳模塊劃分?jǐn)?shù),其模塊劃分明細(xì)為:{(2,29,48,53),(5,14,24,35,41,58,68),(7,50,56,64),(8,10,30,32,51,54,59),(9,16,25,36,52,60,62,67),(11,21,69,70,73,76),(3,13,18,26,33,34,38,40,43,61,63,66,71),(4,6,15,37,42,46),(1,12,17,22,23,28,55,57,65),(19,20,27,31,39,44,45,47,49,72,74,75)} 依次對(duì)應(yīng)伐木機(jī)的右側(cè)送料輪模塊、左側(cè)送料輪模塊、伐木機(jī)切割模塊、伐木機(jī)左側(cè)后夾刀模塊、伐木機(jī)右側(cè)后夾刀模塊、伐木機(jī)主體模塊、伐木機(jī)支撐旋轉(zhuǎn)模塊、伐木機(jī)中部送料輪模塊、伐木機(jī)前夾刀模塊、伐木機(jī)測(cè)距儀模塊。

        4 結(jié)語

        針對(duì)現(xiàn)有模塊劃分方法過于依靠人工賦值,不利于多零件復(fù)雜機(jī)電產(chǎn)品模塊劃分的特點(diǎn),提出了一種產(chǎn)品DSM 自動(dòng)構(gòu)建結(jié)合人工零部件關(guān)聯(lián)度分析的綜合模塊劃分方法,主要研究成果如下:

        (1)對(duì)NX10.0 進(jìn)行二次開發(fā)設(shè)計(jì)了一個(gè)功能模塊,該模塊能夠獲取產(chǎn)品三維裝配模型的空間位置、零部件名稱等相關(guān)信息,通過計(jì)算裝配模型中各零部件間的空間位置距離判斷零部件間是否存在接觸從而構(gòu)建產(chǎn)品的DSM。

        (2)提出了一種基于產(chǎn)品DSM 自動(dòng)構(gòu)建的復(fù)雜產(chǎn)品模塊劃分方法以及基于模塊內(nèi)外聚合度的模塊劃分評(píng)價(jià)指標(biāo),并應(yīng)用此方法將某型號(hào)伐木機(jī)76 個(gè)零部件合理劃分為十個(gè)模塊,驗(yàn)證了該方法的可行性。

        通過對(duì)產(chǎn)品模塊劃分方法的研究可以看出未來產(chǎn)品模塊的劃分必將朝著智能化自動(dòng)化的方向發(fā)展。同時(shí)也希望該研究成果能夠?qū)髞淼难芯坑幸欢ǖ膯l(fā)。

        猜你喜歡
        產(chǎn)品方法
        好產(chǎn)品,可持續(xù)
        從靈感出發(fā),邂逅好產(chǎn)品
        學(xué)習(xí)方法
        可能是方法不對(duì)
        用對(duì)方法才能瘦
        Coco薇(2016年2期)2016-03-22 02:42:52
        四大方法 教你不再“坐以待病”!
        Coco薇(2015年1期)2015-08-13 02:47:34
        2015產(chǎn)品LOOKBOOK直擊
        Coco薇(2015年1期)2015-08-13 02:23:50
        賺錢方法
        捕魚
        新產(chǎn)品
        玩具(2009年10期)2009-11-04 02:33:14
        久久久久人妻精品一区蜜桃| 最新国产精品国产三级国产av | 日韩中文字幕在线观看一区| 最近中文字幕免费完整版| 中文字幕无码精品亚洲资源网久久 | 国产三级精品美女三级| 日本女同av在线播放| 高清毛茸茸的中国少妇| 99久久精品费精品国产一区二区| 欧美日本道免费二区三区| 青青草成人原视频在线播放视频 | 欧美video性欧美熟妇| 欧美 日韩 人妻 高清 中文| 亚洲国产无套无码av电影| 青青草综合在线观看视频| 亚洲精品中文字幕一二三| 加勒比hezyo黑人专区| 国产乱人伦偷精品视频| 欧美人与动牲交片免费| 国产黄色三级一区二区三区四区| 夜夜躁狠狠躁日日躁2022| 玖玖资源站无码专区| 久久精品国产亚洲av热一区| 国产三级不卡一区不卡二区在线| 少妇性饥渴无码a区免费| 国产午夜在线观看视频播放| 国产午夜精品av一区二区三 | 狠狠色综合7777久夜色撩人| 国产乱人视频在线看| 亚洲av中文字字幕乱码| 日本道色综合久久影院| 国产成人精品日本亚洲11| 欧美日韩一区二区三区视频在线观看 | 亚洲精品亚洲人成在线下载| 国产激情小视频在线观看| 国产成人精品一区二区三区| 国产日韩欧美亚洲精品中字| 国产一区二区三区色区| 97精品人妻一区二区三区蜜桃| 国产97在线 | 免费| 亚洲国产成人无码电影|