楊建興 丁治豪，2

（1.西安近代化學(xué)研究所；2.西北大學(xué)化工學(xué)院）

多年來推進(jìn)劑的雙基球形藥制備多采用內(nèi)溶法成球工藝［1～3］，該工藝是將硝化纖維素或吸收藥、水、乙酸乙酯加入球釜內(nèi)，通過攪拌漿葉的強(qiáng)力攪拌，物料被溶解成具有一定黏度的高分子溶液，隨之被攪拌粉碎成細(xì)小的液滴，液滴與水不相溶，在表面張力的作用下，形成球形，再通過升溫，使乙酸乙酯汽化脫溶，液滴脫溶變硬形成球形藥粒［2～4］。但是通過攪拌使高分子溶液在不相溶的水中形成球粒的過程存在球粒粒度分布不均、耗時(shí)較長等問題，成品質(zhì)量不易控制且成型效率較低。

文獻(xiàn)［5，6］在球形藥制備工藝中嘗試將雙基吸收藥與乙酸乙酯混合塑化后采用水下切割的方式造粒，再通過旋流修圓，最后進(jìn)入脫溶釜內(nèi)完成乙酸乙酯的汽化脫溶。 在上述工藝優(yōu)化過程中，塑化后的雙基吸收藥與乙酸乙酯混合物通過螺桿推動(dòng)由?？讛D出后進(jìn)入切粒水箱切割造粒，螺桿擠出裝置對(duì)假塑性混合物料的擠壓效應(yīng)決定了球形藥的粒度均一性。

筆者采用CFD仿真技術(shù)模擬了物料在擠出裝置中的流動(dòng)過程，獲取了不同時(shí)刻物料在裝置中的速度場(chǎng)與壓力場(chǎng)分布，以此分析螺桿旋轉(zhuǎn)過程中對(duì)物料的擠壓效應(yīng)和擠壓孔處的出料效果。

1 數(shù)值模型

1.1 幾何模型及網(wǎng)格劃分

螺桿擠出裝置結(jié)構(gòu)如圖1所示， 螺桿段長度為320 mm，螺紋螺距為26 mm，螺桿段上套筒內(nèi)徑為20 mm，出藥模具安裝于螺桿段出口處，模具的擠壓孔孔徑為1.25 mm，兩環(huán)分布共60個(gè)。 由于幾何結(jié)構(gòu)較為復(fù)雜， 采用四面體+邊界層網(wǎng)格劃分，網(wǎng)格總數(shù)為410萬，螺桿段設(shè)為動(dòng)區(qū)域，其他部分設(shè)為靜區(qū)域（圖2）。

圖1 螺桿擠出裝置幾何模型

圖2 螺桿擠壓裝置不同部位網(wǎng)格劃分

1.2 發(fā)射藥物性參數(shù)

物料混合物由SC-11發(fā)射藥與乙酸乙酯經(jīng)塑化后形成（配比為1:0.9），密度為1 600 kg/m3，采用Thermo HAAKE MARS 60旋轉(zhuǎn)流變儀測(cè)量其流變性能，25 ℃下發(fā)射藥剪切應(yīng)力與剪切速率的關(guān)系曲線如圖3所示， 觀察可知其符合假塑性流體的流動(dòng)規(guī)律。

圖3 發(fā)射藥剪切應(yīng)力與剪切速率的關(guān)系曲線

采用Ostwald-de Wale冪律方程對(duì)圖中SC-11發(fā)射藥剪切應(yīng)力與剪切速率的關(guān)系曲線進(jìn)行擬合，得出物料表觀黏度與剪切速率的變化方程：

式中 ηa——表觀黏度，Pa·s；

γa——剪切速率，s-1。

其中，曲線擬合度R2=0.9998。

1.3 數(shù)學(xué)模型

物料混合物為非牛頓流體，遵循質(zhì)量、動(dòng)量、能量守恒［7］，守恒方程如下：

式中 Cp——定壓比熱容，kJ/（kg·K）；

k——材料導(dǎo)熱系數(shù)，W/（m·K）；

u——流體速度，m/s；

?——哈密頓算子；

ρ——流體密度，kg/m3；

σ——表面應(yīng)力，Pa；

τ——偏應(yīng)力張量，Pa。

模擬中根據(jù)物料特性和螺桿擠出時(shí)的工藝條件，假設(shè)如下：物料物性不隨溫度變化，無化學(xué)反應(yīng)，為不可壓縮層流；壁面無滑移條件；不計(jì)重力和慣性力作用。

1.4 邊界條件及求解方法

流場(chǎng)動(dòng)區(qū)域內(nèi)螺桿轉(zhuǎn)速為20 r/min， 邊界條件設(shè)為質(zhì)量入口、壓力出口條件，物料入口質(zhì)量流量為10 kg/h，出口壓強(qiáng)為大氣壓。 求解方程采用有限容積法， 壓力-速度的耦合求解采用SIMPLE算法，對(duì)流項(xiàng)離散采用二階迎風(fēng)差分。 考慮到物料黏度會(huì)隨剪切速率發(fā)生變化，為了觀測(cè)瞬態(tài)擠出效應(yīng)，模擬中首先采用MRF（多重參考系法）計(jì)算到收斂（收斂準(zhǔn)則為RMS，收斂殘差定為10-3）， 隨后采用滑移網(wǎng)格模型進(jìn)行瞬態(tài)模擬，時(shí)間步長設(shè)為10-3s，每個(gè)時(shí)間步迭代20次，收斂殘差同樣為10-3。

2 結(jié)果及分析

圖4為螺桿擠壓裝置內(nèi)物料的時(shí)均速度場(chǎng)分布，由圖可知，當(dāng)物料進(jìn)入料斗后即開始旋流，直到進(jìn)入擴(kuò)張錐內(nèi)流動(dòng)方向變?yōu)檩S向（z軸方向），在擠壓孔出口處基本是垂直于出口平面流出的（圖5）。 由擠壓孔出口速度分布云圖（圖6）可以看到，每個(gè)擠壓孔出料速度基本一致，可以推知設(shè)計(jì)的出藥模具（出口擠壓孔孔徑為1.25 mm、兩環(huán)分布共60個(gè)）結(jié)構(gòu)合理，符合生產(chǎn)預(yù)期。

圖4 螺桿擠壓裝置內(nèi)部速度場(chǎng)

圖5 擠壓孔出口速度分布云圖

圖6 擠壓孔出口速度分布云圖

圖7為物料在螺桿擠壓裝置內(nèi)的壓強(qiáng)場(chǎng)分布結(jié)果。 物料進(jìn)入料斗后在慣性作用下克服黏性力而螺旋下滑，壓強(qiáng)逐步降低。 當(dāng)物料進(jìn)入螺桿推進(jìn)部分后，在螺桿的推動(dòng)下向擠壓孔出口方向移動(dòng)的過程中逐步升壓， 在前端第1～2個(gè)螺距之間壓強(qiáng)達(dá)到最大值（圖8），隨即加速減壓進(jìn)入擴(kuò)張錐，由出藥模具的擠壓孔擠出。 在料斗后方螺桿根部區(qū)域形成了負(fù)壓區(qū)（圖9），部分物料由料斗被吸入該區(qū)域，沿螺旋導(dǎo)葉反向運(yùn)動(dòng)到螺桿根部再折向移動(dòng)。 在負(fù)壓區(qū)物料有較大的反向速度（圖10）和逆向流量，對(duì)螺桿的功耗和振動(dòng)都會(huì)有一定影響，可以通過調(diào)整料斗位置和料斗內(nèi)螺旋導(dǎo)板的旋向加以改善。

圖7 螺桿擠壓裝置內(nèi)的壓強(qiáng)場(chǎng)分布

圖8 螺桿擠壓裝置內(nèi)的最大壓力位置

圖9 螺桿擠壓裝置內(nèi)的負(fù)壓區(qū)

圖10 中心截面上軸向（z方向）速度分布云圖

圖11為螺桿旋轉(zhuǎn)一周，每旋轉(zhuǎn)60°時(shí)裝置內(nèi)中心截面上物料的瞬時(shí)速度場(chǎng)分布云圖。 圖12、13為相應(yīng)的壓強(qiáng)場(chǎng)云圖和擠壓孔出口的速度分布云圖。 由圖12、13可以看出，相比時(shí)均值，盡管擠壓孔出口的出藥速度與螺桿段定性相似，但在具體數(shù)值上還是會(huì)有一定的變化，特別是在旋轉(zhuǎn)一圈的結(jié)束和開始的階段，裝置內(nèi)速度會(huì)有較大的變化，這應(yīng)該與螺桿的循環(huán)回轉(zhuǎn)有關(guān)。 擠壓孔出口處速度分布雖然會(huì)有一定的變化，但仍保持均勻分布，且變化量很小，因此出口流量基本不變，與時(shí)均值相比偏差在10-5以上（表1），該結(jié)果說明出口質(zhì)量流量基本不隨時(shí)間發(fā)生變化，進(jìn)而說明該裝置可以保證物料連續(xù)、穩(wěn)定地?cái)D出。

表1 螺桿旋轉(zhuǎn)一周出口質(zhì)量流量變化

圖11 不同瞬時(shí)裝置內(nèi)中心截面上的速度分布云圖（螺桿旋轉(zhuǎn)1圈）

圖12 不同瞬時(shí)裝置內(nèi)中心截面上的壓強(qiáng)場(chǎng)云圖（螺桿旋轉(zhuǎn)1圈）

圖13 不同瞬時(shí)擠壓孔出口速度分布云圖（螺桿旋轉(zhuǎn)1圈）

圖14顯示了裝置內(nèi)部螺桿表面的受力情況，由圖14a可以看出，螺桿表面壓強(qiáng)變化較為一致，沿螺桿端部方向壓強(qiáng)逐漸增大，料斗后側(cè)處于負(fù)壓區(qū)；由圖14b可以看出，切應(yīng)力分布較為復(fù)雜，軸段部分與導(dǎo)葉部分差異較大，受力主要在軸端部分。 對(duì)比圖14a、b可以看出，螺桿表面靜壓強(qiáng)值域范圍遠(yuǎn)大于切應(yīng)力，最大值相差2個(gè)數(shù)量級(jí)。 結(jié)合圖8、9可知，螺桿旋轉(zhuǎn)過程中受力不均勻， 可能會(huì)導(dǎo)致軸段和導(dǎo)葉發(fā)生變形及振動(dòng)的情況。

圖14 螺桿表面受力情況

3 結(jié)束語

應(yīng)用FLUENT模擬計(jì)算了雙基發(fā)射藥螺桿擠壓成型過程，模擬結(jié)果顯示設(shè)計(jì)工況下采用文中設(shè)計(jì)的螺桿結(jié)構(gòu)可以保證雙基發(fā)射藥穩(wěn)定地?cái)D壓成型，成型截面較為規(guī)則，不會(huì)影響后續(xù)的切粒工藝。 運(yùn)行中螺桿根部存在負(fù)壓區(qū)，部分發(fā)射藥物料被吸入該區(qū)域， 逆向流動(dòng)再折向返流，對(duì)螺桿的功耗和振動(dòng)都會(huì)有一定影響，可以通過調(diào)整料斗位置和料斗內(nèi)螺旋導(dǎo)板的旋向加以改善。螺桿表面受力以壓力為主，靜壓強(qiáng)的最大值比表面切應(yīng)力大了2個(gè)數(shù)量級(jí)。