999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于有限積分法的機(jī)織物電磁屏蔽效能仿真分析

2016-07-12 13:25:32蘇欽城趙曉明李衛(wèi)斌李建雄
紡織學(xué)報(bào) 2016年2期
關(guān)鍵詞:結(jié)構(gòu)模型

蘇欽城,趙曉明,李衛(wèi)斌,李建雄

(1.天津工業(yè)大學(xué) 紡織學(xué)院,天津 300387;2.天津工業(yè)大學(xué) 電子與信息工程學(xué)院,天津 300387)

基于有限積分法的機(jī)織物電磁屏蔽效能仿真分析

蘇欽城1,趙曉明1,李衛(wèi)斌1,李建雄2

(1.天津工業(yè)大學(xué) 紡織學(xué)院,天津 300387;2.天津工業(yè)大學(xué) 電子與信息工程學(xué)院,天津 300387)

為研究織物對(duì)電磁波的屏蔽機(jī)制,為電磁屏蔽織物開發(fā)提供理論參考,對(duì)織物微觀結(jié)構(gòu)進(jìn)行了建模,然后利用CST 微波工作室建立了電磁屏蔽織物三維電磁仿真模型。采用有限積分法計(jì)算了在1~18 GHz波段的織物電導(dǎo)率、電磁波入射角、織物密度和織物層數(shù)對(duì)屏蔽效能的影響。結(jié)果表明:織物電導(dǎo)率是影響其屏蔽效能的關(guān)鍵因素;屏蔽效能隨織物密度和層數(shù)的改變呈非線性變化;電磁波入射角增加會(huì)引起織物表面電流密度極化效應(yīng)。研究結(jié)果對(duì)電磁屏蔽織物具有較好適用性,可為設(shè)計(jì)與開發(fā)高性能電磁屏蔽織物提供理論參考。

電磁屏蔽織物;微觀結(jié)構(gòu)模型;屏蔽效能;電磁仿真

電磁干擾是當(dāng)今社會(huì)的公共污染之一,電磁屏蔽是防止電磁干擾的重要手段,其主要機(jī)制是通過對(duì)電磁波反射而起到防護(hù)效果。在常用的電磁屏蔽材料中,電磁屏蔽織物因其具有可塑性強(qiáng)、價(jià)格低和良好服用性能等特點(diǎn),已獲得廣泛的應(yīng)用[1-2]。科研人員已對(duì)影響織物電磁屏蔽效能的因素進(jìn)行了探討,包括不銹鋼纖維織物[2-3]、銅絲網(wǎng)格織物[1,4-6]等,分析了織物密度、緊度、導(dǎo)電纖維含量等因素對(duì)屏蔽效能的影響;對(duì)已有的金屬化織物研究包括對(duì)化學(xué)鍍銀[7-8]、磁控濺射[9]等金屬化織物的表面方塊電阻、金屬纖維含量等因素進(jìn)行實(shí)驗(yàn)分析。由于織物結(jié)構(gòu)復(fù)雜,影響因素較多,造成實(shí)驗(yàn)成本高、研發(fā)周期長(zhǎng)等弊端,因此基于電磁波傳播理論來定量研究織物材料與結(jié)構(gòu)參數(shù)對(duì)電磁屏蔽效能的影響,實(shí)現(xiàn)電磁屏蔽織物的數(shù)字化設(shè)計(jì)已迫在眉睫。Chen等[10]運(yùn)用網(wǎng)孔分析法分析了織物密度、電導(dǎo)率及織物層數(shù)等因素對(duì)屏蔽效能的影響;Rybicki等[11-12]用等效電路法將金屬間隔織物交織結(jié)構(gòu)進(jìn)行等效電路分析;李奇軍等[13]考慮不銹鋼纖維織物三維結(jié)構(gòu),假設(shè)該織物為金屬網(wǎng)進(jìn)行仿真分析。以上研究由于對(duì)織物微觀結(jié)構(gòu)考慮不足,對(duì)于屏蔽織物的紗線截面形狀、織物交織結(jié)構(gòu)等因素未曾涉及,因此普遍適用性較差。

本文以平紋織物為原型,分析織物微觀結(jié)構(gòu),進(jìn)行三維結(jié)構(gòu)建模,并建立電磁屏蔽織物仿真模型,采用有限積分法研究電磁屏蔽織物電導(dǎo)率、結(jié)構(gòu)特征及電磁波入射方向等3個(gè)方面與屏蔽效能之間的關(guān)系,以期對(duì)電磁屏蔽織物實(shí)現(xiàn)數(shù)字化設(shè)計(jì)。

1 織物結(jié)構(gòu)建模

織物結(jié)構(gòu)模型的建立是對(duì)紡織品進(jìn)行仿真研究的重要前提,其結(jié)構(gòu)參數(shù)改變會(huì)導(dǎo)致電磁屏蔽效能非線性變化[14]。TexGen軟件是由英國(guó)諾丁漢大學(xué)(University of Nottingham)設(shè)計(jì)開發(fā)的織物幾何結(jié)構(gòu)模擬軟件,可準(zhǔn)確模擬多種紗線結(jié)構(gòu)。本文使用TexGen對(duì)電磁屏蔽織物進(jìn)行微觀結(jié)構(gòu)建模,考慮紗線截面形狀、織物交織規(guī)律等結(jié)構(gòu)參數(shù),提高了織物結(jié)構(gòu)模型的精確度。

紗線是織物交織的結(jié)構(gòu)單元,其幾何結(jié)構(gòu)主要由紗線的軌跡及其截面形狀構(gòu)成。紗線軌跡可考慮為三維空間中紗線截面中心的連線,表達(dá)紗線長(zhǎng)度及其彎曲狀態(tài)。在TexGen建模中,紗線軌跡是由確定的點(diǎn)通過樣條函數(shù)等方法擬合得到。紗線截面形狀可微觀地反映出紗線交織狀態(tài),是紗線結(jié)構(gòu)的重要參數(shù)。由于紗線在織造過程中受擠壓,造成紗線截面不是圓形,Wong等[15]提出紗線截面為橢圓形和凸透鏡形。TexGen中的紗線橫截面有多種描述方式,常用的形狀如圖1所示。

圖1 紗線截面形狀Fig.1 Cross-sectional shape of yarn.(a) Ellipse; (b) Convex lens

織物是紗線的有機(jī)組合體,織造方式、織物組織等因素決定了織物最終結(jié)構(gòu)。機(jī)織物是由一系列的經(jīng)、緯紗交織組成,織物二維結(jié)構(gòu)模型將經(jīng)、緯紗簡(jiǎn)化為平面網(wǎng)格結(jié)構(gòu)[12],而無法考慮織物交織對(duì)其性能的影響。TexGen產(chǎn)生的紗線結(jié)構(gòu)通過適當(dāng)路徑組合,可構(gòu)成織物三維結(jié)構(gòu)模型,本文設(shè)計(jì)的平紋織物結(jié)構(gòu)模型如圖2所示。

圖2 平紋織物結(jié)構(gòu)模型Fig.2 Structure model of plain fabric.(a) Front view; (b) Lateral view; (c) Fabric laminated model

2 電磁仿真分析

本文采用CST軟件定義了織物組織單元周期邊界條件,然后使用頻域求解器的有限積分法對(duì)織物組織單元進(jìn)行仿真分析。對(duì)織物模型的假設(shè)條件為:1)紗線電導(dǎo)率各向同性;2)織物交織點(diǎn)處經(jīng)、緯紗接觸電阻為0;3)織物由無限個(gè)組織單元在同一平面上循環(huán)而成。將織物的三維模型導(dǎo)入CST微波工作室,設(shè)置織物在自由空間,電磁波沿z軸方向傳播,如圖3所示。以平紋織物為研究對(duì)象,設(shè)置電磁波頻率為1~18 GHz,探討織物電導(dǎo)率、電磁波入射角、織物密度和織物層數(shù)對(duì)電磁屏蔽效能的影響。

圖3 織物仿真模型Fig.3 Simulation model of fabric

2.1 織物電導(dǎo)率對(duì)屏蔽效能的影響

為系統(tǒng)研究織物電導(dǎo)率對(duì)屏蔽效能的影響,設(shè)計(jì)仿真模型的織物組織單元長(zhǎng)度為:c=1 mm,L=0.4 mm,探討電導(dǎo)率100~1 000 S/m,頻率1~18 GHz條件下,計(jì)算織物屏蔽效能,結(jié)果如圖4所示。

圖4 電導(dǎo)率、頻率與和屏蔽效能的關(guān)系Fig.4 Relations among conductivity,different frequency and shielding effectiveness

由圖4可知:當(dāng)電導(dǎo)率小于500 S/m時(shí),屏蔽效能隨頻率變化基本保持不變;但大于500 S/m,屏蔽效能隨頻率變大呈逐漸增大的趨勢(shì),且電導(dǎo)率越大其趨勢(shì)愈加顯著。這是由于電導(dǎo)率較小時(shí),電磁波在織物內(nèi)部產(chǎn)生的電場(chǎng)較弱,電場(chǎng)力較小,織物內(nèi)部自由電子受到約束力較大,難以發(fā)生定向移動(dòng),導(dǎo)致電導(dǎo)率小于500 S/m時(shí)織物屏蔽效能未呈現(xiàn)明顯頻散特性。為進(jìn)一步研究電導(dǎo)率對(duì)屏蔽效能的影響,分別提取圖4中電導(dǎo)率為100、400、700和1 000 S/m時(shí)屏蔽效能的平均值,結(jié)果如圖5所示。

圖5 電導(dǎo)率與屏蔽效能的關(guān)系Fig.5 Relations between conductivity and shielding effectiveness

由圖5可看出:織物結(jié)構(gòu)完全一致情況下,屏蔽效能隨電導(dǎo)率增加而顯著提高;電導(dǎo)率在100~1 000 S/m范圍內(nèi),電導(dǎo)率增加300 S/m,屏蔽效能的增量由5.02 dB減至1.03 dB,這說明屏蔽效能隨電導(dǎo)率的增加呈非線性遞減變化。趨膚深度公式為

式中:δ為趨膚深度m;μ0為真空中的磁導(dǎo)率1.257×10-6H/m;μr為相對(duì)磁導(dǎo)率;σ為材料的電導(dǎo)率,S/m;f為電磁波的頻率,Hz。

由上式可知,趨膚深度隨電導(dǎo)率增加呈非線性減小,使得電磁屏蔽效能出現(xiàn)非線性遞減現(xiàn)象,這與仿真結(jié)果與理論分析一致。為深入研究電導(dǎo)率對(duì)屏蔽效能的影響,分別取電導(dǎo)率為100S/m和1 000S/m時(shí)織物截面的感應(yīng)電流密度分布,如圖6所示。

圖6 不同電導(dǎo)率下織物截面的感應(yīng)電流密度Fig.6 Induced current density of fabric cross section at different conductivity

由圖6可知:電導(dǎo)率為100S/m時(shí),織物截面最大電流密度為2.281×105A/m2,自由電子無法定向移動(dòng);當(dāng)電導(dǎo)率為1 000S/m時(shí),織物內(nèi)部電子形成明顯順時(shí)針環(huán)流,最大電流密度達(dá)5.027×105A/m2。說明隨電導(dǎo)率增加,織物內(nèi)部自由電子受到的約束力會(huì)顯著下降,定向移動(dòng)的電子數(shù)量迅速增加,織物內(nèi)部反向電場(chǎng)強(qiáng)度增大,提高了織物的屏蔽效能;隨頻率變大,電磁波波長(zhǎng)縮短,在織物內(nèi)部反射次數(shù)增加,反射損耗的逐漸增大使織物屏蔽效能隨頻率變大而提高。

2.2 電磁波入射角對(duì)屏蔽效能的影響

選用組織單元長(zhǎng)度c=1 mm,厚度L=0.4 mm,電導(dǎo)率σ=1 000 S/m的織物為仿真樣品,探討織物電磁波入射角對(duì)屏蔽效能的影響,計(jì)算結(jié)果如圖7所示。

圖7 入射角、頻率與屏蔽效能的關(guān)系Fig.7 Relations among incidence angle,different frequency and shielding effectiveness

由圖7可知,屏蔽效能隨入射角增大略有提高。當(dāng)電磁波垂直入射,即入射角為0°時(shí),電磁屏蔽效能最小;入射角從0°增至60°時(shí),屏蔽效能增加了1.5 dB,且其增值與電磁波頻率無關(guān)。這是因?yàn)榭椢锏母袘?yīng)電流密度會(huì)隨著電磁波入射角的變大而增加,提高了織物屏蔽效能。為定量分析電磁波入射角變化對(duì)感應(yīng)電流的影響,選取在一定頻率(9.5 GHz)下,入射角為0°和60°時(shí)織物表面的感應(yīng)電流密度分布,結(jié)果如圖8所示。

圖8 不同入射角下織物表面的感應(yīng)電流密度Fig.8 Induced current density of fabric surface at different incidence angle

由圖8可看出:電磁波入射角為0°時(shí),織物表面電流分布均勻,且電流密度低,最大電流密度為3.036×106A/m2;而入射角為60°時(shí),最大電流密度增至7.195×106A/m2,相對(duì)0°時(shí)顯著增加且分布不均勻,具有極化效應(yīng)。

2.3 織物密度對(duì)屏蔽效能的影響

本文模型選用織物材料電導(dǎo)率為1 000 S/m,依然采用周期邊界條件,保持紗線參數(shù)一致,改變織物密度,設(shè)計(jì)5組仿真試樣,其參數(shù)如表1所示。通過CST有限積分計(jì)算,得到頻率1~18 GHz的電磁屏蔽效能,如圖9所示。

表1 織物規(guī)格參數(shù)表Tab.1 Fabric specification parameters

圖9 單元長(zhǎng)度、頻率與與屏蔽效能的關(guān)系Fig.9 Relations among cell length,different frequency and shielding effectiveness

由圖9可知,織物組織單元長(zhǎng)度短,即織物越緊密電磁屏蔽效能越大,且網(wǎng)格大小對(duì)屏蔽效能的頻散性能影響顯著。在相同頻率下,屏蔽效能隨單元長(zhǎng)度增加梯度減小,此結(jié)論與文獻(xiàn)[5]的實(shí)驗(yàn)結(jié)果具有很好的一致性。根據(jù)網(wǎng)孔分析理論,織物孔隙近似為矩形孔,其傳輸系數(shù)為

式中,n為單位面積網(wǎng)孔數(shù)。

由上式可知,c=1時(shí)織物較緊密,屏蔽效能較高,隨組織單元長(zhǎng)度增加網(wǎng)孔面積b2迅速增大,使傳輸系數(shù)增加而降低織物屏蔽效能,該結(jié)論與仿真結(jié)果吻合。

2.4 織物層數(shù)對(duì)屏蔽效能的影響

織物厚度是影響電磁屏蔽效能的重要因素之一,為對(duì)其進(jìn)行定量分析,選用平紋模型(參數(shù)為單元長(zhǎng)度c=1 mm,單元厚度L=0.4 mm,電導(dǎo)率為1 000 S/m)進(jìn)行層疊,對(duì)比單層、雙層及三層織物對(duì)屏蔽效能的影響,經(jīng)過仿真計(jì)算得到的結(jié)果如圖10所示。

圖10 織物層數(shù)、頻率與與屏蔽效能的關(guān)系Fig.10 Relations among layer count of fabric,different frequency and shielding effectiveness

由圖10可知,織物層數(shù)增加,屏蔽效能明顯增大,且隨層數(shù)增加其頻散特性愈加顯著。這是因?yàn)榭椢锖穸仍黾邮闺姶挪ㄔ诳椢飪?nèi)部的反射損耗顯著增大,降低電磁波透射率,提高織物屏蔽效能。選取頻率9.5 GHz時(shí)3層層疊織物的截面電流密度分布如圖11所示。從圖可知,電磁波入射方向的第1層織物內(nèi)部形成明顯電流,第2層、第3層的織物電流密度迅速降低。由于電磁波入射到第1層時(shí)所形成的反向電場(chǎng)場(chǎng)強(qiáng)大,電磁波大部分已被反射到自由空間,因此隨后第2層織物、第3層織物內(nèi)部電流密度顯著減小;電磁波單次反射損耗以及在織物內(nèi)部多次反射損耗,造成了織物內(nèi)部電流密度梯度性降低。

圖11 三層織物截面的感應(yīng)電流密度Fig.11 Induced current density at three layers fabric

3 結(jié) 論

基于電磁屏蔽織物的結(jié)構(gòu)微觀分析與建模,通過有限積分法對(duì)織物模型進(jìn)行電磁仿真,探討織物電導(dǎo)率、電磁波入射角、織物密度以及織物層數(shù)的理論計(jì)算與比較,可得出以下結(jié)論。

1)織物結(jié)構(gòu)三維模型可較好地刻畫出真實(shí)織物的結(jié)構(gòu)形態(tài),通過參數(shù)調(diào)整可適于對(duì)電磁屏蔽織物的建模,具有較好適用性。

2)織物電磁屏蔽效能隨材料電導(dǎo)率增加顯著提高,當(dāng)電導(dǎo)率大于500 S/m后屏蔽效能出現(xiàn)頻散特性,并隨電導(dǎo)率繼續(xù)增加其頻散特性愈加顯著;織物表面電流密度會(huì)隨著電磁波入射角增大而出現(xiàn)極化效應(yīng)。

3)織物密度是影響屏蔽效能的關(guān)鍵因素。織物層疊后第1層對(duì)電磁波反射損耗起主要作用,因此在設(shè)計(jì)開發(fā)屏蔽織物時(shí),織物密度的梯度性設(shè)計(jì)是調(diào)和屏蔽效能與服用性能的有效途徑。

FZXB

[1] OZDEMIR H,U URLU E S,OZKURT A.The electromagnetic shielding of textured steel yarn based woven fabrics used for clothing[J].Journal of Industrial Textiles,DOI: 10.1177/1528083715569369.

[2] PERUMALRAJ R,DASARADHAN B,NALANKILLI G.Copper,stainless steel,glass core yarn,and ply yarn woven fabric composite materials properties[J].Journal of Reinforced Plastics and Composites,2010,29(20):3074-3082.

[3] SHYR T W,SHIE J W.Electromagnetic shielding mechanisms using soft magnetic stainless steel fiber enabled polyester textiles[J].Journal of Magnetism and Magnetic Materials,2012,324(23):4127-4132.

[4] 肖紅,唐章宏,王群,等.電磁屏蔽織物的導(dǎo)電網(wǎng)格結(jié)構(gòu)及其屏蔽效能的一般影響規(guī)律研究[J].紡織學(xué)報(bào),2015,36(2):36-42.XIAO Hong,TANG Zhanghong,WANG Qun,et al.Research on conductive grid structure and general influence factors to shielding effectiveness of electromagnetic shielding fabrics[J].Journal of Textile Research,2015,36(2):36-42.

[5] XIAO H,SHI M W,WANG Q,et al.The electromagnetic shielding and reflective properties of electromagnetic textiles with pores,planar periodic units and space structures[J].Textile Research Journal,2014,84(16):1679-1691.

[6] DAS A,KOTHARI V K,KOTHARI A,et al.Effect of various parameters on electromagnetic shielding effectiveness of textile fabrics[J].Indian Journal of Fibre & Textile Research,2009,34:144-148.

[7] 馮亞麗,張慧茹.滌綸織物化學(xué)鍍銀工藝探究[J].高分子材料科學(xué)與工程,2014,30(6):106-110.FENG Yali,ZHANG Huiru.Electroless silver plating on polyester fabric substrate[J].Polymer Materials Science and Engineering,2014,30(6):106-110.

[8] 徐文龍,熊杰,徐勤.化學(xué)鍍銀滌綸織物的制備及其性能[J].紡織學(xué)報(bào),2011,32(9):42-52.XU Wenlong,XIONG Jie,XU Qin.Preparation and properties of silver electroless plated PET fabrics[J].Journal of Textile Research,2011,32(9):42-52.

[9] 孟靈靈,黃新民,魏取福.滌綸基布磁控濺射銅膜及其電磁屏蔽性能[J].印染,2013(16):1-5.MENG Lingling,HUANG Xinmin,WEI Qufu.Structure of copper films on polyester fabrics by magnetron sputtering and its electromagnetic shielding properties[J].China Dyeing & Printing,2013(16):1-5.

[10] CHEN H C,LEE K C,LIN J H,et al.Fabrication of conductive woven fabric and analysis of electromagnetic shielding via measurement and empirical equation[J].Journal of Materials Processing Technology,2007,184(1/2/3):124-30.

[12] SARTO M S,GRECO S,TAMBURRANO A.Shielding Effectiveness of Protective Metallic Wire Meshes: EMModeling and Validation[J].Ieee Transactions on Electromagnetic Compatibility,2014,56(3):615-621.

[13] 李奇軍,劉長(zhǎng)隆,周明,等.含不銹鋼纖維機(jī)織物電磁屏蔽機(jī)理及其效能的仿真研究[J].功能材料,2013,44(14):2041-2046.LI Qijun,LIU Changlong,ZHOU Ming,et al.Study on principle of the shielding effectiveness of woven fabrics containing stainless steel fibers and simula-tion[J].Function Material,2013,44(14):2041-2046.

[14] LIU Z,WANG X C.Influence of fabric weave type on the effectiveness of electromagnetic shielding woven fabric[J].Journal of Electromagnetic Waves and Applications,2012,26(14):1848-1856.

[15] WONG C,LONG A,SHERBURN M,et al.Comparisons of novel and efficient approaches for permeability prediction based on the fabric architec-ture[J].Composites Part A: Applied Science and Manufacturing,2006,37(6):847-857.

Simulation analysis of woven fabric electromagnetic shielding effectiveness using finite integration technique

SU Qincheng1,ZHAO Xiaoming1,LI Weibin1,LI Jianxiong2

(1.SchoolofTextile,TianjinPolytechnicUniversity,Tianjin300387,China;2.SchoolofElectronicsandInformationEngineering,TianjinPolytechnicUniversity,Tianjin300387,China)

In order to study the shielding mechanism of fabric,and provide theoretical guidance for the development of electromagnetic shielding fabrics.The microstructure of the fabric has been modeled,and a three-dimensional electromagnetic simulation model of the fabric was established by CST microwave studio.The finite integration technique was adopted to compute the shielding effectiveness in 1-18 GHz range,including different conductivities,densities and layer counts of the fabric,and different incidence angles of electromagnetic wave.The result shows that the conductivity is the key factor for the fabric.The shielding effectiveness as the fabric structure parameters change is nonlinear.The electromagnetic wave incidence angle increase would cause the polarization effect of surface current density on the fabric.This result has universal applicability,and provides a theoretical basis for optimizing the design for high-performance shielding fabrics.

electromagnetic shielding fabric; microstructure model; shielding effectiveness; electromagnetic simulation

10.13475/j.fzxb.20150603106

2015-06-14

2015-11-02

國(guó)家自然科學(xué)基金面上項(xiàng)目(61372011)

蘇欽城(1990—),男,碩士生。主要從事電磁紡織品數(shù)值模擬方面的研究。趙曉明,通信作者,E-mail: tex_zhao@163.com。

TP 391.9;TS 101.8

A

猜你喜歡
結(jié)構(gòu)模型
一半模型
《形而上學(xué)》△卷的結(jié)構(gòu)和位置
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
論結(jié)構(gòu)
新型平衡塊結(jié)構(gòu)的應(yīng)用
模具制造(2019年3期)2019-06-06 02:10:54
論《日出》的結(jié)構(gòu)
3D打印中的模型分割與打包
FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長(zhǎng)
主站蜘蛛池模板: av一区二区无码在线| 国产综合另类小说色区色噜噜| 欧美啪啪网| 激情无码视频在线看| 日本欧美中文字幕精品亚洲| 国产第一福利影院| 韩国v欧美v亚洲v日本v| 夜夜操国产| 97视频免费在线观看| 中文字幕无线码一区| 国产午夜一级毛片| 国产手机在线小视频免费观看| 久久美女精品| 免费aa毛片| 亚洲人成网站在线播放2019| 亚洲一区二区日韩欧美gif| 国产精品自在在线午夜| 午夜视频在线观看免费网站| 亚洲精品不卡午夜精品| 亚洲精品第一页不卡| 欧美午夜理伦三级在线观看| 亚洲愉拍一区二区精品| 国产精品护士| 中文字幕人妻无码系列第三区| 手机在线免费不卡一区二| 国产农村精品一级毛片视频| 一级做a爰片久久免费| 精品1区2区3区| 亚洲天堂视频在线观看免费| 999精品色在线观看| 日韩精品成人在线| 久久久久夜色精品波多野结衣| 国产免费a级片| 四虎永久在线视频| 久久久国产精品免费视频| 国产理论一区| 99热最新网址| 免费观看男人免费桶女人视频| 欧美激情首页| 亚洲综合狠狠| 在线99视频| 麻豆AV网站免费进入| 无码中文字幕乱码免费2| 久久香蕉国产线| www.亚洲一区| 国产三级成人| 国产91无码福利在线| 国产精品亚洲日韩AⅤ在线观看| 毛片大全免费观看| 欧美日韩国产在线观看一区二区三区 | 欧美亚洲日韩不卡在线在线观看| 国产无码性爱一区二区三区| 免费看a级毛片| 国产成人一区在线播放| 欧美在线网| 成人免费视频一区| 亚洲人成影视在线观看| 国内熟女少妇一线天| 国产亚洲欧美另类一区二区| 亚洲狼网站狼狼鲁亚洲下载| 久久这里只有精品23| аⅴ资源中文在线天堂| 亚洲欧美自拍中文| 欧美福利在线观看| 亚洲乱伦视频| 91九色视频网| 久久精品娱乐亚洲领先| 老色鬼欧美精品| 最近最新中文字幕在线第一页| 国产91成人| 国产美女久久久久不卡| 日韩成人在线网站| 一级片一区| 伊人久久婷婷| 91蜜芽尤物福利在线观看| 国产成人夜色91| www.国产福利| 无码综合天天久久综合网| 欧美激情第一区| 欧美日本激情| 91在线激情在线观看| 久久特级毛片|