[发明专利]一种基于有限元分析的管路随机振动疲劳寿命分析方法

摘要

本发明涉及一种基于有限元分析的管路动态疲劳寿命分析方法，主要通过有限元分析，获得管路随机振动条件下结构动态应力功率谱响应PSD后，然后结合损伤模型和材料S-N曲线，计算管路结构的疲劳寿命，从而判断管路结构是否处于安全工作状态，若非安全工作状态则对管路结构进行优化设计，直至结果合格，本发明分析方法解决了目前运载火箭型号管路振动疲劳寿命分析预测的问题，为管路振动疲劳寿命的设计与试验提供指导与验证；本发明分析方法可对液体运载火箭管路在随机振动载荷作用下的疲劳寿命进行有效的分析预测，为管路结构疲劳设计与试验提供指导和验证作用，该方法在载人航天运载火箭输送管的计算分析中得到了充分运用，取得了良好的效果。

主权项

一种基于有限元分析的管路随机振动疲劳寿命分析方法，其特征在于包括如下步骤：(1)利用有限元软件建立管路结构的有限元分析模型，分析管路结构在外载荷作用下的应力、应变及位移响应，分析管路结构在预应力模型作用下的自振频率振型，分析管路结构在随机振动载荷作用下的动态应力、应变、位移及加速度响应，并从中提取管路结构在静态外载荷作用下的最大Mises应力σz、随机振动载荷作用下的最大均方根应力S以及此部位的应力功率谱G(f)，并将应力功率谱G(f)输出为TXT或者DAT文件，具体步骤如下：a、利用有限元软件的模型功能根据管路结构图纸建立管路三维模型，由于管路结构为长跨度薄壁结构，模型采用壳单元模型；b、为管路三维模型的各个结构模块赋予各自的材料属性，包括弹性模量、屈服强度、破坏强度及延伸率；c、设置有限元分析步骤，主要分为三步：第一步：静态隐式分析步分析管路结构在外载荷作用下的静态应力、应变和位移响应；第二步：线性摄动分析步分析管路结构的模态振型；第三步：随机振动响应分析步分析管路结构在随机振动载荷作用下的动态应力响应谱；d、为分析模型的三个分析步分别设置边界条件及载荷条件，其中边界条件包括管路法兰及支架部位的固定和滑动约束，载荷条件包括管路内压、补偿位移及温度变化；e、分析模型网格划分：采用三维四节点壳单元对管路三维模型进行网格划分，根据管路模型大小设置单元尺寸；f、将步骤a‑e参量输入有限元软件分析模块，通过有限元软件分析模块进行分析计算，得到结果输出文件，包括管路结构在外载荷作用下的静 态应力、应变和位移响应；管路结构的模态振型；管路结构在随机振动载荷作用下的动态应力响应谱；g、利用有限元软件后处理模块对结果输出文件进行处理，提取管路结构在静态外载荷作用下的最大Mises应力σz、随机振动载荷作用下的最大均方根应力S以及此部位的应力功率谱G(f)，并将应力功率谱G(f)输出为TXT或者DAT文件；(2)建立损伤模型，采用多轴应力状态随机振动疲劳寿命Dirlik经验公式，通过损伤模型准则计算管路结构的疲劳寿命参数DL，即： D L = E ( p ) T C ( σ b σ b - σ z ) m ∫ 0 ∞ S m p ( S ) dS 其中：S为均方根应力；σb为材料的强度极限；σz为管路结构的最大Mises应力；m＝3lg(0.9/k)为疲劳曲线参数，k在弯曲时取0.5，拉压时取0.35，扭转时取0.29；C＝106(0.9σb)m为疲劳曲线参数；T为振动时间；DL值表示在T时间内元件发生疲劳破坏的可能性，1表示材料发生了破坏；E[p]与p(S)定义如下： E [ p ] = [ M 4 M 2 ] 1 2 p ( S ) = D 1 q e - z q + D 2 z R 2 e - z 2 2 R 2 + D 3 ze - z 2 2 2 M 0 式中各参量定义如下： M n = ∫ 0 ∞ G ( f ) · f n df z = S 2 M 0 D 1 = 2 ( x m - γ 2 ) 1 + γ 2 D 2 = 1 - γ - D 1 + D 1 2 1 - R D3＝1‑D1‑D2 x m = M 1 M 0 [ M 2 M 4 ] 1 2 ; R = γ - x m - D 1 2 1 - γ - D 1 + D 1 2 γ = [ M 2 2 M 0 M 4 ] 1 2 q = 1.25 ( γ - D 3 - D 2 R ) D 1 (3)根据管路结构的疲劳寿命参数DL进行判断，若DL值小于0.1，则说明管路结构处于安全工作状态；若DL值大于等于0.1，则对管路结构进行优化设计，之后重复步骤(1)‑(3)进行优化后管路结构的疲劳寿命评估，直至管路结构处于安全工作状态。