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

    基于離散單元法模擬引入AGD技術(shù)COREX豎爐物料運(yùn)動行為

    2016-05-10 03:59:55羅志國鄒宗樹
    材料與冶金學(xué)報(bào) 2016年1期
    關(guān)鍵詞:法向力

    游 洋,周 恒,張 濤,羅志國,鄒宗樹

    (東北大學(xué)多金屬共生礦生態(tài)化冶金教育部重點(diǎn)實(shí)驗(yàn)室,沈陽110819)

    ?

    基于離散單元法模擬引入AGD技術(shù)COREX豎爐物料運(yùn)動行為

    游洋,周恒,張濤,羅志國,鄒宗樹

    (東北大學(xué)多金屬共生礦生態(tài)化冶金教育部重點(diǎn)實(shí)驗(yàn)室,沈陽110819)

    摘要:基于離散單元法,建立引入Areal Gas Distribution( AGD)技術(shù)COREX豎爐物料運(yùn)動行為的DEM模型.模型計(jì)算結(jié)果與文獻(xiàn)物理實(shí)驗(yàn)結(jié)果一致,證明模型的合理性.模擬結(jié)果顯示,AGD豎爐內(nèi)存在三種類型的流動區(qū)域:活塞流區(qū)、準(zhǔn)停滯區(qū)及溝流區(qū).物料運(yùn)動流型呈現(xiàn)Flat→波浪→W的演變過程.AGD梁下方形成的三角形空隙是還原煤氣導(dǎo)入豎爐中心的主要煤氣通道.豎爐圍管區(qū)域AGD梁的安裝會影響該區(qū)域物料的均勻下降,局部區(qū)域的較大法向力可能導(dǎo)致物料的擠壓黏結(jié),誘發(fā)圍管slot堵塞.

    關(guān)鍵詞:COREX豎爐; AGD;離散單元;物料運(yùn)動;法向力

    高爐煉鐵是現(xiàn)代冶金工業(yè)生產(chǎn)的主導(dǎo)流程,但其面臨著環(huán)境污染嚴(yán)重,焦煤資源逐漸短缺的問題.為減少煉焦過程的污染排放及擺脫對冶金焦炭的依賴,近年來諸如COREX,F(xiàn)INEX,Hismelt以及HIsarna等非高爐煉鐵技術(shù)得到迅猛發(fā)展[1,2],其中COREX是首先實(shí)現(xiàn)工業(yè)化運(yùn)行的一種用煤和球團(tuán)(塊礦)生產(chǎn)鐵水的煉鐵新工藝,具有無焦或少焦冶煉的特點(diǎn)[3,4].COREX上部豎爐及下部熔化氣化爐均為典型的氣固逆流反應(yīng)器,離散顆粒的運(yùn)動行為直接影響爐況的穩(wěn)定和順行,因此,COREX過程中固體物料的運(yùn)動現(xiàn)象一直是學(xué)者研究的熱點(diǎn)[5~7].寶鋼投產(chǎn)的第二套COREX系統(tǒng)在第一套的基礎(chǔ)上作了大量改進(jìn),最大變化是引入了Areal Gas Distribution ( AGD)技術(shù),即在豎爐圍管區(qū)域加裝AGD梁,利用物料下行運(yùn)動過程中于梁下方自發(fā)形成的無爐料煤氣通道將還原煤氣引入豎爐中心,改善煤氣流分布.實(shí)際生產(chǎn)表明,AGD梁的安裝有利于豎爐中心煤氣發(fā)展,提高了爐內(nèi)金屬化率,但亦改變了豎爐內(nèi)部結(jié)構(gòu),直接影響爐內(nèi)物料運(yùn)動行為.因此,有必要對引進(jìn)AGD技術(shù)COREX豎爐物料運(yùn)動行為進(jìn)行研究.

    煉鐵反應(yīng)器中物料運(yùn)動行為研究方法可分為兩大類:物理實(shí)驗(yàn)研究和數(shù)值模擬研究.前人建立二維或三維冷態(tài)模型對高爐中的物料運(yùn)動行為進(jìn)行了詳細(xì)研究,但關(guān)于COREX豎爐的研究較少.Lee[8]建立冷態(tài)模型研究了有無死料柱對豎爐物料運(yùn)動流型的影響.Zhou等[9]建立三維半周物理模型對引入AGD技術(shù)COREX豎爐物料運(yùn)動行為進(jìn)行了研究并綜合分析了不同操作條件下有無AGD梁爐內(nèi)物料運(yùn)動行為的變化[10].上述研究成果為了解COREX豎爐中物料運(yùn)動行為提供了有效信息,但冷態(tài)實(shí)驗(yàn)結(jié)果只能反映爐料運(yùn)動的宏觀行為,不能從微觀角度描述物料的運(yùn)動特征.近年來,基于拉格朗日坐標(biāo)系考慮顆粒與顆粒/壁面相互作用的離散單元法( DEM)可從顆粒尺度描述物料運(yùn)動行為,無須過多假設(shè),并可獲得微觀的瞬態(tài)信息,已引起學(xué)者們的廣泛關(guān)注.事實(shí)上,DEM方法已被應(yīng)用到模擬COREX豎爐頂部布料[11,12],下部物料下降速度[13],爐內(nèi)非對稱流動現(xiàn)象[14],受黏性力影響的物料運(yùn)動行為[15]等,上述工作主要為傳統(tǒng)豎爐內(nèi)物料運(yùn)動信息,關(guān)于引入AGD技術(shù)COREX豎爐較少.Hou[15]建立DEM模型模擬了不同爐型結(jié)構(gòu)下物料運(yùn)動狀態(tài),關(guān)注AGD梁對物料宏觀下降速度的影響.

    本文基于離散單元法研究AGD技術(shù)COREX豎爐物料運(yùn)動跡線、時(shí)間線分布以及AGD對爐內(nèi)物料運(yùn)動流型及應(yīng)力分布的影響,為進(jìn)一步了解爐內(nèi)物料運(yùn)動行為提供理論依據(jù).

    1 數(shù)學(xué)模型

    1.1 DEM模型建立

    離散物料運(yùn)動過程中受自身重力及顆粒與顆粒/壁面之間的作用力,同時(shí)顆粒還受到切向力扭矩和滾動摩擦扭矩作用,使顆粒產(chǎn)生平移運(yùn)動和轉(zhuǎn)動.根據(jù)牛頓第二定律,顆粒運(yùn)動的控制方程可描述為

    式中,mi,Ii,vi和ωi分別代表顆粒i的質(zhì)量( kg)、轉(zhuǎn)動慣量( kg·m2)、平動速度( m/s)和轉(zhuǎn)動速度( r/s) ; mig為i顆粒所受重力; Fcn,ij,F(xiàn)ct,ij為顆粒間法向、切向接觸力( N) ; Fdn,ij,F(xiàn)dt,ij為顆粒間法向、切向阻尼力( N) ; Tij,Mij分別為顆粒所受切向力扭矩和滾動摩擦扭矩( N·m) ; ki為與i顆粒相接觸顆粒數(shù)目.本文中顆粒間相互作用力模型采用線性“彈簧-阻尼-摩擦板”模型.顆粒所受接觸力、阻尼力、摩擦力及扭矩詳見文獻(xiàn)[5,7,15].

    1.2模擬條件

    本研究采用自主開發(fā)模擬程序進(jìn)行計(jì)算,整個(gè)代碼采用面向?qū)ο笳Z言C + +實(shí)現(xiàn).為減小計(jì)算負(fù)擔(dān),本研究將AGD技術(shù)COREX-3000豎爐實(shí)際幾何尺寸按1∶20縮小,圖1為用于本研究豎爐計(jì)算區(qū)域及幾何尺寸.本研究二維扁平模型寬度取5個(gè)dp( 0.05 m),排料速度設(shè)定每100時(shí)間步長內(nèi)有4個(gè)顆粒從一側(cè)螺旋排出,排料速度約為0.251 kg/s.實(shí)際生產(chǎn)中,COREX豎爐頂部加入球團(tuán)、塊礦、焦炭和熔劑等物料,而球團(tuán)加入量大于其他物料量,可用單一尺寸顆粒表征物料運(yùn)動流型[14,15].因此,本文選取顆粒為單一球形顆粒,表征爐內(nèi)球團(tuán).具體模擬條件和模擬參數(shù)如表1所示.

    圖1 計(jì)算COREX豎爐模型尺寸圖(單位: mm)Fig.1 Schematic diagram of calculation region of COREX shaft furnace ( Unit: mm)

    表1 模擬條件及參數(shù)Table1 Particle properties and simulation conditions

    2 結(jié)果與討論

    2.1模型驗(yàn)證

    為驗(yàn)證模型的正確性,本文將數(shù)值計(jì)算結(jié)果與文獻(xiàn)[8]中的二維物理實(shí)驗(yàn)結(jié)果進(jìn)行了對比.其中,物理實(shí)驗(yàn)中物料為磁鐵礦,為驗(yàn)證豎爐內(nèi)的流動行為,需對球形顆粒物性參數(shù)進(jìn)行修正,特別是關(guān)乎顆粒運(yùn)動及表征非球形度的摩擦系數(shù)進(jìn)行調(diào)整.驗(yàn)證過程中,取顆粒滑動摩擦系數(shù)μs,pp=0.70,滾動摩擦系數(shù)μr,pp= 0.01d,模擬得到顆粒的安息角為40.89(°),如圖2( a)所示,與文獻(xiàn)[8]中磁鐵礦安息角41(°)相近,可用于模擬物理實(shí)驗(yàn)顆粒.由于本模型為二維扁平模型,模型寬度為5個(gè)dp( 0.05 m),底部排料出口基于螺旋排料特性,顆粒主要從螺旋絲杠頂端排出,排料速度設(shè)定為每100時(shí)間步長內(nèi)有4個(gè)顆粒從一側(cè)螺旋排出,待爐內(nèi)物料運(yùn)動達(dá)到穩(wěn)定狀態(tài)時(shí),顆粒爐內(nèi)運(yùn)動狀態(tài)如圖3所示.計(jì)算所得物料運(yùn)動流型及死料區(qū)等特征流動區(qū)域與實(shí)驗(yàn)測量值基本一致,說明了數(shù)學(xué)模型的合理性.本文模擬COREX -3000豎爐中物料運(yùn)動行為,球團(tuán)設(shè)置滑動摩擦系數(shù)μs,pp= 0.50,滾動摩擦系數(shù)μr,pp= 0.005 d,模擬得到安息角為32.8(°),見圖2( b),近似于實(shí)際球團(tuán)礦休止角( 32(°) ),可用于本表征豎爐內(nèi)物料運(yùn)動行為.

    圖2 安息角Fig.2 Angle of repose( a)—μs,pp=0.70,μr,pp=0.01 d;( b)—μs,pp=0.50,μr,pp=0.005 d

    圖3 模擬結(jié)果與實(shí)驗(yàn)結(jié)果對比Fig.3 Comparison between experimental and numerical results( a)—模擬結(jié)果; ( b)—物理實(shí)驗(yàn)[8]

    2.2物料運(yùn)動流型

    COREX豎爐物料堆積到一定高度后連續(xù)排料,當(dāng)爐內(nèi)運(yùn)動狀態(tài)達(dá)到穩(wěn)定狀態(tài)時(shí),在爐頂加入紅色示蹤料層,表征爐內(nèi)的運(yùn)動模式.圖4為排料速度0.251 kg/s下,引入AGD技術(shù)COREX豎爐物料運(yùn)動流型.由圖可知,達(dá)到穩(wěn)定狀態(tài)時(shí)爐內(nèi)存在三種類型的流動區(qū)域:活塞流區(qū)( I)、溝流區(qū)( II)以及準(zhǔn)停滯區(qū)( III).豎爐實(shí)際生產(chǎn)過程中,為避免底部物料的無滑移運(yùn)動,于底部中心人為添加導(dǎo)流錐,本模型亦在其底部中心區(qū)域設(shè)置導(dǎo)流錐( IV).豎爐上部為典型的活塞流區(qū),物料在徑向上均勻下降,提供良好的氣固還原環(huán)境.溝流區(qū)( II)位于豎爐底部螺旋排料絲杠頂端上方區(qū)域,該區(qū)域物料下降速度較快.而豎爐內(nèi)的準(zhǔn)停滯區(qū)( III)位于豎爐下部螺旋絲桿近端上方.這一流動模式特別是豎爐下部的流動狀態(tài)與冷態(tài)模擬實(shí)驗(yàn)確定的流動特征一致[9,10].其中,準(zhǔn)停滯區(qū)的形成受豎爐下部螺旋排料的影響.爐料在重力作用下向螺旋絲桿頂端空隙運(yùn)動,該處有最大下降速度,隨著螺旋絲杠的機(jī)械作用,頂端物料水平向外輸送,從而使絲杠近端螺齒空隙被頂端傳送過來物料填滿,導(dǎo)致螺旋絲杠近端上方物料下降速度變小.本研究設(shè)定螺旋絲杠頂端為主要排料區(qū)域,因此在螺旋絲杠近端上方形成爐料緩慢移動的準(zhǔn)停滯區(qū).另外,模擬發(fā)現(xiàn),顆粒運(yùn)動通過AGD梁時(shí),會于梁下方自發(fā)形成三角形煤氣通道,該現(xiàn)象在物理實(shí)驗(yàn)[9,10]中亦可觀察得到.還原煤氣可從豎爐圍管區(qū)域的AGD煤氣入口鼓入爐內(nèi)并沿該三角形煤氣通道達(dá)到豎爐中心,改善傳統(tǒng)COREX-3000豎爐由于爐徑擴(kuò)大造成的還原煤氣難以到達(dá)豎爐中心的問題,其工作原理可如圖5所示.

    圖4 物料運(yùn)動流型Fig.4 Snapshots of solid flow pattern from initial to steady state

    圖6為引入AGD技術(shù)COREX豎爐物料運(yùn)動跡線及時(shí)間線.顆粒運(yùn)動跡線顯示,AGD梁正上方示蹤顆粒通過AGD梁時(shí)跡線發(fā)生變化,顆粒繞梁向下運(yùn)動.其他區(qū)域示蹤顆粒則在豎爐上部沿豎直方向向下運(yùn)動,進(jìn)入豎爐下部后向豎爐中間的排料出口偏移.物料運(yùn)動流型顯示,豎爐上部物料均勻下降,運(yùn)動流型呈現(xiàn)‘一’型.物料進(jìn)入圍管區(qū)域后,受AGD梁的阻礙作用,梁上方物料下降速度變慢,爐內(nèi)中心及邊緣區(qū)域下降速度快,物料運(yùn)動流型呈現(xiàn)波浪型.進(jìn)入豎爐下部后,受排料出口作用,運(yùn)動流型呈現(xiàn)‘W’型.需要指出的是,豎爐下部螺旋結(jié)構(gòu)直接影響物料運(yùn)動流型,前人[13]模擬發(fā)現(xiàn),螺旋絲杠頂端螺徑稍有偏大則會引起絲杠頂端物料下降速度加快,螺旋絲杠近端上方物料下降速度減小甚至趨于零,說明下部的流型易出現(xiàn)‘W’型,而優(yōu)化螺旋絲杠頂端螺徑可獲得較均勻的物料下降速度.事實(shí)上,漸增型螺距、錐形螺徑、反錐形螺桿以及它們的相互組合均可增大沿螺桿方向排料能力,獲得均勻的物料下降模式.關(guān)于螺旋結(jié)構(gòu)對COREX豎爐物料運(yùn)動行為影響,特別是螺旋結(jié)構(gòu)與不同AGD梁布置方案的相關(guān)匹配工作正在進(jìn)行中,將在后續(xù)報(bào)道中詳細(xì)呈現(xiàn).

    2.3 AGD對物料運(yùn)動的影響

    圖7為AGD對COREX豎爐物料運(yùn)動流型的影響.由圖可知,無AGD的傳統(tǒng)豎爐內(nèi),圍管及其以上區(qū)域?yàn)槲锪暇鶆蛳陆档幕钊鲄^(qū),物料運(yùn)動流型為‘一’型.隨著物料繼續(xù)向下運(yùn)動,爐墻處物料下降速度變慢,在圍管下方運(yùn)動流型轉(zhuǎn)變?yōu)椤甎’型.進(jìn)入豎爐下部后,受排料出口影響,物料運(yùn)動流型變?yōu)椤甒’型.而加裝AGD梁后,受梁阻礙作用,物料在梁上方流型為波浪型且穿過AGD梁后物料呈現(xiàn)‘W’型的雛形,物料在圍管區(qū)域徑向下降速度發(fā)生變化,這說明AGD梁的安裝會影響豎爐圍管區(qū)域物料的均勻下降.AGD對物料運(yùn)動行為的影響將進(jìn)一步通過顆粒間微觀信息分布進(jìn)行研究.

    圖5 AGD原理Fig.5 Working principle of AGD

    圖6 物料運(yùn)動跡線及時(shí)間線ig.6 Streamlines and timelines of burden

    圖7 AGD對物料運(yùn)動流型影響Fig.7 Influence of AGD on solid flow pattern( a)—Without AGD; ( b)—WithAGD

    圖8 AGD對顆粒法向力分布影響Fig.8 Influence of AGD on normal force distribution( a)—Without AGD; ( b)—WithAGD

    圖8為豎爐內(nèi)顆粒法向力分布.由圖可知,加裝AGD梁后圍管區(qū)域應(yīng)力分布出現(xiàn)兩大變化:一是AGD梁上方區(qū)域顆粒具有較大法向力,主要原因是該區(qū)域顆粒坐落梁頂部支撐著豎爐上部顆粒;二是梁正下方區(qū)域顆粒法向力較小,這主要是由于梁的安裝使其正下方顆粒不受上部顆粒的擠壓所致.圖9為顆粒法向力概率密度函數(shù)分布,其中圖9( a)為圍管區(qū)域顆粒,圖9( b)為圍管以下區(qū)域.由圖9( a)可知,隨著AGD梁的安裝,圍管區(qū)域顆粒法向力概率密度分布曲線頂點(diǎn)向右移動,說明圍管區(qū)域具備較高法向力顆粒數(shù)量增多.圖9( b)顯示,受AGD梁對豎爐上部物料的支撐作用,圍管以下區(qū)域顆粒法向力概率密度分布曲線向左移動,顆粒間的法向力減小.需要指出的是,圍管局部區(qū)域法向應(yīng)力增大易導(dǎo)致物料相互擠壓,在高溫發(fā)黏條件下易導(dǎo)致物料間的黏結(jié)結(jié)塊,這些都是誘發(fā)AGD豎爐圍管區(qū)域slot堵塞的原因之一.

    圖9 顆粒法向力概率密度函數(shù)分布Fig.9 Probability density distribution of the normal contact force in bustle zone( a)—圍管區(qū); ( b)—圍管下部

    3 結(jié)論

    本文采用自編程離散單元程序?qū)σ階real Gas Distribution ( AGD)技術(shù)COREX豎爐物料運(yùn)動行為進(jìn)行模擬.模型計(jì)算結(jié)果與文獻(xiàn)實(shí)驗(yàn)結(jié)果吻合良好,證明了模型的合理性.

    AGD技術(shù)COREX豎爐物料流動區(qū)域可分為活塞流區(qū)、溝流區(qū)以及準(zhǔn)停滯區(qū).爐內(nèi)物料運(yùn)動流型呈現(xiàn)‘一’→‘波浪’→‘W’的演變過程.物料通過AGD梁時(shí)在其下方形成的三角形空隙區(qū)域是還原煤氣導(dǎo)入豎爐中心的主要煤氣通道.與無AGD技術(shù)的傳統(tǒng)豎爐物料運(yùn)動流型對比發(fā)現(xiàn),AGD梁的安裝會破壞圍管區(qū)域物料下降的均勻性.顆粒間法向力分布顯示圍管區(qū)域存在局部的較大應(yīng)力區(qū),易導(dǎo)致該區(qū)物料的擠壓黏結(jié),誘發(fā)AGD豎爐圍管區(qū)域slot的堵塞.

    參考文獻(xiàn):

    [1]周渝生.煤基熔融還原煉鐵新工藝開發(fā)現(xiàn)狀評述[J].鋼鐵,2005,40( 11) : 1-8.( Zhou Yusheng.Review of current development of coalbased smelting reduction Ironmaking process[J].Iron and Steel,2005,40( 11) : 1-8.)

    [2]Qu Y X,Yang Y X,Zou Z S,et al.Thermal decomposition behaviour of fine iron ore particles[J].ISIJ International,2014,54( 10) : 2196-2205.

    [3]Wang N,Xie X M,Zou Z S,et al.Analysis of material and energy consumption of Corex C3000[J].Steel Research international,2008,79( 7) : 547-552.

    [4]Qu Y X,Zou Z S,Xiao Y P.A comprehensive static model for COREX process[J].ISIJ International,2012,52( 12) : 2186-2193.

    [5]Li H F,Luo Z G,Zou Z S,et al.Mathematical Simulation of Burden Distribution in COREX Melter gasifier by Discrete Element Method[J].Journal of Iron and Steel Research,International,2012,19( 9) : 36-42.

    [6]Sun J J,Luo Z G,Zou Z S.Numerical simulation of raceway phenomena in a COREX melter-gasifier[J].Powder Technology,2105,281: 159-166.

    [7]Han L H,Luo Z G,Zhou H,et al.Influence of cohesive zone shape on solid flow in COREX melter gasifier by discrete element method[J].Journal of Iron and Steel Research,International,2015,22( 4) 304-310.

    [8]Lee Y.A scaled model study on the solid flow in a shaft type furnace[J].Powder Technology,1999,102( 2) : 194-201.

    [9]Zhou H,Zou Z S,Luo Z G,et al.Analyses of solid flow in latest design COREX shaft furnace by physical simulation[J].Ironmaking Steelmaking,2015,42( 3) : 209-216.

    [10]Zhou H,Luo Z G,Zou Z S,et al.Experimental study on burden descending behavior in COREX shaft furnace with AGD beams[J].Steel Research international,2015,86( 9) : 1073-1081.

    [11]Li Q,F(xiàn)eng M X,Zou Z S.Validation and calibration approach for discrete element simulation of burden charging in pre-reduction shaft furnace of COREX process[J].ISIJ International,2013,53( 8) : 1365-1371.

    [12]Kou M Y,Wu S L,Wang G,et al.Numerical simulation of burden and gas distribution inside COREX shaft furnace[J].Steel Research International,2015,86 ( 6) 686-694.

    [13]Kou M Y,Wu S L,Shen W,et al.Distribution of particle descending velocity in the COREX shaft furnace with DEM simulation[J].ISIJ International,2013,53 ( 12 ) : 2080 -2089.

    [14]Luo Z G,Zhou H,Zhang T,et al.DEM simulation of solid flow including asymmetric phenomena in a model COREX shaft furnace[J].Journal of Iron and Steel Research,International,2014,In press.

    [15]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( 8) : 1772-1780.

    Modeling of solid flow in COREX shaft furnace with AGD introduced by discrete element method

    You Yang,Zhou Heng,Zhang Tao,Luo Zhiguo,Zou Zongshu
    ( Key Laboratory of Ecological Metallurgy of Multimetallic Minerals of Education Ministry,Northeastern University,Shenyang 110819,China)

    Abstract:Based on the principle of discrete element method ( DEM ),a mathematical model is developed to numerical investigate the solid flow behaviour in COREX shaft furnace ( SF) with a new technique called Areal Gas Distribution ( AGD) introduced.The applicability of the DEM model is validated from its good agreement with the previous experiment in terms of solid flow patterns.The results confirm that the solid flow in COREX SF with AGD beams can be divided into three different flow regions: plug flow zone,quasi-stagnant zone and converging flow region.The basic solid flow profile turns out to be clear Flat→Wave→W type.A triangle shaped free area is observed under the AGD beam,which is the main channel for gas flow into the shaft centre.The effect of AGD beams on solid flow is evaluated and the investigation reveals that the AGD beams affect the uniform descending behaviour of particle in bustle zone.The large normal contact force upstream the beam could possible increase the period of static contacts and the related sticking effect.As a result,the choking of gas slots may be induced.

    Key words:COREX shaft furnace; AGD; discrete element method; solid flow; normal contact force

    通訊作者:鄒宗樹( 1958—),男,教授,博士生導(dǎo)師,E-mail: zouzs@ mail.neu.edu.com.

    作者簡介:游洋( 1989—),男,博士研究生.

    基金項(xiàng)目:國家科技支撐計(jì)劃項(xiàng)目( 2011BAE04B02) ;國家自然科學(xué)基金資助項(xiàng)目( 51174053).

    收稿日期:2015-09-02.

    doi:10.14186/j.cnki.1671-6620.2016.01.002

    中圖分類號:TF 557

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1671-6620( 2016) 01-0006-06

    猜你喜歡
    法向力
    桿式天平軸向力元件測量梁夾角優(yōu)化分析
    削弱永磁直線電機(jī)法向力波動優(yōu)化設(shè)計(jì)
    軌道交通用直線電動機(jī)的電磁力特性分析
    基于動力學(xué)分析的大重合度直齒圓柱齒輪強(qiáng)度計(jì)算
    基于Stick-slip振動試驗(yàn)的Stribeck摩擦模型參數(shù)研究
    振蕩剪切模式下磁流變脂法向力特性分析
    次級非中心位置工況雙邊直線電動機(jī)有限元分析
    旋轉(zhuǎn)穩(wěn)定二維修正彈鴨舵法向力計(jì)算模型研究
    磁極錯位削弱永磁直線伺服電動機(jī)齒槽法向力波動方法
    基于極弧系數(shù)選擇的直線電動機(jī)法向力波動削弱方法
    最新中文字幕久久久久| 免费人成在线观看视频色| 国产无遮挡羞羞视频在线观看| 亚洲精品久久午夜乱码| 久久久久久伊人网av| 免费观看a级毛片全部| 国产精品成人在线| 日本a在线网址| 美女午夜性视频免费| 国产视频一区二区在线看| 看黄色毛片网站| 欧美乱码精品一区二区三区| 老司机深夜福利视频在线观看| 热99国产精品久久久久久7| 水蜜桃什么品种好| 一级毛片精品| 欧美人与性动交α欧美精品济南到| 午夜精品在线福利| 日韩大码丰满熟妇| 视频在线观看一区二区三区| 国产精品久久久久久人妻精品电影| 午夜福利一区二区在线看| 99国产精品一区二区三区| 久久人人爽av亚洲精品天堂| 中文字幕最新亚洲高清| 1024视频免费在线观看| 亚洲av成人不卡在线观看播放网| 757午夜福利合集在线观看| 超碰成人久久| 精品福利观看| 变态另类成人亚洲欧美熟女 | av电影中文网址| 免费久久久久久久精品成人欧美视频| 高清毛片免费观看视频网站 | 精品卡一卡二卡四卡免费| 欧美日韩瑟瑟在线播放| 少妇 在线观看| 久久久久视频综合| 18禁裸乳无遮挡免费网站照片 | 国产欧美亚洲国产| 欧美日本中文国产一区发布| 大陆偷拍与自拍| 性色av乱码一区二区三区2| 美女高潮到喷水免费观看| 久久久精品免费免费高清| 丁香欧美五月| 亚洲五月色婷婷综合| 女同久久另类99精品国产91| 国产一区二区三区视频了| 欧美国产精品va在线观看不卡| 欧美乱色亚洲激情| 亚洲一区二区三区不卡视频| 国产乱人伦免费视频| 日日爽夜夜爽网站| 我的亚洲天堂| 一区二区三区国产精品乱码| 欧美黄色淫秽网站| av网站在线播放免费| av网站在线播放免费| 黄色片一级片一级黄色片| 久久久国产成人免费| 久久这里只有精品19| 桃红色精品国产亚洲av| 精品福利永久在线观看| 1024视频免费在线观看| 老司机影院毛片| 免费观看人在逋| 久久久久久久精品吃奶| 人人妻人人爽人人添夜夜欢视频| 国产精品二区激情视频| 亚洲九九香蕉| 精品人妻熟女毛片av久久网站| 亚洲熟妇熟女久久| 国产精品自产拍在线观看55亚洲 | 十分钟在线观看高清视频www| 极品人妻少妇av视频| 一进一出抽搐动态| 精品国产一区二区久久| 三级毛片av免费| 亚洲综合色网址| 夫妻午夜视频| 国精品久久久久久国模美| 真人做人爱边吃奶动态| 午夜日韩欧美国产| 中文字幕高清在线视频| 成人永久免费在线观看视频| 国产精品98久久久久久宅男小说| 大香蕉久久成人网| 久久久久国内视频| а√天堂www在线а√下载 | 69精品国产乱码久久久| 国产亚洲精品久久久久久毛片 | 又黄又粗又硬又大视频| 女人爽到高潮嗷嗷叫在线视频| e午夜精品久久久久久久| 日本撒尿小便嘘嘘汇集6| 亚洲综合色网址| 啦啦啦在线免费观看视频4| 欧美精品亚洲一区二区| 欧美日韩黄片免| 国产深夜福利视频在线观看| 久久国产精品男人的天堂亚洲| 亚洲精品一二三| 午夜福利在线免费观看网站| 桃红色精品国产亚洲av| 免费在线观看黄色视频的| 伦理电影免费视频| 亚洲在线自拍视频| 热re99久久国产66热| 免费在线观看视频国产中文字幕亚洲| 久久久久视频综合| 国产精华一区二区三区| 一进一出抽搐动态| 在线观看免费高清a一片| 亚洲自偷自拍图片 自拍| 久久香蕉精品热| 一a级毛片在线观看| 露出奶头的视频| 欧美激情极品国产一区二区三区| 精品一品国产午夜福利视频| 欧美精品高潮呻吟av久久| 黑人巨大精品欧美一区二区蜜桃| 欧美亚洲日本最大视频资源| 欧美不卡视频在线免费观看 | 在线永久观看黄色视频| 国产xxxxx性猛交| 99久久人妻综合| 国产99白浆流出| 精品国产亚洲在线| 一二三四在线观看免费中文在| 老司机在亚洲福利影院| 超碰成人久久| 757午夜福利合集在线观看| 欧美 亚洲 国产 日韩一| 亚洲国产欧美网| 国精品久久久久久国模美| xxxhd国产人妻xxx| 欧美大码av| 色尼玛亚洲综合影院| 久久草成人影院| 成人国语在线视频| 国产成人免费观看mmmm| 美女 人体艺术 gogo| 色婷婷av一区二区三区视频| 99久久国产精品久久久| 精品国产一区二区久久| 99精品久久久久人妻精品| 热99re8久久精品国产| av福利片在线| 在线观看日韩欧美| 高清毛片免费观看视频网站 | 亚洲欧美色中文字幕在线| 亚洲第一av免费看| 国产精品久久久av美女十八| 巨乳人妻的诱惑在线观看| 搡老岳熟女国产| 午夜福利乱码中文字幕| 丰满的人妻完整版| 亚洲精品美女久久久久99蜜臀| 亚洲熟妇中文字幕五十中出 | 91麻豆精品激情在线观看国产 | 日韩 欧美 亚洲 中文字幕| 国产真人三级小视频在线观看| 久久狼人影院| 在线观看免费视频网站a站| 免费黄频网站在线观看国产| 99国产极品粉嫩在线观看| 女人高潮潮喷娇喘18禁视频| 欧美+亚洲+日韩+国产| av线在线观看网站| 99国产极品粉嫩在线观看| 精品亚洲成国产av| videos熟女内射| 99久久人妻综合| 国产亚洲精品第一综合不卡| 日韩成人在线观看一区二区三区| 亚洲国产看品久久| 搡老岳熟女国产| 国产成人系列免费观看| 国产xxxxx性猛交| 久久亚洲精品不卡| 国产真人三级小视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 精品亚洲成国产av| 波多野结衣av一区二区av| 村上凉子中文字幕在线| 十八禁高潮呻吟视频| 黑人巨大精品欧美一区二区蜜桃| 欧美亚洲日本最大视频资源| 精品国产亚洲在线| 亚洲精品国产区一区二| 午夜两性在线视频| 国产视频一区二区在线看| 午夜两性在线视频| 亚洲精品中文字幕在线视频| 伊人久久大香线蕉亚洲五| 久久精品aⅴ一区二区三区四区| 久久久久久久午夜电影 | 国产成人系列免费观看| 日韩欧美一区视频在线观看| 天堂√8在线中文| 香蕉丝袜av| 欧美亚洲 丝袜 人妻 在线| 久久青草综合色| 人妻丰满熟妇av一区二区三区 | 侵犯人妻中文字幕一二三四区| 国产亚洲一区二区精品| 18禁观看日本| 欧美av亚洲av综合av国产av| 亚洲精品久久成人aⅴ小说| 欧美激情 高清一区二区三区| 人人妻人人爽人人添夜夜欢视频| 久久精品人人爽人人爽视色| 极品教师在线免费播放| 日韩人妻精品一区2区三区| 十八禁人妻一区二区| 国产熟女午夜一区二区三区| aaaaa片日本免费| 精品欧美一区二区三区在线| 最新美女视频免费是黄的| 亚洲欧美日韩高清在线视频| av国产精品久久久久影院| 人人妻人人添人人爽欧美一区卜| 免费在线观看日本一区| 成人精品一区二区免费| 国产精品国产高清国产av | 成在线人永久免费视频| 18禁美女被吸乳视频| 日韩人妻精品一区2区三区| 很黄的视频免费| 国产成人啪精品午夜网站| 精品第一国产精品| 丁香六月欧美| 不卡av一区二区三区| 身体一侧抽搐| 国产亚洲精品第一综合不卡| 国产精品国产av在线观看| www.999成人在线观看| 午夜福利在线免费观看网站| 51午夜福利影视在线观看| 男人的好看免费观看在线视频 | 好男人电影高清在线观看| 久久精品aⅴ一区二区三区四区| 一本一本久久a久久精品综合妖精| 色94色欧美一区二区| 亚洲av熟女| 欧美日韩一级在线毛片| 老司机午夜福利在线观看视频| 9色porny在线观看| 午夜福利影视在线免费观看| 一二三四在线观看免费中文在| 国产精品综合久久久久久久免费 | 日本黄色日本黄色录像| 国产精品自产拍在线观看55亚洲 | 精品久久久精品久久久| 久热爱精品视频在线9| 一区二区三区激情视频| 精品人妻1区二区| av视频免费观看在线观看| 叶爱在线成人免费视频播放| 嫁个100分男人电影在线观看| 中出人妻视频一区二区| 精品一区二区三区av网在线观看| 欧美人与性动交α欧美软件| 欧美成人午夜精品| 国产亚洲精品第一综合不卡| 欧美黄色片欧美黄色片| 99国产精品一区二区蜜桃av | 成年版毛片免费区| 大型黄色视频在线免费观看| 国产又色又爽无遮挡免费看| 啦啦啦 在线观看视频| 欧美精品人与动牲交sv欧美| 很黄的视频免费| 在线播放国产精品三级| 看黄色毛片网站| 亚洲精品成人av观看孕妇| 天天躁狠狠躁夜夜躁狠狠躁| 99精品欧美一区二区三区四区| 亚洲av日韩在线播放| 波多野结衣一区麻豆| 日本vs欧美在线观看视频| 日韩欧美在线二视频 | 精品人妻熟女毛片av久久网站| 亚洲欧美一区二区三区黑人| 男人操女人黄网站| 午夜久久久在线观看| 亚洲成人手机| 成人特级黄色片久久久久久久| 欧美av亚洲av综合av国产av| 成年人免费黄色播放视频| 亚洲精品自拍成人| 国产成人免费观看mmmm| 人人澡人人妻人| 精品福利观看| 亚洲美女黄片视频| 欧美激情极品国产一区二区三区| 在线观看午夜福利视频| www.熟女人妻精品国产| 91九色精品人成在线观看| 超碰97精品在线观看| 日韩欧美一区视频在线观看| 51午夜福利影视在线观看| 很黄的视频免费| 精品国内亚洲2022精品成人 | 午夜福利视频在线观看免费| 欧美黄色片欧美黄色片| 自线自在国产av| 一级a爱片免费观看的视频| 黄色视频,在线免费观看| 18禁黄网站禁片午夜丰满| 久久人人爽av亚洲精品天堂| 亚洲精品一卡2卡三卡4卡5卡| 欧美av亚洲av综合av国产av| 亚洲精品粉嫩美女一区| 亚洲国产欧美一区二区综合| 午夜久久久在线观看| 国产亚洲精品一区二区www | 男女下面插进去视频免费观看| 国产成人免费观看mmmm| 国产1区2区3区精品| 精品国产亚洲在线| 一进一出好大好爽视频| 日韩人妻精品一区2区三区| 高清av免费在线| 老司机深夜福利视频在线观看| 日韩制服丝袜自拍偷拍| 一二三四社区在线视频社区8| 这个男人来自地球电影免费观看| 亚洲精品国产一区二区精华液| 欧美乱色亚洲激情| 亚洲少妇的诱惑av| 99在线人妻在线中文字幕 | 一二三四在线观看免费中文在| 精品一区二区三卡| 性色av乱码一区二区三区2| 日本wwww免费看| 亚洲色图综合在线观看| 老汉色∧v一级毛片| 亚洲精品一二三| 亚洲精华国产精华精| 久久久国产一区二区| 中文字幕精品免费在线观看视频| 久久婷婷成人综合色麻豆| 亚洲av美国av| 久久精品亚洲av国产电影网| 久久久久久久久免费视频了| 亚洲国产欧美日韩在线播放| 亚洲精品一二三| 亚洲一区二区三区欧美精品| 国产精品永久免费网站| 超碰成人久久| 搡老岳熟女国产| 人成视频在线观看免费观看| 久久久国产成人精品二区 | 成人18禁高潮啪啪吃奶动态图| 又黄又爽又免费观看的视频| 日本欧美视频一区| 天天添夜夜摸| 亚洲免费av在线视频| 一本大道久久a久久精品| 91成人精品电影| 在线看a的网站| 免费av中文字幕在线| 国产精品免费视频内射| 一进一出好大好爽视频| 麻豆av在线久日| 国产片内射在线| 亚洲五月色婷婷综合| 日本精品一区二区三区蜜桃| 一级作爱视频免费观看| 热re99久久国产66热| 国产在线一区二区三区精| 久久中文字幕人妻熟女| 很黄的视频免费| 超碰成人久久| 国产不卡一卡二| 嫁个100分男人电影在线观看| 免费av中文字幕在线| 成年人午夜在线观看视频| 美女午夜性视频免费| 国产精品自产拍在线观看55亚洲 | 建设人人有责人人尽责人人享有的| 欧美黄色片欧美黄色片| 两性夫妻黄色片| 国产精品电影一区二区三区 | 又大又爽又粗| 男女高潮啪啪啪动态图| 日本黄色视频三级网站网址 | 国产精品久久电影中文字幕 | 精品国产一区二区三区四区第35| 高清黄色对白视频在线免费看| 久久久久国产一级毛片高清牌| 久久久久国产精品人妻aⅴ院 | 很黄的视频免费| 国产成+人综合+亚洲专区| 丁香欧美五月| 亚洲五月天丁香| 国产精品一区二区在线不卡| 纯流量卡能插随身wifi吗| 黄色视频不卡| 精品一区二区三区四区五区乱码| 免费一级毛片在线播放高清视频 | 久久香蕉激情| 免费一级毛片在线播放高清视频 | 69精品国产乱码久久久| 国产激情欧美一区二区| 亚洲av成人不卡在线观看播放网| 中文字幕人妻丝袜制服| 色综合欧美亚洲国产小说| 欧美日韩av久久| 两个人免费观看高清视频| 十八禁人妻一区二区| 欧美最黄视频在线播放免费 | 色综合婷婷激情| tube8黄色片| 国产精品 国内视频| 窝窝影院91人妻| 久久精品aⅴ一区二区三区四区| 色精品久久人妻99蜜桃| 法律面前人人平等表现在哪些方面| 一级黄色大片毛片| 女人久久www免费人成看片| 色婷婷av一区二区三区视频| 午夜福利影视在线免费观看| 国产成人系列免费观看| 最新的欧美精品一区二区| 窝窝影院91人妻| av天堂在线播放| 国产精品国产av在线观看| 母亲3免费完整高清在线观看| 亚洲精品国产色婷婷电影| 久久婷婷成人综合色麻豆| 精品国产一区二区三区久久久樱花| 国产一区二区三区在线臀色熟女 | 国产午夜精品久久久久久| 亚洲成a人片在线一区二区| 欧美日韩av久久| 天堂动漫精品| 国产单亲对白刺激| 欧美日韩av久久| 99国产精品免费福利视频| 中文字幕色久视频| 亚洲欧美日韩另类电影网站| 十八禁网站免费在线| 久久久国产成人精品二区 | 国产精品98久久久久久宅男小说| 99re6热这里在线精品视频| 色婷婷av一区二区三区视频| 在线观看免费视频日本深夜| 人人妻人人澡人人爽人人夜夜| 成年女人毛片免费观看观看9 | 午夜福利欧美成人| 最新的欧美精品一区二区| 精品国产国语对白av| 两个人免费观看高清视频| 久99久视频精品免费| 免费日韩欧美在线观看| 亚洲av成人av| 大片电影免费在线观看免费| www.自偷自拍.com| 亚洲国产欧美一区二区综合| 免费一级毛片在线播放高清视频 | 中出人妻视频一区二区| 欧美亚洲日本最大视频资源| 男女之事视频高清在线观看| 免费一级毛片在线播放高清视频 | 国产亚洲欧美精品永久| 国产99白浆流出| 满18在线观看网站| 亚洲色图综合在线观看| 日韩 欧美 亚洲 中文字幕| 中文字幕人妻熟女乱码| 一边摸一边抽搐一进一出视频| 欧美激情极品国产一区二区三区| 国产高清国产精品国产三级| 丝瓜视频免费看黄片| 色精品久久人妻99蜜桃| 精品国产国语对白av| 欧美黄色片欧美黄色片| 777米奇影视久久| 老熟妇仑乱视频hdxx| 在线看a的网站| 在线观看舔阴道视频| 国产精品免费一区二区三区在线 | 免费观看精品视频网站| 91成人精品电影| 一级作爱视频免费观看| 黄色成人免费大全| 国内毛片毛片毛片毛片毛片| 亚洲情色 制服丝袜| 一a级毛片在线观看| 黄色片一级片一级黄色片| 九色亚洲精品在线播放| 丝袜美腿诱惑在线| 淫妇啪啪啪对白视频| 日韩 欧美 亚洲 中文字幕| 美女 人体艺术 gogo| 婷婷成人精品国产| 精品福利永久在线观看| 好看av亚洲va欧美ⅴa在| 亚洲男人天堂网一区| 99久久综合精品五月天人人| 飞空精品影院首页| 国产精品久久视频播放| 无限看片的www在线观看| 老汉色av国产亚洲站长工具| 成熟少妇高潮喷水视频| 日韩成人在线观看一区二区三区| 如日韩欧美国产精品一区二区三区| 欧美国产精品va在线观看不卡| 成人亚洲精品一区在线观看| 正在播放国产对白刺激| 免费黄频网站在线观看国产| 一级作爱视频免费观看| 久热爱精品视频在线9| 天堂中文最新版在线下载| 亚洲中文日韩欧美视频| 天堂动漫精品| 黄色丝袜av网址大全| 啦啦啦免费观看视频1| 最新在线观看一区二区三区| 成人av一区二区三区在线看| a在线观看视频网站| 欧美精品亚洲一区二区| 欧美日韩乱码在线| 亚洲国产欧美网| 女人久久www免费人成看片| 国产精品一区二区免费欧美| 亚洲片人在线观看| 中文字幕色久视频| www.自偷自拍.com| 一级作爱视频免费观看| 亚洲精品中文字幕在线视频| cao死你这个sao货| 国产亚洲精品一区二区www | 在线播放国产精品三级| 在线观看一区二区三区激情| 热99re8久久精品国产| 不卡一级毛片| 18禁裸乳无遮挡免费网站照片 | 亚洲第一欧美日韩一区二区三区| 亚洲全国av大片| 国产极品粉嫩免费观看在线| 精品电影一区二区在线| 超色免费av| 亚洲欧美激情在线| 国产男女超爽视频在线观看| 9色porny在线观看| 欧美老熟妇乱子伦牲交| 99久久综合精品五月天人人| 国产深夜福利视频在线观看| 青草久久国产| 亚洲精品成人av观看孕妇| 久久99一区二区三区| 好看av亚洲va欧美ⅴa在| 丝袜美腿诱惑在线| 9191精品国产免费久久| 女人久久www免费人成看片| 国产成人精品无人区| 国产成人免费无遮挡视频| 十八禁网站免费在线| 黄片播放在线免费| 免费在线观看完整版高清| 日韩欧美在线二视频 | 亚洲人成伊人成综合网2020| 午夜视频精品福利| 成人黄色视频免费在线看| 亚洲精品中文字幕在线视频| 在线观看免费日韩欧美大片| 国产精品99久久99久久久不卡| 日本a在线网址| 欧美日韩瑟瑟在线播放| 性色av乱码一区二区三区2| 中文字幕精品免费在线观看视频| 捣出白浆h1v1| 亚洲综合色网址| 久久国产精品男人的天堂亚洲| 久久久国产成人精品二区 | ponron亚洲| 久久久国产成人免费| 在线观看www视频免费| 国产成人影院久久av| 久久国产精品影院| 色精品久久人妻99蜜桃| 国产精品1区2区在线观看. | 免费不卡黄色视频| 亚洲av成人不卡在线观看播放网| 欧美日韩一级在线毛片| ponron亚洲| 黄色片一级片一级黄色片| 久久人人97超碰香蕉20202| 日韩有码中文字幕| 好男人电影高清在线观看| 精品无人区乱码1区二区| 日韩制服丝袜自拍偷拍| 欧美黄色淫秽网站| 精品久久久久久,| 女人被狂操c到高潮| 涩涩av久久男人的天堂| 99久久精品国产亚洲精品| 精品亚洲成a人片在线观看| 99热网站在线观看| 亚洲成人免费av在线播放| 亚洲精品乱久久久久久| 天天影视国产精品| 日韩欧美在线二视频 | 中出人妻视频一区二区| 日韩 欧美 亚洲 中文字幕| 亚洲国产中文字幕在线视频| 久久久精品免费免费高清| 国产人伦9x9x在线观看| 午夜福利,免费看| 欧美最黄视频在线播放免费 | 久久久久久人人人人人|