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

    基于DEM法氫氣直接還原豎爐內(nèi)物料運(yùn)動(dòng)行為研究

    2018-06-22 07:53:34白明華葛俊禮
    重型機(jī)械 2018年3期
    關(guān)鍵詞:豎爐爐料徑向

    劉 棟, 龍 鵠, 白明華, 葛俊禮, 徐 寬

    (1. 寶武集團(tuán)廣東韶關(guān)鋼鐵有限公司,廣東 韶關(guān) 512123; 2. 燕山大學(xué) 國(guó)家冷軋板帶裝備及工藝工程技術(shù)研究中心,河北 秦皇島 066004; 3. 河北港口集團(tuán)港口機(jī)械有限公司,河北 秦皇島 066004)

    0 前言

    由化石燃料燃燒產(chǎn)生的溫室氣體導(dǎo)致的氣候變暖現(xiàn)象已成為全球共同面對(duì)的問(wèn)題,鋼鐵行業(yè)是工業(yè)當(dāng)中的溫室氣體排放大戶。據(jù)統(tǒng)計(jì)[1,2]鋼鐵廠的噸鋼CO2排放量約為1 832 kg,其中高爐煉鐵占全流程的76.1%。同時(shí)高爐由于大量使用焦炭以及燒結(jié)礦,是SO2、NOx、粉塵的主要來(lái)源[1,2]。以碳還原鐵礦石為高爐技術(shù)核心的傳統(tǒng)長(zhǎng)流程鋼鐵工業(yè)亟待轉(zhuǎn)型發(fā)展[3]。

    氫氣(H2)作為還原劑在豎爐內(nèi)生產(chǎn)直接還原鐵,可降低煉鐵過(guò)程中的碳耗并減少CO2排放。該還原過(guò)程在低于鐵礦石軟化溫度下進(jìn)行,因此氫氣還原豎爐是一種典型的氣固逆流式移動(dòng)床[4-6]。物料在豎爐內(nèi)下降過(guò)程中的運(yùn)動(dòng)行為直接影響爐內(nèi)流場(chǎng)分布,與爐況的穩(wěn)定性和生產(chǎn)率息息相關(guān)。當(dāng)固態(tài)物料呈平推流的狀態(tài)均勻向下移動(dòng)時(shí),可使得物料在爐內(nèi)停留時(shí)間一致,即得到較為均勻的還原,有利于操作穩(wěn)定性。否則,當(dāng)運(yùn)動(dòng)速率沿徑向分布不均,會(huì)導(dǎo)致部分物料未得到有效還原,從而降低生產(chǎn)效率[7,8]。前期通過(guò)實(shí)驗(yàn)研究已知?dú)錃膺€原速率較快[9,10],從反應(yīng)速率角度考慮,其還原段高度與傳統(tǒng)豎爐相比可有效降低,但從爐內(nèi)物料運(yùn)動(dòng)均勻性考慮,高度的降低可能會(huì)帶來(lái)新的問(wèn)題。

    本文通過(guò)離散單元法(DEM)建立豎爐內(nèi)爐料顆粒運(yùn)動(dòng)模型,針對(duì)不同高度條件下?tīng)t內(nèi)物料運(yùn)動(dòng)情況進(jìn)行分析,對(duì)比氫氣豎爐與傳統(tǒng)Midrex豎爐內(nèi)物料沿徑向速率分布的差別,為氫氣直接還原豎爐爐型及設(shè)備設(shè)計(jì)優(yōu)化提供理論指導(dǎo)。

    1 豎爐及排料機(jī)參數(shù)

    傳統(tǒng)Midrex及HYL直接還原豎爐主體部分通常由還原段、過(guò)渡段、冷卻段三段構(gòu)成,爐料自豎爐頂部的料斗、爐料分配器進(jìn)入還原段,在還原段內(nèi)下降,與逆流的還原氣之間發(fā)生反應(yīng),完成從氧化鐵礦石到海綿鐵之間的轉(zhuǎn)變。冷卻段下方由擺式或輪式卸料機(jī)負(fù)責(zé)卸料。Corex熔融煉鐵法是通過(guò)豎爐底部的螺旋排料機(jī)將還原得到的熱態(tài)海綿鐵排出豎爐,并由螺旋機(jī)下方的下降管進(jìn)入下部的熔融氣化爐。

    圖1 直接還原鐵熱送至電爐示意圖

    傳統(tǒng)豎爐冷卻段較高,為了減少能耗,并有效降低豎爐整體高度,本文研究采用熱送工藝,即通過(guò)熱螺旋輸送機(jī)進(jìn)行豎爐內(nèi)直接還原鐵排料,然后熱送至電爐煉鋼廠,流程如圖1所示,并通過(guò)調(diào)整螺旋機(jī)葉片直徑、軸徑,齒數(shù)等方法調(diào)節(jié)爐內(nèi)物料運(yùn)動(dòng)情況。

    1.1 爐主要參數(shù)

    Siderca廠[11]的海綿鐵產(chǎn)量為100 t·h-1,年產(chǎn)量約為80萬(wàn)噸,豎爐還原段直徑為4.88 m,高度為10 m,高徑比約為2,有效容積和還原段利用系數(shù)按式(1)、(2)計(jì)算得到,分別為186.9 m3和11.9 t·m-3·d-1。當(dāng)產(chǎn)能一定時(shí),還原段利用系數(shù)與有效容積呈反比關(guān)系。氫氣還原效率高,還原段利用效率相應(yīng)增大,在豎爐直徑保持不變的條件下,高度可減小。本文研究的氫氣豎爐直徑參照Siderca數(shù)據(jù),取5 m,高度有所減小,分別為8 m,6 m,與Midrex豎爐直徑為5 m,高度為10 m情況進(jìn)行對(duì)比。為簡(jiǎn)化模擬未考慮爐壁傾角。

    Veff=0.785D2·H

    (1)

    ηeff=24QS/Veff

    (2)

    式中,Veff是還原段有效容積,m3;D是還原段直徑,m;H是還原段高度,m;ηeff是還原段利用系數(shù),t·m-3·d-1;QS是每小時(shí)海綿鐵產(chǎn)量,t·h-1。

    1.2 螺旋輸送機(jī)參數(shù)計(jì)算

    直接還原豎爐中的螺旋輸送機(jī)用于熱海綿鐵的卸料,該設(shè)備對(duì)調(diào)節(jié)爐內(nèi)物料運(yùn)動(dòng)行為,控制均勻的下降流流型具有重要的作用,并進(jìn)一步影響到爐內(nèi)氣流分布以及操作的穩(wěn)定性。

    螺旋輸送機(jī)的關(guān)鍵參數(shù)主要包括葉片直徑、螺距、以及螺旋葉片升角,通過(guò)公式[12,13]對(duì)螺旋機(jī)關(guān)鍵參數(shù)進(jìn)行計(jì)算。

    螺旋直徑D為

    (3)

    式中,D為螺旋直徑,m;K為物料綜合系數(shù);Q為輸送量,t·h-1;φ為充滿系數(shù);γ為物料的堆比重,t·m-3;c為傾角系數(shù)。

    螺旋輸送機(jī)螺距、螺旋軸直徑以及螺旋葉片的升角分別由式(4)~(6)計(jì)算得到。

    S=bD

    (4)

    d=(0.2~0.35)D

    (5)

    α=arctan(S/2πr)

    (6)

    式中,S為螺距,m;b為螺距/螺旋直徑系數(shù);R為螺旋軸的半徑,m。

    本次研究一共設(shè)8臺(tái)螺旋輸送機(jī),沿豎爐底部徑向均勻分布。計(jì)算得到基本模型條件下螺旋排料裝置葉片尺寸參數(shù)見(jiàn)表1。

    表1 基本模型螺旋排料裝置葉片尺寸參數(shù)

    2 DEM模型的建立

    離散元法( Discrete Element Method)是將整個(gè)研究區(qū)域離散為多個(gè)獨(dú)立的運(yùn)動(dòng)單元,且每個(gè)單元具有一定幾何特征(包括形狀、大小、排列等)和物理特性以及化學(xué)特性[14, 15]。各單元運(yùn)動(dòng)過(guò)程受經(jīng)典運(yùn)動(dòng)方程控制。本次研究利用EDEM軟件,采用離散單元法建立數(shù)值模型,對(duì)氫氣還原豎爐內(nèi)爐料的運(yùn)動(dòng)行為進(jìn)行仿真模擬,分析料型分布,研究固態(tài)物料運(yùn)動(dòng)特性,并與傳統(tǒng)Midrex豎爐進(jìn)行對(duì)比。

    2.1 模擬方法

    為了能夠描述顆粒間的接觸,DEM運(yùn)用振動(dòng)運(yùn)動(dòng)方程來(lái)模擬。顆粒之間沿法線、切線方向的接觸力采用“彈簧-阻尼器-摩擦板”的Voigt彈簧阻尼模型來(lái)計(jì)算,如圖2所示。

    圖2 振動(dòng)接觸模型

    假設(shè)兩個(gè)球形顆粒之間發(fā)生線彈性接觸,二者有重疊發(fā)生,此時(shí)單元體之間由于接觸變形所產(chǎn)生的法向力Fn以及切向力Ft的計(jì)算如式(7)和(8)所示[14]。

    (7)

    Ftij=min{μ|Fn,ij|tij,[Kt[Δut,ij+Δφij]+

    (8)

    式中,K是顆粒的彈性系數(shù);η是顆粒的阻尼系數(shù);μ是顆粒的滑動(dòng)摩擦系數(shù);Δu是顆粒運(yùn)動(dòng)過(guò)程中的線位移;Δφ是顆粒運(yùn)動(dòng)過(guò)程中的角位移;nij是顆粒i到顆粒j的法向矢量;tij是顆粒i到顆粒j的切向矢量。

    進(jìn)行數(shù)值求解時(shí)采用中心差分法,更新速度用兩次迭代時(shí)間步長(zhǎng)的中間點(diǎn)來(lái)表達(dá)并計(jì)算位移關(guān)系,從而得到進(jìn)一步的作用力。進(jìn)行時(shí)間t內(nèi)的k次循環(huán)計(jì)算求解,即可計(jì)算得在時(shí)間t內(nèi)任意時(shí)間點(diǎn)的任意顆粒的運(yùn)動(dòng)。

    2.2 顆粒模型的建立

    離散元法模擬計(jì)算時(shí)顆粒數(shù)目受到工作站能力限制,而實(shí)際生產(chǎn)過(guò)程中顆粒數(shù)為千萬(wàn)級(jí),為此在建立顆粒模型過(guò)程中采用顆粒粒徑放大的處理方法,減少模擬時(shí)的顆粒數(shù)量,降低計(jì)算量,并對(duì)有關(guān)問(wèn)題及規(guī)律做定性分析,這種放大粒徑的方法的適用性及準(zhǔn)確度在多篇文獻(xiàn)中得到了驗(yàn)證[14-16]。

    顆粒參數(shù)如表2所示,顆粒密度取球團(tuán)密度,顆粒粒徑放大了10倍。

    表2 顆粒參數(shù)

    與爐料顆粒接觸的爐壁參數(shù)如表3所示。

    表3 豎爐爐壁材料參數(shù)

    螺旋排料裝置葉片材料參數(shù)如表4所示。

    表4 螺旋排料裝置葉片材料參數(shù)

    2.3 物理模型的建立及仿真

    氫氣豎爐還原段物料運(yùn)動(dòng)離散元仿真過(guò)程步驟如下:首先利用solidwrks軟件建立豎爐及螺旋輸送機(jī)三維模型;其次,將幾何模型導(dǎo)入EDEM軟件,并設(shè)定所有模擬所需計(jì)算參數(shù);然后開(kāi)始填充物料,當(dāng)豎爐內(nèi)物料填充率達(dá)到90%時(shí),使螺旋輸送機(jī)旋轉(zhuǎn),并將豎爐內(nèi)物料排出;當(dāng)物料不再排出時(shí)停止模擬計(jì)算,收集計(jì)算結(jié)果數(shù)據(jù) 并加以分析。

    使用Solidworks三維建模軟件建立豎爐爐體的離散單元法仿真模型,包括豎爐本體和八臺(tái)螺旋排料裝置,如圖3所示。豎爐內(nèi)徑為5 m,還原段高度為10 m,爐壁未設(shè)傾斜角。

    圖3 豎爐離散元模擬計(jì)算物理模型

    3 豎爐內(nèi)物料運(yùn)動(dòng)模擬結(jié)果

    3.1 基礎(chǔ)模型模擬結(jié)果

    圖4所示為基礎(chǔ)模型在豎爐不同高度下,水平橫截面物料運(yùn)動(dòng)速度分布云圖,深色表示較快的速率,負(fù)值表示與Z軸方向相反。觀察可知高度為2 m處,中心速率較快,由中心向四周速率沿徑向逐漸減小,中心處物料下降速率最快。物料下降速率沿圓周方向均勻分布,因此縱截面不同高度下的物料下降速率分布可衡量整個(gè)豎爐物料運(yùn)動(dòng)分布情況。

    圖4 豎爐高度2 m處橫截面物料運(yùn)動(dòng)速度分布

    圖5所示為基礎(chǔ)模型縱截面物料運(yùn)動(dòng)速度分布云圖,料層總高度為10 m,觀察可知在豎爐高度1/4處,及2.5 m處上方,物料運(yùn)動(dòng)速率分布較為均勻,而在2.5 m處下方,出現(xiàn)了下降速率分布不均的現(xiàn)象,中心處速率最大,中心與邊緣相比,不受到來(lái)自爐壁的摩擦力,且在下方中心處葉片直接將物料旋轉(zhuǎn)輸出,因此物料下降速率快。越靠近豎爐上方,或靠近爐壁端,物料下降速率越小。

    圖5 豎爐縱截面物料運(yùn)動(dòng)速度分布

    為了量化分析速率變化情況,分別取豎爐高度1 m、2 m、3 m處,厚度為0.5 m的物料沿徑向速率分布,如圖6~8所示。觀察可知,不同高度下的物料下降速率均沿中心向爐壁方向逐漸減小,其中高度為1 m處速率差別最大,為0.34 m·s-1,高度3 m處速率差值較小,為0.1 m·s-1。說(shuō)明爐料位置越低,沿徑向速率分布越不均勻。豎爐靠近中心處(r/R=0.1),爐料不同高度下的速率分別為0.38 m·s-1(H=1 m),0.32 m·s-1(H=2 m),0.28 m·s-1(H=3 m),高度為1 m處的物料下降運(yùn)動(dòng)速率與3 m處的相比增大了35.7%;而靠近爐壁處 (r/R=0.7),爐料不同高度下的速率分別為0.14 m·s-1(H=1 m),0.18 m·s-1(H=2 m),0.22 m·s-1(H=3 m),高度為1 m處的速率與3 m處的相比減小了36.3%。

    通過(guò)對(duì)基礎(chǔ)模型物料速率分布分析可知,當(dāng)模擬豎爐高度為10 m,直徑為5 m(傳統(tǒng)Midrex豎爐還原段)時(shí),2.5 m以上物料速率分布沿徑向較為均勻,2.5 m以下速率分布不均,中心速率較高,靠近爐壁處速率較低。

    通過(guò)之前分析,氫氣還原速率快,較傳統(tǒng)Midrex豎爐,氫氣豎爐的還原段可進(jìn)一步降低,因此本節(jié)研究了豎爐高度的減小對(duì)爐內(nèi)物料運(yùn)動(dòng)的影響,分別對(duì)直徑為5 m,爐內(nèi)物料總高度分別為10 m、8 m、6 m條件下的爐料運(yùn)動(dòng)行為進(jìn)行仿真。

    圖6 豎爐縱截面物料運(yùn)動(dòng)速率分布

    3.2 豎爐高度的影響

    圖7所示為不同料層總高條件下縱截面物料運(yùn)動(dòng)速度分布云圖,觀察可知當(dāng)料層總高度從10 m降低到8 m時(shí),物料運(yùn)動(dòng)速率分布與基本模型基本相近,上層較為均勻,下層不均勻度略有增高,中心深色區(qū)域高度略有增加,從2.5 m增大到2.8 m,占總高度比例從25%增大到35%。當(dāng)總高度進(jìn)一步降低到6 m時(shí),中心深色區(qū)域顯著增大,高度約為3.8 m,占據(jù)總料層高度的63.3%。

    圖7 縱截面物料運(yùn)動(dòng)速度分布

    為了量化分析高度為6 m的還原豎爐物料運(yùn)動(dòng)速率變化情況,與基本模型一樣,分別取豎爐高度1 m、2 m、3 m處,厚度為0.5 m的物料沿徑向速率分布,如圖8所示。高度為1 m處速率差別最大,為0.38 m·s-1,高度3 m處速率差值仍較高,為0.33 m·s-1。豎爐靠近中心處(r/R=0.1),爐料不同高度下的速率分別為0.38 m·s-1(H=1 m),0.37 m·s-1(H=2 m),0.36 m·s-1(H=3 m),速率均較高,且變化不大;靠近爐壁處 (r/R=0.7),爐料不同高度下的速率分別為0.05 m·s-1(H=1 m),0.08 m·s-1(H=2 m),0.10 m·s-1(H=3 m),總體靠近爐壁處速率均較小。豎爐3.8 m以下物料運(yùn)動(dòng)速率沿徑向變化均較大。

    圖8 豎爐縱截面物料運(yùn)動(dòng)速率分布(H=6 m)

    3.3 影響效果分析

    為了進(jìn)一步量化豎爐高度對(duì)物料運(yùn)動(dòng)行為的影響,對(duì)不同高度下的模型分別取豎爐高度1 m、2 m、3 m處物料速率沿徑向分布進(jìn)行對(duì)比。

    圖9為高度1 m處不同模型的對(duì)比,三種模型物料運(yùn)動(dòng)速率沿徑向顯著減小,其中10m和8m高豎爐變化趨勢(shì)相近,最大速率差分別為0.34 m·s-1、0.35 m·s-1,兩條曲線接近重合。6 m高豎爐速率差別最大,為0.38 m·s-1。

    圖9 不同豎爐高度縱截面物料運(yùn)動(dòng)速率分布(H=1 m)

    圖10為高度2 m處不同模型的對(duì)比,與1 m處相比,10 m和8 m高豎爐物料運(yùn)動(dòng)速率沿徑向降低幅度顯著減小,其中最大速率差分別為為0.2 m·s-1,0.22 m·s-1,分別降低了41.2%、37.1%。6 m高豎爐速率差仍然較大,為0.36 m·s-1,與1 m處相比僅降低5.3%。

    圖11為高度3 m處不同模型的對(duì)比,與1 m處相比,10 m和8 m高豎爐物料運(yùn)動(dòng)速率沿徑向降低幅度進(jìn)一步減小,其中最大速率差分別為為0.1 m·s-1,0.07 m·s-1,分別降低了72.2%、80%。6 m高豎爐速率差仍然較大,為0.33 m·s-1,與1 m處相比僅降低13.2%。

    圖10 不同豎爐高度縱截面物料運(yùn)動(dòng)分布(H=2 m)

    圖11 不同豎爐高度條件下縱截面物料運(yùn)動(dòng)速率分布(H=3 m)

    通過(guò)上述定性定量分析可知,豎爐高度從10 m降低到8 m時(shí),對(duì)物料運(yùn)動(dòng)變化影響較小,而當(dāng)進(jìn)一步降低到6 m時(shí),影響顯著。豎爐高度為6 m,直徑為5 m,即高徑比為1.2時(shí),物料在高度3.8 m以下運(yùn)動(dòng)速率沿徑向分布極不均勻。為了降低這種不均勻性,獲得更理想的流型分布,可通過(guò)改變螺旋輸送機(jī)葉片尺寸的方法,研究其對(duì)爐料運(yùn)動(dòng)行為的影響。

    4 結(jié)論

    通過(guò)離散單元法建立了豎爐內(nèi)爐料顆粒運(yùn)動(dòng)仿真模型,研究了不同高度及螺旋輸送機(jī)參數(shù)對(duì)爐內(nèi)物料運(yùn)動(dòng)的影響,主要結(jié)論如下:

    (1)研究發(fā)現(xiàn)針對(duì)傳統(tǒng)Midrex豎爐(D=5 m,H=10 m),物料在豎爐上部運(yùn)動(dòng)速率沿徑向分布較為均勻,在豎爐下部速率分布不均,從中心向爐壁處下降速率逐漸降低。

    (2)氫氣豎爐(D=5 m,H=6 m)與傳統(tǒng)Midrex豎爐相比,下層物料運(yùn)動(dòng)速率分布不均勻性更顯著,該層高度占據(jù)總料層高度的63.3%,會(huì)對(duì)實(shí)際生產(chǎn)效率和產(chǎn)品質(zhì)量帶來(lái)影響。

    可在此基礎(chǔ)上進(jìn)一步研究調(diào)整螺旋輸送機(jī)參數(shù)對(duì)改善爐內(nèi)物料運(yùn)動(dòng)分布均勻性的影響。

    [1] 殷瑞鈺. 節(jié)能、清潔生產(chǎn)、綠色制造與鋼鐵工業(yè)的可持續(xù)發(fā)展[J]. 鋼鐵, 2002, 37(08): 1-8.

    [2] 張春霞, 上官方欽, 胡長(zhǎng)慶, 等. 鋼鐵流程結(jié)構(gòu)及對(duì)CO2排放的影響[J]. 鋼鐵, 2010, 45(05): 1-6.

    [3] Sakamoto Y, Tonooka Y, Yanagisawa Y. Estimation of CO2 Emission for Each Process in the Japanese Steel Industry: a Process Analysis[J]. Energy Conveision and Management, 1999, 40(11): 1129-1140.

    [4] 方覺(jué). 非高爐煉鐵工藝與理論[M].北京:冶金工業(yè)出版社, 2002.

    [5] Costa A R D, Wagner D, Patisson F. Modelling a New, Low CO2, Emissions, Hydrogen Steelmaking Process[J]. Journal of Cleaner Production, 2013(46):27-35.

    [6] Warner N A. Towards Zero CO2 Continuous Steelmaking Directly from Ore [J]. Metallurgical and Materials Transactions B, 2014, 45(06): 2080-2096.

    [7] 徐輝, 鄒宗樹(shù), 周渝生. 豎爐生產(chǎn)直接還原鐵過(guò)程的數(shù)值模擬[J]. 材料與冶金學(xué)報(bào), 2009, 8(01):7-11.

    [8] 董雪鋒, 肖興國(guó), 鄒宗樹(shù), 等. COREX預(yù)還原反應(yīng)過(guò)程的數(shù)學(xué)物理仿真[J]. 包頭鋼鐵學(xué)院學(xué)報(bào), 1999, 18(02): 117-120.

    [9] 白明華, 葛俊禮, 龍鵠,等. 淺議我國(guó)直接還原鐵技術(shù)今后的發(fā)展[J]. 重型機(jī)械, 2013(05):1-5.

    [10] Minghua Bai, Hu Long, Ren Subo. Reduction Behavior and Kinetics of Iron Ore Pellets under H2-N2 Atmosphere[J]. ISIJ International, 2018, 58(06).

    [11] Parisi D R, Laborde M A. Modeling of Counter Current Moving Bed Gas-solid Reactor Used in Direct Reduction of Iron Ore[J]. Chemical Engineering Journal, 2004, 104(01):35-43.

    [12] 白明華, 葛俊禮, 王健,等. 淺析還原氣需求量及其與爐內(nèi)流場(chǎng)關(guān)系[J]. 重型機(jī)械, 2014(02):31-34.

    [13] Kou M Y, Wu S L, Du K P, et al. Distribution of Particle Descending Velocity in the COREX Shaft Furnace with DEM Simulation[J]. ISIJ International, 2013, 53(12):2080-2089.

    [14] Hou Q F, Samman M, Li J, et al. Modeling the Gas-solid Flow in the Reduction Shaft of COREX[J]. ISIJ International, 2014, 54(08):1772-1780.

    [15] Zhou H, Luo Z G, Zhang T, et al. Analyses of Solid Flow in COREX Shaft Furnace with AGD by Discrete Element Method[J]. Ironmaking and Steelmaking, 2015, 42(10): 774-784.

    猜你喜歡
    豎爐爐料徑向
    提高豎爐產(chǎn)能的生產(chǎn)實(shí)踐探討
    山西冶金(2022年2期)2022-06-04 10:54:24
    酒鋼高爐爐料結(jié)構(gòu)優(yōu)化系統(tǒng)設(shè)計(jì)及應(yīng)用
    淺探徑向連接體的圓周運(yùn)動(dòng)
    RN上一類Kirchhoff型方程徑向?qū)ΨQ正解的存在性
    基于PID+前饋的3MN徑向鍛造機(jī)控制系統(tǒng)的研究
    高生礦比爐料結(jié)構(gòu)生產(chǎn)實(shí)踐
    四川冶金(2019年5期)2019-12-23 09:04:56
    一類無(wú)窮下級(jí)整函數(shù)的Julia集的徑向分布
    萊鋼1880m3高爐低鎂爐料快速開(kāi)爐高效達(dá)產(chǎn)實(shí)踐
    山東冶金(2018年6期)2019-01-28 08:14:28
    COREX CGD豎爐內(nèi)煤氣流分布的數(shù)值模擬
    COREX豎爐內(nèi)部結(jié)構(gòu)對(duì)物料運(yùn)動(dòng)影響的物理模擬
    亚洲最大成人中文| 国产精品自产拍在线观看55亚洲| 成人av一区二区三区在线看| 老司机深夜福利视频在线观看| 啦啦啦观看免费观看视频高清| 欧美另类亚洲清纯唯美| 精品日产1卡2卡| av女优亚洲男人天堂| 午夜免费激情av| 搡老妇女老女人老熟妇| 性色av乱码一区二区三区2| 亚洲av成人av| 亚洲精品色激情综合| av欧美777| 啦啦啦免费观看视频1| 免费看a级黄色片| 色哟哟哟哟哟哟| 757午夜福利合集在线观看| av黄色大香蕉| 国产亚洲欧美98| 51国产日韩欧美| av女优亚洲男人天堂| 久久午夜亚洲精品久久| 天美传媒精品一区二区| 嫩草影院精品99| 97超视频在线观看视频| 99久久精品热视频| 欧美+亚洲+日韩+国产| 我的老师免费观看完整版| 国产精品亚洲一级av第二区| 亚洲无线观看免费| 午夜精品一区二区三区免费看| 日本一本二区三区精品| 观看免费一级毛片| 中出人妻视频一区二区| 非洲黑人性xxxx精品又粗又长| 久久久久性生活片| 两人在一起打扑克的视频| 国产中年淑女户外野战色| 操出白浆在线播放| 给我免费播放毛片高清在线观看| 免费大片18禁| 露出奶头的视频| 宅男免费午夜| 久久久成人免费电影| 成人av一区二区三区在线看| 欧美精品啪啪一区二区三区| 久久久久久久久中文| 国产高清三级在线| 日韩欧美国产在线观看| 成人性生交大片免费视频hd| 男人的好看免费观看在线视频| 亚洲av熟女| 欧美zozozo另类| 亚洲av成人av| 国产av不卡久久| 中国美女看黄片| 国产成人欧美在线观看| 日韩欧美在线二视频| 中文字幕人妻熟人妻熟丝袜美 | 啦啦啦韩国在线观看视频| 在线国产一区二区在线| 久久精品国产亚洲av涩爱 | 国产极品精品免费视频能看的| 亚洲人成网站在线播放欧美日韩| 国产精品99久久99久久久不卡| 午夜两性在线视频| 国产一级毛片七仙女欲春2| 好看av亚洲va欧美ⅴa在| 中文资源天堂在线| 一区二区三区国产精品乱码| 少妇的逼水好多| 国产精品乱码一区二三区的特点| 亚洲欧美日韩东京热| 老司机在亚洲福利影院| 国产精品亚洲一级av第二区| 成人av一区二区三区在线看| av天堂在线播放| 又粗又爽又猛毛片免费看| 成年女人看的毛片在线观看| 国产99白浆流出| 一区二区三区高清视频在线| 欧美+亚洲+日韩+国产| av天堂在线播放| 亚洲精品日韩av片在线观看 | 99国产极品粉嫩在线观看| 波野结衣二区三区在线 | 人妻夜夜爽99麻豆av| 成人国产一区最新在线观看| 此物有八面人人有两片| 日本与韩国留学比较| 香蕉av资源在线| 国产国拍精品亚洲av在线观看 | 日本一本二区三区精品| 好看av亚洲va欧美ⅴa在| 免费在线观看影片大全网站| 久久精品夜夜夜夜夜久久蜜豆| 宅男免费午夜| 五月玫瑰六月丁香| 欧美成人一区二区免费高清观看| 91麻豆av在线| 国产亚洲精品综合一区在线观看| 精品一区二区三区人妻视频| 国产免费av片在线观看野外av| 国产精品久久久久久人妻精品电影| 老司机午夜十八禁免费视频| 日日摸夜夜添夜夜添小说| 熟女电影av网| 国产精品久久久人人做人人爽| 观看美女的网站| 午夜激情福利司机影院| 国产在线精品亚洲第一网站| 亚洲av五月六月丁香网| 97碰自拍视频| 有码 亚洲区| 欧美色视频一区免费| 99久久成人亚洲精品观看| 成人永久免费在线观看视频| 国产高清视频在线播放一区| 精华霜和精华液先用哪个| 国产精品免费一区二区三区在线| 一个人免费在线观看电影| 夜夜夜夜夜久久久久| av黄色大香蕉| 欧美绝顶高潮抽搐喷水| 精品一区二区三区人妻视频| 国产精品一及| 一边摸一边抽搐一进一小说| 亚洲成a人片在线一区二区| 亚洲avbb在线观看| 亚洲av中文字字幕乱码综合| 欧美日韩综合久久久久久 | 欧美成狂野欧美在线观看| 精品久久久久久成人av| 波多野结衣巨乳人妻| 九九在线视频观看精品| 嫩草影视91久久| 少妇人妻一区二区三区视频| 亚洲精品国产精品久久久不卡| 18美女黄网站色大片免费观看| 97超级碰碰碰精品色视频在线观看| 久久久精品大字幕| 久久国产精品影院| 人人妻人人澡欧美一区二区| 美女黄网站色视频| 18禁美女被吸乳视频| 成人性生交大片免费视频hd| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 母亲3免费完整高清在线观看| 黄色片一级片一级黄色片| 国产综合懂色| 好看av亚洲va欧美ⅴa在| 12—13女人毛片做爰片一| 国产蜜桃级精品一区二区三区| 我要搜黄色片| 热99re8久久精品国产| 免费观看精品视频网站| 一个人看的www免费观看视频| 国产av在哪里看| 亚洲av免费在线观看| 久久6这里有精品| 日本撒尿小便嘘嘘汇集6| 国产色婷婷99| 国产精品1区2区在线观看.| 少妇人妻一区二区三区视频| 综合色av麻豆| av视频在线观看入口| 在线看三级毛片| 国产欧美日韩精品一区二区| 三级男女做爰猛烈吃奶摸视频| 国产午夜精品久久久久久一区二区三区 | 在线观看日韩欧美| 亚洲欧美精品综合久久99| 真人做人爱边吃奶动态| 999久久久精品免费观看国产| 国产一区在线观看成人免费| 成年女人永久免费观看视频| 午夜精品在线福利| 国产精品 欧美亚洲| 久久精品国产亚洲av涩爱 | 免费av观看视频| 女同久久另类99精品国产91| 美女大奶头视频| 噜噜噜噜噜久久久久久91| 99在线视频只有这里精品首页| 国产aⅴ精品一区二区三区波| 岛国在线观看网站| 久久人人精品亚洲av| 国产精品一区二区三区四区久久| 黄色片一级片一级黄色片| 亚洲欧美日韩无卡精品| 黄片小视频在线播放| 久久亚洲真实| 91麻豆av在线| 看黄色毛片网站| 亚洲成人久久性| 99热只有精品国产| 欧美成人a在线观看| 国产69精品久久久久777片| 免费在线观看日本一区| 亚洲电影在线观看av| 国产伦精品一区二区三区视频9 | 校园春色视频在线观看| 久久精品亚洲精品国产色婷小说| 一区二区三区高清视频在线| 黄片小视频在线播放| 一区二区三区免费毛片| 国产精品久久电影中文字幕| 老司机深夜福利视频在线观看| 免费看a级黄色片| 噜噜噜噜噜久久久久久91| 无人区码免费观看不卡| 精品一区二区三区av网在线观看| 午夜两性在线视频| 丰满的人妻完整版| 久久亚洲真实| 成人亚洲精品av一区二区| 国产免费男女视频| 国产三级在线视频| 午夜福利免费观看在线| 一级作爱视频免费观看| 国内精品久久久久精免费| 欧美最新免费一区二区三区 | 日本黄色视频三级网站网址| avwww免费| 两个人视频免费观看高清| 婷婷丁香在线五月| 日本熟妇午夜| 国产成人系列免费观看| 国产精品女同一区二区软件 | 久久久久久久久久黄片| 桃红色精品国产亚洲av| 亚洲人成网站高清观看| 亚洲 欧美 日韩 在线 免费| 91麻豆精品激情在线观看国产| 18美女黄网站色大片免费观看| 国内久久婷婷六月综合欲色啪| 免费无遮挡裸体视频| 日韩 欧美 亚洲 中文字幕| 观看美女的网站| 亚洲欧美精品综合久久99| 亚洲va日本ⅴa欧美va伊人久久| 亚洲乱码一区二区免费版| 国产一区二区三区在线臀色熟女| 无人区码免费观看不卡| 99久久精品一区二区三区| 两个人的视频大全免费| 亚洲熟妇熟女久久| 欧美日韩福利视频一区二区| 啦啦啦观看免费观看视频高清| 日韩欧美 国产精品| 国产一区二区亚洲精品在线观看| 一进一出抽搐动态| 中文字幕久久专区| 色哟哟哟哟哟哟| 九色成人免费人妻av| 99热6这里只有精品| 成人一区二区视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 日韩欧美一区二区三区在线观看| 日本 欧美在线| 久久伊人香网站| 欧美日韩综合久久久久久 | 内地一区二区视频在线| 白带黄色成豆腐渣| 搡女人真爽免费视频火全软件 | 在线观看午夜福利视频| 国产高清videossex| 欧美三级亚洲精品| 五月玫瑰六月丁香| 欧美日韩福利视频一区二区| 一区二区三区国产精品乱码| 亚洲真实伦在线观看| 一夜夜www| 校园春色视频在线观看| 变态另类成人亚洲欧美熟女| 他把我摸到了高潮在线观看| 欧美xxxx黑人xx丫x性爽| 天堂网av新在线| 亚洲成人久久性| 十八禁网站免费在线| 日韩有码中文字幕| 一区二区三区国产精品乱码| 日本五十路高清| 国产一区二区激情短视频| 最近在线观看免费完整版| 一区二区三区国产精品乱码| 午夜影院日韩av| 日韩av在线大香蕉| 真实男女啪啪啪动态图| 国产高清视频在线播放一区| 色综合站精品国产| 夜夜夜夜夜久久久久| 欧美不卡视频在线免费观看| avwww免费| 男人的好看免费观看在线视频| 免费在线观看日本一区| 午夜福利在线在线| 观看美女的网站| 日本在线视频免费播放| 国产精品亚洲一级av第二区| 久久久久久久久久黄片| 午夜精品久久久久久毛片777| a级一级毛片免费在线观看| 啦啦啦免费观看视频1| 性色avwww在线观看| 免费人成在线观看视频色| 99国产极品粉嫩在线观看| 男女之事视频高清在线观看| 国产v大片淫在线免费观看| 欧美日韩综合久久久久久 | 色尼玛亚洲综合影院| 老熟妇乱子伦视频在线观看| 丰满人妻一区二区三区视频av | 亚洲午夜理论影院| 啦啦啦免费观看视频1| 亚洲精品久久国产高清桃花| 久久久精品大字幕| 日本五十路高清| 啪啪无遮挡十八禁网站| 男女那种视频在线观看| 中文字幕av成人在线电影| 一夜夜www| 最新美女视频免费是黄的| 日本免费一区二区三区高清不卡| 极品教师在线免费播放| 亚洲精品乱码久久久v下载方式 | 午夜精品在线福利| 午夜免费男女啪啪视频观看 | 精品久久久久久久毛片微露脸| www国产在线视频色| 亚洲中文字幕日韩| 国产精品爽爽va在线观看网站| 少妇的丰满在线观看| 又紧又爽又黄一区二区| 老汉色av国产亚洲站长工具| 免费在线观看亚洲国产| 国产精品国产高清国产av| 色综合亚洲欧美另类图片| 午夜精品久久久久久毛片777| 国产伦一二天堂av在线观看| 九九久久精品国产亚洲av麻豆| 麻豆国产av国片精品| 日本熟妇午夜| 亚洲国产欧洲综合997久久,| 国产视频内射| 日本成人三级电影网站| 国产精品日韩av在线免费观看| 欧美极品一区二区三区四区| 日本黄大片高清| 欧美性猛交╳xxx乱大交人| 国产精品 欧美亚洲| 国产爱豆传媒在线观看| 成人午夜高清在线视频| 国产 一区 欧美 日韩| 亚洲精品在线观看二区| 男女视频在线观看网站免费| 免费一级毛片在线播放高清视频| 91av网一区二区| 一本久久中文字幕| 欧美乱码精品一区二区三区| 国产精品香港三级国产av潘金莲| 99在线视频只有这里精品首页| 日韩国内少妇激情av| 日韩欧美三级三区| e午夜精品久久久久久久| 国产精品亚洲一级av第二区| 淫秽高清视频在线观看| 亚洲国产精品久久男人天堂| 精品人妻偷拍中文字幕| 丝袜美腿在线中文| 成人欧美大片| 国产精品香港三级国产av潘金莲| 1024手机看黄色片| 一个人看视频在线观看www免费 | 免费观看精品视频网站| 男女下面进入的视频免费午夜| 别揉我奶头~嗯~啊~动态视频| 日韩精品中文字幕看吧| 久久精品综合一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 女人高潮潮喷娇喘18禁视频| 成年人黄色毛片网站| 久久久久国产精品人妻aⅴ院| 国产三级黄色录像| 99久久成人亚洲精品观看| 亚洲乱码一区二区免费版| 男人和女人高潮做爰伦理| 久久久久国产精品人妻aⅴ院| 欧美一区二区国产精品久久精品| 日韩欧美国产在线观看| 欧美成人性av电影在线观看| 亚洲在线自拍视频| 国产主播在线观看一区二区| 别揉我奶头~嗯~啊~动态视频| 熟妇人妻久久中文字幕3abv| 97超视频在线观看视频| 亚洲欧美精品综合久久99| 亚洲成人久久性| 18禁美女被吸乳视频| a级一级毛片免费在线观看| avwww免费| 99热6这里只有精品| 精品一区二区三区人妻视频| 深爱激情五月婷婷| 啦啦啦韩国在线观看视频| 午夜福利欧美成人| 亚洲人与动物交配视频| 少妇高潮的动态图| 深爱激情五月婷婷| 国产aⅴ精品一区二区三区波| 最好的美女福利视频网| 男人的好看免费观看在线视频| 午夜日韩欧美国产| 国产精品一及| 两人在一起打扑克的视频| 国产视频一区二区在线看| 俄罗斯特黄特色一大片| 国产成人av激情在线播放| 在线视频色国产色| 精品欧美国产一区二区三| 日韩国内少妇激情av| 99精品欧美一区二区三区四区| 午夜福利成人在线免费观看| 12—13女人毛片做爰片一| 亚洲国产精品成人综合色| 免费人成在线观看视频色| 国产午夜精品久久久久久一区二区三区 | 亚洲五月婷婷丁香| 一本综合久久免费| 香蕉av资源在线| 性欧美人与动物交配| av在线蜜桃| 91在线精品国自产拍蜜月 | 国产精品一及| 又爽又黄无遮挡网站| 精品人妻偷拍中文字幕| 欧美乱妇无乱码| 国产精品亚洲美女久久久| 国产一区二区亚洲精品在线观看| 老鸭窝网址在线观看| 尤物成人国产欧美一区二区三区| 日韩精品青青久久久久久| 一区二区三区免费毛片| 一个人看视频在线观看www免费 | 免费观看精品视频网站| 久久这里只有精品中国| 美女 人体艺术 gogo| 亚洲精品在线美女| 观看美女的网站| 美女被艹到高潮喷水动态| 日本 av在线| 老司机午夜十八禁免费视频| 三级国产精品欧美在线观看| 观看美女的网站| 日韩欧美在线乱码| 男女床上黄色一级片免费看| 亚洲不卡免费看| 久久人人精品亚洲av| 久9热在线精品视频| 三级男女做爰猛烈吃奶摸视频| 噜噜噜噜噜久久久久久91| 天堂网av新在线| 亚洲美女视频黄频| 国产视频内射| 伊人久久大香线蕉亚洲五| 俺也久久电影网| 欧美激情久久久久久爽电影| 听说在线观看完整版免费高清| 亚洲在线观看片| 欧美午夜高清在线| 亚洲avbb在线观看| 此物有八面人人有两片| 精品一区二区三区视频在线 | 中国美女看黄片| 又黄又粗又硬又大视频| 国产精品爽爽va在线观看网站| 欧美一级a爱片免费观看看| www.色视频.com| 99久国产av精品| 精品人妻1区二区| 国产美女午夜福利| 日本成人三级电影网站| 夜夜爽天天搞| 熟妇人妻久久中文字幕3abv| 露出奶头的视频| 精品国内亚洲2022精品成人| 91麻豆精品激情在线观看国产| 亚洲人成网站在线播| 午夜免费成人在线视频| 国产真实乱freesex| 久久国产精品影院| 91久久精品国产一区二区成人 | 在线视频色国产色| 亚洲欧美激情综合另类| 人妻夜夜爽99麻豆av| 久久性视频一级片| 国产精品一区二区三区四区免费观看 | 久久这里只有精品中国| 成年女人毛片免费观看观看9| 久久欧美精品欧美久久欧美| 操出白浆在线播放| 日本黄色视频三级网站网址| 亚洲精品久久国产高清桃花| 婷婷亚洲欧美| 亚洲欧美一区二区三区黑人| 欧美日韩国产亚洲二区| 亚洲av日韩精品久久久久久密| 99精品在免费线老司机午夜| 中文字幕人妻丝袜一区二区| 尤物成人国产欧美一区二区三区| av欧美777| 欧美极品一区二区三区四区| 色在线成人网| e午夜精品久久久久久久| 在线观看av片永久免费下载| 精品熟女少妇八av免费久了| 成年女人看的毛片在线观看| 一本久久中文字幕| 久久香蕉精品热| 18美女黄网站色大片免费观看| 老熟妇乱子伦视频在线观看| 国产高清有码在线观看视频| 国产精品亚洲美女久久久| 偷拍熟女少妇极品色| 国产激情欧美一区二区| 国产黄a三级三级三级人| 国产又黄又爽又无遮挡在线| 搡老岳熟女国产| 中出人妻视频一区二区| 日韩欧美精品v在线| 国产单亲对白刺激| 国产野战对白在线观看| 成年免费大片在线观看| 午夜福利在线在线| 亚洲欧美激情综合另类| 欧美又色又爽又黄视频| 在线免费观看不下载黄p国产 | xxxwww97欧美| 人人妻人人澡欧美一区二区| av片东京热男人的天堂| 女人高潮潮喷娇喘18禁视频| 国产视频内射| 午夜免费成人在线视频| xxxwww97欧美| 国产野战对白在线观看| 久久草成人影院| 搡老岳熟女国产| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利欧美成人| 色av中文字幕| 搡老岳熟女国产| 欧美丝袜亚洲另类 | 别揉我奶头~嗯~啊~动态视频| 成年女人看的毛片在线观看| tocl精华| 每晚都被弄得嗷嗷叫到高潮| 亚洲av免费高清在线观看| 身体一侧抽搐| 中国美女看黄片| 午夜两性在线视频| 精品欧美国产一区二区三| 免费在线观看影片大全网站| 女警被强在线播放| 伊人久久大香线蕉亚洲五| 国产麻豆成人av免费视频| 成人国产一区最新在线观看| 亚洲人成伊人成综合网2020| 男女午夜视频在线观看| 久久精品国产99精品国产亚洲性色| 亚洲精品日韩av片在线观看 | 国产探花极品一区二区| 亚洲国产精品成人综合色| 内射极品少妇av片p| 免费人成视频x8x8入口观看| 国产精品1区2区在线观看.| 18禁裸乳无遮挡免费网站照片| 夜夜看夜夜爽夜夜摸| 99在线视频只有这里精品首页| 大型黄色视频在线免费观看| 日本撒尿小便嘘嘘汇集6| 日韩成人在线观看一区二区三区| 人妻久久中文字幕网| 午夜a级毛片| 久久人妻av系列| 国产一区二区激情短视频| 18+在线观看网站| 69人妻影院| 欧美激情久久久久久爽电影| 成人无遮挡网站| 久久精品91蜜桃| 18禁黄网站禁片午夜丰满| 男女午夜视频在线观看| 人妻久久中文字幕网| 精品国产美女av久久久久小说| 国产精品一区二区免费欧美| 有码 亚洲区| 亚洲国产日韩欧美精品在线观看 | 国产亚洲精品av在线| 51午夜福利影视在线观看| 色老头精品视频在线观看| 嫁个100分男人电影在线观看| 12—13女人毛片做爰片一| 亚洲欧美激情综合另类| 一级a爱片免费观看的视频| 国产高清videossex| 九色成人免费人妻av| 久久久久免费精品人妻一区二区| 女人十人毛片免费观看3o分钟| 小说图片视频综合网站| 精品熟女少妇八av免费久了| 国产午夜精品久久久久久一区二区三区 | 免费看美女性在线毛片视频| 国内精品一区二区在线观看| 毛片女人毛片|