陳啟東,左志全
(1.常熟理工學(xué)院,江蘇常熟 215500;2.常州大學(xué),江蘇常州 213016)
不同側(cè)向入口旋風(fēng)分離器流場(chǎng)數(shù)值分析
陳啟東1,左志全2
(1.常熟理工學(xué)院,江蘇常熟 215500;2.常州大學(xué),江蘇常州 213016)
利用雷諾應(yīng)力模型(RSM)對(duì)直切單入口、直切雙入口、斜切單入口、斜切雙入口、斜切螺旋面單入口、斜切螺旋面雙入口6種不同側(cè)向入口旋風(fēng)分離器內(nèi)部氣相流場(chǎng)進(jìn)行了計(jì)算分析。結(jié)果表明:雙入口結(jié)構(gòu)旋風(fēng)分離器內(nèi)部壓力場(chǎng)和速度場(chǎng)具有更好的對(duì)稱性與穩(wěn)定性;僅改變?nèi)肟谛鼻薪嵌葘?duì)旋風(fēng)分離器內(nèi)部速度場(chǎng)和壓力場(chǎng)的分布影響不大;當(dāng)本文中6種分離器內(nèi)部具有相近的切向速度徑向分布時(shí),斜切螺旋面入口結(jié)構(gòu)分離器壓力損失減少約25%,入口所需總壓降低17%,處理相同氣體量的能耗約下降17%;斜切螺旋面雙入口(XS-L型)分離器是一種綜合性能比較優(yōu)的旋風(fēng)分離器。
旋風(fēng)分離器;側(cè)向入口;壓降;入口總壓;能耗
旋風(fēng)分離器現(xiàn)在已經(jīng)廣泛應(yīng)用于化工、石油、礦山等工業(yè)部門,旋風(fēng)分離器的結(jié)構(gòu)也呈現(xiàn)出多樣性。結(jié)構(gòu)多樣性體現(xiàn)之一就是入口結(jié)構(gòu)多樣性,很多學(xué)者專家對(duì)不同入口結(jié)構(gòu)的旋風(fēng)分離器進(jìn)行了數(shù)值計(jì)算分析,王振興等[1]的研究成果表明,矩形入口長(zhǎng)寬比為1.5時(shí)壓降最低,對(duì)粒徑小于2.5 μm的顆粒,圓形入口分離效果最差,梯形入口分離效果最好;錢付平等[2]的研究闡述了不同入口截面角旋風(fēng)分離器的分離特性,最后得出入口具有一定截面角則有利于顆粒分離;吳彩金等[3]研究了斜切雙入口旋風(fēng)分離器流場(chǎng),發(fā)現(xiàn)本結(jié)構(gòu)能夠很好地改善單入口旋風(fēng)分離器流場(chǎng)的不對(duì)稱性,同時(shí)減小內(nèi)部局部渦流。
在此基礎(chǔ)上,本文對(duì)入口截面為矩形的直切單入口旋風(fēng)分離器(ZD型)、斜切單入口旋風(fēng)分離器(XD型)、直切雙入口旋風(fēng)分離器(ZS型)、斜切雙入口旋風(fēng)分離器(XS型)、斜切螺旋面單入口旋風(fēng)分離器(XD-L型)、斜切螺旋面雙入口旋風(fēng)分離器(XS-L型)等6種不同側(cè)向入口旋風(fēng)分離器在不同入口流速下的流場(chǎng)進(jìn)行計(jì)算,對(duì)比分析6種分離器內(nèi)壓力分布與速度分布,找出不同側(cè)向入口結(jié)構(gòu)對(duì)流場(chǎng)分布的影響。
根據(jù)以往學(xué)者研究結(jié)果[4~6]可知,由于雷諾應(yīng)力模型摒棄了各向同性的渦黏假定并且更加嚴(yán)格地考慮了流線彎曲、旋渦、旋轉(zhuǎn)和張力快速變化,能夠很好地預(yù)測(cè)分離器內(nèi)部的強(qiáng)制渦與自由渦,計(jì)算結(jié)果與實(shí)驗(yàn)測(cè)量值也最吻合,所以本文計(jì)算采用雷諾應(yīng)力模型(RSM)。RSM模型可以參考文獻(xiàn)[7]。
利用ICEM軟件采用六面體網(wǎng)格對(duì)本文中的模型進(jìn)行網(wǎng)格劃分。借用商業(yè)軟件Fluent對(duì)旋風(fēng)分離器內(nèi)部氣相流場(chǎng)進(jìn)行數(shù)值計(jì)算。求解器選擇基于壓力的求解器,采用非穩(wěn)態(tài)計(jì)算[8],時(shí)間步長(zhǎng)1E-4 s[9],計(jì)算方法采用SIMPLIC算法,壓力采用PRESTO離散格式,動(dòng)量采用二階離散格式,湍動(dòng)能、湍動(dòng)能耗散率采用QUICK離散格式,其他采用默認(rèn)格式。
本文計(jì)算采用的物理模型,雙入口旋風(fēng)分離器的入口高度是單入口旋風(fēng)分離器入口高度的一半,排氣管底端與入口底端對(duì)齊,不同單入口分離器入口高度和寬度分別相同,不同雙入口分離器入口高度和寬度分別相同,其他尺寸都相同。本文中ZD型模型采用參考文獻(xiàn)[10]中模型,結(jié)構(gòu)尺寸如下(單位:m):D=0.19,a=0.038,b=0.095,de=0.064,s=0.095,h1=0.285,h2=0.475,d0=0.072 5。旋風(fēng)分離器的物理模型如圖1所示。
圖1 旋風(fēng)分離器物理模型Fig.1 Geometry of cyclone
1)入口采用速度入口(velocity-inlet)。
2)排氣口采用壓力出口(pressure-outlet)[11],靜壓為0。
3)壁面采用無(wú)滑移邊界,近壁采用標(biāo)準(zhǔn)壁面函數(shù)處理。
為了驗(yàn)證數(shù)值計(jì)算模型和設(shè)置的合理性以及計(jì)算結(jié)果的可靠性,本文利用前述數(shù)值計(jì)算模型和設(shè)置計(jì)算文獻(xiàn)[10]中所給模型3個(gè)橫截面上的切向速度和軸向速度,并與文獻(xiàn)[10]中實(shí)驗(yàn)所測(cè)結(jié)果進(jìn)行比較,如圖2所示。比較結(jié)果顯示,運(yùn)用本文數(shù)值計(jì)算模型和設(shè)置計(jì)算出來(lái)的切向速度和軸向速度能夠較好地吻合文獻(xiàn)[10]中實(shí)驗(yàn)測(cè)得的數(shù)據(jù)。在切向速度徑向分布上捕捉到了切向速度隨徑向距離增加而增大的內(nèi)部強(qiáng)制渦以及切向速度隨徑向距離增加而減小的外部自由渦的Rankin組合渦的分布特點(diǎn),除此,在軸向速度徑向分布上還捕捉到了旋風(fēng)分離器內(nèi)壁面附近下旋流、中心區(qū)域上旋流的雙渦運(yùn)動(dòng)結(jié)構(gòu),并且內(nèi)外旋流的交接面位置與實(shí)驗(yàn)測(cè)得的結(jié)果基本一致。上述比較結(jié)果說(shuō)明,本文采用的數(shù)值計(jì)算模型及相關(guān)設(shè)置適用于后面旋風(fēng)分離器內(nèi)部流場(chǎng)計(jì)算。
對(duì)比了入口斜切角分別為5°、10°、15°、20°、25°的分離器在進(jìn)口為7.5 m/s的流速下的切向速度和軸向速度的徑向分布。隨著入口斜切角的增大,切向速度在強(qiáng)制渦區(qū)域內(nèi)徑向分布幾乎沒有變化,在自由渦區(qū)域,切向速度隨著入口角度增大而減??;軸向速度隨著入口斜切角度的增大而增大,在靠近壁面的下旋流中,軸向速度徑向分布相差不大,在中心區(qū)域的上旋流中,入口角度越大,軸向速度徑向分布相差越大??拷诿娴妮S向速度小一點(diǎn),有利于延長(zhǎng)顆粒的分離時(shí)間,靠近中心附近的上旋流軸向速度對(duì)顆粒分離的影響不明顯,如圖3所示。
圖2 旋風(fēng)分離器內(nèi)不同軸向位置上切向速度、軸向速度CFD計(jì)算結(jié)果與文獻(xiàn)[10]中LDA測(cè)量結(jié)果比較(從上到下橫截面為z=0.26 m,0.56 m,0.635 m)Fig.2 Comparison of the time averaged tangential and axial velocity between the RSM simulations and LDA measurements[10]at different axial sections(from top to bottom:sectionz=0.26 m,0.56 m,0.635 m)
圖3 不同入口斜切角XD型分離器3個(gè)橫截面上切向速度、軸向速度徑向分布(r為徑向位置,R為筒體半徑;從上到下橫截面為z=0.26 m,0.36 m,0.56 m)Fig.3 The radial profile for the time averaged tangential and axial velocity at different sections for different inlet angle cyclone XD(r is the radial distance,R is the cyclone radius;from top to bottom:sectionz=0.26 m,0.36 m,0.56 m)
圖4給出了分離器入口總壓及壓降隨著入口斜切角度增大的分布。從圖4可得出,隨著入口斜切角度的增大,入口總壓和壓降都在減小。與入口角度為5o的分離器相比,入口斜切角為10o的分離器內(nèi)切向速度分布曲線與之相差很小,兩者幾乎重合,入口總壓低0.7%,壓降低1%;入口角度為25°的分離器內(nèi)自由渦區(qū)域切向速度分布曲線比之要低0.5~1 m/s,占入口角度為5°的分離器內(nèi)切向速度峰值的6%左右,入口總壓低7.3%,壓降低8.8%。綜合對(duì)比得出,增大入口斜切角度有利于減阻,降低壓降,但效果有限,與此同時(shí)會(huì)減小自由渦區(qū)域切向速度,減弱分離效果。因此,本文所列的幾個(gè)入口傾斜角度中,入口傾斜角為10°的入口結(jié)構(gòu)相對(duì)較優(yōu)。
圖4 不同入口斜切角XD型分離器入口總壓、壓降Fig.4 The inlet total pressure and the pressure drop for different inlet angle cyclone XD
本節(jié)對(duì)ZD型、XD型(入口斜切角為10°)、ZS型、XS型(入口斜切角為10°)、XD-L型、XS-L型6種不同側(cè)向入口結(jié)構(gòu)分離器內(nèi)壓力場(chǎng)與速度場(chǎng)進(jìn)行分析討論。
4.3.1 壓力分布分析
圖5給出了6種分離器在入口流速為7.5 m/s下的兩個(gè)橫截面(z=0.26 m,0.56 m)上的靜壓徑向分布圖。觀察圖5中6種分離器靜壓徑向分布圖發(fā)現(xiàn),靜壓在徑向上由壁面往中心逐漸減少,由于中心處氣流高速旋轉(zhuǎn),中心區(qū)域形成了半徑很小的負(fù)壓區(qū)柱。從圖5還可看出,沿著徑向方向上壓力梯度最大。單入口結(jié)構(gòu)的分離器內(nèi)靜壓分布出現(xiàn)了擺動(dòng)現(xiàn)象,ZD型擺動(dòng)幅度最小,XD-L型擺動(dòng)幅度最大。靜壓分布的不均衡會(huì)導(dǎo)致分離器內(nèi)氣體流動(dòng)不均勻,出現(xiàn)氣體撞擊壁面現(xiàn)象,從而引發(fā)筒壁震動(dòng)。而雙入口結(jié)構(gòu)的分離器內(nèi)靜壓分布呈現(xiàn)出很好的對(duì)稱性。分離器內(nèi)部對(duì)稱的靜壓分布減少了因氣流擾動(dòng)而損耗的能量。對(duì)比后還可以得出,XD-L型、XS-L型分離內(nèi)部負(fù)壓氣柱要短于其他分離器,并且負(fù)壓區(qū)在排塵口上游區(qū)域就消失了,避免了負(fù)壓氣流將灰斗內(nèi)顆粒重新卷?yè)P(yáng)。在同樣的入口流速下,ZD型、XD型、ZS型、XS型4種分離器內(nèi)部靜壓徑向分布很接近,XD-L型、XS-L型分離器內(nèi)部靜壓徑向分布比較接近,并且在中心區(qū)域靜壓分布曲線要高于前四者,而在外部區(qū)域要低于前四者。
圖5 不同側(cè)向入口分離器不同截面上靜壓分布(從左到右為截面z=0.26 m,0.56 m)Fig.5 The contours plots for the time averaged flow variables of static pressure(from left to right:sectionz=0.26 m,0.56 m)
4.3.2 壓降分析
壓降又叫壓力損失,引起壓力損失的因素有很多[12],但主要是由于氣體與旋風(fēng)分離器內(nèi)壁摩擦所引起的能量損失和氣體因旋轉(zhuǎn)而產(chǎn)生的能量損耗。計(jì)算不同入口流速下流場(chǎng),根據(jù)計(jì)算結(jié)果繪制不同分離器各自壓降、入口總壓隨入口流速平方的分布圖以及各自壓降隨入口總壓分布圖,見圖6、圖7。
圖6 6種分離器壓降、入口總壓與入口流速平方關(guān)系Fig.6 The relationship between the pressure drop and the inlet velocity squared as well as the inlet total pressure and the inlet velocity squared
圖7 6種分離器壓降與入口總壓關(guān)系Fig.7 The relationship between the pressure drop and the inlet total pressure
由圖6a可知,在相同入口流速下,ZD型、ZS型、XD型、XS型4種分離器內(nèi)部壓降高于XD-L型和XS-L型分離器內(nèi)部壓降,并且隨著入口流速的增大,壓降的差值也變得越來(lái)越大;ZD型、ZS型、XD型、XS型4種分離器壓降隨速度平方分布曲線幾乎重合,XD-L型和XS-L型分離器內(nèi)部壓降隨速度平方分布曲線幾乎重合,這說(shuō)明僅僅改變進(jìn)口結(jié)構(gòu)的傾斜角度并不能有效地降低分離器內(nèi)壓降,而旋風(fēng)分離器頂部為螺旋面結(jié)構(gòu)能很好地起到減小壓降作用,當(dāng)入口切入形式確定時(shí),分離器是相差180°的雙入口還是單入口結(jié)構(gòu)對(duì)分離器壓降影響很小。觀察圖6b發(fā)現(xiàn),入口總壓與入口流速平方分布趨勢(shì)和壓降與入口流速平方的分布趨勢(shì)相同,在同樣的入口流速下頂部為螺旋面結(jié)構(gòu)的分離器入口總壓要小于其他4種分離器入口總壓,根據(jù)N=PQ可知,其所消耗的能量也越小。
表1 旋風(fēng)分離器性能參數(shù)表Table 1 Cyclone performance parameters
從圖7中可以看出,旋風(fēng)分離器壓降與入口總壓成線性關(guān)系,把壓降與入口總壓的比值定義為壓降率,用來(lái)衡量壓力損耗程度。4種旋風(fēng)分離器各自壓降率不隨入口流速的變化而變化,為一定值。各自分離器壓降率列于表1。從表1中可以得出,XD-L型、XS-L型兩種分離器壓降率要低于其他4種分離器,差值約為0.07。
4.3.3 速度分析
圖8為6種分離器在入口流速為7.5 m/s下縱截面(Y=0)上的切向、軸向和徑向速度分布云圖。
圖8 從上到下為6種分離器縱截面(Y=0)切向速度、軸向速度、徑向速度云圖Fig.8 The contours plots for the time averaged flow variables at Y=0 from top to bottom:the tangential velocity,axial velocity and radial velocity
從圖8切向速度云圖可以看出在軸向上切向速度徑向分布具有一致性,單入口結(jié)構(gòu)分離器內(nèi)切向速度分布具有擺動(dòng)性,并且在排塵口下方中心軸附近出現(xiàn)了切向速度數(shù)值為負(fù)的現(xiàn)象,說(shuō)明這里出現(xiàn)了二次渦,而二次渦容易引起顆粒返混,相反地,雙入口結(jié)構(gòu)分離器內(nèi)切向速度徑向分布呈現(xiàn)了很好的對(duì)稱性,并且沒有出現(xiàn)二次渦現(xiàn)象。從云圖中也可以直觀地看出切向速度最大值分布在與排氣管直徑相近的圓柱面上。從圖8軸向速度云圖可以看出壁面附近軸向速度為向下的下旋流,中心區(qū)域軸向速度為向上的上旋流,并且軸向速度最值出現(xiàn)在排氣口下端處。下旋流能夠?qū)⒈诿娓浇念w粒輸送到灰斗內(nèi),上旋流將氣體排出分離器。單入口分離器內(nèi)軸向速度分布也呈現(xiàn)擺動(dòng)現(xiàn)象,而雙入口分離器呈現(xiàn)出很好的對(duì)稱性。從圖8徑向速度云圖可以看出單入口分離器在中心軸兩邊出現(xiàn)了正負(fù)對(duì)應(yīng)的徑向速度分布,說(shuō)明氣體在這些地方出現(xiàn)了左右擺動(dòng)的流動(dòng)狀態(tài),而雙入口分離器呈現(xiàn)出較好的對(duì)稱性。
圖9 6種分離器不同橫截面上不同入口流速下切向速度徑向分布Fig.9 The radial profile for the time averaged tangential velocity at different sections for different inlet velocity
根據(jù)之前切向速度云圖分析中得到的旋風(fēng)分離器內(nèi)部切向速度徑向分布在軸向上具有一致性的特點(diǎn),圖9給出了分離器直筒部分一個(gè)橫截面z=0.56 m和錐形管部分一個(gè)橫截面z=0.26 m上在不同入口流速下的切向速度徑向分布,通過比較發(fā)現(xiàn),不同分離器在相同位置的橫截面上切向速度分布趨勢(shì)相近。表2和表3統(tǒng)計(jì)了兩個(gè)橫截面上不同入口流速下切向速度峰值與入口流速的比值以及峰值所在的位置,由表2和表3可以發(fā)現(xiàn),入口流速增加,各個(gè)分離器內(nèi)橫截面上峰值所在的位置幾乎沒有變化,直筒部分切向速度峰值所在圓柱面半徑小于錐形管部分,但差值很小,并且各個(gè)分離器間橫截面上峰值所在位置偏差較小。隨著入口流速的增加,切向速度峰值與入口流速比值變化很小,故可用兩個(gè)橫截面的均值來(lái)表示分離器內(nèi)流場(chǎng)切向速度與入口流速的比值,vZDτmax/v=1.73,vXDτmax/v=1.75 , vZSτmax/v=1.76 , vXSτmax/v=1.71 ,vXD-Lτmax/v=1.44 ,vXS-Lτmax/v=1.48 ,結(jié)合表1,可得到壓降、入口總壓與分離器內(nèi)最大切向速度最大值平方的比值,結(jié)果列于表4。從表4中可以發(fā)現(xiàn),ZD型、XD型、ZS型、XS型4種分離器壓降、入口總壓與切向速度峰值的比值很接近,XD-L型、XS-L型兩種分離器壓降、入口總壓與切向速度峰值的比值很接近。結(jié)合圖10和之前所論述,當(dāng)6種分離器具有相同的切向速度峰值的時(shí)候,切向速度徑向分布也很相近,并且后兩者入口流速約是前四者入口流速的1.2倍。結(jié)合表4可得如下結(jié)論,當(dāng)分離器內(nèi)部具有相近的切向速度徑向分布的時(shí)候,螺旋面入口結(jié)構(gòu)分離器壓力損失減少約25%,入口所需總壓要低17%左右,處理相同氣體量的能耗約減少17%。
表2 橫截面z=0.26 m上切向速度最大值與入口流速比值及最大切向速度位置Table 2 Comparison between the maximum tangential velocity value and ratio of inlet velocity and its position at sectionz=0.26 m
表3 橫截面z=0.56 m上切向速度最大值與入口流速比值及最大切向速度位置Table 3 Comparison between the maximum tangential velocity value and ratio of inlet velocity and its position at sectionz=0.56 m
表4 分離器壓降、入口總壓與切向速度峰值平方比值Table 4 The ratio between the pressure drop and the maximum tangential velocity squared as well as the inlet total pressure and the maximum tangential velocity squared
圖10 具有相同切向速度峰值的切向速度徑向分布Fig.10 The radial profile for the time averaged tangential velocity at different sections with same maximum tangential velocity
1)雙入口結(jié)構(gòu)旋風(fēng)分離器內(nèi)部壓力場(chǎng)和速度場(chǎng)具有更好的對(duì)稱性與穩(wěn)定性。
2)僅改變?nèi)肟谛鼻薪嵌葘?duì)旋風(fēng)分離器內(nèi)部速度場(chǎng)和壓力場(chǎng)的分布影響不大。
3)當(dāng)本文中6種分離器內(nèi)部具有相近的切向速度徑向分布的時(shí)候,螺旋面入口結(jié)構(gòu)分離器壓力損失要減少約25%,入口所需總壓要低17%左右,處理相同氣體量的能耗要下降17%左右。
4)XS-L型分離器是一種綜合性能較優(yōu)的旋風(fēng)分離器。
[1]王振興,畢榮山,李玉剛,等.入口形狀對(duì)旋風(fēng)分離器性能的影響[J].青島科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,33(3):277-281.
[2]錢付平,黃星瑋,紀(jì)國(guó)富.不同入口截面角旋風(fēng)分離器的分離特性[J].鍋爐技術(shù),2010,41(4):14-18.
[3]吳彩金,馬正飛,韓 虹.斜切雙入口旋風(fēng)分離器流場(chǎng)的數(shù)值模擬[J].煤礦機(jī)械,2010,31(9):61-63.
[4]馮喜平,趙勝海,李進(jìn)賢,等.不同湍流模型對(duì)旋渦流動(dòng)的數(shù)值模擬[J].航空動(dòng)力學(xué)報(bào),2011,26(6):1209-1214.
[5]張 雅,劉淑艷,王保國(guó).使用RSM模型對(duì)旋風(fēng)除塵器內(nèi)湍流各向異性的數(shù)值模擬[J].工程熱物理學(xué)報(bào),2005,26(增):41-44.
[6]Mehdi Azadia,Mohsen Azadib,Ali Mohebbia.A CFD study of the effect of cyclone size on its performance parameters[J].Journal of Hazardous Materials,2010,182:835-841.
[7]張師帥.計(jì)算流體動(dòng)力學(xué)及其應(yīng)用[M].武漢:華中科技大學(xué)出版社,2011.
[8]趙 通,楊亞平,劉俊龍.穩(wěn)態(tài)與非穩(wěn)態(tài)下旋風(fēng)分離器氣相流場(chǎng)數(shù)值模擬方法研究[J].動(dòng)力工程學(xué)報(bào),2012,32(8):591-597.
[9]Khairy Elsayed,Chris Lacor.Numerical modeling of the flow field and performance in cyclones of different cone-tip diameters[J].Computers&Fluids,2011,51:48-59.
[10]Fraser S M,Abdel-Razek A M,Abdullah M Z.Computational and experimental investigations in a cyclone dust separator[J].Journal of Process Mechanical Engineering,1997,211:247.
[11]王玉梅,楊淑霞,孫國(guó)剛.旋風(fēng)分離器出口邊界條件對(duì)其流場(chǎng)模擬的影響[J].化工裝備技術(shù),2008,29(2):1-4.
[12]金國(guó)淼.除塵設(shè)備[M].北京:化學(xué)工業(yè)出版社,2002.
Numerical analysis on different lateral inlet cyclone flow field
Chen Qidong1,Zuo Zhiquan2
(1.Changshu Institute of Technology,Changshu,Jiangsu 215500,China;2.Changzhou University,Changzhou,Jiangsu 213016,China)
Numerical simulation and calculation were performed on the gas-phase flow field of the straight tangential cyclone with single entry,straight tangential cyclone with double entries,bias tangential cyclone with single entry,bias tangential cyclone with double entries,bias tangential helicoidal cyclone with single entry,bias tangential helicoidal cyclone with double entries by using RSM model.Results showed that,the pressure field and velocity field of the cyclone with double entries had better symmetry and stability;only changing the inlet angle had a limited influence on the pressure field and velocity field.When the tangential velocity distribution were similar in the six kinds of cyclone studied in this article,compared with other cyclones’,the pressure drop of the bias tangential helicoidal cyclone decreased about 25%;the inlet total pressure needed was about 17%lower;the energy consumption dropped about 17%.Moreover,bias tangential helicoidal cyclone with double entries(XS-L type)had a better energy efficiency ratio,and its performance was the most excellent.
cyclone;lateral inlet;pressure drop;inlet total pressure;energy consumption
O355;TH48
A
1009-1742(2014)02-0058-10
2013-08-29
陳啟東,1962年出生,男,江蘇無(wú)錫市人,教授,長(zhǎng)期從事機(jī)械振動(dòng)、流體機(jī)械的教學(xué)和研究工作;E-mail:cqd666@163.com