CN107884290B - 一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法 - Google Patents
一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法 Download PDFInfo
- Publication number
- CN107884290B CN107884290B CN201710975283.5A CN201710975283A CN107884290B CN 107884290 B CN107884290 B CN 107884290B CN 201710975283 A CN201710975283 A CN 201710975283A CN 107884290 B CN107884290 B CN 107884290B
- Authority
- CN
- China
- Prior art keywords
- crack
- influence
- under
- abrasion
- finite element
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/32—Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0001—Type of application of the stress
- G01N2203/0005—Repeated or cyclic
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0062—Crack or flaws
- G01N2203/0066—Propagation of crack
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0073—Fatigue
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开的一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法,属于材料科学与工程应用技术领域。本发明包括如下步骤:步骤一,建立接触结构在微动载荷下的有限元模型;步骤二,更新受磨损影响的在微动载荷下的有限元模型;步骤三,确定裂纹扩展速率dl/dN和裂纹扩展角度θ,并更新裂纹扩展后的受磨损影响的在微动载荷下的有限元模型;步骤四,完成考虑磨损影响下的微动疲劳裂纹扩展寿命预测;还包括步骤五:根据步骤四得到的考虑磨损影响下的微动疲劳裂纹扩展寿命预测结果,得到考虑磨损影响下的微动疲劳裂纹扩展规律,解决微动疲劳裂纹扩展领域实际工程问题。本发明能够实现考虑磨损影响下的微动疲劳裂纹扩展寿命预测,具有精度高的优点。
Description
技术领域
本发明涉及一种接触机构的微动疲劳裂纹扩展寿命预测方法,属于材料科学与工程应用技术领域。
背景技术
目前,微动疲劳问题广泛存在于机械接触配合结构中,如螺栓连接件、过盈配合件、内燃机机体和主轴承盖的接触面以及隔框、桁条和蒙皮的搭接组合结构都比较容易发生微动疲劳损伤。不同于常规疲劳,微动疲劳裂纹扩展寿命占据了大部分的疲劳寿命。因此,准确预测疲劳裂纹扩展寿命能够大大提高构件的使用容限设计。
当前针对微动疲劳裂纹扩展寿命的预测通常会忽略磨损的影响,微动磨损能够逐渐减少试件表层的材料,当裂纹处于磨损轮廓范围内,其实际裂纹长度往往会受到磨损的影响。当微动接触方式为柱-面接触时,磨损会使接触面积增加,从而改变驱使裂纹扩展的外载荷。磨损的影响在一定程度上会提高微动疲劳裂纹扩展寿命,因此在计算微动疲劳裂纹扩展寿命时将磨损行为考虑在内可以大大提高材料的损伤容限设计,充分发挥材料性能。
发明内容
本发明公开的一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法要解决的技术问题为:实现考虑磨损影响下的微动疲劳裂纹扩展寿命预测,具有精度高的优点。
本发明的目的是通过下述技术方案实现的。
本发明公开的一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法,包括如下步骤:
步骤一,建立接触结构在微动载荷下的有限元模型。
所述的接触结构在微动载荷下的有限元模型包括压头模型与试件模型,试件模型的接触表面包含至少一条初始裂纹,初始裂纹的长度li和初始裂纹角度θ根据实际情况确定。
步骤二,更新受磨损影响的在微动载荷下的有限元模型。
步骤2.1:对压头模型施加法向压力F,对试件模型一端施加循环载荷σB,通过建立的有限元模型计算接触面的接触压力分布与相对滑移值分布。
步骤2.2:通过计算所得的接触压力分布与相对滑移值分布,采用磨损模型计算磨损量体积V。
步骤2.3:为了能够通过有限元方法计算接触面每一个节点位置的磨损深度h,将步骤2.2采用磨损模型转化成计算每一个节点的磨损深度增量Δh的磨损模型。
步骤2.4:考虑到计算每一个循环的磨损量会增加工作量,需采用循环跳跃技术,在一定工作循环ΔN内磨损增量是保持不变的,ΔN为预设的工作循环增量数。因此,对计算ΔN个循环磨损深度增量的磨损模型修改成ΔhΔN(x,t)。
步骤2.5:通过步骤2.4计算所得的磨损深度增量值ΔhΔN(x,t)更新受磨损影响的有限元模型接触面节点在垂直于接触面方向上的位置,以实现接触面轮廓的改变。在对压头施加的法向压力不变的情况下,接触面轮廓的改变会改变接触面每一个节点的接触压力与接触剪切力,同时也会消减预置裂纹的实际长度,即完成更新受磨损影响的在微动载荷下的有限元模型。
作为优选,
步骤二具体实现方法为,
步骤2.1:对压头模型施加法向压力,对试件模型一端施加循环载荷,通过建立的有限元模型计算接触面的接触压力分布与相对滑移值分布。
步骤2.2所述的磨损模型优选如式(1)所示的磨损模型,
其中,V是磨损体积,K是磨损系数,P是接触压强,H是材料硬度,δ是相对滑移值。
步骤2.3:为了能够通过有限元方法计算接触面每一个节点位置的磨损深度增量Δh(x,t),需将式(1)磨损模型修改成:
Δh(x,t)=k·p(x,t)δ(x,t) (2)
其中,Δh(x,t)是每一个工作循环的磨损深度增量,p(x,t)是接触压强,δ(x,t)是相对滑移值,k代替K/H成为磨损系数,x为节点位置,t代表时间。
步骤2.4:考虑到计算每一个循环的磨损增量会增加工作量,因此需采用循环跳跃技术,即在一定工作循环ΔN内磨损增量是保持不变的。因此磨损模型最终修改成
ΔhΔN(x,t)=ΔN·k·p(x,t)δ(x,t) (3)
步骤2.5:通过步骤2.4计算所得的磨损深度增量值ΔhΔN(x,t),更新受磨损影响的有限元模型接触面节点在垂直于接触面方向上的位置,以实现接触面轮廓的改变。在对压头施加的法向压力F不变的情况下,接触面轮廓的改变会改变接触面每一个节点的接触压强P与接触剪切力τ,同时也会消减预置裂纹的实际长度l,即完成更新受磨损影响的在微动载荷下的有限元模型。
步骤三,确定裂纹扩展速率dl/dN和裂纹扩展角度θ,并更新裂纹扩展后的受磨损影响的在微动载荷下的有限元模型。
步骤3.1:确定裂纹扩展速率的表达式。
裂纹尖端应力状态可以由准则f表达,准则f的表达式如下:
其中,KI、KII和T是准则中的变量,其他为相关系数,KI代表I型应力强度因子,KII代表II型应力强度因子,T代表T应力,其他系数的表达式如下表所示
其中v是泊松比,是长度尺度参数。
准则f与裂纹扩展速率建立下列关系:dl/dN=α·(fmax)γ (5)
其中,l为裂纹长度,α、γ为疲劳系数,fmax为一个工作循环内f的最大值,fmax如果小于0代表裂纹尖端的应力场不足以使裂纹扩展,此时计算终止。
步骤3.2:确定裂纹扩展方向的表达式。
利用最大轴向应力准则建立裂纹扩展角度θ与KI、KII和T应力的关系
其中临界距离rc可以表示成
其中,σb为拉伸强度,KIc为断裂韧度。
步骤3.3:根据有限元计算KI、KII和T应力。
步骤3.3优选如下方法实现:
在有限元模型裂纹尖端处划分退化的单元以保证计算所需的精度。利用有限元方法计算裂纹尖端处真实场的J积分,在裂纹尖端处建立由点载荷产生的辅助场并计算辅助场Jaux积分。根据叠加原理,真实场合辅助场的叠加场Js积分表达为:
Js=J+Jaux+M1+M2 (8)
其中,M1和M2为真实场和辅助场的交叉项。对于各向同向材料
其中,E′=E/(1-v2)
步骤3.4:计算裂纹扩展速率和裂纹扩展角度。
在预设Nt个循环内磨损对裂尖应力状态保持不变的条件下,通过受磨损影响有限元模型计算裂纹尖端的KI、KII和T应力,结合准则f与最大轴向应力准则,计算裂纹扩展增量Δlj与裂纹扩展角度θ。所述的Nt为ΔN的倍数。
步骤3.5:通过网格重划分技术实现预置裂纹在计算方向θ与计算长度增量Δlj的扩展,得到裂纹扩展更新后的受磨损影响的在微动载荷下的有限元模型。
步骤四,通过步骤二更新Nt个工作循环内受磨损影响的在微动载荷下的有限元模型,从而改变接触面积以及初始裂纹长度li。在载荷一定的情况下接触面积的变化导致接触应力分布改变,结合初始裂纹长度li的变化从而改变了裂纹尖端的应力场分布。通过步骤三计算在步骤二更新后的微动载荷下的有限元模型计算裂纹扩展长度增量Δlj和裂纹扩展角度θ,并更新裂纹扩展后的受磨损影响的在微动载荷下的有限元模型,使步骤二更新后的裂纹在步骤三计算的裂纹扩展角度方向上实现裂纹扩展,所述的裂纹扩展长度为步骤三计算的裂纹扩展长度增量Δlj,所述的更新后的微动载荷下的有限元模型计算裂纹长度作为下一次迭代的初始裂纹长度l。通过步骤二再次更新Nt个工作循环内在步骤三更新裂纹扩展后的受磨损影响的在微动载荷下的有限元模型。通过重复上述迭代过程,直到裂纹长度l达到预设临界值,完成考虑磨损影响下的微动疲劳裂纹扩展寿命预测。
还包括步骤五:根据步骤四得到的考虑磨损影响下的微动疲劳裂纹扩展寿命预测结果,得到考虑磨损影响下的微动疲劳裂纹扩展规律,解决微动疲劳裂纹扩展领域实际工程问题。
所述的解决微动疲劳裂纹扩展领域实际工程问题包括材料损伤容限设计、材料剩余寿命预测等,通过解决上述工程问题进而充分发挥材料使用寿命。
有益效果:
1、现有技术中微动疲劳裂纹扩展寿命预测方法未考虑在裂纹扩展过程中磨损对裂纹长度的影响,以及在裂纹扩展过程中磨损对接触面应力分布的影响。本发明公开的一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法,将磨损深度的改变嵌入到每一个裂纹扩展增量的计算中,在每一个裂纹扩展增量计算中,利用磨损深度的改变来修正裂纹实际长度以及接触面轮廓,在载荷一定的情况下接触面轮廓的变化改变接触应力分布,从而影响裂尖应力状态,通过裂尖应力状态计算裂纹扩展长度增量和裂纹扩展角度,即在微动疲劳裂纹扩展寿命预测中考虑磨损影响,提高预测精度。
2、本发明公开的一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法,根据得到的考虑磨损影响下的微动疲劳裂纹扩展寿命预测结果,能够进一步得到考虑磨损影响下的微动疲劳裂纹扩展规律,解决微动疲劳裂纹扩展领域实际工程问题。
附图说明
图1为预置初始裂纹的微动有限元模型;
图2为裂纹尖端处的退化单元;
图3为有限元磨损过程说明图;
图4为考虑磨损行为的微动疲劳裂纹扩展寿命预测流程图;
图5有无磨损影响下的裂纹长度与循环数对比图。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图对本发明内容作进一步说明。
本实施例公开的一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法,具体实现步骤如下:
步骤一,建立接触结构在微动载荷下的有限元模型。
所述的接触结构在微动载荷下的有限元模型包括压头模型与试件模型(如图1),模型的材料参数弹性模量为119400MPa,试件模型的接触表面包含至少一条初始裂纹,初始裂纹的长度10微米和初始裂纹角度90度。
步骤二,更新受磨损影响的在微动载荷下的有限元模型。
步骤2.1:对压头模型施加法向压力10MPa,对试件模型一端施加循环载荷700MPa,通过建立的有限元模型计算接触面的接触压力分布与相对滑移值分布。
步骤2.2:通过磨损模型式(1)计算磨损量,
其中,V是磨损体积,K是磨损系数,P是接触压强,H是材料硬度,δ是相对滑移值。
步骤2.3:为了能够通过有限元方法计算接触面每一个节点位置的磨损深度增量Δh(x,t),需将式(1)磨损模型修改成:
Δh(x,t)=k·p(x,t)δ(x,t) (2)
其中,Δh(x,t)是每一个工作循环的磨损深度增量,p(x,t)是接触压强,δ(x,t)是相对滑移值,k代替K/H成为磨损系数取值2.75×10-8MPa,x为节点位置,t代表时间。
步骤2.4:考虑到计算每一个循环的磨损增量会增加工作量,因此需采用循环跳跃技术,即在一定工作循环ΔN内磨损增量是保持不变的,ΔN在这里取值200。因此磨损模型(符号统一)最终修改成
ΔhΔN(x,t)=ΔN·k·p(x,t)δ(x,t) (3)
步骤2.5:如图2所示通过步骤2.4计算所得的磨损深度增量值ΔhΔN(x,t)更新受磨损影响的有限元模型接触面节点在垂直于接触面方向上的位置,以实现接触面轮廓的改变。在对压头施加的法向压力F不变的情况下,接触面轮廓的改变会改变接触面每一个节点的接触压强P与接触剪切力τ,同时也会消减预置裂纹的实际长度l,即完成更新受磨损影响的在微动载荷下的有限元模型。
步骤三,确定裂纹扩展速率dl/dN和裂纹扩展角度θ,并更新裂纹扩展后的受磨损影响的在微动载荷下的有限元模型。
步骤3.1:确定裂纹扩展速率的表达式。
裂纹尖端应力状态可以由准则f表达,准则f的表达式如下:
其中,KI、KII和T是准则中的变量,其他为相关系数,KI代表I型应力强度因子,KII代表II型应力强度因子,T代表T应力,其他系数的表达式如下表所示
其中v取值0.3,取值1。
准则f与裂纹扩展速率建立下列关系:dl/dN=α·(fmax)γ (5)其中,l为裂纹长度,α、γ分别取值9×10-10和1,fmax为一个工作循环内f的最大值。
步骤3.2:确定裂纹扩展方向的表达式。
利用最大轴向应力准则建立裂纹扩展角度θ与KI、KII和T应力的关系
其中临界距离rc可以表示成
其中,σb为拉伸强度,KIc为断裂韧度。
步骤3.3:根据有限元计算KI、KII和T应力。
在有限元模型裂纹尖端处划分退化的单元(如图3)以保证计算所需的精度。利用有限元方法计算裂纹尖端处真实场的J积分,在裂纹尖端处建立由点载荷产生的辅助场并计算辅助场Jaux积分。根据叠加原理,真实场合辅助场的叠加场Js积分表达为:
Js=J+Jaux+M1+M2 (8)其中,M1和M2为真实场和辅助场的交叉项。对于各向同向材料
其中,E′=E/(1-v2)
步骤3.4:计算裂纹扩展速率和裂纹扩展角度。
在预设Nt=2000个循环内磨损对裂尖应力状态保持不变的条件下,通过受磨损影响有限元模型计算裂纹尖端的KI、KII和T应力,结合准则f与最大轴向应力准则,计算裂纹扩展增量Δlj与裂纹扩展角度θ。
步骤3.5:通过网格重划分技术实现预置裂纹在计算方向θ与计算长度增量Δlj的扩展,得到裂纹扩展更新后的受磨损影响的在微动载荷下的有限元模型。
步骤四,如图4所示,通过步骤二更新2000个工作循环内受磨损影响的在微动载荷下的有限元模型,从而改变接触面积以及初始裂纹长度li。在载荷一定的情况下接触面积的变化导致接触应力分布改变,结合初始裂纹长度li的变化从而改变了裂纹尖端的应力场分布。通过步骤三计算在步骤二更新后的微动载荷下的有限元模型计算裂纹扩展长度增量Δlj和裂纹扩展角度θ,并更新裂纹扩展后的受磨损影响的在微动载荷下的有限元模型,使步骤二更新后的裂纹在步骤三计算的裂纹扩展角度方向上实现裂纹扩展,所述的裂纹扩展长度为步骤三计算的裂纹扩展长度增量Δlj,所述的更新后的微动载荷下的有限元模型计算裂纹长度作为下一次迭代的初始裂纹长度l。通过步骤二再次更新Nt个工作循环内在步骤三更新裂纹扩展后的受磨损影响的在微动载荷下的有限元模型。通过重复上述迭代过程,直到裂纹长度l达到预设临界值20微米,完成考虑磨损影响下的微动疲劳裂纹扩展寿命预测。图5为考虑磨损与不考虑磨损的裂纹扩展寿命对比图,可以看到考虑磨损的寿命预测值要大于不考虑磨损的预测值。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法,其特征在于:包括如下步骤:
步骤一,建立接触结构在微动载荷下的有限元模型;
所述的接触结构在微动载荷下的有限元模型包括压头模型与试件模型,试件模型的接触表面包含至少一条初始裂纹,初始裂纹的长度li和初始裂纹角度θ根据实际情况确定;
步骤二,更新受磨损影响的在微动载荷下的有限元模型;
步骤2.1:对压头模型施加法向压力F,对试件模型一端施加循环载荷σB,通过建立的有限元模型计算接触面的接触压力分布与相对滑移值分布;
步骤2.2:通过计算所得的接触压力分布与相对滑移值分布,采用磨损模型计算磨损量体积V;
步骤2.3:为了能够通过有限元方法计算接触面每一个节点位置的磨损深度h,将步骤2.2采用磨损模型转化成计算每一个节点的磨损深度增量Δh(x,t)的磨损模型;
步骤2.4:考虑到计算每一个循环的磨损量会增加工作量,需采用循环跳跃技术,在一定工作循环ΔN内磨损增量是保持不变的,ΔN为预设的工作循环增量数;因此,对计算ΔN个循环磨损深度增量的磨损模型修改成ΔhΔN(x,t),x为节点位置,t代表时间;
步骤2.5:通过步骤2.4计算所得的磨损深度增量值ΔhΔN(x,t)更新受磨损影响的有限元模型接触面节点在垂直于接触面方向上的位置,以实现接触面轮廓的改变;在对压头施加的法向压力不变的情况下,接触面轮廓的改变会改变接触面每一个节点的接触压力与接触剪切力,同时也会消减预置裂纹的实际长度,即完成更新受磨损影响的在微动载荷下的有限元模型;
步骤三,确定裂纹扩展速率dl/dN和裂纹扩展角度θ,并更新裂纹扩展后的受磨损影响的在微动载荷下的有限元模型;
步骤3.1:确定裂纹扩展速率的表达式;
裂纹尖端应力状态由准则f表达,准则f的表达式如下:
其中,KI、KII和T是准则中的变量,其他为相关系数,KI代表I型应力强度因子,KII代表II型应力强度因子,T代表T应力,其他系数的表达式如下表所示
其中v是泊松比,是长度尺度参数;
准则f与裂纹扩展速率建立下列关系:dl/dN=α·(fmax)γ (5)
其中,l为裂纹长度,α、γ为疲劳系数,fmax为一个工作循环内f的最大值,fmax如果小于0代表裂纹尖端的应力场不足以使裂纹扩展,此时计算终止;
步骤3.2:确定裂纹扩展方向的表达式;
利用最大轴向应力准则建立裂纹扩展角度θ与KI、KII和T应力的关系
其中临界距离rc表示成
其中,σb为拉伸强度,KIc为断裂韧度;
步骤3.3:根据有限元计算KI、KII和T应力;
步骤3.4:计算裂纹扩展速率和裂纹扩展角度;
在预设Nt个循环内磨损对裂尖应力状态保持不变的条件下,通过受磨损影响有限元模型计算裂纹尖端的KI、KII和T应力,结合准则f与最大轴向应力准则,计算裂纹扩展长度增量Δlj与裂纹扩展角度θ;所述的Nt为ΔN的倍数;
步骤3.5:通过网格重划分技术实现预置裂纹在计算裂纹扩展角度θ与计算裂纹扩展长度增量Δlj的扩展,得到裂纹扩展更新后的受磨损影响的在微动载荷下的有限元模型;
步骤四,通过步骤二更新Nt个工作循环内受磨损影响的在微动载荷下的有限元模型,从而改变接触面积以及初始裂纹长度li;在载荷一定的情况下接触面积的变化导致接触应力分布改变,结合初始裂纹长度li的变化从而改变了裂纹尖端的应力场分布;通过步骤三计算在步骤二更新后的微动载荷下的有限元模型计算裂纹扩展长度增量Δlj和裂纹扩展角度θ,并更新裂纹扩展后的受磨损影响的在微动载荷下的有限元模型,使步骤二更新后的裂纹在步骤三计算的裂纹扩展角度方向上实现裂纹扩展,所述的裂纹扩展长度为步骤三计算的裂纹扩展长度增量Δlj,所述的更新后的微动载荷下的有限元模型计算裂纹长度作为下一次迭代的初始裂纹长度l;通过步骤二再次更新Nt个工作循环内在步骤三更新裂纹扩展后的受磨损影响的在微动载荷下的有限元模型;通过重复上述迭代过程,直到裂纹长度l达到预设临界值,完成考虑磨损影响下的微动疲劳裂纹扩展寿命预测。
2.如权利要求1所述的一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法,其特征在于:还包括步骤五:根据步骤四得到的考虑磨损影响下的微动疲劳裂纹扩展寿命预测结果,得到考虑磨损影响下的微动疲劳裂纹扩展规律,解决微动疲劳裂纹扩展领域实际工程问题。
3.如权利要求1或2所述的一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法,其特征在于:步骤二具体实现方法为,
步骤2.1:对压头模型施加法向压力,对试件模型一端施加循环载荷,通过建立的有限元模型计算接触面的接触压力分布与相对滑移值分布;
步骤2.2所述的磨损模型如式(1)所示的磨损模型,
其中,V是磨损体积,K是磨损系数,P是接触压强,H是材料硬度,δ是相对滑移值;
步骤2.3:为了能够通过有限元方法计算接触面每一个节点位置的磨损深度增量Δh(x,t),需将式(1)磨损模型修改成:
Δh(x,t)=k·p(x,t)δ(x,t) (2)
其中,Δh(x,t)是每一个工作循环的磨损深度增量,p(x,t)是接触压强,δ(x,t)是相对滑移值,k代替K/H成为磨损系数,x为节点位置,t代表时间;
步骤2.4:考虑到计算每一个循环的磨损增量会增加工作量,因此需采用循环跳跃技术,即在一定工作循环ΔN内磨损增量是保持不变的;因此磨损模型最终修改成
ΔhΔN(x,t)=ΔN·k·p(x,t)δ(x,t) (3)
步骤2.5:通过步骤2.4计算所得的磨损深度增量值ΔhΔN(x,t),更新受磨损影响的有限元模型接触面节点在垂直于接触面方向上的位置,以实现接触面轮廓的改变;在对压头施加的法向压力F不变的情况下,接触面轮廓的改变会改变接触面每一个节点的接触压强P与接触剪切力τ,同时也会消减预置裂纹的实际长度l,即完成更新受磨损影响的在微动载荷下的有限元模型。
4.如权利要求3所述的一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法,其特征在于:步骤3.3选如下方法实现:
在有限元模型裂纹尖端处划分退化的单元以保证计算所需的精度;利用有限元方法计算裂纹尖端处真实场的J积分,在裂纹尖端处建立由点载荷产生的辅助场并计算辅助场Jaux积分;根据叠加原理,真实场合辅助场的叠加场Js积分表达为:
Js=J+Jaux+M1+M2 (8)
其中,M1和M2为真实场和辅助场的交叉项;对于各向同向材料
其中,E′=E/(1-v2)。
5.如权利要求2所述的一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法,其特征在于:所述的解决微动疲劳裂纹扩展领域实际工程问题包括材料损伤容限设计、材料剩余寿命预测,通过解决上述工程问题进而充分发挥材料使用寿命。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710975283.5A CN107884290B (zh) | 2017-10-16 | 2017-10-16 | 一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710975283.5A CN107884290B (zh) | 2017-10-16 | 2017-10-16 | 一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107884290A CN107884290A (zh) | 2018-04-06 |
CN107884290B true CN107884290B (zh) | 2019-11-19 |
Family
ID=61781827
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710975283.5A Active CN107884290B (zh) | 2017-10-16 | 2017-10-16 | 一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107884290B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732035B (zh) * | 2018-04-28 | 2020-07-28 | 南京航空航天大学 | 一种榫连接结构的高温微动疲劳寿命测试方法 |
CN108918261B (zh) * | 2018-05-08 | 2020-12-01 | 中国科学院金属研究所 | 一种用少量试验测定材料构件疲劳寿命规律的方法 |
CN108920883B (zh) * | 2018-08-20 | 2020-08-07 | 南京航空航天大学 | 一种基于疲劳、氧化交互作用的热成形模具使用寿命预测方法 |
CN109408900B (zh) * | 2018-09-29 | 2022-04-05 | 西北工业大学 | 镍基单晶合金涡轮叶片榫头微动疲劳寿命预测方法 |
CN109726471B (zh) * | 2018-12-28 | 2022-11-29 | 合肥通用机械研究院有限公司 | 一种冷凝器在微动磨损作用下的管束疲劳失效预测方法 |
CN110426194B (zh) * | 2019-08-23 | 2020-12-25 | 中国兵器科学研究院 | 考虑磨损和疲劳交互作用的机械部件寿命试验方法 |
CN110823735B (zh) * | 2019-11-22 | 2021-06-18 | 吉林大学 | 基于非均匀应力构筑表面硬度差异性的方法 |
CN110849753B (zh) * | 2019-12-09 | 2020-09-29 | 大连理工大学 | 一种基于微观划痕的金属材料疲劳强度预测方法 |
CN110929344B (zh) * | 2019-12-26 | 2024-02-13 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机结构疲劳裂纹扩展方向的预测方法和装置 |
CN111737901A (zh) * | 2020-06-23 | 2020-10-02 | 石家庄铁道大学 | 一种刀盘疲劳寿命预测方法及其应用 |
CN112182927B (zh) * | 2020-09-15 | 2022-06-07 | 西北工业大学 | 一种模拟高温合金线性摩擦焊接头裂纹相互影响的方法 |
CN113514351B (zh) * | 2020-12-23 | 2024-04-19 | 中国特种设备检测研究院 | 一种考虑预应力重分布的疲劳裂纹扩展行为预测方法 |
CN116356289B (zh) * | 2023-02-18 | 2023-11-17 | 阳光中科(福建)能源股份有限公司 | 一种石墨舟卡点、石墨舟及镀膜方法 |
CN116842778B (zh) * | 2023-05-29 | 2024-01-26 | 天津大学 | 基于紧凑拉伸-剪切试样的疲劳裂纹扩展方向、长度及速率计算方法 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1862256A (zh) * | 2005-05-10 | 2006-11-15 | 上海市市政工程管理处 | 一种钢结构桥梁寿命的断裂力学测量方法 |
CN1869640A (zh) * | 2006-05-24 | 2006-11-29 | 浙江大学 | 一种疲劳裂纹扩展测试的方法 |
CN103217346A (zh) * | 2013-04-01 | 2013-07-24 | 北京航空航天大学 | 一种测量材料高温蠕变裂纹扩展门槛值的方法 |
CN103344511A (zh) * | 2013-07-01 | 2013-10-09 | 北京理工大学 | 一种大功率柴油机铝合金活塞的蠕变-疲劳寿命预测方法 |
CN103344515A (zh) * | 2013-07-05 | 2013-10-09 | 北京航空航天大学 | 一种基于局部应力应变法的低周疲劳和高强度冲击耦合的损伤计算方法 |
CN103616179A (zh) * | 2013-12-05 | 2014-03-05 | 广西大学 | 一种基于缺陷建模的传动齿轮疲劳寿命评估方法 |
CN104034576A (zh) * | 2014-06-16 | 2014-09-10 | 中国飞机强度研究所 | 一种飞机结构裂纹扩展试验用载荷谱简化方法 |
CN104122137A (zh) * | 2014-05-19 | 2014-10-29 | 合肥通用机械研究院 | 一种基于寿命的超高压容器疲劳强度设计方法 |
CN104535438A (zh) * | 2014-12-31 | 2015-04-22 | 北京航空航天大学 | 一种试验件高温高低周复合疲劳裂纹扩展试验系统及测量方法 |
CN104792632A (zh) * | 2015-03-26 | 2015-07-22 | 北京理工大学 | 一种气缸盖模拟件疲劳试验方法 |
CN104792633A (zh) * | 2015-04-17 | 2015-07-22 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种飞机机体裂纹扩展寿命预测方法 |
CN105447316A (zh) * | 2015-12-01 | 2016-03-30 | 北京航空航天大学 | 一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法 |
CN105823677A (zh) * | 2016-06-11 | 2016-08-03 | 贵州大学 | 多载荷疲劳寿命测试过程中磨损和旋转的控制方法及装置 |
CN106194426A (zh) * | 2016-07-26 | 2016-12-07 | 北京理工大学 | 一种用于自由活塞直线发电系统切换过程的控制方法 |
CN106568660A (zh) * | 2016-10-14 | 2017-04-19 | 中国人民解放军海军航空工程学院青岛校区 | 一种复合材料胶接修补结构的剩余疲劳寿命的预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4486420B2 (ja) * | 2004-06-23 | 2010-06-23 | 株式会社ブリヂストン | タイヤ経時変化予測方法、装置、プログラム及び媒体 |
-
2017
- 2017-10-16 CN CN201710975283.5A patent/CN107884290B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1862256A (zh) * | 2005-05-10 | 2006-11-15 | 上海市市政工程管理处 | 一种钢结构桥梁寿命的断裂力学测量方法 |
CN1869640A (zh) * | 2006-05-24 | 2006-11-29 | 浙江大学 | 一种疲劳裂纹扩展测试的方法 |
CN103217346A (zh) * | 2013-04-01 | 2013-07-24 | 北京航空航天大学 | 一种测量材料高温蠕变裂纹扩展门槛值的方法 |
CN103344511A (zh) * | 2013-07-01 | 2013-10-09 | 北京理工大学 | 一种大功率柴油机铝合金活塞的蠕变-疲劳寿命预测方法 |
CN103344515A (zh) * | 2013-07-05 | 2013-10-09 | 北京航空航天大学 | 一种基于局部应力应变法的低周疲劳和高强度冲击耦合的损伤计算方法 |
CN103616179A (zh) * | 2013-12-05 | 2014-03-05 | 广西大学 | 一种基于缺陷建模的传动齿轮疲劳寿命评估方法 |
CN104122137A (zh) * | 2014-05-19 | 2014-10-29 | 合肥通用机械研究院 | 一种基于寿命的超高压容器疲劳强度设计方法 |
CN104034576A (zh) * | 2014-06-16 | 2014-09-10 | 中国飞机强度研究所 | 一种飞机结构裂纹扩展试验用载荷谱简化方法 |
CN104535438A (zh) * | 2014-12-31 | 2015-04-22 | 北京航空航天大学 | 一种试验件高温高低周复合疲劳裂纹扩展试验系统及测量方法 |
CN104792632A (zh) * | 2015-03-26 | 2015-07-22 | 北京理工大学 | 一种气缸盖模拟件疲劳试验方法 |
CN104792633A (zh) * | 2015-04-17 | 2015-07-22 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种飞机机体裂纹扩展寿命预测方法 |
CN105447316A (zh) * | 2015-12-01 | 2016-03-30 | 北京航空航天大学 | 一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法 |
CN105823677A (zh) * | 2016-06-11 | 2016-08-03 | 贵州大学 | 多载荷疲劳寿命测试过程中磨损和旋转的控制方法及装置 |
CN106194426A (zh) * | 2016-07-26 | 2016-12-07 | 北京理工大学 | 一种用于自由活塞直线发电系统切换过程的控制方法 |
CN106568660A (zh) * | 2016-10-14 | 2017-04-19 | 中国人民解放军海军航空工程学院青岛校区 | 一种复合材料胶接修补结构的剩余疲劳寿命的预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107884290A (zh) | 2018-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107884290B (zh) | 一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法 | |
CN104794365B (zh) | 一种基于数学模型预测锚杆极限承载力的计算方法 | |
CN106886663B (zh) | 齿轮弯曲疲劳寿命预测方法及装置 | |
Chatterjee et al. | Numerical simulations of pipe–soil interaction during large lateral movements on clay | |
CN100549660C (zh) | 基于嵌入式有限元建模的三维裂纹扩展模拟方法 | |
CN107506535A (zh) | 一种基于临界应变损伤参量的多轴疲劳寿命预测方法 | |
Herath et al. | Structural strength and laminate optimization of self-twisting composite hydrofoils using a Genetic Algorithm | |
CN111090953A (zh) | 基于材料损伤理论与磨损耦合的接触疲劳失效分析方法 | |
CN110083984A (zh) | 一种涡轮叶片区间损伤容限分析方法 | |
CN109115479A (zh) | 一种基于临界面的涡轮轴多轴疲劳寿命预测方法 | |
CN115083541A (zh) | 基于icvefg法计算正交各向异性材料应力强度因子的方法 | |
CN112651061A (zh) | 一种利用Dynamo建立多跨桥梁上部结构BIM模型的方法 | |
Bittencourt et al. | Simulation of arbitrary, cohesive crack propagation | |
CN110543680B (zh) | 一种局部地基沉陷引起地埋管线变形及内力计算方法 | |
CN110375913B (zh) | 一种盾构法隧道非预埋式结构健康监测方法 | |
Manners | Plastic deformation of a sinusoidal surface | |
Vasovic et al. | Determination of stress intensity factors in low pressure turbine rotor discs | |
CN103321273A (zh) | 一种松土器齿尖、推土机及松土器齿尖模型生成方法 | |
CN112765858B (zh) | 一种考虑间隙影响的复合材料导管螺旋桨优化设计方法 | |
CN109101711A (zh) | 一种基于拓扑优化理论的框梁结构设计方法 | |
Huang et al. | Fluid-structure hydroelastic analysis and hydrodynamic cavitation experiments of composite propeller | |
Roberts et al. | SGBEM modeling of fatigue crack growth in particulate composites | |
CN114969630A (zh) | 预测法向载荷作用下粗糙分形表面摩擦系数的方法及装置 | |
CN112948993A (zh) | 一种谱载荷下考虑磨损影响的微动疲劳寿命预测方法 | |
JP2002310867A (ja) | フレッテイング疲労評価方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |