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

    人體整耳傳聲有限元數(shù)值模擬

    2016-05-19 09:09:05姚海峰劉后廣饒柱石黃新生中國礦業(yè)大學(xué)機(jī)電工程學(xué)院江蘇徐州6復(fù)旦大學(xué)附屬中山醫(yī)院耳鼻喉科上海000上海交通大學(xué)機(jī)械系統(tǒng)與振動(dòng)國家重點(diǎn)實(shí)驗(yàn)室上海0040
    噪聲與振動(dòng)控制 2016年2期
    關(guān)鍵詞:中耳耳蝸

    姚海峰,劉后廣,周 雷,饒柱石,黃新生(.中國礦業(yè)大學(xué)機(jī)電工程學(xué)院,江蘇徐州6;.復(fù)旦大學(xué)附屬中山醫(yī)院耳鼻喉科,上海000;.上海交通大學(xué)機(jī)械系統(tǒng)與振動(dòng)國家重點(diǎn)實(shí)驗(yàn)室,上海0040)

    ?

    人體整耳傳聲有限元數(shù)值模擬

    姚海峰1,劉后廣1,周雷2,饒柱石3,黃新生2
    (1.中國礦業(yè)大學(xué)機(jī)電工程學(xué)院,江蘇徐州221116;2.復(fù)旦大學(xué)附屬中山醫(yī)院耳鼻喉科,上海200032;3.上海交通大學(xué)機(jī)械系統(tǒng)與振動(dòng)國家重點(diǎn)實(shí)驗(yàn)室,上海200240)

    摘要:由于人耳的復(fù)雜幾何形態(tài)、超微結(jié)構(gòu)特性,很難通過實(shí)驗(yàn)測(cè)量其動(dòng)態(tài)響應(yīng),所以建立精確有限元模型對(duì)研究人耳聲音傳導(dǎo)特性有非常重要的作用。而現(xiàn)有人耳模型中耳蝸部分多沒有考慮其內(nèi)部?jī)汕惑w積不對(duì)稱。在原有所建人體外耳、中耳力學(xué)模型基礎(chǔ)上,參考耳蝸尺寸實(shí)驗(yàn)測(cè)量數(shù)據(jù),建立考慮耳蝸兩腔不對(duì)稱的整耳模型。并在外耳道施加90 dB的聲壓,利用有限元軟件對(duì)模型進(jìn)行諧響應(yīng)分析。最后,分別通過鐙骨底板位移頻響、鼓膜臍部位移頻響、鐙骨速度傳遞函數(shù)、耳蝸輸入阻抗的仿真結(jié)果數(shù)據(jù)與國外實(shí)驗(yàn)測(cè)量數(shù)據(jù)進(jìn)行對(duì)比,從而驗(yàn)證所建模型的可靠性。

    關(guān)鍵詞:振動(dòng)與波;中耳;耳蝸;前庭階;鼓階;有限元建模

    耳聾是社會(huì)上最常見的疾病之一,2006年全國第二次殘疾人抽樣調(diào)查報(bào)告顯示,國內(nèi)有2 780萬人患有聽力損傷[1]。這就迫切需要加強(qiáng)對(duì)人耳傳聲特性的研究,針對(duì)性地開發(fā)更有效的助聽裝置[2]。由于人耳內(nèi)部的中耳、耳蝸部分結(jié)構(gòu)復(fù)雜,體積微小,這給實(shí)驗(yàn)測(cè)量研究其傳聲特性增大了難度。針對(duì)該問題,國內(nèi)外很多學(xué)者通過力學(xué)模擬的方法,來輔助分析人耳傳聲特性,常見的有集總參數(shù)法[3]、多體動(dòng)力學(xué)法[4]、等效電路法[5]、有限元法[6]。其中,有限元法更易于模擬人耳這類具有復(fù)雜幾何形態(tài)、超微結(jié)構(gòu)特征、非同質(zhì)性的生物組織。近年來,隨著計(jì)算機(jī)運(yùn)算速度的提升,國內(nèi)外很多學(xué)者都開始這一方面的研究。如GAN等[7]通過對(duì)顳骨的組織進(jìn)行切片并結(jié)合CAD軟件,建立中耳有限元模型;劉迎曦等[8]及姚文娟等[9]都借助于更快捷的CT逆向成型法,建立人體整耳模型。王學(xué)林等[10]也通過該方法,進(jìn)一步建立考慮耳蝸主動(dòng)機(jī)制的整耳模型。與此同時(shí),饒柱石等[11]又利用精度更高的Micro-CT掃描,建立了人體整耳有限元模型。這些模型在人耳傳聲模擬研究中具有重要意義。然而,實(shí)驗(yàn)測(cè)量顯示[12],耳蝸內(nèi)的前庭階和鼓階是不對(duì)稱的兩腔,而現(xiàn)有整耳模型的耳蝸部分多做近似對(duì)稱處理。這種簡(jiǎn)化建模方法,會(huì)影響模擬耳蝸內(nèi)部傳聲的準(zhǔn)確性。針對(duì)上述問題,基于實(shí)驗(yàn)報(bào)道的耳蝸內(nèi)部結(jié)構(gòu)數(shù)據(jù),在原有人體外耳、中耳力學(xué)模型基礎(chǔ)上,建立考慮耳蝸前庭階、鼓階不對(duì)稱結(jié)構(gòu)的人體整耳有限元模型。最終,通過4組仿真數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)的對(duì)比,驗(yàn)證了模型的可靠性。

    1 有限元建模

    1.1模型結(jié)構(gòu)與建模方法

    建立的人體整耳力學(xué)模型主要包括耳道、中耳及耳蝸,如圖1所示。其中,中耳的組成結(jié)構(gòu)主要包括鼓膜、3塊聽小骨(錘骨、砧骨和鐙骨)以及在它們之間起連接作用的關(guān)節(jié)、韌帶和肌肉等。該中耳部分及耳道部分的具體建模過程和所取得材料參數(shù)同文獻(xiàn)[13]一致。

    圖1 人體整耳有限元模型

    整耳力學(xué)模型中的耳蝸部分簡(jiǎn)化為非螺旋的雙腔結(jié)構(gòu)。該耳蝸模型主要結(jié)構(gòu)包括:卵圓窗、圓窗、前庭階(SV)、鼓階(ST)、基底膜、蝸孔。建模過程中,首先利用三維建模軟件建立三維立體模型;再將三維模型保存為SAT格式導(dǎo)入至網(wǎng)格劃分軟件Hypermesh中,使用Hyermesh對(duì)耳蝸設(shè)定材料參數(shù)并且劃分六面體網(wǎng)格。最后將劃分網(wǎng)格后的有限元模型導(dǎo)入至有限元軟件Abaqus中,設(shè)定相應(yīng)邊界條件,并在模型外耳道口施加相應(yīng)的模擬聲壓激勵(lì)。

    1.2耳蝸結(jié)構(gòu)尺寸

    Wysocki利用25個(gè)人體顳骨作為樣本進(jìn)行了研究,得出前庭階與鼓階的精確尺寸數(shù)據(jù),包括其橫截面積與寬度[12]。此外Wever針對(duì)于前庭階橫截面也做了相關(guān)實(shí)驗(yàn),得出了精確的實(shí)驗(yàn)數(shù)據(jù)[14]。根據(jù)以上兩位學(xué)者的實(shí)驗(yàn)數(shù)據(jù)建立耳蝸三維模型,如圖2所示,圖中單位均為mm。前庭階與鼓階沿其長(zhǎng)度方向不斷變化的橫截面積大小如圖3所示。

    圖2 耳蝸有限元模型

    圖3 耳蝸兩腔的橫截面積

    同時(shí),模型中的耳蝸寬度也參考Wysocki的實(shí)驗(yàn)數(shù)據(jù)沿耳蝸縱向變化,如圖4所示。

    圖4 耳蝸寬度

    耳蝸整體長(zhǎng)度取36.5 mm。圓窗面積依據(jù)Atturo等實(shí)驗(yàn)數(shù)據(jù)取2.5 mm2[15],并依據(jù)文獻(xiàn)[16]將圓窗厚度近似取0.06 mm?;啄さ奈恢脧木鄨A窗2.5 mm開始,長(zhǎng)度為32 mm,厚度從基部(靠近鐙骨處)開始由7.5μm線性變化到頂部(靠近蝸孔處)的2.5μm,寬度從基部的100μm線性變化為頂部的500μm[17]。蝸孔的長(zhǎng)度為1mm[7]。

    1.3耳蝸材料參數(shù)

    模型中耳蝸部分材料屬性主要包括圓窗、前庭階、鼓階、基底膜的材料參數(shù)。其中前庭階與鼓階中充滿了淋巴液。圓窗密度取為1 200 kg/m3,阻尼系數(shù)β=5×10-5s,泊松比為0.3,楊氏模量為0.3 MPa?;啄ぶ伟宓拿芏葹?.2×103kg/m3,阻尼系數(shù)β= 5×10-5s,泊松比為0.3,楊氏模量為2.1×109MPa?;啄さ拿芏葹? 200 kg/m3,阻尼系數(shù)β=1×10-4s,泊松比為0.3,其楊氏模量由基部的40 MPa線性減小到頂部的3 MPa[11]。

    1.4耳蝸部分流固耦合數(shù)學(xué)模型

    建立的耳蝸模型中包括了兩腔中淋巴液分別與鐙骨底板、圓窗和基底膜間的流固耦合。假設(shè)兩腔中的淋巴液是非黏性液體,滿足式(1)的控制方程[18],即

    耳蝸模型中的淋巴液材料參數(shù)近似取與水的材料參數(shù)相同,體積模量為2 250 MPa,密度為1000kg/m3[19]。

    在有限元軟件Abaqus中定義流固耦合界面,將流體壓力耦合到結(jié)構(gòu)分析中,該模型中耦合面包括:鐙骨底板與前庭階中淋巴液耦合,基底膜兩側(cè)面分別與兩腔中淋巴液耦合,圓窗與鼓階中淋巴液耦合。

    2 模型可靠性驗(yàn)證

    為了驗(yàn)證建立的人體整耳有限元模型的可靠性,在外耳道入口處施加90 dB聲壓(0.632 N/m2面力),并將其計(jì)算結(jié)果與國外相關(guān)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比。

    2.1鐙骨底板位移頻響曲線

    GAN等人對(duì)10例人耳標(biāo)本進(jìn)行試驗(yàn)研究,得出鐙骨底板與鼓膜臍部在90 dB聲壓激勵(lì)下的位移平均值[20],在模型外耳道入口處施加90 dB聲壓后,進(jìn)行相應(yīng)的諧響應(yīng)仿真分析,從而得出其在250 Hz至8 000 Hz范圍內(nèi)對(duì)應(yīng)的鐙骨底板位移諧響應(yīng)曲線,并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比,如圖5所示。由圖可以看出,模型的仿真結(jié)果整體上與實(shí)驗(yàn)結(jié)果一致。高頻段小于實(shí)驗(yàn)測(cè)量值,出現(xiàn)這一現(xiàn)象是因?yàn)椴捎萌鹄枘釋?duì)中耳阻尼進(jìn)行近似模擬,該方法增大了人耳模型中的高頻阻尼值[21]。此外,在頻率為3 000 Hz左右時(shí)出現(xiàn)了峰值現(xiàn)象,這是因?yàn)橥舛涝诖祟l率段對(duì)聲音有共振放大作用[22];而GAN等試驗(yàn)測(cè)量時(shí),聲音直接施加在靠近鼓膜處,沒有在外耳道口處施加,故所測(cè)結(jié)果不會(huì)出現(xiàn)這種耳道放大的效果。

    2.2鼓膜臍部位移頻響曲線

    將模型在90 dB聲壓激勵(lì)下鼓膜臍部位移頻響曲線與GAN等人的平均實(shí)驗(yàn)數(shù)據(jù)曲線進(jìn)行對(duì)比[20],結(jié)果如圖6所示。模型計(jì)算結(jié)果在趨勢(shì)、幅值上和實(shí)驗(yàn)數(shù)據(jù)都較一致。此外,同鐙骨對(duì)比一樣,該對(duì)比結(jié)果也出現(xiàn)3 000 Hz處放大及高頻段小于實(shí)驗(yàn)值的現(xiàn)象。同理,這些分別由耳道放大及阻尼模擬方法所致。

    圖5 鐙骨底板位移曲線仿真值與實(shí)驗(yàn)值對(duì)比

    圖6 鼓膜臍部位移曲線仿真值與實(shí)驗(yàn)值對(duì)比

    2.3鐙骨底板速度傳遞函數(shù)

    Aibara等多位學(xué)者于2001年運(yùn)用多普勒測(cè)振儀以12例顳骨為樣本進(jìn)行實(shí)驗(yàn)測(cè)量,最終得出鐙骨速度傳遞函數(shù)(STF)曲線[23],用以表示中耳聲傳遞特性。其計(jì)算公式如式(2)所示,式中V為鐙骨底板速度,PTM為鼓膜臍部聲壓值。在外耳道開口處施加90 dB聲壓,對(duì)整耳模型進(jìn)行仿真分析,并得出仿真數(shù)據(jù),結(jié)果如圖7。由圖可見,模型仿真結(jié)果與12例顳骨樣本實(shí)驗(yàn)數(shù)據(jù)中幅值最低的一例曲線比較接近,整體趨勢(shì)相一致。

    圖7 鐙骨底板速度傳遞函數(shù)仿真值與實(shí)驗(yàn)值對(duì)比

    2.4耳蝸輸入阻抗

    耳蝸的輸入阻抗是中耳將聲音能量傳遞給耳蝸的一個(gè)重要的參數(shù)變量,可以將它定義為前庭階內(nèi)靠近卵圓窗處聲壓與靠近卵圓窗淋巴液體積速度之比。如圖8所示,將所建模型仿真分析得出的耳蝸輸入阻抗值與Aibara和Puria的研究數(shù)據(jù)進(jìn)行對(duì)比[23, 24]。圖中Mean與Lower Bound曲線分別表示平均值與下限值。模型的仿真分析數(shù)據(jù)與Puria的實(shí)驗(yàn)數(shù)據(jù)大小比較接近,并且整體趨勢(shì)相一致。仿真數(shù)據(jù)變化范圍基本在實(shí)驗(yàn)數(shù)據(jù)的變化范圍之內(nèi)。

    圖8 耳蝸輸入阻抗曲線仿真值與實(shí)驗(yàn)值對(duì)比

    2.5模型可靠性驗(yàn)證結(jié)果分析

    通過以上四組模型仿真數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)對(duì)比分析可知,模型的仿真結(jié)果數(shù)值在整體趨勢(shì)上與實(shí)驗(yàn)測(cè)量數(shù)據(jù)值保持一致,可以滿足對(duì)于人體整耳實(shí)際結(jié)構(gòu)聲學(xué)特性的研究要求。同時(shí)由于建立模型時(shí)所取的個(gè)體樣本差異,以及在建模時(shí)部分采取簡(jiǎn)單近似處理等原因,最終導(dǎo)致仿真結(jié)果與實(shí)驗(yàn)均值存在一些偏差[25],但這并不影響模型整體的有效性。

    3 結(jié)語

    參考耳蝸幾何尺寸實(shí)驗(yàn)測(cè)量數(shù)據(jù),基于前期建立的外耳、中耳模型,建立了人體整耳力學(xué)模型。為了對(duì)模型的可靠性進(jìn)行驗(yàn)證,針對(duì)其四組仿真數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析比對(duì),即鐙骨底板位移頻響、鼓膜臍部位移頻響、鐙骨底板速度傳遞函數(shù)、耳蝸輸入阻抗。根據(jù)分析比對(duì)的結(jié)果顯示,四組仿真數(shù)據(jù)與對(duì)應(yīng)的實(shí)驗(yàn)數(shù)據(jù)在趨勢(shì)與幅值上均比較接近,滿足人耳傳聲特性的模擬要求。后期將會(huì)進(jìn)一步對(duì)該模型進(jìn)行完善改進(jìn),如對(duì)韌帶等組織進(jìn)行粘彈性模擬、考慮基底膜的各向異性及鼓膜的多層結(jié)構(gòu)等,更深入地研究整耳的結(jié)構(gòu)聲學(xué)特性。

    參考文獻(xiàn):

    [1]第二次全國殘疾人抽樣調(diào)查辦公室.第二次全國殘疾人抽樣調(diào)查資料[M].北京:中國統(tǒng)計(jì)出版社,2007.

    [2]劉后廣,閔小峰,塔娜,等.基于耦合模型的人工中耳壓電振子設(shè)計(jì)[J].噪聲與振動(dòng)控制,2010,30(2):130-133.

    [3] Feng B, GAN R Z. Lumped parametric model of the human ear for sound transmission[J]. Biomechanics and Modeling in Mechanobiology, 2004, 3(1):33-47.

    [4] Volandri G, Puccio F Di, Forte P, et al. Model-oriented review and multi-body simulation of theossicular chain of the human middle ear[J]. Medical Engineering & Physics, 2012, 34(9):1339-1355.

    [5] Kevin N O’Connor, Sunil Puria. Middle-ear circuit model parameters based on a population of human ears[J]. Acoustical Society of America, 2008, 123(1):197-211.

    [6] Zhao Fei, Koike Takuji, Wang Jie, et al. Finite element analysis of the middle ear transfer functions and related pathologies[J]. Medical Engineering & Physics, 2009, 31 (8):907-916.

    [7] GAN R Z, REEVES B P, Wang X L. Modeling of sound transmission from ear canal to cochlea[J]. Ann Biomed Eng, 2007, 35(12):2108-2195.

    [8]劉迎曦,李生,孫秀珍.人耳傳聲數(shù)值模型[J].力學(xué)學(xué)報(bào),2008,40(1):107-113.

    [9]姚文娟,李兵,胡寶琳,等.置換部分聽骨贗復(fù)物后對(duì)人耳聽力恢復(fù)的影響[J].醫(yī)用生物力學(xué),2012,27(1):58-64.

    [10]王學(xué)林,周建軍,凌玲,等.含主動(dòng)耳蝸的人耳傳聲有限元模擬[J].振動(dòng)與沖擊,2012,31(21):41-45.

    [11] Tian J B, Huang X S, Rao Z S, et al. Finite element analysis of the effect of actuator coupling conditions on round window stimulation[J]. Journal of Mechanics in Medicine and Biology, 2015, 15(4):1550048(19pages).

    [12] Wysocki J. Dimensions of the humen vestibular and tympanic scalae[J]. Hearing Research, 1999, 135(1- 2):39-46.

    [13] Zhou L, Feng M L, Wang W, et al. Study on the role of ossicular joint using finite element method[J]. Journal of Mechanics in Medicine and Biology, 2015, 16(2):1650041(20pages).

    [14] Wever E G. Theory of hearing[M]. London, 1949:437-444.

    [15] Atturo F, Barbara M, Rask- Andersen H. Is the human round window really round? an anatomic study with surgical implications[J]. Otology & Neurotology, 2014, 35 (8):1354-1360.

    [16] Kwacz M, Marek P, Borkowski P, et al. A threedimensional finite element model of round window membrane vibration before and after stapedotomy surgery [J]. Biomech Model Mechanobiol, 2013, 12(6):1243-1261.

    [17] Wang X L, Wang L L, Zhou J J, et al. Finite element modeling of human auditory periphery including a feedforward amplification of the cochlea[J], Computer Methods in Biomechanics and Biomedical Engineering, 2014, 17(10):1096-107.

    [18] Tian JB, Ta N, Rao Z S, et al. Finiteelement modeling of sound transmission based on micro- computer tomography for human ear[J]. Applied Mechanics and Materials, 2013, 419:593-601.

    [19] Kwacz M, Marek P, Borkowski P, et al. Effect of different stapes prostheses on the passive vibration of the basilar membrane[J]. Hear Research, 2014, 310(2):13-26.

    [20] GAN R Z, Wood M W, Dormer K J. Human middle ear transfer function measured by double laser interferometry system[J]. Otology & Neurotology, 2004, 25(4):423-435.

    [21] Zhang X, GAN R Z. Finite element modeling of energy absorbance in normal and disordered human ears[J]. Hear Research, 2013, 301(7):146-155.

    [22] Ferris P, Prendergast PJ. Middle-ear dynamicsbeforeand after ossicular replacement[J]. Journal of Biomechanics, 2000, 33(5):581-590.

    [23] Aibara R, Welsh J T, Puria S. Human middle-ear sound transfer function and cochlear input impedance[J]. Hearing Research, 2001,152(1-2):100-109.

    [24] Puria S, Peake W T, Rosowski J J. Sound- pressure measurements in the cochlear vestibule of human-cadaver ears[J].Acoust Soc Am, 1997, 101(5):2754-2770.

    [25]劉后廣,塔娜,饒柱石.人體中耳有限元法數(shù)值仿真[J].系統(tǒng)仿真學(xué)報(bào),2009,21(24):7899-7901.

    Finite Element Simulation of Human Ears

    YAO Hai-feng1, LIU Hou-guang1, ZHOU Lei2, RAO Zhu-shi3, HUANG Xin-sheng2
    ( 1. School of Mechatronic Engineering, ChinaUniversity of Miningand Technology, Xuzhou 221116, Jiangsu China; 2. Department of Otorhinolaryngology, Shanghai Zhongshan Hospital Affiliatedto Fudan University, Shanghai 200032, China; 3. StateKey Laboratory of Mechanical Systemand Vibration, Shanghai Jiaotong University, Shanghai 200240, China)

    Abstract:Due to the complexity and tiny structure of human ears, it is difficult to carry out the experiment for measuring their dynamic response. In thispaper, ahuman ear finiteelement model, which considersthetwo cavity volume’s asymmetry of its cochlea, was established. Then, a sound pressure of 90 dB was applied to the entrance of the external ear canal, and a harmonic analysis was carried out. Finally, the validity of this model was confirmed by comparing the modelpredicted results, which are the cochlear input impedance and the motions of the tympanic membrane and stapes, with published measurement data for human temporal bones. The result shows that the model is reasonable in predicting the biomechanicsbehavior of humanears.

    Key words:vibrationandwave; middleear; cochlea; scalavestibule; scalatympani; finiteelement modeling

    通訊作者:劉后廣,男,博士生導(dǎo)師,副教授。Email:liuhg@cumt.edu.cn

    作者簡(jiǎn)介:姚海峰(1990- ),男,安徽省合肥市人,碩士研究生,主要研究方向?yàn)槲?chuàng)醫(yī)療技術(shù)與器械研究。

    基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(51305442);江蘇省自然科學(xué)基金資助項(xiàng)目(BK20130194);博士后特別資助項(xiàng)目(2015T80597)

    收稿日期:2015-10-27

    文章編號(hào):1006-1355(2016)02-0061-04+83

    中圖分類號(hào):Q62

    文獻(xiàn)標(biāo)識(shí)碼:ADOI編碼:10.3969/j.issn.1006-1335.2016.02.013

    猜你喜歡
    中耳耳蝸
    耳蝸微音器電位臨床操作要點(diǎn)
    開放式鼓室成形術(shù)治療中耳膽脂瘤的效果及鼓膜愈合率觀察
    DR內(nèi)聽道像及多層螺旋CT三維重建對(duì)人工耳蝸的效果評(píng)估
    中耳膽脂瘤并周圍性面癱的臨床分析
    豚鼠耳蝸Hensen細(xì)胞脂滴的性質(zhì)與分布
    基于Gammachirp耳蝸能量譜特征提取的音頻指紋算法
    潛水員中耳氣壓傷的發(fā)病特點(diǎn)與防治
    先天性中耳畸形的臨床分型及其與耳聾的相關(guān)性研究
    內(nèi)聽道狹窄合并耳蝸神經(jīng)管狹窄的影像學(xué)表現(xiàn)
    小耳朵怎么會(huì)有水
    中文字幕精品亚洲无线码一区| 51午夜福利影视在线观看| 中文字幕熟女人妻在线| 成人特级av手机在线观看| 国产蜜桃级精品一区二区三区| 国内揄拍国产精品人妻在线| 国产精品久久久久久亚洲av鲁大| 亚洲电影在线观看av| 亚洲一区二区三区不卡视频| or卡值多少钱| 免费观看人在逋| 天美传媒精品一区二区| 欧美zozozo另类| bbb黄色大片| 日本三级黄在线观看| 在线观看午夜福利视频| 老鸭窝网址在线观看| 欧美zozozo另类| 男女视频在线观看网站免费| 亚洲av熟女| 香蕉丝袜av| 女生性感内裤真人,穿戴方法视频| 国产成年人精品一区二区| 日本精品一区二区三区蜜桃| 亚洲精品成人久久久久久| 日本与韩国留学比较| 亚洲国产精品sss在线观看| 日韩欧美国产一区二区入口| 久久久久久人人人人人| 熟女电影av网| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久久中文| 国产精品野战在线观看| 人人妻人人看人人澡| 国产激情欧美一区二区| 免费在线观看日本一区| 看免费av毛片| 久久久精品大字幕| 国产伦一二天堂av在线观看| 最新在线观看一区二区三区| 又爽又黄无遮挡网站| 精品无人区乱码1区二区| 日韩欧美免费精品| 琪琪午夜伦伦电影理论片6080| 午夜福利免费观看在线| or卡值多少钱| 国产精品 国内视频| 有码 亚洲区| 少妇丰满av| 草草在线视频免费看| 亚洲熟妇中文字幕五十中出| 国产真实伦视频高清在线观看 | 美女 人体艺术 gogo| 久99久视频精品免费| 操出白浆在线播放| 久久精品国产99精品国产亚洲性色| 18禁在线播放成人免费| 国产三级黄色录像| 人人妻,人人澡人人爽秒播| 久久婷婷人人爽人人干人人爱| 国产精品久久久久久久电影 | 男女视频在线观看网站免费| 一夜夜www| 欧美性猛交╳xxx乱大交人| 男女床上黄色一级片免费看| 日本精品一区二区三区蜜桃| 久久精品人妻少妇| 成人三级黄色视频| 内射极品少妇av片p| 脱女人内裤的视频| 国产精品嫩草影院av在线观看 | 国产亚洲精品久久久久久毛片| 国产乱人伦免费视频| 中文字幕高清在线视频| 少妇的丰满在线观看| 亚洲黑人精品在线| 天堂动漫精品| 亚洲激情在线av| 国产精品久久久久久久电影 | 成人三级黄色视频| 国产精品香港三级国产av潘金莲| 3wmmmm亚洲av在线观看| www.999成人在线观看| 国产亚洲欧美在线一区二区| 日本免费a在线| 亚洲av第一区精品v没综合| 在线观看午夜福利视频| 19禁男女啪啪无遮挡网站| 亚洲七黄色美女视频| 成年人黄色毛片网站| 首页视频小说图片口味搜索| 男女那种视频在线观看| 欧美+日韩+精品| 亚洲精品在线美女| 国产真实伦视频高清在线观看 | 一级黄色大片毛片| 男人舔奶头视频| 丁香欧美五月| or卡值多少钱| 欧美成人免费av一区二区三区| 白带黄色成豆腐渣| 高清日韩中文字幕在线| 男插女下体视频免费在线播放| 久久久精品大字幕| 婷婷精品国产亚洲av在线| 俄罗斯特黄特色一大片| 校园春色视频在线观看| 午夜两性在线视频| 亚洲内射少妇av| bbb黄色大片| 欧美日韩国产亚洲二区| 一区福利在线观看| 给我免费播放毛片高清在线观看| 国产在线精品亚洲第一网站| 三级国产精品欧美在线观看| 哪里可以看免费的av片| 热99re8久久精品国产| 亚洲成av人片在线播放无| 少妇人妻一区二区三区视频| 嫩草影院入口| xxx96com| 日韩精品青青久久久久久| 深爱激情五月婷婷| 亚洲男人的天堂狠狠| 变态另类成人亚洲欧美熟女| 亚洲国产欧洲综合997久久,| 日本撒尿小便嘘嘘汇集6| 99久久99久久久精品蜜桃| 国产野战对白在线观看| 一个人观看的视频www高清免费观看| 色尼玛亚洲综合影院| 国产极品精品免费视频能看的| 九色成人免费人妻av| 99热6这里只有精品| 九色国产91popny在线| 亚洲午夜理论影院| 精品不卡国产一区二区三区| 亚洲在线自拍视频| 精品人妻一区二区三区麻豆 | 女同久久另类99精品国产91| av国产免费在线观看| 舔av片在线| 在线观看日韩欧美| 久久久精品大字幕| 国产亚洲精品久久久久久毛片| 日本 av在线| 51午夜福利影视在线观看| 亚洲人与动物交配视频| 成人亚洲精品av一区二区| 免费在线观看日本一区| av福利片在线观看| 午夜福利18| 国产一区二区在线av高清观看| 成人无遮挡网站| 亚洲精品456在线播放app | 天天躁日日操中文字幕| 国产单亲对白刺激| 小蜜桃在线观看免费完整版高清| 999久久久精品免费观看国产| 一区福利在线观看| 搡老妇女老女人老熟妇| 日本在线视频免费播放| 麻豆久久精品国产亚洲av| 乱人视频在线观看| 色尼玛亚洲综合影院| 身体一侧抽搐| 18禁美女被吸乳视频| 法律面前人人平等表现在哪些方面| 国产亚洲精品久久久久久毛片| 成年女人永久免费观看视频| 精品99又大又爽又粗少妇毛片 | 少妇高潮的动态图| 长腿黑丝高跟| 91麻豆av在线| 精品一区二区三区人妻视频| 午夜亚洲福利在线播放| 久久久久久久亚洲中文字幕 | 国产主播在线观看一区二区| 久久久久久久久大av| 香蕉丝袜av| 亚洲国产精品成人综合色| 日韩欧美免费精品| 国产精品三级大全| 高清日韩中文字幕在线| 久久久精品欧美日韩精品| 亚洲av美国av| 精华霜和精华液先用哪个| 婷婷精品国产亚洲av在线| 搡女人真爽免费视频火全软件 | 亚洲av免费在线观看| 亚洲av熟女| 91麻豆av在线| 成人永久免费在线观看视频| 99国产极品粉嫩在线观看| 国产精品三级大全| 久9热在线精品视频| 欧美中文综合在线视频| 无限看片的www在线观看| 此物有八面人人有两片| 哪里可以看免费的av片| 免费观看的影片在线观看| 久久草成人影院| 欧美日韩国产亚洲二区| 精品人妻1区二区| 欧洲精品卡2卡3卡4卡5卡区| 日韩精品中文字幕看吧| 黄色女人牲交| 给我免费播放毛片高清在线观看| 成人欧美大片| 久久欧美精品欧美久久欧美| 极品教师在线免费播放| 欧美黄色淫秽网站| 又黄又爽又免费观看的视频| 久久久久久久久大av| 久9热在线精品视频| 亚洲国产欧美网| 国产真实伦视频高清在线观看 | 老汉色∧v一级毛片| 一区福利在线观看| 99国产综合亚洲精品| 午夜亚洲福利在线播放| 青草久久国产| 十八禁人妻一区二区| 好男人在线观看高清免费视频| 久久久久九九精品影院| 亚洲人成伊人成综合网2020| tocl精华| 久久久久国内视频| 成年女人毛片免费观看观看9| 亚洲第一电影网av| 久久久久久久亚洲中文字幕 | 欧美日韩中文字幕国产精品一区二区三区| 国产精华一区二区三区| 两个人看的免费小视频| 在线天堂最新版资源| 亚洲精品美女久久久久99蜜臀| 成年人黄色毛片网站| 精品久久久久久成人av| 国产免费男女视频| 成人三级黄色视频| 69人妻影院| 久久草成人影院| 成人高潮视频无遮挡免费网站| 丁香欧美五月| 日本a在线网址| 精品日产1卡2卡| 岛国在线观看网站| 搡女人真爽免费视频火全软件 | 最后的刺客免费高清国语| 国产aⅴ精品一区二区三区波| 欧美又色又爽又黄视频| a在线观看视频网站| 欧美日韩国产亚洲二区| 国产三级黄色录像| 日本一本二区三区精品| 草草在线视频免费看| xxxwww97欧美| 12—13女人毛片做爰片一| 91久久精品国产一区二区成人 | 国产精品久久久人人做人人爽| e午夜精品久久久久久久| 亚洲成人免费电影在线观看| 啪啪无遮挡十八禁网站| 内射极品少妇av片p| 麻豆国产97在线/欧美| 法律面前人人平等表现在哪些方面| 神马国产精品三级电影在线观看| 禁无遮挡网站| 免费av不卡在线播放| 女生性感内裤真人,穿戴方法视频| 色精品久久人妻99蜜桃| 久久精品国产清高在天天线| 狠狠狠狠99中文字幕| 一边摸一边抽搐一进一小说| 精品熟女少妇八av免费久了| 在线看三级毛片| 校园春色视频在线观看| 99精品在免费线老司机午夜| 在线观看66精品国产| 91久久精品电影网| 亚洲成人精品中文字幕电影| 婷婷六月久久综合丁香| av在线天堂中文字幕| 亚洲人与动物交配视频| 色综合婷婷激情| 男女午夜视频在线观看| 男女之事视频高清在线观看| 精品99又大又爽又粗少妇毛片 | 免费看十八禁软件| 观看免费一级毛片| 国产高潮美女av| 最近视频中文字幕2019在线8| 免费看日本二区| av福利片在线观看| 亚洲精品美女久久久久99蜜臀| bbb黄色大片| 一个人看视频在线观看www免费 | 久久久久精品国产欧美久久久| 身体一侧抽搐| 国产精品精品国产色婷婷| 国产久久久一区二区三区| 三级国产精品欧美在线观看| 国产精品精品国产色婷婷| 老司机福利观看| 网址你懂的国产日韩在线| 丰满人妻一区二区三区视频av | 一级毛片高清免费大全| 村上凉子中文字幕在线| 啪啪无遮挡十八禁网站| 日本一本二区三区精品| 97碰自拍视频| 亚洲aⅴ乱码一区二区在线播放| 亚洲天堂国产精品一区在线| 老熟妇仑乱视频hdxx| 免费大片18禁| 内地一区二区视频在线| 日本成人三级电影网站| 亚洲av日韩精品久久久久久密| 亚洲久久久久久中文字幕| 99视频精品全部免费 在线| 国产高清视频在线观看网站| 岛国在线免费视频观看| 色av中文字幕| 国产亚洲精品综合一区在线观看| 亚洲精品成人久久久久久| 久久九九热精品免费| 国产亚洲精品av在线| 国产欧美日韩一区二区三| 三级毛片av免费| 男人的好看免费观看在线视频| 99精品在免费线老司机午夜| 中文字幕人成人乱码亚洲影| 久久九九热精品免费| 熟妇人妻久久中文字幕3abv| 精品乱码久久久久久99久播| 亚洲成a人片在线一区二区| 国产精品国产高清国产av| 男人舔女人下体高潮全视频| 十八禁人妻一区二区| 国产激情偷乱视频一区二区| 亚洲自拍偷在线| 51午夜福利影视在线观看| 国产精品永久免费网站| 成人一区二区视频在线观看| 99国产精品一区二区三区| 久久久久国产精品人妻aⅴ院| 日韩欧美在线乱码| 在线a可以看的网站| 日本与韩国留学比较| 国产一区二区三区视频了| 国产av在哪里看| 香蕉av资源在线| 超碰av人人做人人爽久久 | 精品久久久久久久人妻蜜臀av| 97人妻精品一区二区三区麻豆| 欧美xxxx黑人xx丫x性爽| av专区在线播放| 欧美av亚洲av综合av国产av| 久久天躁狠狠躁夜夜2o2o| 国产激情欧美一区二区| 亚洲精品成人久久久久久| 欧美乱妇无乱码| 无限看片的www在线观看| 国产精品一及| 久久久久性生活片| 日本成人三级电影网站| 亚洲国产欧美人成| 欧美另类亚洲清纯唯美| 日韩欧美 国产精品| 69av精品久久久久久| av专区在线播放| 国产免费av片在线观看野外av| 1024手机看黄色片| 无人区码免费观看不卡| 国产中年淑女户外野战色| 日韩中文字幕欧美一区二区| 啦啦啦免费观看视频1| 三级毛片av免费| 亚洲中文字幕一区二区三区有码在线看| or卡值多少钱| 毛片女人毛片| 一个人看视频在线观看www免费 | 国产麻豆成人av免费视频| 亚洲精品一卡2卡三卡4卡5卡| 99热只有精品国产| tocl精华| 全区人妻精品视频| 日本a在线网址| xxx96com| 成人三级黄色视频| 青草久久国产| 日本熟妇午夜| av欧美777| 精品一区二区三区视频在线 | 国产精品爽爽va在线观看网站| 欧美中文综合在线视频| 熟女电影av网| 母亲3免费完整高清在线观看| 午夜日韩欧美国产| 亚洲av美国av| 色av中文字幕| 女人高潮潮喷娇喘18禁视频| 亚洲成人久久爱视频| 国产精品亚洲一级av第二区| 搡老妇女老女人老熟妇| 看片在线看免费视频| 久久精品国产清高在天天线| 欧美中文日本在线观看视频| 午夜福利高清视频| 国产精品久久久久久人妻精品电影| 亚洲欧美日韩高清在线视频| 日韩欧美精品v在线| 一个人观看的视频www高清免费观看| 国产精品野战在线观看| 国产高清激情床上av| 国产日本99.免费观看| 亚洲精品在线美女| 欧美在线一区亚洲| 人人妻人人看人人澡| 黄片小视频在线播放| 首页视频小说图片口味搜索| 亚洲成人久久性| 日本 欧美在线| 无限看片的www在线观看| 亚洲内射少妇av| 黄色丝袜av网址大全| 美女大奶头视频| 女警被强在线播放| 国产精品久久久久久久电影 | 又黄又粗又硬又大视频| 久久精品影院6| 黄色丝袜av网址大全| 91av网一区二区| 狂野欧美白嫩少妇大欣赏| 日韩中文字幕欧美一区二区| 啪啪无遮挡十八禁网站| 久久精品亚洲精品国产色婷小说| 99在线人妻在线中文字幕| 亚洲精品456在线播放app | 女警被强在线播放| 熟妇人妻久久中文字幕3abv| 国产高清视频在线观看网站| 午夜福利成人在线免费观看| 久久久久亚洲av毛片大全| 亚洲最大成人手机在线| h日本视频在线播放| 亚洲七黄色美女视频| 欧美日本视频| 有码 亚洲区| 夜夜躁狠狠躁天天躁| 亚洲精品日韩av片在线观看 | 亚洲性夜色夜夜综合| 尤物成人国产欧美一区二区三区| 亚洲av成人不卡在线观看播放网| 免费人成视频x8x8入口观看| 国产伦在线观看视频一区| 国产三级黄色录像| 可以在线观看毛片的网站| 深夜精品福利| 成人18禁在线播放| 一区二区三区高清视频在线| 亚洲一区二区三区色噜噜| 18禁美女被吸乳视频| 午夜视频国产福利| 久久精品国产亚洲av香蕉五月| 搞女人的毛片| 日韩精品青青久久久久久| 色播亚洲综合网| 欧美最黄视频在线播放免费| 国产真实伦视频高清在线观看 | 在线观看66精品国产| 久久久精品欧美日韩精品| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 免费人成在线观看视频色| 精品国产美女av久久久久小说| 狂野欧美激情性xxxx| 男人舔奶头视频| av女优亚洲男人天堂| 全区人妻精品视频| 午夜福利在线观看吧| www日本在线高清视频| 在线观看美女被高潮喷水网站 | 色av中文字幕| 亚洲国产精品999在线| 在线看三级毛片| 97碰自拍视频| 久99久视频精品免费| 精品国产超薄肉色丝袜足j| 日本 欧美在线| 国产精品99久久99久久久不卡| 69av精品久久久久久| 在线十欧美十亚洲十日本专区| 极品教师在线免费播放| 天堂av国产一区二区熟女人妻| 午夜免费观看网址| 亚洲男人的天堂狠狠| 国产精品久久久久久亚洲av鲁大| 成人无遮挡网站| 日韩精品中文字幕看吧| 女人十人毛片免费观看3o分钟| 伊人久久精品亚洲午夜| 一夜夜www| 亚洲国产精品合色在线| 丝袜美腿在线中文| 亚洲色图av天堂| 色播亚洲综合网| 免费av不卡在线播放| 国模一区二区三区四区视频| av福利片在线观看| www.www免费av| 国产亚洲精品久久久久久毛片| 俺也久久电影网| 精品一区二区三区视频在线观看免费| 久久性视频一级片| 色哟哟哟哟哟哟| 欧美日韩福利视频一区二区| 好男人电影高清在线观看| 亚洲人成电影免费在线| 国产精品久久久久久久电影 | 欧美绝顶高潮抽搐喷水| 黄色视频,在线免费观看| netflix在线观看网站| 一个人免费在线观看的高清视频| 啦啦啦观看免费观看视频高清| 亚洲国产欧洲综合997久久,| 最新美女视频免费是黄的| 国产精品一及| 最新美女视频免费是黄的| 国产91精品成人一区二区三区| 在线播放国产精品三级| 国产精品久久久人人做人人爽| 精品久久久久久久久久久久久| 欧美午夜高清在线| 制服人妻中文乱码| 亚洲av美国av| 久久久精品大字幕| АⅤ资源中文在线天堂| 久久久精品大字幕| 少妇的逼好多水| 一区二区三区国产精品乱码| 国产精品爽爽va在线观看网站| 在线播放无遮挡| 欧美性猛交╳xxx乱大交人| 亚洲av五月六月丁香网| av女优亚洲男人天堂| 午夜福利欧美成人| 99精品在免费线老司机午夜| 91久久精品国产一区二区成人 | 免费看十八禁软件| 日韩欧美国产一区二区入口| av黄色大香蕉| 深爱激情五月婷婷| 亚洲精品在线美女| 国产v大片淫在线免费观看| 亚洲精品色激情综合| 熟妇人妻久久中文字幕3abv| 亚洲成av人片免费观看| 18禁黄网站禁片免费观看直播| 免费观看精品视频网站| 欧美在线一区亚洲| 成人国产一区最新在线观看| 黄色丝袜av网址大全| 亚洲性夜色夜夜综合| 成人特级黄色片久久久久久久| 女人被狂操c到高潮| 亚洲国产欧洲综合997久久,| 日韩精品中文字幕看吧| 一进一出抽搐动态| 久久99热这里只有精品18| 亚洲成av人片在线播放无| 一本久久中文字幕| 亚洲精品日韩av片在线观看 | 香蕉av资源在线| 中文资源天堂在线| 欧美乱妇无乱码| 啦啦啦免费观看视频1| 国产欧美日韩精品一区二区| 最近最新免费中文字幕在线| 国产av在哪里看| 99热这里只有是精品50| 床上黄色一级片| 亚洲国产日韩欧美精品在线观看 | 亚洲精品一区av在线观看| 在线免费观看不下载黄p国产 | 久久国产精品人妻蜜桃| 有码 亚洲区| 又黄又粗又硬又大视频| 深爱激情五月婷婷| 国产成人欧美在线观看| av国产免费在线观看| 国产成人啪精品午夜网站| 美女黄网站色视频| 精品久久久久久久人妻蜜臀av| 欧美性感艳星| 国产精品国产高清国产av| 十八禁人妻一区二区| 国产爱豆传媒在线观看| 性欧美人与动物交配| 一级a爱片免费观看的视频| 国产在视频线在精品| 亚洲,欧美精品.| av黄色大香蕉| 丰满的人妻完整版| 可以在线观看毛片的网站| 制服丝袜大香蕉在线| 亚洲精品一卡2卡三卡4卡5卡| 男人和女人高潮做爰伦理| 国产精品影院久久| 最新美女视频免费是黄的| 亚洲七黄色美女视频| 嫩草影视91久久| 午夜久久久久精精品| 亚洲av第一区精品v没综合| 成人国产一区最新在线观看| 男女那种视频在线观看|