CN105631148B - 应力氧化环境下单向陶瓷基复合材料力学性能分析方法 - Google Patents
应力氧化环境下单向陶瓷基复合材料力学性能分析方法 Download PDFInfo
- Publication number
- CN105631148B CN105631148B CN201511031307.9A CN201511031307A CN105631148B CN 105631148 B CN105631148 B CN 105631148B CN 201511031307 A CN201511031307 A CN 201511031307A CN 105631148 B CN105631148 B CN 105631148B
- Authority
- CN
- China
- Prior art keywords
- stress
- matrix
- ceramic matric
- matric composite
- unidirectional
- 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.)
- Expired - Fee Related
Links
- 238000007254 oxidation reaction Methods 0.000 title claims abstract description 92
- 239000002131 composite material Substances 0.000 title claims abstract description 90
- 230000003647 oxidation Effects 0.000 title claims abstract description 71
- 239000000919 ceramic Substances 0.000 title claims abstract description 68
- HPNSNYBUADCFDR-UHFFFAOYSA-N chromafenozide Chemical compound CC1=CC(C)=CC(C(=O)N(NC(=O)C=2C(=C3CCCOC3=CC=2)C)C(C)(C)C)=C1 HPNSNYBUADCFDR-UHFFFAOYSA-N 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 88
- 238000005336 cracking Methods 0.000 claims abstract description 48
- 230000008859 change Effects 0.000 claims abstract description 40
- 208000037656 Respiratory Sounds Diseases 0.000 claims abstract description 31
- 208000016261 weight loss Diseases 0.000 claims abstract description 29
- 230000004580 weight loss Effects 0.000 claims abstract description 29
- 230000007246 mechanism Effects 0.000 claims abstract description 5
- 230000035882 stress Effects 0.000 claims description 126
- 239000000835 fiber Substances 0.000 claims description 31
- 239000000463 material Substances 0.000 claims description 21
- 230000001590 oxidative effect Effects 0.000 claims description 20
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 13
- 229910052760 oxygen Inorganic materials 0.000 claims description 13
- 239000001301 oxygen Substances 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 12
- 229910052799 carbon Inorganic materials 0.000 claims description 8
- 239000000758 substrate Substances 0.000 claims description 8
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 7
- 239000007789 gas Substances 0.000 claims description 7
- 230000008646 thermal stress Effects 0.000 claims description 7
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- PNEYBMLMFCGWSK-UHFFFAOYSA-N aluminium oxide Inorganic materials [O-2].[O-2].[O-2].[Al+3].[Al+3] PNEYBMLMFCGWSK-UHFFFAOYSA-N 0.000 claims description 3
- 238000009833 condensation Methods 0.000 claims description 3
- 230000005494 condensation Effects 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 11
- 239000011204 carbon fibre-reinforced silicon carbide Substances 0.000 description 25
- 230000007423 decrease Effects 0.000 description 4
- 230000000630 rising effect Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000012856 packing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000035876 healing Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009954 braiding Methods 0.000 description 1
- 239000011153 ceramic matrix composite Substances 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (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
本发明公开了一种应力氧化环境下单向陶瓷基复合材料力学性能分析方法,包括以下步骤:确定单向陶瓷基复合材料基体的裂纹数随应力的变化;确定裂纹宽度变化规律;根据已有的400‑900℃无应力氧化环境下,单向陶瓷基复合材料的氧化机理,确定考虑应力作用下单向陶瓷基复合材料基体开裂对质量失重率λr变化的影响;确定考虑单向陶瓷基复合材料基体裂纹的情况下,单向陶瓷基复合材料剩余模量在应力氧化环境下的变化规律;根据单向陶瓷基复合材料体积分数变化规律,确定考虑基体裂纹情况下,单向陶瓷基复合材料剩余强度在应力氧化环境下的变化情况。本发明为UD‑CMC在应力氧化环境下力学性能的分析提供了相关理论支持。
Description
技术领域
本发明涉及一种分析单向陶瓷基复合材料(unidirectional ceramic matrixcomposite,简称UD-CMC)力学性能的理论方法,主要用于分析单向陶瓷基复合材料在应力氧化环境下材料剩余力学性能的预测。
背景技术
为实现航空发动机的高推比,提高涡轮前温度并降低发动机质量,CMC是航空发动机高温部件的理想选材。由于CMC服役环境十分恶劣,常常受到高温、高应力、热震及腐蚀性气体等的作用,大大降低了材料的性能。应力和高温是CMC服役环境最为典型的两个环境因素。考虑CMC在应力-氧化环境下力学性能的变化十分有必要。
在以往的研究分析中,对CMC在无应力氧化环境下力学性能变化的实验分析及理论分析并不多,对CMC在应力-氧化环境下力学性能变化的实验分析很少,而理论分析更是缺乏。单向陶瓷基复合材料结构简单,便于分析,为其他复杂编织陶瓷基复合材料提供理论参考和研究切入点。
CMC陶瓷基体相对纤维具有脆性,在应力的作用下将发生开裂,裂纹为氧气的进入提供通道。同时拉伸应力作用下,已经产生的裂纹宽度增加,增大了材料氧化的有效面积,促进材料的氧化。氧化使得材料性能降低,在拉伸作用下,材料的失效行为同常温环境下有很大的不同。应力和高温的相互作用,相互促进,形成一种新的陶瓷基复合材料动力学。
曾增,张茂庆,刘伟先等(单向C/SiC陶瓷基复合材料基体失效机制与强度预测.复合材料学报.2015,32(4):1075-1082)对单向C/SiC复合材料基体失效机制,并预测了材料的拉伸强度,但并未考虑氧化因素。孙志刚,牛序铭,王振剑等(Verification andprediction of residual strength of C/SiC composites under non-stressoxidation.Journal of Materials Science.2014,23(49):8192-8203)分析了氧化环境下单向陶瓷基复合材料力学性能的变化规律,但并未考虑拉伸应力的影响。
本发明考虑了基体在应力作用下形成裂纹,对UD-CMC在高温环境下的氧化力学性能变化规律进行分析和预测。
发明内容
为了对CMC在服役环境下力学性能更好的进行了解和预测,本发明的目的是提供一种应力氧化环境下单向陶瓷基复合材料力学性能分析方法,以解决现有技术中对CMC在高温下剩余力学性能变化的分析没有考虑应力的作用的问题。
为实现上述目的,本发明采用以下技术方案:
一种应力氧化环境下单向陶瓷基复合材料力学性能分析方法,包括以下步骤:
步骤一、确定单向陶瓷基复合材料基体的裂纹数随应力的变化:考虑热残余应力,并根据拉伸应力作用下基体承担的应力,推导得出单向陶瓷基复合材料基体裂纹数的变化情况;
步骤二、确定裂纹宽度变化规律:根据裂纹宽度与残余热应力及基体所承担的应力,推导得到应力作用下裂纹宽度的表达式;
步骤三、根据已有的400-900℃无应力氧化环境下,单向陶瓷基复合材料的氧化机理,结合步骤一和步骤二的结果,确定考虑应力作用下单向陶瓷基复合材料基体开裂对质量失重率λr变化的影响;
步骤四、确定考虑单向陶瓷基复合材料基体裂纹的情况下,单向陶瓷基复合材料剩余模量在应力氧化环境下的变化规律;
步骤五、根据步骤四得到的氧化后的单向陶瓷基复合材料体积分数变化规律,确定考虑基体裂纹情况下,单向陶瓷基复合材料剩余强度在应力氧化环境下的变化情况。
所述步骤一的具体步骤是:
采用蒙特卡罗法对基体在拉伸应力作用下的开裂进行分析,当材料受到轴向拉伸应力作用时,假设基体失效概率服从泊松分布,且应力作用下基体产生至少一条裂纹的概率等于1减去基体产生零条裂纹的概率,有:
P(ξ=σ;η=L0)=1-exp{-M(A)},N(A)≥1
式中,M(A)为泊松参数,N(A)为应力作用下产生的裂纹条数,P(ξ;η)表示特征长度为L0,应力为σ时,基体的失效概率;σR为特征应力,σth为热残余应力,σ*为基体的初始开裂应力;
根据蒙特卡罗法,采用计算机对裂纹数进行模拟,得到应力作用下相对裂纹数。
所述步骤二的具体步骤为:
假设裂纹均为贯穿裂纹,裂纹宽度与残余热应力σth及基体所承担的应力大小σm有关:
得到应力作用下裂纹宽度的表达式为:
式中,e为裂纹宽度,e0为初始裂纹宽度,αf、αm分别为纤维和基体的热膨胀系数,Ef、Em分别为纤维、基体的弹性模量,Vf、Vm分别为纤维、基体的体积分数;ΔT表示环境温度和基体开裂温度之间的温度差,基体开裂温度一般取为1030℃。
所述步骤三的具体步骤为:
已知在无应力氧化环境下,单向陶瓷基复合材料的质量失重率为:
400℃<T<700℃
700℃<T<900℃
Seff=nel
其中,Seff为氧气与纤维的有效接触面积,n为裂纹数量,e为裂纹宽度,l为裂纹深度;K0为一个常数,该常数与氧化反应速率相关,本发明中取为1.083×106;R为气体常数,Er为氧化反应所需活化能,W为陶瓷基复合材料的初始质量,T为氧化温度,P为大气压强,t为氧化时间,χO2为环境氧分压,Mc为碳的摩尔质量,Nc为碳的摩尔密度,Tc为基体开裂温度;n为裂纹数量,l为裂纹深度;Rg表示理想气体常数,M0表示空气的摩尔质量;
结合步骤二中得到的裂纹宽度变化规律,得到考虑基体开裂影响时,单向陶瓷基复合材料在应力氧化环境下的质量失重率的变化规律为:
400℃<T<700℃
700℃<T<900
℃
式中出现的参数同首次出现的意义相同。
所述步骤四的具体步骤是:
假设单向陶瓷基复合材料在氧化前后总体积保持不变,假设纤维均匀氧化,Vf和Vf’分别表示氧化前后纤维的体积分数,根据质量守恒有:
式中,ρc和ρf分别为单向陶瓷基复合材料的密度和纤维的密度;
在假设陶瓷基氧化前后体积不变的前提下,单向陶瓷基复合材料质量的减少主要是由于纤维体积分数减小引起的;因此,根据混合定律,氧化后单向陶瓷基复合材料的剩余模量Ec’为:
结合步骤三,将两个温度区间考虑基体裂纹情况下,质量失重率λr的表达式带入上式即可得到400-700℃温度区间和700-900℃温度区间,UD-CMC剩余强度变化:
400℃<T<700℃时:
700℃<T<900℃时:
步骤五的具体步骤是:
考虑多基体裂纹时,单向陶瓷基复合材料的剩余强度σuts的表达式为:
式中,m为weibull模量,σc为特征强度,此特征应力不同于基体失效概率中的特征应力,考虑应力氧化时有:
其中,σ0为无损情况下复合材料的拉伸强度,τi为界面剪切应力,L0为取样长度,rf为纤维半径;
因此,考虑基体裂纹情况下,单向陶瓷基复合材料剩余强度的变化规律为:
将步骤四得到的纤维体积分数带入上式即可得到考虑基体裂纹情况下,单向陶瓷基复合材料剩余强度在应力氧化环境下的变化规律:
400℃<T<700℃时:
700℃<T<900℃时:
本发明的有益效果是:与现有技术相比,本发明具有以下优点:
(1)本发明方法为UD-CMC在应力氧化环境下力学性能的分析提供了相关理论支持;(2)对UD-CMC在应力氧化环境下力学性能的变化有更深的了解,对UD-CMC的安全使用提供参考;(3)采用本发明方法预测的单向C/SiC复合材料在应力氧化环境下剩余模量和剩余强度值与实验数据吻合较好;(4)不考虑基体裂纹,对应力氧化环境下单向C/SiC复合材料力学性能的预测值和实验值之间的误差随着应力的增大而增大;而采用本发明方法预测得到的结果则稳定在某个范围之内,该结果表明本发明方法的合理性。
附图说明
图1为裂纹宽度应力、随氧化时间和氧化温度的变化曲线;
图2为400-700℃,单向C/SiC复合材料质量失重率随应力和时间的变化曲线;
图3为400-700℃,单向C/SiC复合材料质量失重率随应力和温度的变化曲线;
图4为700-900℃,单向C/SiC复合材料质量失重率随应力和时间的变化曲线;
图5为700-900℃,单向C/SiC复合材料质量失重率随应力和温度的变化曲线;
图6为400-700℃,单向C/SiC复合材料剩余模量随应力和时间的变化曲线;
图7为400-700℃,单向C/SiC复合材料剩余模量随应力和温度的变化曲线;
图8为700-900℃,单向C/SiC复合材料剩余模量随应力和时间的变化曲线;
图9为700-900℃,单向C/SiC复合材料剩余模量随应力和温度的变化曲线;
图10为400-700℃,单向C/SiC复合材料剩余强度随应力和时间的变化曲线;
图11为400-700℃,单向C/SiC复合材料剩余强度随应力和温度的变化曲线;
图12为700-900℃,单向C/SiC复合材料剩余强度随应力和时间的变化曲线;
图13为700-900℃,单向C/SiC复合材料剩余强度随应力和温度的变化曲线。
具体实施方式
下面结合附图对本发明作更进一步的说明。
研究发现,由于陶瓷基体具有一定的脆性,在应力作用下UD-CMC基体的失效过程为:1)在应力作用下基体中形成“裂纹源”;2)“裂纹源”发生扩展形成裂纹;3)新的“裂纹源”形成,同时裂纹继续扩展成为贯穿裂纹;4)裂纹数量不再增多,基体开裂达到饱和状态。
在高温环境下,应力的影响对UD-CMC力学性能的影响主要体现在增大了氧化的有效面积。主要有以下两个途径:1)应力增加,裂纹数增多;2)应力增加,裂纹宽度变大。因此,本发明对UD-CMC在应力氧化环境下力学性能的分析主要分为以下几个步骤:
步骤一、确定裂纹数随应力的变化
由于基体的开裂是一个随机过程,采用蒙特卡罗法对基体在拉伸应力作用下的开裂进行分析。实际的UD-CMC基体各处强度并不是均一的,而是具有一定分散性。当材料受到轴向拉伸应力作用时,假设基体失效概率服从泊松分布。且应力作用下基体产生至少一条裂纹的概率等于1减去基体产生零条裂纹的概率,有:
P(ξ=σ;η=L0)=1-exp{-M(A)},N(A)≥1 (1)
其中,M(A)为泊松参数,N(A)为应力作用下产生的裂纹条数,P(ξ;η)表示特征长度为L0,应力为σ时,基体的失效概率。σR为特征应力,σth为热残余应力,σ*为基体的初始开裂应力。
根据蒙特卡罗法,采用计算机对裂纹数进行模拟,可以得到应力作用下相对裂纹数。
根据步骤一给出的公式和方法得到一定应力作用下基体产生的裂纹数,如表1:
表1 应力作用下裂纹数的模拟结果
采用表2给出的相关参数值对单向C/SiC复合材料在应力氧化环境下的力学性能进行分析。
表2 相关参数
步骤二、确定裂纹宽度变化规律
假设裂纹均为贯穿裂纹,裂纹宽度与残余热应力σth及基体所承担的应力大小σm有关:
可以得到应力作用下裂纹宽度的表达式为:
其中,e为裂纹宽度,e0为初始裂纹宽度,αf、αm分别为纤维和基体的热膨胀系数,Ef、Em分别为纤维、基体的弹性模量,Vf、Vm分别为纤维、基体的体积分数;ΔT表示环境温度和基体开裂温度之间的温度差,基体开裂温度一般取为1030℃。
图1为根据步骤二中的公式(5)得到的裂纹宽度随拉伸应力和温度的变化曲线。其中曲线1-曲线5分别表示应力为200MPa、150MPa、100MPa、50MPa和0时裂纹宽度随氧化温度的变化情况,可以看到,基体的裂纹宽度与环境温度、应力大小以及材料各组分的热膨胀系数等因素有关。随温度的上升裂纹宽度逐渐降低;同一温度条件下,应力越大裂纹宽度越大。
步骤三、根据已有的400-900℃无应力氧化环境下,UD-CMC的氧化机理。结合步骤一和步骤二的结果,确定考虑应力作用下基体开裂对质量失重率λr变化的影响。
已知在无应力氧化环境下,UD-CMC的质量失重率为:
400℃<T<700℃ (6)
700℃<T<900℃ (7)
Seff=nel (8)
其中,Seff为氧气与纤维的有效接触面积,K0为一个常数,该常数与氧化反应速率相关,本发明中取1.083×106;R为气体常数,Er为氧化反应所需活化能,n为裂纹数量,l为裂纹深度,W为陶瓷基复合材料的初始质量,T为氧化温度,P为大气压强,t为氧化时间,χO2为环境氧分压,Mc为碳的摩尔质量,Nc为碳的摩尔密度,Tc为基体开裂温度,一般取1303K;Rg表示理想气体常数,M0表示空气的摩尔质量。
式(9)中的λ表示等号右边的关系式,是一个和环境相关的常数。
结合步骤二中得到的裂纹宽度变化规律,可以得到考虑基体开裂影响时,UD-CMC在应力氧化环境下的质量失重率的变化规律为。
400℃<T<700℃ (10)
700℃<T<900
℃ (11)
式中出现的参数同首次出现的意义相同。
图2-图5为根据步骤三中的公式(10)和公式(11)得到的400-700℃和700-900℃环境下,单向C/SiC复合材料质量失重率的变化曲线。其中,编号1-5分别表示应力大小为328MPa、240MPa、220MPa、186MPa和146MPa。图2为400-700℃质量失重率随氧化时间的变化曲线;图3为400-700℃质量失重率随氧化温度的变化曲线;图4为700-900℃质量失重率随氧化时间的变化曲线;图5为700-900℃质量失重率随氧化温度的变化曲线。可以看到,当复合材料的氧化由碳氧反应控制时(400-700℃),应力的作用对氧化具有促进作用,对同一氧化时间和氧化温度,应力越大材料的质量失重率越大;质量失重率随时间的增加呈线性增加,随温度的增加呈非线性增加,并且可以看到曲线的斜率逐渐增大,复合材料的质量失重率增加速度在逐渐增大。而700-900℃,质量失重率随氧化时间的增加不再呈线性变化,随氧化时间的增加材料失重增加,增加速度减缓,同时与400-700℃区间相比,氧化后失重率的变化量很小,应力越大质量失重率增加越多;而质量失重率随氧化温度的上升失重率略有降低,这是由于在该阶段氧气浓度的降低使得氧化程度降低,同时氧化产物对裂纹有一定的封填作用,应力越大,质量失重率越大,材料性能降低越多。
步骤四、确定考虑基体裂纹的情况下,UD-CMC剩余模量在应力氧化环境下的变化规律。
假设UD-CMC在氧化前后总体积保持不变,假设纤维均匀氧化,Vf和Vf’分别表示氧化前后纤维的体积分数,根据质量守恒有:
其中,ρc和ρf分别为UD-CMC的密度和纤维的密度。根据已有的研究知道,在400-900℃温度区间,陶瓷基体不发生氧化反应。在假设陶瓷基氧化前后体积不变的前提下,UD-CMC质量的减少主要是由于纤维体积分数减小引起的。因此,根据混合定律,氧化后UD-CMC的剩余模量Ec’为:
结合步骤三,将两个温度区间考虑基体裂纹情况下,质量失重率λr的表达式带入(13)式即可得到400-700℃温度区间和700-900℃温度区间,UD-CMC剩余强度变化。
400℃<T<700℃时:
700℃<T<900℃时:
图6-附图9为根据步骤四的公式(13)及步骤三的公式(10)和公式(11)得到的400-700℃和700-900℃环境下,单向C/SiC复合材料剩余模量的变化曲线。其中,编号1-5分别表示应力大小为146MPa、186MPa、220MPa、240MPa和328MPa。图6为400-700℃剩余模量随氧化时间的变化曲线;图7为400-700℃剩余模量随氧化温度的变化曲线;图8为700-900℃剩余模量随氧化时间的变化曲线;图9为700-900℃剩余模量随氧化温度的变化曲线。可以看到,温度为400-700℃时,剩余刚度随氧化时间的增加线性下降,应力的作用使得刚度下降较快;剩余刚度随氧化温度的增加呈非线性降低,同时拉伸应力的存在使得裂纹宽度增加,应力越大刚度下降越多。温度为700-900℃时,随着氧化时间的增加陶瓷基复合材料的剩余刚度缓慢降低,而应力的增加使得剩余刚度下降的幅度增大,这主要是由于应力的作用阻碍了该阶段裂纹的愈合;随着氧化温度的增加,陶瓷基复合材料的剩余刚度有略微的上升,这主要是由于该阶段裂纹受到氧化物的封填使得裂纹愈合造成的,应力越大,复合材料的剩余刚度越大,同时,复合材料剩余刚度的变化幅度越大。
步骤五、根据步骤四中得到的氧化后的UD-CMC体积分数变化规律(式(12)),确定考虑基体裂纹情况下,UD-CMC剩余强度在应力氧化环境下的变化情况。
考虑多基体裂纹时,UD-CMC的剩余强度σuts的表达式为:
式中,m为weibull模量,σc为特征强度,此特征应力不同于基体失效概率中的特征应力,考虑应力氧化时有:
其中,σ0为无损情况下复合材料的拉伸强度,τi为界面剪切应力,L0为取样长度,rf为纤维半径。
因此,考虑基体裂纹情况下,CMC剩余强度的变化规律为:
将步骤四得到的纤维体积分数(式(12))带入上式即可得到考虑基体裂纹情况下,UD-CMC剩余强度在应力氧化环境下的变化规律。
400℃<T<700℃时:
700℃<T<900℃时:
图10-附图13为根据步骤五的公式(16)及步骤四的公式(12)得到的400-700℃和700-900℃环境下,单向C/SiC复合材料剩余强度的变化曲线。其中,编号1-5分别表示应力大小为146MPa、186MPa、220MPa、240MPa和328MPa。图10为400-700℃剩余强度随氧化时间的变化曲线;图11为400-700℃剩余强度随氧化温度的变化曲线;图12为700-900℃剩余强度随氧化时间的变化曲线;图13为700-900℃剩余强度随氧化温度的变化曲线。可以看到,400-700℃温度区间,剩余强度随氧化时间、氧化温度及拉伸应力的变化规律与材料剩余刚度的变化规律基本一致。温度为700-900℃时,剩余强度随氧化时间的增加逐渐降低而随氧化温度的增加而有略微的上升,这主要是由于在该温度区间氧化过程由氧气的扩散控制,复合材料的氧化程度有所降低,同时之前的氧化产物对裂纹进行封填,复合材料的质量失重率下降所致。
表3和表4为采用本发明的方法预测的单向C/SiC复合材料的剩余模量(剩余刚度)和剩余强度与实验数据的对比结果。同时对比了不考虑基体裂纹情况下预测的单向C/SiC复合材料的剩余刚度和剩余强度。可以看到,不考虑基体开裂对陶瓷基复合材料剩余强度和剩余刚度的预测值与实验值相比其误差随着应力的增大而增大,而采用本发明的方法预测的结果随应力的增大误差始终稳定在8%-10%的范围内,表明考虑基体开裂的影响对分析单向陶瓷基复合材料在应力氧化条件下力学性能的变化是合理的。
表3 采用本发明预测的单向C/SiC复合材料剩余模量与实验和不考虑基体裂纹预测结果的对比结果
表4 采用本发明预测的单向C/SiC复合材料剩余强度与实验和不考虑基体裂纹预测结果的对比结果
虽然本发明已以较佳实施例公开如上,但他们并不是用来限定本发明的,任何熟习此技艺者,在不脱离本发明之精神和范围内,自当可做各种变化或润饰,因此本发明的保护范围应当以本申请的专利保护范围所界定的为准。
Claims (1)
1.一种应力氧化环境下单向陶瓷基复合材料力学性能分析方法,其特征在于:包括以下步骤:
步骤一、确定单向陶瓷基复合材料基体的裂纹数随应力的变化:考虑残余热应力,并根据拉伸应力作用下基体承担的应力,推导得出单向陶瓷基复合材料基体裂纹数的变化情况;具体步骤是:
采用蒙特卡罗法对基体在拉伸应力作用下的开裂进行分析,当材料受到轴向拉伸应力作用时,假设基体失效概率服从泊松分布,且应力作用下基体产生至少一条裂纹的概率等于1减去基体产生零条裂纹的概率,有:
C(ξ=σ;η=L0)=1-exp{-M(A)},N(A)≥1
式中,M(A)为泊松参数,N(A)为应力作用下产生的裂纹条数,C(ξ;η)表示特征长度为L0,应力为σ时,基体的失效概率;σR为特征应力,σth为残余热应力,σ*为基体的初始开裂应力;
根据蒙特卡罗法,采用计算机对裂纹数进行模拟,得到应力作用下相对裂纹数;
步骤二、确定裂纹宽度变化规律:根据裂纹宽度与残余热应力及基体所承担的应力,推导得到应力作用下裂纹宽度的表达式;具体步骤为:
假设裂纹均为贯穿裂纹,裂纹宽度与残余热应力σth及基体所承担的应力大小σm有关:
得到应力作用下裂纹宽度的表达式为:
式中,e为裂纹宽度,e0为初始裂纹宽度,αf、αm分别为纤维和基体的热膨胀系数,Ef、Em分别为纤维、基体的弹性模量,Vf、Vm分别为纤维、基体的体积分数,ΔT表示环境温度和基体开裂温度之间的温度差;
步骤三、根据已有的400-900℃无应力氧化环境下,单向陶瓷基复合材料的氧化机理,结合步骤一和步骤二的结果,确定考虑应力作用下单向陶瓷基复合材料基体开裂对质量失重率λr变化的影响;具体步骤为:
已知在无应力氧化环境下,单向陶瓷基复合材料的质量失重率为:
Seff=nel
令:
其中,Seff为氧气与纤维的有效接触面积;K0为一个常数,该常数与氧化反应速率相关,R为气体常数,Er为氧化反应所需活化能,W为陶瓷基复合材料的初始质量,T为氧化温度,P为大气压强,t为氧化时间,χO2为环境氧分压,Mc为碳的摩尔质量,Nc为碳的摩尔密度,Tc为基体开裂温度;n为裂纹数量,l为裂纹深度;Rg表示理想气体常数,M0表示空气的摩尔质量;
结合步骤二中得到的裂纹宽度变化规律,得到考虑基体开裂影响时,单向陶瓷基复合材料在应力氧化环境下的质量失重率的变化规律为:
式中出现的参数同首次出现的意义相同;
步骤四、确定考虑单向陶瓷基复合材料基体裂纹的情况下,单向陶瓷基复合材料剩余模量在应力氧化环境下的变化规律;具体步骤是:
假设单向陶瓷基复合材料在氧化前后总体积保持不变,假设纤维均匀氧化,Vf和Vf’分别表示氧化前后纤维的体积分数,根据质量守恒有:
式中,ρc和ρf分别为单向陶瓷基复合材料的密度和纤维的密度;
在假设陶瓷基氧化前后体积不变的前提下,单向陶瓷基复合材料质量的减少主要是由于纤维体积分数减小引起的;因此,根据混合定律,氧化后单向陶瓷基复合材料的剩余模量Ec’为:
结合步骤三,将两个温度区间考虑基体裂纹情况下,质量失重率λr的表达式带入上式即可得到400-700℃温度区间和700-900℃温度区间,单向陶瓷基复合材料剩余模量变化:
400℃<T<700℃时:
700℃<T<900℃时:
步骤五、根据步骤四得到的氧化后的单向陶瓷基复合材料体积分数变化规律,确定考虑基体裂纹情况下,单向陶瓷基复合材料剩余强度在应力氧化环境下的变化情况,具体步骤是:
考虑多基体裂纹时,单向陶瓷基复合材料的剩余强度σuts的表达式为:
式中,m为weibull模量,σc为特征强度,此特征应力不同于基体失效概率中的特征应力,考虑应力氧化时有:
其中,σ0为无损情况下复合材料的拉伸强度,τi为界面剪切应力,L0为取样长度,rf为纤维半径;
因此,考虑基体裂纹情况下,单向陶瓷基复合材料剩余强度的变化规律为:
将步骤四得到的纤维体积分数带入上式即可得到考虑基体裂纹情况下,单向陶瓷基复合材料剩余强度在应力氧化环境下的变化规律:
400℃<T<700℃时:
700℃<T<900℃时:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511031307.9A CN105631148B (zh) | 2015-12-31 | 2015-12-31 | 应力氧化环境下单向陶瓷基复合材料力学性能分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511031307.9A CN105631148B (zh) | 2015-12-31 | 2015-12-31 | 应力氧化环境下单向陶瓷基复合材料力学性能分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105631148A CN105631148A (zh) | 2016-06-01 |
CN105631148B true CN105631148B (zh) | 2019-05-14 |
Family
ID=56046077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511031307.9A Expired - Fee Related CN105631148B (zh) | 2015-12-31 | 2015-12-31 | 应力氧化环境下单向陶瓷基复合材料力学性能分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105631148B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108845035A (zh) * | 2018-07-06 | 2018-11-20 | 北京领示科技有限公司 | 一种基于超声无损定量检出概率的延寿方法 |
CN109583052B (zh) * | 2018-11-15 | 2022-04-05 | 西南交通大学 | 纤维增强陶瓷基复合材料微结构精确建模及模量计算方法 |
CN109781546B (zh) * | 2018-12-29 | 2020-09-22 | 南京航空航天大学 | 一种编织陶瓷基复合材料拉伸强度的预测方法 |
CN109884105B (zh) * | 2019-01-07 | 2020-03-10 | 南京航空航天大学 | 一种陶瓷基复合材料氧化环境碳界面消耗体积确定方法 |
CN109740287B (zh) * | 2019-01-23 | 2020-09-22 | 南京航空航天大学 | 一种编织陶瓷基复合材料模量衰退的预测方法 |
CN110096732B (zh) * | 2019-03-15 | 2020-03-06 | 南京航空航天大学 | 一种陶瓷基复合材料在应力氧化环境下剩余刚度预测方法 |
CN110096731B (zh) * | 2019-03-15 | 2020-03-06 | 南京航空航天大学 | 一种陶瓷基复合材料在应力氧化环境下质量变化预测方法 |
CN109992850B (zh) * | 2019-03-15 | 2020-01-24 | 南京航空航天大学 | 一种陶瓷基复合材料应力氧化环境剩余拉伸强度预测方法 |
CN109815643B (zh) * | 2019-04-04 | 2020-04-14 | 南京航空航天大学 | 单向陶瓷基复合材料任意应变加卸载本构关系预测方法 |
CN110146550B (zh) * | 2019-06-13 | 2020-03-13 | 南京航空航天大学 | 基于电阻抗成像的复合材料高温部件氧化程度监测方法 |
CN110246548B (zh) * | 2019-06-17 | 2020-04-07 | 南京航空航天大学 | 一种应力水蒸气环境下单向陶瓷基复合材料内部氧化形貌预测方法 |
CN110309619B (zh) * | 2019-07-09 | 2020-05-26 | 南京航空航天大学 | 一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法 |
CN110362956B (zh) * | 2019-07-25 | 2021-02-19 | 南京航空航天大学 | 一种陶瓷基复合材料在高温应力环境下剩余刚度计算方法 |
CN111826603B (zh) * | 2020-07-21 | 2022-09-30 | 苏州微创关节医疗科技有限公司 | 制备金属基底表面氧化陶瓷层的方法及应用 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6468915B1 (en) * | 2000-09-21 | 2002-10-22 | Taiwan Semiconductor Manufacturing Company | Method of silicon oxynitride ARC removal after gate etching |
CN102339348A (zh) * | 2011-08-11 | 2012-02-01 | 西北工业大学 | 氧化环境中单向C/SiC复合材料细观应力的检测方法 |
CN103871059A (zh) * | 2014-03-13 | 2014-06-18 | 南京航空航天大学 | 纤维增强复合材料等效弹性参数的计算方法 |
CN104866690A (zh) * | 2015-06-16 | 2015-08-26 | 南京航空航天大学 | 单向陶瓷基复合材料任意加卸载应力应变行为预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100454715B1 (ko) * | 2002-03-14 | 2004-11-05 | 한국과학기술연구원 | MoSi₂―Si₃N₄복합피복층 및 그 제조방법 |
-
2015
- 2015-12-31 CN CN201511031307.9A patent/CN105631148B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6468915B1 (en) * | 2000-09-21 | 2002-10-22 | Taiwan Semiconductor Manufacturing Company | Method of silicon oxynitride ARC removal after gate etching |
CN102339348A (zh) * | 2011-08-11 | 2012-02-01 | 西北工业大学 | 氧化环境中单向C/SiC复合材料细观应力的检测方法 |
CN103871059A (zh) * | 2014-03-13 | 2014-06-18 | 南京航空航天大学 | 纤维增强复合材料等效弹性参数的计算方法 |
CN104866690A (zh) * | 2015-06-16 | 2015-08-26 | 南京航空航天大学 | 单向陶瓷基复合材料任意加卸载应力应变行为预测方法 |
Non-Patent Citations (4)
Title |
---|
C/SiC复合材料应力氧化失效机理;刘小瀛 等;《无机材料学报》;20060930;第21卷(第5期);全文 |
Verification and prediction of residual strength of C/SiC composites under non一stress oxidation;Zhigang Sun et al.;《Journal of Materials Science》;20141130;第49卷(第23期);全文 |
氧化环境下陶瓷基符合材料力学性能研究;王振剑;《中国优秀硕士论文全文数据库工程科技I辑》;20120715(第7期);论文第1-37页和摘要 |
高温氧化气氛下3D C/SiC质量变化率与剩余强度的相关性;李建章 等;《复合材料学报》;20070831;第24卷(第4期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN105631148A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105631148B (zh) | 应力氧化环境下单向陶瓷基复合材料力学性能分析方法 | |
Zok et al. | Degradation of a SiC‐SiC composite in water vapor environments | |
Gauthier et al. | Delayed failure of Hi‐Nicalon and Hi‐Nicalon S multifilament tows and single filaments at intermediate temperatures (500–800 C) | |
Gauthier et al. | Oxidation of silicon carbide fibers during static fatigue in air at intermediate temperatures | |
Patel et al. | Effect of oxidation and residual stress on mechanical properties of SiC seal coated C/SiC composite | |
CN111400906B (zh) | 应力氧化环境下单向陶瓷基复材应力应变曲线预测方法 | |
Mehrman et al. | Influence of hold times on the elevated-temperature fatigue behavior of an oxide–oxide ceramic composite in air and in steam environment | |
Chawla et al. | Monotonic and Cyclic Fatigue Behavior of High‐Performance Ceramic Fibers | |
Ruggles-Wrenn et al. | Effects of steam environment on creep behavior of Nextel™ 720/alumina ceramic composite at elevated temperature | |
Almeida et al. | Thermal exposure effects on the long‐term behavior of a mullite fiber at high temperature | |
Li et al. | Estimate interface frictional coefficient of ceramic matrix composites from hysteresis loops | |
Laforêt et al. | Static fatigue of multifilament tows at high temperatures above 900 C | |
Parthasarathy et al. | Predicted effects of interfacial roughness on the behavior of selected ceramic composites | |
Rugg et al. | Creep of SiC–SiC microcomposites | |
Białas et al. | Numerical simulation of segmentation cracking in thermal barrier coatings by means of cohesive zone elements | |
Longbiao | Time-dependent proportional limit stress of carbon fiber-reinforced silicon carbide ceramic-matrix composites considering interface oxidation | |
Song et al. | Long term stability and mechanical properties of Al2O3–NiAl composites reinforced with partially fragmented long fibers | |
Sun et al. | Verification and prediction of residual strength of C/SiC composites under non-stress oxidation | |
Longbiao | Modeling the effect of oxidation on fatigue life of carbon fiber-reinforced ceramic-matrix composites at elevated temperature | |
Zhang et al. | Comparative Analysis of Low‐Cycle Fatigue Behavior of 2D‐Cf‐PyC/SiC Composites in Different Environments | |
Ochiai et al. | Degradation of SiC/SiC composite due to exposure at high temperatures in vacuum in comparison with that in air | |
Ding et al. | A Progressive Oxidative Damage Model of C/SiC Composites under Stressed Oxidation Environments | |
Li | Effect of Interface Properties on Tensile and Fatigue Behavior of 2D Woven SiC/SiC Fiber‐Reinforced Ceramic‐Matrix Composites | |
Li | Effect of multiple loading sequence on time‐dependent stress rupture of fiber‐reinforced ceramic‐matrix composites | |
Jones et al. | Crack bridging by sic fibers during slow crack growth and the resultant fracture toughness of sicsicf composites |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190514 |
|
CF01 | Termination of patent right due to non-payment of annual fee |