亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        工藝參數(shù)對雙輥薄帶鑄軋中心層偏析影響

        2018-07-09 08:50:26張德全許志強安仲健杜鳳山魏潔平
        中國有色金屬學(xué)報 2018年6期

        張德全,許志強,安仲健,杜鳳山,魏潔平

        ?

        工藝參數(shù)對雙輥薄帶鑄軋中心層偏析影響

        張德全,許志強,安仲健,杜鳳山,魏潔平

        (燕山大學(xué) 國家冷軋板帶裝備及工藝工程技術(shù)研究中心,秦皇島 066004)

        雙輥薄帶鑄軋過程中工藝參數(shù)對中心層偏析有重大影響,針對雙輥薄帶鑄軋機的特點,利用多相流技術(shù)建立完全耦合溶質(zhì)場、流場、溫度場的數(shù)學(xué)模型,對鑄軋熔池中液相?液固兩相?固相同時存在的復(fù)雜凝固過程進(jìn)行模擬仿真,研究分析鑄軋熔池中Mg、Si元素溶質(zhì)場與流場分布特征,對比研究不同工藝參數(shù)(澆鑄溫度、輥縫寬度)對鑄軋薄帶中心層偏析的影響規(guī)律,并結(jié)合鋁合金6061的鑄軋實驗對結(jié)論進(jìn)行驗證。

        雙輥薄帶鑄軋;中心線偏析;流場;Kiss點高度

        雙輥薄帶鑄軋(Twin-roll strip casting)屬于近終型薄帶材生產(chǎn)工藝之一,其有效地促進(jìn)了現(xiàn)代冶金工業(yè)的轉(zhuǎn)型發(fā)展[1]。偏析是在軋制過程中比較常見的現(xiàn)象,中心層偏析屬于宏觀偏析種類之一,其對鑄軋薄帶的質(zhì)量和力學(xué)性能產(chǎn)生惡劣影響[2],嚴(yán)重時薄帶會出現(xiàn)中心分層現(xiàn)象[3],尤其對鋁帶材料的性能影響較重[4]。溶質(zhì)元素在固液兩相中溶解度不同以及固液兩相中的溶質(zhì)輸運是產(chǎn)生中心層偏析的根本原因,但其如何影響中心層偏析以及其影響效果目前尚不清楚[5],目前薄帶連鑄過程中宏觀偏析的數(shù)值模擬較多,但雙輥薄帶鑄軋過程中的數(shù)值模擬往往是脫離溶質(zhì)場孤立地研究熔池內(nèi)的流場分布。由于澆鑄與凝固同步進(jìn)行,澆鑄過程中的動量、溫度差以及濃度差引起的對流會影響鑄軋帶坯的組織和成分,因此在分析熔池的溶質(zhì)場以及中心偏析時,必須將流場與溫度場、溶質(zhì)場耦合建立適用于雙輥薄帶鑄軋的數(shù)學(xué)模型[6]。

        1 多場耦合數(shù)學(xué)模型的構(gòu)建

        1.1 研究方法及基本假設(shè)

        鑄軋熔池內(nèi)同時存在液相、固相以及固液兩相區(qū),在Kiss點以上主要為金屬液的凝固以及坯殼生長,Kiss點以下金屬液已經(jīng)完全凝固,主要發(fā)生帶坯的固態(tài)軋制。中心層偏析的產(chǎn)生除了與溫度場有關(guān),更與熔池內(nèi)部流場相關(guān)。鑄軋力通過凝固坯殼對流場的影響與熔池內(nèi)金屬液自身流動相比影響較小故而忽略,熔池區(qū)域統(tǒng)一處理為廣義流體,忽略軋制力對熔池內(nèi)部流場的影響,采用該方法對熔池內(nèi)溶質(zhì)場求解更加符合實際[7]。

        考慮到幾何模型的對稱性以及熔池內(nèi)復(fù)雜的物理現(xiàn)象,為了簡化數(shù)值模擬模型鑄軋熔池區(qū)域需作如下假設(shè)。

        1) 由實際三維結(jié)構(gòu)簡化為二維軸對稱模型;

        2) 熔池內(nèi)金屬液視為不可壓縮Newton流體;

        3) 凝固過程遵循局部組分平衡;

        4) 忽略鑄軋輥面以及薄帶熱變形,帶坯與鑄軋輥表面無滑移;

        5) 熔池內(nèi)金屬液的流動視為湍流運動[8];

        6) 熔池中液面認(rèn)為水平固定,忽略液面的波 動[9];

        7) 熔池凝固組織全部視為柱狀晶,凝固坯殼以柱狀晶形式沿與熱傳遞相反á001?方向生長,不考慮間隙生成;同時固態(tài)相變潛熱遠(yuǎn)小于凝固潛熱,故忽略固態(tài)相變潛熱。

        1.2 幾何模型及邊界條件的確定

        熔池計算區(qū)域如圖1所示,邊界條件如下。

        1) 入口(Entrance—)邊界條件

        布流器入口采用Velocity-inlet條件,澆鑄速度in根據(jù)出入口質(zhì)量守恒來計算。由于出坯速度恒定,因此入口速度為根據(jù)出口速度變化的恒定值[10]。

        2) 對稱 (Symmtry—)邊界條件

        由于幾何對稱性其熱邊界條件設(shè)置為無熱量傳遞。

        3) 鑄軋輥(Contact surface—)壁面邊界條件

        根據(jù)鑄軋帶坯與鑄軋輥面無滑移的假設(shè),凝固帶坯沿輥面切向方向運動,紊流動能及耗散率均為0。

        4) 出口(Exit—)邊界條件

        出口處鋁液流動狀態(tài)充分發(fā)展,除壓力外的各個變量沿軋制方向的梯度值為零,且與時間無關(guān),瞬態(tài)計算過程中每一時刻都是充分發(fā)展的。

        圖1 鑄軋熔池計算區(qū)域及邊界條件

        2 數(shù)值模擬與結(jié)果分析

        以雙輥鑄軋設(shè)備熔池區(qū)域為研究對象研究不同工藝參數(shù)(輥縫寬度、熔池溫度)對熔池內(nèi)金屬液流動、溶質(zhì)輸運以及薄帶寬度方向合金元素中心層偏析的影響,計算分析的實驗方案如表1和表2所列。

        表1 不同輥縫寬度下數(shù)值模擬實驗方案

        表2 不同澆鑄溫度下數(shù)值模擬實驗方案

        2.1 鑄軋熔池流場與溶質(zhì)場特征

        由于鋁合金液凝固過程存在選分結(jié)晶現(xiàn)象,加之鋁液流動造成的對流擴散,Mg、Si元素在鑄軋熔池內(nèi)的分布極度不均勻。鑄軋熔池內(nèi)流場、溶質(zhì)場分布特征如下:

        鑄軋熔池區(qū)域流場分布如圖2所示,熔池主要存在2個渦流區(qū)(Vortex zone)[11],Vortex 1靠近自由液面,其主要影響鋁液的滯留時間和液面波動,釋放推送至熔池表面鋁液部分潛熱,另一方面將鋁液推送至鑄軋輥面進(jìn)一步對流換流形成初始坯殼;Vortex 2的形成主要是鋁液凝固造成黏度增大以及鑄軋輥的軋制作用造成的,Vortex 2呈細(xì)長楔形并且尖部深入到Kiss point附近區(qū)域。Vortex 2不僅影響熔池固液兩相區(qū)所占比例,而且能使偏聚在Kiss point附近富溶質(zhì)鋁液返回到熔池進(jìn)行混合再分配。渦流的存在不僅可以促使富含溶質(zhì)的鋁液返回到熔池中進(jìn)行再混合均勻化,另一方面該渦流的存在還能夠以較大的速度沖擊枝晶臂,以及芯部較高的溫度能熔斷枝晶,這兩方面都有利于增加形核率,促使Kiss點區(qū)域溫度場均勻化從而有助于芯部等軸晶的形成[12]。

        圖2 鑄軋熔池液相分?jǐn)?shù)及流線

        2.2 澆鑄溫度對中心層偏析的影響

        合金元素的宏觀偏析比定義為各部位局部濃度與初始濃度之比即/0,且/0>1稱為正偏析,/0<1稱為負(fù)偏析[13]。不同澆鑄溫度下薄帶厚度方向Mg、Si元素宏觀偏析比分布如圖3所示,Mg、Si元素含量由鑄軋薄帶邊部向中心線區(qū)域逐漸增高,且越靠近芯部Mg、Si元素含量增長越快,并在中心線處達(dá)到最大值時形成正偏析尖峰,這主要是由于各部位材料凝固時間和位置不同造成的,越靠近中心線區(qū)域鋁液中合金元素含量也就越高。澆鑄溫度由680 ℃上升至700 ℃時Mg、Si元素中心線位置偏析比尖峰值基本不變,但繼續(xù)上升至720 ℃時中心線處Mg、Si元素的宏觀偏析比值變化較大,分別達(dá)到了3.37和2.34。3種澆鑄條件下Si元素正偏析區(qū)域尺寸均為0.32 mm,占鑄軋帶坯帶坯厚度方向16%,680 ℃和700 ℃澆鑄溫度下Mg元素正偏析尺寸基本相同,但是720 ℃澆鑄溫度條件下Mg元素正偏析尺寸由0.84 mm減小至 0.72 mm。

        由圖4和5可以看出,隨著澆鑄溫度的升高,Vortex 1逐漸向熔池液面移動并且渦流區(qū)域范圍減小,而入口鋁液沖擊深度保持17.5 mm左右不變,Vortex 1的退化導(dǎo)致鋁液進(jìn)入熔池后大部分直接沖向輥面,不利于鋁液的散熱和凝固;Vortex 2回流區(qū)速度拐點2下降了5.3 mm,峰值速度也下降了25.3 mm/s,Vortex 2的變化對Kiss點區(qū)域的攪拌強度降低,高濃度鋁液得不到新鮮鋁液的有效稀釋從而加重中心層偏析。720 ℃澆鑄條件下Vortex 3下移并最終與Vortex 2相連導(dǎo)致了流動通道消失(Disappearance of flow passage),從而無低濃度金屬液在Vortex 2與Vortex 3之間沖擊高濃度鋁液,Vortex 2 effect消失后偏析曲線由邊部呈現(xiàn)持續(xù)上升趨勢。

        圖3 不同澆鑄溫度下薄帶厚度方向Mg、Si元素分布

        圖4 不同澆鑄溫度下熔池液相分?jǐn)?shù)及流線

        圖5 不同澆鑄溫度下熔池中心線速度分布

        2.3 輥縫寬度對中心層偏析的影響

        不同輥縫寬度下薄帶厚度方向Mg、Si元素含量分布如圖6所示,隨著輥縫由2 mm增大至4 mm,Mg、Si元素在中心線處宏觀偏析比分別由3.14、2.02上升至5.08、3.05,中心層偏析加重。

        一方面根據(jù)凝固凝殼平面長大平方根定律可知凝固時間與坯殼厚度呈平方關(guān)系,因此隨著輥縫寬度增大,凝固時間激增從而凝固速率降低,在較低的凝固速率條件下固液界面析出的溶質(zhì)元素在鋁液對流作用下被沖刷帶入熔池芯部與新鮮鋁液再混合從而使邊部溶質(zhì)元素含量維持在較低值。另一方面溶質(zhì)元素在輥縫出口處質(zhì)量守恒必須滿足式(1),在輥縫寬度增大條件下進(jìn)入熔池溶質(zhì)含量上升,而鑄軋薄帶邊部溶質(zhì)元素含量基本不變且保持較低值,因此熔池輥縫越寬越會有較多的溶質(zhì)聚集在芯部,中心層偏析程度越嚴(yán)重。

        式中:為熔池輥縫寬度,mm;in為入口澆鑄速度,mm/s;out為出口速度,即軋制速度,mm/s;為澆鑄入口尺寸,mm。

        在相同布流器尺寸和軋制速度條件下,隨著輥縫由2 mm擴大至4 mm,根據(jù)出入口質(zhì)量守恒入口速度由10 mm/s擴大至20 mm/s。由圖7和圖8可知,隨著入口速度由10 mm/s上升至15 mm/s時入口沖擊深度(Impact depth)增加了10.2 mm,而輥縫寬度由3 mm升至4 mm時沖擊深度(Impact depth)基本改變,沖擊深度增加能夠?qū)⑿迈r鋁液盡可能帶入熔池底端沖擊高濃度回流鋁液起到稀釋作用,同時帶動較多的鋁液隨之流動,降低了卷渣的可能性,便于夾雜物和高濃度鋁液上浮。凝固初期靠近輥面的自然對流引起較強的渦流Vortex 1,其能夠?qū)⒛糖把氐娜苜|(zhì)帶入熔池內(nèi)部起到稀釋作用。隨著入口速度升高回流Vortex 1水平方向逐漸向熔池中心線方向靠近,豎直高度逐漸降低向遠(yuǎn)離液面位置移動,同時渦流區(qū)域范圍逐漸增大,初始坯殼附近開始形成Vortex 3,并且隨著輥縫增大Vortex 3回流區(qū)域范圍也逐漸增大。渦流Vortex 1位置的改變不利于初始坯殼固液界面處析出溶質(zhì)元素通過對流沖刷進(jìn)入熔池,邊部溶質(zhì)元素升高速度加快,而Vortex 3的出現(xiàn)能夠形成Vortex 3抑制合金元素的持續(xù)上升。

        圖6 不同輥縫下薄帶厚度方向Mg、Si元素分布

        圖7 輥縫不同時熔池液相分?jǐn)?shù)以及流線

        伴隨著Vortex 1位置的改變,Vortex 2的位置基本沒有太大變化,但是Vortex 2區(qū)域范圍急劇減小,由圖8可以看出,隨著輥縫由2 mm升至3 mm渦流2在中心線上的速度銳減,攪動強度降低。輥縫為4 mm時鋁液回流區(qū)Vortex 2同上方渦流之間的聯(lián)系被截斷(Disappearance of back flow passage),部分注入熔池的鋁液未經(jīng)過出口處回流區(qū)冷卻也未被推送至初始接觸區(qū)進(jìn)行冷卻而直接被卷向熔池出口。

        當(dāng)輥縫較大時,較多的鋁液未經(jīng)過Vortex 1對流換熱冷卻,也未被推送至液面受到初始接觸條件冷卻因而Kiss點位置大幅度降低(見圖9),輥縫由2 mm升至3mm時Kiss點位置下降了7.02 mm,而輥縫繼續(xù)上升至4 mm時Kiss位置只下降了3.40 mm,但是輥縫3 mm和4mm時Mg、Si元素中心偏析差異卻較大,因此輥縫較大時Kiss點位置對中心偏析的影響減小。

        圖8 輥縫不同時熔池中心線速

        圖9 輥縫不同時熔池中心線液相分?jǐn)?shù)

        由圖8可以看出中心線上兩個速度拐點之間的速度峰值現(xiàn)象逐漸消失,Vortex 2造成的速度拐點之間尺寸也由28.2 mm銳減至8.6 mm,Kiss點附近高濃度鋁液的攪拌效果減弱,無法將偏聚的中心線處的高濃度鋁液帶回熔池進(jìn)行再分配,最終隨兩側(cè)坯殼的焊合而封閉在中心線處加重中心偏析。

        3 實驗驗證

        本實驗中通過雙輥鑄軋設(shè)備來驗證鑄軋6061鋁合金薄帶的可行性(見圖10、11),參照6061鋁合金的鑄造工藝參數(shù)特性,經(jīng)過多次反復(fù)的實驗并在實驗過程中對實驗參數(shù)進(jìn)行分析,初步確定了在現(xiàn)有實驗室雙輥鑄軋設(shè)備軋制能力下鑄軋6061鋁合金薄帶的合理工藝參數(shù)[14](見表3),成功軋制的6061鋁合金薄帶。

        圖10 雙輥薄帶鑄軋機照片[15]

        圖11 雙輥薄帶鑄軋工藝流程示意圖[16]

        表3 6061鋁合金薄帶鑄軋實驗工藝參數(shù)

        為驗證所構(gòu)建數(shù)學(xué)模型的正確性,需要將帶寬方向溶質(zhì)元素的分布檢測值與數(shù)值仿真模擬值進(jìn)行對比,本實驗采用斷面取點掃描的方式檢測合金元素含量。圖12、13所示為實驗結(jié)果的金相組織圖和元素掃描圖。

        圖14所示為相同工藝參數(shù)下鑄軋薄帶厚度方向溶質(zhì)元素宏觀偏析比的模擬結(jié)果和實際檢測結(jié)果,模擬結(jié)果與實際檢測得到的宏觀偏析比分布極為相似,在薄帶中心部位出現(xiàn)較為明顯的正偏析尖峰,中心線兩側(cè)也存在由于回流漩渦造成的成分偏低區(qū)域。

        圖12 鑄軋薄帶的顯微組織

        圖13 檢測點分布以及EDS掃描結(jié)果

        同時圖14也表明數(shù)值模擬結(jié)果與實測結(jié)果存在一定的差異,一方面中心位置合金元素含量實測值比模擬值較小,這主要是合金元素在熔煉、鑄軋過程中部分合金元素氧化成為熔渣減少了合金元素的實際含量;另一方面模擬結(jié)果顯示中心正偏析尺寸比實際檢測尺寸要大,Mg、Si元素上升速率相對較慢,這表明實際上盡管凝固末期熔池區(qū)域在Kiss點附近Vortex 2的對流對用已經(jīng)很微弱,但是在鑄軋輥的軋制作用下Kiss點附近糊狀區(qū)高濃度鋁液仍會有一定的流動,從而迫使Kiss點附近部分高濃度鋁液回流返回熔池,降低了合金元素含量上升速率。

        圖14 薄帶厚度方向Mg、Si元素偏析比模擬結(jié)果與實測值

        4 結(jié)論

        1) 不同澆鑄溫度下,Mg、Si正偏析尖峰均出現(xiàn)在薄帶中心線處,但澆鑄溫度由680 ℃上升至700 ℃時Mg、Si元素中心線位置偏析比尖峰值基本不變,但繼續(xù)上升至720 ℃時中心線處Mg、Si元素的宏觀偏析比值變化較大。

        2) 隨著輥縫由2 mm增大至4 mm,Mg、Si元素在中心線處宏觀偏析比分別由3.14、2.02上升至5.08、3.05,中心層偏析加重。

        3) 通過將相同工藝參數(shù)下鑄軋薄帶厚度方向中心層偏析模擬值與實測值對比發(fā)現(xiàn),兩者之間的趨勢基本一致,但由于數(shù)值模擬過程中忽略了鑄軋輥的軋制力作用,因而中心線位置處模擬值比實測值偏高。

        [1] ZAPISKALOV N. Comparison of continuous strip casting with conventional technology[J]. ISIJ International, 2003, 43(8): 1115?1127.

        [2] BIROL Y, GEBZE K. Centerline segregation in a twin-roll cast AA8011 alloy[J]. Aluminium, 1998, 77(5): 318?321.

        [3] HADAZADEH A, WELL M A. Inverse and centerline segregation formation in twin roll cast AZ31 magnesium alloy[J]. Materials Science & Technology, 2015: 1715?1726.

        [4] 徐國棟, 小森重喜, 石洪志. 寶鋼連鑄板坯中心偏析改進(jìn)[C]// 中國金屬學(xué)會. 2001中國鋼鐵年會論文集. 下卷. 北京: 冶金工業(yè)出版社, 2001: 816?817. XU Guo-dong, XIAO Sen-chong-xi, SHI Hong-zhi. Continuous casting slab center segregation improvement of Baosteel[C]// The Chinese Society for Metals. 2001 China’s Steel Conference Proceedings. Book Two. Beijing: Metallurgical Industry Press, 2001: 816?817.

        [5] 張興中, 廖 鵬, 王明林. 雙輥薄帶連鑄技術(shù)狀況的調(diào)查分析[J]. 鋼鐵, 2010, 45(3): 13?17. ZHANG Xing-zhong, LIAO Peng, WANG Ming-lin. Investigation on development of twin roll cast strip[J]. Iron and Steel, 2010, 45(3): 13?17.

        [6] 杜鳳山, 呂 征, 黃華貴, 許志強, 安仲健. 雙輥薄帶鑄軋中心線偏析機理與實驗研究[J]. 中國有色金屬學(xué)報, 2015, 25(10): 2739?2744. DU Feng-shan, Lü Zheng, HUANG Hua-gui, XU Zhi-qiang, AN Zhong-jian. Mechanism and experimental research on centerline segregation of twin-roll strip casting[J]. The Chinese Journal of Nonferrous Metals, 2015, 25(10): 2739?2744.

        [7] 孟哲儒. 振動對雙輥薄帶鑄軋流場溫度場影響的研究[D]. 秦皇島: 燕山大學(xué), 2014. MENG Zhe-ru. Effect of vibration on twin-roll strip roll casting's flow field and temperature field[D]. Qinhuangdao: Yanshan University, 2014.

        [8] GUTHRIE R I L, TAVARES R P. Mathematical and physical modeling of steel flow and solidification in twin-roll/horizontal belt thin-strip casting machines[J]. Applied Mathematical –Modeling, 1998, 22(11): 851?872.

        [9] ABOUTALEBI M R, HASAN M, GUTHRIE R I L. Coupled turbulent flow, heat, and solute transport in continuous casting processes[J]. Metallurgical and Materials Transactions B, 1995, 26(4): 731?744.

        [10] LI Q, ZHANG Y K, LIU L G, ZHANG P, ZHANG Y, FANG Y. Effect of casting parameters on the freezing point position of the 304 stainless steel during twin-roll strip casting process by numerical simulation[J]. Journal of Materials Science, 2012, 47(9): 3953?3960.

        [11] 司開田. 1235 合金鑄軋帶材中心層偏析機制的研究及其預(yù)防措施[J]. 輕合金加工技術(shù), 2007, 35(10): 5?6. SI Kai-tian. Research on segregation and its mechanism of 1235alloy cast-rolling stock and its preventive measures[J]. Light Alloy Fabrication Technology, 2007, 35(10): 5?6.

        [12] ZHAO H, LI P, HE L. Coupled analysis of temperature and flow during twin-roll casting of magnesium alloy strip[J]. Journal of Materials Processing Technology, 2011, 211(6): 1197?1202.

        [13] JING Cai-liang, TIAN Zhi-hong, WANG Xin-hua. A simple micro-segregation model for steel billet[J]. Journal of Iron and Steel Research, 2013, 25(6): 18?23.

        [14] 王 娜, 周志敏, 路貴民, 趙大志. 工藝條件對6061鋁合金近液相線鑄造微觀組織的影響[J]. 特種鑄造及有色合金, 2009, 29(6): 508?511.WANG Na, ZHOU Zhi-min, LU Gui-min, ZHAO Da-zhi. Effect of technological conditions on 6061 aluminum alloy near liquidus casting microstructure[J]. Special Casting & Nonferrous Alloys, 2009, 29(6): 508?511.

        [15] 杜鳳山, 許志強, 路鵬程, 張 沛. 振動式雙輥薄帶鑄軋機: 中國, CN101181719A[P]. 2008?05?21. DU Feng-shan, XU Zhi-qiang, LU Peng-cheng, ZHANG Pei. Vibrating twin-rolls strip rolling machine: China, CN101181719A[P]. 2008?05?21.

        [16] LEE D S, LEE J S, KANG T. Robust molten steel level control in a strip-casting process[J]. ISIJ International, 2005, 45(8): 1165?1172.

        Effect of process parameters on centerline segregation of twin-roll strip casting

        ZHANG De-quan, XU Zhi-qiang, AN Zhong-jian, DU Feng-shan, WEI Jie-ping

        (National Engineering Research Center for Equipment and Technology of Cold Strip Rolling, Yanshan University, Qinhuangdao 066004, China)

        The technical parameters of twin-roll strip casting have significant impact on the centerline segregation. The mathematical model in completely coupled temperature field, flow field and solute field was established by using multiphase flow technology for the characteristics of twin-roll strip casting mill. The solidification process of simultaneous existence of liquid phase, liquid-solid-phase and solid phase in molten pool was simulated based on the coupled mathematical model. The distribution characteristics of Mg and Si solute field and flow field in the molten pool were studied by the simulation results, and the influence of different technical parameters (the width of the roll seam and casting temperature) on the centerline segregation of the thin strip was studied. The research conclusion was verified by the casting and rolling experiments of the aluminum alloy 6061.

        twin-roll strip casting; centerline segregation; flow field; Kiss point height

        Project(51374184) supported by the National Natural Science Foundation of China; Project (E2013203342) supported by the National Natural Science Foundation of Hebei Province—Joint Steel Research Fund Project, China; Project(E2017203043) supported by the Key Project of National Natural Science Foundation of Hebei Province, China

        2017-03-17;

        2017-06-07

        XU Zhi-qiang; Tel: +86-13930357680; E-mail: xzq@ysu.edu.cn

        國家自然科學(xué)基金資助項目(51374184);河北省自然科學(xué)基金—鋼鐵聯(lián)合研究基金資助項目(E2013203342);河北省自然科學(xué)基金重點資助項目(E2017203043)

        2017-03-17;

        2017-06-07

        許志強,教授,博士;電話:13930357680;E-mail: xzq@ysu.edu.cn

        10.19476/j.ysxb.1004.0609.2018.06.06

        1004-0609(2018)-06-1128-08

        TG233.6

        A

        (編輯 王 超)

        精品日韩欧美一区二区三区在线播放| 毛片a级毛片免费观看| 天天看片视频免费观看| 久久精品免费无码区| 在线精品亚洲一区二区三区 | 无套内射无矿码免费看黄| 国产99在线视频| 午夜日本理论片最新片| 开心五月婷婷激情综合网| 亚洲成av人片在线观看www| 国产普通话对白视频二区| 亚洲福利第一页在线观看| 亚洲中文字幕九色日本| 亚洲a∨无码男人的天堂| 国产日韩欧美亚洲精品中字| 在线视频一区二区亚洲| 开心五月骚婷婷综合网| 午夜精品久久久久久久99热| 亚洲深深色噜噜狠狠爱网站| 亚洲人成无码网站十八禁| 一区二区三区免费观看日本 | 国精无码欧精品亚洲一区| 国产精品久久婷婷婷婷| 亚洲精品天堂日本亚洲精品| 真实的国产乱xxxx在线| 亚洲精品无码人妻无码| 国产精品亚洲综合色区丝瓜| 亚洲免费一区二区三区四区| 少妇被粗大的猛烈进出免费视频| 色综合自拍| 亚洲色图视频在线观看,| 色哟哟亚洲色精一区二区| 男女性高爱潮免费网站| 国产精品白浆视频免费观看| 亚洲一本二区偷拍精品| 亚洲人成网网址在线看| 无码视频一区二区三区在线观看| 国产91精品丝袜美腿在线| 日韩人妻精品中文字幕专区| 亚洲av无码av制服另类专区| 国产欧美日韩不卡一区二区三区 |