一种飞机结构焊点疲劳寿命计算分析方法

标题：一种飞机结构焊点疲劳寿命计算分析方法

摘要：本发明属于飞机结构疲劳寿命计算分析领域，涉及一种飞机结构焊点疲劳寿命计算分析方法，其特征在于，包括如下步骤：第一，采用模拟焊核的等效刚性梁单元，计算焊核的结构应力，取绝对值最大的主应力作为损伤参量；第二，通过材料疲劳性能试验，获取焊核和母材的S-N曲线；第三，通过准静态方法对每个计算点的有效应力历程进行计算，确定其应力谱，采用损伤累积法则计算疲劳损伤，得到所有焊点损伤和寿命的分布情况。本方案的优点是：提出适合工程应用的焊点疲劳寿命计算模型与分析方法能在设计的初期对焊点的疲劳寿命通过计算来进行寿命预测，了解整个结构上焊点的分布情况，可指导实际工艺中对焊点的个数、分布方式进行合理的调整。

申请号：CN201410154091.4

申请日：2014/4/17

申请人：中国航空工业集团公司沈阳飞机设计研究所

首项权利要求：1.一种飞机结构焊点疲劳寿命计算分析方法，其特征在于，包括如下步骤：
第一，通过有限元方法来计算结构应力时，采用模拟焊核的等效刚性梁单元，模拟焊核的刚性梁单元长度为0.5(s1+s2)，其中s1和s2分别为板1和板2的厚度，点1、2分别为梁单元在两层壳单元上的端点，点3为焊核中心线与两板连接面的交点；等效刚性梁单元传递的力和力矩被用来计算焊核和连接板周围的结构应力，建立焊接结构的细节有限元模型，通过有限元方法来计算结构应力，从计算所得的数据结果文件中分别提取点1、2、3的力Fx、Fy、Fz和弯矩Mx、My、Mz，依此来计算板1和板2内表面以及焊核在与两板交接点处的结构应力(沿焊核圆周方向每15度取一个点来计算)；点1和点2上的力和弯矩是焊核作用到板上的，而点3上的力和力矩为上层板作用于下层板的；上述结构应力计算如下：
σv1=-σmax(Fx1)cosθ-σmax(Fy1)sinθ+σ(Fz1)+σmax(Mx1)sinθ-σmax(My1)cosθ??(1)
式中：
σmax(Fx1)=Fx1/πds1；σmax(Fy1)=Fy1/πds1；σ(Fz1)=λ(1.744Fz1/s12), 当Fz1＞0时；σ(Fz1)=0，当Fz1≤0时，只有焊核轴向力中的拉伸分量引起损伤，同时：
σmax(Mx1)=λ(1.872Mx1/ds12)；σmax(My1)=(1.872My1/ds12)；
上述应力计算公式计入了通过大量试验得到的经验因子，式中λ=0.6s0.5(作为对弯曲应力梯度效应的补偿)，d为焊核直径；点2的结构应力计算与点1类似；
点3的结构应力计算用绝对最大主应力作为损伤参量，如下：
τ＝τmax(Fx3)sin2θ+τmax(Fy3)cos2θ??(2)
σ=σ(Fz3)cosθ+σmax(NX3)sinθ-σmax(My3)cosθ??(3)
式中：
τmax(Fx3)=16Fx3/(3πd2)；τmax(Fy3)=16Fy3/(3πd2)；σ(Fz3)=4Fz3/(πd2), 当Fz3＞0时；
σ(Fz3)=0，当Fz3≤0时；σmax(Mx3)=32Mx3/(πd3)；
σmax(My3)=32My3/(πd3)
面内的主应力可以从焊核中的剪应力和正应力求得：
式(4)中应力绝对值最大的主应力将作为损伤参量；
第二，焊核和母材具有不同的疲劳特性，通过材料疲劳性能试验，获取焊核和母材的S-N曲线；
第三，结合焊核和母材的材料疲劳特性；通过准静态方法对每个计算点的有效应力历程进行计算，然后，用雨流循环计数，确定其应力谱，采用损伤累积法则计算疲劳损伤，得到所有焊点损伤和寿命的分布情况。

专利类型：发明申请