一种可进行疲劳裂纹扩展速率估算和寿命预测的方法
【专利摘要】本发明公开了一种可进行疲劳裂纹扩展速率估算和寿命预测的方法，步骤：(1)进行疲劳裂纹试验，记录裂纹长度、时间等数据；(2)处理数据估计裂纹长度对应的扩展速率；(3)对Paris公式作最小二乘线合确定金属材料参数；(4)建立裂纹扩展摄动级数模型，得到控制方程与边界条件；(5)选择摄动参数渐进展开控制方程，形成摄动参数渐进序列；(6)合并同类项得到不同幂次下微分方程，确定待定系数；对微分方程迭代求解，得下一次计算时的边界条件；(7)利用计算结果，完成裂纹扩展速率估算，判断裂纹长度是否超临界裂纹，若没超，则转到(5)，若超过，则完成疲劳寿命预测。本发明可以有效估算裂纹扩展速率和预测疲劳寿命，提高计算精度。
【专利说明】
-种可进行疲劳裂纹扩展速率估算和寿命预测的方法
技术领域
[0001] 本发明设及疲劳裂纹扩展速率估算和寿命预测技术方法的研究，特别设及考虑疲 劳裂纹扩展模型中裂纹长度初始扰动量的摄动级数刻画表征W及基于摄动级数法的裂纹 扩展速率估算和寿命预测，W提高裂纹扩展速率估算和寿命预测的计算精度，保证计算结 果的正确性和有效性。
【背景技术】
[0002] 疲劳与断裂是工程中最常见的构件失效原因，在力学、航空、船舶等领域引起了广 泛关注。在高校及研究院所的研究者们大多关注含裂纹结构的特性，包括疲劳裂纹扩展、剩 余强度和剩余寿命。结构出现疲劳裂纹后，其在交变载荷作用下剩余寿命主要由裂纹扩展 速率所控制。因此，对带有初始裂纹构件的裂纹扩展速率和寿命预测的研究具有重要意义。
[0003] 寿命预测的研究，归根结底就是对裂纹扩展速率的研究。断裂力学是研究具有初 始缺陷的材料和结构在各种环境下裂纹的扩展、失稳和止裂的规律，W裂纹尺寸大小和裂 纹的扩展速率为结构损伤的判据，并用来估算疲劳裂纹的扩展寿命。在已有的研究疲劳裂 纹扩展速率的断裂力学的理论方法中，Paris疲劳裂纹扩展公式是研究裂纹扩展最基本的 也是被人们多年来普遍采用的公式，它建立了应力强度因子和裂纹扩展速率之间的关系， 是当今工程应用中预测疲劳裂纹扩展寿命理论的基础。然而需要指出的是，大量研究和实 验结果表明，材料疲劳特性具有很大的分散性，而且对载荷及环境、材料及结构、加工工艺 等多方面相当敏感。也就是说，针对具体工程问题，结构不可避免地存在一些初始缺陷或瑕 疵，导致裂纹扩展公式中的参数存在扰动，特别针对航空航天结构，工作环境和载荷日益复 杂严苛。因此，需要提出或发展一种全新的数值计算方法来考量运些因素或扰动量，进而给 出更为合理精确的疲劳预估寿命。在运方面，相关研究已有一些。然而需要强调的是：与 Paris疲劳裂纹扩展公式相比，运些修正的或发展的模型或方法形式过于复杂，且要求设计 人员具备熟悉精准的疲劳断裂知识。因此，如何利用普遍使用的化ris公式，考量初始扰动 参数的本质存在，在同时保证计算结果的精确有效性W及计算过程的方便适用性，对疲劳 裂纹扩展寿命进行合理评估，是一个具有重要意义和价值的研究课题。
[0004] 为了保留化ris公式的方便适用性，同时考虑疲劳裂纹参数扰动量的存在，利用摄 动级数法对其进行刻画标征，进而提出基于化ris公式和摄动级数法的疲劳裂纹扩展速率 评估和寿命预测公式，是一种实现结构疲劳裂纹寿命预测的一种即简便又有效的方法。
【发明内容】
[0005] 本发明解决的技术问题是:保持普遍使用的化ris疲劳裂纹扩展公式计算方便可 行的特点，同时考虑结构疲劳裂纹长度参数中不可避免的初始扰动量，将Paris公式和摄动 级数法有效结合起来，针对航空航天工程结构的疲劳裂纹扩展及剩余寿命的精确预测问 题，提供一种可有效评估结构疲劳裂纹扩展速率和预测疲劳寿命方法。
[0006] 本发明充分考虑实际工程结构疲劳裂纹长度的初始扰动量，W断裂力学中广泛使 用的化ris裂纹扩展评估公式为基础，引入摄动级数对初始扰动量进行刻画表征，推导出疲 劳裂纹扩展的摄动级数公式，进而对疲劳寿命进行预估。所得到的评估结果不仅可W达到 一定程度的精度，而且计算方便，便于研究人员理解和接受，且由于考虑了初始扰动的存 在，结果更加合理可信。本发明采用的技术方案实现步骤如下：
[0007] 第一步:按照金属材料疲劳裂纹扩展速率试验方法，采用INSTR0N8803-250kN疲劳 试验机，用标准试件在应力变化幅A o = const、应力比R、预制裂纹长度a日的条件下进行疲 劳裂纹扩展试验。标准试件裂纹每经历一定的时间或载荷循环数后，记录裂纹扩展长度ai W及加载时间tiW及循环次数Ni的试验数据，直到标准试件拉断或裂纹超出视野范围；其中 曰1表示第i次记录的裂纹长度、ti表示第i次记录的加载时间，Ni表示第i次记录的循环次数。 读数装置是带精密光栅与步进电机的读数显微镜，通过显示器确定裂纹尖端的位置，并由 坐标值计算裂纹长度。
[0008] 第二步:基于第一步得到的试验数据，包括扩展裂纹长度曰1、加载时间ti及循环次 数Ni,利用常用的疲劳裂纹扩展速率数据处理方法，对试验数据(ai，ti)进行处理，估计裂纹 扩展速率(daAlt)i。如用割线法，W二相邻数据点割线的斜率，作为二相邻数据点平均裂纹 尺寸S =(〇,.+,+ A)/ 2所对应的裂纹扩展速率，故有：
[0009] (da/dt)i= (ai+i-ai)/(ti+i_ti)
[0010]式中（da/dt)i即为裂纹从ai扩展到aw的平均速度。
[0011] 第S步:基于常用的疲劳裂纹扩展化ri S公式da/dt = cab，对其两边取对数运算， 公式重写为：
[0012] lg(da/dt) = lgc+blga
[0013] 其中a(t)表明裂纹长度a为与时间t相关的函数;C与b为金属材料常参数。利用第 二步估计的裂纹扩速率数据[(daAlt)i，ai]，对上式进行最小二乘法线性拟合，利用微积分 极值理论和克拉默法则，确定金属材料参数C与b。
[0014] 第四步:基于第S步得到的疲劳裂纹扩展化ris公式da(t)/dt = c{a(t)}b，引入摄 动级数刻画表征裂纹长度的扰动量A a(t)，即裂纹长度a(t)可表示为：
[0015]
[0016] 具甲e刃小撼功寥数，
a日UhaiW，a2W，…为摄动参数相应幕次下的待定系数， 建立疲劳裂纹扩展的摄动级数数学模型，略去高阶项，得到摄动级数数学模型的控制方程 为：
[0017]
[001引边界条件为：
[0019] a(0) = ao(0)+ai(0)e+a2(0)e^+---
[0020] a(tc) =ac = ao(tc)+ai(tc)e+a2(tc)e2+...
[0021] 其中ao(0)为初始裂纹长度的平均值或名义值，ai(0)为初始扰动量；曰2(0)，曰3 (0 )，…一般取0; a。为临界裂纹长度，t。为最终加载时间。
[0022] 第五步:选择小摄动参数e，利用化ylor级数展开和多项式定理，将第四步得到的 摄动级数数学模型的控制方程右式dtX.WE'丫'渐进展开，形成一个关于摄动参数的渐进 !-亡a 序列，即：
[0023]
[0024] 其中渐进序列中摄动参数的幕次Ei前面的系数待定；h，12,…，Ik为非负整数，分 别为多项式定理展开时ai，曰2，…，ak相对应的幕次。
[0025] 第六步:针对第五步得到的控制方程的渐进序列，合并摄动参数e的相同幕次项， 得到摄动参数不同幕次下的微分方程，即：
[0026；
[0027；
[0028；
[0029；
[0030；
[0031] 其中i表示摄动参数e的幕次为i，为非负整数。并利用第四步给出的初始边界条件 确定渐进序列中的每一项待定系数;利用确定的每一项系数，对微分方程迭代求解，所得结 果为下一次计算时的边界条件。
[0032] 第屯步:利用第六步得到的微分方程迭代求解的计算结果，完成相应加载时间下 的裂纹扩展长度和扩展速率估算，判断当前时刻裂纹长度a(t)是否超过临界裂纹长度a。， 若没有超过临界裂纹长度a。，则转到(5)继续求解，完成了当前时刻的疲劳裂纹扩展速率估 算，若超过临界裂纹长度a。，则终止计算，完成了疲劳裂纹寿命预测。
[0033] 所述步骤一中的一定载荷循环次数为500~1000次。
[0034] 所述步骤一中的一定时间为0.5~1分钟。
[0035] 所述步骤二中常用的疲劳裂纹扩展速率数据处理方法包括割线法W及修正割线 法。
[0036] 所述步骤S中常用的疲劳裂纹扩展公式化ris公式的更通用形式为：
[0037]
[0038] 其中C和m为金属材料常参数，N为循环次数，AK为应力强度因子变化幅；由应力强 度因子呆 I":
[0039]
[0040] 其中F为形状因子，A O为应力变化幅;则得：
[0041]
[0042] 最终得金属材料常参曼
[0043] 所述步骤四中引入的小摄动参数e为标量，在工程实际中取为1;初始扰动量ai(0) 来源于预制裂纹加工误差或测量装置的分辨率。
[0044] 所述步骤五中所述的化ylor级数展开建立在金属材料参数b不为整数的前提下； 针对特殊情况b为整数，仅需多项式定理对控制方程渐进展开。
[0045] 所述步骤六中的下一次计算时的边界条件指下一次计算时相应时刻的初始裂纹 长度及其扰动量。
[0046] 所沐击驢木中的临界裂纹长度ac由下式确定：
[0047]
[0048] 其中Kc为金属材料的断裂强度，F为形状因子，Omax为最大加载应力。
[0049] 本发明与现有技术相比的优点在于:本发明针对工程实际中含初始裂纹长度扰动 的结构提供了一种可考量初始扰动的结构疲劳裂纹扩展速率估算和寿命预测的新思路，在 保持化ris裂纹扩展公式方便实用的前提下，利用摄动级数对扰动进行合理地刻画表征，将 Paris公式和摄动级数方法二者的优势有效结合起来。所构建的疲劳裂纹扩展速率估算和 寿命预测的摄动级数方法，不仅保证了计算过程与化ris公式同样方便简单，且考虑了初始 扰动对裂纹扩展和寿命预测的影响。在对初始扰动存在的结构进行疲劳裂纹寿命估算时， 可W充分考虑裂纹长度初始扰动的影响，在确保结构疲劳寿命预测过程简单实用的前提下 可大大提高计算精度和可信性。
【附图说明】
[0050] 图1是本发明针对可进行疲劳裂纹扩展速率估算和寿命预测方法的流程图；
[0051] 图2是本发明中的单边裂纹侣合金板的几何模型示意图；
[0052] 图3是本发明中的单元裂纹侣合金板的有限云模型示意图；
[0053] 图4是本发明针对单元裂纹侣合金板的几种不同方法估算的时间-裂纹长度(t-a) 曲线比较示意图。
【具体实施方式】
[0054] 如图1所示，本发明提出了一种可进行疲劳裂纹扩展速率估算和寿命预测方法，包 括W下步骤：
[0055] (1)按照金属材料疲劳裂纹扩展速率试验方法，用标准试件在应力变化幅A曰= const、应力比R、预制裂纹长度ao的条件下进行疲劳裂纹扩展试验。标准试件裂纹每经历一 定的时间或载荷循环数后，记录裂纹长度aiW及加载时间tiW及循环次数Ni的试验数据，直 到标准试件拉断为止;其中ai表示第i次记录的裂纹长度、ti表示第i次记录的加载时间，Ni 表示第i次记录的循环次数。读数装置是带精密光栅与步进电机的读数显微镜，通过显示器 确定裂纹尖端的位置，并由坐标值计算裂纹长度。
[0056] (2)基于第一步得到的试验数据，包括扩展裂纹长度曰1、加载时间ti及循环次数Ni, 利用常用的疲劳裂纹扩展速率数据处理方法，包括割线法、修正割线法对试验数据(ai，ti) 进行处理，估计裂纹扩展速率(da/dt)i。如用割线法，W二相邻数据点割线的斜率，作为二 相邻数据点平均裂纹尺寸巧=(〇,,,+3,)/2所对应的裂纹扩展速率，故有：
[0057] (da/dt)i= (ai+i-ai)/(ti+i_ti)
[005引式中（daAlt)i即为裂纹从日1扩展到aw的平均速度。
[0059] 而修正割线法是指任意一个实验数据点（ai，ti)所对应的裂纹扩展速率(da/dt)i 为该数据点上下两段割线斜率的平均值，即有：
[0060]
[0061] 其中曰1-1，曰1，曰1+汾别为第i-1次、i次、i+1次所记录的裂纹长度;分别为 第i-1次、i次、i+1次所记录的加载时间。
[0062] (3)基于常用的疲劳裂纹扩展化ris公式da/dt = cab，对其两边取对数运算，公式 重写为：
[0063] lg(da/dt) = lgc+blga
[0064] 利用第二步估计的裂纹扩速率数据[(daAlt)i，ai]，对上式进行最小二乘法线性拟 合，确定金属材料参数C与b。具体操作过程为：
[0065] 根据微积分极值理论，应有拟合曲线满足：
[0066]
[0067]
[006引
[0069] 法则求解可得：
[0070]
[0071]
[0072] 进而可得到金属材料参数C和b。
[0073] (4)基于第S步得到的疲劳裂纹扩展化ris公式da(t)/dt = c{a(t)}b，引入摄动级 数刻画表征裂纹长度的扰动量Aa(t)，即裂纹长度a(t)可表示为：
[0074]
[0075] 其中e为小摄动参数，a日(*)，曰1(〇，曰2(*)，-'为摄动参数相应幕次下的待定系数；
[0076] 讲而姑欣的微分化古九，
[0077]
[0078] 建立疲劳裂纹扩展的摄动级数数学模型，略去高阶项，得到摄动级数数学模型的 控制方程为：
[0079
[0080]边界条件为：
[0081 ] a(0) = ao(0)+ai(0)e+a2(0)e^+---
[0082] a(tc) =ac = ao(tc)+ai(tc) e+a2(tc)e^+---
[0083] 其中ao(0)为初始裂纹长度的平均值或名义值，ai(0)为初始扰动量；曰2(0)，曰3 (0 )，…一般取0; a。为临界裂纹长度，t。为最终加载时间。
[0084] (5)选择小摄动参数e，利用化ylor级数展开和多项式定理，将第四步得到的摄动 级数数学模型的控制方程右；
^进展开，形成一个关于摄动参数的渐进序列。 由化ylor级数展开，略去高阶项，有：
[0085；
[0086；
[0087；
[0088] 其中n为正整数，满足11《4;11，12，-，，11<为非负整数，分别为多项式定理展开时曰1， 曰2，…，ak相对应的幕次。最终控制方程右式渐进序列得：
[0089]
[0090] 其中渐进序列中摄动参数的幕次Ei前面的系数待定。
[0091] (6)针对第五步得到的控制方程的渐进序列，合并摄动参数e的相同幕次项，得到 摄动参数不同幕次下的微分方程，即：
[0092]
[0093]
[0094]
[0095]
[0096]
[0097] 其中i表示摄动参数e的幕次为i，为非负整数。利用第四步给出的初始边界条件： ao(0)为初始裂纹名义值，ai (0)为初始扰动量，曰2 (0)，曰3 (0)，…均为0，确定渐进序列中的每 一项待定系数;利用确定的每一项系数，对微分方程迭代求解，所得结果为下一次计算时 的边界条件。且体讨超为：针对摄动参猶e的一次巧，有：
[009引
[0099] 得；
[0100] ao(t+ A t) = ao(t)+c{ao(t) }b A t
[0101] 针对摄动参数e的二次项，有：
[0102]
[0103] 得；
[0104] ai(t+A t) =ai(t)+cb{ao(t)}b-iai(t) A t[0105] 依次类推，迭代求解，可得a2(t+At), a3(t+At),…；最终得时刻t+At时的裂纹扩 展长度*?
[0106]
[0107] (7)利用第六步得到的微分方程迭代求解的计算结果，完成相应加载时间下的裂
纹扩展长度和扩展速率估心刻裂纹长度a(t)是否超过临界裂纹长度a。，其中 临界裂纹长度ac可由式子 确定，其中Kc为金属材料的断裂强度，F为形状因 子，Omax为最大加载应力。若没有超过临界裂纹长度a。，则转到（5)继续求解，完成了当前时 刻的疲劳裂纹扩展速率估算，若超过临界裂纹长度a。，则终止计算，终止时间为t。，完成了疲 劳裂纹寿命N。预测。疲劳裂纹寿命N。为：
[010 引 Nc = tcXf
[0109] 其中f为加载频率。
[0110] 实施例；
[0111] 为了更充分地了解该发明的特点及其对工程实际的适用性，本发明针对如图2-4 所示的受均布载荷作用的含单边裂纹的2024-T3侣合金矩形板进行疲劳裂纹扩展速率估算 与寿命预测。该侣合金矩形板长5 = 120mm，W=30mm。断裂强度A：, = 33M/^n/品，形状因子尸为 1.07，金属材料常参数C和m分别为3.868和0.377e-10。初始裂纹长度名义值ao为5mm，初始 扰动来源于测量装置分辨率。
[0112] 利用不考虑初始扰动的化ris公式W及大型有限元通用软件ANSYS求解得到的应 力强度因子，可W分别得到疲劳裂纹扩展的时间-裂纹长度(t-a)曲线；其中，大型有限元 通用软件ANSYS通过建立含单边裂纹的侣合金板的壳有限元模型，通过J积分的方法求解应 力强度因子。同样，利用提出的可进行疲劳裂纹扩展速率估算和寿命预测的摄动级数方法 对含单边侣合金板的寿命进行了预测，=种方法的预测结果如下表，相应的比较曲线图如 图4所示。
[0113]
[0114] 该实施例采用S种方法分别是：不考虑初始扰动的化ris公式、ANSYS有限元求解 W及考虑初始扰动的化ris公式的摄动级数法。从结果可W看出，ANSYS有限元求解方法的 结果与不考虑初始扰动的Paris公式相比，结果更接近于实验结果，但需要指出的是ANSYS 有限元求解过程复杂，时间消耗较大。而另一方面，考虑初始扰动的化ris公式的摄动级数 法，寿命预测结果更接近于实验结果，即精度较高，且保留了化ris公式计算方便可行的优 势。也就是说，本发明提出的摄动级数方法可有效预测结构的疲劳裂纹扩展寿命，工程适用 性较高。
[0115] 综上所述，本发明提出了一种可进行疲劳裂纹扩展速率估算和寿命预测方法。首 先，根据疲劳裂纹扩展试验得到的试验数据，利用割线法、修正割线法对裂纹扩展速率进行 评估，进而借助最小二乘法和克拉默法则，对化ris公式中的金属材料参数进行估计;其次， 结合摄动级数法，对裂纹长度的初始扰动进行合理刻画表征，引入小摄动参数，建立疲劳裂 纹扩展的摄动级数数学模型，并对相应的控制方程进行渐进序列展开;最后，略去高阶项， 合并摄动参数的相同幕次项，得到摄动参数不同幕次下的微分方程;借助边界条件，对微分 方程迭代求解，最终得到当前时刻的疲劳裂纹长度，进而与临界裂纹长度比较，完成疲劳裂 纹寿命预测。
[0116] W上仅是本发明的具体步骤，对本发明的保护范围不构成任何限制;凡采用等同 变换或者等效替换而形成的技术方案，均落在本发明权利保护范围之内。
[0117] 本发明未详细阐述部分属于本领域技术人员的公知技术。
【主权项】
1. 一种可进行疲劳裂纹扩展速率估算和寿命预测的方法，其特征在于实现步骤如下： 第一步：按照金属材料疲劳裂纹扩展速率试验方法，用标准试件在应力变化幅Λσ = const、应力比R、预制裂纹长度ao的条件下进行疲劳裂纹扩展试验，标准试件裂纹每经历一 定的时间或载荷循环数后，记录裂纹长度 &1以及加载时间t以及循环次数K的试验数据，直 到标准试件拉断为止;其中ai表示第i次记录的裂纹长度、^表示第i次记录的加载时间，化 表示第i次记录的循环次数； 第二步:基于第一步得到的试验数据，包括扩展裂纹长度ai、加载时间U及循环次数K， 利用常用的疲劳裂纹扩展速率数据处理方法，对扩展裂纹长度ai、加载时间^进行处理，估 计裂纹扩展速率(da/dth，采用割线法，以二相邻数据点割线的斜率，作为二相邻数据点平 均裂纹尺寸$ = (α,+1 + ?,) / 2所对应的裂纹扩展速率，有： (da/dt)i= (ai+i-ai)/(ti+i-ti) (1) 式中(da/dth即为裂纹从ai扩展到a1+1的平均速度； 第三步:基于常用的疲劳裂纹扩展Paris公式da(t)/dt = c{a(t)}b，对其两边取对数运 算为： lg(da/dt) = lgc+blga (2) 其中a(t)表明裂纹长度a为与时间t相关的函数;c与b为金属材料常参数，利用第二步 估计的裂纹扩速率数据[(da/dOuai]，对上(2)式进行最小二乘法线性拟合，确定金属材料 参数c与b; 第四步:基于第三步得到的疲劳裂纹扩展?31^8公式(^(1：)/(1丨=〇{3(1：)}15，引入摄动级 数刻画表征裂纹长度的扰动量Aa(t)，即裂纹长度a(t)表示为：其中ε为小摄动参数，&〇(〇，&1(〇，&2(〇，一为摄动参数相应幂次下的待定系数，建立 疲劳裂纹扩展的摄动级数数学模型，略去高阶项，得到摄动级数数学模型的控制方程为：边界条件为： a(0) = ao(0)+ai(0)e+a2(0)e2+··· a(tc) =ac = ao(tc)+ai(tc)e+a2(tc)e2+··· 其中ao (0)为初始裂纹长度的平均值或名义值，ai (0)为初始扰动量;a2 (0)，a3 (0)，…一 般取〇; a。为临界裂纹长度，t。为最终加载时间； 第五步:选择小摄动参数ε，利用Taylor级数展开和多项式定理，将第四步得到的摄动 级数数学模型的控制方程右式渐进展开，形成一个关于摄动参数的渐进序列， i-M 即：其中渐进序列中摄动参数的幂次ε1前面的系数待定;llh，···，lk为非负整数，分别为多 项式定理展开时ai，a2，…，ak相对应的幂次； 第六步:针对第五步得到的控制方程的渐进序列，合并摄动参数ε的相同幂次项，得到 摄动参数不同幂次下的微分方程，ΒΡ :其中i表示摄动参数ε的幂次为i，为非负整数。并利用第四步给出的初始边界条件确定 渐进序列中的每一项待定系数;利用确定的每一项系数，对微分方程迭代求解，所得结果为 下一次计算时的边界条件； 第七步:利用第六步得到的微分方程迭代求解的计算结果，完成相应加载时间下的裂 纹扩展长度和扩展速率估算，判断当前时刻裂纹长度a(t)是否超过临界裂纹长度a。，若没 有超过临界裂纹长度a。，则转到（5)继续求解，完成了当前时刻的疲劳裂纹扩展速率估算， 若超过临界裂纹长度a。，则终止计算，完成了疲劳裂纹寿命预测。2. 根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法，其特 征在于:所述步骤一中的一定载荷循环次数为500~1000次。3. 根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法，其特 征在于:所述步骤一中的一定时间为0.5~1分钟。4. 根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法，其特 征在于:所述步骤二中常用的疲劳裂纹扩展速率数据处理方法包括割线法以及修正割线 法。5. 根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法，其特 征在于:所述步骤三中常用的疲劳裂纹扩展公式Paris公式的更通用形式为：其中C和m为金属材料常参数，N为循环次数，ΛΚ为应力强度因子变化幅；由应力强度因 子为1··= >则有：其中F为形状因子，Λσ为应力变化幅;则得：最终得金属材料常参16. 根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法，其特 征在于:所述步骤四中引入的小摄动参数ε为标量，在工程实际中取为1;初始扰动量aKO) 来源于预制裂纹加工误差或测量装置的分辨率。7. 根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法，其特 征在于:所述步骤五中所述的Taylor级数展开建立在金属材料参数b不为整数的前提下;针 对特殊情况b为整数，仅需多项式定理对控制方程渐进展开。8. 根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法，其特 征在于:所述步骤六中的下一次计算时的边界条件指下一次计算时相应时刻的初始裂纹长 度及其扰动量。9. 根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法，其特 征在于:所述步骤七中的临界裂纹长度a。由下式确定：其中K。为金属材料的断裂强度，F为形状因子，〇max为最大加载应力。
【文档编号】G06F17/50GK105956315SQ201610326164
【公开日】2016年9月21日
【申请日】2016年5月17日
【发明人】邱志平, 陈潇, 王晓军, 王磊, 李云龙, 苏欢
【申请人】北京航空航天大学