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

    多模式離子推力器放電室和柵極設(shè)計及其性能實驗研究*

    2022-10-16 09:23:26李建鵬趙以德靳伍銀張興民李娟王彥龍
    物理學(xué)報 2022年19期
    關(guān)鍵詞:推力器束流柵極

    李建鵬 趙以德 靳伍銀 張興民 李娟 王彥龍

    1) (蘭州理工大學(xué)機(jī)電工程學(xué)院,蘭州 730050)

    2) (蘭州空間技術(shù)物理研究所,真空技術(shù)與物理重點實驗室,蘭州 730000)

    針對小天體探測任務(wù)對離子推力器設(shè)計要求,完成總體方案設(shè)計,開展了四極環(huán)形磁鋼會切場放電室和凸面三柵變孔徑柵極設(shè)計,采用試驗研究和理論分析的方法研究了放電室、柵極設(shè)計的合理性及匹配性,試驗結(jié)果表明: 多模式離子推力器實現(xiàn)了寬范圍穩(wěn)定放電和引出,在277—3120 W 功率下,推力從9.9 mN 線性增加到117.2 mN,比沖從1269 s 臺階上升到3492 s,束流發(fā)散角從30.7°下降到26.8°并趨于穩(wěn)定,各工作點矢量偏角小于1.5°,束流平直度大于0.75,柵極穩(wěn)定化處理是控制柵間距一致性的重要方法,采用應(yīng)力釋放鉬安裝環(huán)熱設(shè)計后柵極間距極差得到有效控制,最大下降百分比為90%,本研究為離子推力器多模式化設(shè)計及在軌應(yīng)用工程研制提供參考.

    1 引言

    離子電推進(jìn)因其高比沖、長壽命、推力寬范圍精確可調(diào)等優(yōu)勢成為先進(jìn)衛(wèi)星平臺的標(biāo)志,為航天器提供動力,在星際深空探測領(lǐng)域,多模式離子推力器更是受到各國的普遍重視并得到工程應(yīng)用[1-3].美國DS-1 和DAWN 號深空探測器配備NSTAR-30 推力器,口徑為30 cm,放電室設(shè)計采用三極環(huán)型會切場,離子光學(xué)系統(tǒng)采用凸面鉬雙柵,推力在20—92 mN 寬范圍可調(diào)[4-6],隨后研發(fā)36 cm NEXT推力器,放電室采用多極環(huán)切場,并采取多余物抑制措施,離子光學(xué)系統(tǒng)為鉬雙柵結(jié)構(gòu),推進(jìn)劑采用放電室反向供給,推力在26—236 mN 寬范圍可調(diào),比沖最高4200 s[7].同時針對深空探測美國也開發(fā)了口徑為57 cm的NEXIS 推力器,放電室采用六極環(huán)切場,離子光學(xué)系統(tǒng)采用碳碳復(fù)合材料,功率覆蓋范圍13—28 kW[8].22 cm T6 推力器采用發(fā)散場磁場構(gòu)型放電室,鋇鎢空心陰極,變孔徑鉬屏柵+碳加速柵凹面雙柵離子光學(xué)系統(tǒng),在2.43—4.50 kW 功率下,推力75—145 mN,比沖最高達(dá)到4120 s[9].

    國內(nèi)外學(xué)者就多模式離子推力器等離子體運動、碰撞、激發(fā)、電離、輸送、約束等開展了大量理論研究[10-14],放電室仿真模擬方面,Brophy 和Wilbur[15]建立了描述放電室性能的零維模型,Arakawa 和Wilbur[16]與Mahalingam 和Menart[17]基于PIC-MCC 算法通過二維模型探究了磁場結(jié)構(gòu)對放電室等離子體的約束.Wang 等[18]、陳茂林等[19]、龍建飛等[20]以及Wirz 和Goebel[21]通過仿真研究了多模式離子推力器離子光學(xué)系統(tǒng)設(shè)計及其對推力器的影響,多模式離子推力器在國外已實現(xiàn)在軌應(yīng)用,公開研究成果聚集于在軌表現(xiàn),地面磨損測試等方面,但對離子推力器多模式化研究、詳細(xì)工程設(shè)計報道較少.

    我國研究起步較晚,蘭州空間技術(shù)物理研究所研制的單工況點LIPS-200、雙工況點LIPS-300 推力器適用于地球同步軌道衛(wèi)星,分別于2017 年、2020 年實現(xiàn)在軌工程應(yīng)用[22-24],針對我國首次小天體探測任務(wù)探測器與太陽距離1—3.5 AU 大范圍變化需求,設(shè)計了300—3000 W 多工作點電推進(jìn)系統(tǒng),為小天體探測器從地球向主帶彗星轉(zhuǎn)移段提供必要的推力[25].

    本文基于我國首次小天體探測任務(wù)對推力器的實際需求,提出總體方案設(shè)計,開展放電室寬范圍適應(yīng)性、離子光學(xué)系統(tǒng)屏柵變孔徑及熱穩(wěn)定性設(shè)計,搭建性能測試、熱穩(wěn)定性和束流分布特性測試平臺,驗證推力器放電室、離子光學(xué)系統(tǒng)設(shè)計的合理性及匹配性,本文的研究對離子推力器多模式化設(shè)計及工程化樣機(jī)研制有重要意義.

    2 多模式離子推力器設(shè)計

    離子推力器多模式設(shè)計要求提供原初電子的主陰極具有寬范圍電子發(fā)射能力,放電室兼顧寬范圍適應(yīng)性和小功率放電穩(wěn)定性,離子光學(xué)系統(tǒng)要解決大功率模式下柵極中心束流欠聚焦和小功率模式下柵極邊緣束流過聚焦問題,最終通過放電室與離子光學(xué)系統(tǒng)匹配設(shè)計實現(xiàn)較高的束流平直度,離子推力器在結(jié)構(gòu)、幾何參數(shù)一定的條件下,通過電氣參數(shù)的控制實現(xiàn)多模式調(diào)節(jié)和穩(wěn)定輸出.

    2.1 總體方案設(shè)計

    推力器總體方案設(shè)計:為了提升放電室放電效率、改善放電室下游束流分布特性,放電室采用四極環(huán)形會切場構(gòu)型,永磁鐵為耐高溫低磁損釤鈷材料,空心陰極采用20 A 全石墨觸持極六硼化鑭發(fā)射體方案,離子光學(xué)系統(tǒng)采用高引出、低濺射率三柵極技術(shù)方案以提高工作壽命,并采用應(yīng)力釋放結(jié)構(gòu)及鉬材料安裝環(huán)增加了熱穩(wěn)定性,供電采用高低壓分組及大電流供電內(nèi)引線三線并聯(lián)技術(shù)方案,總體承力結(jié)構(gòu)采用主支撐環(huán)集成式框架設(shè)計.

    2.2 放電室

    放電室是產(chǎn)生和約束等離子體的場所,小天體探測多模式離子推力器放電室的設(shè)計需要重點考慮功率寬范圍適應(yīng)性和低功率下放電穩(wěn)定性.

    束電流口徑設(shè)計是推力器設(shè)計的初始條件,由Child-Langmuir 方程可得柵極組件的最大束電流密度為[26]

    式中,ε0為真空介電常數(shù);e為元電荷量;mi為工質(zhì)氣體離子質(zhì)量;le為柵極有效柵間距,le=為柵極間距,ts為屏柵厚度,ds為屏柵孔直徑;Vt為屏柵與加速柵之間的電壓.

    考慮推力器可靠性設(shè)計裕度,取極限束電流密度的80%為推力器工作最大電流密度,則束電流口徑:

    推力器最大工況下Ibmax=2.1 A,結(jié)合下文柵極設(shè)計參數(shù),束電流口徑需大于24.6 cm,考慮設(shè)計裕度及柵極抗力學(xué)特性,柵極口徑設(shè)計為30 cm.

    其中放電室口徑Dd取值一般較束電流口徑略大[27],取:

    放電室長度設(shè)計需要兼顧工質(zhì)利用率和電離效率,長徑比Ld/Dd一般在0.4—1 范圍內(nèi)分布,放電室?guī)缀螛?gòu)型采用直段陽極筒+錐段陽極筒的方式,這主要是為了增加原初電子的約束長度,提高推力器效率和束流均勻性,減少對柵極的濺射刻蝕,提升推力器壽命.

    多模式離子推力器整體結(jié)構(gòu)設(shè)計參數(shù)來源于最大功率約束條件下束流口徑,極限束電流密度等.最小功率點下,離子推力器放電室內(nèi)工質(zhì)氣體流率較小,為了兼顧最大功率點設(shè)計要求,放電室?guī)缀螀?shù)較大,小流率下放電室中性氣體濃度低,當(dāng)陰極發(fā)射電流一定,放電室等離子體密度低容易引發(fā)放電不穩(wěn)定.

    因此主陰極發(fā)射電流和工質(zhì)流量大小也是影響推力器小功率下放電穩(wěn)定性的重要因素,通過增加放電室氣體工質(zhì)流量,犧牲推力器工質(zhì)利用率,增大主陰極觸持電流犧牲效率可提高多模式離子推力器小功率下放電穩(wěn)定性.但為了提升推力器性能,拓寬推力器多模式調(diào)節(jié)能力,需要開展磁場的優(yōu)化設(shè)計.

    多模式離子推力器放電室磁場由四組磁環(huán)組成,如圖1 所示,分別為上磁極、中磁極、錐段磁極和下磁極,通過偶數(shù)的四極場可以使靠近柵極的區(qū)域趨向為無場區(qū),使得噴口處等離子體密度均勻性得到改善,滿足功率寬范圍多工作點輸出要求.

    圖1 四極環(huán)形磁鐵會切場示意圖Fig.1.Schematic diagram of ring magnets cusp field.

    環(huán)形會切場放電室將擺放永磁鐵的四周打開,南北極方向交替的永磁體組成的磁場結(jié)構(gòu)如圖2所示,從圖可以看出,永磁鐵會形成尖頭形磁力線結(jié)構(gòu),靠近永磁鐵的磁感應(yīng)強(qiáng)度有一個最大值,并沿著腔內(nèi)方向變小,這種結(jié)構(gòu)使得放電室等離子體空間的大部分為無磁場區(qū).

    圖2 四極會切磁場多極邊界剖面圖,包括磁力線及磁場等值線分布Fig.2.Cross section (side) view of a four-ring-cusp magnetic multipole boundary showing the magnetic field lines and examples of contours of constant magnetic field.

    其中放電室磁場強(qiáng)度的空間分布為[26]By是磁極表面的磁場強(qiáng)度,d是磁極間的距離,w是磁極的寬度,x是沿著陽極壁面距離磁極的距離,y是與陽極壁垂直的方向.

    磁鐵間中心線總的磁場強(qiáng)度為

    q是磁偶極子數(shù)量,δ是磁鐵高度的一半,M是每個磁鐵單位長度的磁偶極子磁場強(qiáng)度,等于磁偶極子的數(shù)量乘以磁鐵長度:

    其中Br是磁鐵剩磁,Vm是磁鐵體積.

    放電室磁場設(shè)計要兼顧放電高效性、均勻性、小功率下穩(wěn)定性,磁場強(qiáng)度高可提升推力器效率,但太高會導(dǎo)致陽極吸收面積較小,陽極面積小到一定值時等離子體電位將變?yōu)樨?fù)值,放電電壓一定,負(fù)等離子體電位降低了原初電子能量,電子、原子的碰撞頻率降低,電離率下降,放電室在小功率下放電震蕩變大,甚至導(dǎo)致推力器熄弧.

    離子推力器放電電流:

    其中kB為玻爾茲曼常數(shù),1.3807×10—23J/K ;Te為電子溫度,eV ;m為電子質(zhì)量,kg;e為電子電量,1.602×10—19C;ne為等離子體電子密度,m—3;Aa為陽極表面能吸收電子和離子的面積大小,m2;φ為陽極的等離子電勢,V.

    Aa與電子和離子的拉莫爾半徑有關(guān):

    其中rh為混合拉莫爾半徑,m;Lc為磁極的總長度,m;re為電子拉莫爾半徑,m;ri為離子拉莫爾半徑,m.

    從(9)式可以看出,隨著磁場強(qiáng)度減小,陽極表面能吸收電子和離子的面積Aa增大,放電電流增大,即放電電流與磁場強(qiáng)度成反比.為了避免小功率下推力器出現(xiàn)放電不穩(wěn)定甚至熄弧現(xiàn)象,通過減弱放電室磁場強(qiáng)度可增加磁空區(qū),進(jìn)而提高放電均勻性,同時由于磁場對原初電子運動約束能力的降低,放電電流增加,通過減弱磁場與放電電流的設(shè)計匹配,可滿足小天體多模式離子推力器低功率下穩(wěn)定放電需求.

    基于離子推力器口徑和放電室?guī)缀螛?gòu)型,尺寸,磁場設(shè)計中磁極數(shù)量和強(qiáng)度的選擇要考慮將磁場控制在陽極表面區(qū)域,增大放電室內(nèi)部的磁空區(qū),迭代優(yōu)化后放電室內(nèi)部磁感應(yīng)強(qiáng)度分布如圖3所示,兼顧寬范圍放電室均勻性和小功率下放電穩(wěn)定性,多模式離子推力器放電室設(shè)計參數(shù)如表1 所列.

    圖3 放電室磁感應(yīng)強(qiáng)度等值線Fig.3.Magnetic field contours in the ion thruster discharge chamber.

    表1 離子推力器放電室設(shè)計參數(shù)Table 1.Design parameters of the ion thruster discharge chamber.

    2.3 離子光學(xué)系統(tǒng)

    離子光學(xué)系統(tǒng)其引出性能取決于本身的設(shè)計特性和與上游放電室等離子體分布的匹配性.離子光學(xué)系統(tǒng)引出束流因空間電荷效應(yīng)限制而遵守Child-Langmuir 方程,一定總加速電壓下,用單孔導(dǎo)流系數(shù)p表征離子光學(xué)系統(tǒng)的引出能力

    ib為單孔束流強(qiáng)度,因此,要使離子光學(xué)系統(tǒng)獲得較大的導(dǎo)流系數(shù)就需要柵極的孔徑大于柵間距,理想的離子光學(xué)系統(tǒng)要求是束流密度高、發(fā)散角小、均勻度高,選取合適的柵極孔徑和柵間距離可以獲得適宜的導(dǎo)流系數(shù),從而獲得較好的離子引出,減小對柵極的轟擊.

    凹面柵極會擠占放電室空間,沿柵極表面想要設(shè)計成無場區(qū)較困難,束流平直度較差,進(jìn)而引起柵極溫度梯度增加,濺射加劇,同時選擇球面凸柵可保證柵極在寬范圍引出熱循環(huán)工況下形變方向的一致性,因此選擇面向小天體探測的多模式離子推力器采用凸面三柵結(jié)構(gòu).

    束流均勻性直接關(guān)系著推力器的壽命和可靠性,為了提升推力器束流均勻性,采用變孔徑屏柵孔徑ds(r) 與等離子體密度ni(r) 和電子溫度Te(r)分布相匹配,即屏柵上不同徑向位置的屏柵孔徑按該位置的ni和Te來確定,單孔束流ib與等離子體密度ni、電子溫度Te、損失因子 exp(?2ts/ds) 和屏柵孔徑ds各個參數(shù)的關(guān)系為[27]

    將屏柵進(jìn)行分區(qū),然后計算出各區(qū)域ni和Te的相對平均值,若將(11)式指數(shù)項展開為級數(shù),可得

    式中,下標(biāo)k 代表k 區(qū)域,下標(biāo)j 代表j 區(qū)域.根據(jù)各區(qū)域內(nèi)單孔引出束流相同,即ib,k=ib,j,確定不同區(qū)域的屏柵孔徑ds(r),參數(shù)如表2 所列,通過不同的孔徑引出相同的束流來改善束流密度均勻性,其示意圖如圖4 所示.

    圖4 變孔徑屏柵示意圖Fig.4.Diagram of the variable aperture screen grid.

    表2 離子光學(xué)系統(tǒng)屏柵變孔徑分區(qū)及歸一化后參數(shù)Table 2.Normalized parameters for screen grids with variable aperture zones.

    推力器熱態(tài)點火時,柵極材料受熱膨脹造成屏柵、加速柵和減速柵之間的間距沿著半徑方向發(fā)生變化,由于各電極間厚度、幾何透明度、能量沉積的差異性,三柵極工作時各電極的變形速率和變形量不同,柵間距分布不一致使得導(dǎo)流系數(shù)發(fā)生改變,柵孔對中性差使得離子運動軌跡發(fā)生偏移,影響束流均勻性和推力器壽命.柵極組件初始設(shè)計采用無應(yīng)力釋放鈦合金支撐環(huán)方案,點火中推力器雙柵間束流閃爍較多,在高功率點火后甚至出現(xiàn)柵片邊緣翹曲現(xiàn)象,如圖5 所示.

    圖5 改進(jìn)前后柵極邊緣變形情況對比 (a) 無應(yīng)力釋放鈦安裝環(huán);(b) 應(yīng)力釋放鉬安裝環(huán)Fig.5.Comparison of grid edge deformation before and after improvement: (a) Ti mounting ring without strain relief;(b) Mo mounting ring with strain relief.

    為了滿足多模式離子推力器柵極組件寬范圍束流引出需求,提高束流均勻性,延長壽命.采取以下措施提高柵極組件熱穩(wěn)定性: 1) 將柵極組件安裝環(huán)設(shè)計在推力器放電室外部,以降低安裝環(huán)處溫度,減小放電室對柵極組件邊緣的熱影響;2) 柵極安裝環(huán)采用鉬材料代替之前的鈦合金材料,使工作時安裝環(huán)與柵片間的熱變形協(xié)調(diào)一致;3) 柵極安裝環(huán)設(shè)計應(yīng)力釋放槽,提高柵極在高溫下的形變一致性,圖5 給出了優(yōu)化前采用無應(yīng)力釋放鈦安裝環(huán)和優(yōu)化后采用應(yīng)力釋放鉬安裝環(huán)柵極邊緣變形示意圖,可以看出優(yōu)化后推力器高溫工況下柵片邊緣無翹曲,熱穩(wěn)定性好.

    基于多模式離子推力器離子光學(xué)系統(tǒng)束流寬范圍引出特性,設(shè)計參數(shù)如表3 所列.

    表3 離子光學(xué)系統(tǒng)設(shè)計參數(shù)Table 3.Design parameters of the ion optical system for ion thrusters.

    空心陰極設(shè)計發(fā)射電流20 A,為了滿足長壽命需要,陰極觸持極為全石墨材料,空心陰極發(fā)射體為六硼化鑭材料.圖6為推力器原理樣機(jī)照片,圖7為推力器現(xiàn)場工作圖.

    圖6 離子推力器原理樣機(jī)Fig.6.Ion thruster prototype model.

    圖7 離子推力器點火照片F(xiàn)ig.7.Discharge of the ion thruster.

    3 試驗系統(tǒng)和方法

    離子推力器性能測試實驗是在TS-7 B 電推進(jìn)真空實驗系統(tǒng)上開展的,設(shè)備參數(shù)參考文獻(xiàn)[3],試驗組成如圖8 所示,采用激光干涉推力測量系統(tǒng)對推力進(jìn)行實測,對比沖進(jìn)行計算.

    圖8 試驗組成圖Fig.8.Schematic diagram of experimental principle.

    束流發(fā)散角和推力偏角是束流特性的重要表征量,其測試方法是: 通過探針采集所在測點的電流密度來確定束流面積和束流中心,根據(jù)束流面積計算束發(fā)散角γ,根據(jù)束流中心計算推力偏心角λ,束流發(fā)散角通過測量離子推力器出口下游1000 mm處束流密度分布得到,束流均勻性分布測量采用單法拉第探針在束流下游通過線掃描的方式進(jìn)行測量.

    束發(fā)散角:

    式中,γ為推力器束流發(fā)散全角,一般指束流覆蓋90%的區(qū)域;ФB為推力器束徑(300 mm);ФT為束發(fā)散角測試直徑;L為法拉第筒到推力器出口的距離.

    推力偏角:

    其中d為法拉第探測面處束流中心到推力器中心軸線的距離;L0為束流發(fā)散角所在直線與推力器中心軸線的交點到法拉第探針探測面的距離.

    性能測試試驗: 在額定供電、供氣參數(shù)下進(jìn)行23 個工作點下的測試試驗[25],采用激光干涉推力測量系統(tǒng)進(jìn)行推力實測和比沖計算.

    柵極熱穩(wěn)定性設(shè)計驗證: 在推力器不引束流的冷態(tài)情況下開展柵極間距測試,測試時,在柵片上設(shè)置有13 個間距測量點,其中1—12 測量區(qū)域為柵極半徑OF(200 mm)處所在直徑15 mm的圓心區(qū)域,如上圖3 所示,表征柵極四周區(qū)域間距的變化,13 測量區(qū)域為圓心O所在直徑15 mm 區(qū)域,表征柵極中心區(qū)域的間距變化.

    4 結(jié)果和討論

    4.1 功率與束電流、推力、比沖

    離子推力器輸入功率為工作時所輸入的所有電功率之和:

    式中,Vb,Vd,Va,Vk,Vn分別為屏柵電壓、放電電壓、加速電壓、陰極和中和器觸持電壓,單位為V;Ib,Id,Ia,Ik,In分別為束電流、放電電流、加速電流、陰極和中和器觸持電流,單位為A.

    離子推力器工作參數(shù)與性能的關(guān)系如下[23]:

    推力:

    其中α為發(fā)散角修正系數(shù);β為雙荷離子修正系數(shù);M是氙離子的質(zhì)量,2.18×10—25kg;g為重力加速度,9.8 m/s2;e是電子的電量,1.6×10—19C;Ub為束電壓;θ為90%束流區(qū)域內(nèi)的發(fā)散角半角.

    采用激光干涉推力測量系統(tǒng)對推力進(jìn)行實測并與(17)式對比修正,α取0.972,β取0.977.比沖按(18)式計算:

    從圖9 和圖10 可以看出,多模式離子推力器設(shè)計合理,實現(xiàn)了功率寬范圍穩(wěn)定放電和束流可靠引出,能夠滿足多工作點變推力需求,束電流隨功率增大,從0.3 A 不斷上升到2.1 A,推力隨功率近線性增加,在277—3120 W 功率下能夠?qū)崿F(xiàn)9.9—117.2 mN 變推力寬范圍工作,推力實測值與計算值最大誤差小于2.7%,比沖隨功率增加開始快速增大,臺階上升,到達(dá)一定值后趨于穩(wěn)定.

    圖9 離子束電流隨功率變化曲線Fig.9.Ion beam current as a function of input power.

    圖10 推力、比沖隨功率變化曲線Fig.10.Thrust,specific impulse as a function of input power.

    功率寬范圍調(diào)節(jié)策略分為2 步,首先是保持屏柵電壓一定,精確調(diào)節(jié)陽極電流和放電室流率,在其屏柵電壓下引出不同束流,實現(xiàn)功率小范圍調(diào)節(jié),再通過大步長調(diào)節(jié)屏柵電壓實現(xiàn)功率的寬范圍調(diào)節(jié),離子推力器工質(zhì)流量供給分為3 路,即中和器、主陰極和陽極,一定柵極電壓下,離子推力器的比沖取決于流率,推力器寬范圍調(diào)節(jié)流率設(shè)置原則是: 工作點一定,主陰極、中和器工作在點狀模式且比沖最大化,即中和器、陰極工作在點狀模式的流率最小值為該工作點流率設(shè)計值,陽極流率設(shè)置要求是放電室工質(zhì)利用率不高于放電損耗突增轉(zhuǎn)折點,且最高不超過90%,通過電氣參數(shù)大小步長迭代調(diào)節(jié)能夠?qū)崿F(xiàn)比沖在1269—3492 s 寬范圍工作.

    4.2 柵極熱穩(wěn)定性設(shè)計驗證

    柵間距和柵間距分布一致性是離子推力器多模式調(diào)節(jié)和長壽命實現(xiàn)的重要影響因素,這是由于柵極幾何參數(shù)不變,一定束電流下,受離子光學(xué)系統(tǒng)交叉限制和導(dǎo)流系數(shù)限制,柵極可正常聚焦的總加速電壓范圍一定,即柵極束流正常聚焦引出的限制其實為束電流和總加速電壓之間的限制關(guān)系[26].如果柵間距不一致,則交叉限制屏柵電壓受限于柵極的最小熱態(tài)柵間距,同時導(dǎo)流系數(shù)限制屏柵電壓又受限于柵極的最大熱態(tài)柵間距,柵極可正常聚焦多模式調(diào)節(jié)范圍變小,當(dāng)工作在最優(yōu)導(dǎo)流系數(shù)附近,柵極通過靜電場的作用對離子聚焦、定向和加速,能夠避免柵極對離子的攔截.當(dāng)工作電壓或者是電流超過極限范圍后,由于柵極間的空間電荷飽和效應(yīng)會造成離子發(fā)散角急劇增加,并導(dǎo)致高能離子對加速柵和減速柵直接轟擊,降低推力器工作壽命.

    離子推力器柵極為薄壁多孔弱剛度結(jié)構(gòu),由于受到等離子轟擊存在熱變形,在工程應(yīng)用中為了保證柵極組件在多次熱沖擊下的穩(wěn)定性,通常在柵極初始裝配后開展柵極穩(wěn)定化熱處理,方法是將初裝的柵極組件裝配與離子推力器上按照低功率-高功率,高功率-低功率順序開展束流引出,通過高低溫循環(huán)釋放柵極加工、制造、裝配產(chǎn)生的應(yīng)力,因此在熱穩(wěn)定性驗證實驗中分別測量初始裝配和熱穩(wěn)定化后柵極間距用以評估熱設(shè)計.

    圖11 給出了柵極組件采用無應(yīng)力釋放鈦安裝環(huán)和應(yīng)力釋放鉬安裝環(huán)屏柵-加速柵各測量點間距測量值.優(yōu)化前,柵極組件初始裝配后屏柵-加速柵二者間距為0.912 mm,穩(wěn)定化后釋放了應(yīng)力,二者間距增大到1.140 mm,初始裝配與穩(wěn)定化后間距對比,四周變化百分比25%,中心變化百分比25.8%.優(yōu)化后,柵極組件初始裝配后屏柵-加速柵二者間距為0.892 mm,穩(wěn)定化后變化為1.023 mm,四周前后變化百分比14.6%,中心變化百分比18.5%.

    圖11 不同測量點屏柵-加速柵間距 (a) 優(yōu)化前;(b) 優(yōu)化后Fig.11.Screen grid-acceleration grid spacing at different measurement position: (a) Before optimization;(b) after optimization.

    圖12 給出了優(yōu)化前后加速-減速柵各測量點間距測量值,采用無應(yīng)力釋放鈦安裝環(huán),柵極初始裝配和穩(wěn)定化處理后加速-減速柵二者間距分別為0.798 mm 和1.004 mm,采用應(yīng)力釋放鉬安裝環(huán)柵極初始裝配后加速-減速柵二者間距為0.804 mm,穩(wěn)定化后變?yōu)?.825 mm,四周前后變化百分比從優(yōu)化前的25.8%降低到2.6%,中心前后變化百分比從優(yōu)化前的20.8%降低到7.4%.表4 給出了優(yōu)化前后柵極穩(wěn)定化處理后柵間距變化對比數(shù)據(jù).從表4 可以看出,采用應(yīng)力釋放結(jié)構(gòu)優(yōu)化和鉬安裝環(huán)材料優(yōu)化熱設(shè)計提升了柵極熱穩(wěn)定性,屏柵-加速四周極差從0.1 mm 降低到0.01 mm,中心極差從0.330 mm 降低到0.180 mm,加速-減速四周極差從0.290 mm 減低到0.075 mm,中心極差從0.290 mm降低到0.115 mm,中心極差較四周大,主要是因為離子推力器柵片連接采用四周固定,在柵極組件經(jīng)歷高溫工況后遠(yuǎn)離固定位置的中心發(fā)生較大的偏移.

    表4 優(yōu)化前后柵極穩(wěn)定化處理后柵間距變化量Table 4.Variation in grid gap with stabilization before and after optimization.

    圖12 不同測量點加速柵-減速柵間距 (a) 優(yōu)化前;(b) 優(yōu)化后Fig.12.Accelerator-decelerator spacing at different measurement position: (a) Before optimization;(b) after optimization.

    4.3 束流分布特性

    圖13 給出了推力器不同工作點下的束流發(fā)散角、推力偏角測試結(jié)果,可以看出最小功率工作點277 W 時,束流發(fā)散角為30.7°,其余工作點束流發(fā)散角均小于30°,矢量偏角均在1.5°以內(nèi),符合設(shè)計預(yù)期,有利于長壽命實現(xiàn),且隨著功率的增加,束流發(fā)散角呈現(xiàn)變小趨勢并趨于穩(wěn)定,分析認(rèn)為該現(xiàn)象與柵極組件引束流的聚焦性有關(guān),離子鞘層形狀相對柵極位置并非固定不變,而是受等離子體密度,柵間電壓,柵極材料特性、柵孔孔徑和間距等各種因素的影響,理想工作條件下,鞘層彎曲成球冠狀,凹面對著屏柵孔,使得鞘面發(fā)射的離子全部進(jìn)入屏柵孔內(nèi),離子有理想聚焦的運動軌跡,小功率工作點下等離子體密度低,鞘層與離子加速區(qū)域較遠(yuǎn),部分離子束過聚焦導(dǎo)致束流發(fā)散角大[23].此外,高功率下束流發(fā)散角較小也與導(dǎo)流系數(shù)p有關(guān),隨著功率的增大,等離子體密度增大,沉積在柵極組件的能量變大,柵極有效間距l(xiāng)e 減小促使柵極導(dǎo)流能力增強(qiáng),束流發(fā)散角變小,后續(xù)采用邊緣變孔徑設(shè)計,以解決小功率下邊緣區(qū)域束流欠聚焦問題.

    圖13 推力器不同工作點下的束流發(fā)散角和矢量偏角Fig.13.Beam divergence angle and thrust vector angle at different operating points.

    離子束流的均勻性是離子推力器設(shè)計的關(guān)鍵指標(biāo),是影響離子推力器寬范圍調(diào)節(jié)和長壽命實現(xiàn)的關(guān)鍵因素.推力器束流密度徑向分布測試是準(zhǔn)確評估束流平直度的有效方法,束流密度徑向分布曲線是束流均勻性的表征量,可為推力器性能改善和壽命優(yōu)化提供數(shù)據(jù),圖14 給出了推力器功率最大工況點1.2 倍拉偏下柵極熱穩(wěn)定性設(shè)計優(yōu)化前后束流密度對比曲線,可以看出,改進(jìn)前離子推力器束流密度分布呈雙山峰形,改進(jìn)后密度分布曲線基本為平頂形,束流均勻性較好,束流平直度提升了10%.圖15 給出了優(yōu)化后7 個典型工作點下束流密度徑向分布特性測試結(jié)果,可以看出多模式離子推力器隨著功率增大,束電流密度增大,在不同的工作點下束流密度分布曲線無尖峰,較平緩,束流均勻性較好,束流平直度達(dá)到0.75,說明多模式離子推力器放電室磁場和離子光學(xué)系統(tǒng)設(shè)計合理且匹配良好,能夠滿足電推進(jìn)系統(tǒng)多模式,長壽命設(shè)計要求.

    圖14 柵極熱穩(wěn)定設(shè)計優(yōu)化前后束流密度分布Fig.14.Radial beam current density profile before and after optimization of grid thermal stabilization design.

    圖15 推力器不同工作點下的束流密度分布Fig.15.Radial beam current density profile at different operating points.

    5 結(jié)論

    1) 基于小天體探測任務(wù)對電推進(jìn)系統(tǒng)的設(shè)計需求,提出四極環(huán)形會切場放電室,凸面三柵變孔徑柵極,全石墨觸持極六硼化鑭發(fā)射體空心陰極、集成式主支撐環(huán)承力結(jié)構(gòu)的總體設(shè)計方案,開展了功率寬范圍適應(yīng)性、低功率放電穩(wěn)定性的放電室參數(shù)詳細(xì)設(shè)計;

    2) 通過變孔徑屏柵設(shè)計方案實現(xiàn)了引出束流與放電室電子溫度分布相匹配,束流均勻性好,各工況下束流平直度大于0.75,提出應(yīng)力釋放鉬安裝環(huán)熱設(shè)計并開展試驗驗證,結(jié)果表明: 柵極穩(wěn)定化處理是控制柵間距一致性的重要方法,優(yōu)化后柵極變形得到有效控制,柵極間距極差最大下降百分比為90%,優(yōu)化前后屏柵-加速,加速-減速四周極差分別從0.10 mm,0.290 mm 下降到0.01 mm 和0.075 mm,中心極差分別從0.330 mm,0.290 mm下降到0.180 mm 和0.115 mm;

    3) 搭建了離子推力器性能測試和束流特性測試平臺,試驗結(jié)果表明: 多模式離子推力器推力隨功率近線性增加,比沖隨功率臺階上升,束流發(fā)散角隨功率呈現(xiàn)變小趨勢并趨于穩(wěn)定,在277—3120 W功率下能夠?qū)崿F(xiàn)推力在9.9—117.2 mN,比沖在1269—3492 s 寬范圍穩(wěn)定工作,各工作點束流發(fā)散角小于30°,矢量偏角小于1.5°.

    猜你喜歡
    推力器束流柵極
    單組元推力器倒置安裝多余物控制技術(shù)驗證
    離子推力器三柵極組件熱形變仿真分析及試驗研究
    真空與低溫(2022年2期)2022-03-30 07:11:22
    高能同步輻射光源低能束流輸運線設(shè)計研究
    中國散裂中子源加速器注入束流損失調(diào)節(jié)研究
    柵極液壓成型專用設(shè)備的研制
    ECR中和器束流引出實驗研究
    用于小行星探測的離子推力器技術(shù)研究
    IGBT柵極驅(qū)動電阻的選擇
    離子推力器和霍爾推力器的異同
    太空探索(2015年6期)2015-07-12 12:48:42
    一種無升壓結(jié)構(gòu)的MOSFET柵極驅(qū)動電路
    大陆偷拍与自拍| 亚洲avbb在线观看| 亚洲精品中文字幕一二三四区 | 中文字幕制服av| 九色亚洲精品在线播放| www日本在线高清视频| 中文字幕av电影在线播放| 制服人妻中文乱码| 国产精品一区二区在线不卡| 高清在线国产一区| 高清av免费在线| 男女国产视频网站| 法律面前人人平等表现在哪些方面 | 亚洲精华国产精华精| 黄频高清免费视频| 亚洲国产日韩一区二区| 国产91精品成人一区二区三区 | 热99久久久久精品小说推荐| 国产精品国产av在线观看| 男女无遮挡免费网站观看| 国产日韩欧美在线精品| 国产野战对白在线观看| 丝袜在线中文字幕| 超碰97精品在线观看| 国产日韩欧美视频二区| 丝袜脚勾引网站| 又黄又粗又硬又大视频| 免费日韩欧美在线观看| 高清欧美精品videossex| 一边摸一边抽搐一进一出视频| 欧美 亚洲 国产 日韩一| 五月天丁香电影| 老司机亚洲免费影院| 五月天丁香电影| 精品福利观看| 亚洲国产精品成人久久小说| 狂野欧美激情性bbbbbb| 国产精品1区2区在线观看. | 另类亚洲欧美激情| 欧美午夜高清在线| 国产xxxxx性猛交| 精品亚洲成a人片在线观看| 午夜影院在线不卡| 一个人免费看片子| 精品亚洲成国产av| 精品亚洲成国产av| 国产日韩欧美在线精品| 男女无遮挡免费网站观看| avwww免费| 亚洲精品久久成人aⅴ小说| 日韩精品免费视频一区二区三区| 久热这里只有精品99| 蜜桃在线观看..| 亚洲精品久久成人aⅴ小说| 亚洲国产精品999| 91麻豆精品激情在线观看国产 | 中文字幕人妻丝袜制服| 国产亚洲精品久久久久5区| 国产在线视频一区二区| 国产精品香港三级国产av潘金莲| 999久久久国产精品视频| 亚洲精品在线美女| 97人妻天天添夜夜摸| 国产精品久久久久久精品电影小说| 国产免费视频播放在线视频| 脱女人内裤的视频| 日日爽夜夜爽网站| 中文字幕高清在线视频| 亚洲中文日韩欧美视频| 天堂8中文在线网| 最黄视频免费看| 亚洲熟女精品中文字幕| 香蕉国产在线看| 久久九九热精品免费| 日本av免费视频播放| 波多野结衣一区麻豆| 中国美女看黄片| a在线观看视频网站| 日本av免费视频播放| 亚洲人成电影免费在线| 两性夫妻黄色片| 最黄视频免费看| 日韩一卡2卡3卡4卡2021年| 一本—道久久a久久精品蜜桃钙片| 女警被强在线播放| 女警被强在线播放| 国产成人av教育| 成人国语在线视频| 成年女人毛片免费观看观看9 | 美女主播在线视频| 日本av免费视频播放| 欧美精品一区二区免费开放| 亚洲天堂av无毛| 日韩制服骚丝袜av| 一本—道久久a久久精品蜜桃钙片| 久久精品久久久久久噜噜老黄| 色婷婷av一区二区三区视频| 国产精品亚洲av一区麻豆| 在线观看免费午夜福利视频| 在线观看免费午夜福利视频| 国产一区二区在线观看av| 亚洲综合色网址| 精品第一国产精品| 国产1区2区3区精品| 中文字幕高清在线视频| 王馨瑶露胸无遮挡在线观看| 亚洲精品国产av成人精品| 老熟妇乱子伦视频在线观看 | 国产精品香港三级国产av潘金莲| 在线观看一区二区三区激情| 热99久久久久精品小说推荐| 欧美黑人欧美精品刺激| 亚洲国产精品999| 久久人妻熟女aⅴ| 中文字幕另类日韩欧美亚洲嫩草| 最近中文字幕2019免费版| 欧美日本中文国产一区发布| 黄频高清免费视频| 久久精品国产亚洲av高清一级| 色精品久久人妻99蜜桃| 制服诱惑二区| 别揉我奶头~嗯~啊~动态视频 | 国产精品久久久久久精品古装| 亚洲精品国产区一区二| 午夜精品国产一区二区电影| 国产精品一二三区在线看| 国产高清视频在线播放一区 | 国产精品成人在线| 欧美日韩视频精品一区| 精品国产乱码久久久久久小说| 操美女的视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人欧美在线观看 | 国产欧美日韩一区二区三 | 欧美xxⅹ黑人| 满18在线观看网站| 天天操日日干夜夜撸| 最近最新免费中文字幕在线| 久久99一区二区三区| 欧美精品一区二区大全| 如日韩欧美国产精品一区二区三区| 三级毛片av免费| 菩萨蛮人人尽说江南好唐韦庄| 国产成人欧美在线观看 | videosex国产| 日韩熟女老妇一区二区性免费视频| 啦啦啦 在线观看视频| 国内毛片毛片毛片毛片毛片| 午夜福利在线免费观看网站| 永久免费av网站大全| 这个男人来自地球电影免费观看| 欧美日韩福利视频一区二区| 国产男人的电影天堂91| 久久久久久久久免费视频了| 99国产精品99久久久久| 色婷婷久久久亚洲欧美| 手机成人av网站| 妹子高潮喷水视频| 亚洲伊人色综图| 巨乳人妻的诱惑在线观看| 一二三四社区在线视频社区8| 国产亚洲精品久久久久5区| 亚洲黑人精品在线| 成年人黄色毛片网站| 日本vs欧美在线观看视频| 国产精品 欧美亚洲| 亚洲中文日韩欧美视频| 高清欧美精品videossex| 欧美日韩黄片免| 久久久精品免费免费高清| www.av在线官网国产| 精品国产乱码久久久久久男人| 国产成人影院久久av| 制服诱惑二区| 婷婷丁香在线五月| 窝窝影院91人妻| 99国产精品一区二区三区| 999精品在线视频| 午夜免费成人在线视频| av视频免费观看在线观看| 亚洲男人天堂网一区| 欧美黄色片欧美黄色片| 99久久综合免费| 十分钟在线观看高清视频www| 色播在线永久视频| 伊人久久大香线蕉亚洲五| 天天添夜夜摸| 十分钟在线观看高清视频www| www日本在线高清视频| 国产真人三级小视频在线观看| 亚洲国产欧美一区二区综合| 一级毛片电影观看| 亚洲欧美激情在线| 亚洲久久久国产精品| 欧美日韩一级在线毛片| 黄色毛片三级朝国网站| 久久精品国产a三级三级三级| 黑丝袜美女国产一区| 久久久久久亚洲精品国产蜜桃av| 亚洲自偷自拍图片 自拍| 男女午夜视频在线观看| 黄频高清免费视频| 夫妻午夜视频| 热re99久久精品国产66热6| www.999成人在线观看| 一区二区av电影网| 欧美人与性动交α欧美软件| 男女无遮挡免费网站观看| 日韩大码丰满熟妇| 淫妇啪啪啪对白视频 | 天天添夜夜摸| av网站在线播放免费| 久久人人爽人人片av| 久久天堂一区二区三区四区| 午夜精品久久久久久毛片777| 欧美日韩国产mv在线观看视频| 在线观看人妻少妇| 亚洲黑人精品在线| 人人澡人人妻人| www.自偷自拍.com| 国产免费福利视频在线观看| www.精华液| 精品少妇一区二区三区视频日本电影| 日本91视频免费播放| 啪啪无遮挡十八禁网站| 久久综合国产亚洲精品| 亚洲国产欧美日韩在线播放| 久久女婷五月综合色啪小说| 91麻豆av在线| 日韩中文字幕欧美一区二区| 美女高潮到喷水免费观看| 一本—道久久a久久精品蜜桃钙片| 老司机在亚洲福利影院| 男女床上黄色一级片免费看| 久9热在线精品视频| 丰满少妇做爰视频| 欧美精品高潮呻吟av久久| 亚洲国产欧美一区二区综合| 91大片在线观看| 亚洲 国产 在线| 亚洲一码二码三码区别大吗| 亚洲国产精品一区三区| 18禁观看日本| 高清黄色对白视频在线免费看| 丰满迷人的少妇在线观看| 欧美日韩黄片免| 欧美日韩亚洲综合一区二区三区_| 国产成人啪精品午夜网站| 建设人人有责人人尽责人人享有的| 国产免费现黄频在线看| 久久久久国产一级毛片高清牌| 下体分泌物呈黄色| 国产伦理片在线播放av一区| 国产一区二区在线观看av| 曰老女人黄片| 在线精品无人区一区二区三| 精品国产一区二区久久| 9热在线视频观看99| 免费在线观看完整版高清| 午夜老司机福利片| 老司机亚洲免费影院| 亚洲成人国产一区在线观看| 午夜老司机福利片| 国产精品久久久久久精品电影小说| 亚洲成av片中文字幕在线观看| 老司机午夜十八禁免费视频| 老司机亚洲免费影院| 欧美在线黄色| 超色免费av| 国产成人啪精品午夜网站| xxxhd国产人妻xxx| 日韩 亚洲 欧美在线| 精品亚洲乱码少妇综合久久| 精品乱码久久久久久99久播| 久久中文字幕一级| 国产成人影院久久av| 免费不卡黄色视频| 国产片内射在线| 在线观看免费高清a一片| 国产精品熟女久久久久浪| 老汉色∧v一级毛片| 亚洲国产看品久久| 女性生殖器流出的白浆| 午夜福利在线观看吧| www日本在线高清视频| 男人爽女人下面视频在线观看| 亚洲精品国产区一区二| 国产在线观看jvid| 国产精品九九99| 亚洲一区二区三区欧美精品| 桃花免费在线播放| xxxhd国产人妻xxx| 久热爱精品视频在线9| 亚洲国产成人一精品久久久| 久久亚洲精品不卡| 黑丝袜美女国产一区| 国产一区二区三区av在线| 午夜成年电影在线免费观看| 国产成人精品久久二区二区免费| 交换朋友夫妻互换小说| 各种免费的搞黄视频| 亚洲av成人一区二区三| 激情视频va一区二区三区| 少妇被粗大的猛进出69影院| 老司机影院毛片| 男人操女人黄网站| 丰满迷人的少妇在线观看| 水蜜桃什么品种好| 午夜福利视频在线观看免费| 在线天堂中文资源库| 十八禁人妻一区二区| 国产男人的电影天堂91| 99久久99久久久精品蜜桃| 丰满饥渴人妻一区二区三| h视频一区二区三区| 美女脱内裤让男人舔精品视频| 欧美日韩视频精品一区| 一级片'在线观看视频| 咕卡用的链子| 日韩欧美国产一区二区入口| 精品人妻熟女毛片av久久网站| 汤姆久久久久久久影院中文字幕| 亚洲欧洲日产国产| av电影中文网址| 欧美+亚洲+日韩+国产| 美女高潮喷水抽搐中文字幕| 一级,二级,三级黄色视频| 国产男人的电影天堂91| 女人爽到高潮嗷嗷叫在线视频| 狠狠精品人妻久久久久久综合| 亚洲成av片中文字幕在线观看| 这个男人来自地球电影免费观看| 大片电影免费在线观看免费| 9191精品国产免费久久| av线在线观看网站| 久久亚洲国产成人精品v| 免费看十八禁软件| 一级片'在线观看视频| 男人爽女人下面视频在线观看| 欧美av亚洲av综合av国产av| 美女脱内裤让男人舔精品视频| 精品久久久久久电影网| 精品国产乱子伦一区二区三区 | 欧美精品高潮呻吟av久久| 午夜福利免费观看在线| 一区在线观看完整版| 亚洲情色 制服丝袜| 国产亚洲av高清不卡| 搡老岳熟女国产| 大陆偷拍与自拍| 女警被强在线播放| 高潮久久久久久久久久久不卡| 亚洲精品成人av观看孕妇| av视频免费观看在线观看| 在线看a的网站| 亚洲国产av新网站| 窝窝影院91人妻| 飞空精品影院首页| 99re6热这里在线精品视频| 又大又爽又粗| 伦理电影免费视频| 欧美精品亚洲一区二区| 久久国产精品男人的天堂亚洲| 91av网站免费观看| 欧美黄色淫秽网站| 真人做人爱边吃奶动态| 老熟妇仑乱视频hdxx| 国产欧美亚洲国产| 欧美人与性动交α欧美精品济南到| 青青草视频在线视频观看| 黄网站色视频无遮挡免费观看| 午夜福利一区二区在线看| 永久免费av网站大全| 18禁裸乳无遮挡动漫免费视频| av不卡在线播放| 久久国产精品男人的天堂亚洲| 久久人妻熟女aⅴ| 精品欧美一区二区三区在线| 老司机影院成人| 巨乳人妻的诱惑在线观看| 精品福利永久在线观看| 欧美另类亚洲清纯唯美| 欧美人与性动交α欧美软件| 无遮挡黄片免费观看| 免费女性裸体啪啪无遮挡网站| 精品欧美一区二区三区在线| 亚洲va日本ⅴa欧美va伊人久久 | 日本黄色日本黄色录像| 最黄视频免费看| 亚洲国产欧美网| 精品国产一区二区三区四区第35| 男女免费视频国产| 国产97色在线日韩免费| 夜夜骑夜夜射夜夜干| 日本欧美视频一区| av超薄肉色丝袜交足视频| 丰满迷人的少妇在线观看| 免费在线观看影片大全网站| 国产精品.久久久| 热99re8久久精品国产| 国产亚洲精品一区二区www | av福利片在线| 亚洲国产欧美网| 欧美av亚洲av综合av国产av| 高清黄色对白视频在线免费看| 午夜成年电影在线免费观看| 亚洲精品成人av观看孕妇| 中文字幕人妻丝袜一区二区| 久久久国产成人免费| 黑人操中国人逼视频| 青春草视频在线免费观看| 久久精品成人免费网站| 亚洲色图综合在线观看| 国产野战对白在线观看| 国产福利在线免费观看视频| 高清在线国产一区| 中文字幕制服av| 99精国产麻豆久久婷婷| 免费久久久久久久精品成人欧美视频| 亚洲欧美日韩另类电影网站| 黄色视频不卡| 成人国语在线视频| 国产精品免费视频内射| 日韩,欧美,国产一区二区三区| e午夜精品久久久久久久| 岛国在线观看网站| 亚洲男人天堂网一区| 一进一出抽搐动态| 在线精品无人区一区二区三| 99精品久久久久人妻精品| 一边摸一边抽搐一进一出视频| 80岁老熟妇乱子伦牲交| 大香蕉久久成人网| 久久 成人 亚洲| 欧美在线黄色| 制服诱惑二区| 一区二区三区乱码不卡18| av欧美777| 国产av又大| 一本久久精品| a级片在线免费高清观看视频| 99香蕉大伊视频| 久久精品亚洲熟妇少妇任你| 亚洲男人天堂网一区| 精品久久久久久电影网| 狂野欧美激情性bbbbbb| 国产三级黄色录像| 夜夜夜夜夜久久久久| 少妇 在线观看| 国产精品一区二区在线观看99| 熟女少妇亚洲综合色aaa.| 欧美 日韩 精品 国产| 性色av乱码一区二区三区2| 午夜老司机福利片| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品成人av观看孕妇| 国产亚洲欧美在线一区二区| 国产亚洲欧美精品永久| 亚洲精品国产区一区二| 成年美女黄网站色视频大全免费| 欧美大码av| tube8黄色片| 欧美亚洲 丝袜 人妻 在线| 99九九在线精品视频| 久久中文看片网| 国产精品亚洲av一区麻豆| 不卡一级毛片| 国产免费福利视频在线观看| 高潮久久久久久久久久久不卡| 国产真人三级小视频在线观看| 午夜精品久久久久久毛片777| 老熟妇乱子伦视频在线观看 | 国产黄色免费在线视频| 亚洲国产欧美网| 黄片小视频在线播放| 亚洲免费av在线视频| 国产精品久久久久久精品古装| 90打野战视频偷拍视频| 欧美另类亚洲清纯唯美| 天天躁夜夜躁狠狠躁躁| 纵有疾风起免费观看全集完整版| 国产免费视频播放在线视频| 午夜免费观看性视频| 少妇人妻久久综合中文| 亚洲熟女毛片儿| 最近最新中文字幕大全免费视频| 50天的宝宝边吃奶边哭怎么回事| 美女中出高潮动态图| 日韩有码中文字幕| a级毛片黄视频| 久久女婷五月综合色啪小说| 高清在线国产一区| 新久久久久国产一级毛片| 免费人妻精品一区二区三区视频| 国产精品久久久久成人av| 捣出白浆h1v1| 丰满人妻熟妇乱又伦精品不卡| 人人妻,人人澡人人爽秒播| 久久九九热精品免费| 丰满少妇做爰视频| 亚洲av美国av| 国产亚洲午夜精品一区二区久久| 12—13女人毛片做爰片一| 一级,二级,三级黄色视频| xxxhd国产人妻xxx| 黄色 视频免费看| 精品福利永久在线观看| 男男h啪啪无遮挡| 另类亚洲欧美激情| 亚洲免费av在线视频| 亚洲欧美一区二区三区久久| 精品一品国产午夜福利视频| 国产精品秋霞免费鲁丝片| 高潮久久久久久久久久久不卡| 午夜激情av网站| 日韩视频一区二区在线观看| 亚洲精品国产色婷婷电影| 日本五十路高清| 操出白浆在线播放| 亚洲熟女毛片儿| av网站免费在线观看视频| 王馨瑶露胸无遮挡在线观看| 欧美激情极品国产一区二区三区| 国产片内射在线| 精品国产一区二区三区四区第35| 日本vs欧美在线观看视频| 黄色 视频免费看| 一区二区三区精品91| 老熟妇乱子伦视频在线观看 | 亚洲国产欧美一区二区综合| 黑人操中国人逼视频| 精品人妻熟女毛片av久久网站| 欧美国产精品一级二级三级| 国产野战对白在线观看| 后天国语完整版免费观看| 亚洲精品成人av观看孕妇| 免费观看av网站的网址| 国产免费视频播放在线视频| 91精品三级在线观看| 午夜久久久在线观看| 亚洲国产毛片av蜜桃av| 国产麻豆69| 亚洲专区中文字幕在线| 别揉我奶头~嗯~啊~动态视频 | 天堂中文最新版在线下载| 午夜91福利影院| 国产精品.久久久| 脱女人内裤的视频| 99精品欧美一区二区三区四区| 免费久久久久久久精品成人欧美视频| 日本91视频免费播放| av片东京热男人的天堂| 亚洲国产成人一精品久久久| 久久精品亚洲av国产电影网| av免费在线观看网站| 亚洲七黄色美女视频| 夜夜骑夜夜射夜夜干| 亚洲人成77777在线视频| 亚洲国产欧美网| 人人妻人人澡人人爽人人夜夜| 啦啦啦在线免费观看视频4| 老汉色av国产亚洲站长工具| 久久人人爽av亚洲精品天堂| 国产一卡二卡三卡精品| 欧美日韩亚洲国产一区二区在线观看 | 国产精品 国内视频| 搡老熟女国产l中国老女人| 性高湖久久久久久久久免费观看| 亚洲精品久久午夜乱码| 宅男免费午夜| 国产成人啪精品午夜网站| 国产日韩欧美在线精品| 久久久精品国产亚洲av高清涩受| 老司机在亚洲福利影院| 在线观看一区二区三区激情| 大型av网站在线播放| √禁漫天堂资源中文www| 婷婷成人精品国产| 亚洲九九香蕉| 中文字幕最新亚洲高清| 国产精品影院久久| 汤姆久久久久久久影院中文字幕| 9191精品国产免费久久| 亚洲精品粉嫩美女一区| 99久久综合免费| 亚洲精品美女久久av网站| 亚洲 欧美一区二区三区| 精品人妻熟女毛片av久久网站| 黄片播放在线免费| 亚洲男人天堂网一区| 91老司机精品| 欧美黑人欧美精品刺激| 久久久久久久久免费视频了| 99精品欧美一区二区三区四区| 久久人人爽av亚洲精品天堂| 深夜精品福利| 久久久久久久精品精品| 亚洲精品美女久久av网站| 日日爽夜夜爽网站| 性高湖久久久久久久久免费观看| 91精品国产国语对白视频| 色94色欧美一区二区| 国产免费一区二区三区四区乱码| 日韩视频一区二区在线观看| 精品久久久精品久久久| 国产有黄有色有爽视频| 色婷婷久久久亚洲欧美| 黄色视频,在线免费观看| 王馨瑶露胸无遮挡在线观看| 一本综合久久免费| 中文字幕人妻丝袜制服| av不卡在线播放| 电影成人av| 黑人巨大精品欧美一区二区mp4| 亚洲国产精品一区二区三区在线|