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

    基于SNS-FEM的電動汽車電纜線電磁輻射仿真

    2023-09-04 14:36:04黃湛勇安玉民
    計(jì)算機(jī)仿真 2023年7期
    關(guān)鍵詞:電纜線四面體感應(yīng)電流

    黃湛勇,王 剛,安玉民

    (1. 河北工業(yè)大學(xué)機(jī)械工程學(xué)院,天津 300401;2. 天津市新能源汽車動力傳動與安全技術(shù)重點(diǎn)實(shí)驗(yàn)室,天津 300130;3. 國家技術(shù)創(chuàng)新方法與實(shí)施工具工程技術(shù)研究中心,天津 300401)

    1 引言

    電動汽車車內(nèi)電子設(shè)備眾多,使得電磁環(huán)境更為復(fù)雜。干擾電磁場不僅會使設(shè)備之間相互干涉產(chǎn)生故障,影響行車安全,還會對司乘人員的健康構(gòu)成威脅[1,2]。某電動汽車動力總成布置如圖1所示,動力電池組、電機(jī)控制器位于后方區(qū)域,驅(qū)動電機(jī)居于車身前端。在汽車行駛時(shí),電池組將電能提供給電機(jī)控制器,再通過動力電纜向驅(qū)動電機(jī)供電。受車身尺寸限制,駕駛員距動力電纜線較近,極易受到電纜線低頻輻射電磁場的影響。

    圖1 某電動汽車動力總成布置

    精確高效地預(yù)測電動汽車車內(nèi)干擾電磁場的分布,是電磁屏蔽設(shè)計(jì)的關(guān)鍵。目前的預(yù)測方法主要可分為兩類,即實(shí)驗(yàn)測量與數(shù)值仿真。在實(shí)驗(yàn)測量方面,李志強(qiáng)等[3]采用手持式頻譜分析儀測量確定了電動客車電磁輻射超標(biāo)元器件。劉詳武[4]對車載中控控制模塊進(jìn)行輻射發(fā)射測試,并根據(jù)測試結(jié)果提出了屏蔽建議。黃青云[5]運(yùn)用實(shí)驗(yàn)測量的方式對某電動汽車驅(qū)動電機(jī)產(chǎn)生的電磁場分布進(jìn)行了探究。然而,在研究人體暴露于電纜線低頻電磁場問題時(shí),實(shí)驗(yàn)測量的方式雖容易測得電纜線周圍空間的電磁場強(qiáng)度,但受限于測量方式本身,難以測量人體上的感應(yīng)場。

    數(shù)值仿真方法通過抽象物理建模,擺脫了對實(shí)驗(yàn)設(shè)備及測量位置的限制,彌補(bǔ)了實(shí)驗(yàn)方法的不足。常用的數(shù)值仿真方法包括有限元法、有限差分法(FTDT)、矩量法(MOM)等[6],這其中以有限元法的應(yīng)用最為廣泛。董旭偉等[7]應(yīng)用COMSOL有限元軟件分析了動力電纜線對不同乘坐位置兒童的輻射影響。杜明磊等[8]通過建立車身及動力線纜的電磁輻射模型,研究了線纜長度、端接電阻等因素對線纜輻射強(qiáng)度的影響。陳國強(qiáng)等[9]運(yùn)用矩量法和傳輸線法對驅(qū)動電機(jī)母線的輻射電磁場分布規(guī)律進(jìn)行了探究。然而,應(yīng)用FEM對低頻電磁場人體暴露問題進(jìn)行研究時(shí),仍存在一些不足之處。如采用六面體單元進(jìn)行離散可以提高FEM求解精度,但人體結(jié)構(gòu)復(fù)雜,六面體單元離散會增加網(wǎng)格劃分負(fù)擔(dān),況且車內(nèi)空氣域同樣需要進(jìn)行離散,每個(gè)單元都需要高階插值和多個(gè)積分點(diǎn),導(dǎo)致計(jì)算量巨大。四面體單元可以提升網(wǎng)格劃分效率,降低模型邊界離散誤差,但四面體單元較低的插值階次限制了其求解精度。更重要的是,傳統(tǒng)有限元法由于本身系統(tǒng)剛度“過硬”的缺陷[6],使得其在電磁場數(shù)值計(jì)算中的求解精度并不高。

    適當(dāng)“軟化”離散系統(tǒng)的剛度是提升電磁場數(shù)值仿真精度與效率的關(guān)鍵。受無網(wǎng)格法中應(yīng)變光滑思想的啟發(fā),Liu等提出了G空間與弱-弱形式(W2)的概念[10,11],并進(jìn)一步構(gòu)造了光滑有限元法的求解框架(smoothed finite element methods,S-FEM)[12,13]。其中,基于節(jié)點(diǎn)的光滑有限元模型(node-based smoothed finite element method,NS-FEM)通過“軟化”離散系統(tǒng)的剛度能夠給出真實(shí)解的“上界”[14],在固體力學(xué)問題分析中取得了很好的效果[15]。然而,NS-FEM在時(shí)間域內(nèi)是不穩(wěn)定的,這限制了其在電磁計(jì)算領(lǐng)域的應(yīng)用[16]。為克服該缺陷,Cui等通過泰勒展開補(bǔ)充穩(wěn)定項(xiàng),進(jìn)而提出了穩(wěn)定節(jié)點(diǎn)積分算法(SNS-FEM)[17,18],在不顯著增加計(jì)算量的前提下提升了求解精度,有效地避免了NS-FEM的時(shí)域不穩(wěn)定性,并在聲學(xué)與靜力學(xué)領(lǐng)域得到了廣泛應(yīng)用[19,20]。

    為提高電動汽車電纜線低頻電磁場預(yù)測的精度與效率,本文進(jìn)一步將梯度光滑操作[11](GST)與穩(wěn)定項(xiàng)增強(qiáng)引入到電纜線電磁輻射計(jì)算的控制方程中,構(gòu)建了電纜線對人體輻射的穩(wěn)定節(jié)點(diǎn)積分方法。對車身-電纜線-人體建立幾何建模,采用四面體單元進(jìn)行離散,在背景網(wǎng)格上構(gòu)造點(diǎn)光滑域,基于GST對電勢梯度與磁感應(yīng)強(qiáng)度進(jìn)行光滑處理,通過泰勒展開構(gòu)造積分增強(qiáng)項(xiàng),依托光滑Galerkin弱形式獲得電磁系統(tǒng)離散方程?;谒鶚?gòu)造的方法預(yù)測了人體感應(yīng)電流密度分布,結(jié)果表明:SNS-FEM相比于FEM具有更高的預(yù)測精度和更快的求解效率。

    2 交變電磁場基本理論

    電動汽車動力電纜線通過低頻電流時(shí)會在周圍空間產(chǎn)生呈交流正弦變化的電磁場,并在人體中產(chǎn)生感應(yīng)電渦流,可用如下的麥克斯韋方程來表征

    ?×H=J

    (1)

    (2)

    ?·B=0

    (3)

    式中,H為磁場強(qiáng)度,E為電場強(qiáng)度,B為磁感應(yīng)強(qiáng)度,J為電流密度。

    引入磁阻率ν,電導(dǎo)率σ,矢量磁位A,標(biāo)量電位φ,穩(wěn)態(tài)電渦流問題控制方程可寫為

    (4)

    (5)

    對于本文考慮的電渦流問題,邊界條件分為三類,即Γ=Γ1∪Γ2∪Γ3。其中,Γ1為狄利克雷邊界條件,Γ2為紐曼邊界條件,Γ3為磁導(dǎo)體邊界條件,相關(guān)約束形式如下

    3 電磁場穩(wěn)定節(jié)點(diǎn)積分算法

    3.1 光滑域構(gòu)造與梯度光滑操作

    (7)

    (8)

    在進(jìn)行梯度光滑操作之前,需要基于背景網(wǎng)格進(jìn)一步構(gòu)造光滑域。如圖2所示,圍繞節(jié)點(diǎn)k所形成的光滑域可通過順次連接四面體單元的中心、表面三角形的形心、網(wǎng)格邊中點(diǎn)來構(gòu)建。整個(gè)問題域?qū)⒒诠?jié)點(diǎn)形成Nn個(gè)光滑域,不同光滑域之間滿足Ω=Ω1∪Ω2∪…∪ΩNn,Ωi∪Ωj=?,(i≠j,i=1,2,3…,Nn,j=1,2,3…,Nn)。在光滑域內(nèi),依托梯度光滑技術(shù)[11]對磁感應(yīng)強(qiáng)度與電勢的梯度進(jìn)行光滑處理

    圖2 基于節(jié)點(diǎn)k的光滑域

    (9)

    (10)

    (11)

    將式(7)、(8)、(11)代入(9)、(10)中,可得如下形式的光滑梯度項(xiàng)

    (12)

    (13)

    (14)

    基于光滑Galerkin弱形式[11],穩(wěn)態(tài)電渦流問題控制方程可描述為

    (15)

    經(jīng)整理可得如下矩陣形式

    (16)

    (17)

    上述方式構(gòu)建的電渦流方程在空間域是穩(wěn)定的而在時(shí)間域內(nèi)是不穩(wěn)定的,因此不能直接用于電磁問題地求解。

    3.2 近似光滑域及穩(wěn)定項(xiàng)構(gòu)造

    圖3 近似光滑域的構(gòu)造

    (18)

    (19)

    基于式(18)、(19),圖3虛擬積分點(diǎn)1、2處場函數(shù)梯度值分別為

    (20)

    (21)

    (22)

    (23)

    4 算例分析

    依托上述理論,基于Fortran語言開發(fā)了低頻電磁場問題求解的SNS-FEM源程序代碼。通過如下兩個(gè)電纜線低頻電磁場數(shù)值算例來驗(yàn)證所構(gòu)造方法在精度、收斂性等方面的優(yōu)勢。

    4.1 電纜線低頻電磁場輻射問題

    考慮圖4所示的電纜線輻射模型,長方體空氣域內(nèi)穿過一段電纜線,旁邊放置一矩形體,為便于描述說明,在矩形體上定義P點(diǎn),路徑a與路徑b。

    圖4 電纜線輻射模型

    采用4節(jié)點(diǎn)四面體單元對問題域進(jìn)行離散,單元總數(shù)為227234,節(jié)點(diǎn)總數(shù)為37045??紤]到此問題無解析解,這里采用COMSOL軟件中節(jié)點(diǎn)間距極小的有限元解作為參考解(1032645個(gè)四面體單元,168620個(gè)節(jié)點(diǎn))。

    圖5(a)和(b)分別給出了FEM和SNS-FEM計(jì)算得到的矩形體側(cè)面與底面的感應(yīng)電流密度分布,圖5(c)為參考解??梢钥闯?①對于底面感應(yīng)電流密度的分布,兩種方法所得結(jié)果均接近參考解。但在中線附近,FEM由于系統(tǒng)剛度“過硬”的缺陷,所得結(jié)果與參考解仍存在細(xì)微偏差;②對于側(cè)面的感應(yīng)電流密度分布特別是在靠近電纜線的區(qū)域,FEM數(shù)值解與參考解存在較大程度的偏差,反觀SNS-FEM,由于進(jìn)行了梯度光滑操作和剛度修正,使其在整個(gè)區(qū)域上與參考解基本保持一致。上述結(jié)果反映了SNS-FEM求解電磁場時(shí)精度方面的優(yōu)勢,同時(shí)也驗(yàn)證了所構(gòu)造方法的正確性。

    圖5 矩形體感應(yīng)電流密度分布

    為定量說明所提方法計(jì)算誤差,圖6、圖7分別給出了兩種方法在路徑a與路徑b上計(jì)算的感應(yīng)電流密度分布??梢园l(fā)現(xiàn):①在靠近電纜線的位置,FEM得到了偏小的感應(yīng)電流密度值,與參考解差別較大。如在a上z=0.1m處,與參考解存在近10%的偏差;②SNS-FEM可以提供與參考解十分一致的數(shù)值結(jié)果,吻合程度高。

    圖6 路徑a上感應(yīng)電流密度值

    圖7 路徑b上感應(yīng)電流密度值

    為進(jìn)一步驗(yàn)證SNS-FEM的高效性,取圖4的P點(diǎn)(x=1.3m)的感應(yīng)電流密度進(jìn)行說明。采用不同尺寸的網(wǎng)格對問題域離散,并引入相對誤差Er對計(jì)算結(jié)果進(jìn)行度量。圖8給出了FEM與SNS-FEM兩種方法的收斂特性比較,可以看出:①就收斂速度而言,SNS-FEM收斂速度優(yōu)于FEM;②在相同精度下SNS-FEM所需的自由度更少,如Er=4%時(shí),SNS-FEM所需的自由度僅為FEM的60%。上述分析表明SNS-FEM能有效減少前期網(wǎng)格劃分負(fù)擔(dān),提升建模效率。

    圖8 兩種方法收斂特性比較

    4.2 電纜線低頻電磁場對真實(shí)人體輻射問題

    為進(jìn)一步探究SNS-FEM處理電動汽車電磁場實(shí)際問題的能力,本節(jié)構(gòu)造了真實(shí)人體暴露于電纜線低頻電磁場的數(shù)值模型。如圖9所示,按照國際通用1.75m成年人身高比例建立人體模型,并選用人體駕駛姿態(tài)。過于復(fù)雜的車體外殼會對建模效率產(chǎn)生不利影響,因此這里對實(shí)際模型進(jìn)行了一定程度的簡化,僅留下電動汽車司乘艙、貫穿車體的動力電纜線,同時(shí)保留前擋風(fēng)玻璃車窗等幾何特征。

    圖9 電纜線對真實(shí)人體電磁輻射模型

    將圖9問題域離散成391800個(gè)4節(jié)點(diǎn)四面體單元,節(jié)點(diǎn)總數(shù)為63106??紤]到此問題無解析解,同樣,在COMSOL軟件中采用節(jié)點(diǎn)間距極小的有限元解作為參考解(770605個(gè)四面體單元,128168個(gè)節(jié)點(diǎn))。

    圖10(a)和(b)分別給出了FEM和SNS-FEM計(jì)算得到的人體軀干表面感應(yīng)電流密度分布,圖10(c)為參考解。為便于說明,這里取腳踝處的感應(yīng)電流密度分布進(jìn)行放大。對比可發(fā)現(xiàn):①FEM在右腿腳踝處得到的感應(yīng)電流密度分布與參考解存在較大偏差,幾乎完全偏離了參考解;FEM在左腿腳踝的所得結(jié)果與參考解也存在較大出入,由此可見FEM對于復(fù)雜模型適用性較差;②SNS-FEM得到的數(shù)值結(jié)果精度較好,無論在左腿或右腿處均能獲得與參考解較為一致的感應(yīng)電流密度分布。這是由于梯度光滑技術(shù)與穩(wěn)定項(xiàng)修正對系統(tǒng)剛度進(jìn)行了適當(dāng)?shù)摹败浕?從而降低數(shù)值模型的離散誤差,帶來計(jì)算精度的大幅提高。

    為定量分析兩種方法在求解復(fù)雜低頻電磁場時(shí)的誤差情況,在人體靠近電纜側(cè)的左腿上定義路徑1與路徑2,位置如圖11所示。圖12和圖13分別給出了不同數(shù)值方法計(jì)算得到的路徑1和路徑2上的感應(yīng)電流密度分布??梢钥闯?①FEM計(jì)算結(jié)果在峰值處普遍偏小,在兩條路徑上與參考解均存在偏離,特別是在路徑2上z=0.11~0.15(m)區(qū)域,其結(jié)果與參考解偏離程度較大;②SNS-FEM在路徑1上的數(shù)值解與參考解高度吻合。路徑2上,SNS-FEM雖然在z=0.11~0.15(m)區(qū)域與參考解存在稍許偏離,但其結(jié)果依然優(yōu)于FEM。除此之外,在其余區(qū)域上SNS-FEM均能給出與參考解非常一致的結(jié)果。由此可以認(rèn)為,相比FEM,SNS-FEM更加適用于電動汽車電磁場問題的數(shù)值計(jì)算。

    圖11 路徑1、路徑2位置

    圖12 路徑1上感應(yīng)電流密度值

    圖13 路徑2上感應(yīng)電流密度值

    5 結(jié)論

    針對FEM分析電動汽車電纜線低頻電磁場時(shí)存在的計(jì)算精度低、效率差等數(shù)值缺陷,構(gòu)造了電纜線低頻電磁場仿真的SNS-FEM數(shù)值模型,理論分析與數(shù)值算例表明:

    1)SNS-FEM可以用常規(guī)4節(jié)點(diǎn)四面體單元對問題域進(jìn)行離散,有效突破低階四面體單元對FEM求解精度的限制,減輕了前期網(wǎng)格劃分負(fù)擔(dān),提升了建模效率。

    2)SNS-FEM由于在FEM中引入梯度光滑操作與穩(wěn)定項(xiàng)增強(qiáng),使得在預(yù)測感應(yīng)場時(shí),峰值處精度比FEM提升10%,達(dá)到相同精度所需自由度約為FEM的60%,不用過多加密網(wǎng)格也可獲得十分精準(zhǔn)的結(jié)果。

    3)SNS-FEM對具有復(fù)雜幾何形狀的電磁場工程問題仍然可以提供非常精確的數(shù)值解,表明SNS-FEM在實(shí)際工程應(yīng)用中具有較大潛力。

    猜你喜歡
    電纜線四面體感應(yīng)電流
    楞次定律推論的巧妙應(yīng)用
    四面體小把戲
    金屬軌道的形狀影響感應(yīng)電流的大小嗎
    第23和24太陽活動周高緯地磁感應(yīng)電流分布特性
    R3中四面體的幾個(gè)新Bonnesen型不等式
    R3中四面體的Bonnesen型等周不等式
    一種新型光電纜線敷設(shè)車研究
    通信傳輸中信號衰減的原因及對策
    通信直流電源系統(tǒng)電纜線徑的選擇方式
    找故障
    亚洲精品国产一区二区精华液| 人人妻人人添人人爽欧美一区卜| 久久午夜综合久久蜜桃| av天堂久久9| 啦啦啦中文免费视频观看日本| 在线十欧美十亚洲十日本专区| 老熟妇仑乱视频hdxx| 精品视频人人做人人爽| 亚洲一码二码三码区别大吗| 久久久久久亚洲精品国产蜜桃av| 夜夜骑夜夜射夜夜干| 成年av动漫网址| 在线精品无人区一区二区三| 欧美日韩视频精品一区| 日日爽夜夜爽网站| 老熟妇仑乱视频hdxx| 亚洲国产精品999| 一级,二级,三级黄色视频| 国产成人欧美在线观看 | 在线观看免费日韩欧美大片| 国产成人啪精品午夜网站| 一本大道久久a久久精品| 丝瓜视频免费看黄片| 欧美日本中文国产一区发布| 女警被强在线播放| 99久久综合免费| 每晚都被弄得嗷嗷叫到高潮| 国产欧美亚洲国产| 精品国产一区二区久久| 免费一级毛片在线播放高清视频 | 亚洲精品美女久久av网站| 桃红色精品国产亚洲av| 人人妻人人添人人爽欧美一区卜| 动漫黄色视频在线观看| 99精品欧美一区二区三区四区| 男女边摸边吃奶| 熟女少妇亚洲综合色aaa.| 久久精品亚洲熟妇少妇任你| 波多野结衣av一区二区av| 亚洲人成77777在线视频| 美女脱内裤让男人舔精品视频| 亚洲人成77777在线视频| 日韩欧美免费精品| 成在线人永久免费视频| 精品免费久久久久久久清纯 | 亚洲国产成人一精品久久久| 亚洲成人免费av在线播放| 精品视频人人做人人爽| 欧美日本中文国产一区发布| 欧美日韩亚洲综合一区二区三区_| 交换朋友夫妻互换小说| 波多野结衣一区麻豆| 亚洲人成电影免费在线| 日本wwww免费看| av片东京热男人的天堂| 波多野结衣一区麻豆| 国产精品亚洲av一区麻豆| 91麻豆av在线| 成人18禁高潮啪啪吃奶动态图| 精品国内亚洲2022精品成人 | 成年人免费黄色播放视频| 91av网站免费观看| 十八禁人妻一区二区| 久久久精品国产亚洲av高清涩受| 一本综合久久免费| 欧美97在线视频| 亚洲国产精品一区二区三区在线| 欧美国产精品一级二级三级| 超碰97精品在线观看| 国产精品二区激情视频| 两性夫妻黄色片| 国产精品一区二区免费欧美 | 成人国产av品久久久| 亚洲欧美日韩另类电影网站| 国产在线视频一区二区| 99热网站在线观看| 黄色视频在线播放观看不卡| 深夜精品福利| 人人妻人人澡人人看| 真人做人爱边吃奶动态| 韩国精品一区二区三区| 亚洲人成77777在线视频| 两人在一起打扑克的视频| 午夜91福利影院| 男女午夜视频在线观看| 国产成人精品久久二区二区91| 亚洲国产欧美在线一区| 建设人人有责人人尽责人人享有的| kizo精华| 日韩有码中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲欧美精品自产自拍| 久久久久精品国产欧美久久久 | 狠狠狠狠99中文字幕| 国产深夜福利视频在线观看| 亚洲国产日韩一区二区| 永久免费av网站大全| 亚洲国产欧美一区二区综合| 免费高清在线观看视频在线观看| 亚洲欧美成人综合另类久久久| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧洲日产国产| 国产1区2区3区精品| 9热在线视频观看99| 丝瓜视频免费看黄片| 亚洲国产精品999| 69av精品久久久久久 | 国产一区二区在线观看av| 中国国产av一级| 免费在线观看视频国产中文字幕亚洲 | 99国产精品免费福利视频| 在线观看免费日韩欧美大片| 亚洲精品av麻豆狂野| 日韩制服丝袜自拍偷拍| 永久免费av网站大全| 久久久精品94久久精品| 激情视频va一区二区三区| 十八禁高潮呻吟视频| 国产av精品麻豆| 中国美女看黄片| 欧美少妇被猛烈插入视频| 国产精品国产三级国产专区5o| 99精品欧美一区二区三区四区| 国产欧美日韩综合在线一区二区| 免费日韩欧美在线观看| 亚洲国产日韩一区二区| 午夜激情久久久久久久| 狂野欧美激情性bbbbbb| 国产三级黄色录像| 亚洲第一欧美日韩一区二区三区 | 国产日韩欧美视频二区| 免费高清在线观看日韩| 丰满人妻熟妇乱又伦精品不卡| 性色av一级| 国产免费现黄频在线看| 欧美另类一区| 久久人人97超碰香蕉20202| 香蕉丝袜av| 亚洲精品国产一区二区精华液| 亚洲三区欧美一区| 99精国产麻豆久久婷婷| av天堂久久9| 在线亚洲精品国产二区图片欧美| 天堂中文最新版在线下载| svipshipincom国产片| 动漫黄色视频在线观看| 免费女性裸体啪啪无遮挡网站| 欧美另类一区| 欧美 日韩 精品 国产| 黑丝袜美女国产一区| 国产男人的电影天堂91| 亚洲成人免费电影在线观看| 香蕉丝袜av| av一本久久久久| 国产一级毛片在线| 国产免费福利视频在线观看| 国产99久久九九免费精品| 新久久久久国产一级毛片| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精华国产精华精| 精品久久蜜臀av无| 在线观看舔阴道视频| 少妇粗大呻吟视频| 国产不卡av网站在线观看| 国产精品久久久久久人妻精品电影 | 亚洲av欧美aⅴ国产| 亚洲国产欧美在线一区| 丰满少妇做爰视频| av天堂在线播放| 视频区欧美日本亚洲| 别揉我奶头~嗯~啊~动态视频 | 国产精品自产拍在线观看55亚洲 | 国产精品亚洲av一区麻豆| 久久久久国产精品人妻一区二区| 日韩熟女老妇一区二区性免费视频| 国产男人的电影天堂91| 亚洲专区中文字幕在线| 亚洲激情五月婷婷啪啪| 午夜91福利影院| 国产免费av片在线观看野外av| 日韩有码中文字幕| 女性被躁到高潮视频| 中文字幕人妻丝袜一区二区| 亚洲avbb在线观看| 亚洲伊人色综图| 男女边摸边吃奶| 精品高清国产在线一区| 久久中文看片网| 国产精品秋霞免费鲁丝片| 色老头精品视频在线观看| 亚洲av美国av| 一级片'在线观看视频| 国产亚洲一区二区精品| 久久精品国产亚洲av高清一级| 亚洲情色 制服丝袜| av在线app专区| 中亚洲国语对白在线视频| 国产1区2区3区精品| 亚洲精品久久成人aⅴ小说| 欧美日韩一级在线毛片| 侵犯人妻中文字幕一二三四区| 久久久国产精品麻豆| 2018国产大陆天天弄谢| 99国产极品粉嫩在线观看| 亚洲精华国产精华精| 精品少妇内射三级| 精品人妻在线不人妻| 日本一区二区免费在线视频| 秋霞在线观看毛片| 男女之事视频高清在线观看| 黑人欧美特级aaaaaa片| 精品视频人人做人人爽| 女人精品久久久久毛片| 黄色 视频免费看| 国产成人精品久久二区二区91| 国产99久久九九免费精品| 精品国产一区二区三区四区第35| 在线观看免费午夜福利视频| 亚洲午夜精品一区,二区,三区| 美女高潮到喷水免费观看| 欧美中文综合在线视频| 国产高清视频在线播放一区 | 免费在线观看影片大全网站| 这个男人来自地球电影免费观看| 热99国产精品久久久久久7| 国产伦理片在线播放av一区| 久久久国产一区二区| 91麻豆精品激情在线观看国产 | 国产日韩欧美在线精品| 久久精品久久久久久噜噜老黄| 亚洲午夜精品一区,二区,三区| 国产亚洲欧美在线一区二区| 青春草亚洲视频在线观看| 丝袜喷水一区| bbb黄色大片| 97人妻天天添夜夜摸| 人人妻人人澡人人爽人人夜夜| 国产精品欧美亚洲77777| 好男人电影高清在线观看| 精品福利永久在线观看| 一区二区三区乱码不卡18| 国产精品二区激情视频| 国产高清videossex| 乱人伦中国视频| 久久久久国内视频| 欧美激情久久久久久爽电影 | 操美女的视频在线观看| 女人被躁到高潮嗷嗷叫费观| 91国产中文字幕| 欧美 亚洲 国产 日韩一| 欧美日本中文国产一区发布| 久久久精品94久久精品| avwww免费| 人人妻人人澡人人看| 亚洲中文字幕日韩| 久久免费观看电影| 老司机亚洲免费影院| 亚洲中文av在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲伊人色综图| 国产色视频综合| 亚洲国产精品成人久久小说| 欧美xxⅹ黑人| 欧美乱码精品一区二区三区| 国产一区二区激情短视频 | 国产三级黄色录像| 黄色怎么调成土黄色| 美女中出高潮动态图| 老司机午夜十八禁免费视频| 亚洲专区中文字幕在线| 免费观看a级毛片全部| 久热爱精品视频在线9| 国产伦人伦偷精品视频| 搡老熟女国产l中国老女人| 99国产精品免费福利视频| 国产成人啪精品午夜网站| 国产日韩欧美亚洲二区| 日本一区二区免费在线视频| 午夜视频精品福利| 十八禁人妻一区二区| 久久久久网色| 国产伦人伦偷精品视频| 在线观看免费高清a一片| kizo精华| 国产精品免费大片| 一个人免费看片子| 国产97色在线日韩免费| 极品少妇高潮喷水抽搐| 亚洲伊人久久精品综合| 欧美+亚洲+日韩+国产| 久久精品亚洲av国产电影网| 又紧又爽又黄一区二区| 久久天躁狠狠躁夜夜2o2o| 亚洲精品国产av成人精品| 日韩,欧美,国产一区二区三区| 久久ye,这里只有精品| 国产xxxxx性猛交| 黄色片一级片一级黄色片| 亚洲熟女毛片儿| 99国产精品一区二区蜜桃av | 少妇猛男粗大的猛烈进出视频| 日韩欧美一区二区三区在线观看 | 最近最新免费中文字幕在线| 国产精品久久久久久精品电影小说| 久久中文看片网| 国产国语露脸激情在线看| 国产日韩一区二区三区精品不卡| 一二三四在线观看免费中文在| 老汉色av国产亚洲站长工具| 搡老乐熟女国产| 精品国产一区二区三区四区第35| 91麻豆av在线| 日韩视频在线欧美| 欧美一级毛片孕妇| 男人舔女人的私密视频| 人人妻人人澡人人看| 12—13女人毛片做爰片一| 国产黄色免费在线视频| 国产日韩欧美视频二区| 黄片大片在线免费观看| 性色av一级| 又大又爽又粗| 黄网站色视频无遮挡免费观看| av超薄肉色丝袜交足视频| 亚洲欧美一区二区三区久久| 一区二区日韩欧美中文字幕| 十八禁网站免费在线| 久久天堂一区二区三区四区| 视频区图区小说| 在线 av 中文字幕| 国产欧美亚洲国产| 国产亚洲欧美精品永久| 建设人人有责人人尽责人人享有的| 欧美黑人欧美精品刺激| 中文精品一卡2卡3卡4更新| 99久久人妻综合| 18禁黄网站禁片午夜丰满| 不卡一级毛片| 在线观看免费日韩欧美大片| 999久久久精品免费观看国产| 超色免费av| av免费在线观看网站| 欧美精品一区二区大全| 亚洲欧美精品综合一区二区三区| 一级片免费观看大全| 一边摸一边抽搐一进一出视频| 欧美xxⅹ黑人| 老熟妇乱子伦视频在线观看 | 美女扒开内裤让男人捅视频| 精品一品国产午夜福利视频| 久久久精品94久久精品| 午夜福利一区二区在线看| 国产区一区二久久| 无遮挡黄片免费观看| 日本一区二区免费在线视频| 欧美 亚洲 国产 日韩一| 亚洲av美国av| 精品久久久久久电影网| 中文字幕制服av| 亚洲熟女精品中文字幕| 少妇粗大呻吟视频| 三级毛片av免费| www.自偷自拍.com| 女人精品久久久久毛片| 久久国产精品男人的天堂亚洲| 在线av久久热| 国产又爽黄色视频| 2018国产大陆天天弄谢| 另类亚洲欧美激情| 9色porny在线观看| 欧美激情极品国产一区二区三区| 99久久综合免费| 亚洲人成77777在线视频| 99久久99久久久精品蜜桃| 日韩视频一区二区在线观看| 国产精品久久久久成人av| 51午夜福利影视在线观看| 又黄又粗又硬又大视频| 交换朋友夫妻互换小说| 99热全是精品| 后天国语完整版免费观看| 亚洲精品国产区一区二| 99久久国产精品久久久| 午夜福利视频在线观看免费| 亚洲第一欧美日韩一区二区三区 | 99re6热这里在线精品视频| 午夜免费鲁丝| 成人亚洲精品一区在线观看| 日韩免费高清中文字幕av| 久久人妻熟女aⅴ| 老司机影院成人| 亚洲国产精品一区二区三区在线| 亚洲欧美精品综合一区二区三区| 黄片大片在线免费观看| 老熟妇仑乱视频hdxx| 国产国语露脸激情在线看| 乱人伦中国视频| 极品少妇高潮喷水抽搐| 老司机午夜福利在线观看视频 | 日本黄色日本黄色录像| 一级毛片女人18水好多| 国产亚洲欧美在线一区二区| 国产在视频线精品| 欧美精品人与动牲交sv欧美| 国产精品一二三区在线看| 亚洲天堂av无毛| 久久久久久久大尺度免费视频| 99精品久久久久人妻精品| 成人国产一区最新在线观看| 亚洲成国产人片在线观看| 日本欧美视频一区| 夫妻午夜视频| 精品少妇黑人巨大在线播放| 国产1区2区3区精品| 免费高清在线观看日韩| 欧美激情久久久久久爽电影 | 久久久久久久久免费视频了| 黄色视频不卡| 亚洲欧美清纯卡通| 一级片免费观看大全| 久热爱精品视频在线9| 中文字幕另类日韩欧美亚洲嫩草| 国产免费视频播放在线视频| 肉色欧美久久久久久久蜜桃| 老司机深夜福利视频在线观看 | 在线观看舔阴道视频| 天堂中文最新版在线下载| a 毛片基地| 久久精品成人免费网站| 少妇粗大呻吟视频| 成人国产av品久久久| 啦啦啦 在线观看视频| 久久国产亚洲av麻豆专区| 精品国产一区二区三区久久久樱花| 一本综合久久免费| 久久免费观看电影| 飞空精品影院首页| 在线av久久热| 少妇被粗大的猛进出69影院| 精品国产乱子伦一区二区三区 | 精品国产国语对白av| 久久国产精品大桥未久av| 欧美人与性动交α欧美软件| 久久精品成人免费网站| 纯流量卡能插随身wifi吗| 国产97色在线日韩免费| 欧美变态另类bdsm刘玥| 窝窝影院91人妻| 亚洲人成77777在线视频| 国产精品熟女久久久久浪| 亚洲精品自拍成人| 少妇粗大呻吟视频| 国产成人免费观看mmmm| 欧美变态另类bdsm刘玥| 一级毛片女人18水好多| 中文字幕另类日韩欧美亚洲嫩草| 少妇 在线观看| 免费观看a级毛片全部| 97在线人人人人妻| 一级毛片电影观看| 99re6热这里在线精品视频| 久久久久久久久免费视频了| kizo精华| 国产成人免费观看mmmm| 免费在线观看视频国产中文字幕亚洲 | 国产日韩一区二区三区精品不卡| 丝袜喷水一区| 亚洲成国产人片在线观看| 午夜影院在线不卡| 日韩中文字幕视频在线看片| 久久国产精品男人的天堂亚洲| 99精品久久久久人妻精品| 国产亚洲一区二区精品| 波多野结衣一区麻豆| 国产无遮挡羞羞视频在线观看| 中文欧美无线码| 老熟妇仑乱视频hdxx| 女警被强在线播放| 老熟女久久久| 首页视频小说图片口味搜索| 男女午夜视频在线观看| 三级毛片av免费| 18禁观看日本| 国产极品粉嫩免费观看在线| 女人被躁到高潮嗷嗷叫费观| 一区二区三区精品91| 国产成人欧美| 高清黄色对白视频在线免费看| 国产亚洲精品一区二区www | 少妇粗大呻吟视频| 久9热在线精品视频| 另类精品久久| 热99re8久久精品国产| 亚洲专区中文字幕在线| 久久精品久久久久久噜噜老黄| av天堂在线播放| 久久久久国产一级毛片高清牌| 在线观看免费视频网站a站| av天堂久久9| 欧美日本中文国产一区发布| 日日爽夜夜爽网站| 大香蕉久久网| 日日爽夜夜爽网站| 大香蕉久久网| 国产成人a∨麻豆精品| 免费观看av网站的网址| 我要看黄色一级片免费的| 男女免费视频国产| 亚洲精品成人av观看孕妇| 波多野结衣一区麻豆| 美女视频免费永久观看网站| 99精品久久久久人妻精品| 老司机靠b影院| 免费在线观看黄色视频的| 女人高潮潮喷娇喘18禁视频| 18禁裸乳无遮挡动漫免费视频| 精品国产一区二区三区四区第35| 亚洲精品国产一区二区精华液| 一区福利在线观看| 久久久国产一区二区| kizo精华| 青青草视频在线视频观看| 热re99久久精品国产66热6| 男女边摸边吃奶| 国产精品99久久99久久久不卡| 欧美精品一区二区大全| 国产精品一区二区精品视频观看| 日韩欧美免费精品| 欧美激情久久久久久爽电影 | 久久av网站| 国产成人免费无遮挡视频| 亚洲国产日韩一区二区| 国产精品一区二区免费欧美 | 成人三级做爰电影| 搡老岳熟女国产| 中文字幕最新亚洲高清| 亚洲国产日韩一区二区| 亚洲国产欧美网| 大型av网站在线播放| 午夜两性在线视频| 欧美日韩成人在线一区二区| 国产一级毛片在线| 亚洲人成电影免费在线| 无遮挡黄片免费观看| 国产精品欧美亚洲77777| 精品少妇黑人巨大在线播放| 一级毛片精品| 91成年电影在线观看| 9191精品国产免费久久| 99国产精品一区二区蜜桃av | av片东京热男人的天堂| 巨乳人妻的诱惑在线观看| 青春草亚洲视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 久热这里只有精品99| av一本久久久久| 国产在线免费精品| 两性午夜刺激爽爽歪歪视频在线观看 | 精品福利观看| 1024视频免费在线观看| 建设人人有责人人尽责人人享有的| 香蕉国产在线看| 18禁观看日本| 日韩大片免费观看网站| 日韩精品免费视频一区二区三区| 精品一区在线观看国产| a级片在线免费高清观看视频| 青青草视频在线视频观看| 国产亚洲欧美在线一区二区| 一区二区三区乱码不卡18| 大型av网站在线播放| 欧美成人午夜精品| 亚洲一卡2卡3卡4卡5卡精品中文| 老熟女久久久| 精品欧美一区二区三区在线| 80岁老熟妇乱子伦牲交| 青青草视频在线视频观看| 日韩熟女老妇一区二区性免费视频| 女人爽到高潮嗷嗷叫在线视频| 国产成人精品久久二区二区91| 侵犯人妻中文字幕一二三四区| 亚洲美女黄色视频免费看| 悠悠久久av| av又黄又爽大尺度在线免费看| 亚洲精品粉嫩美女一区| 久久国产精品影院| 日韩大片免费观看网站| 在线天堂中文资源库| 久热爱精品视频在线9| 亚洲免费av在线视频| 欧美黑人欧美精品刺激| 9热在线视频观看99| 亚洲av男天堂| 亚洲国产日韩一区二区| 又紧又爽又黄一区二区| 久久久久视频综合| 国产一区二区三区综合在线观看| 国产精品1区2区在线观看. | 少妇被粗大的猛进出69影院| netflix在线观看网站| 亚洲五月婷婷丁香| 男女无遮挡免费网站观看| 99久久99久久久精品蜜桃| 香蕉国产在线看| 久久精品aⅴ一区二区三区四区| 亚洲欧美一区二区三区久久| 男人操女人黄网站| 午夜福利乱码中文字幕| 国产亚洲精品一区二区www | 国产成人影院久久av| 免费少妇av软件| 亚洲av片天天在线观看| 亚洲色图 男人天堂 中文字幕| 美女午夜性视频免费| netflix在线观看网站|