半軸殼體左右兩面孔加工組合機(jī)床總體設(shè)計
喜歡這套資料就充值下載吧。資源目錄里展示的都可在線預(yù)覽哦。下載后都有,請放心下載,文件全都包含在內(nèi),有疑問咨詢QQ:414951605 或 1304139763
洛陽理工學(xué)院
畢業(yè)設(shè)計(論文)開題報告
院(部): 機(jī)械工程 年 月 日 (學(xué)生填表)
課題名稱
半軸殼體左右兩面孔加工組合機(jī)床的總體設(shè)計
學(xué)生姓名
聶俊超
專業(yè)班級
B110207
課題類型
工程設(shè)計
指導(dǎo)教師
鮑莉
職稱
課題來源
1. 綜述本課題國內(nèi)外研究動態(tài),說明選題的依據(jù)和意義
隨著社會的不斷進(jìn)步,機(jī)械加工技術(shù)的不斷發(fā)展,傳統(tǒng)的機(jī)床已不能完全適應(yīng)新形勢的要求。傳統(tǒng)的機(jī)床只能對一種零件進(jìn)行單刀,單工位,單軸,單面加工,生產(chǎn)效率低且加工精度不穩(wěn)定,為了克服傳統(tǒng)機(jī)床的弊端,工程技術(shù)人員相應(yīng)地設(shè)計出了專用機(jī)床。但由于專用機(jī)床是根據(jù)某一工藝要求專門設(shè)計制造的,且它的組成部件均是專門設(shè)計制造的,因此相對于傳統(tǒng)機(jī)床而言,專用機(jī)床的造價過于昂貴,設(shè)計制造周期長。為了解決傳統(tǒng)機(jī)床與專用機(jī)床之間的矛盾組合機(jī)床便應(yīng)運(yùn)而生了,組合機(jī)床兼有低成本和高效率的優(yōu)點,在大批、大量生產(chǎn)中得到廣泛應(yīng)用,在組合機(jī)床上可以完成鉆孔、擴(kuò)孔、銑削、磨削等工序,生產(chǎn)效率高,加工精度穩(wěn)定,引起了越來越多工程人員的關(guān)注。本課題針對半軸殼體鉆孔組合 機(jī)床設(shè)計,有利于提高大批量生產(chǎn)的生產(chǎn)效率,提高加工精度穩(wěn)定性,節(jié)約各的
資源。
最早的組合機(jī)床于1911年在美國制成,用于加工汽車零件之后便廣泛應(yīng)用于大批量生產(chǎn)的機(jī)械工業(yè)中,并且隨著機(jī)械工業(yè)的發(fā)展而逐步完善。我國的組合機(jī)床的發(fā)展已有31年的歷史,其科研和生產(chǎn)都具有相當(dāng)?shù)幕A(chǔ),應(yīng)用也深入到很多行業(yè),它是提高生產(chǎn)效率和實現(xiàn)高速發(fā)展必不可少的設(shè)備之一。組合機(jī)床及其自動線是集機(jī)電于一體的綜合自動化程度較高的制造技術(shù)和成套工藝裝備。它的特征是高效、高質(zhì)、經(jīng)濟(jì)實用,因而被廣泛應(yīng)用于工程機(jī)械、交通、能源、軍工、輕工、家電等行業(yè)。我國傳統(tǒng)的組合機(jī)床及組合機(jī)床自動線主要采用機(jī)、電、氣、液壓控制,它的加工對象主要是生產(chǎn)批量比較大的大中型箱體類和軸類零件(近年研制的組合機(jī)床加工連桿、板件等也占一定份額),完成鉆孔、擴(kuò)孔、鉸孔,加工各種螺紋、鏜孔、車端面和凸臺,在孔內(nèi)鏜各種形狀槽,以及銑削平面和成形面等。隨著技術(shù)的不斷進(jìn)步,一種新型的組合機(jī)床——柔性組合機(jī)床越來越受到人們的青睞,它應(yīng)用多位主軸箱、可換主軸箱、編碼隨行夾具和刀具的自動更換,配以可編程序控制器(PLC)、數(shù)字控制(NC)等,能任意改變工作循環(huán)控制和驅(qū)動系統(tǒng),并能靈活適應(yīng)多品種加工的可調(diào)可變的組合機(jī)床。另外,近年來組合機(jī)床加工中心、數(shù)控組合機(jī)床、機(jī)床輔機(jī)(清洗機(jī)、裝配機(jī)、綜合測量機(jī)、試驗機(jī)、輸送線)等在組合機(jī)床行業(yè)中所占份額也越來越大。
縱觀世界四大博覽會的概況我國組合機(jī)床裝備的發(fā)展程度與先進(jìn)技術(shù)還有相當(dāng)大的差距,因此我國組合機(jī)床需要繼續(xù)發(fā)展和提高組合機(jī)床自動線水平,
因此我國組合機(jī)床需要繼續(xù)發(fā)展和提高組合機(jī)床自動線水平,這是我國當(dāng)前機(jī)械加工行業(yè)的一項艱巨任務(wù)。在發(fā)展的同時還需進(jìn)一步發(fā)展不同程度的低,中,高檔機(jī)機(jī)電一體化類產(chǎn)品,這有助于滿足不同生產(chǎn)規(guī)模,大中小不同技術(shù)水平的企業(yè)的要求。這對我國加工制造業(yè)提高加工精度,加工效率,自動化程度和綜合制造技術(shù)水平都有著非常深遠(yuǎn)和有利的長遠(yuǎn)趨勢。
2. 研究的基本內(nèi)容,擬解決的主要問題
本課題的基本內(nèi)容是設(shè)計載重汽車主傳動軸萬向節(jié)叉端面鉆孔組合機(jī)床。?
基本內(nèi)容:
?(1)零件分析?
仔細(xì)了解所欲加工零件的加工特點,精度和技術(shù)要求。
?(2)擬定零件工藝方案?
確定零件在組合機(jī)床上完成加工時的定位基準(zhǔn)以及零件的加工工藝方案。?
(3)總體方案設(shè)計?
根據(jù)工藝方案確定機(jī)床的配置形式及總體結(jié)構(gòu)布局。?(4)組合機(jī)床設(shè)計——三圖一卡設(shè)計。
主要解決方法是通過已學(xué)過的基礎(chǔ)理論和專業(yè)知識并參考各類相關(guān)專業(yè)書籍,熟悉組合機(jī)床工作原理和設(shè)計原則,了解各項設(shè)計參數(shù)和技術(shù)要求進(jìn)行設(shè)計計算,同時利用繪圖軟件進(jìn)行初步結(jié)構(gòu)設(shè)計,參考已有的先進(jìn)設(shè)計實例,學(xué)習(xí)其計算過程和其中的先進(jìn)創(chuàng)新理念并運(yùn)用到自己的計算創(chuàng)新設(shè)計中。相關(guān)自已無力解決的問題可以與同學(xué)討論解決,一些自己的創(chuàng)新思想也可與同學(xué)討論其可行性,或者向指導(dǎo)教師請教。
3. 研究步驟、方法及措施
計劃階段:提出任務(wù)—分析零件加工對機(jī)器的需求—確定任務(wù)要求。?
方案設(shè)計:機(jī)床功能分析—提出幾種可能的解決方案—組合幾種可能的方案—對比評價—討論優(yōu)化—得出最佳方案。?
技術(shù)設(shè)計:明確構(gòu)形要求—結(jié)構(gòu)化—選擇材料、計算尺寸—評價—決策確定結(jié)構(gòu)形狀及尺寸—零件設(shè)計—部件設(shè)計—總體設(shè)計。
可行性分析:?
(1)對組合機(jī)床的設(shè)計有一定的掌握。?
(2).運(yùn)用機(jī)械制造技術(shù)知識可以完成“三圖一卡”。
(3)運(yùn)用機(jī)械制造裝備設(shè)計可以解決主軸箱設(shè)計的過程??。
(4)運(yùn)用機(jī)械設(shè)計的知識可以完成傳動零件的校核??。
(5)運(yùn)用機(jī)械制圖的知識和計算機(jī)繪圖軟件的幫助可以完成設(shè)計圖紙的繪制?。
(6)已經(jīng)擁有大量有關(guān)組合機(jī)床的書籍、圖樣,可以到圖書館查閱相關(guān)國家標(biāo)準(zhǔn)?。
(7)參考相關(guān)的設(shè)計的經(jīng)驗.。
(8)與指導(dǎo)老師討論協(xié)商。?
自己已掌握的機(jī)械設(shè)計和組合機(jī)床的相關(guān)專業(yè)知識和多次課程設(shè)計的經(jīng)驗,再結(jié)合與同學(xué)的討論及指導(dǎo)老師的幫助可以完成本課題的設(shè)計。
4. 研究工作進(jìn)度
第1~3周 節(jié)約相關(guān)資料、了解相關(guān)知識,畢業(yè)實習(xí),完成開題報告。
第4~6周 熟悉加工零件,了解用途,繪制零件圖。
第7~13周 確定設(shè)計方案,分析,計算,繪圖。
第14~15周 編寫說明書,定稿,打印,答辯前準(zhǔn)備階段
5. 主要參考文獻(xiàn)
【1】楊可楨、陳光藻主編《機(jī)械設(shè)計基礎(chǔ)》,人民出版社,1990年。
【2】余振華主編,《組合機(jī)床及其調(diào)整使用》機(jī)械工業(yè)出版社,1990年4月。
【3】哈爾濱工業(yè)學(xué)編,《機(jī)床制造工藝設(shè)計手冊》,哈爾濱工業(yè)大學(xué)出版社,1981年5月。
【4】《機(jī)床設(shè)計手冊》編寫組編,《機(jī)床設(shè)計手冊》,化學(xué)工業(yè)出版社,1986年12月
【5】姜作敬主編,《機(jī)械制造工藝學(xué)》,華中理工大學(xué)出版社,1989年6月。
【6】孫靖民主編,《機(jī)械優(yōu)化設(shè)計》,機(jī)械工業(yè)出版社,1990年5月。
【7】楊玉鼎主編,《機(jī)械設(shè)計》,重慶大學(xué)出版社,1989年9月
【8】沈樂年、范珍良等編,《機(jī)械設(shè)計基礎(chǔ)》,清華大學(xué)出版社,1989年9月。
教研室意見
教研室主任簽字: 年 月 日
機(jī)床生產(chǎn)率計算卡
被加工零件
圖號
毛胚種類
鑄 件
名稱
半軸殼體
毛胚重量
材料
HT200
硬 度
170~230HB
工序名稱
半軸殼體左右兩面孔加工
工序號
序
號
工步名稱
被加
工零
件數(shù)
加工
直徑
mm
加工
長度
mm
工作
行程
mm
切削
速度
m/min
轉(zhuǎn)速
r/min
每轉(zhuǎn)
進(jìn)給量
mm/r
每分鐘
進(jìn)給量
mm/min
工 時 min
機(jī)動
時間
輔助
時間
共計
1
裝入工件
1
0.25
2
工件定位夾緊
0.25
3
左動力部件
快進(jìn)
145
5000
0.029
鉆?23.4孔
23.4
35
55
14.69
200
0.375
75
0.582
鉆?17.4孔
17.4
34
55
16.39
300
0.25
75
0.585
快退
200
5000
0.02
4
右動力部件
快進(jìn)
145
10000
0.014
鉆?13孔
13
15
35
16.32
400
0.18
72
0.553
快退
180
10000
0.018
5
松開工件
0.002
6
卸下工件
0.5
備
注
本機(jī)床裝卸取1分
單位工時
1.915
機(jī)床實際生產(chǎn)率Q1
25.53件/h
機(jī)床理想生產(chǎn)率Q
31.32件/h
負(fù)荷率
0.81
被加工零件
圖號
毛胚種類
鑄 件
名稱
半軸殼體
毛胚重量
材料
HT200
硬 度
170~230HB
工序名稱
半軸殼體左右兩面孔加工
工序號
序
號
工步名稱
被加
工零
件數(shù)
加工
直徑
mm
加工
長度
mm
工作
行程
mm
切削
速度
m/min
轉(zhuǎn)速
r/min
每轉(zhuǎn)
進(jìn)給量
mm/r
每分鐘
進(jìn)給量
mm/min
工 時 min
機(jī)動
時間
輔助
時間
共計
1
裝入工件
1
0.25
2
工件定位夾緊
0.25
3
左動力部件
快進(jìn)
145
5000
0.029
鉆?23.4孔
23.4
35
55
14.69
200
0.375
75
0.582
鉆?17.4孔
17.4
34
55
16.39
300
0.25
75
0.585
快退
200
5000
0.02
4
右動力部件
快進(jìn)
145
10000
0.014
鉆?13孔
13
15
35
16.32
400
0.18
72
0.553
快退
180
10000
0.018
5
松開工件
0.002
6
卸下工件
0.5
備
注
本機(jī)床裝卸取1分
單位工時
1.915
機(jī)床實際生產(chǎn)率Q1
25.53件/h
機(jī)床理想生產(chǎn)率Q
31.32件/h
負(fù)荷率
0.81
洛陽理工學(xué)院
畢業(yè)設(shè)計(論文)中期教學(xué)檢查表
系(部): 專業(yè):
指導(dǎo)教師姓名: 檢查時間: 年 月 日 (系(部)填表)
項目
檢查情況
指導(dǎo)教師對設(shè)計任務(wù)準(zhǔn)備和布置情況
畢業(yè)設(shè)計(論文)周數(shù) ( )
畢業(yè)設(shè)計(論文)任務(wù)數(shù)下達(dá)時間 ( )
畢業(yè)設(shè)計(論文)起始時間 ( )
學(xué)生對設(shè)計任務(wù)清楚程度如何? ( )
設(shè)計前教師進(jìn)行的準(zhǔn)備環(huán)節(jié):
1、
2、
3、
教師指導(dǎo)學(xué)生
每周指導(dǎo)學(xué)生 ( )次
一周指導(dǎo)學(xué)生 ( )小時
輔導(dǎo)學(xué)生設(shè)計認(rèn)真程度 ( )
已講課的次數(shù) ( )次
教師對學(xué)生的要求
學(xué)生是否按進(jìn)度計劃執(zhí)行 ( )
目前進(jìn)展情況 ( )
參考文獻(xiàn)要求 ( )本
學(xué)生對指導(dǎo)教師的評價及意見:
其它建議及意見
對該指導(dǎo)教師的指導(dǎo)情況評價
好 ( )
較好 ( )
一般 ( )
差 ( )
檢查人員簽名
洛陽理工學(xué)院
畢業(yè)設(shè)計(論文)任務(wù)書
填表時間: 2015年 2 月 28 日 (指導(dǎo)教師填表)
學(xué)生姓名
聶俊超
專業(yè)班級
B110207
指導(dǎo)教師
鮑莉
課題 類型
工程設(shè)計
題目
半軸殼體左右兩面孔加工組合機(jī)床的總體設(shè)計
主要研究
目標(biāo)
(或研
究內(nèi)容)
1. 組合機(jī)床方案的擬定;
2. 選擇刀具、確定相應(yīng)的切削用量;
3. 繪制被加工零件工序圖、加工示意圖、機(jī)床聯(lián)系尺寸圖;
4. 編制機(jī)床生產(chǎn)率計算卡
課題要求、主要任務(wù)及數(shù)量(指圖紙規(guī)格、張數(shù),說明書頁數(shù)、論文字?jǐn)?shù)等)
獨立完成不少于3張A0圖紙的結(jié)構(gòu)設(shè)計圖、裝配圖和零件圖,其中應(yīng)包含1張以上用計算機(jī)繪制的具有中等難度的A1號圖紙,同時至少有折合A1圖幅以上的圖紙用三維繪制,查閱15篇以上與題目相關(guān)的文獻(xiàn)。并需翻譯一萬印刷符(或譯出3000漢字)以上的有關(guān)技術(shù)資料(并附原文)。畢業(yè)論文的字?jǐn)?shù)一般為1..5萬字(約40~50頁)。
進(jìn)度計劃
第1~3周:借閱相關(guān)資料、了解相關(guān)知識,畢業(yè)實習(xí),完成開題報告
第4~6周:熟悉加工零件、了解用途、繪制零件圖
第7~13周:確定設(shè)計方案、分析、計算、繪圖
第14~15周:編寫說明書、定稿、打印、答辯前準(zhǔn)備階段
主要參
考文獻(xiàn)
[1] 沈陽工業(yè)大學(xué) 大連鐵道學(xué)院 吉林工學(xué)院編.組合機(jī)床設(shè)計.上海:上??茖W(xué)技術(shù)出版社,1985.
[2] 周誦明. 機(jī)械類專業(yè) 畢業(yè)設(shè)計指導(dǎo)書. 武漢:華中理工大學(xué)出社,1993.
[3] 趙雪松,任小中, 于華. 機(jī)械制造裝備設(shè)計.武漢:華中科技大學(xué)出版社, 2009.
[4] 謝家瀛. 組合機(jī)床設(shè)計簡明手冊.北京:機(jī)械工業(yè)出版社,2002.
指導(dǎo)教師簽字: 系主任簽字: 年 月 日
226 2012,24(2):226-234 DOI: 10.1016/S1001-6058(11)60238-2 RESEARCH OF INNER FLOW IN A DOUBLE BLADES PUMP BASED ON OPENFOAM* LIU Hou-lin, REN Yun, WANG Kai, WU Deng-hao, RU Wei-min,TAN Ming-gao Research Center of Fluid Machinery Engineering and Technology, Jiangsu University, Zhengjiang 212013, China, E-mail: (Received September 24, 2011, Revised January 9, 2012) Abstract: The inner flow analysis of centrifugal pumps has gradually become an important issue for the hydraulic design and performance improvement. Nowadays, CFD simulation toolbox of pump inner flow mainly contains commercial tools and open source tools. There are some defects for commercial CFD software for the numerical simulation of 3-D turbulent internal flow in pump, especially in capturing the flow characteristics under the off-design operating conditions. Additionally, it is difficult for researchers to do further investigation because of the undeclared source. Therefore, an open source software like Open Field Operation and Manipulation (OpenFOAM) is increasingly popular with researchers from all over the world. In this paper, a new computational study was implemented based on the original solver and was used to directly simulate the steady-state inner flow in a double blades pump, with the specific speed is 111. In order to disclose the characteristics deeply, three research schemes were conducted. The ratios (/dQ Q) of the flow rate are 0.8, 1.0 and 1.2, respectively. The simulation results were verified with the Particle Imaging Velocimetry (PIV) experimental results, and the numerical calculation results agree well with the experimental data. Meanwhile, the phenomena of flow separation under the off-design operating conditions are well captured by OpenFOAM. The results indicate that OpenFOAM possesses obvious strong predominance in computing the internal flow field of pump. The analysis results can also be used as the basis for the further research and the improvement of centrifugal pump. Key words: numerical simulation, double blades pump, internal flow, Particle Imaging Velocimetry (PIV) Introduction Double blades pumps is a kind of centrifugal pump with two blades. There are two symmetrical curve impeller passages from its inlet to outlet and the impeller outlet is quite wide. Therefore, it usually becomes one of impeller shapes in solid-liquid two- phase centrifugal pump. However, due to the short development history and the imperfect design theory, its impellers are often designed with a combination of experience of the designers so far in practice, so that its performance and stability are not ensured1-3. Determining the pump performance is decided by * Project supported by the National Natural Science Funds for Distinguished Young Scholar (Grant No. 50825902), the National Natural Science Foundation of China (Grant Nos. 51079062, 51179075 and 51109095) and the Natural Science Foundation of Jiangsu Province (Grant Nos. BK2009006, BK2010346). Biography: LIU Hou-lin (1971-), Male, Ph. D., Professor examining its inner flow characteristics is undoubte- dly the best method to improve the performance of pumps4-6. Recently, with the rapid progress in CFD and computer technology, the internal flow simulation has gradually become the important foundation of optimization and design for turbo-machinery7-10. Now, the fluid machinery CFD simulation toolbox mainly contains commercial tools and open source tools. Over the years, the commercial software pack- ages are fashionable in the world by its abundant fun- ction and fine easy-use quality. On the other hand, as a CFD software, commercial tool is not very profe- ssional, the computational results for pumps are less than satisfactory, especially in capturing the flow cha- racteristics under the off-design operating conditions. Furthermore, its undeclared source brings a considera- ble inconvenience for the application of numerical simulation in fluid machineries. Although codes can be added to implementation through user-define-fun- ctions, it has a strong limitation. For example, when the SIMPLE algorithm and two-equation turbulence 227models have to be improved, one need to have a through understanding of the governing equations, discretization method, turbulence models and iteration algorithm. However, the core algorithm code and data processing method can not be acquired because of commercialization and only several options can be chosen. Therefore, many open source CFD software are becoming popular and a high quality open source CFD simulation platform like the Open Field Operation and Manipulation (OpenFOAM) is out- standing due to its powerful function, clear archite- cture, expand feature and unified format. The OpenFOAM CFD toolbox was released as an open source December 10, 2004, which is based on C+ routine and contains many C+ modules which can be freely combined with some other modules such as tensor, vector, turbulence models, numerical algori- thm, discretion modules, automatic control modules and so on. Therefore, it is convenient to employ its solvers in simulating complex physical models in che- mical reaction, turbulent flow and heat conduction, etc. A variety of work on the internal flow in fluid machinery via OpenFOAM was done. Nilsson11 con- ducted the steady and unsteady computation of the flow in the Hllerforsen turbine runner and draft tube, and compared the results obtained by the OpenFOAM with those by the CFX-5 and in experiments. Eve- ntually, the applicability and reliability of the OpenFOAM in a Kaplan water turbine runner and draft tube was verified. Petit et al.12 validated an im- plementation of the General Grid Interface (GGI) in the OpenFOAM by using the frozen rotor steady app- roach and the sliding grid unsteady approach. How- ever, all the simulations were performed for a simpli- fied 2-D model of a centrifugal pump. Li13 simulated boundary layer in wind tunnel by OpenFOAM, which revealed that it was suitable for using OpenFOAM to conduct a Computational Wind Engineering (CWE) research. Currently, there is not much work regarding the study of the inner flow of pumps through compa- rison between the OpenFOAM simulation and experi- ments, and the related reports are seldom found. As an open source code, the OpenFOAM provi- des direct access to models and solver implementation details. However, there are some defects for the OpenFOAM in the numerical simulation of 3-D turbu- lent internal flow in hydroturbines. For the CFD simu- lation of hydroturbines, separate 3-D mesh passages or full geometry are generally connected together in order to simulate the flow of water through a succe- ssion of complex geometries like pumps where the stationary suction and volute are located, along with the rotating impeller. The requirement to fit all the meshes with conformal matching interfaces is often very difficult or leads to geometric compromises that would affect the numerical quality of the simulation results. Therefore, there is a need for a treatment of rotor-stator interfaces, which is necessary for the simulation of the whole hydroturbines inner flow. Although the capability has existed in the OpenFOAM, there is no definition in the case directory. There is also a need for a set of boundary conditions that makes it easy to capture basic features in a similar way as it can be done in some other CFD solvers. Furthermore, relaxation factors which control under- relaxation, have an important influence on improving stability of a computation. However, there is no gui- ding principle about those factors. Therefore, the paper focuses on centrifugal pump with the conside- ration of these factors. The numerical simulation of pump by using the OpenFOAM is achieved and the computational results are verified by Particle Imaging Velocimetry (PIV) experiments. The research also provides the foundation for achieving higher computa- tion accuracy of pump inner flow by improving the CFD method with self-compiling program in the OpenFOAM. In this paper, in order to compute the intera- ctions between rotating and fixed components in pumps, a Multiple Reference Frames (MRF) solver is used. At the same time, the simulation results are also validated by PIV test. The operating system used is SUSE Linux 10.3, and the version number used for the present computations is OpenFOAM 1.5-dev. 1. Numerical method and model 1.1 Governing equations The OpenFOAM toolbox already provides a solver called MRFSimpleFoam for solving the steady- state Reynolds-Averaged Navier-Stokes equations with turbulence models, such as the standard k model. The coupling between velocity and pressure is treated using the SIMPLE method14. The MRFSimpleFoam solver employs the finite-volume technique to discretize the Navier-Stokes equations in the rotating reference frame: =0Ru +2+=RRR uuur +Rp u where Ru is flow velocity in rotating frame, rotating frame speed, r position vector, p fluid pressure, fluid density, and the kinematic vis- cosity. 1.2 General Grid Interface (GGI) Due to the interaction between stator and rotor, 228 how to cope with the grids and information tran- smission of the coupling parts in the computational domain is a key issue to simulate precisely the flow fields15-17. The frozen rotor method in MRFSimpleFOAM solver is a steady-state formula- tion where the relative positions in rotor and stator are fixed. In the same time frame, the rotor and the stator parts will be meshed separately. For non-stationary turbo-machinery simulations, the relative rotation of mesh parts will necessarily produce non-conformal in- terfaces between the fixed and moving sections. A connection between these meshes is needed in order to simplify the mesh complexity in various turbo-machi- nery simulations and hence reduce the computer time cost. The GGI, developed by Beaudoin and Jasak18 can be used for that purpose in OpenFOAM. It is a new coupling interface for OpenFOAM, joining multi- ple non-conformal regions where the patches nodes on each side of the interface do not match. This interface uses weighted interpolation to eva- luate and transmit flow values across a pair of confor- mal or non-conformal coupled patches. The basic GGI interface is similar to a case of “static” sliding inter- faces with the advantage that no remeshing is required for the neighboring cells of the interface. The GGI uses the Sutherland-Hodgman algori- thm12 for computing the master and shadow face in- tersection area. Some quick rejection algorithms based on an axis aligned bounding box have been impleme- nted to speed up the search for potential face neigh- bors. Then, in order to rapidly handle the final non- overlapping filtering test, an efficient Hormann- Agathos point-in-polygon algorithm19 has been inclu- ded into the separating axis theorem algorithm18. Finally, discretization effects are taken into account in order to properly scale the GGI weighting factors to handle the possible presence of non-overlapping faces and hence keep the GGI interface conservative. 1.3 Boundary and initial conditions While the Partial Differential Equation (PDE) is solved with the finite volume method, a suitable inter- polation scheme of values typically from cell centres to face centres has a great effect on the numerical results, especially for the convection term. The conve- ction scheme of the existing solver is specified as default by limited linear differencing, which is a Total Variation Diminishing (TVD) scheme. Although it offers a second-order accurate discretization scheme for the convection term, it creates an unconditionally unstable discretisation practice20. In order to achieve stability, a fist-order accurate upwind differencing scheme has been introduced and the simulation results show that the TVD scheme can result more easily in iteration divergence and computation failure than the upwind differencing scheme in simulating the pump inner flow. Therefore, TVD scheme is not applicable to simulate the pump inner flow in the OpenFOAM. In this paper, Guassian up-wind scheme is used and can get the satisfactory results in numerical tests. In addition to appropriate discretization schemes, under-relaxation is another important technique used for improving the stability of a computation, particu- larly in solving steady-state problems. Under-relaxa- tion works by limiting the amount, in which a variable changes from one iteration to the next, either by limi- ting the solution matrix and source prior to solving for a field or by modifying the field directly. An under- relaxation factor , 01 specifies the amount of under-relaxation, ranging from none at all for =1 and increasing in strength as 0. There- fore, selecting an appropriate relaxation factor has great influence on the efficiency of computation. If the relaxation factor is too large, it will lead to divergence of computation easily. If it is too small, the result will converge slowly. Besides, an appropriate relaxation factor depends on the specific problem itself. Thus, there is no instruction about relaxation factor in simu- lating pump inner flow. At the same time, if the rela- xation factors are used by default in the OpenFOAM to examine the pump inner flow, the result will be un- stable. In this paper, the relaxation factors are suitable for simulating the pump inner flow according to +=1up20. At last, the relaxation factors are determined after many tests and the under-relaxation factors of relevant variables (i.e., pressure, momentum, turbulence kinetic energy and turbulence dissipation rate) are 0.3, 0.7, 0.3 and 0.3, respectively. In order to simulate the flow field in and the whole pump, the GGI method is needed to transmit the information between rotor and stator. 1.4 Model A 3-D model of double-blade pump for simula- tion is generated by Pro/E. The suction chamber is designed by semi-spiral method, while the volute is designed by equal velocity moment method and the cross section is rectangular, and type line is the loga- rithmic spiral. The design parameters of the double blades pump are shown in Table 1. Herein the calculation formulas of sn is 3/43.65=sn QnH Before the simulation, studying the grid indepe- ndency and selecting the turbulent model are nece- ssary21,22. The geometry is meshed in hybrid grid by GAMBIT and the OpenFOAM is used to simulate the inner flow in the double-blade pump. The data for stu- dying grid independency are shown in Table 2. If the head difference is smaller than 0.2%, the grid number is acceptable. According to the computation, Scheme 2 meets the need. So Scheme 2 of grid is adopted. The grids are show in Fig.1. 229 Table 1 The design parameters Parameter Sign Value Flow rate Q (m3/h) 25.86 Head H (m) 2.53 Rotation speed n (r/min) 750 Specific speed sn 111 Suction chamber inlet diametersD (m) 0.08 Impeller inlet diameter iD (m) 0.09 Blade inlet diameter 1D (m) 0.0812 Impeller outlet diameter 2D (m) 0.2 Impeller outlet width 2b(m) 0.047 Blade inlet angle 1(o) 18.3 Blade outlet angle 2 (o) 30 Volute base circle diameter 3D (m) 0.212 Volute inlet width 3b (m) 0.077 Table 2 Data for studying grid independency Grid number No. Impeller Volute Suction Total Head, H (m)1 618 738 220 175 300 103 1 139 016 2.46034 2 577 913 179 340 255 937 1 013 190 2.45974 3 505 999 143 450 202 830 852 279 2.43180 4 378 565 151 341 99 619 629 525 2.40151 Fig.1 Grids Table 3 Comparison of turbulence models Turbulent model Head, H (m) Standard k 2.45974 RNG k 2.44706 Omega SST 2.42225 The standard k, RNG k and Omega SST turbulence models have been used to simulate the inner flow in centrifugal pumps. With the same grid, the comparison among the three turbulence models was made and the results are shown in Table 3. It is found that compared with the experimental data, the head obtained by standard k model is the most accurate. Therefore, the standard k turbulence model is used to perform the simulation in this paper. For the test region, the impeller passage near the volute tongue is selected. In order to analyze the inner flow better, 7 pieces of plane curve are set in the im- peller mid-height, equidistant from the impeller inlet to the impeller outlet, and 12 points are distributed equidistantly on each curve. All the analyses of post processing in this paper are based on those monito- ring points as shown in Fig.2. Fig.2 Sample points In order to disclose the characteristics effectively, three research schemes are presented. The ratios (/dQ Q) of the flow rate are 0.8, 1.0 and 1.2, respe- ctively. 2. Results and analysis 2.1 Relative velocity distribution The relative velocity distributions are illustrated in Fig.3. It can be seen that from inlet to outlet under the same working condition, there is a low-speed zone near the middle of the pressure side at the inlet, along the direction of speed changed, which can be found in the top view of Fig.3. Then, a back flow vortex appears, with a remarkable jet-wake flow model. Furthermore, the velocities of all points in the low speed zone go upward with the increase of radius and the jet-wake flow feature becomes unnoticeable, fina- lly disappears. On the circle of impeller inlet, the rela- tive velocities decline gradually at first, from the suction side nearby to the pressure side nearby, and then, when it comes close to the pressure side, the velocity goes up again, peaking at the pressure side. With the increase of radius, the variations of velocity near the suction side are less regular than that near the pressure side. However, when starting from the middle of impeller passage, the velocity near the suction side rises gradually with the increase of radius. But the velocity near the pressure side always rise with the increase of radius. Moreover, the velocity gradient near the pressure side is larger than that near 230 the suction side. On the outlet circle, the velocity goes up gradually from the suction side. The flow in impeller passages is asymmetric due to the effect of volute. Particularly, the relative velocity of points near the volute tongue rise fastest, peaking at the patch slightly below the volute tongue. After that, they decline. In short, the values of velocity of points near the volute tongue are greater than that far from the tongue. Fig.3 relative velocity distribution in impeller passage and its local top view Under the same working condition, at the edge of the passage on the inlet circle, the velocity near the pressure side is much greater than that near the suction side. Besides, on the same circle near the inlet, the velocity near the suction side are greater than that near the pressure side. With the increase of radius, the gaps between the pressure side and suction side are narrow- wing and the velocities on those two sides tend to equal at the middle passage, and thus an equal speed area is formed. Then, the velocity near the pressure side is far greater than that near the suction side. Fig.4 Static pressure distribution in impeller passage Under different working conditions, the jet-wake model at the same position still holds. But with the increase of flow, the velocity of points in the back- ward flow vortex area also increase, while the scope reduces, and the feature with jet-wake model becomes less and less significant. The whole passage with the low flow conditions, the velocity reaches its minimum at the middle of the pressure side near the inlet, while peaking at the patch slightly below the volute tongue. Meanwhile, with the increase of flow, the greatest velocity in the whole passage decreases, the equal speed area enlarges, and gradually approaches impe- ller outlet. 231 Fig.5 Total pressure distribution in impeller passage 2.2 Static pressure and total pressure distribution The static pressure distribution and total pressure distribution of each point under different working con- ditions are illustrated in Fig.4 and Fig.5 respectively. It can be seen from Fig.4 that under the sa
收藏