CN107103140A - 一种基于双线性累积损伤的时变疲劳可靠性分析方法 - Google Patents
一种基于双线性累积损伤的时变疲劳可靠性分析方法 Download PDFInfo
- Publication number
- CN107103140A CN107103140A CN201710290972.2A CN201710290972A CN107103140A CN 107103140 A CN107103140 A CN 107103140A CN 201710290972 A CN201710290972 A CN 201710290972A CN 107103140 A CN107103140 A CN 107103140A
- Authority
- CN
- China
- Prior art keywords
- fatigue
- life
- damage
- accumulated damage
- bilinearity
- 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.)
- Pending
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于双线性累积损伤的时变疲劳可靠性分析方法,包括:A、确定疲劳寿命数据的概率分布类型及其统计特征参量;B、选取合适的参考寿命建立双线性累积损伤物理模型;C、根据双线性累积损伤物理模型建立概率累积损伤模型;D、根据应力‑强度干涉理论建立时变疲劳可靠度模型;E、绘制可靠度曲线并与试验结果进行比较;综合考虑疲劳失效过程裂纹萌生和裂纹扩展的两阶段特性以及累积损伤过程中的随机性,以双线性损伤理论为建模基础,构建概率累积损伤模型,结合应力‑强度干涉理论,建立基于“累积损伤‑临界损伤”时变可靠度模型,体现疲劳载荷作用的时变特性以及累积损伤过程的动态行为,实现了疲劳全寿命周期内的可靠度预测。
Description
技术领域
本发明属于疲劳可靠性分析与疲劳寿命评估技术领域,具体涉及一种疲劳失效的可靠性分析方法。
背景技术
疲劳破坏目前已成为大多数机械装备结构的主要失效模式之一。据统计,机械零构件由疲劳而引发的断裂破坏已高达50%~90%。疲劳现象已遍布航空航天、轨道交通、船舶海洋、桥梁建筑、机械制造等多个领域,由此引发的疲劳事故已层出不穷。疲劳破坏具有很强的隐蔽性和突发性,破坏前无明显的征兆,对机械装备的安全运行构成了严重威胁,一旦发生破坏易于造成重大事故和生命财产损失。机械装备零构件或结构的疲劳寿命及其可靠性是制约装备整机寿命和系统可靠性水平的关键因素。机械结构的疲劳失效是结构内部损伤不断累积的过程,当累积损伤达到损伤阈值时即发生疲劳破坏。在工程实际中,疲劳损伤的特征表现为累积损伤的不可逆性和随机性,其随机性包含构成材料疲劳属性的内在分散性和与服役载荷、环境等因素相关的外在分散性。由于疲劳失效过程具有随机性和不确定性,疲劳可靠性分析是结构抗疲劳设计的重要内容之一。基于疲劳累积损伤的可靠性建模既要符合疲劳失效的物理机制,还要能正确描述累积损伤过程中的不确定性。
疲劳累积损伤理论为结构寿命预测和可靠性分析提供了有效途径,根据累积损伤的形式,该理论可分为线性、非线性和双线性损伤法则。其中,线性损伤法则或Miner法则因其形式简单在结构疲劳分析中得到了广泛应用,但该法则存在以下两个缺陷:①它对疲劳失效机制进行了简化,未考虑载荷加载的顺序效应;②它是一个确定性的累积损伤理论,不能反映累积损伤的分散性和不确定性。针对第一个缺陷,许多学者提出了各种形式的非线性损伤理论对Miner法则进行修正,这些非线性模型具有特定的物理基础,能够较好地解释载荷顺序效应。但是,该理论同样属于确定性理论,且模型参数较多,可靠性建模过程复杂,难以在工程实际中得到应用。针对第二个缺陷,现有常用的方法是通过引入概率统计理论进行不确定性建模与分析,量化或表征累积损伤过程中的随机性。但这些建模方法大多仍以Miner法则为依据,无法揭示失效过程的物理本质。
从微观物理角度分析,疲劳全寿命周期包含疲劳裂纹萌生和疲劳裂纹扩展两个阶段。双线性累积损伤理论正是以这两个阶段为基础,在每个阶段下均采用线性累积损伤法则进行建模。该理论形式简单,能正确描述失效过程的两阶段特性,物理意义明确。为此,本发明以双线性损伤理论为基础建立概率累积损伤模型。同时,考虑到疲劳载荷作用的时变特性,结合应力-强度干涉模型,构建基于“累积损伤-临界损伤”动态干涉模型,进行时变疲劳可靠性分析,为机械装备的安全评定、维修决策和健康管理提供理论支撑。
发明内容
本发明为解决上述技术问题,提出了一种基于双线性累积损伤的时变疲劳可靠性分析方法,综合考虑疲劳失效过程裂纹萌生和裂纹扩展的两阶段特性以及累积损伤过程中的随机性,构建基于“累积损伤-临界损伤”时变可靠度模型,实现了疲劳全寿命周期内的可靠度预测。
本发明采用的技术方案是:一种基于双线性累积损伤的时变疲劳可靠性分析方法,包括:
A、确定疲劳寿命数据的概率分布类型及其统计特征参量;
B、选取合适的参考寿命建立双线性累积损伤物理模型;
C、根据双线性累积损伤物理模型建立概率累积损伤模型;
D、根据应力-强度干涉理论建立时变疲劳可靠度模型;
E、绘制可靠度曲线并与试验结果进行比较。
进一步地,所述步骤A具体包括:
A1、根据等幅加载疲劳试验,分别收集整理不同应力级别下的若干组疲劳寿命试验数据;
A2、将步骤A1收集到的疲劳寿命试验数据按从小到大进行排序,计算各应力级别对应的疲劳寿命均值和疲劳寿命标准差;
A3、将步骤A1收集到的疲劳寿命试验数据进行对数标准化处理,形成对数寿命数据序列,并计算得到各应力级别对应的疲劳寿命数据的统计特征参量,所述疲劳寿命数据的统计特征参量包括:疲劳寿命的对数均值和疲劳寿命的对数标准差;
A4、根据步骤A3的对数寿命数据序列绘制概率纸,确定疲劳寿命的概率分布类型。
进一步地,所述步骤B具体包括:
B1、将最大应力级别对应的疲劳寿命均值作为参考寿命水平;选取疲劳极限附近的寿命水平107作为基准寿命水平,并绘制双线性损伤中值曲线;
B2、根据双线性累积损伤的拐点坐标公式和插值公式,计算步骤B1所述双线性中值累积损伤曲线中的未知参数,所述未知参数包括:裂纹萌生寿命、裂纹扩展寿命、裂纹萌生寿命分数以及裂纹萌生与裂纹扩展两个阶段下的损伤演化速率;
B3、根据裂纹萌生寿命、裂纹扩展寿命、裂纹萌生寿命分数以及两个阶段下的损伤演化速率,建立双线性累积损伤物理模型。
进一步地,所述步骤C包括:
C1、在等幅疲劳加载下,根据疲劳寿命的概率分布类型、疲劳寿命的对数均值以及疲劳寿命的对数标准差,得到裂纹萌生和裂纹扩展两个阶段下的累积损伤概率统计特性;
C2、根据步骤C1得到的裂纹萌生和裂纹扩展阶段下累积损伤概率统计特性,确定临界损伤的概率分布类型以及其统计特征参量;
C3、根据步骤C2得到的临界损伤的概率统计特征参量,确定累积损伤的概率分布类型以及其统计特征参量。
进一步地,所述步骤D具体为:
D1、在等幅疲劳加载试验条件下,根据临界损伤概率分布类型、临界损伤的统计特征参量、累积损伤的概率分布类型以及累积损伤的统计特征参量,建立结构的安全余量方程;
D2、根据结构安全余量方程,结合裂纹萌生与裂纹扩展两个阶段,根据应力-强度干涉理论构建等幅加载时变疲劳可靠度模型;
D3、收集同类构件在变幅疲劳载荷谱试验下的失效寿命数据,根据应力-强度干涉理论构建变幅加载时变疲劳可靠度模型。
进一步地,步骤E所述绘制可靠度曲线包括:根据等幅加载试验下的疲劳可靠度模型,绘制等幅加载下的可靠度曲线;
以及根据变幅加载时变疲劳可靠度模型,绘制变幅加载下的可靠度变化曲线。
本发明的有益效果:本发明通过综合考虑疲劳失效过程裂纹萌生和裂纹扩展的两阶段特性以及累积损伤过程中的随机性,以双线性累积损伤理论为建模基础,通过分析疲劳寿命数据的概率分布,采用工程中当量化分布的方法确定临界损伤和累积损伤的概率分布类型及其统计特征参量,构建双线性累积损伤物理模型;并在此基础上,应用传统应力-强度干涉理论,构建了基于“累积损伤-临界损伤”时变可靠度模型,该模型能够体现疲劳载荷作用的时变特性以及累积损伤的动态行为,实现了疲劳全寿命周期内的可靠度预测。
附图说明
图1为本发明提供的一种基于双线性累积损伤的时变疲劳可靠性分析方法流程图。
图2为本发明实施例的带切口的45钢摩擦焊接头的标准试件。
图3为本发明实施例的5种不同应力水平下对数疲劳寿命概率纸。
图4为本发明实施例的5种不同应力水平下双线性疲劳累积损伤中值曲线示意图。
图5为本发明实施例的带切口的45钢摩擦焊接头程序块载荷谱。
图6为本发明实施例的恒幅加载下带切口的45钢摩擦焊接头可靠度预测曲线与试验结果比较。
图7为本发明实施例的程序块加载下带切口的45钢摩擦焊接头可靠度预测曲线与试验结果比较。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
以带切口的45钢摩擦焊接头疲劳试验为例对本发明作进一步说明,如图1所示为本发明提供的一种基于双线性累积损伤的时变疲劳可靠性分析方法的流程图,具体包括以下步骤:
A、确定疲劳寿命数据的概率分布类型及其统计特征参量;
B、选取合适的参考寿命建立双线性累积损伤物理模型;
C、根据双线性累积损伤物理模型建立概率累积损伤模型;
D、根据应力-强度干涉理论建立时变疲劳可靠度模型;
E、绘制可靠度曲线并与试验结果进行比较。
所述步骤A确定疲劳寿命数据的概率分布类型及其统计特征参量的具体方法为:
A1、根据等幅加载疲劳试验,分别收集整理不同应力级别下的若干组疲劳寿命试验数据;
以带切口的45钢摩擦焊接头疲劳试验为例,试件母材选用经正火处理的45钢,其化学组分为:0.65%Mn、0.51%C、0.22%Si和余量Fe;力学性能为:屈服强度σy=441MPa,抗拉强度σu=703MPa,延伸率δ10=15.3%,断面收缩率ψ=50.7%;摩擦焊接工艺参数为:摩擦压力为σ摩擦=130MPa,摩擦时间t摩擦=1.4s,顶锻压力σ顶锻=270MPa,保压时间t顶锻=4s。45钢摩擦焊接头经车、磨处理加工成带U型切口的标准试件,如图2所示。试件的应力集中系数为Kt=2。
在等幅疲劳加载下,焊接接头试件在PQ-6型疲劳试验机上进行旋转弯曲加载,加载模式为应力控制,应力比R=-1,加载频率为f=50Hz。试验加载的应力级别分为5种,相应的等效应力幅分别为750MPa、650MPa、630MPa、590MPa、520MPa。每个应力级别分别选取10组样本进行试验,最终收集得到5个不同应力级别下的50组疲劳寿命试验数据,如表1所示。
表1
A2、将步骤A1收集到的疲劳寿命试验数据按从小到大进行排序,计算各应力级别对应的疲劳寿命均值和疲劳寿命标准差;
表1所示的试验寿命数据已按从小到大顺序进行排序,分别采用式(1)和式(2)计算各级应力水平下的疲劳寿命均值和疲劳寿命标准差,计算结果如表2所示。
其中,N为样本量,在本实施例中,各级应力水平下对应的样本量为10,即本申请实施例中N为10,1≤i≤N;Nfi为第i个试验疲劳寿命,为疲劳寿命均值,为疲劳寿命标准差。
表2
A3、将步骤A1收集到的疲劳寿命试验数据进行对数标准化处理,形成对数寿命数据序列,计算得到疲劳寿命数据的统计特征参量,即疲劳寿命的对数均值和疲劳寿命的对数标准差;
对表1中不同应力级别下的试验寿命数据取对数形成对数疲劳寿命lnNf的数据序列,结合式(1)和式(2)计算寿命的对数均值和对数标准差,计算结果如表3所示。
表3
A4、根据步骤A3的对数寿命数据序列绘制概率纸,确定疲劳寿命的概率分布类型。
根据表3对数寿命数据序列,在Matlab运行环境下调用命令“Normal probabilityplot”绘制概率纸判断数据序列是否服从正态分布,最终得到不同应力级别下寿命与失效概率之间的关系,如图3所示。
由图3可知,对数寿命数据点基本都位于一条直线上,表明不同应力级别下的对数寿命比较接近于正态分布。此外,调用Matlab函数“h=lillietest(x)”对数据序列进行正态分布假设检验,计算结果均为“h=0”,接受正态分布。由此可得,疲劳寿命可认为服从对数正态分布。
所述步骤B选取合适的参考寿命建立双线性累积损伤物理模型的具体方法为:
B1、将最大应力级别对应的疲劳寿命均值作为参考寿命水平;在不考虑疲劳极限以下小载荷对损伤贡献的情况下,同时避免外推计算或因改变寿命参考点而重新计算拐点,选取疲劳极限附近的寿命水平107作为基准寿命水平,并绘制双线性损伤中值曲线;
由表1可知,试验加载的最大应力水平为σ1=750MPa,其均值寿命为则将设定为参考寿命水平,即为Nref=32835;基准寿命水平设定为Nfb=107。通常,将最大应力水平下的损伤曲线简化为线性Miner法则,由此,可绘制其他4种应力水平下的中值累积损伤曲线,如图4所示。在图4中,疲劳累积损伤分为两个阶段,即阶段I和阶段II,它们分别对应着裂纹萌生和裂纹扩展两个阶段,两个阶段的过渡处存在拐点,Di为裂纹萌生损伤,Dc为临界损伤,Ni为裂纹萌生寿命。
B2、根据双线性累积损伤的拐点坐标公式和插值公式,计算步骤B1中所绘制双线性损伤中值曲线的未知参数,所述未知参数包括:裂纹萌生寿命、裂纹扩展寿命、裂纹萌生寿命分数以及裂纹萌生与裂纹扩展两个阶段下的损伤演化速率;
根据图4,按照双线性累积损伤理论的计算约定,以Nref=32835为参考寿命、Nfb=107为基准寿命,可采用式(3)计算裂纹萌生损伤:
裂纹萌生寿命采用插值公式计算,描述如下:
其中,α′为裂纹萌生阶段的寿命分数,α′∈[0,1]。和Z的表达式如下:
裂纹扩展寿命的计算公式为:
其中,Np为裂纹扩展寿命。
两个阶段的疲劳损伤演化速率的计算公式分别为:
其中,kI和kII分别表示阶段I和阶段II下的疲劳损伤演化速率。
结合表2中不同应力水平下的疲劳寿命均值数据,对上述疲劳损伤参数进行计算,计算结果如表4所示。
表4
B3、根据裂纹萌生寿命、裂纹扩展寿命、裂纹萌生寿命分数以及两个阶段下的损伤演化速率,建立双线性累积损伤物理模型。
根据双线性累积损伤理论,裂纹萌生和裂纹扩展两个阶段下的累积损伤与载荷循环次数均具有线性关系;以裂纹萌生寿命Ni为分界点,累积损伤可用以下分段函数描述:
其中,n为加载循环次数。
将表4中的疲劳损伤参数代入式(10),由此可构建不同应力水平下的双线性累积损伤物理模型。
所述步骤C根据双线性累积损伤物理模型建立概率累积损伤模型的具体方法为:
C1、在等幅疲劳加载下,根据疲劳寿命的概率分布类型、疲劳寿命的对数均值以及疲劳寿命的对数标准差,得到裂纹萌生和裂纹扩展两个阶段下的累积损伤概率统计特性;
由于疲劳寿命Nf具有统计特性,对于给定应力水平下的使用寿命或加载循环次数n,相应地也是具有某种统计特征的随机变量。在概率建模过程中,如何确定使用寿命的变异性或标准差是求解累积损伤分布的关键。根据之前的分析,在等幅疲劳加载下,疲劳寿命是服从对数正态分布的随机变量,即或对数使用寿命的变异性可用下式描述:
其中,ln n是随机变量,其概率分布与ln Nf一致,都服从正态分布。
根据式(10)得到的双线性累积损伤物理模型,不同阶段下的累积损伤可描述为:
其中,nI和nII分别是阶段I和阶段II下的加载循环次数,D(n)I和D(n)II分别是阶段I和阶段II下的累积损伤。
对式(12)两边取对数可得:
由此可知,ln D(n)与ln n之间具有线性关系。由于ln n服从正态分布,则D(n)I和D(n)II均服从对数正态分布,即和ln D(n)I和ln D(n)II的统计特征参数分别用式(14)和(15)表示。
C2、根据步骤C1得到的裂纹萌生和裂纹扩展阶段下累积损伤概率统计特性,确定临界损伤的概率分布类型以及其统计特征参量;
由式(10)和式(13),将累积损伤公式改写为:
在疲劳失效状态下,临界损伤Dc由两个阶段损伤构成,即裂纹萌生损伤D(Ni)I和裂纹扩展损伤D(Np)II,且D(Ni)I和D(Np)II均是服从对数正态分布的随机变量,因此Dc的概率分布无法用显式描述。在工程中,可对Dc进行当量化处理,使其概率分布仍服从对数正态分布。
当量化处理的条件是要保持变量当量化前后的矩相等,如下式:
其中,E[□]表示变量的一阶矩,即Dc的均值,Var[□]表示变量的二阶矩,即Dc的方差。
根据对数正态分布的性质,D(Ni)I的均值和方差可通过和转换得到,即
同理,可得D(Np)II的均值和方差,如下:
将式(18)和(19)代入式(17),可得:
再次利用对数正态分布的性质,ln Dc的统计特征参数可由和转换得到:
C3、根据步骤C2得到的临界损伤的概率统计特征参量,确定累积损伤的概率分布类型以及统计特征参量。
根据临界损伤概率统计特征参量的分析过程,可推导得到任意使用寿命下的瞬时累积损伤统计特征参量。根据双线性累积损伤的两阶段特性,累积损伤的统计参量可分阶段描述,描述如下:
1)阶段I:裂纹萌生阶段,n≤Ni,ln D(n)的均值和标准差为:
2)阶段II:裂纹扩展阶段,Ni<n≤Nf,有:
对式(23)进行对数转换,可得:
式(22)和式(24)的统计特征参数可根据疲劳损伤演化速率、疲劳寿命的对数均值及其对数标准差求出,并用于疲劳可靠性分析。
所述步骤D根据应力-强度干涉理论建立时变疲劳可靠度模型的具体方法为:
D1、在等幅疲劳加载试验条件下,根据临界损伤概率分布模型、临界损伤的统计特征参量、累积损伤的概率分布类型以及累积损伤的统计特征参量,建立结构的安全余量方程;
在结构疲劳设计中,应力-强度干涉模型是结构疲劳可靠性分析的基础,若结构内部累积的损伤量小于临界损伤,则结构处于安全状态,否则处于失效状态。根据前面的分析,累积损伤D(n)和临界损伤Dc均是服从对数正态分布的随机变量。若不考虑D(n)与Dc之间的相关性,建立结构的安全余量方程为:
G(n)=ln Dc-ln D(n) (25)
由于lnDc和lnD(n)服从正态分布,则G(n)也服从正态分布,得结构可靠度函数:
D2、根据结构安全余量方程,结合裂纹萌生与裂纹扩展两个阶段,根据应力-强度干涉理论构建等幅加载时变疲劳可靠度模型;
在等幅加载下,结构可靠度函数可分以下两个阶段描述。
1)阶段I:裂纹萌生阶段,n≤Ni:
2)阶段II:裂纹扩展阶段,Ni<n≤Nf:
D3、收集同类构件在变幅疲劳载荷谱试验下的失效寿命数据,根据应力-强度干涉理论构建变幅加载时变疲劳可靠度模型。
在变幅加载下,仍以带切口45钢摩擦焊接头为研究对象进行程序块载荷谱疲劳试验,试验载荷谱如图5所示。试验加载参数及失效寿命数据如表5所示。
表5
在变幅加载下,不同载荷阶段下的疲劳损伤均假设服从对数正态分布,因此累积损伤变量是由一系列对数正态分布的随机变量组成,其概率密度函数无法用显式表征。在此,采用当量化分布方法,将累积损伤和∑D(n)简化为服从对数正态分布的随机变量。在疲劳失效状态下,根据当量化处理的条件,可得临界损伤Dc的均值和标准差为:
利用对数正态分布的性质,有:
同理,可得ln(∑D(n))的统计特征参量为:
将式(30)和(31)代入式(26),可得变幅加载下结构可靠度函数为:
等幅加载下式(27)和式(28)以及变幅加载下式(32)所示的可靠度模型均是使用寿命的函数,体现了累积损伤过程的动态行为。在全寿命周期内,使用寿命与可靠度之间具有一一对应关系,从而实现了时变疲劳可靠度估算。
所述步骤E具体为:根据对数正态分布下概率累积损伤建模分析,在等幅加载下利用式(27)和式(28)的等幅加载下结构可靠度函数,绘制5种不同应力水平下的可靠度曲线,并与试验寿命数据计算的经验分布可靠度进行比较,结果如图6所示。
同理,在变幅加载下利用式(32)的变幅加载下结构可靠度函数,绘制带切口的45钢摩擦焊接头在不同程序块加载阶段下的可靠度变化曲线,如图7所示。
在图6中,可靠度函数随载荷循环次数增加单调递减,具有时变特性。在疲劳全寿命周期内,可靠度曲线分为两个阶段,其中阶段I和阶段II分别代表裂纹萌生和裂纹扩展。在裂纹萌生阶段,可靠度水平较高,处于稳定状态,其值与初始寿命状态下的可靠度(R=1)几乎相同;同时,随着应力水平的减小,裂纹萌生寿命增大,此阶段下的可靠度周期也随之延长。在裂纹扩展阶段,随载荷循环次数增加,可靠度曲线急剧下降。上述现象可解释为:在裂纹萌生阶段,疲劳累积损伤速率相对较低,构件内部累积的损伤量较少,发生失效的概率较小,使该阶段下的可靠度仍保持在较高水平状态且变化平稳;随着疲劳载荷继续作用,当进入裂纹扩展阶段时,裂纹增长速率和累积损伤速率迅速增加,导致构件内部的累积损伤快速趋于临界损伤,其发生失效的概率大大增加,因此该阶段对应的可靠度急剧下降。
在图7中,可靠度随载荷循环逐渐衰减,体现了载荷作用的时变特性,预测曲线与试验结果大致吻合。在初始的几个程序块加载下(疲劳失效前期),可靠度水平很高,几乎无变化;在最后两个程序块加载下(疲劳失效后期),可靠度曲线迅速下降。在疲劳失效前期,构件处于裂纹萌生阶段,累积损伤量较少;而在失效后期,构件进入裂纹扩展阶段,此时累积损伤速率大大增加,累积损伤量快速趋于临界值,导致可靠度水平迅速降低。上述程序块谱下的可靠度变化也同时体现了疲劳失效的两阶段特性。
本发明通过综合考虑疲劳失效过程裂纹萌生和裂纹扩展的两阶段特性以及累积损伤过程中的随机性,以双线性累积损伤理论为建模基础,通过分析疲劳寿命数据的概率分布,采用当量化分布的方法确定临界损伤和累积损伤的概率分布类型及其统计特征参量,建立了概率累积损伤模型。在此基础上,应用传统应力-强度干涉理论,构建了基于“累积损伤-临界损伤”时变可靠度模型,该模型能够体现疲劳载荷作用的时变特性以及累积损伤的动态行为,实现了疲劳全寿命周期内的可靠度预测。理论模型的预测结果与试验结果吻合较好,且能正确反映疲劳失效过程裂纹萌生和裂纹扩展的两阶段特性,为现阶段结构疲劳可靠性、寿命分析与评估提供了有效方法和途径,同时,也为后续开展结构的安全评定、维修决策和健康管理奠定基础。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (5)
1.一种基于双线性累积损伤的时变疲劳可靠性分析方法,其特征在于,具体包括:
A、确定疲劳寿命数据的概率分布类型及其统计特征参量;
B、选取合适的参考寿命建立双线性累积损伤物理模型;
C、根据双线性累积损伤物理模型建立概率累积损伤模型;
D、根据应力-强度干涉理论建立时变疲劳可靠度模型;
E、绘制可靠度曲线并与试验结果进行比较。
2.根据权利要求1所述的基于双线性累积损伤的时变疲劳可靠性分析方法,其特征在于,所述步骤A的具体方法为:
A1、根据等幅加载疲劳试验,分别收集整理不同应力级别下的若干组疲劳寿命试验数据;
A2、将步骤A1收集到的疲劳寿命试验数据按从小到大进行排序,计算各应力级别对应的疲劳寿命均值和疲劳寿命标准差;
A3、将步骤A1收集到的疲劳寿命试验数据进行对数标准化处理,形成对数寿命数据序列,并计算得到各应力级别对应的疲劳寿命数据的统计特征参量,所述疲劳寿命数据的统计特征参量包括:疲劳寿命的对数均值和疲劳寿命的对数标准差;
A4、根据步骤A3的对数寿命数据序列绘制概率纸,确定疲劳寿命的概率分布类型。
3.根据权利要求1所述的基于双线性累积损伤的时变疲劳可靠性分析方法,其特征在于,所述步骤B的具体方法为:
B1、将最大应力级别对应的疲劳寿命均值作为参考寿命水平;选取疲劳极限附近的寿命水平107作为基准寿命水平,绘制双线性中值累积损伤曲线;
B2、根据双线性累积损伤的拐点坐标公式和插值公式,计算步骤B1所述双线性中值累积损伤曲线中的未知参数,所述未知参数包括:裂纹萌生寿命、裂纹扩展寿命、裂纹萌生寿命分数以及裂纹萌生与裂纹扩展两个阶段下的损伤演化速率;
B3、根据裂纹萌生寿命、裂纹扩展寿命、裂纹萌生寿命分数以及两个阶段下的损伤演化速率,建立双线性累积损伤物理模型。
4.根据权利要求1所述的基于双线性累积损伤的时变疲劳可靠性分析方法,其特征在于,所述步骤C的具体方法为:
C1、在等幅疲劳加载下,根据疲劳寿命的概率分布类型、疲劳寿命的对数均值以及疲劳寿命的对数标准差,得到裂纹萌生和裂纹扩展两个阶段下的累积损伤概率统计特性;
C2、根据步骤C1得到的裂纹萌生和裂纹扩展阶段下累积损伤概率统计特性,确定临界损伤的概率分布类型以及其统计特征参量;
C3、根据步骤C2得到的临界损伤的概率统计特征参量,确定累积损伤的概率分布类型以及其统计特征参量。
5.根据权利要求1所述的基于双线性累积损伤的时变疲劳可靠性分析方法,其特征在于,所述步骤D的具体方法为:
D1、在等幅疲劳加载试验条件下,根据临界损伤概率分布类型、临界损伤的统计特征参量、累积损伤的概率分布类型以及累积损伤的统计特征参量,建立结构的安全余量方程;
D2、根据结构安全余量方程,结合裂纹萌生与裂纹扩展两个阶段,根据应力-强度干涉理论构建等幅加载时变疲劳可靠度模型;
D3、收集同类构件在变幅疲劳载荷谱试验下的失效寿命数据,根据应力-强度干涉理论构建变幅加载时变疲劳可靠度模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710290972.2A CN107103140A (zh) | 2017-04-28 | 2017-04-28 | 一种基于双线性累积损伤的时变疲劳可靠性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710290972.2A CN107103140A (zh) | 2017-04-28 | 2017-04-28 | 一种基于双线性累积损伤的时变疲劳可靠性分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107103140A true CN107103140A (zh) | 2017-08-29 |
Family
ID=59657915
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710290972.2A Pending CN107103140A (zh) | 2017-04-28 | 2017-04-28 | 一种基于双线性累积损伤的时变疲劳可靠性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107103140A (zh) |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107729647A (zh) * | 2017-10-13 | 2018-02-23 | 电子科技大学 | 一种将时变可靠性转换成静态可靠性的分析方法 |
CN108182327A (zh) * | 2017-12-30 | 2018-06-19 | 北京工业大学 | 一种基于线性损伤累积的多轴热机械疲劳寿命预测方法 |
CN108664716A (zh) * | 2018-04-27 | 2018-10-16 | 佛山科学技术学院 | 一种基于超高周疲劳性能的高强钛合金热处理工艺评价模型 |
CN108875171A (zh) * | 2018-06-05 | 2018-11-23 | 温州大学 | 一种基于概率干涉和弹性梁理论的有轨电车输电线振动可靠性计算和预测方法 |
CN109241618A (zh) * | 2018-09-04 | 2019-01-18 | 北京航空航天大学 | 一种复合材料结构的高周弯曲疲劳寿命预测方法 |
CN109684697A (zh) * | 2018-12-14 | 2019-04-26 | 中国航空工业集团公司西安飞机设计研究所 | 一种当量损伤模型的确定方法 |
CN110083970A (zh) * | 2019-05-08 | 2019-08-02 | 中国兵器工业第五九研究所 | 弹用高分子材料累积损伤的可视化仿真方法 |
CN110175394A (zh) * | 2019-05-23 | 2019-08-27 | 电子科技大学 | 一种涡轮叶片疲劳-蠕变损伤耦合概率寿命预测计算方法 |
CN110208086A (zh) * | 2019-06-05 | 2019-09-06 | 东北石油大学 | 脉动压裂岩石累积损伤求解方法 |
CN110750899A (zh) * | 2019-10-18 | 2020-02-04 | 大连理工大学 | 基于概率服役寿命代理关系的核主泵系统可靠性分析方法 |
CN110929316A (zh) * | 2019-09-16 | 2020-03-27 | 江苏卓然智能重工有限公司 | 一种大型塔吊缆索疲劳评估中的可靠性分析模型 |
CN111159827A (zh) * | 2019-12-27 | 2020-05-15 | 电子科技大学 | 一种高铁轴承概率疲劳可靠性评估方法 |
CN111581835A (zh) * | 2020-05-14 | 2020-08-25 | 内蒙古工业大学 | 一种机械结构体的安全信息获取方法 |
CN111693387A (zh) * | 2020-06-15 | 2020-09-22 | 一汽解放汽车有限公司 | 确定快速疲劳试验最少次数的方法 |
CN111721647A (zh) * | 2020-06-24 | 2020-09-29 | 四川大学 | 一种低周疲劳测试数据处理与内应力评估方法 |
CN111881564A (zh) * | 2020-07-17 | 2020-11-03 | 北京理工大学 | 一种关于机械结构变幅疲劳寿命预测方法 |
CN112036661A (zh) * | 2020-09-10 | 2020-12-04 | 南京工程学院 | 基于刀具力学性能随机分布性的陶瓷刀具可靠度预测方法 |
CN112149242A (zh) * | 2020-08-26 | 2020-12-29 | 北京航空航天大学 | 一种考虑应力松弛和辐照影响的堆内构件压紧弹簧疲劳可靠性评估方法 |
CN112733332A (zh) * | 2020-12-23 | 2021-04-30 | 北京航空航天大学 | 一种考虑拧紧力矩的螺栓疲劳可靠性评估方法 |
CN112836361A (zh) * | 2021-01-21 | 2021-05-25 | 新疆大学 | 多级循环载荷作用下材料疲劳寿命预测方法 |
CN113378310A (zh) * | 2021-04-29 | 2021-09-10 | 北京航空航天大学 | 一种基于不确定微分方程的疲劳裂纹扩展建模方法 |
CN113536490A (zh) * | 2021-08-27 | 2021-10-22 | 西安石油大学 | 一种抽油杆疲劳寿命预测方法及系统 |
CN113591340A (zh) * | 2021-06-24 | 2021-11-02 | 内蒙古工业大学 | 随机载荷下齿轮弯曲疲劳时变可靠性分析方法 |
CN113591268A (zh) * | 2021-06-24 | 2021-11-02 | 内蒙古工业大学 | 变幅载荷下齿轮接触疲劳寿命可靠性评估方法及装置 |
CN113742948A (zh) * | 2021-08-23 | 2021-12-03 | 西安石油大学 | 一种超高强度抽油杆p-s-n曲线拟合新模型及方法 |
CN115408755A (zh) * | 2022-09-05 | 2022-11-29 | 北京交通大学 | 一种考虑时变效应的组合梁桥动力疲劳可靠性评估方法 |
CN116049922A (zh) * | 2023-01-12 | 2023-05-02 | 中国工程物理研究院研究生院 | 应力强度因子门槛值概率建模的金属结构疲劳可靠性预测方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106202913A (zh) * | 2016-07-07 | 2016-12-07 | 华东理工大学 | 时间相关的蠕变疲劳损伤评定方法 |
-
2017
- 2017-04-28 CN CN201710290972.2A patent/CN107103140A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106202913A (zh) * | 2016-07-07 | 2016-12-07 | 华东理工大学 | 时间相关的蠕变疲劳损伤评定方法 |
Non-Patent Citations (1)
Title |
---|
ZHAOCHUN PENG, ETC: "Probabilistic Modeling of Double Linear Damage Accumulation for Fatigue Reliability Analysis", 《PROCEEDINGS OF 2016 INTERNATIONAL CONFERENCE ON QUALITY, RELIABILITY, RISK, MAINTENANCE, AND SAFETY ENGINEERING》 * |
Cited By (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107729647B (zh) * | 2017-10-13 | 2020-09-25 | 电子科技大学 | 一种将时变可靠性转换成静态可靠性的分析方法 |
CN107729647A (zh) * | 2017-10-13 | 2018-02-23 | 电子科技大学 | 一种将时变可靠性转换成静态可靠性的分析方法 |
CN108182327A (zh) * | 2017-12-30 | 2018-06-19 | 北京工业大学 | 一种基于线性损伤累积的多轴热机械疲劳寿命预测方法 |
CN108182327B (zh) * | 2017-12-30 | 2021-06-11 | 北京工业大学 | 一种基于线性损伤累积的多轴热机械疲劳寿命预测方法 |
CN108664716B (zh) * | 2018-04-27 | 2021-11-09 | 佛山科学技术学院 | 一种基于超高周疲劳性能的高强钛合金热处理工艺评价模型 |
CN108664716A (zh) * | 2018-04-27 | 2018-10-16 | 佛山科学技术学院 | 一种基于超高周疲劳性能的高强钛合金热处理工艺评价模型 |
CN108875171A (zh) * | 2018-06-05 | 2018-11-23 | 温州大学 | 一种基于概率干涉和弹性梁理论的有轨电车输电线振动可靠性计算和预测方法 |
CN109241618A (zh) * | 2018-09-04 | 2019-01-18 | 北京航空航天大学 | 一种复合材料结构的高周弯曲疲劳寿命预测方法 |
CN109684697A (zh) * | 2018-12-14 | 2019-04-26 | 中国航空工业集团公司西安飞机设计研究所 | 一种当量损伤模型的确定方法 |
CN109684697B (zh) * | 2018-12-14 | 2023-04-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种当量损伤模型的确定方法 |
CN110083970A (zh) * | 2019-05-08 | 2019-08-02 | 中国兵器工业第五九研究所 | 弹用高分子材料累积损伤的可视化仿真方法 |
CN110083970B (zh) * | 2019-05-08 | 2023-05-02 | 中国兵器工业第五九研究所 | 弹用高分子材料累积损伤的可视化仿真方法 |
CN110175394A (zh) * | 2019-05-23 | 2019-08-27 | 电子科技大学 | 一种涡轮叶片疲劳-蠕变损伤耦合概率寿命预测计算方法 |
CN110208086A (zh) * | 2019-06-05 | 2019-09-06 | 东北石油大学 | 脉动压裂岩石累积损伤求解方法 |
CN110929316A (zh) * | 2019-09-16 | 2020-03-27 | 江苏卓然智能重工有限公司 | 一种大型塔吊缆索疲劳评估中的可靠性分析模型 |
CN110750899A (zh) * | 2019-10-18 | 2020-02-04 | 大连理工大学 | 基于概率服役寿命代理关系的核主泵系统可靠性分析方法 |
CN110750899B (zh) * | 2019-10-18 | 2023-08-01 | 大连理工大学 | 基于概率服役寿命代理关系的核主泵系统可靠性分析方法 |
CN111159827A (zh) * | 2019-12-27 | 2020-05-15 | 电子科技大学 | 一种高铁轴承概率疲劳可靠性评估方法 |
CN111581835A (zh) * | 2020-05-14 | 2020-08-25 | 内蒙古工业大学 | 一种机械结构体的安全信息获取方法 |
CN111693387A (zh) * | 2020-06-15 | 2020-09-22 | 一汽解放汽车有限公司 | 确定快速疲劳试验最少次数的方法 |
CN111721647A (zh) * | 2020-06-24 | 2020-09-29 | 四川大学 | 一种低周疲劳测试数据处理与内应力评估方法 |
CN111721647B (zh) * | 2020-06-24 | 2021-12-28 | 四川大学 | 一种低周疲劳测试数据处理与内应力评估方法 |
CN111881564A (zh) * | 2020-07-17 | 2020-11-03 | 北京理工大学 | 一种关于机械结构变幅疲劳寿命预测方法 |
CN111881564B (zh) * | 2020-07-17 | 2022-09-09 | 北京理工大学 | 一种关于机械结构变幅疲劳寿命预测方法 |
CN112149242A (zh) * | 2020-08-26 | 2020-12-29 | 北京航空航天大学 | 一种考虑应力松弛和辐照影响的堆内构件压紧弹簧疲劳可靠性评估方法 |
CN112036661B (zh) * | 2020-09-10 | 2023-11-21 | 南京工程学院 | 基于刀具力学性能随机分布性的陶瓷刀具可靠度预测方法 |
CN112036661A (zh) * | 2020-09-10 | 2020-12-04 | 南京工程学院 | 基于刀具力学性能随机分布性的陶瓷刀具可靠度预测方法 |
CN112733332B (zh) * | 2020-12-23 | 2022-06-28 | 北京航空航天大学 | 一种考虑拧紧力矩的螺栓疲劳可靠性评估方法 |
CN112733332A (zh) * | 2020-12-23 | 2021-04-30 | 北京航空航天大学 | 一种考虑拧紧力矩的螺栓疲劳可靠性评估方法 |
CN112836361B (zh) * | 2021-01-21 | 2022-10-28 | 新疆大学 | 多级循环载荷作用下材料疲劳寿命预测方法 |
CN112836361A (zh) * | 2021-01-21 | 2021-05-25 | 新疆大学 | 多级循环载荷作用下材料疲劳寿命预测方法 |
CN113378310A (zh) * | 2021-04-29 | 2021-09-10 | 北京航空航天大学 | 一种基于不确定微分方程的疲劳裂纹扩展建模方法 |
CN113591340B (zh) * | 2021-06-24 | 2023-09-01 | 内蒙古工业大学 | 随机载荷下齿轮弯曲疲劳时变可靠性分析方法 |
CN113591268A (zh) * | 2021-06-24 | 2021-11-02 | 内蒙古工业大学 | 变幅载荷下齿轮接触疲劳寿命可靠性评估方法及装置 |
CN113591340A (zh) * | 2021-06-24 | 2021-11-02 | 内蒙古工业大学 | 随机载荷下齿轮弯曲疲劳时变可靠性分析方法 |
CN113742948B (zh) * | 2021-08-23 | 2022-07-22 | 西安石油大学 | 一种超高强度抽油杆p-s-n曲线拟合新模型及方法 |
CN113742948A (zh) * | 2021-08-23 | 2021-12-03 | 西安石油大学 | 一种超高强度抽油杆p-s-n曲线拟合新模型及方法 |
CN113536490B (zh) * | 2021-08-27 | 2022-05-17 | 西安石油大学 | 一种抽油杆疲劳寿命预测方法及系统 |
CN113536490A (zh) * | 2021-08-27 | 2021-10-22 | 西安石油大学 | 一种抽油杆疲劳寿命预测方法及系统 |
CN115408755A (zh) * | 2022-09-05 | 2022-11-29 | 北京交通大学 | 一种考虑时变效应的组合梁桥动力疲劳可靠性评估方法 |
CN115408755B (zh) * | 2022-09-05 | 2024-03-22 | 北京交通大学 | 一种考虑时变效应的组合梁桥动力疲劳可靠性评估方法 |
CN116049922B (zh) * | 2023-01-12 | 2023-08-22 | 中国工程物理研究院研究生院 | 应力强度因子门槛值概率建模的金属结构疲劳可靠性预测方法 |
CN116049922A (zh) * | 2023-01-12 | 2023-05-02 | 中国工程物理研究院研究生院 | 应力强度因子门槛值概率建模的金属结构疲劳可靠性预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107103140A (zh) | 一种基于双线性累积损伤的时变疲劳可靠性分析方法 | |
CN103810328B (zh) | 一种基于混合模型的变压器维修决策方法 | |
CN105868912A (zh) | 基于数据融合的电力变压器状态评价方法及装置 | |
CN109060497A (zh) | 一种考虑低于疲劳极限加载的疲劳损伤分析方法 | |
CN103018063B (zh) | 基于Mittag-Leffler分布的桥梁随机疲劳寿命预测方法 | |
Kwon et al. | Fatigue life assessment and lifetime management of aluminum ships using life-cycle optimization | |
CN109543355A (zh) | 一种基于有限元法的应力-循环次数(s-n)的疲劳失效算法 | |
Berenguer et al. | Inspection and maintenance planning: An application of semi-Markov decision processes | |
Mouhib et al. | Experimental characterization and damage reliability analysis of central core strand extracted from steel wire rope | |
Ku et al. | Structural reliability applications in developing risk-based inspection plans for a floating production installation | |
CN111882238A (zh) | 基于云模型和eahp的门座起重机结构健康评估方法 | |
CN112613191A (zh) | 电缆健康状态评估方法、装置、计算机设备和存储介质 | |
Dissanayake et al. | Reliability index for structural health monitoring of aging bridges | |
CN111639466A (zh) | 一种基于监测数据的斜拉桥拉索使用寿命评估方法 | |
Kepka et al. | Structure service life assessment under combined loading using probability approach | |
Gill et al. | The determination procedure of the onset of the object wear-out period based on monitoring of the empirical failure intensity function | |
Roy et al. | A stochastic model for evolution of creep damage in engineering materials | |
CN113282988B (zh) | 一种大跨度桥梁上部结构状态评估方法及系统 | |
Qin et al. | Fatigue Life Scatter Modeling and Failure Probability Prediction of Multisite Damage Structural Part under Variable Amplitude Load History | |
Suleimenov et al. | Concept of developing an intelligent system for control and operational diagnostics of technological equipment condition | |
CN101482464A (zh) | 确定高温构件持久寿命可靠度的方法 | |
Weller et al. | The role of accelerated testing in reliability prediction | |
CN112255393B (zh) | 一种污水隧道结构损伤风险评级方法及装置 | |
CN117877646B (zh) | 基于不同评价指标的轮轨材料滚动接触疲劳损伤预测方法 | |
Farrow et al. | Development of a New, Correlated FEA Method of Assessing Mooring Chain Fatigue |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170829 |
|
RJ01 | Rejection of invention patent application after publication |