• <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专区在线播放| 联通29元200g的流量卡| 亚洲最大成人中文| 亚洲电影在线观看av| 九九爱精品视频在线观看| 99国产精品一区二区蜜桃av| 免费人成视频x8x8入口观看| 少妇猛男粗大的猛烈进出视频 | 国产探花极品一区二区| 亚洲av二区三区四区| 看片在线看免费视频| 插逼视频在线观看| 成人三级黄色视频| 国产成人91sexporn| 青青草视频在线视频观看| 午夜福利在线观看吧| 1024手机看黄色片| 国产日韩欧美在线精品| 国内精品一区二区在线观看| 国产精品无大码| 久久这里只有精品中国| 国产精品爽爽va在线观看网站| 国产老妇伦熟女老妇高清| 欧美3d第一页| 午夜视频国产福利| 九九在线视频观看精品| 久久久久久伊人网av| a级毛片免费高清观看在线播放| 激情 狠狠 欧美| 亚洲欧美中文字幕日韩二区| 免费黄网站久久成人精品| 久久久欧美国产精品| 免费看光身美女| 欧美日韩一区二区视频在线观看视频在线 | 青春草亚洲视频在线观看| 欧美zozozo另类| 99久国产av精品国产电影| 亚洲一区高清亚洲精品| 91狼人影院| 精品不卡国产一区二区三区| 国产三级中文精品| 神马国产精品三级电影在线观看| 亚洲精品粉嫩美女一区| 日韩欧美精品免费久久| 日韩一区二区三区影片| 欧美+日韩+精品| 欧美日韩综合久久久久久| 国产免费一级a男人的天堂| 亚洲一级一片aⅴ在线观看| 中文字幕熟女人妻在线| 亚洲精品影视一区二区三区av| 99久国产av精品国产电影| 久久久久久久亚洲中文字幕| 国产伦理片在线播放av一区 | 国产一区二区在线av高清观看| 免费看美女性在线毛片视频| 99在线视频只有这里精品首页| 成人二区视频| 日韩一区二区视频免费看| 伦理电影大哥的女人| ponron亚洲| 色噜噜av男人的天堂激情| 狠狠狠狠99中文字幕| 简卡轻食公司| 色尼玛亚洲综合影院| 大又大粗又爽又黄少妇毛片口| 久久久a久久爽久久v久久| 亚洲人成网站高清观看| 亚洲av二区三区四区| 91av网一区二区| 亚洲av中文av极速乱| 日韩欧美精品v在线| 亚洲激情五月婷婷啪啪| 黑人高潮一二区| 一级毛片我不卡| 一级av片app| 99热精品在线国产| 三级经典国产精品| 亚州av有码| 听说在线观看完整版免费高清| 性插视频无遮挡在线免费观看| 精品久久久久久成人av| 国产伦精品一区二区三区四那| 色综合色国产| 亚洲av二区三区四区| 亚洲欧美精品自产自拍| 菩萨蛮人人尽说江南好唐韦庄 | 哪个播放器可以免费观看大片| 国产精品av视频在线免费观看| 亚洲最大成人av| 性色avwww在线观看| 好男人在线观看高清免费视频| 天堂网av新在线| 九草在线视频观看| 亚洲精品日韩av片在线观看| 亚洲七黄色美女视频| 国产美女午夜福利| 国产探花极品一区二区| 国产精品久久久久久精品电影| 欧美成人一区二区免费高清观看| 久久久精品大字幕| 黄色日韩在线| 国产免费一级a男人的天堂| 男女边吃奶边做爰视频| 亚洲婷婷狠狠爱综合网| 非洲黑人性xxxx精品又粗又长| 99热精品在线国产| 熟女电影av网| 草草在线视频免费看| 亚洲精品国产成人久久av| 91久久精品国产一区二区成人| 亚洲一区高清亚洲精品| 亚洲av中文字字幕乱码综合| 精品人妻视频免费看| 少妇的逼好多水| 亚洲婷婷狠狠爱综合网| 欧美最黄视频在线播放免费| 在线天堂最新版资源| 国产精品福利在线免费观看| 国产一区亚洲一区在线观看| 在线观看午夜福利视频| 亚洲美女搞黄在线观看| 欧美日韩综合久久久久久| 久久精品国产自在天天线| 中文字幕熟女人妻在线| 国产成人影院久久av| 国产又黄又爽又无遮挡在线| 变态另类成人亚洲欧美熟女| 黑人高潮一二区| 日韩欧美精品免费久久| 又爽又黄a免费视频| 日本黄色视频三级网站网址| 热99在线观看视频| 亚洲av第一区精品v没综合| 91精品国产九色| 在线播放无遮挡| 亚洲一区二区三区色噜噜| 欧美一区二区国产精品久久精品| 在现免费观看毛片| 又黄又爽又刺激的免费视频.| 男的添女的下面高潮视频| 97超视频在线观看视频| 两个人视频免费观看高清| 少妇人妻精品综合一区二区 | 男女下面进入的视频免费午夜| 久久久久久久午夜电影| 18禁在线无遮挡免费观看视频| 中文字幕人妻熟人妻熟丝袜美| 不卡视频在线观看欧美| 亚洲精品乱码久久久久久按摩| 久久久久久大精品| 亚洲人与动物交配视频| a级一级毛片免费在线观看| 你懂的网址亚洲精品在线观看 | 又粗又硬又长又爽又黄的视频 | 精品99又大又爽又粗少妇毛片| 又爽又黄a免费视频| 一个人看视频在线观看www免费| 不卡一级毛片| 国产成人福利小说| 一本久久中文字幕| 亚洲人成网站在线观看播放| 日本免费一区二区三区高清不卡| 偷拍熟女少妇极品色| 51国产日韩欧美| 99久久人妻综合| 国产免费男女视频| 少妇被粗大猛烈的视频| av卡一久久| 少妇裸体淫交视频免费看高清| 国产真实乱freesex| 在线免费观看不下载黄p国产| 欧美+日韩+精品| 精品熟女少妇av免费看| 乱人视频在线观看| 夜夜夜夜夜久久久久| 人人妻人人澡人人爽人人夜夜 | 爱豆传媒免费全集在线观看| 免费观看的影片在线观看| 天堂av国产一区二区熟女人妻| 久久精品91蜜桃| 国产精品日韩av在线免费观看| 亚洲av免费高清在线观看| 亚洲人成网站在线观看播放| 国产私拍福利视频在线观看| 夜夜夜夜夜久久久久| av福利片在线观看| 国产精品女同一区二区软件| 内地一区二区视频在线| 三级男女做爰猛烈吃奶摸视频| 国产男人的电影天堂91| 成人性生交大片免费视频hd| 尾随美女入室| 成人国产麻豆网| 我要搜黄色片| 亚洲欧美日韩高清在线视频| 永久网站在线| 偷拍熟女少妇极品色| 又爽又黄无遮挡网站| 国产伦一二天堂av在线观看| 91麻豆精品激情在线观看国产| 成人无遮挡网站| 亚洲图色成人| 成人毛片a级毛片在线播放| 五月玫瑰六月丁香| 国产成人aa在线观看| 老熟妇乱子伦视频在线观看| 亚洲一级一片aⅴ在线观看| 国产高清有码在线观看视频| 久久人妻av系列| 99视频精品全部免费 在线| 哪个播放器可以免费观看大片| 久久人人爽人人爽人人片va| 狂野欧美白嫩少妇大欣赏| 亚洲久久久久久中文字幕| 欧美日韩综合久久久久久| 国产私拍福利视频在线观看| 成人二区视频| 99视频精品全部免费 在线| 中国国产av一级| 非洲黑人性xxxx精品又粗又长| 日韩欧美在线乱码| 国产爱豆传媒在线观看| 青青草视频在线视频观看| 久久精品国产99精品国产亚洲性色| 日本成人三级电影网站| 一级毛片电影观看 | 精品久久久久久久人妻蜜臀av| 国产亚洲5aaaaa淫片| 亚洲国产欧美人成| 亚洲人成网站高清观看| 神马国产精品三级电影在线观看| 狠狠狠狠99中文字幕| 欧美最黄视频在线播放免费| 国产精品福利在线免费观看| 国产亚洲av嫩草精品影院| 黄色欧美视频在线观看| 永久网站在线| 最近视频中文字幕2019在线8| 亚洲av不卡在线观看| 在线观看免费视频日本深夜| 国产精品一及| 亚洲av不卡在线观看| 国产黄片美女视频| 久久久久久久午夜电影| 国产v大片淫在线免费观看| 尾随美女入室| 久久精品影院6| 韩国av在线不卡| 波野结衣二区三区在线| 午夜a级毛片| 亚洲欧美清纯卡通| 少妇的逼水好多| 免费看a级黄色片| 日本免费a在线| 欧美色欧美亚洲另类二区| 欧美+日韩+精品| 一边亲一边摸免费视频| 乱系列少妇在线播放| 国产黄片美女视频| 亚洲精品粉嫩美女一区| 麻豆av噜噜一区二区三区| 欧美人与善性xxx| 日韩中字成人| 狂野欧美白嫩少妇大欣赏| 久久久久免费精品人妻一区二区| 日韩中字成人| 久久韩国三级中文字幕| 日韩大尺度精品在线看网址| 岛国毛片在线播放| 免费看美女性在线毛片视频| 高清在线视频一区二区三区 | 欧美一级a爱片免费观看看| 欧美又色又爽又黄视频| 亚洲aⅴ乱码一区二区在线播放| 国产精品一二三区在线看| 国产精品99久久久久久久久| 国产精品免费一区二区三区在线| 久久久久久久久久久丰满| 亚洲欧美精品自产自拍| 国模一区二区三区四区视频| 99久久精品一区二区三区| 我的女老师完整版在线观看| 少妇人妻精品综合一区二区 | 国产激情偷乱视频一区二区| 深夜精品福利| 久99久视频精品免费| 欧美最新免费一区二区三区| 国产私拍福利视频在线观看| 在线播放国产精品三级| 日本黄色视频三级网站网址| 一级毛片我不卡| 老师上课跳d突然被开到最大视频| 天堂中文最新版在线下载 | 91在线精品国自产拍蜜月| 高清在线视频一区二区三区 | 国产三级中文精品| 日韩av在线大香蕉| 两性午夜刺激爽爽歪歪视频在线观看| 日韩制服骚丝袜av| 亚洲图色成人| 日韩一本色道免费dvd| 国产精品久久视频播放| 亚洲丝袜综合中文字幕| 中文精品一卡2卡3卡4更新| 欧美色欧美亚洲另类二区| 91久久精品电影网| av在线播放精品| 国产精品不卡视频一区二区| 国产高清激情床上av| 少妇猛男粗大的猛烈进出视频 | 精品久久久久久久久亚洲| videossex国产| 国产乱人偷精品视频| 国内少妇人妻偷人精品xxx网站| 18禁在线播放成人免费| 欧美高清成人免费视频www| 日本黄大片高清| 18禁裸乳无遮挡免费网站照片| 亚洲av成人精品一区久久| 精品一区二区免费观看| 少妇熟女欧美另类| 亚洲精品亚洲一区二区| 免费搜索国产男女视频| h日本视频在线播放| 久久久久久伊人网av| 亚洲久久久久久中文字幕| 身体一侧抽搐| 成人漫画全彩无遮挡| 亚洲欧美日韩高清在线视频| 黄片wwwwww| 卡戴珊不雅视频在线播放| 久久久a久久爽久久v久久| 蜜臀久久99精品久久宅男| 亚洲国产精品成人综合色| 国产精品麻豆人妻色哟哟久久 | 麻豆国产av国片精品| 人体艺术视频欧美日本| 青青草视频在线视频观看| 国产大屁股一区二区在线视频| 婷婷六月久久综合丁香| 日本撒尿小便嘘嘘汇集6| 欧美最新免费一区二区三区| 男女视频在线观看网站免费| 免费av观看视频| 国产成人福利小说| 直男gayav资源| 中文字幕av成人在线电影| 精品欧美国产一区二区三| av在线亚洲专区| 色综合亚洲欧美另类图片| 亚洲国产精品久久男人天堂| 亚洲真实伦在线观看| 女人十人毛片免费观看3o分钟| 色5月婷婷丁香| 97超视频在线观看视频| 精品日产1卡2卡| 精品人妻视频免费看| 国产老妇伦熟女老妇高清| 国产成人一区二区在线| 久久久精品欧美日韩精品| 亚洲国产精品久久男人天堂| 久久精品人妻少妇| 久久午夜亚洲精品久久| 国内精品宾馆在线| 18+在线观看网站| 啦啦啦观看免费观看视频高清| 欧美激情国产日韩精品一区| 日本与韩国留学比较| 久久6这里有精品| 国产女主播在线喷水免费视频网站 | 免费看光身美女| 美女 人体艺术 gogo| 国产一区二区三区av在线 | 2021天堂中文幕一二区在线观| 欧美区成人在线视频| 日本成人三级电影网站| 国产视频首页在线观看| 午夜福利成人在线免费观看| 男人舔奶头视频| 丝袜美腿在线中文| 午夜视频国产福利| 色哟哟·www| 嘟嘟电影网在线观看| 国产亚洲5aaaaa淫片| 男人和女人高潮做爰伦理| 国语自产精品视频在线第100页| 亚洲七黄色美女视频| 18+在线观看网站| 国产黄a三级三级三级人| 亚洲图色成人| 国产不卡一卡二| 国产中年淑女户外野战色| 国产亚洲精品久久久久久毛片| 欧美性猛交╳xxx乱大交人| 国产伦精品一区二区三区四那| 我的女老师完整版在线观看| 又爽又黄a免费视频| 热99在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 一区福利在线观看| 校园春色视频在线观看| 99久久精品国产国产毛片| 搞女人的毛片| av视频在线观看入口| 国产成人freesex在线| 欧美日韩乱码在线| 精品少妇黑人巨大在线播放 | 亚洲最大成人手机在线| 天堂中文最新版在线下载 | 亚洲性久久影院| 淫秽高清视频在线观看| 看片在线看免费视频| 床上黄色一级片| 免费黄网站久久成人精品| 日本黄大片高清| 亚洲精品国产av成人精品| 亚洲中文字幕日韩| 99久久精品国产国产毛片| 亚洲欧美精品自产自拍| 好男人视频免费观看在线| 一个人看的www免费观看视频| 一级毛片久久久久久久久女| 18+在线观看网站| 黄色视频,在线免费观看| 成人欧美大片| АⅤ资源中文在线天堂| 久久精品国产鲁丝片午夜精品| 婷婷色综合大香蕉| 成人二区视频| 亚洲欧美精品专区久久| 国产单亲对白刺激| 国产精品无大码| 嘟嘟电影网在线观看| 久久国产乱子免费精品| 婷婷色av中文字幕| 亚洲美女视频黄频| 久久精品国产亚洲av天美| 精品久久久久久久久亚洲| 精品久久久久久成人av| 1000部很黄的大片| 亚洲欧美日韩高清专用| 一进一出抽搐gif免费好疼| 日本一二三区视频观看| 99久久久亚洲精品蜜臀av| 久久精品综合一区二区三区| 热99在线观看视频| 一本一本综合久久| 中国美白少妇内射xxxbb| 精华霜和精华液先用哪个| 一区二区三区免费毛片| 国产成人精品一,二区 | 久久精品国产99精品国产亚洲性色| 国产精品一区二区性色av| 18禁裸乳无遮挡免费网站照片| 亚洲成人av在线免费| 国产精品女同一区二区软件| 国产免费男女视频| 日韩成人av中文字幕在线观看| 国产午夜精品论理片| 日韩av在线大香蕉| 欧美最新免费一区二区三区| 夜夜爽天天搞| 永久网站在线| 午夜福利在线观看吧| 免费av观看视频| 亚洲综合色惰| 中文字幕制服av| 22中文网久久字幕| 99riav亚洲国产免费| 黄色一级大片看看| 欧美一区二区国产精品久久精品| 尤物成人国产欧美一区二区三区| 免费看a级黄色片| 男人的好看免费观看在线视频| 欧美高清成人免费视频www| 男女下面进入的视频免费午夜| 久久99热6这里只有精品| 欧美日韩一区二区视频在线观看视频在线 | 一级黄片播放器| 亚洲成人久久性| 夜夜爽天天搞| 搞女人的毛片| 深爱激情五月婷婷| 啦啦啦韩国在线观看视频| av国产免费在线观看| 国产69精品久久久久777片| 色哟哟哟哟哟哟| 日本免费a在线| 免费av观看视频| 黑人高潮一二区| 最近的中文字幕免费完整| 成人毛片60女人毛片免费| 婷婷六月久久综合丁香| 国产成人精品婷婷| 日产精品乱码卡一卡2卡三| 有码 亚洲区| 午夜免费激情av| 国产亚洲欧美98| 看片在线看免费视频| 神马国产精品三级电影在线观看| 亚洲内射少妇av| 国产亚洲欧美98| 美女脱内裤让男人舔精品视频 | 国产一区二区在线av高清观看| 如何舔出高潮| 天天躁夜夜躁狠狠久久av| 欧美成人a在线观看| 日韩一区二区视频免费看| 2021天堂中文幕一二区在线观| 插阴视频在线观看视频| 看十八女毛片水多多多| a级毛片免费高清观看在线播放| 久久精品国产亚洲av天美| 2021天堂中文幕一二区在线观| 一本久久精品| 在线观看免费视频日本深夜| 97热精品久久久久久| 国内少妇人妻偷人精品xxx网站| 最近的中文字幕免费完整| 亚洲欧美精品专区久久| 久久人妻av系列| 99在线人妻在线中文字幕| 国产精品人妻久久久久久| 丝袜美腿在线中文| 夜夜看夜夜爽夜夜摸| 老司机影院成人| 2022亚洲国产成人精品| 久久久国产成人精品二区| 五月伊人婷婷丁香| 欧美极品一区二区三区四区| 亚洲自偷自拍三级| 赤兔流量卡办理| 亚洲最大成人av| 欧美最新免费一区二区三区| 好男人在线观看高清免费视频| 天天躁夜夜躁狠狠久久av| 精品99又大又爽又粗少妇毛片| 美女国产视频在线观看| 青春草视频在线免费观看| 搡老妇女老女人老熟妇| 三级经典国产精品| 日本一二三区视频观看| 亚洲欧美日韩无卡精品| 蜜桃久久精品国产亚洲av| 国产精品一及| 欧美日韩国产亚洲二区| 亚洲国产精品成人综合色| 亚洲国产高清在线一区二区三| 免费在线观看成人毛片| 精华霜和精华液先用哪个| 级片在线观看| 青春草亚洲视频在线观看| www.av在线官网国产| 欧美成人一区二区免费高清观看| 亚洲经典国产精华液单| 给我免费播放毛片高清在线观看| 精品人妻偷拍中文字幕| 少妇人妻精品综合一区二区 | 男人的好看免费观看在线视频| 成人毛片60女人毛片免费| 成人性生交大片免费视频hd| 欧美bdsm另类| 成年免费大片在线观看| 日韩一本色道免费dvd| 丰满的人妻完整版| 夜夜爽天天搞| 久久久久久久久久成人| 国产一区二区激情短视频| 久久久久性生活片| 亚洲中文字幕一区二区三区有码在线看| 嫩草影院新地址| a级一级毛片免费在线观看| 国产精品免费一区二区三区在线| 高清毛片免费看| 亚洲精品成人久久久久久| 午夜激情欧美在线| 久久精品91蜜桃| 能在线免费看毛片的网站| 1000部很黄的大片| 久久国内精品自在自线图片| 老熟妇乱子伦视频在线观看| 校园人妻丝袜中文字幕| 亚洲成av人片在线播放无| 国产一区二区在线av高清观看| 国产亚洲av片在线观看秒播厂 | 自拍偷自拍亚洲精品老妇| 久久久国产成人免费| 国产极品精品免费视频能看的| 免费在线观看成人毛片| 亚洲精品色激情综合| 亚洲美女搞黄在线观看| 2021天堂中文幕一二区在线观| 久99久视频精品免费| 久久这里有精品视频免费| 亚洲av二区三区四区| 91久久精品国产一区二区三区| 国产 一区 欧美 日韩| 国产三级在线视频| 男女下面进入的视频免费午夜| 欧美成人精品欧美一级黄| 22中文网久久字幕| 韩国av在线不卡| 最新中文字幕久久久久| 黑人高潮一二区| 长腿黑丝高跟| 欧美一区二区亚洲| 国产伦一二天堂av在线观看| 深爱激情五月婷婷| 亚洲成人av在线免费| 国产老妇伦熟女老妇高清| 久久99精品国语久久久| 欧美成人免费av一区二区三区| av免费在线看不卡| 美女被艹到高潮喷水动态|