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

    電容儲能型脈沖電源分時分段觸發(fā)策略自動計算方法

    2016-06-29 09:44:36劉旭堃于歆杰劉秀成
    電工技術(shù)學(xué)報 2016年11期

    劉旭堃 于歆杰 劉秀成

    (電力系統(tǒng)國家重點(diǎn)實(shí)驗室(清華大學(xué)) 北京 100084)

    電容儲能型脈沖電源分時分段觸發(fā)策略自動計算方法

    劉旭堃于歆杰劉秀成

    (電力系統(tǒng)國家重點(diǎn)實(shí)驗室(清華大學(xué))北京100084)

    摘要提出一種基于Matlab數(shù)值計算平臺、針對電磁軌道炮發(fā)射過程的快速求解方法,并以此為基礎(chǔ)提出一種兼具準(zhǔn)確性和普適性、針對電容儲能型脈沖電源分時分段觸發(fā)策略的自動計算方法。首先給出電磁軌道炮的電路模型,給出發(fā)射過程中各機(jī)電物理量的表達(dá)式;其次闡述發(fā)射過程求解方法所采用的離散時間步長思想、電源負(fù)載解耦思想和具體實(shí)現(xiàn)方法;進(jìn)而闡述觸發(fā)策略的確定依據(jù)和計算方法的具體實(shí)現(xiàn);最后以一個具有典型代表性、初始儲能為4 MJ的系統(tǒng)為例,通過對比Simplorer仿真結(jié)果,驗證上述方法的可行性和準(zhǔn)確性。結(jié)果表明,自動計算方法得到的觸發(fā)策略能夠使導(dǎo)軌電流波形呈近似平頂波,快速求解方法得到的電樞出膛速度與Simplorer仿真結(jié)果之間的相對誤差低于3%,驗證了該方法的有效性。

    關(guān)鍵詞:觸發(fā)策略發(fā)射過程電容儲能電磁軌道炮

    0引言

    電磁軌道炮是利用電磁發(fā)射技術(shù),以電磁力發(fā)射超高速彈丸,以彈丸動能毀傷目標(biāo)的動能武器系統(tǒng)[1,2],其憑借響應(yīng)快、精確可控等優(yōu)勢成為近年來國內(nèi)外軍事裝備領(lǐng)域的研究熱點(diǎn)[3-5]。電容儲能以發(fā)展成熟、控制簡單等特點(diǎn)成為目前電磁軌道炮脈沖電源的主流儲能形式[6-8]。

    電磁軌道炮在結(jié)構(gòu)和功能上可分為脈沖電源和軌道炮負(fù)載(即發(fā)射裝置)兩部分[9],其結(jié)構(gòu)如圖1所示。脈沖電源的作用是為軌道炮負(fù)載提供短時間、高能量的電流脈沖,其由若干組依次觸發(fā)的脈沖形成單元(Pulse Forming Unit,PFU)并聯(lián)而成。軌道炮負(fù)載的作用是加速電樞,將脈沖電源提供的能量轉(zhuǎn)換為電樞的動能。為避免混淆,特此說明,下述內(nèi)容中“軌道炮負(fù)載”指發(fā)射裝置及其電路模型,“電磁軌道炮”指由脈沖電源和軌道炮負(fù)載構(gòu)成的整體系統(tǒng)。

    圖1 電磁軌道炮的結(jié)構(gòu)Fig.1 Sketch of the electromagnetic railgun

    電磁軌道炮的發(fā)射過程是含有人為控制因素的、復(fù)雜的機(jī)電耦合過程。就軌道炮負(fù)載而言,導(dǎo)軌電流、空間磁場強(qiáng)度等電氣物理量會直接影響電樞加速過程,電樞速度、電樞位移等機(jī)械物理量也會直接影響負(fù)載的電氣特性(如電樞電動勢、導(dǎo)軌電感等)。就脈沖電源而言,所有PFU并聯(lián)向負(fù)載放電,由于存在觸發(fā)先后順序,同一時刻各組PFU的電容電壓和電感電流各不相同,在很大程度上增加了電路的求解難度。Simplorer仿真平臺雖然能夠較為準(zhǔn)確地求解發(fā)射過程,但也存在諸多局限,譬如所有PFU都需要預(yù)先構(gòu)建、所有系統(tǒng)參數(shù)都需要手動設(shè)置、仿真時長(即電樞加速時間)、需要預(yù)先估計上限。因此,基于數(shù)值計算平臺(如Matlab)實(shí)現(xiàn)電磁軌道炮發(fā)射過程的快速求解,是有價值、有意義的。

    對于實(shí)際裝置而言,為了避免軌道燒蝕和電樞過載,導(dǎo)軌電流存在最高閾值約束;為了獲得近似平頂波的導(dǎo)軌電流波形和較高的電樞出膛速度,脈沖電源需要以PFU分時分段觸發(fā)的方式工作[10];為了獲得合理的觸發(fā)策略,以往的方式是通過反復(fù)試驗或仿真,不斷摸索和調(diào)整出大致合理的觸發(fā)策略。但是,這種方式會消耗大量的時間和精力,既繁瑣又低效。因此,實(shí)現(xiàn)脈沖電源分時分段觸發(fā)策略的自動計算,具有很明確、很重要的現(xiàn)實(shí)意義。

    本文的研究工作即圍繞發(fā)射過程的快速求解和觸發(fā)策略的自動計算而展開。

    1系統(tǒng)模型

    1.1PFU的電路模型

    PFU的電路拓?fù)淙鐖D2所示,其中,S為觸發(fā)控制開關(guān),C為儲能電容,L為脈沖形成電感,R為電感電阻與引線電阻之和[11,12]。從脈沖電源側(cè)看,在微小時間段Δt內(nèi),軌道炮負(fù)載可視為恒定電壓源(具體原因?qū)⒃诘?.2節(jié)中論述)。

    圖2 PFU的電路拓?fù)銯ig.2 Topology of the PFU

    根據(jù)Δt開始時刻電容電壓uC0和電感電流iL0的取值,PFU存在以下3種工作狀態(tài):

    1)二階狀態(tài)

    當(dāng)uC0>0、iL0>0時,VD1關(guān)斷,VD2開通,PFU處于RLC二階放電狀態(tài),Δt結(jié)束時刻電容電壓uC和電感電流iL的表達(dá)式分別為

    uC(Δt)=uload+(uC 0-uload)e-αΔtcos(ωΔt)+

    (1)

    e-αΔtsin(ωΔt)

    (2)

    式中,uload為負(fù)載電壓;α和ω為RLC二階電路動態(tài)特性參數(shù),其表達(dá)式為

    (3)

    (4)

    2)一階狀態(tài)

    當(dāng)uC0= 0、iL0>0時,VD1開通,VD2開通,PFU處于RL一階放電狀態(tài),Δt結(jié)束時刻電感電流iL為

    (5)

    式中,τ為RL一階電路動態(tài)特性參數(shù),其表達(dá)式為

    (6)

    3)停止?fàn)顟B(tài)

    當(dāng)uC0= 0、iL0= 0時,VD1關(guān)斷,VD2關(guān)斷,PFU儲能耗盡,放電結(jié)束。

    為敘述簡潔,將PFU電容電壓uC和電感電流iL統(tǒng)稱為PFU狀態(tài)。

    1.2軌道炮負(fù)載的等效電路模型

    軌道炮負(fù)載的等效電路如圖3所示,其中,Lr為導(dǎo)軌電感,Rr為導(dǎo)軌電阻,Rc為導(dǎo)軌與電樞間初始接觸電阻,RVSEC1為導(dǎo)軌接觸速度趨膚效應(yīng)電阻,RVSEC2為電樞接觸速度趨膚效應(yīng)電阻,Ra為電樞電阻,Rp為炮口熄弧電阻,uEMF為電樞動生電動勢。上述電感、電阻和電動勢的表達(dá)式如式(7)~式(16)所示[13-16]。從軌道炮負(fù)載側(cè)看,在微小時間段Δt內(nèi),脈沖電源可視為恒定電流源(具體原因?qū)⒃诘?.2節(jié)中論述)。

    圖3 軌道炮負(fù)載的等效電路Fig.3 Equivalent circuit of the railgun load

    (7)

    (8)

    (9)

    (10)

    (11)

    RVSEC=RVCv1.5

    (12)

    RVSEC1=(1-kVSEC)RVSEC

    (13)

    RVSEC2=kVSECRVSEC

    (14)

    (15)

    uEMF=L′virail

    (16)

    式中,h為導(dǎo)軌高度;b為導(dǎo)軌間距;μ為導(dǎo)軌材料磁導(dǎo)率;L′為導(dǎo)軌電感梯度;ρ10為導(dǎo)軌材料(銅質(zhì))初始電阻率;ρ1為考慮電流擴(kuò)散效應(yīng)的導(dǎo)軌等效電阻率;ρ20為電樞材料(鋁質(zhì))初始電阻率;ρ2為考慮電流擴(kuò)散效應(yīng)的電樞等效電阻率;β為電阻率熱效應(yīng)常數(shù);kc為初始接觸電阻常數(shù);RVC為接觸速度趨膚效應(yīng)電阻常數(shù);kVSEC為接觸速度趨膚效應(yīng)電阻分配比例;irail為導(dǎo)軌電流(其值等于脈沖電源輸出電流isource);ia為電樞電流;t為當(dāng)前時刻;v為電樞速度;x為電樞位移。

    忽略導(dǎo)軌電感Lr,給定脈沖電源輸出電流isource,負(fù)載電壓uload、電樞電流ia、電樞加速度a的表達(dá)式如式(17)~式(19)所示,微小時間段Δt內(nèi)電樞速度增量Δv和電樞位移增量Δx的表達(dá)式如式(20)、式(21)所示。

    (17)

    (18)

    (19)

    Δv=aΔt

    (20)

    Δx=(v0+Δv)Δt

    (21)

    式中,m為電樞質(zhì)量;v0為微小時間段Δt開始時刻電樞速度。

    為敘述簡潔,將導(dǎo)軌電流irail、負(fù)載電壓uload、電樞速度v和電樞位移x統(tǒng)稱為負(fù)載狀態(tài)。

    2發(fā)射過程快速求解方法

    電磁軌道炮發(fā)射過程是非常復(fù)雜的機(jī)電耦合過程,描述這一過程的關(guān)鍵物理量(如導(dǎo)軌電流、電樞速度、電樞位移等)的解析表達(dá)式難以推導(dǎo)。如果希望擺脫仿真平臺限制、基于數(shù)值計算平臺實(shí)現(xiàn)發(fā)射過程的快速求解,引入離散時間步長思想和電源負(fù)載解耦思想是較好的解決方案。

    2.1離散時間步長思想

    將電樞加速時間(即電樞開始加速時刻至最終出膛時刻,毫秒量級)剖分許許多多的微小時間段(微秒量級)。依次計算、串接各微小時間段的PFU狀態(tài)和負(fù)載狀態(tài),即可完成對發(fā)射過程的求解,得到導(dǎo)軌電流、電樞速度、電樞位移等關(guān)鍵物理量的波形及最終的電樞出膛速度和系統(tǒng)發(fā)射效率。

    2.2電源負(fù)載解耦思想

    所有PFU并聯(lián)向負(fù)載放電,由于存在觸發(fā)先后順序,同一時刻各組PFU的狀態(tài)各不相同,直接對電路進(jìn)行求解將非常困難。解耦思想的引入,正是為了實(shí)現(xiàn)微小時間段內(nèi)電路的求解和PFU狀態(tài)、負(fù)載狀態(tài)的計算。

    對于實(shí)際裝置而言,導(dǎo)軌電感感值較小,導(dǎo)軌電流近似平頂波,因此直接忽略導(dǎo)軌電感是可行的(軌道炮負(fù)載隨即可等效為阻值時變的電阻)。進(jìn)而做出如下假設(shè):微小時間段內(nèi)負(fù)載電壓變化幅度很小,近似認(rèn)為其恒定不變。脈沖形成電感的存在,使得各PFU輸出電流不可突變;所有PFU相互并聯(lián)、彼此獨(dú)立,其輸出電流匯聚為導(dǎo)軌電流。因此,在微小時間段內(nèi),導(dǎo)軌電流變化幅度也很小,可以近似認(rèn)為恒定不變?;谏鲜黾僭O(shè)和分析,在微小時間段內(nèi),可以將脈沖電源和軌道炮負(fù)載拆分,分別等效為恒定電流源和恒定電壓源。給定微小時間段開始時刻的PFU狀態(tài)和負(fù)載狀態(tài),根據(jù)第1.1節(jié)和第1.2節(jié)中的電路模型和解析表達(dá)式,即可獲得微小時間段結(jié)束時刻的PFU狀態(tài)和負(fù)載狀態(tài)。

    2.3具體實(shí)現(xiàn)方法

    基于上述思想,將時間剖分、電路解耦,依次計算、串接各微小時間段的PFU狀態(tài)和負(fù)載狀態(tài),即為發(fā)射過程快速求解的實(shí)現(xiàn)方法,其流程圖如圖4所示,具體步驟如下。

    圖4 發(fā)射過程快速求解方法的流程Fig.4 Flow chart of the fast solving method for the launch process

    1)將t時刻負(fù)載電壓代入各組PFU中,根據(jù)t時刻各組PFU的狀態(tài)和式(1)~式(6),求出t+Δt時刻各組PFU的狀態(tài)。

    2)將t+Δt時刻各組PFU的電感電流累加,求出t+Δt時刻脈沖電源輸出電流。

    3)將t+Δt時刻電源輸出電流代入軌道炮負(fù)載中,根據(jù)t時刻負(fù)載狀態(tài)和式(7)~式(21),求出t+Δt時刻負(fù)載狀態(tài)。

    4)重復(fù)步驟1~步驟3,直至電樞出膛。

    上述求解方法是基于“微小時間段內(nèi)負(fù)載電壓不變”的假設(shè)。關(guān)于這一假設(shè)的驗證,將在第4.3節(jié)中給出。

    3觸發(fā)策略自動計算方法

    脈沖電源分時分段觸發(fā)策略的具體含義是指,各組PFU的觸發(fā)單元數(shù)和觸發(fā)時刻。實(shí)現(xiàn)觸發(fā)策略自動計算的具體含義是指,在系統(tǒng)參數(shù)確定的前提下,自動進(jìn)行PFU分時分段設(shè)置,并盡可能獲得較高的電樞出膛速度和系統(tǒng)發(fā)射效率。

    3.1觸發(fā)策略確定依據(jù)

    自動分時分段觸發(fā)策略的確定依據(jù)由以下兩條規(guī)則構(gòu)成:

    1)合適的觸發(fā)單元數(shù)

    能夠使該組PFU觸發(fā)之后的導(dǎo)軌電流峰值在滿足最高閾值約束的前提下達(dá)到最大可能值。

    2)合適的觸發(fā)時刻

    人為給定導(dǎo)軌電流最低閾值約束;導(dǎo)軌電流降至最低閾值的時刻,即作為下一組PFU的觸發(fā)時刻;第1組PFU的觸發(fā)時刻默認(rèn)為0時刻。

    3.2具體實(shí)現(xiàn)方法

    觸發(fā)策略自動計算方法的流程圖如圖5所示,具體步驟如下:

    1)提取第k組PFU觸發(fā)時刻各組PFU狀態(tài)和負(fù)載狀態(tài)。

    2)令第k組PFU的觸發(fā)單元數(shù)從1至nleft(尚未觸發(fā)的PFU的單元數(shù))遍歷。根據(jù)2.3節(jié)給出的PFU狀態(tài)和負(fù)載狀態(tài)計算方法,求解導(dǎo)軌電流波形,尋找導(dǎo)軌電流峰值,根據(jù)3.1節(jié)給出的確定依據(jù),找到合適的觸發(fā)單元數(shù)。

    3)從第k組PFU觸發(fā)之后的導(dǎo)軌電流峰值時刻起,求解導(dǎo)軌電流波形。若當(dāng)前仍有PFU未觸發(fā),則尋找導(dǎo)軌電流幅值降低至最低閾值的時刻或電樞出膛時刻,作為第k+1組PFU的觸發(fā)時刻或發(fā)射過程的結(jié)束時刻;若當(dāng)前已無PFU未觸發(fā),則直接尋找電樞出膛時刻。

    4)重復(fù)步驟1~步驟3,直至電樞出膛。

    圖5 觸發(fā)策略自動計算方法的流程Fig.5 Flow chart of the automatic calculation method for the triggering strategy

    4驗證實(shí)例

    以一個具有典型代表性的、初始儲能為4 MJ的系統(tǒng)為例,對上述方法的可行性和準(zhǔn)確性進(jìn)行驗證。具體方法是將Matlab計算結(jié)果與Simplorer仿真結(jié)果進(jìn)行對比分析。

    需要特別說明的是,Simplorer仿真采用的電路模型和系統(tǒng)參數(shù),已與實(shí)際裝置進(jìn)行了校驗。就電樞出膛速度而言,在不同系統(tǒng)初始儲能下,實(shí)際試驗結(jié)果與仿真結(jié)果之間相對誤差的平均值為2.75%[17-19],說明以Simplorer仿真結(jié)果作為參照對象是合理的。

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

    電磁軌道炮的系統(tǒng)參數(shù)如表1所示。

    表1 電磁軌道炮的系統(tǒng)參數(shù)

    4.2結(jié)果展示

    觸發(fā)策略的計算結(jié)果如表2所示,發(fā)射過程的求解結(jié)果如表3和圖6~圖10所示,其中負(fù)載總等效電阻即負(fù)載電壓與導(dǎo)軌電流的比值。

    分析上述結(jié)果可知,在自動計算方法所得到的觸發(fā)策略的作用下,導(dǎo)軌電流能夠按照預(yù)期在最高閾值和最低閾值之間波動,說明上述觸發(fā)策略計算方法是可行的。

    值得一提的是,給定系統(tǒng)參數(shù),觸發(fā)策略的計算時長約為0.5 s;給定系統(tǒng)參數(shù)和觸發(fā)策略,發(fā)射過程的求解時長約為0.1 s,說明上述計算方法和求解方法是快速的、高效的。

    表2 觸發(fā)策略的計算結(jié)果

    表3 電樞加速時長和電樞出膛速度的計算結(jié)果

    圖6 導(dǎo)軌電流的波形Fig.6 Waveform of the rail current

    圖7 負(fù)載電壓的波形Fig.7 Waveform of the load voltage

    圖8 電樞速度的波形Fig.8 Waveform of the armature velocity

    圖9 電樞位移的波形Fig.9 Waveform of the armature displacement

    圖10 負(fù)載總等效電阻的波形Fig.10 Waveform of the load equivalent resistance

    4.3假設(shè)校驗

    2.3節(jié)中提到,發(fā)射過程快速求解方法是基于“微小時間段內(nèi)負(fù)載電壓不變”的假設(shè)。采用相同的系統(tǒng)模型(即不考慮導(dǎo)軌電感)、系統(tǒng)參數(shù)和觸發(fā)策略進(jìn)行Simplorer仿真,將所得仿真結(jié)果與Matlab計算結(jié)果進(jìn)行對比分析,即可驗證這一假設(shè)的正確性。

    導(dǎo)軌電流的波形對比如圖11所示。Matlab計算結(jié)果和Simplorer仿真結(jié)果之間的吻合度很高,說明“微小時間段內(nèi)負(fù)載電壓不變”的假設(shè)是正確的。

    圖11 導(dǎo)軌電流波形的對比Fig.11 Comparison of the rail current waveforms

    4.4誤差分析

    發(fā)射過程快速求解方法的主要誤差來源是電源負(fù)載解耦時對導(dǎo)軌電感的忽略,這一忽略亦為不得已的舉措。快速求解方法的本質(zhì)是將脈沖電源和軌道炮負(fù)載拆開,分別等效為電流源和電壓源。如果考慮導(dǎo)軌電感,那么在“微小時間段內(nèi)負(fù)載電壓不變”的假設(shè)下,經(jīng)過微小時間段之后,導(dǎo)軌電流會產(chǎn)生很大的跳變,導(dǎo)軌電感電壓也會隨之產(chǎn)生很大的跳變。這將使計算結(jié)果完全錯誤。

    采用相同的系統(tǒng)參數(shù)和觸發(fā)策略,將考慮導(dǎo)軌電感的Simplorer仿真結(jié)果和不考慮導(dǎo)軌電感的Matlab計算結(jié)果進(jìn)行對比分析,即可評估忽略導(dǎo)軌電感所產(chǎn)生的誤差,具體結(jié)果如表4和圖12、圖13所示。

    表4 電樞加速時長和電樞出膛速度的對比

    注:相對誤差的計算,是以Simplorer仿真結(jié)果為真實(shí)值。

    圖12 導(dǎo)軌電流波形的對比Fig.12 Comparison of the rail current waveforms

    圖13 電樞速度波形的對比Fig.13 Comparison of the armature velocity waveforms

    分析上述結(jié)果可知,導(dǎo)軌電感起到了阻礙導(dǎo)軌電流變化的作用,因而使Matlab計算結(jié)果與Simplorer仿真結(jié)果之間出現(xiàn)了一定的誤差。更為重要的是,導(dǎo)軌電感的忽略與否,對電樞出膛速度的影響較小,不足3%。這對于之后將進(jìn)行的、基于上述方法的系統(tǒng)發(fā)射效率研究和系統(tǒng)參數(shù)優(yōu)化而言較為有利。

    5結(jié)論

    本文提出了一種針對電磁軌道炮發(fā)射過程的快速求解方法,并以此為基礎(chǔ)提出了一種電容儲能型脈沖電源分時分段觸發(fā)策略的自動計算方法,同時以一個有典型代表性的系統(tǒng)為例,驗證了上述方法的可行性和準(zhǔn)確性。需要特別強(qiáng)調(diào)的是,實(shí)現(xiàn)發(fā)射過程快速求解和觸發(fā)策略自動計算的更深層次意義在于:基于上述方法能夠進(jìn)行系統(tǒng)發(fā)射效率研究和系統(tǒng)參數(shù)優(yōu)化。

    對實(shí)際裝置而言,影響發(fā)射過程的因素有很多。文中給出的PFU電路模型和軌道炮負(fù)載等效電路模型只是實(shí)際裝置的理論抽象,剔除和忽略了很多次要因素,無法反映所有的物理現(xiàn)象和物理效應(yīng)。但是,本文提出的發(fā)射過程快速求解方法和觸發(fā)策略自動計算方法是具有普適性的,其可行性并不依賴于理論模型的精確性。換言之,更改理論模型和系統(tǒng)參數(shù),上述方法依然能夠適用。

    參考文獻(xiàn)

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

    [2]陳慶國,王永紅,魏新勞,等.電容驅(qū)動型軌道電磁炮電磁過程的計算機(jī)仿真[J].電工技術(shù)學(xué)報,2006,21(4):68-71.

    Chen Qingguo,Wang Yonghong,Wei Xinlao,et al.Computer simulation of electromagnetic process in the capacitor drived rail gun[J].Transactions of China Electrotechnical Society,2006,21(4):68-71.

    [3]Fair H D.Electromagnetic launch science and technology in the United States enters a new era[J].IEEE Transactions on Magnetics,2005,41(1):158-164.

    [4]Thornhill L D,Batteh J H,Littrell D M.Scaling study for the performance of railgun armatures[J].IEEE Transactions on Plasma Science,1989,17(3):409-421.

    [5]李軍,嚴(yán)萍,袁偉群.電磁軌道炮發(fā)射技術(shù)的發(fā)展與現(xiàn)狀[J].高電壓技術(shù),2014,40(4):1052-1064.

    Li Jun,Yan Ping,Yuan Weiqun.Electromagnetic gun technology and its development[J].High Voltage Engineering,2014,40(4):1052-1064.

    [6]馬山剛,于歆杰,李臻.用于電磁發(fā)射的電感儲能型脈沖電源的研究現(xiàn)狀綜述[J].電工技術(shù)學(xué)報,2015,30(24):222-228,236.

    Ma Shangang,Yu Xinjie,Li Zhen.A review of the current research situation of inductive pulsed-power supplies for electromagnetic launch[J].Transactions of China Electrotechnical Society,2015,30(24):222-228,236.

    [7]龔晨,于歆杰,劉秀成.電容儲能型軌道炮連續(xù)發(fā)射系統(tǒng)設(shè)計與仿真[J].電工技術(shù)學(xué)報,2013,28(增2):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.

    [8]李化,章妙,林福昌,等.金屬化膜電容器自愈理論及規(guī)律研究[J].電工技術(shù)學(xué)報,2012,27(9):218-223,230.

    Li Hua,Zhang Miao,Lin Fuchang,et al.Study on theory and influence factors of self-healing in metallized film capacitors[J].Transactions of China Electrotechnical Society,2012,27(9):218-223,230.

    [9]Gong Chen,Yu Xinjie,Liu Xiucheng.Study on the system efficiency of the capacitive pulsed-power supply[J].IEEE Transactions on Plasma Science,2015,43(5):1441-1447.

    [10]范昭楠,于歆杰.基于過程集成的電磁軌道炮脈沖電源多目標(biāo)優(yōu)化[J].電工技術(shù)學(xué)報,2010,25(5):20-24.

    Fan Zhaonan,Yu Xinjie.Process-integration based multi-objective optimization for pulsed power supply of electromagnetic guns[J].Transactions of China Electrotechnical Society,2010,25(5):20-24.

    [11]McFarland J,McNab I R.A long-range naval railgun[J].IEEE Transactions on Magnetics,2003,39(1):289-294.

    [12]Fair H D.Electric launch science and technology in the United States[J].IEEE Transactions on Magnetics,2003,39(1):11-17.

    [13]Deadrick J F,Hawke S R,Scudder D J.MAGRAC-a railgun simulation program[J].IEEE Transactions on Magnetics,1982,18(1):94-104.

    [14]Yu Xinjie,F(xiàn)an Zhaonan.Simulation and two-objective optimization of the electromagnetic-railgun model considering VSEC resistance and contact resistance[J].IEEE Transactions on Plasma Science,2011,39(1):405-410.

    [15]Hsieh K T,Satapathy S,Hsieh M T.Effects of pressure-dependent contact resistivity on contact interfacial conditions[J].IEEE Transactions on Magnetics,2009,45(1):313-318.

    [16]Engel T G,Neri J M,Veracka M J.Characterization of the velocity skin effect in the surface layer of a railgun sliding contact[J].IEEE Transactions on Magnetics,2008,44(7):1837-1844.

    [17]史正軍.高功率脈沖電源的仿真與優(yōu)化設(shè)計[D].北京:清華大學(xué),2008.

    [18]范昭楠.電磁軌道炮仿真系統(tǒng)的驗證與雙目標(biāo)優(yōu)化[D].北京:清華大學(xué),2010.

    [19]龔晨.電容儲能型電磁軌道炮系統(tǒng)效率的研究[D].北京:清華大學(xué),2014.

    An Automatic Calculation Method for the Triggering Strategy of the Capacitive Pulsed-Power Supply

    Liu XukunYu XinjieLiu Xiucheng

    (State Key Laboratory of Power System Department of Electrical Engineering Tsinghua UniversityBeijing100084China)

    AbstractBased on Matlab platform,a fast solving method for the launch process of the electromagnetic railgun is proposed.By using this method,an automatic calculation method for the triggering strategy of the capacitive pulsed-power supply is then proposed,which possesses the advantages of both accuracy and applicability.The circuit model of the electromagnetic railgun system is firstly analyzed and the mathematical expressions of the relevant electrical and mechanical quantities are also given.Then the time discretization idea,the circuit decoupling idea and the specific implementation of the fast solving method are elaborated in detail.And the basic principle and the specific implementation of the automatic calculation method are also discussed.Finally,a typically representative system with 4 MJ initial energy is given.The feasibility and accuracy of these two methods are proven through the comparison with the Simplorer simulation results which are confirmed to be compliant with the real system.Results show that the rail current waveform is approximate to a flat wave under the action of the triggering strategy acquired from the automatic calculation method,and the relative error of the armature launch velocity between the Simplorer simulation result and the Matlab calculation result is within 3%.These results confirm the feasibility and accuracy of the above methods.

    Keywords:Triggering strategy,launch process,capacitive energy-storage,electromagnetic railgun

    收稿日期2015-04-20改稿日期2015-06-29

    作者簡介E-mail:liuxukuncn@126.com E-mail:yuxj@tsinghua.edu.cn(通信作者)

    中圖分類號:TM133

    國家自然科學(xué)基金(50877039)和清華大學(xué)自主科研計劃(20121087927)資助項目。

    劉旭堃男,1993年生,博士研究生,研究方向為脈沖電源和電磁發(fā)射技術(shù)。

    于歆杰男,1973年生,副教授,博士生導(dǎo)師,研究方向為電磁發(fā)射技術(shù)、無線能量傳輸技術(shù)、計算電磁學(xué)。

    久久99精品国语久久久| 亚洲精品第二区| 成人毛片60女人毛片免费| 亚洲精品一区蜜桃| 亚洲高清免费不卡视频| 69精品国产乱码久久久| 亚洲精品国产av蜜桃| 青青草视频在线视频观看| 免费观看a级毛片全部| 夫妻性生交免费视频一级片| 高清毛片免费看| 精品少妇久久久久久888优播| 成人手机av| 午夜影院在线不卡| 一二三四中文在线观看免费高清| 亚洲精品久久久久久婷婷小说| 97超碰精品成人国产| 婷婷色av中文字幕| 美女福利国产在线| 草草在线视频免费看| 人体艺术视频欧美日本| 香蕉精品网在线| 男女边摸边吃奶| 午夜免费男女啪啪视频观看| 蜜桃久久精品国产亚洲av| 丰满少妇做爰视频| 十分钟在线观看高清视频www| 三级国产精品欧美在线观看| 在线免费观看不下载黄p国产| 久久精品久久精品一区二区三区| 亚洲综合色惰| 男女啪啪激烈高潮av片| 国产伦理片在线播放av一区| 简卡轻食公司| 欧美另类一区| 99热国产这里只有精品6| 五月伊人婷婷丁香| 欧美日韩精品成人综合77777| 制服诱惑二区| 亚洲美女黄色视频免费看| 大香蕉久久成人网| 精品人妻偷拍中文字幕| 熟女av电影| 夜夜爽夜夜爽视频| 久久99热6这里只有精品| 黄片播放在线免费| 日韩人妻高清精品专区| kizo精华| 久久国产亚洲av麻豆专区| 久久久久久久久大av| 国产有黄有色有爽视频| 蜜桃在线观看..| 蜜桃久久精品国产亚洲av| 亚洲欧洲国产日韩| 久久久久久久久久久免费av| 国产亚洲av片在线观看秒播厂| 国产成人精品婷婷| 成人亚洲精品一区在线观看| 精品一区二区三区视频在线| 日本猛色少妇xxxxx猛交久久| 日本av免费视频播放| 男人添女人高潮全过程视频| 久久99热6这里只有精品| 王馨瑶露胸无遮挡在线观看| 婷婷色麻豆天堂久久| 日韩电影二区| 免费黄色在线免费观看| 美女xxoo啪啪120秒动态图| 欧美亚洲 丝袜 人妻 在线| 亚洲精品久久成人aⅴ小说 | 国产免费又黄又爽又色| 日韩欧美精品免费久久| 一本色道久久久久久精品综合| 伊人亚洲综合成人网| 高清毛片免费看| 51国产日韩欧美| 日韩欧美精品免费久久| 国产日韩一区二区三区精品不卡 | 欧美xxxx性猛交bbbb| 丰满乱子伦码专区| 亚洲av男天堂| 国产黄频视频在线观看| 欧美最新免费一区二区三区| 91成人精品电影| 亚洲国产精品成人久久小说| 久久人人爽av亚洲精品天堂| 日日撸夜夜添| 欧美日韩成人在线一区二区| 看免费成人av毛片| 精品久久久噜噜| 欧美 日韩 精品 国产| 亚洲精品久久午夜乱码| 97超视频在线观看视频| 九草在线视频观看| 少妇被粗大的猛进出69影院 | 精品人妻在线不人妻| 18禁动态无遮挡网站| 亚洲国产最新在线播放| 国产精品欧美亚洲77777| 亚洲欧洲国产日韩| 伊人亚洲综合成人网| 欧美激情国产日韩精品一区| 自线自在国产av| 亚洲欧美清纯卡通| 9色porny在线观看| 18禁裸乳无遮挡动漫免费视频| 最黄视频免费看| 亚洲av欧美aⅴ国产| 女人久久www免费人成看片| 天天影视国产精品| 午夜91福利影院| 校园人妻丝袜中文字幕| 18禁裸乳无遮挡动漫免费视频| 人妻夜夜爽99麻豆av| 不卡视频在线观看欧美| 久久毛片免费看一区二区三区| 午夜精品国产一区二区电影| 亚洲欧美一区二区三区国产| av卡一久久| 777米奇影视久久| 亚洲欧洲国产日韩| 欧美日韩成人在线一区二区| 插逼视频在线观看| 91精品一卡2卡3卡4卡| 蜜桃久久精品国产亚洲av| 男人爽女人下面视频在线观看| 久久97久久精品| 三级国产精品片| 国产熟女欧美一区二区| 夫妻性生交免费视频一级片| 精品国产露脸久久av麻豆| 美女脱内裤让男人舔精品视频| 日本黄大片高清| 国产精品无大码| 26uuu在线亚洲综合色| 国产精品欧美亚洲77777| 欧美精品国产亚洲| 极品人妻少妇av视频| 国产一区亚洲一区在线观看| 亚洲av男天堂| av一本久久久久| 综合色丁香网| 国产精品久久久久久av不卡| 国产精品久久久久久精品电影小说| 三级国产精品片| 久久精品国产亚洲av天美| 婷婷色av中文字幕| 亚洲美女搞黄在线观看| 精品人妻一区二区三区麻豆| 免费少妇av软件| 欧美激情国产日韩精品一区| 亚洲经典国产精华液单| 大话2 男鬼变身卡| 免费人妻精品一区二区三区视频| 搡女人真爽免费视频火全软件| 国产精品人妻久久久影院| 精品少妇久久久久久888优播| 中文乱码字字幕精品一区二区三区| 亚洲国产欧美日韩在线播放| 免费少妇av软件| 久久久国产精品麻豆| 精品人妻偷拍中文字幕| 亚洲欧洲国产日韩| av天堂久久9| 国产免费一区二区三区四区乱码| 最近的中文字幕免费完整| 久久精品国产a三级三级三级| kizo精华| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品自拍成人| 亚洲国产av影院在线观看| 伦理电影免费视频| 99九九在线精品视频| 久久人人爽人人片av| av在线老鸭窝| 人妻一区二区av| 欧美三级亚洲精品| 91国产中文字幕| 各种免费的搞黄视频| av国产精品久久久久影院| 九色亚洲精品在线播放| 97在线人人人人妻| 欧美激情 高清一区二区三区| 三级国产精品欧美在线观看| 少妇被粗大猛烈的视频| 成人国产av品久久久| 少妇人妻精品综合一区二区| 亚洲国产精品国产精品| 国产精品不卡视频一区二区| 夫妻性生交免费视频一级片| 亚洲美女黄色视频免费看| 一本色道久久久久久精品综合| 在线观看国产h片| 色婷婷久久久亚洲欧美| 久久久久久人妻| 在线播放无遮挡| 狠狠婷婷综合久久久久久88av| 如日韩欧美国产精品一区二区三区 | 色5月婷婷丁香| 在线看a的网站| 精品久久久久久电影网| 欧美少妇被猛烈插入视频| 国产在视频线精品| 91精品国产九色| 久久精品久久久久久久性| 国产亚洲av片在线观看秒播厂| 人妻制服诱惑在线中文字幕| 在线观看免费日韩欧美大片 | 国产一级毛片在线| 日韩视频在线欧美| 免费观看在线日韩| 国产熟女午夜一区二区三区 | 国产淫语在线视频| 亚洲成人av在线免费| 街头女战士在线观看网站| 国产男女内射视频| 人人妻人人添人人爽欧美一区卜| 美女国产高潮福利片在线看| videos熟女内射| 看十八女毛片水多多多| 男女边摸边吃奶| 欧美精品国产亚洲| 欧美 亚洲 国产 日韩一| 亚洲av在线观看美女高潮| 国产在线视频一区二区| 蜜臀久久99精品久久宅男| 一级片'在线观看视频| 午夜福利影视在线免费观看| 中文字幕制服av| 另类亚洲欧美激情| 亚洲精品aⅴ在线观看| 国产伦精品一区二区三区视频9| 亚洲色图综合在线观看| 亚洲国产欧美日韩在线播放| 十八禁网站网址无遮挡| 久久久午夜欧美精品| 一级片'在线观看视频| 一本一本综合久久| 免费少妇av软件| 男女免费视频国产| 在线观看三级黄色| 一本久久精品| 伊人久久国产一区二区| 五月开心婷婷网| 国产熟女午夜一区二区三区 | 精品视频人人做人人爽| tube8黄色片| 夫妻性生交免费视频一级片| 国产伦精品一区二区三区视频9| 免费观看无遮挡的男女| 妹子高潮喷水视频| 亚洲国产精品999| 国产精品无大码| 人人澡人人妻人| 久久久精品区二区三区| 成人午夜精彩视频在线观看| 日本-黄色视频高清免费观看| 十八禁网站网址无遮挡| 国产有黄有色有爽视频| 亚洲欧美成人综合另类久久久| 国产综合精华液| 性色av一级| 亚洲,欧美,日韩| 超色免费av| 丰满迷人的少妇在线观看| 热re99久久精品国产66热6| 国产精品不卡视频一区二区| 超色免费av| 亚洲情色 制服丝袜| 热re99久久精品国产66热6| 国产精品99久久久久久久久| 制服诱惑二区| 欧美日韩视频高清一区二区三区二| 欧美xxxx性猛交bbbb| 赤兔流量卡办理| 只有这里有精品99| 最近中文字幕高清免费大全6| 久热久热在线精品观看| 国产探花极品一区二区| 美女主播在线视频| 日日爽夜夜爽网站| 亚洲av在线观看美女高潮| 99国产综合亚洲精品| 我的女老师完整版在线观看| 久久久久网色| 夜夜爽夜夜爽视频| 中文天堂在线官网| 久久ye,这里只有精品| a级毛色黄片| 色5月婷婷丁香| av免费观看日本| 男男h啪啪无遮挡| 高清毛片免费看| 夜夜看夜夜爽夜夜摸| 自拍欧美九色日韩亚洲蝌蚪91| 狂野欧美白嫩少妇大欣赏| 国产精品国产av在线观看| 伊人亚洲综合成人网| 91久久精品国产一区二区三区| 三级国产精品片| 色5月婷婷丁香| 亚洲熟女精品中文字幕| 啦啦啦中文免费视频观看日本| 天堂俺去俺来也www色官网| 久久精品国产亚洲av涩爱| 精品国产一区二区久久| 18禁在线无遮挡免费观看视频| 自线自在国产av| 人体艺术视频欧美日本| 99热网站在线观看| 一级毛片我不卡| 国产深夜福利视频在线观看| 肉色欧美久久久久久久蜜桃| 一级片'在线观看视频| 亚洲精品第二区| 亚洲国产日韩一区二区| 日韩精品有码人妻一区| 国模一区二区三区四区视频| 久久久久久久亚洲中文字幕| 人妻夜夜爽99麻豆av| av天堂久久9| 久久久久久久国产电影| 你懂的网址亚洲精品在线观看| 国产亚洲欧美精品永久| 97精品久久久久久久久久精品| 97超碰精品成人国产| 男人操女人黄网站| 午夜福利视频精品| av在线app专区| 99re6热这里在线精品视频| 22中文网久久字幕| 日韩强制内射视频| 妹子高潮喷水视频| 黑人巨大精品欧美一区二区蜜桃 | videos熟女内射| 两个人的视频大全免费| 久久午夜福利片| 精品一区二区免费观看| 少妇的逼好多水| 久久 成人 亚洲| av线在线观看网站| 成年美女黄网站色视频大全免费 | 久久久国产欧美日韩av| 天堂中文最新版在线下载| 又粗又硬又长又爽又黄的视频| 国产 精品1| 色5月婷婷丁香| 国产熟女午夜一区二区三区 | 大陆偷拍与自拍| 激情五月婷婷亚洲| 999精品在线视频| 亚洲成人av在线免费| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲五月色婷婷综合| 欧美激情国产日韩精品一区| 国产精品秋霞免费鲁丝片| 国产老妇伦熟女老妇高清| 久久婷婷青草| 欧美亚洲 丝袜 人妻 在线| 免费观看性生交大片5| 久久精品国产亚洲av天美| 日本-黄色视频高清免费观看| 亚洲精品色激情综合| 少妇 在线观看| 黄色配什么色好看| 我的女老师完整版在线观看| 人妻人人澡人人爽人人| 国产淫语在线视频| 国产精品一国产av| videossex国产| 在线观看免费视频网站a站| av在线观看视频网站免费| 午夜福利视频精品| 国产日韩欧美在线精品| 国产精品一国产av| 国产高清国产精品国产三级| 久久精品久久久久久久性| 国产欧美亚洲国产| 男人操女人黄网站| 国产在视频线精品| 亚洲精品久久成人aⅴ小说 | 一区二区av电影网| 99久久精品国产国产毛片| 成年女人在线观看亚洲视频| av播播在线观看一区| 午夜福利影视在线免费观看| 亚洲国产精品999| 精品久久久久久电影网| 久久久久国产精品人妻一区二区| 97超碰精品成人国产| 亚洲美女视频黄频| a 毛片基地| 黄色配什么色好看| 精品久久蜜臀av无| 在线观看www视频免费| 亚洲成人av在线免费| 2018国产大陆天天弄谢| 免费黄色在线免费观看| 亚洲国产日韩一区二区| 中国三级夫妇交换| 狂野欧美激情性bbbbbb| 伦理电影大哥的女人| 中国三级夫妇交换| 久久精品国产亚洲网站| 99视频精品全部免费 在线| 青春草亚洲视频在线观看| 老熟女久久久| 日韩av在线免费看完整版不卡| 久久精品久久精品一区二区三区| 国产色婷婷99| 亚洲精品亚洲一区二区| 日日爽夜夜爽网站| 中文字幕久久专区| 男人操女人黄网站| 在线 av 中文字幕| 久久精品国产亚洲av涩爱| 免费不卡的大黄色大毛片视频在线观看| 99热全是精品| 夫妻性生交免费视频一级片| 激情五月婷婷亚洲| 新久久久久国产一级毛片| 国产精品国产三级专区第一集| 午夜福利影视在线免费观看| 国产av一区二区精品久久| 少妇丰满av| 人妻少妇偷人精品九色| 国产精品一区二区在线观看99| 啦啦啦在线观看免费高清www| 亚洲av不卡在线观看| 亚洲av日韩在线播放| 一区二区三区四区激情视频| 日韩精品免费视频一区二区三区 | 日韩成人伦理影院| 色婷婷av一区二区三区视频| 国产精品国产三级国产专区5o| 九九在线视频观看精品| 丰满乱子伦码专区| 新久久久久国产一级毛片| 伊人亚洲综合成人网| 少妇猛男粗大的猛烈进出视频| 99久久精品一区二区三区| 亚洲av成人精品一二三区| 久久人人爽人人片av| 美女脱内裤让男人舔精品视频| 18禁在线播放成人免费| 黄色欧美视频在线观看| 国产欧美日韩综合在线一区二区| 久久久久网色| 日韩熟女老妇一区二区性免费视频| 精品一区在线观看国产| 免费黄频网站在线观看国产| 国产精品无大码| 亚洲高清免费不卡视频| 久久亚洲国产成人精品v| 亚洲怡红院男人天堂| 一级黄片播放器| 你懂的网址亚洲精品在线观看| 一二三四中文在线观看免费高清| 一区二区日韩欧美中文字幕 | 伊人久久精品亚洲午夜| 97超视频在线观看视频| 超碰97精品在线观看| 天天躁夜夜躁狠狠久久av| 日韩欧美精品免费久久| 黄色配什么色好看| 国产亚洲午夜精品一区二区久久| 中文字幕免费在线视频6| 亚洲av欧美aⅴ国产| 国产成人免费无遮挡视频| 18+在线观看网站| 97超视频在线观看视频| 免费黄网站久久成人精品| 亚洲精品久久久久久婷婷小说| 最近手机中文字幕大全| 国产又色又爽无遮挡免| 日本欧美视频一区| 免费高清在线观看视频在线观看| 18禁在线播放成人免费| 日韩一区二区视频免费看| 我的女老师完整版在线观看| 欧美成人精品欧美一级黄| 亚洲一区二区三区欧美精品| 春色校园在线视频观看| 黄色一级大片看看| 欧美性感艳星| 成人综合一区亚洲| 菩萨蛮人人尽说江南好唐韦庄| 狂野欧美激情性bbbbbb| 我的老师免费观看完整版| 性色avwww在线观看| 国产精品一区二区三区四区免费观看| 亚洲精品国产色婷婷电影| 男人添女人高潮全过程视频| 桃花免费在线播放| 欧美日韩亚洲高清精品| av网站免费在线观看视频| 国产精品久久久久成人av| 国产在线视频一区二区| 99九九线精品视频在线观看视频| 成人午夜精彩视频在线观看| 亚州av有码| 亚洲精品一区蜜桃| 精品国产乱码久久久久久小说| 韩国高清视频一区二区三区| 男女啪啪激烈高潮av片| 日韩亚洲欧美综合| 99精国产麻豆久久婷婷| 美女内射精品一级片tv| 亚洲精品美女久久av网站| 中文字幕最新亚洲高清| 亚洲精品一二三| 国产精品国产三级国产av玫瑰| 高清在线视频一区二区三区| 少妇被粗大猛烈的视频| 妹子高潮喷水视频| 日韩av在线免费看完整版不卡| 女性被躁到高潮视频| 亚洲成色77777| 在线观看免费高清a一片| 欧美激情极品国产一区二区三区 | 国产一区亚洲一区在线观看| 内地一区二区视频在线| 婷婷色综合大香蕉| 热re99久久精品国产66热6| 久久久午夜欧美精品| 亚洲国产色片| 国产av精品麻豆| 丝袜喷水一区| 99热网站在线观看| 亚洲图色成人| 成人国产麻豆网| 中文字幕制服av| 又大又黄又爽视频免费| 丝袜在线中文字幕| 日韩一本色道免费dvd| 成年人午夜在线观看视频| 免费观看在线日韩| 婷婷色综合www| .国产精品久久| videossex国产| 夫妻性生交免费视频一级片| 亚洲精品av麻豆狂野| 久久鲁丝午夜福利片| 亚洲情色 制服丝袜| 午夜福利,免费看| 天美传媒精品一区二区| 国产探花极品一区二区| 午夜免费鲁丝| 欧美+日韩+精品| 欧美性感艳星| 国产伦理片在线播放av一区| 99九九线精品视频在线观看视频| 日日摸夜夜添夜夜添av毛片| 一边亲一边摸免费视频| 亚洲三级黄色毛片| 2018国产大陆天天弄谢| 国产淫语在线视频| 国产男女内射视频| 日韩欧美一区视频在线观看| av一本久久久久| 91精品一卡2卡3卡4卡| 国产亚洲最大av| 免费看不卡的av| 国产日韩一区二区三区精品不卡 | 日韩av不卡免费在线播放| 丰满少妇做爰视频| 人体艺术视频欧美日本| 久久国产亚洲av麻豆专区| 亚洲精品美女久久av网站| 高清不卡的av网站| 黄色配什么色好看| 好男人视频免费观看在线| 久久国产精品大桥未久av| 人妻人人澡人人爽人人| 一级毛片黄色毛片免费观看视频| 日本欧美国产在线视频| 成人毛片a级毛片在线播放| 丝袜脚勾引网站| 这个男人来自地球电影免费观看 | 久久精品国产a三级三级三级| 中国三级夫妇交换| 日韩电影二区| 丰满迷人的少妇在线观看| 一区二区三区乱码不卡18| 久久久久久久久久成人| 亚洲成人手机| 国产免费福利视频在线观看| 男人爽女人下面视频在线观看| 欧美日韩在线观看h| 丁香六月天网| a级毛片免费高清观看在线播放| 免费看光身美女| 夜夜骑夜夜射夜夜干| 午夜福利网站1000一区二区三区| a级毛色黄片| 少妇人妻久久综合中文| 国产淫语在线视频| 免费大片18禁| 九九在线视频观看精品| 99热6这里只有精品| 看非洲黑人一级黄片| 国产精品麻豆人妻色哟哟久久| 久久久久人妻精品一区果冻| av有码第一页| 国产成人aa在线观看| 亚洲色图 男人天堂 中文字幕 | 91精品伊人久久大香线蕉| 国产成人精品无人区| 免费看av在线观看网站| www.色视频.com| 免费人成在线观看视频色| 一级毛片电影观看| 人妻制服诱惑在线中文字幕| 美女脱内裤让男人舔精品视频| 久久国产亚洲av麻豆专区|