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

    點狀熱干擾源誘導(dǎo)的準(zhǔn)各向同性高溫超導(dǎo)股線失超特性

    2024-02-05 09:33:06王起悅劉子秋孫梓源
    電工技術(shù)學(xué)報 2024年2期
    關(guān)鍵詞:股線傳播速度帶材

    王起悅 劉子秋 孫梓源 楊 宇 皮 偉

    點狀熱干擾源誘導(dǎo)的準(zhǔn)各向同性高溫超導(dǎo)股線失超特性

    王起悅1劉子秋2孫梓源1楊 宇1皮 偉1

    (1. 華北電力大學(xué)新能源電力系統(tǒng)國家重點實驗室 北京 102206 2. 國網(wǎng)山東省電力公司濟南供電公司 濟南 250000)

    第二代高溫超導(dǎo)帶材因其較高的臨界電流密度以及良好的機械特性,近年來在超導(dǎo)研究領(lǐng)域得到極大關(guān)注。由REBCO帶材組成的高溫超導(dǎo)導(dǎo)體在液氮溫度下的零電阻特性使其在超導(dǎo)輸電領(lǐng)域具有巨大的潛力。從安全角度看,高溫超導(dǎo)導(dǎo)體的失超保護研究是其獲得廣泛應(yīng)用的一個關(guān)鍵。該文建立了準(zhǔn)各向同性超導(dǎo)股線的三維失超仿真模型,耦合了電場、磁場和溫度場,模擬了不同運行電流以及不同點狀熱干擾能量下超導(dǎo)股線的失超行為,求解出股線的溫度、電流密度以及磁場分布,得到了準(zhǔn)各向同性高溫超導(dǎo)股線的最小失超能(MQE)以及失超傳播速度(QPV),對超導(dǎo)裝置的失超保護具有重要意義。該文對準(zhǔn)各向同性超導(dǎo)股線建立的三維失超模型,也可以應(yīng)用到其他結(jié)構(gòu)的高溫超導(dǎo)導(dǎo)體的失超研究中。

    準(zhǔn)各向同性超導(dǎo)股線 失超行為 最小失超能 失超傳播速度

    0 引言

    近年來,隨著超導(dǎo)電力技術(shù)的研究和發(fā)展[1-6],以REBCO為載流層的第二代高溫超導(dǎo)帶材[7]因其載流能力強、交流損耗低,在超導(dǎo)電力領(lǐng)域獲得了極大關(guān)注,被廣泛應(yīng)用于超導(dǎo)電纜、超導(dǎo)限流器、超導(dǎo)變壓器等超導(dǎo)電力設(shè)備中?;诘诙邷爻瑢?dǎo)帶材的超導(dǎo)電纜導(dǎo)體,例如RACC(roebel assembled coated conductor)導(dǎo)體/電纜[8-10]、CORC(conductor on round core)導(dǎo)體/電纜[11-12],TSTC(twisted stacked-tapes cable)導(dǎo)體/電纜[13]以及準(zhǔn)各向同性超導(dǎo)股線(Quasi-Isotropic Superconducting, Q-IS)導(dǎo)體/電纜[14]等被相繼提出,其中,Q-IS采用REBCO二代高溫超導(dǎo)帶材,在制作時將四股堆疊帶材(每股都是正方形截面)組裝起來,且每股堆疊帶材都相對于相鄰股旋轉(zhuǎn)90°,外覆金屬填充物和金屬包套,這種結(jié)構(gòu)使Q-IS具有通流密度大、受外場影響小的優(yōu)點,在超導(dǎo)電力應(yīng)用中具有廣闊的前景。

    高溫超導(dǎo)裝置在電網(wǎng)中運行時,易受到電力系統(tǒng)故障,如短路故障、電流不平衡的影響,嚴重時可能會使超導(dǎo)導(dǎo)體因大電流或機械應(yīng)力而發(fā)生失超行為。詳細來講,故障處超導(dǎo)體會由超導(dǎo)態(tài)轉(zhuǎn)變?yōu)檎B(tài),發(fā)生局部失超,電流從超導(dǎo)區(qū)域分流到金屬區(qū)域,由于高溫超導(dǎo)體的傳輸電流很高,此時在金屬區(qū)域產(chǎn)生大量焦耳熱,使得超導(dǎo)體溫度持續(xù)上升,當(dāng)焦耳熱的積累量高于導(dǎo)體與液氮熱交換導(dǎo)出的熱量時,失超會繼續(xù)擴散,不僅有可能損壞高溫超導(dǎo)裝置,甚至影響整個電力系統(tǒng)的穩(wěn)定運行,因此超導(dǎo)體的失超特性是目前超導(dǎo)裝置實用化過程中需要重點解決的關(guān)鍵技術(shù)之一。

    目前單根超導(dǎo)帶材失超特性已被大量研究,例如,W. K. Chan等通過建立由多個高縱橫比薄層組成的數(shù)值混合電熱模型(2D-3D)來模擬超導(dǎo)帶材的失超行為[15-16],F(xiàn). Roy等通過建立三維磁熱模型對不同基底的高溫超導(dǎo)帶材的失超特性進行了詳細的研究[17]?;诘诙邷爻瑢?dǎo)帶材的超導(dǎo)導(dǎo)體的失超特性也被研究,如Wang Yawei等建立了基于T-A的三維磁熱模型對CORC導(dǎo)體在不同接觸電阻下的失超特性進行了研究[18]。Kang Rui等建立了熱-電-力模型,得到了影響液氦溫度下TSTC導(dǎo)體的失超特性的因素[19]。

    本文將通過建立三維電磁熱模型對液氮溫度下的準(zhǔn)各向同性超導(dǎo)股線(Q-IS)的失超特性進行研究,得到準(zhǔn)各向同性超導(dǎo)股線在自場下通以不同運行電流以及施加不同熱干擾的磁場、電流密度和溫度分布,進而計算出該股線的最小失超能[20]以及失超傳播速度。

    1 準(zhǔn)各向同性超導(dǎo)股線失超特性模型的建立

    1.1 準(zhǔn)各向同性超導(dǎo)股線結(jié)構(gòu)

    準(zhǔn)各向同性超導(dǎo)股線結(jié)構(gòu)如圖1所示,由超導(dǎo)線芯、鋁箔、鋁填充、銅包套構(gòu)成,其中超導(dǎo)線芯是由四股超導(dǎo)帶材依次旋轉(zhuǎn)90°堆疊而成,超導(dǎo)線芯外側(cè)包裹鋁箔、鋁填充以及銅包套。

    (a)準(zhǔn)各向同性超導(dǎo)股線結(jié)構(gòu)

    (b)股線橫截面圖

    (c)超導(dǎo)帶材放大圖

    圖1 準(zhǔn)各向同性股線示意圖

    Fig.1 Schematic diagram of the quasi-isotropic superconducting strand

    股線的主要參數(shù)見表1,REBCO超導(dǎo)帶材由上海超導(dǎo)有限公司提供,忽略過渡層的超導(dǎo)帶材主要尺寸見表2。

    表1 股線的主要參數(shù)

    Tab.1 Main parameters of the cable

    表2 超導(dǎo)帶材的主要尺寸

    Tab.2 Main parameters of the tape

    1.2 失超模型的建立

    因股線具有對稱性,為減少計算量,根據(jù)對稱性僅對股線的一半(7 cm)進行仿真,將上述準(zhǔn)各向同性股線在有限元分析軟件中建立幾何模型如圖2所示,尺寸見表1,其中股線外部表示液氮浴。對該幾何模型進行網(wǎng)格劃分,劃分結(jié)果如圖3所示。

    圖2 準(zhǔn)各向同性股線幾何模型的建立

    為實現(xiàn)股線電磁場與溫度場的耦合,采用方法[21-22]求解股線的磁場,在有限元軟件中通過PDE模塊與固體傳熱模塊的耦合實現(xiàn)股線電磁熱模型的求解[23]。

    圖3 準(zhǔn)各向同性股線幾何模型的網(wǎng)格劃分

    1.2.1 電磁場仿真模塊的建立

    基于法求解股線的磁場部分,自變量為磁場強度變量,其控制方程為

    式中,0為真空中的磁導(dǎo)率;r為相對磁導(dǎo)率。

    對于三維問題,上述控制方程可以寫為分量形式,即

    第二代超導(dǎo)體的電阻率滿足-指數(shù)規(guī)律,即

    式中,C0為臨界電場強度(10-4V/m);為反映曲線陡度的指數(shù),在本模型中取值為21;C為臨界電流密度,其值的大小受磁場強度以及溫度的影響[24],即

    式中,C為臨界轉(zhuǎn)變溫度,本文取90 K;C0為80根單根超導(dǎo)帶材在自場下77 K時的臨界電流疊加值,即

    式中,為超導(dǎo)帶材的根數(shù)(80);s77K為單根帶材在77 K溫度下的臨界電流(90 A)。

    由于本文模擬股線在零場液氮冷卻條件下,因此本模型的磁場初始條件為

    當(dāng)準(zhǔn)各向同性超導(dǎo)股線通電流時,對股線橫截面采用積分約束的方式施加電流,如式(9)所示??赏ㄟ^有限元軟件中PDE模塊中的逐點約束實現(xiàn)此功能。

    式中,0為準(zhǔn)各向同性超導(dǎo)股線所通的運行電流。

    1.2.2 傳熱模塊的建立

    為了求解股線的溫度分布,需建立股線的傳熱模型。式(5)中的未知變量除了磁場強度還有溫度,通過傳熱模塊與磁場模塊的聯(lián)合求解即可得到股線的2個未知量:磁場強度與溫度。本文對超導(dǎo)帶材的熱力學(xué)參數(shù)采用均一化策略,以點熱源的形式添加熱脈沖,建立仿真模型。熱源的添加如圖4所示,熱源是=0截面處以中心點為圓點,熱脈沖周期為200 ms。

    傳熱模塊的控制方程為

    式中,為密度(kg/m3);和分別為比定壓熱容(J?kg-1?K-1)和熱導(dǎo)率(W?m-1?K-1);d和j分別為外加的熱脈沖(W/m3)以及單位體積的銅每秒產(chǎn)生的焦耳熱(W/m3)。

    圖4 準(zhǔn)各向同性超導(dǎo)股線的熱脈沖施加

    由圖4可以看出,超導(dǎo)線芯的帶材可以分為兩部分:一部分是與平面平行的平行帶材,另一部分是與平面垂直的垂直帶材。由于超導(dǎo)帶材的熱導(dǎo)率為各向異性[25],股線超導(dǎo)線芯部分的熱導(dǎo)率可寫成矩陣形式為

    式中,為單根超導(dǎo)帶材中不同材料的體積占比,見表2;k為超導(dǎo)帶材中不同材料的熱導(dǎo)率。

    由于本文模擬的是股線在零場液氮冷卻條件下的失超特性,傳熱模塊初始條件可寫為

    超導(dǎo)股線的長度為=140 mm,為減小計算量,仿真計算時僅需研究股線的一半(/2=70 mm),股線的中點位置為=0 mm,邊界條件為

    式中,n=[nxnynz]為垂直于銅包套外表面的單位法向量;為熱流密度(J×m-2×s-1);為液氮與銅的換熱系數(shù)[26](W?m-2?K-1);ext為液氮溫度(77 K)。

    2 仿真結(jié)果與分析

    2.1 準(zhǔn)各向同性超導(dǎo)股線溫度分布、電流密度分布與磁場分布

    本節(jié)利用上述失超模型研究了外界熱干擾下股線的失超特性。仿真過程中作了如下假設(shè):

    (1)股線各處的界面接觸熱阻(TCR)為0,即界面熱接觸為理想接觸。

    (2)超導(dǎo)帶材溫度高于臨界溫度(90 K)時其電阻率視為銅的電阻率,且銅的電阻率隨溫度的變化而變化[27]。

    (3)不考慮超導(dǎo)帶材終端接觸電阻對失超特性的影響。

    仿真過程如下:首先,對準(zhǔn)各向同性超導(dǎo)股線施加運行電流(100 ms后增加到3 000 A),經(jīng)過一段時間待股線穩(wěn)定后,在股線的中央以點熱源的方式施加一個熱脈沖。熱脈沖周期為200 ms,熱脈沖期間加熱功率恒定。

    未施加熱干擾時,通過仿真得出該股線臨界電流為4 188 A,之后給股線施加3 000 A的運行電流,待其穩(wěn)定后對股線進行分析,此時股線的磁場強度分布如圖5所示。

    自場下對此超導(dǎo)股線施加熱脈沖,當(dāng)干擾熱量為1.40 J/mm3時,所施加的熱干擾未達到最小失超能,股線發(fā)生失超恢復(fù)行為。準(zhǔn)各向同性超導(dǎo)股線被施加熱干擾后5 s內(nèi)的溫度分布(=0 mm,= 0 mm)如圖6所示,從圖中可以看出,隨著時間的增長,股線的失超區(qū)域逐步縮小直至整根股線恢復(fù)超導(dǎo)狀態(tài)。

    圖5 準(zhǔn)各向同性超導(dǎo)股線通流3 000 A的磁感應(yīng)強度分布

    圖6 超導(dǎo)股線通流3 000 A施加1.40 J/mm3干擾熱量后5 s內(nèi)的溫度分布

    股線通3 000 A電流時在=0 mm處徑向的溫度、電流密度及磁場分布如圖7所示。從圖中可以看出,0.2 s時=0 mm處的股線溫度高于臨界溫度,超導(dǎo)股線發(fā)生失超,電流轉(zhuǎn)移到銅基底層和銅包套,磁場分布也相應(yīng)發(fā)生變化。當(dāng)=5 s時=0 mm處的股線發(fā)生失超恢復(fù)行為,溫度降低至臨界溫度以下,電流重新全部轉(zhuǎn)回超導(dǎo)層內(nèi)。

    逐步增加干擾熱量至1.48 J/mm3時,股線發(fā)生失超傳播行為,12 s內(nèi)的溫度分布(=0 mm,= 0 mm)如圖8所示,失超區(qū)域隨著時間增大不斷擴散,故此時超導(dǎo)股線的最小失超能為1.48 J/mm3。

    由于此時發(fā)生了失超傳播行為,需對除熱源施加位置=0 mm處以外其他位置的溫度、電流密度進行觀測,這里選擇股線模型的中點=35 mm處,以反映股線的失超傳播過程。股線通流3 000 A時在=35 mm處徑向的溫度、電流密度以及磁場分布如圖9所示。

    從圖9中可以看出,0.2 s時股線=35 mm處還未發(fā)生失超,電流集中在超導(dǎo)層。當(dāng)=10 s時,失超區(qū)域已從熱點所在區(qū)域擴散到=35 mm處,股線發(fā)生失超傳播行為。在此,記錄該處發(fā)生失超的開始時刻為=10 s,后續(xù)將以此為基礎(chǔ)計算股線失超傳播速度。

    (a)=0.2 s (b)=5 s

    圖7 施加1.40 J/mm3干擾熱量后超導(dǎo)股線=0 mm處徑向溫度、電流密度及磁場分布

    Fig.7 The temperature field, current distribution and magnetic field of the quasi-isotropic strand under the heat energy of 1.40 J/mm3

    圖8 超導(dǎo)股線通流3 000 A施加1.48 J/mm3干擾熱量后12 s內(nèi)的溫度分布

    2.2 準(zhǔn)各向同性超導(dǎo)股線的最小失超能與失超傳播速度

    超導(dǎo)體在正常運行時,外界施加在超導(dǎo)體上的干擾熱量會誘發(fā)超導(dǎo)體失超,此時正常態(tài)下的超導(dǎo)帶電阻急劇增大,運行電流轉(zhuǎn)移到金屬基體中產(chǎn)生焦耳熱,如果這部分焦耳熱較小可以沿著超導(dǎo)股線進行熱擴散并通過液氮與銅包套的熱交換被直接導(dǎo)出,失超區(qū)域不會擴散。而當(dāng)股線中焦耳熱的積累大于熱擴散散失的熱量時,熱量逐漸積累,進而導(dǎo)致相鄰的超導(dǎo)體溫度隨之升高而發(fā)生失超現(xiàn)象,失超區(qū)域擴大,超導(dǎo)股線出現(xiàn)失超傳播行為。

    (a)=0.2 s (b)=10 s

    圖9 施加1.48 J/mm3干擾熱量后超導(dǎo)股線=35 mm處徑向溫度、電流密度及磁場分布圖

    Fig.9 The temperature field, current distribution and magnetic field of the quasi-isotropic strand under the heat energy of 1.48 J/mm3

    由此可見,外界施加在超導(dǎo)體上的干擾熱量的大小會影響股線是否發(fā)生失超行為。能夠引起超導(dǎo)股線失超的最小干擾熱量即被定義為最小失超能(Minimum Quench Energy, MQE)。

    本文所建立的失超模型中,外界所施加的干擾熱量d為

    式中,heat為熱脈沖功率體密度(W×m-3);h為加熱區(qū)域(半徑為2 mm的球體)的體積;d為熱脈沖的持續(xù)時間(0.2 s)。

    增大干擾熱量d的值直至超導(dǎo)股線失超,MQE可被表示為

    式中,dmin為誘發(fā)超導(dǎo)股線發(fā)生失超的最小干擾 熱量。

    失超傳播速度(QPV)可表示為

    式中,為股線軸上兩不同位置間的距離(cm);+1和分別為兩位置的溫度達到臨界溫度所需的時間。MQE與QPV隨運行電流的變化曲線如圖10所示。

    圖10 股線MQE與QPV隨運行電流的變化曲線

    由圖10可以看出,準(zhǔn)各向同性超導(dǎo)股線在運行電流為2 600~3 400 A時其最小失超能由2.17 J/mm3降至1.04 J/mm3,失超傳播速度由0.61 cm/s增加到1.45 cm/s。股線的最小失超能隨運行電流的增大而減小,失超傳播速度隨運行電流的增大而增大。當(dāng)運行電流為2 600 A時,即歸一化電流為0.59時,MQE可達到2.17 J/mm3,此時股線最高熱點溫度已達190 K,然而在實際應(yīng)用中REBCO超導(dǎo)帶材臨界熱應(yīng)力[28]所對應(yīng)的溫度約為200 K[29],因此股線在低電流運行時盡管不易發(fā)生失超行為,一旦由于外界干擾發(fā)生失超傳播,過高的熱點溫度非??赡軙p壞股線,同時低電流下股線的失超傳播速度極慢,不利于失超保護裝置的及時檢測,因此超導(dǎo)股線的運行電流大小極大地影響失超保護裝置參數(shù)的設(shè)置。為了減小股線的熱點溫度,可采取增加股線中銅的含量[30],利用其高熱導(dǎo)率,將熱量在股線徑向方向上及時將熱量傳導(dǎo)出去,然而銅含量的增加會減緩失超傳播速度,不利于失超檢測。另外,減小股線內(nèi)部的界面熱接觸電阻也會降低熱點溫度,同增加股線中銅含量一樣,會減緩失超傳播速度[31]。

    3 結(jié)論

    本文對點狀熱干源擾誘導(dǎo)的準(zhǔn)各向同性超導(dǎo)股線(Q-IS)的失超特性建立了三維電磁熱模型,仿真了股線的失超行為及恢復(fù)過程,得到了不同運行電流下的最小失超能與失超傳播速度,結(jié)論如下:

    股線的最小失超能隨運行電流的增大而減小,失超傳播速度隨運行電流的增大而增大。銅包套的存在使得股線的熱穩(wěn)定性有了顯著提高,銅包套通過及時分流避免超導(dǎo)體因熱點溫度過高而造成自身不可逆的損壞,矛盾的是,緩慢的失超傳播速度不利于對股線的失超檢測與失超保護,因此設(shè)計股線的銅包套尺寸時不僅要考慮降低超導(dǎo)體熱點的溫度,還要考慮其失超傳播速度不可過低,以避免超過失超保護裝置的失超檢測時間。

    本文所建立的針對準(zhǔn)各向同性超導(dǎo)股線考慮自場的三維磁熱失超模型可推廣應(yīng)用到其他超導(dǎo)導(dǎo)體失超特性的研究中去。

    [1] 黨衛(wèi)軍, 孫奇珍, 薛藝為, 等. 施工缺陷對超導(dǎo)電纜中間接頭內(nèi)電場分布的影響[J]. 江蘇電機工程, 2020, 39(5): 23-29.

    Dang Weijun, Sun Qizhen, Xue Yiwei, et al. Influ- ence of construction defects on electric field distri- bution of superconducting cable intermediate joints[J]. Jiangsu Electrical Engineering, 2020, 39(5): 23-29.

    [2] 李顯皓, 徐穎, 任麗, 等. 非均勻高溫超導(dǎo)帶材對CORC電纜失超特性的影響研究[J]. 電工技術(shù)學(xué)報, 2022, 37(19): 5044-5055.

    Li Xianhao, Xu Ying, Ren Li, et al. Influence of non-uniform high temperature superconducting tapes on quench characteristics of CORC cable[J]. Transa- ctions of China Electrotechnical Society, 2022, 37(19): 5044-5055.

    [3] 張國民, 陳建輝, 邱清泉, 等. 超導(dǎo)直流能源管道的研究進展[J]. 電工技術(shù)學(xué)報, 2021, 36(21): 4389- 4398, 4428.

    Zhang Guomin, Chen Jianhui, Qiu Qingquan, et al. Research progress on the superconducting DC energy pipeline[J]. Transactions of China Electrotechnical Society, 2021, 36(21): 4389-4398, 4428.

    [4] 祝乘風(fēng), 厲彥忠, 譚宏博, 等. 熱擾動沖擊下的高溫超導(dǎo)電纜失超恢復(fù)特性[J]. 電工技術(shù)學(xué)報, 2021, 36(18): 3884-3890.

    Zhu Chengfeng, Li Yanzhong, Tan Hongbo, et al. Numerical analysis on the quench and recovery of the high temperature superconducting cable subjected to thermal disturbance[J]. Transactions of China Elec- trotechnical Society, 2021, 36(18): 3884-3890.

    [5] 秦偉, 馬育華, 呂剛, 等. 一種可用于低真空管道的高溫超導(dǎo)無鐵心直線感應(yīng)磁懸浮電機[J]. 電工技術(shù)學(xué)報, 2022, 37(16): 4038-4046.

    Qin Wei, Ma Yuhua, Lü Gang, et al. Analyzing and designing a novel coreless linear induction maglev motor for low vacuum pipeline[J]. Transactions of China Electrotechnical Society, 2022, 37(16): 4038-4046.

    [6] 龔珺, 諸嘉慧, 方進, 等. 電阻型高溫超導(dǎo)限流器暫態(tài)電阻特性分析[J]. 電工技術(shù)學(xué)報, 2018, 33(9): 2130-2138.

    Gong Jun, Zhu Jiahui, Fang Jin, et al. Transient resistance analysis of resistive high temperature superconducting current limiter[J]. Transactions of China Electrotechnical Society, 2018, 33(9): 2130- 2138.

    [7] 夏芳敏, 王醒東. 實用化高溫超導(dǎo)帶材的應(yīng)用進展[J]. 新材料產(chǎn)業(yè), 2014(4): 48-52.

    Xia Fangmin, Wang Xingdong. Progress in the application of practical HTS tapes[J]. Advanced Materials Industry, 2014(4): 48-52.

    [8] Fleiter J, Ballarino A, Bottura L, et al. Charac- terization of roebel cables for potential use in high-field magnets[J]. IEEE Transactions on Applied Superconductivity, 2015, 25(3): 1-4.

    [9] Goldacker W, Grilli F, Pardo E, et al. Roebel cables from REBCO coated conductors: a one-century-old concept for the superconductivity of the future[J]. Superconductor Science and Technology, 2014, 27(9): 093001.

    [10] Komeda T, Amemiya N, Tsukamoto T, et al. Experimental comparison of AC loss in REBCO roebel cables consisting of six strands and ten strands[J]. IEEE Transactions on Applied Super- conductivity, 2014, 24(3): 1-5.

    [11] van der Laan D C. YBa2Cu3O7-coated conductor cabling for low ac-loss and high-field magnet applications[J]. Superconductor Science and Tech- nology, 2009, 22(6): 065013.

    [12] van der Laan D C, Weiss J D, Noyes P, et al. Record current density of 344 A mm-24.2 K and 17 T in CORC?accelerator magnet cables[J]. Superconductor Science and Technology, 2016, 29(5): 055009.

    [13] Takayasu M, Chiesa L, Bromberg L, et al. Cabling method for high current conductors made of HTS tapes[J]. IEEE Transactions on Applied Super- conductivity, 2011, 21(3): 2340-2344.

    [14] Wang Yinshun, Baasansuren S, Xue Chi, et al. Development of a quasi-isotropic strand stacked by 2G wires[J]. IEEE Transactions on Applied Super- conductivity, 2016, 26(4): 1-6.

    [15] Chan W K, Masson P J, Luongo C, et al. Three- dimensional micrometer-scale modeling of quenching in high-aspect-ratio YBa2Cu3O7-coated conductor tapes-part I: Model development and validation[J]. IEEE Transactions on Applied Superconductivity, 2010, 20(6): 2370-2380.

    [16] Chan W K, Schwartz J. Three-dimensional micrometer- scale modeling of quenching in high-aspect-ratio YBa2Cu3O7-coated conductor tapes-part II: Influ- ence of geometric and material properties and impli- cations for conductor engineering and magnet design[J]. IEEE Transactions on Applied Super- conductivity, 2011, 21(6): 3628-3634.

    [17] Roy F, Therasse M, Dutoit B, et al. Numerical studies of the quench propagation in coated conductors for fault current limiters[J]. IEEE Transactions on Applied Superconductivity, 2009, 19(3): 2496-2499.

    [18] Wang Yawei, Zheng Jinxing, Zhu Zixuan, et al. Quench behavior of high-temperature superconductor (RE) Ba2Cu3OCORC cable[J]. Journal of Physics D: Applied Physics, 2019, 52(34): 345303.

    [19] Kang Rui, Uglietti D, Wesche R, et al. Quench simulation of REBCO cable-in-conduit conductor with twisted stacked-tape cable[J]. IEEE Transactions on Applied Superconductivity, 2020, 30(1): 1-7.

    [20] 林良真. 超導(dǎo)電性及其應(yīng)用[M]. 北京: 北京工業(yè)大學(xué)出版社, 1998.

    [21] Hong Z, Campbell A M, Coombs T A. Numerical solution of critical state in superconductivity by finite element software[J]. Superconductor Science and Technology, 2006, 19(12): 1246-1252.

    [22] Brambilla R, Grilli F, Martini L. Development of an edge-element model for AC loss computation of high-temperature superconductors[J]. Superconductor Science and Technology, 2007, 20(1): 16-24.

    [23] Ma Jun, Geng Jianzhao, Chan Wan kan, et al. A temperature-dependent multilayer model for direct current carrying HTS coated-conductors under per- pendicular AC magnetic fields[J]. Superconductor Science and Technology, 2020, 33(4): 045007.

    [24] Grilli F, Sirois F, Zerme?o V M R, et al. Self- consistent modeling of thecof HTS devices: how accurate do models really need to be?[J]. IEEE Transa- ctions on Applied Superconductivity, 2014, 24(6): 1-8.

    [25] Bagrets N, Goldacker W, Jung A, et al. Thermal properties of ReBCO copper stabilized supercon- ducting tapes[J]. IEEE Transactions on Applied Superconductivity, 2013, 23(3): 6600303.

    [26] Li Tingting, Wang Yinshun, Shi Chenjie, et al. Thermal stability of a quasi-isotropic strand made from coated conductors[J]. IEEE Transactions on Applied Superconductivity, 2016, 26(4): 1-5.

    [27] Casali M, Breschi M, Ribani P L. Two-dimensional anisotropic model of YBCO coated conductors[J]. IEEE Transactions on Applied Superconductivity, 2015, 25(1): 1-12.

    [28] Duan Yujie, Gao Yuanwen. Delamination and current-carrying degradation behavior of epoxy- impregnated superconducting coil winding with 2G HTS tape caused by thermal stress[J]. AIP Advances, 2020, 10(2): 1-13.

    [29] Marinucci C, Bottura L, Calvi M, et al. Quench analysis of a high-current forced-flow HTS conductor model for fusion magnets[J]. IEEE Transactions on Applied Superconductivity, 2011, 21(3): 2445-2448.

    [30] Kang Rui, Uglietti D, Song Yuntao. To the opti- mization of quench performance for a large REBCO cable-in-conduit conductor[J]. IEEE Transactions on Applied Superconductivity, 2021, 31(2): 1-9.

    [31] Kang Rui, Uglietti D, Wesche R, et al. Quench simulation of REBCO cable-in-conduit conductor with twisted stacked-tape cable[J]. IEEE Transactions on Applied Superconductivity, 2020, 30(1): 1-7.

    Quench Characteristics of Quasi-Isotropic Superconducting Strand Triggered by Point Thermal Disturbance

    12111

    (1. The State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources North China Electric Power University Beijing 102206 China 2. State Grid Shandong Electric Power Company Jinan Power Supply Company Jinan 250000 China)

    With the development of second-generation (2G) high-temperature superconductors (HTS) and the progress in cryogenic technology, 2G HTS has been widely applied in superconducting power technology due to its excellent mechanical properties, complete diamagnetism, and zero resistance characteristic at liquid nitrogen temperatures. However, the practical applications of HTS face a key issue: the quench behavior of HTS, mainly embodied in the minimum quench energy (MQE) and quench propagation velocity (QPV). This paper studies the quench characteristics of quasi-isotropic superconducting (Q-IS) strands when subjected to a point thermal disturbance and the influence of different operating currents for quench detection and protection of HTS devices in power systems.

    Firstly, the structure and parameters of the superconducting strand used in the simulation model are described in detail, and the geometric model of a 140 mm long Q-IS strand immersed in liquid nitrogen is established. Considering the effect of the self-magnetic field, a 3D electric-magnetic-thermal simulation model is established to analyze the thermal stability of the strand using the finite element method (FEM). At the beginning of the simulation, the operating current is applied to the HTS strand until the strand operates stably. Then, a thermal pulse is added to the strand’s center to simulate a point thermal disturbance. The thermal disturbance power is constant during the 200 ms duration of the thermal disturbance application.

    The thermal stability simulation results show that when the operating current is 3 000 A and the disturbance energy is 1.40 J/mm3, the strand exhibits an apparent quench recovery behavior, suggesting that 1.40 J/mm3does not reach the MQE of the strand. As the disturbance energy increases to 1.48 J/mm3, the quench region diffuses continuously over time, resulting in quench propagation with the strand. Under this operating current, the MQE of the superconducting strand is 1.48 J/mm3, and the QPV is 0.79 cm/s. Then, MQEs and QPVs under various operating currents are simulated by the same method. The MQE of the Q-IS strand decreases from 2.17 J/mm3to 1.04 J/mm3, and the QPV increases from 0.61 cm/s to 1.45 cm/s as the operating current increases from 2 600 A to 3 400 A.

    It is also shown that the MQE of the Q-IS strand decreases and the QPV increases with the increase of operating current. The current sharing of the copper sheath can significantly improve the thermal stability of the Q-IS strand. However, it is essential to consider the size design of the copper sheath to ensure it is manageable, potentially causing a delay in quench detection by the protection device. The electric-magnetic-thermal simulation model established in this paper can also be extended to study the thermal stability of Q-IS strands in external fields and other superconducting strands.

    Quasi-isotropic superconducting strand, quench behavior, minimum quench energy (MQE), quench propagation velocity (QPV)

    TM26

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

    國家自然科學(xué)基金資助項目(52277025, 51877083)。

    2022-10-18

    2022-12-19

    王起悅 女,1999年生,碩士研究生,研究方向為超導(dǎo)電力技術(shù)。E-mail: superyue990826@163.com

    皮 偉 男,1979年生,副教授,研究方向為超導(dǎo)電力技術(shù)。E-mail: ppiiwei@ncepu.edu.cn(通信作者)

    (編輯 郭麗軍)

    猜你喜歡
    股線傳播速度帶材
    2022年7月板帶材產(chǎn)量匯總表
    2022年6月板帶材產(chǎn)量匯總表
    2022年5月板帶材產(chǎn)量匯總表
    2022年3月板帶材產(chǎn)量匯總表
    代謝綜合征患者臂踝脈搏波傳播速度與頸動脈粥樣硬化的關(guān)系
    大型汽輪發(fā)電機線棒股線斷裂原因分析和預(yù)防
    大型水輪發(fā)電機定子換位線棒匝間短路股線損耗
    新雷
    一類廣義canmassa—Holm方程的無限傳播速度與漸近行為
    多截面場路耦合法和解析法對比分析水輪機定子線棒4 種換位下的環(huán)流
    機電工程(2015年12期)2015-11-18 12:26:30
    嫩草影院精品99| 在线观看一区二区三区| 国产精品国产三级专区第一集| 中文欧美无线码| 黄片无遮挡物在线观看| 男女下面进入的视频免费午夜| 免费av毛片视频| 亚洲精品色激情综合| 色综合站精品国产| 少妇裸体淫交视频免费看高清| 91在线精品国自产拍蜜月| 日日摸夜夜添夜夜爱| 一级毛片我不卡| 欧美xxxx黑人xx丫x性爽| 亚洲综合精品二区| 国产美女午夜福利| 免费看av在线观看网站| 别揉我奶头 嗯啊视频| 极品少妇高潮喷水抽搐| 亚洲在线观看片| 亚洲欧美成人精品一区二区| 一个人看视频在线观看www免费| 亚洲av二区三区四区| 肉色欧美久久久久久久蜜桃 | 亚洲av中文av极速乱| av国产免费在线观看| 99热网站在线观看| 一级毛片电影观看| 一个人免费在线观看电影| 欧美成人一区二区免费高清观看| 成年人午夜在线观看视频 | 久久鲁丝午夜福利片| 九九在线视频观看精品| 亚洲激情五月婷婷啪啪| 五月伊人婷婷丁香| 亚洲熟妇中文字幕五十中出| 狂野欧美激情性xxxx在线观看| 毛片一级片免费看久久久久| 日本欧美国产在线视频| 老司机影院成人| 三级国产精品片| 丝瓜视频免费看黄片| 国产淫语在线视频| 日韩国内少妇激情av| 全区人妻精品视频| 18禁动态无遮挡网站| 内地一区二区视频在线| 十八禁国产超污无遮挡网站| 波多野结衣巨乳人妻| 久久久久国产网址| 亚洲经典国产精华液单| 国产老妇女一区| 亚洲精品亚洲一区二区| 搡女人真爽免费视频火全软件| 精品久久久久久久久av| 校园人妻丝袜中文字幕| 亚洲最大成人手机在线| 国产乱来视频区| 尤物成人国产欧美一区二区三区| 久久这里有精品视频免费| 看免费成人av毛片| 肉色欧美久久久久久久蜜桃 | 男人舔女人下体高潮全视频| 国产成年人精品一区二区| 日韩强制内射视频| 久久韩国三级中文字幕| 99久国产av精品国产电影| 国产精品一二三区在线看| 色尼玛亚洲综合影院| 少妇的逼水好多| 成人亚洲精品一区在线观看 | 成年人午夜在线观看视频 | 亚洲成人av在线免费| 久久99精品国语久久久| 日韩大片免费观看网站| www.色视频.com| 欧美另类一区| 毛片一级片免费看久久久久| 最近的中文字幕免费完整| 91精品国产九色| 免费少妇av软件| 99re6热这里在线精品视频| 欧美日韩一区二区视频在线观看视频在线 | 大香蕉97超碰在线| 少妇熟女欧美另类| 直男gayav资源| 中文资源天堂在线| 欧美区成人在线视频| 亚洲av中文av极速乱| 高清日韩中文字幕在线| 啦啦啦中文免费视频观看日本| 亚洲av免费高清在线观看| 久久精品久久久久久久性| 久久久a久久爽久久v久久| 亚洲无线观看免费| 日本免费在线观看一区| 亚洲精品日韩av片在线观看| 久久久国产一区二区| 久久精品久久精品一区二区三区| 禁无遮挡网站| 亚洲人成网站在线播| 中文字幕av在线有码专区| av天堂中文字幕网| av在线观看视频网站免费| 亚洲一级一片aⅴ在线观看| 国产黄片视频在线免费观看| 午夜精品在线福利| 国产v大片淫在线免费观看| 青春草视频在线免费观看| 黄色一级大片看看| 亚洲精品aⅴ在线观看| av国产免费在线观看| 色吧在线观看| 免费av毛片视频| 床上黄色一级片| 99热网站在线观看| 女人十人毛片免费观看3o分钟| 免费观看的影片在线观看| 久久草成人影院| 天堂av国产一区二区熟女人妻| 白带黄色成豆腐渣| 男女边摸边吃奶| 亚洲国产精品专区欧美| 欧美另类一区| 春色校园在线视频观看| 女的被弄到高潮叫床怎么办| 看黄色毛片网站| 亚洲成色77777| 一区二区三区高清视频在线| 欧美日韩精品成人综合77777| 久久久久久久久久黄片| 国产男人的电影天堂91| 2018国产大陆天天弄谢| 国产视频首页在线观看| 性插视频无遮挡在线免费观看| 日韩 亚洲 欧美在线| 91精品国产九色| 亚洲三级黄色毛片| 人妻一区二区av| 日本免费在线观看一区| freevideosex欧美| 国产 亚洲一区二区三区 | 丰满人妻一区二区三区视频av| 99九九线精品视频在线观看视频| 亚洲国产最新在线播放| 三级毛片av免费| 直男gayav资源| av在线亚洲专区| 国产精品久久久久久久电影| 国产黄色免费在线视频| 国产精品av视频在线免费观看| 九色成人免费人妻av| 麻豆乱淫一区二区| 最近2019中文字幕mv第一页| 男人狂女人下面高潮的视频| 99久久人妻综合| 国产激情偷乱视频一区二区| 国产午夜精品论理片| 国产有黄有色有爽视频| 天天躁夜夜躁狠狠久久av| 亚洲精品成人久久久久久| 国产亚洲91精品色在线| 国产成人a∨麻豆精品| 秋霞在线观看毛片| 亚洲av不卡在线观看| 精品久久久久久久久久久久久| 日韩电影二区| 国产毛片a区久久久久| 夜夜看夜夜爽夜夜摸| 亚洲乱码一区二区免费版| 99久久中文字幕三级久久日本| 国产精品爽爽va在线观看网站| 国产精品久久久久久av不卡| 五月天丁香电影| 成年av动漫网址| 小蜜桃在线观看免费完整版高清| 亚洲色图av天堂| 午夜福利高清视频| 草草在线视频免费看| 色哟哟·www| 午夜免费激情av| 特大巨黑吊av在线直播| 亚洲人成网站在线观看播放| 最近中文字幕2019免费版| 国产不卡一卡二| 嫩草影院精品99| 麻豆成人午夜福利视频| 国产精品久久久久久精品电影小说 | 久久久久久久国产电影| 久久亚洲国产成人精品v| 国产一级毛片七仙女欲春2| 亚洲自拍偷在线| 国产精品久久久久久精品电影| a级毛片免费高清观看在线播放| 国精品久久久久久国模美| 亚洲成人久久爱视频| 伊人久久国产一区二区| 国产探花在线观看一区二区| 夫妻性生交免费视频一级片| 街头女战士在线观看网站| 国产av不卡久久| 一个人观看的视频www高清免费观看| 非洲黑人性xxxx精品又粗又长| 久久久久久久久久久丰满| 久久久欧美国产精品| 欧美激情久久久久久爽电影| 亚洲精品视频女| 亚洲综合色惰| 色播亚洲综合网| 久久精品夜色国产| 国产日韩欧美在线精品| 欧美最新免费一区二区三区| 亚洲欧美精品专区久久| 国产男人的电影天堂91| 毛片一级片免费看久久久久| 99热全是精品| 国产免费又黄又爽又色| 国产老妇女一区| 男女下面进入的视频免费午夜| av在线天堂中文字幕| 午夜精品在线福利| 久久久久久久久久久丰满| 国产探花极品一区二区| 免费高清在线观看视频在线观看| 精品一区二区免费观看| av卡一久久| 亚洲av在线观看美女高潮| www.av在线官网国产| 免费大片18禁| 韩国av在线不卡| 精品久久久久久久人妻蜜臀av| 三级国产精品欧美在线观看| 五月玫瑰六月丁香| 在线 av 中文字幕| 午夜免费观看性视频| 22中文网久久字幕| 激情五月婷婷亚洲| 一级黄片播放器| 亚洲精品一二三| 男女边摸边吃奶| or卡值多少钱| 啦啦啦啦在线视频资源| 久久精品国产亚洲av涩爱| 国产一区二区三区av在线| 又粗又硬又长又爽又黄的视频| 99久国产av精品国产电影| 啦啦啦韩国在线观看视频| 日本黄大片高清| 久久久国产一区二区| 免费黄色在线免费观看| 国产亚洲精品av在线| av播播在线观看一区| 日本免费a在线| 免费看日本二区| 国产黄片美女视频| 日韩大片免费观看网站| 久久久久网色| 青春草视频在线免费观看| 久久精品人妻少妇| 听说在线观看完整版免费高清| 水蜜桃什么品种好| 亚洲精品一二三| 哪个播放器可以免费观看大片| 嫩草影院新地址| 亚洲av电影不卡..在线观看| 日韩,欧美,国产一区二区三区| 少妇丰满av| 少妇猛男粗大的猛烈进出视频 | 天堂影院成人在线观看| 精品熟女少妇av免费看| 久久精品久久久久久噜噜老黄| 美女国产视频在线观看| 久久久久九九精品影院| 国产成人a区在线观看| 国产老妇女一区| 国产成人精品婷婷| 午夜激情欧美在线| 亚洲欧美日韩东京热| 国产成人91sexporn| 久久久久九九精品影院| 免费大片18禁| 国产在线一区二区三区精| 91狼人影院| av免费观看日本| 欧美成人a在线观看| 直男gayav资源| 99视频精品全部免费 在线| 久久久久久久久久黄片| 男人和女人高潮做爰伦理| 大又大粗又爽又黄少妇毛片口| 身体一侧抽搐| 激情五月婷婷亚洲| 嘟嘟电影网在线观看| 国产精品1区2区在线观看.| 国产一区亚洲一区在线观看| 久久人人爽人人片av| 久久鲁丝午夜福利片| 丰满人妻一区二区三区视频av| 亚洲欧美中文字幕日韩二区| 国产成人午夜福利电影在线观看| av免费观看日本| 一区二区三区乱码不卡18| 亚洲三级黄色毛片| 中文欧美无线码| 精品久久久久久久末码| 亚洲人成网站高清观看| 国产成人精品久久久久久| 国产亚洲午夜精品一区二区久久 | 日韩欧美精品免费久久| 一级毛片aaaaaa免费看小| 精品人妻偷拍中文字幕| 天美传媒精品一区二区| 亚洲国产精品国产精品| 三级国产精品欧美在线观看| 乱人视频在线观看| 晚上一个人看的免费电影| 国产高清不卡午夜福利| 天堂√8在线中文| 国产视频内射| 成年版毛片免费区| 精品一区二区三区视频在线| 日韩成人伦理影院| 久久久久精品久久久久真实原创| 午夜福利在线观看吧| 夜夜爽夜夜爽视频| 91aial.com中文字幕在线观看| 性色avwww在线观看| 中文字幕av在线有码专区| 少妇人妻一区二区三区视频| 性色avwww在线观看| 亚洲精华国产精华液的使用体验| 少妇丰满av| 男人爽女人下面视频在线观看| 亚洲欧美成人精品一区二区| 热99在线观看视频| 成人性生交大片免费视频hd| 久久亚洲国产成人精品v| 久久久久精品久久久久真实原创| 国产精品熟女久久久久浪| 国产午夜精品论理片| 久久久久久久久大av| 精品人妻熟女av久视频| 高清视频免费观看一区二区 | 亚洲国产精品国产精品| 久久久亚洲精品成人影院| 欧美激情在线99| 国产色婷婷99| 亚洲精品日韩av片在线观看| 午夜福利视频1000在线观看| 国产又色又爽无遮挡免| 日韩精品有码人妻一区| 国产老妇伦熟女老妇高清| 亚洲性久久影院| 亚洲av.av天堂| 久久久久精品性色| 国内精品一区二区在线观看| 免费不卡的大黄色大毛片视频在线观看 | 日韩一区二区视频免费看| 美女主播在线视频| 国产在视频线精品| 波多野结衣巨乳人妻| av国产免费在线观看| 春色校园在线视频观看| 国产成人一区二区在线| 亚洲av男天堂| 久久久久性生活片| 国产老妇女一区| 99九九线精品视频在线观看视频| 天堂影院成人在线观看| 久久99精品国语久久久| 男插女下体视频免费在线播放| 18+在线观看网站| 热99在线观看视频| 97精品久久久久久久久久精品| 成人鲁丝片一二三区免费| 网址你懂的国产日韩在线| 97在线视频观看| 欧美高清性xxxxhd video| 久久精品人妻少妇| 能在线免费看毛片的网站| 高清视频免费观看一区二区 | 国产淫语在线视频| 又大又黄又爽视频免费| 男人爽女人下面视频在线观看| 国产精品一区二区三区四区免费观看| 夫妻午夜视频| 精品人妻一区二区三区麻豆| 日本与韩国留学比较| 美女cb高潮喷水在线观看| 欧美人与善性xxx| 一级片'在线观看视频| 中文字幕制服av| 久久久亚洲精品成人影院| videossex国产| 欧美人与善性xxx| 九色成人免费人妻av| 国产大屁股一区二区在线视频| 黄片wwwwww| 午夜激情欧美在线| 亚洲国产高清在线一区二区三| 亚洲av一区综合| 免费观看精品视频网站| 少妇丰满av| 99热全是精品| 欧美xxⅹ黑人| 美女内射精品一级片tv| 日韩在线高清观看一区二区三区| 最近最新中文字幕免费大全7| 日本黄色片子视频| 大话2 男鬼变身卡| 性插视频无遮挡在线免费观看| 中文资源天堂在线| 天堂√8在线中文| 日日啪夜夜爽| 日韩伦理黄色片| 人妻系列 视频| 亚洲av日韩在线播放| 免费观看性生交大片5| 免费人成在线观看视频色| 免费播放大片免费观看视频在线观看| 亚洲人与动物交配视频| 国产有黄有色有爽视频| 欧美zozozo另类| 欧美高清成人免费视频www| 91av网一区二区| 国产精品一区www在线观看| 夫妻性生交免费视频一级片| 夫妻午夜视频| 最近最新中文字幕大全电影3| 亚洲熟女精品中文字幕| 亚洲国产精品成人久久小说| 国产老妇女一区| 国精品久久久久久国模美| 小蜜桃在线观看免费完整版高清| 七月丁香在线播放| 欧美日韩亚洲高清精品| 久久热精品热| 亚洲aⅴ乱码一区二区在线播放| 美女被艹到高潮喷水动态| 免费av毛片视频| 18禁动态无遮挡网站| 只有这里有精品99| 日韩欧美精品免费久久| 精品酒店卫生间| 欧美一区二区亚洲| 亚洲av中文字字幕乱码综合| 亚州av有码| 亚洲国产精品国产精品| 看非洲黑人一级黄片| 久久久久性生活片| 国产永久视频网站| 毛片女人毛片| 最近中文字幕2019免费版| 亚洲aⅴ乱码一区二区在线播放| 搞女人的毛片| 中国国产av一级| 国产在线男女| 五月玫瑰六月丁香| 丝袜喷水一区| 搡女人真爽免费视频火全软件| 蜜桃亚洲精品一区二区三区| 亚洲av中文av极速乱| 国产精品久久久久久精品电影小说 | 国产乱人偷精品视频| 日产精品乱码卡一卡2卡三| 在线免费十八禁| 全区人妻精品视频| 一区二区三区免费毛片| 九草在线视频观看| 啦啦啦中文免费视频观看日本| 亚洲精华国产精华液的使用体验| 少妇的逼好多水| 亚洲在线观看片| 美女黄网站色视频| 精品人妻一区二区三区麻豆| 亚洲av成人精品一二三区| 特大巨黑吊av在线直播| 亚洲精品成人久久久久久| 麻豆av噜噜一区二区三区| 日韩一本色道免费dvd| 麻豆成人午夜福利视频| 国产免费福利视频在线观看| 久久久久九九精品影院| 欧美区成人在线视频| 亚洲aⅴ乱码一区二区在线播放| 内射极品少妇av片p| 欧美另类一区| 精品久久久久久久人妻蜜臀av| 亚洲熟妇中文字幕五十中出| 亚洲av成人精品一区久久| 亚洲一级一片aⅴ在线观看| 午夜久久久久精精品| 一级毛片 在线播放| 99久久人妻综合| 国产伦精品一区二区三区视频9| 久久久久久久大尺度免费视频| 国产精品无大码| 又爽又黄无遮挡网站| 人人妻人人澡欧美一区二区| 韩国高清视频一区二区三区| 久久久午夜欧美精品| 亚洲欧美一区二区三区国产| 精品人妻偷拍中文字幕| 精品人妻熟女av久视频| 日韩一本色道免费dvd| 亚洲人与动物交配视频| 91午夜精品亚洲一区二区三区| 久99久视频精品免费| 欧美精品一区二区大全| 在线免费十八禁| 天天躁夜夜躁狠狠久久av| 国产黄色视频一区二区在线观看| 天天躁夜夜躁狠狠久久av| 欧美xxxx性猛交bbbb| 日韩一本色道免费dvd| 久久精品夜色国产| 久久久午夜欧美精品| 啦啦啦啦在线视频资源| 2021天堂中文幕一二区在线观| 韩国高清视频一区二区三区| 亚洲成人精品中文字幕电影| 国产视频内射| 亚洲成人精品中文字幕电影| 免费看日本二区| 久久精品久久久久久噜噜老黄| av线在线观看网站| 美女大奶头视频| 亚洲精品乱码久久久v下载方式| 亚洲国产精品专区欧美| 乱码一卡2卡4卡精品| 亚洲18禁久久av| 精品国内亚洲2022精品成人| 国产亚洲午夜精品一区二区久久 | 国产一区二区亚洲精品在线观看| 日韩一区二区视频免费看| 亚洲av电影在线观看一区二区三区 | 日日啪夜夜爽| 亚洲一级一片aⅴ在线观看| 成人美女网站在线观看视频| 亚洲一级一片aⅴ在线观看| 国产日韩欧美在线精品| 精品午夜福利在线看| 天堂网av新在线| 国产成年人精品一区二区| 国产日韩欧美在线精品| 精品久久久久久久久av| 赤兔流量卡办理| 激情 狠狠 欧美| 欧美激情国产日韩精品一区| 久久韩国三级中文字幕| 国产老妇女一区| 在线 av 中文字幕| 中文天堂在线官网| 国产精品国产三级国产专区5o| 亚洲欧美日韩无卡精品| 中文字幕制服av| 精品一区二区免费观看| 国产一区有黄有色的免费视频 | 一级黄片播放器| 亚洲欧美日韩无卡精品| 国产91av在线免费观看| 久久久久久伊人网av| 国产免费又黄又爽又色| 欧美另类一区| 特级一级黄色大片| 最近中文字幕高清免费大全6| h日本视频在线播放| 一本久久精品| 日韩欧美 国产精品| 三级国产精品片| 亚洲精品日韩av片在线观看| 成人亚洲精品av一区二区| www.色视频.com| 欧美一级a爱片免费观看看| 国国产精品蜜臀av免费| 中国美白少妇内射xxxbb| 亚洲欧美日韩卡通动漫| 国产男女超爽视频在线观看| 免费观看的影片在线观看| 国产爱豆传媒在线观看| 又爽又黄无遮挡网站| 国产成人91sexporn| 精品一区二区三卡| 国产精品精品国产色婷婷| 亚洲在久久综合| 好男人视频免费观看在线| 日本黄大片高清| 人妻系列 视频| 国产不卡一卡二| 国产精品麻豆人妻色哟哟久久 | 欧美高清成人免费视频www| 国产黄色小视频在线观看| 成人高潮视频无遮挡免费网站| videossex国产| 边亲边吃奶的免费视频| 建设人人有责人人尽责人人享有的 | av在线播放精品| 国产有黄有色有爽视频| 亚洲欧美一区二区三区国产| av国产免费在线观看| 十八禁国产超污无遮挡网站| 久久精品综合一区二区三区| 亚洲性久久影院| 美女内射精品一级片tv| 午夜福利网站1000一区二区三区| 三级国产精品片| 亚洲怡红院男人天堂| 国产探花在线观看一区二区| 国产综合精华液| 插逼视频在线观看| 熟女人妻精品中文字幕| av免费观看日本| 国产精品麻豆人妻色哟哟久久 | 久久久久免费精品人妻一区二区| 久久韩国三级中文字幕|