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

    凝聚槽類型對轉(zhuǎn)杯內(nèi)氣流場影響的數(shù)值模擬

    2017-05-24 14:44:16楊瑞華高衛(wèi)東生態(tài)紡織教育部重點實驗室江南大學江蘇無錫214122
    紡織學報 2017年5期
    關(guān)鍵詞:成紗槽內(nèi)紡紗

    劉 超, 楊瑞華, 薛 元, 高衛(wèi)東(生態(tài)紡織教育部重點實驗室(江南大學), 江蘇 無錫 214122)

    凝聚槽類型對轉(zhuǎn)杯內(nèi)氣流場影響的數(shù)值模擬

    劉 超, 楊瑞華, 薛 元, 高衛(wèi)東
    (生態(tài)紡織教育部重點實驗室(江南大學), 江蘇 無錫 214122)

    為探究G、T、U、S型凝聚槽對36 mm直徑轉(zhuǎn)杯內(nèi)氣流場的影響,在軟件Fluent中對三維流場進行數(shù)值計算分析。結(jié)果表明:在相同工藝條件下,凝聚槽在1周(0°~360°) 范圍內(nèi),4種槽型的速度大小為G型>T型>U型>S型,在0°與360°位置處靜壓大小為G型>S型>U型>T型,其余各角度位置處靜壓大小為S型>U型>T型>G型。各槽型內(nèi)氣流靜壓和速度趨勢基本一致,以T型槽為例,輸棉通道內(nèi)的靜壓位于-32 886.15~18 224.56 Pa之間,轉(zhuǎn)杯內(nèi)的靜壓大部分處于-13 719.63~-7 330.80 Pa之間;輸棉通道內(nèi)的氣流隨著管道直徑的減小而加速運動,在出口處達到最大值261.81 m/s。

    轉(zhuǎn)杯紡; 凝聚槽; 數(shù)值模擬; 氣流場; 速度; 靜壓

    轉(zhuǎn)杯紡是一種自由端紡紗方法,成紗過程中自由態(tài)纖維由離心力與氣流力作用在凝聚槽凝聚,隨高速旋轉(zhuǎn)的轉(zhuǎn)杯加捻成紗。轉(zhuǎn)杯紡紗機按照轉(zhuǎn)杯型號的不同,可分為自排風式和抽氣式,二者成紗過程基本一致,但形成負壓的方式不同[1]。

    轉(zhuǎn)杯型號是以凝聚槽的形狀進行分類的,可分為G型、T型、S型、U型等。G型槽適紡中、低線密度機織、針織、牛仔布用紗以及染色布用紗,適紡線密度小于100 tex;T型槽應用面廣,適紡線密度小于60 tex;S型槽適紡高含雜棉及各類化纖,適紡線密度大于28 tex;U型槽適紡粗支紗,紡出的紗線強度較低,波紋較少,適紡線密度大于35 tex[1-2]。成紗質(zhì)量與凝聚槽的槽型、凝聚角度的大小和摩擦因數(shù)等因素有很大關(guān)系。

    曾泳春等[3]采用二維模型對噴氣織機噴嘴中的氣流運動進行數(shù)值計算,結(jié)果表明,噴孔傾角為45°時,噴嘴內(nèi)氣流速度分布均勻,成紗性能好;鄒專勇等[4]建立了噴氣渦流紡噴嘴結(jié)構(gòu)三維模型,表征了噴嘴內(nèi)氣流場的流動情況,分析了流場內(nèi)的速度分布、壓力分布以及流場與成紗機制之間的聯(lián)系;梁海順等[5]應用NUMECA軟件模擬了噴氣織機主噴嘴內(nèi)部流場特征;林惠婷等[6]通過數(shù)值計算方法研究了纖維輸送管和轉(zhuǎn)杯內(nèi)的氣流運動特性,探討了纖維輸送管特征數(shù)和空間位置角對纖維伸直狀態(tài)的影響。

    本文以轉(zhuǎn)杯紡紗通道內(nèi)氣流場作為研究對象,應用Fluent對轉(zhuǎn)杯紡紗通道內(nèi)的氣流流動特征進行模擬,并分析了4種不同型號的凝聚槽對轉(zhuǎn)杯內(nèi)氣流場的影響,為研究凝聚槽型號對紗線質(zhì)量的影響提供理論依據(jù)。

    1 抽氣式轉(zhuǎn)杯紡紗通道流場三維模型

    圖1示出紡紗通道內(nèi)流場的三維模型,包括引紗管、輸棉通道、轉(zhuǎn)杯等。此模型是根據(jù)日發(fā)機械股份有限公司生產(chǎn)的某種型號轉(zhuǎn)杯紡紗機內(nèi)的轉(zhuǎn)杯所建立,轉(zhuǎn)杯直徑為36 mm。從圖1(b)可看出,轉(zhuǎn)杯直徑最大位置處形成一個凝聚纖維的凹槽,即為凝聚槽。為探究凝聚槽類型對轉(zhuǎn)杯內(nèi)氣流場的影響,對直徑為36 mm,槽型為G型、T型、S型、U型的轉(zhuǎn)杯分別進行模擬。

    圖1 紡紗通道內(nèi)流場三維模型Fig.1 3-D Airflow field model of spinning channel (a) and cross sectional view of YZ (b)

    2 紡紗通道內(nèi)流場的數(shù)值計算

    2.1 物理模型

    由于氣流在紡紗通道內(nèi)高速流動,因此,紡紗通道內(nèi)氣流的流動是可壓縮的黏性流動。為簡化數(shù)值計算,忽略紡紗通道流場內(nèi)的纖維運動和熱傳遞?;究刂品匠虨槟芰渴睾惴匠毯蛣恿渴睾惴匠獭?/p>

    能量守恒方程的微分表達式[7]為

    (1)

    式中:ρ為氣體密度;h為熵;ui為氣流速度在xi方向的分量;k為分子傳導率;kt為湍流傳遞引起的傳導率;T為溫度;Sh為被定義的體積源[7]。

    動量守恒方程的微分表達式[7]為

    (2)

    式中:uj為氣流速度在xj方向的分量;p為壓強;τij為應力矢量;gi為沿著i方向的重力分量;Fi為阻力和能源引起的其他源項[7]。

    湍流模型采用標準的k-ε模型方程:

    式中:Gk為層流速度梯度產(chǎn)生的湍動能;Gb為浮力產(chǎn)生的湍動能;YM為可壓縮湍流中脈動膨脹的貢獻;μt為湍流黏性系數(shù);σk為湍動能k對應的普朗特數(shù);σε為耗散能ε對應的普朗特數(shù);Sk和Sε為用戶定義的源項;C1ε=1.42,C2ε=1.68,C3ε=0.09[7]。

    2.2 邊界條件設定

    入口邊界:設置輸棉通道入口為質(zhì)量流量入口,入口質(zhì)量流量=氣體密度(1.29 kg/m3)×入口速度(50 m/s)×入口截面積(83.72 mm2)=0.005 4 kg/s。

    出口邊界:由轉(zhuǎn)杯紡紗工藝可知,出口邊界為壓力出口,出口處的壓強值設為-8 kPa[8]。

    固體壁面:在紡紗過程中,轉(zhuǎn)杯沿著順時針方向高速旋轉(zhuǎn),因此,轉(zhuǎn)杯壁面的運動類型為旋轉(zhuǎn)[9],轉(zhuǎn)速設為12萬r/min。

    3 結(jié)果與討論

    3.1 紡紗通道內(nèi)流場分析(T型)

    在對模擬計算結(jié)果進行后處理的過程中發(fā)現(xiàn),G型、T型、S型和U型槽的紡紗通道內(nèi)三維流場的氣流速度矢量圖及壓強分布圖的趨勢基本保持一致。本文將T型槽的紡紗通道內(nèi)三維流場作為例子分析其流場模擬的結(jié)果。

    圖2示出T型槽靜壓分布圖。由圖2(a)可知,轉(zhuǎn)杯紡紗通道內(nèi)的靜壓處于-13.719 63~-941.96 Pa之間,小于輸棉通道內(nèi)的靜壓。由圖2(b)和圖2(c)可知,輸棉通道入口處壓力達到最大值,為18 224.56 Pa,當纖維由輸送管入口處向出口處運動的過程中,輸棉通道內(nèi)的靜壓逐漸變成了負壓,出口處的靜壓值為-32 886.10 Pa,纖維輸送管內(nèi)的負壓有利于纖維順利脫離分梳輥,向轉(zhuǎn)杯凝聚槽內(nèi)運動。從圖2(c)還可看出:轉(zhuǎn)杯內(nèi)的大部分靜壓為負壓,處于-13 719.63~-7 330.80 Pa之間,有利于纖維沿著滑移面向凝聚槽運動,使紡紗順利進行;轉(zhuǎn)杯內(nèi)的靜壓呈現(xiàn)不均勻的分布,使得轉(zhuǎn)杯由于受到氣流不平衡力的影響,在旋轉(zhuǎn)過程中的穩(wěn)定性降低,磨損程度增大,使用壽命大大縮短。

    圖2 T型槽靜壓分布圖Fig.2 T-type static pressure in spinning channel (a), transfer channel (b) and rotor (c)

    圖3示出T型槽速度分布失量圖。由圖3(a)紡紗通道內(nèi)速度分布矢量圖和圖3(b)輸棉通道內(nèi)速度分布矢量圖可知,隨著輸棉通道的橫截面積逐漸縮小,輸棉通道中的氣流速度值呈梯度遞增,出口處的速度值最大,達到261.81 m/s,氣流在輸棉通道中的加速運動既有利于纖維順利向轉(zhuǎn)杯凝聚槽內(nèi)移動,更有利于輸棉通道中彎鉤纖維的伸直。由圖3(c)轉(zhuǎn)杯內(nèi)速度分布矢量圖可知:輸棉通道出口處的氣流碰撞到轉(zhuǎn)杯滑移面后,向周圍發(fā)散運動;隨著轉(zhuǎn)杯的高速旋轉(zhuǎn),由于受到杯內(nèi)負壓和滑移面摩擦阻力的影響,分成2股不同大小和方向的氣流,一股繞著轉(zhuǎn)杯順時針向上運動,一股繞著轉(zhuǎn)杯逆時針向上運動,在2股氣流交匯處形成渦流,渦流的形成是纖維產(chǎn)生包纏加捻的關(guān)鍵,有助于成紗。

    3.2 凝聚槽類型對轉(zhuǎn)杯內(nèi)氣流場的影響

    在高速旋轉(zhuǎn)的轉(zhuǎn)杯離心力作用下,纖維由滑移面向凝聚槽內(nèi)移動,疊合成環(huán)形的凝聚須條。凝聚的目的是為了連續(xù)紡紗,這是轉(zhuǎn)杯紡紗的必要工藝過程,因此,本文對G型、T型、S型和U型槽內(nèi)的氣流場進行數(shù)值模擬計算,結(jié)果如圖4~5所示。

    圖4示出G型、T型、S型和U型槽的壁面靜壓分布云圖。由圖可知:在相同的工藝條件下,G型槽內(nèi)靜壓處于-43 468.81~22 059.51 Pa之間,T型槽內(nèi)靜壓處于-32 886.15~18 224.56 Pa之間,S型槽內(nèi)靜壓處于-31 944.15~19 432.92 Pa之間,U型槽內(nèi)靜壓處于-32 269.47~19 154.20 Pa之間;在4種槽型中,G型槽高壓值最大,T型槽高壓值最小,S型槽低壓值最大,G型槽低壓值最小。

    圖3 T型槽速度分布矢量圖Fig.3 T-type airflow velocity in spinning channel (a), transfer channel (b) and rotor (c)

    注:單位為Pa。圖4 不同類型凝聚槽壁面靜壓分布Fig.4 Static pressure in different grooves. (a) G type; (b) T type; (c) S type; (d) U type

    圖5示出不同類型凝聚槽壁面的速度分布,圖6示出轉(zhuǎn)杯凝聚槽1周0°~360°各個角度處的壓力和速度分布規(guī)律。0°位置處為輸棉通道中心線延長至出口與凝聚槽交匯處,角度沿著順時針方向遞增。由圖5、6可知,在相同工藝條件下,G型槽內(nèi)速度值最大,T型槽與U型槽次之,且T型槽內(nèi)速度值高于U型槽,S型槽內(nèi)速度值最小。

    注:單位為Pa。圖5 不同類型凝聚槽壁面速度分布Fig.5 Velocity of different types of grooves. (a) G type; (b) T type; (c) S type; (d) U type

    圖6 凝聚槽中壓力和速度分布規(guī)律Fig.6 Pressure distribution of grooves (a) and velocity distribution of grooves (b)

    從圖6(a)可看出,凝聚槽的靜壓在0°位置處達到最大值,隨著轉(zhuǎn)杯的高速旋轉(zhuǎn),靜壓值急速下降,且在0°~20° 區(qū)間內(nèi),靜壓值由正變?yōu)樨摗5?20°~360°區(qū)間內(nèi),靜壓值又由負變?yōu)檎?,?0°~320°區(qū)間內(nèi),壓力相差不大,靜壓比較穩(wěn)定。從圖6(b)可看出,在0°位置處,氣流速度較低,隨著角度遞增,速度急劇增加,在0°~20°區(qū)間內(nèi),速度梯度較大,在20°位置處,速度達到最大值,這是由于輸棉通道出口處的纖維沿著轉(zhuǎn)杯滑移面向凝聚槽中運動的過程中,纖維在該位置處的聚合效果好。沿著順時針方向,氣流速度在凝聚槽120°~180°區(qū)間內(nèi)又減少了1個數(shù)量級,結(jié)合圖5可看出,輸棉通道出口處流出的高速氣流,碰到轉(zhuǎn)杯滑移面后,在滑移面摩擦阻力和負壓作用下,形成2股不同大小、反向運動的氣流,一股沿著轉(zhuǎn)杯順時針方向旋轉(zhuǎn),一股沿著轉(zhuǎn)杯逆時針方向旋轉(zhuǎn),2股氣流在180°位置處相遇,所以此處速度值和壓力值較小,不利于纖維的聚合、加捻成紗。

    由圖6可知:4種槽型在0°和360°位置處靜壓大小為G型>S型>U型>T型,其余角度位置處靜壓大小為S型>U型>T型>G型;4種槽型在各角度位置處的速度大小為G型>T型>U型>S型。由氣流紡回轉(zhuǎn)紗條張力的測試與分析[10]結(jié)果可知,速度越大,凝聚槽內(nèi)的氣流由于受到轉(zhuǎn)杯內(nèi)離心力的作用,越來越易聚集,對纖維須條產(chǎn)生的抱合力較大,纖維與纖維間的緊密度增加,使得纖維須條的加捻阻力減小,導致紗線的捻回傳遞長度減小,紡出的紗線強力相應增加,因此,G型槽紡成的紗線強力最大,大于T型槽和U型槽,而S型凝聚槽中的纖維須條呈現(xiàn)蓬松狀態(tài),使得須條的加捻阻力增大,導致紗線的捻回傳遞長度變大,纖維纏繞增多,成紗強力變小。這又和轉(zhuǎn)杯紡捻度傳遞長度的解析研究[11]中的結(jié)論一致,因為凝聚槽的夾角越大,紡出紗線的捻度傳遞長度就越大,G型凝聚槽的夾角是35°,T型凝聚槽的夾角是45°,U型凝聚槽的夾角是80°,S型凝聚槽的夾角是85°,所以4種槽型紡出的紗線的捻回傳遞長度的大小為S型>U型>T型>G型。捻度傳遞長度纏繞纖維與成紗質(zhì)量[12]的實驗研究表明,紗線的捻回傳遞長度與成紗質(zhì)量的關(guān)系密切。當捻回傳遞長度增加,紗線的強力下降;當捻回傳遞長度減小,紗線的強力增大。4種槽型的成紗強力大小為G型>T型>U型>S型。

    4 結(jié) 論

    采用3-D模型,借助流體計算分析軟件對紡紗通道內(nèi)的氣流場進行模擬計算,并研究了G型、T型、U型和S型槽對轉(zhuǎn)杯內(nèi)氣流流動特征的影響,得到以下結(jié)論。

    1)4種槽型的紡紗通道內(nèi)氣流場靜壓分布云圖與速度分布矢量圖的趨勢基本一致,以T型槽為例,轉(zhuǎn)杯紡紗通道內(nèi)的靜壓低于纖維輸送管內(nèi)的靜壓,纖維在纖維輸送管入口處靜壓值最大為18 224.56 Pa,出口處靜壓值最小為-32 886.10 Pa;轉(zhuǎn)杯內(nèi)的靜壓分布不均勻,且大部分靜壓為負壓。

    2)隨著纖維輸送管管道直徑逐漸減小,纖維由輸送管進入轉(zhuǎn)杯內(nèi)的過程中加速運動,且在纖維輸送管出口處速度達到最大值261.81 m/s,有利于纖維的定向伸直;纖維到達轉(zhuǎn)杯滑移面上時一方面向凝聚槽運動,一方面沿著轉(zhuǎn)杯旋轉(zhuǎn)方向運動。

    3)G型槽內(nèi)靜壓處于-43 468.81~22 059.51 Pa之間,速度處于0~296.27 m/s之間;T型槽內(nèi)靜壓處于-32 886.15~18 224.56 Pa之間,速度處于0~261.81 m/s之間;S型槽內(nèi)靜壓處于-31 944.15~19 432.92 Pa之間,速度處于0~253.31 m/s之間;U型槽內(nèi)靜壓處于-32 269.47~19 154.20 Pa之間,速度處于0~261.20 m/s之間。

    4)以輸棉通道中心線延長至出口與凝聚槽交匯處為0°,角度沿凝聚槽順時針方向遞增至360°;4種槽型在0°與360°位置處靜壓大小為G型>S型>U型>T型,其余角度位置處靜壓大小為S型>U型>T型>G型;4種槽型的速度大小為G型>T型>U型>S型;G型槽紡成的纖維抱合力大,纖維須條緊密,須條加捻阻力小,捻度易于傳遞,所紡紗線強度高于其他3種槽型。

    FZXB

    [1] 張百祥,周慈念.轉(zhuǎn)杯紡紗[M].北京:紡織工業(yè)出版社, 1993:61-65. ZHANG Baixiang,ZHOU Cinian.Rotor Spinning [M]. Beijing:Textile Industry Press,1993:61-65.

    [2] 葉鴻璣,徐潼.轉(zhuǎn)杯紡紗[M].濟南:濟南出版社,1989:8-13. YE Hongji,XU Tong.Rotor Spinning [M].Jinan:Jinan Press,1989:8-13.

    [3] 曾泳春,郁崇文.噴氣紡噴嘴中氣流流動的數(shù)值計算[J].東華大學學報(自然科學版),2002(5):11-16. ZENG Yongchun, YU Chongwen.Value of airflow jet spinning nozzle computing [J].Journal of Donghua University (Natural Science),2002(5):11-16.

    [4] 鄒專勇,俞建勇,薛文良,等.噴氣渦流紡噴嘴內(nèi)部三維流場的數(shù)值研究[J].紡織學報,2008,29(2):86-89. ZOU Zhuanyong, YU Jianyong, XUE Wenliang,et al. Numerical study of jet vortex spinning nozzle inside the three-dimensional flow field [J].Journal of Textile Research,2008,29(2):86-89.

    [5] 梁海順,楊昆,王貫超,等.基于NUMECA技術(shù)的噴氣織機主噴嘴內(nèi)部流場數(shù)值模擬[J].紡織器材,2008(3):12-16. LIANG Haishun, YANG Kun, WANG Guanchao, et al. Main nozzle internal flow field numerical simulation based on NUMECA technology[J].Textile Accessories, 2008(3):12-16.

    [6] 林惠婷,汪軍,曾泳春.輸棉通道幾何參數(shù)對轉(zhuǎn)杯紡氣流場影響的數(shù)值研究[J].紡織學報,2015, 36(2):98-104. LIN Huiting, WANG Jun,ZENG Yongchun.Numerical study cotton spinning channel geometry on the flow field [J]. Journal of Textile Research, 2015,36(2):98-104.

    [7] 王福軍.計算流體動力學分析[M].北京:清華大學出版社,2004:7-9. WANG Fujun.Computational Fluid Dynamics Analy-sis [M].Beijing: Tsinghua University Press,2004:7-9.

    [8] 劉超,楊瑞華,王鴻博,等.轉(zhuǎn)杯紡紗通道三維流場的數(shù)值模擬[J].紡織學報,2016,37(9):145-150. LIU Chao,YANG Ruihua,WANG Hongbo, et al.Numerical simulation of three-dimensional flow field in rotor spinning passage [J].Journal of Textile Research, 2016,37(9):145-150.

    [9] 朱紅均,林遠華,謝龍漢.FLUENT流體分析及仿真實用教程[M].北京:人民郵電出版社,2010:237-238. ZHU Hongjun,LIN Yuanhua,XIE Longhan.FLUENT Fluid Analysis and Simulation Practical Tutorial [M].Beijing:The People′s Posts and Telecommunications Press, 2010:237-238.

    [10] 張百祥,秦洪奇,陳麗珍.氣流紡回轉(zhuǎn)紗條張力的測試與分析[J].東華大學學報(自然科學版),1981(1):8. ZHANG Baixiang, QIN Hongqi,CHEN Lizhen. Rotor spinning yarn tension test and analysis [J].Journal of Donghua University (Natural Science),1981(1): 8.

    [11] 汪軍,黃秀寶.轉(zhuǎn)杯紡捻度傳遞長度的解析研究[J].中國紡織大學學報,2000(1):64-69. WANG Jun, HUANG Xiubao. Analytical studies spinning twist length transfer[J]. Journal of China Textile University, 2000(1): 64-69.

    [12] 黃秀寶,梁金茹.捻度傳遞長度纏繞纖維與成紗質(zhì)量[J]. 紡織學報,1983,5(11):5-11. HUANG Xiubao,LIANG Jinru.Twist transfer length wound fiber and yarn quality[J].Journal of Textile Research,1983,5(11):5-11.

    Numerical simulation of influence of groove type on flow field inside rotor

    LIU Chao, YANG Ruihua, XUE Yuan, GAO Weidong
    (KeyLaboratoryofEco-Textiles(JiangnanUniversity),MinistryofEducation,Wuxi,Jiangsu214122,China)

    The influence of groove type on high speed airflow during rotor spun yarn was investigated. Airflow speed and static pressure in four types of groove, such as G, T, U and S of the 36 mm diameter rotor, were studied by Fluent software. The results show that under the same conditions, speeds in four groove size are G>T>U>S within the range from 0° to 360°. In 0° and 360° position, the static pressures are G>S>U>T. While for the rest of the angle position, the static pressures are S>U> T>G. Taking T groove as example, static pressures within the transfer channel are between -32 886.15 Pa and 18 224.56 Pa, static pressures within the rotor are between -13 719.63 Pa and -7 330.80 Pa. The airstream accelerates from the transfer channel inlet to the outlet with the decrease diameter of the pipe, and reaches the largest value to 261.81 m/s at the outlet.

    rotor spinning; groove; numerical simulation; airflow field; speed; pressure

    2016-03-23

    2017-02-06

    國家自然科學基金項目(51403085);中央高?;究蒲袠I(yè)務費專項資金資助項目(JUSRP51631A);江蘇省自然科學基金項目(BK20130148);江蘇省政策引導類計劃(產(chǎn)學研合作)-前瞻性聯(lián)合研究項目(BY2016022-29)

    劉超(1989—),女,碩士生。主要研究方向為轉(zhuǎn)杯紡復合紗成紗結(jié)構(gòu)與成紗氣流分析。楊瑞華,通信作者,E-mail: yangrh@jiangnan.edu.cn。

    10.13475/j.fzxb.20160304807

    TS 101.2

    A

    猜你喜歡
    成紗槽內(nèi)紡紗
    新型器材專件提高成紗質(zhì)量的應用體會
    紡織器材(2020年1期)2020-03-24 09:14:34
    渠首閘后防沖槽的拋石粒徑與消能試驗研究
    棉卷性能與成紗質(zhì)量關(guān)系模型的建立與驗證
    紡紗設備及技術(shù)的研究進展及發(fā)展前景
    迷你甜甜圈制作機
    KD405:多功能折疊筷子
    紡紗器材的技術(shù)進步與應用
    紡織導報(2015年6期)2015-12-31 19:29:30
    發(fā)熱纖維/棉/黏膠多組分噴氣紡紗線的開發(fā)
    新型紡紗工藝技術(shù)與設備的發(fā)展
    長絨棉纖維性能與成紗質(zhì)量關(guān)系的研究?
    在线观看一区二区三区| 国产成人aa在线观看| 国产精品久久视频播放| 制服丝袜大香蕉在线| 99在线人妻在线中文字幕| 美女免费视频网站| 亚洲 国产 在线| 老司机靠b影院| 老熟妇乱子伦视频在线观看| 中文字幕av在线有码专区| 久久热在线av| 少妇的丰满在线观看| 亚洲五月婷婷丁香| 床上黄色一级片| 国产精品,欧美在线| 亚洲精品一区av在线观看| 欧美三级亚洲精品| 久久久久国产精品人妻aⅴ院| 亚洲专区字幕在线| 国产成人精品久久二区二区免费| 国产探花在线观看一区二区| 观看免费一级毛片| 欧美日韩亚洲国产一区二区在线观看| www国产在线视频色| 欧美中文日本在线观看视频| 夜夜夜夜夜久久久久| 俺也久久电影网| 男女之事视频高清在线观看| 久久久久久久午夜电影| 一二三四社区在线视频社区8| 国产又色又爽无遮挡免费看| 国产精品香港三级国产av潘金莲| 久久久国产欧美日韩av| 午夜福利18| 国产男靠女视频免费网站| 无遮挡黄片免费观看| 午夜福利18| 一区福利在线观看| 国产欧美日韩一区二区三| 久久精品国产清高在天天线| 88av欧美| 国产精品久久久人人做人人爽| 99久久久亚洲精品蜜臀av| 窝窝影院91人妻| 亚洲中文日韩欧美视频| 操出白浆在线播放| 久久亚洲真实| 日本精品一区二区三区蜜桃| 男女午夜视频在线观看| 最好的美女福利视频网| 男女午夜视频在线观看| 国产三级中文精品| 午夜视频精品福利| 久久精品91蜜桃| 中文亚洲av片在线观看爽| 又紧又爽又黄一区二区| 国产视频一区二区在线看| 成人特级黄色片久久久久久久| 欧美日韩黄片免| 日韩中文字幕欧美一区二区| 麻豆成人av在线观看| 一边摸一边做爽爽视频免费| 亚洲中文字幕日韩| 亚洲专区中文字幕在线| 国产三级在线视频| 日韩精品免费视频一区二区三区| 麻豆国产97在线/欧美 | videosex国产| 亚洲男人天堂网一区| 久久精品影院6| 波多野结衣高清作品| 一二三四社区在线视频社区8| 长腿黑丝高跟| 视频区欧美日本亚洲| 国产精品精品国产色婷婷| 国产亚洲精品久久久久5区| 岛国在线免费视频观看| 亚洲av第一区精品v没综合| 国产黄色小视频在线观看| 国产精品自产拍在线观看55亚洲| 国产亚洲精品久久久久5区| 欧美在线一区亚洲| 欧美在线一区亚洲| 久久亚洲真实| 高清毛片免费观看视频网站| 国模一区二区三区四区视频 | 国产主播在线观看一区二区| 亚洲 国产 在线| 每晚都被弄得嗷嗷叫到高潮| 久久久久久九九精品二区国产 | 欧美乱妇无乱码| 亚洲人与动物交配视频| 成年免费大片在线观看| 国产成人欧美在线观看| 欧美中文综合在线视频| 国产成人啪精品午夜网站| 午夜久久久久精精品| 伦理电影免费视频| 精品国产乱子伦一区二区三区| 亚洲欧洲精品一区二区精品久久久| www.精华液| 成人18禁高潮啪啪吃奶动态图| 高潮久久久久久久久久久不卡| 午夜久久久久精精品| 日日摸夜夜添夜夜添小说| 色综合婷婷激情| 老熟妇乱子伦视频在线观看| 亚洲 欧美 日韩 在线 免费| 日韩欧美国产一区二区入口| 日本五十路高清| 色噜噜av男人的天堂激情| 亚洲av成人av| 两个人视频免费观看高清| 嫁个100分男人电影在线观看| 99在线视频只有这里精品首页| 女人被狂操c到高潮| xxx96com| 亚洲一区二区三区不卡视频| 亚洲国产精品合色在线| 全区人妻精品视频| 国产爱豆传媒在线观看 | 久久香蕉激情| 在线看三级毛片| 亚洲欧美精品综合久久99| 99热这里只有是精品50| 国产99久久九九免费精品| 一本大道久久a久久精品| 欧美三级亚洲精品| 村上凉子中文字幕在线| 一区福利在线观看| 亚洲精品美女久久av网站| 久久久精品欧美日韩精品| 亚洲国产精品成人综合色| 午夜精品在线福利| 亚洲免费av在线视频| 制服人妻中文乱码| 日日摸夜夜添夜夜添小说| 亚洲国产精品合色在线| 国产一区二区三区视频了| 在线视频色国产色| 国产精品久久久久久精品电影| 久久久久久久午夜电影| 日韩av在线大香蕉| 欧美成狂野欧美在线观看| 国产97色在线日韩免费| 国产成人啪精品午夜网站| 真人做人爱边吃奶动态| 曰老女人黄片| 动漫黄色视频在线观看| 久久久国产欧美日韩av| 19禁男女啪啪无遮挡网站| 亚洲性夜色夜夜综合| 精品久久久久久久末码| 国产av又大| 国内精品久久久久久久电影| 日本三级黄在线观看| 看免费av毛片| 99国产精品一区二区三区| 免费无遮挡裸体视频| 嫩草影视91久久| 两个人免费观看高清视频| 亚洲人成77777在线视频| 一个人观看的视频www高清免费观看 | 亚洲五月天丁香| 免费人成视频x8x8入口观看| 男插女下体视频免费在线播放| 一级作爱视频免费观看| 99久久精品国产亚洲精品| 国产高清videossex| 高潮久久久久久久久久久不卡| 午夜日韩欧美国产| 午夜激情av网站| 人妻夜夜爽99麻豆av| 日韩成人在线观看一区二区三区| 亚洲美女视频黄频| 日本精品一区二区三区蜜桃| 久久久久久大精品| 亚洲国产精品成人综合色| 成人手机av| 久久精品国产99精品国产亚洲性色| 中文在线观看免费www的网站 | 日韩精品青青久久久久久| 亚洲熟女毛片儿| 亚洲成av人片在线播放无| 欧美成人免费av一区二区三区| 在线观看免费午夜福利视频| 国产精品 国内视频| 久久久久久久午夜电影| 我要搜黄色片| 国产精品一及| 国产精品一及| 国产欧美日韩精品亚洲av| 99精品欧美一区二区三区四区| 丰满人妻一区二区三区视频av | 99在线人妻在线中文字幕| 午夜福利在线观看吧| 婷婷丁香在线五月| 婷婷精品国产亚洲av| 国产高清激情床上av| 一级作爱视频免费观看| 天天一区二区日本电影三级| 午夜福利免费观看在线| 伦理电影免费视频| 18禁黄网站禁片免费观看直播| 禁无遮挡网站| 欧美性长视频在线观看| 久久精品综合一区二区三区| 九色国产91popny在线| 亚洲乱码一区二区免费版| 一边摸一边抽搐一进一小说| 99国产精品一区二区三区| 国产高清视频在线观看网站| 国产v大片淫在线免费观看| 国产一区二区三区视频了| 午夜亚洲福利在线播放| 一区二区三区国产精品乱码| 五月玫瑰六月丁香| 成年免费大片在线观看| 久热爱精品视频在线9| 久久久久久大精品| 啦啦啦观看免费观看视频高清| 国产一区在线观看成人免费| 99热6这里只有精品| 国产精品亚洲一级av第二区| 搡老妇女老女人老熟妇| 一级毛片高清免费大全| 日本精品一区二区三区蜜桃| 成人亚洲精品av一区二区| 日韩免费av在线播放| 亚洲精华国产精华精| 国产亚洲精品综合一区在线观看 | 久久国产精品人妻蜜桃| 国产精品香港三级国产av潘金莲| 亚洲免费av在线视频| 精品国产超薄肉色丝袜足j| 亚洲av第一区精品v没综合| 啦啦啦韩国在线观看视频| 日韩 欧美 亚洲 中文字幕| 亚洲人成电影免费在线| 夜夜看夜夜爽夜夜摸| 亚洲欧美一区二区三区黑人| 午夜日韩欧美国产| 黄色视频不卡| 国产精品永久免费网站| av中文乱码字幕在线| 高清毛片免费观看视频网站| 国产精品一区二区精品视频观看| 在线观看美女被高潮喷水网站 | 欧美午夜高清在线| 亚洲国产中文字幕在线视频| 伊人久久大香线蕉亚洲五| a在线观看视频网站| 好看av亚洲va欧美ⅴa在| 99热只有精品国产| 久久伊人香网站| 日本五十路高清| 亚洲色图 男人天堂 中文字幕| 757午夜福利合集在线观看| 麻豆国产97在线/欧美 | 久久人人精品亚洲av| 听说在线观看完整版免费高清| 免费在线观看成人毛片| 日韩精品中文字幕看吧| 日本精品一区二区三区蜜桃| 99久久精品国产亚洲精品| 久久精品亚洲精品国产色婷小说| 国产蜜桃级精品一区二区三区| 又大又爽又粗| 久久 成人 亚洲| 好男人电影高清在线观看| videosex国产| 99久久综合精品五月天人人| 在线观看66精品国产| 久久精品国产清高在天天线| 操出白浆在线播放| 不卡av一区二区三区| 男女下面进入的视频免费午夜| 一级片免费观看大全| 亚洲av电影不卡..在线观看| 久久久久久免费高清国产稀缺| 亚洲熟妇中文字幕五十中出| av福利片在线| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕人成人乱码亚洲影| 毛片女人毛片| 久久久久免费精品人妻一区二区| 国产精品av视频在线免费观看| 久久天躁狠狠躁夜夜2o2o| 正在播放国产对白刺激| 悠悠久久av| 亚洲av电影在线进入| 狂野欧美激情性xxxx| 亚洲国产精品合色在线| 免费在线观看完整版高清| 国产av又大| 精品久久久久久久人妻蜜臀av| 久久精品国产清高在天天线| 国产片内射在线| 啦啦啦免费观看视频1| 免费在线观看成人毛片| 亚洲av成人不卡在线观看播放网| 欧美成狂野欧美在线观看| av福利片在线观看| 母亲3免费完整高清在线观看| 欧美日本视频| 九九热线精品视视频播放| 久久国产精品人妻蜜桃| 国产主播在线观看一区二区| 99热这里只有精品一区 | 黄色女人牲交| 午夜成年电影在线免费观看| 亚洲精品久久成人aⅴ小说| 欧美丝袜亚洲另类 | 女生性感内裤真人,穿戴方法视频| 国产又黄又爽又无遮挡在线| 91麻豆精品激情在线观看国产| 老司机福利观看| 亚洲美女黄片视频| aaaaa片日本免费| 12—13女人毛片做爰片一| 丝袜人妻中文字幕| 久久国产精品影院| 大型av网站在线播放| 成人永久免费在线观看视频| 国内毛片毛片毛片毛片毛片| 韩国av一区二区三区四区| 国产一区二区激情短视频| 久久午夜综合久久蜜桃| 香蕉av资源在线| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区在线观看日韩 | videosex国产| 色综合婷婷激情| 国产一区二区激情短视频| 亚洲欧美日韩高清在线视频| 国产亚洲精品第一综合不卡| 欧美人与性动交α欧美精品济南到| 亚洲中文字幕日韩| 亚洲精品久久国产高清桃花| 国产亚洲精品久久久久久毛片| 色播亚洲综合网| 757午夜福利合集在线观看| 国内毛片毛片毛片毛片毛片| 又黄又粗又硬又大视频| 一个人免费在线观看电影 | 国产精品久久久久久亚洲av鲁大| 中文字幕精品亚洲无线码一区| 久久久久久九九精品二区国产 | 亚洲成人精品中文字幕电影| www国产在线视频色| 亚洲美女黄片视频| 亚洲第一欧美日韩一区二区三区| 精华霜和精华液先用哪个| 首页视频小说图片口味搜索| 亚洲国产精品999在线| aaaaa片日本免费| 波多野结衣高清无吗| 女警被强在线播放| 欧美+亚洲+日韩+国产| 国产精品美女特级片免费视频播放器 | 亚洲成人久久性| 精品欧美国产一区二区三| 国产亚洲精品久久久久久毛片| 免费看十八禁软件| 久久精品91蜜桃| 午夜免费激情av| 女同久久另类99精品国产91| 午夜免费激情av| 熟女电影av网| 久久九九热精品免费| 在线a可以看的网站| 老鸭窝网址在线观看| 免费无遮挡裸体视频| 欧美在线一区亚洲| 又紧又爽又黄一区二区| 午夜福利视频1000在线观看| 老司机靠b影院| 久久久国产成人免费| 12—13女人毛片做爰片一| 久9热在线精品视频| 999精品在线视频| 久久久国产成人免费| 国产高清视频在线播放一区| 非洲黑人性xxxx精品又粗又长| 日韩精品青青久久久久久| 99热6这里只有精品| av天堂在线播放| 99热6这里只有精品| 午夜免费成人在线视频| 亚洲免费av在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一区高清亚洲精品| 啦啦啦免费观看视频1| 欧美日韩亚洲国产一区二区在线观看| 久久精品国产亚洲av香蕉五月| 成人18禁在线播放| 99热只有精品国产| 日日爽夜夜爽网站| av在线天堂中文字幕| 精品久久久久久久久久免费视频| av国产免费在线观看| 欧美人与性动交α欧美精品济南到| 曰老女人黄片| 香蕉丝袜av| 看免费av毛片| 黄片小视频在线播放| 老司机午夜福利在线观看视频| 黄色 视频免费看| 日本 欧美在线| 黑人操中国人逼视频| 亚洲精品色激情综合| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲精品第一综合不卡| 18禁黄网站禁片午夜丰满| 全区人妻精品视频| 国产精品久久视频播放| 亚洲一区中文字幕在线| 免费观看人在逋| 国产av麻豆久久久久久久| 他把我摸到了高潮在线观看| 熟女少妇亚洲综合色aaa.| 男人舔奶头视频| 成年人黄色毛片网站| 国产黄片美女视频| 亚洲av熟女| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩亚洲综合一区二区三区_| 一本大道久久a久久精品| 日韩av在线大香蕉| 亚洲精品在线美女| 精品一区二区三区四区五区乱码| 五月伊人婷婷丁香| 国产片内射在线| 亚洲熟妇中文字幕五十中出| 亚洲成av人片在线播放无| 日本a在线网址| 三级男女做爰猛烈吃奶摸视频| 久久九九热精品免费| 亚洲 国产 在线| 欧美3d第一页| 少妇被粗大的猛进出69影院| 99热这里只有是精品50| 精品久久久久久久毛片微露脸| 又黄又粗又硬又大视频| 久久精品成人免费网站| 久久久久久九九精品二区国产 | 他把我摸到了高潮在线观看| 国产精品电影一区二区三区| 极品教师在线免费播放| 婷婷精品国产亚洲av| 成年人黄色毛片网站| 此物有八面人人有两片| 欧美另类亚洲清纯唯美| 美女黄网站色视频| 午夜日韩欧美国产| 91成年电影在线观看| 精品国产乱子伦一区二区三区| 床上黄色一级片| 亚洲一区二区三区不卡视频| 国产精品久久视频播放| 午夜福利高清视频| 国产单亲对白刺激| 国产精品久久久人人做人人爽| 一级毛片女人18水好多| 欧美日韩亚洲综合一区二区三区_| 白带黄色成豆腐渣| 真人一进一出gif抽搐免费| 成人国产一区最新在线观看| 老司机福利观看| 亚洲,欧美精品.| 看免费av毛片| 国内毛片毛片毛片毛片毛片| 在线a可以看的网站| 99精品久久久久人妻精品| 色在线成人网| 国产人伦9x9x在线观看| 午夜福利免费观看在线| 国产黄片美女视频| 欧美日韩一级在线毛片| 欧美在线黄色| 精品一区二区三区av网在线观看| 叶爱在线成人免费视频播放| 小说图片视频综合网站| 久久久精品国产亚洲av高清涩受| 国内揄拍国产精品人妻在线| 99精品在免费线老司机午夜| 男女之事视频高清在线观看| 亚洲 欧美 日韩 在线 免费| e午夜精品久久久久久久| 日本 av在线| 脱女人内裤的视频| 亚洲七黄色美女视频| 亚洲人成网站在线播放欧美日韩| 久久久精品国产亚洲av高清涩受| 欧美黑人精品巨大| 日日爽夜夜爽网站| 亚洲aⅴ乱码一区二区在线播放 | 又爽又黄无遮挡网站| 久久精品国产清高在天天线| 国产亚洲精品综合一区在线观看 | 免费在线观看影片大全网站| 狂野欧美激情性xxxx| 久久中文看片网| 色综合婷婷激情| 亚洲国产欧美一区二区综合| 一进一出抽搐gif免费好疼| 久久久久国产一级毛片高清牌| 欧美黄色片欧美黄色片| 国产伦在线观看视频一区| 久久精品91无色码中文字幕| 日韩精品青青久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 国产精品日韩av在线免费观看| 国产精品一区二区三区四区免费观看 | 国产精品久久久久久精品电影| 久久 成人 亚洲| 人成视频在线观看免费观看| 国产爱豆传媒在线观看 | 国产欧美日韩精品亚洲av| 亚洲自偷自拍图片 自拍| 男插女下体视频免费在线播放| 日本五十路高清| 国产人伦9x9x在线观看| 欧美性猛交╳xxx乱大交人| 一级毛片精品| 国产成人精品无人区| 黄色 视频免费看| 欧美一级毛片孕妇| 99在线人妻在线中文字幕| 在线十欧美十亚洲十日本专区| 91国产中文字幕| 亚洲自拍偷在线| 日本一本二区三区精品| 丝袜美腿诱惑在线| 亚洲熟女毛片儿| 欧美丝袜亚洲另类 | 亚洲成av人片免费观看| 欧美av亚洲av综合av国产av| 夜夜躁狠狠躁天天躁| 免费看美女性在线毛片视频| 亚洲成人免费电影在线观看| 国产精品亚洲av一区麻豆| 宅男免费午夜| 特级一级黄色大片| 1024手机看黄色片| 男女下面进入的视频免费午夜| 黄色毛片三级朝国网站| 色精品久久人妻99蜜桃| 国产av在哪里看| 九色国产91popny在线| 欧美日韩乱码在线| 亚洲一区二区三区色噜噜| 无人区码免费观看不卡| 欧美午夜高清在线| 国产亚洲精品一区二区www| 91国产中文字幕| 丁香欧美五月| 看免费av毛片| 久久亚洲真实| 国产高清视频在线播放一区| 精品欧美一区二区三区在线| 宅男免费午夜| 国产亚洲精品综合一区在线观看 | 老熟妇仑乱视频hdxx| 国产激情久久老熟女| 一二三四社区在线视频社区8| 欧美激情久久久久久爽电影| 狂野欧美白嫩少妇大欣赏| 国产蜜桃级精品一区二区三区| 黄色毛片三级朝国网站| 欧美黑人巨大hd| 午夜视频精品福利| 国产激情久久老熟女| 露出奶头的视频| 一区二区三区激情视频| 色综合欧美亚洲国产小说| 久久精品综合一区二区三区| 在线观看免费视频日本深夜| 国产视频内射| 长腿黑丝高跟| 亚洲欧美精品综合一区二区三区| 日本撒尿小便嘘嘘汇集6| 国产成年人精品一区二区| 欧美激情久久久久久爽电影| 午夜激情福利司机影院| 香蕉av资源在线| 国产免费男女视频| 99精品久久久久人妻精品| 国产精品久久久人人做人人爽| 国产伦一二天堂av在线观看| 久久欧美精品欧美久久欧美| 国产精品 欧美亚洲| av在线天堂中文字幕| 我的老师免费观看完整版| 蜜桃久久精品国产亚洲av| 天天添夜夜摸| 亚洲一区二区三区色噜噜| 午夜成年电影在线免费观看| 欧美性猛交╳xxx乱大交人| 我的老师免费观看完整版| 成人18禁在线播放| 少妇裸体淫交视频免费看高清 | 又大又爽又粗| 精华霜和精华液先用哪个| 亚洲中文日韩欧美视频| 青草久久国产| 亚洲片人在线观看| 亚洲一区二区三区不卡视频| 一级作爱视频免费观看| 久久天堂一区二区三区四区| 亚洲va日本ⅴa欧美va伊人久久| 亚洲人与动物交配视频| 岛国视频午夜一区免费看| 老司机在亚洲福利影院| 国产97色在线日韩免费|