多載荷作用下動(dòng)渦旋盤應(yīng)力和變形研究

多載荷作用下動(dòng)渦旋盤應(yīng)力和變形研究*

杜萌1， 張賢明1**， 王立存2,3

(1.重慶工商大學(xué) 廢油資源化技術(shù)與裝備教育部工程研究中心，重慶 400067；2.制造裝備機(jī)構(gòu)設(shè)計(jì)與控制重慶市

重點(diǎn)實(shí)驗(yàn)室，重慶 400067；3.裝備系統(tǒng)服役健康保障重慶市國際聯(lián)合研究中心， 重慶 400067)

摘要：利用NX建立動(dòng)渦旋盤的三維實(shí)體模型，運(yùn)用有限元法研究了動(dòng)渦旋盤在氣體載荷、溫度載荷下的應(yīng)力與變形;并對氣體載荷以及溫度載荷耦合作用下的動(dòng)渦盤應(yīng)力和變形進(jìn)行了研究;仿真結(jié)果表明,溫度載荷作用下，最大應(yīng)力和變形發(fā)生在渦旋齒齒頭部位；氣體載荷作用下，最大變形在嚙合點(diǎn)處，最大應(yīng)力在渦旋齒中心及尾部;耦合作用下在端板與渦旋齒連接處產(chǎn)生最大應(yīng)力，在渦旋齒中心及尾部位產(chǎn)生了最大變形。

關(guān)鍵詞：真空泵；渦旋盤；載荷；耦合；應(yīng)力與變形

doi:10.16055/j.issn.1672-058X.2015.0011.023

收稿日期：2015-06-21；修回日期：2015-07-25.

基金項(xiàng)目：*國家自然科學(xué)

作者簡介：杜萌(1990-)，男，碩士研究生，從事環(huán)保關(guān)鍵技術(shù)與設(shè)備研究.E-mail:1010656355@qq.com.

通訊作者：**張賢明(1955-)，男，重慶市人，教授，博士后導(dǎo)師，從事真空濾油技術(shù)與裝備研究.

中圖分類號：O175.1文獻(xiàn)標(biāo)志碼：A

在渦旋真空泵工作過程中，動(dòng)、靜渦旋盤受到載荷的影響極易產(chǎn)生變形，變形會(huì)造成軸向間隙的增加，同時(shí)也會(huì)影響徑向間隙的密封和徑向嚙合效果，造成氣體泄漏增加，摩擦損耗增加等[1]，因此，對渦旋盤的應(yīng)力和變形研究是必不可少的。由于氣體被壓縮以及渦旋盤之間摩擦產(chǎn)生熱載荷，渦旋盤發(fā)生熱膨脹，進(jìn)而產(chǎn)生了熱應(yīng)力和熱變形。另外渦旋真空泵工作時(shí)，相鄰壓縮腔內(nèi)的氣體壓力差產(chǎn)生的力載荷同樣會(huì)使渦旋盤產(chǎn)生應(yīng)力與變形。目前，大多數(shù)文獻(xiàn)對渦旋盤應(yīng)力和變形研究基于單獨(dú)載荷作用下，并簡化渦旋盤的溫度場為均勻溫度場[2-6]。此處結(jié)合有限元分析軟件對渦旋盤耦合作用下的應(yīng)力和變形進(jìn)行研究，對渦旋盤的加工具有一定的意義。

1動(dòng)渦旋盤載荷

渦旋盤會(huì)受到外界各種力和荷載的作用從而發(fā)生應(yīng)力和變形，這其中以熱載荷和氣體載荷兩種影響最大[7]。

1.1溫度載荷

根據(jù)現(xiàn)有相關(guān)文獻(xiàn)研究成果，知道渦旋真空泵穩(wěn)定工況時(shí)，在一定的誤差范圍內(nèi)，動(dòng)、靜渦旋盤的溫度場為穩(wěn)態(tài)溫度場，即渦旋盤的溫度分布基本不隨時(shí)間發(fā)生變化。利用Workbench環(huán)境中table型表格數(shù)組進(jìn)行加載,溫度載荷沿端板和渦旋齒徑向線性衰減，軸承處沿軸向變化[8]。

1.2氣體載荷

渦旋真空泵的壓縮腔是隨著動(dòng)盤在不斷變化的，壓縮腔容積變化會(huì)影響壓縮腔內(nèi)的壓力分布。渦旋盤受到氣體壓力載荷的作用是渦旋齒內(nèi)、外壁面的壓力差所致。參考現(xiàn)有文獻(xiàn)的研究情況，此處將該問題適當(dāng)簡化，當(dāng)渦旋真空泵開始排氣，渦旋盤內(nèi)外壁表面產(chǎn)生最大壓差，此時(shí)圓盤和渦旋齒受力情況最為惡劣[9-11]。

壓縮腔內(nèi)體積變化的公式如下：

(1)

式(1)中，Vi為各壓縮腔體積編號；p為渦旋節(jié)距；t為渦旋齒厚；h為渦旋齒高。

假設(shè)渦旋真空泵壓縮過程是絕熱進(jìn)行的，則第i個(gè)壓縮腔對應(yīng)于主軸轉(zhuǎn)角的氣體壓力為

(2)

式(2)中，ps為吸氣壓力；k為氣體等熵指數(shù)；Vs為吸氣容積；Vi為第i個(gè)壓縮腔的容積。

2有限元仿真前處理

2.1渦旋盤模型建立

此次仿真通過NX8.5軟件建立渦旋真空泵主要部件的三維仿真模型，模型的徑向?yàn)閤y平面，齒高方向?yàn)閆方向，基圓半徑Rg=2.748 4 mm，Ror=4 mm,圈數(shù)N=3.25 mm，節(jié)距P=17.26 mm，齒厚t=4.71 mm，齒高h(yuǎn)=26.15 mm。

2.2單元網(wǎng)格劃分

渦旋盤材料目前大部分選擇合金或者鑄鐵等(圖1)，此次仿真采用的渦旋盤材料為鑄鐵。該材料的質(zhì)量密度為7.15 kg/cm3，膨脹系數(shù)13.3 E-6/℃，導(dǎo)熱系數(shù)0.052 6 W/mm·K彈性模量為113 Gpa，泊松比為0.26，其他屬性均在材料庫中給出。此次仿真將最小單元尺寸設(shè)置為1.6mm，選擇十節(jié)點(diǎn)四面體單元法對渦旋盤進(jìn)行自動(dòng)網(wǎng)格劃分，最后確定動(dòng)渦旋盤有限元網(wǎng)格中的單元數(shù)為7 588，節(jié)點(diǎn)數(shù)為51 260[12]。

2.3約束條件

端板周圍z方向自由度為零；動(dòng)渦旋軸承孔內(nèi)壁面xy方向的自由度為零；動(dòng)渦旋軸承孔頂部z向自由度為零。

圖1　動(dòng)渦旋盤網(wǎng)格模型

圖2　動(dòng)渦旋盤溫度場分布

3仿真結(jié)果與分析

3.1溫度載荷下的應(yīng)力和變形分析

結(jié)合文獻(xiàn)的計(jì)算結(jié)果，設(shè)計(jì)排氣口溫度為108.18 ℃，進(jìn)氣口溫度為22.2 ℃。選定渦旋盤徑向上的關(guān)鍵點(diǎn)來設(shè)置渦旋盤相應(yīng)面上的線性變化溫度載荷， 將各個(gè)腔計(jì)算的溫度依次施加到動(dòng)渦盤上面。設(shè)置渦旋盤的參考溫度為25 ℃，對施加于渦旋盤的溫度載荷進(jìn)行分析求解如圖2所示。圖3為動(dòng)渦盤在熱載荷作用下的應(yīng)力和變形圖，熱載荷應(yīng)力分布大致從渦旋齒中心向兩側(cè)逐漸降低，渦旋盤的最大等效應(yīng)力位于底座與端板連接位置也就是排氣口處，等效應(yīng)力值大小為225 MPa。渦旋盤變形最大部位在渦旋盤的渦旋齒頂面，變形最大值為16.36 μm，從渦旋齒中心向兩側(cè)遞減。

3.2氣體載荷下的應(yīng)力和變形分析

根據(jù)式(1)(2)計(jì)算結(jié)果，設(shè)定吸氣腔①的壓強(qiáng)為0.631 MPa，第一壓縮腔②的壓強(qiáng)為1.127 MPa，第二壓縮腔③的壓強(qiáng)為1.668 MPa，中心腔④的壓強(qiáng)為2.342 MPa。氣體載荷作用下的應(yīng)力和變形見圖4。

圖3　動(dòng)渦旋盤熱荷載應(yīng)力變形圖

圖4　不同壓縮腔定義

如圖5所示，渦旋真空泵工作時(shí)，處于渦旋盤外側(cè)的壓力為吸氣壓力，中心腔內(nèi)的壓力為排氣壓力。在渦旋齒上，氣體徑向力的作用是由于相鄰壓縮腔之間的壓差所產(chǎn)生，大小從中心的排氣壓力到外側(cè)的吸氣壓力逐漸變化。應(yīng)力最大值發(fā)生在渦旋齒中心，應(yīng)力值大小為83.02 MPa，并且沿渦旋齒向外等效應(yīng)力逐漸減小。動(dòng)渦旋盤氣體載荷變形最大值為9.449 μm。渦旋齒中心的變形比較大為4.04 μm，并且沿渦旋齒向外等效應(yīng)力逐漸減小，但在渦旋齒尾部變形突然變大為9.449 μm。這是由于其壁厚相對于整個(gè)渦旋齒來說較薄，載荷工況類似懸臂梁，因此應(yīng)力產(chǎn)生集中，變形也相應(yīng)變大。

圖5　動(dòng)渦旋盤氣體載荷的應(yīng)力變形圖

3.3動(dòng)渦旋盤多物理場耦合分析

圖6為渦旋盤在氣體載荷與熱載荷耦合作用下的應(yīng)力變形圖，應(yīng)力最大值發(fā)生在渦旋齒中心，為315 MPa，并且沿渦旋齒向外等效應(yīng)力逐漸減小。動(dòng)渦旋盤變形最大值為31.438 μm，位于渦旋齒中心。渦旋齒尾部以及渦旋盤底面處的變形比較小，渦旋盤變形沿渦旋齒向外逐漸遞減。

圖6　動(dòng)渦盤多物理場耦合應(yīng)力變形圖

4結(jié)論

溫度載荷從渦旋齒中心沿渦旋盤向兩邊遞減，熱變形產(chǎn)生的動(dòng)渦旋盤最大變形為16.36 μm。最大等效應(yīng)力為261.26 MPa，發(fā)生在渦旋齒中心。渦旋盤相鄰壓縮腔的渦旋齒壁面兩側(cè)受到的氣體壓力不同，故渦旋齒在兩側(cè)嚙合處存在波動(dòng)差，壓力差越大，波動(dòng)越大，在渦旋齒與端板連接的部分最為明顯。動(dòng)盤等效應(yīng)力最大值為83.02 MPa。渦旋齒中心的齒頭部位發(fā)生最大變形，并且沿齒頂向齒尾變形逐漸減小。動(dòng)盤單元節(jié)點(diǎn)變形最大幅值為6.107 μm。多場耦合作用下等效應(yīng)力值和變形不是線性疊加的，溫度載荷對應(yīng)力及變形趨勢影響較大，渦旋盤的最大應(yīng)力為315.2 MPa，發(fā)生在渦旋齒中心處。經(jīng)粗略計(jì)算可知，鑄鐵的強(qiáng)度可以承受動(dòng)渦旋盤的最大主應(yīng)力，即渦旋盤結(jié)構(gòu)強(qiáng)度是安全的，最大變形為31.438 μm。

參考文獻(xiàn)：

[1] 劉興旺，趙嫚，李超，等. 渦旋壓縮機(jī)的徑向迷宮密封研究[J]. 機(jī)械工程學(xué)報(bào)，2012，48(21)：97-104

[2] 楊廣衍，盛林，張秀麗，等. 無油渦旋真空泵溫度場、應(yīng)力場和熱變形的ANSYS模擬[J]. 真空，2008，45(5):17-19

[3] 楊廣衍，張先鋒，趙元. 渦旋真空泵應(yīng)力與變形的ANSYS模擬與研究[J]. 真空， 2007，44(1):12-13

[4] 殷俊，楊美傳，馮鑒.基于熱應(yīng)力場耦合的渦旋壓縮機(jī)動(dòng)渦旋盤有限元分析[J].壓縮機(jī)技術(shù)，2011(6)：6-9

[5] 金丹，陳旭，趙樹峰.動(dòng)渦盤在氣體力和溫度載荷作用下的應(yīng)力及變形分析[J].壓縮機(jī)技術(shù)，2004(5)：7-10

[6] 張立群，劉振全.渦旋壓縮機(jī)靜渦旋盤實(shí)際工況下的變形分析[J].流體機(jī)械，2000，28(2)：18-20

[7] 李超，謝文君，趙嫚.多場耦合作用下動(dòng)渦旋盤的變形和應(yīng)力研究[J].流體機(jī)械， 2013，41(8):16-20

[8] 黃蕾,唐景春,韓坤.溫度場對動(dòng)渦旋盤應(yīng)力與應(yīng)變的影響[J].低溫與超導(dǎo),2013,41(05):60-63

[9] QIANG J，PENG B，LIU Z. Dynamic Model for the Orbiting Scroll Based on The Pressures in Scroll Chambers — Part II: Investigations on Scroll Compressors snd Model Validation[J]. International Journal of Refrigeration，2013，36:1850-1865

[10] LIU Y G， HUNG C H， CHANG Y C. Study on Involute of Circle with Variable Radii in a Scroll Compressor [J]. Mechanism and Machine Theory，2010，45(11):1520-1536

[11] HUANG H，LI Q，YUAN D，et al. An experimental study on variable air volume operation of ducted air-conditioning with digital scroll compressor and conventional scroll compressor[J]. Applied Thermal Engineering，2008，28(7):761-766

[12] 曹華軍，舒林森，許磊，等.復(fù)雜機(jī)械零件的六面體有限元網(wǎng)格生成方法[J].機(jī)械工程學(xué)報(bào)，2014，50(15)：113-118

[13] 王金梅，劉乘.應(yīng)力-能量法求取泡沫塑料緩沖曲線時(shí)函數(shù)模型的研究[J].包裝工程工程版，2014，35(5)：79-82；101

Research on Stress and Deformation of Dynamic VortexDisk under Multi-load

DU Meng1， ZHANG Xian-ming1*， WANG Li-cun2，3

(1. Engineering Research Center for Waste Oil Recovery Technology and Equipment of Ministry of Education

of China,Chongqing Technology and Business University,Chongqing 400067,China; 2. Chongqing Key Lab

of Manufacturing Equipment Mechanism Design and Control,Chongqing Technology and Business University,

Chongqing 400067,China; 3.Chongqing Research Center of System HealthMaintenance,Chongqing Technology

and Business University,Chongqing 400067,China)

Abstract：The 3D models of dynamic scroll disks in NX software are set up,the stress and deformation of dynamic vortex disks under the condition of the given gas load and constraint conditions are studied by the finite element analysis method,and the stress and deformation of scroll disk under coupled action of gas load and temperature load is also studied. Simulation results showed that,under the action of temperature load,maximum stress and deformation occurred at the tooth head of the vortex teeth,that the maximum deformation occurred at meshing points and that the maximum stress occurred at the center and the tail of the scroll teeth. The maximum stress was generated at the connections between end plate and scroll teeth while the maximum deformation occurred at the center and the tail of the vortex teeth.

Key words： vacuum pump; scroll disk; load; coupling; stress and deformation