一种基于有限元的高压气缸动态热应力分析方法技术

本发明专利技术属机械结构有限元分析相关技术领域，并公开了一种基于有限元的高压气缸动态热应力分析方法。该方法包括下列步骤：S1模拟并获得瞬态运动对应的温度分布模拟值和应力分布模拟值；利用温度分布模拟值和应力分布模拟值预测稳态运动，获得稳态运动中不同周期的温度分布预测值和应力分布预测值；S2计算并获得稳态运动中各个周期对应的温度分布的计算值和应力分布的计算值；S3分别计算稳态运动各个周期的温度分布预测值和计算值之间以及应力分布的预测值和计算值之间的差值，以此获得预测值和计算值之间的误差，实现热应力分析。通过本发明专利技术，实现瞬态运动的热应力的模拟，有效反映实际的瞬态运动，为实践提供指导意义。为实践提供指导意义。为实践提供指导意义。

【技术实现步骤摘要】
一种基于有限元的高压气缸动态热应力分析方法


[0001]本专利技术属机械结构有限元分析相关
，更具体地，涉及一种基于有限元的高压气缸动态热应力分析方法。

技术介绍

[0002]自20世纪以来，天然气在能源结构中所占的比例越来越高，在现代工业中有着十分广阔的发展前景。天然气压缩机广泛应用于排水和天然气出产、集气处置、管道运输、贮存等范畴。在工作中，气缸内部承受压缩气体带来的高温高压，工作环境非常恶劣，高压气缸在运行中受到的热应力已经成为高压气缸故障的重要原因之一。但是目前国内外对高压气缸动态热应力方面的研究甚少，也缺乏气缸动态热应力规律的总结，忽略了高压气缸故障的重要原因之一，从而误判其他原因才是导致高压气缸故障的根本原因，无法从故障中吸取经验教训。
[0003]高压气缸作为压缩机的重要部件，高压气缸的安全运行是压缩机系统稳定的必要条件。高压气缸运行时的高温及产生的热应力会导致高压气缸变形，严重时会直接导致气缸故障，会提高压缩机运行时故障的发生率。因此，准确的计算出高压气缸运行时各个周期的热应力及其作用规律，降低压缩机的故障率，对高压气缸的安全稳定工作起到关键的作用。
[0004]在目前其他关于高压气缸的研究中，更多的关注于气缸的稳态热应力或者与其他稳态载荷的耦合关系，并没有给出高压气缸运行时各个周期的动态热应力计算方法，也缺乏对动态热应力与稳态热应力间的关系总结，因此，本领域存在着发展一种基于有限元的高压气缸动态热应力分析方法的技术需求。

技术实现思路

[0005]针对现有技术的以上缺陷或改进需求，本专利技术提供了一种基于有限元的高压气缸动态热应力分析方法，其中通过利用瞬态运动温度和应力的模拟值来预测稳态运动温度和应力的预测值，然后利用预测值和稳态运动的计算值进行比较获得误差，当误差较大时，说明瞬态运动的模拟不符合实际，当误差较小时，说明瞬态运动的模拟贴合实际，实现对瞬态运动热应力的分析，有效反映真实的瞬态运动，为实践操作提供指导意义。
[0006]为实现上述目的，按照本专利技术，提供了一种基于有限元的高压气缸动态热应力分析方法，其特征在于，该方法包括下列步骤：
[0007]S1瞬态运动分析
[0008]对高压气缸的有限元模型施加动态热边界条件，设定不同的高压气缸模型自由表面初始温度，以此获得高压气缸在不同自由表面初始温度下，瞬态运动中各个周期对应的温度分布模拟值和应力分布模拟值；
[0009]将不同自由表面初始温度下的瞬态运动的温度分布模拟值和应力分布模拟值分别拟合为温度分布函数和应力分布函数，将该温度分布函数和应力分布函数作为预测函数
预测稳态运动周期中的温度分布场和应力分布，以此获得不同自由表面初始温度下的稳态运动中不同周期的温度分布预测值和应力分布预测值；
[0010]S2稳态运动分析
[0011]对高压气缸有限元模型施加稳态热边界条件，按照所瞬态运动设定的自由表面初始温度进行温度和应力的计算，获得在不同自由表面初始温度下的稳态运动中各个周期对应的温度分布的计算值和应力分布的计算值；
[0012]S3计算误差
[0013]在不同自由表面初始温度下，分别计算稳态运动各个周期的温度分布预测值和计算值之间以及应力分布的预测值和计算值之间的差值，以此获得预测值和计算值之间的误差，实现热应力分析。
[0014]进一步优选地，在步骤S1中，根据不同自由表面温度下瞬态运动各个周期对应的应力分布模拟值，确定高压气缸模型中的危险区域，其中，危险区域为每个周期中应力峰值所对应的高压气缸模型中的区域。
[0015]进一步优选地，在步骤S1中，所述动态热边界条件的设定包括：将所述高压缸的有限元模型设定为绝热模型，高压气缸的进排气腔满足第三类强制对流换热边界条件，高压气缸的自由表面满足第三类自然对流换热边界条件，高压气缸的压缩腔满足第三类强制对流换热边界条件和Woschni半经验公式得到的随周期变化的热边界条件。
[0016]进一步优选地，在步骤S1中，求解所述应力场的分布模拟值采用热-机耦合计算的方法。
[0017]进一步优选地，在步骤S1中，求解应力分布模拟值中，设定高压气缸的有限元模型的缸体法兰端面的边界条件是轴向位移分量为0，自由表面满足自由边界条件。
[0018]进一步优选地，在步骤S2中，不同自由表面初始温度下，稳态运动中的温度分布计算值和应力分布的计算值相同。
[0019]进一步优选地，在步骤S3中，在计算误差时，在高压气缸的整个运动周期中选择一个周期i计算误差即可，其中周期i的选择按照下列方式进行：
[0020]首先，分别构建周期与温度和应力的二维平面图，其中，每个周期中选择该周期中的温度最大值或者应力最大值作为二维平面图中该周期对应的温度或应力；
[0021]其次，周期与温度以及周期与应力的二维平面图中曲率变化均为0对应的周期即为所需的周期i。
[0022]总体而言，通过本专利技术所构思的以上技术方案与现有技术相比，具备下列有益效果：
[0023]1.本专利技术中通过构建瞬态运动的边界条件，对瞬态运动进行模拟，再根据模拟值预测稳态运动获得预测值，最后将稳态运动的预测值与计算值进行比较，当误差较小时，说明本专利技术采用的方法对瞬态运动的模拟与实际相符合，当误差较大是，说明本专利技术对瞬态运动的模拟与实际相差较大；采用上述方法能有效对高压气缸刚启动时的瞬态运动进行模拟和分析；
[0024]2.本专利技术中通过对瞬态运动的温度和应力分布的模拟值构建温度和应力的拟合函数，并以此作为稳态运动的预测函数，通过该方式建立了瞬态运动和稳态运动的关系，最终可以通过瞬态运动的预测值与计算值之间的误差反映瞬态运动的模拟值是否准确，方法
简单有效，预测精度高；
[0025]3.本专利技术通过根据瞬态运动中各个周期对应的应力分布的最大值对应的高压缸中的位置，确定高压缸中最危险的区域，根据该最危险的区域获知实际高压缸运动中哪个部位承受的应力最大，进而为后续的防范工作起到指导作用。
附图说明
[0026]图1是按照本专利技术的优选实施例所构建的有限元的动态热应力分析方法的流程示意图；
[0027]图2是按照本专利技术的优选实施例所构建的高压缸的有限元模型；
[0028]图3是按照本专利技术的优选实施例所构建的高压气缸的压缩腔的动态热边界条件曲轴转角和换热系数的对应关系。
具体实施方式
[0029]为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。此外，下面所描述的本专利技术各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0030]高压气缸的运动分为启动初期的瞬态运动和后期的稳态运动，稳态运动由于运动状态稳定，其对应的温度和应力可以直接计算获得，但是对于初期的瞬态运动，其运动没有规律，无法直接获得温度和应力的分布，但是实际本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于有限元的高压气缸动态热应力分析方法，其特征在于，该方法包括下列步骤：S1瞬态运动分析对高压气缸的有限元模型施加动态热边界条件，设定不同的高压气缸模型自由表面初始温度，以此获得高压气缸在不同自由表面初始温度下，瞬态运动中各个周期对应的温度分布模拟值和应力分布模拟值；将不同自由表面初始温度下的瞬态运动的温度分布模拟值和应力分布模拟值分别拟合为温度分布函数和应力分布函数，将该温度分布函数和应力分布函数作为预测函数预测稳态运动周期中的温度分布场和应力分布，以此获得不同自由表面初始温度下的稳态运动中不同周期的温度分布预测值和应力分布预测值；S2稳态运动分析对高压气缸有限元模型施加稳态热边界条件，按照所瞬态运动设定的自由表面初始温度进行温度和应力的计算，获得在不同自由表面初始温度下的稳态运动中各个周期对应的温度分布的计算值和应力分布的计算值；S3计算误差在不同自由表面初始温度下，分别计算稳态运动各个周期的温度分布预测值和计算值之间以及应力分布的预测值和计算值之间的差值，以此获得预测值和计算值之间的误差，实现热应力分析。2.如权利要求1所述的一种基于有限元的高压气缸动态热应力分析方法，其特征在于，在步骤S1中，根据不同自由表面温度下瞬态运动各个周期对应的应力分布模拟值，确定高压气缸模型中的危险区域，其中，危险区域为每个周期中应力峰值所对应的高压气缸模型中的区域。3.如权利要求1或2所述的一种基于有限元的高压气缸动态热应力分析方法，其特征...

【专利技术属性】
技术研发人员：余祖耀，刘彪，魏美容，曾雨昊，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：