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

    固體推進(jìn)劑/襯層粘接界面脫粘失效的數(shù)值模擬

    2020-11-10 01:32:00張志成戴開(kāi)達(dá)陳靜靜
    科學(xué)技術(shù)與工程 2020年28期
    關(guān)鍵詞:襯層細(xì)觀推進(jìn)劑

    張志成, 戴開(kāi)達(dá)*, 陳靜靜

    (1.北京理工大學(xué)機(jī)電學(xué)院, 北京 100081; 2.上海航天動(dòng)力技術(shù)研究所, 上海 201109)

    固體火箭發(fā)動(dòng)機(jī)由于結(jié)構(gòu)簡(jiǎn)單、可靠性高、儲(chǔ)存維護(hù)方便、作戰(zhàn)反應(yīng)快、加速性能好的特點(diǎn),在國(guó)防事業(yè)中備受重視,目前已廣泛應(yīng)用在各類戰(zhàn)術(shù)、戰(zhàn)略導(dǎo)彈及火箭彈中。固體推進(jìn)劑是固體火箭發(fā)動(dòng)機(jī)的能源,發(fā)動(dòng)機(jī)的重要組成部分。固體推進(jìn)劑/襯層粘接界面的損傷研究很大程度上是圍繞載荷和環(huán)境以及相關(guān)性能展開(kāi)的,研究?jī)?nèi)容涉及損傷的表征、實(shí)驗(yàn)?zāi)M及觀測(cè)方法、損傷對(duì)脫粘的影響和損傷本構(gòu)關(guān)系等各個(gè)方面[1-4]。粘接失效是固體火箭發(fā)動(dòng)機(jī)最常見(jiàn)的失效方式之一,其中固體推進(jìn)劑/襯層粘接界面的脫粘是最主要的失效模式。

    眾多學(xué)者對(duì)發(fā)動(dòng)機(jī)界面粘接強(qiáng)度、界面斷裂韌性及失效過(guò)程進(jìn)行了大量的定性研究分析。Liu[5]利用掃描電鏡觀察發(fā)現(xiàn)ANB3066推進(jìn)劑與襯層界面在緩慢拉伸條件下出現(xiàn)了初始損傷區(qū),該區(qū)域隨著拉伸繼續(xù)沿著裂紋尖端擴(kuò)展。吳豐軍等[6]探索了粘接界面細(xì)觀力學(xué)性能、結(jié)構(gòu)與破壞方式的內(nèi)在關(guān)聯(lián),提出了推進(jìn)劑/襯層粘接界面失效模式。邱欣等[7]利用CCD(charge coupled device)光學(xué)顯微鏡觀察HTPB(hydroxyl-terminated poly butadiene)推進(jìn)劑拉伸斷裂過(guò)程,證明靠近襯層附近的高氯酸銨顆粒與襯層脫濕是影響該推進(jìn)劑粘接性能的主要因素。楊明等[8]采用掃描電鏡原位拉伸試驗(yàn)系統(tǒng)對(duì) HTPB推進(jìn)劑/襯層粘接界面試件拉伸破壞過(guò)程進(jìn)行了觀察,實(shí)時(shí)采集了界面變形破壞過(guò)程的SEM圖像,結(jié)合粘接界面的宏觀應(yīng)力-應(yīng)變曲線,分析了其在拉伸過(guò)程中細(xì)觀變破壞機(jī)理。

    為了表征粘接界面力學(xué)特性和對(duì)粘接界面失效進(jìn)行數(shù)值模擬,大量的研究集中在界面力學(xué)模型中的內(nèi)聚力模型[9-12],王廣等[13]基于復(fù)合固體推進(jìn)劑/襯層粘接試件微CT(computed tomography)掃描后重構(gòu)的圖像統(tǒng)計(jì),得到細(xì)觀形貌特征及推進(jìn)劑顆粒典型尺寸,以最小代表體積元為基礎(chǔ)構(gòu)建了二維細(xì)觀數(shù)值模型,使用粘接界面單元表征推進(jìn)劑/襯層界面,分析了不同老化時(shí)間復(fù)合固體推進(jìn)劑/襯層在單向拉伸作用下的脫粘過(guò)程。李高春等[17]依據(jù)粘接界面細(xì)觀結(jié)構(gòu),建立了界面的細(xì)觀模型,并在顆粒/基體以及推進(jìn)劑/襯層之間引入界面元,對(duì)界面細(xì)觀變形和破壞過(guò)程進(jìn)行了數(shù)值模擬研究。目前,研究人員普遍使用分子動(dòng)力學(xué)算法,根據(jù)推進(jìn)劑級(jí)配關(guān)系生成固體推進(jìn)劑隨機(jī)球形顆粒填充的細(xì)觀模型,基于SEM重構(gòu)的細(xì)觀模型的報(bào)道并不多見(jiàn)。

    本文基于SEM圖片重構(gòu)建立了固體推進(jìn)劑/襯層的細(xì)觀結(jié)構(gòu)有限元模型,進(jìn)而研究了界面上不同顆粒含量以及界面不同強(qiáng)化/老化情況下的界面脫粘過(guò)程。

    1 固體推進(jìn)劑/襯層計(jì)算模型

    為了從細(xì)觀角度真實(shí)地模擬出推進(jìn)劑/襯層界面脫粘擴(kuò)展過(guò)程,并考慮界面處顆粒對(duì)推進(jìn)劑/襯層界面脫粘的影響,利用掃描電鏡得到復(fù)合固體推進(jìn)劑的SEM(scanning electron microscope)圖如圖1所示。將圖1導(dǎo)入Auto CAD軟件進(jìn)行圖像處理,然后導(dǎo)入ABAQUS生成有限元模型,見(jiàn)圖2。其中推進(jìn)劑包含顆粒和基體,尺寸為558 μm×380 μm,襯層厚為77 μm。計(jì)算模型采用準(zhǔn)靜態(tài)位移加載,加載速率為0.5 μm/s。模型左右兩側(cè)的邊界為自由邊界,底部約束y方向的位移、邊界條件如圖2所示。網(wǎng)格劃分采用三節(jié)點(diǎn)平面應(yīng)力單元,總單元數(shù)為67 126。

    圖1 推進(jìn)劑掃描電鏡圖Fig.1 Propellent SEM diagram

    圖2 推進(jìn)劑/襯層數(shù)值模型Fig.2 Propellent/liner numerical model

    圖3 顆粒/基體界面和推進(jìn)劑襯層界面Fig.3 Particle/binder interface and propellant/liner interface

    圖4 Cohesive單元模型Fig.4 Cohesive element model

    ABAQUS內(nèi)置了Cohesive單元來(lái)表征界面分離時(shí)力和位移的關(guān)系,分別在顆粒和基體、推進(jìn)劑和襯層之間插入Cohesive單元,如圖3所示。圖4為Cohesive單元模型,1、2、3、4為三角形單元的節(jié)點(diǎn)編號(hào),1、2、3號(hào)節(jié)點(diǎn)組成了顆粒單元,2、3、4號(hào)節(jié)點(diǎn)組成了基體單元。在未插入Cohesive單元之前,基體單元和顆粒單元有兩個(gè)共同的節(jié)點(diǎn)2和3,離散這兩個(gè)單元并插入新的兩個(gè)節(jié)點(diǎn),便形成了由2、3、5、6四個(gè)節(jié)點(diǎn)組成的Cohesive單元。在初始時(shí)刻,Cohesive單元的厚度為0,即單元的上下表面的坐標(biāo)相同,受到拉伸載荷上下表面分離。

    2 模型參數(shù)

    固體推進(jìn)劑的力學(xué)性能由基體、顆粒、襯層和粘接界面的參數(shù)共同決定。基體部分為黏彈性材料,其力學(xué)性能(如剪切模量G、體積模量K)隨著時(shí)間τ和溫度的改變而改變。用ABAQUS內(nèi)置的Prony級(jí)數(shù)來(lái)表示基體的黏彈性,參數(shù)如表1所示[14]?;w的彈性模量E為8 MPa,泊松比υ為0.499[15]。

    模型中的顆粒為AP(ammonium perchlorate),由于顆粒的彈性模量遠(yuǎn)大于基體,將顆粒視為彈性體,具體參數(shù)見(jiàn)表2[13]。

    表1 基體的Prony級(jí)數(shù)

    表2 顆粒參數(shù)

    ABAQUS中的粘接單元是通過(guò)黏聚裂紋模型簡(jiǎn)化而來(lái)的,圖5是雙線性粘接裂紋模型,在該模型中,在初始階段隨著張開(kāi)位移的增加,裂紋表面黏結(jié)力將逐漸增加;當(dāng)張開(kāi)位移達(dá)到一定值時(shí),裂紋面間的粘接力達(dá)到最大值;隨著裂紋表面間張開(kāi)位移的繼續(xù)增加,裂紋面間黏結(jié)力將逐漸降低到零,此時(shí)裂紋面最終形成。在黏結(jié)力達(dá)到開(kāi)裂強(qiáng)度之后,材料表現(xiàn)為線性軟化行為,是材料的軟化階段,也稱為損傷階段。

    圖5 雙線性粘接裂紋模型Fig.5 Bilinear bond crack model

    界面的本構(gòu)關(guān)系式為

    (1)

    (2)

    損傷起始就意味著材料開(kāi)始軟化。當(dāng)應(yīng)力滿足最大名義應(yīng)力準(zhǔn)則損傷起始準(zhǔn)則的時(shí)候,開(kāi)始出現(xiàn)損傷。

    (3)

    式(3)中:σ為名義應(yīng)力;Nmax、Smax和Tmax分別代表當(dāng)變形是法向(對(duì)應(yīng)Ⅰ型裂紋)、剪切1方向(對(duì)應(yīng)Ⅱ型裂紋)以及剪切2方向(對(duì)應(yīng)Ⅲ型裂紋)時(shí)名義應(yīng)力的峰值。

    當(dāng)滿足損傷起始準(zhǔn)則后,單元發(fā)生損傷,即進(jìn)入損傷演化階段。常用的損傷演化規(guī)律包括線性衰減演化與指數(shù)衰減演化,采用線性衰減演化。Cohesive單元的參數(shù)有三個(gè),分別是拉伸強(qiáng)度σ、失效位移δ和初始剛度k,具體參數(shù)見(jiàn)表3[16]。

    表3 Cohesive單元參數(shù)

    3 模擬結(jié)果分析

    3.1 界面顆粒含量的影響

    在推進(jìn)劑/襯層界面含有小顆粒,為了研究不同推進(jìn)劑/襯層界面的顆粒含量對(duì)界面失效過(guò)程的影響,在模型1的基礎(chǔ)上增加和減少了界面的顆粒尺寸含量,但不改變推進(jìn)劑/襯層界面的強(qiáng)度,如圖6所示。界面的顆粒含量Δ用顆粒/襯層界面長(zhǎng)度與推進(jìn)劑/襯層界面長(zhǎng)度的比值表示,具體如表4所示。

    圖6 推進(jìn)劑/襯層界面顆粒分布Fig.6 Particle distribution of matrix/liner interface

    表4 推進(jìn)劑/襯層界面顆粒含量

    圖7為不同時(shí)刻模型1(Δ=17.29%)的SDEG(scalar stiffness degradation)云圖,當(dāng)Cohesive單元的SDEG=1時(shí)表示該單元完全失效。由圖7(a)可以看出,在拉伸載荷作用下推進(jìn)劑內(nèi)部的大顆粒表面和襯層表面都出現(xiàn)了損傷演化,大顆粒首先脫濕產(chǎn)生裂紋,襯層界面損傷趨于臨界值,小顆粒與基體粘接完好。圖7(b)表明當(dāng)應(yīng)變?cè)黾拥?1.6%,大顆粒完全脫粘,由于裂紋尖端的應(yīng)力集中導(dǎo)致大顆粒附近的顆粒也出現(xiàn)了裂紋,但襯層沒(méi)有產(chǎn)生裂紋。由圖7(c)可知當(dāng)應(yīng)變達(dá)到24.77%時(shí),推進(jìn)劑內(nèi)部產(chǎn)生三條和拉伸載荷方向垂直的裂紋,襯層界面處有一個(gè)顆粒與基體脫粘,由于顆粒的脫粘導(dǎo)致該顆粒附近的襯層/基體界面應(yīng)力集中,損傷變量達(dá)到臨界值,此時(shí)模型處于軟化階段,所以隨著位移載荷繼續(xù)增加該推進(jìn)劑/襯層界面SDEG逐漸降低。

    圖7 模型1不同應(yīng)變時(shí)的SDEG云圖Fig.7 Model 1 of different strain SDEG cloud

    圖8為模型2(Δ=1.43%)應(yīng)變?yōu)?2.91%時(shí)的SDEG云圖。模型2在單軸拉伸下的損傷演化趨勢(shì)與模型1基本一致,大顆粒附近首先出現(xiàn)脫濕,襯層損傷有所演化但是沒(méi)有脫粘。在應(yīng)變達(dá)到22.91%時(shí),模型處于軟化階段,推進(jìn)劑內(nèi)部有兩條與拉伸載荷方向垂直的裂紋,右側(cè)裂紋是圖7(c)的中間和右側(cè)裂紋貫通而成。和模型1相比較,軟化階段模型1襯層處發(fā)現(xiàn)顆粒的脫濕[見(jiàn)圖7(c)箭頭所指],但在模型2中沒(méi)有發(fā)生顆粒脫濕。

    圖8 模型2應(yīng)變?yōu)?2.91%時(shí)的SDEG云圖Fig.8 Model 2 strain was 22.91% of SDEG cloud

    圖9 模型3不同應(yīng)變時(shí)的SDEG云圖Fig.9 Model 3 of different strain SDEG cloud

    圖9為模型3(Δ=38.34%)不同時(shí)刻的SDEG云圖,圖9(a)表明應(yīng)變?yōu)?2.91%時(shí),達(dá)到該模型的應(yīng)力峰值,大顆粒界面裂紋有所擴(kuò)展,與模型1和模型2不同的是大顆粒界面沒(méi)有完全脫濕,大顆粒附近的顆粒也沒(méi)有出現(xiàn)裂紋,而襯層附近的顆?;景l(fā)生脫粘。圖9(b)表明應(yīng)變?yōu)?6.4%時(shí),推進(jìn)劑內(nèi)部有三條與單軸拉伸方向垂直的裂紋,但是擴(kuò)展寬度與模型1相比較小。襯層上的顆粒與基體完全脫粘,襯層與基體之間的界面達(dá)到損傷臨界值,位移進(jìn)一步增加后推進(jìn)劑與襯層完全斷開(kāi)。

    圖10 模型整體的平均應(yīng)力-應(yīng)變曲線Fig.10 The average stress-strain curve of the whole model

    圖10為模型1、模型2和模型3的整體平均應(yīng)力-應(yīng)變曲線,可以看到,隨著襯層界面顆粒含量的增加,達(dá)到應(yīng)力峰值需要的應(yīng)變?cè)龃螅@是由顆粒界面的Cohesive單元在達(dá)到應(yīng)力峰值時(shí)貢獻(xiàn)的位移提高導(dǎo)致的。結(jié)合圖7(b)和圖9(a)可以發(fā)現(xiàn),模型1和模型2在達(dá)到應(yīng)力峰值時(shí),僅有推進(jìn)劑內(nèi)部顆粒的脫濕,導(dǎo)致了模型的軟化,而模型3達(dá)到應(yīng)力峰值時(shí),推進(jìn)劑內(nèi)部顆粒的脫濕程度減小,襯層界面的顆粒完全脫粘,導(dǎo)致了模型的軟化。達(dá)到峰值時(shí)的應(yīng)力最大值沒(méi)有呈現(xiàn)一定的趨勢(shì),在Δ=17.29%時(shí),應(yīng)力減小,Δ=38.34%時(shí)應(yīng)力增加。

    3.2 推進(jìn)劑/襯層界面強(qiáng)度的影響

    襯層的老化會(huì)引起界面模量的降低,改變界面所用的黏結(jié)劑也會(huì)引起界面模量的變化。為了研究不同推進(jìn)劑/襯層界面模量對(duì)界面脫粘力學(xué)響應(yīng)的影響,在初始模型的基礎(chǔ)上增強(qiáng)和降低了推進(jìn)劑/襯層界面強(qiáng)度,具體參數(shù)如表5所示。

    表5 不同推進(jìn)劑/襯層界面強(qiáng)度參數(shù)

    圖11為不同強(qiáng)度的推進(jìn)劑/襯層界面在準(zhǔn)靜態(tài)拉伸載荷作用下的應(yīng)力應(yīng)變曲線,初始模型在應(yīng)變?yōu)?1.6%時(shí)達(dá)到應(yīng)力峰值,老化模型1和老化模型2均在應(yīng)變?yōu)?0.2%時(shí)達(dá)到應(yīng)力峰值,隨著老化程度的提高,應(yīng)力峰值不斷減小。強(qiáng)化模型1在應(yīng)變?yōu)?3.46%時(shí)達(dá)到應(yīng)力峰值,強(qiáng)化模型2和強(qiáng)化模型3在應(yīng)變?yōu)?4%時(shí)達(dá)到應(yīng)力峰值,且峰值相差不大。強(qiáng)化模型2和強(qiáng)化模型3的曲線差異不大,在不發(fā)生老化的情況下,推進(jìn)劑/襯層界面強(qiáng)度的提升無(wú)法提高模型的強(qiáng)度。

    圖11 不同強(qiáng)度的推進(jìn)劑/襯層界面在準(zhǔn)靜態(tài)拉伸載荷作用下的應(yīng)力應(yīng)變曲線Fig.11 Stress-strain curves of propellant/liner interfaces of different intensities under quasi-static tensile loading

    圖12(a)為老化模型1(σ=1 MPa,E=150 MPa)達(dá)到應(yīng)力峰值時(shí)的SDEG云圖。由于推進(jìn)劑/襯層界面強(qiáng)度的降低,襯層界面的顆粒與基體開(kāi)始脫粘,推進(jìn)劑內(nèi)部有大顆粒的脫濕。圖12(b)為老化模型2(σ=0.8 MPa,E=120 MPa)達(dá)到應(yīng)力峰值時(shí)的SDEG云圖。由于推進(jìn)劑/襯層界面強(qiáng)度的進(jìn)一步降低,此時(shí)推進(jìn)劑/襯層界面損傷都開(kāi)始演化,推進(jìn)劑內(nèi)部顆粒的損傷演化減小。

    圖12 老化模型達(dá)到應(yīng)力峰值的SDEG分布Fig.12 The SDEG distribution of aging model at the peak stress

    圖13(a)為強(qiáng)化模型1(σ=1.5 MPa,E=300 MPa),達(dá)到應(yīng)力峰值時(shí)推進(jìn)劑內(nèi)部出現(xiàn)多處顆粒的脫濕,襯層界面的顆粒與基體脫粘。圖13(b)為強(qiáng)化模型2(σ=2 MPa,E=500 MPa),此時(shí)僅有推進(jìn)劑內(nèi)部顆粒有多處脫濕,襯層界面及界面上的顆粒損傷演化程度很小。

    圖13 強(qiáng)化模型達(dá)到應(yīng)力峰值時(shí)的SDEG分布Fig.13 The SDEG distribution of strengthening model at the peak stress

    圖13(c)為強(qiáng)化模型3(σ=3 MPa,E=1 000 MPa),損傷變量分布與強(qiáng)化模型2類似,僅有推進(jìn)劑內(nèi)部顆粒的脫濕,可見(jiàn)再增強(qiáng)推進(jìn)劑/襯層界面強(qiáng)度已經(jīng)對(duì)模型達(dá)到應(yīng)力峰值時(shí)的損傷分布沒(méi)有影響。

    4 結(jié)論

    基于細(xì)觀有限元模型,結(jié)合Cohesive單元模擬了推進(jìn)劑/襯層界面在單軸拉伸下的脫粘失效過(guò)程,計(jì)算得到應(yīng)力應(yīng)變曲線,取得如下結(jié)論。

    (1)隨著載荷的增加,首先在推進(jìn)劑內(nèi)部的大顆粒附近界面由于應(yīng)力集中首先發(fā)生損傷,隨后大顆粒開(kāi)始脫濕并引起周圍顆粒的脫濕形成裂紋,模型開(kāi)始軟化。

    (2)襯層界面處顆粒含量的增加不能引起模型應(yīng)力峰值的線性改變。當(dāng)模型達(dá)到應(yīng)力峰值時(shí)若僅有推進(jìn)劑內(nèi)部顆粒的裂紋擴(kuò)展,顆粒含量增加會(huì)導(dǎo)致應(yīng)力峰值減小。若此時(shí)襯層處顆粒同樣產(chǎn)生裂紋,則會(huì)使應(yīng)力集中在襯層界面,由于襯層/基體界面強(qiáng)度更高,模型的應(yīng)力峰值會(huì)相應(yīng)提高。

    (3)推進(jìn)劑/襯層界面的強(qiáng)度的提高無(wú)法使該模型獲得強(qiáng)度更高的應(yīng)力應(yīng)變曲線,界面本身的強(qiáng)度已經(jīng)使推進(jìn)劑/襯層模型在受到單軸拉伸載荷時(shí)破壞的主要形式為推進(jìn)劑內(nèi)部顆粒的脫粘。推進(jìn)劑/襯層界面的老化使模型在承受載荷時(shí)以推進(jìn)劑內(nèi)部顆粒脫粘失效為主逐漸轉(zhuǎn)變?yōu)橥七M(jìn)劑/襯層界面的脫粘失效,應(yīng)力峰值降低,并且達(dá)到應(yīng)力峰值需要的應(yīng)變也降低。

    猜你喜歡
    襯層細(xì)觀推進(jìn)劑
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動(dòng)彈性模量和阻尼比試驗(yàn)研究
    地震研究(2021年1期)2021-04-13 01:05:24
    超磁致伸縮傳感器背襯層參數(shù)優(yōu)化的實(shí)驗(yàn)研究
    電子科技(2021年2期)2021-01-08 02:26:02
    基于四叉樹(shù)網(wǎng)格加密技術(shù)的混凝土細(xì)觀模型
    KNSB推進(jìn)劑最佳配比研究
    固體火箭發(fā)動(dòng)機(jī)襯層固化狀態(tài)超聲波實(shí)時(shí)監(jiān)測(cè)方法
    含LLM-105無(wú)煙CMDB推進(jìn)劑的燃燒性能
    無(wú)鋁低燃速NEPE推進(jìn)劑的燃燒性能
    DNTF-CMDB推進(jìn)劑的燃燒機(jī)理
    PBX炸藥的抗壓強(qiáng)度及抗拉強(qiáng)度細(xì)觀尺度的數(shù)值計(jì)算
    開(kāi)裂混凝土中水分傳輸過(guò)程的細(xì)觀模型
    看非洲黑人一级黄片| 亚洲av.av天堂| 日韩中文字幕视频在线看片| av卡一久久| 国产精品国产三级国产专区5o| 天天影视国产精品| 99精国产麻豆久久婷婷| 亚洲综合色惰| 国产综合精华液| 国产视频首页在线观看| 久久精品熟女亚洲av麻豆精品| 午夜激情av网站| 一区二区三区免费毛片| 国产亚洲精品久久久com| 在现免费观看毛片| 亚洲欧美日韩另类电影网站| 99热6这里只有精品| 日本91视频免费播放| 色94色欧美一区二区| 少妇被粗大猛烈的视频| 天美传媒精品一区二区| 婷婷色av中文字幕| 看免费成人av毛片| 亚洲精品日韩av片在线观看| av在线播放精品| 亚洲精品中文字幕在线视频| 精品人妻熟女毛片av久久网站| 18在线观看网站| 亚洲欧洲日产国产| a级毛色黄片| 少妇的逼水好多| 我的老师免费观看完整版| 女人精品久久久久毛片| 大片免费播放器 马上看| 亚洲av不卡在线观看| 久久av网站| 汤姆久久久久久久影院中文字幕| 国产精品女同一区二区软件| 三级国产精品片| 国产av一区二区精品久久| 我的老师免费观看完整版| 97超视频在线观看视频| 亚洲国产最新在线播放| 天天影视国产精品| 边亲边吃奶的免费视频| 日韩不卡一区二区三区视频在线| 亚洲av日韩在线播放| 在现免费观看毛片| a 毛片基地| 午夜激情福利司机影院| 女性被躁到高潮视频| 亚洲四区av| 韩国av在线不卡| 999精品在线视频| 91成人精品电影| 国产不卡av网站在线观看| 成人亚洲精品一区在线观看| 视频在线观看一区二区三区| 精品人妻偷拍中文字幕| 亚洲精品国产色婷婷电影| 少妇人妻精品综合一区二区| 亚洲欧美成人精品一区二区| 91国产中文字幕| 18禁在线播放成人免费| 久热久热在线精品观看| 夫妻性生交免费视频一级片| 午夜日本视频在线| 91国产中文字幕| 亚洲一级一片aⅴ在线观看| 激情五月婷婷亚洲| 观看美女的网站| 美女内射精品一级片tv| 国产 精品1| 亚洲美女搞黄在线观看| 免费看不卡的av| 91久久精品电影网| 国产精品秋霞免费鲁丝片| 亚洲,一卡二卡三卡| 18禁动态无遮挡网站| 看十八女毛片水多多多| 日本wwww免费看| 亚洲国产精品一区三区| 国内精品宾馆在线| 免费人成在线观看视频色| 国产日韩欧美亚洲二区| 在线播放无遮挡| 久久人妻熟女aⅴ| 欧美日韩精品成人综合77777| √禁漫天堂资源中文www| 男男h啪啪无遮挡| 免费黄色在线免费观看| 狠狠婷婷综合久久久久久88av| 日本免费在线观看一区| 欧美变态另类bdsm刘玥| 男女啪啪激烈高潮av片| 日本av手机在线免费观看| 国产精品久久久久久精品古装| 精品人妻熟女毛片av久久网站| 99热国产这里只有精品6| 毛片一级片免费看久久久久| 少妇精品久久久久久久| 2022亚洲国产成人精品| 日日摸夜夜添夜夜爱| 亚洲精品乱码久久久v下载方式| 久久国产精品大桥未久av| 亚洲国产精品专区欧美| 人妻人人澡人人爽人人| 国产av一区二区精品久久| freevideosex欧美| 天堂8中文在线网| 丝袜在线中文字幕| 九九在线视频观看精品| 啦啦啦在线观看免费高清www| 亚洲国产色片| 久久亚洲国产成人精品v| 成人黄色视频免费在线看| 亚洲欧美日韩卡通动漫| 十八禁高潮呻吟视频| 26uuu在线亚洲综合色| 亚洲国产精品一区三区| 国产欧美日韩一区二区三区在线 | 91国产中文字幕| 亚洲精品久久午夜乱码| 久久久久久久亚洲中文字幕| 久久免费观看电影| 国产黄片视频在线免费观看| 欧美成人午夜免费资源| 纵有疾风起免费观看全集完整版| 亚洲国产精品一区三区| 亚洲精品国产av蜜桃| 亚洲人成77777在线视频| 成人国语在线视频| 亚洲精品成人av观看孕妇| 亚洲在久久综合| 日韩不卡一区二区三区视频在线| a级毛片黄视频| 亚洲精品一二三| 国产午夜精品久久久久久一区二区三区| 久久久久久久精品精品| 精品亚洲成a人片在线观看| 满18在线观看网站| 校园人妻丝袜中文字幕| 高清毛片免费看| 高清午夜精品一区二区三区| 国产高清国产精品国产三级| 亚洲图色成人| 日韩av在线免费看完整版不卡| 午夜免费鲁丝| 国产高清三级在线| 免费高清在线观看视频在线观看| 日韩欧美精品免费久久| 久久毛片免费看一区二区三区| 哪个播放器可以免费观看大片| 精品卡一卡二卡四卡免费| 91精品国产九色| 熟女av电影| 汤姆久久久久久久影院中文字幕| 日本猛色少妇xxxxx猛交久久| 亚洲欧美日韩另类电影网站| 成人亚洲欧美一区二区av| xxx大片免费视频| 欧美日韩av久久| 亚洲国产精品999| 免费播放大片免费观看视频在线观看| 国产 精品1| 免费看光身美女| 亚洲综合色惰| 在线观看www视频免费| 成人亚洲精品一区在线观看| 欧美少妇被猛烈插入视频| 亚洲精品456在线播放app| 大陆偷拍与自拍| 九色成人免费人妻av| 在线观看www视频免费| 老司机影院毛片| 日韩av免费高清视频| 日本黄色片子视频| 精品熟女少妇av免费看| 国产精品女同一区二区软件| 91成人精品电影| 国产黄色免费在线视频| 国精品久久久久久国模美| 丝袜喷水一区| 久久久亚洲精品成人影院| 久久久久久久精品精品| 久久婷婷青草| 国产精品秋霞免费鲁丝片| a 毛片基地| 欧美激情极品国产一区二区三区 | 亚洲欧美清纯卡通| 人妻一区二区av| 亚洲av.av天堂| 99热6这里只有精品| 人人妻人人澡人人爽人人夜夜| 少妇的逼好多水| 亚洲av综合色区一区| 九九在线视频观看精品| 最近中文字幕高清免费大全6| 久久久久久久亚洲中文字幕| 免费高清在线观看视频在线观看| 男女国产视频网站| 日韩一区二区三区影片| 一级毛片黄色毛片免费观看视频| 丰满乱子伦码专区| videos熟女内射| 亚洲精品乱码久久久久久按摩| 91精品国产国语对白视频| 日韩视频在线欧美| 草草在线视频免费看| 简卡轻食公司| 欧美最新免费一区二区三区| 青青草视频在线视频观看| av在线老鸭窝| 五月伊人婷婷丁香| 亚洲欧美日韩另类电影网站| 大陆偷拍与自拍| 一级a做视频免费观看| 曰老女人黄片| 亚洲精品乱码久久久v下载方式| 亚洲丝袜综合中文字幕| 国产熟女欧美一区二区| 亚洲天堂av无毛| 亚洲国产精品一区三区| 99热这里只有是精品在线观看| 欧美变态另类bdsm刘玥| 午夜日本视频在线| 青春草亚洲视频在线观看| 久久久精品区二区三区| 日日摸夜夜添夜夜添av毛片| 只有这里有精品99| 午夜免费鲁丝| 久久久久久久大尺度免费视频| 国产成人freesex在线| av福利片在线| 51国产日韩欧美| 成人亚洲精品一区在线观看| 日韩成人av中文字幕在线观看| 国产一区有黄有色的免费视频| 欧美日韩成人在线一区二区| 色婷婷久久久亚洲欧美| 亚洲成色77777| 久久午夜综合久久蜜桃| av在线老鸭窝| 极品人妻少妇av视频| 99热6这里只有精品| 国产熟女欧美一区二区| 午夜免费鲁丝| a级片在线免费高清观看视频| 成人国产av品久久久| 超色免费av| 高清欧美精品videossex| 国产成人一区二区在线| 国产亚洲最大av| 午夜福利视频在线观看免费| 精品国产一区二区久久| 另类精品久久| 午夜日本视频在线| 下体分泌物呈黄色| 永久网站在线| 人妻人人澡人人爽人人| 亚洲经典国产精华液单| 免费观看性生交大片5| 午夜免费观看性视频| 最新中文字幕久久久久| 亚洲成人一二三区av| 国产女主播在线喷水免费视频网站| 亚洲精品av麻豆狂野| 80岁老熟妇乱子伦牲交| 精品午夜福利在线看| 亚洲国产成人一精品久久久| 久久久午夜欧美精品| 中文字幕精品免费在线观看视频 | 久久精品久久久久久噜噜老黄| 亚洲精品成人av观看孕妇| 亚洲第一区二区三区不卡| av在线播放精品| 亚洲精华国产精华液的使用体验| 五月玫瑰六月丁香| 又黄又爽又刺激的免费视频.| 3wmmmm亚洲av在线观看| 久久99热6这里只有精品| 亚洲国产精品一区二区三区在线| 日韩免费高清中文字幕av| 国产成人freesex在线| 一级毛片电影观看| 男人添女人高潮全过程视频| 五月玫瑰六月丁香| 亚洲欧美色中文字幕在线| 一区二区av电影网| 99热这里只有精品一区| 视频在线观看一区二区三区| 欧美xxxx性猛交bbbb| videosex国产| 国产男人的电影天堂91| 免费看av在线观看网站| videosex国产| 国产亚洲欧美精品永久| 久久久久精品性色| 午夜免费男女啪啪视频观看| 国产精品偷伦视频观看了| 黄色毛片三级朝国网站| 久久久国产欧美日韩av| a级毛片在线看网站| 黄片播放在线免费| 亚洲精品一区蜜桃| 美女大奶头黄色视频| 人人澡人人妻人| 亚洲精品乱码久久久v下载方式| 午夜福利在线观看免费完整高清在| 狂野欧美激情性xxxx在线观看| 91精品国产九色| 精品国产一区二区三区久久久樱花| 久久人人爽人人爽人人片va| 久久久国产欧美日韩av| 国产高清有码在线观看视频| 国产探花极品一区二区| 日本爱情动作片www.在线观看| 国产精品不卡视频一区二区| 久久ye,这里只有精品| 成人综合一区亚洲| 一边亲一边摸免费视频| 中文字幕人妻熟人妻熟丝袜美| 26uuu在线亚洲综合色| 亚洲成人一二三区av| 亚洲欧美精品自产自拍| 国产高清有码在线观看视频| 日本91视频免费播放| 亚洲av中文av极速乱| av在线播放精品| 精品人妻一区二区三区麻豆| 高清黄色对白视频在线免费看| 青春草亚洲视频在线观看| 精品亚洲乱码少妇综合久久| 国产永久视频网站| 亚洲色图综合在线观看| 亚洲欧美色中文字幕在线| 欧美日韩在线观看h| 成人无遮挡网站| 久久av网站| 中国三级夫妇交换| 麻豆乱淫一区二区| 欧美日韩综合久久久久久| 一级毛片电影观看| 自线自在国产av| av网站免费在线观看视频| 最近中文字幕高清免费大全6| 男女无遮挡免费网站观看| 久久精品久久久久久噜噜老黄| 成年人午夜在线观看视频| 免费av不卡在线播放| 91精品伊人久久大香线蕉| av有码第一页| 久久韩国三级中文字幕| 亚洲经典国产精华液单| 性色avwww在线观看| 大陆偷拍与自拍| 国模一区二区三区四区视频| 你懂的网址亚洲精品在线观看| 性色avwww在线观看| 久久久久久久久久人人人人人人| 日韩免费高清中文字幕av| 人体艺术视频欧美日本| 春色校园在线视频观看| 国产免费一级a男人的天堂| 在现免费观看毛片| 久久热精品热| 在线免费观看不下载黄p国产| 亚洲激情五月婷婷啪啪| 青春草国产在线视频| 这个男人来自地球电影免费观看 | 少妇被粗大猛烈的视频| 亚洲图色成人| 97超视频在线观看视频| 亚洲第一av免费看| 免费久久久久久久精品成人欧美视频 | av一本久久久久| 如何舔出高潮| 亚洲人成77777在线视频| 精品久久国产蜜桃| 晚上一个人看的免费电影| 交换朋友夫妻互换小说| 国产黄片视频在线免费观看| 精品一品国产午夜福利视频| 最新中文字幕久久久久| 少妇人妻精品综合一区二区| 久久久久久久久久人人人人人人| 9色porny在线观看| 国产成人精品一,二区| 少妇丰满av| 一边摸一边做爽爽视频免费| 欧美三级亚洲精品| 日日爽夜夜爽网站| 久久久精品免费免费高清| 久久久久久久久大av| 亚洲av综合色区一区| 中文天堂在线官网| 国语对白做爰xxxⅹ性视频网站| 亚洲无线观看免费| 久久ye,这里只有精品| 欧美激情 高清一区二区三区| 亚洲av电影在线观看一区二区三区| 国产成人一区二区在线| 久久青草综合色| 少妇高潮的动态图| 婷婷色综合www| 天堂中文最新版在线下载| 五月天丁香电影| 欧美精品一区二区大全| 麻豆精品久久久久久蜜桃| 午夜av观看不卡| 欧美 日韩 精品 国产| 亚洲av成人精品一二三区| 极品少妇高潮喷水抽搐| 九色成人免费人妻av| 99久久精品一区二区三区| 欧美少妇被猛烈插入视频| 欧美激情 高清一区二区三区| 人人妻人人添人人爽欧美一区卜| 亚洲中文av在线| 精品国产国语对白av| 亚洲第一区二区三区不卡| 欧美日韩av久久| 国产极品天堂在线| 三上悠亚av全集在线观看| 亚洲伊人久久精品综合| 男人添女人高潮全过程视频| 久久影院123| 日韩,欧美,国产一区二区三区| 啦啦啦视频在线资源免费观看| 免费高清在线观看日韩| 欧美 亚洲 国产 日韩一| 久久久久久久久久久久大奶| 超碰97精品在线观看| 3wmmmm亚洲av在线观看| 午夜福利视频在线观看免费| 97在线视频观看| 中国美白少妇内射xxxbb| 国产av国产精品国产| 成人亚洲精品一区在线观看| 美女福利国产在线| 十八禁网站网址无遮挡| 日日爽夜夜爽网站| 亚洲第一区二区三区不卡| 能在线免费看毛片的网站| www.av在线官网国产| 一区二区三区四区激情视频| 九草在线视频观看| 黄色怎么调成土黄色| 国产成人一区二区在线| 狂野欧美激情性xxxx在线观看| 精品一品国产午夜福利视频| 青青草视频在线视频观看| 国产白丝娇喘喷水9色精品| 亚洲精品国产av成人精品| 日本-黄色视频高清免费观看| 国产成人aa在线观看| 如日韩欧美国产精品一区二区三区 | 久久久久久久久久成人| 精品国产国语对白av| 高清黄色对白视频在线免费看| 一级黄片播放器| 欧美bdsm另类| 人人妻人人添人人爽欧美一区卜| 国产探花极品一区二区| 午夜福利影视在线免费观看| 欧美精品一区二区大全| 亚洲成人一二三区av| 菩萨蛮人人尽说江南好唐韦庄| 黄色欧美视频在线观看| 天堂8中文在线网| 夫妻性生交免费视频一级片| 国产 精品1| 亚洲色图 男人天堂 中文字幕 | h视频一区二区三区| 免费大片18禁| 99国产综合亚洲精品| 美女国产高潮福利片在线看| av播播在线观看一区| 如何舔出高潮| 大码成人一级视频| a级片在线免费高清观看视频| 国产免费视频播放在线视频| 人妻夜夜爽99麻豆av| 麻豆精品久久久久久蜜桃| 成人国产麻豆网| 精品一区二区三卡| 亚洲精品久久午夜乱码| 日本-黄色视频高清免费观看| 99热全是精品| 女的被弄到高潮叫床怎么办| 成人国产av品久久久| 午夜激情久久久久久久| 亚洲在久久综合| 插阴视频在线观看视频| 久久久午夜欧美精品| 亚洲美女搞黄在线观看| 国产男女内射视频| 伊人久久精品亚洲午夜| 亚洲av免费高清在线观看| freevideosex欧美| 亚洲精品一区蜜桃| 91久久精品国产一区二区三区| 久久av网站| 成人无遮挡网站| 亚洲精品亚洲一区二区| 飞空精品影院首页| 最新中文字幕久久久久| 99久久人妻综合| 成人国产av品久久久| 国产精品一区二区在线观看99| 最近手机中文字幕大全| 国产极品天堂在线| 极品少妇高潮喷水抽搐| 中国国产av一级| 大陆偷拍与自拍| 少妇熟女欧美另类| 国产一区二区在线观看av| 精品人妻在线不人妻| 久久国产亚洲av麻豆专区| 日本黄大片高清| 亚洲精品中文字幕在线视频| 国产精品人妻久久久影院| 日日摸夜夜添夜夜添av毛片| 国产在视频线精品| 免费人妻精品一区二区三区视频| 国产乱人偷精品视频| 久久ye,这里只有精品| .国产精品久久| 少妇人妻精品综合一区二区| www.色视频.com| av视频免费观看在线观看| 精品少妇内射三级| 久久久午夜欧美精品| av在线app专区| 成人18禁高潮啪啪吃奶动态图 | 91久久精品国产一区二区成人| 亚洲国产精品一区二区三区在线| 欧美 日韩 精品 国产| 大片电影免费在线观看免费| 日本黄色日本黄色录像| 国产精品嫩草影院av在线观看| 免费观看a级毛片全部| 日本欧美视频一区| 欧美性感艳星| 国产在线免费精品| 一区二区三区乱码不卡18| 在线看a的网站| 亚洲欧美日韩另类电影网站| 国产精品偷伦视频观看了| 久久久久久久国产电影| 日韩大片免费观看网站| 一级毛片aaaaaa免费看小| 男女国产视频网站| 精品人妻熟女毛片av久久网站| 91精品国产国语对白视频| 黑丝袜美女国产一区| 婷婷色综合www| 亚洲av日韩在线播放| av女优亚洲男人天堂| 亚洲精品亚洲一区二区| 国产亚洲一区二区精品| 国产黄频视频在线观看| 亚洲精品456在线播放app| 最新中文字幕久久久久| 十八禁高潮呻吟视频| 亚洲国产欧美在线一区| 一级毛片aaaaaa免费看小| 亚洲欧洲精品一区二区精品久久久 | 免费看av在线观看网站| 精品久久久久久久久亚洲| 亚洲图色成人| 黑人高潮一二区| 99久久精品一区二区三区| 亚洲国产精品一区三区| 久久鲁丝午夜福利片| 欧美日韩视频高清一区二区三区二| 伦精品一区二区三区| 久久综合国产亚洲精品| 夫妻性生交免费视频一级片| 日本黄色片子视频| 午夜激情久久久久久久| 国产白丝娇喘喷水9色精品| 久久精品国产亚洲网站| av又黄又爽大尺度在线免费看| 精品人妻熟女毛片av久久网站| 国产精品久久久久久精品电影小说| 搡女人真爽免费视频火全软件| 少妇的逼水好多| 婷婷色综合www| av福利片在线| 18+在线观看网站| 欧美 亚洲 国产 日韩一| 99热国产这里只有精品6| 国产精品人妻久久久久久| 婷婷色综合www| 亚洲精品久久久久久婷婷小说| 三上悠亚av全集在线观看| 国产熟女午夜一区二区三区 | 97在线人人人人妻| 国产亚洲欧美精品永久| 欧美老熟妇乱子伦牲交| 热re99久久国产66热| 91久久精品电影网| 久久久久久久大尺度免费视频| 日韩大片免费观看网站| 少妇被粗大的猛进出69影院 | 久久ye,这里只有精品| 国产精品一区二区在线观看99| 午夜老司机福利剧场| 国产极品天堂在线| 国产在线一区二区三区精| a级毛片免费高清观看在线播放| videosex国产| 国产精品99久久99久久久不卡 | 一本久久精品|