改进的三维虚拟裂纹闭合方法.pdfVIP

维普资讯 2008年2月 西 北 工 业 大 学 学 报 Feb． 2008 第26卷第l期 JournalofNorthwesternPolyteehnicalUniversity Vo1．26No．1 改进的三维虚拟裂纹闭合方法 谢 伟 ，黄其青 ，MasanoriKikuchi。 (1．西北工业大学 航空学院，陕西 西安 710072；2．东京理科大学 机械工程学院，日本 千叶 278—8510) 摘 要：改进 了一种基于有限元技术分析三维裂纹应力强度因子的计算方法——虚拟裂纹闭合方 法，称为改进 的虚拟裂纹 闭合方法。应用传统的虚拟裂纹闭合方法求解三维裂纹体应力强度 因子 时，对有限元模型具有一定的要求，即过裂纹前缘的有限单元裂纹面必须具有相同的面积且对称。 而应用文 中改进的虚拟裂纹闭合方法求解三维裂纹体应力强度 因子时，裂纹前缘的裂纹面可以是 任意形状。文中应用改进虚拟裂纹闭合方法，建立2O结点等参元有限元模型，对三维表面裂纹应力 强度 因子进行 了验证性的分析和讨论。通过与Newman解相比较 ，证明该方法适用性强且计算精 度高。 关 键 词：有限元，改进的虚拟裂纹闭合方法，裂纹面，应力强度因子 中图分类号 ：V131．5 文献标识码：A 文章编号：1000—2758(2008)01—0116—05 在线弹性断裂力学领域[1]，应力强度 因子和能 1 三维虚拟裂纹闭合方法 量释放率被用作断裂准则和表征裂纹尖端应力场的 控制参量 。应用数值计算方法计算这些裂尖控制参 应用 Irwin裂纹闭合积分方法 计算 图1所示 量一直是断裂力学领域重要的研究 内容 。 的裂纹前缘单元的能量释放率，其计算公式为 有一些利用有限元分析结果计算裂纹尖端应力 强度因子的方法，例如：直接利用有限元应力和位移 分析结果来计算应力强度因子的应力法和位移 法 “。；建立在结点力基础上计算应力强度因子的结 点力法 。还有一些利用能量释放率来计算应力强 度因子的方法，其能量释放率可 以通过 积分 方 法或裂纹闭合积分方法[5]获得。 在有限元分析基础上，为了利用裂纹闭合积分 概念，Rybicki和 KanninenE印首先提出虚拟裂纹 闭 合方法 (VCCM)来解决二维裂纹问题。该方法很快 被Shivakumar等[7]推广用来求解三维裂纹问题。但 是普通的三维虚拟裂纹闭合方法需要裂纹前缘的有 限元单元相互垂直且尺寸相等，对三维裂纹建立计 图1 三维虚拟裂纹闭合方法 算模型受到限制，本文发展了一种在裂纹前缘基于 考虑斜的不对称 的有限元单元的虚拟裂纹闭合方 1 广 法，该方法简单可行，对三维有限元网格没有限制， GTotal J ()vl(za—X1)dS (1) 适用性强。 式中， 和 。坐标方 向分别垂直和平行裂纹前缘， 收稿 日期 ：2007—02-06 作者简介：谢 伟(1978--)，西北工业大学博士生，主要从事结构疲劳、断裂与可靠性分析研究。 维普资讯 ． ． ． H ．．． 第 1期 谢 伟等：改进的三维虚拟裂纹 闭合方法 X。轴是裂纹面法线。△是垂直裂纹前缘的单元宽度 ， 2 改进 的虚拟裂纹闭合方法 ∥ 是平行裂纹前缘的单元长度。 表示裂纹前缘单 元编号。。是裂纹前方延伸面的应力。 是裂纹张开