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

    J-TEXT托卡馬克上電磁彈丸注入系統(tǒng)的X型電樞設(shè)計(jì)

    2022-10-14 06:33:42陳忠勇張維康唐俊輝夏勝國(guó)
    電工技術(shù)學(xué)報(bào) 2022年19期

    陳忠勇 張維康 唐俊輝 李 峰 夏勝國(guó)

    J-TEXT托卡馬克上電磁彈丸注入系統(tǒng)的X型電樞設(shè)計(jì)

    陳忠勇1,2張維康1,2唐俊輝3李 峰1,2夏勝國(guó)3

    (1. 磁約束聚變與等離子體國(guó)際合作聯(lián)合實(shí)驗(yàn)室(華中科技大學(xué)電氣與電子工程學(xué)院) 武漢 430074 2. 強(qiáng)電磁工程與新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(華中科技大學(xué)電氣與電子工程學(xué)院) 武漢 430074 3. 脈沖功率技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室(華中科技大學(xué)) 武漢 430074)

    國(guó)際熱核聚變實(shí)驗(yàn)堆(ITER)計(jì)劃是我國(guó)參與規(guī)模最大的國(guó)際科技合作項(xiàng)目,目標(biāo)是驗(yàn)證大型托卡馬克裝置實(shí)現(xiàn)聚變能的可行性。等離子體大破裂是ITER安全運(yùn)行面臨的最大威脅,會(huì)對(duì)裝置造成嚴(yán)重的損壞,因此緩解破裂危害是ITER亟待解決的關(guān)鍵問題。目前緩解破裂危害的基本策略是主動(dòng)注入大量粒子,但現(xiàn)有的破裂緩解系統(tǒng)尚不能完全滿足ITER的需求。該文通過分析現(xiàn)有破裂緩解系統(tǒng)的特點(diǎn),介紹J-TEXT托卡馬克上新一代用于破裂緩解的電磁彈丸注入系統(tǒng)(EMI),該系統(tǒng)使用電磁力發(fā)射彈丸,可有效提高彈丸的注入速度并縮短響應(yīng)時(shí)間,克服了其他系統(tǒng)的局限性。電樞是EMI的核心部件,在發(fā)射中起決定性作用,通過介紹尾翼接觸型電樞的結(jié)構(gòu)及功能,說明這類電樞運(yùn)用于EMI時(shí)在減速階段存在樞軌電接觸性能不足、運(yùn)動(dòng)穩(wěn)定性不夠高等局限性,由此根據(jù)EMI減速階段對(duì)電樞性能產(chǎn)生的特殊需求,設(shè)計(jì)了一款具有X型結(jié)構(gòu)的新型固體電樞。仿真結(jié)果表明,電樞的電磁、機(jī)械性能滿足發(fā)射要求,電接觸特性良好。在發(fā)射性能試驗(yàn)中,電樞完成了加速至520m/s后主動(dòng)減速至0m/s;電樞-彈丸分離性能試驗(yàn)中實(shí)現(xiàn)了電樞、彈丸的穩(wěn)定分離,彈丸的飛行速度為358m/s,且該速度能隨著加速能量的增加進(jìn)一步提高。EMI為托卡馬克等離子體破裂緩解貢獻(xiàn)了一種高效注入雜質(zhì)的先進(jìn)方案,該文提出的X型電樞具有良好的發(fā)射性能,為EMI提供了一種優(yōu)良的新型固體電樞結(jié)構(gòu)及設(shè)計(jì)方法。

    托卡馬克 破裂緩解系統(tǒng) 電磁發(fā)射 軌道炮 固體電樞

    0 引言

    聚變能是被認(rèn)為有可能從根本上解決能源問題的方式之一,它的研發(fā)是我國(guó)的戰(zhàn)略性重大科技問題,被美國(guó)工程院評(píng)為21世紀(jì)十四大科技挑戰(zhàn)之一。托卡馬克是一種使用磁約束等離子體實(shí)現(xiàn)核聚變的裝置,但由于其自身特性或工程技術(shù)因素,在放電過程中可能出現(xiàn)等離子體約束變差并導(dǎo)致放電終止的情況,即“大破裂”(以下簡(jiǎn)稱為“破裂”)[1-2]。等離子體破裂的危害主要包括熱負(fù)荷、電磁負(fù)荷及逃逸電流三個(gè)方面[3-5]。

    通過對(duì)等離子體的參數(shù)進(jìn)行監(jiān)控能夠降低發(fā)生破裂的風(fēng)險(xiǎn),但并不能完全避免破裂的發(fā)生,對(duì)于國(guó)際熱核聚變實(shí)驗(yàn)堆(International Thermonuclear Experimental Reactor, ITER)而言,為了實(shí)現(xiàn)聚變反應(yīng)需在高參數(shù)下運(yùn)行[6],此時(shí)發(fā)生破裂的可能性更大,一旦發(fā)生破裂而不加以緩解,裝置將會(huì)被嚴(yán)重地毀壞。因此一旦確定等離子體將發(fā)生破裂,必須通過破裂緩解系統(tǒng)主動(dòng)將等離子體快速關(guān)斷從而保護(hù)裝置[7]?,F(xiàn)階段ITER破裂緩解系統(tǒng)的基本策略是主動(dòng)注入大量的粒子,產(chǎn)生足夠高的輻射功率以耗散等離子體的能量,由此發(fā)展的系統(tǒng)包括大量雜質(zhì)氣體注入(Massive Gas?Injection, MGI)、散裂彈丸注入(Shattered?cryogenic?Pellet?Injection, SPI)[8]、小型低溫彈丸注入(small?cryogenic?Pellet?Injection, PI)[9]等。

    目前,ITER破裂緩解系統(tǒng)的設(shè)計(jì)研究主要集中在MGI和SPI。MGI的工作原理較為簡(jiǎn)單,儲(chǔ)存在工作腔中的高壓氣體在快速閥門打開后,經(jīng)過氣路注入真空室內(nèi)[10]。SPI則將氣體降溫冷凍成彈丸,再用壓縮氣體將彈丸發(fā)射出去,經(jīng)過彎管破碎后形成碎片進(jìn)入真空室[11]。

    MGI由于本身較簡(jiǎn)單的結(jié)構(gòu),可以在數(shù)ms內(nèi)開啟閥門并釋放較多的雜質(zhì)氣體,但雜質(zhì)粒子主要沉積在等離子體邊緣附近,與等離子體的混合效率并不理想,只能達(dá)到約20%[12]。與MGI相比,SPI雖然能實(shí)現(xiàn)更大的穿透深度,但是對(duì)于ITER尺度的裝置,SPI可能難以充分穿透等離子體[13]。同時(shí),受推進(jìn)氣體聲速的限制,MGI與SPI存在飛行速度上限,其中氖氣的飛行速度約為580m/s[14];使用氦氣作為推進(jìn)氣體的SPI彈丸最高飛行速度約為400m/s[15]。若將現(xiàn)有破裂緩解系統(tǒng)用于ITER,MGI、SPI的響應(yīng)時(shí)間將超過約10ms的最佳破裂預(yù)警時(shí)間[6]。

    因此,為了克服現(xiàn)有破裂緩解系統(tǒng)的局限性,在發(fā)展MGI與SPI的同時(shí),ITER也在鼓勵(lì)發(fā)展新型破裂緩解系統(tǒng)。R. Raman等提出了電磁粒子注入器(ElectromagneticParticleInjector, EPI)的概念:使用軌道炮發(fā)射攜帶彈丸的C型電樞,電樞處于兩導(dǎo)軌間,當(dāng)導(dǎo)軌通以脈沖電流時(shí),電流及產(chǎn)生的磁場(chǎng)形成洛倫茲力驅(qū)動(dòng)電樞[16]。在平臺(tái)的實(shí)驗(yàn)里,EPI證明了將軌道炮與電磁發(fā)射運(yùn)用于托卡馬克的可行性,能在峰值24kA的脈沖電流作用下將4.1g的C型電樞及0.2g彈丸加速至約200m/s,通過撞擊阻攔靶板減速電樞并釋放彈丸,分離后的彈丸以200m/s速度沿原軌跡飛行,而電樞隨后撞擊前方的引導(dǎo)靶板并發(fā)生偏轉(zhuǎn),最終落入下方沙箱中完成回收;在無攜帶彈丸實(shí)驗(yàn)中將1.9g輕質(zhì)C型電樞加速至600m/s[17]。

    EPI使用的C型電樞通過撞擊靶板獲得沖擊力釋放彈丸并發(fā)生偏轉(zhuǎn)[17],但碰撞位置的偏移會(huì)令彈丸獲得橫向速度分量,使飛行軌跡發(fā)散并降低可重復(fù)性,且這一過程損壞了電樞結(jié)構(gòu)并產(chǎn)生較多碎片,碎片將跟隨彈丸進(jìn)入托卡馬克,影響裝置正常運(yùn)行。為了確保在ITER熱猝滅前將足夠多的雜質(zhì)粒子注入等離子體,彈丸的飛行速度需達(dá)到約700m/s,雜質(zhì)粒子數(shù)達(dá)到2×1023[6],若將EPI運(yùn)用于ITER需要提高其發(fā)射參數(shù)以增加電樞及彈丸的飛行速度,其攜帶動(dòng)能也將極大增加,分離彈丸及偏轉(zhuǎn)方向的過程會(huì)變得更具風(fēng)險(xiǎn),即撞擊靶板瞬間產(chǎn)生的巨大沖擊力可能會(huì)產(chǎn)生電樞碎裂、彈丸偏離預(yù)定軌道、碎片進(jìn)入真空室等影響正常放電的嚴(yán)重后果。故為了克服現(xiàn)有系統(tǒng)的局限性,提高系統(tǒng)整體安全性及可靠性,實(shí)現(xiàn)電樞及彈丸的高發(fā)射動(dòng)能,于J-TEXT托卡馬克上設(shè)計(jì)了電磁彈丸注入系統(tǒng)(Electromagnetic Pellet Injection system, EMI)。EMI配備有電磁軌道加速及減速功能,能在加速完成后通過主動(dòng)電磁減速來平穩(wěn)分離電樞-彈丸,并在釋放彈丸后軟回收電樞。

    EPI所使用的C型電樞是尾翼接觸型電樞的代表之一,根據(jù)縱截面形狀可將這類電樞劃分為C型、H型等。尾翼接觸型電樞的結(jié)構(gòu)劃分如圖1和圖2所示,①~④分別代表頭部、喉部、肩部、尾翼四塊區(qū)域。

    電樞頭部指代運(yùn)動(dòng)方向的前端區(qū)域,主要功能是調(diào)整質(zhì)心并攜帶負(fù)載。C型電樞頭部不與軌道接觸,通過與前端非金屬材質(zhì)的導(dǎo)引連接,減小運(yùn)動(dòng)時(shí)的橫向位移。H型電樞頭部起到導(dǎo)引的功能,用于提高運(yùn)動(dòng)穩(wěn)定性,由于無上翹面,此處接觸壓力很小,發(fā)射時(shí)幾乎不通過電流[18]。

    圖1 C型電樞的主要結(jié)構(gòu)劃分

    圖2 H型電樞的主要結(jié)構(gòu)劃分

    電樞喉部指代電樞前端位于頭部后方的區(qū)域,是傳導(dǎo)電流的位置。電樞肩部指代喉部向后與尾部銜接的轉(zhuǎn)角區(qū)域,因此也是載流的關(guān)鍵位置。在發(fā)射中,尾翼與軌道保持過盈配合并受肩部拖動(dòng)向前運(yùn)動(dòng),故肩部為應(yīng)力集中的區(qū)域。

    電樞尾翼指代尾部與軌道直接接觸的區(qū)域。該處設(shè)置有一定過盈量,在裝填入膛后上翹面發(fā)生形變以產(chǎn)生初始接觸壓力。

    C型電樞具有質(zhì)量輕、通流能力強(qiáng)等優(yōu)點(diǎn),雖然這類電樞的研究結(jié)果較為豐富,但以往的研究重點(diǎn)在于加速過程的特性,其受力減速后運(yùn)動(dòng)過程的研究較少,對(duì)于將EMI用于托卡馬克破裂緩解產(chǎn)生的減速及軟回收的特殊需求,現(xiàn)有的尾翼接觸型電樞尚不能很好地滿足,主要原因如下:EPI所使用的C型電樞為具有薄尾翼的輕質(zhì)結(jié)構(gòu),但B探針信號(hào)波形的抖動(dòng)表明電樞在加速時(shí)可能存在不良電接觸的情況,且過薄的尾翼會(huì)降低電樞整體強(qiáng)度,若該電樞直接運(yùn)用于EMI可能發(fā)生尾翼磨損導(dǎo)致電接觸失效的情形。

    C型、H型等尾翼接觸型電樞的結(jié)構(gòu)功能相似,其中H型電樞的結(jié)構(gòu)與功能本質(zhì)上同安裝了導(dǎo)引的C型電樞一致,均通過后端的尾翼饋入電流,由前端的頭部/導(dǎo)引約束橫向位移,故存在電磁減速時(shí)接觸壓力被削弱的問題,以下將以C型電樞為代表,討論尾翼接觸型電樞的局限性。減速階段尾翼接觸型電樞的電磁接觸壓力如圖3所示,在導(dǎo)通減速電流后,尾翼流過的電流將與軌道電流部分反向,削弱原有同向電流時(shí)產(chǎn)生的電磁力,使尾翼與軌道貼合的電磁接觸壓力變小。同時(shí)在較長(zhǎng)距離的運(yùn)動(dòng)后,尾翼已被較大幅度地磨損,機(jī)械過盈帶來的初始接觸壓力大幅度減小,隨著減速電流的提高,電磁接觸壓力進(jìn)一步減小,導(dǎo)致樞軌接觸變差,最終很可能發(fā)生轉(zhuǎn)捩。

    圖3 減速階段尾翼接觸型電樞的電磁接觸壓力

    EMI的電源輸出電流及該電流作用下C型電樞的炮口電壓典型波形如圖4所示。炮口電壓的主要組成是電樞處兩端的電壓,可以反映電樞轉(zhuǎn)捩開始時(shí)刻,一般認(rèn)為炮口電壓幅值大于20V時(shí)開始發(fā)生轉(zhuǎn)捩,該電壓幅值越高表明電接觸性能越差[19-20]。由圖4可以看出在加速階段,炮口電壓幅值不超過10V,電樞與軌道能保持較好的電接觸,在減速電源觸發(fā)后,由于電樞前端的軌道通過電流,炮口電壓的組成由電樞兩端壓降擴(kuò)大為電樞壓降與軌道壓降之和,故此后轉(zhuǎn)捩的發(fā)生時(shí)刻表征為信號(hào)的高頻振蕩時(shí)刻而非幅值單純大于20V。當(dāng)=3.5ms時(shí)減速電源觸發(fā),此時(shí)炮口電壓幅值迅速爬升至350V并伴有約50V的高頻振蕩,表明電樞進(jìn)入減速階段后隨即發(fā)生轉(zhuǎn)捩;在隨后的5ms內(nèi),隨著減速電流上升沿的結(jié)束,炮口電壓平均幅值逐漸減小,但轉(zhuǎn)捩電弧持續(xù)存在。轉(zhuǎn)捩電弧溫度極高(超過3 000℃),將燒蝕組件并極大降低其壽命。產(chǎn)生的雜質(zhì)以等離子體形態(tài)為主,若進(jìn)入真空室則會(huì)引入非必要的雜質(zhì)粒子,嚴(yán)重影響托卡馬克裝置正常放電。

    圖4 C型電樞的發(fā)射結(jié)果

    另一方面,C型、H型電樞其自身結(jié)構(gòu)的重心位于靠前端的位置,減速過程中很可能會(huì)產(chǎn)生偏轉(zhuǎn)力矩使電樞發(fā)生橫向偏移,導(dǎo)致彈丸與電樞分離后同樣擁有橫向速度分量[21],最終沿著非預(yù)設(shè)路徑飛行,與軌道炮內(nèi)膛或真空室壁碰撞。

    基于上述原因,現(xiàn)有C型、H型等尾翼接觸型電樞尚不能完全滿足EMI的發(fā)射需求,需要設(shè)計(jì)能同時(shí)適應(yīng)加速及減速軟回收過程的新型電樞。本文首先介紹J-TEXT(Joint Texas Experimental Tokamak)托卡馬克裝置及EMI系統(tǒng);其次根據(jù)其實(shí)際需求及系統(tǒng)特點(diǎn),設(shè)計(jì)了符合需求的X型電樞,對(duì)該電樞進(jìn)行有限元仿真分析;最后裝載于EMI實(shí)驗(yàn)平臺(tái)上,測(cè)試其性能。

    1 J-TEXT托卡馬克及EMI系統(tǒng)

    1.1 J-TEXT托卡馬克裝置

    J-TEXT托卡馬克是典型的圓截面鐵心托卡馬克裝置,裝置實(shí)物如圖5所示。J-TEXT 托卡馬克的前身是美國(guó)德克薩斯大學(xué)奧斯汀分校聚變研究中心的TEXT-U裝置,后來由美國(guó)的能源部批準(zhǔn)贈(zèng)與了華中科技大學(xué)。J-TEXT是當(dāng)前教育部所有高校中唯一的中型托卡馬克裝置,主要以歐姆加熱為主,2019年實(shí)現(xiàn)了電子回旋加熱[22],裝置的主要運(yùn)行參數(shù)見表1。

    圖5 J-TEXT托卡馬克裝置

    表1 J-TEXT托卡馬克運(yùn)行參數(shù)

    Tab.1 The operating parameters of J-TEXT Tokamak

    1.2 J-TEXT托卡馬克上的EMI

    軌道炮是電磁炮的一種,利用電磁力在短時(shí)間內(nèi)將物體加速到高速度,結(jié)構(gòu)簡(jiǎn)單可靠,容易實(shí)現(xiàn)數(shù)km/s的超高發(fā)射速度。EMI以軌道炮為設(shè)計(jì)基礎(chǔ),具有高發(fā)射速度、低響應(yīng)時(shí)間、注入量大等優(yōu)點(diǎn)。常規(guī)軌道炮發(fā)射時(shí)利用空氣動(dòng)力分離電樞及載荷,并使用打擊靶板的方式將電樞轉(zhuǎn)向或減速[23-25],而EMI將運(yùn)行至托卡馬克的高真空內(nèi),無法實(shí)現(xiàn)氣動(dòng)分離與減速,且擊打靶板很可能產(chǎn)生電樞碎片,這些碎片若進(jìn)入真空室會(huì)嚴(yán)重影響放電,故EMI通過主動(dòng)減速環(huán)節(jié)達(dá)到分離載荷并軟回收電樞的目的。

    EMI的工作原理如圖6所示,由加速脈沖電源、減速脈沖電源、導(dǎo)軌、攜帶彈丸載荷的電樞組成。其中,為電源饋入電流,為電樞受到電磁力,為電場(chǎng)強(qiáng)度,為磁感應(yīng)強(qiáng)度。發(fā)射時(shí)導(dǎo)通加速脈沖電源,電樞經(jīng)過加速達(dá)到速度峰值,隨即觸發(fā)減速脈沖電源,使電樞受到減速力作用而釋放彈丸,彈丸脫離電樞后保持原速度飛行至托卡馬克真空室,而電樞被進(jìn)一步減速直至軟回收,避免進(jìn)入真空室損害裝置。為J-TEXT搭建的EMI試驗(yàn)平臺(tái)如圖7所示,部分參數(shù)見表2。

    圖6 EMI的工作原理

    不計(jì)阻力時(shí),電樞受到電磁力及運(yùn)動(dòng)方程為[26-28]

    圖7 EMI的試驗(yàn)平臺(tái)

    Fig.7 The core component of the EMI

    表2 EMI試驗(yàn)平臺(tái)的主要參數(shù)

    Tab.2 The main parameters of the EMI experimental platform

    2 X型電樞的設(shè)計(jì)

    電樞是EMI的核心部件,是電能-動(dòng)能轉(zhuǎn)換的樞紐,在發(fā)射過程中需承受數(shù)百kA脈沖電流、數(shù)百M(fèi)Pa應(yīng)力的作用,整體溫度也將升高至數(shù)百℃,其設(shè)計(jì)的優(yōu)劣決定EMI發(fā)射性能及破裂緩解效果的好壞,故需要設(shè)計(jì)匹配EMI的電樞。

    電樞需要在強(qiáng)電流、大電磁力、高溫升下保持與軌道的良好電接觸、通流能力和機(jī)械強(qiáng)度,故設(shè)計(jì)要求主要有三點(diǎn):①需具備良好的發(fā)射性能,能在惡劣環(huán)境下保持足夠的機(jī)械強(qiáng)度,發(fā)射過程中保證結(jié)構(gòu)的完整性,同時(shí)電樞與軌道配合能提供良好的滑動(dòng)電接觸性能;②電樞要滿足射向需求,控制結(jié)構(gòu)對(duì)稱度;③電樞要滿足內(nèi)膛壽命技術(shù)要求,能夠抑制和消除轉(zhuǎn)捩電弧,并控制材料的熔化磨損和沉積[17]。

    2.1 X型電樞各部分的參數(shù)設(shè)計(jì)

    X型電樞以現(xiàn)有典型固體電樞為基礎(chǔ),具有前翼、尾翼兩處接觸區(qū)域,提供了更強(qiáng)大的樞軌電接觸性能及運(yùn)動(dòng)穩(wěn)定性,其參數(shù)分布如圖9所示,電樞材料選用高熔化載流特征量的6061Al,以下為X型電樞主要結(jié)構(gòu)參數(shù)設(shè)計(jì)。

    圖9 X型電樞的主要參數(shù)分布

    1)肩部厚度0設(shè)計(jì)

    以電流在肩部處產(chǎn)生的電阻熱恰好熔化電樞材料為依據(jù),根據(jù)電流作用積分方法可以確定尾翼肩部厚度0值[29]。

    式中,為裕度系數(shù),0<<1;為6061Al的熔化載流特征量,=2.524×1016(A/m2)2·s;為絕緣間距,=12mm。裕度系數(shù)可以在0.5~0.9的范圍內(nèi)取值。結(jié)合電流波形,取峰值電流200kA、上升沿0.5ms,可計(jì)算出最小肩部厚度為4.37mm。

    2)前翼、尾翼長(zhǎng)度1設(shè)計(jì)

    以往的尾翼接觸型電樞,其尾翼長(zhǎng)度可與軌道間距保持一致,并根據(jù)接觸壓力分布的結(jié)果對(duì)1進(jìn)行微調(diào),使接觸壓力達(dá)到要求,而且壓力分布的面積要大,分布位置居于尾翼接觸面中心,由于X型電樞具有前后兩端尾翼,考慮到需要控制電樞的質(zhì)量及總長(zhǎng)度,1設(shè)置為12mm。

    3)喉部厚度2設(shè)計(jì)

    考慮到減速階段存在加速電流、減速電源的電流同時(shí)饋入電樞的情形,喉部厚度2應(yīng)大于尾翼肩部厚度0,同時(shí)為了盡量降低電樞整體質(zhì)量,2可取1.3~1.7倍0。

    4)過盈量3設(shè)計(jì)

    過盈量決定了初始接觸壓力的大小,初始接觸壓力的基本要求是大于0.01N/A。根據(jù)200kA的峰值電流可得到最小接觸壓力為2 000N,通過仿真接觸壓力分布、接觸面積大小,最終確定3為0.3mm。

    5)前翼、尾翼末端厚度4設(shè)計(jì)

    前翼、尾翼末端厚度4的取值主要考慮兩個(gè)方面的問題。一是尾翼熔化磨損速度,約0.5mm/m,要求電樞出膛時(shí)翼的末端仍然有一定厚度剩余,且機(jī)械強(qiáng)度能夠承受發(fā)射過程中樞軌界面的接觸壓力。EMI的軌道長(zhǎng)3m,故熔化磨損厚度約1.5mm,根據(jù)接觸壓力的仿真結(jié)果調(diào)整4,可確定4為3mm。

    2.2 X型電樞的有限元仿真結(jié)果

    綜合考慮計(jì)算與設(shè)計(jì)經(jīng)驗(yàn),得到X型電樞的主要參數(shù)見表3,使用有限元軟件的固體力學(xué)模塊仿真無電流饋入電樞時(shí)的初始應(yīng)力及接觸壓力分布,取電樞單邊過盈量為0.3mm,假設(shè)發(fā)生線彈性形變及塑性形變;使用電流模塊仿真其瞬態(tài)電流密度分布,其接觸壓力使用固體力學(xué)模塊仿真結(jié)果,仿真加速階段單獨(dú)饋入200kA加速電流、減速階段剩余50kA加速電流及150kA減速電流共同饋入電樞的情形。

    表3 X型電樞主要參數(shù)的設(shè)計(jì)結(jié)果

    Tab.3 Design results of the main parameters of the X-type armature

    電樞的初始應(yīng)力分布云圖如圖10所示,可以看到應(yīng)力集中在肩部、喉部及翼的前端,對(duì)于該材質(zhì)的電樞,當(dāng)應(yīng)力超過350MPa時(shí)將發(fā)生塑性形變,而結(jié)果顯示最大應(yīng)力值約為440MPa,故肩部、喉部及翼的表面出現(xiàn)少部分塑性形變,電樞主要區(qū)域仍處于線彈性形變的范圍內(nèi),且通過對(duì)電樞應(yīng)力分布進(jìn)行體積積分,除以電樞體積后獲得平均值可得出其平均應(yīng)力值約為230MPa。初始應(yīng)力分布表明電樞結(jié)構(gòu)設(shè)計(jì)合理,整體機(jī)械強(qiáng)度滿足發(fā)射要求,不會(huì)發(fā)生大范圍的塑性形變以減弱電接觸。

    圖10 初始應(yīng)力分布云圖

    電樞的初始接觸壓力分布云圖如圖11所示,可以看到接觸壓力主要分布在樞軌接觸面的中央,主要接觸面呈“沙漏形”。在橫向方向上,壓強(qiáng)由中心向兩側(cè)逐漸增大,在邊沿位置取得最大值,約為150MPa;在縱向方向上,接觸面中央存在壓強(qiáng)最大值,約為80MPa。由于越靠近翼的末端形變程度越大,但越靠近末端其厚度也在減少,導(dǎo)致相同形變下的接觸壓力減少,最終在縱向方向上的中央部位存在最大接觸壓強(qiáng)。另一方面,在橫向方向上,越靠近邊沿電樞形變量越大,故在邊沿處取得最大壓強(qiáng)。翼的總接觸壓力為4 465N,滿足接觸壓力設(shè)計(jì)指標(biāo)中初始機(jī)械過盈接觸壓力不小于2 000N的要求,且單邊過盈量0.3mm的取值能夠保證翼在經(jīng)運(yùn)動(dòng)磨損后仍與軌道有較大的機(jī)械接觸壓力,若取值過大則會(huì)增加接觸壓力,使得電樞與軌道間的摩擦力增加,導(dǎo)致發(fā)射效率降低、磨損增大。

    圖11 初始接觸壓力分布云圖

    電樞的電流密度分布云圖如圖12所示。由于電流將沿著電阻值最小的路徑通過,在加速階段電流主要經(jīng)過路徑更短的電樞尾翼,其電流密度與接觸壓力分布基本一致,電流主要通過接觸面的中央?yún)^(qū)域[30],最高值為4.07×109A/m2,此時(shí)電樞的前翼有少量電流通過,最高電流密度為2.23×109A/m2;在減速階段,由于加速電源仍有部分殘余能量,電樞的尾翼仍有較大電流通過,其電流密度最高值為3.52×109A/m2,而由于加速及減速電源的電流同時(shí)饋入電樞前翼,其最高值較加速階段尾翼的最高值稍有增加,為4.43×109A/m2。仿真表明電流密度分布合理,通流面積較大,滿足加速及減速階段的通流能力要求。

    3 X型電樞的發(fā)射試驗(yàn)

    攜帶彈丸的X型電樞實(shí)物如圖13所示,電樞頭部通過楔形結(jié)構(gòu)與聚碳酸酯材質(zhì)的導(dǎo)引連接,導(dǎo)引用于攜帶彈丸載荷并進(jìn)一步改善電樞運(yùn)動(dòng)軌跡,彈丸設(shè)計(jì)為子彈形以減少飛行時(shí)的空氣阻力。為了測(cè)試X型電樞的性能,使用液壓油缸驅(qū)動(dòng)推桿將電樞裝填入EMI的尾端,試驗(yàn)使用的主要測(cè)量手段見表4,分別開展電樞發(fā)射性能試驗(yàn)、電樞-彈丸分離性能試驗(yàn)。

    圖13 X型電樞及攜帶彈丸的實(shí)物

    表4 測(cè)量手段的設(shè)置

    Tab.4 The measurement method settings

    3.1 電樞的發(fā)射性能試驗(yàn)

    在加速電源饋入EMI后延時(shí)觸發(fā)減速電源,使電樞先加速至預(yù)定速度再減速軟回收,驗(yàn)證其在加速及減速階段的電接觸性與運(yùn)動(dòng)穩(wěn)定性,試驗(yàn)參數(shù)的設(shè)置見表5,本次試驗(yàn)結(jié)果如圖14所示,軟回收后的電樞實(shí)物如圖15所示。

    表5 試驗(yàn)參數(shù)的設(shè)置

    Tab.5 The experimental parameter settings

    圖14 X型電樞的發(fā)射性能試驗(yàn)結(jié)果

    圖15 軟回收后的電樞

    從電源輸出波形可以看到加速電源的輸出電流峰值達(dá)到了165kA,經(jīng)過4ms延時(shí)后觸發(fā)減速電源,其輸出電流的峰值約170kA。

    炮口電壓表征樞軌界面的電接觸狀態(tài),當(dāng)電壓超過約20V時(shí)開始發(fā)生轉(zhuǎn)捩,樞軌界面間產(chǎn)生電弧[19]。本試驗(yàn)的加速階段中,炮口電壓低于10V,電樞與軌道能保持較好的電接觸,當(dāng)=4.0ms時(shí)減速電源觸發(fā),炮口電壓幅值迅速爬升至300V,除觸發(fā)瞬間外其幅值無高頻振蕩,表明電樞進(jìn)入減速階段后沒有發(fā)生轉(zhuǎn)捩;在隨后的3ms內(nèi),炮口電壓平均幅值由峰值逐漸減小到20V以下,且電壓峰值持續(xù)時(shí)間僅為1ms,遠(yuǎn)低于C型電樞約2ms的持續(xù)時(shí)間。炮口電壓信號(hào)表明與C型等尾翼接觸型電樞相比,在減速階段,X型電樞的樞軌電接觸性能更為良好,僅在減速電源觸發(fā)瞬間有電弧產(chǎn)生[31]。

    B探針通過感應(yīng)電樞電流產(chǎn)生的磁場(chǎng)變化來測(cè)量電樞在膛內(nèi)的速度變化,且探針信號(hào)正負(fù)尖峰之間的過零點(diǎn)時(shí)刻為電流脈沖到達(dá)探針正對(duì)位置的時(shí)刻[32],故從圖14可以分析出電樞的位置-時(shí)間關(guān)系曲線,求導(dǎo)后得到電樞的速度-時(shí)間關(guān)系曲線,可以看到電樞經(jīng)過2.8ms的加速后達(dá)到速度峰值,最高速度約520m/s,隨后由于加速電源電流衰減及空氣阻力的提高,電樞速度有所降低,在=4.0ms時(shí)減速電源觸發(fā)開始主動(dòng)減速,最終在=6.6ms時(shí)電樞減速至0m/s。圖15電樞的實(shí)物表明經(jīng)過軟回收后的電樞能夠保持完整形貌,各結(jié)構(gòu)清晰、完整表明電樞各參數(shù)設(shè)計(jì)較合理,但其表面覆蓋有較多熔鋁層,說明電樞經(jīng)加速及減速兩次電流饋入受到了較大的熱負(fù)荷。

    3.2 電樞-彈丸分離性能試驗(yàn)

    在驗(yàn)證X型電樞加速、減速過程的良好電接觸性及受力驅(qū)動(dòng)性能后,需要進(jìn)一步開展電樞-彈丸分離試驗(yàn),測(cè)試其受力減速后能否穩(wěn)定釋放所攜帶的彈丸。

    考慮到過大的減速電源電流及轉(zhuǎn)捩電弧產(chǎn)生的熱負(fù)荷可能會(huì)熔化導(dǎo)引、燒蝕彈丸結(jié)構(gòu),進(jìn)而阻礙彈丸釋放,故在本次試驗(yàn)中采用了較低的減速電源參數(shù),圖16為本試驗(yàn)的電源輸出電流波形及電樞速度-時(shí)間曲線。加速電源的輸出電流峰值達(dá)到175kA,經(jīng)過4.5ms延時(shí)后觸發(fā)減速電源,其輸出電流的峰值約75kA。電樞經(jīng)過2.9ms的加速后達(dá)到405m/s的速度峰值;在=4.5ms時(shí)減速電源觸發(fā)對(duì)電樞進(jìn)行主動(dòng)減速,隨著減速電流的增大,電樞受到的減速電磁力增加;當(dāng)=4.7ms時(shí)減速電磁力超過分離臨界值使得電樞釋放彈丸,此時(shí)電樞、彈丸速度為358m/s;隨后彈丸保持該速度飛行至出膛,而電樞進(jìn)一步減速至239m/s。

    圖16 X型電樞-彈丸分離性能試驗(yàn)結(jié)果

    電樞及彈丸出膛后的飛行畫面如圖17所示,彈丸結(jié)構(gòu)完整、清晰,射向的準(zhǔn)直性良好,能夠在將來的托卡馬克破裂緩解試驗(yàn)中實(shí)現(xiàn)精準(zhǔn)注入粒子;出膛畫面顯示電樞結(jié)構(gòu)完整,未產(chǎn)生嚴(yán)重的轉(zhuǎn)捩電弧燒蝕,進(jìn)一步說明X型電樞的樞軌界面電接觸良好;電樞出膛后,由于電源的能量未被完全耗盡,殘余電流以炮口電弧的形式釋放,在電樞的尾部產(chǎn)生較大電弧光及雜質(zhì)等離子體。炮口電弧會(huì)燒蝕軌道,降低軌道預(yù)期壽命,并產(chǎn)生非必要的雜質(zhì)等離子體,因此需要探索抑制炮口電弧的手段。

    圖17 高速相機(jī)拍攝的電樞、彈丸出膛后的飛行軌跡

    4 結(jié)論

    托卡馬克作為磁約束聚變的裝置,在運(yùn)行過程中可能發(fā)生等離子體破裂的問題,而破裂所釋放的巨大能量會(huì)對(duì)裝置造成嚴(yán)重破壞。為了抑制等離子體破裂的負(fù)面影響,保障人員及裝置的安全,為ITER及J-TEXT托卡馬克設(shè)計(jì)了新一代的破裂緩解系統(tǒng)EMI,與以往的破裂緩解系統(tǒng)MGI和SPI相比,EMI具有低響應(yīng)時(shí)間、高注入速度、注入量大等優(yōu)點(diǎn)。

    本文通過分析C型、H型等尾翼接觸型固體電樞的結(jié)構(gòu)與功能,結(jié)合EMI對(duì)電樞發(fā)射性能的需求,設(shè)計(jì)了具有更強(qiáng)運(yùn)動(dòng)穩(wěn)定性、在加速及減速過程均能保持良好電接觸的X型電樞,克服了常規(guī)固體電樞在減速過程中接觸壓力減小、電接觸性能變差的局限性。

    仿真結(jié)果顯示,X型電樞的接觸壓力達(dá)到要求、壓強(qiáng)分布均勻,在尾翼中部達(dá)到最大值且接觸面較大;電流密度與接觸壓力分布基本一致,主要通流區(qū)域?yàn)榻佑|面中央?yún)^(qū)域,滿足通流能力要求;電樞應(yīng)力集中在肩部及喉部,滿足發(fā)射過程的機(jī)械性能要求。

    EMI的發(fā)射性能試驗(yàn)結(jié)果表明,在峰值165kA的脈沖電流作用下,X型電樞能加速至約520m/s,經(jīng)過4ms延時(shí)后觸發(fā)減速電源,其電流峰值約170kA,在=6.6ms時(shí)減速至0m/s。本試驗(yàn)的炮口電壓表明,僅在減速電源觸發(fā)瞬間有轉(zhuǎn)捩電弧產(chǎn)生,樞軌界面電接觸性能良好,且經(jīng)過軟回收后的電樞能夠保持完整形貌,各結(jié)構(gòu)清晰、完整。EMI的X型電樞-彈丸分離性能試驗(yàn)表明,減速電源觸發(fā)后,電樞能與彈丸穩(wěn)定分離,釋放后的彈丸速度達(dá)358m/s,且該速度能在提高加速電流后進(jìn)一步增加。

    EMI為下一代大型托卡馬克裝置上的破裂緩解提供了新型高效雜質(zhì)注入緩解系統(tǒng),本文提出的X型電樞具有良好的發(fā)射性能,為EMI貢獻(xiàn)了一種優(yōu)良的新型電樞結(jié)構(gòu)及設(shè)計(jì)方法。EMI在發(fā)射過程中仍有部分炮口電弧及額外雜質(zhì)產(chǎn)生,這會(huì)對(duì)組件壽命及托卡馬克正常放電產(chǎn)生影響,故未來將針對(duì)炮口電弧抑制及雜質(zhì)控制展開研究。

    [1] 石秉仁. 磁約束聚變?cè)砼c實(shí)踐[M]. 北京: 原子能出版社, 1999.

    [2] 何開輝, 潘傳紅, 馮開明. 托卡馬克等離子體大破裂及防治綜述[J]. 中國(guó)核科技報(bào)告, 2002, 12: 117-129.

    He Kaihui, Pan Chuanhong, Feng Kaiming. An overview on plasma disruption mitigation and avoidance in Tokamak[J]. China Nuclear Science and Technology Report, 2002, 12: 117-129.

    [3] Arnoux G, Loarte A, Riccardo V, et al. Heat loads on plasma facing components during disruptions on JET[J]. Nuclear Fusion, 2009, 49(8): 5038-5046.

    [4] Arnoux G, Bazylev B, Lehnen M, et al. Heat load measurements on the JET first wall during disruptions[J]. Journal of Nuclear Materials, 2011, 415(1): S817-S820.

    [5] Matthews G F, Bazylev B, Baron-Wiechec A, et al. Melt damage to the JET ITER-like wall and divertor[J]. Physica Scripta, 2016, 167(1): 4070-4078.

    [6] Lehnen M, Campbell D J, Maruyama S, et al. Control, detection and mitigation of disruptions on ITER[C]//2015 IEEE 26th Symposium on Fusion Engineering, Austin, TX, USA, 2015: 1-8.

    [7] Hollmann E M, Aleynikov P B, Fül?p T, et al. Status of research toward the ITER disruption mitigation system[J]. Physics of Plasmas, 2014, 22(2): 1802-1819.

    [8] Taylor P L, Kellman A G, Evans T E, et al. Disruption mitigation studies in DIII-D[J]. Physics of Plasmas, 1999, 6(5): 1872-1879.

    [9] Pautasso G, Buchl K, Fuchs J C, et al. Use of impurity pellets to control energy dissipation during disruption[J]. Nuclear Fusion, 1996, 36(10): 1291-1297.

    [10] Yu Yaowei, Kim Y O, Kim H K, et al. First results on disruption mitigation by massive gas injection in Korea superconducting Tokamak advanced research[J]. Review of Scientific Instruments, 2012, 83(12): 3509-3515.

    [11] Combs S K, Baylor L R, Fisher P W, et al. ORNL mock-up tests of inside launch pellet injection on JET and LHD[J]. Fusion Engineering and Design, 2001, 58(11): 343-347.

    [12] Hollmann E M, Jernigan T C, Parks P B, et al. Measurements of injected impurity assimilation during massive gas injection experiments in DIII-D[J]. Nuclear Fusion, 2008, 48(11): 5007-5019.

    [13] Raman R, Lay W S, Jarboe T R, et al. Electromagnetic particle injector for fast time response disruption mitigation in tokamaks[J]. Nuclear Fusion, 2019, 59(1): 6021-6032.

    [14] 湯熠. J-TEXT上基于高速相機(jī)研究MGI產(chǎn)生的冷前鋒傳播過程[M]. 武漢: 華中科技大學(xué), 2015.

    [15] Baylor L R, Barbier C C, Carmichael J R, et al. Disruption mitigation system developments and design for ITER[J]. Fusion Science and Technology, 2015, 68(2): 211-215.

    [16] Raman R, Jarboe T R, Menard J E, et al. Fast time response electromagnetic disruption mitigation concept[J]. Fusion Science and Technology, 2015, 68(4): 797-805.

    [17] Raman R, Lunsford R, Clauser C F, et al. Prototype tests of the electromagnetic particle injector-2 for fast time response disruption mitigation in tokamaks[J]. Nuclear Fusion, 2021, 61(12): 6034-6048.

    [18] 黨晟罡. 幾種典型固體電樞的形狀設(shè)計(jì)與接觸特性研究[D]. 秦皇島: 燕山大學(xué), 2016.

    [19] (美)理查德·埃斯特里·馬歇爾. 電磁軌道炮的科學(xué)與技術(shù)[M]. 北京: 兵器工業(yè)出版社, 2006.

    [20] 李湘平, 魯軍勇, 馮軍紅, 等. 電磁發(fā)射彈丸飛行彈道仿真[J]. 國(guó)防科技大學(xué)學(xué)報(bào), 2019, 41(4): 25-32.

    Li Xiangping, Lu Junyong, Feng Junhong, et al. Simulation of flight ballistic of electromagnetic launch projectile[J]. Journal of National University of Defense Technology, 2019, 41(4): 25-32.

    [21] 錢杏芳. 導(dǎo)彈飛行力學(xué)[M]. 北京: 北京工業(yè)學(xué)院出版社, 1987.

    [22] 佟瑞海. J-TEXT托卡馬克基于大量雜質(zhì)氣體注入的破裂緩解實(shí)驗(yàn)研究[D]. 武漢: 華中科技大學(xué), 2019.

    [23] 楊啟仁, 徐直軍. 脫殼動(dòng)力學(xué)[M]. 北京: 國(guó)防工業(yè)出版社, 1996.

    [24] 古剛, 李宣. 一體化彈丸彈體和彈托分離特性研究[J]. 艦船科學(xué)技術(shù), 2020, 42(15): 32-37.

    Gu Gang, Li Xuan. Research on separation characteristic of intergration launch projectile's body[J]. Ship Science and Technology, 2020, 42(15): 32-37.

    [25] Holcomb C T, Jarboe T R, Mattick A T, et al. Nonperturbing field profile measurements of a sustained spheromak[J]. Review of Scientific Instruments, 2001, 72(1): 1054-1058.

    [26] 唐波, 徐英桃, 栗保明. 電磁軌道炮被動(dòng)式炮口消弧裝置的理論分析與優(yōu)化設(shè)計(jì)[J]. 彈道學(xué)報(bào), 2016, 28(4): 62-67.

    Tang Bo, Xu Yingtao, Li Baoming. Theoretical analysis and optimization design of passive muzzle arc suppression device in a railgun[J]. Journal of Ballistics, 2016, 28(4): 62-67.

    [27] 劉旭堃, 于歆杰, 劉秀成. 電容儲(chǔ)能型脈沖電源分時(shí)分段觸發(fā)策略自動(dòng)計(jì)算方法[J]. 電工技術(shù)學(xué)報(bào), 2016, 31(11): 186-193.

    Liu Xukun, Yu Xinjie, Liu Xiucheng. An automatic calculation method for the triggering strategy of the capacitive pulsed-power supply[J]. Transactions of China Electrotechnical Society, 2016, 31(11): 186-193.

    [28] 龔晨, 于歆杰, 劉秀成. 電容儲(chǔ)能型軌道炮連續(xù)發(fā)射系統(tǒng)設(shè)計(jì)與仿真[J]. 電工技術(shù)學(xué)報(bào), 2013, 28(S2): 111-115, 121.

    Gong Chen, Yu Xinjie, Liu Xiucheng. Continuous emission scheme and its simulation for capacitor-based railgun system[J]. Transactions of China Electrotechnical Society, 2013, 28(S2): 111-115, 121.

    [29] Chen Lixue, He Junjia, Xia Shengguo, et al. Some key parameters of monolithic C-type armature in rectangular caliber railgun[J]. IEEE Transactions on Plasma Science, 2017, 45(7): 1465-1469.

    [30] 阮景煇, 陳立學(xué), 夏勝國(guó), 等. 電磁軌道炮電流分布特性研究綜述[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(21): 4423-4431.

    Ruan Jinghui, Chen Lixue, Xia Shengguo, et al. A review of current distribution in electromagnetic railguns[J]. Transactions of China Electrotechnical Society, 2020, 35(21): 4423-4431.

    [31] 王振春, 鮑志勇, 曹海要, 等. 增強(qiáng)型電磁軌道炮電樞軌道接觸特性研究[J]. 兵工學(xué)報(bào), 2018, 39(3): 451-456.

    Wang Zhenchun, Bao Zhiyong, Cao Haiyao, et al. Research on contact characteristics of armature and rail in augmented electromagnetic railgun[J]. Acta Armamentarii, 2018, 39(3): 451-456.

    [32] 劉福才, 李歡, 王大正, 等. 基于B探針的膛內(nèi)電樞速度測(cè)量系統(tǒng)研究[J]. 兵工學(xué)報(bào), 2014, 35(6): 762-768.

    Liu Fucai, Li Huan, Wang Dazheng, et al. Research on velocity measurement of solid armature in electromagnetic launch based on B-dot probes[J]. Acta Armamentarii, 2014, 35(6): 762-768.

    Design of X-Type Armature of Electromagnetic Pellet Injection System on J-TEXT Tokamak

    Chen Zhongyong1,2Zhang Weikang1,2Tang Junhui3Li Feng1,2Xia Shengguo3

    (1. International Joint Research Laboratory of Magnetic Confinement Fusion and Plasma Physics School of Electrical and Electronic Engineering Huazhong University of Science and Technology Wuhan 430074 China 2. State Key Laboratory of Advanced Electromagnetic Engineering and Technology School of Electrical and Electronic Engineering Huazhong University of Science and Technology Wuhan 430074 China 3. Key Laboratory of Pulsed Power Technology Ministry of Education Huazhong University of Science and Technology Wuhan 430074 China)

    The International Thermonuclear Experimental Reactor (ITER) program is the largest international scientific and technological cooperation project that China has participated in, which goal is to verify the scientific and technical feasibility of using large tokamak devices to obtain fusion energy. The major plasma disruption is the biggest threat to the safe operation of ITER, resulting in the damage of the device. Therefore, the disruption mitigation is a key scientific and technical problem to be solved urgently. At present, the basic strategy for disruption mitigation is to actively inject a large number of particles, but the existing disruption mitigation systems cannot fully meet the requirements of ITER. By analyzing the characteristics of existing systems, this paper introduces a new generation of Electromagnetic Pellet Injection system (EMI) for disruption mitigation on J-TEXT tokamak. The system uses electromagnetic force to launch pellets, which can effectively improve the injection speed and reduced response time, overcoming the limitations of other systems. The armature is the core component of EMI. By introducing the structure and function of the tail contact armature, it is shown that when this type of armature is used in EMI, there are insufficient armature-rail electrical contact performance and motion stability in the deceleration stage. According to the special requirements of the armature performance in the EMI deceleration stage, a new solid armature with an X-shaped structure is designed. The simulation results show that the electromagnetic and mechanical properties of the armature meet the launch requirements, and the electrical contact characteristics are brilliant. In the launch performance test, the armature is accelerated to 520m/s and then actively decelerated to 0m/s; in the armature-pellet separation test, the stable separation of the armature and pellet is achieved, and the flight velocity of pellet is 358m/s, and this velocity can be further increased with the increase of acceleration energy. EMI provides an advanced scheme for efficient disruption mitigation on tokamak, and the X-type armature proposed in this paper has good launch performance which provides an excellent new type of armature structure and design method for EMI.

    Tokamak, disruption mitigation system, electromagnetic launch, railgun, solid armature

    10.19595/j.cnki.1000-6753.tces.221242

    TM33; TL67

    國(guó)家重點(diǎn)研發(fā)計(jì)劃(2019YFE03010004)和國(guó)家自然科學(xué)基金(51821005)資助項(xiàng)目。

    2022-06-29

    2022-08-08

    陳忠勇 男,1978年生,教授,博士生導(dǎo)師,研究方向?yàn)橥锌R克等離子體破裂及破裂緩解。E-mail:zychen@hust.edu.cn(通信作者)

    張維康 男,1998年生,碩士研究生,研究方向?yàn)橥锌R克破裂緩解及電磁發(fā)射。E-mail:m202071656@hust.edu.cn

    (編輯 赫蕾)

    美女xxoo啪啪120秒动态图| 特大巨黑吊av在线直播| 亚洲人成网站在线观看播放| 交换朋友夫妻互换小说| 一本色道久久久久久精品综合| 色视频在线一区二区三区| 久久久久精品性色| 国产亚洲av片在线观看秒播厂| 国产精品.久久久| 亚洲精品自拍成人| 男女边吃奶边做爰视频| 久久精品国产亚洲网站| 人妻一区二区av| 午夜福利在线观看免费完整高清在| 亚洲欧洲国产日韩| av在线播放精品| a级毛色黄片| 男人添女人高潮全过程视频| 天天操日日干夜夜撸| 日韩三级伦理在线观看| 老司机影院成人| 国产高清三级在线| 亚洲av日韩在线播放| 亚洲一级一片aⅴ在线观看| 亚洲成人手机| 国产爽快片一区二区三区| 精品久久久噜噜| 超碰97精品在线观看| 亚洲av免费高清在线观看| 久久久国产一区二区| 国产国拍精品亚洲av在线观看| 免费观看的影片在线观看| 丝瓜视频免费看黄片| 一个人看视频在线观看www免费| 日产精品乱码卡一卡2卡三| 男女边吃奶边做爰视频| 国产乱来视频区| 久久毛片免费看一区二区三区| 男的添女的下面高潮视频| 国产欧美另类精品又又久久亚洲欧美| 日韩免费高清中文字幕av| 嘟嘟电影网在线观看| 日日摸夜夜添夜夜添av毛片| 少妇人妻一区二区三区视频| 观看美女的网站| 免费黄色在线免费观看| 国产成人一区二区在线| 99热国产这里只有精品6| 嫩草影院新地址| 妹子高潮喷水视频| 日韩,欧美,国产一区二区三区| 中文欧美无线码| 亚洲精品国产色婷婷电影| 国产精品99久久99久久久不卡 | av在线app专区| 亚洲成色77777| 亚洲欧美成人综合另类久久久| av国产精品久久久久影院| 另类精品久久| 九草在线视频观看| 日日撸夜夜添| 麻豆成人午夜福利视频| 亚洲精品一二三| 一级,二级,三级黄色视频| 亚洲欧洲日产国产| 欧美丝袜亚洲另类| 亚洲中文av在线| 日韩伦理黄色片| 免费观看无遮挡的男女| 纯流量卡能插随身wifi吗| 十八禁网站网址无遮挡 | 久久精品国产a三级三级三级| 女性生殖器流出的白浆| 99久久中文字幕三级久久日本| 成人无遮挡网站| 高清在线视频一区二区三区| 高清黄色对白视频在线免费看 | 欧美老熟妇乱子伦牲交| 69精品国产乱码久久久| videossex国产| 久久久国产欧美日韩av| 国内少妇人妻偷人精品xxx网站| 色哟哟·www| 亚洲国产色片| 亚洲在久久综合| av福利片在线| 亚洲精品国产av成人精品| 久久久a久久爽久久v久久| 精品熟女少妇av免费看| 一边亲一边摸免费视频| 亚洲人与动物交配视频| 精品少妇黑人巨大在线播放| 亚洲av.av天堂| 99九九线精品视频在线观看视频| 综合色丁香网| 午夜老司机福利剧场| 久久精品国产亚洲网站| 老熟女久久久| 99久久精品一区二区三区| 黄色一级大片看看| 国产免费福利视频在线观看| 日本av免费视频播放| 国产成人免费观看mmmm| 欧美少妇被猛烈插入视频| 亚洲国产最新在线播放| 国产乱来视频区| 免费少妇av软件| 国产欧美亚洲国产| 一级,二级,三级黄色视频| 日韩中字成人| 亚洲欧洲国产日韩| www.色视频.com| 亚洲内射少妇av| 精品一区在线观看国产| 9色porny在线观看| 国产成人免费无遮挡视频| 国产片特级美女逼逼视频| 黄色视频在线播放观看不卡| 青春草国产在线视频| 国产乱人偷精品视频| 欧美老熟妇乱子伦牲交| 成人免费观看视频高清| 汤姆久久久久久久影院中文字幕| 亚洲av男天堂| 9色porny在线观看| 亚洲精品乱码久久久v下载方式| 人人妻人人澡人人爽人人夜夜| 女人精品久久久久毛片| 亚洲人与动物交配视频| av天堂中文字幕网| 国产亚洲午夜精品一区二区久久| 色视频www国产| 国产欧美亚洲国产| 少妇的逼好多水| 麻豆精品久久久久久蜜桃| 不卡视频在线观看欧美| 国产一区亚洲一区在线观看| 另类精品久久| 人妻夜夜爽99麻豆av| 亚洲精品国产色婷婷电影| 婷婷色av中文字幕| 国产免费一级a男人的天堂| 激情五月婷婷亚洲| 99久久人妻综合| 精品熟女少妇av免费看| 少妇裸体淫交视频免费看高清| 亚洲av在线观看美女高潮| 亚洲电影在线观看av| 欧美精品国产亚洲| 久久久久久久国产电影| 在线观看av片永久免费下载| 大片免费播放器 马上看| 在线观看人妻少妇| 午夜日本视频在线| 国产免费一区二区三区四区乱码| 亚洲国产毛片av蜜桃av| 日韩av在线免费看完整版不卡| 在线观看美女被高潮喷水网站| 精品99又大又爽又粗少妇毛片| 日本免费在线观看一区| 少妇裸体淫交视频免费看高清| 老熟女久久久| 性高湖久久久久久久久免费观看| 亚洲精品,欧美精品| 综合色丁香网| 国产一级毛片在线| 国产男女超爽视频在线观看| 欧美性感艳星| 一个人免费看片子| 日本色播在线视频| 欧美日韩av久久| 欧美激情国产日韩精品一区| 亚洲国产欧美在线一区| 亚洲精品中文字幕在线视频 | 亚洲精品国产成人久久av| 亚洲av综合色区一区| 免费av不卡在线播放| 久久人人爽人人爽人人片va| 国产91av在线免费观看| 最新的欧美精品一区二区| 久久国产亚洲av麻豆专区| 国产日韩欧美在线精品| 免费观看a级毛片全部| 一区二区三区乱码不卡18| 少妇的逼水好多| 久久精品久久久久久噜噜老黄| 热99国产精品久久久久久7| 欧美老熟妇乱子伦牲交| 精品久久久久久电影网| 婷婷色av中文字幕| 日本-黄色视频高清免费观看| 精品一区在线观看国产| 国产在线一区二区三区精| 在线观看av片永久免费下载| 国模一区二区三区四区视频| 免费看不卡的av| 国产一区二区在线观看日韩| 99热6这里只有精品| 成人免费观看视频高清| 久久影院123| 乱人伦中国视频| 亚洲精品视频女| 黄色毛片三级朝国网站 | 人人澡人人妻人| 十八禁高潮呻吟视频 | 国产精品久久久久久久久免| 亚洲国产精品一区二区三区在线| 国产黄片视频在线免费观看| 久久这里有精品视频免费| 最黄视频免费看| 99久久精品国产国产毛片| 欧美老熟妇乱子伦牲交| 一边亲一边摸免费视频| 91成人精品电影| 国产黄色免费在线视频| 欧美 亚洲 国产 日韩一| 色5月婷婷丁香| 丰满人妻一区二区三区视频av| 免费在线观看成人毛片| 国产精品一区二区性色av| 久久人人爽av亚洲精品天堂| 美女中出高潮动态图| 观看av在线不卡| 我要看黄色一级片免费的| 在线精品无人区一区二区三| 日本猛色少妇xxxxx猛交久久| 在线亚洲精品国产二区图片欧美 | 国产亚洲91精品色在线| 国产一级毛片在线| 亚洲精品久久久久久婷婷小说| 亚洲成人手机| 日本黄色日本黄色录像| 在线观看国产h片| 新久久久久国产一级毛片| 日韩人妻高清精品专区| 内射极品少妇av片p| 国产精品无大码| 久久99热这里只频精品6学生| 国产精品久久久久久精品古装| 欧美激情国产日韩精品一区| 精品久久久久久久久亚洲| 日日摸夜夜添夜夜爱| 久久精品国产a三级三级三级| 插逼视频在线观看| 国产国拍精品亚洲av在线观看| 亚洲欧美日韩另类电影网站| 在线播放无遮挡| 免费人妻精品一区二区三区视频| 亚洲av电影在线观看一区二区三区| 国产精品一二三区在线看| kizo精华| 国产精品无大码| 99久久综合免费| av黄色大香蕉| 亚洲av日韩在线播放| 我要看日韩黄色一级片| 嘟嘟电影网在线观看| 国内少妇人妻偷人精品xxx网站| 在线观看免费高清a一片| 久久这里有精品视频免费| 亚洲伊人久久精品综合| 欧美3d第一页| 女人久久www免费人成看片| 亚洲美女黄色视频免费看| 亚洲一级一片aⅴ在线观看| 毛片一级片免费看久久久久| 99久久人妻综合| 欧美激情国产日韩精品一区| 亚洲欧美精品自产自拍| 国产成人91sexporn| 男男h啪啪无遮挡| 久久精品国产亚洲av天美| 久久久久久久久久人人人人人人| av免费在线看不卡| 2021少妇久久久久久久久久久| 欧美人与善性xxx| 大片电影免费在线观看免费| 欧美精品高潮呻吟av久久| 最新中文字幕久久久久| 国产精品久久久久久精品古装| 在线观看av片永久免费下载| 亚洲国产色片| 高清黄色对白视频在线免费看 | 亚洲国产精品一区二区三区在线| 91精品国产国语对白视频| 亚洲精华国产精华液的使用体验| 日本猛色少妇xxxxx猛交久久| videos熟女内射| 成人黄色视频免费在线看| 人妻少妇偷人精品九色| 欧美xxⅹ黑人| 国产在视频线精品| 亚洲成人av在线免费| 久久久久久久精品精品| 国产又色又爽无遮挡免| 欧美少妇被猛烈插入视频| 午夜影院在线不卡| 日韩大片免费观看网站| 成年av动漫网址| 国产午夜精品久久久久久一区二区三区| 亚洲欧美中文字幕日韩二区| 天天操日日干夜夜撸| 国产在视频线精品| 国产熟女欧美一区二区| 91精品国产国语对白视频| 亚洲,欧美,日韩| 欧美国产精品一级二级三级 | 国产综合精华液| 日韩精品免费视频一区二区三区 | 国产有黄有色有爽视频| videossex国产| 黄色一级大片看看| 黄色怎么调成土黄色| 久久这里有精品视频免费| 国产精品人妻久久久影院| 十八禁高潮呻吟视频 | 一级爰片在线观看| 爱豆传媒免费全集在线观看| 一个人免费看片子| 日韩三级伦理在线观看| 亚洲第一av免费看| 亚州av有码| av黄色大香蕉| 十分钟在线观看高清视频www | 国产精品无大码| 性色avwww在线观看| 久久这里有精品视频免费| 久久精品久久久久久噜噜老黄| 精品人妻熟女av久视频| 国产成人精品福利久久| 精品久久国产蜜桃| 纯流量卡能插随身wifi吗| 一级,二级,三级黄色视频| 一级毛片我不卡| 在线精品无人区一区二区三| 黑人高潮一二区| 91精品国产九色| 人妻夜夜爽99麻豆av| 亚洲,一卡二卡三卡| 夜夜爽夜夜爽视频| av播播在线观看一区| 久久 成人 亚洲| 乱人伦中国视频| 亚洲精品久久久久久婷婷小说| 国语对白做爰xxxⅹ性视频网站| 午夜免费男女啪啪视频观看| www.色视频.com| 三级国产精品片| 国产黄片美女视频| 夫妻午夜视频| 欧美丝袜亚洲另类| 在线天堂最新版资源| 国产亚洲午夜精品一区二区久久| 超碰97精品在线观看| 欧美日韩综合久久久久久| 免费av中文字幕在线| 男人狂女人下面高潮的视频| 2022亚洲国产成人精品| 韩国高清视频一区二区三区| 国产精品偷伦视频观看了| a 毛片基地| av网站免费在线观看视频| 一本一本综合久久| 纵有疾风起免费观看全集完整版| 欧美xxxx性猛交bbbb| 亚洲欧美精品自产自拍| 欧美日韩av久久| 美女大奶头黄色视频| 成年人午夜在线观看视频| 我的老师免费观看完整版| 亚洲综合精品二区| 亚洲国产精品一区二区三区在线| 男人爽女人下面视频在线观看| 日本免费在线观看一区| 一本一本综合久久| 在线观看av片永久免费下载| 青春草国产在线视频| 自拍偷自拍亚洲精品老妇| 欧美精品亚洲一区二区| 国产黄色免费在线视频| 国产精品久久久久久精品古装| 深夜a级毛片| 一级av片app| 婷婷色综合大香蕉| 99久久人妻综合| 99热网站在线观看| 日韩视频在线欧美| 激情五月婷婷亚洲| 国产精品99久久久久久久久| 一区二区三区四区激情视频| 国产成人一区二区在线| 啦啦啦在线观看免费高清www| 极品人妻少妇av视频| 毛片一级片免费看久久久久| 丁香六月天网| 久久久久久久久久人人人人人人| 日韩大片免费观看网站| 国产中年淑女户外野战色| 欧美bdsm另类| 国产美女午夜福利| 国产成人精品福利久久| 久久久久人妻精品一区果冻| av在线app专区| 亚洲精品乱码久久久久久按摩| 色视频在线一区二区三区| 精品亚洲成a人片在线观看| 欧美日韩视频精品一区| 欧美日韩国产mv在线观看视频| 国产一区二区三区综合在线观看 | 久久久久久久久久久免费av| av一本久久久久| 十八禁高潮呻吟视频 | 亚洲精品456在线播放app| 国产女主播在线喷水免费视频网站| av国产精品久久久久影院| 高清毛片免费看| 精华霜和精华液先用哪个| 成人二区视频| 国产黄频视频在线观看| 国产亚洲一区二区精品| 日韩欧美 国产精品| 欧美 亚洲 国产 日韩一| 美女脱内裤让男人舔精品视频| 一二三四中文在线观看免费高清| 国产成人一区二区在线| 丰满迷人的少妇在线观看| 国产免费又黄又爽又色| 亚洲综合色惰| 啦啦啦啦在线视频资源| 在线天堂最新版资源| av天堂久久9| 久久这里有精品视频免费| 在线观看三级黄色| tube8黄色片| 精品一区二区免费观看| 国产美女午夜福利| 我的老师免费观看完整版| 日本欧美国产在线视频| 插逼视频在线观看| 在线观看免费高清a一片| 国产精品不卡视频一区二区| 一区二区av电影网| 三上悠亚av全集在线观看 | 亚洲国产成人一精品久久久| 黄片无遮挡物在线观看| 欧美丝袜亚洲另类| 午夜福利在线观看免费完整高清在| 欧美性感艳星| 日韩欧美一区视频在线观看 | 涩涩av久久男人的天堂| 如日韩欧美国产精品一区二区三区 | 尾随美女入室| 91午夜精品亚洲一区二区三区| 国产精品人妻久久久影院| 亚洲图色成人| 看非洲黑人一级黄片| 特大巨黑吊av在线直播| 国产一区二区三区av在线| 一级av片app| 免费高清在线观看视频在线观看| 午夜福利影视在线免费观看| 乱人伦中国视频| 亚洲真实伦在线观看| 少妇熟女欧美另类| 亚洲激情五月婷婷啪啪| a 毛片基地| 纵有疾风起免费观看全集完整版| 国产精品久久久久成人av| 夫妻性生交免费视频一级片| 少妇的逼好多水| av专区在线播放| 日本猛色少妇xxxxx猛交久久| 成年女人在线观看亚洲视频| 一本—道久久a久久精品蜜桃钙片| 亚洲av综合色区一区| a级毛片在线看网站| 一级毛片我不卡| 国产一区二区三区综合在线观看 | 欧美人与善性xxx| 日本av手机在线免费观看| 久久精品熟女亚洲av麻豆精品| 高清av免费在线| 久热这里只有精品99| 欧美另类一区| 在线播放无遮挡| 国产视频首页在线观看| 一级毛片aaaaaa免费看小| 777米奇影视久久| 老女人水多毛片| 日韩三级伦理在线观看| 国产高清三级在线| 黑人高潮一二区| 久久久久久久亚洲中文字幕| 亚洲av福利一区| 少妇被粗大的猛进出69影院 | 大又大粗又爽又黄少妇毛片口| 国产精品人妻久久久影院| 中文字幕制服av| 日韩不卡一区二区三区视频在线| 久久久久久久大尺度免费视频| 精华霜和精华液先用哪个| 丝袜在线中文字幕| 欧美日韩精品成人综合77777| 中文字幕免费在线视频6| 极品人妻少妇av视频| 一个人免费看片子| 亚洲精品乱久久久久久| 午夜影院在线不卡| 精品久久国产蜜桃| a级一级毛片免费在线观看| 日韩制服骚丝袜av| 亚洲国产欧美日韩在线播放 | 超碰97精品在线观看| 国产真实伦视频高清在线观看| 国产精品久久久久久久久免| 日韩中字成人| 亚洲国产成人一精品久久久| 国产伦精品一区二区三区视频9| 久久久久久久久久久丰满| 亚洲av男天堂| 成人美女网站在线观看视频| 日本猛色少妇xxxxx猛交久久| 国产午夜精品一二区理论片| 一级二级三级毛片免费看| 日本黄色日本黄色录像| 热re99久久国产66热| 插逼视频在线观看| 亚洲国产色片| freevideosex欧美| 国产精品久久久久久精品古装| 看非洲黑人一级黄片| av线在线观看网站| 成人国产av品久久久| 国产国拍精品亚洲av在线观看| 在线观看免费视频网站a站| 国产精品三级大全| 91精品一卡2卡3卡4卡| 精品少妇久久久久久888优播| 久久久久国产网址| 热re99久久国产66热| 黑丝袜美女国产一区| 最近手机中文字幕大全| 人妻夜夜爽99麻豆av| 少妇人妻精品综合一区二区| 亚洲真实伦在线观看| 特大巨黑吊av在线直播| 免费黄频网站在线观看国产| 精品一区二区免费观看| 亚洲国产毛片av蜜桃av| 有码 亚洲区| 蜜桃在线观看..| 久久青草综合色| 夫妻性生交免费视频一级片| 高清欧美精品videossex| 欧美精品一区二区大全| 秋霞在线观看毛片| 中文欧美无线码| 99热这里只有精品一区| 亚洲av.av天堂| 日韩电影二区| 国产精品一区二区在线不卡| 久久鲁丝午夜福利片| 国产免费一区二区三区四区乱码| 日本-黄色视频高清免费观看| 国产 精品1| 黄片无遮挡物在线观看| 国产爽快片一区二区三区| 日韩大片免费观看网站| 国产欧美日韩精品一区二区| 亚洲欧洲精品一区二区精品久久久 | 黑人高潮一二区| 80岁老熟妇乱子伦牲交| 18禁在线播放成人免费| 看十八女毛片水多多多| av不卡在线播放| 亚洲一区二区三区欧美精品| 91精品国产九色| 亚洲综合色惰| 国产欧美日韩精品一区二区| 观看美女的网站| av天堂中文字幕网| 曰老女人黄片| 2022亚洲国产成人精品| 国产91av在线免费观看| 日韩三级伦理在线观看| 国产又色又爽无遮挡免| 亚洲欧美日韩卡通动漫| 日本免费在线观看一区| 少妇 在线观看| 观看av在线不卡| 久久久a久久爽久久v久久| 久久久精品94久久精品| 免费大片黄手机在线观看| 亚洲欧美清纯卡通| 久久免费观看电影| 如何舔出高潮| 黑人高潮一二区| 日日摸夜夜添夜夜添av毛片| 男男h啪啪无遮挡| 国产女主播在线喷水免费视频网站| 久久精品国产亚洲网站| av福利片在线| 精品久久久噜噜| 一级毛片久久久久久久久女| 久久女婷五月综合色啪小说| 久久久久精品性色| 亚洲熟女精品中文字幕| 91成人精品电影| 搡老乐熟女国产| 久久久久久久久久久免费av| 成年人午夜在线观看视频| 国产一区二区三区av在线| 91成人精品电影| 国产高清国产精品国产三级| 亚洲av二区三区四区|