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

    基于自由場(chǎng)爆炸的豬鼓膜破裂規(guī)律實(shí)驗(yàn)研究

    2024-12-13 00:00:00向書(shū)毅薛松波杜智博趙楊王興皓田旭高志強(qiáng)馮國(guó)棟費(fèi)舟莊茁柳占立
    爆炸與沖擊 2024年12期
    關(guān)鍵詞:沖擊波

    摘要: 聽(tīng)覺(jué)系統(tǒng)各組成部分的機(jī)械損傷是爆炸后造成聽(tīng)力損失的主要原因,強(qiáng)脈沖聲致聽(tīng)覺(jué)損害風(fēng)險(xiǎn)準(zhǔn)則仍然存在許多爭(zhēng)議,例如:指標(biāo)選擇沖量還是超壓峰值,正壓持續(xù)時(shí)間是否重要等。本研究基于自由場(chǎng)實(shí)爆條件,設(shè)計(jì)并搭建了大動(dòng)物爆炸致傷平臺(tái),探究了不同爆炸參數(shù)對(duì)鼓膜破裂的影響規(guī)律,并建立了基于自由場(chǎng)超壓峰值和正壓持續(xù)時(shí)間的鼓膜創(chuàng)傷量效關(guān)系。通過(guò)筆形壓力傳感器測(cè)量自由場(chǎng)超壓,通過(guò)Friedlander 公式擬合超壓時(shí)程曲線,確定沖擊波超壓峰值和正壓持續(xù)時(shí)間,并對(duì)時(shí)域中記錄的波形進(jìn)行歸一化能量頻譜分析,以確定沖擊波在頻域上的信號(hào)能量分布。對(duì)爆炸后的小型豬進(jìn)行解剖,記錄不同爆炸參數(shù)下鼓膜創(chuàng)傷程度。以超壓峰值和正壓持續(xù)時(shí)間為自變量,對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行二元邏輯回歸分析,并給出鼓膜破裂風(fēng)險(xiǎn)曲線。研究發(fā)現(xiàn),當(dāng)自由場(chǎng)超壓峰值低于170 kPa 時(shí),鼓膜無(wú)明顯損傷;當(dāng)自由場(chǎng)超壓峰值高于237 kPa 時(shí),部分鼓膜出現(xiàn)不同程度的破裂和充血。距爆心越近,超壓峰值越大,但鼓膜創(chuàng)傷的嚴(yán)重程度并未隨之單調(diào)增加。在8.0 kg TNT 當(dāng)量的爆炸實(shí)驗(yàn)中,鼓膜破裂的嚴(yán)重程度隨爆心距的減小呈現(xiàn)先提高再降低的趨勢(shì)。通過(guò)對(duì)沖擊波載荷特征的分析可知,距爆心越近,正壓持續(xù)時(shí)間越短,高頻段能量占比相對(duì)更大,小型豬鼓膜破裂的概率可能反而降低,此時(shí)仍然出現(xiàn)顯著的聽(tīng)力損失和耳蝸損傷。鼓膜作為通過(guò)振動(dòng)傳遞聲信號(hào)的黏彈性薄膜結(jié)構(gòu),其動(dòng)力學(xué)響應(yīng)可能與載荷頻率成分密切相關(guān)。除了超壓峰值,沖擊波波形頻譜分布對(duì)鼓膜破裂程度影響顯著。

    關(guān)鍵詞: 自由場(chǎng)爆炸;沖擊波;鼓膜破裂;聽(tīng)力損失;致傷風(fēng)險(xiǎn)

    中圖分類號(hào): O389 國(guó)標(biāo)學(xué)科代碼: 13035 文獻(xiàn)標(biāo)志碼: A

    聽(tīng)覺(jué)器官是爆炸沖擊波最敏感的靶器官之一,其爆炸沖擊傷發(fā)生率高,影響專業(yè)人員的作業(yè)能力,開(kāi)展爆炸沖擊波致聽(tīng)覺(jué)功能損害評(píng)估具有重要意義,目前爆炸致聽(tīng)力損失的形成和預(yù)防機(jī)制尚不清晰。鼓膜(tympanic membrane, TM)破裂是最常見(jiàn)的耳部爆炸傷,Gan 等[1] 通過(guò)激波管動(dòng)物實(shí)驗(yàn)測(cè)量研究了鼓膜破裂閾值和超壓波形之間的關(guān)系。Gan 等[2] 通過(guò)激波管人體顳骨實(shí)驗(yàn)研究了沖擊波方向與鼓膜破裂閾值的關(guān)系,基于人耳三維有限元建模計(jì)算,提出了鼓膜應(yīng)力相對(duì)于超壓峰值的變化率可以表征鼓膜在沖擊波下的力學(xué)損傷的觀點(diǎn)。盡管閾值指標(biāo)僅選擇了超壓峰值,基于人耳顳骨測(cè)量得到的鼓膜破裂閾值數(shù)據(jù)仍具有臨床診斷價(jià)值。Gan 等[3] 通過(guò)激波管人體顳骨實(shí)驗(yàn)和人耳三維有限元建模評(píng)估了耳塞針對(duì)于沖擊波的防護(hù)機(jī)制,將傳遞到鼓膜處的沖擊波超壓峰值和超壓波形作為評(píng)價(jià)耳塞防護(hù)效果的2 個(gè)指標(biāo)。

    鼓膜由于其低破裂閾值且易觀察的特性,被視為一種爆炸傷標(biāo)記物,用于判斷聽(tīng)力的損害程度。隨著醫(yī)學(xué)研究的深入,研究者提出,鼓膜破裂程度和聽(tīng)力損失程度并不一致,不能準(zhǔn)確地判斷聽(tīng)覺(jué)功能的損害程度。爆炸波如何與鼓膜相互作用并損害外周聽(tīng)覺(jué)系統(tǒng)的機(jī)制仍然不清楚,Leckness 等[4] 建立了人耳三維有限元模型來(lái)預(yù)測(cè)爆炸超壓通過(guò)外耳中耳的傳輸。Brown 等[5-6] 和Bradshaw 等[7] 建立了從外耳到耳蝸的爆炸波傳輸三維有限元模型,并計(jì)算了鼓膜、中耳聽(tīng)小骨和耳蝸基底膜的位移響應(yīng)。有限元模型能模擬沖擊波加載時(shí)中耳和耳蝸的響應(yīng),為闡明爆炸致聽(tīng)力損失機(jī)制提供參考,可應(yīng)用于協(xié)助預(yù)防、診斷和治療戰(zhàn)場(chǎng)人員由于爆炸引起的聽(tīng)力損失。Jiang 等[8-9] 通過(guò)激波管人耳顳骨實(shí)驗(yàn)對(duì)爆炸暴露的人耳進(jìn)行了生物力學(xué)測(cè)量,他們使用雙激光多普勒測(cè)振儀來(lái)測(cè)量爆炸導(dǎo)致的鼓膜以及鐙骨底板的運(yùn)動(dòng)。在此基礎(chǔ)上,Bien 等[10] 同時(shí)測(cè)量了沖擊波導(dǎo)致的鐙骨底板的運(yùn)動(dòng)和耳蝸內(nèi)流體的壓力。這些研究提供了一種新的方法來(lái)定量表征爆炸導(dǎo)致的中耳力學(xué)響應(yīng)并量化進(jìn)入耳蝸的能量通量,實(shí)驗(yàn)數(shù)據(jù)可用于幫助闡明外周聽(tīng)覺(jué)系統(tǒng)損傷的形成機(jī)制。Gan 等[11] 通過(guò)掃描激光多普勒測(cè)振技術(shù)對(duì)正常和爆炸后的人顳骨鼓膜進(jìn)行了全場(chǎng)表面運(yùn)動(dòng)測(cè)量,該方法可描述爆炸后鼓膜傳聲運(yùn)動(dòng)的改變以及可能的局部受損位置。Luo 等[12-13] 通過(guò)霍普金森拉桿在300~2 000 s?1 的高應(yīng)變率下對(duì)人耳鼓膜條帶試樣進(jìn)行了拉伸實(shí)驗(yàn),測(cè)試了爆炸前后鼓膜的高應(yīng)變率力學(xué)性能。Engles 等[14] 通過(guò)激光多普勒測(cè)振儀在200~8 000 Hz 的頻率范圍測(cè)試了聲加載導(dǎo)致的鼓膜振動(dòng),通過(guò)有限元建模的逆問(wèn)題求解確定了爆炸對(duì)鼓膜存儲(chǔ)模量和損耗模量的影響。Liang 等[15-17] 通過(guò)微條紋投影法表征了爆炸沖擊波對(duì)鼓膜非線性靜態(tài)力學(xué)性能的影響,并且表征了人耳圓窗膜的力學(xué)性質(zhì)。

    研究者也通過(guò)激波管動(dòng)物實(shí)驗(yàn)對(duì)聽(tīng)覺(jué)功能損害和聽(tīng)覺(jué)系統(tǒng)病理性變化進(jìn)行了直接表征。Chen 等[18]通過(guò)龍貓爆炸模型量化了重復(fù)低強(qiáng)度爆炸導(dǎo)致的急性和進(jìn)行性聽(tīng)覺(jué)損害。Smith 等[19] 研究了龍貓反復(fù)暴露在引起輕度創(chuàng)傷性腦損傷(mild traumatic brain injury, mTBI)的爆炸載荷下所產(chǎn)生的進(jìn)行性聽(tīng)覺(jué)損害,在出現(xiàn)永久聽(tīng)力損失的動(dòng)物中觀察到了聽(tīng)覺(jué)皮層的損傷,結(jié)果表明引起mTBI 的爆炸次數(shù)嚴(yán)重影響聽(tīng)力損失。Jiang 等[20] 提出了一種新方法來(lái)對(duì)反復(fù)爆炸所導(dǎo)致的外周聽(tīng)覺(jué)系統(tǒng)(peripheral auditorysystem, PAS)和中央聽(tīng)覺(jué)系統(tǒng)(central audirory system, CAS)損傷進(jìn)行造模,通過(guò)使用頭盔以及耳塞的動(dòng)物模型和生物力學(xué)模型對(duì)這2 種損傷進(jìn)行區(qū)分。Shao 等[21] 在此基礎(chǔ)上審查了PAS 和CAS 的改變,發(fā)現(xiàn)中樞聽(tīng)覺(jué)反應(yīng)獨(dú)立于外周聽(tīng)覺(jué),耳塞可以保護(hù)鼓膜和中耳免受結(jié)構(gòu)損傷,但聽(tīng)力水平、耳蝸外毛細(xì)胞和CAS(興奮性和抑制性神經(jīng)遞質(zhì)受體水平)僅在測(cè)試的爆炸超壓水平上得到部分保護(hù)。耳蝸核持續(xù)性的神經(jīng)退化可能在長(zhǎng)期的聽(tīng)力損害中起主要作用。

    然而,沖擊波致聽(tīng)力損失的形成機(jī)制和預(yù)防機(jī)制并不清晰。爆炸致聽(tīng)覺(jué)損害風(fēng)險(xiǎn)準(zhǔn)則仍然存在很多爭(zhēng)議,例如:指標(biāo)選擇沖量還是超壓峰值,正壓持續(xù)時(shí)間是否重要,短脈沖還是長(zhǎng)脈沖更危險(xiǎn)等。關(guān)于閾值的討論組成了一個(gè)龐大的數(shù)據(jù)庫(kù),這個(gè)數(shù)據(jù)庫(kù)由看似不相容的數(shù)據(jù)集和相互矛盾的觀點(diǎn)組成。聽(tīng)力是一個(gè)復(fù)雜的物理過(guò)程,爆炸致聽(tīng)覺(jué)損害風(fēng)險(xiǎn)標(biāo)準(zhǔn)需要足夠復(fù)雜、能夠處理這些數(shù)據(jù)庫(kù)中各種看似矛盾的效應(yīng)。此外,鼓膜作為多層黏彈性膜組織,其對(duì)于爆炸沖擊波和高強(qiáng)聲的失效準(zhǔn)則如何制定仍然面臨挑戰(zhàn)。盡管明確鼓膜破裂與沖擊波波形相關(guān),但沒(méi)有建立相應(yīng)的量效關(guān)系。現(xiàn)有研究絕大部分都是基于小動(dòng)物激波管實(shí)驗(yàn)平臺(tái)開(kāi)展,且鼓膜損傷與耳蝸損傷的關(guān)系沒(méi)有進(jìn)一步表征。

    本研究基于自由場(chǎng)實(shí)爆條件,設(shè)計(jì)并搭建大動(dòng)物爆炸致傷平臺(tái),在不同爆心距離下得到一系列具有超壓峰值梯度和正壓持續(xù)時(shí)間梯度的沖擊波載荷工況,探究不同爆炸參數(shù)對(duì)于鼓膜破裂的影響規(guī)律,并建立基于自由場(chǎng)超壓峰值和正壓持續(xù)時(shí)間的鼓膜創(chuàng)傷量效關(guān)系。

    1 方 法

    1.1 搭建大動(dòng)物致傷平臺(tái)

    采用自由場(chǎng)實(shí)爆條件對(duì)爆炸致動(dòng)物耳創(chuàng)傷進(jìn)行研究,炸藥當(dāng)量分別為1.9 kg TNT 以及8.0 kg TNT,爆心距離地面高度為1.8 m。選用14 頭質(zhì)量均約15 kg 的健康雄性小型豬(巴馬小香豬)作為致傷豬,用于觀察沖擊波傷情,布放方法如圖1 所示。共有7 種致傷工況,每個(gè)工況布置2 頭小型豬。7 種致傷工況的爆炸條件分別為1.9 kg TNT 當(dāng)量爆炸,距離爆心1.8、2.6 和3.2 m 以及8.0 kg TNT 當(dāng)量爆炸,距離爆心2.6、2.9、3.2 和3.8 m,得到一系列具有超壓峰值梯度和正壓時(shí)間梯度的沖擊波載荷工況。

    將焊制的半圓柱形鐵質(zhì)艙室固定在增高支架上用于實(shí)驗(yàn),使艙室所處高度與爆心一致。采用2% 戊巴比妥鈉作為麻醉劑,以20 mg/kg 的劑量對(duì)小型豬進(jìn)行誘導(dǎo)麻醉。將麻醉后的小型豬軀干放置在艙室中,注意頭部面對(duì)爆心,并在縫隙處填充泡沫膠保證密封性,泡沫膠固化后就實(shí)現(xiàn)了小型豬的固定。給固定后的小型豬佩戴防火面罩以保護(hù)口鼻,如圖2 所示。通過(guò)筆形壓力傳感器(型號(hào)PCB137B22/PCB137B23)測(cè)量自由場(chǎng)壓力,采樣頻率為5 MHz。壓力傳感器安裝位置離致傷豬較近,且與致傷豬保持在同一半徑,探頭離地面的高度與爆心一致,且方向與沖擊波傳播方向垂直。

    爆炸前后,均采用聽(tīng)性腦干反應(yīng)(auditory brainstem response, ABR)對(duì)小型豬的聽(tīng)力功能進(jìn)行測(cè)試,記錄爆炸前后ABR 閾值的變化。將小型豬麻醉后轉(zhuǎn)至簡(jiǎn)易隔聲屏蔽室,置于電熱毯上維持體溫。安裝完電極后,以短聲( click)及1、2、4、6 和8 kHz 的短純音作為刺激聲音,最大刺激強(qiáng)度為90 dB,按10 dB 遞減,當(dāng)出現(xiàn)無(wú)規(guī)律難以辨認(rèn)的波形時(shí),遞增5 dB,將誘發(fā)出可重復(fù)規(guī)律波形的最低刺激強(qiáng)度記為ABR 閾值。爆炸后記錄每個(gè)布放距離的即刻死亡率。將麻醉狀態(tài)下的小型豬放置在解剖臺(tái)上,以大動(dòng)脈放血的方式將其處死。沿顱頂中線進(jìn)行開(kāi)顱,取出腦組織后,暴露鼓膜和耳蝸。摘取耳蝸放置于電鏡固定液中保存,通過(guò)掃描電子顯微鏡(scanning electron microscope, SEM)觀察。使用咬骨鉗仔細(xì)剝離鼓膜周圍骨質(zhì),從中耳向外耳道方向暴露鼓膜,觀察記錄鼓膜的創(chuàng)傷情況。

    1.2 自由場(chǎng)爆炸載荷條件

    1.2.1 超壓時(shí)程曲線

    沖擊波超壓時(shí)程(p-t)曲線具有一個(gè)特征形狀,當(dāng)沖擊波到達(dá)時(shí),壓力幾乎是瞬間升高的,然后呈現(xiàn)指數(shù)型衰減,在一段時(shí)間t+后,壓力返回到環(huán)境值并進(jìn)入負(fù)壓階段[22]。Friedlander[23] 提出這種特征形狀可以描述為:

    p = pse-γt/ t+(1-t/ t+) (1)

    式中:ps 為超壓峰值;t+為正壓持續(xù)時(shí)間;γ 為無(wú)量綱常數(shù),這里設(shè)置為4[24]。式(1) 對(duì)較寬峰值超壓范圍的爆炸沖擊波超壓時(shí)間歷史提供了較好的描述,是獲得超壓峰值和正壓持續(xù)時(shí)間的較好方法[22]。

    1.2.2 沖擊波波形分析

    對(duì)超壓波形進(jìn)行脈沖壓力能量頻譜分析,以確定信號(hào)在頻率上的能量分布。通過(guò)傅里葉變換將記錄的超壓波形轉(zhuǎn)化為頻率(25 Hz~5 MHz)上的壓力分布。脈沖能量通量即單位面積的能量,由下式[1]給出:

    式中:E*為脈沖能量通量;p(t) 為時(shí)域上的脈沖壓力信號(hào);T 為信號(hào)總時(shí)長(zhǎng),本文中T=40 ms;ρ 為空氣密度;c 為空氣中的聲速;本文中假設(shè)ρ 和c 是與壓力無(wú)關(guān)的常數(shù)。

    采用中心頻率為62.5、125、250 和500 Hz,1、2、4、8 和16 kHz 的倍頻帶通濾波器,以及低通濾波器L62.5 和高通濾波器H16k,對(duì)不同頻段的信號(hào)進(jìn)行篩選和過(guò)濾。濾波信號(hào)的能量通量也通過(guò)式(2) 來(lái)計(jì)算。由于脈沖信號(hào)在時(shí)域上壓力平方的積分在頻域上保持恒等[25],將各頻段的能量通量除以總脈沖能量通量,得到各個(gè)頻段的歸一化能量通量。

    1.3 基于邏輯回歸模型識(shí)別創(chuàng)傷量效關(guān)系

    與爆炸沖擊波致顱腦創(chuàng)傷(blast-induced traumatic brain injury, bTBI)類似,鼓膜創(chuàng)傷風(fēng)險(xiǎn)曲線可以表示為超壓峰值ps 與正壓持續(xù)時(shí)間t+的函數(shù)。根據(jù)Rafaels 等[26] 的描述,bTBI 耐受性曲線可以通過(guò)冪律關(guān)系進(jìn)行現(xiàn)象學(xué)描述:

    ps = p0[1+α(t+)-β](3)

    式中:p0、α 和β 為需要通過(guò)實(shí)驗(yàn)確定的參數(shù),且β>0。本文中通過(guò)式(3) 對(duì)小型豬耳鼓膜致傷載荷的特點(diǎn)進(jìn)行描述。采用Sigmoid 函數(shù)表示損傷概率:

    Y (A) =1/1+e-A (4)

    式中:Y 為損傷概率,其值限制在0~1 之間;A 為超壓峰值ps 和正壓持續(xù)時(shí)間t+共同貢獻(xiàn)的度量。A 定義

    為如下形式[27]:

    A = C1 +C2 ln{p0[1+α(t+)-β]=ps}(5)

    式中:除p0、α、β 外,C1 和C2 也為待定常數(shù)。通過(guò)最大似然估計(jì)對(duì)實(shí)驗(yàn)中測(cè)得的創(chuàng)傷數(shù)據(jù)進(jìn)行二元邏輯回歸分析,確定鼓膜破裂風(fēng)險(xiǎn)曲線,并繪制在ps-t+圖中。

    2 結(jié)果與討論

    2.1 沖擊波載荷特征

    2.1.1 超壓峰值和正壓持續(xù)時(shí)間

    圖3 展示了在1.9 kg TNT 當(dāng)量爆炸距離爆心1.8 m 處的實(shí)測(cè)自由場(chǎng)超壓時(shí)程曲線以及用Friedlander 方程擬合的結(jié)果,可以看出,F(xiàn)riedlander 方程提供了一個(gè)較好的超壓時(shí)間歷史的描述。由此確定沖擊波超壓峰值和正壓持續(xù)時(shí)間,如表1 所示。由表1 可以看出,在爆炸沖擊波作用范圍1.8~5.5 m 內(nèi),實(shí)驗(yàn)測(cè)得超壓峰值為96.3~628.3 kPa,沖擊波正壓持續(xù)時(shí)間為1.4~6.6 ms。由距離擬合曲線(圖4~5)可見(jiàn),隨著爆心距的增大,爆炸沖擊波超壓峰值減小,正壓持續(xù)時(shí)間增長(zhǎng),符合沖擊波衰減規(guī)律。

    2.1.2 沖擊波能量頻譜分析

    不同當(dāng)量爆炸條件下,不同位置處所測(cè)沖擊波信號(hào)在不同頻帶的歸一化能量分布如圖6 所示。可以看出,沖擊波大部分能量集中在小于500 Hz 的較低頻率下;距離爆心越近,正壓持續(xù)時(shí)間越短,能量分布向高頻移動(dòng)。這表明,不同正壓持續(xù)時(shí)間的沖擊波表現(xiàn)出不同的能量分布特征,正壓持續(xù)時(shí)間越短,高頻段能量占比相對(duì)越高。

    2.2 爆炸致鼓膜破裂量效關(guān)系

    2.2.1 鼓膜創(chuàng)傷形式

    實(shí)驗(yàn)中記錄到的鼓膜損傷情況如圖7 所示,可分為3 類,分別是鼓膜無(wú)明顯損傷(圖7(a))、鼓膜片狀出血(圖7(b))和鼓膜穿孔(圖7(c))?;诖朔诸惤y(tǒng)計(jì)的鼓膜損傷情況如表1 所示,可以看出,當(dāng)自由場(chǎng)峰值超壓小于170 kPa 時(shí),鼓膜無(wú)明顯損傷,為半透明狀薄膜,呈現(xiàn)張緊狀態(tài)。當(dāng)峰值超壓大于237 kPa時(shí),鼓膜會(huì)出現(xiàn)不同程度的充血甚至破裂。本研究中觀察到的鼓膜破裂呈現(xiàn)圓孔狀,占完整鼓膜面積的1/3~2/3,甚至全部消失。沒(méi)有觀察到Koike 等[28] 在臨界破裂壓力觀察到的沿徑向方向的撕裂以及靜壓作用下鼓膜的撕裂[29]。隨著離爆心距離越近,超壓峰值越大,但鼓膜創(chuàng)傷嚴(yán)重程度并未隨之單調(diào)增加。在8.0 kg TNT 當(dāng)量爆炸實(shí)驗(yàn)中,鼓膜破裂嚴(yán)重程度隨著爆心距的減小呈現(xiàn)先提高再降低的趨勢(shì)。

    2.2.2 沖擊波致鼓膜破裂風(fēng)險(xiǎn)曲線

    利用式(4)~(5),以超壓峰值和正壓持續(xù)時(shí)間為自變量,對(duì)表1 中實(shí)驗(yàn)測(cè)得的鼓膜破裂率進(jìn)行二元邏輯回歸分析,結(jié)果如圖8 所示。參數(shù)值確定如下:α=20 ms1.2, β=1.2, p0=67 kPa, C1=2, C2=?8.12。

    對(duì)99%、50% 和1% 的致傷概率,繪制了沖擊波對(duì)小型豬鼓膜破裂風(fēng)險(xiǎn)曲線,同時(shí)也將實(shí)驗(yàn)測(cè)得的數(shù)據(jù)繪制其中,如圖9 所示。隨爆心距的增大,超壓峰值降低,正壓持續(xù)時(shí)間增長(zhǎng)。從建立的量效關(guān)系可以看出,1.9 kg TNT 當(dāng)量爆炸實(shí)驗(yàn)中,鼓膜致傷概率隨爆心距的增大而降低,而在8.0 kg TNT 當(dāng)量爆炸實(shí)驗(yàn)中,鼓膜致傷概率隨著爆心距的增大呈現(xiàn)先提高再降低的趨勢(shì)。和實(shí)驗(yàn)觀察到的規(guī)律一致,回歸模型能捕捉到鼓膜破裂量效關(guān)系主要趨勢(shì)和特征。

    2.3 討論

    自由場(chǎng)沖擊波的特點(diǎn)是上升沿陡峭、超壓峰值高、信號(hào)帶寬和持續(xù)時(shí)間短(微秒-毫秒級(jí))[30]。常用的沖擊波生物毀傷準(zhǔn)則包括超壓峰值準(zhǔn)則、結(jié)合持續(xù)時(shí)間的超壓峰值準(zhǔn)則、沖量準(zhǔn)則和超壓峰值-沖量準(zhǔn)則,其中超壓峰值準(zhǔn)則使用最廣泛[30]。然而,聽(tīng)力損失是一個(gè)復(fù)雜的物理過(guò)程,爆炸致聽(tīng)覺(jué)損害風(fēng)險(xiǎn)準(zhǔn)則仍然存在很多爭(zhēng)議,例如:指標(biāo)選擇沖量還是超壓峰值、正壓持續(xù)時(shí)間是否重要和短脈沖還是長(zhǎng)脈沖更危險(xiǎn)等[31]。本研究中就出現(xiàn)了與直覺(jué)相違背的現(xiàn)象。在8 kg TNT 當(dāng)量的自由場(chǎng)爆炸中,鼓膜創(chuàng)傷程度沒(méi)有隨著超壓峰值的變大(爆心距的減?。┒鴨握{(diào)提高。在離爆心2.6 m 的位置,即使致死率為50%,鼓膜破裂率為0%,充血率為25%。此時(shí)超壓峰值為628.28 kPa,正壓持續(xù)時(shí)間為1.3 ms,這一沖擊波強(qiáng)度足以使小型豬致死,卻沒(méi)有出現(xiàn)鼓膜破裂。隨著爆心距的增大,沒(méi)有出現(xiàn)小型豬死亡的情況。在爆心距2.9 m 處,鼓膜破裂率為0%,充血率為100%,此時(shí)超壓峰值為528.74 kPa,正壓持續(xù)時(shí)間為2.11 ms。在爆心距3.2 m 處,鼓膜破裂率提高到100%,在爆心距3.8 m 處,鼓膜破裂率降低到50%,充血率為25%。

    鼓膜作為一個(gè)黏彈性薄膜結(jié)構(gòu),通過(guò)自身的振動(dòng)將聲音信號(hào)傳遞進(jìn)中耳及內(nèi)耳。Cullis[32] 討論了沖擊波與結(jié)構(gòu)相互作用的影響因素,認(rèn)為爆炸形成的沖擊波將來(lái)自爆炸流場(chǎng)的能量耦合到結(jié)構(gòu)中而與結(jié)構(gòu)相互作用。隨著沖擊波的傳播,其超壓峰值降低,正壓區(qū)被不斷拉寬,形狀發(fā)生改變,相應(yīng)地沖擊波載荷的能量頻率分布也發(fā)生了變化。圖6 也能提供直接的依據(jù),離爆心越遠(yuǎn),低頻能量占比越高。然而,受載結(jié)構(gòu)具有一系列固定的特征頻率及相應(yīng)的振動(dòng)模式,結(jié)構(gòu)的響應(yīng)依賴于沖擊波如何將能量耦合到這些振動(dòng)模式中[32]。因此,除了超壓峰值和沖量,沖擊波的整體形狀也對(duì)結(jié)構(gòu)的響應(yīng)產(chǎn)生顯著影響。此外,鼓膜及其附著的聽(tīng)骨鏈、肌肉和韌帶都有黏彈性,具有一定的壓力緩沖能力[33]。當(dāng)沖擊波的正壓持續(xù)時(shí)間與結(jié)構(gòu)振蕩周期相比較短時(shí),載荷的部分能量就會(huì)被結(jié)構(gòu)的慣性吸收,從而減小結(jié)構(gòu)變形[32]。這些機(jī)制提供了一定的角度去理解鼓膜的創(chuàng)傷程度是如何隨爆心距的改變而改變的。然而,確切的致傷機(jī)制需要進(jìn)一步研究。一方面,可以通過(guò)實(shí)驗(yàn)測(cè)量或者有限元仿真去定量研究鼓膜在不同類型沖擊波作用下的動(dòng)態(tài)響應(yīng)特性;另一方面,需要發(fā)展能反映鼓膜材料動(dòng)態(tài)損傷特征的失效準(zhǔn)則。

    實(shí)驗(yàn)中發(fā)現(xiàn),盡管沒(méi)有發(fā)生鼓膜穿孔,卻仍然出現(xiàn)了嚴(yán)重的聽(tīng)力損失和內(nèi)耳耳蝸損傷,如圖10 所示。圖10 為8.0 kg TNT 當(dāng)量爆炸條件下距爆心2.9 m 處的小型豬在爆炸前后的ABR 聽(tīng)力閾值變化以及耳蝸毛細(xì)胞SEM 圖像。可以看出,短聲ABR 閾值上升40 dB,高頻聽(tīng)力損失嚴(yán)重,2 kHz 以上純音ABR 閾值上升40~70 dB,同時(shí)也觀測(cè)到耳蝸內(nèi)毛細(xì)胞(inner hair cell, IHC)與外毛細(xì)胞(outer hair cell,OHC)的顯著損傷。這一現(xiàn)象也印證了Xie 等[34] 和Littlefield 等[35] 提出的鼓膜破裂和聽(tīng)力功能損失不一致的現(xiàn)象,鼓膜的破裂程度不能直接作為聽(tīng)覺(jué)爆炸傷標(biāo)記物。Amrein 等[31] 指出,即使正壓持續(xù)時(shí)間較短,也會(huì)立即出現(xiàn)嚴(yán)重的聽(tīng)力損害。Koike 等[28] 提出鼓膜破裂可以被視為耳部的一種合適的保護(hù)機(jī)制。鼓膜與圓窗膜的面積比能夠放大鼓膜所接收到的聲壓,實(shí)現(xiàn)外耳空氣聲與內(nèi)耳流體聲的阻抗匹配,使得大部分聲能可以傳遞進(jìn)內(nèi)耳[36]。如果鼓膜是非常頑強(qiáng)的,只會(huì)在顯著超壓下破裂,那么突然的非生理性超壓會(huì)導(dǎo)致壓力波的能量直接傳遞進(jìn)入聽(tīng)骨鏈以及耳蝸。這會(huì)導(dǎo)致更加嚴(yán)重的問(wèn)題,例如:聽(tīng)骨鏈的脫位和圓窗膜的破裂等。但是,鼓膜破裂是一種機(jī)械損傷,并不是所有高強(qiáng)度的沖擊波都會(huì)導(dǎo)致鼓膜破裂,例如:在本研究8 kg TNT 當(dāng)量的實(shí)爆實(shí)驗(yàn)中,距爆心2.6 和2.9 m 處,正壓持續(xù)時(shí)間短,即使超壓峰值高,也沒(méi)有出現(xiàn)鼓膜破裂。因此,具有高聲強(qiáng)、短正壓時(shí)間特征的載荷(例如小口徑武器產(chǎn)生的高強(qiáng)度短脈沖[37])對(duì)于聽(tīng)覺(jué)功能來(lái)說(shuō)是危險(xiǎn)的,在防護(hù)的過(guò)程中需要特別注意。這種載荷一般不會(huì)引起鼓膜破裂,容易受到忽略,類似的載荷還有低強(qiáng)度沖擊波的反復(fù)暴露(低于鼓膜破裂閾值)[18]。盡管鼓膜破裂程度不能作為長(zhǎng)期聽(tīng)覺(jué)功能障礙的標(biāo)記物,但鼓膜仍然是沖擊波與耳部作用的第一步。沖擊波如何與鼓膜相互作用并將能量傳遞到外周聽(tīng)覺(jué)系統(tǒng)的機(jī)制仍然需要研究,這對(duì)于闡明強(qiáng)脈沖聲致聽(tīng)力損失的形成機(jī)制和預(yù)防機(jī)制至關(guān)重要。

    3 結(jié) 論

    基于自由場(chǎng)實(shí)爆條件,設(shè)計(jì)并搭建了大動(dòng)物爆炸致傷平臺(tái),在不同爆心距離下得到了一系列具有超壓峰值梯度和正壓持續(xù)時(shí)間梯度的沖擊波載荷工況?;诒緦?shí)驗(yàn)中小型豬耳部傷情的觀測(cè)結(jié)果,探究了不同爆炸參數(shù)對(duì)鼓膜破裂的影響規(guī)律,通過(guò)邏輯回歸分析建立了基于自由場(chǎng)超壓峰值和正壓持續(xù)時(shí)間的鼓膜創(chuàng)傷量效關(guān)系。研究發(fā)現(xiàn),當(dāng)自由場(chǎng)超壓峰值低于170 kPa 時(shí),鼓膜無(wú)明顯損傷;當(dāng)自由場(chǎng)超壓峰值大于237 kPa 時(shí),部分鼓膜出現(xiàn)不同程度的破裂和充血。隨著離爆心越近,超壓峰值越大,但鼓膜創(chuàng)傷的嚴(yán)重程度并未隨之單調(diào)提高。通過(guò)分析沖擊波載荷特征可知,離爆心越近正壓持續(xù)時(shí)間越短,高頻段能量占比相對(duì)越高,小型豬鼓膜破裂的概率可能反而降低,此時(shí),仍然出現(xiàn)顯著的聽(tīng)力損失和耳蝸損傷。除了超壓峰值,沖擊波波形頻譜分布可能對(duì)鼓膜破裂程度影響顯著。爆炸波如何與鼓膜相互作用并損害聽(tīng)覺(jué)系統(tǒng)的機(jī)制仍需進(jìn)一步研究。本研究可為爆炸性武器沖擊波生物損傷效應(yīng)標(biāo)準(zhǔn)的制定提供參考。

    參考文獻(xiàn):

    [1]GAN R Z, NAKMALI"D, JI X D, et al. Mechanical damage of tympanic membrane in relation to impulse pressure waveform:a study in chinchillas [J]. Hearing Research, 2016, 340: 25–34. DOI: 10.1016/j.heares.2016.01.004.

    [2]GAN R Z, LECKNESS K, NAKMALI D, et al. Biomechanical measurement and modeling of human eardrum injury inrelation to blast wave direction [J]. Military Medicine, 2018, 183(S1): 245–251. DOI: 10.1093/milmed/usx149.

    [3]GAN R Z, LECKNESS K, SMITH K, et al. Characterization of protection mechanisms to blast overpressure for personalhearing protection devices: biomechanical measurement and computational modeling [J]. Military Medicine, 2019, 184(S1):251–260. DOI: 10.1093/milmed/usy299.

    [4]LECKNESS K, NAKMALI D, GAN R Z. Computational modeling of blast wave transmission through human ear [J].Military Medicine, 2018, 183(S1): 262–268. DOI: 10.1093/milmed/usx226.

    [5]BROWN M A, JI X D, GAN R Z. 3D finite element modeling of blast wave transmission from the external ear to cochlea [J].Annals of Biomedical Engineering, 2021, 49(2): 757–768. DOI: 10.1007/s10439-020-02612-y.

    [6]BROWN M A, BRADSHAW J J, GAN R Z. Three-dimensional finite element modeling of blast wave transmission from theexternal ear to a spiral cochlea [J]. Journal of Biomechanical Engineering, 2022, 144(1): 014503. DOI: 10.1115/1.4051925.

    [7]BRADSHAW J J, BROWN M A, JIANG S Y, et al. 3D finite element model of human ear with 3-chamber spiral cochlea forblast wave transmission from the ear canal to cochlea [J]. Annals of Biomedical Engineering, 2023, 51(5): 1106–1118. DOI:10.1007/s10439-023-03200-6.

    [8]JIANG S Y, SMITH K, GAN R Z. Dual-laser measurement and finite element modeling of human tympanic membranemotion under blast exposure [J]. Hearing Research, 2019, 378: 43–52. DOI: 10.1016/j.heares.2018.12.003.

    [9]JIANG S Y, DAI C K, GAN R Z. Dual-laser measurement of human stapes footplate motion under blast exposure [J]. HearingResearch, 2021, 403: 108177. DOI: 10.1016/j.heares.2021.108177.

    [10]BIEN A G, JIANG S Y, GAN R Z. Real-time measurement of stapes motion and intracochlear pressure during blast exposure [J].Hearing Research, 2023, 429: 108702. DOI: 10.1016/j.heares.2023.108702.

    [11]GAN R Z, JIANG S Y. Surface motion changes of tympanic membrane damaged by blast waves [J]. Journal of BiomechanicalEngineering, 2019, 141(9): 091009. DOI: 10.1115/1.4044052.

    [12]LUO H Y, DAI C K, GAN R Z, et al. Measurement of Young’s modulus of human tympanic membrane at high strain rates [J].Journal of Biomechanical Engineering, 2009, 131(6): 064501. DOI: 10.1115/1.3118770.

    [13]LUO H Y, JIANG S Y, NAKMALI D U, et al. Mechanical properties of a human eardrum at high strain rates after exposure toblast waves [J]. Journal of Dynamic Behavior of Materials, 2016, 2(1): 59–73. DOI: 10.1007/s40870-015-0041-3.

    [14]ENGLES W G, WANG X L, GAN R Z. Dynamic properties of human tympanic membrane after exposure to blast waves [J].Annals of Biomedical Engineering, 2017, 45(10): 2383–2394. DOI: 10.1007/s10439-017-1870-0.

    [15]LIANG J F, LUO H Y, YOKELL Z, et al. Characterization of the nonlinear elastic behavior of chinchilla tympanic membraneusing micro-fringe projection [J]. Hearing Research, 2016, 339: 1–11. DOI: 10.1016/j.heares.2016.05.012.

    [16]LIANG J F, YOKELL Z A, NAKMAILI D U, et al. The effect of blast overpressure on the mechanical properties of achinchilla tympanic membrane [J]. Hearing Research, 2017, 354: 48–55. DOI: 10.1016/j.heares.2017.08.003.

    [17]LIANG J F, SMITH K D, GAN R Z, et al. The effect of blast overpressure on the mechanical properties of the humantympanic membrane [J]. Journal of the Mechanical Behavior of Biomedical Materials, 2019, 100: 103368. DOI: 10.1016/j.jmbbm.2019.07.026.

    [18]CHEN T, SMITH K, JIANG S Y, et al. Progressive hearing damage after exposure to repeated low-intensity blasts inchinchillas [J]. Hearing Research, 2019, 378: 33–42. DOI: 10.1016/j.heares.2019.01.010.

    [19]SMITH K D, CHEN T, GAN R Z. Hearing damage induced by blast overpressure at mild TBI level in a chinchilla model [J].Military Medicine, 2020, 185(S1): 248–255. DOI: 10.1093/milmed/usz309.

    [20]JIANG S Y, GANNON A N, SMITH K D, et al. Prevention of blast-induced auditory injury using 3D printed helmet andhearing protection device: a preliminary study on biomechanical modeling and animal [J]. Military Medicine, 2021, 186(S1):537–545. DOI: 10.1093/milmed/usaa317.

    [21]SHAO N N, JIANG S Y, YOUNGER D, et al. Central and peripheral auditory abnormalities in chinchilla animal model ofblast-injury [J]. Hearing Research, 2021, 407: 108273. DOI: 10.1016/j.heares.2021.108273.

    [22]DEWEY J M. The shape of the blast wave: studies of the Friedlander equation [C]//Proceedings of the 21st InternationalSymposium on Military Aspects of Blast and Shock. 2010: 1–9.

    [23]FRIEDLANDER F G. The diffraction [23] of sound pulses Ⅰ: diffraction by a semi-infinite plane [J]. Proceedings of the Royal Society of London Series A: Mathematical and Physical Sciences, 1946, 186(1006): 322–344. DOI: 10.1098/rspa.1946.0046.

    [24]IYOHO A E, HO K, CHAN P. The development of a tympanic membrane model and probabilistic dose-response riskassessment of rupture because of blast [J]. Military Medicine, 2020, 185(S1): 234–242. DOI: 10.1093/milmed/usz215.

    [25]YOUNG R W. On the energy transported with a sound pulse [J]. The Journal of the Acoustical Society of America, 1970,47(2A): 441–442. DOI: 10.1121/1.1911547.

    [26]RAFAELS K, BASS C, SALZAR R S, et al. Survival risk assessment for primary blast exposures to the head [J]. Journal ofNeurotrauma, 2011, 28: 2319–2328. DOI: 10.1089/neu.2009.1207.

    [27]ZHU F, CHOU C C, YANG K H, et al. Some considerations on the threshold and inter-species scaling law for primary blastinducedtraumatic brain injury: a semi-analytical approach [J]. Journal of Mechanics in Medicine and Biology, 2013, 13(4):1350065. DOI: 10.1142/S0219519413500656.

    [28]KOIKE T, WADA H, ITO F, et al. High-speed video observation of tympanic membrane rupture in guinea pigs [J]. JSMEInternational Journal Series C: Mechanical Systems, Machine Elements and Manufacturing, 2003, 46(4): 1434–1440. DOI:10.1299/jsmec.46.1434.

    [29]JENSEN J H, BONDING P. Experimental pressure induced rupture of the tympanic membrane in man [J]. Acta Oto-Laryngologica, 1993, 113(1): 62–67. DOI: 10.3109/00016489309135768.

    [30]張潔元. 爆炸沖擊波致豚鼠聽(tīng)覺(jué)和前庭功能損害評(píng)估技術(shù)與標(biāo)準(zhǔn)研究 [D]. 重慶: 陸軍軍醫(yī)大學(xué), 2022. DOI:10.27001/d.cnki.gtjyu.2022.000239.

    ZHANG J Y. Evaluation techniques and criteria for auditory and vestibular function impairment caused by blast shock wavesin guinea pigs [D]. Chongqing: Army Medical University, 2022. DOI: 10.27001/d.cnki.gtjyu.2022.000239.

    [31]AMREIN B E, LETOWSKI T R. High level impulse sounds and human hearing: standards, physiology, quantification: ARLTR-6017 [R]. Affiliation: U. S. Army Research Laboratory, 2012.

    [32]CULLIS I G. Blast waves and how they interact with structures [J]. BMJ Military Health, 2001, 147(1): 16–26. DOI:10.1136/jramc-147-01-02.

    [33]PADURARIU S, DE GREEF D, JACOBSEN H, et al. Pressure buffering by the tympanic membrane: in vivo measurementsof middle ear pressure fluctuations during elevator motion [J]. Hearing Research, 2016, 340: 113–120. DOI: 10.1016/j.heares.2015.12.004.

    [34]XIE P P, PENG Y, HU J J, et al. Assessment of hearing loss induced by tympanic membrane perforations under blastenvironment [J]. European Archives of Oto-Rhino-Laryngology, 2020, 277(2): 453–461. DOI: 10.1007/s00405-019-05710-3.

    [35]LITTLEFIELD P D, BRUNGART D S. Long-term sensorineural hearing loss in patients with blast-induced tympanicmembrane perforations [J]. Ear and Hearing, 2020, 41(1): 165–172. DOI: 10.1097/AUD.0000000000000751.

    [36]FAY J, PURIA S, DECRAEMER W F, et al. Three approaches for estimating the elastic modulus of the tympanic membrane [J].Journal of Biomechanics, 2005, 38(9): 1807–1815. DOI: 10.1016/j.jbiomech.2004.08.022.

    [37]PRICE G R, KIM H N, LIM D J, et al. Hazard from weapons impulses: histological and electrophysiological evidence [J]. TheJournal of the Acoustical Society of America, 1989, 85(3): 1245–1254. DOI: 10.1121/1.397455.

    (責(zé)任編輯 張凌云)

    基金項(xiàng)目: 國(guó)家重點(diǎn)研發(fā)計(jì)劃(2020-JCJQ-ZD-254,2022YFC3320502)

    猜你喜歡
    沖擊波
    水下爆破沖擊波對(duì)混凝土塢門(mén)影響試驗(yàn)研究
    經(jīng)筋排刺法結(jié)合體外沖擊波治療中老年肩周炎的臨床觀察
    防護(hù)裝置粘接強(qiáng)度對(duì)爆炸切割沖擊波的影響
    武漢沖擊波
    能源物聯(lián)網(wǎng)沖擊波
    能源(2018年10期)2018-12-08 08:02:34
    醫(yī)生集團(tuán)沖擊波
    放散式體外沖擊波治療鈣化性岡上肌肌腱炎的臨床療效
    放散式?jīng)_擊波治療早中期股骨頭缺血性壞死臨床療效觀察
    牽引定位下沖擊波治療股骨頭缺血性壞死72例
    超聲雙探頭聯(lián)合定位法在體外沖擊波碎石術(shù)中的應(yīng)用
    成人美女网站在线观看视频| 亚洲成人一二三区av| 中文字幕av成人在线电影| 免费大片黄手机在线观看| av国产久精品久网站免费入址| 日日摸夜夜添夜夜添av毛片| 美女xxoo啪啪120秒动态图| 国产淫片久久久久久久久| 2021少妇久久久久久久久久久| 中文资源天堂在线| 国产成人福利小说| 亚洲电影在线观看av| 久久久久网色| 99热6这里只有精品| 神马国产精品三级电影在线观看| 大香蕉久久网| 91精品伊人久久大香线蕉| 美女国产视频在线观看| 麻豆久久精品国产亚洲av| xxx大片免费视频| 国产美女午夜福利| 99精国产麻豆久久婷婷| 丝瓜视频免费看黄片| 一本久久精品| 一级a做视频免费观看| 九九久久精品国产亚洲av麻豆| 中文资源天堂在线| 波多野结衣巨乳人妻| 欧美成人精品欧美一级黄| 在线看a的网站| 亚洲精华国产精华液的使用体验| 国产在视频线精品| 亚洲国产欧美人成| 婷婷色综合大香蕉| 97超碰精品成人国产| 亚洲国产成人一精品久久久| 激情五月婷婷亚洲| 国产极品天堂在线| 王馨瑶露胸无遮挡在线观看| 欧美日韩视频精品一区| 国内精品美女久久久久久| 成年女人在线观看亚洲视频 | 国产亚洲av嫩草精品影院| 亚洲国产色片| 一本一本综合久久| 免费大片18禁| 欧美日本视频| 国产 一区 欧美 日韩| 一级毛片久久久久久久久女| 色哟哟·www| 一本色道久久久久久精品综合| 欧美激情久久久久久爽电影| 国产精品伦人一区二区| 日本猛色少妇xxxxx猛交久久| 深夜a级毛片| 亚洲欧美成人精品一区二区| 99热这里只有精品一区| 免费看日本二区| 日韩 亚洲 欧美在线| 国产男女内射视频| 国产精品成人在线| 国产精品偷伦视频观看了| 一个人看视频在线观看www免费| 真实男女啪啪啪动态图| 麻豆精品久久久久久蜜桃| 少妇猛男粗大的猛烈进出视频 | 亚洲自偷自拍三级| 国产一区二区三区av在线| 亚洲av国产av综合av卡| 午夜视频国产福利| 伊人久久国产一区二区| 成人二区视频| 日韩免费高清中文字幕av| 久久人人爽人人爽人人片va| 少妇人妻 视频| 精品人妻视频免费看| 大陆偷拍与自拍| 欧美潮喷喷水| 亚洲精品456在线播放app| 亚洲欧美清纯卡通| 综合色av麻豆| 丰满人妻一区二区三区视频av| 欧美日韩在线观看h| 亚洲av中文字字幕乱码综合| 成人综合一区亚洲| 精品少妇黑人巨大在线播放| 午夜免费男女啪啪视频观看| 亚洲av男天堂| 久久99蜜桃精品久久| 日本黄大片高清| 成人亚洲精品av一区二区| 国产精品一二三区在线看| 国产精品国产三级国产专区5o| av国产精品久久久久影院| 欧美性感艳星| 国产成人精品婷婷| 精华霜和精华液先用哪个| 1000部很黄的大片| 看免费成人av毛片| 中国美白少妇内射xxxbb| 尤物成人国产欧美一区二区三区| 中文字幕免费在线视频6| 国产亚洲av嫩草精品影院| 久久久久久伊人网av| 交换朋友夫妻互换小说| 别揉我奶头 嗯啊视频| 欧美 日韩 精品 国产| 美女内射精品一级片tv| eeuss影院久久| 91久久精品电影网| 欧美精品人与动牲交sv欧美| 欧美 日韩 精品 国产| 在线免费观看不下载黄p国产| 久久久精品欧美日韩精品| 日日啪夜夜爽| 亚洲人与动物交配视频| 亚洲精品一二三| 91狼人影院| 国产亚洲一区二区精品| av在线播放精品| 女人十人毛片免费观看3o分钟| 啦啦啦中文免费视频观看日本| 丝袜喷水一区| 久久国内精品自在自线图片| 在线看a的网站| av又黄又爽大尺度在线免费看| 蜜臀久久99精品久久宅男| 在线观看一区二区三区激情| 伦理电影大哥的女人| 91午夜精品亚洲一区二区三区| 99热国产这里只有精品6| 九色成人免费人妻av| 春色校园在线视频观看| 美女视频免费永久观看网站| 人妻一区二区av| 22中文网久久字幕| 日韩三级伦理在线观看| 性插视频无遮挡在线免费观看| 免费黄色在线免费观看| 蜜桃久久精品国产亚洲av| 99久久精品一区二区三区| 国产探花极品一区二区| 男女啪啪激烈高潮av片| 亚洲av男天堂| av在线播放精品| 亚洲欧美日韩另类电影网站 | 免费电影在线观看免费观看| 91精品伊人久久大香线蕉| 久久久久久久大尺度免费视频| 小蜜桃在线观看免费完整版高清| 一级毛片黄色毛片免费观看视频| 99久久精品热视频| 亚洲一区二区三区欧美精品 | 亚洲国产高清在线一区二区三| 草草在线视频免费看| 国产精品爽爽va在线观看网站| 精品国产露脸久久av麻豆| 一个人看视频在线观看www免费| 天堂网av新在线| 国产精品久久久久久久电影| 亚洲国产高清在线一区二区三| 亚洲av成人精品一二三区| 久久久久久九九精品二区国产| 熟女av电影| 久久99热这里只频精品6学生| 亚洲电影在线观看av| 十八禁网站网址无遮挡 | 自拍欧美九色日韩亚洲蝌蚪91 | 神马国产精品三级电影在线观看| 国产精品人妻久久久久久| 黄色怎么调成土黄色| 欧美日韩精品成人综合77777| 日本-黄色视频高清免费观看| 国产伦在线观看视频一区| 超碰97精品在线观看| 国产亚洲精品久久久com| 国产精品一区www在线观看| 美女高潮的动态| 国产精品伦人一区二区| 欧美成人午夜免费资源| 亚洲av免费在线观看| 欧美日韩视频精品一区| 久久久久性生活片| 嘟嘟电影网在线观看| 久久久久国产精品人妻一区二区| 人妻少妇偷人精品九色| 插阴视频在线观看视频| 亚洲美女视频黄频| 人妻一区二区av| 午夜福利在线观看免费完整高清在| 日韩电影二区| 少妇 在线观看| 国产熟女欧美一区二区| 伦精品一区二区三区| 制服丝袜香蕉在线| 久久久久久久午夜电影| 三级男女做爰猛烈吃奶摸视频| 亚洲人与动物交配视频| 国国产精品蜜臀av免费| 五月开心婷婷网| 亚洲婷婷狠狠爱综合网| 国产精品三级大全| 国产精品一区二区性色av| 国产午夜精品久久久久久一区二区三区| 97在线视频观看| 亚洲精品第二区| 亚洲精华国产精华液的使用体验| av播播在线观看一区| 精品少妇久久久久久888优播| 久久99热6这里只有精品| www.av在线官网国产| 国产成人免费观看mmmm| 下体分泌物呈黄色| 国产成人午夜福利电影在线观看| 亚洲av国产av综合av卡| 在线观看国产h片| 免费看av在线观看网站| 国产伦理片在线播放av一区| av免费观看日本| 亚洲精品国产成人久久av| 韩国av在线不卡| 另类亚洲欧美激情| 国产熟女欧美一区二区| 大片电影免费在线观看免费| 狂野欧美激情性bbbbbb| 日韩成人伦理影院| 国产真实伦视频高清在线观看| 成人国产av品久久久| 成人一区二区视频在线观看| 亚洲av国产av综合av卡| 在线观看国产h片| 男人狂女人下面高潮的视频| 日本一本二区三区精品| 高清av免费在线| 插逼视频在线观看| 国产精品一区二区三区四区免费观看| 亚洲av中文字字幕乱码综合| 白带黄色成豆腐渣| 欧美老熟妇乱子伦牲交| 免费黄色在线免费观看| 国产精品国产三级国产av玫瑰| 国产极品天堂在线| 日本黄大片高清| 五月伊人婷婷丁香| 亚洲一区二区三区欧美精品 | 日韩一区二区三区影片| 国产av不卡久久| 99精国产麻豆久久婷婷| 一区二区三区乱码不卡18| 超碰av人人做人人爽久久| 天天躁日日操中文字幕| 一级毛片黄色毛片免费观看视频| 水蜜桃什么品种好| 亚洲,一卡二卡三卡| 22中文网久久字幕| 国产精品国产av在线观看| 久久久久久久久久成人| 亚洲精品久久久久久婷婷小说| 久久久精品94久久精品| 特级一级黄色大片| 97超视频在线观看视频| 大陆偷拍与自拍| 亚洲欧美成人精品一区二区| 国产成人精品福利久久| 看非洲黑人一级黄片| 日韩成人伦理影院| 亚洲精品456在线播放app| 免费观看在线日韩| 国模一区二区三区四区视频| 看非洲黑人一级黄片| 人人妻人人看人人澡| 国产探花在线观看一区二区| 自拍欧美九色日韩亚洲蝌蚪91 | 免费看光身美女| 女的被弄到高潮叫床怎么办| 国产大屁股一区二区在线视频| 欧美人与善性xxx| 91狼人影院| 少妇裸体淫交视频免费看高清| videossex国产| 国产精品嫩草影院av在线观看| 3wmmmm亚洲av在线观看| 亚洲美女视频黄频| 蜜臀久久99精品久久宅男| 亚洲不卡免费看| 九九爱精品视频在线观看| 国产精品99久久久久久久久| 国产黄片美女视频| 久久精品国产a三级三级三级| 免费看不卡的av| 亚洲在久久综合| 黄色视频在线播放观看不卡| 97超视频在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品自拍成人| 天天躁夜夜躁狠狠久久av| h日本视频在线播放| 人体艺术视频欧美日本| 国产成人freesex在线| 建设人人有责人人尽责人人享有的 | 久久这里有精品视频免费| 欧美最新免费一区二区三区| 亚洲天堂国产精品一区在线| 男女那种视频在线观看| 三级经典国产精品| 麻豆成人午夜福利视频| 亚洲人成网站在线播| 又爽又黄a免费视频| 久久亚洲国产成人精品v| 亚洲国产高清在线一区二区三| 两个人的视频大全免费| 亚洲国产精品999| 亚洲天堂国产精品一区在线| 久热这里只有精品99| 国产探花在线观看一区二区| 亚洲久久久久久中文字幕| 久久久久久久大尺度免费视频| 老司机影院成人| 特大巨黑吊av在线直播| 久久久久九九精品影院| 欧美亚洲 丝袜 人妻 在线| 国产成人91sexporn| 国内精品美女久久久久久| 久热这里只有精品99| 日韩欧美精品v在线| 老司机影院毛片| 男人舔奶头视频| 国产精品无大码| 高清av免费在线| 国产一区亚洲一区在线观看| 伦精品一区二区三区| 亚洲人成网站在线观看播放| 丝袜美腿在线中文| 蜜桃亚洲精品一区二区三区| 18禁在线无遮挡免费观看视频| 在线观看美女被高潮喷水网站| 白带黄色成豆腐渣| 亚洲精品国产av蜜桃| 精品一区二区三区视频在线| 中国国产av一级| 国产男女内射视频| 久久人人爽av亚洲精品天堂 | 久久久久久久久大av| 一级av片app| 久久女婷五月综合色啪小说 | 国产成人精品婷婷| 寂寞人妻少妇视频99o| 直男gayav资源| 免费观看av网站的网址| av在线天堂中文字幕| 亚洲天堂国产精品一区在线| 国产精品福利在线免费观看| 2021少妇久久久久久久久久久| 日日摸夜夜添夜夜添av毛片| 少妇人妻久久综合中文| 日韩一区二区三区影片| 国产有黄有色有爽视频| 国产亚洲最大av| 久久精品国产a三级三级三级| 午夜福利网站1000一区二区三区| 深爱激情五月婷婷| 久久久a久久爽久久v久久| 交换朋友夫妻互换小说| 好男人视频免费观看在线| 波多野结衣巨乳人妻| 精品久久久久久电影网| 国产精品久久久久久精品电影小说 | 亚洲va在线va天堂va国产| 91精品一卡2卡3卡4卡| 婷婷色综合大香蕉| 91久久精品国产一区二区成人| 精品视频人人做人人爽| 成人国产麻豆网| 一级黄片播放器| 国产精品.久久久| 欧美少妇被猛烈插入视频| 欧美 日韩 精品 国产| 亚洲激情五月婷婷啪啪| 国产免费一区二区三区四区乱码| 联通29元200g的流量卡| 日韩伦理黄色片| 99热全是精品| 国产成人精品久久久久久| 69人妻影院| 亚洲av中文字字幕乱码综合| 日本wwww免费看| 韩国高清视频一区二区三区| av.在线天堂| 精品久久国产蜜桃| 欧美 日韩 精品 国产| 交换朋友夫妻互换小说| 午夜亚洲福利在线播放| 日韩欧美一区视频在线观看 | 午夜免费男女啪啪视频观看| 亚洲美女视频黄频| 久久人人爽人人片av| 色网站视频免费| 观看美女的网站| 美女被艹到高潮喷水动态| 婷婷色综合大香蕉| 午夜激情福利司机影院| 免费观看在线日韩| 国产成人免费无遮挡视频| 国产精品人妻久久久久久| 亚洲一级一片aⅴ在线观看| 在线亚洲精品国产二区图片欧美 | 久久综合国产亚洲精品| 亚洲三级黄色毛片| 美女国产视频在线观看| 热re99久久精品国产66热6| 交换朋友夫妻互换小说| av卡一久久| 在线精品无人区一区二区三 | 亚洲真实伦在线观看| 亚洲伊人久久精品综合| 亚洲av免费在线观看| 国产欧美另类精品又又久久亚洲欧美| 最近的中文字幕免费完整| 人体艺术视频欧美日本| 综合色av麻豆| 亚洲三级黄色毛片| 亚洲av.av天堂| 高清在线视频一区二区三区| 亚洲精品成人久久久久久| 韩国av在线不卡| 国产免费视频播放在线视频| 日本一本二区三区精品| 日本wwww免费看| 欧美 日韩 精品 国产| 免费高清在线观看视频在线观看| 毛片女人毛片| 神马国产精品三级电影在线观看| av在线播放精品| 可以在线观看毛片的网站| 日韩精品有码人妻一区| 亚洲精品第二区| 纵有疾风起免费观看全集完整版| 18禁在线无遮挡免费观看视频| 久久精品熟女亚洲av麻豆精品| 国产精品秋霞免费鲁丝片| 狂野欧美激情性bbbbbb| 久久久精品免费免费高清| 色吧在线观看| 免费在线观看成人毛片| 夜夜看夜夜爽夜夜摸| 午夜福利视频精品| 国产又色又爽无遮挡免| 久久久成人免费电影| 99热这里只有精品一区| 国产黄色免费在线视频| 久久97久久精品| 国产淫语在线视频| av免费观看日本| 国产免费一区二区三区四区乱码| 国产在线男女| 国产真实伦视频高清在线观看| 国产国拍精品亚洲av在线观看| 男人爽女人下面视频在线观看| 看免费成人av毛片| 大片电影免费在线观看免费| 97精品久久久久久久久久精品| 精品国产露脸久久av麻豆| 久久久久久久精品精品| 欧美变态另类bdsm刘玥| 久久久国产一区二区| 老女人水多毛片| 亚洲av成人精品一二三区| 成人二区视频| 国产亚洲5aaaaa淫片| 热99国产精品久久久久久7| 国产精品一区二区在线观看99| 久久99热这里只频精品6学生| 午夜免费男女啪啪视频观看| 国国产精品蜜臀av免费| 内地一区二区视频在线| 黄片无遮挡物在线观看| 国产精品一区www在线观看| 国产男女超爽视频在线观看| 啦啦啦中文免费视频观看日本| 日韩欧美一区视频在线观看 | 国产视频内射| 肉色欧美久久久久久久蜜桃 | 欧美日韩综合久久久久久| 观看美女的网站| av在线亚洲专区| 欧美少妇被猛烈插入视频| 看黄色毛片网站| 在线观看一区二区三区| 性色av一级| 深夜a级毛片| 免费看av在线观看网站| 自拍偷自拍亚洲精品老妇| 一级黄片播放器| 国产极品天堂在线| 99热全是精品| 天美传媒精品一区二区| 观看免费一级毛片| 视频区图区小说| 在线免费十八禁| 黄色日韩在线| 国产成人精品婷婷| 午夜精品一区二区三区免费看| 永久网站在线| 在线 av 中文字幕| 色网站视频免费| 亚洲va在线va天堂va国产| 熟女电影av网| 高清欧美精品videossex| 国产熟女欧美一区二区| 国产亚洲精品久久久com| 日本-黄色视频高清免费观看| 亚洲四区av| 黄色配什么色好看| 国产综合精华液| 美女被艹到高潮喷水动态| 国产成人福利小说| 日韩电影二区| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲,一卡二卡三卡| 国产成人精品福利久久| 国精品久久久久久国模美| 午夜免费观看性视频| 欧美国产精品一级二级三级 | 久久人人爽av亚洲精品天堂 | 看免费成人av毛片| 亚洲精品久久午夜乱码| 欧美激情久久久久久爽电影| 丝袜美腿在线中文| 国产老妇女一区| 国精品久久久久久国模美| 波野结衣二区三区在线| 女人十人毛片免费观看3o分钟| 国产成人精品福利久久| 91精品伊人久久大香线蕉| 99精国产麻豆久久婷婷| 欧美日韩国产mv在线观看视频 | 国产精品久久久久久久电影| 久久久色成人| 在线看a的网站| 国产黄色视频一区二区在线观看| 日韩av在线免费看完整版不卡| 国产精品一二三区在线看| 亚洲精品日韩av片在线观看| 韩国高清视频一区二区三区| 亚洲真实伦在线观看| 干丝袜人妻中文字幕| 亚洲欧美清纯卡通| 网址你懂的国产日韩在线| 大片免费播放器 马上看| 美女脱内裤让男人舔精品视频| 午夜老司机福利剧场| 国产精品一区二区性色av| 水蜜桃什么品种好| 国产极品天堂在线| 国产v大片淫在线免费观看| 亚洲伊人久久精品综合| 亚洲欧美日韩另类电影网站 | 亚洲av日韩在线播放| 最近的中文字幕免费完整| 亚洲精品乱码久久久久久按摩| 亚洲精品国产av成人精品| .国产精品久久| 午夜激情久久久久久久| 亚洲在线观看片| 少妇人妻久久综合中文| 日本一本二区三区精品| 九色成人免费人妻av| 欧美日本视频| 精品久久久久久久久av| 国产在线一区二区三区精| 日本与韩国留学比较| 亚洲精品成人av观看孕妇| 日本色播在线视频| 如何舔出高潮| 1000部很黄的大片| 99热6这里只有精品| 亚洲精品亚洲一区二区| 干丝袜人妻中文字幕| 看非洲黑人一级黄片| 新久久久久国产一级毛片| 日韩制服骚丝袜av| 国产欧美日韩精品一区二区| 日韩伦理黄色片| 成年av动漫网址| 亚洲自拍偷在线| 青春草国产在线视频| 久久精品久久久久久噜噜老黄| 久久久久性生活片| 日韩精品有码人妻一区| 极品少妇高潮喷水抽搐| eeuss影院久久| 亚洲人与动物交配视频| 五月伊人婷婷丁香| 最近最新中文字幕大全电影3| 久久久久久伊人网av| 国产在线一区二区三区精| 国产成人精品婷婷| 亚洲无线观看免费| 卡戴珊不雅视频在线播放| 国产av码专区亚洲av| 国内精品宾馆在线| 久久久亚洲精品成人影院| 国产 精品1| 欧美日韩视频精品一区| 蜜桃亚洲精品一区二区三区| 国产精品久久久久久精品电影| 国产一区有黄有色的免费视频| 一级片'在线观看视频| 熟女av电影| 国产免费一区二区三区四区乱码| 天堂中文最新版在线下载 | 国产日韩欧美在线精品| 国产视频首页在线观看| 日日摸夜夜添夜夜添av毛片| 亚洲婷婷狠狠爱综合网|