• <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ī)法向力波動削弱方法
    婷婷丁香在线五月| 变态另类丝袜制服| 宅男免费午夜| 国产成人福利小说| 一个人免费在线观看的高清视频| 亚洲狠狠婷婷综合久久图片| 最近中文字幕高清免费大全6 | 少妇人妻精品综合一区二区 | 给我免费播放毛片高清在线观看| 嫁个100分男人电影在线观看| 夜夜爽天天搞| 国产高清有码在线观看视频| 在线免费观看的www视频| 床上黄色一级片| 成人特级av手机在线观看| www.999成人在线观看| 一二三四社区在线视频社区8| 国产精品免费一区二区三区在线| 婷婷色综合大香蕉| 久久国产乱子免费精品| 国产不卡一卡二| 国产免费av片在线观看野外av| 男人舔女人下体高潮全视频| 91麻豆精品激情在线观看国产| 搡女人真爽免费视频火全软件 | 禁无遮挡网站| 一个人免费在线观看电影| 有码 亚洲区| 久久久久久大精品| 他把我摸到了高潮在线观看| 成人国产一区最新在线观看| 中文字幕人成人乱码亚洲影| 丰满的人妻完整版| 国产精品一及| or卡值多少钱| 日本免费a在线| 国内毛片毛片毛片毛片毛片| 老鸭窝网址在线观看| av在线老鸭窝| 亚洲在线自拍视频| 搡老熟女国产l中国老女人| 精品国产亚洲在线| 美女大奶头视频| 久久中文看片网| 亚洲欧美日韩卡通动漫| 18禁黄网站禁片午夜丰满| 中亚洲国语对白在线视频| xxxwww97欧美| 高清在线国产一区| 看黄色毛片网站| 久久国产精品影院| 亚洲自拍偷在线| 精品乱码久久久久久99久播| 国产精品女同一区二区软件 | 久久精品夜夜夜夜夜久久蜜豆| 亚洲 欧美 日韩 在线 免费| 国产精品人妻久久久久久| 亚洲精品乱码久久久v下载方式| 精品无人区乱码1区二区| 亚洲精华国产精华精| 国产精品久久久久久人妻精品电影| 国产午夜精品久久久久久一区二区三区 | 欧洲精品卡2卡3卡4卡5卡区| 日本三级黄在线观看| 精品日产1卡2卡| 一级黄片播放器| 国产日本99.免费观看| 麻豆成人午夜福利视频| 日韩欧美精品免费久久 | 国产在线男女| 国产高清三级在线| 午夜激情欧美在线| 亚洲欧美激情综合另类| 最近在线观看免费完整版| 一二三四社区在线视频社区8| 黄片小视频在线播放| 久久99热这里只有精品18| 亚洲av.av天堂| 亚洲国产欧洲综合997久久,| 制服丝袜大香蕉在线| 欧美成狂野欧美在线观看| 国产视频内射| 最近最新中文字幕大全电影3| 午夜日韩欧美国产| 国产免费男女视频| 一区二区三区激情视频| 国产三级中文精品| 深夜a级毛片| 久久久久久久久大av| 一级毛片久久久久久久久女| 最后的刺客免费高清国语| 精品无人区乱码1区二区| 日本撒尿小便嘘嘘汇集6| 亚洲av电影在线进入| 欧美黑人巨大hd| 99热6这里只有精品| 一级av片app| 国产精品久久视频播放| 很黄的视频免费| 99热精品在线国产| 老熟妇乱子伦视频在线观看| 亚洲精品在线观看二区| 免费电影在线观看免费观看| 欧美性猛交╳xxx乱大交人| 男插女下体视频免费在线播放| 两人在一起打扑克的视频| 久久久国产成人精品二区| 如何舔出高潮| 国产av一区在线观看免费| 在线播放无遮挡| 国产白丝娇喘喷水9色精品| 日韩中文字幕欧美一区二区| 久久热精品热| 亚洲自拍偷在线| 变态另类丝袜制服| 午夜福利在线观看吧| 一级作爱视频免费观看| 成年女人毛片免费观看观看9| 日韩大尺度精品在线看网址| 特级一级黄色大片| 黄色日韩在线| 亚洲精品亚洲一区二区| 国产高清激情床上av| 一二三四社区在线视频社区8| 欧美在线一区亚洲| 亚洲精品影视一区二区三区av| 九九久久精品国产亚洲av麻豆| 国产精品日韩av在线免费观看| a在线观看视频网站| 男人狂女人下面高潮的视频| 亚洲欧美精品综合久久99| 亚洲av成人精品一区久久| 欧美日韩瑟瑟在线播放| 91午夜精品亚洲一区二区三区 | 成年女人看的毛片在线观看| 午夜福利在线在线| 亚洲精品色激情综合| 我要搜黄色片| 成人特级黄色片久久久久久久| 精品久久久久久,| 国产成人福利小说| 性色av乱码一区二区三区2| 最近最新免费中文字幕在线| 成人一区二区视频在线观看| 一a级毛片在线观看| 99riav亚洲国产免费| 精品国内亚洲2022精品成人| 国产一区二区在线观看日韩| 久99久视频精品免费| 成人永久免费在线观看视频| 亚洲美女视频黄频| 能在线免费观看的黄片| 免费无遮挡裸体视频| 日韩欧美三级三区| 欧美成人免费av一区二区三区| 一级黄色大片毛片| 露出奶头的视频| 99久国产av精品| 变态另类成人亚洲欧美熟女| 俄罗斯特黄特色一大片| 蜜桃久久精品国产亚洲av| 国产极品精品免费视频能看的| 久久久久久大精品| 国产 一区 欧美 日韩| 欧美+亚洲+日韩+国产| 欧美成人性av电影在线观看| 精品久久久久久久久久免费视频| 亚洲自拍偷在线| www.www免费av| 国产在线男女| 国产高清激情床上av| 久久久国产成人免费| 国产成+人综合+亚洲专区| 久久99热这里只有精品18| 可以在线观看的亚洲视频| 伊人久久精品亚洲午夜| 丰满人妻一区二区三区视频av| 国产成人av教育| 又黄又爽又免费观看的视频| 综合色av麻豆| 国产在线精品亚洲第一网站| av国产免费在线观看| 国产亚洲av嫩草精品影院| 成人欧美大片| 久久香蕉精品热| 日韩欧美 国产精品| 色播亚洲综合网| 婷婷丁香在线五月| 国产毛片a区久久久久| 小说图片视频综合网站| 又粗又爽又猛毛片免费看| 十八禁网站免费在线| 亚洲av免费高清在线观看| 看免费av毛片| 一级黄色大片毛片| 中文字幕人成人乱码亚洲影| 嫩草影院精品99| 无人区码免费观看不卡| 亚洲午夜理论影院| 精品福利观看| 久久久久亚洲av毛片大全| 久久久久久大精品| 国产亚洲欧美98| 日本一本二区三区精品| 成人毛片a级毛片在线播放| 免费大片18禁| 午夜激情欧美在线| 亚洲三级黄色毛片| 国语自产精品视频在线第100页| 亚洲真实伦在线观看| 亚洲18禁久久av| 国产精品久久久久久人妻精品电影| 日本黄色片子视频| 亚洲在线自拍视频| 日韩欧美精品v在线| 少妇的逼好多水| 亚洲成av人片免费观看| 亚洲欧美日韩高清在线视频| 又黄又爽又免费观看的视频| 精品一区二区免费观看| 最近视频中文字幕2019在线8| 99精品久久久久人妻精品| 又爽又黄a免费视频| 99久久99久久久精品蜜桃| 国产精品影院久久| 国产精品一区二区三区四区久久| 中国美女看黄片| 亚洲av成人精品一区久久| 日韩中文字幕欧美一区二区| 热99在线观看视频| 69人妻影院| 亚洲第一欧美日韩一区二区三区| 老司机深夜福利视频在线观看| 精品久久久久久,| 国产一区二区激情短视频| 午夜老司机福利剧场| 中文字幕人妻熟人妻熟丝袜美| 国产毛片a区久久久久| 精品人妻视频免费看| 亚洲av一区综合| 99久久精品国产亚洲精品| 韩国av一区二区三区四区| 丁香六月欧美| 久久国产乱子免费精品| 免费看美女性在线毛片视频| 99精品久久久久人妻精品| av专区在线播放| 免费看a级黄色片| 两个人视频免费观看高清| 国产aⅴ精品一区二区三区波| 亚洲一区二区三区色噜噜| 国产精品自产拍在线观看55亚洲| 国产日本99.免费观看| 亚洲国产日韩欧美精品在线观看| 搡老岳熟女国产| 久久亚洲真实| 国产大屁股一区二区在线视频| 一进一出抽搐gif免费好疼| 听说在线观看完整版免费高清| 国产精品嫩草影院av在线观看 | 久久人人爽人人爽人人片va | 搞女人的毛片| 一a级毛片在线观看| 老鸭窝网址在线观看| 国产色爽女视频免费观看| 亚洲成人久久性| 久久久久久大精品| 99热这里只有是精品在线观看 | 亚洲av成人av| 激情在线观看视频在线高清| 亚洲av成人不卡在线观看播放网| 免费高清视频大片| 亚洲欧美激情综合另类| 亚洲国产精品合色在线| 精品人妻视频免费看| 精品久久国产蜜桃| 美女高潮的动态| 国产精品不卡视频一区二区 | 午夜激情福利司机影院| 亚洲无线观看免费| 自拍偷自拍亚洲精品老妇| 最近中文字幕高清免费大全6 | 成人无遮挡网站| 日本在线视频免费播放| 女同久久另类99精品国产91| 亚洲18禁久久av| 午夜精品久久久久久毛片777| 韩国av一区二区三区四区| 搡老岳熟女国产| 在线观看午夜福利视频| 成人亚洲精品av一区二区| 久久精品国产自在天天线| 99久久无色码亚洲精品果冻| 久久久久九九精品影院| 91在线观看av| 亚洲精品成人久久久久久| 亚洲狠狠婷婷综合久久图片| 国产精品永久免费网站| 在线观看免费视频日本深夜| 99国产综合亚洲精品| 国产又黄又爽又无遮挡在线| 国产极品精品免费视频能看的| 久久精品国产亚洲av香蕉五月| 亚洲av第一区精品v没综合| 亚洲一区高清亚洲精品| 99久久99久久久精品蜜桃| 男女下面进入的视频免费午夜| 男人的好看免费观看在线视频| 国产精品野战在线观看| 精品久久久久久久久亚洲 | 欧美成人一区二区免费高清观看| 午夜精品一区二区三区免费看| av天堂在线播放| 欧美xxxx黑人xx丫x性爽| 看免费av毛片| 丁香欧美五月| 男插女下体视频免费在线播放| 18+在线观看网站| 热99在线观看视频| 别揉我奶头~嗯~啊~动态视频| 好看av亚洲va欧美ⅴa在| 丰满的人妻完整版| 热99re8久久精品国产| 99久久九九国产精品国产免费| 精品午夜福利在线看| 熟女人妻精品中文字幕| 成人高潮视频无遮挡免费网站| 亚洲熟妇熟女久久| 欧美高清性xxxxhd video| 亚洲专区中文字幕在线| 丰满乱子伦码专区| 99热这里只有是精品在线观看 | 午夜亚洲福利在线播放| 欧美成人免费av一区二区三区| 全区人妻精品视频| 变态另类丝袜制服| 亚洲av不卡在线观看| 日本黄大片高清| 久久久精品欧美日韩精品| 狂野欧美白嫩少妇大欣赏| 狂野欧美白嫩少妇大欣赏| 伊人久久精品亚洲午夜| 日本精品一区二区三区蜜桃| 国产精品亚洲av一区麻豆| 久久久久久九九精品二区国产| 在线观看美女被高潮喷水网站 | 人人妻人人看人人澡| 精品一区二区免费观看| 久久精品国产99精品国产亚洲性色| 日本免费一区二区三区高清不卡| 一区二区三区四区激情视频 | 色噜噜av男人的天堂激情| 国产成人av教育| 女同久久另类99精品国产91| 亚洲国产欧美人成| 女人被狂操c到高潮| 小说图片视频综合网站| 美女黄网站色视频| 美女 人体艺术 gogo| 成人国产综合亚洲| 黄片小视频在线播放| 99精品久久久久人妻精品| 欧美激情国产日韩精品一区| 深夜a级毛片| 日本撒尿小便嘘嘘汇集6| 精品人妻1区二区| x7x7x7水蜜桃| 好男人电影高清在线观看| 一个人观看的视频www高清免费观看| 日韩欧美在线乱码| 欧美性猛交╳xxx乱大交人| 看黄色毛片网站| 久久久久久国产a免费观看| 久久久精品大字幕| 美女大奶头视频| 能在线免费观看的黄片| 久久久久亚洲av毛片大全| a在线观看视频网站| 亚洲一区二区三区不卡视频| 免费看美女性在线毛片视频| 欧美日韩瑟瑟在线播放| 久久天躁狠狠躁夜夜2o2o| 亚洲专区国产一区二区| 高清毛片免费观看视频网站| 精品一区二区三区视频在线观看免费| 久久性视频一级片| ponron亚洲| 热99re8久久精品国产| 久久欧美精品欧美久久欧美| 最新在线观看一区二区三区| eeuss影院久久| 亚洲av电影不卡..在线观看| 九九久久精品国产亚洲av麻豆| 在线播放国产精品三级| 97超视频在线观看视频| 亚洲三级黄色毛片| or卡值多少钱| 搞女人的毛片| 日本 欧美在线| 欧美区成人在线视频| 免费观看的影片在线观看| 精品一区二区三区人妻视频| 18禁黄网站禁片午夜丰满| 国产精品久久久久久亚洲av鲁大| 很黄的视频免费| 亚洲美女黄片视频| 欧美中文日本在线观看视频| 日本 av在线| 色尼玛亚洲综合影院| 宅男免费午夜| 日日摸夜夜添夜夜添小说| 露出奶头的视频| 成人永久免费在线观看视频| 男女视频在线观看网站免费| 又紧又爽又黄一区二区| 麻豆国产av国片精品| 日本免费一区二区三区高清不卡| 桃色一区二区三区在线观看| avwww免费| 99久久成人亚洲精品观看| 久久中文看片网| 久久久久久久午夜电影| 天堂av国产一区二区熟女人妻| 有码 亚洲区| 免费av观看视频| 最近视频中文字幕2019在线8| 熟妇人妻久久中文字幕3abv| 最近最新免费中文字幕在线| 身体一侧抽搐| xxxwww97欧美| 五月伊人婷婷丁香| 中文在线观看免费www的网站| 亚洲18禁久久av| 波多野结衣巨乳人妻| 国产精品女同一区二区软件 | 欧美bdsm另类| 国产大屁股一区二区在线视频| 国产精品一区二区三区四区久久| 亚洲午夜理论影院| 日韩av在线大香蕉| 亚洲精品在线美女| 美女 人体艺术 gogo| 久久国产乱子伦精品免费另类| 校园春色视频在线观看| 亚洲av.av天堂| 国产探花在线观看一区二区| 偷拍熟女少妇极品色| 婷婷亚洲欧美| 欧美乱妇无乱码| 91午夜精品亚洲一区二区三区 | 丝袜美腿在线中文| 51午夜福利影视在线观看| 夜夜躁狠狠躁天天躁| 亚洲综合色惰| 校园春色视频在线观看| 国产精品影院久久| 久9热在线精品视频| 亚洲 欧美 日韩 在线 免费| 中亚洲国语对白在线视频| 成人高潮视频无遮挡免费网站| 国产精品自产拍在线观看55亚洲| 欧美潮喷喷水| 男女床上黄色一级片免费看| www.熟女人妻精品国产| 婷婷色综合大香蕉| 国产成人欧美在线观看| 嫩草影视91久久| 亚洲成人久久爱视频| 精品一区二区三区人妻视频| 男人的好看免费观看在线视频| 免费av不卡在线播放| 国产精品久久久久久人妻精品电影| 搞女人的毛片| 欧美性感艳星| 赤兔流量卡办理| 床上黄色一级片| 十八禁国产超污无遮挡网站| 国产欧美日韩精品亚洲av| 9191精品国产免费久久| 嫩草影院入口| 欧美三级亚洲精品| 桃色一区二区三区在线观看| 日韩欧美一区二区三区在线观看| 啦啦啦韩国在线观看视频| 欧美日韩乱码在线| 国产精品久久久久久久电影| 国内精品久久久久久久电影| avwww免费| 久久久精品大字幕| 亚洲av一区综合| 欧美日韩综合久久久久久 | 国产亚洲精品久久久com| 此物有八面人人有两片| 亚洲男人的天堂狠狠| 亚洲熟妇中文字幕五十中出| 亚洲,欧美,日韩| 精品久久久久久久末码| 久久热精品热| 99热这里只有精品一区| 亚洲内射少妇av| 欧美乱色亚洲激情| 免费观看的影片在线观看| 婷婷精品国产亚洲av在线| 蜜桃久久精品国产亚洲av| 天堂√8在线中文| 一个人看视频在线观看www免费| 欧美成人一区二区免费高清观看| 大型黄色视频在线免费观看| 国产精品久久电影中文字幕| 天堂√8在线中文| 99国产精品一区二区三区| 俺也久久电影网| 婷婷精品国产亚洲av在线| 久久伊人香网站| 国产精品三级大全| 国产精品一区二区三区四区免费观看 | 一a级毛片在线观看| 午夜精品一区二区三区免费看| 超碰av人人做人人爽久久| 国产亚洲精品久久久久久毛片| 免费在线观看亚洲国产| 香蕉av资源在线| 在线国产一区二区在线| 毛片一级片免费看久久久久 | 97超视频在线观看视频| 亚洲人与动物交配视频| 亚洲中文字幕日韩| 搡老妇女老女人老熟妇| 国产免费男女视频| 亚洲黑人精品在线| 在线天堂最新版资源| 国产成人aa在线观看| 国产午夜精品久久久久久一区二区三区 | 亚洲一区二区三区不卡视频| 十八禁网站免费在线| 一本精品99久久精品77| 一a级毛片在线观看| 国产69精品久久久久777片| 丁香欧美五月| 高清在线国产一区| 精品久久久久久久久av| 国产av麻豆久久久久久久| 舔av片在线| 两性午夜刺激爽爽歪歪视频在线观看| 一区二区三区四区激情视频 | 一a级毛片在线观看| 日韩欧美国产在线观看| 十八禁网站免费在线| 9191精品国产免费久久| 亚洲黑人精品在线| 我的老师免费观看完整版| 男插女下体视频免费在线播放| 夜夜看夜夜爽夜夜摸| 欧美色视频一区免费| 91午夜精品亚洲一区二区三区 | 桃红色精品国产亚洲av| 亚洲av二区三区四区| 精品久久久久久久末码| 全区人妻精品视频| 成人av在线播放网站| 国产免费一级a男人的天堂| 毛片一级片免费看久久久久 | 丁香六月欧美| 久久久久久久午夜电影| 亚洲国产日韩欧美精品在线观看| 搡老妇女老女人老熟妇| 国产白丝娇喘喷水9色精品| 性色avwww在线观看| 三级国产精品欧美在线观看| 51国产日韩欧美| 日韩成人在线观看一区二区三区| 欧美极品一区二区三区四区| 国产伦精品一区二区三区视频9| 老女人水多毛片| 精品久久久久久久久av| 在线观看舔阴道视频| 女人十人毛片免费观看3o分钟| 一个人看的www免费观看视频| 免费观看的影片在线观看| 久久久久久九九精品二区国产| 99热这里只有是精品在线观看 | 欧美极品一区二区三区四区| 香蕉av资源在线| 欧美激情国产日韩精品一区| 国产真实乱freesex| 日韩欧美一区二区三区在线观看| 色5月婷婷丁香| 好男人在线观看高清免费视频| 国产高清激情床上av| 亚洲国产日韩欧美精品在线观看| 男人和女人高潮做爰伦理| 日韩中字成人| 免费av观看视频| 成年人黄色毛片网站| 午夜福利欧美成人| 9191精品国产免费久久| 精品欧美国产一区二区三| 亚洲美女视频黄频| 国产伦在线观看视频一区| 最好的美女福利视频网| 国产在线精品亚洲第一网站| 欧美在线一区亚洲| 国产男靠女视频免费网站| 亚洲在线观看片| 夜夜爽天天搞| 国产午夜精品久久久久久一区二区三区 | 90打野战视频偷拍视频| 动漫黄色视频在线观看| 蜜桃久久精品国产亚洲av| 极品教师在线免费播放| 麻豆久久精品国产亚洲av| 亚洲美女搞黄在线观看 | 久久精品人妻少妇| 久久精品综合一区二区三区| 亚洲久久久久久中文字幕|