基于MATLAB的攪拌機軸的強度計算

高 路，高勝寒

(1.吉林化工學(xué)院 機電工程學(xué)院，吉林 吉林 132022；2.中南大學(xué) 數(shù)學(xué)與統(tǒng)計學(xué)院,湖南 長沙 410083)

攪拌機是攪拌反應(yīng)器的重要組成部分.攪拌機中攪拌軸的主要作用是將扭矩從傳動裝置傳遞給攪拌器的葉片.攪拌軸的設(shè)計直接關(guān)系到攪拌反應(yīng)器的振動、軸封性能等，因此攪拌軸的設(shè)計尤其重要.作為攪拌裝置的主要零件之一的攪拌軸由電機驅(qū)動，通過減速器減速以一定的轉(zhuǎn)速實現(xiàn)攪拌功能[1].攪拌軸既與攪拌器連接,又要與軸封裝置以及軸承、聯(lián)軸器等組成軸系,且須保證以一定轉(zhuǎn)速回轉(zhuǎn),其設(shè)計計算需考慮諸多因素.

1 攪拌軸的力學(xué)模型

作用在攪拌軸上的力相當復(fù)雜，對于一個從容器頂部垂直懸掛的攪拌軸來說，主要有攪拌器及軸上其它零件的重力、慣性力、流體作用力、軸的轉(zhuǎn)矩、流體的徑向力、攪拌器和攪拌軸在組合重心處質(zhì)量偏心引起的離心力等[2].

攪拌軸通過剛性聯(lián)軸器同減速機輸出軸聯(lián)接，輸出軸用30215滾動軸承支撐，攪拌軸材料為0Cr18Ni9，分別用普通平鍵與兩個攪拌漿和剛性聯(lián)軸器聯(lián)接.本攪拌軸為懸臂軸，其力學(xué)模型如圖1所示[3].

圖1 力學(xué)模型

2 攪拌軸的扭轉(zhuǎn)變形計算

2.1 按扭轉(zhuǎn)變形計算

由于攪拌軸受轉(zhuǎn)矩和彎矩的聯(lián)合作用，保證攪拌軸單位長度的最大扭轉(zhuǎn)角γ在許用范圍之內(nèi)，避免由于扭轉(zhuǎn)變形過大造成攪拌軸振動，使軸封失效[4,5].

已知電機功率Pn=7.5 KW，攪拌轉(zhuǎn)數(shù)n=52 rpm，攪拌軸的直徑d=φ75 mm，如圖2所示[3].

圖2 軸零件圖

(1)

=1.35×106N·mm

式(1)中Mmax為攪拌軸傳遞的最大扭矩，N·m；η1為傳動側(cè)軸承前的傳動裝置效率，取η1=98%；Pn為電機功率，KW；n為攪拌軸的轉(zhuǎn)速，rpm；攪拌軸產(chǎn)生扭轉(zhuǎn)變形時，扭轉(zhuǎn)角γ.

(2)

=0.125°/m[γ]

式中[γ]：軸的許用轉(zhuǎn)角，°/m，對于懸臂軸[γ]=0.35°/m；G：軸材料剪切彈性模量，MPa，G=199×103MPa；N0：空心軸內(nèi)徑與外徑比，N0=0.經(jīng)計算扭轉(zhuǎn)角合格.

2.2 按扭轉(zhuǎn)變形的MATLAB計算

M函數(shù)文件Torque1(x)的建立:

function[Mmax,Gamma]=Torque1(lota,Pn,n,G,d,N0)

Mmax=9.55*10^6*(Pn/n)*lota；

Gamma=5836*Mmax*10^(-3)/(G*d^4*(1-N0^4))*10^5；

在命令窗口中輸入實參lota、Pn、n、G、d、N0的具體數(shù)值，然后調(diào)用計算最大扭矩與許用轉(zhuǎn)角的函數(shù)文件Torque1(lota,Pn,n,G,d,N0)，便可得出計算結(jié)果[6-8].

函數(shù)程序截圖如圖3所示，運算結(jié)果截圖如圖4所示.

圖3 函數(shù)程序截圖

圖4 運算結(jié)果截圖

相應(yīng)的程序運算精確結(jié)果為：Mmax=1.349 9×106N·mm

γ=0.1251°/mm

3 攪拌軸的強度計算

由于攪拌釜體內(nèi)存在不規(guī)則流動，故攪拌機徑向受力，使攪拌軸在承擔扭矩的同時還承受彎矩作用.因此，對于常壓容器中攪拌軸的強度計算，應(yīng)綜合考慮扭矩和彎矩聯(lián)合作用下的受力情況[4,5].

3.1 彎扭聯(lián)合強度計算

在彎扭組合的當量扭矩作用下，攪拌軸的強度條件為：

M=MR+MA

式中：Mte為彎扭組合的當量扭矩，N·mm；M為軸上彎矩總和，N·mm；MA為由軸上推力引起作用于軸上的彎矩，N·mm；MR為由軸上推力引起作用于軸上的彎矩，N·mm；Mn為軸上扭矩，N·mm；PN為電機額定功率，kw；η1為傳動裝置的效率，η1=95%；n為軸的轉(zhuǎn)速，r/min.

上述計算參數(shù)見表1.

表1 強度計算參數(shù)

0Cr18Ni9(304)材料的 [σ]=137 MPa，[τ]=(0.5～0.6)，[σ]=(68.5～82.2)MPa

計算得τ=49MPa[τ]

3.2 彎扭聯(lián)合強度的MATLAB計算

相應(yīng)內(nèi)容的M函數(shù)文件建立與計算過程[6-8]:

第i個攪拌器上流體徑向力：

function[Fhi]=force1(K1,Mn,Dj)

Fhi=K1*Mn*10^3/(3/8*Dj)；

函數(shù)程序截圖如圖5所示，運算結(jié)果截圖如圖6所示.

圖5 函數(shù)程序截圖

圖6 運算結(jié)果截圖

相應(yīng)的程序運算精確結(jié)果為：

Fhi=1.0908×103N

在命令窗口中輸入實參K1、Mn、Dj的具體數(shù)值，然后調(diào)用計算第i個攪拌器上流體徑向力的函數(shù)文件force1(K1,Mn,Dj)，便可得出計算結(jié)果.

攪拌軸與各層圓盤組合質(zhì)量偏心引起的離心力：

function[F]=force2(nk,Mw1,n)

F=((pi^2)/9)*Mw1*n^2*floor(exp(1))*10^(-5)/(1-(n/nk)^2)；

函數(shù)程序截圖如圖7所示，運算結(jié)果截圖如圖8所示.

相應(yīng)的程序運算精確結(jié)果為：

Fhi=1.0908×103N

在命令窗口中輸入實參nk、Mw1、n的具體數(shù)值，然后調(diào)用計算攪拌軸與各層圓盤組合質(zhì)量偏心引起的離心力的函數(shù)文件force2(nk,Mw1,n)，便可得出計算結(jié)果.

攪拌軸及攪拌器組合重心至傳動側(cè)軸承的距離：

function[Le]=distance1(ML1,L1,Mw)

s=sum(Li.*Mi)；

Le=(s+(ML1*L1)/2/Mw)；

圖7 函數(shù)程序截圖

圖8 運算結(jié)果截圖

在命令窗口中輸入每個函數(shù)相應(yīng)實參ML1，L1，Mw的具體數(shù)值，然后調(diào)用計算攪拌軸及攪拌器組合重心至傳動側(cè)軸承的距離的函數(shù)文件function[Le]，按矩陣格式輸入Li，Mi，便可得出計算結(jié)果.

函數(shù)程序截圖如圖9所示，運算結(jié)果截圖如圖10所示.

圖9 函數(shù)程序截圖

圖10 運算結(jié)果截圖

相應(yīng)的程序運算精確結(jié)果為：

按彎扭聯(lián)合強度，相應(yīng)MATLAB計算.函數(shù)程序截圖如圖11所示，運算結(jié)果截圖如圖12所示.

圖11 函數(shù)程序截圖

圖12 運算結(jié)果截圖

相應(yīng)的程序運算精確結(jié)果為：τ=49.046 MPa

4 結(jié) 論

基于MATLAB的機械設(shè)計方法利于機械設(shè)計的優(yōu)化，MATLAB軟件進行機械零件的強度計算，不用編寫大量的算法程序，為驗證攪拌軸強度計算的精確程度，采用了計算機模擬測試，測試出的結(jié)果與計算結(jié)果完全吻合，提高了計算效率和計算的精度.同時，根據(jù)實際的設(shè)計和計算經(jīng)驗證明，用MATLAB計算得出的結(jié)果與實際完全一致，是完全可以信賴的優(yōu)化計算方法，對工程實際具有一定的指導(dǎo)作用.