• <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ì)觀模型
    欧美性感艳星| 久热这里只有精品99| 国产成人一区二区在线| 男的添女的下面高潮视频| 久久久成人免费电影| 亚洲国产高清在线一区二区三| 亚洲第一区二区三区不卡| 国产午夜精品久久久久久一区二区三区| 亚洲国产av新网站| 国产免费视频播放在线视频| 最近的中文字幕免费完整| 九色成人免费人妻av| 久久久久久久久久久丰满| 一区二区av电影网| 久久久精品免费免费高清| 久久这里有精品视频免费| 日韩欧美一区视频在线观看 | 制服丝袜香蕉在线| 国产精品不卡视频一区二区| 成人午夜精彩视频在线观看| 成人欧美大片| 久久精品国产亚洲av天美| 亚洲av成人精品一区久久| 一级爰片在线观看| 成人一区二区视频在线观看| 免费看不卡的av| 欧美人与善性xxx| 特级一级黄色大片| 欧美日韩精品成人综合77777| 免费大片18禁| 晚上一个人看的免费电影| 内射极品少妇av片p| 91aial.com中文字幕在线观看| 亚洲精品国产成人久久av| 久久久久精品性色| 国产精品国产av在线观看| 两个人的视频大全免费| 国产精品三级大全| 日本一二三区视频观看| 色视频在线一区二区三区| 日本三级黄在线观看| 国产亚洲一区二区精品| 日韩欧美精品免费久久| 日韩强制内射视频| 黑人高潮一二区| 深夜a级毛片| 日产精品乱码卡一卡2卡三| 国产成人精品久久久久久| 亚洲成色77777| 国产精品99久久久久久久久| 真实男女啪啪啪动态图| 超碰97精品在线观看| 蜜桃久久精品国产亚洲av| 国产老妇伦熟女老妇高清| 成人无遮挡网站| 国产亚洲最大av| 99久久中文字幕三级久久日本| 只有这里有精品99| 99视频精品全部免费 在线| 欧美激情国产日韩精品一区| 不卡视频在线观看欧美| 国产精品.久久久| 肉色欧美久久久久久久蜜桃 | 永久免费av网站大全| 在线天堂最新版资源| 黄色日韩在线| 制服丝袜香蕉在线| freevideosex欧美| 日本午夜av视频| 九九久久精品国产亚洲av麻豆| 亚洲av福利一区| 久久久久久久精品精品| 91久久精品国产一区二区成人| 日韩欧美精品免费久久| 一级黄片播放器| 熟妇人妻不卡中文字幕| 秋霞伦理黄片| 久久女婷五月综合色啪小说 | a级毛片免费高清观看在线播放| 我要看日韩黄色一级片| 我要看日韩黄色一级片| 国产精品爽爽va在线观看网站| 亚洲综合色惰| av播播在线观看一区| 国产爱豆传媒在线观看| 亚洲欧美精品自产自拍| 老司机影院毛片| 成人高潮视频无遮挡免费网站| 赤兔流量卡办理| 亚洲无线观看免费| 国产精品国产av在线观看| 国产免费又黄又爽又色| av线在线观看网站| 蜜桃亚洲精品一区二区三区| 麻豆乱淫一区二区| 国产69精品久久久久777片| 免费看av在线观看网站| 男女无遮挡免费网站观看| 伦精品一区二区三区| 五月开心婷婷网| 看十八女毛片水多多多| 亚洲av一区综合| 国产免费一区二区三区四区乱码| 国产成人aa在线观看| 99热这里只有精品一区| 国产免费又黄又爽又色| 日韩,欧美,国产一区二区三区| av网站免费在线观看视频| 搡女人真爽免费视频火全软件| 一级毛片电影观看| 在现免费观看毛片| 国产乱人偷精品视频| 一个人看视频在线观看www免费| 天天躁夜夜躁狠狠久久av| av免费在线看不卡| 一个人观看的视频www高清免费观看| 男人爽女人下面视频在线观看| 少妇 在线观看| 亚洲欧洲日产国产| 亚洲av.av天堂| xxx大片免费视频| 亚洲av福利一区| 天美传媒精品一区二区| 欧美日韩亚洲高清精品| 国产精品人妻久久久久久| 一级av片app| 综合色丁香网| 亚洲四区av| 又爽又黄无遮挡网站| 老师上课跳d突然被开到最大视频| 夜夜爽夜夜爽视频| 日韩国内少妇激情av| 久久精品熟女亚洲av麻豆精品| 丝瓜视频免费看黄片| 国产在线男女| 久热久热在线精品观看| 亚洲精品久久午夜乱码| 自拍偷自拍亚洲精品老妇| 日韩人妻高清精品专区| 网址你懂的国产日韩在线| 久热这里只有精品99| 欧美成人午夜免费资源| 女的被弄到高潮叫床怎么办| 麻豆乱淫一区二区| 真实男女啪啪啪动态图| kizo精华| 久久精品人妻少妇| 啦啦啦在线观看免费高清www| 一级毛片我不卡| 亚洲av不卡在线观看| 小蜜桃在线观看免费完整版高清| 国产精品一区二区性色av| 熟妇人妻不卡中文字幕| 久久99热6这里只有精品| 亚洲天堂国产精品一区在线| 欧美高清成人免费视频www| 伊人久久精品亚洲午夜| 在线看a的网站| 国产伦精品一区二区三区视频9| 中国国产av一级| 51国产日韩欧美| 亚洲天堂国产精品一区在线| 亚洲成色77777| 日本午夜av视频| 久久久久久国产a免费观看| 亚洲精品日本国产第一区| 欧美成人精品欧美一级黄| 91狼人影院| 亚洲精品自拍成人| 国产男女内射视频| 能在线免费看毛片的网站| 精品亚洲乱码少妇综合久久| 人人妻人人澡人人爽人人夜夜| 久久99蜜桃精品久久| 丝袜美腿在线中文| 欧美另类一区| 日韩av免费高清视频| 深爱激情五月婷婷| 建设人人有责人人尽责人人享有的 | 国产成人精品福利久久| 少妇人妻一区二区三区视频| 国产综合精华液| 亚洲精品中文字幕在线视频 | 天美传媒精品一区二区| 成人二区视频| 日韩一区二区视频免费看| 亚洲av二区三区四区| 深夜a级毛片| 少妇裸体淫交视频免费看高清| 亚洲人与动物交配视频| 日韩欧美一区视频在线观看 | 久久精品国产a三级三级三级| 精品酒店卫生间| 欧美bdsm另类| 91久久精品国产一区二区成人| 晚上一个人看的免费电影| 天天一区二区日本电影三级| 国产真实伦视频高清在线观看| 亚洲av电影在线观看一区二区三区 | 大又大粗又爽又黄少妇毛片口| 高清午夜精品一区二区三区| 综合色av麻豆| 国产色婷婷99| 一区二区av电影网| 成人二区视频| 中文天堂在线官网| 校园人妻丝袜中文字幕| 男女无遮挡免费网站观看| av免费观看日本| 2021少妇久久久久久久久久久| 国产黄片美女视频| 免费大片18禁| 午夜视频国产福利| 九九久久精品国产亚洲av麻豆| 少妇 在线观看| 色网站视频免费| 激情 狠狠 欧美| 又爽又黄无遮挡网站| 免费观看av网站的网址| 免费观看无遮挡的男女| 亚洲成人一二三区av| 国产一区有黄有色的免费视频| 最近最新中文字幕大全电影3| 久久精品久久久久久噜噜老黄| 国产精品国产三级国产av玫瑰| 久久久久久久精品精品| 香蕉精品网在线| 色综合色国产| 在线观看一区二区三区| 免费观看无遮挡的男女| 国产久久久一区二区三区| 一级毛片黄色毛片免费观看视频| 99久久中文字幕三级久久日本| 国产亚洲av嫩草精品影院| 大又大粗又爽又黄少妇毛片口| 久久久亚洲精品成人影院| 激情五月婷婷亚洲| 人妻系列 视频| 久久鲁丝午夜福利片| 欧美日韩在线观看h| av在线播放精品| 久久久久久久久久久丰满| 国内精品美女久久久久久| 国产一区二区亚洲精品在线观看| 蜜桃亚洲精品一区二区三区| 亚洲精品乱码久久久久久按摩| 美女内射精品一级片tv| 日韩,欧美,国产一区二区三区| 色哟哟·www| 亚洲美女搞黄在线观看| 亚洲成人中文字幕在线播放| 99热全是精品| 国产高潮美女av| 亚洲精品影视一区二区三区av| 18+在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 欧美成人一区二区免费高清观看| 国产av不卡久久| 中国三级夫妇交换| 国语对白做爰xxxⅹ性视频网站| 精品一区二区三区视频在线| 嘟嘟电影网在线观看| 丝袜喷水一区| 91狼人影院| 国产老妇女一区| 国产乱人视频| 大片电影免费在线观看免费| 色哟哟·www| 国精品久久久久久国模美| 亚洲婷婷狠狠爱综合网| 国产老妇伦熟女老妇高清| 青春草国产在线视频| 精品亚洲乱码少妇综合久久| 亚洲色图av天堂| 亚洲国产av新网站| 五月天丁香电影| 亚洲av中文字字幕乱码综合| 亚洲自拍偷在线| 美女内射精品一级片tv| 国产成人精品一,二区| 国内少妇人妻偷人精品xxx网站| 在线观看国产h片| 欧美成人精品欧美一级黄| 久久久久久久久久人人人人人人| 国产亚洲91精品色在线| 秋霞在线观看毛片| 国产欧美亚洲国产| 日韩三级伦理在线观看| 性色av一级| 大陆偷拍与自拍| 九色成人免费人妻av| 建设人人有责人人尽责人人享有的 | 亚洲av在线观看美女高潮| 国产精品无大码| 97精品久久久久久久久久精品| 国产精品久久久久久久电影| 国模一区二区三区四区视频| 免费av不卡在线播放| av在线播放精品| 国产伦理片在线播放av一区| 一级二级三级毛片免费看| 青春草视频在线免费观看| 高清日韩中文字幕在线| 成年av动漫网址| 国产一区二区三区av在线| 亚洲精品第二区| 成人毛片a级毛片在线播放| 最近手机中文字幕大全| 国产精品久久久久久久电影| 午夜福利视频1000在线观看| 天堂中文最新版在线下载 | 日韩国内少妇激情av| 国产精品爽爽va在线观看网站| 日韩电影二区| 精品一区二区免费观看| 人妻少妇偷人精品九色| 国产又色又爽无遮挡免| 22中文网久久字幕| 少妇人妻一区二区三区视频| 精品久久久久久久久亚洲| 美女国产视频在线观看| 爱豆传媒免费全集在线观看| 久久人人爽人人爽人人片va| 男男h啪啪无遮挡| 国产美女午夜福利| 新久久久久国产一级毛片| 国产淫语在线视频| 久久久欧美国产精品| 亚洲av二区三区四区| 精品国产三级普通话版| 插阴视频在线观看视频| 看黄色毛片网站| 国产精品秋霞免费鲁丝片| 久久久久精品性色| 亚洲国产日韩一区二区| 欧美精品国产亚洲| 男人狂女人下面高潮的视频| 精品国产露脸久久av麻豆| 听说在线观看完整版免费高清| 尤物成人国产欧美一区二区三区| 国产黄色免费在线视频| 插阴视频在线观看视频| 亚洲性久久影院| 免费电影在线观看免费观看| 国产精品国产三级国产专区5o| 网址你懂的国产日韩在线| 国国产精品蜜臀av免费| av专区在线播放| 亚洲内射少妇av| 精品久久久噜噜| 国产免费又黄又爽又色| 97超视频在线观看视频| 国产色婷婷99| 久久久久久国产a免费观看| 亚洲三级黄色毛片| 日韩av免费高清视频| 成人综合一区亚洲| 国国产精品蜜臀av免费| 亚洲精品一二三| 精品一区在线观看国产| 久久久精品免费免费高清| 国国产精品蜜臀av免费| 久久亚洲国产成人精品v| 男女国产视频网站| 日本一本二区三区精品| 深夜a级毛片| 久久久久精品性色| 热99国产精品久久久久久7| 建设人人有责人人尽责人人享有的 | 国产黄色视频一区二区在线观看| av线在线观看网站| 国产永久视频网站| 久久精品久久精品一区二区三区| 日日啪夜夜爽| 黄片无遮挡物在线观看| 免费观看的影片在线观看| 久久久久久久久久成人| 大片电影免费在线观看免费| 高清毛片免费看| 午夜福利在线在线| 国产亚洲精品久久久com| 欧美老熟妇乱子伦牲交| 亚洲av免费在线观看| 欧美成人精品欧美一级黄| 亚洲精品自拍成人| 禁无遮挡网站| 黄色一级大片看看| 一级a做视频免费观看| 亚洲国产成人一精品久久久| 亚洲av中文av极速乱| 亚洲av免费高清在线观看| 国产黄片美女视频| 国精品久久久久久国模美| 国产成人一区二区在线| 国产黄a三级三级三级人| 亚洲天堂国产精品一区在线| 亚洲av福利一区| 欧美潮喷喷水| 嫩草影院精品99| 99久久精品一区二区三区| 夫妻午夜视频| 搞女人的毛片| 一区二区三区免费毛片| 亚洲第一区二区三区不卡| 国产老妇伦熟女老妇高清| 在线观看av片永久免费下载| 久久精品久久久久久久性| 日本午夜av视频| 高清毛片免费看| 七月丁香在线播放| 亚洲三级黄色毛片| av在线天堂中文字幕| 日韩国内少妇激情av| 卡戴珊不雅视频在线播放| 亚洲成人av在线免费| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠久久av| 麻豆精品久久久久久蜜桃| 亚洲成人中文字幕在线播放| 97在线人人人人妻| 中文精品一卡2卡3卡4更新| 日韩一区二区视频免费看| 我要看日韩黄色一级片| 五月开心婷婷网| 一区二区三区乱码不卡18| 亚洲欧美一区二区三区黑人 | videossex国产| 成年女人看的毛片在线观看| 少妇人妻精品综合一区二区| 亚洲成人中文字幕在线播放| 国产日韩欧美亚洲二区| 少妇丰满av| 晚上一个人看的免费电影| 亚洲av中文av极速乱| 久久久久国产精品人妻一区二区| 亚洲天堂国产精品一区在线| 欧美97在线视频| 久久久久久久久大av| 蜜桃久久精品国产亚洲av| 精品国产三级普通话版| 干丝袜人妻中文字幕| 免费高清在线观看视频在线观看| 最近的中文字幕免费完整| 人妻夜夜爽99麻豆av| 国产老妇伦熟女老妇高清| 麻豆久久精品国产亚洲av| 最近手机中文字幕大全| 人体艺术视频欧美日本| 性色avwww在线观看| 久久女婷五月综合色啪小说 | 少妇熟女欧美另类| 日韩一本色道免费dvd| 丰满少妇做爰视频| 精品国产露脸久久av麻豆| 国产亚洲精品久久久com| 高清午夜精品一区二区三区| 国产亚洲精品久久久com| 99精国产麻豆久久婷婷| 亚洲综合精品二区| 最近中文字幕2019免费版| 一二三四中文在线观看免费高清| 久久精品国产亚洲av涩爱| 一级毛片久久久久久久久女| 亚洲精品一二三| 我要看日韩黄色一级片| 免费看不卡的av| 午夜激情久久久久久久| 精品亚洲乱码少妇综合久久| 26uuu在线亚洲综合色| 特大巨黑吊av在线直播| 久久99蜜桃精品久久| 久久久久性生活片| 久久热精品热| 精品一区二区三区视频在线| 日本黄大片高清| 国产亚洲精品久久久com| 国产av码专区亚洲av| 六月丁香七月| 人人妻人人爽人人添夜夜欢视频 | 男人添女人高潮全过程视频| 亚洲av.av天堂| 极品少妇高潮喷水抽搐| 国产精品av视频在线免费观看| 欧美成人午夜免费资源| 搡女人真爽免费视频火全软件| 国产成人精品一,二区| 搡老乐熟女国产| 成人亚洲精品一区在线观看 | 蜜桃亚洲精品一区二区三区| 97超视频在线观看视频| 国产伦理片在线播放av一区| 丰满少妇做爰视频| 亚洲av国产av综合av卡| 蜜桃久久精品国产亚洲av| 国产av不卡久久| 内射极品少妇av片p| 一本久久精品| 国产黄色免费在线视频| 伊人久久国产一区二区| 日韩强制内射视频| .国产精品久久| 下体分泌物呈黄色| 日韩伦理黄色片| 大又大粗又爽又黄少妇毛片口| 国产综合精华液| 高清午夜精品一区二区三区| 久久99热6这里只有精品| 亚洲伊人久久精品综合| 18+在线观看网站| 亚洲精品影视一区二区三区av| av网站免费在线观看视频| 亚洲av不卡在线观看| 亚洲av福利一区| 美女高潮的动态| 久久久久久久久大av| 国产精品国产av在线观看| 卡戴珊不雅视频在线播放| 日韩欧美 国产精品| 女的被弄到高潮叫床怎么办| 久久久国产一区二区| 天天一区二区日本电影三级| 久久久久网色| 日本色播在线视频| 青春草视频在线免费观看| 啦啦啦中文免费视频观看日本| 日本黄大片高清| 91精品国产九色| 亚州av有码| 久久精品熟女亚洲av麻豆精品| 男女国产视频网站| 91精品一卡2卡3卡4卡| 欧美性感艳星| 国产色婷婷99| 国产一区二区三区综合在线观看 | 欧美日韩一区二区视频在线观看视频在线 | 王馨瑶露胸无遮挡在线观看| 高清视频免费观看一区二区| 99九九线精品视频在线观看视频| 国产成人一区二区在线| 精品少妇久久久久久888优播| 免费看光身美女| 一区二区三区四区激情视频| 亚洲经典国产精华液单| 最近手机中文字幕大全| 人妻夜夜爽99麻豆av| 国产美女午夜福利| 国产精品国产三级专区第一集| 国产69精品久久久久777片| 99热国产这里只有精品6| 日日啪夜夜爽| 自拍欧美九色日韩亚洲蝌蚪91 | videossex国产| 日日啪夜夜撸| 男女那种视频在线观看| 免费观看性生交大片5| 交换朋友夫妻互换小说| 九九爱精品视频在线观看| 简卡轻食公司| 成年女人看的毛片在线观看| 嘟嘟电影网在线观看| 国产av不卡久久| 亚洲色图av天堂| 一级黄片播放器| 丰满人妻一区二区三区视频av| av在线亚洲专区| 大香蕉97超碰在线| 精品国产乱码久久久久久小说| 午夜免费观看性视频| 国精品久久久久久国模美| 中文字幕av成人在线电影| 高清在线视频一区二区三区| 一级毛片久久久久久久久女| 1000部很黄的大片| 日本三级黄在线观看| 亚洲av电影在线观看一区二区三区 | 成人漫画全彩无遮挡| 国产午夜福利久久久久久| 在线亚洲精品国产二区图片欧美 | 女人久久www免费人成看片| 乱系列少妇在线播放| 国产成人91sexporn| 国产伦理片在线播放av一区| 老司机影院成人| 26uuu在线亚洲综合色| 18+在线观看网站| 2021少妇久久久久久久久久久| 日韩中字成人| a级一级毛片免费在线观看| 18禁裸乳无遮挡免费网站照片| 日韩欧美精品v在线| 美女被艹到高潮喷水动态| 天天躁日日操中文字幕| 自拍偷自拍亚洲精品老妇| 亚洲三级黄色毛片| av国产久精品久网站免费入址| 亚洲欧美精品专区久久| 成人一区二区视频在线观看| 18禁在线无遮挡免费观看视频| 国产亚洲5aaaaa淫片| 青春草视频在线免费观看| 好男人视频免费观看在线| 一区二区三区精品91| 精品视频人人做人人爽| 97热精品久久久久久| 国产精品无大码| 人妻 亚洲 视频| 高清欧美精品videossex| 免费高清在线观看视频在线观看| 欧美一级a爱片免费观看看| 69av精品久久久久久| 婷婷色麻豆天堂久久| 观看免费一级毛片| 丝袜美腿在线中文| 国产黄频视频在线观看| 精品一区二区三卡|