輪胎微波硫化圓柱諧振腔的設(shè)計及數(shù)值模擬研究

李 濤，孫 斌，陳海龍，張茂東，胡冰濤，李慶領(lǐng)

（青島科技大學(xué) 機(jī)電工程學(xué)院，山東 青島 266061）

硫化是輪胎等橡膠制品在制造加工過程中必不可少的一道工序，對橡膠制品的實際應(yīng)用具有重要意義。微波硫化是一種新型橡膠硫化方式，其硫化快、效率高[1-2]。微波加熱最大的缺陷是橡膠內(nèi)外溫度不均勻，因此，設(shè)計一種專用于輪胎硫化的微波加熱裝置具有重要的現(xiàn)實意義。

國內(nèi)外學(xué)者對諧振腔腔體設(shè)計進(jìn)行了相關(guān)研究。巨漢基等[3]對影響微波爐性能的各個因素進(jìn)行了仿真研究。W.Cha-um等[4]利用數(shù)值和試驗方法對比研究了樣品位置、尺寸和微波功率對加熱效率的影響。孫鵬等[5]利用有限元分析軟件仿真了饋口位置、樣品大小、樣品位置和多饋口激勵時對微波利用率的影響。Y.Huang等[6]運(yùn)用實驗輔助理論的方法分析了規(guī)則矩形諧振腔的場均勻性問題。V.Sebera等[7]通過有限元法模擬模式攪拌器對多模腔內(nèi)微波加熱均勻性的影響。楊繼孔等[8]對高功率矩形微波反應(yīng)器加熱效率和均勻性仿真進(jìn)行了研究。曹湘琪等[9-10]對圓柱形微波加熱器的效率和均勻性仿真進(jìn)行了研究，并利用HFSS軟件仿真圓柱形微波加熱器饋口位置和長度、內(nèi)筒半徑和高度以及負(fù)載厚度對微波吸熱效率的影響。

迄今為止，國內(nèi)外學(xué)者對用于輪胎微波硫化諧振腔的相關(guān)研究較少，僅有少數(shù)關(guān)于輪胎微波硫化矩形諧振腔的研究，而圓柱諧振腔具有品質(zhì)因數(shù)高、調(diào)諧方便、結(jié)構(gòu)堅固和易于加工制作等優(yōu)點。本工作對輪胎微波硫化圓柱諧振腔進(jìn)行設(shè)計，并對圓柱諧振腔內(nèi)輪胎微波硫化過程進(jìn)行數(shù)值模擬研究。

1 諧振腔尺寸理論計算

1.1 設(shè)計原理

對于圓柱形諧振腔，其諧振頻率的計算公式[11]如下：

式中，fr為諧振腔的諧振頻率；R為圓柱形諧振腔半徑；l為圓柱形諧振腔的高度；m，n，p為模式標(biāo)號，取值為正整數(shù)（0，1，2，3，4…），每一種組合稱為一種工作模式，其組合數(shù)稱為模式數(shù)；μmn為m階貝塞爾函數(shù)的導(dǎo)函數(shù)第n個根植。以上參數(shù)均采用國際單位制。

諧振頻率只有滿足公式（2）[12]時的模式才能在腔體內(nèi)生成。

式中，f0為工作頻率；Δf為頻寬。

確定工作頻率和頻寬后，對于一定半徑和高度的圓柱形諧振腔，可以計算該圓柱形諧振腔的模式數(shù)、簡并比。通過計算不同尺寸圓柱形諧振腔的模式數(shù)、簡并比，并從中選擇模式數(shù)較大、簡并比為0的圓柱形諧振腔（此時諧振腔內(nèi)電磁場分布更為均勻且損耗低）。

1.2 設(shè)計過程

（1）確定工作頻率和頻寬。目前，中國用于工業(yè)加熱設(shè)備的微波固定專用頻率為L波段[（915±25） MHz]和S波段[（2 450±50） MHz]。915 MHz型的設(shè)備正受到重視，因為橡膠在915和2 450 MHz下的吸收性能并無明顯差異，但915 MHz卻有透射能力強(qiáng)、投資成本低、功率利用率高、磁控管壽命長、可達(dá)功率高和維修較方便等優(yōu)點。

（2）確定圓柱形諧振腔的尺寸。模擬選用輪胎為7.50R20內(nèi)支撐外胎，確定圓柱形諧振腔的最小半徑為472 mm，最小高度為230 mm。設(shè)定R為472～682 mm，ΔR為15 mm；l為230～530 mm，Δl為30 mm。

（3）計算各尺寸腔體的模式數(shù)、簡并比。根據(jù)公式（2）計算不同尺寸腔體的模式數(shù)、簡并比。

（4）初步確定圓柱形諧振腔尺寸。對計算結(jié)果進(jìn)行篩選，選擇模式數(shù)大于或者等于9且簡并比為0的圓柱形諧振腔，初步確定圓柱形諧振腔尺寸及其對應(yīng)的模式數(shù)和簡并比如表1所示。

表1 模式數(shù)大于8的圓柱形諧振腔尺寸及其對應(yīng)的模式數(shù)和簡并比

2 數(shù)值模型與物理模型

2.1 模型簡化

為節(jié)省計算時間和降低CPU使用率，對模型進(jìn)行合理假設(shè)，具體如下：

假設(shè)1：構(gòu)成輪胎的橡膠材料均勻且各向同性；

假設(shè)2：介電性能不變；

假設(shè)3：輪胎和空氣的初始溫度均勻；

假設(shè)4：傳質(zhì)過程忽略不計；

假設(shè)5：空氣與輪胎之間的化學(xué)反應(yīng)忽略不計；

假設(shè)6：空氣的介電損耗忽略不計；

假設(shè)7：波導(dǎo)和腔體的材質(zhì)為銅；

假設(shè)8：空氣中熱傳導(dǎo)忽略不計；

假設(shè)9：輪胎邊界完全絕緣；

假設(shè)10：端口由Y方向的橫向電磁駐波場激勵。

2.2 物理模型

將7.50R20內(nèi)支撐外胎簡化為5層，分別為胎面、胎肩、胎側(cè)、氣密層和帶束層，輪胎輪廓如圖1所示。根據(jù)工作頻率（915 MHz）選擇BJ8型號波導(dǎo)，波導(dǎo)長度定為100 mm，放置在圓柱形諧振腔側(cè)面正中位置。輪胎、諧振腔和波導(dǎo)的相對位置如圖2所示。

圖1 7.50R20輪胎輪廓示意

圖2 輪胎、諧振腔和波導(dǎo)相對位置示意

2.3 控制方程

COMSOL Multiphysics軟件用于輪胎微波硫化過程數(shù)值模擬，該數(shù)值模擬過程耦合了電磁方程和傳熱方程。COMSOL Multiphysics軟件使用了有限元方法（FEM）的數(shù)值計算模型，采用麥克斯韋方程計算數(shù)值模擬過程中的電磁傳播。電場波控制方程如下：

式中，μr為相對磁導(dǎo)率；k0為自由空間的波數(shù)；εr為相對介電常數(shù)；σ為電導(dǎo)率；ω為角頻率；ε0為真空介電常數(shù)（8.85×10-12F·m-1）；E為電場強(qiáng)度，V·m-1。

k0根據(jù)公式（4）計算：

式中，c0為真空光速。

介電常數(shù)為復(fù)數(shù)形式，表達(dá)式如下：

式中，ε*為復(fù)介電常數(shù)，F(xiàn)·m-1；ε′為介電常數(shù)實部；ε″為介電常數(shù)虛部，F(xiàn)·m-1；j=（-1）1/2。

微波場傳熱方程用傅里葉能量平衡公式來耦合，表達(dá)式如下：

式 中，ρ為 密 度，kg·m-3；cp為 常 壓 比 熱 容，J·kg-1·K-1；T為熱力學(xué)溫度，K；k為熱導(dǎo)率，W·m-1·K-1；f為頻率；δ為介電損耗角；Q為體積熱源。

2.4 網(wǎng)格質(zhì)量

采用物理控制網(wǎng)格構(gòu)建整個幾何體。網(wǎng)格單元尺寸設(shè)置為細(xì)化。為得到準(zhǔn)確的結(jié)果，最大網(wǎng)格單元尺寸精細(xì)到到微波波長的1/10。輪胎構(gòu)建網(wǎng)格單元（見圖3）數(shù)目為153 633。網(wǎng)格單元質(zhì)量評價如圖4所示。根據(jù)經(jīng)驗，當(dāng)單元質(zhì)量值大于0.6時，數(shù)值模擬結(jié)果是可靠的。

圖3 網(wǎng)格單元

圖4 網(wǎng)格單元質(zhì)量評價

3 試驗研究

3.1 微波加熱裝置

試驗選用Milestone Microwave Laboratory Systems系列中的High Performance Microwave Digestion Unit MLS 1200 mega作為加熱橡膠的微波加熱器（工作頻率為2 450 MHz），微波加熱器如圖5所示。

圖5 微波加熱器示意

微波加熱器的主要部件包括諧振腔（微波加熱腔體）、磁控管、功率放大器、波導(dǎo)、變壓器、電源、爐門、控制面板和通風(fēng)裝置。微波加熱器接通電源后，磁控管開始運(yùn)轉(zhuǎn)，其輸出的能量經(jīng)功率放大器放大后，再通過波導(dǎo)傳輸，最后由耦合窗耦合到諧振腔中，微波在加熱腔體內(nèi)來回反射從而起到加熱物料的作用。該微波加熱器的優(yōu)勢在于能夠在一個較大范圍內(nèi)對加熱功率和時間進(jìn)行設(shè)定，且可以進(jìn)行功率、時間組合加熱條件的設(shè)定。

3.2 溫度測量裝置

通常情況下，采用熱電偶、溫度計等對物料進(jìn)行測溫，但無論是熱電偶還是溫度計中都含有金屬，將金屬放置在微波環(huán)境中，會發(fā)生打火現(xiàn)象，這使得傳統(tǒng)的測溫方式束手無措。測溫儀器在微波加熱橡膠環(huán)境中必須滿足以下兩點要求：

（1）測溫元件不受微波強(qiáng)電磁場的影響，能夠精確無誤測量物料的溫度；

（2）測溫元件中不含金屬材料，并且使用的材料能夠耐200 ℃以上溫度。

經(jīng)過篩選，光纖溫度傳感器和紅外熱成像儀滿足以上兩點要求。試驗中選用光纖溫度傳感器實時測量膠片在微波加熱時的溫度，選用紅外熱成像儀測量膠片加熱后的最終表面溫度。試驗中所用光纖溫度在線檢測裝置如圖6所示，其通過數(shù)據(jù)線與筆記本電腦連接，測溫端與橡膠試樣連接，另一端與光纖溫度在線監(jiān)測裝置連接。傳感器技術(shù)參數(shù)如下：

圖6 光纖溫度在線監(jiān)測裝置示意

電源電壓 220 V，測量范圍 －50～＋200 ℃，測量精度 ±0.5 ℃FS，分辨率 0.1 ℃，采樣頻率 1 Hz，數(shù)字接口 RS485，波特率 9 600 bps，工作溫度 －40～＋75 ℃。

3.3 橡膠塊測溫點布置

對不同尺寸的橡膠塊（如表2所示）微波加熱。為得到微波加熱過程中橡膠塊的溫度分布規(guī)律，在橡膠塊內(nèi)布置9個測溫點并對測溫點標(biāo)號，如圖7所示。

表2 微波諧振腔、波導(dǎo)、聚四氟乙烯板和橡膠塊尺寸 mm

圖7 橡膠塊測溫點布置示意

3.4 試驗驗證

使用微波加熱裝置對橡膠塊加熱并采用光纖溫度在線檢測裝置測量橡膠塊內(nèi)測溫點的溫度。使用COMSOL Multiphysics模擬橡膠塊的微波加熱過程。將幾個測溫點的試驗結(jié)果和數(shù)值模擬結(jié)果進(jìn)行對比，如圖8所示。數(shù)據(jù)分析表明：試驗結(jié)果和模擬結(jié)果中測溫點的溫升規(guī)律相似，模擬數(shù)據(jù)與試驗數(shù)據(jù)密切相符；所有測溫點的模擬數(shù)據(jù)和試驗數(shù)據(jù)之間的偏差在5%以內(nèi)，驗證了數(shù)值模擬模型的可靠性。

圖8 測溫點的試驗結(jié)果和數(shù)值模擬結(jié)果

4 結(jié)果與分析

微波激發(fā)頻率為915 MHz，模式為TE10，微波輸出功率為10 kW。受熱的輪胎橡膠初始溫度為20 ℃，加熱時間為1 800 s。為使7.50R20輪胎加熱均勻，讓其在諧振腔內(nèi)進(jìn)行旋轉(zhuǎn)。對選定的10組尺寸圓柱形諧振腔內(nèi)的輪胎微波硫化過程進(jìn)行數(shù)值模擬，得到微波作用（加熱）720 s時的輪胎溫度云圖（見圖9）。

從圖9可以看出，隨著模式數(shù)增加，輪胎微波硫化的均勻性有所不同，即使同一模式數(shù)的諧振腔，輪胎微波硫化的均勻性也有所不同。其中5#諧振腔和6#諧振腔內(nèi)輪胎溫度分布較均勻，最大溫差分別為36.3和48.7 ℃；1#，4#和7#諧振腔內(nèi)輪胎溫度分布相對不均勻，最大溫差分別為76.9，70.0和105.1 ℃。

圖9 微波加熱720 s時的輪胎溫度分布云圖

模擬得到微波加熱720 s時1#—10#諧振腔內(nèi)輪胎體平均溫度、回波損耗S參數(shù)和最大溫差，通過S=10log（P反射/P入射）對S參數(shù)計算，得到各諧振腔的微波利用率。將微波加熱720 s時1#—10#諧振腔內(nèi)輪胎體平均溫度、微波利用率、回波損耗S參數(shù)和最大溫差統(tǒng)計于表3中。

根據(jù)表3數(shù)據(jù)作圖，探索微波利用率和輪胎體平均溫度的關(guān)系，如圖10所示。

從表3和圖10可以看出，輪胎體平均溫度變化趨勢與微波利用率變化趨勢相同，即隨著微波利用率增大而增大，反之亦然，即微波利用率越大，輪胎體平均溫度越高。這是因為微波利用率越高，則在相同的加熱時間內(nèi)，輪胎獲得的能量越多，輪胎的加熱體積就越大，輪胎體平均溫度就越高。為了更好地得到微波利用率和輪胎體平均溫度的關(guān)系，對輪胎體平均溫度和微波利用率進(jìn)行線性擬合，結(jié)果如圖11所示。

圖10 圓柱形諧振腔的微波利用率及其加熱720 s時的輪胎體平均溫度

表3 模式數(shù)大于8的圓柱形諧振腔性能參數(shù)

從圖11可以看出，擬合曲線與輪胎體平均溫度和微波利用率兩組數(shù)據(jù)切合程度較高，輪胎體平均溫度隨著微波利用率增大而增大。通過擬合得到y(tǒng)＝80.107x＋22.047，可決系數(shù)R2＝0.990 52，這表明擬合曲線能夠較好地表征輪胎體平均溫度和微波利用率之間的關(guān)系，且輪胎體平均溫度與微波利用率呈現(xiàn)良好的一次函數(shù)關(guān)系。

圖11 輪胎體平均溫度和微波利用率關(guān)系的擬合曲線

根據(jù)表3中輪胎體平均溫度和最大溫差數(shù)據(jù)作圖，結(jié)果如圖12所示。

圖12 圓柱諧形振腔的最大溫差及其加熱720 s時的輪胎體平均溫度

從圖12可以看出：1#，2#，3#，4#，7#和10#諧振腔內(nèi)輪胎體平均溫度均大于90 ℃，微波利用率較大；1#，3#，4#，7#，9#和10#諧振腔內(nèi)輪胎最大溫差大于60 ℃且加熱相對不均勻；2#諧振腔的輪胎最大溫差較小，輪胎微波硫化不均勻；5#，6#，8#和9#諧振腔內(nèi)輪胎體平均溫度在80～90 ℃，低于其他諧振腔，微波利用率較小；8#和9#諧振腔的輪胎最大溫差大于5#和6#諧振腔且輪胎微波硫化不均勻；5#諧振腔內(nèi)輪胎體平均溫度大于6#諧振腔內(nèi)輪胎體平均溫度且輪胎最大溫差小于6#諧振腔內(nèi)的輪胎最大溫差，雖然5#諧振腔內(nèi)輪胎體平均溫度較小、微波利用率較低，但輪胎最大溫差最小，輪胎微波硫化也最為均勻。為取得更好的加熱效果，確定5#諧振腔的尺寸為最終圓柱形諧振腔的尺寸，即圓柱形諧振腔半徑為652 mm、高度為470 mm。

5 結(jié)論

（1）隨著圓柱形諧振腔模式數(shù)增大，輪胎微波硫化的均勻性有所不同，即使同一模式數(shù)諧振腔，輪胎微波硫化的均勻性也有所不同。

（2）輪胎體平均溫度變化趨勢與微波利用率變化趨勢相同，呈現(xiàn)良好的一次函數(shù)關(guān)系。

（3）確定5#諧振腔的尺寸為最終圓柱形諧振腔尺寸，即圓柱形諧振腔半徑為652 mm，高度為470 mm。