CN109885874A - 一种基于abaqus的多轴蠕变疲劳预测方法 - Google Patents
一种基于abaqus的多轴蠕变疲劳预测方法 Download PDFInfo
- Publication number
- CN109885874A CN109885874A CN201910026871.3A CN201910026871A CN109885874A CN 109885874 A CN109885874 A CN 109885874A CN 201910026871 A CN201910026871 A CN 201910026871A CN 109885874 A CN109885874 A CN 109885874A
- Authority
- CN
- China
- Prior art keywords
- fatigue
- creep
- stress
- strain
- tensor
- 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
- 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/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/02—Knowledge representation; Symbolic representation
-
- 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/0071—Creep
-
- 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
-
- 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/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0214—Calculations a priori without experimental data
-
- 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
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Immunology (AREA)
- Artificial Intelligence (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Chemical & Material Sciences (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于ABAQUS的多轴蠕变疲劳预测方法,其包括步骤:S1:建立ABAQUS有限元模型,通过用户子程序UMAT定义待测试材料的粘塑性本构方程;S2:确定粘塑性本构方程所需的模型参数;S3:建立待测试材料的多轴应力应变状态的疲劳损伤计算模型和蠕变损伤计算模型;S4:建立多轴应力应变状态下的ABAQUS有限元模型,基于所定义的粘塑性本构方程和模型参数,计算得到每个循环周次的应力应变张量;S5:通过用户子程序USDFLD计算等效应力和等效塑性应变,基于疲劳损伤计算模型和蠕变损伤计算模型并结合应力应变张量,通过线性累计损伤准则叠加每个循环周次的疲劳损伤和蠕变损伤,得到裂纹萌生寿命。
Description
技术领域
本发明涉及数值模拟领域,尤其涉及一种多轴蠕变疲劳预测方法。
背景技术
随着高温旋转部件大型化、高性能、长寿命的需求日益提升,包含几何不连续的结构在复杂苛刻环境下的结构完整性评估成为亟需解决的关键技术瓶颈之一。由几何不连续和复杂加载历史导致的多轴应力应变状态不可避免地成为限制了此类部件的使用寿命。近些年,有限元软件的发展可以很好地满足人们对于复杂应力应变行为的理解以及提供了该状态下精确寿命预测的可行性。
现今针对复杂结构的蠕变疲劳分析及寿命预测的主要可以分为三类。第一类是近年基于微观分析方法的晶体塑性有限元理论,其通过累积塑性滑移带和存储的能量耗散等疲劳指示因子表征蠕变疲劳演化过程。第二类是能够描述包含裂纹萌生和扩展阶段的连续损伤力学理论,其通过损伤变量引入统一的蠕变疲劳本构的方式描述材料在循环载荷下损伤累积直至断裂的过程。第三类流行于当前的设计准则,其通过非统一本构模型描述蠕变疲劳行为,并通过分别计算的蠕变损伤和疲劳损伤方程以及唯象的包络线来预测蠕变寿命。但是这三类方法有各自的缺点:第一类方法仅适用于描述微观层面的应力应变行为,对于宏观层面的大型高温部件并不适用;第二类方法侧重于描述裂纹扩展阶段的蠕变疲劳行为,其编程复杂、收敛性差、计算成本高等特点决定了这类方法不具有很强的普适性;第三类方法虽然具有操作性强的特点,但是其多用于分析稳态情况下单轴应力应变行为,对于复杂应力应变状态和复杂加载历史情况下的蠕变疲劳分析和预测并不准确。
基于此,期望获得一种新的蠕变疲劳预测方法,以更好地实现几何不连续结构在多轴应力应变状态下的蠕变疲劳分析,获得更为直观、精确度高的结果,且该蠕变疲劳预测方法实用性强。
发明内容
本发明的目的在于提供一种基于ABAQUS的多轴蠕变疲劳预测方法,该多轴蠕变疲劳预测方法可以更好地实现几何不连续结构在多轴应力应变状态下的蠕变疲劳分析,获得更为直观、精确度高的结果,且该蠕变疲劳预测方法实用性强。
根据上述发明目的,本发明提出了一种基于ABAQUS的多轴蠕变疲劳预测方法,其包括步骤:
S1:建立ABAQUS有限元模型,通过用户子程序UMAT定义待测试样材料在循环载荷过程中的粘塑性本构方程;
S2:确定所述粘塑性本构方程所需的模型参数;
S3:建立待测试材料的多轴应力应变状态的疲劳损伤计算模型和蠕变损伤计算模型;
S4:建立多轴应力应变状态下的ABAQUS有限元模型,基于步骤S1中用户子程序UMAT所定义的粘塑性本构方程和步骤S2中的模型参数,计算得到每个循环周次的应力应变张量;
S5:通过用户子程序USDFLD计算等效应力和等效塑性应变,基于步骤S3中的疲劳损伤计算模型和蠕变损伤计算模型并结合步骤S4中得到的应力应变张量,通过线性累计损伤准则叠加每个循环周次的疲劳损伤和蠕变损伤,得到待测试材料的裂纹萌生寿命。
在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中,步骤S4中的用户子程序UMAT可以计算每个节点的应力应变张量。步骤S5中的用户子程序USDFLD设置的作用为:(1)提取步骤S4的应力应变张量,并进行标量化运算,计算应力应变等效数值;(2)将步骤S3的疲劳损伤计算模型和蠕变损伤计算模型输入用户子程序USDFLD,以获得每个循环周次的疲劳损伤和蠕变损伤,上述损伤的计算可以基于等效应力应变的;(3)对每个循环周次的损伤进行线性累计损伤的叠加,以得到裂纹萌生寿命。
进一步地,在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中,在步骤S2中,在给定温度下对待测试材料进行单轴拉伸试验以及该给定温度下不同应变幅和保载时间的单轴蠕变疲劳试验,以得到高温拉伸曲线、循环软化曲线、应力松弛曲线和磁滞回线,用以确定所述粘塑性本构方程所需的模型参数。
进一步地,在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中,在步骤S2中,通过试参法模拟ABAQUS有限元模型的高温拉伸曲线、循环软化曲线、应力松弛曲线和磁滞回线,使之与实验结果吻合,也就是说,通过试参法所获得的曲线与实验所获得的曲线具有较好的拟合程度,以获取粘塑性本构方程所需的模型参数。
进一步地,在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中,在步骤S1中,粘塑性本构方程包括:粘塑性本构的主控方程、粘塑性本构的粘塑性方程、粘塑性本构的背应力张量非弹性随动强化方程和粘塑性本构的各项同性强化方程。
更进一步地,在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中,S1步骤进一步包括步骤:
S11:采用下述公式(1)和公式(2)描述粘塑性本构的主控方程:
εt=εe+εin (1);
其中,εt为总应变张量,εe为弹性应变张量,εin为非弹性应变张量,E为弹性模量,v为泊松比,σ为应力张量,trσ为应力张量的迹,I为二阶单位张量;
S12:采用下述公式(3)、公式(4)和公式(5)描述所述粘塑性本构的粘塑性方程:
其中,为非弹性应变率张量,为累积非弹性应变率,s为应力张量的偏量,a为背应力张量的偏量,J(σ-α)为Von-Mises应力空间距离,α为背应力张量,K和n为率相关的材料参数,R为各向同性变形抗性,κ为弹性区域初始尺寸,“:”表示张量的内积运算;
S13:采用下述公式(6)和公式(7)描述所述粘塑性本构的背应力张量非弹性随动强化方程:
m(q)=φ1e-q/ω+φ2 (7);
其中,αi表示若干个背应力张量部分中的每一个部分,ζi和ri为背应力张量每个部分的材料参数,γ为描述静力回复项的材料参数,m(q)为描述静力回复项的指数方程,q为塑性应变幅值,φ1,φ2和ω为指数方程中三个材料参数,J(αi)表示背应力第二不变量,e表示以自然常数为底的指数函数,表示背应力张量每个部分的应力变化率;
S14:采用下述公式(8)描述所述粘塑性本构的各项同性强化方程:
其中,Q为在第一阶段快速软化的各项同性抗力渐近值,b为接近渐近值的速度参数,H为第二阶段线性软化的斜率相关的参数,p为累积非弹性应变,表示各项同性强化率。
进一步地,在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中,其中背应力张量分为8个部分,即
进一步地,在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中,在所述S3步骤中,多轴应力应变状态的所述纯疲劳损伤计算模型为:
其中,“()max”表示临界平面上的最大疲劳损伤因子,τmax为临界平面上的最大剪切应力,τ′f为剪切疲劳强度常数,Δγ/2为临界平面上的剪切应变幅,σn,max为临界平面上的最大正应力,σ′f为疲劳强度常数,Δεn/2为临界平面上的正应变幅,G为剪切模量,df为一个循环周次的疲劳损伤,b0为疲劳强度指数,γ′f为剪切疲劳延性常数,c0为疲劳延性指数;
多轴应力应变状态的所述蠕变损伤计算模型为:
其中,dc为一个循环周次的蠕变损伤,th为一个循环周次的保载时间,Z为弹性跟随因子,t表示一个循环周次内从保载开始的时间,为蠕变损伤第一线性回归参数,MDF为多轴延性因子,n1为蠕变损伤第二线性回归参数,wf,trans为失效应变能密度平台值,为一个循环周次保载前最大等效应力,A为松弛第一参数,B为松弛第二参数,为等效弹性模量,为一个循环周次由疲劳产生的等效塑性应变范围,为一个循环周次的等效平均应力,n2为稳态蠕变指数,σH为静水应力,表示等效应力。
进一步地,在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中,在步骤S4中,建立多轴应力应变状态下的ABAQUS有限元模型,施加边界条件和外部载荷,划分模型网格,以得到各个积分点每个循环周次的应力应变张量。
进一步地,在本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法中,步骤S5进一步包括步骤:
S51:通过用户子程序USDFLD提取ABAQUS模型中每个节点的应力张量和应变张量;
S52:通过用户子程序USDFLD,并结合所述步骤S51,通过标量化运算获得每一个时刻的等效应力、等效塑性应变和每一个循环周次保载时间内的弹性跟随因子,最终得到每一个循环周次的疲劳损伤和蠕变损伤;
S53:通过用户子程序USDFLD采用下述公式(11)计算多轴蠕变疲劳条件下的总损伤:
其中,D(n)为前n个循环周次的累积总损伤,为第i个循环周次产生的疲劳损伤,为第i个循环周次产生的蠕变损伤;
其中,当某个节点的总损伤叠加率先到失效数值1时,即可定义为最危险节点并确定裂纹萌生寿命ni。
需要说明的是,在本发明所述的技术方案中,裂纹萌生寿命采用循环周次ni进行表征。
本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法具有如下所述的优点和有益效果:
(1)本发明所述的多轴蠕变疲劳预测方法利用基于用户子程序UMAT定义待测试材料的粘塑性本构方程,从而得到多轴应力应变状态下的蠕变疲劳行为;
(2)本发明所述的多轴蠕变疲劳预测方法利用基于用户子程序USDFLD计算等效应力和等效塑性应变,从而获得每个循环周次中各个积分点的蠕变损伤、疲劳损伤和总损伤数值;
(3)本发明所述的多轴蠕变疲劳预测方法具有很强的直观性,可以直观地获得几何不连续结构的裂纹萌生位置和该位置的裂纹萌生寿命。
附图说明
图1为本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的流程图。
图2示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴拉伸试验和模拟曲线的拟合结果图。
图3示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验循环软化数据和模拟曲线的拟合结果图。
图4示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验应力松弛数据和模拟曲线的拟合结果图。
图5示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验静力回复项的指数方程拟合结果图。
图6示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的两个潜在危险点每周次疲劳损伤和每周次蠕变损伤随循环周次的轨迹图。
图7示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的某一个潜在危险点前100个循环周次的磁滞回线。
图8示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的另一个潜在危险点前100个循环周次的磁滞回线。
图9示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的缺口根部次表面最危险积分点的蠕变疲劳损伤轨迹图和裂纹萌生寿命预测。
图10示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在另一种实施方式下的两个潜在危险点每周次疲劳损伤和每周次蠕变损伤随循环周次的轨迹图。
图11示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在另一种实施方式下的某一个潜在危险点前100个循环周次的磁滞回线。
图12示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在另一种实施方式下的另一个潜在危险点前100个循环周次的磁滞回线。
图13示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在另一种实施方式下的缺口根部表面的蠕变疲劳损伤轨迹图和裂纹萌生寿命预测。
具体实施方式
下面将结合说明书附图和具体的实施例对本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法做进一步的详细说明。
图1为本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的流程图。
如图1所示,在本实施方式中,基于ABAQUS的多轴蠕变疲劳预测方法包括步骤:
S1:建立ABAQUS有限元模型,通过用户子程序UMAT定义待测试材料在循环载荷过程中的粘塑性本构方程;
S2:确定所述粘塑性本构方程所需的模型参数;
S3:建立待测试材料的多轴应力应变状态的疲劳损伤计算模型和蠕变损伤计算模型;
S4:建立多轴应力应变状态下的ABAQUS有限元模型,基于步骤S1中用户子程序UMAT所定义的粘塑性本构方程和步骤S2中的模型参数,计算得到每个循环周次的应力应变张量;
S5:通过用户子程序USDFLD计算等效应力和等效塑性应变,基于步骤S3中的疲劳损伤计算模型和蠕变损伤计算模型并结合步骤S4中得到的应力应变张量,通过线性累计损伤准则叠加每个循环周次的疲劳损伤和蠕变损伤,得到待测试材料的裂纹萌生寿命。
其中,在步骤S2中,在给定温度下对待测试材料进行单轴拉伸试验以及该给定温度下不同应变幅和保载时间的单轴蠕变疲劳试验,以得到高温拉伸曲线、循环软化曲线、应力松弛曲线和磁滞回线,用以确定所述粘塑性本构方程所需的模型参数,并且可以通过试参法模拟ABAQUS有限元模型的高温拉伸曲线、循环软化曲线、应力松弛曲线和磁滞回线,以获取粘塑性本构方程所需的模型参数。
其中,在步骤S1中,粘塑性本构方程包括:粘塑性本构的主控方程、粘塑性本构的粘塑性方程、粘塑性本构的背应力张量非弹性随动强化方程和粘塑性本构的各项同性强化方程。其可以通过下述步骤获得:
S11:采用下述公式(1)和公式(2)描述所述粘塑性本构的主控方程:
εt=εe+εin (1);
其中,εt为总应变张量,εe为弹性应变张量,εin为非弹性应变张量,E为弹性模量,v为泊松比,σ为应力张量,trσ为应力张量的迹,I为二阶单位张量;
S12:采用下述公式(3)、公式(4)和公式(5)描述所述粘塑性本构的粘塑性方程:
其中,为非弹性应变率张量,为累积非弹性应变率,s为应力张量的偏量,a为背应力张量的偏量,J(σ-α)为Von-Mises应力空间距离,α为背应力张量,K和n为率相关的材料参数,R为各向同性变形抗性,κ为弹性区域初始尺寸,“:”表示张量的内积运算;
S13:采用下述公式(6)和公式(7)描述所述粘塑性本构的背应力张量非弹性随动强化方程:
m(q)=φ1e-q/ω+φ2 (7);
其中,αi表示若干个背应力张量部分中的每一个部分,ζi和ri为背应力张量每个部分的材料参数,γ为描述静力回复项的材料参数,m(q)为描述静力回复项的指数方程,q为塑性应变幅值,φ1,φ2和ω为指数方程中三个材料参数,J(αi)表示背应力第二不变量,e表示以自然常数为底的指数函数,表示背应力张量每个部分的应力变化率;
S14:采用下述公式(8)描述所述粘塑性本构的各项同性强化方程:
其中,Q为在第一阶段快速软化的各项同性抗力渐近值,b为接近渐近值的速度参数,H为第二阶段线性软化的斜率相关的参数,p为累积非弹性应变,表示各项同性强化率。
需要说明的是,上述步骤中,背应力张量分为8个部分,即
而在S3步骤中,多轴应力应变状态的所述纯疲劳损伤计算模型为:
其中,“( )max”表示临界平面上的最大疲劳损伤因子,τmax为临界平面上的最大剪切应力,τ′f为剪切疲劳强度常数,Δγ/2为临界平面上的剪切应变幅,σn,max为临界平面上的最大正应力,σ′f为疲劳强度常数,Δεn/2为临界平面上的正应变幅,G为剪切模量,df为一个循环周次的疲劳损伤,b0为疲劳强度指数,γ′f为剪切疲劳延性常数,c0为疲劳延性指数;
多轴应力应变状态的所述蠕变损伤计算模型为:
其中,dc为一个循环周次的蠕变损伤,th为一个循环周次的保载时间,Z为弹性跟随因子,t表示一个循环周次内从保载开始的时间,为蠕变损伤第一线性回归参数,MDF为多轴延性因子,n1为蠕变损伤第二线性回归参数,wf,trans为失效应变能密度平台值,为一个循环周次保载前最大等效应力,A为松弛第一参数,B为松弛第二参数,为等效弹性模量,为一个循环周次由疲劳产生的等效塑性应变范围,为一个循环周次的等效平均应力,n2为稳态蠕变指数,σH为静水应力,表示等效应力。
而在步骤S4中,建立多轴应力应变状态下的ABAQUS有限元模型,可以施加边界条件和外部载荷,划分模型网格,以得到各个积分点每个循环周次的应力应变张量
在本实施方式中,步骤S5可以进一步包括步骤
S51:通过用户子程序USDFLD提取ABAQUS模型中每个节点的应力张量和应变张量;
S52:通过用户子程序USDFLD,并结合所述步骤S51,通过标量化运算获得每一个时刻的等效应力、等效塑性应变和每一个循环周次保载时间内的弹性跟随因子,最终得到每一个循环周次的疲劳损伤和蠕变损伤;
S53:通过用户子程序USDFLD采用下述公式(11)计算多轴蠕变疲劳条件下的总损伤:
其中,D(n)为前n个循环周次的累积总损伤,为第i个循环周次产生的疲劳损伤,为第i个循环周次产生的蠕变损伤;
其中,当某个节点的总损伤叠加率先到失效数值1时,即可定义为最危险节点并确定裂纹萌生寿命ni。
为了更好地阐述本案的基于AVAQUS的多轴蠕变疲劳预测方法的预测效果,将采用缺口半径为8毫米的单边缺口试样进行验证。
验证所采用的单边缺口试样的材料为高温镍基GH4169超合金,蠕变疲劳试验在650℃空气环境下进行。试验时,试样两端施加的外载和为整体的应变控制,由于单边缺口试样几何不连续性的影响,其缺口根本最薄弱的部位处于多轴应力应变状态。在此之前,需要进行650℃空气环境下的单轴拉伸试验以及该环境下不同应变幅和保载时间的单轴蠕变疲劳试验,获得的试验结果用于确定公式(1)~(8)粘塑性本构方程所需的材料参数。
通过试参法调整单轴拉伸试验的模拟结果,使之与单轴拉伸试验能够较好地吻合。图2示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴拉伸试验和模拟曲线的拟合结果图。
如图2所示,I表示单轴拉伸试验数据,II表示采用试参法模拟得到的单轴拉伸试验曲线,通过试参法调整单轴拉伸试验曲线II的模拟结果,使之与单轴蠕变疲劳试验所得的单轴拉伸试验数据能够较好地吻合,其中,弹性模量E=177GPa,泊松比v=0.33,弹性区域初始尺寸κ=815MPa,各个背应力张量每个部分的随动强化材料参数:ζ1=6130,ζ2=1807,ζ3=892,ζ4=352,ζ5=150,ζ6=88.2,ζ7=75.0,ζ8=28.4,r1=23.4,r2=68.0,r3=75.9,r4=48.0,r5=43.4,r6=25.4,r7=54.5,r8=28.0,率相关的粘塑性材料参数:K=400,n=2.0。
图3示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验循环软化数据和模拟曲线的拟合结果图。
如图3所示,III表示循环软化试验数据,IV表示采用试参法模拟得到的循环软化曲线,通过试参法调整循环软化曲线IV的模拟结果,使之与单轴蠕变疲劳试验所得的循环软化数据能够较好地吻合,其中,各项同性强化的材料参数:第二阶段线性软化的斜率相关的参数H=-8.5,接近渐近值的速度参数b=4.1,第一阶段快速软化的各项同性抗力渐近值Q=618MPa。
图4示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验应力松弛数据和模拟曲线的拟合结果图。
如图4所示,V表示单轴蠕变疲劳试验应力松弛数据,VI表示采用试参法模拟得到的应力松弛曲线,通过试参法调整第一个循环周次的应力松弛曲线的模拟结果,使之与单轴蠕变疲劳试验所得的应力松弛数据能够较好地吻合,其中,静力回复项中独立材料参数:γ=4.0×10-7。
与此同时,拟合公式(7)中的静力回复项的指数方程,拟合结果图如图5所示。图5示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的单轴蠕变疲劳试验静力回复项的指数方程拟合结果图。
如图5所示,VII表示单轴蠕变疲劳试验数据,VIII表示模拟结果曲线,其中,静力回复项中指数方程的材料参数:φ1=0.37,φ2=2.82,ω=6.6×10-4。
在本案的一个实施方式中,选取数据为总应变范围为1.0%,每个循环周次应变最大处的保载时间为3600s,650℃空气环境下的单边缺口蠕变疲劳试验。其蠕变疲劳裂纹萌生于缺口根部次表面的位置,裂纹萌生寿命为76个循环周次。
图6给出了该实施方式下两个典型位置每周次疲劳损伤和蠕变损伤的轨迹图。
如图6所示,曲线IX表示缺口根部表面每周次疲劳损伤曲线,曲线X表示缺口根部次表面每周次疲劳损伤曲线,曲线XI表示缺口根部表面每周次蠕变损伤曲线,曲线XII表示缺口根部次表面每周次蠕变损伤曲线,其中,选取缺口根部表面的积分点,因其在循环载荷的过程中受到应力集中的影响,疲劳裂纹通常萌生于表面,所以缺口根部表面的积分点通常被认为是疲劳过程中的潜在危险点;而选取缺口根部次表面的积分点,是由于次表面通常具有相比于表面更高的应力三轴度,蠕变裂纹通常萌生于内部,所以缺口根部次表面的积分点通常被认为是蠕变过程中的潜在危险点。而从图6中可以看出:此实施方式下由于保载时间较长,循环载荷的过程中蠕变损伤占主导地位,所计算最大总损伤的位置在缺口根部次表面,与实验所观测的现象完全一致。
图7和图8给出了图6所示的缺口根部表面和次表面两个位置前100个循环周次的磁滞回线。其中,图7表示缺口根部表面积分点处的磁滞回线,图8表示缺口根部次表面积分点处的磁滞回线,图7和图8中的A1表示第一个循环周次,A2表示第一百个循环周次。
结合图7和图8中可以看出:虽然缺口根部表面积分点具有更大的非弹性应变范围,但是次表面积分点具有更加明显的非弹性应变的累积。这是由于受到多轴应力应变状态的影响,几何不连续处在保载阶段处于应力应变混合控制的模式。结合图7和图8还可以看出,次表面积分点更加接近应力控制模式,导致其蠕变损伤更加明显,进而解释了蠕变主导的循环载荷其裂纹萌生位置出现在缺口根部次表面。
图9示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在一种实施方式下的缺口根部次表面最危险积分点的蠕变疲劳损伤轨迹图和裂纹萌生寿命预测。
借助蠕变疲劳损伤交互图,追踪缺口根部附近最危险的次表面积分点,如图9所示,可以获得该实施方式下的裂纹萌生寿命为105个循环周次,这与实验测得的76个循环周次较为接近,在1.5倍误差带范围之内证明此数值模拟方法具有较高的可靠性。
需要说明的是,图9中的曲线XIII表示
在本案的另一个实施方式中,选取数据为总应变范围为1.0%,每个循环周次应变最大处的保载时间为60s,650℃空气环境下的单边缺口蠕变疲劳试验。其蠕变疲劳裂纹萌生于缺口根部表面的位置,裂纹萌生寿命为480个循环周次。
图10给出了该实施方式下两个典型位置每周次疲劳损伤和蠕变损伤的轨迹图。
如图10所示,曲线XIV表示缺口根部表面每周次疲劳损伤曲线,曲线XV表示缺口根部次表面每周次疲劳损伤曲线,曲线XVI表示缺口根部表面每周次蠕变损伤曲线,曲线XVII表示缺口根部次表面每周次蠕变损伤曲线,其中,选取缺口根部表面的积分点,是因其在循环载荷的过程中受到应力集中的影响,疲劳裂纹通常萌生于表面,所以缺口根部表面的积分点通常被认为是疲劳过程中的潜在危险点;选取缺口根部次表面的积分点,是由于次表面通常具有相比于表面更高的应力三轴度,蠕变裂纹通常萌生于内部,所以缺口根部次表面的积分点通常被认为是蠕变过程中的潜在危险点。而从图10中可以看出:在此实施方式下由于保载时间较短,循环载荷的过程中疲劳损伤占主导地位,所计算最大总损伤的位置在缺口根部表面,与实验所观测的现象完全一致。
图11和图12给出了图10所示的缺口根部表面和次表面两个位置前100个循环周次的磁滞回线。其中,图11表示缺口根部表面积分点处的磁滞回线,图12表示缺口根部次表面积分点处的磁滞回线,图11和图12中的A1表示第一个循环周次,A2表示第一百个循环周次。
结合图11和图12可以看出,缺口根部表面积分点具有更大的非弹性应变范围,而且缺口根本表面和次表面两个位置的积分点具有几乎没有非弹性应变的累积。结合图11和图12可以看出,该实施方式下保载时间较短导致的蠕变/松弛的现象几乎可以忽略不计,在这种情况下,由非弹性应变范围导致的疲劳损伤使得裂纹萌生位置出现在缺口根部表面。
图13示意性地显示了本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法在另一种实施方式下的缺口根部表面的蠕变疲劳损伤轨迹图和裂纹萌生寿命预测。
借助蠕变疲劳损伤交互图,追踪缺口根部表面最危险的积分点,如图13所示,可以获得该实施方式下的裂纹萌生寿命为346个循环周次,这与实验测得的480个循环周次较为接近,在1.5倍误差带范围之内证明此数值模拟方法具有较高的可靠性。
需要说明的是,图13中的曲线XVIII表示
综上所述可以看出,本发明所述的基于ABAQUS的多轴蠕变疲劳预测方法利用基于用户子程序UMAT定义待测试材料的粘塑性本构方程,从而得到多轴应力应变状态下的蠕变疲劳行为。
此外,本发明所述的多轴蠕变疲劳预测方法利用基于用户子程序USDFLD计算等效应力和等效塑性应变,从而获得每个循环周次中各个积分点的蠕变损伤、疲劳损伤和总损伤数值;
另外,本发明所述的多轴蠕变疲劳预测方法具有很强的直观性,可以直观地获得几何不连续结构的裂纹萌生位置和该位置的裂纹萌生寿命。
需要说明的是,本发明的保护范围中现有技术部分并不局限于本申请文件所给出的实施例,所有不与本发明的方案相矛盾的现有技术,包括但不局限于在先专利文献、在先公开出版物,在先公开使用等等,都可纳入本发明的保护范围。
此外,本案中各技术特征的组合方式并不限本案权利要求中所记载的组合方式或是具体实施例所记载的组合方式,本案记载的所有技术特征可以以任何方式进行自由组合或结合,除非相互之间产生矛盾。
还需要注意的是,以上所列举的实施例仅为本发明的具体实施例。显然本发明不局限于以上实施例,随之做出的类似变化或变形是本领域技术人员能从本发明公开的内容直接得出或者很容易便联想到的,均应属于本发明的保护范围。
Claims (9)
1.一种基于ABAQUS的多轴蠕变疲劳预测方法,其特征在于,包括步骤:
S1:建立ABAQUS有限元模型,通过用户子程序UMAT定义待测试材料在循环载荷过程中的粘塑性本构方程;
S2:确定所述粘塑性本构方程所需的模型参数;
S3:建立待测试材料的多轴应力应变状态的疲劳损伤计算模型和蠕变损伤计算模型;
S4:建立多轴应力应变状态下的ABAQUS有限元模型,基于步骤S1中用户子程序UMAT所定义的粘塑性本构方程和步骤S2中的模型参数,计算得到每个循环周次的应力应变张量;
S5:通过用户子程序USDFLD计算等效应力和等效塑性应变,基于步骤S3中的疲劳损伤计算模型和蠕变损伤计算模型并结合步骤S4中得到的应力应变张量,通过线性累计损伤准则叠加每个循环周次的疲劳损伤和蠕变损伤,得到待测试材料的裂纹萌生寿命。
2.如权利要求1所述的基于ABAQUS的多轴蠕变疲劳预测方法,其特征在于,在所述步骤S2中,在给定温度下对待测试材料进行单轴拉伸试验以及该给定温度下不同应变幅和保载时间的单轴蠕变疲劳试验,以得到高温拉伸曲线、循环软化曲线、应力松弛曲线和磁滞回线,用以确定所述粘塑性本构方程所需的模型参数。
3.如权利要求2所述的基于ABAQUS的多轴蠕变疲劳预测方法,其特征在于,在所述步骤S2中,通过试参法模拟ABAQUS有限元模型的高温拉伸曲线、循环软化曲线、应力松弛曲线和磁滞回线。
4.如权利要求1所述的基于ABAQUS的多轴蠕变疲劳预测方法,其特征在于,在步骤S1中,所述粘塑性本构方程包括:粘塑性本构的主控方程、粘塑性本构的粘塑性方程、粘塑性本构的背应力张量非弹性随动强化方程和粘塑性本构的各项同性强化方程。
5.如权利要求4所述的基于ABAQUS的多轴蠕变疲劳预测方法,其特征在于,所述S1步骤进一步包括步骤:
S11:采用下述公式(1)和公式(2)描述所述粘塑性本构的主控方程:
εt=εe+εin (1);
其中,εt为总应变张量,εe为弹性应变张量,εin为非弹性应变张量,E为弹性模量,ν为泊松比,σ为应力张量,trσ为应力张量的迹,I为二阶单位张量;
S12:采用下述公式(3)、公式(4)和公式(5)描述所述粘塑性本构的粘塑性方程:
其中,为非弹性应变率张量,为累积非弹性应变率,s为应力张量的偏量,a为背应力张量的偏量,J(σ-α)为Von-Mises应力空间距离,α为背应力张量,K和n为率相关的材料参数,R为各向同性变形抗性,K为弹性区域初始尺寸,“:”表示张量的内积运算;
S13:采用下述公式(6)和公式(7)描述所述粘塑性本构的背应力张量非弹性随动强化方程:
m(q)=φ1e-q/ω+φ2 (7);
其中,αi表示若干个背应力张量部分中的每一个部分,ζi和ri为背应力张量每个部分的材料参数,γ为描述静力回复项的材料参数,m(q)为描述静力回复项的指数方程,q为塑性应变幅值,φ1,φ2和ω为指数方程中三个材料参数,J(αi)表示背应力第二不变量,e表示以自然常数为底的指数函数,表示背应力张量每个部分的应力变化率;
S14:采用下述公式(8)描述所述粘塑性本构的各项同性强化方程:
其中,Q为在第一阶段快速软化的各项同性抗力渐近值,b为接近渐近值的速度参数,H为第二阶段线性软化的斜率相关的参数,p为累积非弹性应变,表示各项同性强化率。
6.如权利要求5所述的基于ABAQUS的多轴蠕变疲劳预测方法,其特征在于,其中背应力张量分为8个部分,即
7.如权利要求1所述的基于ABAQUS的多轴蠕变疲劳预测方法,其特征在于,在所述S3步骤中,多轴应力应变状态的所述纯疲劳损伤计算模型为:
其中,“()max”表示临界平面上的最大疲劳损伤因子,τmmax为临界平面上的最大剪切应力,τ′f为剪切疲劳强度常数,Δγ/2为临界平面上的剪切应变幅,σn,max为临界平面上的最大正应力,σ′f为疲劳强度常数,Δεn/2为临界平面上的正应变幅,G为剪切模量,df为一个循环周次的疲劳损伤,b0为疲劳强度指数,γ′f为剪切疲劳延性常数,c0为疲劳延性指数;
多轴应力应变状态的所述蠕变损伤计算模型为:
其中,dc为一个循环周次的蠕变损伤,th为一个循环周次的保载时间,Z为弹性跟随因子,t表示一个循环周次内从保载开始的时间,为蠕变损伤第一线性回归参数,MDF为多轴延性因子,n1为蠕变损伤第二线性回归参数,wf,trans为失效应变能密度平台值,为一个循环周次保载前最大等效应力,A为松弛第一参数,B为松弛第二参数,E为等效弹性模量,为一个循环周次由疲劳产生的等效塑性应变范围,为一个循环周次的等效平均应力,n2为稳态蠕变指数,σH为静水应力,表示等效应力。
8.如权利要求1所述的基于ABAQUS的多轴蠕变疲劳预测方法,其特征在于,在所述步骤S4中,建立多轴应力应变状态下的ABAQUS有限元模型,施加边界条件和外部载荷,划分模型网格,以得到各个积分点每个循环周次的应力应变张量。
9.如权利要求1所述的基于ABAQUS的多轴蠕变疲劳预测方法,其特征在于,
所述步骤S5进一步包括步骤:
S51:通过用户子程序USDFLD提取ABAQUS模型中每个节点的应力张量和应变张量;
S52:通过用户子程序USDFLD,并结合所述步骤S51,通过标量化运算获得每一个时刻的等效应力、等效塑性应变和每一个循环周次保载时间内的弹性跟随因子,最终得到每一个循环周次的疲劳损伤和蠕变损伤;
S53:通过用户子程序USDFLD采用下述公式(11)计算多轴蠕变疲劳条件下的总损伤:
其中,D(n)为前n个循环周次的累积总损伤,为第i个循环周次产生的疲劳损伤,为第i个循环周次产生的蠕变损伤;
其中,当某个节点的总损伤叠加率先到失效数值1时,即可定义为最危险节点,并确定裂纹萌生寿命ni。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910026871.3A CN109885874B (zh) | 2019-01-11 | 2019-01-11 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
PCT/CN2019/114718 WO2020143284A1 (zh) | 2019-01-11 | 2019-10-31 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
US17/421,929 US20230103274A9 (en) | 2019-01-11 | 2019-10-31 | A Multiaxial Creep-Fatigue Prediction Method Based On ABAQUS |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910026871.3A CN109885874B (zh) | 2019-01-11 | 2019-01-11 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109885874A true CN109885874A (zh) | 2019-06-14 |
CN109885874B CN109885874B (zh) | 2022-12-23 |
Family
ID=66925793
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910026871.3A Active CN109885874B (zh) | 2019-01-11 | 2019-01-11 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20230103274A9 (zh) |
CN (1) | CN109885874B (zh) |
WO (1) | WO2020143284A1 (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110411851A (zh) * | 2019-07-22 | 2019-11-05 | 北京科技大学 | 一种高温合金涡轮叶片服役损伤评价及蠕变寿命预测方法 |
CN110826285A (zh) * | 2019-11-22 | 2020-02-21 | 华东理工大学 | 一种几何不连续结构的疲劳寿命预测方法 |
CN110967267A (zh) * | 2019-11-25 | 2020-04-07 | 中国民用航空飞行学院 | 一种判定疲劳裂纹萌生寿命的试验方法 |
WO2020143284A1 (zh) * | 2019-01-11 | 2020-07-16 | 华东理工大学 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
CN111460583A (zh) * | 2020-04-14 | 2020-07-28 | 华东理工大学 | 一种针对复杂几何结构件的蠕变-疲劳寿命设计方法 |
CN111504779A (zh) * | 2020-04-07 | 2020-08-07 | 武汉大学 | 利用脆性指数确定岩石软化曲线的方法及装置 |
CN112364535A (zh) * | 2020-11-06 | 2021-02-12 | 华东理工大学 | 一种基于晶体塑性的蠕变疲劳寿命预测方法 |
CN112711901A (zh) * | 2020-12-10 | 2021-04-27 | 华南理工大学 | 基于泛协同Kriging模型的机构疲劳寿命预测方法 |
CN112836342A (zh) * | 2021-01-05 | 2021-05-25 | 北京航空航天大学 | 一种基于连续损伤力学的微米烧结银芯片粘接层疲劳失效物理模型建模与验证方法 |
CN112948941A (zh) * | 2021-03-16 | 2021-06-11 | 宁波大学 | 受剪构件的高周疲劳损伤数值计算方法、装置 |
CN113049376A (zh) * | 2021-03-10 | 2021-06-29 | 华东理工大学 | 一种用于过热器管板的蠕变疲劳损伤评估方法 |
CN113514343A (zh) * | 2021-07-13 | 2021-10-19 | 华东理工大学 | 一种多轴蠕变性能参数的测试方法 |
CN113611377A (zh) * | 2021-07-05 | 2021-11-05 | 南京工业大学 | 一种利用晶体塑性模型模拟混合控制蠕变疲劳变形的方法 |
CN113627054A (zh) * | 2021-08-02 | 2021-11-09 | 东南大学 | 一种优化的延性金属损伤参数确定方法 |
CN113722946A (zh) * | 2021-07-30 | 2021-11-30 | 淮浙电力有限责任公司凤台发电分公司 | 一种汽轮机转子的蠕变-疲劳寿命预测方法及预测系统 |
CN113742914A (zh) * | 2021-09-02 | 2021-12-03 | 南京工业大学 | 适用于预测多种控制模式循环载荷变形行为的方法 |
CN114462147A (zh) * | 2022-01-28 | 2022-05-10 | 中国人民解放军陆军工程大学 | 含损伤的推进剂蠕变型本构模型的构建与有限元应用方法 |
CN114493012A (zh) * | 2022-01-28 | 2022-05-13 | 武汉科技大学 | 考虑硬化损伤协同效应的胶结充填体蠕变失稳预测方法及装置 |
CN114622907A (zh) * | 2022-05-12 | 2022-06-14 | 中国科学院地质与地球物理研究所 | 油气储层最小水平主应力剖面模型预测方法、系统和设备 |
CN115964919A (zh) * | 2023-02-06 | 2023-04-14 | 南京工业大学 | 熔盐腐蚀-力交互作用下晶间裂纹萌生及扩展的预测方法 |
CN115982958A (zh) * | 2022-12-07 | 2023-04-18 | 南京工业大学 | 一种基于工程损伤力学的材料蠕变疲劳寿命预测方法 |
CN116052813A (zh) * | 2022-12-30 | 2023-05-02 | 南京工业大学 | 一种全数值驱动的实时蠕变疲劳损伤监测方法 |
Families Citing this family (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111950101A (zh) * | 2020-08-07 | 2020-11-17 | 中铁轨道交通装备有限公司 | 一种基于张量的多轴疲劳分析方法 |
CN112163352B (zh) * | 2020-08-13 | 2024-04-19 | 西安建筑科技大学 | 中碳钢超细晶棒材3d-spd轧制损伤预测方法及模型 |
USD997992S1 (en) | 2020-08-21 | 2023-09-05 | Vulcan Industrial Holdings, LLC | Fluid end for a pumping system |
CN112115540B (zh) * | 2020-09-27 | 2022-12-20 | 武汉第二船舶设计研究所(中国船舶重工集团公司第七一九研究所) | 一种融合实测数据的海工支撑结构疲劳分析方法及系统 |
CN112520064B (zh) * | 2020-12-04 | 2022-09-20 | 中国航空工业集团公司成都飞机设计研究所 | 一种基于应变监测的损伤自动识别方法 |
CN112711835B (zh) * | 2020-12-16 | 2022-10-21 | 北京理工大学 | 一种基于修正塑性应变能的金属材料疲劳寿命预测方法 |
CN112668167A (zh) * | 2020-12-21 | 2021-04-16 | 合图智造科技(西安)有限公司 | 基于少量实验数据的材料参数构建方法 |
CN113297754B (zh) * | 2021-06-30 | 2024-02-06 | 华中科技大学 | 一种太阳能高温吸热器疲劳-蠕变损伤评估方法及系统 |
CN113378327B (zh) * | 2021-07-02 | 2022-09-02 | 吉林重通成飞新材料股份公司 | 风电叶片后缘抗开裂设计方法、装置、设备及存储介质 |
CN113642207A (zh) * | 2021-07-09 | 2021-11-12 | 恒大新能源技术(深圳)有限公司 | 金属失效模型构建方法、装置、终端设备及存储介质 |
US11846356B1 (en) | 2021-08-18 | 2023-12-19 | Vulcan Industrial Holdings, LLC | Self-locking plug |
CN113884371B (zh) * | 2021-09-30 | 2024-04-02 | 西南石油大学 | 一种基于温度效应计算围岩蠕变时间的方法及其计算装置 |
CN114112676B (zh) * | 2021-12-03 | 2023-12-19 | 中国林业科学研究院木材工业研究所 | 一种木材横纹抗压全时程本构关系的构建方法 |
CN114459941B (zh) * | 2021-12-29 | 2023-11-21 | 宁波职业技术学院 | 一种排气系统中静子件开裂风险的预测方法与系统 |
CN114065558B (zh) * | 2022-01-18 | 2022-04-15 | 中国科学院力学研究所 | 一种浮动摩擦片全流程损伤识别、寿命预测方法及设备 |
CN114462146B (zh) * | 2022-01-28 | 2023-02-03 | 中国人民解放军陆军工程大学 | 考虑老化损伤的推进剂蠕变本构模型的构建与有限元应用方法 |
CN114580098A (zh) | 2022-02-09 | 2022-06-03 | 南京航空航天大学 | 基于虚拟应变能的金属材料多轴疲劳寿命预测方法与系统 |
CN114547753B (zh) * | 2022-03-03 | 2024-03-22 | 中交基础设施养护集团有限公司 | 考虑地层空间变异性的长条形地下结构分段长度确定方法 |
CN114611389A (zh) * | 2022-03-04 | 2022-06-10 | 北京航空航天大学 | 一种基于人工智能的复合材料失效高效模拟方法 |
CN114741922B (zh) * | 2022-04-11 | 2024-02-23 | 西安交通大学 | 一种基于Attention机制的透平叶片蠕变-疲劳寿命预测方法 |
CN114943121B (zh) * | 2022-04-12 | 2024-02-23 | 西北工业大学 | 一种引入残余应力的摩擦组焊增材件蠕变孕育期预测方法 |
US11434900B1 (en) | 2022-04-25 | 2022-09-06 | Vulcan Industrial Holdings, LLC | Spring controlling valve |
US11920684B1 (en) | 2022-05-17 | 2024-03-05 | Vulcan Industrial Holdings, LLC | Mechanically or hybrid mounted valve seat |
CN115183965B (zh) * | 2022-05-17 | 2023-08-08 | 中铁西北科学研究院有限公司 | 一种适用于振动台试验的隧道衬砌地震累积损伤评价方法 |
CN115204003B (zh) * | 2022-06-15 | 2023-09-29 | 中国人民解放军火箭军工程大学 | 一种聚合物材料非线性蠕变型粘弹性本构模型的构建方法 |
CN115310311B (zh) * | 2022-07-01 | 2023-08-01 | 南京航空航天大学 | 一种金属材料在谱载荷下疲劳寿命计算方法 |
CN115640666B (zh) * | 2022-07-25 | 2023-03-28 | 南京航空航天大学 | 一种基于损伤等效的航空发动机加速任务试车谱编制方法 |
CN115358113B (zh) * | 2022-07-29 | 2023-04-07 | 中国航发沈阳发动机研究所 | 一种航空发动机用聚四氟乙烯软管脉冲寿命计算方法 |
CN115017745B (zh) * | 2022-08-09 | 2022-10-21 | 天津大学 | 一种用于火电机组结构件的蠕变疲劳寿命预测方法及系统 |
CN115640683B (zh) * | 2022-10-20 | 2024-02-23 | 青岛理工大学 | 一种岩土体蠕变预测方法及系统 |
CN115392807B (zh) * | 2022-10-31 | 2023-04-18 | 南京复创大数据产业发展有限公司 | 一种基于过程工艺数据的大数据质量分析方法和系统 |
CN115906583B (zh) * | 2022-12-16 | 2023-08-01 | 中国人民解放军陆军工程大学 | 一种基于虚单元法的药柱结构完整性仿真分析方法及系统 |
CN116046980B (zh) * | 2022-12-16 | 2023-09-01 | 北京航空航天大学 | 一种基于应变监测的结构疲劳损伤诊断方法 |
CN116029180B (zh) * | 2023-03-30 | 2023-07-14 | 湖南云箭科技有限公司 | 机载外挂物疲劳仿真方法及系统 |
CN116296820B (zh) * | 2023-03-31 | 2024-04-05 | 四川大学 | 双向逼近确定岩石弹性参数与起裂应力的迭代方法及产品 |
CN116842778B (zh) * | 2023-05-29 | 2024-01-26 | 天津大学 | 基于紧凑拉伸-剪切试样的疲劳裂纹扩展方向、长度及速率计算方法 |
CN116818292B (zh) * | 2023-06-01 | 2024-02-23 | 南京航空航天大学 | 一种航空发动机安全起落次数确定方法 |
CN116895350B (zh) * | 2023-08-04 | 2024-01-16 | 辽宁工业大学 | 一种在复合位移加载下波纹管的多轴疲劳寿命预测方法 |
CN117350138B (zh) * | 2023-12-06 | 2024-02-23 | 西北工业大学 | 一种芯片粘接结构热循环疲劳失效物理模型建模方法 |
CN117558381B (zh) * | 2024-01-12 | 2024-03-22 | 四川大学 | 一种金属材料温度与应变速率相关塑性硬化模型计算方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101853317A (zh) * | 2010-04-20 | 2010-10-06 | 北京航空航天大学 | 一种涡轮盘结构概率设计系统的构建方法 |
JP2013148540A (ja) * | 2012-01-23 | 2013-08-01 | Babcock Hitachi Kk | 高Cr鋼の応力集中部に適したクリープ損傷評価方法 |
CN103926152A (zh) * | 2014-04-09 | 2014-07-16 | 北京工业大学 | 一种高温多轴谱载下低周蠕变-疲劳寿命评估方法 |
US20140336954A1 (en) * | 2013-05-13 | 2014-11-13 | Michael Bruyneel | Calculating Fatigue and Fatigue Failure of Structures |
CN105158084A (zh) * | 2015-09-15 | 2015-12-16 | 华东理工大学 | 一种材料的蠕变-疲劳寿命预测方法 |
CN108931448A (zh) * | 2018-05-07 | 2018-12-04 | 华南理工大学 | 一种高铬钢材料热力学响应及疲劳-蠕变损伤的预测方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150168262A1 (en) * | 2013-10-11 | 2015-06-18 | Solar Turbines Incorporated | Single crystal turbine blade lifing process and system |
CN106202913B (zh) * | 2016-07-07 | 2018-03-06 | 华东理工大学 | 时间相关的蠕变疲劳损伤评定方法 |
CN106934168B (zh) * | 2017-03-21 | 2018-08-28 | 中国石油大学(华东) | 一种材料多轴蠕变失效应变预测方法 |
CN109885874B (zh) * | 2019-01-11 | 2022-12-23 | 华东理工大学 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
-
2019
- 2019-01-11 CN CN201910026871.3A patent/CN109885874B/zh active Active
- 2019-10-31 WO PCT/CN2019/114718 patent/WO2020143284A1/zh active Application Filing
- 2019-10-31 US US17/421,929 patent/US20230103274A9/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101853317A (zh) * | 2010-04-20 | 2010-10-06 | 北京航空航天大学 | 一种涡轮盘结构概率设计系统的构建方法 |
JP2013148540A (ja) * | 2012-01-23 | 2013-08-01 | Babcock Hitachi Kk | 高Cr鋼の応力集中部に適したクリープ損傷評価方法 |
US20140336954A1 (en) * | 2013-05-13 | 2014-11-13 | Michael Bruyneel | Calculating Fatigue and Fatigue Failure of Structures |
CN103926152A (zh) * | 2014-04-09 | 2014-07-16 | 北京工业大学 | 一种高温多轴谱载下低周蠕变-疲劳寿命评估方法 |
CN105158084A (zh) * | 2015-09-15 | 2015-12-16 | 华东理工大学 | 一种材料的蠕变-疲劳寿命预测方法 |
CN108931448A (zh) * | 2018-05-07 | 2018-12-04 | 华南理工大学 | 一种高铬钢材料热力学响应及疲劳-蠕变损伤的预测方法 |
Cited By (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020143284A1 (zh) * | 2019-01-11 | 2020-07-16 | 华东理工大学 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
CN110411851A (zh) * | 2019-07-22 | 2019-11-05 | 北京科技大学 | 一种高温合金涡轮叶片服役损伤评价及蠕变寿命预测方法 |
WO2021098169A1 (zh) * | 2019-11-22 | 2021-05-27 | 华东理工大学 | 一种几何不连续结构的疲劳寿命预测方法 |
CN110826285A (zh) * | 2019-11-22 | 2020-02-21 | 华东理工大学 | 一种几何不连续结构的疲劳寿命预测方法 |
CN110967267A (zh) * | 2019-11-25 | 2020-04-07 | 中国民用航空飞行学院 | 一种判定疲劳裂纹萌生寿命的试验方法 |
CN111504779A (zh) * | 2020-04-07 | 2020-08-07 | 武汉大学 | 利用脆性指数确定岩石软化曲线的方法及装置 |
CN111460583A (zh) * | 2020-04-14 | 2020-07-28 | 华东理工大学 | 一种针对复杂几何结构件的蠕变-疲劳寿命设计方法 |
CN112364535A (zh) * | 2020-11-06 | 2021-02-12 | 华东理工大学 | 一种基于晶体塑性的蠕变疲劳寿命预测方法 |
CN112364535B (zh) * | 2020-11-06 | 2024-03-29 | 华东理工大学 | 一种基于晶体塑性的蠕变疲劳寿命预测方法 |
CN112711901A (zh) * | 2020-12-10 | 2021-04-27 | 华南理工大学 | 基于泛协同Kriging模型的机构疲劳寿命预测方法 |
CN112711901B (zh) * | 2020-12-10 | 2023-08-22 | 华南理工大学 | 基于泛协同Kriging模型的机构疲劳寿命预测方法 |
CN112836342A (zh) * | 2021-01-05 | 2021-05-25 | 北京航空航天大学 | 一种基于连续损伤力学的微米烧结银芯片粘接层疲劳失效物理模型建模与验证方法 |
CN113049376A (zh) * | 2021-03-10 | 2021-06-29 | 华东理工大学 | 一种用于过热器管板的蠕变疲劳损伤评估方法 |
CN113049376B (zh) * | 2021-03-10 | 2022-05-17 | 华东理工大学 | 一种用于过热器管板的蠕变疲劳损伤评估方法 |
CN112948941A (zh) * | 2021-03-16 | 2021-06-11 | 宁波大学 | 受剪构件的高周疲劳损伤数值计算方法、装置 |
CN112948941B (zh) * | 2021-03-16 | 2023-12-12 | 宁波大学 | 受剪构件的高周疲劳损伤数值计算方法、装置 |
CN113611377A (zh) * | 2021-07-05 | 2021-11-05 | 南京工业大学 | 一种利用晶体塑性模型模拟混合控制蠕变疲劳变形的方法 |
CN113611377B (zh) * | 2021-07-05 | 2023-08-22 | 南京工业大学 | 一种利用晶体塑性模型模拟混合控制蠕变疲劳变形的方法 |
CN113514343A (zh) * | 2021-07-13 | 2021-10-19 | 华东理工大学 | 一种多轴蠕变性能参数的测试方法 |
CN113514343B (zh) * | 2021-07-13 | 2022-05-06 | 华东理工大学 | 一种多轴蠕变性能参数的测试方法 |
CN113722946A (zh) * | 2021-07-30 | 2021-11-30 | 淮浙电力有限责任公司凤台发电分公司 | 一种汽轮机转子的蠕变-疲劳寿命预测方法及预测系统 |
CN113627054B (zh) * | 2021-08-02 | 2022-11-22 | 东南大学 | 一种优化的延性金属损伤参数确定方法 |
CN113627054A (zh) * | 2021-08-02 | 2021-11-09 | 东南大学 | 一种优化的延性金属损伤参数确定方法 |
CN113742914B (zh) * | 2021-09-02 | 2023-08-08 | 南京工业大学 | 适用于预测多种控制模式循环载荷变形行为的方法 |
CN113742914A (zh) * | 2021-09-02 | 2021-12-03 | 南京工业大学 | 适用于预测多种控制模式循环载荷变形行为的方法 |
CN114462147A (zh) * | 2022-01-28 | 2022-05-10 | 中国人民解放军陆军工程大学 | 含损伤的推进剂蠕变型本构模型的构建与有限元应用方法 |
CN114493012B (zh) * | 2022-01-28 | 2023-08-08 | 武汉科技大学 | 考虑硬化损伤协同效应的胶结充填体蠕变失稳预测方法及装置 |
CN114462147B (zh) * | 2022-01-28 | 2023-02-03 | 中国人民解放军陆军工程大学 | 含损伤的推进剂蠕变型本构模型的构建与有限元应用方法 |
CN114493012A (zh) * | 2022-01-28 | 2022-05-13 | 武汉科技大学 | 考虑硬化损伤协同效应的胶结充填体蠕变失稳预测方法及装置 |
CN114622907B (zh) * | 2022-05-12 | 2022-07-19 | 中国科学院地质与地球物理研究所 | 油气储层最小水平主应力剖面模型预测方法、系统和设备 |
CN114622907A (zh) * | 2022-05-12 | 2022-06-14 | 中国科学院地质与地球物理研究所 | 油气储层最小水平主应力剖面模型预测方法、系统和设备 |
CN115982958A (zh) * | 2022-12-07 | 2023-04-18 | 南京工业大学 | 一种基于工程损伤力学的材料蠕变疲劳寿命预测方法 |
CN115982958B (zh) * | 2022-12-07 | 2023-10-13 | 南京工业大学 | 一种基于工程损伤力学的材料蠕变疲劳寿命预测方法 |
CN116052813A (zh) * | 2022-12-30 | 2023-05-02 | 南京工业大学 | 一种全数值驱动的实时蠕变疲劳损伤监测方法 |
CN116052813B (zh) * | 2022-12-30 | 2023-10-13 | 南京工业大学 | 一种全数值驱动的实时蠕变疲劳损伤监测方法 |
CN115964919A (zh) * | 2023-02-06 | 2023-04-14 | 南京工业大学 | 熔盐腐蚀-力交互作用下晶间裂纹萌生及扩展的预测方法 |
CN115964919B (zh) * | 2023-02-06 | 2023-11-10 | 南京工业大学 | 熔盐腐蚀-力交互作用下晶间裂纹萌生及扩展的预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109885874B (zh) | 2022-12-23 |
US20220026326A1 (en) | 2022-01-27 |
WO2020143284A1 (zh) | 2020-07-16 |
US20230103274A9 (en) | 2023-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109885874A (zh) | 一种基于abaqus的多轴蠕变疲劳预测方法 | |
Philippidis et al. | Life prediction methodology for GFRP laminates under spectrum loading | |
CN106840877B (zh) | 一种基于应力的多轴小裂纹全寿命预测方法 | |
Goldberg et al. | Implementation of an associative flow rule including hydrostatic stress effects into the high strain rate deformation analysis of polymer matrix composites | |
Hoc et al. | A procedure for identifying the plastic behavior of single crystals from the local response of polycrystals | |
Passipoularidis et al. | A study of factors affecting life prediction of composites under spectrum loading | |
Huang et al. | Analysis of fatigue life of PMMA at Different frequencies based on a new damage mechanics model | |
Sarfaraz et al. | Variable amplitude fatigue of adhesively-bonded pultruded GFRP joints | |
Mandegarian et al. | A general energy based fatigue failure criterion for the carbon epoxy composites | |
Harris | A parametric constant-life model for prediction of the fatigue lives of fibre-reinforced plastics | |
Kim et al. | A time‐integration method for the viscoelastic–viscoplastic analyses of polymers and finite element implementation | |
Mars et al. | Computing tire component durability via critical plane analysis | |
Ge et al. | Experimental investigation and multiscale modeling of VE damper considering chain network and ambient temperature influence | |
Muliana et al. | A time-integration algorithm for thermo-rheologically complex polymers | |
Erber et al. | Hysteresis and fatigue | |
Ha et al. | Multiscale fatigue damage model for CFRP laminates considering the effect of progressive interface debonding | |
Granados et al. | Fatigue prediction of composite materials, based on serial/parallel mixing theory | |
CN109388833A (zh) | 一种基于疲劳寿命的弹性元件结构优化设计方法 | |
Carrella-Payan et al. | Implementation of fatigue model for unidirectional laminate based on finite element analysis: theory and practice | |
CN108520158B (zh) | 一种循环双曲线弹塑性土体本构模型及其应用方法 | |
Ding et al. | A study of fatigue crack growth with changing loading direction | |
Pineda et al. | Multiscale model for progressive damage and failure of laminated composites using an explicit finite element method | |
Sabounchi et al. | Microplane model of cylindrical geometry for transversely isotropic polymer composites | |
Ha et al. | Fatigue life prediction of CFRP laminates with stress concentration lamina level failure criteria | |
CN116052813B (zh) | 一种全数值驱动的实时蠕变疲劳损伤监测方法 |
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 |