1. 奇异单元简介
在利用ABAQUS计算裂纹体的应力强度因子时,为了满足裂尖奇异性的要求,通常会在裂纹尖端附近的区域采用奇异单元进行划分。对于平面单元,通过将四边形单元的其中一条单元边折叠,并调节二阶四边形单元中间节点的位置,可以获得满足不同奇异性要求的奇异单元。
(1)将二阶四边形单元的其中一条单元边折叠,同时将折叠边的三个节点绑定到一起,同时将中间节点移动到单元边的1/4位置处,则可以在裂纹尖端获得r-1/2的奇异性:
为了将折叠边的三个节点绑定到一起,可以对三个节点使用线性约束方程或者多点约束,一种更为简便的方法为对折叠边的三个节点使用相同的节点号(需要通过修改inp文件实现),如图1.1所示。
图1.1 奇异单元的构建(r-1/2奇异性)
(2)如果将四边形单元的其中一条单元边折叠,但折叠边的三个节点可以自由移动,并且中间节点仍然保持在单元边的中间位置,如图1.2所示,则可在裂纹尖端获得r-1的奇异性:
图1.2 奇异单元的构建(r-1奇异性)
(3)如果将四边形单元的其中一条单元边折叠,但压缩边的三个节点可以自由移动,并且将中间节点移动到1/4的位置,如图1.3所示,则可在裂纹尖端获得上述奇异性的组合:
图1.3 奇异单元的构建(组合奇异性)
由弹性理论可知,裂纹尖端附近的应力场可统一表示为如下形式:
其中σij为裂尖附近某点处的应力分量,K为应力强度因子,r为该点到裂纹尖端的距离,θ为该点与裂尖构成的连线与裂纹面所呈的夹角,可以看到裂纹尖端的应力分布具有r-1的奇异性,因此预计通过第一种方法构造的奇异单元计算得到的应力强度因子相比于其它两种方法是最为精确的。
在二维中心裂纹板的应力强度因子分析中,笔者已经详细介绍了具有不同奇异性的奇异单元的构建方法。由于在前文中,奇异单元是利用ABAQUS GUI中构建的,因此尽管裂纹尖端采用三角形单元划分,但ABAQUS会自动通过调整单元节点完成奇异单