超臨界二氧化碳介質(zhì)箔片軸承彈流耦合研究
發(fā)布時(shí)間：2018-09-19 14:35
【摘要】：超臨界二氧化碳由于其特殊的傳熱以及壓縮特性,非常適合作為熱力循環(huán)的工質(zhì),近年來關(guān)于超臨界二氧化碳閉式布雷頓循環(huán)的研究正逐漸成為熱點(diǎn),而由于高轉(zhuǎn)速與較高的溫度和壓力,該循環(huán)系統(tǒng)中的渦輪-電機(jī)-壓縮機(jī)共軸轉(zhuǎn)子的支承也面臨著一系列的問題。本文針對(duì)超臨界二氧化碳閉式布雷頓循環(huán)的支承系統(tǒng),采用理論推導(dǎo)與數(shù)值計(jì)算相結(jié)合的方法,對(duì)以超臨界二氧化碳為潤(rùn)滑介質(zhì)的徑向箔片軸承和推力箔片軸承進(jìn)行仿真研究。全文共分為三部分:第一部分,使用理論方法,從連續(xù)性方程和能量方程出發(fā),得到兩組適用于變粘度變密度工況的雷諾方程,在雷諾方程中加入隨雷諾數(shù)變化的半經(jīng)驗(yàn)常數(shù),以此來考慮高轉(zhuǎn)速大雷諾數(shù)下潤(rùn)滑介質(zhì)對(duì)軸承特性的影響。對(duì)兩個(gè)雷諾方程進(jìn)行了無量綱處理,并建立合適的箔片變形模型,將潤(rùn)滑膜厚度與壓力及箔片的變形相互耦合在一起,推導(dǎo)了基于有限差分法的顯式松弛迭代的穩(wěn)態(tài)雷諾方程求解推進(jìn)格式。第二部分,使用基于MATLAB編寫的程序計(jì)算空氣介質(zhì)的箔片徑向和推力軸承的承載力與偏位角,并與文獻(xiàn)結(jié)果進(jìn)行對(duì)比,結(jié)果表明數(shù)值計(jì)算準(zhǔn)確可靠。接著結(jié)合物性計(jì)算工具使用該求解器仿真超臨界二氧化碳潤(rùn)滑下的兩種軸承的靜特性,主要包括壓力、潤(rùn)滑膜厚度和湍流潤(rùn)滑系數(shù)分布、壓力峰值、偏位角、承載力和摩擦力矩等,并討論其與相同參數(shù)傳統(tǒng)空氣相比靜特性的區(qū)別,對(duì)于推力軸承還確定了最佳箔片數(shù)。第三部分,將變粘度變密度的瞬態(tài)雷諾方程中的壓力和潤(rùn)滑膜厚度在平衡位置處基于小擾動(dòng)法線性展開,按動(dòng)態(tài)擾動(dòng)量歸并后得到和靜參數(shù)相關(guān)的剛度和阻尼相互耦合的偏微分方程,將該方程與箔片變形方程耦合,求解軸承的動(dòng)力學(xué)特性,并分析動(dòng)力學(xué)系數(shù)隨軸承偏心率、間隙量、軸承數(shù)、箔片柔度、激振頻率的變化規(guī)律。本研究工作中針對(duì)兩種軸承不同雷諾方程所共有的結(jié)構(gòu)特點(diǎn),提出一種顯式迭代求解雷諾方程的新思路,即將雷諾方程展開并化為標(biāo)準(zhǔn)形式,建立標(biāo)準(zhǔn)形式的通用迭代求解模塊,之后通過將不同雷諾方程化為標(biāo)準(zhǔn)型并調(diào)用通用模塊,可以大大簡(jiǎn)化數(shù)值編程,并增加程序可靠性。在對(duì)超臨界二氧化碳物性進(jìn)行計(jì)算時(shí),通過建立物性映射表,并輔以等距分割快速查找算法,可以快速插值求取高精度的物性參數(shù),確保計(jì)算結(jié)果的準(zhǔn)確。在軸承的動(dòng)特性計(jì)算之中,基于穩(wěn)態(tài)雷諾方程,提出了一種適用于擾動(dòng)頻率為0工況的新型動(dòng)特性計(jì)算方法,該方法與靜壓計(jì)算程序結(jié)合能夠計(jì)算變物性的軸承特性。本文所使用的方法對(duì)于其它軸承仿真研究具有一定的借鑒作用,研究得到的一些規(guī)律和結(jié)論可用于超臨界二氧化碳箔片軸承的潤(rùn)滑特性預(yù)測(cè),對(duì)箔片軸承的設(shè)計(jì)具有一定的指導(dǎo)作用。
[Abstract]:Because of its special heat transfer and compression characteristics, supercritical carbon dioxide is very suitable as a working medium for thermal cycle. In recent years, the research on closed Brayton cycle of supercritical carbon dioxide is becoming a hot topic. Because of the high speed and high temperature and pressure, the support of the coaxial rotor of the turbo-motor-compressor is faced with a series of problems. In this paper, the radial foil bearing and thrust foil bearing with supercritical carbon dioxide as lubricating medium are simulated by combining theoretical derivation and numerical calculation for the supporting system of supercritical carbon dioxide closed Brayton cycle. The paper is divided into three parts: in the first part, two sets of Reynolds equations suitable for variable viscosity and variable density conditions are obtained by using the theoretical method and the energy equation. The semi-empirical constant with Reynolds number is added to the Reynolds equation. The influence of lubricating medium on bearing characteristics at high speed and large Reynolds number is considered. The two Reynolds equations are dimensionless and a suitable foil deformation model is established. The lubrication film thickness is coupled with the pressure and the deformation of the foil. Based on the finite difference method, an explicit relaxation iteration scheme for steady-state Reynolds equation is derived. In the second part, a program based on MATLAB is used to calculate the bearing capacity and offset angle of foil radial and thrust bearings in air media, and the results are compared with the results in literature. The results show that the numerical calculation is accurate and reliable. Then the static characteristics of two kinds of bearings under supercritical carbon dioxide lubrication are simulated by using the solver combined with physical property calculation tools, including pressure, thickness of lubricating film and distribution of turbulent lubricating coefficient, peak pressure, offset angle, and so on. The bearing capacity and friction moment are discussed, and the difference between static characteristics of thrust bearing and conventional air with the same parameters is discussed. The optimum foil number is also determined for thrust bearing. In the third part, the pressure and the thickness of the lubricating film in the transient Reynolds equation with variable viscosity and density are linearly expanded at the equilibrium position based on the small perturbation method. According to the partial differential equation of stiffness and damping coupling with static parameters, the dynamic characteristic of bearing is solved by coupling the equation with foil deformation equation, and the dynamic coefficient is analyzed with the eccentricity of bearing. Change law of clearance, bearing number, foil flexibility and exciting frequency. In this work, a new idea of solving Reynolds equation by explicit iteration is proposed, which is to expand the Reynolds equation and transform it into a standard form, according to the common structural characteristics of two different Reynolds equations of bearings. The general iterative solution module in standard form is established, and the numerical programming can be greatly simplified and the program reliability can be increased by converting the different Reynolds equations into the standard form and calling the general module. In the calculation of the physical properties of supercritical carbon dioxide, by establishing the physical property mapping table and adding the fast searching algorithm of equidistant segmentation, the high precision physical parameters can be obtained quickly, and the accuracy of the calculation results can be ensured. Based on the steady state Reynolds equation, a new method for calculating the dynamic characteristics of bearings is presented in this paper. The method is combined with the static pressure calculation program to calculate the characteristics of the bearings with variable physical properties. The method used in this paper can be used as a reference for the simulation of other bearings. The obtained laws and conclusions can be used to predict the lubrication characteristics of supercritical carbon dioxide foil bearings. It can guide the design of foil bearing.
【學(xué)位授予單位】：哈爾濱工業(yè)大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2017
【分類號(hào)】：TH133.3
本文編號(hào)：2250419
[Abstract]:Because of its special heat transfer and compression characteristics, supercritical carbon dioxide is very suitable as a working medium for thermal cycle. In recent years, the research on closed Brayton cycle of supercritical carbon dioxide is becoming a hot topic. Because of the high speed and high temperature and pressure, the support of the coaxial rotor of the turbo-motor-compressor is faced with a series of problems. In this paper, the radial foil bearing and thrust foil bearing with supercritical carbon dioxide as lubricating medium are simulated by combining theoretical derivation and numerical calculation for the supporting system of supercritical carbon dioxide closed Brayton cycle. The paper is divided into three parts: in the first part, two sets of Reynolds equations suitable for variable viscosity and variable density conditions are obtained by using the theoretical method and the energy equation. The semi-empirical constant with Reynolds number is added to the Reynolds equation. The influence of lubricating medium on bearing characteristics at high speed and large Reynolds number is considered. The two Reynolds equations are dimensionless and a suitable foil deformation model is established. The lubrication film thickness is coupled with the pressure and the deformation of the foil. Based on the finite difference method, an explicit relaxation iteration scheme for steady-state Reynolds equation is derived. In the second part, a program based on MATLAB is used to calculate the bearing capacity and offset angle of foil radial and thrust bearings in air media, and the results are compared with the results in literature. The results show that the numerical calculation is accurate and reliable. Then the static characteristics of two kinds of bearings under supercritical carbon dioxide lubrication are simulated by using the solver combined with physical property calculation tools, including pressure, thickness of lubricating film and distribution of turbulent lubricating coefficient, peak pressure, offset angle, and so on. The bearing capacity and friction moment are discussed, and the difference between static characteristics of thrust bearing and conventional air with the same parameters is discussed. The optimum foil number is also determined for thrust bearing. In the third part, the pressure and the thickness of the lubricating film in the transient Reynolds equation with variable viscosity and density are linearly expanded at the equilibrium position based on the small perturbation method. According to the partial differential equation of stiffness and damping coupling with static parameters, the dynamic characteristic of bearing is solved by coupling the equation with foil deformation equation, and the dynamic coefficient is analyzed with the eccentricity of bearing. Change law of clearance, bearing number, foil flexibility and exciting frequency. In this work, a new idea of solving Reynolds equation by explicit iteration is proposed, which is to expand the Reynolds equation and transform it into a standard form, according to the common structural characteristics of two different Reynolds equations of bearings. The general iterative solution module in standard form is established, and the numerical programming can be greatly simplified and the program reliability can be increased by converting the different Reynolds equations into the standard form and calling the general module. In the calculation of the physical properties of supercritical carbon dioxide, by establishing the physical property mapping table and adding the fast searching algorithm of equidistant segmentation, the high precision physical parameters can be obtained quickly, and the accuracy of the calculation results can be ensured. Based on the steady state Reynolds equation, a new method for calculating the dynamic characteristics of bearings is presented in this paper. The method is combined with the static pressure calculation program to calculate the characteristics of the bearings with variable physical properties. The method used in this paper can be used as a reference for the simulation of other bearings. The obtained laws and conclusions can be used to predict the lubrication characteristics of supercritical carbon dioxide foil bearings. It can guide the design of foil bearing.
【學(xué)位授予單位】：哈爾濱工業(yè)大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2017
【分類號(hào)】：TH133.3
【參考文獻(xiàn)】
相關(guān)期刊論文 前10條
1 閆佳佳;劉占生;張廣輝;韋繼勤;;考慮稀薄氣體效應(yīng)的止推箔片軸承靜特性分析[J];航空動(dòng)力學(xué)報(bào);2016年06期
2 吳毅;王佳瑩;王明坤;戴義平;;基于超臨界CO_2布雷頓循環(huán)的塔式太陽能集熱發(fā)電系統(tǒng)[J];西安交通大學(xué)學(xué)報(bào);2016年05期
3 郭偉;戴俊;李國(guó)平;章藝;王小靜;;徑向可傾瓦軸承動(dòng)特性特征參數(shù)影響研究[J];船舶工程;2014年S1期
4 陳汝剛;陳韜;龔超;;箔片動(dòng)壓止推氣體軸承流固耦合數(shù)值模擬[J];西安交通大學(xué)學(xué)報(bào);2014年05期
5 黃彥平;王俊峰;;超臨界二氧化碳在核反應(yīng)堆系統(tǒng)中的應(yīng)用[J];核動(dòng)力工程;2012年03期
6 段承杰;王捷;楊小勇;;反應(yīng)堆超臨界CO_2 Brayton循環(huán)特性[J];原子能科學(xué)技術(shù);2010年11期
7 陳汝剛;周權(quán);劉燁;侯予;;新型箔片動(dòng)壓止推氣體軸承承載特性的試驗(yàn)研究[J];西安交通大學(xué)學(xué)報(bào);2010年09期
8 許懷錦;劉占生;張廣輝;王永亮;;波箔片氣體動(dòng)壓徑向軸承氣膜壓力分布特點(diǎn)[J];軸承;2008年06期
9 虞烈,戚社苗,耿海鵬;彈性箔片空氣動(dòng)壓軸承的完全氣彈潤(rùn)滑解[J];中國(guó)科學(xué)E輯:工程科學(xué) 材料科學(xué);2005年07期
10 熊聯(lián)友,王 瑾,劉井龍,陳純正;平箔式箔片止推氣體軸承靜特性的理論研究[J];潤(rùn)滑與密封;2002年02期
相關(guān)博士學(xué)位論文 前2條
1 徐方程;基于庫(kù)倫摩擦效應(yīng)的徑向氣體箔片軸承-轉(zhuǎn)子系統(tǒng)振動(dòng)特性[D];哈爾濱工業(yè)大學(xué);2013年
2 張廣輝;高速動(dòng)靜壓混合氣體軸承轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)特性研究[D];哈爾濱工業(yè)大學(xué);2010年
相關(guān)碩士學(xué)位論文 前3條
1 胡錦;嵌入彈簧式氣體箔片軸承理論分析與實(shí)驗(yàn)研究[D];湖南大學(xué);2016年
2 浮翔;新型跨臨界CO_2冷熱電聯(lián)合循環(huán)性能研究與分析[D];重慶大學(xué);2014年
3 閆佳佳;彈性箔片動(dòng)壓氣體止推軸承-轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)特性研究[D];哈爾濱工業(yè)大學(xué);2012年
，本文編號(hào)：2250419
本文鏈接：http://www.sikaile.net/jixiegongchenglunwen/2250419.html
最近更新
教材專著
