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

    殘余應(yīng)力與拘束對 GH4169合金疲勞及蠕變-疲勞裂紋擴展速率影響研究

    2024-04-12 19:27:40雷雪雪劉芳崔元元王驍曉楊杰陳浩峰
    上海理工大學(xué)學(xué)報 2024年1期

    雷雪雪 劉芳 崔元元 王驍曉 楊杰 陳浩峰

    摘要:為深入理解殘余應(yīng)力與拘束在疲勞載荷下的交互作用,以鎳基高溫合金GH4169為研究對象,選用緊湊拉伸(CT)試樣,對不同拘束CT試樣的上方施加不同大小的預(yù)加載荷從而在裂尖產(chǎn)生不同應(yīng)力。將該應(yīng)力作為殘余應(yīng)力,系統(tǒng)研究殘余應(yīng)力和拘束交互作用下GH4169合金的疲勞和蠕變—疲勞裂紋擴展速率。結(jié)果表明:隨著裂尖殘余壓應(yīng)力的增加,不同拘束下GH4169合金的疲勞和蠕變—疲勞裂紋擴展速率均降低。與低拘束試樣相比,高拘束試樣的疲勞和蠕變—疲勞裂紋擴展速率對殘余應(yīng)力的變化更加敏感,這主要與裂尖Mises應(yīng)力和垂直于裂紋擴展方向的正應(yīng)力有關(guān)。與疲勞裂紋擴展速率相比,蠕變—疲勞裂紋擴展速率對殘余應(yīng)力的變化更加敏感。

    關(guān)鍵詞:拘束;殘余應(yīng)力;疲勞裂紋擴展;蠕變—疲勞裂紋擴展;GH4169合金

    中圖分類號:TH 114;O 346 ?文獻標(biāo)志碼:A

    Study on the effect of residual stress and constraints on fatigue and creep-fatigue crack propagation rates

    LEI Xuexue1, LIU Fang2, CUI Yuanyuan2, WANG Xiaoxiao3, YANG Jie1, CHEN Haofeng3

    (1. School of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China;?2. School of Mechanical Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China;3. School of?Mechanical and Power Engineering, East China University of Science and Technology, Shanghai 200237, China)

    Abstract: In order to deeply understand the interaction between residual stress and constraint under fatigue, taking GH4169 superalloy as the research object, compact tension (CT) specimens with different constraints were selected. Different loads were applied to CT specimens to produce different stresses at the crack tip, and the stresses were introduced into the CT specimen as residual stress. And then, the fatigue and creep-fatigue crack propagation rates of GH4169 superalloy were systematically studied under the interaction of residual stress and constraint. The results showed that the fatigue and creep-fatigue crack propagation rates of GH4169 alloy decreased with the increase of residual compressive stress at crack tip. In contrast to the specimen with low constraint, the fatigue and creep-fatigue crack propagation rates of the specimen with high constraint were more sensitive to the change in residual stress, and it was mainly related to the Mises stress at the crack tip and the normal stress perpendicular to the crack propagation direction. Compared with the fatigue crack propagation rate, the creep-fatigue crack propagation rate was more sensitive to the change of residual stress.

    Keywords:constraint; residual stress; fatigue crack propagation; creep-fatigue crack propagation; GH4169 alloy

    鎳基高溫合金 GH4169因其高強度、良好的抗氧化性、穩(wěn)定的熱加工性及優(yōu)良的可焊接性等特點,在航空航天領(lǐng)域被用作制造渦輪盤的材料[1]。然而,在惡劣的服役環(huán)境和長時間工作下,渦輪盤不可避免地發(fā)生疲勞或蠕變–疲勞失效[2]。為了保障航空發(fā)動機的安全性和可靠性,需對 GH4169合金的疲勞行為和蠕變–疲勞行為展開深入研究。

    疲勞和蠕變–疲勞裂紋擴展速率是疲勞和蠕變–疲勞行為的主要指標(biāo),也是結(jié)構(gòu)長壽命設(shè)計中壽命估算的重要依據(jù),因此受到了學(xué)者們的廣泛關(guān)注。近年來,在疲勞裂紋擴展速率方面,Ye 等[3] 研究了 GH4169合金在最大應(yīng)力700 MPa、不同應(yīng)力比下的疲勞裂紋擴展行為,提出了使用裂紋尖端可逆塑性區(qū)尺寸的線性函數(shù)來預(yù)測 GH4169合金的裂紋擴展速率。Brynk 等[4]發(fā)現(xiàn)適量的稀土元素可以使疲勞裂紋擴展速率降低。在蠕變–疲勞裂紋擴展速率方面, Chen 等[5]研究了熱壓縮保載對 GH4169合金蠕變–疲勞裂紋擴展行為的影響。發(fā)現(xiàn)在高溫條件下,壓縮保載導(dǎo)致高溫下蠕變變形的積累,從而顯著加快了裂紋擴展速度。Wang 等[6] 研究了 GH4169合金的蠕變–疲勞交互行為,并用裂紋尖端張開位移表征蠕變–疲勞小裂紋擴展速率。

    殘余應(yīng)力是影響疲勞和蠕變–疲勞裂紋擴展速率的重要因素之一。為了研究和理解殘余應(yīng)力對疲勞和蠕變–疲勞裂紋擴展行為的影響, Sun 等[7]研究了深冷處理、激光沖擊處理和激光沖擊加工與后加工的混合處理方法在 GH4169合金表面產(chǎn)生的殘余應(yīng)力,結(jié)果表明,激光沖擊使試樣表面產(chǎn)生殘余壓應(yīng)力,使得裂紋擴展速度減慢,疲勞壽命提高,混合處理方法使疲勞壽命增加了2倍以上。 Zhang 等[8]對感應(yīng)淬火后的疲勞裂紋擴展行為進行了實驗研究,認為在高周疲勞下,壽命的提

    高歸結(jié)于殘余壓應(yīng)力對裂紋萌生及擴展的抑制作用。李萍萍等[9]發(fā)現(xiàn)低載強化能夠在不影響原有殘余壓應(yīng)力的基礎(chǔ)上,提高零件疲勞壽命。榮劍英[10]對緊湊拉伸(CT)試樣通過有限元模擬,制備了拉伸殘余應(yīng)力,使得試樣受拉應(yīng)力,結(jié)果表明,拉伸殘余應(yīng)力加速了蠕變–疲勞裂紋擴展速率。

    與殘余應(yīng)力一樣,裂尖拘束也是影響材料疲勞和蠕變–疲勞裂紋擴展速率的重要因素。拘束是指試樣或結(jié)構(gòu)對裂紋尖端材料塑性變形的阻礙,可分為面內(nèi)拘束和面外拘束。面內(nèi)拘束受裂紋擴展方向上試樣或結(jié)構(gòu)尺寸(如試樣寬度)的影響;面外拘束則受與裂紋前端面相平行方向上試樣或結(jié)構(gòu)尺寸(如試樣厚度)的影響。不同幾何形狀、不同厚度(即不同拘束)試樣的裂紋擴展速率也會存在差別[11-13]。一些學(xué)者觀察到高拘束會降低疲勞裂紋擴展速率[11-12],但是,也有學(xué)者觀察到了相反的現(xiàn)象[13]。在蠕變–疲勞條件下,拘束的影響同樣明顯。 Jing 等[14]研究了 CT 試樣的面內(nèi)、面外拘束和保載時間對蠕變–疲勞裂紋擴展速率的影響。結(jié)果表明,面內(nèi)和面外拘束的增加會導(dǎo)致蠕變–疲勞裂紋擴展速率加快。Lu 等[15]研究了拘束對 GH4169合金蠕變–疲勞裂紋擴展速率的影響,并使用統(tǒng)一拘束參數(shù) Ap 建立了蠕變–疲勞裂紋擴展速率與拘束的關(guān)聯(lián)。

    因機加工、熱加工等過程均會產(chǎn)生殘余應(yīng)力,而一旦有了裂紋,拘束就會存在,所以,在渦輪盤等實際結(jié)構(gòu)中,殘余應(yīng)力和拘束往往同時存在。為了更準(zhǔn)確地考察材料的疲勞與蠕變–疲勞裂紋擴展速率,需要同時納入殘余應(yīng)力與拘束的影響。然而,在目前已有的研究中,同時考慮拘束與殘余應(yīng)力對疲勞和蠕變–疲勞裂紋擴展速率影響的研究還不多。 Ma 等[16]研究了不同拘束 CT 試樣和偏心加載單邊缺口(ESET)試樣中殘余應(yīng)力的分布;榮劍英[10]研究了不同拘束 CT 試樣在拉伸殘余應(yīng)力下的蠕變–疲勞裂紋擴展速率。吳進等[17]研究了拘束與殘余應(yīng)力交互作用下異種金屬焊接接頭的斷裂行為,但在研究中尚未考慮疲勞與蠕變–疲勞的影響。

    鑒于此,本文在之前研究基礎(chǔ)上,選用渦輪盤常用材料 GH4169合金為研究對象,針對殘余應(yīng)力和拘束交互作用下的疲勞和蠕變–疲勞裂紋擴展速率進行研究,以期深入理解殘余應(yīng)力與拘束在疲勞載荷下的交互作用,得到不同拘束下疲勞和蠕變–疲勞裂紋擴展速率隨殘余應(yīng)力的變化規(guī)律,從而為航空發(fā)動機的安全可靠運行提供保障。

    1 裂紋擴展模型

    1.1 疲勞裂紋擴展模型

    采用?Yang 等[18]提出的疲勞裂紋擴展模型對疲勞裂紋擴展速率進行計算。在該模型中,采用修正的?Chaboche非線性隨動強化準(zhǔn)則和各向同性硬化規(guī)則計算瞬時應(yīng)力?應(yīng)變響應(yīng),采用修正廣義應(yīng)變幅(MGSA)疲勞損傷參數(shù)計算疲勞損傷,并基于逐周次損傷累積概念對累積疲勞損傷進行計算,前?i 個循環(huán)周次的累積疲勞損傷為

    D i)= d j) ???????(1)

    式中, d j)為第j 個循環(huán)周次產(chǎn)生的疲勞損傷。

    當(dāng)某一個單元的累積疲勞損傷達到0.99[14],該單元材料的楊氏模量降為1,此單元喪失承載能力并判定為失效,裂紋向前擴展一個單元的長度。

    1.2 蠕變–疲勞裂紋擴展模型

    采用 Lu等[15]提出的蠕變–疲勞裂紋擴展模型對蠕變–疲勞裂紋擴展速率進行計算。該模型基于逐周次概念,通過引入平均應(yīng)力,結(jié)合把應(yīng)力松弛與塑性應(yīng)變范圍及保載時間相關(guān)聯(lián)的概念,對應(yīng)變能密度耗散模型進行修正。并結(jié)合瞬時應(yīng)力和應(yīng)變響應(yīng),計算每次循環(huán)材料的蠕變損傷和疲勞損傷。將疲勞損傷和蠕變損傷相累積,可得到材料的總損傷。

    前 i 個循環(huán)周次的累積蠕變損傷

    式中: d j)為第 j 周產(chǎn)生的蠕變損傷; wc;crit 為無蠕變區(qū)的臨界應(yīng)變能密度,它反映了給定溫度下保載階段產(chǎn)生蠕變損傷的門檻值; j)為瞬態(tài)非彈性應(yīng)變能密度耗散率; t )為第 j 周的保載時間;λc 和 nc 為材料參數(shù)。

    前 i 個循環(huán)周次的總損傷

    D = D i)+ D i) ???????(3)

    與疲勞裂紋擴展模型[18]相同,當(dāng)某一個單元的累積總損傷達到0.99,判定該單元失效,裂紋向前擴展一個單元的長度。

    2 方案設(shè)計與有限元數(shù)值計算

    2.1 材料參數(shù)

    GH4169合金化學(xué)成分如表1所示[19],疲勞裂紋擴展模型和蠕變–疲勞裂紋擴展模型中所需材料參數(shù)如表2所示[19],這些參數(shù)均通過實驗進行了標(biāo)定,并進行了進一步的驗證[15, 18]。表2中:E 為彈性模量;ν為泊松比; Q0為初始屈服應(yīng)力; K 和 n 是代表材料黏性特性的材料參數(shù);ζ i 和 ri 為每一個背應(yīng)力部分的材料參數(shù);φ1,φ2和ω是由不同載荷工況下蠕變–疲勞試驗確定的溫度相關(guān)參數(shù); Qsa 為快速軟化第一階段各項同性變形抗力的漸近值; b 表示逼近漸近值的速度參數(shù);H為線性軟化第二階段的斜率;τf(′)為剪切疲勞強度;σf(′)為疲勞強度系數(shù);γf(′)為剪切延性強度; b0,c0表示關(guān)于疲勞強度和延性的指數(shù); G 為剪切模量;φ1和 n1分別為材料和溫度相關(guān)的模型參數(shù); n2為穩(wěn)態(tài)蠕變指數(shù); A , B是與材料有關(guān)的常數(shù); wf,crit 為臨界失效能應(yīng)變密度。

    2.2 試樣設(shè)計

    選擇 CT 試樣在不同殘余應(yīng)力、不同拘束下的疲勞和蠕變–疲勞裂紋擴展行為進行研究。 CT 試樣如圖1所示,其中,試樣寬度 W=32 mm ,初始裂紋長度為 a。通過應(yīng)變ε控制加載,疲勞和蠕變–疲勞加載波形如圖2所示。加載應(yīng)變幅值均為1%,蠕變–疲勞加載時保載時間為100 s。

    通過改變 CT 試樣的初始裂紋長度(a=9.6 mm, a/ W=0.3; a=16 mm , a/ W=0.5)改變拘束;為了在裂紋尖端得到不同的殘余應(yīng)力,在疲勞或蠕變–疲勞加載前,分別在試樣上下兩面施加不同預(yù)加載荷(0,10,20,50,100,150 MPa),以在裂尖產(chǎn)生不同的應(yīng)力。將該應(yīng)力作為殘余應(yīng)力引入 CT 試樣(僅引入應(yīng)力,不引入應(yīng)變),而后再進行疲勞或蠕變–疲勞加載。

    2.3 殘余應(yīng)力的引入及有限元數(shù)值計算

    采用有限元軟件 ABAQUS 對 CT 試樣建模,由于加載和 CT 試樣幾何的對稱性,僅建立二分之一模型。在試樣加載孔中心建立參考點,參考點和試樣內(nèi)表面剛性耦合,對參考點施加疲勞或蠕變–疲勞載荷,對裂尖前端的韌帶部分施加對稱邊界條件,如圖1所示。模型網(wǎng)格采用8節(jié)點雙線性減縮積分平面應(yīng)變單元(CPE8R),最小網(wǎng)格尺寸為0.1 mm×0.1 mm,網(wǎng)格無關(guān)性已在之前的工作中[18]進行了驗證。

    將建立好的 CT 試樣模型 Model 復(fù)制為模型 Model-copy,通過在 Model 上施加不同大小的均布載荷(0,10,20,50,100,150 MPa)而得到裂尖應(yīng)力場,將此應(yīng)力場導(dǎo)入 Model-copy,作為試樣受疲勞載荷前存在的殘余應(yīng)力。最后,對含有殘余應(yīng)力的模型 Model-copy 進行疲勞或蠕變–疲勞加載。在加載過程中,疲勞裂紋擴展模型和蠕變疲勞裂紋擴展模型將分別通過子程序的方式嵌入模型 Model-copy,從而最終獲得納入不同殘余應(yīng)力時不同拘束下的疲勞和蠕變–疲勞裂紋擴展行為。

    3 殘余應(yīng)力對不同拘束下 GH4169合金疲勞裂紋擴展速率的影響

    在不同預(yù)加載荷下,不同拘束 CT 試樣在第一周循環(huán)時的滯后回線如圖3所示。由圖可見,隨著預(yù)加載荷的增加,滯后回線起始段呈現(xiàn)逐漸下移的趨勢,尤其是滯后回線的起始點越來越低。在沒有預(yù)加載荷時,起始點所對應(yīng)的應(yīng)力為0 MPa,隨著預(yù)加載荷的增加,起始點所對應(yīng)的應(yīng)力逐漸降低。這說明預(yù)加載荷所產(chǎn)生的應(yīng)力已成功地作為殘余應(yīng)力引入到模型中,將在后續(xù)對 GH4169合金的疲勞裂紋擴展行為產(chǎn)生影響。在不同的預(yù)加載荷下,在裂尖產(chǎn)生了殘余壓應(yīng)力(應(yīng)力為負值),殘余壓應(yīng)力與預(yù)加載荷呈正相關(guān)。與低拘束(a/W=0.3)相比,在高拘束(a/W =0.5)下,隨著預(yù)加載荷的增加,起始點所對應(yīng)的應(yīng)力降低趨勢更加明顯。此外,滯后回線圍繞區(qū)域的面積隨著預(yù)加載荷的增加而減小,這說明隨著殘余壓應(yīng)力的增加,每次循環(huán)所消耗的能量有所減少。

    在不同預(yù)加載荷下,不同拘束 CT 試樣疲勞裂紋擴展長度 a 與疲勞循環(huán)周次 Nf 的關(guān)系如圖4所 示。由圖可見,無論是在低拘束還是在高拘束下,隨著預(yù)加載荷的增加,裂紋擴展至同一長度所經(jīng)歷的循環(huán)周次均明顯增多。在低拘束下,預(yù)加載荷為0,10,20,50,100,150 MPa 時裂紋擴展至1 mm所需的循環(huán)周次分別為520,524,529,548,602,665周;在高拘束下,預(yù)加載荷為0,10,20,50,100,150 MPa 時裂紋擴展至1 mm 所需的循環(huán)周次分別為1011,1059,1120,1392,1525,2076周。

    這說明隨著預(yù)加載荷的增加,不同拘束試樣的疲勞裂紋擴展速率呈現(xiàn)降低的趨勢,且在高拘束下疲勞裂紋擴展速率的變化更加明顯。以裂紋擴展至1 mm 為例,在低拘束下施加150 MPa 預(yù)加載荷時較無預(yù)加載荷時疲勞裂紋擴展速率降低27.89%,在高拘束下施加150 MPa 預(yù)加載荷時較無預(yù)加載荷時疲勞裂紋擴展速率降低105.34%。由圖5擴展至1 mm 時的云圖也可以發(fā)現(xiàn),無預(yù)加載荷時試樣的損傷范圍明顯大于預(yù)加載荷為150 MPa 時試樣的損傷范圍,且在同樣的預(yù)加載荷下,低拘束試樣損傷范圍比高拘束試樣更大。這些均與圖3中滯后回線表現(xiàn)出來的規(guī)律一致。

    提取不同拘束試樣中不同預(yù)加載荷所引起的裂尖 Mises應(yīng)力,該應(yīng)力也是疲勞裂紋擴展前裂尖的殘余應(yīng)力,如圖6所示。由圖可見,隨著預(yù)加載荷的增加,裂尖 Mises應(yīng)力呈現(xiàn)逐漸增加的趨勢,即疲勞裂紋擴展前試樣中的殘余應(yīng)力逐漸增加。且 a/W=0.5高拘束試樣裂尖的 Mises應(yīng)力較 a/W =0.3低拘束試樣更大。圖7的裂尖應(yīng)力場更加清晰地展示了這一現(xiàn)象,當(dāng)預(yù)加載荷為50 MPa 時,高拘束試樣裂尖的 Mises應(yīng)力明顯高于低拘束試樣。這與圖5所呈現(xiàn)的變化規(guī)律一致,是高拘束試樣損傷范圍低于低拘束試樣的根本原因。

    為進一步分析殘余應(yīng)力在裂尖附近的分布情況,提取了沿裂紋擴展方向上裂尖前方各點的壓應(yīng)力 S22,如圖8所示??梢园l(fā)現(xiàn),最大殘余壓應(yīng)力分布在裂尖附近,從裂尖開始,沿裂紋擴展方向, S22呈現(xiàn)遞減的趨勢,在距離裂尖較遠處逐漸趨于平緩,預(yù)加載荷所起的作用減小。在高預(yù)加載荷下,裂紋擴展路徑上分布的殘余應(yīng)力比低預(yù)加載荷下更高。

    4 殘余應(yīng)力對不同拘束下 GH4169合金蠕變–疲勞裂紋擴展速率的影響

    承受蠕變–疲勞載荷時,不同拘束 CT 試樣在第一周循環(huán)時的滯后回線如圖9所示。由圖可見,與疲勞加載時(圖3)一致,隨著預(yù)加載荷的增加,滯后回線起始段逐漸下移,起始點越來越低,所對應(yīng)的應(yīng)力逐漸降低。這是由于無論是疲勞加載還是蠕變–疲勞加載,所引入的殘余應(yīng)力是一致的,圖6~8所展示的數(shù)據(jù)和規(guī)律同樣存在于蠕變–疲勞載荷加載時,此處不再贅述。

    值得一提的是,在承受蠕變–疲勞載荷時,滯后回線圍繞區(qū)域的面積同樣隨著預(yù)加載荷的增加而減小,即殘余壓應(yīng)力的存在同樣降低了蠕變–疲勞載荷下每次循環(huán)所消耗的能量。

    在不同預(yù)加載荷下,不同拘束 CT 試樣蠕變–疲勞裂紋擴展長度 a 與蠕變–疲勞循環(huán)周次 Nc-f 的關(guān)系如圖10所示。由圖可見,與圖4一致,在不同拘束下,隨著預(yù)加載荷的增加,裂紋擴展至同一長度所經(jīng)歷的循環(huán)周次均增多。在低拘束下,預(yù)加載荷為0,10,20,50,100和150 MPa 時裂紋擴展至1 mm所需的循環(huán)周次分別為24,24,24,28,31周;在高拘束下,預(yù)加載荷為0,10,20,50,100,150 MPa 時裂紋擴展至1 mm 所需的循環(huán)周次分別為36,38,40,60,67,74周。即試樣的蠕變–疲勞裂紋擴展速率隨著殘余壓應(yīng)力的增加而降低,且在高拘束下降低更加明顯。同樣以裂紋擴展至1 mm ,150 MPa 預(yù)加載荷為例,在低拘束下蠕變–疲勞裂紋擴展速率較無預(yù)加載荷時降低29.17%,在高拘束下蠕變–疲勞裂紋擴展速率較無預(yù)加載荷時降低105.56%。由圖11裂紋擴展至1 mm 時的損傷云圖同樣可見,在殘余壓應(yīng)力下,試樣損傷范圍明顯減小,且在高拘束下,損傷范圍減小更多。

    綜合比較疲勞裂紋擴展行為和蠕變–疲勞裂紋擴展行為,可以發(fā)現(xiàn),其相同點在于在殘余壓應(yīng)力的作用下,不同拘束試樣的疲勞與蠕變–疲勞裂紋擴展速率均明顯下降;其不同點在于在同一預(yù)加載荷下,裂紋擴展至同一長度時,蠕變–疲勞裂紋擴展所需要的循環(huán)周次遠少于疲勞裂紋擴展。蠕變–疲勞時滯后回線所圍繞區(qū)域面積更大, GH4169合金的裂紋擴展速度更快,壽命更低,對殘余應(yīng)力的變化更加敏感。與低拘束相比,高拘束下疲勞和蠕變–疲勞裂紋擴展速率對殘余應(yīng)力的變化均更加敏感。

    5 結(jié) 論

    a.隨著裂尖殘余壓應(yīng)力的增加,在不同拘束下 GH4169合金的疲勞和蠕變–疲勞裂紋擴展速率均降低,每次疲勞循環(huán)和每次蠕變–疲勞循環(huán)所消耗的能量均減少。

    b.與疲勞載荷相比,在蠕變–疲勞載荷下, GH4169合金的裂紋擴展速率更快,壽命更低,且對殘余應(yīng)力的變化更加敏感。

    c.與低拘束相比,高拘束下 GH4169合金的疲勞和蠕變–疲勞裂紋擴展速率對殘余應(yīng)力的變化更加敏感,這主要與裂尖 Mises應(yīng)力和垂直子裂紋擴展方向的正應(yīng)力有關(guān)。

    參考文獻:

    [1] AZARBARMAS M, AGHAIE-KHAFRI M, CABRERA J M, et al. Dynamic recrystallization mechanisms and twining evolution during hot deformation of Inconel 718[J]. Materials Science and Engineering: A, 2016, 678:137–152.

    [2] CINI A, IRVING P E. Development of fatigue cracks from mechanically machined scratches on 2024-T351 aluminium alloy —part I: experimentation and fractographic analysis [J]. Fatigue & Fracture of Engineering Materials & Structures, 2017, 40(5):776–789.

    [3] YE S, GONG J G, ZHANG X C, et al. Effect of stress ratio on the fatigue crack propagation behavior of the nickel- based GH4169 alloy[J]. Acta Metallurgica Sinica (English Letters), 2017, 30(9):809–821.

    [4] BRYNK T, PAKIELA Z, LUDWICHOWSKA K, et al. Fatigue crack growth rate and tensile strength of Re modified Inconel 718 produced by means of selective laser melting[J]. Materials Science and Engineering: A, 2017, 698:289–301.

    [5] CHEN X, PETTIT R G, DUDZINSKI D, et al. On the role of crack tip creep deformation in hot compressive dwell fatigue crack growth acceleration in aluminum and nickel engine alloys[J]. International Journal of Fatigue, 2020, 145:106082.

    [6] WANG Z, WU W W, LIANG J C, et al. Creep –fatigue interaction behavior of nickel-based single ?crystal superalloy at high temperature by in-situ SEM observation[J]. International Journal of Fatigue, 2020, 141:105879.

    [7] SUN J F, SU A P, WANG T M, et al. Effect of laser shock processing with post-machining and deep cryogenic treatment on fatigue life of GH4169 super alloy[J]. International Journal of Fatigue, 2019, 119:261–267.

    [8] ZHANG S J, XIE J J, JIANG Q Q, et al. Fatigue crack growth behavior in gradient microstructure of hardened surface layer for an axle steel[J]. Materials Science and Engineering: A, 2017, 700:66–74.

    [9]李萍萍, 盧曦, 鄭松林.低載強化對圓柱齒輪殘余應(yīng)力影響的試驗研究[J].上海理工大學(xué)學(xué)報 , 2008, 30(2):139–142.

    [10]榮劍英.殘余應(yīng)力對?G115鋼650℃蠕變?疲勞裂紋擴展速率的作用[D].天津:天津大學(xué), 2018.

    [11] VARFOLOMEEV I, LUKE M, MOROZ S. Experimental and numerical investigations of fatigue crack growth in various specimen geometries[J]. Procedia Engineering, 2010, 2(1):1829–1837.

    [12] HUTA? P, SEITL S, KN?SL Z. Effect of constraint on fatigue crack propagation near threshold in medium carbon steel[J]. Computational Materials Science, 2006, 37(1/2):51–57.

    [13] TONG J. T-stress and its implications for crack growth[J]. Engineering ?Fracture ?Mechanics, ?2002, ?69(12):1325–1337.

    [14] JING H Y, SU D B, XU L Y, et al. Finite element simulation of creep-fatigue crack growth behavior for P91 steel at 625°C considering creep-fatigue interaction[J]. International Journal of Fatigue, 2017, 98:41–52.

    [15] LU R S, TAN J P, YANG J, et al. A new creep-fatigue crack growth model and a correlation of the creep-fatigue crack growth rate with unified constraint parameter[J]. International Journal of Fatigue, 2023, 166:107248.

    [16] MA Y E, STARON P, FISCHER T, et al. Size effects on residual stress and fatigue crack growth in friction stir welded 2195-T8 aluminium—Part I: Experiments[J]. International Journal of Fatigue, 2011, 33(11):1417–1425. [17]吳進 , 楊杰 , 陳浩峰.納入殘余應(yīng)力時不同拘束下 DMWJ 的斷裂行為[J].金屬學(xué)報, 2022, 58(7):956–964.

    [18] YANG J, GUO G L, LU R S, et al. Numerical modelling of a new FCP model and a correlation of the FCP rate with the constraint[J]. International Journal of Fatigue, 2022, 163:107036.

    [19] WANG R Z, GUO S J, CHEN H F, et al. Multi-axial creep- fatigue life prediction considering history-dependent damage evolution: a new numerical procedure and experimental validation[J]. Journal of the Mechanics and Physics of Solids, 2019, 131:313–336.

    (編輯:石 瑛)

    成人18禁在线播放| 国产成人精品无人区| 18禁国产床啪视频网站| av天堂在线播放| 中文亚洲av片在线观看爽| 亚洲av片天天在线观看| 脱女人内裤的视频| 人妻久久中文字幕网| 99riav亚洲国产免费| 中文字幕另类日韩欧美亚洲嫩草| 宅男免费午夜| 99re在线观看精品视频| 一级a爱片免费观看的视频| 免费在线观看影片大全网站| 精品福利观看| 久久久久免费精品人妻一区二区 | 亚洲av片天天在线观看| 最近最新中文字幕大全电影3 | 怎么达到女性高潮| 亚洲成国产人片在线观看| 麻豆国产av国片精品| 欧美黑人巨大hd| 露出奶头的视频| 亚洲熟妇中文字幕五十中出| 日韩欧美一区二区三区在线观看| 精华霜和精华液先用哪个| 此物有八面人人有两片| 欧美一区二区精品小视频在线| 亚洲最大成人中文| 国产成人啪精品午夜网站| 亚洲成a人片在线一区二区| 日韩高清综合在线| 99riav亚洲国产免费| 丝袜美腿诱惑在线| 成人亚洲精品一区在线观看| 国产精品一区二区三区四区久久 | 亚洲无线在线观看| 18禁黄网站禁片免费观看直播| 少妇粗大呻吟视频| 每晚都被弄得嗷嗷叫到高潮| 成人av一区二区三区在线看| 妹子高潮喷水视频| 夜夜看夜夜爽夜夜摸| 999久久久精品免费观看国产| 色老头精品视频在线观看| av免费在线观看网站| 国产免费男女视频| 国产精品乱码一区二三区的特点| 热re99久久国产66热| 怎么达到女性高潮| 淫妇啪啪啪对白视频| 日韩欧美国产在线观看| 欧美久久黑人一区二区| 首页视频小说图片口味搜索| 好男人在线观看高清免费视频 | 亚洲av片天天在线观看| 国产精品野战在线观看| 日日爽夜夜爽网站| 看片在线看免费视频| 怎么达到女性高潮| 一级毛片高清免费大全| 啦啦啦韩国在线观看视频| 免费观看精品视频网站| 久久 成人 亚洲| 精品高清国产在线一区| 亚洲av片天天在线观看| 又黄又粗又硬又大视频| 欧美日本视频| 国产精品野战在线观看| 精品国产国语对白av| 波多野结衣巨乳人妻| 久久久久精品国产欧美久久久| √禁漫天堂资源中文www| 成人国语在线视频| 亚洲精品国产一区二区精华液| 伦理电影免费视频| 天堂影院成人在线观看| 国产aⅴ精品一区二区三区波| 校园春色视频在线观看| 精品久久久久久久久久久久久 | 亚洲av五月六月丁香网| 激情在线观看视频在线高清| 国产精品1区2区在线观看.| 激情在线观看视频在线高清| 欧美不卡视频在线免费观看 | 国语自产精品视频在线第100页| 两性午夜刺激爽爽歪歪视频在线观看 | 免费在线观看完整版高清| 91麻豆精品激情在线观看国产| 国产主播在线观看一区二区| 亚洲成人精品中文字幕电影| 亚洲欧美日韩无卡精品| 亚洲第一av免费看| 国产91精品成人一区二区三区| 88av欧美| 国产一区二区激情短视频| 欧美午夜高清在线| 免费在线观看完整版高清| 亚洲熟妇中文字幕五十中出| 亚洲国产欧洲综合997久久, | 亚洲avbb在线观看| 亚洲人成网站高清观看| 日本熟妇午夜| 欧美亚洲日本最大视频资源| 欧美午夜高清在线| 变态另类丝袜制服| 久久中文字幕人妻熟女| 母亲3免费完整高清在线观看| 1024视频免费在线观看| 国产成人精品无人区| 日本 欧美在线| 亚洲精华国产精华精| 亚洲五月天丁香| 午夜福利在线观看吧| 日韩大码丰满熟妇| 成年人黄色毛片网站| 国产一卡二卡三卡精品| 久久国产精品男人的天堂亚洲| ponron亚洲| 2021天堂中文幕一二区在线观 | 欧美日韩一级在线毛片| 日韩一卡2卡3卡4卡2021年| svipshipincom国产片| 午夜久久久久精精品| 50天的宝宝边吃奶边哭怎么回事| 午夜免费鲁丝| 成人一区二区视频在线观看| netflix在线观看网站| 久久精品国产综合久久久| 好看av亚洲va欧美ⅴa在| 精品福利观看| ponron亚洲| 中文字幕人成人乱码亚洲影| 一区二区三区精品91| cao死你这个sao货| 最近最新免费中文字幕在线| 欧美日韩福利视频一区二区| 人人妻人人澡人人看| 国产亚洲欧美精品永久| 免费一级毛片在线播放高清视频| 一本大道久久a久久精品| 久久香蕉精品热| 亚洲精品在线美女| 久久久久久人人人人人| 97人妻精品一区二区三区麻豆 | 国产亚洲精品久久久久5区| 高清毛片免费观看视频网站| e午夜精品久久久久久久| 午夜免费鲁丝| 又黄又爽又免费观看的视频| 在线观看一区二区三区| 国产黄色小视频在线观看| 很黄的视频免费| 久久久国产成人精品二区| 黄色丝袜av网址大全| 日日爽夜夜爽网站| 免费无遮挡裸体视频| 一进一出抽搐gif免费好疼| 亚洲自拍偷在线| avwww免费| 国产亚洲欧美在线一区二区| 最近最新中文字幕大全免费视频| 一本一本综合久久| 成熟少妇高潮喷水视频| 淫妇啪啪啪对白视频| 久久精品亚洲精品国产色婷小说| 亚洲熟女毛片儿| 午夜免费成人在线视频| 色综合亚洲欧美另类图片| 视频区欧美日本亚洲| 中国美女看黄片| 成年免费大片在线观看| 日本三级黄在线观看| 国产激情偷乱视频一区二区| 色老头精品视频在线观看| 91字幕亚洲| 精品国产美女av久久久久小说| 黑丝袜美女国产一区| 亚洲五月天丁香| xxxwww97欧美| 亚洲人成网站高清观看| 国产精品乱码一区二三区的特点| 国产一级毛片七仙女欲春2 | 成人国语在线视频| 久久国产亚洲av麻豆专区| 亚洲一区中文字幕在线| 国产男靠女视频免费网站| avwww免费| 丝袜美腿诱惑在线| 午夜精品在线福利| 97超级碰碰碰精品色视频在线观看| 看免费av毛片| 啦啦啦韩国在线观看视频| 又紧又爽又黄一区二区| 精品人妻1区二区| 亚洲激情在线av| 性色av乱码一区二区三区2| 亚洲五月婷婷丁香| 亚洲第一欧美日韩一区二区三区| 精品久久久久久久久久免费视频| 成人永久免费在线观看视频| 黄网站色视频无遮挡免费观看| 久久精品人妻少妇| 宅男免费午夜| 国产精品精品国产色婷婷| 最近在线观看免费完整版| 久久久久精品国产欧美久久久| 女人爽到高潮嗷嗷叫在线视频| 欧美不卡视频在线免费观看 | 亚洲国产欧美日韩在线播放| 国产精品久久久av美女十八| 女生性感内裤真人,穿戴方法视频| 国产黄色小视频在线观看| 久9热在线精品视频| 成人国语在线视频| 中文字幕久久专区| 91九色精品人成在线观看| 亚洲人成网站在线播放欧美日韩| 无遮挡黄片免费观看| 亚洲精品国产精品久久久不卡| 大型av网站在线播放| 精品无人区乱码1区二区| 50天的宝宝边吃奶边哭怎么回事| 国产精品国产高清国产av| 男人舔奶头视频| 色综合亚洲欧美另类图片| 免费观看人在逋| 真人一进一出gif抽搐免费| 可以免费在线观看a视频的电影网站| 欧美久久黑人一区二区| 亚洲av五月六月丁香网| 国产av在哪里看| 亚洲国产欧洲综合997久久, | 精品卡一卡二卡四卡免费| 国产亚洲av高清不卡| 久久人人精品亚洲av| 亚洲精品国产精品久久久不卡| 亚洲五月天丁香| av片东京热男人的天堂| 日日夜夜操网爽| 国产精品一区二区精品视频观看| 亚洲精品一区av在线观看| 国产成年人精品一区二区| 精品免费久久久久久久清纯| 好男人在线观看高清免费视频 | 久久国产精品男人的天堂亚洲| 国产av不卡久久| 搡老妇女老女人老熟妇| 欧美日韩一级在线毛片| 久久热在线av| 又紧又爽又黄一区二区| 亚洲欧美日韩无卡精品| 91成年电影在线观看| 这个男人来自地球电影免费观看| 亚洲精品国产一区二区精华液| 日韩精品中文字幕看吧| 亚洲国产精品sss在线观看| 琪琪午夜伦伦电影理论片6080| 18禁裸乳无遮挡免费网站照片 | 欧美精品啪啪一区二区三区| 999久久久国产精品视频| 成人亚洲精品一区在线观看| 日本撒尿小便嘘嘘汇集6| 欧美乱色亚洲激情| 午夜免费观看网址| 午夜老司机福利片| 亚洲一区中文字幕在线| 亚洲av美国av| av片东京热男人的天堂| 国产乱人伦免费视频| 久久精品国产亚洲av香蕉五月| 成在线人永久免费视频| 大香蕉久久成人网| 亚洲天堂国产精品一区在线| 两性午夜刺激爽爽歪歪视频在线观看 | 99热6这里只有精品| 国产精品野战在线观看| 欧美亚洲日本最大视频资源| 两个人免费观看高清视频| 亚洲aⅴ乱码一区二区在线播放 | 国产97色在线日韩免费| 久久久久久久久中文| 天天躁夜夜躁狠狠躁躁| 国产午夜精品久久久久久| 黄片大片在线免费观看| 亚洲av电影不卡..在线观看| 久久久久国产精品人妻aⅴ院| 国产成人系列免费观看| 一个人免费在线观看的高清视频| 国产麻豆成人av免费视频| 黑人操中国人逼视频| 欧美av亚洲av综合av国产av| 久久久国产精品麻豆| 精品熟女少妇八av免费久了| 国产激情欧美一区二区| 少妇粗大呻吟视频| 久久香蕉国产精品| 国产精品美女特级片免费视频播放器 | 好看av亚洲va欧美ⅴa在| 999久久久精品免费观看国产| 欧美午夜高清在线| 男人的好看免费观看在线视频 | 国产v大片淫在线免费观看| 在线观看舔阴道视频| 国产精品综合久久久久久久免费| 欧美日本亚洲视频在线播放| 亚洲性夜色夜夜综合| 很黄的视频免费| 超碰成人久久| 99re在线观看精品视频| 国产亚洲精品久久久久5区| 久久久国产成人精品二区| 高潮久久久久久久久久久不卡| 国产精品综合久久久久久久免费| 亚洲av日韩精品久久久久久密| 欧美日韩福利视频一区二区| 一级毛片女人18水好多| 黑人巨大精品欧美一区二区mp4| 久久99热这里只有精品18| 老鸭窝网址在线观看| 国产又黄又爽又无遮挡在线| 欧美三级亚洲精品| 成人免费观看视频高清| 国产亚洲欧美98| 国产成+人综合+亚洲专区| 国产又爽黄色视频| 日日干狠狠操夜夜爽| 日韩欧美在线二视频| 欧美三级亚洲精品| 变态另类丝袜制服| 亚洲人成伊人成综合网2020| 国产成人精品久久二区二区免费| 国产男靠女视频免费网站| 香蕉丝袜av| 欧洲精品卡2卡3卡4卡5卡区| av免费在线观看网站| 亚洲精品久久成人aⅴ小说| 高清毛片免费观看视频网站| 人妻丰满熟妇av一区二区三区| 午夜两性在线视频| 十八禁网站免费在线| 十分钟在线观看高清视频www| 日日夜夜操网爽| 成人三级黄色视频| 一本久久中文字幕| 精品国内亚洲2022精品成人| 欧美亚洲日本最大视频资源| 亚洲久久久国产精品| av福利片在线| 一本大道久久a久久精品| 男女做爰动态图高潮gif福利片| 母亲3免费完整高清在线观看| 日本在线视频免费播放| 757午夜福利合集在线观看| 亚洲第一青青草原| 后天国语完整版免费观看| 日韩大尺度精品在线看网址| 女人高潮潮喷娇喘18禁视频| 久久久久久久精品吃奶| 精品午夜福利视频在线观看一区| 国产精品1区2区在线观看.| 久久久久久久久久黄片| 国产国语露脸激情在线看| 亚洲国产日韩欧美精品在线观看 | 国产高清激情床上av| 亚洲精华国产精华精| 国产成人av激情在线播放| 国产精品日韩av在线免费观看| 免费看十八禁软件| 免费一级毛片在线播放高清视频| av天堂在线播放| 满18在线观看网站| 成人18禁高潮啪啪吃奶动态图| 此物有八面人人有两片| 99国产精品一区二区蜜桃av| 日韩欧美国产在线观看| 悠悠久久av| 叶爱在线成人免费视频播放| 91麻豆av在线| 动漫黄色视频在线观看| 一个人免费在线观看的高清视频| 老司机午夜福利在线观看视频| 国产精品99久久99久久久不卡| 日本精品一区二区三区蜜桃| 久久久精品国产亚洲av高清涩受| 熟妇人妻久久中文字幕3abv| 两个人免费观看高清视频| 一级作爱视频免费观看| 久久久久免费精品人妻一区二区 | а√天堂www在线а√下载| 人妻久久中文字幕网| 麻豆成人av在线观看| 变态另类丝袜制服| 欧美黄色淫秽网站| 侵犯人妻中文字幕一二三四区| 亚洲熟妇中文字幕五十中出| 老司机午夜福利在线观看视频| 久久久久久免费高清国产稀缺| 视频区欧美日本亚洲| 亚洲第一电影网av| 好男人电影高清在线观看| 亚洲av熟女| 一级毛片女人18水好多| 精品人妻1区二区| 丝袜人妻中文字幕| 久久国产亚洲av麻豆专区| 国产一区在线观看成人免费| 亚洲色图 男人天堂 中文字幕| 中文资源天堂在线| 亚洲男人天堂网一区| 在线av久久热| 国产午夜福利久久久久久| 天天添夜夜摸| 法律面前人人平等表现在哪些方面| 不卡一级毛片| 91av网站免费观看| 中文字幕最新亚洲高清| 国产三级在线视频| 最近最新免费中文字幕在线| 99国产精品99久久久久| 午夜福利欧美成人| 国产极品粉嫩免费观看在线| 免费无遮挡裸体视频| АⅤ资源中文在线天堂| 色av中文字幕| 99久久久亚洲精品蜜臀av| 男女之事视频高清在线观看| 人妻丰满熟妇av一区二区三区| 国产精品久久久av美女十八| 又大又爽又粗| 国产成+人综合+亚洲专区| 免费在线观看完整版高清| 亚洲第一欧美日韩一区二区三区| 国产亚洲精品久久久久久毛片| 久久久久久久午夜电影| 女警被强在线播放| 美国免费a级毛片| 搡老岳熟女国产| 少妇裸体淫交视频免费看高清 | 久久人妻av系列| xxxwww97欧美| 日韩中文字幕欧美一区二区| а√天堂www在线а√下载| 国产精品电影一区二区三区| 日韩欧美一区视频在线观看| 中文字幕人妻熟女乱码| 精品少妇一区二区三区视频日本电影| 国产精品av久久久久免费| 久久热在线av| 亚洲片人在线观看| 老司机午夜十八禁免费视频| 999精品在线视频| 色播亚洲综合网| 成人永久免费在线观看视频| 欧美黑人欧美精品刺激| 国产片内射在线| 国产精品香港三级国产av潘金莲| 99久久精品国产亚洲精品| 亚洲精品中文字幕在线视频| 欧美日韩亚洲综合一区二区三区_| 日本一本二区三区精品| 搡老岳熟女国产| 一进一出抽搐动态| 久久精品91蜜桃| 亚洲一码二码三码区别大吗| 91在线观看av| 免费高清视频大片| www.精华液| 午夜福利在线在线| 亚洲成av人片免费观看| 一二三四在线观看免费中文在| 日韩视频一区二区在线观看| 久久热在线av| 日本 av在线| 久久久国产成人精品二区| 久久婷婷成人综合色麻豆| 国产成人精品久久二区二区91| 国产成人精品久久二区二区免费| 91麻豆精品激情在线观看国产| 精品一区二区三区四区五区乱码| 一区二区三区激情视频| av在线播放免费不卡| 中文字幕精品亚洲无线码一区 | 999精品在线视频| 狠狠狠狠99中文字幕| 女人爽到高潮嗷嗷叫在线视频| 国产精品二区激情视频| 精品午夜福利视频在线观看一区| 老司机在亚洲福利影院| 国产一区二区激情短视频| svipshipincom国产片| 亚洲精品中文字幕在线视频| 欧美性猛交╳xxx乱大交人| 国产激情欧美一区二区| 国产午夜精品久久久久久| 久久久久久久久免费视频了| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美98| 日韩免费av在线播放| 桃红色精品国产亚洲av| 日本成人三级电影网站| 搡老妇女老女人老熟妇| 国产激情偷乱视频一区二区| 中文字幕人妻丝袜一区二区| 国产人伦9x9x在线观看| 午夜福利欧美成人| 最新在线观看一区二区三区| 国产精品影院久久| 国产国语露脸激情在线看| 一区二区三区高清视频在线| 18美女黄网站色大片免费观看| 亚洲av中文字字幕乱码综合 | 精品欧美一区二区三区在线| 国产高清视频在线播放一区| 日本a在线网址| 亚洲国产精品合色在线| 久久中文字幕人妻熟女| 亚洲午夜理论影院| 美女高潮到喷水免费观看| 成人一区二区视频在线观看| 国产精品一区二区三区四区久久 | 侵犯人妻中文字幕一二三四区| 老鸭窝网址在线观看| 亚洲国产精品合色在线| 国产伦人伦偷精品视频| 一本综合久久免费| 久久中文字幕一级| 99久久99久久久精品蜜桃| a级毛片a级免费在线| 真人做人爱边吃奶动态| videosex国产| 国产1区2区3区精品| 免费高清视频大片| 国内精品久久久久精免费| 亚洲av美国av| 黑人欧美特级aaaaaa片| 黄色女人牲交| 精品欧美一区二区三区在线| 2021天堂中文幕一二区在线观 | 日韩精品青青久久久久久| 日韩有码中文字幕| 亚洲国产欧美日韩在线播放| tocl精华| 丝袜在线中文字幕| 一级作爱视频免费观看| 麻豆久久精品国产亚洲av| 两性夫妻黄色片| 身体一侧抽搐| 俺也久久电影网| 亚洲中文字幕日韩| 男人的好看免费观看在线视频 | 99久久久亚洲精品蜜臀av| 久久香蕉激情| 日本免费一区二区三区高清不卡| 麻豆成人午夜福利视频| 12—13女人毛片做爰片一| 久久久国产成人免费| 老司机靠b影院| 操出白浆在线播放| 自线自在国产av| 国产视频内射| 精华霜和精华液先用哪个| 欧美成人免费av一区二区三区| 一a级毛片在线观看| 国产激情久久老熟女| 日韩欧美一区视频在线观看| 法律面前人人平等表现在哪些方面| 日韩成人在线观看一区二区三区| 可以免费在线观看a视频的电影网站| 亚洲欧美日韩无卡精品| 色综合欧美亚洲国产小说| 日韩精品免费视频一区二区三区| 日本在线视频免费播放| 精品国产美女av久久久久小说| 两个人看的免费小视频| 久久久久久久久久黄片| 午夜免费鲁丝| 亚洲精品粉嫩美女一区| 色综合婷婷激情| 91av网站免费观看| 亚洲va日本ⅴa欧美va伊人久久| 满18在线观看网站| 一级a爱视频在线免费观看| 大香蕉久久成人网| 欧美一区二区精品小视频在线| 亚洲人成伊人成综合网2020| 在线视频色国产色| 欧美中文综合在线视频| 香蕉av资源在线| 国内精品久久久久精免费| 视频区欧美日本亚洲| 熟妇人妻久久中文字幕3abv| 母亲3免费完整高清在线观看| 亚洲国产高清在线一区二区三 | 国产人伦9x9x在线观看| 国产91精品成人一区二区三区| 久久青草综合色| 免费高清在线观看日韩| 在线观看舔阴道视频| 亚洲aⅴ乱码一区二区在线播放 | 夜夜看夜夜爽夜夜摸| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品影院久久| 精品久久蜜臀av无| 久久婷婷人人爽人人干人人爱| 91老司机精品| 亚洲人成77777在线视频| 熟妇人妻久久中文字幕3abv| 男女视频在线观看网站免费 | 老熟妇乱子伦视频在线观看| 桃红色精品国产亚洲av| 91字幕亚洲| 性欧美人与动物交配| 午夜a级毛片| 欧美久久黑人一区二区|