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

    物性參數(shù)溫度變化下激光熔覆多場(chǎng)耦合模擬與實(shí)驗(yàn)

    2019-08-06 01:50:30李昌于志斌高敬翔李云飛韓興
    兵工學(xué)報(bào) 2019年6期
    關(guān)鍵詞:熔池物性溫度場(chǎng)

    李昌, 于志斌, 高敬翔, 李云飛, 韓興

    (1.遼寧科技大學(xué) 機(jī)械工程與自動(dòng)化學(xué)院, 遼寧 鞍山 114051;2.中國(guó)能源建設(shè)集團(tuán) 東北電力第一工程有限公司, 遼寧 沈陽(yáng) 110179)

    0 引言

    激光熔覆涉及物理學(xué)、化學(xué)、冶金及材料科學(xué)等多學(xué)科交叉,熔覆過(guò)程是以不同送粉方式將熔覆粉置于基體表面,經(jīng)激光輻照使其與基體層同時(shí)熔化,并快速凝固成稀釋率極低、與基體成冶金結(jié)合的熔覆層,從而改善基體表面的耐磨、耐蝕、耐熱特性。激光熔覆材料具有稀釋率小、組織致密、涂層與基體結(jié)合性好、效率高、速度快、綠色環(huán)保等特點(diǎn),在航空航天、汽車、海洋、石油、化工等工業(yè)領(lǐng)域具有廣泛的應(yīng)用前景。

    20世紀(jì)60年代,隨著大功率激光器的發(fā)展,激光開(kāi)始進(jìn)入工業(yè)表面處理領(lǐng)域。早在1974年,Gnanamuthu[1]在基體上熔覆了一層金屬,取得了激光熔覆專利。側(cè)重于溫度場(chǎng)計(jì)算的研究有:Mazumder等[2]建立了激光熔覆溫度場(chǎng)準(zhǔn)穩(wěn)態(tài)系統(tǒng),為溫度場(chǎng)模擬提供了理論基礎(chǔ),但沒(méi)有考慮熔池內(nèi)相變潛熱對(duì)溫度場(chǎng)的影響; Brucker等[3]考慮粉末流與激光的相互作用,分析了熔覆溫度場(chǎng);Kou等[4]建立了準(zhǔn)穩(wěn)態(tài)激光熔凝流動(dòng)和導(dǎo)熱模型,計(jì)算了熔池形貌;Ehsan等[5]研究了沉積模式對(duì)激光熔覆熱應(yīng)力場(chǎng)的影響;Parisa等[6]對(duì)激光熔覆單層多道的溫度場(chǎng)、應(yīng)力場(chǎng)進(jìn)行了研究;Gao等[7]模擬激光熔覆過(guò)程中的溫度場(chǎng),尤其是溫度變化、冷卻速率和固體- 液體(簡(jiǎn)稱固液)界面凝固速率,研究了不同工藝下的溫度場(chǎng)變化規(guī)律,但未考慮Marangoni對(duì)流對(duì)熔池的作用。側(cè)重于流場(chǎng)計(jì)算的研究有:Picasso等[8]建立了激光熔覆二維數(shù)學(xué)模型,研究了熔體池中的流體運(yùn)動(dòng)和液體- 氣體(簡(jiǎn)稱液氣)界面的變形,但熔覆層高度需預(yù)先給定;Toyserkani等[9]考慮熔覆過(guò)程粉末顆粒對(duì)激光能量的衰減作用,建立了顆粒受熱溫升及熔池內(nèi)對(duì)流和熱傳導(dǎo)模型;Gan等[10]利用多物理場(chǎng)仿真軟件COMSOL建立了激光熔覆數(shù)學(xué)模型,包括溫度場(chǎng)、流場(chǎng),考慮到可移動(dòng)熔池的對(duì)流和擴(kuò)散。但這些研究均未考慮熔覆基體﹑粉材物性參數(shù)的溫度變化(簡(jiǎn)稱溫變)影響以及粉末對(duì)激光能量的衰減作用。

    激光熔覆過(guò)程涉及激光、粉末、基體間的交互作用。熔池內(nèi)的溫度場(chǎng)和流場(chǎng)將影響對(duì)流、傳熱和凝固,直接影響熔覆質(zhì)量。激光熔覆多場(chǎng)耦合相互影響演變規(guī)律如圖1所示。由于熔池體積小、溫變快、極強(qiáng)瞬時(shí)性特點(diǎn),難以用實(shí)驗(yàn)法跟蹤多場(chǎng)耦合動(dòng)態(tài)演變規(guī)律,而有限元法為研究熔覆機(jī)理提供了有效途徑。目前對(duì)激光熔覆有限元建模均未考慮基體和熔覆粉物性參數(shù)的溫變影響。而激光熔覆過(guò)程是熱- 彈性- 塑性- 流體(簡(jiǎn)稱熱- 彈- 塑- 流)多場(chǎng)耦合演化過(guò)程,建模中忽略基體和熔覆粉物性參數(shù)溫變影響將造成誤差。因此,本文建立了熔覆過(guò)程熱- 彈- 塑- 流多場(chǎng)耦合三維數(shù)學(xué)模型,以CALPHAD相圖法計(jì)算基體和粉材的溫變物性參數(shù),對(duì)模型求解,得出了熔覆過(guò)程中溫度場(chǎng)、流速場(chǎng)、凝固行為與演變規(guī)律。通過(guò)溫度梯度G和凝固生長(zhǎng)速率S預(yù)測(cè)其形態(tài)和凝固組織,為優(yōu)化激光熔覆工藝參數(shù)奠定重要理論基礎(chǔ)。

    1 實(shí)驗(yàn)材料與實(shí)驗(yàn)設(shè)備

    實(shí)驗(yàn)基體為45號(hào)鋼,退火處理,成分見(jiàn)表1,余量成分為Fe. 熔覆粉為Fe60,粒度為53-150目,成分見(jiàn)表2. 激光熔覆實(shí)驗(yàn)試件如圖2所示,激光熔覆過(guò)程原理如圖3所示。實(shí)驗(yàn)選取德國(guó)TRUMPF公司生產(chǎn)的TruDisk4002激光器,配以6自由度機(jī)器人系統(tǒng)實(shí)現(xiàn)激光熔覆,碟光激光器與6自由度機(jī)器人組合的激光熔覆系統(tǒng)如圖4所示。碟片激光器優(yōu)勢(shì)在于將固體激光器的棒狀晶改為碟片晶,碟片晶薄且直徑與厚度比大,可及時(shí)有效冷卻,一維熱傳導(dǎo)使晶體內(nèi)溫度分布均勻,解決了熱透鏡問(wèn)題,改善了光束質(zhì)量、轉(zhuǎn)換效率及功率穩(wěn)定性,比傳統(tǒng)激光器具有優(yōu)勢(shì)。所用激光器波長(zhǎng)為1 030 nm,激光功率80~4 000 W,輸出光束直徑0.2 mm,光束質(zhì)量為8 mm·mard,其轉(zhuǎn)換功率可達(dá)30%,比棒式激光器提高近10倍。將試件沿垂直掃描方向縱向切開(kāi),對(duì)剖面打磨、拋光,采用4%的硝酸酒精溶液腐蝕,用德國(guó)Carl Zeiss公司生產(chǎn)的Axioskop2電子顯微鏡和Zeiss-ΣIGMA HD場(chǎng)發(fā)射電子顯微鏡進(jìn)行金相實(shí)驗(yàn),觀察熔覆層金相組織形貌。

    表1 45號(hào)鋼的元素組成

    表2 Fe60的元素組成

    圖2 激光熔覆實(shí)驗(yàn)試件Fig.2 Experimental specimens for laser cladding

    圖3 激光熔覆過(guò)程原理示意圖Fig.3 Schematic diagram of laser cladding principle

    圖4 碟片激光器與6自由度機(jī)器人組合的激光熔覆系統(tǒng)Fig.4 Laser cladding system combined with a disk laser and a six degrees of freedom robot

    2 激光熔覆理論建模與求解

    激光熔覆過(guò)程建模基于以下假設(shè):

    1) 熔池金屬流假定為層流,且為不可壓縮牛頓流體;

    2) 激光束能量在光斑內(nèi)呈高斯分布,且功率恒定;

    3) 材料為各向同性;

    4) 粉末流濃度服從高斯分布,掉落到熔池中的粉末立即熔化。

    2.1 激光熔覆過(guò)程總體控制方程

    激光熔覆過(guò)程總體控制方程[11-13]包括連續(xù)質(zhì)量方程((1)式)和Navier-Stokes動(dòng)量方程((2)式)。

    (1)

    式中:ρ為密度;t為時(shí)間;u為熔池內(nèi)金屬流動(dòng)速度。

    (2)

    式中:μ為流體動(dòng)力黏度;p為壓強(qiáng);I為單位矩陣;K0為由多孔形態(tài)確定的常數(shù),取值為2×107;B為避免分母為0的極小參數(shù),取值為10-3.fl為液體質(zhì)量分?jǐn)?shù),可計(jì)算[14-15]為

    (3)

    T為溫度,Ts、Tl分別為固相、液相溫度。

    (2)式中:左邊為慣性力;右邊第1項(xiàng)為壓力,第2項(xiàng)為黏性力,最后1項(xiàng)為熔池與基體交接處糊狀區(qū)的動(dòng)量耗散,根據(jù)Carman-Kozeny方程進(jìn)行量化。

    (4)

    式中:cp為比熱容;k為導(dǎo)熱系數(shù);H為金屬熔化潛焓,即ΔH=Lfl,L為潛熱。 (4)式是能量方程:左邊第1項(xiàng)為熱累積,第2項(xiàng)為熱對(duì)流;右邊第1項(xiàng)為熱傳導(dǎo)。

    2.2 熔覆粉與激光束的相互作用

    本文采取同軸送粉方式,熔覆粉在光束中具有流動(dòng)發(fā)散性,送粉速率和載氣流量對(duì)發(fā)散性有直接影響,需要實(shí)現(xiàn)粉流束與激光束間的良好匹配。影響粉- 光(簡(jiǎn)稱粉光)匹配的參數(shù)有激光功率、離焦量、激光束與工件表面的相對(duì)移動(dòng)速度、送粉速率、載氣流量、熔覆顆粒大小、熔覆材料的物性參數(shù)。粉光匹配原則是熔覆粉在激光束的有效直徑內(nèi)加熱,并落到基體表面激光束作用有效區(qū)域內(nèi),如圖5所示[16]。工作表面處于3位置為最佳匹配狀態(tài);在2、3位置之間,熔覆層外觀良好,但光束有效直徑外的材料不能與基體實(shí)現(xiàn)冶金結(jié)合,降低了結(jié)合質(zhì)量;處于2位置以下不能實(shí)現(xiàn)熔覆。處于3、4位置間,靠近4位置會(huì)降低光束能量利用率,增大基體熱影響范圍;由4位置向5位置逼近,接近噴嘴,由于激光反射、被加熱粉的高溫輻射使噴嘴溫升,導(dǎo)致熔覆粉軟化粘連堵死噴嘴,移動(dòng)空間嚴(yán)重受限。

    圖5 激光熔覆粉光匹配示意圖Fig.5 Schematic diagram of matching of powder flow and laser light in the process of laser cladding

    激光熔覆中,當(dāng)光束穿過(guò)粉末流時(shí),其部分能量被粉末吸收、反射和散射,剩余能量穿過(guò)粉末流到達(dá)熔池表面。同時(shí)粉末吸收激光能量,在落入熔池前溫度已經(jīng)升高。根據(jù)Lambert-Beer光透射定理和米氏理論,粉末流對(duì)激光束的衰減作用與激光功率、波長(zhǎng)、粉末顆粒特性、粉末流濃度以及激光束穿越的粉末流長(zhǎng)度有關(guān),可表示[17-18]為

    q′(r,φ)=qexp(-δNφ),

    (5)

    式中:q′(r,φ)為距離噴嘴φ處平面上激光中心點(diǎn)的功率密度;r為消光面積;q為激光熱通量;δ為吸收率;N為單位體積顆粒數(shù)。

    粉末落入熔池前的溫升由熱平衡式得出:

    (6)

    式中:Radd為粉末顆粒半徑;ap為粉末對(duì)激光能量的吸收率;s為噴嘴到熔池的距離;vf為顆粒速度;ρm為粉末顆粒密度;cpm為粉末比熱容;ΔT′為粉末的溫升,

    (7)

    2.3 模型求解初始條件與邊界設(shè)置

    2.3.1 激光光源

    采用移動(dòng)高斯熱源,激光熱通量

    (8)

    式中:a為能量吸收率;P為激光功率;R為光斑半徑;v為掃描速度;x、y分別為光斑中心的瞬時(shí)坐標(biāo);hc為熱傳遞系數(shù);ε為發(fā)射率;σb為Stefan-Boltzmann常量;T0為初始環(huán)境溫度。

    綜上所述,本研究成功建立了新西蘭白兔減壓病模型,該模型中動(dòng)物肢體癱瘓率較高,可用于脊髓型減壓病研究。后續(xù)研究可根據(jù)實(shí)驗(yàn)?zāi)康?,調(diào)整高氣壓暴露的壓力、時(shí)間及減壓速率以制備出預(yù)期嚴(yán)重程度的減壓病動(dòng)物模型。建立的指標(biāo)評(píng)價(jià)體系能比較全面地反映動(dòng)物模型的發(fā)病情況,可根據(jù)不同的實(shí)驗(yàn)需要選擇部分或全部指標(biāo)評(píng)價(jià)干預(yù)措施的效能。

    2.3.2 邊界設(shè)置

    液氣動(dòng)量方程的邊界條件為

    (9)

    式中:Fl-g為液氣動(dòng)量;等號(hào)右側(cè)分別代表毛細(xì)力和熱毛細(xì)力;σ為表面張力;n*為曲面法向參數(shù);κ為曲面曲率。

    采取同軸送粉方式,考慮到液氣界面移動(dòng),邊界移動(dòng)速度可表述為

    vl-g=ul-g·n*+vp·n*,

    (10)

    式中:ul-g為在液氣界面的移動(dòng)速度;vp為因熔覆粉添加導(dǎo)致的液體和氣體界面的移動(dòng)速度,

    (11)

    式中:mf為送粉量;ηm為粉末流率;ρm為粉末密度;Rp為粉末流半徑;z為z軸方向單位向量。

    3 激光熔覆過(guò)程多場(chǎng)耦合數(shù)值模擬

    3.1 有限元模型的建立與參數(shù)設(shè)置

    圖6 激光熔覆有限元網(wǎng)格劃分Fig.6 Finite element mesh division of laser cladding

    基體為45號(hào)鋼,粉末為Fe60粉?;贑OMSOL Multiphysics平臺(tái)搭建熔覆過(guò)程熱- 彈- 塑- 流多場(chǎng)耦合有限元模型??紤]熔覆件左右對(duì)稱,建立1/2模型,尺寸為20 mm×10 mm×6 mm,采用自由四面體劃分網(wǎng)格,打開(kāi)網(wǎng)格自動(dòng)加密功能,網(wǎng)格包含315 988 個(gè)域單元、11 324 個(gè)邊界元和360個(gè)邊單元,如圖6所示。計(jì)算過(guò)程中選用熱應(yīng)力、層流、兩相流水平集、稀物質(zhì)傳遞模塊計(jì)算熔覆過(guò)程多場(chǎng)耦合變化規(guī)律,采用基于任意拉格朗日- 歐拉法的動(dòng)網(wǎng)格描述熔池動(dòng)態(tài)形狀變化,仿真模型與實(shí)驗(yàn)選取的具體工藝參數(shù)[19]如表3所示。

    表3 激光熔覆的計(jì)算參數(shù)

    注:λ0為激光波長(zhǎng)。

    工作初始溫度為環(huán)境溫度,熔池初速為0 m/s. 基體對(duì)光束能量吸收率與激光波長(zhǎng)、基體粗糙度有關(guān),且隨溫度變化,由Hagen-Rubens公式得

    (12)

    熔池的物性參數(shù)為

    Wbath=βWm+(1-β)Wp,

    (13)

    式中:Wbath為熔池?zé)嵛镄詤?shù);Wm和Wp分別為基體和熔覆粉的熱物性參數(shù);β為混合分?jǐn)?shù),取決于二者在熔覆層橫截面積占比。用CALPHAD相圖法計(jì)算基體、粉材的溫變物性參數(shù),分別如圖7、圖8所示。計(jì)算結(jié)果表明:基體和熔覆粉的物性參數(shù)受溫度影響較大,熔覆過(guò)程中忽略物性參數(shù)的溫變影響將給計(jì)算帶來(lái)較大誤差。用插值函數(shù)導(dǎo)入COMSOL中,設(shè)置求解和邊界條件,利用SIMPLE求解器求解。

    3.2 模擬計(jì)算與結(jié)果分析

    利用DELL T5600工作站對(duì)多場(chǎng)耦合模型進(jìn)行求解,得出熔覆過(guò)程多場(chǎng)耦合變化規(guī)律。

    圖7 基體物性參數(shù)隨溫度變化曲線Fig.7 Variations of physical properties of substrate with temperature

    圖8 熔覆粉Fe60的物性參數(shù)隨溫度變化曲線Fig.8 Variations of physical properties of cladding powder Fe60 with temperature

    3.2.1 溫度場(chǎng)模擬計(jì)算

    熔池溫度直接影響熔池內(nèi)相變凝固、晶體生長(zhǎng)、形核率以及微觀組織形成。模擬計(jì)算得到從0~700 ms熔覆溫度場(chǎng)變化規(guī)律如圖9所示。繪制熔覆過(guò)程中不同位置溫度變化曲線如圖10所示。圖9的計(jì)算表明:在初始階段,由于激光輻照產(chǎn)生急熱效應(yīng),加熱區(qū)溫度迅速升高,到700 ms熔池溫度達(dá)最高值2 590 K. 形成近似2.0 mm×1.5 mm×1.0 mm的橢球體熔池,最高溫度出現(xiàn)在移動(dòng)光斑中心偏后位置。當(dāng)最高溫度超過(guò)固相線后,基體發(fā)生固液相轉(zhuǎn)變,開(kāi)始形成熔池。各時(shí)刻熔池前方溫度等高線呈密集分布,溫度變化較為劇烈,后方溫度等高線分布較為稀疏,溫度變化較為平緩。圖9(c)是對(duì)熔池沿掃描方向等距剖分得出的熔池內(nèi)部溫度分布,結(jié)果表明:隨著粉末的不斷添加,導(dǎo)致熔覆層逐漸變厚,但每一時(shí)刻熔池最高溫度始終處于熔覆層頂面,不同剖面溫度圍繞光斑呈橢圓形分布,溫度由內(nèi)向外逐漸遞減。圖9(d)是將熔池沿掃描方向等距剖分觀察到的溫度分布,結(jié)果表明:隨著光斑的移動(dòng),縱向剖面內(nèi)溫度影響區(qū)同步前移,熱影響區(qū)剖面溫度呈橢圓分布,熔池表面溫度較高,由表及里沿熔池縱深方向溫度逐漸遞減。圖10的結(jié)果表明:沿著熔池深度的1號(hào)軌跡方向,靠近光斑溫度呈遞增趨勢(shì),變化梯度隨時(shí)間逐漸減小,最終趨于一致;沿著熔池y向的2號(hào)軌跡,靠近光斑溫度呈遞增趨勢(shì),增長(zhǎng)梯度逐漸減小,最終趨于一致;沿著掃描方向的3號(hào)軌跡,不同時(shí)刻的溫度先增大、后減小,光斑偏后位置溫度最高,增長(zhǎng)梯度高于減小梯度,最高溫度隨時(shí)間逐漸增大。

    圖9 激光熔覆溫度場(chǎng)演變規(guī)律Fig.9 Temperature field evolution of laser cladding

    圖10 激光熔覆不同位置溫度變化曲線Fig.10 Temperature change curves at different positions of laser cladding

    3.2.2 流速場(chǎng)模擬計(jì)算

    熔池液態(tài)金屬流對(duì)傳熱的影響,可由對(duì)流速率與擴(kuò)散速率之比PeT表示:

    (14)

    式中:L為熔池半徑的特征長(zhǎng)度;α1為熱擴(kuò)散系數(shù),

    (15)

    kr為熔池導(dǎo)熱系數(shù),ρr為熔池密度,cpr為熔池比熱容。隨著PeT的增大,熱傳遞中擴(kuò)散傳遞的比例減少,對(duì)流輸運(yùn)的比例增大。

    通過(guò)計(jì)算的流速場(chǎng)、PeT數(shù),可判定熔池內(nèi)熱傳導(dǎo)和熱對(duì)流哪個(gè)起主導(dǎo)作用。Marangoni效應(yīng)是熔池對(duì)流現(xiàn)象,如圖11所示為激光熔覆熔池Marangomi效應(yīng)示意圖。該效應(yīng)將剪切應(yīng)力的法向分量與溫度的切向?qū)?shù)相關(guān)聯(lián),用Marangoni數(shù)表征流動(dòng)特性,取決于幾何形狀。熔池內(nèi)Marangoni力會(huì)影響流體流動(dòng)和溫度分布,進(jìn)而改變?nèi)鄢匦螤睢?/p>

    圖11 激光熔覆熔池Marangoni效應(yīng)示意圖Fig.11 Schematic diagram of Marangoni effect in a laser cladding pool

    圖12 激光熔覆不同時(shí)刻的熔池速度流場(chǎng)Fig.12 Velocity fields in laser cladding pool at different times

    計(jì)算得出不同時(shí)刻的熔池速度流場(chǎng)如圖12所示。圖12的結(jié)果表明:由于表面張力的溫度系數(shù)為負(fù),熔池金屬流從激光束中心向熔池邊緣移動(dòng)。在t=10 ms接近初始階段,熔池金屬流速最大為6.64×10-3m/s,計(jì)算得到PeT數(shù)小于5,說(shuō)明熔池內(nèi)熱傳導(dǎo)起重要作用,形成半球形熔池邊界。隨著對(duì)流的增強(qiáng),當(dāng)t=700 ms時(shí),熔池內(nèi)金屬流速最大值達(dá)到0.3 m/s,計(jì)算得到PeT大于200,熔池內(nèi)熱對(duì)流起主要作用,熔池被強(qiáng)烈的Marangoni流所控制。圖13所示為激光熔覆不同位置的流速瞬時(shí)變化。圖13(a)表明:沿著1號(hào)軌跡線,遠(yuǎn)離熔池區(qū)流速為0 m/s,越接近光斑,流速呈現(xiàn)先增大、后減小的變化,最大流速逐漸變大。圖13(b)表明:沿著2號(hào)軌跡線,熔池內(nèi)部流速曲線直到300 ms以前,熔池內(nèi)流速變化趨勢(shì)相近,流速最大值不斷升高,在300 ms時(shí)刻最大流速達(dá)到0.19 m/s,呈現(xiàn)微小變動(dòng)。300 ms以后流速先增大、后減小,然后繼續(xù)增大,最大流速仍隨著時(shí)間逐漸變大,在700 ms時(shí)刻最大流速達(dá)到0.243 m/s. 圖13(c)表明:沿著3號(hào)軌跡即掃描方向,熔池流速呈現(xiàn)兩邊高、中間低的趨勢(shì),最高值出現(xiàn)在熔池前端,不同位置流速隨著時(shí)間逐漸增大。

    圖13 激光熔覆不同位置的流速瞬時(shí)變化Fig.13 Instantaneous change in velocity at different positions of laser cladding

    3.2.3 熔覆過(guò)程物性參數(shù)溫變影響分析

    物性參數(shù)隨溫度變化對(duì)熔覆過(guò)程溫度場(chǎng)、流速場(chǎng)的影響結(jié)果如圖14所示。圖14(a)~圖14(c)表明:考慮物性參數(shù)溫變影響下的溫度比未考慮物性參數(shù)溫變影響下的計(jì)算溫度偏低,二者變化趨勢(shì)一致。圖14(d)~圖14(f)表明:考慮物性參數(shù)溫變影響下的流速比未考慮物性參數(shù)溫變影響下的流速偏低,二者變化趨勢(shì)一致。物性參數(shù)溫變對(duì)計(jì)算結(jié)果有一定的影響,計(jì)算中忽略物性參數(shù)的溫變影響將導(dǎo)致計(jì)算產(chǎn)生誤差。

    3.3 激光熔覆金相實(shí)驗(yàn)

    圖15描述了冷卻速率、溫度梯度和凝固速度對(duì)凝固微觀組織結(jié)構(gòu)和尺寸的影響。由圖15可見(jiàn):影響凝固微觀組織的兩個(gè)關(guān)鍵因素是凝固前溫度梯度G和固/液界面推移速率S;冷卻速率G×S影響凝固組織的尺寸大小,較高的冷卻速率導(dǎo)致組織更加精細(xì);G/S決定了微觀組織從平面晶、胞晶、柱狀枝晶、等軸晶的形態(tài)變化。

    (16)

    式中:ΔTd為單位距離上溫度變化;Δd為距離變化;θ為固液界面推移速率(凝固生長(zhǎng)速率)與激光掃描速度方向夾角。接近熔池與基體相交底部,S與v近似垂直,越靠近熔池頂端,θ角越小,S越大。G由模擬計(jì)算得出,如圖15所示。試件沿垂直掃描方向剖開(kāi),用Axioskop2 SEM電子顯微鏡、Zeiss-ΣIGMA HD場(chǎng)發(fā)射電子顯微鏡觀察熔覆層顯微組織形貌,如圖16所示。根據(jù)G與計(jì)算得出的S,可導(dǎo)出冷卻速率G×S和形狀控制因子G/S:A點(diǎn)G×S=126,G/S=5 099;B點(diǎn)G×S=1 157,G/S=864;C點(diǎn)G×S=1 996,G/S=698;進(jìn)而可預(yù)測(cè)熔覆中熔池凝固組織的形態(tài)和尺寸。

    模擬與實(shí)驗(yàn)結(jié)果對(duì)比表明:A點(diǎn)位置處于熔池底部,具有較高的G/R和較低的G×R,處于平面結(jié)晶區(qū);C點(diǎn)位置處于熔池頂部,具有較低的G/R和較高的G×R,為等軸晶與枝晶共存區(qū),觀察可知為等軸晶和較細(xì)的枝晶分布;B點(diǎn)位置處于枝晶區(qū),觀察可知為柱狀和粗枝狀枝晶;實(shí)驗(yàn)結(jié)果與仿真趨勢(shì)一致。

    圖14 激光熔覆中物性參數(shù)溫變對(duì)溫度場(chǎng)和流場(chǎng)的影響Fig.14 Influences of physical parameters with temperature on temperature field and flow field in the process of laser cladding

    圖15 激光熔覆中G和S對(duì)熔池凝固組織的影響Fig.15 Effects of G and S on the morphology of solidifled microstructure in the process of laser cladding

    圖16 激光熔覆過(guò)程中G變化Fig.16 Cloud charts of G during laser cladding

    圖17 熔覆層的實(shí)驗(yàn)剖面與A、B、C點(diǎn)位置的微觀組織Fig.17 Microstructures at the positions A, B and C on experimental cross-section of cladding layer

    4 結(jié)論

    1) 采用碟片激光器在45號(hào)鋼基體熔覆Fe60粉,加熱區(qū)溫度迅速升高,形成近似2.0 mm×1.5 mm×1.0 mm的橢球體熔池,最高溫度出現(xiàn)在移動(dòng)光斑中心偏后位置;熔覆初期,熔池金屬流速最大值為6.64×10-3m/s,熔池內(nèi)熱傳導(dǎo)起主要作用。700 ms時(shí)最大金屬流速為0.3 m/s,熔池?zé)釋?duì)流起主要作用。

    2) 激光熔覆熔池底部處于平面結(jié)晶區(qū),頂部為等軸晶與較細(xì)枝晶的共存區(qū),中部為粗大的枝晶。模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)吻合良好,驗(yàn)證了本文方法的有效性。

    猜你喜歡
    熔池物性溫度場(chǎng)
    R1234ze PVTx熱物性模擬計(jì)算
    能源工程(2022年1期)2022-03-29 01:06:26
    中韓天氣預(yù)報(bào)語(yǔ)篇的及物性分析
    LKP狀態(tài)方程在天然氣熱物性參數(shù)計(jì)算的應(yīng)用
    煤氣與熱力(2021年6期)2021-07-28 07:21:30
    鋁合金加筋板焊接溫度場(chǎng)和殘余應(yīng)力數(shù)值模擬
    基于紋影法的溫度場(chǎng)分布測(cè)量方法
    MJS工法與凍結(jié)法結(jié)合加固區(qū)溫度場(chǎng)研究
    建筑科技(2018年6期)2018-08-30 03:41:08
    低孔低滲儲(chǔ)層物性下限確定方法及其適用性
    電弧焊熔池表征與熔透狀態(tài)映射研究
    焊接(2016年10期)2016-02-27 13:05:32
    一種基于頻域變換的熔池運(yùn)動(dòng)模糊圖像恢復(fù)算法
    焊接(2015年11期)2015-07-18 11:12:46
    MIG焊接熔池形成與凝固過(guò)程數(shù)值模擬
    焊接(2015年9期)2015-07-18 11:03:51
    精品国产一区二区三区久久久樱花| 亚洲 欧美一区二区三区| 一边摸一边抽搐一进一出视频| 自线自在国产av| 国产精品二区激情视频| 欧美精品一区二区免费开放| 50天的宝宝边吃奶边哭怎么回事| 精品国产超薄肉色丝袜足j| 另类亚洲欧美激情| 亚洲精华国产精华精| 在线av久久热| av电影中文网址| 亚洲av成人不卡在线观看播放网 | 中国国产av一级| 久久精品亚洲熟妇少妇任你| 欧美在线黄色| 天天操日日干夜夜撸| 一本色道久久久久久精品综合| 一级毛片精品| 亚洲av成人不卡在线观看播放网 | 亚洲va日本ⅴa欧美va伊人久久 | 亚洲天堂av无毛| 高清在线国产一区| 超碰97精品在线观看| 欧美激情 高清一区二区三区| 久久久久久久久久久久大奶| 亚洲av电影在线观看一区二区三区| 黄色视频在线播放观看不卡| 成年动漫av网址| 亚洲国产精品一区二区三区在线| 夫妻午夜视频| 97人妻天天添夜夜摸| 日韩视频在线欧美| 亚洲av国产av综合av卡| 精品少妇久久久久久888优播| 亚洲成人免费电影在线观看| 亚洲欧美成人综合另类久久久| 各种免费的搞黄视频| 欧美另类亚洲清纯唯美| 久久精品亚洲熟妇少妇任你| 另类亚洲欧美激情| 十八禁人妻一区二区| 亚洲精品国产色婷婷电影| 国产精品免费视频内射| 一本一本久久a久久精品综合妖精| 又大又爽又粗| 精品一区二区三区av网在线观看 | 国产在线一区二区三区精| 国产免费av片在线观看野外av| 亚洲视频免费观看视频| 不卡av一区二区三区| 又紧又爽又黄一区二区| 人人妻,人人澡人人爽秒播| 久久天堂一区二区三区四区| 成人免费观看视频高清| 亚洲国产精品一区二区三区在线| 国产精品偷伦视频观看了| 日韩有码中文字幕| 亚洲第一欧美日韩一区二区三区 | 免费不卡黄色视频| 日日摸夜夜添夜夜添小说| 久久久久久人人人人人| 日韩视频一区二区在线观看| 成人三级做爰电影| 亚洲欧美激情在线| 日韩熟女老妇一区二区性免费视频| 操出白浆在线播放| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人免费av在线播放| 免费观看av网站的网址| 老司机午夜福利在线观看视频 | 国产亚洲一区二区精品| 国产欧美日韩精品亚洲av| 亚洲色图 男人天堂 中文字幕| 淫妇啪啪啪对白视频 | 亚洲精品自拍成人| 亚洲欧美日韩高清在线视频 | 十八禁网站免费在线| 亚洲熟女精品中文字幕| 在线观看免费高清a一片| 好男人电影高清在线观看| 久久女婷五月综合色啪小说| 一级片'在线观看视频| 亚洲一区中文字幕在线| 亚洲av美国av| 亚洲国产av新网站| 天天躁日日躁夜夜躁夜夜| 丝袜喷水一区| av网站免费在线观看视频| 国产成人啪精品午夜网站| 一级a爱视频在线免费观看| 欧美在线黄色| 90打野战视频偷拍视频| 亚洲 国产 在线| 国产一级毛片在线| 久久久国产欧美日韩av| 男女床上黄色一级片免费看| 色94色欧美一区二区| e午夜精品久久久久久久| 中文字幕精品免费在线观看视频| 午夜福利一区二区在线看| 亚洲精品日韩在线中文字幕| 啦啦啦在线免费观看视频4| 成人国产一区最新在线观看| 国产精品av久久久久免费| 91麻豆av在线| 亚洲熟女毛片儿| 欧美黑人欧美精品刺激| 国产精品秋霞免费鲁丝片| 午夜91福利影院| av天堂在线播放| 国产伦人伦偷精品视频| 国产精品一区二区在线不卡| 黑人巨大精品欧美一区二区蜜桃| 国内毛片毛片毛片毛片毛片| 美女扒开内裤让男人捅视频| 婷婷色av中文字幕| 国产有黄有色有爽视频| 高清在线国产一区| 嫩草影视91久久| 少妇裸体淫交视频免费看高清 | 国产视频一区二区在线看| 波多野结衣av一区二区av| 国产日韩一区二区三区精品不卡| 亚洲精品久久午夜乱码| 午夜福利视频精品| 国产亚洲av高清不卡| 自线自在国产av| 午夜影院在线不卡| 捣出白浆h1v1| 亚洲国产欧美一区二区综合| 日本vs欧美在线观看视频| 日日夜夜操网爽| 午夜精品国产一区二区电影| 国产一区二区三区综合在线观看| 亚洲七黄色美女视频| 波多野结衣一区麻豆| 亚洲第一欧美日韩一区二区三区 | 制服诱惑二区| 国产片内射在线| av电影中文网址| 纯流量卡能插随身wifi吗| 亚洲国产精品999| 亚洲av电影在线进入| 大陆偷拍与自拍| 国内毛片毛片毛片毛片毛片| 一本—道久久a久久精品蜜桃钙片| 欧美人与性动交α欧美软件| 国产在线一区二区三区精| 国内毛片毛片毛片毛片毛片| 久久精品成人免费网站| 久久午夜综合久久蜜桃| 另类精品久久| 日本精品一区二区三区蜜桃| 久久国产精品大桥未久av| 性少妇av在线| 亚洲自偷自拍图片 自拍| 18禁观看日本| 午夜福利,免费看| 汤姆久久久久久久影院中文字幕| 日本猛色少妇xxxxx猛交久久| 狂野欧美激情性bbbbbb| 中文字幕人妻丝袜一区二区| 人人妻人人澡人人看| 精品少妇内射三级| 国产日韩欧美在线精品| 人成视频在线观看免费观看| 亚洲国产av新网站| 大片免费播放器 马上看| 免费在线观看影片大全网站| 1024香蕉在线观看| 51午夜福利影视在线观看| 69精品国产乱码久久久| 啦啦啦中文免费视频观看日本| 热re99久久精品国产66热6| 久久天堂一区二区三区四区| 一区二区三区精品91| 伦理电影免费视频| 成人影院久久| 我要看黄色一级片免费的| 午夜免费观看性视频| 美女大奶头黄色视频| 成人国产av品久久久| 亚洲 欧美一区二区三区| 嫩草影视91久久| 亚洲九九香蕉| 亚洲成国产人片在线观看| 777久久人妻少妇嫩草av网站| 一区二区日韩欧美中文字幕| 国产免费福利视频在线观看| 国产成人av教育| 男男h啪啪无遮挡| 欧美日韩亚洲高清精品| 人妻一区二区av| 狠狠狠狠99中文字幕| 秋霞在线观看毛片| 伊人久久大香线蕉亚洲五| 亚洲久久久国产精品| 在线观看免费午夜福利视频| 色视频在线一区二区三区| 日韩中文字幕欧美一区二区| 国产免费av片在线观看野外av| 免费观看人在逋| 菩萨蛮人人尽说江南好唐韦庄| 免费一级毛片在线播放高清视频 | 男女高潮啪啪啪动态图| 欧美日韩av久久| tube8黄色片| 午夜福利在线免费观看网站| 精品熟女少妇八av免费久了| 精品少妇久久久久久888优播| 美女视频免费永久观看网站| 精品人妻熟女毛片av久久网站| 亚洲色图 男人天堂 中文字幕| 久久 成人 亚洲| 久久亚洲精品不卡| 久久这里只有精品19| 九色亚洲精品在线播放| 少妇裸体淫交视频免费看高清 | 欧美黑人精品巨大| 欧美黄色淫秽网站| 老汉色∧v一级毛片| 在线看a的网站| 久久青草综合色| 国产精品香港三级国产av潘金莲| 亚洲人成电影免费在线| 在线观看www视频免费| 中文字幕另类日韩欧美亚洲嫩草| 国产成人a∨麻豆精品| 久久久欧美国产精品| 亚洲一区二区三区欧美精品| 久久青草综合色| 日韩电影二区| a在线观看视频网站| 久久99热这里只频精品6学生| 日韩视频在线欧美| 首页视频小说图片口味搜索| 婷婷色av中文字幕| 十八禁人妻一区二区| 免费在线观看视频国产中文字幕亚洲 | 国产精品成人在线| 黄色视频在线播放观看不卡| 男女之事视频高清在线观看| 啦啦啦在线免费观看视频4| 久久久精品免费免费高清| 性少妇av在线| 999久久久精品免费观看国产| 成年美女黄网站色视频大全免费| 国精品久久久久久国模美| 久久精品aⅴ一区二区三区四区| 国产精品秋霞免费鲁丝片| 亚洲熟女毛片儿| 日本a在线网址| 中文字幕高清在线视频| 亚洲精品美女久久久久99蜜臀| 老司机福利观看| 狠狠精品人妻久久久久久综合| 视频在线观看一区二区三区| a 毛片基地| 最近中文字幕2019免费版| 久久人人爽av亚洲精品天堂| 亚洲精品一二三| 国产精品亚洲av一区麻豆| 国产一级毛片在线| 中文欧美无线码| 精品福利观看| 天天操日日干夜夜撸| 久久精品国产综合久久久| 99re6热这里在线精品视频| www.精华液| 无遮挡黄片免费观看| 欧美人与性动交α欧美精品济南到| 欧美精品亚洲一区二区| 男女下面插进去视频免费观看| 国产成人精品久久二区二区免费| 亚洲精品美女久久av网站| 男女国产视频网站| 精品人妻在线不人妻| 夜夜夜夜夜久久久久| 欧美 日韩 精品 国产| 又大又爽又粗| av在线播放精品| 国产精品国产av在线观看| 狂野欧美激情性xxxx| 久久人人97超碰香蕉20202| 亚洲五月婷婷丁香| 国产亚洲av高清不卡| 99国产精品一区二区蜜桃av | 淫妇啪啪啪对白视频 | 少妇的丰满在线观看| 亚洲专区中文字幕在线| 最黄视频免费看| 国产日韩欧美视频二区| 天天操日日干夜夜撸| 天天躁日日躁夜夜躁夜夜| 亚洲精品粉嫩美女一区| 精品国内亚洲2022精品成人 | 日韩一区二区三区影片| 1024视频免费在线观看| 亚洲中文日韩欧美视频| 欧美老熟妇乱子伦牲交| 91大片在线观看| 男人添女人高潮全过程视频| 色综合欧美亚洲国产小说| 欧美+亚洲+日韩+国产| 男人操女人黄网站| 亚洲精品在线美女| 男女之事视频高清在线观看| 精品一区在线观看国产| 另类亚洲欧美激情| 一本色道久久久久久精品综合| 9热在线视频观看99| 亚洲人成77777在线视频| 脱女人内裤的视频| 人妻一区二区av| 精品一区在线观看国产| 蜜桃在线观看..| 法律面前人人平等表现在哪些方面 | 香蕉国产在线看| 精品少妇一区二区三区视频日本电影| 久久免费观看电影| 一二三四在线观看免费中文在| 在线看a的网站| 午夜福利影视在线免费观看| 国产亚洲精品久久久久5区| 午夜91福利影院| 久久久国产精品麻豆| av国产精品久久久久影院| 高清视频免费观看一区二区| 精品一区二区三区av网在线观看 | 国产精品一区二区在线观看99| 精品少妇久久久久久888优播| 国产91精品成人一区二区三区 | 国产精品免费视频内射| 亚洲精品日韩在线中文字幕| 亚洲熟女毛片儿| 又大又爽又粗| 精品国产超薄肉色丝袜足j| 国产成人影院久久av| 97精品久久久久久久久久精品| 69av精品久久久久久 | 男人添女人高潮全过程视频| 国产男女内射视频| 久久亚洲精品不卡| 久久久久精品国产欧美久久久 | 日本vs欧美在线观看视频| 久久人人97超碰香蕉20202| 欧美日韩一级在线毛片| 嫩草影视91久久| 三级毛片av免费| 狠狠精品人妻久久久久久综合| 久久香蕉激情| 成人国产一区最新在线观看| 99久久精品国产亚洲精品| 国产精品影院久久| 日韩大片免费观看网站| 国产又爽黄色视频| √禁漫天堂资源中文www| 亚洲国产看品久久| 一级毛片电影观看| 欧美激情 高清一区二区三区| 国产精品一区二区在线不卡| 深夜精品福利| 啦啦啦中文免费视频观看日本| 精品人妻熟女毛片av久久网站| 我的亚洲天堂| www.av在线官网国产| 国产99久久九九免费精品| 日韩,欧美,国产一区二区三区| 黄色a级毛片大全视频| 人人妻人人爽人人添夜夜欢视频| 黄色怎么调成土黄色| 国产欧美日韩精品亚洲av| 久久99热这里只频精品6学生| 精品少妇黑人巨大在线播放| 国产99久久九九免费精品| 国产欧美亚洲国产| 9191精品国产免费久久| 中文字幕人妻丝袜制服| 亚洲专区字幕在线| 交换朋友夫妻互换小说| 一边摸一边做爽爽视频免费| 青草久久国产| 999久久久国产精品视频| 99久久国产精品久久久| 久久久国产欧美日韩av| 午夜激情久久久久久久| 一区二区日韩欧美中文字幕| 国产日韩一区二区三区精品不卡| 国产一区二区三区综合在线观看| 99久久综合免费| 一个人免费在线观看的高清视频 | 日本黄色日本黄色录像| 大片免费播放器 马上看| 黑丝袜美女国产一区| 国产成+人综合+亚洲专区| 少妇裸体淫交视频免费看高清 | 国产精品一区二区在线观看99| 国产免费一区二区三区四区乱码| 成年美女黄网站色视频大全免费| 九色亚洲精品在线播放| 中文欧美无线码| 少妇粗大呻吟视频| 亚洲国产av新网站| 国产又爽黄色视频| 正在播放国产对白刺激| 国产在线视频一区二区| 欧美日韩国产mv在线观看视频| 久久人妻福利社区极品人妻图片| 老司机亚洲免费影院| 美女大奶头黄色视频| 99热全是精品| 欧美日韩福利视频一区二区| 男女午夜视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 国产日韩一区二区三区精品不卡| 日本黄色日本黄色录像| 精品久久久精品久久久| tocl精华| 日韩制服丝袜自拍偷拍| 一本—道久久a久久精品蜜桃钙片| 国产又爽黄色视频| 国产精品一区二区在线不卡| 亚洲精品粉嫩美女一区| 亚洲精品国产av成人精品| 国产成人精品无人区| 又紧又爽又黄一区二区| 别揉我奶头~嗯~啊~动态视频 | 一级片'在线观看视频| av视频免费观看在线观看| 欧美久久黑人一区二区| 色94色欧美一区二区| 国产成+人综合+亚洲专区| 永久免费av网站大全| 侵犯人妻中文字幕一二三四区| 欧美精品一区二区免费开放| 99久久人妻综合| 精品久久蜜臀av无| 欧美变态另类bdsm刘玥| av国产精品久久久久影院| 老司机影院毛片| 国产精品久久久久成人av| 一级片免费观看大全| 久久久精品区二区三区| 日本wwww免费看| 波多野结衣av一区二区av| 午夜福利视频精品| 女性生殖器流出的白浆| 欧美人与性动交α欧美软件| 最近最新免费中文字幕在线| 动漫黄色视频在线观看| 一进一出抽搐动态| av在线老鸭窝| 国产黄色免费在线视频| 黄色 视频免费看| 国产精品秋霞免费鲁丝片| 搡老岳熟女国产| 国产高清videossex| 99re6热这里在线精品视频| 美女福利国产在线| 亚洲av日韩在线播放| 久久人人爽人人片av| 男人操女人黄网站| 国产野战对白在线观看| 夫妻午夜视频| 亚洲国产日韩一区二区| 我要看黄色一级片免费的| 男女边摸边吃奶| 老汉色∧v一级毛片| 99热网站在线观看| 日韩精品免费视频一区二区三区| 淫妇啪啪啪对白视频 | 高潮久久久久久久久久久不卡| 国产精品.久久久| 国产成+人综合+亚洲专区| 最近最新免费中文字幕在线| 国产高清视频在线播放一区 | 热99久久久久精品小说推荐| 丁香六月欧美| 精品国产一区二区三区久久久樱花| 曰老女人黄片| 国产片内射在线| 考比视频在线观看| 啦啦啦视频在线资源免费观看| 狠狠婷婷综合久久久久久88av| 国产精品免费视频内射| 欧美日韩福利视频一区二区| 亚洲第一av免费看| 91老司机精品| 欧美日韩亚洲综合一区二区三区_| 国产精品一区二区精品视频观看| av福利片在线| 精品一区二区三卡| xxxhd国产人妻xxx| 男女之事视频高清在线观看| 天堂俺去俺来也www色官网| 亚洲欧美成人综合另类久久久| 男人添女人高潮全过程视频| 高潮久久久久久久久久久不卡| 国产有黄有色有爽视频| 国产亚洲av高清不卡| 国产国语露脸激情在线看| 亚洲第一av免费看| 我的亚洲天堂| 性色av乱码一区二区三区2| 久久国产精品影院| 日本wwww免费看| 欧美黑人精品巨大| 91精品国产国语对白视频| 亚洲欧美一区二区三区黑人| 婷婷成人精品国产| 少妇的丰满在线观看| 国产伦理片在线播放av一区| 午夜福利影视在线免费观看| 一区二区三区精品91| 极品少妇高潮喷水抽搐| av一本久久久久| 90打野战视频偷拍视频| 亚洲专区国产一区二区| 亚洲国产av影院在线观看| 精品人妻在线不人妻| 国产1区2区3区精品| 亚洲av欧美aⅴ国产| 三上悠亚av全集在线观看| 久久久国产一区二区| 国产一区二区三区综合在线观看| www日本在线高清视频| 中文欧美无线码| 一个人免费看片子| 色视频在线一区二区三区| 天堂中文最新版在线下载| 777久久人妻少妇嫩草av网站| 精品国产乱码久久久久久小说| 日韩 欧美 亚洲 中文字幕| 脱女人内裤的视频| 交换朋友夫妻互换小说| 91字幕亚洲| 免费av中文字幕在线| av有码第一页| 一二三四在线观看免费中文在| 国产成人一区二区三区免费视频网站| 国产免费av片在线观看野外av| 丝袜美足系列| 汤姆久久久久久久影院中文字幕| 国产亚洲av片在线观看秒播厂| 久久热在线av| 51午夜福利影视在线观看| 黄色 视频免费看| 人妻久久中文字幕网| 亚洲精品美女久久久久99蜜臀| 热99国产精品久久久久久7| 免费观看a级毛片全部| 男人添女人高潮全过程视频| 免费一级毛片在线播放高清视频 | 欧美少妇被猛烈插入视频| 女性生殖器流出的白浆| 精品少妇一区二区三区视频日本电影| e午夜精品久久久久久久| 精品乱码久久久久久99久播| 两性午夜刺激爽爽歪歪视频在线观看 | 高清欧美精品videossex| 青草久久国产| 欧美av亚洲av综合av国产av| 国产成人av激情在线播放| 最近最新中文字幕大全免费视频| 少妇粗大呻吟视频| 国产精品 欧美亚洲| 人人妻人人澡人人爽人人夜夜| 国产亚洲av高清不卡| 一本—道久久a久久精品蜜桃钙片| 99国产精品免费福利视频| www日本在线高清视频| 久久精品久久久久久噜噜老黄| 久久久精品免费免费高清| 在线观看舔阴道视频| 超色免费av| 亚洲第一av免费看| 在线观看免费高清a一片| 最新的欧美精品一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲 欧美一区二区三区| 久久免费观看电影| 精品亚洲成国产av| 久久青草综合色| 亚洲精华国产精华精| 精品亚洲成a人片在线观看| 黄色怎么调成土黄色| 久久久久久久久久久久大奶| 亚洲少妇的诱惑av| 亚洲精品久久午夜乱码| 日韩人妻精品一区2区三区| 男女午夜视频在线观看| 日本av手机在线免费观看| 电影成人av| 丁香六月天网| 窝窝影院91人妻| 久久久久国内视频| 天天影视国产精品| 免费一级毛片在线播放高清视频 | 两个人看的免费小视频| 国产成人精品久久二区二区91| 黑人巨大精品欧美一区二区蜜桃| 9色porny在线观看| 纵有疾风起免费观看全集完整版| 成人国产一区最新在线观看| 国产一区二区 视频在线| 大型av网站在线播放| 成人国产一区最新在线观看| 9色porny在线观看| 啪啪无遮挡十八禁网站| 亚洲伊人色综图| 亚洲精品久久午夜乱码| 视频在线观看一区二区三区| 午夜影院在线不卡| 大片免费播放器 马上看|