CN109632887B - 一种监测编织陶瓷基复合材料高温疲劳损伤累积的方法 - Google Patents

一种监测编织陶瓷基复合材料高温疲劳损伤累积的方法 Download PDF

Info

Publication number
CN109632887B
CN109632887B CN201811642683.5A CN201811642683A CN109632887B CN 109632887 B CN109632887 B CN 109632887B CN 201811642683 A CN201811642683 A CN 201811642683A CN 109632887 B CN109632887 B CN 109632887B
Authority
CN
China
Prior art keywords
fiber
equation
matrix
stress
matrix interface
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
Application number
CN201811642683.5A
Other languages
English (en)
Other versions
CN109632887A (zh
Inventor
李龙彪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201811642683.5A priority Critical patent/CN109632887B/zh
Publication of CN109632887A publication Critical patent/CN109632887A/zh
Application granted granted Critical
Publication of CN109632887B publication Critical patent/CN109632887B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00Investigating or analyzing materials by the use of thermal means
    • G01N25/72Investigating presence of flaws

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

本发明属于复合材料高温疲劳损伤监测技术领域,具体涉及一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法。本发明利用与温度和循环数相关的纤维/基体界面剪应力建立编织陶瓷基复合材料的纤维/基体界面脱粘长度方程,并以此为基础,得到编织陶瓷基复合材料的疲劳耗散能方程,再通过疲劳耗散能方程得到疲劳损伤参数,用于监测编织陶瓷基复合材料的高温疲劳损伤累积。本发明提供的上述方法,充分考虑温度和循环次数对复合材料的基体和纤维纤维/基体界面的影响,所得复合材料的高温疲劳损伤累积更加准确。

Description

一种监测编织陶瓷基复合材料高温疲劳损伤累积的方法
技术领域
本发明属于复合材料疲劳寿命预测技术领域,具体涉及一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法。
背景技术
编织陶瓷基复合材料具有耐高温、耐腐蚀、低密度、高比强、高比模等优点,相比高温合金,能够承受更高的温度,减少冷却气流,提高涡轮效率,目前已经应用于航空发动机燃烧室、涡轮导向叶片、涡轮壳环、尾喷管等。由CFM公司研制的LEAP(Leading EdgeAviation Propulsion,LEAP)系列发动机,高压涡轮采用了编织陶瓷基复合材料部件,LEAP-1B发动机为空客A320和波音737MAX提供动力,LEAP-X1C发动机也是我国大型飞机C919选用的唯一动力装置。
为了保证编织陶瓷基复合材料在飞机和航空发动机结构中使用的可靠性与安全性,国内外研究学者均将陶瓷基复合材料性能评估、损伤演化、强度与寿命预测工具的开发作为陶瓷基复合材料结构部件适航取证的关键。编织陶瓷基复合材料的疲劳寿命对材料使用的安全性有直接影响,如何准确预测编织陶瓷基复合材料的高温疲劳损伤累积,是保证其使用可靠性与安全性的关键。
发明内容
本发明的目的在于提供一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法,本发明提供的方法将高温环境和循环次数的影响融入监测过程,提高了编织陶瓷基复合材料高温疲劳损伤累积监测结果的准确性。
为了实现上述目的,本发明提供如下技术方案:
本发明提供了一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法,包括以下步骤:
(1)根据断裂力学脱粘准则,利用与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力建立纤维/基体界面脱粘长度方程;
(2)根据断裂力学脱粘准则、纤维/基体界面滑移机理和所述步骤(1)得到的纤维/基体界面脱粘长度方程,建立卸载纤维/基体界面反向滑移长度方程;
(3)根据断裂力学脱粘准则、纤维/基体界面滑移机理、所述步骤(1)得到的纤维/基体界面脱粘长度方程和所述步骤(2)得到的卸载纤维/基体界面反向滑移长度方程,建立重新加载纤维/基体界面新滑移长度方程;
(4)根据所述步骤(1)得到的纤维/基体界面脱粘长度方程和所述步骤(2)得到的卸载纤维/基体界面反向滑移长度方程,结合总体载荷承担准则,建立卸载应力-应变方程;
根据所述步骤(1)得到的纤维/基体界面脱粘长度方程、所述步骤(2)得到的卸载纤维/基体界面反向滑移长度方程和所述步骤(3)得到的重新加载纤维/基体界面新滑移长度方程,结合总体载荷承担准则,建立重新加载应力-应变方程;
(5)根据所述步骤(4)得到的卸载应力-应变方程和重新加载应力-应变方程,建立疲劳迟滞耗散能方程,监测编织陶瓷基复合材料的高温疲劳损伤累积。
优选的,所述步骤(1)中纤维/基体界面脱粘长度方程如式1所示:
Figure GDA0002290954250000021
式1中,ld表示纤维/基体界面脱粘长度;
τi(T,N)表示与温度和循环数相关的纤维/基体界面脱粘区磨擦剪应力;
rf表示纤维半径;
Vm表示基体体积;
Em表示基体弹性模量;
σ表示应力;
χ表示沿应力加载方向纤维有效体积含量系数;
Vf表示编织陶瓷基复合材料中纤维体积含量;
Ec表示编织陶瓷基复合材料弹性模量;
ρ表示剪滞模型参数;
Ef表示纤维弹性模量;
ζd表示纤维/基体界面脱粘能。
优选的,所述步骤(1)中断裂力学脱粘准则如式1-1所示:
Figure GDA0002290954250000031
纤维轴向位移如式1-2所示:
Figure GDA0002290954250000032
所述纤维相对基体轴向位移如式1-3所示:
Figure GDA0002290954250000033
式1-1、1-2和1-3中,
F表示基体裂纹平面纤维承担载荷;
Figure GDA0002290954250000036
表示纤维轴向位移中的x为0时,对纤维/基体界面脱粘长度求导;
Figure GDA0002290954250000037
表示纤维相对基体轴向位移中的x为0时,对纤维/基体界面脱粘长度求导;
wf(x)表示纤维轴向位移;
v(x)表示纤维相对基体轴向位移;
lc表示基体裂纹间距。
优选的,所述卸载纤维/基体界面反向滑移长度方程如式2所示:
Figure GDA0002290954250000034
式2中,y表示卸载纤维/基体界面反向滑移长度;
σmax表示疲劳峰值应力。
优选的,所述重新加载纤维/基体界面新滑移长度方程如式3所示:
Figure GDA0002290954250000035
式3中,z表示重新加载纤维/基体界面新滑移长度;
σmin表示疲劳谷值应力。
优选的,纤维/基体界面部分脱粘时,所述卸载应力-应变方程如式4-1所示:
Figure GDA0002290954250000041
所述重新加载应力-应变方程优选如式4-2所示:
Figure GDA0002290954250000042
纤维/基体界面完全脱粘时,所述卸载应力-应变方程如式4-3所示:
Figure GDA0002290954250000043
所述重新加载应力-应变方程如式4-4所示:
Figure GDA0002290954250000044
式4-1、4-2、4-3和4-4中,
εunloading_p表示纤维/基体界面部分脱粘时卸载应力对应的应变;
εreloading_p表示纤维/基体界面部分脱粘时重新加载应力对应的应变;
εunloading_f表示纤维/基体界面完全脱粘时卸载应力对应的应变;
εreloading_f表示纤维/基体界面完全脱粘时重新加载应力对应的应变;
ac表示复合材料热膨胀系数;
af表示纤维热膨胀系数;
△T表示编织陶瓷基复合材料的测试温度与制备温度的差。
优选的,所述步骤(5)中疲劳迟滞耗散能方程如式5所示:
Figure GDA0002290954250000045
式5中,Un表示不同循环数的疲劳迟滞耗散能;
εunloading(σ)表示编织陶瓷基复合材料的卸载应变;
εreloading(σ)表示编织陶瓷剂复合材料的重新加载应变;
σmax表示疲劳峰值应力;
σmin表示疲劳谷值应力。
优选的,还包括根据所述步骤(5)得到的疲劳迟滞耗散能方程得到初始循环的疲劳迟滞耗散能和不同循环数的疲劳迟滞耗散能,结合编织陶瓷基复合材料的弹性应变能,得到编织陶瓷基复合材料的损伤参数方程,监测编织陶瓷基复合材料的高温疲劳损伤累积。
优选的,所述疲劳迟滞耗散能的损伤参数方程如式6所示:
Figure GDA0002290954250000051
式6中,Un表示不同循环数的疲劳迟滞耗散能;
Ue表示编织陶瓷基复合材料的弹性应变能;
Uinitial表示初始循环的疲劳迟滞耗散能;
Ф表示编织陶瓷基复合材料的损伤参数。
本发明利用与温度和循环数相关的纤维/基体界面剪应力建立编织陶瓷基复合材料的纤维/基体界面脱粘长度方程,并以此为基础,得到编织陶瓷基复合材料的疲劳耗散能方程,监测编织陶瓷基复合材料的高温疲劳损伤累积。本发明提供的上述方法,充分考虑温度和循环次数对复合材料的基体和纤维纤维/基体界面的影响,所得复合材料的高温疲劳损伤累积的监测结果更加准确。
附图说明
图1为本发明提供的编织陶瓷基复合材料应力状态下纤维开裂的剪滞单胞模型图;
图2为本发明提供的编织陶瓷基复合材料高温疲劳损伤累积关系曲线。
具体实施方式
本发明提供的通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法中涉及的各项符号、含义及其获取方法汇总于表1中,以下具体实施方式中,除特殊说明外,各个方程或者关系式中的符号含义、获取方法均以表1的内容为准。
在本发明中,所述高温指600℃以上的温度;纤维/基体界面指编织陶瓷基复合材料受损区域纤维和基体的纤维/基体界面;轴向指纤维轴向方向。
表1编织陶瓷基复合材料高温疲劳损伤累积监测方法的参数说明
Figure GDA0002290954250000052
Figure GDA0002290954250000061
Figure GDA0002290954250000071
为进一步清楚描述本发明所述的通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法,本发明优选提供编织陶瓷基复合材料的剪滞单胞模型图(图1所示),以对本发明出现的参数含义做进一步说明。
图1中,fober为纤维;matrix为基体;σ表示应力;ld表示纤维/基体界面脱粘长度,该长度对应的区域为纤维/基体界面脱粘区;lc表示基体裂纹间距,lc/2-ld对应的区域为纤维/基体界面粘结区;x表示轴向方向,w(x)表示沿轴向的位移,v(0)表示x为0时,纤维相对基体位移;σ/Vf表示基体裂纹平面纤维承担应力。
在应力作用下,编织陶瓷基复合材料受损区域的纤维、基体相对移动产生纤维/基体界面摩擦剪应力,以τi(x)表示;纤维/基体界面脱粘区对应纤维/基体界面脱粘区摩擦剪应力,以τi表示;编织陶瓷基复合材料纤维/基体界面摩擦剪应力受温度的影响,将与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力以τi(T,N)表示。
基于表1和图1的说明,对本发明提供的方法进行如下的说明:
本发明提供了一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法,包括以下步骤:
(1)根据断裂力学脱粘准则,利用与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力建立纤维/基体界面脱粘长度方程;
(2)根据断裂力学脱粘准则、纤维/基体界面滑移机理和所述步骤(1)得到的纤维/基体界面脱粘长度方程,建立卸载纤维/基体界面反向滑移长度方程;
(3)根据断裂力学脱粘准则、纤维/基体界面滑移机理、所述步骤(1)得到的纤维/基体界面脱粘长度方程和所述步骤(2)得到的卸载纤维/基体界面反向滑移长度方程,建立重新加载纤维/基体界面新滑移长度方程;
(4)根据所述步骤(1)得到的纤维/基体界面脱粘长度方程和所述步骤(2)得到的卸载纤维/基体界面反向滑移长度方程,结合总体载荷承担准则,建立卸载应力-应变方程;
根据所述步骤(1)得到的纤维/基体界面脱粘长度方程、所述步骤(2)得到的卸载纤维/基体界面反向滑移长度方程和所述步骤(3)得到的重新加载纤维/基体界面新滑移长度方程,结合总体载荷承担准则,建立重新加载应力-应变方程;
(5)根据所述步骤(4)得到的卸载应力-应变方程和重新加载应力-应变方程,建立疲劳迟滞耗散能方程,监测编织陶瓷基复合材料的高温疲劳损伤累积。
本发明根据断裂力学脱粘准则,利用与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力建立纤维/基体界面脱粘长度方程。
在本发明中,所述纤维/基体界面脱粘长度方程优选如式1所示:
Figure GDA0002290954250000081
式1中,ld表示纤维/基体界面脱粘长度;
τi(T,N)表示与温度和循环数相关的纤维/基体界面脱粘区磨擦剪应力;
rf表示纤维半径;
Vm表示基体体积;
Em表示基体弹性模量;
σ表示应力;
χ表示沿应力加载方向纤维有效体积含量系数;
Vf表示编织陶瓷基复合材料中纤维体积含量;
Ec表示编织陶瓷基复合材料弹性模量;
ρ表示剪滞模型参数;
Ef表示纤维弹性模量;
ζd表示纤维/基体界面脱粘能。
在本发明中,所述纤维/基体界面脱粘长度方程为纤维/基体界面脱粘长度、纤维/基体界面脱粘区摩擦剪应力和应力的关系式,其中,纤维/基体界面脱粘区摩擦剪应力为与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力,将循环因素和温度因素融入脱粘长度方程中,使监测过程接近复合材料实际应用条件,进而提高了监测结果的准确性。
在本发明中,构建脱粘长度方程所用断裂力学脱粘准则优选如式1-1所示:
Figure GDA0002290954250000091
式1-1中,
F表示基体裂纹平面纤维承担荷载;
Figure GDA0002290954250000093
表示纤维轴向位移中的x为0时,对纤维/基体界面脱粘长度求导;
Figure GDA0002290954250000094
表示纤维相对基体轴向位移中的x为0时,对纤维/基体界面脱粘长度求导;
在本发明中,纤维轴向位移优选如式1-2表示;
Figure GDA0002290954250000092
wf(x)表示纤维轴向位移;
τi(T,N)表示与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力;
x表示轴向坐标;
lc表示基体裂纹间距。
在本发明中,纤维相对基体轴向位移优选如式1-3所示:
Figure GDA0002290954250000101
式1-3中,v(x)表示纤维相对基体轴向位移。
在本发明中,编织陶瓷基复合材料受损区域在应力作用下,纤维和基体均会产生相对移动,其中纤维的移动距离以纤维轴向位移(wf(x))表示,基体的移动距离以基体轴向位移(wm(x))表示;所述基体轴向位移优选如式1-4所示:
Figure GDA0002290954250000102
在本发明中,纤维轴向位移和基体轴向位移差的绝对值即为纤维相对基体轴向位移。本发明优选通过式1-2和1-4得到如式1-3所示的纤维相对基体轴向位移的表达式。本发明优选将式1-2和1-3与式1-1结合,得到式1所示的纤维/基体界面脱粘长度方程。
在本发明中,式1中的沿应力加载方向纤维有效体积含量系数(χ)与编织陶瓷基复合材料中纤维的编织维度相关:
所述编织陶瓷基复合材料的编织维度为2时,χ为0.5;
所述编织陶瓷基复合材料的编织维度为2.5时,χ为0.75;
所述编织陶瓷基复合材料的编织维度为3时,χ为0.93。
在本发明具体实施例中,所述编织陶瓷基复合材料的维度优选为2。
在本发明中,式1中的剪滞模型参数(ρ)优选通过剪滞模型计算得到,所述剪滞模型优选为BHE剪滞模型。本发明对所述计算方式没有特殊要求,采用本领域技术人员熟知的即可。
在本发明中,与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力τi(T,N)优选通过式1-5得到;
Figure GDA0002290954250000103
式1-5中,μ为编织陶瓷基复合材料纤维/基体界面摩擦系数;
αrf为纤维径向热膨胀系数;
αrm为基体径向热膨胀系数;
T0为编织陶瓷基复合材料制备温度;
T为编织陶瓷基复合材料使用温度;
A表示编织陶瓷基复合材料弹性常数;
τi(T,N)表示与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力;
τ0_i(N)为纤维/基体界面脱粘区摩擦剪应力。
如式1-5,本发明所述与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力优选通过不同循环数的纤维/基体界面脱粘区摩擦剪应力、编织陶瓷基复合材料纤维/基体界面摩擦系数、纤维径向热膨胀系数、基体径向热膨胀系数、编织陶瓷基复合材料的制备温度和使用温度、编织陶瓷基复合材料弹性常数得到。
在本发明中,式1中的基体体积含量、纤维体积含量、纤维半径、基体弹性模量、纤维弹性模量、纤维/基体界面脱粘能、复合材料弹性模量优选按照表1提供的方式得到,本发明对上述参数的测量方法或者计算方法没有特殊要求,采用本领域技术人员熟知的即可。在本发明中,式1中的应力为加载在编织陶瓷基复合材料上的应力,优选通过测量得到。
在本发明中,所述纤维/基体界面脱粘长度优选是利用包括与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力的参数得到,上述参量受温度和/或循环数的影响,因此,纤维/基体界面脱粘长度是纳入了温度和循环次数影响因素的一个参量,再以纤维/基体界面脱粘长度这个参量为基础,构建后续的迟滞耗散能方程,使编织陶瓷基复合材料疲劳损伤参数的结果更接近实际情况,对提高编织陶瓷基复合材料的高温疲劳损伤累积监测值的准确性有利。
得到纤维/基体界面脱粘长度方程后,本发明根据断裂力学脱粘准则、纤维/基体界面滑移机理和所述的纤维/基体界面脱粘长度方程,建立卸载纤维/基体界面反向滑移长度方程。在本发明中,所述卸载纤维/基体界面反向滑移长度方程优选如式2所示:
Figure GDA0002290954250000111
式2中,y表示卸载纤维/基体界面反向滑移长度方程;
σmax表示疲劳峰值应力。
得到卸载纤维/基体界面反向滑移长度方程后,本发明根据断裂力学脱粘准则、纤维/基体界面滑移机理、所述的纤维/基体界面脱粘长度方程和所述的卸载纤维/基体界面反向滑移长度方程,建立重新加载纤维/基体界面新滑移长度方程。在本发明中,所述重新加载纤维/基体界面新滑移长度方程优选如式3所示:
Figure GDA0002290954250000121
式3中,z表示重新加载纤维/基体界面新滑移长度;
σmin表示疲劳谷值应力。
在本发明中,所述重新加载纤维/基体界面新滑移长度方程优选为包括重新加载纤维/基体界面新滑移长度、卸载纤维/基体界面反向滑移长度、纤维/基体界面脱粘区摩擦剪应力、脱粘长度和应力的关系式,为建立应力-应变方程提供基础。
得到重新加载纤维/基体界面新滑移长度方程后,本发明根据所述的纤维/基体界面脱粘长度方程和所述的卸载纤维/基体界面反向滑移长度方程,结合总体载荷承担准则,建立卸载应力-应变方程;
根据所述的纤维/基体界面脱粘长度方程、所述的卸载纤维/基体界面反向滑移长度方程和所述的重新加载纤维/基体界面新滑移长度方程,结合总体载荷承担准则,建立重新加载应力-应变方程。
在本发明中,所述卸载应力-应变方程和重新加载应力-应变方程优选分为纤维/基体界面部分脱粘和纤维/基体界面完全脱粘两种情况:
当纤维/基体界面部分脱粘时,所述卸载应力-应变方程如式4-1所示:
Figure GDA0002290954250000122
所述重新加载应力-应变方程优选如式4-2所示:
Figure GDA0002290954250000131
当纤维/基体界面完全脱粘时,所述卸载应力-应变方程如式4-3所示:
Figure GDA0002290954250000132
所述重新加载应力-应变方程优选如式4-4所示:
Figure GDA0002290954250000133
式4-1、4-2、4-3和4-4中,
εunloading_p表示纤维/基体界面部分脱粘时卸载应力对应的应变;
εreloading_p表示纤维/基体界面部分脱粘时重新加载应力对应的应变;
εunloading_f表示纤维/基体界面完全脱粘时卸载应力对应的应变;
εreloading_f表示纤维/基体界面完全脱粘时重新加载应力对应的应变;
ac表示复合材料热膨胀系数;
af表示纤维热膨胀系数;
△T表示编织陶瓷基复合材料的测试温度与制备温度的差。
得到应力-应变方程后,本发明根据所述的卸载应力-应变方程和重新加载应力-应变方程,建立疲劳迟滞耗散能方程,监测编织陶瓷基复合材料的高温疲劳损伤累积。在本发明中,所述疲劳迟滞耗散能方程优选如式5所示:
Figure GDA0002290954250000134
式5中,Un表示不同循环数疲劳迟滞耗散能,n表示循环数;
εunloading(σ)表示编织陶瓷基复合材料的卸载应变;
εreloading(σ)表示编织陶瓷剂复合材料的重新加载应变;
σmax表示疲劳峰值应力;
σmin表示疲劳谷值应力。
在本发明中,所述重新加载应变包括纤维/基体界面部分脱粘的重新加载应力对应的应变和纤维/基体界面完全脱粘的重新加载应力对应的应变,将部分脱粘和完全脱粘区域对应的应变参数分别代入式5,可得到准确度更高的疲劳迟滞耗散能;同理,卸载应变包括纤维/基体界面部分脱粘的卸载应力对应的应变和纤维/基体界面完全脱粘的卸载应力对应的应变,代入方式同上,不再重复说明。
在本发明中,通过所述疲劳迟滞耗散能方程可得到编织陶瓷基复合材料不同循环数下的疲劳迟滞耗散能,再根据疲劳迟滞耗散能确定编织陶瓷基界面脱粘和界面滑移的范围及长度,进而确定编织陶瓷基复合材料疲劳损伤积累的程度,以达到根据疲劳迟滞耗散能-循环数关系监测编织陶瓷基复合材料的高温疲劳损伤累积情况的目的。
得到疲劳迟滞耗散能方程后,本发明还优选根据所述的疲劳迟滞耗散能方程得到初始循环的疲劳迟滞耗散能和不同循环数的疲劳迟滞耗散能,结合编织陶瓷基复合材料的弹性应变能,得到编织陶瓷基复合材料的损伤参数方程,监测编织陶瓷基复合材料的高温疲劳损伤累积。
在本发明中,所述疲劳迟滞耗散能的损伤参数方程优选如式6所示:
Figure GDA0002290954250000141
式6中,Un表示不同循环数的疲劳迟滞耗散能;
Ue表示编织陶瓷基复合材料的弹性应变能;
Uinitial表示初始循环的疲劳迟滞耗散能;
Ф表示编织陶瓷基复合材料的损伤参数。
在本发明中,所述初始循环的疲劳迟滞耗散能优选为N=1时,根据疲劳迟滞耗散能方程得到的值。
在本发明中,所述弹性应变能优选通过疲劳峰值应力、疲劳谷值应力、疲劳峰值应力对应的应变和疲劳谷值应力对应的应变得到;更优选通过如式7所示的关系式得到:
Figure GDA0002290954250000142
式7中,Ue表示编织陶瓷基复合材料的弹性应变能;
σmax表示疲劳峰值应力;
σmin表示疲劳谷值应力;
εmax表示疲劳峰值应力对应的应变;
εmin表示疲劳谷值应力对应的应变。
本发明利用初始循环的疲劳迟滞耗散能、编织陶瓷基复合材料的弹性应变能,不同循环数对应的疲劳迟滞耗散能得到编织陶瓷基复合材料损伤参数,再根据损伤参数得到编织陶瓷基复合材料的疲劳损伤累积信息。在本发明中,所述损伤参数除包括疲劳迟滞耗散能外,还包括编织陶瓷基复合材料的弹性应变能,考虑的因素更接近材料实际使用的情况,因此获取的编织陶瓷基复合材料高温疲劳损伤积累信息也更为准确。
为了进一步说明本发明,下面结合附图和实施例对本发明提供的一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法进行详细地描述,但不能将它们理解为对本发明保护范围的限定。
实施例1
以2D SiC/SiC编织陶瓷基复合材料为测试样品,测试样品在高温600℃惰性气体环境,然后对测试样品的疲劳损伤进行监测:
提供参数:Ef=230GPa,Em=350GPa,Vf=40%,rf=7.5μm,αf=3.9×10-6/℃,ξd=1.0J/m2;T0=1000℃,χ=0.5,αm=2.0×10-6/℃,Vm=60%,σmax=130MPa,σmin=0MPa。
将上述参数代入建立的各项方程中,可得到编织陶瓷基复合材料不同循环数的疲劳迟滞耗散能,用于建立如图2所示的疲劳迟滞耗散能-循环数关系曲线。由图2可知,本发明提供方法所得疲劳迟滞耗散能的预测值与实验测试结果的吻合度极高,说明本发明提供的方法能够准确获取编织陶瓷基复合材料的疲劳迟滞耗散能,因此,根据该预测方法监测编织陶瓷基复合材料的疲劳损伤积累信息的可靠性较高。
尽管上述实施例对本发明做出了详尽的描述,但它仅仅是本发明一部分实施例,而不是全部实施例,人们还可以根据本实施例在不经创造性前提下获得其他实施例,这些实施例都属于本发明保护范围。

Claims (8)

1.一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法,包括以下步骤:
(1)根据断裂力学脱粘准则,利用与温度和循环数相关的纤维/基体界面脱粘区摩擦剪应力建立纤维/基体界面脱粘长度方程;
(2)根据断裂力学脱粘准则、纤维/基体界面滑移机理和所述步骤(1)得到的纤维/基体界面脱粘长度方程,建立卸载纤维/基体界面反向滑移长度方程;
(3)根据断裂力学脱粘准则、纤维/基体界面滑移机理、所述步骤(1)得到的纤维/基体界面脱粘长度方程和所述步骤(2)得到的卸载纤维/基体界面反向滑移长度方程,建立重新加载纤维/基体界面新滑移长度方程;
(4)根据所述步骤(1)得到的纤维/基体界面脱粘长度方程和所述步骤(2)得到的卸载纤维/基体界面反向滑移长度方程,结合总体载荷承担准则,建立卸载应力-应变方程;
根据所述步骤(1)得到的纤维/基体界面脱粘长度方程、所述步骤(2)得到的卸载纤维/基体界面反向滑移长度方程和所述步骤(3)得到的重新加载纤维/基体界面新滑移长度方程,结合总体载荷承担准则,建立重新加载应力-应变方程;
(5)根据所述步骤(4)得到的卸载应力-应变方程和重新加载应力-应变方程,建立疲劳迟滞耗散能方程,监测编织陶瓷基复合材料的高温疲劳损伤累积;
所述步骤(1)中纤维/基体界面脱粘长度方程如式1所示:
Figure FDA0002495485760000011
式1中,ld表示纤维/基体界面脱粘长度;
τi(T,N)表示与温度和循环数相关的纤维/基体界面脱粘区磨擦剪应力;
rf表示纤维半径;
Vm表示基体体积;
Em表示基体弹性模量;
σ表示应力;
χ表示沿应力加载方向纤维有效体积含量系数;
Vf表示编织陶瓷基复合材料中纤维体积含量;
Ec表示编织陶瓷基复合材料弹性模量;
ρ表示剪滞模型参数;
Ef表示纤维弹性模量;
ζd表示纤维/基体界面脱粘能。
2.如权利要求1所述的方法,其特征在于,所述步骤(1)中断裂力学脱粘准则如式1-1所示:
Figure FDA0002495485760000021
纤维轴向位移如式1-2所示:
Figure FDA0002495485760000022
所述纤维相对基体轴向位移如式1-3所示:
Figure FDA0002495485760000023
式1-1、1-2和1-3中,
F表示基体裂纹平面纤维承担载荷;
Figure FDA0002495485760000024
表示纤维轴向位移中的x为0时,对纤维/基体界面脱粘长度求导;
Figure FDA0002495485760000025
表示纤维相对基体轴向位移中的x为0时,对纤维/基体界面脱粘长度求导;
wf(x)表示纤维轴向位移;
v(x)表示纤维相对基体轴向位移;
lc表示基体裂纹间距。
3.如权利要求1所述的方法,其特征在于,所述卸载纤维/基体界面反向滑移长度方程如式2所示:
Figure FDA0002495485760000026
式2中,y表示卸载纤维/基体界面反向滑移长度;
σmax表示疲劳峰值应力。
4.如权利要求3所述的方法,其特征在于,所述重新加载纤维/基体界面新滑移长度方程如式3所示:
Figure FDA0002495485760000031
式3中,z表示重新加载纤维/基体界面新滑移长度;
Figure FDA0002495485760000032
σmin表示疲劳谷值应力。
5.如权利要求4所述的方法,其特征在于,纤维/基体界面部分脱粘时,所述卸载应力-应变方程如式4-1所示:
Figure FDA0002495485760000033
所述重新加载应力-应变方程如式4-2所示:
Figure FDA0002495485760000034
纤维/基体界面完全脱粘时,所述卸载应力-应变方程如式4-3所示:
Figure FDA0002495485760000035
所述重新加载应力-应变方程如式4-4所示:
Figure FDA0002495485760000036
式4-1、4-2、4-3和4-4中,
lc表示基体裂纹间距;
εunloading_p表示纤维/基体界面部分脱粘时卸载应力对应的应变;
εreloading_p表示纤维/基体界面部分脱粘时重新加载应力对应的应变;
εunloading_f表示纤维/基体界面完全脱粘时卸载应力对应的应变;
εreloading_f表示纤维/基体界面完全脱粘时重新加载应力对应的应变;
ac表示复合材料热膨胀系数;
af表示纤维热膨胀系数;
△T表示编织陶瓷基复合材料的测试温度与制备温度的差。
6.如权利要求1所述的方法,其特征在于,所述步骤(5)中疲劳迟滞耗散能方程如式5所示:
Figure FDA0002495485760000041
式5中,Un表示不同循环数的疲劳迟滞耗散能;
εunloading(σ)表示编织陶瓷基复合材料的卸载应变;
εreloading(σ)表示编织陶瓷剂复合材料的重新加载应变;
σmax表示疲劳峰值应力;
σmin表示疲劳谷值应力。
7.如权利要求1所述的方法,其特征在于,还包括根据所述步骤(5)得到的疲劳迟滞耗散能方程得到初始循环的疲劳迟滞耗散能和不同循环数的疲劳迟滞耗散能,结合编织陶瓷基复合材料的弹性应变能,得到编织陶瓷基复合材料的损伤参数方程,监测编织陶瓷基复合材料的高温疲劳损伤累积。
8.如权利要求7所述的方法,其特征在于,所述疲劳迟滞耗散能的损伤参数方程如式6所示:
Figure FDA0002495485760000042
式6中,Un表示不同循环数的疲劳迟滞耗散能;
Ue表示编织陶瓷基复合材料的弹性应变能;
Uinitial表示初始循环的疲劳迟滞耗散能;
Ф表示编织陶瓷基复合材料的损伤参数。
CN201811642683.5A 2018-12-29 2018-12-29 一种监测编织陶瓷基复合材料高温疲劳损伤累积的方法 Active CN109632887B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811642683.5A CN109632887B (zh) 2018-12-29 2018-12-29 一种监测编织陶瓷基复合材料高温疲劳损伤累积的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811642683.5A CN109632887B (zh) 2018-12-29 2018-12-29 一种监测编织陶瓷基复合材料高温疲劳损伤累积的方法

Publications (2)

Publication Number Publication Date
CN109632887A CN109632887A (zh) 2019-04-16
CN109632887B true CN109632887B (zh) 2020-08-14

Family

ID=66055113

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811642683.5A Active CN109632887B (zh) 2018-12-29 2018-12-29 一种监测编织陶瓷基复合材料高温疲劳损伤累积的方法

Country Status (1)

Country Link
CN (1) CN109632887B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110208086B (zh) * 2019-06-05 2021-07-09 东北石油大学 脉动压裂岩石累积损伤求解方法
CN111474062B (zh) * 2020-04-23 2021-05-14 南京航空航天大学 一种编织陶瓷基复合材料高温静疲劳损伤演化预测方法
CN111474044B (zh) * 2020-04-24 2021-03-19 南京航空航天大学 一种考虑多级疲劳损伤的编织陶瓷基复合材料疲劳迟滞回线的预测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09329531A (ja) * 1996-06-07 1997-12-22 Hitachi Ltd 高温構造部材の寿命診断素子
US8409491B1 (en) * 2007-09-28 2013-04-02 The United States of America as represented by the Administrator of the National Aeronautics & Space Administration (NASA) In-situ formation of reinforcement phases in ultra high temperature ceramic composites
CN104866690A (zh) * 2015-06-16 2015-08-26 南京航空航天大学 单向陶瓷基复合材料任意加卸载应力应变行为预测方法
CN105701312A (zh) * 2015-12-17 2016-06-22 南京航空航天大学 复杂编织结构陶瓷基复合材料疲劳迟滞行为预测方法
CN105760605A (zh) * 2015-12-17 2016-07-13 南京航空航天大学 复杂编织结构陶瓷基复合材料疲劳寿命预测方法
CN108195706A (zh) * 2017-12-25 2018-06-22 北京航空航天大学 一种陶瓷基复合材料结构件的热疲劳试验系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09329531A (ja) * 1996-06-07 1997-12-22 Hitachi Ltd 高温構造部材の寿命診断素子
US8409491B1 (en) * 2007-09-28 2013-04-02 The United States of America as represented by the Administrator of the National Aeronautics & Space Administration (NASA) In-situ formation of reinforcement phases in ultra high temperature ceramic composites
CN104866690A (zh) * 2015-06-16 2015-08-26 南京航空航天大学 单向陶瓷基复合材料任意加卸载应力应变行为预测方法
CN105701312A (zh) * 2015-12-17 2016-06-22 南京航空航天大学 复杂编织结构陶瓷基复合材料疲劳迟滞行为预测方法
CN105760605A (zh) * 2015-12-17 2016-07-13 南京航空航天大学 复杂编织结构陶瓷基复合材料疲劳寿命预测方法
CN108195706A (zh) * 2017-12-25 2018-06-22 北京航空航天大学 一种陶瓷基复合材料结构件的热疲劳试验系统

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
"Damage development and lifetime prediction of fiber-reinforced ceramic-matrix composites subjected to dwell-fatigue loading at elevated temperatures in oxidizing atmosphere";Li LONGBIAO;《Journal of the Ceramic Society of Japan》;20180731;第126卷(第7期);第516-528页 *
"Damage monitor and life prediction of carbon fiber-reinforced ceramic-matrix composites at room and elevated temperatures using hysteresis dissipated energy-based damage parameter";Li Longbiao;《Composite Interfaces》;20180212;第25卷(第4期);第335-356页 *
"MODELING IN-PHASE THERMOMECHANICAL FATIGUE HYSTERESIS LOOPS OF FIBER-REINFORCED CERAMIC-MATRIX COMPOSITES";LI LONGBIAO;《Ceramics-Silikáty》;20180830;第62卷(第3期);第240-247页 *
"Modeling out-of-phase thermomechanical fatigue hysteresis of long fiber-reinforced ceramic-matrix composites subjected to different loading sequence";Li LONGBIAO;《Journal of the Ceramic Society of Japan》;20180531;第126卷(第5期);第346-353页 *
"基于迟滞耗散能的纤维增强陶瓷基复合材料疲劳寿命预测方法";李龙彪;《复合材料学报》;20160430;第33卷(第4期);第841-851页 *
"纤维增强陶瓷基复合材料疲劳迟滞回线模型研究";李龙彪;《力学学报》;20140930;第46卷(第5期);第710-729页 *
"陶瓷基复合材料疲劳寿命预测方法";胡殿印 等;《航空发动机》;20170831;第43卷(第41期);第17-22页 *

Also Published As

Publication number Publication date
CN109632887A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN109614755B (zh) 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳纤维/基体界面剪应力的方法
CN109632530B (zh) 一种编织陶瓷基复合材料热机械疲劳迟滞回线的预测方法
CN109684598B (zh) 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳寿命的方法
CN109670272B (zh) 一种编织陶瓷基复合材料过载疲劳迟滞回线的预测方法
CN109781546B (zh) 一种编织陶瓷基复合材料拉伸强度的预测方法
CN109632887B (zh) 一种监测编织陶瓷基复合材料高温疲劳损伤累积的方法
CN111523237B (zh) 一种编织陶瓷基复合材料蠕变断裂行为的预测方法
CN111474044B (zh) 一种考虑多级疲劳损伤的编织陶瓷基复合材料疲劳迟滞回线的预测方法
CN109711100B (zh) 一种编织陶瓷基复合材料疲劳寿命的预测方法
CN111523232B (zh) 一种编织陶瓷基复合材料拉伸行为预测方法
CN111474043B (zh) 一种考虑多级疲劳损伤的编织陶瓷基复合材料剩余强度的预测方法
Longbiao A hysteresis dissipated energy-based parameter for damage monitoring of carbon fiber-reinforced ceramic–matrix composites under fatigue loading
WO2015053949A1 (en) Single crystal turbine blade lifing process and system
CN113032982B (zh) 一种考虑基体与纤维碎断的编织陶瓷基复合材料疲劳迟滞回线的预测方法
CN109740287B (zh) 一种编织陶瓷基复合材料模量衰退的预测方法
CN111474062B (zh) 一种编织陶瓷基复合材料高温静疲劳损伤演化预测方法
CN109598098A (zh) 一种考虑环境影响的编织陶瓷基复合材料稳态基体开裂应力的预测方法
CN109616163B (zh) 一种考虑环境影响的编织陶瓷基复合材料多基体开裂的预测方法
CN113032983B (zh) 一种考虑纤维断裂影响的编织陶瓷基复合材料振动阻尼的预测方法
CN114936474A (zh) 一种通过迟滞回线宽度预测纤维增强陶瓷基复合材料界面剪应力的方法
CN114943148A (zh) 一种通过切线模量预测纤维增强陶瓷基复合材料界面属性的方法
CN114139388A (zh) 一种预测纤维增强陶瓷基复合材料非闭合迟滞回线的方法
Abdi et al. Development of ASTM Test Standards for the Mode I Interlaminar Fracture Toughness (GIc--Crack Growth Resistance) of Ceramic Matrix Composites
CN113051719A (zh) 一种编织陶瓷基复合材料拉压疲劳迟滞回线的预测方法
Li Comparisons of the Thermomechanical Fatigue Behavior of C/SIC and SIC/SIC Ceramic-Matrix Composites Subjected to Different Phase Angle

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