999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

Ti60合金疲勞蠕變交互作用下蠕變應力門檻值預測

2025-06-08 00:00:00趙翰博張建宇劉浩龔愉
重慶大學學報 2025年4期
關鍵詞:模型

中圖分類號:V252 文獻標志碼:A 文章編號:1000-582X(2025)04-029-11

doi:10.11835/j.issn.1000-582X.2024.257

引用格式:,等.Ti60合金疲勞蠕變交互作用下蠕變應力門檻值預測[J].重慶大學學報,2025,48(4): 29-39.

Prediction of threshold value of creep stress in Ti60 alloy under fatigue-creep interaction

ZHAO Hanbo, ZHANG Jianyu, LIU Hao, GONG Yu (College of Aerospace Engineering, Chongqing University, Chongqing , P. R. China)

Abstract: Fatigue-creep interaction tests were conducted on Ti60 alloy specimens at . Under test conditions of maximum fatigue stress MPa and stress ratio R=0.1 , the influence of different creep stresses on the fatigue-creep behavior of thealloy was investigated.Based on thetest data and the Norton model,anovel method was proposed to improve the understanding of fatigue-creep behavior with a limited number of samples. This approach establishes a prediction model capable of estimating the threshold value of creep stress over extended periods through short-term experimental data. The model's predictions were compared with results from the maximum axial stress method and the step-loading method, yielding a relative error of less than 2% . The results show that the proposed prediction model can accurately determine the creep stress threshold under fatigue-creep interaction conditions.

Keywords: Ti60 alloy; fatigue-creep interaction; creep stress threshold; stress prediction model

近年來隨著航空航天技術的不斷發展,航空飛行器的速度取得了重大突破,更高飛行速度的服役環境對飛機結構材料的要求也隨之增加。超音速飛機在高速飛行情況下會產生大量的熱,并且在短時間內無法及時有效地散熱,同時伴隨著機械載荷,在此種工況下,傳統的鋁合金已經無法滿足高溫環境下強度需要,而高溫鈦合金因其耐腐蝕、耐高溫等特點被用于機身上溫度較高、載荷較大、工況較惡劣的構件中,在此環境中構件會受到蠕變、疲勞及二者交互作用影響,這嚴重影響其服役壽命;而疲勞蠕變二者交互作用可能會比單獨的疲勞或蠕變作用造成更大的損傷,因此高溫環境下金屬的疲勞、蠕變及二者的交互效應研究一直是理論和工程研究的關注重點。

對疲勞蠕變交互作用的研究,目前主要基于損傷力學和斷裂力學理論,通過損傷參量和損傷累積研究材料的失效機理以及裂紋萌生壽命。目前疲勞蠕變壽命的預測方法可分為三大類,即基于Manson-Coffin方程發展的方法、基于應變區分的方法(strain range partitioning,SRP),以及基于線性損傷累積準則的方法[13]。按照試驗加載方式,可分為應變控制試驗和應力控制試驗,Saad等]對P91合金進行應變控疲勞蠕變試驗,得出該材料疲勞蠕變分為非線性軟化、軟化平穩及裂紋快速擴展至斷裂的3個階段,Cristalli等認為疲勞蠕變是一種蠕變損傷增強疲勞損傷的機制,加劇了循環軟化效應,導致疲勞蠕變壽命縮短。魏峰郝玉龍]和劉洪杰對應力控試驗條件進行了研究,發現交互作用中,蠕變損傷有時會被疲勞損傷抑制;在應力控制的試驗模式下, ΔXu 等[基于Narasimhachary等[的試驗結果,建立了疲勞蠕變損傷模型,認為隨著保載時間的延長,蠕變損傷在總損傷中比例增大。學者對高溫下合金疲勞蠕變交互壽命預測開展了大量研究,但對材料疲勞-蠕變交互作用的影響因素研究較少,Zhang等研究表明,不同試驗條件下影響FGH97高溫疲勞蠕變失效的主要因素不同,并分析了與應力水平和保載時間相關的損傷發展主導因素;Liu等[23-4研究了FGH97在不同溫度下疲勞蠕變交互試驗中的主要失效因素,發現蠕變效應發揮作用存在一個溫度門檻值,在該門檻溫度之下可忽略蠕變效應。大量的研究表明,疲勞-蠕變試驗中,蠕變損傷與疲勞損傷之間存在非常復雜的關系,然而,不管是相互競爭、相互促進還是相互抑制,均是在特定的試驗條件之下,因此,也可能存在可忽略端變效應的門檻應力。鑒于此,筆者通過在 550°"下對Ti60合金試件進行疲勞端變試驗,探究疲勞變試驗中不同蠕變應力的影響,以應變量 1% 為判據判斷其是否發生蠕變效應,得出疲勞蠕變試驗中蠕變應力門檻值的預測模型,并結合升降法試驗結果進行驗證。

1疲勞蠕變試驗方案

疲勞-蠕變試驗采用圖1所示Ti60合金增材制造試驗件,根據《金屬材料蠕變-疲勞試驗方法》(GB/T38822—2020)[,在INSTRON8801疲勞試驗機上進行。試件與試驗機通過銷釘連接,試驗時保證加載裝置與試件有良好的對中性,試驗溫度為 "550° ,通過高溫爐對試件進行加熱,試驗前保溫 30min 確保試件溫度分布均勻。采用應力控制,試件的應變由高溫引伸計測得,根據飛機的飛-續-飛載荷工況,飛機構件將承受疲勞-蠕變載荷,設計試驗加載波形如圖2所示。

圖1疲勞蠕變試件示意圖 Fig.1 Fatigue creep specimen

圖2中, 為疲勞載荷的峰值與谷值, 為蠕變應力, 為疲勞載荷循環數, 為蠕變保載時間。本試驗中,保持疲勞載荷的最大應力水平為 450MPa 以及應力比 R=0.1 不變,在環境溫度 "550° 下研究蠕變應力對試件疲勞蠕變的影響。為模擬高速飛行器薄壁結構的實際工況,由循環載荷與蠕變保載組成塊譜,如圖2所示,循環加載 次,頻率為 ,每一塊譜中保載時間 0

2 試驗結果

下疲勞最大應力水平 蠕變應力 為例,試驗結果如圖3所示。從圖3中可以看出,疲勞蠕變試驗所得應變-時間曲線與純蠕變試驗中的蠕變曲線相似,其應變-時間關系有明顯的3個階段,即應變速率由快變慢的第1階段;應變速率基本保持不變的線性增長第2階段以及應變速率持續增加直至破壞的第3階段。此外,不同蠕變應力、蠕變時長為 50h 和 100h 的蠕變應變-時間曲線如圖4和圖5所示。

從圖4中可以看出,當蠕變時長為 50h ,蠕變應力不高于 390MPa 情況下的應變速率基本一致;蠕變應力高于 405MPa 時,應變速率隨蠕變應力增加而增大。據此推斷,蠕變應力在 390~405MPa 之間可能存在一個影響應變速率改變的轉折點。

從圖5中可以看到,Ti60合金在 、蠕變應力為 450MPa 蠕變時長 100h 時的應變量在 6% 以內,蠕變第2階段應變速率基本不變、第3階段應變速率緩慢穩定增長,與 等合金材料相關試驗[27-30相比,應變速率增長較緩慢,表明Ti60材料蠕變性能良好。

圖5疲勞-蠕變應變-時間曲線 (100h. 】Fig.5 Strain-time curve in fatigue-creep test( 100h

3蠕變應力門檻值預測模型

3.1基于最大軸向應力預測模型

2009年陳學東等在延性耗竭理論及損傷力學有效應力基礎上提出最大軸向應力壽命預測模型,其表達式如式(1)所示。

式中: 為疲勞蠕變壽命; c 和 m 為材料常數; 為最大軸向應力。

在等溫(給定的相同溫度)及相同保載時間的試驗條件下,利用式(1)得出基于最大軸向應力的蠕變應力門檻值預測模型,如式(2)所示。

式中: 為蠕變應力; 為試驗時長,單位為h; 為材料參數。

應用式(2)對蠕變應力門檻值的預測步驟如圖6所示。

圖6基于最大軸向應力計算蠕變應力流程圖Fig.6 Flow chart of creep stress calculation based on maximum axial stress model

按照上述流程,計算達到不同應變量下,進行 50h 試驗的蠕變應力門檻值的各材料參數(見表1),計算結果如表2所示。

3.2基于平均應變速率預測模型

3.1節所提到的預測模型步驟較復雜,計算不同應變量下的門檻值時均需要重新通過試驗數據確定相關參數。下面以平均應變速率為依據,提出步驟更為簡潔的蠕變應力門檻值預測模型。

試驗中觀察到無論試驗時長為 50h 或 ,蠕變曲線第2階段的線性應變增加速率隨蠕變應力 的變化而改變,且隨著蠕變應力 的增加,蠕變應變速率也隨之增加,二者具有明顯正相關性,在實際應用中,蠕變第1階段與第3階段時間占比小,應變速率變化快且不穩定,而蠕變第2階段在整個蠕變過程中時間占比較大,速率穩定,且不同的 對應的應變速率有明顯差異,根據疲勞蠕變過程中不同蠕變應力 對應的平均應變變化率提出疲勞蠕變試驗中蠕變應力門檻值的預測模型。

根據目前試驗所得數據,不同蠕變應力下試驗的應變速率不同,且隨著蠕變應力的增大,蠕變速率也隨之增加,故基于Norton模型[提出。

式中:t為疲勞-蠕變試驗時長; 為蠕變應力; 為擬合參數,在雙對數坐標系下進行線性擬合確定。

在式(3)運用時,假設應變速率與設計達到應變 和設計壽命 有式(4)所示關系,將式(4)帶入式(3)中可得

式中, 為蠕變應力門檻值。以式(5)模型為依據,通過試驗驗證其疲勞 蠕變過程中,蠕變應力 與疲勞-蠕變時長t的關系,具體流程如圖7所示。本試驗中,蠕變應力 與應變速率 dε/dt 在雙對數坐標下的擬合曲線如圖8所示。

圖7基于平均應變速率計算蠕變應力流程圖Fig.7 Flow chart of creep stress calculation based on average strain rate model
圖8蠕變應力-應變速率曲線 Fig.8 Creep stress-strain rate curve

在雙對數坐標下,根據式(5)對試驗數據進行擬合,可得

根據式(8)計算不同應變量下,進行 50h 試驗的蠕變應力門檻值,結果如表3所示。

表3基于平均應變速率蠕變應力門檻值計算結果Table 3 Creep stress threshold value based on average strain rate model

3.3 模型預測與試驗結果對比

為驗證3.1和3.2節中模型的有效性,用升降法測得 $550~\mathrm{{^\circC}}$ 下,疲勞最大應力水平為 450MPa ,經歷試驗時間為 50h 和 ,使試件應變達到 1% 的蠕變應力 ,用試驗實際測得的2組蠕變應力 與式(3)及式(8)所得結果進行對比,驗證其準確性和有效性。 50h 和 100h 疲勞蠕變試驗運用升降法得到的升降圖如圖9~10所示。

試驗結果與模型預測結果對比及誤差如表4所示。

由表4可知,基于最大軸向應力預測模型與本文所提出的模型均有良好的預測效果,但是基于最大軸向應力預測模型在預測蠕變門檻應力時由于其所需的試驗數據較多,無法預測本文中蠕變時長為 100h 的蠕變門檻應力。

按照3.1和3.2節的2種預測蠕變應力門檻值的模型,計算不同應變量下試驗時長 50h 的蠕變應力門檻值,結果對比如表5所示。

從圖11中可以看出,應用3.2節所提出模型計算結果與3.1節模型結果相一致。由表5可知,隨著試驗應變量的增加,2種模型的相對誤差逐漸減小,計算應變量較大的情況二者計算結果趨于一致。應用3.1節的最大軸向應力模型計算端變應力門檻值時,需要大量的數據來確定參數,且如需計算不同應變量下對應的蠕變應力

門檻值時,需重新確認公式中的常數,計算量較大;而3.2節的基于平均應變速率模型僅需對試驗數據進行一次參數確定即可獲得不同應變量下的蠕變應力門檻值,且能保證預測的準確性,減少大量計算。

圖112種模型蠕變應力門檻值計算結果對比Fig. 11 Comparison of calculation results of creep stress threshold values of two model

4結論

1)Ti60合金在 下進行疲勞端變試驗結果顯示,在試驗過程中,試件有明顯的蠕變損傷過程,即第1階段應變速率由快變慢,第2階段應變速率保持不變,以及應變速率持續增加的第3階段;以圖3的疲勞最大應力 450MPa 和蠕變應力 450MPa 為例,雖然應變-時間曲線有明顯的蠕變3個階段,但第1階段持續時間短,所占比重較小;第2階段所占比重較大且應變速率穩定;雖然有明顯的第3階段,但在第3階段中,與其他材料的疲勞蠕變試驗不同,應變速率增長比較緩慢,在 6% 應變以內無應變快速增加直至斷裂現象,說明Ti60合金具有良好且穩定的蠕變性能。

2)當固定試驗溫度、疲勞-蠕變試驗中的疲勞載荷、保載時間等因素時,材料的應變速率隨著蠕變應力的不同而不同,且蠕變應力越大,應變速率越大。根據此現象提出的疲勞蠕變預測模型的操作簡單,精度較高。

3)所提出蠕變應力門檻值的預測模型僅僅考慮蠕變應力單個因素的影響,但金屬在高溫環境下的疲勞蠕變交互作用影響因素十分復雜,從試驗結果看,疲勞和蠕變2種失效形式存在明顯的交互作用,未來的研究可以不僅僅局限于蠕變應力的影響,也將疲勞失效考慮其中,使得疲勞蠕變預測模型體系更加完善、準確。

參考文獻

[1]李欣,趙軍,劉時兵,等.航空用高溫鈦合金的研究進展[C]//2020中國鑄造活動周論文集.合肥:中國機械工程學會鑄造分會,鑄造行業生產力促進中心,2020:109-113.LiX,Zhao J,LiuSB,et al.Research progressofhigh temperature titaniumallysfor aviation[C]//22ChinaFoundryCongress.Hefei:Foundry Institution of Chinese MechanicalEnginering Society,Foundry ProductivityPromotion Center,2020:109-113. (in Chinese)

[2]何春艷,張利軍.國內外高溫鈦合金的發展與應用[J].世界有色金屬,2016(1):21-25.HeCY,Zhang LJ.Thedevelopment andapplicationofhightemperature titaniumaloyatdomestic andabroad[]WorldNonferrous Metals, 2016(1): 21-25. (in Chinese)

[3]張鴻漸.高溫鈦合金的發展與應用[J].技術與市場,2015,22(12):208.ZhangHJ.Development andapplicationofhightemperaturetitaniumaloy[J]TechnologyandMarket,2015,22(12):208.(in Chinese)

「4丨涂善東.高溫結構完整性原理IM1北京:科學出版社.2003.

Tu S D. High temperature structural integrity[M]. Beijing: Science Press, 2003. (in Chinese)

[5]張俊善.材料的高溫變形與斷裂[M].北京:科學出版社,2007. Zhang JS.Deformation and fracture ofmaterials at high temperature[M].Beijing:Science Press,207.(in Chinese)

[6]ViswanathanR,Stringer J.Failuremechanismsofhightemperaturecomponentsin powerplants[].Journalof Engineering Materials and Technology, 2000, 122(3): 246-255.

[7]龔偉忠.316H不銹鋼缺口蠕變-疲勞行為及其壽命預測方法研究[D].上海:華東理工大學,2022. Gong W Z.Researchonnotchcreep-fatiguebehaviorand life prediction methodof 316H stainless stee[D].Shanghai: East China University of Science and Technology, 2022. (in Chinese)

[8]王潤梓,廖鼎,張顯程,等.高溫結構蠕變疲勞壽命設計方法:從材料到結構[J].機械工程學報,2021,57(16):66-86,105. Wang RZ,LiaoD, Zhang XC,etal.Creep-fatigue life design methods inhigh-temperature structures:frommaterials to components[J]. Journal of Mechanical Engineering, 2021, 57(16): 66-86, 105. (in Chinese)

[9]王潤梓.基于能量密度耗散準則的蠕變-疲勞壽命預測模型及應用[D].上海:華東理工大學,2019. Wang RZ.Acree-fatiguelifepredictionmodelbasedonstrainenergy densityexhaustioncriterionanditsapplicationonaeroengine turbine discs[D]. Shanghai: East China University of Science and Technology, 2019.(in Chinese)

[10]Zhang SDTakahashiYCrep-fatiguelifeanddamageevaluationunderarousstrainwaveformsforNi-basedAloy74H[]. International Journal of Fatigue, 2023, 176: 107833.

[11]TakazawaSKangJAbeMetal.Demonstrationofsingle-framecoerentX-aydiffactionimagingusingtriangularapeue: towards dynamic nanoimaging of extended objects[J].Optics Express, 2021, 29(10):14394.

[12]TakahashiYModellngofruptureductilityofmetallcmaterialsforiderangesoftemperaturesandloadingcoditionsartI: development of basic model[J]. Materials at High Temperatures, 2020, 37(6): 357-369.

[13]TakahashiYModellingofruptureductlitofmetallcmaterialsoveriderangesoftemperaturesandloadingonditionspart II:comparison with strain energy-based approach[J].Materials atHigh Temperatures,2020,37(5):340-350.

[14]王家璇,李夢陽,鄭澤邦.蠕變-疲勞交互作用下P91鋼變形行為的研究進展[J].熱加工工藝,2024,53(9):1-7. Wang JX,LiMY,Zheng ZB.Advancesindeformation behaviour of P91 steelundercreep-fatigue interaction[J].Hot Working Technology,2024,53(9):1-7.(in Chinese)

[15]SaadAA, Bachok Z,SunW.A studyonthedamage evolutionof P9lsteelundercyclic loading athightemperature[]. International Journal of Automotive and Mechanical Engineering, 2016, 13(3): 3564-3573.

[16]CristallCAgostiniPBerardDetal.Lowcyclefatigue,crep-fatigueandrelaxatiofatiguetestsoPl]ualof Physical Science and Application, 2017, 7(2): 18-26.

[17]魏峰.P91鋼蠕變-疲勞交互作用損傷模型研究及壽命評估[D].成都:西南交通大學,2009. WeiF.Researchoncreep-fatigueinterationdamage modelandlife assessmentsofP9lsteel[D.Chengdu:Southwest Jiaotong University, 2009. (in Chinese)

[18]郝玉龍.P91鋼蠕變特性及蠕變疲勞交互作用研究[D].成都:西南交通大學,2005. Hao Y L. Study on creep and creep-fatigue interaction of P91 stee[D]. Chengdu: Southwest Jiaotong University, 2005.(in Chinese)

[19]劉洪杰.電站鍋爐用P91鋼蠕變/疲勞交互作用的試驗研究[J].動力工程,2007,27(6):990-995. Liu HJ.Experimental studyoncreep-fatigue interactionbehaviorof stelP9lforpower plant boilers[J].Journal of Power Engineering, 2007, 27(6): 990-995.(in Chinese)

[20]XuL,ZhaoLGaoZtalAnvelcreep-fatgueinteractindamagemdelwiththetressefecttosimulatethecrfatigue crack growth behavior[J]. International Journal of Mechanical Sciences, 2017,130:143-153.

[21]Narasimhachary SB,Saxena A.ResultsoftheASTMroundrobinoncreep-fatiguecrack growthtestingofaP91steel[]. Materials Performance and Characterization, 2019, 8(1): 20190125.

[22]ZhangMZhangYX,LiuHetal.Judgmentcriterionofthedominantfactorofcrep-fatiguecrackgrowthinanickelbase superalloy at elevated temperature[J]. International Journal of Fatigue, 2019, l18: 176-184.

[23]LiuH,BaoR,FeiBJ.Determinationofcreecrack growththresholdbyexperiments underelevatedtemperaturewith prestressed specimens[J]. Advanced Materials Research, 2014, 891/892: 371-376.

[24]LiuH,BaoR,LeiWMetal.Evaluatingthecriticaemperatureofcreep-fatigueinteractianickelasedpowdertallugy superalloy[J].Key Engineering Materials, 2013, 577/578: 625-628.

[25]李舜酪.機械疲勞與可靠性設計[M].北京:科學出版社,2006. Li S M.Mechanical fatigue and reliability design[M]. Beijing: Science Press, 20o6. (in Chinese)

[26]國家市場監督管理總局,國家標準化管理委員會.金屬材料蠕變-疲勞試驗方法:GB/T38822——2020[S].北京:中國標準 出版社,2020. State Administrationfor Market Regulation,National Standardizaton Administration.Metallicmaterials:creep fatigue test method: GB/T 38822—2020[S]. Beijing: Standards Press of China, 2020. (in Chinese)

[27] Liu Z, Gong J G, Zhao P, et al. Creep-fatigue interaction and damage behavior in 9-12%Cr steel under stress-controlled cycling atelevatedtemperature:effctsofholdingtimeandloadingrate[J]InternationalJournalofFatigue2256:06684.

[28]TakahashiY.Studyoncreep-fatigueevaluationproceduresforhighchromium stels:partIsensitivitytocalculated deformation[J]. International Journal of Pressure Vessels and Piping, 2008, 85(6): 423-440.

[29] Dong C L, Yu H C, Jiao ZH ,et al.Low cycle fatigue, creep and creep-fatigue interaction behavior of a TiAl alloy at high temperatures[J]. Scripta Materialia, 2018, 144: 60-63.

[30]陳凌,張賢明,歐陽平.一種疲勞-蟠變交互作用壽命預測模型及試驗驗證[J].中國機械工程,2015,26(10):1356-1361. Chen L,Zhang XMOuyang PAmodeloflifepredictionforfatigue-creepinteractionanditsexperimentalverification[]. China Mechanical Engineering, 2015, 26(10): 1356-1361. (in Chinese)

[31]陳學東,范志超,江慧豐,等.復雜加載條件下壓力容器典型用鋼疲勞蠕變壽命預測方法[J].機械工程學報,2009,45(2): 81-87. Chen XD,Fan ZC,Jiang HFetal.Cree-fatiguelifepredictionmethodsofpressureesseltypicalsteelsundercomplicated loading conditions[J]. Chinese Journal of Mechanical Engineering,2009, 45(2): 81-87.(in Chinese)

[32]姜禮益,宋世杰,闕前華.T91鋼高溫蠕變實驗和耦合損傷模型研究[J].四川輕化工大學學報(自然科學版),2023,36(1): 1-10. JiangLY,Song SJ,KanQH.Studyoncreeexperiment anddamage-coupledconstitutivemodelofT91steelatelevated temperature[J]Journal of SichuanUniversityof Scienceamp;Engineering (Natural Science Edition)2023,36(1):1-10. (in Chinese)

[33]中國航發北京航空材料研究院,郭廣平,丁傳富.航空材料力學性能檢測[M].北京:機械工業出版社,2018. AECC Beijing Institute of Aeronautical Materials, Guo GP,Ding CF.Mechanical testing of Aeronautical Materials[M]. Beijing: China Machine Press, 2018. (in Chinese)

(編輯鄭潔)

猜你喜歡
模型
一半模型
一種去中心化的域名服務本地化模型
適用于BDS-3 PPP的隨機模型
提煉模型 突破難點
函數模型及應用
p150Glued在帕金森病模型中的表達及分布
函數模型及應用
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
3D打印中的模型分割與打包
主站蜘蛛池模板: 91破解版在线亚洲| 国产福利小视频在线播放观看| 欧美区在线播放| a在线观看免费| 五月天综合网亚洲综合天堂网| 国产又粗又猛又爽| AV不卡国产在线观看| aa级毛片毛片免费观看久| 国产SUV精品一区二区| 91亚瑟视频| 亚洲成在人线av品善网好看| 亚洲黄色高清| 久久免费看片| 欧美全免费aaaaaa特黄在线| 九月婷婷亚洲综合在线| 欧美日韩一区二区在线播放| 亚洲婷婷在线视频| 亚洲成人免费看| 国产a v无码专区亚洲av| 青青青草国产| av免费在线观看美女叉开腿| 一区二区自拍| 久久狠狠色噜噜狠狠狠狠97视色 | 超清无码一区二区三区| 在线观看无码av五月花| 久久亚洲国产视频| jijzzizz老师出水喷水喷出| 国产精品免费电影| 国产精品v欧美| 国产高清在线丝袜精品一区| 欧美中文字幕第一页线路一| 一级做a爰片久久免费| 国产国产人成免费视频77777 | 国产丝袜一区二区三区视频免下载| 最近最新中文字幕在线第一页| 国产真实自在自线免费精品| 这里只有精品国产| 国产女人爽到高潮的免费视频| 制服无码网站| 欧美日韩成人| 好久久免费视频高清| 在线观看国产小视频| 亚洲国产天堂在线观看| 女高中生自慰污污网站| 国产97视频在线观看| 熟女成人国产精品视频| 国产成人欧美| 国产jizzjizz视频| 99这里只有精品在线| 亚洲伊人久久精品影院| 91在线激情在线观看| 日本午夜三级| 午夜一级做a爰片久久毛片| 亚洲成aⅴ人片在线影院八| 超级碰免费视频91| 亚洲欧美日韩成人在线| 国产成人精品综合| 国产亚洲精品yxsp| 青草视频在线观看国产| 久久久久人妻一区精品色奶水| 亚洲色图欧美一区| 欧美高清日韩| 波多野结衣的av一区二区三区| 日本人妻丰满熟妇区| 精品偷拍一区二区| 美女毛片在线| 国产欧美综合在线观看第七页| 亚洲国产AV无码综合原创| 亚洲大尺度在线| 亚洲福利视频一区二区| 亚洲成a人在线观看| 十八禁美女裸体网站| 国产人妖视频一区在线观看| 456亚洲人成高清在线| 欧美精品1区2区| 国产理论一区| 亚洲综合精品第一页| 日本精品一在线观看视频| 99一级毛片| 成人午夜视频在线| 激情综合五月网| 精品国产乱码久久久久久一区二区|