離心攪拌機(jī)葉片流固耦合應(yīng)力場分析

李紅軍，楊 康，張 弛，陳 偉

離心攪拌機(jī)葉片流固耦合應(yīng)力場分析

李紅軍1,2，楊 康1，張 弛1，陳 偉1,2

（1.武漢紡織大學(xué) 機(jī)械工程與自動(dòng)化學(xué)院，湖北 武漢 430200；2.工業(yè)雷管智能裝配湖北工程研究中心，湖北 武漢 430200）

為研究離心攪拌機(jī)葉片工作時(shí)應(yīng)力分布問題，本文建立了一種離心攪拌機(jī)葉片流固耦合應(yīng)力場分析模型，利用ANSYS-CFX流體分析軟件將流體分析與結(jié)構(gòu)分析進(jìn)行耦合，對(duì)離心攪拌機(jī)葉片在高轉(zhuǎn)速工況下的應(yīng)力場進(jìn)行分析，得到流體載荷在葉片內(nèi)部產(chǎn)生的等效應(yīng)力分布圖及整體葉片耦合后的最大形變量和工作應(yīng)力場。最終通過分析應(yīng)力場結(jié)果得到了離心攪拌機(jī)合適的工作轉(zhuǎn)速，為工業(yè)生產(chǎn)提供參考建議。

攪拌葉片；CFX；流固耦合；應(yīng)力分析

攪拌葉片作為離心攪拌機(jī)的核心部件，其強(qiáng)度和可靠性直接影響著攪拌機(jī)的性能。分析攪拌機(jī)葉片受力情況的傳統(tǒng)方法通常是將攪拌機(jī)功率轉(zhuǎn)化為作用在葉片上的均布載荷，但在實(shí)際攪拌過程中葉片與流體之間作用載荷并不是均勻分布的，從而造成較大誤差。所以流固耦合的分析方法受到越來越多學(xué)者關(guān)注，它的優(yōu)點(diǎn)在于能夠分析出兩種不同工作場下的模擬實(shí)驗(yàn)數(shù)據(jù)。張衎[1]等分析了流固耦合情況下葉輪流體壓力數(shù)值和離心載荷對(duì)葉輪強(qiáng)度的影響；蔡琳[2]等利用CFX-ANSYS耦合分析了JDY型攪拌機(jī)葉片在水中的形變與受力情況。趙悟[3]等分析了雙立軸攪拌機(jī)在不同工況下的受力，求出攪拌機(jī)構(gòu)的形變、應(yīng)力和固有頻率。但是在攪拌機(jī)葉片受不同轉(zhuǎn)速工況下的流固耦合應(yīng)力場的分析特性較少，為此，本文使用了單向流固耦合的方法研究了在不同轉(zhuǎn)速下葉片的形變與受力。給出了攪拌機(jī)工作轉(zhuǎn)速適宜范圍，得出了應(yīng)力集中部位，為攪拌機(jī)葉片的疲勞分析和設(shè)計(jì)優(yōu)化提供參考。

1 建立攪拌機(jī)葉片整體模型

利用Inventor三維建模軟件建立了葉片的整體式結(jié)構(gòu)模型。由上至下帶有三層攪拌刮刀，作用是處理物料與攪拌腔內(nèi)壁粘連部分；中部有圓錐形散料裝置，作用是使物料能夠快速進(jìn)入到攪拌腔內(nèi)。同時(shí)底部三個(gè)攪拌葉片呈120°均勻分布在攪拌盤周圍，主要起到攪拌混合物料的作用。攪拌葉片工作半徑為346mm，豎直方向最大高度為218.37mm。具體分布方式如圖1所示。

1散料裝置；2第一層刮刀；3第二層刮刀；4第三層刮刀；5攪拌葉片；6攪拌軸

2 流固耦合分析理論方法

流固耦合是一種基于流體力學(xué)和固體力學(xué)相結(jié)合的研究方法，可分為單向耦合和雙向耦合兩種。單向耦合只考慮流場內(nèi)壓力變化對(duì)固體場的影響，而忽略固體場內(nèi)的形變反過來對(duì)流場的影響；雙向耦合既要考慮流場對(duì)固體場的影響，也不能忽視固體形變對(duì)流體的影響[4]。由于本文主要研究流場內(nèi)壓力對(duì)攪拌機(jī)葉片的影響，且考慮到材料特性，葉片形變量較小，對(duì)流場的反作用可忽略,因此考慮到以上因素本文選擇使用單向流固耦合的方法來對(duì)攪拌機(jī)葉片強(qiáng)度進(jìn)行分析計(jì)算。

2.1 流體的控制方程

任何流動(dòng)問題都必須滿足守恒定律，即對(duì)于此問題中流體的控制方程同樣需要遵循基本的質(zhì)量、動(dòng)量、能量守恒方程。攪拌過程屬于一種復(fù)雜的高速旋轉(zhuǎn)流動(dòng)，由雷諾數(shù)方程可知對(duì)一般管道流體而言當(dāng)雷諾數(shù)大于2300時(shí)流體會(huì)發(fā)生湍流現(xiàn)象，而處于湍流模式下的不可壓縮流體只需要考慮質(zhì)量守恒和動(dòng)量守恒方程[5],不考慮能量之間的傳遞，具體質(zhì)量方程、動(dòng)量方程如下所示。

（1）質(zhì)量守恒方程：

式中ρ—流體密度；v—流體速度；t—時(shí)間；?—哈密頓算子

（2）動(dòng)量守恒方程：

式中τ—切向應(yīng)力；p—二階應(yīng)力張量；F—質(zhì)量力；g—重力加速度

2.2 單向流固耦合方程

由于流體運(yùn)動(dòng)所引發(fā)的固體振動(dòng)和位移方程為：

對(duì)于流體和固體作用的接觸面需要滿足流體與固體的應(yīng)力、位移、溫度和熱量交換等參數(shù)守恒，具體守恒方程如下：

式中τ—所受應(yīng)力；r—位移；T—溫度；q—熱量；下標(biāo)f—流體；s—固體。

3 攪拌葉片流固耦合仿真分析

3.1 仿真前處理

將Inventor內(nèi)建立好的離心攪拌機(jī)整體三維模型以STP格式保存導(dǎo)入到ANSYS CFX中進(jìn)行攪拌腔流體域與攪拌葉片旋轉(zhuǎn)域的獨(dú)立網(wǎng)格劃分。生成四面體網(wǎng)格劃分，得出結(jié)果流體域劃分網(wǎng)格單元為405519個(gè)；葉片旋轉(zhuǎn)域劃分網(wǎng)格單元為85270個(gè)，網(wǎng)格單元?jiǎng)澐执笮?0mm。對(duì)攪拌葉片、刮刀與攪拌盤連接處進(jìn)行sizing控制細(xì)化網(wǎng)格，得到較為精確結(jié)果。具體網(wǎng)格劃分情況如圖2和圖3所示。

使用CFX進(jìn)行流體仿真模擬時(shí)，首先對(duì)流體和固體基本材料屬性進(jìn)行設(shè)置，本文中選用脫硫石膏、粉煤灰與少量水混合成流體[7]，該混合流體密度約為2000kg/m3，動(dòng)力粘度為3pa/s-1。本文忽略溫度對(duì)結(jié)果的影響，不考慮比熱容與熱傳遞系數(shù)。攪拌葉片選用不銹鋼，其密度為7850kg/m3，彈性模量2.0E+05，泊松比為0.3，材料屈服極限為250Mpa。

圖2 流體域網(wǎng)格劃分

圖3 旋轉(zhuǎn)域葉片網(wǎng)格劃分

同時(shí)定義攪拌過程出口和入口邊界條件。本文中設(shè)置進(jìn)口邊界條件為入口速度2m/s，出口邊界條件為靜壓力0Pa，初步設(shè)定攪拌葉片工作轉(zhuǎn)速為400r/min，攪拌機(jī)腔體壁面設(shè)定為無滑移邊界，考慮到攪拌出口處無回流要求即對(duì)出口進(jìn)行一定長度的外延。設(shè)定邊界如圖4所示。

圖4 流體邊界設(shè)置

圖5 400r/min攪拌葉片流場壓力分布

根據(jù)流場壓力云圖可以看出最大壓力分布在第一層刮刀外側(cè)，最大值約為0.16Mpa。原因?yàn)榇颂幮D(zhuǎn)半徑最大，所以切向速度相較于其他刮刀和葉片也最大。造成流體流速較快反饋到刮刀表面引起壓力增大。

3.2 流固耦合場分析

求解出上文中整體攪拌機(jī)葉片的表面壓力值，流體壓力作用在攪拌葉片上造成葉片形變，同時(shí)葉片的形變也會(huì)引起流場的變化，但在本文中葉片的形變較小，對(duì)流場的影響可以忽略，故只需要考慮單向耦合過程。本文中引起葉片的應(yīng)力變化的載荷主要有兩種，分別為流場中流體作用產(chǎn)生的壓力載荷和離心攪拌機(jī)旋轉(zhuǎn)產(chǎn)生的離心載荷。因此將以上兩種載荷綜合考慮到結(jié)構(gòu)靜力學(xué)分析中進(jìn)行求解，同時(shí)在靜力學(xué)模塊中根據(jù)攪拌葉片的放置方向設(shè)定電機(jī)軸逆時(shí)針旋轉(zhuǎn)轉(zhuǎn)速為400r/min。最終求得攪拌機(jī)葉片整體的最大形變量和最大等效壓力云圖如圖6、圖7所示。

圖6 400r/min葉片形變量分布

圖7 400r/min葉片等效應(yīng)力分布

根據(jù)模擬分析結(jié)果可知最大形變量約為0.0254mm，發(fā)生在攪拌葉片底端；最大等效應(yīng)力為6.0201Mpa，發(fā)生在攪拌葉片與攪拌盤相連接處內(nèi)側(cè)。這說明在攪拌過程中，攪拌葉片與流體主要接觸區(qū)域?yàn)閿嚢枨惑w底部靠邊側(cè)位置。從而引起最大形變位置在攪拌葉片底部，最大應(yīng)力集中位置靠近葉片與攪拌盤連接處內(nèi)側(cè)。

3.3 耦合結(jié)果分析

按照以上相同的分析方法和步驟，在CFX流體分析中對(duì)離心攪拌機(jī)轉(zhuǎn)速從100r/min至800r/min進(jìn)行模擬計(jì)算，將流場中得出的結(jié)果耦合到結(jié)構(gòu)分析中，最終得出最大形變量分布情況如圖8所示與最大等效應(yīng)力。通過圖示結(jié)果可以看出在攪拌機(jī)不同工況下攪拌葉片的最大形變量相對(duì)較小，極限轉(zhuǎn)速800r/min時(shí)最大形變量不超過0.1mm。這也就是在選擇流固分析類型時(shí)選擇單向流固耦合分析的原因，攪拌葉片的形變所引起的流場變化對(duì)模擬計(jì)算的結(jié)果影響較小，可以忽略。

圖8 不同轉(zhuǎn)速下葉片最大形變量分布圖

得到流固耦合條件下的等效應(yīng)力后，將在流場中得到的葉片表面壓力分布抑制，即單獨(dú)在結(jié)構(gòu)分析環(huán)境中求出耦合前僅受離心載荷下葉片的受力情況。得到最大等效應(yīng)力的位置與耦合后一致。耦合前后對(duì)比分析如圖9所示。從圖中可以得出，隨著轉(zhuǎn)速的增加，耦合前后的等效應(yīng)力差值呈現(xiàn)遞增狀態(tài)。100r/min至300r/min時(shí)差值變化較小，但是低速攪拌狀態(tài)下攪拌效率較低；轉(zhuǎn)速到500r/min以上時(shí)等效應(yīng)力差值較大，流體壓力在關(guān)鍵部位的影響已經(jīng)達(dá)到總應(yīng)力約18%已不可忽略。由此可見為了增加攪拌葉片壽命，減少流體對(duì)葉片的磨損情況，可以確定攪拌機(jī)的工作轉(zhuǎn)速規(guī)定在300r/min至500r/min時(shí)較為合適，攪拌效率和使用壽命都能得到提高。

圖9 耦合前后應(yīng)力對(duì)比圖

4 結(jié)語

（1）建立以一種用于生產(chǎn)環(huán)保墻板材料的離心攪拌機(jī)模型。為了驗(yàn)證其在工作時(shí)的葉片的應(yīng)力分布情況，使用了Workbench環(huán)境中流體分析工具CFX和結(jié)構(gòu)靜力學(xué)分析工具Static-Structural二者進(jìn)行單向流固耦合分析。在CFX中通過添加流體屬性、定義邊界條件、求解設(shè)定等得出流場中葉片表面壓力分布，將數(shù)值導(dǎo)入到結(jié)構(gòu)分析進(jìn)行求解得到最終結(jié)果。

（2）以攪拌軸轉(zhuǎn)速為400r/min為例，模擬得出攪拌機(jī)葉片在流場中的壓力分布數(shù)值，以及耦合到固體場中的最大形變和最大等效應(yīng)力。同時(shí)使用相同的分析方法對(duì)比100r/min至800r/min時(shí)模擬數(shù)據(jù)，耦合后最大形變量為0.0957mm，最大等效應(yīng)力為23.894Mpa，確定了流固耦合分析的類型，以及最佳的離心攪拌機(jī)工作轉(zhuǎn)速。通過分析最大形變、等效應(yīng)力等數(shù)據(jù)，此流固耦合模型可以為同類型攪拌葉片的疲勞分析提供依據(jù)，同時(shí)為工業(yè)生產(chǎn)提供參考建議。

[1] 張衎, 梁尚明. 基于Workbench整體葉輪流固耦合強(qiáng)度分析[J]. 機(jī)械, 2014, 41(04):46-49.

[2] 蔡琳, 江寒冰. 攪拌機(jī)葉片的流固耦合計(jì)算機(jī)仿真分析[J]. 制造業(yè)自動(dòng)化, 2012, 34(02):69-71.

[3] 趙悟, 岳超, 沈威威. 雙立軸振動(dòng)攪拌機(jī)攪拌機(jī)構(gòu)有限元分析與結(jié)構(gòu)優(yōu)化[J]. 機(jī)械強(qiáng)度, 2020, 42(06): 1496- 1502.

[4] 黎義斌, 梁開一, 李正貴. 基于流固耦合的斜軸式攪拌器水力性能數(shù)值分析[J]. 過程工程學(xué)報(bào), 2020, 20(12): 1424-1431.

[5] 高耀東, 胡鑫. 基于Workbench混凝土攪拌車葉片應(yīng)力分析[J]. 機(jī)械強(qiáng)度, 2016, 38(04): 887-891.

[6] 熊雄, 朱目成, 黃尚坤. 基于正交試驗(yàn)的槳式攪拌器流固耦合分析[J]. 化工機(jī)械, 2018, 45(06): 752-756.

[7] 李紅軍, 周嘯, 李旭, 等. 基于離散元素法的離心攪拌機(jī)混合特性研究[J]. 武漢紡織大學(xué)學(xué)報(bào), 2021, 34(04): 23-26.

Analysis of Fluid-solid Coupling Stress Field of Centrifugal Mixer Blade

LI Hong-jun1,2,YANG Kang1, ZHANG Chi1, CHEN Wei1,2

(1.College of Mechanical Engineering and Automation,Wuhan Textile University, Wuhan Hubei 430200, China; 2.Hubei Engineering Research Center for Intelligent Assembly of Industrial Detonators, Wuhan Hubei, 430200, China)

In order to study the stress distribution of centrifugal mixer blades during operation, this paper establishes a centrifugal mixer blade fluid-solid coupling stress field analysis model, using ANSYS-CFX The fluid analysis software couples fluid analysis and structural analysis, analyzes the stress field of the centrifugal mixer blade under high-speed conditions, and obtains the equivalent stress distribution diagram generated by the fluid load inside the blade and the maximum value after the overall blade is coupled. Deformation and working stress field. Finally, the appropriate working speed of the centrifugal mixer was obtained by analyzing the results of the stress field, which provided reference suggestions for industrial production.

mixer blade; CFX; Fluid-solid coupling; stress analysis

李紅軍（1973-），男，教授，研究方向：工業(yè)自動(dòng)化與圖形學(xué).

TP391.1

A

2095－414X(2022)01－0026－04