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

        ?

        封頭補強仿形織物的曲面幾何映射研究

        2024-01-03 02:37:48馬文鎖許泳豪馬振宇漫睿東
        產(chǎn)業(yè)用紡織品 2023年9期

        姬 智 馬文鎖 許泳豪 馬振宇 漫睿東

        河南科技大學(xué) 機電工程學(xué)院,河南 洛陽 471003

        以碳纖維復(fù)合材料代替金屬材料作為殼體的壓力容器具有質(zhì)輕、比強度和比模量高等優(yōu)點,在航空航天及儲氣領(lǐng)域得到廣泛的應(yīng)用[1-5]。但在水壓爆破試驗過程中發(fā)現(xiàn),碳纖維殼體壓力容器的封頭部分,尤其是接近金屬件邊緣處,容易產(chǎn)生應(yīng)力集中區(qū)域,從而導(dǎo)致碳纖維殼體發(fā)生低壓爆破。目前,常用的解決方法是對殼體封頭部分進行補強[6]。二維織物生產(chǎn)工藝成熟、成本低、可織造的尺寸范圍廣,且經(jīng)精確設(shè)計的織物結(jié)構(gòu)及其變形范圍能夠基本滿足復(fù)合材料構(gòu)件的織物設(shè)計要求,故被用作復(fù)合材料壓力容器封頭區(qū)域的補強織物[7-10]。

        傳統(tǒng)經(jīng)緯正交織造方法分段加工的平面二維織物在覆蓋封頭曲面時,紗線無法精確滿足封頭曲面的母線與環(huán)向正交的工程需求,且由于覆蓋過程中二維織物受到載荷的作用,紗線及纖維產(chǎn)生局部滑移,聚集狀態(tài)發(fā)生改變,纖維體積百分含量、力學(xué)性能及其復(fù)合材料構(gòu)件性能等發(fā)生非均勻變化[11-13]。仿形織物是按部件的外形輪廓設(shè)計織制的單層機織物,織物可呈曲面結(jié)構(gòu),在成型時可將織物纏繞或?qū)雍掀饋?使紗線連續(xù)的同時還可以確保經(jīng)緯紗之間的夾角[14]。仿形織造技術(shù)為高性能復(fù)合材料的低成本制造開辟了新的途徑。仿形織物與目標(biāo)曲面貼合的關(guān)鍵在于目標(biāo)曲面的精確展開。鄭鵬飛等[15]根據(jù)測地線的短程性與線上各點測地曲率為零的幾何特性,采用在不可展的封頭曲面上構(gòu)建由測地線和法截線交織成線網(wǎng)的方法,獲得橢圓形和球形封頭展開圖的邊界。梁堰波等[16]基于能量模型將待展開曲面以三角網(wǎng)格形式表示,將曲面簡化為三角網(wǎng)格頂點組成的質(zhì)點系統(tǒng),使曲面網(wǎng)格的劃分更具有操作性。常用的曲面展開方法和技術(shù)還有幾何逼近法[17]、力學(xué)模型法[18]、能量方程法[19-20]和三維軟件法[21]等。

        合理的二維織物與曲面幾何覆蓋映射之間的關(guān)系,有利于解決無缺陷覆蓋的織物結(jié)構(gòu)設(shè)計難題。本文基于封頭空間曲面的剖分式理想展開方法,研究二維織物與封頭曲面之間的幾何映射關(guān)系,提出可覆蓋封頭曲面的二維織物的加工工藝參數(shù)理論計算方法,為封頭補強仿形織物的成型加工提供理論參考。

        1 封頭補強織物成型方法及二維織物的單元類型

        1.1 封頭補強織物成型方法

        與封頭曲面外形特征相吻合的仿形織物的成型方法有兩種,如圖1所示:一種是對平面的二維織物進行裁剪,再通過織物變形進行覆蓋,成型為所需要的封頭補強織物,即平面覆蓋法;另一種是在滿足封頭外形要求的前提下,設(shè)計加工工藝參數(shù),在線一次仿形加工出完整的曲面補強織物,使用時將織物纏繞或?qū)雍掀饋?即仿形覆蓋法。平面覆蓋法補強織物產(chǎn)品雖然可以較好地貼合封頭曲面,但部分紗線的裁斷將導(dǎo)致補強織物的力學(xué)性能降低,影響補強效率和質(zhì)量可靠性;仿形覆蓋法加工的織物,其環(huán)向紗線連續(xù)且與經(jīng)紗方向一致,可以很好地貼合封頭曲面表面,同時可避免織物由于紗線裁斷而導(dǎo)致的力學(xué)性能損失[22]。

        圖1 曲面封頭補強織物成型方法

        1.2 二維織物的單元類型

        用于覆蓋封頭曲面的二維織物的結(jié)構(gòu)組織主要有平紋組織、斜紋組織和緞紋組織[23]。這3種組織的二維織物的一般單元類型為正交單元,其中包括正方單元和矩形單元。正交單元在外力作用下可變形為可逆平行四邊形單元和極限剪切單元。可逆平行四邊形單元在二維平面內(nèi)結(jié)構(gòu)不穩(wěn)定,在覆蓋為封頭曲面時,需實現(xiàn)經(jīng)緯紗正交,從而使單元結(jié)構(gòu)達到穩(wěn)定狀態(tài)。極限剪切單元為二維織物單元的最大剪切狀態(tài),具有極限穩(wěn)定性,在外力作用下很難發(fā)生變形。3種組織的二維織物的各單元類型如表1所示。

        表1 二維織物的一般單元類型

        利用二維織物覆蓋封頭曲面時,考慮紗線截面及幾何形狀在一定范圍內(nèi)的變形,將織物單元的變形分為拉伸、壓縮和剪切3種,即可求得可覆蓋封頭表面的織物最大剪切單元。該單元為表面覆蓋的最小可重復(fù)單元。使用尺寸不小于最小可重復(fù)單元加工的二維織物,若其變形后仍能夠滿足工程需求,并且可無缺陷覆蓋封頭曲面,則稱這種單元為該二維織物的適配單元(不唯一)。適配單元的幾何尺寸決定封頭曲面的覆蓋程度。二維織物及其適配單元示意圖如圖2所示。

        圖2 二維織物適配單元

        2 仿形覆蓋加工封頭曲面的二維展開方案

        封頭織物結(jié)構(gòu)如圖3所示,封頭曲面由球頭面與圓柱面組合而成。經(jīng)紗與徑圓重合、緯紗與回轉(zhuǎn)曲面母線重合,被認(rèn)為是封頭曲面最理想的覆蓋方式。為降低織造難度,在封頭補強織物實際加工過程中可采用分段織造、間斷卷曲的加工工藝,使補強織物的局部仍表現(xiàn)為平面織物。針對封頭補強織物與封頭曲面貼合的部分,將回轉(zhuǎn)封頭曲面等分,采用每份曲面上覆蓋的織物緯紗與母線重合、經(jīng)紗與徑圓重合的織物近似展平方法。

        圖3 封頭織物結(jié)構(gòu)

        為考察經(jīng)、緯紗展平后的軌跡,將封頭曲面沿母線和環(huán)向進行分割,再將分割后的經(jīng)紗層逐條展開,其展開面的形狀為扇形,封頭圓柱部分的展開面形狀為矩形,如圖4b)所示。封頭曲面每一層展開扇形的圓弧半徑不同,但圓心均分布在中心線上。剖分式理想展開方案使得每一經(jīng)紗層均能精確展開,但每一層與相鄰層之間存在間隙,且間隙的大小均不相同,層與層之間的最大間隙為Emax。并且由圖4c)可知,在剖分式理想展開過程中,平行于經(jīng)紗方向的紗線連續(xù),但緯紗方向的網(wǎng)格不連續(xù)。二維織物的經(jīng)紗在織造過程中與緯紗正交,可沿經(jīng)紗方向?qū)暭唵卧W(wǎng)格進行近似展開。單元網(wǎng)格的長度(又稱環(huán)向單元長)記為w,單元網(wǎng)格的寬度記為h。

        圖4 剖分式理想展開

        二維織物的織造過程中,經(jīng)紗是相互平行的,理想展開經(jīng)紗的角度較難控制,因為受織物加工工藝的影響。將理想展開的經(jīng)紗投影在相互平行的位置,如圖5所示,即為平行經(jīng)紗(徑圓)的近似展開方案。取短徑長為Lmin,經(jīng)紗被分為5份,短徑和長徑在平行經(jīng)紗上的投影長度不變,交織點的長度不發(fā)生改變,相鄰層經(jīng)紗展開時的短圓弧長和長圓弧長相等,則在相鄰平行經(jīng)紗的展開過程中,緯紗為傾斜的紗線。具體過程如圖6所示,緯紗與經(jīng)紗的交織點均分布在平行經(jīng)紗上,將經(jīng)紗逐層理想展開后,網(wǎng)格劃分后的封頭曲面如圖6a)所示,將理想展開的經(jīng)紗投影在相互平行的位置后,連接同一緯紗的交織點,形成具有一定傾斜角度的緯紗,獲得近似展開的二維織物,即經(jīng)紗相互平行,緯紗呈弧線分布,如圖6b)所示。隨著交織點數(shù)量的增加,即經(jīng)密和緯密的增大,緯紗的弧線近似平滑,當(dāng)平行經(jīng)紗的數(shù)量較多時,單元尺寸的長度較小,緯紗可以通過一定弧度進行打緯,形成近似展開的二維織物,覆蓋封頭曲面。

        圖5 平行經(jīng)紗近似展開

        圖6 平行經(jīng)紗近似展開二維織物

        由于單元網(wǎng)格的尺寸不同,剖分式理想展開后織物與曲面幾何覆蓋映射之間的誤差也不一樣。單元尺寸取決于紗線參數(shù)和工藝參數(shù),進而影響層與層之間的最大間隙Emax。封頭曲面的等分?jǐn)?shù)同樣影響二維織物的覆蓋精度,設(shè)xoy坐標(biāo)面上封頭曲面的等分?jǐn)?shù)為m,則封頭曲面的幾何關(guān)系如圖7所示。

        圖7 封頭曲面的幾何關(guān)系

        封頭半徑為R的圓柱曲面可精確展開,封頭的局部球面本身為不可展曲面,其小端半徑為r,大端半徑為R。被等分的封頭曲面每一份的大端弧長為Li,小端弧長為li??傻闷鋷缀侮P(guān)系為

        (1)

        設(shè)回轉(zhuǎn)曲面r(θ,φ)由xoz坐標(biāo)面上的曲線L1和L2繞z軸旋轉(zhuǎn)而成,則曲面以θ、φ為參數(shù)的曲面方程為

        r(θ,φ)={f(φ)cosθ,f(φ)sinθ,g(φ)}

        (2)

        曲面r(θ,φ)中的一點映射為切柱面中同緯度直線上的一點,取曲線的θ向參數(shù)為Θ,對應(yīng)的可展切柱面的參數(shù)方程為

        (3)

        封頭不可展曲面的參數(shù)方程為

        x=f1(θ,φ),y=f2(θ,φ),z=f3(θ,φ)

        (4)

        回轉(zhuǎn)母線L1:

        x=f(φ)=Rcosφ,z=g(φ)=Rsinφ

        (5)

        由式(2)、式(3)和式(5)得,球面片:

        r′(θ,φ)={Rcosφcosθ,Rcosφsinθ,Rsinφ}

        (6)

        切柱面片:

        (7)

        結(jié)合式(1)~式(7)即可獲得封頭曲面等分后的曲面幾何關(guān)系。

        3 二維織物和封頭曲面的幾何映射

        封頭曲面覆蓋的織物與相應(yīng)的二維織物存在幾何對應(yīng)關(guān)系。假設(shè)紗線橫截面為橢圓形且面積恒定為S=πab(a和b分別為橢圓的長半軸和短半軸);單元中紗線的長度不變,且織物在變形前后單元的數(shù)量不變。對單個不可展曲面,通過剖分式理想展開和平行經(jīng)紗近似展開,將空間不可展曲面上的單元映射在二維平面,如圖8a)和圖8b)所示;再通過研究網(wǎng)格單元的變形,得到二維織物的單元合適形狀和尺寸范圍,如圖8c)所示。封頭曲面是由球頭面與圓柱面形成的組合曲面。Ⅰ區(qū)域?qū)?yīng)圓柱曲面的展開,單元僅發(fā)生拉伸變形或壓縮變形,單元的4個節(jié)點組成的形狀仍為矩形。Ⅱ區(qū)域?qū)?yīng)球頭面的展開,單元發(fā)生剪切變形,可通過單元節(jié)點的變化,探究曲母線回轉(zhuǎn)面與二維織物的幾何映射關(guān)系。

        圖8 適配單元的近似展開

        圖9 曲母線回轉(zhuǎn)面

        在xoz平面內(nèi),曲母線z=-px2+H上切線L1的斜率為K1(K1=-2px),將過坐標(biāo)原點o的曲母線切線line 1的法線line 2斜率用K2表示,由K1K2=-1可求得line 1與line 2的交點A(x1,z1)。點A與x軸上的空間曲面底面半徑BC構(gòu)成一個平面(z1/x1)x+z=0,該平面與空間曲面相交構(gòu)成空間曲面的剖切面ABC,相交曲線f(x,y,z)為剖切面ABC的輪廓線,如圖10所示。曲面可對稱2次完成覆蓋及幾何運算,覆蓋半回轉(zhuǎn)曲面時,經(jīng)緯向紗線分別與曲線f(x,y,z)和f(y,z)吻合并均勻映射。原始平紋織物單元邊長均勻變形均布曲線段后,再向Ⅰ、Ⅱ、Ⅲ和Ⅳ區(qū)域覆蓋映射,如圖11所示。

        圖10 回轉(zhuǎn)曲面幾何關(guān)系

        圖11 回轉(zhuǎn)曲面分區(qū)及幾何映射方法

        紗線在不同區(qū)域內(nèi)任意位置發(fā)生不同程度的對角協(xié)調(diào)變形,且圖11中的Ⅰ、Ⅲ區(qū)域和Ⅱ、Ⅳ區(qū)域關(guān)于平面yoz對稱。二維織物在區(qū)域內(nèi)部發(fā)生對角協(xié)調(diào)變形,一個二維織物單元用4個節(jié)點描述,連續(xù)變化的四節(jié)點單元構(gòu)成空間區(qū)域。以點A為基點,對應(yīng)區(qū)域Ⅰ紗線節(jié)點為W(i,j),i表示沿曲線f4(x,y,z)方向的節(jié)點序列,j表示沿曲線f1(x,z)方向的節(jié)點序列。在W(1,j)和W(i,1)上單元發(fā)生內(nèi)稟變形。從基點A開始在曲線f(x,y,z)上和曲線f(y,z)上以二維織物內(nèi)稟變形后的單元邊長L0等分曲線段,尋找到一系列幾何節(jié)點W(1,k+1)和W(n+1,1)并作為覆蓋初始點,其中k和n分別表示需要覆蓋曲線f(x,y,z)和f(y,z)的單元數(shù)量,xi,j,yi,j,zi,j為節(jié)點W(i,j)在曲面上的空間坐標(biāo)。

        用基點A、點W(1,2)和點W(2,1)這3個節(jié)點尋找第一個單元的第四個節(jié)點,并確定可能的單元初始紗線夾角α1,1,換算對應(yīng)夾角的單元長度L1,1=g(α1,1)。依次疊加獲得區(qū)域I內(nèi)所有單元經(jīng)緯紗線的夾角αi,j,并換算成對應(yīng)的變形后單元邊長Li,j=g(αi,j)。αi,j和Li,j之間滿足如下關(guān)系:

        (8)

        以點A為基點在區(qū)域Ⅲ內(nèi)以相同的方法疊加運算,即可獲得完整幾何覆蓋算法??筛采w二維織物展開幾何輪廓如圖12所示,用設(shè)計的紗線股數(shù)和經(jīng)緯密加工二維織物,保證經(jīng)緯紗數(shù)量Q1和Q2,并按文中所述覆蓋方法即可精確覆蓋所述空間曲面。

        圖12 半球覆蓋映射

        4 二維可織織物經(jīng)緯密的理論計算方法

        二維織物的結(jié)構(gòu)參數(shù)及外輪廓形狀取決于其經(jīng)緯密的大小。為區(qū)別于傳統(tǒng)的二維織物,后文將與封頭曲面相對應(yīng)的仿形織物稱作二維可織織物。二維可織織物織造時經(jīng)紗的排布如圖13所示,經(jīng)紗沿著封頭曲面母線均勻排布。隨著封頭曲面大端到小端圓弧半徑的逐漸減小,緯紗數(shù)量逐漸減少。沿封頭曲面的回轉(zhuǎn)母線,對封頭曲面按照第2節(jié)描述的剖分式理想展開法進行精確展開。取m等分后的一個曲面,曲面圓弧對應(yīng)封頭曲面內(nèi)的角度為2π/m,不可展曲面部分沿平行經(jīng)紗層逐條展開,展開面形狀為扇形;可展曲面部分的展開形狀為矩形,如圖14a)所示。

        圖13 織造過程中的經(jīng)紗排列

        圖14 剖分式理想展開的經(jīng)紗層間關(guān)系

        精確展開的經(jīng)紗層與層之間的關(guān)系如圖14b)所示,第j層與相鄰層之間的公用圓弧長度相等,但由于每一層的圓弧展開角度不同,導(dǎo)致各層之間存在間隙。經(jīng)紗展開后的空間弧線段長度為

        (9)

        假設(shè)最長空間弧線段具有K個節(jié)點,其將弧線等分為K-1段,各段的長度與紗線的參數(shù)相關(guān),且均勻分布在經(jīng)紗上。隨著封頭曲面大端到小端圓弧半徑的逐漸減小,空間弧線段上的節(jié)點數(shù)減少。由于單元的長度w是固定的,且空間弧線段上的節(jié)點數(shù)是確定的,最終得到展開圓弧段的長度為

        (10)

        在經(jīng)紗層的展開過程中,短圓弧網(wǎng)格展開后的理論長度和實際長度存在差異,將該誤差記為E1;同理,長圓弧網(wǎng)格展開后的理論長度與實際長度存在差異,將該誤差記為E2。經(jīng)紗方向的單元節(jié)點數(shù)為Kj。由經(jīng)紗方向的單元節(jié)點數(shù)即可知緯紗數(shù)量。設(shè)P為間隙占單元長度的百分比,其值越小,表明展開的精度越高。則各個參數(shù)之間存在式(11)的關(guān)系。

        (11)

        式中:Ej為第j層經(jīng)紗展開圓弧段的長度誤差。

        在理想展開法的精確度分析過程中,理想單元與實際單元的重合度決定了封頭曲面的展開精度。假設(shè)網(wǎng)格劃分過程中,曲面等分?jǐn)?shù)m=8,織物結(jié)構(gòu)單元為正方單元,即w=h;當(dāng)單元尺寸分別取10、8、6和4時,織物單元理想分布(矩形方框表示)和實際分布(菱形方框)如圖15所示。

        圖15 單元尺寸對單元分布的影響

        由圖15可知,在封頭曲面結(jié)構(gòu)參數(shù)相同的情況下,單元理想分布和實際單元分布的重合度隨著單元尺寸的減小而增大,即封頭曲面的展開精度隨著單元尺寸的減小而增大;且單元的尺寸大小影響短圓弧和長圓弧網(wǎng)格展開后的理論長度與實際長度的誤差,短圓弧誤差E1和單元節(jié)點數(shù)量K隨著單元尺寸的減小而增大,長圓弧誤差E2和層與層之間的最大間隙Emax隨著單元尺寸的減小而減小。則當(dāng)單元尺寸設(shè)計得足夠小時,E1和E2的值較小,即長短圓弧展開的理論長度和實際長度相差較小,同時Emax值也較小,即層與層之間的最大間隙較小。當(dāng)長短圓弧展開的理論長度和實際長度的誤差和層與層之間的最大間隙小到可以忽略時,即可實現(xiàn)圖6所示的平行經(jīng)紗近似展開。在網(wǎng)格劃分過程中,等分?jǐn)?shù)的大小同樣影響單元理想分布與實際分布的誤差值,本文不做過多討論。

        5 試驗驗證

        5.1 試驗材料與裝置

        被覆蓋的封頭曲面如圖16所示,其大端半徑R=43 mm,小端半徑r=16 mm。在網(wǎng)格劃分過程中,將回轉(zhuǎn)封頭曲面等分,等分?jǐn)?shù)m=4,w=h=4 mm。手工織制出與封頭曲面相對應(yīng)的二維可織織物,所用國產(chǎn)碳纖維紗線的規(guī)格如表2所示??椫频亩S可織織物(圖17)與等分后的部分封頭曲面相對應(yīng)。

        表2 碳纖維紗線的規(guī)格

        圖16 被覆蓋的封頭曲面

        圖17 織物試樣

        5.2 結(jié)果與討論

        織物鋪覆效果如圖18所示。由圖18可以看出,二維可織織物覆蓋后,封頭曲面部分上的經(jīng)紗與封頭曲面的徑圓基本重合,緯紗傾斜一定角度,經(jīng)紗和緯紗基本處于極限剪切穩(wěn)定狀態(tài),且緯紗密度從小端半徑至大端半徑逐漸增大,經(jīng)緯紗的位置形狀與本文理論分析部分相吻合,驗證了封頭曲面剖分式理想展開法的可行性,以及封頭曲面與二維可織織物之間幾何映射關(guān)系的正確性。

        圖18 鋪覆后的封頭曲面

        然而,由于鋪覆所用的二維可織織物采用手工織制,鋪覆后的封頭曲面部分紗線間存在較大間隙,尤其是曲面母線附近。出現(xiàn)這種現(xiàn)象的原因可能為:一是手工織制的二維可織織物本身就存在較大間隙,織物在外力作用下鋪覆在封頭曲面表面上,為適應(yīng)曲率不斷變化的封頭曲面,紗線會產(chǎn)生旋轉(zhuǎn)和滑移,從而造成間隙的擴大,采用專業(yè)的織機織造可以盡可能避免此類問題;二是在劃分網(wǎng)格過程中,曲面等分?jǐn)?shù)和單元尺寸選取不合理,導(dǎo)致實際單元與理想單元之間的分布誤差較大。兩種因素均會導(dǎo)致織造出的織物與封頭曲面不能完全貼合。碳纖維織物在與基體材料進行復(fù)合時,紗線間隙的存在會形成復(fù)合材料富樹脂區(qū),進而導(dǎo)致封頭壓力容器的力學(xué)性能下降,因此在實際生產(chǎn)過程中應(yīng)盡可能避免紗線間隙的形成[24-25]。

        6 結(jié)論

        本文以三元組織二維織物覆蓋封頭曲面的展開方法為背景,提出了可確保覆蓋二維織物經(jīng)緯紗線連續(xù)和正交的封頭曲面剖分式理想展開法;探討了封頭曲面覆蓋織物與二維可織織物存在的幾何映射關(guān)系,并對二維可織織物的織造工藝參數(shù)進行設(shè)計。研究發(fā)現(xiàn)理想展開后的封頭曲面在結(jié)構(gòu)參數(shù)相同的情況下,曲面展開精度隨著單元尺寸的減小而增大,且各層經(jīng)紗的短圓弧誤差和單元節(jié)點數(shù)量隨著單元尺寸的減小而增大,長圓弧誤差和層與層之間的最大間隙隨著單元尺寸的減小而減小。即當(dāng)單元尺寸設(shè)計得足夠小時,緯紗的弧線近似平滑,緯紗可以通過一定弧度進行打緯,利用平行經(jīng)紗近似展開形成相應(yīng)的封頭補強二維可織織物。最后,本文進行了二維可織織物的織制與封頭覆蓋試驗,結(jié)果表明,使用相應(yīng)的二維可織織物覆蓋封頭后,經(jīng)緯紗的位置形狀與理論分析部分基本相吻合,驗證了封頭曲面剖分式理想展開法的可行性及封頭曲面與二維可織織物之間幾何映射關(guān)系的正確性。

        亚洲av乱码二区三区涩涩屋| 无码国产精品第100页| 丝袜人妻无码中文字幕综合网| 国产一区二区熟女精品免费| 婷婷四虎东京热无码群交双飞视频| 国产成人精品电影在线观看| 国产成人精品精品欧美| 中国老太老肥熟女视频| 97超碰精品成人国产| 亚洲欧美日韩成人高清在线一区| 婷婷综合久久中文字幕蜜桃三电影| 亚洲一区二区精品久久岳| 蜜桃精品视频一二三区| 伊人久久大香线蕉av色| 亚洲av无码成人黄网站在线观看| 91网红福利精品区一区二| 人妻少妇偷人精品一区二区| 无码人妻一区二区三区免费视频 | 久久精品一区二区熟女| 国产婷婷色一区二区三区在线| 久久久噜噜噜久久中文字幕色伊伊| 亚洲午夜久久久精品国产| 激情五月我也去也色婷婷| 国产xxxxx在线观看| 中文字幕第八页| 女女同性av一区二区三区| 日本韩国男男作爱gaywww| 丰满少妇被猛男猛烈进入久久| 福利网在线| 高清不卡av一区二区| 亚洲av美国av产亚洲av图片| 久久AV中文一区二区三区| 国产亚洲精品一区二区在线播放 | 精品女厕偷拍视频一区二区区| 99久久婷婷国产综合亚洲| 日本高清www午色夜高清视频| 亚洲欧美日韩精品久久| 色妺妺视频网| 亚洲精品二区三区在线观看| 日本一区二区三区视频在线观看| 人妻av中文字幕无码专区|