CN105956315A - 一种可进行疲劳裂纹扩展速率估算和寿命预测的方法 - Google Patents
一种可进行疲劳裂纹扩展速率估算和寿命预测的方法 Download PDFInfo
- Publication number
- CN105956315A CN105956315A CN201610326164.2A CN201610326164A CN105956315A CN 105956315 A CN105956315 A CN 105956315A CN 201610326164 A CN201610326164 A CN 201610326164A CN 105956315 A CN105956315 A CN 105956315A
- Authority
- CN
- China
- Prior art keywords
- crack
- growth rate
- perturbation
- sigma
- fatigue crack
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种可进行疲劳裂纹扩展速率估算和寿命预测的方法,步骤:(1)进行疲劳裂纹试验,记录裂纹长度、时间等数据;(2)处理数据估计裂纹长度对应的扩展速率;(3)对Paris公式作最小二乘线合确定金属材料参数;(4)建立裂纹扩展摄动级数模型,得到控制方程与边界条件;(5)选择摄动参数渐进展开控制方程,形成摄动参数渐进序列;(6)合并同类项得到不同幂次下微分方程,确定待定系数;对微分方程迭代求解,得下一次计算时的边界条件;(7)利用计算结果,完成裂纹扩展速率估算,判断裂纹长度是否超临界裂纹,若没超,则转到(5),若超过,则完成疲劳寿命预测。本发明可以有效估算裂纹扩展速率和预测疲劳寿命,提高计算精度。
Description
技术领域
本发明涉及疲劳裂纹扩展速率估算和寿命预测技术方法的研究,特别涉及考虑疲劳裂纹扩展模型中裂纹长度初始扰动量的摄动级数刻画表征以及基于摄动级数法的裂纹扩展速率估算和寿命预测,以提高裂纹扩展速率估算和寿命预测的计算精度,保证计算结果的正确性和有效性。
背景技术
疲劳与断裂是工程中最常见的构件失效原因,在力学、航空、船舶等领域引起了广泛关注。在高校及研究院所的研究者们大多关注含裂纹结构的特性,包括疲劳裂纹扩展、剩余强度和剩余寿命。结构出现疲劳裂纹后,其在交变载荷作用下剩余寿命主要由裂纹扩展速率所控制。因此,对带有初始裂纹构件的裂纹扩展速率和寿命预测的研究具有重要意义。
寿命预测的研究,归根结底就是对裂纹扩展速率的研究。断裂力学是研究具有初始缺陷的材料和结构在各种环境下裂纹的扩展、失稳和止裂的规律,以裂纹尺寸大小和裂纹的扩展速率为结构损伤的判据,并用来估算疲劳裂纹的扩展寿命。在已有的研究疲劳裂纹扩展速率的断裂力学的理论方法中,Paris疲劳裂纹扩展公式是研究裂纹扩展最基本的也是被人们多年来普遍采用的公式,它建立了应力强度因子和裂纹扩展速率之间的关系,是当今工程应用中预测疲劳裂纹扩展寿命理论的基础。然而需要指出的是,大量研究和实验结果表明,材料疲劳特性具有很大的分散性,而且对载荷及环境、材料及结构、加工工艺等多方面相当敏感。也就是说,针对具体工程问题,结构不可避免地存在一些初始缺陷或瑕疵,导致裂纹扩展公式中的参数存在扰动,特别针对航空航天结构,工作环境和载荷日益复杂严苛。因此,需要提出或发展一种全新的数值计算方法来考量这些因素或扰动量,进而给出更为合理精确的疲劳预估寿命。在这方面,相关研究已有一些。然而需要强调的是:与Paris疲劳裂纹扩展公式相比,这些修正的或发展的模型或方法形式过于复杂,且要求设计人员具备熟悉精准的疲劳断裂知识。因此,如何利用普遍使用的Paris公式,考量初始扰动参数的本质存在,在同时保证计算结果的精确有效性以及计算过程的方便适用性,对疲劳裂纹扩展寿命进行合理评估,是一个具有重要意义和 价值的研究课题。
为了保留Paris公式的方便适用性,同时考虑疲劳裂纹参数扰动量的存在,利用摄动级数法对其进行刻画标征,进而提出基于Paris公式和摄动级数法的疲劳裂纹扩展速率评估和寿命预测公式,是一种实现结构疲劳裂纹寿命预测的一种即简便又有效的方法。
发明内容
本发明解决的技术问题是:保持普遍使用的Paris疲劳裂纹扩展公式计算方便可行的特点,同时考虑结构疲劳裂纹长度参数中不可避免的初始扰动量,将Paris公式和摄动级数法有效结合起来,针对航空航天工程结构的疲劳裂纹扩展及剩余寿命的精确预测问题,提供一种可有效评估结构疲劳裂纹扩展速率和预测疲劳寿命方法。
本发明充分考虑实际工程结构疲劳裂纹长度的初始扰动量,以断裂力学中广泛使用的Paris裂纹扩展评估公式为基础,引入摄动级数对初始扰动量进行刻画表征,推导出疲劳裂纹扩展的摄动级数公式,进而对疲劳寿命进行预估。所得到的评估结果不仅可以达到一定程度的精度,而且计算方便,便于研究人员理解和接受,且由于考虑了初始扰动的存在,结果更加合理可信。本发明采用的技术方案实现步骤如下:
第一步:按照金属材料疲劳裂纹扩展速率试验方法,采用INSTRON8803-250kN疲劳试验机,用标准试件在应力变化幅Δσ=const、应力比R、预制裂纹长度a0的条件下进行疲劳裂纹扩展试验。标准试件裂纹每经历一定的时间或载荷循环数后,记录裂纹扩展长度ai以及加载时间ti以及循环次数Ni的试验数据,直到标准试件拉断或裂纹超出视野范围;其中ai表示第i次记录的裂纹长度、ti表示第i次记录的加载时间,Ni表示第i次记录的循环次数。读数装置是带精密光栅与步进电机的读数显微镜,通过显示器确定裂纹尖端的位置,并由坐标值计算裂纹长度。
第二步:基于第一步得到的试验数据,包括扩展裂纹长度ai、加载时间ti及循环次数Ni,利用常用的疲劳裂纹扩展速率数据处理方法,对试验数据(ai,ti)进行处理,估计裂纹扩展速率(da/dt)i。如用割线法,以二相邻数据点割线的斜率,作为二相邻数据点平均裂纹尺寸所对应的裂纹扩展速率,故有:
(da/dt)i=(ai+1-ai)/(ti+1-ti)
式中(da/dt)i即为裂纹从ai扩展到ai+1的平均速度。
第三步:基于常用的疲劳裂纹扩展Paris公式da/dt=cab,对其两边取对数运算,公式重写为:
lg(da/dt)=lgc+blga
其中a(t)表明裂纹长度a为与时间t相关的函数;c与b为金属材料常参数。利用第二步估计的裂纹扩速率数据[(da/dt)i,ai],对上式进行最小二乘法线性拟合,利用微积分极值理论和克拉默法则,确定金属材料参数c与b。
第四步:基于第三步得到的疲劳裂纹扩展Paris公式da(t)/dt=c{a(t)}b,引入摄动级数刻画表征裂纹长度的扰动量Δa(t),即裂纹长度a(t)可表示为:
其中ε为小摄动参数,a0(t),a1(t),a2(t),…为摄动参数相应幂次下的待定系数,建立疲劳裂纹扩展的摄动级数数学模型,略去高阶项,得到摄动级数数学模型的控制方程为:
边界条件为:
a(0)=a0(0)+a1(0)ε+a2(0)ε2+…
a(tc)=ac=a0(tc)+a1(tc)ε+a2(tc)ε2+…
其中a0(0)为初始裂纹长度的平均值或名义值,a1(0)为初始扰动量;a2(0),a3(0),…一般取0;ac为临界裂纹长度,tc为最终加载时间。
第五步:选择小摄动参数ε,利用Taylor级数展开和多项式定理,将第四步得到的摄动级数数学模型的控制方程右式渐进展开,形成一个关于摄动参数的渐进序列,即:
其中渐进序列中摄动参数的幂次εi前面的系数待定;l1,l2,…,lk为非负整数,分别为多项式定理展开时a1,a2,…,ak相对应的幂次。
第六步:针对第五步得到的控制方程的渐进序列,合并摄动参数ε的相同幂次项,得到摄动参数不同幂次下的微分方程,即:
其中i表示摄动参数ε的幂次为i,为非负整数。并利用第四步给出的初始边界条件确定渐进序列中的每一项待定系数;利用确定的每一项系数,对微分方程迭代求解,所得结果为下一次计算时的边界条件。
第七步:利用第六步得到的微分方程迭代求解的计算结果,完成相应加载时间下的裂纹扩展长度和扩展速率估算,判断当前时刻裂纹长度a(t)是否超过临界裂纹长度ac,若没有超过临界裂纹长度ac,则转到(5)继续求解,完成了当前时刻的疲劳裂纹扩展速率估算,若超过临界裂纹长度ac,则终止计算,完成了疲劳裂纹寿命预测。
所述步骤一中的一定载荷循环次数为500~1000次。
所述步骤一中的一定时间为0.5~1分钟。
所述步骤二中常用的疲劳裂纹扩展速率数据处理方法包括割线法以及修正割线法。
所述步骤三中常用的疲劳裂纹扩展公式Paris公式的更通用形式为:
其中C和m为金属材料常参数,N为循环次数,ΔK为应力强度因子变化幅;由应力强度因子为则有:
其中F为形状因子,Δσ为应力变化幅;则得:
最终得金属材料常参数和b=m/2。
所述步骤四中引入的小摄动参数ε为标量,在工程实际中取为1;初始扰动量a1(0)来 源于预制裂纹加工误差或测量装置的分辨率。
所述步骤五中所述的Taylor级数展开建立在金属材料参数b不为整数的前提下;针对特殊情况b为整数,仅需多项式定理对控制方程渐进展开。
所述步骤六中的下一次计算时的边界条件指下一次计算时相应时刻的初始裂纹长度及其扰动量。
所述步骤七中的临界裂纹长度ac由下式确定:
其中Kc为金属材料的断裂强度,F为形状因子,σmax为最大加载应力。
本发明与现有技术相比的优点在于:本发明针对工程实际中含初始裂纹长度扰动的结构提供了一种可考量初始扰动的结构疲劳裂纹扩展速率估算和寿命预测的新思路,在保持Paris裂纹扩展公式方便实用的前提下,利用摄动级数对扰动进行合理地刻画表征,将Paris公式和摄动级数方法二者的优势有效结合起来。所构建的疲劳裂纹扩展速率估算和寿命预测的摄动级数方法,不仅保证了计算过程与Paris公式同样方便简单,且考虑了初始扰动对裂纹扩展和寿命预测的影响。在对初始扰动存在的结构进行疲劳裂纹寿命估算时,可以充分考虑裂纹长度初始扰动的影响,在确保结构疲劳寿命预测过程简单实用的前提下可大大提高计算精度和可信性。
附图说明
图1是本发明针对可进行疲劳裂纹扩展速率估算和寿命预测方法的流程图;
图2是本发明中的单边裂纹铝合金板的几何模型示意图;
图3是本发明中的单元裂纹铝合金板的有限云模型示意图;
图4是本发明针对单元裂纹铝合金板的几种不同方法估算的时间-裂纹长度(t-a)曲线比较示意图。
具体实施方式
如图1所示,本发明提出了一种可进行疲劳裂纹扩展速率估算和寿命预测方法,包括以下步骤:
(1)按照金属材料疲劳裂纹扩展速率试验方法,用标准试件在应力变化幅Δσ=const、应力比R、预制裂纹长度a0的条件下进行疲劳裂纹扩展试验。标准试件裂纹每经历一定的时间或载荷循环数后,记录裂纹长度ai以及加载时间ti以及循环次数Ni的试验数据,直到标准试件拉断为止;其中ai表示第i次记录的裂纹长度、ti表示第i次记 录的加载时间,Ni表示第i次记录的循环次数。读数装置是带精密光栅与步进电机的读数显微镜,通过显示器确定裂纹尖端的位置,并由坐标值计算裂纹长度。
(2)基于第一步得到的试验数据,包括扩展裂纹长度ai、加载时间ti及循环次数Ni,利用常用的疲劳裂纹扩展速率数据处理方法,包括割线法、修正割线法对试验数据(ai,ti)进行处理,估计裂纹扩展速率(da/dt)i。如用割线法,以二相邻数据点割线的斜率,作为二相邻数据点平均裂纹尺寸所对应的裂纹扩展速率,故有:
(da/dt)i=(ai+1-ai)/(ti+1-ti)
式中(da/dt)i即为裂纹从ai扩展到ai+1的平均速度。
而修正割线法是指任意一个实验数据点(ai,ti)所对应的裂纹扩展速率(da/dt)i为该数据点上下两段割线斜率的平均值,即有:
其中ai-1,ai,ai+1分别为第i-1次、i次、i+1次所记录的裂纹长度;ti-1,ti,ti+1分别为第i-1次、i次、i+1次所记录的加载时间。
(3)基于常用的疲劳裂纹扩展Paris公式da/dt=cab,对其两边取对数运算,公式重写为:
lg(da/dt)=lgc+blga
利用第二步估计的裂纹扩速率数据[(da/dt)i,ai],对上式进行最小二乘法线性拟合,确定金属材料参数c与b。具体操作过程为:
根据微积分极值理论,应有拟合曲线满足:
或
其中m为估计的裂纹扩展速率的个数。利用消元法或克拉默法则求解可得:
进而可得到金属材料参数c和b。
(4)基于第三步得到的疲劳裂纹扩展Paris公式da(t)/dt=c{a(t)}b,引入摄动级数刻画表征裂纹长度的扰动量Δa(t),即裂纹长度a(t)可表示为:
其中ε为小摄动参数,a0(t),a1(t),a2(t),…为摄动参数相应幂次下的待定系数;
进而相应的微分形式为:
建立疲劳裂纹扩展的摄动级数数学模型,略去高阶项,得到摄动级数数学模型的控制方程为:
边界条件为:
a(0)=a0(0)+a1(0)ε+a2(0)ε2+…
a(tc)=ac=a0(tc)+a1(tc)ε+a2(tc)ε2+…
其中a0(0)为初始裂纹长度的平均值或名义值,a1(0)为初始扰动量;a2(0),a3(0),…一般取0;ac为临界裂纹长度,tc为最终加载时间。
(5)选择小摄动参数ε,利用Taylor级数展开和多项式定理,将第四步得到的摄动级数数学模型的控制方程右式渐进展开,形成一个关于摄动参数的渐进序列。由Taylor级数展开,略去高阶项,有:
由多项式定理有:
其中n为正整数,满足n≤k;l1,l2,…,lk为非负整数,分别为多项式定理展开时a1,a2,…,ak相对应的幂次。最终控制方程右式渐进序列得:
其中渐进序列中摄动参数的幂次εi前面的系数待定。
(6)针对第五步得到的控制方程的渐进序列,合并摄动参数ε的相同幂次项,得到摄动参数不同幂次下的微分方程,即:
其中i表示摄动参数ε的幂次为i,为非负整数。利用第四步给出的初始边界条件:a0(0)为初始裂纹名义值,a1(0)为初始扰动量,a2(0),a3(0),…均为0,确定渐进序列中的每一 项待定系数;利用确定的每一项系数,对微分方程迭代求解,所得结果为下一次计算时的边界条件。具体过程为:针对摄动参数ε的一次项,有:
得:
a0(t+Δt)=a0(t)+c{a0(t)}bΔt
针对摄动参数ε的二次项,有:
得:
a1(t+Δt)=a1(t)+cb{a0(t)}b-1a1(t)Δt
依次类推,迭代求解,可得a2(t+Δt),a3(t+Δt),…;最终得时刻t+Δt时的裂纹扩展长度为:
(7)利用第六步得到的微分方程迭代求解的计算结果,完成相应加载时间下的裂纹扩展长度和扩展速率估算,判断当前时刻裂纹长度a(t)是否超过临界裂纹长度ac,其中临界裂纹长度ac可由式子确定,其中Kc为金属材料的断裂强度,F为形状因子,σmax为最大加载应力。若没有超过临界裂纹长度ac,则转到(5)继续求解,完成了当前时刻的疲劳裂纹扩展速率估算,若超过临界裂纹长度ac,则终止计算,终止时间为tc,完成了疲劳裂纹寿命Nc预测。疲劳裂纹寿命Nc为:
Nc=tc×f
其中f为加载频率。
实施例:
为了更充分地了解该发明的特点及其对工程实际的适用性,本发明针对如图2-4所示的受均布载荷作用的含单边裂纹的2024-T3铝合金矩形板进行疲劳裂纹扩展速率估算与寿命预测。该铝合金矩形板长S=120mm,W=30mm。断裂强度形状因子F为1.07,金属材料常参数C和m分别为3.868和0.377e-10。初始裂纹长度名义值a0为5mm,初始扰动来源于测量装置分辨率。
利用不考虑初始扰动的Paris公式以及大型有限元通用软件ANSYS求解得到的应力强度因子,可以分别得到疲劳裂纹扩展的时间-裂纹长度(t-a)曲线;其中,大型有限 元通用软件ANSYS通过建立含单边裂纹的铝合金板的壳有限元模型,通过J积分的方法求解应力强度因子。同样,利用提出的可进行疲劳裂纹扩展速率估算和寿命预测的摄动级数方法对含单边铝合金板的寿命进行了预测,三种方法的预测结果如下表,相应的比较曲线图如图4所示。
该实施例采用三种方法分别是:不考虑初始扰动的Paris公式、ANSYS有限元求解以及考虑初始扰动的Paris公式的摄动级数法。从结果可以看出,ANSYS有限元求解方法的结果与不考虑初始扰动的Paris公式相比,结果更接近于实验结果,但需要指出的是ANSYS有限元求解过程复杂,时间消耗较大。而另一方面,考虑初始扰动的Paris公式的摄动级数法,寿命预测结果更接近于实验结果,即精度较高,且保留了Paris公式计算方便可行的优势。也就是说,本发明提出的摄动级数方法可有效预测结构的疲劳裂纹扩展寿命,工程适用性较高。
综上所述,本发明提出了一种可进行疲劳裂纹扩展速率估算和寿命预测方法。首先,根据疲劳裂纹扩展试验得到的试验数据,利用割线法、修正割线法对裂纹扩展速率进行评估,进而借助最小二乘法和克拉默法则,对Paris公式中的金属材料参数进行估计;其次,结合摄动级数法,对裂纹长度的初始扰动进行合理刻画表征,引入小摄动参数,建立疲劳裂纹扩展的摄动级数数学模型,并对相应的控制方程进行渐进序列展开;最后,略去高阶项,合并摄动参数的相同幂次项,得到摄动参数不同幂次下的微分方程;借助边界条件,对微分方程迭代求解,最终得到当前时刻的疲劳裂纹长度,进而与临界裂纹长度比较,完成疲劳裂纹寿命预测。
以上仅是本发明的具体步骤,对本发明的保护范围不构成任何限制;凡采用等同变换或者等效替换而形成的技术方案,均落在本发明权利保护范围之内。
本发明未详细阐述部分属于本领域技术人员的公知技术。
Claims (9)
1.一种可进行疲劳裂纹扩展速率估算和寿命预测的方法,其特征在于实现步骤如下:
第一步:按照金属材料疲劳裂纹扩展速率试验方法,用标准试件在应力变化幅△σ=const、应力比R、预制裂纹长度a0的条件下进行疲劳裂纹扩展试验,标准试件裂纹每经历一定的时间或载荷循环数后,记录裂纹长度ai以及加载时间ti以及循环次数Ni的试验数据,直到标准试件拉断为止;其中ai表示第i次记录的裂纹长度、ti表示第i次记录的加载时间,Ni表示第i次记录的循环次数;
第二步:基于第一步得到的试验数据,包括扩展裂纹长度ai、加载时间ti及循环次数Ni,利用常用的疲劳裂纹扩展速率数据处理方法,对扩展裂纹长度ai、加载时间ti进行处理,估计裂纹扩展速率(da/dt)i,采用割线法,以二相邻数据点割线的斜率,作为二相邻数据点平均裂纹尺寸所对应的裂纹扩展速率,有:
(da/dt)i=(ai+1-ai)/(ti+1-ti) (1)
式中(da/dt)i即为裂纹从ai扩展到ai+1的平均速度;
第三步:基于常用的疲劳裂纹扩展Paris公式da(t)/dt=c{a(t)}b,对其两边取对数运算为:
lg(da/dt)=lgc+blga (2)
其中a(t)表明裂纹长度a为与时间t相关的函数;c与b为金属材料常参数,利用第二步估计的裂纹扩速率数据[(da/dt)i,ai],对上(2)式进行最小二乘法线性拟合,确定金属材料参数c与b;
第四步:基于第三步得到的疲劳裂纹扩展Paris公式da(t)/dt=c{a(t)}b,引入摄动级数刻画表征裂纹长度的扰动量△a(t),即裂纹长度a(t)表示为:
其中ε为小摄动参数,a0(t),a1(t),a2(t),…为摄动参数相应幂次下的待定系数,建立疲劳裂纹扩展的摄动级数数学模型,略去高阶项,得到摄动级数数学模型的控制方程为:
边界条件为:
a(0)=a0(0)+a1(0)ε+a2(0)ε2+…
a(tc)=ac=a0(tc)+a1(tc)ε+a2(tc)ε2+…
其中a0(0)为初始裂纹长度的平均值或名义值,a1(0)为初始扰动量;a2(0),a3(0),…一般取0;ac为临界裂纹长度,tc为最终加载时间;
第五步:选择小摄动参数ε,利用Taylor级数展开和多项式定理,将第四步得到的摄动级数数学模型的控制方程右式渐进展开,形成一个关于摄动参数的渐进序列,即:
其中渐进序列中摄动参数的幂次εi前面的系数待定;l1,l2,…,lk为非负整数,分别为多项式定理展开时a1,a2,…,ak相对应的幂次;
第六步:针对第五步得到的控制方程的渐进序列,合并摄动参数ε的相同幂次项,得到摄动参数不同幂次下的微分方程,即:
ε0:
ε1:
ε2:
ε3:
εi:
其中i表示摄动参数ε的幂次为i,为非负整数。并利用第四步给出的初始边界条件确定渐进序列中的每一项待定系数;利用确定的每一项系数,对微分方程迭代求解,所得结果为下一次计算时的边界条件;
第七步:利用第六步得到的微分方程迭代求解的计算结果,完成相应加载时间下的裂纹扩展长度和扩展速率估算,判断当前时刻裂纹长度a(t)是否超过临界裂纹长度ac,若没有超过临界裂纹长度ac,则转到(5)继续求解,完成了当前时刻的疲劳裂纹扩展速率估算,若超过临界裂纹长度ac,则终止计算,完成了疲劳裂纹寿命预测。
2.根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法,其特征在于:所述步骤一中的一定载荷循环次数为500~1000次。
3.根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法,其特征在于:所述步骤一中的一定时间为0.5~1分钟。
4.根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法,其特征在于:所述步骤二中常用的疲劳裂纹扩展速率数据处理方法包括割线法以及修正割线法。
5.根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法,其特征在于:所述步骤三中常用的疲劳裂纹扩展公式Paris公式的更通用形式为:
其中C和m为金属材料常参数,N为循环次数,△K为应力强度因子变化幅;由应力强度因子为则有:
其中F为形状因子,△σ为应力变化幅;则得:
最终得金属材料常参数和b=m/2。
6.根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法,其特征在于:所述步骤四中引入的小摄动参数ε为标量,在工程实际中取为1;初始扰动量a1(0)来源于预制裂纹加工误差或测量装置的分辨率。
7.根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法,其特征在于:所述步骤五中所述的Taylor级数展开建立在金属材料参数b不为整数的前提下;针对特殊情况b为整数,仅需多项式定理对控制方程渐进展开。
8.根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法,其特征在于:所述步骤六中的下一次计算时的边界条件指下一次计算时相应时刻的初始裂纹长度及其扰动量。
9.根据权利要求1所述的一种可进行疲劳裂纹扩展速率估算和寿命预测的方法,其特征在于:所述步骤七中的临界裂纹长度ac由下式确定:
其中Kc为金属材料的断裂强度,F为形状因子,σmax为最大加载应力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610326164.2A CN105956315B (zh) | 2016-05-17 | 2016-05-17 | 一种可进行疲劳裂纹扩展速率估算和寿命预测的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610326164.2A CN105956315B (zh) | 2016-05-17 | 2016-05-17 | 一种可进行疲劳裂纹扩展速率估算和寿命预测的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105956315A true CN105956315A (zh) | 2016-09-21 |
CN105956315B CN105956315B (zh) | 2018-03-09 |
Family
ID=56912696
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610326164.2A Expired - Fee Related CN105956315B (zh) | 2016-05-17 | 2016-05-17 | 一种可进行疲劳裂纹扩展速率估算和寿命预测的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105956315B (zh) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106525365A (zh) * | 2016-10-28 | 2017-03-22 | 中航动力股份有限公司 | 一种航空发动机轮盘裂纹予制试验系统及试验方法 |
CN106568660A (zh) * | 2016-10-14 | 2017-04-19 | 中国人民解放军海军航空工程学院青岛校区 | 一种复合材料胶接修补结构的剩余疲劳寿命的预测方法 |
CN107389471A (zh) * | 2017-07-13 | 2017-11-24 | 电子科技大学 | 裂纹失效模式判定方法及基于该方法的疲劳寿命预测方法 |
CN107449894A (zh) * | 2017-08-23 | 2017-12-08 | 中国人民解放军海军工程大学 | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 |
CN108334666A (zh) * | 2018-01-08 | 2018-07-27 | 中车青岛四方机车车辆股份有限公司 | 一种车轴抗疲劳断裂分析方法 |
CN108828069A (zh) * | 2018-06-19 | 2018-11-16 | 北京领示科技有限公司 | 一种基于超声定量数据的关键件疲劳寿命的定寿方法 |
CN109142529A (zh) * | 2018-08-27 | 2019-01-04 | 佛山科学技术学院 | 一种高强钛合金电子束焊接接头超高周疲劳寿命预测方法 |
CN109708966A (zh) * | 2018-12-04 | 2019-05-03 | 中国大唐集团科学技术研究院有限公司火力发电技术研究院 | 一种电站高温承压管道焊接接头多缺陷安全性评定方法 |
CN109918776A (zh) * | 2019-03-05 | 2019-06-21 | 北京航空航天大学 | 基于双步最小二乘法的机械产品的疲劳裂纹预测方法 |
CN110489848A (zh) * | 2019-08-12 | 2019-11-22 | 武汉钢铁有限公司 | 一种不同海水流速腐蚀疲劳裂纹扩展速率预测方法 |
CN110907270A (zh) * | 2019-11-29 | 2020-03-24 | 浙江大学宁波理工学院 | 利用铁磁性材料的弱磁信号预测疲劳寿命的方法 |
CN111044367A (zh) * | 2019-12-31 | 2020-04-21 | 成都理工大学 | 一种基于三轴应力-应变曲线的岩石裂缝亚临界扩展速率实验测试方法 |
CN111815694A (zh) * | 2020-09-14 | 2020-10-23 | 湖南大学 | 一种疲劳裂纹扩展寿命预测方法、装置、设备及存储介质 |
CN112380656A (zh) * | 2020-11-20 | 2021-02-19 | 西安热工研究院有限公司 | 一种燃气轮机燃烧室部件裂纹扩展寿命评估方法 |
CN112464572A (zh) * | 2020-12-10 | 2021-03-09 | 北京航空航天大学 | 一种传动轴加工性能评价及基于关键工艺的控制方法 |
CN112699479A (zh) * | 2020-12-29 | 2021-04-23 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机结构裂纹形成寿命理论值分散系数的确定方法 |
CN112765854A (zh) * | 2021-01-21 | 2021-05-07 | 东南大学 | 一种路面内部裂缝数量预测方法 |
CN112784495A (zh) * | 2021-01-28 | 2021-05-11 | 郑州轻工业大学 | 一种基于数据驱动的机械结构实时疲劳寿命预测方法 |
CN112834368A (zh) * | 2020-10-27 | 2021-05-25 | 西南交通大学 | 一种受电弓疲劳裂纹扩展寿命的预测方法 |
CN113051768A (zh) * | 2021-04-08 | 2021-06-29 | 广州海关技术中心 | 一种金属疲劳寿命评价方法、系统、设备及介质 |
CN114117652A (zh) * | 2020-08-31 | 2022-03-01 | 重庆大学 | 基于工业ct图像的工件内部裂纹提取和扩展分析的方法 |
CN114117670A (zh) * | 2021-11-24 | 2022-03-01 | 江苏科技大学 | 一种磨削凸轮轴疲劳辉纹间距的估算方法 |
CN114218661A (zh) * | 2022-02-21 | 2022-03-22 | 中国海洋大学 | 一种基于疲劳裂纹扩展的疲劳寿命预测方法 |
CN116773374A (zh) * | 2023-06-15 | 2023-09-19 | 上海发电设备成套设计研究院有限责任公司 | 核电汽轮机的汽缸应力腐蚀与低周疲劳长寿命监控方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104392122A (zh) * | 2014-11-17 | 2015-03-04 | 北京航空航天大学 | 基于裂纹检出概率模型的概率寿命评估方法 |
-
2016
- 2016-05-17 CN CN201610326164.2A patent/CN105956315B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104392122A (zh) * | 2014-11-17 | 2015-03-04 | 北京航空航天大学 | 基于裂纹检出概率模型的概率寿命评估方法 |
Non-Patent Citations (8)
Title |
---|
HAMED HADDAD KHODAPARAST等: ""Perturbation methods for the estimation of parameter variability in stochastic model updating"", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》 * |
PETER A. HESS等: ""Mechanisms of elevated temperature fatigue crack growth in Zr–Ti–Cu–Ni–Be bulk metallic glass"", 《ACTA MATERIALIA》 * |
ZHIPING QIU等: ""Fatigue crack growth modeling and prediction with uncertainties via stochastic perturbation series expansion method"", 《INTERNATIONAL JOURNAL OF MECHANICAL SCIENCES》 * |
ZHIPING QIU等: ""Predicting fatigue crack growth evolution via perturbation series expansion method based on the generalized multinomial theorem"", 《THEORETICAL AND APPLIED FRACTURE MECHANICS》 * |
潘灏: ""裂纹扩展问题渐近场分析的非一致有效性"", 《钢铁研究学报》 * |
王军等: ""疲劳裂纹扩展的不确定理论"", 《北京航空航天大学学报》 * |
苏长青等: ""转轴裂纹扩展的可靠性灵敏度分析"", 《航空动力学报》 * |
金雅娟等: ""任意分布参数的涡轮盘裂纹扩展寿命可靠性分析"", 《工程设计学报》 * |
Cited By (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106568660B (zh) * | 2016-10-14 | 2019-06-25 | 中国人民解放军海军航空工程学院青岛校区 | 一种复合材料胶接修补结构的剩余疲劳寿命的预测方法 |
CN106568660A (zh) * | 2016-10-14 | 2017-04-19 | 中国人民解放军海军航空工程学院青岛校区 | 一种复合材料胶接修补结构的剩余疲劳寿命的预测方法 |
CN106525365A (zh) * | 2016-10-28 | 2017-03-22 | 中航动力股份有限公司 | 一种航空发动机轮盘裂纹予制试验系统及试验方法 |
CN106525365B (zh) * | 2016-10-28 | 2019-01-04 | 中航动力股份有限公司 | 一种航空发动机轮盘裂纹予制试验系统及试验方法 |
CN107389471A (zh) * | 2017-07-13 | 2017-11-24 | 电子科技大学 | 裂纹失效模式判定方法及基于该方法的疲劳寿命预测方法 |
CN107449894A (zh) * | 2017-08-23 | 2017-12-08 | 中国人民解放军海军工程大学 | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 |
CN107449894B (zh) * | 2017-08-23 | 2019-08-20 | 中国人民解放军海军工程大学 | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 |
CN108334666A (zh) * | 2018-01-08 | 2018-07-27 | 中车青岛四方机车车辆股份有限公司 | 一种车轴抗疲劳断裂分析方法 |
CN108828069A (zh) * | 2018-06-19 | 2018-11-16 | 北京领示科技有限公司 | 一种基于超声定量数据的关键件疲劳寿命的定寿方法 |
CN109142529A (zh) * | 2018-08-27 | 2019-01-04 | 佛山科学技术学院 | 一种高强钛合金电子束焊接接头超高周疲劳寿命预测方法 |
CN109142529B (zh) * | 2018-08-27 | 2021-08-03 | 佛山科学技术学院 | 一种高强钛合金电子束焊接接头超高周疲劳寿命预测方法 |
CN109708966A (zh) * | 2018-12-04 | 2019-05-03 | 中国大唐集团科学技术研究院有限公司火力发电技术研究院 | 一种电站高温承压管道焊接接头多缺陷安全性评定方法 |
CN109918776B (zh) * | 2019-03-05 | 2020-09-11 | 北京航空航天大学 | 基于双步最小二乘法的机械产品的疲劳裂纹预测方法 |
CN109918776A (zh) * | 2019-03-05 | 2019-06-21 | 北京航空航天大学 | 基于双步最小二乘法的机械产品的疲劳裂纹预测方法 |
CN110489848A (zh) * | 2019-08-12 | 2019-11-22 | 武汉钢铁有限公司 | 一种不同海水流速腐蚀疲劳裂纹扩展速率预测方法 |
CN110907270B (zh) * | 2019-11-29 | 2022-05-27 | 浙大宁波理工学院 | 利用铁磁性材料的弱磁信号预测疲劳寿命的方法 |
CN110907270A (zh) * | 2019-11-29 | 2020-03-24 | 浙江大学宁波理工学院 | 利用铁磁性材料的弱磁信号预测疲劳寿命的方法 |
CN111044367A (zh) * | 2019-12-31 | 2020-04-21 | 成都理工大学 | 一种基于三轴应力-应变曲线的岩石裂缝亚临界扩展速率实验测试方法 |
CN111044367B (zh) * | 2019-12-31 | 2020-12-18 | 成都理工大学 | 一种基于三轴应力-应变曲线的岩石裂缝亚临界扩展速率实验测试方法 |
US10989639B1 (en) | 2019-12-31 | 2021-04-27 | Chengdu University Of Technology | Experimental test method for subcritical propagation rate of rock fractures based on triaxial stress—strain curve |
CN114117652A (zh) * | 2020-08-31 | 2022-03-01 | 重庆大学 | 基于工业ct图像的工件内部裂纹提取和扩展分析的方法 |
CN111815694A (zh) * | 2020-09-14 | 2020-10-23 | 湖南大学 | 一种疲劳裂纹扩展寿命预测方法、装置、设备及存储介质 |
CN112834368A (zh) * | 2020-10-27 | 2021-05-25 | 西南交通大学 | 一种受电弓疲劳裂纹扩展寿命的预测方法 |
CN112380656A (zh) * | 2020-11-20 | 2021-02-19 | 西安热工研究院有限公司 | 一种燃气轮机燃烧室部件裂纹扩展寿命评估方法 |
CN112464572A (zh) * | 2020-12-10 | 2021-03-09 | 北京航空航天大学 | 一种传动轴加工性能评价及基于关键工艺的控制方法 |
CN112699479A (zh) * | 2020-12-29 | 2021-04-23 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机结构裂纹形成寿命理论值分散系数的确定方法 |
CN112699479B (zh) * | 2020-12-29 | 2024-04-09 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机结构裂纹形成寿命理论值分散系数的确定方法 |
CN112765854B (zh) * | 2021-01-21 | 2024-01-09 | 东南大学 | 一种路面内部裂缝数量预测方法 |
CN112765854A (zh) * | 2021-01-21 | 2021-05-07 | 东南大学 | 一种路面内部裂缝数量预测方法 |
CN112784495B (zh) * | 2021-01-28 | 2021-09-24 | 郑州轻工业大学 | 一种基于数据驱动的机械结构实时疲劳寿命预测方法 |
CN112784495A (zh) * | 2021-01-28 | 2021-05-11 | 郑州轻工业大学 | 一种基于数据驱动的机械结构实时疲劳寿命预测方法 |
CN113051768A (zh) * | 2021-04-08 | 2021-06-29 | 广州海关技术中心 | 一种金属疲劳寿命评价方法、系统、设备及介质 |
CN113051768B (zh) * | 2021-04-08 | 2023-10-17 | 广州海关技术中心 | 一种金属疲劳寿命评价方法、系统、设备及介质 |
CN114117670A (zh) * | 2021-11-24 | 2022-03-01 | 江苏科技大学 | 一种磨削凸轮轴疲劳辉纹间距的估算方法 |
CN114117670B (zh) * | 2021-11-24 | 2024-05-31 | 江苏科技大学 | 一种磨削凸轮轴疲劳辉纹间距的估算方法 |
CN114218661A (zh) * | 2022-02-21 | 2022-03-22 | 中国海洋大学 | 一种基于疲劳裂纹扩展的疲劳寿命预测方法 |
CN114218661B (zh) * | 2022-02-21 | 2022-06-03 | 中国海洋大学 | 一种基于疲劳裂纹扩展的疲劳寿命预测方法 |
CN116773374A (zh) * | 2023-06-15 | 2023-09-19 | 上海发电设备成套设计研究院有限责任公司 | 核电汽轮机的汽缸应力腐蚀与低周疲劳长寿命监控方法 |
CN116773374B (zh) * | 2023-06-15 | 2024-05-31 | 上海发电设备成套设计研究院有限责任公司 | 核电汽轮机的汽缸应力腐蚀与低周疲劳长寿命监控方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105956315B (zh) | 2018-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105956315A (zh) | 一种可进行疲劳裂纹扩展速率估算和寿命预测的方法 | |
Rege et al. | A one-parameter nonlinear fatigue damage accumulation model | |
Madge et al. | A combined wear and crack nucleation–propagation methodology for fretting fatigue prediction | |
Wang et al. | Multi-axial creep-fatigue life prediction considering history-dependent damage evolution: A new numerical procedure and experimental validation | |
Saad et al. | Characterization of viscoplasticity behaviour of P91 and P92 power plant steels | |
Božić et al. | Multiscale fatigue crack growth modelling for welded stiffened panels | |
Wang et al. | Effect of constraint induced by crack depth on creep crack-tip stress field in CT specimens | |
Aid et al. | An equivalent stress process for fatigue life estimation under multiaxial loadings based on a new non linear damage model | |
Baxevanis et al. | On the fracture toughness of pseudoelastic shape memory alloys | |
Xiong et al. | A load history generation approach for full-scale accelerated fatigue tests | |
Sun et al. | A multi-scale damage model for fatigue accumulation due to short cracks nucleation and growth | |
Wen et al. | Creep crack growth by grain boundary cavitation under monotonic and cyclic loading | |
Levkovitch et al. | Simulation of deformation and lifetime behavior of a fcc single crystal superalloy at high temperature under low-cycle fatigue loading | |
Kamaya | Assessment of thermal fatigue damage caused by local fluid temperature fluctuation (part I: characteristics of constraint and stress caused by thermal striation and stratification) | |
Wang et al. | An improved viscoplastic constitutive model and its application to creep behavior of turbine blade | |
Hou et al. | Low cycle fatigue behavior of single crystal superalloy with temperature gradient | |
Chen et al. | Multiaxial thermo-mechanical fatigue life prediction based on notch local stress-strain estimation considering temperature change | |
Singh et al. | Effect of direction of approach of test temperature on fracture toughness of Zr–2.5 Nb pressure tube material | |
Zuo et al. | Stochastic fatigue life and reliability prediction based on residual strength | |
Fan et al. | Modeling of fatigue crack growth of stainless steel 304L | |
Hurley et al. | A new method for predicting fatigue crack propagation rates | |
Voothaluru et al. | Determination of lattice level energy efficiency for fatigue crack initiation | |
CN116306083A (zh) | 基于位错攀移机制的晶体塑性有限元模拟方法 | |
CN104881705A (zh) | 一种基于最优组合模型预测一次风机状态趋势的方法 | |
Hellebrand et al. | A multiscale approach to investigate residual stresses due to targeted cooling of hot bulk formed parts |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180309 |
|
CF01 | Termination of patent right due to non-payment of annual fee |