首页 理论教育裂纹扩展过程中的多种形态模拟

裂纹扩展过程中的多种形态模拟

【摘要】:在裂纹扩展、聚合过程中,裂纹扩展的方式及其形态有时会比较复杂。可以看出,一旦上述两裂纹相交,裂纹4尖端的应力奇异性将消失,裂纹4即会停止继续往前扩展。图3.16裂纹扩展的三种基本形态利用上述模型进行裂纹扩展过程模拟的一个困难是如何实现裂纹聚合的模拟。locked状态意味着该裂纹尖端没有应力奇异性,而unlocked状态意味着该尖端含有奇异性。图3.17数字化裂纹描述图3.18给出了几种典型的裂纹及其可能的扩展形态。

在裂纹扩展、聚合过程中,裂纹扩展的方式及其形态有时会比较复杂。但是,所有的扩展形态都可以归结于三种基本形态或这三种基本形态的组合。这三种基本形态如图3.16所示,在图3.16(a)中,裂纹1与2在某端点相遇,因而实际上融合成为一条裂纹;图3.16(b)显示了两个裂纹在扩展时交叉的情形,此时裂纹4扩展并与裂纹3的裂纹面相交。可以看出,一旦上述两裂纹相交,裂纹4尖端的应力奇异性将消失,裂纹4即会停止继续往前扩展。

图3.16 裂纹扩展的三种基本形态

利用上述模型进行裂纹扩展过程模拟的一个困难是如何实现裂纹聚合的模拟。如上所说,这里的聚合主要包括裂纹的融合和交叉。借助传统有限元实现上述裂纹的聚合模拟,其优点在于可以在商用软件如ABAQUS通过面向对象的Python语言接口进行二次开发来实现,从而为较大规模的混凝土工程结构损伤分析提供便利。

在裂纹扩展过程数值分析中,裂纹面以向量方式近似描述,选取裂纹面上合适距离的几何点(Geometric Points,GP)序列即可较为精确地模拟裂纹构形。以大写的斜体字母下标,如I,J等将用来表示不同裂纹,而斜体下标如i,j等于表示裂纹面上的不同几何点。此外,显然每个裂纹都含有两个端点,这里分别记为头部(Head,H)和尾部(Tail,T)。为了实现裂纹的聚合计算,在每个端点(H或T)上引入两种状态,即锁定的(locked)和自由的(unlocked)。locked状态意味着该裂纹尖端没有应力奇异性,而unlocked状态意味着该尖端含有奇异性。如图3.17给出了上述裂纹在模拟时数字化表示方法的示意图。图中裂纹I通过n+1个连续的几何点(从0到n)来表示。该裂纹含有一个锁住(locked)的头部(用圆圈标注,下同)和一个自由(unlocked)的尾部(用三角形标注,下同)。下文中,斜体大写字符,如I,J,用于描述裂纹,而斜体小写字符,如i和j,用于描述裂纹表面的几何点。

图3.17 数字化裂纹描述

图3.18给出了几种典型的裂纹及其可能的扩展形态。对于两端均位于介质内部的普通裂纹而言,两个端点都被设置成unlocked状态(图3.18(a))。而单边裂纹可以通过锁定其中的一端来实现(图3.18(b))。如果两个裂纹融合,那么融合的端点同时被设定为locked,其余两个端点保持原有状态。

图3.18 几种典型裂纹及其扩展形态

(a)普通裂纹;(b)边裂纹;(c)裂纹融合;(d)裂纹交叉;(e)3裂纹交汇;(f)一条裂纹与另两条裂纹交叉

如图3.18(c)中,裂纹K的T端与裂纹L的T端将被锁定。如果两个裂纹交叉,则交叉点所在的端点将被锁定,其余端点维持原有状态不变。如图3.18(d)中,裂纹N的H端扩展并与裂纹M交叉,此时裂纹N的H端将被设置成锁定状态。还存在其他比较复杂的裂纹型态,同样可以用上述这四种基本形式描述出来。如图3.18(e)所示,三条裂纹在同一点融合。相应地,需要将裂纹R的H端,裂纹S的T端,裂纹T的H端锁定。其余维持原有状态即可。此时实质上形成与图3.18(d)情形相同的裂纹交叉状态,当然描述存在不同。图3.18(f)给出了一个更为复杂的形态,即裂纹V与裂纹W同时与U裂纹交叉。

每一步状态分析之后,会得到最危险裂纹尖端T0(位于裂纹I上)的扩展方向信息。沿此方向做一个射线,探测沿射线方向距离最近的边界(边界包括外边界、内边界以及所有裂纹面),记为LJ。其中J为该边界线段的编号。若大于程序所设置的默认裂纹扩展步长Δadefault,则裂纹按照既定方式扩展;否则若小于Δadefault,则将该裂纹I以及T0与边界线J进行连接。上述探测过程的示意如图3.19所示,其中图3.19(a),(b)是裂纹扩展抵达边界的情形,图3.19(c)、(d)是与其他裂纹聚合的情形。在实际操作中,若该裂纹尖端本身就距离某边界很近(这里指圆心距离而非沿扩展射线探测距离),这时候会直接将该裂纹尖端与该边界连接。

图3.19 裂纹聚合过程的探测

程序中专门发展了一个模块来处理不同裂纹之间的连接问题:比如连接信息的存储,连接信息的更新,节点插入,裂纹聚合探测,边界探测等操作,此处不再赘述。有兴趣的读者可参阅文献[16]。