• <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
    国产美女午夜福利| 1024手机看黄色片| 欧美一区二区亚洲| 观看美女的网站| 我要搜黄色片| 日韩av在线大香蕉| 久久亚洲真实| 亚洲黑人精品在线| 夜夜爽天天搞| 日韩欧美精品免费久久 | 听说在线观看完整版免费高清| 99国产综合亚洲精品| 搡老妇女老女人老熟妇| 成人午夜高清在线视频| 男女视频在线观看网站免费| 男插女下体视频免费在线播放| 宅男免费午夜| 免费电影在线观看免费观看| 亚洲精品色激情综合| 国产一区二区在线av高清观看| 听说在线观看完整版免费高清| 啪啪无遮挡十八禁网站| 午夜免费激情av| 久久九九热精品免费| 麻豆国产av国片精品| bbb黄色大片| 国产成人欧美在线观看| 蜜桃久久精品国产亚洲av| 最新美女视频免费是黄的| 欧美最黄视频在线播放免费| netflix在线观看网站| 免费一级毛片在线播放高清视频| 国产激情欧美一区二区| 国产亚洲精品av在线| 岛国在线观看网站| 欧美不卡视频在线免费观看| 天天一区二区日本电影三级| av福利片在线观看| 身体一侧抽搐| 又紧又爽又黄一区二区| 少妇高潮的动态图| 亚洲成人免费电影在线观看| 亚洲电影在线观看av| 中文字幕av成人在线电影| 亚洲国产精品999在线| 免费av毛片视频| 19禁男女啪啪无遮挡网站| 免费观看人在逋| 性色av乱码一区二区三区2| 一边摸一边抽搐一进一小说| 极品教师在线免费播放| 亚洲熟妇中文字幕五十中出| 日韩欧美三级三区| 亚洲精品粉嫩美女一区| 国产欧美日韩精品亚洲av| 乱人视频在线观看| 91在线精品国自产拍蜜月 | 免费av不卡在线播放| 亚洲成人免费电影在线观看| 久久欧美精品欧美久久欧美| 欧美av亚洲av综合av国产av| 色吧在线观看| 九九热线精品视视频播放| 18禁美女被吸乳视频| 日韩欧美精品免费久久 | 免费看日本二区| 最近最新中文字幕大全电影3| 18+在线观看网站| 岛国在线观看网站| 久久久久久九九精品二区国产| 男人舔女人下体高潮全视频| 国产三级黄色录像| 91久久精品国产一区二区成人 | 在线观看舔阴道视频| av天堂中文字幕网| 男人舔奶头视频| 久久这里只有精品中国| 国产一区二区激情短视频| 久久6这里有精品| av女优亚洲男人天堂| 色综合欧美亚洲国产小说| 高潮久久久久久久久久久不卡| 国产精品99久久久久久久久| 亚洲精品在线美女| 日韩成人在线观看一区二区三区| 精品乱码久久久久久99久播| 亚洲狠狠婷婷综合久久图片| 18禁裸乳无遮挡免费网站照片| 日本黄色片子视频| 99热这里只有是精品50| 内射极品少妇av片p| 人妻夜夜爽99麻豆av| 亚洲精品成人久久久久久| 国产欧美日韩精品亚洲av| 色播亚洲综合网| bbb黄色大片| 午夜久久久久精精品| 亚洲最大成人中文| 亚洲成人久久性| 啦啦啦免费观看视频1| 99在线人妻在线中文字幕| 色播亚洲综合网| 欧美黄色淫秽网站| 亚洲av不卡在线观看| 久久伊人香网站| 午夜免费男女啪啪视频观看 | 国产高清videossex| 免费人成在线观看视频色| 18+在线观看网站| 国产极品精品免费视频能看的| 日韩欧美国产一区二区入口| 精品国产美女av久久久久小说| 长腿黑丝高跟| 欧美日韩乱码在线| 欧美xxxx黑人xx丫x性爽| 熟女人妻精品中文字幕| 校园春色视频在线观看| 美女被艹到高潮喷水动态| 少妇高潮的动态图| 中文字幕av成人在线电影| 久久久成人免费电影| 亚洲中文字幕日韩| 波多野结衣高清无吗| 色综合亚洲欧美另类图片| 极品教师在线免费播放| 久久久国产成人免费| 欧美3d第一页| 亚洲欧美日韩高清在线视频| 国产毛片a区久久久久| 亚洲精品亚洲一区二区| 国产精品久久久久久精品电影| 国产精品精品国产色婷婷| 99riav亚洲国产免费| 色视频www国产| 人妻丰满熟妇av一区二区三区| 岛国在线观看网站| 亚洲专区中文字幕在线| 亚洲,欧美精品.| 日韩欧美免费精品| 日韩亚洲欧美综合| 日本在线视频免费播放| 搡老熟女国产l中国老女人| 全区人妻精品视频| 村上凉子中文字幕在线| 老汉色av国产亚洲站长工具| 日韩av在线大香蕉| 欧美黄色淫秽网站| 丁香六月欧美| 国产亚洲精品久久久久久毛片| 丁香欧美五月| 琪琪午夜伦伦电影理论片6080| 啪啪无遮挡十八禁网站| 日韩有码中文字幕| 天堂影院成人在线观看| 免费看光身美女| 大型黄色视频在线免费观看| 91麻豆精品激情在线观看国产| 久久久久久大精品| 俺也久久电影网| 欧美成人免费av一区二区三区| 99久久精品一区二区三区| 欧美区成人在线视频| 国产野战对白在线观看| 3wmmmm亚洲av在线观看| 90打野战视频偷拍视频| 亚洲av免费高清在线观看| 欧洲精品卡2卡3卡4卡5卡区| 午夜日韩欧美国产| 日韩欧美 国产精品| 久久这里只有精品中国| 一夜夜www| 午夜激情福利司机影院| 啦啦啦韩国在线观看视频| 国产老妇女一区| 国产一区二区在线观看日韩 | 日本撒尿小便嘘嘘汇集6| 国产伦精品一区二区三区四那| 午夜两性在线视频| 久久亚洲精品不卡| av视频在线观看入口| 国产成人欧美在线观看| 9191精品国产免费久久| 国产高清视频在线播放一区| 欧美黑人欧美精品刺激| 亚洲aⅴ乱码一区二区在线播放| 97超级碰碰碰精品色视频在线观看| 99久久精品一区二区三区| 日本黄大片高清| 国产精品,欧美在线| a级毛片a级免费在线| а√天堂www在线а√下载| 亚洲成人久久爱视频| 中文字幕久久专区| 国产极品精品免费视频能看的| 午夜免费观看网址| 两个人视频免费观看高清| 99在线视频只有这里精品首页| 国产av一区在线观看免费| 国产午夜精品论理片| 色哟哟哟哟哟哟| 一a级毛片在线观看| 长腿黑丝高跟| 色精品久久人妻99蜜桃| 欧美一级毛片孕妇| 欧美乱码精品一区二区三区| 精品电影一区二区在线| 免费高清视频大片| 亚洲精品久久国产高清桃花| 色av中文字幕| 欧美乱色亚洲激情| 日本五十路高清| 最近视频中文字幕2019在线8| aaaaa片日本免费| 欧美日本视频| 久久久久久九九精品二区国产| 国产精品三级大全| 欧美日韩一级在线毛片| 99热精品在线国产| 男人舔奶头视频| 久久久精品欧美日韩精品| 三级毛片av免费| 波多野结衣巨乳人妻| 中文字幕av成人在线电影| 亚洲av日韩精品久久久久久密| 成人性生交大片免费视频hd| 看片在线看免费视频| 国内少妇人妻偷人精品xxx网站| 19禁男女啪啪无遮挡网站| 欧美日韩福利视频一区二区| 国产探花在线观看一区二区| 精品一区二区三区人妻视频| 国产高清视频在线观看网站| 久久九九热精品免费| 国内精品久久久久精免费| 精品一区二区三区人妻视频| 19禁男女啪啪无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 免费在线观看日本一区| 精品一区二区三区人妻视频| 亚洲精品粉嫩美女一区| 非洲黑人性xxxx精品又粗又长| 一区二区三区国产精品乱码| 97碰自拍视频| 国产亚洲欧美在线一区二区| 无人区码免费观看不卡| 韩国av一区二区三区四区| 黄色日韩在线| 高清日韩中文字幕在线| a在线观看视频网站| 九色国产91popny在线| 亚洲电影在线观看av| 99在线人妻在线中文字幕| 国产麻豆成人av免费视频| 看黄色毛片网站| 久久人人精品亚洲av| 岛国在线免费视频观看| 国产精品永久免费网站| 亚洲精品在线美女| 久久香蕉国产精品| 午夜激情福利司机影院| 国产免费男女视频| 国产一区在线观看成人免费| 色播亚洲综合网| 国产又黄又爽又无遮挡在线| 老熟妇仑乱视频hdxx| 叶爱在线成人免费视频播放| 老司机午夜福利在线观看视频| 成年版毛片免费区| 国产高清视频在线观看网站| 51国产日韩欧美| 成人永久免费在线观看视频| 亚洲自拍偷在线| 99国产精品一区二区蜜桃av| 丰满人妻熟妇乱又伦精品不卡| 久久香蕉精品热| 少妇的逼好多水| 亚洲国产色片| 午夜日韩欧美国产| 免费在线观看影片大全网站| 午夜福利欧美成人| 亚洲七黄色美女视频| 亚洲不卡免费看| 国产99白浆流出| 国产一区二区激情短视频| 美女高潮的动态| 婷婷精品国产亚洲av| 丰满人妻一区二区三区视频av | 中文字幕av在线有码专区| 国产真人三级小视频在线观看| 久久亚洲真实| 久久久精品大字幕| 国产精品影院久久| 最近最新中文字幕大全免费视频| svipshipincom国产片| 久久精品91蜜桃| 国产一区二区在线av高清观看| 国产69精品久久久久777片| 久久久久国产精品人妻aⅴ院| 宅男免费午夜| 少妇人妻精品综合一区二区 | a在线观看视频网站| 成人一区二区视频在线观看| 99国产极品粉嫩在线观看| 亚洲欧美日韩卡通动漫| 午夜久久久久精精品| 精品电影一区二区在线| 一夜夜www| 九色国产91popny在线| 亚洲精品亚洲一区二区| 69av精品久久久久久| 国产aⅴ精品一区二区三区波| av在线蜜桃| 天天添夜夜摸| 欧美另类亚洲清纯唯美| 亚洲中文字幕一区二区三区有码在线看| 国产黄a三级三级三级人| 欧美中文日本在线观看视频| 99视频精品全部免费 在线| 十八禁网站免费在线| 夜夜夜夜夜久久久久| 国产亚洲精品久久久久久毛片| 久久精品综合一区二区三区| 久久精品国产99精品国产亚洲性色| av天堂中文字幕网| 亚洲精品国产精品久久久不卡| 女人被狂操c到高潮| 国产视频一区二区在线看| 国产欧美日韩精品亚洲av| 国产视频内射| 99国产精品一区二区三区| 久久亚洲精品不卡| 嫩草影院精品99| 国产v大片淫在线免费观看| 黄色成人免费大全| 亚洲国产欧洲综合997久久,| 午夜激情欧美在线| 色老头精品视频在线观看| 国产高清激情床上av| 成年人黄色毛片网站| 久久久久久久久中文| 国产精品久久久久久亚洲av鲁大| 国产一区二区激情短视频| 十八禁人妻一区二区| 久久久久精品国产欧美久久久| 99久久久亚洲精品蜜臀av| 在线观看美女被高潮喷水网站 | 中文字幕精品亚洲无线码一区| 成人国产综合亚洲| 丰满人妻一区二区三区视频av | 麻豆成人av在线观看| 欧美性猛交╳xxx乱大交人| 久久国产精品人妻蜜桃| 国产真实伦视频高清在线观看 | 午夜两性在线视频| 日韩欧美免费精品| 看免费av毛片| 国产高清激情床上av| 国产麻豆成人av免费视频| 又紧又爽又黄一区二区| 欧美日韩中文字幕国产精品一区二区三区| 精品无人区乱码1区二区| 久久国产精品影院| 听说在线观看完整版免费高清| av黄色大香蕉| 神马国产精品三级电影在线观看| 免费看十八禁软件| 91九色精品人成在线观看| 美女黄网站色视频| 性欧美人与动物交配| 亚洲 欧美 日韩 在线 免费| 欧美日韩乱码在线| 国产色婷婷99| avwww免费| 可以在线观看的亚洲视频| 婷婷亚洲欧美| 午夜福利免费观看在线| 日本在线视频免费播放| av天堂在线播放| 男插女下体视频免费在线播放| 久久精品亚洲精品国产色婷小说| 丰满人妻一区二区三区视频av | 一个人看的www免费观看视频| 亚洲av日韩精品久久久久久密| 日韩精品青青久久久久久| 最近最新免费中文字幕在线| 亚洲av中文字字幕乱码综合| 久久亚洲精品不卡| 很黄的视频免费| 欧美三级亚洲精品| 精品免费久久久久久久清纯| 国产激情欧美一区二区| 国产伦精品一区二区三区视频9 | 老司机深夜福利视频在线观看| 久久精品国产清高在天天线| x7x7x7水蜜桃| 国产亚洲精品久久久久久毛片| 一个人免费在线观看电影| 亚洲不卡免费看| 在线a可以看的网站| 超碰av人人做人人爽久久 | 亚洲av成人不卡在线观看播放网| 精品久久久久久久末码| 久久精品国产综合久久久| 久久中文看片网| 精品午夜福利视频在线观看一区| 亚洲专区中文字幕在线| 最后的刺客免费高清国语| 内射极品少妇av片p| 麻豆国产97在线/欧美| 欧洲精品卡2卡3卡4卡5卡区| 女人高潮潮喷娇喘18禁视频| 男人和女人高潮做爰伦理| 白带黄色成豆腐渣| 国产精品精品国产色婷婷| 国内久久婷婷六月综合欲色啪| 午夜福利在线在线| 免费观看的影片在线观看| 成人国产综合亚洲| 亚洲狠狠婷婷综合久久图片| av黄色大香蕉| 2021天堂中文幕一二区在线观| 久久性视频一级片| 国产精品国产高清国产av| 最近视频中文字幕2019在线8| 欧美一级a爱片免费观看看| 婷婷六月久久综合丁香| 成人鲁丝片一二三区免费| 精品人妻偷拍中文字幕| 嫁个100分男人电影在线观看| 国产成+人综合+亚洲专区| 久久精品国产清高在天天线| 亚洲av不卡在线观看| 欧美+亚洲+日韩+国产| 亚洲精华国产精华精| 熟妇人妻久久中文字幕3abv| 天堂√8在线中文| 精品久久久久久久久久久久久| 最好的美女福利视频网| 亚洲人与动物交配视频| 日本成人三级电影网站| 观看美女的网站| 丁香六月欧美| 18禁黄网站禁片免费观看直播| 午夜精品久久久久久毛片777| 精品不卡国产一区二区三区| 亚洲七黄色美女视频| 欧美高清成人免费视频www| 少妇丰满av| av片东京热男人的天堂| 国产亚洲精品久久久久久毛片| 国产黄色小视频在线观看| 啦啦啦观看免费观看视频高清| 热99在线观看视频| 女人被狂操c到高潮| 精品人妻偷拍中文字幕| 国产精华一区二区三区| 国产精品98久久久久久宅男小说| 亚洲男人的天堂狠狠| 三级国产精品欧美在线观看| 97超级碰碰碰精品色视频在线观看| 国产一区二区在线观看日韩 | 欧美日本亚洲视频在线播放| 精品99又大又爽又粗少妇毛片 | 热99在线观看视频| 国产色婷婷99| 香蕉丝袜av| 亚洲精品一卡2卡三卡4卡5卡| 成年人黄色毛片网站| 神马国产精品三级电影在线观看| 精品99又大又爽又粗少妇毛片 | 男女做爰动态图高潮gif福利片| 淫妇啪啪啪对白视频| 久久久久久久午夜电影| 又黄又爽又免费观看的视频| 亚洲av免费高清在线观看| 波多野结衣高清无吗| 色av中文字幕| 免费高清视频大片| 搞女人的毛片| 国产亚洲av嫩草精品影院| 欧美一区二区精品小视频在线| 首页视频小说图片口味搜索| 国产伦精品一区二区三区视频9 | 亚洲av免费在线观看| 免费一级毛片在线播放高清视频| 岛国在线免费视频观看| 身体一侧抽搐| 久久精品国产99精品国产亚洲性色| 他把我摸到了高潮在线观看| 中文字幕av在线有码专区| 欧美激情在线99| 国产一区二区三区视频了| 夜夜夜夜夜久久久久| 99久久精品热视频| 国产av一区在线观看免费| 丰满的人妻完整版| 亚洲第一电影网av| 亚洲美女视频黄频| 久久久久精品国产欧美久久久| 久久精品夜夜夜夜夜久久蜜豆| 99热这里只有精品一区| 草草在线视频免费看| 熟女少妇亚洲综合色aaa.| 婷婷六月久久综合丁香| www.熟女人妻精品国产| 一进一出好大好爽视频| 国产精品 欧美亚洲| 搡老岳熟女国产| 两人在一起打扑克的视频| 在线天堂最新版资源| 搡老妇女老女人老熟妇| 成年人黄色毛片网站| 国产视频一区二区在线看| 免费看十八禁软件| 久久精品亚洲精品国产色婷小说| 在线观看美女被高潮喷水网站 | 久久久久久九九精品二区国产| 国产精品久久视频播放| av专区在线播放| 午夜福利欧美成人| 成熟少妇高潮喷水视频| 国产乱人视频| 别揉我奶头~嗯~啊~动态视频| 色老头精品视频在线观看| 国产成人福利小说| 美女被艹到高潮喷水动态| 亚洲激情在线av| 级片在线观看| 国产午夜精品久久久久久一区二区三区 | 三级男女做爰猛烈吃奶摸视频| 国产精品女同一区二区软件 | 色综合婷婷激情| 久久久久久国产a免费观看| 欧美性猛交黑人性爽| a级一级毛片免费在线观看| 国产精品av视频在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 在线观看av片永久免费下载| 国产一区二区激情短视频| 99热精品在线国产| 国产一区在线观看成人免费| 久久午夜亚洲精品久久| 日本黄色视频三级网站网址| 69人妻影院| 国产精品一区二区三区四区免费观看 | 亚洲avbb在线观看| 午夜福利欧美成人| 一级黄片播放器| 久久人人精品亚洲av| 欧美zozozo另类| 亚洲av五月六月丁香网| 亚洲aⅴ乱码一区二区在线播放| 啦啦啦观看免费观看视频高清| 国产午夜精品久久久久久一区二区三区 | 白带黄色成豆腐渣| 欧美黑人巨大hd| 亚洲av二区三区四区| 90打野战视频偷拍视频| 无限看片的www在线观看| 91在线精品国自产拍蜜月 | 精品国产美女av久久久久小说| 国产中年淑女户外野战色| 身体一侧抽搐| 亚洲成人久久爱视频| 夜夜看夜夜爽夜夜摸| 好男人在线观看高清免费视频| 国产亚洲精品av在线| 18禁在线播放成人免费| 国产高清videossex| 一区福利在线观看| 性欧美人与动物交配| 搡老岳熟女国产| 一级黄片播放器| 免费无遮挡裸体视频| 久久精品91无色码中文字幕| 免费看光身美女| 最近最新免费中文字幕在线| 欧美乱码精品一区二区三区| 亚洲国产欧洲综合997久久,| 久久伊人香网站| 亚洲真实伦在线观看| 99久久精品国产亚洲精品| 亚洲电影在线观看av| 69av精品久久久久久| 白带黄色成豆腐渣| 99在线人妻在线中文字幕| 天堂影院成人在线观看| 88av欧美| 最近最新中文字幕大全电影3| 五月玫瑰六月丁香| 国产色爽女视频免费观看| www.999成人在线观看| 日韩欧美在线二视频| 岛国视频午夜一区免费看| 99久久精品一区二区三区| www.999成人在线观看| 中文字幕av在线有码专区| 免费高清视频大片| 亚洲午夜理论影院| 亚洲熟妇中文字幕五十中出| 天堂av国产一区二区熟女人妻| 国内毛片毛片毛片毛片毛片| 午夜日韩欧美国产| 狠狠狠狠99中文字幕| 欧美日本视频| 亚洲国产精品成人综合色| 内射极品少妇av片p| 国产精品 国内视频| 熟妇人妻久久中文字幕3abv| 丰满乱子伦码专区| 一本综合久久免费| 99国产综合亚洲精品| 757午夜福利合集在线观看| 叶爱在线成人免费视频播放| 欧美不卡视频在线免费观看|