[发明专利]一种多轴非比例载荷下剪应力幅值计算方法

说明书

本发明公开了一种多轴非比例载荷下剪应力幅值计算方法，即最大非比例系数法(MNFM)。多轴非比例载荷下，材料平面上的剪应力矢量方向和大小都是随着时间不断变化的，在一个多轴非比例循环载荷作用下，剪应力矢量端将描绘出一条迹线。因此，在多轴非比例载荷下，确定剪应力幅值是常困难的。本发明通过坐标转换，得到平面上的剪应力时间历程；用离散的点描述平面上剪应力矢量端所描绘的迹线；定义面上的q方向，将剪应力矢量端迹线向q方向投影，得到面上的最大投影LP；引入非比例度系数；应用MNFM公式得到面上的剪应力幅值；对经过坐标原点的所有平面进行上述计算，所有材料平面上的剪应力幅值中的最大值即为多轴非比例载荷下的剪应力幅值。

技术领域

本发明属于航空系统技术领域，具体指代一种多轴非比例载荷下剪应力幅值计算方法。

背景技术

随着航空事业的发展，新型航空器越来越多地呈现出飞行空域辽阔、飞行速度提高、飞行寿命延长等特点，对于安全性和经济性的要求也越来越高。研究表明，疲劳破坏是结构失效的主要原因之一，尤其是多轴疲劳失效。在机械零件失效中大约有80％以上属于疲劳破坏。由于疲劳破坏前往往没有明显的变形征兆，因此往往容易导致重大事故的发生。目前工程应用中的多轴疲劳问题往往采用将其简化为单轴疲劳的形式来进行处理，这样大大简化了计算步骤，但也导致疲劳寿命预测与真实寿命之间相差过大，降低了疲劳寿命预测的参考价值。因此，研究结构件在多轴加载下的疲劳寿命，寻找更准确更精简的疲劳寿命预测方法仍然是当前工程中主要课题。

在多轴疲劳损伤预测中，计算材料平面上的剪应力幅值是疲劳研究最基本，同时也是非常关键的问题之一。在真实的多轴疲劳加载过程中，部件所受到的载荷往往是变幅载荷，甚至是随机载荷。如图1所示，在一个循环的复杂加载过程中，应力向量S_n描绘了一个封闭的曲线Ψ′，其中应力向量S_n可以分解为正应力向量σ_n和剪应力向量τ,正应力σ_n只改变其大小，不改变方向。因此它的幅值和均值都较容易确定。然而，剪应力向量τ的大小和方向往往是随时间变化的，τ的尖端描绘了一个封闭的曲线Ψ，由于该曲线的形状随着外载荷的变化而变化，因此确定剪应力的幅值和均值都是一个复杂的问题。

发明内容

针对于上述现状，本发明的目的在于通过引入一个反映剪应力曲线非比例程度的因子，对于传统的最大投影法进行修正，使之适用于非比例载荷的情况。

为达到上述目的，本发明采用的技术方案如下：

本发明的一种多轴非比例载荷下剪应力幅值计算方法，包括步骤如下：

(1)定义疲劳危险点O，并将其设为坐标原点；定义直角坐标系Oxyz；定义平面Δ；定义面Δ上的局部坐标系Ouv；

(2)输入多轴应力加载历程；

(3)对于给定的平面经过坐标转换，将O点上的应力的空间矢量投影到所研究的面Δ₀内，得到在多轴非比例载荷下材料在面Δ₀上的剪应力矢量随时间变化历程；

(4)用离散的点描述这一材料在面Δ₀上剪应力矢量端在一个多轴非比例循环载荷下所描绘的迹线；

(5)在面Δ₀上定义方向q，设方向q与u轴正方向夹角为γ_i，γ_i∈[0°，180°]，将上述迹线上各个离散的点投影到q方向上，所有投影点中最远的两个点的距离设为LP_γi，使γ_i的值在0°到180°之间逐渐变化，得到Δ₀面上与局部坐标轴u轴之间夹角为0°到180°之间任意角度方向上的LP_γi，取最大投影为平面Δ₀的剪应力投影；