基于格心型有限体积法的喷管热应力分析  

Nozzle thermal stress analysis based on cell-centred finite volume method

作　　者：马源 董瑞 陈晨 谭志城 梁冰 陈亚楼 明平剑 MA Yuan;DONG Rui;CHEN Chen;TAN Zhicheng;LIANG Bin;CHEN Yalou;MING Pingjian(Sino-French Institute of Nuclear Engineering and Technology,Sun Yat-sen University,Zhuhai 519082,China;National Key Laboratory of Solid Rocket Propulsion,The Institute of Xi'an Aerospace Solid Propulsion Technology,Xi'an 710025,China;College of Power and Energy Engineering,Harbin Engineering University,Harbin 150001,China)

机构地区：[1]中山大学中法核工程与技术学院,珠海519082 [2]西安航天动力技术研究所,固体推进全国重点实验室,西安710025 [3]哈尔滨工程大学动力与能源工程学院,哈尔滨150001

出　　处：《固体火箭技术》2025年第2期299-309,共11页Journal of Solid Rocket Technology

摘　　要：基于格心型有限体积法(CC-FVM),以位移作为求解变量,建立CC-FVM弹性方程求解方法;针对非线性接触问题,采用位移插值方法计算接触穿透量,根据罚函数法和库伦摩擦模型,推导出法向应力和摩擦力的计算方法,计算得出接触应力和界面摩擦力;采用迭代方法实现非线性问题求解,并通过二维滑块和三维滑块算例验证了新接触应力计算方法的正确性。基于该方法,对某复合材料喷管进行冷态和热负荷工况的应力和接触应力分析,分析温度载荷对喷管应力分布的影响,计算得出喷管的应力分布和接触应力的集中区域,并提出了降低喷管内部接触应力的措施。结果表明,金属层和隔热层在收敛段几何折角处接触应力最大,正常工况下喷管最大应力位于固定层与螺栓连接的倒几何折角处,小于材料极限强度,热负荷使得喷管固定层的应力增大37.5%,非金属层应力增大40%;喷管部件连接处的平滑处理可有效降低这些区域的接触应力。Based on the cell-centred finite volume method(CC-FVM),analogous to the finite volume method(FVM)used in computational fluid dynamics(CFD),the CC-FVM elasticity equations were established according to the principle of force equilibrium by taking the displacements as the solution variables.Aiming at the nonlinear contact problem,the displacement interpolation method was used to calculate the contact penetration.According to the penalty function method and Coulomb friction model,the calculation method of normal stress and friction force was deduced,and the contact stress and interface friction force were calculated.The iterative method was used to solve the nonlinear problem,and the correctness of the new contact stress calculation method was verified by two-dimensional slider and three-dimensional slider examples.Based on this method,the stress and contact stress of a composite material nozzle under cold and hot load conditions were analyzed;the influence of temperature load on the stress distribution of the nozzle was analyzed;the stress distribution of the nozzle and the concentration area of the contact stress were calculated,and the measures to reduce the contact stress in the nozzle were put forward.The results show that the maximum contact stress between the metal layer and the heat insulation layer occurs at the chamfer of the convergence section.Under normal conditions,the maximum stress of the nozzle is located at the chamfer of the fixed layer and the bolt connection,which is less than the ultimate strength of the material.The thermal load increases the stress of the nozzle fixed layer by 37.5%,and the stress of the non-metal layer increases by 40%.The smooth treatment of the nozzle component joints can effectively reduce the contact stress in these areas.

关 键 词：有限体积法 固体火箭发动机 喷管 接触应力 弹性力学 

分 类 号：V414[航空宇航科学与技术—航空宇航推进理论与工程]