随机响应疲劳分析
随机载荷下结构的疲劳寿命研究
随机响应分析的功率谱密度 (PSD) 结果用于计算力矩 mn ,这些力矩用于生成周期数与应力范围的概率密度函数。
根据随机响应分析中产生的 PSD 应力计算 PSD 力矩。
功率谱密度 (PSD) 力矩

力矩的计算依据为:
- fk
- 频率值。
- Gk
- 特定频率下的 PSD 响应值 fk 。
计算应力范围出现概率
针对每个循环段内初始和最终应力范围值之间的应力范围出现概率的计算是用户定义的。
P(ΔSi) 在 (ΔSi − δS/2) 和(ΔSi + δS/2) 之间的应力范围出现的概率为:
概率密度函数(周期数与应力范围的概率密度)
以如上方式计算的 PSD 力矩用于应力范围概率密度函数 f(mn) 的生成。此函数基于指定的损伤模型。DIRLIK、LALANNE、NARROW 和 Steinberg 三波段 (THREE) 选项用于定义损伤模型。
- DIRLIK
- DIRLIK 提出了用于确定概率密度函数的封闭形式的解,即:p(S) = D1Qe−ZQ + D2ZR2e−Z22R2 + D3Ze−Z222√m0
- LALANNE
- LALANNE 随机疲劳损伤模型将概率密度函数描述为:p(S) = 1√m0 √1 − γ2√2πe−S28m0(1 − γ2) + Sγ4√m0(1 + erf(Sγ2√2m0(1 − γ2)))
- NARROW
- 窄带随机疲劳损伤模型使用以下概率函数:p(S) = (S4m0e−(S28m0))
- THREE
- Steinberg 三波段随机疲劳损伤模型使用以下概率函数:P(S) = {0.683 at 2√m00.271 at 4√m00.043 at 6√m0
上曲面应力范围因子
应力范围上限的计算方式为:
应力范围上限 = 2*RMS 应力*上曲面应力范围因子
RMS 应力是由随机响应子工况输出的。目标应力范围受限于应力范围的上限。在疲劳损伤计算中,将不考虑任何超出上限的应力。
上曲面应力范围
直接指定上曲面应力范围。
循环次数
计算应力范围的宽度 DS = δS ,并据此计算概率(参见Figure 2)。默认值为 100,第一个循环从 0.0 开始,至δS 。应力范围宽度计算为:DS=上曲面应力范围/循环次数
应力范围宽度
直接定义应力范围的宽度 ( δS )。
计算应力范围出现概率
- DIRLIK, LALANNE, NARROW
-
P(Si)
在
(ΔSi − δS/2)
和
(ΔSi + δS/2)
之间的应力范围出现的概率为: P(Si) = pi(Si)δS
选择损伤模型
- 应力的 PSD 力矩用于计算相应的力矩,这些力矩用于确定应力范围的概率密度函数。
- DIRLIK 和 LALANNE 模型在更广泛的应力范围谱系分布中产生概率。因此,当输入的随机信号由多个频率的各种应力范围组成时,应使用这些模型。如果使用 DIRLIK 和 LALANNE,概率密度函数中的信息能更好地捕捉到应力范围分布的更广范围。
- NARROW 模型适用于应力范围预计与特定应力范围分布的高概率密切相关的随机信号。因此,如果您知道输入的随机数据没有广泛的应力范围分布,而是主要集中在某个特定的应力范围,则应该选择 NARROW。这个模型预计最高的应力范围概率位于这个特定的应力范围或其周围。
- THREE 模型与 NARROW 模型相同,只是它预计随机信号的分布除了包含与 1*RMS 的关联外,还包含与 2*RMS 和 3*RMS 的关联(尽管关联较小)。因此,如果您的输入随机数据主要集中在 1*RMS 的应力范围内,其次是 2*RMS 和 3*RMS,那么您应该选择 THREE。
峰值数和零交叉数
- NARROW, THREE
- 原始时域随机载荷(基于频率的随机 PSD 载荷由其产生)中的每秒零交叉数确定为:nzcross = √m2m0其中,
- mn
- 计算出的相应力矩。
- DIRLIK, LALANNE
- 原始时域随机载荷(基于频率的随机 PSD 载荷由其产生)中的每秒峰值数确定为:npeaks = √m4m2其中,
- mn
- 计算出的相应力矩。
周期数
- 窄带 (NARROW),三波段 (THREE)
- 总周期数计算为:NT = nzcrossT其中,
- T
- 总曝光时间。
- DIRLIK, LALANNE
- 总周期数计算为:NT = npeaksT其中,
- T
- 总曝光时间。
- 特定应力范围的总周期数
- 应力范围
ΔSi
的总周期数计算方式为:Ni = P(ΔSi)NT
疲劳寿命和损伤
疲劳寿命(材料在失效前的特定应力范围 Si 的最大周期数)根据材料 SN 曲线计算为:
作为施加的随机载荷的结果,总疲劳损伤计算为:
为了将任何导致平均应力不为零的载荷的平均应力修正考虑在内,您可以定义一个由这种载荷(通常是重力载荷)组成的静态子工况。这个静态子工况可以在随机疲劳分析事件设置中引用。