• <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)用
    真人一进一出gif抽搐免费| 俄罗斯特黄特色一大片| 美女高潮喷水抽搐中文字幕| 九色成人免费人妻av| 亚洲最大成人av| 乱人视频在线观看| 亚洲五月天丁香| 给我免费播放毛片高清在线观看| 欧美丝袜亚洲另类 | 久久欧美精品欧美久久欧美| 国产一区二区激情短视频| 亚洲无线在线观看| 亚洲精华国产精华精| 久久久色成人| 毛片女人毛片| 国产v大片淫在线免费观看| 校园春色视频在线观看| 国产v大片淫在线免费观看| 色综合婷婷激情| 久久久久国产精品人妻aⅴ院| 亚洲人成网站在线播放欧美日韩| av在线观看视频网站免费| 国产国拍精品亚洲av在线观看| 我要搜黄色片| 亚洲欧美清纯卡通| 中出人妻视频一区二区| 高清在线国产一区| av天堂中文字幕网| 亚洲精品日韩av片在线观看| 欧美黑人欧美精品刺激| 久久6这里有精品| 亚洲最大成人手机在线| 又黄又爽又刺激的免费视频.| 搡老熟女国产l中国老女人| 一个人看的www免费观看视频| 精品一区二区三区人妻视频| 国产精品久久久久久亚洲av鲁大| av在线天堂中文字幕| 久久天躁狠狠躁夜夜2o2o| 国产亚洲av嫩草精品影院| 久久久久性生活片| 又爽又黄a免费视频| 日韩欧美在线乱码| 亚洲精品在线美女| 波多野结衣高清无吗| 国产av在哪里看| 亚洲美女视频黄频| 日本熟妇午夜| 亚洲三级黄色毛片| 成熟少妇高潮喷水视频| 观看美女的网站| 精品一区二区免费观看| 日韩成人在线观看一区二区三区| 亚洲中文字幕日韩| 精品不卡国产一区二区三区| 中出人妻视频一区二区| 99精品久久久久人妻精品| 国产精品嫩草影院av在线观看 | 一个人观看的视频www高清免费观看| 亚洲美女黄片视频| 三级国产精品欧美在线观看| 欧美最新免费一区二区三区 | 日本熟妇午夜| 午夜日韩欧美国产| 夜夜夜夜夜久久久久| 俺也久久电影网| 五月玫瑰六月丁香| 日本免费一区二区三区高清不卡| 婷婷精品国产亚洲av| 亚洲成人免费电影在线观看| 亚洲成人精品中文字幕电影| 88av欧美| 亚洲电影在线观看av| 欧美最黄视频在线播放免费| 熟女电影av网| 亚洲自拍偷在线| 哪里可以看免费的av片| 欧美一区二区亚洲| 亚洲人成网站在线播| 波多野结衣巨乳人妻| 小说图片视频综合网站| 欧美在线黄色| 最后的刺客免费高清国语| aaaaa片日本免费| 男女那种视频在线观看| 亚洲av五月六月丁香网| 丰满人妻一区二区三区视频av| 中文字幕熟女人妻在线| 亚洲av熟女| 在线观看午夜福利视频| 精品人妻熟女av久视频| 熟女人妻精品中文字幕| 亚洲最大成人手机在线| 日韩亚洲欧美综合| 欧美乱妇无乱码| 18禁裸乳无遮挡免费网站照片| avwww免费| 大型黄色视频在线免费观看| 九色国产91popny在线| a在线观看视频网站| 有码 亚洲区| 久久精品夜夜夜夜夜久久蜜豆| 欧美最新免费一区二区三区 | 两性午夜刺激爽爽歪歪视频在线观看| 一本综合久久免费| 久久人人精品亚洲av| 最新中文字幕久久久久| 波多野结衣高清作品| 欧美性猛交╳xxx乱大交人| 久久久久久久精品吃奶| xxxwww97欧美| 欧美乱色亚洲激情| 精品国产亚洲在线| 成年女人永久免费观看视频| 高清毛片免费观看视频网站| 国内久久婷婷六月综合欲色啪| 天堂av国产一区二区熟女人妻| 搡老熟女国产l中国老女人| 亚洲精品一区av在线观看| 亚洲自偷自拍三级| 午夜老司机福利剧场| 国产大屁股一区二区在线视频| 日韩国内少妇激情av| 成人国产一区最新在线观看| 我的女老师完整版在线观看| 别揉我奶头~嗯~啊~动态视频| av黄色大香蕉| 成人永久免费在线观看视频| 国产精品久久久久久亚洲av鲁大| 精品乱码久久久久久99久播| 天美传媒精品一区二区| 很黄的视频免费| 国产成人a区在线观看| 欧美精品啪啪一区二区三区| 欧美又色又爽又黄视频| 男女视频在线观看网站免费| 国产视频内射| 国产伦一二天堂av在线观看| 成人av一区二区三区在线看| 精品人妻偷拍中文字幕| 亚洲成av人片在线播放无| 成人av在线播放网站| 国产日本99.免费观看| 色精品久久人妻99蜜桃| 欧美黑人欧美精品刺激| 国产精品一区二区三区四区久久| 草草在线视频免费看| 少妇丰满av| www日本黄色视频网| 2021天堂中文幕一二区在线观| 亚洲美女黄片视频| 亚洲国产日韩欧美精品在线观看| 中文字幕久久专区| 亚洲男人的天堂狠狠| 国产成人啪精品午夜网站| 日本三级黄在线观看| 亚洲成av人片在线播放无| 国产伦一二天堂av在线观看| 国产在视频线在精品| 国产亚洲欧美98| 亚洲欧美日韩高清专用| 久久久色成人| av女优亚洲男人天堂| 午夜福利在线观看吧| 美女cb高潮喷水在线观看| 午夜免费激情av| 一级作爱视频免费观看| 国产日本99.免费观看| 精品久久久久久久人妻蜜臀av| 国产午夜福利久久久久久| 国产伦精品一区二区三区视频9| 少妇裸体淫交视频免费看高清| 久久久久久久精品吃奶| 99热6这里只有精品| 亚洲在线观看片| 一区二区三区激情视频| 成年女人永久免费观看视频| 一进一出抽搐动态| 日本三级黄在线观看| 免费黄网站久久成人精品 | 青草久久国产| 亚洲av免费高清在线观看| 桃红色精品国产亚洲av| 性插视频无遮挡在线免费观看| 精品久久久久久久久av| 精品久久久久久久久亚洲 | 欧美日韩乱码在线| 欧美一区二区亚洲| 丝袜美腿在线中文| 久99久视频精品免费| bbb黄色大片| 91麻豆精品激情在线观看国产| 国产精品99久久久久久久久| 亚洲经典国产精华液单 | 亚洲三级黄色毛片| 国产在线男女| 亚洲无线观看免费| 国产在线男女| 午夜精品一区二区三区免费看| 好看av亚洲va欧美ⅴa在| 国产午夜精品论理片| 国语自产精品视频在线第100页| 国产精品1区2区在线观看.| 一区二区三区激情视频| 亚洲,欧美,日韩| 俄罗斯特黄特色一大片| 亚洲狠狠婷婷综合久久图片| 男女之事视频高清在线观看| 成人性生交大片免费视频hd| 床上黄色一级片| 老司机福利观看| 亚洲成av人片免费观看| 国产淫片久久久久久久久 | 精品久久久久久久久久久久久| 国产69精品久久久久777片| 欧美另类亚洲清纯唯美| 国产精品久久电影中文字幕| 欧美3d第一页| 蜜桃久久精品国产亚洲av| 天堂影院成人在线观看| 亚洲一区高清亚洲精品| 成年人黄色毛片网站| 国产av在哪里看| 中文亚洲av片在线观看爽| 国产精品野战在线观看| 亚洲专区国产一区二区| 欧美日韩黄片免| 国产精品久久久久久精品电影| 欧美性猛交黑人性爽| 色精品久久人妻99蜜桃| 欧美xxxx性猛交bbbb| 99热这里只有精品一区| 如何舔出高潮| 麻豆国产97在线/欧美| 我要搜黄色片| 欧美成人性av电影在线观看| 一个人观看的视频www高清免费观看| 亚洲五月婷婷丁香| 男人的好看免费观看在线视频| 热99re8久久精品国产| 亚洲五月天丁香| www.色视频.com| 丰满乱子伦码专区| 偷拍熟女少妇极品色| 性色avwww在线观看| 欧美另类亚洲清纯唯美| 搡老熟女国产l中国老女人| 国产极品精品免费视频能看的| 亚洲美女视频黄频| 精品久久国产蜜桃| 久久精品久久久久久噜噜老黄 | 精品久久久久久成人av| 乱码一卡2卡4卡精品| 色av中文字幕| 色尼玛亚洲综合影院| 中文字幕久久专区| 日本熟妇午夜| 成人高潮视频无遮挡免费网站| 国产亚洲精品久久久久久毛片| 伦理电影大哥的女人| 男人舔奶头视频| 亚洲av免费在线观看| 亚洲国产精品成人综合色| 在线免费观看不下载黄p国产 | 我要看日韩黄色一级片| 亚洲五月天丁香| 亚洲黑人精品在线| 日韩欧美精品v在线| 久久久久久九九精品二区国产| 久久久久久久久大av| 一二三四社区在线视频社区8| 亚洲自偷自拍三级| 亚洲真实伦在线观看| 搡女人真爽免费视频火全软件 | 国产日本99.免费观看| 亚洲一区二区三区不卡视频| www日本黄色视频网| 亚洲黑人精品在线| 欧美一级a爱片免费观看看| 国产精品野战在线观看| 国产精品久久久久久亚洲av鲁大| 欧美日韩中文字幕国产精品一区二区三区| 中出人妻视频一区二区| 国产亚洲精品久久久久久毛片| 国产私拍福利视频在线观看| 欧美乱妇无乱码| 久久伊人香网站| 日本在线视频免费播放| 国产白丝娇喘喷水9色精品| 国产男靠女视频免费网站| 亚洲七黄色美女视频| 亚洲不卡免费看| .国产精品久久| 美女黄网站色视频| 免费看美女性在线毛片视频| 国产人妻一区二区三区在| 国产精品亚洲美女久久久| 婷婷精品国产亚洲av| 精品一区二区三区视频在线观看免费| 免费在线观看影片大全网站| 美女大奶头视频| 亚洲国产高清在线一区二区三| 丁香欧美五月| 深爱激情五月婷婷| 免费搜索国产男女视频| 亚洲黑人精品在线| 亚洲五月婷婷丁香| 国产精品亚洲av一区麻豆| 久久久久久久亚洲中文字幕 | 天天躁日日操中文字幕| 中出人妻视频一区二区| 亚洲电影在线观看av| 国产精品爽爽va在线观看网站| 91麻豆av在线| 午夜福利在线观看免费完整高清在 | 高清日韩中文字幕在线| 成人国产综合亚洲| 亚洲成人免费电影在线观看| 国产三级中文精品| 一级a爱片免费观看的视频| 精品人妻视频免费看| 九九热线精品视视频播放| 日韩有码中文字幕| 国产成人啪精品午夜网站| 九色国产91popny在线| 日日干狠狠操夜夜爽| 黄片小视频在线播放| 最近中文字幕高清免费大全6 | 亚洲国产欧洲综合997久久,| 国产精品,欧美在线| 男女下面进入的视频免费午夜| avwww免费| 97超视频在线观看视频| 亚洲18禁久久av| 宅男免费午夜| 搡老岳熟女国产| 少妇的逼水好多| 国产伦在线观看视频一区| 亚洲 欧美 日韩 在线 免费| 久久这里只有精品中国| 性色av乱码一区二区三区2| 三级毛片av免费| 亚洲不卡免费看| 在线十欧美十亚洲十日本专区| 国产白丝娇喘喷水9色精品| 精品一区二区免费观看| 无人区码免费观看不卡| 亚洲av二区三区四区| 久久精品影院6| 亚洲五月婷婷丁香| 国产精品一区二区三区四区免费观看 | 精品人妻偷拍中文字幕| 一区二区三区免费毛片| 国产探花在线观看一区二区| eeuss影院久久| 日本 欧美在线| 99久久99久久久精品蜜桃| 我的女老师完整版在线观看| 久久草成人影院| 身体一侧抽搐| 天堂√8在线中文| 非洲黑人性xxxx精品又粗又长| h日本视频在线播放| 丰满乱子伦码专区| 国产av一区在线观看免费| 久久欧美精品欧美久久欧美| 国产亚洲精品久久久com| 久久久久久久久久成人| 少妇被粗大猛烈的视频| 久久久久久久精品吃奶| 99久久精品一区二区三区| 亚洲熟妇熟女久久| 精品人妻一区二区三区麻豆 | 亚洲欧美清纯卡通| 九九久久精品国产亚洲av麻豆| 日韩欧美三级三区| 国产精品久久视频播放| 国产成人aa在线观看| 欧美潮喷喷水| 97超级碰碰碰精品色视频在线观看| 国产av麻豆久久久久久久| 日本免费一区二区三区高清不卡| 一二三四社区在线视频社区8| 国模一区二区三区四区视频| 欧洲精品卡2卡3卡4卡5卡区| 婷婷精品国产亚洲av| 亚洲精品在线观看二区| 深爱激情五月婷婷| 欧美最黄视频在线播放免费| 精品国产三级普通话版| 日本撒尿小便嘘嘘汇集6| 婷婷精品国产亚洲av| 久久人人爽人人爽人人片va | 18禁黄网站禁片午夜丰满| 网址你懂的国产日韩在线| 国产一区二区激情短视频| 欧美高清性xxxxhd video| 日本黄色视频三级网站网址| 99在线人妻在线中文字幕| 99久久无色码亚洲精品果冻| 中文在线观看免费www的网站| 激情在线观看视频在线高清| av在线观看视频网站免费| 成人特级av手机在线观看| 色精品久久人妻99蜜桃| 热99re8久久精品国产| 9191精品国产免费久久| 三级男女做爰猛烈吃奶摸视频| 午夜福利18| 观看美女的网站| 亚洲久久久久久中文字幕| 久久草成人影院| 国产成人影院久久av| 高潮久久久久久久久久久不卡| 成人高潮视频无遮挡免费网站| 国产精品一及| 日韩有码中文字幕| 国产精品人妻久久久久久| 2021天堂中文幕一二区在线观| 观看美女的网站| 毛片女人毛片| 精品久久久久久久人妻蜜臀av| 2021天堂中文幕一二区在线观| 69av精品久久久久久| 午夜福利成人在线免费观看| 亚洲一区高清亚洲精品| 日韩欧美在线二视频| 欧美日韩综合久久久久久 | 高清毛片免费观看视频网站| 99久久成人亚洲精品观看| 国产一级毛片七仙女欲春2| 久久99热这里只有精品18| 国语自产精品视频在线第100页| 亚洲自偷自拍三级| 色播亚洲综合网| 久久精品国产亚洲av天美| 亚洲国产高清在线一区二区三| 亚洲av第一区精品v没综合| 99久久成人亚洲精品观看| 网址你懂的国产日韩在线| 欧美一区二区国产精品久久精品| 国内久久婷婷六月综合欲色啪| 999久久久精品免费观看国产| 十八禁人妻一区二区| 麻豆成人av在线观看| 亚洲国产高清在线一区二区三| 亚洲精品久久国产高清桃花| 亚洲一区二区三区不卡视频| 亚洲自偷自拍三级| av黄色大香蕉| 麻豆成人av在线观看| 久久久久国产精品人妻aⅴ院| 欧美精品国产亚洲| 亚洲人成网站在线播| 国产男靠女视频免费网站| 成人av在线播放网站| 欧美日韩黄片免| 亚洲精华国产精华精| 少妇熟女aⅴ在线视频| 国产亚洲欧美在线一区二区| 十八禁国产超污无遮挡网站| 国产午夜精品久久久久久一区二区三区 | www.熟女人妻精品国产| 床上黄色一级片| 欧美日韩国产亚洲二区| 久久午夜福利片| 日本三级黄在线观看| 久久精品夜夜夜夜夜久久蜜豆| 最近最新中文字幕大全电影3| 国产一区二区亚洲精品在线观看| 少妇熟女aⅴ在线视频| a在线观看视频网站| 精品人妻偷拍中文字幕| 国产视频内射| 亚洲美女视频黄频| 国产精品电影一区二区三区| 性色av乱码一区二区三区2| 一级黄色大片毛片| 国产成人av教育| 日韩中字成人| 久久国产精品人妻蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩瑟瑟在线播放| 国产精品嫩草影院av在线观看 | 国内揄拍国产精品人妻在线| 中文字幕久久专区| 男女之事视频高清在线观看| 亚洲电影在线观看av| 很黄的视频免费| 国内揄拍国产精品人妻在线| 免费人成在线观看视频色| 观看美女的网站| 亚洲人成网站在线播放欧美日韩| 99热6这里只有精品| 午夜a级毛片| 欧美最黄视频在线播放免费| 日日干狠狠操夜夜爽| 亚洲第一区二区三区不卡| 欧美色欧美亚洲另类二区| 草草在线视频免费看| 精品久久久久久久久久免费视频| 美女xxoo啪啪120秒动态图 | 国产麻豆成人av免费视频| 婷婷色综合大香蕉| 高潮久久久久久久久久久不卡| 我要搜黄色片| 精品国产亚洲在线| 成人鲁丝片一二三区免费| 国产91精品成人一区二区三区| 色综合婷婷激情| 欧美成人a在线观看| 午夜福利高清视频| 99久国产av精品| 国产精品免费一区二区三区在线| 91九色精品人成在线观看| 国模一区二区三区四区视频| 久久久久亚洲av毛片大全| 别揉我奶头 嗯啊视频| 久久久久久国产a免费观看| 黄色配什么色好看| 99精品久久久久人妻精品| 五月伊人婷婷丁香| 脱女人内裤的视频| 欧美zozozo另类| 国产高潮美女av| 国产精品伦人一区二区| 精品午夜福利在线看| 搡老熟女国产l中国老女人| 18禁裸乳无遮挡免费网站照片| av天堂在线播放| 精品人妻一区二区三区麻豆 | 村上凉子中文字幕在线| 九色国产91popny在线| 深爱激情五月婷婷| 日韩国内少妇激情av| 99国产精品一区二区蜜桃av| 成年女人毛片免费观看观看9| 久久香蕉精品热| 99riav亚洲国产免费| 日韩中文字幕欧美一区二区| 禁无遮挡网站| 91久久精品电影网| 18禁黄网站禁片午夜丰满| 成人国产一区最新在线观看| 在线看三级毛片| 亚洲精品一区av在线观看| 桃色一区二区三区在线观看| 亚洲av成人不卡在线观看播放网| 久99久视频精品免费| 亚洲精品在线观看二区| 欧美中文日本在线观看视频| 国产一区二区在线观看日韩| 国产精品久久电影中文字幕| 久久精品久久久久久噜噜老黄 | 88av欧美| 高潮久久久久久久久久久不卡| 国产一区二区在线av高清观看| 99热精品在线国产| 欧美黄色片欧美黄色片| 老女人水多毛片| 长腿黑丝高跟| 中文字幕高清在线视频| 免费看美女性在线毛片视频| 少妇熟女aⅴ在线视频| 国产精品久久视频播放| 亚洲国产高清在线一区二区三| 麻豆国产av国片精品| 一进一出抽搐gif免费好疼| 日韩 亚洲 欧美在线| 大型黄色视频在线免费观看| 久久久久性生活片| 五月伊人婷婷丁香| 国产在视频线在精品| 哪里可以看免费的av片| 在线播放无遮挡| 欧美zozozo另类| 全区人妻精品视频| 别揉我奶头~嗯~啊~动态视频| 日本五十路高清| 蜜桃久久精品国产亚洲av| 免费av毛片视频| 国产色婷婷99| 极品教师在线免费播放| 欧美日韩乱码在线| 91在线观看av| 色播亚洲综合网| 成人特级黄色片久久久久久久| 亚洲欧美清纯卡通| ponron亚洲| 女人被狂操c到高潮| 亚洲国产精品999在线| 国产69精品久久久久777片| 最新中文字幕久久久久| 国产aⅴ精品一区二区三区波| 一级a爱片免费观看的视频| 尤物成人国产欧美一区二区三区| 国产国拍精品亚洲av在线观看| 免费观看精品视频网站| 听说在线观看完整版免费高清| 黄色配什么色好看| 午夜免费激情av| 又粗又爽又猛毛片免费看| 国产单亲对白刺激| 成人国产综合亚洲| 国产成年人精品一区二区| 可以在线观看毛片的网站| 少妇的逼好多水| 9191精品国产免费久久| 国产高清有码在线观看视频| 久久香蕉精品热| 少妇人妻精品综合一区二区 | 亚洲一区二区三区色噜噜| 美女 人体艺术 gogo| av欧美777| 色5月婷婷丁香|