• <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
    国产毛片a区久久久久| 丰满乱子伦码专区| 大香蕉97超碰在线| 尤物成人国产欧美一区二区三区| 欧美+日韩+精品| av.在线天堂| 26uuu在线亚洲综合色| 综合色av麻豆| 久久久久性生活片| 中文字幕久久专区| 久久精品91蜜桃| 国产一区二区三区av在线| 国产亚洲91精品色在线| 久久久色成人| 联通29元200g的流量卡| 久久精品91蜜桃| 久久99热这里只有精品18| 最近中文字幕高清免费大全6| 中文字幕精品亚洲无线码一区| 国产日韩欧美在线精品| 亚洲中文字幕一区二区三区有码在线看| 日本wwww免费看| 久久综合国产亚洲精品| 成人鲁丝片一二三区免费| 看片在线看免费视频| 身体一侧抽搐| 国产黄a三级三级三级人| 一区二区三区免费毛片| 又爽又黄a免费视频| 看免费成人av毛片| 亚洲成人中文字幕在线播放| 免费观看人在逋| 久久精品久久久久久久性| 国产精品野战在线观看| 久久久亚洲精品成人影院| 免费不卡的大黄色大毛片视频在线观看 | 97在线视频观看| 美女cb高潮喷水在线观看| 亚洲图色成人| 插逼视频在线观看| 日韩欧美三级三区| 精品欧美国产一区二区三| 直男gayav资源| 精品酒店卫生间| 日韩欧美精品v在线| 我要看日韩黄色一级片| 干丝袜人妻中文字幕| 国产亚洲精品av在线| 久久久久久九九精品二区国产| 亚州av有码| 一区二区三区免费毛片| 国产一区有黄有色的免费视频 | 纵有疾风起免费观看全集完整版 | 18禁在线播放成人免费| 女人被狂操c到高潮| 欧美高清成人免费视频www| 欧美一区二区精品小视频在线| 国产高清视频在线观看网站| 日韩亚洲欧美综合| 国产精品野战在线观看| 2022亚洲国产成人精品| 综合色丁香网| 久久久久久久国产电影| 网址你懂的国产日韩在线| 日韩中字成人| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲av天美| 看免费成人av毛片| 看非洲黑人一级黄片| 女人十人毛片免费观看3o分钟| 国产亚洲一区二区精品| 黄色配什么色好看| 国产伦精品一区二区三区视频9| 欧美zozozo另类| 如何舔出高潮| 人体艺术视频欧美日本| 久久人人爽人人爽人人片va| 狂野欧美白嫩少妇大欣赏| 精品国内亚洲2022精品成人| 日韩精品青青久久久久久| 男女边吃奶边做爰视频| 久久久久性生活片| 国产伦理片在线播放av一区| 婷婷色综合大香蕉| 午夜老司机福利剧场| av黄色大香蕉| 国产极品精品免费视频能看的| 日韩三级伦理在线观看| 成人性生交大片免费视频hd| 波多野结衣巨乳人妻| 国产乱人偷精品视频| 欧美日韩综合久久久久久| 国产精品av视频在线免费观看| 女人十人毛片免费观看3o分钟| 午夜福利成人在线免费观看| 欧美不卡视频在线免费观看| 久久99热这里只有精品18| 性色avwww在线观看| 久久人妻av系列| 亚洲电影在线观看av| 99久久中文字幕三级久久日本| 男人的好看免费观看在线视频| 日韩欧美国产在线观看| 两个人的视频大全免费| 淫秽高清视频在线观看| 特大巨黑吊av在线直播| 国产精品一区二区在线观看99 | 免费播放大片免费观看视频在线观看 | 国产单亲对白刺激| 老司机影院毛片| 国产成人免费观看mmmm| 简卡轻食公司| 欧美激情国产日韩精品一区| 高清日韩中文字幕在线| 99在线视频只有这里精品首页| 免费观看精品视频网站| 黄片无遮挡物在线观看| 99久久无色码亚洲精品果冻| 观看免费一级毛片| 身体一侧抽搐| 亚洲av一区综合| 久久草成人影院| 看十八女毛片水多多多| 长腿黑丝高跟| 一级黄片播放器| 色5月婷婷丁香| 国产真实伦视频高清在线观看| 日本免费a在线| 成人无遮挡网站| 久久久色成人| 嫩草影院入口| 99久久无色码亚洲精品果冻| 亚洲欧美日韩卡通动漫| 国产午夜精品久久久久久一区二区三区| 免费av观看视频| 非洲黑人性xxxx精品又粗又长| 午夜亚洲福利在线播放| 91久久精品国产一区二区成人| 六月丁香七月| 久久久色成人| 欧美成人精品欧美一级黄| 成人综合一区亚洲| 久久久久国产网址| 1000部很黄的大片| 亚洲欧美成人综合另类久久久 | 欧美97在线视频| 国产免费男女视频| 国产精品综合久久久久久久免费| 91aial.com中文字幕在线观看| 国产精品人妻久久久影院| 亚洲av电影在线观看一区二区三区 | 中国美白少妇内射xxxbb| 看免费成人av毛片| 精品国内亚洲2022精品成人| 人妻制服诱惑在线中文字幕| 亚洲天堂国产精品一区在线| 美女高潮的动态| 亚洲aⅴ乱码一区二区在线播放| a级毛片免费高清观看在线播放| 色哟哟·www| 欧美一区二区国产精品久久精品| 一区二区三区免费毛片| 热99re8久久精品国产| 精品久久久久久久久久久久久| av在线亚洲专区| 99视频精品全部免费 在线| 日韩欧美精品v在线| 97超视频在线观看视频| av视频在线观看入口| 在现免费观看毛片| 国产精品伦人一区二区| 久久精品熟女亚洲av麻豆精品 | 亚洲中文字幕一区二区三区有码在线看| 久久久精品94久久精品| 久久久久久久久久久丰满| 我的老师免费观看完整版| 美女xxoo啪啪120秒动态图| 国产精品蜜桃在线观看| 欧美三级亚洲精品| 免费观看人在逋| 国产在视频线在精品| 久久热精品热| 亚洲精品影视一区二区三区av| 亚洲人成网站高清观看| 黄色日韩在线| 亚洲av中文av极速乱| 亚洲18禁久久av| 搡女人真爽免费视频火全软件| 国产女主播在线喷水免费视频网站 | 在线天堂最新版资源| 91狼人影院| 国产极品精品免费视频能看的| 日韩av在线免费看完整版不卡| 成年女人永久免费观看视频| 亚洲va在线va天堂va国产| 国产一级毛片在线| 一级爰片在线观看| 黄色配什么色好看| 免费看a级黄色片| 亚洲国产欧美人成| 免费观看的影片在线观看| 晚上一个人看的免费电影| kizo精华| 亚洲中文字幕一区二区三区有码在线看| 国产伦精品一区二区三区视频9| 嫩草影院新地址| 亚洲精品456在线播放app| 日本免费在线观看一区| 最近视频中文字幕2019在线8| 色尼玛亚洲综合影院| 久久热精品热| 亚洲高清免费不卡视频| 国产av一区在线观看免费| 亚洲欧美成人精品一区二区| 大话2 男鬼变身卡| 永久免费av网站大全| 色网站视频免费| 国内精品一区二区在线观看| 国产淫语在线视频| 国产精品无大码| 91av网一区二区| 精品久久久久久电影网 | 啦啦啦韩国在线观看视频| 日本-黄色视频高清免费观看| 国产成人91sexporn| 一区二区三区四区激情视频| 免费看av在线观看网站| 欧美日本视频| 日韩人妻高清精品专区| 国产黄a三级三级三级人| 精品一区二区三区视频在线| 99热这里只有是精品50| 亚洲国产高清在线一区二区三| 亚洲av免费在线观看| 亚洲va在线va天堂va国产| 亚洲欧美成人精品一区二区| 中文精品一卡2卡3卡4更新| 又爽又黄无遮挡网站| 精品一区二区三区人妻视频| 欧美一区二区亚洲| 男女那种视频在线观看| 看非洲黑人一级黄片| 国模一区二区三区四区视频| 男人狂女人下面高潮的视频| 亚洲在线观看片| 亚洲天堂国产精品一区在线| 免费看av在线观看网站| 亚洲婷婷狠狠爱综合网| 午夜福利在线观看吧| 日韩欧美精品免费久久| 亚洲av不卡在线观看| 国产探花在线观看一区二区| 精品一区二区三区人妻视频| 一级二级三级毛片免费看| 又黄又爽又刺激的免费视频.| 亚洲精品自拍成人| 久久国内精品自在自线图片| 亚洲成色77777| 我要看日韩黄色一级片| 可以在线观看毛片的网站| 色噜噜av男人的天堂激情| 能在线免费观看的黄片| 夜夜看夜夜爽夜夜摸| 中国国产av一级| 午夜久久久久精精品| 亚洲成人精品中文字幕电影| 最近的中文字幕免费完整| 在线a可以看的网站| 国产精品精品国产色婷婷| 久久人妻av系列| 哪个播放器可以免费观看大片| 大香蕉97超碰在线| 在线观看66精品国产| 一级二级三级毛片免费看| 国产高清三级在线| 日韩制服骚丝袜av| 女人被狂操c到高潮| 亚洲成人av在线免费| 色吧在线观看| 最近中文字幕高清免费大全6| 亚洲精品456在线播放app| 国产综合懂色| 国产日韩欧美在线精品| 91精品国产九色| 啦啦啦啦在线视频资源| 亚洲人成网站在线播| 搞女人的毛片| 九九在线视频观看精品| 午夜a级毛片| 成人无遮挡网站| 免费观看精品视频网站| 22中文网久久字幕| 欧美一区二区国产精品久久精品| 国产精品国产三级国产专区5o | 国产高清国产精品国产三级 | 中文字幕久久专区| 在线观看av片永久免费下载| 国产亚洲最大av| 国产单亲对白刺激| 国产精品一区二区在线观看99 | 亚洲高清免费不卡视频| 久久久久国产网址| 观看美女的网站| 免费一级毛片在线播放高清视频| 日本午夜av视频| 男女视频在线观看网站免费| 美女高潮的动态| 亚洲精华国产精华液的使用体验| 大香蕉久久网| 国产熟女欧美一区二区| 成人亚洲欧美一区二区av| 纵有疾风起免费观看全集完整版 | 欧美高清成人免费视频www| 久久欧美精品欧美久久欧美| 国产高清国产精品国产三级 | 女人久久www免费人成看片 | 久久精品夜色国产| 搞女人的毛片| 久久精品影院6| 婷婷色av中文字幕| 免费无遮挡裸体视频| 少妇人妻一区二区三区视频| 视频中文字幕在线观看| 国产单亲对白刺激| 大香蕉久久网| 色吧在线观看| 禁无遮挡网站| 成人综合一区亚洲| 亚洲中文字幕日韩| 国产一区二区在线av高清观看| 国产av一区在线观看免费| 国产熟女欧美一区二区| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品日韩av片在线观看| 男人和女人高潮做爰伦理| 国产午夜精品论理片| 少妇丰满av| 国产精品一区二区三区四区久久| 哪个播放器可以免费观看大片| 欧美zozozo另类| 少妇人妻精品综合一区二区| 69人妻影院| 日本三级黄在线观看| 国产一级毛片七仙女欲春2| 超碰av人人做人人爽久久| 99久久成人亚洲精品观看| 久久久午夜欧美精品| 性插视频无遮挡在线免费观看| 亚洲三级黄色毛片| 少妇高潮的动态图| 久久久午夜欧美精品| 欧美一区二区国产精品久久精品| 亚洲高清免费不卡视频| 热99re8久久精品国产| 国产伦在线观看视频一区| 国产爱豆传媒在线观看| 观看美女的网站| 国产精品麻豆人妻色哟哟久久 | 精品少妇黑人巨大在线播放 | 午夜福利在线观看吧| 晚上一个人看的免费电影| 可以在线观看毛片的网站| 国产黄片美女视频| 亚洲精品乱码久久久久久按摩| 日韩大片免费观看网站 | 国产在视频线在精品| 亚洲av日韩在线播放| 偷拍熟女少妇极品色| 久久久久久久亚洲中文字幕| 国产真实伦视频高清在线观看| 亚洲av免费高清在线观看| 精品少妇黑人巨大在线播放 | 成人综合一区亚洲| 3wmmmm亚洲av在线观看| 能在线免费观看的黄片| 国产精品一区二区在线观看99 | 国产成人91sexporn| 久久精品久久久久久久性| 永久网站在线| 欧美zozozo另类| 国产成人午夜福利电影在线观看| 高清日韩中文字幕在线| 干丝袜人妻中文字幕| 国产免费一级a男人的天堂| 97热精品久久久久久| 色网站视频免费| 美女cb高潮喷水在线观看| 国产一级毛片七仙女欲春2| 日韩亚洲欧美综合| 亚洲精品自拍成人| 99热这里只有是精品在线观看| 神马国产精品三级电影在线观看| 桃色一区二区三区在线观看| 在线观看av片永久免费下载| 国内精品宾馆在线| 精品国产露脸久久av麻豆 | 亚洲一区高清亚洲精品| 日韩在线高清观看一区二区三区| 久久久久久九九精品二区国产| 亚洲成人中文字幕在线播放| 美女国产视频在线观看| 淫秽高清视频在线观看| 一本久久精品| 99热全是精品| 国产亚洲精品av在线| 亚洲精品亚洲一区二区| a级毛片免费高清观看在线播放| 久久精品国产99精品国产亚洲性色| 欧美三级亚洲精品| 国国产精品蜜臀av免费| 国产黄色视频一区二区在线观看 | 日韩中字成人| 中文乱码字字幕精品一区二区三区 | 免费av不卡在线播放| 国产精品人妻久久久影院| 国产精品国产三级国产专区5o | 国产精品女同一区二区软件| 欧美另类亚洲清纯唯美| 亚洲精品aⅴ在线观看| 水蜜桃什么品种好| 在线a可以看的网站| 色综合站精品国产| 精品99又大又爽又粗少妇毛片| www.av在线官网国产| 边亲边吃奶的免费视频| 美女黄网站色视频| 日本av手机在线免费观看| 非洲黑人性xxxx精品又粗又长| 国产亚洲av片在线观看秒播厂 | 色综合站精品国产| 能在线免费看毛片的网站| 免费人成在线观看视频色| 亚洲欧美精品专区久久| 一二三四中文在线观看免费高清| 午夜精品一区二区三区免费看| 97人妻精品一区二区三区麻豆| 最后的刺客免费高清国语| 亚洲欧美成人精品一区二区| 国产高清视频在线观看网站| 国产亚洲5aaaaa淫片| 国产午夜福利久久久久久| 亚洲av成人av| 熟女人妻精品中文字幕| 午夜精品国产一区二区电影 | 国产成年人精品一区二区| 一级毛片电影观看 | 亚洲精品国产成人久久av| 99久久精品国产国产毛片| 小蜜桃在线观看免费完整版高清| 亚洲av成人精品一二三区| 综合色丁香网| 久久精品国产亚洲网站| 又黄又爽又刺激的免费视频.| 少妇高潮的动态图| 免费大片18禁| 丝袜喷水一区| 精品人妻视频免费看| 国产一级毛片在线| 色吧在线观看| 久久久a久久爽久久v久久| 亚洲精品成人久久久久久| 日韩欧美国产在线观看| 99热网站在线观看| 国产淫语在线视频| 国国产精品蜜臀av免费| 免费一级毛片在线播放高清视频| 日本免费a在线| 特级一级黄色大片| 日韩高清综合在线| 成年女人看的毛片在线观看| 狂野欧美白嫩少妇大欣赏| 天天躁夜夜躁狠狠久久av| 久久久久网色| 91aial.com中文字幕在线观看| 久久久久性生活片| 九九在线视频观看精品| 嫩草影院入口| 国产精品.久久久| 国产高清有码在线观看视频| 天堂中文最新版在线下载 | 国产精品国产三级专区第一集| 能在线免费观看的黄片| 亚洲欧美精品自产自拍| 久久午夜福利片| 18+在线观看网站| 国产精品国产高清国产av| 九草在线视频观看| 国产私拍福利视频在线观看| 久99久视频精品免费| 亚洲成色77777| 好男人视频免费观看在线| 欧美+日韩+精品| 搡老妇女老女人老熟妇| 99视频精品全部免费 在线| 18禁裸乳无遮挡免费网站照片| 国内少妇人妻偷人精品xxx网站| 国语对白做爰xxxⅹ性视频网站| 熟女电影av网| 人妻系列 视频| 神马国产精品三级电影在线观看| 别揉我奶头 嗯啊视频| 色综合亚洲欧美另类图片| 亚洲在线自拍视频| 日本猛色少妇xxxxx猛交久久| 久久精品熟女亚洲av麻豆精品 | 性色avwww在线观看| 精品久久久噜噜| 久久久成人免费电影| 国产视频内射| 免费av毛片视频| 在线播放国产精品三级| 综合色av麻豆| kizo精华| 丰满少妇做爰视频| 亚洲国产高清在线一区二区三| 亚洲国产最新在线播放| 蜜臀久久99精品久久宅男| 色视频www国产| 99久久中文字幕三级久久日本| 国产成人精品一,二区| 深爱激情五月婷婷| 免费看光身美女| 性色avwww在线观看| 国产免费福利视频在线观看| 99久国产av精品国产电影| 亚洲乱码一区二区免费版| 久久精品国产99精品国产亚洲性色| 在线免费观看的www视频| 国产淫语在线视频| 一区二区三区四区激情视频| 看十八女毛片水多多多| 最近的中文字幕免费完整| 精品国产露脸久久av麻豆 | 午夜福利在线观看吧| 99视频精品全部免费 在线| 国产精品国产三级国产专区5o | 国产91av在线免费观看| 久久这里只有精品中国| 欧美另类亚洲清纯唯美| 97在线视频观看| 亚洲人成网站在线观看播放| 欧美成人一区二区免费高清观看| 日本黄色视频三级网站网址| 韩国高清视频一区二区三区| 精品久久久久久久久亚洲| 久久久精品大字幕| 色综合站精品国产| .国产精品久久| 久久国内精品自在自线图片| 精品熟女少妇av免费看| 成人av在线播放网站| 男的添女的下面高潮视频| 亚洲国产成人一精品久久久| 亚洲激情五月婷婷啪啪| 国产在视频线在精品| 久久精品影院6| 国产精品.久久久| 99久久精品国产国产毛片| 哪个播放器可以免费观看大片| 亚洲不卡免费看| 激情 狠狠 欧美| www.av在线官网国产| 最近2019中文字幕mv第一页| 国产精品福利在线免费观看| 国产日韩欧美在线精品| 久久久精品94久久精品| 亚洲成人av在线免费| 亚洲欧美中文字幕日韩二区| 日本wwww免费看| videos熟女内射| 视频中文字幕在线观看| 99久久九九国产精品国产免费| 亚洲va在线va天堂va国产| 一个人看的www免费观看视频| 国产美女午夜福利| 国产 一区精品| 精品人妻一区二区三区麻豆| 久久婷婷人人爽人人干人人爱| 国产免费又黄又爽又色| 国产精品日韩av在线免费观看| 在线a可以看的网站| 国产淫片久久久久久久久| 亚洲美女视频黄频| 黄色欧美视频在线观看| av卡一久久| 欧美三级亚洲精品| 久久久久久久国产电影| 亚州av有码| 黄色一级大片看看| 美女cb高潮喷水在线观看| 日韩中字成人| 国产色爽女视频免费观看| av视频在线观看入口| 高清午夜精品一区二区三区| 晚上一个人看的免费电影| 美女国产视频在线观看| 不卡视频在线观看欧美| 美女内射精品一级片tv| 亚洲真实伦在线观看| 中文字幕亚洲精品专区| 中文字幕熟女人妻在线| 国语自产精品视频在线第100页| 建设人人有责人人尽责人人享有的 | 91aial.com中文字幕在线观看| 国产精品久久电影中文字幕| 成人三级黄色视频| 波多野结衣巨乳人妻| 噜噜噜噜噜久久久久久91| 大话2 男鬼变身卡| a级毛色黄片| 国产精品一及| 波多野结衣巨乳人妻| 毛片女人毛片| 国产免费福利视频在线观看| 国产亚洲一区二区精品|