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

    瞬態(tài)條件下電磁軌道發(fā)射裝置絕緣體熱損傷分析

    2023-11-11 06:11:20胡鑫凱魯軍勇張嘉煒
    電工技術(shù)學(xué)報(bào) 2023年21期
    關(guān)鍵詞:模型

    胡鑫凱 魯軍勇 李 白 譚 賽 張嘉煒

    (海軍工程大學(xué)電磁能技術(shù)全國(guó)重點(diǎn)實(shí)驗(yàn)室 武漢 430033)

    0 引言

    電磁軌道發(fā)射裝置( Electromagnetic Rail Launcher, EMRL)是一種在高脈沖強(qiáng)磁場(chǎng)環(huán)境下利用洛倫茲力發(fā)射彈丸的直線加速器[1-6]。絕緣體作為EMRL 的核心部件放置在兩條導(dǎo)軌之間,起到了絕緣和支撐的雙重作用。由于電磁發(fā)射過程復(fù)雜,膛內(nèi)惡劣的熱學(xué)環(huán)境會(huì)造成絕緣體的熱損傷,導(dǎo)致絕緣體支撐強(qiáng)度下降和絕緣失效。絕緣體熱損傷主要分為兩個(gè)階段:一是電樞膛內(nèi)運(yùn)動(dòng)階段,二是電樞出膛后炮口電弧回流階段。研究這兩個(gè)階段絕緣體的熱損傷情況對(duì)提升絕緣體支撐和絕緣性能具有重要的實(shí)際意義。

    電樞膛內(nèi)運(yùn)動(dòng)階段產(chǎn)生的氣動(dòng)熱[7]對(duì)絕緣體作用時(shí)間短,熱量傳遞可忽略,主要考慮導(dǎo)軌傳遞給絕緣體的熱量,包括焦耳熱、樞-軌接觸熱及摩擦熱。李丹等[8-9]通過二維瞬態(tài)仿真發(fā)現(xiàn)絕緣體最高溫度和最大熱應(yīng)力均出現(xiàn)在與導(dǎo)軌的接觸面上,導(dǎo)軌的熱效應(yīng)對(duì)絕緣體影響不大,但連發(fā)造成的熱量累積會(huì)導(dǎo)致絕緣體強(qiáng)度降低和絕緣老化。姜遠(yuǎn)志等[10]通過多物理場(chǎng)耦合仿真,發(fā)現(xiàn)連發(fā)條件下每次發(fā)射后都會(huì)出現(xiàn)300℃的瞬態(tài)高溫,G10 絕緣體在高溫環(huán)境下反復(fù)工作,容易導(dǎo)致樹脂基體的破壞和層間結(jié)合力的減弱。但這些研究都沒有考慮電樞運(yùn)動(dòng)帶來的導(dǎo)軌熱量的時(shí)空變化對(duì)絕緣體的熱影響。

    炮口電弧回流階段絕緣體主要受到高溫電弧等離子體的熱侵蝕。電樞出膛后,發(fā)射系統(tǒng)內(nèi)大量的剩余能量會(huì)以電弧的形式通過炮口引弧裝置釋放出來[11]。由于電樞的運(yùn)動(dòng)速度遠(yuǎn)大于聲速,導(dǎo)致內(nèi)膛氣壓遠(yuǎn)低于標(biāo)準(zhǔn)大氣壓,炮口內(nèi)外的壓差是炮口電弧向膛尾回流的主要原因。J.J.Weimer 等[12]使用高速攝像機(jī)和光譜儀分析了回流電弧的組分、流動(dòng)以及溫度信息,發(fā)現(xiàn)回流溫度最高可達(dá)4 000 K,為鋁顆粒在1 atm(1 atm=1.013×105Pa)中的燃燒溫度,回流持續(xù)了60 ms,速度約600 m/s,與文獻(xiàn)[13-16]通過磁流體仿真獲取的回流溫度和速度相近,高溫電弧導(dǎo)致了絕緣體的燒蝕。P.J.Cote 等[17-18]通過激光脈沖加熱試驗(yàn),復(fù)現(xiàn)了試驗(yàn)中觀察到的絕緣體熱損傷,包括環(huán)氧樹脂熱分解以及纖維尖端的軟化。中國(guó)科學(xué)院電工研究所進(jìn)行發(fā)射試驗(yàn)后,發(fā)現(xiàn)炮口區(qū)域絕緣體燒蝕炭化最嚴(yán)重[19],表面材料因高溫而分解脫落,燒蝕增大了絕緣體表面粗糙度、提高了氧化物含量并帶來了金屬污染[20-21],導(dǎo)致絕緣體表面電阻率和閃絡(luò)電壓大幅降低[22-24]。綜上可知,目前對(duì)電弧回流引起的絕緣體熱損傷的研究主要集中在實(shí)驗(yàn)現(xiàn)象的描述和分析,缺乏對(duì)流-固-熱耦合環(huán)境下絕緣體熱侵蝕機(jī)理的研究。

    因此,本文采用三維瞬態(tài)電-磁-熱耦合模型分析了電樞膛內(nèi)運(yùn)動(dòng)過程中隨時(shí)間以及電樞位置變化的導(dǎo)軌熱量對(duì)絕緣體的影響。建立考慮G10 絕緣體細(xì)觀結(jié)構(gòu)的流-固-熱耦合模型解釋了電弧回流對(duì)絕緣體的熱侵蝕機(jī)理,并得到了試驗(yàn)驗(yàn)證。

    1 電-磁-熱耦合模型

    為研究電樞膛內(nèi)運(yùn)動(dòng)階段導(dǎo)軌對(duì)絕緣體的熱傳導(dǎo)作用,建立了考慮電樞運(yùn)動(dòng)的三維瞬態(tài)電-磁-熱耦合模型。

    1.1 控制方程

    基于準(zhǔn)靜態(tài)麥克斯韋方程和本構(gòu)方程,通過引入矢量磁位和標(biāo)量電位,得到EMRL 瞬態(tài)電磁場(chǎng)的控制方程統(tǒng)一形式為[25]

    式中,μ為相對(duì)磁導(dǎo)率;A為矢量磁位;φ為標(biāo)量電位;σ為電導(dǎo)率;t為時(shí)間;v為電樞運(yùn)動(dòng)速度。

    EMRL 的熱擴(kuò)散方程可以寫成[26]

    式中,ρ為材料的密度;cp為比定壓熱容;T為溫度;κ為導(dǎo)熱系數(shù);Q為熱源,主要由三部分組成:導(dǎo)軌焦耳熱、樞-軌接觸熱及摩擦熱,各項(xiàng)熱源分別用Qe、Qc、Qf表示,方程為

    式中,J為電流密度;ρa(bǔ)為導(dǎo)軌和電樞的平均電阻率;Hsoft為電樞硬度;P為接觸壓強(qiáng);δ為摩擦系數(shù);經(jīng)驗(yàn)接觸常數(shù)c=9.45×10-4,m=0.63[27]。

    假設(shè)電樞為剛體,其在洛倫茲力和摩擦力的共同作用下沿著導(dǎo)軌x軸運(yùn)動(dòng),運(yùn)動(dòng)方程可以寫成

    式中,M為電樞質(zhì)量,M=160 g;x為電樞的位置,起始位置為0.2 m;B為磁感應(yīng)強(qiáng)度;f為摩擦阻力;V為電樞運(yùn)動(dòng)區(qū)域。

    1.2 建立模型

    本文的研究對(duì)象為50 mm 方口徑的EMRL。身管長(zhǎng)4 m,主要由銅合金導(dǎo)軌、鋁合金電樞和G10絕緣體組成,材料物性參數(shù)見表1。在EMRL 周圍設(shè)置空氣域以求解電磁場(chǎng),鑒于EMRL 結(jié)構(gòu)的對(duì)稱特性,為提高計(jì)算效率,建立的1/4 幾何模型如圖1 所示。電樞運(yùn)動(dòng)通過改變電樞運(yùn)動(dòng)區(qū)域材料屬性的方法實(shí)現(xiàn)[28-29],即在圖1 所示的電樞運(yùn)動(dòng)區(qū)域使用鋁合金屬性代表電樞,其他區(qū)域定義為空氣屬性。在電樞和空氣的交界處,電導(dǎo)率采用平滑邏輯函數(shù)的形式,以避免不收斂。

    圖1 1/4 EMRL 幾何模型Fig.1 Geometric model of the quarter EMRL

    表1 材料物性參數(shù)Tab.1 Parameters of material properties

    1.3 定解條件

    磁場(chǎng)邊界條件:邊界S1 垂直于磁感應(yīng)強(qiáng)度方向,n×H=0;邊界S2 相切于磁感應(yīng)強(qiáng)度方向,n×A=0;空氣域邊界S3 視作無窮遠(yuǎn),n×A=0;邊界S4、S5上A的方向與電流方向一致,n×A=0。A的初始值為0。

    電場(chǎng)邊界條件:邊界S6 為導(dǎo)軌端面,給時(shí)變電壓輸入;邊界S7 為電樞運(yùn)動(dòng)區(qū)域的對(duì)稱面,設(shè)為零電位,φ=0 V;電絕緣條件適用于除S6 和S7 以外的所有邊界,n·J=0。φ的初始值為0 V。

    溫度場(chǎng)邊界條件:由于電樞運(yùn)動(dòng)時(shí)間為ms 量級(jí),所以邊界S1、S2、S3、S4 和S5 均施加絕熱邊界條件,n·(κT? )=0。T初始值為293.15 K。

    2 流-固-熱耦合模型

    電樞出膛后,由于高速摩擦,轉(zhuǎn)捩及炮口擊穿放電,炮口區(qū)域?qū)⒊霈F(xiàn)復(fù)雜的高溫空氣電弧等離子體流場(chǎng),并在氣壓差的影響下向膛尾回流。為研究電弧回流持續(xù)時(shí)間內(nèi)電弧對(duì)絕緣體增強(qiáng)體(紗線)以及基體(樹脂)的熱侵蝕作用,建立了考慮G10絕緣體細(xì)觀結(jié)構(gòu)的流-固-熱耦合三維模型。

    2.1 控制方程

    為分析絕緣體在回流電弧作用下的流固耦合與傳熱問題,引入流體狀態(tài)無量綱參數(shù)——雷諾數(shù)Re來判定回流電弧在膛內(nèi)的流動(dòng)狀態(tài),當(dāng)Re>2 320時(shí),膛內(nèi)為湍流[30]。經(jīng)計(jì)算,電弧在4 000 K 下的雷諾數(shù)為18 113.4,可知電弧的流動(dòng)狀態(tài)為湍流,故采用適用性和魯棒性高的k-ε模型[31]。

    將空氣電弧等離子體視作物性參數(shù)隨溫度變化的流體,該過程涉及流體的流動(dòng)與傳熱、流體與固體界面之間的熱對(duì)流以及固體自身的熱傳導(dǎo)。模型涉及質(zhì)量、動(dòng)量和能量三大守恒定律以及k-ε湍流方程為

    式中,fρ為流體密度;vf為流體微元的速度矢量;k和ε分別為流體的湍流動(dòng)能及耗散率;μf為流體動(dòng)力粘度;tμ為湍流粘度;kσ和εσ為湍流的普朗特常數(shù);Pk為粘性項(xiàng);Cε1和Cε2為模型常數(shù);Cμ為湍流粘度常數(shù)。

    2.2 建立模型

    G10 絕緣體是一種玻璃纖維層壓材料,用環(huán)氧樹脂粘合在一起。G10 采用玻璃纖維平面布作為增強(qiáng)相,鋪層厚度為[(0°,90°)]120,50%的纖維分布在x方向上,其余50%分布在y方向上,z方向?yàn)槠矫娌嫉亩询B方向。為準(zhǔn)確地計(jì)算電弧回流對(duì)G10 絕緣體的熱侵蝕作用,構(gòu)建了層合結(jié)構(gòu)復(fù)合材料G10 的細(xì)觀模型如圖2 所示,其中紗線寬度為0.7 mm,厚度為0.125 mm,密度為10 根/cm。細(xì)觀模型可以對(duì)紗線尺度進(jìn)行分析,相較于宏觀模型可以更好地反映紗線和樹脂的溫度分布,也不會(huì)像微觀(纖維尺度)模型那樣帶來巨量的網(wǎng)格。紗線和樹脂的物性參數(shù)以及隨溫度變化的空氣電弧等離子體物性參數(shù)[32-34]見表2 和圖3 所示。

    圖2 G10 細(xì)觀模型Fig.2 Meso-model of G10

    圖3 電弧等離子體物性參數(shù)Fig.3 Physical parameters of arc plasma

    表2 材料物性參數(shù)Tab.2 Parameters of material properties

    2.3 定解條件

    在COMSOL Multiphysics 有限元軟件中建立幾何模型并賦予物性參數(shù)后,使用了“湍流,k-ε”和“固體和流體傳熱”兩個(gè)物理場(chǎng)模塊,并通過多物理場(chǎng)中的“非等溫流動(dòng)”模塊將兩個(gè)物理場(chǎng)耦合起來。

    根據(jù)文獻(xiàn)[12-16],假設(shè)炮口回流電弧的溫度為4 000 K,流速為600 m/s,電樞出膛時(shí)內(nèi)膛氣壓接近于真空狀態(tài)?;亓鞒掷m(xù)時(shí)間與發(fā)射能級(jí)相關(guān),本文選取試驗(yàn)工況下實(shí)測(cè)的10 ms,得到定解條件如下:

    內(nèi)膛區(qū)域?yàn)榱黧w,流體邊界條件:在“湍流,k-ε”模塊中選擇邊界S8 為流體入口,對(duì)其賦予“充分發(fā)展的流動(dòng)”邊界條件,平均流速為600 m/s。邊界S9 為流體出口,設(shè)置出口壓力為0 Pa。邊界S10、S11 設(shè)置為對(duì)稱面,其余面設(shè)置為壁面。內(nèi)膛流體區(qū)域的初始速度及初始?jí)毫Χ紴?。

    G10 絕緣體區(qū)域?yàn)楣腆w,固體和流體傳熱邊界條件:邊界S8 為高溫流體入口,溫度為4 000 K;邊界S9 為出口;邊界S10、S11 設(shè)置為對(duì)稱面;由于計(jì)算時(shí)間僅為10 ms,所以其他邊界設(shè)置為絕熱條件。整個(gè)區(qū)域的初始溫度都為293.15 K。

    3 結(jié)果與討論

    3.1 電樞膛內(nèi)運(yùn)動(dòng)階段

    仿真計(jì)算得到的驅(qū)動(dòng)電流波形如圖4 所示。驅(qū)動(dòng)電流為脈沖大電流,上升沿約為1 ms,峰值約為0.442 MA,電樞出膛時(shí)刻為5.7 ms。

    圖4 驅(qū)動(dòng)電流曲線Fig.4 Driving current curve

    在電流趨膚效應(yīng)、速度趨膚效應(yīng)以及鄰近效應(yīng)的影響下,電流密度會(huì)趨于導(dǎo)軌內(nèi)表面分布,圖5為電流上升沿時(shí)刻(0.7 ms)導(dǎo)軌的電流密度分布。可以看出導(dǎo)軌內(nèi)表面的電流密度明顯高于導(dǎo)軌內(nèi)部,且更趨向于集中在導(dǎo)軌兩側(cè)邊緣,即與絕緣體的接觸面。再加上樞-軌接觸電阻以及滑動(dòng)摩擦的影響,集中于導(dǎo)軌內(nèi)表面的熱量會(huì)通過與絕緣體的接觸面?zhèn)鬟f到絕緣體上。

    圖5 0.7 ms 時(shí)刻導(dǎo)軌的電流密度分布Fig.5 Current density distribution of the rail at 0.7 ms

    1.2 ms(最高溫度時(shí)刻)和5.7 ms(電樞出膛時(shí)刻)導(dǎo)軌和絕緣體的溫度分布分別如圖6 和圖7 所示。由圖6 和圖7 可知隨著電樞的運(yùn)動(dòng),導(dǎo)軌和絕緣體的高溫區(qū)域不斷向口部方向延伸,但尾部至口部的溫升呈逐漸減小的趨勢(shì),最高溫度始終出現(xiàn)在電樞起始位置附近。原因如下:首先,電樞在起始位置附近速度較低,加熱時(shí)間較長(zhǎng);其次,該段時(shí)間為電流急劇上升段,渦流效應(yīng)更加明顯,對(duì)應(yīng)的焦耳熱效應(yīng)更為顯著;最后,該階段電樞與導(dǎo)軌的接觸狀態(tài)較差,接觸電阻很大。

    圖6 1.2 ms 時(shí)刻導(dǎo)軌和絕緣體的溫度分布Fig.6 Temperature distributions of the rail and insulator at 1.2 ms

    圖7 5.7 ms 時(shí)刻導(dǎo)軌和絕緣體的溫度分布Fig.7 Temperature distributions of the rail and insulator at 5.7 ms

    由于G10 絕緣體的絕緣等級(jí)為B 級(jí),這意味著G10 絕緣體最高允許工作溫度為403 K[10]。單次發(fā)射過程中,在1.2 ms 時(shí)刻與導(dǎo)軌接觸的絕緣體內(nèi)表面出現(xiàn)最高溫度401 K,接近最高允許工作溫度,如圖6 中放大的A 區(qū)所示??紤]連續(xù)發(fā)射過程中絕緣體溫度的累積,需對(duì)導(dǎo)軌進(jìn)行冷卻系統(tǒng)設(shè)計(jì)[35],尤其是導(dǎo)軌尾部、電樞起始位置附近以及導(dǎo)軌內(nèi)表面的有效冷卻,以避免絕緣材料老化而影響支撐性能。

    3.2 炮口電弧回流階段

    電樞出膛后,炮口區(qū)域的高溫電弧等離子體會(huì)在短時(shí)間內(nèi)流回內(nèi)膛,炮口電弧回流示意圖如圖8所示。圖9 和圖10 分別為電弧持續(xù)回流10 ms 后G10 絕緣體紗線和樹脂的溫度分布。其中,玻璃纖維紗線的熱分解溫度約1 000 K,環(huán)氧樹脂的熱分解溫度約450 K。

    圖8 炮口電弧回流示意圖Fig.8 Schematic of muzzle arc blowback

    圖9 紗線的溫度分布Fig.9 Temperature distribution of yarn

    圖10 樹脂的溫度分布Fig.10 Temperature distribution of resin

    從圖9 可以看出,得益于纖維紗線較高的分解溫度,僅距離絕緣體口部1.3 mm 范圍內(nèi)的紗線達(dá)到熱分解點(diǎn)。熱侵蝕深度從口部0 mm 處的0.125 mm逐漸降至口部1.3 mm 處的0 mm,如圖中放大的B區(qū)所示??芍亓麟娀H會(huì)對(duì)絕緣體口部1.3 mm 范圍內(nèi)的紗線造成“倒角”狀熱損傷,即回流電弧對(duì)絕緣體紗線的影響很小。

    從圖10 可以看出,由于環(huán)氧樹脂分解溫度較低,距離絕緣體口部10 mm 范圍內(nèi)的樹脂均達(dá)到熱分解點(diǎn)。熱侵蝕深度沿尾部方向略呈下降趨勢(shì),口部0 mm 處最深為0.125 mm,如圖中放大的C 區(qū)所示,平均侵蝕深度約0.07 mm?;谟?jì)算效率的考慮,本文僅建立了10 mm 長(zhǎng)模型,但通過對(duì)仿真結(jié)果的合理推測(cè),可知回流電弧對(duì)絕緣體內(nèi)表面電弧流經(jīng)區(qū)域的樹脂均會(huì)造成一定程度的熱分解,樹脂的大范圍熱分解是電弧回流帶來的主要熱損傷。

    3.3 試驗(yàn)結(jié)果及分析

    在實(shí)際的發(fā)射過程中,質(zhì)量為160 g 的電樞在電磁力的作用下向炮口運(yùn)動(dòng)?;贓MRL 動(dòng)態(tài)特性測(cè)試平臺(tái)[36-37],對(duì)電樞通過B 點(diǎn)探頭位置的時(shí)間值進(jìn)行插值擬合,獲得發(fā)射過程中電樞的實(shí)時(shí)速度及位置。電樞實(shí)測(cè)和仿真的運(yùn)動(dòng)速度及位移曲線如圖11 所示,兩者一致性較好。

    圖11 電樞運(yùn)動(dòng)曲線Fig.11 Motion curve of the armature

    對(duì)試驗(yàn)后的EMRL 樣機(jī)進(jìn)行拆解,觀察膛內(nèi)G10 絕緣體的熱損傷情況。圖12 為試驗(yàn)后距離炮口500 mm 范圍內(nèi)的絕緣體的典型宏觀狀態(tài),其中放大區(qū)域?yàn)榻^緣體口部局部位置。區(qū)域D 出現(xiàn)了“倒角”狀損傷,判斷為該處的纖維紗線和樹脂均被熱分解所致。膛內(nèi)電弧流經(jīng)區(qū)域的絕緣體內(nèi)表面纖維紗線之間出現(xiàn)了大量的孔洞,并附著炭化分解產(chǎn)物以及銀色金屬污染物,原因是填充在紗線之間的樹脂被大范圍熱分解炭化,同時(shí)回流電弧中的熔融鋁顆粒會(huì)附著在絕緣體表面??芍?,回流電弧僅會(huì)對(duì)絕緣體口部很小范圍紗線造成“倒角”狀熱損傷,樹脂的熱分解損傷占主導(dǎo)地位,試驗(yàn)現(xiàn)象與仿真得到的熱侵蝕規(guī)律一致。

    圖12 試驗(yàn)后絕緣體的典型宏觀狀態(tài)Fig.12 Typical macroscopic state of the insulator after test

    炮口電弧回流帶來的絕緣體內(nèi)表面樹脂熱分解炭化和金屬污染,會(huì)導(dǎo)致絕緣體表面電阻率和閃絡(luò)電壓大幅降低。抑制炮口電弧回流以及提高絕緣材料的抗熱侵蝕能力有助于避免絕緣失效。采用炮口消弧裝置是一種抑制炮口電弧回流的有效手段;提高纖維紗線密度以減小單位體積內(nèi)易被熱分解的樹脂含量,或者選用耐熱性能優(yōu)異的改性樹脂都可以很好地提高絕緣材料的抗熱侵蝕能力。

    4 結(jié)論

    本文建立了三維瞬態(tài)電-磁-熱耦合模型以及考慮G10 絕緣體細(xì)觀結(jié)構(gòu)的流-固-熱耦合模型來分析絕緣體的熱損傷情況,得到以下結(jié)論:

    1)電樞膛內(nèi)運(yùn)動(dòng)階段,絕緣體的高溫區(qū)域隨電樞運(yùn)動(dòng)而向口部延伸,但尾部至口部的溫升逐漸降低,最高溫度始終出現(xiàn)在電樞起始位置附近,接近絕緣體最高允許工作溫度。

    2)炮口電弧回流階段,回流電弧對(duì)絕緣體紗線的影響很小,僅口部1.3 mm 范圍內(nèi)的紗線形成“倒角”狀熱損傷。但對(duì)絕緣體內(nèi)表面電弧流經(jīng)區(qū)域的樹脂造成了大范圍熱分解,樹脂的熱分解損傷占主導(dǎo)地位。試驗(yàn)現(xiàn)象與仿真得到的熱侵蝕規(guī)律一致。

    3)連續(xù)發(fā)射會(huì)帶來絕緣體溫度的累積,對(duì)導(dǎo)軌進(jìn)行冷卻系統(tǒng)設(shè)計(jì)以避免絕緣材料老化而影響支撐性能。電弧回流會(huì)帶來絕緣體內(nèi)表面炭化和金屬污染,采用炮口消弧裝置以及提高樹脂的抗熱侵蝕能力有助于避免絕緣失效。

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    性高湖久久久久久久久免费观看| 欧美精品一区二区免费开放| 成人三级做爰电影| 中文字幕精品免费在线观看视频| 少妇粗大呻吟视频| 美女高潮到喷水免费观看| 一级,二级,三级黄色视频| 国产97色在线日韩免费| 超碰成人久久| 午夜福利视频在线观看免费| 久久精品国产亚洲av高清一级| 国产欧美日韩综合在线一区二区| 亚洲精品av麻豆狂野| 老司机在亚洲福利影院| 在线永久观看黄色视频| 一级,二级,三级黄色视频| 亚洲avbb在线观看| 国产精品99久久99久久久不卡| 丁香欧美五月| 18禁美女被吸乳视频| 亚洲人成伊人成综合网2020| 国产成+人综合+亚洲专区| 好男人电影高清在线观看| 十八禁网站网址无遮挡| 操出白浆在线播放| 在线观看免费高清a一片| 欧美 日韩 精品 国产| 亚洲第一青青草原| tocl精华| 99精品在免费线老司机午夜| 丝瓜视频免费看黄片| 国产午夜精品久久久久久| 欧美成狂野欧美在线观看| 成人18禁高潮啪啪吃奶动态图| 777久久人妻少妇嫩草av网站| 男女午夜视频在线观看| 狠狠婷婷综合久久久久久88av| 亚洲国产av影院在线观看| 真人做人爱边吃奶动态| 国产男靠女视频免费网站| 性高湖久久久久久久久免费观看| 亚洲成a人片在线一区二区| 久久国产亚洲av麻豆专区| 不卡一级毛片| 国产精品.久久久| 午夜精品久久久久久毛片777| 精品一区二区三区av网在线观看 | 欧美精品高潮呻吟av久久| 久久国产精品大桥未久av| bbb黄色大片| 99精品在免费线老司机午夜| 亚洲精品在线观看二区| 男女之事视频高清在线观看| 国产亚洲午夜精品一区二区久久| 亚洲成人免费电影在线观看| 欧美日本中文国产一区发布| 久9热在线精品视频| 青青草视频在线视频观看| 亚洲中文字幕日韩| 国产亚洲欧美在线一区二区| 免费在线观看黄色视频的| 国产免费视频播放在线视频| 国产精品影院久久| 久久精品国产亚洲av高清一级| 婷婷丁香在线五月| 亚洲专区国产一区二区| 国产精品二区激情视频| 久久久精品国产亚洲av高清涩受| 亚洲欧美一区二区三区久久| 久久亚洲精品不卡| 电影成人av| 亚洲欧洲日产国产| 黄网站色视频无遮挡免费观看| 美女高潮喷水抽搐中文字幕| 99久久国产精品久久久| 老司机深夜福利视频在线观看| 日韩成人在线观看一区二区三区| 1024视频免费在线观看| av国产精品久久久久影院| 欧美日韩亚洲高清精品| 欧美日韩亚洲综合一区二区三区_| 操美女的视频在线观看| 黄色视频在线播放观看不卡| 午夜福利,免费看| av超薄肉色丝袜交足视频| 最新的欧美精品一区二区| h视频一区二区三区| av在线播放免费不卡| 国产男女超爽视频在线观看| 777米奇影视久久| 下体分泌物呈黄色| 亚洲成av片中文字幕在线观看| 免费人妻精品一区二区三区视频| 久久国产亚洲av麻豆专区| 国产在线一区二区三区精| 久久国产精品大桥未久av| 成人国语在线视频| 国产精品98久久久久久宅男小说| 99re6热这里在线精品视频| 一区在线观看完整版| 少妇猛男粗大的猛烈进出视频| 嫩草影视91久久| 香蕉丝袜av| 高清av免费在线| 亚洲五月色婷婷综合| 久热这里只有精品99| 1024视频免费在线观看| 国产99久久九九免费精品| 伊人久久大香线蕉亚洲五| 国产一区二区三区在线臀色熟女 | 青草久久国产| 国产成人精品久久二区二区免费| 久久精品熟女亚洲av麻豆精品| 亚洲av国产av综合av卡| 国产1区2区3区精品| 亚洲伊人色综图| 亚洲精品中文字幕一二三四区 | 国产欧美亚洲国产| 老司机午夜福利在线观看视频 | 久久中文字幕人妻熟女| 美女扒开内裤让男人捅视频| 亚洲国产欧美日韩在线播放| 久久精品国产亚洲av高清一级| 色尼玛亚洲综合影院| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜| 女人精品久久久久毛片| 亚洲欧美激情在线| 国产成人精品在线电影| 午夜成年电影在线免费观看| 日韩免费av在线播放| 亚洲九九香蕉| 中文字幕最新亚洲高清| 亚洲少妇的诱惑av| 1024视频免费在线观看| 男男h啪啪无遮挡| 丁香六月天网| 久久性视频一级片| 老司机深夜福利视频在线观看| 80岁老熟妇乱子伦牲交| 热99re8久久精品国产| 亚洲va日本ⅴa欧美va伊人久久| 精品人妻在线不人妻| 国产av又大| 热re99久久国产66热| 日韩成人在线观看一区二区三区| 汤姆久久久久久久影院中文字幕| 亚洲午夜理论影院| 久久久欧美国产精品| 别揉我奶头~嗯~啊~动态视频| 精品国产国语对白av| 国产精品香港三级国产av潘金莲| 国产欧美日韩一区二区精品| 手机成人av网站| 深夜精品福利| 狂野欧美激情性xxxx| 少妇 在线观看| 亚洲精品国产精品久久久不卡| 精品久久久久久电影网| 欧美成人午夜精品| 免费观看av网站的网址| 国产成人av激情在线播放| 亚洲伊人色综图| 欧美精品人与动牲交sv欧美| 99精品在免费线老司机午夜| 两个人免费观看高清视频| 国产欧美日韩一区二区精品| 欧美日韩精品网址| 国产免费现黄频在线看| 伦理电影免费视频| 新久久久久国产一级毛片| 天天操日日干夜夜撸| 人人妻人人添人人爽欧美一区卜| 涩涩av久久男人的天堂| 欧美亚洲 丝袜 人妻 在线| 91精品国产国语对白视频| 啦啦啦免费观看视频1| 中文欧美无线码| 国产亚洲精品一区二区www | 一区在线观看完整版| 无限看片的www在线观看| 另类亚洲欧美激情| 夜夜夜夜夜久久久久| 国产精品98久久久久久宅男小说| 91成年电影在线观看| 熟女少妇亚洲综合色aaa.| 精品国产超薄肉色丝袜足j| 久久精品人人爽人人爽视色| 成人影院久久| 国产男靠女视频免费网站| 中国美女看黄片| 国产免费视频播放在线视频| 亚洲成人国产一区在线观看| 777米奇影视久久| 蜜桃国产av成人99| 亚洲av第一区精品v没综合| 91大片在线观看| av天堂在线播放| 国产精品 国内视频| av国产精品久久久久影院| 涩涩av久久男人的天堂| 最黄视频免费看| 亚洲va日本ⅴa欧美va伊人久久| 一本一本久久a久久精品综合妖精| 免费不卡黄色视频| 老司机靠b影院| 久久99一区二区三区| 精品一区二区三区av网在线观看 | 精品国产一区二区久久| 激情视频va一区二区三区| 久久中文字幕人妻熟女| 日日爽夜夜爽网站| 黄片小视频在线播放| 精品久久久久久电影网| 老司机靠b影院| 99国产精品一区二区蜜桃av | 中文字幕制服av| 最近最新中文字幕大全免费视频| 国产成人精品久久二区二区免费| av视频免费观看在线观看| 热99re8久久精品国产| 久久国产精品影院| 高清黄色对白视频在线免费看| 亚洲欧美日韩高清在线视频 | 久久人人爽av亚洲精品天堂| 一区福利在线观看| 亚洲精品乱久久久久久| 国产高清视频在线播放一区| 免费人妻精品一区二区三区视频| 老司机在亚洲福利影院| 久久人妻熟女aⅴ| 高潮久久久久久久久久久不卡| 亚洲专区中文字幕在线| 国产精品二区激情视频| 精品福利观看| 老司机靠b影院| 99精品欧美一区二区三区四区| 桃花免费在线播放| 成人av一区二区三区在线看| 免费人妻精品一区二区三区视频| 欧美精品亚洲一区二区| 国产亚洲精品久久久久5区| 欧美精品av麻豆av| 天天影视国产精品| 国产97色在线日韩免费| www.自偷自拍.com| 91国产中文字幕| 色老头精品视频在线观看| 视频区欧美日本亚洲| 制服人妻中文乱码| avwww免费| aaaaa片日本免费| 精品国产亚洲在线| 这个男人来自地球电影免费观看| 在线十欧美十亚洲十日本专区| 一二三四社区在线视频社区8| 午夜久久久在线观看| 99国产极品粉嫩在线观看| 青草久久国产| 日韩一区二区三区影片| 国产精品 国内视频| 午夜日韩欧美国产| 19禁男女啪啪无遮挡网站| 成人18禁在线播放| 巨乳人妻的诱惑在线观看| 国产国语露脸激情在线看| 在线观看66精品国产| 少妇裸体淫交视频免费看高清 | 国产精品美女特级片免费视频播放器 | 欧美精品高潮呻吟av久久| 99精品久久久久人妻精品| 97人妻天天添夜夜摸| 90打野战视频偷拍视频| 免费看十八禁软件| 免费观看人在逋| 免费不卡黄色视频| 精品福利观看| 欧美黑人欧美精品刺激| 久久久久久亚洲精品国产蜜桃av| 国产精品98久久久久久宅男小说| 国产伦理片在线播放av一区| 制服诱惑二区| 9191精品国产免费久久| 他把我摸到了高潮在线观看 | 9191精品国产免费久久| 天天躁狠狠躁夜夜躁狠狠躁| 18在线观看网站| 99九九在线精品视频| 色在线成人网| 日韩人妻精品一区2区三区| av天堂久久9| 中文字幕制服av| 少妇裸体淫交视频免费看高清 | 久久精品亚洲av国产电影网| av免费在线观看网站| 大陆偷拍与自拍| 丁香六月欧美| 国产三级黄色录像| 黑人欧美特级aaaaaa片| 丁香欧美五月| 亚洲第一欧美日韩一区二区三区 | 久久香蕉激情| 欧美 日韩 精品 国产| 嫁个100分男人电影在线观看| www日本在线高清视频| 国产成人精品在线电影| 国产成人精品久久二区二区免费| videosex国产| 色播在线永久视频| 91成年电影在线观看| 99精品在免费线老司机午夜| 19禁男女啪啪无遮挡网站| 国产精品一区二区免费欧美| 欧美黑人欧美精品刺激| 我的亚洲天堂| 一本一本久久a久久精品综合妖精| 丰满人妻熟妇乱又伦精品不卡| 啦啦啦免费观看视频1| 最新在线观看一区二区三区| 丁香六月天网| 菩萨蛮人人尽说江南好唐韦庄| 三上悠亚av全集在线观看| 国产精品偷伦视频观看了| 侵犯人妻中文字幕一二三四区| 亚洲av欧美aⅴ国产| 在线十欧美十亚洲十日本专区| 久久中文字幕人妻熟女| 又黄又粗又硬又大视频| 亚洲精品中文字幕一二三四区 | 十八禁人妻一区二区| 成年版毛片免费区| 人成视频在线观看免费观看| 激情在线观看视频在线高清 | 亚洲精品国产一区二区精华液| 国产成人系列免费观看| 老司机午夜福利在线观看视频 | 99热国产这里只有精品6| 欧美日韩黄片免| 久久性视频一级片| 老汉色∧v一级毛片| 大香蕉久久成人网| 999精品在线视频| 中文欧美无线码| 建设人人有责人人尽责人人享有的| 电影成人av| 老司机在亚洲福利影院| 久久精品91无色码中文字幕| 在线观看66精品国产| 一级片免费观看大全| 国产老妇伦熟女老妇高清| 狠狠狠狠99中文字幕| 中文字幕精品免费在线观看视频| 狠狠狠狠99中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 久久人人97超碰香蕉20202| 美女主播在线视频| 亚洲av成人不卡在线观看播放网| 日韩一卡2卡3卡4卡2021年| 大型黄色视频在线免费观看| 国产一区二区在线观看av| 黑人巨大精品欧美一区二区蜜桃| 操美女的视频在线观看| 亚洲欧美色中文字幕在线| 日韩中文字幕视频在线看片| av在线播放免费不卡| 久久久国产一区二区| 美女午夜性视频免费| av网站在线播放免费| av在线播放免费不卡| 国产激情久久老熟女| 亚洲精品久久成人aⅴ小说| 国产精品久久电影中文字幕 | 岛国毛片在线播放| 在线观看免费高清a一片| 久久久国产精品麻豆| 中文字幕av电影在线播放| 国产三级黄色录像| 在线观看舔阴道视频| 免费高清在线观看日韩| 香蕉久久夜色| 大型av网站在线播放| 免费日韩欧美在线观看| 亚洲精华国产精华精| 麻豆成人av在线观看| 一级,二级,三级黄色视频| 一级黄色大片毛片| 久久精品亚洲av国产电影网| 中文字幕制服av| 国产成人av教育| 12—13女人毛片做爰片一| 肉色欧美久久久久久久蜜桃| 日韩有码中文字幕| 99久久精品国产亚洲精品| 一二三四在线观看免费中文在| 精品一区二区三区视频在线观看免费 | a级毛片黄视频| 精品福利永久在线观看| 亚洲欧美激情在线| av福利片在线| 国产野战对白在线观看| 熟女少妇亚洲综合色aaa.| 久久久精品免费免费高清| 成人影院久久| 在线天堂中文资源库| 精品高清国产在线一区| 一夜夜www| 欧美精品啪啪一区二区三区| 视频在线观看一区二区三区| 亚洲精品中文字幕在线视频| 日本wwww免费看| 久久久国产欧美日韩av| 法律面前人人平等表现在哪些方面| 一区福利在线观看| 国产在线精品亚洲第一网站| 性高湖久久久久久久久免费观看| 我要看黄色一级片免费的| 99国产精品一区二区蜜桃av | a级片在线免费高清观看视频| 久久久久久久久免费视频了| 久久久国产精品麻豆| av电影中文网址| 亚洲一码二码三码区别大吗| 一级黄色大片毛片| 别揉我奶头~嗯~啊~动态视频| av线在线观看网站| 亚洲精品乱久久久久久| 欧美变态另类bdsm刘玥| 精品少妇久久久久久888优播| 日韩熟女老妇一区二区性免费视频| 丰满少妇做爰视频| 欧美精品人与动牲交sv欧美| 妹子高潮喷水视频| 日韩一卡2卡3卡4卡2021年| 亚洲国产毛片av蜜桃av| 国产高清videossex| 99国产精品免费福利视频| 少妇猛男粗大的猛烈进出视频| 亚洲九九香蕉| 久热爱精品视频在线9| 欧美成狂野欧美在线观看| 999久久久国产精品视频| 国产区一区二久久| 一区二区日韩欧美中文字幕| 亚洲欧美日韩高清在线视频 | 欧美日韩精品网址| 国产黄频视频在线观看| 日韩制服丝袜自拍偷拍| 午夜福利在线观看吧| 手机成人av网站| 在线观看免费日韩欧美大片| 韩国精品一区二区三区| 亚洲成国产人片在线观看| 日韩精品免费视频一区二区三区| 国产人伦9x9x在线观看| 美女高潮到喷水免费观看| 无限看片的www在线观看| 老鸭窝网址在线观看| 日韩中文字幕欧美一区二区| 国产成人av教育| 国产精品久久久av美女十八| 午夜成年电影在线免费观看| 国产一区二区三区综合在线观看| 天天躁日日躁夜夜躁夜夜| 狠狠狠狠99中文字幕| 亚洲欧美日韩另类电影网站| 国产欧美日韩综合在线一区二区| 久久国产精品男人的天堂亚洲| 国产高清国产精品国产三级| 成人亚洲精品一区在线观看| 激情在线观看视频在线高清 | 啦啦啦中文免费视频观看日本| 99国产综合亚洲精品| 手机成人av网站| 香蕉国产在线看| 国产精品影院久久| 国产精品免费大片| 成年人免费黄色播放视频| 国产在视频线精品| 免费看十八禁软件| 成人国产一区最新在线观看| 免费久久久久久久精品成人欧美视频| 老汉色av国产亚洲站长工具| 久久精品亚洲精品国产色婷小说| 最近最新中文字幕大全免费视频| 999久久久精品免费观看国产| 亚洲伊人久久精品综合| avwww免费| 精品少妇久久久久久888优播| 叶爱在线成人免费视频播放| 黑人巨大精品欧美一区二区mp4| 欧美日韩一级在线毛片| 天天躁狠狠躁夜夜躁狠狠躁| 成人18禁高潮啪啪吃奶动态图| 色94色欧美一区二区| 曰老女人黄片| 精品熟女少妇八av免费久了| 欧美精品高潮呻吟av久久| 91大片在线观看| 亚洲人成77777在线视频| 精品福利观看| 国产欧美日韩综合在线一区二区| av在线播放免费不卡| 高潮久久久久久久久久久不卡| 丝袜美腿诱惑在线| 天堂俺去俺来也www色官网| 视频区图区小说| 国产真人三级小视频在线观看| 久久热在线av| 黄色视频,在线免费观看| 老司机在亚洲福利影院| 国产极品粉嫩免费观看在线| h视频一区二区三区| 亚洲 欧美一区二区三区| 51午夜福利影视在线观看| 国产日韩欧美视频二区| 两人在一起打扑克的视频| 久久天堂一区二区三区四区| 大片免费播放器 马上看| 母亲3免费完整高清在线观看| 纯流量卡能插随身wifi吗| 91麻豆av在线| 成年人免费黄色播放视频| 色婷婷久久久亚洲欧美| 久久久国产成人免费| 免费少妇av软件| xxxhd国产人妻xxx| 久久99一区二区三区| 飞空精品影院首页| 大码成人一级视频| 99在线人妻在线中文字幕 | 精品国产一区二区三区四区第35| 看免费av毛片| 在线看a的网站| 女人被躁到高潮嗷嗷叫费观| 国产男女内射视频| 亚洲av日韩精品久久久久久密| 黄色视频在线播放观看不卡| 三级毛片av免费| 一区福利在线观看| 久久影院123| 高清视频免费观看一区二区| 无限看片的www在线观看| 这个男人来自地球电影免费观看| 国产日韩欧美亚洲二区| 午夜久久久在线观看| 男女床上黄色一级片免费看| 亚洲中文av在线| 老熟妇乱子伦视频在线观看| 国产成人av教育| 国产精品国产av在线观看| 亚洲中文日韩欧美视频| 又大又爽又粗| 免费看a级黄色片| 国产精品一区二区免费欧美| 亚洲精品中文字幕一二三四区 | 日韩人妻精品一区2区三区| 成人特级黄色片久久久久久久 | 精品国产超薄肉色丝袜足j| 一区在线观看完整版| 亚洲伊人色综图| 久久国产精品人妻蜜桃| 黄频高清免费视频| 国产高清视频在线播放一区| 国产精品 欧美亚洲| 国产男女超爽视频在线观看| 丰满饥渴人妻一区二区三| 可以免费在线观看a视频的电影网站| 精品人妻熟女毛片av久久网站| 国产精品久久久久久精品电影小说| 国产国语露脸激情在线看| 国产精品国产av在线观看| 精品久久久精品久久久| 后天国语完整版免费观看| 一个人免费看片子| 国精品久久久久久国模美| 中亚洲国语对白在线视频| 精品亚洲成a人片在线观看| 我的亚洲天堂| av网站免费在线观看视频| 亚洲精品在线观看二区| 青青草视频在线视频观看| 大型av网站在线播放| 亚洲av成人不卡在线观看播放网| 成人18禁高潮啪啪吃奶动态图| 精品免费久久久久久久清纯 | 五月开心婷婷网| 国产视频一区二区在线看| 少妇的丰满在线观看| 国产精品九九99| 国产亚洲av高清不卡| 精品久久蜜臀av无| 性高湖久久久久久久久免费观看| 午夜两性在线视频| 亚洲综合色网址| 亚洲av电影在线进入| 妹子高潮喷水视频| 欧美中文综合在线视频| 国产午夜精品久久久久久| 亚洲五月色婷婷综合| 成年女人毛片免费观看观看9 | 亚洲av欧美aⅴ国产| 老熟妇仑乱视频hdxx| 国产97色在线日韩免费| 精品一品国产午夜福利视频| 丰满迷人的少妇在线观看| 啦啦啦 在线观看视频| 国产成人影院久久av| 久久久久久人人人人人| 性少妇av在线| 在线观看免费视频网站a站| 成人18禁高潮啪啪吃奶动态图| 亚洲中文av在线| 国产老妇伦熟女老妇高清| 老司机午夜十八禁免费视频| 十八禁网站免费在线|