文章编号: 1004-0609(2006)02-0253-07
纤维拔出时界面分离能释放率的应力函数分析
刘鹏飞, 陶伟明, 郭乙木
(浙江大学 工程力学系, 杭州 310027)
摘 要: 当引入应力平衡方程、 边界条件和连续性条件时, 提出了一个描述纤维拔出过程中、 分离界面伴随摩擦时的弹性应力传递理论模型。 根据最小势能原理, 获得了纤维和基体中所有应力的求解和考察界面断裂特性的分离能释放率G的表达式, 且当引入一个界面分离准则G≥Γi时, 提供了一个确定临界分离长度的方法。 对纤维增强复合材料SiC/Ti-6Al-4V进行了数值计算, 且和shear-lag(剪滞)理论结果作比较。
关键词: SiC纤维增强复合材料; 纤维拔出; 应力分析; 能量释放率
中图分类号: TG332 文献标识码: A
Analysis of energy release rate for interfacial debonding in fiber pull-out by using stress functions
LIU Peng-fei, TAO Wei-ming, GUO Yi-mu
(Department of Engineering Mechanics, Zhejiang University, Hangzhou 310027, China)
Abstract: By introducing stress equilibrium equations, boundary-and continuity-conditions, a theoretical model on the elastic stress transfer with friction at the debonded interface in fiber pull-out was established. Based on the minimum complementary energy principle, all stress solutions in the fiber and matrix, an expression for the energy release rate G which explores the interfacial fracture properties were obtained. When an interfacial debonding criterion G≥Γi was introduced, a method for determining the critical debond length was provided. Numerical calculations were conducted for the fiber-reinforced composite SiC/Ti-6Al-4V and the results were compared with those obtained by using the shear-lag theories.
Key words: SiC fiber-reinforced composites; fiber pull-out; stress analysis; energy release rate;
SiC纤维增强金属基复合材料由于具有高强度、 耐高温和高断裂韧性等优点, 在航空航天等领域具有广泛的应用价值, 而其中一个重要的增韧机理是纤维桥联基体裂纹效应。 桥联过程中分离界面上伴随着复杂的摩擦和滑移, 很大程度上影响复合材料的力学特性[1-4]。
由于纤维施加在基体裂纹上的桥联载荷对于纤维来说是拔出载荷, 因此对纤维拔出的微观力学测试已成为考察界面断裂和失效特性的广泛使用的方法[5]。 为分析增韧效应, 在纤维拔出过程中, 建立一个有效的界面分离准则是必要的。 当前使用较多的界面分离准则是以界面剪切强度作为参考, 通过计算最大界面剪切应力来判断界面是否进一步分离[6]。 按照基于能量的界面分离准则, 当界面分离能释放率G超过界面分离韧性Γi时, 界面裂纹扩展。 已经证明, 基于能量的分离准则在预测界面失效时较强度准则更可靠[7]。
显然, 准则地计算纤维和基体中的应力和应变场对于获得一个准确的能量释放率的求解是很关键的。 分别运用剪滞理论和拉梅方法, Hsueh[8] 和Ochiai等[9]获得了分离能释放率的解, 然而由于忽略了纤维中的剪切应力和应变能、 界面径向应力、 基体轴向应力随径向位置的变化及泊松效应等, 具有一定缺陷。 当基体轴向应力被集中于一个有效半径上的一个等效轴向应力代替时, Chiang[10]推导出包括纤维轴向应变能、 基体轴向和剪切应变能在内的分离能释放率的表达式。 Rauchs等[11]运用有限元法, 获得了能量释放率的数值解, 然而过度的简化导致较大误差。
本文作者采用应力平衡方程, 引入应力函数来分析界面分离准则, 在最小势能原理的基础上, 获得了界面分离且伴随摩擦滑移时、 纤维和基体中的各项应力和能量释放率的求解。 数值结果和Hsueh[8] 、 Chiang[10]结果作比较。
1 界面分离准则
1.1 应力分析
图1所示为一个半径为r1、 嵌入外半径为r2的基体中的单纤维。 纤维体积分数φf=r21/r22, 基体体积分数φm=1-φf。 嵌入基体的纤维长度L, 界面分离长度Ld。 纤维和基体均考虑成线弹性。 定义一个柱坐标系(r, θ, z), z轴为纤维轴向, 纤维加载端和嵌入端分别为z=0和z=L。 纤维拔出应力σb平行于z轴。
图1 拔出应力作用下分离界面上的摩擦滑移
Fig.1 Frictional sliding at debonded interface under pull-out stress
纤维轴向应力σzf和界面剪切应力τi(z)满足
式中 σzf(z)考虑为纤维横截面上的平均轴向应力。
纤维和基体轴向应力之间的平衡要求为
在分离界面上, 假设剪切应力τi(z)退化为摩擦应力τs, 则联合式(1)和(2)获得
在柱坐标系下, 对于一个空间轴对称问题, 忽略体积力时的应力平衡方程可表达为
按Wu等[12], 当引入一个应力函数Φi=fj(f)×gj(z)时, 满足式(5)的应力解可表达为
式中 i=f, j=1和i=m, j=2分别表示纤维和基体。 r, θ, z分别表示径向、 环向和轴向。
因此, 界面结合区域Ld≤z〈L内的应力解可表达为
1.2 应力边界条件
式中 式(8a)和(8b)分别表示纤维和基体轴向应力在分离端z=Ld处的连续性; 式(8c)和(8d)分别表示纤维和基体轴向应力在z=L处的边界条件[10]; 式(8e)和(8f)分别表示纤维和基体径向和剪切应力在界面结合处的连续性; 式(8g)和(8h)分别表示基体径向和剪切应力在基体外表面r=r2处的边界条件。 σd为分离端纤维轴向应力[5]。 Ef, Em, Ec分别为纤维、 基体和复合材料弹性模量。
1.3 应力求解
求解式(5b), 纤维剪切应力τrzf(r, z)计算为
根据式(1), (2), (5a), (8f)和(8h), 基体剪切应力τrzm(r, z)计算为
将式(7a)代入(8a), 函数f1(r)和f2(r)可表达为
式中 Ci(i=1-4)为常系数。
联合式(2)、 (7a)和(11), 则函数g1(z)和g2(z)之间的关系表达为
联合式(5), (8e)~(8h), (11)和(12), 则系数Ci和纤维、 基体中的所有应力表达为
根据式(13), 纤维和基体中的所有应力均由函数g1(z)来决定。
1.4 最小势能原理
对于一个稳定的平衡系统, 最小势能原理可表述为: 在所有几何可能位移中, 真实位移总是使总势能Π取最小值
式中 Ue为纤维和基体总应变能; Ued和Ueb分别为分离界面和结合界面的应变能; Uw为拔出应力σb作功; Uf为摩擦应力τs作功。
界面分离区域0≤z〈Ld内的应变能Ued可分为由远场拉伸应力σb和由摩擦应力τs引起的纤维和基体的应变能之和
式中 G为剪切模量, 且G=E/[2(1+v)]。 忽略了径向和环向应力的影响。
界面结合区域Ld≤z〈L内的应变能Ueb表达为
纤维和基体的应力—应变关系表达为
式中 v为泊松比。
联合式(16)~(18), 应变能Ueb可计算为
对于总势能Π, 欧拉-拉格朗日方程可表达为
式中 F为式(19)中的积分表达式
将表达式F代入式(20)中, 获得
求解式(21), 获得
根据式(8a)~(8d)和(24), 系数A和B分别计算为
当L→∞时, 系数A进一步计算为
因此, 界面结合区域Ld≤z〈L内的纤维轴向应力σzf(z)获得
联合式(1)和(27), 界面结合区域Ld≤z〈L内的剪切应力τi(z)获得
1.5 能量释放率
纤维和基体轴向位移分别计算为
式中 wf(z)和wm(z)方向均与z轴相反。 当L→∞时, 式(29)和(30)忽略了幂次项的影响。
纤维和基体之间的相对轴向位移v(z)计算为
摩擦应力τs作功Uf计算为
式中 ushear=τsr1/(2Gm)为基体剪切扭曲对裂纹展开位移的贡献; 为一个无量纲的参数, 可表达为[5]
拔出应力σb作功Uw可计算为
式中 Udebond为由界面分离导致的复合材料附加位移wf(z), 且可定义为纤维轴向位移在纤维加载端z=0的值与界面无分离时的复合材料位移wc之差
当界面裂纹扩展dLd时, 分离界面上的应变能以dA=2πr1dLd的因子变化。 联合式(14), (16), (19), (32)和(34), 能量释放率G的表达式获得
上式表明, 当材料和几何参数确定时, 能量释放率L是分离长度Ld的二阶函数。
与式(36)比较, Hsueh[8] 和Chiang[10]分别获得如下结果
2 结果与讨论
SiC纤维增强复合材料SiC/Ti-6Al-4V[13] 用于数值计算。 材料参数: Ef=400GPa, Em=115GPa, φf=0.17, φm=0.3[13]。 纤维拉伸强度σs=4.19GPa[14]。 纤维半径r1=10μm。
图2 不同条件下能量释放率随标准化分离长度的分布
Fig.2 Distributions of energy release rate with normalized debond length under different conditions
图2所示为随标准化分离长度Ld/r1的增加, 能量释放率G的分布。 理论上, 在式(36)~(38)中, 摩擦应力τs以Ecτ2s/(r1EfφmEm)的形式贡献于G, 因此, 随Ld/r1的增加, 产生了一条先下降后上升的曲线。 然而, 由于界面分离只发生在G>Γi这里, 界面分离韧性假设为Γi=1J/m2, 当条件G=Γi满足时, 界面分离停止, 因此G-Ld/r1曲线后上升部分实际上是无意义的。
对于G-Ld/r1曲线下降部分, 当G>Γi时, 本文获得的G值比Hsueh[8]和Chiang[10]稍小, 这可解释为: Hsueh忽略了纤维和基体中的剪切应力和应变能, 纤维和基体径向应力, 基体轴向应力随径向位置的变化及泊松效应。 Chiang只进一步考虑到基体剪切应力和应变能的影响。
比较图2(a)和(b), 随摩擦应力τs增加, 本文获得G值和Hsueh、 Chiang获得的G值之间的差距变大。 结论也能从式(36)获得解释: 由参数β′和远场纤维轴向应力σ∞f决定的纤维、 基体剪切项和泊松效应的强度均正比于摩擦应力τs。 因此, 随摩擦应力τs增加, 剪滞理论中忽略的纤维、 基体应变能和泊松效应, 抑制界面失效作用增大。
当引入一个界面分离准则G≥Γi时, 通过求解式(36), 临界分离长度能确定
由于G-Ld/r1曲线只有下降部分有意义, 因此上式中只有较小的Ld1(Ld1〈Ld2)是有意义的。 结论也由Liu等[15]证实, 区别在于: Liu获得的结论是建立在拉梅求解的基础上, 和剪滞理论具有同样的缺陷。 当时Γi=1J/m2, 本文获得的Ld1值较Hsueh和Chiang稍小, 在摩擦应力分别为τs=10, 20MPa时, 临界分离长度分别为Ld1=275, 135μm。
图3所示为摩擦应力τs对能量释放率G的影响。 对于G-Ld/r1曲线下降部分, 在界面分离停止之前, 随摩擦应力τs增加, 出现一段更长的分离长度, 且当Ld/r1相同时, 获得更小的G值。 当摩擦应力τs减小且逐渐接近零时, G-Ld/r1曲线逐渐接近水平, 在τs=0时, 分离能释放率获得最大值G=6.03J/m2, 此时界面抵抗失效的能力最弱。
图3 当σb=1GPa, φf=1%时, 不同摩擦应力下能量释放率随标准化分离长度的分布
Fig.3 Distributions of energy release rate with normalized debond length under σb=1GPa, φf=1%, and different friction stresses
图4 当σb=1GPa, τs=10MPa时,
不同纤维体积分数下能量释放率随标准化分离长度的分布
Fig.4 Distributions of energy release rate G with normalized debond length under σb=1GPa, τs=10MPa, and different fiber volume fractions
图4所示为纤维体积分数φf对能量释放率G的影响。 对于G-Ld/r1曲线下降部分, 当φf较小时, G值相差不大; 而当φf较大时, G值较大程度的减小, 抑制界面失效的能力增强。
3 结论
假设分离界面上的摩擦应力为常数, 建立了纤维拔出过程中伴随界面分离时, 界面弹性应力传递的应力函数求解模型。 基于最小势能原理, 分别获得了满足所有平衡条件、 界面连续性和边界条件的纤维和基体的应力解和界面分离能释放率的表达式。 对SiC纤维增强复合材料SCS-6/Ti-6Al-4V作数值计算, 分析了各种参数对界面分离特性的影响, 且结果和剪滞理论结果作比较。
1) 引入一个界面分离准则后, 一个求解临界界面分离长度的方法获得。
2) 当摩擦应力τs值一直减小且接近零时, G-Ld/r1曲线趋于水平, 直到界面无摩擦时, G值为一个独立于界面分离长度Ld的常数。
3) 理论上, G-Ld/r1曲线呈现先下降后上升的趋势, 然而只有在G>Γi时才发生界面分离, 当条件G=Γi满足时界面分离停止, 因此曲线后上升部分是无意义的。
4) 随摩擦应力τs增加, 纤维和基体剪切效应、 泊松效应抑制界面失效的作用增强。
REFERENCES
[1]张强, 陈国钦, 武高辉, 等.含高体积分数SiC的铝基复合材料制备与性能[J].中国有色金属学报, 2003, 13(5): 1180-1183.
ZHANG Qiang, CHEN Guo-qin, WU Gao-hui, et a1. Fabrication and property of SiCp/Al composites with high content of SiC[J]. The Chinese Journal of Non-ferrous Metals, 2003, 13(5): 1180-1183.
[2]Hampe A, Kalinka G, Meretz S, et al. An advanced equipment for single-fibre pull-out test designed to monitor the fracture process[J]. Composites, 1995, 26 (1): 40-46.
[3]Piggott M R. A new model for interface failure in fibre-reinforced polymers[J]. Composites Science and Technology, 1995, 55 (3): 269-276.
[4]Begley M R, McMeeking R M. Numerical analysis of fibre bridging and fatigue crack growth in metal matrix composite materials[J]. Mater Sci Eng A, 1995, A200(1-2): 12-20.
[5]Budiansky B, Cui Y L. Toughening of ceramics by short aligned fibers[J]. Mechanics of Materials, 1995, 21 (2): 139-146.
[6]Lacroix T H, Keunings R, Desaeger M, et al. A new data reduction scheme for the fragmentation testing of polymer composites[J]. Journal of Materials Science, 1995, 30 (1): 683-692.
[7]Hon da K, Kagawa Y. Debonding criterion in the pushout process of fiber-reinforced ceramics[J]. Acta Materialia, 1996, 44 (8): 3267-3277.
[8]Hsueh C H. Crack-wake interfacial debonding criteria for fiber-reinforced ceramic composites[J]. Acta Materialia, 1996, 44 (6): 2211-2216.
[9]Ochiai S, Hojo M, Inoue T. Shear-lag simulation of the progress of interfacial debonding in unindirectional composites[J]. Composites Science and Technology, 1999, 59 (1): 77-88.
[10]Chiang Y C. On fiber debonding and matrix cracking in fiber-reinforced ceramics[J]. Composites Science and Technology, 2001, 61(12): 1743-175.
[11]Rauchs G, Withers P J. Computational assessment of the influence of load ratio on fatigue crack growth in fibre-reinforced metal matrix composites[J]. International Journal of Fatigue, 2002, 24 (12): 1205-1211.
[12]Wu W, Verpoest I, Varna J. A novel axisymmetric variational analysis of stress transfer into fibres through a partially debonded interface[J]. Composites Science and Technology, 1998, 58 (12): 1863-1877.
[13]Preuss M, Rauchs G, Doel T J A, et al. Measurements of fibre bridging during fatigue crack growth in Ti/SiC fibre metal matrix composites[J]. Acta Materialia, 2003, 51(4): 1045-1057.
[14]Warrier S G, Maruyama B, Majumdar B S, et al. Behavior of several interfaces during fatigue crack growth in SiC/Ti-6Al-4V composites[J]. Mater Sci Eng A, 1999, A259(2): 189-200.
[15]Liu Y F, Kagawa Y. The energy release rate for an interfacial debond crack in a fiber pull-out model[J]. Composites Science and Technology, 2000, 60 (2): 167-171.
(编辑陈爱华)
基金项目: 浙江省自然科学基金资助项目(M503095)
收稿日期: 2005-06-23; 修订日期: 2005-11-08
作者简介: 刘鹏飞(1980-), 男, 博士研究生
通讯作者: 陶伟明, 教授; 电话: 0571-87951769; E-mail: taowm@zju.edu.cn