脆性材料在腐蚀环境中表面裂缝扩展的数值研究方法.docVIP

脆性材料在腐蚀环境中表面裂缝扩展的数值研究方法 唐志波 谢永和 楼林燕 浙江海洋学院 船舶与建筑工程学院，316000，舟山 电子信箱: Zhibo_tang@zjou.edu.cn 传真: 0580-2550401 电话: 0580-2550445 摘要:作者采用了一定的数学力学模型，用有限元计算方法模拟了脆性材料表面裂缝在应力和腐蚀耦合作用下的动态扩展过程。我们用front-tracking图形技术开发了有限元网格剖分软件,从根本上解决了该问题数值模拟的关键技术。 关键词:应力腐蚀耦合；裂缝动态扩展；有限元方法 陶瓷材料等高强度脆性材料在切割、打磨等加工制作工程中，难免表面会形成很多微裂缝。在持续应力作用下并且在腐蚀环境工作时，这些微裂缝会向内扩展，裂缝长度不断增加，元器件的有效受力截面逐渐变小，最终导致断裂破坏，这种现象称为应力腐蚀疲劳问题。这种方式的裂缝扩展主要是因为裂缝尖端附近应力集中的区域材料被长时间快速腐蚀引起，其中腐蚀的速度跟裂缝附近的应力状态、材料常数、裂缝的几何形状以及化学反应环境有关。实验表明对于给定的材料和化学反应环境，裂缝扩展最终导致材料破坏的时间长短取决于外荷载的大小。但特别要提到的是，当裂尖初始应力小于某个阀值的时候，由于裂缝的钝化现象，这种应力腐蚀引起的断裂破坏现象可以避免。近几十年来有很多学者对这种现象发生了兴趣，纷纷提出了各种解释该现象的机理和数学力学模型[2-10]，其中Charles-Hillig 在1965年提出的模型[1]被广为接受。1991年，美国国家计量局的学者Chuang 和Fuller[11]在该模型的基础上做出了杰出的工作，他们假设了规则形状的微裂缝，并根据在应力和腐蚀耦合作用下裂缝表面被腐蚀掉（表现为裂缝扩展）的初始条件，从理论上证实了应力阀值的存在，并预测了裂缝扩展过程中的四种可能形态。但因Charles-Hillig模型的复杂性，在裂缝扩展数值模拟上一直没有大的突破。在本文中，作者报道了应用自己开发的有限元程序进行的裂缝扩展全过程数值模拟工作，计算结果和文献中的试验结果有较好的相似性。 1 理论模型 图1应力腐蚀模型示意图 图 2 数值计算模型示意图 Fig.1 Sketch map for stress corrosion problem Fig.2 Idealized geometry and six nodded triangle elements Charles 和 Hillig 提出的模型描述的裂缝表面（包括裂缝尖端）材料被腐蚀的速度乃取决于裂缝附近局部的“热力学驱动力”（thermodynamic driving force）和能障（energy barrier）的变化，而这两个因素又取决于局部的应力状态和裂缝表面的曲率分布，其模型见图1所示，图中，分别为弹性模量和泊桑比，是表面曲率，σ是裂缝表面的切向应力。为研究方便起见，作者采用了平面(平面应变)假设，裂缝表面包括裂缝尖端的法向倒退速度（也即材料表面被腐蚀的速度）可以表示为[1] （1） 表示在无应力状态下表面材料被腐蚀的速度，一般在 ( 取负值是因为表面总是因物质被腐蚀而后退)。另外, （2） 其中是裂缝表面的切向应力；是材料的摩尔体积；和是两个无量纲的常数（其物理意义及相关热力学推导见[1]）；是平面应变假设的弹性模量；和分别是无应力状态下裂缝单位面积上的表面能和表面曲率，表示凹平面。 我们假设一个边缘包含一条微裂缝的脆性体，裂缝处于化学腐蚀环境，脆性体在远端受垂直于裂缝长度方向的力。典型的有限元单元剖分图见图2，其中我们根据对称性只取一半进行计算，图中剖分的单元均为6节点的三角形单元，微裂缝初始形状为扁平的半椭圆，两半轴长度之比为0.01，远端有垂直于裂缝长轴的载荷。这样，我们就可以得到相应的位移场和应力场，其本构关系为： （3） 同时应力场必须满足： （4） 在接下来的讨论中，我们假设弹性体在荷载作用下的弹性变形忽略不计，所以裂缝形状和长度的变化被认为完全由裂缝表面的材料被腐蚀所致。 2 数值模拟方法 我们应用有限元方法来计算上述跟时间有关的应力腐蚀耦合问题，大致分为三个步骤。首先计算时的应力场，得到裂缝表面各点的切向应力分量值。然后计算在微小时间步后的裂缝形状，得到新的计算图形。最后在新的计算域内剖分单元，重新计算应力场。再重复前两步工作，循环计算，直到指定的时间为止，同时记录并用TechPlot显示裂缝随时间扩展的全过程。 我们应用了标