CN113051719B - 一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法 - Google Patents
一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法 Download PDFInfo
- Publication number
- CN113051719B CN113051719B CN202110259040.8A CN202110259040A CN113051719B CN 113051719 B CN113051719 B CN 113051719B CN 202110259040 A CN202110259040 A CN 202110259040A CN 113051719 B CN113051719 B CN 113051719B
- Authority
- CN
- China
- Prior art keywords
- stress
- stage
- equation
- composite material
- matrix
- 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
- 239000011153 ceramic matrix composite Substances 0.000 title claims abstract description 121
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000000463 material Substances 0.000 claims abstract description 85
- 238000007906 compression Methods 0.000 claims abstract description 56
- 230000006835 compression Effects 0.000 claims abstract description 51
- 238000003825 pressing Methods 0.000 claims abstract description 28
- 239000002131 composite material Substances 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims description 128
- 239000000835 fiber Substances 0.000 claims description 95
- 238000009826 distribution Methods 0.000 claims description 52
- 238000005336 cracking Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 4
- 238000002360 preparation method Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 abstract description 11
- 238000004364 calculation method Methods 0.000 description 13
- 238000011160 research Methods 0.000 description 5
- 238000009940 knitting Methods 0.000 description 4
- 238000010008 shearing Methods 0.000 description 4
- 230000006978 adaptation Effects 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011226 reinforced ceramic Substances 0.000 description 1
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
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- 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
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)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法,属于复合材料疲劳迟滞回线预测技术领域。本发明首先将拉压过程的卸载阶段以及重新加载阶段分别划分为三个阶段,确定基体裂纹闭合应力以及完全压缩应力,分析编织陶瓷基复合材料在上述三个阶段的损伤细观应力场,根据断裂力学界面脱粘准则建立界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程以及界面新滑移长度方程,最终获得拉压过程三个阶段的应力‑应变关系方程,以此预测编织陶瓷基复合材料拉压疲劳迟滞回线。本发明提供的方法能够准确地预测拉压疲劳载荷对编织陶瓷基复合材料造成的损伤问题,提高了编织陶瓷基复合材料迟滞回线预测的准确性。
Description
技术领域
本发明涉及复合材料疲劳迟滞回线预测技术领域,尤其涉及一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法。
背景技术
编织陶瓷基复合材料具有耐高温、耐腐蚀、低密度、高比强、高比模等优点,相比高温合金,能够承受更高的温度,减少冷却气流,提高涡轮效率,目前已经应用于航空发动机燃烧室、涡轮导向叶片、涡轮壳环、尾喷管等。由CFM公司研制的LEAP(LeadingEdgeAviation Propulsion)系列发动机,高压涡轮采用了编织陶瓷基复合材料部件,LEAP-1B发动机为空客A320和波音737MAX提供动力,LEAP-X1C发动机为大型飞机C919提供动力。
为了保证编织陶瓷基复合材料在飞机和航空发动机结构中使用的可靠性与安全性,国内外研究人员将陶瓷基复合材料性能评估、损伤演化、强度与寿命预测工具的开发作为陶瓷基复合材料结构部件适航取证的关键。在拉压疲劳载荷作用下,编织陶瓷基复合材料出现基体多开裂、纤维/基体界面脱粘与滑移等多重损伤机制,使得卸载与重新加载过程中应力-应变曲线出现明显的迟滞现象。
目前针对编织陶瓷基复合材料疲劳迟滞回线的研究,未考虑压缩载荷对迟滞回线的影响(李龙彪,纤维增强陶瓷基复合材料疲劳迟滞回线模型研究[J],力学学报,2014,5:710-729)。如何考虑压缩载荷对编织陶瓷基复合材料疲劳迟滞回线的影响,监测压缩载荷对复合材料造成的损伤,是编织陶瓷基复合材料结构实际工程应用需要解决的关键技术问题。
发明内容
本发明的目的在于提供一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法,本发明提供的方法考虑了拉压疲劳载荷因素对编织陶瓷基复合材料界面脱粘以及界面滑移的影响,能够准确地预测拉压疲劳载荷对编织陶瓷基复合材料造成的损伤问题,提高了编织陶瓷基复合材料迟滞回线预测的准确性。
为了实现上述发明目的,本发明提供以下技术方案:
本发明提供了一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法,包括如下步骤:
(1)将编织陶瓷基复合材料的卸载阶段划分为第一卸载阶段、第二卸载阶段和第三卸载阶段,其中,所述第一卸载阶段为界面反向滑移阶段,第二卸载阶段为基体裂纹闭合阶段,第三卸载阶段为基体完全压缩阶段;
将编织陶瓷基复合材料的重新加载阶段划分为第一重新加载阶段、第二重新加载阶段和第三重新加载阶段,其中,所述第一重新加载阶段为基体完全压缩阶段,第二重新加载阶段为基体裂纹闭合阶段,第三重新加载阶段为界面新滑移阶段;
根据卸载与重新加载纤维相对基体在界面处的滑移机理,得到编织陶瓷基复合材料拉压疲劳载荷条件下的疲劳损伤细观应力场,基于所述编织陶瓷基复合材料拉压疲劳载荷条件下的损伤细观应力场,分别建立第一卸载阶段、第二卸载阶段、第三卸载阶段、第一重新加载阶段、第二重新加载阶段和第三重新加载阶段时纤维轴向应力分布方程和基体轴向应力分布方程;
(2)根据断裂力学界面脱粘准则,分别建立界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程和界面新滑移长度方程;
(3)根据步骤(1)中第一卸载阶段、第二卸载阶段、第三卸载阶段、第一重新加载阶段、第二重新加载阶段和第三重新加载阶段时纤维轴向应力分布方程和基体轴向应力分布方程以及步骤(2)中界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程和界面新滑移长度方程,建立编织陶瓷基复合材料拉压疲劳载荷条件下的应力-应变关系方程,以此预测编织陶瓷基复合材料拉压疲劳迟滞回线。
优选地,所述步骤(1)中,在第一卸载阶段时,编织陶瓷基复合材料的纤维轴向应力分布方程如式1所示,基体轴向应力分布方程如式2所示:
其中,σf(x)为纤维轴向应力,σunloading为卸载应力,τi为界面脱粘区摩擦剪应力,rf为纤维半径,Lcounter_slip为界面反向滑移长度,Ldebonding为界面脱粘长度,Lcracking为基体裂纹间距,ρ为剪滞模型参数,χ为复合材料沿应力加载方向纤维有效体积含量系数,Vf为复合材料中纤维体积含量;σfo为界面粘结区纤维轴向应力,σmo为界面粘结区基体轴向应力,x为轴向取值,Vm为复合材料中基体体积含量,σm(x)为基体轴向应力;
定义基体裂纹闭合应力为σclosure,在所述基体裂纹闭合应力作用下,基体闭合应力时纤维应变和基体闭合应力时基体应变满足式3所示方程:
εf(σclosure)=εm(σclosure) 式3;
其中,εf(σclosure)为基体闭合应力时纤维应变,εm(σclosure)为基体闭合应力时基体应变;
当卸载应力低于所述基体裂纹闭合应力时,进入第二卸载阶段,此时编织陶瓷基复合材料的纤维轴向应力分布方程如式4所示,基体轴向应力分布方程如式5所示:
其中,Lcompressive为界面压缩滑移长度;
定义完全压缩应力为σcomplete_compressive,当卸载应力低于所述完全压缩应力时,进入第三卸载阶段,此时编织陶瓷基复合材料的纤维轴向应力分布方程如式6-1所示,基体轴向应力分布方程如式6-2所示:
σf(x)=σfo 式6-1;
σm(x)=σmo 式6-2。
优选地,所述步骤(1)中,当重新加载应力低于所述完全压缩应力时,进入第一重新加载阶段,此时编织陶瓷基复合材料的纤维轴向应力由式6-1确定,基体轴向应力由式6-2确定;
当重新加载应力高于所述完全压缩应力且低于所述基体裂纹闭合应力时,进入所述第二重新加载阶段,此时编织陶瓷基复合材料的纤维轴向应力由式4确定,基体轴向应力由式5确定;
当重新加载应力高于所述基体裂纹闭合应力时,进入所述第三重新加载阶段,此时编织陶瓷基复合材料的纤维轴向应力分布方程如式7所示,基体轴向应力分布方程如式8所示:
其中,σreloading为重新加载应力,Lnew_slip为界面新滑移长度。
优选地,所述步骤(2)中,断裂力学界面脱粘准则满足式9所示方程:
其中,Γd为界面脱粘能,wf(0)为基体裂纹平面纤维轴向位移,v(x)为纤维相对基体轴向位移,F为基体裂纹平面纤维承担载荷。
优选地,所述步骤(2)中,界面脱粘长度方程如式10所示:
其中,Em为基体弹性模量,σmax为疲劳峰值应力,Ec为复合材料弹性模量,Ef为纤维弹性模量;
界面反向滑移长度方程如式11所示:
界面压缩滑移长度方程如式12所示:
其中,σclosure为基体裂纹闭合应力;
界面新滑移长度方程如式13所示:
优选地,所述步骤(3)中,在第一卸载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程如式14所示:
其中,εunloading为卸载应变,αc为复合材料热膨胀系数,αf为纤维热膨胀系数,ΔT为测试温度与制备温度差值,η为界面脱粘比,γ为界面反向滑移比;
所述η由式15-1确定,所述γ由式15-2确定:
在第二卸载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程如式16所示:
其中,δ为界面压缩滑移比;
所述δ由式17确定:
在第三卸载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程如式18所示:
优选地,所述步骤(3)中,在第一重新加载阶段时,所述编织陶瓷基复合材料的应力-应变关系由式18确定;
在第二重新加载阶段时,所述编织陶瓷基复合材料的应力-应变关系由式16确定;
在第三重新加载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程如式19所示:
其中,εreloading为重新加载应变,σ为应力,为界面新滑移比;
所述由式20确定:
优选地,所述编织陶瓷基复合材料沿应力加载方向纤维有效体积含量系数满足式21所示公式:
其中,Vf_loading为复合材料沿应力加载方向纤维体积含量。
本发明提供了一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法,本发明首先将拉压过程的卸载阶段以及重新加载阶段分别划分为三个阶段,确定基体裂纹闭合应力以及完全压缩应力,分析编织陶瓷基复合材料在上述三个阶段的损伤细观应力场,根据断裂力学界面脱粘准则建立界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程以及界面新滑移长度方程,最终获得拉压过程三个阶段的应力-应变关系方程,以此预测编织陶瓷基复合材料拉压疲劳迟滞回线。本发明提供的方法考虑了拉压疲劳载荷因素对编织陶瓷基复合材料界面脱粘以及界面滑移的影响,能够准确地预测拉压疲劳载荷对编织陶瓷基复合材料造成的损伤问题,提高了编织陶瓷基复合材料迟滞回线预测的准确性。
附图说明
图1为本发明中编织陶瓷基复合材料拉压疲劳迟滞回线;
图2为本发明中试验和理论预测的编织陶瓷基复合材料拉压疲劳迟滞回线。
具体实施方式
本发明提供的编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法中涉及的各项符号、含义及其获取方法汇总于表1中,以下具体实施方式中,除特殊说明外,各个方程或者关系式中的符号含义、获取方法均以表1的内容为准,不再一一赘述。
表1编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法中参数说明
注:表1中复合材料表示编织陶瓷基复合材料,纤维表示编织陶瓷基复合材料中的纤维,基体表示编织陶瓷复合材料中的基体,轴向指的是应力加载方向,界面指的是基体/纤维界面。
基于表1的说明,对本发明提供的编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法的具体实施过程进行如下说明:
本发明提供了一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法,包括如下步骤:
(1)将编织陶瓷基复合材料的卸载阶段划分为第一卸载阶段、第二卸载阶段和第三卸载阶段,其中,所述第一卸载阶段为界面反向滑移阶段,第二卸载阶段为基体裂纹闭合阶段,第三卸载阶段为基体完全压缩阶段;
将编织陶瓷基复合材料的重新加载阶段划分为第一重新加载阶段、第二重新加载阶段和第三重新加载阶段,其中,所述第一重新加载阶段为基体完全压缩阶段,第二重新加载阶段为基体裂纹闭合阶段,第三重新加载阶段为界面新滑移阶段;
根据卸载与重新加载纤维相对基体在界面处的滑移机理,得到编织陶瓷基复合材料拉压疲劳载荷条件下的损伤细观应力场,基于所述编织陶瓷基复合材料拉压疲劳载荷条件下的损伤细观应力场,分别建立第一卸载阶段、第二卸载阶段、第三卸载阶段、第一重新加载阶段、第二重新加载阶段和第三重新加载阶段时纤维轴向应力分布方程和基体轴向应力分布方程;
(2)根据断裂力学界面脱粘准则,分别建立界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程和界面新滑移长度方程;
(3)根据步骤(1)中第一卸载阶段、第二卸载阶段、第三卸载阶段、第一重新加载阶段、第二重新加载阶段和第三重新加载阶段时纤维轴向应力分布方程和基体轴向应力分布方程以及步骤(2)中界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程和界面新滑移长度方程,建立编织陶瓷基复合材料拉压疲劳载荷条件下的应力-应变关系方程,以此预测编织陶瓷基复合材料拉压疲劳迟滞回线。
本发明考虑拉压疲劳载荷(即包括卸载阶段以及重新加载阶段)因素的影响,首先将拉压过程的卸载阶段以及重新加载阶段分别划分为三个阶段,确定基体裂纹闭合应力以及完全压缩应力,分析编织陶瓷基复合材料在上述三个阶段的损伤细观应力场,根据断裂力学界面脱粘准则建立界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程以及界面新滑移长度方程,最终获得拉压过程三个阶段的应力-应变关系方程,以此预测编织陶瓷基复合材料拉压疲劳迟滞回线。
本发明将编织陶瓷基复合材料的卸载阶段划分为第一卸载阶段、第二卸载阶段和第三卸载阶段,其中,所述第一卸载阶段为界面反向滑移阶段,第二卸载阶段为基体裂纹闭合阶段,第三卸载阶段为基体完全压缩阶段;将编织陶瓷基复合材料的重新加载阶段划分为第一重新加载阶段、第二重新加载阶段和第三重新加载阶段,其中,所述第一重新加载阶段为基体完全压缩阶段,第二重新加载阶段为基体裂纹闭合阶段,第三重新加载阶段为界面新滑移阶段;根据卸载与重新加载纤维相对基体在界面处的滑移机理,得到编织陶瓷基复合材料拉压疲劳载荷条件下的损伤细观应力场,基于所述编织陶瓷基复合材料拉压疲劳载荷条件下的损伤细观应力场,分别建立第一卸载阶段、第二卸载阶段、第三卸载阶段、第一重新加载阶段、第二重新加载阶段和第三重新加载阶段时纤维轴向应力分布方程和基体轴向应力分布方程。
在本发明中,在第一卸载阶段时,编织陶瓷基复合材料的纤维轴向应力分布方程优选如式1所示,基体轴向应力分布方程优选如式2所示:
其中,σf(x)为纤维轴向应力,σunloading为卸载应力,τi为界面脱粘区摩擦剪应力,rf为纤维半径,Lcounter_slip为界面反向滑移长度,Ldebonding为界面脱粘长度,Lcracking为基体裂纹间距,ρ为剪滞模型参数,χ为复合材料沿应力加载方向纤维有效体积含量系数,Vf为复合材料中纤维体积含量;σfo为界面粘结区纤维轴向应力,σmo为界面粘结区基体轴向应力,x为轴向取值,Vm为复合材料中基体体积含量,σm(x)为基体轴向应力。
如式1和式2所示,本发明在研究编织陶瓷基复合材料在第一卸载阶段时的纤维轴向应力分布以及基体轴向应力分布时,优选研究基体裂纹区域段(0~Lcracking/2),更优选将这一区域段划分为界面反向滑移区(0~Lcounter_slip)、界面滑移区(Lcounter_slip~Ldebonding)和界面粘结区(Ldebonding~Lcracking/2),并针对不同区域,提出不同的计算方法,以提高第一卸载阶段时纤维轴向应力分布以及基体轴向应力分布预测结果的准确性。当x=Lcounter_slip时,纤维轴向应力以及基体轴向应力优选代入针对界面反向滑移区的公式中进行计算;当x=Ldebonding时,纤维轴向应力以及基体轴向应力优选代入针对界面滑移区的公式中进行计算。
本发明优选定义基体裂纹闭合应力为σclosure,在所述基体裂纹闭合应力作用下,基体闭合应力时纤维应变和基体闭合应力时基体应变满足式3所示方程:
εf(σclosure)=εm(σclosure) 式3;
其中,εf(σclosure)为基体闭合应力时纤维应变,εm(σclosure)为基体闭合应力时基体应变;
当卸载应力低于所述基体裂纹闭合应力时,进入第二卸载阶段,此时编织陶瓷基复合材料的纤维轴向应力分布方程优选如式4所示,基体轴向应力分布方程优选如式5所示:
其中,Lcompressive为界面压缩滑移长度。
如式4和式5所示,本发明在研究编织陶瓷基复合材料在第二卸载阶段时的纤维轴向应力分布以及基体轴向应力分布时,优选研究基体裂纹区域段(0~Lcracking/2),更优选将这一区域段划分为界面压缩区(0~Lcompressive)、界面反向滑移区(Lcompressive~Lcounter_slip)、界面滑移区(Lcounter_slip~Ldebonding)和界面粘结区(Ldebonding~Lcracking/2),并针对不同区域,提出不同的计算方法,以提高第二卸载阶段时纤维轴向应力分布以及基体轴向应力分布预测结果的准确性。当x=Lcompressive时,纤维轴向应力以及基体轴向应力优选代入针对界面压缩区的公式中进行计算;当x=Lcounter_slip时,纤维轴向应力以及基体轴向应力优选带入针对界面反向滑移区的公式中进行计算;当x=Ldebonding时,纤维轴向应力以及基体轴向应力优选带入针对界面脱粘区的公式中进行计算。
本发明优选定义完全压缩应力为σcomplete_compressive,当卸载应力低于所述完全压缩应力时(此时的卸载应力低于基体裂纹闭合应力),进入第三卸载阶段,此时编织陶瓷基复合材料的纤维轴向应力分布方程如式6-1所示,基体轴向应力分布方程如式6-2所示:
σf(x)=σfo 式6-1;
σm(x)=σmo 式6-2。
如式6-1和式6-2所示,在本发明中,在第三卸载阶段,编织陶瓷基复合材料的纤维轴向应力为界面粘结区纤维轴向应力,基体轴向应力为界面粘结区基体轴向应力。
在本发明中,当重新加载应力低于所述完全压缩应力时,进入第一重新加载阶段,此时编织陶瓷基复合材料的纤维轴向应力优选由式6-1确定,基体轴向应力由优选式6-2确定,在此不再赘述。
在本发明中,当重新加载应力高于所述完全压缩应力且低于所述基体裂纹闭合应力时,进入所述第二重新加载阶段,此时编织陶瓷基复合材料的纤维轴向应力由优选式4确定,基体轴向应力优选由式5确定,在此不再赘述。
在本发明中,当重新加载应力高于所述基体裂纹闭合应力时,进入所述第三重新加载阶段,此时编织陶瓷基复合材料的纤维轴向应力分布方程优选如式7所示,基体轴向应力分布方程优选如式8所示:
其中,σreloading为重新加载应力,Lnew_slip为界面新滑移长度。
如式7和式8所示,本发明在研究编织陶瓷基复合材料在第三重新加载阶段时的纤维轴向应力分布以及基体轴向应力分布时,优选研究基体裂纹区域段(0~Lcracking/2),更优选将这一区域段划分为界面新滑移区(0~Lnew_slip)、界面反向滑移区(Lnew_slip~Lcounter_slip)、界面滑移区(Lcounter_slip~Ldebonding)和界面粘结区(Ldebonding~Lcracking/2),并针对不同区域,提出不同的计算方法,以提高第三重新加载阶段时纤维轴向应力分布以及基体轴向应力分布预测结果的准确性。当x=Lnew_slip时,纤维轴向应力以及基体轴向应力优选代入针对界面新滑移区的公式中进行计算;当x=Lcounter_slip时,纤维轴向应力以及基体轴向应力优选代入针对界面反向滑移区的公式中进行计算;当x=Ldebonding时,纤维轴向应力以及基体轴向应力优选代入针对界面滑移区的公式中进行计算。
本发明根据断裂力学界面脱粘准则,分别建立界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程和界面新滑移长度方程。在本发明中,断裂力学界面脱粘准则优选满足式9所示方程:
/>
其中,Γd为界面脱粘能,wf(0)为基体裂纹平面纤维轴向位移,v(x)为纤维相对基体轴向位移,F为基体裂纹平面纤维承担载荷。
在本发明中,所述界面脱粘长度方程优选如式10所示:
其中,Em为基体弹性模量,σmax为疲劳峰值应力,Ec为复合材料弹性模量,Ef为纤维弹性模量;
界面反向滑移长度方程优选如式11所示:
界面压缩滑移长度方程如式12所示:
其中,σclosure为基体裂纹闭合应力;
界面新滑移长度方程优选如式13所示:
本发明通过式10~式13所示界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程以及界面新滑移长度方程,揭示了界面卸载与重新加载过程中的脱粘与滑移范围。
本发明根据所述第一卸载阶段、第二卸载阶段、第三卸载阶段、第一重新加载阶段、第二重新加载阶段和第三重新加载阶段时纤维轴向应力分布方程和基体轴向应力分布方程以及所述界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程和界面新滑移长度方程,建立编织陶瓷基复合材料拉压疲劳载荷条件下的应力-应变关系方程,以此预测编织陶瓷基复合材料拉压疲劳迟滞回线。
在本发明中,在第一卸载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程优选如式14所示:
其中,εunloading为卸载应变,αc为复合材料热膨胀系数,αf为纤维热膨胀系数,ΔT为测试温度与制备温度差值,η为界面脱粘比,γ为界面反向滑移比;
所述η优选由式15-1确定,所述γ优选由式15-2确定:
在本发明中,在第二卸载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程优选如式16所示:
其中,δ为界面压缩滑移比;
所述δ优选由式17确定:
在本发明中,在第三卸载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程优选如式18所示:
在本发明中,在第一重新加载阶段时,所述编织陶瓷基复合材料的应力-应变关系由式18确定,在此不再赘述。
在本发明中,在第二重新加载阶段时,所述编织陶瓷基复合材料的应力-应变关系由式16确定,在此不再赘述。
在本发明中,在第三重新加载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程优选如式19所示:
其中,εreloading为重新加载应变,σ为应力,为界面新滑移比;
所述优选由式20确定:
/>
在本发明中,上述技术方案涉及的公式中,编织陶瓷基复合材料沿应力加载方向纤维有效体积含量系数(χ)优选满足式21所示公式:
其中,Vf_loading为复合材料沿应力加载方向纤维体积含量。
在本发明中,所述复合材料沿应力加载方向纤维有效体积含量系数(χ)与编织陶瓷基复合材料中纤维的编织维度相关:
所述编织陶瓷基复合材料的编织维度为2时,χ为0.5;
所述编织陶瓷基复合材料的编织维度为2.5时,χ为0.75;
所述编织陶瓷基复合材料的编织维度为3时,χ为0.93。
在本发明中,上述技术方案涉及的公式中,剪滞模型参数(ρ)优选通过剪滞模型计算得到,所述剪滞模型优选为BHE剪滞模型。本发明对所述计算方式没有特殊要求,采用本领域技术人员熟知的方式即可。
本发明建立编织陶瓷基复合材料拉压疲劳载荷条件下的应力-应变关系方程,得到第一卸载阶段、第二卸载阶段、第三卸载阶段、第一重新加载阶段、第二重新加载阶段和第三重新加载阶段时应力随应变的变化曲线,以此得到编织陶瓷基复合材料拉压疲劳载荷条件下的迟滞回线,可以监测拉压疲劳载荷条件对编织陶瓷基复合材料造成的损伤,提高编织陶瓷基复合材料结构实际工程应用过程中的安全性。具体的,图1为本发明中编织陶瓷基复合材料拉压疲劳迟滞回线,按图1所示方式将卸载阶段以及重新加载阶段各划分为三个阶段,通过确定基体裂纹闭合应力以及完全压缩应力,分析编织陶瓷基复合材料在上述三个阶段的损伤细观应力场,根据断裂力学界面脱粘准则建立界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程以及界面新滑移长度方程,最终获得拉压过程三个阶段的应力-应变关系方程,以此能够实现编织陶瓷基复合材料拉压疲劳迟滞回线的准确预测。
下面将结合本发明中的实施例,对本发明中的技术方案进行清楚、完整地描述。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
按照本发明提供的预测方法建立所需的应力-应变关系方程,以编织陶瓷基复合材料(SiC/CAS)为测试样品,对测试样品进行拉压疲劳载荷试验,并对其拉压疲劳载荷迟滞回线进行预测:
提供参数:Ef=220GPa,Em=90GPa,rf=7.5μm,αf=4×10-6/℃,αm=4.6×10-6/℃,τi=20MPa,Γd=0.1J/m2,Vf=35%;
该复合材料的制备温度为1020℃,测试温度为20℃,则ΔT=-1000℃;
然后根据式1、式2、式4、式5、式6-1、式6-2、式7、式8、式10、式11、式12和式13所示方程,建立编织陶瓷基复合材料拉压疲劳载荷条件下的应力-应变关系方程,以此得到应力-应变关系,进而构建如图2所示的应力-应变关系曲线,得到编织陶瓷基复合材料拉压疲劳载荷迟滞回线。
图2中实线为采用本方案构建的应力-应变关系曲线,不同的点为实际测试数据,可见采用本发明提供的方法预测的迟滞回线形状和位置与实验数据相吻合,说明本发明提供的方法能够准确预测编织陶瓷基复合材料拉压疲劳载荷迟滞回线。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (6)
1.一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法,包括如下步骤:
(1)将编织陶瓷基复合材料的卸载阶段划分为第一卸载阶段、第二卸载阶段和第三卸载阶段,其中,所述第一卸载阶段为界面反向滑移阶段,第二卸载阶段为基体裂纹闭合阶段,第三卸载阶段为基体完全压缩阶段;
将编织陶瓷基复合材料的重新加载阶段划分为第一重新加载阶段、第二重新加载阶段和第三重新加载阶段,其中,所述第一重新加载阶段为基体完全压缩阶段,第二重新加载阶段为基体裂纹闭合阶段,第三重新加载阶段为界面新滑移阶段;
根据卸载与重新加载纤维相对基体在界面处的滑移机理,得到编织陶瓷基复合材料拉压疲劳载荷条件下的疲劳损伤细观应力场,基于所述编织陶瓷基复合材料拉压疲劳载荷条件下的损伤细观应力场,分别建立第一卸载阶段、第二卸载阶段、第三卸载阶段、第一重新加载阶段、第二重新加载阶段和第三重新加载阶段时纤维轴向应力分布方程和基体轴向应力分布方程;
(2)根据断裂力学界面脱粘准则,分别建立界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程和界面新滑移长度方程;
(3)根据步骤(1)中第一卸载阶段、第二卸载阶段、第三卸载阶段、第一重新加载阶段、第二重新加载阶段和第三重新加载阶段时纤维轴向应力分布方程和基体轴向应力分布方程以及步骤(2)中界面脱粘长度方程、界面反向滑移长度方程、界面压缩滑移长度方程和界面新滑移长度方程,建立编织陶瓷基复合材料拉压疲劳载荷条件下的应力-应变关系方程,以此预测编织陶瓷基复合材料拉压疲劳迟滞回线;
所述步骤(3)中,在第一卸载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程如式14所示:
其中,εunloading为卸载应变,αc为复合材料热膨胀系数,αf为纤维热膨胀系数,ΔT为测试温度与制备温度差值,η为界面脱粘比,γ为界面反向滑移比,σunloading为卸载应力,Vf为复合材料中纤维体积含量,Ef为纤维弹性模量,τi为界面脱粘区摩擦剪应力,Ldebonding为界面脱粘长度,rf为纤维半径,σfo为界面粘结区纤维轴向应力,ρ为剪滞模型参数,σmo为界面粘结区基体轴向应力,Vm为复合材料中基体体积含量,Lcracking为基体裂纹间距;
所述η由式15-1确定,所述γ由式15-2确定:
其中,Lcounter_slip为界面反向滑移长度;
在第二卸载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程如式16所示:
其中,δ为界面压缩滑移比,σclosure为基体裂纹闭合应力,
所述δ由式17确定:
其中,Lcompressive为界面压缩滑移长度;
在第三卸载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程如式18所示:
其中,Ec为复合材料弹性模量;
所述步骤(3)中,在第一重新加载阶段时,所述编织陶瓷基复合材料的应力-应变关系由式18确定;
在第二重新加载阶段时,所述编织陶瓷基复合材料的应力-应变关系由式16确定;
在第三重新加载阶段时,所述编织陶瓷基复合材料的应力-应变关系方程如式19所示:
其中,εreloading为重新加载应变,σ为应力,为界面新滑移比;
所述由式20确定:
2.根据权利要求1所述的预测方法,其特征在于,所述步骤(1)中,在第一卸载阶段时,编织陶瓷基复合材料的纤维轴向应力分布方程如式1所示,基体轴向应力分布方程如式2所示:
其中,σf(x)为纤维轴向应力,χ为复合材料沿应力加载方向纤维有效体积含量系数,x为轴向取值,σm(x)为基体轴向应力;
在所述基体裂纹闭合应力作用下,基体闭合应力时纤维应变和基体闭合应力时基体应变满足式3所示方程:
εf(σclosure)=εm(σclosure) 式3;
其中,εf(σclosure)为基体闭合应力时纤维应变,εm(σclosure)为基体闭合应力时基体应变;
当卸载应力低于所述基体裂纹闭合应力时,进入第二卸载阶段,此时编织陶瓷基复合材料的纤维轴向应力分布方程如式4所示,基体轴向应力分布方程如式5所示:
定义完全压缩应力为σcomplete_compressive,当卸载应力低于所述完全压缩应力时,进入第三卸载阶段,此时编织陶瓷基复合材料的纤维轴向应力分布方程如式6-1所示,基体轴向应力分布方程如式6-2所示:
σf(x)=σfo 式6-1;
σm(x)=σmo 式6-2。
3.根据权利要求2所述的预测方法,其特征在于,所述步骤(1)中,当重新加载应力低于所述完全压缩应力时,进入第一重新加载阶段,此时编织陶瓷基复合材料的纤维轴向应力由式6-1确定,基体轴向应力由式6-2确定;
当重新加载应力高于所述完全压缩应力且低于所述基体裂纹闭合应力时,进入所述第二重新加载阶段,此时编织陶瓷基复合材料的纤维轴向应力由式4确定,基体轴向应力由式5确定;
当重新加载应力高于所述基体裂纹闭合应力时,进入所述第三重新加载阶段,此时编织陶瓷基复合材料的纤维轴向应力分布方程如式7所示,基体轴向应力分布方程如式8所示:
其中,σreloading为重新加载应力,Lnew_slip为界面新滑移长度。
4.根据权利要求3所述的预测方法,其特征在于,所述步骤(2)中,断裂力学界面脱粘准则满足式9所示方程:
其中,Γd为界面脱粘能,wf(0)为基体裂纹平面纤维轴向位移,v(x)为纤维相对基体轴向位移,F为基体裂纹平面纤维承担载荷。
5.根据权利要求4所述的预测方法,其特征在于,所述步骤(2)中,界面脱粘长度方程如式10所示:
其中,Em为基体弹性模量,σmax为疲劳峰值应力;
界面反向滑移长度方程如式11所示:
界面压缩滑移长度方程如式12所示:
其中,σclosure为基体裂纹闭合应力;
界面新滑移长度方程如式13所示:
6.根据权利要求2~5任一项所述的预测方法,其特征在于,所述编织陶瓷基复合材料沿应力加载方向纤维有效体积含量系数满足式21所示公式:
其中,Vf_loading为复合材料沿应力加载方向纤维体积含量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110259040.8A CN113051719B (zh) | 2021-03-10 | 2021-03-10 | 一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110259040.8A CN113051719B (zh) | 2021-03-10 | 2021-03-10 | 一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113051719A CN113051719A (zh) | 2021-06-29 |
CN113051719B true CN113051719B (zh) | 2024-05-31 |
Family
ID=76511245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110259040.8A Active CN113051719B (zh) | 2021-03-10 | 2021-03-10 | 一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113051719B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114139388A (zh) * | 2021-12-03 | 2022-03-04 | 南京航空航天大学 | 一种预测纤维增强陶瓷基复合材料非闭合迟滞回线的方法 |
CN115031948B (zh) * | 2022-08-10 | 2022-11-11 | 常州市宏发纵横新材料科技股份有限公司 | 模块化叶片疲劳寿命预测方法、装置、控制装置和存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109632887A (zh) * | 2018-12-29 | 2019-04-16 | 南京航空航天大学 | 一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法 |
CN109632530A (zh) * | 2018-12-29 | 2019-04-16 | 南京航空航天大学 | 一种编织陶瓷基复合材料热机械疲劳迟滞回线的预测方法 |
CN109670272A (zh) * | 2019-01-23 | 2019-04-23 | 南京航空航天大学 | 一种编织陶瓷基复合材料过载疲劳迟滞回线的预测方法 |
CN109684598A (zh) * | 2018-12-29 | 2019-04-26 | 南京航空航天大学 | 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳寿命的方法 |
CN110196996A (zh) * | 2019-04-30 | 2019-09-03 | 南京航空航天大学 | 一种金属基复合材料拉压疲劳迟滞回线预测方法 |
CN111474062A (zh) * | 2020-04-23 | 2020-07-31 | 南京航空航天大学 | 一种编织陶瓷基复合材料高温静疲劳损伤演化预测方法 |
CN111474044A (zh) * | 2020-04-24 | 2020-07-31 | 南京航空航天大学 | 一种考虑多级疲劳损伤的编织陶瓷基复合材料疲劳迟滞回线的预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109614755B (zh) * | 2018-12-29 | 2023-04-07 | 南京航空航天大学 | 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳纤维/基体界面剪应力的方法 |
-
2021
- 2021-03-10 CN CN202110259040.8A patent/CN113051719B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109632887A (zh) * | 2018-12-29 | 2019-04-16 | 南京航空航天大学 | 一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法 |
CN109632530A (zh) * | 2018-12-29 | 2019-04-16 | 南京航空航天大学 | 一种编织陶瓷基复合材料热机械疲劳迟滞回线的预测方法 |
CN109684598A (zh) * | 2018-12-29 | 2019-04-26 | 南京航空航天大学 | 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳寿命的方法 |
CN109670272A (zh) * | 2019-01-23 | 2019-04-23 | 南京航空航天大学 | 一种编织陶瓷基复合材料过载疲劳迟滞回线的预测方法 |
CN110196996A (zh) * | 2019-04-30 | 2019-09-03 | 南京航空航天大学 | 一种金属基复合材料拉压疲劳迟滞回线预测方法 |
CN111474062A (zh) * | 2020-04-23 | 2020-07-31 | 南京航空航天大学 | 一种编织陶瓷基复合材料高温静疲劳损伤演化预测方法 |
CN111474044A (zh) * | 2020-04-24 | 2020-07-31 | 南京航空航天大学 | 一种考虑多级疲劳损伤的编织陶瓷基复合材料疲劳迟滞回线的预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113051719A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109614755B (zh) | 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳纤维/基体界面剪应力的方法 | |
Hirano | Application of eutectic composites to gas turbine system and fundamental fracture properties up to 1700° C | |
CN113051719B (zh) | 一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法 | |
CN109670272B (zh) | 一种编织陶瓷基复合材料过载疲劳迟滞回线的预测方法 | |
CN111523237B (zh) | 一种编织陶瓷基复合材料蠕变断裂行为的预测方法 | |
WO2018171096A1 (zh) | 一种设计考核高温旋转构件寿命与强度的双判据图方法 | |
CN113032982B (zh) | 一种考虑基体与纤维碎断的编织陶瓷基复合材料疲劳迟滞回线的预测方法 | |
CN111024486B (zh) | 一种单向陶瓷基复合材料蠕变行为预测方法 | |
CN109684598B (zh) | 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳寿命的方法 | |
CN111474044B (zh) | 一种考虑多级疲劳损伤的编织陶瓷基复合材料疲劳迟滞回线的预测方法 | |
CN109781546B (zh) | 一种编织陶瓷基复合材料拉伸强度的预测方法 | |
CN109711100B (zh) | 一种编织陶瓷基复合材料疲劳寿命的预测方法 | |
CN111474043B (zh) | 一种考虑多级疲劳损伤的编织陶瓷基复合材料剩余强度的预测方法 | |
Li | High temperature mechanical behavior of ceramic-matrix composites | |
CN113051720B (zh) | 一种考虑基体碎断的编织陶瓷基复合材料裂纹张开位移的预测方法 | |
Holmes et al. | Elevated temperature creep behavior of continuous fiber-reinforced ceramics | |
Ye et al. | A new creep model and its application in the evaluation of creep properties of a titanium alloy at 500 C | |
CN111474062A (zh) | 一种编织陶瓷基复合材料高温静疲劳损伤演化预测方法 | |
CN113032983B (zh) | 一种考虑纤维断裂影响的编织陶瓷基复合材料振动阻尼的预测方法 | |
CN109616163B (zh) | 一种考虑环境影响的编织陶瓷基复合材料多基体开裂的预测方法 | |
CN114970158B (zh) | 一种编织陶瓷基复合材料主裂纹张开位移预测方法及*** | |
Sanders et al. | Longitudinal fatigue response of a metal matrix composite under strain controlled mode at elevated temperature | |
Carrere et al. | Multi-scale modelling of silicon carbide reinforced titanium MMCs: Application to advanced compressor design | |
Mahadevan et al. | Probabilistic fatigue–creep life prediction of composites | |
Siadkowska et al. | Selecting a material for an aircraft diesel engine block |
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 |