CT试样三维疲劳裂纹扩展数值模拟.docVIP

作者简介:代鹏(1984-,男,宁夏中卫人,博士研究生,研究方向:疲劳裂纹。 通迅作者:冯淼林,男,副教授,(E-mail: mlfeng@sjtu.edu.cn CT 试样三维疲劳裂纹扩展数值模拟 (代鹏,冯淼林* (上海交通大学 船舶海洋与建筑工程学院,上海,200240 摘要:在循环载荷下疲劳裂纹的裂纹形貌在稳定扩展区近似为半椭圆形状,因此通过Paris 方法根据疲劳裂纹表面尖端点应力强度因子的变化幅值(ΔK得到扩展速率与真实的裂纹速率会有误差。为了更好的研究疲劳裂纹的性质,本文通过分析紧凑拉伸(CT 试样的疲劳裂纹扩展后的三维形貌,采用Jiang 一Sehitoglu 循环塑性模型和疲劳准则以及Paris 公式对裂纹扩展速率进行了有限元数值模拟,对扩展形貌的沿厚度方向裂纹的扩展速率进行了分析,发现在裂纹中间区域,疲劳裂纹扩展速率几乎相等。 关键词:循环塑性,疲劳准则,裂纹扩展,Paris 方法 中图分类号:XXX 文献标识码:A 1 引言 对疲劳裂纹的研究一直是固体力学学科中的研究热点之一,因此,为了更好的研究疲劳裂纹的性质,在近来的疲劳裂纹扩展的研究中,对疲劳裂纹的扩展形貌的研究逐渐成为热点。然而,对疲劳裂纹速率的计算方法通常使用Paris 公式来得到,这种方法在计算表面裂纹的扩展速率总被广泛的使用。众所周知,Paris 公式需要得到的疲劳裂纹表面尖端点应力强度因子的变化幅值(ΔK,对于疲劳裂纹面多数为椭圆形貌,裂纹面中间的ΔK 值明显不同于裂纹表面点的ΔK 值。正是如此,在疲劳裂纹扩展速率的研究中,通过Paris 公式得到扩展速率与真实的裂纹速率有误差。所以,为了更好的研究疲劳裂纹的性质,对于半椭圆形状的裂纹形貌下的裂纹扩展研究就显得很有必要。在循环载荷下疲劳裂纹的裂纹形貌在稳定扩展区近似为半椭 圆形状。Putra 和 Schijve [1] 通过电腐蚀方法来研究7075-T6铝板中的疲劳裂纹扩展,观察到疲劳裂 纹的形貌近似于半椭圆形。McFadyen [2] 等通过试验方法研究疲劳裂纹的半椭圆形状在不同应力比下的扩展趋势。 在试验研究的同时,大量数值模拟方法也被应 用到疲劳裂纹扩展的研究当中。Newman 和 Raju [3, 4] 应用双耦合Paris 疲劳裂纹扩展公式,设定一种 特定的半椭圆形貌来模拟疲劳裂纹。Smith 和 Cooper [5] 建立了一种通用的有限元模型,在裂纹形状的变化的时候,自动形成新的疲劳裂纹面及更新 程序。Lin 和 Smith [6] 发展了一种多自由度模拟程序,来预测表面疲劳的裂纹扩展速率及计算由于半椭圆面形状;应力比等的不同而产生的差异。 本文通过分析紧凑拉伸(CT试样的疲劳裂纹扩展后的三维形貌,运用数值模拟方法对扩展形貌的沿厚度方向裂纹的扩展速率进行了分析,得到了在裂纹中间区域,疲劳裂纹扩展速率几乎相等;同时对于CT 试样,由于裂纹扩展面假设满足椭圆方程,在边界处裂纹扩展的法线方向前缘没有单元,或者只存在极少数单元,较难得到理想的裂纹扩展的数值结果。 2 CT 疲劳形貌观察[7] 图1是疲劳试验试件,材料为16MnR 钢,其材料常数如下:220E GPa =; 0.31ν=。试验为载荷控制,其幅值为26.67KN。R (最小载荷与最大载荷之比值为0.1。裂纹形貌确定采用半幅降载勾线法,即经过一定的加载循环次数后,保持勾线阶段的平均载荷与常幅加载的平均载荷一致,其加载幅度降低到常幅加载幅度的一半,使得裂纹扩展速率显著降低,获得一条细而清晰的亮条纹,即裂纹扩展形貌。 图1紧凑拉伸试件(mm Fig.1 Compact tension specimen (mm 图2为试件在横幅载荷作用下的疲劳裂纹扩展结果,从图中可以很明显看出在不同裂纹长度下的疲劳裂纹面具有相同的扩展趋势,在裂纹的扩展方向上每个裂纹面中部都有一定的凸起,从裂纹起裂开始后,每条裂纹面形貌都为近似的半椭圆形状。 图2 0.1R =的常幅加载下的紧凑拉伸试样的裂纹 扩展形貌演化 Fig. 2 Evolution of the fatigue crack growth shape at 1.0=R under constant amplitude loading 选取裂纹长度16.24a mm =和 22.52a mm =的两条裂纹分别拟合出这两条裂纹 的形状曲线,假设曲线的形状为半椭圆形,分别满 足: 1裂纹长度16.24a mm =: 22 11.4212.5z x ????+=???????? (1 2裂纹长度22.52a mm =: 22 11.4412.5z x ???? +=???????? (2 图 3 为试件的裂纹形貌的试验数据形貌与拟合出的椭圆曲线的比较,图中,从图中的比较结果来看,拟合出的裂纹