表面裂纹疲劳扩展的数值模拟.pdfVIP

第! 卷# 第$ 期 应 用 力 学 学 报 %’( !# )( $ !**+ 年,! 月 !#$%% ’()*$+, (- +..,#%/ 0%!+$#! -./( !**+ 文章编号：,***0$11（!**+ ）*$0*2+0*2 表面裂纹疲劳扩展的数值模拟! 吴志学 （扬州大学# 扬州 !!2**1 ） 摘要：建立了一种无形状约束的模拟表面裂纹在线弹性断裂力学条件下疲劳扩展的数值方法，并 研究了表面疲劳裂纹形状演化和裂纹尖端应力强度因子（345 ）的分布特征。该方法以三维有限单 元技术和6789: 疲劳裂纹扩展规律为基础，并在裂纹扩展增量计算中考虑了裂纹闭合影响。本文第 一部分主要介绍模拟三维疲劳裂纹扩展的数值方法的理论背景和相关的技术细节。着重分析和讨 论基于三维有限单元法计算裂纹345 所涉及的几个主要问题：裂纹尖端单元网格密度对估算精度 的影响；自由表面的影响及其修正方法；裂纹尖端非正交单元网格的影响及修正方法。 关键词：应力强度因子；线弹性断裂力学；表面裂纹；疲劳扩展；数值模拟；三维有限元分析 中图分类号： ;$+ ；=,,2# # # 文献标识码：# 则!*G 的345 估算误差将会引起一倍以上的疲劳寿 ,# 引# 言 命偏差。提高实际裂纹的345 估算精度主要有两条 途径：提高数值方法的估算精度和尽可能准确的描 表面裂纹是很多工程结构中常见的缺陷。无论 述实际裂纹形状。显然，后者尤为重要，因为用假设 在损伤容限设计还是在缺陷评定阶段，人们需要尽 的裂纹形状来描述实际裂纹几何势必造成误差，并 可能准确地估算构件中裂纹的345 ，因为这一参量 且由于这一误差的存在，采用高精度的针对假设裂 是预测含裂纹构件的断裂以及估算剩余疲劳寿命的 纹形状的345 解来估算实际裂纹的345 将变得没有 主要依据。一般情况下，只要已知裂纹形状，其345 实际意义。一般情况下，受测试手段限制，准确测量 可以通过多种数值方法估算得到，如,? $ 节点位移 构件中的裂纹形状是困难的，加上345 解对裂纹形 ［,0! ］ ［ ］ ［$ ］ 法 、虚拟裂纹扩展法 、节点力法 等。对于简 状的高度敏感性，近年来，人们开始借助数值模拟技 ［+0,* ］ 单裂纹形状（如半椭圆形裂纹）已有许多经验解可 术来预测疲劳裂纹扩展 ，以期望提高裂纹形状 供使用，如针对有限厚度和宽度平板的表面裂纹的 描述及345 估算精度。H9B 和3A9IJ［+0K ］曾对表面疲 ［2 ］ 劳裂纹形状演化进行了模拟，并就裂纹形状和 ).@A7B 和C7DE 经验公式 。尽管如此，对于结构 345 复杂或受复杂载荷作用的结构件中的非规则形状裂 分布与假设半椭圆形状的裂纹结果进行了比较分 纹，其345 仍然需要通过数值方法估算得到。 析，发现在某些情况下，扩展裂纹形状及其相应的 提高345 的估算精度具有重要的实际意义。例 345 与假设半椭圆形状裂纹的结果之间有较大的偏 如，对于在线弹性条件下扩展的疲劳裂纹，假设裂纹 差。但在他们的数值模型中并未考虑裂纹闭合可能 扩展速率满足6789: