Study on the Influence of Crack Growth Increment on Crack Growth Stability

裂纹扩展增量对裂纹扩展稳定性影响的研究 Study on the Influence of Crack Growth Increment on Crack Growth Stability first-author 朱凯林（1998—），男，硕士研究生，主要从事轨道车辆动力学结构强度及可靠性方面的研究。E-mail：1612504315@qq.com。 朱凯林（1998—），男，硕士研究生，主要从事轨道车辆动力学结构强度及可靠性方面的研究。E-mail：1612504315@qq.com。 1612504315@qq.com 西南交通大学 机械工程学院，四川 成都 610031 西南交通大学 机械工程学院，四川 成都 610031 School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China 在实际的裂纹扩展仿真分析中需要确定三个主要参数，包括局部裂纹扭转角度、局部裂纹扩展增量以及扩展后裂纹尖端的拟合外伸，而对于裂纹扩展增量往往根据计算结果精度的需求来确定其大小，忽略了裂纹扩展增量对裂纹扩展稳定性的影响。针对这一问题，利用FRANC3D和Abaqus软件联合仿真分析计算了Q235钢CT试样在不同裂纹扩展增量情况下的裂纹扩展特点，并根据裂纹尖端应力强度因子的分布情况得出两点结论：（1）裂纹两端接近结构自由表面，两端趋向于平面应力状态，中间内部裂纹处于平面应变状态，导致应力强度因子SIF呈现中间大、两端小的对称分布，使得在较大裂纹扩展增量的情况下裂纹中部与两端的扩展值差异较大，最终导致裂纹凸凹扩展甚至失稳；（2）在进行裂纹扩展分析时，为保证计算精度和节约计算时间，并且保证裂纹扩展的稳定性，裂纹扩展增量选取为当前裂纹尖端单元特征尺寸的15%~30%比较合适。 In an actual crack extension needs to determine three main parameters in the simulation analysis, including local crack torsion Angle, local crack growth increment and expanded fitting overhang of the crack tip, and for the crack growth increment is often according to the requirements of the precision of the calculation results to determine its size, ignoring the effect of crack growth increment on the stability of crack propagation. To solve this problem, FRANC3D and Abaqus software were used to analyze and calculate the crack growth characteristics of Q235 steel CT sample under different crack growth increments, and two conclusions were drawn according to the distribution of stress intensity factor at the crack tip: (1) on both ends of the crack near free surface structure, both ends of trend surface stress state, intermediate internal crack in the plane strain state, the stress intensity factor SIF presents small symmetry distribution on both ends of the large, middle, made in the case of large crack growth increment at the ends of the central crack and extension value differences, eventually lead to crack the concavo-convex extension even instability; (2) In the analysis of crack growth, in order to ensure the calculation accuracy, save calculation time and ensure the stability of crack growth, it is appropriate to select 15%~30% of the characteristic size of the current crack tip element for the increment of crack growth. FRANC3D CT试样 应力强度因子 裂纹扩展增量 FRANC3D computed tomography (CT) specimen stress intensity factor crack growth increment 前言 1 工程上构件疲劳失效的主要原因是构件危险部位裂纹的萌生及扩展，目前机械结构疲劳裂纹剩余扩展寿命计算都是基于Paris公式，计算出裂纹尖端的应力分布情况和应力强度因子（stress intensity factor，SIF）大小，结合相应工程判据计算出裂纹尖端扭转角度和材料断裂韧性来确定裂纹扩展速率及其扩展寿命 ［ 1 1 ］ 。 利用裂纹尖端应力计算SIF需要裂纹尖端具有非常细密的网格，而常用的有限元分析软件如ANSYS、Abaqus和Nastran等难以对局部网格进行理想的细化。FRANC3D是FAC公司开发的新一代裂纹分析软件，不仅可以对局部裂纹尖端进行网格细化，还可以在裂纹尖端生成三圈环状奇异单元，提高计算精度，并且提供了ANSYS、Abaqus和Nastran等软件接口 ［ 2 2 ］ 。 目前有关动态裂纹的研究最经典的案例是利用拉格朗日有限差分法分析阶跃载荷作用下中心裂纹板的动态应力强度因子。彭英 ［ 3 3 ］ 应用相互作用积分法计算动态应力强度因子，考虑了裂纹面接触作用、冲击载荷形式、裂纹长度、裂纹角度及裂纹位置等参数的影响；姜伟 ［ 4 4 ］ 采用将动态有限元分析过程和相互作用积分相结合的方法，在有限元分析软件ANSYS中求解动态应力强度因子。尽管动态求解方法更趋近于实际，但其计算过程极其复杂，不适用于复杂模型，故常用的疲劳裂纹扩展分析方法是通过设置不同的静态裂纹分析步进行迭代，在相邻两分析步中，很重要的一个参数是裂纹扩展增量。扩展增量越大，裂纹扩展分析步较少，计算时间短；增量越小，裂纹扩展分析步较多，计算精度高，但计算时间较长。故裂纹扩展增量的确定对裂纹扩展仿真分析有重要意义。而在实际的裂纹扩展仿真分析中，诸多研究学者往往根据结果精度的需求来确定裂纹扩展增量的大小，忽略了扩展增量对裂纹扩展稳定性的影响，为探究这一问题，本文结合FRANC3D软件，分析了在不同扩展增量下紧凑拉伸（CT）试样的裂纹扩展稳定性情况，并根据贯穿裂纹SIF的分布特点分析了裂纹凸凹扩展的原因及裂纹扩展增量的确定方法，为裂纹扩展剩余寿命计算问题提供参考依据。 疲劳裂纹扩展分析方法 1 对结构进行疲劳裂纹扩展分析首先是建立模型，本文采用HyperMesh对模型进行网格划分并将其导入Abaqus定义材料属性、载荷以及约束，并提交分析。然后导入FRANC3D中在裂纹扩展区域引入裂纹并重新划分网格，得到新的完整有限元模型后再提交到Abaqus中进行应力分析，计算出裂纹尖端的SIF值，根据SIF大小计算出裂纹尖端节点的局部扭转角度，结合给定的裂纹扩展距离或载荷循环次数确定裂纹扩展几何并再次重新划分裂纹网格，如此反复，便可得到裂纹每一步扩展的应力、SIF以及扩展寿命的详细数据。具体的疲劳裂纹扩展仿真模拟流程如 图1 图1 所示。 疲劳裂纹扩展过程 1 对裂纹进行疲劳裂纹扩展分析时需要确定三个重要参数，分别为局部裂纹扭转角度、局部裂纹扩展增量以及扩展后裂纹尖端的拟合外伸，如 图2 图2 所示。 计算裂纹尖端局部扭转角度 2 通过Abaqus应力分析后，根据 式（1） 式（1） 计算出裂纹尖端的SIF大小。 <math display="block" id="M1"><msub><mrow><mi>K</mi></mrow><mrow><mi>m</mi></mrow></msub><mo>=</mo><mi>F</mi><msup><mrow><mi>σ</mi></mrow><mrow><mi>m</mi></mrow></msup><mroot><mrow><mi>π</mi><mi mathvariant="normal">a</mi></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034221&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034202&type= 20.74333191 4.99533367 式中　 m 取Ⅰ、Ⅱ、Ⅲ，分别表示三种裂纹类型； a 为裂纹尺寸； F 为与裂纹尺寸相关的系数。结合工程判据计算出裂纹尖端各个节点的偏转角度（见 图2 图2 中的 <math id="M2"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034224&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034223&type= 1.86266661 2.62466669 角）。其工程判据有最大张应力准则、局部对称原则、最大能量释放率准则以及最小应变能密度准则等，这些判据本质上相差不大，都是倾向于最小化Ⅱ型裂纹应力强度因子的准则，均可得到准确且相似的结果，本文采用最大张应力准则 ［ 5 5 ］ ，如 图3 图3 所示，其裂纹偏转角 <math id="M3"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034224&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034223&type= 1.86266661 2.62466669 计算公式见 式（2） 式（2） 。 <math display="block" id="M4"><mi>θ</mi><mo>=</mo><mfrac><mrow><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">I</mi></mrow></msub></mrow><mrow><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">I</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034206&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034205&type= 10.58333302 8.63599968 最大张应力准则的基本假设如下： （1）裂纹沿最大周向应力 <math id="M5"><msub><mrow><mi>σ</mi></mrow><mrow><mi>θ</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034231&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034229&type= 6.34999990 2.62466669 的方向开裂； （2）当此方向的周向应力到达临界值时，裂纹失稳扩展。 计算裂纹尖端局部扩展距离 2 得到裂纹尖端的扭转角后，还需确定裂纹尖端节点的扩展距离（见 图2 图2 中的 <math id="M6"><mo>∆</mo><mi>a</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034354&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034353&type= 3.97933316 2.70933342 ），从而确定扩展后裂纹尖端的位置。但裂纹尖端各个节点的SIF大小不一致，裂纹扩展时各个节点的扩展距离也不一致，且裂纹尖端局部扩展距离是针对一个节点而言，为了确定所有节点的扩展距离，需要指定一个节点的扩展距离，其余节点扩展距离根据其SIF值与指定节点SIF的比值大小来确定。由于裂纹两端接近表面，局部不再保持平面应变条件 ［ 6 6 ］ ，计算出的SIF精度较差，为了减小误差，通常指定的节点都选取为裂纹中部且SIF值接近于平均应力强度因子的节点，其扩展距离为 <math id="M7"><mo>∆</mo><msub><mrow><mi>a</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">d</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034246&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034244&type= 7.19666624 2.70933342 ，其余节点 <math id="M8"><mo>∆</mo><mi>a</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034354&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034353&type= 3.97933316 2.70933342 由 式（1） 式（1） 确定。 <math display="block" id="M9"><mo>∆</mo><mi>a</mi><mo>=</mo><mo>∆</mo><msub><mrow><mi>a</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">d</mi></mrow></msub><mfrac><mrow><mo>∆</mo><mi>k</mi></mrow><mrow><mo>∆</mo><msub><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">d</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034255&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034252&type= 22.43666649 8.38199997 扩展后裂纹的拟合及外伸 2 扩展后的裂纹尖端由各节点扩展后的位置连接确定，但直接相连接会形成锯齿形裂纹尖端（如 图2 图2 中的蓝虚线所示），形成数值“噪音”，这种非光顺的裂纹在下一步裂纹扩展计算出的SIF曲线也是锯齿形的，导致后续裂纹扩展极其不稳定。为了消除数值“噪音”的影响，需要对每一步扩展后的裂纹尖端进行光顺化处理（如 图2 图2 中红实线所示）。 裂纹扩展增量对裂纹扩展稳定性的影响 1 在指定裂纹尖端局部扩展距离时，有以下两种确定方法。 （1）指定载荷循环次数来推算。 根据指定的载荷循环次数来推算出裂纹扩展距离，但指定的载荷循环次数在裂纹扩展初期可能偏小，导致扩展距离步过小，计算时间加长，效率低。而对于裂纹扩展后期而言又可能过大，裂纹扩展过快导致失稳断裂。 （2）指定扩展增量计算。 指定的扩展增量若过小，裂纹扩展过程会较为缓慢，计算时间虽然加长，但裂纹前缘相对比较稳定。若裂纹扩展增量过大，由于裂纹两端的SIF值精度较低，会加剧裂纹两端与中间部分扩展增量的差值，造成裂纹的不稳定扩展。 本文以CT拉伸试样为例，分析不同裂纹扩展增量对裂纹扩展稳定性的影响及其裂纹扩展特点，根据GB/T 6398—2017《金属材料疲劳试验疲劳裂纹扩展方法》 ［ 7 7 ］ ，确定CT试样的基本外形尺寸参数如 图4 图4 所示，其厚度为7.5 mm。 在Hypermesh中构建有限元网格模型如 图5 图5 所示，随后导入Abaqus中定义载荷及约束条件，分别将上、下加载孔心与其孔表面耦合在一起，对上加载孔心施加6 kN向上的拉力，对下加载孔孔心的 XYZ 位移及其转动全部约束住。 材料选取为低碳钢Q235，参考董达善 ［ 8 8 ］ 和赵章焰 ［ 9 9 ］ 等用三点弯曲试验测定的Q235钢材料的疲劳裂纹扩展性能参数，得到材料断裂参数 C 、 m 、应力强度因子阀值 <math id="M10"><mo>∆</mo><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034234&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034267&type= 6.26533365 2.70933342 以及断裂韧性 <math id="M11"><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">C</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034238&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034236&type= 3.97933316 2.62466669 如 表1 表1 所示。 屈服强度 R eL /MPa 抗拉强度 R m /MPa C m SIF阀值 <math id="M12"><mo>∆</mo><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034242&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034240&type= 4.65666676 1.94733346 /MPa·mm 1/2 材料断裂韧性 <math id="M13"><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">C</mi></mrow></msub><mo>/</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034276&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034274&type= 3.97933316 1.86266661 MPa·mm 1/2 弹性模量 D /MPa 235 400 4.552 <math id="M14"><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>-</mo><mn>17</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034281&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034279&type= 6.60400009 1.86266661 4.198 7 99.93 6052.64 2.08 <math id="M15"><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>5</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034285&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034283&type= 4.91066647 1.86266661 在Abaqus中提交应力分析后导入FRANC3D软件中，在CT试样缺口处预设一条长8 mm的贯穿裂纹并重新绘制局部裂纹网格，得到如 图6 图6 所示的初始裂纹尖端应力分布。 根据GB/T 6398-2017标准规定，设置其应力比 R =0.1，分别在0.2 mm、0.5 mm、1.0 mm、1.5 mm以及2.0 mm的裂纹扩展增量（ <math id="M16"><mo>∆</mo><mi>a</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034354&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034353&type= 3.97933316 2.70933342 ）情况下对裂纹进行扩展分析，得到每种情况下的裂纹扩展如 图7 图7 所示。 growth increments 1;2;3;4;5; 为了比较在不同 <math id="M17"><mo>∆</mo><mi>a</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034354&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034353&type= 3.97933316 2.70933342 情况下的裂纹扩展状况的好坏，将裂纹扩展过程中裂纹尖端中部SIF值的稳定性作为评估标准，GB/T 6398-2017标准中给出了CT试样SIF值的计算公式 ［ 10 10 ］ 如 式（2） 式（2） 所示，将其计算结果作为理论参考值。 <math display="block" id="M18"><mi>K</mi><mo>=</mo><mfrac><mrow><mi>F</mi></mrow><mrow><mi>B</mi><mroot><mrow><mi>W</mi></mrow><mrow/></mroot></mrow></mfrac><mi mathvariant="normal">g</mi><mfenced separators="|"><mrow><mfrac><mrow><mi>a</mi></mrow><mrow><mi>w</mi></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034344&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034343&type= 26.41600037 9.31333351 式中　 F 为施加在上加载孔的拉力； B 为CT试样厚度； W 为试样宽度，即从上下加载孔心连线到试样边缘的距离； a 为从上下加载孔心连线到裂纹尖端的垂直距离； <math id="M19"><mi mathvariant="normal">g</mi><mfenced separators="|"><mrow><mfrac><mrow><mi>a</mi></mrow><mrow><mi>w</mi></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034345&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034362&type= 8.80533314 8.29733276 为应力强度因子形状参数，其大小由 式（3） 式（3） 确定。 <math display="block" id="M20"><mi>g</mi><mfenced separators="|"><mrow><mi>α</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><mn mathvariant="normal">0.886</mn><mo>+</mo><mn mathvariant="normal">4.64</mn><mi>α</mi><mo>-</mo><mn mathvariant="normal">13.32</mn><msup><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mn mathvariant="normal">14.72</mn><msup><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5.6</mn><msup><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>-</mo><mi>α</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">1.5</mn></mrow></msup></mrow></mfrac><mo>×</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034368&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034366&type= 70.52732849 8.55133343 <math display="block" id="M21"><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo>+</mo><mi>α</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034371&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034378&type= 9.05933285 3.38666677 需要注意的是，标准中公式计算出的理论SIF值代表整个裂纹尖端的值，即认为CT试样贯穿裂纹的裂纹尖端SIF值处处相等。而在实际仿真计算中，贯穿裂纹两端为自由表面，靠近两端的裂纹尖端逐渐趋向于平面应力状态，中间内部裂纹处于平面应变状态，使得实际仿真计算出的CT试样贯穿裂纹SIF值会出现两端小、中间大的情况，长8 mm的预制裂纹尖端SIF值及在 <math id="M22"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">2.0</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034464&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034445&type= 18.37266731 2.79399991 时后续4步裂纹扩展过程与理论参考值的比较如 图8 图8 所示。 由 图7 图7 可知， <math id="M23"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">0.2</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034450&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034449&type= 18.37266731 2.79399991 、 <math id="M24"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">0.5</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034491&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034477&type= 18.37266731 2.79399991 以及 <math id="M25"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">1.0</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034488&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034475&type= 18.37266731 2.79399991 时，每一步扩展拟合后的裂纹尖端形状都比较平稳，裂纹扩展过程稳定； <math id="M26"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">1.5</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034444&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034460&type= 18.37266731 2.79399991 以及 <math id="M27"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">2.0</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034464&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034445&type= 18.37266731 2.79399991 时，裂纹扩展尖端明显趋向凸凹状，扩展极其不稳定，甚至在扩展后期由于裂纹尖端形状异常导致重新引入裂纹网格的困难，无法进行下一步扩展分析，需要修正裂纹尖端的拟合形状。由 图8 图8 可知，贯穿竖直裂纹尖端中部SIF大、两端SIF小，导致在扩展时中部扩展距离长、两端短，形成凸状裂纹尖端。而凸状裂纹尖端两端SIF大、中间SIF小，导致在扩展时两端距离长、中间短，若扩展增量合适，裂纹尖端会重新形成竖直贯穿裂纹（见 图7 图7 中 <math id="M28"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">1.0</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034488&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034475&type= 18.37266731 2.79399991 ）；若扩展增量较大，则会形成凹状裂纹尖端（见 图7 图7 中 <math id="M29"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">2.0</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034464&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034445&type= 18.37266731 2.79399991 ）。 为弱化裂纹两端SIF值的影响，选取裂纹尖端中部计算结果为评估值，其与理论计算结果的误差为7.59%，主要来源于标准给出的理论计算公式中假定了CT试样为刚性体，未考虑外加载荷和边界条件带来的影响 ［ 11 11 - 12 12 11-12 ］ 。而仿真计算时为了提高计算精度，在裂纹尖端还设置了三圈奇异单元环，使得理论计算结果要小一些，与实际结果存在一定偏差。将不同扩展增量下每一扩展步的SIF结果与理论计算结果进行比较，如 图9 图9 所示。 由 图9 图9 可知，在小的裂纹扩展增量情况下，裂纹中部尖端SIF值与理论参考值接近重合，存在一定可接受范围内的误差，随着扩展增量的加大，尤其是在 <math id="M30"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">1.5</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034444&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034460&type= 18.37266731 2.79399991 和 <math id="M31"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">2.0</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034464&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034445&type= 18.37266731 2.79399991 时，裂纹扩展变得逐渐不稳定，呈现凸凹状，且SIF值也在理论参考值附近波动；在 <math id="M32"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">0.2</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034450&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034449&type= 18.37266731 2.79399991 和 <math id="M33"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">0.5</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034491&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034477&type= 18.37266731 2.79399991 时，裂纹扩展平稳，计算结果与理论值也相差不大； <math id="M34"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">1.0</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034488&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034475&type= 18.37266731 2.79399991 的裂纹扩展后期逐渐出现不稳定扩展趋势。为了在保证计算精度的前提下需尽可能减少计算时间，根据本文数据可知，裂纹扩展增量 <math id="M35"><mo>∆</mo><mi>a</mi><mo>=</mo><mn mathvariant="normal">0.5</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034491&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034477&type= 18.37266731 2.79399991 时裂纹扩展稳定、计算时间短并且计算精度高。本文模型的裂纹尖端单元特征尺寸为2 mm，结合管建东 ［ 12 12 ］ 的网格参数对FRANC3D软件计算SIF精度的影响研究，裂纹扩展增量选取为当前裂纹尖端单元特征尺寸的15%~30%比较合适。 结论与展望 1 本文通过利用FRANC3D和Abaqus软件联合仿真分析了Q235钢CT试样在不同裂纹扩展增量情况下的裂纹扩展稳定性，得出了以下结论： （1）在实际仿真计算中，贯穿裂纹两端为自由表面，由于靠近两端的裂纹尖端逐渐趋向于平面应力状态，中间内部裂纹处于平面应变状态，导致应力强度因子SIF呈现中间大、两端小的对称分布，进而使得裂纹扩展时中间部分裂纹尖端扩展距离长，形成凸状裂纹尖端。凸状裂纹尖端中间SIF小、两端大，在大的裂纹扩展增量情况下裂纹形成凹状裂纹尖端，凸凹交替的裂纹扩展极其不稳定，甚至会影响下一扩展步的网格划分。 （2）在进行裂纹扩展分析时，为保证计算精度和节约计算时间，并且保证裂纹扩展的稳定性，裂纹扩展增量选取为当前裂纹尖端单元特征尺寸的15%~30%比较合适。 本文研究的对象是标准拉伸试样，裂纹为贯穿裂纹，属于Ⅰ型裂纹扩展，而裂纹扩展增量对复合型裂纹扩展结构的影响未曾分析，还有待研究。 参考文献 孙远韬 ， 王少梅 . 基于SIF幅值拟合的机械结构裂纹扩展寿命分析方法研究 ［J］. 武汉理工大学学报：交通科学与工程版 ， 2008 ， 32 （ 5 ）： 4 . 基于SIF幅值拟合的机械结构裂纹扩展寿命分析方法研究 Sun Yuantao ， WANG Shaomei . Analysis of Crack Growth Life of Mechanical Structures Based on SIF Amplitude Fitting ［J］. Journal of Wuhan University of Technology： Transportation Science & Engineering ， 2008 ， 32 （ 5 ）： 4 . Analysis of Crack Growth Life of Mechanical Structures Based on SIF Amplitude Fitting 艾书民 ， 于明 ， 成晓鸣 ， 等 . 基于Franc3D软件的三维裂纹扩展分析与应用 ［J］. 机械强度 ， 2018 ， 40 （ 1 ）： 251 - 254 . 基于Franc3D软件的三维裂纹扩展分析与应用 AI Shumin ， YU Ming ， Cheng Xiaoming ， et al . Analysis and application of three-dimensional crack propagation based on Franc3D software ［J］. Journal of mechanical strength ， 2018 ， 40 （ 1 ）： 251 - 254 . Analysis and application of three-dimensional crack propagation based on Franc3D software 彭英 ， 杨平 ， 姜伟 . 考虑接触的斜裂纹动态应力强度因子分析 ［J］. 哈尔滨工程大学学报 ， 2019 ， 40 （ 8 ）： 7 . 考虑接触的斜裂纹动态应力强度因子分析 Peng Ying ， Yang Ping ， Jiang Wei . Dynamic Stress intensity factor analysis of Oblique Crack considering Contact ［J］. Journal of Harbin Engineering University ， 2019 ， 40 （ 8 ）： 7 . Dynamic Stress intensity factor analysis of Oblique Crack considering Contact 姜伟 ， 杨平 ， 董琴 . 平板穿透裂纹尖端动态应力强度因子研究 ［J］. 武汉理工大学学报（交通科学与工程版） ， 2016 ， 40 （ 5 ）： 820 - 825 . DOI： 10.3963/j.issn.2095-3844.2016.05.012 http://dx.doi.org/10.3963/j.issn.2095-3844.2016.05.012 . 平板穿透裂纹尖端动态应力强度因子研究 Jiang Wei ， Yang Ping ， DONG Qin . Journal of wuhan university of technology （transportation science and engineering） ， 2016 ， 40 （ 5 ）： 820 - 825 . DOI： 10.3963/j.issn http://dx.doi.org/10.3963/j.issn . Journal of wuhan university of technology （transportation science and engineering） 2095 - 3844 . 2016 . 05 . 012 . 2016 王铎 . 断裂力学（上）基本部分 ［M］. 哈尔滨 ： 哈尔滨工业大学出版社 ， 1989 . 断裂力学（上）基本部分 JRP Corrêa ， Junior J S ， Carrasco J ， et al . Structural integrity assessment of a 5CT P110 steel riser pipe according to BS 7910：2013 standard ［J］. International Journal of Pressure Vessels and Piping ， 2020 ， 188 （ 3 ）： 104206 . Structural integrity assessment of a 5CT P110 steel riser pipe according to BS 7910：2013 standard 钢铁研究总院 . 金属材料疲劳试验疲劳裂纹扩展方法 ： GB/T6398-2017 ［S］. 北京 ： 人民交通出版社 ， 2017 . 董达善 ， 梅潇 . Q235材料的Paris常数C，m及 <math id="M36"><mo>∆</mo><msub><mrow><mi mathvariant="normal">K</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034480&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034478&type= 5.58799982 2.28600001 值测试 ［J］. 机械强度 ， 2003 ， 25 （ 2 ）： 215 - 218 . Q235材料的Paris常数C，m及∆Kth值测试 Dong DA-shan ， MEI Xiao . Paris Constant C， m and <math id="M37"><mo>∆</mo><msub><mrow><mi mathvariant="normal">K</mi></mrow><mrow><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034500&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034497&type= 5.58799982 2.28600001 values of Q235 material ［J］. Journal of mechanical strength ， 2003 ， 25 （ 2 ）： 215 - 218 . Paris Constant C， m and ∆Kth values of Q235 material 赵章焰 ， 吕运冰 ， 孙国正 . J积分法测量低碳钢Q235的断裂韧性 <math id="M38"><msub><mrow><mi mathvariant="normal">K</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">C</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034521&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034518&type= 4.23333359 2.20133328 ［J］. 武汉理工大学学报 ， 2002 ， 24 （ 4 ）： 111 - 112 . J积分法测量低碳钢Q235的断裂韧性KIC Zhao Zhangyan ， Lv Yunbing ， Sun Guozheng . Measurement of fracture toughness <math id="M39"><msub><mrow><mi mathvariant="normal">K</mi></mrow><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">C</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034521&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034518&type= 4.23333359 2.20133328 of Low Carbon Steel Q235 by J Integral Method ［J］. Journal of wuhan university of technology ， 2002 ， 24 （ 4 ）： 111 - 112 . Measurement of fracture toughness KIC of Low Carbon Steel Q235 by J Integral Method Barrionuevo G O ， Guerrero B . ANLISIS NUMRICO DE LA MECNICA DE FRACTURA EN TUBERAS FORZADAS APLICANDO EL ESTNDAR BS7910 Y SIMULACIN POR ELEMENTOS FINITOS ［C］. III Congreso de Estudiantes de Ingeniería UC-USM ， 2018 . ANLISIS NUMRICO DE LA MECNICA DE FRACTURA EN TUBERAS FORZADAS APLICANDO EL ESTNDAR BS7910 Y SIMULACIN POR ELEMENTOS FINITOS 熊勋 ， 杨莹 ， 汪舟 ， 等 . 基于FRANC3D和ABAQUS联合仿真三维疲劳裂纹扩展分析及寿命预测 ［J］. 武汉理工大学学报（交通科学与工程版） ， 2020 ， 44 （ 03 ）： 506 - 512 . 基于FRANC3D和ABAQUS联合仿真三维疲劳裂纹扩展分析及寿命预测 Xiong Xun ， Yang Ying ， Wang Zhou ， et al . Three-dimensional Fatigue Crack Propagation Analysis and Life Prediction Based on FRANC3D and ABAQUS Co-Simulation ［J］. Journal of Wuhan University of Technology （Transportation Science and Engineering） ， 2020 ， 44 （ 03 ）： 506 - 512 . Three-dimensional Fatigue Crack Propagation Analysis and Life Prediction Based on FRANC3D and ABAQUS Co-Simulation 管建东 . FRANC3D计算应力强度因子的可靠性及网格参数对精度影响的研究 ［J］. 建筑工程技术与设计 ， 2018 （ 04 ）： 2473 . FRANC3D计算应力强度因子的可靠性及网格参数对精度影响的研究 Guan Jiandong . Reliability of FRANC3D Calculation of Stress intensity factor and Influence of Grid Parameters on accuracy ［J］. Building Engineering Technology and Design ， 2018 （ 04 ）： 2473 . Reliability of FRANC3D Calculation of Stress intensity factor and Influence of Grid Parameters on accuracy 编辑部网址： http：//www.71dhj.com http://www.71dhj.com 编辑部网址： 10.7512/j.issn.1001-2303.2022.03.09 TG405 A 1001-2303（2022）03-0065-07 朱凯林，米彩盈.裂纹扩展增量对裂纹扩展稳定性影响的研究［J］.电焊机，2022，52（3）：65-71. 朱凯林，米彩盈.裂纹扩展增量对裂纹扩展稳定性影响的研究［J］.电焊机，2022，52（3）：65-71. ZHU Kailin, MI Caiying.Study on the Influence of Crack Growth Increment on Crack Growth Stability[J].Electric Welding Machine, 2022, 52(3): 65-71. ZHU Kailin, MI Caiying.Study on the Influence of Crack Growth Increment on Crack Growth Stability[J].Electric Welding Machine, 2022, 52(3): 65-71. 2021-09-03 2021-11-12 2022-03-30 2022-04-26 电焊机 3 52