999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

快速陡脈沖重復(fù)電場下高頻變壓器絕緣介質(zhì)損耗與沖擊能量積聚特性

2023-03-11 09:40:04王威望李睿喆何杰峰張曉彤李盛濤
電工技術(shù)學(xué)報(bào) 2023年5期
關(guān)鍵詞:變壓器

王威望 李睿喆 何杰峰 張曉彤 李盛濤

快速陡脈沖重復(fù)電場下高頻變壓器絕緣介質(zhì)損耗與沖擊能量積聚特性

王威望 李睿喆 何杰峰 張曉彤 李盛濤

(電力設(shè)備電氣絕緣國家重點(diǎn)實(shí)驗(yàn)室(西安交通大學(xué)) 西安 710049)

電力電子變壓器中高頻變壓器須承受高頻非正弦與高溫的復(fù)雜應(yīng)力。隨著SiC等大功率半導(dǎo)體器件的使用,高頻變壓器中非正弦方波電壓上升與下降沿出現(xiàn)陡脈沖應(yīng)力,其d/d大于10 kV/ms,給高頻變壓器絕緣可靠性帶來巨大挑戰(zhàn)。高d/d包含大量高頻諧波,從而導(dǎo)致絕緣介質(zhì)損耗明顯增加;高頻、高d/d下絕緣遭受累積性沖擊電應(yīng)力,導(dǎo)致局部電場畸變,引發(fā)絕緣局部放電與損傷。該文圍繞d/d對(duì)高頻變壓器絕緣介質(zhì)損耗與沖擊能量積聚的影響展開研究。采用階躍響應(yīng)函數(shù)模擬高d/d脈沖電壓,結(jié)合絕緣介質(zhì)損耗分解與計(jì)算,提出了d/d與方波電壓疊加下絕緣介質(zhì)損耗的計(jì)算方法。結(jié)果表明,絕緣介質(zhì)損耗隨d/d的增大而增加,絕緣熱效應(yīng)明顯增加。根據(jù)絕緣在方波交變電場下能量儲(chǔ)存和釋放特性,分析了高d/d下絕緣能量積聚密度d。采用有限元仿真分析10 kW、10 kHz、1 000 V/750 V高頻變壓器絕緣電場、位移電流與沖擊能量密度。結(jié)果發(fā)現(xiàn),d隨著d/d增加而增大。結(jié)合介質(zhì)損耗計(jì)算,研究指出在方波電壓高d/d處(電壓極性反轉(zhuǎn)時(shí)),絕緣遭受累積電-熱沖擊應(yīng)力,造成絕緣損傷。研究結(jié)果為大容量高頻變壓器絕緣失效與設(shè)計(jì)提供指導(dǎo)。

高頻變壓器 d/d介質(zhì)損耗 絕緣損傷 能量沖擊功率密度 位移電流

0 引言

電力電子變壓器(Power Electronic Transformer,PET)是實(shí)現(xiàn)可再生與分布式能源接入、并網(wǎng)和靈活交直流輸配電的關(guān)鍵設(shè)備[1-4]。高頻變壓器是PET的核心部件,發(fā)揮著電壓變換、電氣隔離與能量傳輸?shù)戎匾饔肹5-6]。高頻變壓器的體積小、質(zhì)量輕、功率密度大、電-磁-熱-力耦合強(qiáng)[7-8]。與傳統(tǒng)工頻變壓器相比,高頻變壓器承受高頻非正弦電壓,電壓波形通常為雙極性方波。功率半導(dǎo)體器件的快速開通和關(guān)斷造成快速的電壓上升沿,即高d/d,高次電壓諧波也會(huì)疊加在瞬態(tài)繞組電壓上,導(dǎo)致高頻變壓器中出現(xiàn)多頻電應(yīng)力[9]。隨著碳化硅、氮化鎵等高頻大功率寬禁帶半導(dǎo)體器件的應(yīng)用,高頻變壓器中出現(xiàn)的暫態(tài)陡脈沖d/d可達(dá)10 kV/ms。高頻、高d/d脈沖與多諧波電壓給高頻變壓器絕緣可靠性帶來巨大挑戰(zhàn)。

電壓波形與頻率影響絕緣擊穿、局部放電、劣化與老化性能。趙義焜等研究了聚酰亞胺、nomex紙等四種絕緣材料的高頻擊穿特性,發(fā)現(xiàn)頻率從1 kHz上升到20 kHz時(shí),絕緣擊穿場強(qiáng)下降到原值的30%~40%;頻率為1~5 kHz時(shí),下降速率最快,達(dá)到了1.01 kV/kHz[10]。M. Khanali等對(duì)兩臺(tái)相同的變壓器分別施加工頻正弦電壓和高頻脈沖電壓 500 h進(jìn)行老化測試,并定期測量兩臺(tái)變壓器的局部放電值和絕緣電阻等介電特性。結(jié)果發(fā)現(xiàn)變壓器處于高頻脈沖電壓下時(shí)其絕緣性能明顯劣化[11]。郝春艷等采用針-板電極對(duì)交聯(lián)聚乙烯施加5~7kHz的交流電壓探究其電樹枝的生長規(guī)律,發(fā)現(xiàn)隨著頻率升高,電樹枝顏色加深,樹枝通道中碳素增多,電樹枝平均生長速率也顯著加快[12]。本課題組前期研究表明,溫度為80 ℃時(shí),環(huán)氧樹脂絕緣擊穿電場強(qiáng)度從0.5 kHz到10 kHz下降了34%[13-15]。

高頻諧波電壓影響絕緣的局部放電與擊穿性能。方田等對(duì)全膜電容器的絕緣介質(zhì)分別施加工頻電壓、工頻疊加3次諧波電壓、工頻疊加5次諧波電壓進(jìn)行短時(shí)擊穿實(shí)驗(yàn),發(fā)現(xiàn)擊穿電場強(qiáng)度隨著諧波次數(shù)增大而降低[16]。J.Kridsananont等研究了在疊加不同次數(shù)和不同百分比諧波電壓的高壓交流電壓下油紙絕緣的局部放電特性,發(fā)現(xiàn)局部放電起始電壓(Partial Discharge Inception Voltage, PDIV)隨諧波次數(shù)和百分比增大而減小,平均放電量隨諧波次數(shù)和百分比增大而增加[17]。研究發(fā)現(xiàn)4 kHz PWM諧波下低密度聚乙烯的擊穿強(qiáng)度比工頻50 Hz降低了33%[18]。Li Xiaonan等發(fā)現(xiàn)疊加諧波時(shí)油紙絕緣電壓的局部放電起始電壓和擊穿電壓均低于純交流電壓下的結(jié)果[19]。

高d/d脈沖電壓上升或下降沿變化快,并包含了多諧波電壓,從而影響絕緣電場、局部放電、老化與擊穿性能。梁成軍等采用棒-板電極和針-板電極,對(duì)不同電壓上升速率下的SF6間隙放電電壓測試發(fā)現(xiàn),間隙放電電壓隨電壓上升速率增大呈現(xiàn)先減小后增大的趨勢[20]。這主要是由于電壓上升速率增大之后“電暈穩(wěn)定化”作用減弱。Wang Peng等在1~5 kHz頻率范圍和50 ns~16 μs電壓上升時(shí)間范圍內(nèi)的重復(fù)方波電壓下對(duì)聚酰胺酰亞胺漆包線進(jìn)行了耐久性試驗(yàn),發(fā)現(xiàn)聚酰胺酰亞胺漆包線的絕緣壽命隨著頻率的增大和電壓上升的時(shí)間減?。磀/d增大)而減少[21]。于超凡等對(duì)變頻電極匝間絕緣材料試樣施加了電壓上升時(shí)間為60~1 000 ns的重復(fù)方波電壓進(jìn)行局部放電實(shí)驗(yàn),發(fā)現(xiàn)隨著電壓上升時(shí)間減小(即d/d增大),放電頻域能量呈現(xiàn)上升趨勢[22]。T. Hammarstrom等在4~40 μs電壓上升時(shí)間范圍的方波電壓下對(duì)帶有絕緣介質(zhì)電極的腔體進(jìn)行了局部放電實(shí)驗(yàn),發(fā)現(xiàn)上升時(shí)間越短(上升速率越大),局部放電振幅越大,且局部放電上升時(shí)間越少,會(huì)導(dǎo)致絕緣快速劣化并且使用壽命降低[23]。A. Komuro等對(duì)針-板電極施加了不同電壓上升速率的脈沖電壓,并拍攝了放電過程中的流注形狀,發(fā)現(xiàn)隨著電壓上升速率增大,放電電流、主流注的速度、流注通道的直徑和次流注的發(fā)射長度都會(huì)增大[24]。

研究指出,高頻脈沖電壓下局部放電是影響絕緣老化與劣化的關(guān)鍵[25-29]。已有研究指出,隨頻率增加局放數(shù)量和幅值先增加后減小,局部放電隨頻率增加出現(xiàn)一個(gè)拐點(diǎn)。研究認(rèn)為高頻時(shí)空間電荷極化需要一定時(shí)間,此時(shí)極化效應(yīng)不占主導(dǎo),而放電后的電荷復(fù)合作用增加,殘留同極性電荷,導(dǎo)致局放幅值和數(shù)量減小[30]。但電-熱效應(yīng)對(duì)空間電荷極化與復(fù)合過程并不清晰。局部放電起始電壓隨頻率增加而增加也難以解釋[28]。另外有研究表明,局部放電幅值和數(shù)量隨溫度增加而增加,而與頻率呈負(fù)相關(guān)[25]。頻率的增加一方面增加了單位時(shí)間雪崩放電積累破壞作用,另一方面增加了電荷注入/抽出次數(shù),產(chǎn)生熱電子,造成極化損耗與熱效應(yīng)。同時(shí),高頻電場在絕緣介質(zhì)中產(chǎn)生了較大位移電流,沖擊作用下導(dǎo)致局部過熱失效[25]。但很難解釋頻變對(duì)局部放電的作用弱于溫度的影響。重復(fù)脈沖電壓下固體絕緣老化機(jī)理復(fù)雜,研究表明高頻電場下分子基團(tuán)振動(dòng)造成的偶極子極化損耗增加了絕緣內(nèi)部熱效應(yīng),從而導(dǎo)致絕緣電-熱老化[31]。雖然高頻熱效應(yīng)是導(dǎo)致絕緣破壞的重要原因,但研究指出快速暫態(tài)脈沖下單純的介電損耗難以積累足夠能量導(dǎo)致分子鏈斷裂[14],因此,暫態(tài)脈沖下熱效應(yīng)可能不占主導(dǎo),高頻脈沖下電-熱耦合對(duì)極化、熱電子的本質(zhì)影響仍然需要深入研究。目前,這些研究大多集中在高頻正弦電壓下,高頻非正弦、d/d陡脈沖對(duì)絕緣電-熱耦合下極化損耗與沖擊能量積聚特性的影響仍鮮有報(bào)道。

本文針對(duì)高頻變壓器絕緣材料在高頻、高d/d電壓下的電-熱失效問題,開展d/d與方波電壓疊加下高頻變壓器固體絕緣環(huán)氧樹脂介質(zhì)損耗計(jì)算與分析,以及高d/d對(duì)絕緣沖擊能量損傷的影響。通過有限元物理場仿真研究獲得了高頻變壓器絕緣溫升與沖擊能量特性。研究結(jié)果對(duì)大容量高頻變壓器絕緣失效與可靠性具有重要意義。

1 高頻非正弦電壓下絕緣介質(zhì)損耗分析

1.1 絕緣介質(zhì)損耗

絕緣介質(zhì)在交變電壓下產(chǎn)生的損耗主要分為電導(dǎo)損耗和松弛損耗。電導(dǎo)損耗主要來源于介質(zhì)在外加電場作用下形成的傳導(dǎo)電流,松弛損耗主要來源于介質(zhì)在周期性極化過程中偶極子轉(zhuǎn)向和界面極化所消耗的能量。圖1為絕緣介質(zhì)等效電路與相量圖。在交變電場下,電子位移極化和離子位移極化時(shí)間短,不產(chǎn)生損耗。這兩種極化等效為純電容支路,相應(yīng)電流為瞬時(shí)充電電流i1。松弛極化可等效為電容和電阻串聯(lián),等效為吸收電流i2。而電導(dǎo)損耗對(duì)應(yīng)純電阻支路,相應(yīng)電流稱為電導(dǎo)電流i1。此時(shí)電流與外加電壓間有一個(gè)相位差,絕緣介電常數(shù)為一個(gè)復(fù)數(shù),介電常數(shù)虛部″可用Debye理論表示[13]為

式中,ε∞ 為高頻介電常數(shù),F(xiàn)/m;τ為松弛時(shí)間,s;γ為電導(dǎo)率,S/m;εs 為靜態(tài)介電常數(shù),F(xiàn)/m。

絕緣介質(zhì)的復(fù)介電常數(shù)可以通過介電譜測試獲得。圖2為環(huán)氧樹脂不同溫度下的復(fù)介電常數(shù)虛部結(jié)果??梢钥闯鼋殡姵?shù)虛部在玻璃化轉(zhuǎn)變溫度(g=120 ℃)之前變化不大,而在之后隨溫度迅速增加。介電常數(shù)虛部隨頻率的增大而增加。此結(jié)果可用于后續(xù)的介質(zhì)損耗計(jì)算。因此,與傳統(tǒng)介質(zhì)損耗計(jì)算不同,本文考慮了介電常數(shù)虛部與頻率和溫度的非線性關(guān)系進(jìn)行計(jì)算。

圖2 環(huán)氧樹脂絕緣介電常數(shù)虛部與頻率和溫度的關(guān)系

1.2 高頻方波下絕緣介質(zhì)損耗計(jì)算

當(dāng)絕緣介質(zhì)承受的外加電壓為非正弦電壓時(shí),絕緣介質(zhì)損耗為不同諧波電壓下介質(zhì)損耗的疊加,可以表示為

式中,0為真空電容(由0計(jì)算得到);V,RMS為所施加電壓第次諧波的有效值;為諧波次數(shù)。

由式(2)可知,諧波電壓幅值越大且諧波次數(shù)越多,絕緣介質(zhì)損耗也會(huì)隨之升高,基于這樣的結(jié)果,對(duì)d/d=0.5 kV/μs、d/d=1 kV/μs、d/d=5 kV/μs、d/d=10 kV/μs四種方波電壓進(jìn)行頻域分析。圖3為獲得的幅頻特性結(jié)果。隨著d/d增大,方波電壓的諧波次數(shù)和幅值增加,絕緣介質(zhì)損耗增大。為了將d/d對(duì)介質(zhì)損耗的影響進(jìn)行量化,通過模擬分析不同d/d方波電壓下絕緣介質(zhì)損耗。

圖3 不同dv/dt方波的諧波分量

圖4為采用低通濾波器的階躍響應(yīng)曲線模擬理想方波。低通濾波器的階躍響應(yīng)函數(shù)為

將式(3)與傅里葉級(jí)數(shù)結(jié)合可以計(jì)算得到方波

圖4 階躍函數(shù)擬合方波電壓

電壓的諧波分量有效值為

式中,sq為方波幅值,V;s為方波頻率,Hz。

將式(4)代入式(2)可得到高頻方波激勵(lì)下絕緣介質(zhì)損耗表達(dá)式為

將式(6)進(jìn)行等效變換,可以表示為

由于方波電壓的諧波次數(shù)為無窮大,因此可以將近似為

由于傅里葉分量的離散疊加值略高于連續(xù)積分值,因此式(7)和式(8)的近似計(jì)算過程會(huì)帶來一些誤差。為消除這種誤差,引入校正系數(shù)進(jìn)行校正,則可以表示為

式中,為校正系數(shù),=1.16。

由于高頻變壓器內(nèi)部電場分布并不均勻,因此式(9)可以改寫為

式中,sq為高頻變壓器內(nèi)部的絕緣材料所承受方波電場,V/m。

由此得到高頻變壓器內(nèi)部的絕緣介質(zhì)損耗功率密度為

1.3 dv/dt對(duì)絕緣介質(zhì)損耗的定量影響計(jì)算

由式(11)可以看出,絕緣介質(zhì)損耗密度除了與外加電壓頻率和所在區(qū)域電場強(qiáng)度有關(guān)之外,還和所使用的低通階躍響應(yīng)函數(shù)的截止頻率c有關(guān),而截止頻率c由電壓上升時(shí)間r計(jì)算得到。且r是與電壓上升速率d/d密切相關(guān)的值,即

設(shè)不同dv/dt方波電壓下電場強(qiáng)度均為1 kV/mm,可以計(jì)算得到不同dv/dt方波電壓下絕緣介質(zhì)損耗功率密度如圖5所示。結(jié)果表明,絕緣介質(zhì)損耗功率密度隨電壓上升速率dv/dt增加而增大,且呈現(xiàn)非線性增加。當(dāng)dv/dt為10 kV/μs時(shí),損耗功率最高達(dá)到了38.2 kW/m3,比dv/dt為0.5 kV/μs時(shí)的18.7 kW/m3增大了兩倍。

隨著大功率半導(dǎo)體器件的發(fā)展,高頻變壓器絕緣所承受的d/d電壓將不斷升高,相同電場強(qiáng)度下的介質(zhì)損耗將顯著增加。此外,本文計(jì)算所用的電場強(qiáng)度僅為1 kV/mm,隨著高頻變壓器進(jìn)一步小型化,結(jié)構(gòu)更加緊湊,其內(nèi)部電場強(qiáng)度進(jìn)一步增大,從而增加絕緣介質(zhì)高頻損耗。因此,高d/d陡脈沖與方波下絕緣介質(zhì)損耗明顯增加,從而增大了高頻變壓器溫升。

2 高頻變壓器物理場仿真模型

2.1 幾何模型與材料參數(shù)

本文所研究的對(duì)象是10 kW、10 kHz、1 000 V/ 750 V的高頻變壓器。變壓器磁心材料為鐵氧體,結(jié)構(gòu)為E型磁心。為了減小高頻趨膚與鄰近效應(yīng)的影響,高頻變壓器繞組采用多股利茲線。采用環(huán)氧樹脂作為高頻變壓器的主絕緣。環(huán)氧樹脂整體灌封于鐵心窗口,實(shí)現(xiàn)繞組絕緣。根據(jù)變化器設(shè)計(jì)尺寸和絕緣裕度,構(gòu)建高頻變壓器的幾何模型。圖6為10 kW、10 kHz、1000 V/750 V高頻變壓器二維幾何結(jié)構(gòu)與磁心窗口放大圖。為了保證計(jì)算精度,按照實(shí)際繞組而不是線圈域進(jìn)行構(gòu)建。同時(shí)為了保證計(jì)算效率,需要在仿真中進(jìn)行簡化,因此忽略導(dǎo)線絕緣層和變壓器外殼等部件,并且使用二維對(duì)稱幾何模型,可以提高仿真速度。表1列出了仿真中使用的材料參數(shù)。

圖6 高頻變壓器仿真幾何模型

表1 材料參數(shù)

Tab.1 Material specification

2.2 電磁-熱物理場計(jì)算方法

基于麥克斯韋方程,可以得到變壓器電磁場控制方程為

式中,為磁導(dǎo)率,H/m;為電導(dǎo)率,s/m。

在渦流區(qū)

在非渦流區(qū)

溫度場中,變壓器鐵心、繞組和環(huán)氧樹脂這些相互接觸的物體之間傳熱形式為熱傳導(dǎo)。而變壓器邊界處與空氣之間的傳熱形式為熱對(duì)流和熱輻射。

式中,λ、λ、λ分別為變壓器內(nèi)部不同方向的熱導(dǎo)率;C為恒壓熱容;為密度;為時(shí)間;為變壓器溫度;為熱源單位體積發(fā)熱功率。采用第三類邊界條件,即

式中,為流體的導(dǎo)熱系數(shù);為流體的表面熱系數(shù);a為流體溫度。

3 陡脈沖下絕緣電場與沖擊能量分析

3.1 絕緣電場分布計(jì)算

根據(jù)已有的高頻變壓器幾何模型,施加一次和二次側(cè)的輸入電壓,分別為1 000 V和750 V方波電壓,電壓的d/d為0.5 kV/μs。計(jì)算一個(gè)周期內(nèi)環(huán)氧樹脂絕緣的電場強(qiáng)度。圖7為計(jì)算獲得的隨時(shí)間變化的電場強(qiáng)度分布結(jié)果。結(jié)果表明電場集中在繞組端部到磁心之間的部位。電場分布隨時(shí)間而變化。在一個(gè)周期內(nèi),絕緣電場呈現(xiàn)周期性變化。在整個(gè)d/d時(shí)間段內(nèi),電場強(qiáng)度逐漸增加。當(dāng)= 0.002 ms時(shí)方波電壓上升到最大值,此時(shí)電場強(qiáng)度最大,約為1.2 kV/mm。最大電場集中在繞組端部絕緣處。

最大電場強(qiáng)度max隨時(shí)間變化規(guī)律(一個(gè)周期)如圖8所示。由圖8結(jié)果可知,最大電場值隨時(shí)間變化規(guī)律與電壓變化規(guī)律一致,則到達(dá)最大電場值所需時(shí)間與電壓上升速率d/d有關(guān),再分別輸入d/d為1 kV/μs、5 kV/μs、10 kV/μs的方波電壓進(jìn)行計(jì)算,得到最大電場強(qiáng)度所需時(shí)間隨d/d變化規(guī)律如圖9所示。隨著電壓上升速率增加,最大電場出現(xiàn)的時(shí)間更早。當(dāng)d/d=10 kV/ms時(shí),加壓0.1ms絕緣電場便可以達(dá)到最大值。這說明高d/d可在短時(shí)間內(nèi)導(dǎo)致絕緣局部電場畸變,從而引發(fā)絕緣劣化。

圖8 最大電場Emax值隨時(shí)間變化規(guī)律(一個(gè)周期)

圖9 最大電場值所需時(shí)間隨電壓上升速率變化規(guī)律

圖10為高頻變壓器溫度場仿真結(jié)果。高頻變壓器溫升主要來源于鐵心損耗、繞組損耗和絕緣介質(zhì)損耗。本文主要考慮環(huán)氧固體絕緣介電參數(shù),采用上述高頻方波下絕緣介質(zhì)損耗計(jì)算方法,計(jì)算獲得絕緣損耗,同時(shí)考慮d/d的影響(10 kV/μs)。鐵心與繞組損耗根據(jù)文獻(xiàn)[13]進(jìn)行計(jì)算。損耗作為熱源計(jì)算高頻變壓器的溫度分布。結(jié)果表明高頻變壓器熱點(diǎn)溫度出現(xiàn)在中部繞組中,可達(dá)56.4 ℃。鐵心外部溫度在45 ℃附近。

圖10 考慮介質(zhì)損耗的10 kW, 10 kHz高頻變壓器溫度場仿真結(jié)構(gòu)

3.2 絕緣沖擊能量積聚特性

由電場分布的結(jié)果可知,同樣施加電壓下,最大電場強(qiáng)度隨d/d的變化很小,但最大電場強(qiáng)度出現(xiàn)的時(shí)間縮短,表明高d/d時(shí)短時(shí)間絕緣會(huì)出現(xiàn)較強(qiáng)的電場畸變。為探究d/d對(duì)絕緣損傷的影響,根據(jù)絕緣在交變電場下的能量儲(chǔ)存和釋放特性,分析高d/d時(shí)絕緣暫態(tài)沖擊能量特性。此能量與絕緣介質(zhì)損耗和位移電流有關(guān)。當(dāng)絕緣兩端承受幅值為sq的方波電壓時(shí),可以將絕緣整體看作是一個(gè)復(fù)合電容器。考慮到絕緣材料內(nèi)部電場分布并不均勻,使用等效電容計(jì)算所儲(chǔ)存的能量無法反映電場集中處的情況,而這些電場畸變區(qū)域更容易發(fā)生擊穿失效。因此需要根據(jù)絕緣電場分布結(jié)果,計(jì)算每一點(diǎn)處的能量密度,即

式中,為絕緣材料介電常數(shù),F(xiàn)/m;sq為絕緣材料所承受方波電場強(qiáng)度值,V/m。

電場強(qiáng)度越大,相應(yīng)區(qū)域的能量密度就越大,能量密度與電場強(qiáng)度值的二次方成正比。介質(zhì)中的能量密度與電位移矢量有關(guān),取決于交變電場下介質(zhì)內(nèi)部位移電流。絕緣介質(zhì)中總的電流密度可表示為

電導(dǎo)率與及介電常數(shù)的關(guān)系為

因此,總的電流密度包含了電導(dǎo)電流與極化電流。式(20)中后一項(xiàng)會(huì)產(chǎn)生介質(zhì)損耗,前一項(xiàng)表示介質(zhì)存儲(chǔ)電場能量,不產(chǎn)生損耗。在電壓上升期間,每一點(diǎn)處的電場值也在持續(xù)增大,將電場能量密度對(duì)時(shí)間求導(dǎo)可以得到每一點(diǎn)處的功率密度,即

式中,d為絕緣材料中位移電流密度,A/m2。

提取上述一個(gè)周期內(nèi)繞組端部與磁心之間的最大強(qiáng)度計(jì)算結(jié)果,并采用式(20)進(jìn)行計(jì)算和處理,得到絕緣材料所承受的能量沖擊功率密度d在一個(gè)周期內(nèi)的變化規(guī)律,如圖11所示。

圖11 方波電壓下能量沖擊功率密度隨時(shí)間變化規(guī)律(dv/dt=0.5 kV/μs)

由圖11可知,在一個(gè)周期內(nèi)絕緣材料所承受的能量沖擊功率密度d在電壓上升和下降階段(d/d)變化最迅速,并且在電場最大時(shí)達(dá)到峰值,為17.1 MW/m3。電場能量密度遠(yuǎn)大于絕緣介質(zhì)損耗密度。與介質(zhì)損耗發(fā)熱不同,這部分能量會(huì)累積在絕緣材料內(nèi)部,在高頻電壓與高d/d不斷的作用下進(jìn)行充放電。在電壓上升、下降沿或極性反轉(zhuǎn)時(shí)能量積聚明顯增加。累積能量在短時(shí)間難以釋放,造成極化能量漲落較大,從而影響介質(zhì)材料松弛和損耗特性。此外,高d/d作用下絕緣介質(zhì)損耗增加,熱效應(yīng)會(huì)導(dǎo)致絕緣介電常數(shù)和電導(dǎo)增加,高頻分量和諧波作用下絕緣介質(zhì)熱效應(yīng)在短時(shí)間內(nèi)累積,以致局部放電、空間電荷積聚等現(xiàn)象發(fā)生,從而導(dǎo)致絕緣過早失效。

3.3 dv/dt對(duì)沖擊能量積聚的影響

將計(jì)算獲得的d/d=0.5 kV/μs、d/d=1 kV/μs、d/d=5 kV/μs、d/d=10 kV/μs四種方波電壓激勵(lì)下的電場分布作為已知量,根據(jù)式(21)計(jì)算獲得不同d/d方波電壓下沖擊能量密度d分布。圖12為不同d/d下繞組端部d的分布結(jié)果。結(jié)果表明,能量沖擊功率密度d的最大值集中在繞組端部,與電場分布一致,且隨著d/d增大d也明顯增大,表明d/d越大由電壓極性反轉(zhuǎn)帶來的沖擊能量積聚越大。通過能量沖擊功率密度d能夠反映出絕緣在高頻非正弦、高d/d電應(yīng)力下的絕緣性能。

圖12 不同dv/dt下繞組端部功率密度

4 高頻變壓器實(shí)驗(yàn)測試驗(yàn)證

通過雙向LLC諧振DC-DC功率變換器測量10 kW、10 kHz、1 kV/750 V的高頻變壓器,測試平臺(tái)如圖13所示。輸入電壓為1 kV,輸出電壓為750 V。該變換器的最大功率可達(dá)50 kW。制造的10 kW高頻變壓器位于中間。基于Si-IGBT模塊分別用于高壓側(cè)和低壓側(cè)。開關(guān)頻率與諧振頻率相等(10 kHz)。諧振式DC-DC變換器可實(shí)現(xiàn)軟開關(guān),使整個(gè)變換器效率最大化。為了實(shí)現(xiàn)簡單的控制,在一次側(cè)采用開環(huán)恒頻控制方法,在二次側(cè)采用非控制整流器。

圖14為高頻變壓器在10 kW運(yùn)行30 min后的溫度測試結(jié)果。高頻變壓器外部表面熱點(diǎn)溫度約為47.6 ℃,需要注意的是,鐵心和繞組熱點(diǎn)溫度高于外部溫度(47.6 ℃)。溫度測試結(jié)果與圖10中的溫度仿真(45 ℃)結(jié)果類似,誤差較小。因此,本文考慮d/d的高頻變壓器絕緣介質(zhì)損耗計(jì)算可得到驗(yàn)證。研究結(jié)果有利于指導(dǎo)大容量高頻變壓器分析與設(shè)計(jì)。

圖14 紅外熱像儀測試高頻變壓器溫度分布

5 結(jié)論

本文圍繞高頻、高d/d電壓下高頻變壓器固體絕緣環(huán)氧樹脂介質(zhì)損耗計(jì)算與分析,以及高d/d對(duì)絕緣沖擊能量影響展開研究。提出了考慮d/d的絕緣介質(zhì)損耗計(jì)算方法,獲得了不同d/d方波電壓對(duì)絕緣介質(zhì)損耗的影響。分析了高d/d對(duì)絕緣暫態(tài)沖擊能量密度的影響,獲得的主要結(jié)論如下:

1)根據(jù)絕緣介質(zhì)損耗機(jī)理,采用階躍函數(shù)模擬d/d脈沖,通過頻域疊加,推導(dǎo)得到了不同d/d方波電壓下絕緣介質(zhì)損耗計(jì)算公式。計(jì)算表明高d/d時(shí)絕緣介質(zhì)損耗可增加兩倍以上。

2)采用有限元物理場仿真發(fā)現(xiàn)絕緣電場畸變發(fā)生在繞組端部,且隨d/d周期變化很小,最大電場變化很小,但出現(xiàn)時(shí)間縮短。表明高d/d時(shí),絕緣內(nèi)部出現(xiàn)多次電場集中。

3)從絕緣材料在交變電場下能量儲(chǔ)存和釋放的角度分析了能量密度d。其與絕緣位移電流和電場畸變密切相關(guān)。計(jì)算表明d隨d/d增加而顯著增大,表明極性反轉(zhuǎn)時(shí)的沖擊能量顯著增加,d與介質(zhì)損耗協(xié)同作用導(dǎo)致高d/d下絕緣電-熱累積效應(yīng)增加,從而影響絕緣失效特性。10 kW、10 kHz高頻變壓器實(shí)驗(yàn)結(jié)果表明變壓器溫度與仿真一致,驗(yàn)證了本文的介質(zhì)損耗與能量沖擊的分析。

通過陡脈沖電場下介質(zhì)損耗與沖擊能量分析能夠反映出絕緣在高頻非正弦、高d/d電應(yīng)力下的絕緣性能,為大容量高頻變壓器絕緣可靠性與設(shè)計(jì)提供指導(dǎo)。

[1] 李凱, 趙爭鳴, 袁立強(qiáng), 等. 面向交直流混合配電系統(tǒng)的多端口電力電子變壓器研究綜述[J]. 高電壓技術(shù), 2021, 47(4): 1233-1250. Li Kai, Zhao Zhengming, Yuan Liqiang, et al. Overview on research of multi-port power electronic transformer oriented for AC/DC hybrid distribution grid[J]. High Voltage Engineering, 2021, 47(4): 1233-1250.

[2] 王威望, 劉瑩, 何杰峰, 等. 高壓大容量電力電子變壓器中高頻變壓器研究現(xiàn)狀和發(fā)展趨勢[J]. 高電壓技術(shù), 2020, 46(10): 3362-3373. Wang Weiwang, Liu Ying, He Jiefeng, et al. Research status and development of high frequency transformer used in high voltage and large capacity power electronic transformer[J]. High Voltage Engineering, 2020, 46(10): 3362-3373.

[3] 胡鈺杰, 李子欣, 羅龍, 等. 串聯(lián)諧振間接矩陣型電力電子變壓器高頻電流特性分析及開關(guān)頻率設(shè)計(jì)[J]. 電工技術(shù)學(xué)報(bào), 2022, 37(6): 1442-1454. Hu Yujie, Li Zixin, Luo Long, et al. Characteristic analysis of high-frequency-link current of series resonant indirect matrix type power electronics transformer and switching frequency design[J]. Transactions of China Electrotechnical Society, 2022, 37(6): 1442-1454.

[4] Wang Weiwang, He Jiefeng, Wang Xin, et al. Analysis of electric field stress and dielectric loss in insulation of magnetic component for cascaded power electronic transformer[C]//2020 4th International Conference on HVDC (HVDC), Xi'an, China, 2020: 1078-1083.

[5] Fauri M. Harmonic modelling of non-linear load by means of crossed frequency admittance matrix[J]. IEEE Transactions on Power Systems, 1997, 12(4): 1632-1638.

[6] Guillod T, F?rber R, Krismer F, et al. Computation and analysis of dielectric losses in MV power electronic converter insulation[C]//2016 IEEE Energy Conversion Congress and Exposition (ECCE), Milwaukee, WI, USA, 2017: 1-8.

[7] 韓繼業(yè), 李勇, 曹一家, 等. 基于模塊化多電平型固態(tài)變壓器的新型直流微網(wǎng)架構(gòu)及其控制策略[J]. 電網(wǎng)技術(shù), 2016, 40(3): 733-740. Han Jiye, Li Yong, Cao Yijia, et al. A new DC microgrid architecture based on MMC-SST and its control strategy[J]. Power System Technology, 2016, 40(3): 733-740.

[8] 王迎迎, 程紅. 應(yīng)用于功率變換器的多繞組高頻變壓器模型[J]. 電工技術(shù)學(xué)報(bào), 2021, 36(19): 4140-4147. Wang Yingying, Cheng Hong. Dual multi-winding high-frequency transformer equivalent circuit for power converter applications[J]. Transactions of China Electrotechnical Society, 2021, 36(19): 4140-4147.

[9] 汪濤, 駱仁松, 文繼峰, 等. 基于輔助繞組的高頻變壓器繞組損耗測量方法[J]. 電工技術(shù)學(xué)報(bào), 2022, 37(10): 2622-2630, 2655. Wang Tao, Luo Rensong, Wen Jifeng, et al. A measurement method of winding loss for high-frequency transformer based on auxiliary winding[J]. Transactions of China Electrotechnical Society, 2022, 37(10): 2622-2630, 2655.

[10] 趙義焜, 張國強(qiáng), 郭潤睿, 等. 高頻變壓器用耐高溫型匝間絕緣材料的擊穿特性[J]. 高電壓技術(shù), 2020, 46(2): 657-665. Zhao Yikun, Zhang Guoqiang, Guo Runrui, et al. Breakdown characteristics of high-temperature resistant turn-to-turn insulation materials for high-frequency transformers[J]. High Voltage Engineering, 2020, 46(2): 657-665.

[11] Khanali M, Jayaram S, Cheng J. Effects of voltages with high-frequency contents on the transformer insulation properties[C]//2013 IEEE Electrical Insulation Conference (EIC), Ottawa, ON, Canada, 2013: 235-238.

[12] 郝春艷, 郭軍科, 賀欣, 等. 高頻交流電壓下XLPE的電樹枝形貌特征及生長機(jī)理分析[J]. 絕緣材料, 2018, 51(3): 53-57. Hao Chunyan, Guo Junke, He Xin, et al. Analysis of morphology properties and growth mechanism of electrical trees in XLPE under high frequency AC voltage[J]. Insulating Materials, 2018, 51(3): 53-57.

[13] Wang Weiwang, He Jiefeng, Liu Ying, et al. Effects of spike voltages coupling with high d/dsquare wave on dielectric loss and electric-thermal field of high-frequency transformer[J]. IEEE Access, 2021, 9: 137733-137743.

[14] Wang Weiwang, Wang Xin, He Jiefeng, et al. Electric stress and dielectric breakdown characteristics under high-frequency voltages with multi-harmonics in a solid-state transformer[J]. International Journal of Electrical Power & Energy Systems, 2021, 129: 106861.

[15] Wang Xin, Wang Weiwang, Liu Y, et al. Influence of high frequency voltage with harmonic contents and temperature on breakdown of epoxy resin used in power electronic transformer[C]//International Symposium on Electrical Insulating Materials, Tokyo, Japan, 2020: 462-465.

[16] 方田, 李化, 黃想, 等. 工頻疊加諧波電壓下溫度對(duì)全膜電容器絕緣介質(zhì)擊穿特性的影響[J/OL]. 高壓電器, 2022: 1-7. (2022-04-25). https://kns.cnki.net/ kcms/detail/ 61.1127.TM.20220424.1602.002.html. Fang Tian, Li Hua, Huang Xiang, et al. Influence of temperature on the breakdown characteristics of foil-film capacitor insulation under AC superimposed harmonic voltage[J/OL]. High Voltage Apparatus, 2022: 1-7. (2022-04-25). https://kns.cnki.net/kcms/ detail/ 61.1127.TM.20220424.1602.002.html.

[17] Kridsananont J, Li Xining, Cui Yanjie, et al. Impact of harmonics on needle-plane partial discharge in oil-paper insulation[C]//2018 12th International Conference on the Properties and Applications of Dielectric Materials (ICPADM), Xi'an, China, 2018: 512-516.

[18] Suzuki H, Mukai S, Ohki Y, et al. Dielectric breakdown of low-density polyethylene under simulated inverter voltages[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 1997, 4(2): 238-240.

[19] Li Xiaonan, Wu Guangning, Yang Yan, et al. Partial discharge characteristics of oil-paper insulation for on-board traction transformers under superposed inter-harmonic AC voltages[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2020, 27(1): 240-248.

[20] 梁成軍, 劉軒東, 張玲俐, 等. 電壓上升率對(duì)GIS絕緣缺陷擊穿特性的影響[J]. 電網(wǎng)技術(shù), 2021, 45(3): 1195-1200.Liang Chengjun, Liu Xuandong, Zhang Lingli, et al. Effect of voltage rising rate on breakdown characteristics of GIS insulation defects[J]. Power System Technology, 2021, 45(3): 1195-1200.

[21] Wang Peng, Cavallini A, Montanari G C. The influence of repetitive square wave voltage parameters on enameled wire endurance[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2014, 21(3): 1276-1284.

[22] 于超凡, 王鵬, 馬世金, 等. 重復(fù)方波參數(shù)對(duì)變頻電機(jī)絕緣放電頻域能量分布影響研究[J]. 絕緣材料, 2022, 55(2): 78-83. Yu Chaofan, Wang Peng, Ma Shijin, et al. Influence of repetitive square wave voltage parameters on PD frequency domain energy distribution of inverter-fed motor insulation[J]. Insulating Materials, 2022, 55(2): 78-83.

[23] Hammarstrom T, Bengtsson T, Blennow J, et al. Evidence for changing PD properties at short voltage rise times[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2011, 18(5): 1686-1692.

[24] Komuro A, Ono R, Oda T. Effects of pulse voltage rise rate on velocity, diameter and radical production of an atmospheric-pressure streamer discharge[J]. Plasma Sources Science and Technology, 2013, 22(4): 045002.

[25] 韓帥, 李慶民, 劉偉杰, 等. 溫-頻耦合效應(yīng)對(duì)高頻固態(tài)變壓器絕緣局部放電特性的影響[J]. 電工技術(shù)學(xué)報(bào), 2015, 30(2): 204-210. Han Shuai, Li Qingmin, Liu Weijie, et al. Impacts of coupled temperature-frequency effects on partial discharge characteristics of high frequency solid state transformer insulation[J]. Transactions of China Electrotechnical Society, 2015, 30(2): 204-210.

[26] 劉濤,董國靜,李慶民,等. 高頻脈沖下電?熱應(yīng)力對(duì)聚酰亞胺絕緣壽命的耦合作用分析[J]. 高電壓技術(shù),2020, 46(7): 2504-2510.Liu Tao, Dong Guojing, Li Qingmin, et al. Coupling effect of electrical and thermal stresses on insulation life of polyimide under high frequency impulses[J]. High Voltage Engineering, 2020, 46(7): 2504-2510.

[27] 趙義焜, 張國強(qiáng), 韓冬, 等. 基于材料絕緣壽命的高頻變壓器繞組絕緣試驗(yàn)電壓確定方法[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(13): 2932-2939. Zhao Yikun, Zhang Guoqiang, Han Dong, et al. An insulation test method based on the dielectric insulation life for windings in high-frequency transformer[J]. Transactions of China Electrotechnical Society, 2020, 35(13): 2932-2939.

[28] Ghassemi M. Accelerated insulation aging due to fast, repetitive voltages: a review identifying challenges and future research needs[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2019, 26(5): 1558-1568.

[29] Zuo Zhou, Yao Chenguo, Dissado L A, et al. Simulation of electro-thermal ageing and breakdown in polymeric insulation under high frequency trapezoidal-wave pulses[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2017, 24(6): 3766-3775.

[30] Jiang Jun, Zhao Mingxin, Zhang Chaohai, et al. Partial discharge analysis in high-frequency transformer based on high-frequency current transducer[J]. Energies, 2018, 11(8): 1997.

[31] Niayesh K, Gockenbach E. On the aging mechanism of solid insulating materials exposed to repetitive high voltage pulses[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2014, 21(1): 304-310.

Dielectric Loss and Impact Energy Accumulation of High Frequency Transformer Insulation under Rapidly Repetitive Pulsed Voltages

Wang Weiwang Li Ruizhe He Jiefeng Zhang Xiaotong Li Shengtao

(State Key Laboratory of Electrical Insulation and Power Equipment Xi'an Jiaotong University Xi'an 710049 China)

As a key component, the high frequency transformer (HFT) plays an increasingly important role in voltage isolation and power transfer in solid state transformer (SST). Unlike the power frequency transformer, HFT faces the complex stresses of non-sinusoidal high frequency and high temperature. The fast and steep pulse voltage occurs on the rising and falling edges of a non-sinusoidal square wave voltage in HFTs. Its voltage change rate in the rise period (d/d) is higher than 10 kV/ms, which poses great challenges to the insulation of HFT. A High d/dvoltage contains a large number of high frequency harmonics, resulting in a significant increase in dielectric loss. The insulation suffers from cumulative impact stress at high frequency and high d/dvoltages, leading to local electric field distortion and partial discharge in the insulation. However, insulation degradation and breakdown is difficult to understand owing to the rapidly repetitive pulsed voltages. This paper focuses on the influence of d/don HFT dielectric losses and impact energy accumulation.

Firstly, the step response function was used to simulate a high d/dsquare-wave voltage. It combines with the decomposition and calculation of the dielectric loss. This paper used the complex dielectric parameter, including the¢and¢¢of epoxy resin. It depends on frequency and temperature, which can be used in the dielectric loss calculation. The actual square wave voltage can be simulated and treated by frequency domain response. Then a calculation method for dielectric losses with superposition of d/dand square-wave voltage was proposed. This equivalent technique simplified the complex multi-frequency dielectric loss calculation. The results indicated that the dielectric loss increases with increasing d/dof the square wave voltage, apparently enhancing the insulation heating. The dielectric loss reaches 38.2 kW/m3at 10 kV/μs, which is two times higher than at 0.5 kV/μs.

Secondly, a finite element simulation (FEM) of 10 kW, 10 kHz, 1 000 V/750 V HFT was performed, considering the dielectric parameters of epoxy resin and calculating the dielectric loss. The insulation electric field distortion occurs at the end of the winding, and its variation within the d/dperiod is small. However, the occurrence time ofmaxis short. It shows that the electric field distroation in the insulation occurs several times at high d/dperiod. For example, the insulation electric field can reach the maximum value at 0.1s (d/d=10 kV/ms). A high d/dcan causea local electric field distortion in a short time, resulting in insulation deterioration. In this case, the simulated maximum temperature on the winding can reach up to 56.4 ℃.

Finally, according to the energy storage of insulation under square-wave field, the impact power density of the insulationdat high d/dis studied. The FEM results showed that the electric field, displacement current and energy impact power density depend on d/d. The impact energy density induced by the square-wave voltage (dat 1 kV, 10 kHz) increases with d/d. The impact energy accumulation increases due to the high d/dcaused by the voltage polarity reversal. Thedcan reflect the insulation performance under high frequency non-sinusoidal voltages with high d/d. The experimental results of the 10 kW, 10 kHz, 1 000 V/750 V HFT were performed by a dual active bridge (DAB) platform. The temperature measured after 30 minutes operation under load agrees with the results of the FEM simulation. It can verify the dielectric loss calculation and thedanalysis. In combination with the dielectric loss, it indicates that the insulation is subjected to severe electro-thermal stresses, leading to insulation degradation at high d/dvoltages (during voltage polarity reversal). The research results provide guidelines for insulation failure and the design of high capacity HFTs.

High frequency transformer, d/ddielectric loss, insulation degradation, energy impact power density, displacement current

國家自然科學(xué)基金面上項(xiàng)目(52177025)、陜西省自然科學(xué)基礎(chǔ)研究計(jì)劃青年項(xiàng)目(2020JQ-045)和電力設(shè)備電氣絕緣國家重點(diǎn)實(shí)驗(yàn)室中青年基礎(chǔ)研究創(chuàng)新基金(EIPE21314)資助。

2022-09-26

2022-10-07

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

TM433

王威望 男,1987 年生,博士,副教授,研究方向?yàn)榻^緣介質(zhì)理論與應(yīng)用、高頻磁件設(shè)計(jì)與絕緣可靠性。E-mail:weiwwang@xjtu.edu.cn

李盛濤 男,1963 年生,博士,教授,研究方向?yàn)殡娊橘|(zhì)理論與應(yīng)用、新型電介質(zhì)材料與器件。E-mail:sli@mail.xjtu.edu.cn(通信作者)

(編輯 郭麗軍)

猜你喜歡
變壓器
10KV變壓器配電安裝柜中常見問題分析
理想變壓器的“三個(gè)不變”與“三個(gè)變”
超導(dǎo)變壓器技術(shù)展望
開關(guān)電源中高頻變壓器的設(shè)計(jì)
變壓器光纖測溫探頭的安裝固定
電子制作(2017年8期)2017-06-05 09:36:15
探究供配電變壓器的推廣運(yùn)行
電子制作(2017年8期)2017-06-05 09:36:15
變壓器運(yùn)行中的異?,F(xiàn)象與故障處理
電子制作(2017年7期)2017-06-05 09:36:14
變壓器差動(dòng)保護(hù)誤動(dòng)原因探討
電子制作(2017年10期)2017-04-18 07:23:16
一種不停電更換變壓器的帶電作業(yè)法
變壓器免維護(hù)吸濕器的開發(fā)與應(yīng)用
主站蜘蛛池模板: 国产真实二区一区在线亚洲| 国产精品吹潮在线观看中文| 一本无码在线观看| 欧美一区中文字幕| 激情网址在线观看| 欧美日韩国产在线人| 欧美、日韩、国产综合一区| 青青草国产在线视频| 国产哺乳奶水91在线播放| 999福利激情视频| yjizz国产在线视频网| 狠狠色丁香婷婷| 欧美成人午夜影院| 2021国产v亚洲v天堂无码| 54pao国产成人免费视频| 欧美A级V片在线观看| 欧美五月婷婷| 亚洲av无码片一区二区三区| 拍国产真实乱人偷精品| 美女被操91视频| 91精品专区国产盗摄| 亚洲V日韩V无码一区二区| 国产福利2021最新在线观看| 亚洲av片在线免费观看| 一区二区欧美日韩高清免费| 在线观看91香蕉国产免费| 欧美日韩导航| 亚洲一区二区三区香蕉| 国产麻豆另类AV| 免费国产高清视频| 女人18一级毛片免费观看| 国产精品露脸视频| 欧美在线黄| 亚洲第一成年网| 在线国产资源| 国产制服丝袜无码视频| 久久久成年黄色视频| 亚洲一区二区成人| 久久网欧美| 四虎综合网| 国产理论一区| 国产成人毛片| 最新亚洲人成无码网站欣赏网 | a级毛片免费看| 久久综合丝袜日本网| 中文字幕亚洲乱码熟女1区2区| 丁香综合在线| 精品国产乱码久久久久久一区二区| 亚洲精品成人福利在线电影| 国产菊爆视频在线观看| 国产在线精品99一区不卡| 国产91av在线| 国产自视频| 澳门av无码| 中文一级毛片| 精品国产一区91在线| 亚洲视频免费在线看| 最新国产麻豆aⅴ精品无| 精品国产香蕉在线播出| 国产美女免费| 成人欧美日韩| 色首页AV在线| 99在线视频精品| 在线a视频免费观看| 亚洲一区波多野结衣二区三区| 91毛片网| 91av成人日本不卡三区| 亚洲精品无码成人片在线观看| 女人一级毛片| 最新无码专区超级碰碰碰| 国产网友愉拍精品| 国产xx在线观看| 国产白浆视频| 亚洲男人的天堂久久精品| 亚洲国产日韩一区| 日本三区视频| 国外欧美一区另类中文字幕| 亚洲成a人在线观看| 亚洲成人播放| 美女免费精品高清毛片在线视| 欧美精品黑人粗大| 一级高清毛片免费a级高清毛片|