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

    激光選區(qū)熔化TC11鈦合金的斷裂機(jī)理及異種載荷下裂紋擴(kuò)展研究

    2022-04-15 05:54:00黃鵬于海洋郭寧周慶軍唐炳濤趙偉李新生
    精密成形工程 2022年4期
    關(guān)鍵詞:相場柱狀晶鈦合金

    黃鵬,于海洋,郭寧,周慶軍,唐炳濤,趙偉,李新生

    激光選區(qū)熔化TC11鈦合金的斷裂機(jī)理及異種載荷下裂紋擴(kuò)展研究

    黃鵬1,2,于海洋1,2,郭寧1,2,周慶軍3,唐炳濤1,2,趙偉1,2,李新生4

    (1. 齊魯工業(yè)大學(xué)(山東省科學(xué)院) 機(jī)械工程學(xué)院,濟(jì)南 250353;2. 山東省機(jī)械設(shè)計(jì)研究院,濟(jì)南 250031;3. 首都航天機(jī)械有限公司,北京 100076;4. 金雷科技股份公司,濟(jì)南 271105)

    基于航天航空領(lǐng)域損傷容限性鈦合金的設(shè)計(jì)理念,研究沉積態(tài)TC11鈦合金拉伸速率相關(guān)的斷裂機(jī)理及異種載荷形式下裂紋擴(kuò)展路徑。采用電子背散射衍射(EBSD)表征由激光選區(qū)熔化(SLM)制備的沉積態(tài)TC11鈦合金的晶粒形貌,利用掃描電子顯微鏡分析不同應(yīng)變速率下TC11鈦合金的斷口形貌和缺陷分布,針對無法原位實(shí)時(shí)追蹤裂紋擴(kuò)展路徑的問題,利用Abaqus相場UEL對裂紋擴(kuò)展進(jìn)行研究。沉積態(tài)TC11鈦合金中,∥BD截面上的晶粒主要為柱狀晶并有少部分等軸晶,⊥BD截面上的晶粒全部為等軸晶,組織為由密排立方結(jié)構(gòu)(HCP)的針狀α相構(gòu)成的馬氏體,無明顯晶粒取向;瞬時(shí)載荷下,載荷為1、1.5、2 MPa時(shí),切向裂紋寬度分別為0、7.69、14 mm,角度分別為108°、92°、82°;周期載荷下,載荷為1、1.1、1.2、1.35、1.5 MPa時(shí),裂紋寬度分別為15.56、11.87、26.23、20.51、20.92 mm。隨拉伸速率的增加,斷裂機(jī)理由韌性斷裂轉(zhuǎn)化為脆性斷裂;瞬時(shí)加載時(shí),裂紋萌生并切向擴(kuò)展,進(jìn)而裂紋分叉直至完全斷裂,且載荷的增加對沉積態(tài)TC11鈦合金切向裂紋有促進(jìn)作用,同時(shí)可抑制裂紋的分叉;周期加載時(shí),載荷的增加對裂紋寬度先是起促進(jìn)作用,裂紋達(dá)到最大值后穩(wěn)定在25~20 mm。

    激光選區(qū)熔化;TC11鈦合金;相場斷裂;裂紋擴(kuò)展;斷口形貌

    雙相鈦合金具有密度低、比強(qiáng)度高、耐腐蝕、耐高溫、耐低溫等優(yōu)異性能,廣泛應(yīng)用于航天航空領(lǐng)域的主承力構(gòu)件。與傳統(tǒng)的制造技術(shù)相比,SLM具有近凈成形、材料利用率高、可以生產(chǎn)復(fù)雜幾何形狀的變截面結(jié)構(gòu)件等巨大優(yōu)勢[1]。SLM工藝存在復(fù)雜的熱循環(huán)過程,在SLM成形的鈦合金內(nèi)部會存在孔隙等缺陷,而這些缺陷周圍常常伴隨著應(yīng)力集中,這會引起裂紋的萌生及擴(kuò)展?;谀壳昂教旌娇疹I(lǐng)域損傷容限型鈦合金設(shè)計(jì)理念,由裂紋誘發(fā)的材料失效及裂紋的擴(kuò)展演化是制約鈦合金在航天航空領(lǐng)域主承力結(jié)構(gòu)件上應(yīng)用的內(nèi)在因素。

    許多學(xué)者針對SLM成形、LMD成形的TC11鈦合金斷裂機(jī)理及組織與性能進(jìn)行了相關(guān)研究,針對退火溫度和保溫時(shí)間對SLM成形的TC11鈦合金斷裂機(jī)理的影響,竇恩惠等[2]研究了不同熱處理制度下的沉積態(tài)TC11鈦合金斷裂機(jī)制,為通過熱處理的手段改善SLM成形的鈦合金的組織性能提供了參考,并且為鈦合金的工程應(yīng)用提供了一定的參考。針對復(fù)雜結(jié)構(gòu)件實(shí)際成形熱過程中SLM成形鈦合金的斷裂機(jī)理,周慶軍等[3]研究了沉積態(tài)TC11鈦合金的組織及力學(xué)性能,為調(diào)控鈦合金激光熔化沉積(LMD)成形的大型復(fù)雜構(gòu)件的組織性能、抑制成形過程中出現(xiàn)的開裂提供了重要依據(jù)。針對不同溫度和時(shí)間下熱暴露對TC11鈦合金斷裂機(jī)理的影響,張衛(wèi)方等[4]研究了使用溫度附近TC11鈦合金在不同溫度和時(shí)間下熱暴露的斷裂機(jī)理,揭示了熱暴露對TC11鈦合金斷裂機(jī)理的影響。由于通過實(shí)驗(yàn)的手段無法原位實(shí)時(shí)追蹤裂紋擴(kuò)展路徑,對沉積態(tài)TC11鈦合金裂紋誘發(fā)的斷裂機(jī)理尚缺乏深入的研究,因此,利用數(shù)值模擬來實(shí)現(xiàn)實(shí)時(shí)追蹤裂紋擴(kuò)展的路徑尤為關(guān)鍵。

    許多學(xué)者提出了多種模擬固體材料裂紋擴(kuò)展的方法,如Guo等[5]提出了一種能夠模擬任意形狀裂紋表面的二維裂紋擴(kuò)展方法,這種方法與XFEM相比,可直接用網(wǎng)格描述裂紋,不需要引入額外的技術(shù)來追蹤界面,但其很難擴(kuò)展應(yīng)用到復(fù)雜的三維問題上。Zhou等[6]開發(fā)了一種動態(tài)嵌入內(nèi)聚元素的三維有限元模型處理網(wǎng)格依賴關(guān)系,但內(nèi)聚單元的裂紋一定是沿每個(gè)單元邊界傳播的,所以必然會有一定程度的網(wǎng)格依賴無法消除。Azevedo等[7]提出了一種混合方法,將粒子法應(yīng)用于較大結(jié)構(gòu)的分析,即僅在斷裂過程中使用剛性的圓形顆粒,而其余部分使用有限元建模。Gergely等[8-9]提出了一種新的相場UEL來實(shí)現(xiàn)對裂紋擴(kuò)展路徑的追蹤,這種方法是基于Belytschko團(tuán)隊(duì)[10-11]提出的XFEM方法以及Gürses等[12]提出的配置力驅(qū)動的鋒利裂縫前緣思想構(gòu)建的,這種方法不需要預(yù)先定義裂縫,可以模擬出曲折路線的裂紋路徑、裂紋分叉及合并,但由于采用了交錯時(shí)的集成算法,所以需要一個(gè)相對較小的時(shí)間步長,故無法模擬長時(shí)間載荷狀態(tài)下的裂紋擴(kuò)展路徑。

    文中通過EBSD表征了沉積態(tài)TC11鈦合金晶粒和組織形貌及取向信息,通過SEM分析了不同應(yīng)變速率下TC11鈦合金拉伸試棒斷口形貌并揭示了其斷裂機(jī)理,最后通過Abaqus相場UEL預(yù)測了瞬時(shí)與周期加載的條件下沉積態(tài)TC11鈦合金裂紋的擴(kuò)展路徑,揭示了不同類載荷下裂紋擴(kuò)展的演化規(guī)律,為航天航空領(lǐng)域損傷容限型鈦合金的設(shè)計(jì)提供了一定的參考,并為TC11鈦合金工程應(yīng)用提供了依據(jù)。

    1 TC11鈦合金組織表征及力學(xué)性能實(shí)驗(yàn)方法

    1.1 SLM工藝

    采用湖南華曙高科技股份有限公司生產(chǎn)的FS271型激光選區(qū)熔化成形設(shè)備,原材料采用西安歐中材料科技有限公司生產(chǎn)的粒徑為15~53 μm的TC11粉末,其具體化學(xué)成分見表1。SLM具體工藝參數(shù)如下:激光功率為400 W,掃描速度為1 400 mm/s,掃描間距為0.12 mm,分層厚度為60 μm。

    表1 TC11粉末成分(質(zhì)量分?jǐn)?shù))

    1.2 不同應(yīng)變速率的TC11鈦合金拉伸實(shí)驗(yàn)方法

    參照GB/T 228.1—2010[13],利用激光加工獲得的拉伸試棒具體尺寸如圖1所示。采用濟(jì)南鑫光實(shí)驗(yàn)機(jī)制造有限公司生產(chǎn)的WDW-200E電子萬能實(shí)驗(yàn)機(jī)測試沉積態(tài)TC11鈦合金室溫拉伸性能,拉伸速率分別為0.1、1、2、5 mm/min,最大拉伸力為27~28 kN,拉伸方向垂直于SLM的沉積方向。

    1.3 TC11鈦合金斷口形貌表征方法

    采用北京日立高新科學(xué)儀器有限公司生產(chǎn)的SU3500掃描電子顯微鏡對拉斷件的斷口形貌進(jìn)行觀察。實(shí)驗(yàn)采用15 000 V的加速電壓、148 000 mA左右的發(fā)射電流。

    1.4 TC11鈦合金晶粒形貌及組織表征方法

    采用EBSD來表征沉積態(tài)TC11鈦合金晶粒形貌及組織。EBSD樣品制備過程如下:①利用金相砂紙打磨至800目,保證無劃痕或者有極少量的淺劃痕;②使用粒度為1~3 μm的氧化鋁拋光劑粗拋,使用粒度為0.2~0.5 μm的氧化鋁拋光劑細(xì)拋;③電解拋光采用體積比為1︰6︰3的高氯酸、甲醇、正丁醇的混合溶液,陰極材料為不銹鋼,電壓為30 V,電流密度為100~ 150 A/dm2,在溫度≤?10 ℃的條件下電解30~60 s。EBSD采用日本電子JSM-7900F掃描電鏡,采用edax Hikari XP探頭,最后采用EDAX OIM軟件處理數(shù)據(jù)。掃描區(qū)域?yàn)?1 μm× 48 μm,步長為0.1 μm,放大倍數(shù)為1 000。

    圖1 TC11鈦合金拉伸試樣尺寸

    2 結(jié)果與分析

    2.1 沉積態(tài)TC11鈦合金晶粒及組織表征

    將沉積態(tài)TC11鈦合金(拉伸方向⊥BD)長方體試棒按如圖2所示的方向切開,在其截面的9個(gè)區(qū)域內(nèi),分別通過EBSD表征其晶粒取向及組織形貌。

    圖2 原始態(tài)樣品切口位置及EBSD序號分布

    通過圖3a—c可以觀察到EBSD晶粒取向,圖2中1—3區(qū)域內(nèi)主要為柱狀晶,在某些區(qū)域內(nèi)會出現(xiàn)等軸晶與柱狀晶交替生長的情況,如圖3c所示。比較圖3a—c可以得出,1區(qū)域的柱狀晶最為粗大,而2區(qū)域的柱狀晶尺寸較小。柱狀晶內(nèi)部主要為由α相構(gòu)成的細(xì)針狀馬氏體組織,如圖3a所示,同時(shí)通過立方鈦與六方鈦的分布(見圖3e)可知,組織內(nèi)部多為由HCP結(jié)構(gòu)的六方鈦構(gòu)成的α相,同時(shí)由于SLM工藝熔池冷卻速度較快,還存在極少量由BCC結(jié)構(gòu)的立方鈦構(gòu)成的β相,六方鈦體積分?jǐn)?shù)為99.4%,而立方鈦的體積分?jǐn)?shù)僅為0.6%,如圖3f所示。部分區(qū)域出現(xiàn)了粗大的針狀馬氏體組織,如圖3b所示。這種細(xì)針狀的馬氏體組織對提高沉積態(tài)TC11鈦合金的綜合力學(xué)性能有益。從圖3可以觀察到,SLM制成的TC11鈦合金無明顯的晶粒取向,各晶粒之間取向不同。

    圖2中6—8區(qū)域的EBSD晶粒取向分別見圖4a—c,可以觀察到在6—8區(qū)域內(nèi),晶粒為等軸晶晶粒。晶粒內(nèi)部仍為由細(xì)針狀α相組成的馬氏體,同時(shí)與1—3區(qū)域相同,組織內(nèi)六方鈦占絕大多數(shù),體積分?jǐn)?shù)為99.8%,而立方鈦占比很小,體積分?jǐn)?shù)僅為0.2%,見圖4d。

    ∥BD方向截面與⊥BD方向截面的微觀組織完全一致,均為由六方鈦組成的α相組成的馬氏體組織,但2個(gè)區(qū)域晶粒形狀有明顯差異,垂直于熔道面的1—3區(qū)域多為長條狀柱狀晶,而在熔道面上的6—8區(qū)域則基本為等軸晶。

    2.2 斷口形貌及機(jī)理分析

    不同拉伸速率下SLM成形斷口形貌見圖5,可以看到隨著應(yīng)變速率的增加,邊緣延展區(qū)(塑性區(qū))的面積逐漸減小,而中間瞬斷區(qū)(脆性區(qū))的面積逐漸增大,在拉伸速率為5 mm/min時(shí)邊緣延展區(qū)面積非常小,斷裂面絕大部分為瞬斷區(qū)。從宏觀來看,隨著拉伸速率的增加,縮頸現(xiàn)象越來越不明顯。從高倍SEM中可以看到斷面上中心瞬斷區(qū)的孔隙,同時(shí)在部分孔隙的邊緣觀察到細(xì)小的微裂紋,這是由孔隙周圍的應(yīng)力集中造成的,如圖5f和g所示,并且隨著拉伸速率的增加,斷面的韌窩數(shù)量減少,韌窩尺寸變小,說明隨著拉伸速率的增加,TC11鈦合金的斷裂由韌性斷裂逐漸向脆性斷裂轉(zhuǎn)化。

    圖4 沉積態(tài)TC11鈦合金晶粒及組織形貌

    圖5 不同拉伸速率下SLM成形斷口形貌

    3 裂紋擴(kuò)展的相場斷裂模擬

    3.1 斷裂表面的相場描述

    相場法主要是通過建立內(nèi)部長度尺度c與損傷變量的關(guān)系研究材料的斷裂問題,將材料內(nèi)部的尖銳裂紋擴(kuò)散到擁有彈性或塑性的固體體積中,通過監(jiān)測損傷變量來判斷局部的材料是否發(fā)生斷裂,若=1表示材料發(fā)生斷裂,若<1表示材料未發(fā)生斷裂。損傷變量與內(nèi)部長度尺度c的關(guān)系見式(1)。

    (1)

    斷裂面多維度裂紋表面密度與裂紋長度參數(shù)相關(guān),由相場函數(shù)引入裂紋表面密度,裂紋表面密度可以由式(2)表示[14]:

    (2)

    理論斷裂面是由式(2)的整體體積積分獲得的:

    3.2 能量函數(shù)

    彈塑性動態(tài)問題的能量函數(shù)包括式(4)的拉格朗日函數(shù):

    由式(8)可以看出材料在拉伸時(shí)會發(fā)生退化而在壓縮時(shí)材料不發(fā)生變化[15]。退化函數(shù)()為:

    式中:表示退化函數(shù)整體的數(shù)值穩(wěn)定性,為10?8級。

    塑性應(yīng)變能為:

    相場法的最關(guān)鍵思想就是將本身不連續(xù)的點(diǎn)作為已發(fā)生斷裂與未發(fā)生斷裂材料之間的一個(gè)連續(xù)場,故斷裂能量可由式(2)改寫為損傷變量的一個(gè)函數(shù),即:

    式中:表示斷裂能的能量密度;c表示產(chǎn)生一個(gè)單位的斷裂表面所需要的表面能;c表示長度尺度參數(shù),主要表征了損傷擴(kuò)散的尺度。

    3.3 Abaqus相場UEL

    Abaqus相場UEL的實(shí)現(xiàn)過程如下:以分層的方式使用2種元素類型,同時(shí)這2層中每一層都與同一個(gè)節(jié)點(diǎn)相連接,但不同的層之間對同一節(jié)點(diǎn)不同自由度(DOF)的剛度有著不同的貢獻(xiàn)。圖6a為用二維表示的三維有限元結(jié)構(gòu)在Abaqus中的框架圖,具體操作流程如圖6b所示。

    圖6 Abaqus相場UEL框架圖及實(shí)現(xiàn)方式

    3.4 裂紋擴(kuò)展的相場斷裂模擬

    3.4.1 有限元模型建立

    對沉積態(tài)TC11鈦合金進(jìn)行裂紋擴(kuò)展的相場模擬。在厚度為1 mm、長度為100 mm、寬度為40 mm的二維矩形板的上下端分別施加瞬時(shí)載荷與周期載荷,瞬時(shí)載荷分別為1、1.5、2 MPa,周期載荷分別為1.5、1.2、1.1、1 MPa,周期加載方式為正弦加載,如圖7a所示。網(wǎng)格單元類型為CPS4R,單元數(shù)量為1 000,如圖7b所示。材料屬性如下:彈性模量為115 GPa,泊松比為0.33,密度為4 600 kg/m3,斷裂能量密度c為3 J/m3,內(nèi)部長度尺度c為0.004 m,不施加硬化和應(yīng)變懲罰,即硬化常數(shù)=0、應(yīng)變懲罰ε=0。瞬時(shí)載荷中分析步長為1×10?7s,總分析時(shí)間為0.000 1 s。周期載荷中分析步長為0.000 2 s,總分析時(shí)間為1 s。

    圖7 數(shù)值模擬相關(guān)參數(shù)

    3.4.2 不同加載方式下的TC11鈦合金裂紋擴(kuò)展演化研究

    瞬時(shí)載荷相場UEL模擬結(jié)果如圖8所示??梢钥闯鲭S著施加載荷的增大,裂紋擴(kuò)展路徑也隨之發(fā)生變化。當(dāng)載荷為1MPa時(shí),切向裂紋的長度為0,在預(yù)制裂紋的前端就發(fā)生裂紋分叉,由于載荷較小其裂紋沒有擴(kuò)展到材料邊緣發(fā)生斷裂。當(dāng)載荷為1.5 MPa時(shí),切向裂紋長度明顯增加,當(dāng)切向裂紋擴(kuò)展一段距離后,裂紋開始發(fā)生分叉現(xiàn)象,此時(shí)切向裂紋的長度為7.69 mm,同時(shí)裂紋即將擴(kuò)展到材料邊緣發(fā)生斷裂。當(dāng)載荷增加到2 MPa時(shí),切向裂紋長度進(jìn)一步增加到14 mm,此時(shí)大載荷使材料發(fā)生斷裂。同時(shí),發(fā)現(xiàn)裂紋分叉以切向裂紋為軸呈軸對稱分布,隨著載荷的增加,裂紋分叉的角度逐漸減小,在加載載荷為1 MPa時(shí),裂紋分叉角度最大約為108°,在載荷為2 MPa時(shí),裂紋分叉角度最小約為82°。

    圖8 瞬時(shí)載荷下裂紋擴(kuò)展路徑云圖

    施加周期載荷的相場UEL結(jié)果如圖9所示。由于施加周期載荷的時(shí)間足夠長,所以均發(fā)生斷裂,并且均無切向裂紋產(chǎn)生,說明在周期載荷加載的條件下,TC11鈦合金很難產(chǎn)生切向裂紋,主要裂紋方向?yàn)檠乩旆较虻牧鸭y。其次,分別在沿拉伸方向上取上、中、下3個(gè)位置測量其裂紋寬度,經(jīng)計(jì)算得到均值進(jìn)行比較,可以發(fā)現(xiàn)隨著周期載荷的增加,裂紋的切向?qū)挾戎饾u增加。在載荷為1.2 MPa時(shí),裂紋寬度最大為26.23 mm,隨后逐漸降低,穩(wěn)定在20 mm左右,在1、1.1、1.2、1.35、1.5 MPa的載荷下,裂紋寬度分別為15.56、11.87、26.23、20.51、20.92 mm,說明在載荷為1~1.2 MPa時(shí),載荷的增加會促進(jìn)裂紋的擴(kuò)展,在1.2 MPa時(shí)促進(jìn)效果最明顯,在1.2~1.5 MPa時(shí)載荷的增加對裂紋寬度的影響較小,整體保持在25~20 mm。

    圖9 周期載荷下裂紋擴(kuò)展路徑云圖

    4 結(jié)論

    通過EBSD、SEM和Abaqus相場UEL研究了SLM制備的沉積態(tài)TC11鈦合金的組織、斷口形貌和裂紋擴(kuò)展情況,揭示了沉積態(tài)TC11鈦合金應(yīng)變速率相關(guān)的斷裂機(jī)理,預(yù)測了瞬時(shí)與周期加載下不同加載載荷的裂紋擴(kuò)展演化,得到以下結(jié)論。

    1)沉積態(tài)TC11鈦合金∥BD截面方向上,微觀組織多為長條狀柱狀晶,部分區(qū)域出現(xiàn)柱狀晶與等軸晶交替生長的情況,晶粒內(nèi)部大部分為由HCP結(jié)構(gòu)六方鈦構(gòu)成的針狀α相和極少量的β相?!虰D截面方向上的微觀組織與∥BD截面方向上的組織一致,均為由α相構(gòu)成的馬氏體組織,但晶粒形狀為等軸晶,無明顯晶粒取向。

    2)隨著拉伸速率的增加,沉積態(tài)TC11鈦合金斷面的韌窩數(shù)量減少,韌窩尺寸變小,邊緣延展區(qū)(韌性區(qū))的面積減小,斷裂由韌性斷裂逐漸向脆性斷裂轉(zhuǎn)化。

    3)對沉積態(tài)TC11鈦合金施加瞬時(shí)載荷時(shí),裂紋出現(xiàn)分叉并以切向裂紋為軸呈軸對稱分布,隨施加載荷的增加,裂紋分叉的角度逐漸減小,隨施加載荷的增加,切向裂紋長度逐漸增加,在載荷為1 MPa和1.5 MPa時(shí),材料未發(fā)生斷裂,在載荷為2 MPa時(shí)材料發(fā)生斷裂。沉積態(tài)TC11鈦合金施加周期載荷時(shí),裂紋無明顯切向的裂紋擴(kuò)展,在載荷為1~1.2 MPa時(shí)隨加載載荷的增加裂紋寬度增加,1.2 MPa時(shí)裂紋寬度最大,在載荷為1.2~1.5 MPa時(shí),隨載荷的增加,裂紋寬度穩(wěn)定在25~20 mm。

    [1] LORE T, FREDERIK V, TOM C, et al. A Study of the Microstructural Evolution during Selective Laser Melting of Ti-6Al-4V[J]. Acta Materialia, 2010, 58: 3303-3312.

    [2] 竇恩惠, 肖美立, 柯林達(dá), 等. 熱處理對激光選區(qū)熔化成形TC11鈦合金組織性能的影響[J]. 中國激光, 2021(6): 201-209.

    DOU En-hui, XIAO Mei-li, KE Lin-da, et al. Effect of Heat Treatment on Microstructure and Mechanical Properties of Selective-Laser-Melted TC11 Titanium Alloys[J]. Chinese Journal of Lasers, 2021(6): 201-209.

    [3] 周慶軍, 嚴(yán)振宇, 韓旭, 等. 激光熔化沉積TC11鈦合金的組織與力學(xué)性能[J]. 中國激光, 2018, 45(11): 1-6.

    ZHOU Qing-jun, YAN Zhen-yu, HAN Xu, et al. Microstructure and Mechanical Properties of Laser Melting Deposited TC11 Titanium Alloys[J]. Chinese Journal of Lasers, 2018, 45(11): 1-6.

    [4] 張衛(wèi)方, 高威, 劉高遠(yuǎn), 等. 熱暴露對TC11鈦合金斷裂行為的影響[J]. 航空材料學(xué)報(bào), 2004, 24(2): 1-5.

    ZHANG Wei-fang, GAO Wei, LIU Gao-yuan, et al. The Influence of High Temperature Exposure on Fracture Behaviors of TC11 Titanium Alloy[J]. Journal of Aeronautical Materials, 2004, 24(2): 1-5.

    [5] GUO Liang-peng, WANG Yu-heng. A Node Split Method for Crack Growth Problem[J]. Applied Mechanics and Materials. 2012, 182/183: 1524-1528.

    [6] ZHOU F, MOLINARI J F. Dynamic Crack Propagation with Cohesive Elements: a Methodology to Address Mesh Dependency[J]. International Journal for Numerical Meyhods in Engineering, 2004, 59(1): 1-24.

    [7] AZEVEDO N M, LEMOS J V. Hybrid Discrete Element/Finite Element Method for Fracture Analysis[J]. Computer Methods in Applied Mechanics and Engineering, 2006, 195(33/34/35/36): 4579-4593.

    [8] GERGELY M, ANTHONY G. 2D and 3D Abaqus Implementation of a Robust Staggered Phase-Field Solution for Modeling Brittle Fracture[J]. Finite Elements in Analysis and Design, 2017, 130: 27-38.

    [9] GERGELY M, ANTHONY G, RIAN S, et al. An Open-Source Abaqus Implementation of the Phase-Field Method to Study the Effect of Plasticity on the Instantaneous Fracture Toughnessin Dynamic Crack Propagation[J]. Computer Methods in Applied Mechanics and Engineering, 2020, 365: 1-27.

    [10] MO?S N, DOLBOW J, BELYTSCHKO T. Afinite Element Method for Crack Growth without Remeshing[J]. International Journal for Numerical Meyhods in Engineering, 1999, 46(1): 131-150.

    [11] MO?S N, GRAVOUIL A, BELYTSCHKO T. Non-Planar 3D Crack Growth by the Extended Finite Element and Level Sets-Part I: Mechanical Model[J]. International Journal for Numerical Meyhods in Engineering, 2002, 53(11):2549-2568

    [12] GüRSES E, MIEHE C. A Computational Framework of Three-Dimensional Configurational-Force-Driven Brittle Crack Propagation[J]. Computer Methods in Applied Mechanics and Engineering, 2009, 198(15/16): 1413- 1428.

    [13] GB/T 228.1—2010, 金屬材料拉伸試驗(yàn)第1部分: 室溫試驗(yàn)方法[S].

    GB/T 228.1—2010, Tensile Testing of Metallic Materials Part 1: Room Temperature Test Method[S].

    [14] MIEHE C, WELSCHINGER F, HOFACKKER M. Thermodynamically Consistent Phase-Field Models of Fracture: Variational Principles and Multi-Field FE Implementations[J]. International Journal for Numerical Meyhods in Engineering, 2010, 83(10): 1273-1311.

    [15] BERNARD P E, MO?S N, CHEVAUGEON N. Damage Growth Modeling Using the Thick Level Set (TLS) Approach: Efficient Discretization for Quasi-Static Loadings, Comput[J]. Computer Methods in Applied Mechanics and Engineering, 2012, 233: 11-27.

    [16] 周英豪. 選區(qū)激光熔化增材制造Ti-22Al-25Nb高溫合金的組織與性能研究[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2020: 47-55.

    ZHOU Ying-hao. Study on Microstructure and Properties of Ti-22Al-25Nb Superalloy Produced by Selective Laser Melting Additive[D]. Harbin: Harbin Institute of Technology, 2020: 47-55.

    [17] 尹衍軍. 選區(qū)激光熔化成型316L不銹鋼流動規(guī)律及組織、性能研究[D]. 北京: 北京科技大學(xué), 2019: 49-119.

    YIN Yan-jun. Study on Flow Law, Microstructure and Properties of 316L Stainless Steel by Selective Laser Melting[D]. Beijing: Beijing University of Science and Technology, 2019: 49-119.

    [18] 王震. 激光選區(qū)熔化Ti6Al4V合金及其表面超聲滾壓加工的組織與性能研究[D]. 廣州: 華南理工大學(xué), 2019: 39-41.

    WANG Zhen. Laser Selective Melting Ti6Al4V Study on Microstructure and Properties of Alloy and Its Surface by Ultrasonic Rolling[D]. Guangzhou: South China University of Technology, 2019: 39-41.

    Fracture Mechanism and Crack Propagation of TC11 Titanium Alloy Melted by Laser Selective Melting under Different Loads

    HUANG Peng1,2, YU Hai-yang1,2, GUO Ning1,2, ZHOU Qing-jun3, TANG Bing-tao1,2, ZHAO Wei1,2, LI Xin-sheng4

    (1. School of Mechanical Engineering, Qilu University of Technology (Shandong Academy of Sciences), Jinan 250353, China; 2. Shandong Mechanical Design and Research Institute, Jinan 250031, China; 3. Capital Aerospace Machinery Corporation Limited, Beijing 100076, China; 4. Jinlei Technology Co., Ltd., Jinan 271105, China)

    The work aims to study the fracture law and crack propagation path related to the rate of deposited TC11 titanium alloy according to the current design concept of damage tolerant titanium alloy in aerospace field. In this paper, the grain morphology of deposited TC11 titanium alloy prepared by laser selective melting (SLM) was characterized by electron backscattered diffraction (EBSD). The fracture morphology and defect distribution of TC11 titanium alloy under different strain rates were analyzed with scanning electron microscope (SEM). Aiming at the problem that the crack propagation path cannot be tracked in situ and in real time, the Abaqus phase field UEL was used to study the crack propagation. In the deposited TC11 titanium alloy, the grains on the//BD section were mainly columnar crystals and a small part of equiaxed crystals. All the grains on the^BD section were equiaxed crystals and the structure was needle like with dense cubic structure (HCP) α phase martensite tissue. No obvious grain orientation was observed; Under instantaneous load, when the load was 1, 1.5, 2 MPa, the tangential crack width was 0, 7.69, 14 mm respectively. The angles were 108°, 92°, 82° respectively; Under cyclic load, when the load was 1, 1.1, 1.2, 1.35, 1.5 MPa, the crack width was 15.56, 11.87, 26.23, 20.51, 20.92 mm respectively. With the increase of tensile rate, the fracture mechanism changes from ductile fracture to brittle fracture; Under transient loading, the crack initiates and propagates in a tangential sway, and then the crack bifurcates until complete fractures. The increase of load can promote the tangential crack of deposited TC11 titanium alloy and inhibit the crack bifurcation at the same time; During periodic loading, the increase of load promotes the crack width, and then the crack stabilizes at 25~20 mm after it reaches the maximum value.

    selective laser melting (SLM); TC11 titanium alloy; phase field fracture; crack propagation; fracture morphology

    10.3969/j.issn.1674-6457.2022.04.011

    TG146.2+3

    A

    1674-6457(2022)04-0094-08

    2021-12-14

    國家自然科學(xué)基金(52105370);山東省自然科學(xué)基金(ZR2020QE168)

    黃鵬(1998—),男,碩士生,主要研究方向?yàn)榧す庠霾闹圃焘伜辖鸬牧W(xué)性能。

    郭寧(1989—),男,博士,講師,主要研究方向?yàn)榧す庠霾闹圃旃に?組織-力學(xué)性能模擬。

    責(zé)任編輯:蔣紅晨

    猜你喜歡
    相場柱狀晶鈦合金
    基于子單元光滑有限元的混凝土相場損傷模型研究
    鑄件凝固微觀組織仿真程序開發(fā)
    “神的金屬”鈦合金SHINE YOUR LIFE
    中國自行車(2018年8期)2018-09-26 06:53:32
    基于相場理論的瀝青自愈合微觀進(jìn)程與機(jī)理研究進(jìn)展
    石油瀝青(2018年1期)2018-04-12 07:31:51
    鈦合金板鍛造的工藝實(shí)踐
    四川冶金(2017年6期)2017-09-21 00:52:30
    基于COMSOL的相場模擬研究
    科技視界(2017年8期)2017-07-31 10:31:17
    對流擴(kuò)散-多相相變體系內(nèi)柱狀晶/等軸晶形成過程的數(shù)值模擬
    材料工程(2017年6期)2017-06-22 13:33:24
    PtCo23.2永磁合金加工開裂原因分析
    消防射水對導(dǎo)線二次短路熔痕的影響分析
    醫(yī)用鈦合金的研究與應(yīng)用
    国产视频内射| 母亲3免费完整高清在线观看| 97人妻精品一区二区三区麻豆| 人人妻,人人澡人人爽秒播| 我的老师免费观看完整版| 12—13女人毛片做爰片一| 精品国产亚洲在线| 亚洲中文字幕一区二区三区有码在线看 | 亚洲中文av在线| 欧美日韩乱码在线| 男人的好看免费观看在线视频 | 亚洲成人久久性| 亚洲美女视频黄频| 久久精品影院6| 琪琪午夜伦伦电影理论片6080| 99国产综合亚洲精品| 夜夜爽天天搞| 国产高清videossex| 看片在线看免费视频| 亚洲国产精品合色在线| 一二三四在线观看免费中文在| 小说图片视频综合网站| 亚洲 国产 在线| 亚洲av成人不卡在线观看播放网| 中文字幕人成人乱码亚洲影| 少妇被粗大的猛进出69影院| 一区福利在线观看| 黄色女人牲交| 在线观看一区二区三区| videosex国产| 国产精品亚洲av一区麻豆| 国产亚洲精品av在线| 日本免费一区二区三区高清不卡| 两个人免费观看高清视频| 一夜夜www| 最新在线观看一区二区三区| 岛国在线观看网站| 搡老岳熟女国产| 免费观看人在逋| av片东京热男人的天堂| 亚洲国产欧美网| 欧美三级亚洲精品| 人人妻人人看人人澡| 亚洲成人久久爱视频| 很黄的视频免费| 人妻久久中文字幕网| 中文字幕最新亚洲高清| 午夜久久久久精精品| netflix在线观看网站| 亚洲av片天天在线观看| 色老头精品视频在线观看| 熟女电影av网| 无遮挡黄片免费观看| 国产亚洲精品第一综合不卡| 国产亚洲精品第一综合不卡| 国产91精品成人一区二区三区| 90打野战视频偷拍视频| 久久精品成人免费网站| 国产v大片淫在线免费观看| 亚洲专区字幕在线| 亚洲片人在线观看| 日韩大码丰满熟妇| 亚洲人成网站高清观看| 999久久久国产精品视频| 色综合站精品国产| 免费看日本二区| 不卡av一区二区三区| 亚洲国产精品999在线| 欧美日韩乱码在线| 正在播放国产对白刺激| 我的老师免费观看完整版| 国产av不卡久久| 禁无遮挡网站| 美女扒开内裤让男人捅视频| 国产高清videossex| 亚洲人成网站高清观看| 黄色丝袜av网址大全| 丁香六月欧美| 国产午夜精品论理片| 夜夜看夜夜爽夜夜摸| 宅男免费午夜| 亚洲一卡2卡3卡4卡5卡精品中文| 精品高清国产在线一区| 精华霜和精华液先用哪个| av视频在线观看入口| 国产区一区二久久| 日日摸夜夜添夜夜添小说| 精品高清国产在线一区| 国产精品久久视频播放| 亚洲精品久久国产高清桃花| 久久久国产成人免费| 91麻豆av在线| 中文字幕高清在线视频| 小说图片视频综合网站| 美女大奶头视频| 两个人视频免费观看高清| 欧洲精品卡2卡3卡4卡5卡区| 悠悠久久av| 90打野战视频偷拍视频| 蜜桃久久精品国产亚洲av| 一级黄色大片毛片| 国产91精品成人一区二区三区| 一本一本综合久久| 国产成+人综合+亚洲专区| 97人妻精品一区二区三区麻豆| 人成视频在线观看免费观看| 99热这里只有精品一区 | 99在线人妻在线中文字幕| 夜夜爽天天搞| 男插女下体视频免费在线播放| 欧美色欧美亚洲另类二区| 人妻久久中文字幕网| 两个人看的免费小视频| 9191精品国产免费久久| 久久久久免费精品人妻一区二区| 国产乱人伦免费视频| 午夜福利在线在线| 18禁国产床啪视频网站| 国产伦人伦偷精品视频| 午夜精品在线福利| 9191精品国产免费久久| 他把我摸到了高潮在线观看| 午夜激情福利司机影院| 又黄又粗又硬又大视频| 亚洲18禁久久av| 国产成人一区二区三区免费视频网站| 日韩国内少妇激情av| 一区福利在线观看| 在线观看日韩欧美| 欧美精品啪啪一区二区三区| АⅤ资源中文在线天堂| 精品久久蜜臀av无| 12—13女人毛片做爰片一| 99国产综合亚洲精品| 亚洲av电影不卡..在线观看| 日本黄色视频三级网站网址| 叶爱在线成人免费视频播放| 成人三级做爰电影| 不卡一级毛片| 精品免费久久久久久久清纯| 91大片在线观看| 成人av一区二区三区在线看| 亚洲七黄色美女视频| 亚洲第一电影网av| 亚洲真实伦在线观看| 日韩av在线大香蕉| 制服人妻中文乱码| 亚洲乱码一区二区免费版| 成人午夜高清在线视频| 亚洲精品av麻豆狂野| 亚洲国产精品999在线| 1024视频免费在线观看| 国产亚洲精品综合一区在线观看 | 国内久久婷婷六月综合欲色啪| 国产一区二区三区视频了| 亚洲狠狠婷婷综合久久图片| 伦理电影免费视频| av中文乱码字幕在线| 亚洲国产精品成人综合色| 少妇熟女aⅴ在线视频| 国产亚洲精品第一综合不卡| 久久午夜亚洲精品久久| 一本大道久久a久久精品| 一区福利在线观看| 亚洲人成伊人成综合网2020| 国产精品久久久久久久电影 | 伊人久久大香线蕉亚洲五| 亚洲激情在线av| 99国产综合亚洲精品| 99精品久久久久人妻精品| 国产人伦9x9x在线观看| 激情在线观看视频在线高清| 悠悠久久av| 精品国产美女av久久久久小说| 亚洲七黄色美女视频| av片东京热男人的天堂| 91老司机精品| 午夜福利18| 中亚洲国语对白在线视频| 国产av麻豆久久久久久久| 欧美av亚洲av综合av国产av| 好男人在线观看高清免费视频| 一本大道久久a久久精品| 91麻豆精品激情在线观看国产| 免费在线观看黄色视频的| 日本熟妇午夜| 麻豆国产97在线/欧美 | 村上凉子中文字幕在线| 国产亚洲av嫩草精品影院| 国产熟女xx| 国产又色又爽无遮挡免费看| 精品一区二区三区四区五区乱码| 全区人妻精品视频| 亚洲精品中文字幕一二三四区| 欧美大码av| 欧美+亚洲+日韩+国产| cao死你这个sao货| 午夜日韩欧美国产| 两个人免费观看高清视频| 国产成人aa在线观看| 国产精品一区二区三区四区免费观看 | 777久久人妻少妇嫩草av网站| 日韩三级视频一区二区三区| 午夜免费激情av| 99热这里只有精品一区 | 在线看三级毛片| 中文字幕久久专区| 亚洲欧美精品综合一区二区三区| 亚洲精品在线观看二区| 久久久久久久久久黄片| 一夜夜www| 桃色一区二区三区在线观看| 最近最新中文字幕大全电影3| 国产乱人伦免费视频| 欧美在线黄色| 久久久久久久精品吃奶| 人妻久久中文字幕网| 国产激情欧美一区二区| 三级男女做爰猛烈吃奶摸视频| 波多野结衣高清无吗| 19禁男女啪啪无遮挡网站| 久久久久九九精品影院| 国产精品一区二区免费欧美| 国产一区在线观看成人免费| 免费在线观看成人毛片| 99久久久亚洲精品蜜臀av| 99国产精品一区二区蜜桃av| 国产伦在线观看视频一区| 五月玫瑰六月丁香| svipshipincom国产片| 国产精品影院久久| av福利片在线| 亚洲成人免费电影在线观看| 免费观看人在逋| 午夜精品一区二区三区免费看| 精华霜和精华液先用哪个| av福利片在线| 妹子高潮喷水视频| 国模一区二区三区四区视频 | 国产爱豆传媒在线观看 | 国产高清视频在线播放一区| 一本久久中文字幕| av免费在线观看网站| 听说在线观看完整版免费高清| 中文字幕最新亚洲高清| 精品久久久久久,| 麻豆国产av国片精品| 国产精品国产高清国产av| 久久久久久国产a免费观看| 成人18禁在线播放| 久久婷婷成人综合色麻豆| 久久香蕉国产精品| 亚洲人成电影免费在线| 精品国产乱子伦一区二区三区| 黄色 视频免费看| 麻豆国产97在线/欧美 | 午夜成年电影在线免费观看| 国语自产精品视频在线第100页| 久久婷婷人人爽人人干人人爱| 69av精品久久久久久| 丰满人妻熟妇乱又伦精品不卡| 在线观看舔阴道视频| 欧美国产日韩亚洲一区| 曰老女人黄片| 宅男免费午夜| 亚洲精品久久国产高清桃花| 国产人伦9x9x在线观看| 亚洲欧美精品综合一区二区三区| 99久久精品热视频| 国模一区二区三区四区视频 | 校园春色视频在线观看| 日韩欧美精品v在线| 精品欧美国产一区二区三| 欧美不卡视频在线免费观看 | 黄色成人免费大全| 少妇裸体淫交视频免费看高清 | 美女黄网站色视频| 天堂影院成人在线观看| 午夜a级毛片| 最新在线观看一区二区三区| 午夜两性在线视频| avwww免费| 国产爱豆传媒在线观看 | av福利片在线| 香蕉久久夜色| www日本黄色视频网| 欧美日韩乱码在线| 99久久综合精品五月天人人| 麻豆成人av在线观看| 国产精品一及| 熟女少妇亚洲综合色aaa.| 91九色精品人成在线观看| 亚洲真实伦在线观看| 啪啪无遮挡十八禁网站| 欧美成人免费av一区二区三区| 中文资源天堂在线| 在线观看一区二区三区| 制服诱惑二区| 好看av亚洲va欧美ⅴa在| 老汉色av国产亚洲站长工具| 国产在线精品亚洲第一网站| 免费在线观看完整版高清| 日韩欧美在线二视频| 日韩欧美精品v在线| 午夜激情福利司机影院| 此物有八面人人有两片| 国产成人一区二区三区免费视频网站| 天天躁夜夜躁狠狠躁躁| 国产麻豆成人av免费视频| 色哟哟哟哟哟哟| 亚洲精品久久国产高清桃花| 国产视频一区二区在线看| 亚洲一码二码三码区别大吗| 日韩精品中文字幕看吧| 最近最新中文字幕大全免费视频| 亚洲最大成人中文| 亚洲av熟女| 女人爽到高潮嗷嗷叫在线视频| 色综合站精品国产| 在线国产一区二区在线| 久久久久久大精品| 亚洲午夜理论影院| 久久久精品欧美日韩精品| 天堂av国产一区二区熟女人妻 | 色av中文字幕| 村上凉子中文字幕在线| 很黄的视频免费| 999久久久国产精品视频| 好男人在线观看高清免费视频| 午夜两性在线视频| 给我免费播放毛片高清在线观看| 亚洲国产高清在线一区二区三| 国模一区二区三区四区视频 | 久久99热这里只有精品18| 午夜福利成人在线免费观看| 亚洲精品久久成人aⅴ小说| 成人三级做爰电影| 色噜噜av男人的天堂激情| 女人被狂操c到高潮| 午夜免费成人在线视频| av片东京热男人的天堂| 久久国产乱子伦精品免费另类| 狂野欧美白嫩少妇大欣赏| ponron亚洲| 久久香蕉激情| 亚洲全国av大片| 欧美日韩福利视频一区二区| 精品电影一区二区在线| 18禁黄网站禁片午夜丰满| 欧美精品亚洲一区二区| or卡值多少钱| 麻豆国产97在线/欧美 | 久久久久九九精品影院| 成人精品一区二区免费| 巨乳人妻的诱惑在线观看| 又大又爽又粗| 丁香六月欧美| 欧美黑人精品巨大| 麻豆国产av国片精品| 国产精品自产拍在线观看55亚洲| www国产在线视频色| 麻豆成人av在线观看| a级毛片a级免费在线| 最新美女视频免费是黄的| 色在线成人网| 欧美日韩一级在线毛片| 午夜激情av网站| 亚洲18禁久久av| 美女扒开内裤让男人捅视频| 两个人的视频大全免费| 久久久久国内视频| 亚洲av五月六月丁香网| 国产av不卡久久| 久久精品国产99精品国产亚洲性色| 一边摸一边做爽爽视频免费| 国产精品亚洲一级av第二区| 最好的美女福利视频网| 免费搜索国产男女视频| 国内精品久久久久精免费| 亚洲av成人不卡在线观看播放网| 女生性感内裤真人,穿戴方法视频| 熟女少妇亚洲综合色aaa.| 日韩欧美 国产精品| 亚洲av熟女| 91老司机精品| 亚洲18禁久久av| xxx96com| 亚洲熟妇熟女久久| 91在线观看av| 国产成人精品久久二区二区91| 婷婷亚洲欧美| 午夜免费激情av| 不卡一级毛片| 美女 人体艺术 gogo| 婷婷精品国产亚洲av在线| 日本撒尿小便嘘嘘汇集6| 午夜福利18| 亚洲欧美精品综合久久99| av国产免费在线观看| 亚洲av电影不卡..在线观看| 欧美午夜高清在线| 国产99久久九九免费精品| 久久久久久国产a免费观看| 午夜日韩欧美国产| 狠狠狠狠99中文字幕| 在线十欧美十亚洲十日本专区| 亚洲一码二码三码区别大吗| 久久香蕉精品热| 午夜福利视频1000在线观看| 12—13女人毛片做爰片一| 无人区码免费观看不卡| 久久久国产欧美日韩av| 亚洲国产精品久久男人天堂| 亚洲精品色激情综合| 一级毛片精品| 欧美黄色片欧美黄色片| aaaaa片日本免费| 国产99白浆流出| 舔av片在线| 夜夜躁狠狠躁天天躁| 国产一区二区在线av高清观看| 亚洲中文日韩欧美视频| 国产亚洲精品久久久久5区| 亚洲午夜精品一区,二区,三区| 欧美成人一区二区免费高清观看 | 欧美又色又爽又黄视频| 特大巨黑吊av在线直播| 国产一区在线观看成人免费| 国产视频内射| 少妇粗大呻吟视频| 18禁美女被吸乳视频| 成人欧美大片| 国产精品一及| 校园春色视频在线观看| 波多野结衣高清无吗| 最近最新免费中文字幕在线| 国产熟女xx| 亚洲av美国av| 免费看日本二区| 两个人的视频大全免费| 一边摸一边抽搐一进一小说| 国产乱人伦免费视频| 一本久久中文字幕| 大型黄色视频在线免费观看| 中出人妻视频一区二区| 18禁国产床啪视频网站| 波多野结衣巨乳人妻| 日本免费a在线| 国产精品一区二区精品视频观看| 精品欧美国产一区二区三| 神马国产精品三级电影在线观看 | 亚洲,欧美精品.| 可以在线观看毛片的网站| 国产精品精品国产色婷婷| 亚洲一区中文字幕在线| 青草久久国产| av中文乱码字幕在线| 熟女少妇亚洲综合色aaa.| 黄片小视频在线播放| 精品久久久久久久人妻蜜臀av| 丝袜美腿诱惑在线| 一级黄色大片毛片| 女警被强在线播放| 欧美日本亚洲视频在线播放| 九九热线精品视视频播放| 日日摸夜夜添夜夜添小说| 免费观看人在逋| 久热爱精品视频在线9| 久久九九热精品免费| 亚洲全国av大片| 精品久久蜜臀av无| 日韩欧美三级三区| 9191精品国产免费久久| 国产精品 国内视频| 国产一级毛片七仙女欲春2| 中文字幕精品亚洲无线码一区| 亚洲国产精品999在线| 国产av一区在线观看免费| 亚洲精品美女久久av网站| 午夜日韩欧美国产| av超薄肉色丝袜交足视频| 我的老师免费观看完整版| 亚洲免费av在线视频| 久久婷婷人人爽人人干人人爱| 脱女人内裤的视频| 亚洲欧美一区二区三区黑人| 在线观看www视频免费| 国产成人一区二区三区免费视频网站| 男女下面进入的视频免费午夜| 亚洲在线自拍视频| 香蕉丝袜av| 国产亚洲av高清不卡| 18禁黄网站禁片午夜丰满| 国产真人三级小视频在线观看| 欧美久久黑人一区二区| 丝袜美腿诱惑在线| 亚洲熟女毛片儿| x7x7x7水蜜桃| 久久精品aⅴ一区二区三区四区| 精品国内亚洲2022精品成人| 国产伦人伦偷精品视频| 精品国产超薄肉色丝袜足j| aaaaa片日本免费| 国产av一区二区精品久久| 日日摸夜夜添夜夜添小说| 欧美成人性av电影在线观看| 人妻夜夜爽99麻豆av| 国产亚洲欧美98| 久久天堂一区二区三区四区| 亚洲午夜理论影院| 久久久久国内视频| 婷婷精品国产亚洲av在线| 色尼玛亚洲综合影院| 亚洲中文av在线| 曰老女人黄片| 免费在线观看黄色视频的| 国产不卡一卡二| 看免费av毛片| 狠狠狠狠99中文字幕| 不卡一级毛片| 99久久综合精品五月天人人| 亚洲av成人av| 两性午夜刺激爽爽歪歪视频在线观看 | 搡老熟女国产l中国老女人| 中亚洲国语对白在线视频| 又粗又爽又猛毛片免费看| 法律面前人人平等表现在哪些方面| 国产午夜精品久久久久久| 日韩精品中文字幕看吧| 操出白浆在线播放| 亚洲精品久久成人aⅴ小说| 国产av麻豆久久久久久久| 国产高清有码在线观看视频 | 久久天躁狠狠躁夜夜2o2o| www.自偷自拍.com| 级片在线观看| 嫩草影院精品99| 亚洲人成网站在线播放欧美日韩| 免费在线观看亚洲国产| 18禁国产床啪视频网站| 久久精品亚洲精品国产色婷小说| 亚洲成av人片免费观看| 巨乳人妻的诱惑在线观看| 精品国内亚洲2022精品成人| 国产精品久久久久久人妻精品电影| 国产爱豆传媒在线观看 | 国产三级黄色录像| 777久久人妻少妇嫩草av网站| 黄色片一级片一级黄色片| 黄色女人牲交| 欧美乱码精品一区二区三区| 国产亚洲av嫩草精品影院| 丁香欧美五月| 丁香六月欧美| 久久精品成人免费网站| 国产亚洲精品第一综合不卡| 91老司机精品| 99国产精品一区二区蜜桃av| 伊人久久大香线蕉亚洲五| 伦理电影免费视频| 我要搜黄色片| 99久久久亚洲精品蜜臀av| 国产激情久久老熟女| 午夜福利免费观看在线| 在线观看免费视频日本深夜| 操出白浆在线播放| 精品久久久久久,| 精品日产1卡2卡| 国产高清激情床上av| 成人18禁高潮啪啪吃奶动态图| 十八禁网站免费在线| 99国产精品99久久久久| 欧美性猛交╳xxx乱大交人| 日本熟妇午夜| 这个男人来自地球电影免费观看| 久热爱精品视频在线9| 在线观看免费日韩欧美大片| 黄色 视频免费看| 精品无人区乱码1区二区| 亚洲国产精品sss在线观看| 日韩欧美在线二视频| 男女床上黄色一级片免费看| 中文字幕高清在线视频| 免费观看人在逋| 一二三四在线观看免费中文在| 香蕉国产在线看| 国产精品影院久久| 欧美另类亚洲清纯唯美| 精品乱码久久久久久99久播| 熟妇人妻久久中文字幕3abv| 国产真实乱freesex| 婷婷六月久久综合丁香| 国产午夜精品论理片| 一级作爱视频免费观看| 国产精品久久久久久精品电影| 欧美乱色亚洲激情| 日韩有码中文字幕| 亚洲激情在线av| 少妇熟女aⅴ在线视频| 午夜老司机福利片| 亚洲精品色激情综合| 亚洲av五月六月丁香网| 桃色一区二区三区在线观看| АⅤ资源中文在线天堂| 日本成人三级电影网站| 欧美中文日本在线观看视频| 亚洲国产欧美人成| 夜夜看夜夜爽夜夜摸| 亚洲avbb在线观看| 九色成人免费人妻av| 一边摸一边做爽爽视频免费| 特大巨黑吊av在线直播| 亚洲狠狠婷婷综合久久图片| 老汉色∧v一级毛片| 国产片内射在线|