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

    一種澆注PBX固化過程的實驗與數(shù)值模擬

    2015-05-10 02:40:30張百磊常雙君袁俊明歐亞鵬
    含能材料 2015年3期
    關(guān)鍵詞:藥柱應力場光柵

    張百磊,常雙君,袁俊明,歐亞鵬

    (中北大學化工與環(huán)境學院,山西 太原 030051)

    隨著炸藥新材料的不斷應用及炸藥配方與工藝的優(yōu)化和發(fā)展,澆注高聚物粘結(jié)炸藥(PBX)由于其高能量、低感度和高安全性等多種優(yōu)點,正在被越來越廣泛地應用在高技術(shù)武器系統(tǒng)中[1-2]。

    在澆注和固化過程中由于工藝條件不當會導致成型藥柱內(nèi)部熱應力和氣孔的出現(xiàn),因此需要對炸藥傳統(tǒng)固化工藝進行優(yōu)化,但此方面的研究主要靠反復實驗的方法,周期較長且缺乏科學理論指導。近年來凝固過程的數(shù)值模擬由金屬材料推廣至含能材料領(lǐng)域,國內(nèi)專家學者已經(jīng)對熔鑄類炸藥凝固過程的數(shù)值模擬進行了大量的探索和研究,李敬明等[3]用ProCAST軟件對熔黑梯(RHT)炸藥凝固過程建立了數(shù)學模型并進行了計算,結(jié)果表明,該軟件適合模擬炸藥的凝固過程。馬松等[4]通過水浴和自然條件下熔鑄炸藥凝固過程的數(shù)值計算,得出水浴條件下熔鑄炸藥質(zhì)量更好。黃勇[5]通過對熔鑄炸藥澆注過程進行數(shù)值模擬,發(fā)現(xiàn)澆注系統(tǒng)與環(huán)境間的換熱對凝固過程有著決定性的作用,特別是冒口部分,如果換熱系數(shù)過大,根本起不到保溫補縮的作用; 另外加壓能顯著減少縮松缺陷。馬松等[6]對2,4-二硝基苯甲醚(DNAN)炸藥凝固過程進行了實驗與數(shù)值模擬,結(jié)果表明,DNAN在96 ℃附近出現(xiàn)溫度滯后現(xiàn)象,中心部位最明顯,凝固后最大孔隙率可達0.38。然而關(guān)于澆注PBX固化工藝過程的數(shù)值模擬在國內(nèi)外報道較少,而澆注PBX的固化過程和熔鑄炸藥的凝固過程又有著本質(zhì)的區(qū)別:前者主要是靠發(fā)生化學交聯(lián)反應放出熱量完成,后者是一個物理降溫過程?;诖嗽跀?shù)值模擬過程中考慮化學反應的放熱量,完成對澆注PBX固化技術(shù)中實現(xiàn)計算機模擬與仿真,將對提高該炸藥的質(zhì)量和確保工藝安全性起到重要作用。

    本研究對澆注PBX固化過程中溫度場和應力場進行了數(shù)值模擬,得出了溫度場和應力場的變化規(guī)律,并與實驗測試結(jié)果對比,以期為澆注PBX固化工藝參數(shù)的優(yōu)化、方案的改進等提供參考。

    2 實驗部分

    2.1 實驗裝置

    FBG-3000型干涉解調(diào)儀;多點分布式Bragg光柵組,3根光柵之間距離24 mm,每根光柵上的測試點間距60 mm,測試點3、6、9距離模具內(nèi)底部60 mm,測試點1、4、7距離上表面60mm,測試點2、5、8距離上下測試點均為60 mm。實驗前對光柵進行標定;模具尺寸為Φ110 mm×245 mm,壁厚5 mm;AHX-8型水浴烘箱。系統(tǒng)原理示意圖如圖1所示。

    2.2 實驗原理

    澆注PBX的固化原理一般認為[7]是澆注PBX粘結(jié)劑體系的反應,即羥基與異氰酸酯基反應生成氨基甲酸酯。反應機理如式(a)、(b)、(c)。

    圖1 Bragg光柵測量炸藥溫度原理圖
    Fig.1 Principal diagram of measuring the temperature of explosive by Bragg grating

    測量溫度前需要將Bragg光柵固定在模具中,其測量溫度的基本原理[8]是:通過外界溫度微小的變化引起光柵折射率的改變,從而引起光柵的Bragg波長產(chǎn)生微小的移位,通過測量Bragg波長的變化,利用其反射波長與溫度之間的線性關(guān)系,經(jīng)過標定,實時顯示溫度。實驗采用FBG-3000型干涉解調(diào)儀,實現(xiàn)Bragg光柵波長的測量。

    2.3 實驗過程

    本研究采用的澆注PBX配方為64%RDX、20%Al和16%粘結(jié)劑體系。澆注PBX采用捏合-澆注-加熱固化技術(shù)制備,將3根光柵按圖1所示位置固定在模具中,在加熱保溫條件下將藥漿澆注到模具內(nèi),放入水浴烘箱在60 ℃加熱固化,開啟FBG-3000型干涉解調(diào)儀實時記錄藥柱溫度的變化,溫度采集時間間隔設(shè)置為600 s,在藥柱溫度降到60 ℃穩(wěn)定后,關(guān)閉烘箱,當藥柱冷卻至室溫(約25 ℃)后停止實驗。

    3 數(shù)值模擬

    3.1 模型建立

    由圖1可得模擬對象是圓柱型固化系統(tǒng),其具有旋轉(zhuǎn)對稱的特點,故在前處理軟件PRO/E中建立1/4模型,再導入ProCAST進行網(wǎng)格劃分,總節(jié)點數(shù)8058,總單元數(shù)36800,如圖2所示。

    圖2 藥柱固化模型
    Fig.2 Grain curing model

    3.2 溫度場求解設(shè)置

    采用有限單元法[9-10]對澆注PBX的固化過程進行了數(shù)值模擬,其中溫度場分布采用不穩(wěn)定熱傳導的偏微分方程進行計算[11],其數(shù)值計算依據(jù)瞬態(tài)傳熱偏微分方程,如(1)式所示。由澆注PBX的固化原理可知,該模擬溫度場模擬的邊界條件取第三類邊界條件,即外界環(huán)境溫度Tf和換熱系數(shù)α已知,如(2)式所示。表1給出了數(shù)值模擬所采用的材料性能參數(shù)、初始條件及邊界條件參數(shù)[12-13]。

    (1)

    (2)

    式中,方程(1)左邊表示熱積蓄項; 右邊第一、二、三項表示導熱項; 第四項為材料內(nèi)部熱源強度; 其中k為炸藥導熱系數(shù),W·m-1·K-1;cp為比熱,kJ·kg-1·K-1;ρ為密度,kg·m-3;T為溫度,℃;t為時間,s;Q為炸藥內(nèi)部熱源強度,kJ·m-3·s-1; 式(2)中,Tf為環(huán)境溫度,℃。

    表1 澆注PBX及模具的性能參數(shù)
    Table 1 Property parameters of cast PBX and steel mould

    itemsk/W·m-1·K-1cp/kJ·kg-1·K-1ρ/kg·m-3explosive0.301.301656mould43.530.487850

    澆注之前要對模具進行預熱,所以模具和藥漿的初始溫度都設(shè)為60 ℃。澆注PBX的固化過程會發(fā)生交聯(lián)反應而放出熱量,該熱量即為式(1)中的Q。該反應熱的計算是多英全等[14]以Van Krevelen和 Chermin等人的基團估算方法及數(shù)據(jù)為基礎(chǔ),結(jié)合高分子內(nèi)聚能的概念,由(3)式得出:

    (3)

    參照聚氨酯類彈性體(TPU)結(jié)構(gòu)單元的形式及Van Krevelen和Chermin的數(shù)據(jù),可得到1mol TPU 結(jié)構(gòu)單元為假想氣態(tài)時的生成焓與硬段含量H(即固化劑的含量)之間具有表達式(4):

    (4)

    依據(jù)TPU結(jié)構(gòu)可將其分為3個部分:其中,A代表硬段的結(jié)構(gòu)單元,該結(jié)構(gòu)單元是由異氰酸酯基和羥基之間發(fā)生反應生成的絡(luò)合物; A與B構(gòu)成硬段,B代表固化劑,即該硬段是由上述絡(luò)合物和固化劑發(fā)生反應生成的另一種絡(luò)合物; C為軟段部分,代表粘結(jié)劑的結(jié)構(gòu)單元。

    同理,可得1 mol TPU結(jié)構(gòu)單元中硬段與軟段內(nèi)能總和與硬段含量H之間可表示如下:

    (5)

    (6)

    3.3 應力場求解設(shè)置

    熱應力數(shù)值模擬是通過對藥柱固化過程中溫度場與熱應力場的耦合計算來得到應力場模擬結(jié)果。應力場計算模型包括線彈性模型,彈塑性模型及粘彈塑性模型等[15],考慮到澆注PBX的力學性能,模擬采用彈塑性雙線性硬化模型對其成型過程中的應力場進行計算,圖3給出了該模型的應力應變曲線示意圖[15]。隨著應力的增大,模型由彈性模量計算轉(zhuǎn)為了硬度計算,其轉(zhuǎn)折點是屈服應力。應變也隨著溫度T的升高而增大。由于炸藥澆注成型時采用金屬模具,其模量和強度均遠大于炸藥的模量和強度,且澆注溫度(60 ℃)比金屬的熔點低,因此在進行應力場計算時將金屬簡化為剛性模型。表2給出了應力場計算所需參數(shù)[9,16-17]。

    圖3 彈塑性雙線性硬化模型[15]
    Fig.3 Model of elastic-plastic bilinear hardening[15]

    表2 應力場計算所需性能參數(shù)
    Table2 Property parameters needed in calculation of the stress fields

    itemsα/W·m-2·K-1E/MPaσ/MPavexplosive121087.60(T=23℃)7.18(T=40℃)6.10(T=60℃)12(T=25℃)9(T=40℃)3(T=60℃)0.50mould---

    Note:Eis elastic modulus,σis yield stress,vis Poisson′s ratio.αis coefficient of heat transfer.

    模具選用剛性模具,所以相對炸藥來說模具的E、σ、v、都沒有數(shù)值,但是其換熱系數(shù)α是需要設(shè)置的。模擬中的換熱系數(shù)α表征的是藥柱與模具之間傳熱的快慢,實際生產(chǎn)中模具的材料、厚度及有無絕熱層等因素都會影響到換熱系數(shù)。實驗中模具的換熱系數(shù)為10 W·m-2·K-1,故本模擬換熱系數(shù)分別取8,10,12 W·m-2·K-1進行計算。

    4 計算結(jié)果與討論

    4.1 溫度場數(shù)值模擬結(jié)果及實驗驗證

    依據(jù)1/4模型的建立和參數(shù)的設(shè)置,在3.88×105s藥柱進入降溫階段,藥柱內(nèi)部的溫度梯度最大,該時間點藥柱內(nèi)部的溫度分布,如圖4所示。由于藥柱的熱導率數(shù)值(熱導率表征熱量在一種材料中的導熱速率)很小,所以溫度在藥柱內(nèi)部的分布由中心向外跨度較大。另外藥柱冷卻過程放熱相對較少,所以該過程也可被認為是由中心向外部傳熱降溫的過程。圖5、6、7分別給出了第1、2、3、4、5、6、7、8、9點溫度變化曲線的數(shù)值模擬結(jié)果和實驗測試結(jié)果,即每根光柵所在藥柱位置縱軸方向上的溫度分布。

    圖4 3.88×105s時藥柱的溫度場分布
    Fig.4 Temperature field distribution in grain at 3.88×105s

    從圖5~圖7可見,澆注PBX中不同位置溫度變化趨勢的數(shù)值模擬結(jié)果和實驗測試結(jié)果基本一致。由于澆注PBX在固化過程中會發(fā)生反應放出熱量,故所有的測試點都會出現(xiàn)一段溫度維持在60 ℃以上的時間平臺,最高可達65 ℃;大約96 h后固化反應結(jié)束,體系熱量失衡開始降溫。當藥柱溫度降到60 ℃穩(wěn)定后,關(guān)閉烘箱,整個體系自然降溫。其中7、8、9三個測試點因離模具內(nèi)壁較近,向外部環(huán)境傳熱相對較快,所以這三個點所在光柵上測得的溫度會比其它兩個光柵上測得的溫度低。同理4、5、6三個點所在光柵上測得的溫度會比1、2、3點所在光柵上測得的溫度低。

    另外由圖5~圖7可分別觀察到,藥柱內(nèi)部縱向溫差較大,可能是由于測試條件及方法的限制,所以在后續(xù)的研究中可以通過在模具內(nèi)壁增加一層絕熱層或者改變模具材料及厚度來降低溫差。

    a.tested curves

    b.simulated curves

    圖5 第1、2、3點溫度變化曲線
    Fig.5 Temperature change curves of points 1,2,3

    a.tested curves

    b.simulated curves

    圖6 第4、5、6點溫度變化曲線
    Fig.6 Temperature change curves of points 4,5,6

    a.tested curves

    b.simulated curves

    圖7 第7、8、9點溫度變化曲線
    Fig.7 Temperature change curves of points 7,8,9

    4.2 應力場數(shù)值模擬結(jié)果

    依據(jù)設(shè)置的3種換熱系數(shù)8,10,12 W·m-2·K-1,計算得到藥柱在3.88×105s有效應力分布如圖8所示。應力集中區(qū)域主要分布在藥柱下半部分,隨著換熱系數(shù)的增加,藥柱內(nèi)部應力集中的區(qū)域及最大有效應力也隨之增大。這是因為隨著換熱系數(shù)的增大,藥柱通過模具向外界環(huán)境傳遞出熱量的速度就會相應加快,導致藥柱內(nèi)部溫度梯度過大,從而影響到藥柱內(nèi)部的應力分布。3種換熱系數(shù)下的最大有效應力見表3。

    表3 澆注PBX在3種換熱系數(shù)下的最大有效應力
    Table 3 Maximum effective stress of cast PBX under three kinds of heat transfer coefficient

    α/W·m-2·K-181012maximumeffectivestress/kPa1.955.666.69

    a.8 W·m-2·K-1b.10 W·m-2·K-1c.12 W·m-2·K-1
    圖8 澆注PBX在不同換熱系數(shù)下的有效應力分布
    Fig.8 Effective stress distribution of cast PBX under different heat transfer coefficient

    實驗采用的澆注類PBX的屈服應力[18]是0.06 GPa,所以藥柱內(nèi)部不會因為應力集中的原因而出現(xiàn)熱裂痕。但是藥柱內(nèi)部應力場會隨著炸藥的配方、環(huán)境溫度的變化率、模具材料及厚度等因素的變化而變化,例如在藥柱冷卻階段,如果環(huán)境溫度變化過快,就會導致藥柱內(nèi)部溫度梯度過大,產(chǎn)生較大的熱應力,這將會對藥柱內(nèi)部質(zhì)量產(chǎn)生很大影響。

    5 結(jié) 論

    (1) 結(jié)合澆注PBX性能參數(shù),利用有限單元法對其固化過程中的溫度場和應力場進行了計算,并與實驗相結(jié)合。環(huán)境溫度是60 ℃時,固化初期藥柱溫度有一定程度的升高,其中藥柱中心溫度最高為65 ℃,并由藥柱中心向模具內(nèi)壁逐漸降低;在3.88×105s時藥柱內(nèi)部溫度梯度較大。

    (2) 應力集中區(qū)域主要分布在藥柱下半部分;在3種換熱系數(shù)下,藥柱內(nèi)部在3.88×105s時呈現(xiàn)出不同的應力集中,其最大有效應力會隨著換熱系數(shù)的升高而增大。當換熱系數(shù)為12 W·m-2·K-1時,最大有效應力是6.69 kPa。

    參考文獻:

    [1] 孫業(yè)斌,惠君明,曹欣茂,等.軍用混合炸藥[M].兵器工業(yè)出版社,1995:349-358.

    SUN Ye-bin,HUI Jun-ming,CAO Xin-mao,et al.Military composite explosive[M].Weapon Industry Press,1995:349-358.

    [2] 李媛嬡,南海.國外澆注PBX炸藥在硬目標侵徹武器中的應用[J].飛航導彈,2012 (11):88-91.

    LI Yuan-yuan,NAN Hai.The application of casting PBX explosives in the hard target penetration weapon at foreign[J].AerodynamicMissileJournal,2012 (11):88-91.

    [3] 李敬明,田勇,張明,等.熔黑梯炸藥固化過程的數(shù)值模擬與實驗驗證[J].含能材料,2009,17(4):428-430.

    LI Jing-ming,TIAN Yong,ZHANG Ming,et al.Numerical simulation and experimental validation of RHT solidification process[J].ChineseJournalofEnergeticMaterials(HannengCailiao).2009,17(4):428-430.

    [4] 馬松,袁俊明,劉玉存,等.水浴和自然條件下熔鑄炸藥固化過程的數(shù)值計算[J].火炸藥學報,2013,36(3):54-58.

    MA Song,YUAN Jun-ming,LIU Yu-cun,et al.Numerical calculation of solidification process of melt-cast explosive in water bath and natural conditions[J].ChineseJournalofExplosives&Propellants,2013,36(3):54-58.

    [5] 黃勇.炸藥熔鑄過程數(shù)值模擬研究[J].中國工程物理研究院科技年報,2012:106-107.

    [6] 馬松,袁俊明,劉玉存,等.DNAN炸藥凝固過程的實驗與數(shù)值模擬[J].含能材料,2014,22(2):240-244.

    MA Song,YUAN Jun-ming,LIU Yu-cun,et al.Experiment and numerical simulation of DNAN solidification process[J].ChineseJournalofEnergeticMaterials(HannengCailiao).2014,22(2):240-244.

    [7] 陳春燕,王曉峰,高立龍,等.固化參數(shù)對HTPB固化反應和力學性能的影響[J].化學研究與應用,2013,25(10):1381-1385.

    CHEN Chun-yan,WANG Xiao-feng,GAO Li-long,et al.Effects of NCO/OH molar ratio on the cure reaction and mechanical property of HTPB[J].ChemicalResearchAndApplication,2013,25(10):1381-1385.

    [8] 姜志剛.FBG光纖光柵的原理和應用[J].中國水運,2008(5):128-129.

    JIANG Zhi-gang.Principle and typical current applications of fiber grating [J].ChinaWaterTransport,2008(5):128-129.

    [9] 田勇.炸藥熔鑄成型過程監(jiān)測評價及數(shù)值模擬研究[D].北京:中國科學院研究生院(工程熱物理研究所),2010.

    TIAN Yong.Process monitoring/evaluation and numerical simulation during casting explosive solidif-cation[D].Beijing:Graduate University of China Academy of Sciences(Institute of Engineering Thermophysics),2010.

    [10] Mudryy R,Singh S.Modeling and simulation of melt cast explosives[C]∥Insensitive Munitions and Energetic Materials Technology Symposium.Arlington 2007.

    [11] 曹堅,錢蘇翔,焦衛(wèi)東,等.基于有限單元法的鑄件三維溫度場和應力場耦合計算研究[J].機械制造,2007,45(1):15-17.

    CAO Jian,QIAN Su-xiang,JIAO Wei-dong,et al.Numerical simulation of finite elements in 3-dimension temperature field and mechanical effects[J].Machinery,2007,45(1):15-17.

    [12] Maienschein J L,Nichols III A L,Reaugh J E,et al.Integrated experimental and modeling studies to predict the impact response of explosives and propellants[R].Lawrence Livermore National Laboratory (LLNL),Livermore,CA,2005.

    [13] McClelland M A,Tran T D,Cunningham B J,et al.Cook off response of PBXN-109:material characterization and ALE3D model [C]∥Proceedings of JANNAF 19th Propulsion Systems Hazards Subcommittee Meeting.2000:191.

    [14] 多英全,陳福泰.熱塑性聚氨酯彈性體及推進劑生成焓的估算[J].推進技術(shù),2000,21(6):79-82.

    DUO Ying-quan,CHEN Fu-tai.Estimation of formation enthalpy of thermoplastic polyurethane binder for novel thermoplastic propellant[J].JournalofPropulsionTechnology,2000,21(6):79-82.

    [15] 李日,馬軍賢,崔啟玉.鑄造工藝仿真ProCAST從入門到精通[M].中國水利水電出版社,2010:265-266.

    LI Ri,MA Jun-xian,CUI Qi-yu.Casting process simulation ProCAST from abecedarian to proficiency[M].ChinaWaterPowerPress,2010:265-266.

    [16] Hamshere B L,Lochert I J,Dexter R M.Evaluation of PBXN-109:the explosive fill for the Penguin Anti-Ship Missile Warhead[R].Defence Science And Technology Organisation Salisbury(Australia)Systems Sciences Lab,2003.

    [18] McClelland M A,Tran T D,Cunningham B J,et al.Cookoff response of PBXN-109:material characterization and ALE3D thermal predictions[C]∥Proceedings of 50th JANNAF Propulsion Meeting,Salt Lake City,UT,CPIA.2001.

    猜你喜歡
    藥柱應力場光柵
    高聚物黏結(jié)炸藥沖擊波感度試驗方法
    四川化工(2022年6期)2023-01-15 10:54:54
    更 正
    含能材料(2017年6期)2017-03-07 06:32:07
    CDIO教學模式在超聲光柵實驗教學中的實踐
    基于LabView的光柵衍射虛擬實驗研究
    鋁合金多層多道窄間隙TIG焊接頭應力場研究
    焊接(2016年9期)2016-02-27 13:05:22
    減面燃燒規(guī)律的藥柱選用準則*
    考慮斷裂破碎帶的丹江口庫區(qū)地應力場與水壓應力場耦合反演及地震預測
    固體推進劑組合藥柱的界面力學性能
    火炸藥學報(2014年1期)2014-03-20 13:17:28
    基于位移相關(guān)法的重復壓裂裂縫尖端應力場研究
    斷塊油氣田(2014年5期)2014-03-11 15:33:49
    岸坡應力場及卸荷帶劃分量化指標研究
    熟妇人妻久久中文字幕3abv| 久久久久九九精品影院| 久热爱精品视频在线9| 国产精品免费视频内射| 成年版毛片免费区| 亚洲av中文字字幕乱码综合| 国内精品久久久久精免费| 久久人妻av系列| 欧洲精品卡2卡3卡4卡5卡区| 国产69精品久久久久777片 | 91av网站免费观看| 床上黄色一级片| 亚洲av电影在线进入| 免费看美女性在线毛片视频| 黄频高清免费视频| 中文字幕熟女人妻在线| 国产午夜精品久久久久久| 九九热线精品视视频播放| 亚洲人成77777在线视频| 国产精品精品国产色婷婷| 在线观看美女被高潮喷水网站 | 18禁黄网站禁片免费观看直播| 美女免费视频网站| 黄色女人牲交| 麻豆国产97在线/欧美 | 麻豆成人午夜福利视频| 亚洲成av人片免费观看| 亚洲成av人片免费观看| 国产熟女午夜一区二区三区| 最近在线观看免费完整版| 中文字幕精品亚洲无线码一区| 免费在线观看黄色视频的| 国产精品影院久久| 国产久久久一区二区三区| 成人18禁在线播放| 国内精品久久久久精免费| 最近最新中文字幕大全电影3| 91大片在线观看| 50天的宝宝边吃奶边哭怎么回事| 日本一本二区三区精品| 亚洲中文字幕日韩| 在线国产一区二区在线| netflix在线观看网站| 不卡av一区二区三区| 手机成人av网站| 国产精品1区2区在线观看.| 久久精品91蜜桃| 国产精品久久久久久亚洲av鲁大| 亚洲精品一区av在线观看| 在线永久观看黄色视频| 国产免费av片在线观看野外av| 国产精品久久久人人做人人爽| 日韩欧美一区二区三区在线观看| 久久这里只有精品中国| 久久天堂一区二区三区四区| 91麻豆精品激情在线观看国产| 亚洲精品粉嫩美女一区| 久久久久久久午夜电影| 国产欧美日韩精品亚洲av| 听说在线观看完整版免费高清| 69av精品久久久久久| 啪啪无遮挡十八禁网站| 亚洲精品久久国产高清桃花| 日韩高清综合在线| 给我免费播放毛片高清在线观看| 18禁黄网站禁片午夜丰满| 亚洲免费av在线视频| 欧美日韩亚洲国产一区二区在线观看| 国产精品 欧美亚洲| 午夜精品一区二区三区免费看| 黄色毛片三级朝国网站| netflix在线观看网站| 午夜两性在线视频| 国产精品久久久久久久电影 | 黄频高清免费视频| 最近最新免费中文字幕在线| 久9热在线精品视频| 成人国语在线视频| 久久久久精品国产欧美久久久| 波多野结衣巨乳人妻| 午夜久久久久精精品| 最新美女视频免费是黄的| 免费一级毛片在线播放高清视频| 99久久99久久久精品蜜桃| 国产亚洲精品一区二区www| 中文资源天堂在线| 亚洲自偷自拍图片 自拍| 又黄又爽又免费观看的视频| 国产亚洲精品久久久久久毛片| 国产99白浆流出| av中文乱码字幕在线| www日本黄色视频网| 久久精品aⅴ一区二区三区四区| 黄色女人牲交| 99在线人妻在线中文字幕| 少妇裸体淫交视频免费看高清 | 国产激情久久老熟女| 免费一级毛片在线播放高清视频| 亚洲一区中文字幕在线| 夜夜看夜夜爽夜夜摸| 此物有八面人人有两片| 午夜a级毛片| 变态另类丝袜制服| 日本a在线网址| 99国产精品一区二区蜜桃av| 这个男人来自地球电影免费观看| 久久九九热精品免费| 黄色片一级片一级黄色片| 欧美极品一区二区三区四区| 一本精品99久久精品77| 亚洲av片天天在线观看| 白带黄色成豆腐渣| 免费在线观看影片大全网站| 久久人人精品亚洲av| 在线永久观看黄色视频| or卡值多少钱| 中亚洲国语对白在线视频| 日本一二三区视频观看| 中文字幕精品亚洲无线码一区| 亚洲免费av在线视频| 中文在线观看免费www的网站 | 亚洲男人天堂网一区| 黄色视频,在线免费观看| 这个男人来自地球电影免费观看| 亚洲人成伊人成综合网2020| 美女扒开内裤让男人捅视频| 人妻丰满熟妇av一区二区三区| 久久久久免费精品人妻一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜成年电影在线免费观看| 怎么达到女性高潮| 午夜日韩欧美国产| 久久久久国产精品人妻aⅴ院| 亚洲电影在线观看av| 国产精品影院久久| 国产欧美日韩一区二区三| 视频区欧美日本亚洲| 一级片免费观看大全| 午夜精品久久久久久毛片777| 国产单亲对白刺激| 国产在线精品亚洲第一网站| 亚洲欧美日韩东京热| 亚洲熟女毛片儿| 久久久久久国产a免费观看| 俄罗斯特黄特色一大片| 欧美黑人巨大hd| 操出白浆在线播放| 人妻久久中文字幕网| 男男h啪啪无遮挡| 欧美成人午夜精品| 亚洲av成人不卡在线观看播放网| 国产真人三级小视频在线观看| 久久久久久大精品| 国产精品久久视频播放| 国产v大片淫在线免费观看| 可以在线观看毛片的网站| 国产精品免费视频内射| 成人三级黄色视频| 久久精品综合一区二区三区| 黄片大片在线免费观看| 狂野欧美激情性xxxx| 国产激情偷乱视频一区二区| 黄色成人免费大全| 午夜影院日韩av| 舔av片在线| 99久久国产精品久久久| 国产成人av激情在线播放| 欧美精品亚洲一区二区| 欧美又色又爽又黄视频| 亚洲激情在线av| 最近最新中文字幕大全免费视频| 99久久国产精品久久久| 国产区一区二久久| 成人av在线播放网站| 在线观看午夜福利视频| 亚洲av第一区精品v没综合| 国产精品九九99| 中文字幕人妻丝袜一区二区| 99re在线观看精品视频| 午夜福利高清视频| 两个人的视频大全免费| 黄色成人免费大全| 超碰成人久久| 亚洲中文日韩欧美视频| 夜夜躁狠狠躁天天躁| 日韩欧美在线二视频| 国产欧美日韩一区二区三| 午夜亚洲福利在线播放| 免费搜索国产男女视频| av超薄肉色丝袜交足视频| 亚洲精品久久国产高清桃花| 丝袜美腿诱惑在线| 在线看三级毛片| 50天的宝宝边吃奶边哭怎么回事| 熟妇人妻久久中文字幕3abv| 色哟哟哟哟哟哟| 免费av毛片视频| 一级片免费观看大全| 亚洲avbb在线观看| 国产亚洲精品综合一区在线观看 | 性色av乱码一区二区三区2| 午夜激情福利司机影院| 亚洲 欧美一区二区三区| 色哟哟哟哟哟哟| 久久天堂一区二区三区四区| 日韩av在线大香蕉| 久久久久久亚洲精品国产蜜桃av| 脱女人内裤的视频| 国产亚洲精品久久久久久毛片| 国产免费av片在线观看野外av| 日日干狠狠操夜夜爽| 国产亚洲欧美98| 国产激情偷乱视频一区二区| 国产视频一区二区在线看| 日韩 欧美 亚洲 中文字幕| av在线天堂中文字幕| 香蕉国产在线看| 黄色丝袜av网址大全| 色老头精品视频在线观看| 精华霜和精华液先用哪个| 嫩草影院精品99| 麻豆国产av国片精品| 日韩欧美国产在线观看| 国语自产精品视频在线第100页| 男女午夜视频在线观看| 久久久精品欧美日韩精品| 国产成人一区二区三区免费视频网站| 国产成人av激情在线播放| 小说图片视频综合网站| 精品久久久久久,| 俄罗斯特黄特色一大片| 免费在线观看亚洲国产| 免费一级毛片在线播放高清视频| 三级国产精品欧美在线观看 | 国产精品免费一区二区三区在线| 免费看美女性在线毛片视频| 亚洲人成伊人成综合网2020| 精品高清国产在线一区| 久久久久久九九精品二区国产 | 九色国产91popny在线| 黄色成人免费大全| 欧美黑人精品巨大| 最近最新中文字幕大全免费视频| 免费看日本二区| 国产成+人综合+亚洲专区| 亚洲自拍偷在线| 91av网站免费观看| xxxwww97欧美| 久久精品国产综合久久久| 人人妻,人人澡人人爽秒播| 日日摸夜夜添夜夜添小说| 国产免费男女视频| 日本五十路高清| 亚洲 欧美 日韩 在线 免费| 国产亚洲精品综合一区在线观看 | 亚洲五月婷婷丁香| 欧美一级a爱片免费观看看 | 久久精品国产亚洲av高清一级| 久久久久九九精品影院| 91老司机精品| 两个人的视频大全免费| 亚洲av成人av| 久久久久九九精品影院| 免费一级毛片在线播放高清视频| 婷婷精品国产亚洲av在线| 18禁美女被吸乳视频| 最近在线观看免费完整版| 精品一区二区三区av网在线观看| 在线观看66精品国产| 国产精品久久久久久精品电影| 久久九九热精品免费| 国产99久久九九免费精品| 黄色 视频免费看| 久久久久久国产a免费观看| 日本一区二区免费在线视频| 99国产精品99久久久久| www.自偷自拍.com| 欧美日韩乱码在线| 国产区一区二久久| 欧美绝顶高潮抽搐喷水| 美女大奶头视频| 国内毛片毛片毛片毛片毛片| 国产v大片淫在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 国产成人系列免费观看| 亚洲av第一区精品v没综合| 特大巨黑吊av在线直播| 欧美激情久久久久久爽电影| 国产精品香港三级国产av潘金莲| 又大又爽又粗| 久久国产精品人妻蜜桃| 久久久久久九九精品二区国产 | 可以在线观看毛片的网站| 人人妻,人人澡人人爽秒播| 亚洲九九香蕉| 两人在一起打扑克的视频| 欧美绝顶高潮抽搐喷水| 天天躁夜夜躁狠狠躁躁| 舔av片在线| 中文字幕精品亚洲无线码一区| 久久精品国产综合久久久| 免费看a级黄色片| 免费人成视频x8x8入口观看| 国产成人啪精品午夜网站| 久久国产精品影院| 天堂动漫精品| 久久这里只有精品19| 搡老熟女国产l中国老女人| 一级a爱片免费观看的视频| 色精品久久人妻99蜜桃| 中文字幕最新亚洲高清| 精品一区二区三区视频在线观看免费| 丰满人妻一区二区三区视频av | 色综合婷婷激情| 巨乳人妻的诱惑在线观看| 人成视频在线观看免费观看| 蜜桃久久精品国产亚洲av| 久久久久国产精品人妻aⅴ院| 在线永久观看黄色视频| 免费看十八禁软件| 国产探花在线观看一区二区| 女人高潮潮喷娇喘18禁视频| 亚洲七黄色美女视频| 国产av一区在线观看免费| 一区二区三区高清视频在线| 亚洲无线在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产乱人伦免费视频| 岛国视频午夜一区免费看| 国内精品久久久久久久电影| 亚洲色图av天堂| 岛国视频午夜一区免费看| 午夜激情福利司机影院| 亚洲va日本ⅴa欧美va伊人久久| 国产私拍福利视频在线观看| 国产麻豆成人av免费视频| 免费在线观看视频国产中文字幕亚洲| 一a级毛片在线观看| 在线观看免费日韩欧美大片| 午夜精品一区二区三区免费看| 日韩免费av在线播放| 亚洲人成77777在线视频| 欧美日韩亚洲国产一区二区在线观看| 日本成人三级电影网站| 亚洲熟妇中文字幕五十中出| 可以在线观看的亚洲视频| 欧美绝顶高潮抽搐喷水| 宅男免费午夜| 一级作爱视频免费观看| 日日夜夜操网爽| 亚洲美女黄片视频| 中国美女看黄片| 脱女人内裤的视频| 亚洲真实伦在线观看| 免费搜索国产男女视频| 怎么达到女性高潮| 男女午夜视频在线观看| 国产黄色小视频在线观看| 日韩欧美国产在线观看| 日韩欧美 国产精品| 国产亚洲精品久久久久5区| 夜夜夜夜夜久久久久| 一级毛片精品| 久久久国产欧美日韩av| 亚洲五月天丁香| 国产又黄又爽又无遮挡在线| 神马国产精品三级电影在线观看 | 国产精品乱码一区二三区的特点| 美女大奶头视频| 欧美成狂野欧美在线观看| 91九色精品人成在线观看| 国产精品久久久久久亚洲av鲁大| 精品国产乱码久久久久久男人| 中文在线观看免费www的网站 | 三级毛片av免费| 久久精品亚洲精品国产色婷小说| 一进一出好大好爽视频| 精品熟女少妇八av免费久了| 国产亚洲精品一区二区www| 免费av毛片视频| 狂野欧美白嫩少妇大欣赏| 中文亚洲av片在线观看爽| 亚洲精品粉嫩美女一区| 亚洲一区高清亚洲精品| 国产精品乱码一区二三区的特点| 香蕉av资源在线| 麻豆国产av国片精品| 成人高潮视频无遮挡免费网站| 国产私拍福利视频在线观看| 中文字幕精品亚洲无线码一区| 国语自产精品视频在线第100页| 国产黄a三级三级三级人| 午夜福利成人在线免费观看| 手机成人av网站| 18禁国产床啪视频网站| 国内精品久久久久久久电影| 国产精品98久久久久久宅男小说| 国产一区二区激情短视频| 精品久久蜜臀av无| 色哟哟哟哟哟哟| 久久精品国产亚洲av高清一级| 久久国产精品人妻蜜桃| 2021天堂中文幕一二区在线观| 欧美成人一区二区免费高清观看 | 91在线观看av| 国产欧美日韩一区二区三| 成人av在线播放网站| 精品福利观看| 精品久久久久久久毛片微露脸| 成人三级做爰电影| 日本黄大片高清| 午夜亚洲福利在线播放| 50天的宝宝边吃奶边哭怎么回事| 巨乳人妻的诱惑在线观看| 在线观看日韩欧美| 欧美日韩精品网址| 黄色视频,在线免费观看| 制服人妻中文乱码| 亚洲国产精品sss在线观看| 久久性视频一级片| 亚洲成av人片在线播放无| www.自偷自拍.com| 露出奶头的视频| 亚洲人成电影免费在线| 91大片在线观看| 在线观看66精品国产| aaaaa片日本免费| 狠狠狠狠99中文字幕| 久久久精品国产亚洲av高清涩受| 欧美大码av| 午夜激情av网站| 国产一区二区在线观看日韩 | 一区二区三区高清视频在线| 国产精品久久久久久精品电影| 91国产中文字幕| 亚洲专区字幕在线| 亚洲国产欧洲综合997久久,| 首页视频小说图片口味搜索| 动漫黄色视频在线观看| 51午夜福利影视在线观看| 国产精品野战在线观看| 久久这里只有精品中国| 国产亚洲欧美98| 手机成人av网站| 国产精品日韩av在线免费观看| 又爽又黄无遮挡网站| 精品少妇一区二区三区视频日本电影| 日韩三级视频一区二区三区| 国产av一区在线观看免费| 最好的美女福利视频网| 在线观看免费日韩欧美大片| 757午夜福利合集在线观看| 亚洲av电影在线进入| 欧美久久黑人一区二区| 国产一区在线观看成人免费| 麻豆成人午夜福利视频| 欧美乱色亚洲激情| 色综合站精品国产| 天堂影院成人在线观看| 国产精品久久久av美女十八| 很黄的视频免费| 中文资源天堂在线| 久久精品国产亚洲av香蕉五月| 国产麻豆成人av免费视频| 黑人巨大精品欧美一区二区mp4| 欧美日韩亚洲综合一区二区三区_| 久久国产精品人妻蜜桃| 少妇人妻一区二区三区视频| 国产黄色小视频在线观看| 黄色a级毛片大全视频| 欧美黑人精品巨大| 亚洲一区高清亚洲精品| 亚洲全国av大片| 亚洲欧美日韩高清专用| av福利片在线观看| 日本a在线网址| 看黄色毛片网站| 亚洲国产欧洲综合997久久,| 国产精品 欧美亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人aa在线观看| 久久精品91蜜桃| 欧美在线一区亚洲| 亚洲国产欧美网| 少妇被粗大的猛进出69影院| 精品一区二区三区四区五区乱码| 久久精品91蜜桃| 久久精品91无色码中文字幕| 又黄又爽又免费观看的视频| 精品久久久久久久毛片微露脸| 精品国内亚洲2022精品成人| 一区福利在线观看| 好男人在线观看高清免费视频| 午夜激情av网站| 欧美日韩乱码在线| 欧美久久黑人一区二区| av欧美777| 少妇的丰满在线观看| 老汉色∧v一级毛片| 色精品久久人妻99蜜桃| 无限看片的www在线观看| 色综合站精品国产| 亚洲人成网站高清观看| 亚洲激情在线av| 99热这里只有精品一区 | 熟妇人妻久久中文字幕3abv| 中文字幕熟女人妻在线| 一级毛片女人18水好多| 99精品欧美一区二区三区四区| 成人三级做爰电影| 国产亚洲欧美在线一区二区| 伦理电影免费视频| 两个人免费观看高清视频| 人妻久久中文字幕网| 国产av在哪里看| 午夜老司机福利片| 不卡一级毛片| 成人手机av| 午夜两性在线视频| 国模一区二区三区四区视频 | 国产精品免费视频内射| 国产成人精品无人区| 久久久久国产精品人妻aⅴ院| 日韩欧美免费精品| 亚洲av片天天在线观看| 一进一出好大好爽视频| 欧美一级毛片孕妇| 亚洲国产看品久久| 狂野欧美激情性xxxx| 精品久久久久久,| 欧美黄色淫秽网站| 色哟哟哟哟哟哟| 国产一区二区三区在线臀色熟女| 午夜成年电影在线免费观看| 妹子高潮喷水视频| 亚洲欧美一区二区三区黑人| 久久久久久久久免费视频了| 无遮挡黄片免费观看| 男插女下体视频免费在线播放| 一本大道久久a久久精品| 国产91精品成人一区二区三区| 国产麻豆成人av免费视频| 欧美精品亚洲一区二区| av有码第一页| 欧美3d第一页| 国产精品亚洲一级av第二区| 久久中文字幕人妻熟女| 日本在线视频免费播放| 好男人电影高清在线观看| 亚洲国产精品成人综合色| 波多野结衣巨乳人妻| 国产精品久久久人人做人人爽| 精品一区二区三区四区五区乱码| 在线观看66精品国产| 91av网站免费观看| 精品国产乱子伦一区二区三区| 国产又色又爽无遮挡免费看| 黄色片一级片一级黄色片| av福利片在线观看| 99久久无色码亚洲精品果冻| 小说图片视频综合网站| 91国产中文字幕| 每晚都被弄得嗷嗷叫到高潮| 国产真实乱freesex| 97碰自拍视频| 亚洲一区高清亚洲精品| 五月伊人婷婷丁香| 国产成人啪精品午夜网站| 久久天堂一区二区三区四区| 91老司机精品| 女生性感内裤真人,穿戴方法视频| 男男h啪啪无遮挡| 免费观看人在逋| 亚洲国产精品合色在线| 麻豆av在线久日| 国产99白浆流出| 国产麻豆成人av免费视频| 国产av在哪里看| 欧美一级毛片孕妇| 国产av一区二区精品久久| 日本一本二区三区精品| 国产一区二区激情短视频| 在线国产一区二区在线| 日韩欧美免费精品| 午夜激情av网站| 高清毛片免费观看视频网站| 在线a可以看的网站| 99在线视频只有这里精品首页| 国产精品久久久久久亚洲av鲁大| 老司机午夜福利在线观看视频| 一区福利在线观看| 欧美大码av| 欧美性猛交黑人性爽| 亚洲自偷自拍图片 自拍| 久久久久久久精品吃奶| 午夜两性在线视频| 久99久视频精品免费| 日韩成人在线观看一区二区三区| 美女扒开内裤让男人捅视频| 99热只有精品国产| 成人精品一区二区免费| 亚洲七黄色美女视频| 国产亚洲av高清不卡| 99精品久久久久人妻精品| 人人妻,人人澡人人爽秒播| 精品久久蜜臀av无| 精品国产美女av久久久久小说| 中文字幕久久专区| 露出奶头的视频| 亚洲午夜精品一区,二区,三区| 俺也久久电影网| 日韩免费av在线播放| 亚洲第一欧美日韩一区二区三区| 老熟妇仑乱视频hdxx| 成人国语在线视频|