CN104020062A - 一种编制疲劳加速试验载荷谱的方法 - Google Patents

Abstract

一种编制疲劳加速试验载荷谱的方法，该方法有三大步骤：步骤一、疲劳载荷循环的判别；步骤二、低载截除；步骤三、归并二级波，直至无法再次合并。本发明简单实用、操作方便、损伤等效精度高、时间压缩效果显著。本发明在测试技术领域里具有较好的实用价值和广阔的应用前景。

Description

一种编制疲劳加速试验载荷谱的方法

技术领域

本发明提供一种编制疲劳加速试验载荷谱的方法，属于试验测试技术领域。

背景技术

编制疲劳加速试验载荷谱是结构疲劳验证试验的重要前提，由于实测载荷谱往往包含大量的高频低载，施加此类载荷对零部件进行全尺寸试验所需的时间常常占总试验时间的大部分，造成极大的时间、人力、物力资源浪费，因此，编制疲劳加速试验载荷谱，缩短试验时间，能极大降低试验成本。本发明提出一种编制疲劳加速试验载荷谱的方法，在保持实测载荷谱中大载荷的同时，截除大量不造成疲劳损伤的高阶小载循环，并按疲劳损伤等效原理归并大量造成一定损伤的中等大小的载荷循环，这样编制的疲劳加速试验载荷谱既能保持载荷顺序效应，损伤等效精度高，又能有效地缩短试验时间，具有重要的学术意义和工程应用价值。

发明内容

1、目的：本发明的目的是提供一种编制疲劳加速试验载荷谱的方法，该方法具有简单实用、操作方便、损伤等效精度高、时间压缩效果显著等优点。

2、技术方案：本发明一种编制疲劳加速试验载荷谱的方法，该方法具体步骤如下：

步骤一、疲劳载荷循环的判别

从工程观点出发，通常取N_f＝10⁷次循环和N_T＝10⁴次循环对应的疲劳应力循环的水平分别作为三级波和主波的判别门槛值，即裂纹形成寿命长于N_f所对应的应力循环可判别为三级波，而裂纹形成寿命短于N_T的应力循环则判别为主波。根据S-N曲面公式

${(\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-S_m}{S}_a-S_0)}^{m}N=C---\left(1\right)$

可以得出

${\mathrm{\σ}}_b{s}_a+[{\left(\frac{C}{N_f}\right)}^{\frac{1}{m}}+S_0]s_m\≤[{\left(\frac{C}{N_f}\right)}^{\frac{1}{m}}+S_0]{\mathrm{\σ}}_b---\left(2\right)$

${\mathrm{\σ}}_b{s}_a+[{\left(\frac{C}{N_T}\right)}^{\frac{1}{m}}+S_0]s_m\≥[{\left(\frac{C}{N_T}\right)}^{\frac{1}{m}}+S_0]{\mathrm{\σ}}_b---\left(3\right)$

式中符号说明如下：σ_b为材料强度极限，s_a为应力幅值，s_m为应力均值，S₀为材料拟合疲劳极限，N为疲劳寿命，C和m为常数。

当应力循环满足式(2)时，则判别为三级波；而当应力循环满足式(3)时，则判别为主波；当式(2)和式(3)都不满足时，应力循环判别为二级波。

步骤二、低载截除

保留主波应力循环不变，并截除三级波应力循环。

步骤三、归并二级波，直至无法再次合并。

根据Miner线性累积损伤理论作为损伤当量折算的基础，一个应力循环(s_a,s_m)造成的损伤可写成

$D(s_a,s_m)=\frac{1}{N(s_a,s_m)}---\left(4\right)$

式中，N(s_a,s_m)由式(1)确定。为了将两个相邻的连续应力循环(s_a1,s_m1)和(s_a2,s_m2)归并为新二级波，根据损伤当量原则，由式(4)和式(1)可建立疲劳裂纹形成的损伤当量折算公式：

${[\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-{\left(s_m\right)}_{\mathrm{eq}}}{\left(s_a\right)}_{\mathrm{eq}}-S_0]}^{-m}={[\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m1}}{s}_{a1}-S_0]}^{-m}+{[\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m2}}{s}_{a2}-S_0]}^{-m}---\left(5\right)$

令(s_min)_eq＝s_min2＝s_m2-s_a2， $M={[{(\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m1}}{s}_{a1}-S_0)}^{-m}+{(\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m2}}{s}_{a2}-S_0)}^{-m}]}^{-\frac{1}{m}}+S_0,$ 则由式(5)可分别得到归并后新二级波的应力幅值和均值分别为：

${\left(s_a\right)}_{\mathrm{eq}}=\frac{({\mathrm{\σ}}_b+s_{a2}-s_{m2})M}{{\mathrm{\σ}}_b+M}---\left(6\right)$

${\left(s_m\right)}_{\mathrm{eq}}=\frac{{\mathrm{\σ}}_b(s_{m2}-s_{a2}-M)+2M(s_{m2}-s_{a2})}{{\mathrm{\σ}}_b+M}---\left(7\right)$

式中D表示当量损伤。

3、优点及功效：本发明一种编制疲劳加速试验载荷谱的方法，其特点是简单实用、操作方便、损伤等效精度高、时间压缩效果显著，

附图说明

图1为是本发明所述方法的流程框图。

具体实施方式

图1为本发明所述方法的流程框图，本发明分三步实现，具体为：

步骤一、疲劳载荷循环的判别

从工程观点出发，通常取N_f＝10⁷次循环和N_T＝10⁴次循环分别作为三级波和主波的判别门槛值，即裂纹形成寿命长于N_f所对应的应力循环可判别为三级波，而裂纹形成寿命短于N_T的应力循环则判别为主波。根据S-N曲面公式

${(\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-S_m}{S}_a-S_0)}^{m}N=C---\left(8\right)$

可以得出：

${\mathrm{\σ}}_b{s}_a+[{\left(\frac{C}{N_f}\right)}^{\frac{1}{m}}+S_0]s_m\≤[{\left(\frac{C}{N_f}\right)}^{\frac{1}{m}}+S_0]{\mathrm{\σ}}_b---\left(9\right)$

${\mathrm{\σ}}_b{s}_a+[{\left(\frac{C}{N_T}\right)}^{\frac{1}{m}}+S_0]s_m\≥[{\left(\frac{C}{N_T}\right)}^{\frac{1}{m}}+S_0]{\mathrm{\σ}}_b---\left(10\right)$

式中符号说明如下：σ_b为材料强度极限，s_a为应力幅值，s_m为应力均值，S₀为材料拟合疲劳极限，N为疲劳寿命，C和m为常数。

当应力循环满足式(9)时，则判别为三级波；而当应力循环满足式(10)时，则判别为主波；当式(9)和式(10)都不满足时，应力循环判别为二级波。

步骤二、低载截除

保留主波应力循环不变，并截除三级波应力循环。

步骤三、归并二级波，直至无法再次合并。

根据Miner线性累积损伤理论作为损伤当量折算的基础，一个应力循环(s_a,s_m)造成的损伤可写成

$D(s_a,s_m)=\frac{1}{N(s_a,s_m)}---\left(11\right)$

式中，N(s_a,s_m)由式(8)确定。为了将两个相邻的连续应力循环(s_a1,s_m1)和(s_a2,s_m2)归并为新二级波，根据损伤当量原则，由式(11)和式(8)可建立疲劳裂纹形成的损伤当量折算公式：

${[\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-{\left(s_m\right)}_{\mathrm{eq}}}{\left(s_a\right)}_{\mathrm{eq}}-S_0]}^{-m}={[\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m1}}{s}_{a1}-S_0]}^{-m}+{[\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m2}}{s}_{a2}-S_0]}^{-m}---\left(12\right)$

令(s_min)_eq＝s_min2＝s_m2-s_a2， $M={[{(\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m1}}{s}_{a1}-S_0)}^{-m}+{(\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m2}}{s}_{a2}-S_0)}^{-m}]}^{-\frac{1}{m}}+S_0,$ 则由式(12)可分别得到归并后新二级波的应力幅值和均值分别为：

${\left(s_a\right)}_{\mathrm{eq}}=\frac{({\mathrm{\σ}}_b+s_{a2}-s_{m2})M}{{\mathrm{\σ}}_b+M}---\left(13\right)$

${\left(s_m\right)}_{\mathrm{eq}}=\frac{{\mathrm{\σ}}_b(s_{m2}-s_{a2}-M)+2M(s_{m2}-s_{a2})}{{\mathrm{\σ}}_b+M}---\left(14\right)$

式中D为当量损伤。

Claims (1)

1.一种编制疲劳加速试验载荷谱的方法，其特征在于：该方法具体步骤如下：

步骤一、疲劳载荷循环的判别

从工程观点出发，通常取N_f＝10⁷次循环和N_T＝10⁴次循环对应的疲劳应力循环的水平分别作为三级波和主波的判别门槛值，即裂纹形成寿命长于N_f所对应的应力循环判别为三级波，而裂纹形成寿命短于N_T的应力循环则判别为主波；根据S-N曲面公式

${(\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-S_m}{S}_a-S_0)}^{m}N=C---\left(1\right)$

得出

${\mathrm{\σ}}_b{s}_a+[{\left(\frac{C}{N_f}\right)}^{\frac{1}{m}}+S_0]s_m\≤[{\left(\frac{C}{N_f}\right)}^{\frac{1}{m}}+S_0]{\mathrm{\σ}}_b---\left(2\right)$

${\mathrm{\σ}}_b{s}_a+[{\left(\frac{C}{N_T}\right)}^{\frac{1}{m}}+S_0]s_m\≥[{\left(\frac{C}{N_T}\right)}^{\frac{1}{m}}+S_0]{\mathrm{\σ}}_b---\left(3\right)$

式中符号说明如下：σ_b为材料强度极限，s_a为应力幅值，s_m为应力均值，S₀为材料拟合疲劳极限，N为疲劳寿命，C和m为常数；

当应力循环满足式(2)时，则判别为三级波；而当应力循环满足式(3)时，则判别为主波；当式(2)和式(3)都不满足时，应力循环判别为二级波；

步骤二、低载截除

保留主波应力循环不变，并截除三级波应力循环；

步骤三、归并二级波，直至无法再次合并；

根据Miner线性累积损伤理论作为损伤当量折算的基础，一个应力循环(s_a,s_m)造成的损伤写成

$D(s_a,s_m)=\frac{1}{N(s_a,s_m)}---\left(4\right)$

式中，N(s_a,s_m)由式(1)确定，为了将两个相邻的连续应力循环(s_a1,s_m1)和(s_a2,s_m2)归并为新二级波，根据损伤当量原则，由式(4)和式(1)建立疲劳裂纹形成的损伤当量折算公式：

${[\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-{\left(s_m\right)}_{\mathrm{eq}}}{\left(s_a\right)}_{\mathrm{eq}}-S_0]}^{-m}={[\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m1}}{s}_{a1}-S_0]}^{-m}+{[\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m2}}{s}_{a2}-S_0]}^{-m}---\left(5\right)$

令(s_min)_eq＝s_min2＝s_m2-s_a2， $M={[{(\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m1}}{s}_{a1}-S_0)}^{-m}+{(\frac{{\mathrm{\σ}}_b}{{\mathrm{\σ}}_b-s_{m2}}{s}_{a2}-S_0)}^{-m}]}^{-\frac{1}{m}}+S_0,$ 则由式(5)分别得到归并后新二级波的应力幅值和均值分别为：

${\left(s_a\right)}_{\mathrm{eq}}=\frac{({\mathrm{\σ}}_b+s_{a2}-s_{m2})M}{{\mathrm{\σ}}_b+M}---\left(6\right)$

${\left(s_m\right)}_{\mathrm{eq}}=\frac{{\mathrm{\σ}}_b(s_{m2}-s_{a2}-M)+2M(s_{m2}-s_{a2})}{{\mathrm{\σ}}_b+M}---\left(7\right)$

式中D表示当量损伤。