亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于多端口波導(dǎo)結(jié)構(gòu)的寬頻帶聲觸發(fā)器*

        2023-09-06 02:15:28龐乃琦王垠葛勇施斌杰袁壽其孫宏祥2
        物理學(xué)報(bào) 2023年16期
        關(guān)鍵詞:觸發(fā)器聲壓聲學(xué)

        龐乃琦 王垠 葛勇 施斌杰 袁壽其 孫宏祥2)?

        1) (江蘇大學(xué)物理與電子工程學(xué)院,流體機(jī)械工程技術(shù)研究中心,鎮(zhèn)江 212013)

        2) (中國(guó)科學(xué)院聲學(xué)研究所,聲場(chǎng)聲信息國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100190)

        基于線性相干及相位調(diào)控機(jī)制設(shè)計(jì)制備了兩類聲學(xué)觸發(fā)器,所設(shè)計(jì)的聲學(xué)觸發(fā)器由相控單元和多端口波導(dǎo)結(jié)構(gòu)組成,其寬度與長(zhǎng)度分別為0.32λ 和0.82λ (λ 為聲波波長(zhǎng)),具有亞波長(zhǎng)結(jié)構(gòu)特征.基于相控單元的相位調(diào)制及聲波的線性相干機(jī)制,分別實(shí)現(xiàn)了T 觸發(fā)器和D 觸發(fā)器的聲學(xué)邏輯功能,且相對(duì)帶寬(工作帶寬與工作頻帶的中心頻率之比)分別可以達(dá)到0.23 和0.22.實(shí)驗(yàn)測(cè)量與數(shù)值模擬的結(jié)果吻合很好.本文所提出的聲觸發(fā)器具有寬頻帶、亞波長(zhǎng)尺寸及結(jié)構(gòu)簡(jiǎn)單等特點(diǎn),可為設(shè)計(jì)新型聲觸發(fā)器及聲邏輯門提供理論方案與原理性器件.

        1 引言

        聲邏輯基礎(chǔ)元器件由于其在邏輯運(yùn)算、信息處理和集成聲學(xué)等領(lǐng)域的廣泛應(yīng)用前景[1-7],受到研究學(xué)者越來(lái)越多的關(guān)注,并已發(fā)展出多種不同機(jī)制的聲邏輯門.如基于非線性機(jī)制,Li 等[8]基于球形粒子驅(qū)動(dòng)鏈的非線性動(dòng)力學(xué)效應(yīng)實(shí)現(xiàn)了與門(AND)和或門(OR)的聲邏輯功能.近年來(lái),聲子晶體[9-19]的迅猛發(fā)展為設(shè)計(jì)實(shí)現(xiàn)高性能聲邏輯門器件提供了可行性.例如Bringuier 等[11]基于聲子晶體中的輸入與控制信號(hào)的線性相干機(jī)制設(shè)計(jì)實(shí)現(xiàn)了與非門(NAND),異或門(XOR)和非門(NOT)等類型聲邏輯門.張婷等[12]基于包含線缺陷聲子晶體的自準(zhǔn)直聲束線性干涉,設(shè)計(jì)制備了AND,OR 和NOT等基礎(chǔ)聲邏輯門器件.夏建平等[14]提出了具有可編程功能的聲拓?fù)浣^緣體,實(shí)現(xiàn)了OR 和XOR 的聲邏輯功能.陸雨靜等[19]基于谷聲子晶體設(shè)計(jì)了一種X 型拓?fù)洳▽?dǎo)結(jié)構(gòu),基于聲波的線性相干及谷守恒機(jī)制,實(shí)現(xiàn)了具有魯棒性的OR 和XOR 雙功能聲邏輯門器件.然而,基于聲子晶體的聲邏輯門器件普遍存在著尺寸大問(wèn)題,為了克服該不足,研究人員基于聲學(xué)超材料[20-31]設(shè)計(jì)實(shí)現(xiàn)了各類具有寬頻帶和超薄特性的聲邏輯門器件.如張婷等[27]基于蜷曲空間單元結(jié)構(gòu)設(shè)計(jì)制備了具有超薄平面結(jié)構(gòu)的聲邏輯門器件,通過(guò)控制聲信號(hào)的初始相位實(shí)現(xiàn)各類基礎(chǔ)聲邏輯門功能.左承毅等[28]提出一種基于多端口圓波導(dǎo)的寬頻帶聲邏輯門,通過(guò)使用具有一定相位差的兩個(gè)輸入信號(hào)即可實(shí)現(xiàn)相關(guān)的聲邏輯功能,且所實(shí)現(xiàn)的工作帶寬大于5 kHz.在上述聲邏輯門中,研究人員均通過(guò)調(diào)制聲信號(hào)初始相位和幅值實(shí)現(xiàn)各類聲學(xué)邏輯功能.此外,王垠等[29]將相控單元置于多端口波導(dǎo)結(jié)構(gòu)中,基于被動(dòng)單元結(jié)構(gòu)調(diào)制聲信號(hào)的初始相位,同樣也可以實(shí)現(xiàn)各類聲學(xué)邏輯功能.上述工作均有效實(shí)現(xiàn)了各類聲學(xué)邏輯門器件,且有效地推動(dòng)了其應(yīng)用推廣.然而,對(duì)于邏輯功能更為復(fù)雜的聲觸發(fā)器基礎(chǔ)元器件,相關(guān)工作仍有待深入展開(kāi).

        本文提出了兩類基于多端口波導(dǎo)結(jié)構(gòu)的聲學(xué)觸發(fā)器.通過(guò)在波導(dǎo)結(jié)構(gòu)中嵌入相控單元,利用單元的相位調(diào)控及聲波線性相干機(jī)制,在同一閾值情形下實(shí)現(xiàn)了T 觸發(fā)器和D 觸發(fā)器的聲邏輯功能.在此基礎(chǔ)上,進(jìn)一步實(shí)驗(yàn)驗(yàn)證了T 觸發(fā)器和D 觸發(fā)器的性能,實(shí)驗(yàn)測(cè)量與數(shù)值模擬結(jié)果吻合較好.最后,實(shí)驗(yàn)測(cè)量了兩類聲觸發(fā)器的工作帶寬,所設(shè)計(jì)的T 觸發(fā)器和D 觸發(fā)器的相對(duì)帶寬(工作帶寬與工作頻帶的中心頻率之比)分別可以達(dá)到0.23和0.22,具有典型的寬頻帶特性.

        2 聲觸發(fā)器設(shè)計(jì)與性能

        2.1 數(shù)值模型

        為了設(shè)計(jì)聲學(xué)觸發(fā)器,引入文獻(xiàn)[32]中的相控單元結(jié)構(gòu).如圖1(a)所示,單元(長(zhǎng)和寬分別為l和h)由直空氣通道與10 組對(duì)稱的斜擋板(角度為θ)組成,壁厚為w,其結(jié)構(gòu)參數(shù)如下:l=λ/2,h=λ/10 及w=0.005λ,其中λ為聲波波長(zhǎng).單元的固體部分由環(huán)氧樹(shù)脂3 D 打印制備,可以滿足硬邊界條件.本文采用有限元多物理場(chǎng)耦合軟件COMSOL Multiphysics 中的壓力聲學(xué)模塊建立數(shù)值模型,模擬單元及相關(guān)的器件性能.在數(shù)值模型中,輸入與輸出端口的邊界設(shè)置為平面波輻射邊界,其他邊界均設(shè)置為硬聲場(chǎng)邊界,采用的空氣密度ρ=1.21 kg/m3,聲速c=343 m/s.這里,設(shè)入射聲波波長(zhǎng)λ為10 cm (對(duì)應(yīng)空氣中的頻率f0=3.43 kHz).

        圖1 (a) 相控單元示意圖;(b) 頻率為3.43 kHz 的聲波通過(guò)具有不同參數(shù)θ 的單元產(chǎn)生的相位延遲與透射率;(c) 聲觸發(fā)器示意圖.T 和Qn 端口的紅色箭頭表示輸入聲信號(hào),Qn+1 端口的藍(lán)色箭頭表示輸出聲信號(hào)Fig.1.(a) Schematic of a phased unit cell;(b) phase delays (blue solid line) and transmissions (red dashed line) of sound wave with frequency of 3.43 kHz caused by the phased unit cells with different values of θ;(c) schematic of an acoustic trigger.The red arrows at the ports T and Qn represent input sound signals,and the blue arrows at the port Qn+1 are output sound signals.

        圖1(b)顯示頻率為3.43 kHz 的聲波通過(guò)具有不同參數(shù)θ的單元產(chǎn)生的相位延遲及透射率,變化參數(shù)θ對(duì)應(yīng)的相位延遲(藍(lán)實(shí)線)可以覆蓋整個(gè)2π 區(qū)間,且在15° <θ< 90°范圍,聲透射率(紅虛線)在0.8 以上.這里,選擇θ=19°和57°對(duì)應(yīng)的相控單元I 和II (黑色空心點(diǎn))設(shè)計(jì)聲觸發(fā)器,其聲透射率分別為TI=0.92 和TII=0.96,相位差φII—φI=π.圖1(c)顯示所設(shè)計(jì)的聲觸發(fā)器示意圖,聲觸發(fā)器左右兩側(cè)輸入端T 和Qn對(duì)應(yīng)的直通道分別放置相控單元,中間區(qū)域填充環(huán)氧樹(shù)脂剛性固體,其中參數(shù)d1=0.075λ及d2=0.1λ.值得注意的是,聲觸發(fā)器長(zhǎng)度和寬度分別為0.82λ和0.32λ,具有亞波長(zhǎng)結(jié)構(gòu)特征.此外,當(dāng)T 和Qn端口對(duì)應(yīng)的通道放置相同相控單元時(shí),Qn+1端口的聲信號(hào)為干涉增強(qiáng);當(dāng)T 和Qn端口對(duì)應(yīng)的通道放置反相相控單元時(shí)(圖1(c)中單元I 和II),Qn+1端口的聲信號(hào)則為干涉相消.因此,基于線性聲干涉機(jī)制,可以實(shí)現(xiàn)多種聲觸發(fā)器功能.

        2.2 T 觸發(fā)器結(jié)構(gòu)與性能

        如圖1(c),T 觸發(fā)器的T 和Qn端口對(duì)應(yīng)的通道分別放置相控單元I 和II,令T 和Qn端口的輸入聲信號(hào)初始相位和振幅相同,聲觸發(fā)器的輸入態(tài)表示為(Q1,Q2),Qi表示T 和Qn端口有無(wú)輸入信號(hào),分別定義為“1”和“0”.T 觸發(fā)器主要用于實(shí)現(xiàn)保持與反轉(zhuǎn)功能[33,34],其特性方程滿足Qn+1=T⊕Qn.圖2(a)顯示數(shù)值模擬不同輸入態(tài)激發(fā)T 觸發(fā)器產(chǎn)生的聲壓幅值場(chǎng)分布,當(dāng)輸入態(tài)為(1,1)時(shí),Qn+1端口的聲壓幅值接近于0.這主要由于所采用的是一對(duì)反相相控單元,聲信號(hào)在Qn+1端口產(chǎn)生干涉相消.然而,當(dāng)輸入態(tài)為(1,0)和(0,1)時(shí),Qn+1端口的聲壓幅值較強(qiáng),明顯高于輸入態(tài)(1,1)對(duì)應(yīng)的輸出聲信號(hào).

        圖2 (a) 數(shù)值模擬頻率為3.43 kHz 不同輸入態(tài)激發(fā)T 觸發(fā)器產(chǎn)生的聲壓幅值場(chǎng)分布;(b),(c) 對(duì)應(yīng)的輸出端Qn+1 的聲能級(jí)和真值表Fig.2.(a) Simulated pressure amplitude distributions caused by the T-type trigger with different input states at 3.43 kHz;(b),(c) simulated acoustic intensity levels at the output port Qn+1 and truth table.

        為了量化聲觸發(fā)器的輸出性能,引入聲能級(jí)定義

        式中p1和p2分別表示有無(wú)聲觸發(fā)器結(jié)構(gòu)對(duì)應(yīng)的出射端區(qū)域S中的聲壓.這里,選取閾值It=10 dB,當(dāng)I>It時(shí),其輸出態(tài)為“1”;而I<It時(shí),其輸出態(tài)為“0”.圖2(b)和圖2(c)分別顯示數(shù)值模擬T 觸發(fā)器輸出端Qn+1的聲能級(jí)及其真值表.基于所選擇的閾值It=10 dB,可以得到,輸入態(tài)(1,1),(1,0),(0,1)及(0,0)對(duì)應(yīng)的Qn+1端口輸出態(tài)分別為(1),(1),(1)及(0).如圖2(c)真值表所示,當(dāng)T 端口的輸入態(tài)為(0)時(shí),Qn+1端口輸出態(tài)與Qn端口輸入態(tài)保持一致;而當(dāng)T 端口的輸入態(tài)為(1)時(shí),Qn+1端口輸出態(tài)與Qn端口輸入態(tài)相反,實(shí)現(xiàn)了T 觸發(fā)器的聲邏輯功能.

        2.3 D 觸發(fā)器結(jié)構(gòu)與性能

        在T 觸發(fā)器結(jié)構(gòu)的中間固體區(qū)域引入控制端A,可以設(shè)計(jì)實(shí)現(xiàn)D 觸發(fā)器.D 觸發(fā)器主要用于實(shí)現(xiàn)置0 和置1 功能[33,34],其特性方程滿足 Qn+1=D .下面為了實(shí)現(xiàn)控制端A 的功能,數(shù)值模擬聲波通過(guò)具有不同參數(shù)θ的單元(斜擋板寬度為3w)產(chǎn)生的相位延遲 (藍(lán)實(shí)線)及透射率(紅虛線),如圖3 所示,選取相控單元III 放置控制端A,其擋板角度θ=52°,聲透射率TIII=0.7,相位延遲φIII=φII— π/5.

        圖3 數(shù)值模擬聲波通過(guò)具有不同參數(shù)θ 的單元(斜擋板寬度為3w)產(chǎn)生的相位延遲 (藍(lán)色實(shí)線)及透射率(紅色虛線)Fig.3.Simulated phase delays (blue solid line) and transmissions (red dashed line) caused by the phased unit cells with different values of θ.

        圖4(a)給出了數(shù)值模擬不同輸入態(tài)激發(fā)D 觸發(fā)器產(chǎn)生的聲壓幅值場(chǎng)分布.可以看出,輸入態(tài)為(1,1)時(shí),Qn+1端口的輸出聲壓幅值較強(qiáng),與圖2(a)中輸入態(tài)(1,1)對(duì)應(yīng)的結(jié)果相反,這主要由于控制端A 的聲信號(hào)的存在.此外,輸入態(tài)(0,1)對(duì)應(yīng)的輸出聲壓幅值高于輸入態(tài)(1,1)對(duì)應(yīng)的結(jié)果,這主要由于D 端口的輸入態(tài)為(0),因此無(wú)信號(hào)與輸入端Qn的信號(hào)產(chǎn)生干涉相消.然而,當(dāng)輸入態(tài)為(1,0)和(0,0)時(shí),Qn+1端口的聲壓幅值較弱.前者是由于控制端A 和輸入端D 的聲信號(hào)在輸出端處產(chǎn)生弱干涉相消;而后者是由于D 與Qn端口無(wú)聲信號(hào)輸入,而控制端A 通道中的單元III 聲透射率相對(duì)較低(約為0.7),且透射聲信號(hào)能量會(huì)均勻分散到其他三個(gè)端口,因此輸出端信號(hào)的聲能級(jí)低于閾值.圖4(b)和圖4(c)分別顯示輸出端Qn+1的聲能級(jí)及其真值表.基于統(tǒng)一的閾值,輸入態(tài)(1,1),(1,0),(0,1)及(0,0)對(duì)應(yīng)的Qn+1端口輸出態(tài)分別為(1),(0),(1)及(0).如圖4(c)所示,Qn+1端口的輸出態(tài)始終與D 端口的輸入態(tài)保持一致,而與Qn端口的輸入態(tài)無(wú)關(guān),從而實(shí)現(xiàn)了D 觸發(fā)器相關(guān)的聲邏輯功能.

        圖4 (a) 數(shù)值模擬頻率為3.43 kHz 不同輸入態(tài)激發(fā)D 觸發(fā)器產(chǎn)生的聲壓幅值場(chǎng)分布;(b),(c)對(duì)應(yīng)的輸出端Qn+1 的聲能級(jí)和真值表Fig.4.(a) Simulated pressure amplitude distributions caused by the D-type trigger with different input states at 3.43 kHz;(b),(c) simulated acoustic intensity levels at the output port Qn+1 and truth table.

        3 實(shí)驗(yàn)結(jié)果與討論

        為了驗(yàn)證T 觸發(fā)器和D 觸發(fā)器的聲學(xué)性能,實(shí)驗(yàn)測(cè)量?jī)深惵曈|發(fā)器輸出端Qn+1的聲能級(jí)強(qiáng)度.圖5(a)為實(shí)驗(yàn)測(cè)量裝置,圖5(b)和圖5(c)為樣品照片.為了實(shí)現(xiàn)相同輸入聲信號(hào),在左側(cè)兩個(gè)輸入端放置相同的聲源,相控單元與聲源的間距為5 mm.輸入聲信號(hào)由功率放大器驅(qū)動(dòng)聲源產(chǎn)生.采用0.25 in (1 in=2.54 cm)麥克風(fēng)(Brüel&Kj?r-4961 型) 在右側(cè)輸出端測(cè)量輸出聲信號(hào).實(shí)驗(yàn)測(cè)量數(shù)據(jù)由Brüel&Kj?r3160-A-022 模塊進(jìn)行采集記錄,并通過(guò)Pulse Labshop 軟件進(jìn)行數(shù)據(jù)分析.圖5(d),(e)顯示實(shí)驗(yàn)測(cè)量的頻率為3.430 kHz 的不同輸入態(tài)聲波激發(fā)兩類聲觸發(fā)器對(duì)應(yīng)的輸出端Qn+1的聲能級(jí),其中參考信號(hào)為實(shí)驗(yàn)測(cè)量的環(huán)境噪聲,與圖2(b)和圖4(b)中的數(shù)值模擬結(jié)果相比,實(shí)驗(yàn)測(cè)量與數(shù)值模擬結(jié)果吻合較好,從而驗(yàn)證了所設(shè)計(jì)的T 觸發(fā)器和D 觸發(fā)器可行性.

        圖5 (a) 實(shí)驗(yàn)裝置示意圖;(b),(c) T 型與D 型觸發(fā)器樣品照片;(d),(e) 實(shí)驗(yàn)測(cè)量的頻率為3.43 kHz 的不同輸入態(tài)聲波激發(fā)T 觸發(fā)器和D 觸發(fā)器對(duì)應(yīng)輸出端Qn+1 的聲能級(jí)Fig.5.(a) Schematic of experimental set-up;(b),(c) photographs of the T-type trigger and D-type trigger;(d),(e) experimental measurement of the acoustic intensity levels at the output port Qn+1 of T-type trigger and D-type trigger at a frequency of 3.43 kHz.

        為了進(jìn)一步展示所設(shè)計(jì)的T 觸發(fā)器和D 觸發(fā)器的帶寬,實(shí)驗(yàn)測(cè)量?jī)深愑|發(fā)器輸出端Qn+1對(duì)應(yīng)的聲能級(jí)譜,分別見(jiàn)圖6(a)和圖6(b).可以看出,基于統(tǒng)一閾值It=10 dB,在陰影區(qū)域范圍中均可實(shí)現(xiàn)T 觸發(fā)器和D 觸發(fā)器的功能,對(duì)應(yīng)的工作頻帶分別為3.293—4.069 kHz 和3.400—4.138 kHz,相對(duì)帶寬分別可以達(dá)到0.23 和0.22,從而實(shí)驗(yàn)驗(yàn)證了兩類聲觸發(fā)器的寬頻帶特征.除此之外,如圖6(a)和圖6(b)所示,兩類觸發(fā)器不同輸入態(tài)所對(duì)應(yīng)的聲能級(jí)譜變化幅度較大,這主要由于所設(shè)計(jì)的單元相位差與頻率密切相關(guān),同時(shí)也限制了觸發(fā)器的工作帶寬.在此基礎(chǔ)上,可知兩類觸發(fā)器的工作帶寬不一致,這主要由于D 觸發(fā)器控制端A 中的單元III 相位延遲及聲透射率同樣與頻率相關(guān),從而影響了D 觸發(fā)器的工作帶寬.

        圖6 T 觸發(fā)器(a)和D 觸發(fā)器(b)輸出端處的不同輸入態(tài)對(duì)應(yīng)的聲能級(jí)譜.黑色陰影區(qū)域范圍分別為(a) 3.293—4.069 kHz,(b) 3.400—4.138 kHzFig.6.Measured intensity level spectra at the output ports of the T-type trigger (a) and D-type trigger (b) for different input states.Black shaded regions cover the ranges of 3.293—4.069 kHz in panel (a) and 3.400—4.138 kHz in panel(b).

        4 結(jié)論

        本文基于多端口波導(dǎo)結(jié)構(gòu)及相控單元,實(shí)驗(yàn)設(shè)計(jì)了兩類聲學(xué)觸發(fā)器,所設(shè)計(jì)的聲學(xué)觸發(fā)器寬度和長(zhǎng)度分別為0.32λ和0.82λ,具有亞波長(zhǎng)結(jié)構(gòu)特征.研究結(jié)果表明: 在多端口波導(dǎo)結(jié)構(gòu)中嵌入相控單元,利用單元的相位調(diào)控及聲波的線性相干機(jī)制,可以實(shí)現(xiàn)T 觸發(fā)器和D 觸發(fā)器對(duì)應(yīng)的聲邏輯功能,且對(duì)應(yīng)的相對(duì)帶寬分別可以達(dá)到0.23 和0.22,具有寬頻帶特征.實(shí)驗(yàn)測(cè)量與數(shù)值模擬結(jié)果吻合較好.與其他類型的聲邏輯器件相比,所設(shè)計(jì)的聲觸發(fā)器具有具有小型化、寬頻帶、相同閾值及波導(dǎo)結(jié)構(gòu)易連接等優(yōu)點(diǎn),在聲開(kāi)關(guān)、聲通信及聲二極管等領(lǐng)域具有一定的潛在應(yīng)用,同時(shí)也為設(shè)計(jì)新型聲邏輯器件提供了理論方案與原理性器件.

        猜你喜歡
        觸發(fā)器聲壓聲學(xué)
        基于嘴唇處的聲壓數(shù)據(jù)確定人體聲道半徑
        愛(ài)的就是這股Hi-Fi味 Davis Acoustics(戴維斯聲學(xué))Balthus 70
        車輛結(jié)構(gòu)噪聲傳遞特性及其峰值噪聲成因的分析
        汽車工程(2018年12期)2019-01-29 06:46:36
        Acoustical Treatment Primer:Diffusion談?wù)劼晫W(xué)處理中的“擴(kuò)散”
        Acoustical Treatment Primer:Absorption談?wù)劼晫W(xué)處理中的“吸聲”(二)
        Acoustical Treatment Primer:Absorption 談?wù)劼晫W(xué)處理中的“吸聲”
        主從JK觸發(fā)器邏輯功能分析
        電子世界(2017年22期)2017-12-02 03:03:45
        使用觸發(fā)器,強(qiáng)化安全性
        基于GIS內(nèi)部放電聲壓特性進(jìn)行閃絡(luò)定位的研究
        基于聲壓原理的柴油發(fā)動(dòng)機(jī)檢測(cè)室噪聲的測(cè)量、分析與治理
        日本成人在线不卡一区二区三区| 亚洲乱码日产精品一二三| 亚洲av无码成人网站在线观看| 中国农村熟妇性视频| 色一乱一伦一图一区二区精品| 久久久精品电影| 国产毛片一区二区日韩| 日本高清不卡二区三区| 日韩精品视频免费网站| 国产精品 无码专区| 午夜福利试看120秒体验区| 国产肉丝袜在线观看| 亚洲av无码专区在线亚| 看全色黄大黄大色免费久久 | 国产精品亚洲A∨无码遮挡| 国产目拍亚洲精品二区| 美女脱掉内裤扒开下面让人插| 日韩一区二区超清视频| 天天鲁一鲁摸一摸爽一爽| 国产成人精品无码播放| 亚洲欧美久久婷婷爱综合一区天堂| 性视频毛茸茸女性一区二区| 最近更新中文字幕一区二区| 久久久久亚洲av成人人电影| 国产亚洲2021成人乱码| 巨熟乳波霸若妻在线播放| 日韩人妻无码精品二专区| 一区二区三区日韩毛片| 国产夫妻自拍视频在线播放 | 亚洲深深色噜噜狠狠网站| 国产熟女内射oooo| 国产福利姬喷水福利在线观看| 国产午夜亚洲精品一级在线| 色婷婷一区二区三区四| 那有一级内射黄片可以免费看| 国产电影无码午夜在线播放| 国产乱妇乱子在线播视频播放网站| 久久精品—区二区三区无码伊人色| 久久精品国产亚洲精品色婷婷| 精品精品国产三级av在线| 精人妻无码一区二区三区|