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

    熱聲發(fā)電系統(tǒng)最大聲電效率阻抗匹配分析

    2017-05-16 01:08:11夏加寬李文瑞
    電工技術(shù)學(xué)報 2017年9期
    關(guān)鍵詞:熱聲負載電阻阻抗匹配

    何 新 夏加寬 蘇 浩 李文瑞 張 健

    (1.沈陽工業(yè)大學(xué)電氣工程學(xué)院 沈陽 110870 2.沈陽師范大學(xué)數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院 沈陽 110034)

    ?

    熱聲發(fā)電系統(tǒng)最大聲電效率阻抗匹配分析

    何 新1,2夏加寬1蘇 浩1李文瑞1張 健1

    (1.沈陽工業(yè)大學(xué)電氣工程學(xué)院 沈陽 110870 2.沈陽師范大學(xué)數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院 沈陽 110034)

    為提高熱聲發(fā)電系統(tǒng)的聲電轉(zhuǎn)換效率,根據(jù)相似理論,結(jié)合聲學(xué)阻抗流源機理進行分析,提出一種永磁直線發(fā)電機彈簧機械阻抗和負載電阻抗的聲學(xué)阻抗匹配設(shè)計方法。該方法以提高捕獲聲功為約束條件,實現(xiàn)聲電轉(zhuǎn)換效率最大。進一步地,在彈簧機械阻抗和負載容抗聲學(xué)阻抗匹配時,給出發(fā)電機成為熱聲發(fā)動機聲學(xué)阻性負載、聲學(xué)感性負載和聲學(xué)容性負載的閾值規(guī)律。最后搭建實驗平臺,通過建模、求解及分析設(shè)計,驗證了發(fā)電機聲學(xué)阻抗匹配方法的有效性,為熱聲發(fā)電系統(tǒng)(TAEGS)的設(shè)計提供理論依據(jù)。

    熱聲發(fā)電系統(tǒng) 永磁直線發(fā)電機 熱聲發(fā)動機 聲學(xué)阻抗匹配 最大聲電效率

    0 引言

    熱聲發(fā)電系統(tǒng)(Thermoacoustic Electric Generation System,TAEGS)主要由熱聲發(fā)動機和直線發(fā)電機組成[1]。熱聲發(fā)動機由于回?zé)崞鲀啥舜嬖跍囟忍荻龋鶕?jù)熱聲效應(yīng)原理在活塞表面產(chǎn)生聲功。永磁直線發(fā)電機(Permanent Magnet Linear Generator,PMLG)嵌入到熱聲發(fā)動機的聲學(xué)流道內(nèi),其永磁體動子在活塞帶動下做直線往復(fù)運動,將捕獲的聲功轉(zhuǎn)換為電功[2,3]。熱聲發(fā)電技術(shù)可以將汽車尾氣、太陽能、工業(yè)余熱等低品質(zhì)熱源的熱能轉(zhuǎn)換成電能,可以如風(fēng)力發(fā)電技術(shù)[4-6],成為微電網(wǎng)[7]及能源互聯(lián)網(wǎng)的重要組成部分[8,9]。

    美國Los Alamos國家實驗室利用熱聲發(fā)電技術(shù),在熱聲發(fā)動機輸入500 W熱量時,發(fā)電機捕獲聲功143 W,輸出電功116 W,其聲電效率達到81%[10,11]。中國科學(xué)院理化技術(shù)研究所設(shè)計了聚能型百瓦級TAEGS和千瓦級太陽能TAEGS,其中百瓦級TAEGS在熱聲發(fā)動機輸入2 800 W熱量時,發(fā)電機捕獲聲功126 W,輸出電功97 W,其聲電效率為76%;千瓦級TAEGS在熱聲發(fā)動機輸入3 800 W熱量時,發(fā)電機捕獲聲功700 W,輸出電功481 W,其聲電效率為68.7%[12,13]。馬來西亞諾明漢大學(xué)設(shè)計了為偏遠山區(qū)提供電力支持的TAEGS,利用兩種不同燃料燃燒做飯火爐余熱而設(shè)計低成本TAEGS,在熱聲發(fā)動機輸入熱量2.54 kW時,發(fā)電機捕獲聲功460 W,其中由丙烷驅(qū)動火爐余熱的TAEGS產(chǎn)生約15 W電能,其聲電轉(zhuǎn)換效率為3.3%;而由燃燒木頭驅(qū)動火爐余熱的TAEGS產(chǎn)生約22.7 W電能,其聲電轉(zhuǎn)換效率為4.9%[14]。為了提高熱聲發(fā)電系統(tǒng)聲學(xué)阻抗匹配特性,S.Backhaus等[15]提出TAEGS聲學(xué)流源分析方法,優(yōu)化設(shè)計熱聲發(fā)動機使之與發(fā)電機匹配,提高系統(tǒng)聲電轉(zhuǎn)換效率。Z.B.Yu等[16,17]通過聲力電類比方法分析彈簧剛度系數(shù)和聲波壓力等對TAEGS輸出電功的影響。文獻[18-21]分析了負載電阻、活塞直徑、聲波壓力頻率和彈簧剛度系數(shù)對TAEGS的影響,并得出機械諧振和電路諧振可提高系統(tǒng)輸出電功。

    發(fā)電機嵌入到熱聲發(fā)動機的聲學(xué)流道后,由于氣體工質(zhì)聲學(xué)柔性特征,發(fā)電機聲學(xué)阻抗嵌入熱聲發(fā)動機作為其聲學(xué)負載直接參與工作循環(huán)。若直線發(fā)電機的聲學(xué)阻抗設(shè)計不合理,會導(dǎo)致耦合處聲波體積流率減小、聲波壓力與體積流率的聲學(xué)相位差增大,系統(tǒng)自激蕩起振難。兩機聲學(xué)阻抗不匹配使得發(fā)電機捕獲聲功能力減弱,最終影響TAEGS輸出的電功率和聲電轉(zhuǎn)換效率。而系統(tǒng)地分析直線發(fā)電機聲學(xué)阻抗與熱聲發(fā)動機聲學(xué)阻抗匹配機理,將發(fā)電機作為熱聲發(fā)動機聲學(xué)負載控制器,通過對發(fā)電機聲學(xué)阻抗進行優(yōu)化設(shè)計,提高發(fā)電機捕獲聲功,最終提高系統(tǒng)聲電轉(zhuǎn)換效率的發(fā)電機聲學(xué)阻抗設(shè)計,還沒有成熟的理論和方法。

    本文基于相似理論,在交變流體網(wǎng)絡(luò)里根據(jù)TAEGS聲學(xué)阻抗流源機理,采用類電路的相量法建立系統(tǒng)聲功捕獲和聲電效率的數(shù)學(xué)模型,給出以提高捕獲聲功為約束條件,實現(xiàn)TAEGS聲電效率最大的發(fā)電機聲學(xué)阻抗匹配設(shè)計方法,同時得到在發(fā)電機彈簧機械阻抗和外接負載容抗聲學(xué)阻抗匹配時,發(fā)電機成為熱聲發(fā)動機聲學(xué)阻性負載、聲學(xué)感性負載和聲學(xué)容性負載的負載電容設(shè)計的閾值規(guī)律。最后通過Matlab數(shù)值計算和實驗相結(jié)合,驗證了發(fā)電機聲學(xué)阻抗匹配方法的正確性,為TAEGS機械結(jié)構(gòu)和控制策略設(shè)計提供理論依據(jù)。

    1 TAEGS數(shù)學(xué)模型

    帶有可控電負載的TAEGS主要由熱聲發(fā)動機、發(fā)電機及其電負載組成,系統(tǒng)結(jié)構(gòu)如圖1所示。熱聲發(fā)動機由主冷卻器、回?zé)崞?、加熱器、熱緩沖管、次冷卻器、聲容管和慣性管等組成。PMLG動子軸通過機械彈簧連接活塞直接嵌入熱聲發(fā)動機的聲學(xué)流道。

    圖1 帶有可控電負載的TAEGS結(jié)構(gòu)Fig.1 Structure of TAEGS with controllable electric loads

    根據(jù)熱聲流源學(xué)說和聲學(xué)線性理論[16],回?zé)崞魇菬崧曅?yīng)的主要場所,且內(nèi)部粘滯損耗比熱馳豫損耗重要得多,因此回?zé)崞骱图訜崞鞫x為流阻和流源組合。主冷卻器和次冷卻器在聲學(xué)流道局部引入一個長度很長但橫截面積很小的聲學(xué)結(jié)構(gòu),相當(dāng)于引入聲阻和聲感。聲容管和慣性管為氣體介質(zhì)運行提供行波相位反饋通路,用聲感和聲容表示。熱緩沖管實現(xiàn)加熱器和次冷卻器熱隔離,表示成聲阻和聲容。帶有可控電負載的PMLG等效成聲阻、聲感和聲容。因此在氣體介質(zhì)運行小振幅線性分割假設(shè)條件下,熱聲發(fā)動機看成由低品質(zhì)熱源提供源動力,聲阻、聲感和聲容元件組成的聲學(xué)網(wǎng)絡(luò)系統(tǒng);而嵌入到聲學(xué)流道里的帶有電負載的發(fā)電機是其聲學(xué)負載。根據(jù)相似理論,采用類電學(xué)相量法,得到帶有可控電負載的TAEGS聲學(xué)阻抗流源模型如圖2所示。其中熱聲發(fā)動機聲學(xué)回路通過能量轉(zhuǎn)換器連接到發(fā)電機機械回路,機械回路再通過能量轉(zhuǎn)換器連接到帶有可變電負載的電氣回路。

    圖2 帶有可控電負載的TAEGS聲學(xué)阻抗流源模型Fig.2 Acoustic impedance current source model of TAEGS with controllable electric loads

    外接可控電負載的聲學(xué)阻抗與發(fā)電機電氣結(jié)構(gòu)的聲學(xué)阻抗互聯(lián);電氣結(jié)構(gòu)的聲學(xué)阻抗與其機械結(jié)構(gòu)的聲學(xué)阻抗互聯(lián);由可控電負載阻抗、電氣結(jié)構(gòu)阻抗和機械結(jié)構(gòu)阻抗組成的發(fā)電機聲學(xué)阻抗與熱聲發(fā)動機聲學(xué)阻抗互聯(lián),如圖2所示。根據(jù)聲學(xué)兩相流定理[22]和介觀熱力循環(huán)理論[23],與可控電負載阻抗互聯(lián)的發(fā)電機聲學(xué)阻抗變化會引起熱聲發(fā)動機聲學(xué)流道耦合處氣體微團的聲波壓力和體積流率之間相位發(fā)生變化,然后通過熱聲發(fā)動機機械結(jié)構(gòu)分布,氣體微團的能量傳遞會導(dǎo)致回?zé)崞骼餁怏w微團的聲波壓力和體積流率之間的相位也相應(yīng)發(fā)生變化,從而使熱聲發(fā)動機產(chǎn)生的聲功發(fā)生變化。

    通過聲力電類比方法,根據(jù)牛頓運動定律和基爾霍夫電壓定律,TAEGS運動方程和電壓方程為

    (1)

    結(jié)合特勒根定理和戴維南定理,求解式(1)可得定子繞組電流和兩機耦合處聲波體積流率為

    (2)

    (3)

    其中

    (4)

    (5)

    式中,ω=2πf。

    2 TAEGS聲學(xué)阻抗模型

    根據(jù)電路阻抗原理,帶有可控電負載的發(fā)電機的電阻抗為

    (6)

    帶有可控電負載的發(fā)電機的電阻抗模為

    (7)

    帶有可控電負載的發(fā)電機的電阻抗角為

    (8)

    根據(jù)相似理論,采用類電路的相量法,得到帶有電負載的直線發(fā)電機在兩機耦合處的聲學(xué)阻抗為

    Za=Ra+jXa

    (9)

    式中

    發(fā)電機在兩機耦合處的聲學(xué)阻抗模為

    (10)

    發(fā)電機在兩機耦合處的聲學(xué)阻抗角為

    ζa=arctan

    (11)

    3 TAEGS最大聲電效率聲學(xué)阻抗匹配設(shè)計

    嵌入熱聲發(fā)動機聲學(xué)流道的發(fā)電機捕獲聲功為

    (12)

    嵌入熱聲發(fā)動機聲學(xué)流道的發(fā)電機輸出電功為

    (13)

    則TAEGS的聲電效率為

    (14)

    從TAEGS聲電效率數(shù)學(xué)模型推導(dǎo)過程知,聲電效率與發(fā)電機捕獲聲功密切相關(guān)。由式(12)知,捕獲聲功既受發(fā)電機電氣參數(shù)(外接負載阻抗、定子繞組內(nèi)阻、定子繞組電感、機電常數(shù))和機械參數(shù)(彈簧剛度系數(shù)、阻尼系數(shù)、動子質(zhì)量、背腔容積、活塞面積和質(zhì)量)的影響,又受熱聲發(fā)動機參數(shù)(聲波壓力、氣體介質(zhì)運行頻率、氣體常數(shù)、氣體平均壓力)的影響。

    (15)

    (16)

    在提高捕獲聲功條件下,TAEGS聲電效率為

    (17)

    接下來通過優(yōu)化負載電阻實現(xiàn)聲電效率最大,即由dηae-better/dRL=0,得到負載電阻最優(yōu)值時TAEGS聲電效率最大為

    (18)

    因此以提高捕獲聲功為約束條件,以TAEGS聲電效率最大為優(yōu)化設(shè)計目標的發(fā)電機彈簧機械阻抗和負載電阻抗的聲學(xué)阻抗匹配設(shè)計方法,即

    (19)

    (20)

    (21)

    式中,f0為TAEGS穩(wěn)態(tài)時氣體運行頻率。

    當(dāng)機械彈簧阻抗采用最大聲電效率的聲學(xué)阻抗匹配方法進行設(shè)計時,改變外接容抗可改變發(fā)電機聲學(xué)阻抗,使其成為熱聲發(fā)動機的聲學(xué)容性負載、聲學(xué)感性負載和聲學(xué)阻性負載。當(dāng)電容按Ce>1/[(2πf0)2Lg]設(shè)計時,發(fā)電機呈現(xiàn)聲容特征,為熱聲發(fā)動機的聲學(xué)容性負載;當(dāng)電容按Ce=1/[(2πf0)2Lg]設(shè)計時,發(fā)電機呈現(xiàn)聲阻特征,為熱聲發(fā)動機的聲學(xué)阻性負載;當(dāng)電容按Ce<1/[(2πf0)2Lg]設(shè)計時,發(fā)電機呈現(xiàn)聲感特征,為熱聲發(fā)動機聲學(xué)感性負載。由此得到在機械彈簧實現(xiàn)聲學(xué)阻抗匹配時,改變電容使發(fā)電機成為熱聲發(fā)動機聲學(xué)容性負載、聲學(xué)感性負載和聲學(xué)阻性負載的閾值規(guī)律。

    4 TAEGS實例分析

    為驗證以提高捕獲聲功為約束條件、以TAEGS聲電效率最大為優(yōu)化設(shè)計目標的發(fā)電機聲學(xué)阻抗匹配設(shè)計方法,首先采用Matlab軟件進行數(shù)值計算,分析發(fā)電機聲學(xué)阻抗對捕獲聲功的影響。表1列出了數(shù)值計算過程中TAEGS平臺的主要參數(shù)。

    表1 熱聲發(fā)電系統(tǒng)平臺參數(shù)

    當(dāng)系統(tǒng)接入80 μF電容時,發(fā)電機在熱聲發(fā)動機的聲學(xué)流道里呈現(xiàn)聲容特征。數(shù)值計算分析彈簧機械阻抗和負載電阻變化對聲學(xué)容性發(fā)電機捕獲聲功的影響,如圖3所示,當(dāng)發(fā)電機作為熱聲發(fā)動機聲學(xué)容性負載時,連接動子和活塞的彈簧機械阻抗對捕獲聲功起關(guān)鍵作用。當(dāng)彈簧機械阻抗位于TAEGS機械諧振頻寬內(nèi)可迅速提高捕獲的聲功。發(fā)電機外接負載電阻變化也影響捕獲的聲功,但其影響權(quán)重遠低于彈簧機械阻抗。當(dāng)負載電容小于38 μF時,發(fā)電機呈現(xiàn)聲感特征,彈簧機械阻抗和外接負載電阻對聲功捕獲的影響與其作為聲學(xué)容性負載相似。

    圖3 彈簧機械阻抗和負載電阻對聲學(xué)容性發(fā)電機捕獲聲功的影響Fig.3 Influence on acoustic capacitive generator capturing acoustic power of spring stiffness and resistance

    當(dāng)負載電容按照最大聲電效率聲學(xué)阻抗匹配方法式(20)進行設(shè)計,即TAEGS接入電容38 μF,發(fā)電機呈現(xiàn)聲阻特征,數(shù)值計算分析彈簧機械阻抗和負載電阻變化對聲學(xué)阻性發(fā)電機捕獲聲功的影響,如圖4所示,發(fā)電機作為熱聲發(fā)動機聲學(xué)阻性負載時,彈簧機械阻抗和負載電阻對聲功捕獲都起到重要作用。當(dāng)彈簧機械阻抗位于TAEGS機械諧振頻寬內(nèi),系統(tǒng)捕獲聲功隨著負載電阻增大而增大。

    圖4 彈簧機械阻抗和負載電阻對聲學(xué)阻性發(fā)電機捕獲聲功的影響Fig.4 Influence on acoustic resistance generator capturing acoustic power of spring stiffness and resistance

    圖3和圖4說明發(fā)電機不僅是TAEGS聲電能量轉(zhuǎn)換裝置,而且直接參與其聲學(xué)循環(huán)。由機械阻抗和電氣阻抗所組成的發(fā)電機聲學(xué)阻抗必須與熱聲發(fā)動機的聲學(xué)阻抗相匹配,才能提高捕獲的聲功。無論發(fā)電機因為外接電負載改變而呈現(xiàn)聲容特征、聲感特征或者聲阻特征參與到熱聲發(fā)動機的聲學(xué)循環(huán)中,連接發(fā)電機動子和活塞的彈簧機械阻抗對捕獲聲功都起到關(guān)鍵作用。當(dāng)彈簧機械阻抗處于非諧振頻寬時,系統(tǒng)捕獲的聲功迅速降低。由于TAEGS的聲學(xué)特性,發(fā)電機捕獲熱聲發(fā)動機聲學(xué)流道里最大聲功的頻寬有限,因此彈簧機械阻抗設(shè)計尤其重要。

    按照TAEGS最大聲電效率聲學(xué)阻抗匹配設(shè)計方法進行彈簧機械阻抗的設(shè)計。根據(jù)式(19)彈簧剛度系數(shù)為285.5 kN/m,TAEGS在機械諧振頻寬內(nèi)運行。當(dāng)TAEGS機械結(jié)構(gòu)確定后,通過改變電負載改變發(fā)電機聲學(xué)阻抗,使其成為熱聲發(fā)動機聲學(xué)容性負載、聲學(xué)感性負載和聲學(xué)阻性負載。圖5和圖6分別為機械彈簧實現(xiàn)聲學(xué)阻抗匹配設(shè)計時,電負載變化對發(fā)電機聲學(xué)阻抗虛部和實部的影響。

    圖5 可控電容負載對發(fā)電機聲學(xué)阻抗虛部的影響Fig.5 Influence of the controllable capacitive load on acoustic impedance imaginary part

    圖6 可控電容負載對發(fā)電機聲學(xué)阻抗實部的影響Fig.6 Influence of the controllable capacitive load on acoustic impedance real part

    由圖5可知,當(dāng)外接電阻分別為8 Ω、20 Ω和50 Ω,電容為38 μF時,熱聲發(fā)動機的聲波壓力和體積流率的聲學(xué)阻抗角趨近于零,發(fā)電機成為熱聲發(fā)動機的聲學(xué)阻性負載;當(dāng)電容小于38 μF時,發(fā)電機為聲學(xué)感性負載;當(dāng)電容大于38 μF時,發(fā)電機為聲學(xué)容性負載。實際運行中熱聲發(fā)動機由于低品質(zhì)熱源加熱,不能處于理想的聲功輸出狀態(tài),可通過優(yōu)化設(shè)計發(fā)電機外接電容負載對兩機耦合處聲波壓力和體積流率的聲學(xué)阻抗角進行補償,提高TAEGS整機性能。

    由圖6可知,外接電容在10~80 μF范圍進行選值:當(dāng)外接電阻為8 Ω時,發(fā)電機聲學(xué)阻抗的實部范圍為(0.2×108~2.75×108) Pa·s/m3,其聲學(xué)阻抗的虛部范圍為(-1.45×108~1.45×108) Pa·s/m3;當(dāng)外接電阻為20 Ω時,發(fā)電機聲學(xué)阻抗的實部范圍為(0.2×108~1.45×108) Pa·s/m3,其聲學(xué)阻抗的虛部范圍為(-0.6×108~0.5×108) Pa·s/m3;當(dāng)外接電阻為50 Ω時,發(fā)電機聲學(xué)阻抗的實部范圍為(0.2×108~0.6×108) Pa·s/m3,其聲學(xué)阻抗的虛部范圍為(-0.2×108~0.2×108) Pa·s/m3。

    搭建熱聲發(fā)電系統(tǒng)實驗平臺,其原理如圖7所示,實驗平臺主要參數(shù)見表1。

    圖7 熱聲發(fā)電系統(tǒng)實驗平臺原理Fig.7 Principle diagram of TAEGS experimental platform

    熱聲發(fā)電系統(tǒng)實驗平臺主要由熱聲發(fā)動機、直線發(fā)電機、檢測系統(tǒng)和控制系統(tǒng)組成,如圖7所示。實驗過程中,熱聲發(fā)動機內(nèi)充入5.5 MPa的氦氣,采用丁烷氣體燃燒模擬低品質(zhì)熱源為TAEGS供熱。檢測系統(tǒng)主要由熱電偶、壓力傳感器、流量傳感器和數(shù)字示波器組成,它們分別采集熱聲發(fā)動機回?zé)崞鳒囟萒1、加熱器溫度T2和水冷卻器溫度T3、丁烷體積流速V1、活塞表面和背腔的聲波壓力Fa1和Fa2、發(fā)電機背腔的平均壓力F0、活塞的位移xg、定子繞組電流ig和電壓Ug等參數(shù)。檢測系統(tǒng)將檢測信號提供給控制系統(tǒng),控制系統(tǒng)控制丁烷氣體流速、水泵電機速度、直線發(fā)電機的起振和可控電負載的切換等功能。實驗中根據(jù)文獻[24]所提出的方法進行聲功的測量,根據(jù)文獻[25]所提出的方法進行TAEGS起振控制。在系統(tǒng)起振的初始階段,控制直線發(fā)電機電動運行拖動熱聲發(fā)動機產(chǎn)生振蕩。熱聲發(fā)動機在丁烷加熱條件下,回?zé)崞骱庾约ふ袷幃a(chǎn)生聲功,當(dāng)熱聲發(fā)動機和直線發(fā)電機達到諧振運行狀態(tài)后,控制器切換成整流狀態(tài),電機由電動拖動模式轉(zhuǎn)換成熱聲發(fā)動機的聲學(xué)負載發(fā)電模式運行。當(dāng)加熱到2 159W時,TAEGS進入穩(wěn)態(tài)運行狀態(tài)。采用燈泡作為可視負載觀測系統(tǒng)輸出電功狀態(tài),同時采用示波器實時監(jiān)測系統(tǒng)輸出的電參數(shù)。

    實驗平臺根據(jù)式(19)設(shè)計彈簧剛度系數(shù)為285.5kN/m,根據(jù)式(20)設(shè)計負載電容為38μF。圖8 給出了彈簧機械阻抗匹配和負載容抗匹配實驗條件下,發(fā)電機捕獲的聲功和系統(tǒng)輸出電功及聲電轉(zhuǎn)換效率隨著外接負載電阻變化的影響。

    實驗過程中考慮到定子繞組電流和動子位移的限制,電阻設(shè)定在16~50Ω之間進行測試。當(dāng)外接電阻16Ω時,捕獲聲功、輸出電功和聲電效率的仿真結(jié)果分別為124W、100W和80.5%;當(dāng)外接電阻50Ω時,捕獲聲功、輸出電功和聲電效率的仿真結(jié)果分別為305W、227W和74.5%。仿真計算中系統(tǒng)聲電效率最大值為80.8%,此時外接電阻為20Ω,捕獲聲功和輸出電功的仿真結(jié)果分別為148W和120W。在相同條件下,當(dāng)外接電阻16Ω時,實驗測得捕獲聲功、輸出電功和聲電效率分別為120W、91W和75%;當(dāng)外接電阻50Ω時,實驗測得捕獲聲功、輸出電功和聲電效率分別為278W、189W和68.2%。實驗中得到當(dāng)捕獲聲功和輸出電功分別為142W和108W時,系統(tǒng)聲電效率最大,其值為75.8%,此時對應(yīng)的電阻為20Ω,恰好與最大聲電效率聲學(xué)阻抗匹配方法中式(21)所確定的外接電阻優(yōu)化設(shè)計值相同。

    圖8 負載電阻對聲功、電功和聲電效率的影響Fig.8 Influence of resistance on acoustic power,electric power and acoustic-to-electric efficiency

    由于仿真計算時只考慮了壓力為小振幅振蕩、熱聲發(fā)動機氣體工作為層流狀態(tài),忽略了氣體介質(zhì)的紊流效應(yīng)及其導(dǎo)致的粘性耗散;另外由于熱聲發(fā)電系統(tǒng)的聲學(xué)柔性,外接電阻及加熱溫度變化導(dǎo)致兩機耦合處聲波壓力幅值發(fā)生變化;再有發(fā)電機端部漏磁及測量誤差等因素,實驗中測量得到的捕獲聲功、輸出電功和聲電效率數(shù)據(jù)都低于仿真計算結(jié)果。但是系統(tǒng)聲功捕獲和聲電效率變化趨勢相同,且聲電效率出現(xiàn)最大值所對應(yīng)的阻值相同。

    5 結(jié)論

    結(jié)合TAEGS聲學(xué)特性,提出以提高聲功捕獲為約束條件,實現(xiàn)最大聲電轉(zhuǎn)換效率的發(fā)電機聲學(xué)阻抗匹配優(yōu)化設(shè)計方法,并搭建TAEGS實驗平臺驗證方法的正確性。具體結(jié)論如下:

    1)彈簧剛度系數(shù)變化會導(dǎo)致發(fā)電機由機械阻抗組成的聲學(xué)阻抗特性發(fā)生變化。

    2)負載電阻抗變化會導(dǎo)致發(fā)電機由負載電阻抗組成的聲學(xué)阻抗特性發(fā)生變化。

    3)優(yōu)化設(shè)計彈簧剛度系數(shù)、負載電阻抗可以提高兩機耦合處聲學(xué)阻抗匹配特性,提高聲電轉(zhuǎn)換效率。

    4)彈簧機械阻抗和負載電阻抗直接影響TAEGS輸出電功率和熱聲發(fā)動機輸出聲功率,不同的TAEGS其機械彈簧、負載電阻和負載電容存在不同的優(yōu)化閾值,進行實際設(shè)計時必須結(jié)合設(shè)備具體參數(shù)進行其雙功率的優(yōu)化設(shè)計。

    本文所提聲學(xué)阻抗匹配設(shè)計方法可為TAEGS機械結(jié)構(gòu)和控制策略設(shè)計提供理論依據(jù)和技術(shù)支撐。

    [1]WuZhanghua,ManMan,LuoEC,etal.Experimentalinvestigationofa500Wtravelingwavethermoacousticelectricitygenerator[J].ChineseScienceBulletin,2011,56(19):1975-1977.

    [2] 夏加寬,何新,賈海源,等.混合磁路熱聲直線發(fā)電機:中國,201010539540.9[P].2012-11-14.

    [3] 夏加寬,李文瑞,何新,等.各向異性橫向磁通永磁直線電機穿片漏磁分析[J].電工技術(shù)學(xué)報,2015,30(14):518-524.XiaJiakuan,LiWenrui,HeXin,etal.Analysisofanisotropicpermeateleakagefluxontransversefluxpermanentmagnetlinearmachine[J].TransactionsofChinaElectrotechnicalSociety,2015,30(14):518-524.

    [4] 楊睿茜,陳紅坤,陶玉波,等.雙饋式風(fēng)機控制器參數(shù)整定的實用方法[J].電力系統(tǒng)保護與控制,2015,43(2):63-69.YangRuixi,ChenHongkun,TaoYubo,etal.Apracticalmethodondoubly-fedinductiongeneratorcontrollerparameterdetermination[J].PowerSystemProtectionandControl,2015,43(2):63-69.

    [5] 王愷,關(guān)少卿,汪令祥,等.基于模糊信息?;妥钚《酥С窒蛄繖C的風(fēng)電功率聯(lián)合預(yù)測建模[J].電力系統(tǒng)保護與控制,2015,43(3):26-32.WangKai,GuanShaoqing,WangLingyang,etal.Acombinedforcastingmodelforwindpowerpredictionbasedonfuzzyinformationgranulationandleastsquaressupportvectormachine[J].PowerSystemProtectionandControl,2015,43(3):26-32.

    [6] 馮江霞,梁軍,馮益坤.基于風(fēng)電功率min級分量波動特性的風(fēng)電場儲能容量優(yōu)化計算[J].電力系統(tǒng)保護與控制,2015,43(3):1-6.FengJiangxia,LiangJun,F(xiàn)engYikun.Optimizationcalculationofwindfarmstoragecapacitybasedonminute-scalewindpowercharacteristic[J].PowerSystemProtectionandControl,2015,43(3):1-6.

    [7] 張建華,于雷,劉念,等.含風(fēng)/光/柴/蓄及海水淡化負荷的微電網(wǎng)容量優(yōu)化配置[J].電工技術(shù)學(xué)報,2015,30(2):102-112.ZhangJianhua,YuLei,LiuNian,etal.CapacityconfigurationoptimizationforIslandmicrogridwithwind/photovoltaic/diesel/storageandseawaterdesalinationload[J].TransactionsofChinaElectrotechnicalSociety,2015,30(2):102-112.

    [8] 趙海,蔡巍,王進法,等.能源互聯(lián)網(wǎng)架構(gòu)設(shè)計與拓撲模型[J].電工技術(shù)學(xué)報,2015,30(11):30-36.ZhaoHai,CaiWei,WangJinfa,etal.Analysisofanisotropicpermeateleakagefluxontransversefluxpermanentmagnetlinearmachine[J].TransactionsofChinaElectrotechnicalSociety,2015,30(11):30-36.

    [9] 余曉丹,徐憲東,陳碩翼,等.綜合能源系統(tǒng)與能源互聯(lián)網(wǎng)簡述[J].電工技術(shù)學(xué)報,2016,31(1):1-13.YuXiaodan,XuXiandong,ChenShuoyi,etal.Abriefreviewtointegratedenergysystemandenergyinternet[J].TransactionsofChinaElectrotechnicalSociety,2016,31(1):1-13.

    [10]BackhausS,SwiftGW.Athermoacoustic-stirlingheatengine:detailedstudy[J].JournaloftheAcousticalSocietyofAmerica,2000,107(6):3148-3166.

    [11]PetachM,TwardE,BackhausS.Designofahighefficiencypowersource(HEPS)basedonthermoacoustictechnology[R].Cleveland,OH:NASAGlennResearchCenter,2004.

    [12]WuZhanghua,ZhangLimin,DaiWei,etal.Investigationona1kWtravelingwavethermoacousticelectricalgenerator[J].AppliedEnergy,2014,124(3):140-147.

    [13]WuZhanghua,ManMan,DaiWei,etal.Asolar-poweredtraveling-wavethermoacousticelectricitygenerator[J].SolarEnergy,2012,86(9):2376-2382.

    [14]BaimanC,AbdallaAY,PaulHR,etal.Developmentandassessmentofthermoacousticgeneratorsoperatingbywasteheatfromcookingstove[J].2012,4(12):894-902.

    [15]BackhausS,TwardE,PetachM.Traveling-wavethermoacousticelectricgenerator[J].AppliedPhysicsLetters,2004,85(6):1085-1087.

    [16]YuZB,SaechanP,JaworskiAJ.Amethodofcharacteringperformanceofaudioloudspeakersforlinearalternatorapplicationsinlow-costthermoacousticelectricitygenerators[J].AppliedAcoustics,2011,72(5):260-267.

    [17]YuZB,JaworskiAJ,BackhausS.Travelling-wavethermoacousticelectricitygeneratorusinganultra-compliantalternatorforutilizationoflow-gradethermalenergy[J].AppliedEnergy,2012(99):135-145.

    [18]邱利民,趙益濤,王凱,等.百瓦級熱聲斯特林發(fā)電系統(tǒng)輸出特性研究[J].中國電機工程學(xué)報,2013,33(8):101-106.QiuLimin,ZhaoYitao,WangKai,etal.Studyonoutputcharacteristicsof100W-classthermoacoustic-stirlingelectricgenerators[J].ProceedingsoftheCSEE,2013,33(8):101-106.

    [19]夏加寬,李文瑞,何新,等.熱聲發(fā)電系統(tǒng)最大電能輸出阻抗匹配及機械諧振特性分析[J].中國電機工程學(xué)報,2015,35(22):5908-5914.XiaJiakuan,LiWenrui,HeXin,etal.Analysisonimpedancematchingandmechanicalresonancecharacteristicsofmaximumelectricpoweroutputforthermoacousticelectricgenerationsystems[J].Proceed-ingsoftheCSEE,2015,35(22):5908-5914.

    [20]王凱,孫大明,郭軼楠,等.行波熱聲發(fā)電系統(tǒng)的諧振特性研究[J].中國電機工程學(xué)報,2014,34(9):1423-1428.WangKai,SunDaming,GuoYinan,etal.Studyontheresonancecharacteristicsofatraveling-wavethermoacousticelectricgenerationsystem[J].Proceed-ingsoftheCSEE,2014,34(9):1423-1428.

    [21]吳張華,羅二倉,戴巍.熱聲發(fā)電之直線發(fā)電機的理論研究[J].太陽能學(xué)報,2008,29(4):493-497.WuZhanghua,LuoErcang,DaiWei.Theoreticalinvestigationonlinearalternatorinthermoacousticpowergenerationsystem[J].ActaEnergiaeSolarisSinica,2008,29(4):493-497.

    [22]李正宇.熱動力學(xué)相似性——受熱流動系統(tǒng)的時間特征參量研究[D].武漢:華中科技大學(xué),2004.

    [23]羅二倉,戴巍,吳張華,等.交變流動熱機的介觀熱力循環(huán)理論,第二部分:發(fā)動機回?zé)崞鞯慕橛^熱力循環(huán)模型及分析[J].低溫工程,2004(2):1-10.LuoErcang,DaiWei,WuZhanghua,etal.Meso-scopethermodynamictheoryforcyclicflowengines,partⅡ,meso-scopethermodynamicmodelforacyclicflowregeneratorofprimemovers[J].Cryogenics,2004(2):1-10.

    [24]趙瑞東,吳張華,戴巍,等.回?zé)崞鳠峁D(zhuǎn)換特性的實驗研究[J].太陽能學(xué)報,2014,35(11):2171-2176.ZhaoRuidong,WuZhanghua,DaiWei,etal.Studyonthethermoacousticconversioncharacteristicoftheregenerator[J].ActaEnergiaeSolarisSinica,2014,35(11):2171-2176.

    [25]夏加寬,何新,趙志強,等.熱聲熱聲發(fā)動機起振控制方法:中國,201010581563.6[P].2012-01-14.

    (編輯 張洪霞)

    Analysis on Impedance Matching of Thermoacoustic Electric Generation System for the Maximum Acoustic-to-Electric Efficiency

    HeXin1,2XiaJiakuan1SuHao1LiWenrui1ZhangJian1

    (1.School of Electrical Engineering Shenyang University of Technology Shenyang 110870 China 2.School of Mathematics and System Science Shenyang Normal University Shenyang 110034 China)

    In order to improve the acoustic-to-electric efficiency for thermoacoustic electric generation system (TAEGS), an acoustic impedance matching design on permanent magnet linear generator (PMLG) is proposed. Based on the similarity theory, analysis and calculations are carried out to study the spring mechanical and load impedance, combining with acoustic impedance current source mechanism. The method obtains the maximum acoustic-to-electric efficiency with the constraints to improve capturing acoustic power for PMLG. Furthermore, the threshold law of external capacitance is achieved for PMLG acting as acoustic impedance, acoustic capacitive and acoustic inductive load on thermoacoustic engine, under the acoustic impedance matching of the spring mechanical impedance and capacitance reactance. Finally, the experimental platform is built. The results show that the effectiveness of the proposed acoustic impedance matching design, providing the theoretical and technical basis for TAEGS design.

    Thermoacoustic electric generation system,permanent magnet linear generator,thermoacoustic engine,acoustic impedance matching,the maximum acoustic-to-electric efficiency

    國家自然科學(xué)基金(51377108)和遼寧省博士啟動基金(20141073)資助項目。

    2016-01-20 改稿日期2016-06-28

    TM615

    何 新 女,1974年生,博士,副教授,研究方向為熱聲發(fā)電系統(tǒng)設(shè)計及其控制等。

    E-mail:hexin9898@163.com(通信作者)

    夏加寬 男,1962年生,博士,教授,博士生導(dǎo)師,研究方向為直線電機及其控制等。

    E-mail:xiajk_mail@163.com

    猜你喜歡
    熱聲負載電阻阻抗匹配
    熱驅(qū)動熱聲制冷技術(shù)發(fā)展現(xiàn)狀與展望
    基于LCC-S補償拓樸MCR-WPT系統(tǒng)的寬范圍阻抗匹配研究
    非對稱雙級環(huán)路行波熱聲熱機的實驗研究*
    新能源進展(2022年2期)2022-05-06 04:43:38
    太陽能電池板加DC-DC前后負載功率隨負載電阻變化規(guī)律的實驗研究
    熱聲效應(yīng)及其應(yīng)用研究進展
    多階阻抗匹配與寬頻帶天線的設(shè)計分析
    基于E類功率放大器的非接觸感應(yīng)耦合電能傳輸系統(tǒng)
    高頻熱聲發(fā)動機的聲耦合特性
    巧解原線圈有負載電阻的變壓器問題
    一種阻抗匹配層的截斷邊界
    日本爱情动作片www.在线观看 | 免费看光身美女| 免费无遮挡裸体视频| av中文乱码字幕在线| 美女黄网站色视频| 国产精品久久久久久av不卡| 成人二区视频| 真人做人爱边吃奶动态| 动漫黄色视频在线观看| 午夜a级毛片| 制服丝袜大香蕉在线| 国产欧美日韩一区二区精品| 真实男女啪啪啪动态图| 熟女电影av网| 国产伦精品一区二区三区视频9| 在线观看美女被高潮喷水网站| 国产美女午夜福利| 亚洲欧美日韩高清在线视频| 不卡视频在线观看欧美| 直男gayav资源| 精品日产1卡2卡| 99热网站在线观看| 日本一本二区三区精品| 久久人人爽人人爽人人片va| 欧美高清成人免费视频www| 哪里可以看免费的av片| 天堂动漫精品| 日本 欧美在线| 亚洲欧美精品综合久久99| 国产亚洲精品av在线| 国产爱豆传媒在线观看| 久久人人爽人人爽人人片va| 熟女人妻精品中文字幕| 免费一级毛片在线播放高清视频| 精品人妻偷拍中文字幕| 成人特级av手机在线观看| 日日撸夜夜添| 狠狠狠狠99中文字幕| 成年免费大片在线观看| 国产精品亚洲美女久久久| 国产乱人视频| 久久精品影院6| 国产亚洲精品久久久com| 久久精品国产亚洲av天美| 亚洲最大成人中文| 亚洲美女视频黄频| 成人无遮挡网站| 亚洲男人的天堂狠狠| 国产精品国产高清国产av| 国产极品精品免费视频能看的| 又紧又爽又黄一区二区| 亚洲美女搞黄在线观看 | 国产精品99久久久久久久久| 午夜亚洲福利在线播放| 国产熟女欧美一区二区| 久久久色成人| 久久久色成人| 亚洲欧美日韩高清在线视频| 成人av一区二区三区在线看| 欧美性猛交黑人性爽| 乱系列少妇在线播放| 国产色婷婷99| 悠悠久久av| 波多野结衣高清作品| 1024手机看黄色片| 日韩在线高清观看一区二区三区 | 国模一区二区三区四区视频| 乱码一卡2卡4卡精品| 国产高清不卡午夜福利| 五月玫瑰六月丁香| 高清毛片免费观看视频网站| 欧美国产日韩亚洲一区| 不卡一级毛片| 麻豆国产97在线/欧美| 欧美xxxx黑人xx丫x性爽| 亚洲最大成人手机在线| 亚洲黑人精品在线| 成人永久免费在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产高清在线一区二区三| 欧美国产日韩亚洲一区| 日韩高清综合在线| 亚洲va在线va天堂va国产| 蜜桃久久精品国产亚洲av| 国产精品1区2区在线观看.| 男女做爰动态图高潮gif福利片| 99久久无色码亚洲精品果冻| 国产大屁股一区二区在线视频| 国产老妇女一区| 免费av毛片视频| 性色avwww在线观看| 观看免费一级毛片| 国内少妇人妻偷人精品xxx网站| 麻豆国产97在线/欧美| 午夜激情欧美在线| 日韩国内少妇激情av| 午夜福利欧美成人| 亚洲va在线va天堂va国产| 久久久久久久亚洲中文字幕| 国内少妇人妻偷人精品xxx网站| 欧美日韩中文字幕国产精品一区二区三区| 乱系列少妇在线播放| 成人亚洲精品av一区二区| 干丝袜人妻中文字幕| 久久精品综合一区二区三区| av.在线天堂| 一本精品99久久精品77| 欧美成人性av电影在线观看| 噜噜噜噜噜久久久久久91| 美女免费视频网站| 国产精品无大码| 国产精品久久久久久亚洲av鲁大| 午夜影院日韩av| 精品久久久久久久久亚洲 | 欧美xxxx性猛交bbbb| 国产麻豆成人av免费视频| 免费在线观看成人毛片| 日本黄色片子视频| 国产精品久久久久久久久免| 成年女人毛片免费观看观看9| 亚洲avbb在线观看| 国产亚洲av嫩草精品影院| 婷婷六月久久综合丁香| 日本a在线网址| 欧美丝袜亚洲另类 | 大型黄色视频在线免费观看| 少妇丰满av| 国产一区二区在线av高清观看| 亚洲最大成人手机在线| 九九爱精品视频在线观看| 亚洲午夜理论影院| 欧美不卡视频在线免费观看| 午夜福利高清视频| 91在线精品国自产拍蜜月| 在线播放国产精品三级| 午夜福利高清视频| 日日干狠狠操夜夜爽| 91麻豆av在线| 国产精品亚洲一级av第二区| a级毛片免费高清观看在线播放| 国产单亲对白刺激| 夜夜看夜夜爽夜夜摸| 国产一区二区亚洲精品在线观看| 亚洲人成网站在线播| 又紧又爽又黄一区二区| 久久精品国产自在天天线| 老熟妇仑乱视频hdxx| 精品一区二区三区视频在线观看免费| 村上凉子中文字幕在线| 男女边吃奶边做爰视频| 亚洲熟妇熟女久久| 亚洲精品粉嫩美女一区| 国产真实乱freesex| 午夜日韩欧美国产| 久久久久久伊人网av| 男人和女人高潮做爰伦理| 国产精品人妻久久久久久| 日韩中字成人| 成人无遮挡网站| 十八禁国产超污无遮挡网站| 欧美绝顶高潮抽搐喷水| 国内精品一区二区在线观看| 国产在线男女| 天天一区二区日本电影三级| 亚洲综合色惰| 亚洲成人久久爱视频| 国产亚洲精品综合一区在线观看| 国产高清激情床上av| 禁无遮挡网站| 国产精品人妻久久久久久| 少妇丰满av| 久久午夜亚洲精品久久| 亚洲天堂国产精品一区在线| 国产精品久久电影中文字幕| 国产精品永久免费网站| 干丝袜人妻中文字幕| 色视频www国产| 亚洲美女搞黄在线观看 | 看片在线看免费视频| 小说图片视频综合网站| 男人舔奶头视频| 长腿黑丝高跟| 日本 欧美在线| 欧美色欧美亚洲另类二区| 欧美一级a爱片免费观看看| 国产精品国产高清国产av| 亚洲成人精品中文字幕电影| 国产国拍精品亚洲av在线观看| 久久久国产成人免费| 成人精品一区二区免费| 亚洲专区国产一区二区| 村上凉子中文字幕在线| 精品国产三级普通话版| 淫秽高清视频在线观看| 久久久久性生活片| 中文字幕免费在线视频6| 成人一区二区视频在线观看| 久久精品国产亚洲av涩爱 | 俄罗斯特黄特色一大片| 丰满乱子伦码专区| av在线天堂中文字幕| 亚洲最大成人手机在线| 亚洲精品乱码久久久v下载方式| 日韩欧美在线乱码| 永久网站在线| 免费av毛片视频| 亚洲aⅴ乱码一区二区在线播放| 日本黄色片子视频| 人妻丰满熟妇av一区二区三区| 看十八女毛片水多多多| 成人三级黄色视频| 久久久久久国产a免费观看| 18禁裸乳无遮挡免费网站照片| 能在线免费观看的黄片| 亚洲人成网站在线播| 久久人人爽人人爽人人片va| 免费不卡的大黄色大毛片视频在线观看 | 免费人成视频x8x8入口观看| 国产亚洲欧美98| 欧美日本视频| 午夜福利在线观看免费完整高清在 | 成年人黄色毛片网站| 直男gayav资源| 国产伦在线观看视频一区| 麻豆久久精品国产亚洲av| 久久婷婷人人爽人人干人人爱| 九色成人免费人妻av| 天美传媒精品一区二区| 男人的好看免费观看在线视频| 免费看av在线观看网站| 少妇的逼水好多| 永久网站在线| 99热只有精品国产| 亚洲欧美日韩高清专用| 国产色爽女视频免费观看| 无遮挡黄片免费观看| 午夜精品在线福利| 国产精品1区2区在线观看.| 噜噜噜噜噜久久久久久91| 天堂影院成人在线观看| 日本色播在线视频| 午夜激情欧美在线| 麻豆久久精品国产亚洲av| 69av精品久久久久久| 成人一区二区视频在线观看| 香蕉av资源在线| 麻豆精品久久久久久蜜桃| 99视频精品全部免费 在线| 免费大片18禁| 成人欧美大片| 我的女老师完整版在线观看| 特大巨黑吊av在线直播| 神马国产精品三级电影在线观看| av女优亚洲男人天堂| 国产精品国产高清国产av| 久久精品综合一区二区三区| 性欧美人与动物交配| 熟妇人妻久久中文字幕3abv| 永久网站在线| 波多野结衣巨乳人妻| 久久精品91蜜桃| 成人三级黄色视频| 国产精品av视频在线免费观看| 欧美xxxx黑人xx丫x性爽| 又黄又爽又刺激的免费视频.| 中文亚洲av片在线观看爽| 天堂√8在线中文| 欧美xxxx性猛交bbbb| 老司机深夜福利视频在线观看| 免费电影在线观看免费观看| 亚洲精品影视一区二区三区av| 国产探花在线观看一区二区| 午夜福利在线观看免费完整高清在 | 亚洲av日韩精品久久久久久密| 韩国av在线不卡| 真人做人爱边吃奶动态| 此物有八面人人有两片| 亚洲一级一片aⅴ在线观看| 一区二区三区四区激情视频 | 搡女人真爽免费视频火全软件 | 久久精品国产自在天天线| 亚洲精品影视一区二区三区av| 亚洲国产精品久久男人天堂| 亚洲性久久影院| 我的老师免费观看完整版| 亚洲熟妇熟女久久| 成人无遮挡网站| 亚洲精品乱码久久久v下载方式| 网址你懂的国产日韩在线| 国产 一区精品| 搡老妇女老女人老熟妇| 国产精品一区www在线观看 | 国产成年人精品一区二区| 一区二区三区免费毛片| 91狼人影院| 久久人人爽人人爽人人片va| 久久久精品大字幕| 亚洲精品色激情综合| 精品午夜福利视频在线观看一区| 亚洲精品一区av在线观看| 国产真实乱freesex| 国产熟女欧美一区二区| 久久久久国内视频| 欧美日韩综合久久久久久 | 日日摸夜夜添夜夜添小说| 99热网站在线观看| 又爽又黄a免费视频| 少妇人妻一区二区三区视频| 97超级碰碰碰精品色视频在线观看| 不卡一级毛片| 美女cb高潮喷水在线观看| 国产毛片a区久久久久| 久久久精品欧美日韩精品| 啦啦啦韩国在线观看视频| 丰满人妻一区二区三区视频av| 免费黄网站久久成人精品| 免费看a级黄色片| 国产精品美女特级片免费视频播放器| 亚洲第一区二区三区不卡| 午夜福利在线观看吧| 日韩精品中文字幕看吧| 欧美激情在线99| 免费av观看视频| 国产精品三级大全| 天堂av国产一区二区熟女人妻| 五月玫瑰六月丁香| 国内精品久久久久精免费| www.色视频.com| 亚洲久久久久久中文字幕| 亚洲精品日韩av片在线观看| 亚洲精品456在线播放app | 日韩欧美在线二视频| 亚洲精品日韩av片在线观看| 99久国产av精品| 午夜影院日韩av| 很黄的视频免费| 大型黄色视频在线免费观看| 午夜激情欧美在线| 婷婷精品国产亚洲av在线| bbb黄色大片| 欧美日韩黄片免| 久久久久国内视频| 在线观看66精品国产| 岛国在线免费视频观看| 国产精品一区二区三区四区免费观看 | 69人妻影院| 欧美日韩国产亚洲二区| 国内毛片毛片毛片毛片毛片| 内射极品少妇av片p| 欧美性感艳星| 亚洲国产日韩欧美精品在线观看| 久久久精品欧美日韩精品| 少妇的逼水好多| 久久久久久国产a免费观看| 精品国内亚洲2022精品成人| 亚洲国产欧美人成| 18+在线观看网站| 女生性感内裤真人,穿戴方法视频| 日韩大尺度精品在线看网址| 精品久久国产蜜桃| 国产激情偷乱视频一区二区| 成年女人永久免费观看视频| 真人做人爱边吃奶动态| 久久久精品欧美日韩精品| 午夜日韩欧美国产| 深爱激情五月婷婷| 欧美人与善性xxx| 黄色视频,在线免费观看| 国产精品不卡视频一区二区| 午夜亚洲福利在线播放| 色综合色国产| 久久天躁狠狠躁夜夜2o2o| .国产精品久久| 色精品久久人妻99蜜桃| 日韩中文字幕欧美一区二区| 国产亚洲av嫩草精品影院| 亚洲电影在线观看av| 99久久精品国产国产毛片| av.在线天堂| 99久久精品国产国产毛片| 动漫黄色视频在线观看| 简卡轻食公司| 久久欧美精品欧美久久欧美| 麻豆精品久久久久久蜜桃| 最好的美女福利视频网| 久久久国产成人精品二区| 2021天堂中文幕一二区在线观| 韩国av在线不卡| av中文乱码字幕在线| 夜夜夜夜夜久久久久| 99国产极品粉嫩在线观看| 99热这里只有是精品在线观看| 国产三级在线视频| 久久草成人影院| 久久精品国产自在天天线| 亚洲男人的天堂狠狠| 婷婷六月久久综合丁香| av在线老鸭窝| 有码 亚洲区| 亚洲中文日韩欧美视频| 制服丝袜大香蕉在线| 免费看光身美女| 丝袜美腿在线中文| 国产精品98久久久久久宅男小说| 中国美女看黄片| .国产精品久久| 18禁裸乳无遮挡免费网站照片| 欧美三级亚洲精品| 麻豆成人午夜福利视频| 精品人妻1区二区| 97超级碰碰碰精品色视频在线观看| 午夜福利18| 性欧美人与动物交配| 欧美日韩精品成人综合77777| 亚洲久久久久久中文字幕| 久久精品国产亚洲av香蕉五月| 久久久久久久久久久丰满 | 精品99又大又爽又粗少妇毛片 | 日韩一区二区视频免费看| 国产一区二区激情短视频| 欧美高清成人免费视频www| 亚洲av一区综合| 一个人看视频在线观看www免费| 精品免费久久久久久久清纯| 色综合色国产| 免费观看的影片在线观看| 国产一区二区三区视频了| 久99久视频精品免费| 一本久久中文字幕| 两个人的视频大全免费| 中文字幕精品亚洲无线码一区| 人人妻,人人澡人人爽秒播| 男人狂女人下面高潮的视频| 午夜a级毛片| 亚洲天堂国产精品一区在线| 国产成人福利小说| 伦精品一区二区三区| 亚洲av第一区精品v没综合| 黄色一级大片看看| 午夜老司机福利剧场| 欧美日韩综合久久久久久 | netflix在线观看网站| 美女xxoo啪啪120秒动态图| 国产免费男女视频| 在线观看午夜福利视频| 国产精品一区二区性色av| 三级毛片av免费| 国产精品福利在线免费观看| 亚洲不卡免费看| 国产91精品成人一区二区三区| 免费在线观看日本一区| 舔av片在线| 亚洲性久久影院| 国产一区二区激情短视频| 国产成人a区在线观看| 最新在线观看一区二区三区| 中文字幕久久专区| 88av欧美| 高清在线国产一区| 久久精品国产清高在天天线| 精品久久久久久成人av| 亚洲精品成人久久久久久| 国产69精品久久久久777片| av在线观看视频网站免费| 夜夜夜夜夜久久久久| 2021天堂中文幕一二区在线观| 一进一出抽搐动态| 免费在线观看影片大全网站| 嫁个100分男人电影在线观看| h日本视频在线播放| 亚洲欧美日韩卡通动漫| 黄色配什么色好看| 我要搜黄色片| 日日摸夜夜添夜夜添av毛片 | 国内精品久久久久久久电影| 亚洲国产高清在线一区二区三| 久久久久久久精品吃奶| 欧美高清成人免费视频www| 亚洲熟妇熟女久久| 国产成人a区在线观看| 免费av毛片视频| 九九爱精品视频在线观看| 3wmmmm亚洲av在线观看| 国产亚洲欧美98| 国产男靠女视频免费网站| 99精品久久久久人妻精品| 亚洲第一区二区三区不卡| www.www免费av| 婷婷精品国产亚洲av| 亚洲国产精品成人综合色| 搞女人的毛片| 欧美日韩乱码在线| 色5月婷婷丁香| 久久久精品欧美日韩精品| 国产av麻豆久久久久久久| 亚洲狠狠婷婷综合久久图片| 国产精品1区2区在线观看.| 麻豆成人午夜福利视频| 亚洲成人中文字幕在线播放| 国产又黄又爽又无遮挡在线| 色综合站精品国产| 国内精品一区二区在线观看| 美女cb高潮喷水在线观看| 人人妻,人人澡人人爽秒播| 欧美高清性xxxxhd video| 久久精品影院6| 内射极品少妇av片p| 国产成人影院久久av| 日韩精品青青久久久久久| 特大巨黑吊av在线直播| 亚洲av成人av| 看免费成人av毛片| 久久九九热精品免费| 国产精品福利在线免费观看| 老女人水多毛片| 嫩草影院精品99| 变态另类成人亚洲欧美熟女| 少妇人妻精品综合一区二区 | 久久久久精品国产欧美久久久| 精品人妻偷拍中文字幕| 午夜a级毛片| 欧美极品一区二区三区四区| 搡老岳熟女国产| 色av中文字幕| 国产主播在线观看一区二区| 精品日产1卡2卡| 久久精品国产自在天天线| 国产v大片淫在线免费观看| 国产精品久久久久久亚洲av鲁大| 毛片女人毛片| 成人av一区二区三区在线看| 欧美xxxx性猛交bbbb| 国产精品福利在线免费观看| 免费看日本二区| www.色视频.com| 国产视频内射| 在线观看av片永久免费下载| 1024手机看黄色片| 成人毛片a级毛片在线播放| 别揉我奶头~嗯~啊~动态视频| 国产老妇女一区| 最近在线观看免费完整版| 亚洲国产高清在线一区二区三| 丝袜美腿在线中文| 人妻制服诱惑在线中文字幕| 国产高清有码在线观看视频| 老司机深夜福利视频在线观看| 嫩草影院新地址| 亚洲国产精品合色在线| 日本 欧美在线| 欧美区成人在线视频| 亚洲最大成人手机在线| 亚洲成人久久性| 99精品久久久久人妻精品| 久久精品人妻少妇| 99热6这里只有精品| 国产精品无大码| 蜜桃久久精品国产亚洲av| 日日撸夜夜添| 国产色婷婷99| 色综合站精品国产| 中文资源天堂在线| 亚洲va在线va天堂va国产| 午夜影院日韩av| 淫妇啪啪啪对白视频| 少妇裸体淫交视频免费看高清| 日本黄大片高清| 欧美日韩亚洲国产一区二区在线观看| 琪琪午夜伦伦电影理论片6080| 看十八女毛片水多多多| 国内精品久久久久久久电影| 国产精品免费一区二区三区在线| 两性午夜刺激爽爽歪歪视频在线观看| 久久久国产成人免费| 久久亚洲精品不卡| 在线播放无遮挡| 给我免费播放毛片高清在线观看| 国产 一区精品| 欧美潮喷喷水| videossex国产| 午夜影院日韩av| 亚洲欧美日韩高清在线视频| а√天堂www在线а√下载| 国内少妇人妻偷人精品xxx网站| 国产精品永久免费网站| 国产精品女同一区二区软件 | 午夜亚洲福利在线播放| 国产女主播在线喷水免费视频网站 | 婷婷丁香在线五月| 成年女人永久免费观看视频| 国产老妇女一区| 国产精品伦人一区二区| 精品久久国产蜜桃| 女同久久另类99精品国产91| 又爽又黄无遮挡网站| 国产一区二区在线av高清观看| 一级黄色大片毛片| 男人狂女人下面高潮的视频| 女生性感内裤真人,穿戴方法视频| 国产高清不卡午夜福利| 成人二区视频| 国国产精品蜜臀av免费| 亚洲专区国产一区二区| 免费看av在线观看网站| 国产一区二区在线av高清观看| 麻豆成人av在线观看| 蜜桃亚洲精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 日韩中文字幕欧美一区二区| 午夜福利成人在线免费观看| 久久精品综合一区二区三区| 日本黄色视频三级网站网址| 亚洲人成网站在线播| av在线天堂中文字幕| 国产免费一级a男人的天堂|