High temperature mechanical behavior and damage mechanism of 2.5D woven carbon fiber-glass fiber/bismaleimide resin composite
-
摘要: 采用三维纺织技术与树脂传递模塑工艺(RTM)制备了2.5D机织碳纤维-玻璃纤维/双马来酰亚胺树脂复合材料,分别在室温(25℃)和高温(150℃、240℃、300℃)环境下对复合材料进行了三点弯曲和层间剪切力学性能测试,探究了温度对复合材料力学行为及损伤机制的影响。结果表明:温度对2.5D机织碳纤维-玻璃纤维/双马来酰亚胺树脂复合材料的力学特性和损伤模式均有显著影响。温度上升导致纤维/基体界面结合力减弱,复合材料在300℃下的弯曲强度、弯曲模量和层间剪切强度相比室温环境分别降低了23.06%、70.01%和18.93%。弯曲载荷下,2.5D机织混杂复合材料的室温破坏模式主要以局部的纤维断裂和基体开裂为主,而高温损伤则以纤维/基体界面脱粘为主导。剪切载荷下,2.5D机织混杂复合材料的室温破坏模式主要为分层破坏,而随着温度升高,复合材料因基体软化出现塑性变形,基体开裂、界面脱粘及分层破坏决定了材料的最终失效。
-
关键词:
- 2.5D机织复合材料 /
- 混杂纤维 /
- 温度 /
- 力学行为 /
- 损伤机制
Abstract: 2.5D woven carbon fiber-glass fiber/bismaleimide resin composite was prepared by three-dimensional textile technology and resin transfer molding (RTM), and the mechanical properties of three-point bending and interlaminar shear were tested at room temperature (25℃) and high temperature (150℃, 240℃, 300℃), respectively, and the influence of temperature on the mechanical behavior and damage mechanism of the composites was explored. The results show that temperature has a significant effect on the mechanical properties and damage mode of 2.5D woven carbon fiber-glass fiber/bismaleimide resin composites. Temperature rise leads to weakening of fiber/matrix interface adhesion. Compared with the room temperature environment, the flexural strength, flexural modulus and interlaminar shear strength of the composite at 300℃ decrease by 23.06%, 70.01% and 18.93%. Under bending load, the failure mode of 2.5D woven hybrid composites at room temperature is mainly dominated by local fiber fracture and matrix cracking, while high temperature damage is dominated by fiber/matrix interface debonding. Under shear load, the failure mode of 2.5D woven hybrid composites at room temperature is mainly stratified failure. With the increase of temperature, the composites appear plastic deformation due to the softening of matrix, matrix cracking, interface debonding and stratified failure determine the final failure of the material. -
连续碳化硅纤维增韧碳化硅陶瓷基复合材料(SiC/SiC)具备高强度、耐高温、低密度等特点,可作为新型航空发动机的优选材料之一[1-5]。随着SiC/SiC复合材料应用范围愈加广泛,不仅应用于航空领域,在化工领域中也是制造密封环的理想材料,在核聚变反应堆中可用于包层结构、包层流道内衬和转换器等方面[6],这些尖端领域均对复合材料的稳定性提出要求。化学气相渗透法(CVI)工艺导致陶瓷基复合材料(CMC)内部普遍存在一些缺陷(束内外孔隙、分层、预制体铺层错位和基体开裂等),故其力学性能具有较大的离散性,强度作为结构设计的一项重要参数势必会受到其离散性的影响,从而限制结构的高可靠性设计。
CMC成型编织和致密化工艺复杂,易形成损伤源,有必要综合考虑制备工艺中的损伤对复合材料力学性能的影响,建立概率论可靠性模型[7]。大量研究表明CMC拉伸强度服从Weibull分布[8-13]。陶永强等[14]通过渐进损伤的简化剪滞理论得到了各层应力分布及90°层裂纹密度与施加应力之间的关系;Curtin等[15]基于基体裂纹随机开裂、纤维随机破坏理论,对单向增强CMC的应力-应变行为进行了模拟;孟志新等[16]探究了不同纤维束下SiC CMC的拉伸行为;Vagaggini等[17]计算了界面性能参数对复合材料拉伸行为的影响,表明增大界面脱粘能和界面滑移应力利于提高拉伸强度。目前CMC已实现工程化应用,评估工艺放大对复合材料性能分散性的影响成为主要的技术问题,分析单向纤维束复合材料拉伸强度的可靠性,评估结构单元和微结构特征对分散性的影响,也是该材料由小试走向中试的重要环节,小试指实验室开发优化阶段,中试指在小试的基础上继续放大、量产化。关于各组元对CMC性能影响的研究已有很多报道,但是关于小试炉与中试炉对CMC力学性能及其各组元影响的研究鲜有报道。
本文采用CVI技术分别在小试炉与中试炉中制备了单向纤维束SiC/SiC复合材料(Mini-SiC/SiC)的两批试样,探究工艺放大制备对复合材料力学性能离散的影响,并讨论造成复合材料力学性能离散的来源。其中小试炉为实验室制备,一台CVI沉积炉,中试为90台沉积炉,小试炉膛直径为600 mm,中试炉膛直径为1200 mm。通过双参数Weibull分布模型对Mini-SiC/SiC复合材料强度分布和可靠性进行分析,并基于双参数Weibull分布损伤模型建立拉伸载荷-位移关系,最终采用ORS Dragonfly软件分割各组元,结合基体间距裂纹公式分析分散性对可靠性的影响。
1. 材料制备与测试分析
1.1 Mini-SiC/SiC复合材料的制备
采用CVI法制备Mini-SiC/SiC复合材料,见图1。首先,均匀地将SiC纤维束缠绕到石墨框上,获得纤维预制体。然后,在纤维预制体表面沉积BN界面层,其中以BCl3作为硼源,NH3作为氨源,Ar为稀释气体,沉积温度600℃[18]。最后,将沉积BN界面层的SiC纤维预制体在化学气相沉积炉中浸渗SiC基体,选择三氯甲基硅烷(MTS)为SiC基体的源物质,以H2为载气,Ar为稀释气体,MTS与H2摩尔比1∶10,沉积温度1150℃。Mini-SiC/SiC复合材料制备工艺参数如表1所示,Mini-SiC/SiC A为小试炉制备(炉膛直径为600 mm),Mini-SiC/SiC B为中试炉制备(炉膛直径为1200 mm),小试炉制备复合材料平均密度为3.09 g/cm2,中试炉制备复合材料平均密度为3.07 g/cm2。SiC纤维耐温为1350℃,故制备温度对其影响不大。
表 1 Mini-SiC/SiC复合材料的制备工艺参数Table 1. Process parameters of Mini-SiC/SiC compositesMaterial name Deposition time/h Inside diameter/mm BN interface SiC matrix Mini-SiC/SiC A 100 160 600 Mini-SiC/SiC B 100 160 1200 Notes: Mini-SiC/SiC A—Preparation in a small test furnace (Furnace diameter of 600 mm); Mini-SiC/SiC B—Preparation in a pilot test furnace (Furnace diameter of 1200 mm). 1.2 力学测试与分析方法
拉伸测试在美特斯(MTS)公司生产的CMT6203型微机控制电子万能试验机上进行。实验执行标准为GB/T 1040.5—2008[19],标距50 mm,加载速度3 mm/min。其拉伸测试特定夹具组件见图2。
Mini-SiC/SiC复合材料性能分析,采用以概率论为根据的结构可靠性模型,基于双参数Weibull分布模型对Mini-SiC/SiC复合材料的强度分布和可靠性进行分析。
1.3 表征与分析
采用扫描电子显微镜(SEM,SU3800,日本日立)观察Mini-SiC/SiC复合材料的显微结构和微观形貌。ORS Dragonfly软件在阈值分割的基础上,对Mini-SiC/SiC复合材料断层扫描(CT,AX-2000,宁波奥影检测)图像使用深度学习工具分割出灰度值相近的部分进行处理[20]。
2. 结果与讨论
2.1 Mini-SiC/SiC复合材料非线性拉伸行为
图3为典型Mini-SiC/SiC复合材料载荷-位移曲线。该图由两组曲线组成,图3(a)、3(b)分别为小试炉与中试炉制备,可见曲线规律一致,故各选取其中一条曲线分析,曲线分三阶段:第一阶段加载初期,复合材料整体发生弹性变形,拉伸曲线斜率保持不变;第二阶段曲线呈非线性,且斜率逐渐减小,SiC基体产生新裂纹,随应力增加而增多;当达到最大载荷时曲线进入第三阶段,该阶段曲线未表现出韧性现象,呈现当基体裂纹达到饱和时,直接发生脆性断裂,其断口主要以纤维单丝形式拔出,未出现纤维簇形式拔出。
2.2 Mini-SiC/SiC复合材料拉伸强度统计分布规律
采用中位估计法获取材料强度分布,计算每个强度所对应的失效概率,该失效概率为试验估计值。中位估计法失效概率估计值Fn可表示为
Fn(σ)=i−0.5N (1) 式中:i为强度从小到大排列的序数,即强度σ所对应的序号,i的最大值为 N;N为样本总量。
双参数Weibull分布模型累积失效分布函数F(σ) 数学表达式为
F(σ)=1−exp[−(σσ0)m] (2) 式中:
σ0 为尺度参数,当 σ0=σ时,Weibull函数取最大值;m为形状参数,表示材料强度的分布规则,m越大,强度分散性越小。将式(2)整理得:
ln{ln[11−F(σ)]}=mlnσ−mlnσ0 (3) 由式(3)将
ln{ln[1/(1−F(σ))]} 作为y轴,lnσ作为x轴,可得:y=mx−mlnσ0 (4) 将式(4)进行作图,并线性拟合分析,由此可得m和σ0值。
2.2.1 参数估计方法
Weibull分布参数求解及检验主要数据(见附录表4),线性拟合如图4所示。以Mini-SiC/SiC A复合材料为例,进行双参数Weibull分布模型的参数估计。
由图4(a)可得,Mini-SiC/SiC A线性拟合曲线函数为
y=20.59x−122.03 (5) 由此可知,m=20.59,σ0=374.79 MPa。同理可得,Mini-SiC/SiC B的m=5.01,σ0=400.74 MPa,前者m值大于后者,表明小试炉制备稳定性更优。此外,Weibull分布参数和线性相关系数如表2所示。由表可得,线性相关系数r都接近于1,变量之间存在高度相关性,线性拟合合理,因此m与σ0接近真实值。
表 2 两种Mini-SiC/SiC复合材料拉伸强度双参数Weibull分布的参数和线性相关系数Table 2. Parameters and linear correlation coefficient of two-parameter Weibull distribution for tensile strength of two Mini-SiC/SiC compositesName Shape parameter b Scale parameter a/mm r Mini-SiC/SiC A 20.59 374.79 0.94 Mini-SiC/SiC B 5.01 400.74 0.98 将上述得到的m和σ0代入式(2),Mini-SiC/SiC A复合材料的双参数Weibull累积失效分布函数可表示为
F(σ)=1−exp[−(σ374.79)20.59] (6) 同理可得,Mini-SiC/SiC B的双参数Weibull累积失效分布函数可表示为
F(σ)=1−exp[−(σ400.74)5.01] (7) 图5为两种Mini-SiC/SiC复合材料双参数Weibull累积失效分布函数曲线与中位估计数据曲线对比图。分析可知,分别采用式(6)、式(7)计算的曲线与中位估计法得到的曲线能够较好吻合。因此,拉伸强度分布规律能够采用双参数Weibull分布模型表征,且效果较好。
为检验Mini-SiC/SiC A复合材料拉伸强度与双参数Weibull分布的拟合优度。本文引入Kolmogorov-Smirnov(K-S)检验,可通过样本观测值的失效概率估计值Fn(σ)与假设的双参数Weibull累积失效概率分布F(σ)的差值进行对比分析,从而确定该失效概率的分布类型。
Fn(σ)与F(σ)的最大差值DN由下式得出:
DN=max|Fn(σ)−F(σ)| (8) Fn(σ)和F(σ)在上文中已完成计算,前者由中位估计法获得,后者则通过双参数Weibull分布函数计算得到,都代表材料在强度达到σ时发生失效的概率。DN数据见附录表4。
工程上一般采取显著水平α=0.05或α=0.01。经查阅DN临界值表,当N大于35,DN临界值等于D35与
√N 的比值,如附录表1所示。当复合材料Mini-SiC/SiC A拉伸强度为356.48 MPa时(附录表4),Kolmogorov距离取最大值DN=0.1252<D60,0.05=0.1756<D60,0.01=0.2104,可知DN小于显著水平α=0.05和α=0.01的DN临界标准值。
采用相同方法可检验Mini-SiC/SiC B拉伸强度是否服从双参数Weibull分布,检验结果数据如附录表2所示。由表可知,DN分别为0.0995均小于D60,0.01和D60,0.05,故拉伸强度均服从双参数Weibull分布,且拟合优度良好。
根据K-S检验可知,DN小于显著水平α=0.05和α=0.01的DN临界标准值。故拉伸强度均服从双参数Weibull分布,且拟合优度良好。
可靠度R是度量可靠性的标准,其与失效概率关系可表示为
R(σ)=1−F(σ) (9) 该式意味着给定任意强度,都可计算出其可靠度为多大,并随着强度增大,失效概率会增大,可靠度会减小。由于给出任意强度数据便可计算其可靠度,因此对于每个给定的可靠度,也可以计算出相应的可靠拉伸强度。联立式(2)和式(9)并求反函数,可获得给定可靠度下的可靠拉伸强度σ如下式:
σ(R)=σ0(ln1R)1m (10) 由此就可获得任意可靠度所对应的可靠强度值。在工程上,一般将5%、25%、50%、75%、95%可靠度作为可靠质量标准。5个可靠度下的可靠拉伸强度值如附录表3所示。同时,不同可靠度下的可靠拉伸强度如附录图1所示。
研究表明,可靠度为50%时的可靠拉伸强度与平均拉伸强度较为接近,表3给出了可靠度为50%的可靠拉伸强度与平均拉伸强度。由表中数据通过计算绝对偏差和相对偏差可知:最大绝对偏差仅4.69 MPa,最大相对偏差仅1.2%。因此,采用平均拉伸强度可准确预测50%可靠度所对应的可靠强度,且误差较小。
表 3 两种Mini-SiC/SiC复合材料在可靠度50%时的可靠拉伸强度与平均拉伸强度Table 3. Reliable tensile strength and average tensile strength of two Mini-SiC/SiC composites at 50% reliabilityName Mean tensile strength/MPa Tensile strength of 50% reliability/MPa Absolute deviation/MPa Relative deviation/% Mini-SiC/SiC A 365.08 368.18 3.10 0.8 Mini-SiC/SiC B 367.78 372.47 4.69 1.2 2.3 Mini-SiC/SiC复合材料拉伸载荷与位移关系的确定
复合材料拉伸应力σ与载荷P关系为
σ=PS (11) 复合材料拉伸应变ε与伸长量∆l关系为
ε=Δll0 (12) 式中:S为材料截面积;
l0 为标距,值为50 mm。材料宏观损伤变量服从双参数Weibull分布,由式(2)与拉伸强度建立联系。在拉伸加载时,损伤贯穿整个过程,故累积损伤量和位移x同样服从双参数Weibull分布。因此累积损伤因子Da随着位移x变化的分布函数可表示为类似于式(2)形式,如下所示:
Da=1−exp[−(xa)b] (13) 由于SiC/SiC复合材料为脆性材料,为便于计算参数a、b,获取Weibull分布模型,因此将x视为断裂位移xmax。
将式(13)整理得:
ln{ln[11−Da]}=blnxmax−blna (14) 可将
ln{ln[1/(1−Da)]} 作为y轴,lnxmax作为x轴,得:y=bxmax−blna (15) SiC/SiC复合材料中,SiC纤维和SiC基体均为脆性材料,各组分材料达到破坏极限前几乎未发生塑性变形。因此,基于双参数Weibull分布的弹性损伤模型表征Mini-SiC/SiC复合材料拉伸加载过程:
P=k(1−Da)x (16) 式中,k为抗拉刚度。
2.3.1 双参数Weibull分布模型参数估计
Weibull分布参数求解及检验的详细数据见附录表5,线性拟合如图6所示。以Mini-SiC/SiC A复合材料为例,进行双参数Weibull分布模型参数估计。
Mini-SiC/SiC A复合材料拟合曲线函数为
y=3.09x+0.42 (17) 由此可得,b=3.09,a=0.87 mm。同理可得,Mini-SiC/SiC B的b=2.50,a=1.16 mm。此外,Weibull分布参数和线性相关系数r如表4所示。可知,线性相关系数都接近于1,变量之间存在高度相关关系,线性拟合合理,a与b接近真实数据。
表 4 两种Mini-SiC/SiC复合材料断裂位移双参数Weibull分布参数、线性相关系数Table 4. Weibull distribution parameters and linear correlation coefficients of two Mini-SiC/SiC composites with two-parameter fracture displacementName Shape parameter b Scale parameter a/mm r Mini-SiC/SiC A 3.09 0.87 0.97 Mini-SiC/SiC B 2.50 1.16 0.99 代入参数后,双参数Weibull累积损失因子损伤函数可分别表示为
Da=1−exp[−(x0.87)3.09] (18) Da=1−exp[−(x1.16)2.50] (19) 图7为拉伸强度的Weibull累积损伤因子函数拟合结果与试验结果对比图。分析可知,中位估计数据分别和式(18)、式(19)函数曲线吻合较好,因此,双参数Weibull模型能较好地表征累积损伤因子分布规律。
DN为Kolmogorov距离;D60,0.01为样本总量为60,显著水平α取0.01;D60,0.05为样本总量为60,显著水平α取0.05。
采用K-S检验进行拟合优度分析,计算估计值D和Da最大差异值可获得DN。如表5所示,经比对发现,两种材料DN值均小于D60,0.01、D60,0.05,故断裂位移均服从双参数Weibull分布,且拟合优度良好。
表 5 Mini SiC/SiC复合材料的DN与D60,0.01、D60,0.05Table 5. DN and D60,0.01, D60,0.05 of Mini SiC/SiC compositesName DN D60,0.01 D60,0.05 Mini-SiC/SiC A 0.1139 0.2104 0.1756 Mini-SiC/SiC B 0.1036 0.2104 0.1756 Notes: DN—Kolmogorov distance; D60,0.01—Total sample size is 60 and the significance level α is taken as 0.01; D60,0.05—Total sample size is 60 and α is taken as 0.05. 2.3.2 拉伸载荷与位移关系
将式(13)代入式(16)得拉伸载荷-位移关系模型:
P=k(exp[−(xa)b])x (20) 为求解k值,假设拉伸载荷-位移为线性关系。将每个试样断裂位移和最大载荷数据代入式(20),可得每个试样值,取平均值,即为k值。经计算,两种材料k值分别为370.89 N/mm、251.09 N/mm。将a、b、k值代入式(20),得Mini-SiC/SiC A、Mini-SiC/SiC B拉伸载荷-位移模型分别为
P=370.89(exp[−(x0.87)3.09])x (21) P=251.09(exp[−(x1.16)2.50])x (22) 将拉伸应力和载荷关系(式(11))与拉伸应变和加载位移之间的关系(式(12))及拉伸载荷-位移模型(式(20))联立进行求解,得拉伸应力-应变关系如下式:
σ=kexp[−(εl0/a)b]εl0S (23) 代入相关数据得上述两种材料拉伸应力-应变关系分别如下式所示:
σ=370.89exp[−(50ε/0.87)3.09]50εS (24) σ=251.09exp[−(50ε/1.16)2.50]50εS (25) 为验证拉伸应力-应变关系合理性,抽取拉伸强度接近于平均值的试样进行检验。Mini-SiC/SiC A中选取最大载荷128.79 N,断裂位移0.83 mm的试样。Mini-SiC/SiC B中选取最大载荷132.60 N,断裂位移为1.03 mm的试样。求出试样断裂应变εmax后,将该值分别代入式(24)、式(25)可获得采用拉伸应力-应变关系计算强度值(表6)。比对发现,与真实强度误差不超过7%。因此,该关系建立合理,可预测应力-应变关系。
表 6 Mini-SiC/SiC复合材料拉伸应力-应变关系检验过程主要数据Table 6. Main data of the test process of tensile stress-strain relationship of Mini-SiC/SiC compositesName Mini-SiC/
SiC AMini-SiC/
SiC BMaximum load P/N 128.79 132.60 Fracture displacement xmax/mm 0.83 1.03 True strength value σ/MPa 368.72 393.48 Fracture strain εmax 0.0166 0.0205 Strength calculation value of con-
stitutive model/MPa373.50 366.40 Error between calculated value
and real value+1.2% –6.8% 2.4 基体对复合材料强度分散性的影响
源于CVI工艺的复杂性,故复合材料内部存在缺陷,这些缺陷导致复合材料宏观力学性能具有分散性。图8为Mini-SiC/SiC复合材料初始态与断口SEM图像。图8(a)为初始态单向纤维,可看出纤维束表面沉积基体较粗糙,存在参差不齐菜花状β-SiC,表明基体沉积状态具有不稳定性。由图8(b)、8(c)可见,Mini-SiC/SiC复合材料断口纤维呈纤维单丝形式拔出,基体裂纹未出现大范围偏转,表现瞬间断裂。当纤维束表面存在缺陷时,拉伸过程中会由于应力集中而成为断裂源,造成纤维早期破环,故降低纤维表面缺陷可提高纤维束强度,从而降低强度的分散性。
图9采用ORS Dragonfly软件对Mini-SiC/SiC复合材料拉断CT图进行处理,图9(a)为初始状态下的正视断口,在阈值分割基础上采用深度学习方法提取拉伸断口裂纹分布(图9(b))情况,其中近似垂直于纤维且相互平行的线条为复合材料周期性裂纹。接着提取SiC纤维(图9(c)),可看出纤维分布并不均匀,纤维少的地方为基体与孔隙填充。随后,对断口俯视角度也进行学习处理(图9(d)),先提取SiC基体(图9(e)),可清楚看到纤维与基体分布状况,证实纤维少的地方由基体填充。最后,提取纤维表面BN界面相(图9(f)),可以看出界面相均匀覆盖纤维表面,形成一根根圆柱管。
由图8和图9从Mini-SiC/SiC复合材料宏观表面和分离复合材料内部各组元来看,制备态表面状态与复合材料内部纤维分布情况,会影响界面相与基体分布状态 ,导致复合材料分散性大,影响其可靠性。
复合材料拉伸断裂强度σu如下式所示:
σu=σfuVf (26) 式中:
σfu 为纤维断裂强度;Vf 为纤维体积分数。由式(26)求得理论拉伸强度为499.45 MPa,但由于单向复合材料断裂时,连续纤维(l>30lc)并非同时断裂,下式中载荷传递因子β为0.7,其纤维承载效率大于99%:
σfu=σ∗fu(1−1−βl/lc) (27) 式中:
σ∗fu 为纤维的原位断裂强度;l/lc 为比长度;l为纤维长度;lc为临界纤维长度。此时求得理论拉伸强度为349.61 MPa,该值与实测拉伸强度365.08 MPa相差仅4.24%。表明纤维基本全部承载,纤维对复合材料分散性影响较小。
图10对Mini-SiC/SiC复合材料拉伸产生的周期性裂纹间距进行测量,最小裂纹间距为83.2 μm,最大裂纹为107.8 μm。根据图9(f)可知界面相均匀分布,故视界面滑移应力(τ)为定值,不影响复合材料的分散性。由式(28)可计算出(附录表8、表9) Mini-SiC/SiC复合材料基体裂纹宽度为89.43 μm,该宽度处于实际裂纹宽度范围内。验证了该式可以较好地描述Mini-SiC/SiC复合材料力学断裂行为。
将最小与最大裂纹间距代入式(25),分别求出其对应拉伸强度为376.16 MPa和406.14 MPa,该值同样处于该批试样实际拉伸强度(331.02 MPa,407.82 MPa)内,表明基体裂纹间距具有离散性,基体是影响Mini-SiC/SiC复合材料可靠性的主要来源。
ς=b2(1−a1Vf)2Rσp24dτ0EmVf2 (28) 式中:Vf为纤维体积分数;
σp 为峰值应力;R为纤维半径;d为裂纹间距;Em为基体模量。通过双参数Weibull分布模型可知,小试炉与中试炉拟合良好,故裂纹间距仍选用图10。此时中试炉试样拉伸强度(附录表6)范围为(161.09 MPa,540.95 MPa),该范围大于小试炉试样拉伸强度范围,这是由于中试炉沉积基体工艺的放大化所导致,这种现象与m值对应,表明中试炉具有更大离散性(图3)。最大和最小裂纹间距对应的拉伸强度仍处于该范围内。
图11为基体影响复合材料形状参数m与力学性能离散的机制图。根据图9(f)与式(26)、式(27)可知,纤维与界面对复合材料性能离散的影响不大,故基体为影响复合材料可靠性的重要原因。m值影响应力-应变曲线非线性段的平稳性[21],m越大非线性段越平稳,表明复合材料在达到基体最大开裂后曲线越规律,离散性越小。图3与表2可明显看出小试炉制备的复合材料更加稳定,这是由于在小试炉中CVI沉积复合材料基体更加均匀,复合材料孔隙分布集中在内部中心区域(图11(a)),该区域的形成主要原因是基体在沉积时会形成沉积梯度,靠近表层区域的部分复合材料密度大、孔隙被充分填充,致使复合材料内部孔隙无法填充。反观中试炉制备,CVI工艺得到放大,同批次制备试样大大增多,造成复合材料内部孔隙发生不均匀分布现象(图11(b)),进而导致制备出的试样离散性增大、m值减小,影响其可靠性。
3. 结 论
(1) 双参数Weibull分布模型可准确表征小试炉制备单向纤维束 SiC/SiC (Mini-SiC/SiC A)、中试炉制备Mini-SiC/SiC B拉伸强度分布,其拉伸强度分别为365.08 MPa和367.78 MPa,形状参数分别为20.59和5.01。
(2) 基于双参数Weibull分布模型建立的可靠度分析模型可准确预测不同可靠度下的可靠强度。其中,当可靠度为50%时,可靠强度与拉伸强度最大误差为1.2%,故在工程中可采用材料的平均拉伸强度预测50%可靠拉伸强度。
(3) 基于双参数Weibull分布模型建立的拉伸载荷-位移关系可准确预测拉伸应力-应变关系,与真实强度误差不超过7%。
(4) 复合材料内部界面相均匀包裹纤维,界面相滑移应力视为定值;连续纤维承载大于99%,纤维对复合材料性能离散性影响较小;通过基体裂纹间距公式可知,最小和最大裂纹间距所对应的拉伸强度为376.16 MPa和406.14 MPa,该值均处于Mini-SiC/SiC A复合材料拉伸强度(331.02 MPa,407.82 MPa)和Mini-SiC/SiC B复合材料拉伸强度(161.09 MPa,540.95 MPa)范围内。然而,Mini-SiC/SiC B复合材料拉伸强度范围大于Mini-SiC/SiC A,且前者 Weibull模数5.01比后者 Weibull模数 20.59低75.7%,中试分散性更大是由于中试炉尺寸放大所致,基体非均匀性是影响复合材料可靠性的主要原因。
-
表 1 复合材料各组分性能参数
Table 1 Performance parameters of each component of composite materials
Material Type Tensile strength/MPa Tensile modulus/GPa Elongation at break/% Density/(g·cm−3) Carbon fiber (CF) T300-3 K 3530 230 2.1 1.8 Glass fiber (GF) E 3400 73 4.8 2.54 Bismaleimide R801 90 4.2 2.5 1.25 表 2 2.5D机织混杂复合材料参数
Table 2 2.5D woven hybrid composite material parameters
Yarn Fiber linear
density/texLayer Preformed unit
density/(yarn·cm−1)Composite
thickness/mmFiber volume
fraction/vol%Warp (CF) 200 2 8 1.85 49.87 Warp (GF) 144 4 8 Weft (CF) 200 3 4 Weft (GF) 144 4 8 -
[1] 李昱霖, 安庆升, 杨坤好, 等. 防热承载一体化复合材料电缆罩分析及验证[J]. 空天防御, 2019, 2(3):1-7. DOI: 10.3969/j.issn.2096-4641.2019.03.001 LI Yulin, AN Qingsheng, YANG Kunhao, et al. Analysis and validation of integrated thermal protection and load bearing for composite material cable cover[J]. Air&Space Defense,2019,2(3):1-7(in Chinese). DOI: 10.3969/j.issn.2096-4641.2019.03.001
[2] 周世豪, 倪楠楠, 刘彬, 等. 结构热防护一体化复合材料研究进展[J]. 航空材料学报, 2022, 42(4):1-15. DOI: 10.11868/j.issn.1005-5053.2020.000053 ZHOU Shihao, NI Nannan, LIU Bin, et al. Research progress on structural thermal protection integration composite[J]. Journal of Aeronautical Materials,2022,42(4):1-15(in Chinese). DOI: 10.11868/j.issn.1005-5053.2020.000053
[3] LI W J, ZHANG Z W, ZHU M D, et al. Novel strategy and multi-scale modelling of integrated multifunctional composite for thermal protection under extreme environment[J]. Applied Thermal Engineering,2022,209:118313. DOI: 10.1016/j.applthermaleng.2022.118313
[4] 杨智勇, 张东, 顾春辉, 等. 国外空天往返飞行器用先进树脂基复合材料研究与应用进展[J]. 复合材料学报, 2022, 39(7):3029-3043. YANG Zhiyong, ZHANG Dong, GU Chunhui, et al. Research and application of advanced resin matrix composites for aerospace shuttle vehicles abroad[J]. Acta Materiae Compositae Sinica,2022,39(7):3029-3043(in Chinese).
[5] 高军鹏, 白江波, 邓华, 等. 间隙率对平纹及三轴向织物复合材料弹性性能的影响[J]. 宇航材料工艺, 2014, 44(5):20-24,35. DOI: 10.3969/j.issn.1007-2330.2014.05.005 GAO Junpeng, BAI Jiangbo DENG Hua, et al. Effect of gap ratio on elastic properties of plain weave fabric and laminate with triaxial weave fabric composites[J]. Aerospace Materials & Technology,2014,44(5):20-24,35(in Chinese). DOI: 10.3969/j.issn.1007-2330.2014.05.005
[6] FAN W, DANG W, LIU T, et al. Fatigue behavior of the 3D orthogonal carbon/glass fibers hybrid composite under three-point bending load[J]. Materials & Design,2019,183:108112.
[7] PANKOW M, RIOSBASS M, JUSTUSSON B, et al. Effect of fiber hybridization on tensile fracture of 3D woven textile composites[J]. Composites Part A: Applied Science and Manufacturing,2023,165:107306. DOI: 10.1016/j.compositesa.2022.107306
[8] 龚亮, 薛利利. 热氧老化对三向正交碳/玻璃纤维/双马复合材料力学性能的影响[J]. 材料工程, 2021, 49(10):144-155. DOI: 10.11868/j.issn.1001-4381.2021.000045 GONG Liang, XUE Lili. Influence of thermo-oxidative aging on mechanical performance of three-dimensional orthogonal carbon/glass fiber/bismaleimide composites[J]. Journal of Materials Engineering,2021,49(10):144-155(in Chinese). DOI: 10.11868/j.issn.1001-4381.2021.000045
[9] PAPA I, BOCCARUSSO L, LANGELLA A, et al. Carbon/glass hybrid composite laminates in vinylester resin: Bending and low velocity impact tests[J]. Composite Structures,2020,232:111571. DOI: 10.1016/j.compstruct.2019.111571
[10] 王忠远, 蔡长春, 王振军, 等. 三维角联锁机织铝基复合材料面内拉伸力学行为与失效机制[J]. 复合材料学报, 2021, 38(9):2989-2999. DOI: 10.13801/j.cnki.fhclxb.20201116.007 WANG Zhongyuan, CAI Changchun, WANG Zhenjun, et al. In-plane tensile mechanical behavior and failure mechanism of three-dimensional angle interlocking woven aluminum matrix composites[J]. Acta Materiae Compositae Sinica,2021,38(9):2989-2999(in Chinese). DOI: 10.13801/j.cnki.fhclxb.20201116.007
[11] LI T, CHEN L, SHEN W, et al. Quasi-static compression response of carbon fiber reinforced 2.5D woven composites at different loading directions[J]. Materials,2022,15(11):3953. DOI: 10.3390/ma15113953
[12] 翟军军, 王露晨, 孔祥霞. 三维多向编织复合材料温度效应综述:热传导、热膨胀性质和力学响应[J]. 复合材料学报, 2021, 38(8):2459-2478. ZHAI Junjun, WANG Luchen, KONG Xiangxia. A review on temperature effect of 3D multi-directional braided composites: Heat conduction, thermal expansion properties and mechanical response[J]. Acta Materiae Compositae Sinica,2021,38(8):2459-2478(in Chinese).
[13] DONG F, WANG X, ZHANG C, et al. Ablation behavior and damage mechanisms of carbon/boron-modified phenolic 2.5D woven composite[J]. Polymer Degradation and Stability,2023,209:110279. DOI: 10.1016/j.polymdegradstab.2023.110279
[14] RUGGLES-WRENN M B, ALNATIFAT S A. Fully-reversed tension-compression fatigue of 2D and 3D woven polymer matrix composites at elevated temperature[J]. Polymer Testing,2021,97:107179. DOI: 10.1016/j.polymertesting.2021.107179
[15] ZHAO Y R, JIAN S, WEN W D, et al. Thermo-mechanical behaviors of 2.5D shallow straight-link-shaped woven composites under the warp direction fatigue loading at room and elevated temperatures[J]. Composite Structures,2022,289:115489. DOI: 10.1016/j.compstruct.2022.115489
[16] SONG J, WEN W, CUI H. Experimental and numerical investigation of mechanical behaviors of 2.5D woven composites at ambient and un-ambient temperatures[J]. Composite Structures,2018,201:699-720. DOI: 10.1016/j.compstruct.2018.06.054
[17] DANG M G, LI D S, LEI J. Temperature effects on mechanical response and failure mechanism of 3D angle-interlock woven carbon/epoxy composites[J]. Composites Communications,2020,18:37-42. DOI: 10.1016/j.coco.2020.01.001
[18] 曹淼. 三维角联锁机织碳纤维/环氧树脂复合材料热氧老化效应及冲击后剩余强度[D]. 上海: 东华大学, 2020. CAO Miao. Residual strength of 3-D angle-interlock woven carbon fiber/epoxy resin composites after thermo-oxidative ageing and impact[D]. Shanghai: Donghua University, 2020(in Chinese).
[19] RATHORE D K, PRUSTY R K, MOHANTY S C, et al. In-situ elevated temperature flexural and creep response of inter-ply glass/carbon hybrid FRP composites[J]. Mechanics of Materials,2017,105:99-111. DOI: 10.1016/j.mechmat.2016.11.013
[20] 于洋, 樊威, 薛利利, 等. 热氧老化对三维编织碳纤维-玻璃纤维/双马来酰亚胺树脂复合材料力学性能的影响[J]. 复合材料学报, 2021, 38(12):4060-4072. DOI: 10.13801/j.cnki.fhclxb.20210311.004 YU Yang, FAN Wei, XUE Lili, et al. Influence of thermo-oxidative aging on the mechanical performance of three-dimensional braided carbonfiber-glass fiber/bismaleimide composites[J]. Acta Materiae Compositae Sinica,2021,38(12):4060-4072(in Chinese). DOI: 10.13801/j.cnki.fhclxb.20210311.004
[21] LIANG Q, LIU J, WANG X, et al. Flexural progressive failure mechanism of hybrid 3D woven composites: Combination of X-ray tomography, acoustic emission and digital image correlation[J]. Composite Structures,2022,280:114894. DOI: 10.1016/j.compstruct.2021.114894
[22] 张雷, 李姗姗, 庄毅等. 碳纤维与玻-碳层间混杂2.5维机织复合材料的力学性能对比研究[J]. 材料导报, 2022, 36(S1):199-203. ZHANG Lei, LI Shanshan, ZHUANG Yi, et al. Comparative study on mechanical properties of carbon fiber and glass-carbon interlaminar hybrid 2.5D woven composites[J]. Materials Reports,2022,36(S1):199-203(in Chinese).
[23] 中国国家标准化管理委员会. 纤维增强塑料弯曲性能试验方法: GB/T 1449—2005[S]. 北京: 中国标准出版社, 2005. Standardization Administration of the People's Republic of China. Fibre-reinforced plastic composites—Dertermination of flexural properties: GB/T 1449—2005[S]. Beijing: China Standards Press, 2005(in Chinese).
[24] 中国建筑材料联合会. 纤维增强塑料 短梁法测定层间剪切强度: JC/T 773—2010[S]. 北京: 中国建材工业出版社, 2010. China Building Materials Federation. Fibre-reinforced plastics composites—Determination of apparent interlaminar shear strength by short-beam method: JC/T 773—2010[S]. Beijing: China Building Materials Industry Press, 2010(in Chinese).
[25] 李紫伦, 杨安坤, 覃小红, 等. 三维编织玻璃纤维/环氧树脂复合材料薄壁管轴向压缩性能的温度效应[J]. 复合材料学报, 2023, 40(10): 5588-5600. LI Zilun, YANG Ankun, QIN Xiaohong, et al. Temperature effect on axial compressive properties of three-dimensional glass fiber/epoxy resin braided composite thin-walled tubes[J]. Acta Materiae Compositae Sinica, 2023, 40(10): 5588-5600(in Chinese).
[26] LI J, FAN W, TAO L, et al. The temperature effect on the inter-laminar shear properties and failure mechanism of 3D orthogonal woven composites[J]. Textile Research Journal,2020,90(23-24):2806-2817. DOI: 10.1177/0040517520927009
[27] YANG B, YUE Z, GENG X, et al. Temperature effects on transverse failure modes of carbon fiber/bismaleimides composites[J]. Journal of Composite Materials,2017,51(2):261-272. DOI: 10.1177/0021998316639122
[28] 双超, 刘璐璐, 赵振华, 等. 湿热老化后碳纤维复合材料层间剪切强度实验方法对比研究[J]. 航空材料学报, 2017, 37(5):90-98. DOI: 10.11868/j.issn.1001-4381.2016.000226 SHUANG Chao, LIU Lulu, ZHAO Zhenhua, et al. Comparison study on interlaminar shear strength testing methods of CFRP under hygrothermal aging conditions[J]. Journal of Aeronautical Materials,2017,37(5):90-98(in Chinese). DOI: 10.11868/j.issn.1001-4381.2016.000226
[29] HE Y, MEI M, WEI K, et al. Interlaminar shear behaviour and meso damage suppression mechanism of stitched composite under short beam shear using X-ray CT[J]. Composites Science and Technology,2022,218:109189. DOI: 10.1016/j.compscitech.2021.109189
-
期刊类型引用(1)
1. 李新娅,王宁,卢佳浩,张鹏,夏兆鹏,侯耒. 基于改进Weibull模型的高强缝合锚钉缝线强度预测. 现代纺织技术. 2024(06): 52-60 . 百度学术
其他类型引用(1)
-
目的
随着超高速飞行器马赫数的提升,对飞行器材料结构件在缩小体积、减轻重量、提高力学承载性能、抗静电性能和防热性能等方面提出了更为苛刻的要求。传统的飞行器耐高温与承载部件大多采用金属材料,并在表面涂覆涂层来实现防热,这种结构存在质量大、隔热差以及耐用性难保证等缺点,难以满足结构-功能一体化的设计需求。基于此,本文设计并制备了2.5D机织碳纤维-玻璃纤维/双马来酰亚胺树脂复合材料(2.5D机织混杂复合材料),并对其高温力学性能展开了研究。
方法采用三维纺织技术与树脂传递模塑工艺(RTM)制备了2.5D机织碳纤维-玻璃纤维/双马来酰亚胺树脂复合材料,并根据复合材料的动态热机械性能设定了力学性能测试温度上限,开展了不同温度场(25℃、150℃、240℃、300℃)下2.5D机织混杂复合材料的三点弯曲和层间剪切力学性能测试。借助光学显微镜和扫描电镜对试样的断口形貌进行了观测,以探究温度对复合材料力学行为及损伤机制的影响。
结果动态热机械测试(DMA)结果表明2.5D机织混杂复合材料具有优异的耐高温性能。通过DMA测得该复合材料玻璃化转变温度为304℃,当测试温度升至300℃时,其储能模量保留率仍有72.32%,表现出良好的热稳定性。此外,2.5D机织混杂复合材料具有明显的温度效应,复合材料的弯曲强度、弯曲模量和层间剪切强度随温度上升逐渐下降。2.5D机织混杂复合材料的三点弯曲失效为压缩应力、拉伸应力和剪切应力耦合作用下的结果,整个破坏是一个从受压面外侧到受拉面内侧的渐进过程。试样在150℃、240℃和300℃下的平均弯曲强度相较于25℃下分别降低了15.31%、18.42%和23.06%;平均弯曲模量相较于25℃下分别降低了10.06%、32.05%和70.01%。从破坏形貌来看,随着温度升高,试样受压面和受拉面上的损伤逐渐减弱,从侧面观察到试样的损伤主要集中在试样的上半部分并沿经纱方向扩展,主要的损伤模式包括纱线断裂、基体裂纹、纤维/基体脱粘和分层等。2.5D机织混杂复合材料层间剪切失效为层间剪切应力、弯曲应力、横向剪切应力和局部挤压应力多力耦合作用下的结果。试样在150℃、240℃和300℃下试样的平均弯曲强度相较于25℃下分别降低了4.92%、7.47%和18.93%,可以看出,试样在各温度场下的剪切强度变化不大,表现出优异的层间剪切性能。从破坏形貌来看,短梁试样上半部分受压出现纱线与基体界面开裂,下半部分呈拉伸破坏并存在明显的分层现象。随着温度升高,界面退化引起损伤范围增大,导致试样最终失效。主要的损伤模式为,纱线断裂,基体开裂,界面脱粘及分层等。
结论温度对2.5D机织碳纤维-玻璃纤维/双马来酰亚胺树脂复合材料的力学特性和损伤模式均有着显著影响。温度上升导致纤维/基体界面结合力减弱,复合材料在300℃下的弯曲强度、弯曲模量和层间剪切强度相比于室温环境分别降低了23.06%、70.01%和18.93%。弯曲载荷下,2.5D机织混杂复合材料的室温破坏模式主要以局部的纤维断裂和基体开裂为主,而高温损伤则以纤维/基体界面脱粘为主导。剪切载荷下,2.5D机织混杂复合材料的室温破坏模式主要为分层破坏,而随着温度升高,复合材料因基体软化出现塑性变形,基体开裂、界面脱粘及分层破坏决定了材料的最终失效。
-
随着超高速飞行器马赫数的提升,对飞行器材料结构件在缩小体积、减轻重量、提高力学承载性能、抗静电性能和防热性能等方面提出了更为苛刻的要求,碳/玻璃纤维混杂不仅在力学性能上可以做到纤维间性能取长补短,而且还可以对混杂纤维进行合理的排布和编织,实现复合材料的结构-功能一体化。目前混杂结构多为二维层合结构,在受到载荷时容易发生分层破坏,2.5D机织结构复合材料采用一体化成型技术,具有层间断裂性能优异、可设计性强和近净成形等优势。在复合材料部件服役过程中,不可避免面临高温环境的威胁,研究复合材料在高温环境下的力学行为具有重要的理论和实践意义,但2.5D机织混杂复合材料的高温力学行为研究鲜有报道。
基于此,本文设计制备了2.5D机织碳纤维-玻璃纤维/双马来酰亚胺树脂复合材料,并根据复合材料的动态热机械性能设定了力学性能测试温度上限,开展了不同温度场(25℃、150℃、240℃、300℃)下复合材料三点弯曲和层间剪切力学性能测试,研究了温度效应下复合材料的力学行为及损伤机制。研究结果表明:2.5D机织混杂复合材料具有明显的温度效应,随着温度升高,基体和纤维结合力减弱,复合材料的三点弯曲和层间剪切力学性能逐渐下降,但在300℃下仍具有较高性能水平,损伤模式从局部的纤维断裂和基体破坏转变为大面积的纤维/基体界面脱粘。
不同温度下2.5D机织混杂复合材料层间剪切力学性能 (a)载荷-位移曲线;(b)剪切强度