• <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)用
    99热6这里只有精品| 高清欧美精品videossex| 欧美精品人与动牲交sv欧美| 国产精品一二三区在线看| 一本—道久久a久久精品蜜桃钙片| 黄色一级大片看看| 少妇人妻精品综合一区二区| 国产精品欧美亚洲77777| 日韩欧美一区视频在线观看 | 亚洲色图综合在线观看| 久久精品熟女亚洲av麻豆精品| 毛片一级片免费看久久久久| 伊人久久精品亚洲午夜| 啦啦啦中文免费视频观看日本| 亚洲精品久久久久久婷婷小说| 观看免费一级毛片| 夫妻午夜视频| 男的添女的下面高潮视频| 亚洲精品国产av蜜桃| 欧美xxxx黑人xx丫x性爽| 日韩不卡一区二区三区视频在线| 久久久久久久久久久免费av| 少妇人妻精品综合一区二区| 深爱激情五月婷婷| 精华霜和精华液先用哪个| 九草在线视频观看| 成年人午夜在线观看视频| 日本猛色少妇xxxxx猛交久久| 午夜日本视频在线| 丝瓜视频免费看黄片| 王馨瑶露胸无遮挡在线观看| 最新中文字幕久久久久| 久久精品熟女亚洲av麻豆精品| 18禁裸乳无遮挡免费网站照片| 久久久久久久久大av| 久久久久久久久久久丰满| 观看美女的网站| 亚洲欧洲国产日韩| 精品酒店卫生间| 欧美zozozo另类| a级毛片免费高清观看在线播放| 国产精品福利在线免费观看| 亚洲国产精品999| av卡一久久| 亚洲经典国产精华液单| 久久久精品94久久精品| 看免费成人av毛片| 久久99热6这里只有精品| 免费av中文字幕在线| 亚洲精品乱码久久久久久按摩| 国产午夜精品一二区理论片| 欧美高清性xxxxhd video| av女优亚洲男人天堂| 各种免费的搞黄视频| 午夜福利高清视频| 波野结衣二区三区在线| 97超碰精品成人国产| 国产成人免费无遮挡视频| 国产亚洲5aaaaa淫片| 欧美97在线视频| 纯流量卡能插随身wifi吗| 国产av国产精品国产| 18禁裸乳无遮挡动漫免费视频| 亚洲精品乱码久久久久久按摩| 大陆偷拍与自拍| 赤兔流量卡办理| 91精品国产九色| 国产成人aa在线观看| 国产免费一区二区三区四区乱码| 中国国产av一级| xxx大片免费视频| 久久热精品热| 全区人妻精品视频| 亚洲色图av天堂| 一级毛片aaaaaa免费看小| 久久精品国产鲁丝片午夜精品| 亚洲精品国产av成人精品| 日韩av免费高清视频| 亚洲精品,欧美精品| 三级国产精品欧美在线观看| 亚洲性久久影院| 99热这里只有是精品50| 国产无遮挡羞羞视频在线观看| 春色校园在线视频观看| 啦啦啦视频在线资源免费观看| 国产精品国产三级国产av玫瑰| av国产免费在线观看| 亚洲av中文av极速乱| 国产av码专区亚洲av| 欧美高清成人免费视频www| 下体分泌物呈黄色| 青春草亚洲视频在线观看| 日本色播在线视频| 网址你懂的国产日韩在线| 国产色婷婷99| 天堂俺去俺来也www色官网| 一本久久精品| 国产精品爽爽va在线观看网站| 亚洲国产精品专区欧美| 日日啪夜夜撸| 成年美女黄网站色视频大全免费 | 男女边吃奶边做爰视频| 黑人猛操日本美女一级片| 亚洲高清免费不卡视频| 97精品久久久久久久久久精品| 免费观看性生交大片5| 丰满人妻一区二区三区视频av| 国产精品偷伦视频观看了| 国产精品一区二区性色av| 亚洲av二区三区四区| 一个人看的www免费观看视频| 18+在线观看网站| 国产在线男女| videos熟女内射| 涩涩av久久男人的天堂| 97在线人人人人妻| 蜜桃久久精品国产亚洲av| 精品视频人人做人人爽| 欧美激情极品国产一区二区三区 | 3wmmmm亚洲av在线观看| 亚洲国产精品一区三区| videossex国产| 亚洲精品第二区| 亚洲av成人精品一区久久| 伦理电影大哥的女人| 国产午夜精品一二区理论片| 日韩视频在线欧美| 九九爱精品视频在线观看| 亚洲精品456在线播放app| 啦啦啦中文免费视频观看日本| 亚洲av男天堂| 国内揄拍国产精品人妻在线| 亚洲精品视频女| 亚洲精品视频女| 中文欧美无线码| 国产精品福利在线免费观看| 国产永久视频网站| 欧美精品一区二区大全| 亚洲激情五月婷婷啪啪| 国产一区二区三区av在线| 欧美三级亚洲精品| 国产精品.久久久| 国产探花极品一区二区| 免费看不卡的av| 99re6热这里在线精品视频| 久久久久久久亚洲中文字幕| 日韩一区二区视频免费看| 女的被弄到高潮叫床怎么办| 国产久久久一区二区三区| 美女视频免费永久观看网站| 亚洲国产成人一精品久久久| 啦啦啦啦在线视频资源| 婷婷色av中文字幕| 亚洲国产成人一精品久久久| 高清毛片免费看| 极品少妇高潮喷水抽搐| 久久久a久久爽久久v久久| 国产久久久一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 97在线人人人人妻| 成人二区视频| 热re99久久精品国产66热6| 欧美极品一区二区三区四区| 久久99热6这里只有精品| 麻豆精品久久久久久蜜桃| 大香蕉97超碰在线| 日韩免费高清中文字幕av| a 毛片基地| 免费黄网站久久成人精品| 蜜桃久久精品国产亚洲av| 亚洲av在线观看美女高潮| 极品少妇高潮喷水抽搐| 99久久精品一区二区三区| 久久久精品免费免费高清| 赤兔流量卡办理| kizo精华| 亚洲欧美日韩东京热| 黑丝袜美女国产一区| 寂寞人妻少妇视频99o| 人妻夜夜爽99麻豆av| 免费观看a级毛片全部| 七月丁香在线播放| 久久亚洲国产成人精品v| 亚洲精品乱码久久久v下载方式| 下体分泌物呈黄色| 男女免费视频国产| 成年av动漫网址| 国语对白做爰xxxⅹ性视频网站| 老司机影院成人| 一本一本综合久久| 天堂俺去俺来也www色官网| 免费黄色在线免费观看| 纯流量卡能插随身wifi吗| 九九爱精品视频在线观看| 伊人久久国产一区二区| 国产日韩欧美亚洲二区| 日韩视频在线欧美| 18+在线观看网站| 日韩成人伦理影院| 中国美白少妇内射xxxbb| 欧美人与善性xxx| 女性生殖器流出的白浆| 亚洲经典国产精华液单| 色综合色国产| 亚洲色图av天堂| 国产在线视频一区二区| 亚洲在久久综合| 国产欧美亚洲国产| 欧美人与善性xxx| 一个人看的www免费观看视频| 成年美女黄网站色视频大全免费 | 欧美日韩视频精品一区| 日本免费在线观看一区| 婷婷色麻豆天堂久久| 亚洲欧美精品专区久久| 欧美日韩视频精品一区| 欧美bdsm另类| 国产精品无大码| 亚洲av成人精品一二三区| 少妇的逼好多水| 自拍欧美九色日韩亚洲蝌蚪91 | av天堂中文字幕网| 91久久精品电影网| 亚洲欧美日韩东京热| av.在线天堂| 精品熟女少妇av免费看| 久久久午夜欧美精品| 亚洲熟女精品中文字幕| 国产成人免费观看mmmm| 噜噜噜噜噜久久久久久91| 视频中文字幕在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲av综合色区一区| 18禁动态无遮挡网站| 亚洲第一av免费看| 国产又色又爽无遮挡免| 国产成人freesex在线| 国产黄片美女视频| 男的添女的下面高潮视频| 少妇人妻一区二区三区视频| 一级a做视频免费观看| 校园人妻丝袜中文字幕| 高清不卡的av网站| 91在线精品国自产拍蜜月| 国内精品宾馆在线| 久久久久精品性色| 极品教师在线视频| 久久国产精品男人的天堂亚洲 | 国产人妻一区二区三区在| 老司机影院成人| 久久精品久久久久久久性| 99久久人妻综合| 欧美+日韩+精品| 国产精品久久久久久久电影| 亚洲精品乱码久久久久久按摩| 免费观看a级毛片全部| 精品一区在线观看国产| 黑人猛操日本美女一级片| 少妇裸体淫交视频免费看高清| 在线观看免费高清a一片| 人妻系列 视频| 中文字幕制服av| 国产又色又爽无遮挡免| 美女脱内裤让男人舔精品视频| 男女边摸边吃奶| 99久久精品热视频| 国产精品爽爽va在线观看网站| av一本久久久久| 舔av片在线| 日本一二三区视频观看| 午夜激情福利司机影院| 亚洲人与动物交配视频| 熟妇人妻不卡中文字幕| 国产在线男女| 欧美最新免费一区二区三区| 日本av免费视频播放| 日本wwww免费看| 亚洲精品,欧美精品| 久久99热这里只有精品18| 久久国产乱子免费精品| 日韩av免费高清视频| 美女内射精品一级片tv| 少妇 在线观看| 免费看av在线观看网站| 日本欧美视频一区| 成人午夜精彩视频在线观看| 亚洲欧美日韩东京热| 日韩欧美 国产精品| 成人黄色视频免费在线看| 欧美国产精品一级二级三级 | 亚洲精品色激情综合| 少妇裸体淫交视频免费看高清| 五月天丁香电影| 如何舔出高潮| 国产淫语在线视频| 最近中文字幕2019免费版| 亚洲欧洲国产日韩| 国产亚洲最大av| 久久久久精品性色| 国产免费又黄又爽又色| 在线播放无遮挡| 全区人妻精品视频| 街头女战士在线观看网站| 日韩三级伦理在线观看| 久久99热这里只频精品6学生| 亚洲精品456在线播放app| 国产高清三级在线| 亚洲,欧美,日韩| 我的老师免费观看完整版| 久久精品人妻少妇| 日韩精品有码人妻一区| 狂野欧美激情性xxxx在线观看| 亚洲va在线va天堂va国产| 男人狂女人下面高潮的视频| 亚洲精品一区蜜桃| 大又大粗又爽又黄少妇毛片口| 久久久久久久久久人人人人人人| 亚洲美女视频黄频| 性色avwww在线观看| 久久久久久久亚洲中文字幕| 国产成人精品久久久久久| 亚洲精品乱码久久久久久按摩| 国产精品av视频在线免费观看| 国产精品爽爽va在线观看网站| 2021少妇久久久久久久久久久| 国产精品久久久久久精品电影小说 | 久久ye,这里只有精品| 久久6这里有精品| 熟妇人妻不卡中文字幕| 波野结衣二区三区在线| 亚洲综合精品二区| 国产乱来视频区| 国产精品国产av在线观看| 婷婷色综合大香蕉| 日韩欧美精品免费久久| 亚洲精品乱码久久久v下载方式| 网址你懂的国产日韩在线| 高清视频免费观看一区二区| 一个人看的www免费观看视频| 亚洲精品视频女| 色吧在线观看| 精品久久久久久久末码| 成年免费大片在线观看| 国产精品一二三区在线看| 国产v大片淫在线免费观看| 夜夜爽夜夜爽视频| 成人亚洲精品一区在线观看 | 日日啪夜夜撸| av免费在线看不卡| 搡老乐熟女国产| av在线播放精品| 久久久精品94久久精品| 99久国产av精品国产电影| 五月天丁香电影| 男人添女人高潮全过程视频| 在线天堂最新版资源| xxx大片免费视频| 五月天丁香电影| 久久这里有精品视频免费| 寂寞人妻少妇视频99o| 国产成人免费观看mmmm| 亚洲欧美日韩卡通动漫| 狂野欧美白嫩少妇大欣赏| 在线播放无遮挡| 哪个播放器可以免费观看大片| 亚洲va在线va天堂va国产| av线在线观看网站| 一个人看的www免费观看视频| 久久久久久伊人网av| 日本欧美国产在线视频| videos熟女内射| 亚洲成色77777| 国产成人免费观看mmmm| 欧美xxxx黑人xx丫x性爽| 高清毛片免费看| 国产成人精品久久久久久| 免费不卡的大黄色大毛片视频在线观看| 亚洲综合精品二区| 老女人水多毛片| 亚洲精品aⅴ在线观看| 国产精品人妻久久久影院| 少妇人妻精品综合一区二区| 日韩,欧美,国产一区二区三区| 国产精品无大码| 王馨瑶露胸无遮挡在线观看| 免费少妇av软件| 欧美精品一区二区大全| 国产高清国产精品国产三级 | 国产在线免费精品| 免费观看a级毛片全部| 男男h啪啪无遮挡| 国产精品久久久久久久久免| 不卡视频在线观看欧美| 国产精品秋霞免费鲁丝片| 性色avwww在线观看| 一个人免费看片子| 22中文网久久字幕| 欧美国产精品一级二级三级 | 国产午夜精品久久久久久一区二区三区| 日本wwww免费看| 美女cb高潮喷水在线观看| 成人美女网站在线观看视频| 制服丝袜香蕉在线| 啦啦啦在线观看免费高清www| 99久国产av精品国产电影| 亚洲人成网站高清观看| 天堂中文最新版在线下载| 最新中文字幕久久久久| 亚洲成人手机| 老女人水多毛片| 成人毛片a级毛片在线播放| 97在线视频观看| 国产中年淑女户外野战色| 十八禁网站网址无遮挡 | 国产精品国产三级专区第一集| av一本久久久久| 国产老妇伦熟女老妇高清| 小蜜桃在线观看免费完整版高清| 国产黄频视频在线观看| 色视频在线一区二区三区| 欧美极品一区二区三区四区| 全区人妻精品视频| 有码 亚洲区| 精品99又大又爽又粗少妇毛片| 久久国产精品大桥未久av | 国产极品天堂在线| 久久99热这里只频精品6学生| 久久精品国产亚洲av天美| 纵有疾风起免费观看全集完整版| 黑人猛操日本美女一级片| 午夜精品国产一区二区电影| 国产一区二区三区综合在线观看 | av免费观看日本| 一区二区三区乱码不卡18| 久久亚洲国产成人精品v| 国产精品一及| 欧美精品亚洲一区二区| 久久毛片免费看一区二区三区| 国产精品av视频在线免费观看| 日日摸夜夜添夜夜爱| 如何舔出高潮| 久久精品久久精品一区二区三区| 18禁裸乳无遮挡免费网站照片| 婷婷色av中文字幕| av视频免费观看在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产精品一二三区在线看| 国产高清三级在线| 国产探花极品一区二区| 成人国产av品久久久| 欧美变态另类bdsm刘玥| 中文资源天堂在线| 亚洲av欧美aⅴ国产| 国产精品偷伦视频观看了| av免费在线看不卡| 久久精品久久精品一区二区三区| 一本一本综合久久| 99热这里只有是精品50| 日韩人妻高清精品专区| 最近中文字幕2019免费版| 亚洲av综合色区一区| 亚洲,欧美,日韩| 国内揄拍国产精品人妻在线| 王馨瑶露胸无遮挡在线观看| 国产日韩欧美亚洲二区| 国产精品国产三级专区第一集| 午夜免费鲁丝| 亚洲精品国产av蜜桃| 亚洲真实伦在线观看| 久久韩国三级中文字幕| 国产精品久久久久成人av| 伊人久久精品亚洲午夜| 国产精品国产三级专区第一集| 男女边吃奶边做爰视频| 在线免费十八禁| 国产色婷婷99| 亚洲激情五月婷婷啪啪| 亚洲av福利一区| 日本欧美国产在线视频| 国产成人午夜福利电影在线观看| 欧美日韩在线观看h| av卡一久久| 在线精品无人区一区二区三 | 国产 一区精品| 一级毛片我不卡| 我要看黄色一级片免费的| 午夜激情福利司机影院| 久久久欧美国产精品| 永久免费av网站大全| 国产欧美另类精品又又久久亚洲欧美| 成人无遮挡网站| 成年免费大片在线观看| 中国美白少妇内射xxxbb| 欧美精品人与动牲交sv欧美| 午夜福利网站1000一区二区三区| 3wmmmm亚洲av在线观看| 一区二区三区四区激情视频| 婷婷色av中文字幕| 寂寞人妻少妇视频99o| 97在线视频观看| 国产精品一区www在线观看| 少妇精品久久久久久久| 亚洲成人av在线免费| 国产精品一及| 香蕉精品网在线| 免费观看av网站的网址| 国产精品三级大全| 国产 一区 欧美 日韩| 在线天堂最新版资源| 精品午夜福利在线看| 免费观看的影片在线观看| 在线观看一区二区三区激情| 插逼视频在线观看| 日日撸夜夜添| 韩国av在线不卡| 国产精品不卡视频一区二区| 国产精品福利在线免费观看| 91在线精品国自产拍蜜月| 国产91av在线免费观看| 国产精品.久久久| 一级毛片电影观看| 久久久久性生活片| 欧美人与善性xxx| 国产一区二区三区av在线| 大片免费播放器 马上看| 国产久久久一区二区三区| 在线精品无人区一区二区三 | 亚洲av电影在线观看一区二区三区| 欧美3d第一页| 人人妻人人爽人人添夜夜欢视频 | 国产精品国产三级国产av玫瑰| 91精品伊人久久大香线蕉| 老熟女久久久| 久久久久久久久久人人人人人人| 亚洲国产av新网站| 国产综合精华液| 99久久精品一区二区三区| 国产真实伦视频高清在线观看| 校园人妻丝袜中文字幕| 一级av片app| 国产伦精品一区二区三区四那| 欧美激情国产日韩精品一区| 亚洲国产高清在线一区二区三| 午夜福利在线在线| 久久国产乱子免费精品| 国产精品人妻久久久影院| 久久综合国产亚洲精品| 天堂俺去俺来也www色官网| 亚洲图色成人| 寂寞人妻少妇视频99o| 激情五月婷婷亚洲| 在线亚洲精品国产二区图片欧美 | 蜜臀久久99精品久久宅男| www.色视频.com| 久久久久网色| 国产成人午夜福利电影在线观看| 国产黄频视频在线观看| 最近中文字幕高清免费大全6| 亚洲精品乱码久久久v下载方式| 伦精品一区二区三区| 91午夜精品亚洲一区二区三区| 18禁在线播放成人免费| 国产欧美日韩一区二区三区在线 | 亚洲经典国产精华液单| 久久久久视频综合| 韩国高清视频一区二区三区| 黄色日韩在线| 搡老乐熟女国产| kizo精华| 99热国产这里只有精品6| 狂野欧美白嫩少妇大欣赏| 国产在线男女| 成人特级av手机在线观看| 干丝袜人妻中文字幕| 大话2 男鬼变身卡| 精品一区二区三区视频在线| 国产精品99久久久久久久久| 精华霜和精华液先用哪个| 丰满少妇做爰视频| 又爽又黄a免费视频| 日本vs欧美在线观看视频 | 欧美精品亚洲一区二区| 三级经典国产精品| 午夜老司机福利剧场| 国产男女超爽视频在线观看| 青春草亚洲视频在线观看| 欧美高清成人免费视频www| 美女xxoo啪啪120秒动态图| 精品亚洲成a人片在线观看 | 国产av国产精品国产| 精品亚洲成a人片在线观看 | 亚洲aⅴ乱码一区二区在线播放| 久久鲁丝午夜福利片| av播播在线观看一区| 免费观看无遮挡的男女| 日韩电影二区| 亚洲精品,欧美精品| 欧美3d第一页| 男女边摸边吃奶| 青春草国产在线视频| 永久免费av网站大全| 成年人午夜在线观看视频| 黄色一级大片看看| 狂野欧美激情性xxxx在线观看| 日韩欧美一区视频在线观看 | 人人妻人人澡人人爽人人夜夜| 久久ye,这里只有精品| 蜜桃在线观看..| 少妇的逼好多水| 小蜜桃在线观看免费完整版高清| 我的女老师完整版在线观看| 一个人免费看片子| 最新中文字幕久久久久| 亚洲av在线观看美女高潮| 久久热精品热|