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

    聲學高溫計拓撲結(jié)構(gòu)優(yōu)化研究

    2023-12-07 08:21:34趙乂鋆趙儉陳鑫虎
    計測技術(shù) 2023年5期
    關(guān)鍵詞:聲壓級聲學聲波

    趙乂鋆,趙儉,陳鑫虎

    (航空工業(yè)北京長城計量測試技術(shù)研究所,北京 100095)

    0 引言

    隨著對極端環(huán)境測溫需求的提高,測量諸如高溫風洞試驗段和航空發(fā)動機燃燒室等高溫、高流速的氣流溫度場環(huán)境的需求也日漸攀升[1]。這類測溫環(huán)境中,有著氣流場復(fù)雜、溫度場變化速度快、環(huán)境干擾強等特點;甚至部分環(huán)境會因密封和防輻射等問題,無法將敏感元件直接浸入到待測氣流場中。復(fù)雜多變的氣流溫度場使得采用傳統(tǒng)的熱電偶等接觸式測溫方式時間常數(shù)較大、不耐腐蝕和輻射等問題被放大,不能很好地滿足高溫、高流速的氣流溫度場環(huán)境的測溫需求[2]。而聲學測溫作為一種較新的測溫方法,具有抗腐蝕、抗振動、測溫范圍廣等特點,逐漸受到各國研究人員的重視。本文的主要研究內(nèi)容是非接觸式的聲學測溫方法。

    從原理上來說,非接觸式的聲學測溫方法大致有聲速法測溫、聲共鳴法測溫以及聲頻率法測溫等[3-5]。在氣流溫度測量中,聲速法測溫與聲共鳴法測溫的研究與應(yīng)用更成熟。這兩種方法最基本的原理都是通過測量聲音在氣流場中的傳播時間,進而反算出聲波傳播路徑上氣流溫度場的平均溫度,再根據(jù)解耦和插值等算法得出溫度數(shù)值或是其他相關(guān)參數(shù)。其中,聲共鳴法測溫更加關(guān)注溫標的建設(shè)與準確度;聲速法測溫更加關(guān)注氣流平面溫度場的還原。因此,在對氣流溫度場進行測量與還原時,首要應(yīng)考慮運用聲速法測溫,進行聲學高溫計的研究。

    國內(nèi)的聲學高溫計相關(guān)研究始于20世紀末[6]。時至今日,聲學高溫計在爐膛等溫度較高、氣流速度不快的封閉空間溫度測量上取得了較為成熟的研究成果[7-12]。當溫度升高時,聲探頭振動的感受部分會受到高溫影響會導(dǎo)致性能下降,直接影響到發(fā)射和接收聲信號的完整程度,進而影響到溫度的測量;此外,聲波在高流速的氣流場傳播時,會在矢量合成的作用下明顯改變傳播方向,使得預(yù)計的聲波傳播距離和實際的聲波傳播距離產(chǎn)生差異,增大聲學高溫計的測溫誤差。當氣流溫度進一步升高,速度進一步加快時,氣流溫度場的測溫環(huán)境將會變得更加惡劣。根據(jù)之前的研究[10],目前現(xiàn)有的一套聲學高溫計在600 ℃,0.2 Ma 的情況下,現(xiàn)有聲學高溫計的中心區(qū)誤差為2%。但當溫度和氣流馬赫數(shù)上升時,中心區(qū)的誤差可達20%以上,且一段時間內(nèi),還原的中心溫度差值可達理論數(shù)值的40%,嚴重影響了溫度場還原的精度和準確性。

    針對以上聲學高溫計在高溫、高流速等極端環(huán)境測溫中存在的普遍問題,采用仿真的方式對聲波傳播路徑在極端環(huán)境中的變化情況進行模擬。以現(xiàn)有的一套聲學高溫計裝置為基礎(chǔ)模型,設(shè)置與變化情況匹配的聲探頭拓撲結(jié)構(gòu)并在模擬環(huán)境中對發(fā)射與接收的聲信號參數(shù)進行仿真優(yōu)化,提高測溫穩(wěn)定性和精度。

    1 聲學高溫計基本原理

    現(xiàn)有聲學高溫計裝置的核心測溫原理是非接觸式聲學測溫中的聲速法測溫。溫場中的聲速由一維微元運動公式[13]推導(dǎo)得出。

    式中:c為溫場中的聲速;γ為被測介質(zhì)的絕熱狀態(tài)指數(shù);R、M分別為摩爾氣體常數(shù)和氣體摩爾質(zhì)量;T為熱力學溫度;Z 為常數(shù),若被測介質(zhì)是空氣,一般值取為20.05。

    在已知聲波穿過待測溫場路徑長度的情況下,為了得到溫場中的聲速c,需要獲取聲波在溫場中的傳播時間。聲波傳播時間可由互相關(guān)法計算得到。在直接進行互相關(guān)計算之前,需要對接收到的聲信號進行包絡(luò)提取,目前常用的方法為希爾伯特變換法。希爾伯特變換法即通過希爾伯特變換求得獲取到的聲波信號包絡(luò)[14],即

    式中:?(t)為連續(xù)時間信號x(t)的希爾伯特變換;t為時間;*為卷積運算符號;j為復(fù)數(shù)符號;y(t)為x(t)的信號包絡(luò)。

    獲取包絡(luò)之后,便可以通過計算機軟件準確計算出聲波的傳播時間,進而獲得聲速。當被測介質(zhì)的濃度和成分穩(wěn)定時,氣體摩爾質(zhì)量為常數(shù),Z 為定值。從式(1)中可得,聲速c和溫度T在理想情況下是單值對應(yīng)關(guān)系。以此得到該聲波路徑上的平均溫度。

    當多條聲波路徑上的平均溫度已知時,即可通過矩陣解耦算法計算得出特定區(qū)域的平均溫度。常用的矩陣解耦算法是最小二乘法。最小二乘法的基本原理是將整個待測平面劃分為若干相同的像素,再根據(jù)聲速與溫度的關(guān)系計算出每一個像素的平均溫度,最后用插值的方式還原整個溫度場[15]?,F(xiàn)有聲學高溫計裝置中使用的最小二乘法基本原理如下:

    首先,每一條聲波沿著其傳播路徑傳播時間t都可以表示為

    式中:a為空間狀態(tài)因子,表示穿越氣流溫場的平均聲速c的倒數(shù),ds表示聲波傳播路徑長度的微分值。

    接著,給每一條路徑和每一個像素編號,令ΔSk,i表示第k條路徑通過第i個像素的長度,ai表示第i個像素的聲波平均傳播速度c的倒數(shù),根據(jù)式(3)可知理論上聲波在第k條路徑的傳播時間tk為

    通過實驗測得的聲波傳播時間為tc,則可計算出聲波傳播時間的誤差εk為

    然后,利用極值求解的方法使誤差εk的平方和最小,得到一個含參數(shù)的正則矩陣。當聲探頭的位置和數(shù)量固定時,正則矩陣參數(shù)為常數(shù)。當矩陣可逆時,可得空間狀態(tài)因子矩陣A。又由于空間狀態(tài)因子是聲波傳播速度的倒數(shù),則可以得到每個像素內(nèi)的區(qū)域平均溫度為

    最后,利用適當?shù)牟逯邓惴?,就可以重建出整個待測平面溫度場的分布。

    2 仿真參數(shù)設(shè)置

    2.1 仿真分析背景

    在利用聲波傳播時間計算氣流溫度時,變化的氣流速度會導(dǎo)致聲信號的聲壓級降低。一般來說,聲壓級越高,聲信號的能量就越強,聲壓的幅值也越大,受到環(huán)境干擾的影響就越小。加上環(huán)境干擾的影響,聲探頭獲取的聲源信號會產(chǎn)生變形。變形聲信號會影響聲波在氣流場中傳播時間獲取的穩(wěn)定性和精度,進一步使氣流溫度場還原的精度以及可靠性降低。根據(jù)矢量疊加的原理可知,聲波在氣流中傳播時,實際的傳播方向為聲波矢量與氣流矢量的疊加方向?,F(xiàn)有聲學高溫計裝置的聲探頭在同一水平面上,傾斜的聲波傳播方向勢必會導(dǎo)致接收到的聲信號聲壓級減小,更易受到環(huán)境干擾。因此,需要從硬件角度尋找并設(shè)計一種適用于當前研究指標的聲探頭拓撲框架,提高聲學高溫計測溫的精度和可靠程度。

    本文設(shè)置不同分布半徑和不同朝向夾角的聲探頭支架作為仿真分析的基礎(chǔ)模型。以Virtual.Lab為仿真軟件,通過設(shè)置恒溫度不同馬赫數(shù)以及恒馬赫數(shù)不同溫度的兩類條件分別對不同結(jié)構(gòu)的聲探頭支架進行聲壓級降低的仿真分析。通過仿真的方式尋找出最佳的聲探頭支架拓撲結(jié)構(gòu)。

    2.2 聲探頭結(jié)構(gòu)的仿真建模及環(huán)境設(shè)置

    為了分析聲探頭朝向?qū)邮章曅盘柕挠绊?,在仿真中分別設(shè)置了以15 cm 和20 cm 為直徑的兩種聲探頭支架。每一種直徑的支架又分別設(shè)計了發(fā)射聲探頭與接收聲探頭朝向正對、探頭朝向與支架平面夾角15°、30°以及45°的情況,總計8 種聲探頭的支架。以直徑和探頭夾角的方式命名,中間用短橫線進行連接。如“15-dir”表示15 cm直徑探頭朝向正對的聲探頭支架,“15-15”表示15 cm 直徑探頭朝向與支架平面夾角15°的聲探頭支架。以15 cm 直徑的聲探頭支架為例,在Virtual.Lab 中建模并簡化其他非必要構(gòu)件后,部分15 cm直徑聲探頭支架Virtual.Lab模型如圖1所示。

    圖1 部分15 cm直徑聲探頭支架Virtual.Lab模型Fig.1 Partial 15 cm diameter acoustic probe support model in Virtual.Lab

    仿真采用有限元仿真的方式進行。將聲探頭模型放入到圓柱形分析環(huán)境中模擬氣流通過部分并網(wǎng)格劃分,得到需要進行分析的有限元分析三維網(wǎng)格劃分,如圖2所示。

    圖2 有限元分析三維網(wǎng)格劃分剖面圖Fig.2 3D mesh division profile of finite element analysis

    最后對建立的網(wǎng)格進行聲環(huán)境以及邊界條件設(shè)置,增加監(jiān)控聲壓變化的場點,如圖3所示。至此完成不同聲探頭結(jié)構(gòu)的仿真建模及環(huán)境設(shè)置。

    圖3 增加場點及邊界條件的有限元分析模型Fig.3 Finite element analysis model with field points and boundary conditions

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

    3.1 恒溫度不同氣流速度下各聲探頭支架仿真情況

    為了分析各聲探頭支架在不同氣流速度下的各不同方向上聲壓級降低情況,設(shè)置恒溫度不同氣流速度的仿真條件。其中,設(shè)置環(huán)境溫度為常溫,氣流速度以0.1 Ma為間隔,最大值為0.3 Ma。由于同一種探頭支架內(nèi)還存在多個聲波接收點。故在分析時分別將正對方向、偏斜22.5°、偏斜45°的聲壓級變化情況命名為“正對方向變化情況”、“側(cè)對方向變化情況”、“邊緣方向變化情況”,以此為分類計算不同聲波接收探頭處的聲壓級占比聲波發(fā)射處聲壓級的比值。其中,兩種聲探頭支架在0.2 Ma氣流速度下的聲壓級分布示意圖如圖4所示。

    圖4 部分聲探頭支架0.2 Ma流速下聲壓級分布示意圖Fig.4 Schematic diagram of sound pressure level distribution at a flow rate of 0.2 Ma for some sound probe supports

    以聲探頭發(fā)射點的聲壓級為標準,仿真計算并記錄不同聲探頭支架不同方向上的聲壓級數(shù)值,令傳播后的聲壓級數(shù)值與發(fā)射點的聲壓級數(shù)值相比作為聲壓級衰減的評判標準,將各個聲探頭支架在恒溫不同流速情況下的聲壓級的衰減情況進行匯總,得到的數(shù)據(jù)如表1所示。

    表1 各聲探頭支架在不同氣流馬赫數(shù)下接收聲壓級占聲源聲壓級比值Tab.1 Ratio of received sound pressure level to source sound pressure level of each acoustic probe support at different airflow Mach numbers

    表1中記錄在不同聲探頭支架內(nèi)徑中,不同支架類型常溫環(huán)境下分別在0.1、0.2、0.3 Ma 時不同聲探頭偏斜角度的聲壓級減小情況。由表1分析可知,當氣流速度逐漸增大時,正對朝向的聲探頭支架的聲壓級分布的均勻性和對稱性越差,且聲壓級減小的程度越大。說明在氣流速度增大時,正對朝向的聲壓級會受到氣流速度變化帶來的較為明顯的影響。

    此外,當氣流速度逐漸增大時,傾斜角度的聲探頭支架的聲壓級并未像正對的聲探頭支架一樣逐漸下降,而是存在上升后又下降的趨勢。對于直徑15 cm 的聲探頭支架組而言,除了15-45 聲探頭支架之外,總體上滿足正對方向上的聲壓級分布情況較側(cè)對方向以及邊緣方向較好;側(cè)對方向上的聲壓分布較邊緣方向較好的特征。三個方向上,聲壓級變化基本符合先增大后減小的趨勢。

    經(jīng)由先前的研究[10]可知,聲探頭的聲波接收范圍存在角度限制。當探頭的距離相同時,聲探頭發(fā)射與接收口之間對得越正,聲壓級的減小情況就越小。反映在15 cm 直徑上的聲探頭支架上時,正對方向上的聲探頭之間相對距離最遠,但因為朝向正對,使得聲探頭的聲源信號有較大區(qū)域的接收口。因此正對方向上的聲壓級減小得少。與之相對的,邊緣方向上的聲壓級減小要多得多。大體上與仿真結(jié)果相符。

    以15 cm 直徑聲探頭支架為例,若以矢量合成的計算方法可知,當聲速恒定的情況下,氣流在0.1~0.3 Ma 時分別傳播了1.5、3.0、4.5 cm。將氣流傳播的距離與探頭朝向的距離進行對比可知,在0.1 Ma 時,正對、側(cè)對和邊緣方向上的氣流傳播距離與15-15 和15-30 聲探頭支架的偏斜角度較為靠近,體現(xiàn)為兩聲探頭支架三個方向上聲壓級衰減較小,但總體上相差并不大;當氣流馬赫數(shù)升高時,氣流傳播的矢量距離變長,有一定傾斜角度的聲探頭支架與合成后的聲波實際傳播方向更為接近。具體表現(xiàn)為在氣流馬赫數(shù)上升,環(huán)境干擾更強,整體聲壓級衰減增大的情況下,15-15和15-30 聲探頭支架的整體聲壓級減小相比其他聲探頭支架要小得多,這與仿真結(jié)果是相符合的。值得注意的是,在邊緣方向上,由于聲波傳播的距離過短,其聲矢量和氣流矢量合成后的角度大致為40°,與仿真結(jié)果中在15-30 和15-45 聲探頭支架邊緣方向上的聲壓級衰減情況相符。至于15-45聲探頭支架在側(cè)對方向上相較于正對和邊緣聲壓級均較小的情況,推測可能是由于傾斜角度過大,導(dǎo)致仿真并未準確解算出實際矢量合成后的聲壓級減小情況。

    總體來說,在直徑15 cm 的聲探頭支架組內(nèi)部,15-30 的聲探頭支架在各個流速情況下均有較好的聲壓級分布和較小的聲壓級減小,綜合各個流場下的聲壓級分布較優(yōu)。

    對于直徑20 cm 的聲探頭支架組而言,正對方向上的聲壓級占比情況較側(cè)對方向略高,較邊緣方向更高。呈現(xiàn)出較為明顯的3 層分布趨勢。同時,3 個方向上的聲壓級隨著氣流溫度改變而波動降低。

    以相同的方法對20 cm 的聲探頭支架組分析可知,當聲探頭支架的直徑增大時,氣流在不同馬赫數(shù)下傳播矢量長度不變。因此在聲探頭直徑增大時,聲矢量與氣流矢量合成的角度相較于15 cm聲探頭支架而言會變小,體現(xiàn)在聲壓級變化上則表現(xiàn)為20-15 與20-30 聲探頭支架的聲壓級之差比15-15 與15-30 聲探頭支架的聲壓級之差小,符合仿真結(jié)果的規(guī)律。

    綜合對比直徑20 cm 的聲探頭支架組的聲壓級減小情況,20-30 的聲探頭支架綜合各個流場下的聲壓級分布較優(yōu)。

    綜合比較15-30 以及20-30 兩個聲探頭支架的聲壓級分布情況。15-30 的聲探頭支架有著更好的聲壓級分布以及較少的聲壓級減小。當氣流速度較快時,20-30 聲探頭支架的聲壓級受氣流速度影響產(chǎn)生的變化相較于15-30 聲探頭支架更為明顯。經(jīng)過之前的計算分析,推測在相同的聲探頭傾斜角度下時,較長的傳播路徑會增大對聲壓級較小的影響。此外,20-30 聲探頭支架側(cè)對方向以及邊緣方向的聲壓級減小情況明顯比15-30 聲探頭支架大,這也有可能是受到相同因素的影響。故綜合比較可知,當恒定氣流溫度改變馬赫數(shù)時,15-30聲探頭支架的結(jié)構(gòu)更優(yōu)。

    3.2 恒馬赫數(shù)不同氣流溫度下各聲探頭支架仿真情況

    為了分析各聲探頭支架在不同氣流溫度下的各不同方向上聲壓級降低情況,按照同樣的方式,設(shè)置了恒馬赫數(shù)不同氣流溫度的仿真條件。其中,設(shè)置馬赫數(shù)恒為0.3 Ma,氣流溫度分別為500、700、900 ℃,對各聲探頭支架的聲壓級降低情況進行分析,部分結(jié)果如圖5所示。將各個聲探頭支架在恒馬赫數(shù)不同溫度情況下的聲壓級的衰減情況進行匯總,得到的結(jié)果如表2所示。

    表2 各聲探頭支架在不同氣流溫度下接收聲壓級占聲源聲壓級比值Tab.2 Ratio of received sound pressure level to source sound pressure level of each acoustic probe support at different temperatures

    表2中記錄在在不同聲探頭支架內(nèi)徑中,不同支架類型分別在恒流速時,500、700、900 ℃下不同聲探頭偏斜角度的聲壓級減小情況。

    當氣流馬赫數(shù)恒定時,聲速隨著溫度的升高而變大,而氣流的實際傳播速度也隨著溫度的增大而增大。實際上,當氣流馬赫數(shù)恒定時,聲矢量與氣流矢量的合成矢量方向是確定的。對15 cm的聲探頭支架組而言,經(jīng)過矢量合成的計算可知,正對方向、側(cè)對方向和邊緣方向上的總體傾斜角度在10° ~ 30°之間。也即是從理論計算上來說,聲壓級分布情況較優(yōu)的支架應(yīng)為15°或30°角傾斜的聲探頭支架。

    對表2 中的數(shù)據(jù)進行歸納可知,對于直徑15 cm 的聲探頭支架組而言,總體上滿足正對方向上的聲壓級分布情況較側(cè)對方向以及邊緣方向較好;側(cè)對方向上的聲壓分布較邊緣方向較好的特征。正對方向以及側(cè)對方向聲壓級變化基本符合先增大后減小的趨勢,而邊緣方向的聲壓級符合先減小后增大的波動趨勢。15-15 和15-30 兩聲探頭支架確實具有相對較小的聲壓級減小情況,與理論分析的結(jié)果大致符合??傮w來說,在直徑15 cm 的聲探頭支架組內(nèi)部,15-30 的聲探頭支架綜合各個流場下的聲壓級均在0.820及以上,對于不同溫度情況下的氣流場適應(yīng)性較好。

    對于直徑20 cm 的聲探頭支架組而言,正對方向上的聲壓級分布情況與側(cè)對方向相當,較邊緣方向較好。正對方向并無明顯的隨著氣流溫度改變而變化的規(guī)律,而邊緣方向的聲壓級符合總體上升的趨勢。

    同理,根據(jù)上面20 cm 聲探頭支架的分析以及3.1 節(jié)中的總結(jié)??梢杂嬎愠鍪噶亢铣傻慕嵌仍谡龑Ψ较颉?cè)對方向和邊緣方向上的總體傾斜角度6°~20°。與3.1 節(jié)中20 cm 聲探頭支架在傾斜角度的分析上相符。也即從理論分析上來說,20 cm 聲探頭支架組中,20-15 聲探頭支架在理論上會獲得更好的聲壓級分布情況。通過對比表2 中20 cm 聲探頭支架組在各方向上的聲壓級減小情況可知,理論分析與仿真的結(jié)果基本相符。20-15 的聲探頭支架綜合各個流場下的聲壓級在0.83 及以上,其聲壓級相較于同組其他聲探頭支架而言分布較優(yōu)。

    綜合比較15-30 以及20-15 兩個聲探頭支架的聲壓級分布情況。15-30 的聲探頭支架有著更好的聲壓級分布以及較少的聲壓級減小,但其受溫度影響產(chǎn)生的變化相較于20-15 聲探頭支架更為明顯。當氣流溫度較高時,20-15 聲探頭支架的聲壓級減小情況明顯比15-30 聲探頭支架小。即當恒定馬赫數(shù)改變氣流溫度時,20-15聲探頭支架的結(jié)構(gòu)更優(yōu)。

    由3.1 與3.2 中的分析可知,15-30 和20-15 分別在恒溫變馬赫數(shù)及恒馬赫數(shù)變溫的測溫環(huán)境中表現(xiàn)最好。對比兩種聲探頭拓撲結(jié)構(gòu),雖然20-15 聲探頭支架在恒定馬赫數(shù)情況下時要比15-30 聲探頭支架更好,但15-30 聲探頭支架在恒定溫度時得到聲壓級減小以及分布情況要比20-15 聲探頭支架好得多。因此,從仿真的聲壓級減小結(jié)果而言,15-30 聲探頭支架相對來說是最優(yōu)解。若從給聲探頭支架增加冷卻管或防止聲探頭的端口處因高溫融化等加工或?qū)嶋H試驗的可靠性角度考慮,也可將20-15聲探頭支架作為聲探頭拓撲結(jié)構(gòu)的備用選項。

    4 試驗測試

    為了實際測試不同聲探頭拓撲結(jié)構(gòu)對測溫精度和穩(wěn)定性的影響,將試驗環(huán)境設(shè)置為300 ℃,0.2 Ma 的氣流,設(shè)計了15-dir、15-30、20-15 三個聲探頭支架作為對照試驗組進行對比。由于試驗環(huán)境中的干擾過多,對聲信號的聲壓級測量會受到較大的氣流環(huán)境噪聲干擾,因此以1 min 內(nèi)三個方向上測得的聲波飛渡時間數(shù)據(jù)為依據(jù),計算其樣本標準差、平均值與理論數(shù)值的差值和理論數(shù)值的比值(相對誤差比值),分別考察不同聲探頭拓撲結(jié)構(gòu)在三個方向上的測溫穩(wěn)定性和精度。樣本標準差越低,穩(wěn)定性越好;比值越低,精度越高。得到的數(shù)據(jù)結(jié)果如表3及表4所示。

    表3 三組聲探頭結(jié)構(gòu)采集數(shù)據(jù)標準差Tab.3 Standard deviation of data collected from three acoustic probe structures

    表4 三組聲探頭結(jié)構(gòu)采集數(shù)據(jù)相對誤差比值Tab.4 Relative error ratio of data collected from three acoustic probe structures

    根據(jù)第2節(jié)中的分析可知,聲壓級的變化反映了當前聲信號傳播能力在氣流溫度場的變化情況。當聲信號幅值受環(huán)境影響降低時,接收聲信號波形越光滑,包絡(luò)提取和互相關(guān)計算也越準確,因此聲壓級的減小越少,采集數(shù)據(jù)的標準差和相對誤差比值也越小。故通過仿真獲取的聲壓級變化趨勢在理論上應(yīng)當與試驗時采集數(shù)據(jù)的標準差和相對誤差比值的變化趨勢相同。從表3 及表4 中可以看出,反映數(shù)據(jù)波動情況的標準差和數(shù)據(jù)精度的相對誤差比值的變化趨勢大致相當。在正對方向上,采集獲得數(shù)據(jù)的精度和穩(wěn)定性最好。側(cè)對方向上要稍差一些,邊緣方向上的情況最差。

    在三個對照試驗組中,正對方向上數(shù)據(jù)精度和穩(wěn)定性最好的是20-15 聲探頭拓撲結(jié)構(gòu),其次是15-30 結(jié)構(gòu),該特征與仿真計算的結(jié)果一致。根據(jù)第3 節(jié)的分析,應(yīng)當是20-15 聲探頭拓撲結(jié)構(gòu)的角度更加靠近實際矢量合成角度的緣故;側(cè)對和邊緣方向上,15-30 結(jié)構(gòu)比20-15 結(jié)構(gòu)有更小的標準差和相對誤差比值,二者比現(xiàn)有15-dir結(jié)構(gòu)均有明顯提升。經(jīng)分析可知,當聲信號傳播的水平距離變短時,氣流矢量的傳播距離不變,因此矢量合成后的角度變大,更加靠近15-30 的聲探頭拓撲結(jié)構(gòu)。此外,20-15 聲探頭拓撲結(jié)構(gòu)有更長的內(nèi)徑和更多的聲信號傳播距離,聲信號在氣流場中的衰減更嚴重,聲信號傳播的穩(wěn)定性也會受到影響而降低。

    以標準差的數(shù)值作為評判標準,15-30 在三個方向上的標準差數(shù)值總和降低22.4%,20-15 在三個方向上的總和降低17.5%,可知15-30 聲探頭拓撲結(jié)構(gòu)的綜合性能更好。

    5 結(jié)論

    從現(xiàn)有的聲學高溫計出發(fā),分析了聲學測溫裝置的測溫原理,及現(xiàn)有聲學高溫計存在的測溫精度和穩(wěn)定性問題的原因。提出并構(gòu)建了將聲學高溫計中的聲探頭朝向進行相應(yīng)的調(diào)整,使聲波和氣流矢量合成之后的實際聲波傳播方向落在接收聲探頭敏感區(qū)域上的聲探頭拓撲結(jié)構(gòu)模型。

    以LMS Virtual.Lab 聲學仿真軟件為仿真環(huán)境,構(gòu)建了數(shù)個不同的分布直徑、不同聲探頭朝向的拓撲結(jié)構(gòu)。設(shè)置好邊界條件和流體參數(shù)后,對各聲探頭拓撲結(jié)構(gòu)進行恒馬赫數(shù)變溫度以及恒溫度變馬赫數(shù)仿真計算。總結(jié)不同聲探頭拓撲結(jié)構(gòu)各方向上的聲壓級減小情況,通過比較分析出15-30和20-15 聲探頭拓撲結(jié)構(gòu)聲壓級綜合減小均在0.2以下,接收到的聲信號強度較高。

    根據(jù)仿真計算的結(jié)果,設(shè)計了同工況下現(xiàn)有聲學高溫計支架與15-30、20-15兩種聲探頭拓撲結(jié)構(gòu)支架的試驗。根據(jù)試驗結(jié)果的數(shù)據(jù)分析可知,在常溫至900 ℃、0~0.3 Ma的氣流溫度場環(huán)境中,15-30 聲探頭拓撲結(jié)構(gòu)的采集數(shù)據(jù)標準差最小,比現(xiàn)有結(jié)構(gòu)降低了22.4%。明顯提升了測溫穩(wěn)定性和精度。

    目前的研究僅對氣流溫度不大于900 ℃,氣流速度不高于0.3 Ma 的環(huán)境進行了分析和試驗,在之后的研究中還需要考慮到隨著溫度上升以及流速加快帶來的壓力及聲波傳播方向變化,進一步提升聲學高溫計的測溫性能。

    猜你喜歡
    聲壓級聲學聲波
    機器噪聲平均聲壓級計算方法差異性實證研究
    電動工具(2024年1期)2024-02-29 01:40:24
    愛的就是這股Hi-Fi味 Davis Acoustics(戴維斯聲學)Balthus 70
    一種計算消聲室聲壓級的新方法
    計量學報(2020年2期)2020-04-11 04:33:22
    全新DXR mkll有源揚聲器
    演藝科技(2019年4期)2019-03-30 03:21:46
    Acoustical Treatment Primer:Diffusion談?wù)劼晫W處理中的“擴散”
    Acoustical Treatment Primer:Absorption談?wù)劼晫W處理中的“吸聲”(二)
    Acoustical Treatment Primer:Absorption 談?wù)劼晫W處理中的“吸聲”
    愛的聲波 將愛留在她身邊
    中國寶玉石(2018年3期)2018-07-09 03:13:58
    聲波殺手
    自適應(yīng)BPSK在井下鉆柱聲波傳輸中的應(yīng)用
    全区人妻精品视频| 亚洲在久久综合| 少妇熟女欧美另类| 精品亚洲乱码少妇综合久久| 三级国产精品片| 少妇人妻久久综合中文| 纯流量卡能插随身wifi吗| 人妻人人澡人人爽人人| 午夜免费鲁丝| 亚洲欧洲精品一区二区精品久久久 | 国产又色又爽无遮挡免| 在线看a的网站| 欧美日韩视频精品一区| 最近中文字幕高清免费大全6| 亚洲欧美中文字幕日韩二区| 亚洲国产欧美日韩在线播放| 夜夜骑夜夜射夜夜干| 亚洲第一区二区三区不卡| 18在线观看网站| 女的被弄到高潮叫床怎么办| 久久精品久久精品一区二区三区| av在线播放精品| 热99国产精品久久久久久7| 啦啦啦视频在线资源免费观看| 精品第一国产精品| 王馨瑶露胸无遮挡在线观看| 9热在线视频观看99| 黑人巨大精品欧美一区二区蜜桃 | videosex国产| 久久久久久久久久成人| 蜜桃在线观看..| 国产成人一区二区在线| 九九在线视频观看精品| 如日韩欧美国产精品一区二区三区| 日韩不卡一区二区三区视频在线| 国国产精品蜜臀av免费| 9191精品国产免费久久| 婷婷成人精品国产| 欧美97在线视频| tube8黄色片| 黄色 视频免费看| 国产乱来视频区| videos熟女内射| 精品人妻熟女毛片av久久网站| 亚洲成人av在线免费| 母亲3免费完整高清在线观看 | 国产日韩欧美亚洲二区| 亚洲欧美成人精品一区二区| 美女脱内裤让男人舔精品视频| 中文字幕另类日韩欧美亚洲嫩草| 国产欧美亚洲国产| 国产有黄有色有爽视频| 国产极品天堂在线| 视频区图区小说| xxx大片免费视频| 91在线精品国自产拍蜜月| 卡戴珊不雅视频在线播放| 久久久久国产精品人妻一区二区| 国产色爽女视频免费观看| 国产在线免费精品| 久久99精品国语久久久| 久久久国产一区二区| 国产亚洲欧美精品永久| 国产亚洲午夜精品一区二区久久| 青青草视频在线视频观看| 久久精品久久精品一区二区三区| 亚洲欧美日韩卡通动漫| 九色亚洲精品在线播放| 欧美精品高潮呻吟av久久| 久久久久久久久久成人| 亚洲人成77777在线视频| 日本与韩国留学比较| 一级片'在线观看视频| 国产av国产精品国产| 极品人妻少妇av视频| 国产视频首页在线观看| 毛片一级片免费看久久久久| 好男人视频免费观看在线| 最近2019中文字幕mv第一页| 九九爱精品视频在线观看| av视频免费观看在线观看| 香蕉丝袜av| 母亲3免费完整高清在线观看 | 永久免费av网站大全| 久久久久久伊人网av| 日韩三级伦理在线观看| 午夜老司机福利剧场| 国产老妇伦熟女老妇高清| 久久久久久久国产电影| 精品人妻熟女毛片av久久网站| 午夜日本视频在线| 少妇的逼好多水| 国产精品嫩草影院av在线观看| 精品亚洲成a人片在线观看| 国产成人精品在线电影| 亚洲一码二码三码区别大吗| 老司机影院毛片| 男女下面插进去视频免费观看 | 久久影院123| 亚洲欧美日韩卡通动漫| 久久99热6这里只有精品| 午夜免费观看性视频| 全区人妻精品视频| www.色视频.com| 啦啦啦在线观看免费高清www| 人人澡人人妻人| 成年人午夜在线观看视频| 18禁观看日本| 亚洲欧美成人精品一区二区| 在线天堂中文资源库| 香蕉丝袜av| 中文字幕另类日韩欧美亚洲嫩草| 日韩中字成人| 国产男人的电影天堂91| 日本爱情动作片www.在线观看| 丰满少妇做爰视频| 激情视频va一区二区三区| 最近的中文字幕免费完整| 亚洲成人av在线免费| 最近最新中文字幕大全免费视频 | 亚洲av成人精品一二三区| 国精品久久久久久国模美| 午夜精品国产一区二区电影| 午夜福利在线观看免费完整高清在| 日本vs欧美在线观看视频| av电影中文网址| 丝袜喷水一区| 国产免费一级a男人的天堂| 夫妻午夜视频| 蜜桃国产av成人99| 欧美日韩一区二区视频在线观看视频在线| 狂野欧美激情性xxxx在线观看| 亚洲精品日本国产第一区| 免费少妇av软件| 男女免费视频国产| 国产精品免费大片| 久久99热6这里只有精品| 不卡视频在线观看欧美| 精品一区二区免费观看| 九色亚洲精品在线播放| 精品酒店卫生间| 国产成人91sexporn| 一本久久精品| 欧美激情 高清一区二区三区| 一区在线观看完整版| 91精品伊人久久大香线蕉| 91成人精品电影| 久久99一区二区三区| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区国产| 久久久久久久久久成人| 成人国产av品久久久| av一本久久久久| 啦啦啦啦在线视频资源| 在线天堂最新版资源| 自拍欧美九色日韩亚洲蝌蚪91| 热99国产精品久久久久久7| 青春草视频在线免费观看| 婷婷色麻豆天堂久久| 国产综合精华液| 欧美亚洲日本最大视频资源| 波野结衣二区三区在线| 久久久精品94久久精品| 老熟女久久久| 视频在线观看一区二区三区| 中文字幕av电影在线播放| 精品99又大又爽又粗少妇毛片| 女性生殖器流出的白浆| 男女高潮啪啪啪动态图| 国产在线一区二区三区精| 纵有疾风起免费观看全集完整版| 搡老乐熟女国产| 香蕉精品网在线| 少妇高潮的动态图| 国国产精品蜜臀av免费| 亚洲国产精品国产精品| 日韩av在线免费看完整版不卡| 久久ye,这里只有精品| 黄色 视频免费看| 久久女婷五月综合色啪小说| 成人18禁高潮啪啪吃奶动态图| 国产69精品久久久久777片| 国产免费现黄频在线看| 少妇人妻久久综合中文| 色哟哟·www| 中文精品一卡2卡3卡4更新| 中文天堂在线官网| 日本与韩国留学比较| 免费人妻精品一区二区三区视频| 人成视频在线观看免费观看| 男人操女人黄网站| xxx大片免费视频| 国产精品久久久久久av不卡| 99热这里只有是精品在线观看| 国产xxxxx性猛交| 香蕉丝袜av| 国产一区亚洲一区在线观看| 啦啦啦视频在线资源免费观看| 女人被躁到高潮嗷嗷叫费观| 国产永久视频网站| 亚洲 欧美一区二区三区| tube8黄色片| 亚洲激情五月婷婷啪啪| 亚洲综合精品二区| 亚洲,欧美精品.| 亚洲伊人久久精品综合| 久久99热6这里只有精品| 99国产综合亚洲精品| 激情五月婷婷亚洲| 十八禁高潮呻吟视频| 日韩一本色道免费dvd| av在线播放精品| 亚洲欧美色中文字幕在线| 天天躁夜夜躁狠狠久久av| 观看美女的网站| 十八禁高潮呻吟视频| 亚洲精品日韩在线中文字幕| 日日摸夜夜添夜夜爱| 久久这里有精品视频免费| av福利片在线| 在线 av 中文字幕| 九草在线视频观看| 美女xxoo啪啪120秒动态图| 国产日韩欧美亚洲二区| 国产 一区精品| 日韩免费高清中文字幕av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲精品久久久久久婷婷小说| 精品国产国语对白av| 蜜桃在线观看..| 日日撸夜夜添| 18禁动态无遮挡网站| 国产免费一级a男人的天堂| 日日爽夜夜爽网站| 国产伦理片在线播放av一区| 国产精品国产三级国产av玫瑰| 啦啦啦啦在线视频资源| 制服丝袜香蕉在线| 日日摸夜夜添夜夜爱| 精品视频人人做人人爽| 十分钟在线观看高清视频www| 性色av一级| 精品熟女少妇av免费看| 丝袜美足系列| 九草在线视频观看| 综合色丁香网| 国产精品国产av在线观看| 精品一区二区免费观看| 99久久精品国产国产毛片| 一本—道久久a久久精品蜜桃钙片| 亚洲av电影在线观看一区二区三区| 97精品久久久久久久久久精品| 免费在线观看完整版高清| 国产 一区精品| 亚洲欧美日韩卡通动漫| 久久久久久人人人人人| 黄色配什么色好看| 如日韩欧美国产精品一区二区三区| 国产高清国产精品国产三级| 免费黄色在线免费观看| 日韩在线高清观看一区二区三区| 免费人妻精品一区二区三区视频| 91精品伊人久久大香线蕉| 天天影视国产精品| 人成视频在线观看免费观看| 久久人人爽人人爽人人片va| 亚洲人与动物交配视频| 天天影视国产精品| 精品人妻偷拍中文字幕| 精品午夜福利在线看| 日韩欧美精品免费久久| 国产精品不卡视频一区二区| 成年美女黄网站色视频大全免费| 久久99精品国语久久久| 欧美性感艳星| 欧美成人午夜免费资源| 日本欧美视频一区| 美女中出高潮动态图| 看十八女毛片水多多多| 99精国产麻豆久久婷婷| 国产欧美另类精品又又久久亚洲欧美| 女性被躁到高潮视频| 久久婷婷青草| 国产精品成人在线| 黄色怎么调成土黄色| 丝袜在线中文字幕| 亚洲成色77777| 一本色道久久久久久精品综合| 色婷婷久久久亚洲欧美| 久久毛片免费看一区二区三区| 国产永久视频网站| 黑人猛操日本美女一级片| 啦啦啦啦在线视频资源| 精品国产一区二区三区久久久樱花| 九色成人免费人妻av| 亚洲av综合色区一区| 十八禁网站网址无遮挡| 国产黄频视频在线观看| 久久青草综合色| 精品一区在线观看国产| 最近手机中文字幕大全| a级片在线免费高清观看视频| 成人黄色视频免费在线看| 国语对白做爰xxxⅹ性视频网站| 久久精品国产亚洲av天美| 黄网站色视频无遮挡免费观看| 精品亚洲成a人片在线观看| 亚洲一码二码三码区别大吗| 国产男女内射视频| 97精品久久久久久久久久精品| 90打野战视频偷拍视频| 一边亲一边摸免费视频| 亚洲欧美中文字幕日韩二区| 亚洲国产毛片av蜜桃av| 寂寞人妻少妇视频99o| 精品午夜福利在线看| 亚洲欧美日韩另类电影网站| 免费在线观看完整版高清| 日韩中字成人| 在现免费观看毛片| 男人添女人高潮全过程视频| 男女边摸边吃奶| 日韩成人av中文字幕在线观看| 十八禁高潮呻吟视频| 国产永久视频网站| 黄色怎么调成土黄色| 国产日韩欧美在线精品| 亚洲精品视频女| 卡戴珊不雅视频在线播放| 免费看光身美女| 中国国产av一级| 久久精品aⅴ一区二区三区四区 | av免费在线看不卡| 丝袜喷水一区| 五月玫瑰六月丁香| 久久国产精品男人的天堂亚洲 | 少妇被粗大的猛进出69影院 | 十八禁高潮呻吟视频| 日韩不卡一区二区三区视频在线| 如何舔出高潮| 日韩精品有码人妻一区| av线在线观看网站| 极品人妻少妇av视频| 日韩制服丝袜自拍偷拍| 美女国产视频在线观看| 国产色爽女视频免费观看| 色婷婷av一区二区三区视频| 国产亚洲av片在线观看秒播厂| 少妇人妻 视频| 成年av动漫网址| 久久精品国产自在天天线| 蜜桃国产av成人99| 国产一区二区在线观看日韩| 纯流量卡能插随身wifi吗| 欧美日韩精品成人综合77777| 亚洲精品中文字幕在线视频| 观看av在线不卡| 婷婷色av中文字幕| 欧美人与性动交α欧美精品济南到 | 黑人欧美特级aaaaaa片| 亚洲国产精品999| 成人漫画全彩无遮挡| 国产成人午夜福利电影在线观看| 2018国产大陆天天弄谢| 免费黄网站久久成人精品| 男女下面插进去视频免费观看 | 在线天堂中文资源库| 天堂8中文在线网| 日韩成人伦理影院| 青青草视频在线视频观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 人人妻人人添人人爽欧美一区卜| 日本黄色日本黄色录像| av卡一久久| 国产成人精品在线电影| 一本久久精品| 成人毛片60女人毛片免费| 在现免费观看毛片| 久久久亚洲精品成人影院| 侵犯人妻中文字幕一二三四区| 一级片免费观看大全| 下体分泌物呈黄色| 亚洲精品乱码久久久久久按摩| 在线免费观看不下载黄p国产| 国产永久视频网站| 黄网站色视频无遮挡免费观看| 大香蕉久久成人网| 久热久热在线精品观看| 免费播放大片免费观看视频在线观看| 亚洲少妇的诱惑av| 久久久久视频综合| 久久久精品免费免费高清| 国产精品人妻久久久影院| 久久国产精品男人的天堂亚洲 | 十八禁网站网址无遮挡| av卡一久久| 大陆偷拍与自拍| 丝袜人妻中文字幕| 少妇熟女欧美另类| 日本色播在线视频| 99久久人妻综合| 亚洲av男天堂| 国产又爽黄色视频| 丰满乱子伦码专区| 十八禁网站网址无遮挡| 男男h啪啪无遮挡| 国产成人精品在线电影| 亚洲婷婷狠狠爱综合网| 人体艺术视频欧美日本| www日本在线高清视频| 边亲边吃奶的免费视频| 韩国av在线不卡| 国产精品一国产av| 99久久人妻综合| 我的女老师完整版在线观看| 久久精品熟女亚洲av麻豆精品| 丰满乱子伦码专区| 国产在线视频一区二区| 日本av手机在线免费观看| 成年人免费黄色播放视频| 91精品伊人久久大香线蕉| 久久久久久久久久久免费av| 欧美 亚洲 国产 日韩一| 久久99一区二区三区| 亚洲av免费高清在线观看| 女人久久www免费人成看片| 在线观看美女被高潮喷水网站| 九草在线视频观看| 夜夜骑夜夜射夜夜干| 精品福利永久在线观看| 国产老妇伦熟女老妇高清| videosex国产| 亚洲色图综合在线观看| kizo精华| 欧美另类一区| 国产xxxxx性猛交| 成人国产av品久久久| 亚洲国产精品999| 国产亚洲午夜精品一区二区久久| 高清在线视频一区二区三区| 欧美精品国产亚洲| 国产伦理片在线播放av一区| 亚洲成人手机| 国产女主播在线喷水免费视频网站| 一级黄片播放器| 91午夜精品亚洲一区二区三区| 日本午夜av视频| 日本av手机在线免费观看| 80岁老熟妇乱子伦牲交| 久久久久视频综合| 精品99又大又爽又粗少妇毛片| 欧美日韩成人在线一区二区| 97人妻天天添夜夜摸| 亚洲国产精品一区二区三区在线| 最后的刺客免费高清国语| 亚洲中文av在线| 天美传媒精品一区二区| 亚洲国产精品999| 国精品久久久久久国模美| 最近中文字幕2019免费版| 成年人午夜在线观看视频| 欧美日韩成人在线一区二区| 高清毛片免费看| 亚洲av福利一区| 免费看不卡的av| 亚洲精品av麻豆狂野| 国产av码专区亚洲av| 午夜激情av网站| 中文欧美无线码| 自线自在国产av| 在线观看三级黄色| 久久av网站| 亚洲婷婷狠狠爱综合网| 成人毛片60女人毛片免费| 青春草国产在线视频| 成人国产麻豆网| 亚洲欧美色中文字幕在线| 狂野欧美激情性xxxx在线观看| 亚洲国产日韩一区二区| 2021少妇久久久久久久久久久| 久久久亚洲精品成人影院| 欧美成人精品欧美一级黄| 我的女老师完整版在线观看| 国产精品人妻久久久久久| 国产成人免费无遮挡视频| 亚洲欧美一区二区三区黑人 | 亚洲丝袜综合中文字幕| 日韩成人av中文字幕在线观看| 99国产综合亚洲精品| 欧美xxxx性猛交bbbb| 欧美日韩精品成人综合77777| 久久久a久久爽久久v久久| 亚洲欧洲精品一区二区精品久久久 | 精品卡一卡二卡四卡免费| 狠狠婷婷综合久久久久久88av| 性色avwww在线观看| 国产熟女欧美一区二区| 美国免费a级毛片| 大陆偷拍与自拍| 亚洲婷婷狠狠爱综合网| 久久女婷五月综合色啪小说| 久久韩国三级中文字幕| 精品酒店卫生间| 91精品伊人久久大香线蕉| 看十八女毛片水多多多| 国产成人a∨麻豆精品| 最近中文字幕2019免费版| 亚洲情色 制服丝袜| 欧美97在线视频| 又粗又硬又长又爽又黄的视频| av在线app专区| 亚洲国产精品专区欧美| 97在线人人人人妻| 草草在线视频免费看| 欧美成人精品欧美一级黄| 日韩伦理黄色片| 免费黄网站久久成人精品| 日韩不卡一区二区三区视频在线| 亚洲美女搞黄在线观看| 久久精品久久久久久噜噜老黄| 国产爽快片一区二区三区| 少妇人妻 视频| 极品少妇高潮喷水抽搐| 精品久久国产蜜桃| 日韩熟女老妇一区二区性免费视频| 精品一区二区三区视频在线| 中文字幕人妻熟女乱码| 天天影视国产精品| 最近中文字幕高清免费大全6| 国产免费福利视频在线观看| 亚洲国产欧美日韩在线播放| 久久久国产一区二区| 看免费成人av毛片| www日本在线高清视频| 哪个播放器可以免费观看大片| av在线老鸭窝| 午夜免费鲁丝| 国产成人精品无人区| 18禁动态无遮挡网站| 久久久国产一区二区| 国产一区二区激情短视频 | 成人综合一区亚洲| 免费观看无遮挡的男女| 天天躁夜夜躁狠狠久久av| 日韩制服丝袜自拍偷拍| 狂野欧美激情性xxxx在线观看| 99精国产麻豆久久婷婷| 热re99久久精品国产66热6| 十分钟在线观看高清视频www| 国产成人一区二区在线| 免费日韩欧美在线观看| 免费观看a级毛片全部| 美女中出高潮动态图| 黄片播放在线免费| videosex国产| av不卡在线播放| 22中文网久久字幕| 亚洲av综合色区一区| 人人妻人人澡人人看| 伊人久久国产一区二区| 一本色道久久久久久精品综合| 美女中出高潮动态图| 搡老乐熟女国产| 精品人妻熟女毛片av久久网站| av卡一久久| 狠狠婷婷综合久久久久久88av| 亚洲国产精品成人久久小说| 最新的欧美精品一区二区| 欧美日韩亚洲高清精品| 日本wwww免费看| 精品一区二区三卡| 欧美精品av麻豆av| 少妇高潮的动态图| 国产在视频线精品| 如何舔出高潮| 国产精品一二三区在线看| 午夜av观看不卡| 国产精品久久久久久久久免| 中国三级夫妇交换| 少妇的逼水好多| 卡戴珊不雅视频在线播放| 亚洲成人av在线免费| 久久久久国产精品人妻一区二区| 五月伊人婷婷丁香| 又黄又爽又刺激的免费视频.| 成人影院久久| 亚洲成av片中文字幕在线观看 | 亚洲熟女精品中文字幕| 精品亚洲成a人片在线观看| 国产亚洲欧美精品永久| 欧美激情国产日韩精品一区| 久久 成人 亚洲| 狠狠婷婷综合久久久久久88av| a级毛片黄视频| 亚洲天堂av无毛| 少妇精品久久久久久久| 国产精品免费大片| 91精品三级在线观看| 最后的刺客免费高清国语| 边亲边吃奶的免费视频| 妹子高潮喷水视频| 成人国产av品久久久| 成人毛片60女人毛片免费| 麻豆乱淫一区二区| 亚洲欧美成人综合另类久久久| 91成人精品电影| 精品久久国产蜜桃| 亚洲久久久国产精品| 最后的刺客免费高清国语| a级毛片在线看网站| 晚上一个人看的免费电影| 色吧在线观看| 欧美变态另类bdsm刘玥| 国精品久久久久久国模美| 日本黄大片高清|