摘要：空间核电系统采用高温气冷堆闭式布雷顿循环方案，使动力转换系统可高效地利用核能发电。动力转换系统采用止推箔片轴承，省略了供油和供气系统。建立柱坐标系下可压缩气体的雷诺方程并耦合求解锥度和不对中误差影响的气膜厚度方程，分析了止推箔片轴承在生产和装配过程中产生的锥度和不对中误差对轴承承载特性的影响。结果表明：在设计间隙不变的情况下，锥度增大导致轴承承载力减小；在轴向载荷恒定的情况下，锥度绝对值越大，最小气膜厚度越小，最大气膜压力则相反；随着不对中角度增加，轴承承载力和摩擦功快速增大，单瓦的承载力出现严重不均衡情况，不对中误差对轴承刚度系数和阻尼系数的影响较明显，直接影响转子系统的稳定性。

　　人类探测太空的要受空间推进技术的限制,针对高速度增量的航天任务,传统化学推进技术的比冲较低,太阳能电池阵—蓄电池技术严重依赖太阳光照,深空生存能力较差[1]。空间核电源系统具有能量密度大,容易实现大功率(数千瓦至数兆瓦)供电,机动性和隐蔽性好,不依赖太阳光照等优势,是深空探测不可替代的空间电源[1]。动态转换的放射性同位素发电系统在空间核电源系统中转换效率最高,其中研究较多且已应用的是闭式布雷顿循环的动力转换系统(图1),旋转单元是该系统的核心组件,其结构和性能直接影响系统的优劣。止推箔片轴承作为一种柔性自适应动压轴承,在旋转单元中承受轴向载荷,与滚动轴承相比,省略了复杂的冷却润滑系统,使系统更加清洁和高效。(a)闭式布雷顿循环系统原理(b)旋转单元图1闭式布雷顿循环的动力转化系统及核心组件Fig.1PowerconversionsystemandcorecomponentsofclosedBraytoncycle1854年,法国科学家提出将气体作为润滑剂的可能,1886年法国科学家雷诺推导出润滑膜内压力分布的雷诺方程[2],对流体润滑的认识上升到理论高度。止推箔片轴承作为阶梯型动压轴承的衍生结构轴承,具有波箔刚度非线性的特点,其理论研究晚于径向箔片轴承且相关研究较少。1983年,文献[3]对止推箔片轴承的静态特性进行了详细的计算分析,计算模型考虑了波箔刚度的非线次使用柔度系数α分析箔片轴承静特性,该方法一直沿用至今。2000年,文献[4]采用有限差分法与有限元相结合的方法,将可压缩流体的雷诺方程求解与弹性变形耦合起来进行止推箔片轴承分析。2017年,文献[5]采用非等温雷诺方程评估了转子静态和动态角度不对中对止推箔片轴承性能的影响,并深入分析了轴承的承载力、摩擦功、气膜刚度和阻尼系数的变化情况。20世纪80年代,文献[6]应用圆柱薄壳线性理论推导出波箔变形方程,这是国内最早关于箔片轴承中波箔刚度的计算模型。文献[7]采用有限单元法对平箔式止推箔片轴承进行求解,并分析了箔片轴承静态特性的影响因素。文献[8]对止推箔片轴承进行了启停试验,分析了箔片材料和表面涂层对轴承启停磨损的影响。文献[9]采用有限差分法得出,推力盘倾斜状态下轴承的轴向承载力和黏性摩擦力矩均随倾角的增大而增大,且倾角越大静态特性越大。文献[10]基于牛顿迭代法建立偏角不对中时止推箔片轴承的润滑模型,与线性替换法的比较说明牛顿迭代法具有同样的适用性。文献[11]基于CFX的有限元计算与MATLAB的有限差分求解动压气体轴承气膜压力,二者计算结果一致,随着偏心率的增加CFX有限元的计算结果比有限差分法计算的略小。文献[12]采用ANSYS和CFX软件计算了超临界二氧化碳止推箔片轴承动态特性,得出动态刚度和阻尼大小主要取决于气膜厚度和箔片结构,且随气膜厚度的增大而减小。随着更多类型箔片轴承的研究,箔片轴承的性能逐渐被熟知,但仍缺乏精准的理论分析模型和实用的设计规范[13]。综上可知,目前在止推箔片轴承的误差研究方面,主要着重于轴承不对中误差对其静态性能的影响,而锥度误差对轴承承载特性的影响以及不对中误差下轴承动态特性变化的研究较少。在生产和装配过程中,锥度和不对中误差均难以避免,不对中误差主要有平行不对中、角不对中和综合不对中。由于旋转单元转子系统的径向设计间隙约束及推力盘径向尺寸大于止推轴承高度,分析过程不考虑平行不对中的影响,本文所提及的不对中特指角不对中或主轴倾斜。本文在柱坐标系下建立止推箔片轴承的雷诺方程,并采用牛顿迭代法和有限差分法对雷诺方程和气膜厚度方程进行耦合求解,分析生产和装配过程引起的径向锥度误差和不对中误差对止推箔片轴承承载特性的影响。

　　止推箔片轴承主要由顶箔、波箔和底板组成。波箔和顶箔组成多个扇形瓦,每个瓦的度数和倾斜面占比相同。止推箔片轴承结构和箔片参数示意图如图2所示,β为扇形瓦的瓦角,βp为扇形瓦节线位置角,σh为扇形瓦倾斜面高度,ω为推力盘角速度,h2为设计间隙。(a)止推箔片轴承(b)箔片参数图2止推箔片轴承结构和箔片参数示意图Fig.2Diagramofthrustfoilbearingstructureandfoilparameters

　　推导等温条件下的雷诺方程时进行如下假设:1)忽略气体惯性力;2)润滑气体介质视为牛顿流体,符合牛顿黏性定律;3)气膜间隙内为层流;4)忽略气膜法向流动;5)润滑介质视为理想气体。经推导得到柱坐标下润滑介质为理想气体的雷诺方程,即(1)式中:r为径向坐标;h为气膜厚度;P为气膜压力;θ为周向坐标;μ为气体动力黏度。对上述参数进行量纲一化处理,令对于恒温定常状态雷诺方程,(1)式可简化为(2)式中:Pa为环境压力;R2为轴承外半径;Λ为轴承数。

　　每个扇形瓦的气膜由倾斜面气膜和水平面气膜组成,根据图2的结构参数得到单个扇形瓦的气膜厚度方程,即h=h2+h(θ)+w(r,θ),(3)式中:w(r,θ)为箔片变形后构建的变形区域气膜厚度;h(θ)为扇形瓦倾斜面气膜厚度。定义b为斜面占比,b=βP/β,气膜厚度的表达式为(4)在工程应用过程中存在加工误差,推力盘不可避免地产生一定的径向锥度βz,将直接影响气膜厚度分布。锥度为负数时表示推力盘出现内锥状态,为正数时出现外锥状态,按内锥度和外锥度这2种形式构建径向锥度模型(图3),Ri为锥度半径,R1为轴承内半径,Δh为锥度引起的气膜厚度变化。锥度βz的定义为(a)内锥度(b)外锥度图3径向锥度模型Fig.3Radialtapermodel(5)同样,对相关参数进行量钢一化处理,令建立考虑推力盘径向锥度的气膜厚度方程,即h0=h2+h(θ)+w(r,θ)+Δh。(6)尽管设计及生产过程要严格控制形位公差,但在装配过程仍会存在不对中现象[14-15],而不对中误差会导致止推箔片轴承气膜间隙发生变化,从而引起轴承承载特性的改变。如图4所示,建立不对中模型和空间坐标系。(a)坐标系及扇形瓦标记(b)坐标系及不对中模型图4空间坐标系和不对中模型Fig.4Spatialcoordinatesystemandmisalignmentmodel假设转子与止推箔片轴承的夹角(不对中角度)为φ,在空间坐标下可分解为绕x轴的θx和绕y轴的θy。考虑转子不对中的气膜厚度方程为h1=h2+h(θ)+w(r,θ)-rcos?θtan?φ?。(7)

　　波箔单位面积刚度Kb的表达式为(8)式中:Eb为波箔的弹性模量;t为波箔厚度;l为半波箔长度;s为波箔单位长度;ν为波箔的泊松比。

　　将单个扇形瓦按图5进行网格划分,沿径向划分为m份,沿周向划分为n份。扇形瓦周围与环境接触,其边界条件为环境边界。量纲一的边界条件为图5单个扇形瓦的网格划分Fig.5Gridpartitionofsinglesegmentalpad

　　为验证止推箔片轴承数学模型的准确性,按文献[3]的算例进行量纲一化处理和计算。文献[3]中量纲一的承载力为(10)式中:W为承载力。本文程序计算与文献[3]的对比结果如图7所示,最大误差为2.53%,两者吻合度较好,说明本文止推箔片轴承数学模型的准确性较高。图7本文程序与文献[3]计算结果对比Fig.7ComparisonofcalculationresultsbetweenthisprogramandRef.[3]

　　锥度一定时,各扇形瓦气膜厚度和气膜压力在同一半径和周向角度上相同,为提高计算精度并减少计算时间,选取单个扇形瓦进行分析。在额定设计间隙下,内锥(Ri=R1)和外锥(Ri=R2)状态下量纲一的气膜压力和气膜厚度分布如图8和图9所示:当内锥度为-4×10-4时,最大气膜压力增大18.4%,对应轴向承载力提高47.8%;当外锥度为4×10-4时,最大气膜压力减小12.7%,对应轴向承载力降低33.8%。以图9b中气膜厚度A和B位置为基准,内锥度为-4×10-4时B区域的气膜厚度减小至10.6μm,外锥度为4×10-4时A区域的气膜厚度增大至29.4μm。综上可知:内锥度减小了平均气膜厚度,提高了承载力;外锥度增大了平均气膜厚度,降低了承载力。(a)-4×10-4

　　图9不同锥度下止推箔片轴承的气膜厚度分布Fig.9Gasfilmthicknessdistributionforthrustfoilbearingunderdifferenttapers设计间隙不变时轴向承载力和摩擦功随锥度的变化以及轴向载荷恒定(200N)时最小气膜厚度和量纲一的最大气膜压力随锥度的变化如图10所示:设计间隙不变时,轴向承载力和摩擦功随着锥度的增加逐渐减小,由于锥度的增加增大了平均气膜厚度,削弱了动压效应进而减小了轴向承载力和摩擦功。在轴向载荷恒定时,随着锥度绝对值的增加,最小气膜厚度逐渐减小,最小气膜厚度达到0.67μm,若继续减小至零时程序计算很难收敛,由于模型未考虑表面粗糙度,工程应用时可认为此时轴承已处于临界碰磨状态。说明恒定载荷工况时只有允许一定范围的锥度误差才能保证轴承正常运行。最大气膜压力与最小气膜厚度正好相反,说明轴向载荷恒定时只要锥度绝对值增加就会增大最大气膜压力。此外,内锥度和外锥度引起的最大气膜压力和最小气膜厚度均关于锥度βz=0对称。?图10止推箔片轴承静态参数随锥度的变化Fig.10Changeofstaticparametersforthrustfoilbearingwithtaper

　　为简化分析模型,掌握不对中对轴承性能的影响,只选取绕某一轴倾斜产生的不对中进行分析。对中和不对中时止推箔片轴承的气膜压力、气膜厚度和箔片变形如图11所示:由于不对中的产生,轴承各瓦的气膜压力和气膜厚度分布不均匀,箔片产生的变形量也不同;2#瓦最大气膜压力(2.458)最大,高于对中状态下各瓦的最大气膜压力(1.805),5#瓦的最大气膜压力最小(1.494);由于气膜压力作用,最大变形量为(1.205μm)同样出现在2#瓦上,是5#瓦最大变形量(0.412μm)的2.925倍;由于箔片变形构建了新的气膜承载区域,气膜压力越高时箔片变形量越大,体现了止推箔片轴承较好的自适应性。??

　　图11止推箔片轴承的气膜压力、气膜厚度和箔片变形量Fig.11Gasfilmpressure,gasfilmthicknessandfoildeformationofthrustfoilbearing止推箔片轴承和各瓦的承载力和摩擦功随不对中角度的变化如图12所示:轴承的承载力和摩擦功随着不对中角度的增大而增加,但各瓦不同;当推力盘按图4坐标系倾斜时,x轴上方的3个扇形瓦平均气膜厚度减小,动压效应增强,气膜压力升高,承载力和摩擦功随不对中角度的增大而增加,2#瓦平均气膜厚度最小,其承载力最大且增速最快,1#瓦次之,3#瓦增速最慢,x轴下方3个扇形瓦则相反;不对中造成各瓦的承载力和摩擦发热严重不均衡,不对中角度越大时各瓦的差异越大。

　　(b)各瓦的承载力(c)各瓦的摩擦功图12止推箔片轴承和各瓦的承载力和摩擦功随不对中角度的变化Fig.12Change

　　2、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。

　　3、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。

　　4、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

　　2025年一级造价师之建设工程技术与计量（安装）题库【精选题】.docx

　　2025年一级造价师之建设工程技术与计量（安装）题库含完整答案【典优】.docx

　　2025年一级造价师之建设工程技术与计量（安装）题库【研优卷】.docx

　　2025年企业人力资源管理师之一级人力资源管理师考试题库附参【综合卷】.docx

　　2025年一级造价师之建设工程技术与计量（安装）题库附参【综合题】.docx

　　2025年一级造价师之建设工程技术与计量（安装）题库及答案参考.docx

　　2026届吉林省梅河口五中等联谊校生物高三第一学期期末质量跟踪监视试题含解析.doc

　　市场营销学（第三版）教学课件第7章 4Ps营销组合策略（上）PPT.pptx

　　原创力文档创建于2008年，本站为文档C2C交易模式，即用户上传的文档直接分享给其他用户（可下载、阅读），本站只是中间服务平台，本站所有文档下载所得的收益归上传人所有。原创力文档是网络服务平台方，若您的权利被侵害，请发链接和相关诉求至 电线) ，上传者