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

    基于廣義波阻抗梯度飛片的準(zhǔn)等熵壓縮技術(shù)*

    2019-06-05 08:05:30陳子博謝普初劉東升王永剛
    爆炸與沖擊 2019年4期

    陳子博,謝普初,劉東升,陳 偉,王永剛

    (寧波大學(xué)沖擊與安全工程教育部重點(diǎn)實(shí)驗(yàn)室,浙江 寧波 315211)

    梯度功能材料是兩種或多種材料復(fù)合且成分和結(jié)構(gòu)呈連續(xù)梯度變化的一種新型復(fù)合材料。波阻抗梯度功能材料(wave impedance functionally graded material,WIFGM)是一類較特殊的梯度功能材料,該材料引起人們的關(guān)注始于材料高壓本構(gòu)和狀態(tài)方程研究[1-2]。在輕氣炮實(shí)驗(yàn)中,波阻抗梯度功能材料制作的飛片在高速碰撞靶材后,可以在靶材中產(chǎn)生具有緩慢上升前沿的壓縮加載波,實(shí)現(xiàn)對(duì)靶材的準(zhǔn)等熵壓縮加載,熵增和溫升遠(yuǎn)小于沖擊絕熱壓縮情況。準(zhǔn)等熵壓縮加載填補(bǔ)了準(zhǔn)靜態(tài)等熵加載和沖擊絕熱加載之間空白,提供了一種考察物質(zhì)壓縮形態(tài)(尤其是偏離Hugoniot狀態(tài))的新實(shí)驗(yàn)途徑,可以得到范圍更廣、壓縮度更高、精度更高的材料物態(tài)方程參數(shù)。另外,把波阻抗梯度飛片的準(zhǔn)等熵壓縮性和二級(jí)輕氣炮相結(jié)合可以發(fā)展三級(jí)炮超高速發(fā)射技術(shù),彈丸速度可以達(dá)到10 km/s以上[3],這在航天器空間碎片防護(hù)領(lǐng)域具有重要的應(yīng)用背景。

    20世紀(jì)80年代,美國圣地亞核武器國家實(shí)驗(yàn)室首先開展了波阻抗梯度材料的制備技術(shù)研究。Barker等[1]采用顆粒共沉降成型結(jié)合高壓粘接方法制備了Ta/Cu/Al/TPX金屬-聚合物體系的波阻抗梯度飛片,可以在靶材中產(chǎn)生準(zhǔn)等熵壓縮波,但該波阻抗梯度飛片的力學(xué)強(qiáng)度較差,在發(fā)射過程中飛片易發(fā)生破壞。Chhabildas等[2]采用焊接方法制備了Ta/Cu/Ti/Al/Mg/TPX體系疊層飛片,但該飛片的波阻抗沿厚度方向是以階躍方式變化的,導(dǎo)致其產(chǎn)生的壓力剖面平順性較差。后來,美國勞倫斯利弗莫爾國家實(shí)驗(yàn)室Martin等[4]采用流延成型法制備了7種組分和19種組分的Mg/Cu體系波阻抗梯度材料,獲得了平滑的壓力剖面。近來,Yep等[5]進(jìn)一步發(fā)展了流延成型技術(shù),制備了近100種組分的Mg/Cu/W體系梯度飛片,其具有更寬廣波阻抗變化范圍,可以實(shí)現(xiàn)500 GPa壓力值的準(zhǔn)等熵壓縮。從20世紀(jì)90年代開始,武漢理工大學(xué)和中國工程物理研究院流體物理研究所合作開展波阻抗梯度材料制備技術(shù)研究[6-8],發(fā)展了放電等離子快速粉末燒結(jié)法,研制出的梯度飛片已成功應(yīng)用于準(zhǔn)等熵壓縮加載、超高速驅(qū)動(dòng)、太帕量級(jí)壓力的高精度平面加載以及復(fù)雜加、卸載等動(dòng)高壓實(shí)驗(yàn)技術(shù)領(lǐng)域。傳統(tǒng)的波阻抗梯度功能材料制備多是選擇多種不同密度不同波速的金屬材料進(jìn)行復(fù)合。為了得到材料波阻抗準(zhǔn)連續(xù)變化,需要進(jìn)行非常多組分配比,鋪層和燒結(jié)制作工序異常繁雜,質(zhì)量也難以精確控制,成品率低,價(jià)格昂貴,這大大制約了波阻抗梯度功能材料實(shí)際應(yīng)用,迫切需要發(fā)展新的制備工藝。最近,Taylar等[9]采用激光立體光刻技術(shù)制備了波阻抗梯度陶瓷飛片,實(shí)現(xiàn)對(duì)平板撞擊實(shí)驗(yàn)中沖擊加載波形的有效調(diào)控。

    增材制造(additive manufacturing, AM)技術(shù)是基于計(jì)算機(jī)輔助設(shè)計(jì)數(shù)據(jù),采用材料逐層累積的方法來制造零件[10-11]。材料結(jié)構(gòu)的計(jì)算機(jī)三維自由設(shè)計(jì)和增材制造技術(shù)的完美結(jié)合為新型功能材料或結(jié)構(gòu)材料的研發(fā)與制備提供新的途徑。本文中擬采用激光選區(qū)熔化(selective laser melting,SLM)金屬增材制造技術(shù),通過幾何結(jié)構(gòu)設(shè)計(jì)來研發(fā)新型的廣義波阻抗梯度飛片材料。首先,基于應(yīng)力波傳播理論來探討波阻抗梯度飛片的幾何構(gòu)型設(shè)計(jì);隨后對(duì)“結(jié)構(gòu)設(shè)計(jì)型”波阻抗梯度飛片的準(zhǔn)等熵壓縮加載特性進(jìn)行數(shù)值模擬,重點(diǎn)討論結(jié)構(gòu)設(shè)計(jì)的關(guān)鍵參數(shù)和碰撞速度對(duì)加載波形的影響規(guī)律;最后,采用增材制造技術(shù)制備波阻抗梯度飛片,開展氣炮實(shí)驗(yàn),對(duì)比分析了實(shí)驗(yàn)與數(shù)值計(jì)算結(jié)果,驗(yàn)證新型波阻抗梯度飛片的設(shè)計(jì)可行性和適用性。

    1 波阻抗梯度飛片的結(jié)構(gòu)設(shè)計(jì)

    材料密度ρ0與波速c的乘積定義為材料波阻抗。在研究變截面桿中波傳播問題時(shí)人們引入廣義波阻抗的概念,定義為材料密度ρ0、波速c以及截面面積A的乘積。由此來看,除了改變材料密度和波速之外,還可以通過改變截面面積來實(shí)現(xiàn)廣義波阻抗的梯度變化。王禮立等[12]通過引入廣義波阻抗把錐形桿中的縱波傳播問題近似化為一維應(yīng)力平面縱波問題來處理,通過特征線方法探討了應(yīng)力波在錐形桿中的傳播規(guī)律,結(jié)果表明,錐形桿具有調(diào)節(jié)應(yīng)力波波形和強(qiáng)弱的特性。羅鑫等[13]對(duì)錐形桿共軸撞擊進(jìn)行數(shù)值分析,結(jié)果表明,撞擊所產(chǎn)生的加載波陣面具有較長的前沿升時(shí),同時(shí)還消除了波形振蕩。陶俊林等[14]在SHPB實(shí)驗(yàn)中將圓柱形子彈改變?yōu)殄F形子彈,實(shí)現(xiàn)對(duì)一維應(yīng)力條件下加載波形的調(diào)整,上升沿時(shí)間顯著增大。受此啟發(fā),在輕氣炮實(shí)驗(yàn)中也采用類似的設(shè)計(jì)思想。輕氣炮實(shí)驗(yàn)系統(tǒng)中采用薄板飛片作為加載手段,以便在樣品中產(chǎn)生一維沖擊波,因此不可以直接把薄板飛片改變?yōu)殄F形飛片。這里,我們把飛片分為兩部分,一部分作為飛片基座,然后在基座上密排疊加許多小正四棱錐臺(tái),猶如基座上長出一根根變截面針刺,簡稱為“針床型”波阻抗梯度飛片,其幾何結(jié)構(gòu)示意圖如圖1所示,圖中e為基座高,h為四棱錐臺(tái)高,m為正四棱錐臺(tái)下底邊寬,n為正四棱錐臺(tái)上頂邊寬。當(dāng)“針床型”波阻抗梯度飛片高速撞擊樣品時(shí),每一個(gè)小正四棱錐臺(tái)都可以看作“點(diǎn)”加載脈沖源,會(huì)產(chǎn)生一系列具有緩慢上升前沿的類球面波,隨著傳播距離增大球面波將相互疊加,最終可以得到近似平面的加載波形,從而滿足一維應(yīng)變加載條件。

    圖1 “針床形”廣義波阻抗梯度飛片幾何結(jié)構(gòu)示意圖Fig.1 Schematic diagram of the “needle-bed” generalized wave impedance gradient flyer

    2 計(jì)算模型

    為了驗(yàn)證上述設(shè)計(jì)思想的可行性,首先開展“針床型”波阻抗梯度飛片高速撞擊靶板的數(shù)值模擬。在數(shù)值計(jì)算之前,還需要選擇合適的算法。如果采用拉格朗日算法,在“針床型”波阻抗梯度飛片高速碰撞靶板時(shí),飛片和靶板都會(huì)產(chǎn)生由于局部接觸引起局部大變形,從而導(dǎo)致局部網(wǎng)格單元發(fā)生畸變問題,使得計(jì)算提前中止。通常的解決方案是采用侵蝕算法,把嚴(yán)重畸變的網(wǎng)格單元從模型中刪除,但網(wǎng)格單元?jiǎng)h除降低了計(jì)算精度。為了解決大變形網(wǎng)格畸變問題,人們常采用歐拉算法或者光滑粒子流體動(dòng)力學(xué)(SPH)算法。相比拉格朗日算法,雖然歐拉算法解決了網(wǎng)格畸變問題,但更加耗費(fèi)機(jī)時(shí)。本文中擬采用SPH算法來計(jì)算“針床型”波阻抗梯度飛片高速碰撞靶板問題。SPH算法是將連續(xù)的固體離散為一系列具有質(zhì)量、速度和能量的粒子,通過求解質(zhì)點(diǎn)組的動(dòng)力學(xué)方程和跟蹤每個(gè)質(zhì)點(diǎn)的運(yùn)動(dòng)軌道求得整個(gè)系統(tǒng)的力學(xué)行為。該方法是一種無網(wǎng)格拉格朗日數(shù)值計(jì)算方法,廣泛應(yīng)用于處理大變形問題。

    “針床型”波阻抗梯度飛片上有非常多的小正四棱錐臺(tái),周期性排列。小正四棱錐臺(tái)頂面尺寸非常小,為了保證足夠的計(jì)算精度,建模時(shí)需要?jiǎng)澐衷S多尺寸非常小粒子,整個(gè)模型的粒子數(shù)龐大,計(jì)算耗費(fèi)時(shí)間非常長。為此,這里采用了簡化模型,即選擇一個(gè)典型單胞來計(jì)算,單胞SPH計(jì)算模型如圖2所示,通過在x軸和y軸方向設(shè)置固定邊界條件來保證僅僅z軸方向發(fā)生應(yīng)變,以滿足一維應(yīng)變加載條件。典型計(jì)算單胞的邊長為2 mm,基座高2 mm,正四棱錐臺(tái)的底邊寬1 mm、頂邊寬0.2 mm、高6 mm,靶板邊長2 mm、高7 mm。模型中SPH粒子的直徑為100 μm。由于高速碰撞的物理過程,目前計(jì)算中常采用Johnson-Cook本構(gòu)模型,具體包含線性熱彈性、Mises屈服準(zhǔn)則、各向同性應(yīng)變硬化、應(yīng)變率硬化、各向同性損傷引起的軟化,其等效Mises屈服應(yīng)力表示如下:

    式中:A、B、n、C和m是模型常數(shù); εeq是等效塑性應(yīng)變是無量綱化等效塑性應(yīng)變率,為參考應(yīng)變率;T?=(T-T0)/(Tm-T0)為無量綱溫度,其中T0是室溫,Tm是材料的熔點(diǎn)。計(jì)算時(shí)波阻抗梯度飛片和靶板材料均為鋼,表1給出了數(shù)值模擬中采用的材料本構(gòu)模型參數(shù)。

    圖2 波阻抗梯度飛片氣炮實(shí)驗(yàn)的SPH單胞計(jì)算模型Fig.2 SPH cell simulation model of impedance gradient flyer in gas gun experiments

    表1 飛片和靶板材料本構(gòu)和狀態(tài)方程參數(shù)Table 1 Flyer and target material parameters

    3 計(jì)算結(jié)果討論與分析

    3.1 平面波產(chǎn)生和傳播特性分析

    圖3 靶板內(nèi)不同時(shí)刻的應(yīng)力波傳播云圖Fig.3 Stress wave propagation contours in the specimen at different times

    首先通過一個(gè)算例來討論一下“針床型”梯度飛片高速撞擊靶板后靶板中應(yīng)力波傳播特性。圖3給出了撞擊速度350 m/s下靶板中沿z軸方向傳播的不同時(shí)刻的壓力云圖。初始階段可以清楚地看到每一個(gè)小四棱錐臺(tái)接觸處產(chǎn)生了球面加載波,彼此獨(dú)立往前傳播;隨著傳播距離增大,球面加載波開始相互作用;隨著傳播距離進(jìn)一步增大,最后可以看到開始的球面波逐漸轉(zhuǎn)變?yōu)槠矫娌?,從而?shí)現(xiàn)對(duì)靶板平面加載。在實(shí)際應(yīng)用中,為了保證靶板承受的是平面加載,需要在靶板前增加背板。

    3.2 撞擊速度和幾何結(jié)構(gòu)參數(shù)的影響

    3.2.1 撞擊速度

    首先討論飛片撞擊速度對(duì)準(zhǔn)等熵壓縮加載特性的影響。計(jì)算中,保持典型單胞幾何結(jié)構(gòu)參數(shù)不變,僅僅改變飛片撞擊速度。圖4給出了350 m/s速度下不同時(shí)刻飛片和靶板的變形過程和塑性應(yīng)變?cè)茍D,從圖中可以看到:(1) 飛片上的小四棱錐臺(tái)與靶板的接觸處首先出現(xiàn)了比較嚴(yán)重的塑性變形,塑性變形區(qū)隨著時(shí)間逐漸增大;(2) 小四棱錐臺(tái)之間間隙不斷被由于嚴(yán)重塑性變形而破損的物質(zhì)填充,小四棱錐臺(tái)的高度不斷減小,其整體壓縮量超過50%。圖5給出了200、350、500、800 m/s撞擊條件下靶板自由面速度時(shí)程曲線,從圖中看到:(1) 靶板中產(chǎn)生了緩慢上升的加載波形,在初始階段自由面速度波形出現(xiàn)了一些振蕩,但隨著速度增大,振蕩現(xiàn)象逐漸轉(zhuǎn)弱,速度時(shí)程曲線趨于光滑,同時(shí)還觀察到初始的振蕩現(xiàn)象還依賴于飛片撞擊速度,高速碰撞時(shí)振蕩現(xiàn)象比較微弱,且持續(xù)時(shí)間比較短;(2) 加載波形的上升沿時(shí)間隨著飛片撞擊速度的增大而不斷減小。自由面速度初始振蕩現(xiàn)象和上升沿時(shí)間對(duì)撞擊速度的依賴性主要?dú)w因于飛片的塑性變形過程及塑性波傳播。圖6給出了不同撞擊速度下飛片最終變形圖,隨著撞擊速度增大,小四棱錐臺(tái)壓縮量不斷增加。

    圖7給出了自由面速度波形的上升沿時(shí)間隨著撞擊速度的變化曲線,結(jié)果顯示兩者之間呈指數(shù)衰減關(guān)系,擬合方程為:

    圖4 不同時(shí)刻飛片和靶板的塑性應(yīng)變?cè)茍DFig.4 Plastic strain contours in the flyer and specimen at different times

    式中:t為自由面速度波形的上升沿時(shí)間,v0為飛片的撞擊速度。另外,靶板自由面速度峰值也依賴于飛片撞擊速度。圖8給出了自由面速度的峰值與撞擊速度的比值隨撞擊速度的變化曲線,隨著撞擊速度增大,自由面速度的峰值更接近于撞擊速度,比值趨近于1,兩者之間呈指數(shù)增長關(guān)系,擬合方程為:

    式中:vmax/v0為靶板自由面速度峰值與飛片撞擊速度比,v0為飛片的撞擊速度。

    圖5 不同撞擊速度下靶板自由面速度時(shí)程曲線Fig.5 Free surface velocity profiles of specimen at different impact velocities

    圖6 不同撞擊速度下飛片變形的最終形態(tài)Fig.6 Final deformation of the flyer at different impact velocities

    圖7 壓縮波上升沿時(shí)間隨著撞擊速度的變化曲線Fig.7 Compression wave front time vs.impact velocities

    圖8 自由面速度峰值與撞擊速度比值隨著撞擊速度的變化曲線Fig.8 Ratios between peak velocity and impact velocity vs.impact velocity

    3.2.2 幾何結(jié)構(gòu)參數(shù)

    除了上述的飛片撞擊速度,“針床型”飛片的幾何結(jié)構(gòu)參數(shù)也對(duì)準(zhǔn)等熵壓縮加載特性有顯著的影響。首先討論小四棱錐臺(tái)高度的影響。其他幾何參數(shù)不變條件下,四棱錐臺(tái)高(h)分別選擇6、8、10 mm,以相同撞擊速度(350 m/s)撞擊靶板。圖9給出了靶板自由表面速度時(shí)程曲線,從中看到加載波形的上升沿時(shí)間隨著錐臺(tái)高度增加而增大,但速度峰值隨著錐臺(tái)高度增加而減小。接著,討論小四棱錐臺(tái)錐角的影響,這里定義正四棱錐傾斜面與底面夾角稱為錐角,如圖10所示,保持小四棱錐臺(tái)上底邊寬度不變,通過改變下底邊邊長來實(shí)現(xiàn)。下底邊邊長分別取 1、1.5、2和 3 mm,對(duì)應(yīng)錐角分別為 86.2°、83.8°、81.4°、76.9°。圖11給出了相同撞擊速度(350 m/s)不同錐角條件下靶板自由面速度時(shí)程曲線的對(duì)比,圖中顯示:上升沿時(shí)間隨著錐角增大而減小,但峰值速度隨著錐角增大而增大。

    圖9 不同四棱錐臺(tái)高條件下靶板自由面速度時(shí)程曲線對(duì)比Fig.9 Comparison of the free surface velocity profiles for the specimen with different heights of the pyramid

    圖12(a)~(b)分別給出了上升沿時(shí)間和峰值速度隨錐角的變化關(guān)系曲線,上升沿時(shí)間隨錐角增大呈線性減小關(guān)系,擬合關(guān)系式為:

    圖10 四棱錐臺(tái)的錐角定義示意圖Fig.10 Schematic diagram of the cone angle of the pyramid

    圖11 不同錐角條件下靶板自由面速度時(shí)程曲線對(duì)比Fig.11 Comparison of free surface velocity profiles for the pyramid specimen with different cone angles

    圖12 壓縮波上升沿時(shí)間和峰值速度隨著正四棱錐臺(tái)錐角的變化曲線Fig.12 Compression wave front time and peak velocity vs.the cone angle of the pyramid

    而峰值速度與錐角之間呈線性增長關(guān)系,擬合關(guān)系式為:

    式中:t表示波陣面前沿上升沿時(shí)間,μs;β表示梯度飛片正四棱錐傾斜面與底面夾角形成錐角;vmax為靶板自由面速度峰值,m/s。

    最后討論正四棱錐臺(tái)的上底邊寬的影響。保持下底邊和錐臺(tái)高度不變,上底邊寬分別為0.2、0.4和0.6 mm,以相同初始速度(350 m/s)沖擊靶板。圖13給出了不同上底邊寬帶條件下的靶板自由面速度的時(shí)程曲線,從圖中可知:自由面速度曲線初始階段出現(xiàn)了明顯的臺(tái)階,隨著上底邊寬度增大,臺(tái)階高度也增大,為了獲得比較平順的準(zhǔn)等熵壓縮加載條件,四棱錐臺(tái)上底邊寬度不易過大;同時(shí)上底邊寬度也影響著上升沿時(shí)間,寬度增大,上升沿時(shí)間變小。

    圖13 不同上底邊寬度條件下靶板自由面速度時(shí)程曲線對(duì)比Fig.13 Comparison of free surface velocity profiles for the specimen with different widths of the upper edge

    4 實(shí) 驗(yàn)

    4.1 廣義波阻抗梯度飛片制備

    利用Pro-E軟件建立了廣義波阻抗梯度飛片三維模型,如圖14(a)所示,飛片直徑53 mm,厚8 mm,基座高2 mm,四棱錐臺(tái)高6 mm,正四棱錐臺(tái)下底邊寬2 mm。采用德國EOS—M280激光選區(qū)熔化金屬3D打印機(jī)進(jìn)行制備,粉末材料采用GP1不銹鋼粉末,制備完成的飛片如圖14(b)所示。

    圖14 廣義波阻抗梯度飛片的3D幾何模型和樣品Fig.14 3D model and product of generalized wave impedance gradient flyer

    4.2 氣炮實(shí)驗(yàn)

    實(shí)驗(yàn)在一級(jí)輕氣炮上進(jìn)行,實(shí)驗(yàn)裝置示意圖如圖15所示,為了保證擊靶平面性,飛片安裝于鋁合金彈托上,通過高壓氣體驅(qū)動(dòng)彈托和飛片高速撞擊靶板,飛片擊靶速度由測(cè)速探針來測(cè)量,采用全光纖激光位移干涉系統(tǒng)(DISAR)來監(jiān)測(cè)靶板自由面速度時(shí)程曲線。實(shí)驗(yàn)中測(cè)得飛片撞擊速度為348 m/s。圖16給出了實(shí)測(cè)的靶板自由面速度時(shí)程曲線。按照實(shí)驗(yàn)條件開展了相關(guān)的數(shù)值模擬工作,計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比也如圖16所示,兩者基本一致,從而驗(yàn)證了設(shè)計(jì)思想和數(shù)值計(jì)算的有效性。

    圖15 輕氣炮實(shí)驗(yàn)裝置示意圖Fig.15 Schematic of the gas gun device

    圖16 實(shí)驗(yàn)與數(shù)值計(jì)算得到靶板自由面速度時(shí)程曲線對(duì)比Fig.16 Experimental and numerical free surface velocities of the target

    5 結(jié) 論

    (1) “針床型”小四棱錐臺(tái)高速碰撞靶板后,在靶板中產(chǎn)生密集分布的球面加載波,隨著傳播距離的增大,球面加載波相互疊加形成緩慢上升前沿的平面加載波,從而實(shí)現(xiàn)對(duì)靶板的準(zhǔn)等熵壓縮加載;

    (2) 準(zhǔn)等熵壓縮加載的特性強(qiáng)烈依賴于飛片擊靶的速度,隨著速度的增加,加載波上升沿時(shí)間減小,而峰值速度增大;

    (3) 廣義波阻抗飛片的關(guān)鍵幾何參數(shù)也對(duì)準(zhǔn)等熵壓縮加載特性影響顯著,隨著四棱錐臺(tái)高增大,加載波上升沿時(shí)間增長,而峰值速度減小;相反,保持四棱錐臺(tái)高不變,只改變錐角的大小,隨著錐角的增大加載波上升沿時(shí)間縮短,而峰值速度增大;小四棱錐臺(tái)變鈍(即上底邊寬度增大)會(huì)導(dǎo)致加載波初始階段出現(xiàn)明顯的臺(tái)階,為了獲得比較平順的準(zhǔn)等熵壓縮加載條件,四棱錐臺(tái)需要保持尖銳度。

    最新的欧美精品一区二区| 欧美日韩精品网址| 国产精品香港三级国产av潘金莲 | 黄色一级大片看看| 超色免费av| www.精华液| 我要看黄色一级片免费的| 久久久久久人人人人人| 国产在线视频一区二区| 日本午夜av视频| 亚洲av电影在线观看一区二区三区| av女优亚洲男人天堂| 极品少妇高潮喷水抽搐| 欧美最新免费一区二区三区| 亚洲一级一片aⅴ在线观看| 看免费av毛片| 熟妇人妻不卡中文字幕| 女性生殖器流出的白浆| 日韩精品有码人妻一区| 午夜91福利影院| 天天躁夜夜躁狠狠躁躁| 黄色一级大片看看| 日本vs欧美在线观看视频| 欧美日韩视频高清一区二区三区二| 狂野欧美激情性bbbbbb| 亚洲,欧美精品.| 视频在线观看一区二区三区| 欧美久久黑人一区二区| 黄色 视频免费看| 另类亚洲欧美激情| 亚洲精品,欧美精品| av一本久久久久| 9热在线视频观看99| 美女午夜性视频免费| 免费女性裸体啪啪无遮挡网站| 性少妇av在线| 欧美乱码精品一区二区三区| 美女中出高潮动态图| 国产有黄有色有爽视频| 中文精品一卡2卡3卡4更新| 丝袜在线中文字幕| 伦理电影大哥的女人| 欧美精品一区二区免费开放| av在线播放精品| 亚洲精品国产区一区二| 国产成人精品在线电影| 亚洲成国产人片在线观看| 人妻人人澡人人爽人人| 久久婷婷青草| 中文欧美无线码| 亚洲av电影在线观看一区二区三区| 麻豆乱淫一区二区| 国产在线视频一区二区| 一区二区av电影网| av卡一久久| 一级a爱视频在线免费观看| 少妇 在线观看| 中文字幕亚洲精品专区| 国产色婷婷99| 香蕉丝袜av| 操美女的视频在线观看| 男男h啪啪无遮挡| 18禁观看日本| 欧美日韩一区二区视频在线观看视频在线| 国产成人免费观看mmmm| 男人爽女人下面视频在线观看| 国产日韩欧美在线精品| 国产极品天堂在线| 久久综合国产亚洲精品| 老汉色∧v一级毛片| 91精品三级在线观看| 国产一区二区三区av在线| 色视频在线一区二区三区| 热re99久久精品国产66热6| 人人妻人人爽人人添夜夜欢视频| 男人操女人黄网站| 久久精品国产综合久久久| 亚洲国产精品国产精品| 国产精品久久久久久人妻精品电影 | 色婷婷久久久亚洲欧美| 久久久欧美国产精品| 亚洲精品自拍成人| 黄片播放在线免费| 国产欧美日韩一区二区三区在线| 欧美最新免费一区二区三区| 免费久久久久久久精品成人欧美视频| 水蜜桃什么品种好| 亚洲久久久国产精品| 69精品国产乱码久久久| 一级毛片我不卡| 亚洲国产毛片av蜜桃av| 男的添女的下面高潮视频| svipshipincom国产片| 精品少妇一区二区三区视频日本电影 | 日韩人妻精品一区2区三区| 免费av中文字幕在线| 亚洲欧美日韩另类电影网站| 欧美激情极品国产一区二区三区| 最新在线观看一区二区三区 | 国产亚洲一区二区精品| videosex国产| 亚洲精品国产一区二区精华液| 国产 精品1| 亚洲精华国产精华液的使用体验| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产av成人精品| 欧美精品高潮呻吟av久久| 高清黄色对白视频在线免费看| 国产成人啪精品午夜网站| 国产黄色视频一区二区在线观看| 看免费成人av毛片| 爱豆传媒免费全集在线观看| 欧美日韩av久久| 久久天堂一区二区三区四区| 欧美日韩福利视频一区二区| 热99国产精品久久久久久7| 秋霞伦理黄片| 99精品久久久久人妻精品| 精品少妇内射三级| 一区福利在线观看| 街头女战士在线观看网站| 久久99热这里只频精品6学生| 国产爽快片一区二区三区| 亚洲成色77777| 成人国产av品久久久| 你懂的网址亚洲精品在线观看| 成年女人毛片免费观看观看9 | 成年av动漫网址| 欧美最新免费一区二区三区| 男女国产视频网站| 老司机靠b影院| 欧美最新免费一区二区三区| 国产片内射在线| 久久久久久免费高清国产稀缺| av有码第一页| 午夜久久久在线观看| 97精品久久久久久久久久精品| 欧美黄色片欧美黄色片| 国产成人欧美在线观看 | 久久亚洲国产成人精品v| 一区在线观看完整版| 91精品伊人久久大香线蕉| 在线观看人妻少妇| 国产精品成人在线| 日韩免费高清中文字幕av| 午夜免费观看性视频| 中文欧美无线码| 亚洲成人手机| 一区福利在线观看| 亚洲欧美成人综合另类久久久| 我要看黄色一级片免费的| 亚洲人成网站在线观看播放| 亚洲av成人精品一二三区| 色婷婷久久久亚洲欧美| 国产精品.久久久| av又黄又爽大尺度在线免费看| 亚洲五月色婷婷综合| 捣出白浆h1v1| 欧美中文综合在线视频| 亚洲专区中文字幕在线 | 麻豆乱淫一区二区| 一级毛片黄色毛片免费观看视频| 中国国产av一级| 亚洲精品国产av蜜桃| 亚洲第一区二区三区不卡| 欧美精品av麻豆av| 久久久国产欧美日韩av| 制服诱惑二区| 欧美日韩av久久| 久久精品国产亚洲av涩爱| 国产成人免费观看mmmm| 蜜桃国产av成人99| 18在线观看网站| 一二三四中文在线观看免费高清| 国产av一区二区精品久久| 看非洲黑人一级黄片| 夜夜骑夜夜射夜夜干| 国产精品久久久人人做人人爽| 最黄视频免费看| 亚洲四区av| 亚洲国产精品999| 日韩制服骚丝袜av| 韩国高清视频一区二区三区| 你懂的网址亚洲精品在线观看| 亚洲欧洲精品一区二区精品久久久 | 国产一区二区三区综合在线观看| 午夜老司机福利片| 亚洲精品久久午夜乱码| 高清不卡的av网站| 久久97久久精品| 免费在线观看视频国产中文字幕亚洲 | 十八禁网站网址无遮挡| 在线看a的网站| 色吧在线观看| 卡戴珊不雅视频在线播放| 国产无遮挡羞羞视频在线观看| 日本欧美视频一区| 国产av码专区亚洲av| 亚洲国产av新网站| 卡戴珊不雅视频在线播放| 亚洲精品国产区一区二| 宅男免费午夜| 欧美在线一区亚洲| 精品久久蜜臀av无| 欧美 亚洲 国产 日韩一| 日本91视频免费播放| 欧美日韩综合久久久久久| 99热网站在线观看| 99久久综合免费| 亚洲精品一区蜜桃| videosex国产| 性少妇av在线| 一本大道久久a久久精品| 波多野结衣av一区二区av| 国产有黄有色有爽视频| 亚洲av成人精品一二三区| 国产免费视频播放在线视频| 人人澡人人妻人| 色播在线永久视频| 国产日韩欧美亚洲二区| 波野结衣二区三区在线| 久久久精品国产亚洲av高清涩受| 日日摸夜夜添夜夜爱| 亚洲国产成人一精品久久久| 免费看不卡的av| 最新的欧美精品一区二区| 青春草亚洲视频在线观看| 天堂俺去俺来也www色官网| 免费看不卡的av| 亚洲美女黄色视频免费看| 综合色丁香网| 婷婷成人精品国产| 精品人妻熟女毛片av久久网站| 超色免费av| 久久影院123| 黄网站色视频无遮挡免费观看| 天天影视国产精品| 成年美女黄网站色视频大全免费| 成人国产av品久久久| 久久久久久久精品精品| 美女视频免费永久观看网站| 亚洲男人天堂网一区| 午夜日本视频在线| 中文字幕最新亚洲高清| 丰满饥渴人妻一区二区三| 91精品三级在线观看| 日韩中文字幕视频在线看片| 高清在线视频一区二区三区| 国产精品无大码| 中文乱码字字幕精品一区二区三区| 欧美人与善性xxx| 大片电影免费在线观看免费| 久久久国产欧美日韩av| 少妇人妻 视频| 国产日韩欧美视频二区| 纵有疾风起免费观看全集完整版| 午夜福利视频精品| 国产一区二区三区综合在线观看| 亚洲情色 制服丝袜| 青春草视频在线免费观看| 亚洲欧美成人综合另类久久久| 王馨瑶露胸无遮挡在线观看| 老熟女久久久| 人人妻人人澡人人爽人人夜夜| 亚洲熟女精品中文字幕| 国产成人a∨麻豆精品| 十八禁高潮呻吟视频| 免费女性裸体啪啪无遮挡网站| 又粗又硬又长又爽又黄的视频| 水蜜桃什么品种好| 亚洲国产精品999| 国产国语露脸激情在线看| 一区二区三区乱码不卡18| 亚洲成人免费av在线播放| 伊人亚洲综合成人网| 中文字幕最新亚洲高清| 在线亚洲精品国产二区图片欧美| 亚洲成人av在线免费| 啦啦啦在线观看免费高清www| 少妇 在线观看| 亚洲精品美女久久av网站| 国产成人一区二区在线| 国产成人系列免费观看| 国产av精品麻豆| 日韩大片免费观看网站| 一边摸一边抽搐一进一出视频| 欧美在线黄色| 亚洲人成77777在线视频| 午夜久久久在线观看| 欧美日韩亚洲高清精品| 美女脱内裤让男人舔精品视频| 又粗又硬又长又爽又黄的视频| 少妇 在线观看| 尾随美女入室| 成人国产麻豆网| 在线 av 中文字幕| 黄片小视频在线播放| 久久久久精品久久久久真实原创| 人妻人人澡人人爽人人| 成人国语在线视频| 黄色视频在线播放观看不卡| 欧美老熟妇乱子伦牲交| 亚洲免费av在线视频| 国产 精品1| 国产精品二区激情视频| 亚洲av男天堂| 久热爱精品视频在线9| 精品人妻熟女毛片av久久网站| 精品一区二区免费观看| 曰老女人黄片| 国产毛片在线视频| 久久久精品国产亚洲av高清涩受| av女优亚洲男人天堂| 久久av网站| 性少妇av在线| 亚洲熟女毛片儿| 欧美激情极品国产一区二区三区| 高清av免费在线| 女的被弄到高潮叫床怎么办| 肉色欧美久久久久久久蜜桃| 精品一区在线观看国产| 成年av动漫网址| 电影成人av| 午夜日韩欧美国产| 在线观看免费视频网站a站| 69精品国产乱码久久久| av不卡在线播放| 1024香蕉在线观看| 菩萨蛮人人尽说江南好唐韦庄| 九色亚洲精品在线播放| 看十八女毛片水多多多| 色综合欧美亚洲国产小说| 婷婷色av中文字幕| 男女国产视频网站| 久久久久久久国产电影| 亚洲欧美一区二区三区黑人| av一本久久久久| 狠狠精品人妻久久久久久综合| 精品福利永久在线观看| 9191精品国产免费久久| 人人妻人人爽人人添夜夜欢视频| 国产激情久久老熟女| 亚洲少妇的诱惑av| 观看av在线不卡| av不卡在线播放| 亚洲av电影在线进入| 91老司机精品| 人妻 亚洲 视频| 国产精品亚洲av一区麻豆 | 国产成人精品在线电影| 丰满乱子伦码专区| 人人妻,人人澡人人爽秒播 | 视频在线观看一区二区三区| 久久精品亚洲熟妇少妇任你| 午夜福利免费观看在线| 亚洲av电影在线观看一区二区三区| 青春草视频在线免费观看| 最新在线观看一区二区三区 | 婷婷色av中文字幕| 精品福利永久在线观看| 丰满少妇做爰视频| 妹子高潮喷水视频| 在线观看人妻少妇| 亚洲一区二区三区欧美精品| 国产人伦9x9x在线观看| 飞空精品影院首页| 大话2 男鬼变身卡| 青青草视频在线视频观看| 亚洲精品一区蜜桃| 男女边摸边吃奶| av电影中文网址| avwww免费| 成人手机av| 成年人午夜在线观看视频| 精品人妻熟女毛片av久久网站| 国产男人的电影天堂91| 亚洲国产欧美日韩在线播放| 大话2 男鬼变身卡| 啦啦啦在线免费观看视频4| 精品国产一区二区三区四区第35| 91aial.com中文字幕在线观看| 婷婷色综合www| 久久精品国产亚洲av涩爱| 国产亚洲欧美精品永久| 欧美成人精品欧美一级黄| 成人毛片60女人毛片免费| 中文字幕av电影在线播放| 两个人免费观看高清视频| 日日撸夜夜添| 亚洲三区欧美一区| 亚洲欧美一区二区三区黑人| 亚洲精品av麻豆狂野| 两个人免费观看高清视频| 99久久精品国产亚洲精品| 飞空精品影院首页| 永久免费av网站大全| 丁香六月欧美| 午夜免费观看性视频| 日本vs欧美在线观看视频| 欧美日韩av久久| 久久女婷五月综合色啪小说| 国产一区二区 视频在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲熟女精品中文字幕| 9热在线视频观看99| av有码第一页| 中文字幕最新亚洲高清| 国产精品麻豆人妻色哟哟久久| 国产 一区精品| 最近的中文字幕免费完整| 免费观看人在逋| 卡戴珊不雅视频在线播放| 国产成人一区二区在线| 中文精品一卡2卡3卡4更新| 男人舔女人的私密视频| 在线观看免费日韩欧美大片| 亚洲国产精品成人久久小说| 久久国产精品大桥未久av| 国产精品免费大片| 欧美久久黑人一区二区| 波多野结衣av一区二区av| 成人黄色视频免费在线看| 亚洲精品日本国产第一区| 国产在线免费精品| 中国国产av一级| 秋霞伦理黄片| 精品午夜福利在线看| 日日啪夜夜爽| 欧美日韩成人在线一区二区| 啦啦啦中文免费视频观看日本| 少妇精品久久久久久久| 亚洲天堂av无毛| 一边摸一边抽搐一进一出视频| 免费久久久久久久精品成人欧美视频| 亚洲国产成人一精品久久久| 日韩一本色道免费dvd| 国产成人免费无遮挡视频| 欧美97在线视频| 午夜老司机福利片| 久久这里只有精品19| 久久精品久久久久久久性| 欧美亚洲 丝袜 人妻 在线| 亚洲av电影在线进入| 精品少妇内射三级| 精品国产露脸久久av麻豆| 亚洲欧美清纯卡通| 亚洲av欧美aⅴ国产| 国产黄色视频一区二区在线观看| 亚洲人成电影观看| 热re99久久精品国产66热6| 制服诱惑二区| 熟妇人妻不卡中文字幕| 99热网站在线观看| 久久国产精品男人的天堂亚洲| 又粗又硬又长又爽又黄的视频| 久久精品熟女亚洲av麻豆精品| 男女边吃奶边做爰视频| 亚洲视频免费观看视频| 久久97久久精品| 久久国产精品大桥未久av| 亚洲综合精品二区| 男女下面插进去视频免费观看| 黄色视频在线播放观看不卡| 日韩伦理黄色片| 自线自在国产av| 操美女的视频在线观看| 精品视频人人做人人爽| 亚洲综合色网址| 人妻 亚洲 视频| 人人澡人人妻人| 在线天堂中文资源库| 日日摸夜夜添夜夜爱| 亚洲欧美成人精品一区二区| 黑人猛操日本美女一级片| 热99久久久久精品小说推荐| 爱豆传媒免费全集在线观看| av线在线观看网站| 精品一区二区免费观看| 高清欧美精品videossex| 伊人亚洲综合成人网| 麻豆乱淫一区二区| 国产伦理片在线播放av一区| xxxhd国产人妻xxx| 在线天堂最新版资源| 亚洲一区中文字幕在线| 亚洲精品久久久久久婷婷小说| 欧美 日韩 精品 国产| 日韩精品有码人妻一区| av在线老鸭窝| 男女高潮啪啪啪动态图| 国产探花极品一区二区| 美国免费a级毛片| 国产免费福利视频在线观看| netflix在线观看网站| 亚洲欧美精品自产自拍| 少妇 在线观看| 欧美日韩视频精品一区| 精品亚洲成国产av| 2018国产大陆天天弄谢| 午夜福利,免费看| 中文字幕av电影在线播放| 肉色欧美久久久久久久蜜桃| 国产精品免费视频内射| av免费观看日本| 亚洲av在线观看美女高潮| 777米奇影视久久| 亚洲视频免费观看视频| 亚洲精品av麻豆狂野| 精品国产露脸久久av麻豆| 国产精品久久久久久精品电影小说| 国产极品天堂在线| 天天躁夜夜躁狠狠躁躁| 狂野欧美激情性xxxx| 18禁国产床啪视频网站| 啦啦啦视频在线资源免费观看| 精品国产露脸久久av麻豆| 一二三四在线观看免费中文在| 日本av免费视频播放| 久久精品国产a三级三级三级| 亚洲综合色网址| 不卡视频在线观看欧美| 老司机影院成人| 日本色播在线视频| 在线观看一区二区三区激情| 亚洲情色 制服丝袜| 国产熟女欧美一区二区| 日韩一区二区三区影片| 色吧在线观看| 亚洲成av片中文字幕在线观看| 高清欧美精品videossex| 中文欧美无线码| 国产一区有黄有色的免费视频| 久久综合国产亚洲精品| 天堂中文最新版在线下载| 婷婷色av中文字幕| 亚洲av福利一区| 国产探花极品一区二区| 在线 av 中文字幕| 操出白浆在线播放| 午夜免费观看性视频| 毛片一级片免费看久久久久| 亚洲精品自拍成人| 女人高潮潮喷娇喘18禁视频| 亚洲欧洲精品一区二区精品久久久 | 美女脱内裤让男人舔精品视频| 亚洲欧美激情在线| 免费不卡黄色视频| 国产毛片在线视频| 亚洲四区av| 欧美人与性动交α欧美精品济南到| 天天躁狠狠躁夜夜躁狠狠躁| 精品少妇久久久久久888优播| 国产免费福利视频在线观看| 女的被弄到高潮叫床怎么办| 国产亚洲午夜精品一区二区久久| 久久女婷五月综合色啪小说| 国产乱来视频区| 一本色道久久久久久精品综合| 免费在线观看完整版高清| 欧美人与性动交α欧美软件| 国产 一区精品| 久久狼人影院| 嫩草影院入口| 国产黄色免费在线视频| 日本91视频免费播放| 日韩成人av中文字幕在线观看| 亚洲av日韩在线播放| 18禁观看日本| 搡老乐熟女国产| 中文字幕人妻丝袜制服| 悠悠久久av| 国产精品欧美亚洲77777| 天天添夜夜摸| 一区二区三区激情视频| 在线观看三级黄色| 看十八女毛片水多多多| 日韩一区二区三区影片| 国产一区亚洲一区在线观看| 日韩,欧美,国产一区二区三区| 午夜激情av网站| 99热网站在线观看| 十八禁网站网址无遮挡| 国产精品一区二区精品视频观看| 国产成人91sexporn| 亚洲,欧美精品.| 老司机影院成人| 亚洲美女视频黄频| 久久亚洲国产成人精品v| 国产精品久久久久成人av| 波多野结衣一区麻豆| 久久精品aⅴ一区二区三区四区| av国产久精品久网站免费入址| 在线观看免费日韩欧美大片| 欧美97在线视频| 大片免费播放器 马上看| 街头女战士在线观看网站| avwww免费| 咕卡用的链子| 制服诱惑二区| av视频免费观看在线观看| 一本久久精品| 精品国产乱码久久久久久小说| 热99久久久久精品小说推荐| 亚洲四区av| 少妇人妻 视频| 一本一本久久a久久精品综合妖精| 老司机亚洲免费影院| 又大又黄又爽视频免费| 黄片小视频在线播放| a级毛片黄视频| e午夜精品久久久久久久| 在线天堂中文资源库| 欧美成人午夜精品| 一级片'在线观看视频| 国产高清不卡午夜福利|