[发明专利]一种多轴变幅加载下基于短裂纹扩展的全寿命预测方法

摘要

本发明公开了一种多轴变幅加载下基于短裂纹扩展的全寿命预测方法，该方法步骤为：在通过W‑B计数确定的每个半循环中选取最大剪应力、较大正应力所在平面为临界面，并利用临界面上的损伤参量来表征短裂纹扩展驱动力；基于拉伸型多轴疲劳损伤参量，建立适用于多轴变幅应力状态下的有效J积分计算公式；通过单轴恒幅拉压裂纹扩展试验中获得的恒幅下试件寿命和与之对应的裂纹长度数据，拟合出裂纹扩展常数c和m；利用S‑N曲线来拟合虚拟的裂纹初始尺寸a₀；利用上述的所提的J积分计算模型并基于Paris公式，可得到每个载荷块的裂纹扩展值，然后进行迭代循环累加，从而确定试件的疲劳寿命。本方法可以很好的描述多轴变幅应力加载对裂纹扩展速率的影响。

主权项

1.一种多轴变幅加载下基于短裂纹扩展的全寿命预测方法，其特征在于：具体步骤如下，步骤1)：薄壁管件在多轴变幅应力加载状态下，裂纹主要萌生于最大剪应力幅、较大正应力所在的平面；由于变幅情况下，每个循环的临界面各不相同，通过W‑B多轴循环计数确定多个半循环，根据每个半循环内最大剪应力幅值、较大正应力所在的平面，确定为每个半循环的临界面，并在该临界面上计算计算损伤参量来表征短裂纹扩展驱动力；步骤2)：以此临界面为基础，用Newman裂纹闭合公式来考虑裂纹闭合；利用上述临界面上的数据并配合下面的J积分公式计算相对应的有效J积分即为有效小裂纹扩展驱动力：其中，a为半裂纹长度，E为弹性模量，n′为循环硬化指数，Δσ_n，eff为考虑了闭合后的有效正应力，Δτ_n为剪应力幅值，Δε_pl.eff为Δσ_n，eff所对应的塑性应变，由Ramberg‑Osgood公式求出；步骤3)：通过单轴拉压试验中获得的恒幅下试件寿命和与之对应的裂纹长度数据，通过计算得到有效J积分与裂纹扩展速度之间的双对数曲线，并由此拟合出裂纹扩展常数c和m，并以此为基线进行下一步计算；Paris形式的裂纹扩展曲线公式如下：其中，是裂纹扩展速率，c，m为单轴Paris常数；步骤4)：确定裂纹初始尺寸，选用Michael Vormwald提出的用材料的S‑N曲线来拟合虚拟的裂纹初始尺寸来考虑材料缺陷因素；具体方法是在S‑N曲线中选取尽可能多的点(σ_i，0，N_i)，对每一个点用公式计算一次a_i，0，最后取所有a_i，0的算术平均值，即可得虚拟的裂纹初始尺寸，具体公式如下：其中，a_i，0是在S‑N曲线上所取第i个点所算的裂纹初始尺寸，a_f为最终失效尺寸，a₀是所有取的n个点所得初始尺寸的平均值；步骤5)：在临界面上进行分析，将多轴有效J积分与Paris公式相结合，便可以得到基于临界面法的短裂纹扩展模型，公式如下：经过Wang‑Brown多轴循环计数，得到该多轴变幅载荷块下多个半循环，计算每个半循环的裂纹扩展量公式如下：多轴变幅载荷块通过Wang‑Brown多轴循环计数，得到该载荷块下在临界面上多个半循环，每个半循环相应产生的裂纹扩展量就为Δa₁，Δa₂，Δa₃，......，Δa_k，一个载荷谱块作用下所产生的裂纹总的增加量则为：Δa＝Δa₁+Δa₂+Δa₃……+Δa_kk为该载荷块多轴计数后半循环的个数；得到了每个载荷块的裂纹扩展值计算方法，进行短裂纹的扩展量迭代循环累加，直到裂纹长度达到短裂纹的临界长度截止；这样便能计算得到试件的疲劳寿命。