基于非线性疲劳累积损伤机理退化-冲击模型的建模方法技术

本发明专利技术提供一种基于非线性疲劳累积损伤机理退化‑冲击模型的建模方法，步骤如下：一：恒幅载荷下的疲劳退化模型；二：冲击载荷模型；三：冲击过载条件下的迟滞效应模型；四：基于分段确定性马尔科夫过程(PDMP)的迟滞效应相关模型；通过以上步骤，本发明专利技术模拟实际使用环境，还考虑到了偶发情况下造成的偶然冲击过程，具有真实性；这对于实际工程应用中，对裂纹扩展会有新的认识和理解，有着实际意义；本发明专利技术应用了科学的方法将疲劳退化过程和随机冲击过程结合起来，给出了明确的数学表达，用以计算样件的可靠性，不仅有学术价值，同样对于工程上也为其提供了一条解决该问题的途径，具有推广应用价值。

Modeling method of degradation impact model based on nonlinear fatigue cumulative damage mechanism

The invention provides a modeling method based on the nonlinear fatigue cumulative damage mechanism. The following steps are as follows: the fatigue degradation model under constant amplitude loading; two: impact load model; three: hysteresis effect model under impact overload conditions; four: Based on the piecewise deterministic Markov process (PDMP) hysteresis Through the above steps, the invention simulates the actual use of the environment and takes into account the accidental impact process caused by occasional cases, which is true. In practical engineering applications, there will be a new understanding and understanding of crack growth and practical significance in the practical engineering application; the invention applies a scientific method to fatigue. Combined with the process of degradation and random impact, a clear mathematical expression is given to calculate the reliability of the sample. It is not only of academic value, but also provides a way to solve the problem in engineering. It has the value of popularization and application.

【技术实现步骤摘要】
基于非线性疲劳累积损伤机理退化-冲击模型的建模方法
本专利技术提供一种基于非线性疲劳累积损伤机理退化-冲击模型的建模方法，它涉及一种随机退化过程和一种随机冲击过程共同作用时的疲劳累积损伤计算方法，它是一种非线性疲劳损伤叠加模型，针对冲击过大时对退化过程造成迟滞和退化过程同时会影响冲击损伤大小的问题，提出了退化和冲击相互作用的非线性疲劳损伤建模方法。这属于可靠性
中的可靠性建模问题。
技术介绍
机身是飞机最重要的部件之一，通常采用损伤容限的概念设计。在服役过程中，机身通常受到由气流引起的疲劳载荷和由湍流引起的偶发冲击而产生裂纹。这种现实载荷条件可以被转化为疲劳加载序列和一些尖峰过载组合的科学问题。研究人员通常用疲劳退化过程和随机冲击过程来模拟这种实际的物理过程。针对这样的现实问题，在国内外已有的研究中，出现了这样的两类方向：1.将疲劳退化过程和随机冲击过程看成是两个相互独立的过程，并非考虑了两者相关。2.一些研究者在某种程度上考虑了退化和冲击两个过程之间的相关性，但主要是针对冲击过大造成系统直接崩溃和冲击不大造成损伤增量进行建模研究。并非考虑了冲击和退化二者之间的相互作用。因此，考虑随机退化过程和随机冲击过程的相互作用，并提出能够计算该种相互作用的建模方法，是亟待解决的问题。
技术实现思路
本专利技术的目的在于提供一种基于疲劳损伤机理随机退化和随机冲击模型的建模方法，它用以解决受到退化和冲击作用时的裂纹扩展问题，考虑冲击载荷过载时造成退化过程迟滞和退化过程影响冲击损伤大小，建立非线性疲劳损伤模型，建立起两者之间的相互作用模型，以解决裂纹扩展整个过程的算法问题。本专利技术是采用以下技术方案实现的，本专利技术是一种基于疲劳损伤机理随机退化和随机冲击模型的建模方法，其步骤如下：步骤一：恒幅载荷下的疲劳退化模型在疲劳损伤阶段，采用巴黎(Paris)模型描述恒幅载荷下的裂纹扩展；Paris模型在当前的裂纹扩展研究中被广泛使用，其裂纹长度a的微分表达式内容如下：Δσ＝σmax-σmin(3)校正因子C受到材料的分散性影响，其被假定为一个服从正态分布的随机变量，即Y是几何因子，几何因子在特殊情况下得以简化，计算一块铝板、铁板中心孔裂纹时，几何因子Y等于1.12；计算一块铝板、铁板表面裂纹时，几何因子Y的值近似为步骤二：冲击载荷模型实际上，恒幅载荷伴随着随机冲击；裂纹增长在恒幅载荷和冲击载荷下示意图如图2所示；在图2中，TN是第N个冲击的出现时间，KN是冲击强度；在本专利技术中，所有用来计算迟滞效应的冲击称为过载；冲击出现的次数和冲击的大小均是需要给定某种随机分布；在第一次冲击到达之前，裂纹的扩展只取决于恒幅载荷；当第N个冲击出现时，裂纹在过载引起的塑性区内增长，其增长速率很低，直至塑性区消失；后来，裂纹增长率恢复到正常在第(N+1)个冲击出现之前；冲击过程中相当于一个跳跃过程；恒幅载荷下的退化过程是一个随时间不断更新的确定过程；X(n)被假定为退化和冲击过程中裂纹扩展的状态，表达式如下：X(n)＝Xn,n∈[TN,TN+1)(5)步骤三：冲击过载条件下的迟滞效应模型在实际中，冲击荷载没有理想的恒幅载荷，而是带有随机性的环境冲击；对于某些材料，如铝和铁，过载会延迟表面裂纹扩展，而不是加速其增长；这种延迟是由冲击过载造成的迟滞效应；维林博格模型(即Willenborg模型)描述了这种过载迟滞现象；在冲击载荷作用下，裂纹尖端附近出现一个大的单调塑性区，裂纹扩展速率达到最小值；在下一个载荷周期中，当塑性区超过上边界值时，迟滞效应逐渐减小，然后完全消失；在迟滞区域中，疲劳裂纹扩展可以用一下一公式进行描述：σcomp＝σreq-σmax(8)σreq是无迟滞效应时的应力,σreq的表达可通过欧文(Irwin)函数和几何准则推导：ρreq＝ρres(9)ρres＝ρol-(a(n)-aol)(10)基于公式(9)～(12),σreq的最终推导结果为：ρreq是由σreq产生。α是Irwin函数的系数，对于平面板的裂纹扩展问题，其取值为1；迟滞效应的持续时间取决于塑性区的相应大小；从公式(7)和(8)，可以找到持续时间，当其应力满足以下关系时，σmin＞σreq-σmax(14)即，当满足等式(13)的条件时，过载迟滞效应结束；通过公式(12)和(13)中的应力强度来判断迟滞区域大小可以转换成通过裂纹长度来判断迟滞区域大小，其转换后的公式表达为：从疲劳模型的公式(1)和(2)可知，未来的裂纹扩展速率取决于当前裂纹长度；因此，传统的线性累积损伤法则不再适用于计算该迟滞区域的裂缝长度；式中，r＝σmin/σmax；推导出的公式(15)能够判定以下的两种情况：能够判断每个冲击到来时，是否会发生迟滞效应；当时发生迟滞效应；能够判断迟滞效应的结束；当累积裂纹长度满足公式(15)时，迟滞效应结束；步骤四：基于分段确定性马尔科夫过程(PDMP)的迟滞效应相关模型该步骤建立了由过载引起的迟滞裂纹扩展模型；在迟滞过程中，裂纹扩展的路径受冲击过程中许多随机因素的影响，如冲击到达时间、冲击尺寸、累积裂纹长度以及材料的性能等；Willenborg模型是用来描述随机冲击载荷和裂纹扩展速率之间的关系，而Paris模型描述了不存在迟滞效应情况下的确定性过程，是恒幅载荷加载下的裂纹扩展；分段确定性马尔科夫过程是可以将随机冲击过程和确定性裂纹扩展模型联系在一起；一般来说，分段确定性马尔科夫过程很适合描述确定性过程加上随机时间：确定性过程是指疲劳退化，随机事件是随机冲击。因此，分段确定性马尔科夫过程被用来描述疲劳退化和随机冲击之间的相互作用关系。分段确定性马尔可夫过程具有三个特征：①确定性过程②跳跃过程(XN,TN,N)，N≥0；③状态转移，疲劳退化过程和随机冲击过程之间的状态转移；分段确定性马尔科夫过程具体算法实施可以通过以下三步进行：1.选择行为向量构造了一个具有环境和应力变化的多维行为向量Z(n)来描述试样的裂纹状态：2.推导流函数这一步的目的是找出两个随机冲击之间的确定路径；假设两个随机冲击的发生时间是TN和TN+1，而流函数表示为：式中，η(TN)表示第N个冲击发生时已积累的裂纹长度，迟滞效应发生的边界条件是：3.推导跳跃过程这一步包括跳跃时间和冲击载荷大小的推导。冲击的大小服从正态分布其均值和标准方差为μol和δol；两个相邻随机冲击之间的间隔时间服从强度为λ的指数分布，TN+1-TN～E(λ)(19)若系统不失效，必须满足其最大应力强度因子Kmax不能超过阈值Kthreshold，并且累积的裂纹长度a(n)不能唱过最大阈值athreshold。可靠度表达式为：R(n)＝P{Kmax(n)≤Kthreshold∩a(n)≤athreshold}(20)通过以上步骤，首先，可以看出该模型考虑了横幅载荷下的疲劳退化和突发性的冲击载荷，全面分析了整个寿命周期经历过程；其次，通过考虑冲击过载情况下造成裂纹扩展发生迟滞效应，扩展速率减缓，推导出发生迟滞效应和结束迟滞效应的边界条件；最后，提出了分段确定性马尔科夫过程，来描述疲劳退化和随机冲击载荷相结合的过程，考虑到了两者之间的相互作用，建立了耦合模型，给出了可靠度的计算方法；基于以上三点明显的优势，本专利技术对于裂纹扩展相关本文档来自技高网...

【技术保护点】
1.基于非线性疲劳累积损伤机理退化‑冲击模型的建模方法，其特征在于：其步骤如下：步骤一：恒幅载荷下的疲劳退化模型在疲劳损伤阶段，采用巴黎即Paris模型描述恒幅载荷下的裂纹扩展；该Paris模型在当前的裂纹扩展研究中被广泛使用，其裂纹长度a的微分表达式

【技术特征摘要】
1.基于非线性疲劳累积损伤机理退化-冲击模型的建模方法，其特征在于：其步骤如下：步骤一：恒幅载荷下的疲劳退化模型在疲劳损伤阶段，采用巴黎即Paris模型描述恒幅载荷下的裂纹扩展；该Paris模型在当前的裂纹扩展研究中被广泛使用，其裂纹长度a的微分表达式内容如下：Δσ＝σmax-σmin(3)校正因子C受到材料的分散性影响，其被假定为一个服从正态分布的随机变量，即Y是几何因子，几何因子在特殊情况下得以简化，计算一块铝板、铁板中心孔裂纹时，几何因子Y等于1.12；计算一块铝板、铁板表面裂纹时，几何因子Y的值近似为步骤二：冲击载荷模型实际上，恒幅载荷伴随着随机冲击；在本发明中，所有用来计算迟滞效应的冲击称为过载；冲击出现的次数和冲击的大小均是需要给定一种随机分布；在第一次冲击到达之前，裂纹的扩展只取决于恒幅载荷；当第N个冲击出现时，裂纹在过载引起的塑性区内增长，其增长速率很低，直至塑性区消失；后来，裂纹增长率恢复到正常在第(N+1)个冲击出现之前；冲击过程中相当于一个跳跃过程；恒幅载荷下的退化过程是一个随时间不断更新的确定过程；X(n)被假定为退化和冲击过程中裂纹扩展的状态，表达式如下：X(n)＝Xn,n∈[TN,TN+1)步骤三：冲击过载条件下的迟滞效应模型在实际中，冲击荷载没有理想的恒幅载荷，而是带有随机性的环境冲击；对于一些材料，如铝和铁，过载会延迟表面裂纹扩展，而不是加速其增长；这种延迟是由冲击过载造成的迟滞效应；维林博格模型即Willenborg模型描述了这种过载迟滞现象；在冲击载荷作用下，裂纹尖端附近出现一个大的单调塑性区，裂纹扩展速率达到最小值；在下一个载荷周期中，当塑性区超过上边界值时，迟滞效应逐渐减小，然后完全消失；在迟滞区域中，疲劳裂纹扩展能用一下一公式进行描述：σcomp＝σreq-σmax(8)σreq是无迟滞效应时的应力,σreq的表达通过欧文即Irwin函数和几何准则推导：ρreq＝ρres(9)ρres＝ρol-(a(n)-aol)(10)基于公式(9)～(12),σreq的最终推导结果为：ρreq是由σreq产生；α是Irwin函数的系数，对于平面板的裂纹扩展问题，其取值为1；迟滞效应的持续时间取决于塑性区的相应大小；从公式(7)和(8)，能找到持续时间，当其应力满足以下关系时，σmin＞σreq-σmax(14)即，当满足等式(13)的条件时，过载迟滞效应结束；通过公式(12)和(13)中的应力强度来判断迟滞区域大小能转换成通过裂纹长度来判断迟滞区域大小，其转换后的公式表达为：从疲劳模型的公式(1)和(2)能知，未来的裂纹扩展速率取决于当前裂纹长度；因此，传统的线性累积损伤法则不再适用于计算该迟滞区域的裂缝长度；式中，r＝σmin/σmax；推导出的公式(15)能够判定以下的两种情况：能够判断每个冲击到来时，是否会发生迟滞效应；当时发生迟滞效应；能够判断迟滞效应的结束；当累积裂纹长度满足公式(15)时，迟滞效应结束；步骤四：基于分段确定性马尔科夫过程(PDMP)的迟滞效应相关模型该步骤建立了由过载引起的迟滞裂纹扩展模型；在迟滞过程中，裂纹扩展的路径受冲击过程中许多随机因素的影响，如冲击到达时间、冲击尺寸、累积裂纹长度以及材料的性能等；Willenborg模...

【专利技术属性】
技术研发人员：许丹，何娇兰，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京,11