复合材料反平面切口端部奇性研究.pdfVIP

第30卷第5期 计 算 物 理 V01．30．No．5 CHINESE OFCOMPUTATIONALPHYSICS 2013年9月 JOURNAL Sep．，2013 文章编号：1001-246X【2013)05-0700-06 复合材料反平面切口端部奇性研究 程长征， 葛仁余， 薛伟伟， 牛忠荣 (合肥工业大学土木与水利工程学院，安徽合肥230009) 摘要：基于复合材料切121尖端位移场的渐近展开，将切121的反平面平衡控制方程转化为关于切口奇性指数的特 征微分方程，采用一种变换将其化为线性特征微分方程组，引入插值矩阵法计算相应边界条件下方程组的特征值 以获取切121尖端的应力奇性指数．研究单相材料切口、双相材料切121及止于异质界面切121的奇异性，算例表明该方 法可以一次性计算出多阶奇性指数．对所取得的非奇异指数尽管切口不表现出奇性状态，但却是描述切口尖端完 整应力场必不可少的参量， 关键词：复合材料；反平面切口；奇性分析；渐近展开 中图分类号：0343．6 文献标识码：A 0 引言 复合材料由于其优异的力学性能被广泛应用于航天、国防等高技术工业领域，但在其加工和使用过程 中，不可避免会出现接头、界面端甚至裂纹等几何或材料不连续区域，这里统称为切口问题．切口端部由于几 何形状或材料特性的突变导致尖端区域发生应力奇异¨1．李航等旧。通过对含穿透切口的复合材料加筋板轴 向压缩试验，发现切口前端出现了非常明显的应变集中且最先发生破坏．Liu等1研究发现对于一些特殊的 复合材料切口角度，通过选择合适的纤维方向，应力奇异性的程度可以减小甚至消失．也就是说，开展切口奇 性研究对确保复合材料结构服役的可靠性和安全性具有重要意义． 对于复合材料反平面切13奇异性的研究，采用复变函数和Mellin变换方法的有李俊林等”3采用断裂复 变方法，研究了正交异性双材料反平面平板搭接界面端问题，显示反平面搭接界面端只有一阶奇异性． 强度的表达式。Lin等¨o采用联合Mellin变换和镜像法分析了二维异质材料切口在反平面荷载和螺旋位错 了不同径向边界条件下切口尖端位移和应力场的表达式；采用有限元法的有：Pageau等旧1采用有限元特征 法研究了平面各向异性复合材料尖劈在反平面载荷作用下的奇性应力状况．平学成等1提出了一个基于位 移的分析尖劈端部奇性位移场和应力场反平面问题的非协调元特征法；采用特征函数法的有：钱俊等¨州先 用特征函数法导出了反平面复合材料切Ill的特征方程，后采用MUller法计算了解的序列．戴耀等¨u将均匀 材料的特征函数结构推广到非均匀材料，并以反平面问题为例，给出了该问题任意阶特征函数的定解方程． 尽管上述有些方法能给出切口奇性指数的表达式旧川，但那是一组超越方程，还需要迭代求解．本文提 出基于渐近展开，将切口尖端的位移场表达成关于奇性指数和位移角函数的组合形式，取其中的典型项代人 复合材料结构反平面问题的控制方程，化切口奇性指数的计算为求一组常微分特征方程组的解，这样可以一 次性计算出所有的奇性指数，而不需要一次一次地迭代计算¨。． 1 单相复合材料反平面切口 图1所示为单相复合材料反平面切口问题．在整体坐标系oxy下反平面应力分量or。or。和反平面应变 收稿日期：2012—11—19；修回日期：2013—04一05 基金项目：国家自然科学基金11372094)和安徽省自然科学基金38)资助项目 万方数据 第5期 程长征等：复合材料反平面切口端部奇性研究 701 分量y。7，：之间的本构关系可表示为 盼R挑≯ ㈩ 其中G。，和G：，为复合材料的切变模量．经过坐标变换后，可以获 得极坐标系orO下的本构关系 盼暖跚12]fr。rz)， (2) 其中 Q1l=COS2(0—00)G13+sin2(日一00)G23， (3) Q1