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

    地球靜止軌道LIPS-300離子電推力器熱設(shè)計與優(yōu)化

    2016-02-13 07:43:52胡幗杰李健劉百麟
    中國空間科學技術(shù) 2016年1期
    關(guān)鍵詞:磁鋼推力器外殼

    胡幗杰,李健,劉百麟

    中國空間技術(shù)研究院通信衛(wèi)星事業(yè)部,北京100094

    地球靜止軌道LIPS-300離子電推力器熱設(shè)計與優(yōu)化

    胡幗杰*,李健,劉百麟

    中國空間技術(shù)研究院通信衛(wèi)星事業(yè)部,北京100094

    LIPS-300離子電推力器工作時產(chǎn)生的大量熱耗會嚴重影響其工作性能,因此推力器的熱控設(shè)計和在軌溫度預(yù)計是推力器設(shè)計的關(guān)鍵技術(shù)之一。文章首先提出了地球靜止軌道LIPS-300離子電推力器熱設(shè)計,建立了離子電推力器熱分析有限元模型。然后,針對典型的地球靜止軌道位保工況開展了LIPS-300離子電推力器的在軌熱分析,得到了LIPS-300離子電推力器重點部件在軌工作時的溫度預(yù)計:下極靴、陽極筒、屏柵、下磁鋼為主要的高溫位置,且下磁鋼溫度可能超出材料溫度上限。最后,基于熱分析的結(jié)果結(jié)合離子電推力器的特點提出適當增加外殼散熱面積是離子電推力器優(yōu)化設(shè)計的有效方向,而增強陽極筒、極靴與屏柵筒間的內(nèi)部換熱是無效的設(shè)計優(yōu)化方向。

    地球靜止軌道;離子電推力器;熱設(shè)計;在軌熱分析;優(yōu)化設(shè)計

    Key words:geosynchronous orbit;ion thruster;thermal design;thermal analysis in orbit; optimization design

    離子電推進技術(shù)是先進的空間推進技術(shù),其比沖比常規(guī)化學燃料發(fā)動機高出數(shù)倍,大大提高了衛(wèi)星的有效載荷質(zhì)量,延長了衛(wèi)星在軌工作壽命。隨著航天器在軌應(yīng)用需求的不斷增加,對離子電推進系統(tǒng)也有了更高的要求。大推力、高比沖、長壽命及高可靠性已經(jīng)成為目前離子電推力器技術(shù)發(fā)展的趨勢[1-2]。

    中國蘭州空間技術(shù)物理研究所研制的20 cm口徑LIPS-200離子電推力器已經(jīng)于2012年進行了飛行演示試驗[3],得到了成功驗證。為了獲得更大的推力和比沖,在LIPS-200離子電推力器的基礎(chǔ)上又研制了推力更大的30 cm口徑LIPS-300離子電推力器[4]。不過,LIPS-300離子電推力器在推力增加的同時產(chǎn)生的熱量也大大增加,從而嚴重影響電推力器的工作性能和使用壽命[5-7],因此電推力器的熱控設(shè)計和在軌工作溫度預(yù)計是電推力器設(shè)計的關(guān)鍵技術(shù)之一。

    關(guān)于電推力器的設(shè)計技術(shù),公開發(fā)表的可借鑒資料甚少,尤其缺少詳細的設(shè)計技術(shù)資料。從可見的報道中推測,熱控設(shè)計中采用的主要散熱措施包括表面涂覆散熱涂層和外殼開多孔兩種。外殼開多孔可以增加推力器放電室與冷空間的換熱,但是也會增強外部熱環(huán)境對推力器內(nèi)部環(huán)境的影響,同時也大大增加了機械加工的難度,因此不是優(yōu)選的設(shè)計方案。而表面涂覆散熱涂層能有效地提高相應(yīng)表面的散熱效率,改善推力器內(nèi)部熱量的散熱條件,是有效的設(shè)計方案。不過,國外資料中沒有給出散熱涂層的涂覆位置、散熱路徑和散熱面設(shè)計等具體熱控設(shè)計方案。國內(nèi)研究中,孫明明等[8]僅根據(jù)對LIPS-300離子電推力器的地面穩(wěn)態(tài)分析討論了表面散熱涂層的優(yōu)化作用,但沒有進行在軌熱分析,也沒有給出系統(tǒng)的熱控設(shè)計和優(yōu)化分析。

    為了提高LIPS-300離子電推力器的性能和可靠性,需對其進行系統(tǒng)的熱控設(shè)計和在軌熱分析,并對電推力器提出設(shè)計優(yōu)化建議。本文即針對典型的地球同步軌道位保工況對LIPS-300離子電推力器開展熱設(shè)計和在軌熱分析,并基于熱分析結(jié)果提出設(shè)計優(yōu)化方案。

    1 LIPS-300離子電推力器

    1.1 電推力器構(gòu)成和星上布局

    LIPS-300離子電推力器主要由主陰極、放電室、柵極組件、中和器、外殼(含安裝板)等部分組成。其中,放電室由陽極筒、屏柵筒、磁鐵和極靴組成,主陰極安裝在放電室底部屏柵筒中央,中和器陰極安裝在外殼外壁上。離子推力器內(nèi)部結(jié)構(gòu)如圖1所示。

    圖1 LIPS-300離子電推力器內(nèi)部結(jié)構(gòu)示意Fig.1 Schematic of the LIPS-300 ion thruster structure

    電推力器在星上的安裝位置如圖2所示。四臺LIPS-300離子電推力器對稱安裝于衛(wèi)星背地板外表面上。該布局電推力器包絡(luò)要求:電推力器在星上安裝后,距離星體南北板外表面的橫向尺寸和縱向距離背地板外表面的尺寸的最大包絡(luò)處均應(yīng)滿足運載包絡(luò)要求。

    LIPS-300離子電推力器通過柱段外殼安裝法蘭與矢量調(diào)節(jié)機構(gòu)的連接板隔熱安裝,矢量調(diào)節(jié)機構(gòu)通過安裝底板上的安裝孔與衛(wèi)星艙板隔熱安裝。

    1.2 內(nèi)部熱耗

    離子電推力器依靠氣體放電并中和產(chǎn)生推力,放電腔內(nèi)部各類帶有能量的粒子在沉積過程中會產(chǎn)生大量熱耗,并由等離子體云團向放電腔不同內(nèi)表面輻射[9-11]。此外,在離子電推力器工作時發(fā)射體需要維持較高的溫度才能克服逸出功并產(chǎn)生電子。

    孫明明等[12-13]基于能量沉積數(shù)學模型計算得到了LIPS-300離子電推力器在不同工作模式下放電室各個部件相應(yīng)位置的內(nèi)部熱耗以及發(fā)射體的溫度。在5.0 k W高功率模式下,總熱耗共585.47 W,其中陽極筒占393.35 W,主陰極、屏柵及加速柵分別加載剩余熱耗,發(fā)射體維持約1740℃高溫;在3.0 kW低功率模式下,總熱耗共396.45 W,其中陽極筒占294.40 W,發(fā)射體溫度約1734℃;在2.4 kW低功率模式下,總熱耗共326.12 W,其中陽極筒占242.81 W,發(fā)射體溫度約為1715℃。

    圖2 星上離子電推力器構(gòu)型布局Fig.2 Arrangement of the ion thruster on the satellite

    2 熱設(shè)計

    由于LIPS-300離子電推力器在軌工作期間放電室產(chǎn)生的高額熱流是制約熱設(shè)計的關(guān)鍵,因此由放電室至深冷空間設(shè)計有效的散熱路徑是實現(xiàn)LIPS-300離子電推力器溫度控制的主要技術(shù)途徑。根據(jù)上述設(shè)計原理,采取屏柵筒外表面、外殼內(nèi)表面高發(fā)射率處理,強化放電室到離子電推力器外殼的輻射換熱;在外殼外表面進行散熱設(shè)計,將工作期間放電室產(chǎn)生的廢熱排散到深冷空間。

    LIPS-300離子電推力器具體的熱控設(shè)計如下:

    1)上、下屏柵筒外表面(朝向外殼)進行本色陽極化處理(紅外發(fā)射率εh為0.82),強化屏柵筒與磁鋼、外殼之間的輻射換熱。

    2)前、后外殼內(nèi)表面(朝向屏柵筒外表面)進行本色陽極化處理(εh為0.82),強化外殼與屏柵筒、磁鋼之間的輻射換熱。

    3)前外殼外表面(朝向冷空間)含鎖緊耳片均噴涂CCAg(Silver Composite Coating)型復(fù)合熱控薄膜[14-15](太陽吸收率αS為0.3,εh為0.75),提高外殼外表面的散熱效率。

    4)后外殼外表面(朝向冷空間)包覆10單元多層隔熱組件,削弱電推力器與星本體之間的熱耦合。

    5)除此之外,電推力器其他表面為原始表面狀態(tài)。

    圖3為基于上述熱控設(shè)計的離子電推力器散熱路徑示意。從圖3中可以看到,本文開展的熱設(shè)計中在屏柵筒和外殼以及外殼和深冷空間之間相應(yīng)采用了本色陽極化處理和噴涂CCAg熱控涂層強化散熱,增強了LIPS-300離子電推力器大量熱耗的排散。此外,如圖3所示,上述熱控設(shè)計中沒有采用強化陽極筒和屏柵筒之間換熱的措施,后續(xù)將詳細分析此部分措施的影響。

    圖3 LIPS-300離子電推力器散熱路徑設(shè)計示意Fig.3 Schematic of the designed heat diffusion path of the LIPS-300 ion thruster

    3 有限元建模

    為便于分析計算,同時考慮到傳熱關(guān)系,在本文的有限元建模中采用了如下簡化與假設(shè):

    1)電推力器安裝在星外,星體內(nèi)部溫度對電推力器的影響只通過電推力器安裝底板與星體接觸面來實現(xiàn),本文模型中把星體安裝板作定壁溫邊界處理,忽略衛(wèi)星本體艙內(nèi)溫度波動對電推力器的影響。

    2)矢量調(diào)節(jié)機構(gòu)按定溫邊界處理,外表面按包覆10單元多層隔熱組件處理。

    3)忽略離子電推力器與矢量調(diào)節(jié)機構(gòu)之間的接觸換熱,按絕熱處理。

    4)電推力器放電室熱流分布按照理論分析計算賦值[12-13],中和器陰極和主陰極發(fā)射體按定溫邊界處理。

    5)屏柵和加速柵為帶孔拋面,由于柵極開孔孔徑較小無法精確建模,因此建模時簡化為實體,然后根據(jù)開孔率折算等效的太陽吸收比和紅外發(fā)射率。

    6)忽略東西板和對地板上的星外凸出部件以及太陽翼對電推進分系統(tǒng)星外部件的溫度影響。

    7)忽略LIPS-300離子電推力器與星體之間連接電纜、管路的發(fā)熱和漏熱。

    8)有限元模型中簡化了與陰極相連的管路、電連接導(dǎo)線、表面連接處的螺孔螺釘?shù)燃毿⊥蛊鸩考?/p>

    對于在軌熱分析,需要對電推力器和星體進行整體有限元建模。其中,主陰極處部件結(jié)構(gòu)和熱耦合關(guān)系復(fù)雜,且發(fā)射體有高溫邊界條件,為了保證有限元模型的準確性,建模中采用自動和人工結(jié)合的網(wǎng)格劃分方法,分別建立了陰極發(fā)射體、陰極管、內(nèi)外陰極屏套等有限元網(wǎng)格及熱耦合關(guān)系。此外,上極靴、中極靴、下極靴是熱分析的關(guān)鍵部件且結(jié)構(gòu)形式復(fù)雜,建模中同樣采用自動和人工結(jié)合的網(wǎng)格劃分方法。

    星本體東西艙板及推力器后外殼表面包覆的多層隔熱組件,采用不均勻多層殼的建模方法。即直接在實體表面繪制面網(wǎng)格,并賦予不均勻多層殼屬性,不均勻多層殼的頂層和底層的輻射屬性分別按多層組件表面(即黃膜或黑膜)和實體的實際材料賦值,頂層和底層之間按多層隔熱組件的當量輻射率賦輻射熱耦合關(guān)系。此方法未新增單元和節(jié)點,計算效率較高。

    圖4為本文建立的LIPS-300離子電推力器和星體的有限元模型,模型共劃分為111 697個面(殼)網(wǎng)格單元。

    圖5給出了LIPS-300離子電推力器的有限元模型,總共劃分為27 250個面(殼)網(wǎng)格單元。

    圖4 整星有限元模型Fig.4 Finite element model of the whole-spacecraft

    圖5 LIPS-300離子電推力器有限元模型Fig.5 Finite element model of the LIPS-300 ion thruster

    由于LIPS-300離子電推力器還未開展熱平衡試驗,本文僅能通過與孫明明等[8]開展的LIPS-300離子電推力器的地面穩(wěn)態(tài)分析進行比對,驗證文中建立的熱分析模型。孫明明等使用ANSYS有限元軟件建模,本文使用NX/ TMG有限元軟件建模,兩者得到的結(jié)果可以相互驗證。不過,由于孫明明等建模的關(guān)注點在熱變形分析,與本文建模的關(guān)注點和簡化假設(shè)有所區(qū)別,因此本文熱分析結(jié)果與其分析結(jié)果的對比會存在一定的差異。

    對比驗證工況取與孫明明等[8]地面穩(wěn)態(tài)分析相同的邊界條件,同樣設(shè)置推力器工作條件為地面環(huán)境,溫度為22℃。圖6給出了本文熱分析模型和孫明明等[8]的地面穩(wěn)態(tài)分析模型得到的下磁鋼及下極靴的溫度分布,從圖6中對比可以看出兩者的溫度分布規(guī)律相同且差異小于10℃,驗證了本文熱分析模型的合理性。

    圖6 下磁鋼及下極靴溫度分布結(jié)果Fig.6 Temperature distribution of the downstream magnet and downstream anode

    4 在軌熱分析

    4.1 在軌熱環(huán)境

    對于LIPS-300離子電推力器典型地球同步軌道的位保工況,工作軌道參數(shù)為:軌道高度35 786 km,軌道周期24 h,軌道傾角0°,最長地球陰影時間72 min。由于軌道高度遠大于10 000 km,衛(wèi)星受到的地球輻射和地球反照可忽略不計,衛(wèi)星外熱流主要為太陽輻射,太陽輻照平均值為1 353 W/m2,春分、夏至、秋分、冬至是外熱流變化的典型特征日,太陽輻照值有所不同,在夏至為1 309 W/m2,冬至為1 399 W/m2,分點為1 361 W/m2。

    4.2 參數(shù)設(shè)置

    在熱分析中選取的電推力器主要部件的熱物性參數(shù)如表1所示。

    表1 LIPS-300離子電推力器材料屬性Table 1 Material properties of the LIPS-300 ion thruster

    4.3 工況設(shè)置

    本文計算的衛(wèi)星處于地球同步軌道,影響衛(wèi)星溫度變化的空間環(huán)境熱源是太陽。經(jīng)分析,離子電推力器極端低溫工況出現(xiàn)在轉(zhuǎn)移軌道,無外熱流的情況。由于冬至太陽光照射南板且太陽輻照度最大,再加上壽命末期熱控涂層退化,太陽吸收比升高,因此壽命末期的冬至靠近南板安裝的離子電推力器將出現(xiàn)高溫工況。

    綜上所述,本文進行以下2個典型高溫、低溫工況的熱分析計算如下。

    工況1:壽命末期,冬至軌道外熱流,瞬態(tài)計算;+X側(cè)2臺離子電推力器24 h工作,離子電推力器內(nèi)部功耗按照5.0 k W高功率模式的熱流分布賦值,其余2臺離子電推力器均不工作;衛(wèi)星本體艙板定溫為40℃,矢量調(diào)節(jié)機構(gòu)定溫為60℃。

    工況2:轉(zhuǎn)移軌道,無外熱流,穩(wěn)態(tài)計算; 4個離子電推力器均不工作;衛(wèi)星本體艙板定溫為0℃,矢量調(diào)節(jié)機構(gòu)定溫為0℃。

    4.4 計算結(jié)果及分析

    選取安裝衛(wèi)星背地板+X+Y側(cè)的離子電推力器為主分析對象,其他3個離子電推力器為輔助分析對象,考慮其對主分析對象的遮擋影響。

    工況1和工況2主要部件溫度平均值Tave(部件各個位置溫度的平均值)見表2,其中工況1還列出了部件溫度的最大值Tmax和最小值Tmin(部件各個位置溫度的最大值和最小值)。

    表2 工況1、工況2主要部件溫度計算結(jié)果Table 2 Simulation results of the main components temperature in Case 1 and Case 2 ℃

    從計算結(jié)果可以看到,在低溫工況,電推力器整體溫度范圍約為-191~-176℃,其中加速柵和支撐座由于與其他部件隔熱設(shè)計,因此溫度較電推力器其他部件低約10℃,電推力器整體溫度在-179℃左右。

    在離子電推力器在軌工作的高溫工況下,下極靴、下陽極筒、上陽極筒、屏柵、下磁鋼和下屏柵筒為主要的高溫位置,具體分析如下:

    1)主陰極有大量內(nèi)部熱耗,并且工作中主陰極維持約1 740℃恒定高溫,受其影響下極靴、下陽極筒溫度較高。此外,對于下陽極筒還有內(nèi)部熱耗影響,因此下極靴、下陽極筒為高溫位置。

    2)上陽極筒、屏柵受內(nèi)部熱耗以及下陽極筒和主陰極的熱輻射影響溫度也較高。

    3)下磁鋼和下屏柵筒沒有內(nèi)部熱耗,不過它們與下極靴和下陽極筒有較強的熱耦合(包括熱傳導(dǎo)和熱輻射),因此下磁鋼和下屏柵筒也是明顯的高溫位置。

    考慮到磁鋼的選材和工作性能要求,磁鋼溫度上限為380℃,下磁鋼的溫度最大值已經(jīng)超過了溫度上限要求。為了更清晰地描述磁鋼在離子電推力器在軌工作高溫工況下的溫度變化和溫度分布情況,圖7和圖8分別給出了工況1下磁鋼在兩個周期內(nèi)溫度T隨時間t的變化曲線和磁鋼在當?shù)卣鐣r刻的溫度云圖。

    圖7 工況1下磁鋼溫度隨時間的變化曲線Fig.7 Temperature evolutions of the downstream magnet in case 1

    從圖7可以看到,在離子電推力器在軌工作的高溫工況下,下磁鋼溫度經(jīng)歷約4 h的快速上升后,穩(wěn)定小幅波動,波動幅度約20℃。磁鋼初始快速升溫主要是由于電推力器點火工作后內(nèi)部熱耗的影響,約4 h后內(nèi)部熱耗對溫度的影響已趨于穩(wěn)態(tài),隨后的溫度波動為外熱流的影響。由于離子電推力器在軌工作的內(nèi)部熱耗較大,引起的溫升遠大于外熱流引起的溫度波動,因此溫度隨時間的變化曲線表現(xiàn)為如圖7所示的溫度快速升高后趨于穩(wěn)定。計算得到的其他部件溫度隨時間變化曲線也有類似的規(guī)律。上述模擬結(jié)果表明在5.0 k W高功率模式下外熱流對電推力器內(nèi)部部件溫度影響較小,內(nèi)熱源是內(nèi)部部件溫度水平的主導(dǎo)因素。

    圖8 工況1磁鋼當?shù)卣鐣r刻溫度云圖Fig.8 Temperature contour of the magnets at the local noon in case 1

    另外,從圖8可以看到靠近下極靴的磁鋼部位是磁鋼的高溫區(qū)域,也是超過磁鋼溫度上限要求的位置。因此,從熱分析的角度設(shè)計優(yōu)化主要有兩個方向:一是增強后外殼向深冷空間散熱從而直接降低下磁鋼溫度;二是增強下極靴和下陽極筒向外的傳熱通道(包括熱輻射和熱傳導(dǎo))使它們的溫度降低,進而減弱它們對下磁鋼的影響使下磁鋼溫度降低。不過,兩個設(shè)計優(yōu)化方向的效果還需要通過模擬計算進行驗證。下節(jié)將從增加散熱面積、增強內(nèi)部換熱兩個方面展開設(shè)計優(yōu)化分析。

    5 設(shè)計優(yōu)化分析

    5.1 增加散熱面積

    在第2節(jié)所述熱控設(shè)計基礎(chǔ)上,增加設(shè)置錐段后外殼外表面為散熱面,即錐段后外殼外表面熱控狀態(tài)由包覆多層隔熱組件改為涂覆CCAg型復(fù)合熱控薄膜[14-15],作為設(shè)計優(yōu)化狀態(tài)1。不過,增加散熱面會影響低溫工況部件的溫度,因此也需要對低溫工況進行核算。這里設(shè)計優(yōu)化狀態(tài)1同樣設(shè)置高溫工況和低溫工況,高溫工況與工況1設(shè)置一致,低溫工況與工況2設(shè)置一致。設(shè)計優(yōu)化狀態(tài)1高溫工況、低溫工況各部件溫度平均值見表3,其中高溫工況還列出了部件溫度的最大值和最小值。

    表3 設(shè)計優(yōu)化狀態(tài)1高溫工況、低溫工況主要部件溫度計算結(jié)果Table 3 Simulation results of the main components temperature for optimization design scheme 1 in the worst hot and cold cases ℃

    從計算結(jié)果可以看到,相比工況1,設(shè)計優(yōu)化狀態(tài)1高溫工況電推力器整體溫度明顯降低,磁鋼溫度降幅約13%~19%,上、下磁鋼溫度均低于380℃。因此,增加散熱面積對降低磁鋼溫度的優(yōu)化效果顯著。

    此外,相比工況2,設(shè)計優(yōu)化狀態(tài)1低溫工況電推力器整體溫度并未降低而是提高約19℃。這是由于錐段后外殼外表面涂覆CCAg型復(fù)合熱控薄膜在增強與深冷空間換熱的同時也會增強電推力器與星本體的熱交換,而星本體艙板溫度受星內(nèi)熱控措施控制,低溫下限為0℃,對應(yīng)多層外表面約為-100℃,高于低溫工況(工況2)電推力器部件溫度,因此對于設(shè)計優(yōu)化狀態(tài)1,低溫工況下,電推力器會從星本體吸熱從而使電推力器部件溫度有所升高,這對于電推力器部件性能也是有利的。

    綜上所述,適當增加散熱面積是有效的設(shè)計優(yōu)化方向。

    5.2 增加內(nèi)部換熱

    下極靴和下陽極筒向外的傳熱通道主要包括熱輻射通道和熱傳導(dǎo)通道,因此可以通過增強極靴、陽極筒與屏柵筒之間的內(nèi)部輻射換熱和傳導(dǎo)換熱來增強內(nèi)部換熱,從而增強下極靴和下陽極筒向外的傳熱通道,實現(xiàn)降低其溫度的效果。

    (1)增強內(nèi)部輻射換熱

    在第2節(jié)所述熱控設(shè)計基礎(chǔ)上,增加將屏柵筒內(nèi)表面與陽極筒外表面做陽極化處理,作為設(shè)計優(yōu)化狀態(tài)2,即加強陽極筒與屏柵筒之間輻射換熱。計算工況設(shè)置高溫工況,與工況1設(shè)置一致,關(guān)注設(shè)計優(yōu)化狀態(tài)2對高溫區(qū)是否有降溫效果。設(shè)計優(yōu)化狀態(tài)2高溫工況各部件溫度平均值、最大值和最小值見表4。

    表4 設(shè)計優(yōu)化狀態(tài)2高溫工況主要部件溫度計算結(jié)果Table 4 Simulation results of the main components temperature for optimization design scheme 2 in the worst hot case ℃

    從計算結(jié)果可以看到,相比工況1,設(shè)計優(yōu)化狀態(tài)2高溫工況雖然陽極筒溫度明顯降低,降幅約16%~20%,但屏柵筒溫度升高,升幅約5%~6%。下磁鋼受屏柵筒影響溫度升高約10℃,下磁鋼工作溫度范圍為359.91~398.91℃,溫度最大值超過380℃的溫度上限要求,因此設(shè)計優(yōu)化狀態(tài)2沒有實現(xiàn)磁鋼降溫的優(yōu)化效果,增強陽極筒與屏柵筒之間的內(nèi)部輻射換熱不是有效的設(shè)計優(yōu)化方向。

    (2)增強內(nèi)部傳導(dǎo)換熱

    在第2節(jié)所述熱控狀態(tài)設(shè)計基礎(chǔ)上,將屏柵筒與極靴間的接觸面積增加一倍,作為設(shè)計優(yōu)化狀態(tài)3,即加強極靴與屏柵筒之間傳導(dǎo)換熱。計算工況同樣設(shè)置高溫工況,與工況1設(shè)置一致,關(guān)注設(shè)計優(yōu)化狀態(tài)3對高溫區(qū)是否有降溫效果。設(shè)計優(yōu)化狀態(tài)3高溫工況各部件溫度平均值、最大值和最小值見表5。

    表5 設(shè)計優(yōu)化狀態(tài)3高溫工況主要部件溫度計算結(jié)果Table 5 Simulation results of the main components temperature for optimization design scheme 3 in the worst hot case ℃

    從計算結(jié)果可以看到,相比工況1,設(shè)計優(yōu)化狀態(tài)3高溫工況電推力器整體溫度水平變化較小,磁鋼溫度降低不超過5℃,溫度最大值還是超過380℃的溫度上限要求。這是由于電推力器內(nèi)部部件之間換熱以熱輻射為主,熱傳導(dǎo)對各部件溫度水平的影響較小。因此,增強極靴與屏柵筒之間的內(nèi)部傳導(dǎo)換熱也不是有效的設(shè)計優(yōu)化方向。

    綜上所述,增強陽極筒、極靴與屏柵筒間的內(nèi)部換熱不是有效的設(shè)計優(yōu)化方向。

    6 結(jié)束語

    本文對地球同步軌道LIPS-300離子電推力器開展了熱設(shè)計和在軌熱分析。熱設(shè)計采用了本色陽極化處理和噴涂CCAg熱控涂層強化屏柵筒和外殼以及外殼和深冷空間之間換熱,從而改善了電推力器內(nèi)部熱耗的散熱路徑。在軌熱分析結(jié)果表明:在極端高溫工況下,下極靴、陽極筒、屏柵、下磁鋼為主要的高溫位置,且下磁鋼最高溫度會超出磁鋼工作溫度上限。此外,電推力器主要部件溫度隨時間的變化規(guī)律均為在初始快速升溫后穩(wěn)定小幅波動,表明電推力器內(nèi)熱源是決定電推力器整體溫度水平的主導(dǎo)因素。基于在軌熱分析,本文還驗證了兩個設(shè)計優(yōu)化方向的有效性,結(jié)果表明適當增加散熱面積是有效的設(shè)計優(yōu)化方向,而增強陽極筒、極靴與屏柵筒間的內(nèi)部換熱是無效的設(shè)計優(yōu)化方向。

    References)

    [1] 溫正,鐘凌偉,王一白,等.離子推力器加速柵極離子運動規(guī)律的數(shù)值研究[J].強激光與離子束,2011,23(6): 1640-1645. WEN Z,ZHONG L W,WANG Y B,et al.Three dimensional numerical study on motion lows of ions in ion thruster optics[J].High Power Laser and Particle Beams,2011,23(6):1640-1645(in Chinese).

    [2] 陳娟娟,張?zhí)炱?賈艷輝,等.20 cm氙離子推力器放電室性能優(yōu)化[J].強激光與離子束,2012,24(10): 2469-2474. CHEN J J,ZHANG T P,JIA Y H,et al. Performance optimization of 20 cm xenon ion thruster discharge chamber[J].High Power Laser and Particle Beams,2012,24(10):2469-2474(in Chinese).

    [3] ZHANG T P,WANG X Y,JIANG H C.Initial flight test results of the LIPS-200 electric propulsion system on SJ-9A satellite[C]∥International Electric Propulsion Conference,2013:47-54.

    [4] 張?zhí)炱?蘭州空間技術(shù)物理研究所電推進新進展[J].火箭推進,2015,41(2):7-12. ZHANG T P.New progress of electric propulsion technology in LIP[J].Journal of Rocket Propulsion, 2015,41(2):7-12(in Chinese).

    [5] WEN L,CROTTY J D,PAWLIK E V.Ion thruster thermal characteristics and performance[J].Journal of Spacecraft and Rockets,1973,10(1):35-41.

    [6] OGLEBAY J C.Thermal analytic model of a 30 cm engineering model mercury ion thruster[C]∥International Electric Propulsion Conferenee,1975:344-360.

    [7] 陳茂林,夏廣慶,徐宗琦,等.柵極熱變形對離子推力器工作過程影響分析[J].物理學報,2015(9):094104. CHEN M L,XIA G Q,XU Z Q,et al.Analysis on the effects of optics thermal deformation on the ion thruster op eration[J].Acta Physica Sinica,2015(9): 094104(in Chinese).

    [8] 孫明明,張?zhí)炱?王亮.30 cm口徑離子推力器熱特性模擬分析[J].真空與低溫,2014,20(3):158-162. SUN M M,ZHANG T P,WANG L.Thermal analysis of 30 cm diameter thruster[J].Vacuum and Cryogenics,2014,20(3):158-162(in Chinese).

    [9] 賀武生,孫安邦,毛根旺,等.離子推力器放電腔數(shù)值模擬[J].強激光與離子束,2010,22(12):3020-3024. HE W S,SUN A B,MAO G W,et al.Numerical simulation of ion thruster discharge chamber[J].High Power Laser and Particle Beams,2010,22(12):3020-3024(in Chinese).

    [10] 張?zhí)炱?陳繼巍,李小平,等.地面壽命試驗中離子推力器表面的濺射沉積量計算[J].航天器環(huán)境工程, 2011,28(5):436-439. ZHANG T P,CHEN J W,LI X P,et al. Calculation of backsputter depositon on LIPS-200 ion thruster in the life test[J].Spacecraft Environment Engineering,2011,28(5):436-439(in Chinese).

    [11] 孫安邦,毛根旺,夏廣慶,等.離子推力器放電腔內(nèi)等離子體流動規(guī)律的全粒子模型[J].推進技術(shù), 2012,33(1):143-149. SUN A B,MAO G W,XIA G Q,et al.Full particles model of plasma flow for ion thruster discharge chamber[J].Journal of Propulsion Technology,2012, 33(1):143-149(in Chinese).

    [12] 孫明明,顧左,郭寧,等.離子推力器空心陰極熱特性模擬分析[J].強激光與離子束,2010,22(5): 1149-1152. SUN M M,GU Z,GUO N,et al.Thermal analysis of hollow cathodes for ion thruster[J].High Power Laser and Particle Beams,2010,22(5): 1149-1152(in Chinese).

    [13] 孫明明,張?zhí)炱?陳娟娟,等.LIPS-200離子推力器熱特性模擬分析研究[J].強激光與離子束,2014,26 (8):084002. SUN M M,ZHANG T P,CHEN J J,et al.Thermal analysis of LIPS-200 ion thruster[J].High Power Laser and Particle Beams,2014,26(8):084002(in Chinese).

    [14] DEVER J,MILLER S.Exposure of polymer film thermal control materials on the materials international space station experiment(MISSE),AIAA-2002-4924[R]. Reston:AIAA,2002.

    [15] 羅列超,趙榮根,孟佳,等.航天器用鋁光亮陽極氧化涂層特性研究[J].無機材料學報,2002,17(6): 1269-1276. LUO L C,ZHAO R G,MENG J,et al.White anodic coating on aluminium surface for spacecraft[J]. Journal of Inorganic Materials,2002,17(6):1269-1276 (in Chinese).

    (編輯:范真真)

    Thermal design and optimization of LlPS-300 ion thruster in geosynchronous orbit

    HU Guojie*,LI Jian,LIU Bailin
    Institute of Telecommunication Satellite,China Academy of Space Technology,Beijing 100094,China

    The Large heat dissipation of LIPS-300 ion thruster seriously affects its performance.Therefore,the thermal design and analysis of LIPS-300 ion thruster in geosynchronous orbit were carried out.The simulation results show that under the hightemperature working condition,the downstream harness,downstream and upstream anode, screen grid,and downstream magnet are the positions with the highest temperature,and the temperature of the downstream magnet exceeds the temperature upper limit.In order to decrease the temperature of the downstream magnet,two optimization design measures were presented and tried.The analysis results indicate that increasing heat radiation area is an effective optimization design measure while enhancing the inner heat exchange is noneffective.

    V439+.4

    :A

    10.3780/j.issn.1000-758X.2016.0015

    2015-11-12;

    :2015-12-10;錄用日期:2015-12-30;< class="emphasis_bold">網(wǎng)絡(luò)出版時間

    時間:2016-02-24 13:41:18

    http:∥www.cnki.net/kcms/detail/11.1859.V.20160224.1341.013.html

    *

    :胡幗杰(1985-),女,博士,hgj919@163.com,主要研究方向為航天器熱設(shè)計及換熱流動優(yōu)化

    胡幗杰,李健,劉百麟.地球靜止軌道LIPS-300離子推力器熱設(shè)計與優(yōu)化[J].中國空間科學技術(shù),2016,36(1):85-93.HU G J,LI J,LIU B L.Thermal design and optimization of LIPS-300 ion thruster in geosynchronous orbit[J]. Chinese Space Science and Technology,2016,36(1):85-93(in Chinese).

    http:∥zgkj.cast.cn

    猜你喜歡
    磁鋼推力器外殼
    單組元推力器倒置安裝多余物控制技術(shù)驗證
    小型電機環(huán)形磁鋼粘接及其粘接強度提高
    基于Halbach陣列磁鋼的PMSM氣隙磁密波形優(yōu)化
    防爆電機(2020年4期)2020-12-14 03:11:08
    U盾外殼組件注塑模具設(shè)計
    塑料外殼注射模設(shè)計
    模具制造(2019年7期)2019-09-25 07:30:00
    基于有限元分析的表貼式磁性轉(zhuǎn)子工藝技術(shù)研究
    微特電機(2019年1期)2019-02-25 02:09:44
    隔爆外殼水壓試驗工藝探討
    用于小行星探測的離子推力器技術(shù)研究
    離子推力器和霍爾推力器的異同
    太空探索(2015年6期)2015-07-12 12:48:42
    一種新的降低渦流損耗的磁鋼分段方式
    久久国产乱子伦精品免费另类| 亚洲精品国产区一区二| 欧美久久黑人一区二区| 久久久久久亚洲精品国产蜜桃av| 香蕉国产在线看| 在线看a的网站| 国产精品爽爽va在线观看网站 | 国产一区二区三区综合在线观看| 男女午夜视频在线观看| 精品国产亚洲在线| 久久 成人 亚洲| 高清av免费在线| 满18在线观看网站| 啦啦啦免费观看视频1| 亚洲欧美激情在线| 不卡av一区二区三区| 美女福利国产在线| 可以在线观看毛片的网站| 天堂影院成人在线观看| 日本wwww免费看| 久久九九热精品免费| 国产成人影院久久av| 国产欧美日韩精品亚洲av| 色婷婷久久久亚洲欧美| 高清在线国产一区| 久久久久久亚洲精品国产蜜桃av| 久99久视频精品免费| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品久久成人aⅴ小说| 精品久久久久久,| 中文字幕人妻丝袜一区二区| 日本 av在线| 欧美+亚洲+日韩+国产| 黑人巨大精品欧美一区二区mp4| 色播在线永久视频| 国产精品 国内视频| 久久精品亚洲精品国产色婷小说| 国产成人av激情在线播放| 丰满人妻熟妇乱又伦精品不卡| 欧洲精品卡2卡3卡4卡5卡区| 搡老乐熟女国产| 国产欧美日韩一区二区精品| 十八禁人妻一区二区| 最近最新免费中文字幕在线| 日本一区二区免费在线视频| 可以免费在线观看a视频的电影网站| 老鸭窝网址在线观看| 中文字幕另类日韩欧美亚洲嫩草| 免费观看人在逋| 国产精品爽爽va在线观看网站 | 欧美午夜高清在线| 久久国产乱子伦精品免费另类| 日本五十路高清| 日本免费一区二区三区高清不卡 | 亚洲国产精品合色在线| 国产精品乱码一区二三区的特点 | 亚洲男人的天堂狠狠| 91麻豆av在线| 欧美中文综合在线视频| 国产精品野战在线观看 | 免费少妇av软件| 一进一出好大好爽视频| 99国产综合亚洲精品| 麻豆av在线久日| 女人被狂操c到高潮| 在线观看免费日韩欧美大片| 99香蕉大伊视频| 婷婷丁香在线五月| 精品久久久久久成人av| 久久影院123| 久久青草综合色| 美女福利国产在线| 看免费av毛片| 超碰成人久久| 久久人人97超碰香蕉20202| 国产欧美日韩一区二区三区在线| 女人高潮潮喷娇喘18禁视频| 村上凉子中文字幕在线| 国产亚洲精品第一综合不卡| 黑人巨大精品欧美一区二区蜜桃| 日韩欧美国产一区二区入口| 国产又爽黄色视频| 欧美黄色片欧美黄色片| 国产精品av久久久久免费| 国产成人精品无人区| 国产欧美日韩综合在线一区二区| 国产精品美女特级片免费视频播放器 | 91av网站免费观看| 嫩草影院精品99| 久久中文字幕人妻熟女| 亚洲国产精品一区二区三区在线| 亚洲人成77777在线视频| 黄色女人牲交| 日韩欧美一区视频在线观看| 国产高清激情床上av| av视频免费观看在线观看| 神马国产精品三级电影在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 国产精品1区2区在线观看.| 丝袜美足系列| 久久国产精品人妻蜜桃| 中文字幕最新亚洲高清| 国产99久久九九免费精品| 亚洲精品av麻豆狂野| 91九色精品人成在线观看| 久久青草综合色| 成人精品一区二区免费| 亚洲av成人一区二区三| 丝袜美足系列| 好看av亚洲va欧美ⅴa在| 真人做人爱边吃奶动态| 脱女人内裤的视频| 国产精品电影一区二区三区| 级片在线观看| 久久热在线av| 精品国产乱子伦一区二区三区| 黄色a级毛片大全视频| 国产精品电影一区二区三区| 久久久久久久久中文| 国产亚洲精品一区二区www| 色播在线永久视频| 两性午夜刺激爽爽歪歪视频在线观看 | 日日干狠狠操夜夜爽| 在线看a的网站| 母亲3免费完整高清在线观看| 国内毛片毛片毛片毛片毛片| 日韩欧美免费精品| 亚洲av第一区精品v没综合| 国产精品久久久久久人妻精品电影| 欧美激情极品国产一区二区三区| 人人妻人人添人人爽欧美一区卜| 99香蕉大伊视频| 18美女黄网站色大片免费观看| a在线观看视频网站| 一级作爱视频免费观看| 69精品国产乱码久久久| 波多野结衣一区麻豆| 青草久久国产| 久久久久久久久久久久大奶| 亚洲五月婷婷丁香| 青草久久国产| 大香蕉久久成人网| 久久午夜亚洲精品久久| 在线观看免费日韩欧美大片| 女警被强在线播放| 亚洲第一欧美日韩一区二区三区| www.精华液| 国产在线精品亚洲第一网站| xxxhd国产人妻xxx| 看片在线看免费视频| 丁香六月欧美| 搡老熟女国产l中国老女人| 国产免费男女视频| 国产精品免费视频内射| 曰老女人黄片| 中文字幕精品免费在线观看视频| 日韩大尺度精品在线看网址 | 日本五十路高清| 乱人伦中国视频| 国产亚洲精品第一综合不卡| 欧美日韩黄片免| 久久久国产欧美日韩av| 精品一品国产午夜福利视频| a级片在线免费高清观看视频| 欧美日本中文国产一区发布| 久久精品国产亚洲av高清一级| 日本免费a在线| 国产av精品麻豆| 免费搜索国产男女视频| 一区二区三区精品91| 麻豆国产av国片精品| 国产区一区二久久| 精品第一国产精品| 麻豆一二三区av精品| 岛国在线观看网站| 身体一侧抽搐| 韩国精品一区二区三区| 亚洲少妇的诱惑av| 在线视频色国产色| 亚洲精品粉嫩美女一区| 免费女性裸体啪啪无遮挡网站| 国产精华一区二区三区| 国产成人欧美在线观看| 国产激情久久老熟女| 人妻久久中文字幕网| 最新美女视频免费是黄的| 最近最新中文字幕大全电影3 | 亚洲午夜精品一区,二区,三区| 欧美日韩av久久| 女人被躁到高潮嗷嗷叫费观| 丝袜在线中文字幕| 国产av在哪里看| 精品人妻在线不人妻| 亚洲av成人av| 精品免费久久久久久久清纯| 大香蕉久久成人网| 99久久国产精品久久久| 高清在线国产一区| 两个人看的免费小视频| 免费一级毛片在线播放高清视频 | 丰满迷人的少妇在线观看| 国产精品久久久久成人av| 婷婷精品国产亚洲av在线| 亚洲性夜色夜夜综合| 国产精品美女特级片免费视频播放器 | 热re99久久精品国产66热6| 成人国产一区最新在线观看| 久久精品影院6| 成人手机av| 人人妻人人澡人人看| 在线观看免费高清a一片| 亚洲一区中文字幕在线| 免费少妇av软件| 母亲3免费完整高清在线观看| 两个人看的免费小视频| 999久久久精品免费观看国产| 巨乳人妻的诱惑在线观看| videosex国产| 日本vs欧美在线观看视频| 国产熟女xx| 亚洲成人免费电影在线观看| 欧美黑人精品巨大| 99riav亚洲国产免费| 在线av久久热| 久久精品亚洲熟妇少妇任你| 国产无遮挡羞羞视频在线观看| 老汉色∧v一级毛片| 欧美激情 高清一区二区三区| 日韩高清综合在线| 成人18禁高潮啪啪吃奶动态图| 夜夜看夜夜爽夜夜摸 | 亚洲黑人精品在线| av网站在线播放免费| 久久国产精品人妻蜜桃| 国产精品亚洲一级av第二区| 成人国产一区最新在线观看| 午夜福利,免费看| 久久人妻福利社区极品人妻图片| 久久久久久久久久久久大奶| 黄色丝袜av网址大全| 国产乱人伦免费视频| 日韩成人在线观看一区二区三区| 久久精品影院6| 久久精品91无色码中文字幕| 久久久久国内视频| 成年人黄色毛片网站| 国产区一区二久久| 黄色怎么调成土黄色| 国产精品1区2区在线观看.| 丁香欧美五月| 美女高潮到喷水免费观看| 欧美日韩视频精品一区| 黄色怎么调成土黄色| 午夜视频精品福利| 亚洲精品美女久久av网站| 精品熟女少妇八av免费久了| 丰满迷人的少妇在线观看| 国产一区二区三区视频了| 91麻豆精品激情在线观看国产 | 亚洲av成人一区二区三| 成年版毛片免费区| 侵犯人妻中文字幕一二三四区| 国产免费av片在线观看野外av| 精品免费久久久久久久清纯| 国产精品秋霞免费鲁丝片| 久久精品国产亚洲av香蕉五月| 国产一区二区三区在线臀色熟女 | av天堂在线播放| 又大又爽又粗| 欧美精品亚洲一区二区| 一个人免费在线观看的高清视频| 脱女人内裤的视频| 亚洲成人国产一区在线观看| 国产av一区在线观看免费| 色在线成人网| 亚洲欧美激情综合另类| 免费在线观看影片大全网站| 久久亚洲真实| 久久精品亚洲熟妇少妇任你| 国产有黄有色有爽视频| 男男h啪啪无遮挡| 欧美午夜高清在线| 日韩欧美在线二视频| 免费看a级黄色片| 亚洲色图 男人天堂 中文字幕| 欧美日韩精品网址| 成人亚洲精品一区在线观看| 日日夜夜操网爽| 男女做爰动态图高潮gif福利片 | 色综合站精品国产| 最新美女视频免费是黄的| 级片在线观看| 精品国产国语对白av| 黄网站色视频无遮挡免费观看| 欧美一区二区精品小视频在线| 国产精品综合久久久久久久免费 | 国产欧美日韩一区二区精品| 国产成人精品久久二区二区91| 精品国产乱码久久久久久男人| 日韩三级视频一区二区三区| av在线播放免费不卡| 在线视频色国产色| 日韩av在线大香蕉| 人人妻人人添人人爽欧美一区卜| 精品电影一区二区在线| 黄色视频,在线免费观看| 午夜久久久在线观看| 亚洲欧美日韩另类电影网站| 无遮挡黄片免费观看| 国产亚洲欧美在线一区二区| 最近最新免费中文字幕在线| 在线观看午夜福利视频| 亚洲欧美一区二区三区黑人| 亚洲精品美女久久久久99蜜臀| 日本精品一区二区三区蜜桃| 国产亚洲精品久久久久5区| 国产免费现黄频在线看| 19禁男女啪啪无遮挡网站| 免费人成视频x8x8入口观看| 性欧美人与动物交配| 曰老女人黄片| 亚洲久久久国产精品| 香蕉丝袜av| 大香蕉久久成人网| 交换朋友夫妻互换小说| 成人影院久久| 欧美日韩av久久| 欧美黑人精品巨大| 亚洲aⅴ乱码一区二区在线播放 | 免费av中文字幕在线| 欧美日韩精品网址| 亚洲成av片中文字幕在线观看| 精品欧美一区二区三区在线| 欧美日本亚洲视频在线播放| 最新在线观看一区二区三区| 男女午夜视频在线观看| 极品人妻少妇av视频| 伊人久久大香线蕉亚洲五| 在线观看日韩欧美| 国产精品综合久久久久久久免费 | 久久久久久亚洲精品国产蜜桃av| 国产精品影院久久| 性少妇av在线| 免费看a级黄色片| 黄色片一级片一级黄色片| 纯流量卡能插随身wifi吗| 久久香蕉激情| 成人国产一区最新在线观看| 操美女的视频在线观看| 成人三级黄色视频| 狂野欧美激情性xxxx| 91成人精品电影| aaaaa片日本免费| ponron亚洲| 在线观看免费日韩欧美大片| 亚洲午夜理论影院| 成人av一区二区三区在线看| 久久久久精品国产欧美久久久| 91在线观看av| 国产高清激情床上av| 桃红色精品国产亚洲av| 国产成人精品久久二区二区免费| 国产在线观看jvid| 日韩大码丰满熟妇| 亚洲av五月六月丁香网| 69av精品久久久久久| www日本在线高清视频| 午夜亚洲福利在线播放| 久久久久国产精品人妻aⅴ院| 日日干狠狠操夜夜爽| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩乱码在线| 久久精品亚洲精品国产色婷小说| 天堂中文最新版在线下载| 欧美乱色亚洲激情| 大型av网站在线播放| 精品福利观看| 99riav亚洲国产免费| 国产视频一区二区在线看| 丝袜人妻中文字幕| 亚洲精品国产区一区二| 国产国语露脸激情在线看| 777久久人妻少妇嫩草av网站| 久久久久久大精品| 亚洲国产中文字幕在线视频| ponron亚洲| 一夜夜www| 夜夜夜夜夜久久久久| 日日干狠狠操夜夜爽| 久久人妻熟女aⅴ| 亚洲欧美激情综合另类| 免费在线观看影片大全网站| 免费搜索国产男女视频| 亚洲人成网站在线播放欧美日韩| 中文字幕人妻丝袜制服| 久久久久国产精品人妻aⅴ院| 天天添夜夜摸| 国产亚洲精品一区二区www| 亚洲av美国av| 人妻丰满熟妇av一区二区三区| 国产在线观看jvid| 18禁黄网站禁片午夜丰满| 国产成人欧美| 手机成人av网站| 久久久久亚洲av毛片大全| 亚洲视频免费观看视频| 久久久久久久精品吃奶| 久久伊人香网站| 天天躁夜夜躁狠狠躁躁| 老汉色av国产亚洲站长工具| 亚洲人成电影免费在线| 国产aⅴ精品一区二区三区波| 一级黄色大片毛片| 日韩视频一区二区在线观看| 18禁黄网站禁片午夜丰满| 黄片播放在线免费| 又黄又粗又硬又大视频| 久久久久亚洲av毛片大全| 成人亚洲精品一区在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 人妻丰满熟妇av一区二区三区| 精品卡一卡二卡四卡免费| 黑人巨大精品欧美一区二区mp4| 亚洲国产毛片av蜜桃av| 亚洲狠狠婷婷综合久久图片| 午夜老司机福利片| 亚洲在线自拍视频| 777久久人妻少妇嫩草av网站| 免费av毛片视频| 99热国产这里只有精品6| 美国免费a级毛片| av中文乱码字幕在线| 久久久国产成人免费| 99国产精品99久久久久| 精品一品国产午夜福利视频| 久久久久九九精品影院| www国产在线视频色| 免费在线观看影片大全网站| 天堂√8在线中文| 一边摸一边抽搐一进一小说| 国产单亲对白刺激| 人人妻人人添人人爽欧美一区卜| 欧美国产精品va在线观看不卡| 丁香欧美五月| av在线天堂中文字幕 | 最近最新免费中文字幕在线| 麻豆久久精品国产亚洲av | 50天的宝宝边吃奶边哭怎么回事| 久久久水蜜桃国产精品网| 色尼玛亚洲综合影院| 亚洲伊人色综图| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品在线美女| 精品一区二区三区四区五区乱码| 亚洲成人久久性| 极品人妻少妇av视频| 丝袜在线中文字幕| 一级片免费观看大全| 黄网站色视频无遮挡免费观看| 日韩人妻精品一区2区三区| 成人特级黄色片久久久久久久| 亚洲国产欧美一区二区综合| 一级a爱视频在线免费观看| 亚洲av熟女| 欧美人与性动交α欧美软件| 久久99一区二区三区| 欧美乱色亚洲激情| 精品一品国产午夜福利视频| 精品国产乱子伦一区二区三区| av视频免费观看在线观看| 国产成人av激情在线播放| 岛国视频午夜一区免费看| 免费在线观看日本一区| 一进一出好大好爽视频| 桃红色精品国产亚洲av| 成人国产一区最新在线观看| 真人一进一出gif抽搐免费| 亚洲少妇的诱惑av| 法律面前人人平等表现在哪些方面| 日本精品一区二区三区蜜桃| 亚洲一区二区三区不卡视频| 亚洲一区二区三区欧美精品| 亚洲专区中文字幕在线| 久久久久国产一级毛片高清牌| 91麻豆精品激情在线观看国产 | 国产精品影院久久| 久久精品亚洲av国产电影网| 女同久久另类99精品国产91| 午夜福利一区二区在线看| 成人影院久久| 免费少妇av软件| 法律面前人人平等表现在哪些方面| 日韩成人在线观看一区二区三区| 日韩一卡2卡3卡4卡2021年| 国产在线观看jvid| 大码成人一级视频| 在线永久观看黄色视频| 国产麻豆69| 日韩精品青青久久久久久| 麻豆国产av国片精品| 欧美精品啪啪一区二区三区| 美女大奶头视频| 亚洲一区二区三区欧美精品| 久久国产精品人妻蜜桃| 亚洲国产精品sss在线观看 | 夜夜夜夜夜久久久久| 精品少妇一区二区三区视频日本电影| 亚洲av电影在线进入| 正在播放国产对白刺激| 在线观看一区二区三区| 中文欧美无线码| 91精品国产国语对白视频| 国产伦人伦偷精品视频| 成人精品一区二区免费| 亚洲国产精品一区二区三区在线| 亚洲av成人一区二区三| 日韩人妻精品一区2区三区| 麻豆久久精品国产亚洲av | 法律面前人人平等表现在哪些方面| 亚洲欧美一区二区三区久久| 午夜福利欧美成人| 欧美久久黑人一区二区| 亚洲第一青青草原| 男男h啪啪无遮挡| 国产三级黄色录像| 91老司机精品| 国产成年人精品一区二区 | 久久精品国产综合久久久| 黄片小视频在线播放| 黄色怎么调成土黄色| 精品人妻1区二区| √禁漫天堂资源中文www| 手机成人av网站| 成人黄色视频免费在线看| 色哟哟哟哟哟哟| e午夜精品久久久久久久| 精品国产超薄肉色丝袜足j| 又大又爽又粗| 在线观看免费视频日本深夜| 日韩欧美一区视频在线观看| 青草久久国产| 亚洲精品国产区一区二| 国产精品偷伦视频观看了| av中文乱码字幕在线| 国产极品粉嫩免费观看在线| 日韩精品青青久久久久久| aaaaa片日本免费| 亚洲精品一二三| 人人澡人人妻人| 中文字幕另类日韩欧美亚洲嫩草| 高清欧美精品videossex| 人成视频在线观看免费观看| 少妇 在线观看| 啦啦啦免费观看视频1| 91大片在线观看| 丁香欧美五月| 午夜久久久在线观看| 日韩欧美一区二区三区在线观看| 在线观看舔阴道视频| 亚洲九九香蕉| 欧美日韩av久久| 一边摸一边抽搐一进一小说| 国产野战对白在线观看| 中文字幕人妻丝袜制服| 999久久久精品免费观看国产| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av片天天在线观看| 香蕉丝袜av| 免费在线观看影片大全网站| 交换朋友夫妻互换小说| 日韩欧美国产一区二区入口| 级片在线观看| 美女高潮喷水抽搐中文字幕| 超碰成人久久| 国产精品美女特级片免费视频播放器 | 亚洲成av片中文字幕在线观看| 欧美另类亚洲清纯唯美| 午夜福利免费观看在线| 国产精品一区二区在线不卡| 一二三四社区在线视频社区8| 亚洲一区中文字幕在线| 99在线人妻在线中文字幕| 国产乱人伦免费视频| 免费高清视频大片| 亚洲午夜理论影院| 欧美成狂野欧美在线观看| 亚洲专区中文字幕在线| 国产成人影院久久av| 每晚都被弄得嗷嗷叫到高潮| svipshipincom国产片| 男女午夜视频在线观看| 老司机午夜福利在线观看视频| 黄色视频不卡| 欧美乱码精品一区二区三区| 国产精品久久久久久人妻精品电影| 午夜日韩欧美国产| av在线天堂中文字幕 | 亚洲欧美精品综合久久99| 久久精品aⅴ一区二区三区四区| 丰满迷人的少妇在线观看| 亚洲av熟女| 美女国产高潮福利片在线看| 老司机午夜十八禁免费视频| 新久久久久国产一级毛片| 欧美色视频一区免费| 亚洲伊人色综图| 日韩大尺度精品在线看网址 | 亚洲精品一区av在线观看| 国产精品久久久久久人妻精品电影| 少妇 在线观看| 国产单亲对白刺激| 国产一区二区激情短视频| 校园春色视频在线观看| 男女做爰动态图高潮gif福利片 | 免费少妇av软件|