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

    彈體超高速侵徹石灰?guī)r靶體地沖擊的數(shù)值模擬研究*

    2022-02-11 09:35:40張山豹孔祥振
    爆炸與沖擊 2022年1期
    關(guān)鍵詞:靶體超高速石灰?guī)r

    張山豹,孔祥振,2,方 秦,洪 建

    (1. 陸軍工程大學(xué)爆炸沖擊防災(zāi)減災(zāi)國家重點實驗室,江蘇 南京 210007;2. 東南大學(xué)爆炸安全防護教育部工程研究中心,江蘇 南京 211189)

    近年來,以“上帝之杖”為代表的超高速動能武器迅速發(fā)展,其終點毀傷效應(yīng)成為武器研發(fā)及工程防護領(lǐng)域關(guān)注的熱點問題。彈體對地打擊的破壞模式主要表現(xiàn)為兩個方面:(1) 彈體沖擊靶體,侵入靶體直接造成破壞;(2) 侵徹沖擊產(chǎn)生的瞬時高壓以應(yīng)力波的形式在靶體中傳播,造成靶體內(nèi)部的間接破壞效應(yīng)(如開裂,震塌等破壞現(xiàn)象),即為侵徹的地沖擊效應(yīng)。已有實驗表明,彈體侵徹混凝土和巖石等脆性靶體時會呈現(xiàn)剛性侵徹、變形不侵蝕侵徹和侵蝕侵徹3 個典型階段,且會出現(xiàn)侵徹深度隨初速度增加而減小的情況。彈體超高速侵徹的直接侵徹深度相較于中低速侵徹沒有明顯增加,地沖擊效應(yīng)可能成為其對地破壞的主要途徑。因此,研究彈體超高速侵徹的地沖擊規(guī)律,對綜合評價超高速動能武器的對地毀傷效應(yīng)具有重要意義。

    “上帝之杖”動能彈的入地侵徹速度可超過3 km/s,但受限于實驗條件,對彈體在該速度區(qū)間侵徹巖石類脆性材料的研究較少。Poelchau 等、李干等、錢秉文等分別開展了不同彈體超高速侵徹巖石類脆性材料的實驗研究,主要針對靶體開坑形態(tài)、侵徹深度的變化規(guī)律等進行研究??紫檎竦葘楏w超高速侵徹半無限金屬靶的6 種典型的理論模型進行了對比分析,發(fā)現(xiàn)其均不適用于混凝土等脆性靶體的侵徹深度預(yù)測。特別是針對彈體超高速侵徹巖石類材料的地沖擊研究,目前還較匱乏且主要以數(shù)值模擬方法為主?;贕EODYN 流體動力學(xué)計算軟件和自定義的材料動態(tài)本構(gòu)關(guān)系,Antoun 等開展了鎢合金彈體超高速侵徹石灰?guī)r靶體的數(shù)值模擬研究,重點探討了初速度為4 和6 km/s 情況下的地沖擊效應(yīng)。鄧國強等以流體彈塑性模型為基礎(chǔ),采用無網(wǎng)格的光滑粒子流體動力學(xué)(smooth particle hydrodynamics, SPH)方法,對彈體超高速侵徹花崗巖靶體進行了數(shù)值模擬研究,分析了超高速打擊的破壞現(xiàn)象及靶體中應(yīng)力波的傳播規(guī)律,并基于破壞特征和應(yīng)力波傳播特征提出了超高速侵徹和裝藥淺埋爆炸的等效方法?;诔咚偾謴爻煽优c淺埋爆炸成坑的對比分析,王明洋等提出了超高速侵徹地沖擊的等效計算方法,并結(jié)合開展的合金鋼彈體超高速侵徹花崗巖靶體實驗對等效方法進行了驗證。可以發(fā)現(xiàn),不同學(xué)者選用的彈靶材料、本構(gòu)模型參數(shù)以及數(shù)值模擬方法均不相同,且對地沖擊的影響因素尚未有系統(tǒng)的分析。

    本文中,以鎢合金長桿彈和石灰?guī)r靶體為研究對象,對超高速侵徹的地沖擊進行數(shù)值模擬。由于超高速侵徹時石灰?guī)r靶體處于高圍壓、高應(yīng)變率狀態(tài),因此為保證數(shù)值模擬的可靠性,首先基于石灰?guī)r材料的靜動態(tài)力學(xué)實驗數(shù)據(jù)對材料模型參數(shù)進行詳細標定,并結(jié)合彈體大范圍著速侵徹石灰?guī)r靶體實驗驗證模型參數(shù)的準確性;然后基于上述驗證的數(shù)值模型,對彈體超高速侵徹石灰?guī)r靶體進行數(shù)值模擬,研究不同初速度下的地沖擊規(guī)律,并進一步詳細探討彈靶參數(shù)對地沖擊的影響規(guī)律。

    1 石灰?guī)r材料模型參數(shù)的標定

    石灰?guī)r是一種多孔、中等強度的巖石材料,抗壓強度約為60 MPa,初始密度和壓實密度分別為2 280 和2 700 kg/m3。已有石灰?guī)r的靜動態(tài)力學(xué)性能實驗數(shù)據(jù)較多,因此靶體選用石灰?guī)r,數(shù)值模擬中材料模型選用前期自主開發(fā)的Kong-Fang 模型。本節(jié)中首先簡要介紹Kong-Fang 材料模型,然后基于已有的靜動態(tài)力學(xué)性能實驗數(shù)據(jù)對模型參數(shù)進行標定。

    1.1 Kong-Fang 模型

    在對已有脆性材料模型(HJC、RHT 和K&C 模型等)進行精細化評估的基礎(chǔ)上,Kong 等提出了適用于爆炸沖擊作用下的Kong-Fang 材料模型,并將其植入LS-DYNA 軟件中。Zhang 等結(jié)合高速侵徹下脆性靶體材料的高圍壓屈服和動態(tài)拉伸特征,對Kong-Fang 模型進行了改進,并將改進的Kong-Fang 模型植入AUTODYN 商業(yè)軟件中,經(jīng)驗證能夠較好地模擬脆性材料靶體的損傷破壞效應(yīng)。本文中,模型介紹及數(shù)值模擬均基于改進后的Kong-Fang 模型。

    Kong-Fang 模型的強度面包含最大強度面σ、屈服強度面σ和殘余強度面σ,并通過損傷變量和η 插值得到當前破壞強度面σ :

    式中:為當前子午面與壓縮子午面的比值;為壓力,受壓為正;λ 為修正的等效塑性應(yīng)變;λ為峰值應(yīng)力對應(yīng)的λ。

    Kong-Fang 模型能夠區(qū)分拉壓損傷,總損傷由拉壓損傷耦合而成:

    式中:為壓縮損傷,為應(yīng)力偏量造成的拉伸損傷,Δ為應(yīng)力球量造成的拉伸損傷。強度面及損傷相關(guān)的參數(shù)標定詳見1.2 節(jié)。

    Kong-Fang 模型使用徑向放大破壞強度面的方法考慮應(yīng)變率效應(yīng),最終強度面按照下式確定:

    式中:為材料的動態(tài)強度增強因子,在Kong-Fang 模型中采用Xu 等提出的經(jīng)驗公式,具體形式及參數(shù)標定詳見1.3 節(jié)。

    1.2 強度面參數(shù)的標定

    Kong-Fang 模型中最大強度面σ需通過圍壓壓縮實驗數(shù)據(jù)標定,屈服強度面σ和殘余強度面σ分別通過材料的應(yīng)變硬化和應(yīng)變軟化曲線特征進行標定,各強度面具體形式為:

    式中:、和分別為材料的單軸靜態(tài)拉伸強度、壓縮強度和屈服強度。、、、和為強度面中的經(jīng)驗參數(shù),需根據(jù)實驗數(shù)據(jù)進行標定;ψ 為拉壓子午面之比,需要根據(jù)圍壓壓縮(σ<σ=σ,σ、σ、σ分別為第一、第二和第三主應(yīng)力,以拉為正)和圍壓拉伸(σ>σ=σ)兩種狀態(tài)的強度面關(guān)系確定。

    利用Fossum 等開展的三軸實驗數(shù)據(jù)對原Kong-Fang 模型的拉壓子午面比ψ 標定為:

    根據(jù)Fossum 等和Walton 等開展的不同圍壓下的石灰?guī)r圍壓壓縮實驗數(shù)據(jù),擬合得到最大強度面曲線,如圖1 所示。同時考慮巖石材料的剪切飽和特性,限定強度面最大值為20。

    圖1 最大強度面模型與實驗數(shù)據(jù)的擬合Fig. 1 Fitting of the maximum strength surface model to the experimental data

    根據(jù)石灰?guī)r壓縮實驗的應(yīng)力-應(yīng)變曲線,取屈服強度=0.8,并對屈服強度面、殘余強度面及損傷相關(guān)的參數(shù)(見表1)進行擬合。基于標定的參數(shù)開展單軸壓縮(UUC)及圍壓壓縮(TXC)數(shù)值模擬,其應(yīng)力-應(yīng)變曲線與實驗結(jié)果的對比如圖2 所示。可以看出,Kong-Fang 模型能夠較好地模擬出石灰?guī)r材料的應(yīng)變硬化和應(yīng)變軟化特征,特別是石灰?guī)r低圍壓下呈現(xiàn)脆性破壞而高圍壓下呈現(xiàn)塑性破壞這一特性。

    表1 石灰?guī)r的強度模型參數(shù)Table 1 Parameters of the strength surface models for limestone

    圖2 石灰?guī)r應(yīng)力-應(yīng)變曲線的數(shù)值模擬結(jié)果與實驗數(shù)據(jù)的對比Fig. 2 Comparison of stress-strain curves of limestone between experimental data and numerical simulation

    1.3 應(yīng)變率效應(yīng)的參數(shù)標定

    Kong-Fang 模型采用Xu 等提出的應(yīng)變率增強因子計算公式:

    式中:、、和為相關(guān)經(jīng)驗參數(shù)。利用已有的實驗數(shù)據(jù)對石灰?guī)r的壓縮動態(tài)強度增強因子(dynamic increase factor of compression, DIFC)和拉伸動態(tài)強度增強因子(dynamic increase factor of tension, DIFT)分別進行標定,其參數(shù)如表2 所示,應(yīng)變率增強因子的曲線與相關(guān)實驗結(jié)果的對比如圖3 所示。

    圖3 石灰?guī)r動態(tài)強度增強因子隨應(yīng)變率的變化曲線Fig. 3 Changes of the dynamic increase factors with strain rate for limestone

    表2 石灰?guī)r的動態(tài)強度增強因子參數(shù)Table 2 Parameters for dynamic increase factors of limestone

    1.4 狀態(tài)方程參數(shù)的標定

    -α 狀態(tài)方程是由Herrmann提出的一種多孔材料的狀態(tài)方程(equation of state, EOS),能夠較好地描述孔隙壓實過程中的材料性質(zhì)變化,因此本文中的石灰?guī)r采用該狀態(tài)方程。-α 狀態(tài)方程根據(jù)材料孔隙壓實過程分為3 個階段。第一階段為壓力小于孔隙坍塌壓力時,采用線彈性狀態(tài)方程:

    式中:為材料的體積模量,為體積應(yīng)變。

    第二階段為孔隙壓實的過程,階段的起始分別由孔隙坍塌壓力和完全壓實壓力控制:

    式中:α 為材料的孔隙率,α為對應(yīng)于的材料孔隙率;()為材料在第三階段密實段所采用的方程;為壓縮系數(shù)。第三階段材料完全密實,可采用多項式狀態(tài)方程描述:

    式中:、、為壓縮段參數(shù),一般取=;為拉伸段參數(shù),拉伸時按照線彈性考慮,即=。

    對于線彈性段和孔隙壓實段,可根據(jù)石灰?guī)r的靜水壓力實驗結(jié)果標定,而對于高壓密實段主要通過平板撞擊實驗數(shù)據(jù)進行參數(shù)標定,擬合得到的狀態(tài)方程參數(shù)及曲線分別如表3 和圖4 所示。

    表3 石灰?guī)r的狀態(tài)方程參數(shù)Table 3 Equation of state parameters for limestone

    圖4 石灰?guī)r狀態(tài)方程曲線與實驗數(shù)據(jù)的擬合Fig. 4 Equation of state of limestone fitted to experimental data

    1.5 模型驗證

    為了驗證上述材料模型及參數(shù)的準確性,對Frew 等開展的4340 鋼彈體大范圍打擊速度侵徹石灰?guī)r靶體的實驗進行數(shù)值模擬,并對結(jié)果進行比較。實驗中彈體直徑分別為7.1 和12.7 mm,長徑比均為10,靶體尺寸根據(jù)彈體初始速度的不同而變化?;谇捌诮?jīng)驗,數(shù)值模擬采用有限元與SPH 相結(jié)合的方法,其中彈體采用有限元,靶體采用SPH 粒子,并利用二維軸對稱算法以提高計算效率。實驗后彈體變形較小,因此視為剛性材料,石灰?guī)r靶體由Kong-Fang 模型描述且參數(shù)按1.4 節(jié)中標定的選取。

    數(shù)值模擬得到的不同直徑彈體侵徹石灰?guī)r靶體的侵徹深度()與實驗數(shù)據(jù)的對比如圖5 所示,圖中展示了彈體直徑分別為7.1 和12.7 mm 工況的結(jié)果對比。可以看出,數(shù)值模擬得到的侵徹深度與實驗數(shù)據(jù)較接近,驗證了石灰?guī)r材料模型及其參數(shù)的可靠性。

    圖5 數(shù)值模擬的侵徹深度與實驗數(shù)據(jù)的對比Fig. 5 Simulated depths of penetration at different initial projectile velocities compared with experimental data

    2 彈體超高速侵徹石灰?guī)r靶體的地沖擊

    當彈體撞擊到靶體表面時,彈靶交界面處會產(chǎn)生瞬時高壓力。一方面,彈體繼續(xù)向前運動,隧道區(qū)逐漸擴展,侵徹深度不斷增大,且彈體的速度也逐漸降低,并可能發(fā)生變形、侵蝕等破壞現(xiàn)象;另一方面,在彈靶交界處產(chǎn)生的瞬時高壓會在靶體中以應(yīng)力波的形式傳播,并對內(nèi)部靶體產(chǎn)生不同程度的影響,如圖6 所示。應(yīng)力波在靶體中衰減、反射并發(fā)生相互作用,會造成石灰?guī)r材料的壓實、破碎或斷裂,并可能進一步在靶體中形成開裂、開坑及震塌等破壞現(xiàn)象。

    圖6 應(yīng)力波造成的靶體破壞分區(qū)Fig. 6 Damage regions of the target caused by stress waves

    當彈體侵徹速度較低時,彈體對靶體的破壞主要表現(xiàn)為直接侵徹(隧道區(qū));而當彈體侵徹速度較高時,彈體會發(fā)生較嚴重的磨蝕甚至完全侵蝕,直接侵徹隧道區(qū)深度不會顯著增大,而較高的侵徹速度會導(dǎo)致彈靶交界面處的壓力升高,因此應(yīng)力波對靶體的破壞效應(yīng)凸顯?;谏鲜鲵炞C的數(shù)值模擬方法及材料模型參數(shù),本節(jié)對鎢合金彈體超高速侵徹石灰?guī)r靶體的地沖擊進行研究。

    2.1 數(shù)值模型

    采用AUTODYN-2D 軸對稱算法,其中圓柱形鎢合金彈體長2 m,直徑為0.2 m,初速度區(qū)間考慮為1.5~5.0 km/s,采用有限元方法建模;石灰?guī)r靶體長20 m,直徑為10 m,采用SPH 方法建模。彈體有限元網(wǎng)格大小按前期經(jīng)驗確定為10 mm×10 mm,同時對靶體的粒子間距大小進行了敏感性分析,分別開展了靶體粒子間距為50、40 和30 mm 的數(shù)值模擬,發(fā)現(xiàn)數(shù)值預(yù)測的侵徹深度隨粒子間距減小趨于穩(wěn)定。因此,同時考慮計算的準確性和計算效率,粒子間距確定為40 mm。在彈體著靶點正下方6、7、8、9、10、12、14、16、18 和20 m 處的靶體內(nèi)部分別布置壓力測點。

    鎢合金彈體采用Shock 狀態(tài)方程結(jié)合Steinberg 強度模型,能夠較好地體現(xiàn)鎢合金材料在高速侵徹下的高應(yīng)變率高壓力學(xué)性能,參數(shù)取AUTODYN 中的默認值。為了較好地描述彈體的變形及侵蝕,采用不同的單元刪除算法對初速度為5.0 km/s 的工況進行了收斂性論證。為了盡可能減少單元刪除的影響,且綜合考慮彈體的破壞模式、計算效率等因素,采用基于應(yīng)變的單元刪除算法,并確定單元刪除閾值為0.8。石灰?guī)r靶體采用1.1 節(jié)中介紹的Kong-Fang 模型及所標定的材料參數(shù)。

    從模擬得到的彈體破壞程度來看,當彈體初速度低于2.5 km/s 時,侵徹完成時彈體尚有殘余長度;而當彈體速度高于3.0 km/s 時,在侵徹過程中彈體被完全侵蝕,之后隧道區(qū)底部的靶體材料由于慣性會繼續(xù)運動,造成侵徹深度進一步增大。不同撞擊速度下模擬得到的最終侵徹深度如圖7 所示,其中黑色虛線為根據(jù)模擬結(jié)果擬合得到的趨勢線??梢钥闯?,在研究的速度區(qū)間內(nèi),侵徹可大致分為侵蝕侵徹階段和流體侵徹階段,侵徹深度隨彈體初速度的變化整體呈先減小再增大的趨勢,并最終趨于穩(wěn)定,與前期的研究結(jié)論相符。并且最終的穩(wěn)定侵徹深度大于流體動力學(xué)極限,這是由靶體材料的慣性運動造成的。

    圖7 數(shù)值模擬得到不同撞擊速度下的的侵徹深度Fig. 7 Simulated depths of penetration at different initial projectile velocities

    2.2 超高速侵徹的應(yīng)力波傳播

    以彈體初速度為5.0 km/s 的工況為例,對彈體超高速侵徹石灰?guī)r靶體的過程進行分析,其破壞現(xiàn)象和壓力波傳播過程如圖8 所示。彈體撞擊靶體時,彈靶交界面產(chǎn)生瞬時壓力(高達30 GPa)并開始向周圍傳播,侵徹過程中彈靶接觸面的高壓力狀態(tài)會持續(xù)以應(yīng)力波的形式傳播,如圖8(a)所示;1.4 ms 時,彈體已被完全磨蝕,直接侵徹深度約為4.45 m,如圖8(b)所示;彈體完全磨蝕后,坑底靶體材料在慣性作用下繼續(xù)運動,彈坑繼續(xù)擴張,在3.0 ms 時,侵徹深度約為4.81 m;最終計算到100.0 ms 時,侵徹深度為6.20 m且趨于穩(wěn)定。從圖8(c)~(d)可以看出,壓力波基本呈半球形向靶體內(nèi)傳播,與地面淺埋爆炸具有一定的相似性;但由于侵徹過程具有方向性,沿深度方向的壓力波峰值大于水平方向同等距離的壓力波峰值。圖8(d)中,當主要的地沖擊已經(jīng)傳播到靶體內(nèi)部較深處時,靶體靠近彈坑處仍然有壓力分布,這是由于彈靶之間的相互作用會持續(xù)在彈靶交界面處產(chǎn)生新的高壓狀態(tài),并持續(xù)向靶體中傳播。

    圖8 超高速侵徹過程中的破壞現(xiàn)象及壓力波傳播Fig. 8 Damage and pressure wave propagation in the target during hypervelocity penetration

    其他彈體初速度下的壓力波傳播過程與上述算例基本一致。超高速侵徹造成的地沖擊:一方面,會在隧道區(qū)周圍造成高壓區(qū)域,導(dǎo)致材料發(fā)生壓實破碎甚至相變;另一方面,可能傳入地下防護結(jié)構(gòu)頂板,在頂板自由表面反射產(chǎn)生拉伸波,造成開裂、震塌等破壞效應(yīng)。

    2.3 超高速侵徹的地沖擊分析

    圖9 給出了不同侵徹速度下靶體中各測點沿深度方向的徑向應(yīng)力時程變化??梢钥闯霾煌謴厮俣认碌膽?yīng)力波波形存在一定差異:當侵徹速度較低時,其升壓時間較長,且可能出現(xiàn)壓力進入平臺段或開始下降、后再次升高的情況,這是由于侵徹完成時彈體尚有殘留,會持續(xù)在彈靶接觸面產(chǎn)生壓力并向靶體內(nèi)部傳播;而當侵徹速度較高時,其應(yīng)力迅速達到峰值并降低,這可能是由于彈體磨蝕較快,地沖擊主要是由于初始侵徹階段引起的瞬時高壓造成的。當侵徹速度較高時,其升壓時間較短,特別是彈體初速度為5.0 km/s 時,其升壓時間均在0.5 ms 以內(nèi)。需要說明的是,在彈體侵徹過程中,彈靶接觸面上的應(yīng)力峰值較大。對于布置深度小于或略大于侵徹深度的壓力測點,會明顯受到彈體的干擾,因此在總結(jié)地沖擊規(guī)律時只選取了地沖擊傳播穩(wěn)定后的測點數(shù)據(jù),且不同初速度下的直接侵徹深度不同,因而地沖擊傳播穩(wěn)定后的測點數(shù)量也不同。

    圖9 不同深度的應(yīng)力時程曲線Fig. 9 Time histories of stress waves at different depths

    圖10 給出了不同侵徹速度下應(yīng)力峰值隨深度變化的曲線,可以看出:隨著深度的增大,應(yīng)力峰值均不斷降低,應(yīng)力波在傳播過程中不斷衰減;隨著侵徹速度的升高,同一深度處的應(yīng)力峰值也會隨之升高;當初速度高于3.0 km/s 時,不同深度處的應(yīng)力峰值升高較快??傮w來看,彈體初速度較低時,侵徹產(chǎn)生的地沖擊較小,對靶體內(nèi)部的影響也較?。欢敵跛俣戎饾u升高時,地沖擊的影響范圍逐漸加深。

    圖10 不同彈體初速度下應(yīng)力峰值隨深度的變化趨勢Fig. 10 Change of peak stress with depth at different initial projectile velocities

    3 地沖擊的影響因素

    對于中、低速侵徹,彈靶材料、尺寸等因素對最終的破壞效應(yīng)有顯著的影響。對于超高速侵徹,已有研究表明:彈靶的強度不再起主導(dǎo)作用,破壞效應(yīng)主要與彈體長徑比(/)、彈靶密度比等因素有關(guān)。而對應(yīng)力波的傳播,靶體材料的孔隙率會有較大影響。因此,本節(jié)以彈體初速度5.0 km/s 工況為基準,探究彈體長徑比、彈靶密度比和靶體孔隙率等因素對地沖擊的影響。

    3.1 彈體長徑比

    彈體長徑比,特別是彈體長度是影響侵徹深度的重要因素,本節(jié)采用改變彈體長度、保持彈體直徑不變的方式探討彈體長徑比對侵徹深度的影響。數(shù)值模擬得到的侵徹深度如圖11所示,可以看出,在彈徑不變的情況下,隨著彈體長度的增加,侵徹深度也不斷增大。

    圖11 侵徹深度與彈體長徑比的關(guān)系(彈徑不變)Fig. 11 Depth of penetration versus length-to-diameter ratio at a constant projectile diameter

    圖12(a)給出了不同長徑比下應(yīng)力峰值隨深度的變化趨勢,顯然,隨著侵徹深度的增大,地沖擊影響的范圍也會相應(yīng)增大。因此引入相對深度(深度與最終侵徹深度之比)作為評判地沖擊影響范圍的依據(jù),如圖12(b)所示。可以看出,當橫坐標為相對深度時,不同長徑比地沖擊應(yīng)力峰值變化趨勢和數(shù)值基本一致。綜合來看,在彈徑不變的情況下,長徑比對侵徹深度影響較大,而對地沖擊規(guī)律影響不大。

    圖12 彈體長徑比對地沖擊的影響Fig. 12 Effect of length-to-diameter ratio on ground shock wave

    3.2 彈體密度

    在超高速侵徹階段,彈靶密度比也是侵徹深度的重要影響因素之一。考慮到地質(zhì)類材料的密度一般較為接近,因此只改變彈體的材料密度以考慮彈靶密度比對地沖擊的影響。本節(jié)中保持彈體的尺寸不變,彈體材料選擇密度差異較大的鋁、鋼作為對照,其計算條件及結(jié)果如表4所示,表中ρ為彈體密度,ρ為靶體材料密度。

    表4 不同彈體密度情況下的侵徹深度Table 4 Depths of penetration at different projectile densities

    3 種彈體材料下靶體中不同深度的應(yīng)力峰值變化趨勢如圖13 所示。從圖13(a)可以看出,彈體材料密度越高,地沖擊影響的深度也越大,實際中超高速動能武器主要采用鎢合金等高密度金屬材料;而從圖13(b)來看,不同彈體材料下地沖擊影響的相對深度差別不大??傮w來看,不同的彈體密度對侵徹深度有較大影響,但對地沖擊規(guī)律的影響不大。

    圖13 彈體密度對地沖擊的影響Fig. 13 Effects of projectile density on ground shock wave

    3.3 靶體孔隙率

    地質(zhì)材料性質(zhì)差異大,對彈體的侵徹過程及應(yīng)力波的傳播都存在一定程度的影響??紤]到地質(zhì)材料性質(zhì)差異主要由孔隙率引起,本節(jié)中僅改變材料的孔隙率,保持靶體強度等參數(shù)不變,探討孔隙率對地沖擊的影響。

    原石灰?guī)r材料的孔隙率約為15%,在數(shù)值模擬時采用-α 狀態(tài)方程考慮其孔隙壓實過程的性質(zhì)變化?;诂F(xiàn)有的實驗數(shù)據(jù),采用多項式狀態(tài)方程,擬合得到無孔隙對照組巖石材料的狀態(tài)方程曲線,如圖14所示。

    圖14 不同孔隙率巖石材料的狀態(tài)方程Fig. 14 Equations of state for limestones with different porosities

    數(shù)值模擬得到的不同孔隙率靶體內(nèi)不同深度的應(yīng)力峰值變化趨勢如圖15 所示??梢钥闯觯煌紫堵实膽?yīng)力峰值差異較大,且孔隙率越高,其應(yīng)力峰值越低。這主要是由于材料的孔隙率對應(yīng)力波的傳播有較大的影響,孔隙率越高,應(yīng)力波衰減越快,因此孔隙率對地沖擊影響較大,這與Antoun的研究結(jié)論一致。進一步分析可以發(fā)現(xiàn),已有對地沖擊的實驗和數(shù)值模擬研究中采用的巖石不同,得到的地沖擊數(shù)值差異較大,可能是由于孔隙率不同導(dǎo)致的。

    圖15 靶體孔隙率對地沖擊的影響Fig. 15 Effect of target porosity on ground shock wave

    基于上述討論可以看出,彈體的參數(shù)變化(彈體的長徑比、密度等)主要影響其侵徹深度,并會影響地沖擊的破壞深度;但若從相對深度衡量,彈體參數(shù)變化對地沖擊規(guī)律的影響不大。而靶體孔隙率影響靶體中應(yīng)力波的傳播,因此對地沖擊有較大的影響。

    4 結(jié) 論

    為探究超高速動能武器的對地破壞效應(yīng)及其影響因素,采用數(shù)值模擬方法對彈體超高速侵徹的地沖擊進行了研究。利用現(xiàn)有石灰?guī)r靜動態(tài)力學(xué)性能實驗數(shù)據(jù),對Kong-Fang 材料模型進行了參數(shù)標定,并驗證了材料模型及其參數(shù)的準確性。開展了鎢合金彈體超高速侵徹石灰?guī)r靶體的數(shù)值模擬,對不同速度侵徹下的靶體破壞過程和破壞機理進行了分析,得到以下結(jié)論。

    (1)彈體超高速侵徹靶體時,彈靶接觸時會產(chǎn)生瞬時高壓,并以應(yīng)力波的形式在靶體中傳播,造成靶體內(nèi)部不同程度的破壞。

    (2)從靶體內(nèi)部不同測點的地沖擊應(yīng)力時程曲線來看,當侵徹速度較低時,其升壓和降壓時間均較長;而當侵徹速度較高時,升壓時間較短,可近似視為沖擊波。

    (3)從地沖擊沿深度方向的應(yīng)力峰值來看,隨著深度的增大,應(yīng)力峰值不斷降低,地沖擊不斷衰減;隨著侵徹速度的提高,同一深度處的應(yīng)力峰值升高,地沖擊顯著增強;本文中的算例表明,當彈體侵徹速度高于3.0 km/s 時,地沖擊的影響較大。

    (4)從相對深度來看,彈體參數(shù)(彈體長徑比、密度)對地沖擊規(guī)律的影響不大;靶體的特征特別是孔隙率對地沖擊規(guī)律具有較大的影響。

    猜你喜歡
    靶體超高速石灰?guī)r
    靶體結(jié)構(gòu)對前混合水射流噴丸強化應(yīng)力特性的影響
    某石灰?guī)r礦區(qū)巖溶涌水治理處理方法
    Dyeing a “Blue Lagoon”black
    石灰?guī)r
    一種提升石灰?guī)r瀝青混合料路用性能的新型有機添加劑
    石油瀝青(2019年4期)2019-09-02 01:41:52
    中國風(fēng)投行業(yè)迎來超高速發(fā)展
    金橋(2018年1期)2018-09-28 02:24:38
    采用超高速離心鑄造技術(shù)開發(fā)耐磨薄壁鑄管
    四川冶金(2018年1期)2018-09-25 02:39:26
    基于ST-SRCKF的超高速強機動目標跟蹤算法
    彈丸對預(yù)開孔混凝土靶體侵徹的實驗研究
    IR推出堅固可靠的超高速1400V IGBT為感應(yīng)加熱和軟開關(guān)應(yīng)用作出優(yōu)化
    国产精品秋霞免费鲁丝片| 青春草亚洲视频在线观看| 又黄又爽又刺激的免费视频.| 丰满乱子伦码专区| 久久久久国产精品人妻一区二区| 欧美日韩精品成人综合77777| 又黄又爽又刺激的免费视频.| 亚洲欧洲日产国产| 成人影院久久| 亚洲成人一二三区av| kizo精华| 国精品久久久久久国模美| 亚洲三级黄色毛片| 久久影院123| 美女福利国产在线| 青春草视频在线免费观看| 午夜福利网站1000一区二区三区| av网站免费在线观看视频| 久久av网站| 2021少妇久久久久久久久久久| 视频在线观看一区二区三区| 亚洲精品久久久久久婷婷小说| 国产熟女欧美一区二区| 国语对白做爰xxxⅹ性视频网站| 亚洲av欧美aⅴ国产| 亚洲人成网站在线播| 观看美女的网站| 国产极品天堂在线| 亚洲精华国产精华液的使用体验| 男人爽女人下面视频在线观看| 日本猛色少妇xxxxx猛交久久| 22中文网久久字幕| 日韩精品免费视频一区二区三区 | 亚洲国产精品成人久久小说| 国产一区二区在线观看av| 色94色欧美一区二区| 亚洲少妇的诱惑av| videossex国产| 成人国语在线视频| 国产有黄有色有爽视频| 国产色爽女视频免费观看| 中文字幕精品免费在线观看视频 | 久久久久久久久久人人人人人人| 乱人伦中国视频| av女优亚洲男人天堂| 免费久久久久久久精品成人欧美视频 | 乱人伦中国视频| 国精品久久久久久国模美| 日韩制服骚丝袜av| 国产免费一级a男人的天堂| 午夜免费观看性视频| 少妇熟女欧美另类| 91在线精品国自产拍蜜月| 老司机影院毛片| 丝袜喷水一区| 午夜激情av网站| 成人影院久久| 如日韩欧美国产精品一区二区三区 | 久久精品久久久久久久性| 久久久久久久精品精品| 最新中文字幕久久久久| 国产成人精品在线电影| 国产爽快片一区二区三区| 人妻人人澡人人爽人人| 亚洲精品色激情综合| 我的老师免费观看完整版| 黑人巨大精品欧美一区二区蜜桃 | 十分钟在线观看高清视频www| 国产熟女欧美一区二区| 男女无遮挡免费网站观看| 美女cb高潮喷水在线观看| 久久久国产一区二区| av播播在线观看一区| 免费大片黄手机在线观看| 乱码一卡2卡4卡精品| 一本大道久久a久久精品| 国产亚洲午夜精品一区二区久久| 制服诱惑二区| 亚洲成人av在线免费| 80岁老熟妇乱子伦牲交| 晚上一个人看的免费电影| 五月天丁香电影| 亚洲欧美日韩卡通动漫| 最近最新中文字幕免费大全7| 高清午夜精品一区二区三区| 国产不卡av网站在线观看| 日韩一区二区三区影片| 欧美日韩国产mv在线观看视频| 国产精品国产三级国产av玫瑰| 在线看a的网站| 91精品一卡2卡3卡4卡| 国产69精品久久久久777片| 国精品久久久久久国模美| 国产欧美另类精品又又久久亚洲欧美| 久久精品国产亚洲av涩爱| 777米奇影视久久| 久久久精品免费免费高清| 美女视频免费永久观看网站| 大片免费播放器 马上看| 国产亚洲午夜精品一区二区久久| 只有这里有精品99| 夜夜看夜夜爽夜夜摸| 热99久久久久精品小说推荐| 亚洲国产成人一精品久久久| 国产精品久久久久久精品电影小说| 天堂8中文在线网| 日本wwww免费看| 欧美精品人与动牲交sv欧美| 亚洲性久久影院| 极品少妇高潮喷水抽搐| 少妇人妻精品综合一区二区| 久久青草综合色| 亚洲av.av天堂| 久久99精品国语久久久| 一个人看视频在线观看www免费| 在线观看一区二区三区激情| 一级黄片播放器| 中文字幕av电影在线播放| 日韩三级伦理在线观看| 观看av在线不卡| 欧美日韩综合久久久久久| 午夜视频国产福利| 99久久精品国产国产毛片| 97超视频在线观看视频| 国产熟女欧美一区二区| 老女人水多毛片| 99久久精品国产国产毛片| 国产精品人妻久久久影院| 纵有疾风起免费观看全集完整版| 在线观看免费视频网站a站| 建设人人有责人人尽责人人享有的| 黑人欧美特级aaaaaa片| 日韩熟女老妇一区二区性免费视频| 欧美丝袜亚洲另类| 夫妻午夜视频| 黑丝袜美女国产一区| 久久久久久人妻| 视频区图区小说| 午夜免费男女啪啪视频观看| 五月天丁香电影| av天堂久久9| 亚洲综合色惰| 男女边摸边吃奶| 一个人看视频在线观看www免费| 中文字幕亚洲精品专区| 日韩一区二区三区影片| 亚洲国产精品成人久久小说| 成人无遮挡网站| 国产av国产精品国产| 久久ye,这里只有精品| 午夜激情久久久久久久| 久久久久久久久大av| 伊人亚洲综合成人网| 老司机影院毛片| 久久这里有精品视频免费| 97精品久久久久久久久久精品| 欧美性感艳星| 亚洲欧洲精品一区二区精品久久久 | 国产精品国产三级国产av玫瑰| 黄片无遮挡物在线观看| 亚洲美女黄色视频免费看| 国产免费一级a男人的天堂| 日本黄大片高清| 亚洲欧美精品自产自拍| 高清午夜精品一区二区三区| 久久女婷五月综合色啪小说| 免费看av在线观看网站| 十八禁高潮呻吟视频| a 毛片基地| 久久精品久久精品一区二区三区| 女的被弄到高潮叫床怎么办| av在线播放精品| 成人18禁高潮啪啪吃奶动态图 | a级毛片在线看网站| 国产成人freesex在线| videos熟女内射| 狠狠精品人妻久久久久久综合| 国内精品宾馆在线| 国产亚洲午夜精品一区二区久久| 人人妻人人澡人人爽人人夜夜| 精品一区二区三卡| 国产精品久久久久久精品古装| av免费在线看不卡| 国产熟女欧美一区二区| 97在线人人人人妻| 秋霞在线观看毛片| av卡一久久| 欧美国产精品一级二级三级| 国产精品无大码| 日本-黄色视频高清免费观看| 亚洲精品美女久久av网站| 日韩视频在线欧美| av线在线观看网站| 久久精品久久久久久久性| 曰老女人黄片| 99久久精品国产国产毛片| 丰满少妇做爰视频| 国产免费视频播放在线视频| 久久97久久精品| 校园人妻丝袜中文字幕| 国产精品一区二区在线不卡| 五月玫瑰六月丁香| 18禁裸乳无遮挡动漫免费视频| 只有这里有精品99| 2022亚洲国产成人精品| 母亲3免费完整高清在线观看 | 久久精品人人爽人人爽视色| 秋霞伦理黄片| h视频一区二区三区| 大又大粗又爽又黄少妇毛片口| av在线app专区| 久热久热在线精品观看| 欧美亚洲日本最大视频资源| 成人二区视频| 欧美 日韩 精品 国产| 女的被弄到高潮叫床怎么办| 久久精品久久久久久久性| 欧美丝袜亚洲另类| 亚洲成人一二三区av| 在线免费观看不下载黄p国产| 中文字幕精品免费在线观看视频 | 在线观看三级黄色| 亚洲综合色惰| 免费不卡的大黄色大毛片视频在线观看| 亚洲一级一片aⅴ在线观看| 亚洲综合色网址| 男女边吃奶边做爰视频| 中文字幕最新亚洲高清| 中文字幕最新亚洲高清| 啦啦啦中文免费视频观看日本| 丰满迷人的少妇在线观看| 国产亚洲最大av| 久久影院123| 综合色丁香网| 午夜老司机福利剧场| 夜夜骑夜夜射夜夜干| 色婷婷久久久亚洲欧美| 免费观看的影片在线观看| 在线观看国产h片| 草草在线视频免费看| 亚洲av男天堂| 国产精品国产三级专区第一集| 亚洲一区二区三区欧美精品| 91国产中文字幕| 色94色欧美一区二区| 国产亚洲av片在线观看秒播厂| 一区二区三区免费毛片| 少妇的逼好多水| 男女边摸边吃奶| 精品酒店卫生间| 在线观看一区二区三区激情| 亚洲精品第二区| 一区二区三区免费毛片| 免费观看无遮挡的男女| 成人国语在线视频| av一本久久久久| 桃花免费在线播放| 国产男人的电影天堂91| 久久国产亚洲av麻豆专区| 搡女人真爽免费视频火全软件| 免费观看性生交大片5| 亚洲婷婷狠狠爱综合网| 免费观看的影片在线观看| a 毛片基地| 香蕉精品网在线| 又黄又爽又刺激的免费视频.| 少妇的逼水好多| 少妇高潮的动态图| 成人毛片a级毛片在线播放| 日本色播在线视频| 18禁裸乳无遮挡动漫免费视频| 777米奇影视久久| 久久午夜综合久久蜜桃| 欧美日韩国产mv在线观看视频| 下体分泌物呈黄色| 亚洲精品成人av观看孕妇| 纯流量卡能插随身wifi吗| 男人添女人高潮全过程视频| h视频一区二区三区| av专区在线播放| 伦精品一区二区三区| 国产白丝娇喘喷水9色精品| 国产成人精品久久久久久| 久久久午夜欧美精品| 亚洲欧美成人综合另类久久久| 大片电影免费在线观看免费| 日本av手机在线免费观看| 人妻人人澡人人爽人人| 欧美3d第一页| 午夜福利视频在线观看免费| 两个人免费观看高清视频| 国产成人freesex在线| 美女大奶头黄色视频| 人人妻人人澡人人看| 日韩强制内射视频| 国产免费一区二区三区四区乱码| 国产精品99久久99久久久不卡 | 男女高潮啪啪啪动态图| av在线播放精品| 国产成人aa在线观看| 97在线人人人人妻| 亚洲在久久综合| 日韩人妻高清精品专区| 国产国语露脸激情在线看| 国产片内射在线| 午夜福利影视在线免费观看| 亚洲精品第二区| 丝袜喷水一区| 九草在线视频观看| 亚洲av欧美aⅴ国产| 视频区图区小说| 国产高清国产精品国产三级| 国产免费一区二区三区四区乱码| 韩国高清视频一区二区三区| tube8黄色片| 日本黄色日本黄色录像| 精品人妻偷拍中文字幕| www.色视频.com| 国产日韩一区二区三区精品不卡 | 麻豆乱淫一区二区| 日日啪夜夜爽| 伦精品一区二区三区| 高清黄色对白视频在线免费看| av免费观看日本| 有码 亚洲区| 成人无遮挡网站| 香蕉精品网在线| 少妇被粗大的猛进出69影院 | 成年人午夜在线观看视频| 如何舔出高潮| 国产不卡av网站在线观看| 你懂的网址亚洲精品在线观看| 国产黄色免费在线视频| 中国三级夫妇交换| 国产黄片视频在线免费观看| 在线 av 中文字幕| 亚洲av不卡在线观看| 中文字幕久久专区| 免费av不卡在线播放| 欧美日韩在线观看h| 亚洲综合色网址| 性色avwww在线观看| 91精品国产九色| 亚洲人成网站在线观看播放| 制服丝袜香蕉在线| 搡女人真爽免费视频火全软件| 日本猛色少妇xxxxx猛交久久| 能在线免费看毛片的网站| 亚洲av综合色区一区| 两个人免费观看高清视频| 亚洲精品aⅴ在线观看| 精品久久久久久电影网| 一级毛片我不卡| 在线观看www视频免费| 99久久精品一区二区三区| 日本免费在线观看一区| 美女国产视频在线观看| 啦啦啦视频在线资源免费观看| 妹子高潮喷水视频| 亚洲无线观看免费| 91午夜精品亚洲一区二区三区| 国产亚洲欧美精品永久| 99久国产av精品国产电影| 久久精品国产a三级三级三级| 日本黄色日本黄色录像| 欧美日韩精品成人综合77777| av在线播放精品| 婷婷成人精品国产| 久久精品国产亚洲av涩爱| 成人手机av| 精品久久蜜臀av无| videossex国产| 亚洲av在线观看美女高潮| 久久精品国产亚洲av涩爱| 这个男人来自地球电影免费观看 | 久久人妻熟女aⅴ| 桃花免费在线播放| 卡戴珊不雅视频在线播放| 天堂俺去俺来也www色官网| 午夜视频国产福利| 久久热精品热| 有码 亚洲区| 国产成人aa在线观看| 亚洲精品久久成人aⅴ小说 | 亚洲欧美成人综合另类久久久| 日日摸夜夜添夜夜爱| 美女视频免费永久观看网站| 亚洲久久久国产精品| 在线观看免费日韩欧美大片 | 亚洲三级黄色毛片| 国产免费福利视频在线观看| 久久人人爽人人爽人人片va| 男人爽女人下面视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 日本免费在线观看一区| 亚洲综合色网址| 丝袜脚勾引网站| 99久久精品一区二区三区| 日韩一本色道免费dvd| 国产片内射在线| 成人二区视频| 亚洲中文av在线| 亚洲怡红院男人天堂| 老司机影院成人| 久久毛片免费看一区二区三区| 视频在线观看一区二区三区| 99九九线精品视频在线观看视频| 丝袜在线中文字幕| 国产精品一区二区在线观看99| 国产一区二区在线观看av| 一本大道久久a久久精品| 美女中出高潮动态图| 婷婷色麻豆天堂久久| 日本爱情动作片www.在线观看| 自线自在国产av| 夜夜看夜夜爽夜夜摸| 国产亚洲欧美精品永久| 99热全是精品| 国产在线视频一区二区| 九九久久精品国产亚洲av麻豆| av卡一久久| 亚洲精品一区蜜桃| 高清黄色对白视频在线免费看| 各种免费的搞黄视频| 一本大道久久a久久精品| 简卡轻食公司| 高清不卡的av网站| 精品久久久久久电影网| 欧美精品一区二区大全| 91精品一卡2卡3卡4卡| 18禁观看日本| 又黄又爽又刺激的免费视频.| 成人亚洲欧美一区二区av| 色吧在线观看| 曰老女人黄片| 久久精品国产亚洲av涩爱| 亚洲欧美成人综合另类久久久| 精品卡一卡二卡四卡免费| 18禁观看日本| 男女免费视频国产| 精品久久国产蜜桃| 国产亚洲一区二区精品| 伦理电影大哥的女人| 久久久久视频综合| 午夜激情av网站| 国产精品99久久久久久久久| 九九久久精品国产亚洲av麻豆| 大片免费播放器 马上看| 涩涩av久久男人的天堂| 国产女主播在线喷水免费视频网站| 少妇的逼好多水| 97精品久久久久久久久久精品| 久久久久久久久久成人| 国产国语露脸激情在线看| 国产精品久久久久成人av| 日韩精品免费视频一区二区三区 | 永久网站在线| 少妇的逼水好多| 中文欧美无线码| 毛片一级片免费看久久久久| 中文字幕免费在线视频6| 人人妻人人澡人人爽人人夜夜| 国产午夜精品久久久久久一区二区三区| 久久 成人 亚洲| 乱码一卡2卡4卡精品| 日韩av不卡免费在线播放| 午夜免费鲁丝| 久久久久久人妻| 国产免费福利视频在线观看| 狂野欧美激情性xxxx在线观看| 亚洲精品456在线播放app| 99精国产麻豆久久婷婷| 激情五月婷婷亚洲| 飞空精品影院首页| 精品国产露脸久久av麻豆| 免费av中文字幕在线| 狂野欧美激情性xxxx在线观看| 波野结衣二区三区在线| 中国三级夫妇交换| 啦啦啦中文免费视频观看日本| 亚洲精品国产av蜜桃| 一级毛片黄色毛片免费观看视频| 大码成人一级视频| 97精品久久久久久久久久精品| 亚洲国产日韩一区二区| 日韩大片免费观看网站| 少妇 在线观看| 国产精品一区二区在线观看99| 18禁裸乳无遮挡动漫免费视频| 下体分泌物呈黄色| 免费看光身美女| 免费不卡的大黄色大毛片视频在线观看| 国产成人aa在线观看| 精品亚洲成a人片在线观看| 大话2 男鬼变身卡| 色视频在线一区二区三区| av在线app专区| 最新的欧美精品一区二区| 在线观看免费高清a一片| 久久精品熟女亚洲av麻豆精品| 亚洲精品日韩av片在线观看| 我要看黄色一级片免费的| 高清视频免费观看一区二区| 少妇人妻 视频| 国产片特级美女逼逼视频| 国产一区二区三区综合在线观看 | 天天躁夜夜躁狠狠久久av| 国产亚洲午夜精品一区二区久久| 欧美成人精品欧美一级黄| 在线看a的网站| 色视频在线一区二区三区| 一级爰片在线观看| 2022亚洲国产成人精品| 人妻人人澡人人爽人人| 亚洲无线观看免费| 免费不卡的大黄色大毛片视频在线观看| 午夜激情福利司机影院| 美女脱内裤让男人舔精品视频| 亚洲情色 制服丝袜| 国产色爽女视频免费观看| 国产一区二区三区av在线| 国产精品熟女久久久久浪| www.av在线官网国产| 亚洲综合色网址| 18禁观看日本| 最新的欧美精品一区二区| 中文字幕精品免费在线观看视频 | 男女无遮挡免费网站观看| 国产成人91sexporn| 妹子高潮喷水视频| 永久免费av网站大全| 嫩草影院入口| 91精品三级在线观看| 大码成人一级视频| 中文字幕亚洲精品专区| 亚洲精品乱码久久久久久按摩| 毛片一级片免费看久久久久| 国产精品不卡视频一区二区| 国产欧美另类精品又又久久亚洲欧美| 久久精品久久久久久久性| 久久99一区二区三区| 91久久精品国产一区二区成人| 熟女av电影| 黑人高潮一二区| 久久久久久久国产电影| av网站免费在线观看视频| 一区二区三区免费毛片| 国产在线视频一区二区| 狂野欧美激情性xxxx在线观看| tube8黄色片| 最新中文字幕久久久久| 国产视频内射| 人成视频在线观看免费观看| 18禁在线无遮挡免费观看视频| 日韩成人伦理影院| 一边亲一边摸免费视频| 十八禁高潮呻吟视频| 少妇被粗大猛烈的视频| 欧美日韩av久久| 插逼视频在线观看| 成人漫画全彩无遮挡| 日韩视频在线欧美| 七月丁香在线播放| 高清毛片免费看| 男人操女人黄网站| 久久久久精品性色| 精品国产一区二区三区久久久樱花| 精品一区在线观看国产| 国产一区有黄有色的免费视频| 亚洲美女视频黄频| 亚洲国产欧美日韩在线播放| 亚洲激情五月婷婷啪啪| 美女国产高潮福利片在线看| 全区人妻精品视频| 亚洲婷婷狠狠爱综合网| 大香蕉久久网| 午夜免费鲁丝| 久久av网站| 中文字幕最新亚洲高清| 另类亚洲欧美激情| 国产精品三级大全| 久久久a久久爽久久v久久| 亚洲国产av新网站| 18禁观看日本| 久久久久精品久久久久真实原创| 亚洲成人一二三区av| 久久99蜜桃精品久久| 欧美xxxx性猛交bbbb| 大话2 男鬼变身卡| 日韩亚洲欧美综合| 久久人妻熟女aⅴ| 性色avwww在线观看| 久久女婷五月综合色啪小说| 男男h啪啪无遮挡| 黄色一级大片看看| xxxhd国产人妻xxx| 99国产精品免费福利视频| 青青草视频在线视频观看| 尾随美女入室| 国产高清国产精品国产三级| 久久这里有精品视频免费| 日本猛色少妇xxxxx猛交久久| 久久久久久久亚洲中文字幕| 777米奇影视久久| 日韩电影二区| 2018国产大陆天天弄谢| 日韩一区二区视频免费看| 欧美少妇被猛烈插入视频| 亚洲人成网站在线观看播放| 男女高潮啪啪啪动态图| 蜜桃在线观看..| 国产精品一区二区在线不卡| 欧美老熟妇乱子伦牲交| 国产精品女同一区二区软件| 十八禁网站网址无遮挡|