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

    基于多尺度建模與響應(yīng)面法的55NiCrMoV7模具鋼調(diào)質(zhì)熱處理工藝優(yōu)化

    2022-02-22 06:43:26趙繼全曲周德武川石如星禹興勝
    精密成形工程 2022年2期
    關(guān)鍵詞:棒料模具鋼調(diào)質(zhì)

    趙繼全,曲周德,武川,,石如星,禹興勝

    基于多尺度建模與響應(yīng)面法的55NiCrMoV7模具鋼調(diào)質(zhì)熱處理工藝優(yōu)化

    趙繼全1,曲周德1,武川1,2,3,石如星2,3,禹興勝2,3

    (1. 天津職業(yè)技術(shù)師范大學 汽車模具智能制造國家地方聯(lián)合工程實驗室,天津 300222;2. 洛陽中重鑄鍛有限責任公司,河南 洛陽 471039;3. 中信重工機械股份有限公司,河南 洛陽 471039)

    探究大型55NiCrMoV7模具鋼鍛件(1000 mm×1700 mm)調(diào)質(zhì)熱處理的最優(yōu)工藝方案。建立大型55NiCrMoV7模具鋼鍛件的宏微觀耦合的多尺度有限元模型,基于現(xiàn)場實驗測試不同位置溫度場和組織場的演化規(guī)律驗證該模型。模擬研究此鍛件在不同淬火與回火工藝下的應(yīng)力分布、HRC硬度和屈服強度變化。以殘余應(yīng)力均方差、硬度和屈服強度為目標函數(shù),建立淬、回火工藝參數(shù)與力學性能之間的回歸模型。結(jié)合所建立的響應(yīng)面模型分析與遺傳算法,對該模型進行參量約束條件下的多目標優(yōu)化。通過響應(yīng)面模型分析得出最優(yōu)工藝范圍為:淬火溫度為850~910 ℃,回火溫度為400~550 ℃。多目標優(yōu)化后得到Pareto最優(yōu)解集:淬火907.4 ℃+回火443.1 ℃、淬火907.4 ℃+回火411.5 ℃、淬火903.4 ℃+回火485.1 ℃、淬火895.2 ℃+回火535.8 ℃和淬火892.1 ℃+回火507.1 ℃。建立的淬、回火溫度對應(yīng)力均方差、HRC硬度和屈服強度的回歸模型,為大型55NiCrMoV7模具鋼鍛件調(diào)質(zhì)熱處理工藝優(yōu)化提供了理論依據(jù)。

    55NiCrMoV7模具鋼;多目標優(yōu)化;響應(yīng)面法;遺傳算法;多尺度模型

    55NiCrMoV7模具鋼作為一種熱作模具鋼,因其高溫下具有良好的綜合力學性能,被廣泛應(yīng)用于各種熱作模具沖壓制造領(lǐng)域[1]。目前,Jebaraj M等[2-3]研究了LN2和CO2冷卻劑在55NiCrMoV7鋼銑削中的作用,對55NiCrMoV7模具鋼在低溫環(huán)境下端銑加工表面完整性進行了探究。S. Hoja等[4]對55NiCrMoV7模具鋼碳氮共滲處理的發(fā)展進行了研究。Weifang Zhang等[5]對大型模鍛壓力機上的模具座(由55NiCrMoV7鋼制成)的脆性斷裂行為進行了研究,指出在使用過程中產(chǎn)生的大尺寸碳化物是導致材料脆化和沖擊韌性退化,進而導致模座脆性斷裂的主要原因。在熱處理工藝方面,張占平等[6-7]提出了55NiCrMoV7鋼在復(fù)雜回火工藝條件下的動力學模型,建立了回火組織特征馬氏體(211)峰物理半高寬-時間-回火溫度的單一模型和分段模型。Mengting Fan等[8-9]對55NiCrMoV7鋼進行了不同速率的淬火處理,并進行了600,620,650 ℃回火探究。禹興勝等[10]模擬研究了65 t的55NiCrMoV7鋼錠均質(zhì)熱處理過程中宏觀偏析擴散規(guī)律,得出鉻元素偏析是影響鋼錠均勻化的主要因素,為制定均質(zhì)熱處理溫度和時間提供了科學指導。

    對大型鋼錠進行恰當?shù)恼{(diào)質(zhì)熱處理是提高其綜合質(zhì)量的關(guān)鍵,目前針對大型55NiCrMoV7模具鋼鍛件的調(diào)質(zhì)工藝優(yōu)化還鮮有報道[11]。因此文中將對55NiCrMoV7鋼大型鍛件進行調(diào)質(zhì)工藝優(yōu)化的探究,以改善其力學性能。為了得到調(diào)質(zhì)工藝與力學性能之間的關(guān)系,文中擬采用響應(yīng)面法建立響應(yīng)面模型,并利用Matlab遺傳算法工具箱進行多目標優(yōu)化,相關(guān)研究結(jié)論為大型鍛件熱處理工藝的優(yōu)化提供了理論依據(jù)和科學指導。

    1 材料與實驗

    實驗采用55NiCrMoV7鋼錠鍛制成規(guī)格為1000 mm×1700 mm的棒料,質(zhì)量為10 625 kg,具體化學成分如表1所示。

    實驗熱處理流程:對55NiCrMoV7模具鋼進行890 ℃×10 h退火+740 ℃×24 h高溫回火的鍛后熱處理工藝,以消除內(nèi)應(yīng)力,改善材料內(nèi)部組織,為調(diào)質(zhì)熱處理做好組織準備,在鍛后熱處理過程中每隔7 min記錄一次溫度變化。調(diào)質(zhì)熱處理工藝采用860 ℃×16 h淬火+560 ℃×40 h回火,并對鍛件進行金相組織觀察、硬度及屈服強度測試,為大型55NiCrMoV7模具鋼鍛件的模型驗證提供基礎(chǔ)。

    模擬調(diào)質(zhì)熱處理工藝可提供大量數(shù)據(jù)支撐,避免了純實驗所帶來的巨大工作量,提高了工作效率。其工藝流程如下:≤200 ℃起溫加熱至淬火溫度(790~ 910 ℃)后,空氣預(yù)冷3 min+油冷200 min+上下串動10 min,淬火前油溫≤60 ℃,回火(100~650 ℃)并出爐空冷至室溫。控制升溫速度為60 ℃/h,且采用階梯加熱的方式,具體調(diào)質(zhì)工藝如圖1所示。在模擬調(diào)質(zhì)熱處理結(jié)果中提取殘余應(yīng)力值、HRC硬度和屈服強度值,并計算出殘余應(yīng)力的均方差,以衡量殘余應(yīng)力在大型鍛件內(nèi)部分布的均勻性。

    表1 實驗材料55NiCrMoV7鋼化學成分(質(zhì)量分數(shù))

    Tab.1 Chemical component of 55NiCrMoV7 steel (mass fraction) %

    2 模型的建立

    2.1 大型55NiCrMoV7鋼調(diào)質(zhì)熱處理的多尺度模型

    基于DEFORM對材料構(gòu)建的模型如圖2所示,網(wǎng)格劃分以六面體為主,邊角處采用四面體網(wǎng)格細化局部特征,熱處理過程中選用第3類邊界條件,試件與外部環(huán)境通過熱對流和熱輻射進行熱交換,空淬取恒溫25 ℃,油冷取恒溫50 ℃,考慮到試件的受熱情況,升溫階段熱交換系數(shù)設(shè)為0.1,降溫階段的換熱系數(shù)通過實驗測得,從試樣心部到表層均勻選取32個點以提取殘余應(yīng)力值。

    圖1 模擬調(diào)質(zhì)工藝流程

    圖2 55NiCrMoV7鋼模型的建立

    2.2 有限元模型的驗證

    2.2.1 溫度場驗證

    為了驗證有限元模型的可靠性,在鍛后熱處理過程中采用熱電偶對55NiCrMoV7模具鋼進行了溫度測量,如圖3a和3b所示,在55NiCrMoV7模具鋼的外部至內(nèi)部均勻選取了11個測溫點,溫度數(shù)據(jù)每間隔7 min提取一次,選?、俸?位置的提取溫度與同步模擬溫度數(shù)據(jù)進行對比,結(jié)果如圖4所示,可以看到在圖4a中,①位置(即表層溫度)模擬表面溫度曲線與實測溫度曲線基本重合,這證明了選取的網(wǎng)格劃分方式以及熱交換方式的正確性。在?位置(即心部溫度)差別最大,如圖4b所示,這是由鍛件體積過于龐大以及在實際測溫中不可避免的操作誤差導致,但這些誤差都在合理的范圍之內(nèi),綜合考量該模型具有一定的可靠性。

    圖3 不同部位溫度測量

    圖4 模擬溫度與實際溫度對比

    2.2.2 微觀組織

    55NiCrMoV7模具鋼同步調(diào)質(zhì)熱處理模擬結(jié)果如圖5所示,在熱處理過程中棒料邊角處吸熱、散熱能力較強,溫度變化較快,故馬氏體組織集中在棒料兩端的邊角處,如圖5c所示,體積分數(shù)最大為23.9%。棒料內(nèi)部在高溫回火的條件下,滲碳體以顆粒狀分布于鐵素體的基體上,呈現(xiàn)粒狀珠光體組織,由圖5a可以看出珠光體在體積分數(shù)最小處都達到了52.1%。由圖5b可以看出貝氏體主要分布在棒料的外層,體積分數(shù)占到了30.7%。綜上分析可知,整個棒料的組織成分分布與溫度場變化表現(xiàn)一致,這進一步印證了模型的可靠性。

    圖5 模擬調(diào)質(zhì)55NiCrMoV7鋼組織分布

    對調(diào)質(zhì)熱處理后的55NiCrMoV7模具鋼棒料取樣進行金相組織觀察,如圖6所示,1號取樣位置距端部50 mm,2號取樣位置距端部425 mm(即1/4處),3號取樣位置為棒料正中間,4號取樣位置為棒料另一端且距端部50 mm處,所有取樣寬度均為20 mm。每個取樣位置由表及里分別取樣并編號為:1-1,1-2,1-3,2-1,2-2,2-3,3-1,3-2,3-3,4-1,4-2,4-3,取樣位置詳見圖6。

    55NiCrMoV7鋼棒料表層組織成分變化如圖7a,d,g所示。由圖7a可以看出,基體以珠光體為主,針狀馬氏體組織附著在其表面,極少量的粒狀碳化物存在于馬氏體之間,大量碳化物都溶入了馬氏體之中,且馬氏體組織明顯多于圖7d和圖7g。在圖7d中貝氏體組織顯現(xiàn),與馬氏體組織夾雜分布,組織間碳化物增多且分布雜亂。圖7g中幾乎看不到馬氏體的蹤跡,組織以珠光體和貝氏體為主,粒狀碳化物較圖7d略小。隨著向棒料內(nèi)部推進組織成分的變化,如圖7b,e,h所示,從圖7b可以看出,部分碳化物在馬氏體中析出,馬氏體組織含量減少。圖7e中馬氏體消失,碳化物不斷長大。在圖7h中碳化物體積略有減小,且趨于均勻化分布。棒料心部組織變化如圖7c,f,i所示,可以看到心部共同特點是碳化物顆粒比較細小且分布均勻,其中圖7i棒料中心部位,碳化物顆粒細小而彌散均勻分布于鐵素體基體上,呈現(xiàn)為粒狀珠光體組織,該部位組織均勻性最好。圖7j—l取樣位置與圖7a—c取樣位置呈對稱分布,觀察其組織分布可以發(fā)現(xiàn)二者一致,故證明了該實驗結(jié)果的正確性。模擬組織分布與實驗結(jié)果觀察一致,也證明了大鍛件55NiCrMoV7鋼有限元模型的可靠性。

    圖6 金相取樣部位示意

    3 結(jié)果與討論

    3.1 淬火溫度影響

    3.1.1 殘余應(yīng)力變化

    790~910 ℃不同溫度淬火+510 ℃回火模擬后的殘余應(yīng)力分布情況如圖8所示,從圖8a—c可以看出,殘余應(yīng)力最小區(qū)域主要集中在棒料中心兩端且成對稱分布,棒料側(cè)面應(yīng)力較小區(qū)域隨著淬火溫度的增加及淬透性的增強,而逐漸消失,這符合熱應(yīng)力產(chǎn)生規(guī)律[12]。從圖8d可知,隨著淬火溫度的升高,殘余應(yīng)力最大值略微提高,且不超過800 MPa,最小值略有浮動。綜上所述,在模擬過程中隨著淬火溫度的上升,雖然殘余應(yīng)力最小區(qū)域不斷縮小但高低應(yīng)力區(qū)域劃分明顯,故淬火溫度對改善殘余應(yīng)力分布影響較小。

    圖7 55NiCrMoV7鋼調(diào)質(zhì)后不同部位金相組織

    圖8 不同溫度淬火+510 ℃回火應(yīng)力分布

    3.1.2 模擬硬度和屈服強度變化

    由圖9a可知,790 ℃淬火時硬度僅為42.8HRC,但在820 ℃及以上淬火后,硬度均大于44HRC,其中820 ℃淬火硬度最大為45.1HRC,在860 ℃淬火后實驗測得硬度為44.7HRC,與預(yù)測值接近,故不同淬火溫度對硬度的影響很小。屈服強度隨淬火溫度變化如圖9b所示,在淬火溫度為880 ℃時,屈服強度最大為1335 MPa,860 ℃淬火實驗值為1310 MPa。綜上所述,為了兼顧材料硬度和屈服強度,淬火溫度應(yīng)保持在850 ℃以上。

    3.2 回火溫度影響

    3.2.1 殘余應(yīng)力變化

    850 ℃淬火+不同溫度回火模擬后的殘余應(yīng)力分布如圖10所示,由圖10a—c可以看出,殘余應(yīng)力最大區(qū)域總是集中在棒料表層,這是由于表層溫度變化較快,塑性變形不均勻而產(chǎn)生熱應(yīng)力,殘余應(yīng)力最小區(qū)域呈杯狀分布在棒料兩端及側(cè)面,且隨著回火溫度的升高側(cè)面低應(yīng)力區(qū)域逐漸發(fā)散最后消失,端部杯狀區(qū)域也有所減小,高、低應(yīng)力區(qū)劃分愈發(fā)明顯。其中380 ℃回火后,如圖10b所示,低應(yīng)力區(qū)域發(fā)散明顯,高低應(yīng)力區(qū)域界線模糊,應(yīng)力均勻性較好。由圖10d可以看出,隨著回火溫度的升高,殘余應(yīng)力最大值不斷降低,最低可達476 MPa。綜上所述,回火溫度對殘余應(yīng)力分布影響較大,其中380 ℃回火可以有效提高大型鍛件殘余應(yīng)力分布均勻性。

    3.2.2 模擬硬度和屈服強度變化

    由上述可知,回火溫度的升高有利于消除殘余應(yīng)力,提高應(yīng)力分布均勻性,但是過高的回火溫度會造成其他力學性能下降,如圖11所示,由圖11a可知,當回火溫度高于510 ℃時,材料的硬度會急劇下降,這是因為500 ℃以上高溫回火,碳化物不斷在馬氏體中析出,導致馬氏體的含碳量大大降低,硬度呈現(xiàn)斷崖式下跌[13]。如圖11b所示,屈服強度在380 ℃時達到最大值1450 MPa,此后回火溫度越高,屈服強度下降越明顯,這與回火過程中馬氏體組織軟化有關(guān)。在560 ℃回火實驗后測得硬度為41.8HRC,屈服強度為1130 MPa,預(yù)測值與實驗值在誤差允許的范圍之內(nèi)。

    圖9 淬火溫度對硬度和屈服強度的影響

    圖10 850 ℃淬火不同溫度回火應(yīng)力分布

    圖11 回火溫度對硬度和屈服強度的影響

    3.3 熱處理工藝優(yōu)化

    3.3.1 響應(yīng)面模型的建立

    基于調(diào)質(zhì)模擬實驗設(shè)計及結(jié)果(如表2所示)對調(diào)質(zhì)熱處理工藝進行殘余應(yīng)力分布、硬度和屈服強度的多目標優(yōu)化。實驗采用二因素、五水平中心組合[14](Central Composite Design)的設(shè)計方法,獲得了調(diào)質(zhì)工藝與力學性能的響應(yīng)關(guān)系,即將淬火溫度和回火溫度作為自變量,應(yīng)力均方差洛氏硬度和屈服強度s作為因變量,根據(jù)表2中模擬數(shù)據(jù)構(gòu)建函數(shù)模型,經(jīng)擬合得到回歸方程如下:

    (2)

    式中:1為淬火溫度(℃);2為回火溫度(℃)。

    3.3.2 響應(yīng)面模型分析

    根據(jù)調(diào)質(zhì)模擬實驗結(jié)果繪制淬火溫度和回火溫度對應(yīng)力均方差、硬度和屈服強度的3D響應(yīng)曲面模型,如圖12所示。在400~550 ℃區(qū)間回火時,850~ 910 ℃的淬火溫度會使應(yīng)力均方差略低,即鍛件應(yīng)力分布均勻性較好。同時由圖13—14可知,在850~ 910 ℃隨著淬火溫度的升高,硬度和屈服強度的值也在增大,因此淬火溫度應(yīng)控制在850~910 ℃。由圖12可以看出400 ℃以下回火,會導致殘余應(yīng)力均方差的值大大增加,即大型鍛件的應(yīng)力分布均勻性變差,隨著回火溫度的升高,應(yīng)力均方差值逐漸減小,有利于提高殘余應(yīng)力分布的均勻性,因此應(yīng)控制回火溫度不低于400 ℃。但是回火溫度太高會導致材料的硬度和屈服強度大幅降低,如圖13所示,超過550 ℃的高溫回火,隨著回火溫度的升高,會使材料的硬度急劇下降,同時對比圖14可知,屈服強度的變化趨勢與硬度變化基本一致。因此,為了保證各項指標的合格性,回火溫度應(yīng)控制在400~550 ℃。

    表2 調(diào)質(zhì)模擬實驗設(shè)計及結(jié)果

    Tab.2 Experiment design and results of simulated tempering and quenching

    3.3.3 基于遺傳算法的多目標優(yōu)化

    3D響應(yīng)曲面模型可以直觀反映出調(diào)質(zhì)工藝對某一性能指標的影響規(guī)律,但不能對應(yīng)力均方差、硬度和屈服強度3個目標協(xié)調(diào)求解。為了獲得匹配性較好的調(diào)質(zhì)工藝,文中基于Matlab遺傳算法工具箱(GA工具箱)進一步對大型鍛件55NiCrMoV7模具鋼的力

    圖12 3D響應(yīng)曲面模型

    學性能進行多目標優(yōu)化。在多目標的優(yōu)化過程中,所有的目標函數(shù)必須同時被最大或最小化,GA工具箱默認函數(shù)值越小越優(yōu)[15]。故多目標函數(shù)可表示為:

    式中:1()為應(yīng)力均方差模型;2()為硬度模型;3()為屈服強度模型s。根據(jù)式(4)建立優(yōu)化模型,設(shè)置初始種群為50,最優(yōu)個體系數(shù)為0.3,進化代數(shù)為200,交叉概率為0.95,變異概率為0.05,運行結(jié)果如表3所示。

    表3 多目標優(yōu)化模型運行結(jié)果

    Tab.3 Running results of multi-objective optimization model

    通過多目標優(yōu)化模型運行結(jié)果可知,并不存在同時使應(yīng)力均方差最小、硬度和屈服強度最大的解,但是可以找到符合各項指標要求的解。要求殘余應(yīng)力均方差不大于60,硬度不低于45HRC,屈服強度不小于1300 MPa,如圖13所示,最終符合綜合力學性能的點為:1,2,3,8,13。對應(yīng)的工藝參數(shù)分別為:淬火907.4 ℃+回火443.1 ℃、淬火907.4 ℃+回火411.5 ℃、淬火903.4 ℃+回火485.1 ℃、淬火895.2 ℃+回火535.8 ℃和淬火892.1 ℃+回火507.1 ℃。在Pareto最優(yōu)解求解的過程中,每次計算過程進行的交叉進化和變異不同,導致Pareto最優(yōu)解求解運行的結(jié)果不同,但每次結(jié)果的改變都是通過犧牲其他目標函數(shù)實現(xiàn)的,每個Pareto最優(yōu)解都同時滿足多目標優(yōu)化問題[16]。這就意味著以上5組工藝屬于平行關(guān)系,無優(yōu)劣之分。

    圖13 多目標優(yōu)化Pareto最優(yōu)解

    4 結(jié)語

    1)建立了大型55NiCrMoV7模具鋼鍛件的有限元模型,并通過實驗測溫和微觀組織分析驗證了該模型。

    2)建立了淬、回火溫度對應(yīng)力均方差、HRC硬度和屈服強度的回歸模型。

    3)確定最優(yōu)工藝范圍:淬火溫度為880~910 ℃,回火溫度為400~550 ℃。得到Pareto最優(yōu)解集:淬火907.4 ℃+回火443.1 ℃、淬火907.4 ℃+回火411.5 ℃、淬火903.4 ℃+回火485.1 ℃、淬火895.2 ℃+回火535.8 ℃和淬火892.1 ℃+回火507.1 ℃。

    [1] 史正良, 劉金龍. 熱作模具鋼的選擇與應(yīng)用[J]. 煤礦機械, 2007(12): 126-129.

    SHI Zheng-liang, LIU Jin-long. Choice and Application of Hot Die Steels[J]. Coal Mine Machinery, 2007(12): 126-129.

    [2] JEBARAJ M, PRADEEP K M, ANBURAJ R. Effect of LN2and CO2Coolants in Milling of 55NiCrMoV7 Steel[J]. Journal of Manufacturing Processes, 2020, 53: 318-327.

    [3] JEBARAJ M, PRADEEP K M, YUVARAJ N, et al. Investigation of Surface Integrity in End Milling of 55NiCrMoV7 Die Steel under the Cryogenic Environments[J]. Machining Science and Technology, 2020, 24(3): 465-488.

    [4] HOJA S, SKALECKI M G, KLüMPER-WESTKAMP H, et al. Carbonitrieren von Warmarbeitsst?hlen[J]. HTM- Journal of Heat Treatment and Materials, 2017, 72(4): 187-198.

    [5] ZHANG W, WANG H, ZHANG J, et al. Brittle Fracture Behaviors of Large Die Holders Used in Hot Die Forging[J]. Metals, 2017, 7(6): 198.

    [6] 張占平, 齊育紅, DELAGNES D, 等. 鋼的回火時間-溫度-硬度動力學關(guān)系[J]. 金屬熱處理學報, 2004(1): 41-45.

    ZHANG Zhan-ping, QI Yu-hong, DELAGNES D, et al. Time-Temperature-Hardness Kinetics Relation of Tempering Steels[J]. Transactions of Materials and Heat Treatment, 2004(1): 41-45.

    [7] 張占平, 齊育紅, DELAGNES D, 等. 55NiCrMoV7鋼的回火組織與性能[J]. 金屬熱處理, 2004(3): 27-31.

    ZHANG Zhan-ping, QI Yu-hong, DELAGNES D, et al. Microstructures and Properties of the Tempered Steel 55NiCrMoV7[J]. Heat Treatment of Metals, 2004(3): 27-31.

    [8] FAN Meng-ting, SUN Ming-yue, LI Dian-zhong. Effect of Microstructural Evolution on Mechanical Properties of 55NiCrMoV7 Steel[J]. Advanced Materials Research, 2011, 1334: 848-852.

    [9] 樊夢婷, 孫明月, 李殿中. 大型壓力機模座熱處理過程模擬及工藝優(yōu)化[J]. 材料工程, 2011(11): 44-50.

    FAN Meng-ting, SUN Ming-yue, LI Dian-zhong. Simulation and Optimization of Heat Treatment Process for Large Press Die Holder[J]. Journal of Materials Engineering, 2011(11): 44-50.

    [10] 禹興勝, 武川, 李昌義, 等. 55NiCrMoV7鋼錠均質(zhì)化處理宏觀偏析熱擴散元胞自動機模擬研究[J]. 科技創(chuàng)新與應(yīng)用, 2020(31): 18-19.

    YU Xing-sheng, WU Chuan, LI Chang-yi, et al. Simulation of Macrosegregation Thermal Diffusion Cellular Automata for Homogenization of 55NiCrMoV7 Ingot[J]. Technology Innovation and Application, 2020(31): 18- 19.

    [11] 趙繼全, 曲周德, 禹興勝, 等. 調(diào)質(zhì)熱處理對55NiCr-MoV7模具鋼組織與力學性能的影響[J]. 特殊鋼, 2021, 42(4): 60-65.

    ZHAO Ji-quan, QU Zhou-de, YU Xing-sheng, et al. Micro-structure Evolution and Mechanical Properties Study of 55NiCrMoV7 Die Steel after Quenching and Tempering[J]. Special Steel, 2021, 42(4): 60-65.

    [12] WU Chuan, MENG Qing-ling. Microstructural Evolution of a Steam Turbine Rotor Subjected to a Water Quenching Process: Numerical Simulation and Experimental Verification[J]. Advances in Manufacturing, 2019, 7(1): 84-104.

    [13] 李爽, 時彥林, 楊曉彩, 等. 鉬鎢釩合金化熱作模具鋼高溫回火組織演變[J]. 工程科學學報, 2020, 42(7): 902-911.

    LI Shuang, SHI Yan-lin, YANG Xiao-cai, et al. Mic-ros-tr-uc-tural Evolution of Mo-W-V Alloyed Hot-Work Die Steel during High-Temperature Temp-ering[J]. Chinese Journal of Engineering, 2020, 42(7): 902-911.

    [14] 肖瑞, 楊明, 黃朝文. 基于響應(yīng)面法對211ZX新型高強鋁合金固溶工藝的優(yōu)化設(shè)計[J]. 稀有金屬, 2019, 43(10): 1040-1046.

    XIAO Rui, YANG Ming, HUANG Chao-wen. Optimal Design of 211ZX High Strength Aluminum Alloy Solid Solution Process Based on Response Surface Method[J]. Chinese Journal of Rare Metals, 2019, 43(10): 1040- 1046.

    [15] PHAN-NGUYEN Huu. Multi-Objective Optimization in Titanium Powder Mixed Electrical Discharge Machining Process Parameters for Die Steels[J]. Alexandria Engineering Journal, 2020, 59(6): 4063-4079.

    [16] 盧秋霞, 劉潔, 王志鵬. 鋁合金板件自沖鉚接性能的響應(yīng)面-遺傳算法優(yōu)化[J]. 鍛壓技術(shù), 2021(5): 57-63.

    LU Qiu-xia, LIU Jie, WANG Zhi-peng. Optimization on Self-Piercing Riveted Performance of Aluminum Alloy Plate Part Based on Response Surface and Genetic Algorithm[J]. Forging & Stamping Technology, 2021(5): 57-63.

    Optimization of Quenching and Tempering Heat Treatment Process for 55NiCrMoV7 Die Steel Based on Multi-Scale Modeling and Response Surface Methodology

    ZHAO Ji-quan1, QU Zhou-de1, WU Chuan1,2,3, SHI Ru-xing2,3, YU Xing-sheng2,3

    (1. National-local Joint Engineering Laboratory of Intelligent Manufacturing Oriented Automobile Die & Mould, Tianjin University of Technology and Education, Tianjin 300222, China; 2. Luoyang CITIC Casting and Forging Co., Ltd., Luoyang 471039, China; 3. CITIC Heavy Industries Co., Ltd., Luoyang 471039, China)

    The work aims to investigate the optimal process scheme of quenching and tempering heat treatment of large 55NiCrMoV7 die steel forgings (1000 mm × 1700 mm). A multi-scale finite element model of macro-micro coupling for 55NiCrMoV7 die steel forgings was established to test the temperature field and microstructure field evolution at different locations based on field experiments and then the model was verified. The stress distribution, HRC hardness and yield strength of the forgings were simulated under different quenching and tempering processes. The regression model between the parameters of quenching and tempering and mechanical properties was established with the mean square error of residual stress, hardness and yield strength as the objective function. Combined with the response surface model analysis and genetic algorithm, the multi-objective optimization of the model under parameter constraints was carried out. According to the analysis of response surface model, the optimal process range was as follows: quenching temperature of 850~910 ℃, tempering temperature of 400~550 ℃. After multi-objective optimization, the optimal solution set of Pareto was obtained: quenching of 907.4 ℃+ tempering of 443.1 ℃, quenching of 907.4 ℃+ tempering of 411.5 ℃, quenching of 903.4 ℃+ tempering of 485.1 ℃, quenching of 895.2 ℃+ tempering of 535.8 ℃ and quenching of 892.1 ℃+ tempering of 507.1 ℃. The established regression model of quenching and tempering temperature on stress mean square error, HRC hardness and yield strength provides a theoretical basis for the optimization of quenching and tempering heat treatment process for large 55NiCrMoV7 die steel forgings.

    55NiCrMoV7 die steel; multi-objective optimization; response surface method; genetic algorithm; multi-scale model

    10.3969/j.issn.1674-6457.2022.02.012

    TG316

    A

    1674-6457(2022)02-0074-09

    2021-05-20

    中國博士后第67批面上項目(2020M672309)

    趙繼全(1996—),男,碩士生,主要研究方向為模具鋼熱處理。

    曲周德(1972—),男,博士,教授,主要研究方向為模具鋼熱處理。

    猜你喜歡
    棒料模具鋼調(diào)質(zhì)
    CaO調(diào)質(zhì)污泥自然失水特性及數(shù)學模型研究
    LN17模具鋼預(yù)硬化熱處理正交試驗
    模具制造(2019年7期)2019-09-25 07:30:06
    基于隨機采樣一致性的規(guī)則棒料位姿識別*
    調(diào)質(zhì)型管線鋼關(guān)鍵工藝技術(shù)研究
    山東冶金(2019年3期)2019-07-10 00:54:00
    車床自動送料機構(gòu)的棒料振動變形分析
    重型機械(2018年6期)2019-01-07 11:13:48
    預(yù)硬型塑料模具鋼1.2312生產(chǎn)實踐
    四川冶金(2018年1期)2018-09-25 02:39:28
    直讀光譜法測定熱作模具鋼中硫的不確定度評定
    Cr12MoV冷作模具鋼滲鉻層界面能譜分析
    水產(chǎn)膨化常見的幾種調(diào)質(zhì)器的結(jié)構(gòu)特征與性能比較
    酸調(diào)質(zhì)—氧化降解法提高含油污泥的沉降性能
    一边摸一边抽搐一进一小说| 国产精品爽爽va在线观看网站| 日本黄大片高清| 一二三四在线观看免费中文在| 夜夜看夜夜爽夜夜摸| 动漫黄色视频在线观看| 午夜福利高清视频| 亚洲最大成人中文| 久久精品91蜜桃| 精品久久久久久,| 五月伊人婷婷丁香| 99re在线观看精品视频| 变态另类丝袜制服| 在线观看美女被高潮喷水网站 | 丰满人妻一区二区三区视频av | 正在播放国产对白刺激| 国内久久婷婷六月综合欲色啪| 亚洲 欧美一区二区三区| 男男h啪啪无遮挡| 看片在线看免费视频| 一本综合久久免费| 国产高清视频在线播放一区| 午夜免费成人在线视频| 精品乱码久久久久久99久播| av国产免费在线观看| 国产不卡一卡二| 久久这里只有精品中国| 久久精品综合一区二区三区| 亚洲九九香蕉| 十八禁人妻一区二区| 亚洲人成伊人成综合网2020| 成年女人毛片免费观看观看9| 一卡2卡三卡四卡精品乱码亚洲| 久久久久国产一级毛片高清牌| 久久久久精品国产欧美久久久| 可以在线观看毛片的网站| 国产aⅴ精品一区二区三区波| 久久 成人 亚洲| 日日摸夜夜添夜夜添小说| 亚洲中文日韩欧美视频| 国产午夜精品论理片| 欧美不卡视频在线免费观看 | videosex国产| 国产精品,欧美在线| 成在线人永久免费视频| 制服丝袜大香蕉在线| 日韩av在线大香蕉| 日韩高清综合在线| 禁无遮挡网站| 18禁黄网站禁片午夜丰满| 午夜福利成人在线免费观看| 两个人看的免费小视频| 亚洲av中文字字幕乱码综合| 亚洲最大成人中文| 午夜两性在线视频| 哪里可以看免费的av片| 最新在线观看一区二区三区| 午夜激情av网站| 亚洲精品久久成人aⅴ小说| 国产高清有码在线观看视频 | 亚洲av成人不卡在线观看播放网| 久久精品综合一区二区三区| 中文字幕精品亚洲无线码一区| 成人一区二区视频在线观看| 国产男靠女视频免费网站| 久久中文字幕一级| 欧美极品一区二区三区四区| 亚洲熟女毛片儿| 两个人免费观看高清视频| 久久精品国产亚洲av高清一级| 午夜福利成人在线免费观看| 一二三四在线观看免费中文在| 丰满的人妻完整版| 最近在线观看免费完整版| 国内精品久久久久精免费| 成人av一区二区三区在线看| 欧美日韩中文字幕国产精品一区二区三区| 国产成人一区二区三区免费视频网站| 欧美 亚洲 国产 日韩一| 在线视频色国产色| 午夜福利在线观看吧| 久久久久精品国产欧美久久久| 亚洲熟女毛片儿| 天堂动漫精品| 18禁黄网站禁片午夜丰满| 国产成人精品久久二区二区免费| 18禁黄网站禁片免费观看直播| 狠狠狠狠99中文字幕| 成熟少妇高潮喷水视频| 国产亚洲精品久久久久久毛片| 国产精品一及| 亚洲成人久久爱视频| 亚洲成av人片在线播放无| 精品久久久久久,| 亚洲av熟女| 日本在线视频免费播放| 国产高清视频在线播放一区| 国产黄色小视频在线观看| 亚洲avbb在线观看| 国产亚洲精品久久久久久毛片| 国产成年人精品一区二区| 午夜精品在线福利| 日日摸夜夜添夜夜添小说| av天堂在线播放| 在线观看美女被高潮喷水网站 | 国产野战对白在线观看| 51午夜福利影视在线观看| 亚洲av电影不卡..在线观看| 小说图片视频综合网站| 久久香蕉国产精品| www国产在线视频色| 成人高潮视频无遮挡免费网站| av在线播放免费不卡| 麻豆av在线久日| 九色国产91popny在线| 悠悠久久av| 久久久久久国产a免费观看| av片东京热男人的天堂| 曰老女人黄片| 免费电影在线观看免费观看| 99热6这里只有精品| 欧美日本视频| 国产91精品成人一区二区三区| 香蕉国产在线看| 一区二区三区国产精品乱码| 99re在线观看精品视频| 成人手机av| 午夜久久久久精精品| 久久精品国产综合久久久| 欧美日韩黄片免| 亚洲成人中文字幕在线播放| 亚洲人成77777在线视频| 两性夫妻黄色片| 国产亚洲欧美98| 国产免费男女视频| 丰满的人妻完整版| 97超级碰碰碰精品色视频在线观看| 欧美绝顶高潮抽搐喷水| 免费在线观看完整版高清| 可以在线观看的亚洲视频| 午夜日韩欧美国产| 五月伊人婷婷丁香| 久99久视频精品免费| 又黄又粗又硬又大视频| av天堂在线播放| 久久久久久国产a免费观看| 操出白浆在线播放| 白带黄色成豆腐渣| 久久久久久大精品| 亚洲天堂国产精品一区在线| 人人妻,人人澡人人爽秒播| 最近最新中文字幕大全免费视频| 精品国产乱子伦一区二区三区| 日韩精品青青久久久久久| 国产又黄又爽又无遮挡在线| 久久 成人 亚洲| 日韩av在线大香蕉| 亚洲成a人片在线一区二区| 成人一区二区视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产爱豆传媒在线观看 | 免费看十八禁软件| 欧美一级毛片孕妇| 久久久久久久久中文| 亚洲一区二区三区色噜噜| 男女做爰动态图高潮gif福利片| 精品免费久久久久久久清纯| 99热6这里只有精品| 99精品在免费线老司机午夜| 国产亚洲av高清不卡| 国产精品香港三级国产av潘金莲| 极品教师在线免费播放| 亚洲成人久久爱视频| 怎么达到女性高潮| 妹子高潮喷水视频| 嫩草影院精品99| 日本三级黄在线观看| 久久久久久人人人人人| 亚洲精品中文字幕一二三四区| 观看免费一级毛片| www日本黄色视频网| 午夜福利视频1000在线观看| 亚洲午夜理论影院| 在线永久观看黄色视频| 久久久久久亚洲精品国产蜜桃av| 12—13女人毛片做爰片一| 日韩欧美免费精品| 国产99久久九九免费精品| 国产成人系列免费观看| 操出白浆在线播放| 国产精品久久久久久久电影 | 亚洲熟妇中文字幕五十中出| 欧美中文日本在线观看视频| 日韩免费av在线播放| 亚洲av熟女| 免费观看精品视频网站| 国产精品香港三级国产av潘金莲| 一级作爱视频免费观看| 黑人巨大精品欧美一区二区mp4| 国产三级黄色录像| 中文字幕最新亚洲高清| 成人三级黄色视频| 黄色片一级片一级黄色片| 欧美激情久久久久久爽电影| 亚洲熟妇熟女久久| 高清在线国产一区| 在线视频色国产色| 国产成人aa在线观看| 亚洲国产精品合色在线| 麻豆国产97在线/欧美 | 中文字幕人成人乱码亚洲影| 精品无人区乱码1区二区| 黄色成人免费大全| 制服丝袜大香蕉在线| 亚洲av成人精品一区久久| 久久久国产欧美日韩av| av福利片在线| 熟女电影av网| cao死你这个sao货| 久久精品亚洲精品国产色婷小说| 日韩三级视频一区二区三区| 91大片在线观看| 无遮挡黄片免费观看| 国产亚洲av嫩草精品影院| 两性午夜刺激爽爽歪歪视频在线观看 | 成人午夜高清在线视频| 看黄色毛片网站| 免费av毛片视频| 久久久久性生活片| 12—13女人毛片做爰片一| 人成视频在线观看免费观看| 韩国av一区二区三区四区| 免费电影在线观看免费观看| 国产成人啪精品午夜网站| 波多野结衣高清无吗| 日韩高清综合在线| 亚洲av五月六月丁香网| 在线观看美女被高潮喷水网站 | 91麻豆精品激情在线观看国产| 亚洲欧美激情综合另类| 日韩免费av在线播放| 18美女黄网站色大片免费观看| 日韩欧美在线乱码| 亚洲熟妇熟女久久| 日本一本二区三区精品| 18禁黄网站禁片免费观看直播| 国产真实乱freesex| 精品福利观看| 最近最新中文字幕大全免费视频| 一本综合久久免费| or卡值多少钱| 日日夜夜操网爽| 亚洲精品久久成人aⅴ小说| 精品一区二区三区四区五区乱码| 国产高清视频在线播放一区| 成年人黄色毛片网站| 国产精品国产高清国产av| 亚洲熟妇中文字幕五十中出| 日本在线视频免费播放| 国产精品香港三级国产av潘金莲| 日韩有码中文字幕| 无限看片的www在线观看| 国产野战对白在线观看| 精品免费久久久久久久清纯| 成人亚洲精品av一区二区| 国产爱豆传媒在线观看 | 午夜福利高清视频| 99在线视频只有这里精品首页| 色老头精品视频在线观看| 久久久久久久午夜电影| 国产探花在线观看一区二区| 色综合亚洲欧美另类图片| 国产单亲对白刺激| 亚洲乱码一区二区免费版| 精品欧美国产一区二区三| av免费在线观看网站| av在线天堂中文字幕| 免费观看精品视频网站| 国产高清视频在线观看网站| 日本一本二区三区精品| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲全国av大片| 国产成人影院久久av| 欧美日韩瑟瑟在线播放| 性色av乱码一区二区三区2| 亚洲中文日韩欧美视频| 亚洲中文av在线| 黑人欧美特级aaaaaa片| 免费av毛片视频| 老司机福利观看| 三级国产精品欧美在线观看 | 亚洲欧洲精品一区二区精品久久久| 少妇人妻一区二区三区视频| 久久久精品国产亚洲av高清涩受| 老汉色∧v一级毛片| a级毛片在线看网站| 搞女人的毛片| 黄频高清免费视频| 免费在线观看视频国产中文字幕亚洲| 亚洲人成77777在线视频| 国产又黄又爽又无遮挡在线| 一区二区三区高清视频在线| 露出奶头的视频| 国产激情久久老熟女| 男女午夜视频在线观看| 两人在一起打扑克的视频| 一级片免费观看大全| 韩国av一区二区三区四区| 亚洲精品国产精品久久久不卡| 天堂影院成人在线观看| 国产精品久久久人人做人人爽| 精品久久蜜臀av无| 日韩成人在线观看一区二区三区| 欧美黄色淫秽网站| 欧美日本亚洲视频在线播放| 免费在线观看成人毛片| 亚洲全国av大片| 好男人在线观看高清免费视频| 黑人巨大精品欧美一区二区mp4| 久久久久久久精品吃奶| 国产免费男女视频| 我要搜黄色片| 欧美 亚洲 国产 日韩一| 欧美性猛交╳xxx乱大交人| 黑人欧美特级aaaaaa片| 又粗又爽又猛毛片免费看| 国产av不卡久久| 亚洲精品色激情综合| 国产69精品久久久久777片 | 久久伊人香网站| 欧美一区二区精品小视频在线| 91av网站免费观看| 国产精品久久电影中文字幕| 久久久久久九九精品二区国产 | 老司机在亚洲福利影院| 久久婷婷人人爽人人干人人爱| 国模一区二区三区四区视频 | 岛国在线观看网站| 一进一出好大好爽视频| 一本久久中文字幕| 国产99久久九九免费精品| 免费看a级黄色片| 狠狠狠狠99中文字幕| 黄色a级毛片大全视频| 男人舔女人的私密视频| 男女视频在线观看网站免费 | 国产一区二区在线av高清观看| 国产97色在线日韩免费| 黑人巨大精品欧美一区二区mp4| 琪琪午夜伦伦电影理论片6080| 色综合欧美亚洲国产小说| 悠悠久久av| 精品久久久久久久人妻蜜臀av| 日本 欧美在线| 欧美性长视频在线观看| 极品教师在线免费播放| 特级一级黄色大片| 亚洲七黄色美女视频| 国产免费男女视频| 成人午夜高清在线视频| 欧美精品亚洲一区二区| 又大又爽又粗| 精品不卡国产一区二区三区| 青草久久国产| 国产高清视频在线播放一区| 中文字幕精品亚洲无线码一区| 亚洲 国产 在线| 黄色片一级片一级黄色片| 日韩免费av在线播放| 欧美黑人巨大hd| 国产亚洲精品综合一区在线观看 | 日本三级黄在线观看| 国产av在哪里看| 禁无遮挡网站| 欧美日韩瑟瑟在线播放| 五月玫瑰六月丁香| 亚洲第一欧美日韩一区二区三区| 亚洲国产高清在线一区二区三| 亚洲av电影在线进入| svipshipincom国产片| 亚洲精品在线观看二区| 久久久久九九精品影院| 女人爽到高潮嗷嗷叫在线视频| 国产成年人精品一区二区| 特大巨黑吊av在线直播| 久久久久久九九精品二区国产 | 国产精品1区2区在线观看.| 国产亚洲欧美在线一区二区| 国产成人啪精品午夜网站| 久久久久亚洲av毛片大全| 丝袜人妻中文字幕| 成人永久免费在线观看视频| 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| www.自偷自拍.com| 两性夫妻黄色片| 精品久久久久久久人妻蜜臀av| 欧美一区二区国产精品久久精品 | 88av欧美| 久久 成人 亚洲| a级毛片a级免费在线| 在线观看舔阴道视频| 免费看美女性在线毛片视频| 亚洲一区中文字幕在线| 麻豆av在线久日| 亚洲 欧美 日韩 在线 免费| 久9热在线精品视频| 999久久久国产精品视频| 很黄的视频免费| 日韩欧美精品v在线| 人妻久久中文字幕网| 国产精品av视频在线免费观看| 正在播放国产对白刺激| 天堂√8在线中文| 麻豆一二三区av精品| 亚洲人成77777在线视频| 国产伦一二天堂av在线观看| 夜夜爽天天搞| 不卡一级毛片| 亚洲欧美激情综合另类| 亚洲国产看品久久| 91大片在线观看| 免费高清视频大片| 国产亚洲精品av在线| 欧美日韩亚洲综合一区二区三区_| 老司机深夜福利视频在线观看| cao死你这个sao货| 精品国产超薄肉色丝袜足j| 欧美日韩一级在线毛片| 午夜影院日韩av| 精品欧美一区二区三区在线| 黄色女人牲交| 免费在线观看日本一区| 亚洲九九香蕉| 亚洲精品一卡2卡三卡4卡5卡| 一本久久中文字幕| 又粗又爽又猛毛片免费看| 国产69精品久久久久777片 | 免费看十八禁软件| 一进一出抽搐gif免费好疼| 三级毛片av免费| 久久人人精品亚洲av| 欧美成人免费av一区二区三区| 国产精品久久久人人做人人爽| 亚洲av成人精品一区久久| 他把我摸到了高潮在线观看| 亚洲avbb在线观看| 男女那种视频在线观看| 亚洲熟妇熟女久久| 国产精品亚洲美女久久久| cao死你这个sao货| 久久久久久免费高清国产稀缺| 日韩有码中文字幕| 18禁黄网站禁片免费观看直播| 成人18禁在线播放| 亚洲熟妇熟女久久| 国产免费男女视频| a在线观看视频网站| 午夜福利高清视频| 婷婷丁香在线五月| 美女高潮喷水抽搐中文字幕| 成人精品一区二区免费| 99国产精品一区二区蜜桃av| 一区二区三区激情视频| 亚洲欧美激情综合另类| 草草在线视频免费看| 高潮久久久久久久久久久不卡| 国产精品爽爽va在线观看网站| 亚洲欧美日韩东京热| 18美女黄网站色大片免费观看| 两个人的视频大全免费| 亚洲免费av在线视频| 久久精品国产亚洲av香蕉五月| 国产精品一区二区免费欧美| 天天添夜夜摸| 欧美日韩精品网址| 一级作爱视频免费观看| 亚洲成人中文字幕在线播放| 亚洲人成网站高清观看| 丝袜美腿诱惑在线| 女同久久另类99精品国产91| 国产麻豆成人av免费视频| 啦啦啦免费观看视频1| 黄色丝袜av网址大全| 亚洲精品国产精品久久久不卡| 国产免费男女视频| 午夜福利视频1000在线观看| 日韩精品青青久久久久久| 欧美极品一区二区三区四区| 一二三四社区在线视频社区8| 女人爽到高潮嗷嗷叫在线视频| a级毛片在线看网站| 亚洲专区国产一区二区| 亚洲五月天丁香| 在线观看66精品国产| 在线观看午夜福利视频| 午夜成年电影在线免费观看| 久久国产精品影院| 黄频高清免费视频| 特级一级黄色大片| 51午夜福利影视在线观看| 真人做人爱边吃奶动态| 天堂√8在线中文| 精品一区二区三区av网在线观看| 午夜a级毛片| 成年女人毛片免费观看观看9| 亚洲国产精品999在线| 在线看三级毛片| 99国产极品粉嫩在线观看| 中文字幕精品亚洲无线码一区| 每晚都被弄得嗷嗷叫到高潮| 一本一本综合久久| 亚洲 欧美 日韩 在线 免费| 午夜免费激情av| av福利片在线观看| 久久天躁狠狠躁夜夜2o2o| 手机成人av网站| 欧美zozozo另类| 亚洲欧美日韩东京热| 性色av乱码一区二区三区2| 日韩国内少妇激情av| 校园人妻丝袜中文字幕| 成人高潮视频无遮挡免费网站| 最近中文字幕高清免费大全6| 内射极品少妇av片p| 美女xxoo啪啪120秒动态图| 欧美高清成人免费视频www| 99在线人妻在线中文字幕| 成人综合一区亚洲| АⅤ资源中文在线天堂| 大又大粗又爽又黄少妇毛片口| 成年女人永久免费观看视频| 99精品在免费线老司机午夜| 岛国在线免费视频观看| 欧美一区二区精品小视频在线| 午夜福利成人在线免费观看| 亚洲欧洲日产国产| 99久久九九国产精品国产免费| 国产精品一区www在线观看| 男女下面进入的视频免费午夜| 一级毛片aaaaaa免费看小| а√天堂www在线а√下载| 亚洲电影在线观看av| 久久精品久久久久久久性| 日韩欧美三级三区| 国产成人freesex在线| 男人和女人高潮做爰伦理| 国产一区二区三区av在线 | 亚洲成人中文字幕在线播放| 国产探花在线观看一区二区| 91久久精品国产一区二区成人| 国产爱豆传媒在线观看| 国产精品伦人一区二区| 在线观看一区二区三区| 日本三级黄在线观看| 日韩三级伦理在线观看| 赤兔流量卡办理| 校园人妻丝袜中文字幕| 12—13女人毛片做爰片一| 国产精品久久电影中文字幕| 淫秽高清视频在线观看| 色综合色国产| 最近2019中文字幕mv第一页| 婷婷六月久久综合丁香| 亚洲欧美成人精品一区二区| 禁无遮挡网站| 日韩中字成人| a级一级毛片免费在线观看| 中国美女看黄片| 国产极品天堂在线| 99热这里只有是精品50| 国产高清有码在线观看视频| 日本免费一区二区三区高清不卡| 在线播放国产精品三级| 超碰av人人做人人爽久久| 色哟哟哟哟哟哟| 亚洲人成网站在线播| 国产视频首页在线观看| 免费观看在线日韩| 校园春色视频在线观看| 国产一区二区三区在线臀色熟女| 亚洲电影在线观看av| 老司机福利观看| 亚洲精品国产成人久久av| 全区人妻精品视频| 精品日产1卡2卡| 91久久精品国产一区二区三区| 在线天堂最新版资源| 嘟嘟电影网在线观看| 国产一区二区亚洲精品在线观看| 国产老妇女一区| 男人狂女人下面高潮的视频| 成年女人永久免费观看视频| 99视频精品全部免费 在线| 99热全是精品| 亚洲美女视频黄频| 欧美最黄视频在线播放免费| 免费大片18禁| 在线天堂最新版资源| 青春草国产在线视频 | 国产乱人偷精品视频| 欧美日韩乱码在线| 日韩人妻高清精品专区| 国产亚洲av片在线观看秒播厂 | 男女做爰动态图高潮gif福利片| 国产亚洲精品av在线| 麻豆国产av国片精品| 久久久成人免费电影| 日韩国内少妇激情av| 午夜激情欧美在线| av又黄又爽大尺度在线免费看 | 日韩欧美 国产精品| 天堂影院成人在线观看| 两个人视频免费观看高清| 日产精品乱码卡一卡2卡三|