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

    邊光滑有限元-邊界元耦合法計(jì)算二維瞬態(tài)渦流場(chǎng)

    2021-03-17 03:24:24王洋洋蔣興良
    關(guān)鍵詞:計(jì)算精度鋁板有限元法

    王洋洋,蔣興良

    (輸配電裝備及系統(tǒng)安全與新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(重慶大學(xué)),重慶 400044)

    在瞬態(tài)渦流場(chǎng)研究中,常用的數(shù)值計(jì)算方法有積分法[1-2]、有限元法[3]、邊界元法以及有限元邊界元耦合法[4-7]. 雖然傳統(tǒng)有限元法(FEM)[8]作為最成功的數(shù)值方法之一,已被廣泛用于解決工程和科學(xué)的各種實(shí)際問題,但是因?yàn)槠溆泄逃械娜秉c(diǎn),如剖分網(wǎng)格密度以及計(jì)算階數(shù)的復(fù)雜程度都對(duì)有限元計(jì)算精度影響非常大. 單元數(shù)量的增加,及其插值函數(shù)精度的提高,才能得到更精確的數(shù)值解. 特別是在某些復(fù)雜的設(shè)計(jì)優(yōu)化實(shí)例中,有限元法需要重復(fù)多次進(jìn)行網(wǎng)格剖分以及重要部分網(wǎng)格加密操作,這不僅增加了設(shè)計(jì)的難度,而且還降低了設(shè)計(jì)計(jì)算的效率. 因此國(guó)內(nèi)外學(xué)者一直致力于研究如何克服傳統(tǒng)有限元法的缺點(diǎn).

    近年來(lái),在努力開發(fā)先進(jìn)數(shù)值方法時(shí),文獻(xiàn)[9]將無(wú)網(wǎng)格中的應(yīng)變光滑技術(shù)與傳統(tǒng)有限元法相結(jié)合,引入梯度光滑的操作概念,提出一系列廣義梯度光滑有限元法. 其包括用梯度光滑技術(shù)分別對(duì)構(gòu)造的節(jié)點(diǎn)光滑域、邊光滑域以及面光滑域進(jìn)行光滑操作而形成的點(diǎn)光滑有限元法(NS-FEM)[10]、邊光滑有限元法(ES-FEM)[11]以及面光滑有限元法(FS-FEM)[12]. 3種方法各有優(yōu)缺點(diǎn). 相對(duì)于傳統(tǒng)有限元法,NS-FEM在精度和效率方面有所提高. 但是NS-FEM在空間上是穩(wěn)定的在時(shí)間上可能是不穩(wěn)定的,并且由于“過軟”特性而不能直接應(yīng)用于解決動(dòng)態(tài)問題[13]. 然而,ES-FEM可以克服時(shí)間不穩(wěn)定性[14],由于本文研究的是時(shí)變渦流場(chǎng),因此克服時(shí)間不穩(wěn)定性對(duì)于時(shí)變渦流場(chǎng)分析尤其重要. 此外,對(duì)于固體力學(xué)問題,具有線性積分的ES-FEM比線性有限元法更加準(zhǔn)確和有效[15]. 至此,ES-FEM近年來(lái)得到了廣泛的應(yīng)用. 對(duì)于FS-FEM更適用于三維問題.

    本文提出一種用ES-FEM與邊界元耦合法計(jì)算二維瞬態(tài)渦流場(chǎng). 該耦合法同時(shí)具有ES-FEM計(jì)算精度高以及邊界元法占用計(jì)算機(jī)內(nèi)存少等優(yōu)點(diǎn). 首先將問題域離散化為一組三角形單元,然后進(jìn)一步構(gòu)造與三角形單元邊相關(guān)的積分域,采用伽遼金加權(quán)余量法,推導(dǎo)出基于邊光滑有限元-邊界元耦合算法的離散公式,在導(dǎo)體區(qū)域內(nèi)采用ES-FEM,非導(dǎo)體區(qū)域內(nèi)采用邊界元法. 研究導(dǎo)體區(qū)域的磁場(chǎng)分布,并與測(cè)量結(jié)果和有限元-邊界元耦合法計(jì)算結(jié)果進(jìn)行對(duì)比. 計(jì)算結(jié)果表明,在相同網(wǎng)格密度條件下,基于ES-FEM邊界元耦合法計(jì)算精度更高,該方法具有更好的應(yīng)用前景.

    1 電磁場(chǎng)的分析模型

    電磁分析的問題實(shí)際上是在給定邊界條件下求解一組麥克斯韋方程組的問題. 如圖1所示,為開域渦流場(chǎng)求解區(qū)域示意圖. 整個(gè)電磁場(chǎng)求解域分為渦流區(qū)V1和非渦流區(qū)V2,Γ1是V1和V2的內(nèi)部邊界,非渦流區(qū)V2包括無(wú)源電流的非渦流區(qū)和有源電流的非渦流區(qū).V1區(qū)域內(nèi)采用ES-FEM,V2區(qū)域內(nèi)采用BEM, 在內(nèi)部邊界上滿足的邊界條件為

    (1)

    (2)

    ν

    (3)

    圖1 開域渦流場(chǎng)求解示意圖

    2 邊光滑有限元-邊界元耦合法介紹

    2.1 邊光滑有限元法

    邊光滑有限元法是在求解有限元的基礎(chǔ)上對(duì)其求解域進(jìn)一步光滑,形成多個(gè)光滑子單元,然后在新形成的光滑子單元內(nèi)引入梯度光滑操作,梯度光滑操作是指在新形成的光滑域上對(duì)有限元求解的系數(shù)矩陣進(jìn)行光滑平均運(yùn)算,其數(shù)學(xué)表達(dá)式為

    其中Ax為光滑域的面積. 圖2為邊光滑有限元法的計(jì)算流程圖.

    圖2 計(jì)算流程圖

    對(duì)于二維求解域,先用劃分軟件Hypermesh對(duì)其進(jìn)行三角形劃分,其中生成N個(gè)節(jié)點(diǎn)和p條邊. 傳統(tǒng)有限元的離散矩陣[16]為

    (4)

    (5)

    式中AL為光滑域的面積,在二維求解域單元內(nèi),AL的表達(dá)式可根據(jù)參考文獻(xiàn)[9]計(jì)算得到.

    圖3 基于邊的光滑域

    對(duì)式(7)進(jìn)行格林公式換算,那么對(duì)光滑域內(nèi)位函數(shù)的面積分可以轉(zhuǎn)化為光滑域邊界上的線積分,即

    (6)

    將式(6)代入式(5)中,可得

    (7)

    其中:

    (8)

    (9)

    對(duì)光滑域ΩL采用高斯積分,式(12)的光滑域邊界積分形式為

    (10)

    式中:NA為電磁場(chǎng)求解域邊界L的個(gè)數(shù),NB為高斯點(diǎn)的個(gè)數(shù),δq為高斯點(diǎn)積分的權(quán)函數(shù).

    由式(10)可知,dij是由光滑域各個(gè)邊中點(diǎn)處的形函數(shù)求得,而各個(gè)中點(diǎn)處的形函數(shù)可通過場(chǎng)節(jié)點(diǎn)插值得到,如圖4所示,點(diǎn)1、2、3、4分別為場(chǎng)節(jié)點(diǎn),A、B分別為中心節(jié)點(diǎn),g1、g2、g3、g4分別為高斯積分點(diǎn)[17].

    基于邊L的光滑有限元系統(tǒng)的局部光滑域的系數(shù)矩陣可以寫成:

    (11)

    對(duì)上述局部系數(shù)矩陣進(jìn)行整體組裝,得到域內(nèi)整體系數(shù)矩陣為

    (12)

    圖4 ES-FEM形函數(shù)分布

    2.2 邊界元法

    對(duì)于非渦流區(qū)應(yīng)用格林公式,可得對(duì)應(yīng)邊界積分方程[5]為

    (13)

    采用伽遼金余量法對(duì)式(13)離散,得

    (14)

    2.3 基于邊光滑有限元邊界元耦合算法

    將式(4)和式(14)聯(lián)立,可得

    (15)

    對(duì)式(15)中的導(dǎo)數(shù)項(xiàng),采用向后差分,即

    (16)

    (17)

    (18)

    3 計(jì)算算例與分析

    3.1 脈沖線圈-鋁板表面磁場(chǎng)分布

    脈沖線圈-鋁板的結(jié)構(gòu)參數(shù):矩形截面為0.6 mm×4.8 mm的銅絲纏繞成內(nèi)徑、外徑分別為3.2、25.5 mm,匝數(shù)為30的圓柱形線圈[19]. 脈沖線圈流過的瞬時(shí)電流波形,如圖5所示. 在線圈正上方放置0.8 mm厚的鋁板,電導(dǎo)率δ=3.48×107s/m. 線圈與鋁板的間隙為2 mm. 采用本文所提的基于邊有限元邊界元耦合法求得鋁板表面的磁感應(yīng)強(qiáng)度,其結(jié)果分別與測(cè)量值以及有限元-邊界元法的計(jì)算結(jié)果進(jìn)行對(duì)比,如圖6~8所示.

    圖5 脈沖電流隨時(shí)間變化

    圖6 磁矢量位A線的分布云圖

    圖7 切向磁感應(yīng)強(qiáng)度隨時(shí)間的變化(r=10 mm)

    由圖6可知,離線圈越近,磁矢量位A值越大,與實(shí)際情況一致,說(shuō)明本文的計(jì)算方法合理. 為了驗(yàn)證本文所提方法的計(jì)算精度,設(shè)相對(duì)誤差Bi為

    (19)

    式中:B0為測(cè)量值,B1是計(jì)算值,i=x、y分別代表磁場(chǎng)沿切向、垂直方向的磁場(chǎng)分量.

    圖8 垂直磁感應(yīng)強(qiáng)度隨時(shí)間的變化(r=7.6 mm)

    由圖7、8可知,基于邊光滑有限元邊界元耦合法與測(cè)量結(jié)果一致,因此可用本文所提的方法計(jì)算渦流場(chǎng). 在網(wǎng)格單元?jiǎng)澐置芏认嗤臈l件下,運(yùn)行兩種方法效率對(duì)比見表1,雖然本文所提方法的運(yùn)行時(shí)間稍長(zhǎng)于有限元邊界元耦合法,但其計(jì)算精度遠(yuǎn)遠(yuǎn)高于有限元邊界元耦合法,另外,對(duì)于目前計(jì)算機(jī)的計(jì)算能力而言,本文通過改變模型的大小,對(duì)比兩種方法的運(yùn)行時(shí)間. 由表2可知,兩種方法運(yùn)行時(shí)間的差別可以忽略. 由此可知,在計(jì)算渦流場(chǎng)時(shí)基于邊光滑有限元邊界元法能大大提高計(jì)算精度.

    表1 效率對(duì)比

    表2 運(yùn)行時(shí)間對(duì)比

    綜上所述,在進(jìn)行渦流場(chǎng)分析時(shí),應(yīng)優(yōu)先考慮基于邊光滑有限元邊界元耦合法.

    3.2 鋁板表面磁場(chǎng)分布的影響因素

    3.2.1 電流密度變化的影響

    脈沖電流-鋁板模型參數(shù)如上所述,在進(jìn)行電流密度對(duì)其磁場(chǎng)影響分析時(shí),保持模型的其他參數(shù)不變,改變流過脈沖線圈的電流密度Js從10×106A/m2到180×106A/m2變化,間隔為30×106A/m2. 則上述電流密度變化的情況下鋁板表面磁場(chǎng)分量分布如圖9所示.

    圖9 磁場(chǎng)強(qiáng)度隨渦流密度的變化

    如圖9所示,磁場(chǎng)強(qiáng)度隨渦流密度的增大而增大,這是因?yàn)闇u流密度Js與系數(shù)矩陣F成正相關(guān),即當(dāng)Js增大時(shí),F(xiàn)增大. 由式(20)可知,隨著F增大,所求磁矢量位A增大,因此磁場(chǎng)強(qiáng)度增大.

    3.2.2 脈沖線圈-鋁板間隙變化的影響

    脈沖電流-鋁板模型參數(shù)如上所述,在進(jìn)行脈沖線圈-鋁板間隙對(duì)其磁場(chǎng)影響分析時(shí),保持模型的其他參數(shù)不變,改變間隙h從2 mm到8 mm之間變化,間隔為1 mm. 則上述間隙變化的情況下鋁板表面磁場(chǎng)分量分布如圖10所示.

    圖10 磁場(chǎng)強(qiáng)度隨間隙的變化

    由圖10可知,隨著脈沖線圈與鋁板之間間隙h的增大,磁場(chǎng)強(qiáng)度逐漸減小. 這說(shuō)明如果為了增大鋁板中的渦流作用,即產(chǎn)生大的脈沖力,那么脈沖線圈需盡可能近的接近鋁板.

    4 結(jié) 論

    1)提出了一種計(jì)算渦流場(chǎng)分布的基于邊光滑有限元邊界元耦合法,該方法即具有ES-FEM計(jì)算精度高和邊界元法占用內(nèi)存少的優(yōu)點(diǎn),并詳細(xì)給出了該方法的推導(dǎo)過程. 同時(shí),利用本文所提的基于邊光滑有限元邊界元耦合法計(jì)算分析脈沖線圈上方鋁板表面的磁場(chǎng)分布,與測(cè)量結(jié)果對(duì)比表明,該方法是正確合理的.

    2)在單元網(wǎng)格密度劃分相同的條件下,將本文得到的計(jì)算結(jié)果與傳統(tǒng)有限元邊界元耦合法計(jì)算得到的結(jié)果對(duì)比,雖然本文所提方法的運(yùn)行時(shí)間稍長(zhǎng)于有限元邊界元耦合法,但其計(jì)算精度遠(yuǎn)遠(yuǎn)高于有限元邊界元耦合法,另外,對(duì)于目前計(jì)算機(jī)的計(jì)算能力而言,兩種方法運(yùn)行時(shí)間的差別可以忽略.

    3)在計(jì)算模型的其他參數(shù)保持不變,只改變渦流密度的情況下,磁場(chǎng)強(qiáng)度隨渦流密度的增大而增大且曲線增大率幾乎保持不變,同理,磁感應(yīng)強(qiáng)度隨著脈沖線圈與鋁板之間間隙h的增大而減小,且其曲線降低率也幾乎保持不變.

    猜你喜歡
    計(jì)算精度鋁板有限元法
    哪塊邊角料的面積大
    大型鋁板拉伸機(jī)液壓底座的設(shè)計(jì)計(jì)算
    正交各向異性材料裂紋疲勞擴(kuò)展的擴(kuò)展有限元法研究
    雙曲弧形純鋁板內(nèi)幕墻的施工探討
    智能城市(2018年7期)2018-07-10 08:30:24
    基于SHIPFLOW軟件的某集裝箱船的阻力計(jì)算分析
    廣東造船(2018年1期)2018-03-19 15:50:50
    長(zhǎng)脈沖激光與連續(xù)激光對(duì)鋁板熱破壞仿真對(duì)比
    單元類型和尺寸對(duì)拱壩壩體應(yīng)力和計(jì)算精度的影響
    三維有限元法在口腔正畸生物力學(xué)研究中發(fā)揮的作用
    鋼箱計(jì)算失效應(yīng)變的沖擊試驗(yàn)
    集成對(duì)稱模糊數(shù)及有限元法的切削力預(yù)測(cè)
    久久精品国产鲁丝片午夜精品| 三级经典国产精品| 亚洲真实伦在线观看| 国产蜜桃级精品一区二区三区| 一级毛片电影观看 | 91麻豆精品激情在线观看国产| 18禁在线播放成人免费| 三级毛片av免费| av黄色大香蕉| av.在线天堂| 2021天堂中文幕一二区在线观| 精品久久久久久久久av| 亚洲av免费高清在线观看| 久久久久久久久久久丰满| 日韩国内少妇激情av| 一进一出抽搐gif免费好疼| av免费在线看不卡| 亚洲精品久久国产高清桃花| 久久久久精品国产欧美久久久| 亚洲欧美成人精品一区二区| 国产三级在线视频| 久久久久久大精品| 亚洲成人久久爱视频| 精品一区二区三区av网在线观看| 亚洲成人av在线免费| 3wmmmm亚洲av在线观看| 国产黄色小视频在线观看| 成人av在线播放网站| 国产在线精品亚洲第一网站| 人妻夜夜爽99麻豆av| 色吧在线观看| 一本一本综合久久| 欧美潮喷喷水| 日韩精品青青久久久久久| 国产国拍精品亚洲av在线观看| 最近手机中文字幕大全| 国产 一区 欧美 日韩| 波野结衣二区三区在线| 精品乱码久久久久久99久播| 亚洲av熟女| 干丝袜人妻中文字幕| 欧美高清成人免费视频www| 国内少妇人妻偷人精品xxx网站| av免费在线看不卡| 99在线视频只有这里精品首页| 十八禁国产超污无遮挡网站| 午夜福利高清视频| 色综合站精品国产| 国产日本99.免费观看| 亚洲中文日韩欧美视频| 亚洲美女黄片视频| 日日干狠狠操夜夜爽| 国产精品久久久久久久久免| 网址你懂的国产日韩在线| 国产不卡一卡二| 91在线精品国自产拍蜜月| 桃色一区二区三区在线观看| 久久久久国产网址| 精品国产三级普通话版| 日韩高清综合在线| 十八禁国产超污无遮挡网站| 天天躁夜夜躁狠狠久久av| 亚洲av成人av| 51国产日韩欧美| 成人欧美大片| 亚洲激情五月婷婷啪啪| 国产精品爽爽va在线观看网站| 狂野欧美白嫩少妇大欣赏| 亚洲中文日韩欧美视频| 亚洲精品成人久久久久久| 国产黄色小视频在线观看| 蜜桃亚洲精品一区二区三区| 国产免费男女视频| 久久这里只有精品中国| 午夜福利视频1000在线观看| 人妻丰满熟妇av一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 免费电影在线观看免费观看| 在线免费十八禁| 成人一区二区视频在线观看| 嫩草影视91久久| 一级毛片电影观看 | 噜噜噜噜噜久久久久久91| 国产中年淑女户外野战色| 午夜免费男女啪啪视频观看 | 成年av动漫网址| 熟妇人妻久久中文字幕3abv| 看非洲黑人一级黄片| 69人妻影院| 极品教师在线视频| 国产精品亚洲一级av第二区| 国产一区二区激情短视频| 国产精品美女特级片免费视频播放器| 波多野结衣巨乳人妻| 日韩av不卡免费在线播放| 亚洲无线在线观看| 亚洲欧美精品综合久久99| 欧美激情久久久久久爽电影| 99久久中文字幕三级久久日本| 日韩高清综合在线| 少妇的逼水好多| av中文乱码字幕在线| 深爱激情五月婷婷| 亚洲不卡免费看| .国产精品久久| 菩萨蛮人人尽说江南好唐韦庄 | 国产黄色视频一区二区在线观看 | 久久久午夜欧美精品| 亚洲国产精品成人久久小说 | 免费看光身美女| 亚洲国产日韩欧美精品在线观看| av福利片在线观看| 国产不卡一卡二| 亚洲av熟女| 久久久国产成人精品二区| 可以在线观看毛片的网站| 亚洲电影在线观看av| 99久久精品热视频| 99久国产av精品国产电影| 最近最新中文字幕大全电影3| 不卡视频在线观看欧美| 少妇被粗大猛烈的视频| 日本a在线网址| 国产精品伦人一区二区| 日韩国内少妇激情av| 国产真实伦视频高清在线观看| 国产一区二区激情短视频| 久99久视频精品免费| 亚洲av中文字字幕乱码综合| 婷婷精品国产亚洲av在线| 国产一区二区亚洲精品在线观看| 麻豆精品久久久久久蜜桃| 中文字幕久久专区| 国产男靠女视频免费网站| 激情 狠狠 欧美| 日本撒尿小便嘘嘘汇集6| 日日摸夜夜添夜夜添av毛片| 看片在线看免费视频| 欧美一区二区亚洲| 国产一区二区亚洲精品在线观看| 国产在视频线在精品| 久久6这里有精品| 亚洲国产精品sss在线观看| 亚洲av五月六月丁香网| 日本三级黄在线观看| 亚洲av免费高清在线观看| 露出奶头的视频| 欧美日本亚洲视频在线播放| 99热网站在线观看| 简卡轻食公司| 一个人看视频在线观看www免费| 18禁黄网站禁片免费观看直播| 69av精品久久久久久| 高清毛片免费看| 你懂的网址亚洲精品在线观看 | 两个人的视频大全免费| 听说在线观看完整版免费高清| 成人鲁丝片一二三区免费| 不卡视频在线观看欧美| 亚洲18禁久久av| 乱系列少妇在线播放| 中国美白少妇内射xxxbb| 成人三级黄色视频| 春色校园在线视频观看| 日韩在线高清观看一区二区三区| 亚洲欧美成人综合另类久久久 | 日日撸夜夜添| 99热6这里只有精品| 国产伦在线观看视频一区| 精品久久久久久久久亚洲| 一个人观看的视频www高清免费观看| 又爽又黄a免费视频| 日本-黄色视频高清免费观看| АⅤ资源中文在线天堂| 精品午夜福利视频在线观看一区| 亚洲第一区二区三区不卡| 国产一区二区激情短视频| 美女高潮的动态| 99久久九九国产精品国产免费| 亚洲一级一片aⅴ在线观看| 美女内射精品一级片tv| 观看美女的网站| 18+在线观看网站| 日韩欧美精品v在线| 麻豆一二三区av精品| 99热这里只有是精品50| 国产精品永久免费网站| 久久久久性生活片| 一本久久中文字幕| 一区二区三区免费毛片| 一个人免费在线观看电影| 激情 狠狠 欧美| avwww免费| 国产黄a三级三级三级人| 在线免费观看不下载黄p国产| 熟妇人妻久久中文字幕3abv| 日韩三级伦理在线观看| 最新在线观看一区二区三区| 日本黄色视频三级网站网址| 嫩草影院新地址| 国产综合懂色| 亚洲自拍偷在线| 亚洲婷婷狠狠爱综合网| 在线播放无遮挡| 一区福利在线观看| 国产黄色视频一区二区在线观看 | 午夜精品国产一区二区电影 | 国产精品久久久久久亚洲av鲁大| 午夜福利在线观看吧| 国产探花在线观看一区二区| 国产精品av视频在线免费观看| 男人的好看免费观看在线视频| 久久久欧美国产精品| 国产一区二区在线观看日韩| 国产美女午夜福利| 免费av不卡在线播放| 波多野结衣高清作品| 久久99热6这里只有精品| 国产黄片美女视频| 国产一区二区三区在线臀色熟女| 国产 一区 欧美 日韩| 日韩人妻高清精品专区| 精品久久久久久久末码| 亚洲性久久影院| 丰满的人妻完整版| 非洲黑人性xxxx精品又粗又长| 婷婷亚洲欧美| 国产精品爽爽va在线观看网站| 午夜精品国产一区二区电影 | 美女免费视频网站| 国产真实乱freesex| 久久久久国产网址| 久久久精品94久久精品| 精品人妻一区二区三区麻豆 | av女优亚洲男人天堂| 国产免费男女视频| 有码 亚洲区| 亚洲国产精品久久男人天堂| 波野结衣二区三区在线| 精品一区二区免费观看| 日韩制服骚丝袜av| 亚洲欧美日韩无卡精品| 午夜免费男女啪啪视频观看 | 久久久久久久午夜电影| 欧美xxxx黑人xx丫x性爽| 中文字幕精品亚洲无线码一区| 级片在线观看| av在线老鸭窝| 午夜精品一区二区三区免费看| 亚洲精华国产精华液的使用体验 | 精品久久久久久久久久久久久| 欧美国产日韩亚洲一区| 成人高潮视频无遮挡免费网站| 国产精品亚洲一级av第二区| 国产午夜精品久久久久久一区二区三区 | 国产伦精品一区二区三区四那| 99热这里只有是精品50| 免费观看人在逋| 精品国产三级普通话版| 神马国产精品三级电影在线观看| 成人毛片a级毛片在线播放| 九九久久精品国产亚洲av麻豆| 欧美日韩精品成人综合77777| 一a级毛片在线观看| 99久久精品热视频| 日韩制服骚丝袜av| 国产精品人妻久久久影院| 十八禁网站免费在线| 国产精品久久视频播放| 欧美激情国产日韩精品一区| 久久久久久国产a免费观看| 免费黄网站久久成人精品| 国产爱豆传媒在线观看| 美女cb高潮喷水在线观看| 欧美人与善性xxx| 女同久久另类99精品国产91| 国产成人91sexporn| 麻豆一二三区av精品| 亚洲,欧美,日韩| 身体一侧抽搐| 亚洲自拍偷在线| 1000部很黄的大片| 在线观看66精品国产| 亚洲成av人片在线播放无| 国产久久久一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 一本一本综合久久| 性插视频无遮挡在线免费观看| 男女视频在线观看网站免费| 久久久久久大精品| 激情 狠狠 欧美| 热99在线观看视频| 长腿黑丝高跟| 午夜a级毛片| 成人国产麻豆网| 久久精品久久久久久噜噜老黄 | 美女大奶头视频| 国产精品一区www在线观看| 国产黄片美女视频| 亚洲欧美日韩卡通动漫| 一个人看视频在线观看www免费| 国产精品三级大全| 日本一本二区三区精品| 高清日韩中文字幕在线| 在线免费十八禁| 亚洲av二区三区四区| 黄色视频,在线免费观看| 日韩大尺度精品在线看网址| 精品久久久久久久久亚洲| 99热全是精品| 91久久精品电影网| 日本欧美国产在线视频| 国产视频一区二区在线看| 亚洲av熟女| 中国美白少妇内射xxxbb| 九色成人免费人妻av| 老女人水多毛片| 国产一区二区三区av在线 | 淫妇啪啪啪对白视频| 色噜噜av男人的天堂激情| 国产精品一二三区在线看| ponron亚洲| 亚洲av免费高清在线观看| 在线a可以看的网站| 亚洲欧美日韩卡通动漫| 日韩一区二区视频免费看| 国产精品一及| 国产精品永久免费网站| 久久精品国产亚洲av香蕉五月| 久久久国产成人精品二区| 日韩精品青青久久久久久| 99riav亚洲国产免费| 看免费成人av毛片| 天美传媒精品一区二区| 亚洲性久久影院| 淫秽高清视频在线观看| 国产精品人妻久久久久久| 五月玫瑰六月丁香| 级片在线观看| 国产精品久久久久久精品电影| 欧美成人精品欧美一级黄| 麻豆久久精品国产亚洲av| 久久久久久久午夜电影| 精品久久久噜噜| avwww免费| 国产欧美日韩一区二区精品| 老熟妇乱子伦视频在线观看| 亚洲av熟女| 欧美绝顶高潮抽搐喷水| 在线观看美女被高潮喷水网站| 国产激情偷乱视频一区二区| 国产一区二区激情短视频| 亚洲丝袜综合中文字幕| 久久精品综合一区二区三区| 真人做人爱边吃奶动态| 九九热线精品视视频播放| 丝袜喷水一区| 国产一区亚洲一区在线观看| 亚洲在线观看片| 国产精品免费一区二区三区在线| 久久人妻av系列| 一进一出好大好爽视频| 国产又黄又爽又无遮挡在线| 日韩欧美三级三区| 欧美最黄视频在线播放免费| 国产国拍精品亚洲av在线观看| 大又大粗又爽又黄少妇毛片口| 国产成年人精品一区二区| 老师上课跳d突然被开到最大视频| 草草在线视频免费看| 亚洲第一电影网av| 色吧在线观看| 亚洲在线观看片| 国产欧美日韩精品亚洲av| 99热精品在线国产| 亚洲最大成人av| 亚洲中文日韩欧美视频| 中文字幕熟女人妻在线| 国产午夜精品论理片| 听说在线观看完整版免费高清| 欧美日本视频| 观看美女的网站| 亚洲欧美日韩东京热| 久久6这里有精品| 看十八女毛片水多多多| 一级毛片我不卡| 精品不卡国产一区二区三区| 99久久精品国产国产毛片| 亚洲欧美精品自产自拍| 在现免费观看毛片| 老熟妇仑乱视频hdxx| 简卡轻食公司| 最好的美女福利视频网| 女的被弄到高潮叫床怎么办| 热99在线观看视频| 国产乱人偷精品视频| 自拍偷自拍亚洲精品老妇| 春色校园在线视频观看| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品国产清高在天天线| 国产中年淑女户外野战色| 中国美白少妇内射xxxbb| 国产成年人精品一区二区| 日本黄色片子视频| 亚洲国产高清在线一区二区三| 国产精品爽爽va在线观看网站| 精品人妻一区二区三区麻豆 | 91在线精品国自产拍蜜月| 国产精品电影一区二区三区| 国产不卡一卡二| 美女 人体艺术 gogo| 黄色配什么色好看| 亚洲av成人av| 在线看三级毛片| 亚洲成人久久性| 成人av在线播放网站| 午夜福利高清视频| 非洲黑人性xxxx精品又粗又长| 免费看a级黄色片| 99久国产av精品| 久久精品久久久久久噜噜老黄 | 如何舔出高潮| 国产精品久久久久久精品电影| 国产精品一区二区免费欧美| 有码 亚洲区| 国产亚洲精品av在线| 久久综合国产亚洲精品| 大香蕉久久网| 给我免费播放毛片高清在线观看| 男女啪啪激烈高潮av片| 国产v大片淫在线免费观看| 国产高清激情床上av| 国产大屁股一区二区在线视频| 日韩三级伦理在线观看| 久久久久久久久中文| 免费观看人在逋| 看十八女毛片水多多多| 一本一本综合久久| 亚洲真实伦在线观看| 亚洲精品成人久久久久久| 一级av片app| 久久精品久久久久久噜噜老黄 | 色av中文字幕| 一进一出好大好爽视频| 国产激情偷乱视频一区二区| 男女之事视频高清在线观看| 亚洲经典国产精华液单| 亚洲最大成人手机在线| 国产av在哪里看| 久久草成人影院| 精品一区二区免费观看| 可以在线观看毛片的网站| 美女xxoo啪啪120秒动态图| 日韩,欧美,国产一区二区三区 | 嫩草影院新地址| 亚洲专区国产一区二区| 人妻少妇偷人精品九色| 深夜a级毛片| 老司机午夜福利在线观看视频| 久久精品人妻少妇| av在线播放精品| 久久久精品大字幕| 噜噜噜噜噜久久久久久91| 熟妇人妻久久中文字幕3abv| 久久久久国内视频| 我的老师免费观看完整版| 久久精品国产亚洲网站| 97人妻精品一区二区三区麻豆| 亚洲美女搞黄在线观看 | 欧美高清成人免费视频www| av视频在线观看入口| 黄色日韩在线| 麻豆国产av国片精品| 久久欧美精品欧美久久欧美| 亚洲av不卡在线观看| 亚洲一级一片aⅴ在线观看| 麻豆一二三区av精品| 嫩草影院新地址| 国产精品亚洲一级av第二区| 午夜激情福利司机影院| 一区福利在线观看| 欧美最黄视频在线播放免费| 超碰av人人做人人爽久久| 国内精品一区二区在线观看| 男人和女人高潮做爰伦理| 嫩草影视91久久| 蜜桃亚洲精品一区二区三区| 亚州av有码| 午夜视频国产福利| av在线播放精品| 亚洲第一区二区三区不卡| 久久精品国产亚洲网站| 国产蜜桃级精品一区二区三区| 免费人成在线观看视频色| 欧美性感艳星| 亚洲人成网站在线播放欧美日韩| 99热6这里只有精品| 国内揄拍国产精品人妻在线| 啦啦啦韩国在线观看视频| 在线观看美女被高潮喷水网站| 日本成人三级电影网站| 一级黄片播放器| 色噜噜av男人的天堂激情| 久久久精品大字幕| 观看美女的网站| 全区人妻精品视频| 99热这里只有是精品在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 美女 人体艺术 gogo| 男人舔女人下体高潮全视频| 一级毛片aaaaaa免费看小| 亚洲av一区综合| 国产熟女欧美一区二区| 秋霞在线观看毛片| 色综合亚洲欧美另类图片| 少妇的逼水好多| 一个人观看的视频www高清免费观看| av在线蜜桃| 老司机午夜福利在线观看视频| 亚洲av美国av| 久久久久免费精品人妻一区二区| 成人无遮挡网站| 日韩一本色道免费dvd| 亚洲欧美精品综合久久99| 国产精品免费一区二区三区在线| av在线亚洲专区| 高清毛片免费观看视频网站| 久久久午夜欧美精品| 成年av动漫网址| 天堂√8在线中文| 欧美激情国产日韩精品一区| 99久国产av精品国产电影| 深夜a级毛片| 亚洲av成人精品一区久久| 观看美女的网站| 日本-黄色视频高清免费观看| 最新在线观看一区二区三区| 亚州av有码| 精品少妇黑人巨大在线播放 | 亚洲国产精品成人综合色| 寂寞人妻少妇视频99o| 欧美3d第一页| 久久久久久久久久黄片| 日韩国内少妇激情av| 麻豆国产av国片精品| 国产一区二区在线观看日韩| 99久久精品一区二区三区| 一区二区三区四区激情视频 | 床上黄色一级片| 伦理电影大哥的女人| 久久久国产成人免费| 精华霜和精华液先用哪个| 成人美女网站在线观看视频| 狠狠狠狠99中文字幕| 成人毛片a级毛片在线播放| 亚洲av不卡在线观看| 女生性感内裤真人,穿戴方法视频| 白带黄色成豆腐渣| 精品无人区乱码1区二区| 国产伦精品一区二区三区视频9| 人人妻人人看人人澡| 18禁在线无遮挡免费观看视频 | 国产久久久一区二区三区| 综合色av麻豆| 亚洲经典国产精华液单| 久久久久性生活片| 免费不卡的大黄色大毛片视频在线观看 | 国产蜜桃级精品一区二区三区| 18+在线观看网站| 高清毛片免费观看视频网站| 亚洲精品日韩av片在线观看| 中文字幕久久专区| 毛片一级片免费看久久久久| 久久欧美精品欧美久久欧美| 久久人妻av系列| 69人妻影院| 22中文网久久字幕| 别揉我奶头 嗯啊视频| 国产精品爽爽va在线观看网站| 你懂的网址亚洲精品在线观看 | 在线免费观看不下载黄p国产| 久久精品国产清高在天天线| 欧美绝顶高潮抽搐喷水| 国产精品一区www在线观看| 成人国产麻豆网| 亚洲成av人片在线播放无| 国国产精品蜜臀av免费| 日日摸夜夜添夜夜添av毛片| 狠狠狠狠99中文字幕| 亚洲人成网站在线播放欧美日韩| 最近在线观看免费完整版| 国产爱豆传媒在线观看| 免费看av在线观看网站| 男女下面进入的视频免费午夜| 别揉我奶头 嗯啊视频| 国产精品人妻久久久影院| 亚洲丝袜综合中文字幕| 深爱激情五月婷婷| 一个人看的www免费观看视频| 亚洲va在线va天堂va国产| av在线蜜桃| 欧美色欧美亚洲另类二区| 成人三级黄色视频| 欧美高清性xxxxhd video| 99热只有精品国产| 99久久精品一区二区三区| 欧美性感艳星| 亚洲av五月六月丁香网| 俺也久久电影网| 亚洲成人久久性| 久久久久性生活片| 一a级毛片在线观看| 国产成人a区在线观看| 一个人免费在线观看电影| 国产乱人视频|