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

        ?

        基于模糊多準(zhǔn)則的滌綸低彈絲生產(chǎn)工藝參數(shù)優(yōu)化

        2021-04-06 06:02:18邵景峰蔡再生
        紡織學(xué)報(bào) 2021年1期
        關(guān)鍵詞:收縮率滌綸賦權(quán)

        邵景峰, 李 寧, 蔡再生

        (1. 西安工程大學(xué) 管理學(xué)院, 陜西 西安 710048; 2. 東華大學(xué)化學(xué)化工與生物工程學(xué)院, 上海 201620)

        滌綸低彈絲(DTY)作為滌綸的一種變形絲類型,具有很高的取向度,不僅手感柔軟順滑,且具有很好的質(zhì)地,已成為化纖行業(yè)實(shí)現(xiàn)綠色制造的理想原料。然而,在滌綸低彈絲生產(chǎn)過程中,因工藝參數(shù)(如加工速度、熱箱溫度等)之間存在強(qiáng)耦合性而導(dǎo)致整個(gè)生產(chǎn)過程能耗過高,從而影響生產(chǎn)過程難以實(shí)現(xiàn)綠色制造;因此,如何對(duì)滌綸低彈絲工藝參數(shù)之間的強(qiáng)耦合關(guān)系進(jìn)行優(yōu)化,已成為亟待解決的問題。

        在國外,學(xué)者們對(duì)該問題的研究,主要集中在工藝?yán)碚搩?yōu)化與技術(shù)應(yīng)用兩方面。其中,在理論構(gòu)筑方面,利用生命周期理論、碳核算理論、能耗測量技術(shù)探討了滌綸低彈絲環(huán)境負(fù)荷與工藝參數(shù)之間關(guān)系,還有能耗實(shí)時(shí)測量方法,以及碳足跡核算與評(píng)價(jià)方法[1-3]等,目前已經(jīng)明晰了各參數(shù)之間的關(guān)聯(lián)關(guān)系,實(shí)現(xiàn)了滌綸低彈絲在給定范圍的生產(chǎn)質(zhì)量預(yù)測[4-6],為滌綸低碳生產(chǎn)工藝參數(shù)的優(yōu)化奠定了理論基礎(chǔ)。在技術(shù)應(yīng)用方面,主要集中于工藝參數(shù)優(yōu)化算法、優(yōu)化方法等方面,目前已解決了假捻變形參數(shù)對(duì)滌綸紗線結(jié)構(gòu)和卷曲性能的影響、滌綸低彈絲加捻變形參數(shù)之間的關(guān)系、以及基于VIKOR的滌綸低彈絲性能評(píng)價(jià)方法[7-8]。

        在我國,眾多學(xué)者也對(duì)工藝參數(shù)之間強(qiáng)耦合性問題進(jìn)行了研究,其研究焦點(diǎn)在于工藝參數(shù)優(yōu)化和生產(chǎn)實(shí)踐方面。其中,在工藝參數(shù)優(yōu)化方面,通過構(gòu)建滌綸低彈絲生產(chǎn)過程多目標(biāo)優(yōu)化調(diào)度模型,開展工藝參數(shù)的多目標(biāo)優(yōu)化[9-11],同時(shí)借助參數(shù)優(yōu)化方法(比如熱定型過程參數(shù)優(yōu)化)開展溫度控制、碳足跡核算,但還存在工藝參數(shù)強(qiáng)耦合性問題探討尚不深入的問題。在生產(chǎn)實(shí)踐方面,通過建立生產(chǎn)過程圖、能源消費(fèi)碳排放因素分解、優(yōu)化工藝方案等[12-14],開展面向滌綸低彈絲生產(chǎn)過程碳足跡核算、影響因素以及能耗損失分析,并借助試驗(yàn)和應(yīng)用驗(yàn)證的方法對(duì)工藝參數(shù)相關(guān)關(guān)系進(jìn)行改進(jìn)。

        目前,國內(nèi)外學(xué)者們對(duì)滌綸低彈絲生產(chǎn)工藝參數(shù)之間強(qiáng)耦合問題的研究,已明晰了各參數(shù)之間的關(guān)聯(lián)關(guān)系,解決了假捻、加捻變形參數(shù)關(guān)系、多目標(biāo)優(yōu)化、碳排放因素分解問題,但還存在工藝參數(shù)之間的強(qiáng)耦合性問題尚未徹底解決。該問題的解決,將對(duì)提升滌綸低彈絲成品性能、實(shí)現(xiàn)綠色制造有重要意義。為此,本文借助碳足跡理論從滌綸低彈絲生產(chǎn)工藝流程入手,提取關(guān)鍵工藝參數(shù)作為研究對(duì)象,構(gòu)建基于模糊多準(zhǔn)則的工藝參數(shù)優(yōu)化方法,以實(shí)現(xiàn)關(guān)鍵工藝參數(shù)間的解耦。

        1 滌綸低彈絲工藝參數(shù)關(guān)系表示

        1.1 滌綸低彈絲生產(chǎn)過程碳排放分析

        滌綸低彈絲生產(chǎn)主要是由高速紡制滌綸預(yù)取向絲(POY)通過牽伸加捻工藝而形成[15-17],首先將POY絲掛絲喂入拉伸變形機(jī)中,通過加熱箱對(duì)POY絲進(jìn)行加熱變形拉伸處理,然后將加熱變形后的絲經(jīng)過冷卻板冷卻后進(jìn)行加捻,同時(shí)對(duì)加捻處理后的絲進(jìn)行二次加熱拉伸,并通過上油處理減少成形后紡紗過程的靜電效應(yīng),最后經(jīng)卷繞、落絲以及分級(jí)質(zhì)檢等工序后出廠。

        在此基礎(chǔ)上,利用碳足跡理論對(duì)滌綸低彈絲生產(chǎn)過程中的碳排放信息進(jìn)行分析。具體而言,影響滌綸低彈絲品質(zhì)的關(guān)鍵工藝參數(shù)主要有加工速度、牽伸比、第一熱箱溫度、第二熱箱溫度、油輪轉(zhuǎn)速、油尺高度。在整個(gè)滌綸低彈絲生產(chǎn)過程中,碳排放信息主要通過電能的消耗進(jìn)行轉(zhuǎn)化,拉伸變形機(jī)中的變形拉伸以及各設(shè)備零器件之間的傳動(dòng)等是電能消耗的主要途徑[18]。

        為此,利用碳足跡理論,構(gòu)建如圖1所示的滌綸低彈絲生產(chǎn)過程碳排放分析模型。

        通過該模型,建立碳排放信息量Ce。

        (1)

        式中:Ei為生產(chǎn)過程中各項(xiàng)能源(i)消耗轉(zhuǎn)化的碳排放量;fi為各能源消耗對(duì)應(yīng)的碳排放系數(shù)。

        1.2 滌綸低彈絲工藝參數(shù)關(guān)系

        在式(1)的基礎(chǔ)上,選擇油輪轉(zhuǎn)速(x1)、加工速度(x2)、牽伸比(x3)、油尺高度(x4)、第一熱箱溫度(x5)、第二熱箱溫度(x6)為自變量,碳排放值Ce。為因變量,利用碳足跡理論建立因變量與自變量之間的函數(shù)關(guān)系。進(jìn)而以83 dtex(36 f)低彈絲為例,在FK6V-700型拉伸變形機(jī)上獲取的參數(shù)及指標(biāo)如表1~3所示。

        表1 DTY加工參數(shù)及性能指標(biāo)

        表2 不同含油率下DTY的線密度及斷裂強(qiáng)度

        表3 加工不同線密度DTY的油輪轉(zhuǎn)速及含油率

        基于表1~3中的工藝參數(shù),利用MatLab R2018b軟件對(duì)選擇的自變量與因變量關(guān)系進(jìn)行仿真試驗(yàn),得到的二次回歸方程為

        (2)

        2 滌綸低彈絲工藝參數(shù)關(guān)系優(yōu)化

        2.1 基于信噪比的工藝參數(shù)正交試驗(yàn)設(shè)計(jì)

        信噪比具有望大特性與望小特性2種屬性,通常作為衡量產(chǎn)品質(zhì)量穩(wěn)健性的重要指標(biāo)[19],將DTY指標(biāo)轉(zhuǎn)換為信噪比值,進(jìn)而評(píng)估DTY質(zhì)量品質(zhì)的優(yōu)劣。為此,在式(2)的基礎(chǔ)上,為表征DTY的工藝參數(shù)關(guān)系,還需將工藝參數(shù)和正交設(shè)計(jì)相互關(guān)聯(lián)。為此,根據(jù)信息論理論[20-21]將信噪比引入工藝參數(shù)關(guān)系正交試驗(yàn)設(shè)計(jì)。對(duì)DTY工藝參數(shù)而言,相互之間關(guān)聯(lián)度越高,意味著滌綸低彈絲成品的質(zhì)量指標(biāo)越好。DTY指標(biāo)轉(zhuǎn)換為信噪比值的計(jì)算過程如下。

        斷裂強(qiáng)度(Z1)、斷裂伸長率(Z2)、卷曲收縮率(Z3)、沸水收縮率(Z4)、卷曲穩(wěn)定度(Z5)作為DTY質(zhì)量指標(biāo),在一定范圍內(nèi)越大越好,因此采用望大特性的信噪比進(jìn)行數(shù)據(jù)轉(zhuǎn)換如式(3)所示。

        (3)

        式中:aij為轉(zhuǎn)換后的第i行第j列的DTY質(zhì)量指標(biāo)信噪比值;n為每組試驗(yàn)的重復(fù)次數(shù);yij為第i行第j列的DTY質(zhì)量指標(biāo)值。

        DTY碳排放量(Ce)越小越好,因此采用望小特性的信噪比進(jìn)行數(shù)據(jù)轉(zhuǎn)換如式(4)所示。

        (4)

        式中:bij為轉(zhuǎn)換后的第i行第j列的DTY碳排放量信噪比值;n為每組試驗(yàn)的重復(fù)次數(shù);yij為第i行第j列的DTY碳排放量值。

        由此,具體的正交試驗(yàn)方案設(shè)計(jì)如下。選擇油輪轉(zhuǎn)速(Vo)、加工速度(Vl)、牽伸比(B)、油尺高度(Ho)、第一熱箱溫度(T1)及第二熱箱溫度(T2)作為試驗(yàn)優(yōu)化指標(biāo),其他參數(shù)為設(shè)定標(biāo)準(zhǔn)值。反映DTY質(zhì)量的斷裂強(qiáng)度(Z1)、斷裂伸長率(Z2)、卷曲收縮率(Z3)、沸水收縮率(Z4)、卷曲穩(wěn)定度(Z5)以及碳排放量(Ce)作為綜合工藝優(yōu)化,設(shè)計(jì)6因素5水平的正交試驗(yàn)方案L25(56),試驗(yàn)因素及水平如表4所示。

        2.2 基于模糊多準(zhǔn)則的工藝參數(shù)優(yōu)化

        2.2.1 數(shù)據(jù)預(yù)處理

        為消除因數(shù)據(jù)量綱不同帶來的影響,還需要對(duì)轉(zhuǎn)換后的DTY指標(biāo)信噪比值進(jìn)行規(guī)范化處理,具體的處理過程如下。

        表4 試驗(yàn)因素及水平

        轉(zhuǎn)換后的DTY斷裂強(qiáng)度(Z1)、斷裂伸長率(Z2)、卷曲收縮率(Z3)、沸水收縮率(Z4)、卷曲穩(wěn)定度(Z5)的信噪比值處理如式(5)所示。

        (5)

        式中:Aij為aij經(jīng)規(guī)范化處理后第i行第j列的結(jié)果值。

        轉(zhuǎn)換后的DTY碳排放量(Ce)的信噪比值的計(jì)算公式如式(6)所示。

        (6)

        式中:Bij為bij經(jīng)規(guī)范化處理后第i行第j列的結(jié)果值。

        與此同時(shí),將Aij與Bij進(jìn)行數(shù)據(jù)組合,形成DTY指標(biāo)數(shù)據(jù)規(guī)范化矩陣G,并令Gij為G的第i行第j列元素。

        2.2.2 滌綸低彈絲質(zhì)量指標(biāo)組合賦權(quán)

        在數(shù)據(jù)處理的基礎(chǔ)上,為充分優(yōu)化工藝參數(shù)之間的強(qiáng)關(guān)聯(lián)關(guān)系,結(jié)合熵值賦權(quán)法與客觀權(quán)重賦值法CRITIC賦權(quán)法,形成一種適合滌綸低彈絲工藝參數(shù)優(yōu)化的組合賦權(quán)法,以待通過熵值賦權(quán)法,彌補(bǔ)CRITIC賦權(quán)法在數(shù)據(jù)離散性方面不足的問題。其中:用ej表示熵值賦權(quán)法中第j個(gè)指標(biāo)的權(quán)重值,如式(7)所示;用dj表示CRITIC賦權(quán)法中第j個(gè)指標(biāo)的權(quán)重值,如式(8)所示。

        (7)

        (8)

        式中:σj為標(biāo)準(zhǔn)差;rij為指標(biāo)之間的相關(guān)系數(shù)。

        在式(8)基礎(chǔ)上,利用乘除法將dj與ej進(jìn)行組合計(jì)算,形成滌綸低彈絲質(zhì)量指標(biāo)的組合權(quán)重kj,如式(9)所示。

        (9)

        2.2.3 計(jì)算各試驗(yàn)方案的利益比率

        (10)

        (11)

        2.2.3.2計(jì)算群體效益值 令Si為滌綸低彈絲各工藝參數(shù)試驗(yàn)方案中的群體效益值,則具體的計(jì)算過程如式(12)所示。

        (12)

        2.2.3.3計(jì)算個(gè)體遺憾值 令Wi為滌綸低彈絲各工藝參數(shù)試驗(yàn)方案中的個(gè)體遺憾值,具體計(jì)算過程如式(13)所示。令Qi為滌綸低彈絲各工藝參數(shù)試驗(yàn)方案中的利益比率,對(duì)應(yīng)的計(jì)算過程如式(14)所示。

        (13)

        (14)

        式中:S+、S-分別表示試驗(yàn)方案的群體最大、最小效益值;W+、W-分別表示決策方案的最大、最小個(gè)體遺憾值;v表示決策機(jī)制系數(shù),用于平衡Si與Wi的影響差異,取0.5。

        3 試驗(yàn)驗(yàn)證

        3.1 基于信噪比的正交試驗(yàn)數(shù)據(jù)規(guī)范化

        根據(jù)表4中的試驗(yàn)因素及水平開展正交試驗(yàn),得到DTY的斷裂強(qiáng)度(Z1)、斷裂伸長率(Z2)、卷曲收縮率(Z3)、沸水收縮率(Z4)、卷曲穩(wěn)定度(Z5),并核算各試驗(yàn)方案的碳排放量(Ce)。將具有望大特性的DTY質(zhì)量指標(biāo)值代入式(3),并將具有望小特性的碳排放量(Ce)代入式(4),經(jīng)計(jì)算得到的信噪比值如表5所示。

        利用VIKOR方法在處理綜合評(píng)價(jià)問題方面的優(yōu)勢[22-23],將多工藝參數(shù)強(qiáng)耦合問題通過矩陣歸一確定正理想解和負(fù)理想解,計(jì)算每組試驗(yàn)方案與理想解之間的差距,并比較與理想解的接近度。根據(jù)表5中的數(shù)據(jù),將DTY質(zhì)量指標(biāo)轉(zhuǎn)化后的信噪比值代入式(5)進(jìn)行處理,并將碳排放量(Ce)轉(zhuǎn)化后的信噪比值代入式(6)進(jìn)行處理,得到的規(guī)范化處理后結(jié)果如表6所示。

        表5 各試驗(yàn)方案指標(biāo)數(shù)據(jù)計(jì)算后的信噪比值

        根據(jù)2.2.2節(jié)中DTY指標(biāo)組合賦權(quán)方法,通過式(7)~(9)的數(shù)據(jù)處理,最終得到滌綸低彈絲生產(chǎn)過程中碳排放量及質(zhì)量指標(biāo)的組合權(quán)重值,如表7所示。

        3.2 試驗(yàn)方案的利益比率確定

        在得到試驗(yàn)數(shù)據(jù)規(guī)范化處理結(jié)果的基礎(chǔ)上,結(jié)合式(12)~(14),分別對(duì)每組數(shù)據(jù)進(jìn)行計(jì)算,得到每組試驗(yàn)對(duì)應(yīng)的群體效益值Si、個(gè)體遺憾值Wi和決策方案的利益比率Qi,如表8所示。

        表6 信噪比值規(guī)范化

        表7 滌綸低彈絲指標(biāo)組合權(quán)重

        表8 每組試驗(yàn)的Si、Wi、Qi值

        在此基礎(chǔ)上,利用影響DTY成品品質(zhì)的5個(gè)關(guān)鍵工藝參數(shù)水平對(duì)應(yīng)的Qi值對(duì)決策方案進(jìn)行評(píng)價(jià),主要包括油輪轉(zhuǎn)速(Vo)、加工速度(Vl)、牽伸比(R)、油尺高度(Ho)、第一熱箱溫度(T1)及第二熱箱溫度(T2),經(jīng)計(jì)算得到DTY關(guān)鍵工藝參數(shù)平均Qi值(即決策方案的利益比率)如表9所示。其中,平均Qi值越小越接近正理想解,則說明工藝參數(shù)水平越優(yōu)。

        由表9可知,對(duì)工藝參數(shù)加工速度、牽伸比、第一熱箱溫度、第二熱箱溫度、油輪轉(zhuǎn)速、油尺高度而言,其Qi值之間分別存在如下關(guān)系,即:

        表9 工藝參數(shù)平均Qi值

        1)對(duì)于加工速度,由于Q3

        2)對(duì)于牽伸比,由于Q2

        3)對(duì)于第一熱箱溫度,由于Q3

        4)對(duì)于第二熱箱溫度,由于Q2

        5)對(duì)于油輪轉(zhuǎn)速,由于Q4

        6)對(duì)于油尺高度,由于Q1

        分別對(duì)應(yīng)最優(yōu)工藝水平下的滌綸低彈絲關(guān)鍵工藝參數(shù)值,可見,DTY關(guān)鍵工藝參數(shù)的最優(yōu)組合為:加工速度600 m/min、牽伸比1.55、第一熱箱溫度195 ℃、第二熱箱溫度150 ℃、油輪轉(zhuǎn)速0.7 r/min、油尺高度180 mm。

        由表9可知,在不同工藝參數(shù)水平下,熱箱溫度、油輪轉(zhuǎn)速及油尺高度的平均Qi值的變化最為顯著,說明三者的交互作用對(duì)優(yōu)化綜合目標(biāo)的影響程度也最為顯著。在此基礎(chǔ)上,為進(jìn)一步驗(yàn)證工藝參數(shù)之間的優(yōu)化效果,還需對(duì)關(guān)鍵工藝參數(shù)之間的強(qiáng)耦合性關(guān)系進(jìn)行殘差值分析。為此,借助MatLab R2018b對(duì)油輪轉(zhuǎn)速、加工速度、牽伸比、油尺高度、第一熱箱溫度、第二熱箱溫度之間的耦合關(guān)系進(jìn)行擬合,得到的二次回歸方程殘差分布圖如圖2所示。

        圖2 關(guān)鍵工藝參數(shù)的殘差分布圖

        由圖2可知,關(guān)鍵工藝參數(shù)的殘差值的正態(tài)分布概率逼近平衡線,且分布均勻,說明工藝參數(shù)之間的優(yōu)化效果良好。

        3.3 工藝優(yōu)化前后滌綸低彈絲性能對(duì)比

        以83 dtex(36 f)低彈絲加彈工藝為例,結(jié)合表1中所示的FK6V-700型拉伸變形機(jī)加工工藝參數(shù),對(duì)比分析工藝優(yōu)化前和優(yōu)化后滌綸低彈絲生產(chǎn)工藝參數(shù)的變化,具體結(jié)果如表10所示??芍?,在83 dtex(36 f)低彈絲的加彈工藝中,熱箱溫度、油輪轉(zhuǎn)速及油尺高度3個(gè)工藝參數(shù)的強(qiáng)耦合性得到了較好的解耦。經(jīng)工藝優(yōu)化后,滌綸低彈絲的質(zhì)量指標(biāo)變化情況,如表11所示。

        表10 工藝參數(shù)優(yōu)化前后對(duì)比值

        表11 質(zhì)量指標(biāo)優(yōu)化前后對(duì)比值

        由表11可知,較優(yōu)化前,優(yōu)化后83 dtex(36 f)低彈絲的平均斷裂強(qiáng)度達(dá)到了3.75 cN/dtex,提高了15.84%;平均斷裂伸長率達(dá)到了26.49%,提高了4.04%;平均卷曲收縮率以及平均卷曲穩(wěn)定度分別提高了15.53%和1.24%,后加工過程中不易產(chǎn)生毛絲與斷頭現(xiàn)象。然而經(jīng)工藝優(yōu)化后,83 dtex(36 f)低彈絲的平均沸水收縮率有所下降,這其中的主要原因在于加工速度的下降,減少了絲束的熱定型時(shí)間,從而增加了纖維的內(nèi)應(yīng)力松弛度,造成了絲條在沸水中的應(yīng)力松弛變小,從而使得沸水收縮率下降。

        在此基礎(chǔ)上,將表10中優(yōu)化后的工藝參數(shù)代入公式(2),在滌綸低彈絲關(guān)鍵工藝參數(shù)最優(yōu)組合下得到的滌綸低彈絲生產(chǎn)過程中的碳排放量為191.534 6 kg(CO2),與優(yōu)化前工藝條件下的碳排放量相比,其下降了4.58%。

        4 結(jié) 論

        針對(duì)滌綸低彈絲生產(chǎn)工藝參數(shù)間存在強(qiáng)耦合的問題,在明晰生產(chǎn)流程的基礎(chǔ)上,對(duì)生產(chǎn)過程的碳排放進(jìn)行了分析,并對(duì)滌綸低彈絲工藝參數(shù)之間關(guān)系進(jìn)行了表示;進(jìn)而設(shè)計(jì)了基于信噪比的工藝參數(shù)關(guān)系正交試驗(yàn)方案,開展了基于組合賦權(quán)VIKOR方法的工藝參數(shù)優(yōu)化。

        經(jīng)試驗(yàn)驗(yàn)證表明,當(dāng)工藝優(yōu)化參數(shù)組合為油輪轉(zhuǎn)速0.7 r/min,加工速度600 m/min,牽伸比1.55,油尺高度180 mm,第一熱箱溫度195 ℃,第二熱箱溫度150 ℃時(shí),可降低DTY生產(chǎn)過程中的碳排放量,同時(shí)提高滌綸低彈絲產(chǎn)品質(zhì)量特征值。工藝參數(shù)優(yōu)化后,滌綸低彈絲的平均斷裂強(qiáng)度提高了15.84%、平均斷裂伸長率提高了4.04%、碳排放量降低了4.58%,充分說明構(gòu)建的基于綜合賦權(quán)的工藝參數(shù)優(yōu)化模型有利于解決工藝參數(shù)之間的強(qiáng)耦合問題。

        猜你喜歡
        收縮率滌綸賦權(quán)
        滌綸短纖維干熱收縮率測量不確定度分析
        影響滌綸短纖維干熱收縮率因素探討
        論鄉(xiāng)村治理的有效賦權(quán)——以A縣扶貧項(xiàng)目為例
        中國西部(2022年2期)2022-05-23 13:28:20
        聚丙烯材料收縮率影響因素研究
        上海塑料(2021年3期)2022-01-06 14:05:02
        企業(yè)數(shù)據(jù)賦權(quán)保護(hù)的反思與求解
        試論新媒體賦權(quán)
        活力(2019年15期)2019-09-25 07:22:12
        基于改進(jìn)AHP熵博弈賦權(quán)的輸變電工程評(píng)價(jià)
        母線生產(chǎn)中端面線收縮率與寬厚比間關(guān)系的統(tǒng)計(jì)規(guī)律
        山東冶金(2018年5期)2018-11-22 05:12:20
        茶碳滌綸網(wǎng)布填補(bǔ)國內(nèi)空白
        11 dtex/6 f滌綸超細(xì)低彈絲的研制
        国产精品一区二区久久久av | 中文字幕人妻熟女人妻| 50岁熟妇的呻吟声对白| 精品国产av无码一道| 一区二区三区岛国av毛片 | 国产av国片精品有毛| 五月婷婷俺也去开心| 欧美性xxx久久| 99精品人妻少妇一区二区三区| 日本成本人片视频免费| 国产成人综合在线视频| 97久久综合区小说区图片专区| 中文字幕一区乱码在线观看| 成人国成人国产suv| 比比资源先锋影音网| 国产精品va在线观看一| 亚洲岛国一区二区三区| 永久免费人禽av在线观看| 中文字幕av日韩精品一区二区| 国内精品九九久久精品小草| 最新中文字幕日韩精品| 中文字幕无码毛片免费看| 在线欧美不卡| 亚洲不卡毛片在线观看| 精品无码一区二区三区的天堂| 无遮无挡爽爽免费视频| 亚州AV成人无码久久精品| 日本频道一区二区三区| 午夜免费福利小电影| 国产精品视频一区国模私拍| 女同性恋亚洲一区二区| 在线视频国产91自拍| 国产熟女露脸大叫高潮| 久久精品国产亚洲av大全相关| 精品老熟女一区二区三区在线| 亚洲av成人精品日韩在线播放| 久久网视频中文字幕综合| 天堂av中文在线官网| 无码精品人妻一区二区三区漫画| a亚洲va欧美va国产综合| 精品999无码在线观看|